Python素数列表生成器错误

时间:2019-02-19 09:22:14

标签: python loops primes

此代码从5开始,并列出以下素数供您选择,在本例中为以下17个素数。当我运行它时,将打印25和49。为什么没有将它们过滤掉?

start = 5
number = 1
divisor = 3
upper = start - 2
doc = open("hey.txt", "w")
while number <= 17:
    if start % divisor == 0:
        start = start + 2
        divisor = 3
    elif divisor == upper:
        doc.write(str(start))
        doc.write(", ")
        number = number + 1
        start = start + 2
        divisor = 3
    else:
        divisor = divisor + 2


hey.txt: 5, 7, 11, 13, 17, 19, 23, 25, 29, 31, 35, 37, 41, 43, 47, 49, 53, 

4 个答案:

答案 0 :(得分:1)

您需要更新您的upper变量,我将解释:

start编号写入文件时,这意味着您发现此编号为质数,因此需要将upper变量更新为新的{{1} }值,因为您增加了start -2。因此您的函数应如下所示:

start

答案 1 :(得分:0)

可能是因为您没有增加鞋帮。因此,您测试的最大除数为3。

答案 2 :(得分:0)

您已经准备好回答您的问题。但是看看这个问题。这是一种进行素数Fast primes with python的更快方法:

def primes2(n):
    """ Input n>=6, Returns a list of primes, 2 <= p < n """
    n, correction = n-n%6+6, 2-(n%6>1)
    sieve = [True] * (n//3)
    for i in range(1,int(n**0.5)//3+1):
      if sieve[i]:
        k=3*i+1|1
        sieve[      k*k//3      ::2*k] = [False] * ((n//6-k*k//6-1)//k+1)
        sieve[k*(k-2*(i&1)+4)//3::2*k] = [False] * ((n//6-k*(k-2*(i&1)+4)//6-1)//k+1)
    return [2,3] + [3*i+1|1 for i in range(1,n//3-correction) if sieve[i]]

是的,它不是很可读。但是首先,奥尔威特先看看别人是否已经做了您想做的事情。 Python是一种具有庞大社区的语言,它的好处是拥有大量的库和大量的开发的算法/程序,它们可能比您自己更干净,更快速地完成您想做的事情。因此,请享受社区的利益。)

答案 3 :(得分:0)

更短更短:

lower = 5
upper = 50

print("Prime numbers between {} and {} are: ".format(lower, upper))

doc = open("hey.txt", "w")
for num in range(lower,upper + 1):
   # prime numbers are greater than 1
   if num > 1:
       for i in range(2,num):
           if (num % i) == 0:
               break
       else:
           print(num)
           doc.write(str(num))
           doc.write(",")

输出:

Prime numbers between 5 and 50 are: 
5,7,11,13,17,19,23,29,31,37,41,43,47,