OFFSET
1,2
FORMULA
a(p) = p for p = primes (A000040).
a(n) = n^c(n) where c(n) only depends on the prime signature of n. - David A. Corneth, May 05 2020
EXAMPLE
For n = 6; divisors d of 6: {1, 2, 3, 6}; pod(d): {1, 2, 3, 36}; lcm_{d|6} pod(d) = pod(6) = 36; a(6) = 36/1 * 36/2 * 36/3 * 36/36 = 7776.
MATHEMATICA
pod[n_] := Times @@ Divisors[n]; a[n_] := pod[n]^Length[(d = Divisors[n])]/Times @@ (pod /@ d); Array[a, 30] (* Amiram Eldar, May 03 2020 *)
PROG
(Magma) [&*[ LCM([&*Divisors(d): d in Divisors(n)]) / &*Divisors(d): d in Divisors(n)]: n in [1..100]]
(PARI) pod(n) = vecprod(divisors(n));
a(n) = my(d=divisors(n), podn = pod(n)); prod(k=1, #d, podn/pod(d[k])); \\ Michel Marcus, May 03-11 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, May 03 2020
STATUS
approved