扩展欧几里德算法

• 遇见这道题，对求最大公约数又有了新的认识

欧几里得问题

int gcd(int a, int b) {
    if (b == 0) 
        return a;
    else 
        return gcd(b, a%b);
}

扩展欧几里得问题

• 对于ax + by = c; 在 c % gcd(a,b) = 0的时候，我们一般先求 a/gcd(a,b) x + b/gcd(a,b) y = 1的解，然后在扩大c倍，就是原问题的解。当然c % gcd(a,b) ！= 0 的时候，原问题无解。
• 根据下面参考文献中的原理，我们可以得到下面的程序：

#include<iostream>
using namespace std;
void gcd(int a, int b, int &d, int &x, int &y) {
    if (!b) {
        d = a;
        x = 1;
        y = 0;
    }else {
        gcd(b, a%b, d, y, x);
        y -= x*(a/b);   //***
    }
}
int main() {
    int a, b, d, x, y;
    a = 6;
    b = 15; 
    gcd(a, b, d, x ,y);
    cout << x << " " << y << " " << d;
}

• 对于程序中加*号的地方，我的理解：
  
  • 由参考文献可知，x1 = y2; y1 = x2 - a/b * y2; 因为经过了一步，gcd(b, a%b, d, y, x); 所以 此时的 y = x2; x = y2; 所以根据将 y1 公式翻过去就是 y = y - x *(a/b);

参考

• http://blog.csdn.net/rising_fallmoon/article/details/10724239（这个是推导过程）