qpow快速幂乘法
思路:
形如a^b,当b很大时可以对指数进行拆分。
若把b写成二进制,则b有 [log_2b]+1 个二进制位
b=b_0⋅20+b_1⋅21+…+b_n⋅2^n,其中b_i=0或1,n= [log_2b]
则ab=ab_0∗2^0* ab_1∗21… ab_n∗2n,
总共只需要乘(n+1)次,而(n+1)远远小于b。
思路:构造循环,对于b的二进制形式,每次都检验最低位,若为1则获得该位的权值,获得乘数
ab_n∗2n;后对b右移运算,把本次检验位舍弃,上一位移至最低位再度检验…直至b的值变为零
(所有位均被舍弃),终止循环。
代码如下
long long qpow(long long i,long long k,long long mod){//加上了取mod的运算
long long sum_i_k = 1;
while(k){
if(b&1){
sum_i_k=sum_i_k*i%mod;
}
i=i*i%mod;//真正的运算简便在于i的变大
b>>1;//b只是开关
}
}

浙公网安备 33010602011771号