乘法逆元的一些求法

逆元

• 【逆元素-百度一下】

广义的来讲，对于任何域中的元素，有乘法运算和单位元$I$，如果对于该域中的元素$A$，存在另一个元素$A'$，且满足$A imes A'=I$，那么$A'$就是$A$的逆元。

这里我们只讨论在整数域里的逆元，也就是当$Ain m Z$且$I=1$，其实这里的逆元$A'=frac{1}{A}$，但是我们要在模的意义下讨论它的求法。

---

在取模意义下，我们只需求出一个数$inv$，是的这个数与$frac{1}{A}$同余即可，那么这个数$inv$就是$A$在取模意义下的逆元。

费马小定理

内容：对于$x$，在模$p,pin mathbb{P} ext{质数}$的意义下，有$x^{p-1}equiv 1(mod p)$

我们变换一个形式，左右同时除以$x$，就是$x^{p-2}equiv frac{1}{x}(mod p)$，那么我们就发现$x^{p-2}$和$frac{1}{x}$是同余的，那么$x^{p-2}$就是$x$的逆元了。

费马小定理证明：

1. 用欧拉定理直接证明，这个就后面再说。
2. 我们用剩余系来看：

引理1：假设$a$为$1sim p-1$内的一个数，那么$a,2a,3a,cdots,(p-1)a$可以在模$p$的意义下不重复的取遍$1sim p-1$的值。引理1的证明-这里面有

那么我们将$a,2a,cdots,(p-1)a$乘起来，在模$p$的意义下也就相当于把$1,2,cdots,(p-1)$乘起来，那么可以得到$a^{p-1}(p-1)!equiv (p-1)!(mod p)$，然后两端同时除以$(p-1)!$就可以得到$a^{p-1}equiv 1(mod p)$了，于是得证。

所以用费马小定理求取逆元的代码如下：

int inv(int a,int b=p-2){
	int ans=1;
	for(;b;b>>=1,a=(a*a)%p)if(b&1)ans=(ans*a)%p;
	return ans%p;
}

---

欧拉定理

• 百度百科

内容：对于$x$，在模$p$（$p$为任意与$x$互质的数）意义下，有$x^{varphi(p)}equiv 1(mod p)$，其中$varphi(p)$为欧拉函数，表示$1sim p-1$中与$p$互质的个数。

那么我们知道，当$p$为质数时就有$varphi(p)=p-1$，那么这个就是费马小定理了。

我们同样进行变形，可以得到$x^{varphi(p)-1}equiv frac{1}{x}(mod p)$，那么这里$x^{varphi(p)-1}$就为$x$的逆元了。

我们看，这里只要求$x,p$互质，并没有要求$p$为质数，所以欧拉定理求逆元比费马小定理求逆元应用更加广泛。

其他应用：当快速幂取模时，指数较大，但我们发现$x^{varphi(p)}equiv 1(mod p)$，所以我们可以将指数模了$varphi(p)$，相当于除以很多个$1$，然后再来快速幂。

欧拉定理的证明：

我们先令$x_1sim x_{varphi(p)}$为$p$的简化剩余系，再令$p_i=xx_i$，其中$gcd(x_i,p)=1$

引理1. $p_i$之间两两在模$p$的意义下不同余，同样$x_i$也是。

• 证明引理1. ：

假设有$p_iequiv p_j(mod p)$，那么$p_i-p_jequiv 0(mod p)$，我们将其变形得$xx_i-xx_jequiv 0(mod p)$，也就是$x(x_i-x_j)equiv 0(mod p)$，由于$gcd(x,p)=1$，所以只有在$x_i-x_j$为$p$的倍数时，这个式子才满足，而$x_i,x_jleq p$，那么它们相减就更加不可能为$p$的倍数了，所以假设不成立，得证。

引理2. 每个$p_i$模$p$的结果都与$p$互质。

• 证明引理2. ：

假设$p_iequiv r(mod p),gcd(p,r)>1$
那么变形就有$xx_i=kp+r$，也就是$xx_i-kp=r$，我们可以得知$xx_i$肯定与$p$互质（因为$x,x_i$分别与$p$互质），那么我们令$d=gcd(p,r)$，那么变形可知$xx_i=r+kp$，也就是$xx_i=d(w_1+kw_2) ext{其中}(dw_1=r,dw_2=p)$，那么$xx_i$肯定有因子$d$，那么它与$p$就不互质了，所以假设不成立，得证。

• 证明欧拉定理：

由1,2引理可得，$(xx_1 imes xx_2cdots imes xx_{varphi(p)})equiv (x_1 imes x_2cdots imes x_{varphi(p)})(mod p)$

所以变形得$x^{varphi(p)}prod x_iequiv prod x_i(mod p)$，两端同除就得$x^{varphi(p)}equiv 1(mod p)$

扩展欧拉定理

内容：
$x ^ { k } = left{ egin{array} { l l } { x ^ { k mod varphi ( p ) } } & { g c d ( x , p ) = 1 } \ { x ^ { k } } & { k < varphi ( p ) } \ { x ^ { k mod varphi ( p ) + varphi ( p ) } } & { g c d ( x , p ) eq 1 , k geq varphi ( p ) } end{array} ight.$

---

扩展欧几里得算法

内容 : exgcd(a,b,x,y) 【见详解于这里】

这个也就是求方程$ax+by=c$的解。
那么我们将$axequiv 1(mod y)$看作方程$ax+by=1$，那么当$gcd(a,b)=1$的时候，这个肯定有解，所以我们用扩展欧几里得定理求出$x$，那么可以发现$axequiv 1(mod y)$的，所以此时的$x$就是在$mod y$意义下的逆元。

---

线性递推求法

首先我们知道$1$的逆元一定为$1$，所以考虑$i$的逆元($i>1$)。

对于模数为$p$，我们令$p=i imes k+d(d<i)$，那么对于$k,d$则有：$k=lfloorfrac{p}{i} floor,d=p mod i$

那么原式可以写成$i imes k+dequiv 0(mod p)$

两边同时乘以$i^{-1} imes d^{-1}$可得

$d^{-1} imes k+i^{i-1}equiv 0(mod p)$

那么$i^{-1}equiv (-k imes d^{-1})(mod p)$，所以将$k,d$带回即有：

$i^{-1}equiv -lfloorfrac{p}{i} floor imes(p mod i)^{-1}(mod p)$

由于$i$以内的逆元已经求出，所以递推即可。

所以$inv(i)=(-frac{p}{i}) imes inv(p mod i)$

---

参考大佬文章IN-Luogu