OFFSET
1,2
COMMENTS
a(9) > 30000. - Michel Marcus, Sep 02 2019
a(13) > 200000. - Giovanni Resta, Sep 02 2019
MATHEMATICA
Do[s=PrimePi[n*Prime[n]]/n; If[IntegerQ[s], Print[n]], {n, 1, 100000}]
Select[Range[16000], IntegerQ[PrimePi[#*Prime[#]]/#]&] (* Harvey P. Dale, Jul 23 2015 *)
PROG
(PARI) isok(n) = denominator(primepi(n*prime(n))/n) == 1; \\ Michel Marcus, Sep 02 2019
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Labos Elemer, May 27 2003
EXTENSIONS
Corrected by Harvey P. Dale, Jul 23 2015
a(9)-a(12) from Giovanni Resta, Sep 02 2019
a(13) from Chai Wah Wu, May 14 2020
STATUS
approved