欧几里得算法及拓展欧几里得算法

---

前言

本文主要介绍的是求$\gcd (a,b)$及方程$ax+by=\gcd (a,b)$的解的方法，即欧几里得算法和拓展欧几里得算法。

---

欧几里得算法

定理：$\forall a,b\in \mathbb{N},b\ne 0,\gcd (a,b)=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$。
证明：
当 $a<b$ ，则$\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)=\gcd (b,a)=\gcd (a,b)$，命题成立。
当$a\ge b$，设$a=qb+r$，其中$0\le r<b$。显然$r=a\mathrm{m}\mathrm{o}\mathrm{d}b$。而对于$a,b$的任意公约数$d|a,d|q\times b$，故$d|(a-qb)$，即$d|r$。因此$d$也是$b,r$的公约数。反之亦成立。故$a,b$的公约数集合与$b,a\mathrm{m}\mathrm{o}\mathrm{d}b$的公约数集合相同，他们的最大公约数自然也相等。
证毕。

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

---

拓展欧几里得算法

求方程$ax+by=\gcd (a,b)$的整数解。
设方程$ax+by=\gcd (a,b)$的一组解为$\{\begin{array}{l}x=x_1\\ y=y_1\end{array}$。由欧几里得算法得知，方程$ax+by=\gcd (a,b)$可转化为$bx+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y=\gcd (a,b)$，则有另一组解为$\{\begin{array}{l}x=x_2\\ y=y_2\end{array}$。此时把方程变形一下就可以变成：
$$\begin{gathered} a{x}_1+b{y}_1=b{x}_2+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_2 \\ a{x}_1+b{y}_1=b{x}_2+\left(a-\lfloor \frac{a}{b}\rfloor \times b\right)y_2 \\ a{x}_1+b{y}_1=a{y}_2+b{x}_2-\lfloor \frac{a}{b}\rfloor \times b{y}_2 \end{gathered}$$
根据恒等定理，有一组解为$\{\begin{array}{l}x=y_2\\ y=x_2-\lfloor \frac{a}{b}\rfloor \times y_2\end{array}$。可以发现，$x,y$基于$x_2,y_2$，于是在求最大公约数的时候，可以顺便求出。

int x,y;
int exgcd(int a,int b)
{
	if(b==0)
	{
		x=1,y=0;
		return a;
	}
	int temp=exgcd(b,a%b);
	int tx=x;
	x=y;
	y=tx-a/b*y;
	return temp;
}