用欧拉筛法(线性筛法),更简洁了。
import time
def EchoTime(info):
# 格式化成2016-03-20 11:45:39形式
print(info)
print(time.strftime("%Y-%m-%d %H:%M:%S", time.localtime()))
def FindPrime(p_yield):
p=[2,3]
t=5
#EchoTime('开始时间:')
while p[-1]<p_yield:
k=0
for i in p:
if t%i!=0:
k+=1
else:
break
if k==len(p):
p.append(t)
t=p[-1]+2
else:
t+=2
#if len(p)%1000==0:
# print(p[-1])
# print(time.strftime("%Y-%m-%d %H:%M:%S", time.localtime()))
print(p)
#EchoTime('运算结束时间:')
if __name__ == '__main__':
FindPrime(1000000)
#import profile
#profile.run("FindPrime(1000000)")
以下是旧算法。
_____________________________________________________________________
def FindPrime(PrimeSeek,p_yield):
p=[2,3,5,7,11,13]
x=PrimeSeek
while x<p_yield:
temp=int(sqrt(x))
IsPrime=True #初始化判定标记
k=len(p)
while(temp<p[k-1:k][0]):
k-=1
for i in p[:k]:
if x%i!=0:
continue
else:
IsPrime=False
break
if(IsPrime==True): #如果到循环结束仍然为真则为质数
p.append(x)
x+=1
print(p)
在Idle内逐行输入
from math import sqrt
FindPrime(17,200)
17是种子值,需要是个质数。
200表示生成200以内的质数。
Python 3.5.1测试通过。