欧几里得算法(gcd)及其扩展

欧几里得算法(gcd)即辗转相除法，求两个数的最大公约数
辗转相除法的关键在于恒等式gcd ( a,b ) = gcd ( b,a % b ),以及边界条件gcd ( a ,0 ) = a
gcd函数代码如下

int gcd(int a,int b){
    return b == 0 ? a : gcd(b,a % b);
}

扩展的欧几里得算法用于解决如下问题：
已知a，b，求解一组x，y，使其满足 ax+by=gcd(a,b) （该等式称贝祖等式）
证明：
假设a>b
(1)当b=0，gcd(a,0)=a,则方程变为ax=a，所以一组解为x=1，y=0
(2)当b!=0，假设
ax1 + by1 = gcd(a,b)      
bx2 + (a%b)y2 = gcd(b,a%b)
因为gcd(a,b)  = gcd(b,a%b)，所以ax1 + by1 = bx2 + (a%b)y2    
从编程的角度上来说，a%b = a - (a/b)*b
所以上式等价变形为
ax1 + by1 = bx2 + [a - (a/b)*b]y2
                 = bx2 + ay2 - (a/b)*by2
                 = ay2 + b[x2 - (a/b)y2]
将a，b看做变量，根据对应项系数相等，有
x1