OFFSET
1,1
EXAMPLE
a(3)=rad(3*4*5*6)=30.
PROG
(PARI) rad(n)=local(p, i); p=factor(n)[, 1]; prod(i=1, length(p), p[i]) for (k=1, 100, print1(rad(k*(k+1)*(k+2)*(k+3))", "))
(PARI) a(n)=my(v=[]); for(i=0, 3, v=concat(v, factor(n+i)[, 1])); v=vecsort(v, , 8); prod(i=1, #v, v[i]) \\ Charles R Greathouse IV, Jun 20 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Jon Perry, Dec 12 2002
STATUS
approved