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.
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