OFFSET
1,2
COMMENTS
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The z-density of a multiset S of positive integers is Sum_{s in S} (omega(s) - 1) - omega(lcm(S)) where omega = A001221 is number of distinct prime factors.
MATHEMATICA
zens[n_]:=If[n==1, 0, Total@Cases[FactorInteger[n], {p_, k_}:>k*(PrimeNu[PrimePi[p]]-1)]-PrimeNu[LCM@@Cases[FactorInteger[n], {p_, k_}:>PrimePi[p]]]];
Select[Range[1000], zens[#]==0&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, May 24 2018
STATUS
approved