|
%I #9 May 27 2018 10:33:45
%S 1,1,2,5,15,52,203,877,4140,21147,115975,678570,4213596,27643917,
%T 190873448,1382311874,10468191626,82676728870,679376431959,
%U 5796059673820,51241781683988,468639956185787,4426928995172738,43131913020157751,432877526792148982
%N Number of set partitions of [n] such that j is member of block b only if b = 1 or at least one of j-1, ..., j-9 is member of a block >= b-1.
%H Alois P. Heinz, <a href="/A287672/b287672.txt">Table of n, a(n) for n = 0..30</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Partition_of_a_set">Partition of a set</a>
%F a(n) = A287641(n,9).
%F a(n) = A000110(n) for n <= 11.
%e a(12) = 4213596 = 4213597 - 1 = A000110(12) - 1 counts all set partitions of [12] except: 13456789(10)(11)|2|(12).
%p b:= proc(n, l) option remember; `if`(n=0, 1, add(b(n-1,
%p [seq(max(l[i], j), i=2..nops(l)), j]), j=1..l[1]+1))
%p end:
%p a:= n-> b(n, [0$9]):
%p seq(a(n), n=0..20);
%t b[n_, l_] := b[n, l] = If[n == 0, 1, Sum[b[n - 1, Append[Table[Max[l[[i]], j], {i, 2, Length[l]}], j]], {j, 1, l[[1]] + 1}]];
%t a[n_] := b[n, Table[0, 9]];
%t Table[a[n], {n, 0, 20}] (* _Jean-François Alcover_, May 27 2018, from Maple *)
%Y Column k=9 of A287641.
%Y Cf. A000110.
%K nonn
%O 0,3
%A _Alois P. Heinz_, May 29 2017
|
