扩展欧几里得

前言

划水选手来把以前学的算法写博客上，防止以后忘了
然后先定义$gcd(a,0)=a$

欧几里得算法

求$gcd(a,b)$即最大公约数，有什么用，能求$lcm(a,b)=a\ast b/gcd(a,b)$，）大雾
方法：$gcd(a,b)=gcd(b,a\%b)$递归即可
复杂度：$\mathcal{O}(logn)$

扩展欧几里得

求$ax+by=gcd(a,b)$的一组整数解
）当然保证一定有解
证明？因为我能求出来
方法：
我们发现，如果$b=0$，那么显然答案显然是$\{\begin{array}{l}x=1\\ y=0\end{array}$
既然是扩展欧几里得，那么我们就用欧几里得算法来
假设我们求出了$a^{\prime }=b,b^{\prime }=a\%b$时的解$a^{\prime }{x}^{\prime }+b^{\prime }{y}^{\prime }=gcd(a^{\prime },b^{\prime })=gcd(a,b)$
我们来化简一下$(a\%b)y^{\prime }+b{x}^{\prime }=gcd(a,b)$
我们发现现在离我们的目标$ax+by=gcd(a,b)$还有一点距离，就是第一项
我们向目标靠近
$$\begin{array}{rl}a{y}^{\prime }+b{x}^{\prime }& =gcd(a,b)+(a-a\%b)y^{\prime }\\ a{y}^{\prime }+b{x}^{\prime }& =gcd(a,b)+\lfloor \frac{a}{b}\rfloor b{y}^{\prime }\\ a{y}^{\prime }+b(x^{\prime }-\lfloor \frac{a}{b}\rfloor y^{\prime })& =gcd(a,b)\end{array}$$
这样我们只要递归就能找到解了
随便搞个代码

int EXgcd(const int a,const int b,int &x,int &y)
{
    if(!b)
	{
		x=1,y=0;
		return a;
	}
    const int res=EXgcd(b,a%b,y,x);
    y-=a/b*x;
    return res;
}

返回值是gcd(a,b)，然后x,y定义一下变量直接传进去（都不用预赋值），函数做完里面就是求的x,y
然后的话，至于通解
我们发现$a\ast b/gcd(a,b)=b\ast a/gcd(a,b)$
所以如果获取一组$\{\begin{array}{l}{x}^{\prime }\\ y^{\prime }\end{array}$
那么通解就是$\{\begin{array}{l}x=x^{\prime }+k\ast b/gcd(a,b)\\ y=y^{\prime }+k\ast a/gcd(a,b)\end{array},k\in Z$

额外的应用

然后如果要求$ax+by=c$的整数解咋办咧？
若$c\%gcd(a,b)\ne 0$无解
显然在x,y都是整数的时候$ax+by$一定是$gcd(a,b)$的整数倍
所以相当于变成求$ax+by=l\ast gcd(a,b)$的整数解
假设我们求出了$ax+by=gcd(a,b)$的解，貌似x、y各乘$l$就好了，通解也是原本的通解分别乘$l$

结语

这篇博客巨水，感兴趣的可以去看看裴蜀定理，那个看起来更牛逼bushi一点