OFFSET
1,1
MATHEMATICA
Module[{nn=2500, sqf}, sqf=Select[Range[nn], SquareFreeQ]; Select[Prime[ Range[ PrimePi[Length[sqf]]]], PrimeQ[sqf[[#]]]&]] (* Harvey P. Dale, Jul 28 2021 *)
PROG
(PARI) lista(nn) = my(sqfs = select(n->issquarefree(n), vector(nn, i, i))); forprime(p = 2, #sqfs, if (isprime(sqfs[p]), print1(p, ", "))) \\ Michel Marcus, Sep 11 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Irina Gerasimova, Sep 10 2013
EXTENSIONS
Corrected by Harvey P. Dale, Jul 28 2021
STATUS
approved