递推打表求逆元

递推打表求乘法逆元

• 目的：得到在(mod m)的意义下的[2, n)区间内每个数分别对应的逆元
• 推导出递推表达式：
  • 任取一个同余式$ax\equiv c(modm)$消去a与c，得到x关于m的式子即为递推表达式
  • 对于这个同余式，一定存在一个数y使得$ax-c=my$
  • 根据定义$c\ast (c的逆元)\equiv 1(modm)$
  • 我们把c的逆元放在inv数组里，所以$c\ast inv[c]\equiv 1(modm)$
  • 将 c 替换为$ax-my$，$(ax-my)\ast inv[c]\equiv 1(modm)$
  • $ax\ast inv[c]-my\ast inv[c]\equiv 1(modm)$
    • 由于 $my\%m=0$
    • 而原式等价于$(ax\ast inv[c])\%m-(my\ast inv[c])\%m=1$
    • 所以等价于$(ax\ast inv[c])\%m=1$
  • 所以$ax\ast inv[c]\equiv 1(modm)$
    • 在 $ax-c=my中，a=my/x；-c=my\%x$
  • $my/x\ast x\ast inv[-my\%x]\equiv 1(modm)$
  • 对于任意的y都满足此式，不妨取y = -1
  • $-m/x\ast inv[m\%x]\ast x\equiv 1(modm)$
  • 在(mod m)的意义下，$-m/x\ast inv[m\%x]$是x的逆元
  • 为保证求得结果都为正数，$(m-m/x)\ast inv[m\%x]$
  • 所以在(mod m)的意义下，$(m-m/x)\ast inv[m\%x]$是x在区间[0, m)内的逆元

#include<cstdio>

int inv[100001];
int mod = 9973;

//在(mod 9973)的意义下，[2, 10000)分别对应的逆元
void Inv() {
    inv[1] = 1;

    for (int i = 2; i < 10000; i++)
        inv[i] = (mod - mod/i) * inv[mod % i] % mod;
}

int main()
{
    Inv();
    for (int i = 0; i < 20; i++)
        printf("%d ", inv[i]);

    return 0;
}