2022年11月29日 星期二

Prime

 import time

n = 2*10**6


stTime = time.time()

p = [2,3,5]

for i in range(6,n,6):

    pflag = True

    for j in p:

        if (i+1)%j==0 :

            pflag = False

            break

        if j*j>i+1:break

    if pflag:p.append(i+1)

    pflag = True

    for j in p:

        if (i+5)%j==0 :

            pflag = False

            break

        if j*j>i+5:break

    if pflag:p.append(i+5)

print(p[-5:])

edTime = time.time()

print(edTime-stTime,len(p))

沒有留言:

張貼留言