OFFSET
1,3
COMMENTS
A020754 probably gives the records in this sequence. - Charles R Greathouse IV, Aug 02 2013
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
MATHEMATICA
rad[n_] := Times @@ FactorInteger[n][[All, 1]]; a[1] = 1; a[n_] := For[k = 1, True, k++, If[ rad[n+k] > rad[n], Return[k]]]; Table[ a[n], {n, 1, 105}] (* Jean-François Alcover, Oct 05 2012 *)
PROG
(Haskell)
a079229 n = a079229_list !! (n-1)
a079229_list = f a007947_list where
f (x:xs) = ((+ 1) . length $ takeWhile (<= x) xs) : f xs
-- Reinhard Zumkeller, Oct 06 2012
(PARI) rad(n)=my(f=factor(n)[, 1]); prod(i=1, #f, f[i])
a(n)=my(r=rad(n), k); while(r>=rad(n+k++), ); k \\ Charles R Greathouse IV, Aug 02 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jan 02 2003
STATUS
approved