The fastest implementation I could come up with:
isprime = [True]*Nisprime[0] = isprime[1] = Falsefor i in range(4, N, 2): isprime[i] = Falsefor i in range(3, N, 2): if isprime[i]: for j in range(i*i, N, 2*i): isprime[j] = False