快速幂
typedef long long ll;
ll quick(ll a,ll b){
ll s=1;
while(b > 0){
if( b&1 ) s *= a; //如果b不是2的倍数,即b二进制此位为1,s则*=a
b /= 2; //b往下一位二进制数走
a *= a; //b走往下一位,a同理*=a;
}
return s;
}
//快速幂取模算法
//取模的一些公式
(a*b) % m = (a % m) * (b % m) % m
(a/b) % M = (a * b^(M-2)) % M; //费马小定理, M 为素数
(a+b) % M = (a % M+b % M) % M;
快速乘
举个栗子,你就知道快速乘的作用。。。假如现在两个数 a 和 b 都是 1e18 范围内的数,求 (a * b) % 9999999999999999(16个9),你会发现直接算 a * b 会存在爆 long long 的情况,就算你先取模也没用
#define LL long long
LL multiply(LL a, LL b, LL mod){
LL ans = 0;
while (b){
if (b & 1){
ans = (ans + a) % mod;
}
a = (a + a) % mod;
b >>= 1;
}
return ans;
}