|
%I #13 Nov 22 2023 01:44:38
%S 1,3,8,14,25,40,62,89,127,173,233,304,393,497,624,769,942,1139,1369,
%T 1627,1925,2257,2635,3053,3524,4042,4621,5253,5954,6717,7557,8466,
%U 9462,10536,11706,12963,14326,15786,17363,19046,20857,22786,24854
%N Number of asymmetric (identity) trees with n nodes and 4 leaves.
%H G. C. Greubel, <a href="/A055335/b055335.txt">Table of n, a(n) for n = 9..1000</a>
%H <a href="/index/Tra#trees">Index entries for sequences related to trees</a>
%H <a href="/index/Rec#order_11">Index entries for linear recurrences with constant coefficients</a>, signature (2,0,-1,0,-2,2,0,1,0,-2,1).
%F G.f.: x^9*(1+x+2*x^2-x^3)/((1-x)^2*(1-x^2)^2*(1+x^2)*(1-x^3)).
%F a(n) = (-225 -762*n +516*n^2 -100*n^3 +6*n^4)/1152 -(3/128)*(-1)^n*(2*n -11) -(1/16)*(2 -(-1)^n)*(-1)^binomial(n,2) -(1/9)*ChebyshevU(n-1, -1/2) + [n=1]. - _G. C. Greubel_, Nov 10 2023
%t Drop[CoefficientList[Series[x^9*(1+x+2*x^2-x^3)/((1-x)*Product[1-x^j, {j,4}]), {x,0,50}], x], 9] (* _G. C. Greubel_, Nov 10 2023 *)
%o (Magma) R<x>:=PowerSeriesRing(Integers(), 50); Coefficients(R!( x^9*(1+x+2*x^2-x^3)/((1-x)*(&*[1-x^j: j in [1..4]])) )); // _G. C. Greubel_, Nov 10 2023
%o (SageMath)
%o def A055335_list(prec):
%o P.<x> = PowerSeriesRing(ZZ, prec)
%o return P( x^9*(1+x+2*x^2-x^3)/((1-x)*product(1-x^j for j in range(1,5))) ).list()
%o a=A055335_list(50); a[9:] # _G. C. Greubel_, Nov 10 2023
%Y Column 4 of A055334.
%K nonn
%O 9,2
%A _Christian G. Bower_, May 12 2000
|
