viết chương trình đếm số các số nguyên tố < 10000 . python

Ta có công thức tính số lượng số nguyên tố $< n$ là $\dfrac{n}{\ln(n)}$

=> Để sàng $10^4$ số nguyên tố thì $n = 10^6$ nên ta cần sàng $\dfrac{10^6}{\ln(10^6)}\approx72382$ số, đủ vì $72382 > 10000$

from math import sqrt
from math import trunc

prime = [True for i in range(1000000 + 1)]

def sieve():
    prime[0] = prime[1] = False
    for i in range(2, 1001):
        if prime[i]:
            for j in range(i * i, 1000000 + 1, i):
                prime[j] = False

sieve()

cnt = 0
for i in range(2, 10000):
    if prime[i]:
        cnt += 1

print(cnt)