扩展欧几里得

首先接触的就是欧几里得求最大公约数（GCD）

int gcd(int a, int b){
    // b < a
    while(b){
        int t = a%b;
        a = b;
        b = t;
    }
    return a;
}

递归的写法：

int gcd(int a, int b){
   if(!b) return a;
   return gcd(b, a%b);
}

 

在扩展欧几里得中，设法求 ax + by = gcd(a, b) 中的一组 x 和 y 的解

当 b = 0 的时候， gcd(a, b) = a,此时的 x = 1, y = 0为其中的一组的解

当 b ≠ 0 的时候　， 设 

　　ax1 + by1 = gcd(a, b)

又辗转相除法得：

　　gcd(a, b) = gcd(b, a%b)

设 bx2 + (a%b)y2 = gcd(b, a%b) = gcd(a, b)

所以综合得到一个等式：　　ax1 + by1 = bx2 + (a%b)y2　　(想得到x1, y1)

又 a%b = a - (a/b)*b ，带入上式中

得　　bx2 + [ a - (a/b)*b]y2 = ax1 + by1

　　　ay2 + bx2 - (a/b)*b*y2 = ax1 + by1

所以在知道x2, y2 之后，便可以推出

　　x1 = y2

　　y1 = x2 - (a/b)*y2

然后推回去就行了

int ex_gcd(int a, int b, int &x, int &y){
    if(b == 0){
        x = 1;
        y = 0;
        return a;
    }
    int r = ex_gcd(b, a%b, x, y);
    int t = x;
    x = y;
    y = t - floor(a/b)*y;
    return r;
}