int gcd(int a,int b){return a%b==0? b:gcd(b,a%b);}
首先对于基础gcd,基于性质gcd(a,b)=gcd(b,a%b)
int x,y;
int exgcd(int a,int b)
{
if(b == 0)
{
x=1;
y=0;
return a;
}
int ans = exgcd(b,a%b);
int t=y;
y = x - a/b*t;
x = t;
return ans;
}
对于exgcd,基于以下证明
关于求逆元
求a的逆元相当于ax=1(mod m),即ax-my = 1,调用exgcd(a,m)解出x即为a的逆元
如果求m=(a/b)(mod p)
即求( a / b ) % p =a * inv ( b , p ) %p =( a%p * inv ( b , p )%p ) %p