我正在尝试用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列表呢?
答案 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)