OFFSET
1,12
COMMENTS
By convention, mu() = 0.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
LINKS
FORMULA
EXAMPLE
Heinz number of (2,1,1) is 12, so mu(2,1,1) = a(12) = 2.
MATHEMATICA
nn=120;
ptns=Table[If[n===1, {}, Join@@Cases[FactorInteger[n]//Reverse, {p_, k_}:>Table[PrimePi[p], {k}]]], {n, nn}];
tris=Join@@Map[Tuples[IntegerPartitions/@#]&, ptns];
mu[y_]:=mu[y]=If[Length[y]===1, 1, -Sum[Times@@mu/@t, {t, Select[tris, And[Length[#]>1, Sort[Join@@#, Greater]===y]&]}]];
mu/@ptns
CROSSREFS
KEYWORD
sign
AUTHOR
Gus Wiseman, Feb 05 2018
STATUS
approved