OFFSET
1,4
COMMENTS
a(1) = 1 (empty product).
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
EXAMPLE
45 = 3^2 * 5^1. So a(45) = 2^2 * 1^1 = 4.
72 = 2^3 * 3^2. So a(72) = 3^3 * 2^2 = 108.
MATHEMATICA
f[n_] := Block[{fi = Last@Transpose@FactorInteger@n}, Times @@ (fi^fi)]; Rest@Array[f, 93] (* Robert G. Wilson v *)
PROG
(PARI) a(n)=local(v, r, i); v=factorint(n); r=1; for(i=1, matsize(v)[1], r*=v[i, 2]^v[i, 2]); r (Herrgesell)
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Leroy Quet, Dec 25 2005
EXTENSIONS
More terms from Robert G. Wilson v and Lambert Herrgesell (zero815(AT)googlemail.com), Dec 27 2005
Corrected the starting offset, data section extended to 105 terms - Antti Karttunen, May 28 2017
STATUS
approved