OFFSET
1,2
COMMENTS
LINKS
R. Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
a(12)=a(2*6)=a(8)=a(2*4)=a(6)=a(2*3)=a(5)=a(5*1)=5.
MATHEMATICA
lpf[n_] := FactorInteger[n][[1, 1]]; a[n_] := a[n] = If[n <= lpf[n]^2, lpf[n], a[lpf[n] + n/lpf[n]]]; Table[a[n], {n, 1, 74}](* Jean-François Alcover, Dec 21 2011 *)
PROG
(Haskell)
a068319 n = if n <= spf ^ 2 then spf else a068319 $ spf + div n spf
where spf = a020639 n
-- Reinhard Zumkeller, Jun 24 2013
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Reinhard Zumkeller, Feb 27 2002, Jul 13 2007
STATUS
approved