OFFSET
0,3
COMMENTS
The maximum value is 56899, which first occurs at a(275). The maximum next occurs at a(977). T. D. Noe verified that the terms around a(275) and a(977) are the same. Hence the period is 977 - 275 = 702. The actual period starts at a(24) with the interesting terms 349, 45, 45, 349, 349. For some different initial conditions, the period is different. The point of the SORT operation here is that it "mixes" the sequence and the questions are, considering cycles as orbits, all about ergodicity. To turn this into the sorted Fibonacci sequence (A069638), use a(0)=0, a(1)=1, a(2)=1. This is a "base" sequence, but has analogs in other bases; for instance, SORT(base 2)[n] means count the 1's in the binary, call that k and output 2^(k-1). How does this sequence depend on SORT(base M)[n] for various M? Are there any initial values such that the sequence us unbounded? If not, how does cycle length depend upon initial values?
LINKS
Richard I. Hess, Problem 920: sorted Fibonacci sequence, Pi Mu Epsilon Journal, Vol. 10 (Fall 1998) No. 9, pp. 754-755.
FORMULA
a(0) = 1, a(1) = 1, a(2) = 2 and for n>1: a(n+1) = SORT[a(n) + a(n-1) + a(n-2)] where SORT places digits in ascending order and deletes 0.
EXAMPLE
a(8) = 18 because a(5) + a(6) + a(7) = 13 + 24 + 44 = 81 and SORT(81) = 18.
MATHEMATICA
nxt[{a_, b_, c_}]:=Module[{d=FromDigits[Sort[IntegerDigits[a+b+c]]]}, {b, c, d}]; Transpose[NestList[nxt, {1, 1, 2}, 65]][[1]] (* Harvey P. Dale, Feb 07 2011 *)
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
Jonathan Vos Post, Jun 08 2005
STATUS
approved