A337010 a(n) = exp(-1/2) * Sum_{k>=0} (2*k + 3)^n / (2^k * k!).
1, 4, 18, 92, 532, 3440, 24552, 191280, 1612304, 14597952, 141123872, 1449324992, 15743376704, 180203389696, 2166381979264, 27274611880704, 358690234163456, 4916123783848960, 70076765972288000, 1036967662211324928, 15902394743591408640
Offset: 0
Keywords
Programs
-
Mathematica
nmax = 20; CoefficientList[Series[Exp[3 x + (Exp[2 x] - 1)/2], {x, 0, nmax}], x] Range[0, nmax]! a[0] = 1; a[n_] := a[n] = 4 a[n - 1] + Sum[Binomial[n - 1, k - 1] 2^(k - 1) a[n - k], {k, 2, n}]; Table[a[n], {n, 0, 20}]