A106593 Six-symbol substitution based on doubling the Rauzy substitution : n=1 characteristic polynomial: x^6-3*x^4+3*x^2-1.

1, 4, 5, 6, 2, 3, 1, 2, 3, 4, 5, 4, 5, 6, 4, 5, 2, 3, 2, 3, 1, 2, 3, 2, 3, 4, 5, 4, 5, 4, 5, 6, 4, 5, 4, 5, 2, 3, 2, 3, 2, 3, 1, 2, 3, 2, 3, 2, 3, 4, 5, 4, 5, 4, 5, 4, 5, 6, 4, 5, 4, 5, 4, 5, 2, 3, 2, 3, 2, 3, 2, 3, 1, 2, 3, 2, 3, 2, 3, 2, 3, 4, 5, 4, 5, 4, 5, 4, 5, 4, 5, 6, 4, 5, 4, 5, 4, 5, 4, 5, 2, 3, 2, 3, 2

Offset

0,2

Comments

Triangular form: {1}, {4, 5, 6}, {2, 3, 1, 2, 3}, {4, 5, 4, 5, 6, 4, 5}, {2, 3, 2, 3, 1, 2, 3, 2, 3}, {4, 5, 4, 5, 4, 5, 6, 4, 5, 4, 5},

References

Curtis McMullen, Prym varieties and Teichmuller curves.

Links

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

Formula

1->{4.5.6}, 2->{4}*n, 3->{5}, 4->{2}, 5->{3}.n 6->{1, 2, 3}

Mathematica

n0=6 n=1 s[1] = {4, 5, 6}; s[2] = Table[If[i <= n, 4, {}], {i, 1, n0}]; s[3] = Table[If[i <= n, 5, {}], {i, 1, n0}]; s[4] = Table[If[i <= n, 2, {}], {i, 1, n0}]; s[5] = Table[If[i <= n, 3, {}], {i, 1, n0}]; s[6] = {1, 2, 3}; t[a_] := Flatten[s /@ a]; p[0] = {1}; p[1] = t[p[0]]; p[n_] := t[p[n - 1]] aa = Table[p[i], {i, 0, 10}]; MatrixForm[aa] aaa = Flatten[aa]

Crossrefs

Sequence in context: A082486 A106591 A106592 * A276036 A305006 A338688

Adjacent sequences: A106590 A106591 A106592 * A106594 A106595 A106596

Keyword

nonn,uned

Author

Roger L. Bagula, May 10 2005

Status

approved