OFFSET
1,1
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = rad(n)*rad(n+1)*rad(n+2) if n is odd; or rad(n/2)*rad(n+1)*rad(n/2+1) if n is even. - C. Ronaldo (aga_new_ac(AT)hotmail.com), Dec 13 2004
EXAMPLE
a(3) = rad(3*4*5) = 30.
MAPLE
with(numtheory):rad:=proc(n) local s, i: s:=ifactors(n)[2]: RETURN(mul(s[i][1], i=1..nops(s))): end; seq(rad(n*(n+1)*(n+2)), n=1..60); seq(piecewise(n mod 2=0, rad(n/2)*rad(n+1)*rad(n/2+1), rad(n)*rad(n+1)*rad(n+2)), n=1..60); (C. Ronaldo)
MATHEMATICA
lsf[n_]:=Max[Select[Divisors[n], SquareFreeQ]]; lsf/@Table[n(n+1)(n+2), {n, 50}] (* Harvey P. Dale, Oct 18 2020 *)
a[n_] := Times @@ Union @@ (FactorInteger[#][[;; , 1]] & /@ (n + {0, 1, 2})); Array[a, 50] (* Amiram Eldar, Jun 30 2022 *)
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))", "))
(Haskell)
a078637 n = a007947 $ product [n..n+2] -- Reinhard Zumkeller, Jul 04 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Jon Perry, Dec 12 2002
STATUS
approved