求最大公约数的奇妙技巧

最大公约数

题目：
求正整数a和b的最大公约数gcd(a, b)

普通方法：
假设a < b,遍历i从a到1，如果a和b均能整除i，则i即为所求。

def gcd(a, b):
	if a > b: a, b = b, a
	for i in range(a, 0, -1):
		if a % i == 0 and b % i == 0: return a

奇妙技巧：
对于两个数a和b，必定存在整数k和r，使得a = k * b + r。
设b和r的最大公约数为gcd(b, r)，令b = k₁gcd(b, r)，r = k₂gcd(b, r)，带入上式，可得:
a = k * k₁gcd(b, r) + k₂gcd(b, r) = (kk₁+k₂)gcd(b, r)，从而gcd(b, r)也是a的约数。
又gcd(b, r)是b的约数，因此gcd(b, r)是a和b的公约数。
设gcd(a, b)是a和b的最大公约数，从而必定有gcd(a, b) % gcd(b, r) = 0

将a = k * b + r改写成r = a - k * b，同理可知, gcd(a, b) 一定是r的约数，从而gcd(a, b)是a和r的公约数。
从而必定有gcd(b, r) % gcd(a, b) = 0

可得：gcd(a, b) % gcd(b, r) = gcd(b, r) % gcd(a, b)，
因此：gcd(a, b) = gcd(b, r) = gcd(b, a % b)
注意到：当r = 0, 即a % b = 0时，b即为最大公约数，很明显，可以用递归求解。

def gcd(a, b):
	return a if b == 0 else gcd(b, a % b)