OFFSET
1,7
COMMENTS
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
EXAMPLE
The a(21) = 5 maximal chains are the rows:
(111111)<(21111)<(2211)<(222)<(42)
(111111)<(21111)<(2211)<(411)<(42)
(111111)<(21111)<(2211)<(321)<(42)
(111111)<(21111)<(3111)<(411)<(42)
(111111)<(21111)<(3111)<(321)<(42)
MATHEMATICA
pcovs[ptn_]:=Select[Union[Reverse/@Sort/@Join@@@Tuples[IntegerPartitions/@ptn]], Length[#]===Length[ptn]+1&];
coc[ptn_]:=coc[ptn]=If[Max[ptn]===1, 1, Total[coc/@pcovs[ptn]]];
primeMS[n_]:=If[n===1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Table[coc[Reverse[primeMS[n]]], {n, 50}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 04 2018
STATUS
approved