cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 152 results. Next

A240926 a(n) = 2 + L(2*n) = 2 + A005248(n), n >= 0, with the Lucas numbers (A000032).

Original entry on oeis.org

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
Offset: 0

Views

Author

Kival Ngaokrajang, Aug 03 2014

Keywords

Comments

This sequence also gives the curvature of touching circles inscribed in a special way in the smaller segment of a circle of radius 5/4 cut 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 so that the larger sagitta also has 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 conditions that C_n touches (i) the circle C, (ii) the chord and (iii) the circle C_(n-1). The curvature of the n-th circle, C_n = 1/R_n, n >= 0, is 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.
Essentially a duplicate of A092387. - R. J. Mathar, Jul 07 2023

Crossrefs

Programs

  • Magma
    [2+Lucas(2*n): n in [0..40]]; // Vincenzo Librandi, Oct 08 2015
    
  • Mathematica
    Table[2 + LucasL[2 n], {n, 0, 50}] (* Vincenzo Librandi, Oct 08 2015 *)
  • PARI
    vector(100, n, n--; 2 + fibonacci(2*n-1) + fibonacci(2*n+1)) \\ Altug Alkan, Oct 08 2015

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
a(n) = (phi^n + phi^(-n))^2, where phi = A001622 = (1 + sqrt(5))/2. - Diego Rattaggi, Jun 10 2020
Sum_{k>=0} 1/a(k) = A338303. - Amiram Eldar, Oct 22 2020

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

A339129 Odd composite integers m such that A005248(m-J(m,5)) == 2 (mod m), where J(m,5) is the Jacobi symbol.

Original entry on oeis.org

9, 49, 63, 121, 169, 289, 323, 361, 377, 441, 529, 841, 961, 1127, 1369, 1681, 1849, 1891, 2209, 2303, 2809, 2961, 3481, 3721, 3751, 3827, 4181, 4489, 4901, 4961, 5041, 5329, 5491, 5777, 6137, 6241, 6601, 6721, 6889, 7381, 7921, 8149, 9409, 10201, 10609, 10877, 10933, 11449, 11663
Offset: 1

Views

Author

Ovidiu Bagdasar, Nov 24 2020

Keywords

Comments

The generalized Pell-Lucas sequences of integer parameters (a,b) defined by V(m+2)=a*V(m+1)-b*V(m) and V(0)=2, V(1)=a, satisfy the identity
V(p-J(p,D)) == 2 (mod p) when p is prime, b=1 and D=a^2-4.
This sequence contains the odd composite integers with V(m-J(m,D)) == 2 (mod m).
For a=3 and b=1, we have D=5 and V(m) recovers A005248(m).

References

  • D. Andrica, O. Bagdasar, Recurrent Sequences: Key Results, Applications and Problems. Springer, 2020.
  • D. Andrica, O. Bagdasar, On some new arithmetic properties of the generalized Lucas sequences, Mediterr. J. Math. (to appear, 2021)
  • D. Andrica, O. Bagdasar, On generalized pseudoprimality of level k (submitted)

Crossrefs

Cf. A005248.
Cf. A339125 (a=1, b=-1), A339126 (a=3, b=-1), A339127 (a=5, b=-1), A339128 (a=7, b=-1), A339130 (a=5, b=1), A339131 (a=7, b=1).

Programs

  • Mathematica
    Select[Range[3, 12000, 2], CompositeQ[#] && Divisible[LucasL[2*(# - JacobiSymbol[#, 5])] - 2, #] &] (* Amiram Eldar, Nov 26 2020 *)

A188378 Partial sums of A005248.

Original entry on oeis.org

2, 5, 12, 30, 77, 200, 522, 1365, 3572, 9350, 24477, 64080, 167762, 439205, 1149852, 3010350, 7881197, 20633240, 54018522, 141422325, 370248452, 969323030, 2537720637, 6643838880, 17393796002, 45537549125, 119218851372, 312119004990, 817138163597, 2139295485800
Offset: 0

Views

Author

Gabriele Fici, Mar 29 2011

Keywords

Comments

Different from A024851.
Luo proves that these integers cannot be uniquely decomposed as the sum of distinct and nonconsecutive terms of the Lucas number sequence. - Michel Marcus, Apr 20 2020

Crossrefs

Programs

  • Magma
    [5*Fibonacci(n)*Fibonacci(n+1)+1+(-1)^n: n in [0..40]]; // Vincenzo Librandi, Jan 24 2016
  • Maple
    f:= gfun:-rectoproc({a(n+3)-4*a(n+2)+4*a(n+1)-a(n), a(0) = 2, a(1) = 5, a(2) = 12}, a(n), remember):
    map(f, [$0..60]); # Robert Israel, Feb 02 2016
  • Mathematica
    LinearRecurrence[{4,-4,1},{2,5,12},30] (* Harvey P. Dale, Oct 05 2015 *)
    Accumulate@ LucasL@ Range[0, 58, 2] (* Michael De Vlieger, Jan 24 2016 *)
  • PARI
    a(n) = 5*fibonacci(n)*fibonacci(n+1) + 1 + (-1)^n; \\ Michel Marcus, Aug 26 2013
    
  • PARI
    Vec((-2+3*x)/((x-1)*(x^2-3*x+1)) + O(x^100)) \\ Altug Alkan, Jan 24 2016
    

Formula

a(n) = A000032(2n+1)+1 = A002878(n)+1 = 2*A027941(n+1)-3*A027941(n).
G.f.: ( -2+3*x ) / ( (x-1)*(x^2-3*x+1) ). - R. J. Mathar, Mar 30 2011
a(n) = 5*A001654(n) + 1 + (-1)^n, n>=0. [Wolfdieter Lang, Jul 23 2012]
(a(n)^3 + (a(n)-2)^3) / 2 = A000032(A016945(n)) = Lucas(6*n+3) = A267797(n), for n>0. - Altug Alkan, Jan 31 2016
a(n) = 2^(-1-n)*(2^(1+n)-(3-sqrt(5))^n*(-1+sqrt(5))+(1+sqrt(5))*(3+sqrt(5))^n). - Colin Barker, Nov 02 2016

A335672 Odd composite integers m such that A005248(m) == 3 (mod m).

Original entry on oeis.org

15, 105, 195, 231, 323, 377, 435, 665, 705, 1443, 1551, 1891, 2465, 2737, 2849, 3289, 3689, 3745, 3827, 4181, 4465, 4879, 5655, 5777, 6479, 6601, 6721, 7055, 7743, 8149, 9879, 10815, 10877, 11305, 11395, 11663, 12935, 13201, 13981, 15251, 15301, 17119, 17261, 17711, 18407, 18915, 19043, 20999
Offset: 1

Views

Author

Ovidiu Bagdasar, Jun 17 2020

Keywords

Comments

If p is a prime, then A005248(p)==3 (mod p).
This sequence contains the odd composite integers for which the congruence holds.
The generalized Pell-Lucas sequences of integer parameters (a,b) defined by V(n+2)=a*V(n+1)-b*V(n) and V(0)=2, V(1)=a, satisfy the identity V(p)==a (mod p) whenever p is prime and b=-1,1.
For a=3, b=1, V(n) recovers A005248(n) (bisection of Fibonacci numbers).

Examples

			15 is the first odd composite integer for which A005248(15)=18604984==3 (mod 15).
		

References

  • D. Andrica, O. Bagdasar, Recurrent Sequences: Key Results, Applications and Problems. Springer (to appear, 2020).

Crossrefs

Cf. A005248, A335669 (a=3,b=-1), A335673 (a=4,b=1), A335674 (a=5,b=1).

Programs

  • Mathematica
    Select[Range[3, 10000, 2], CompositeQ[#] && Divisible[LucasL[2#] - 3, #] &] (* Amiram Eldar, Jun 18 2020 *)

A337777 Even composite integers m such that U(m)^2 == 1 (mod m) and V(m) == 3 (mod m), where U(m)=A001906(m) and V(m)=A005248(m) are the m-th generalized Lucas and Pell-Lucas numbers of parameters a=3 and b=1, respectively.

Original entry on oeis.org

4, 44, 836, 1364, 2204, 7676, 7964, 9164, 11476, 12524, 23804, 31124, 32642, 39556, 73124, 80476, 99644, 110564, 128876, 156484, 192676, 199924, 287804, 295196, 315524, 398924, 542242, 715604, 780044, 934876, 987524, 1050524, 1339516, 1390724, 1891124, 1996796
Offset: 1

Views

Author

Ovidiu Bagdasar, Sep 20 2020

Keywords

Comments

For a, b integers, the following sequences are defined:
generalized Lucas sequences by U(n+2)=a*U(n+1)-b*U(n) and U(0)=0, U(1)=1;
generalized Pell-Lucas sequences by V(n+2)=a*V(n+1)-b*V(n) and V(0)=2, V(1)=a.
These satisfy the identities U(p)^2 == 1 and V(p)==a (mod p) for p prime and b=1,-1.
These numbers may be called weak generalized Lucas-Bruckner pseudoprimes of parameters a and b. The current sequence is defined for a=3 and b=1.

Crossrefs

Cf. A337626.

Programs

  • Mathematica
    Select[Range[2, 20000, 2], CompositeQ[#] && Divisible[LucasL[2#] - 3, #] && Divisible[ChebyshevU[#-1, 3/2]*ChebyshevU[#-1, 3/2] - 1, #] &]

Extensions

More terms from Amiram Eldar, Sep 21 2020

A339521 Odd composite integers m such that A005248(2*m-J(m,5)) == 3 (mod m), where J(m,5) is the Jacobi symbol.

Original entry on oeis.org

21, 203, 323, 329, 377, 451, 609, 861, 1001, 1081, 1183, 1547, 1729, 1819, 1891, 2033, 2211, 2821, 3081, 3549, 3653, 3827, 4089, 4181, 4669, 5671, 5777, 5887, 6601, 6721, 8149, 8557, 8841, 10877, 11309, 11663, 13201, 13601, 13861, 13981, 14701, 15051, 15251
Offset: 1

Views

Author

Ovidiu Bagdasar, Dec 07 2020

Keywords

Comments

The generalized Pell-Lucas sequences of integer parameters (a,b) defined by V(m+2)=a*V(m+1)-b*V(m) and V(0)=2, V(1)=a, satisfy V(k*p-J(p,D)) == V(k-1) (mod p) whenever p is prime, k is a positive integer, b=1 and D=a^2-4.
The composite integers m with the property V(k*m-J(m,D)) == V(k-1) (mod m) are called generalized Pell-Lucas pseudoprimes of level k+ and parameter a.
Here b=1, a=3, D=5 and k=2, while V(m) recovers A005248(m).

References

  • D. Andrica, O. Bagdasar, Recurrent Sequences: Key Results, Applications and Problems. Springer, 2020.
  • D. Andrica, O. Bagdasar, On some new arithmetic properties of the generalized Lucas sequences, Mediterr. J. Math. (to appear, 2021).
  • D. Andrica, O. Bagdasar, On generalized pseudoprimality of level k (submitted).

Crossrefs

Cf. A005248, A071904, A339129 (a=3, b=1, k=1).
Cf. A339522 (a=5, b=1), A339523 (a=7, b=1).

Programs

  • Mathematica
    Select[Range[3, 25000, 2], CoprimeQ[#, 5] && CompositeQ[#] && Divisible[LucasL[2*(2*# - JacobiSymbol[#, 5])] - 3, #] &]

A339728 Odd composite integers m such that A005248(3*m-J(m,5)) == 7 (mod m), where J(m,5) is the Jacobi symbol.

Original entry on oeis.org

9, 21, 27, 63, 161, 189, 207, 261, 287, 323, 341, 377, 671, 783, 861, 901, 987, 1007, 1107, 1269, 1281, 1287, 1449, 1853, 1891, 2071, 2241, 2407, 2431, 2501, 2529, 2567, 2743, 2961, 3201, 3827, 4181, 4623, 5029, 5473, 5611, 5777, 5781, 6119, 6601, 6721, 7161
Offset: 1

Views

Author

Ovidiu Bagdasar, Dec 14 2020

Keywords

Comments

The generalized Pell-Lucas sequences of integer parameters (a,b) defined by V(m+2)=a*V(m+1)-b*V(m) and V(0)=2, V(1)=a, satisfy V(k*p-J(p,D)) == V(k-1) (mod p) whenever p is prime, k is a positive integer, b=1 and D=a^2-4.
The composite integers m with the property V(k*m-J(m,D)) == V(k-1) (mod m) are called generalized Pell-Lucas pseudoprimes of level k+ and parameter a.
Here b=1, a=3, D=5 and k=3, while V(m) recovers A005248(m), with V(2)=7.

References

  • D. Andrica, O. Bagdasar, Recurrent Sequences: Key Results, Applications and Problems. Springer, 2020.
  • D. Andrica, O. Bagdasar, On some new arithmetic properties of the generalized Lucas sequences, Mediterr. J. Math. (to appear, 2021).
  • D. Andrica, O. Bagdasar, On generalized pseudoprimality of level k (submitted).

Crossrefs

Cf. A005248, A071904, A339129 (a=3, b=1, k=1), A339521 (a=3, b=1, k=2).
Cf. A339729 (a=5, b=1), A339730 (a=7, b=1).

Programs

  • Mathematica
    Select[Range[3, 7500, 2], CoprimeQ[#, 5] && CompositeQ[#] && Divisible[LucasL[2*(3*# - JacobiSymbol[#, 5])] - 7, #] &]

A345379 Number of terms m <= n, where m is a term in the bisection of Lucas numbers (A005248).

Original entry on oeis.org

0, 0, 1, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5
Offset: 0

Views

Author

Ovidiu Bagdasar, Jun 16 2021

Keywords

Examples

			a(0)=a(1)=0, since the least term in A005248 is 2.
a(2)=1 since A005248(0) = 2 is followed in that sequence by 3.
a(k)=3 for 3 <= k <= 6 since the first terms of A005248 are {0, 2, 3, 7}.
		

Crossrefs

Cf. A005248, A108852 (Fibonacci), A130245 (Lucas), A130260.

Programs

  • Mathematica
    Block[{a = 3, b = 1, nn = 105, u, v = {}}, u = {2, a}; Do[AppendTo[u, Total[{-b, a} u[[-2 ;; -1]]]]; AppendTo[v, Count[u, _?(# <= i &)]], {i, nn}]; {Boole[First[u] <= 0]}~Join~v] (* or *)
    {0}~Join~Accumulate@ ReplacePart[ConstantArray[0, Last[#]], Map[# -> 1 &, #]] &@ LucasL@ Range[0, 10, 2] (* Michael De Vlieger, Jun 16 2021 *)

A033550 a(n) = A005248(n) - n.

Original entry on oeis.org

2, 2, 5, 15, 43, 118, 316, 836, 2199, 5769, 15117, 39592, 103670, 271430, 710633, 1860483, 4870831, 12752026, 33385264, 87403784, 228826107, 599074557, 1568397585, 4106118220, 10749957098, 28143753098, 73681302221, 192900153591
Offset: 0

Views

Author

Keywords

Comments

Also distinct compositions of the wheel graph W_n. - Ralf Stephan, Jan 02 2003

Crossrefs

Programs

  • GAP
    List([0..50], n-> Lucas(1,-1, 2*n)[2] - n ); # G. C. Greubel, Oct 12 2019
  • Magma
    [Lucas(2*n) - n: n in [0..50]]; // G. C. Greubel, Oct 12 2019
    
  • Maple
    with(combinat); seq(fibonacci(2*n+1)+fibonacci(2*n-1)-n, n=0..50); # G. C. Greubel, Oct 12 2019
  • Mathematica
    Table[LucasL[2*n]-n, {n,0,50}] (* G. C. Greubel, Oct 12 2019 *)
  • PARI
    a(n)=fibonacci(2*n+1)+fibonacci(2*n-1)-n
    
  • Sage
    [lucas_number2(2*n, 1,-1) - n for n in range(50)] # G. C. Greubel, Oct 12 2019
    

Formula

a(n) = 3*a(n-1) - a(n-2) + n - 1.
G.f.: (2 - 8*x + 11*x^2 - 4*x^3)/((1-3*x+x^2)*(1-x)^2).
a(n) = Lucas(2*n) - n. - G. C. Greubel, Oct 12 2019
E.g.f.: 2*exp(3*x/2)*cosh(sqrt(5)*x/2) - exp(x)*x. - Stefano Spezia, Oct 14 2024

A000032 Lucas numbers beginning at 2: L(n) = L(n-1) + L(n-2), L(0) = 2, L(1) = 1.

Original entry on oeis.org

2, 1, 3, 4, 7, 11, 18, 29, 47, 76, 123, 199, 322, 521, 843, 1364, 2207, 3571, 5778, 9349, 15127, 24476, 39603, 64079, 103682, 167761, 271443, 439204, 710647, 1149851, 1860498, 3010349, 4870847, 7881196, 12752043, 20633239, 33385282, 54018521, 87403803
Offset: 0

Views

Author

N. J. A. Sloane, May 24 1994

Keywords

Comments

Cf. A000204 for Lucas numbers beginning with 1.
Also the number of independent vertex sets and vertex covers for the cycle graph C_n for n >= 2. - Eric W. Weisstein, Jan 04 2014
Also the number of matchings in the n-cycle graph C_n for n >= 3. - Eric W. Weisstein, Oct 01 2017
Also the number of maximal independent vertex sets (and maximal vertex covers) for the n-helm graph for n >= 3. - Eric W. Weisstein, May 27 2017
Also the number of maximal independent vertex sets (and maximal vertex covers) for the n-sunlet graph for n >= 3. - Eric W. Weisstein, Aug 07 2017
This is also the Horadam sequence (2, 1, 1, 1). - Ross La Haye, Aug 18 2003
For distinct primes p, q, L(p) is congruent to 1 mod p, L(2p) is congruent to 3 mod p and L(pq) is congruent 1 + q(L(q) - 1) mod p. Also, L(m) divides F(2km) and L((2k + 1)m), k, m >= 0.
a(n) = Sum_{k=0..ceiling((n - 1)/2)} P(3; n - 1 - k, k), n >= 1, with a(0) = 2. These are the sums over the SW-NE diagonals in P(3; n, k), the (3, 1) Pascal triangle A093560. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs. Also SW-NE diagonal sums of the (1, 2) Pascal triangle A029635 (with T(0, 0) replaced by 2).
Suppose psi = log(phi) = A002390. We get the representation L(n) = 2*cosh(n*psi) if n is even; L(n) = 2*sinh(n*psi) if n is odd. There is a similar representation for Fibonacci numbers (A000045). Many Lucas formulas now easily follow from appropriate sinh- and cosh-formulas. For example: the identity cosh^2(x) - sinh^2(x) = 1 implies L(n)^2 - 5*F(n)^2 = 4*(-1)^n (setting x = n*psi). - Hieronymus Fischer, Apr 18 2007
From John Blythe Dobson, Oct 02 2007, Oct 11 2007: (Start)
The parity of L(n) follows easily from its definition, which shows that L(n) is even when n is a multiple of 3 and odd otherwise.
The first six multiplication formulas are:
L(2n) = L(n)^2 - 2*(-1)^n;
L(3n) = L(n)^3 - 3*(-1)^n*L(n);
L(4n) = L(n)^4 - 4*(-1)^n*L(n)^2 + 2;
L(5n) = L(n)^5 - 5*(-1)^n*L(n)^3 + 5*L(n);
L(6n) = L(n)^6 - 6*(-1)^n*L(n)^4 + 9*L(n)^2 - 2*(-1)^n.
Generally, L(n) | L(mn) if and only if m is odd.
In the expansion of L(mn), where m represents the multiplier and n the index of a known value of L(n), the absolute values of the coefficients are the terms in the m-th row of the triangle A034807. When m = 1 and n = 1, L(n) = 1 and all the terms are positive and so the row sums of A034807 are simply the Lucas numbers. (End)
From John Blythe Dobson, Nov 15 2007: (Start)
The comments submitted by Miklos Kristof on Mar 19 2007 for the Fibonacci numbers (A000045) contain four important identities that have close analogs in the Lucas numbers:
For a >= b and odd b, L(a + b) + L(a - b) = 5*F(a)*F(b).
For a >= b and even b, L(a + b) + L(a - b) = L(a)*L(b).
For a >= b and odd b, L(a + b) - L(a - b) = L(a)*L(b).
For a >= b and even b, L(a + b) - L(a - b) = 5*F(a)*F(b).
A particularly interesting instance of the difference identity for even b is L(a + 30) - L(a - 30) = 5*F(a)*832040, since 5*832040 is divisible by 100, proving that the last two digits of Lucas numbers repeat in a cycle of length 60 (see A106291(100)). (End)
From John Blythe Dobson, Nov 15 2007: (Start)
The Lucas numbers satisfy remarkable difference equations, in some cases best expressed using Fibonacci numbers, of which representative examples are the following:
L(n) - L(n - 3) = 2*L(n - 2);
L(n) - L(n - 4) = 5*F(n - 2);
L(n) - L(n - 6) = 4*L(n - 3);
L(n) - L(n - 12) = 40*F(n - 6);
L(n) - L(n - 60) = 4160200*F(n - 30).
These formulas establish, respectively, that the Lucas numbers form a cyclic residue system of length 3 (mod 2), of length 4 (mod 5), of length 6 (mod 4), of length 12 (mod 40) and of length 60 (mod 4160200). The divisibility of the last modulus by 100 accounts for the fact that the last two digits of the Lucas numbers begin to repeat at L(60).
The divisibility properties of the Lucas numbers are very complex and still not fully understood, but several important criteria are established in Zhi-Hong Sun's 2003 survey of congruences for Fibonacci numbers. (End)
Sum_{n>0} a(n)/(n*2^n) = 2*log(2). - Jaume Oliver Lafont, Oct 11 2009
A010888(a(n)) = A030133(n). - Reinhard Zumkeller, Aug 20 2011
The powers of phi, the golden ratio, approach the values of the Lucas numbers, the odd powers from above and the even powers from below. - Geoffrey Caveney, Apr 18 2014
Inverse binomial transform is (-1)^n * a(n). - Michael Somos, Jun 03 2014
Lucas numbers are invariant to the following transformation for all values of the integers j and n, including negative values, thus: L(n) = (L(j+n) + (-1)^n * L(j-n))/L(j). The same transformation applied to all sequences of the form G(n+1) = m * G(n) + G(n-1) yields Lucas numbers for m = 1, except where G(j) = 0, regardless of initial values which may be nonintegers. The corresponding sequences for other values of m are: for m = 2, 2*A001333; for m = 3, A006497; for m = 4, 2*A001077; for m = 5, A087130; for m = 6, 2*A005667; for m = 7, A086902. The invariant ones all have G(0) = 2, G(1) = m. A related family of sequences is discussed at A059100. - Richard R. Forberg, Nov 23 2014
If x=a(n), y=a(n+1), z=a(n+2), then -x^2 - z*x - 3*y*x - y^2 + y*z + z^2 = 5*(-1)^(n+1). - Alexander Samokrutov, Jul 04 2015
A conjecture on the divisibility of infinite subsequences of Lucas numbers by prime(n)^m, m >= 1, is given in A266587, together with the prime "entry points". - Richard R. Forberg, Dec 31 2015
A trapezoid has three lengths of sides in order L(n-1), L(n+1), L(n-1). For increasing n a very close approximation to the maximum area will have the fourth side equal to 2*L(n). For a trapezoid with sides L(n-1), L(n-3), L(n-1), the fourth side will be L(n). - J. M. Bergot, Mar 17 2016
Satisfies Benford's law [Brown-Duncan, 1970; Berger-Hill, 2017]. - N. J. A. Sloane, Feb 08 2017
Lucas numbers L(n) and Fibonacci numbers F(n), being related by the formulas F(n) = (F(n-1) + L(n-1))/2 and L(n) = 2 F(n+1) - F(n), are a typical pair of "autosequences" (see the link to OEIS Wiki). - Jean-François Alcover, Jun 09 2017
For n >= 3, the Lucas number L(n) is the dimension of a commutative Hecke algebra of affine type A_n with independent parameters. See Theorem 1.4, Corollary 1.5, and the table on page 524 in the link "Hecke algebras with independent parameters". - Jia Huang, Jan 20 2019
From Klaus Purath, Apr 19 2019: (Start)
While all prime numbers appear as factors in the Fibonacci numbers, this is not the case with the Lucas numbers. For example, L(n) is never divisible by the following prime numbers < 150: 5, 13, 17, 37, 53, 61, 73, 89, 97, 109, 113, 137, 149 ... See A053028. Conjecture: Three properties can be determined for these prime numbers:
First observation: The prime factors > 3 occur in the Fibonacci numbers with an odd index.
Second observation: These are the prime numbers p congruent to 2, 3 (modulo 5), which occur both in Fibonacci(p+1) and in Fibonacci((p+1)/2) as prime factors, or the prime numbers p congruent to 1, 4 (modulo 5), which occur both in Fibonacci((p-1)/2) and in Fibonacci((p-1)/(2^k)) with k >= 2.
Third observation: The Pisano period lengths of these prime numbers, given in A001175, are always divisible by 4, but not by 8. In contrast, those of the prime factors of Lucas numbers are divisible either by 2, but not by 4, or by 8. (See also comment in A053028 by N. J. A. Sloane, Feb 21 2004). (End)
L(n) is the sum of 4*k consecutive terms of the Fibonacci sequence (A000045) divided by Fibonacci(2*k): (Sum_{i=0..4*k-1, k>=1} F(n+i))/F(2*k) = L(n+2*k+1). Sequences extended to negative indices, following the rule a(n-1) = a(n+1) - a(n). - Klaus Purath, Sep 15 2019
If one forms a sequence (A) of the Fibonacci type with the initial values A(0) = A022095(n) and A(1) = A000285(n), then A(n+1) = L(n+1)^2 always applies. - Klaus Purath, Sep 29 2019
From Kai Wang, Dec 18 2019: (Start)
L((2*m+1)k)/L(k) = Sum_{i=0..m-1} (-1)^(i*(k+1))*L((2*m-2*i)*k) + (-1)^(m*k).
Example: k=5, m=2, L(5)=11, L(10)=123, L(20)=15127, L(25)=167761. L(25)/L(5) = 15251, L(20) + L(10) + 1 = 15127 + 123 + 1 = 15251. (End)
From Peter Bala, Dec 23 2021: (Start)
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) ( mod p^k ) hold for all prime p and positive integers n and k.
For a positive integer k, the sequence (a(n))n>=1 taken modulo k becomes a purely periodic sequence. For example, taken modulo 11, the sequence becomes [1, 3, 4, 7, 0, 7, 7, 3, 10, 2, 1, 3, 4, 7, 0, 7, 7, 3, 10, 2, ...], a periodic sequence with period 10. (End)
For any sequence with recurrence relation b(n) = b(n-1) + b(n-2), it can be shown that the recurrence relation for every k-th term is given by: b(n) = A000032(k) * b(n-k) + (-1)^(k+1) * b(n-2k), extending to negative indices as necessary. - Nick Hobson, Jan 19 2024
For n >= 3, L(n) is the number of (n-1)-digit numbers where all consecutive pairs of digits have a difference of at least 8. - Edwin Hermann, Apr 19 2025

Examples

			G.f. = 2 + x + 3*x^2 + 4*x^3 + 7*x^4 + 11*x^5 + 18*x^6 + 29*x^7 + ...
		

References

  • P. Bachmann, Niedere Zahlentheorie (1902, 1910), reprinted Chelsea, NY, 1968, vol. 2, p. 69.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 32,50.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 499.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 46.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 112, 202-203.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.5 The Fibonacci and Related Sequences, pp. 287-288.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 148.
  • Silvia Heubach and Toufik Mansour, Combinatorics of Compositions and Words, CRC Press, 2010.
  • V. E. Hoggatt, Jr., Fibonacci and Lucas Numbers. Houghton, Boston, MA, 1969.
  • Thomas Koshy, Fibonacci and Lucas Numbers with Applications, John Wiley and Sons, 2001.
  • C. N. Menhinick, The Fibonacci Resonance and other new Golden Ratio discoveries, Onperson, (2015), pages 200-206.
  • Paulo Ribenboim, My Numbers, My Friends: Popular Lectures on Number Theory, Springer-Verlag, NY, 2000, p. 3.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 45-46, 59.
  • Michel Rigo, Formal Languages, Automata and Numeration Systems, 2 vols., Wiley, 2014. Mentions this sequence - see "List of Sequences" in Vol. 2.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • S. Vajda, Fibonacci and Lucas numbers and the Golden Section, Ellis Horwood Ltd., Chichester, 1989.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See pp. 83-84.

Crossrefs

Cf. A000204. A000045(n) = (2*L(n + 1) - L(n))/5.
First row of array A103324.
a(n) = A101220(2, 0, n), for n > 0.
a(k) = A090888(1, k) = A109754(2, k) = A118654(2, k - 1), for k > 0.
Cf. A131774, A001622, A002878 (L(2n+1)), A005248 (L(2n)), A006497, A080039, A049684 (summation of Fibonacci(4n+2)), A106291 (Pisano periods), A057854 (complement), A354265 (generalized Lucas numbers).
Cf. sequences with formula Fibonacci(n+k)+Fibonacci(n-k) listed in A280154.
Subsequence of A047201.

Programs

  • Haskell
    a000032 n = a000032_list !! n
    a000032_list = 2 : 1 : zipWith (+) a000032_list (tail a000032_list)
    -- Reinhard Zumkeller, Aug 20 2011
    
  • Magma
    [Lucas(n): n in [0..120]];
    
  • Maple
    with(combinat): A000032 := n->fibonacci(n+1)+fibonacci(n-1);
    seq(simplify(2^n*(cos(Pi/5)^n+cos(3*Pi/5)^n)), n=0..36)
  • Mathematica
    a[0] := 2; a[n] := Nest[{Last[#], First[#] + Last[#]} &, {2, 1}, n] // Last
    Array[2 Fibonacci[# + 1] - Fibonacci[#] &, 50, 0] (* Joseph Biberstine (jrbibers(AT)indiana.edu), Dec 26 2006 *)
    Table[LucasL[n], {n, 0, 36}] (* Zerinvary Lajos, Jul 09 2009 *)
    LinearRecurrence[{1, 1}, {2, 1}, 40] (* Harvey P. Dale, Sep 07 2013 *)
    LucasL[Range[0, 20]] (* Eric W. Weisstein, Aug 07 2017 *)
    CoefficientList[Series[(-2 + x)/(-1 + x + x^2), {x, 0, 20}], x] (* Eric W. Weisstein, Sep 21 2017 *)
  • PARI
    {a(n) = if(n<0, (-1)^n * a(-n), if( n<2, 2-n, a(n-1) + a(n-2)))};
    
  • PARI
    {a(n) = if(n<0, (-1)^n * a(-n), polsym(x^2 - x - 1, n)[n+1])};
    
  • PARI
    {a(n) = real((2 + quadgen(5)) * quadgen(5)^n)};
    
  • PARI
    a(n)=fibonacci(n+1)+fibonacci(n-1) \\ Charles R Greathouse IV, Jun 11 2011
    
  • PARI
    polsym(1+x-x^2, 50) \\ Charles R Greathouse IV, Jun 11 2011
    
  • Python
    def A000032_gen(): # generator of terms
        a, b = 2, 1
        while True:
            yield a
            a, b = b, a+b
    it = A000032_gen()
    A000032_list = [next(it) for  in range(50)] # _Cole Dykstra, Aug 02 2022
    
  • Python
    from sympy import lucas
    def A000032(n): return lucas(n) # Chai Wah Wu, Sep 23 2023
    
  • Python
    [(i:=3)+(j:=-1)] + [(j:=i+j)+(i:=j-i) for  in range(100)] # _Jwalin Bhatt, Apr 02 2025
  • Sage
    [lucas_number2(n,1,-1) for n in range(37)] # Zerinvary Lajos, Jun 25 2008
    

Formula

G.f.: (2 - x)/(1 - x - x^2).
L(n) = ((1 + sqrt(5))/2)^n + ((1 - sqrt(5))/2)^n = phi^n + (1-phi)^n.
L(n) = L(n - 1) + L(n - 2) = (-1)^n * L( - n).
L(n) = Fibonacci(2*n)/Fibonacci(n) for n > 0. - Jeff Burch, Dec 11 1999
E.g.f.: 2*exp(x/2)*cosh(sqrt(5)*x/2). - Len Smiley, Nov 30 2001
L(n) = F(n) + 2*F(n - 1) = F(n + 1) + F(n - 1). - Henry Bottomley, Apr 12 2000
a(n) = sqrt(F(n)^2 + 4*F(n + 1)*F(n - 1)). - Benoit Cloitre, Jan 06 2003 [Corrected by Gary Detlefs, Jan 21 2011]
a(n) = 2^(1 - n)*Sum_{k=0..floor(n/2)} C(n, 2k)*5^k. a(n) = 2T(n, i/2)( - i)^n with T(n, x) Chebyshev's polynomials of the first kind (see A053120) and i^2 = - 1. - Paul Barry, Nov 15 2003
L(n) = 2*F(n + 1) - F(n). - Paul Barry, Mar 22 2004
a(n) = (phi)^n + ( - phi)^( - n). - Paul Barry, Mar 12 2005
From Miklos Kristof, Mar 19 2007: (Start)
Let F(n) = A000045 = Fibonacci numbers, L(n) = a(n) = Lucas numbers:
L(n + m) + (-1)^m*L(n - m) = L(n)*L(m).
L(n + m) - (-1)^m*L(n - m) = 8*F(n)*F(m).
L(n + m + k) + (-1)^k*L(n + m - k) + (-1)^m*(L(n - m + k) + (-1)^k*L(n - m - k)) = L(n)*L(m)*L(k).
L(n + m + k) - (-1)^k*L(n + m - k) + (-1)^m*(L(n - m + k) - (-1)^k*L(n - m - k)) = 5*F(n)*L(m)*F(k).
L(n + m + k) + (-1)^k*L(n + m - k) - (-1)^m*(L(n - m + k) + (-1)^k*L(n - m - k)) = 5*F(n)*F(m)*L(k).
L(n + m + k) - (-1)^k*L(n + m - k) - (-1)^m*(L(n - m + k) - (-1)^k*L(n - m - k)) = 5*L(n)*F(m)*F(k). (End)
Inverse: floor(log_phi(a(n)) + 1/2) = n, for n>1. Also for n >= 0, floor((1/2)*log_phi(a(n)*a(n+1))) = n. Extension valid for all integers n: floor((1/2)*sign(a(n)*a(n+1))*log_phi|a(n)*a(n+1)|) = n {where sign(x) = sign of x}. - Hieronymus Fischer, May 02 2007
Let f(n) = phi^n + phi^(-n), then L(2n) = f(2n) and L(2n + 1) = f(2n + 1) - 2*Sum_{k>=0} C(k)/f(2n + 1)^(2k + 1) where C(n) are Catalan numbers (A000108). - Gerald McGarvey, Dec 21 2007, modified by Davide Colazingari, Jul 01 2016
Starting (1, 3, 4, 7, 11, ...) = row sums of triangle A131774. - Gary W. Adamson, Jul 14 2007
a(n) = trace of the 2 X 2 matrix [0,1; 1,1]^n. - Gary W. Adamson, Mar 02 2008
From Hieronymus Fischer, Jan 02 2009: (Start)
For odd n: a(n) = floor(1/(fract(phi^n))); for even n>0: a(n) = ceiling(1/(1 - fract(phi^n))). This follows from the basic property of the golden ratio phi, which is phi - phi^(-1) = 1 (see general formula described in A001622).
a(n) = round(1/min(fract(phi^n), 1 - fract(phi^n))), for n>1, where fract(x) = x - floor(x). (End)
E.g.f.: exp(phi*x) + exp(-x/phi) with phi: = (1 + sqrt(5))/2 (golden section). 1/phi = phi - 1. See another form given in the Smiley e.g.f. comment. - Wolfdieter Lang, May 15 2010
L(n)/L(n - 1) -> A001622. - Vincenzo Librandi, Jul 17 2010
a(n) = 2*a(n-2) + a(n-3), n>2. - Gary Detlefs, Sep 09 2010
L(n) = floor(1/fract(Fibonacci(n)*phi)), for n odd. - Hieronymus Fischer, Oct 20 2010
L(n) = ceiling(1/(1 - fract(Fibonacci(n)*phi))), for n even. - Hieronymus Fischer, Oct 20 2010
L(n) = 2^n * (cos(Pi/5)^n + cos(3*Pi/5)^n). - Gary Detlefs, Nov 29 2010
L(n) = (Fibonacci(2*n - 1)*Fibonacci(2*n + 1) - 1)/(Fibonacci(n)*Fibonacci(2*n)), n != 0. - Gary Detlefs, Dec 13 2010
L(n) = sqrt(A001254(n)) = sqrt(5*Fibonacci(n)^2 - 4*(-1)^(n+1)). - Gary Detlefs, Dec 26 2010
L(n) = floor(phi^n) + ((-1)^n + 1)/2 = A014217(n) +((-1)^n+1)/2, where phi = A001622. - Gary Detlefs, Jan 20 2011
L(n) = Fibonacci(n + 6) mod Fibonacci(n + 2), n>2. - Gary Detlefs, May 19 2011
For n >= 2, a(n) = round(phi^n) where phi is the golden ratio. - Arkadiusz Wesolowski, Jul 20 2012
a(p*k) == a(k) (mod p) for primes p. a(2^s*n) == a(n)^(2^s) (mod 2) for s = 0,1,2.. a(2^k) == - 1 (mod 2^k). a(p^2*k) == a(k) (mod p) for primes p and s = 0,1,2,3.. [Hoggatt and Bicknell]. - R. J. Mathar, Jul 24 2012
From Gary Detlefs, Dec 21 2012: (Start)
L(k*n) = (F(k)*phi + F(k - 1))^n + (F(k + 1) - F(k)*phi)^n.
L(k*n) = (F(n)*phi + F(n - 1))^k + (F(n + 1) - F(n)*phi)^k.
where phi = (1 + sqrt(5))/2, F(n) = A000045(n).
(End)
L(n) = n * Sum_{k=0..floor(n/2)} binomial(n - k,k)/(n - k), n>0 [H. W. Gould]. - Gary Detlefs, Jan 20 2013
G.f.: G(0), where G(k) = 1 + 1/(1 - (x*(5*k-1))/((x*(5*k+4)) - 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 15 2013
L(n) = F(n) + F(n-1) + F(n-2) + F(n-3). - Bob Selcoe, Jun 17 2013
L(n) = round(sqrt(L(2n-1) + L(2n-2))). - Richard R. Forberg, Jun 24 2014
L(n) = (F(n+1)^2 - F(n-1)^2)/F(n) for n>0. - Richard R. Forberg, Nov 17 2014
L(n+2) = 1 + A001610(n+1) = 1 + Sum_{k=0..n} L(k). - Tom Edgar, Apr 15 2015
L(i+j+1) = L(i)*F(j) + L(i+1)*F(j+1) with F(n)=A000045(n). - J. M. Bergot, Feb 12 2016
a(n) = (L(n+1)^2 + 5*(-1)^n)/L(n+2). - J. M. Bergot, Apr 06 2016
Dirichlet g.f.: PolyLog(s,-1/phi) + PolyLog(s,phi), where phi is the golden ratio. - Ilya Gutkovskiy, Jul 01 2016
L(n) = F(n+2) - F(n-2). - Yuchun Ji, Feb 14 2016
L(n+1) = A087131(n+1)/2^(n+1) = 2^(-n)*Sum_{k=0..n} binomial(n,k)*5^floor((k+1)/2). - Tony Foster III, Oct 14 2017
L(2*n) = (F(k+2*n) + F(k-2*n))/F(k); n >= 1, k >= 2*n. - David James Sycamore, May 04 2018
From Greg Dresden and Shaoxiong Yuan, Jul 16 2019: (Start)
L(3n + 4)/L(3n + 1) has continued fraction: n 4's followed by a single 7.
L(3n + 3)/L(3n) has continued fraction: n 4's followed by a single 2.
L(3n + 2)/L(3n - 1) has continued fraction: n 4's followed by a single -3. (End)
From Klaus Purath, Sep 15 2019: (Start)
All involved sequences extended to negative indices, following the rule a(n-1) = a(n+1) - a(n).
L(n) = (2*L(n+2) - L(n-3))/5.
L(n) = (2*L(n-2) + L(n+3))/5.
L(n) = F(n-3) + 2*F(n).
L(n) = 2*F(n+2) - 3*F(n).
L(n) = (3*F(n-1) + F(n+2))/2.
L(n) = 3*F(n-3) + 4*F(n-2).
L(n) = 4*F(n+1) - F(n+3).
L(n) = (F(n-k) + F(n+k))/F(k) with odd k>0.
L(n) = (F(n+k) - F(n-k))/F(k) with even k>0.
L(n) = A001060(n-1) - F(n+1).
L(n) = (A022121(n-1) - F(n+1))/2.
L(n) = (A022131(n-1) - F(n+1))/3.
L(n) = (A022139(n-1) - F(n+1))/4.
L(n) = (A166025(n-1) - F(n+1))/5.
The following two formulas apply for all sequences of the Fibonacci type.
(a(n-2*k) + a(n+2*k))/a(n) = L(2*k).
(a(n+2*k+1) - a(n-2*k-1))/a(n) = L(2*k+1). (End)
L(n) = F(n-k)*L(k+1) + F(n-k-1)*L(k), for all k >= 0, where F(n) = A000045(n). - Michael Tulskikh, Dec 06 2019
F(n+2*m) = L(m)*F(n+m) + (-1)^(m-1)*F(n) for all n >= 0 and m >= 0. - Alexander Burstein, Mar 31 2022
a(n) = i^(n-1)*cos(n*c)/cos(c) = i^(n-1)*cos(c*n)*sec(c), where c = Pi/2 + i*arccsch(2). - Peter Luschny, May 23 2022
From Yike Li and Greg Dresden, Aug 25 2022: (Start)
L(2*n) = 5*binomial(2*n-1,n) - 2^(2*n-1) + 5*Sum_{j=1..n/5} binomial(2*n,n+5*j) for n>0.
L(2*n+1) = 2^(2n) - 5*Sum_{j=0..n/5} binomial(2*n+1,n+5*j+3). (End)
From Andrea Pinos, Jul 04 2023: (Start)
L(n) ~ Gamma(1/phi^n) + gamma.
L(n) = Re(phi^n + e^(i*Pi*n)/phi^n). (End)
L(n) = ((Sum_{i=0..n-1} L(i)^2) - 2)/L(n-1). - Jules Beauchamp, May 03 2025
From Peter Bala, Jul 09 2025: (Start)
The following series telescope:
For k >= 1, Sum_{n >= 1} (-1)^((k+1)*(n+1)) * a(2*n*k)/(a((2*n-1)*k)*a((2*n+1)*k)) = 1/a(k)^2.
For positive even k, Sum_{n >= 1} 1/(a(k*n) - (a(k) + 2)/a(k*n)) = 1/(a(k) - 2) and
Sum_{n >= 1} (-1)^(n+1)/(a(k*n) + (a(k) - 2)/a(k*n)) = 1/(a(k) + 2).
For positive odd k, Sum_{n >= 1} 1/(a(k*n) - (-1)^n*(a(2*k) + 2)/a(k*n)) = (a(k) + 2)/(2*(a(2*k) - 2)) and
Sum_{n >= 1} (-1)^(n+1)/(a(k*n) - (-1)^n*(a(2*k) + 2)/a(k*n)) = (a(k) - 2)/(2*(a(2*k) - 2)). (End)
Showing 1-10 of 152 results. Next