OFFSET
1,4
COMMENTS
a(n) = if n=p^k (p prime, k>0) then p^(k-1) else n/(spf(n)*gpf(n)).
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
MATHEMATICA
gc[n_]:=Module[{fi=Transpose[FactorInteger[n]][[1]]}, GCD[n/First[fi], n/Last[ fi]]]; Array[gc, 110] (* Harvey P. Dale, Jun 17 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Aug 07 2002
STATUS
approved