A032033 -id:A032033 - cp's OEIS frontend

A354291 Expansion of e.g.f. exp(f(x) - 1) where f(x) = 1/(4 - 3*exp(x)) = e.g.f. for A032033.

1, 3, 30, 435, 8211, 190056, 5196099, 163541055, 5815620696, 230350071189, 10048990989747, 478467217544322, 24678559536271581, 1370217125170670367, 81457311857722336614, 5160975525978898855143, 347090708803947931122807, 24690132231344937537382560

Offset: 0

Seiichi Manyama, May 23 2022

nonn

Cf. A075729, A354290.

Cf. A032033, A354287, A354289.

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

a_vector(n) = my(v=vector(n+1)); v[1]=1; for(i=1, n, v[i+1]=sum(j=1, i, sum(k=0, j, 3^k*k!*stirling(j, k, 2))*binomial(i-1, j-1)*v[i-j+1])); v;

a(0) = 1; a(n) = Sum_{k=1..n} A032033(k) * binomial(n-1,k-1) * a(n-k).
a(n) = Sum_{k=0..n} 3^k * A000262(k) * Stirling2(n,k).
a(n) ~ exp(-7/8 - n + 1/(8*log(4/3)) + sqrt(n/log(4/3))) * n^(n - 1/4) / (2*log(4/3)^(n + 1/4)). - Vaclav Kotesovec, May 23 2022

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.

Original entry on oeis.org

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

A004123 Number of generalized weak orders on n points.

Original entry on oeis.org

1, 2, 10, 74, 730, 9002, 133210, 2299754, 45375130, 1007179562, 24840104410, 673895590634, 19944372341530, 639455369290922, 22079273878443610, 816812844197444714, 32232133532123179930, 1351401783010933015082

Offset: 1

N. J. A. Sloane

Number of bipartitional relations on a set of cardinality n. - Ralf Stephan, Apr 27 2003
From Peter Bala, Jul 08 2022: (Start)
Conjecture: Let k be a positive integer. The sequence obtained by reducing a(n) modulo k is eventually periodic with the period dividing phi(k) = A000010(k). For example, modulo 7 we obtain the sequence [1, 2, 3, 4, 2, 0, 0, 2, 3, 4, 2, 0, 0, 2, 3, 4, 2, 0, 0, ...] with an apparent period of 6 = phi(7) starting at a(2). Cf. A000670.
More generally, we conjecture that the same property holds for integer sequences having an e.g.f. of the form G(exp(x) - 1), where G(x) is an integral power series. (End)

L Santocanale, F Wehrung, G Grätzer, F Wehrung, Generalizations of the Permutohedron, in Grätzer G., Wehrung F. (eds) Lattice Theory: Special Topics and Applications. Birkhäuser, Cham, pp. 287-397; DOI https://doi.org/10.1007/978-3-319-44236-5_8
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 1..100
Paul Barry, Eulerian-Dowling Polynomials as Moments, Using Riordan Arrays, arXiv:1702.04007 [math.CO], 2017.
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
D. Foata and C. Krattenthaler, Graphical Major Indices, II, Seminaire Lotharingien de Combinatoire, B34k, 16 pp., 1995.
D. Foata and D. Zeilberger, The Graphical Major Index, arXiv:math/9406220 [math.CO], 1994.
Jacob Sprittulla, On Colored Factorizations, arXiv:2008.09984 [math.CO], 2020.
Carl G. Wagner, Enumeration of generalized weak orders, Arch. Math. (Basel) 39 (1982), no. 2, 147-152.
C. G. Wagner, Enumeration of generalized weak orders, Preprint, 1980. [Annotated scanned copy]
C. G. Wagner and N. J. A. Sloane, Correspondence, 1980

Cf. A004121, A004122, A000165, A000670, A032033.
Second row of array A094416 (generalized ordered Bell numbers).
Equals 2 * A050351(n) for n>0.

a[n_] := (1/3)*PolyLog[-n + 1, 2/3]; a[1]=1; Table[a[n], {n, 1, 18}] (* Jean-François Alcover, Jun 11 2012 *)
CoefficientList[Series[1/(3-2*Exp[x]), {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, Aug 07 2013 *)

{a(n)=polcoeff(sum(m=0, n, 2^m*m!*x^(m+1)/prod(k=1, m, 1-k*x+x*O(x^n))), n)} /* Paul D. Hanna, Jul 20 2011 */

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

A004123 = lambda n: sum(stirling_number2(n-1,k)*(2^k)*factorial(k) for k in (0..n-1))
[A004123(n) for n in (1..18)] # Peter Luschny, Jan 18 2016

E.g.f. for sequence with offset 0: 1/(3-2*exp(x)).
a(n) = 2^n*A(n,3/2); A(n,x) the Eulerian polynomials. - Peter Luschny, Aug 03 2010
O.g.f.: Sum_{n>=0} 2^n*n!*x^(n+1)/Product_{k=0..n} (1-k*x). - Paul D. Hanna, Jul 20 2011
a(n) = Sum_{k>=0} k^n*(2/3)^k/3.
a(n) = Sum_{k=0..n} Stirling2(n, k)*(2^k)*k!.
Stirling transform of A000165. - Karol A. Penson, Jan 25 2002
"AIJ" (ordered, indistinct, labeled) transform of 2, 2, 2, 2, ...
Recurrence: a(n) = 2*Sum_{k=1..n} binomial(n, k)*a(n-k), a(0)=1. - Vladeta Jovovic, Mar 27 2003
a(n) ~ (n-1)!/(3*(log(3/2))^n). - Vaclav Kotesovec, Aug 07 2013
a(n) = log(3/2)*Integral_{x>=0} floor(x)^n * (3/2)^(-x) dx. - Peter Bala, Feb 14 2015
E.g.f.: (x - log(3 - 2*exp(x)))/3. - Ilya Gutkovskiy, May 31 2018
Conjectural o.g.f. as a continued fraction of Stieltjes type:  1/(1 - 2*x/(1 - 3*x/(1 - 4*x/(1 - 6*x/(1 - ... - 2*n*x/(1 - 3*n*x/(1 - ...))))))). - Peter Bala, Jul 08 2022

More terms from Christian G. Bower

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)

A094419 Generalized ordered Bell numbers Bo(6,n).

Original entry on oeis.org

1, 6, 78, 1518, 39390, 1277646, 49729758, 2258233998, 117196187550, 6842432930766, 443879517004638, 31674687990494478, 2465744921215207710, 207943837884583262286, 18885506918597311159518, 1837699347783655374914958, 190743171535070652261555870, 21035482423625416328497024206

Offset: 0

Ralf Stephan, May 02 2004

nonn

Sixth row of array A094416, which has more information.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Wikipedia, Polylogarithm.

Cf. A032033, A094416, A094417, A094418, A131689.
Cf. A346985, A354252, A365555, A365556, A365557.
Cf. A384514, A384521, A384522, A384523.

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

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

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

a(n) = (-1)^(n+1)*polylog(-n, 7/6)/7; \\ Seiichi Manyama, Jun 01 2025

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

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

A201354 Expansion of e.g.f. exp(x) / (4 - 3*exp(x)).

Original entry on oeis.org

1, 4, 28, 292, 4060, 70564, 1471708, 35810212, 995827420, 31153998244, 1082931514588, 41407678132132, 1727226633730780, 78051253062575524, 3798351192214837468, 198049421007186054052, 11014905131587945490140, 650903915009792820650404, 40726453234725158535472348

Offset: 0

Paul D. Hanna, Nov 30 2011

			E.g.f.: E(x) = 1 + 4*x + 28*x^2/2! + 292*x^3/3! + 4060*x^4/4! + 70564*x^5/5! + ...
O.g.f.: A(x) = 1 + 4*x + 28*x^2 + 292*x^3 + 4060*x^4 + 70564*x^5 + ...
where A(x) = 1 + 4*x/(1+x) + 2!*4^2*x^2/((1+x)*(1+2*x)) + 3!*4^3*x^3/((1+x)*(1+2*x)*(1+3*x)) + 4!*4^4*x^4/((1+x)*(1+2*x)*(1+3*x)*(1+4*x)) + ...

G. C. Greubel, Table of n, a(n) for n = 0..370

Cf. A000629, A032033, A050352, A123227, A136728.

R:=PowerSeriesRing(Rationals(), 20); Coefficients(R!(Laplace( 1/(4*Exp(-x) -3) ))); // G. C. Greubel, Jun 08 2020

seq(coeff(series(1/(4*exp(-x) -3), x, n+1)*n!, x, n), n = 0..20); # G. C. Greubel, Jun 08 2020

Table[Sum[(-1)^(n-k)*4^k*StirlingS2[n,k]*k!,{k,0,n}],{n,0,20}] (* Vaclav Kotesovec, Jun 13 2013 *)

{a(n)=n!*polcoeff(exp(x+x*O(x^n))/(4 - 3*exp(x+x*O(x^n))), n)}

{a(n)=polcoeff(sum(m=0, n, 4^m*m!*x^m/prod(k=1, m, 1+k*x+x*O(x^n))), n)}

{Stirling2(n, k)=if(k<0||k>n, 0, sum(i=0, k, (-1)^i*binomial(k, i)/k!*(k-i)^n))}
{a(n)=sum(k=0, n, (-1)^(n-k)*4^k*Stirling2(n, k)*k!)}

[sum( (-1)^(n-j)*4^j*factorial(j)*stirling_number2(n,j) for j in (0..n)) for n in (0..20)] # G. C. Greubel, Jun 08 2020

O.g.f.: A(x) = Sum_{n>=0} n! * 4^n*x^n / Product_{k=0..n} (1+k*x).
O.g.f.: A(x) = 1/(1 - 4*x/(1-3*x/(1 - 8*x/(1-6*x/(1 - 12*x/(1-9*x/(1 - 16*x/(1-12*x/(1 - 20*x/(1-15*x/(1 - ...))))))))))), a continued fraction.
a(n) = Sum_{k=0..n} (-1)^(n-k) * 4^k * Stirling2(n,k) * k!.
a(n) = 4*A050352(n) for n>0.
a(n) = Sum_{k=0..n} A123125(n,k)*4^k*3^(n-k). - Philippe Deléham, Nov 30 2011
a(n) = log(4/3) * Integral_{x = 0..oo} (ceiling(x))^n * (4/3)^(-x) dx. - Peter Bala, Feb 06 2015
G.f.: 2/G(0), where G(k) = 1 + 1/(1 - 8*x*(k+1)/(8*x*(k+1) - 1 + 6*x*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
a(n) ~ n! / (3*(log(4/3))^(n+1)). - Vaclav Kotesovec, Jun 13 2013
a(n) = 1 + 3 * Sum_{k=0..n-1} binomial(n,k) * a(k). - Ilya Gutkovskiy, Jun 08 2020
From Seiichi Manyama, Nov 15 2023: (Start)
a(0) = 1; a(n) = -4*Sum_{k=1..n} (-1)^k * binomial(n,k) * a(n-k).
a(0) = 1; a(n) = 4*a(n-1) + 3*Sum_{k=1..n-1} binomial(n-1,k) * a(n-k). (End)
a(n) = (4/3)*A032033(n) - (1/3)*0^n. - Seiichi Manyama, Dec 21 2023

A050352 Number of 4-level labeled linear rooted trees with n leaves.

Original entry on oeis.org

1, 1, 7, 73, 1015, 17641, 367927, 8952553, 248956855, 7788499561, 270732878647, 10351919533033, 431806658432695, 19512813265643881, 949587798053709367, 49512355251796513513, 2753726282896986372535, 162725978752448205162601

Offset: 0

Christian G. Bower, Oct 15 1999

nonn

G. C. Greubel, Table of n, a(n) for n = 0..370
Marian Muresan, A concrete approach to classical analysis, CMS Books in Mathematics (2009) Table 10.2.
Norihiro Nakashima, Shuhei Tsujie, Enumeration of Flats of the Extended Catalan and Shi Arrangements with Species, arXiv:1904.09748 [math.CO], 2019.
Index entries for sequences related to rooted trees

Cf. A000670, A050351 - A050359.

Equals 1/3 * A032033(n) for n>0.

[1] cat [(&+[3^(j-1)*Factorial(j)*StirlingSecond(n,j): j in [1..n]]) : n in [0..20]]; // G. C. Greubel, Jun 08 2020

seq(coeff(series( (3-2*exp(x))/(4-3*exp(x)), x, n+1)*n!, x, n), n = 0..20); # G. C. Greubel, Jun 08 2020

With[{nn=20}, CoefficientList[Series[(3-2Exp[x])/(4-3Exp[x]),{x,0,nn}], x]*Range[0,nn]!] (* Harvey P. Dale, Aug 16 2012 *)

a(n)=n!*if(n<0,0,polcoeff((3-2*exp(x))/(4-3*exp(x))+O(x^(n+1)),n))

{a(n)=if(n==0, 1, (1/12)*round(suminf(k=1, k^n * (3/4)^k *1.)))} \\ Paul D. Hanna, Nov 28 2014

[1]+[sum( 3^(j-1)*factorial(j)*stirling_number2(n,j) for j in (1..n)) for n in (1..20)] # G. C. Greubel, Jun 08 2020

E.g.f.: (3 - 2*exp(x))/(4 - 3*exp(x)).
a(n) is asymptotic to (1/12)*n!/log(4/3)^(n+1). - Benoit Cloitre, Jan 30 2003
For m-level trees (m>1), e.g.f. is (m-1-(m-2)*e^x)/(m-(m-1)*e^x) and number of trees is 1/(m*(m-1))*sum(k>=0, (1-1/m)^k*k^n). Here m=4, so a(n)=(1/12)*sum(k>=0, (3/4)^k*k^n) (for n>0). - Benoit Cloitre, Jan 30 2003
Let f(x) = (1+x)*(1+2*x). Let D be the operator g(x) -> d/dx(f(x)*g(x)). Then for n>=1, a(n) = D^(n-1)(1) evaluated at x = 1. Compare with the result A000670(n) = D^(n-1)(1) at x = 0. See also A194649. - Peter Bala, Sep 05 2011
E.g.f.: 1 + x/(G(0)-4*x) where G(k)= x + k + 1 - x*(k+1)/G(k+1); (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Jul 11 2012
a(n) = (1/12) * Sum_{k>=1} k^n * (3/4)^k for n>0. - Paul D. Hanna, Nov 28 2014
a(n) = Sum_{k=1..n} Stirling2(n, k) * k! * 3^(k-1). - Paul D. Hanna, Nov 28 2014, after Vladeta Jovovic in A050351
a(n) = 1 + 3 * Sum_{k=1..n-1} binomial(n,k) * a(k). - Ilya Gutkovskiy, Jun 08 2020

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

A255927 a(n) = (3/4) * Sum_{k>=0} (3*k)^n/4^k.

Original entry on oeis.org

1, 1, 5, 33, 285, 3081, 40005, 606033, 10491885, 204343641, 4422082005, 105265315233, 2733583519485, 76902684021801, 2329889536156005, 75629701786875633, 2618654297178083085, 96336948993312237561, 3752590641305604502005, 154294551397830418471233, 6677999524135208461382685

Offset: 0

Karol A. Penson, Sep 03 2015

nonn

			a(5) = 729*hypergeom([2,2,2,2,2],[1,1,1,1],1/4)/16 = 3081.

Alois P. Heinz, Table of n, a(n) for n = 0..400
P. Blasiak, K. A. Penson and A. I. Solomon, Dobinski-type relations and the Log-normal distribution, arXiv:quant-ph/0303030, 2003.
P. Blasiak, K. A. Penson and A. I. Solomon, Dobinski-type relations and the Log-normal distribution, J. Phys. A: Math. Gen. 36, (2003), L273.
Eric Weisstein's World of Mathematics, Lerch Transcendent

Cf. A000670, A094418, A004123 A032033, A094417, A122704, A326323.

S:= series(3/(4-exp(3*x)), x, 51):
seq(coeff(S,x,n)*n!, n=0..50); # Robert Israel, Sep 03 2015
seq(add(combinat:-eulerian1(n,k)*4^k, k=0..n), n=0..20); # Peter Luschny, Jun 27 2019

a[n_] := 3^(n+1)/4 HurwitzLerchPhi[1/4, -n, 0];
Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Sep 18 2018 *)
Eulerian1[0, 0] = 1; Eulerian1[n_, k_] := Sum[(-1)^j (k-j+1)^n Binomial[n+1, j], {j, 0, k+1}]; Table[Sum[Eulerian1[n, k] 4^k, {k, 0, n}], {n, 0, 20}] (* Jean-François Alcover, Jul 13 2019, after Peter Luschny *)

a(n) = sum(k=0, n, stirling(n,k,2)*k!*3^(n-k)); \\ Michel Marcus, Sep 03 2015

a(n) = Sum_{k>=0} Stirling2(n,k)*k!*3^(n-k).
E.g.f.: 3/(4-exp(3*x)).
Special values of the generalized hypergeometric function n_F_(n-1):
a(n) = (3^(n+1)/16) * hypergeom([2,2,..2],[1,1,..1],1/4), where the sequence in the first square bracket ("upper" parameters) has n elements all equal to 2 whereas the sequence in the second square bracket ("lower" parameters) has n-1 elements all equal to 1.
Example: a(5) = 729 * hypergeom([2,2,2,2,2],[1,1,1,1],1/4)/16 = 3081.
a(n) is the n-th moment of the discrete weight function W(x) = (3/4)*sum(k>=0, Dirac(x-3*k)/4^k), n>=1.
a(n) ~ n! * 3^(n+1) / ((log(2))^(n+1) * 2^(n+3)). - Vaclav Kotesovec, Jul 09 2018
G.f.: Sum_{j>=0} j!*x^j / Product_{k=1..j} (1 - 3*k*x). - Ilya Gutkovskiy, Apr 04 2019
a(n) = A_{4}(n) where A_{n}(x) are the Eulerian polynomials as defined in A326323. - Peter Luschny, Jun 27 2019

a(0)=1 prepended by Alois P. Heinz, Sep 18 2018

cp's OEIS Frontend

A354291 Expansion of e.g.f. exp(f(x) - 1) where f(x) = 1/(4 - 3*exp(x)) = e.g.f. for A032033.

Views

Author

Keywords

Crossrefs

Programs

PARI

PARI

Formula

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

A004123 Number of generalized weak orders on n points.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Sage

Formula

Extensions

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

A094419 Generalized ordered Bell numbers Bo(6,n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

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

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.