乘法逆元

乘法逆元：若a*x ≡1（mod p），且gcd（a,p）=1，则称a关于模p的乘法逆元为x。

求解方法

1.拓展欧几里得

exgcd用以求解ax+by=gcd(a,b)式的问题，当gcd（a,b)=1时， 令p=b，原式化为ax+py=1，即a*x ≡1（mod p）

#include<iostream>
#include<cstring>
using namespace std;

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

int inv(int a,int p)
{
	int x,y;
	exgcd(a,p,x,y);
	return (x%p+p)%p;
}

int a,p;

int main()
{
	cin>>a>>p;//a*b=1(modp)
	cout<<inv(a,p)<<endl;
	return 0;
}

2.费马小定理
∀a,p∈N,若p为质数，且p ∤ a，有a^(p-1)与1同余p

p需为质数

#include<iostream>
#include<cstring>
using namespace std;

int fpow(int a,int b,int p)//a^b % p 
{
	if(!b)return 1;
	int ans=1;
	if(b&1){b--;ans=a;}
	a=((a%p)*(a%p))%p;
	b>>=1;
	return ((fpow(a,b,p)%p)*(ans%p))%p;
}

int inv(int a,int p)
{
	return fpow(a,p-2,p);//a^(p-2) % p;
}

int a,p;

int main()
{
	cin>>a>>p;//a^(p-1)=1(modp)--->a^(p-2)=a^(-1) (mod p)
	cout<<inv(a,p)<<endl;
	return 0;
}

3.欧拉定理

#include<iostream>
#include<cstring>
using namespace std;

int fpow(int a,int b,int p)//a^b % p 
{
	if(!b)return 1;
	int ans=1;
	if(b&1){b--;ans=a;}
	a=((a%p)*(a%p))%p;
	b>>=1;
	return ((fpow(a,b,p)%p)*(ans%p))%p;
}

int phi(int x)
{
	int ans=x;
	for(int i=2;i*i<=x;i++)
	{
		if(x%i==0)ans=ans-ans/i;
		while(x%i==0)x/=i;
	}
	if(x>1)ans=ans-ans/x;
	return ans;
}

int inv(int a,int p)
{
	return fpow(a,phi(p)-1,p);//a^(phi(n)-1) % p;
}

int a,n;

int main()
{
	cin>>a>>n;//a^(phi(n))=1 (mod p )
	cout<<inv(a,n)<<endl;
	return 0;
}

4.线性递推

初始化inv[0]=inv[1]=1;则有递推式inv[i]=(p-p/i)*inv[p%i]%p。
p需为质数

#include<iostream>
#include<cstring>
using namespace std;

int n,p,inv[100010];

int main()
{
	cin>>n>>p;
	inv[1]=1;
	for(int i=2;i<=n;i++)
	{
		inv[i]=(p-p/i)*inv[p%i]%p;
	}
	cout<<inv[n]; 
	return 0;
}