实施筛网以筛检Eratosthenes时出现问题

时间:2018-10-29 17:47:15

标签: python primes sieve-of-eratosthenes sieve wheel-factorization

我在进一步优化主要计算功能方面有些挣扎。

到目前为止,我最终得到了Eratosthenes筛子。

我在https://primesieve.org/上找到了一个提示,可以通过实施滚轮和与本文的链接来进一步优化此提示: ftp://ftp.cs.wisc.edu/pub/techreports/1991/TR1028.pdf

我试图将此伪代码转换为python,但无法正常工作。我感到步骤B中的迭代不正确。计算prime_sieve_fast(100, 3)时,不会删除 91 。这是合乎逻辑的,因为运行中的变量永远不会达到7*1313*7。我怎么了?

def prime_sieve(n):
    prime_list=[0,0]
    for i in range(2,n+1):
        prime_list.append(1)
    for p in range(2,int(n**(1/2))+1):
        for j in range(p**2,n+1,p):
            prime_list[j]=0
    primenumbers=[]
    for i in range(n):
        if prime_list[i+1]==1:
            primenumbers.append(i+1)
    return prime_list,primenumbers    

def prime_sieve_faster(n,n_wheel):
    primes=prime_sieve(100)[1][:n_wheel+1]
    w=wheel(n_wheel,primes[:-1])
    M=len(w)
    prime_list=[1]*(n+1)
    for r in range(M):
        if w[r%M]==0:
            b=0
        else:
            b=1
        for i in range(r,n+1,M):
            prime_list[i]=b
    for i in range(n_wheel):
        prime_list[primes[i]]=1
    prime_list[1]=0
    for p in range(primes[n_wheel],int(n**(1/2))+1):
        print(p)
        step=w[p%M]
        if step==0:
            prime_list[p]=0
        else:
            for f in range(p,p+M+1,step):
                for x in range(p*f,n+1,M*p):
                    prime_list[x]=0
                    print(p,f,x,M*p)
    primenumbers=[]
    for i in range(n):
        if prime_list[i+1]==1:
            primenumbers.append(i+1)
    return prime_list,primenumbers


def wheel(k,primes):
    M=1
    for prime in primes:
        M*=prime
    W=[1]*M
    for prime in primes:
        for x in range(0,M,prime):
            W[x]=0
    W[M-1]=2
    prev=M-1
    for x in range (M-2,0,-1):
        if W[x]!=0:
            W[x]=prev-x
            prev=x
    return W

0 个答案:

没有答案