用反证法证明辗转相除法

辗转相除法是一种快速求两个数最大公约数的一种算法，其核心的思想是 gcd(a,b) = gcd(b, a mod b)。至于为什么这样，本文就给出证明过程。

辗转相除法定义

设 a, b 为正整数，不妨令 a >= b，gcd(m, n)是求m, n的最大公约数，那么以下算法：

1. 令 c = a mod b，若 c = 0，则b是最大公约数；
2. 若 c ≠ 0，求 gcd(b, a mod b)，回到第1步

流程图如下所示：

注1*：此操作的含义是将b和c的值作为新的a和b的值再代入gcd(a, b)。举例来说，假设输入a = 27, b = 18，则所要求的为 gcd(27,18)。先计算 c = 27 mod 18 = 9，此时 c ≠ 0，那么开始第二次循环，此时 a = b = 18, b = c = 9，所以 gcd(27, 18) 变为 gcd(18, 9)。

分析

辗转相除法的核心要求是证明 gcd(a, b) = gcd(b, a mod b) ，其中 $a,b$ 是正整数，函数 gcd 用于求 a, b 的最大公约数。设正整数 g = gcd(a, b)，那么就是要证明gcd(b, a mod b)=g。

证明

设 r = a mod b，则 a = k * a + r，其中 k 为整数。设 g = gcd(a, b)，那么：
① b mod g = 0，乘以 k 后可得：(k1 * b) mod g = 0；
② a mod g = 0，即 (k * b + r) mod g = 0；
联合①，②式可得 r mod g = 0，也就是说 r 也是 g 的整数倍，
所以 a，b，r 可以表示成 k1 * g, k2 * g 和 k3 * g，其中 k1 >= k2 > k3
因此，gcd(b, a mod b) = gcd(b, r) = gcd(k2 * g, k3 * g)，所以 g 是其公约数。

但是g是不是最大公约数还要证明，这里使用反证法，即假设 g 不是最大的公约数，那么存在 G = gcd(b, a mod b) G > g。G 是 b 和 a mod b 的公约数，而 a = k * b + (a mod b)，所以 G 也是 a 的公约数，即 G = gcd(a, b)，这与 g = gcd(a, b) 矛盾，所以 g 是 gcd(b, a mod b) 的最大公约数。

证毕。

小结

本文对辗转相除法利用反证法进行了证明，同时笔者还发现，还有一种通过两边范围相夹的方法来证明，参见：辗转相除法求得最大公约数的证明。基本证明方法见附录的总结。

附：另一种证明方法

设 g1 = gcd(a, b) 和 g2 = gcd(b, a mod b)，其中 a, b, c 均为正整数，那么：

• g1 $\le$ g2
  因为g1是a, b的公约数，所以 g1 也是 a mod b 的约数（参见本文前面的结论），所以 g1 是 b, c 的公约数，而 g2 是b, a mod b 的最大公约数，所以 g1 $\le$ g2。
• g2 $\le$ g1
  证明方法类似，可以先证明 g2 是 a, b 的公约数，而 g1 是a, b 的最大公约数，所以g2 $\le$ g1。

基于1和2可得 g1 = g2，即 gcd(a, b) = gcd(b, a mod b)，请毕。