STATUS
proposed
approved
proposed
approved
editing
proposed
allocated for Peter Luschny
a(n) = Sum_{k=0..n} (-1)^(n-k)*(n!/k!) * [x^n] (1 - exp(-LambertW(x*exp(-x))))^k.
1, 1, 6, 53, 647, 10092, 191915, 4309769, 111682044, 3281731611, 107860953795, 3921762633846, 156322429050397, 6779458454252941, 317841794915501862, 16020304439710056785, 863955306007083830051, 49641711131738762890764, 3027776406780183894833791, 195382900651186641677702197
0,3
egf := k -> (1 - exp(-LambertW(x*exp(-x))))^k / k!:
ser := k -> series(egf(k), x, 22):
T := (n, k) -> (-1)^(n-k)*n!*coeff(ser(k), x, n):
seq(add(T(n, k), k = 0..n), n = 0..19);
allocated
nonn
Peter Luschny, Jan 29 2023
approved
editing
allocated for Peter Luschny
allocated
approved