OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
FORMULA
If n = Product (p_j^k_j) then a(n) = Product A000041(p_j^k_j).
MAPLE
a:= n-> mul(combinat[numbpart](i[1]^i[2]), i=ifactors(n)[2]):
seq(a(n), n=1..60); # Alois P. Heinz, Apr 11 2025
MATHEMATICA
Table[Times @@ (PartitionsP[#[[1]]^#[[2]]] & /@ FactorInteger[n]), {n, 1, 60}]
PROG
(PARI) a(n) = my(f=factor(n)); prod(i=1, #f~, numbpart(f[i, 1]^f[i, 2])); \\ Michel Marcus, Apr 11 2025
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Ilya Gutkovskiy, Apr 10 2025
STATUS
approved
