cp's OEIS Frontend

Limit a(n+1)/a(n) ~ 11.08...

Limit is equal to A001622^5 = A244593 = 11.090169943749474241... - Vaclav Kotesovec, Aug 31 2025

Limit a(n+1)/a(n) ~ 3.3829757679..., real root of cubic (1+x+3*x^2-x^3).

Compare to the corresponding array A108553 of crystal ball sequences for D_n lattice.

From Peter Bala, Jul 18 2008: (Start)

Row reverse of A099608.

This array has a remarkable relationship with the constant zeta(2). The row, column and diagonal entries of the array occur in series acceleration formulas for zeta(2).

For the entries in row n we have zeta(2) = 2*(1 - 1/2^2 + 1/3^2 - ... + (-1)^(n+1)/n^2) + (-1)^n*Sum_{k >= 1} 1/(k^2*T(n,k-1)*T(n,k)). For example, n = 4 gives zeta(2) = 2*(1 - 1/4 + 1/9 - 1/16) + 1/(1*21) + 1/(4*21*131) + 1/(9*131*471) + ... . See A142995 for further details.

For the entries in column k we have zeta(2) = (1 + 1/4 + 1/9 + ... + 1/k^2) + 2*Sum_{n >= 1} (-1)^(n+1)/(n^2*T(n-1,k)*T(n,k)). For example, k = 4 gives zeta(2) = (1 + 1/4 + 1/9 + 1/16) + 2*(1/(1*9) - 1/(4*9*61) + 1/(9*61*309) - ... ). See A142999 for further details.

Also, as consequence of Apery's proof of the irrationality of zeta(2), we have a series acceleration formula along the main diagonal of the table: zeta(2) = 5 * Sum_{n >= 1} (-1)^(n+1)/(n^2*T(n,n)*T(n-1,n-1)) = 5*(1/3 - 1/(2^2*3*19) + 1/(3^2*19*147) - ...).

There also appear to be series acceleration results along other diagonals. For example, for the main subdiagonal, calculation supports the result zeta(2) = 2 - Sum_{n >= 1} (-1)^(n+1)*(n^2+(2*n+1)^2)/(n^2*(n+1)^2*T(n,n-1)*T(n+1,n)) = 2 - 10/(2^2*7) + 29/(6^2*7*55) - 58/(12^2*55*471) + ..., while for the main superdiagonal we appear to have zeta(2) = 1 + Sum_{n >= 1} (-1)^(n+1)*((n+1)^2 + (2*n+1)^2)/(n^2*(n+1)^2*T(n-1,n)*T(n,n+1)) = 1 + 13/(2^2*5) - 34/(6^2*5*37) + 65/(12^2*37*309) - ... .

Similar series acceleration results hold for Apery's constant zeta(3) involving the crystal ball sequences for the product lattices A_n x A_n; see A143007 for further details. Similar results also hold between the constant log(2) and the crystal ball sequences of the hypercubic lattices A_1 x...x A_1 and between log(2) and the crystal ball sequences for lattices of type C_n ; see A008288 and A142992 respectively for further details. (End)

This array is the Hilbert transform of triangle A008459 (see A145905 for the definition of the Hilbert transform). - Peter Bala, Oct 28 2008

If each row started with an initial 0 (i.e., a(n,k) = (n+1)^k - n^k) then each row would be the binomial transform of the preceding row. - Henry Bottomley, May 31 2001

a(n-1, k-1) is the number of ordered k-tuples of positive integers such that the largest of these integers is n. - Alford Arnold, Sep 07 2005

From Alford Arnold, Jul 21 2006: (Start)

The sequences in A047969 can also be calculated using the Eulerian Array (A008292) and Pascal's Triangle (A007318) as illustrated below: (cf. A101095).

1 1 1 1 1 1

1 1 1 1 1 1

-----------------------------------------

1 2 3 4 5 6

1 2 3 4 5

1 3 5 7 9 11

-----------------------------------------

1 3 6 10 15 21

4 12 24 40 60

1 3 6 10

1 7 19 37 61 91

-----------------------------------------

1 4 10 20 35 56

11 44 110 220 385

11 44 110 220

1 4 10

1 15 65 175 369 671

----------------------------------------- (End)

From Peter Bala, Oct 26 2008: (Start)

The above remarks of Alford Arnold may be summarized by saying that (the transpose of) this array is the Hilbert transform of the triangle of Eulerian numbers A008292 (see A145905 for the definition of the Hilbert transform). In this context, A008292 is best viewed as the array of h-vectors of permutohedra of type A. See A108553 for the Hilbert transform of the array of h-vectors of type D permutohedra. Compare this array with A009998.

The polynomials n^k - (n-1)^k, k = 1,2,3,..., which give the nonzero entries in the columns of this array, satisfy a Riemann hypothesis: their zeros lie on the vertical line Re s = 1/2 in the complex plane. See A019538 for the connection between the polynomials n^k - (n-1)^k and the Stirling polynomials of the simplicial complexes dual to the type A permutohedra.

(End)

Empirical: (n+1)^(k+1) - n^(k+1) is the number of first differences of length k+1 arrays of numbers in 0..n, k > 0. - R. H. Hardin, Jun 30 2013

a(n-1, k-1) is the number of bargraphs of width k and height n. Examples: a(1,2) = 7 because we have [1,1,2], [1,2,1], [2,1,1], [1,2,2], [2,1,2], [2,2,1], and [2,2,2]; a(2,1) = 5 because we have [1,3], [2,3], [3,1], [3,2], and [3,3] (bargraphs are given as compositions). This comment is equivalent to A. Arnold's Sep 2005 comment. - Emeric Deutsch, Jan 30 2017

The lattice C_n consists of all integer lattice points v = (x_1,...,x_n) in Z^n such that the sum x_1 + ... + x_n is even. Let ||v|| = 1/2 * Sum_{i = 1..n} |x_i|; this defines a norm on C_n. The k-th term of the crystal ball sequence of C_n gives the number of lattice points v in C_n with ||v|| <= k [Bacher et al.]. The case n = 2 is illustrated in the Example section below.

This array has a remarkable relationship with the constant log(2). The row, column and (conjecturally) the diagonal entries of the array occur in series acceleration formulas for log(2) (see the Formula section below for some examples).

See A103884 for the table of coordination sequences of the C_n lattices. For the crystal ball sequences for the A_n and D_n lattices see A108625 and A108553 respectively. For the crystal ball sequences for the product lattices A_1 x ... x A_1(n copies) and A_n x A_n see A008288 and A143007 respectively.

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.

Row n equals the (n+1)-th differences of row n of the square array A108553. G.f. of row n equals: (1-x)^(n+1)*CBD_n(x), where CBD_n denotes the g.f. of the crystal ball sequence for D_n lattice.

From Peter Bala, Oct 23 2008: (Start)

Let D_n be the root lattice generated as a monoid by {+-e_i +- e_j: 1 <= i not equal to j <= n}. Let P(D_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(D_n) [Ardila et al.]. See A108556 for the corresponding array of f-vectors for these type D_n polytopes. See A008459 for the array of h-vectors for type A_n polytopes and A086645 for the array of h-vectors associated with type C_n polytopes.

The Hilbert transform of this array (as defined in A145905) equals A108553.

(End)

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

Compare with A108553, where row n equals the crystal ball sequence for D_n lattice.

Compare to diagonal A108554 of square array A108553, in which row n equals the crystal ball sequence for D_n lattice.