（4）线性方程与最大公因数

线性方程定理

已知两个整数a，b。我们观察a的倍数加上b的倍数得到的所有可能的整数。也就是$ax+by$得到的所有整数，其中x与y可以为任意整数，注意x，y可以取到负值。

例如我们列出a=42，b=30的所有整数，通过观察发现所有的数都是6的倍数，这很好理解，因为42与30都可以被6整除，$42x+30y=6(7x+5y)$的每个数都是6的倍数。更一般的，显然形如$ax+by$的每一个数都可以被$gcd(x,y)$整除，因为$gcd(x,y)$可以整除a，b。

第二个观察是42与30最大公因数6出现在了表中。

$42\ast (-2)+30\ast 3=6=gcd(42,30)$

我们这时可以假设一个断言：

形如 $ax+by$的最小正整数等于gcd(a,b)。

有许多方法可以证明这是成立的，我们可以利用欧几里得算法构造出合适的x，y，这里将计算过程列成表格，顺便与欧几里得算法过程进行对比。

$a=2\ast b+16$	$16=a-2b$
$b=1\ast 16+6$	$6=b-1\ast 16$
\,	$=b-1\ast (a-2b)$
\,	$=-a+3b$
$16=2\ast 6+4$	$4=16-2\ast 6$
\,	$=(a-2b)-2\ast (-a+3b)$
\,	$=3a-8b$
$6=1\ast 4+2$	$2=6-1\ast 4$
\,	$=(-a+3a)-1\ast (3a-8b)$
\,	$=-4a+11b$
$4=2\ast 2+0$	\,

按这种方法我们每次都会获得形如

最新余数=a的倍数+b的倍数

的等式。最终得到最后的非零余数，这就是$gcd(x,y)$，得解。

现在我们知道$ax+by=gcd(a,b)$总有整数解（x，y），我们可以进一步讨论方程有多少个解以及怎么表示所有的解。
从特殊到一般，我们首先从互素的a，b开始，即让$gcd(a，b)=1$，假设（x₁，y₁）是

$ax+by=1$

方程的一个解。通过x₁减去一个b的倍数和y₁加上一个a的相同的倍数，可得到其他解，换句话说，对于任何整数k，我们得到了新解$(x_1+kb,y_1-ka)$，通过计算

$a(x_1+kb)+b(y_1-ka)=a{x}_1+akb+b{y}_1-bka=a{x}_1+b{y}_1=1$

可验证这确实是解。

同样我们也可以通过解方程推导出答案，仍然观察$gcd(a,b)=1$的情况，假设(x₁,y₁)，（x₂,y₂）是方程的两个解，则有：

$a{x}_1+b{y}_1=1a{x}_2+b{y}_2=1$

我们用y₂乘第一个方程，用y₁乘第二个方程，相减就消去了b，重新整理后得到

$a{x}_1{y}_2=a{x}_2{y}_1=y_2-y_1$

类似地，如果用x₂乘以第一个方程，用x₁乘以第二个方程，再相减便得到

$b{x}_2{y}_1-b{x}_1{y}_2=x_2-x_1$

如果令$k=x_2{y}_1-x_1{y}_2$，则得到

$x_2=x_1+kb与y_2=y_1+ka$

这表明第二个解是由第一个解得到的，通过取不同的k可以得到每一个解$(x_1+kb,y_1-ka)$。

如果$gcd(a,b)>1$的情况呢？我们令$g=gcd(a,b)$由欧几里得一定可以得出

$ax+by=g$

至少有一个解（x₁，y₂）而g整除a，b，故（x₁，y₁）也是简单方程

$\frac{a}{g}x+\frac{b}{g}y=1$

的解，于是可以代入前面的方法得到解为$(x_1+k\ast \frac{b}{g},y_1-k\ast \frac{a}{g})$。

线性方程定理设a，b是非零整数，$g=gcd(a,b)$。方程

$ax+by=g$

总是有一个整数解（x₁，y₁），它可由前面的欧几里得算法得出，则方程的每一个解可由

$(x_1+k\ast \frac{b}{g},y_1-k\ast \frac{a}{g})$

得出，其中k为任意整数。

这节我们证明了方程$ax+by=g$总是有解，这个事实很重要，我们在接下来的数论中将多次用到。

算法代码

模拟上面的证明过程

int Eucild(int a1,int b1,int &x,int &y)
{
	int a=a1,b=b1;
	int r=a%b;
	x=1;
	int n=0,m=1,q;
	while(r)
	{		
		q=a/b;
		x=m-q*n;
		m=n;n=x;
		a=b;b=r;
		r=a%b;
	}
	y=(b-x*a1)/b1;
	return b;
}

但这并不是一个优秀的算法。（是我自己编的，太菜了）

接下来给出一个标准算法
对于一个方程$ax+by=gcd(a,b)$来说，我们可以做一下推导：
设有$a{x}_1+b{y}_1=gcd(a,b)$；
$b{x}_2+(a\%b)y_2=gcd(b,a\%b)$
对于这个方程组我们需要知道x₁，x₂，y₁，y₂的关系，才能递归的解决这个问题。
观察$b{x}_2+(a\%b)y_2=gcd(b,a\%b)$这个式子其中$(a\%b)$可以写为$a-\lfloor \frac{a}{b}\rfloor \ast b$。将括号打开合并a，b得

$a{y}_2+b(x_2+\lfloor \frac{a}{b}\rfloor y_2)$

由于欧几里得算法的原理$gcd(a,b)==gcd(b,a\%b)$,我们将两式子联立,对比系数即可得到$x_1=y_2,y_1=x_2-\lfloor ab\rfloor \ast y_2$
这个递归的边界是什么呢?我们知道,当朴素欧几里得到达边界时,$gcd(a,0)=a$,那么边界条件就是对$a\ast x_0+b\ast y_0=a$求解,很显然,此时x₀=1,y₀=0

int eucild(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;
}