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.

A034807 Triangle T(n,k) of coefficients of Lucas (or Cardan) polynomials.

Original entry on oeis.org

2, 1, 1, 2, 1, 3, 1, 4, 2, 1, 5, 5, 1, 6, 9, 2, 1, 7, 14, 7, 1, 8, 20, 16, 2, 1, 9, 27, 30, 9, 1, 10, 35, 50, 25, 2, 1, 11, 44, 77, 55, 11, 1, 12, 54, 112, 105, 36, 2, 1, 13, 65, 156, 182, 91, 13, 1, 14, 77, 210, 294, 196, 49, 2, 1, 15, 90, 275, 450, 378, 140, 15, 1, 16, 104
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

These polynomials arise in the following setup. Suppose G and H are power series satisfying G + H = G*H = 1/x. Then G^n + H^n = (1/x^n)*L_n(-x).
Apart from signs, triangle of coefficients when 2*cos(nt) is expanded in terms of x = 2*cos(t). For example, 2*cos(2t) = x^2 - 2, 2*cos(3t) = x^3 - 3x and 2*cos(4t) = x^4 - 4x^2 + 2. - Anthony C Robin, Jun 02 2004
Triangle of coefficients of expansion of Z_{nk} in terms of Z_k.
Row n has 1 + floor(n/2) terms. - Emeric Deutsch, Dec 25 2004
T(n,k) = number of k-matchings of the cycle C_n (n > 1). Example: T(6,2)=9 because the 2-matchings of the hexagon with edges a, b, c, d, e, f are ac, ad, ae, bd, be, bf, ce, cf and df. - Emeric Deutsch, Dec 25 2004
An example for the first comment: G=c(x), H=1/(x*c(x)) with c(x) the o.g.f. Catalan numbers A000108: (x*c(x))^n + (1/c(x))^n = L(n,-x)= Sum_{k=0..floor(n/2)} T(n,k)*(-x)^k.
This triangle also supplies the absolute values of the coefficients in the multiplication formulas for the Lucas numbers A000032.
From L. Edson Jeffery, Mar 19 2011: (Start)
This sequence is related to rhombus substitution tilings. A signed version of it (see A132460), formed as a triangle with interlaced zeros extending each row to n terms, begins as
{2}
{1, 0}
{1, 0, -2}
{1, 0, -3, 0}
{1, 0, -4, 0, 2}
{1, 0, -5, 0, 5, 0}
....
For the n X n tridiagonal unit-primitive matrix G_(n,1) (n >= 2) (see the L. E. Jeffery link below), defined by
G_(n,1) =
(0 1 0 ... 0)
(1 0 1 0 ... 0)
(0 1 0 1 0 ... 0)
...
(0 ... 0 1 0 1)
(0 ... 0 2 0),
Row n (i.e., {T(n,k)}, k=0..n) of the signed table gives the coefficients of its characteristic function: c_n(x) = Sum_{k=0..n} T(n,k)*x^(n-k) = 0. For example, let n=3. Then
G_(3,1) =
(0 1 0)
(1 0 1)
(0 2 0),
and row 3 of the table is {1,0,-3,0}. Hence c_3(x) = x^3 - 3*x = 0. G_(n,1) has n distinct eigenvalues (the solutions of c_n(x) = 0), given by w_j = 2*cos((2*j-1)*Pi/(2*n)), j=1..n. (End)
For n > 0, T(n,k) is the number of k-subsets of {1,2,...,n} which contain neither consecutive integers nor both 1 and n. Equivalently, T(n,k) is the number of k-subsets without neighbors of a set of n points on a circle. - José H. Nieto S., Jan 17 2012
With the first column omitted, this gives A157000. - Philippe Deléham, Mar 17 2013
The number of necklaces of k black and n - k white beads with no adjacent black beads (Kaplansky 1943). Coefficients of the Dickson polynomials D(n,x,-a). - Peter Bala, Mar 09 2014
From Tom Copeland, Nov 07 2015: (Start)
This triangular array is composed of interleaved rows of reversed, unsigned A127677 (cf. A156308, A217476, A263916) and reversed A111125 (cf. A127672).
See also A113279 for another connection to symmetric and Faber polynomials.
The difference of consecutive rows gives the previous row shifted.
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). (End)
Diagonals are related to multiplicities of eigenvalues of the Laplacian on hyperspheres through A029635. - Tom Copeland, Jan 10 2016
For n>=3, also the independence and matching polynomials of the n-cycle graph C_n. See also A284966. - Eric W. Weisstein, Apr 06 2017
Apparently, with the rows aerated and then the 2s on the diagonal removed, this matrix becomes the reverse, or mirror, of unsigned A117179. See also A114525 - Tom Copeland, May 30 2017
Briggs's (1633) table with an additional column of 2s on the right can be used to generate this table. See p. 69 of the Newton reference. - Tom Copeland, Jun 03 2017
From Liam Solus, Aug 23 2018: (Start)
For n>3 and k>0, T(n,k) equals the number of Markov equivalence classes with skeleton the cycle on n nodes having exactly k immoralities. See Theorem 2.1 of the article by A. Radhakrishnan et al. below.
For n>2 odd and r = floor(n/2)-1, the n-th row is the coefficient vector of the Ehrhart h*-polynomial of the r-stable (n,2)-hypersimplex. See Theorem 4.14 in the article by B. Braun and L. Solus below.
(End)
Conjecture: If a(n) = H(a,b,c,d,n) is a second-order linear recurrence with constant coefficients defined as a(0) = a, a(1)= b, a(n) = c*a(n-1) + d*a(n-2) then a(m*n) = H(a, H(a,b,c,d,m), Sum_{k=0..floor(m/2)} T(m,k)*c^(m-2*k)*d^k, (-1)^(m+1)*d^m, n) (Wolfdieter Lang). - Gary Detlefs, Feb 06 2023
For the proof of the preceding conjecture see the Detlefs and Lang link. There also proofs for several properties of this table are found. - Wolfdieter Lang, Apr 25 2023
From Mohammed Yaseen, Nov 09 2024: (Start)
Let m - 1/m = x, then
m^2 + 1/m^2 = x^2 + 2,
m^3 - 1/m^3 = x^3 + 3*x,
m^4 + 1/m^4 = x^4 + 4*x^2 + 2,
m^5 - 1/m^5 = x^5 + 5*x^3 + 5*x,
m^6 + 1/m^6 = x^6 + 6*x^4 + 9*x^2 + 2,
m^7 - 1/m^7 = x^7 + 7*x^5 + 14*x^3 + 7*x, etc. (End)

Examples

			I have seen two versions of these polynomials: One version begins L_0 = 2, L_1 = 1, L_2 = 1 + 2*x, L_3 = 1 + 3*x, L_4 = 1 + 4*x + 2*x^2, L_5 = 1 + 5*x + 5*x^2, L_6 = 1 + 6*x + 9*x^2 + 2*x^3, L_7 = 1 + 7*x + 14*x^2 + 7*x^3, L_8 = 1 + 8*x + 20*x^2 + 16*x^3 + 2*x^4, L_9 = 1 + 9*x + 27*x^2 + 30*x^3 + 9*x^4, ...
The other version (probably the more official one) begins L_0(x) = 2, L_1(x) = x, L_2(x) = 2 + x^2, L_3(x) = 3*x + x^3, L_4(x) = 2 + 4*x^2 + x^4, L_5(x) = 5*x + 5*x^3 + x^5, L_6(x) = 2 + 9*x^2 + 6*x^4 + x^6, L_7(x) = 7*x + 14*x^3 + 7*x^5 + x^7, L_8(x) = 2 + 16*x^2 + 20*x^4 + 8*x^6 + x^8, L_9(x) = 9*x + 30*x^3 + 27*x^5 + 9*x^7 + x^9.
From _John Blythe Dobson_, Oct 11 2007: (Start)
Triangle begins:
  2;
  1;
  1,  2;
  1,  3;
  1,  4,  2;
  1,  5,  5;
  1,  6,  9,   2;
  1,  7, 14,   7;
  1,  8, 20,  16,   2;
  1,  9, 27,  30,   9;
  1, 10, 35,  50,  25,   2;
  1, 11, 44,  77,  55,  11;
  1, 12, 54, 112, 105,  36,   2;
  1, 13, 65, 156, 182,  91,  13;
  1, 14, 77, 210, 294, 196,  49,  2;
  1, 15, 90, 275, 450, 378, 140, 15;
(End)
From _Peter Bala_, Mar 20 2025: (Start)
Let S = x + y and M = -x*y. Then the triangle gives the coefficients when expressing the symmetric polynomial x^n + y^n as a polynomial in S and M. For example,
x^2 + y^2 = S^2 + 2*M; x^3 + y^3 = S^3 + 3*S*M; x^4 + y^4 = S^4 + 4*(S^2)*M + 2*M^2;
x^5 + y^5 = S^5 + 5*(S^3)*M + 5*S*M^2; x^6 + y^6 = S^6 + 6*(S^4)*M + 9*(S^2)*M^2 + 2*M^3. See Woko. In general x^n + y^n = 2*(-i)^n *(sqrt(M))^n * T(n, i*S/(2*sqrt(M))), where T(n, x) denotes the n-th Chebyshev polynomial of the first kind. (End)

References

  • A. Brousseau, Fibonacci and Related Number Theoretic Tables. Fibonacci Association, San Jose, CA, 1972, p. 148.
  • C. D. Godsil, Algebraic Combinatorics, Chapman and Hall, New York, 1993.
  • Thomas Koshy, Fibonacci and Lucas Numbers with Applications. New York, etc.: John Wiley & Sons, 2001. (Chapter 13, "Pascal-like Triangles," is devoted to the present triangle.)
  • The Royal Society Newton Tercentenary Celebrations, Cambridge Univ. Press, 1947.

Links

Crossrefs

Cf. A029635, A061896, A111125, A113279, A114525, A127672, A127677, A156308, A217476, A263916.
Cf. A117179.

Programs

  • Maple
    T:= proc(n,k) if n=0 and k=0 then 2 elif k>floor(n/2) then 0 else n*binomial(n-k,k)/(n-k) fi end: for n from 0 to 15 do seq(T(n,k), k=0..floor(n/2)) od; # yields sequence in triangular form # Emeric Deutsch, Dec 25 2004
  • Mathematica
    t[0, 0] = 2; t[n_, k_] := Binomial[n-k, k] + Binomial[n-k-1, k-1]; Table[t[n, k], {n, 0, 16}, {k, 0, Floor[n/2]}] // Flatten (* Jean-François Alcover, Dec 30 2013 *)
CoefficientList[Table[x^(n/2) LucasL[n, 1/Sqrt[x]], {n, 0, 15}], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
Table[Select[Reverse[CoefficientList[LucasL[n, x], x]], 0 < # &], {n, 0, 16}] // Flatten (* Robert G. Wilson v, May 03 2017 *)
CoefficientList[FunctionExpand @ Table[2 (-x)^(n/2) Cos[n ArcSec[2 Sqrt[-x]]], {n, 0, 15}], x] // Flatten (* Eric W. Weisstein, Apr 03 2018 *)
CoefficientList[Table[2 (-x)^(n/2) ChebyshevT[n, 1/(2 Sqrt[-x])], {n, 0, 15}], x] // Flatten (* Eric W. Weisstein, Apr 03 2018 *)
  • PARI
    {T(n, k) = if( k<0 || 2*k>n, 0, binomial(n-k, k) + binomial(n-k-1, k-1) + (n==0))}; /* Michael Somos, Jul 15 2003 */

Formula

Row sums = A000032. T(2n, n-1) = A000290(n), T(2n+1, n-1) = A000330(n), T(2n, n-2) = A002415(n). T(n, k) = A029635(n-k, k), if n>0. - Michael Somos, Apr 02 1999
Lucas polynomial coefficients: 1, -n, n*(n-3)/2!, -n*(n-4)*(n-5)/3!, n*(n-5)*(n-6)*(n-7)/4!, - n*(n-6)*(n-7)*(n-8)*(n-9)/5!, ... - Herb Conn and Gary W. Adamson, May 28 2003
G.f.: (2-x)/(1-x-x^2*y). - Vladeta Jovovic, May 31 2003
T(n, k) = T(n-1, k) + T(n-2, k-1), n>1. T(n, 0) = 1, n>0. T(n, k) = binomial(n-k, k) + binomial(n-k-1, k-1) = n*binomial(n-k-1, k-1)/k, 0 <= 2*k <= n except T(0, 0) = 2. - Michael Somos, Apr 02 1999
T(n,k) = (n*(n-1-k)!)/(k!*(n-2*k)!), n>0, k>=0. - Alexander Elkins (alexander_elkins(AT)hotmail.com), Jun 09 2007
O.g.f.: 2-(2xt+1)xt/(-t+xt+(xt)^2). (Cf. A113279.) - Tom Copeland, Nov 07 2015
T(n,k) = A011973(n-1,k) + A011973(n-3,k-1) = A011973(n,k) - A011973(n-4,k-2) except for T(0,0)=T(2,1)=2. - Xiangyu Chen, Dec 24 2020
L_n(x) = ((x+sqrt(x^2+4))/2)^n + (-((x+sqrt(x^2+4))/2))^(-n). See metallic means. - William Krier, Sep 01 2023

Extensions

Improved description, more terms, etc., from Michael Somos

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

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