我在进一步优化主要计算功能方面有些挣扎。
到目前为止,我最终得到了Eratosthenes筛子。
我在https://primesieve.org/上找到了一个提示,可以通过实施滚轮和与本文的链接来进一步优化此提示: ftp://ftp.cs.wisc.edu/pub/techreports/1991/TR1028.pdf
我试图将此伪代码转换为python,但无法正常工作。我感到步骤B中的迭代不正确。计算prime_sieve_fast(100, 3)
时,不会删除 91 。这是合乎逻辑的,因为运行中的变量永远不会达到7*13
或13*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