A036040 revisited

wouter meeussen wouter.meeussen at pandora.be
Sun Mar 23 14:41:00 CET 2003


rule1: Abramowitz & Stegun is always right
rule 2:  even in case of FUBAR, rule 1 applies.

A036040 reads:
1,
1,1,
1,3,1,
1,4,3,6,1,
1,5,10,10,15,10,1,
1,6,15,10,15,60,15,20,45,15,1
           ** **      ** **
1,6,15,15,10,60,20,15,45,15,1 is what Mathematica gives.

Mathematica orders the partitions of 6 as
{6},{5,1},{4,2},{4,1,1},{3,3},{3,2,1},{3,1,1,1},{2,2,2},{2,2,1,1},{2,1,1,1,1
},{1,1,1,1,1,1}
                            ****     **                  ******    ****
The A036040 as it now stands is given by :
    runs[li:{__Integer}]:=((Length/@ Split[#]))&[Sort@ li]
    Table[temp=Map[Reverse,Sort@ (Sort/@ Partitions[w]),{1}];
    Apply[Multinomial,temp,{1}]/Apply[Times,(runs/@ temp)!,{1}],{w,6}]

while the 'natural' (??) Mathematica ordering is produced by
    runs[li:{__Integer}]:=((Length/@ Split[#]))&[Sort@ li]
    Table[Apply[Multinomial,Partitions[w],{1}]/Apply[Times,(runs/@
Partitions[w])!,{1}],{w,6}]

Maybe both should be in the OEIS?

An alternate description would be:
T[n,m]=count of set partitions of n with block lengths given by the m-th
partition of n.

The current description "Triangle of multinomial coefficients" gives the
impression that a calculation
is possible without generating the partitions of n. Did I look over it ?

fumblingly yours,

Wouter.







More information about the SeqFan mailing list