|
|
A367785
|
|
Expansion of e.g.f. exp(exp(3*x) - x - 1).
|
|
3
|
|
|
1, 2, 13, 89, 772, 7745, 87949, 1109288, 15332539, 229840361, 3706130914, 63857565095, 1169261937973, 22646779177898, 462143532144937, 9902312863237637, 222119823632283628, 5202170552214520637, 126914730275907871201, 3218552632981994910248, 84686139239808135094879, 2307953474037054591248501
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
a(n) = exp(-1) * Sum_{k>=0} (3*k-1)^n / k!.
a(0) = 1; a(n) = -a(n-1) + Sum_{k=1..n} binomial(n-1,k-1) * 3^k * a(n-k).
a(n) = Sum_{k=0..n} (-1)^(n-k) * binomial(n,k) * 3^k * Bell(k).
|
|
MATHEMATICA
|
nmax = 21; CoefficientList[Series[Exp[Exp[3 x] - x - 1], {x, 0, nmax}], x] Range[0, nmax]!
a[0] = 1; a[n_] := a[n] = -a[n - 1] + Sum[Binomial[n - 1, k - 1] 3^k a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 21}]
Table[Sum[(-1)^(n - k) Binomial[n, k] 3^k BellB[k], {k, 0, n}], {n, 0, 21}]
|
|
PROG
|
(PARI) my(x='x+O('x^30)); Vec(serlaplace(exp(exp(3*x) - x - 1))) \\ Michel Marcus, Nov 30 2023
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|