为什么辗转相除法能求出最大公约数（转载+修改）

为什么辗转相除法能求出最大公约数

public static int getMaxPact1(int m, int n) {
if (m < n) {
int temp = m;
m = n;
n = temp;
}
// m>=n
int k = 0;
while (m % n != 0) {
k = m % n;
m = n;
n = k;
}
return n;
}

 

上边这段代码就是辗转相除法求最大公约数的程序，很多人都在用这种方法求最大公约数，但是知道为什么这样能求最大公约数的人就不多了。其实这是一个数论的定理：

y = ax + b => gcd(x, y) = gcd(b, x) (a, b, x, y ∈ Z)

只要证明了这个式子成立，只需要证明x、y的公约数和b、x的公约数相同即可，因为如果他们的约数都相同了，那么其中一定有一个是最大的。

 

首先介绍一些简单的数学知识：

y = ax(a, x, y ∈ Z)，就说x能整除y，记作x|y

x|y, x|z => x|ax + bz（这个定理很简单就不做证明了）

算法的证明如下

设d|x, d|y

d|x, d|y => d|y-ax 即 d|b

所以有d|x, d|y => d|x, d|b

即x、y的约数都是b、x的约数

 

假定x,y ∈ Z,x >= y, 令r[0] = x , r[1] = y 用辗转相除法有：

1 2 3 4 5 6 7 8 9 10 11 12 13

r[0] = r[1] * q[1] + r[2]    0<=r[2]<r[1]        r[1] = r[2] * q[2] + r[3]    0<=r[3]<r[2]                     .                     .                     .        r[n-2] = r[n-1] * q[n-1] + r[n]    0<=r[n]<r[n-2]        r[n-1] = r[n]q[n]

 

在序列x=r[0]>r[1]>r[2]...>=0中最多有a项，最终余数为0

由y = ax + b => gcd(x, y) = gcd(b, x) (a, b, x, y ∈ Z)

gcd(x, y) = gcd(r[0], r[1]) = gcd(r[1], r[2])=...=r[n]

即：gcd(x, y )为辗转相除法中组后一个非0的余数

转载自：http://lamoop.com/post/2012-08-09/40033469363