OFFSET
1,2
MAPLE
a:= n-> (l-> mul(ithprime(i+1)^l[i][2], i=1..nops(l)))
(sort(ifactors(n)[2], (x, y)->x[2]>y[2])):
seq(a(n), n=1..100); # Alois P. Heinz, Jan 12 2023
PROG
(PARI) a(n) = { my(f=vecsort(factor(n)[, 2], , 4)); prod(i=1, #f, prime(i+1)^f[i]) } \\ Andrew Howroyd, Jan 12 2023
(Python)
from math import prod
from sympy import prime, factorint
def A359600(n): return prod(prime(i)**e for i, e in enumerate(sorted(factorint(n).values(), reverse=True), 2)) # Chai Wah Wu, Jan 12 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 12 2023
STATUS
approved