快速幂

#include<iostream>
#include<cstdio>
using namespace std;
//公式 
//(ab)mod c = ((a mod c) * (b mod c))mod c
//(a^b) mod c = ((a mod c) ^ b ) mod c
int ggt(int x,int y,int z)
{
	/*利用二进制进行， 
		对每一个数思考他对应的二进制，则ai = 1 || 0 
		 y = a1 + a2 *2^1 + a3 *2^2 + ... + an * 2^(n-1)
		 则 x^y = x^(a1 + a2 *2^1 + a3 *2^2 + ... + an * 2^(n-1))
		 =x^(a1)*x^(a2 *2^1)*...*x^(an *2^n-1)
		 设 A1 = x^(a1) A2 =  x^(a2 *2^1)  --> 可知 An = An-1 * 2 
		 那么对上面的每一项思考，如果等于零，则不运算 
	*/ 
	int ans = 1;
	x = x % z;
	while(y != 0)
	{
		if(y&1) ans = (ans*x)%z;//判断当前位二进制是否为零 
		y >>= 1; //二进制向右移一位，即丢弃一位 
		x = (x*x)%z;//计算An+1 
	}
	return ans;
}
int pm(int x,int y,int z)
{
	/*
		因为a^b = (a^2)^ b/2 
		依次进行循环，直到 b == 0 
		计算以x为主，对x进行累乘 
	*/ 
	int ans = 1;
	x = x % z;
	while(y != 0)
	{
		if(y%2==1) ans = (ans * x) % z;//如果y是奇数，则多出来一位，单独乘 
		y = y / 2;//递推 
		x = (x * x)%z;//一层一层计算 
	}
	return ans;
}
int main()
{
	int a , b , c;
	c = 1000; //取模运算 
	while(scanf("%d %d",&a,&b))
	{
		if(a==0&&b==0) break;
		int sum , sum_1 ,sum_0;
		sum = ggt(a,b,c);
		sum_0 = pm(a,b,c);
		printf("%d\n",sum);
	}
	return 0;
 }