cnspp-ch04.pptVIP

* Compute successive instances of GCD(a,b) = GCD(b,a mod b). Note this MUST always terminate since will eventually get a mod b = 0 (ie no remainder left). * Stallings Table 4.6 * See text for additional discussion. Reduction comes from observation that if in GF(2n) then irreducible poly g(x) has highest term xn , and if compute xn mod g(x) answer is g(x)- xn 如果gcd(a, n)=1, 一定存在整数x, 0xn，满足 ax mod n=1 构造乘法表 例如: 可以用计算最大公约数的Euclid算法的扩展来求逆。 在GF(p)中求乘法逆元 例: 求gcd(1970, 1066) 1970 = 1 x 1066 + 904 gcd(1066, 904) 1066 = 1 x 904 + 162 gcd(904, 162) 904 = 5 x 162 + 94 gcd(162, 94) 162 = 1 x 94 + 68 gcd(94, 68) 94 = 1 x 68 + 26 gcd(68, 26) 68 = 2 x 26 + 16 gcd(26, 16) 26 = 1 x 16 + 10 gcd(16, 10) 16 = 1 x 10 + 6 gcd(10, 6) 10 = 1 x 6 + 4 gcd(6, 4) 6 = 1 x 4 + 2 gcd(4, 2) 4 = 2 x 2 + 0 gcd(2, 0) 因此Gcd(1970,1066)=2 在GF(p)中求乘法逆元 如果gcd(m, b)=1，可使用扩展欧几里德算法求出b的乘法逆元 EXTENDED EUCLID(m, b) 1.(A1, A2, A3)=(1, 0, m); (B1, B2, B3)=(0, 1, b) 2. if B3 = 0 return A3 = gcd(m, b); no inverse 3. if B3 = 1 return B3 = gcd(m, b); B2 = b–1 mod m 4. Q = A3 div B3 5. (T1, T2, T3)=(A1 – Q B1, A2 – Q B2, A3 – Q B3) 6. (A1, A2, A3)=(B1, B2, B3) 7. (B1, B2, B3)=(T1, T2, T3) 8. goto 2 多项式运算 三种多项式运算 使用代数基本规则的普通多项式运算 系数运算是模p运算的多项式运算，即系数在GF(p)中 系数在GF(p)中，且多项式被定义为模一个n次多项式m(x)的多项式运算 普通多项式运算 普通多项式运算 一个n次多项式(n=0)的表达形式如下 f(x) = anxn + an-1xn-1 + … + a1x + a0 = ∑ aixi 其中ai是某个指定数集S中的元素，该数集称为系数集，且an≠ 0，f(x)是定义在系数集S上的多项式 零次多项式称为常数多项式，是系数集里的一个元素 如果an=1，对应的n次多项式就称为首1多项式 普通多项式运算 加或减就是相应系数的加减，乘除则要用到所有系数 Eg. 系数在 GF(p)中的多项式运算 在计算每个系数的值时需要做模运算 可以模任何素数p 模2的运算, 系数在 GF(2)中的多项式运算 也就是说所有的系数不是0就是1 比如，令 f(x) = x3 + x2 , g(x) = x2 + x + 1 则 f(x) + g(x) = x3 + x + 1 f(x) - g(x) = x3 + x + 1 f(x) g(x) = x5 + x2 f(x) /g(x) =x……x 既约多项式 多项式可以写成如下形式: f(x) = q(x) m(x) + r(x) 其中，r(x)就可被看作是余数 r(x) = f(x) mod m(x)