整除性与最大公因数（欧几里得算法）

整除性与最大公因数

1 基本概念

  设$m,n$是整数，$m\ne 0$。$m$整除$n$指：$n$是$m$的倍数，即存在整数$k$使得$n=km$。整除$n$的数称为$n$的因数。

  对于两个已知整数，它们的公因数指：可以整除它们两个的数。

  两个整数$a,b$（不全为0）的最大公因数是整除它们两个的最大数，记为$gcd(a,b)$。如果$gcd(a,b)=1$，则称$a$与$b$互素。

  最大公因数是一个很重要的量，经常在数论中见到。

2 一些求最大公因数的方法

• 因数分解，例如求：
  $$gcd(225,120)=15$$
  通过因数分解$225=3^2\ast 5^2$，$120=2^3\ast 3\ast 5$可以验证答案是正确的。

• 打表，我们可以把$a,b$的所有公因数都找出来，然后找出同时出现在两个表中的最大数

3 欧几里得算法

  欧几里得算法是求两个数最大公因数的最有效方法，算法由直到余数为0的一系列带余除法组成的。

3.1 一个例子

  计算：
$$gcd(36,132)$$
第一步：
$$132=3\ast 36+24$$
第二步取$36$，用前一步的余数$24$除$36$得：
$$36=1\ast 24+12$$
下一步同理：
$$24=2\ast 12+0$$
欧几里得算法说明当得到余数0时，则前一步的余数就是最初两个数的最大公因数。因此$gcd(36,163)=12$。

3.2 形式化地表示

  欧几里得算法可以用下式来形式化地表示：
$$A=B\ast Q+R$$
每一步用$B,R$替换原来的$A,B$，重复此过程知道得到余数0为止。

  对于一般情形给出$a,b$有：
$$\begin{array}{rl}& a=q_1\ast b+r_1\\ & b=q_2\ast r_1+r_2\\ & r_1=q_3\ast r_2+r_3\\ & r_2=q_4\ast r_3+r_4\\ & \cdots \cdots \\ & r_{n-3}=q_{n-1}\ast r_{n-2}+r_{n-1}\\ & r_{n-2}=q_n\ast r_{n-1}+r_n\\ & r_{n-1}=q_{n+1}\ast r_n+0\end{array}$$
如果我们令$a=r_{-1},b=r_0$则上述的每一行都形如为：
$$r_{i-1}=q_{i+1}\ast r_i+r_{i+1}$$

$r_n$就是$a,b$的最大公因数

3.3 证明

  首先证明：$r_n$是$a,b$的公因数。自底向上分析。最后一行$r_{n-1}=q_{n+1}\ast r_n$说明$r_n$整除$r_{n-1}$，同理倒数第二行$r_{n-2}=q_n\ast r_{n-1}+r_n$表明$r_n$整除$r_{n-2}$，再往上依旧成立，直到第$1,2$行可以得到$r_n$整除$a$和$b$，故$r_n$是$a,b$的公因数得证。

  再证明：$r_n$是$a,b$的最大公因数。设$d$是$a,b$的任意公因数。自顶向下分析，我们从第一个等式$a=q_1\ast b+r_1$出发，由于$d$整除$a$和$b$，可知$d$整除$r_1$，到第二行式子$b=q_2\ast r_1+r_2$，由于$d$整除$b$和$r_1$，可知$d$整除$r_2$。再往下也是同理，对于每一步，$d$整除前两个余数$r_i$和$r_{i-1}$。最后我们能够得到$d$整除$r_n$的结论。既然$a,b$的任意公因数$d$都整除$r_n$，$r_n$也是$a,b$的公因数。那么$r_n$一定是$a$与$b$的最大公因数。

3.4 定理（欧几里得算法）

  要计算两个整数$a$与$b$的最大公因数，先令$r_{-1}=a$且$r_0=b$，然后计算相继的商和余数：
$$r_{i-1}=q_{i+1}\ast r_i+r_{i+1}(i=0,1,2,\cdots )$$
直到某余数$r_{n+1}$为$0$。最后的非零余数$r_n$就是$a,b$的最大公因数。

4 算法实现

  python实现$gcd$并计算$gcd(316258250,1160718174)$

def gcd(a, b):
    while a % b:
        a, b = b, a % b
    return b
    
print(gcd(316258250, 1160718174))

递归实现

def gcd(a, b):
    if not b:
        return a
    return gcd(b, a % b)

print(gcd(316258250, 1160718174))

结果：

1078

5 参考资料

[1] 《数论概论》P19-P22