裴蜀定理（扩展欧几里得算法）

一、内容

$\forall a,b\in Z$, 则$\exists x,y\in Z$, 满足方程
$$ax+by=gcd(a,b)$$

ps : 该方程有多组解x,y。

---

二、通解推导

对于方程
$$a{x}_1+b{y}_1=gcd(a,b),a{x}_2+b{y}_2=gcd(a,b)$$
将两式相减：
$$a(x_1-x_2)=b(y_2-y_1)$$
等式两边同时除以$gcd(a,b)$：
$$a^{\prime }(x_1-x_2)=b^{\prime }(y_2-y_1)$$
从而$b^{\prime }|a^{\prime }(x_1-x_2)$。又因为$a^{\prime },b^{\prime }$互质，所以$b^{\prime }|(x_1-x_2)$。令$(x_1-x_2)/b^{\prime }=k$，则
$$x_1=x_2+k\ast b^{\prime }=x_2+k\ast \frac{b}{gcd(a,b)}$$
同理可以得到：
$$y_1=y_2-k\ast \frac{a}{gcd(a,b)}$$

---

所以最后得到通解为：
$$\{\begin{array}{l}X=x+k_1\ast \frac{b}{gcd(a,b)}\\ Y=y-k_2\ast \frac{a}{gcd(a,b)}\end{array}$$
其中x, y为扩展欧几里得算法求出来的一组解。

---

三、算法推导

若
$$a{x}_1+b{y}_1=gcd(a,b)\text{\hspace{1em}(1)}$$
则用$b$替换$a$, $a\%b$替换$b$, 得到
$$b{x}_2+(a\%b)y_2=gcd(b,a\%b)\text{\hspace{1em}(2)}$$
如果找到了$x_1,x_2$的关系以及$y_1,y_2$的关系, 就可以利用递归来计算出原方程中的解, 即$x_1,y_1.$

我们将$(2)$进行变形：
$$b{x}_2+(a-\lfloor \frac{a}{b}\rfloor b)y_2=gcd(b,a\%b)$$
合并同类项得到：
$$a{y}_2+b(x_2-\lfloor \frac{a}{b}\rfloor y_2)=gcd(b,a\%b)\text{\hspace{1em}(3)}$$

由于$gcd(a,b)=gcd(b,a\%b)$, 所以由$(2),(3)$得：
$$a{x}_1+b{y}_1=a{y}_2+b(x_2-\lfloor \frac{a}{b}\rfloor y_2)\text{\hspace{1em}(4)}$$
对比系数可以得到：

$$\{\begin{array}{l}{x}_1=y_2\\ y_1=x_2-\lfloor \frac{a}{b}\rfloor y_2\end{array}\text{\hspace{1em}(5)}$$

因此，我们可以先计算出$x_2,y_2$, 然后推出原方程的解$x_1,y_1$。显然，这是一个递归的过程。

---

如果我们用$a\%b$替换a, $b$替换$b$, 则$(2)$变成了
$$b{y}_2+(a\%b)x_2=gcd(a\%b,b)\text{\hspace{1em}(6)}$$
则最后的推导结果为：
$$\{\begin{array}{l}{x}_1=x_2\\ y_1=y_2-\lfloor \frac{a}{b}\rfloor x_2\end{array}\text{\hspace{1em}(7)}$$

---

四、代码部分

/* 采用第二种推导结果 */
int exgcd(int a, int b, int &x, int &y) { // 返回两个值x y，因此传引用
	/* 由辗转相除法可知：递归到最后一定是 ax + 0y = gcd(a, 0) = a
	   所以递归结束的条件为 b == 0，而此时x = 1，y = 0 */
	if(b == 0) {
		x = 1;
		y = 0;
		return a;
	}
	int d = exgcd(b, a % b, y, x); // 参数对应方程 : by + (a % b)x = gcd(b, a%b)
	y -= a / b * x;
	return d;
}

五、拓展

判断下列不定方程解的情况
$$ax+by=c\text{\hspace{1em}(8)}$$
先说结论：若$gcd(a,b)|c$，则方程一定有解；否则一定无解

---

1. $c\%gcd(a,b)!=0$

由于等式左边
$$gcd(a,b)|a,gcd(a,b)|b$$
所以方程一定无解

2. $c\%gcd(a,b)==0$

将方程$(1)$两边同时乘上$\frac{c}{gcd(a,b)}$，即可得到方程$(8)$的解。所以此时方程一定有解。