Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- N = 10 ** 7
- def task(N):
- lp = [0] * N
- pr = []
- for i in xrange(2, N):
- if not lp[i]:
- lp[i] = i
- pr.append(i)
- for j, pr_j in enumerate(pr):
- if pr_j <= lp[i] and i * pr_j <= N - 1:
- lp[i * pr_j] = pr_j
- else:
- break
- return pr, lp
- pr, lp = task(N)
- print pr[:10]
Advertisement
Add Comment
Please, Sign In to add comment