A293037 E.g.f.: exp(1 + x - exp(x)).

1, 0, -1, -1, 2, 9, 9, -50, -267, -413, 2180, 17731, 50533, -110176, -1966797, -9938669, -8638718, 278475061, 2540956509, 9816860358, -27172288399, -725503033401, -5592543175252, -15823587507881, 168392610536153, 2848115497132448, 20819319685262839

Offset: 0

Author

Seiichi Manyama, Sep 28 2017

Keywords

• sign

Links

• Seiichi Manyama, Table of n, a(n) for n = 0..593

Crossrefs

Column k=1 of A293051.

Column k=1 of A335977.

Cf. A000587 (k=0), this sequence (k=1), A293038 (k=2), A293039 (k=3), A293040 (k=4).

Cf. A000296, A014182, A193683, A193684, A196835.

Programs

Maple

f:= series(exp(1 + x - exp(x)), x= 0, 101): seq(factorial(n) * coeff(f, x, n), n = 0..30); # Muniru A Asiru, Oct 31 2017
# second Maple program:
b:= proc(n, t) option remember; `if`(n=0, 1-2*t,
      add(b(n-j, 1-t)*binomial(n-1, j-1), j=1..n))
    end:
a:= n-> b(n+1, 1):
seq(a(n), n=0..35);  # Alois P. Heinz, Dec 01 2021

Mathematica

m = 26; Range[0, m]! * CoefficientList[Series[Exp[1 + x - Exp[x]], {x, 0, m}], x] (* Amiram Eldar, Jul 06 2020 *)
Table[Sum[Binomial[n, k] * BellB[k, -1], {k, 0, n}], {n, 0, 30}] (* Vaclav Kotesovec, Jul 06 2020 *)

PARI

my(N=40, x='x+O('x^N)); Vec(serlaplace(exp(-exp(x)+1+x)))

PARI

a(n) = if(n==0, 1, -sum(k=0, n-2, binomial(n-1, k)*a(k))); \\ Seiichi Manyama, Aug 02 2021

Formula

a(n) = exp(1) * Sum_{k>=0} (-1)^k*(k + 1)^n/k!. - Ilya Gutkovskiy, Jun 13 2019

a(n) = Sum_{k=0..n} binomial(n,k) * Bell(k, -1). - Vaclav Kotesovec, Jul 06 2020

a(0) = 1; a(n) = - Sum_{k=0..n-2} binomial(n-1,k) * a(k). - Seiichi Manyama, Aug 02 2021