OFFSET
1,2
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = A066636(n) * n. - Amiram Eldar, Sep 15 2023
EXAMPLE
a(40)=10^3 > 40 > 10=rad(40).
MATHEMATICA
rad[n_] := Times @@ First /@ FactorInteger[n]; mpe[n_] := Max @@ Last /@ FactorInteger[n]; mpe[1] = 0; a[n_] := rad[n]^mpe[n]; Table[a[n], {n, 1, 59}] (* Jean-François Alcover, Mar 27 2013 *)
PROG
(PARI) a(n)=if(n==1, return(1)); my(f=factor(n)); prod(i=1, #f~, f[i, 1])^ vecmax(f[, 2]) \\ Charles R Greathouse IV, Aug 21 2013
(Haskell)
a066638 n = a007947 n ^ a051903 n -- Reinhard Zumkeller, Jun 17 2015
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Reinhard Zumkeller, Jan 09 2002
STATUS
approved