#include <iostream>
using namespace std;
const int MOD = 1000;//像这样的一个常量就应该专门定义一下
int PowMod(int a, int n)//a^n%MOD
 {
       int ret = 1;
       while(n) {
             if(n & 1) ret = ret * a % MOD;
             //变为二进制,然后就可以专门进行分解计算,十分的方便,要求是结合位运算一同使用
             a = a * a % MOD;
             //这里要求特别的注意,因为是二进制,所以每个步骤都是平方次处理,如果是三进制就会发现要求有^3的形式了,所以在这里要求平方
             n >>= 1;
       }
             return ret;
 } 	//可以将上面的看成一个专门的模板
int main()
{
    int a,b;
     cin>>a>>b;
    while(a!=0 && b!=0||a!=0 && b==0||a==0 && b!=0)
    {
        cout<<PowMod(a,b)<<endl;
         cin>>a>>b;
}
    return 0;
}
//位运算的基本操作已经基本搞懂了,后期就是处理问题和自己解决问题的时间,当然,这个是最基本的一道题
矩阵快速幂AC代码HDU 2035,布布扣,bubuko.com
原文地址:http://www.cnblogs.com/tianxia2s/p/3853025.html