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

A273167 Numerators of coefficient triangle for expansion of x^(2*n) in terms of Chebyshev polynomials of the first kind T(2*m, x) (A127674).

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 5, 15, 3, 1, 35, 7, 7, 1, 1, 63, 105, 15, 45, 5, 1, 231, 99, 495, 55, 33, 3, 1, 429, 3003, 1001, 1001, 91, 91, 7, 1, 6435, 715, 1001, 273, 455, 35, 15, 1, 1, 12155, 21879, 1989, 4641, 1071, 765, 51, 153, 9, 1, 46189, 20995, 62985, 4845, 4845, 969, 4845, 285, 95, 5, 1
Offset: 0

Views

Author

Wolfdieter Lang, Jun 12 2016

Keywords

Comments

The denominator triangle is given in A273168.
The expansion is x^(2*n) = Sum_{m=0..n} R(n, m)*Tnx(2*m, x), n >= 0, with the rational triangle R(n, m) = a(n, m)/A273168(n, m).
Compare this with A127673.
This is equivalent to the expansion cos(x)^(2n) = Sum_{m=0..n} R(n, m)*cos(2*m*x), n >= 0. Compare this with the even numbered rows of A273496.
See A244420/A244421 for the expansion of x^(2*n+1) in terms of odd-indexed Chebyshev polynomials of the first kind.
The signed rational triangle S(n, m) = R(n, m) * (-1)^m appears in the expansion sin(x)^(2n) = Sum_{m=0..n} S(n, m) * cos(2*m*x), n >= 0. This is equivalent to the identity (1-x^2)^n = Sum_{m=0..n} S(n, m) * T(2*m, x).

Examples

			The triangle a(n, m) begins:
n\m     0     1    2    3    4   5  6   7 8 9
0:      1
1:      1     1
2:      3     1    1
3:      5    15    3    1
4:     35     7    7    1    1
5:     63   105   15   45    5   1
6:    231    99  495   55   33   3  1
7:    429  3003 1001 1001   91  91  7   1
8:   6435   715 1001  273  455  35 15   1 1
9:  12155 21879 1989 4641 1071 765 51 153 9 1
...
The rational triangle R(n, m) begins:
n\m  0       1      2     3      4     5  ...
0:   1
1:  1/2     1/2
2:  3/8     1/2    1/8
3:  5/16   15/32   3/16  1/32
4: 35/128   7/16   7/32  1/16  1/128
5: 63/256 105/256 15/64 45/512 5/256 1/512
...
row 6: 231/1024 99/256 495/2048 55/512 33/1024 3/512 1/2048,
row 7: 429/2048 3003/8192 1001/4096 1001/8192 91/2048 91/8192 7/4096 1/8192,
row 8: 6435/32768 715/2048 1001/4096 273/2048 455/8192 35/2048 15/4096 1/2048 1/32768,
row 9: 12155/65536 21879/65536 1989/8192 4641/32768 1071/16384 765/32768 51/8192 153/131072 9/65536 1/131072,
...
n=3: x^6 = (5/16)*T(0, x) + (15/32)*T(2, x)
  +(3/16)*T(4, x) + (1/32)*T(6,x).
  cos^6(x) = (5/16) + (15/32)*cos(2*x) +
    (3/16)*cos(4*x) + (1/32)*cos(6*x).
  sin^6(x) = (5/16) - (15/32)*cos(2*x) +
    (3/16)*cos(4*x) - (1/32)*cos(6*x).
		

Crossrefs

Programs

  • Mathematica
    T[MaxN_] := Function[{n}, With[{exp = Expand[(1/2)^(2 n) (Exp[I x] + Exp[-I x])^(2 n)]}, Prepend[ 2 Coefficient[exp, Exp[I 2 # x]] & /@ Range[1, n], exp /. {Exp[_] -> 0}]]][#] & /@ Range[0, MaxN];
    T[5] // ColumnForm
    T2[MaxN_] := Table[Inverse[Outer[Coefficient[#1, x, #2] &, Prepend[ChebyshevT[#, x] & /@ Range[2 MaxN], 1], Range[0, 2 MaxN]]][[n, m]], {n, 1, 2 MaxN, 2}, {m, 1, n, 2}]
    T2[6] // ColumnForm (* Bradley Klee, Jun 14 2016 *)
  • PARI
    a(n, m) = if (m == 0, numerator((1/2^(2*n-1)) * binomial(2*n,n)/2), numerator((1/2^(2*n-1))*binomial(2*n, n-m)));
    tabl(nn) = for (n=0, nn, for (k=0, n, print1(a(n,k), ", ")); print()); \\ Michel Marcus, Jun 19 2016

Formula

a(n, m) = numerator(R(n, m)), n >= 0, m = 1, ..., n, with the rationals R(n, m) given by R(n, 0) = (1/2^(2*n-1))*binomial(2*n,n)/2 and R(n ,m) = (1/2^(2*n-1))*binomial(2*n, n-m) for m =1..n, n >= 0.

A273168 Denominators of coefficient triangle for expansion of x^(2*n) in terms of Chebyshev polynomials of the first kind T(2*m, x) (A127674).

Original entry on oeis.org

1, 2, 2, 8, 2, 8, 16, 32, 16, 32, 128, 16, 32, 16, 128, 256, 256, 64, 512, 256, 512, 1024, 256, 2048, 512, 1024, 512, 2048, 2048, 8192, 4096, 8192, 2048, 8192, 4096, 8192, 32768, 2048, 4096, 2048, 8192, 2048, 4096, 2048, 32768, 65536, 65536, 8192, 32768, 16384, 32768, 8192, 131072, 65536, 131072, 262144, 65536, 262144, 32768, 65536, 32768, 524288, 131072, 262144, 131072, 524288
Offset: 0

Views

Author

Wolfdieter Lang, Jun 12 2016

Keywords

Comments

The numerator sequence is given in A273167, where details are given.

Examples

			The triangle a(n, m) begins:
n\m     0    1    2    3    4    5    6    7
0:      1
1:      2    2
2:      8    2    8
3:     16   32   16   32
4:    128   16   32   16  128
5:    256  256   64  512  256  512
6:   1024  256 2048  512 1024  512 2048
7:   2048 8192 4096 8192 2048 8192 4096 8192
...
row 8: 32768 2048 4096 2048 8192 2048 4096 2048 32768,
row 9: 65536 65536 8192 32768 16384 32768 8192 131072 65536 131072,
...
		

Crossrefs

Cf. A273167.

Programs

  • PARI
    a(n, m) = if (m == 0, denominator((1/2^(2*n-1)) * binomial(2*n,n)/2), denominator((1/2^(2*n-1))*binomial(2*n, n-m)));
    tabl(nn) = for (n=0, nn, for (k=0, n, print1(a(n,k), ", ")); print()); \\ Michel Marcus, Jun 19 2016

Formula

a(n, m) = denominator(R(n, m)), n >= 0, m = 1, ..., n, with the rationals R(n, m) given by R(n, 0) = (1/2^(2*n-1)) * binomial(2*n,n)/2 and R(n ,m) = (1/2^(2*n-1))*binomial(2*n, n-m) for m =1..n, n >= 0.

A063007 T(n,k) = binomial(n,k)*binomial(n+k,k), 0 <= k <= n, triangle read by rows.

Original entry on oeis.org

1, 1, 2, 1, 6, 6, 1, 12, 30, 20, 1, 20, 90, 140, 70, 1, 30, 210, 560, 630, 252, 1, 42, 420, 1680, 3150, 2772, 924, 1, 56, 756, 4200, 11550, 16632, 12012, 3432, 1, 72, 1260, 9240, 34650, 72072, 84084, 51480, 12870, 1, 90, 1980, 18480, 90090, 252252, 420420, 411840, 218790, 48620
Offset: 0

Views

Author

Henry Bottomley, Jul 02 2001

Keywords

Comments

T(n,k) is the number of compatible k-sets of cluster variables in Fomin and Zelevinsky's Cluster algebra of finite type B_n. Take a row of this triangle regarded as a polynomial in x and rewrite as a polynomial in y := x+1. The coefficients of the polynomial in y give a row of triangle A008459 (squares of binomial coefficients). For example, x^2+6*x+6 = y^2+4*y+1. - Paul Boddington, Mar 07 2003
T(n,k) is the number of lattice paths from (0,0) to (n,n) using steps E=(1,0), N=(0,1) and D=(1,1) (i.e., bilateral Schroeder paths), having k N=(0,1) steps. E.g. T(2,0)=1 because we have DD; T(2,1) = 6 because we have NED, NDE, EDN, END, DEN and DNE; T(2,2)=6 because we have NNEE, NENE, NEEN, EENN, ENEN and ENNE. - Emeric Deutsch, Apr 20 2004
Another version of [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] DELTA [0, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...] = 1; 1, 0; 1, 2, 0; 1, 6, 6, 0; 1, 12, 30, 20, 0; ..., where DELTA is the operator defined in A084938. - Philippe Deléham Apr 15 2005
Terms in row n are the coefficients of the Legendre polynomial P(n,2x+1) with increasing powers of x.
From Peter Bala, Oct 28 2008: (Start)
Row n of this triangle is the f-vector of the simplicial complex dual to an associahedron of type B_n (a cyclohedron) [Fomin & Reading, p.60]. See A008459 for the corresponding h-vectors for associahedra of type B_n and A001263 and A033282 respectively for the h-vectors and f-vectors for associahedra of type A_n.
An alternative description of this triangle in terms of f-vectors is as follows. Let A_n be the root lattice generated as a monoid by {e_i - e_j: 0 <= i,j <= n+1}. Let P(A_n) be the polytope formed by the convex hull of this generating set. Then the rows of this array are the f-vectors of a unimodular triangulation of P(A_n) [Ardila et al.]. A008459 is the corresponding array of h-vectors for these type A_n polytopes. See A127674 (without the signs) for the array of f-vectors for type C_n polytopes and A108556 for the array of f-vectors associated with type D_n polytopes.
The S-transform on the ring of polynomials is the linear transformation of polynomials that is defined on the basis monomials x^k by S(x^k) = binomial(x,k) = x(x-1)...(x-k+1)/k!. Let P_n(x) denote the S-transform of the n-th row polynomial of this array. In the notation of [Hetyei] these are the Stirling polynomials of the type B associahedra. The first few values are P_1(x) = 2*x + 1, P_2(x) = 3*x^2 + 3*x + 1 and P_3(x) = (10*x^3 + 15*x^2 + 11*x + 3)/3. These polynomials have their zeros on the vertical line Re x = -1/2 in the complex plane, that is, the polynomials P_n(-x) satisfy a Riemann hypothesis. See A142995 for further details. The sequence of values P_n(k) for k = 0,1,2,3, ... produces the n-th row of A108625. (End)
This is the row reversed version of triangle A104684. - Wolfdieter Lang, Sep 12 2016
T(n, k) is also the number of (n-k)-dimensional faces of a convex n-dimensional Lipschitz polytope of real functions f defined on the set X = {1, 2, ..., n+1} which satisfy the condition f(n+1) = 0 (see Gordon and Petrov). - Stefano Spezia, Sep 25 2021
The rows seem to give (up to sign) the coefficients in the expansion of the integer-valued polynomial ((x+1)*(x+2)*(x+3)*...*(x+n) / n!)^2 in the basis made of the binomial(x+i,i). - F. Chapoton, Oct 09 2022
Chapoton's observation above is correct: the precise expansion is ((x+1)*(x+2)*(x+3)*...*(x+n)/ n!)^2 = Sum_{k = 0..n} (-1)^k*T(n,n-k)*binomial(x+2*n-k, 2*n-k), as can be verified using the WZ algorithm. For example, n = 3 gives ((x+1)*(x+2)*(x+3)/3!)^2 = 20*binomial(x+6,6) - 30*binomial(x+5,5) + 12*binomial(x+4,4) - binomial(x+3,3). - Peter Bala, Jun 24 2023

Examples

			The triangle T(n, k) starts:
  n\k 0  1    2     3     4      5      6      7      8     9
  0:  1
  1:  1  2
  2:  1  6    6
  3:  1 12   30    20
  4:  1 20   90   140    70
  5:  1 30  210   560   630    252
  6:  1 42  420  1680  3150   2772    924
  7:  1 56  756  4200 11550  16632  12012   3432
  8:  1 72 1260  9240 34650  72072  84084  51480  12870
  9:  1 90 1980 18480 90090 252252 420420 411840 218790 48620
... reformatted by _Wolfdieter Lang_, Sep 12 2016
From _Petros Hadjicostas_, Jul 11 2020: (Start)
Its inverse (from Table II, p. 92, in Ser's book) is
   1;
  -1/2,  1/2;
   1/3, -1/2,    1/6;
  -1/4,  9/20,  -1/4,   1/20;
   1/5, -2/5,    2/7,  -1/10,  1/70;
  -1/6,  5/14, -25/84,  5/36, -1/28,  1/252;
   1/7, -9/28,  25/84, -1/6,   9/154, -1/84, 1/924;
   ... (End)
		

References

  • J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p. 366.
  • J. Ser, Les Calculs Formels des Séries de Factorielles. Gauthier-Villars, Paris, 1933, Table I, p. 92.
  • D. Zagier, Integral solutions of Apery-like recurrence equations, in: Groups and Symmetries: from Neolithic Scots to John McKay, CRM Proc. Lecture Notes 47, Amer. Math. Soc., Providence, RI, 2009, pp. 349-366.

Crossrefs

See A331430 for an essentially identical triangle, except with signed entries.
Columns include A000012, A002378, A033487 on the left and A000984, A002457, A002544 on the right.
Main diagonal is A006480.
Row sums are A001850. Alternating row sums are A033999.
Cf. A033282 (f-vectors type A associahedra), A108625, A080721 (f-vectors type D associahedra).
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692,A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Haskell
    a063007 n k = a063007_tabl !! n !! k
    a063007_row n = a063007_tabl !! n
    a063007_tabl = zipWith (zipWith (*)) a007318_tabl a046899_tabl
    -- Reinhard Zumkeller, Nov 18 2014
    
  • Magma
    /* As triangle: */ [[Binomial(n,k)*Binomial(n+k,k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 03 2015
  • Maple
    p := (n,x) -> orthopoly[P](n,1+2*x): seq(seq(coeff(p(n,x),x,k), k=0..n), n=0..9);
  • Mathematica
    Flatten[Table[Binomial[n, k]Binomial[n + k, k], {n, 0, 10}, {k, 0, n}]] (* Harvey P. Dale, Dec 24 2011 *)
    Table[CoefficientList[Hypergeometric2F1[-n, n + 1, 1, -x], x], {n, 0, 9}] // Flatten
    (* Peter Luschny, Mar 09 2018 *)
  • PARI
    {T(n, k) = local(t); if( n<0, 0, t = (x + x^2)^n; for( k=1, n, t=t'); polcoeff(t, k) / n!)} /* Michael Somos, Dec 19 2002 */
    
  • PARI
    {T(n, k) = binomial(n, k) * binomial(n+k, k)} /* Michael Somos, Sep 22 2013 */
    
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, (n+k)! / (k!^2 * (n-k)!))} /* Michael Somos, Sep 22 2013 */
    

Formula

T(n, k) = (n+k)!/(k!^2*(n-k)!) = T(n-1, k)*(n+k)/(n-k) = T(n, k-1)*(n+k)*(n-k+1)/k^2 = T(n-1, k-1)*(n+k)*(n+k-1)/k^2.
binomial(x, n)^2 = Sum_{k>=0} T(n,k) * binomial(x, n+k). - Michael Somos, May 11 2012
T(n, k) = A109983(n, k+n). - Michael Somos, Sep 22 2013
G.f.: G(t, z) = 1/sqrt(1-2*z-4*t*z+z^2). Row generating polynomials = P_n(1+2z), i.e., T(n, k) = [z^k] P_n(1+2*z), where P_n are the Legendre polynomials. - Emeric Deutsch, Apr 20 2004
Sum_{k>=0} T(n, k)*A000172(k) = Sum_{k>=0} T(n, k)^2 = A005259(n). - Philippe Deléham, Jun 08 2005
1 + z*d/dz(log(G(t,z))) = 1 + (1 + 2*t)*z + (1 + 8*t + 8*t^2)*z^2 + ... is the o.g.f. for a signed version of A127674. - Peter Bala, Sep 02 2015
If R(n,t) denotes the n-th row polynomial then x^3 * exp( Sum_{n >= 1} R(n,t)*x^n/n ) = x^3 + (1 + 2*t)*x^4 + (1 + 5*t + 5*t^2)*x^5 + (1 + 9*t + 21*t^2 + 14*t^3)*x^6 + ... is an o.g.f for A033282. - Peter Bala, Oct 19 2015
P(n,x) := 1/(1 + x)*Integral_{t = 0..x} R(n,t) dt are (modulo differences of offset) the row polynomials of A033282. - Peter Bala, Jun 23 2016
From Peter Bala, Mar 09 2018: (Start)
R(n,x) = Sum_{k = 0..n} binomial(2*k,k)*binomial(n+k,n-k)*x^k.
R(n,x) = Sum_{k = 0..n} binomial(n,k)^2*x^k*(1 + x)^(n-k).
n*R(n,x) = (1 + 2*x)*(2*n - 1)*R(n-1,x) - (n - 1)*R(n-2,x).
R(n,x) = (-1)^n*R(n,-1 - x).
R(n,x) = 1/n! * (d/dx)^n ((x^2 + x)^n). (End)
The row polynomials are R(n,x) = hypergeom([-n, n + 1], [1], -x). - Peter Luschny, Mar 09 2018
T(n,k) = C(n+1,k)*A009766(n,k). - Bob Selcoe, Jan 18 2020 (Connects this triangle with the Catalan triangle. - N. J. A. Sloane, Jan 18 2020)
If we let A(n,k) = (-1)^(n+k)*(2*k+1)*(n*(n-1)*...*(n-(k-1)))/((n+1)*...*(n+(k+1))) for n >= 0 and k = 0..n, and we consider both T(n,k) and A(n,k) as infinite lower triangular arrays, then they are inverses of one another. (Empty products are by definition 1.) See the example below. The rational numbers |A(n,k)| appear in Table II on p. 92 in Ser's (1933) book. - Petros Hadjicostas, Jul 11 2020
From Peter Bala, Nov 28 2021: (Start)
Row polynomial R(n,x) = Sum_{k >= n} binomial(k,n)^2 * x^(k-n)/(1+x)^(k+1) for x > -1/2.
R(n,x) = 1/(1 + x)^(n+1) * hypergeom([n+1, n+1], [1], x/(1 + x)).
R(n,x) = (1 + x)^n * hypergeom([-n, -n], [1], x/(1 + x)).
R(n,x) = hypergeom([(n+1)/2, -n/2], [1], -4*x*(1 + x)).
If we set R(-1,x) = 1, we can run the recurrence n*R(n,x) = (1 + 2*x)*(2*n - 1)*R(n-1,x) - (n - 1)*R(n-2,x) backwards to give R(-n,x) = R(n-1,x).
R(n,x) = [t^n] ( (1 + t)*(1 + x*(1 + t)) )^n. (End)
n*T(n,k) = (2*n-1)*T(n-1,k) + (4*n-2)*T(n-1,k-1) - (n-1)*T(n-2,k). - Fabián Pereyra, Jun 30 2022
From Peter Bala, Oct 07 2024: (Start)
n-th row polynomial R(n,x) = Sum_{k = 0..n} binomial(n, k) * x^k o (1 + x)^(n-k), where o denotes the black diamond product of power series as defined by Dukes and White (see Bala, Section 4.4, exercise 3).
Denote this triangle by T. Then T * transpose(T) = A143007, the square array of crystal ball sequences for the A_n X A_n lattices.
Let S denote the triangle ((-1)^(n+k)*T(n, k))n,k >= 0, a signed version of this triangle. Then S^(-1) * T = A007318, Pascal's triangle; it appears that T * S^(-1) = A110098.
T = A007318 * A115951. (End)

A086645 Triangle read by rows: T(n, k) = binomial(2n, 2k).

Original entry on oeis.org

1, 1, 1, 1, 6, 1, 1, 15, 15, 1, 1, 28, 70, 28, 1, 1, 45, 210, 210, 45, 1, 1, 66, 495, 924, 495, 66, 1, 1, 91, 1001, 3003, 3003, 1001, 91, 1, 1, 120, 1820, 8008, 12870, 8008, 1820, 120, 1, 1, 153, 3060, 18564, 43758, 43758, 18564, 3060, 153, 1, 1, 190, 4845, 38760
Offset: 0

Views

Author

Philippe Deléham, Jul 26 2003

Keywords

Comments

Terms have the same parity as those of Pascal's triangle.
Coefficients of polynomials (1/2)*((1 + x^(1/2))^(2n) + (1 - x^(1/2))^(2n)).
Number of compositions of 2n having k parts greater than 1; example: T(3, 2) = 15 because we have 4+2, 2+4, 3+2+1, 3+1+2, 2+3+1, 2+1+3, 1+3+2, 1+2+3, 2+2+1+1, 2+1+2+1, 2+1+1+2, 1+2+2+1, 1+2+1+2, 1+1+2+2, 3+3. - Philippe Deléham, May 18 2005
Number of binary words of length 2n - 1 having k runs of consecutive 1's; example: T(3,2) = 15 because we have 00101, 01001, 01010, 01011, 01101, 10001, 10010, 10011, 10100, 10110, 10111, 11001, 11010, 11011, 11101. - Philippe Deléham, May 18 2005
Let M_n be the n X n matrix M_n(i, j) = T(i, j-1); then for n > 0, det(M_n) = A000364(n), Euler numbers; example: det([1, 1, 0, 0; 1, 6, 1, 0; 1, 15, 15, 1; 1, 28, 70, 28 ]) = 1385 = A000364(4). - Philippe Deléham, Sep 04 2005
Equals ConvOffsStoT transform of the hexagonal numbers, A000384: (1, 6, 15, 28, 45, ...); e.g., ConvOffs transform of (1, 6, 15, 28) = (1, 28, 70, 28, 1). - Gary W. Adamson, Apr 22 2008
From Peter Bala, Oct 23 2008: (Start)
Let C_n be the root lattice generated as a monoid by {+-2*e_i: 1 <= i <= n; +-e_i +- e_j: 1 <= i not equal to j <= n}. Let P(C_n) be the polytope formed by the convex hull of this generating set. Then the rows of this array are the h-vectors of a unimodular triangulation of P(C_n) [Ardila et al.]. See A127674 for (a signed version of) the corresponding array of f-vectors for these type C_n polytopes. See A008459 for the array of h-vectors for type A_n polytopes and A108558 for the array of h-vectors associated with type D_n polytopes.
The Hilbert transform of this triangle is A142992 (see A145905 for the definition of this term).
(End)
Diagonal sums: A108479. - Philippe Deléham, Sep 08 2009
Coefficients of Product_{k=1..n} (cot(k*Pi/(2n+1))^2 - x) = Sum_{k=0..n} (-1)^k*binomial(2n,2k)*x^k/(2n+1-2k). - David Ingerman (daviddavifree(AT)gmail.com), Mar 30 2010
Generalized Narayana triangle for 4^n (or cosh(2x)). - Paul Barry, Sep 28 2010
Coefficients of the matrix inverse appear to be T^(-1)(n,k) = (-1)^(n+k)*A086646(n,k). - R. J. Mathar, Mar 12 2013
Let E(y) = Sum_{n>=0} y^n/(2*n)! = cosh(sqrt(y)). Then this triangle is the generalized Riordan array (E(y), y) with respect to the sequence (2*n)! as defined in Wang and Wang. Cf. A103327. - Peter Bala, Aug 06 2013
Row 6, (1,66,495,924,495,66,1), plays a role in expansions of powers of the Dedekind eta function. See the Chan link, p. 534, and A034839. - Tom Copeland, Dec 12 2016

Examples

			From _Peter Bala_, Oct 23 2008: (Start)
The triangle begins
n\k|..0.....1.....2.....3.....4.....5.....6
===========================================
0..|..1
1..|..1.....1
2..|..1.....6.....1
3..|..1....15....15.....1
4..|..1....28....70....28.....1
5..|..1....45...210...210....45.....1
6..|..1....66...495...924...495....66.....1
...
(End)
From _Peter Bala_, Aug 06 2013: (Start)
Viewed as the generalized Riordan array (cosh(sqrt(y)), y) with respect to the sequence (2*n)! the column generating functions begin
1st col: cosh(sqrt(y)) = 1 + y/2! + y^2/4! + y^3/6! + y^4/8! + ....
2nd col: 1/2!*y*cosh(sqrt(y)) = y/2! + 6*y^2/4! + 15*y^3/6! + 28*y^4/8! + ....
3rd col: 1/4!*y^2*cosh(sqrt(y)) = y^2/4! + 15*y^3/6! + 70*y^4/8! + 210*y^5/10! + .... (End)
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 224.

Crossrefs

Cf. A008459, A108558, A127674, A142992. - Peter Bala, Oct 23 2008
Cf. A103327 (binomial(2n+1, 2k+1)), A103328 (binomial(2n, 2k+1)), A091042 (binomial(2n+1, 2k)). -Wolfdieter Lang, Jan 06 2013
Cf. A086646 (unsigned matrix inverse), A103327.
Cf. A034839.

Programs

  • Magma
    /* As triangle: */ [[Binomial(2*n, 2*k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Dec 14 2016
  • Maple
    A086645:=(n,k)->binomial(2*n,2*k): seq(seq(A086645(n,k),k=0..n),n=0..12);
  • Mathematica
    Table[Binomial[2 n, 2 k], {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Dec 13 2016 *)
  • Maxima
    create_list(binomial(2*n,2*k),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
    
  • PARI
    {T(n, k) = binomial(2*n, 2*k)};
    
  • PARI
    {T(n, k) = sum( i=0, min(k, n-k), 4^i * binomial(n, 2*i) * binomial(n - 2*i, k-i))}; /* Michael Somos, May 26 2005 */
    

Formula

T(n, k) = (2*n)!/((2*(n-k))!*(2*k)!) row sums = A081294. COLUMNS: A000012, A000384
Sum_{k>=0} T(n, k)*A000364(k) = A000795(n) = (2^n)*A005647(n).
Sum_{k>=0} T(n, k)*2^k = A001541(n). Sum_{k>=0} T(n, k)*3^k = 2^n*A001075(n). Sum_{k>=0} T(n, k)*4^k = A083884(n). - Philippe Deléham, Feb 29 2004
O.g.f.: (1 - z*(1+x))/(x^2*z^2 - 2*x*z*(1+z) + (1-z)^2) = 1 + (1 + x)*z +(1 + 6*x + x^2)*z^2 + ... . - Peter Bala, Oct 23 2008
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A081294(n), A001541(n), A090965(n), A083884(n), A099140(n), A099141(n), A099142(n), A165224(n), A026244(n) for x = 0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Sep 08 2009
Product_{k=1..n} (cot(k*Pi/(2n+1))^2 - x) = Sum_{k=0..n} (-1)^k*binomial(2n,2k)*x^k/(2n+1-2k). - David Ingerman (daviddavifree(AT)gmail.com), Mar 30 2010
From Paul Barry, Sep 28 2010: (Start)
G.f.: 1/(1-x-x*y-4*x^2*y/(1-x-x*y)) = (1-x*(1+y))/(1-2*x*(1+y)+x^2*(1-y)^2);
E.g.f.: exp((1+y)*x)*cosh(2*sqrt(y)*x);
T(n,k) = Sum_{j=0..n} C(n,j)*C(n-j,2*(k-j))*4^(k-j). (End)
T(n,k) = 2*T(n-1,k) + 2*T(n-1,k-1) + 2*T(n-2,k-1) - T(n-2,k) - T(n-2,k-2), with T(0,0)=T(1,0)=T(1,1)=1, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Nov 26 2013
From Peter Bala, Sep 22 2021: (Start)
n-th row polynomial R(n,x) = (1-x)^n*T(n,(1+x)/(1-x)), where T(n,x) is the n-th Chebyshev polynomial of the first kind. Cf. A008459.
R(n,x) = Sum_{k = 0..n} binomial(n,2*k)*(4*x)^k*(1 + x)^(n-2*k).
R(n,x) = n*Sum_{k = 0..n} (n+k-1)!/((n-k)!*(2*k)!)*(4*x)^k*(1-x)^(n-k) for n >= 1. (End)

A084930 Triangle of coefficients of Chebyshev polynomials T_{2n+1} (x).

Original entry on oeis.org

1, -3, 4, 5, -20, 16, -7, 56, -112, 64, 9, -120, 432, -576, 256, -11, 220, -1232, 2816, -2816, 1024, 13, -364, 2912, -9984, 16640, -13312, 4096, -15, 560, -6048, 28800, -70400, 92160, -61440, 16384, 17, -816, 11424, -71808, 239360, -452608, 487424, -278528, 65536, -19, 1140, -20064, 160512, -695552
Offset: 0

Views

Author

Gary W. Adamson, Jun 12 2003

Keywords

Comments

From Herb Conn, Jan 28 2005: (Start)
"Letting x = 2 Cos 2A, we have the following trigonometric identities:
"Sin 3A = 3*Sin A - 4*Sin^3 A
"Sin 5A = 5*Sin A - 20*Sin^3 A + 16*Sin^5 A
"Sin 7A = 7*Sin A - 56*Sin^3 A + 112*Sin^5 A - 64*Sin^7 A
"Sin 9A = 9*Sin A - 120*Sin^3 A + 432*Sin^5 A - 576*Sin^7 A + 256*Sin^9 A, etc." (End)
Cayley (1876) states "Write sin u = x, then we have sin u = x, [...] sin 3u = 3x - 4x^3, [...] sin 5u = 5x - 20x^3 + 16 x^5, [...]". Since T_n(cos(u)) = cos(nu) for all integer n, sin(u) = cos(u - Pi/2), and sin(u + k*Pi) = (-1)^k sin(u) it follows that T_n(sin(u)) = (-1)^((n-1)/2) sin(nu) for all odd integer n. - Michael Somos, Jun 19 2012
From Wolfdieter Lang, Aug 05 2014: (Start)
The coefficient triangle t(n,k) for the row polynomials Todd(n, x) := T_{2*n+1}(sqrt(x))/sqrt(x) = sum(t(n,k)*x^k, k=0..n) is the Riordan triangle ((1-z)/(1+z)^2, 4*z/(1+z)^2) (rewrite the g.f. for the present triangle a(n,k) given in the formula section). The triangle entries t(n,k) = a(n,k), but the interpretation of the row polynomials is different for both cases.
From the relation Todd(n, x) = S(n, 2*(2*x-1)) - S(n-1, 2*(2*x-1)) with the Chebyshev S-polynomials (see A049310 and the formula section of A130777) follows the recurrence: Todd(n, x) = 2*(-1)^n*(1-x)*Todd(n-1, 1-x) + (2x-1)*Todd(n-1, x), n >= 1, Todd(0, x) = 1.
This corresponds to the triangle recurrence t(n,k) = (2*(k+1)*(-1)^(n-k) - 1)*t(n-1,k) + 2*(1 +(-1)^(n-k))*t(n-1,k-1) + 2*(-1)^(n-k)*sum(binomial(l+1,k)*t(n-1,l), l=k+1..n-1), n >= k >= 1, t(n,k) = 0 if n < k, t(n,0) = (-1)^n*(2*n+1). Compare this with the shorter recurrence involving the rational A-sequence for this Riordan triangle which has g.f. x^2/(2-x-2*sqrt(1-x)). t(n,k) = sum(A(j)*t(n-1,k-1+j), j=0..n-k), n >= k >= 1. The Z-sequence has g.f. -(1 + 2/sqrt(1-x)). For the A- and Z-sequence see a link under A006232. (End)

Examples

			The triangle a(n,k):
n   2n+1\k 0     1      2       3       4        5        6         7        8        9      10 ...
0    1:    1
1    3:   -3     4
2    5:    5   -20     16
3    7:   -7    56   -112      64
4    9:    9  -120    432    -576     256
5   11:  -11   220  -1232    2816   -2816     1024
6   13:   13  -364   2912   -9984   16640   -13312     4096
7   15:  -15   560  -6048   28800  -70400    92160   -61440     16384
8   17:   17  -816  11424  -71808  239360  -452608   487424   -278528    65536
9   19:  -19  1140 -20064  160512 -695552  1770496 -2723840   2490368 -1245184   262144
10  21:   21 -1540  33264 -329472 1793792 -5870592 12042240 -15597568 12386304 -5505024 1048576
... formatted and extended by _Wolfdieter Lang_, Aug 02 2014.
---------------------------------------------------------------------------------------------------
First few polynomials T_{2n+1}(x) are
1*x - 3*x + 4*x^3
5*x - 20*x^3 + 16*x^5
- 7*x + 56*x^3 - 112*x^5 + 64*x^7
9*x - 120*x^3 + 432*x^5 - 576*x^7  + 256*x^9
		

References

  • A. Cayley, On an Expression for 1 +- sin(2p+1)u in Terms of sin u, Messenger of Mathematics, 5 (1876), pp. 7-8 = Mathematical Papers Vol. 10, n. 630, pp. 1-2.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2nd ed., Wiley, New York, 1990. p. 37, eq. (1.96) and p. 4, eq. (1.10).

Crossrefs

Cf. A053120 (coefficient triangle of T-polynomials), A127674 (even-indexed T polynomials).
Cf. A127675 (row reversed triangle with different signs).

Programs

  • Mathematica
    row[n_] := (T = ChebyshevT[2*n+1, x]; Coefficient[T, x, #]& /@ Range[1, Exponent[T, x], 2]); Table[row[n], {n, 0, 9} ] // Flatten (* Jean-François Alcover, Aug 06 2014 *)

Formula

Alternate rows of A008310.
a(n,k)=((-1)^(n-k))*(2^(2*k))*binomial(n+1+k,2*k+1)*(2*n+1)/(n+1+k) if n>=k>=0 else 0.
From Wolfdieter Lang, Aug 02 2014: (Start)
a(n,k) = [x^(2*k+1)] T_{2*n+1}(x), n >= k >= 0.
G.f. for row polynomials: x*(1-z)/(1 + 2*(1- 2*x^2)*z + z^2). (End)
The first column sequences are: A157142, 4*(-1)^(n+1)*A000330(n), 16*(-1)^n*A005585(n-1), 64*(-1)^(n+1)*A050486(n-3), 256*(-1)^n*A054333(n-4), ... - Wolfdieter Lang, Aug 05 2014

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)lycos.com), Jun 26 2003
Edited; example rewritten (to conform with the triangle) by Wolfdieter Lang, Aug 02 2014

A084534 Triangle read by rows: row #n has n+1 terms. T(n,0)=1, T(n,n)=2, T(n,m) = T(n-1,m-1) + Sum_{k=0..m} T(n-1-k,m-k).

Original entry on oeis.org

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

Views

Author

Gary W. Adamson, May 29 2003

Keywords

Comments

Sum of row #n = A000204(2n). (But sum of row #0 = 1.)
Row #n has the unsigned coefficients of the monic polynomial whose roots are 2 cos(Pi*(2k-1)/(4n)) for k=1..2n. [Comment corrected by Barry Brent, Jan 03 2006]
The positive roots are some diagonal lengths of a regular (4n)-gon, inscribed in the unit circle.
Polynomial of row #n = Sum_{m=0..n} (-1)^m * T(n,m) x^(2*n-2*m).
This is the unsigned version of the coefficient table for scaled Chebyshev T(2*n,x) polynomials. - Wolfdieter Lang, Mar 07 2007
Reversed A127677 (cf. A156308, A217476, A263916). - Tom Copeland, Nov 07 2015

Examples

			First few Chebyshev T(2*n,x) polynomials:
  T(2*0,x) = 1;
  T(2*1,x) = x^2 -   2;
  T(2*2,x) = x^4 -   4*x^2 +  2;
  T(2*3,x) = x^6 -   6*x^4 +  9*x^2 -  2;
  T(2*4,x) = x^8 -   8*x^6 + 20*x^4 - 16*x^2 +  2;
  T(2*5,x) = x^10 - 10*x^8 + 35*x^6 - 50*x^4 + 25*x^2 - 2;
Triangle begins as:
  1;
  1,  2;
  1,  4,  2;
  1,  6,  9,   2;
  1,  8, 20,  16,   2;
  1, 10, 35,  50,  25,  2;
  1, 12, 54, 112, 105, 36, 2;
		

References

  • I. Kaplansky and J. Riordan, The problème des ménages, Scripta Math. 12, (1946), 113-124. See p. 118.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2. ed., Wiley, New York, 1990. p. 37, eq.(1.96) and p. 4. eq.(1.10).

Crossrefs

Row sums are A005248 for n > 0.
Companion triangle A082985.
Cf. A082985 (unsigned scaled coefficient table for Chebyshev's T(2*n+1, x) polynomials).

Programs

  • Magma
    A084534:= func< n,k | k eq 0 select 1 else 2*(n/k)*Binomial(2*n-k-1, k-1) >;
    [A084534(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Feb 02 2022
    
  • Maple
    T := proc(n, m): if n=0 then 1 else binomial(2*n-m, m)*2*n/(2*n-m) fi: end: seq(seq(T(n,m),m=0..n),n=0..10); # Johannes W. Meijer, May 31 2018
  • Mathematica
    a[n_, m_] := Binomial[2n-m, m]*2n/(2n-m); a[0, 0] = 1; Table[a[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (* Jean-François Alcover, Apr 12 2016, after Wolfdieter Lang *)
  • PARI
    T(n,m) = if(n==0, m==0, binomial(2*n-m, m)*2*n/(2*n-m)) \\ Andrew Howroyd, Dec 18 2017
    
  • Sage
    def A084534(n,k): return 1 if (k==0) else 2*(n/k)*binomial(2*n-k-1, k-1)
    flatten([[A084534(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Feb 02 2022

Formula

T(n,m) = binomial(2*n-m, m)*2*n/(2*n-m) for n > 0. - Andrew Howroyd, Dec 18 2017
Signed version from Wolfdieter Lang, Mar 07 2007: (Start)
a(n,m)=0 if n
a(n,m)=0 if n
a(n,m)=0 if nA127674(n,n-m)/2^(2*(n-m)-1) (scaled coefficients of Chebyshev's T(2*n,x), decreasing even powers). [corrected by Johannes W. Meijer, May 31 2018] (End)

Extensions

Edited by Don Reble, Nov 12 2005

A108556 Triangle, read by rows, where row n equals the inverse binomial transform of the crystal ball sequence for D_n lattice.

Original entry on oeis.org

1, 1, 2, 1, 4, 4, 1, 12, 30, 20, 1, 24, 120, 192, 96, 1, 40, 330, 940, 1080, 432, 1, 60, 732, 3200, 6240, 5568, 1856, 1, 84, 1414, 8708, 25200, 37184, 27104, 7744, 1, 112, 2480, 20352, 80960, 173824, 206080, 126976, 31744, 1, 144, 4050, 42588, 221544, 643824, 1096032, 1085760, 579456, 128768
Offset: 0

Author

Paul D. Hanna, Jun 10 2005

Keywords

Comments

Row n equals the inverse binomial transform of row n of the square array A108553.
Array of f-vectors for type D root polytopes [Ardila et al.]. See A063007 and A127674 for the arrays of f-vectors for type A and type C root polytopes respectively. - Peter Bala, Oct 23 2008

Examples

			Triangle begins:
1;
1,2;
1,4,4;
1,12,30,20;
1,24,120,192,96;
1,40,330,940,1080,432;
1,60,732,3200,6240,5568,1856;
1,84,1414,8708,25200,37184,27104,7744;
1,112,2480,20352,80960,173824,206080,126976,31744; ...
		

Crossrefs

Cf. A108553, A108557 (row sums), A108558, Rows are inverse binomial transforms of: A001844 (row 2), A005902 (row 3), A007204 (row 4), A008356 (row 5), A008358 (row 6), A008360 (row 7), A008362 (row 8), A008377 (row 9), A008379 (row 10).

Programs

  • Mathematica
    T[n_, k_] := Module[{A}, A = Table[Table[If[r - 1 == 0 || c - 1 == 0, 1, If[r - 1 == 1, 2c - 1, Sum[Binomial[r + c - j - 2, c - j - 1] (Binomial[2r - 2, 2j] - 2(r - 1) Binomial[r - 3, j - 1]), {j, 0, c - 1}]]], {c, 1, n + 1}], {r, 1, n + 1}]; SeriesCoefficient[((A[[n + 1]]. x^Range[0, n]) /. x -> x/(1 + x))/(1 + x), {x, 0, k}]];
    Table[T[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 26 2018, from PARI *)
  • PARI
    T(n,k)=local(A=vector(n+1,r,vector(n+1,c,if(r-1==0 || c-1==0,1,if(r-1==1,2*c-1, sum(j=0,c-1,binomial(r+c-j-2,c-j-1)*(binomial(2*r-2,2*j)-2*(r-1)*binomial(r-3,j-1)))))))); polcoeff(subst(Ser(A[n+1]),x,x/(1+x))/(1+x),k)

Formula

Main diagonal equals A008353: 2^(n-1)*(2^n-n) for n>1.
O.g.f. : rational function N(x,z)/D(x,z), where N(x,z) = 1 - 3*(1 + 2*x)*z + (3 + 8*x + 8*x^2)*z^2 - (1 + 2*x)*(1 - 6*x - 6*x^2)z^3 - 8*x*(1 + x)(1 + 2*x + 2*x^2)*z^4 + 2*x*(1 + x)*(1 + 2*x)*z^5 and D(x,z) = ((1-z)^2 - 4*x*z)*(1 - z*(1 + 2*x))^2. - Peter Bala, Oct 23 2008

A075733 Array of coefficients in Zagier's polynomials P_(n,0)(x).

Original entry on oeis.org

1, -2, 1, -8, 8, 1, -18, 48, -32, 1, -32, 160, -256, 128, 1, -50, 400, -1120, 1280, -512, 1, -72, 840, -3584, 6912, -6144, 2048, 1, -98, 1568, -9408, 26880, -39424, 28672, -8192, 1, -128, 2688, -21504, 84480, -180224, 212992, -131072, 32768, 1, -162, 4320, -44352, 228096, -658944, 1118208
Offset: 1

Author

Benoit Cloitre, Oct 18 2002

Keywords

Comments

Zagier's polynomial P_{n,m} (m<=n and m>=0) are used under PARI/GP (for algorithmic acceleration of alternating or positive series)
Up to signs, coefficients of Chebyshev's T-polynomials for even index. See A127674.

Crossrefs

Cf. A000332.

Formula

P_(n, 0)(x) = (-1)^n*2^(2*n-1)*Product_{i=1..n} (x - cos(Pi*(2*i-1)/(4*n))^2).

A004141 Norm of a matrix.

Original entry on oeis.org

1, 2, 8, 48, 256, 1280, 6912, 39424, 212992, 1118208, 6553600, 36765696, 199229440, 1133117440, 6499598336, 36175872000, 200655503360, 1167945891840, 6620826304512, 36681168191488, 212364657950720, 1219998345330688, 6864598984556544, 38958828003262464, 226089827240509440, 1287455960675123200
Offset: 1

Keywords

Comments

Maximum in the row n-1 of the absolute values of the triangle A127674. - R. J. Mathar, Jul 15 2015

References

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

Crossrefs

Cf. A259868.

Programs

  • Maple
    A := proc(n,k)
        2*n/(n+k)*binomial(n+k,n-k)*2^(2*k-1) ;
    end proc:
    A004141 := proc(n)
        seq(abs(A(n,k)),k=0..n-1) ;
        max(%) ;
    end proc:
    seq(A004141(n),n=1..30) ; # R. J. Mathar, Jul 15 2015
  • Mathematica
    A[n_, k_] := 2n/(n+k) Binomial[n+k, n-k] 2^(2k-1);
    row[n_] := Table[A[n, k], {k, 0, n-1}] // Abs;
    a[n_] := If[n < 3, n, row[n-1] // Max];
    Array[a, 30] (* Jean-François Alcover, Mar 30 2020, from Maple *)

Formula

The scanned page from Kammler and McGlinn (page 843 of the article) gives a fairly explicit way to calculate a(n). - N. J. A. Sloane, Jul 09 2015

A305402 A number triangle T(n,k) read by rows for 0<=k<=n, related to the Taylor expansion of f(u, p) = (1/2)*(1+1/(sqrt(1-u^2)))*exp(p*sqrt(1-u^2)).

Original entry on oeis.org

1, 1, -2, 3, -4, 2, 15, -18, 9, -2, 105, -120, 60, -16, 2, 945, -1050, 525, -150, 25, -2, 10395, -11340, 5670, -1680, 315, -36, 2, 135135, -145530, 72765, -22050, 4410, -588, 49, -2, 2027025, -2162160, 1081080, -332640, 69300, -10080, 1008, -64, 2
Offset: 0

Author

Johannes W. Meijer, May 31 2018

Keywords

Comments

The function f(u, p) = (1/2)*(1+1/(sqrt(1-u^2))) * exp(p*sqrt(1-u^2)) was found while studying the Fresnel-Kirchhoff and the Rayleigh-Sommerfeld theories of diffraction, see the Meijer link.
The Taylor expansion of f(u, p) leads to the number triangle T(n, k), see the example section.
Normalization of the triangle terms, dividing the T(n, k) by T(n-k, 0), leads to A084534.
The row sums equal A003436, n >= 2, respectively A231622, n >= 1.

Examples

			The first few terms of the Taylor expansion of f(u; p) are:
f(u, p) = exp(p) * (1 + (1-2*p) * u^2/4 + (3-4*p+2*p^2) * u^4/16 + (15-18*p+9*p^2-2*p^3) * u^6/96 + (105-120*p+60*p^2-16*p^3+2*p^4) * u^8/768 + ... )
The first few rows of the T(n, k) triangle are:
n=0:     1
n=1:     1,     -2
n=2:     3,     -4,    2
n=3:    15,    -18,    9,    -2
n=4:   105,   -120,   60,   -16,   2
n=5:   945,  -1050,  525,  -150,  25,  -2
n=6: 10395, -11340, 5670, -1680, 315, -36, 2
		

References

  • J. W. Goodman, Introduction to Fourier Optics, 1996.
  • A. Papoulis, Systems and Transforms with Applications in Optics, 1968.

Crossrefs

Cf. Related to the left hand columns: A001147, A001193, A261065.
Cf. Related to the right hand columns: A280560, A162395, A006011, A040977, A053347, A054334, A266561.

Programs

  • Magma
    [[n le 0 select 1 else (-1)^k*2^(k-n+1)*Factorial(2*n-k-1)*Binomial(n, k)/Factorial(n-1): k in [0..n]]: n in [1..10]]; // G. C. Greubel, Nov 08 2018
  • Maple
    T := proc(n, k): if n=0 then 1 else (-1)^k*2^(k-n+1)*n*(2*n-k-1)!/(k!*(n-k)!) fi: end: seq(seq(T(n, k), k=0..n), n=0..8);
  • Mathematica
    Table[If[n==0 && k==0,1, (-1)^k*2^(k-n+1)*n*(2*n-k-1)!/(k!*(n-k)!)], {n, 0, 10}, {k,0,n}]//Flatten (* G. C. Greubel, Nov 08 2018 *)
  • PARI
    T(n,k) = {if(n==0, 1, (-1)^k*2^(k-n+1)*n*(2*n-k-1)!/(k!*(n-k)!))}
    for(n=0, 10, for(k=0, n, print1(T(n, k), ", ")); print); \\ Andrew Howroyd, Nov 08 2018
    

Formula

T(n, k) = (-1)^k*2^(k-n+1)*n*(2*n-k-1)!/(k!*(n-k)!), n > 0 and 0 <= k <= n, T(0, 0) = 1.
T(n, k) = (-1)^k*A001147(n-k)*A084534(n, k), n >= 0 and 0 <= k <= n.
T(n, k) = 2^(2*(k-n)+1)*A001147(n-k)*A127674(n, n-k), n > 0 and 0 <= k <= n, T(0, 0) = 1.
T(n, k) = (-1)^k*(A001497(n, k) + A132062(n, k)), n >= 1, T(0,0) = 1.
Showing 1-10 of 11 results. Next