login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A305197 Number of set partitions of [n] with symmetric block size list of length A004525(n). 4

%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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 20:05 EDT 2024. Contains 371254 sequences. (Running on oeis4.)