快速幂
功能
快速计算 a b a^b ab
思路
若将 b b b表示为 ∑ p i × 2 i \sum p_i \times 2^i ∑pi×2i,则 a b a^b ab可以表示为 ∏ ( a 2 i ) p i \prod (a^{2^i})^{p_i} ∏(a2i)pi,其中 p i p_i pi表示 b b b的二进制从右往左第 i i i位数字。
时间复杂度
O ( log n ) O(\log n) O(logn)
测试
HDU:1061
模板
typedef long long LL;
/**
* @param a: the base of power
* @param b: the exponent of power
* @return: a^b
* @other: b >= 0
*/
LL FP(LL a, LL b) {
LL ans = 1;
while (b > 0) {
if (b & 1 == 1) ans *= a; // pi
a = a*a; // a^(2^i)
b >>= 1;
}
return ans;
}
扩展
对于较大的数需要取模,由 ( a m o d    c ) × ( b m o d    c ) ≡ ( a × b ) m o d    c (a \mod c) \times (b \mod c) \equiv (a \times b) \mod c (amodc)×(bmodc)≡(a×b)modc得 ∏ ( a 2 i ) p i m o d    c ≡ ( ∏ ( ( a m o d    c ) 2 i ) p i ) m o d    c \prod (a^{2^i})^{p_i} \mod c \equiv (\prod ((a \mod c)^{2^i})^{p_i}) \mod c ∏(a2i)pimodc≡(∏((amodc)2i)pi)modc
模板
typedef long long LL;
const LL MOD = 1e9+7; // the divisor of answer
/**
* @param a: the base of power
* @param b: the exponent of power
* @return: a^b
* @other: b >= 0
*/
LL FP(LL a, LL b) {
a %= MOD;
LL ans = 1;
while (b > 0) {
if (b & 1 == 1) ans = (ans * a) % MOD; // pi
a = (a * a) % MOD; // (a mod c)^(2^i)
b >>= 1;
}
return ans;
}