打印连续的素数

时间:2018-08-24 15:16:12

标签: python primes

我正在尝试用Python创建一个程序,以打印连续的素数。

estPrimes = [2, 3]
num = 3
run = True

while run == True:
    for i in range(estPrimes[0], estPrimes[(len(estPrimes) - 1)]):
        if i % num == 0:
            num += 2
            sleep(1)
        else:
            estPrimes.append(num)
            print(num)
            sleep(1)
            num += 2

它只会以两倍的打印量结束,我知道它与num += 2有关,但是为什么不除以estPrimes列表呢?

1 个答案:

答案 0 :(得分:1)

从哪里开始?

  • 您正在迭代一个列表。那通常是病 建议-可能会出现奇怪的错误。

  • 无论情况如何,您都可以增加num,这是您的主要问题

  • 您正在检查找到的第一个和最后一个质数之间的每个数字,而不仅仅是质数


这是您代码的修正

def print_primes_to(n):

    primes = [2,3]
    for prime in primes:
        print(prime)

    for num in range(5, n, 2):
        isprime = True
        for prime in primes:
            if num % prime == 0: #not prime
                isprime = False
                break
        if isprime:  # could replace with for/else
            primes.append(num)
            print(num)

测试用例

>>> print_primes_to(30)
2
3
5
7
11
13
17
19
23
29

有更好的算法,这是我闲逛的算法:

import itertools as _itertools

def primelist(n):
    """Returns a list of all primes less than n."""
    if n <= 2:
        return []
    nums = [True]*(n - 2)
    for num in range(2, _rootp1(n)):
        for counter in _itertools.count(num):
            product = counter * num
            if product >= n:
                break
            nums[product - 2] = False
    return [index for index, val in enumerate(nums, 2) if val]

def _rootp1(n: float) -> int:
    """returns the truncated square root of n plus 1 for use in range"""
    return int(n**0.5) + 1

看看生成小于100万的所有素数需要多长时间

%timeit primelist(1_000_000)  
1.33 s ± 30.6 ms per loop (mean ± std. dev. of 7 runs, 1 loop each)