扩展欧几里得算法

裴蜀定理

裴蜀定理：对任意的整数a、b，关于未知数x、y的线性丢番图方程ax+by=m，当且仅当m是gcd(a,b)的整数倍时方程有解。
该方程称为裴蜀等式。

证明$ax+by=gcd(a,b)$时有解：
设$d=gcd(a,b)$，则有$d\mid a$，$d\mid b$，因此$\forall x,y\in \mathbb{Z},d\mid ax+by$，即$d\mid m$
设$s$是最小的$m$，则有$d\mid s$
令$q=\lfloor \frac{a}{s}\rfloor ,r=a\%s=a-qs=a-q(ax+by)=a(1-qx)+b(-qy)$，可见$r$也是$x,y$的线性组合，因为$s$是$x,y$线性组合的最小值，$0\le r<s$，所以$r=0$，则有$s\mid a$
同理有$s\mid b$，因此$s$是$a,b$的公约数，而$d$是$a,b$的最大公约数，因此有$s\mid d$
综上，$s=d$

证明$ax+by=m$有解的充要条件是$gcd(a,b)\mid m$：
充分性：若$gcd(a,b)\mid m$，则$m=k\ast gcd(a,b)$，设$ax+by=gcd(a,b)$的一组解为$x_0,y_0$，显然$k{x}_0,k{y}_0$是$ax+by=m$的一组解。
必要性：设$d=gcd(a,b)$，则有$d\mid a$，$d\mid b$，因此$\forall x,y\in \mathbb{Z},d\mid ax+by$，即$d\mid m$

扩展欧几里得算法

显然，当m等于k倍的gcd(a,b)时，裴蜀等式的解是m恰好等于gcd(a,b)时解的k倍。
因此，对任意的裴蜀等式，我们只需考察m恰好等于gcd(a,b)时裴蜀等式的解。
扩展欧几里得算法就是用来求当m恰好等于gcd(a,b)时裴蜀等式的解的一种算法。

设$a\ast x_1+b\ast y_1=gcd(a,b)，b\ast x_2+(a\%b)\ast y_2=gcd(b,a\%b)$
由欧几里得算法（辗转相除法）可知，$gcd(a,b)=gcd(b,a\%b)$
$\therefore a\ast x_1+b\ast y_1=b\ast x_2+(a\%b)\ast y_2$
又$\because a\%b=a-\lfloor \frac{a}{b}\rfloor \ast b$
$\therefore a\ast x_1+b\ast y_1=b\ast x_2+a\ast y_2-\lfloor \frac{a}{b}\rfloor \ast b\ast y_2$
移项得，$a\ast (x_1-y_2)+b\ast (y_1-x_2+\lfloor \frac{a}{b}\rfloor \ast y_2)=0$
因此，$x_1=y_2,y_1=x_2-\lfloor \frac{a}{b}\rfloor \ast y_2$就是一组解

根据上面的推导可知：对于一组线性方程ax+by=gcd(a,b)，它的一组解$x_1,y_1$可以通过另一组线性方程bx+(a%b)y=gcd(b, a%b)的一组解$x_2,y_2$计算得出。这恰好满足递归的条件。
例如，当求(a,b)为(15,8)的解时，我们需要先计算(8,7)的解，进而需要先计算(7,1)的解，再进而需要先计算(1,0)的解。
显然，递归的边界是b=0
当b等于0时，方程变为ax=gcd(a,0)，又因为gcd(a,0)=a，所以方程变为ax=a，解为x=1
这样，我们运用递归算法就可以求出任意一组线性方程ax+by=gcd(a,b)的一组解：
先不断递推，一定能推至b=0，b=0时解是已知的，回溯时用已经知道的$x_2,y_2$，计算上一层的$x_1,y_1$,这样就可以算出最开始的解。
代码如下：

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

该算法在递归中巧妙地运用了引用类型，使得代码无比简洁。
由于x，y参数是引用类型，因此在若干层递归的过程中，我们改变的一直是同一块内存中存储的数值。

exgcd(b, a % b, y, x);

这行代码是整个算法的核心。根据前文的推导，将a、b变为b、a%b很好理解，那么为什么要颠倒x，y的次序呢？
还是根据前文的推导$x_1=y_2$，也就是说当前层的x和下一层的y相等，所以我们只需要将当前层的x引用传给下一层的y引用。引用可以理解为一个变量的别名，颠倒x，y的次序可以理解成：当前层的x和下一层的y是同一个变量的不同名称，他们其实是一个变量。 这样就做到了递归回溯后$x_1=y_2$
那么$y_1=x_2-\lfloor \frac{a}{b}\rfloor \ast y_2$是如何做到的呢？
由于x，y在递归的上下层中次序颠倒，使得回溯后$x_1=y_2$，同样也有$y_1=x_2$
即在执行这行代码后$x_1=y_2$，$y_1=x_2$
那么$y_1=x_2-\lfloor \frac{a}{b}\rfloor \ast y_2$就变成了$y_1=y_1-\lfloor \frac{a}{b}\rfloor \ast x_1$
于是就有了下面一行代码

y -= a / b * x; 

要注意的是，裴蜀等式有解时一定有无穷多的解，这个算法求得的只是其中一组解。
代码应用如下：

#include<iostream>
using namespace std;
void exgcd(int a, int b, int& x, int& y)
{
	if (b == 0) x = 1;
	else
	{
		exgcd(b, a % b, y, x);
		y -= a / b * x;
	}
}
int main()
{
	int a, b, x, y;
	x = 0;
	y = 0;
	cin >> a >> b;
	exgcd(a, b, x, y);
	cout << x << " " << y << endl;
	return 0;
}

最小正整数解

int min_x = (x % b + b) % b;//x的最小正整数解
int min_y = (y % a + a) % a;//y的最小正整数解
cout << min_x << " " << (x - min_x) / b * a + y << endl;//x的最小正整数解与其对应的y
cout << (y - min_y) / a * b + x << " " << min_y;//y的最小正整数解与其对应的x

证明：
不妨设$x_0$为$x$的最小正整数解，$y_0$为$x_0$对应的$y$
$\because a\ast x+b\ast y=gcd(a,b)$
$\therefore a\ast x+a\ast b\ast k+b\ast y-a\ast b\ast k=gcd(a,b)$
合并得，$a\ast (x+k\ast b)+b\ast (y-k\ast a)=gcd(a,b)$
因此，所有整数解都可写为$x=x_0+k\ast b,y=y_0-k\ast a$的形式，即x多一个b，y就要少一个a
显然，若$x$为正数，$x\%b=(x_0+k\ast b)\%b=x_0$
若x为负数，$x\%b=[(x_0-b)+(k+1)\ast b]\%b=x_0-b$
$\therefore x\%b+b=x_0或x_0+b$
$\therefore (x\%b+b)\%b=x_0$
综上所述$x$的最小正整数解为$(x\%b+b)\%b$
同理可证$y$的最小正整数解为$(y\%a+a)\%a$
注意：一般情况下，x的最小正整数解与y的最小正整数解并不是一组解。