扩展欧几里得

我在很久很久以前写过这篇文章但是因为自己特别菜当时不会$\text{LATEX}$，而且一些地方理解不够深入现在来做一下简单修改

扩展欧几里得算法($exgcd$)可以用来求下面这个关于$x,y$的不定方程的最小整数解
$$ax+by=\gcd (a,b)$$

我们回忆一下，我们之前使用辗转相除法求$\gcd$的过程
当我们辗转到$b_0=0$的时候，那么这个时候的$a_0=\gcd (a,b)$
那么由于满足
$$a_0{x}_0+b_0{y}_0=\gcd (a,b)$$
显然$x_0=1,y_0=0$
那么初值我们就可以确定了
接下来考虑回溯的时候如何更新$x,y$

根据迭代的方法，我们应该有一个等式
$$\begin{array}{rl}{a}_2{x}_2+b_2{y}_2& =a_1{x}_1+b_1{y}_1=\gcd (a,b)\\ a_2{x}_2+b_2{y}_2& =b_2{x}_1+(a_2\mathrm{m}\mathrm{o}\mathrm{d}{b}_2)y_1\\ a_2{x}_2+b_2{y}_2& =b_2{x}_1+(a_2-\lfloor \frac{a_2}{b_2}\rfloor b_2)y_1\\ a_2{x}_2+b_2{y}_2& =b_2{x}_1+a_2{y}_1-\lfloor \frac{a_2}{b_2}\rfloor b_2{y}_1\\ (x_2-y_1)a_2& +(y_2-x_1+\lfloor \frac{a_2}{b_2}\rfloor y_1)b_2=0\\ \because a_2,b_2& \ne 0\\ \therefore x_2-y_1=0,& y_2-x_1+\lfloor \frac{a_2}{b_2}\rfloor y_1=0\\ \therefore x_2=y_1,& y_2=x_1-\lfloor \frac{a_2}{b_2}\rfloor y_1\end{array}$$
然后就很显然了，我们在之前求$\gcd$的代码上简单修改就可以了好像还不少
然后注意我们最后求的是$x$的最小整数解，也就是说 $x$不保证是一个正数
要加回来
下面给出$exgcd$的核心代码

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