Contact
CoCalc Logo Icon
StoreFeaturesDocsShareSupport News AboutSign UpSign In
| Download

Play around with sieves and benchmarking...

Views: 276
Image: ubuntu2004
1
import numpy as np
2
from numba import jit
3
4
5
@jit(nopython=True)
6
def pi():
7
number = 1000000
8
9
#primes = list(range(number + 1))
10
primes = np.arange(number + 1)
11
12
i = 2
13
while i * i <= number:
14
if primes[i] != 0:
15
for j in range(2, number):
16
if primes[i] * j > number:
17
break
18
else:
19
primes[primes[i] * j] = 0
20
i += 1
21
22
cnt = 0
23
for i in range(2, number + 1):
24
if primes[i] != 0:
25
cnt += 1
26
27
return cnt
28
29
30
@jit(nopython=True)
31
def main():
32
j = 0
33
for i in range(100):
34
j += pi()
35
print(j)
36
37
38
main()
39
40