|
| |
|
|
A063007
|
|
T(n,k) = binomial(n,k)*binomial(n+k,k), 0 <= k <= n, triangle read by rows.
|
|
61
|
|
|
|
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
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
0,3
|
|
|
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.
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)
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
|
|
|
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.
|
|
|
LINKS
|
|
|
|
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
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
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
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
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
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
|
|
|
EXAMPLE
|
The triangle T(n, k) starts:
n\k 0 1 2 3 4 5 6 7
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
row n = 8: 1 72 1260 9240 34650 72072 84084 51480 12870,
row n = 9: 1 90 1980 18480 90090 252252 420420 411840 218790 48620,
row n = 10: 1 110 2970 34320 210210 756756 1681680 2333760 1969110 923780 184756.
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)
|
|
|
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
|
|
|
PROG
|
(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 */
(Haskell)
a063007 n k = a063007_tabl !! n !! k
a063007_row n = a063007_tabl !! n
a063007_tabl = zipWith (zipWith (*)) a007318_tabl a046899_tabl
(Magma) /* As triangle: */ [[Binomial(n, k)*Binomial(n+k, k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 03 2015
|
|
|
CROSSREFS
|
See A331430 for an essentially identical triangle, except with signed entries.
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.)
|
|
|
KEYWORD
|
|
|
|
AUTHOR
|
|
|
|
STATUS
|
approved
|
| |
|
|
