计算逆元的三种方法

$$计算逆元的三种方法$$

一、首先是两种O(log)时间求单个数的逆元

我们要求$a$在模$p$意义下的逆元
其中$a$是小于$p$的正整数，$p$是一个素数

① 费马小定理：

由费马小定理$a^{p-1}\equiv 1(modp)$可以得到：
$a\cdot a^{p-2}\equiv 1(modp)$
所以可以直接使用快速幂求出$a^{p-2}$即为$a$在模$p$意义下的逆元

② 扩展欧几里得：

先假设我们要求的$a$在模$p$意义下的逆元为$x$
可以得到方程$a\cdot x\equiv 1(modp)$
变成一般等式：$a\cdot x+p\cdot y=1$
现在要求的就是$x$
$p$是素数，所以$gcd(a,p)=1$，可以使用扩展欧几里得计算得到$x$。

二、线性求逆元：

显然$1^{-1}=1$
假设现在已经计算出$i$之前的所有正整数的逆元
现在要求$i$的逆元
先写出等式：$i\cdot k+r=p$
变成模$p$下的恒等式：$i\cdot k+r\equiv 0(modp)$
两边同乘$i^{-1}\cdot r^{-1}$得到：$k\cdot j^{-1}+i^{-1}\equiv 0(modp)$
移相：$i^{-1}\equiv -k\cdot j^{-1}(modp)$
$k$和$j$可以分别用$p/i$和$p\%i$来代替：$i^{-1}\equiv -(p/i)\cdot (p\%i{)}^{-1}(modp)$
$$取正数：i^{-1}=(p-p/i)\cdot (p\%i{)}^{-1}\%p$$
$p\%i<i$ 所以$p\%i$的逆元已经算出来了，右端项没有未知数。所以就可以$O(n)$的时间内算出$1-n$之内所有数的逆元了

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