算法导论-第31章-数论算法

本章将处理大整数。“大输入”通常指包含“大整数”的输入，而不是包含“很多整数“的输入。因此，我们利用输入所需的位数来度量输入的大小，而不仅仅是输入中整数的数目。

给定 $k$ 个整数输入 $a_1,a_2,\cdots ,a_n$，如果算法可以在关于 $\log a_1,\log a_2,\cdots ,\log a_n$ 的多项式时间内完成，即算法在关于二进制编码后的输入长度的多项式时间内完成，则称该算法为多项式时间算法。 （不是很懂🤣）

31.1 基础数论概念

整数集 Z = { ⋯   , − 2 , − 1 , 0 , 1 , 2 , ⋯   } Z=\{\cdots, -2, -1, 0, 1, 2, \cdots\} Z={⋯,−2,−1,0,1,2,⋯}。

自然数集 N = { 0 , 1 , 2 , ⋯   } N=\{0, 1, 2, \cdots\} N={0,1,2,⋯}。

整除性与约数

**符号 $d|a$ （读作” $d$ 整除 $a$“）的含义是，存在某个整数 $k$，使得 $\ast \ast a=kd$。

任何整数均可整除0。

如果 $d|a$ 且 $d\ge 0$，则称 $d$ 是 $a$ 的约数。

不失一般性，规定约数为非负数，非零整数 $a$ 的约数至少为1，且不大于 $|a|$。例如，24的约数是1，2，3，4，6，8，12，24。

任何正整数 $a$ 均可被平凡数1和其自身 $a$ 整除。整数 $a$ 的非平凡约数称为 $a$ 的因子。例如，20的因子是2， 4， 5， 10。

素数与合数

如果一个整数 $a>1$ 且只能被平凡约数1和它自身所整除，则这个数是素数。

前20个素数排列如下：2， 3， 5， 7， 11， 13， 17， 19， 23， 29， 31， 37， 41， 43， 47， 53， 59， 61， 67， 71.

如果一个整数 $a>1$ 且不是素数，则这个数是合数。

0， 1和所有负整数都既不是素数也不是合数。

除法定理、余数和模

除法定理：对于任何整数 $a$ 和任何正整数 $n$，存在唯一整数 $q$ 和 $r$，满足 $0\le r<n$ 且 $a=qn+r$。

称 $q=\lfloor a/n\rfloor$ 为除法的商，值 $r=amodn$ 为除法的余数。$n|a$ 当且仅当 $amodn=0$。

模 $n$ 的等价类： [ a ] n = { a + k n ∣ k ∈ Z } [a]_n=\{a+kn|k \in Z\} [a]n​={a+kn∣k∈Z}

• $a\in [b{]}_n$，则$a\equiv b(modn)$。

公约数与最大公约数

$a,b$ 不全为零的两个整数，则最大公约数 $gcd(a,b)$ 是 { a x + b y ∣ x , y ∈ Z } \{ax+by|x,y \in Z\} {ax+by∣x,y∈Z} 中最小的正整数。

• 对所有整数 $a$ 和 $b$，如果 $d|a$，$d|b$，则 $d|gcd(a,b)$。
• 对所有整数 $a$ 和 $b$，和非负整数 $n$，有 $gcd(an,bn)=ngcd(a,b)$。
• 对所有正整数 $n$、$a$ 和 $b$，如果 $n|ab$ 且 $gcd(a,n)=1$，则 $n|b$。

互素数和素数唯一性分解定理

如果 $gcd(a,b)=1$，称 $a$ 和 $b$ 为互素数。

$\forall a,b,p\in Z$，如果 $gcd(a,p)=1$ 和 $gcd(b,p)=1$，则 $gcd(ab,p)=1$。

唯一分解定理：一个合数能被唯一写成形式 $a=p_1^{e_1}{p}_2^{e_2}\cdots p_r^{e_r}$，这里 $p_i$ 是素数，$p_1<p_2<\cdots <p_r$，$e_i$ 是正整数。例如，$6000=2^4\times 3\times 5^3$

31.2 最大公约数

一种直观解GCD

$a=p_1^{e_1}{p}_2^{e_2}\cdots p_r^{e_r}$ 和 $b=p_1^{f_1}{p}_2^{f_2}\cdots p_r^{f_r}$ $\Rightarrow$ g c d ( a , b ) = p 1 m i n { e 1 , f 1 } p 2 m i n { e 2 , f 2 } ⋯ p r m i n { e r , f r } gcd(a,b) = p_{1}^{min\{e_1,f_1\}} p_{2}^{min\{e_2,f_2\}} \cdots p_{r}^{min\{e_r,f_r\}} gcd(a,b)=p1min{e1​,f1​}​p2min{e2​,f2​}​⋯prmin{er​,fr​}​。

注：这种解法需要整数的素因子分解，而素因子分解是一个很难的问题（NP-hard）。

Euclid‘s Algorithm

EUCLID(a, b)
    if b == 0
        return a
    else 
        return EUCLID(b, a mod b)

以 $gcd(30,21)$ 为例，EUCLID(30, 21) = EUCLID(21, 9) = EUCLID(9, 3) = EUCLID(3, 0) = 3

Euclid‘s Algorithm的递归调用次数为 $\mathrm{O}(\log b)$。

Extended-Euclid‘s Algorithm

31.3 求解线性模方程

求 $x$ 满足下面的模方程 $ax\equiv b(modn)$，其中 $a>0,n>0$。

当且仅当 $gcd(a,n)|b$ 时，方程 $ax\equiv b(modn)$ 有解。

例题，判断 $4x\equiv 2(mod6)$ 和 $4x\equiv 3(mod6)$ 有无解。

• 因为 $gcd(4,6)|2$，所以 $4x\equiv 2(mod6)$ 有解。
• 因为 $gcd(4,6)不整除3$，所以 $4x\equiv 3(mod6)$ 无解。

设 $d=gcd(a,n)$ 且 $d=a{x}^{\prime }+n{y}^{\prime }$，对某个整数 $x^{\prime }$ 和 $y^{\prime }$。如果 $d|b$，则有特解 $x_0=x^{\prime }(b/d)modn$。

设 $x_0$ 为 $ax\equiv b(modn)$ 的一个解，则有 $d$ 个不同的解：$x_i=x_0+i(n/d),i=0,1,2,\cdots ,d-1$。

方程 $ax\equiv b(modn)$ 要么对模 $n$ 有 $gcd(a,n)$ 个不同的解，要么无解。

例题，求 $14x\equiv 30(mod100)$。

解：调用$ExtendedEuclid(14,100)=(d,x^{\prime },y^{\prime })=(2,-7,1)$

∵2|30

∴特解$x_0=-7\times 30/2(mod100)=95$

调用求全部解算法：二个解95,45

$x_0=95+0\times 100/2=95$

$x_1=95+1\times 100/2=145=45(mod100)$

31.4 中国剩余定理

模逆: $a{a}^{-1}\equiv 1(modn)$。

定义 $m_i=n/n_i$，于是可得 $m_i=n_1{n}_2\cdots n_{i-1}{n}_{i+1}\cdots n_k$。接着定义 $c_i=m_i(m_i^{-1}mod{n}_i)$。
$$a\equiv (a_1{c}_1+a_2{c}_2+\cdots +a_k{c}_k)(modn)$$
例题，已知 $a\equiv 2(mod5),a\equiv 3(mod13)$。

解：
$$a_1=2,a_2=3$$

$$n_1=5,n_2=13$$

$$m_1=n_2=13,m_2=n_1=5$$

$$m_1{m}_1^{-1}=1(mod5),m_1^{-1}=2$$

$$m_2{m}_2^{-1}=1(mod13),m_2^{-1}=8$$

$$c_1=13(2mod5)=26$$

$$c_2=5(8mod13)=40$$

$$\begin{gathered} a\equiv 2\cdot 26+3\cdot 40(mod65) \\ \equiv (52+120)(mod65) \\ \equiv 42(mod65) \end{gathered}$$

例题，已知 $x\equiv 4(mod5)$ 和 $x\equiv 5(mod11)$。

解：
$$a_1=4,a_2=5$$

$$n_1=5,n_2=11$$

$$m_1=n_2=11$$

$$m_2=n_1=5$$

$$\begin{gathered} c_1=m_1(m_1^{-1}mod{n}_1) \\ =11(11^{-1}mod5),其中11\cdot 11^{-1}=1(mod5),得11^{-1}=1 \\ =11(1mod5) \\ =11 \end{gathered}$$

$$\begin{gathered} c_2=m_2(m_2^{-1}mod{n}_2) \\ =5(5^{-1}mod11),其中5\cdot 5^{-1}=1(mod11),得5^{-1}=9 \\ =5(9mod11) \\ =45 \end{gathered}$$

$$\begin{gathered} a=(a_1{c}_1+a_2{c}_2)(modn),其中n=n_1{n}_2=5\times 11=55 \\ =(4\times 11+5\times 45)(mod55) \\ =269mod55 \\ =49 \end{gathered}$$

关键在于求模逆。其他套公式。

31.5 RSA

构造过程：

1. 随机选择两个大素数 $p$ 和 $q$， $p<>q$；
2. 计算 $n=pq$；
3. 选择小的奇整数 $e$，使得 $e$ 与 $(p-1)(q-1)$ 互素；
4. 计算 $e$ 模 $(p-1)(q-1)$ 的逆 $d$，即 $ed=1(mod(p-1)(q-1))$；
5. 公布 $e,n$为RSA的公钥；
6. 保存 $d,p,q$ 为RSA的私钥。

加解密过程：

• 加密：$P(M)=M^{e}modn$; (=C)
• 解密：$S(C)=C^{d}modn$;

RSA（Rivest-Shamir-Adleman）是一种非对称加密算法，其正确性基础基于数论中两个关键性质：质因数分解困难性和模幂运算的可逆性。

1. 质因数分解困难性：RSA的安全性依赖于大整数的质因数分解问题的困难性。给定一个大的合数，将其分解为质数的乘积是一个计算上非常困难的问题，特别是当质数很大时。RSA利用了这一困难性，通过使用大素数进行加密和解密，使得破解者难以找到私钥（包含了质因数）。

2. 模幂运算的可逆性：RSA算法中使用了模幂运算，即对一个数进行指数幂运算后再取模的操作。模幂运算具有可逆性，即给定模数、底数和指数，可以通过模幂运算计算出唯一的结果。在RSA中，公钥加密和私钥解密过程都涉及到模幂运算，确保了正确的加密和解密操作。

31.6 素数判定

简单素数测试算法

伪素数测试算法

31.7 MR算法的改进措施和算法复杂性

• 措施1：随机选择若干个基，以减少计算量
• 措施2：在计算 $a^{n-1}modn$ 过程中，加入合数证人（Witness）的查找

设 $n$ 为 $\beta -bit$ 整数，算法需要 $\mathrm{O}(s\beta )$ 算数运算，$\mathrm{O}({\beta }^3)$ 位运算。