扩展欧几里得算法 推导 + 代码

1. 目的
   求解ax + by = c 的整数解

2. 推导过程
     先求解方程$a\ast x+b\ast y=gcd(a,b)$的一组特解 $$\begin{gathered} a\ast x_1+b\ast y_1=gcd(a,b)(1) \\ b\ast x_2+(a\%b)\ast y_2=gcd(b,a\%b)(2) \end{gathered}$$  我们对公式（2）根据$a\%b=a-\lfloor \frac{a}{b}\rfloor \ast b$做展开得到$$a\ast y_2+b\ast (x_2-\lfloor \frac{a}{b}\rfloor \ast y_2)=gcd(b,a\%b)(3)$$  根据欧几里得算法$gcd(a,b)==gcd(b,a\%b)$对方程（1），（2）进行合并可得$$x_1=y_2,y_1=x_2-\lfloor \frac{a}{b}\rfloor \ast y_2(4)$$  如此反复进入递归k - 1次，边界为$b_k=0,gcd(a_k,0)=a_k$边界方程为$a_k\ast x_k=a_k$得到一特解$$x_k=1,y_k=0$$  根据方程（4）我们可以在回溯过程中求出一组特解$x_1,y_1$

//这里为什么要加&取值符，因为在递归调用中，父函数依赖于子函数的解（x，y）
//但是子函数不方便直接return给父函数解（当然也可以借助结构体，较为麻烦）
//所以直接借助&取值符，这样子函数中的x1，y1本质就是父函数中的x2，y2（因为传入的是地址，子函数与父函数公用，不产生新的复制参数）
//这样子函数x1，y1有解了，父函数中x2，y2也有值了（本来就是同一对变量）
int exgcd(int a,int b,int &x1,int &y1)
{
	if(!b)
	{
		x1=1,y1=0;
		return a;
	}
	int x2,y2;
	int d = exgcd(b,a%b,x2,y2);
	x1=y2,y1=x2-(a/b)*y2;
	return d;
}

3. 求通解
     方程$a\ast x+b\ast y=gcd(a,b)$中,我们对特解$x$加上常数$k_1$，对特解$y_1$加上常数$k_2$，若果$k_1$，$k_2$造成的影响相互抵消，那么$x+k_1，y+k_2$就是一组新的解
     我们对方程$(x+k_1)\ast a+(y+k_2)\ast b=gcd(a,b)$做展开
   $$a\ast x+b\ast y+k_1\ast a+k_2\ast b=gcd(a,b)(5)$$  易知只要满足$k_1\ast a=k_2\ast b=t\ast lcm(a,b)$  t为整数， 方程（5）就会成立
   $$k_1=t\ast lcm(a,b)/a=t\ast b/gcd(a,b)$$ $$k_2=t\ast lcm(a,b)/a=t\ast a/gcd(a,b)$$  即可得到任意整数解$(x+k_1,y+k_2)$
     回到方程ax + by = c，有整数解的前提为$gcd(a,b)\ast q==c$,即为$c\%gcd(a,b)==0$
   解集为$a\ast x+b\ast y=gcd(a,b)$的解集 $\ast c/gcd(a,b)$