A346843 -id:A346843 - cp's OEIS frontend

A346842 E.g.f.: exp(exp(x) - 1) * (exp(x) - 1)^3 / 3!.

1, 10, 75, 520, 3556, 24626, 174805, 1279240, 9677151, 75750752, 613656836, 5142797660, 44557627661, 398786697398, 3683575764083, 35084121263136, 344242894197456, 3476490965903174, 36104281709286841, 385257741260565844, 4220537246457019687, 47432055430482106880

Offset: 3

Ilya Gutkovskiy, Aug 05 2021

nonn

Cf. A000110, A000292, A000392, A003128, A005493, A049020, A189924, A346843, A346844.

b:= proc(n, m) option remember;
     `if`(n=0, binomial(m, 3), m*b(n-1, m)+b(n-1, m+1))
    end:
a:= n-> b(n, 0):
seq(a(n), n=3..24);  # Alois P. Heinz, Aug 05 2021

nmax = 24; CoefficientList[Series[Exp[Exp[x] - 1] (Exp[x] - 1)^3/3!, {x, 0, nmax}], x] Range[0, nmax]! // Drop[#, 3] &
Table[Sum[StirlingS2[n, k] Binomial[k, 3], {k, 0, n}], {n, 3, 24}]
Table[Sum[Binomial[n, k] StirlingS2[k, 3] BellB[n - k], {k, 0, n}], {n, 3, 24}]
Table[(BellB[n+3] - 6*BellB[n+2] + 8*BellB[n+1] - BellB[n])/6, {n, 3, 24}] (* Vaclav Kotesovec, Aug 06 2021 *)

my(x='x+O('x^25)); Vec(serlaplace(exp(exp(x)-1)*(exp(x)-1)^3/3!)) \\ Michel Marcus, Aug 06 2021

a(n) = Sum_{k=0..n} Stirling2(n,k) * binomial(k,3).
a(n) = Sum_{k=0..n} binomial(n,k) * Stirling2(k,3) * Bell(n-k).
a(n) = (Bell(n+3) - 6*Bell(n+2) + 8*Bell(n+1) - Bell(n))/6. - Vaclav Kotesovec, Aug 06 2021
a(n) ~ exp(-1 - n + n/LambertW(n)) * (n - LambertW(n))^3 * n^n / (6 * sqrt(1 + LambertW(n)) * LambertW(n)^(n+3)). - Vaclav Kotesovec, Jun 28 2022

A346844 E.g.f.: exp(exp(x) - 1) * (exp(x) - 1)^5 / 5!.

Original entry on oeis.org

1, 21, 287, 3290, 34671, 350889, 3492511, 34669734, 346231886, 3497726232, 35872743270, 374387203190, 3982122624117, 43207791878715, 478532965417529, 5411213661200830, 62482405993313229, 736696756305382411, 8868148033487285103, 108969560832001750716

Offset: 5

Ilya Gutkovskiy, Aug 05 2021

nonn

Cf. A000110, A000389, A000481, A003128, A005493, A049020, A094816, A346842, A346843.

b:= proc(n, m) option remember;
     `if`(n=0, binomial(m, 5), m*b(n-1, m)+b(n-1, m+1))
    end:
a:= n-> b(n, 0):
seq(a(n), n=5..24);  # Alois P. Heinz, Aug 05 2021

nmax = 24; CoefficientList[Series[Exp[Exp[x] - 1] (Exp[x] - 1)^5/5!, {x, 0, nmax}], x] Range[0, nmax]! // Drop[#, 5] &
Table[Sum[StirlingS2[n, k] Binomial[k, 5], {k, 0, n}], {n, 5, 24}]
Table[Sum[Binomial[n, k] StirlingS2[k, 5] BellB[n - k], {k, 0, n}], {n, 5, 24}]
Table[(-BellB[n] + 89*BellB[n+1] - 145*BellB[n+2] + 75*BellB[n+3] - 15*BellB[n+4] + BellB[n+5])/120, {n, 5, 24}] (* Vaclav Kotesovec, Aug 06 2021 *)

my(x='x+O('x^25)); Vec(serlaplace(exp(exp(x)-1)*(exp(x)-1)^5/5!)) \\ Michel Marcus, Aug 06 2021

a(n) = Sum_{k=0..n} Stirling2(n,k) * binomial(k,5).
a(n) = Sum_{k=0..n} binomial(n,k) * Stirling2(k,5) * Bell(n-k).
a(n) = (-Bell(n) + 89*Bell(n+1) - 145*Bell(n+2) + 75*Bell(n+3) - 15*Bell(n+4) + Bell(n+5))/120. - Vaclav Kotesovec, Aug 06 2021

cp's OEIS Frontend

A346842 E.g.f.: exp(exp(x) - 1) * (exp(x) - 1)^3 / 3!.

Views

Author

Keywords

Crossrefs

Programs

Maple

Mathematica

PARI

Formula