擅长:python、mysql、java
<p>更快的素数列表是这样的:</p>
<pre><code>import numpy as np
def primesfrom2to(n):
""" Returns a array of primes, p < n """
assert n>=2
sieve = np.ones(n/2, dtype=np.bool)
for i in xrange(3,int(n**0.5)+1,2):
if sieve[i/2]:
sieve[i*i/2::i] = False
return np.r_[2, 2*np.nonzero(sieve)[0][1::]+1]
</code></pre>
<p>那么就是这些的总和:</p>
^{pr2}$
<p>如果要使用函数,建议重写:</p>
^{3}$