OFFSET
1,1
COMMENTS
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}]]]; P = Select[ Range[2m], PrimeQ]; lp = L[[Select[P, # <= Length[L] &]]]; pl = P[[Select[L, # <= Length[P] &]]]; pp = P[[Select[P, # <= Length[P] &]]]; ll = L[[Select[L, # <= Length[L] &]]]; Intersection[lp, pl, pp, ll] (* after Jean-François Alcover at A000959 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 19 2019
STATUS
approved