质数:一个大于1的自然数,除了1和它自身外,不能被其他自然数整除的数叫做质数;否则称为合数
最近在一个做题网站中遇到的一个问题,编写一个函数,通过传入的参数判断这个数是不是质数。然后发现了三种解决办法,下面按照效率高低依次进行分享。欢迎批评指正。
方法一:
根据解释本身,传入一个参数num,依次判断num - 1中是否有num的约数,如果没有就是质数,否则就是合数。
def is_prime1(num):
i = 2
while i < num-1:
if num % i == 0:
return False
i += 1
return True
第一种办法需要对小于num 的每一个数进行判断,效率最低。
方法二:
对传入的参数num进行因式分解,只判断左半部分,如果左边没有可以被整除的,右边肯定也没有。
from cmath import sqrt
def is_prime2(num):
i = 2
tmp = sqrt(num).real
while i <= tmp:
if num % i == 0:
return False
i += 1
return True
方法三:
关于质数的分布规律,大于等于5的质数一定和6的倍数相邻,例如5和7 11和13 17和19 等等。证明:x >= 1 将x >= 5 的自然数表示如下:6x -2 ,6x -1,6x ,6x+1 ,6x+2 ,6x+3 ,6x+4 ,6x+5 ,6(x+1) ,6(x+1)+1 .....可以看到,不在6的倍数两侧的可以提去一个公约数:6x -2、6x+2、6x+3、6x+4 等于 2(3x-1) 2(3x+1) 3(2x+1) 2(2x+2)。所以他们一定不是素数,在除去6x本身,显然,素数要出现只可能出现在6x的两侧。当然也存在6x两侧的也可能不是素数,还需要在判断一次。而且我们判断的步长可以改成6,
from cmath import sqrt
def is_prime3(num):
if num in [2,3]:
return True
# 不在6的倍数两侧绝对不是质数
elif (num % 6) != 1 and (num % 6) != 5 or num <= 1:
return False
tmp = sqrt(num).real
i = 5
# 在6的倍数两侧也可能不是质数
while i <= tmp:
if num % i == 0 or num % (i+2) == 0:
return False
i += 6
return True
最后对这三个方法进行性能测试,同时处理100000以内的数
通过耗时对比可以看出来,方法二和方法三相差不是很多,但数据量很多的话差距就会明显。方法一简直了。
个人拙见,有更好的解法,欢迎批评指正