[go: up one dir, main page]

login
Search: a000910 -id:a000910
     Sort: relevance | references | number | modified | created      Format: long | short | data
Triangle read by rows: T(n,k) = C(n+k,n)*C(n,k)/(k+1), for n >= 0, k = 0..n.
+10
38
1, 1, 1, 1, 3, 2, 1, 6, 10, 5, 1, 10, 30, 35, 14, 1, 15, 70, 140, 126, 42, 1, 21, 140, 420, 630, 462, 132, 1, 28, 252, 1050, 2310, 2772, 1716, 429, 1, 36, 420, 2310, 6930, 12012, 12012, 6435, 1430, 1, 45, 660, 4620, 18018, 42042, 60060, 51480, 24310, 4862
OFFSET
0,5
COMMENTS
Row sums: A006318 (Schroeder numbers). Essentially same as triangle A060693 transposed.
T(n,k) is number of Schroeder paths (i.e., consisting of steps U=(1,1), D=(1,-1), H=(2,0) and never going below the x-axis) from (0,0) to (2n,0), having k U's. E.g., T(2,1)=3 because we have UHD, UDH and HUD. - Emeric Deutsch, Dec 06 2003
Little Schroeder numbers A001003 have a(n) = Sum_{k=0..n} A088617(n,k)*(-1)^(n-k)*2^k. - Paul Barry, May 24 2005
Conjecture: The expected number of U's in a Schroeder n-path is asymptotically Sqrt[1/2]*n for large n. - David Callan, Jul 25 2008
T(n, k) is also the number of order-preserving and order-decreasing partial transformations (of an n-chain) of width k (width(alpha) = |Dom(alpha)|). - Abdullahi Umar, Oct 02 2008
The antidiagonals of this lower triangular matrix are the rows of A055151. - Tom Copeland, Jun 17 2015
REFERENCES
Charles Jordan, Calculus of Finite Differences, Chelsea 1965, p. 449.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 0..11475 (rows 0 <= n <= 150)
Anwar Al Ghabra, K. Gopala Krishna, Patrick Labelle, and Vasilisa Shramchenko, Enumeration of multi-rooted plane trees, arXiv:2301.09765 [math.CO], 2023.
Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.4.
Paul Barry, Generalized Catalan Numbers Associated with a Family of Pascal-like Triangles, J. Int. Seq., Vol. 22 (2019), Article 19.5.8.
Paul Barry, On the inversion of Riordan arrays, arXiv:2101.06713 [math.CO], 2021.
Manosij Ghosh Dastidar and Michael Wallner, Bijections and congruences involving lattice paths and integer compositions, arXiv:2402.17849 [math.CO], 2024. See p. 16.
Samuele Giraudo, Tree series and pattern avoidance in syntax trees, arXiv:1903.00677 [math.CO], 2019.
Hsien-Kuei Hwang and Satoshi Kuriki, Integrated empirical measures and generalizations of classical goodness-of-fit statistics, arXiv:2404.06040 [math.ST], 2024. See p. 11.
C. Jordan, Calculus of Finite Differences, Budapest, 1939. [Annotated scans of pages 448-450 only]
M. Klazar, On numbers of Davenport-Schinzel sequences, Discr. Math., 185 (1998), 77-87.
Paul W. Lapey and Aaron Williams, A Shift Gray Code for Fixed-Content Łukasiewicz Words, Williams College, 2022.
A. Laradji and A. Umar, A. Combinatorial results for semigroups of order-preserving partial transformations, Journal of Algebra 278, (2004), 342-359.
A. Laradji and A. Umar, Combinatorial results for semigroups of order-decreasing partial transformations, J. Integer Seq. 7 (2004), 04.3.8.
Jason P. Smith, The poset of graphs ordered by induced containment, arXiv:1806.01821 [math.CO], 2018.
FORMULA
Triangle T(n, k) read by rows; given by [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...] DELTA [[1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...] where DELTA is Deléham's operator defined in A084938.
T(n, k) = A085478(n, k)*A000108(k); A000108 = Catalan numbers. - Philippe Deléham, Dec 05 2003
Sum_{k=0..n} T(n, k)*x^k*(1-x)^(n-k) = A000108(n), A001003(n), A007564(n), A059231(n), A078009(n), A078018(n), A081178(n), A082147(n), A082181(n), A082148(n), A082173(n) for x = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11. - Philippe Deléham, Aug 18 2005
Sum_{k=0..n} T(n,k)*x^k = (-1)^n*A107841(n), A080243(n), A000007(n), A000012(n), A006318(n), A103210(n), A103211(n), A133305(n), A133306(n), A133307(n), A133308(n), A133309(n) for x = -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively. - Philippe Deléham, Oct 18 2007
O.g.f. (with initial 1 excluded) is the series reversion with respect to x of (1-t*x)*x/(1+x). Cf. A062991 and A089434. - Peter Bala, Jul 31 2012
G.f.: 1 + (1 - x - T(0))/y, where T(k) = 1 - x*(1+y)/( 1 - x*y/T(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Nov 03 2013
From Peter Bala, Jul 20 2015: (Start)
O.g.f. A(x,t) = ( 1 - x - sqrt((1 - x)^2 - 4*x*t) )/(2*x*t) = 1 + (1 + t)*x + (1 + 3*t + 2*t^2)*x^2 + ....
1 + x*(dA(x,t)/dx)/A(x,t) = 1 + (1 + t)*x + (1 + 4*t + 3*t^2)*x^2 + ... is the o.g.f. for A123160.
For n >= 1, the n-th row polynomial equals (1 + t)/(n+1)*Jacobi_P(n-1,1,1,2*t+1). Removing a factor of 1 + t from the row polynomials gives the row polynomials of A033282. (End)
From Tom Copeland, Jan 22 2016: (Start)
The o.g.f. G(x,t) = {1 - (2t+1) x - sqrt[1 - (2t+1) 2x + x^2]}/2x = (t + t^2) x + (t + 3t^2 + 2t^3) x^2 + (t + 6t^2 + 10t^3 + 5t^3) x^3 + ... generating shifted rows of this entry, excluding the first, was given in my 2008 formulas for A033282 with an o.g.f. f1(x,t) = G(x,t)/(1+t) for A033282. Simple transformations presented there of f1(x,t) are related to A060693 and A001263, the Narayana numbers. See also A086810.
The inverse of G(x,t) is essentially given in A033282 by x1, the inverse of f1(x,t): Ginv(x,t) = x [1/(t+x) - 1/(1+t+x)] = [((1+t) - t) / (t(1+t))] x - [((1+t)^2 - t^2) / (t(1+t))^2] x^2 + [((1+t)^3 - t^3) / (t(1+t))^3] x^3 - ... . The coefficients in t of Ginv(xt,t) are the o.g.f.s of the diagonals of the Pascal triangle A007318 with signed rows and an extra initial column of ones. The numerators give the row o.g.f.s of signed A074909.
Rows of A088617 are shifted columns of A107131, whose reversed rows are the Motzkin polynomials of A055151, related to A011973. The diagonals of A055151 give the rows of A088671, and the antidiagonals (top to bottom) of A088617 give the rows of A107131 and reversed rows of A055151. The diagonals of A107131 give the columns of A055151. The antidiagonals of A088617 (bottom to top) give the rows of A055151.
(End)
T(n, k) = [x^k] hypergeom([-n, 1 + n], [2], -x). - Peter Luschny, Apr 26 2022
EXAMPLE
Triangle begins:
[0] 1;
[1] 1, 1;
[2] 1, 3, 2;
[3] 1, 6, 10, 5;
[4] 1, 10, 30, 35, 14;
[5] 1, 15, 70, 140, 126, 42;
[6] 1, 21, 140, 420, 630, 462, 132;
[7] 1, 28, 252, 1050, 2310, 2772, 1716, 429;
[8] 1, 36, 420, 2310, 6930, 12012, 12012, 6435, 1430;
[9] 1, 45, 660, 4620, 18018, 42042, 60060, 51480, 24310, 4862;
MAPLE
R := n -> simplify(hypergeom([-n, n + 1], [2], -x)):
Trow := n -> seq(coeff(R(n, x), x, k), k = 0..n):
seq(print(Trow(n)), n = 0..9); # Peter Luschny, Apr 26 2022
MATHEMATICA
Table[Binomial[n+k, n] Binomial[n, k]/(k+1), {n, 0, 10}, {k, 0, n}]//Flatten (* Michael De Vlieger, Aug 10 2017 *)
PROG
(PARI) {T(n, k)= if(k+1, binomial(n+k, n)*binomial(n, k)/(k+1))}
(Magma) [[Binomial(n+k, n)*Binomial(n, k)/(k+1): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Jun 18 2015
(SageMath) flatten([[binomial(n+k, 2*k)*catalan_number(k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 22 2022
KEYWORD
nonn,tabl,easy
AUTHOR
N. J. A. Sloane, Nov 23 2003
STATUS
approved
Triangle (0 <= k <= n) read by rows: T(n, k) is the number of Schröder paths from (0,0) to (2n,0) having k peaks.
+10
25
1, 1, 1, 2, 3, 1, 5, 10, 6, 1, 14, 35, 30, 10, 1, 42, 126, 140, 70, 15, 1, 132, 462, 630, 420, 140, 21, 1, 429, 1716, 2772, 2310, 1050, 252, 28, 1, 1430, 6435, 12012, 12012, 6930, 2310, 420, 36, 1, 4862, 24310, 51480, 60060, 42042, 18018, 4620, 660, 45, 1, 16796
OFFSET
0,4
COMMENTS
The rows sum to A006318 (Schroeder numbers), the left column is A000108 (Catalan numbers); the next-to-left column is A001700, the alternating sum in each row but the first is 0.
T(n,k) is the number of Schroeder paths (i.e., consisting of steps U=(1,1), D=(1,-1), H=(2,0) and never going below the x-axis) from (0,0) to (2n,0), having k peaks. Example: T(2,1)=3 because we have UU*DD, U*DH and HU*D, the peaks being shown by *. E.g., T(n,k) = binomial(n,k)*binomial(2n-k,n-1)/n for n>0. - Emeric Deutsch, Dec 06 2003
A090181*A007318 as infinite lower triangular matrices. - Philippe Deléham, Oct 14 2008
T(n,k) is also the number of rooted plane trees with maximal degree 3 and k vertices of degree 2 (a node may have at most 2 children, and there are exactly k nodes with 1 child). Equivalently, T(n,k) is the number of syntactically different expressions that can be formed that use a unary operation k times, a binary operation n-k times, and nothing else (sequence of operands is fixed). - Lars Hellstrom (Lars.Hellstrom(AT)residenset.net), Dec 08 2009
LINKS
Vincenzo Librandi, Rows n = 0..100, flattened
J. Agapito, A. Mestre, P. Petrullo, and M. Torres, Counting noncrossing partitions via Catalan triangles, CEAFEL Seminar, June 30, 2015
Jean-Christophe Aval and François Bergeron, Rectangular Schröder Parking Functions Combinatorics, arXiv:1603.09487 [math.CO], 2016.
Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, J. Integer Sequ., Vol. 9 (2006), Article 06.2.4.
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
Paul Barry, Riordan Pseudo-Involutions, Continued Fractions and Somos 4 Sequences, arXiv:1807.05794 [math.CO], 2018.
Paul Barry, The Central Coefficients of a Family of Pascal-like Triangles and Colored Lattice Paths, J. Int. Seq., Vol. 22 (2019), Article 19.1.3.
Paul Barry, Generalized Catalan Numbers Associated with a Family of Pascal-like Triangles, J. Int. Seq., Vol. 22 (2019), Article 19.5.8.
Paul Barry, On the inversion of Riordan arrays, arXiv:2101.06713 [math.CO], 2021.
Paul Barry, On Motzkin-Schröder Paths, Riordan Arrays, and Somos-4 Sequences, J. Int. Seq. (2023) Vol. 26, Art. 23.4.7.
David Callan and Toufik Mansour, Five subsets of permutations enumerated as weak sorting permutations, arXiv:1602.05182 [math.CO], 2016.
Samuele Giraudo, Tree series and pattern avoidance in syntax trees, arXiv:1903.00677 [math.CO], 2019.
Nate Kube and Frank Ruskey, Sequences That Satisfy a(n-a(n))=0, Journal of Integer Sequences, Vol. 8 (2005), Article 05.5.5.
Krishna Menon and Anurag Singh, Grassmannian permutations avoiding identity, arXiv:2212.13794 [math.CO], 2022.
Jean-Christophe Novelli and Jean-Yves Thibon, Duplicial algebras and Lagrange inversion, arXiv preprint arXiv:1209.5959 [math.CO], 2012.
FORMULA
Triangle T(n, k) (0 <= k <= n) read by rows; given by [1, 1, 1, 1, 1, ...] DELTA [1, 0, 1, 0, 1, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Aug 12 2003
If C_n(x) is the g.f. of row n of the Narayana numbers (A001263), C_n(x) = Sum_{k=1..n} binomial(n,k-1)*(binomial(n-1,k-1)/k) * x^k and T_n(x) is the g.f. of row n of T(n,k), then T_n(x) = C_n(x+1), or T(n,k) = [x^n]Sum_{k=1..n}(A001263(n,k)*(x+1)^k). - Mitch Harris, Jan 16 2007, Jan 31 2007
G.f.: (1 - t*y - sqrt((1-y*t)^2 - 4*y)) / 2.
T(n, k) = binomial(2n-k, n)*binomial(n, k)/(n-k+1). - Philippe Deléham, Dec 07 2003
A060693(n, k) = binomial(2*n-k, k)*A000108(n-k); A000108: Catalan numbers. - Philippe Deléham, Dec 30 2003
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000108(n), A006318(n), A047891(n+1), A082298(n), A082301(n), A082302(n), A082305(n), A082366(n), A082367(n), for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, respectively. - Philippe Deléham, Apr 01 2007
T(n,k) = Sum_{j>=0} A090181(n,j)*binomial(j,k). - Philippe Deléham, May 04 2007
Sum_{k=0..n} T(n,k)*x^(n-k) = (-1)^n*A107841(n), A080243(n), A000007(n), A000012(n), A006318(n), A103210(n), A103211(n), A133305(n), A133306(n), A133307(n), A133308(n), A133309(n) for x = -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, respectively. - Philippe Deléham, Oct 18 2007
From Paul Barry, Jan 29 2009: (Start)
G.f.: 1/(1-xy-x/(1-xy-x/(1-xy-x/(1-xy-x/(1-xy-x/(1-.... (continued fraction);
G.f.: 1/(1-(x+xy)/(1-x/(1-(x+xy)/(1-x/(1-(x+xy)/(1-.... (continued fraction). (End)
T(n,k) = [k<=n]*(Sum_{j=0..n} binomial(n,j)^2*binomial(j,k))/(n-k+1). - Paul Barry, May 28 2009
T(n,k) = A104684(n,k)/(n-k+1). - Peter Luschny, May 17 2011
From Tom Copeland, Sep 21 2011: (Start)
With F(x,t) = (1-(2+t)*x-sqrt(1-2*(2+t)*x+(t*x)^2))/(2*x) an o.g.f. (nulling the n=0 term) in x for the A060693 polynomials in t,
G(x,t) = x/(1+t+(2+t)*x+x^2) is the compositional inverse in x.
Consequently, with H(x,t) = 1/(dG(x,t)/dx) = (1+t+(2+t)*x+x^2)^2 / (1+t-x^2), the n-th A060693 polynomial in t is given by (1/n!)*((H(x,t)*d/dx)^n) x evaluated at x=0, i.e., F(x,t) = exp(x*H(u,t)*d/d) u, evaluated at u = 0.
Also, dF(x,t)/dx = H(F(x,t),t). (End)
See my 2008 formulas in A033282 to relate this entry to A088617, A001263, A086810, and other matrices. - Tom Copeland, Jan 22 2016
Rows of this entry are non-vanishing antidiagonals of A097610. See p. 14 of Agapito et al. for a bivariate generating function and its inverse. - Tom Copeland, Feb 03 2016
From Werner Schulte, Jan 09 2017: (Start)
T(n,k) = A126216(n,k-1) + A126216(n,k) for 0 < k < n;
Sum_{k=0..n} (-1)^k*(1+x*(n-k))*T(n,k) = x + (1-x)*A000007(n).
(End)
Conjecture: Sum_{k=0..n} (-1)^k*T(n,k)*(n+1-k)^2 = 1+n+n^2. - Werner Schulte, Jan 11 2017
EXAMPLE
Triangle begins:
00: [ 1]
01: [ 1, 1]
02: [ 2, 3, 1]
03: [ 5, 10, 6, 1]
04: [ 14, 35, 30, 10, 1]
05: [ 42, 126, 140, 70, 15, 1]
06: [ 132, 462, 630, 420, 140, 21, 1]
07: [ 429, 1716, 2772, 2310, 1050, 252, 28, 1]
08: [ 1430, 6435, 12012, 12012, 6930, 2310, 420, 36, 1]
09: [ 4862, 24310, 51480, 60060, 42042, 18018, 4620, 660, 45, 1]
10: [16796, 92378, 218790, 291720, 240240, 126126, 42042, 8580, 990, 55, 1]
...
MAPLE
A060693 := (n, k) -> binomial(n, k)*binomial(2*n-k, n)/(n-k+1); # Peter Luschny, May 17 2011
MATHEMATICA
t[n_, k_] := Binomial[n, k]*Binomial[2 n - k, n]/(n - k + 1); Flatten[Table[t[n, k], {n, 0, 9}, {k, 0, n}]] (* Robert G. Wilson v, May 30 2011 *)
PROG
(PARI) T(n, k) = binomial(n, k)*binomial(2*n - k, n)/(n - k + 1);
for(n=0, 10, for(k=0, n, print1(T(n, k), ", ")); print); \\ Indranil Ghosh, Jul 28 2017
(Python)
from sympy import binomial
def T(n, k): return binomial(n, k) * binomial(2 * n - k, n) / (n - k + 1)
for n in range(11): print([T(n, k) for k in range(n + 1)]) # Indranil Ghosh, Jul 28 2017
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
F. Chapoton, Apr 20 2001
EXTENSIONS
More terms from Vladeta Jovovic, Apr 21 2001
New description from Philippe Deléham, Aug 12 2003
New name using a comment by Emeric Deutsch from Peter Luschny, Jul 26 2017
STATUS
approved
a(n) = (n-3)*(n-2)*(n-1)*n*(n+1)/30.
+10
6
0, 0, 0, 0, 4, 24, 84, 224, 504, 1008, 1848, 3168, 5148, 8008, 12012, 17472, 24752, 34272, 46512, 62016, 81396, 105336, 134596, 170016, 212520, 263120, 322920, 393120, 475020, 570024, 679644, 805504, 949344, 1113024, 1298528, 1507968, 1743588, 2007768, 2303028
OFFSET
0,5
COMMENTS
The following sequences are provided by the formula n*binomial(n,k) - binomial(n,k+1) = k*binomial(n+1,k+1):
. A000217(n) for k=1,
. A007290(n+1) for k=2,
. A050534(n) for k=3,
. a(n) for k=4,
. A000910(n+1) for k=5.
Sum of reciprocals of a(n), for n>3: 5/16.
From a(2), convolution of oblong numbers (A002378) with themselves. - Bruno Berselli, Oct 24 2016
LINKS
C. P. Neuman and D. I. Schonbach, Evaluation of sums of convolved powers using Bernoulli numbers, SIAM Rev. 19 (1977), no. 1, 90--99. MR0428678 (55 #1698). See Table 3. - N. J. A. Sloane, Mar 23 2014
FORMULA
G.f.: 4*x^4/(1-x)^6.
a(n) = n*binomial(n,4)-binomial(n,5) = 4*binomial(n+1,5) = 4*A000389(n+1).
a(n) = 2*A177747(2*n-7), with A177747(-7) = A177747(-5) = A177747(-3) = A177747(-1) = 0.
(n-4)*a(n) = (n+1)*a(n-1).
E.g.f.: (1/30)*x^4*(5+x)*exp(x). - G. C. Greubel, May 23 2022
Sum_{n>=4} (-1)^n/a(n) = 20*log(2) - 655/48. - Amiram Eldar, Jun 02 2022
MAPLE
f:=n->(n^5-5*n^3+4*n)/30;
[seq(f(n), n=0..50)]; # N. J. A. Sloane, Mar 23 2014
MATHEMATICA
LinearRecurrence[{6, -15, 20, -15, 6, -1}, {0, 0, 0, 0, 4, 24}, 39]
CoefficientList[Series[4x^4/(1-x)^6, {x, 0, 40}], x] (* Vincenzo Librandi, Mar 24 2014 *)
Times@@@Partition[Range[-3, 40], 5, 1]/30 (* Harvey P. Dale, Sep 19 2020 *)
PROG
(Magma) [4*Binomial(n+1, 5): n in [0..38]];
(Maxima) makelist(coeff(taylor(4*x^4/(1-x)^6, x, 0, n), x, n), n, 0, 38);
(PARI) a(n)=(n-3)*(n-2)*(n-1)*n*(n+1)/30 \\ Charles R Greathouse IV, Oct 07 2015
(SageMath) [4*binomial(n+1, 5) for n in (0..40)] # G. C. Greubel, May 23 2022
CROSSREFS
First differences are in A033488.
KEYWORD
nonn,easy
AUTHOR
Bruno Berselli, Mar 23 2012
STATUS
approved
Triangular array of ways of drawing k non-intersecting chords between n points on a circle; i.e., Motzkin polynomial coefficients.
+10
4
1, 1, 0, 1, 1, 0, 1, 3, 0, 0, 1, 6, 2, 0, 0, 1, 10, 10, 0, 0, 0, 1, 15, 30, 5, 0, 0, 0, 1, 21, 70, 35, 0, 0, 0, 0, 1, 28, 140, 140, 14, 0, 0, 0, 0, 1, 36, 252, 420, 126, 0, 0, 0, 0, 0, 1, 45, 420, 1050, 630, 42, 0, 0, 0, 0, 0, 1, 55, 660, 2310, 2310, 462, 0, 0, 0, 0, 0, 0, 1, 66, 990, 4620
OFFSET
0,8
FORMULA
For n >= 2k: T(n, k) = n!/((n-2k)!k!(k+1)!) = A007318(n, 2k)*A000108(k).
T(n,k) = A055151(n,k).
EXAMPLE
Rows start: 1; 1,0; 1,1,0; 1,3,0,0; 1,6,2,0,0; 1,10,10,0,0,0; 1,15,30,5,0,0,0; etc.
CROSSREFS
Visible version of A055151. Row sums are A001006 (Motzkin numbers). Columns include A000012, A000217, A034827 and perhaps A000910.
KEYWORD
nonn,tabl
AUTHOR
Henry Bottomley, Jan 31 2003
STATUS
approved
Triangle read by rows. The Motzkin triangle, the coefficients of the Motzkin polynomials. M(n, k) = binomial(n, k) * CatalanNumber(k/2) if k is even, otherwise 0.
+10
3
1, 1, 0, 1, 0, 1, 1, 0, 3, 0, 1, 0, 6, 0, 2, 1, 0, 10, 0, 10, 0, 1, 0, 15, 0, 30, 0, 5, 1, 0, 21, 0, 70, 0, 35, 0, 1, 0, 28, 0, 140, 0, 140, 0, 14, 1, 0, 36, 0, 252, 0, 420, 0, 126, 0, 1, 0, 45, 0, 420, 0, 1050, 0, 630, 0, 42, 1, 0, 55, 0, 660, 0, 2310, 0, 2310, 0, 462, 0
OFFSET
0,9
COMMENTS
The generalized Motzkin numbers M(n, k) are a refinement of the Motzkin numbers M(n) (A001006) in the sense that they are coefficients of polynomials M(n, x) = Sum_{n..k} M(n, k) * x^k that take the value M(n) at x = 1. The coefficients of x^n are the aerated Catalan numbers A126120.
Variants are the irregular triangle A055151 with zeros deleted, A097610 with reversed rows, A107131 and A080159.
In the literature the name 'Motzkin triangle' is also used for the triangle A026300, which is generated from the powers of the generating function of the Motzkin numbers.
LINKS
M. Artioli, G. Dattoli, S. Licciardi, and S. Pagnutti, Motzkin Numbers: an Operational Point of View, arXiv:1703.07262 [math.CO], 2017, (Table 1 on p. 3).
FORMULA
Let p(n, x) = hypergeom([(1 - n)/2, -n/2], [2], (2*x)^2).
p(n, 1) = A001006(n); p(n, sqrt(2)) = A025235(n); p(n, 2) = A091147(n).
p(2, n) = A002522(n); p(3, n) = A056107(n).
p(n, n) = A359649(n); 2^n*p(n, 1/2) = A000108(n+1).
M(n, k) = [x^k] p(n, x).
M(n, k) = [t^k] (n! * [x^n] exp(x) * BesselI(1, 2*t*x) / (t*x)).
M(n, k) = [t^k][x^n] ((1 - x - sqrt((x-1)^2 - (2*t*x)^2)) / (2*(t*x)^2)).
M(n, n) = A126120(n).
M(n, n-1) = A138364(n), the number of Motzkin n-paths with exactly one flat step.
M(2*n, 2*n) = A000108(n), the number of peakless Motzkin paths having a total of n up and level steps.
M(4*n, 2*n) = A359647(n), the central terms without zeros.
M(2*n+2, 2*n) = A002457(n) = (-4)^n * binomial(-3/2, n).
Sum_{k=0..n} M(n - k, k) = A023426(n).
Sum_{k=0..n} k * M(n, k) = 2*A014531(n-1) = 2*GegenbauerC(n - 2, -n, -1/2).
Sum_{k=0..n} i^k*M(n, k) = A343773(n), (i the imaginary unit), is the excess of the number of even Motzkin n-paths (A107587) over the odd ones (A343386).
Sum_{k=0..n} Sum_{j=0..k} M(n, j) = A189912(n).
Sum_{k=0..n} Sum_{j=0..k} M(n, n-j) = modified A025179(n).
For a recursion see the Python program.
EXAMPLE
Triangle M(n, k) starts:
[0] 1;
[1] 1, 0;
[2] 1, 0, 1;
[3] 1, 0, 3, 0;
[4] 1, 0, 6, 0, 2;
[5] 1, 0, 10, 0, 10, 0;
[6] 1, 0, 15, 0, 30, 0, 5;
[7] 1, 0, 21, 0, 70, 0, 35, 0;
[8] 1, 0, 28, 0, 140, 0, 140, 0, 14;
[9] 1, 0, 36, 0, 252, 0, 420, 0, 126, 0;
MAPLE
CatalanNumber := n -> binomial(2*n, n)/(n + 1):
M := (n, k) -> ifelse(irem(k, 2) = 1, 0, CatalanNumber(k/2)*binomial(n, k)):
for n from 0 to 9 do seq(M(n, k), k = 0..n) od;
# Alternative, as coefficients of polynomials:
p := n -> hypergeom([(1 - n)/2, -n/2], [2], (2*x)^2):
seq(print(seq(coeff(simplify(p(n)), x, k), k = 0..n)), n = 0..9);
# Using the exponential generating function:
egf := exp(x)*BesselI(1, 2*x*t)/(x*t): ser:= series(egf, x, 11):
seq(print(seq(coeff(simplify(n!*coeff(ser, x, n)), t, k), k = 0..n)), n = 0..9);
PROG
(Python)
from functools import cache
@cache
def M(n: int, k: int) -> int:
if k % 2: return 0
if n < 3: return 1
if n == k: return (2 * (n - 1) * M(n - 2, n - 2)) // (n // 2 + 1)
return (M(n - 1, k) * n) // (n - k)
for n in range(10): print([M(n, k) for k in range(n + 1)])
CROSSREFS
Cf. A001006 (Motzkin numbers), A026300 (Motzkin gf. triangle), A126120 (aerated Catalan), A000108 (Catalan).
KEYWORD
nonn,tabl
AUTHOR
Peter Luschny, Jan 09 2023
STATUS
approved
a(n) = 132*binomial(n,12).
(Formerly N2345)
+10
1
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 132, 1716, 12012, 60060, 240240, 816816, 2450448, 6651216, 16628040, 38798760, 85357272, 178474296, 356948592, 686439600, 1274816400, 2294669520, 4015671660, 6850263420, 11417105700, 18627909300, 29804654880, 46835886240, 72382733280, 110147637600, 165221456400, 244527755472, 357386719536, 516225261552, 737464659360
OFFSET
0,13
REFERENCES
Charles Jordan, Calculus of Finite Differences, Chelsea 1965, p. 450.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
LINKS
Index entries for linear recurrences with constant coefficients, signature (13,-78,286,-715,1287,-1716,1716,-1287,715,-286,78,-13,1).
FORMULA
G.f.: -132*x^12 / (x-1)^13. - Colin Barker, Jul 22 2013
a(n) = 132*A010965(n). - R. J. Mathar, Jul 15 2015
MATHEMATICA
132*Binomial[Range[0, 40], 12] (* Harvey P. Dale, Feb 16 2020 *)
PROG
(Magma) [132*Binomial(n, 12):n in [0..40]]; // Vincenzo Librandi, Feb 17 2020
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Jun 10 2012. This sequence was in the 1973 "Handbook", but was later dropped from the database. I have now reinstated it. - N. J. A. Sloane, Jun 10 2012
STATUS
approved

Search completed in 0.014 seconds