一、正常求逆元(扩栈欧几里得)
//a对b求逆元
ll exgcd(ll a,ll b,ll &x,ll &y)
{
ll t,gcd;
if(b==0)
{
x=1,y=0;
return a;
}
gcd=exgcd(b,a%b,x,y);
t=x,x=y,y=t-a/b*y;
return gcd;
}
ll inv(ll a,ll b)
{
ll x,y;
exgcd(a,b,x,y);
x%=b;
x=(x+b)%b;
return x;
}
二、以下方法要求模为素数
1.
//a对mod的逆元,mod为质数
ll inv(ll a,ll mod)
{
return a == 1 ? 1 : (mod - mod / a) * inv(mod % a,mod) % mod;
}
2.
ll inv(ll n,ll p)
{
ll m=p-2,ans=1;
while(m)
{
if(m&1) ans=ans*n%p;
m>>=1;
n=n*n%p;
}
return ans;
}
3.
ll inv[MAXN];
inv[1] = 1;
for (int i = 2; i<MAXN; i++)
inv[i] = inv[MOD%i]*(MOD-MOD/i)%MOD;