|
%I #56 Apr 07 2023 10:52:19
%S 1,1,3,1,15,10,1,105,105,25,1,945,1260,490,56,1,10395,17325,9450,1918,
%T 119,1,135135,270270,190575,56980,6825,246,1,2027025,4729725,4099095,
%U 1636635,302995,22935,501,1,34459425,91891800,94594500,47507460,12122110,1487200
%N Triangle of Ward numbers T(n,k) (n>=0, k=0 if n=0, otherwise 0 <= k <= n-1) read by rows.
%C It appears that the sum of row(n) is A000311(n+1). - _Michel Marcus_, Feb 07 2013
%C Conjecture on row sums was proved in the first paragraph of the formula section of the reverse matrix A134991 in 2008 (e.g.f. evaluated at t=1). - _Tom Copeland_, Jan 03 2016
%D Charles Jordan, Calculus of Finite Differences, Chelsea 1950, p. 172, Table C_{m, i}.
%H L. Carlitz, <a href="http://projecteuclid.org/euclid.dmj/1077377459">The coefficients in an asymptotic expansion and certain related numbers</a>, Duke Math. J., vol 35 (1968), p. 83-90. See page 85.
%H José L. Cereceda, <a href="https://arxiv.org/abs/2001.03208">Figurate numbers and sums of powers of integers</a>, arXiv:2001.03208 [math.NT], 2020. See Table 8 p. 11.
%H Lane Clark, <a href="http://dx.doi.org/10.1016/S0012-365X(99)00019-9">Asymptotic normality of the Ward numbers</a>, Discrete Math. 203 (1999), no. 1-3, 41-48. MR1696232 (2000d:11101)
%H T. Copeland, <a href="http://tcjpn.wordpress.com/2015/12/21/generators-inversion-and-matrix-binomial-and-integral-transforms/">Generators, Inversion, and Matrix, Binomial, and Integral Transforms</a>
%H Andrew Elvey Price and Alan D. Sokal, <a href="https://arxiv.org/abs/2001.01468">Phylogenetic trees, augmented perfect matchings, and a Thron-type continued fraction (T-fraction) for the Ward polynomials</a>, arXiv:2001.01468 [math.CO], 2020.
%F T(n, k) = Sum_{m = 0..n-k} (-1)^(n - k + m)*C(2*n - k, n + m)*Stirling2(n + m, m). - _Peter Luschny_, Feb 19 2021
%e Triangle begins:
%e 1
%e 1
%e 3 1
%e 15 10 1
%e 105 105 25 1
%e 945 1260 490 56 1
%e 10395 17325 9450 1918 119 1 ...
%p A181996 := (n,k) -> add((-1)^(n - k + m)*binomial(2*n - k, n + m)*Stirling2(n + m, m), m = 0..n-k):
%p seq(seq(A181996(n, k), k = 0..n-1+0^n), n=0..8); # _Peter Luschny_, Feb 19 2021
%o (PARI) T(n,k) = {if (n == 0, return(1)); if (k == 0, return (prod(x=2,n, 2*x-1))); if (k == n, return (0)); return((2*n-1-k)*T(n-1,k) + (n-k)*T(n-1, k-1));} \\ _Michel Marcus_, Feb 07 2013
%Y See A134991, which is the mirror image and is the main entry for this triangle, for further information.
%Y Cf. A000311.
%K nonn,tabf
%O 0,3
%A _N. J. A. Sloane_, Apr 05 2012
%E More terms from _Michel Marcus_, Feb 07 2013
|
