模运算------x的y次方再模n

问题描述：ax mod n （a的x次方再模n）求结果，对小数量的求模一般数据类型可以应付，但是当x很大时就会产生很大的中间结果。我们希望加速运算并且不产生大结果。

From: http://blog.csdn.net/sqx2011/article/details/8241713

我们知道模运算满足分配率，（a*b）mod n = ((a mod n)*(b mod n)) mod n,所以我们将指数运算同模运算结合起来，将指数运算作为一系列的乘法运算，对于每次指数运算做一次模运算，（a*a*a*a*a*a*）mod n = (((a*a)mod n)((a*a)mod n)((a*a)mod n)*(x mod n))mod n.代码实现如下：

[cpp] view plain copy

1. #include <iostream>  
2. using namespace std;  
3.   
4. int pow(int x, int y, int mod)  
5. {  
6.     int res = 1;  
7.     while(y)  
8.     {  
9.         if(y & 1) res =  res * x % mod;  
10.         x = x * x % mod;  
11.         y = y / 2;  
12.     }  
13.     return res;  
14. }  
15.   
16. int main()  
17. {  
18.     int x, y, z;  
19.     while(1){  
20.     cout << "请输入底数 指数 模，空格隔开：" ;  
21.     cin >> x >> y >> z;  
22.     cout <<"结果为：" << pow(x, y, z) << endl;  
23.     }  
24.     return 0;  
25. }  

其中（y & 1）判断mod n 的个数是否是奇数，如果为1说明有单独的x mod n 没有进行指数运算,计算它保存到临时res,并且到最后只剩一个指数运算的时候将结果和临时res相乘取模。如果y是2倍数就不存在这种情况了。