OFFSET
1,1
COMMENTS
Iteration g[x] applied in A031176 is slightly modified to obtain actual function to iterate here: f[x]=1+g[x].Initial values resulting in fixed points are collected.
EXAMPLE
n=35 is here because list={36,46,53,[35],35,...} with transient t=3,
c=1 cycle-length;
MATHEMATICA
ed[x_] :=IntegerDigits[x]; f[x_] :=Apply[Plus, ed[x]^2]+1; itef[x_, ho_] :=NestList[f, x, ho]; tmc=Table[Length[Union[itef[w, 100], {w, 1, 256}]; c1=Table[Min[Flatten[Position[itef[w, Length[Union[itef[w, 100]]]] -Last[itef[w, Length[Union[itef[w, 100]]]]], 0]]], {w, 1, 256}]; Flatten[Position[c1, 1]]
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Labos Elemer, Nov 11 2004
STATUS
approved