HDU 2035人见人爱A^B(数字快速幂法)

人见人爱A^B

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 21906    Accepted Submission(s): 15264

Problem Description

求A^B的最后三位数表示的整数。
说明：A^B的含义是“A的B次方”

 

Input

输入数据包含多个测试实例，每个实例占一行，由两个正整数A和B组成（1<=A,B<=10000），如果A=0, B=0，则表示输入数据的结束，不做处理。

 

Output

对于每个测试实例，请输出A^B的最后三位表示的整数，每个输出占一行。

 

Sample Input

  
  

   
   2 3
12 6
6789 10000
0 0
  
  

 

Sample Output

  
  

   
   8
984
1
  
  

 

/*假如题目中的b由上限10^4变为10^9，则原来的暴力法会超时

但是由于10^9小于2^30，我们可以利用递推法来求出*/
/*我们先将B转换为二进制整数
例如B=(45)10=(101101)2
也就意味着45=2^5+2^3+2^2+2^0
                           =32+8+4+1


那么x^45 = x^32 * x^8 * x^4 * x^1


计算出x,x^2,x^4,x^8 …… 很快吧。
就算B有10^9这么大，Log2(10^9)<30
不超过30次就搞定了。
*/
#include<iostream>
using namespace std;
const int maxn=32;
const int MOD=1000;
int data[maxn];
int a,b,ans;
int main()
{
int i;
while(cin>>a>>b)
{
if(a==0&&b==0){break;}
data[0]=a;//数组的第一项应该为a，与b=1时相呼应
for(i=1;i<maxn;i++)
{data[i]=data[i-1]*data[i-1]%MOD;}/*把a的2,4,8...次幂的结果依次算出来*/
ans=1;
for(i=0;i<maxn;i++)
{
if(b&1)/*如果这一位存在的话要乘上*/
{ans=ans*data[i]%MOD;}
b>>=1;/*无论是否存在这一位，b都要右移一位*/
}
cout<<ans<<endl;
}
return 0;
}