演算法93试除因数分解法的虚拟码范例930.PPTVIP

9.3 因數分解 因數分解（factorization）的問題從過去一直被研究至今，而這樣的研究在未來似乎會持續下去。因數分解在公開金鑰密碼系統（參見第十章）的安全性扮演一個非常重要的角色。 本節討論主題 算術的基本定理 因數分解的方法 費瑪分解法 Pollard p – 1分解法 Pollard rho分解法 更快速的分解法 9.3.1 算術的基本定理 最大公因數 最小公倍數 演算法 9.3 試除因數分解法的虛擬碼 範例9.30 使用試除因數法來找出 1523357784 的因數。 解法：我們執行一個程式並得到以下的結果。 9.3.3 費瑪分解法 費瑪因數分解法（Fermat factorization method，演算法 9.4）可以將一個數值 n 分解成 兩個正整數 a 和 b（不一定為質數），使得 n = a × b。 演算法 9.4 費瑪因數分解法的虛擬碼 9.3.4 Pollard p – 1 分解法 此方法植基於一個條 件：p – 1 不能有大於預定邊界 B 的因數。Pollard 證明在這種條件之下： 演算法 9.5 Pollard p – 1 因數分解法的虛擬碼 範例9.31 使用 Pollard p – 1 分解法來求出的因數，其中邊界 B = 8。 解法：我們依這個演算法執行程式，並找出 p = 421。事實上= 421 × 135979。注意421 是質數，而且 p – 1 沒有任何大於 8 的因數（421 – 1 = 22 × 3 × 5 × 7）。 圖 9.3 Pollard rho 的連續數值 演算法 9.6 Pollard rho 分解法的虛擬碼 範例9.32 假設有一台電腦 1 秒可以執行 230（將近 10 億）次位元運算，其分解下列大小的整數 大約需要多少時間？ a. 60 位數。 b. 100 位數。 範例9.32 (續) 解法： a. 一個 60 位數將近有 200 個位元，其複雜度為 2nb / 4 或 250。當每秒可以執行 230 次運算時，此演算法可以在 220 秒完成，大約為 12 天。 b. 一個 100 位數將近有 300 個位元，其複雜度為 275。當每秒可以執行 230 次運算時，此演算法可以在 245 秒完成，需要好幾年。 範例9.33 我們寫出一個程式來分解 434617，其結果為 709（434617 = 709 × 613）。表 9.2 顯 示執行此程式時的數對（x 和 y）以及 p 的值。 表 9.2 9.3.6 更快速的分解法 二次篩選法 這個方法 使用篩選方式來找出 x2 mod n 的值。 數體篩選法 這個方法在代數的環結構中使用篩選方式來找出 x2 ? y2 mod n。 範例9.34 假設有一台電腦 1 秒可以執行 230（將近 10 億）次位元運算。使用下列方法來分解一 個 100 位數的整數，大約各需要多少時間？ a. 二次篩選法。 b. 數體篩選法。 解法：一個 100 位數將近有 300 個位元（n = 2300）。ln(2300) = 207，而 lnln (2300) = 5。 9.4 中國餘數定理 中國餘數定理（Chinese remainder theorem, CRT）是用來求解模數兩兩相異且互質之單變數的同餘方程組。此方程組如下所述： 範例9.35 以下範例為模數均相異的同餘方程組： 我們將在下一節說明這個方程組的解法；目前，先注意到這個方程組的解為 x = 23。這個值可以滿 足所有的方程式：23 ≡ 2 (mod 3)、23 ≡ 3 (mod 5) 和 23 ≡ 2 (mod 7)。 範例9.35 (續) 解法：我們可以依照以下的步驟來解這個方程組： 求出 M = m1 × m2 × … × mk。這是所有方程式共同的模數。 求出 M1 = M/m1, M2 = M/m2, …, Mk = M/mk。 在相對應的模數 (m1, m2, …, mk) 下，求出 M1, M2, …, Mk的乘法反元素 M1?1, M2?1, …, Mk ?1 。 這些方程式共同的解為： 範例9.38 假設我們需要計算 z = x + y，其