login
Weigh transform of the Bell numbers (A000110).
4

%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