扩展欧几里得算法

如果事先不了解 欧几里得算法 ，请点击。

扩展欧几里得算法：对于不完全为0的非负整数 a,b，必然存在整数对 X,Y，使得 aX + bY = gcd(a,b)。

解法见注释。

/*    How to solve "aX + bY = gcd(a,b)" ?1、if b=0,     gcd(a,b) = a,     X = 1 , Y is any number.     For convenience,we make Y equals 0.2、if b!=0,    gcd(a,b)   = aX1 + bY1;    gcd(b,a%b) = bX2 + (a%b)Y2;                 ----> (a-a/b*b)Y2;--> aX1 + bY1  = aY2 +  b(X2-a/b*Y2);--> X1 = Y2 , Y1 = X2-a/b*Y2;    So we can use recursion to sovle it.*/#include <stdio.h>int GCD;void ExGcd(int a,int b,int &x,int &y){    if(b == 0)    {        x = 1;        y = 0;        GCD = a;        return ;    }    ExGcd(b,a%b,x,y);    int t = x;    x = y;    y = t - a/b*y;}int main(){    int a,b,x,y;    scanf("%d%d",&a,&b);    ExGcd(a,b,x,y);    printf("gcd(%d,%d) -> %d * %d + %d * %d\n",a,b,a,x,b,y);return 0;}

下面讨论如何用扩展欧几里得算法解决不定方程 aX + bY = c。

1、当 c mod gcd(a,b) = 0 的时候，该方程才有整数解。

2、方程可转化为 aX + bY = gcd(a,b) * c/gcd(a,b)。

等价于 aX/(c/gcd(a,b)) + bY/(c/gcd(a,b)) = gcd(a,b)。

所以将方程 aX + bY = gcd(a,b) 的解乘以 c/gcd(a,b) 即可得到解 X0,Y0。

对于上述两个方程的其他解，我们可以用方程 a(X0+△X) + b(Y0+△Y) = c 表示。

于是得到 a△X + b△Y = 0，即 △X / △Y = - b/a = - ( b/gcd(a,b) ) / ( a/gcd(a,b) ) = t。( t 为与X、Y无关的参数)

于是方程的所有解可以用 X = X0 + b/gcd(a,b)*z ， Y = Y0 - a/gcd(a,b)*z 来表示。（ z 为任意整数）

另外有一个没有证明过的猜想：利用上述方法求出的 (X0,Y0) 是所有 |X| + |Y| 的值中最小的（并不一定唯一），从而有X0 < |b|。（假设a < b）