OFFSET
1,2
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
a(24) = 1536, since 24 = 2^3 * 3^1 and 1536 = 2^(3^2) * 3^(1^2).
MAPLE
a:= n-> mul(i[1]^(i[2]^2), i=ifactors(n)[2]):
seq(a(n), n=1..50); # Alois P. Heinz, Jun 09 2014
MATHEMATICA
f[p_, e_] := p^(e^2); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Aug 31 2023 *)
PROG
(Haskell)
a054496 n = product $
zipWith (^) (a027748_row n) (map a000290 $ a124010_row n)
-- Reinhard Zumkeller, Apr 27 2013
(PARI) a(n) = my(f=factor(n)); for (i=1, #f~, f[i, 2] = f[i, 2]^2); factorback(f); \\ Michel Marcus, Jun 09 2014
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
Leroy Quet, May 14 2000
STATUS
approved