login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A145904
Square array read by antidiagonals: Hilbert transform of the Narayana numbers A001263.
5
1, 1, 1, 1, 3, 1, 1, 6, 5, 1, 1, 10, 16, 7, 1, 1, 15, 40, 31, 9, 1, 1, 21, 85, 105, 51, 11, 1, 1, 28, 161, 295, 219, 76, 13, 1, 1, 36, 280, 721, 771, 396, 106, 15, 1, 1, 45, 456, 1582, 2331, 1681, 650, 141, 17, 1
OFFSET
0,5
COMMENTS
Refer to A145905 for the definition of the Hilbert transform of a lower triangular array. For the Hilbert transform of A008459, the array of type B Narayana numbers, see A108625.
This seems to be a duplicate of A273350. - Alois P. Heinz, Jun 04 2016. This could probably be proved by showing that the g.f.s are the same. - N. J. A. Sloane, Jul 02 2016
FORMULA
Row n generating function: 1/(n+1) * 1/(1-x) * Jacobi_P(n,1,1,(1+x)/(1-x)) = N_n(x)/(1-x)^n where N_n(x) denotes the shifted Narayana polynomial N_n(x) = sum{k = 1..n} A001263(k)*x^(k-1) of degree n-1.
Conjectural column n generating function: N_n(x^2)/(1-x)^(2n+1).
The entries in row n are given by the values of a polynomial function p_n(x) at x = 0,1,2,... . The first few are p_1(x) = 2x + 1, p_2(x) = (5x^2 + 5x + 2)/2, p_3(x) = (2x + 1)*(7x^2 + 7x + 6)/6 and p_4(x) = (7x^4 + 14x^3 + 21x^2 + 14x + 4)/4. These polynomials appear to have their zeros on the line Re x = -1/2; that is, the polynomials p_n(-x) appear to satisfy a Riemann hypothesis. The corresponding result for A108625 is true (see A142995 for details).
Contribution from Paul Barry, Jan 06 2009: (Start)
The g.f. for the corresponding number triangle is:
1/(1-x-xy-x^2y/(1-x-x^2y/(1-x-xy-x^2y/(1-x-x^2y/(1-x-xy-x^2y.... (a continued fraction). (End)
This g.f. satisfies x^2*y*g^2 - (1-x-x*y)*g + 1 = 0. - R. J. Mathar, Jun 16 2016
G.f.: ((y-1)*sqrt(((x^2+2*x+1)*y-x^2+2*x-1)/(y-1))+(-x-1)*y-x+1)/(2*x^2*y). - Vladimir Kruchinin, Jan 15 2018
T(n,m) = 1/(n+1)*Sum_{i=0..m+1} C(n+1,i-1)*C(n+1,i)*C(n+m-i+1,m+1-i). - Vladimir Kruchinin, Jan 15 2018
EXAMPLE
The array begins
n\k|..0.....1.....2.....3.....4.....5
=====================================
0..|..1.....1.....1.....1.....1.....1
1..|..1.....3.....5.....7.....9....11
2..|..1.....6....16....31....51....76
3..|..1....10....40...105...219...396
4..|..1....15....85...295...771..1681
5..|..1....21...161...721..2331..6083
...
Row 2: (1 + 3x + x^2)/(1 - x)^3 = 1 + 6x + 16x^2 + 31x^3 + ... .
Row 3: (1 + 6x + 6x^2 + x^3)/(1 - x)^4 = 1 + 10x + 40x^2 + 105x^3 + ... .
MATHEMATICA
Table[1/(# + 1)*Sum[Binomial[# + 1, i - 1] Binomial[# + 1, i] Binomial[# + k - i + 1, k + 1 - i], {i, 0, k + 1}] &[m - k], {m, 0, 9}, {k, 0, m}] // Flatten (* Michael De Vlieger, Jan 15 2018 *)
PROG
(Maxima)
taylor(((y-1)*sqrt(((x^2+2*x+1)*y-x^2+2*x-1)/(y-1))+(-x-1)*y-x+1)/(2*x^2*y), x, 0, 10, y, 0, 10);
T(n, m, k):=1/(n+1)*sum(binomial(n+1, i-1)*binomial(n+1, i)*binomial(n+m-i+1, m+1-i), i, 0, m+1); /* Vladimir Kruchinin, Jan 15 2018 */
CROSSREFS
Cf. A001263, A005891 (row 2), A063490 (row 3), A108625 (Hilbert transform of h-vectors of type B associahedra).
Cf. also A273350.
Sequence in context: A129818 A123970 A055898 * A273350 A328083 A370174
KEYWORD
easy,nonn,tabl
AUTHOR
Peter Bala, Oct 31 2008
STATUS
approved