prime?=[True?for?_?in?range(limit?+?1)]
p?=?2
while?(p?*?p
#如果prime[p]沒有被改變，那麼它是一個素數
if?(prime[p]==?True):
#更新所有p的倍數為非素數】
可在超算中心運行的時候發現，
當生成小於或等於10000000的所有素數時，即limit?=?10000000，這個函數會消耗大量內存和計算時間。
先前使用自己的電腦模擬並沒有出現這種情況，可使用超算，問題就比較明顯。
那種感覺，就好像頂配的電腦裡邊，出現了一張1050顯卡。
許青舟輕輕轉動手中簽字筆，正在考慮該怎麼修改一下。
時間一分一秒過去。
晚上22：30，許青舟吐出一口氣，看著屏幕，滿意地點了點頭。
【primes?=?np.array([2，?3，?5，?7，?11，?13，?1
7，?19，?23，?29，?31
#計算素數間隔
gaps?=?np.diff(primes)
#計算均值和標準差】
素數定理的近似表達式，它給出了小於或等於某個數$?x?$的素數數量\pi(x)π(x)的近似值：[\pi(x)\approx?\frac{x}{\ln?x}]
這個公式在$?x?$很大時非常精確。
在此基礎上，直接使用primerange函數，與之前的Python函數有同樣的效果，可整個運算過程簡化了很多。
伸了個懶腰，許青舟收拾東西回宿舍。
許青舟推門進去，才發現宿舍的氛圍非常詭異，孔先波坐在椅子上，雙目無神，翟自強看了看他，欲言又止。
就徐正洋還好點，但此時也是直勾勾地盯著他。
許青舟放下東西，擔憂地問：“伱們都沒事吧？”
“我沒事，但老孔他們道心崩潰了。”
徐正洋眼神那叫一個幽怨，咬著牙，說道：“老許，這究竟是tm怎麼回事！”
孔先波和翟自強兩個人也看了過來，想聽許青舟怎麼說。
“就是兩篇論文，不是什麼大事，本來打算等晚點再告訴你們的。”許青舟謙虛地笑了笑。
就兩篇論文？不是什麼大事？
徐正洋沉默了，這個時候，他突然覺得，先前老孔那些裝逼行為都很容易接受，甚至感到親切。
媽的，這貨才是真正的逼王！
孔先波捂著胸口，有一種心肌梗塞的感覺，翟自強吐出一口氣，伸出了一個大拇指。
“這兩篇我暑假的時候就開始構思了，顧誌鐘教授也幫我改了改。”許青舟笑著解釋。
顧教授幫著改過？
聽到這話，孔先波心裡終於好受一點，如果顧教授幫自己改，那.
但是，顧教授為什麼沒署名？連通訊作者的名字都沒有.
孔先波還沒問，翟自強已經發現了重點：“可顧教授是教數學的，怎麼幫你改物理論文？”
“哦，這篇啊，是我自己寫的，顧教授隻幫我改了數學那篇，不過因為改動不大，所以他說不署名了。”
“其實，上次我沒說謊，比起數學，我真更擅長物理。”
(本章完)