打印连续素数

Printing consecutive 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列表?

从哪里开始?

  • 您正在编辑一个列表,同时遍历它。那一般是病了 建议 - 可能会出现奇怪的错误。

  • 无论您的情况如何,您都会增加 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)