快速幂

思路：

​ 1.每一步都把指数分成两半，而相应的底数做平方运算。这样不仅能把非常大的指数给不断变小，所需要执行的循环次数也变小，而最后表示的结果却一直不会变。

​ 例：$3^{10}=3\times 3\times 3\times 3\times 3\times 3\times 3\times 3\times 3\times 3$ ，尽量想办法把指数变小来，这里的指数为10。

​ $3^{10}=(3\times 3)\times (3\times 3)\times (3\times 3)\times (3\times 3)\times (3\times 3)$

​ $3^{10}=(3\times 3{)}^5$

​ $3^{1}0=9^5$

​ 2.此时指数由 $10$ 缩减一半变成了 $5$ ，而底数变成了原来的平方，求 $3^{10}$ 原本需要执行 $10$ 次循环操作，求 $9^5$ 却只需要执行 $5$ 次循环操作，但是 $3^{10}$ 却等于 $9^5$，用一次（底数做平方操作）的操作减少了原本一半的循环量，特别是在幂特别大的时候效果非常好，例如 $2^{10000}=4^{5000}$，底数只是做了一个小小的平方操作，而指数就从 $10000$ 变成了$5000$ ，减少了$5000$ 次的循环操作。

​ 3.现在问题是如何把指数 $5$ 变成原来的一半，$5$ 是一个奇数，$5$ 的一半是 $2.5$ ，但是指数不能为小数，因此不能简单粗暴地直接执行 $5/2$，然而，这里还有另一种方法能表示 $9^5$，$9^5=9^4\times 9^1$

​ 4.此时抽出了一个底数的一次方，这里即为 $9^1$，这个 $9^1$ 先单独移出来,剩下的 $9^4$ 又能够在执行“缩指数”操作了，把指数缩小一半，底数执行平方操作。$9^5=81^2\times 9^1$

​ 5.把指数缩小一半，底数执行平方操作，$9^5=6561^1\times 9^1$

​ 6.此时，发现指数又变成了一个奇数1，按照上面对指数为奇数的操作方法，应该抽出了一个底数的一次方，这里即为$6561^1$，这个 $6561^1$ 先单独移出来，但是此时指数却变成了 $0$ ，也就意味着我们无法再进行“缩指数”操作了。

​ $9^5=(6561^0)\times (9^1)\times (6561^1)=1\times (9^1)\times (6561^1)=(9^1)\times (6561^1)=9\times 6561=59049$

​ 7.能够发现，最后的结果是 $9\times 6561$。所以能发现一个规律：最后求出的幂结果实际上就是在变化过程中所有当指数为奇数时底数的乘积。

​ 8.继续优化：$b$ % $2$==$1$可以用更快的“位运算”来代替，例如：$b$ & $1$。因为如果 $b$ 为偶数，则其二进制表示的最后一位一定是 $0$ ；如果 $b$ 是奇数，则其二进制表示的最后一位一定是 $1$ 。将他们分别与 $1$ 的二进制做“与”运算，得到的就是 $b$ 二进制最后一位的数字了，是 $0$ 则为偶数，是 $1$ 则为奇数。例如 $9$ 是奇数，则 $9$ &$1=1$；而 $8$是偶数，则 $8$ &$1=0$；因此奇偶数的判断就可以用“位运算”来替换了。

​ 9.$m=m/2$ 也可以用更快的移位操作来代替，例如：$6$ 的四位二进制为 $0110$ ，而 $6/2=3$，$3$ 的四位二进制为 $0011$，可以发现，$a$ 的一半，结果为 $a$ 的二进制码向右移一位，即 $m>>=1$。

public static long num(long n, long m, long p) {
    long result = 1;
    while (m > 0) {
        if ((m & 1 ) == 1) {
            result = result * n % p;
        }
        m >>= 1;
        n = (n * n) % p;
    }
    return result;
}