%I #10 May 27 2018 10:33:04

%S 1,1,2,5,15,52,203,876,4116,20827,112538,645045,3900512,24769152,

%T 164546915,1139818861,8209631792,61331709492,474221335902,

%U 3787741281763,31199052157724,264605708064825,2307562757319104,20666169125398768,189855243829576499

%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-4 is member of a block >= b-1.

%H Alois P. Heinz, <a href="/A287667/b287667.txt">Table of n, a(n) for n = 0..90</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Partition_of_a_set">Partition of a set</a>

%F a(n) = A287641(n,4).

%F a(n) = A000110(n) for n <= 6.

%e a(7) = 876 = 877 - 1 = A000110(7) - 1 counts all set partitions of [7] except: 13456|2|7.

%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$4]):

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

%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, 4]];

%t Table[a[n], {n, 0, 26}] (* _Jean-François Alcover_, May 27 2018, from Maple *)

%Y Column k=4 of A287641.

%Y Cf. A000110.

%K nonn

%O 0,3

%A _Alois P. Heinz_, May 29 2017