扩展欧几里得算法

裴蜀定理

证明

若 $a,b\in \mathbb{Z}$，那么对于 $\forall x,y\in \mathbb{Z}$，$\gcd (a,b)\mid a\times x+b\times y$。

此外，一定 $\exists x,y\in \mathbb{Z}$，使得 $a\times x+b\times y=\gcd (a,b)$ 成立。

证明

证明第一点：

$$\begin{array}{l}\because \gcd (a,b)\mid a,\gcd (a,b)\mid b\\ \therefore \gcd (a,b)\mid ax,\gcd (a,b)\mid by(x,y\in \mathbb{Z})\\ \therefore \gcd (a,b)=a\times x+b\times y\end{array}$$

证明第二点：

设 $a\times x+b\times y$ 的最小正整数值为 $r$，考虑证明 $a\mathrm{mod}r=0$。

设 $p=\lfloor \frac{a}{r}\rfloor$，则可得 $a\mathrm{mod}r=a-p\times r$​。

代入 $r$，得到 $a\mathrm{mod}r=a-p\times (a\times x+b\times y)$。

进一步可得 $a\times (1-p\times x)+b\times (-p\times y)$，这就是 $a{x}^{\prime }+b{y}^{\prime }$​ 的形式。

因为 $0\le a\mathrm{mod}r<r$，且 $r$ 为 $a\times x+b\times y$ 的最小正整数值，

所以 $a\mathrm{mod}r=0$，即 $r$ 是 $a$ 的因数。

同理可得 $r$ 也是 $b$ 的因数。

所以 $r$ 是 $a,b$ 的公因数，即 $r\mid \gcd (a,b)$。

又因为 $a\times x+b\times y$ 都是 $\gcd (a,b)$ 的倍数，所以 $\gcd (a,b)\mid r$。

故 $r=\gcd (a,b)$。

$\text{Exgcd}$

首先考虑欧几里得算法是如何求解 $\gcd (a,b)$ 的。

1. 当 $b=0$ 时，$a$ 即为 $a,b$ 的最大公约数。
2. 求出 $\gcd (b,a\mathrm{mod}b)$，设为 $g_0$。
3. 令 $g=g_0$，$g$ 就是 $a,b$​ 的最大公约数。

---

令 $d=\gcd (a,b)$ 于是可以仿照这个过程，求解不定方程：

1. 当 $b=0$ 时，$ax+by=d$ 的解是 $x=\frac{d}{a}$，$y$ 任意。

2. 求出 $bx+(a\mathrm{mod}b)y=d$ 的一组解，设为 $(x_0,y_0)$。

3. 令 $(x,y)=f((x_0,y_0))$，那么 $(x,y)$ 就是 $ax+by=d$ 的一组解。

   其中 $f((x_0,y_0))$ 表示对 $(x_0,y_0)$ 的一种变换。

现在的问题就是，如何通过 $(x_0,y_0)$ 得到 $(x,y)$​。

PS在欧几里得算法中，$\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$，故下面的式子中 $d$ 是相等的。

已知 $b{x}_0+(a\mathrm{mod}b)y_0=d$，转化形式可得：
$$b{x}_0+(a-\lfloor \frac{a}{b}\rfloor b)y_0=d$$
将 $a,b$ 提出来，可得
$$a{y}_0+b(x_0-\lfloor \frac{a}{b}\rfloor y_0)=d=ax+by$$
通过待定系数法，可得
$$\{\begin{array}{l}x=y_0\\ y=x_0-\lfloor \frac{a}{b}\rfloor y_0\end{array}$$

---

具体实现中，在数学上，当最后 $b=0$ 时 $y$ 可以取任意值，但是如果 $y$ 取得很大的话，中途中可能会超过 int 和 long long 的存储范围。

但是如果最后 $b=0$ 时 $y=0$，那么对于其它的迭代步骤，求解 $ax+by=\gcd (a,b)$ 得到的一组解 $(x,y)$ 一定满足：
$$\{\begin{array}{l}|x|\le \left|\frac{b}{2\times \gcd (a,b)}\right|\\ |y|\le \left|\frac{a}{2\times \gcd (a,b)}\right|\end{array}$$

证明：

PS：这里只考虑 $a>b$ 的情况，显然 $a\le b$ 在依次迭代后会变成 $a>b$ 的情况。

首先，当 $a\mathrm{mod}b=0$ 时，因为 $a>b$ 且 $a$ 是 $b$ 的倍数，所以 $a\ge 2\times b$。迭代得到的解
$$\{\begin{array}{l}x=y_0=0<\left|\frac{b}{2\times \gcd (a,b)}\right|\\ y=x_0-\lfloor \frac{a}{b}\rfloor y_0=x_0=\left|\frac{2\times b}{2\times b}\right|\le \left|\frac{a}{2\times \gcd (a,b)}\right|\end{array}$$
其次，设 $bx+(a\mathrm{mod}b)y=\gcd (b,a\mathrm{mod}b)$ 得到的解为 $(x_0,y_0)$，满足：
$$\{\begin{array}{l}{x}_0\le \left|\frac{a\mathrm{mod}b}{2\times \gcd (b,a\mathrm{mod}b)}\right|\\ y_0\le \left|\frac{b}{2\times \gcd (b,a\mathrm{mod}b)}\right|\end{array}$$
则 $ax+by=\gcd (a,b)$ 得到的解 $(x,y)$ 也满足：

1. $|x|=|y_0|\le \left|\frac{b}{2\times \gcd (b,a\mathrm{mod}b)}\right|=\left|\frac{b}{2\times \gcd (a,b)}\right|$；

2. 由于 $|x_0|\le \left|\frac{a\mathrm{mod}b}{2\times \gcd (b,a\mathrm{mod}b)}\right|$，$|\lfloor \frac{a}{b}\rfloor y_0|\le \left|\frac{\lfloor \frac{a}{b}\rfloor }{2\times \gcd (a,b)}\right|$，

   所以
   $$|y|=|x_0-\lfloor \frac{a}{b}\rfloor y_0|\le |x_0|+|\lfloor \frac{a}{b}\rfloor y_0|=\left|\frac{\lfloor \frac{a}{b}\rfloor b+a\mathrm{mod}b}{2\times \gcd (a,b)}\right|=\left|\frac{a}{2\times \gcd (a,b)}\right|$$

又可以发现，设 $ax+by=d=k\times \gcd (a,b)$，那么每次得到的解都是在求解 $ax+by=\gcd (a,b)$ 时的解的 $k$ 倍，所以不管 $d$ 取多少，在最后的迭代中取 $x=k,y=0$ 同样是没有问题的。

---

求出一对解 $x_0,y_0$ 后，可求出通解为：
$$\{\begin{array}{l}x=x_0+\frac{b}{\gcd (a,b)}\times t\\ y=y_0-\frac{a}{\gcd (a,b)}\times t\end{array}t\in \mathbb{Z}$$

---

时间复杂度与欧几里得算法相同，故时间复杂度为 O ( log ⁡ max ⁡ { a , b } ) O(\log\max\{a,b\}) O(logmax{a,b})。

$\text{Code}$

using lint = long long;

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