您的位置:首页 > 其它

(A^B)%C 快速幂乘

2015-10-24 16:52 225 查看
<span style="font-size:18px;"><strong>简单数论的快速幂乘问题</strong></span>
#include<iostream>
using namespace std;
int main()
{
int a,b,c;
while (scanf("%d%d%d",&a,&b,&c)!=EOF&&a!=0&&b!=0&&c!=0)
{
int ans = 1;
a=a%c;
while (b>0)
{
if(b%2==1)     // 作用1:当 b为奇数,则先单独乘一个a
// 作用:当 b=1时,即 已经乘了 b=b/2=1后,将值赋给k
ans = (ans * a) % c;
a = (a * a) % c;
b = b / 2;       //a1=a%c                             1=2^0
// a2=((a%c)*(a%c))%c                2=2^1
// a3=((a%c*a%c)%c*(a%c*a%c)%c)      4=2^2
//由此,可知 b=b/2 ,每次的 a 的个数为上一次的 2倍
}
printf("%d\n",ans);
}
return 0;
}


引用快速幂乘的经典实例uva 11582

思路:设g(i)=f(i) mod n,当二元组(g(i)、g(i+1))出现重复时,整个序列就开始重复

余数最多n种,所以最多n^2项就会出现重复。设周期为m,只需计算出g(0)~g(n^2)项,然后计算g(a^b)等于其中的哪一项即可

但是,n<=1000,那么n方的规模会达到10^6,显然f(n^2)有些庞大

我们可以利用(a+b)mod n=((a mod n)+(b mod n))mod n 优化

那么g(i+2)=g(i+1)+g(i)=(f(i+1)mod n + f(i)mod n )mod n=(g(i+1)+g(i))mod n



#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 2005001;
int g[maxn];
unsigned long long a,b;
int PowerMod(unsigned long long a,unsigned long long b,int c)
{
int ans = 1;
a = a % c;
while(b>0)
{
if(b % 2 == 1)
ans = (ans * a) % c;
a = (a * a) % c;
b = b / 2;
}
return ans;
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
int n;
cin >> a >> b >> n;
if(n == 1)
{
puts("0");
continue;
}
int m = 1;
g[1] = 1,g[2] = 1;
int i = 3;
while(1)
{
g[i] = (g[i - 2] + g[i - 1]) % n;
if(g[1] == g[i-1] && g[2] == g[i])
break;
i++;
}
m = i - 2;
int s = PowerMod(a,b,m);
cout<<g[s]<<endl;
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: