【#初始化一个布尔数组,所有值都设为True,表示所有数都可能是素数
prime=[Truefor_inrange(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)approxfrac{x}{lnx}]
这个公式在$x$很大时非常精确。
在此基础上,直接使用primerange函数,与之前的Python函数有同样的效果,可整个运算过程简化了很多。
伸了个懒腰,许青舟收拾东西回宿舍。
许青舟推门进去,才发现宿舍的氛围非常诡异,孔先波坐在椅子上,双目无神,翟自强看了看他,欲言又止。
就徐正洋还好点,但此时也是直勾勾地盯着他。
许青舟放下东西,担忧地问:「伱们都没事吧?」
「我没事,但老孔他们道心崩溃了。」
徐正洋眼神那叫一个幽怨,咬着牙,说道:「老许,这究竟是tm怎麽回事!」
孔先波和翟自强两个人也看了过来,想听许青舟怎麽说。
「就是两篇论文,不是什麽大事,本来打算等晚点再告诉你们的。」许青舟谦虚地笑了笑。
就两篇论文?不是什麽大事?
徐正洋沉默了,这个时候,他突然觉得,先前老孔那些装逼行为都很容易接受,甚至感到亲切。
妈的,这货才是真正的逼王!
孔先波捂着胸口,有一种心肌梗塞的感觉,翟自强吐出一口气,伸出了一个大拇指。
「这两篇我暑假的时候就开始构思了,顾志锺教授也帮我改了改。」许青舟笑着解释。
顾教授帮着改过?
听到这话,孔先波心里终于好受一点,如果顾教授帮自己改,那。
但是,顾教授为什麽没署名?连通讯作者的名字都没有。
孔先波还没问,翟自强已经发现了重点:「可顾教授是教数学的,怎麽帮你改物理论文?」
「哦,这篇啊,是我自己写的,顾教授只帮我改了数学那篇,不过因为改动不大,所以他说不署名了。」
「其实,上次我没说谎,比起数学,我真更擅长物理。」
(本章完)