数学

primes below 1000000

import math

def find_primes(n):

    sieve = [True] * (n + 1)

    sieve[0], sieve[1] = False, False

    for i in range(2, int(math.sqrt(n)) + 1):

        if sieve[i]:

            sieve[i*i : n+1 : i] = [False] * len(sieve[i*i : n+1 : i])

    primes = [i for i, is_prime in enumerate(sieve) if is_prime]

    return primes

primes = find_primes(1000000)

print(primes)