A^B mod C(高速べき乗+高速乗+型取り)問題解
1192 ワード
A^B mod C
Given A,B,C, You should quickly calculate the result of A^B mod C. (1<=A,B,C<2^63).
Input
There are multiply testcases. Each testcase, there is one line contains three integers A, B and C, separated by a single space.
Output
For each testcase, output an integer, denotes the result of A^B mod C.
Sample Input
私が使っているunsigned_int 64データを開きます.高速べき乗と高速乗算で、注意:%を使わないでください.タイムアウトします(大物が言っています.そして、私は検証しました.233)
コード:
転載先:https://www.cnblogs.com/KirinSB/p/9409142.html
Given A,B,C, You should quickly calculate the result of A^B mod C. (1<=A,B,C<2^63).
Input
There are multiply testcases. Each testcase, there is one line contains three integers A, B and C, separated by a single space.
Output
For each testcase, output an integer, denotes the result of A^B mod C.
Sample Input
3 2 4
2 10 1000
Sample Output 1
24
私が使っているunsigned_int 64データを開きます.高速べき乗と高速乗算で、注意:%を使わないでください.タイムアウトします(大物が言っています.そして、私は検証しました.233)
コード:
#include
#define ull unsigned __int64
ull mul(ull a,ull b,ull c){
ull res=0;
a=a%c;
while(b){
if(b & 1) res+=a;
if(res>=c) res-=c; // %
b/=2;
a+=a;
if(a>=c) a-=c; // %
}
return res;
}
int main(){
ull a,b,c,ans;
while(scanf("%I64u%I64u%I64u",&a,&b,&c)!=EOF){
ans=1;
a=a%c;
while(b){
if(b%2==1) ans=mul(ans,a,c);
a=mul(a,a,c);
b/=2;
}
printf("%I64u
",ans);
}
return 0;
}
転載先:https://www.cnblogs.com/KirinSB/p/9409142.html