Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import numpy as np
- n = int(input("Insert n: "))
- if n > 1:
- A = [True for x in range(n)]
- for i in range(2, int(np.sqrt(n) + 1)):
- if A[i] == True:
- print("i: " + str(i))
- j = np.power(i, 2) + i * (i-2)
- for k in range(n):
- j = np.power(i, 2) + i*k
- if j >= n: break
- print("j: " + str(j))
- A[j] = False
- print("A: " + str(A))
- print("Primes: " + str([i for i, v in enumerate(A) if v and i > 1]))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement