This site is supported by donations to The OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A131133 a(0)=1; for n > 0, a(n) = (1/n + 1/a(n-1))*lcm(n, a(n-1)). 2
 1, 2, 2, 5, 9, 14, 10, 17, 25, 34, 22, 3, 5, 18, 16, 31, 47, 64, 41, 60, 4, 25, 47, 70, 47, 72, 49, 76, 26, 55, 17, 48, 5, 38, 36, 71, 107, 144, 91, 10, 5, 46, 44, 87, 131, 176, 111, 158, 103, 152, 101, 152, 51, 104, 79, 134, 95, 8, 33, 92, 38, 99, 161, 32, 3, 68, 67, 2, 35, 104 (list; graph; refs; listen; history; text; internal format)
 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 Cf. A131134. Sequence in context: A278386 A262614 A184296 * A096569 A076607 A054229 Adjacent sequences:  A131130 A131131 A131132 * A131134 A131135 A131136 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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified August 14 06:01 EDT 2018. Contains 313748 sequences. (Running on oeis4.)