OFFSET
0,4
FORMULA
a(n) = n! * Sum_{k=0..n-1} (-3)^(n-1-k) / ((n-k) * k!).
a(0) = 0, a(1) = 1, a(n) = (-3 * n + 4) * a(n-1) + 3 * (n-1) * a(n-2) + 1.
a(n) ~ -(-1)^n * (n-1)! * 3^(n-1) / exp(1/3). - Vaclav Kotesovec, Jun 08 2022
PROG
(PARI) my(N=20, x='x+O('x^N)); concat(0, Vec(serlaplace(log(1+3*x)*exp(x)/3)))
(PARI) a(n) = n!*sum(k=0, n-1, (-3)^(n-1-k)/((n-k)*k!));
(PARI) a_vector(n) = my(v=vector(n+1)); v[1]=0; v[2]=1; for(i=2, n, v[i+1]=(-3*i+4)*v[i]+3*(i-1)*v[i-1]+1); v;
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, May 27 2022
STATUS
approved