Recipe 576596 revision 3

def eratosthenes(n):
    N=range(n+1)
    z=[0]*(n/2)
    for i in range(2, int(n**.5)+1):
        if N[i]:
            N[i*i::i] = z[:(n/i)-i+1]
    return filter(None, N[2:])

• « Back to Recipe 576596

History

• revision 3 (15 years ago)
• previous revisions are not available