OFFSET
1,1
COMMENTS
The largest term is less than 10^12 because if m>12 then 1^9+2^9+...+n^9 < pi(10^(m-1)). There is no further term up to 41*10^7.
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..33 (full sequence)
EXAMPLE
326730783 is in the sequence because pi(326730783) = 1^3 + 2^2 + 3^6 + 4^7 + 5^3 + 6^0 + 7^7 + 8^8 + 9^3 = 17618732.
MATHEMATICA
Do[d=IntegerDigits[n]; k=Length[d]; If[PrimePi[n]==Sum[j^d[[j]], {j, k}], Print[n]], {n, 410000000}]
CROSSREFS
KEYWORD
base,fini,full,nonn
AUTHOR
Farideh Firoozbakht, Sep 17 2005
EXTENSIONS
a(19)-a(25) from Donovan Johnson, Nov 09 2010
STATUS
approved