快速幂求模

C - 人见人爱A^B

Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u

Submit  Status

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 
     
     

 

这道题A,B可以取到10000，暴力运算会溢出。可以采用快速幂求模运算。求最后三位数，即对1000求模：

#include<iostream>
#include<algorithm>
using namespace std;

int main()
{
	int ans,a,b;
	while(cin>>a>>b)
	{
		if(a==0&&b==0) break;
		ans=1;
		while(b>0)
		{
			if(b % 2 == 1)
			{
				ans = (ans * a) % 1000;
				b=b-1;
			}
			if(b%2==0)
			{
				b = b/2;
				a = (a*a)%1000;
			}
			
		}
		cout<<ans<<endl;
	}
	
	return 0;
}

比如。对2的5次方对10求模。

将2的5次方拆成2的4次方与2相乘，分别对两边求模。

对2求模，可以得到ans=2；对2的4次方求模，根据程序，会变成：a=（2*2）%10=4,b=4/2=2 ；a=（4*4）%10=6.b=2/2=1;

这时b=1，可以得到ans=(ans*a)%10=2;程序结束。