OFFSET
1,5
COMMENTS
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
Problem: What is the image? In the nonnegative case it appears to start: 0, 1, 2, 3, 5, 7, 9, ...
FORMULA
a(n) = n - A122111(n).
EXAMPLE
The partition (4,4,1,1) has Heinz number 196 and its conjugate (4,2,2,2) has Heinz number 189, so a(196) = 196 - 189 = 7.
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
conj[y_]:=If[Length[y]==0, y, Table[Length[Select[y, #>=k&]], {k, 1, Max[y]}]];
Table[n-Times@@Prime/@conj[primeMS[n]], {n, 30}]
CROSSREFS
A similar sequence is A175508.
KEYWORD
sign
AUTHOR
Gus Wiseman, Mar 20 2022
STATUS
approved