btc address: 1FmWXNJT3jVKaHBQs2gAs6PLGVWx1zPPHf
eth address: 0xd91c747b4a76B8013Aa336Cbc52FD95a7a9BD3D9

上一篇文章中，ECC椭圆曲线加密(一) 介绍了椭圆曲线的加法，乘法。

同余运算

同余就是有相同的余数，两个整数 a、 b，若它们除以正整数 m所得的余数相等，则称 a， b对于模m同余，用表示。
在数论中叫 “时钟运算”。 这个 “时钟运算” 跟CPU时钟周期没有关系。

python中用% 计算同余，如何计算 呢?

2**137%71 
#或者
pow(2,137,71) 
12%5
27%5

有限域

上一篇中的椭圆曲线，对坐标（x，y）没有任何限制，只要符合曲线方程就可以，坐标可以是整数、负数、有理数，即在实数范围内，实数用 表示。

椭圆曲线是连续的，并不适合用于加密；所以，我们必须把椭圆曲线变成离散的点，我们要把椭圆曲线定义在有限域上。
而椭圆曲线密码所使用的椭圆曲线是定义在有限域内，有限域用 表示。
有限域最常见的例子是，当元素为质数时的有限域（用表示）所组成的整数集合。

假设在 中的椭圆曲线 是什么计算规则呢？
，用到同余的概念。
我们找到3个坐标 (192, 105) (17, 56), (1, 193)都在有限域中。

#计算他们的余数是否相同
(192**3+7)%223 == 105**2%223
(17**3+7)%223 == 56**2%223
(1**3+7)%223 == 193**2%223

椭圆曲线上的离散对数

在椭圆曲线密码中，我们首先定义一条椭圆曲线，然后对椭圆曲线上的某一点之间的运算进行定义，并用这些运算来进行密码技术的相关计算，这就是椭圆曲线加密算法的数学依据。

如果椭圆曲线上一点P，我们计算nP，显然点的分布与顺序都是杂乱无章的。

乘法逆元

在模7乘法中：

• 1的逆元为1 (1*1)%7=1
• 2的逆元为4 (2*4)%7=1
• 3的逆元为5 (3*5)%7=1
• 4的逆元为2 (4*2)%7=1
• 5的逆元为3 (5*3)%7=1
• 6的逆元为6 (6*6)%7=1

扩展欧几里得算法用来求乘法逆元

在 mod p 的意义下我们把x的乘法逆元写作 ，乘法逆元有如下的性质：

乘法逆元的一大应用是模意义下的除法，除法在模意义下并不是封闭的，但我们可以根据上述公式，将其转化为乘法。
假设需要1/4 mod 23，可以转化为，又可以转化为1*(4和23的乘法逆元) mod 23。

def ext_euclid(a, b):
     if b == 0:
         return 1, 0, a
     else:
         x, y, q = ext_euclid(b, a % b) # q = gcd(a, b) = gcd(b, a%b)
         x, y = y, (x - (a // b) * y)
         return x, y, q

标量乘法

除了加法，我们定义另外一种运算：标量乘法，也即数乘

写成如上形式的话，nP 的计算看上去需要 n 次加法。如果 n 有 k 位二进制位的话，即位，那我们的算法复杂度就是, 计算量有点大，但是其实存在更快速的方案。

其中一个就是先做倍数再做加法。假设n=151，其对应的二进制是10010111。而二进制数字可以转化为：

我们可以这么写：

所以，该运算过程是这样的：

• 获取P
• 取P的2倍，得到2P
• 2P加上P
• 把2P再取2倍，得到4P
• 4P加上2P加上P
• 4P再取2倍，得到8P
• 不取8P做运算
• 8P取2倍，得到16P
• 16P加上4P加上2P加上P
• ……

最终，要得到151P我们只是做了一些简单的倍数以及加法。

如果我们计算 6P，只需要

• P
• 通过P得到2P
• 通过2P得到4P
• 2P+4P得到6P

可以看看下面的Python代码实现

def bits(n):
    '''
    bits(151) => 1, 1, 1, 0, 1, 0, 0, 1
    '''
    while n:
        yield n & 1
        n >>= 1

def double_and_add(n, x):
    result = 0
    addend = x

    for bit in bits(n):
        if bit == 1:
            result += addend
        addend *= 2

    return result

参考：
https://www.cnblogs.com/Kalafinaian/p/7392505.html
https://oi.men.ci/mul-inverse/
https://eng.paxos.com/blockchain-101-foundational-math
https://eng.paxos.com/blockchain-101-elliptic-curve-cryptography
https://andrea.corbellini.name/2015/05/23/elliptic-curve-cryptography-finite-fields-and-discrete-logarithms/