拓展欧几里得算法

拓展欧几里得算法，可以用来求解一个一元二次方程的特解。
形如：$ax+by=gcd(a,b)$及其变式

为什么是拓展欧几里得呢，因为这个特解是从欧几里得算法(辗转相除法)的基础上来求得的，即$gcd->exgcd$。
我们知道求$gcd$时，最后$return$出来的结果是$b$，如果我们此时再往深递归一层，也就是$b=0$的时候，此时返回的结果值就是$gcd=a$，由此我们得到了在此情况下的一个特殊解$x=1,y=0$，由这个解我们就可以回溯一层层的求出其它的解了。

由$gcd(a,b)=gcd(b,a\%b)$，我们设当前层为$gcd(b,a\%b)$的话，设特解为$x_0和y_0$，其上一层的解为$x_1和y_1$，那么
$a{x}_1+b{y}_1=gcd(a,b)=gcd(b,a\%b)=b{x}_0+(a\%b)y_0$，由可知
$a\%b=a-\frac{a}{b}\ast b$，则
$a{x}_1+b{y}_1=b{x}_0+(a\%b)y_0=b{x}_0+(a-\frac{a}{b}\ast b)y_0=a{y}_0+b(x_0-\frac{a}{b}{y}_0)$
这样通过对应系数相等，我们就可以比较出来
$$\{\begin{array}{l}{x}_1=y_0\\ y_1=x_0-\frac{a}{b}\ast b\end{array}$$
这样我们就可以得到回溯得到我们需要的一组解了。

一
如果题目要我们求得是
$ax+by=c$这样的不定方程的解，如果$c\%gcd(a,b)\ne 0$，那么就无解，否则的话只需要把$c$先除到等式左边再乘上$gcd(a,b)$就变成了一般我们熟悉的形式了。
同理这样的我们求出来的$x和y$，也应改逆向操作，先除$gcd$再乘$c$，就是我们要求的解了。
二
如果求$ax+by=gcd(a,b)$的多组解，我们可以构造
$x=x^{\ast }+k\ast b$，$y=y^{\ast }-k\ast a$，形如这样的解带回原方程还是只剩下$ax+by=gcd(a,b)$这样的形式，注意这里有个细节就是防止$a$和$b$不互质这里要用$\frac{a}{gcd(a,b)}$和$\frac{b}{gcd(a,b)}$操作。
三
如果求得是最小的正整数解，求出来的$x$只需要先令$t=\frac{b}{gcd(a,b)}$（保证其为正，后面要取模），然后$x=(x\%t+t)\%t$即可。
证明的话，这篇大佬博客写的不错
四
还可以用扩展欧几里得算法来求解逆元,$ax\equiv 1(modp)$，即说明我们知道逆元存在条件即为$gcd(a,p)=1$，这样同余式其实就是可化为$ax+py=1$，这样逆元就可以用拓欧来直接求出来了。

未完待续~

void exgcd(int a,int b,int &x,int &y) {//拓欧求逆元
	if(b == 0) {//最后一层 满足 ax + by = gcd(a,b) 的特解 x = 1,y = 0
		x = 1; y = 0;
		return ;
	}
	exgcd(b,a%b,x,y);//上一层的特解x
	int t = x;  //x = y1,y = x1 - (a/b)*y1
	x = y;
	y = t - a / b * y;
}