OFFSET
1,2
COMMENTS
LINKS
FORMULA
MATHEMATICA
terms = 72;
np = Ceiling[terms/3] + 1;
nc = Ceiling[(2/3) terms];
pp = Prime[Range[np]];
cc = Partition[Select[Range[FindRoot[n == nc + PrimePi[n] + 1, {n, nc, 2nc}][[1, 2]] // Floor], CompositeQ], 2];
Join[{1}, Riffle[pp, cc] // Flatten][[1 ;; terms]] (* Jean-François Alcover, Nov 15 2021 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Reinhard Zumkeller, Jan 20 2006
STATUS
approved