%I #11 Nov 18 2020 06:50:50
%S 0,1,3,10,40,166,749,3683,19275,107806,640970,4024912,26653653,
%T 185401581,1350624721,10282222002,81592209580,673535269054,
%U 5773214891137,51291776763863,471617190143567,4481375500319334,43947651280912186,444258975094335440
%N Number of partitions of n-set with at least one even block.
%H Alois P. Heinz, <a href="/A089005/b089005.txt">Table of n, a(n) for n = 1..500</a>
%F E.g.f.: exp(sinh(x))*(exp(cosh(x)-1)-1).
%p with(combinat):
%p b:= proc(n, i, t) option remember; `if`(n=0, t, `if`(i<1, 0,
%p add(multinomial(n, n-i*j, i$j)/j!*b(n-i*j, i-1,
%p max(t, `if`(j=0, 0, 1-irem(i, 2)))), j=0..n/i)))
%p end:
%p a:= n-> b(n$2, 0):
%p seq(a(n), n=1..30);
%t multinomial[n_, k_List] := n!/Times @@ (k!);
%t b[n_, i_, t_] := b[n, i, t] = If[n == 0, t, If[i<1, 0, Sum[multinomial[n, {n - i j} ~Join~ Table[i, {j}]]/j! b[n - i j, i - 1, Max[t, If[j == 0, 0, 1 - Mod[i, 2]]]], {j, 0, n/i}]]];
%t a[n_] := b[n, n, 0];
%t Array[a, 30] (* _Jean-François Alcover_, Nov 18 2020, after Maple *)
%o (PARI) my(x='x+O('x^30)); concat(0, Vec(serlaplace(exp(sinh(x))*(exp(cosh(x)-1)-1))))
%Y Cf. A000110, A005046, A003724, A059838, A047967.
%K nonn
%O 1,3
%A _Vladeta Jovovic_, Nov 02 2003
|