OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..451
FORMULA
a(n) = n! * Sum_{k=0..n-1} (-3)^k / ((n-k) * k!).
a(0) = 0, a(1) = 1, a(n) = (n-4) * a(n-1) + 3 * (n-1) * a(n-2) + (-3)^(n-1). - Seiichi Manyama, May 27 2022
a(n) ~ exp(-3) * (n-1)!. - Vaclav Kotesovec, Jun 08 2022
MATHEMATICA
nmax = 24; CoefficientList[Series[-Log[1 - x] Exp[-3 x], {x, 0, nmax}], x] Range[0, nmax]!
Table[n! Sum[(-3)^k/((n - k) k!), {k, 0, n - 1}], {n, 0, 24}]
PROG
(PARI) a_vector(n) = my(v=vector(n+1)); v[1]=0; v[2]=1; for(i=2, n, v[i+1]=(i-4)*v[i]+3*(i-1)*v[i-1]+(-3)^(i-1)); v; \\ Seiichi Manyama, May 27 2022
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Jul 15 2021
STATUS
approved