OFFSET
1,2
COMMENTS
rad(n) is the product of the primes dividing n (A007947).
LINKS
T. D. Noe, Table of n, a(n) for n = 1..100
FORMULA
a(n)= n^A007947(n).
EXAMPLE
a(1) = 1^1 = 1 ; a(2) = 2^2 = 4 ; a(3) = 3^3 = 27 ; a(4) = 4^2 = 16.
MAPLE
with(numtheory):for n from 1 to 200 do:t1:= ifactors(n)[2] : p :=mul(t1[i][1], i=1..nops(t1)):x:=n^p:print(x):od:
MATHEMATICA
Table[n^(Times @@ Transpose[FactorInteger[n]][[1]]), {n, 20}]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michel Lagneau, Apr 06 2010
STATUS
approved