拓展欧几里德算法

        拓展欧几里德定理：对于不完全为0的非负整数a，b，必然存在整数对x，y，使得ax+by=gcd(a,b)。

        利用拓展欧几里德定理，可以对形如ax+by=c的式子进行是否存在整数解的判断，以及存在时解为何值的运算，整个过程称为拓展欧几里德算法。

        首先，对于ax+by=c这一式子，根据拓展欧几里得定理ax+by=gcd(a,b)，若c%gcd(x,y)=0成立，则表示式子存在整数解，此时，我们有:

                ax₁+by₁=gcd(a,b)=gcd(b,a%b)=bx₂+(a%b)y₂=bx₂+(a-a/b*b)y₂=ay₂+bx₂-b(a/b)y₂

        根据恒等定理，显然有x₁=y₂,y₁=x₂-(a/b)y₂。从中可以看出x₁和y₁的值取决于x₂和y₂，以此类推，可知x₂和y₂的值取决于x₃和y₃……直到x_n和y_n，此时b=0，使得ax_n=gcd(a,0)=a，故x_n=1,y_n=0，再往上逆推x_n-1和y_n-1等等，直到x₁和y₁的值，整个过程可以用递归来实现。此时，需要注意到ax₁+by₁=gcd(a,b)，(c/gcd(a,b))x₁+(c/gcd(a,b))y₁=c，而且这只是这个方程的某一对解而已，因为我们令y_n=0了。

         另外，可以留意到的是，求解x和y的过程中也是可以顺便求出gcd(a,b)的，至于代码则可以参考下面贴出来的：

//对于ax+by=gcd(a,b)，求x和y的值
void exEuclid( int a, int b, int& x, int& y, int& gcd ) {
    if ( b == 0 ) {
        x = 1;
        y = 0;
        gcd = a;
    }
    else {
        exEuclid( b, a % b, y, x, gcd );
        y -= a / b * x;
    }
}