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-8 of 8 results.

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

A086645 Triangle read by rows: T(n, k) = binomial(2n, 2k).

Original entry on oeis.org

1, 1, 1, 1, 6, 1, 1, 15, 15, 1, 1, 28, 70, 28, 1, 1, 45, 210, 210, 45, 1, 1, 66, 495, 924, 495, 66, 1, 1, 91, 1001, 3003, 3003, 1001, 91, 1, 1, 120, 1820, 8008, 12870, 8008, 1820, 120, 1, 1, 153, 3060, 18564, 43758, 43758, 18564, 3060, 153, 1, 1, 190, 4845, 38760
Offset: 0

Views

Author

Philippe Deléham, Jul 26 2003

Keywords

Comments

Terms have the same parity as those of Pascal's triangle.
Coefficients of polynomials (1/2)*((1 + x^(1/2))^(2n) + (1 - x^(1/2))^(2n)).
Number of compositions of 2n having k parts greater than 1; example: T(3, 2) = 15 because we have 4+2, 2+4, 3+2+1, 3+1+2, 2+3+1, 2+1+3, 1+3+2, 1+2+3, 2+2+1+1, 2+1+2+1, 2+1+1+2, 1+2+2+1, 1+2+1+2, 1+1+2+2, 3+3. - Philippe Deléham, May 18 2005
Number of binary words of length 2n - 1 having k runs of consecutive 1's; example: T(3,2) = 15 because we have 00101, 01001, 01010, 01011, 01101, 10001, 10010, 10011, 10100, 10110, 10111, 11001, 11010, 11011, 11101. - Philippe Deléham, May 18 2005
Let M_n be the n X n matrix M_n(i, j) = T(i, j-1); then for n > 0, det(M_n) = A000364(n), Euler numbers; example: det([1, 1, 0, 0; 1, 6, 1, 0; 1, 15, 15, 1; 1, 28, 70, 28 ]) = 1385 = A000364(4). - Philippe Deléham, Sep 04 2005
Equals ConvOffsStoT transform of the hexagonal numbers, A000384: (1, 6, 15, 28, 45, ...); e.g., ConvOffs transform of (1, 6, 15, 28) = (1, 28, 70, 28, 1). - Gary W. Adamson, Apr 22 2008
From Peter Bala, Oct 23 2008: (Start)
Let C_n be the root lattice generated as a monoid by {+-2*e_i: 1 <= i <= n; +-e_i +- e_j: 1 <= i not equal to j <= n}. Let P(C_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(C_n) [Ardila et al.]. See A127674 for (a signed version of) the corresponding array of f-vectors for these type C_n polytopes. See A008459 for the array of h-vectors for type A_n polytopes and A108558 for the array of h-vectors associated with type D_n polytopes.
The Hilbert transform of this triangle is A142992 (see A145905 for the definition of this term).
(End)
Diagonal sums: A108479. - Philippe Deléham, Sep 08 2009
Coefficients of Product_{k=1..n} (cot(k*Pi/(2n+1))^2 - x) = Sum_{k=0..n} (-1)^k*binomial(2n,2k)*x^k/(2n+1-2k). - David Ingerman (daviddavifree(AT)gmail.com), Mar 30 2010
Generalized Narayana triangle for 4^n (or cosh(2x)). - Paul Barry, Sep 28 2010
Coefficients of the matrix inverse appear to be T^(-1)(n,k) = (-1)^(n+k)*A086646(n,k). - R. J. Mathar, Mar 12 2013
Let E(y) = Sum_{n>=0} y^n/(2*n)! = cosh(sqrt(y)). Then this triangle is the generalized Riordan array (E(y), y) with respect to the sequence (2*n)! as defined in Wang and Wang. Cf. A103327. - Peter Bala, Aug 06 2013
Row 6, (1,66,495,924,495,66,1), plays a role in expansions of powers of the Dedekind eta function. See the Chan link, p. 534, and A034839. - Tom Copeland, Dec 12 2016

Examples

			From _Peter Bala_, Oct 23 2008: (Start)
The triangle begins
n\k|..0.....1.....2.....3.....4.....5.....6
===========================================
0..|..1
1..|..1.....1
2..|..1.....6.....1
3..|..1....15....15.....1
4..|..1....28....70....28.....1
5..|..1....45...210...210....45.....1
6..|..1....66...495...924...495....66.....1
...
(End)
From _Peter Bala_, Aug 06 2013: (Start)
Viewed as the generalized Riordan array (cosh(sqrt(y)), y) with respect to the sequence (2*n)! the column generating functions begin
1st col: cosh(sqrt(y)) = 1 + y/2! + y^2/4! + y^3/6! + y^4/8! + ....
2nd col: 1/2!*y*cosh(sqrt(y)) = y/2! + 6*y^2/4! + 15*y^3/6! + 28*y^4/8! + ....
3rd col: 1/4!*y^2*cosh(sqrt(y)) = y^2/4! + 15*y^3/6! + 70*y^4/8! + 210*y^5/10! + .... (End)
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 224.

Crossrefs

Cf. A008459, A108558, A127674, A142992. - Peter Bala, Oct 23 2008
Cf. A103327 (binomial(2n+1, 2k+1)), A103328 (binomial(2n, 2k+1)), A091042 (binomial(2n+1, 2k)). -Wolfdieter Lang, Jan 06 2013
Cf. A086646 (unsigned matrix inverse), A103327.
Cf. A034839.

Programs

  • Magma
    /* As triangle: */ [[Binomial(2*n, 2*k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Dec 14 2016
  • Maple
    A086645:=(n,k)->binomial(2*n,2*k): seq(seq(A086645(n,k),k=0..n),n=0..12);
  • Mathematica
    Table[Binomial[2 n, 2 k], {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Dec 13 2016 *)
  • Maxima
    create_list(binomial(2*n,2*k),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
    
  • PARI
    {T(n, k) = binomial(2*n, 2*k)};
    
  • PARI
    {T(n, k) = sum( i=0, min(k, n-k), 4^i * binomial(n, 2*i) * binomial(n - 2*i, k-i))}; /* Michael Somos, May 26 2005 */
    

Formula

T(n, k) = (2*n)!/((2*(n-k))!*(2*k)!) row sums = A081294. COLUMNS: A000012, A000384
Sum_{k>=0} T(n, k)*A000364(k) = A000795(n) = (2^n)*A005647(n).
Sum_{k>=0} T(n, k)*2^k = A001541(n). Sum_{k>=0} T(n, k)*3^k = 2^n*A001075(n). Sum_{k>=0} T(n, k)*4^k = A083884(n). - Philippe Deléham, Feb 29 2004
O.g.f.: (1 - z*(1+x))/(x^2*z^2 - 2*x*z*(1+z) + (1-z)^2) = 1 + (1 + x)*z +(1 + 6*x + x^2)*z^2 + ... . - Peter Bala, Oct 23 2008
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A081294(n), A001541(n), A090965(n), A083884(n), A099140(n), A099141(n), A099142(n), A165224(n), A026244(n) for x = 0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Sep 08 2009
Product_{k=1..n} (cot(k*Pi/(2n+1))^2 - x) = Sum_{k=0..n} (-1)^k*binomial(2n,2k)*x^k/(2n+1-2k). - David Ingerman (daviddavifree(AT)gmail.com), Mar 30 2010
From Paul Barry, Sep 28 2010: (Start)
G.f.: 1/(1-x-x*y-4*x^2*y/(1-x-x*y)) = (1-x*(1+y))/(1-2*x*(1+y)+x^2*(1-y)^2);
E.g.f.: exp((1+y)*x)*cosh(2*sqrt(y)*x);
T(n,k) = Sum_{j=0..n} C(n,j)*C(n-j,2*(k-j))*4^(k-j). (End)
T(n,k) = 2*T(n-1,k) + 2*T(n-1,k-1) + 2*T(n-2,k-1) - T(n-2,k) - T(n-2,k-2), with T(0,0)=T(1,0)=T(1,1)=1, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Nov 26 2013
From Peter Bala, Sep 22 2021: (Start)
n-th row polynomial R(n,x) = (1-x)^n*T(n,(1+x)/(1-x)), where T(n,x) is the n-th Chebyshev polynomial of the first kind. Cf. A008459.
R(n,x) = Sum_{k = 0..n} binomial(n,2*k)*(4*x)^k*(1 + x)^(n-2*k).
R(n,x) = n*Sum_{k = 0..n} (n+k-1)!/((n-k)!*(2*k)!)*(4*x)^k*(1-x)^(n-k) for n >= 1. (End)

A108553 Square array, read by antidiagonals, where row n equals the crystal ball sequence for D_n lattice.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 5, 1, 1, 13, 13, 7, 1, 1, 25, 55, 25, 9, 1, 1, 41, 169, 147, 41, 11, 1, 1, 61, 411, 625, 309, 61, 13, 1, 1, 85, 853, 2051, 1681, 561, 85, 15, 1, 1, 113, 1583, 5577, 6981, 3721, 923, 113, 17, 1, 1, 145, 2705, 13203, 23673, 18733, 7225, 1415, 145, 19, 1
Offset: 0

Views

Author

Paul D. Hanna, Jun 10 2005

Keywords

Comments

Rows 0 and 2 are included by extension since they fit the formula. Row 1 equals the odd numbers in order that triangle A108556 maintains that A108556(n,n-1) = (n/2)*A108556(n,n) for all n>=1, where row n of triangle A108556 equals the inverse binomial transform of row n of this square array.

Examples

			Square array begins:
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,...
1,3,5,7,9,11,13,15,17,19,21,23,25,27,...
1,5,13,25,41,61,85,113,145,181,221,265,...
1,13,55,147,309,561,923,1415,2057,2869,...
1,25,169,625,1681,3721,7225,12769,21025,...
1,41,411,2051,6981,18733,42783,86983,...
1,61,853,5577,23673,76389,204205,476113,...
1,85,1583,13203,68853,264825,824083,...
Inverse binomial transform of rows gives
rows of triangle A108556:
1,
1,2,
1,4,4,
1,12,30,20,
1,24,120,192,96,
1,40,330,940,1080,432, ...
Product of the g.f. of row n and (1-x)^(n+1)
generates the symmetric triangle A108558:
1;
1,1;
1,2,1;
1,9,9,1;
1,20,54,20,1;
1,35,180,180,35,1; ...
The row sums of triangle A108558 equals the
main diagonal of triangle A108556.
		

Crossrefs

Cf. A108554 (diagonal), A108555 (antidiagonal sums), A108556, A108558, A001844 (row 2), A005902 (row 3), A007204 (row 4), A008356 (row 5), A008358 (row 6), A008360 (row 7), A008362 (row 8), A008377 (row 9), A008379 (row 10).

Programs

  • PARI
    T(n,k)=if(n<0 || k<0,0,if(n==0 || k==0,1,if(n==1,2*k+1, sum(j=0,k,binomial(n+k-j,k-j)*(binomial(2*n,2*j)-2*n*binomial(n-2,j-1))))))

Formula

T(n, k) = Sum_{j=0..n} C(n+k-j, k-j)*[C(2*n, 2*j) - 2*j*(n-j)*C(n, j)/(n-1)] for n>1, with T(0, k)=1, T(1, k)=2*k+1.

A145905 Square array read by antidiagonals: Hilbert transform of triangle A060187.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 9, 5, 1, 1, 27, 25, 7, 1, 1, 81, 125, 49, 9, 1, 1, 243, 625, 343, 81, 11, 1, 1, 729, 3125, 2401, 729, 121, 13, 1, 1, 2187, 15625, 16807, 6561, 1331, 169, 15, 1, 1, 6561, 78125, 117649, 59049, 14641, 2197, 225, 17, 1, 1, 19683, 390625, 823543
Offset: 0

Views

Author

Peter Bala, Oct 27 2008

Keywords

Comments

Definition of the Hilbert transform of a triangular array:
For many square arrays in the database the entries in a row are polynomial in the column index, of degree d say and hence the row generating function has the form P(x)/(1-x)^(d+1), where P is some polynomial function. Often the array whose rows are formed from the coefficients of these P polynomials is of independent interest. This suggests the following definition.
Let [L(n,k)]n,k>=0 be a lower triangular array and let R(n,x) := sum {k = 0 .. n} L(n,k)*x^k, denote the n-th row generating polynomial of L. Then we define the Hilbert transform of L, denoted Hilb(L), to be the square array whose n-th row, n >= 0, has the generating function R(n,x)/(1-x)^(n+1).
In this particular case, L is the array A060187, the array of Eulerian numbers of type B, whose row polynomials are the h-polynomials for permutohedra of type B. The Hilbert transform is an infinite Vandermonde matrix V(1,3,5,...).
We illustrate the Hilbert transform with a few examples:
(1) The Delannoy number array A008288 is the Hilbert transform of Pascal's triangle A007318 (view as the array of coefficients of h-polynomials of n-dimensional cross polytopes).
(2) The transpose of the array of nexus numbers A047969 is the Hilbert transform of the triangle of Eulerian numbers A008292 (best viewed in this context as the coefficients of h-polynomials of n-dimensional permutohedra of type A).
(3) The sequence of Eulerian polynomials begins [1, x, x + x^2, x + 4*x^2 + x^3, ...]. The coefficients of these polynomials are recorded in triangle A123125, whose Hilbert transform is A004248 read as square array.
(4) A108625, the array of crystal ball sequences for the A_n lattices, is the Hilbert transform of A008459 (viewed as the triangle of coefficients of h-polynomials of n-dimensional associahedra of type B).
(5) A142992, the array of crystal ball sequences for the C_n lattices, is the Hilbert transform of A086645, the array of h-vectors for type C root polytopes.
(6) A108553, the array of crystal ball sequences for the D_n lattices, is the Hilbert transform of A108558, the array of h-vectors for type D root polytopes.
(7) A086764, read as a square array, is the Hilbert transform of the rencontres numbers A008290.
(8) A143409 is the Hilbert transform of triangle A073107.

Examples

			Triangle A060187 (with an offset of 0) begins
1;
1, 1;
1, 6, 1;
so the entries in the first three rows of the Hilbert transform of
A060187 come from the expansions:
Row 0: 1/(1-x) = 1 + x + x^2 + x^3 + ...;
Row 1: (1+x)/(1-x)^2 = 1 + 3*x + 5*x^2 + 7*x^3 + ...;
Row 2: (1+6*x+x^2)/(1-x)^3 = 1 + 9*x + 25*x^2 + 49*x^3 + ...;
The array begins
n\k|..0....1.....2.....3......4
================================
0..|..1....1.....1.....1......1
1..|..1....3.....5.....7......9
2..|..1....9....25....49.....81
3..|..1...27...125...343....729
4..|..1...81...625..2401...6561
5..|..1..243..3125.16807..59049
...
		

Crossrefs

Cf. A008292, A039755, A052750 (first superdiagonal), A060187, A114172, A145901.

Programs

  • Maple
    T:=(n,k) -> (2*k + 1)^n: seq(seq(T(n-k,k),k = 0..n),n = 0..10);

Formula

T(n,k) = (2*k + 1)^n, (see equation 4.10 in [Franssens]). This array is the infinite Vandermonde matrix V(1,3,5,7, ....) having a LDU factorization equal to A039755 * diag(2^n*n!) * transpose(A007318).

A108556 Triangle, read by rows, where row n equals the inverse binomial transform of the crystal ball sequence for D_n lattice.

Original entry on oeis.org

1, 1, 2, 1, 4, 4, 1, 12, 30, 20, 1, 24, 120, 192, 96, 1, 40, 330, 940, 1080, 432, 1, 60, 732, 3200, 6240, 5568, 1856, 1, 84, 1414, 8708, 25200, 37184, 27104, 7744, 1, 112, 2480, 20352, 80960, 173824, 206080, 126976, 31744, 1, 144, 4050, 42588, 221544, 643824, 1096032, 1085760, 579456, 128768
Offset: 0

Views

Author

Paul D. Hanna, Jun 10 2005

Keywords

Comments

Row n equals the inverse binomial transform of row n of the square array A108553.
Array of f-vectors for type D root polytopes [Ardila et al.]. See A063007 and A127674 for the arrays of f-vectors for type A and type C root polytopes respectively. - Peter Bala, Oct 23 2008

Examples

			Triangle begins:
1;
1,2;
1,4,4;
1,12,30,20;
1,24,120,192,96;
1,40,330,940,1080,432;
1,60,732,3200,6240,5568,1856;
1,84,1414,8708,25200,37184,27104,7744;
1,112,2480,20352,80960,173824,206080,126976,31744; ...
		

Crossrefs

Cf. A108553, A108557 (row sums), A108558, Rows are inverse binomial transforms of: A001844 (row 2), A005902 (row 3), A007204 (row 4), A008356 (row 5), A008358 (row 6), A008360 (row 7), A008362 (row 8), A008377 (row 9), A008379 (row 10).

Programs

  • Mathematica
    T[n_, k_] := Module[{A}, A = Table[Table[If[r - 1 == 0 || c - 1 == 0, 1, If[r - 1 == 1, 2c - 1, Sum[Binomial[r + c - j - 2, c - j - 1] (Binomial[2r - 2, 2j] - 2(r - 1) Binomial[r - 3, j - 1]), {j, 0, c - 1}]]], {c, 1, n + 1}], {r, 1, n + 1}]; SeriesCoefficient[((A[[n + 1]]. x^Range[0, n]) /. x -> x/(1 + x))/(1 + x), {x, 0, k}]];
    Table[T[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 26 2018, from PARI *)
  • PARI
    T(n,k)=local(A=vector(n+1,r,vector(n+1,c,if(r-1==0 || c-1==0,1,if(r-1==1,2*c-1, sum(j=0,c-1,binomial(r+c-j-2,c-j-1)*(binomial(2*r-2,2*j)-2*(r-1)*binomial(r-3,j-1)))))))); polcoeff(subst(Ser(A[n+1]),x,x/(1+x))/(1+x),k)

Formula

Main diagonal equals A008353: 2^(n-1)*(2^n-n) for n>1.
O.g.f. : rational function N(x,z)/D(x,z), where N(x,z) = 1 - 3*(1 + 2*x)*z + (3 + 8*x + 8*x^2)*z^2 - (1 + 2*x)*(1 - 6*x - 6*x^2)z^3 - 8*x*(1 + x)(1 + 2*x + 2*x^2)*z^4 + 2*x*(1 + x)*(1 + 2*x)*z^5 and D(x,z) = ((1-z)^2 - 4*x*z)*(1 - z*(1 + 2*x))^2. - Peter Bala, Oct 23 2008

A109001 Triangle, read by rows, where g.f. of row n equals the product of (1-x)^n and the g.f. of the coordination sequence for root lattice B_n, for n >= 0.

Original entry on oeis.org

1, 1, 1, 1, 6, 1, 1, 15, 23, 1, 1, 28, 102, 60, 1, 1, 45, 290, 402, 125, 1, 1, 66, 655, 1596, 1167, 226, 1, 1, 91, 1281, 4795, 6155, 2793, 371, 1, 1, 120, 2268, 12040, 23750, 18888, 5852, 568, 1, 1, 153, 3732, 26628, 74574, 91118, 49380, 11124, 825, 1, 1, 190, 5805, 53544, 201810, 350196, 291410, 114600, 19629, 1150, 1
Offset: 0

Views

Author

Paul D. Hanna, Jun 17 2005

Keywords

Comments

Compare to triangle A108558, where row n equals the (n+1)-th differences of the crystal ball sequence for D_n lattice.

Examples

			G.f.s of initial rows of square array A108998 are:
  (1),
  (1 + x)/(1-x),
  (1 + 6*x + x^2)/(1-x)^2;
  (1 + 15*x + 23*x^2 + x^3)/(1-x)^3;
  (1 + 28*x + 102*x^2 + 60*x^3 + x^4)/(1-x)^4.
Triangle begins:
  1;
  1,   1;
  1,   6,    1;
  1,  15,   23,     1;
  1,  28,  102,    60,     1;
  1,  45,  290,   402,   125,     1;
  1,  66,  655,  1596,  1167,   226,     1;
  1,  91, 1281,  4795,  6155,  2793,   371,     1;
  1, 120, 2268, 12040, 23750, 18888,  5852,   568,   1;
  1, 153, 3732, 26628, 74574, 91118, 49380, 11124, 825, 1;
		

Crossrefs

Cf. A108998, A108999, A109000, A022144 (row 2), A022145 (row 3), A022146 (row 4), A022147 (row 5), A022148 (row 6), A022149 (row 7), A022150 (row 8), A022151 (row 9), A022152 (row 10), A022153 (row 11), A022154 (row 12).

Programs

  • GAP
    Flat(List([0..10],n->List([0..n],k->Binomial(2*n+1,2*k)-2*n*Binomial(n-1,k-1)))); # Muniru A Asiru, Dec 14 2018
  • Mathematica
    T[n_, k_] := Binomial[2n+1, 2k] - 2n * Binomial[n-1, k-1]; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Amiram Eldar, Dec 14 2018 *)
  • PARI
    T(n,k)=binomial(2*n+1,2*k)-2*n*binomial(n-1,k-1)
    

Formula

T(n, k) = C(2*n+1, 2*k) - 2*n*C(n-1, k-1).
Row sums are 2^n*(2^n - n) for n >= 0.
G.f. for coordination sequence of B_n lattice: ((Sum_{i=0..n} binomial(2*n+1, 2*i)*z^i) - 2*n*z*(1+z)^(n-1))/(1-z)^n. [Bacher et al.]

A204621 Triangle read by rows: coordinator triangle for lattice A*_n.

Original entry on oeis.org

1, 1, 1, 1, 4, 1, 1, 5, 5, 1, 1, 6, 16, 6, 1, 1, 7, 22, 22, 7, 1, 1, 8, 29, 64, 29, 8, 1, 1, 9, 37, 93, 93, 37, 9, 1, 1, 10, 46, 130, 256, 130, 46, 10, 1, 1, 11, 56, 176, 386, 386, 176, 56, 11, 1, 1, 12, 67, 232, 562, 1024, 562, 232, 67, 12, 1
Offset: 0

Views

Author

N. J. A. Sloane, Jan 17 2012

Keywords

Examples

			Triangle begins:
                   1
                1    1
              1    4    1
            1    5    5    1
          1    6    16    6    1
        1    7    22    22    7    1
      1    8    29    64    29    8    1
    1    9    37    93    93    37    9    1
  1    10    46    130    256    130    46    10    1
1     11    56    176    386    386    176    56    11    1
...
		

Crossrefs

The triangle for Z^n is A007318, A_n is A008459, D_n is A108558, D*_n is A008518.
T(2n,n) gives A000302.

Programs

  • GAP
    Flat(List([0..10],n->List([0..n],k->Sum([0..Minimum(k,n-k)],i->Binomial(n+1,i))))); # Muniru A Asiru, Dec 14 2018
  • Mathematica
    T[n_, k_] := Sum[Binomial[n+1, i] , {i, 0, Min[k, n-k]}]; Table[T[n,k], {n,0,10}, {k,0,n}] // Flatten (* Amiram Eldar, Dec 14 2018 *)

Formula

T(n, k) = Sum_{i=0..min(k,n-k)} binomial(n+1,i). [Wang and Yu, Theorem 4.1] - Eric M. Schmidt, Dec 07 2017

A305693 a(n) = binomial(4*n, 2*n) - 4*n*binomial(2*n-2, n-1).

Original entry on oeis.org

1, 2, 54, 852, 12550, 183356, 2698108, 40090728, 600970566, 9074671980, 137844584020, 2104090834456, 32247569822364, 495918392331992, 7648690018326840, 118264579157865424, 1832624131015069254, 28453041434367110220, 442512540108817131364
Offset: 0

Views

Author

Seiichi Manyama, Jun 08 2018

Keywords

Crossrefs

Programs

  • PARI
    {a(n) = binomial(4*n, 2*n)-4*n*binomial(2*n-2, n-1)}

Formula

a(n) = A000984(2*n) - 4*n*A000984(n-1) for n > 0.
G.f.: sqrt(1 + sqrt(1 - 16*x))/sqrt(2*(1 - 16*x)) - 4*x*(1 - 2*x)/(1 - 4*x)^(3/2). - Ilya Gutkovskiy, Jun 08 2018
Showing 1-8 of 8 results.