必存在一对整数(x,y) 使得a*x+b*y=gcd(a,b);

gcd(a,b)= gcd(b,a%b) (辗转相除法)

又 a*x1+b*y1=gcd(a,b)=gcd(b,a%b)=b*x2+(a%b)*y2;

又b*x2+(a%b)*y2=b*x2+(a-[a/b]*b)*y2=a*y2-b*(x2-[a/b]*y2);

对比系数得 x1=y2,y1=x2-[a/b]*y2;

按以上递归得最后必有x*a=a(这个a表示上面的a和b的最大公约数); 即x=1
void exgcd(int a,int b,int &d,int& x,int &y){
    if(!b){
        d=a;x=1;y=0;
        return ;
    }
    exgcd(b,a%b,d,y,x);
    y-=a/b*x;
}


对于上诉代码 调换了x,y 也就导致了 y1=x2,并且y2=x1; 所以有 y-=x*(a/b);

用d记录最小公倍数,如果a*x+b*y=c,c != k*d 则无整数解