快速幂算法

$S$表示$p$的二进制权位序列：$2^n，...，2^3，2^2，2^1，2^0$。
$B$表示$p$的二进制序列。
$b^p$可以表示成
$$b^p=\prod _{k=0}^n{b}^{S_k\ast B_k}$$
记 $T_k=b^{S_k\ast B_k}$，
$$T_k=\{\begin{array}{ll}1& B_k=0\\ b^{2^k}& B_k=1\end{array}$$

可以观察发现，
（1）$B_k$决定了第$k$项 $T_k$ 是否应该放入累乘中。
（2）存在递推关系$T_{k+1}=T_k\ast T_k$。

以$3^{13}$为例。13的二进制序列为1101，从右开始数，0号位为1，1号位为0，2号位为1，3号位为1。则
$$B=1,1,0,1$$$$S=2^3,2^2,2^1,2^0$$$$T=3^8,3^4,3^2,3^1$$,故$3^{13}=3\ast 1\ast 3^4\ast 3^8$。

1，快速幂的迭代法：

long long binaryPow(long long b, long long p)
{
    long long r = 1;
    while(p>0)
    {
        if(p & 1)   // 判断p的二进制末位是不是1
        {
            r = r * b;
        }
        p = p >> 1;    // p的二进制右移一位
        b = b * b;
    }
    return r;
}

2，快速幂取余：
涉及到同余的性质$$abmodn=[(amodn)(bmodn)]modn$$即乘积的模等于乘数模的乘积的模。

long long binaryPowMod(long long b, long long p, long long m)
{
    long long r = 1;
    while(p>0)
    {
        if(p & 1)   // 判断p的二进制末位是不是1
        {
            r = r * b % m;
        }
        p = p >> 1;    // p的二进制右移一位
        b = b * b % m;
    }
    return r;
}

推荐一篇通俗易懂的博客可以参考：快速幂算法从零开始一步一步优化

3.例题

剑指 Offer 16. 数值的整数次方