A052318 Number of labeled rooted trimmed trees with n nodes.
1, 2, 3, 16, 145, 1536, 19579, 290816, 4942305, 94689280, 2020278931, 47523053568, 1222147737265, 34117226135552, 1027550555918475, 33213871550365696, 1146891651823112641, 42135941698113503232, 1641164216596258397347, 67550839668807638712320
Offset: 1
Links
Programs
-
Maple
A:= proc(n) option remember; if n<=1 then x else convert(series(x* exp(A(n-1)-x^2), x,n), polynom) fi end: a:= n-> coeff(A(n+1), x,n)*n!: seq(a(n), n=1..25); # Alois P. Heinz, Aug 23 2008
-
Mathematica
a[n_] := Sum[ Boole[ EvenQ[n-m]]*(m^((n+m)/2-2)/((n-m)/2)!)*((-1)^((n-m)/2)/(m-1)!), {m, 1, n}]*n!; Table[a[n], {n, 1, 18}] (* Jean-François Alcover, Sep 10 2012, after Vladimir Kruchinin *) Rest[CoefficientList[Series[-LambertW[-x/E^(x^2)],{x,0,20}],x]*Range[0,20]!] (* Vaclav Kotesovec, Jan 08 2014 *)
-
Maxima
a(n):=sum((if mod(n-m,2)=0 then m^((n+m)/2-2)/((n-m)/2)!*(-1)^((n-m)/2) else 0)/(m-1)!,m,1,n); /* Vladimir Kruchinin, Aug 07 2012 */
Formula
E.g.f. satisfies A(x) = x*exp(A(x) - x^2).
E.g.f.: -LambertW(-x/exp(x^2)). - Vaclav Kotesovec, Jan 08 2014
a(n) ~ sqrt(1 + LambertW(-2*exp(-2))) * 2^(n/2) * n^(n-1) / (exp(n) * (-LambertW(-2*exp(-2)))^(n/2)). - Vaclav Kotesovec, Jan 08 2014
Comments