OFFSET
1,1
LINKS
EXAMPLE
Start with
2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83 89 97 101 103 ... and delete every second term, giving
2 5 11 17 23 31 41 47 59 67 73 83 97 103 ... and delete every 3rd term, giving
2 5 17 23 41 47 67 73 97 103 ... and delete every 5th term, giving
.... Continue forever and what's left is the sequence.
MAPLE
S[1]:={seq(ithprime(i), i=1..322)}: for n from 2 to 322 do S[n]:=S[n-1] minus {seq(S[n-1][ithprime(n-1)*i], i=1..nops(S[n-1])/ithprime(n-1))} od: S[322]; # Emeric Deutsch, Nov 17 2004
MATHEMATICA
alle[0]=Table[Prime[i], {i, 1, 10000}]; alle[i_]:=alle[i]= Module[{zuloeschen= Table[alle[i-1][[j]], {j, Prime[i], Length[alle[i-1]], Prime[i]}]}, Complement[alle[i-1], zuloeschen]] (* alle[i] gives the sequence after the i-th iteration and here the first Prime[i] elements are fixed and will not chang in later iterations. So to get the first Prime[10]=29 terms, type *) Take[alle[10], Prime[10]] (* Michael Taktikos, Nov 16 2004 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Nov 16 2004
EXTENSIONS
More terms from Michael Taktikos and Ray Chandler, Nov 16 2004
STATUS
approved