OFFSET
1,3
COMMENTS
FORMULA
a(prime(n)) = A000041(n).
a(x * y) <= a(x) * a(y).
EXAMPLE
The a(30) = 5 partitions are (321), (2211), (3111), (21111), (111111), with corresponding Heinz numbers: 30, 36, 40, 48, 64.
MATHEMATICA
primeMS[n_]:=If[n===1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Table[Length[Union[Sort/@Join@@@Tuples[IntegerPartitions/@primeMS[n]]]], {n, 50}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 04 2018
STATUS
approved