扩展欧几里德算法

在学习扩展欧几里德算法之前，先来回顾一下欧几里德算法。欧几里德算法又称辗转相除法，它是用来求两个数的最大公约数其依赖于一个定理：gcd(a,b) = gcd(b,a%b);

欧几里德算法：

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

扩展欧几里德算法是用来在已知a,b求解一组x，y使得a*x+b*y=gcd(a,b)(根据数论中不定方程解的定理，可知解一定存在)，扩展欧几里得算法常用在求解模线性方程及方程组中。

下面给出证明过程：

根据欧几里德定理：gcd(a,b) = gcd(b,a%b);我们将gcd(b,a%b)中的b换为a'，a%b换为b'，即a'=b，b' = a%b使得a'*x+b'*y=gcd(a',b')；由于 b' = a% b = a - a / b * b，那么可得到如下变换：a'*x+b'*y=gcd(a',b')  ===》bx + (a - a / b * b)y = gcd(a', b') = gcd(a, b)  ===》ay +b(x - a / b*y) = gcd(a, b),因此对于a、b而言，他们对应未知量的解分别为：y和(x - a / b*y) 

下面是不定方程的求解过程，求解a*x+b*y=n；

1、先计算gcd(a,b)，若n不能被gcd(a,b)整除，则方程无整数解；否则，在方程两边同时除以gcd(a,b)，得到新的不定方程a'* x + b' * y = n'，此时gcd(a',b')=1;

2、利用上面所说的欧几里德算法求出方程a' * x + b' * y = 1的一组整数解x0,y0，则n' * x0,n' *y0是方程a' * x + b' * y = n'的一组整数解；

3、根据数论中的相关定理，可得方程a'* x + b' * y = n'的所有整数解为：x = n' * x0 + b' * t          y = n' * y0 - a' * t,其中n' = n/gcd(a,b);  b' = b/gcd(a,b), a' = a/gcd(a,b);

故上式可化为：x = n' * x0 + b' * t = n/gcd(a,b) *x0 + b/gcd(a,b) *t；

                       y = n' * y0 - a' * t = n/gcd(a,b) * y0 - a/gcd(a,b) * t；（t为整数）

上式就是不定方程a*x+b*y=n的所有整数解；

下面给出扩展欧几里得算法代码：

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

返回值为最大公约数

求解不定方程的一组解：

//d为a，b的最大公约数，不必自行给出，x，y也不必自行给出，只需在主函数中定义就行了
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);
    }
}

如果说以上写的还是没看懂，那么这里给出一个个人认为比较好的博客链接看下大牛们的也许就能理解了：

http://www.cnblogs.com/frog112111/archive/2012/08/19/2646012.html