OFFSET
1,4
FORMULA
a(n) = n for n = 1 and squarefree semiprimes (A006881).
EXAMPLE
For n = 6; a(6) = pod(1)/1 * pod(2)/2 * pod(3)/3 * pod(6)/6 = 1/1 * 2/2 * 3/3 * 36/6 = 6.
MATHEMATICA
Array[Product[Apply[Times, Divisors@ d]/d, {d, Divisors@ #}] &, 51] (* Michael De Vlieger, Jan 19 2019 *)
PROG
(Magma) [&*[&*[c: c in Divisors(d)] / d: d in Divisors(n)]: n in [1..100]]
(PARI) a(n) = my(x=1); fordiv(n, d, x*=vecprod(divisors(d))/d); x; \\ Michel Marcus, Dec 23 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Dec 23 2018
STATUS
approved