1,1

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000

(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

a(n)=A007947(n)+A003557(n).

Sequence in context: A193768 A205791 A039696 * A245092 A092601 A162906

Adjacent sequences: A076329 A076330 A076331 * A076333 A076334 A076335

nonn

Reinhard Zumkeller, Nov 06 2002

approved