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 (ctibor.zizka(AT)seznam.cz), 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: A011273 A223387 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 | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .