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

A008459 Square the entries of Pascal's triangle.

Original entry on oeis.org

1, 1, 1, 1, 4, 1, 1, 9, 9, 1, 1, 16, 36, 16, 1, 1, 25, 100, 100, 25, 1, 1, 36, 225, 400, 225, 36, 1, 1, 49, 441, 1225, 1225, 441, 49, 1, 1, 64, 784, 3136, 4900, 3136, 784, 64, 1, 1, 81, 1296, 7056, 15876, 15876, 7056, 1296, 81, 1, 1, 100, 2025, 14400, 44100, 63504, 44100, 14400, 2025, 100, 1
Offset: 0

Views

Author

Keywords

Comments

Number of lattice paths from (0, 0) to (n, n) with steps (1, 0) and (0, 1), having k right turns. - Emeric Deutsch, Nov 23 2003
Product of A007318 and A105868. - Paul Barry, Nov 15 2005
Number of partitions that fit in an n X n box with Durfee square k. - Franklin T. Adams-Watters, Feb 20 2006
From Peter Bala, Oct 23 2008: (Start)
Narayana numbers of type B. Row n of this triangle is the h-vector of the simplicial complex dual to an associahedron of type B_n (a cyclohedron) [Fomin & Reading, p. 60]. See A063007 for the corresponding f-vectors for associahedra of type B_n. See A001263 for the h-vectors for associahedra of type A_n. The Hilbert transform of this triangular array is A108625 (see A145905 for the definition of this term).
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 h-vectors of a unimodular triangulation of P(A_n) [Ardila et al.]. A063007 is the corresponding array of f-vectors for these type A_n polytopes. See A086645 for the array of h-vectors for type C_n polytopes and A108558 for the array of h-vectors associated with type D_n polytopes.
(End)
The n-th row consists of the coefficients of the polynomial P_n(t) = Integral_{s = 0..2*Pi} (1 + t^2 - 2*t*cos(s))^n/Pi/2 ds. For example, when n = 3, we get P_3(t) = t^6 + 9*t^4 + 9*t^2 + 1; the coefficients are 1, 9, 9, 1. - Theodore Kolokolnikov, Oct 26 2010
Let E(y) = Sum_{n >= 0} y^n/n!^2 = BesselJ(0, 2*sqrt(-y)). Then this triangle is the generalized Riordan array (E(y), y) with respect to the sequence n!^2 as defined in Wang and Wang. - Peter Bala, Jul 24 2013
From Colin Defant, Sep 16 2018: (Start)
Let s denote West's stack-sorting map. T(n,k) is the number of permutations pi of [n+1] with k descents such that s(pi) avoids the patterns 132, 231, and 321. T(n,k) is also the number of permutations pi of [n+1] with k descents such that s(pi) avoids the patterns 132, 312, and 321.
T(n,k) is the number of permutations of [n+1] with k descents that avoid the patterns 1342, 3142, 3412, and 3421. (End)
The number of convex polyominoes whose smallest bounding rectangle has size (k+1)*(n+1-k) and which contain the lower left corner of the bounding rectangle (directed convex polyominoes). - Günter Rote, Feb 27 2019
Let P be the poset [n] X [n] ordered by the product order. T(n,k) is the number of antichains in P containing exactly k elements. Cf. A063746. - Geoffrey Critzer, Mar 28 2020

Examples

			Pascal's triangle begins
  1
  1  1
  1  2   1
  1  3   3   1
  1  4   6   4   1
  1  5  10  10   5   1
  1  6  15  20  15   6   1
  1  7  21  35  35  21   7   1
...
so the present triangle begins
  1
  1   1
  1   4    1
  1   9    9     1
  1  16   36    16     1
  1  25  100   100    25    1
  1  36  225   400   225   36   1
  1  49  441  1225  1225  441  49   1
...
		

References

  • T. K. Petersen, Eulerian Numbers, Birkhauser, 2015, Chapter 12.
  • J. Riordan, An introduction to combinatorial analysis, Dover Publications, Mineola, NY, 2002, page 191, Problem 15. MR1949650
  • P. G. Tait, On the Linear Differential Equation of the Second Order, Proceedings of the Royal Society of Edinburgh, 9 (1876), 93-98 (see p. 97) [From Tom Copeland, Sep 09 2010, vol number corrected Sep 10 2010]

Crossrefs

Row sums are in A000984. Columns 0-3 are A000012, A000290, A000537, A001249.
Family of polynomials (see A062145): this sequence (c=1), A132813 (c=2), A062196 (c=3), A062145 (c=4), A062264 (c=5), A062190 (c=6).
Cf. A007318, A055133, A116647, A001263, A086645, A063007, A108558, A108625 (Hilbert transform), A145903, A181543, A086645 (logarithmic derivative), A105868 (inverse binomial transform), A093118.

Programs

  • GAP
    Flat(List([0..10],n->List([0..n],k->Binomial(n,k)^2))); # Muniru A Asiru, Mar 30 2018
    
  • Magma
    /* As triangle */ [[Binomial(n, k)^2: k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Dec 15 2016
    
  • Maple
    seq(seq(binomial(n, k)^2, k=0..n), n=0..10);
  • Mathematica
    Table[Binomial[n, k]^2, {n, 0, 11}, {k, 0, n}]//Flatten (* Alonso del Arte, Dec 08 2013 *)
  • Maxima
    create_list(binomial(n,k)^2,n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
    
  • Maxima
    T(n,k):=if n=k then 1 else if k=0 then 1 else T(n-1,k)*(n+k)/(n-k)+T(n-1,k-1); /* Vladimir Kruchinin, Oct 18 2014 */
    
  • Maxima
    A(x,y):=1/sqrt(1-2*x-2*x*y+x^2-2*x^2*y+x^2*y^2);
    taylor(x*A(x,y)+x*y*A(x,y)+sqrt(1+4*x^2*y*A(x,y)^2),x,0,7,y,0,7); /* Vladimir Kruchinin, Oct 23 2020 */
    
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, binomial(n, k)^2)}; /* Michael Somos, May 03 2004 */
    
  • PARI
    {T(n,k)=polcoeff(polcoeff(sum(m=0,n,(2*m)!/m!^2*x^(2*m)*y^m/(1-x-x*y+x*O(x^n))^(2*m+1)),n,x),k,y)} \\ Paul D. Hanna, Oct 31 2010
    
  • Python
    def A008459(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),n-comb(r+1,2))**2 # Chai Wah Wu, Nov 12 2024

Formula

T(n,k) = A007318(n,k)^2. - Sean A. Irvine, Mar 29 2018
E.g.f.: exp((1+y)*x)*BesselI(0, 2*sqrt(y)*x). - Vladeta Jovovic, Nov 17 2003
G.f.: 1/sqrt(1-2*x-2*x*y+x^2-2*x^2*y+x^2*y^2); g.f. for row n: (1-t)^n P_n[(1+t)/(1-t)] where the P_n's are the Legendre polynomials. - Emeric Deutsch, Nov 23 2003 [The original version of the bivariate g.f. has been modified with the roles of x and y interchanged so that now x corresponds to n and y to k. - Petros Hadjicostas, Oct 22 2017]
G.f. for column k is Sum_{j = 0..k} C(k, j)^2*x^(k+j)/(1 - x)^(2*k+1). - Paul Barry, Nov 15 2005
Column k has g.f. (x^k)*Legendre_P(k, (1+x)/(1-x))/(1 - x)^(k+1) = (x^k)*Sum_{j = 0..k} C(k, j)^2*x^j/(1 - x)^(2*k+1). - Paul Barry, Nov 19 2005
Let E be the operator D*x*D, where D denotes the derivative operator d/dx. Then (1/n!^2) * E^n(1/(1 - x)) = (row n generating polynomial)/(1 - x)^(2*n+1) = Sum_{k >= 0} binomial(n+k, k)^2*x^k. For example, when n = 3 we have (1/3!)^2*E^3(1/(1 - x)) = (1 + 9*x + 9*x^2 + x^3)/(1 - x)^7 = (1/3!)^2 * Sum_{k >= 0} ((k+1)*(k+2)*(k+3))^2*x^k. - Peter Bala, Oct 23 2008
G.f.: A(x, y) = Sum_{n >= 0} (2*n)!/n!^2 * x^(2*n)*y^n/(1 - x - x*y)^(2*n+1). - Paul D. Hanna, Oct 31 2010
From Peter Bala, Jul 24 2013: (Start)
Let E(y) = Sum_{n >= 0} y^n/n!^2 = BesselJ(0, 2*sqrt(-y)). Generating function: E(y)*E(x*y) = 1 + (1 + x)*y + (1 + 4*x + x^2)*y^2/2!^2 + (1 + 9*x + 9*x^2 + x^3)*y^3/3!^2 + .... Cf. the unsigned version of A021009 with generating function exp(y)*E(x*y).
The n-th power of this array has the generating function E(y)^n*E(x*y). In particular, the matrix inverse A055133 has the generating function E(x*y)/E(y). (End)
T(n,k) = T(n-1,k)*(n+k)/(n-k) + T(n-1,k-1), T(n,0) = T(n,n) = 1. - Vladimir Kruchinin, Oct 18 2014
Observe that the recurrence T(n,k) = T(n-1,k)*(n+k)/(n-k) - T(n-1,k-1), for n >= 2 and 1 <= k < n, with boundary conditions T(n,0) = T(n,n) = 1 gives Pascal's triangle A007318. - Peter Bala, Dec 21 2014
n-th row polynomial R(n, x) = [z^n] (1 + (1 + x)*z + x*z^2)^n. Note that 1/n*[z^(n-1)] (1 + (1 + x)*z + x*z^2)^n gives the row polynomials of A001263. - Peter Bala, Jun 24 2015
Binomial transform of A105868. If G(x,t) = 1/sqrt(1 - 2*(1 + t)*x + (1 - t)^2*x^2) denotes the o.g.f. of this array then 1 + x*d/dx log(G(x,t)) = 1 + (1 + t)*x + (1 + 6*t + t^2)*x^2 + ... is the o.g.f. for A086645. - Peter Bala, Sep 06 2015
T(n,k) = Sum_{i=0..n} C(n-i,k)*C(n,i)*C(n+i,i)*(-1)^(n-i-k). - Vladimir Kruchinin, Jan 14 2018
G.f. satisfies A(x,y) = x*A(x,y)+x*y*A(x,y)+sqrt(1+4*x^2*y*A(x,y)^2). - Vladimir Kruchinin, Oct 23 2020
G.f. satisfies the differential equation y * d^2(A(x,y))/dy^2 - x^2 * d^2(x*A(x,y))/dx^2 + 2*x^2* A(x,y)^3 = 0. - Sergii Voloshyn, Mar 07 2025
T(n,k) = Sum_{i=0..n} C(2*n+1,i)*C(n+k-i,n)^2*(-1)^i. - Natalia L. Skirrow, Apr 14 2025

A062196 Triangle read by rows, T(n, k) = binomial(n, k)*binomial(n + 2, k).

Original entry on oeis.org

1, 1, 3, 1, 8, 6, 1, 15, 30, 10, 1, 24, 90, 80, 15, 1, 35, 210, 350, 175, 21, 1, 48, 420, 1120, 1050, 336, 28, 1, 63, 756, 2940, 4410, 2646, 588, 36, 1, 80, 1260, 6720, 14700, 14112, 5880, 960, 45, 1, 99, 1980, 13860, 41580, 58212, 38808, 11880, 1485, 55
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Comments

Also the coefficient triangle of certain polynomials N(2; m,x) := Sum_{k=0..m} T(m,k)*x^k. The e.g.f. of the m-th (unsigned) column sequence without leading zeros of the generalized (a=2) Laguerre triangle L(2; n+m,m) = A062139(n+m,m), n >= 0, is N(2; m,x)/(1-x)^(3+2*m), with the row polynomials N(2; m,x).

Examples

			Triangle starts:
  n\k 0...1.....2......3..... 4.....;
  [0] 1;
  [1] 1,  3;
  [2] 1,  8,    6;
  [3] 1, 15,   30,    10;
  [4] 1, 24,   90,    80,    15;
  [5] 1, 35,  210,   350,   175,    21;
  [6] 1, 48,  420,  1120,  1050,   336,    28;
  [7] 1, 63,  756,  2940,  4410,  2646,   588,    36;
  [8] 1, 80, 1260,  6720, 14700, 14112,  5880,   960,   45;
  [9] 1, 99, 1980, 13860, 41580, 58212, 38808, 11880, 1485, 55.
		

Crossrefs

Family of polynomials (see A062145): A008459 (c=1), A132813 (c=2), this sequence (c=3), A062145 (c=4), A062264 (c=5), A062190 (c=6).
Sums include: A001791 (row), (-1)^n*A089849(n+1) (alternating sign row).
Diagonals: A000217 (k=n), A002417 (k=n-1), A001297 (k=n-2), A105946 (k=n-3), A105947 (k=n-4), A105948 (k=n-5), A107319 (k=n-6).
Columns: A005563 (k=1), A033487 (k=2), A027790 (k=3), A107395 (k=4), A107396 (k=5), A107397 (k=6), A107398 (k=7), A107399 (k=8).

Programs

  • Magma
    A062196:= func;
    [A062196(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Feb 21 2025
    
  • Maple
    T := (n, k) -> binomial(n, k)*binomial(n + 2, k);
    seq(seq(T(n, k), k=0..n), n=0..9); # Peter Luschny, Sep 30 2021
  • Mathematica
    A062196[n_, k_]:= Binomial[n, k]*Binomial[n+2, k];
    Table[A062196[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Feb 21 2025 *)
  • SageMath
    def A062196(n,k): return binomial(n,k)*binomial(n+2,k)
    print(flatten([[A062196(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Feb 21 2025

Formula

T(m, k) = [x^k] N(2; m, x), where N(2; m, x) = ((1-x)^(3+2*m))*(d^m/dx^m)(x^m/(m!*(1-x)^(m+3))).
N(2; m, x) = Sum_{j=0..m} ((binomial(m, j)*(2*m+2-j)!/((m+2)!*(m-j)!)*(x^(m-j)))*(1-x)^j).
T(n,m) = binomial(n, m)*(binomial(n+1, m) + binomial(n+1, m-1)). - Vladimir Kruchinin, Apr 06 2018
From G. C. Greubel, Feb 21 2025: (Start)
T(2*n, n) = (n+1)^2*A000108(n)*A000108(n+1).
T(2*n-1, n) = (4*n^2 - 1)*A000108(n-1)*A000108(n), n >= 1.
T(2*n+1, n) = (1/2)*binomial(n+2,2)*A000108(n+1)*A000108(n+2). (End)

Extensions

New name by Peter Luschny, Sep 30 2021

A062145 Triangle read by rows: T(n, k) = [z^k] P(n, z) where P(n, z) = Sum_{k=0..n} binomial(n, k) * Pochhammer(n - k + c, k) * z^k / k! and c = 4.

Original entry on oeis.org

1, 1, 4, 1, 10, 10, 1, 18, 45, 20, 1, 28, 126, 140, 35, 1, 40, 280, 560, 350, 56, 1, 54, 540, 1680, 1890, 756, 84, 1, 70, 945, 4200, 7350, 5292, 1470, 120, 1, 88, 1540, 9240, 23100, 25872, 12936, 2640, 165, 1, 108, 2376, 18480, 62370, 99792, 77616, 28512, 4455, 220
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Comments

Coefficient triangle of certain polynomials N(3; m,x).

Examples

			As a square array:
    1,    1,     1,     1,     1,     1,    1,  1, ... A000012;
    4,   10,    18,    28,    40,    54,   70, 88, ... A028552;
   10,   45,   126,   280,   540,   945, 1540, ....... A105938;
   20,  140,   560,  1680,  4200,  9240, ............. A105939;
   35,  350,  1890,  7350, 23100, 62370, ............. A027803;
   56,  756,  5292, 25872, 99792, .................... A105940;
   84, 1470, 12936, 77616, ........................... A105942;
  120, 2640, 28512, .................................. A105943;
  165, 4455, 57015, .................................. A105944;
  ....;
As a triangle:
  1;
  1,   4;
  1,  10,   10;
  1,  18,   45,    20;
  1,  28,  126,   140,    35;
  1,  40,  280,   560,   350,    56;
  1,  54,  540,  1680,  1890,   756,    84;
  1,  70,  945,  4200,  7350,  5292,  1470,   120;
  1,  88, 1540,  9240, 23100, 25872, 12936,  2640,  165;
  1, 108, 2376, 18480, 62370, 99792, 77616, 28512, 4455, 220;
  ....;
		

Crossrefs

Family of polynomials: A008459 (c=1), A132813 (c=2), A062196 (c=3), this sequence (c=4), A062264 (c=5), A062190 (c=6).
Columns: A028552 (k=1), A105938 (k=2), A105939 (k=3), A027803 (k=4), A105940 (k=5), A105942 (k=6), A105943 (k=7), A105944 (k=8).
Diagonals: A000292 (k=n), A027800 (k=n-1), A107417 (k=n-2), A107418 (k=n-3), A107419 (k=n-4), A107420 (k=n-5), A107421 (k=n-6), A107422 (k=n-7).
Sums: A002054 (row).

Programs

  • Magma
    A062145:= func< n,k | Binomial(n,k)*Binomial(n+3,k) >;
    [A062145(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Mar 07 2025
    
  • Mathematica
    NN[3, m_, x_] := x^m*(2*m+3)!*Hypergeometric2F1[-m, -m, -2*m-3, (x-1)/x]/( (m+3)!*m!); Table[CoefficientList[NN[3, m, x], x], {m, 0, 9}] // Flatten (* Jean-François Alcover, Sep 18 2013 *)
    P[c_, n_, z_] := Sum[Binomial[n, k] Pochhammer[n-k+c, k] z^k /k!, {k,0,n}];
    CL[c_] := Table[CoefficientList[P[c, n, z], z], {n, 0, 5}] // TableForm
    CL[4]  (* Peter Luschny, Feb 12 2024 *)
    A062145[n_,k_]:= Binomial[n,k]*Binomial[n+3,k];
    Table[A062145[n,k], {n,0,12},{k,0,n}]//Flatten (* G. C. Greubel, Mar 07 2025 *)
  • SageMath
    def A062145(n,k): return binomial(n,k)*binomial(n+3,k)
    print(flatten([[A062145(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Mar 07 2025

Formula

The e.g.f. of the m-th (unsigned) column sequence without leading zeros of the generalized (a=3) Laguerre triangle L(3; n+m, m) = A062137(n+m, m), n >= 0, is N(3; m, x)/(1-x)^(2*(m+2)), with the row polynomials N(3; m, x) := Sum_{k=0..m} a(m, k)*x^k.
N(3; m, x) := ((1-x)^(2*(m+2)))*(d^m/dx^m)(x^m/(m!*(1-x)^(m+4))); a(m, k) = [x^k]N(3; m, x).
N(3; m, x) = Sum_{j=0..m} ((binomial(m, j)*(2*m+3-j)!/((m+3)!*(m-j)!))*(x^(m-j))*(1-x)^j).
N(3; m, x)= x^m*(2*m+3)! * 2F1(-m, -m; -2*m-3; (x-1)/x)/((m+3)!*m!). - Jean-François Alcover, Sep 18 2013
From G. C. Greubel, Mar 07 2025 : (Start)
T(n, k) = binomial(n, k)*binomial(n+3, k).
T(2*n, n) = (1/2)*(n+1)^2*A000108(n)*A000108(n+2).
Sum_{k=0..n} (-1)^k*T(n, k) = (-1)^floor((n+2)/2)*(A047074(n+3) - A047074(n+ 2)). (End)

Extensions

New name by Peter Luschny, Feb 12 2024
More terms from G. C. Greubel, Mar 07 2025

A062264 Coefficient triangle of certain polynomials N(4; m,x).

Original entry on oeis.org

1, 1, 5, 1, 12, 15, 1, 21, 63, 35, 1, 32, 168, 224, 70, 1, 45, 360, 840, 630, 126, 1, 60, 675, 2400, 3150, 1512, 210, 1, 77, 1155, 5775, 11550, 9702, 3234, 330, 1, 96, 1848, 12320, 34650, 44352, 25872, 6336, 495, 1, 117, 2808, 24024, 90090, 162162, 144144, 61776, 11583, 715
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Comments

The e.g.f. of the m-th (unsigned) column sequence without leading zeros of the generalized (a=4) Laguerre triangle L(4; n+m,m) = A062140(n+m,m), n >= 0, is N(4; m,x)/(1-x)^(5+2*m), with the row polynomials N(4; m,x) := Sum_{k=0..m} T(m,k)*x^k.

Examples

			Triangle begins as:
  1;
  1,   5;
  1,  12,   15;
  1,  21,   63,    35;
  1,  32,  168,   224,     70;
  1,  45,  360,   840,    630,    126;
  1,  60,  675,  2400,   3150,   1512,    210;
  1,  77, 1155,  5775,  11550,   9702,   3234,    330;
  1,  96, 1848, 12320,  34650,  44352,  25872,   6336,    495;
  1, 117, 2808, 24024,  90090, 162162, 144144,  61776,  11583,   715;
  1, 140, 4095, 43680, 210210, 504504, 630630, 411840, 135135, 20020, 1001;
		

Crossrefs

Family of polynomials (see A062145): A008459 (c=1), A132813 (c=2), A062196 (c=3), A062145 (c=4), this sequence (c=5), A062190 (c=6).
Columns: A028347 (k=2), A104473 (k=3), A104474 (k=4), A104475 (k=5), A027814 (k=6), A103604 (k=7), A104476 (k=8), A104478 (k=9).
Diagonals: A000332 (k=n), A027810 (k=n-1), A105249 (k=n-2), A105250 (k=n-3), A105251 (k=n-4), A105252 (k=n-5), A105253 (k=n-6), A105254 (k=n-7).
Sums: A002694 (row).

Programs

  • Magma
    A062264:= func< n,k | Binomial(n,k)*Binomial(n+4,k) >;
    [A062264(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Mar 03 2025
    
  • Mathematica
    A062264[n_, k_]:= Binomial[n,k]*Binomial[n+4,k];
    Table[A062264[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Mar 03 2025 *)
  • SageMath
    def A062264(n,k): return binomial(n,k)*binomial(n+4,k)
    print(flatten([[A062264(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Mar 03 2025

Formula

T(m, k) = [x^k] N(4; m, x), with N(4; m, x) = ((1-x)^(2*m+5))*(d^m/dx^m)((x^m)/(m!*(1-x)^(m+5))).
N(4; m, x) = Sum_{j=0..m} (binomial(m, j)*(2*m+4-j)!/((m+4)!*(m-j)!)*(x^(m-j))*(1-x)^j).
From G. C. Greubel, Mar 03 2025: (Start)
T(n, k) = binomial(n,k)*binomial(n+4,k).
Sum_{k=0..n} (-1)^k*T(n, k) = (1/4)*( (1+(-1)^n)*(-1)^((n+2)/2)*(n^2 + 5*n - 2)*Catalan((n+2)/2)/(n+1) + 8*(1-(-1)^n)*(-1)^((n+1)/2)*Catalan((n+1)/2) ). (End)

A132813 Triangle read by rows: A001263 * A127648 as infinite lower triangular matrices.

Original entry on oeis.org

1, 1, 2, 1, 6, 3, 1, 12, 18, 4, 1, 20, 60, 40, 5, 1, 30, 150, 200, 75, 6, 1, 42, 315, 700, 525, 126, 7, 1, 56, 588, 1960, 2450, 1176, 196, 8, 1, 72, 1008, 4704, 8820, 7056, 2352, 288, 9, 1, 90, 1620, 10080, 26460, 31752, 17640, 4320, 405, 10
Offset: 0

Views

Author

Gary W. Adamson, Sep 01 2007

Keywords

Comments

Also T(n,k) = binomial(n-1, k-1)*binomial(n, k-1), related to Narayana polynomials (see Sulanke reference). - Roger L. Bagula, Apr 09 2008
h-vector for cluster complex associated to the root system B_n. See p. 8, Athanasiadis and C. Savvidou. - Tom Copeland, Oct 19 2014

Examples

			First few rows of the triangle are:
  1;
  1,  2;
  1,  6,   3;
  1, 12,  18,   4;
  1, 20,  60,  40,   5;
  1, 30, 150, 200,  75,   6;
  1, 42, 315, 700, 525, 126, 7;
  ...
		

Crossrefs

Family of polynomials (see A062145): A008459 (c=1), this sequence (c=2), A062196 (c=3), A062145 (c=4), A062264 (c=5), A062190 (c=6).
Columns: A000012 (k=0), A002378 (k=1), A006011 (k=2), 4*A006542 (k=3), 5*A006857 (k=4), 6*A108679 (k=5), 7*A134288 (k=6), 8*A134289 (k=7), 9*A134290 (k=8), 10*A134291 (k=9).
Diagonals: A000027 (k=n), A002411 (k=n-1), A004302 (k=n-2), A108647 (k=n-3), A134287 (k=n-4).
Main diagonal: A000894.
Sums: (-1)^floor((n+1)/2)*A001405 (signed row), A001700 (row), A203611 (diagonal).
Cf. A103371 (mirrored).

Programs

  • GAP
    Flat(List([0..10],n->List([0..n], k->(k+1)*Binomial(n+1,k+1)*Binomial(n+1,k)/(n+1)))); # Muniru A Asiru, Feb 26 2019
    
  • Haskell
    a132813 n k = a132813_tabl !! n !! k
    a132813_row n = a132813_tabl !! n
    a132813_tabl = zipWith (zipWith (*)) a007318_tabl $ tail a007318_tabl
    -- Reinhard Zumkeller, Apr 04 2014
    
  • Magma
    /* triangle */ [[(k+1)*Binomial(n+1,k+1)*Binomial(n+1,k)/(n+1): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Oct 19 2014
    
  • Maple
    P := (n, x) -> hypergeom([1-n, -n], [1], x): for n from 1 to 9 do PolynomialTools:-CoefficientList(simplify(P(n,x)),x) od; # Peter Luschny, Nov 26 2014
  • Mathematica
    T[n_,k_]=Binomial[n-1,k-1]*Binomial[n,k-1]; Table[Table[T[n,k],{k,1,n}],{n,1,11}]; Flatten[%] (* Roger L. Bagula, Apr 09 2008 *)
    P[n_, x_] := HypergeometricPFQ[{1-n, -n}, {1}, x]; Table[CoefficientList[P[n, x], x], {n, 1, 10}] // Flatten (* Jean-François Alcover, Nov 27 2014, after Peter Luschny *)
  • PARI
    tabl(nn) = {for (n = 1, nn, for (k = 1, n, print1(binomial(n-1, k-1)*binomial(n, k-1) , ", ");););} \\ Michel Marcus, Feb 12 2014
    
  • SageMath
    def A132813(n,k): return binomial(n,k)*binomial(n+1,k)
    print(flatten([[A132813(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Mar 12 2025

Formula

T(n,k) = (k+1)*binomial(n+1,k+1)*binomial(n+1,k)/(n+1), n >= k >= 0.
From Roger L. Bagula, May 14 2010: (Start)
T(n, m) = coefficients(p(x,n)), where
p(x,n) = (1-x)^(2*n)*Sum_{k >= 0} binomial(k+n-1, k)*binomial(n+k, k)*x^k,
or p(x,n) = (1-x)^(2*n)*Hypergeometric2F1([n, n+1], [1], x). (End)
T(n,k) = binomial(n,k) * binomial(n+1,k). - Reinhard Zumkeller, Apr 04 2014
These are the coefficients of the polynomials Hypergeometric2F1([1-n,-n], [1], x). - Peter Luschny, Nov 26 2014
G.f.: A(x,y) = A281260(x,y)/(1-A281260(x,y))/x. - Vladimir Kruchinin, Oct 10 2020

A027818 a(n) = (n+1)*binomial(n+6,6).

Original entry on oeis.org

1, 14, 84, 336, 1050, 2772, 6468, 13728, 27027, 50050, 88088, 148512, 241332, 379848, 581400, 868224, 1268421, 1817046, 2557324, 3542000, 4834830, 6512220, 8665020, 11400480, 14844375, 19143306, 24467184, 31011904, 39002216, 48694800, 60381552, 74393088, 91102473
Offset: 0

Views

Author

Thi Ngoc Dinh (via R. K. Guy)

Keywords

Comments

Number of 13-subsequences of [ 1, n ] with just 6 contiguous pairs.

Crossrefs

Cf. A093564 ((7, 1) Pascal, column m=7). Partial sums of A050403.

Programs

  • GAP
    List([0..30], n-> (n+1)*Binomial(n+6,6)); # G. C. Greubel, Aug 29 2019
  • Haskell
    a027818 n = (n + 1) * a007318' (n + 6) 6
    -- Reinhard Zumkeller, Aug 31 2014
    
  • Magma
    [(n+1)*Binomial(n+6,6): n in [0..30]]; // G. C. Greubel, Aug 29 2019
    
  • Maple
    seq((n+1)*binomial(n+6,6), n=0..30); # Zerinvary Lajos, Oct 19 2006
  • Mathematica
    Table[(n+1)*Binomial[n+6,6], {n,0,30}] (* G. C. Greubel, Aug 29 2019 *)
  • PARI
    a(n) = (n+1)*binomial(n+6,6) \\ Charles R Greathouse IV, Jun 11 2015
    
  • Sage
    [(n+1)*binomial(n+6,6) for n in (0..30)] # G. C. Greubel, Aug 29 2019
    

Formula

G.f.: (1+6*x)/(1-x)^8.
a(n) = A245334(n+6,6)/A000142(6). - Reinhard Zumkeller, Aug 31 2014
E.g.f.: (7! +9360*x +20520*x^2 +15000*x^3 +4650*x^4 +666*x^5 +43*x^6 + x^7)*exp(x)/7!. - G. C. Greubel, Aug 29 2019
From Amiram Eldar, Jan 28 2022: (Start)
Sum_{n>=0} 1/a(n) = Pi^2 - 5269/600.
Sum_{n>=0} (-1)^n/a(n) = Pi^2/2 - 512*log(2)/5 + 40189/600. (End)

A062152 Sixth (unsigned) column of triangle A062138 (generalized a=5 Laguerre).

Original entry on oeis.org

1, 66, 2772, 96096, 3027024, 90810720, 2663781120, 77630192640, 2270683134720, 67111301537280, 2013339046118400, 61498356317798400, 1916698771904716800, 61039483966811750400, 1988143192061868441600
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Examples

			a(2) = (2+5)! * binomial(2+10,10) / 5! = (5040 * 66) / 120 = 2772. - _Indranil Ghosh_, Feb 24 2017
		

Crossrefs

Programs

  • Magma
    [Factorial(n+5)*Binomial(n+10,10)/Factorial(5): n in [0..20]]; // G. C. Greubel, May 11 2018
  • Mathematica
    Table[(n+5)!*Binomial[n+10,10]/5!,{n,0,14}] (* Indranil Ghosh, Feb 24 2017 *)
  • PARI
    a(n) = (n+5)!*binomial(n+10,10)/5! \\ Indranil Ghosh, Feb 24 2017
    
  • Python
    import math
    f=math.factorial
    def C(n, r):return f(n)/f(r)/f(n-r)
    def A062152(n): return f(n+5)*C(n+10, 10)/f(5) # Indranil Ghosh, Feb 24 2017
    

Formula

a(n) = A062138(n+5, 5).
a(n) = (n+5)!*binomial(n+10, 10)/5!.
E.g.f.: N(5;5, x)/(1-x)^16 with N(5;5, x) := Sum_{k=0..5} A062190(5, k)* x^k = 1 + 50*x + 450*x^2 + 1200*x^3 + 1050*x^4 + 252*x^5.
If we define f(n,i,x) = Sum_{k=i..n} Sum_{j=i..k} binomial(k,j)* Stirling1(n,k)*Stirling2(j,i)*x^(k-j) then a(n-10) = (-1)^n*f(n,10,-6), (n>=10). - Milan Janjic, Mar 01 2009

A104673 a(n) = binomial(n+5, n)*binomial(n+10, n).

Original entry on oeis.org

1, 66, 1386, 16016, 126126, 756756, 3699696, 15402816, 56316546, 184940756, 554822268, 1540663488, 4001445448, 9802357488, 22805484768, 50678855040, 108088495515, 222161129190, 441579528390, 851355545040, 1596291646950, 2917485413700, 5208073135200
Offset: 0

Views

Author

Zerinvary Lajos, Apr 22 2005

Keywords

Examples

			If n=0 then C(0+5,5)*C(0+10,0+0) = C(5,5)*C(10,0) = 1*1 = 1.
If n=6 then C(6+5,5)*C(6+10,6+0) = C(11,5)*C(16,6) = 462*8008 = 3699696.
		

Crossrefs

Cf. A062190.

Programs

  • Magma
    A104673:= func< n | Binomial(n+5,n)*Binomial(n+10,n) >;
    [A104673(n): n in [0..30]]; // G. C. Greubel, Mar 01 2025
    
  • Mathematica
    Table[Binomial[n+5,5]Binomial[n+10,n],{n,0,30}]
  • SageMath
    def A104673(n): return binomial(n+5,n)*binomial(n+10,n)
    print([A104673(n) for n in range(31)]) # G. C. Greubel, Mar 01 2025

Formula

G.f.: (1 + 50*x + 450*x^2 + 1200*x^3 + 1050*x^4 + 252*x^5)/(1-x)^16. - Harvey P. Dale, Nov 24 2011
From Amiram Eldar, Aug 30 2022: (Start)
Sum_{n>=0} 1/a(n) = 17875*Pi^2/3 - 7468753315/127008.
Sum_{n>=0} (-1)^n/a(n) = 208828891/127008 - 725*Pi^2/6 - 40960*log(2)/63. (End)

Extensions

Corrected and extended by Harvey P. Dale, Nov 24 2011

A104678 a(n) = binomial(n+4,4) * binomial(n+9,4).

Original entry on oeis.org

126, 1050, 4950, 17325, 50050, 126126, 286650, 600600, 1178100, 2187900, 3879876, 6613425, 10892700, 17409700, 27096300, 41186376, 61289250, 89475750, 128378250, 181306125, 252378126, 346673250, 470401750, 631098000, 837837000, 1101476376, 1434925800
Offset: 0

Views

Author

Zerinvary Lajos, Apr 22 2005

Keywords

Examples

			If n=0 then C(0+4,0+0)*C(0+9,4) = C(4,0)*C(9,4) = 1*126 = 126.
If n=5 then C(5+4,5+0)*C(5+9,4) = C(9,5)*C(14,4) = 126*1001 = 126126.
		

Crossrefs

Programs

Formula

From R. J. Mathar, Nov 29 2015: (Start)
a(n) = A000332(n+4)*A000332(n+9).
G.f.: ( -126+84*x-36*x^2+9*x^3-x^4 ) / (x-1)^9. (End)
From Amiram Eldar, Aug 30 2022: (Start)
Sum_{n>=0} 1/a(n) = 9/980.
Sum_{n>=0} (-1)^n/a(n) = 1/140. (End)

A104670 a(n) = binomial(n+2, 2)*binomial(n+7, n).

Original entry on oeis.org

1, 24, 216, 1200, 4950, 16632, 48048, 123552, 289575, 629200, 1283568, 2482272, 4585308, 8139600, 13953600, 23193984, 37509021, 59183784, 91333000, 138138000, 205134930, 299562120, 430775280, 610740000, 854611875, 1181415456, 1614834144, 2184124096, 2925166200
Offset: 0

Views

Author

Zerinvary Lajos, Apr 22 2005

Keywords

Examples

			If n=0 then C(2+0,2)*C(7+0,0+0) = C(2,2)*C(7,0) = 1*1 = 1;
if n=6 then C(2+6,2)*C(7+6,0+6) = C(8,2)*C(13,6) = 28*1716 = 48048.
		

Crossrefs

Programs

  • Magma
    A104670:= func< n | Binomial(n+2,n)*Binomial(n+7,n) >;
    [A104670(n): n in [0..30]]; // G. C. Greubel, Mar 01 2025
    
  • Maple
    [seq(stirling2(n+1,n)*binomial(n+6,7),n=1..25)]; # Zerinvary Lajos, Dec 06 2006
  • Mathematica
    a[n_] := Binomial[n + 2, 2] * Binomial[n + 7, 7]; Array[a, 25, 0] (* Amiram Eldar, Aug 30 2022 *)
  • SageMath
    def A104670(n): return binomial(n+2,n)*binomial(n+7,n)
    print([A104670(n) for n in range(31)]) # G. C. Greubel, Mar 01 2025

Formula

G.f.: (1 + 14*x + 21*x^2)/(1-x)^10. - Colin Barker, Mar 18 2012
From Amiram Eldar, Aug 30 2022: (Start)
Sum_{n>=0} 1/a(n) = 49*Pi^2/3 - 288281/1800.
Sum_{n>=0} (-1)^n/a(n) = 448*log(2)/3 - 35*Pi^2/6 - 1799/40. (End)

Extensions

Corrected and extended by Don Reble, Nov 21 2006
Showing 1-10 of 28 results. Next