快速幂
题目描述:
给出三个整数 a , b , m a,b,m a,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
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 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因为: 5 2 4 ∗ 0 = 5 2 1 ∗ 0 = 1 5^{2^4*0}=5^{2^1*0}=1 524∗0=521∗0=1
根据指数的运算性质: n 2 m − 1 ∗ n 2 m − 1 = n 2 m n^{2^{m-1}}*n^{2^{m-1}}=n^{2m} n2m−1∗n2m−1=n2m 可以得到 5 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 5^{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;
}