扩展欧几里得

推理过程：

$ax+by=gcd(a,b)$
求解 x，y的方法的理解
设 a>b。
①显然当 $b=0$，$gcd(a,b)=a$。此时 x=1，y=0；
②$a>b>0$ 时
设
$ax_1+ by_1= gcd(a,b)$
$bx_2+ (a\mod b)y_2= gcd(b,a\mod b)$;
根据朴素的欧几里德原理有
$gcd(a,b) = gcd(b,a\mod b)$;
则:$ax_1+ by_1= bx_2+ (a mod b) y_2$;
即:$ax_1+ by_1= bx_2+ (a - [\frac{a}{b}] * b)y_2=ay_2+ bx_2- [\frac{a}{b}] * by_2$
也就是$ax_1+ by_1 == ay_2+ b(x_2- [\frac{a}{b} y_2)$;
根据恒等定理得：$x1=y2$ ; $y1=x2- [a / b] *y2$;
这样我们就得到了求解 $x_1,y_1$ 的方法：
$x_1，y_1$ 的值基于 $x_2，y_2$.
上面的思想是以递归定义的
∵ gcd 不断的递归求解一定会有个时候 $b=0$
∴递归可以结束。

代码：

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