%I #4 Mar 12 2014 16:36:46
%S 1,1,1,2,1,1,2,1,2,1,2,1,3,1,1,2,1,2,1,2,1,3,1,2,1,2,1,3,1,2,1,3,1,2,
%T 1,3,1,2,6,1,1,2,1,2,1,2,1,3,1,2,1,2,1,3,1,2,1,3,1,2,1,3,1,2,6,1,2,1,
%U 2,1,3,1,2,1,3,1,2,1,3,1,2,6,1,2,1,3,1,2,1,3,1,2,6,1,2,1,3,1,2,6,1,2,1,3,1
%N A bi-Rauzy Pisot substitution.
%F 1->{1, 2} 2->{1, 3} 3->6 4->{4, 5} 5->{4, 6) 6->1
%t s[1] = {1, 2}; s[2] = {1, 3}; s[3] = {6}; s[4] = {4, 5}; s[5] = {4, 6}; s[6] = {1}; t[a_] := Join[a, Flatten[s /@ a]]; p[0] = {1}; p[1] = t[{1}]; p[n_] := t[p[n - 1]] aa = Flatten[Table[p[n], {n, 0, 6}]]
%Y Cf. A073058, A103684.
%K nonn,uned
%O 0,4
%A _Roger L. Bagula_, Apr 07 2005
|