C++快速幂函数

该快速幂函数用于求 $a^b\%p$ 的结果，既用来计算$a$的$b$次幂取模$p$的结果。
函数原理：
(利用费马小定理)¹
依次判断指数$b$的二进制位有哪些是$1$。如果二进制最后一位是$1$,就将最终答案$ans$乘上$a$，结束判断后再将$a$平方，同时将$b$右移$1$位，直到$b$等于$0$为止。
注意：在每一次计算之后，都要取模 $p$。

---

函数代码：

long long fast_power(long long a, long long b, long long p)
{ // 快速幂函数，求解a^b%p
    long long ans = 1;
    while (b) // 直到b等于0停止循环
    {
        if (b & 1) // 如果二进制最后一位是1
            ans = ans * a % p; // 将最终答案ans乘上a
        a = a * a % p; // 将a平方
        b >>= 1; // 将b右移1位
    }
    return ans;
}

---

1. 费马小定理(Fermat’s little theorem)是数论中的一个重要定理，在$1636$年提出。如果$p$是一个质数，而整数$a$不是$p$的倍数，则有$a^{(p-1)}\equiv 1(modp)$。
   跳转百度百科 ↩︎