快速幂
01
02
03
04
05
06
07
08
09
10
11
12
|
int
pow
(
int
a ,
int
k)
{
int
rec = 1;
while
( k )
{
if
(k & 1)
rec *= a;
a *= a;
k >>= 1;
}
return
rec;
}
|
举个例子,3^13;13=1*2^3+1*2^2+0*2^1+1*2^0;那么第一次K=13;这个时候K&1表示探查2^0的这一位是否为1;为1,则进行rec*=a; k>>=1表示将13右移一位,即13/2;可以在之前的13的上面直接表示出来.6=1*2^2+1*2^1+0*2^0;这个时候,k&1则为0了。因为2^0的系数是0;
快速幂取模
01
02
03
04
05
06
07
08
09
10
11
12
|
int
ModPow(
int
a,
int
b,
int
n)
{
int
rec=1;
while
(b)
{
if
(b & 1)
rec = (rec * a) % n;
a = (a * a) % n;
b >>= 1;
}
return
rec % n;
}
|
模运算规则:
模运算与基本四则运算有些相似,但是除法例外。其规则如下:
(a + b) % p = (a % p + b % p) % p
(a – b) % p = (a % p – b % p) % p
(a * b) % p = (a % p * b % p) % p
ab % p = ((a % p)b) % p
结合率:
((a+b) % p + c) % p = (a + (b+c) % p) % p
((a*b) % p * c)% p = (a * (b*c) % p) % p