A360176 - OEIS

%I #7 Jan 29 2023 21:02:26

%S 1,0,1,0,-5,1,0,37,-15,1,0,-393,223,-30,1,0,5481,-3815,745,-50,1,0,

%T -95053,76051,-18870,1865,-75,1,0,1975821,-1749811,514381,-65730,3920,

%U -105,1,0,-47939601,45876335,-15316854,2358181,-183610,7322,-140,1

%N Triangle read by rows. T(n, k) = Sum_{j=k..n} binomial(n, j) * (-j)^(n - j) * (-1)^(j - k)* A360177(j, k).

%F E.g.f. of column k: (1 - exp(-LambertW(x*exp(-x))))^k / k!.

%e Triangle T(n, k) starts:

%e [0] 1;

%e [1] 0,         1;

%e [2] 0,        -5,        1;

%e [3] 0,        37,      -15,         1;

%e [4] 0,      -393,      223,       -30,       1;

%e [5] 0,      5481,    -3815,       745,     -50,       1;

%e [6] 0,    -95053,    76051,    -18870,    1865,     -75,    1;

%e [7] 0,   1975821, -1749811,    514381,  -65730,    3920, -105,    1;

%e [8] 0, -47939601, 45876335, -15316854, 2358181, -183610, 7322, -140, 1;

%p T := (n, k) -> add(binomial(n, j) * (-j)^(n - j) * (-1)^(j - k) * A360177(j, k), j = k..n): for n from 0 to 9 do seq(T(n, k), k = 0..n) od;

%p # Alternative:

%p egf := k -> (1 - exp(-LambertW(x*exp(-x))))^k / k!:

%p ser := k -> series(egf(k), x, 22): T := (n, k) -> n!*coeff(ser(k), x, n):

%p for n from 0 to 8 do seq(T(n, k), k = 0..n) od;

%Y Cf. A360177, A273954 (column 1), A028895 (subdiagonal).

%K sign,tabl

%O 0,5

%A _Peter Luschny_, Jan 28 2023