问题:
给定三个整数a, b, m(a<1000000000, b<1000000, 1<m<1000000000),求(aEb)%m。
求余运算 (ab)%m=(a%m)(b%m);
快速幂的二分思想:
*如果b是奇数, 那么aEb=a * aE(b-1);
如果b是偶数, 那么aEb=aE(b/2) aE(b/2);
b是奇数的情况总是可以在下一步转换为b是偶数的情况,b是偶数总可以化解为b/2,经过logb之后就可以将b变成零。
typedef long long LL;
LL binaryPow(LL a, LL b, LL m)
{
if(b==0)
{
return 1;
}
if(b%2==1)
{
return (a*binaryPow(a, b-1, m))%m;
}
else
{
LL mul=binaryPow(a, b/2, m);
return (mul*mul)%m;
}
}
快速幂的迭代写法:
对于a的13次方, 13的二进制表示就是1101, 即1+4+8, 按照这个规律我们可以把任意b都写成二进制形式然后表示成a, a^2, a^4, a^8 a^2k的乘积的形式
代码实现如下:
1.初始令ans等于一, 用来存放累乘的结果
2.判断b的二进制末尾是否为一, 如果是的话, 令b乘上a
3.令a平方, 将b右移一位(即向前找一个二进制数字,1101 110 11 1)
4.b如果b大于零,返回2
typedef long long LL;
LL binaryPow(LL a, LL b , LL m)
{
LL ans=1;
while(b>0)
{
if(b&1)//也可以写为if(b%2)
{
ans=(ans*a)%m;
}
a=(a*a)%m;
b>>=1;
}
return ans;
}