A105264 Theta(1) Pisot substitution level 7 : characteristic polynomial x^4-x^3-1=0.

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

Offset

0,2

Comments

Program for getting Polynomial: s[1] = {2, 0, 0, 0}; s[2] = {3, 0, 0, 0}; s[3] = {4, 0, 0, 0}; s[4] = {4, 1, 0, 0}; M = Table[Table[Count[s[j], i], {i, 1, n0}], {j, 1, n0}] Det[M - x*IdentityMatrix[n0]]

Links

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

Eric Weisstein's World of Mathematics, Pisot Number

Formula

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

Mathematica

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

Crossrefs

Sequence in context: A348190 A211100 A329326 * A063787 A307092 A335458

Adjacent sequences: A105261 A105262 A105263 * A105265 A105266 A105267

Keyword

nonn,uned

Author

Roger L. Bagula, Apr 15 2005

Status

approved