A094419 -id:A094419 - cp's OEIS frontend

A131689 Triangle of numbers T(n,k) = k!Stirling2(n,k) = A000142(k)A048993(n,k) read by rows, T(n, k) for 0 <= k <= n.

1, 0, 1, 0, 1, 2, 0, 1, 6, 6, 0, 1, 14, 36, 24, 0, 1, 30, 150, 240, 120, 0, 1, 62, 540, 1560, 1800, 720, 0, 1, 126, 1806, 8400, 16800, 15120, 5040, 0, 1, 254, 5796, 40824, 126000, 191520, 141120, 40320, 0, 1, 510, 18150, 186480, 834120, 1905120, 2328480, 1451520, 362880

Offset: 0

Philippe Deléham, Sep 14 2007

Triangle T(n,k), 0 <= k <= n, read by rows given by [0,1,0,2,0,3,0,4,0,5,0,6,0,7,0,...] DELTA [1,1,2,2,3,3,4,4,5,5,6,6,...] where DELTA is the operator defined in A084938; another version of A019538.
See also A019538: version with n > 0 and k > 0. - Philippe Deléham, Nov 03 2008
From Peter Bala, Jul 21 2014: (Start)
T(n,k) gives the number of (k-1)-dimensional faces in the interior of the first barycentric subdivision of the standard (n-1)-dimensional simplex. For example, the barycentric subdivision of the 1-simplex is o--o--o, with 1 interior vertex and 2 interior edges, giving T(2,1) = 1 and T(2,2) = 2.
This triangle is used when calculating the face vectors of the barycentric subdivision of a simplicial complex. Let S be an n-dimensional simplicial complex and write f_k for the number of k-dimensional faces of S, with the usual convention that f_(-1) = 1, so that F := (f_(-1), f_0, f_1,...,f_n) is the f-vector of S. If M(n) denotes the square matrix formed from the first n+1 rows and n+1 columns of the present triangle, then the vector F*M(n) is the f-vector of the first barycentric subdivision of the simplicial complex S (Brenti and Welker, Lemma 2.1). For example, the rows of Pascal's triangle A007318 (but with row and column indexing starting at -1) are the f-vectors for the standard n-simplexes. It follows that A007318*A131689, which equals A028246, is the array of f-vectors of the first barycentric subdivision of standard n-simplexes. (End)
This triangle T(n, k) appears in the o.g.f. G(n, x) = Sum_{m>=0} S(n, m)*x^m with S(n, m) = Sum_{j=0..m} j^n for n >= 1 as G(n, x) = Sum_{k=1..n} (x^k/(1 - x)^(k+2))*T(n, k). See also the Eulerian triangle A008292 with a Mar 31 2017 comment for a rewritten form. For the e.g.f. see A028246 with a Mar 13 2017 comment. - Wolfdieter Lang, Mar 31 2017
T(n,k) = the number of alignments of length k of n strings each of length 1. See Slowinski. An example is given below. Cf. A122193 (alignments of strings of length 2) and A299041 (alignments of strings of length 3). - Peter Bala, Feb 04 2018
The row polynomials R(n,x) are the Fubini polynomials. - Emanuele Munarini, Dec 05 2020
From Gus Wiseman, Feb 18 2022: (Start)
Also the number of patterns of length n with k distinct parts (or with maximum part k), where we define a pattern to be a finite sequence covering an initial interval of positive integers. For example, row n = 3 counts the following patterns:
  (1,1,1)  (1,2,2)  (1,2,3)
           (2,1,2)  (1,3,2)
           (2,2,1)  (2,1,3)
           (1,1,2)  (2,3,1)
           (1,2,1)  (3,1,2)
           (2,1,1)  (3,2,1)
(End)
Regard A048994 as a lower-triangular matrix and divide each term A048994(n,k) by n!, then this is the matrix inverse. Because Sum_{k=0..n} (A048994(n,k) * x^n / n!) = A007318(x,n), Sum_{k=0..n} (A131689(n,k) * A007318(x,k)) = x^n. - Natalia L. Skirrow, Mar 23 2023
T(n,k) is the number of ordered partitions of [n] into k blocks. - Alois P. Heinz, Feb 21 2025

			The triangle T(n,k) begins:
  n\k 0 1    2     3      4       5        6        7        8        9      10 ...
  0:  1
  1:  0 1
  2:  0 1    2
  3:  0 1    6     6
  4:  0 1   14    36     24
  5:  0 1   30   150    240     120
  6:  0 1   62   540   1560    1800      720
  7:  0 1  126  1806   8400   16800    15120     5040
  8:  0 1  254  5796  40824  126000   191520   141120    40320
  9:  0 1  510 18150 186480  834120  1905120  2328480  1451520   362880
  10: 0 1 1022 55980 818520 5103000 16435440 29635200 30240000 16329600 3628800
  ... reformatted and extended. - _Wolfdieter Lang_, Mar 31 2017
From _Peter Bala_, Feb 04 2018: (Start)
T(4,2) = 14 alignments of length 2 of 4 strings of length 1. Examples include
  (i) A -    (ii) A -    (iii) A -
      B -         B -          - B
      C -         - C          - C
      - D         - D          - D
There are C(4,1) = 4 alignments of type (i) with a single gap character - in column 1, C(4,2) = 6 alignments of type (ii) with two gap characters in column 1 and C(4,3) = 4 alignments of type (iii) with three gap characters in column 1, giving a total of 4 + 6 + 4 = 14 alignments. (End)

Vincenzo Librandi, Rows n = 0..100, flattened
Peter Bala, Deformations of the Hadamard product of power series
F. Brenti and V. Welker, f-vectors of barycentric subdivisions, arXiv:math/0606356 [math.CO], Math. Z., 259(4), 849-865, 2008.
M. Dukes and C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
Germain Kreweras, Une dualité élémentaire souvent utile dans les problèmes combinatoires, Mathématiques et Sciences Humaines 3 (1963): 31-41.
Jerry Metzger and Thomas Richards, A Prisoner Problem Variation, Journal of Integer Sequences, Vol. 18 (2015), Article 15.2.7.
Massimo Nocentini, An algebraic and combinatorial study of some infinite sequences of numbers supported by symbolic and logic computation, PhD Thesis, University of Florence, 2019. See Ex. 36.
Mircea Dan Rus, Yet another note on notation, arXiv:2501.08762 [math.HO], 2025. See p. 6.
J. B. Slowinski, The Number of Multiple Alignments, Molecular Phylogenetics and Evolution 10:2 (1998), 264-266. doi:10.1006/mpev.1998.0522
M. Z. Spivey, On Solutions to a General Combinatorial Recurrence, J. Int. Seq. 14 (2011) # 11.9.7.
Wikipedia, Barycentric subdivision
Wikipedia, Simplicial complex
Wikipedia, Simplex
Gus Wiseman, Sequences counting and ranking compositions by the patterns they match or avoid.

Columns k=0..10 are A000007, A000012, A000918, A001117, A000919, A001118, A000920, A135456, A133068, A133360, A133132,
Case m=1 of the polynomials defined in A278073.
Cf. A000142 (diagonal), A000670 (row sums), A000012 (alternating row sums), A210029 (central terms).
Cf. A001286, A037960, A037961, A037962, A037963.
Cf. A008292, A028246 (o.g.f. and e.g.f. of sums of powers).
Cf. A019538, A122193, A299041.
A version for partitions is A116608, or by maximum A008284.
A version for compositions is A235998, or by maximum A048004.
Classes of patterns:
- A000142 = strict
- A005649 = anti-run, complement A069321
- A019536 = necklace
- A032011 = distinct multiplicities
- A060223 = Lyndon
- A226316 = (1,2,3)-avoiding, weakly A052709, complement A335515
- A296975 = aperiodic
- A345194 = alternating, up/down A350354, complement A350252
- A349058 = weakly alternating
- A351200 = distinct runs
- A351292 = distinct run-lengths
Cf. A007318, A097805, A335456, A335457, A344605, A349050.

function T(n, k)
    if k < 0 || k > n return 0 end
    if n == 0 && k == 0 return 1 end
    k*(T(n-1, k-1) + T(n-1, k))
end
for n in 0:7
    println([T(n, k) for k in 0:n])
end
# Peter Luschny, Mar 26 2020

A131689 := (n,k) -> Stirling2(n,k)*k!: # Peter Luschny, Sep 17 2011
# Alternatively:
A131689_row := proc(n) 1/(1-t*(exp(x)-1)); expand(series(%,x,n+1)); n!*coeff(%,x,n); PolynomialTools:-CoefficientList(%,t) end:
for n from 0 to 9 do A131689_row(n) od; # Peter Luschny, Jan 23 2017

t[n_, k_] := k!*StirlingS2[n, k]; Table[t[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 25 2014 *)
T[n_, k_] := If[n <= 0 || k <= 0, Boole[n == 0 && k == 0], Sum[(-1)^(i + k) Binomial[k, i] i^(n + k), {i, 0, k}]]; (* Michael Somos, Jul 08 2018 *)

{T(n, k) = if( n<0, 0, sum(i=0, k, (-1)^(k + i) * binomial(k, i) * i^n))};
/* Michael Somos, Jul 08 2018 */

@cached_function
def F(n): # Fubini polynomial
    R. = PolynomialRing(ZZ)
    if n == 0: return R(1)
    return R(sum(binomial(n, k)*F(n - k)*x for k in (1..n)))
for n in (0..9): print(F(n).list()) # Peter Luschny, May 21 2021

T(n,k) = k*(T(n-1,k-1) + T(n-1,k)) with T(0,0)=1. Sum_{k=0..n} T(n,k)*x^k = (-1)^n*A000629(n), A033999(n), A000007(n), A000670(n), A004123(n+1), A032033(n), A094417(n), A094418(n), A094419(n) for x = -2, -1, 0, 1, 2, 3, 4, 5, 6 respectively. [corrected by Philippe Deléham, Feb 11 2013]
Sum_{k=0..n} T(n,k)*x^(n-k) = A000012(n), A000142(n), A000670(n), A122704(n) for x=-1, 0, 1, 2 respectively. - Philippe Deléham, Oct 09 2007
Sum_{k=0..n} (-1)^k*T(n,k)/(k+1) = Bernoulli numbers A027641(n)/A027642(n). - Peter Luschny, Sep 17 2011
G.f.: F(x,t) = 1 + x*t + (x+x^2)*t^2/2! + (x+6*x^2+6*x^3)*t^3/3! + ... = Sum_{n>=0} R(n,x)*t^n/n!.
The row polynomials R(n,x) satisfy the recursion R(n+1,x) = (x+x^2)*R'(n,x) + x*R(n,x) where ' indicates differentiation with respect to x. - Philippe Deléham, Feb 11 2013
T(n,k) = [t^k] (n! [x^n] (1/(1-t*(exp(x)-1)))). - Peter Luschny, Jan 23 2017
The n-th row polynomial has the form x o x o ... o x (n factors), where o denotes the black diamond multiplication operator of Dukes and White. See also Bala, Example E8. - Peter Bala, Jan 08 2018

A032033 Stirling transform of A032031.

Original entry on oeis.org

1, 3, 21, 219, 3045, 52923, 1103781, 26857659, 746870565, 23365498683, 812198635941, 31055758599099, 1295419975298085, 58538439796931643, 2848763394161128101, 148537065755389540539, 8261178848690959117605, 488177936257344615487803, 30544839926043868901604261

Offset: 0

Christian G. Bower

Also "AIJ" (ordered, indistinct, labeled) transform of 3,3,3,3...

Third row of array A094416 (generalized ordered Bell numbers).

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.
P. Blasiak, K. A. Penson and A. I. Solomon, Dobinski-type relations and the log-normal distribution, arXiv:quant-ph/0303030, J. Phys. A.: Math. Gen 36 (2003) L273.
Christian G. Bower, Transforms (2)
Jacob Sprittulla, On Colored Factorizations, arXiv:2008.09984 [math.CO], 2020.

Cf. A032031, A094416, A094418, A094419.

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

a[n_] := PolyLog[-n, 3/4]/4; a[0] = 1; Table[a[n], {n, 0, 16}] (* Jean-François Alcover, Nov 14 2011 *)
t = 30; Range[0, t]! CoefficientList[Series[1/(4 - 3 Exp[x]), {x, 0, t}], x] (* Vincenzo Librandi, Mar 16 2014 *)

a(n)=ceil(polylog(-n,3/4)/4) \\ Charles R Greathouse IV, Jul 14 2014

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

E.g.f.: 1/(4-3*exp(x)).
a(n) = 3 * A050352(n), n > 0.
a(n) = Sum_{k=0..n} Stirling2(n, k) * (3^k) * k!.
a(n) = (1/4) * Sum_{k>=0} k^n * (3/4)^k. - Karol A. Penson, Jan 25 2002
a(n) = Sum_{k=0..n} A131689(n,k)*3^k. - Philippe Deléham, Nov 03 2008
G.f. A(x)=B(x)/x, where B(x)=x+3*x^2+21*x^3+... = Sum_{n>=1} b(n)*x^n satisfies 4*B(x)-x = 3*B(x/(1-x)), and b(n)=3*Sum_{k=1..n-1} binomial(n-1,k-1)*b(k), b(1)=1. - Vladimir Kruchinin, Jan 27 2011
a(n) = log(4/3)*Integral_{x = 0..inf} (floor(x))^n * (4/3)^(-x) dx. - Peter Bala, Feb 14 2015
a(0) = 1; a(n) = 3 * Sum_{k=1..n} binomial(n,k) * a(n-k). - Ilya Gutkovskiy, Jan 17 2020
a(0) = 1; a(n) = 3 * a(n-1) - 4 * Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 16 2023
a(n) = (3/4) * Sum_{k=0..n} 4^k * (-1)^(n-k) * k! * Stirling2(n,k) for n > 0. - Seiichi Manyama, Jun 01 2025

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

Original entry on oeis.org

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)

A094418 Generalized ordered Bell numbers Bo(5,n).

Original entry on oeis.org

1, 5, 55, 905, 19855, 544505, 17919055, 687978905, 30187495855, 1490155456505, 81732269223055, 4931150091426905, 324557348772511855, 23141780973332248505, 1776997406800302687055, 146197529083891406394905, 12829862285488250150167855, 1196280147496701351115120505

Offset: 0

Ralf Stephan, May 02 2004

nonn

Fifth 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. A094416, A094417, A094419, A094422, A131689.

Cf. A346984, A365568, A365569, A365570.

A094416:= func< n,k | (&+[Factorial(j)*n^j*StirlingSecond(k,j): j in [0..k]]) >;
A094418:= func< k | A094416(5,k) >;
[A094418(n): n in [0..30]]; // G. C. Greubel, Jan 12 2024

t = 30; Range[0, t]! CoefficientList[Series[1/(6 - 5 Exp[x]), {x, 0, t}], x] (* Vincenzo Librandi, Mar 16 2014 *)

my(N=25,x='x+O('x^N)); Vec(serlaplace(1/(6 - 5*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 A094418(k): return A094416(5,k)
[A094418(n) for n in range(31)] # G. C. Greubel, Jan 12 2024

E.g.f.: 1/(6 - 5*exp(x)).
a(n) = Sum_{k=0..n} A131689(n,k) * 5^k. - Philippe Deléham, Nov 03 2008
a(n) ~ n! / (6*(log(6/5))^(n+1)). - Vaclav Kotesovec, Mar 14 2014
a(0) = 1; a(n) = 5 * Sum_{k=1..n} binomial(n,k) * a(n-k). - Ilya Gutkovskiy, Jan 17 2020
a(0) = 1; a(n) = 5 * a(n-1) - 6 * 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)/6 * Li_{-n}(6/5), where Li_{n}(x) is the polylogarithm function.
a(n) = (1/6) * Sum_{k>=0} k^n * (5/6)^k.
a(n) = (5/6) * Sum_{k=0..n} 6^k * (-1)^(n-k) * A131689(n,k) for n > 0. (End)

A094416 Array read by antidiagonals: generalized ordered Bell numbers Bo(r,n).

Original entry on oeis.org

1, 2, 3, 3, 10, 13, 4, 21, 74, 75, 5, 36, 219, 730, 541, 6, 55, 484, 3045, 9002, 4683, 7, 78, 905, 8676, 52923, 133210, 47293, 8, 105, 1518, 19855, 194404, 1103781, 2299754, 545835, 9, 136, 2359, 39390, 544505, 5227236, 26857659, 45375130, 7087261

Offset: 1

Ralf Stephan, May 02 2004

Also, r times the number of (r+1)-level labeled linear rooted trees with n leaves.
"AIJ" (ordered, indistinct, labeled) transform of {r,r,r,...}.
Stirling transform of r^n*n!, i.e. of e.g.f. 1/(1-r*x).
Also, Bo(r,s) is ((x*d/dx)^n)(1/(1+r-r*x)) evaluated at x=1.
r-th ordered Bell polynomial (A019538) evaluated at n.
Bo(r,n) is the n-th moment of a geometric distribution with probability parameter = 1/(r+1).  Here, geometric distribution is the number of failures prior to the first success. - Geoffrey Critzer, Jan 01 2019
Row r (starting at r=0), Bo(r+1, n), is the Akiyama-Tanigawa algorithm applied to the powers of r+1. See Python program below. - Shel Kaphan, May 03 2024

			Array begins as:
  1,  3,   13,    75,     541,     4683,      47293, ...
  2, 10,   74,   730,    9002,   133210,    2299754, ...
  3, 21,  219,  3045,   52923,  1103781,   26857659, ...
  4, 36,  484,  8676,  194404,  5227236,  163978084, ...
  5, 55,  905, 19855,  544505, 17919055,  687978905, ...
  6, 78, 1518, 39390, 1277646, 49729758, 2258233998, ...

G. C. Greubel, Antidiagonals n = 1..50, flattened
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
P. Blasiak, K. A. Penson and A. I. Solomon, Dobinski-type relations and the log-normal distribution, arXiv:quant-ph/0303030, 2003.
C. G. Bower, Transforms

Rows 1-10 are A000670, A004123, A032033, A094417, A094418, A094419, A238464, A238465, A238466, A238467.
Columns include A014105, A094421.
Main diagonal is A094420.
Antidiagonal sums are A094422.
Cf. A019538, A131689, A344499.

A094416:= func< n,k | (&+[Factorial(j)*n^j*StirlingSecond(k,j): j in [0..k]]) >;
[A094416(n-k+1,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Jan 12 2024

Bo[, 0]=1; Bo[r, n_]:= Bo[r, n]= r*Sum[Binomial[n,k] Bo[r,n-k], {k, n}];
Table[Bo[r-n+1, n], {r, 10}, {n, r}] // Flatten (* Jean-François Alcover, Nov 03 2018 *)

# The Akiyama-Tanigawa algorithm applied to the powers of r + 1
# generates the rows. Adds one row (r=0) and one column (n=0).
# Adapted from Peter Luschny on A371568.
def f(n, r): return (r + 1)**n
def ATtransform(r, len, f):
  A = [0] * len
  R = [0] * len
  for n in range(len):
      R[n] = f(n, r)
      for j in range(n, 0, -1):
          R[j - 1] = j * (R[j] - R[j - 1])
      A[n] = R[0]
  return A
for r in range(8): print([r], ATtransform(r, 8, f)) # Shel Kaphan, May 03 2024

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

E.g.f.: 1/(1 + r*(1 - exp(x))).
Bo(r, n) = Sum_{k=0..n} k!*r^k*Stirling2(n, k) = 1/(r+1) * Sum_{k>=1} k^n * (r/(r+1))^k, for r>0, n>0.
Recurrence: Bo(r, n) = r * Sum_{k=1..n} C(n, k)*Bo(r, n-k), with Bo(r, 0) = 1.
Bo(r,0) = 1, Bo(r,n) = r*Bo(r,n-1) - (r+1)*Sum_{j=1..n-1} (-1)^j * binomial(n-1,j) * Bo(r,n-j). - Seiichi Manyama, Nov 17 2023

Offset corrected by Geoffrey Critzer, Jan 01 2019

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

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

Original entry on oeis.org

1, 3, 30, 489, 11127, 325218, 11612595, 489926559, 23846152332, 1315294430043, 81078316924035, 5523729981650004, 412148874577007037, 33425421047034028743, 2927620572178735480350, 275410244285003264624949, 27695140477706524122414867

Offset: 0

Seiichi Manyama, May 21 2022

nonn

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

Cf. A305404, A354242, A354253.
Cf. A094419, A346985, A354252, A365556, A365557.
Cf. A011781.

my(N=20, x='x+O('x^N)); Vec(serlaplace(1/sqrt(7-6*exp(x))))

my(N=20, x='x+O('x^N)); Vec(serlaplace(sum(k=0, N, binomial(2*k, k)*(3*(exp(x)-1)/2)^k)))

a(n) = sum(k=0, n, (3/2)^k*(2*k)!*stirling(n, k, 2)/k!);

E.g.f.: Sum_{k>=0} binomial(2*k,k) * (3 * (exp(x) - 1)/2)^k.
a(n) = Sum_{k=0..n} (3/2)^k * (2*k)! * Stirling2(n,k)/k!.
a(n) ~ sqrt(2/7) * n^n / (exp(n) * log(7/6)^(n + 1/2)). - Vaclav Kotesovec, Jun 04 2022
a(0) = 1; a(n) = Sum_{k=1..n} (6 - 3*k/n) * binomial(n,k) * a(n-k). - Seiichi Manyama, Sep 09 2023
a(0) = 1; a(n) = 3*a(n-1) - 7*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 17 2023

A238464 Generalized ordered Bell numbers Bo(7,n).

Original entry on oeis.org

1, 7, 105, 2359, 70665, 2646007, 118893705, 6232661239, 373405001865, 25167452766967, 1884759251911305, 155262005162499319, 13952854271421949065, 1358385484966283220727, 142418920493123648992905, 15998363870912950298468599

Offset: 0

Vincenzo Librandi, Mar 17 2014

nonn

Row 7 of array A094416, which has more information.

Vincenzo Librandi, Table of n, a(n) for n = 0..250

Cf. A000670, A004123, A032033, A094417, A094418, A094419.

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

t=30; Range[0, t]! CoefficientList[Series[1/(8 - 7 Exp[x]), {x, 0, t}], x]

x='x+O('x^66); Vec(serlaplace(1/(8 - 7*exp(x)))) \\ Joerg Arndt, Mar 17 2014

E.g.f.: 1/(8 - 7*exp(x)).
a(n) ~ n! / (8*(log(8/7))^(n+1)). - Vaclav Kotesovec, Mar 20 2014
a(0) = 1; a(n) = 7*a(n-1) - 8*Sum_{k=1..n-1} (-1)^k * binomial(n-1,k) * a(n-k). - Seiichi Manyama, Nov 17 2023

A345077 a(0) = 1; a(n) = 6 * Sum_{k=1..n} binomial(n,k) * a(k-1).

Original entry on oeis.org

1, 6, 48, 414, 3876, 38946, 416808, 4722774, 56379756, 706236426, 9250945008, 126342991614, 1794459834036, 26445918969906, 403610795535288, 6367606516836774, 103683034842399996, 1739933892930544986, 30052751213767045248, 533635421576480845134, 9730601644306627161156

Offset: 0

Ilya Gutkovskiy, Jun 07 2021

nonn

Cf. A040027, A094419, A144223, A343523, A343975, A344735, A344840, A345078, A345081.

a[0] = 1; a[n_] := a[n] = 6 Sum[Binomial[n, k] a[k - 1], {k, 1, n}]; Table[a[n], {n, 0, 20}]
nmax = 20; A[] = 0; Do[A[x] = 1 + 6 x A[x/(1 - x)]/(1 - x)^2 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]

G.f. A(x) satisfies: A(x) = 1 + 6 * x * A(x/(1 - x)) / (1 - x)^2.

A384514 Expansion of e.g.f. 6/(7 - exp(6*x)).

Original entry on oeis.org

1, 1, 8, 78, 960, 14736, 272448, 5881968, 145105920, 4026744576, 124159039488, 4211132779008, 155814875873280, 6245695887446016, 269610827961212928, 12469729905669224448, 615184657168540631040, 32246522356406129197056, 1789714914567248392224768

Offset: 0

Seiichi Manyama, Jun 01 2025

nonn

Wikipedia, Polylogarithm.

Cf. A326323.

Cf. A094419, A384521, A384522, A384523, A384524.

PARI
```
a(n) = (-6)^(n+1)*polylog(-n, 7)/7;
```

a(n) = (-6)^(n+1)/7 * Li_{-n}(7), where Li_{n}(x) is the polylogarithm function.
a(n) = 6^(n+1) * Sum_{k>=0} k^n * (1/7)^(k+1).
a(n) = Sum_{k=0..n} 6^(n-k) * k! * Stirling2(n,k).
a(n) = (1/7) * Sum_{k=0..n} 7^k * (-6)^(n-k) * k! * Stirling2(n,k) for n > 0.
a(0) = 1; a(n) = Sum_{k=1..n} 6^(k-1) * binomial(n,k) * a(n-k).
a(0) = 1; a(n) = a(n-1) + 7 * Sum_{k=1..n-1} (-6)^(k-1) * binomial(n-1,k) * a(n-k).

cp's OEIS Frontend

A131689 Triangle of numbers T(n,k) = k!*Stirling2(n,k) = A000142(k)*A048993(n,k) read by rows, T(n, k) for 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Julia

Maple

Mathematica

PARI

SageMath

Formula

A032033 Stirling transform of A032031.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

A094418 Generalized ordered Bell numbers Bo(5,n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

SageMath

Formula

A094416 Array read by antidiagonals: generalized ordered Bell numbers Bo(r,n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

Python

SageMath

Formula

Extensions

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

A131689 Triangle of numbers T(n,k) = k!Stirling2(n,k) = A000142(k)A048993(n,k) read by rows, T(n, k) for 0 <= k <= n.