A103371 Number triangle T(n,k) = C(n,n-k)*C(n+1,n-k).
1, 2, 1, 3, 6, 1, 4, 18, 12, 1, 5, 40, 60, 20, 1, 6, 75, 200, 150, 30, 1, 7, 126, 525, 700, 315, 42, 1, 8, 196, 1176, 2450, 1960, 588, 56, 1, 9, 288, 2352, 7056, 8820, 4704, 1008, 72, 1, 10, 405, 4320, 17640, 31752, 26460, 10080, 1620, 90, 1, 11, 550, 7425, 39600, 97020
Offset: 0
Examples
The triangle T(n, k) begins: n\k 0 1 2 3 4 5 6 7 8 9 ... 0: 1 1: 2 1 2: 3 6 1 3: 4 18 12 1 4: 5 40 60 20 1 5: 6 75 200 150 30 1 6: 7 126 525 700 315 42 1 7: 8 196 1176 2450 1960 588 56 1 8: 9 288 2352 7056 8820 4704 1008 72 1 9: 10 405 4320 17640 31752 26460 10080 1620 90 1 ... reformatted. - _Wolfdieter Lang_, Jul 31 2017 From _R. J. Mathar_, Mar 29 2013: (Start) The matrix inverse starts 1; -2, 1; 9, -6, 1; -76, 54, -12, 1; 1055, -760, 180, -20, 1; -21906, 15825, -3800, 450, -30, 1; 636447, -460026, 110775, -13300, 945, -42, 1; (End) O.g.f. of 4th diagonal [4, 40,200, ...] is G(3, x) = 4*(1 + 3*x + x^2)/(1 - x)^7, from the n = 3 row [1, 3, 1] of A001263. See a comment above. - _Wolfdieter Lang_, Jul 31 2017
Links
- Reinhard Zumkeller, Rows n = 0..125 of table, flattened
- Per Alexandersson, Svante Linusson, Samu Potka, and Joakim Uhlin, Refined Catalan and Narayana cyclic sieving, arXiv:2010.11157 [math.CO], 2020.
- Paul Barry and A. Hennessy, A Note on Narayana Triangles and Related Polynomials, Riordan Arrays, and MIMO Capacity Calculations, J. Int. Seq. 14 (2011), Article 11.3.8.
- Robert Coquereaux and Jean-Bernard Zuber, Counting partitions by genus. II. A compendium of results, arXiv:2305.01100 [math.CO], 2023. See p. 16.
- R. Cori and G. Hetyei, Counting genus one partitions and permutations, arXiv:1306.4628 [math.CO], 2013.
- Wolfdieter Lang, On Generating functions of Diagonals Sequences of Sheffer and Riordan Number Triangles, arXiv:1708.01421 [math.NT], August 2017.
- A. Laradji and A. Umar, Combinatorial results for semigroups of order-preserving full transformations, Semigroup Forum 72 (2006), 51-62.
Crossrefs
Programs
-
Haskell
a103371 n k = a103371_tabl !! n !! k a103371_row n = a103371_tabl !! n a103371_tabl = map reverse a132813_tabl -- Reinhard Zumkeller, Apr 04 2014
-
Magma
/* As triangle */ [[Binomial(n,n-k)*Binomial(n+1,n-k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Aug 01 2017
-
Maple
A103371 := (n,k) -> binomial(n,k)^2*(n+1)/(k+1); seq(print(seq(A103371(n, k), k=0..n)), n=0..7); # Peter Luschny, Oct 19 2011
-
Mathematica
Flatten[Table[Binomial[n,n-k]Binomial[n+1,n-k],{n,0,10},{k,0,n}]] (* Harvey P. Dale, May 26 2014 *) CoefficientList[Series[Series[E^(x(1+y))(BesselI[0,2*x*Sqrt[y]]+BesselI[1,2*x*Sqrt[y]]/Sqrt[y]),{x,0,8}],{y,0,8}],{x,y}]*Range[0,8]! (* Natalia L. Skirrow, Apr 14 2025 *)
-
Maxima
create_list(binomial(n,k)*binomial(n+1,k+1),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
-
PARI
for(n=0,10, for(k=0,n, print1(binomial(n,k)*binomial(n+1,k+1), ", "))) \\ G. C. Greubel, Nov 09 2018
Formula
Number triangle T(n, k) = C(n, n-k)*C(n+1, n-k) = C(n, k)*C(n+1, k+1); Column k of this triangle has g.f. Sum_{j=0..k} (C(k, j)*C(k+1, j) * x^(k+j))/(1-x)^(2*k+2); coefficients of the numerators are the rows of the reverse triangle C(n, k)*C(n+1, k).
T(n,k) = C(n, k)*Sum_{j=0..(n-k)} C(n-j, k). - Paul Barry, Jan 12 2006
T(n,k) = (n+1-k)*N(n+1,k+1), with N(n,k):=A001263(n,k), the Narayana triangle (with offset [1,1]).
O.g.f.: ((1-(1-y)*x)/sqrt((1-(1+y)*x)^2-4*x^2*y) -1)/2, (from o.g.f. of A001263, Narayana triangle). - Wolfdieter Lang, Nov 13 2007
From Peter Bala, Jan 24 2008: (Start)
O.g.f. for row n: (1-x)^n*P(n,1,0,(1+x)/(1-x)) = 1/(2*x)*(1-x)^(n+1)*( Legendre_P(n+1,(1+x)/(1-x)) - Legendre_P(n,(1+x)/(1-x)) ), where P(n,a,b,x) denotes the Jacobi polynomial.
O.g.f. for column k: x^k/(1-x)^(k+2)*P(k,0,1,(1+x)/(1-x)). Compare with A008459. (End)
Let S(n,k) = binomial(2*n,n)^(k+1)*((n+1)^(k+1)-n^(k+1))/(n+1)^k. Then T(2*n,n) = S(n,1). (Cf. A194595, A197653, A197654). - Peter Luschny, Oct 20 2011
T(n,k) = A003056(n+1,k+1)*C(n,k)^2/(k+1). - Peter Luschny, Oct 29 2011
From Natalia L. Skirrow, Apr 14 2025: (Start)
T(n,k) = A008459(n,k) + n*N(n,k+1).
E.g.f.: e^(x*(1+y))*(I_0(2*x*sqrt(y)) + I_1(2*x*sqrt(y))/sqrt(y)), where I_n is the modified Bessel function of the first kind. (The I_0 contributes A008459(n,k), the I_1 contributes n*N(n,k+1))
O.g.f. for row n: (n+1)*2F1(-n,-n;2;y) = (n+1)*2F1(2+n,2+n;2;y)*(1-y)^(2*(n+1)) (by Euler's hypergeometric transformation); (n+1)*2F1(2+n,2+n;2;y) is the o.g.f. for row n of (k+n+1)!^2/(k!*(k+1)!*n!*(n+1)!), which is column n+1 of A132812.
O.g.f. for column k: 2F1(1+k,2+k;1;x)*x^k = 2F1(-k,-1-k;1;x)*x^k/(1-x)^(2+2*k). 2F1(-k,-1-k;1;x) is the kth row of A132813, the reflection of the kth row of this triangle.
O.g.f. for diagonal d (beginning at a(d,0)): (d+1)*x^d*2F1(d+1,d+2;2;x*y). 2F1(d+1,d+2;2;x) = 2F1(1-d,-d;2;x)/(1-x)^(2*d+1), numerator being the o.g.f. of row d of the Narayana triangle.
These respectively yield:
T(n,k) = Sum_{i=0..n+k} C(2*(n+1),i)*(-1)^i*A132812(n+1+k-i,n+1),
T(d+k,k) = Sum_{i=0..k} C(d-i+1+2*k,d-i)*T(k,k-i),
T(d+k,k) = Sum_{i=0..d} C(k-i + 2*d,k-i)*N(d,i+1)*(d+1).
E.g.f. for column k: 1F1(2+k;1;x)*x^k/k!.
E.g.f. for diagonal d: (d+1)*x^d*1F1(d+2;2;x*y)/d!. (End)
Comments