A051373 a(n+1) = a(n) + sum of digits of (a(n)^3).

1, 2, 10, 11, 19, 47, 64, 83, 118, 137, 163, 191, 235, 281, 307, 344, 379, 425, 469, 497, 532, 578, 604, 641, 676, 722, 766, 812, 856, 884, 919, 965, 1009, 1046, 1081, 1109, 1153, 1199, 1252, 1289, 1324, 1352, 1387, 1433, 1486, 1532, 1585, 1640, 1666, 1712

Offset: 1

Author

Miklos SZABO (mike(AT)ludens.elte.hu)

Keywords

• easy
• nonn
• base

Links

• Harvey P. Dale, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A033298.

Programs

Mathematica

NestList[#+Total[IntegerDigits[#^3]]&,1,50] (* Harvey P. Dale, Dec 25 2015 *)

PARI

sod(n) = d = digits(n); sum(i=1, #d, d[i])
lista(nn) = {na = 1; for (i=1, nn, print1(na, ", "); na += sod(na^3););} \\ Michel Marcus, Jun 22 2013

Extensions

More terms from James Sellers, May 01 2000