3.7 Python基础之综合练习一
3.7.1 最小公倍数算法
最大公约数(GCD, Greatest Common Divisor
),也称最大公因数、最大公因子,指两个或多个整数共有约数中最大的一个。
比如数12
和数18
的最大公约数是6
,因为12
的约数有1、2、3、4、6、12
,而18
的约数有1、2、3、6、9、18
,通过比较,显然6
是数12
和数18
的最大公约数。
通过上述过程,显然我们可以通过枚举这两个数的所有约数,考虑这两个数共有的约数,然后选择最大的就是这两个数的最大公约数。因为一个数的约数必然是不大于该数的,所以我们可以通过枚举不超过这两个数中的最大者的正整数,来达到上述效果,具体代码如下述所示:
def gcd_1(x, y):
ed = max(x, y)+1
divisor = 1
for i in range(2, ed):
if x % i == 0 and y % i == 0:
divisor = i
return divisor
其实在古代就有能求解出最大公约数的算法了,《九章算术》是中国古代的数学专著,其中的“更相减损术”就可以用来求两个数的最大公约数,原文是:“可半者半之,不可半者,副置分母、子之数,以少减多,更相减损,求其等也。以等数约之。”大致所描述的算法步骤是:
- 任意给定两个正整数;判断它们是否都是偶数。若是,则用
2
约简;若不是则执行第二步; - 以较大的数减较小的数,接着把所得的差与较小的数比较,并以大数减小数。继续这个操作,直到所得的减数和差相等为止;
- 第一步中约掉的若干个
2
与第二步中最后得到的差(或减数)的乘积就是所求的最大公约数。
实际编程中,我们可以省略第一步,这样第二步最后得到的差(或减数)就是这两个数的最大公约数,其具体实现如下述代码所示:
def gcd_2(x, y):
while True:
if x < y:
x, y = y, x
elif x == y:
return x
x -= y
尽管前面已经介绍了两种求最大公约数的方法,但实际生活中,我们更倾向于使用辗转相除法,来求解任意两个正整数的最大公约数,以求解30
和12
的最大公约数为例,按gcd_2
代码,其过程为:
30 - 12 = 18 -> 18 - 12 = 6 -> 12 - 6 = 6
最后因为减数和差相等,即6 - 6 = 0
,故6
就是30
和12
的最大公约数,仔细观察上述过程,我们可以发现第一步和第二步实际上就是被减数30
减了2
次12
,然后在第三步,用上次计算的余数6
继续与12
进行比较。显然,我们可以通过整数求余运算,直接一步求得30
和12
的余数6
,此时余数绝对是比除数小的,那么则将除数代替被除数的位置,余数代替除数的位置,然后重复上述过程,直至余数为0
,那么此时的除数就是原来两个数的最大公约数了。上述过程用递归方式实现的话,代码是非常简短的,具体代码如下:
def gcd(x, y):
return x if y == 0 else gcd(y, x%y)
如何求任意两个正整数的最小公倍数
几个数共有的倍数叫做这几个数的公倍数,其中除0
以外最小的一个公倍数,叫做这几个数的最小公倍数(LCM,Least Common Multiple
)。
如3
和7
的最小公倍数是21
,因为不存在一个比21
还小的正整数,既是3
的倍数,也是7
的倍数。
显然对任意两个正整数a
和b
,a*b
必是他们的公倍数。假设g
为a
和b
的最大公约数,那么a
、b
可以分别写成一个正整数与他们最大公约数的乘积的形式,即a = p * g
,b = q * g
。那么显然c = p * q * g
,是a
和b
的一个公倍数,而且是最小公倍数。因为p
和q
必定不共有大于1
的公约数,所以若减小p、q、g
这三个任意一个数的话,都不能使其乘积还是a
和b
的倍数。
3.7.2 输出指定范围内的素数
如何判断一个正整数是否是素数
素数(Prime Number
),又称质数,一个大于1
的自然数,除了1
和它自身外,不能整除其他自然数的数叫做质数;否则,称为合数(Composite Number
)。1
既不是素数,也不是合数。
如2、3、5、7、11
都是素数,因为找不到除了1
和其本身之外的约数;而4、6、8
都是合数,因为4
可以整除2
,6
可以整除2
和3
,8
可以整除2
和4
。
根据上述定义,我们很容易写出判断一个素数是否是素数的代码:
def is_prime_1(x):
if x == 1:
return False
for i in range(2, x):
if x % i == 0:
return False
return True
假设一个正整数a
,则其可以被写成任意两个正整数之积,即a = p * q
,假设p < q
,那么正整数p
和q
都是a
的约数,注意到,如果我们知道p
是a
的约数,那么可以通过q = a / p
快速求得另外一个约数q
。所以,我们在判断质数的时候,只需要枚举2
到不大于sqrt(a)
的正整数即可。
虽然通过上述方法,已经能让我们在根号级别的复杂度内,判断一个正整数是否是素数,但如果我们要判断很多个数是否为素数呢?是否每次都需要枚举int(sqrt(a)+1)
个数呢?回到我们最初的起点,我们之所以要枚举这些数,就是想找出原数的约数。然后除1
外,任何一个正整数都能写成多个素数的乘积的形式,那么我们枚举特定范围内的所有素数,也能达到相同的效果,而且在判断多个正整数是否是素数的时候,我们只需要枚举更少的质因数与其比较。大家可以看下面不同区间内的素数统计结果:
从上图的统计结果我们可以发现,当区间越来越大,里面的素数个数和区间内所有数字的个数差距也越来越大。所以,我们用区间内的素数,去判断一个整数是不是素数,比较的次数将更少。
而求不超过某个正整数x
内的所有素数,有一个著名的算法——埃拉托斯特尼筛法。其算法描述为:
- 先用一个数组
vis
,把不大于该正整数x
的所有正整数标记为0
,表示没有访问; - 然后从第一个素数
2
开始遍历整个区间,如果当前访问的数没有访问过,则可以认为它是一个素数,那么就将它在该区间内所有的倍数全部标记为已访问,这样就保证外部的循环发现的没有访问过的数都是素数。
其具体实现如下述代码所示:
def sieve(x):
vis = [0 for i in range(x+1)]
prime_table = []
for i in range(2, x+1):
if vis[i] == 0:
prime_table.append(i)
for j in range(i*2, x+1, i):
vis[j] = 1
return prime_table
然而,除了上述筛法,还有其他高效的筛法,比如欧拉筛法。
def ouler(x):
vis = [0 for i in range(x+1)]
prime_table = []
ln = 0
for num in range(2, x+1):
if vis[num] == 0:
prime_table.append(num)
ln += 1
for j in range(ln):
if num * prime_table[j] > x:
break
vis[num * prime_table[j]] = 1
if num % prime_table[j] == 0:
break
return prime_table