算法_1： 数论

除法

• 记：a | b 表示 a 整除 b
• q=a div d； r=a mod d —– a=dq+r
  101 = 11 * 9+2————— 9=101 div 11； 2=101 mod 11

模算术

• m整除 a-b 时 称 a 模 m 同余 b。

• 记： a $\equiv$ b (mod m)——- 同余式
  17是否模6同余5———17 $\equiv$ 5(mod 6) : 6 整除 17-5

• a $\equiv$ b (mod m) 当且仅当 a mod m=b mod m

模m算术

• $a+_mb = (a+b) \ mod \ m$
• $a \cdot_m b = (a \cdot b) \ mod \ m$

---

整数的表示和算法

• 整数的表示
  $n=a_kb^k+a_{k-1}b^{k-1}+......+a_1b+a_0$

• 进制转换

241=2*120 +1

120=2*60 + 0
60 =2*30 + 0
30= 2*15 + 0
15= 2*7  + 0
7=  2*3  + 1
3=  2*1  + 1
1=  2*0  + 1

(241)10=(11100001)

---

算法 : 构造b进制展开式

$procedure \ base \ b \ expansion$
$q := n$
$k := 0$

$while \ \ q \neq 0$
$\\ a_k:=1 \ mod \ b$
$\\ q:= q \ div \ b$
$\\ k:= k + 1$

$return (a_{k-1}......a_1a_0)$

---

整数加法运算

• 原理：
  $a=(a_{n-1}....a_1a_0)_2 和 \ b=(b_{n-1}.....b_1b_0)$

  $a_0+b_0=c_0*2+s_0$——$c_0是进位$
  $a_1+b_1+c_0=c_1*2+s_0$
  $s_n=c_n-1$是首位

• 例子 ： a=(1110) 和 b=(1011)

  $a_0+b_0=0+1=0*2+1$—–$c_0=0,s_0=1$
  $a_1+b_1+c_0=1+1+0=1*2+0$——$c_1=1,s_1=0$
  $a_2+b_2+c_1=1+0+1=1*2+0$—–$c_2=1,s_2=0$
  $a_3+b_3+c_2=1+1+1=1*2+1$—–$c_3=1,s3=1$
  $s_4=c_3=1$

a+b=(11001)

• 算法 : procedure add

  ${a和b分别是整数展开式：(a_{n-1}....a_1a_0) 和 (b_{n-1}.....b_1b_0)}$
  $c:=0$

  $for \ j:=0 \ to \ n-1$
  {
  $\\ d:=\lfloor(a_j+b_j+c)/2\rfloor$
  $\\ s_j := a_j+b_j+c -2d$
  $c:=d$
  $s_n:=c$
  }
  return ($s_ns_{n-1}......s_0$)

• 复杂度=O(n)

---

整数乘法运算

• 原理：
  $ab = a(b_02^0+b_12^1+......+b_{n-1}2^{n-1})$
  $=a(b_02^0)+a(b_12^1)+.....+a(b_{n-1}2^{n-1})$

• 例子：a=(110) 和 b=(101)

$ab_0*2^0=(110)_2*1*2^0=(110)_2$
$ab_1*2^1=(110)_2*0*2^1=(0000)_2$
$ab_2*2^2=(110)_2*1*2^2=(11000)_2$

相加$(110)+(0000)+(11000)=(11110)$

• 算法

for j:= 0 to n-1
     if b_j=1 then c_j := a 移位j位
     else c_j:=0
p:=0
for j:=0 to n-1
    p:=p+c_j

return p {p=ab}

• 复杂度： $O(n^2)$

---

模指数运算

• 原理：
  $n=(a_{k-1}.....a_1a_0)$ 计算 $b^n$
  $b^n=b^{a_{k-1}2^{k-1}+.....+a_1*2+a_0}=b^{a_{k-1}2^{k-1}}.....b^{a_1*2}b^{a_0}$

• 例子：$3^{644}$ mod 645

644=(1010000100)
x=1 and power=3 mod 645=3
i=0: 因为 $a_0=0$, 所以 x=1 和 power=$3^2$ mod 645= 9 mod 645=9;
i=1: 因为 $a_1=0$, 所以 x=1 和 power=$9^2$ mod 645= 81 mod 645=81;
i=2: 因为 $a_2=1$, 所以 x=1*81 mod 645=81 and power=$81^2$ mod 645= 6561 mod 645 =111;
i=3: 因为 $a_3=0$, 所以 x=81 and power=$111^2$ mod 645=12321 mod 645=66;
i=4: 因为 $a_4=0$, 所以 x=81 and power=$66^2$ mod 645=4356 mod 645=486;
i=5: 因为 $a_5=0$, 所以 x=81 and power=$486^2$ mod 645=236196 mod 645=126;
i=6: 因为 $a_6=0$, 所以 x=81 and power=$126^2$ mod 645=15876 mod 645=396;
i=7: 因为 $a_7=1$, 所以 x=(81*396) mod 645=471 和 power=$396^2$ mod 645=1568176 mod 645=81;
i=8: 因为 $a_8=0$, 所以 471 和 power=$81^2$ mod 645=6561 mod 645=111;
i=9: 因为 $a_9=1$, 所以有x=(471*111) mod 645=36

• 算法：

x：=1
power := b mod m
for i := 0 to k-1

 if $a_i$=1 then x:=(x*power) mod m
 power :=(power*power) mod m

return x {x等于 $b^n$ mod m}

• 复杂度：$O((logm)^2logn)$

---

---

---

素数

• 素数： 每个大于1的整数，只能被1或者本身整除，称为素数

*试除法

• n是一个合数，则n必有一个素因子小于或者等于$\sqrt{n}$
• 若整数不能被小于或者等于其平方根的素数整除，则就是素数
  101—– $\sqrt{101}的素数是 2,3,5,7$；因为101不能被2,3,5,7整除，故101是素数

*埃拉托斯特尼筛选

*素数定理

• 当x无限增长的时候，不超过$x$的 素数个数与$x/lnx$之比趋于1

---

最大公约数

• 能整除两个整数的最大整数为这两个整数的最大公约数
• 记： gcd ( a, b )
• 两个数互素则最大公约数是1

*素因数分解法求最大公约数

• 正整数： a 和 b的素因子分解是

• $a=p_1^{a_1}p_2^{a_2}.....p_n^{a_n}$ 和 $b=p_1^{b_1}p_2^{b_2}.....p_n^{b_n}$

• gcd（a,b）=$p_1^{min(a_1,b_1)}p_2^{min(a_2,b_2)}......p_n^{min(a_n,b_n)}$

*素因数分解法求最小公倍数

• lcm=$p_1^{max(a_1,b_1)}p_2^{max(a_2,b_2)}......p_n^{max(a_n,b_n)}$

*欧几里得算法

• 算法

$a=bq+r$, 其中 a, b, q 和 r 均为整数，则 gcd(a,b)=gcd(b,r)

• 例子

662 = 414*1+248
414 = 248*1+166
248 = 166*1+ 82
166= 82*2 +2
82=2*41

所以 $gcd(414, 662)=2$

• 算法

x:=a
y:=b

while y 不等于 0
r:= x mod y
x:= y
y:= r
return x （gcd(a ,b)是x）

---

gcd的线性组合

• 贝祖定理： $gcd (a, b)=sa+tb$
• 例：gcd(6,14)=2 ,而 2=(-2)*6+ 1*14

---

---

---

求解同余方程

*线性同余方程

• $ax \equiv b (mod \ m)$

• 如果 a 和 m 为互素的整数，且m>1，则 a模m的逆存在

• $a \overline {a} \equiv 1 (mod \ m)$ 成立的 $\overline{a}$ 是逆

*中国剩余定理

• 令$m_1, m_2, .... , m_n$ 为大于1的两两互素的正整数， $a_1,a_2,....a_n$ 是任意数。则
  $x \equiv a_1 ( mod \ m_1)$
  $x \equiv a_2 ( mod \ m_2)$
  ……
  $x \equiv a_n ( mod \ m_n)$

有唯一的模 $m=m_1m_2.....m_n$的解

*费马小定理

• 如果p是素数，a是一个不能被p整除的整数，则

$a^{p-1} \equiv 1 (mod \ p)$ ,且 对于每个 a有， $a^p \equiv a(mod \ p)$

• 例子：计算$7^{222}$ mod 11

$7^{10} \equiv$ 1(mod 11)
所以对于每个正整数k, 有
$(7^{10})^k \equiv 1(mod \ 11)$
222=22*10+2

—>$7^{222}=7^{22*10+2} = (7^{10})^{22}7^2 \equiv (1)^{22}*49 \equiv 5(mod \ 11)$

---

---

---

同余应用

*散列函数

*伪随机数

• 线性同余法：模数 m, 倍数 a, 增量 c 和种子 $x_0$；
  $2 <=a <m$
  $0 <=c <m$
  $0 <=x_0 <m$
  生成伪随机序列
  $x_{n+1}=(ax_n+c) mod m$

• 例子： $x_{n+1}=(7x_n+4) mod \ 9$

$x_1=7x_0+4 \ mod \ 9 =7*3+4 \ mod \ 9=25 \ mod \ 9=7$
$x_2=7x_1+4 \ mod \ 9 =7*7+4 \ mod \ 9=53 \ mod \ 9=8$
……

• c=0 则为纯倍增生成器

*校验码