OFFSET
1,2
COMMENTS
a(28) > 2*10^6. [Donovan Johnson, May 30 2010]
EXAMPLE
n = 14: a(14) = 35950, Pi(35950^2) = Pi(1292402500) = 64866916 = 8054^2
MATHEMATICA
Do[ If[ IntegerQ[ Sqrt[ PrimePi[n^2]]], Print[n]], {n, 1, 400000} ]
lst = {}; Do[ If[ IntegerQ@ Sqrt@ PrimePi[n^2], AppendTo[lst, n]; Print@n], {n, 520000}]; lst (* Robert G. Wilson v, Feb 03 2006 *)
PROG
(PARI) pi(x, c=0) = forprime(p=2, x, c++); c for(n=1, 10^6, if(issquare(pi(n^2)), print1(n, " ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Oct 07 2001
EXTENSIONS
More terms from Robert G. Wilson v and Labos Elemer, Oct 08 2001
Further terms from Robert G. Wilson v, Oct 16 2001
a(20) from Robert G. Wilson v, Feb 03 2006
Added missing terms 134141, 142413, 194711 and a(24)-a(27) from Donovan Johnson, May 30 2010
a(28)-a(31) from Chai Wah Wu, Aug 27 2019
STATUS
approved