erato

vector<int> erato(int start, int stop)
{
        int i, j;
        vector<int> res;
        vector<int> tab(stop);
        for(i=start; i<=stop; i++)
                tab[i]=i;
        for(i=start; i<=stop; i++){
                if(tab[i]>1){
                        res.push_back(tab[i]);
                        for(j=i; j<=stop; j=j+i)
                                tab[j]=0;
                }
        }
        return res;
}