第一百一十三章 比起数学，我真更擅长物理 (2/3)

2014年11月5日。

    很多期刊都会提前印刷，这样读者在10月收到11月的期刊时，并不会感觉是过期的。

    还有，群里被消息轰炸了。

    当然，主要是徐正洋的消息，孔先波和翟自强两个人陷入了诡异的安静。

    他发了一个晚上说，收起手机，去吃饭了。

    他正在学三食堂吃饭的时候，季云石就打来电话，说第二专业的手续已经办妥了，他平时上课就跟着物理学专业1班。

    至于要不要去听课，就看许青舟这边有没有时间。

    吃完饭，许青舟找了间空教室，如今构造出来的素数分析模型，只是一个底层模型，说起来的话算是地基，接下来需要添砖加瓦。

    不过，下午在超算中心模拟数值的时候，许青舟发现一个小问题。

    在构造底层模型的时候，为了更好的筛选素数，查找规律，许青舟编写了一个Python函数。

    【#初始化一个布尔数组，所有值都设为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 

本章未完，请点击下一页继续阅读