南邮 OJ 1037 模运算

模运算

时间限制(普通/Java) :  1000 MS/ 3000 MS          运行内存限制 : 65536 KByte
总提交 : 1092            测试通过 : 202 

比赛描述

给定整数a,b,n,要求计算(a^b)mod n

输入

多组数据,每组数据一行,为三个用空格隔开的整数a,b,n

1<=a<=5000,0<=b<=10^8,1<=n<=5000000

输出

每组数据输出一行,为所求值

样例输入

2 3 5
2 2 4

样例输出

3
0

题目来源

NUAA

#include <iostream>
using namespace std;
/*
*功能：快速幂取模
       a*b%c=( (a%c)*b )%c
       b = p(n)*2^n+p(n-1)*2^(n-1)+…+p(1)*2+p(0)
       a^b = a^(p(n)*2^n+p(n-1)*2^(n-1)+...+p(1)*2+p(0))
           = a^(p(n)*2^n)* a^(p(n-1)*2^(n-1))*...*a^(p(1)*2)*a^p(0)
	   用temp暂存a^(p(i)*2^i)%n
	   a^(2^i)%c = ( (a^(2^(i-1))%c)^2 )  %c
*输入：a,b,n
*输出：a*b%c
*作者：陈汝军
*时间：2014-8-24 16:49:38
*/
__int64 expMod(__int64 a,__int64 b,__int64 n){
	__int64 temp=a%n;
	__int64 ret=1;
	while(b){
		if(b&1)
			ret = ret*temp%n;
		temp = temp*temp%n;
		b >>= 1;
	}
	return ret;
}

int main(void){
	__int64 a=0,b=0,n=0;
	while(cin>>a>>b>>n){
		cout<<expMod(a,b,n)<<endl;	
	}
	return 0;
}