login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A138056 Levels of substitution A103684 (based on the morphism f: 1->{1,2}, 2->{1,3}, 3->{3}) like Markov substitution taken as polynomials p(x,n)]and coefficients of the differential polynomials returned as q(x,n) =dp(x,n)dx coefficients (first zero omitted). 0

%I #4 Sep 09 2023 19:36:39

%S 2,2,2,9,2,2,9,4,10,6,2,2,9,4,10,6,7,16,9,30,11,24,2,2,9,4,10,6,7,16,

%T 9,30,11,24,13,28,15,48,17,36,19,20,42,22,69,2,2,9,4,10,6,7,16,9,30,

%U 11,24,13,28,15,48,17,36,19,20,42,22,69,24,50,26,81,28,58,30,31,64,33,102,35

%N Levels of substitution A103684 (based on the morphism f: 1->{1,2}, 2->{1,3}, 3->{3}) like Markov substitution taken as polynomials p(x,n)]and coefficients of the differential polynomials returned as q(x,n) =dp(x,n)dx coefficients (first zero omitted).

%C Row sums with zero: {0, 2, 13, 33, 130, 459, 1533, 5266, 17884, 60532, 205129, ...};

%C This sequence uses the Markov substitution form that I have been using in my chord-geometry/graph sequences.

%C This method of differentiating a substitution appears to be new.

%F f: 1->{1,2}, 2->{1,3}, 3->{3}); p(x,n)=Sum[Substitution[n,m]*t(m-1),{m,1,n}]; q(x,n)=dp(x,n)dx; out_n,m=Coefficients(q(x,n).

%e {2},

%e {2, 2, 9},

%e {2, 2, 9, 4, 10, 6},

%e {2, 2, 9, 4, 10, 6, 7, 16, 9, 30, 11, 24},

%e {2, 2, 9, 4, 10, 6, 7, 16, 9, 30, 11, 24, 13, 28, 15, 48, 17, 36, 19, 20, 42, 22, 69},

%e {2, 2, 9, 4, 10, 6, 7, 16, 9, 30, 11, 24, 13, 28, 15, 48, 17, 36, 19, 20, 42, 22, 69, 24, 50, 26, 81, 28, 58, 30, 31, 64, 33, 102, 35, 72, 37, 76, 39, 120, 41, 84, 43}

%t Clear[a, s, p, t, m, n]; (* substitution *); s[1] = {1, 2}; s[2] = {1, 3}; s[3] = {1}; t[a_] := Flatten[s /@ a]; p[0] = {1}; p[1] = t[p[0]]; p[n_] := t[p[n - 1]]; (*A103684*); a = Table[p[n], {n, 0, 10}]; Flatten[a]; b = Table[CoefficientList[D[Apply[Plus, Table[a[[n]][[m]]*x^( m - 1), {m, 1, Length[a[[n]]]}]], x], x], {n, 1, 11}]; Flatten[b] Table[Apply[Plus, CoefficientList[D[Apply[Plus, Table[a[[n]][[m]]* x^(m - 1), {m, 1, Length[a[[n]]]}]], x], x]], {n, 1, 11}];

%Y Cf. A103684.

%K nonn,uned,tabf

%O 1,1

%A _Roger L. Bagula_, May 02 2008

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 07:33 EDT 2024. Contains 371235 sequences. (Running on oeis4.)