A229413 - OEIS

%I #9 Dec 10 2020 17:38:35

%S 1,1,76,12644,3305017,1245131903,654277037674,467728049807348,

%T 443694809361207824,544852927413901502514,846359710104516310431744,

%U 1629392161877794034658847500,3819592516111353522143561652540,10738740219595085951726635839975852

%N Number of set partitions of {1,...,3n} into sets of size at most n.

%H Alois P. Heinz, <a href="/A229413/b229413.txt">Table of n, a(n) for n = 0..150</a>

%F a(n) = (3n)! * [x^(3n)] exp(Sum_{j=1..n} x^j/j!).

%F a(n) = A229223(3n,n).

%p G:= proc(n, k) option remember; local j; if k>n then G(n, n)

%p       elif n=0 then 1 elif k<1 then 0 else G(n-k, k);

%p       for j from k-1 to 1 by -1 do %*(n-j)/j +G(n-j, k) od; % fi

%p     end:

%p a:= n-> G(3*n, n):

%p seq(a(n), n=0..20);

%t G[n_, k_] := G[n, k] = Module[{j, g}, Which[k > n, G[n, n], n == 0, 1, k < 1, 0, True, g = G[n - k, k]; For[j = k - 1, j >= 1, j--, g = g(n-j)/j + G[n - j, k]]; g]];

%t a[n_] := G[3n, n];

%t a /@ Range[0, 20] (* _Jean-François Alcover_, Dec 10 2020, after _Alois P. Heinz_ *)

%Y Column k=3 of A229243.

%Y Cf. A229223.

%K nonn

%O 0,3

%A _Alois P. Heinz_, Sep 22 2013