Python学习日志——filter函数埃氏筛法

本文仅代表个人理解，如有错误还请不吝赐教

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) # 构造新序列
# 打印1000以内的素数:
for n in primes():
    if n < 1000:
        print(n)
    else:
        break

计算素数的方法是埃氏筛法，题目参照廖雪峰学习平台

算法思想：

首先，列出从2开始的所有自然数，构造一个序列：

2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, ...

取序列的第一个数2，它一定是素数，然后用2把序列的2的倍数筛掉：

3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, ...

取新序列的第一个数3，它一定是素数，然后用3把序列的3的倍数筛掉：

5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, ...

取新序列的第一个数5，然后用5把序列的5的倍数筛掉：

7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, ...

不断筛下去，就可以得到所有的素数。

首先我们先看最开始运行的部分

for n in primes():

在这里，我们会进入到primers()函数里，并且根据这行语句，我们知道：primers()函数得出的一定是可迭代的，这也就是yield的作用

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

接着运行到n=next(it)语句后，开始运行_odd_iter()函数，返回的n为3.

        it = filter(_not_divisible(n), it) # 构造新序列

紧接着是这句语句，值得注意的是，此时并没有真正运行这行语句，而是等待着it函数的运行，也就是_odd_iter()的运行，因为只有next(it)时才会运行，在运行_odd_iter()后会进行这行语句，也就是跳转到_not_divisible(n)函数中

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

而在这里，我们要注意的是，n是指代不算刚刚产生的所有n，而x是刚刚产生的n，而这个代码就是将x与之前产生过的所有素数进行取余以确定x是否是素数的倍数。

yield n
        it = filter(_not_divisible(n), it) # 构造新序列
# 打印1000以内的素数:
for n in primes():
    if n < 1000:
        print(n)
    else:
        break

之后就再度输出可迭代的n,进行判断，之后再循环直到完成为止。

总结：

1、需要明白在yield-next下，函数从next开始运行，到yield结束，而牢记这点，这个代码理解并不困难。

2、在需要调用匿名函数lambda时，要注意"lambda x:"中的x是需要外界输入的，而如果不需要输入是没有这个变量的。如下图所示：

    def f(j):
        return lambda :j*j

3、要明白yield产生的是迭代，n在不同函数中代表也不相同