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-2 of 2 results.

A053120 Triangle of coefficients of Chebyshev's T(n,x) polynomials (powers of x in increasing order).

Original entry on oeis.org

1, 0, 1, -1, 0, 2, 0, -3, 0, 4, 1, 0, -8, 0, 8, 0, 5, 0, -20, 0, 16, -1, 0, 18, 0, -48, 0, 32, 0, -7, 0, 56, 0, -112, 0, 64, 1, 0, -32, 0, 160, 0, -256, 0, 128, 0, 9, 0, -120, 0, 432, 0, -576, 0, 256, -1, 0, 50, 0, -400, 0, 1120, 0, -1280, 0, 512, 0, -11, 0, 220, 0, -1232, 0, 2816, 0, -2816, 0, 1024
Offset: 0

Views

Author

Keywords

Comments

Row sums (signed triangle): A000012 (powers of 1). Row sums (unsigned triangle): A001333(n).
From Wolfdieter Lang, Oct 21 2013: (Start)
The row polynomials T(n,x) equal (S(n,2*x) - S(n-2,2*x))/2, n >= 0, with the row polynomials S from A049310, with S(-1,x) = 0, and S(-2,x) = -1.
The zeros of T(n,x) are x(n,k) = cos((2*k+1)*Pi/(2*n)), k = 0, 1, ..., n-1, n >= 1. (End)
From Wolfdieter Lang, Jan 03 2020 and Paul Weisenhorn: (Start)
The (sub)diagonal sequences {D_{2*k}(m)}{m >= 0}, for k >= 0, have o.g.f. GD{2*k}(x) = (-1)^k*(1-x)/(1-2*x)^(k+1), for k >= 0, and GD_{2*k+1}(x) = 0, for k >= 0. This follows from their o.g.f. GGD(z, x) := Sum_{k>=0} GD_k(x)*z^n which is obtained from the o.g.f. of the T-triangle GT(z, x) = (1-x*z)/(1 - 2*x + z^2) (see the formula section) by GGD(z, x) = GT(z, x/z).
The explicit form is then D_{2*k}(m) = (-1)^k, for m = 0, and
(-1)^k*(2*k+m)*2^(m-1)*risefac(k+1, m-1)/m!, for m >= 1, with the rising factorial risefac(x, n). (End)

Examples

			The triangle a(n,m) begins:
n\m  0  1   2    3     4    5     6     7      8    9   10...
0:   1
1:   0  1
2:  -1  0   2
3:   0 -3   0    4
4:   1  0  -8    0     8
5:   0  5   0  -20     0   16
6:  -1  0  18    0   -48    0    32
7:   0 -7   0   56     0 -112     0    64
8:   1  0 -32    0   160    0  -256     0    128
9:   0  9   0 -120     0  432     0  -576      0  256
10: -1  0  50    0  -400    0  1120     0  -1280    0  512
... Reformatted and extended - _Wolfdieter Lang_, Oct 21 2013
E.g., the fourth row (n=3) corresponds to the polynomial T(3,x) = -3*x + 4*x^3.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964. Tenth printing, Wiley, 2002 (also electronically available), p. 795.
  • F. Hirzebruch et al., Manifolds and Modular Forms, Vieweg 1994 pp. 77, 105.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2. ed., Wiley, New York, 1990.
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 22, page 196.
  • TableCurve 2D, Automated curve fitting and equation discovery, Version 5.01 for Windows, User's Manual, Chebyshev Series Polynomials and Rationals, pages 12-21 - 12-24, SYSTAT Software, Inc., Richmond, WA, 2002.

Crossrefs

The first nonzero (sub)diagonal sequences are A011782, -A001792, A001793(n+1), -A001794, A006974, -A006975, A006976, -A209404.

Programs

  • Julia
    using Nemo
    function A053120Row(n)
        R, x = PolynomialRing(ZZ, "x")
        p = chebyshev_t(n, x)
        [coeff(p, j) for j in 0:n] end
    for n in 0:6 A053120Row(n) |> println end # Peter Luschny, Mar 13 2018
    
  • Magma
    &cat[ Coefficients(ChebyshevT(n)): n in [0..11] ]; // Klaus Brockhaus, Mar 08 2008
    
  • Maple
    with(orthopoly) ;
    A053120 := proc(n,k)
        T(n,x) ;
        coeftayl(%,x=0,k) ;
    end proc: # R. J. Mathar, Jun 30 2013
    T := (n, x) -> `if`(n = 0, 1, add((-1)^(n - k) * (n/(2*k))*binomial(k, n - k) *(2*x)^(2*k - n), k = 1 ..n)):
    seq(seq(coeff(T(n, x), x, k), k = 0..n), n = 0..11); # Peter Luschny, Sep 20 2022
  • Mathematica
    t[n_, k_] := Coefficient[ ChebyshevT[n, x], x, k]; Flatten[ Table[ t[n, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, Jan 16 2012 *)
  • PARI
    for(n=0,5,P=polchebyshev(n);for(k=0,n,print1(polcoeff(P,k)", "))) \\ Charles R Greathouse IV, Jan 16 2012
    
  • SageMath
    def f(n,k): # f = A039991
        if (n<2 and k==0): return 1
        elif (k<0 or k>n): return 0
        else: return 2*f(n-1, k) - f(n-2, k-2)
    def A053120(n,k): return f(n, n-k)
    flatten([[A053120(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Aug 10 2022

Formula

T(n, m) = A039991(n, n-m).
G.f. for row polynomials T(n,x) (signed triangle): (1-x*z)/(1-2*x*z+z^2). If unsigned: (1-x*z)/(1-2*x*z-z^2).
T(n, m) := 0 if n < m or n+m odd; T(n, m) = (-1)^(n/2) if m=0 (n even); otherwise T(n, m) = ((-1)^((n+m)/2 + m))*(2^(m-1))*n*binomial((n+m)/2-1, m-1)/m.
Recursion for n >= 2: T(n, m) = T*a(n-1, m-1) - T(n-2, m), T(n, m)=0 if n < m, T(n, -1) := 0, T(0, 0) = T(1, 1) = 1.
G.f. for m-th column (signed triangle): 1/(1+x^2) if m=0, otherwise (2^(m-1))*(x^m)*(1-x^2)/(1+x^2)^(m+1).
From G. C. Greubel, Aug 10 2022: (Start)
Sum_{k=0..floor(n/2)} T(n-k, k) = A000007(n).
T(2*n, n) = i^n * A036909(n/2) * (1+(-1)^n)/2 + [n=0]/3. (End)
T(n, k) = [x^k] T(n, x) for n >= 1, where T(n, x) = Sum_{k=1..n}(-1)^(n - k)*(n/ (2*k))*binomial(k, n - k)*(2*x)^(2*k - n). - Peter Luschny, Sep 20 2022

A370469 Triangle read by columns where T(n,k) is the number of points in Z^n such that |x1| + ... + |xn| = k, |x1|, ..., |xn| > 0.

Original entry on oeis.org

2, 2, 4, 2, 8, 8, 2, 12, 24, 16, 2, 16, 48, 64, 32, 2, 20, 80, 160, 160, 64, 2, 24, 120, 320, 480, 384, 128, 2, 28, 168, 560, 1120, 1344, 896, 256, 2, 32, 224, 896, 2240, 3584, 3584, 2048, 512, 2, 36, 288, 1344, 4032, 8064, 10752, 9216, 4608, 1024
Offset: 1

Views

Author

Shel Kaphan, Mar 30 2024

Keywords

Comments

T(n,k) is the number of points on the n-dimensional cross polytope with facets at distance k from the origin which have no coordinate equal to 0.
T(n,n) = 2^n. The (n-1)-dimensional simplex at distance n from the origin in Z^n has exactly 1 point with no zero coordinates, at (1,1,...,1). There are 2^n (n-1)-dimensional simplexes at distance n from the origin as part of the cross polytope in Z^n. (The lower dimensional facets do not count as they have at least one 0 coordinate.)
T(2*n,3*n) = T(2*n+1,3*n), and this is A036909.

Examples

			 n\k 1 2 3  4  5   6   7    8    9    10    11    12     13     14      15
   -----------------------------------------------------------------------
 1 | 2 2 2  2  2   2   2    2    2     2     2     2      2      2       2
 2 |   4 8 12 16  20  24   28   32    36    40    44     48     52      56
 3 |     8 24 48  80 120  168  224   288   360   440    528    624     728
 4 |       16 64 160 320  560  896  1344  1920  2640   3520   4576    5824
 5 |          32 160 480 1120 2240  4032  6720 10560  15840  22880   32032
 6 |              64 384 1344 3584  8064 16128 29568  50688  82368  128128
 7 |                 128  896 3584 10752 26880 59136 118272 219648  384384
 8 |                      256 2048  9216 30720 84480 202752 439296  878592
 9 |                           512  4608 23040 84480 253440 658944 1537536
10 |                                1024 10240 56320 225280 732160 2050048
11 |                                      2048 22528 135168 585728 2050048
12 |                                            4096  49152 319488 1490944
13 |                                                   8192 106496  745472
14 |                                                         16384  229376
15 |                                                                 32768
The cross polytope in Z^3 (the octahedron) with points at distance 3 from the origin has 8 triangle facets, each with edge length 4. There is one point in the center of each triangle with coordinates (+-1,+-1,+-1).
		

Crossrefs

Cf. A033996, A333714 (n=3)
Cf. A102860 (n=4).
Cf. A036289, A097064, A134401 (+1-diagonal).
Cf. A001815 (+2-diagonal).
Cf. A371064.
Cf. A036909.
2 * A013609.

Programs

  • Mathematica
    T[n_,k_]:=Binomial[k-1,n-1]*2^n; Table[T[n,k],{k,10},{n,k}]//Flatten
  • Python
    from math import comb
    def A370469_T(n,k): return comb(k-1,n-1)<Chai Wah Wu, Apr 25 2024

Formula

T(n,k) = binomial(k-1,n-1)*2^n.
G.f.: 2*x*y/(1 - y - 2*x*y). - Stefano Spezia, Apr 27 2024
Showing 1-2 of 2 results.