算法 - 数论 - powmod 快速幂取模

理论

• 定理1：$(a+b)\mathrm{m}\mathrm{o}\mathrm{d}c=(a\mathrm{m}\mathrm{o}\mathrm{d}c+b\mathrm{m}\mathrm{o}\mathrm{d}c)\mathrm{m}\mathrm{o}\mathrm{d}c$
  证明：
  设$a=xc+x_0,b=yc+y_0$
  则$(a+b)÷c=(xc+x_0+yc+y_0)÷c$
  其中$xc÷c=x...0$
  $yc÷c=y...0$
  显然有 $(a+b)\mathrm{m}\mathrm{o}\mathrm{d}c=(xc+x_0+yc+y_0)\mathrm{m}\mathrm{o}\mathrm{d}c=(x_0+y_0)\mathrm{m}\mathrm{o}\mathrm{d}c$
  即：$(a+b)\mathrm{m}\mathrm{o}\mathrm{d}c=(a\mathrm{m}\mathrm{o}\mathrm{d}c+b\mathrm{m}\mathrm{o}\mathrm{d}c)\mathrm{m}\mathrm{o}\mathrm{d}c$

• 定理2：$(a\times b)\mathrm{m}\mathrm{o}\mathrm{d}c=((a\mathrm{m}\mathrm{o}\mathrm{d}c)\times (b\mathrm{m}\mathrm{o}\mathrm{d}c))\mathrm{m}\mathrm{o}\mathrm{d}c$
  证明:
  设$a=xc+x_0,b=yc+y_0$
  则
  $$(a\times b)\mathrm{m}\mathrm{o}\mathrm{d}c=((xc+x_0)\times (yc+y_0))\mathrm{m}\mathrm{o}\mathrm{d}c=(xy{c}^2+x_0{y}_0+x_{0}yc+y_{0}xc)\mathrm{m}\mathrm{o}\mathrm{d}c=(x_0{y}_0)\mathrm{m}\mathrm{o}\mathrm{d}c=((a\mathrm{m}\mathrm{o}\mathrm{d}c)\times (b\mathrm{m}\mathrm{o}\mathrm{d}c))\mathrm{m}\mathrm{o}\mathrm{d}c$$

• 推理1: $x^{N_0}\mathrm{m}\mathrm{o}\mathrm{d}c=((mx\mathrm{m}\mathrm{o}\mathrm{d}c)\times (x^{N_1}\mathrm{m}\mathrm{o}\mathrm{d}c{)}^2)\mathrm{m}\mathrm{o}\mathrm{d}c;N_0=m+2\times N_1$

实现

• 递归

#include <cinttypes>

uintmax_t pow_mod(const uintmax_t x,const uintmax_t n, const uintmax_t c = INT_MAX)
{
    if (n == 0) return 1;
    if (n == 1) return x%c;
    if (n == 2) return ((x%c) * (x%c)) % c;
    const auto y = pow_mod(x, n >> 1,c) % c;
    if (n % 2) return ((x % c) * ((y * y) % c)) % c;
    return (y * y) % c;
}

• 循环

#include <cinttypes>

uintmax_t pow_mod(const uintmax_t x,const uintmax_t n, const uintmax_t c = INT_MAX)
{
    uintmax_t res = 1;
    uintmax_t y = x;
    for (uintmax_t i = n; i != 0;i >>= 1) {
        if (i & 1) res = (res * y) % c;
        y = (y * y) % c;
    }
    return res;
}