OFFSET
0,2
MATHEMATICA
f[n_] := Plus @@ (IntegerDigits[n]^2); t = Table[0, {25}]; k = 1; While[k < 150000001, a = Length@ NestWhileList[f, k, UnsameQ@## &, All] - 1; If[a < 25 && t[[a]] == 0, t[[a]] = k; Print[{a, k}]]; k++ ]
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Robert G. Wilson v, Apr 25 2010
STATUS
approved