OFFSET
1,4
COMMENTS
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
The only fixed point is 1 because otherwise the sequence decreases omega (A001222) by one.
EXAMPLE
The partition (3,2,2,1) with Heinz number 90 has differences plus one (2,1,2) with Heinz number 18, so a(90) = 18.
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
db[n_]:=Times@@Prime/@(1+Differences[primeMS[n]]);
Table[db[n], {n, 100}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Apr 23 2019
STATUS
approved