OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
PROG
(PARI) rad(n)=my(f=factor(n)[, 1]); prod(i=1, #f, f[i])
a(n)=my(r=rad(n)); r+n/r \\ Charles R Greathouse IV, Aug 08 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Nov 06 2002
STATUS
approved