The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.



(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A116379 Number of ternary rooted identity (distinct subtrees) trees with n nodes. 4


%S 1,1,1,2,3,6,12,25,52,113,245,542,1205,2707,6113,13907,31780,73010,

%T 168399,389991,906231,2112742,4939689,11580640,27216387,64110091,

%U 151334814,357938832,848153045,2013190671,4786210412,11396004660,27172368314,64875527649

%N Number of ternary rooted identity (distinct subtrees) trees with n nodes.

%C It is not known if these trees have the asymptotic form C rho^{-n} n^{-3/2}, whereas the identity binary trees, A063895, do, see the Jason P. Bell et al. reference.

%H Alois P. Heinz, <a href="/A116379/b116379.txt">Table of n, a(n) for n = 1..1000</a>

%H Jason P. Bell, Stanley N. Burris and Karen A. Yeats, <a href="http://www.arXiv.org/abs/math.CO/0512432">Counting Rooted Trees: The Universal Law t(n) ~ C rho^{-n} n^{-3/2}</a>

%F G.f. satisfies A(x) = x(1+A(x)+A(x)^2/2-A(x^2)/2+A(x)^3/6-A(x)A(x^2)/2+A(x^3)/3), that is A(x) = x(1+Set_{<=3}(A)(x))

%p A:= proc(n) option remember; local T; if n<=1 then x else T:= unapply(A(n-1), x); convert(series(x* (1+T(x)+ T(x)^2/2- T(x^2)/2+ T(x)^3/6- T(x)*T(x^2)/2+ T(x^3)/3), x, n+1), polynom) fi end: a:= n-> coeff(A(n), x, n): seq(a(n), n=1..40); # _Alois P. Heinz_, Aug 22 2008

%t A[n_] := A[n] = If[n <= 1, x, T[y_] = A[n-1] /. x -> y; Normal[Series[y*(1+T[y]+T[y]^2/2-T[y^2]/2+T[y]^3/6-T[y]*T[y^2]/2+T[y^3]/3), {y, 0, n+1}]] /. y -> x] ; a[n_] := Coefficient[A[n], x, n]; Table[a[n], {n, 1, 40}] (* _Jean-Fran├žois Alcover_, Feb 13 2014, after Maple *)

%o (C) #include <ginac/ginac.h> using namespace GiNaC; int main(){ int i, order = 40; symbol x("x"); ex T = x; for (i=0; i<order; i++) T = (x+x*(T + pow(T,2)/2 - T.subs(x==pow(x,2))/2 + pow(T,3)/6 - T*T.subs(x==pow(x,2))/2 + T.subs(x==pow(x,3))/3)).series(x, i+2); for (i=1; i<=order; i++) std::cout << T.coeff(x, i) << ", "; }

%Y Cf. A004111, A063895, A116380.

%K nonn

%O 1,4

%A Karen Yeats (kayeats(AT)bu.edu), Feb 06 2006

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 18 21:20 EST 2020. Contains 332028 sequences. (Running on oeis4.)