gcd(a,b)=gcd(b,a%b)
gcd(a,b)中要求a,b是非负数

//xa+yb=gcd
int exgcd(int a,int b,int &x,int &y){
    if(!b){
        x=1; y=0;
        return a;
    }
    int d=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
ax+by=d通过exgcd求出了x0,y0(ax+by=gcd的解)后求出x和y的通解:
x=(x0-b/gcd)*k
y=(y0+a/gcd)*k

另类的方法(不常用,不实用)

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