Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A099880
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A099880 Number of preferential arrangements (or simple hierarchies) of 2*n labeled elements with two kinds of elements (where each kind has n elements). +0
1
0, 2, 18, 260, 5250, 136332, 4327092, 162309576, 7024896450, 344582629820, 18890850749628, 1144656941236536, 75963981061424820, 5479642938171428600, 426894499408073653800, 35720957482170932284560, 3195135789350678836128450 (list; graph; listen)
OFFSET

0,2

COMMENT

The unlabeled case seems to be given by A003480, which can be generated by the following combstruct command: SeqUnionU := [S, {S=Sequence(Set(U,card>=1),card>=1), U=Union(a,b),a=Atom, b=Atom},unlabeled]; [seq(count(SeqUnionU, size=n), n=0..20)];

Apart from the first term these numbers appear to show up in the middle of the even-numbered rows in triangle A154921. [From Mats Granvik (mats.granvik(AT)abo.fi), Feb 07 2009]

LINKS

Thomas Wieder, Home Page.

Thomas Wieder, (Old) Home Page.

FORMULA

a(n) = binomial(2*n, n) * Sum_{k=1..n} k! * Stirling2(n, k) = binomial(2*n, n)*A000670(n).

Note: a_0(n) = Sum_{k=0..n} k! * Stirling2(n, k) * binomial(2*n, n) gives 1, 2, 18, 260, 5250, 136332, 4327092, 162309576, ...

a(0) = 1, a(n) = Sum_i=1^n Sum_k=1^i ((2*i)!/i!)*(1/k!)*((i-k)!/(2*(i-k))!) * a(i-k).

EXAMPLE

Let a[1], a[2],...,a[n] and b[1],b[2],...,b[n] denote two kinds "a" and "b" of labeled elements where each kind as n elements in total.

Let : denote a level, e.g. if the elements a[1]) and a[2] are on level l=1 and the element b[1] is on level l=2 then a[1]a[2]:b[1] is a preferrential arrangement (a simple hierarchy) with two levels.

Then for n=2 we have a(2) = 18 arrangements:

a[1]a[2]; a[1]:a[2]; a[2]:a[1]; a[1]b[1]; a[1]:b[1]; b[1]:a[1]; a[1]b[2];

a[1]:b[2]; b[2]:a[1]; a[2]b[1]; a[2]:b[1]; b[1]:a[2]; a[2]b[2]; a[2]:b[2];

b[2]:a[2]; b[1]b[2]; b[1]:b[2]; b[2]:b[1].

MAPLE

main := proc(n::integer) local a, level, Summe; with(combinat): Summe := 0; for level from 1 to n do Summe := Summe + level! * stirling2(n, level) * binomial(2*n, n); ; end do; a := Summe; print("n, a(n):", n, a); end proc;

A099880 := proc(n::integer) local i, k, a, Summe; with(combinat): a[0] := 1; for i from 1 to n do Summe := 0; for k from 1 to i do Summe := Summe + ((2*i)!/i!)*(1/k!)*((i-k)!/(2*(i-k))!) * a[i-k]; end do; a[i] := Summe; end do; a[n] := Summe; print("A099880(n) = ", a[n]); end proc;

a:=n->sum(binomial(2*n, n)*(stirling2(n, k))*k!, k=1..n):seq(a(n), n=0..16); - Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Oct 19 2006

MATHEMATICA

f[n_] := Sum[l! StirlingS2[n, l] Binomial[2n, n], {l, n}]; Table[ f[n], {n, 0, 16}] (from Robert G. Wilson v Nov 04 2004)

CROSSREFS

Cf. A000670, A075729, A003480.

Sequence in context: A109517 A143138 A151362 this_sequence A141009 A143154 A032037

Adjacent sequences: A099877 A099878 A099879 this_sequence A099881 A099882 A099883

KEYWORD

nonn

AUTHOR

Thomas Wieder (wieder.thomas(AT)t-online.de), Nov 02 2004

EXTENSIONS

More terms from Robert G. Wilson v (rgwv(AT)rgwv.com), Nov 04 2004

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 25 20:09 EST 2009. Contains 167514 sequences.


AT&T Labs Research