2 3 12 6 6789 10000 0 0
8 984 1第一发 数论,问题描述:求 a^n%m最简单的思想是用循环分部取余,在n很小的情况下是可以过的,如下#include <iostream> #include <cstring> #include <cstdio> #include <cctype> #include <cstdlib> #include <algorithm> #include <set> #include <vector> #include <string> #include <map> #include <queue> using namespace std; int main() { int a,b; while(cin>>a>>b) { if(!a&&!b)break; int p=1; for(int i=1;i<=b;i++) p=p*a%1000; cout<<p<<endl; } return 0; }
但是当n很大的时候这个循环就要挂掉了,所以要用到数学方法。对指数n进行二分,就会有a^n%m= a^(n/2)%m n为偶数a^(n/2)%m * a%m n为奇数非递归 代码如下:#include <iostream> #include <cstring> #include <cstdio> #include <cctype> #include <cstdlib> #include <algorithm> #include <set> #include <vector> #include <string> #include <map> #include <queue> using namespace std; #define LL long long LL multimod(LL a,LL n,LL m) { LL ans=1,tem=a; while(n) { if(n&1)//判断是否为奇数 ans=ans*tem%m; tem=tem*tem%m; n/=2; } return ans; } int main() { LL a,b; while(cin>>a>>b) { if(!a&&!b)break; cout<<multimod(a,b,1000)<<endl; } return 0; }
原文地址:http://blog.csdn.net/qq_16255321/article/details/39022201