OFFSET
2,3
COMMENTS
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
LINKS
MAPLE
a:= n-> (l-> add(i, i=l)/igcd(l[]))(map(i->
numtheory[pi](i[1])$i[2], ifactors(n)[2])):
seq(a(n), n=2..100); # Alois P. Heinz, Jul 03 2018
MATHEMATICA
Table[With[{pms=If[n===1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]]}, Total[pms]/GCD@@pms], {n, 2, 100}]
PROG
(PARI) A316436(n) = { my(f = factor(n), pis = apply(p -> primepi(p), f[, 1]~), es = f[, 2]~, g = gcd(pis)); sum(i=1, #f~, pis[i]*es[i])/g; }; \\ Antti Karttunen, Sep 10 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jul 03 2018
EXTENSIONS
More terms from Antti Karttunen, Sep 10 2018
STATUS
approved