%I #16 Mar 12 2022 11:33:26
%S 1,1,2,7,21,78,305,1304,6007,29854,159012,904986,5479078,35150263,
%T 238033523,1695554145,12663533586,98881246850,805128085616,
%U 6820302066048,59983405937707,546690232627480,5154757226832625,50208266917662433,504482106565647708
%N Weigh transform of the Bell numbers (A000110).
%H Alois P. Heinz, <a href="/A305850/b305850.txt">Table of n, a(n) for n = 0..576</a>
%F G.f.: Product_{k>=1} (1+x^k)^Bell(k).
%p g:= proc(n) option remember; `if`(n=0, 1,
%p add(binomial(n-1, j-1)*g(n-j), j=1..n))
%p end:
%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
%p add(binomial(g(i), j)*b(n-i*j, i-1), j=0..n/i)))
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..30);
%t g[n_] := g[n] = If[n == 0, 1,
%t Sum[Binomial[n - 1, j - 1]*g[n - j], {j, 1, n}]];
%t b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0,
%t Sum[Binomial[g[i], j]*b[n - i*j, i - 1], {j, 0, n/i}]]];
%t a[n_] := b[n, n];
%t Table[a[n], {n, 0, 30}] (* _Jean-François Alcover_, Mar 12 2022, after _Alois P. Heinz_ *)
%Y Cf. A000110, A290351, A305846, A305852.
%K nonn
%O 0,3
%A _Alois P. Heinz_, Jun 11 2018