﻿from math import *

def Erato(n):
    P=list(range(2,n+1))
    j=2
    while j <= sqrt(n):
        if P[j-2]>=1:
            for k in range(2,floor(n/j)+1):
               P[j*k-2]=0
        j=j+1
    L=[]
    for s in P:
       if s!=0:
          L=L+[s]
    return L

