扩展GCD gcd解法

扩展gcd要解决的问题：已知ax+by=gcd(a,b)，求一组可行的x,y。

①由gcd(a,b)gcd(b,a%b)得ax+bybx+(a%b)y，继续化下去，化成ax1+by1的形式即为递归转移式，因为递归过程中a,b不能改变，改的是x,y。

bx+(a%b)y == bx+ (a-(a/b)b) y == bx+ay-(a/b)by == ay + b(x-a/by)

即exgcd(ax+by)== exgcd( ay + b(x-a/b*y) )

然后就可以递归求解了，共用相同的x,y，可以全局，下面代码是用引用实现。

②递归终止：当b0时，式子为axgcd(a,b)，那么显然a=gcd(a,b),x=1,y=0;

//扩展gcd求解x,y同时还把gcd(a,b)求出来了，exgcd的返回值就是gcd(a,b)

C++参考代码：

ll exgcd(ll a, ll b, ll &x, ll &y)  
{
    if(b == 0)
    {
        x = 1;
        y = 0;
        return a;
    }
    ll d = exgcd(b, a % b, x, y);
    ll t = x;
    x = y;
    y = t - (a / b) * y;
    return d;
｝

扩展gcd求出一组解(x0,y0)后，那ax+by=d的其他解为x=x0+kb/d; y=y0-ka/d , 其中d=gcd(a,b)

扩展gcd求出一组解(x0,y0)后，那ax+by=cd的一组解为x=cx0，y=c*y0

ax+by=c无整数解的情况：c不能被gcd(a,b)整除（也即c不是gcd(a,b)的倍数）。

理由：因为ax与by均是gcd(a,b)的倍数,所以其和ax+by也是gcd(a,b)的倍数。