A053122 Triangle of coefficients of Chebyshev's S(n,x-2) = U(n,x/2-1) polynomials (exponents of x in increasing order).
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
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.
Links
- T. D. Noe, Rows n=0..50 of triangle, flattened
- Wolfdieter Lang, First rows of the triangle.
- M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
- P. Barry, Symmetric Third-Order Recurring Sequences, Chebyshev Polynomials, and Riordan Arrays, JIS 12 (2009) 09.8.6.
- Naiomi T. Cameron and Asamoah Nkwanta, On some (pseudo) involutions in the Riordan group, J. of Integer Sequences, 8(2005), 1-16.
- P. Damianou, On the characteristic polynomials of Cartan matrices and Chebyshev polynomials, arXiv preprint arXiv:1110.6620 [math.RT], 2014 (p. 10). - _Tom Copeland_, Oct 11 2014
- Pentti Haukkanen, Jorma Merikoski, Seppo Mustonen, Some polynomials associated with regular polygons, Acta Univ. Sapientiae, Mathematica, 6, 2 (2014) 178-193.
- Eric Weisstein's World of Mathematics, Cartan Matrix
- Eric Weisstein's World of Mathematics, Dynkin Diagram
- Index entries for sequences related to Chebyshev polynomials.
Crossrefs
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)
A052179 Triangle of numbers arising in enumeration of walks on cubic lattice.
1, 4, 1, 17, 8, 1, 76, 50, 12, 1, 354, 288, 99, 16, 1, 1704, 1605, 700, 164, 20, 1, 8421, 8824, 4569, 1376, 245, 24, 1, 42508, 48286, 28476, 10318, 2380, 342, 28, 1, 218318, 264128, 172508, 72128, 20180, 3776, 455, 32, 1, 1137400, 1447338
Offset: 0
Comments
Triangle T(n,k), 0 <= k <= n, read by rows given by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = 4*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + 4*T(n-1,k) + T(n-1,k+1) for k >= 1. - Philippe Deléham, Mar 27 2007
Triangle read by rows: T(n,k) = number of lattice paths from (0,0) to (n,k) that do not go below the line y=0 and consist of steps U=(1,1), D=(1,-1) and four types of steps H=(1,0); example: T(3,1)=50 because we have UDU, UUD, 16 HHU paths, 16 HUH paths and 16 UHH paths. - Philippe Deléham, Sep 25 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Riordan array ((1-4x-sqrt(1-8x+12x^2))/(2x^2), (1-4x-sqrt(1-8x+12x^2))/(2x)). Inverse of A159764. - Paul Barry, Apr 21 2009
6^n = (n-th row terms) dot (first n+1 terms in (1,2,3,...)). Example: 6^3 = 216 = (76, 50, 12, 1) dot (1, 2, 3, 4) = (76 + 100 + 36 + 4) = 216. - Gary W. Adamson, Jun 15 2011
A subset of the "family of triangles" (Deléham comment of Sep 25 2007) is the succession of binomial transforms beginning with triangle A053121, (0,0); giving -> A064189, (1,1); -> A039598, (2,2); -> A091965, (3,3); -> A052179, (4,4); -> A125906, (5,5) ->, etc.; generally the binomial transform of the triangle generated from (n,n) = that generated from ((n+1),(n+1)). - Gary W. Adamson, Aug 03 2011
Examples
Triangle begins: 1; 4, 1; 17, 8, 1; 76, 50, 12, 1; 354, 288, 99, 16, 1; ... Production matrix begins: 4, 1; 1, 4, 1; 0, 1, 4, 1; 0, 0, 1, 4, 1; 0, 0, 0, 1, 4, 1; 0, 0, 0, 0, 1, 4, 1; 0, 0, 0, 0, 0, 1, 4, 1; - _Philippe Deléham_, Nov 04 2011
Links
- G. C. Greubel, Table of n, a(n) for the first 100 rows, flattened
- Rigoberto Flórez, Leandro Junes, José L. Ramírez, Further Results on Paths in an n-Dimensional Cubic Lattice, Journal of Integer Sequences, Vol. 21 (2018), Article 18.1.2.
- R. K. Guy, Catwalks, Sandsteps and Pascal Pyramids, J. Integer Seqs., Vol. 3 (2000), #00.1.6.
Programs
-
Maple
T:= proc(n, k) option remember; `if`(min(n, k)<0, 0, `if`(max(n, k)=0, 1, T(n-1, k-1)+4*T(n-1, k)+T(n-1, k+1))) end: seq(seq(T(n,k), k=0..n), n=0..10); # Alois P. Heinz, Oct 28 2021
-
Mathematica
t[0, 0] = 1; t[n_, k_] /; k < 0 || k > n = 0; t[n_, 0] := t[n, 0] = 4*t[n-1, 0] + t[n-1, 1]; t[n_, k_] := t[n, k] = t[n-1, k-1] + 4*t[n-1, k] + t[n-1, k+1]; Flatten[ Table[t[n, k], {n, 0, 9}, {k, 0, n}]] (* Jean-François Alcover, Oct 10 2011, after Philippe Deleham *)
Formula
Sum_{k>=0} T(m, k)*T(n, k) = T(m+n, 0) = A005572(m+n). - Philippe Deléham, Sep 15 2005
n-th row = M^n * V, where M = the infinite tridiagonal matrix with all 1's in the super and subdiagonals and (4,4,4,...) in the main diagonal. E.g., Row 3 = (76, 50, 12, 1) since M^3 * V = [76, 50, 12, 1, 0, 0, 0, ...]. - Gary W. Adamson, Nov 04 2006
Sum_{k=0..n} T(n,k) = A005573(n). - Philippe Deléham, Feb 04 2007
Sum_{k=0..n} T(n,k)*(k+1) = 6^n. - Philippe Deléham, Mar 27 2007
Sum_{k=0..n} T(n,k)*x^k = A033543(n), A064613(n), A005572(n), A005573(n) for x = -2, -1, 0, 1 respectively. - Philippe Deléham, Nov 28 2009
As an infinite lower triangular matrix = the binomial transform of A091965 and 4th binomial transform of A053121. - Gary W. Adamson, Aug 03 2011
G.f.: 2/(1 - 4*x - 2*x*y + sqrt(1 - 8*x + 12*x^2)). - Daniel Checa, Aug 17 2022
G.f. for the m-th column: x^m*(A(x))^(m+1), where A(x) is the g.f. of the sequence counting the walks on the cubic lattice starting and finishing on the xy plane and never going below it (A005572). Explicitly, the g.f. is x^m*((1 - 4*x - sqrt(1 - 8*x + 12*x^2))/(2*x^2))^(m+1). - Daniel Checa, Aug 28 2022
A008315 Catalan triangle read by rows. Also triangle of expansions of powers of x in terms of Chebyshev polynomials U_n(x).
1, 1, 1, 1, 1, 2, 1, 3, 2, 1, 4, 5, 1, 5, 9, 5, 1, 6, 14, 14, 1, 7, 20, 28, 14, 1, 8, 27, 48, 42, 1, 9, 35, 75, 90, 42, 1, 10, 44, 110, 165, 132, 1, 11, 54, 154, 275, 297, 132, 1, 12, 65, 208, 429, 572, 429, 1, 13, 77, 273, 637, 1001, 1001, 429, 1, 14, 90, 350, 910, 1638, 2002, 1430, 1, 15, 104
Offset: 0
Comments
Number of standard tableaux of shape (n-k,k) (0<=k<=floor(n/2)). Example: T(4,1)=3 because in th top row we can have 124, 134, or 123 (but not 234). - Emeric Deutsch, May 23 2004
T(n,k) is the number of n-digit binary words (length n sequences on {0,1}) containing k 1's such that no initial segment of the sequence has more 1's than 0's. - Geoffrey Critzer, Jul 31 2009
T(n,k) is the number of dispersed Dyck paths (i.e. Motzkin paths with no (1,0) steps at positive heights) of length n and having k (1,1)-steps. Example: T(5,1)=4 because, denoting U=(1,1), D=(1,-1), H=(1,0), we have HHHUD, HHUDH, HUDHH, and UDHHH. - Emeric Deutsch, May 30 2011
T(n,k) is the number of length n left factors of Dyck paths having k (1,-1)-steps. Example: T(5,1)=4 because, denoting U=(1,1), D=(1,-1), we have UUUUD, UUUDU, UUDUU, and UDUUU. There is a simple bijection between length n left factors of Dyck paths and dispersed Dyck paths of length n, that takes D steps into D steps. - Emeric Deutsch, Jun 19 2011
Triangle, with zeros omitted, given by (1, 0, 0, -1, 1, 0, 0, -1, 1, 0, 0, -1, 1, ...) DELTA (0, 1, -1, 0, 0, 1, -1, 0, 0, 1, -1, 0, 0, 1, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 12 2011
T(n,k) are rational multiples of A055151(n,k). - Peter Luschny, Oct 16 2015
T(2*n,n) = Sum_{k>=0} T(n,k)^2 = A000108(n), T(2*n+1,n) = A000108(n+1). - Michael Somos, Jun 08 2020
Examples
Triangle begins: 1; 1; 1, 1; 1, 2; 1, 3, 2; 1, 4, 5; 1, 5, 9, 5; 1, 6, 14, 14; 1, 7, 20, 28, 14; ... T(5,2) = 5 because there are 5 such sequences: {0, 0, 0, 1, 1}, {0, 0, 1, 0, 1}, {0, 0, 1, 1, 0}, {0, 1, 0, 0, 1}, {0, 1, 0, 1, 0}. - _Geoffrey Critzer_, Jul 31 2009
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. 796.
- P. J. Larcombe, A question of proof..., Bull. Inst. Math. Applic. (IMA), 30, Nos. 3/4, 1994, 52-54.
Links
- T. D. Noe, Rows n=0..100 of triangle, flattened
- M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
- Tewodros Amdeberhan, Moa Apagodu, and Doron Zeilberger, Wilf's "Snake Oil" Method Proves an Identity in The Motzkin Triangle, arXiv:1507.07660 [math.CO], 2015.
- Nantel Bergeron, Kelvin Chan, Yohana Solomon, Farhad Soltani, and Mike Zabrocki, Quasisymmetric harmonics of the exterior algebra, arXiv:2206.02065 [math.CO], 2022.
- Chassidy Bozeman, Christine Cheng, Pamela E. Harris, Stephen Lasinis, and Shanise Walker, The Pinnacle Sets of a Graph, arXiv:2406.19562 [math.CO], 2024. See pp. 9-10.
- Suyoung Choi and Hanchul Park, A new graph invariant arises in toric topology, arXiv preprint arXiv:1210.3776 [math.AT], 2012.
- Suyuong Choi and Younghan Yoon, A decomposition of graph a-numbers, arXiv:2508.06855 [math.CO], 2025. See p. 14.
- C. Kenneth Fan, Structure of a Hecke algebra quotient, J. Amer. Math. Soc. 10 (1997), no. 1, 139-167.
- R. K. Guy, Catwalks, sandsteps and Pascal pyramids, J. Integer Sequences, Vol. 3 (2000), Article #00.1.6.
- Lin Jiu, Victor H. Moll, and C. Vignat, Identities for generalized Euler polynomials, arXiv:1401.8037 [math.PR], 2014.
- Nik Lygeros and Oliver Rozier, A new solution to the equation tau(rho) == 0 (mod p), J. Int. Seq. 13 (2010) # 10.7.4.
- Mustafa A. A. Obaid, S. Khalid Nauman, Wafaa M. Fakieh, and Claus Michael Ringel, The numbers of support-tilting modules for a Dynkin algebra, 2014 and J. Int. Seq. 18 (2015) 15.10.6.
- Alon Regev, The central component of a triangulation, Journal of Integer Sequences, Vol. 16 (2013), Article 13.4.1, p. 7.
- J. Riordan, The distribution of crossings of chords joining pairs of 2n points on a circle, Math. Comp., 29 (1975), 215-222.
- J. Riordan, The distribution of crossings of chords joining pairs of 2n points on a circle, Math. Comp., 29 (1975), 215-222. [Annotated scanned copy]
- L. W. Shapiro, A Catalan triangle, Discrete Math. 14 (1976), no. 1, 83-90. [Annotated scanned copy]
- Zheng Shi, Impurity entropy of junctions of multiple quantum wires, arXiv preprint arXiv:1602.00068 [cond-mat.str-el], 2016 (See Appendix A).
- Index entries for sequences related to Chebyshev polynomials.
Crossrefs
Programs
-
Haskell
a008315 n k = a008315_tabf !! n !! k a008315_row n = a008315_tabf !! n a008315_tabf = map reverse a008313_tabf -- Reinhard Zumkeller, Nov 14 2013
-
Maple
b:= proc(x, y) option remember; `if`(y<0 or y>x, 0, `if`(x=0, 1, add(b(x-1, y+j), j=[-1, 1]))) end: T:= (n, k)-> b(n, n-2*k): seq(seq(T(n, k), k=0..n/2), n=0..16); # Alois P. Heinz, Oct 14 2022
-
Mathematica
Table[Binomial[k, i]*(k - 2 i + 1)/(k - i + 1), {k, 0, 20}, {i, 0, Floor[k/2]}] // Grid (* Geoffrey Critzer, Jul 31 2009 *)
-
PARI
{T(n, k) = if( k<0 || k>n\2, 0, if( n==0, 1, T(n-1, k-1) + T(n-1, k)))}; /* Michael Somos, Aug 17 1999 */
Formula
T(n, 0) = 1 if n >= 0; T(2*k, k) = T(2*k-1, k-1) if k>0; T(n, k) = T(n-1, k-1) + T(n-1, k) if k=1, 2, ..., floor(n/2). - Michael Somos, Aug 17 1999
T(n, k) = binomial(n, k) - binomial(n, k-1). - Michael Somos, Aug 17 1999
Rows of Catalan triangle A008313 read backwards. Sum_{k>=0} T(n, k)^2 = A000108(n); A000108 : Catalan numbers. - Philippe Deléham, Feb 15 2004
T(n,k) = C(n,k)*(n-2*k+1)/(n-k+1). - Geoffrey Critzer, Jul 31 2009
Extensions
Expanded description from Clark Kimberling, Jun 15 1997
A089942 Inverse binomial matrix applied to A039599.
1, 0, 1, 1, 1, 1, 1, 3, 2, 1, 3, 6, 6, 3, 1, 6, 15, 15, 10, 4, 1, 15, 36, 40, 29, 15, 5, 1, 36, 91, 105, 84, 49, 21, 6, 1, 91, 232, 280, 238, 154, 76, 28, 7, 1, 232, 603, 750, 672, 468, 258, 111, 36, 8, 1, 603, 1585, 2025, 1890, 1398, 837, 405, 155, 45, 9, 1, 1585, 4213, 5500
Offset: 0
Comments
Triangle T(n,k), 0 <= k <= n, defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = T(n-1,1), T(n,k) = T(n-1,k-1) + T(n-1,k) + T(n-1,k+1) for k >= 1. - Philippe Deléham, Feb 27 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise from choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Riordan array (f(x),x*g(x)), where f(x)is the o.g.f. of A005043 and g(x)is the o.g.f. of A001006. - Philippe Deléham, Nov 22 2009
Riordan array ((1+x-sqrt(1-2x-3x^2))/(2x(1+x)), (1-x-sqrt(1-2x-3x^2))/(2x)). Inverse of Riordan array ((1+x)/(1+x+x^2),x/(1+x+x^2)). E.g.f. of column k is exp(x)*(Bessel_I(k,2x)-Bessel_I(k+1,2x)).
Diagonal sums are A187306.
Simultaneous equations using the first n rows solve for diagonal lengths of odd N = (2n+1) regular polygons, with constants c^0, c^1, c^2, ...; where c = 1 + 2*cos( 2*Pi/N) = sin(3*Pi/N)/sin(Pi/N) = the third longest diagonal of N>5. By way of example, take the first 4 rows relating to the 9-gon (nonagon), N=(2*4 + 1), with c = 1 + 2*cos(2*Pi/9) = 2.5320888.... The simultaneous equations are (1,0,0,0) = 1; (0,1,0,0) = c; (1,1,1,0) = c^2, (1,3,2,1) = c^3. The answers are 1, 2.532..., 2.879..., and 1.879...; the four distinct diagonal lengths of the 9-gon (nonagon) with edge = 1. - Gary W. Adamson, Sep 07 2011
Number of linearly independent irreducible representations of a given weight j in a tensor of given rank n. - Mikkel N. Schmidt, Aug 20 2025
Examples
Triangle begins 1, 0, 1, 1, 1, 1, 1, 3, 2, 1, 3, 6, 6, 3, 1, 6, 15, 15, 10, 4, 1, 15, 36, 40, 29, 15, 5, 1, 36, 91, 105, 84, 49, 21, 6, 1, 91, 232, 280, 238, 154, 76, 28, 7, 1 Production matrix is 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1
Links
- G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
- Paul Barry and Aoife Hennessy, Four-term Recurrences, Orthogonal Polynomials and Riordan Arrays, Journal of Integer Sequences, 2012, article 12.4.2. - From _N. J. A. Sloane_, Sep 21 2012
- J. A. R. Coope, R. F. Snider, and F. R. McCourt, Irreducible Cartesian Tensors, The Journal of Chemical Physics, vol 43, no 7, 1965.
- Emeric Deutsch, Luca Ferrari, and Simone Rinaldi, Production Matrices, Advances in Applied Mathematics, 34 (2005) pp. 101-122.
- Toufik Mansour and Mark Shattuck, Enumeration of Catalan and smooth words according to capacity, Integers (2025) Vol. 25, Art. No. A5. See p. 29.
- Donatella Merlini, Douglas G. Rogers, Renzo Sprugnoli, and M. Cecilia Verri, On some alternative characterizations of Riordan arrays, Canad J. Math., 49 (1997), 301-320.
- Yidong Sun and Luping Ma, Minors of a class of Riordan arrays related to weighted partial Motzkin paths, Eur. J. Comb. (2014) Vol. 39, 157-169. See Table 2.2.
Crossrefs
Row sums give A002426 (central trinomial coefficients).
Programs
-
Maple
T:= (n,k) -> simplify(GegenbauerC(n-k,-n+1,-1/2)-GegenbauerC(n-k-1,-n+1,-1/2)): for n from 1 to 9 do seq(T(n,k), k=1..n) od; # Peter Luschny, May 12 2016 # Or by recurrence: T := proc(n, k) option remember; if n = k then 1 elif k < 0 or n < 0 or k > n then 0 elif k = 0 then T(n-1, 1) else T(n-1, k-1) + T(n-1, k) + T(n-1, k+1) fi end: for n from 0 to 9 do seq(T(n, k), k = 0..n) od; # Peter Luschny, May 25 2021
-
Mathematica
T[n_, k_] := GegenbauerC[n - k, -n + 1, -1/2] - GegenbauerC[n - k - 1, -n + 1, -1/2]; Table[T[n, k], {n,1,10}, {k,1,n}] // Flatten (* G. C. Greubel, Feb 28 2017 *)
Formula
G.f.: (1+z-q)/[(1+z)(2z-t+tz+tq)], where q = sqrt(1-2z-3z^2).
Sum_{k>=0} T(m,k)*T(n,k) = T(m+n,0) = A005043(m+n). - Philippe Deléham, Mar 22 2007
Sum_{k=0..n} T(n,k)*(2k+1) = 3^n. - Philippe Deléham, Mar 22 2007
Sum_{k=0..n} T(n,k)*2^k = A112657(n). - Philippe Deléham, Apr 01 2007
T(n,2k) + T(n,2k+1) = A109195(n,k). - Philippe Deléham, Nov 11 2008
T(n,k) = GegenbauerC(n-k,-n+1,-1/2) - GegenbauerC(n-k-1,-n+1,-1/2) for 1 <= k <= n. - Peter Luschny, May 12 2016
Extensions
Edited by Emeric Deutsch, Mar 04 2004
A091965 Triangle read by rows: T(n,k) = number of lattice paths from (0,0) to (n,k) that do not go below the line y=0 and consist of steps U=(1,1), D=(1,-1) and three types of steps H=(1,0) (left factors of 3-Motzkin steps).
1, 3, 1, 10, 6, 1, 36, 29, 9, 1, 137, 132, 57, 12, 1, 543, 590, 315, 94, 15, 1, 2219, 2628, 1629, 612, 140, 18, 1, 9285, 11732, 8127, 3605, 1050, 195, 21, 1, 39587, 52608, 39718, 19992, 6950, 1656, 259, 24, 1, 171369, 237129, 191754, 106644, 42498, 12177, 2457
Offset: 0
Comments
Reversal of A084536. - Philippe Deléham, Mar 23 2007
Triangle T(n,k), 0 <= k <= n, read by rows given by T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = 3*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + 3*T(n-1,k) + T(n-1,k+1) for k >= 1. - Philippe Deléham, Mar 27 2007
This triangle belongs to the family of triangles defined by T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
5^n = (n-th row terms) dot (first n+1 terms in (1,2,3,...)). Example for row 4: 5^4 = 625 = (137, 132, 57, 12, 1) dot (1, 2, 3, 4, 5) = (137 + 264 + 171 + 48 + 5) = 625. - Gary W. Adamson, Jun 15 2011
Riordan array ((1-3*x-sqrt(1-6*x+5*x^2))/(2*x^2), (1-3*x-sqrt(1-6*x+5*x^2))/(2*x)). - Philippe Deléham, Feb 19 2012
Examples
Triangle begins: 1; 3, 1; 10, 6, 1; 36, 29, 9, 1; 137, 132, 57, 12, 1; 543, 590, 315, 94, 15, 1; 2219, 2628, 1629, 612, 140, 18, 1; T(3,1)=29 because we have UDU, UUD, 9 HHU paths, 9 HUH paths and 9 UHH paths. Production matrix begins 3, 1; 1, 3, 1; 0, 1, 3, 1; 0, 0, 1, 3, 1; 0, 0, 0, 1, 3, 1; 0, 0, 0, 0, 1, 3, 1; 0, 0, 0, 0, 0, 1, 3, 1; 0, 0, 0, 0, 0, 0, 1, 3, 1; 0, 0, 0, 0, 0, 0, 0, 1, 3, 1; 0, 0, 0, 0, 0, 0, 0, 0, 1, 3, 1; - _Philippe Deléham_, Nov 07 2011
References
- A. Nkwanta, Lattice paths and RNA secondary structures, DIMACS Series in Discrete Math. and Theoretical Computer Science, 34, 1997, 137-147.
Links
- Vincenzo Librandi, Rows n = 0..100, flattened
- Shu-Chiuan Chang and Robert Shrock, Structure of the Partition Function and Transfer Matrices for the Potts Model in a Magnetic Field on Lattice Strips, J. Stat. Physics 137 (2009) 667, table 5.
- Helmut Prodinger, The amplitude of Motzkin paths, arXiv:2104.07596 [math.CO], 2021. Mentions this sequence.
- Helmut Prodinger, Multi-edge trees and 3-coloured Motzkin paths: bijective issues, arXiv:2105.03350 [math.CO], 2021.
- Helmut Prodinger, Weighted unary-binary trees, Hex-trees, marked ordered trees, and related structures, arXiv:2106.14782 [math.CO], 2021.
Programs
-
Mathematica
nmax = 9; t[n_, k_] := ((k+1)*n!*Hypergeometric2F1[k+3/2, k-n, 2k+3, -4]) / ((k+1)!*(n-k)!); Flatten[ Table[ t[n, k], {n, 0, nmax}, {k, 0, n}]] (* Jean-François Alcover, Nov 14 2011, after Vladimir Kruchinin *) T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0, T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]]; Table[T[n, k, 3, 3], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, May 22 2017 *)
-
Maxima
T(n,k):=(k+1)*sum((binomial(2*(m+1),m-k)*binomial(n,m))/(m+1),m,k,n); /* Vladimir Kruchinin, Oct 08 2011 */
-
Sage
@CachedFunction def A091965(n,k): if n==0 and k==0: return 1 if k<0 or k>n: return 0 if k==0: return 3*A091965(n-1,0)+A091965(n-1,1) return A091965(n-1,k-1)+3*A091965(n-1,k)+A091965(n-1,k+1) for n in (0..7): [A091965(n,k) for k in (0..n)] # Peter Luschny, Nov 05 2012
Formula
G.f.: G = 2/(1 - 3*z - 2*t*z + sqrt(1-6*z+5*z^2)). Alternatively, G = M/(1 - t*z*M), where M = 1 + 3*z*M + z^2*M^2.
Sum_{k>=0} T(m, k)*T(n, k) = T(m+n, 0) = A002212(m+n+1). - Philippe Deléham, Sep 14 2005
The triangle may also be generated from M^n * [1,0,0,0,...], where M = an infinite tridiagonal matrix with 1's in the super and subdiagonals and [3,3,3,...] in the main diagonal. - Gary W. Adamson, Dec 17 2006
Sum_{k=0..n} T(n,k)*(k+1) = 5^n. - Philippe Deléham, Mar 27 2007
Sum_{k=0..n} T(n,k)*x^k = A117641(n), A033321(n), A007317(n), A002212(n+1), A026378(n+1) for x = -3, -2, -1, 0, 1 respectively. - Philippe Deléham, Nov 28 2009
T(n,k) = (k+1)*Sum_{m=k..n} binomial(2*(m+1),m-k)*binomial(n,m)/(m+1). - Vladimir Kruchinin, Oct 08 2011
The n-th row polynomial R(n,x) equals the n-th degree Taylor polynomial of the function (1 - x^2)*(1 + 3*x + x^2)^n expanded about the point x = 0. - Peter Bala, Sep 06 2022
A038622 Triangular array that counts rooted polyominoes.
1, 2, 1, 5, 3, 1, 13, 9, 4, 1, 35, 26, 14, 5, 1, 96, 75, 45, 20, 6, 1, 267, 216, 140, 71, 27, 7, 1, 750, 623, 427, 238, 105, 35, 8, 1, 2123, 1800, 1288, 770, 378, 148, 44, 9, 1, 6046, 5211, 3858, 2436, 1296, 570, 201, 54, 10, 1, 17303, 15115, 11505, 7590, 4302, 2067, 825, 265
Offset: 0
Comments
The PARI program gives any row k and any n-th term for this triangular array in square or right triangle array format. - Randall L Rathbun, Jan 20 2002
Triangle T(n,k), 0 <= k <= n, read by rows given by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = 2*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + T(n-1,k) + T(n-1,k+1) for k >= 1. - Philippe Deléham, Mar 27 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Triangle read by rows = partial sums of A064189 terms starting from the right. - Gary W. Adamson, Oct 25 2008
Column k has e.g.f. exp(x)*(Bessel_I(k,2x)+Bessel_I(k+1,2x)). - Paul Barry, Mar 08 2011
Examples
From _Paul Barry_, Mar 08 2011: (Start) Triangle begins 1; 2, 1; 5, 3, 1; 13, 9, 4, 1; 35, 26, 14, 5, 1; 96, 75, 45, 20, 6, 1; 267, 216, 140, 71, 27, 7, 1; 750, 623, 427, 238, 105, 35, 8, 1; 2123, 1800, 1288, 770, 378, 148, 44, 9, 1; Production matrix is 2, 1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1 (End)
Links
- Reinhard Zumkeller, Rows n = 0..120 of triangle, flattened
- D. Gouyou-Beauchamps and G. Viennot, Equivalence of the two-dimensional directed animal problem to a one-dimensional path problem, Adv. in Appl. Math. 9 (1988), no. 3, 334-357. See Table 1 on page 340.
Crossrefs
Programs
-
Haskell
import Data.List (transpose) a038622 n k = a038622_tabl !! n !! k a038622_row n = a038622_tabl !! n a038622_tabl = iterate (\row -> map sum $ transpose [tail row ++ [0,0], row ++ [0], [head row] ++ row]) [1] -- Reinhard Zumkeller, Feb 26 2013
-
Maple
T := (n,k) -> simplify(GegenbauerC(n-k,-n+1,-1/2)+GegenbauerC(n-k-1,-n+1,-1/2)): for n from 1 to 9 do seq(T(n,k),k=1..n) od; # Peter Luschny, May 12 2016
-
Mathematica
nmax = 10; t[n_ /; n > 0, k_ /; k >= 1] := t[n, k] = t[n-1, k-1] + t[n-1, k] + t[n-1, k+1]; t[0, 0] = 1; t[0, ] = 0; t[?Negative, ?Negative] = 0; t[n, 0] := 2 t[n-1, 0] + t[n-1, 1]; Flatten[ Table[ t[n, k], {n, 0, nmax}, {k, 0, n}]](* Jean-François Alcover, Nov 09 2011 *)
-
PARI
s=[0,1]; {A038622(n,k)=if(n==0,1,t=(2*(n+k)*(n+k-1)*s[2]+3*(n+k-1)*(n+k-2)*s[1])/((n+2*k-1)*n); s[1]=s[2]; s[2]=t; t)}
Formula
a(n, k) = a(n-1, k-1) + a(n-1, k) + a(n-1, k+1) for k>0, a(n, k) = 2*a(n-1, k) + a(n-1, k+1) for k=0.
Riordan array ((sqrt(1-2x-3x^2)+3x-1)/(2x(1-3x)),(1-x-sqrt(1-2x-3x^2))/(2x)). Inverse of Riordan array ((1-x)/(1+x+x^2),x/(1+x+x^2)). First column is A005773(n+1). Row sums are 3^n (A000244). If L=A038622, then L*L' is the Hankel matrix for A005773(n+1), where L' is the transpose of L. - Paul Barry, Sep 18 2006
T(n,k) = GegenbauerC(n-k,-n+1,-1/2) + GegenbauerC(n-k-1,-n+1,-1/2). In this form also the missing first column of the triangle 1,1,1,3,7,19,... (cf. A002426) can be computed. - Peter Luschny, May 12 2016
From Peter Bala, Jul 12 2021: (Start)
T(n,k) = Sum_{j = k..n} binomial(n,j)*binomial(j,floor((j-k)/2)).
Matrix product of Riordan arrays ( 1/(1 - x), x/(1 - x) ) * ( (1 - x*c(x^2))/(1 - 2*x), x*c(x^2) ) = A007318 * A061554 (triangle version), where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108.
The n-th row polynomial R(n,x) equals the n-th degree Taylor polynomial of the function (1 + x)*(1 + x + x^2)^n expanded about the point x = 0. - Peter Bala, Sep 06 2022
Extensions
More terms from David W. Wilson
A111418 Right-hand side of odd-numbered rows of Pascal's triangle.
1, 3, 1, 10, 5, 1, 35, 21, 7, 1, 126, 84, 36, 9, 1, 462, 330, 165, 55, 11, 1, 1716, 1287, 715, 286, 78, 13, 1, 6435, 5005, 3003, 1365, 455, 105, 15, 1, 24310, 19448, 12376, 6188, 2380, 680, 136, 17, 1, 92378, 75582, 50388
Offset: 0
Comments
Riordan array (c(x)/sqrt(1-4*x),x*c(x)^2) where c(x) is g.f. of A000108. Unsigned version of A113187. Diagonal sums are A014301(n+1).
Triangle T(n,k),0<=k<=n, read by rows defined by :T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=3*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+2*T(n-1,k)+T(n-1,k+1) for k>=1. - Philippe Deléham, Mar 22 2007
Reversal of A122366. - Philippe Deléham, Mar 22 2007
Column k has e.g.f. exp(2x)(Bessel_I(k,2x)+Bessel_I(k+1,2x)). - Paul Barry, Jun 06 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1 . Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Diagonal sums are A014301(n+1). - Paul Barry, Mar 08 2011
This triangle T(n,k) appears in the expansion of odd powers of Fibonacci numbers F=A000045 in terms of F-numbers with multiples of odd numbers as indices. See the Ozeki reference, p. 108, Lemma 2. The formula is: F_l^(2*n+1) = sum(T(n,k)*(-1)^((n-k)*(l+1))* F_{(2*k+1)*l}, k=0..n)/5^n, n >= 0, l >= 0. - Wolfdieter Lang, Aug 24 2012
Central terms give A052203. - Reinhard Zumkeller, Mar 14 2014
This triangle appears in the expansion of (4*x)^n in terms of the polynomials Todd(n, x):= T(2*n+1, sqrt(x))/sqrt(x) = sum(A084930(n,m)*x^m), n >= 0. This follows from the inversion of the lower triangular Riordan matrix built from A084930 and comparing the g.f. of the row polynomials. - Wolfdieter Lang, Aug 05 2014
From Wolfdieter Lang, Aug 15 2014: (Start)
This triangle is the inverse of the signed Riordan triangle (-1)^(n-m)*A111125(n,m).
This triangle T(n,k) appears in the expansion of x^n in terms of the polynomials todd(k, x):= T(2*k+1, sqrt(x)/2)/(sqrt(x)/2) = S(k, x-2) - S(k-1, x-2) with the row polynomials T and S for the triangles A053120 and A049310, respectively: x^n = sum(T(n,k)*todd(k, x), k=0..n). Compare this with the preceding comment.
The A- and Z-sequences for this Riordan triangle are [1, 2, 1, repeated 0] and [3, 1, repeated 0]. For A- and Z-sequences for Riordan triangles see the W. Lang link under A006232. This corresponds to the recurrences given in the Philippe Deléham, Mar 22 2007 comment above. (End)
Examples
From _Wolfdieter Lang_, Aug 05 2014: (Start) The triangle T(n,k) begins: n\k 0 1 2 3 4 5 6 7 8 9 10 ... 0: 1 1: 3 1 2: 10 5 1 3: 35 21 7 1 4: 126 84 36 9 1 5: 462 330 165 55 11 1 6: 1716 1287 715 286 78 13 1 7: 6435 5005 3003 1365 455 105 15 1 8: 24310 19448 12376 6188 2380 680 136 17 1 9: 92378 75582 50388 27132 11628 3876 969 171 19 1 10: 352716 293930 203490 116280 54264 20349 5985 1330 210 21 1 ... Expansion examples (for the Todd polynomials see A084930 and a comment above): (4*x)^2 = 10*Todd(n, 0) + 5*Todd(n, 1) + 1*Todd(n, 2) = 10*1 + 5*(-3 + 4*x) + 1*(5 - 20*x + 16*x^2). (4*x)^3 = 35*1 + 21*(-3 + 4*x) + 7*(5 - 20*x + 16*x^2) + (-7 + 56*x - 112*x^2 +64*x^3)*1. (End) --------------------------------------------------------------------- Production matrix is 3, 1, 1, 2, 1, 0, 1, 2, 1, 0, 0, 1, 2, 1, 0, 0, 0, 1, 2, 1, 0, 0, 0, 0, 1, 2, 1, 0, 0, 0, 0, 0, 1, 2, 1, 0, 0, 0, 0, 0, 0, 1, 2, 1, 0, 0, 0, 0, 0, 0, 0, 1, 2, 1 - _Paul Barry_, Mar 08 2011 Application to odd powers of Fibonacci numbers F, row n=2: F_l^5 = (10*(-1)^(2*(l+1))*F_l + 5*(-1)^(1*(l+1))*F_{3*l} + 1*F_{5*l})/5^2, l >= 0. - _Wolfdieter Lang_, Aug 24 2012
Links
- Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
- Paul Barry, On the Hurwitz Transform of Sequences, Journal of Integer Sequences, Vol. 15 (2012), #12.8.7.
- E. Deutsch, L. Ferrari and S. Rinaldi, Production Matrices, Advances in Applied Mathematics, 34 (2005) pp. 101-122.
- Asamoah Nkwanta and Earl R. Barnes, Two Catalan-type Riordan Arrays and their Connections to the Chebyshev Polynomials of the First Kind, Journal of Integer Sequences, Article 12.3.3, 2012.
- A. Nkwanta, A. Tefera, Curious Relations and Identities Involving the Catalan Generating Function and Numbers, Journal of Integer Sequences, 16 (2013), #13.9.5.
- K. Ozeki, On Melham's sum, The Fibonacci Quart. 46/47 (2008/2009), no. 2, 107-110.
- Sun, Yidong; Ma, Luping Minors of a class of Riordan arrays related to weighted partial Motzkin paths. Eur. J. Comb. 39, 157-169 (2014), Table 2.2.
Crossrefs
Programs
-
Haskell
a111418 n k = a111418_tabl !! n !! k a111418_row n = a111418_tabl !! n a111418_tabl = map reverse a122366_tabl -- Reinhard Zumkeller, Mar 14 2014
-
Mathematica
Table[Binomial[2*n+1, n-k], {n,0,10}, {k,0,n}] (* G. C. Greubel, May 22 2017 *) T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0, T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]]; Table[T[n, k, 3, 2], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, May 22 2017 *)
Formula
T(n, k) = C(2*n+1, n-k).
Sum_{k=0..n} T(n, k) = 4^n.
Sum_{k, 0<=k<=n}(-1)^k *T(n,k) = binomial(2*n,n) = A000984(n). - Philippe Deléham, Mar 22 2007
T(n,k) = sum{j=k..n, C(n,j)*2^(n-j)*C(j,floor((j-k)/2))}. - Paul Barry, Jun 06 2007
Sum_{k, k>=0} T(m,k)*T(n,k) = T(m+n,0)= A001700(m+n). - Philippe Deléham, Nov 22 2009
G.f. row polynomials: ((1+x) - (1-x)/sqrt(1-4*z))/(2*(x - (1+x)^2*z))
(see the Riordan property mentioned in a comment above). - Wolfdieter Lang, Aug 05 2014
A055248 Triangle of partial row sums of triangle A007318(n,m) (Pascal's triangle). Triangle A008949 read backwards. Riordan (1/(1-2x), x/(1-x)).
1, 2, 1, 4, 3, 1, 8, 7, 4, 1, 16, 15, 11, 5, 1, 32, 31, 26, 16, 6, 1, 64, 63, 57, 42, 22, 7, 1, 128, 127, 120, 99, 64, 29, 8, 1, 256, 255, 247, 219, 163, 93, 37, 9, 1, 512, 511, 502, 466, 382, 256, 130, 46, 10, 1, 1024, 1023, 1013, 968, 848, 638, 386, 176, 56, 11, 1
Offset: 0
Comments
In the language of the Shapiro et al. reference (also given in A053121) such a lower triangular (ordinary) convolution array, considered as matrix, belongs to the Riordan-group. The g.f. for the row polynomials p(n,x) (increasing powers of x) is 1/((1-2*z)*(1-x*z/(1-z))).
Binomial transform of the all 1's triangle: as a Riordan array, it factors to give (1/(1-x),x/(1-x))(1/(1-x),x). Viewed as a number square read by antidiagonals, it has T(n,k) = Sum_{j=0..n} binomial(n+k,n-j) and is then the binomial transform of the Whitney square A004070. - Paul Barry, Feb 03 2005
Riordan array (1/(1-2x), x/(1-x)). Antidiagonal sums are A027934(n+1), n >= 0. - Paul Barry, Jan 30 2005; edited by Wolfdieter Lang, Jan 09 2015
Eigensequence of the triangle = A005493: (1, 3, 10, 37, 151, 674, ...); row sums of triangles A011971 and A159573. - Gary W. Adamson, Apr 16 2009
Read as a square array, this is the generalized Riordan array ( 1/(1 - 2*x), 1/(1 - x) ) as defined in the Bala link (p. 5), which factorizes as ( 1/(1 - x), x/(1 - x) )*( 1/(1 - x), x )*( 1, 1 + x ) = P*U*transpose(P), where P denotes Pascal's triangle, A007318, and U is the lower unit triangular array with 1's on or below the main diagonal. - Peter Bala, Jan 13 2016
Examples
The triangle a(n,m) begins: n\m 0 1 2 3 4 5 6 7 8 9 10 ... 0: 1 1: 2 1 2: 4 3 1 3: 8 7 4 1 4: 16 15 11 5 1 5: 32 31 26 16 6 1 6: 64 63 57 42 22 7 1 7: 128 127 120 99 64 29 8 1 8: 256 255 247 219 163 93 37 9 1 9: 512 511 502 466 382 256 130 46 10 1 10: 1024 1023 1013 968 848 638 386 176 56 11 1 ... Reformatted. - _Wolfdieter Lang_, Jan 09 2015 Fourth row polynomial (n=3): p(3,x)= 8 + 7*x + 4*x^2 + x^3. The matrix inverse starts 1; -2, 1; 2, -3, 1; -2, 5, -4, 1; 2, -7, 9, -5, 1; -2, 9, -16, 14, -6, 1; 2, -11, 25,- 30, 20, -7, 1; -2, 13, -36, 55, -50, 27, -8, 1; 2, -15, 49, -91, 105, -77, 35, -9, 1; -2, 17, -64, 140, -196, 182, -112, 44, -10, 1; 2, -19, 81, -204, 336, -378, 294, -156, 54, -11, 1; ... which may be related to A029653. - _R. J. Mathar_, Mar 29 2013 From _Peter Bala_, Dec 23 2014: (Start) With the array M(k) as defined in the Formula section, the infinite product M(0)*M(1)*M(2)*... begins /1 \ /1 \ /1 \ /1 \ |2 1 ||0 1 ||0 1 | |2 1 | |4 3 1 ||0 2 1 ||0 0 1 |... = |4 5 1 | |8 7 4 1 ||0 4 3 1 ||0 0 2 1 | |8 19 9 1 | |... ||0 8 7 4 1 ||0 0 4 3 1| |... | |... ||... ||... | | | = A143494. (End) Matrix factorization of square array as P*U*transpose(P): /1 \ /1 \ /1 1 1 1 ...\ /1 1 1 1 ...\ |1 1 ||1 1 ||0 1 2 3 ... | |2 3 4 5 ... | |1 2 1 ||1 1 1 ||0 0 1 3 ... | = |4 7 11 16 ... | |1 3 3 1 ||1 1 1 1 ||0 0 0 1 ... | |8 15 26 42 ... | |... ||... ||... | |... | - _Peter Bala_, Jan 13 2016
Links
- Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
- Peter Bala, Notes on generalized Riordan arrays
- Peter Bala, A055248: Rapidly converging series for log(2) and Pi
- Jean-Luc Baril, Javier F. González, and José L. Ramírez, Last symbol distribution in pattern avoiding Catalan words, Univ. Bourgogne (France, 2022).
- Paul Barry, On a Central Transform of Integer Sequences, arXiv:2004.04577 [math.CO], 2020.
- Norman Lindquist and Gerard Sierksma, Extensions of set partitions, Journal of Combinatorial Theory, Series A 31.2 (1981): 190-198. See Table I.
- L. W. Shapiro, S. Getu, Wen-Jin Woan and L. C. Woodson, The Riordan Group, Discrete Appl. Maths. 34 (1991) 229-239.
Crossrefs
Programs
-
Haskell
a055248 n k = a055248_tabl !! n !! k a055248_row n = a055248_tabl !! n a055248_tabl = map reverse a008949_tabl -- Reinhard Zumkeller, Jun 20 2015
-
Maple
T := (n,k) -> 2^n - (1/2)*binomial(n, k-1)*hypergeom([1, n + 1], [n-k + 2], 1/2). seq(seq(simplify(T(n,k)), k=0..n),n=0..10); # Peter Luschny, Oct 10 2019
-
Mathematica
a[n_, m_] := Sum[ Binomial[n, m + j], {j, 0, n}]; Table[a[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (* Jean-François Alcover, Jul 05 2013, after Paul Barry *) T[n_, k_] := Binomial[n, k] * Hypergeometric2F1[1, k - n, k + 1, -1]; Flatten[Table[T[n, k], {n, 0, 7}, {k, 0, n}]] (* Peter Luschny, Oct 06 2023 *)
Formula
a(n, m) = A008949(n, n-m), if n > m >= 0.
a(n, m) = Sum_{k=m..n} A007318(n, k) (partial row sums in columns m).
Column m recursion: a(n, m) = Sum_{j=m..n-1} a(j, m) + A007318(n, m) if n >= m >= 0, a(n, m) := 0 if n
G.f. for column m: (1/(1-2*x))*(x/(1-x))^m, m >= 0.
a(n, m) = Sum_{j=0..n} binomial(n, m+j). - Paul Barry, Feb 03 2005
Inverse binomial transform (by columns) of A112626. - Ross La Haye, Dec 31 2006
T(2n,n) = A032443(n). - Philippe Deléham, Sep 16 2009
From Peter Bala, Dec 23 2014: (Start)
Exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(8 + 7*x + 4*x^2/2! + x^3/3!) = 8 + 15*x + 26*x^2/2! + 42*x^3/3! + 64*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ).
Let M denote the present triangle. For k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0 M/ having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. The infinite product M(0)*M(1)*M(2)*..., which is clearly well-defined, is equal to A143494 (but with a different offset). See the Example section. Cf. A106516. (End)
a(n,m) = Sum_{p=m..n} 2^(n-p)*binomial(p-1,m-1), n >= m >= 0, else 0. - Wolfdieter Lang, Jan 09 2015
T(n, k) = 2^n - (1/2)*binomial(n, k-1)*hypergeom([1, n+1], [n-k+2], 1/2). - Peter Luschny, Oct 10 2019
T(n, k) = binomial(n, k)*hypergeom([1, k - n], [k + 1], -1). - Peter Luschny, Oct 06 2023
n-th row polynomial R(n, x) = (2^n - x*(1 + x)^n)/(1 - x). These polynomials can be used to find series acceleration formulas for the constants log(2) and Pi. - Peter Bala, Mar 03 2025
A037952 a(n) = binomial(n, floor((n-1)/2)).
0, 1, 1, 3, 4, 10, 15, 35, 56, 126, 210, 462, 792, 1716, 3003, 6435, 11440, 24310, 43758, 92378, 167960, 352716, 646646, 1352078, 2496144, 5200300, 9657700, 20058300, 37442160, 77558760, 145422675, 300540195, 565722720, 1166803110, 2203961430, 4537567650
Offset: 0
Keywords
Comments
The maximum size of an intersecting (or proper) antichain on an n-set. - Vladeta Jovovic, Dec 27 2000
Number of ordered trees with n+1 edges, having root of degree at least 2 and nonroot nodes of outdegree 0 or 2. - Emeric Deutsch, Aug 02 2002
a(n)=number of Dyck (n+1)-paths that are symmetric but not prime. A prime Dyck path is one that returns to the x-axis only at its terminal point. For example a(3)=3 counts UDUUDDUD, UUDDUUDD, UDUDUDUD. - David Callan, Dec 09 2004
Number of involutions of [n+2] containing the pattern 132 exactly once. For example, a(3)=3 because we have 1'3'2'45, 42'5'13' and 52'4'3'1 (the entries corresponding to the pattern 132 are "primed"). - Emeric Deutsch, Nov 17 2005
Also number of ways to put n eggs in floor(n/2) baskets where order of the baskets matters and all baskets have at least 1 egg. - Ben Paul Thurston, Sep 30 2006
For n >= 1 the number of standard Young tableaux with shapes corresponding to partitions into at most 2 distinct parts. - Joerg Arndt, Oct 25 2012
It seems that 3, 4, 10, ... are Colbourn's Covering Array Numbers CAN(2,k,2). - Ryan Dougherty, May 27 2015
Essentially the same as A007007. - Georg Fischer, Oct 02 2018
a(n) is the number of subsets of {1,2,...,n} that contain exactly 1 more odd than even elements. For example, for n = 6, a(6) = 15 and the 15 sets are {1}, {3}, {5}, {1,2,3}, {1,2,5}, {1,3,4}, {1,3,6}, {1,4,5}, {1,5,6}, {2,3,5}, {3,4,5}, {3,5,6}, {1,2,3,4,5}, {1,2,3,5,6}, {1,3,4,5,6}. - Enrique Navarrete, Dec 21 2019
a(n) is the number of lattice paths of n steps taken from the step set {U=(1,1), D=(1,-1)} that start at the origin, never go below the x-axis, and end strictly above the x-axis; more succinctly, proper left factors of Dyck paths. For example, a(3)=3 counts UUU, UUD, UDU, and a(4)=4 counts UUUU, UUUD, UUDU, UDUU. - David Callan and Emeric Deutsch, Jan 25 2021
For n >= 3, a(n) is also the number of pinnacle sets in the (n-2)-Plummer-Toft graph. - Eric W. Weisstein, Sep 11 2024
Links
- Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
- Cyril Banderier and Michael Wallner, Lattice paths with catastrophes, arXiv:1707.01931 [math.CO], 2017.
- J.-L. Baril, Classical sequences revisited with permutations avoiding dotted pattern, Electronic Journal of Combinatorics, 18 (2011), #P178.
- Jean-Luc Baril, Sergey Kirgizov, and Armen Petrossian, Dyck Paths with catastrophes modulo the positions of a given pattern, Australasian J. Comb. (2022) Vol. 84, No. 2, 398-418.
- Jean-Luc Baril and José L. Ramírez, Fibonacci and Catalan paths in a wall, 2023.
- Paul Barry, d-orthogonal polynomials, Fuss-Catalan matrices and lattice paths, arXiv:2505.16718 [math.CO], 2025. See p. 25.
- C. J. Colbourn, Table of CAN(2, k, 2)
- Emeric Deutsch, Ordered trees with prescribed root degrees, node degrees and branch lengths, Discrete Math., 282, 2004, 89-94.
- O. Guibert and T. Mansour, Restricted 132-involutions, Sem. Lotharingien de Combinatoire, 48, 2002, Article B48a (Corollary 4.2).
- M. Miyakawa, A. Nozaki, G. Pogosyan, and I. G. Rosenberg, A map from the lower-half of the n-Cube onto the (n-1)-Cube which preserves intersecting antichains, Discr. Appl. Math. 92 (2-3) (1999) 223-228.
- M. van de Vel, Determination of msd(L^n), J. Algebraic Combin., 9 (1999), 161-171.
- Eric Weisstein's World of Mathematics, Pinnacle Set.
- Eric Weisstein's World of Mathematics, Plummer-Toft Graph.
Crossrefs
Programs
-
Haskell
a037952 n = a037952_list !! n a037952_list = zipWith (-) (tail a001405_list) a001405_list -- Reinhard Zumkeller, Mar 04 2012
-
Magma
[Binomial(n, Floor((n-1)/2)): n in [0..40]]; // G. C. Greubel, Jun 21 2022
-
Maple
a:= n-> binomial(n, floor((n-1)/2)): seq(a(n), n=0..35); # Alois P. Heinz, Sep 19 2017
-
Mathematica
Table[ Binomial[n, Floor[n/2]], {n, 0, 35}]//Differences (* Jean-François Alcover, Jun 10 2013 *) f[n_] := Binomial[n, Floor[(n-1)/2]]; Array[f, 35, 0] (* Robert G. Wilson v, Nov 13 2014 *)
-
PARI
a(n) = binomial(n, (n-1)\2); \\ Altug Alkan, Oct 03 2018
-
SageMath
[binomial(n, (n-1)//2) for n in (0..40)] # G. C. Greubel, Jun 21 2022
Formula
E.g.f.: BesselI(1, 2*x) + BesselI(2, 2*x). - Vladeta Jovovic, Apr 28 2003
O.g.f.: (1-sqrt(1-4x^2))/(x - 2x^2 + x*sqrt(1-4x^2)).
Convolution of A001405 and A126120 shifted right: g001405(x)*g126120(x) = g037952(x)/x. - Philippe Deléham, Mar 17 2007
D-finite with recurrence: (n+2)*a(n) + (-n-2)*a(n-1) + 2*(-2*n+1)*a(n-2) + 4*(n-2)*a(n-3) = 0. - R. J. Mathar, Jan 25 2013. Proved by Robert Israel, Nov 13 2014
For n > 0: a(n) = A265848(n,0). - Reinhard Zumkeller, Dec 24 2015
a(n) = binomial(n, (n-2)/2) = A001791(n/2), n even; a(n) = binomial(n, (n+1)/2) = A001700((n-1)/2), n odd. - Enrique Navarrete, Dec 21 2019
From R. J. Mathar, Sep 23 2021: (Start)
a(n) = Sum_{m=1..n} A053121(n,m) [comment Callan-Deutsch].
a(2n+1) = A000984(n+1)/2. (End)
a(n) = Sum_{k=2..n} A143359(n,k). [Callan's 2004 comment]. - R. J. Mathar, Sep 24 2021
From Amiram Eldar, Sep 27 2024: (Start)
Sum_{n>=1} 1/a(n) = 1 + Pi/sqrt(3).
Sum_{n>=1} (-1)^(n+1)/a(n) = (3 - Pi/sqrt(3))/9. (End)
A059365 Another version of the Catalan triangle: T(r,s) = binomial(2*r-s-1,r-1) - binomial(2*r-s-1,r), r>=0, 0 <= s <= r.
0, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 5, 5, 3, 1, 0, 14, 14, 9, 4, 1, 0, 42, 42, 28, 14, 5, 1, 0, 132, 132, 90, 48, 20, 6, 1, 0, 429, 429, 297, 165, 75, 27, 7, 1, 0, 1430, 1430, 1001, 572, 275, 110, 35, 8, 1, 0, 4862, 4862, 3432, 2002, 1001, 429, 154, 44
Offset: 0
Examples
Triangle starts 0; 0, 1; 0, 1, 1; 0, 2, 2, 1; 0, 5, 5, 3, 1; 0, 14, 14, 9, 4, 1; 0, 42, 42, 28, 14, 5, 1; 0, 132, 132, 90, 48, 20, 6, 1; 0, 429, 429, 297, 165, 75, 27, 7, 1; 0, 1430, 1430, 1001, 572, 275, 110, 35, 8, 1; 0, 4862, 4862, 3432, 2002, 1001, 429, 154, 44, 9, 1; ...
Links
- G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
- F. R. Bernhart, Catalan, Motzkin and Riordan numbers, Discr. Math., 204 (1999), 73-112.
- D. Callan, A recursive bijective approach to counting permutations containing 3-letter patterns, arXiv:math/0211380 [math.CO], 2002.
- FindStat - Combinatorial Statistic Finder, The number of touch points of a Dyck path., The number of initial rises of a Dyck paths., The number of nodes on the left branch of the tree., The number of subtrees.
- A. Robertson, D. Saracino and D. Zeilberger, Refined restricted permutations, arXiv:math/0203033 [math.CO], 2002.
- Index to sequences related to Catalan
Crossrefs
The three triangles A059365, A106566 and A099039 are the same except for signs and the leading term.
Essentially the same as A033184.
Programs
-
Magma
/* as triangle */ [[[0] cat [Binomial(2*r-s-1, r-1)- Binomial(2*r-s-1, r): s in [1..r]]: r in [0..10]]]; // Vincenzo Librandi, Jan 09 2017
-
Mathematica
Table[Binomial[2*r - s - 1, r - 1] - Binomial[2*r - s - 1, r], {r, 0, 10}, {s, 0, r}] // Flatten (* G. C. Greubel, Jan 08 2017 *)
-
PARI
tabl(nn) = { print(0); for (r=1, nn, for (s=0, r, print1(binomial(2*r-s-1,r-1)-binomial(2*r-s-1,r), ", ");); print(););} \\ Michel Marcus, Nov 01 2013
Comments