扩展欧几里得算法 exgcd 求逆元（适用于模数不为质数的情况）

原理

不打算自己懂。。。

代码

ull exgcd(ull a,ull b,ull &x,ull &y)//扩展欧几里得求模b意义下a的逆元
//返回的d是a和b的最大公约数,而最终的x是a在模b意义下的逆元 
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }    
    ull d=exgcd(b,a%b,y,x);
    y=y-a/b*x;
    return d;
}

exgcd(a,b,x,y);//注意最终x可能返回负数,要加上b变成正数再取模 
        cout<<(x+b)%b<

注意最终的x可能返回的是负数，所以要加上b就成正数再取模。