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.

Previous Showing 21-30 of 152 results. Next

A056854 a(n) = Lucas(4*n).

Original entry on oeis.org

2, 7, 47, 322, 2207, 15127, 103682, 710647, 4870847, 33385282, 228826127, 1568397607, 10749957122, 73681302247, 505019158607, 3461452808002, 23725150497407, 162614600673847, 1114577054219522, 7639424778862807, 52361396397820127, 358890350005878082, 2459871053643326447
Offset: 0

Views

Author

Barry E. Williams, Aug 29 2000

Keywords

Comments

a(n) and b(n) := A004187(n) are the nonnegative proper and improper solutions of the Pell equation a(n)^2 - 5*(3*b(n))^2 = +4. See the cross-reference to A004187 below. - Wolfdieter Lang, Jun 26 2013
Lucas numbers of the form n^2-2. - Michel Lagneau, Aug 11 2014

Examples

			Pell equation: n = 0, 2^2 - 45*0^2 = +4 (improper);  n = 1, 7^2 - 5*(3*1)^2 = +4; n=2, 47^2 - 5*(3*7)^2 = +4. - _Wolfdieter Lang_, Jun 26 2013
		

References

  • R. P. Stanley. Enumerative combinatorics. Vol. 2, volume 62 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 1999.

Crossrefs

Cf. quadrisection of A000032: this sequence (first), A056914 (second), A246453 (third, without 11), A288913 (fourth).
Cf. Lucas(k*n): A000032 (k = 1), A005248 (k = 2), A014448 (k = 3), A001946 (k = 5), A087215 (k = 6), A087281 (k = 7), A087265 (k = 8), A087287 (k = 9), A065705 (k = 10), A089772 (k = 11), A089775 (k = 12).

Programs

  • Magma
    [Lucas(4*n): n in [0..100]]; // Vincenzo Librandi, Apr 14 2011
  • Mathematica
    a[0] = 2; a[1] = 7; a[n_] := 7a[n - 1] - a[n - 2]; Table[ a[n], {n, 0, 19}] (* Robert G. Wilson v, Jan 30 2004 *)
    LinearRecurrence[{7,-1},{2,7},25] (* or *) LucasL[4*Range[0,25]] (* Harvey P. Dale, Aug 08 2011 *)
  • PARI
    a(n)=if(n<0,0,polsym(1-7*x+x^2,n)[n+1])
    
  • PARI
    a(n)=if(n<0,0,2*subst(poltchebi(n),x,7/2))
    
  • Sage
    [lucas_number2(n,7,1) for n in range(27)] #Zerinvary Lajos, Jun 25 2008
    

Formula

a(n) = 7*a(n-1) - a(n-2) with a(0)=2, a(1)=7.
a(n) = A000032(4*n), where A000032 = Lucas numbers.
a(n) = 7*S(n-1, 7) - 2*S(n-2, 7) = S(n, 7) - S(n-2, 7) = 2*T(n, 7/2), with S(n, x) := U(n, x/2), S(-1, x) := 0, S(-2, x) := -1. U(n, x), resp. T(n, x), are Chebyshev's polynomials of the second, resp. first, kind. S(n-1, 7) = A004187(n), n>=0. See A049310 and A053120.
a(n) = ((7+sqrt(45))/2)^n + ((7-sqrt(45))/2)^n.
G.f.: (2-7x)/(1-7x+x^2).
a(n) = A005248(2*n); bisection of A005248.
a(n) = Fibonacci(8*n)/Fibonacci(4*n), n>0. - Gary Detlefs, Dec 26 2010
a(n) = 2 + 5*Fibonacci(2*n)^2 = 2 + 5*A049684(n), n >= 0. This is in Koshy's book (reference under A065563) 15. on p. 88. Compare with the above Chebyshev T formula. - Wolfdieter Lang, Aug 27 2012
From Peter Bala, Jan 06 2013: (Start)
Let F(x) = Product_{n = 0..inf} (1 + x^(4*n+1))/(1 + x^(4*n+3)). Let alpha = 1/2*(7 - 3*sqrt(5)). This sequence gives the simple continued fraction expansion of 1 + F(alpha) = 2.14242 42709 40138 85949 ... = 2 + 1/(7 + 1/(47 + 1/(322 + ...))).
Also F(-alpha) = 0.85670 72882 04563 14901 ... has the continued fraction representation 1 - 1/(7 - 1/(47 - 1/(322 - ...))) and the simple continued fraction expansion 1/(1 + 1/((7-2) + 1/(1 + 1/((47-2) + 1/(1 + 1/((322-2) + 1/(1 + ...))))))). Cf. A005248.
F(alpha)*F(-alpha) has the simple continued fraction expansion 1/(1 + 1/((7^2-4) + 1/(1 + 1/((47^2-4) + 1/(1 + 1/((322^2-4) + 1/(1 + ...))))))).
Added Oct 13 2019: 1/2 + (1/2)*F(alpha)/F(-alpha) = 1.16675297774947414828... has the simple continued fraction expansion 1 + 1/((7 - 2) + 1/(1 + 1/((322 - 2) + 1/(1 + 1/(15127 - 2) + 1/(1 + ...))))). (End)
a(n) = Fibonacci(4*n+2) - Fibonacci(4*n-2), where Fibonacci(-2) = -1. - Bruno Berselli, May 25 2015
a(n) = sqrt(45*(A004187(n))^2+4).
From Peter Bala, Oct 13 2019: (Start)
a(n) = F(4*n+4)/F(4) - F(4*n-4)/F(4) = A004187(n+1) - A004187(n-1).
a(n) = trace(M^n), where M is the 2 X 2 matrix [0, 1; 1, 1]^4 = [2, 3; 3, 5].
Consequently the Gauss congruences hold: a(n*p^k) = a(n*p^(k-1)) (mod p^k) for all prime p and positive integers n and k. See Zarelua and also Stanley (Ch. 5, Ex. 5.2(a) and its solution).
5*Sum_{n >= 1} 1/(a(n) - 9/a(n)) = 1: (9 = Lucas(4)+2 and 5 = Lucas(4)-2)
9*Sum_{n >= 1} (-1)^(n+1)/(a(n) + 5/a(n)) = 1.
Sum_{n >= 1} 1/a(n) = (1/4)*( theta_3((7-3*sqrt(5))/2)^2 - 1 ), where theta_3(q) = 1 + 2*Sum_{n >= 1} q^n^2. Cf. A153415.
Sum_{n >= 1} (-1)^(n+1)/a(n) = (1/4)*( 1 - theta_3((3*sqrt(5)-7)/2)^2 ).
x*exp(Sum_{n >= 1} a(n)*x^/n) = x + 7*x^2 + 48*x^3 + 329*x^4 + ... is the o.g.f. for A004187. (End)
E.g.f.: 2*exp(7*x/2)*cosh(3*sqrt(5)*x/2). - Stefano Spezia, Oct 18 2019
a(2k+1)/7 is the numerator of the continued fraction [3*sqrt(5), 3*sqrt(5), ..., 3*sqrt(5)] with 2k copies of 3*sqrt(5), for k>0. - Greg Dresden and Tracy Z. Wu, Sep 10 2020
a(n) = Sum_{k>=1} Lucas(2*n*k)/(Lucas(2*n)^k). - Diego Rattaggi, Jan 20 2025

Extensions

More terms from James Sellers, Aug 31 2000
Chebyshev comments from Wolfdieter Lang, Oct 31 2002

A004146 Alternate Lucas numbers - 2.

Original entry on oeis.org

0, 1, 5, 16, 45, 121, 320, 841, 2205, 5776, 15125, 39601, 103680, 271441, 710645, 1860496, 4870845, 12752041, 33385280, 87403801, 228826125, 599074576, 1568397605, 4106118241, 10749957120, 28143753121, 73681302245, 192900153616, 505019158605, 1322157322201
Offset: 0

Views

Author

Keywords

Comments

This is the r=5 member in the r-family of sequences S_r(n) defined in A092184 where more information can be found.
Number of spanning trees of the wheel W_n on n+1 vertices. - Emeric Deutsch, Mar 27 2005
Also number of spanning trees of the n-helm graph. - Eric W. Weisstein, Jul 16 2011
a(n) is the smallest number requiring n terms when expressed as a sum of Lucas numbers (A000204). - David W. Wilson, Jan 10 2006
This sequence has a primitive prime divisor for all terms beyond the twelfth. - Anthony Flatters (Anthony.Flatters(AT)uea.ac.uk), Aug 17 2007
From Giorgio Balzarotti, Mar 11 2009: (Start)
Determinant of power series of gamma matrix with determinant 1:
a(n) = Determinant(A + A^2 + A^3 + A^4 + A^5 + ... + A^n)
where A is the submatrix A(1..2,1..2) of the matrix with factorial determinant
A = [[1,1,1,1,1,1,...],[1,2,1,2,1,2,...],[1,2,3,1,2,3,...],[1,2,3,4,1,2,...],
[1,2,3,4,5,1,...],[1,2,3,4,5,6,...],...]. Note: Determinant A(1..n,1..n)= (n-1)!.
See A158039, A158040, A158041, A158042, A158043, A158044, for sequences of matrix 2!,3!,... (End)
The previous comment could be rephrased as: a(n) = -det(A^n - I) where I is the 2 X 2 identity matrix and A = [1, 1; 1, 2]. - Peter Bala, Mar 20 2015
a(n) is also the number of points of Arnold's "cat map" that are on orbits of period n-1. This is a map of the two-torus T^2 into itself. If we regard T^2 as R^2 / Z^2, the action of this map on a two vector in R^2 is multiplication by the unit-determinant matrix A = [2, 1;1, 1], with the vector components taken modulo one. As such, an explicit formula for the n-th entry of this sequence is -det(I-A^n). - Bruce Boghosian, Apr 26 2009
7*a(n) gives the total number of vertices in a heptagonal hyperbolic lattice {7,3} with n total levels, in which an open heptagon is centered at the origin. - Robert M. Ziff, Apr 10 2011
The sequence is the case P1 = 5, P2 = 6, Q = 1 of the 3 parameter family of 4th-order linear divisibility sequences found by Williams and Guy. - Peter Bala, Apr 03 2014
Determinants of the spiral knots S(3,k,(1,-1)). a(k) = det(S(3,k,(1,-1))). These knots are also the weaving knots W(k,3) and the Turk's Head Links THK(3,k). - Ryan Stees, Dec 14 2014
Even-indexed Fibonacci numbers (1, 3, 8, 21, ...) convolved with (1, 2, 2, 2, 2, ...). - Gary W. Adamson, Aug 09 2016
a(n) is the number of ways to tile a bracelet of length n with 1-color square, 2-color dominos, 3-color trominos, etc. - Yu Xiao, May 23 2020
a(n) is the number of face-labeled unfoldings of a pyramid whose base is a simple n-gon. Cf. A103536. - Rick Mabry, Apr 17 2023

Examples

			For k=3, b(3) = sqrt(5)*b(2) - b(1) = 5 - 1 = 4, so det(S(3,3,(1,-1))) = 4^2 = 16.
G.f. = x + 5*x^2 + 16*x^3 + 45*x^4 + 121*x^5 + 320*x^4 + 841*x^5 + ... - _Michael Somos_, Feb 10 2023
		

References

  • I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Wiley, N.Y., 1983, (p. 193, Problem 3.3.40 (a)).
  • N. Hartsfield and G. Ringel, Pearls in Graph Theory, p. 102. Academic Press: 1990.
  • B. Hasselblatt and A. Katok, "Introduction to the Modern Theory of Dynamical Systems," Cambridge University Press, 1997.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

This is the r=5 member of the family S_r(n) defined in A092184.
Cf. A005248. Partial sums of A002878. Pairwise sums of A027941. Bisection of A074392.
Sequence A032170, the Möbius transform of this sequence, is then the number of prime periodic orbits of Arnold's cat map. - Bruce Boghosian, Apr 26 2009
Cf. A103536 for the number of geometrically distinct edge-unfoldings of the regular pyramid. - Rick Mabry, Apr 17 2023

Programs

  • Magma
    [Lucas(n)-2: n in [0..60 by 2]]; // Vincenzo Librandi, Mar 20 2015
  • Mathematica
    Table[LucasL[2*n] - 2, {n, 0, 20}]
    (* Second program: *)
    LinearRecurrence[{4, -4, 1}, {0, 1, 5}, 30] (* Jean-François Alcover, Jan 08 2019 *)
  • PARI
    a(n) = { we = quadgen(5);((1+we)^n) + ((2-we)^n) - 2;} /* Michel Marcus, Aug 18 2012 */
    

Formula

a(n) = A005248(n) - 2 = A000032(2*n) - 2.
a(n+1) = 3*a(n) - a(n-1) + 2.
G.f.: x*(1+x)/(1-4*x+4*x^2-x^3) = x*(1+x)/((1-x)*(1-3*x+x^2)).
a(n) = 2*(T(n, 3/2)-1) with Chebyshev's polynomials T(n, x) of the first kind. See their coefficient triangle A053120.
a(n) = 4*a(n-1) - 4*a(n-2) + a(n-3), n>=3, a(0)=0, a(1)=1, a(2)=5.
a(n) = 2*T(n, 3/2) - 2, with twice the Chebyshev polynomials of the first kind, 2*T(n, x=3/2) = A005248(n).
a(n) = b(n) + b(n-1), n>=1, with b(n):=A027941(n-1), n>=1, b(-1):=0, the partial sums of S(n, 3) = U(n, 3/2) = A001906(n+1), with S(n, x) = U(n, x/2) Chebyshev's polynomials of the second kind.
a(2n) = A000204(2n)^2 - 4 = 5*A000045(2n)^2; a(2n+1) = A000204(2n+1)^2. - David W. Wilson, Jan 10 2006
a(n) = ((3+sqrt(5))/2)^n + ((3-sqrt(5))/2)^n - 2. - Felix Goldberg (felixg(AT)tx.technion.ac.il), Jun 09 2001
a(n) = b(n-1) + b(n-2), n>=1, with b(n):=A027941(n), b(-1):=0, partial sums of S(n, 3) = U(n, 3/2) = A001906(n+1), Chebyshev's polynomials of the second kind.
a(n) = n*Sum_{k=1..n} binomial(n+k-1,2*k-1)/k, n > 0. - Vladimir Kruchinin, Sep 03 2010
a(n) = floor(tau^(2*n)*(tau^(2*n) - floor(tau^(2*n)))), where tau = (1+sqrt(5))/2. - L. Edson Jeffery, Aug 26 2013
From Peter Bala, Apr 03 2014: (Start)
a(n) = U(n-1,sqrt(5)/2)^2, for n >= 1, where U(n,x) denotes the Chebyshev polynomial of the second kind.
a(n) = the bottom left entry of the 2 X 2 matrix T(n, M), where M is the 2 X 2 matrix [0, -3/2; 1, 5/2] and T(n,x) denotes the Chebyshev polynomial of the first kind.
See the remarks in A100047 for the general connection between Chebyshev polynomials of the first kind and 4th-order linear divisibility sequences. (End)
a(k) = det(S(3,k,(1,-1))) = b(k)^2, where b(1)=1, b(2)=sqrt(5), b(k)=sqrt(5)*b(k-1) - b(k-2) = b(2)*b(k-1) - b(k-2). - Ryan Stees, Dec 14 2014
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + Sum_{n >= 1} Fibonacci(2*n)*x^n. Cf. A001350. - Peter Bala, Mar 19 2015
E.g.f.: exp(phi^2*x) + exp(x/phi^2) - 2*exp(x), where phi = (1 + sqrt(5))/2. - G. C. Greubel, Aug 24 2015
a(n) = a(-n) for all n in Z. - Michael Somos, Aug 27 2015
From Peter Bala, Jun 03 2016: (Start)
a(n) = Lucas(2*n) - Lucas(0*n);
a(n)^2 = Lucas(4*n) - 3*Lucas(2*n) + 3*Lucas(0*n) - Lucas(-2*n);
a(n)^3 = Lucas(6*n) - 5*Lucas(4*n) + 10*Lucas(2*n) - 10*Lucas(0*n) + 5*Lucas(-2*n) - Lucas(-4*n) and so on (follows from Binet's formula for Lucas(2*n) and the algebraic identity (x + 1/x - 2)^m = f(x) + f(1/x) where f(x) = (x - 1)^(2*m - 1)/x^(m-1) ). (End)
Limit_{n->infinity} a(n+1)/a(n) = (3 + sqrt(5))/2 = A104457. - Ilya Gutkovskiy, Jun 03 2016
a(n) = (phi^n - phi^(-n))^2, where phi = A001622 = (1 + sqrt(5))/2. - Diego Rattaggi, Jun 10 2020
a(n) = 4*sinh(n*A002390)^2, where A002390 = arcsinh(1/2). - Gleb Koloskov, Sep 18 2021
a(n) = 5*F(n)^2 = L(n)^2 - 4 if n even and a(n) = L(n)^2 = 5*F(n)^2 - 4 if n odd. - Michael Somos, Feb 10 2023
a(n) = n^2 + Sum_{i=1..n-1} i*a(n-i). - Fern Gossow, Dec 03 2024

Extensions

Correction to formula from Nephi Noble (nephi(AT)math.byu.edu), Apr 09 2002
Chebyshev comments from Wolfdieter Lang, Sep 10 2004

A053122 Triangle of coefficients of Chebyshev's S(n,x-2) = U(n,x/2-1) polynomials (exponents of x in increasing order).

Original entry on oeis.org

1, -2, 1, 3, -4, 1, -4, 10, -6, 1, 5, -20, 21, -8, 1, -6, 35, -56, 36, -10, 1, 7, -56, 126, -120, 55, -12, 1, -8, 84, -252, 330, -220, 78, -14, 1, 9, -120, 462, -792, 715, -364, 105, -16, 1, -10, 165, -792, 1716, -2002, 1365, -560, 136, -18, 1, 11, -220, 1287, -3432, 5005, -4368, 2380, -816, 171, -20
Offset: 0

Views

Author

Keywords

Comments

Apart from signs, identical to A078812.
Another version with row-leading 0's and differing signs is given by A285072.
G.f. for row polynomials S(n,x-2) (signed triangle): 1/(1+(2-x)*z+z^2). Unsigned triangle |a(n,m)| has g.f. 1/(1-(2+x)*z+z^2) for row polynomials.
Row sums (signed triangle) A049347(n) (periodic(1,-1,0)). Row sums (unsigned triangle) A001906(n+1)=F(2*(n+1)) (even-indexed Fibonacci).
In the language of Shapiro et al. (see A053121 for the reference) such a lower triangular (ordinary) convolution array, considered as a matrix, belongs to the Bell-subgroup of the Riordan-group.
The (unsigned) column sequences are A000027, A000292, A000389, A000580, A000582, A001288 for m=0..5, resp. For m=6..23 they are A010966..(+2)..A011000 and for m=24..49 they are A017713..(+2)..A017763.
Riordan array (1/(1+x)^2,x/(1+x)^2). Inverse array is A039598. Diagonal sums have g.f. 1/(1+x^2). - Paul Barry, Mar 17 2005. Corrected by Wolfdieter Lang, Nov 13 2012.
Unsigned version is in A078812. - Philippe Deléham, Nov 05 2006
Also row n gives (except for an overall sign) coefficients of characteristic polynomial of the Cartan matrix for the root system A_n. - Roger L. Bagula, May 23 2007
From Wolfdieter Lang, Nov 13 2012: (Start)
The A-sequence for this Riordan triangle is A115141, and the Z-sequence is A115141(n+1), n>=0. For A- and Z-sequences for Riordan matrices see the W. Lang link under A006232 with details and references.
S(n,x^2-2) = sum(r(j,x^2),j=0..n) with Chebyshev's S-polynomials and r(j,x^2) := R(2*j+1,x)/x, where R(n,x) are the monic integer Chebyshv T-polynomials with coefficients given in A127672. Proof from comparing the o.g.f. of the partial sum of the r(j,x^2) polynomials (see a comment on the signed Riordan triangle A111125) with the present Riordan type o.g.f. for the row polynomials with x -> x^2. (End)
S(n,x^2-2) = S(2*n+1,x)/x, n >= 0, from the odd part of the bisection of the o.g.f. - Wolfdieter Lang, Dec 17 2012
For a relation to a generator for the Narayana numbers A001263, see A119900, whose columns are unsigned shifted rows (or antidiagonals) of this array, referring to the tables in the example sections. - Tom Copeland, Oct 29 2014
The unsigned rows of this array are alternating rows of a mirrored A011973 and alternating shifted rows of A030528 for the Fibonacci polynomials. - Tom Copeland, Nov 04 2014
Boas-Buck type recurrence for column k >= 0 (see Aug 10 2017 comment in A046521 with references): a(n, m) = (2*(m + 1)/(n - m))*Sum_{k = m..n-1} (-1)^(n-k)*a(k, m), with input a(n, n) = 1, and a(n,k) = 0 for n < k. - Wolfdieter Lang, Jun 03 2020
Row n gives the characteristic polynomial of the (n X n)-matrix M where M[i,j] = 2 if i = j, -1 if |i-j| = 1 and 0 otherwise. The matrix M is positive definite and has 2-condition number (cot(Pi/(2*n+2)))^2. - Jianing Song, Jun 21 2022
Also the convolution triangle of (-1)^(n+1)*n. - Peter Luschny, Oct 07 2022

Examples

			The triangle a(n,m) begins:
n\m   0    1    2     3     4     5     6    7    8  9
0:    1
1:   -2    1
2:    3   -4    1
3:   -4   10   -6     1
4:    5  -20   21    -8     1
5:   -6   35  -56    36   -10     1
6:    7  -56  126  -120    55   -12     1
7:   -8   84 -252   330  -220    78   -14    1
8:    9 -120  462  -792   715  -364   105  -16    1
9:  -10  165 -792  1716 -2002  1365  -560  136  -18  1
... Reformatted and extended by _Wolfdieter Lang_, Nov 13 2012
E.g., fourth row (n=3) {-4,10,-6,1} corresponds to the polynomial S(3,x-2) = -4+10*x-6*x^2+x^3.
From _Wolfdieter Lang_, Nov 13 2012: (Start)
Recurrence: a(5,1) = 35 = 1*5 + (-2)*(-20) -1*(10).
Recurrence from Z-sequence [-2,-1,-2,-5,...]: a(5,0) = -6 = (-2)*5 + (-1)*(-20) + (-2)*21 + (-5)*(-8) + (-14)*1.
Recurrence from A-sequence [1,-2,-1,-2,-5,...]: a(5,1) = 35 = 1*5  + (-2)*(-20) + (-1)*21 + (-2)*(-8) + (-5)*1.
(End)
E.g., the fourth row (n=3) {-4,10,-6,1} corresponds also to the polynomial S(7,x)/x = -4 + 10*x^2 - 6*x^4 + x^6. - _Wolfdieter Lang_, Dec 17 2012
Boas-Buck type recurrence: -56 = a(5, 2) = 2*(-1*1 + 1*(-6) - 1*21) = -2*28 = -56. - _Wolfdieter Lang_, Jun 03 2020
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2. ed., Wiley, New York, 1990.
  • R. N. Cahn, Semi-Simple Lie Algebras and Their Representations, Dover, NY, 2006, ISBN 0-486-44999-8, p. 62.
  • Sigurdur Helgasson, Differential Geometry, Lie Groups and Symmetric Spaces, Graduate Studies in Mathematics, volume 34. A. M. S.: ISBN 0-8218-2848-7, 1978, p. 463.

Crossrefs

Cf. A285072 (version with row-leading 0's and differing signs). - Eric W. Weisstein, Apr 09 2017

Programs

  • Maple
    seq(seq((-1)^(n+m)*binomial(n+m+1,2*m+1),m=0..n),n=0..10); # Robert Israel, Oct 15 2014
    # Uses function PMatrix from A357368. Adds a row above and a column to the left.
    PMatrix(10, n -> -(-1)^n*n); # Peter Luschny, Oct 07 2022
  • Mathematica
    T[n_, m_, d_] := If[ n == m, 2, If[n == m - 1 || n == m + 1, -1, 0]]; M[d_] := Table[T[n, m, d], {n, 1, d}, {m, 1, d}]; a = Join[M[1], Table[CoefficientList[Det[M[d] - x*IdentityMatrix[d]], x], {d, 1, 10}]]; Flatten[a] (* Roger L. Bagula, May 23 2007 *)
    (* Alternative code for the matrices from MathWorld: *)
    sln[n_] := 2IdentityMatrix[n] - PadLeft[PadRight[IdentityMatrix[n - 1], {n, n - 1}], {n, n}] - PadLeft[PadRight[IdentityMatrix[n - 1], {n - 1, n}], {n, n}] (* Roger L. Bagula, May 23 2007 *)
  • Sage
    @CachedFunction
    def A053122(n,k):
        if n< 0: return 0
        if n==0: return 1 if k == 0 else 0
        return A053122(n-1,k-1)-A053122(n-2,k)-2*A053122(n-1,k)
    for n in (0..9): [A053122(n,k) for k in (0..n)] # Peter Luschny, Nov 20 2012

Formula

a(n, m) := 0 if n
a(n, m) = -2*a(n-1, m) + a(n-1, m-1) - a(n-2, m), a(n, -1) := 0 =: a(-1, m), a(0, 0)=1, a(n, m) := 0 if n
O.g.f. for m-th column (signed triangle): ((x/(1+x)^2)^m)/(1+x)^2.
From Jianing Song, Jun 21 2022: (Start)
T(n,k) = [x^k]f_n(x), where f_{-1}(x) = 0, f_0(x) = 1, f_n(x) = (x-2)*f_{n-1}(x) - f_{n-2}(x) for n >= 2.
f_n(x) = (((x-2+sqrt(x^2-4*x))/2)^(n+1) - ((x-2-sqrt(x^2-4*x))/2)^(n+1))/sqrt(x^2-4x).
The roots of f_n(x) are 2 + 2*cos(k*Pi/(n+1)) = 4*(cos(k*Pi/(2*n+2)))^2 for 1 <= k <= n. (End)

A064831 Partial sums of A001654, or sum of the areas of the first n Fibonacci rectangles.

Original entry on oeis.org

0, 1, 3, 9, 24, 64, 168, 441, 1155, 3025, 7920, 20736, 54288, 142129, 372099, 974169, 2550408, 6677056, 17480760, 45765225, 119814915, 313679521, 821223648, 2149991424, 5628750624, 14736260449, 38580030723, 101003831721
Offset: 0

Author

Howard Stern (hsstern(AT)mindspring.com), Oct 23 2001

Keywords

Comments

The n-th rectangle is F(n)*F(n+1), where F(n) = n-th Fibonacci number (F(1)=1, F(2)=1, F(3)=2, etc.), A000045.
If 2*T(a_n) = the oblong number formed by substituting a(n) in the product formula x(x+1), then 2*T(a_n) = F(n-1)*F(n) * F(n)*F(n+1). Thus a(n) equals the integer part of the square root of the right hand side of the given equation. - Kenneth J Ramsey, Dec 19 2006
Contribution from Johannes W. Meijer, Sep 22 2010: (Start)
The a(n) represent several triangle sums of the Golden Triangle A180662: Kn11 (terms doubled), Kn12(n+1) (terms doubled), Kn4, Ca1 (terms tripled), Ca4, Gi1 (terms quadrupled) and Gi4. See A180662 for the definitions of these sums.
(End)
Define a 2 X (n+1) matrix with elements T(r,0)=A000032(r) and T(r,1) = Fibonacci(r), r=0,1,..,n. The matrix times its transposed is a 2 X 2 matrix with one diagonal element A001654(n+1), the other A216243(n), and A027941(n+1) on both outer diagonals. The determinant of this 2 X 2 matrix is 4*a(n). Example: For n=3 the matrix is 2 X 4 with rows 2 1 3 4; 0 1 1 2 to give as a product the 2 X 2 matrix with rows 30 12; 12 6 and determinant 180-144 = 36 =4*a(3). - J. M. Bergot, Feb 13 2013
a(n+1) is equal to the number of ternary strings of length n without any substring of the form 0x1, where x is in {0,1,2}. - John M. Campbell, Apr 03 2016

Crossrefs

Odd terms of A097083.
Partial sums of A001654.

Programs

  • GAP
    a:=[0,1,3,9];; for n in [5..30] do a[n]:=3*a[n-1]-3*a[n-3]+a[n-4]; od; a; # G. C. Greubel, Jan 09 2019
  • Magma
    m:=30; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!( x/((1-x^2)*(1-3*x+x^2)) )); // G. C. Greubel, Jan 09 2019
    
  • Mathematica
    Table[ Sum[ Fibonacci[k]*Fibonacci[k + 1], {k, n} ], {n, 0, 30}]
    f[n_] := Floor[GoldenRatio^(2 n + 2)/5]; Array[f, 28, 0] (* Robert G. Wilson v, Oct 25 2001 *)
    a[0]= 0; a[1]= 1; a[2]= 3; a[3]= 9; a[n_]:= a[n]= 3a[n-1] - 3a[n-3] + a[n-4]; Table[a[n], {n, 0, 20}] (* Vladimir Reshetnikov, Oct 28 2015 *)
  • PARI
    a(n)=if(n<0,0,fibonacci(n+1)^2-1+n%2)
    
  • PARI
    { for (n=0, 200, a=fibonacci(n+1)^2 - 1 + n%2; write("b064831.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 27 2009
    
  • PARI
    my(x='x+O('x^30)); concat([0], Vec(x/((1-x^2)*(1-3*x+x^2)))) \\ G. C. Greubel, Jan 09 2019
    
  • Sage
    (x/((1-x^2)*(1-3*x+x^2))).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Jan 09 2019
    

Formula

a(n) = F(n+1)^2 - 1 if n is even, or F(n+1)^2 if n is odd.
a(n) = A005313(n+1) - n.
G.f.: x/((1-x^2)*(1-3*x+x^2)). - N. J. A. Sloane Jul 15 2002
a(n) = Sum_{k=0..floor(n/2)} U(n-2k-1, 3/2). - Paul Barry, Nov 15 2003
Let M_n denote the n X n Hankel matrix M_n(i, j)=F(i+j-1) where F = A000045 is Fibonacci numbers, then the characteristic polynomial of M_n is x^n - F(2n)x^(n-1) + a(n-1)x^(n-2) . - Michael Somos, Nov 14 2002
a(n) = a(n-1) + A001654(n) with a(0)=0. (Partial sums of A001654). - Johannes W. Meijer, Sep 22 2010
a(n) = floor(phi^(2*n+2)/5), where phi =(1+sqrt(5))/2. - Gary Detlefs Mar 12 2011
a(n) = (A027941(n) + A001654(n))/2, n>=0. - Wolfdieter Lang, Jul 23 2012
a(n) = A005248(n+1)/5 -1/2 -(-1)^n/10. - R. J. Mathar, Feb 21 2013
Recurrence: a(0) = 0, a(1) = 1, a(2) = 3, a(3) = 9, a(n) = 3*a(n-1) - 3*a(n-3) + a(n-4). - Vladimir Reshetnikov, Oct 28 2015
a(n) = Sum_{i=0..n} (n+1-i)*Fibonacci(i)^2. - Bruno Berselli, Feb 20 2017

Extensions

More terms from Robert G. Wilson v, Oct 25 2001

A023039 a(n) = 18*a(n-1) - a(n-2).

Original entry on oeis.org

1, 9, 161, 2889, 51841, 930249, 16692641, 299537289, 5374978561, 96450076809, 1730726404001, 31056625195209, 557288527109761, 10000136862780489, 179445175002939041, 3220013013190122249, 57780789062419261441
Offset: 0

Keywords

Comments

The primitive Heronian triangle 3*a(n) +- 2, 4*a(n) has the latter side cut into 2*a(n) +- 3 by the corresponding altitude and has area 10*a(n)*A060645(n). - Lekraj Beedassy, Jun 25 2002
Chebyshev polynomials T(n,x) evaluated at x=9.
{a(n)} gives all (unsigned, integer) solutions of Pell equation a(n)^2 - 80*b(n)^2 = +1 with b(n) = A049660(n), n >= 0.
{a(n)} gives all possible solutions for x in Pell equation x^2 - D*y^2 = 1 for D=5, D=20 and D=80. The corresponding values for y are A060645 (D=5), A207832 (D=20) and A049660 (D=80). - Herbert Kociemba, Jun 05 2022
Also gives solutions to the equation x^2 - 1 = floor(x*r*floor(x/r)) where r=sqrt(5). - Benoit Cloitre, Feb 14 2004
Appears to give all solutions > 1 to the equation: x^2 = ceiling(x*r*floor(x/r)) where r=sqrt(5). - Benoit Cloitre, Feb 24 2004
For all terms x of the sequence, 5*x^2 - 5 is a square, A004292(n)^2.
The a(n) are the x-values in the nonnegative integer solutions of x^2 - 5y^2 = 1, see A060645(n) for the corresponding y-values. - Sture Sjöstedt, Nov 29 2011
Rightmost digits alternate repeatedly: 1 and 9 in fact, a(2) = 18*9 - 1 == 1 (mod 10); a(3) = 18*1 - 9 == 9 (mod 10) therefore a(2n) == 1 (mod 10), a(2n+1) == 9 (mod 10). - Carmine Suriano, Oct 03 2013

Examples

			G.f. = 1 + 9*x + 161*x^2 + 2889*x^3 + 51841*x4 + 930249*x^5 + 16692641*x^6 + ...
		

Crossrefs

Row 2 of array A188645.
Row 4 of A322790.

Programs

  • Magma
    I:=[1, 9]; [n le 2 select I[n] else 18*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Feb 13 2012
    
  • Maple
    a := n -> hypergeom([n, -n], [1/2], -4):
    seq(simplify(a(n)), n=0..16); # Peter Luschny, Jul 26 2020
  • Mathematica
    LinearRecurrence[{18, -1}, {1, 9}, 50] (* Sture Sjöstedt, Nov 29 2011 *)
    CoefficientList[Series[(1-9*x)/(1-18*x+x^2), {x, 0, 50}], x] (* G. C. Greubel, Dec 19 2017 *)
  • PARI
    {a(n) = fibonacci(6*n) / 2 + fibonacci(6*n - 1)}; /* Michael Somos, Aug 11 2009 */
    
  • PARI
    x='x+O('x^30); Vec((1-9*x)/(1-18*x+x^2)) \\ G. C. Greubel, Dec 19 2017

Formula

a(n) ~ (1/2)*(sqrt(5) + 2)^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
Limit_{n->infinity} a(n)/a(n-1) = phi^6 = 9 + 4*sqrt(5). - Gregory V. Richardson, Oct 13 2002
a(n) = T(n, 9) = (S(n, 18) - S(n-2, 18))/2, with S(n, x) := U(n, x/2) and T(n, x), resp. U(n, x), are Chebyshev's polynomials of the first, resp. second, kind. See A053120 and A049310. S(-2, x) := -1, S(-1, x) := 0, S(n, 18)=A049660(n+1).
a(n) = sqrt(80*A049660(n)^2 + 1) (cf. Richardson comment).
a(n) = ((9 + 4*sqrt(5))^n + (9 - 4*sqrt(5))^n)/2.
G.f.: (1 - 9*x)/(1 - 18*x + x^2).
a(n) = cosh(2*n*arcsinh(2)). - Herbert Kociemba, Apr 24 2008
a(n) = A001077(2*n). - Michael Somos, Aug 11 2009
From Johannes W. Meijer, Jul 01 2010: (Start)
a(n) = 2*A167808(6*n+1) - A167808(6*n+3).
Limit_{k->infinity} a(n+k)/a(k) = a(n) + A060645(n)*sqrt(5).
Limit_{n->infinity} a(n)/A060645(n) = sqrt(5).
(End)
a(n) = (1/2)*A087215(n) = (1/2)*(sqrt(5) + 2)^(2*n) + (1/2)*(sqrt(5) - 2)^(2*n).
Sum_{n >= 1} 1/( a(n) - 5/a(n) ) = 1/8. Compare with A005248, A002878 and A075796. - Peter Bala, Nov 29 2013
a(n) = 2*A115032(n-1) - 1 = S(n, 18) - 9*S(n-1, 18), with A115032(-1) = 1, and see the above formula with S(n, 18) using its recurrence. - Wolfdieter Lang, Aug 22 2014
a(n) = A128052(3n). - A.H.M. Smeets, Oct 02 2017
a(n) = A049660(n+1) - 9*A049660(n). - R. J. Mathar, May 24 2018
a(n) = hypergeom([n, -n], [1/2], -4). - Peter Luschny, Jul 26 2020
a(n) = L(6*n)/2 for L(n) the Lucas sequence A000032(n). - Greg Dresden, Dec 07 2021
a(n) = cosh(6*n*arccsch(2)). - Peter Luschny, May 25 2022

Extensions

Chebyshev and Pell comments from Wolfdieter Lang, Nov 08 2002
Sture Sjöstedt's comment corrected and reformulated by Wolfdieter Lang, Aug 24 2014

A003501 a(n) = 5*a(n-1) - a(n-2), with a(0) = 2, a(1) = 5.

Original entry on oeis.org

2, 5, 23, 110, 527, 2525, 12098, 57965, 277727, 1330670, 6375623, 30547445, 146361602, 701260565, 3359941223, 16098445550, 77132286527, 369562987085, 1770682648898, 8483850257405, 40648568638127, 194758992933230, 933146396028023, 4470972987206885
Offset: 0

Keywords

Comments

Positive values of x satisfying x^2 - 21*y^2 = 4; values of y are in A004254. - Wolfdieter Lang, Nov 29 2002
Except for the first term, positive values of x (or y) satisfying x^2 - 5xy + y^2 + 21 = 0. - Colin Barker, Feb 08 2014

Examples

			G.f. = 2 + 5*x + 23*x^2 + 110*x^3 + 527*x^4 + 2525*x^5 + ... - _Michael Somos_, Oct 25 2022
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • GAP
    a:=[2,5];; for n in [4..30] do a[n]:=5*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Jan 16 2020
  • Magma
    I:=[2,5]; [n le 2 select I[n] else 5*Self(n-1) -Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 16 2020
    
  • Magma
    R:=PowerSeriesRing(Integers(), 25); Coefficients(R!((2-5*x)/(1-5*x+x^2))); // Marius A. Burtea, Jan 16 2020
    
  • Maple
    seq( simplify(2*ChebyshevT(n, 5/2)), n=0..30); # G. C. Greubel, Jan 16 2020
  • Mathematica
    a[0]=2; a[1]=5; a[n_]:= 5a[n-1] -a[n-2]; Table[a[n], {n,0,30}] (* Robert G. Wilson v, Jan 30 2004 *)
    LinearRecurrence[{5,-1},{2,5},30] (* Harvey P. Dale, May 12 2019 *)
    2*ChebyshevT[Range[0, 30], 5/2] (* G. C. Greubel, Jan 16 2020 *)
    a[ n_] := LucasL[n, 5*I]/I^n; (* Michael Somos, Oct 25 2022 *)
  • PARI
    {a(n) = subst(poltchebi(n),x,5/2)*2};
    
  • PARI
    {a(n) = polchebyshev(n,1,5/2)*2 }; /* Michael Somos, Oct 25 2022 */
    
  • Sage
    [lucas_number2(n,5,1) for n in range(37)] # Zerinvary Lajos, Jun 25 2008
    

Formula

a(n) = 5*S(n-1, 5) - 2*S(n-2, 5) = S(n, 5) - S(n-2, 5) = 2*T(n, 5/2), with S(n, x)=U(n, x/2), S(-1, x)=0, S(-2, x)=-1. U(n, x), resp. T(n, x), are Chebyshev's polynomials of the second, resp. first, kind. S(n-1, 5) = A004254(n), n>=0.
G.f.: (2-5*x)/(1-5*x+x^2). - Simon Plouffe in his 1992 dissertation.
a(n) ~ (1/2*(5 + sqrt(21)))^n. - Joe Keane (jgk(AT)jgk.org), May 16 2002
a(n) = ap^n + am^n, with ap=(5+sqrt(21))/2 and am=(5-sqrt(21))/2.
a(n) = sqrt(4 + 21*A004254(n)^2).
From Peter Bala, Jan 06 2013: (Start)
Let F(x) = Product_{n=0..inf} (1 + x^(4*n+1))/(1 + x^(4*n+3)). Let alpha = 1/2*(5 - sqrt(21)). This sequence gives the simple continued fraction expansion of 1 + F(alpha) = 2.19827 65373 95327 17782 ... = 2 + 1/(5 + 1/(23 + 1/(110 + ...))).
Also F(-alpha) = 0.79824 49142 28050 93561 ... has the continued fraction representation 1 - 1/(5 - 1/(23 - 1/(110 - ...))) and the simple continued fraction expansion 1/(1 + 1/((5-2) + 1/(1 + 1/((23-2) + 1/(1 + 1/((110-2) + 1/(1 + ...))))))).
F(alpha)*F(-alpha) has the simple continued fraction expansion 1/(1 + 1/((5^2-4) + 1/(1 + 1/((23^2-4) + 1/(1 + 1/((110^2-4) + 1/(1 + ...))))))).
(End)
a(n) = (A217787(k+3n) + A217787(k-3n))/A217787(k) for k>=3n. - Bruno Berselli, Mar 25 2013

Extensions

Chebyshev comments from Wolfdieter Lang, Oct 31 2002

A027011 Triangular array T read by rows: T(n,k) = t(n, 2k+1) for 0 <= k <= floor((2n-1)/2), t given by A027960, n >= 0.

Original entry on oeis.org

3, 3, 4, 3, 7, 5, 3, 7, 15, 6, 3, 7, 18, 28, 7, 3, 7, 18, 44, 47, 8, 3, 7, 18, 47, 98, 73, 9, 3, 7, 18, 47, 120, 199, 107, 10, 3, 7, 18, 47, 123, 291, 373, 150, 11, 3, 7, 18, 47, 123, 319, 661, 654, 203, 12, 3, 7, 18, 47, 123, 322, 806, 1404, 1085, 267, 13, 3, 7, 18, 47
Offset: 1

Keywords

Comments

Right-edge columns are polynomials approximating Lucas(2n).

Examples

			                                          3
                                     3,   4
                                3,   7,   5
                           3,   7,  15,   6
                      3,   7,  18,  28,   7
                 3,   7,  18,  44,  47,   8
            3,   7,  18,  47,  98,  73,   9
       3,   7,  18,  47, 120, 199, 107,  10
  3,   7,  18,  47, 123, 291, 373, 150,  11
		

Crossrefs

This is a bisection of the "Lucas array " A027960, see A026998 for the other bisection.
Right-edge columns include A027965, A027967, A027969, A027971.
An earlier version of this entry had (unjustifiably) each row starting with 1.

Programs

  • Maple
    t:=proc(n,k)option remember:if(k=0 or k=2*n)then return 1:elif(k=1)then return 3:else return t(n-1,k-2) + t(n-1,k-1):fi:end:
    T:=proc(n,k)return t(n,2*k+1):end:
    for n from 0 to 8 do for k from 0 to floor((2*n-1)/2) do print(T(n,k));od:od: # Nathaniel Johnston, Apr 18 2011
  • Mathematica
    t[n_, k_] := t[n, k] = If[k == 0 || k == 2*n, 1, If[k == 1, 3, t[n-1, k-2] + t[n-1, k-1]]]; T[n_, k_] := t[n, 2*k+1]; Table[T[n, k], {n, 1, 12}, {k, 0, (2*n-1)/2}] // Flatten (* Jean-François Alcover, Nov 18 2013, after Nathaniel Johnston *)

Formula

T(n, k) = Lucas(2n) = A005248(n) for 2k+1 <= n, otherwise the (2n-2k+1)-th coefficient of the power series for (1+2x)/((1-x-x^2)(1-x)^(2k-n+1)).

Extensions

Edited by Ralf Stephan, May 05 2005

A127677 Scaled coefficient table for Chebyshev polynomials 2*T(2*n, sqrt(x)/2) (increasing even scaled powers, without zero entries).

Original entry on oeis.org

2, -2, 1, 2, -4, 1, -2, 9, -6, 1, 2, -16, 20, -8, 1, -2, 25, -50, 35, -10, 1, 2, -36, 105, -112, 54, -12, 1, -2, 49, -196, 294, -210, 77, -14, 1, 2, -64, 336, -672, 660, -352, 104, -16, 1, -2, 81, -540, 1386, -1782, 1287, -546, 135, -18, 1, 2, -100, 825, -2640, 4290, -4004, 2275, -800, 170, -20, 1
Offset: 0

Author

Wolfdieter Lang, Mar 07 2007

Keywords

Comments

2*T(2*n,x) = Sum_{m=0..n} a(n,m)*(2*x)^(2*m).
Closely related to A284982, which has opposite signs and rows begin with 0 of alternating signs instead of +/2. - Eric W. Weisstein, Apr 07 2017
Bisection triangle of A127672 (without zero entries, even part). The odd part is ((-1)^(n-m))*A111125(n,m).
If the leading 2 is replaced by a 1 we get the essentially identical sequence A110162. - N. J. A. Sloane, Jun 09 2007
Also row n gives coefficients of characteristic polynomial of the Cartan matrix for the root system B_n (or, equally, C_n). - Roger L. Bagula, May 23 2007
From Wolfdieter Lang, Oct 04 2013: (Start)
This triangle a(n,m) is used to express the length ratio side/R given by s(4*n+2) = 2*sin(Pi/(4*n+2)) = 2*cos(2*n*Pi/(4*n+2)) in a regular (4*n+2)-gon, inscribed in a circle with radius R, in terms of rho(4*n+2) = 2*cos(Pi/4*n+2), the length ratio of (the smallest diagonal)/side (for n=2 there is no such diagonal).
s(4*n+2) = Sum_{m=0..n}a(n,m)*rho(4*n+2)^(2*m). This formula is needed to show that the total sum of all length ratios in a (4*n+2)-gon is an integer in the algebraic number field Q(rho(4*n+2)). Note that rho(4*n+2) has degree delta(4*n+2) = A055034(4*n+2). Therefore one has to take s(4*n+2) modulo C(4*n+2, x=rho(4*n+2)), the minimal polynomial of rho(4*n+2) (see A187360). Thanks go to Seppo Mustonen for asking me to look into this problem. See ((-1)^(n-m))*A111125(n,m) for the (4*n)-gon situation. (End)

Examples

			The triangle a(n,m) starts:
n\m  0    1    2     3     4     5     6     7    8   9  10 ...
0:   2
1:  -2    1
2:   2   -4    1
3:  -2    9   -6     1
4:   2  -16   20    -8     1
5:  -2   25  -50    35   -10     1
6:   2  -36  105  -112    54   -12     1
7:  -2   49 -196   294  -210    77   -14     1
8:   2  -64  336  -672   660  -352   104   -16    1
9:  -2   81 -540  1386 -1782  1287  -546   135  -18   1
10:  2 -100  825 -2640  4290 -4004  2275  -800  170 -20  1
... Reformatted and extended by _Wolfdieter Lang_, Nov 21 2012.
n=3: [-2,9,-6,1] stands for -2*1 + 9*(2*x)^2 -6*(2*x)^4 +1*(2*x)^6 = 2*(1+18*x^2-48*x^4+32*x^6) = 2*T(6,x).
(4*n+2)-gon side/radius s(4*n+2) as polynomial in rho(4*n+2) = smallest diagonal/side: n=0: s(2) = 2 (rho(2)=0); n=1: s(6) = -2 + rho(6)^2 = -2 + 3 = 1, (C(6,x) = x^2 - 3); n=2: s(10) = 2 - 4*rho(10)^2 + 1*rho(10)^4 = 2 - 4*rho(10)^2 + (5*rho(10)^2 - 5) = -3 + rho(10)^2, (C(10,x) = x^4 - 5*x^2 + 5). - _Wolfdieter Lang_, Oct 04 2013
		

References

  • R. N. Cahn, Semi-Simple Lie Algebras and Their Representations, Dover, NY, 2006, ISBN 0-486-44999-8, p. 62
  • Sigurdur Helgasson,Differential Geometry, Lie Groups and Symmetric Spaces, Graduate Studies in Mathematics, volume 34. A. M. S. :ISBN 0-8218-2848-7, 1978,p. 463.

Crossrefs

Cf. A284982 (opposite signs and rows begin with 0).
Row sums (signed): -A061347(n+3) for n>=0.
Row sums (unsigned): A005248(n) = L(2*n), where L=Lucas.

Programs

  • Mathematica
    T[n_, m_, d_] := If[ n == m, 2, If[n == d && m == d - 1, -2, If[(n == m - 1 || n == m + 1), -1, 0]]] M[d_] := Table[T[n, m, d], {n, 1, d}, {m, 1, d}] a = Join[M[1], Table[CoefficientList[CharacteristicPolynomial[M[d], x], x], {d, 1, 10} ]] (* Roger L. Bagula, May 23 2007 *)
    CoefficientList[2 ChebyshevT[2 Range[0, 10], Sqrt[x]/2], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
    CoefficientList[Table[(-1)^n LucasL[2 n, Sqrt[-x]], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
  • PARI
    a(n,m) = {if(n>=2, -2*a(n-1,m)+a(n-1,m-1)-a(n-2,m), if(n==0, if(m!=0,0,2), if(m==0,-2, if(m==1,1,0))))};
    for(n=0,10,for(m=0,n,print1(a(n,m),", "))) \\ Hugo Pfoertner, Jul 19 2020

Formula

a(n,m) = 0 if n < m; a(n,0) = 2*(-1)^n; a(n,m) = ((-1)^(n+m))*n*binomial(n+m-1, 2*m-1)/m.
a(n,m) = 0 if n < m, a(0,0) = 2, a(n,m) = (-1)^(n-m)*(2*n/(n+m))*binomial(n+m, n-m), n >= 1. From Waring's formula applied to Chebyshev's T-polynomials. See also A110162. - Wolfdieter Lang, Nov 21 2012
The o.g.f. Sum_{n>=0} p(n,x)*z^n, n>=0, for the row polynomials p(n,x) := Sum_{m=0..n} a(n,m)*x^m is (2 + z*(2-x))/((z+1)^2 - z*x). Here p(n,x) = R(2*n,sqrt(x)) := 2*T(2*n,sqrt(x)/2) with Chebyshev's T-polynomials. For the R-polynomials see A127672. - Wolfdieter Lang, Nov 28 2012
From Tom Copeland, Nov 07 2015: (Start)
A logarithmic generator is 2*(1-log(1+x))-log(1-t*x/(1+x)^2) = 2 - log(1+(2-t)*x+x^2) = 2 + (-2 + t)*x + (2 - 4*t + t^2) x^2/2 + (-2 + 9*t - 6*t^2 + t^3) x^3/3 + ..., so a number of relations to the Faber polynomials of A263916 hold with p(0,x) = 2:
1) p(n,x) = F(n,(2-x),1,0,0,..)
2) p(n,x) = (-1)^n 2 + F(n,-x,2x,-3x,...,(-1)^n n*x)
3) p(n,x) = (-1)^n [2 + F(n,x,2x,3x,...,n*x)].
The unsigned array contains the partial sums of A111125 modified by appending a column of zeros, except for an initial two, to A111125. Then the difference of consecutive rows of unsigned A127677, further modified by appending an initial rows of zeros, generates the modified A111125. Cf. A208513 and A034807.
For relations among the characteristic polynomials of Cartan matrices of the Coxeter root groups, Chebyshev polynomials, cyclotomic polynomials, and the polynomials of this entry, see Damianou (p. 12, 20, and 21) and Damianou and Evripidou (p. 7).
See A111125 for a relation to the squares of the odd row polynomials here with the constant removed.
p(n,x)^2 = 2 + p(2*n,x). See also A127672. (End)
a(n,m) = -2*a(n-1,m) + a(n-1,m-1) - a(n-2,m) for n >= 2 with initial conditions a(0,0) = 2, a(1,0) = -2, a(1,1) = 1, a(0,m) = 0 for m != 0, a(1,m) = 0 for m != 0,1. - William P. Orrick, Jun 09 2020
p(n,x) = (x-2)*p(n-1,x) - p(n-2,x) for n >= 2. - William P. Orrick, Jun 09 2020

Extensions

Definition corrected by Eric W. Weisstein, Apr 06 2017

A164102 Decimal expansion of 2*Pi^2.

Original entry on oeis.org

1, 9, 7, 3, 9, 2, 0, 8, 8, 0, 2, 1, 7, 8, 7, 1, 7, 2, 3, 7, 6, 6, 8, 9, 8, 1, 9, 9, 9, 7, 5, 2, 3, 0, 2, 2, 7, 0, 6, 2, 7, 3, 9, 8, 8, 1, 4, 4, 8, 1, 5, 8, 1, 2, 5, 2, 8, 2, 6, 6, 9, 8, 7, 5, 2, 4, 4, 0, 0, 8, 9, 6, 4, 4, 8, 3, 8, 4, 1, 0, 4, 8, 6, 0, 0, 3, 5, 4, 6, 8, 0, 7, 4, 3, 7, 1, 0, 4, 4, 6, 3, 6, 4, 8, 0
Offset: 2

Author

R. J. Mathar, Aug 10 2009

Keywords

Comments

Surface area of the 4-dimensional unit sphere. The volume of the 4-dimensional unit sphere is a fourth of this, A102753.
Also decimal expansion of Pi^2/5 = 1.973920..., with offset 1. - Omar E. Pol, Oct 04 2011

Examples

			19.739208802178717237668981...
		

References

  • L. A. Santalo, Integral Geometry and Geometric Probability, Addison-Wesley, 1976, see p. 15.

Programs

Formula

Equals 2*A002388 = 4*A102753.
Pi^2/5 = Sum_{k>=1} Lucas(2*k)/(k^2*binomial(2*k,k)) = Sum_{k>=1} A005248(k)/A002736(k) (Seiffert, 1991). - Amiram Eldar, Jan 17 2022

A020876 a(n) = ((5+sqrt(5))/2)^n + ((5-sqrt(5))/2)^n.

Original entry on oeis.org

2, 5, 15, 50, 175, 625, 2250, 8125, 29375, 106250, 384375, 1390625, 5031250, 18203125, 65859375, 238281250, 862109375, 3119140625, 11285156250, 40830078125, 147724609375, 534472656250, 1933740234375, 6996337890625, 25312988281250, 91583251953125
Offset: 0

Keywords

Comments

Number of no-leaf edge-subgraphs in Moebius ladder M_n.

Examples

			G.f. = 2 + 5*x + 15*x^2 + 50*x^3 + 175*x^4 + 625*x^5 + 2250*x^6 + ...
		

Crossrefs

Appears in A109106. - Johannes W. Meijer, Jul 01 2010

Programs

  • Magma
    [Floor(((5+Sqrt(5))/2)^n+((5-Sqrt(5))/2)^n): n in [0..30]]; // Vincenzo Librandi, Aug 08 2014
  • Maple
    G:=(x,n)-> cos(x)^n+cos(3*x)^n:
    seq(simplify(4^n*G(Pi/10,2*n)), n=0..22); # Gary Detlefs, Dec 05 2010
  • Mathematica
    Table[Sum[LucasL[2*i] Binomial[n, i], {i, 0, n}], {n, 0, 50}] (* T. D. Noe, Sep 10 2011 *)
    CoefficientList[Series[(2 - 5 x)/(1 - 5 x + 5 x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Aug 08 2014 *)
    LinearRecurrence[{5,-5},{2,5},30] (* Harvey P. Dale, Mar 13 2016 *)
  • Sage
    [lucas_number2(n,5,5) for n in range(0,24)] # Zerinvary Lajos, Jul 08 2008
    

Formula

Also, a(n) = (sqrt(5)*phi)^n + (sqrt(5)/phi)^n, where phi = golden ratio. - N. J. A. Sloane, Aug 08 2014
Let S(n, m)=sum(k=0, n, binomial(n, k)*fibonacci(m*k)), then for n>0 a(n)= S(2*n, 2)/S(n, 2). - Benoit Cloitre, Oct 22 2003
From R. J. Mathar, Feb 06 2010: (Start)
a(n)= 5*a(n-1) - 5*a(n-2).
G.f.: (2-5*x)/(1-5*x+5*x^2). (End)
From Johannes W. Meijer, Jul 01 2010: (Start)
Lim_{k->infinity} a(n+k)/a(k) = (A020876(n) + A093131(n)*sqrt(5))/2.
Lim_{k->infinity} A020876(n)/A093131(n) = sqrt(5). (End)
Binomial transform of A005248. - Carl Najafi, Sep 10 2011
a(n) = 2*A030191(n) - 5*A030191(n-1). - R. J. Mathar, Mar 02 2012
From Kai Wang, Dec 22 2019: (Start)
a((2*m+1)*k)/a(k) = Sum_{i=0..m-1} (-1)^(i*(k+1))*a(2*(m-i)*k) + 5^(m*k).
A093131(m+r)*A093131(n+s) + A093131(m+s)*A093131(n+r) = (2*a(m+n+r+s) - 5^(n+s)*a(m-n)*a(r-s))/5.
a(m+r)*a(n+s) - a(m+s)*a(n+r) = 5^(n+s+1)*A093131(m-n)*A093131(r-s).
a(m+r)*a(n+s) + a(m+s)*a(n+r) = 2*a(m+n+r+s) + 5^(n+s)*a(m-n)*a(r-s).
a(m+r)*a(n+s) - 5*A093131(m+s)*A093131(n+r) = 5^(n+s)*a(m-n)*a(r-s).
a(m+r)*a(n+s) + 5*A093131(m+s)*A093131(n+r) = 2*a(m+n+r+s)+ 5^(n+s+1)*A093131(m-n)*A093131(r-s).
A093131(m-n) = (A093131(m)*a(n) - a(m)*A093131(n))/(2*5^n).
A093131(m+n) = (A093131(m)*a(n) + a(m)*A093131(n))/2.
a(n)^2 - a(n+1)*a(n-1) = -5^n.
a(n)^2 - a(n+r)*a(n-r) = -5^(n-r+1)*A093131(r)^2.
a(m)*a(n+1) - a(m+1)*a(n) = -5^(n+1)*A093131(m-n).
a(m+n) - 5^(n)*a(m-n) = 5*A093131(m)*A093131(n).
a(m+n) + 5^(n)*a(m-n) = a(m)*a(n).
a(m-n) = (a(m)*a(n) - 5*A093131(m)*A093131(n))/(2*5^n).
a(m+n) = (a(m)*a(n) + 5*A093131(m)*A093131(n))/2. (End)
E.g.f.: 2*exp(5*x/2)*cosh(sqrt(5)*x/2). - Stefano Spezia, Dec 27 2019

Extensions

Definition simplified by N. J. A. Sloane, Aug 08 2014
Previous Showing 21-30 of 152 results. Next