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