OFFSET
1,1
MATHEMATICA
Do[k=Round[Log[n]]; If[PrimePi[n^2]==PrimePi[(n-k)^2]+n, Print[n]], {n, 1, 10^5, 1}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Manuel Valdivia, Mar 13 2007
STATUS
approved