A106589 Three-symbol substitution based on Rauzy: adding one 3 to second substitution: characteristic polynomial: x^3-x^2-2*x-2.

2, 3, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 3, 3, 1, 2, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 3, 3, 1, 2, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 1, 2, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 3, 3, 1, 2, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 3, 3, 1, 2, 3, 1, 2, 3, 2, 3, 3, 1, 2, 3, 3, 3, 1, 2, 3, 1, 2

Offset

0,1

Links

Table of n, a(n) for n=0..104.

Curtis McMullen, Prym varieties and Teichmüller curves, Duke Math. J., Volume 133, Number 3 (2006), 405-616

Formula

1->{2}, 2->{3, 3} 3->{1, 2, 3}.

Mathematica

s[1] = {2}; s[2] = {3, 3}; s[3] = {1, 2, 3}; t[a_] := Flatten[s /@ a]; p[0] = {1}; p[1] = t[p[0]]; p[n_] := t[p[n - 1]] aa = p[7]

Prog

(PARI) {a(n)=my(A); if(n<1, 0, A=[2]; while(length(A)<n||A[1]!=2, A=concat(vector(length(A), k, [[2], [3, 3], [1, 2, 3]][A[k]]))); A[n])} /* Michael Somos, May 16 2005 */

Crossrefs

Sequence in context: A371888 A066517 A108132 * A335490 A334290 A294101

Adjacent sequences: A106586 A106587 A106588 * A106590 A106591 A106592

Keyword

nonn,uned

Author

Roger L. Bagula, May 10 2005

Status

approved