Contact
CoCalc Logo Icon
StoreFeaturesDocsShareSupport News AboutSign UpSign In
| Download

Play around with sieves and benchmarking...

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