深入浅出讲扩展欧几里得-exgcd

扩展欧几里得(exgcd)

扩展欧几里得的典型应用是解决形如

$$a\times x+b\times y=c$$

的二元一次方程的解的存在性问题以及求出特解与通解。

解的存在性问题

其实我们可以通过辗转相减法观察出来，$a,b$辗转相减的时候相当于$a\times x+b\times y$中的$x,y$不断变化的过程。我们通过前面的知识知道，这一个过程后，一定能得出最大公约数，所以$a\times x+b\times y=gcd(a,b)$一定有解，那么是否$gcd(x,y)\ne c$时就无解呢？

首先我们知道$c$是$gcd(a,b)$的倍数才可能有解，因为左边一定含有质因子$gcd(a,b)$，左边等于右边，则右边也一定含有$gcd(a,b)$，因此我们可以通过两边同时除以$gcd(a,b)$得到一个新的式子：

$$k_1\times x+k_2\times y=c/g,gcd(k1,k2)=1$$

由辗转相减可得：$k_1\times x+k_2\times y=1$一定有解。所以，$k_1\times x+k_2\times y=c/g$一定有解。

所以，当$c$是$gcd(a,b)$的倍数时，方程一定有解。

特解与通解

所以我们只要求形如：

$$a\times x+b\times y=gcd(a,b)$$

然后，将解扩大$c/gcd(a,b)$倍就可以得出方程的解了。

观察式子本质就是$a$和$b$辗转相减的过程。假设已经求出了

$$(a-b)\times x_1+b\times y_1=gcd(a,b)$$

的解$x_1,y_1$，那么变换一下就可以得到

$$a\times x_1+b\times (y_1-x_1)=gcd(a,b)$$

我们可以得出

$$x=x_1,y=y_1-x_1$$

得到了原方程的解，所以我们可以将原问题转化为一个规模更小的子问题。然后根据子问题的答案推出原问题的答案。考虑到这个辗转相减可以用辗转相除来替代，我们可以将求解原问题转化为求解

$$b\times x_1+a\%b\times y_1=gcd(a,b)$$

由：

$$a\%b=a-a/b\times b$$

代入得：

$$b\times x_1+(a-a/b\times b)\times y_1=gcd(a,b)$$

等价于：

$$a\times y_1+b\times (x_1-a/b\ast y_1)=gcd(a,b)$$

得到：

$$x=y_1,y=x_1-a/b\ast y_1$$

因此，我们只需要一直缩小问题的规模，直到变成$gcd(a,b)\times x+0\times y=gcd(a,b)$，然后得到一组解$(1,0)$，在这组解反推回去，得到一组原方程的特解。该过程可以用一个函数来实现。

代码：

$$exgcd$$

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

函数执行完毕后，代码里的$(x,y)$就是$a\times x+b\times y=gcd(a,b)$的一组特解。通解的变化规律就是$x,y$的值往相反的方向变化，比如$x$变大一点，$y$变小一点，使得答案不变，设这个变化的最小单位为$d_1,d_2$。

$$a\times (x+d_1)+b\times (y-d_2)=gcd(a,b)$$

$$a\times d_1=b\times d_2$$

两边同除$gcd(a,b)$令$k_1=a/gcd(a,b),k_2=b/gcd(a,b)$

$$k_1\times d_1=k_2\times d_2$$

此时，$gcd(k_1,k_2)=1$，显然$d_1=k_2,d_2=k_1$，可得通解为：

$$(x+k\times d_1,y-k\times a/gcd(a,b))$$

特解的绝对值大小

我们通过递归函数的实现可以观察到，解的绝对值是跟$a,b$的绝对值同一个级别的。（即传入与传出的数据不会相差太大，如传入的很小，传出的爆炸）。