login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 59th year, we have over 358,000 sequences, and we’ve crossed 10,300 citations (which often say “discovered thanks to the OEIS”).

Other ways to Give
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A145051 Numerator of the first convergent to sqrt(n) using the recursion x = (n/x + x)/2. 9
1, 3, 2, 5, 3, 7, 4, 9, 5, 11, 6, 13, 7, 15, 8, 17, 9, 19, 10, 21, 11, 23, 12, 25, 13, 27, 14, 29, 15, 31, 16, 33, 17, 35, 18, 37, 19, 39, 20, 41, 21, 43, 22, 45, 23, 47, 24, 49, 25, 51, 26, 53, 27, 55, 28, 57, 29, 59, 30, 61, 31, 63, 32, 65, 33, 67, 34, 69, 35, 71, 36, 73, 37, 75 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

This is the same as A026741 without the first 2 terms in A026741. The link describes the experimental derivation of the generating function.

From Jaroslav Krizek, May 28 2010: (Start)

Numerators of arithmetic means of the first n positive integers for n >= 1.

See A040001 - denominators of arithmetic means of the first n positive integers.

a(n) = A026741(n+1) = A000217(n) * A040001(n) / n. (End)

Minimum number of line segments to draw into a circle to partition the circle into n+1 congruent circular sectors, i.e., minimum number of straight cuts required to cut a circular cake into n+1 equal slices. - Felix Fröhlich, Sep 01 2015

LINKS

Table of n, a(n) for n=1..74.

Cino Hilliard, Roots by Recursion [Broken link]

Index entries for linear recurrences with constant coefficients, signature (0,2,0,-1).

FORMULA

From Paul Barry, Nov 22 2009: (Start)

G.f.: x*(1 + 3*x - x^3)/(1 - x^2)^2.

a(n+1) = (n + 2)*(3 - (-1)^n)/4;

a(n+1) = Sum_{k=0..n, if(k=floor(n/2) or k=floor((n+1)/2),1,0)*(k+1)}. (End)

EXAMPLE

n=1, x=1; x = (1/1+1)/2 = 1/1;

n=2, x=1; x = (2/1+1)/2 = 3/2;

n=3, x=1; x = (3/1+1)/2 = 2/1.

G.f.: x + 3*x^2 + 2*x^3 + 5*x^4 + 3*x^5 + 7*x^6 + 4*x^7 + 9*x^8 + 5*x^9 + ...

MATHEMATICA

lst={}; Do[a=n^2+n; b=n^2-n; c=a/b; AppendTo[lst, Denominator[c]], {n, 2, 5!}]; lst (* Vladimir Joseph Stephan Orlovsky, Oct 20 2009 *)

PROG

(PARI) g(n, p) = x=1; for(j=1, p, x=(n/x+x)/2; if(j==1, print1(numerator(x), ", ")))

for(k=1, 100, g(k, 1))

(Magma) [(n+1)*(3 - (-1)^(n-1))/4: n in [1..100]]; // Vincenzo Librandi, Sep 02 2015

CROSSREFS

Cf. A026741.

Sequence in context: A323462 A030640 A176447 * A026741 A105658 A083242

Adjacent sequences: A145048 A145049 A145050 * A145052 A145053 A145054

KEYWORD

nonn,frac,easy

AUTHOR

Cino Hilliard, Sep 30 2008

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified December 5 19:04 EST 2022. Contains 358588 sequences. (Running on oeis4.)