【算法】乘法逆元

定义

逆元是广义化的倒数，在除法运算当中，我们经常会遇到 $÷a$，那么 $\frac{1}{a}$ 或 $a^{-1}$ 就是它的逆元，通常用 $inv(a)$ 来表达。

换句话说，在同余意义下，除以一个数等于乘上一个数的逆元。

若某一实数 $x$ 满足

$$ax\equiv 1(\mathrm{mod}p)$$

我们就称 $x$ 是 $a$ 在模 $p$ 意义下的逆元。

任何整数 $a\mathrm{mod}p$ 都有逆元 $a^{-1}\mathrm{mod}p$

同时，逆元还有一些性质

• 对于一个数在模 $p$ 意义下的逆元，必定唯一

• $inv(a)\times inv(b)=inv(a\times b)$

扩展欧几里得求逆元

我们将此同余方程化简一下

$$ax-py=1$$

我们发现，当前不定方程与扩欧十分相像，如果令 $a=a$，$b=p$ 似乎可以直接求出 $x$ 即逆元，可是等号右边却不同。

根据扩欧的性质可以得出（不知道的去看博客），若 $\gcd (a,b)\mid c$，则方程一定有解。那么这时，如果 $a$ 和 $p$ 不互质，也就没有逆元了，所以 $a$ 和 $p$ 一定互质。

随后就可以正常的求扩欧了，得出 $x$ 直接输出即可。

代码

#include<bits/stdc++.h>
using namespace std;
int a,p,x,y;
int exgcd(int a,int b,int &x,int &y)
{
	if(!b)
	{
		x=1;y=0;
		return a;
	}
	int gcd=exgcd(b,a%b,x,y);
	int t=x;
	x=y;
	y=t-(a/b)*y;
	return gcd;
}
int main()
{
	scanf("%d %d",&a,&p);
	int gcd=exgcd(a,p,x,y);
	if(gcd!=1)	
		puts("impossible");
	else printf("%d",(x%p+p)%p);
	return 0;
}

快速幂求逆元

我们知道，费马小定理所证明的式子为

$$a^{p-1}\equiv 1(\mathrm{mod}p)$$

（不会的看博客）

我们可以将它化简一下，即：

$$a^{p-2}\times a\equiv 1(\mathrm{mod}p)$$

显然，此时 $a^{p-2}$ 是 $a$ 在模 $p$ 意义下的逆元。

代码

#include<bits/stdc++.h>
using namespace std;
int a,p;
int gcd(int a,int b){return b?gcd(b,a%b):a;}
int Pow(int n,int k)
{
	int res=1;
	while(k)
	{
		if(k&1)
			(res*=n)%=p;
		(n*=n)%=p;
		k>>=1;
	}
	return (res%p+p)%p;
}
int main()
{
	scanf("%d %d",&a,&p);
	int g=gcd(a,p);
	if(g!=1)	
		puts("impossible");
	else printf("%d",Pow(a,p-2));
	return 0;
}

递推求逆元

我们设模数 $p=k\times i+r$，其中 $i$ 为除数，$k$ 为 $p÷i$ 的商，$r$ 为余数。再将其乘上 $i$ 和 $r$ 的逆元，即:

$$k\times r^{-1}+i^{-1}\equiv 0(\mathrm{mod}p)$$

$$i^{-1}\equiv -k\times r^{-1}(\mathrm{mod}p)$$

再将其替换成我们已知的变量

$$i^{-1}\equiv -\lfloor \frac{p}{i}\rfloor \times (p\mathrm{mod}i{)}^{-1}(\mathrm{mod}p)$$

我们即可从小到大枚举 $i$，递推求出逆元。

代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
int n,p;
int inv[13000001];
signed main()
{
	scanf("%lld %lld",&n,&p);
	inv[1]=1;
	puts("1");
	for(int i=2;i<=n;i++)
	{
	    inv[i]=((((-p/i))*inv[p%i])%p+p)%p;
	    printf("%lld\n",inv[i]);
	}
	return 0;
}