|
%I #13 Feb 12 2023 15:27:23
%S 1,5,13,22,38,63,90,122,158,207,271,352,452,573,698,826,970,1139,1335,
%T 1551,1776,2019,2275,2564,2888,3231,3592,3992,4433,4917,5429,5958,
%U 6534,7159,7835,8564,9348,10189,11089,12050,13050,14074,15163,16319,17544
%N Sum of first n perfect powers.
%H Hugo Pfoertner, <a href="/A076408/b076408.txt">Table of n, a(n) for n = 1..10000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PerfectPower.html">Perfect Powers</a>.
%e a(8) = 1+4+8+9+16+25+27+32 = 122.
%t Accumulate[Join[{1},Select[Range[1500],GCD@@FactorInteger[#][[All,2]]>1&]]] (* _Harvey P. Dale_, Feb 12 2023 *)
%o (PARI) print1(s=1,", ");for(k=2,1225,if(ispower(k),print1(s+=k,", "))) \\ _Hugo Pfoertner_, Jan 01 2019
%Y Cf. A001597, A025479, A076407, A322969.
%K nonn
%O 1,2
%A _Reinhard Zumkeller_, Oct 09 2002
|
