A271841 Number of set partitions of [2n] having exactly n pairs (m,m+1) such that m is in some block b and m+1 is in block b+1.

1, 1, 6, 61, 891, 17081, 404275, 11364373, 368982178, 13564841773, 556179920807, 25136678260282, 1240530238800284, 66339010440041817, 3819462133549622416, 235473674234358044731, 15472450628591543437233, 1079168872840695090981865, 79613621745613390178188361

Offset

0,3

Links

Alois P. Heinz, Table of n, a(n) for n = 0..100

Wikipedia, Partition of a set

Formula

a(n) = A185982(2n,n).

Maple

b:= proc(n, i, m, k) option remember; `if`(k>n, 0, `if`(n=0, 1,
       add(`if`(j=i+1 and k=0, 0, b(n-1, j, max(m, j), k-
      `if`(j=i+1, 1, 0))), j=1..m+1)))
    end:
a:= n-> b(2*n, 1, 0, n):
seq(a(n), n=0..18);

Mathematica

b[n_, i_, m_, k_] := b[n, i, m, k] = If[k > n, 0, If[n == 0, 1, Sum[If[j == i + 1 && k == 0, 0, b[n - 1, j, Max[m, j], k - If[j == i + 1, 1, 0]]], {j, 1, m + 1}]]];
a[n_] := b[2*n, 1, 0, n];
Table[a[n], {n, 0, 18}] (* Jean-François Alcover, May 27 2018, translated from Maple *)

Crossrefs

Cf. A185982.

Sequence in context: A086403 A049120 A346983 * A361526 A056546 A127695

Adjacent sequences: A271838 A271839 A271840 * A271842 A271843 A271844

Keyword

nonn

Author

Alois P. Heinz, Apr 15 2016

Status

approved