1

我在这里写的筛算法遇到了问题。我已经尝试修复它总共大约 10 个小时。我在这里四处寻找类似的问题,但似乎找不到遇到此问题的人。我对 python 比较陌生,在阅读了很多生成器文档之后,我设法编写了可以工作的代码。但是,我仍然不知道为什么我的第一次尝试失败了。

我想出的只是在每个连续的筛分步骤中似乎 gen1 实际上并没有被清空。因此,我尝试在名称 gen1 和 gen2 之间交替,删除每个名称以避免此问题。那也没有用。

我真的很感激对此的一些见解,以及任何改进我现在所拥有的建议。

这是失败的代码:

def primes(n):
    "yields primes up to n. For use with large n"
    q = 0
    yield 2
    gen1 = (x for x in range(3,n,2))
    while q*q < n:
        q = next(gen1)
        gen1 = (x for x in gen1 if x%q != 0)
        yield q
    else:
        while 1:
            try:
                yield next(gen1)
            except:
                StopIteration
                break

这是我当前的代码:

import math
global gen1
global gen
def gen1(x):
    for i in range(3,x,2):
        yield i
def gen(generator,n):
    "Input generator and current starting 'index' for the generator"
    # Recursively defines new generator for sieve of Eratosthenes
    for i in range(n+1):
        predicate = next(generator)
        yield predicate
    for i in generator:
        if i % predicate != 0:
            yield i
def primes(n):
    yield 2
    a = gen1(n)
    for i in range(math.ceil(math.sqrt(n))):
        a = gen(a,i)
    yield from a
4

1 回答 1

4

基本问题是范围之一。在这一行:

  gen1 = (x for x in gen1 if x%q != 0)

qused的值不是生成器表达式创建q时绑定的值,而是生成器表达式执行时的(不断变化的!)值。这与在任何类型的嵌套函数中引用非局部变量的方式相同。q

要在创建时捕获绑定,显而易见且简单的方法是编写一个函数,并传递您希望它使用的值。例如,这种重写在这方面和其他几个方面更加 Pythonic:

def primes(n):
    "yields primes up to n. For use with large n"

    def gen(gen1, q):
        for x in gen1:
            if x % q:
                yield x

    q = 0
    yield 2
    gen1 = iter(range(3, n, 2))
    while q*q < n:
        q = next(gen1)
        gen1 = gen(gen1, q)
        yield q
    yield from gen1

现在,主体中使用的值gen()正是传递给它的值,所以它是显而易见的,而不是神秘的 ;-)

于 2018-12-17T00:59:26.167 回答