|
|
A051373
|
|
a(n+1) = a(n) + sum of digits of (a(n)^3).
|
|
9
|
|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
MATHEMATICA
|
NestList[#+Total[IntegerDigits[#^3]]&, 1, 50] (* Harvey P. Dale, Dec 25 2015 *)
|
|
PROG
|
(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
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn,base
|
|
AUTHOR
|
Miklos SZABO (mike(AT)ludens.elte.hu)
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|