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