【#初始化一个布尔数组，所有值都设为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怎么回事！”
孔先波和翟自强两个人也看了过来，想听许青舟怎么说。
“就是两篇论文，不是什么大事，本来打算等晚点再告诉你们的。”许青舟谦虚地笑了笑。
就两篇论文？不是什么大事？
徐正洋沉默了，这个时候，他突然觉得，先前老孔那些装逼行为都很容易接受，甚至感到亲切。
妈的，这货才是真正的逼王！
孔先波捂着胸口，有一种心肌梗塞的感觉，翟自强吐出一口气，伸出了一个大拇指。
“这两篇我暑假的时候就开始构思了，顾志钟教授也帮我改了改。”许青舟笑着解释。
顾教授帮着改过？
听到这话，孔先波心里终于好受一点，如果顾教授帮自己改，那.
但是，顾教授为什么没署名？连通讯作者的名字都没有.
孔先波还没问，翟自强已经发现了重点：“可顾教授是教数学的，怎么帮你改物理论文？”
“哦，这篇啊，是我自己写的，顾教授只帮我改了数学那篇，不过因为改动不大，所以他说不署名了。”
“其实，上次我没说谎，比起数学，我真更擅长物理。”
(本章完)
本章已阅读完毕(请点击下一章继续阅读!)