|
%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
|
