【#初始化一个布尔数组，所有值都设为True，表示所有数都可能是素数

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， 17， 19， 23， 29， 31

#计算素数间隔

gaps = np.diff(primes)

#计算均值和标准差】

素数定理的近似表达式，它给出了小于或等于某个数$ x $的素数数量\pi(x)π(x)的近似值：[\pi(x)\approx \frac{x}{\ln x}]

这个公式在$ x $很大时非常精确。

在此基础上，直接使用primerange函数，与之前的Python函数有同样的效果，可整个运算过程简化了很多。

伸了个懒腰，许青舟收拾东西回宿舍。

许青舟推门进去，才发现宿舍的氛围非常诡异，孔先波坐在椅子上，双目无神，翟自强看了看他，欲言又止。

就徐正洋还好点，但此时也是直勾勾地盯着他。

许青舟放下东西，担忧地问：「伱们都没事吧？」

「我没事，但老孔他们道心崩溃了。」

徐正洋眼神那叫一个幽怨，咬着牙，说道：「老许，这究竟是tm怎麽回事！」

孔先波和翟自强两个人也看了过来，想听许青舟怎麽说。

「就是两篇论文，不是什麽大事，本来打算等晚点再告诉你们的。」许青舟谦虚地笑了笑。

就两篇论文？不是什麽大事？

徐正洋沉默了，这个时候，他突然觉得，先前老孔那些装逼行为都很容易接受，甚至感到亲切。

妈的，这货才是真正的逼王！

孔先波捂着胸口，有一种心肌梗塞的感觉，翟自强吐出一口气，伸出了一个大拇指。

「这两篇我暑假的时候就开始构思了，顾志锺教授也帮我改了改。」许青舟笑着解释。

顾教授帮着改过？

听到这话，孔先波心里终于好受一点，如果顾教授帮自己改，那.

但是，顾教授为什麽没署名？连通讯作者的名字都没有.

孔先波还没问，翟自强已经发现了重点：「可顾教授是教数学的，怎麽帮你改物理论文？」

「哦，这篇啊，是我自己写的，顾教授只帮我改了数学那篇，不过因为改动不大，所以他说不署名了。」

「其实，上次我没说谎，比起数学，我真更擅长物理。」

(本章完)