# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a112718 Showing 1-1 of 1 %I A112718 #11 Mar 19 2013 16:01:51 %S A112718 2,12,23,113,151,5924,14254,106545,1915765,2798136,31749441,35282317, %T A112718 35389065,35389165,105227821,141291863,193789064,326730783,839512048, %U A112718 882012907,884676937,2780026914,2997751947,8493184690,8493955191 %N A112718 Numbers n such that pi(n)=1^d_1+2^d_2+...+k^d_k where d_1 d_2 ... d_k is the decimal expansion of n. %C A112718 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. %H A112718 Giovanni Resta, Table of n, a(n) for n = 1..33 (full sequence) %e A112718 326730783 is in the sequence because pi(326730783)=1^3+2^2+3^6+ %e A112718 4^7+5^3+6^0+7^7+8^8+9^3. %t A112718 Do[d=IntegerDigits[n];k=Length[d];If[PrimePi[n]==Sum[j^d[[j]], {j, k}], Print[n]], {n, 410000000}] %Y A112718 Cf. A035138, A112719, A112720. %K A112718 base,fini,full,nonn %O A112718 1,1 %A A112718 _Farideh Firoozbakht_, Sep 17 2005 %E A112718 a(19)-a(25) from _Donovan Johnson_, Nov 09 2010 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE