|
%I #10 Dec 14 2020 05:13:10
%S 1,6,42,330,2970,30240,345240,4377240,61122600,933055200,15470254800,
%T 277005128400,5329454130000,109681187616000,2404894892400000,
%U 55977698400624000,1378748676601296000,35829233832135744000,979763376201049440000,28124715476056399200000
%N Number of ordered set partitions of [n] where the maximal block size equals two.
%H Alois P. Heinz, <a href="/A320758/b320758.txt">Table of n, a(n) for n = 2..428</a>
%F E.g.f.: 1/(1-Sum_{i=1..2} x^i/i!) - 1/(1-x).
%F A(n) = A080599(n) - A000142(n).
%p b:= proc(n, k) option remember; `if`(n=0, 1, add(
%p b(n-i, k)*binomial(n, i), i=1..min(n, k)))
%p end:
%p a:= n-> (k-> b(n, k) -b(n, k-1))(2):
%p seq(a(n), n=2..25);
%t b[n_, k_] := b[n, k] = If[n == 0, 1, Sum[b[n - i, k] Binomial[n, i], {i, 1, Min[n, k]}]];
%t a[n_] := With[{k = 2}, b[n, k] - b[n, k-1]];
%t a /@ Range[2, 25] (* _Jean-François Alcover_, Dec 14 2020, after _Alois P. Heinz_ *)
%Y Column k=2 of A276922.
%Y Cf. A000142, A080599.
%K nonn
%O 2,2
%A _Alois P. Heinz_, Oct 20 2018
|
