【数学】乘法逆元进阶

乘法逆元基础

可以使用扩展欧几里得算法或费马小定理求得。

乘法逆元进阶

有一种乘法逆元的线性递推算法：
显然 $1^{-1}\equiv 1(\mathrm{mod}m)$；
对于 $i^{-1}$，我们令
$k=\lfloor \frac{m}{i}\rfloor$，$j=m\mathrm{mod}i$，有 $m=ki+j$，即 $ki+j\equiv 0(\mathrm{mod}m)$；
两边同时乘 $i^{-1}\times j^{-1}$ 得：$k{j}^{-1}+i^{-1}\equiv 0(\mathrm{mod}m)$
即 $i^{-1}\equiv -k{j}^{-1}(\mathrm{mod}m)$
代入 $m=ki+j$ 得：

$$i^{-1}\equiv -\lfloor \frac{m}{i}\rfloor (m\mathrm{mod}i{)}^{-1}(\mathrm{mod}m)$$
故利用迭代易知。

还有一个求 $i{!}^{-1}$ 时的小技巧如下：
假设我们要求 $i\in [1,N]$ 时 $i!$ 在模 $m$ 意义下的逆元。
不妨先使用费马小定理求出 $N{!}^{-1}$
接着我们发现 $i{!}^{-1}\cdot i\equiv (i-1{)}^{-1}$
故可以简单地线性递推。

代码

• $i^{-1}$

inv[1]=1;
for (int i=2;i<=n;i++) inv[i]=(long long)(M-M/i)*inv[M%i]%M;

• $i{!}^{-1}$

fac[0]=1;
for (int i=1;i<=n;i++) fac[i]=fac[i-1]*i%M;
inv[n]=pow(fac[n],M-2);
for (int i=K-1;i>=0;i--) inv[i]=inv[i+1]*(i+1)%M;