A071421 a(n) = a(n-1) + sum of decimal digits of n^n.
1, 5, 14, 27, 38, 65, 90, 127, 172, 173, 214, 268, 326, 378, 477, 565, 663, 771, 898, 929, 1046, 1194, 1340, 1493, 1644, 1798, 1987, 2150, 2317, 2380, 2564, 2769, 2976, 3190, 3450, 3720, 3991, 4256, 4562, 4674, 4982, 5297, 5610, 5935, 6241, 6593, 6967
Offset: 1
Links
- Harvey P. Dale, Table of n, a(n) for n = 1..1000
Programs
-
Mathematica
s=0; Do[s=s+Apply[Plus, IntegerDigits[n^n]]; Print[s], {n, 1, 128}] nxt[{n_,a_}]:={n+1,a+Total[IntegerDigits[(n+1)^(n+1)]]}; NestList[nxt,{1,1},50][[All,2]] (* Harvey P. Dale, Dec 11 2016 *)