|
%I #21 Apr 21 2023 11:23:34
%S 1,1,1,1,1,1,1,1,2,1,1,1,3,7,1,1,1,4,13,34,1,1,1,5,19,85,216,1,1,1,6,
%T 25,154,701,1696,1,1,1,7,31,241,1456,7261,15898,1,1,1,8,37,346,2481,
%U 18136,89125,173468,1,1,1,9,43,469,3776,35761,260002,1277865,2161036,1
%N Square array T(n,k), n >= 0, k >= 0, read by antidiagonals downwards, where T(n,k) = n! * Sum_{j=0..floor(n/2)} (k/2)^j * (j+1)^(n-j-1) / (j! * (n-2*j)!).
%H Seiichi Manyama, <a href="/A362377/b362377.txt">Antidiagonals n = 0..139, flattened</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/LambertW-Function.html">Lambert W-Function</a>.
%F E.g.f. A_k(x) of column k satisfies A_k(x) = exp(x + k*x^2/2 * A_k(x)).
%F A_k(x) = exp(x - LambertW(-k*x^2/2 * exp(x))).
%F A_k(x) = -2 * LambertW(-k*x^2/2 * exp(x))/(k*x^2) for k > 0.
%e Square array begins:
%e 1, 1, 1, 1, 1, 1, 1, ...
%e 1, 1, 1, 1, 1, 1, 1, ...
%e 1, 2, 3, 4, 5, 6, 7, ...
%e 1, 7, 13, 19, 25, 31, 37, ...
%e 1, 34, 85, 154, 241, 346, 469, ...
%e 1, 216, 701, 1456, 2481, 3776, 5341, ...
%e 1, 1696, 7261, 18136, 35761, 61576, 97021, ...
%o (PARI) T(n, k) = n! * sum(j=0, n\2, (k/2)^j*(j+1)^(n-j-1)/(j!*(n-2*j)!));
%Y Columns k=0..3 give A000012, A143740, A125500, A362380.
%Y Cf. A362378, A362394.
%K nonn,tabl
%O 0,9
%A _Seiichi Manyama_, Apr 20 2023
|
