OFFSET
1,3
COMMENTS
Sequence is defined for all n, as A075423(n) < n.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
PROG
(PARI) rad(n)=vecprod(factor(n)[, 1])
a(n)=my(k); while(n>1, n=rad(n)-1; k++); k \\ Charles R Greathouse IV, Jul 09 2013
(Haskell)
a075425 n = snd $ until ((== 1) . fst)
(\(x, i) -> (a075423 x, i + 1)) (n, 0)
-- Reinhard Zumkeller, Aug 14 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Sep 15 2002
STATUS
approved