Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import numpy as np
- def sieve_of_eratosthenes(limit):
- is_prime = np.ones(limit + 1, dtype=np.bool)
- is_prime[[0, 1]] = False
- for i in range(2, limit + 1):
- if is_prime[i]:
- is_prime[np.arange(i * i, limit + 1, i)] = False
- return np.where(is_prime)[0]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement