OFFSET
1,2
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
FORMULA
For n > 1: a(n) = (n + a(n-1)) / GCD(n, a(n-1)). - Reinhard Zumkeller, Nov 26 2015
MAPLE
A131134 := proc(n) option remember ; if n = 1 then 1; else (1/n+1/A131134(n-1))*lcm(n, A131134(n-1)) ; fi ; end: seq(A131134(n), n=1..80) ; # R. J. Mathar, Oct 17 2007
MATHEMATICA
nxt[{n_, a_}]:={n+1, (1/(n+1)+1/a)LCM[n+1, a]}; Transpose[NestList[nxt, {1, 1}, 70]][[2]] (* Harvey P. Dale, Aug 13 2013 *)
PROG
(Haskell)
a131134 n = a131134_list !! (n-1)
a131134_list = 1 : zipWith (\v w -> (v+w) `div` gcd v w) [2..] a131134_list
-- Reinhard Zumkeller, Nov 26 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Jun 17 2007
EXTENSIONS
More terms from R. J. Mathar, Oct 17 2007
STATUS
approved