对于给定的a,b,如何在O(log b)时间内求出a^b mod M?
首先我们来看看3^14 mod M,我们可以将其转化为3^2 * 3^4 * 3^8 mod M,然后写一个while循环,先计算3^2 mod M ,随后用该值计算3^4 mod M,……
那么,3^2 * 3^4 * 3^8是怎么得到的呢?
我们可以用二进制作一个中间转换:
1. 将14化为2进制,得1110
2. 由于幂运算关系在二进制下亦成立,则有3^1110=3^1000+3^100+3^10
3. 再化成10进制,即得到3^2 * 3^4 * 3^8
模板:
typedef long long ll;
const ll mod = ...;
ll pow(ll a, ll r) ///a^r % mod
{
ll ans = 1;
while (r)
{
if (r & 1) ans = ans * a % mod;
a = a * a % mod;
r >>= 1;
}
return ans;
}