第一百一十三章 比起数学，我真更擅长物理（3 / 4）

e(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函数有同样的效果，可整个运算过程简化了很多。

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

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

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