%I #13 Jun 13 2018 14:11:15
%S 1,1,1,1,3,7,19,56,171,470,2066,10299,31346,91925,559987,3939653,
%T 11954993,36298007,282835456,2571177913,7785919355,24158837489,
%U 229359684137,2557117944391,7731656573016,24350208829581,272633076900991,3601150175699409,10876116332074739
%N Number of set partitions of [n] with symmetric block size list of length A004525(n).
%H Alois P. Heinz, <a href="/A305197/b305197.txt">Table of n, a(n) for n = 0..400</a>
%F a(n) = A275281(n,(n+sin(n*Pi/2))/2).
%p b:= proc(n, s) option remember; expand(`if`(n>s,
%p binomial(n-1, n-s-1)*x, 1)+add(binomial(n-1, j-1)*
%p b(n-j, s+j)*binomial(s+j-1, j-1), j=1..(n-s)/2)*x^2)
%p end:
%p a:= n-> coeff(b(n, 0), x, (n+sin(n*Pi/2))/2):
%p seq(a(n), n=0..30);
%t b[n_, s_] := b[n, s] = Expand[If[n > s, Binomial[n - 1, n - s - 1]*x, 1] + Sum[Binomial[n - 1, j - 1]*b[n - j, s + j]*Binomial[s + j - 1, j - 1], {j, 1, (n - s)/2}]*x^2];
%t a[n_] := Coefficient[b[n, 0], x, (n + Sin[n*Pi/2])/2];
%t Table[a[n], {n, 0, 30}] (* _Jean-François Alcover_, Jun 13 2018, from Maple *)
%Y Bisections give A275283 (even part), A305198 (odd part).
%Y Cf. A004525, A275281.
%K nonn
%O 0,5
%A _Alois P. Heinz_, May 27 2018
|