login
rad(n(n+1)(n+2)(n+3)).
1

%I #5 Jun 20 2013 18:51:19

%S 6,30,30,210,210,42,210,330,330,4290,6006,2730,2730,3570,510,1938,

%T 9690,3990,43890,53130,10626,7590,8970,390,2730,15834,6090,188790,

%U 26970,10230,34782,39270,39270,132090,147630,54834,274170,303810,111930

%N rad(n(n+1)(n+2)(n+3)).

%e a(3)=rad(3*4*5*6)=30.

%o (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))","))

%o (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

%Y Cf. A007947.

%K nonn

%O 1,1

%A _Jon Perry_, Dec 12 2002