乘法逆元应用及求法

乘法逆元的应用以及原理

乘法逆元的定义：若$ax\equiv 1(\mathrm{mod}p)$，就称$x$为$a$模$p$的乘法逆元。

乘法逆元的应用及原理：若$a\mid b$，则$\frac{b}{a}\mathrm{mod}p=b\cdot x\mathrm{mod}p$¹。证明：设$b=ak$，即$k=\frac{b}{a}$，就有$b\cdot x\mathrm{mod}p=ak\cdot x\mathrm{mod}p=k\cdot (ax)\mathrm{mod}p$。因为$ax\equiv 1(\mathrm{mod}p)$，由同余定理得：$k\cdot (ax)\equiv k\equiv \frac{b}{a}(\mathrm{mod}p)$，于是得证。

求乘法逆元

费马小定理²

对于任意的质数$p$和正整数$a$，若满足$\gcd (a,p)=1$，则$a^{p-1}\equiv 1(\mathrm{mod}p)$，那么$a\cdot a^{p-2}\equiv 1(\mathrm{mod}p)$，此时$a^{p-2}$就是$a$模$p$的逆元，快速幂求即可。复杂度$O({\log }_{2}p)$。

扩展欧几里得算法

由于费马小定理对$p$和$a$有限制，所以我们引入另一种方法：求解线性同余方程，也即将原同余方程转化为二元一次不定方程$ax+pk=1$的整数解。复杂度$O({\log }_{2}n)$，其中$n$为$a,p$的取值范围。

线性递推

如果要求出$[1,n]$中每个数模$p$的乘法逆元（不存在则为$0$），我们就可以采用$O(n)$的线性递推方法。首先有$1^{-1}\equiv 1(\mathrm{mod}p)$（$1$就是$1$的乘法逆元）（$x^{-1}$表示$x$的乘法逆元）。对于大于$1$的正整数$x$，我们假设$p=kx+b$，其中$k=\lfloor \frac{p}{x}\rfloor$，$b=p\mathrm{mod}x$。于是可以得到：
$$kx+b\equiv 0(\mathrm{mod}p)$$
两边同时乘以$x^{-1}{b}^{-1}$，得到：
$$k{b}^{-1}+x^{-1}\equiv 0(\mathrm{mod}p)$$
那么可以得到：
$$x^{-1}\equiv -k{b}^{-1}(\mathrm{mod}p)$$
将$k$和$b$的值带入，得到：
$$x^{-1}\equiv -\lfloor \frac{p}{x}\rfloor \cdot (p\mathrm{mod}x{)}^{-1}(\mathrm{mod}p)$$
如果说$x$有一个逆元$r$，那么$r+mp$一定都是$x$的逆元，因为$x(r+mp)\equiv xr\equiv 1(\mathrm{mod}p)$。所以说所有满足$t\equiv -\lfloor \frac{p}{x}\rfloor \cdot (p\mathrm{mod}x{)}^{-1}(\mathrm{mod}p)$的$t$都是$x$的逆元，也就是说(-(p/x)*inv[p%x]%p+p)%p就是$x^{-1}$的最小整数解。由于$0\le p\mathrm{mod}x<x$，所以$x^{-1}$可以由前面的数的逆元求得。如果$p\mathrm{mod}x=0$，那么$x\mid p$，对于$x>1$和任意的$k$，一定有$xk\not\equiv 1(\mathrm{mod}p)$，那么$x$的乘法逆元不存在，因此设为$0$。而我们将$0^{-1}$设为$0$，则求逆元的式子仍然成立。由此式递推即可。

---

1. 对于$a\nmid b$，$\lfloor \frac{b}{a}\rfloor \mathrm{mod}p=\frac{\lfloor \frac{b}{a}\rfloor \cdot a}{a}\mathrm{mod}p$，因此只要将$b$转化为$\lfloor \frac{b}{a}\rfloor \cdot a$即可。 ↩︎

2. 关于费马小定理的证明，可以看看我的这篇文章。 ↩︎