login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A240926 a(n) = 2 + L(2*n) = 2 + A005248(n), n >= 0, with the Lucas numbers (A000032). 15
4, 5, 9, 20, 49, 125, 324, 845, 2209, 5780, 15129, 39605, 103684, 271445, 710649, 1860500, 4870849, 12752045, 33385284, 87403805, 228826129, 599074580, 1568397609, 4106118245, 10749957124, 28143753125, 73681302249, 192900153620, 505019158609 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,1

COMMENTS

This sequence gives also the curvature of touching circles inscribed in a special way in the smaller segment of circle of radius 5/4 divided by a chord of length 2.

Consider a circle C of radius 5/4 (in some length units) with a chord of length 2. This has been chosen such that the larger sagitta has also length 2. The smaller sagitta has length 1/2. The input, besides the circle C is the circle C_0 with radius R_0 = 1/4, touching the chord and circle C. The following sequence of circles C_n with radii R_n, n >= 1, is obtained from the condition that C_n touches i) the circle C, ii) the chord and iii) the circle C)(n-1). The circle curvatures C_n = 1/R_n, n >= 0, are conjectured to be a(n). See an illustration given in the link. As found by Wolfdieter Lang (see part II of the proof given by W. Lang in the link), this circle problem is related to the nonnegative solutions of the Pell equation X^2 - 5*Y^2 = 4: a(n) = 2 + X(n) = 2 + A005248(n). For the larger segment below the chord (with sagitta length 2) the sequence would be A115032, see W. Lang's proof given in part I of the link.

  If the circle radius and the sagitta length were both equal to 1, the curvature sequence would be A099938.

LINKS

G. C. Greubel, Table of n, a(n) for n = 0..2375

Kival Ngaokrajang, Illustration of initial terms

Wolfdieter Lang, Proof of the coincidence of a(n) with the touching circle problem (part II).

Wolfdieter Lang, Figures for various touching circle problems.

Index entries for sequences related to Chebyshev polynomials.

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

FORMULA

Conjectures (proved in the next entry) from Colin Barker, Aug 25 2014 (and Aug 27 2014): (Start)

a(n) = (2+(1/2*(3-sqrt(5)))^n+(1/2*(3+sqrt(5)))^n).

a(n) = 4*a(n-1) - 4*a(n-2) + a(n-3).

G.f.: -(5*x^2-11*x+4) / ((x-1)*(x^2-3*x+1)). (End)

From Wolfdieter Lang, Aug 26 2014 (Start)

a(n) = 2 + S(n, 3) - S(n-2, 3) = 2 + 2*S(n, 3)  - 3*S(n-1, 3).

a(n) = 3*a(n-1) - a(n-2) - 2, n >= 1, with a(-1)= 5 and a(0) = 4 (from the S(n, 3) recurrence or from A005248).

The first of the Colin Barker conjectures above is true because of the Binet-de Moivre formula for L(2*n) (see the Jul 24 2003 Dennis P. Walsh comment on A005248). With phi = (1+sqrt(5))/2, use 1/phi = phi-1, phi^2 = phi+1, (phi-1)^2 = 2 - phi.

His third conjecture (the g.f.) follows from the g.f. of A005248 by adding 2/(1-x).

His second conjecture (recurrence) with input a(-3) = 20, a(-2) = 9 and a(-1) = 5 (from the above given recurrence) leads to his g.f. with the expanded denominator. Thus all three conjectures are true.  (End)

a(n) = A005592(n) + 3, with n>0. - Zino Magri, Feb 16 2015

MATHEMATICA

Table[2 + LucasL[2 n], {n, 0, 50}] (* Vincenzo Librandi, Oct 08 2015 *)

PROG

(MAGMA) [2+Lucas(2*n): n in [0..40]]; // Vincenzo Librandi, Oct 8 2015

(PARI) vector(100, n, n--; 2 + fibonacci(2*n-1) + fibonacci(2*n+1)) \\ Altug Alkan, Oct 08 2015

CROSSREFS

Cf. A000032, A005248, A005592, A099938, A115032, A246638, A246639, A246640, A246641, A246642.

Sequence in context: A216225 A276645 A293282 * A041032 A041993 A153068

Adjacent sequences:  A240923 A240924 A240925 * A240927 A240928 A240929

KEYWORD

nonn,easy

AUTHOR

Kival Ngaokrajang, Aug 03 2014

EXTENSIONS

Edited: name changed (after proof has been given in part II of the W. Lang link), comments rewritten, cross refs. and link to Chebyshev index added. - Wolfdieter Lang, Aug 26 2014

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 18 18:58 EDT 2018. Contains 316323 sequences. (Running on oeis4.)