OFFSET
1,2
COMMENTS
Or replace prime(i) in n by prime(n*i).
All terms are odd.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 1 because 1 is the empty product.
a(2) = 3 = prime(2) = prime(2*1) because 2 = prime(1).
a(3) = 13 = prime(6) = prime(3*2) because 3 = prime(2).
a(4) = 49 = 7^2 = prime(4)^2 = prime(4*1)^2 because 4 = prime(1)^2.
MAPLE
a:= n-> mul(ithprime(n*numtheory[pi](i[1]))^i[2], i=ifactors(n)[2]):
seq(a(n), n=1..45);
PROG
(PARI) a(n) = my(f=factor(n)); for (k=1, #f~, f[k, 1] = prime(n*primepi(f[k, 1]))); factorback(f); \\ Michel Marcus, Mar 02 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Mar 01 2022
STATUS
approved