欧几里得算法&扩展欧几里得算法&线性同余方程学习笔记

一、最大公约数

一组整数的公约数，是指能够同时整除这组数中每一个数的整数。
一组整数的最大公约数（Greatest Common Divisor，GCD），是指这组数所有公约数中最大的一个。
例如， $6,12,18,30,52$ 的最大公约数是 $2$ 。

二、欧几里得算法

下面的式子是欧几里得算法的核心：
$$\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$$

接下来给出证明：
① $a,b$ 的公约数必然为 $b,a\mathrm{mod}b$ 的公约数。
设 $d$ 为 $a,b$ 的公约数，$a=kb+c$，则 $d\mid a,d\mid b$，则 $a\mathrm{mod}b=c$，由于 $d\mid a,d\mid (kb)$，故 $d\mid c$ ，因此 $d\mid (a\mathrm{mod}b)$，又因为 $d\mid b$，所以 $a,b$ 的公约数必然为 $b,a\mathrm{mod}b$ 的公约数。
② $b,a\mathrm{mod}b$ 的公约数必然为 $a,b$ 的公约数。
设 $d$ 为 $b,a\mathrm{mod}b$ 的公约数，$a=kb+c$ 则 $d\mid b,d\mid (a\mathrm{mod}b)$ ，则 $a\mathrm{mod}b=c$，由于 $d\mid (kb),d\mid c$，故 $d\mid a$，又因为 $d\mid b$，所以 $b,a\mathrm{mod}b$ 的公约数必然为 $a,b$ 的公约数。

由①②得，$a,b$ 的公约数和 $b,a\mathrm{mod}b$ 的公约数完全相同，因此，它们的最大公约数也相同，即 $\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$ 。

在计算两个数的最大公约数时，可以根据这个等式，不断递归，直到 $b=0$ 。代码如下。

int gcd(int a, int b){
	return b == 0 ? a : gcd(b, a % b);
}

递归过程中 $a\mathrm{mod}b$ 至少将 $a$ 的规模减半，所以递归层数为 $\log n$ 层，复杂度为 $O(\log n)$ 。

三、多个数的最大公约数

计算多个整数的最大公约数时可以根据：
$$\gcd (a_1,a_2,a_3,...,a_n)=\gcd (gcd(a_1,a_2),a_3,a_4,...,a_n)$$

不断减少数字的个数，最终简化成两个整数的情况。

四、扩展欧几里得算法

扩展欧几里得算法常用于求解类似 $ax+by=gcd(a,b)$ 形式的二元一次方程组的一组可行解。
顾名思义，扩展欧几里得算法与欧几里得算法有相似之处，其实，这两种算法都使用了 $\gcd (a,b)=\gcd (b,a\mathrm{mod}b)$ 这条性质不断缩小问题规模，而后回溯求解。因此，不妨构造一个与 $ax+by=\gcd (a,b)$ （以下称作方程①） 类似的方程 $bx+(a\mathrm{mod}b)y=\gcd (b,a\mathrm{mod}b)$ （以下称作方程②） 并且假设已知后面方程组的一组解为 $x=x_0,y=y_0$。因此有：
$$ax+by=\gcd (a,b)=\gcd (b,a\mathrm{mod}b)=b{x}_0+(a\mathrm{mod}b)y_0$$

$$ax+by=b{x}_0+(a-\lfloor \frac{a}{b}\rfloor \times b)y_0$$

$$a(x-y_0)+b(y-x_0+\lfloor \frac{a}{b}\rfloor y_0)=0$$

故可以得到原方程的解为 $x=y_0,y=x_0-\lfloor \frac{a}{b}\rfloor y_0$。这样，在得知方程②解的情况下可以推得方程①的解。那么如何求出方程②的解呢？
通过观察两个方程的系数可以发现，从① $\to$ ②的过程中，系数会通过取模减小的，那么同欧几里得算法一样，不断递归，使系数减小，直到 $b=0$ ，此时根据欧几里得算法有 $a=\gcd (a,b)$，求出此时方程 $\gcd (a,b)x+0y=\gcd (a,b)$ 的解为 $x=1,y$ 为任意值，不妨设 $y=0$，不断回溯得到原方程 $ax+by=\gcd (a,b)$ 的解。代码如下：

// 其中 d 为最大公约数，可以不写，将 exgcd 函数定义为 void 类型
int exgcd(int a, int b, int &x, int &y){
	if (!b){
		x = 1, y = 0;
		return a;
	}
	int d = exgcd(b, a % b, x, y);
	int tmp = x;
	x = y;
	y = tmp - (a / b) * y;
	return d;
}

递归过程同欧几里得算法，因此复杂度也为 $O(\log n)$ 。

五、扩展欧几里得算法求出的特解性质

使用扩展欧几里得算法求出的特解 $x=x_0,y=y_0$ 满足性质 $|x_0|<b,|y_0|<a$ 。采用归纳法证明如下：
不妨假设共进行 $n$ 层递归，原方程为 $ax+by=\gcd (a,b)$，第 $i$ 层递归得到的方程为 $a_{i}x+b_{i}y=\gcd (a,b)$，该层递归得到的解为 $x=x_i,y=y_i$ ，递归出口 $a_{n}x+b_{n}y=\gcd (a,b)$，其中 $b_n=0$。
① 当 $\gcd (a,b)=b_i$ 时：
有 $a\mathrm{mod}b=0$，因此递归层数 $i=n-1$ ，得到解为 $x=x_{n-1}=1,y=y_{n-1}=0$，此时满足$|x_{n-1}|<|b_{n-1}|,|y_{n-1}|<|a_{n-1}|$ 。
② 当 $\gcd (a,b)\ne b_i$ 时：
不妨假设第 $i+1$ 层递归得到的解 $x=x_{i+1},y=y_{i+1}$ 满足 $|x_{i+1}|<|b_{i+1}|,|y_{i+1}|<|a_{i+1}|$ 。根据系数递归的公式有 $a_{i+1}=b_i,b_{i+1}=a_i\mathrm{mod}{b}_i$ 。
则对于 $x_i$ ：
$$|x_i|=|y_{i+1}|<|a_{i+1}|=|b_i|$$

对于 $y_i$ ：
$$|y_i|=|x_{i+1}-\lfloor \frac{a_i}{b_i}\rfloor y_{i+1}|⩽|x_{i+1}|+|\lfloor \frac{a_i}{b_i}\rfloor y_{i+1}|=|x_{i+1}|+\lfloor \frac{a_i}{b_i}\rfloor |y_{i+1}|⩽|b_{i+1}|-\lfloor \frac{a_i}{b_i}\rfloor |y_{i+1}|$$

$$⩽a_i\mathrm{mod}{b}_i+\lfloor \frac{a_i}{b_i}\rfloor |y_{i+1}|⩽a_i-\lfloor \frac{a_i}{b_i}\rfloor b_i+\lfloor \frac{a_i}{b_i}\rfloor |y_{i+1}|⩽a_i-\lfloor \frac{a_i}{b_i}\rfloor (b_i-|y_{i+1}|)$$

又因为 $|y_{i+1}|<|a_{i+1}|=|b_i|$ ，所以 $b_i-|y_{i+1}|>0$，故 $|y_i|<|a_i|$ ，满足性质。
根据数学归纳法可得，原方程的解亦满足性质。

六、求解二元一次方程不定方程

接下来我们要求解二元一次方程 $ax+by=c$ 。
① $\gcd (a,b)\mid c$ 时，有无穷多组解。
使用扩展欧几里得算法求出方程 $ax+by=\gcd (a,b)$ 的特解 $x=x_0,y=y_0$ 后，方程 $ax+by=c$ 的一组特解可以表示为：
$$x=\frac{c}{\gcd (a,b)}{x}_0,y=\frac{c}{\gcd (a,b)}{y}_0$$

因此，$ax+by=c$ 的所有解可表示为：
$$x=\frac{c}{\gcd (a,b)}{x}_0+\frac{b}{\gcd (a,b)}k,y=\frac{c}{\gcd (a,b)}{y}_0-\frac{a}{\gcd (a,b)}k$$

其中 $k$ 为整数。

② $\gcd (a,b)\nmid c$ 时，无解。根据整除性质易证。

七、线性同余方程

形如 $ax\equiv c(\mathrm{mod}b)$ 的方程被称为线性同余方程。

八、求解线性同余方程

线性同余方程 $ax\equiv c(\mathrm{mod}b)$ 的解与 $ax+by=c$ 的解等价。同时，这两个方程解的存在性也相同。即当 $\gcd (a,b)\mid c$ 时，方程存在无数解，当 $\gcd (a,b)\nmid c$ 时，方程无解。

因此采用扩展欧几里得算法求出二元一次不定方程的解，那么这组解也是线性同余方程的解（此时无需关注 $y$ 的解是多少，只关注 $x$ 的解即可）。