CoefficientList[Series[Exp[2*x]/(1+LambertW[-x]), {x, 0, 20}], x]*Range[0, 20]!
Table[1 + Sum[Binomial[n, m]*(1 + Sum[Binomial[m, k]*k^k, {k, 1, m}]), {m, 1, n}], {n, 0, 20}]
Table[2^n + Sum[Binomial[n, k]*2^(n-k)*k^k, {k, 1, n}], {n, 0, 20}] (*
Vaclav Kotesovec, Oct 28 2016 *)