OFFSET
1,2
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.
Also numbers of which it is possible to choose a different constant integer partition of each prime index. The multiset of prime indices of n is row n of A112798.
EXAMPLE
The prime indices of 60 are {1,1,2,3}, and we have the following 4 choices of constant partitions:
((1),(1),(2),(3))
((1),(1),(2),(1,1,1))
((1),(1),(1,1),(3))
((1),(1),(1,1),(1,1,1))
Since none of these is strict, 60 is not in the sequence.
The prime indices of 90 are {1,2,2,3}, and we have the following 4 strict choices:
((1),(2),(1,1),(3))
((1),(2),(1,1),(1,1,1))
((1),(1,1),(2),(3))
((1),(1,1),(2),(1,1,1))
So 90 is in the sequence.
MATHEMATICA
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[100], Select[Tuples[Select[IntegerPartitions[#], SameQ@@#&]&/@prix[#]], UnsameQ@@#&]!={}&]
- or -
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[100], And@@Table[Count[prix[#], k]<=DivisorSigma[0, k], {k, Union[prix[#]]}]&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Sep 04 2025
STATUS
approved
