关于扩展欧几里得算法（Extended Euclidean Algorithm），我是在做青蛙的约会这一经典题目才接触到这个算法的。后面也有关于这一题的AC代码和解题思路。
内容：已知a, b，求解一组x，y，使它们满足贝祖等式： ax+by =gcd(a, b)
扩展欧几里得算法，就和它的名字一样是对欧几里得算法的扩展。何为扩展？一是，该算法保留了欧几里得算法的本质，可以求a与b的最大公约数。二是，已知a, b求解二元一次方程ax+by =gcd(a, b)的一组解（x，y）。
证明：
假设 a>b,
(1) b=0 gcd(a,b) = a , ax = a , 则x=1，y=0;（这里我还是推荐不把gcd(a,0)理解成最大公约数，而是一个计算机求出来的值）
(2) 假设 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，即ax1+by1=bx2+(a-a/b*b)y2 ax1+by1=ay2+b(x2-a/b*y2)
在根据多项式恒等定理（把a,b看成变量），x1=y2; y1=x2-a/b*y2;
（表面上看，就是已知方程一的一组解，可以得到方程二的一组解，已知方程二的一组解，就可以得到方程一的一组解，但是实际情况是，不可能先知道方程一的解（x1,y1）。）上述思想是递归定义的，不断地利用gcd(a,b) =gcd(b,a%b)，到b=0（y的系数为0）时，由(1)的解，根据解之间的关系，最终可以得到方程ax+by =gcd(a, b)的解。
递归形式代码：
1 #include<iostream> 2 using namespace std; 3 4 int exgcd(int a,int b,int &x,int &y) 5 { 6 if(b==0) 7 { 8 x=1; 9 y=0; 10 return a; 11 } 12 int gcd=exgcd(b,a%b,x,y); 13 int x2=x,y2=y; 14 x=y2; 15 y=x2-(a/b)*y2; 16 return gcd; 17 } 18 19 int main() 20 { 21 int x,y,a,b; 22 cout<<"请输入a和b:"<<endl; 23 cin>>a>>b; 24 cout<<"a和b的最大公约数:"<<endl; 25 cout<<exgcd(a,b,x,y)<<endl; 26 cout<<"ax+by=gcd(a,b) 的一组解是:"<<endl; 27 cout<<x<<" "<<y<<endl; 28 return 0; 29 }
非递归形式代码：
1 #include<iostream> 2 using namespace std; 3 4 int exgcd(int a,int b,int &x,int &y) 5 { 6 int x1,y1,x0,y0; 7 x0=1; y0=0; 8 x1=0; y1=1; 9 x=0; y=1; 10 int r=a%b; 11 int q=(a-r)/b; 12 while(r) 13 { 14 x=x0-q*x1; y=y0-q*y1; 15 x0=x1; y0=y1; 16 x1=x; y1=y; 17 a=b; b=r; r=a%b; 18 q=(a-r)/b; 19 } 20 return b; 21 } 22 23 int main() 24 { 25 int x,y,a,b; 26 cout<<"请输入a和b:"<<endl; 27 cin>>a>>b; 28 cout<<"a和b的最大公约数:"<<endl; 29 cout<<exgcd(a,b,x,y)<<endl; 30 cout<<"ax+by=gcd(a,b) 的一组解是:"<<endl; 31 cout<<x<<" "<<y<<endl; 32 return 0; 33 }
运行截图：
同样有两点想说明：
1.扩展欧几里得算法是对欧几里得算法的扩展，可以求出gcd(a,b)，好多人都没意识到这一点。
2.x,y可以用全局变量，参数传递就不用传引用了。