OFFSET
1,3
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
MATHEMATICA
a[1]=a[2]=1; a[n_] := a[n] = a[n-1]+a[Divisors[n-1][[-2]]]; Table[a[n], {n, 1, 58}] (* Jean-François Alcover, Mar 20 2011 *)
PROG
(PARI) a(n) = if(n==1, 1, if(n==2, 1, a(n-1)+a((n-1)/divisors(n-1)[2]))) /* j(m) = (largest divisor) = m / (smallest divisor) */ \\ Michael B. Porter, Oct 06 2009
(Haskell)
a080655 n = a080655_list !! (n-1)
a080655_list = 1 : 1 : f 3 where
f n = (a080655 (n-1) + a080655 (a032742 (n-1))) : f (n+1)
-- Reinhard Zumkeller, Aug 21 2011
CROSSREFS
KEYWORD
nonn,easy,nice
AUTHOR
Leroy Quet, Feb 28 2003
EXTENSIONS
More terms from Naohiro Nomoto, Mar 07 2003
STATUS
approved