|
%I #14 Jun 12 2023 11:35:17
%S 1,1,1,3,6,20,55,223,761,3595,14532,77818,361605,2155525,11274781,
%T 73822175,428004750,3046519516,19348533739,148493347507,1023481273549,
%U 8412534272415,62450994058052,546699337652602,4343869829492281,40308548641909593,340994681344324137
%N Number of partitions of [n] such that in each block the smallest element has the same parity as the largest element.
%H Alois P. Heinz, <a href="/A363592/b363592.txt">Table of n, a(n) for n = 0..150</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Partition_of_a_set">Partition of a set</a>
%F a(n) mod 2 = A131719(n+1).
%e a(0) = 1: () the empty partition.
%e a(1) = 1: 1.
%e a(2) = 1: 1|2.
%e a(3) = 3: 123, 13|2, 1|2|3.
%e a(4) = 6: 123|4, 13|24, 13|2|4, 1|234, 1|24|3, 1|2|3|4.
%e a(5) = 20: 12345, 1235|4, 123|4|5, 1245|3, 125|3|4, 1345|2, 135|24, 13|24|5, 135|2|4, 13|2|4|5, 15|234, 1|234|5, 145|2|3, 15|24|3, 1|24|35, 1|24|3|5, 1|2|345, 15|2|3|4, 1|2|35|4, 1|2|3|4|5.
%p b:= proc(n, x, y, u, v) option remember; `if`(y+u>n, 0, `if`(n=0, 1,
%p `if`(y=0, 0, b(n-1, v, u, y-1, x+1)*y)+b(n-1, v, u, y, x+1)+
%p `if`(v=0, 0, b(n-1, v-1, u+1, y, x)*v)+b(n-1, v, u, y, x)*(u+x)))
%p end:
%p a:= n-> b(n, 0$4):
%p seq(a(n), n=0..30);
%Y Cf. A000110, A131719, A361084.
%K nonn
%O 0,4
%A _Alois P. Heinz_, Jun 10 2023
|
