【ACM数论】快速幂

快速幂

通常地，我们计算$a^x$是将$x$从1一直累乘到$x$，这样是最简单的计算方法，但是这种做法的时间复杂度是$o(n)$，一旦我们计算的幂次方较大或者计算次数较多的时候就会花费很长的时间；例如：计算$10^{100000000}$组类似$999^{100000000}$这样的计算式，你可能算一辈子都算不出来。

​ 这个时候引入二进制对$a^x$进行优化：

​ 考虑$a^{2x}$=$a^x\times a^x$=$(a^x)$${}^2$于是可以利用二进制表示分割成更小的一部分

​ 例如：$9^{13}$=$9^{(1101{)}_2}$=$9^8$×$9^4$×$9^1$

​ 因为$n$有$[lo{g}_{2}x]+1$个二进制位，因此当我们知道了 $a^1,a^2,a^4...a^{2^k}$后，我们只用计算$[lo{g}_{2}x]+1$次乘法就可以计算出$a^x$

​ 因此计算$a^x$只需要将$x$对应的二进制下的数位为1的整系数幂乘起来即可

​ 即：$a^x$=$a^{(x_0{x}_1{x}_2..x_n{)}_2}$=$a^{x_0{2}^0}$×$a^{x_1{2}^1}$×$a^{x_2{2}^2}$…×$a^{x_n{2}^n}$

​ 根据上式发现，原问题被转化成了形式相同的子问题的乘积，并且我们可以在常数时间内从 $2^i$ 项推出$2^{i+1}$项。

​ 通过二进制处理后时间复杂度降为$O(logn)$

​ 下面是我的快速幂代码模版：

long long fastpow(long long a, long long n) //快速幂
{
    long long ans = 1;
    while (n)
    {
        if (n & 1)
        {
            ans *= a;
        }
        a *= a;
        n >>= 1; //位运算
    }
    return ans;
}