我在堆栈溢出中找到了针对Python的Eratosthenes筛的高度优化的实现。我对它的功能有一个大概的了解，但是我必须承认它的运作细节让我难以理解。

我仍然想将其用于一个小项目（我知道有一些库可以做到这一点，但我想使用此功能）。

这是原始的：

'''
    Sieve of Eratosthenes 
    Implementation by Robert William Hanks      
    https://stackoverflow.com/questions/2068372/fastest-way-to-list-all-primes-below-n/3035188
'''

def sieve(n):
    """Return an array of the primes below n."""
    prime = numpy.ones(n//3 + (n%6==2), dtype=numpy.bool)
    for i in range(3, int(n**.5) + 1, 3):
        if prime[i // 3]:
            p = (i + 1) | 1
            prime[       p*p//3     ::2*p] = False
            prime[p*(p-2*(i&1)+4)//3::2*p] = False
    result = (3 * prime.nonzero()[0] + 1) | 1
    result[0] = 3
    return numpy.r_[2,result]

我想要实现的是对其进行修改，以返回n 从 以下 开始的 所有素数 ，x以便：

primes = sieve(50, 100)

会返回介于50和100之间的质数。这 似乎 很容易，我尝试替换了这两行：

def sieve(x, n):
    ...
    for i in range(x, int(n**.5) + 1, 3):
    ...

但是由于我无法解释的原因，上述的值对x返回的numpy数组没有影响！

如何修改sieve()为 仅 在x和之间返回素数n