乘法逆元【解法大全】

定义

若$a\ast x\equiv 1(modp)$，且 $a$ 与 $p$ 互质，那么我们就能定义: $x$ 为 $a$ 的逆元，记为 $a^{-1}$，所以我们也可以称 $x$ 为 $a$ 的倒数( $mod$ $p$ 意义下)。

所以对于 $\frac{a}{b}modp$ ，我们就可以求出 b 在 $modp$意义下的逆元，然后乘上 $a$ ，再 $modp$ ，就是这个乘法逆元的值了。

---

模板题

原题链接》》乘法逆元

---

法一

两种证明：

（1）费马小定理：

定理内容：如果 $a,p$ 互质，那么 $a^{p-1}\equiv 1(modp)$

把上式转换成 $a\ast a^{p-2}\equiv 1(modp)$，结合逆元方程 $a\ast x\equiv 1(modp)$，得到 $a\ast a^{p-2}\equiv a\ast x(modp)$

所以 $a^{p-2}\equiv x(modp)$，即 $x=a^{p-2}modp$, 用快速幂求解即可

（2）欧拉定理
定理内容：如果a,p互质，那么$a^{\phi (p)}\equiv 1(modp)$，当p为质数时，$\phi (p)=p-1$。
所以$a^{p-1}\equiv 1(modp)$，跟费马小定理的式子长得一模一样，按刚才的方式推导即可得 $x=a^{p-2}modp$

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define ll long long

using namespace std;

ll n,p;

ll ksm(ll a,ll b)
{
	ll ans=1,base=a;
	while(b)
	{
		if(b&1)
			ans=(ans*base)%p;
		base=(base*base)%p;
		b>>=1;
	}
	return ans;
}

int main(){
	scanf("%lld%lld",&n,&p);
	for(ll i=1;i<=n;i++)
	{
		printf("%lld\n",ksm(i,p-2)%p);
	}
}

---

法二

解不定方程（扩欧）

解同余方程 $ax\equiv 1(modp)$等价于解不定方程 $ax+py=1$

根据裴蜀定理，这个不定方程有整数解的前提是gcd(a,p)=1，即a，p互质。（而保证互质的一个方法是令a或b为一个大质数，比如998244353或者1e9+7之类的，所以如果题目莫名给了一个大质数或者a，b互质的条件，很可能就是在保证扩欧方程有解的前提）

$Exgcd$求解即可（不会 请左转 Exgcd ）

code

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define ll long long

using namespace std;

ll n,p,x,y;

void exgcd(int a,int b)
{
	if(b==0)
	{
		x=1,y=0;
		return;
	}
	exgcd(b,a%b);
	int z=x;
	x=y,y=z-(a/b)*y;
}

int main(){
	scanf("%lld%lld",&n,&p);
	for(ll i=1;i<=n;i++)
	{
		exgcd(i,p);
		x=(x%p+p)%p;
		printf("%lld\n",x); 
	}
}

---

法三

线性递推：复杂度 $O(n)$！！！

过程：
令 $p=k\ast i+r$ ;$(k=p/i，r=pmodi)$，其中$(i<p,k<p,r<i)$，则有 ：

• $k\ast i+r\equiv 0(modp)$①

①式左右同乘$i^{-1}\ast r^{-1}$得：

• $k\ast r^{-1}+i^{-1}\equiv 0(modp)$移项 得

移项 得

• $i^{-1}\equiv -k\ast r^{-1}(modp)$

带入$k=p/i，r=pmodi$

得

• $i^{-1}\equiv (-\frac{p}{i})\ast (pmodi{)}^{-1}(modp)$②

由于 $(pmodi)<i$ 所以,在求出 $i^{-1}$之前,我们早已求出 $(pmodi{)}^{-1}$因此用数组 $inv[i]记录i^{-1}$ ( $i$ 的逆元)，则

• $inv[i]\equiv (-\frac{p}{i})\ast inv[pmodi]modp$

不要以为到这里就结束了

因为我们需要保证$i^{-1}>0$所以,我们在②式右边$+p(pmodp=0$ , 答案不变)

即 $inv[i]=p-\frac{p}{i}\ast inv[pmodi]modp$

当然 $inv[1]=1,inv[0]=tan90$(赋值为0);

而且$for$循环 要从$2$开始，防止改变 $inv[1]$ 的值;

code:

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define ll long long

using namespace std;

ll n,p,x,y,inv[3000010];

int main(){
	scanf("%lld%lld",&n,&p);
	inv[1]=1;
	printf("%lld\n",inv[1]);
	for(ll i=2;i<=n;i++)
	{
		inv[i]=p-(p/i)*inv[p%i]%p;
		printf("%lld\n",inv[i]);
	}
}

---

应用

一般用于求 $\frac{a}{b}modp$ ，把它转换为求a*（b的逆元）