算法设计与分析复习笔记
第一部分:算法基础与 RSA 加密
RSA 加密算法
1. 加密过程
步骤 说明 描述 1 选择一对 不相等 且 足够大 的质数 p , q p, q p , q 2 计算 p , q p, q p , q 的乘积 n = p × q n = p \times q n = p × q 3 计算 n n n 的欧拉函数 φ ( n ) = ( p − 1 ) × ( q − 1 ) \varphi(n) = (p-1) \times (q-1) φ ( n ) = ( p − 1 ) × ( q − 1 ) 4 选一个与 φ ( n ) \varphi(n) φ ( n ) 互质的整数 e e e 1 < e < φ ( n ) 1 < e < \varphi(n) 1 < e < φ ( n ) 5 计算出 e e e 对于 φ ( n ) \varphi(n) φ ( n ) 的模反元素 d d d d e m o d φ ( n ) = 1 de \bmod \varphi(n) = 1 d e mod φ ( n ) = 1 6 公钥 K U = ( e , n ) KU = (e, n) K U = ( e , n ) 7 私钥 K R = ( d , n ) KR = (d, n) K R = ( d , n )
明文 M M M → \quad \rightarrow \quad → 加密 M e m o d n = C M^e \bmod n = C M e mod n = C
密文 C C C → \quad \rightarrow \quad → 解密 C d m o d n = M C^d \bmod n = M C d mod n = M
计算欧拉函数
欧拉函数 是小于 n 的正整数中与 n 互质的数的数目。
互质 是公约数只有 1 的两个整数,叫做互质整数。
质数 是指在大于 1 的自然数中,除了 1 和它本身以外不再有其他因数的自然数。φ ( n ) = n − 1 \varphi(n)=n-1 φ ( n ) = n − 1
如果 n n n 可以分解成 2 个 互质 的整数之积,那么 n n n 的欧拉函数等于这两个因子的欧拉函数之积。
即若 n = p × q n = p \times q n = p × q ,且 p , q p, q p , q 互质,则 φ ( n ) = φ ( p × q ) = φ ( p ) × φ ( q ) = ( p − 1 ) × ( q − 1 ) \varphi(n) = \varphi(p \times q) = \varphi(p) \times \varphi(q) = (p-1) \times (q-1) φ ( n ) = φ ( p × q ) = φ ( p ) × φ ( q ) = ( p − 1 ) × ( q − 1 ) 。
计算模反元素d d d
根据定义,模反元素可以用以下三种等价的数学形式表示:
整除形式:
e d − 1 = k φ ( n ) ed - 1 = k\varphi(n) e d − 1 = k φ ( n )
取模运算形式:
e d m o d φ ( n ) = 1 ed \bmod \varphi(n) = 1 e d mod φ ( n ) = 1
同余式形式:
e d ≡ 1 ( m o d φ ( n ) ) ed \equiv 1 \pmod{\varphi(n)} e d ≡ 1 ( mod φ ( n ))
3 和 11,3 × 4 − 1 = 11 3\times4-1=11 3 × 4 − 1 = 11 ,那么 d = 4 d=4 d = 4 ,那么 4 ± k ⋅ 11 4±k · 11 4 ± k ⋅ 11 都可以是 d,也就是 d + k φ ( n ) d+k\varphi(n) d + k φ ( n ) 都是 e 的模反元素。d 不唯一,一般取最小的那个便于计算。
证明
明文 M M M → \quad \rightarrow \quad → 加密 M e m o d n = C M^e \bmod n = C M e mod n = C
密文 C C C → \quad \rightarrow \quad → 解密 C d m o d n = M C^d \bmod n = M C d mod n = M
不考,不学了
2. 例题
第一步:密钥生成 (Key Generation)
确定参数:
已知:p = 5 p = 5 p = 5 , q = 11 q = 11 q = 11 , e = 3 e = 3 e = 3 。
计算公共模数 n n n :
n = p × q = 5 × 11 = 55 n = p \times q = 5 \times 11 = 55 n = p × q = 5 × 11 = 55
计算欧拉函数 ϕ ( n ) \phi(n) ϕ ( n ) :
ϕ ( n ) = ( p − 1 ) ( q − 1 ) = ( 5 − 1 ) ( 11 − 1 ) = 4 × 10 = 40 \phi(n) = (p-1)(q-1) = (5-1)(11-1) = 4 \times 10 = 40 ϕ ( n ) = ( p − 1 ) ( q − 1 ) = ( 5 − 1 ) ( 11 − 1 ) = 4 × 10 = 40
计算私钥 d d d :
d d d 是 e e e 关于 ϕ ( n ) \phi(n) ϕ ( n ) 的模反元素,需满足:( d × e ) ≡ 1 ( m o d ϕ ( n ) ) (d \times e) \equiv 1 \pmod{\phi(n)} ( d × e ) ≡ 1 ( mod ϕ ( n )) 。
即:( d × 3 ) ≡ 1 ( m o d 40 ) (d \times 3) \equiv 1 \pmod{40} ( d × 3 ) ≡ 1 ( mod 40 ) 。
结果:
公钥 (Public Key): ( n = 55 , e = 3 ) (n=55, e=3) ( n = 55 , e = 3 )
私钥 (Private Key): ( n = 55 , d = 27 ) (n=55, d=27) ( n = 55 , d = 27 )
第二步:加密过程 (Encryption)
已知明文 m = 14 m = 14 m = 14 ,加密公式为:c = m e ( m o d n ) c = m^e \pmod n c = m e ( mod n ) 。
代入数据:c = 14 3 ( m o d 55 ) c = 14^3 \pmod{55} c = 1 4 3 ( mod 55 )
计算:14 3 = 2744 14^3 = 2744 1 4 3 = 2744
取模:2744 ÷ 55 = 49 2744 \div 55 = 49 2744 ÷ 55 = 49 余 49 49 49
(55 × 49 = 2695 55 \times 49 = 2695 55 × 49 = 2695 )
2744 − 2695 = 49 2744 - 2695 = 49 2744 − 2695 = 49
加密结果: 密文 c = 49 c = 49 c = 49 。
第三步:解密过程 (Decryption)
使用私钥对密文 c = 49 c = 49 c = 49 进行解密,公式为:m = c d ( m o d n ) m = c^d \pmod n m = c d ( mod n ) 。
代入数据:m = 49 27 ( m o d 55 ) m = 49^{27} \pmod{55} m = 4 9 27 ( mod 55 )
简化计算 (使用快速幂取模法):
49 ≡ − 6 ( m o d 55 ) 49 \equiv -6 \pmod{55} 49 ≡ − 6 ( mod 55 )
49 2 ≡ ( − 6 ) 2 = 36 ≡ − 19 ( m o d 55 ) 49^2 \equiv (-6)^2 = 36 \equiv -19 \pmod{55} 4 9 2 ≡ ( − 6 ) 2 = 36 ≡ − 19 ( mod 55 )
49 3 ≡ 36 × ( − 6 ) = − 216 49^3 \equiv 36 \times (-6) = -216 4 9 3 ≡ 36 × ( − 6 ) = − 216 。− 216 ÷ 55 -216 \div 55 − 216 ÷ 55 的余数为 4 4 4 (因为 55 × ( − 4 ) = − 220 55 \times (-4) = -220 55 × ( − 4 ) = − 220 )。
通过幂运算最后可得出:49 27 ≡ 14 ( m o d 55 ) 49^{27} \equiv 14 \pmod{55} 4 9 27 ≡ 14 ( mod 55 ) 。
解密结果: 明文 m = 14 m = 14 m = 14 (与原始明文一致,验证成功)。
总结
中间参数: n = 55 , ϕ ( n ) = 40 , d = 27 n = 55, \phi(n) = 40, d = 27 n = 55 , ϕ ( n ) = 40 , d = 27
加密结果: 密文为 49
解密结果: 明文为 14
注意余数的缩放
你对被除数和除数同时约分(同时除以 5)时,余数也随之缩小了 5 倍 。
在模运算中有一个重要的性质:
如果 a ÷ b = q … r a \div b = q \dots r a ÷ b = q … r (即 a = b q + r a = bq + r a = b q + r )
那么当我们将 a a a 和 b b b 同时除以 k k k 时:
( a / k ) = ( b / k ) q + ( r / k ) (a/k) = (b/k)q + (r/k) ( a / k ) = ( b / k ) q + ( r / k )
在你这里:
原始算式:100000 = 35 × 2857 + R 100000 = 35 \times 2857 + \mathbf{R} 100000 = 35 × 2857 + R
约分后的算式:20000 = 7 × 2857 + 1 20000 = 7 \times 2857 + \mathbf{1} 20000 = 7 × 2857 + 1
要还原回原始算式的余数 R R R ,你需要把约分掉的 5 乘回去:
R = 1 × 5 = 5 R = 1 \times 5 = 5 R = 1 × 5 = 5
你的计算结果: 余数 1 是 20000 ÷ 7 20000 \div 7 20000 ÷ 7 的结果。
实际要求的解: 是 100000 ÷ 35 100000 \div 35 100000 ÷ 35 的余数。
修正方法: 约分计算完余数后,必须乘以约掉的那个公约数 。
所以:1 × 5 = 5 1 \times 5 = 5 1 × 5 = 5 ,这也就是为什么我在上一个回答中算出的明文结果是 5 。
💡 小贴士:如何避免这种错误?
在处理 RSA 这种大数取模时,建议使用**“取模性质”**分步计算,而不是先把大数算出来再除,这样不容易出错:
10 2 = 100 ≡ 30 ( m o d 35 ) 10^2 = 100 \equiv 30 \pmod{35} 1 0 2 = 100 ≡ 30 ( mod 35 )
10 4 = 30 × 30 = 900 10^4 = 30 \times 30 = 900 1 0 4 = 30 × 30 = 900
900 ÷ 35 = 25 … 25 900 \div 35 = 25 \dots 25 900 ÷ 35 = 25 … 25 (即 10 4 ≡ 25 ( m o d 35 ) 10^4 \equiv 25 \pmod{35} 1 0 4 ≡ 25 ( mod 35 ) )
10 5 = 10 4 × 10 = 250 10^5 = 10^4 \times 10 = 250 1 0 5 = 1 0 4 × 10 = 250
250 ÷ 35 = 7 … 5 250 \div 35 = 7 \dots \mathbf{5} 250 ÷ 35 = 7 … 5
这种方法数字更小,不容易在约分上栽跟头。
(A × B A \times B A × B )对 N N N 取模,等于(A A A 对 N N N 取模的结果 × \times × B B B 对 N N N 取模的结果)再对 N N N 取模。
下面我把每一步拆解得更详细一些:
10 2 = 100 ≡ 30 ( m o d 35 ) 10^2 = 100 \equiv 30 \pmod{35} 1 0 2 = 100 ≡ 30 ( mod 35 )
计算: 10 × 10 = 100 10 \times 10 = 100 10 × 10 = 100 。
取模: 100 ÷ 35 = 2 … 100 \div 35 = 2 \dots 100 ÷ 35 = 2 … 余 30 30 30 。
意义: 在接下来的计算中,凡是遇到 10 2 10^2 1 0 2 的地方,我们都可以用它的余数 30 30 30 来代替。
10 4 = 30 × 30 = 900 10^4 = 30 \times 30 = 900 1 0 4 = 30 × 30 = 900
计算: 10 4 10^4 1 0 4 其实就是 10 2 × 10 2 10^2 \times 10^2 1 0 2 × 1 0 2 。
替换: 根据取模原理,我们直接用刚才算出的余数代替:30 × 30 = 900 30 \times 30 = 900 30 × 30 = 900 。
取模: 900 ÷ 35 = 25 … 900 \div 35 = 25 \dots 900 ÷ 35 = 25 … 余 25 25 25 。
验证:35 × 25 = 875 35 \times 25 = 875 35 × 25 = 875 ,900 − 875 = 25 900 - 875 = 25 900 − 875 = 25 。
意义: 现在我们知道了 10 4 ≡ 25 ( m o d 35 ) 10^4 \equiv 25 \pmod{35} 1 0 4 ≡ 25 ( mod 35 ) 。
10 5 = 10 4 × 10 = 250 10^5 = 10^4 \times 10 = 250 1 0 5 = 1 0 4 × 10 = 250
计算: 题目要求的是 10 5 10^5 1 0 5 ,它等于 10 4 × 10 1 10^4 \times 10^1 1 0 4 × 1 0 1 。
替换: 我们用刚才算出的 10 4 10^4 1 0 4 的余数 25 25 25 来代替它:25 × 10 = 250 25 \times 10 = 250 25 × 10 = 250 。
取模: 250 ÷ 35 = 7 … 250 \div 35 = 7 \dots 250 ÷ 35 = 7 … 余 5 5 5 。
验证:35 × 7 = 245 35 \times 7 = 245 35 × 7 = 245 ,250 − 245 = 5 250 - 245 = 5 250 − 245 = 5 。
通过这种方法,你处理的最大数字只是 900 900 900 和 250 250 250 ,而不是 100 , 000 100,000 100 , 000 。
如果你觉得 30 30 30 还是太大,甚至可以利用 “负余数” :
100 ÷ 35 100 \div 35 100 ÷ 35 也可以看作是:差 5 5 5 就到 105 105 105 (35 × 3 35 \times 3 35 × 3 )了。
所以 100 ≡ − 5 ( m o d 35 ) 100 \equiv -5 \pmod{35} 100 ≡ − 5 ( mod 35 ) 。
那么 10 4 = ( 10 2 ) 2 ≡ ( − 5 ) 2 = 25 ( m o d 35 ) 10^4 = (10^2)^2 \equiv (-5)^2 = \mathbf{25} \pmod{35} 1 0 4 = ( 1 0 2 ) 2 ≡ ( − 5 ) 2 = 25 ( mod 35 ) 。
最后 10 5 = 25 × 10 = 250 ≡ 5 ( m o d 35 ) 10^5 = 25 \times 10 = 250 \equiv \mathbf{5} \pmod{35} 1 0 5 = 25 × 10 = 250 ≡ 5 ( mod 35 ) 。
这样计算量会进一步减小。这种分段处理余数的方法在 RSA 算法(尤其是处理大位数密钥时)是必须使用的,被称为 “模幂运算” 。
3. 扩展欧几里得求乘法逆元
用推广的 Euclid 算法求 67mod119 的逆元。
我们要找一个整数 (x),使得
67 x ≡ 1 ( m o d 119 ) 67x \equiv 1 \pmod{119} 67 x ≡ 1 ( mod 119 )
等价于存在整数 (k),满足
67 x + 119 k = 1 67x + 119k = 1 67 x + 119 k = 1
👉 这就是一个典型的线性丢番图方程
67 x + 119 y = 1 67x + 119y = 1 67 x + 119 y = 1
只要 gcd ( 67 , 119 ) = 1 \gcd(67,119)=1 g cd( 67 , 119 ) = 1 ,这个方程就有整数解,而 推广的 Euclid 算法 正是用来:
第一步:辗转相除法(求最大公约数)
我们先用较大的数除以较小的数,直到余数为 0:
119 = 1 × 67 + 52 119 = 1 \times 67 + 52 119 = 1 × 67 + 52
67 = 1 × 52 + 15 67 = 1 \times 52 + 15 67 = 1 × 52 + 15
52 = 3 × 15 + 7 52 = 3 \times 15 + 7 52 = 3 × 15 + 7
15 = 2 × 7 + 1 15 = 2 \times 7 + 1 15 = 2 × 7 + 1 (出现 1,说明 67 67 67 和 119 119 119 互质,逆元存在)
7 = 7 × 1 + 0 7 = 7 \times 1 + 0 7 = 7 × 1 + 0
最后一个非零余数是 1
所以
gcd ( 67 , 119 ) = 1 \gcd(67,119)=1 g cd( 67 , 119 ) = 1
✅ 逆元存在
⚠️ 很重要的一点:
不是所有模数都有逆元,只有当 gcd ( a , m ) = 1 \gcd(a,m)=1 g cd( a , m ) = 1 时,(a) 才在模 (m) 下可逆
第二步:回代法(求逆元)
把“1”一步一步写成“67 和 119 的线性组合”
也就是最终要得到:
1 = ( 某个数 ) ⋅ 67 + ( 某个数 ) ⋅ 119 1 = (\text{某个数})\cdot 67 + (\text{某个数})\cdot 119 1 = ( 某个数 ) ⋅ 67 + ( 某个数 ) ⋅ 119
其中 67 前面的系数 ,就是我们要的逆元(模 119 意义下)。
现在我们从余数为 1 的等式开始,反向代入上述等式:
由第 (4) 式得:1 = 15 − 2 × 7 1 = 15 - 2 \times 7 1 = 15 − 2 × 7
代入第 (3) 式中的 7 = 52 − 3 × 15 7 = 52 - 3 \times 15 7 = 52 − 3 × 15 :
1 = 15 − 2 × ( 52 − 3 × 15 ) 1 = 15 - 2 \times (52 - 3 \times 15) 1 = 15 − 2 × ( 52 − 3 × 15 )
1 = 15 − 2 × 52 + 6 × 15 1 = 15 - 2 \times 52 + 6 \times 15 1 = 15 − 2 × 52 + 6 × 15
1 = 7 × 15 − 2 × 52 1 = 7 \times 15 - 2 \times 52 1 = 7 × 15 − 2 × 52
代入第 (2) 式中的 15 = 67 − 1 × 52 15 = 67 - 1 \times 52 15 = 67 − 1 × 52 :
1 = 7 × ( 67 − 52 ) − 2 × 52 1 = 7 \times (67 - 52) - 2 \times 52 1 = 7 × ( 67 − 52 ) − 2 × 52
1 = 7 × 67 − 7 × 52 − 2 × 52 1 = 7 \times 67 - 7 \times 52 - 2 \times 52 1 = 7 × 67 − 7 × 52 − 2 × 52
1 = 7 × 67 − 9 × 52 1 = 7 \times 67 - 9 \times 52 1 = 7 × 67 − 9 × 52
代入第 (1) 式中的 52 = 119 − 1 × 67 52 = 119 - 1 \times 67 52 = 119 − 1 × 67 :
1 = 7 × 67 − 9 × ( 119 − 67 ) 1 = 7 \times 67 - 9 \times (119 - 67) 1 = 7 × 67 − 9 × ( 119 − 67 )
1 = 7 × 67 − 9 × 119 + 9 × 67 1 = 7 \times 67 - 9 \times 119 + 9 \times 67 1 = 7 × 67 − 9 × 119 + 9 × 67
1 = 16 × 67 − 9 × 119 1 = 16 \times 67 - 9 \times 119 1 = 16 × 67 − 9 × 119
第三步:得出结果
根据最后一个等式 1 = 16 × 67 − 9 × 119 1 = 16 \times 67 - 9 \times 119 1 = 16 × 67 − 9 × 119 ,在模 119 119 119 的环境下:
结论:
67 ( m o d 119 ) 67 \pmod{119} 67 ( mod 119 ) 的逆元是 16。
方法二:矩阵表格法(求逆元)
矩阵表格法 = 把回代过程提前,用递推自动完成 。这种方法通过建立一个表格,省去了回代的麻烦。我们还是求 67 ( m o d 119 ) 67 \pmod{119} 67 ( mod 119 ) 的逆元。
回代法里你做的是:
一层一层代
每次更新「当前余数 = 前两个余数的线性组合」
表格法只是把:
1 = a i ⋅ 119 + b i ⋅ 67 1 = a_i\cdot 119 + b_i\cdot 67 1 = a i ⋅ 119 + b i ⋅ 67
中的 b i b_i b i (67 的系数)提前算出来。
口诀: 先排好商 q q q ,从 0 , 1 0, 1 0 , 1 开始,计算公式为 x i = x i − 2 − q i − 1 × x i − 1 x_i = x_{i-2} - q_{i-1} \times x_{i-1} x i = x i − 2 − q i − 1 × x i − 1 。
先做辗转相除得到商 q q q :
119 ÷ 67 = 1 119 \div 67 = \mathbf{1} 119 ÷ 67 = 1 余 52 52 52
67 ÷ 52 = 1 67 \div 52 = \mathbf{1} 67 ÷ 52 = 1 余 15 15 15
52 ÷ 15 = 3 52 \div 15 = \mathbf{3} 52 ÷ 15 = 3 余 7 7 7
15 ÷ 7 = 2 15 \div 7 = \mathbf{2} 15 ÷ 7 = 2 余 1 1 1
真正被用到的不是余数,而是“商” 。表格法 只吃商,不管余数
q = [ 1 , ; 1 , ; 3 , ; 2 ] \boxed{q = [1,;1,;3,;2]} q = [ 1 ,; 1 ,; 3 ,; 2 ]
填表:
步骤 商 q 序列 x (逆元所在列) 计算过程 -2 0 (初始值) -1 1 (初始值) 1 1 -1 0 − ( 1 × 1 ) = − 1 0 - (1 \times 1) = -1 0 − ( 1 × 1 ) = − 1 2 1 2 1 − ( 1 × − 1 ) = 2 1 - (1 \times -1) = 2 1 − ( 1 × − 1 ) = 2 3 3 -7 − 1 − ( 3 × 2 ) = − 7 -1 - (3 \times 2) = -7 − 1 − ( 3 × 2 ) = − 7 4 2 16 2 − ( 2 × − 7 ) = 16 2 - (2 \times -7) = 16 2 − ( 2 × − 7 ) = 16
结果: 最后一行的 16 就是逆元。
4. 反复平方算法
反复平方算法的目标:
用 O ( log n ) O (\log n) O ( log n ) 次运算
代替 O ( n ) O (n) O ( n ) 次连乘
把指数写成二进制,把幂次拆成若干个平方
也就是:
a n = a 2 1 k ⋅ a 2 2 k ⋯ a^n = a^{2^k_1}\cdot a^{2^k_2}\cdots a n = a 2 1 k ⋅ a 2 2 k ⋯
而每个 a 2 k a^{2^k} a 2 k 都可以通过“平方前一个”得到。
1️⃣ 指数拆分
n = ∑ b i 2 i ( b i ∈ 0 , 1 ) n = \sum b_i 2^i \quad (b_i\in{0,1}) n = ∑ b i 2 i ( b i ∈ 0 , 1 )
2️⃣ 模运算与乘法兼容
( a b ) m o d m = ( ( a m o d m ) ( b m o d m ) ) m o d m (ab)\bmod m = ((a\bmod m)(b\bmod m))\bmod m ( ab ) mod m = (( a mod m ) ( b mod m )) mod m
👉 所以我们可以:
算法步骤
Step 1:指数二进制化
把 n 写成二进制。
Step 2:反复平方(预处理)
依次计算:
\begin{aligned} a^{2^0} &\equiv a \pmod m \ a^{2^1} &\equiv a^2 \pmod m \ a^{2^2} &\equiv $a^2$^2 \pmod m \ a^{2^3} &\equiv $a^2$^4 \pmod m \ &\vdots \end{aligned}
Step 3:按二进制位相乘
只乘 二进制中为 1 的位 对应的幂。
示例
计算 14 3 ( m o d 55 ) 14^3 \pmod{55} 1 4 3 ( mod 55 )
虽然指数小,但流程完全一样 。
Step 1:指数二进制化
3 = ( 11 ) 2 = 2 1 + 2 0 3 = (11)_2 = 2^1 + 2^0 3 = ( 11 ) 2 = 2 1 + 2 0
Step 2:反复平方(每一步都取模)
第 1 项:
14 2 0 ≡ 14 ( m o d 55 ) 14^{2^0} \equiv 14 \pmod{55} 1 4 2 0 ≡ 14 ( mod 55 )
第 2 项:
14 2 1 = 14 2 = 196 14^{2^1} = 14^2 = 196 1 4 2 1 = 1 4 2 = 196
196 m o d 55 = 196 − 3 × 55 = 31 196 \bmod 55 = 196 - 3\times55 = 31 196 mod 55 = 196 − 3 × 55 = 31
所以:
14 2 ≡ 31 ( m o d 55 ) 14^2 \equiv 31 \pmod{55} 1 4 2 ≡ 31 ( mod 55 )
Step 3:选中二进制为 1 的幂相乘
因为:
14 3 = 14 2 1 ⋅ 14 2 0 14^3 = 14^{2^1}\cdot14^{2^0} 1 4 3 = 1 4 2 1 ⋅ 1 4 2 0
所以:
14 3 ≡ 31 × 14 ( m o d 55 ) 14^3 \equiv 31 \times 14 \pmod{55} 1 4 3 ≡ 31 × 14 ( mod 55 )
计算:
31 × 14 = 434 31\times14 = 434 31 × 14 = 434
434 m o d 55 = 434 − 7 × 55 = 49 434 \bmod 55 = 434 - 7\times55 = 49 434 mod 55 = 434 − 7 × 55 = 49
🎯 最终答案
14 3 ≡ 49 ( m o d 55 ) \boxed{14^3 \equiv 49 \pmod{55}} 1 4 3 ≡ 49 ( mod 55 )
表格解法
计算
10 27 ( m o d 55 ) 10^{27} \pmod{55} 1 0 27 ( mod 55 )
Step 1:指数二进制化
27 = ( 11011 ) 2 = 16 + 8 + 2 + 1 27 = (11011)_2 = 16+8+2+1 27 = ( 11011 ) 2 = 16 + 8 + 2 + 1
Step 2:表格法
i e i e_i e i 10 2 i m o d 55 10^{2^i} \bmod 55 1 0 2 i mod 55 累积结果 0 1 10 10 10 1 × 10 = 10 1\times10=10 1 × 10 = 10 1 1 10 2 = 100 ≡ 45 10^2=100\equiv45 1 0 2 = 100 ≡ 45 10 × 45 ≡ 10 10\times45\equiv10 10 × 45 ≡ 10 2 0 45 2 = 2025 ≡ 45 45^2=2025\equiv45 4 5 2 = 2025 ≡ 45 不变 3 1 45 2 ≡ 45 45^2\equiv45 4 5 2 ≡ 45 10 × 45 ≡ 10 10\times45\equiv10 10 × 45 ≡ 10 4 1 45 2 ≡ 45 45^2\equiv45 4 5 2 ≡ 45 10 × 45 ≡ 10 10\times45\equiv10 10 × 45 ≡ 10
最终结果
10 27 ≡ 10 ( m o d 55 ) \boxed{10^{27} \equiv 10 \pmod{55}} 1 0 27 ≡ 10 ( mod 55 )
5. 时间复杂度
RSA 加密与解密的主要计算为模指数运算,采用快速幂算法,其时间复杂度为 O ( log n ) O(\log n) O ( log n ) 次模乘运算。
6. RSA 的计算可行性和不可攻破性
可行性
已知 p,q,就能:轻松算φ(n),用扩展 Euclid 算 d
模指数运算可以用快速幂
RSA 的密钥生成和加解密过程均可通过多项式时间算法完成,尤其是模指数运算可采用反复平方算法,因此在计算上是可行的。
不可攻破性
已知 n,不知道 p,q,无法得到μ(n),无法求 d
在现有计算能力和已知算法条件下,大整数因子分解问题是计算上不可行的,因此 RSA 具有不可攻破性。
0-1 背包问题
1️⃣ 子问题定义(非常关键)
令 dp[i][j] 表示:
从前 i 个物品中选择,背包容量为 j 时可获得的最大价值。
2️⃣ 边界条件
dp[0][j] = 0(没有物品)
dp[i][0] = 0(容量为 0)
例题
dp[i][j] =
在「只考虑前 i 个物品」的情况下,
背包容量为 j 时,能得到的最大价值
⚠️ 注意三点:
只考虑前 i 个物品 (后面的物品当不存在)
背包容量固定是 j
这个格子一定是最优解
也就是说:
“每一个单元格,都是当前条件下的最优解”
看图中第 0 行(i = 0):
不考虑任何物品
那不管背包多大:能装的最大价值 = 0
每填一个格子,你只纠结一件事:
👉 “当前这个物品,我要不要放?”
只有两种选择,没有第三种。
规则 1:如果放不下:只能不放
条件:当前物品重量 > 当前背包容量
那就:dp[i][j] = dp[i-1][j]
👉 直接抄上一行
例子:
2 > 1,放不下,所以:dp[1][1] = dp[0][1] = 0
规则 2:如果放得下:比较「放 vs 不放」
假设:
当前是第 i 个物品
重量 = w
价值 = v
当前背包容量 = j
两种选择
✅ 选择 1:不放
那就是:价值 = dp[i-1][j]
✅ 选择 2:放
放了之后:
占用 w 的容量
剩余容量 = j - w
剩余容量里还能放什么?
👉 只能从前 i-1 个物品里选
所以:价值 = v + dp[i-1][j-w],取最大值
dp[i][j] = max(
dp[i-1][j], // 不放
v + dp[i-1][j-w] // 放
)
也就是说, “放进去划算还是不放划算”
带你手算一个格子(非常重要)
我们看图中这一格,第 2 行、第 5 列
1️⃣ 放得下吗?3 ≤ 5 ✔ 放得下
2️⃣ 不放的价值:dp[1][5] = 3(只放紫球)
3️⃣ 放的价值
剩余容量:5 - 3 = 2
查上一行第 2 列:dp[1][2] = 3
加上水瓶的价值:5 + 3 = 8
4️⃣ 取最大值:max(3, 8) = 8,✔ 所以这一格是 8
为什么一定“查上一行”?因为这是 0/1 背包 :
每个物品只能用一次
如果你查的是当前行 ,那就相当于允许一个物品用多次
那就变成了完全背包 (另一个问题)
最后一格为什么是 9?右下角:
背包容量 6
可以放:
水瓶(3,5)
西瓜(4,6) ❌ 超重
紫球(2,3)
最优组合是:紫球 + 西瓜 = 2 + 4 = 6 价值 = 3 + 6 = 9
✔ 所以答案是 9
DP 表格回溯
目标:
不仅知道最大价值是多少,还要知道“装了什么”
1️⃣ 回溯的核心思想(一句话)
看这一格的值,是不是“来自上面”,还是“来自左上 + 当前物品”
换句话说:
如果 dp[i][j] == dp[i-1][j] 👉 第 i 个物品没选
如果 dp[i][j] != dp[i-1][j] 👉 第 i 个物品选了
⚠️ 注意:
前提是 weight[i] ≤ j,否则根本不可能选
2️⃣ 回溯从哪里开始?
从右下角开始
dp[n][W]
也就是:
3️⃣ 回溯规则(非常重要)
假设你现在在:
dp[i][j]
情况 A:dp[i][j] == dp[i-1][j]
说明: 不选第 i 个物品,也能达到当前最优值
👉 所以:第 i 个物品没选 i = i - 1 j 不变
情况 B:dp[i][j] != dp[i-1][j]
说明: 当前最优解必须选了第 i 个物品
👉 所以:第 i 个物品选了 j = j - weight[i] i = i - 1
4️⃣ 用你那张图完整走一遍回溯
最终答案是:dp[3][6] = 9
第一步:看 dp[3][6]
dp[3][6] = 9
dp[2][6] = 8
不相等 ⇒ 选了第 3 个物品(西瓜)
j = 6 - 4 = 2 i = 2
第二步:看 dp[2][2]
dp[2][2] = 3
dp[1][2] = 3
相等 ⇒ 没选第 2 个物品(水瓶)
i = 1 j = 2
第三步:看 dp[1][2]
dp[1][2] = 3
dp[0][2] = 0
不相等 ⇒ 选了第 1 个物品(紫球)
j = 2 - 2 = 0 i = 0
结束
最终选了:
✅ 物品 1(紫球)
✅ 物品 3(西瓜)
DP 方程的真正含义
你现在看到的方程是:
dp[i][j] = max(
dp[i-1][j],
dp[i-1][j - w[i]] + v[i]
)
我们来逐句拆解人话含义 。
1️⃣ dp[i][j] 在“回答什么问题”?
在前 i 个物品里选,
用容量 j 的背包,
最多能赚多少钱?
2️⃣ 为什么只有两种情况?
因为 0/1 背包的本质限制 :
每个物品
不存在:
3️⃣ 为什么选了就一定是 dp[i-1][j-w]?
因为:
你已经用了第 i 个物品一次
剩下只能从 前 i-1 个物品里选
剩余容量是 j - w[i]
👉 这不是技巧,是问题定义强制的
4️⃣ 为什么取 max 就一定对?
因为,最优解一定来自“所有合法方案中价值最大的那个” ,而:
所有合法方案 = 放 / 不放 两种情况
max 就是在这两类里选最优
三、怎么证明 0/1 背包具有「最优子结构」?
这是理论核心,我们来严谨但不装逼地讲清楚 。
1️⃣ 什么叫最优子结构?
一个问题的最优解
可以由其子问题的最优解组合得到
2️⃣ 对背包问题,子问题是什么?
主问题:前 i 个物品,容量 j 的最优解
子问题:
前 i-1 个物品,容量 j
前 i-1 个物品,容量 j - w[i]
3️⃣ 正式证明思路(关键)
我们用 反证法
假设:
S 是问题 (i, j) 的一个最优解,且 选了第 i 个物品
那么:
S = {第 i 个物品} ∪ S’
S’ 是在 前 i-1 个物品 、容量 j-w[i] 下的一个解
关键一步(最重要)
如果:
S’ 不是 (i-1, j-w[i]) 的最优解
那么:
存在另一个解 S”,价值更大
用 S” + 第 i 个物品
👉 会得到比 S 更优的解
这 与 S 是最优解矛盾
所以:
S’ 必须是 (i-1, j-w[i]) 的最优解
✔ 这就证明了:
最优解包含最优子解
4️⃣ 不选第 i 个物品的情况同理
如果最优解没选第 i 个物品:
那它就是 (i-1, j) 的最优解
否则同样可以构造更优解,产生矛盾
5️⃣ 结论
0/1 背包问题具有最优子结构
所以可以用动态规划解决
排序
快速排序
快速排序采用分而治之的思想,首先从待排序序列中选取一个元素作为基准(pivot),通过一次划分操作将序列分为两部分,使得基准左侧的元素均不大于基准,右侧的元素均不小于基准;然后分别对左右两个子序列递归地进行快速排序,直到子序列规模为 1 或 0 为止。
QuickSort(a, low, high):
if low >= high then return
pivot = Partition(a, low, high)
QuickSort(a, low, pivot - 1)
QuickSort(a, pivot + 1, high)
【数据结构合集 - 快速排序(算法过程, 效率分析, 稳定性分析)】 https://www.bilibili.com/video/BV1y4421Z7hK/
最好:折半查找
最坏:pivot 取最小/最大值
Partition (划分) 函数的伪代码是必考的。你只会在纸上画表是不够的,必须能默写出类似下面的逻辑:
输入 :数组 A [ p . . r ] A[p..r] A [ p .. r ]
逻辑 :
选主元 x = A [ r ] x = A[r] x = A [ r ] (通常选最后一个)
i = p − 1 i = p - 1 i = p − 1
遍历 j j j 从 p p p 到 r − 1 r-1 r − 1 :
如果 A [ j ] ≤ x A[j] \le x A [ j ] ≤ x ,则 i + + i++ i + + ,交换 A [ i ] A[i] A [ i ] 和 A [ j ] A[j] A [ j ]
最后交换 A [ i + 1 ] A[i+1] A [ i + 1 ] 和 A [ r ] A[r] A [ r ]
返回 i + 1 i+1 i + 1
建议: 不要死记硬背具体的 C++ 代码,但要背熟上面这个伪代码的逻辑步骤 ,因为试卷是按行给分的。
2. 细节缺失:复杂度的“递推式”
往年试卷不仅问你复杂度是多少,还要求写出 “复杂度的递推式及其闭式解” 4。
你现在的状态 :记住了 O ( n log n ) O(n \log n) O ( n log n ) 。
考试要求 :你可能需要写出:
最好情况递推式: T ( n ) = 2 T ( n / 2 ) + Θ ( n ) T(n) = 2T(n/2) + \Theta(n) T ( n ) = 2 T ( n /2 ) + Θ ( n )
最坏情况递推式: T ( n ) = T ( n − 1 ) + Θ ( n ) T(n) = T(n-1) + \Theta(n) T ( n ) = T ( n − 1 ) + Θ ( n ) (这是快排特有的,归并排序没有这个问题)
3. 场景补充:最坏情况的触发条件
除了记忆“最坏是 O ( n 2 ) O(n^2) O ( n 2 ) ”,你必须知道什么时候 发生最坏情况。
常考简答题: “在什么情况下快速排序会退化为冒泡排序/效率最低?”
答案: 当每次选择的 pivot 都是当前数组中的最小或最大元素(如已排序或逆序数组)时, 快速排序退化为 T(n)=T(n−1)+Θ(n),时间复杂度为 Θ(n²)。
归并排序
核心思想:分而治之(Divide and Conquer)
将数组一分为二
对左右两部分分别进行归并排序
合并两个有序数组
堆排序
【数据结构合集 - 堆与堆排序(算法过程, 效率分析, 稳定性分析)】 https://www.bilibili.com/video/BV1HYtseiEQ8/
主定理
时间复杂度
蛮力法
分治法
算法基础
一、算法基础(共 25 分)
一(1)给出 Knuth 的算法定义及其特征(5 分)
【标准定义(必背)】
Knuth 对算法的定义:
算法是一个有限的、确定的、可行的指令序列 ,该序列对给定的输入产生相应的输出,并在有限步骤内终止。
【算法的五个特征(每点都要写)】
输入(Input)
算法有 0 个或多个输入。
输出(Output)
算法至少有 1 个输出。
确定性(Definiteness)
每一条指令含义明确、无二义性。
有限性(Finiteness)
算法在有限步骤后终止。
可行性(Effectiveness)
每一步操作都可以在有限时间内完成。
一(2)给出渐近精确界 Θ 的定义及极限判定准则(5 分)
【Θ 的定义(2–3 行最标准)】
若存在正常数 c 1 , c 2 , n 0 c_1, c_2, n_0 c 1 , c 2 , n 0 ,使得对所有 n ≥ n 0 n \ge n_0 n ≥ n 0 ,都有
c 1 g ( n ) ≤ f ( n ) ≤ c 2 g ( n ) c_1 g(n) \le f(n) \le c_2 g(n) c 1 g ( n ) ≤ f ( n ) ≤ c 2 g ( n ) 则称 f ( n ) f(n) f ( n ) 的渐近精确界为 Θ ( g ( n ) ) \Theta(g(n)) Θ ( g ( n )) 。
【极限判定准则(必考)】
若
lim n → ∞ f ( n ) g ( n ) = c , 0 < c < ∞ \lim_{n \to \infty} \frac{f(n)}{g(n)} = c,\quad 0 < c < \infty n → ∞ lim g ( n ) f ( n ) = c , 0 < c < ∞ 则
f ( n ) = Θ ( g ( n ) ) f(n) = \Theta(g(n)) f ( n ) = Θ ( g ( n ))
📌 给分点提示
一定要出现 “常数 c₁、c₂、n₀”
极限要写成 正的有限常数
一(3)画出基本 NP 完全问题归约树,并写出问题中文名(5 分)
【经典 NP 完全归约链(你可以“画成文字”)】
SAT
↓
3-SAT
↓
CLIQUE(团问题)
↓
VERTEX COVER(顶点覆盖)
↓
INDEPENDENT SET(独立集)
【中文名称对应(一定要写)】
英文 中文 SAT 布尔可满足性问题 3-SAT 3-子句可满足性问题 CLIQUE 团问题 VERTEX COVER 顶点覆盖问题 INDEPENDENT SET 独立集问题
一(4)列出常见渐近函数的阶(含中文描述)(5 分)
【标准表(建议直接照写)】
函数 阶 中文描述 $1 O(1) 常数阶 log n O(log n) 对数阶 n O(n) 线性阶 n log n O(n log n) 线性对数阶 n² O(n²) 平方阶 n³ O(n³) 立方阶 2ⁿ O(2ⁿ) 指数阶 n! O(n!) 阶乘阶
一(5)手写多项式合并函数的大师定理(5 分)
【Master 定理标准版】
设递推式为:
T ( n ) = a T ( n b ) + f ( n ) , a ≥ 1 , b > 1 T(n) = aT\left(\frac{n}{b}\right) + f(n),\quad a \ge 1, b > 1 T ( n ) = a T ( b n ) + f ( n ) , a ≥ 1 , b > 1
三种情况(一定要分条写)
若 f ( n ) = O ( n log b a − ε ) f(n) = O(n^{\log_b a - \varepsilon}) f ( n ) = O ( n l o g b a − ε )
T ( n ) = Θ ( n log b a ) T(n) = \Theta(n^{\log_b a}) T ( n ) = Θ ( n l o g b a )
若 f ( n ) = Θ ( n log b a ) f(n) = \Theta(n^{\log_b a}) f ( n ) = Θ ( n l o g b a )
T ( n ) = Θ ( n log b a log n ) T(n) = \Theta(n^{\log_b a}\log n) T ( n ) = Θ ( n l o g b a log n )
若 f ( n ) = Ω ( n log b a + ε ) f(n) = \Omega(n^{\log_b a + \varepsilon}) f ( n ) = Ω ( n l o g b a + ε )
T ( n ) = Θ ( f ( n ) ) T(n) = \Theta(f(n)) T ( n ) = Θ ( f ( n ))
二、算法分析(共 25 分)
二(1)动态规划适用问题的两个特征(5 分)
【标准答案】
最优子结构
原问题的最优解可以由子问题的最优解构成。
重叠子问题
在递归求解过程中,不同子问题会重复出现。
二(2)说明贪心法与动态规划的不同(5 分)
【标准对比答案】
对比 贪心 动态规划 决策方式 局部最优 全局最优 是否回溯 否 否(全面考虑) 子问题 一般不重叠 有重叠 结果 不一定最优 一定最优
【一句话总结(加分)】
贪心算法通过局部最优选择期望得到全局最优解,而动态规划通过系统地求解并保存子问题最优解来保证全局最优。
二(3)Levenshtein 编辑距离定义及公式(5 分)
【定义】
编辑距离是将一个字符串转换为另一个字符串所需的最少编辑操作次数,允许的操作包括插入、删除和替换。
【状态转移方程】
令 dp[i][j] 表示将 X 前 i 个字符变为 Y 前 j 个字符的最小操作数:
d p [ i ] [ j ] = min { d p [ i − 1 ] [ j ] + 1 ( 删除 ) d p [ i ] [ j − 1 ] + 1 ( 插入 ) d p [ i − 1 ] [ j − 1 ] + c o s t ( 替换 ) dp[i][j] =
\min
\begin{cases}
dp[i-1][j] + 1 & (\text{删除}) \
dp[i][j-1] + 1 & (\text{插入}) \
dp[i-1][j-1] + cost & (\text{替换})
\end{cases} d p [ i ] [ j ] = min { d p [ i − 1 ] [ j ] + 1 ( 删除 ) d p [ i ] [ j − 1 ] + 1 ( 插入 ) d p [ i − 1 ] [ j − 1 ] + cos t ( 替换 )
二(4)编辑距离 DP 伪代码及复杂度(5 分)
【伪代码(简洁考试版)】
for i = 0 to m:
dp[i][0] = i
for j = 0 to n:
dp[0][j] = j
for i = 1 to m:
for j = 1 to n:
if X[i] == Y[j]:
dp[i][j] = dp[i-1][j-1]
else
dp[i][j] = 1 + min(dp[i-1][j],
dp[i][j-1],
dp[i-1][j-1])
【复杂度】
时间复杂度:O ( m n ) O(mn) O ( mn )
空间复杂度:O ( m n ) O(mn) O ( mn )
二(5)矩阵链乘法 DP 方程及伪代码(5 分)
【DP 方程】
m [ i ] [ j ] = { 0 , i = j min i ≤ k < j m [ i ] [ k ] + m [ k + 1 ] [ j ] + p i − 1 p k p j , i < j m[i][j] =
\begin{cases}
0, & i=j \
\min_{i\le k<j} {m[i][k] + m[k+1][j] + p_{i-1}p_kp_j}, & i<j
\end{cases} m [ i ] [ j ] = { 0 , i = j min i ≤ k < j m [ i ] [ k ] + m [ k + 1 ] [ j ] + p i − 1 p k p j , i < j
【伪代码(简版)】
for i = 1 to n:
m[i][i] = 0
for l = 2 to n:
for i = 1 to n-l+1:
j = i + l - 1
m[i][j] = ∞
for k = i to j-1:
m[i][j] = min(m[i][j],
m[i][k] + m[k+1][j] + p[i-1]*p[k]*p[j])
第二部分:核心概念精要
1. 大 O 记法定义及其极限判别准则
对于给定的复杂度函数 T ( n ) T(n) T ( n ) ,如果存在两个正数 C C C 和 n 0 n_0 n 0 ,以及函数 τ ( n ) \tau(n) τ ( n ) ,使得当 n ≥ n 0 n \ge n_0 n ≥ n 0 时有 T ( n ) ≤ C τ ( n ) T(n) \le C \tau(n) T ( n ) ≤ C τ ( n ) ,则称 T ( n ) = O ( τ ( n ) ) T(n) = O(\tau(n)) T ( n ) = O ( τ ( n )) 。
T ( n ) = O ( τ ( n ) ) ,当且仅当 0 ≤ lim n → ∞ T ( n ) τ ( n ) < + ∞ T(n) = O(\tau(n)) \text{,当且仅当 } 0 \le \lim_{n \to \infty} \frac{T(n)}{\tau(n)} < +\infty T ( n ) = O ( τ ( n )) ,当且仅当 0 ≤ n → ∞ lim τ ( n ) T ( n ) < + ∞
即 lim n → ∞ T ( n ) τ ( n ) ∈ [ 0 , + ∞ ) \lim_{n \to \infty} \frac{T(n)}{\tau(n)} \in [0, +\infty) lim n → ∞ τ ( n ) T ( n ) ∈ [ 0 , + ∞ ) 成立。
2. Karatsuba 乘法算法基本思想、伪代码、复杂度、闭式解,递推式
他采用与高斯复数乘法技巧相同的方法。
将 x L y R + x R y L x_L y_R + x_R y_L x L y R + x R y L 表述为 ( x L + x R ) ( y L + y R ) − x L y L − x R y R (x_L + x_R)(y_L + y_R) - x_L y_L - x_R y_R ( x L + x R ) ( y L + y R ) − x L y L − x R y R 。
其中 x L y L , x R y R x_L y_L, x_R y_R x L y L , x R y R 是重复计算,使乘法运算由 4 次降低为 3 次,尽管 add 了 3 次加法算法,但加法运算的复杂度是线性的,整体上降低了复杂度。
伪代码:
输入 :两个 n n n 位的正整数 X X X 和 Y Y Y
输出 :X X X 和 Y Y Y 的乘积
multiplyK(x, y)
n ← n \leftarrow n ← X X X 和 Y Y Y 的位数
if n = 1 return x * y
x L , x R = X x_L, x_R = X x L , x R = X 的高 n / 2 n/2 n /2 位和低 n / 2 n/2 n /2 位
y L , y R = Y y_L, y_R = Y y L , y R = Y 的高 n / 2 n/2 n /2 位和低 n / 2 n/2 n /2 位
P 1 = multiplyK ( x L , y L ) P_1 = \text{multiplyK}(x_L, y_L) P 1 = multiplyK ( x L , y L )
P 2 = multiplyK ( x R , y R ) P_2 = \text{multiplyK}(x_R, y_R) P 2 = multiplyK ( x R , y R )
P 3 = multiplyK ( x L + x R , y L + y R ) P_3 = \text{multiplyK}(x_L + x_R, y_L + y_R) P 3 = multiplyK ( x L + x R , y L + y R )
return P 1 × 10 n + ( P 3 − P 1 − P 2 ) × 10 n / 2 + P 2 P_1 \times 10^n + (P_3 - P_1 - P_2) \times 10^{n/2} + P_2 P 1 × 1 0 n + ( P 3 − P 1 − P 2 ) × 1 0 n /2 + P 2
复杂度分析:
T ( n ) = { O ( 1 ) n = 1 3 T ( n / 2 ) + O ( n ) n > 1 T(n) =
\begin{cases}
O(1) & n=1 \\
3T(n/2) + O(n) & n>1
\end{cases} T ( n ) = { O ( 1 ) 3 T ( n /2 ) + O ( n ) n = 1 n > 1
闭式解: T ( n ) = O ( n log 2 3 ) ≈ O ( n 1.58 ) T(n) = O(n^{\log_2 3}) \approx O(n^{1.58}) T ( n ) = O ( n l o g 2 3 ) ≈ O ( n 1.58 )
3. 快速排序基本思想、复杂度
基本思想:
从待排序的数列中选定一个作为枢纽元素的数字,然后采用一种分区算法,将枢纽元素排到正确位置上,使其前元素比它小,其后元素比它大。其他元素划分为两个区,然后再用分治策略利用上述过程求解前后两个分区,其直到分区中仅剩不超过 1 个元素。
复杂度:
最坏复杂度:O ( n 2 ) O(n^2) O ( n 2 )
平均/最好:O ( n log n ) O(n \log n) O ( n log n )
4. 多项式合并函数下的大师定理
当 T ( n ) = a T ( n / b ) + f ( n ) T(n) = aT(n/b) + f(n) T ( n ) = a T ( n / b ) + f ( n ) 中 f ( n ) = O ( n d ) f(n) = O(n^d) f ( n ) = O ( n d ) 时,大师定理可简化:
对常数 a ≥ 1 , b > 1 , d ≥ 0 a \ge 1, b > 1, d \ge 0 a ≥ 1 , b > 1 , d ≥ 0 ,有 T ( n ) = a T ( n / b ) + O ( n d ) T(n) = aT(n/b) + O(n^d) T ( n ) = a T ( n / b ) + O ( n d ) ,则 T ( n ) T(n) T ( n ) 有如下闭式解:
T ( n ) = { O ( n d ) d > log b a O ( n d log n ) d = log b a O ( n log b a ) d < log b a T(n) =
\begin{cases}
O(n^d) & d > \log_b a \\
O(n^d \log n) & d = \log_b a \\
O(n^{\log_b a}) & d < \log_b a
\end{cases} T ( n ) = ⎩ ⎨ ⎧ O ( n d ) O ( n d log n ) O ( n l o g b a ) d > log b a d = log b a d < log b a
5. Knuth 的算法定义及特性
算法是一组有穷规则的集合,这些规则给出了求解特定类型问题的运算序列。它具有 5 个重要特性,即:
有穷性 (Finiteness)
确定性 (Definiteness)
有效性 (Effectiveness)
输入 (Input)
输出 (Output)
6. Euclid GCD 手算表格
No a b q=⌊ a / b ⌋ \lfloor a/b \rfloor ⌊ a / b ⌋ r 1 98 56 1 42 2 56 42 1 14 3 42 14 3 0 4 14 0
7. 适应于动态规划方法求解的问题的 2 个特征
最优子结构性质 (Optimal Substructure)
子问题重叠性 (Overlapping Subproblems)
8. Levenshtein 编辑距离的定义及公式表达
定义:
两个字符串间的编辑距离是指将一个字符串变换为另一个字符串的基本编辑操作数。这里的编辑操作指的是对单个字符的插入 、删除 和替换 3 种基本编辑操作。
公式:
l e v a , b ( i , j ) = { max ( i , j ) min ( i , j ) = 0 min { l e v a , b ( i − 1 , j ) + 1 l e v a , b ( i , j − 1 ) + 1 l e v a , b ( i − 1 , j − 1 ) + δ ( a i , b j ) else lev_{a,b}(i, j) =
\begin{cases}
\max(i, j) & \min(i, j) = 0 \\
\min \begin{cases}
lev_{a,b}(i-1, j) + 1 \\
lev_{a,b}(i, j-1) + 1 \\
lev_{a,b}(i-1, j-1) + \delta(a_i, b_j)
\end{cases} & \text{else}
\end{cases} l e v a , b ( i , j ) = ⎩ ⎨ ⎧ max ( i , j ) min ⎩ ⎨ ⎧ l e v a , b ( i − 1 , j ) + 1 l e v a , b ( i , j − 1 ) + 1 l e v a , b ( i − 1 , j − 1 ) + δ ( a i , b j ) min ( i , j ) = 0 else
(注:公式中若 a i ≠ b j a_i \neq b_j a i = b j 则 δ = 1 \delta=1 δ = 1 ,否则 δ = 0 \delta=0 δ = 0 )
9. 给出 Levenshtein 编辑距离的 DP 算法伪代码、复杂度
复杂度: O ( m n ) O(mn) O ( mn )
伪代码:
Function LevenshteinDistance(S1, S2)
m = length(S1), n = length(S2)
Create DP[m+1][n+1]
// Initialization
From For i = 0 to m: DP[i][0] = i
For j = 0 to n: DP[0][j] = j
// Main Loop
For i = 1 to m:
For j = 1 to n:
If S1[i-1] == S2[j-1]:
cost = 0
else:
cost = 1
DP[i][j] = min(
DP[i-1][j] + 1, // Deletion
DP[i][j-1] + 1, // Insertion
DP[i-1][j-1] + cost // Substitution
)
return DP[m][n]
10. 0-1 背包问题的形式化定义
给定 n n n 个重量分别为 w 1 , w 2 , … , w n w_1, w_2, \dots, w_n w 1 , w 2 , … , w n ,价值分别为 v 1 , v 2 , … , v n v_1, v_2, \dots, v_n v 1 , v 2 , … , v n ,且不可分割的物品及一个容量 W W W 的背包,问装入哪些物品可以获得最大的价值?
11. 0-1 背包问题 DP 算法子问题的定义及 DP 方程
子问题定义:
将前 i i i (0 ≤ i ≤ n 0 \le i \le n 0 ≤ i ≤ n ) 个物品装入容量为 j j j (0 ≤ j ≤ W 0 \le j \le W 0 ≤ j ≤ W ) 的背包,问最大可以获得多少价值?
DP 方程:
V ( i , j ) = { 0 i = 0 or j = 0 V ( i − 1 , j ) w i > j , 1 ≤ i ≤ n max { V ( i − 1 , j ) , V ( i − 1 , j − w i ) + v i } w i ≤ j , 1 ≤ i ≤ n V(i, j) =
\begin{cases}
0 & i=0 \text{ or } j=0 \\
V(i-1, j) & w_i > j, 1 \le i \le n \\
\max \{ V(i-1, j), V(i-1, j-w_i) + v_i \} & w_i \le j, 1 \le i \le n
\end{cases} V ( i , j ) = ⎩ ⎨ ⎧ 0 V ( i − 1 , j ) max { V ( i − 1 , j ) , V ( i − 1 , j − w i ) + v i } i = 0 or j = 0 w i > j , 1 ≤ i ≤ n w i ≤ j , 1 ≤ i ≤ n
复杂度: T ( n , W ) = O ( n W ) T(n, W) = O(nW) T ( n , W ) = O ( nW )
0-1 背包问题的动态规划算法的伪代码
1. int Dp0-1Knapsack(int n, int W, int *w, int *v)
2. {
3. for (int i = 1; i <= n; i++)
4. for (int j = 1; j <= W; j++)
5. if (j < w[i-1])
6. V[i][j] = V[i-1][j];
7. else if (V[i-1][j] >= V[i-1][j-w[i-1]] + v[i-1]) // 这里笔记有一处判断逻辑
8. V[i][j] = V[i-1][j];
9. else
10. V[i][j] = V[i-1][j-w[i-1]] + v[i-1];
11.
11. // Backtracking to find items
12. int j = W;
13. for (int i = n; i > 0; i--)
14. if (V[i][j] == V[i-1][j])
15. x[i] = 0;
16. else {
17. x[i] = 1; j -= w[i-1];
18. }
19. return V[n][W];
20. }
12. 设有价值为 5, 8, 4, 3,重量为 4, 6, 5, 4 和容量为 10 的背包,设计并给出 DP 计算的手写表
物品/容量 (W = 10 W=10 W = 10 ):
i = 1 : w = 4 , v = 5 i=1: w=4, v=5 i = 1 : w = 4 , v = 5
i = 2 : w = 6 , v = 8 i=2: w=6, v=8 i = 2 : w = 6 , v = 8
i = 3 : w = 5 , v = 4 i=3: w=5, v=4 i = 3 : w = 5 , v = 4
i = 4 : w = 4 , v = 3 i=4: w=4, v=3 i = 4 : w = 4 , v = 3
DP 表格 (部分还原):
i i i \ j j j 0 1 2 3 4 5 6 7 8 9 10 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 5 5 5 5 5 5 5 2 0 0 0 0 5 5 8 8 8 8 13 3 0 0 0 0 5 5 8 8 8 9 13 4 0 0 0 0 5 5 8 8 8 9 13
结论: 最优解是物品 2 和 4 (注:根据价值 13 推算,应为物品 1 和 2,但笔记写的是 2 和 4),最优解值 13。
13. 给出 RSA 算法的主流程代码,说明其不可攻破性和可行性,复杂度
不可攻破性: 基于大整数分解的困难性。
可行性: RSA 的安全性基于大整数分解的困难性。在已知公钥 ( e , n ) (e, n) ( e , n ) 的情况下,想要推导出私钥 d d d ,必须对 n n n 进行因式分解,这在计算上是不可行的(Computationally Infeasible)。而 RSA 的加密和解密操作利用模幂运算,在计算上是高效可行的。 ”
复杂度: O ( ln n ) O(\ln n) O ( ln n ) (注:笔记手写标注)
流程:
取两个素数 p p p 和 q q q 。Miller-Rabin 素性测试算法。
计算 n = p q n = pq n = pq , ϕ ( n ) = ( p − 1 ) ( q − 1 ) \phi(n) = (p-1)(q-1) ϕ ( n ) = ( p − 1 ) ( q − 1 ) 。大整数乘法。
取一个与偶数 ϕ ( n ) \phi(n) ϕ ( n ) 互素的数 e e e 。Euclid GCD 算法。
计算 e ( m o d ϕ ( n ) ) e \pmod{\phi(n)} e ( mod ϕ ( n )) 的乘法逆元 d d d 。扩展欧几里得算法。
置 P = ( e , n ) P=(e, n) P = ( e , n ) 为 RSA 公钥, S = ( d , n ) S=(d, n) S = ( d , n ) 为 RSA 私钥。
对 0 ≤ P ≤ n − 1 0 \le P \le n-1 0 ≤ P ≤ n − 1 加密,C = P e m o d n C = P^e \mod n C = P e mod n 。模幂运算的反复平方算法。
对 C C C 解密,P = C d m o d n P = C^d \mod n P = C d mod n 。模幂运算的反复平方算法。
14. 设 p=3, q=11, 进行 RSA 运算
计算 n , ϕ ( n ) n, \phi(n) n , ϕ ( n ) ,确定与 ϕ ( n ) \phi(n) ϕ ( n ) 互素的 e e e ,计算乘法逆元 d d d ,然后用公钥 ( e , n ) (e, n) ( e , n ) 和私钥 ( d , n ) (d, n) ( d , n ) 对 30 进行加密和解密。
n = p × q = 33 n = p \times q = 33 n = p × q = 33
ϕ ( n ) = ( p − 1 ) ( q − 1 ) = 20 \phi(n) = (p-1)(q-1) = 20 ϕ ( n ) = ( p − 1 ) ( q − 1 ) = 20
令 e = 7 e = 7 e = 7
gcd ( e , ϕ ( n ) ) = 1 \gcd(e, \phi(n)) = 1 g cd( e , ϕ ( n )) = 1
扩展欧几里得求逆元 d d d :
No a = e a=e a = e b = ϕ ( n ) b=\phi(n) b = ϕ ( n ) q = ⌊ a / b ⌋ q=\lfloor a/b \rfloor q = ⌊ a / b ⌋ r r r 1 7 20 0 7 2 20 7 2 6 3 7 6 1 1 4 6 1 6 0 5 1 0
逆元推导过程:
∴ a 5 x 1 + b 5 x 0 = 1 b 4 x 1 + ( a 4 − q 4 b 4 ) x 0 = 1 ( a 3 − q 4 b 3 ) x 1 = 1 ∴ 1 b 2 − 6 ( a 2 − q 2 b 2 ) = 1 − 6 a 2 + 3 b 2 = 1 − 6 b 1 + 3 ( a 1 − q 1 b 1 ) = 1 3 a 1 − 6 b 1 = 1 \begin{aligned}
\therefore & a_5 x_1 + b_5 x_0 = 1 \\
& b_4 x_1 + (a_4 - q_4 b_4) x_0 = 1 \\
& (a_3 - q_4 b_3) x_1 = 1 \\
\therefore & 1b_2 - 6(a_2 - q_2 b_2) = 1 \\
& -6a_2 + 3b_2 = 1 \\
& -6b_1 + 3(a_1 - q_1 b_1) = 1 \\
& 3a_1 - 6b_1 = 1
\end{aligned} ∴ ∴ a 5 x 1 + b 5 x 0 = 1 b 4 x 1 + ( a 4 − q 4 b 4 ) x 0 = 1 ( a 3 − q 4 b 3 ) x 1 = 1 1 b 2 − 6 ( a 2 − q 2 b 2 ) = 1 − 6 a 2 + 3 b 2 = 1 − 6 b 1 + 3 ( a 1 − q 1 b 1 ) = 1 3 a 1 − 6 b 1 = 1
即乘法逆元 d = 3 d = 3 d = 3 。
加密与解密运算:
公钥 ( 7 , 33 ) (7, 33) ( 7 , 33 ) ,私钥 ( 3 , 33 ) (3, 33) ( 3 , 33 )
加密: C = m e m o d n = 30 7 m o d 33 C = m^e \mod n = 30^7 \mod 33 C = m e mod n = 3 0 7 mod 33
i i i 2 1 0 b b b 1 1 1 d 2 d^2 d 2 1 900 → 9 900 \to 9 900 → 9 36 → 3 36 \to 3 36 → 3 d 2 m o d n d^2 \mod n d 2 mod n 1 9 3 d ⋅ a d \cdot a d ⋅ a 30 270 90 d ⋅ a m o d n d \cdot a \mod n d ⋅ a mod n 30 6 24
即 C = 24 C=24 C = 24
解密: m = C d m o d n = 24 3 m o d 33 m = C^d \mod n = 24^3 \mod 33 m = C d mod n = 2 4 3 mod 33 (a = 24 , b = 3 , d = 1 , n = 33 a=24, b=3, d=1, n=33 a = 24 , b = 3 , d = 1 , n = 33 )
i i i 1 0 b b b 1 1 d 2 d^2 d 2 1 576 → 15 576 \to 15 576 → 15 d ⋅ a d \cdot a d ⋅ a 24 360 d ⋅ a m o d n d \cdot a \mod n d ⋅ a mod n 24 30
即 m = 30 m=30 m = 30
15. 归并算法基本思想、伪代码
Two Way Merge 基本思想:
顺次比较左右两个子数组中的元素
将比较小者放入临时数组
最后将左或右两个子数组中剩余元素复制到临时数组
将临时数组中的元素复制回原数组
伪代码 (TwoWayMerge):
输入 :数组 a a a , p ∼ q p \sim q p ∼ q 和 q + 1 ∼ r q+1 \sim r q + 1 ∼ r 的元素已排好序的数组
输出 :p ∼ r p \sim r p ∼ r 的元素已排好的数组 a a a
1. TwoWayMerge(a, p, q, r, b)
2. i = p, j = q + 1, k = p
3. for k = p to r
4. if (i <= q and (j > r or a[i] <= a[j]))
5. b[k] = a[i]; i = i + 1
6. else
7. b[k] = a[j]; j = j + 1
8. for k = p to r
9. a[k] = b[k]
复杂度 O ( n ) O(n) O ( n ) 。
归并排序 (MergeSort) 基本思想:
分解: 将排序数组分解为左右两个大小不重叠的子数组。
解决: 当子问题规模大于 1 时递归;当为 1 时则为平凡有序数组,直接返回。
合并: 将两个有序数组合并为一个有序数组。
伪代码 (MergeSort):
输入 :n n n 个元素的数组 a a a
输出 :排序的 a a a
1. MergeSort(a, low, high)
2. if low = high then return
3. mid = (low + high) / 2
4. MergeSort(a, low, mid)
5. MergeSort(a, mid+1, high)
6. TwoWayMerge(a, low, mid, high)
16. NP 难题
令 π \pi π 是一个判定性问题,如果对 NP 问题中的每个问题 π ′ ∈ N P \pi' \in NP π ′ ∈ N P 都有 π ′ ≤ P π \pi' \le_P \pi π ′ ≤ P π ,则称判定性问题 π \pi π 是一个 NP 难题 (NP-Hard)。
17. NP 完全问题
令 π \pi π 是一个判定性问题,如果 π ∈ N P \pi \in NP π ∈ N P ,且对 NP 中的每个问题 π ′ ∈ N P \pi' \in NP π ′ ∈ N P 都有 π ′ ≤ P π \pi' \le_P \pi π ′ ≤ P π ,则称判定性问题 π \pi π 是一个 NP 完全问题 (NPC)。
18. NP 完全问题归约树
graph TD
A [ 布尔电路可满足性问题 ] --> B [ 布尔可满足性问题 SAT ]
B --> C [ 3 - CNF 可满足性问题 ]
C --> D [ 团问题 ]
D --> E [ 顶点覆盖问题 ]
E --> F [ 子集和问题 ]
C --> G [ 哈密尔顿回路问题 ]
G --> H [ 旅行商问题 TSP ]
%% 样式调整(可选,为了更美观)
style A fill :# f9f , stroke :# 333 , stroke-width : 2px
style C fill :# bbf , stroke :# 333 , stroke-width : 2px
第三部分:专题详解
算法设计补考速成 考点
RSA 算法
基本原理:时间复杂度、可行性、不可攻破性
密钥生成
扩展欧几里得算法(含表格)
反复平方算法(手画表格)
算法基础
渐近时间复杂度:Θ \Theta Θ 精确界、O O O 上界、Ω \Omega Ω 下界的定义及极限判定;常见渐近函数的阶和中文描述
大师定理:手写多项式合并函数的大师定理公式
NP 和 NP 完全定理:P、NP、NP-Hard、NPC 的定义关系,归约的意义,常见 NPC 中文名称;经典归约链、归约树大致结构
算法定义 / Knuth 对算法的定义
分治排序算法
Karatsuba 大整数乘法
推导公式、复杂度分析。
核心思想
将原本 4 次的乘法降低到 3 次:
X = A ⋅ 10 n / 2 + B , Y = C ⋅ 10 n / 2 + D X = A \cdot 10^{n/2} + B, \qquad Y = C \cdot 10^{n/2} + D X = A ⋅ 1 0 n /2 + B , Y = C ⋅ 1 0 n /2 + D
X ⋅ Y = A C ⋅ 10 n + ( A D + B C ) ⋅ 10 n / 2 + B D X \cdot Y = AC \cdot 10^n + (AD + BC) \cdot 10^{n/2} + BD X ⋅ Y = A C ⋅ 1 0 n + ( A D + B C ) ⋅ 1 0 n /2 + B D
其中共有四次乘法:A C , A D , B C , B D AC, AD, BC, BD A C , A D , B C , B D 。
将
A D + B C AD + BC A D + B C
替换为
( A + B ) ( C + D ) − A C − B D (A+B)(C+D) - AC - BD ( A + B ) ( C + D ) − A C − B D
于是只需要三次乘法。
递推式与时间复杂度
T ( n ) = 3 T ( n / 2 ) + O ( n ) T(n) = 3T(n/2) + O(n) T ( n ) = 3 T ( n /2 ) + O ( n )
根据大师定理:a = 3 , b = 2 , d = 1 a=3, b=2, d=1 a = 3 , b = 2 , d = 1 ,且 3 > 2 1 3 > 2^1 3 > 2 1 ,故
T ( n ) = O ( n log 2 3 ) ≈ O ( n 1.585 ) T(n)=O\left(n^{\log_2 3}\right) \approx O(n^{1.585}) T ( n ) = O ( n l o g 2 3 ) ≈ O ( n 1.585 )
对比朴素乘法 O ( n 2 ) O(n^2) O ( n 2 ) ,属于典型分治优化。
归并 / 快速排序
要求:会写伪代码、会分析时间复杂度。
归并排序的 Merge 合并逻辑及伪代码
手算执行过程、闭解式
时间复杂度递推式
归并排序 / 快速排序常写为:
T ( n ) = 2 T ( n / 2 ) + Θ ( n ) T(n) = 2T(n/2) + \Theta(n) T ( n ) = 2 T ( n /2 ) + Θ ( n )
二叉最大堆
节点上移操作(shift up)的基本思想
伪代码
复杂度:O ( log n ) O(\log n) O ( log n )
会画树形演变过程
回溯法 VS 分支限界法
动态规划与贪心算法
0/1 背包与动态规划
概念与证明:0/1 背包形式化定义、DP 子问题及动态规划方程,证明最优子结构的正确性
手算表:给定一组物品重量、价值和背包容量,画出 DP 计算矩阵表,并通过表格回推最优解(即哪些物品被装入)
算法对比:说明贪心算法与动态规划的区别;贪心具有贪心选择策略,DP 具有重叠子问题。说明为什么 0/1 背包不能用贪心算法求解(需举例)
0/1 背包的 DP 方程、伪代码、复杂度、手算表
Levenshtein 编辑距离与矩阵链相乘
Levenshtein 编辑距离:定义、公式表达、DP 算法伪代码、时间复杂度
矩阵链相乘:动态规划方程及伪代码
动态规划与贪心法的不同
RSA 加密算法
加密过程
选择一对不相等且足够大的质数 p , q p, q p , q
计算
n = p × q n = p \times q n = p × q
计算 n n n 的欧拉函数
φ ( n ) = ( p − 1 ) ( q − 1 ) \varphi(n) = (p-1)(q-1) φ ( n ) = ( p − 1 ) ( q − 1 )
选一个与 φ ( n ) \varphi(n) φ ( n ) 互质的整数 e e e ,满足
1 < e < φ ( n ) 1 < e < \varphi(n) 1 < e < φ ( n )
计算出 e e e 对于 φ ( n ) \varphi(n) φ ( n ) 的模反元素 d d d ,使得
d e m o d φ ( n ) = 1 de \bmod \varphi(n) = 1 d e mod φ ( n ) = 1
公钥
K U = ( e , n ) KU = (e, n) K U = ( e , n )
私钥
K R = ( d , n ) KR = (d, n) K R = ( d , n )
模反元素 d d d
如果 e e e 和 φ ( n ) \varphi(n) φ ( n ) 互质,那么一定存在整数 d d d ,使得 e d − 1 ed - 1 e d − 1 能被 φ ( n ) \varphi(n) φ ( n ) 整除,也即:
e d ≡ 1 ( m o d φ ( n ) ) ed \equiv 1 \pmod{\varphi(n)} e d ≡ 1 ( mod φ ( n ))
等价表达:
e d − 1 = k φ ( n ) ed - 1 = k\varphi(n) e d − 1 = k φ ( n )
e d m o d φ ( n ) = 1 ed \bmod \varphi(n) = 1 e d mod φ ( n ) = 1
e d ≡ 1 ( m o d φ ( n ) ) ed \equiv 1 \pmod{\varphi(n)} e d ≡ 1 ( mod φ ( n ))
例如:3 3 3 和 11 11 11 ,有
3 × 4 − 1 = 11 3 \times 4 - 1 = 11 3 × 4 − 1 = 11
故 d = 4 d=4 d = 4 。更一般地,4 + k ⋅ 11 4 + k \cdot 11 4 + k ⋅ 11 都可以视为一个解,因此 d + k φ ( n ) d + k\varphi(n) d + k φ ( n ) 都是 e e e 的模反元素。通常取最小正整数解,便于计算。
明文 M M M 、密文 C C C 的加解密公式:
M e m o d n = C M^e \bmod n = C M e mod n = C
C d m o d n = M C^d \bmod n = M C d mod n = M
RSA 加密算法(例题)
密钥生成与加解密
设:
p = 5 , q = 11 , e = 3 p=5, \quad q=11, \quad e=3 p = 5 , q = 11 , e = 3
则
n = p q = 5 × 11 = 55 n=pq=5\times 11=55 n = pq = 5 × 11 = 55
φ ( n ) = ( p − 1 ) ( q − 1 ) = ( 5 − 1 ) ( 11 − 1 ) = 4 × 10 = 40 \varphi(n)=(p-1)(q-1)=(5-1)(11-1)=4\times 10=40 φ ( n ) = ( p − 1 ) ( q − 1 ) = ( 5 − 1 ) ( 11 − 1 ) = 4 × 10 = 40
求 d d d :
3 d ≡ 1 ( m o d 40 ) 3d \equiv 1 \pmod{40} 3 d ≡ 1 ( mod 40 )
可得
d = 27 d=27 d = 27
于是:
公钥:( n , e ) = ( 55 , 3 ) (n,e)=(55,3) ( n , e ) = ( 55 , 3 )
私钥:( n , d ) = ( 55 , 27 ) (n,d)=(55,27) ( n , d ) = ( 55 , 27 )
给定明文
m = 14 m=14 m = 14
加密公式:
C = m e m o d n = 14 3 m o d 55 C = m^e \bmod n = 14^3 \bmod 55 C = m e mod n = 1 4 3 mod 55
原稿后续用快速幂 / 取模化简演算:
14 2 = 196 ≡ 31 ( m o d 55 ) 14^2 = 196 \equiv 31 \pmod{55} 1 4 2 = 196 ≡ 31 ( mod 55 )
14 3 ≡ 31 × 14 = 434 ≡ 49 ( m o d 55 ) 14^3 \equiv 31 \times 14 = 434 \equiv 49 \pmod{55} 1 4 3 ≡ 31 × 14 = 434 ≡ 49 ( mod 55 )
故
C = 49 C = 49 C = 49
再解密:
49 27 m o d 55 = 14 49^{27} \bmod 55 = 14 4 9 27 mod 55 = 14
手稿中给出了连乘取模表,最终得到
m = 14 m=14 m = 14
扩展欧几里得算法求逆元
使用推广 Euclid 算法求 67 m o d 119 67 \bmod 119 67 mod 119 的逆元。
即求整数 x x x ,使得
67 x ≡ 1 ( m o d 119 ) 67x \equiv 1 \pmod{119} 67 x ≡ 1 ( mod 119 )
等价于求整数解:
67 x + 119 k = 1 67x + 119k = 1 67 x + 119 k = 1
只要
gcd ( 67 , 119 ) = 1 \gcd(67,119)=1 g cd( 67 , 119 ) = 1
就有整数解。扩展 Euclid 不仅判断两个数是否互质,还能构造出这样的 ( x , y ) (x,y) ( x , y ) 。
辗转相除法求最大公约数
119 = 1 × 67 + 52 119 = 1\times 67 + 52 119 = 1 × 67 + 52
67 = 1 × 52 + 15 67 = 1\times 52 + 15 67 = 1 × 52 + 15
52 = 3 × 15 + 7 52 = 3\times 15 + 7 52 = 3 × 15 + 7
15 = 2 × 7 + 1 15 = 2\times 7 + 1 15 = 2 × 7 + 1
7 = 1 × 1 + 0 7 = 1\times 1 + 0 7 = 1 × 1 + 0
最后一个非零余数是 1 1 1 ,因此
gcd ( 67 , 119 ) = 1 \gcd(67,119)=1 g cd( 67 , 119 ) = 1
所以逆元存在。
注意:不是所有的模数都有逆元。只有当
gcd ( a , m ) = 1 \gcd(a,m)=1 g cd( a , m ) = 1
时,a a a 才在模 m m m 下有逆元。
回代法求逆元
目标是把 1 1 1 表示成 67 67 67 和 119 119 119 的线性组合。
由
1 = 15 − 2 × 7 1 = 15 - 2\times 7 1 = 15 − 2 × 7
又因为
7 = 52 − 3 × 15 7 = 52 - 3\times 15 7 = 52 − 3 × 15
代入得
1 = 15 − 2 ( 52 − 3 × 15 ) 1 = 15 - 2(52 - 3\times 15) 1 = 15 − 2 ( 52 − 3 × 15 )
= 7 × 15 − 2 × 52 = 7\times 15 - 2\times 52 = 7 × 15 − 2 × 52
再由
15 = 67 − 1 × 52 15 = 67 - 1\times 52 15 = 67 − 1 × 52
得
1 = 7 ( 67 − 52 ) − 2 × 52 1 = 7(67 - 52) - 2\times 52 1 = 7 ( 67 − 52 ) − 2 × 52
= 7 × 67 − 9 × 52 = 7\times 67 - 9\times 52 = 7 × 67 − 9 × 52
再由
52 = 119 − 67 52 = 119 - 67 52 = 119 − 67
得
1 = 7 × 67 − 9 ( 119 − 67 ) 1 = 7\times 67 - 9(119 - 67) 1 = 7 × 67 − 9 ( 119 − 67 )
= 16 × 67 − 9 × 119 = 16\times 67 - 9\times 119 = 16 × 67 − 9 × 119
于是
1 = 16 × 67 + ( − 9 ) × 119 1 = 16\times 67 + (-9)\times 119 1 = 16 × 67 + ( − 9 ) × 119
所以
16 × 67 ≡ 1 ( m o d 119 ) 16\times 67 \equiv 1 \pmod{119} 16 × 67 ≡ 1 ( mod 119 )
故 67 67 67 在模 119 119 119 下的逆元为:
16 16 16
反复平方算法
目的:用 O ( log n ) O(\log n) O ( log n ) 次运算代替 O ( n ) O(n) O ( n ) 次连乘。
计算:
14 3 m o d 55 14^3 \bmod 55 1 4 3 mod 55
指数二进制化
3 = ( 11 ) 2 = 2 1 + 2 0 3=(11)_2=2^1+2^0 3 = ( 11 ) 2 = 2 1 + 2 0
反复平方(每一步取模)
14 2 0 ≡ 14 ( m o d 55 ) 14^{2^0} \equiv 14 \pmod{55} 1 4 2 0 ≡ 14 ( mod 55 )
14 2 = 196 ≡ 31 ( m o d 55 ) 14^2 = 196 \equiv 31 \pmod{55} 1 4 2 = 196 ≡ 31 ( mod 55 )
即:
14 2 1 ≡ 31 ( m o d 55 ) 14^{2^1} \equiv 31 \pmod{55} 1 4 2 1 ≡ 31 ( mod 55 )
按二进制位为 1 的幂相乘
14 3 = 14 2 1 ⋅ 14 2 0 14^3 = 14^{2^1} \cdot 14^{2^0} 1 4 3 = 1 4 2 1 ⋅ 1 4 2 0
14 3 ≡ 31 × 14 = 434 ≡ 49 ( m o d 55 ) 14^3 \equiv 31\times 14 = 434 \equiv 49 \pmod{55} 1 4 3 ≡ 31 × 14 = 434 ≡ 49 ( mod 55 )
原理模板
设
n = ∑ b i 2 i , b i ∈ { 0 , 1 } n = \sum b_i 2^i, \qquad b_i \in \{0,1\} n = ∑ b i 2 i , b i ∈ { 0 , 1 }
则
a n = ∏ i a b i 2 i a^n = \prod_i a^{b_i 2^i} a n = i ∏ a b i 2 i
并利用模运算性质:
( a b m o d m ) = ( ( a m o d m ) ( b m o d m ) ) m o d m (ab \bmod m)=((a\bmod m)(b\bmod m))\bmod m ( ab mod m ) = (( a mod m ) ( b mod m )) mod m
因此每一步都可以先取模,防止数字过大。
反复平方算法(表格解法)
计算:
10 27 m o d 55 10^{27} \bmod 55 1 0 27 mod 55
因为
27 = ( 11011 ) 2 = 16 + 8 + 2 + 1 27=(11011)_2=16+8+2+1 27 = ( 11011 ) 2 = 16 + 8 + 2 + 1
可整理为下表:
i i i e i e_i e i 10 2 i m o d 55 10^{2^i} \bmod 55 1 0 2 i mod 55 累积结果 0 1 10 1 × 10 = 10 1\times 10=10 1 × 10 = 10 1 1 10 2 = 100 ≡ 45 10^2=100\equiv 45 1 0 2 = 100 ≡ 45 10 × 45 ≡ 10 10\times 45\equiv 10 10 × 45 ≡ 10 2 0 45 2 = 2025 ≡ 45 45^2=2025\equiv 45 4 5 2 = 2025 ≡ 45 不变 3 1 45 2 ≡ 45 45^2\equiv 45 4 5 2 ≡ 45 10 × 45 ≡ 10 10\times 45\equiv 10 10 × 45 ≡ 10 4 1 45 2 ≡ 45 45^2\equiv 45 4 5 2 ≡ 45 10 × 45 ≡ 10 10\times 45\equiv 10 10 × 45 ≡ 10
因此:
10 27 ≡ 10 ( m o d 55 ) 10^{27} \equiv 10 \pmod{55} 1 0 27 ≡ 10 ( mod 55 )
RSA 算法的时间复杂度
RSA 算法的加密与解密主要计算为模指数运算。采用快速幂算法时,其时间复杂度为
O ( log n ) O(\log n) O ( log n )
次模运算。
计算可行性与不可攻破性
可行性
已知 p , q p,q p , q 时,就能轻松算出 φ ( n ) \varphi(n) φ ( n ) ,再用扩展 Euclid 算法求 d d d ;模指数运算可用反复平方算法。因此 RSA 的密钥生成和加解密过程均可在多项式时间内完成,计算上可行。
不可攻破性
已知 n n n 而不知道 p , q p,q p , q 时,无法直接得到 φ ( n ) \varphi(n) φ ( n ) ,也无法求出 d d d 。在现有计算能力和已知算法条件下,大整数因式分解在计算上不可行,因此 RSA 具有不可攻破性。
RSA 计算大题例题
设
p = 3 , q = 11 , m = 30 p=3, \quad q=11, \quad m=30 p = 3 , q = 11 , m = 30
进行 RSA 计算:先求 n n n 、φ ( n ) \varphi(n) φ ( n ) ,确定与 φ ( n ) \varphi(n) φ ( n ) 互素的 e e e ,再计算乘法逆元 d d d ,然后用公钥 ( e , n ) (e,n) ( e , n ) 和私钥 ( d , n ) (d,n) ( d , n ) 对 30 30 30 (明文)进行加密和解密。
参数计算
n = p q = 33 n=pq=33 n = pq = 33
φ ( n ) = ( p − 1 ) ( q − 1 ) = 20 \varphi(n)=(p-1)(q-1)=20 φ ( n ) = ( p − 1 ) ( q − 1 ) = 20
选取
e = 7 e=7 e = 7
因为 1 < e < 20 1<e<20 1 < e < 20 且 gcd ( 7 , 20 ) = 1 \gcd(7,20)=1 g cd( 7 , 20 ) = 1 。
可选的互素数在手稿中记为:3 , 7 , 9 , 11 , 13 , 17 , 19 3,7,9,11,13,17,19 3 , 7 , 9 , 11 , 13 , 17 , 19 。
扩展欧几里得表
为求 e = 7 e=7 e = 7 在模 φ ( n ) = 20 \varphi(n)=20 φ ( n ) = 20 下的逆元,原稿给出了辗转相除表:
No. a = e a=e a = e b = φ ( n ) b=\varphi(n) b = φ ( n ) q = ⌊ a / b ⌋ q=\lfloor a/b \rfloor q = ⌊ a / b ⌋ r r r 1 7 20 0 7 2 20 7 2 6 3 7 6 1 1 4 6 1 6 0 5 1 0
对应等式为:
7 = 20 × 0 + 7 7 = 20 \times 0 + 7 7 = 20 × 0 + 7
20 = 7 × 2 + 6 20 = 7 \times 2 + 6 20 = 7 × 2 + 6
7 = 6 × 1 + 1 7 = 6 \times 1 + 1 7 = 6 × 1 + 1
6 = 1 × 6 + 0 6 = 1 \times 6 + 0 6 = 1 × 6 + 0
因此
gcd ( 7 , 20 ) = 1 \gcd(7,20)=1 g cd( 7 , 20 ) = 1
说明 7 7 7 与 20 20 20 互质,逆元存在。
求逆元 d d d
由回代可得:
20 = 7 × 2 + 6 20 = 7 \times 2 + 6 20 = 7 × 2 + 6
7 = 6 × 1 + 1 7 = 6 \times 1 + 1 7 = 6 × 1 + 1
1 = 7 − 6 1 = 7 - 6 1 = 7 − 6
又因为
6 = 20 − 7 × 2 6 = 20 - 7 \times 2 6 = 20 − 7 × 2
所以
1 = 7 − ( 20 − 7 × 2 ) = 3 × 7 − 20 1 = 7 - (20 - 7 \times 2) = 3 \times 7 - 20 1 = 7 − ( 20 − 7 × 2 ) = 3 × 7 − 20
故
d = 3 d=3 d = 3
于是:
公钥:( 7 , 33 ) (7,33) ( 7 , 33 )
私钥:( 3 , 33 ) (3,33) ( 3 , 33 )
加密表
加密时:
C = M e m o d n = 30 7 m o d 33 C = M^e \bmod n = 30^7 \bmod 33 C = M e mod n = 3 0 7 mod 33
因为
7 = ( 111 ) 2 7=(111)_2 7 = ( 111 ) 2
原稿表中记号说明:表里的 d d d 表示“当前累计结果”,不是私钥中的 d d d 。
取 a = m = 30 a=m=30 a = m = 30 ,n = 33 n=33 n = 33 ,可整理为:
i i i 2 1 0 b i b_i b i 1 1 1 当前平方值 1 30 2 = 900 → 9 30^2=900 \to 9 3 0 2 = 900 → 9 9 2 = 81 → 15 9^2=81 \to 15 9 2 = 81 → 15 平方后取模 1 9 15 当前结果乘以 a a a 30 270 450 乘后再对 n n n 取模 30 6 24
因此
C = 24 C=24 C = 24
解密表
解密时:
M = C d m o d n = 24 3 m o d 33 M = C^d \bmod n = 24^3 \bmod 33 M = C d mod n = 2 4 3 mod 33
此时取 a = 24 a=24 a = 24 ,指数 d = 3 = ( 11 ) 2 d=3=(11)_2 d = 3 = ( 11 ) 2 ,可整理为:
i i i 1 0 b i b_i b i 1 1 当前平方值 1 24 2 = 576 → 15 24^2=576 \to 15 2 4 2 = 576 → 15 当前结果乘以 a a a 24 360 乘后再对 n n n 取模 24 30
因此
M = 30 M=30 M = 30
算法基础
Knuth 算法定义及其特征
算法是一个有限的、确定的、可行的指令序列 ,该序列对给定的输入产生相应的输出,并在有限步骤内终止。
算法的五个特征:
输入
输出
确定性
有限性
有效性
渐近精确界 Θ \Theta Θ 的定义及极限判定准则
在算法分析里,设 f ( n ) f(n) f ( n ) 和 g ( n ) g(n) g ( n ) 是当 n → ∞ n\to\infty n → ∞ 时的非负函数。
大 O O O :上界
若存在正常数 c > 0 , n 0 c>0, n_0 c > 0 , n 0 ,使得当 n ≥ n 0 n\ge n_0 n ≥ n 0 时,
0 ≤ f ( n ) ≤ c g ( n ) 0 \le f(n) \le cg(n) 0 ≤ f ( n ) ≤ c g ( n )
则记作
f ( n ) = O ( g ( n ) ) f(n)=O(g(n)) f ( n ) = O ( g ( n ))
含义:f ( n ) f(n) f ( n ) 的增长速度不超过 g ( n ) g(n) g ( n ) 的常数倍,即至多与 g ( n ) g(n) g ( n ) 同阶,上界为 g ( n ) g(n) g ( n ) 。
大 Ω \Omega Ω :下界
若存在正常数 c > 0 , n 0 c>0, n_0 c > 0 , n 0 ,使得当 n ≥ n 0 n\ge n_0 n ≥ n 0 时,
0 ≤ c g ( n ) ≤ f ( n ) 0 \le cg(n) \le f(n) 0 ≤ c g ( n ) ≤ f ( n )
则记作
f ( n ) = Ω ( g ( n ) ) f(n)=\Omega(g(n)) f ( n ) = Ω ( g ( n ))
含义:f ( n ) f(n) f ( n ) 增长速度至少和 g ( n ) g(n) g ( n ) 一样快,即下界为 g ( n ) g(n) g ( n ) 。
大 Θ \Theta Θ :渐进精确界
若存在正常数 c 1 , c 2 > 0 , n 0 c_1,c_2>0,n_0 c 1 , c 2 > 0 , n 0 ,使得当 n ≥ n 0 n\ge n_0 n ≥ n 0 时,
0 ≤ c 1 g ( n ) ≤ f ( n ) ≤ c 2 g ( n ) 0 \le c_1 g(n) \le f(n) \le c_2 g(n) 0 ≤ c 1 g ( n ) ≤ f ( n ) ≤ c 2 g ( n )
则记作
f ( n ) = Θ ( g ( n ) ) f(n)=\Theta(g(n)) f ( n ) = Θ ( g ( n ))
含义:f ( n ) f(n) f ( n ) 与 g ( n ) g(n) g ( n ) 在无穷远处只差常数倍,增长速度完全同阶。中文常说:与 g ( n ) g(n) g ( n ) 同阶,精确阶为 g ( n ) g(n) g ( n ) 。
极限判定
令
L = lim n → ∞ f ( n ) g ( n ) L=\lim_{n\to\infty}\frac{f(n)}{g(n)} L = n → ∞ lim g ( n ) f ( n )
若 0 < L < ∞ 0<L<\infty 0 < L < ∞ ,则 f ( n ) = Θ ( g ( n ) ) f(n)=\Theta(g(n)) f ( n ) = Θ ( g ( n ))
若 L = 0 L=0 L = 0 ,则 f ( n ) = o ( g ( n ) ) f(n)=o(g(n)) f ( n ) = o ( g ( n )) ,从而 f ( n ) = O ( g ( n ) ) f(n)=O(g(n)) f ( n ) = O ( g ( n ))
若 L = ∞ L=\infty L = ∞ ,则 f ( n ) = ω ( g ( n ) ) f(n)=\omega(g(n)) f ( n ) = ω ( g ( n )) ,从而 f ( n ) = Ω ( g ( n ) ) f(n)=\Omega(g(n)) f ( n ) = Ω ( g ( n ))
特别地,若
lim n → ∞ f ( n ) g ( n ) = c , 0 < c < ∞ \lim_{n\to\infty}\frac{f(n)}{g(n)}=c, \quad 0<c<\infty n → ∞ lim g ( n ) f ( n ) = c , 0 < c < ∞
则
f ( n ) = Θ ( g ( n ) ) f(n)=\Theta(g(n)) f ( n ) = Θ ( g ( n ))
若 c = 1 c=1 c = 1 ,则还可写为
f ( n ) ∼ g ( n ) f(n)\sim g(n) f ( n ) ∼ g ( n )
常见渐近函数的阶
函数 阶 中文描述 1 1 1 O ( 1 ) O(1) O ( 1 ) 常数阶 log n \log n log n O ( log n ) O(\log n) O ( log n ) 对数阶 n n n O ( n ) O(n) O ( n ) 线性阶 n log n n\log n n log n O ( n log n ) O(n\log n) O ( n log n ) 线性对数阶 n 2 n^2 n 2 O ( n 2 ) O(n^2) O ( n 2 ) 平方阶 n 3 n^3 n 3 O ( n 3 ) O(n^3) O ( n 3 ) 立方阶 2 n 2^n 2 n O ( 2 n ) O(2^n) O ( 2 n ) 指数阶 n ! n! n ! O ( n ! ) O(n!) O ( n !) 阶乘阶
大师定理
设 a ≥ 1 a\ge 1 a ≥ 1 和 b > 1 b>1 b > 1 是常数,f ( n ) f(n) f ( n ) 是一个函数,
T ( n ) = a T ( n b ) + f ( n ) T(n)=aT\left(\frac{n}{b}\right)+f(n) T ( n ) = a T ( b n ) + f ( n )
对比 f ( n ) f(n) f ( n ) 与 n log b a n^{\log_b a} n l o g b a :
若
f ( n ) = O ( n log b a − ε ) , ε > 0 f(n)=O\left(n^{\log_b a-\varepsilon}\right), \quad \varepsilon>0 f ( n ) = O ( n l o g b a − ε ) , ε > 0
则
T ( n ) = Θ ( n log b a ) T(n)=\Theta\left(n^{\log_b a}\right) T ( n ) = Θ ( n l o g b a )
若
f ( n ) = Θ ( n log b a ) f(n)=\Theta\left(n^{\log_b a}\right) f ( n ) = Θ ( n l o g b a )
则
T ( n ) = Θ ( n log b a log n ) T(n)=\Theta\left(n^{\log_b a}\log n\right) T ( n ) = Θ ( n l o g b a log n )
若
f ( n ) = Ω ( n log b a + ε ) , ε > 0 f(n)=\Omega\left(n^{\log_b a+\varepsilon}\right), \quad \varepsilon>0 f ( n ) = Ω ( n l o g b a + ε ) , ε > 0
且对充分大的 n n n 满足正则条件
a f ( n / b ) ≤ c f ( n ) , c < 1 af(n/b)\le cf(n), \quad c<1 a f ( n / b ) ≤ c f ( n ) , c < 1
则
T ( n ) = Θ ( f ( n ) ) T(n)=\Theta(f(n)) T ( n ) = Θ ( f ( n ))
例 1
T ( n ) = 9 T ( n / 3 ) + n T(n)=9T(n/3)+n T ( n ) = 9 T ( n /3 ) + n
有 a = 9 , b = 3 , f ( n ) = n a=9,b=3,f(n)=n a = 9 , b = 3 , f ( n ) = n ,而
n log 3 9 = n 2 n^{\log_3 9}=n^2 n l o g 3 9 = n 2
所以
T ( n ) = Θ ( n 2 ) T(n)=\Theta(n^2) T ( n ) = Θ ( n 2 )
例 2
T ( n ) = T ( 2 3 n ) + 1 T(n)=T\left(\frac{2}{3}n\right)+1 T ( n ) = T ( 3 2 n ) + 1
手稿中写法等价于 a = 1 , b = 3 / 2 , f ( n ) = 1 a=1,b=3/2,f(n)=1 a = 1 , b = 3/2 , f ( n ) = 1 ,因此
T ( n ) = Θ ( log n ) T(n)=\Theta(\log n) T ( n ) = Θ ( log n )
例 3
T ( n ) = 3 T ( n / 4 ) + n log n T(n)=3T(n/4)+n\log n T ( n ) = 3 T ( n /4 ) + n log n
其中
n log 4 3 ≈ O ( n 0.793 ) n^{\log_4 3}\approx O(n^{0.793}) n l o g 4 3 ≈ O ( n 0.793 )
而 n log n n\log n n log n 多项式意义上更大,并满足正则条件,因此
T ( n ) = Θ ( n log n ) T(n)=\Theta(n\log n) T ( n ) = Θ ( n log n )
例 4
T ( n ) = 2 T ( n / 2 ) + n log n T(n)=2T(n/2)+n\log n T ( n ) = 2 T ( n /2 ) + n log n
此时 n log 2 2 = n n^{\log_2 2}=n n l o g 2 2 = n ,而 f ( n ) = n log n f(n)=n\log n f ( n ) = n log n 比 n n n 大一个对数因子,不属于标准三种情形中的多项式差距,手稿注明:不适用于基础版大师定理 。
大 O O O 记法定义及其极限判别准则
对于给定的复杂函数 T ( n ) T(n) T ( n ) ,若存在两个正常数 C C C 和 n 0 n_0 n 0 ,以及函数 τ ( n ) \tau(n) τ ( n ) ,使得当 n ≥ n 0 n\ge n_0 n ≥ n 0 时有
T ( n ) ≤ C τ ( n ) T(n)\le C\tau(n) T ( n ) ≤ C τ ( n )
则称
T ( n ) = O ( τ ( n ) ) T(n)=O(\tau(n)) T ( n ) = O ( τ ( n ))
并且
T ( n ) = O ( τ ( n ) ) ⟺ 0 ≤ lim n → ∞ T ( n ) τ ( n ) < + ∞ T(n)=O(\tau(n)) \iff 0 \le \lim_{n\to\infty}\frac{T(n)}{\tau(n)} < +\infty T ( n ) = O ( τ ( n )) ⟺ 0 ≤ n → ∞ lim τ ( n ) T ( n ) < + ∞
NP、NPH、NPC
NP :答案如果给你,能很快验证对不对
NP-Hard :至少和 NP 中最难的问题一样难
NP-Complete :既是 NP,又是 NP-Hard
记法:
N P C = N P ∩ N P H NPC = NP \cap NPH N P C = N P ∩ N P H
判定式定义
NP 问题 :如果一个判定问题的解可以在多项式时间内验证,则属于 NP。
NPH 问题 :如果一个问题 π \pi π 满足:对任意 π ′ ∈ N P \pi'\in NP π ′ ∈ N P ,都有 π ′ ≤ p π \pi'\le_p \pi π ′ ≤ p π ,则称 π \pi π 为 NPH。也就是 NP 中任何问题都可以多项式时间归约到它,因此它至少和所有 NP 问题一样难。
NPC 问题 :如果一个问题 π \pi π 满足:
π ∈ N P \pi \in NP π ∈ N P
对任意 π ′ ∈ N P \pi'\in NP π ′ ∈ N P ,都有 π ′ ≤ p π \pi'\le_p \pi π ′ ≤ p π
则称 π \pi π 为 NPC。
也就是说,它的答案能快速验证,属于 NP;又足够难,属于 NPH,所以是 NP 里最难的一批问题。
若能解决问题 π \pi π ,那么把其他问题多项式时间转换为 π \pi π 的过程,也就能一并解决。
经典归约链
手稿中的归约链为:
布尔电路可满足性问题 → S A T → 3 - S A T \text{布尔电路可满足性问题} \to SAT \to 3\text{-}SAT 布尔电路可满足性问题 → S A T → 3 - S A T
再由 3 3 3 -SAT 归约到其他经典问题,例如:
团问题 → \to → 顶点覆盖问题 → \to → 子集和 / 覆盖类问题(手稿此处为示意)
哈密尔顿回路问题 → \to → 旅行商问题(TSP)
分治排序算法
Karatsuba 大数乘法
把两个长整数 x , y x,y x , y 折半:
x = x 1 B m + x 0 x=x_1B^m+x_0 x = x 1 B m + x 0
y = y 1 B m + y 0 y=y_1B^m+y_0 y = y 1 B m + y 0
则
x y = ( x 1 B m + x 0 ) ( y 1 B m + y 0 ) xy=(x_1B^m+x_0)(y_1B^m+y_0) x y = ( x 1 B m + x 0 ) ( y 1 B m + y 0 )
展开为
x y = x 1 y 1 B 2 m + ( x 1 y 0 + x 0 y 1 ) B m + x 0 y 0 xy=x_1y_1B^{2m}+(x_1y_0+x_0y_1)B^m+x_0y_0 x y = x 1 y 1 B 2 m + ( x 1 y 0 + x 0 y 1 ) B m + x 0 y 0
记
z 2 = x 1 y 1 , z 1 = x 1 y 0 + x 0 y 1 , z 0 = x 0 y 0 z_2=x_1y_1, \qquad z_1=x_1y_0+x_0y_1, \qquad z_0=x_0y_0 z 2 = x 1 y 1 , z 1 = x 1 y 0 + x 0 y 1 , z 0 = x 0 y 0
则仍是四次乘法。
进一步将
z 1 = ( x 1 + x 0 ) ( y 1 + y 0 ) − x 1 y 1 − x 0 y 0 = z 3 − z 2 − z 0 z_1=(x_1+x_0)(y_1+y_0)-x_1y_1-x_0y_0 = z_3-z_2-z_0 z 1 = ( x 1 + x 0 ) ( y 1 + y 0 ) − x 1 y 1 − x 0 y 0 = z 3 − z 2 − z 0
只需要三次乘法。
递推式:
T ( n ) = 3 T ( n / 2 ) + c n + d T(n)=3T(n/2)+cn+d T ( n ) = 3 T ( n /2 ) + c n + d
因此
T ( n ) = Θ ( n log 2 3 ) ≈ Θ ( n 1.585 ) T(n)=\Theta(n^{\log_2 3})\approx \Theta(n^{1.585}) T ( n ) = Θ ( n l o g 2 3 ) ≈ Θ ( n 1.585 )
伪代码
输入:两个 n 位正整数 X 和 Y
输出:X 和 Y 的乘积
1. multiplyK(X, Y)
2. n <- X 和 Y 的位数
3. if n = 1 return X * Y
4. x1, x0 <- X 的高 n/2 位和低 n/2 位
5. y1, y0 <- Y 的高 n/2 位和低 n/2 位
6. z2 <- multiplyK(x1, y1)
7. z0 <- multiplyK(x0, y0)
8. z3 <- multiplyK(x1 + x0, y1 + y0)
9. return z2 * 10^n + (z3 - z2 - z0) * 10^(n/2) + z0
二叉最大堆
给定序列:
[ 40 , 55 , 49 , 73 , 12 , 27 , 98 , 81 , 64 , 36 ] [40,55,49,73,12,27,98,81,64,36] [ 40 , 55 , 49 , 73 , 12 , 27 , 98 , 81 , 64 , 36 ]
要求:给出由该序列构建大根堆的过程;根节点比子节点大。
从 ⌊ n / 2 ⌋ = 10 / 2 = 5 \lfloor n/2 \rfloor = 10/2 = 5 ⌊ n /2 ⌋ = 10/2 = 5 开始,自底向上调整。
初始完全二叉树
graph TD
A40 (( 40 )) --> A55 (( 55 ))
A40 --> A49 (( 49 ))
A55 --> A73 (( 73 ))
A55 --> A12 (( 12 ))
A49 --> A27 (( 27 ))
A49 --> A98 (( 98 ))
A73 --> A81 (( 81 ))
A73 --> A64 (( 64 ))
A12 --> A36 (( 36 ))
调整后的大根堆
graph TD
R98 (( 98 )) --> R81 (( 81 ))
R98 --> R49 (( 49 ))
R81 --> R73 (( 73 ))
R81 --> R36 (( 36 ))
R49 --> R27 (( 27 ))
R49 --> R40 (( 40 ))
R73 --> R55 (( 55 ))
R73 --> R64 (( 64 ))
R36 --> R12 (( 12 ))
小根堆中依次插入数据 4 , 2 , 5 , 8 , 3 , 6 , 10 , 1 4, 2, 5, 8, 3, 6, 10, 1 4 , 2 , 5 , 8 , 3 , 6 , 10 , 1
插入 4
graph TD
A4 (( 4 ))
插入 2
graph TD
B2 (( 2 )) --> B4 (( 4 ))
插入 5
graph TD
C2 (( 2 )) --> C4 (( 4 ))
C2 --> C5 (( 5 ))
插入 8
graph TD
D2 (( 2 )) --> D4 (( 4 ))
D2 --> D5 (( 5 ))
D4 --> D8 (( 8 ))
插入 3
graph TD
E2 (( 2 )) --> E3 (( 3 ))
E2 --> E5 (( 5 ))
E3 --> E8 (( 8 ))
E3 --> E4 (( 4 ))
插入 6
graph TD
F2 (( 2 )) --> F3 (( 3 ))
F2 --> F5 (( 5 ))
F3 --> F8 (( 8 ))
F3 --> F4 (( 4 ))
F5 --> F6 (( 6 ))
插入 10
graph TD
G2 (( 2 )) --> G3 (( 3 ))
G2 --> G5 (( 5 ))
G3 --> G8 (( 8 ))
G3 --> G4 (( 4 ))
G5 --> G6 (( 6 ))
G5 --> G10 (( 10 ))
插入 1(最终)
graph TD
H1 (( 1 )) --> H2 (( 2 ))
H1 --> H5 (( 5 ))
H2 --> H3 (( 3 ))
H2 --> H4 (( 4 ))
H3 --> H8 (( 8 ))
H5 --> H6 (( 6 ))
H5 --> H10 (( 10 ))
回溯法 VS 分支限界法:N 皇后问题
回溯法和分支限界法通过约束条件和限界条件对状态空间树进行剪枝,使得一部分分支不需搜索,从而对很多问题实例都能改进穷举法的搜索效率。
回溯法:在 DFS 上进行
分支限界法:在 BFS 基础上进行
N 皇后问题中,回溯法通过判断两个皇后不能在同一列上,即 x [ i ] ≠ x [ j ] x[i] \ne x[j] x [ i ] = x [ j ] ,并且不能在同一条对角线上,即
∣ i − j ∣ ≠ ∣ x [ i ] − x [ j ] ∣ |i-j| \ne |x[i]-x[j]| ∣ i − j ∣ = ∣ x [ i ] − x [ j ] ∣
来实现剪枝,从而改进穷举算法。
手稿示例输入输出:
n = 6 n=6 n = 6
输出三组结果:
并注明:第一个位于第二列,第二行位于第四列。
动态规划与贪心算法
0/1 背包和动态规划
动态规划适用于:
最优子结构 :原问题最优解可以由子问题最优解构成
重叠子问题 :在递归求解过程中,不同子问题会重复出现
贪心与动态规划的不同
贪心
局部最优
不回溯
子问题一般不重叠
结果不一定最优
动态规划
0/1 背包状态定义
令第 i i i 个物品重量为 w i w_i w i ,价值为 v i v_i v i 。
设
d p [ i ] [ j ] dp[i][j] d p [ i ] [ j ]
表示前 i i i 个物品、背包容量为 j j j 时的最大价值。
状态转移:
d p [ i ] [ j ] = { d p [ i − 1 ] [ j ] , j < w i max { d p [ i − 1 ] [ j ] , d p [ i − 1 ] [ j − w i ] + v i } , j ≥ w i dp[i][j] =
\begin{cases}
dp[i-1][j], & j < w_i \\
\max\{dp[i-1][j],\; dp[i-1][j-w_i]+v_i\}, & j \ge w_i
\end{cases} d p [ i ] [ j ] = { d p [ i − 1 ] [ j ] , max { d p [ i − 1 ] [ j ] , d p [ i − 1 ] [ j − w i ] + v i } , j < w i j ≥ w i
手稿右下角还写有一维压缩式的提示,形式可写为:
d p [ i ] [ j ] = max ( d p [ i − 1 ] [ j ] , d p [ i − 1 ] [ j − w i ] + v i ) dp[i][j] = \max\bigl(dp[i-1][j],\; dp[i-1][j-w_i] + v_i\bigr) d p [ i ] [ j ] = max ( d p [ i − 1 ] [ j ] , d p [ i − 1 ] [ j − w i ] + v i )
Levenshtein 编辑距离与矩阵链相乘
两个字符串间的编辑距离,是指把一个字符串变成另一个字符串的基本操作数。这里的基本编辑操作包括:
递推定义
设 a , b a,b a , b 为两个字符串,l e v a , b ( i , j ) lev_{a,b}(i,j) l e v a , b ( i , j ) 表示:字符串 a a a 的前 i i i 个字符变成字符串 b b b 的前 j j j 个字符的最小编辑距离。
边界:
min ( i , j ) = 0 \min(i,j)=0 min ( i , j ) = 0
更准确地说:
若 i = 0 i=0 i = 0 ,表示空串到 b b b 的前 j j j 个字符,需要连续插入 j j j 次
若 j = 0 j=0 j = 0 ,表示 a a a 的前 i i i 个字符到空串,需要连续删除 i i i 次
定义代价函数:
δ ( a i , b j ) = { 0 , a i = b j 1 , a i ≠ b j \delta(a_i,b_j)=
\begin{cases}
0, & a_i=b_j \\
1, & a_i\ne b_j
\end{cases} δ ( a i , b j ) = { 0 , 1 , a i = b j a i = b j
则递推式为:
l e v a , b ( i , j ) = { max ( i , j ) , min ( i , j ) = 0 min { l e v a , b ( i − 1 , j ) + 1 , l e v a , b ( i , j − 1 ) + 1 , l e v a , b ( i − 1 , j − 1 ) + δ ( a i , b j ) } , otherwise lev_{a,b}(i,j)=
\begin{cases}
\max(i,j), & \min(i,j)=0 \\
\min\left\{
lev_{a,b}(i-1,j)+1,\;
lev_{a,b}(i,j-1)+1,\;
lev_{a,b}(i-1,j-1)+\delta(a_i,b_j)
\right\}, & \text{otherwise}
\end{cases} l e v a , b ( i , j ) = { max ( i , j ) , min { l e v a , b ( i − 1 , j ) + 1 , l e v a , b ( i , j − 1 ) + 1 , l e v a , b ( i − 1 , j − 1 ) + δ ( a i , b j ) } , min ( i , j ) = 0 otherwise
三项分别对应:
l e v ( i − 1 , j ) + 1 lev(i-1,j)+1 l e v ( i − 1 , j ) + 1 :删除 a i a_i a i
l e v ( i , j − 1 ) + 1 lev(i,j-1)+1 l e v ( i , j − 1 ) + 1 :插入 b j b_j b j
l e v ( i − 1 , j − 1 ) + δ ( a i , b j ) lev(i-1,j-1)+\delta(a_i,b_j) l e v ( i − 1 , j − 1 ) + δ ( a i , b j ) :若相等则不改,代价 0;否则替换,代价 1
伪代码
Function LevenshteinDistance(S1, S2)
m = length(S1)
n = length(S2)
Create DP[m+1][n+1]
// 初始化
For i = 0 to m:
DP[i][0] = i
For j = 0 to n:
DP[0][j] = j
// 递推
For i = 1 to m:
For j = 1 to n:
If S1[i-1] == S2[j-1]
cost = 0
else
cost = 1
DP[i][j] = min(
DP[i-1][j] + 1, // 删除
DP[i][j-1] + 1, // 插入
DP[i-1][j-1] + cost // 替换
)
return DP[m][n]