clearly is clearly
Wouter Meeussen
eu000949 at pophost.eunet.be
Tue Oct 12 00:40:48 CEST 1999
more where that came from:
res=Table[ Length[#]& /@
Split[Sort[Apply[Plus,KSubsets[Range[2n], n],1]-n(n+1)/2]], {n,0,9}]
{{1},
{1,1},
{1,1,2,1,1},
{1,1,2,3,3,3,3,2,1,1},
{1,1,2,3,5,5,7,7,8,7,7,5,5,3,2,1,1},
{1,1,2,3,5,7,9,11,14,16,18,19,20,20,19,...
...
if we could produce the n^2+1 values in the n'th row of "res"
then we could properly continue this sequence:
Length/@(Union /@ res)
{1,1,2,3,6,12,17,24,32,40}
-------------------------------------------------------
we can put Humpty-Dumpty together again:
multi=Table[Range[0,n^2]+n(n+1)/2,{n,0,9}];
Thread[d[res,multi]]/.d->Dot
{0,3,30,210,1260,6930,36036,180180,875160,4157010}
the above is the sum of all the (2n choose n) sums.
I am sorry, but the terms
0,3,30,210,1260,6930,36036,180180,875160,4157010
do not match anything in the table,
because it's too simple:
Table[Binomial[2n,n] n(2n+1)/2,{n,0,9}]
{0,3,30,210,1260,6930,36036,180180,875160,4157010}
and the average outcome is
{0, 3/2, 5, 21/2, 18, 55/2, 39, 105/2, 68, 171/2}
or n(2n+1)/2
squeezing the (2n choose n) trials into the n^2+1
different outcome-categories *should* yield something
familiar, and it does.
Dr. Wouter L. J. MEEUSSEN
w.meeussen.vdmcc at vandemoortele.be
eu000949 at pophost.eunet.be
More information about the SeqFan
mailing list