OFFSET
0,2
COMMENTS
Equivalently, a(0)=1; a(n) = (a(n-1)+n)/gcd(n,a(n-1)). - Ctibor O. Zizka, Feb 15 2008
LINKS
Harvey P. Dale, Table of n, a(n) for n = 0..1000
MAPLE
a[0] := 1; for n to 70 do a[n] := (1/n+1/a[n-1])*lcm(n, a[n-1]) end do; seq(a[n], n = 0 .. 70) # Emeric Deutsch, Jun 21 2007
MATHEMATICA
nxt[{n_, a_}]:={n+1, (1/(n+1)+1/a)LCM[n+1, a]}; Transpose[NestList[nxt, {0, 1}, 70]][[2]] (* Harvey P. Dale, Mar 04 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Jun 17 2007
EXTENSIONS
More terms from Emeric Deutsch, Jun 21 2007
Edited by N. J. A. Sloane, May 16 2008 at the suggestion of R. J. Mathar
STATUS
approved