扩展欧几里得定理 exgcd

扩展欧几里得求形如 $ax+by=c$ 的方程。

$$ax\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)⟺ax+my=b$$

• 对于不为 $0$ 的整数 $a,b$，存在整数 $x,y$，使得 $ax+by=gcd(a,b)$。

• 方程 $ax+by=c$ 有解，当且仅当 $gcd(a,b)|c$。

特殊情况：对于 $ax+by=gcd(a,b)$，当 $b=0$ 时，令 $x=1,y=0$ 即求得一组解。

一般情况：借鉴辗转相除法 $gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$，假设已经求得了 $(b\mathrm{m}\mathrm{o}\mathrm{d}a)x_0+a{y}_0=gcd(a,b)$ 的一组解 $(x_0,y_0)$，方程可以根据 $a\mathrm{m}\mathrm{o}\mathrm{d}m=a-m\lfloor \frac{a}{m}\rfloor$ 变形为

$$(b-a\lfloor \frac{b}{a}\rfloor )x_0+a{y}_0=gcd(a,b)$$

$$\lfloor \frac{b}{a}\rfloor x_0+a{y}_0+b{x}_0=gcd(a,b)$$

$$a(y_0-\lfloor \frac{b}{a}\rfloor x_0)+b{x}_0=gcd(a,b)$$

交换 $x_0$ 与 $y_0$
$$a(x_0-\lfloor \frac{b}{a}\rfloor y_0)+b{y}_0=gcd(a,b)$$

$$\{\begin{array}{l}ax+by=gcd(a,b)\\ a(x_0-\lfloor \frac{b}{a}\rfloor y_0)+b{y}_0=gcd(a,b)\end{array}⟺\{\begin{array}{l}x=x_0-y_0\lfloor \frac{b}{a}\rfloor \\ y=y_0\end{array}$$

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

如何求最小正整数解。假设我们求出了
$$a{x}_0+b{y}_0=gcd(a,b)$$
的一组解 $(x_0,y_0)$。有解的必要条件为 $gcd(a,b)|c$，所以令 $k=\frac{c}{gcd(a,b)}$。
$$ax+by=c$$

$$ax+by=k\times gcd(a,b)$$

$$\{\begin{array}{l}a\frac{x}{k}+b\frac{y}{k}=gcd(a,b)\\ a{x}_0+b{y}_0=gcd(a,b)\end{array}⟺\{\begin{array}{l}x=k\times x_0\\ y=k\times y_0\end{array}⟺\{\begin{array}{l}x=\frac{c}{gcd(a,b)}\times x_0\\ y=\frac{c}{gcd(a,b)}\times y_0\end{array}$$

$$a\times \frac{c}{gcd(a,b)}\times x_0+b\times \frac{c}{gcd(a,b)}\times y_0=c$$

为了求最小整数解，需要把上式的 $x_0$ 尽可能地减小，如果 $x_0$ 减小，对应的 $y_0$ 会增加。所以设 $x_0$ 减小了 $t_0$，$y_0$ 增加了 $t_1$。
$$\begin{gathered} \{\begin{array}{l}x=x_0\times \frac{c}{gcd(a,b)}-\frac{c}{gcd(a,b)}\times t_0\\ y=y_0\times \frac{c}{gcd(a,b)}+\frac{c}{gcd(a,b)}\times t_1\end{array} \\ \Downarrow \\ a\times x_0\times \frac{c}{gcd(a,b)}-a\times \frac{c}{gcd(a,b)}\times t_0+b\times y_0\times \frac{c}{gcd(a,b)}+b\times \frac{c}{gcd(a,b)}\times t_1=c \\ \because a\times \frac{c}{gcd(a,b)}\times x_0+b\times \frac{c}{gcd(a,b)}\times y_0=c \\ \therefore a\times \frac{c}{gcd(a,b)}\times t_0=b\times \frac{c}{gcd(a,b)}\times t_1 \\ \frac{t_0}{t_1}=\frac{b\times \frac{c}{gcd(a,b)}}{a\times \frac{c}{gcd(a,b)}} \\ \frac{t_0}{t_1}=\frac{b}{a} \end{gathered}$$

$$a(x_0-t_0)+b(y_0+t_1)=c⟺a(x_0-t_0)+b(y_0+\frac{a}{b}\times t_0)$$

让 $t_0=\frac{b}{gcd(a,b)}$，将 $x_0$ 减到不能再减为止，并考虑 $\frac{b}{gcd(a,b)}<0$ 的情况需要加上模数再取模。总结公式为
$$\{\begin{array}{l}{x}_{min}=(x_0\mathrm{m}\mathrm{o}\mathrm{d}\frac{b}{gcd(a,b)}+\frac{b}{gcd(a,b)})\mathrm{m}\mathrm{o}\mathrm{d}\frac{b}{gcd(a,b)}\\ y_{min}=\frac{c–a\times x_0}{b}\end{array}$$