原理:
我们用分治思想是比一个一个乘快的
即比如我们求a的8次方 :a1a1 = a2 ,那么我们直接a2a2 = a4,a4*a4 = a8
参数就是几次幂。返回值是对应参数的幂 (这里对p取余了)(一般也把a当参数)
long long a, p;//a^c mod p=s
long long fp(int c)
{
if(c==0)return 1;
if (c == 1)return a % p;
long long tmp = fp(c / 2);
if (c % 2 == 0)
return tmp * tmp % p;
else
return tmp * tmp % p * a % p;
}
long long fp(long long a,int c)
{
if(c==0)return 1;
if (c == 1)return a % MOD;
long long tmp = fp(a,c / 2);
if (c % 2 == 0)
return tmp * tmp % MOD;
else
return tmp * tmp % MOD * a % MOD;
}
细节解释:
c == 1就是1次幂。
tmp就是a的c/2次幂。我们要返回c次幂,整数除法是向下取整的。
(比如5次幂,5/2==2,那么需要额外乘一个使得为c次幂)
——————————
非递归写法:(二进制拆分)
#define ll long long
ll mod = 1e9+7;
ll q(ll a,ll b)
{
ll s = 1;
while(b)
{
if(b&1) s = s * a %mod;
a = a * a % mod;
b >> 1;
}
return s;
}
———————
快速幂求逆元:
前提:
费马小定理:
1
≡
a
p
−
1
(
m
o
d
p
)
1 ≡ a^{p-1} ~(mod ~p)
1≡ap−1 (mod p) ( 1和 a^(p-1)在mod质数p下同余 )
逆元也就是倒数,若我们求a的逆元x,则满足:
a x ≡ 1 ( m o d p ) ax ≡ 1~ (mod ~p) ax≡1 (mod p)
而由费马小定理可得
a
x
≡
a
p
−
1
(
m
o
d
p
)
ax ≡ a^{p-1}~ (mod ~p)
ax≡ap−1 (mod p)
即
x
≡
a
p
−
2
(
m
o
d
p
)
x ≡ a^{p-2}~ (mod ~p)
x≡ap−2 (mod p)