OFFSET
1,1
MATHEMATICA
m = 10^4; L = Table[2*i + 1, {i, 0, m}]; For[n = 2, n < Length[L], r = L[[n++]]; L = ReplacePart[L, Table[r*i -> Nothing, {i, 1, Length[L]/r}]]]; ll = L[[Select[L, # <= Length[L] &]]]; pp = Prime@ Prime@ Range@ PrimePi@ PrimePi@ (2m); Intersection[pp, ll] (* after Jean-François Alcover at A000959 and Giovanni Resta at A303403 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 19 2019
STATUS
approved