A346983 -id:A346983 - cp's OEIS frontend

A094417 Generalized ordered Bell numbers Bo(4,n).

1, 4, 36, 484, 8676, 194404, 5227236, 163978084, 5878837476, 237109864804, 10625889182436, 523809809059684, 28168941794178276, 1641079211868751204, 102961115527874385636, 6921180217049667005284, 496267460209336700111076, 37807710659221213027893604

Offset: 0

Ralf Stephan, May 02 2004

nonn

Fourth row of array A094416, which has more information.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.

Cf. A000670, A004123, A032033, A050353, A094418, A094419, A131689.

Cf. A346983, A354242, A365567.

m:=20; R:=LaurentSeriesRing(RationalField(), m); b:=Coefficients(R!(1/(5 - 4*Exp(x)))); [Factorial(n-1)*b[n]: n in [1..m]]; // Bruno Berselli, Mar 17 2014

a:= proc(n) option remember;
      `if`(n=0, 1, 4* add(binomial(n, k) *a(k), k=0..n-1))
    end:
seq(a(n), n=0..20);

max = 16; f[x_] := 1/(5-4*E^x); CoefficientList[Series[f[x], {x, 0, max}], x]*Range[0, max]! (* Jean-François Alcover, Nov 14 2011, after g.f. *)

my(N=25,x='x+O('x^N)); Vec(serlaplace(1/(5 - 4*exp(x)))) \\ Joerg Arndt, Jan 15 2024

def A094416(n,k): return sum(factorial(j)*n^j*stirling_number2(k,j) for j in range(k+1)) # array
def A094417(k): return A094416(4,k)
[A094417(n) for n in range(31)] # G. C. Greubel, Jan 12 2024

E.g.f.: 1/(5 - 4*exp(x)).
a(n) = 4 * A050353(n) for n>0.
a(n) = Sum_{k=0..n} A131689(n,k) * 4^k. - Philippe Deléham, Nov 03 2008
E.g.f.: A(x) with A_n = 4 * Sum_{k=0..n-1} C(n,k) * A_k; A_0 = 1. - Vladimir Kruchinin, Jan 27 2011
G.f.: 2/G(0), where G(k)= 1 + 1/(1 - 8*x*(k+1)/(8*x*(k+1) - 1 + 10*x*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
a(n) = log(5/4)*int {x = 0..inf} (floor(x))^n * (5/4)^(-x) dx. - Peter Bala, Feb 14 2015
a(0) = 1; a(n) = 4 * a(n-1) - 5 * Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 16 2023
From Seiichi Manyama, Jun 01 2025: (Start)
a(n) = (-1)^(n+1)/5 * Li_{-n}(5/4), where Li_{n}(x) is the polylogarithm function.
a(n) = (1/5) * Sum_{k>=0} k^n * (4/5)^k.
a(n) = (4/5) * Sum_{k=0..n} 5^k * (-1)^(n-k) * A131689(n,k) for n > 0. (End)

A346982 Expansion of e.g.f. 1 / (4 - 3 * exp(x))^(1/3).

Original entry on oeis.org

1, 1, 5, 41, 477, 7201, 133685, 2945881, 75145677, 2177900241, 70687244965, 2539879312521, 100086803174077, 4291845333310081, 198954892070938645, 9914294755149067961, 528504758009562261677, 30010032597449931644721, 1808359960001658961070725

Offset: 0

Ilya Gutkovskiy, Aug 09 2021

Stirling transform of A007559.

Seiichi Manyama, Table of n, a(n) for n = 0..373

Cf. A000670, A007559, A305404, A346983, A346984, A346985, A352117, A352118, A352119.

Cf. A032033, A365558.

g:= proc(n) option remember; `if`(n<2, 1, (3*n-2)*g(n-1)) end:
b:= proc(n, m) option remember;
     `if`(n=0, g(m), m*b(n-1, m)+b(n-1, m+1))
    end:
a:= n-> b(n, 0):
seq(a(n), n=0..18);  # Alois P. Heinz, Aug 09 2021

nmax = 18; CoefficientList[Series[1/(4 - 3 Exp[x])^(1/3), {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[StirlingS2[n, k] 3^k Pochhammer[1/3, k], {k, 0, n}], {n, 0, 18}]

a(n) = Sum_{k=0..n} Stirling2(n,k) * A007559(k).
a(n) ~ n! / (Gamma(1/3) * 2^(2/3) * n^(2/3) * log(4/3)^(n + 1/3)). - Vaclav Kotesovec, Aug 14 2021
From  Peter Bala, Aug 22 2023: (Start)
O.g.f. (conjectural): 1/(1 - x/(1 - 4*x/(1 - 4*x/(1 - 8*x/(1 - 7*x/(1 - 12*x/(1 - ... - (3*n-2)*x/(1 - 4*n*x/(1 - ... ))))))))) - a continued fraction of Stieltjes-type (S-fraction).
More generally, it appears that the o.g.f. of the sequence whose e.g.f. is equal to 1/(r+1 - r*exp(s*x))^(m/s) corresponds to the S-fraction 1/(1 - r*m*x/(1 - s*(r+1)*x/(1 - r*(m+s)*x/(1 - 2*s(r+1)*x/(1 - r*(m+2*s)*x/(1 - 3*s(r+1)*x/( 1 - ... ))))))). This is the case r = 3, s = 1, m = 1/3. (End)
a(0) = 1; a(n) = Sum_{k=1..n} (3 - 2*k/n) * binomial(n,k) * a(n-k). - Seiichi Manyama, Sep 09 2023
a(0) = 1; a(n) = a(n-1) - 4*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 16 2023

A346985 Expansion of e.g.f. 1 / (7 - 6 * exp(x))^(1/6).

Original entry on oeis.org

1, 1, 8, 113, 2325, 62896, 2109143, 84403033, 3924963750, 207976793991, 12369246804853, 815880360117978, 59107920881218525, 4665585774576259261, 398534278371999103888, 36627974592437584634573, 3603954453161886215458025, 377983931878997401821759456, 42095013846928585982896180123

Offset: 0

Ilya Gutkovskiy, Aug 09 2021

nonn

Stirling transform of A008542.

In general, for k >= 1, if e.g.f. = 1 / (k + 1 - k*exp(x))^(1/k), then a(n) ~ n! / (Gamma(1/k) * (k+1)^(1/k) * n^(1 - 1/k) * log(1 + 1/k)^(n + 1/k)). - Vaclav Kotesovec, Aug 14 2021

Alois P. Heinz, Table of n, a(n) for n = 0..343

Cf. A000670, A008542, A094419, A305404, A346982, A346983, A346984, A352117, A352118, A352119.

Cf. A094419, A354252, A365555, A365556, A365557.

g:= proc(n) option remember; `if`(n<2, 1, (6*n-5)*g(n-1)) end:
b:= proc(n, m) option remember;
     `if`(n=0, g(m), m*b(n-1, m)+b(n-1, m+1))
    end:
a:= n-> b(n, 0):
seq(a(n), n=0..18);  # Alois P. Heinz, Aug 09 2021

nmax = 18; CoefficientList[Series[1/(7 - 6 Exp[x])^(1/6), {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[StirlingS2[n, k] 6^k Pochhammer[1/6, k], {k, 0, n}], {n, 0, 18}]

a[n]:=if n=0 then 1 else (1/n)*sum(binomial(n,k)*(n+5*k)*a[k],k,0,n-1);
makelist(a[n],n,0,50); /* Tani Akinari, Aug 22 2023 */

a(n) = Sum_{k=0..n} Stirling2(n,k) * A008542(k).
a(n) ~ n! / (Gamma(1/6) * 7^(1/6) * n^(5/6) * log(7/6)^(n + 1/6)). - Vaclav Kotesovec, Aug 14 2021
For n > 0, a(n) = (1/n)*Sum_{k=0..n-1} binomial(n,k)*(n+5*k)*a(k). - Tani Akinari, Aug 22 2023
O.g.f. (conjectural): 1/(1 - x/(1 - 7*x/(1 - 7*x/(1 - 14*x/(1 - 13*x/(1 - 21*x/(1 - ... - (6*n-5)*x/(1 - 7*n*x/(1 - ... ))))))))) - a continued fraction of Stieltjes-type (S-fraction). - Peter Bala, Aug 25 2023
a(0) = 1; a(n) = a(n-1) - 7*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 17 2023

A346984 Expansion of e.g.f. 1 / (6 - 5 * exp(x))^(1/5).

Original entry on oeis.org

1, 1, 7, 85, 1495, 34477, 983983, 33476437, 1322441575, 59492222077, 3002578396255, 168005805229285, 10321907081030167, 690761732852321677, 50015387402165694607, 3895721046926471861365, 324805103526730206129607, 28861947117644330678207389, 2722944810091827410698112959

Offset: 0

Ilya Gutkovskiy, Aug 09 2021

Stirling transform of A008548.

Seiichi Manyama, Table of n, a(n) for n = 0..351

Cf. A000670, A008548, A094418, A305404, A346982, A346983, A346985, A352117, A352118, A352119.

g:= proc(n) option remember; `if`(n<2, 1, (5*n-4)*g(n-1)) end:
b:= proc(n, m) option remember;
     `if`(n=0, g(m), m*b(n-1, m)+b(n-1, m+1))
    end:
a:= n-> b(n, 0):
seq(a(n), n=0..18);  # Alois P. Heinz, Aug 09 2021

nmax = 18; CoefficientList[Series[1/(6 - 5 Exp[x])^(1/5), {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[StirlingS2[n, k] 5^k Pochhammer[1/5, k], {k, 0, n}], {n, 0, 18}]

a(n) = Sum_{k=0..n} Stirling2(n,k) * A008548(k).
a(n) ~ n! / (Gamma(1/5) * 6^(1/5) * n^(4/5) * log(6/5)^(n + 1/5)). - Vaclav Kotesovec, Aug 14 2021
O.g.f. (conjectural): 1/(1 - x/(1 - 6*x/(1 - 6*x/(1 - 12*x/(1 - 11*x/(1 - 18*x/(1 - ... - (5*n-4)*x/(1 - 6*n*x/(1 - ... ))))))))) - a continued fraction of Stieltjes-type. - Peter Bala, Aug 22 2023
a(0) = 1; a(n) = Sum_{k=1..n} (5 - 4*k/n) * binomial(n,k) * a(n-k). - Seiichi Manyama, Sep 09 2023
a(0) = 1; a(n) = a(n-1) - 6*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 16 2023

A347016 Expansion of e.g.f. 1 / (1 + 4 * log(1 - x))^(1/4).

Original entry on oeis.org

1, 1, 6, 62, 916, 17644, 419360, 11859840, 388965600, 14514046560, 607165485120, 28143329181120, 1431690475207680, 79302863940387840, 4751108622148907520, 306118435580577146880, 21107196651940518551040, 1550773243761690603179520, 120947288498720390755353600

Offset: 0

Ilya Gutkovskiy, Aug 10 2021

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..360

Cf. A007696, A007840, A346978, A346983, A347015.

g:= proc(n) option remember; `if`(n<2, 1, (4*n-3)*g(n-1)) end:
a:= n-> add(abs(Stirling1(n, k))*g(k), k=0..n):
seq(a(n), n=0..18);  # Alois P. Heinz, Aug 10 2021

nmax = 18; CoefficientList[Series[1/(1 + 4 Log[1 - x])^(1/4), {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[Abs[StirlingS1[n, k]] 4^k Pochhammer[1/4, k], {k, 0, n}], {n, 0, 18}]

a(n) = Sum_{k=0..n} |Stirling1(n,k)| * A007696(k).
a(n) ~ n! * exp(n/4) / (Gamma(1/4) * 2^(1/2) * n^(3/4) * (exp(1/4) - 1)^(n + 1/4)). - Vaclav Kotesovec, Aug 14 2021
a(0) = 1; a(n) = Sum_{k=1..n} (4 - 3*k/n) * (k-1)! * binomial(n,k) * a(n-k). - Seiichi Manyama, Sep 09 2023

A347021 Expansion of e.g.f. 1 / (1 - 4 * log(1 + x))^(1/4).

Original entry on oeis.org

1, 1, 4, 32, 364, 5444, 100520, 2210760, 56406240, 1637877600, 53327583360, 1924096475520, 76198487927040, 3285955396558080, 153273199794071040, 7689131281851770880, 412809183978447306240, 23616192920003184176640, 1434201753814306170808320

Offset: 0

Ilya Gutkovskiy, Aug 11 2021

nonn

Cf. A006252, A007696, A320343, A346983, A347016, A347020, A347022, A347023.

nmax = 18; CoefficientList[Series[1/(1 - 4 Log[1 + x])^(1/4), {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[StirlingS1[n, k] 4^k Pochhammer[1/4, k], {k, 0, n}], {n, 0, 18}]

a(n) = Sum_{k=0..n} Stirling1(n,k) * A007696(k).
a(n) ~ n! * exp(1/16) / (Gamma(1/4) * 2^(1/2) * n^(3/4) * (exp(1/4) - 1)^(n + 1/4)). - Vaclav Kotesovec, Aug 14 2021
a(0) = 1; a(n) = Sum_{k=1..n} (-1)^(k-1) * (4 - 3*k/n) * (k-1)! * binomial(n,k) * a(n-k). - Seiichi Manyama, Sep 11 2023

A365567 Expansion of e.g.f. 1 / (5 - 4 * exp(x))^(3/4).

Original entry on oeis.org

1, 3, 24, 297, 5001, 106578, 2748399, 83182347, 2890153626, 113364686403, 4954547485149, 238734066994272, 12573018414279501, 718498413957515103, 44278797576715884024, 2927171415480872824197, 206625238881832412874501, 15511299587628626891270178

Offset: 0

Seiichi Manyama, Sep 09 2023

nonn

Cf. A094417, A346983, A354242.

a[n_] := Sum[Product[4*j + 3, {j, 0, k - 1}] * StirlingS2[n, k], {k, 0, n}]; Array[a, 18, 0] (* Amiram Eldar, Sep 11 2023 *)

a(n) = sum(k=0, n, prod(j=0, k-1, 4*j+3)*stirling(n, k, 2));

a(n) = Sum_{k=0..n} (Product_{j=0..k-1} (4*j+3)) * Stirling2(n,k).
a(0) = 1; a(n) = Sum_{k=1..n} (4 - k/n) * binomial(n,k) * a(n-k).
a(n) ~ Gamma(1/4) * n^(n + 1/4) / (5^(3/4) * sqrt(Pi) * exp(n) * log(5/4)^(n + 3/4)). - Vaclav Kotesovec, Nov 11 2023
a(0) = 1; a(n) = 3*a(n-1) - 5*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 16 2023

cp's OEIS Frontend

A094417 Generalized ordered Bell numbers Bo(4,n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

A346982 Expansion of e.g.f. 1 / (4 - 3 * exp(x))^(1/3).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A346985 Expansion of e.g.f. 1 / (7 - 6 * exp(x))^(1/6).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

Formula

A346984 Expansion of e.g.f. 1 / (6 - 5 * exp(x))^(1/5).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A347016 Expansion of e.g.f. 1 / (1 + 4 * log(1 - x))^(1/4).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A347021 Expansion of e.g.f. 1 / (1 - 4 * log(1 + x))^(1/4).

Views

Author

Keywords

Crossrefs

Programs

Mathematica

Formula

A365567 Expansion of e.g.f. 1 / (5 - 4 * exp(x))^(3/4).

Views

Author

Keywords

Crossrefs

Programs

Mathematica

PARI

Formula