快速指数算法 和 求逆元算法

源地址：http://www.cnblogs.com/yinghuochong/archive/2011/10/08/2203107.html

快速指数算法 和 求逆元算法

快速指数算法 和 求逆元 的算法是加密中常用到的两个算法。

这两个算法主要都是涉及到的 模运算 ，对于模运算的性质总结如下：

（1）     (a + b) % n = (a % n + b % n) % n

（2）     (a - b) % n = (a % n - b % n) % n 
（3）     (a * b) % n = (a % n * b % n) % n 
（4）      a^b % n = ((a % n)^b) % n 
（1）式证明
∵ a = k1*n + r1

b = k2*n + r2

a % n = r1

b % n = r2

∴(a+b) % n = ((k1+k2)*n + (r1+r2)) % n = (r1+r2) % n = (a % n + b % n)% n
得证
（2）式证明同上
（3）式证明
a = k1*n + r1
b = k2*n + r2
(a*b) % n = (k1k2n² + (k1r2+k2r1)n + r1r2) % n = r1r2 % n = (a %n * b %n ) % n
(4)式证明
设 a % n = r
a^b %n= (a * a * a * a…*a) %n = (a %n * a %n * a %n * … * a %n) %n = r^b % n = ((a % n) ^b) % n

 

快速指数算法： 要计算x的e次方对m取余的值，如果直接计算会数值会比较大，可以利用上面的模的性质进行降幂计算。

要计算 x^e%m 的值 如 62⁶⁵% 133  可以用如下的方法：

62⁶⁵ % 133
= 62 * 62⁶⁴ % 133
= 62 * (62²)³² % 133
= 62 * 3844³² % 133
= 62 * (3844 % 133)³² % 133
= 62 * 120³²% 133

= 62 * 36¹⁶ % 133
= 62 * 99⁸ % 133
= 62 * 92⁴ % 133
= 62 * 85² % 133
= 62 * 43 % 133
= 2666 % 133
= 6

C# 中实现如下：

输入62、65、133测试数据输出结果如上图。

 

求逆元算法：

定义 如果ab≡1(mod m), 则称b是a的模m逆,记作a的模m逆是方程ax≡1(mod m)的解. 两个数互质一定有逆元。

求逆元可以使用辗转相除法，但是只有两个数都是质数的时候才有逆元，举例如下：

例：求5的模7逆

做辗转相除法, 求得整数b,k使得 5b+7k=1, 则b是5的模7逆.

计算如下:

     7=5+2,  5=2×2+1.

回代    1=5-2×2=5-2×(7-5)= 3×5-2×7,

得 5 ^-1≡3(mod7).

例：求21的模73逆

做辗转相除法, 求得整数b,k使得 21b+73k=1, 则b是21的模73逆.

计算如下:

     73=21*3+10

21=10*2+1

回代    1=21-10*2

1=21-（73-21*3）*2

=21-73*2+6*21

=7*21-73*2

得 21 ^-1≡7(mod73).

例：求7的模96逆

做辗转相除法, 求得整数b,k使得 7b+96k=1, 则b是7的模96逆.

计算如下:

96=7*13+5

7=5*1+2

5=2*2+1 

回代 1=5-2*2

1=5-（7-5*1）*2

=5*3-7*2

=(96-7*13)*3-7*2

=96*3-41*7

-41 mod 96 =55 所以55就是7关于96的逆元。

得 7 ^-1≡55(mod96).

 

C#中实现如下：

测试程序如下：

?

int num, modeNum;

           while (true)

           {

               num = Int32.Parse(Console.ReadLine());

               modeNum = Int32.Parse(Console.ReadLine());

               int gcd=Gcd(num, modeNum);

               //不互为质数，有公约数

               if (gcd!= 1)

               {                   

                   Console.Write("有公约数，公约数为：" + gcd.ToString() + "\n");

                   Console.ReadKey();

                   continue;

               }

               int opposite = GetOpposite(num, modeNum);

               Console.Write("opposite=" + opposite.ToString() + "\n");

               Console.ReadKey();

           }

　　

测试结果如下：

分类: 算法与数据结构