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 17 results. Next

A127672 Monic integer version of Chebyshev T-polynomials (increasing powers).

Original entry on oeis.org

2, 0, 1, -2, 0, 1, 0, -3, 0, 1, 2, 0, -4, 0, 1, 0, 5, 0, -5, 0, 1, -2, 0, 9, 0, -6, 0, 1, 0, -7, 0, 14, 0, -7, 0, 1, 2, 0, -16, 0, 20, 0, -8, 0, 1, 0, 9, 0, -30, 0, 27, 0, -9, 0, 1, -2, 0, 25, 0, -50, 0, 35, 0, -10, 0, 1, 0, -11, 0, 55, 0, -77, 0, 44, 0, -11, 0, 1, 2, 0, -36, 0, 105, 0, -112, 0, 54, 0, -12, 0, 1, 0, 13, 0, -91
Offset: 0

Views

Author

Wolfdieter Lang, Mar 07 2007

Keywords

Comments

The row polynomials R(n,x) := Sum_{m=0..n} a(n,m)*x^m have been called Chebyshev C_n(x) polynomials in the Abramowitz-Stegun handbook, p. 778, 22.5.11 (see A049310 for the reference, and note that on p. 774 the S and C polynomials have been mixed up in older printings). - Wolfdieter Lang, Jun 03 2011
This is a signed version of triangle A114525.
The unsigned column sequences (without zeros) are, for m=1..11: A005408, A000290, A000330, A002415, A005585, A040977, A050486, A053347, A054333, A054334, A057788.
The row polynomials R(n,x) := Sum_{m=0..n} a(n,m)*x*m, give for n=2,3,...,floor(N/2) the positive zeros of the Chebyshev S(N-1,x)-polynomial (see A049310) in terms of its largest zero rho(N):= 2*cos(Pi/N) by putting x=rho(N). The order of the positive zeros is falling: n=1 corresponds to the largest zero rho(N) and n=floor(N/2) to the smallest positive zero. Example N=5: rho(5)=phi (golden section), R(2,phi)= phi^2-2 = phi-1, the second largest (and smallest) positive zero of S(4,x). - Wolfdieter Lang, Dec 01 2010
The row polynomial R(n,x), for n >= 1, factorizes into minimal polynomials of 2*cos(Pi/k), called C(k,x), with coefficients given in A187360, as follows.
R(n,x) = Product_{d|oddpart(n)} C(2*n/d,x)
= Product_{d|oddpart(n)} C(2^(k+1)*d,x),
with oddpart(n)=A000265(n), and 2^k is the largest power of 2 dividing n, where k=0,1,2,...
(Proof: R and C are monic, the degree on both sides coincides, and the zeros of R(n,x) appear all on the r.h.s.) - Wolfdieter Lang, Jul 31 2011 [Theorem 1B, eq. (43) in the W. Lang link. - Wolfdieter Lang, Apr 13 2018]
The zeros of the row polynomials R(n,x) are 2*cos(Pi*(2*k+1)/(2*n)), k=0,1, ..., n-1; n>=1 (from those of the Chebyshev T-polynomials). - Wolfdieter Lang, Sep 17 2011
The discriminants of the row polynomials R(n,x) are found under A193678. - Wolfdieter Lang, Aug 27 2011
The determinant of the N X N matrix M(N) with entries M(N;n,m) = R(m-1,x[n]), 1 <= n,m <= N, N>=1, and any x[n], is identical with twice the Vandermondian Det(V(N)) with matrix entries V(N;n,m) = x[n]^(m-1). This is an instance of the general theorem given in the Vein-Dale reference on p. 59. Note that R(0,x) = 2 (not 1). See also the comments from Aug 26 2013 under A049310 and from Aug 27 2013 under A000178. - Wolfdieter Lang, Aug 27 2013
This triangle a(n,m) is also used to express in the regular (2*(n+1))-gon, inscribed in a circle of radius R, the length ratio side/R, called s(2*(n+1)), as a polynomial in rho(2*(n+1)), the length ratio (smallest diagonal)/side. See the bisections ((-1)^(k-s))*A111125(k,s) and A127677 for comments and examples. - Wolfdieter Lang, Oct 05 2013
From Tom Copeland, Nov 08 2015: (Start)
These are the characteristic polynomials a_n(x) = 2*T_n(x/2) for the adjacency matrix of the Coxeter simple Lie algebra B_n, related to the Cheybshev polynomials of the first kind, T_n(x) = cos(n*q) with x = cos(q) (see p. 20 of Damianou). Given the polynomial (x - t)*(x - 1/t) = 1 - (t + 1/t)*x + x^2 = e2 - e1*x + x^2, the symmetric power sums p_n(t,1/t) = t^n + t^(-n) of the zeros of this polynomial may be expressed in terms of the elementary symmetric polynomials e1 = t + 1/t = y and e2 = t*1/t = 1 as p_n(t,1/t) = a_n(y) = F(n,-y,1,0,0,...), where F(n,b1,b2,...,bn) are the Faber polynomials of A263916.
The partial sum of the first n+1 rows given t and y = t + 1/t is PS(n,t) = Sum_{k=0..n} a_n(y) = (t^(n/2) + t^(-n/2))*(t^((n+1)/2) - t^(-(n+1)/2)) / (t^(1/2) - t^(-1/2)). (For n prime, this is related simply to the cyclotomic polynomials.)
Then a_n(y) = PS(n,t) - PS(n-1,t), and for t = e^(iq), y = 2*cos(q), and, therefore, a_n(2*cos(q)) = PS(n,e^(iq)) - PS(n-1,e^(iq)) = 2*cos(nq) = 2*T_n(cos(q)) with PS(n,e^(iq)) = 2*cos(nq/2)*sin((n+1)q/2) / sin(q/2).
(End)
R(45, x) is the famous polynomial used by Adriaan van Roomen (Adrianus Romanus) in his Ideae mathematicae from 1593 to pose four problems, solved by Viète. See, e.g., the Havil reference, pp. 69-74. - Wolfdieter Lang, Apr 28 2018
From Wolfdieter Lang, May 05 2018: (Start)
Some identities for the row polynomials R(n, x) following from the known ones for Chebyshev T-polynomials (A053120) are:
(1) R(-n, x) = R(n, x).
(2) R(n*m, x) = R(n, R(m, x)) = R(m, R(n, x)).
(3) R(2*k+1, x) = (-1)^k*x*S(2*k, sqrt(4-x^2)), k >= 0, with the S row polynomials of A049310.
(4) R(2*k, x) = R(k, x^2-2), k >= 0.
(End)
For y = z^n + z^(-n) and x = z + z^(-1), Hirzebruch notes that y(z) = R(n,x) for the row polynomial of this entry. - Tom Copeland, Nov 09 2019

Examples

			Row n=4: [2,0,-4,0,1] stands for the polynomial 2*y^0 - 4*y^2 + 1*y^4. With y^m replaced by 2^(m-1)*x^m this becomes T(4,x) = 1 - 8*x^2 + 8*x^4.
Triangle begins:
n\m   0   1   2   3   4   5   6   7   8   9  10 ...
0:    2
1:    0   1
2:   -2   0   1
3:    0  -3   0   1
4:    2   0  -4   0   1
5:    0   5   0  -5   0   1
6:   -2   0   9   0  -6   0   1
7:    0  -7   0  14   0  -7   0   1
8:    2   0 -16   0  20   0  -8   0   1
9:    0   9   0 -30   0  27   0  -9   0   1
10:  -2   0  25   0 -50   0  35   0 -10   0   1 ...
Factorization into minimal C-polynomials:
R(12,x) = R((2^2)*3,x) = C(24,x)*C(8,x) = C((2^3)*1,x)*C((2^3)*3,x). - _Wolfdieter Lang_, Jul 31 2011
		

References

  • Julian Havil, The Irrationals, A Story of the Numbers You Can't Count On, Princeton University Press, Princeton and Oxford, 2012, pp. 69-74.
  • F. Hirzebruch et al., Manifolds and Modular Forms, Vieweg 1994 pp. 77, 105.
  • R. Vein and P. Dale, Determinants and Their Applications in Mathematical Physics, Springer, 1999.

Crossrefs

Row sums (signed): A057079(n-1). Row sums (unsigned): A000032(n) (Lucas numbers). Alternating row sums: A099837(n+3).
Bisection: A127677 (even n triangle, without zero entries), ((-1)^(n-m))*A111125(n, m) (odd n triangle, without zero entries).

Programs

  • Maple
    seq(seq(coeff(2*orthopoly[T](n,x/2),x,j),j=0..n),n=0..20); # Robert Israel, Aug 04 2015
  • Mathematica
    a[n_, k_] := SeriesCoefficient[(2 - t*x)/(1 - t*x + x^2), {x, 0, n}, {t, 0, k}]; Flatten[Table[a[n, k], {n, 0, 12}, {k, 0, n}]] (* L. Edson Jeffery, Nov 02 2017 *)

Formula

a(n,0) = 0 if n is odd, a(n,0) = 2*(-1)^(n/2) if n is even, else a(n,m) = t(n,m)/2^(m-1) with t(n,m):=A053120(n,m) (coefficients of Chebyshev T-polynomials).
G.f. for m-th column (signed triangle): 2/(1+x^2) if m=0 else (x^m)*(1-x^2)/(1+x^2)^(m+1).
Riordan type matrix ((1-x^2)/(1+x^2),x/(1+x^2)) if one puts a(0,0)=1 (instead of 2).
O.g.f. for row polynomials: R(x,z) := Sum_{n>=0} R(n,x)*z^n = (2-x*z)*S(x,z), with the o.g.f. S(x,z) = 1/(1 - x*z + z^2) for the S-polynomials (see A049310).
Note that R(n,x) = R(2*n,sqrt(2+x)), n>=0 (from the o.g.f.s of both sides). - Wolfdieter Lang, Jun 03 2011
a(n,m) := 0 if n < m or n+m odd; a(n,0) = 2*(-1)^(n/2) (n even); else a(n,m) = ((-1)^((n+m)/2 + m))*n*binomial((n+m)/2-1,m-1)/m.
Recursion for n >= 2 and m >= 2: a(n,m) = a(n-1,m-1) - a(n-2,m), a(n,m) = 0 if n < m, a(2*k,1) = 0, a(2*k+1,1) = (2*k+1)*(-1)^k. In addition, for column m=0: a(2*k,0) = 2*(-1)^k, a(2*k+1,0) = 0, k>=0.
Chebyshev T(n,x) = Sum{m=0..n} a(n,m)*2^(m-1)*x^m. - Wolfdieter Lang, Jun 03 2011
R(n,x) = 2*T(n,x/2) = S(n,x) - S(n-2,x), n>=0, with Chebyshev's T- and S-polynomials, showing that they are integer and monic polynomials. - Wolfdieter Lang, Nov 08 2011
From Tom Copeland, Nov 08 2015: (Start)
a(n,x) = sqrt(2 + a(2n,x)), or 2 + a(2n,x) = a(n,x)^2, is a reflection of the relation of the Chebyshev polynomials of the first kind to the cosine and the half-angle formula, cos(q/2)^2 = (1 + cos(q))/2.
Examples: For n = 2, -2 + x^2 = sqrt(2 + 2 - 4*x^2 + x^4).
For n = 3, -3*x + x^3 = sqrt(2 - 2 + 9*x^2 - 6*x^4 + x^6).
(End)
L(x,h1,h2) = -log(1 - h1*x + h2*x^2) = Sum_{n>0} F(n,-h1,h2,0,...,0) x^n/n = h1*x + (-2*h2 + h1^2) x^2/2 + (-3*h1*h2 + h1^3) x^3/3 + ... is a log series generator of the bivariate row polynomials where T(0,0) = 0 and F(n,b1,b2,...,bn) are the Faber polynomials of A263916. exp(L(x,h1,h2)) = 1 / (1 - h1*x + h2*x^2) is the o.g.f. of A049310. - Tom Copeland, Feb 15 2016

Extensions

Name changed and table rewritten by Wolfdieter Lang, Nov 08 2011

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

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.

Crossrefs

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

A111125 Triangle read by rows: T(k,s) = ((2*k+1)/(2*s+1))*binomial(k+s,2*s), 0 <= s <= k.

Original entry on oeis.org

1, 3, 1, 5, 5, 1, 7, 14, 7, 1, 9, 30, 27, 9, 1, 11, 55, 77, 44, 11, 1, 13, 91, 182, 156, 65, 13, 1, 15, 140, 378, 450, 275, 90, 15, 1, 17, 204, 714, 1122, 935, 442, 119, 17, 1, 19, 285, 1254, 2508, 2717, 1729, 665, 152, 19, 1, 21, 385, 2079, 5148, 7007, 5733, 2940, 952, 189, 21, 1
Offset: 0

Views

Author

N. J. A. Sloane, Oct 16 2005

Keywords

Comments

Riordan array ((1+x)/(1-x)^2, x/(1-x)^2). Row sums are A002878. Diagonal sums are A003945. Inverse is A113187. An interesting factorization is (1/(1-x), x/(1-x))(1+2*x, x*(1+x)). - Paul Barry, Oct 17 2005
Central coefficients of rows with odd numbers of term are A052227.
From Wolfdieter Lang, Jun 26 2011: (Start)
T(k,s) appears as T_s(k) in the Knuth reference, p. 285.
This triangle is related to triangle A156308(n,m), appearing in this reference as U_m(n) on p. 285, by T(k,s) - T(k-1,s) = A156308(k,s), k>=s>=1 (identity on p. 286). T(k,s) = A156308(k+1,s+1) - A156308(k,s+1), k>=s>=0 (identity on p. 286).
(End)
A111125 is jointly generated with A208513 as an array of coefficients of polynomials v(n,x): initially, u(1,x)= v(1,x)= 1; for n>1, u(n,x)= u(n-1,x) +x*(x+1)*v(n-1) and v(n,x)= u(n-1,x) +x*v(n-1,x) +1. See the Mathematica section. The columns of A111125 are identical to those of A208508. Here, however, the alternating row sums are periodic (with period 1,2,1,-1,-2,-1). - Clark Kimberling, Feb 28 2012
This triangle T(k,s) (with signs and columns scaled with powers of 5) appears in the expansion of Fibonacci numbers F=A000045 with multiples of odd numbers as indices in terms of odd powers of F-numbers. See the Jennings reference, p. 108, Theorem 1. Quoted as Lemma 3 in the Ozeki reference given in A111418. The formula is: F_{(2*k+1)*n} = Sum_{s=0..k} ( T(k,s)*(-1)^((k+s)*n)*5^s*F_{n}^(2*s+1) ), k >= 0, n >= 0. - Wolfdieter Lang, Aug 24 2012
From Wolfdieter Lang, Oct 18 2012: (Start)
This triangle T(k,s) appears in the formula x^(2*k+1) - x^(-(2*k+1)) = Sum_{s=0..k} ( T(k,s)*(x-x^(-1))^(2*s+1) ), k>=0. Prove the inverse formula (due to the Riordan property this will suffice) with the binomial theorem. Motivated to look into this by the quoted paper of Wang and Zhang, eq. (1.4).
Alternating row sums are A057079.
The Z-sequence of this Riordan array is A217477, and the A-sequence is (-1)^n*A115141(n). For the notion of A- and Z-sequences for Riordan triangles see a W. Lang link under A006232. (End)
The signed triangle ((-1)^(k-s))*T(k,s) gives the coefficients of (x^2)^s of the polynomials C(2*k+1,x)/x, with C the monic integer Chebyshev T-polynomials whose coefficients are given in A127672 (C is there called R). See the odd numbered rows there. This signed triangle is the Riordan array ((1-x)/(1+x)^2, x/(1+x)^2). Proof by comparing the o.g.f. of the row polynomials where x is replaced by x^2 with the odd part of the bisection of the o.g.f. for C(n,x)/x. - Wolfdieter Lang, Oct 23 2012
From Wolfdieter Lang, Oct 04 2013: (Start)
The signed triangle S(k,s) := ((-1)^(k-s))*T(k,s) (see the preceding comment) is used to express in a (4*(k+1))-gon the length ratio s(4*(k+1)) = 2*sin(Pi/4*(k+1)) = 2*cos((2*k+1)*Pi/(4*(k+1))) of a side/radius as a polynomial in rho(4*(k+1)) = 2*cos(Pi/4*(k+1)), the length ratio (smallest diagonal)/side:
s(4*(k+1)) = Sum_{s=0..k} ( S(k,s)*rho(4*(k+1))^(2*s+1) ).
This is to be computed modulo C(4*(k+1), rho(4*(k+1)) = 0, the minimal polynomial (see A187360) in order to obtain s(4*(k+1)) as an integer in the algebraic number field Q(rho(4*(k+1))) of degree delta(4*(k+1)) (see A055034). Thanks go to Seppo Mustonen for asking me to look into the problem of the square of the total length in a regular n-gon, where this formula is used in the even n case. See A127677 for the formula in the (4*k+2)-gon. (End)
From Wolfdieter Lang, Aug 14 2014: (Start)
The row polynomials for the signed triangle (see the Oct 23 2012 comment above), call them todd(k,x) = Sum_{s=0..k} ( (-1)^(k-s)*T(k,s)*x^s ) = S(k, x-2) - S(k-1, x-2), k >= 0, with the Chebyshev S-polynomials (see their coefficient triangle (A049310) and S(-1, x) = 0), satisfy the recurrence todd(k, x) = (-1)^(k-1)*((x-4)/2)*todd(k-1, 4-x) + ((x-2)/2)*todd(k-1, x), k >= 1, todd(0, x) = 1. From the Aug 03 2014 comment on A130777.
This leads to a recurrence for the signed triangle, call it S(k,s) as in the Oct 04 2013 comment: S(k,s) = (1/2)*(1 + (-1)^(k-s))*S(k-1,s-1) + (2*(s+1)*(-1)^(k-s) - 1)*S(k-1,s) + (1/2)*(-1)^(k-s)*Sum_{j=0..k-s-2} ( binomial(j+s+2,s)*4^(j+2)* S(k-1, s+1+j) ) for k >= s >= 1, and S(k,s) = 0 if k < s and S(k,0) = (-1)^k*(2*k+1). Note that the recurrence derived from the Riordan A-sequence A115141 is similar but has simpler coefficients: S(k,s) = sum(A115141(j)*S(k-1,s-1+j), j=0..k-s), k >= s >=1.
(End)
From Tom Copeland, Nov 07 2015: (Start)
Rephrasing notes here: Append an initial column of zeros, except for a 1 at the top, to A111125 here. Then the partial sums of the columns of this modified entry are contained in A208513. Append an initial row of zeros to A208513. Then the difference of consecutive pairs of rows of the modified A208513 generates the modified A111125. Cf. A034807 and A127677.
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. 20 and 21) and Damianou and Evripidou (p. 7).
As suggested by the equations on p. 7 of Damianou and Evripidou, the signed row polynomials of this entry are given by (p(n,x))^2 = (A(2*n+1, x) + 2)/x = (F(2*n+1, (2-x), 1, 0, 0, ... ) + 2)/x = F(2*n+1, -x, 2*x, -3*x, ..., (-1)^n n*x)/x = -F(2*n+1, x, 2*x, 3*x, ..., n*x)/x, where A(n,x) are the polynomials of A127677 and F(n, ...) are the Faber polynomials of A263196. Cf. A127672 and A127677.
(End)
The row polynomials P(k, x) of the signed triangle S(k, s) = ((-1)^(k-s))*T(k, s) are given from the row polynomials R(2*k+1, x) of triangle A127672 by
P(k, x) = R(2*k+1, sqrt(x))/sqrt(x). - Wolfdieter Lang, May 02 2021

Examples

			Triangle T(k,s) begins:
k\s  0    1     2     3     4     5     6    7    8   9 10
0:   1
1:   3    1
2:   5    5     1
3:   7   14     7     1
4:   9   30    27     9     1
5:  11   55    77    44    11     1
6:  13   91   182   156    65    13     1
7:  15  140   378   450   275    90    15    1
8:  17  204   714  1122   935   442   119   17    1
9:  19  285  1254  2508  2717  1729   665  152   19   1
10: 21  385  2079  5148  7007  5733  2940  952  189  21  1
... Extended and reformatted by _Wolfdieter Lang_, Oct 18 2012
Application for Fibonacci numbers F_{(2*k+1)*n}, row k=3:
F_{7*n} = 7*(-1)^(3*n)*F_n + 14*(-1)^(4*n)*5*F_n^3 + 7*(-1)^(5*n)*5^2*F_n^5 + 1*(-1)^(6*n)*5^3*F_n^7, n>=0. - _Wolfdieter Lang_, Aug 24 2012
Example for the  Z- and A-sequence recurrences  of this Riordan triangle: Z = A217477 = [3,-4,12,-40,...]; T(4,0) = 3*7 -4*14 +12*7 -40*1 = 9. A =  [1, 2, -1, 2, -5, 14, ..]; T(5,2) = 1*30 + 2*27 - 1*9 + 2*1= 77. _Wolfdieter Lang_, Oct 18 2012
Example for the (4*(k+1))-gon length ratio s(4*(k+1))(side/radius) as polynomial in the ratio rho(4*(k+1)) ((smallest diagonal)/side): k=0, s(4) = 1*rho(4) = sqrt(2); k=1, s(8) = -3*rho(8) + rho(8)^3 = sqrt(2-sqrt(2)); k=2, s(12) = 5*rho(12) - 5*rho(12)^3 + rho(12)^5, and C(12,x) = x^4 - 4*x^2 + 1, hence rho(12)^5 = 4*rho(12)^3 - rho(12), and s(12) = 4*rho(12) - rho(12)^3 = sqrt(2 - sqrt(3)). - _Wolfdieter Lang_, Oct 04 2013
Example for the recurrence for the signed triangle S(k,s)= ((-1)^(k-s))*T(k,s) (see the Aug 14 2014 comment above):
S(4,1) = 0 + (-2*2 - 1)*S(3,1) - (1/2)*(3*4^2*S(3,2) + 4*4^3*S(3,3)) = - 5*14 - 3*8*(-7) - 128*1 = -30. The recurrence from the Riordan A-sequence A115141 is S(4,1) = -7 -2*14 -(-7) -2*1 = -30. - _Wolfdieter Lang_, Aug 14 2014
		

Crossrefs

Mirror image of A082985, which see for further references, etc.
Also closely related to triangles in A098599 and A100218.

Programs

  • Magma
    [((2*n+1)/(n+k+1))*Binomial(n+k+1, 2*k+1): k in [0..n], n in [0..12]];  // G. C. Greubel, Feb 01 2022
  • Mathematica
    (* First program *)
    u[1, x_]:=1; v[1, x_]:=1; z=16;
    u[n_, x_]:= u[n-1, x] + x*v[n-1, x];
    v[n_, x_]:= u[n-1, x] + (x+1)*v[n-1, x] + 1;
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]  (* A208513 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]  (* A111125 *) (* Clark Kimberling, Feb 28 2012 *)
    (* Second program *)
    T[n_, k_]:= ((2*n+1)/(2*k+1))*Binomial[n+k, 2*k];
    Table[T[n, k], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Feb 01 2022 *)
  • Sage
    @CachedFunction
    def T(n,k):
        if n< 0: return 0
        if n==0: return 1 if k == 0 else 0
        h = 3*T(n-1,k) if n==1 else 2*T(n-1,k)
        return T(n-1,k-1) - T(n-2,k) - h
    A111125 = lambda n,k: (-1)^(n-k)*T(n,k)
    for n in (0..9): [A111125(n,k) for k in (0..n)] # Peter Luschny, Nov 20 2012
    

Formula

T(k,s) = ((2*k+1)/(2*s+1))*binomial(k+s,2*s), 0 <= s <= k.
From Peter Bala, Apr 30 2012: (Start)
T(n,k) = binomial(n+k,2*k) + 2*binomial(n+k,2*k+1).
The row generating polynomials P(n,x) are a generalization of the Morgan-Voyce polynomials b(n,x) and B(n,x). They satisfy the recurrence equation P(n,x) = (x+2)*P(n-1,x) - P(n-2,x) for n >= 2, with initial conditions P(0,x) = 1, P(1,x) = x+r+1 and with r = 2. The cases r = 0 and r = 1 give the Morgan-Voyce polynomials A085478 and A078812 respectively. Andre-Jeannin has considered the case of general r.
P(n,x) = U(n+1,1+x/2) + U(n,1+x/2), where U(n,x) denotes the Chebyshev polynomial of the second kind - see A053117. P(n,x) = (2/x)*(T(2*n+2,u)-T(2*n,u)), where u = sqrt((x+4)/4) and T(n,x) denotes the Chebyshev polynomial of the first kind - see A053120. P(n,x) = Product_{k = 1..n} ( x + 4*(sin(k*Pi/(2*n+1)))^2 ). P(n,x) = 1/x*(b(n+1,x) - b(n-1,x)) and P(n,x) = 1/x*{(b(2*n+2,x)+1)/b(n+1,x) - (b(2*n,x)+1)/b(n,x)}, where b(n,x) := Sum_{k = 0..n} binomial(n+k,2*k)*x^k are the Morgan-Voyce polynomials of A085478. Cf. A211957.
(End)
From Wolfdieter Lang, Oct 18 2012 (Start)
O.g.f. column No. s: ((1+x)/(1-x)^2)*(x/(1-x)^2)^s, s >= 0. (from the Riordan data given in a comment above).
O.g.f. of the row polynomials R(k,x):= Sum_{s=0..k} ( T(k,s)*x^s ), k>=0: (1+z)/(1-(2+x)*z+z^2) (from the Riordan property).
(End)
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k), T(0,0) = 1, T(1,0) = 3, T(1,1) = 1, T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Nov 12 2013

Extensions

More terms from Paul Barry, Oct 17 2005

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)

A263916 Coefficients of the Faber partition polynomials.

Original entry on oeis.org

-1, -2, 1, -3, 3, -1, -4, 4, 2, -4, 1, -5, 5, 5, -5, -5, 5, -1, -6, 6, 6, -6, 3, -12, 6, -2, 9, -6, 1, -7, 7, 7, -7, 7, -14, 7, -7, -7, 21, -7, 7, -14, 7, -1, -8, 8, 8, -8, 8, -16, 8, 4, -16, -8, 24, -8, -8, 12, 24, -32, 8, 2, -16, 20, -8, 1
Offset: 1

Author

Tom Copeland, Oct 29 2015

Keywords

Comments

The coefficients of the Faber polynomials F(n,b(1),b(2),...,b(n)) (Bouali, p. 52) in the order of the partitions of Abramowitz and Stegun. Compare with A115131 and A210258.
These polynomials occur in discussions of the Virasoro algebra, univalent function spaces and the Schwarzian derivative, symmetric functions, and free probability theory. They are intimately related to symmetric functions, free probability, and Appell sequences through the raising operator R = x - d log(H(D))/dD for the Appell sequence inverse pair associated to the e.g.f.s H(t)e^(xt) (cf. A094587) and (1/H(t))e^(xt) with H(0)=1.
Instances of the Faber polynomials occur in discussions of modular invariants and modular functions in the papers by Asai, Kaneko, and Ninomiya, by Ono and Rolen, and by Zagier. - Tom Copeland, Aug 13 2019
The Faber polynomials, denoted by s_n(a(t)) where a(t) is a formal power series defined by a product formula, are implicitly defined by equation 13.4 on p. 62 of Hazewinkel so as to extract the power sums of the reciprocals of the zeros of a(t). This is the Newton identity expressing the power sum symmetric polynomials in terms of the elementary symmetric polynomials/functions. - Tom Copeland, Jun 06 2020
From Tom Copeland, Oct 15 2020: (Start)
With a_n = n! * b_n = (n-1)! * c_n for n > 0, represent a function with f(0) = a_0 = b_0 = 1 as an
A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!
B) ordinary generating function (o.g.f.), or formal power series: f(x) = 1/(1-b.x) = 1 + Sum_{n > 0} b_n * x^n
C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0} c_n * x^n /n.
Expansions of log(f(x)) are given in
I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] = e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials
II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] = log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.
Expansions of exp(f(x)-1) are given in
III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind
IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials
V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] = e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.
Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)

Examples

			F(1,b1) = - b1
F(2,b1,b2) = -2 b2 + b1^2
F(3,b1,b2,b3) = -3 b3 + 3 b1 b2 - b1^3
F(4,b1,...) = -4 b4 + 4 b1 b3 + 2 b2^2  - 4 b1^2 b2 + b1^4
F(5,...) = -5 b5 + 5 b1 b4 + 5 b2 b3 - 5 b1^2 b3 - 5 b1 b2^2 + 5 b1^3 b2 - b1^5
------------------------------
IF(1,b1) = -b1
IF(2,b1,,b2) = -b2 + b1^2
IF(3,b1,b2,b3) = -2 b3 + 3 b1 b2 - b1^3
IF(4,b1,...) = -6 b4 + 8 b1 b3 + 3 b2^2  - 6 b1^2 b2 + b1^4
IF(5,...) = -24 b5 + 30 b1 b4 + 20 b2 b3 - 20 b1^2 b3 - 15 b1 b2^2 + 10 b1^3 b2 - b1^5
------------------------------
For 1/(1+x)^2 = 1- 2x + 3x^2 - 4x^3 + 5x^4 - ..., F(n,-2,3,-4,...) = (-1)^(n+1) 2.
------------------------------
F(n,x,2x,...,nx), F(n,-x,2x,-3x,...,(-1)^n n*x), and F(n,(2-x),1,0,0,...) are related to the Chebyshev polynomials through A127677 and A111125. See also A110162, A156308, A208513, A217476, and A220668.
------------------------------
For b1 = p, b2 = q, and all other indeterminates 0, see A113279 and A034807.
For b1 = -y, b2 = 1 and all other indeterminates 0, see A127672.
		

References

  • H. Airault, "Symmetric sums associated to the factorization of Grunsky coefficients," in Groups and Symmetries: From Neolithic Scots to John McKay, CRM Proceedings and Lecture Notes: Vol. 47, edited by J. Harnad and P. Winternitz, American Mathematical Society, 2009.
  • D. Bleeker and B. Booss, Index Theory with Applications to Mathematics and Physics, International Press, 2013, (see section 16.7 Characteristic Classes and Curvature).
  • M. Hazewinkel, Formal Groups and Applications, Academic Press, New York San Francisco London, 1978, p. 120.
  • F. Hirzebruch, Topological methods in algebraic geometry. Second, corrected printing of the third edition. Die Grundlehren der Mathematischen Wissenschaften, Band 131 Springer-Verlag, Berlin Heidelberg New York, 1978, p. 11 and 92.
  • D. Knutson, λ-Rings and the Representation Theory of the Symmetric Group, Lect. Notes in Math. 308, Springer-Verlag, 1973, p. 35.
  • D. Yau, Lambda-Rings, World Scientific Publishing Co., Singapore, 2010, p. 45.

Crossrefs

Programs

  • Mathematica
    F[0] = 1; F[1] = -b[1]; F[2] = b[1]^2 - 2 b[2]; F[n_] := F[n] = -b[1] F[n - 1] - Sum[b[n - k] F[k], {k, 1, n - 2}] - n b[n] // Expand;
    row[n_] := (List @@ F[n]) /. b[_] -> 1 // Reverse;
    Table[row[n], {n, 1, 8}] // Flatten // Rest (* Jean-François Alcover, Jun 12 2017 *)

Formula

-log(1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n.
-d(1 + b(1) x + b(2) x^2 + ...)/dx / (1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) x^(n-1).
F(n,b(1),...,b(n)) = -n*b(n) - Sum_{k=1..n-1} b(n-k)*F(k,b(1),...,b(k)).
Umbrally, with B(x) = 1 + b(1) x + b(2) x^2 + ..., B(x) = exp[log(1-F.x)] and 1/B(x) = exp[-log(1-F.x)], establishing a connection to the e.g.f. of A036039 and the symmetric polynomials.
The Stirling partition polynomials of the first kind St1(n,b1,b2,...,bn;-1) = IF(n,b1,b2,...,bn) (cf. the Copeland link Lagrange a la Lah, signed A036039, and p. 184 of Airault and Bouali), i.e., the cyclic partition polynomials for the symmetric groups, and the Faber polynomials form an inverse pair for isolating the indeterminates in their definition, for example, F(3,IF(1,b1),IF(2,b1,b2)/2!,IF(3,b1,b2,b3)/3!)= b3, with bk = b(k), and IF(3,F(1,b1),F(2,b1,b2),F(3,b1,b2,b3))/3!= b3.
The polynomials specialize to F(n,t,t,...) = (1-t)^n - 1.
See Newton Identities on Wikipedia on relation between the power sum symmetric polynomials and the complete homogeneous and elementary symmetric polynomials for an expression in multinomials for the coefficients of the Faber polynomials.
(n-1)! F(n,x[1],x[2]/2!,...,x[n]/n!) = - p_n(x[1],...,x[n]), where p_n are the cumulants of A127671 expressed in terms of the moments x[n]. - Tom Copeland, Nov 17 2015
-(n-1)! F(n,B(1,x[1]),B(2,x[1],x[2])/2!,...,B(n,x[1],...,x[n])/n!) = x[n] provides an extraction of the indeterminates of the complete Bell partition polynomials B(n,x[1],...,x[n]) of A036040. Conversely, IF(n,-x[1],-x[2],-x[3]/2!,...,-x[n]/(n-1)!) = B(n,x[1],...,x[n]). - Tom Copeland, Nov 29 2015
For a square matrix M, determinant(I - x M) = exp[-Sum_{k>0} (trace(M^k) x^k / k)] = Sum_{n>0} [ P_n(-trace(M),-trace(M^2),...,-trace(M^n)) x^n/n! ] = 1 + Sum_{n>0} (d[n] x^n), where P_n(x[1],...,x[n]) are the cycle index partition polynomials of A036039 and d[n] = P_n(-trace(M),-trace(M^2),...,-trace(M^n)) / n!. Umbrally, det(I - x M)= exp[log(1 - b. x)] = exp[P.(-b_1,..,-b_n)x] = 1 / (1-d.x), where b_k = tr(M^k). Then F(n,d[1],...,d[n]) = tr[M^n]. - Tom Copeland, Dec 04 2015
Given f(x) = -log(g(x)) = -log(1 + b(1) x + b(2) x^2 + ...) = Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n, action on u_n = F(n,b(1),...,b(n)) with A133932 gives the compositional inverse finv(x) of f(x), with F(1,b(1)) not equal to zero, and f(g(finv(x))) = f(e^(-x)). Note also that exp(f(x)) = 1 / g(x) = exp[Sum_{n>=1} F(n,b(1),...,b(n)) * x^n/n] implies relations among A036040, A133314, A036039, and the Faber polynomials. - Tom Copeland, Dec 16 2015
The Dress and Siebeneicher paper gives combinatorial interpretations and various relations that the Faber polynomials must satisfy for integral values of its arguments. E.g., Eqn. (1.2) p. 2 implies [2 * F(1,-1) + F(2,-1,b2) + F(4,-1,b2,b3,b4)] mod(4) = 0. This equation implies that [F(n,b1,b2,...,bn)-(-b1)^n] mod(n) = 0 for n prime. - Tom Copeland, Feb 01 2016
With the elementary Schur polynomials S(n,a_1,a_2,...,a_n) = Lah(n,a_1,a_2,...,a_n) / n!, where Lah(n,...) are the refined Lah polynomials of A130561, F(n,S(1,a_1),S(2,a_1,a_2),...,S(n,a_1,...,a_n)) = -n * a_n since sum_{n > 0} a_n x^n = log[sum{n >= 0} S(n,a_1,...,a_n) x^n]. Conversely, S(n,-F(1,a_1),-F(2,a_1,a_2)/2,...,-F(n,a_1,...,a_n)/n) = a_n. - Tom Copeland, Sep 07 2016
See Corollary 3.1.3 on p. 38 of Ardila and Copeland's two MathOverflow links to relate the Faber polynomials, with arguments being the signed elementary symmetric polynomials, to the logarithm of determinants, traces of powers of an adjacency matrix, and number of walks on graphs. - Tom Copeland, Jan 02 2017
The umbral inverse polynomials IF appear on p. 19 of Konopelchenko as partial differential operators. - Tom Copeland, Nov 19 2018

Extensions

More terms from Jean-François Alcover, Jun 12 2017

A110162 Riordan array ((1-x)/(1+x), x/(1+x)^2).

Original entry on oeis.org

1, -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

Paul Barry, Jul 14 2005

Keywords

Comments

Inverse of Riordan array A094527. Rows sums are A099837. Diagonal sums are A110164. Product of Riordan array A102587 and inverse binomial transform (1/(1+x), x/(1+x)).
Coefficients of polynomials related to Cartan matrices of types C_n and B_n: p(x, n) = (-2 + x)*p(x, n - 1) - p(x, n - 2), with p(x,0) = 1; p(x,1) = 2-x; p(x,2) = x^2-4*x-2. - Roger L. Bagula, Apr 12 2008
From Wolfdieter Lang, Nov 16 2012: (Start)
The alternating row sums are given in A219233.
For n >= 1 the row polynomials in the variable x^2 are R(2*n,x):=2*T(2*n,x/2) with Chebyshev's T-polynomials. See A127672 and also the triangle A127677.
(End)
From Peter Bala, Jun 29 2015: (Start)
Riordan array has the form ( x*h'(x)/h(x), h(x) ) with h(x) = x/(1 + x)^2 and so belongs to the hitting time subgroup H of the Riordan group (see Peart and Woan).
T(n,k) = [x^(n-k)] f(x)^n with f(x) = (1 - 2*x + sqrt(1 - 4*x))/2. In general the (n,k)th entry of the hitting time array ( x*h'(x)/h(x), h(x) ) has the form [x^(n-k)] f(x)^n, where f(x) = x/( series reversion of h(x) ). (End)

Examples

			Triangle T(n,k) begins:
m\k  0    1    2     3     4     5     6    7    8   9 10 ...
0:   1
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 16 2012
Row polynomial n=2: P(2,x) = 2 - 4*x + x^2. R(4,x):= 2*T(4,x/2) = 2 - 4*x^2 + x^4. For P and R see a comment above. - _Wolfdieter Lang_, Nov 16 2012.
		

Crossrefs

Cf. A128411. See A127677 for an almost identical triangle.

Programs

  • Magma
    /* As triangle */ [[(-1)^(n-k)*(Binomial(n+k,n-k) + Binomial(n+k-1,n-k-1)): k in [0..n]]: n in [0.. 12]]; // Vincenzo Librandi, Jun 30 2015
    
  • Mathematica
    Table[If[n==0 && k==0, 1, (-1)^(n-k)*(Binomial[n+k, n-k] + Binomial[n+k-1, n-k-1])], {n, 0, 15}, {k, 0, n}]//Flatten (* G. C. Greubel, Dec 16 2018 *)
  • PARI
    {T(n,k) = (-1)^(n-k)*(binomial(n+k,n-k) + binomial(n+k-1,n-k-1))};
    for(n=0, 12, for(k=0, n, print1(T(n,k), ", "))) \\ G. C. Greubel, Dec 16 2018
    
  • Sage
    [[(-1)^(n-k)*(binomial(n+k,n-k) + binomial(n+k-1,n-k-1)) for k in range(n+1)] for n in range(12)] # G. C. Greubel, Dec 16 2018

Formula

T(n,k) = (-1)^(n-k)*(C(n+k,n-k) + C(n+k-1,n-k-1)), with T(0,0) = 1. - Paul Barry, Mar 22 2007
From Wolfdieter Lang, Nov 16 2012: (Start)
O.g.f. row polynomials P(n,x) := Sum(T(n,k)*x^k, k=0..n): (1-z^2)/(1+(x-2)*z+z^2) (from the Riordan property).
O.g.f. column No. k: ((1-x)/(1+x))*(x/(1+x)^2)^k, k >= 0.
T(0,0) = 1, T(n,k) = (-1)^(n-k)*(2*n/(n+k))*binomial(n+k,n-k), n>=1, and T(n,k) = 0 if n < k. (From the Chebyshev T-polynomial formula due to Waring's formula.)
(End)
T(n,k) = -2*T(n-1,k) + T(n-1,k-1) - T(n-2,k), T(0,0)=1, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Nov 29 2013

A100051 A Chebyshev transform of 1,1,1,...

Original entry on oeis.org

1, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1
Offset: 0

Author

Paul Barry, Oct 31 2004

Keywords

Comments

1, followed by period 6: repeat [1, -1, -2, -1, 1, 2]. - Joerg Arndt, Aug 28 2024
A Chebyshev transform of 1/(1-x): if A(x) is the g.f. of a sequence, map it to ((1-x^2)/(1+x^2))A(x/(1+x^2)).
Transform of 1/(1+x) under the mapping g(x)->((1+x)/(1-x))g(x/(1-x)^2). - Paul Barry, Dec 01 2004
Multiplicative with a(p^e) = -1 if p = 2; -2 if p = 3; 1 otherwise. - David W. Wilson, Jun 10 2005

Examples

			G.f. = 1 + x - x^2 - 2*x^3 - x^4 + x^5 + 2*x^6 + x^7 - x^8 - 2*x^9 - x^10 + ...
		

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[(1 - x^2)/(1 - x + x^2), {x,0,50}], x] (* G. C. Greubel, May 03 2017 *)
    LinearRecurrence[{1,-1},{1,1,-1},80] (* Harvey P. Dale, Mar 25 2019 *)
  • PARI
    {a(n) = - (n == 0) + [2, 1, -1, -2, -1, 1][n%6 + 1]}; /* Michael Somos, Mar 21 2011 */

Formula

From Paul Barry, Dec 01 2004: (Start)
G.f.: (1-x^2)/(1-x+x^2).
a(n) = a(n-1) - a(n-2), n>2.
a(n) = n*Sum_{k=0..floor(n/2)} (-1)^k*binomial(n-k, k)/(n-k).
a(n) = Sum_{k=0..n} binomial(n+k, 2k)*(2n/(n+k))*(-1)^k, n>1. (End)
Moebius transform is length 6 sequence [1, -2, -3, 0, 0, 6].
Euler transform of length 6 sequence [1, -2, -1, 0, 0, 1].
a(n) = a(-n). a(n) = c_6(n) if n>1 where c_k(n) is Ramanujan's sum. - Michael Somos, Mar 21 2011
a(n) = A087204(n), n>0. - R. J. Mathar, Sep 02 2008
a(n) = A057079(n+1), n>0. Dirichlet g.f. zeta(s) *(1-2^(1-s)-3^(1-s)+6^(1-s)). - R. J. Mathar, Apr 11 2011

A132460 Irregular triangle read by rows of the initial floor(n/2) + 1 coefficients of 1/C(x)^n, where C(x) is the g.f. of the Catalan sequence (A000108).

Original entry on oeis.org

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

Author

Paul D. Hanna, Aug 21 2007

Keywords

Comments

The length of row n is A008619(n).
Essentially equals a signed version of A034807, the triangle of Lucas polynomials. The initial n coefficients of 1/C(x)^n consist of row n followed by floor((n-1)/2) zeros for n > 0.
For the following formula for 1/C(x)^n see the W. Lang reference, proposition 1 on p. 411:
1/C(x)^n = (sqrt(x))^n*(S(n,1/sqrt(x)) - sqrt(x)*S(n-1,1/sqrt(x))*C(x)), n >= 0, with the Chebyshev polynomials S(n,x) with coefficients given in A049310. See also the coefficient array A115139 for P(n,x) = (sqrt(x)^(n-1))*S(n-1, 1/sqrt(x)). - Wolfdieter Lang, Sep 14 2013
This triangular array is composed of interleaved rows of reversed, A127677 (cf. A156308, A217476, A263916) and reversed, signed A111125. - Tom Copeland, Nov 07 2015
It seems that the n-th row lists the coefficients of the HOMFLYPT (HOMFLY) polynomial reduced to one variable for link family n, see Jablan's slide 38. - Andrey Zabolotskiy, Jan 16 2018
For n >= 1 row n gives the coefficients of the Girard-Waring formula for the sum of x1^n + x2^n in terms of the elementary symmetric functions e_1(x1,x2) = x1 + x2 and e_2(x1,x2) = x1*x2. This is an array using the partitions of n, in the reverse Abramowitz-Stegun order, with all partitions with parts larger than 2 eliminated. E.g., n = 4: x1^4 + x2^4 = 1*e1^4 - 4*e1^3*e2 + 2*e1*e2^2. See also A115131, row n = 4, with the mentioned partitions omitted. - Wolfdieter Lang, May 03 2019
Row n lists the coefficients of the n-th Faber polynomial for the replicable function given in A154272 with offset -1. - Ben Toomey, May 12 2020

Examples

			The irregular triangle T(n,k) begins:
n\k 0    1    2    3    4    5    6   7 ...
-------------------------------------------------
0:  1
1:  1
2:  1   -2
3:  1   -3
4:  1   -4    2
5:  1   -5    5
6:  1   -6    9   -2
7:  1   -7   14   -7
8:  1   -8   20  -16    2
9:  1   -9   27  -30    9
10: 1  -10   35  -50   25   -2
11: 1  -11   44  -77   55  -11
12: 1  -12   54 -112  105  -36    2
13: 1  -13   65 -156  182  -91   13
14: 1  -14   77 -210  294 -196   49  -2
... (reformatted - _Wolfdieter Lang_, May 03 2019)
		

Crossrefs

Cf. A000108, A008619, A034807 (Lucas polynomials), A111125, A115131 (Waring numbers), A127677, A132461 (row squared sums), A156308, A217476, A263916.

Programs

  • Mathematica
    T[0, 0] = 1; T[n_, k_] := (-1)^k (Binomial[n-k, k] + Binomial[n-k-1, k-1]);
    Table[T[n, k], {n, 0, 14}, {k, 0, n/2}] // Flatten (* Jean-François Alcover, Jun 04 2018 *)
  • PARI
    {T(n,k)=if(k>n\2,0,(-1)^k*(binomial(n-k, k)+binomial(n-k-1, k-1)))}

Formula

T(n,k) = (-1)^k*( C(n-k,k) + C(n-k-1,k-1) ) for n >= 0, 0 <= k <= floor(n/2).
T(0,0) = 1; T(n,k) = (-1)^k*n*binomial(n-k,k)/(n-k), k = 0..floor(n/2). - Wolfdieter Lang, May 03 2019

A156308 Inverse of triangle S(n,m) defined by sequence A156290, n >= 1, 1 <= m <= n.

Original entry on oeis.org

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

Author

Hartmut F. W. Hoft, Feb 07 2009

Keywords

Comments

From Wolfdieter Lang, Jun 26 2011: (Start)
This triangle S(n,m) appears as U_m(n) in the Knuth reference on p. 285. It is related to the Riordan triangle T_m(n) = A111125(n,m) by S(n,m) = A111125(n,m) - A111125(n-1,m), n >= m >= 1 (identity on p. 286).
Also, S(n,m)-S(n-1,m) = A111125(n-1,m-1), n >= 2, m >= 1 (identity on p. 286). (End)
These polynomials may be expressed in terms of the Faber polynomials of A263916 and are embedded in A127677 and A208513. - Tom Copeland, Nov 06 2015

Examples

			Triangle starts:
  n=1:  1;
  n=2:  4,  1;
  n=3:  9,  6,  1;
  n=4: 16, 20,  8,  1;
  ...
		

Crossrefs

Same as triangle A208513 with the first column truncated.
Columns: A000290 (m=1), A002415 (m=2), A040977 (m=3), A053347 (m=4), A054334 (m=5).

Programs

  • Magma
    [(n/k)*Binomial(n+k-1, 2*k-1): k in [1..n], n in [1..12]]; // G. C. Greubel, Feb 01 2022
  • Mathematica
    S[m_] := Flatten[Table[k/j Binomial[k + j - 1, 2 j - 1], {k, 1, m}, {j, 1, k}]]
  • Sage
    flatten([[(n/k)*binomial(n+k-1, 2*k-1) for k in (1..n)] for n in (1..12)]) # G. C. Greubel, Feb 01 2022
    

Formula

S(n, m) = (n/m) * binomial(n + m - 1, 2*m - 1).
From Peter Bala, May 01 2014: (Start)
The n-th row o.g.f. is polynomial R(n,x) = 2/x*( T(n,(x + 2)/2) - 1 ), where T(n,x) is Chebyshev polynomial of the first kind. They form a divisibility sequence: if n divides m then R(n,x) divides R(m,x) in the ring Z[x].
R(2*n,x) = (x + 4)*U(n-1,(x + 2)/2)^2;
R(2*n + 1,x) = ( U(n,(x + 2)/2) + U(n-1,(x + 2)/2) )^2.
O.g.f.: Sum_{n >= 0} R(n,x)*z^n = z*(1 + z)/( (1 - z)*(1 - (x + 2)*z + z^2) ). (End)
The polynomial R(n,x) defined above satisfies (x + 1/x - 2) * R(n, x + 1/x - 2) = x^n + 1/x^n - 2. - Alexander Burstein, May 23 2021

Extensions

Edited by Max Alekseyev, Mar 05 2018

A208513 Triangle of coefficients of polynomials u(n,x) jointly generated with A111125; see the Formula section.

Original entry on oeis.org

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

Author

Clark Kimberling, Feb 28 2012

Keywords

Comments

The columns of A208513 are identical to those of A208509. Here, however, the alternating row sums are periodic (with period 1,0,-2,-3,-2,0).
From Tom Copeland, Nov 07 2015: (Start)
These polynomials may be expressed in terms of the Faber polynomials of A263916, similar to A127677.
Rephrasing notes in A111125: Append an initial column of zeros except for a 1 at the top to A111125. Then the rows of this entry contain the partial sums of the column sequences of modified A111125; therefore, the difference of consecutive pairs of rows of this entry, modified by appending an initial row of zeros to it, generates the modified A111125. (End)

Examples

			First five rows:
  1;
  1,  1;
  1,  4,  1;
  1,  9,  6, 1;
  1, 16, 20, 8, 1;
First five polynomials u(n,x):
  u(1,x) = 1;
  u(2,x) = 1 +    x;
  u(3,x) = 1 +  4*x +    x^2;
  u(4,x) = 1 +  9*x +  6*x^2 +   x^3;
  u(5,x) = 1 + 16*x + 20*x^2 + 8*x^3 + x^4;
		

Programs

  • Magma
    A208513:= func< n,k | k eq 1 select 1 else (2*(n-1)/(n+k-2))*Binomial(n+k-2, 2*k-2) >;
    [A208513(n,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Feb 02 2022
    
  • Mathematica
    (* First program *)
    u[1, x_]:=1; v[1, x_]:=1; z=16;
    u[n_, x_]:= u[n-1, x] + x*v[n-1, x];
    v[n_, x_]:= u[n-1, x] + (x+1)*v[n-1, x] + 1;
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n,z}];
    TableForm[cu]
    Flatten[%]  (* A208513 *)
    Table[Expand[v[n, x]], {n,z}]
    cv = Table[CoefficientList[v[n, x], x], {n,z}];
    TableForm[cv]
    Flatten[%]  (* A111125 *)
    (* Second program *)
    T[n_, k_]:= If[k==1, 1, ((n-1)/(k-1))*Binomial[n+k-3, 2*k-3]];
    Table[T[n, k], {n,15}, {k,n}]//Flatten (* G. C. Greubel, Feb 02 2022 *)
  • Sage
    def A208513(n,k): return 1 if (k==1) else ((n-1)/(k-1))*binomial(n+k-3, 2*k-3)
    flatten([[A208513(n,k) for k in (1..n)] for n in (1..12)]) # G. C. Greubel, Feb 02 2022

Formula

Coefficients of u(n, x) from the mixed recurrence relations:
u(n,x) = u(n-1,x) + x*v(n-1,x),
v(n,x) = u(n-1,x) + (x+1)*v(n-1,x) + 1,
where u(1,x) = 1, u(2,x) = 1+x, v(1,x) = 1, v(2,x) = 3+x.
From Peter Bala, May 01 2012: (Start)
Working with an offset of 0: T(n,0) = 1; T(n,k) = (n/k)*binomial(n+k-1,2*k-1) = (n/k)*A078812(n,k) for k > 0. Cf. A156308.
O.g.f.: ((1-t)^2 + t^2*x)/((1-t)*((1-t)^2-t*x)) = 1 + (1+x)*t + (1+4*x+x^2)*t^2 + ....
u(n+1,x) = -1 + (b(2*n,x) + 1)/b(n,x), where b(n,x) = Sum_{k = 0..n} binomial(n+k, 2*k)*x^k are the Morgan-Voyce polynomials of A085478.
This triangle is formed from the even numbered rows of A211956 with a factor of 2^(k-1) removed from the k-th column entries.
(End)
T(n, k) = (2*(n-1)/(n+k-2))*binomial(n+k-2, 2*k-2). - G. C. Greubel, Feb 02 2022
Showing 1-10 of 17 results. Next