您好,登錄后才能下訂單哦!
腦袋一熱,想看一下300以內的最大素數是多少,就寫了一個。
注意:對正整數n,如果用2到n的平方根之間的所有整數去除,均無法整除,則n為質數(素數)。
# -*- coding:utf-8 -*-
import math
import time
ss = [] # 放可能是素數的列表
fss = [] # 放可能是非素數的列表
result = [] # 最終結果
x = 300
print("start!" )
start = time.clock()
# 遍歷所有小于X,大于2的數
for xx in range(2, x+1):
# 只要xx的數,不能被2至xx的平方根的所有數整除,就是素數
for i in range(2, int(math.sqrt(xx)+1)):
if (xx % i) != 0:
i = i + 1
# print("素數: ",xx)
ss.append(xx)
else:
# print("非素數", xx)
fss.append(xx)
# 只要x中的數沒有出現在非素數列表中,則它就是素數
for j in range(2, x+1):
if j not in fss:
result.append(j)
print("result: ", result)
end = time.clock()
times = ''.join("%f s" % (end - start)) # 看用時多久
print("end!\n耗時:", times)
結果:
start!
result: [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, 293]
end!
耗時:0.002066 s
百度百科里有一個比較簡單:而且速度也快。
# 判斷是否為素數
def is_prime(n):
if n == 1:
return False
for i in range(2, int(math.sqrt(n)+1)):
if n % i == 0:
return False
return True
start = time.clock()
print("start:" )
for n in range(2, x+1):
if is_prime(n) is True:
result.append(n)
print(result)
end = time.clock()
times = ''.join("%f s" % (end - start))
print("end!\n耗時:", times)
結果
start:
[2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, 293]
end!
耗時: 0.000434 s
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。