python3 题解（4）最大公约数

最大公约数

【问题】整数域上，a的约数指的是：所有能都够整除a的数字。a，b的公约数指：所有能同时整除a，b两个数字的数。a，b的最大公约数，毫无疑问就是这些数字中最大的一个。
比如： 5 8 最大公约数为1, 15 40 最大公约数是5。
编程求 a，b的最大公约数。

这个问题最笨的解法是从较小的一数开始，一直往下测试每一个数字的整除可能性。找到第一个能同时整除a，b的即是解。

当然，数字一大，上法就太慢了。可以用中国古代就有的辗转相除法，西方叫欧几里德定理。
其原理是：如果a==b， 则a为所求。否则不妨设a>b, 则 b, a-b 的最大公约数为所求。

用python 表述为：

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

if __name__ == '__main__':
	print(gcd(30, 125))
	print(gcd(0, 8))
	print(gcd(8, 0))

代码工作固然不错。但如果考虑整个整数域，仍有问题。
-5 与 0 的最大公约数？
按定义，能同时除尽它们的最大数字是：5，而不是 -5。
0 与 0 的最大公约数呢？ 多大的数字都不够大，因而为无穷大。也可以说无意义。
0.3 与 5 呢？ 那就过分了。总要是整数才讨论公约数啊。

python 是解释型的脚本语言。它的好处是灵活，类型动态。但参数与返回值的类型就没有什么保证。为了安全都要自己在运行时检查。所以函数的设计要仔细。

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

	if not (isinstance(a,int) and isinstance(b,int)): 
		raise TypeError('a,b must be int')
	if a==0 and b==0:
		raise ValueError('a, b can not all be zero')
	return f(abs(a),abs(b))

if __name__ == '__main__':
	print(gcd(30, 125))
	print(gcd(0, 8))
	print(gcd(8, 0))
	print(gcd(-30, 20))
	print(gcd(30, -35))
	print(gcd(-28, -35))

除了最大公约数问题，还有个它的孪生问题：最小公倍数。
只要利用下面的公式就好了。

$$a，b的最小公倍数=\frac{ab}{gcd(a,b)}$$

注意，这是只对正整数成立的。最小公倍数应该定义为正数，否则，负数可以无穷大，找不到一个所谓最小的。