扩展欧几里得算法的由来和C++实现

1、简述因子的概念

假设存在两个整数 m、n，若 n 除以 m 的结果是整数，即没有余数，则表示 m 是 n 的因子，记作 m | n。

2、欧几里得辗转相除法的证明

证：假设存在 a > b > 0，且 a , b∈Z，有 $\mathbf{a}=\mathbf{q}\cdot \mathbf{b}+\mathbf{r}$，即 $\mathbf{q}=[\mathbf{a}/\mathbf{b}]$，$0\le \mathbf{r}<\mathbf{b}$

  令 $\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{x},\mathbf{y})$ 表示 x 和 y 的最大公约数

  ① 当 r = 0 时，$\mathbf{a}=\mathbf{q}\cdot \mathbf{b}$，即 $\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})=\mathbf{b}$

  ② 当 r ≠ 0 时，假设存在一个数 d 使得 d|b 且 d|b，即有：

     $\mathbf{d}|\mathbf{a}\wedge \mathbf{d}|\mathbf{b}⟶\mathbf{d}|\mathbf{b}\wedge \mathbf{d}|(\mathbf{a}-\mathbf{q}\cdot \mathbf{b})⟶\mathbf{d}|\mathbf{b}\wedge \mathbf{d}|\mathbf{r}$
     $\mathbf{d}|\mathbf{b}\wedge \mathbf{d}|\mathbf{r}⟶\mathbf{d}|(\mathbf{q}\cdot \mathbf{b}+\mathbf{r})\wedge \mathbf{d}|\mathbf{b}⟶\mathbf{d}|\mathbf{a}\wedge \mathbf{d}|\mathbf{b}$

  于是得到：$\mathbf{d}|\mathbf{a}\wedge \mathbf{d}|\mathbf{b}⟷\mathbf{d}|\mathbf{b}\wedge \mathbf{d}|\mathbf{r}$
  令   d 1 = {   d   ∣   d ∣ a ∧ d ∣ b   } 、 d 2 = {   d   ∣   d ∣ b ∧ d ∣ r   } \ \bm{d_1=\{\ d\ \bigg| \ d|a\land d|b \ \}}、 \bm{d_2=\{\ d\ \bigg| \ d|b\land d|r \ \}}  d1​={ d ∣∣∣∣​ d∣a∧d∣b }、d2​={ d ∣∣∣∣​ d∣b∧d∣r }
  即 $\mathbf{d}|\mathbf{a}\wedge \mathbf{d}|\mathbf{b}⟷\mathbf{d}|\mathbf{b}\wedge \mathbf{d}|\mathbf{r}⟹{\mathbf{d}}_1={\mathbf{d}}_2⟹{\mathbf{d}}_{1\mathbf{m}\mathbf{a}\mathbf{x}}={\mathbf{d}}_{2\mathbf{m}\mathbf{a}\mathbf{x}}$

  其中${\mathbf{d}}_1、{\mathbf{d}}_2$是集合，${\mathbf{d}}_{1\mathbf{m}\mathbf{a}\mathbf{x}}、{\mathbf{d}}_{2\mathbf{m}\mathbf{a}\mathbf{x}}$是具体的值

  $\therefore \mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})=\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{b},\mathbf{r})=\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{b},\mathbf{a}\mathbf{m}\mathbf{o}\mathbf{d}\mathbf{b})$

3、最小公倍数和最大公约数的关系

两个数的乘积等于其最大公约数和最小公倍数的乘积

证：假设存在 a > b > 0 ，且 a , b∈Z

  令 $\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{x},\mathbf{y})$ 表示 x 和 y 的最大公约数，$\mathbf{l}\mathbf{c}\mathbf{m}(\mathbf{x},\mathbf{y})$ 表示 x 和 y 的最小公倍数

  又假设 $\mathbf{t}=\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})$，所以 ${\mathbf{a}}_0=\mathbf{a}/\mathbf{t}$ 和 ${\mathbf{b}}_0=\mathbf{b}/\mathbf{t}$互质

  $\therefore \mathbf{g}\mathbf{c}\mathbf{d}({\mathbf{a}}_0,{\mathbf{b}}_0)=1$，$\mathbf{l}\mathbf{c}\mathbf{m}({\mathbf{a}}_0,{\mathbf{b}}_0)={\mathbf{a}}_0\cdot {\mathbf{b}}_0$

  $\because \mathbf{l}\mathbf{c}\mathbf{m}(\mathbf{a},\mathbf{b})=\mathbf{l}\mathbf{c}\mathbf{m}({\mathbf{a}}_0\cdot \mathbf{t},{\mathbf{b}}_0\cdot \mathbf{t})=\mathbf{t}\cdot \mathbf{l}\mathbf{c}\mathbf{m}({\mathbf{a}}_0,{\mathbf{b}}_0)=\mathbf{t}\cdot {\mathbf{a}}_0\cdot {\mathbf{b}}_0$

  $\therefore \mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})\cdot \mathbf{l}\mathbf{c}\mathbf{m}(\mathbf{a},\mathbf{b})={\mathbf{t}}^2\cdot \mathbf{a}\cdot \mathbf{b}=(\mathbf{t}\cdot \mathbf{a})(\mathbf{t}\cdot \mathbf{b})={\mathbf{a}}_0\cdot {\mathbf{b}}_0$

  故两个数的乘积等于其最大公约数和最小公倍数的乘积得证

4、贝祖定理

关于 x , y（x，y 都是整数）的不定方程，一定存在一组解 （x , y） 使得 $\mathbf{a}\mathbf{x}+\mathbf{b}\mathbf{y}=\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})$ 成立，其中 a , b 都是整数，$\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})$表示 a 和 b 的最大公约数。

证：假设存在 a > b > 0 ，且 a , b∈Z，令$\mathbf{d}=\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})$

  $\because \mathbf{a}=\mathbf{b}\cdot \mathbf{q}+\mathbf{r}(0\le \mathbf{r}<\mathbf{b})$，且 $\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})=\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{b},\mathbf{r})$

  $\therefore \mathbf{b}=\mathbf{r}\cdot {\mathbf{q}}_1+{\mathbf{r}}_1(0\le {\mathbf{r}}_1<\mathbf{r})$

  又 $\because \mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{b},\mathbf{r})=\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{r},{\mathbf{r}}_1)$

   $\therefore \mathbf{r}={\mathbf{r}}_1\cdot {\mathbf{q}}_2+{\mathbf{r}}_2(0\le {\mathbf{r}}_2<{\mathbf{r}}_1)$

  以此类推：$\dots \dots$

  ${\mathbf{r}}_{\mathbf{n}-3}={\mathbf{r}}_{\mathbf{n}-2}\cdot {\mathbf{q}}_{\mathbf{n}-1}+{\mathbf{r}}_{\mathbf{n}-1}(0\le {\mathbf{r}}_{\mathbf{n}-1}<{\mathbf{r}}_{\mathbf{n}-2})$

  ${\mathbf{r}}_{\mathbf{n}-2}={\mathbf{r}}_{\mathbf{n}-1}\cdot {\mathbf{q}}_{\mathbf{n}}+{\mathbf{r}}_{\mathbf{n}}(0\le {\mathbf{r}}_{\mathbf{n}}<{\mathbf{r}}_{\mathbf{n}-1})$

  ${\mathbf{r}}_{\mathbf{n}-1}={\mathbf{r}}_{\mathbf{n}}\cdot {\mathbf{q}}_{\mathbf{n}+1}({\mathbf{r}}_{\mathbf{n}+1}=0)$ ， 即 $\mathbf{d}={\mathbf{r}}_{\mathbf{n}}$

  此时 $\mathbf{b}>\mathbf{r}>{\mathbf{r}}_1>{\mathbf{r}}_2>\cdot \cdot \cdot >{\mathbf{r}}_{\mathbf{n}-1}>{\mathbf{r}}_{\mathbf{n}}>{\mathbf{r}}_{\mathbf{n}+1}=0$

  $\therefore \mathbf{d}={\mathbf{r}}_{\mathbf{n}}={\mathbf{r}}_{\mathbf{n}-2}-{\mathbf{q}}_{\mathbf{n}}\cdot {\mathbf{r}}_{\mathbf{n}-1}$
     $={\mathbf{r}}_{\mathbf{n}-2}-{\mathbf{q}}_{\mathbf{n}}\cdot ({\mathbf{r}}_{\mathbf{n}-3}-{\mathbf{q}}_{\mathbf{n}-1}\cdot {\mathbf{r}}_{\mathbf{n}-2})$
     $=-{\mathbf{q}}_{\mathbf{n}}\cdot {\mathbf{r}}_{\mathbf{n}-3}+(1+{\mathbf{q}}_{\mathbf{n}}\cdot {\mathbf{q}}_{\mathbf{n}-1})\cdot {\mathbf{r}}_{\mathbf{n}-2}$
     $\cdots \cdots$
     $=\mathbf{x}\cdot \mathbf{a}+\mathbf{y}\cdot \mathbf{b}$

  由于$\mathbf{q},{\mathbf{q}}_1,{\mathbf{q}}_2,...,{\mathbf{q}}_{\mathbf{n}-1},{\mathbf{q}}_{\mathbf{n}}$都是整数

  所以 x 和 y 必是整数

  故一定存在一组解 （x , y） 使得 $\mathbf{a}\mathbf{x}+\mathbf{b}\mathbf{y}=\mathbf{g}\mathbf{c}\mathbf{d}(\mathbf{a},\mathbf{b})$ 成立

5、通过扩展欧几里得算法求解贝祖等式

解：令 ${\mathbf{a}}_0{\mathbf{x}}_0+{\mathbf{b}}_0{\mathbf{y}}_0=\mathbf{g}\mathbf{c}\mathbf{d}({\mathbf{a}}_0,{\mathbf{b}}_0)---(1)$

  由欧几里得辗转相除公式可得：

  $\mathbf{g}\mathbf{c}\mathbf{d}({\mathbf{b}}_0,{\mathbf{a}}_0\mathbf{m}\mathbf{o}\mathbf{d}{\mathbf{b}}_0)=\mathbf{g}\mathbf{c}\mathbf{d}({\mathbf{a}}_0,{\mathbf{b}}_0)---(2)$

  $\because \mathbf{x}\mathbf{m}\mathbf{o}\mathbf{d}\mathbf{y}=\mathbf{x}-\lfloor \mathbf{x}/\mathbf{y}\rfloor \cdot \mathbf{y}$ ，其中 $\lfloor \rfloor$ 为向下取整符号

  $\therefore \mathbf{g}\mathbf{c}\mathbf{d}({\mathbf{b}}_0,{\mathbf{a}}_0\mathbf{m}\mathbf{o}\mathbf{d}{\mathbf{b}}_0)={\mathbf{b}}_0{\mathbf{x}}_1+({\mathbf{a}}_0-\lfloor {\mathbf{a}}_0/{\mathbf{b}}_0\rfloor \cdot {\mathbf{b}}_0)\cdot {\mathbf{y}}_1---(3)$

  将 (1) (3) 式代入 (2) 可得：

  ${\mathbf{a}}_0{\mathbf{x}}_0+{\mathbf{b}}_0{\mathbf{y}}_0={\mathbf{b}}_0{\mathbf{x}}_1+({\mathbf{a}}_0-\lfloor {\mathbf{a}}_0/{\mathbf{b}}_0\rfloor \cdot {\mathbf{b}}_0)\cdot {\mathbf{y}}_1---(4)$

  整理 (4) 得：${\mathbf{a}}_0{\mathbf{x}}_0+{\mathbf{b}}_0{\mathbf{y}}_0={\mathbf{a}}_0{\mathbf{y}}_1+{\mathbf{b}}_0({\mathbf{x}}_1-\lfloor {\mathbf{a}}_0/{\mathbf{b}}_0\rfloor \cdot {\mathbf{y}}_1)---(5)$

  $\therefore ({\mathbf{x}}_0,{\mathbf{y}}_0)=({\mathbf{y}}_1,{\mathbf{x}}_1-\lfloor {\mathbf{a}}_0/{\mathbf{b}}_0\rfloor \cdot {\mathbf{y}}_1)$

  为了进一步求得 x1 和 y1 的值，可以按上述步骤求解

  故$({\mathbf{x}}_1,{\mathbf{y}}_1)=({\mathbf{y}}_2,{\mathbf{x}}_2-\lfloor {\mathbf{a}}_1/{\mathbf{b}}_1\rfloor \cdot {\mathbf{y}}_2)$

  其中 $\mathbf{g}\mathbf{c}\mathbf{d}({\mathbf{a}}_1,{\mathbf{b}}_1)=\mathbf{g}\mathbf{c}\mathbf{d}({\mathbf{b}}_0,{\mathbf{a}}_0\mathbf{m}\mathbf{o}\mathbf{d}{\mathbf{b}}_0)$，即${\mathbf{a}}_1={\mathbf{b}}_0,{\mathbf{b}}_1={\mathbf{a}}_0\mathbf{m}\mathbf{o}\mathbf{d}{\mathbf{b}}_0$

  $\cdots \cdots$

  以此类推，直到递归至第 n 次使${\mathbf{b}}_{\mathbf{n}}=0,{\mathbf{a}}_{\mathbf{n}}=\mathbf{g}\mathbf{c}\mathbf{d}({\mathbf{a}}_{\mathbf{n}},{\mathbf{b}}_{\mathbf{n}})$，即$({\mathbf{x}}_{\mathbf{n}},{\mathbf{y}}_{\mathbf{n}})=(1,0)$

  然后将该组解逐一回代能够求得：$({\mathbf{x}}_{\mathbf{n}-1},{\mathbf{y}}_{\mathbf{n}-1}),({\mathbf{x}}_{\mathbf{n}-2},{\mathbf{y}}_{\mathbf{n}-2}),\dots ,({\mathbf{x}}_0,{\mathbf{y}}_0)$

  由于原方程为不定方程 ，于是求得的$({\mathbf{x}}_0,{\mathbf{y}}_0)$只是一组特解

  通解为：$(\mathbf{x},\mathbf{y})=({\mathbf{x}}_0\pm \mathbf{k}{\mathbf{b}}_0,{\mathbf{y}}_0\mp \mathbf{k}{\mathbf{a}}_0)$，$\mathbf{k}\in \mathbf{Z}$

  $(\mathbf{a}\mathbf{x}+\mathbf{b}\mathbf{y}=\mathbf{c}⟹\mathbf{a}(\mathbf{x}+\mathbf{k}\mathbf{b})+\mathbf{b}(\mathbf{y}-\mathbf{k}\mathbf{a})=\mathbf{a}\mathbf{x}+\mathbf{b}\mathbf{y}=\mathbf{c})$

C++算法实现

void extend_gcd(int a, int b,int &x, int &y){
	if(b==0) x=1,y=0;
	else {
		extend_gcd(b,a%b,x,y);
		int temp = y;
		y = x- a/b*y;
		x = temp;
	}
	return;
}