mirror of
https://github.com/python/cpython.git
synced 2025-07-27 21:24:32 +00:00
Improve the sieve() recipe in the itertools docs (gh-109199)
Lazier sieve
This commit is contained in:
parent
75cd86599b
commit
d3ed9921cd
1 changed files with 7 additions and 4 deletions
|
@ -1030,13 +1030,16 @@ The following recipes have a more mathematical flavor:
|
||||||
def sieve(n):
|
def sieve(n):
|
||||||
"Primes less than n."
|
"Primes less than n."
|
||||||
# sieve(30) --> 2 3 5 7 11 13 17 19 23 29
|
# sieve(30) --> 2 3 5 7 11 13 17 19 23 29
|
||||||
|
if n > 2:
|
||||||
|
yield 2
|
||||||
|
start = 3
|
||||||
data = bytearray((0, 1)) * (n // 2)
|
data = bytearray((0, 1)) * (n // 2)
|
||||||
data[:3] = 0, 0, 0
|
|
||||||
limit = math.isqrt(n) + 1
|
limit = math.isqrt(n) + 1
|
||||||
for p in compress(range(limit), data):
|
for p in iter_index(data, 1, start, limit):
|
||||||
|
yield from iter_index(data, 1, start, p*p)
|
||||||
data[p*p : n : p+p] = bytes(len(range(p*p, n, p+p)))
|
data[p*p : n : p+p] = bytes(len(range(p*p, n, p+p)))
|
||||||
data[2] = 1
|
start = p*p
|
||||||
return iter_index(data, 1) if n > 2 else iter([])
|
yield from iter_index(data, 1, start)
|
||||||
|
|
||||||
def factor(n):
|
def factor(n):
|
||||||
"Prime factors of n."
|
"Prime factors of n."
|
||||||
|
|
Loading…
Add table
Add a link
Reference in a new issue