login
Number of forests of different trees.
2

%I #23 Nov 15 2014 17:51:32

%S 1,1,1,2,3,6,12,24,49,105,231,517,1188,2783,6643,16101,39606,98605,

%T 248287,631214,1618878,4183964,10889305,28517954,75111521,198851386,

%U 528929895,1412993746,3789733399,10201625514,27555373561,74664487653,202908119046,552939614498

%N Number of forests of different trees.

%H Alois P. Heinz, <a href="/A035055/b035055.txt">Table of n, a(n) for n = 0..1000</a>

%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>

%F Weigh transform of A000055.

%F a(n) ~ c * d^n / n^(5/2), where d = A051491 = 2.9557652856519949747148175..., c = 0.89246007934060351292465521837... . - _Vaclav Kotesovec_, Aug 25 2014

%p with(numtheory):

%p b:= proc(n) option remember; `if`(n<2, n,

%p (add(add(d*b(d), d=divisors(j))*b(n-j), j=1..n-1))/(n-1))

%p end:

%p h:= proc(n) option remember; `if`(n=0, 1, b(n)-(add(b(k)*b(n-k),

%p k=0..n) -`if`(irem(n, 2)=0, b(n/2), 0))/2)

%p end:

%p g:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,

%p add(binomial(h(i), j)*g(n-i*j, i-1), j=0..n/i)))

%p end:

%p a:= n-> g(n, n):

%p seq(a(n), n=0..40); # _Alois P. Heinz_, May 19 2013

%t nn = 20; t[x_] := Sum[a[n] x^n, {n, 1, nn}]; a[0] = 0;

%t b = Flatten[

%t sol = SolveAlways[

%t 0 == Series[

%t t[x] - x Product[1/(1 - x^i)^ a[i], {i, 1, nn}], {x, 0, nn}],

%t x]; Table[a[n], {n, 0, nn}] /. sol];

%t r[x_] := Sum[b[[n]] x^(n - 1), {n, 1, nn + 1}]; c =

%t Drop[CoefficientList[

%t Series[r[x] - (r[x]^2/2 - r[x^2]/2), {x, 0, nn}], x],

%t 1]; CoefficientList[

%t Series[Product[(1 + x^i)^c[[i]], {i, 1, nn}], {x, 0, nn}], x] (* _Geoffrey Critzer_, Nov 15 2014 *)

%Y Cf. A005195.

%K nonn

%O 0,4

%A _Christian G. Bower_, Oct 15 1998