筛选法求全体素数

计算素数的一个方法是埃氏筛法，它的算法理解起来非常简单,参考《廖雪峰的python教程》:

def _odd_iter():
n = 1
while True:
n = n + 2
yield n

def _not_divisible(n):
return lambda x: x % n > 0

def primes():
yield 2
it = _odd_iter() # 初始序列
while True:
n = next(it) # 返回序列的第一个数
yield n
it = filter(_not_divisible(n), it) # 构造新序列

def print_prime(num):
# 打印num以内的素数:
for n in primes():
if n < num:
print(n)
else:
break

if __name__ == '__main__':
n=int(input("please input an integer:"))
print_prime(n)

注意到Iterator是惰性计算的序列，所以我们可以用Python表示“全体自然数”，“全体素数”这样的序列，而代码非常简洁。

埃氏筛选法算法描述