快速幂
题目描述:
给出三个整数 a,b,ma,b,ma,b,m ,求 ab mod m 的值。
输入格式:
一行三个整数 a、b、m。
输出格式:
一个整数,表示 ab mod m 的值。
输入样例:
2 100 1007
输出样例:
169
提示:
对于全部数据,1≤a,b,m≤109。
思路:
首先需要知道一个原理:
两数之积取模等于对两数分别取模然后对积取模,即 (a * b) % p = (a%p) * (b%p)%p
对于指数较大的情况,为了解决时间复杂度高的问题,引入快速幂。
快速幂的核心思想:利用二进制来加速计算
例如:计算 545
545(10)=5101101(2)=525∗1+24∗0+23∗1+22∗1+21∗0+20∗1=525∗1∗524∗0∗523∗1∗522∗1∗521∗0∗520∗1 5^{45(10)} = 5^{101101(2)} = 5^{2^5*1+2^4*0+2^3*1+2^2*1+2^1*0+2^0*1}=5^{2^5*1}*5^{2^4*0}*5^{2^3*1}*5^{2^2*1}*5^{2^1*0}*5^{2^0*1} 545(10)=5101101(2)=525∗1+24∗0+23∗1+22∗1+21∗0+20∗1=525∗1∗524∗0∗523∗1∗522∗1∗521∗0∗520∗1因为:524∗0=521∗0=15^{2^4*0}=5^{2^1*0}=1524∗0=521∗0=1
根据指数的运算性质:n2m−1∗n2m−1=n2mn^{2^{m-1}}*n^{2^{m-1}}=n^{2m}n2m−1∗n2m−1=n2m 可以得到 525=524∗524、524=523∗523、523、522、521、5205^{2^5}=5^{2^4}*5^{2^4}、5^{2^4}=5^{2^3}*5^{2^3}、5^{2^3}、5^{2^2}、5^{2^1}、5^{2^0}525=524∗524、524=523∗523、523、522、521、520
综上:
如果一个指数的二进制形式最低位为 1 ,则累乘项数,为 0 ,则跳过。
再计算下一项
将指数二进制右移一位(使用位运算符)
代码:
#include<stdio.h>
int main()
{
int b,m;
long long a=0,s=1;
scanf("%lld %d %d",&a,&b,&m);
while(b)
{
if(b&1) //判断指数的二进制最低位是否为 0
s=s*a%m; //如果为 0 ,则和前面的项累乘并存储
a=a*a%m; //计算下一项(例如5^4=5^2*5^2)
b>>=1; //二进制右移一位
}
printf("%lld",s);
return 0;
}