A074059 -id:A074059 - cp's OEIS frontend

A134685 Irregular triangle read by rows: coefficients C(j,k) of a partition transform for direct Lagrange inversion.

1, -1, 3, -1, -15, 10, -1, 105, -105, 10, 15, -1, -945, 1260, -280, -210, 35, 21, -1, 10395, -17325, 6300, 3150, -280, -1260, -378, 35, 56, 28, -1, -135135, 270270, -138600, -51975, 15400, 34650, 6930, -2100, -1575, -2520, -630, 126, 84, 36, -1

Offset: 1

Tom Copeland, Jan 26 2008, Sep 13 2008

Let f(t) = u(t) - u(0) = Ev[exp(u.* t) - u(0)] = log{Ev[(exp(z.* t))/z_0]} = Ev[-log(1- a.* t)], where the operator Ev denotes umbral evaluation of the umbral variables u., z. or a., e.g., Ev[a.^n + a.^m] = a_n + a_m . The relation between z_n and u_n is given in reference in A127671 and u_n = (n-1)! * a_n .
If u_1 is not equal to 0, then the compositional inverse for these expressions is given by g(t) = Sum_{j>=1} P(j,t) where, with u_n denoted by (n') for brevity,
P(1,t) = (1')^(-1) * [ 1 ] * t
P(2,t) = (1')^(-3) * [ -(2') ] * t^2 / 2!
P(3,t) = (1')^(-5) * [ 3 (2')^2 - (1')(3') ] * t^3 / 3!
P(4,t) = (1')^(-7) * [ -15 (2')^3 + 10 (1')(2')(3') - (1')^2 (4') ] * t^4 / 4!
P(5,t) = (1')^(-9) * [ 105 (2')^4 - 105 (1') (2')^2 (3') + 10 (1')^2 (3')^2 + 15 (1')^2 (2') (4') - (1')^3 (5') ] * t^5 / 5!
P(6,t) = (1')^(-11) * [ -945 (2')^5 + 1260 (1') (2')^3 (3') - 280 (1')^2 (2') (3')^2 - 210 (1')^2 (2')^2 (4') + 35 (1')^3 (3')(4') + 21 (1')^3 (2')(5') - (1')^4 (6') ] * t^6 / 6!
P(7,t) = (1')^(-13) * [ 10395 (2')^6 - 17325 (1') (2')^4 (3') +  (1')^2 [ 6300 (2')^2 (3')^2 + 3150  (2')^3 (4')] - (1')^3 [280 (3')^3 + 1260  (2')(3')(4') + 378  (2')^2(5')] + (1')^4 [35 (4')^2 + 56 (3')(5') + 28 (2')(6')]  - (1')^5 (7') ] * t^7 / 7!
P(8,t) = (1')^(-15) * [ -135135 (2')^7 + 270270 (1') (2')^5 (3') -  (1')^2 [ 138600 (2')^3 (3')^2 +  51975 (2')^4 (4')] + (1')^3 [15400 (2')(3')^3 + 34650  (2')^2(3')(4') + 6930  (2')^3(5')] - (1')^4 [2100 (3')^2(4') + 1575 (2')(4')^2 + 2520 (2')(3')(5') + 630 (2')^2(6') ]  + (1')^5 [126 (4')(5') + 84 (3')(6') + 36 (2')(7')] - (1')^6 (8') ] * t^8 / 8!
...
Substituting ((m-1)') for (m') in each partition and ignoring the (0') factors, the partitions in the brackets of P(n,t) become those of n-1 listed in Abramowitz and Stegun on page 831 (in the reversed order) and the number of partitions in P(n,t) is given by A000041(n-1).
Combinatorial interpretations are given in the link.
From Tom Copeland, Jul 10 2018: (Start)
Coefficients occurring in prolongation for the special Euclidean group SE(2) and special affine group SA(2) in the Olver presentation on moving frames (MFP) in slides 33 and 42. These are a result of applying an iterated derivative of the form h(x)d/dx = d/dy as in this entry (more generally as g(x) d/dx as discussed in A145271). See also p. 6 of Olver's paper on contact forms, but note that the 12 should be a 15 in the formula for the compositional inverse of S(t).
Change variables in the MFP to obtain connections to the partition polynomials Prt_n = n! * P(n,1) above. Let delta and beta in the formulas for the equi-affine curves in MFP be L and B, respectively, and D_y = (1/(L-B*u_x)) d/dx = (1/w_x) d/dx. Then v_(yy) = (1/B) [-w_(xx)/(w_x)^3] in MFP (there is an overall sign error in MFP for v_(yy) and higher derivatives w.r.t. y), and (d/dy)^n v = v_n = (1/B)* [(1/w_1)*(d/dx)]^(n-2) [-w_2/(w_1)^3] for n > 1, with w_n = (d/dx)^n w. Consequently, in the partition polynomials Prt_n for n > 1 here substitute (n') = -B*u_n = w_n for n > 1 and (1')  = L-B*u_1 = w_1, where u_n = (d/dx)^n u, and then divide by B. For example, v_4 = (1/B)*Prt_4 = (1/B)*4!*P(4,1) = (1/B) (L-B*u_n)^(-7) [-15*(-B*u_2)^3 + 10 (L-B*u_1)(-B*u_2)(-B*u_3) - (L-B*u_1)^2 (-B*u_4)], agreeing with v_4 in MFP except for the overall sign.
For the SE(2) transformation formulas in MFP, let w_x = cos(phi) + sin(phi)*u_x, and then the same transformations apply as above with cos(phi) and sin(phi) substituted for L and -B, respectively. (End)

			Examples and checks:
1) Let u_1 = -1 and u_n = 1 for n>1,
then f(t) = exp(u.*t) - u(0) = exp(t)-2t-1
and g(t) = [e.g.f. of signed A000311];
therefore, the row sums of unsigned [C(j,k)] are A000311 =
(0,1,1,4,26,236,2752,...) = (0,-P(1,1),2!*P(2,1),-3!*P(3,1),4!*P(4,1),...).
2) Let u_1 = -1 and u_n = (n-1)! for n>1,
then f(t) = -log(1-t)-2t
and g(t) = [e.g.f. of signed (0,A032188)]
with (0,A032188) = (0,1,1,5,41,469,6889,...) = (0,-P(1,1),2!*P(2,1),-3!P(3,1),...).
3) Let u_1 = -1 and u_n = (-1)^n (n-2)! for n>1, then
f(t) = (1+t) log(1+t) - 2t
and g(t) = [e.g.f. of signed (0,A074059)]
with (0,A074059) = (0,1,1,2,7,34,213,...) = (0,-P(1,1),2!*P(2,1),-3!*P(3,1),...).
4) Let u_1 = 1, u_2 = -1 and u_n = 0 for n>2,
then f(t) = t(1-t/2)
and g(t) = [e.g.f. of (0,A001147)] = 1 - (1-2t)^(1/2)
with (0,A001147) = (0,1,1,3,15,105,945...) =(0,P(1,1),2!*P(2,1),3!*P(3,1),...).
5) Let u_1 = 1, u_2 = -2 and u_n = 0 for n>2,
then f(t)= t(1-t)
and g(t) = t * [o.g.f. of A000108] = [1 - (1-4t)^(1/2)] / 2
with (0,A000108) = (0,1,1,2,5,14,42,...) = (0,P(1,1),P(2,1),P(3,1),...).
.
From _Peter Luschny_, Feb 19 2021: (Start)
Triangle starts:
 [1]  1;
 [2] -1;
 [3]  3,     -1;
 [4] -15,     10,    -1;
 [5]  105,   -105,   [10, 15],  -1;
 [6] -945,    1260,  [-280, -210], [35, 21],  -1;
 [7]  10395, -17325, [6300, 3150], [-280, -1260, -378], [35, 56, 28], -1;
 [8] -135135, 270270, [-138600, -51975], [15400, 34650, 6930], [-2100, -1575, -2520, -630], [126, 84, 36], -1
The coefficients can be seen as a refinement of the Ward numbers: Let R(n, k) = Sum T(n, k), where the sum collects adjacent terms with equal sign, as indicated by the square brackets in the table, then R(n+1, k+1) = (-1)^(n-k)*W(n, k), where W(n, k) are the Ward numbers A181996, for n >= 0 and 0 <= k <= n-1.  (End)

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, p. 831.
D. S. Alexander, A History of Complex Dynamics: From Schröder to Fatou to Julia, Friedrich Vieweg & Sohn, 1994, p. 10.
J. Riordan, Combinatorial Identities, Robert E. Krieger Pub. Co., 1979, (unsigned partition polynomials in Table 5.2 on p. 181, but may have errors).

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
O. Arnaldsson, Élie Cartan’s Theory of Moving Frames, slide presentation, 2014.
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms, 2015.
Tom Copeland, Important formulas in combinatorics, MathOverflow answer, 2015.
Tom Copeland, Compositional inversion and generating functions in algebraic geometry, MathOverflow question, 2014.
Tom Copeland, Compositional inverse pairs, the inviscid Burgers-Hopf equation, and the Stasheff associahedra, 2014.
Tom Copeland, Lagrange a la Lah, 2011.
Tom Copeland, Short Note on Lagrange Inversion, 2008.
Tom Copeland, Formal group laws and binomial Sheffer sequences, 2018.
Bishal Deb and Alan D. Sokal, Higher-order Stirling cycle and subset triangles: Total positivity, continued fractions and real-rootedness, arXiv:2507.18959 [math.CO], 2025. See p. 32.
Hector Figueroa and Jose M. Gracia-Bondia, Combinatorial Hopf algebras in quantum field theory I, arXiv:0408145 [hep-th], 2005, (p. 38).
Ezra Getzler, The semi-classical approximation for modular operads, arXiv:alg-geom/9612005, 1996 (see p. 2).
Peter J. Olver, The canonical contact form.
Peter J. Olver, Moving Frames, slide presentation, 2009.
Andrew Elvey Price and Alan D. Sokal, Phylogenetic trees, augmented perfect matchings, and a Thron-type continued fraction (T-fraction) for the Ward polynomials, arXiv:2001.01468 [math.CO], 2020.
Jair Patrick Taylor, Formal group laws and hypergraph colorings, doctoral thesis, Univ. of Wash., 2016, p. 95.

Cf. A145271, (A134991, A019538) = (reduced array, associated g(x)).
Cf. A000108, A000311, A001147, A032188, A049019, A074059, A129185, A133314, A139605, A036040.
Cf. A181996 (Ward numbers).

rows[n_] := {{1}}~Join~Module[{h = 1/(1 + Sum[u[k] y^k/k!, {k, n-1}] + O[y]^n), g = y, r}, r = Reap[Do[g = h D[g, y]; Sow[Expand[Normal@g /. {y -> 0}]], {k, n}]][[2, 1, ;;]]; Table[Coefficient[r[[k]], Product[u[t], {t, p}]], {k, 2, n}, {p, Reverse@Sort[Sort /@ IntegerPartitions[k-1]]}]];
rows[8] // Flatten (* Andrei Zabolotskii, Feb 19 2024 *)

The bracketed partitions of P(n,t) are of the form (u_1)^e(1) (u_2)^e(2) ... (u_n)^e(n) with coefficients given by (-1)^(n-1+e(1)) * [2*(n-1)-e(1)]! / [2!^e(2)*e(2)!*3!^e(3)*e(3)! ... n!^e(n)*e(n)! ].
From Tom Copeland, Sep 05 2011: (Start)
Let h(t) = 1/(df(t)/dt)
  = 1/Ev[u.*exp(u.*t)]
  = 1/(u_1+(u_2)*t+(u_3)*t^2/2!+(u_4)*t^3/3!+...),
  an e.g.f. for the partition polynomials of A133314
  (signed A049019) with an index shift.
  Then for the partition polynomials of A134685,
  n!*P(n,t) = ((t*h(y)*d/dy)^n) y evaluated at y=0,
  and the compositional inverse of f(t) is
  g(t) = exp(t*h(y)*d/dy) y evaluated at y=0.
  Also, dg(t)/dt = h(g(t)). (Cf. A000311 and A134991)(End)
From Tom Copeland, Oct 30 2011: (Start)
With exp[x* PS(.,t)] = exp[t*g(x)]=exp[x*h(y)d/dy] exp(t*y) eval. at y=0, the raising/creation and lowering/annihilation operators
defined by R PS(n,t)=PS(n+1,t) and L PS(n,t)= n*PS(n-1,t) are
R = t*h(d/dt) = t * 1/[u_1+(u_2)*d/dt+(u_3)*(d/dt)^2/2!+...],  and
L = f(d/dt)=(u_1)*d/dt+(u_2)*(d/dt)^2/2!+(u_3)*(d/dt)^3/3!+....
Then  P(n,t) = (t^n/n!) dPS(n,z)/dz  eval. at z=0. (Cf. A139605, A145271, and link therein to Mathemagical Forests for relation to planted trees on p. 13.) (End)
The bracketed partition polynomials of P(n,t) are also given by (d/dx)^(n-1) 1/[u_1 + u_2 * x/2! + u_3 * x^2/3! + ... + u_n * x^(n-1)/n!]^n evaluated at x=0. - Tom Copeland, Jul 07 2015
Equivalent matrix computation: Multiply the m-th diagonal (with m=1 the index of the main diagonal) of the lower triangular Pascal matrix by u_m = (d/dx)^m f(x) evaluated at x=0 to obtain the matrix UP with UP(n,k) = binomial(n,k) u_{n+1-k}. Then P(n,t) = (1, 0, 0, 0, ...) [UP^(-1) * S]^(n-1) FC * t^n/n!, where S is the shift matrix A129185, representing differentiation in the basis x^n//n!, and FC is the first column of UP^(-1), the inverse matrix of UP. These results follow from A145271 and A133314. - Tom Copeland, Jul 15 2016
Also, P(n,t) = (1, 0, 0, 0, ...) [UP^(-1) * S]^n (0, 1, 0, ..)^T * t^n/n! in agreement with A139605. - Tom Copeland, Aug 27 2016
From Tom Copeland, Sep 20 2016: (Start)
Let PS(n,u1,u2,...,un) = P(n,t) / (t^n/n!), i.e., the square-bracketed part of the partition polynomials in the expansion for the inverse in the comment section, with u_k = uk.
Also let PS(n,u1=1,u2,...,un) = PB(n,b1,b2,...,bK,...) where each bK represents the partitions of PS, with u1 = 1, that have K components or blocks, e.g., PS(5,1,u2,...,u5) = PB(5,b1,b2,b3,b4) = b1 + b2 + b3 + b4 with b1 = -u5, b2 = 15 u2 u4 + 10 u3^2, b3 = -105 u2^2 u3,  and b4 = 105 u2^4.
The relation between solutions of the inviscid Burgers' equation and compositional inverse pairs (cf. link and A086810) implies, for n > 2,  PB(n, 0 * b1, 1 * b2,..., (K-1) * bK, ...) = (1/2) * Sum_{k = 2..n-1} binomial(n+1,k) * PS(n-k+1,u_1=1,u_2,...,u_(n-k+1)) * PS(k,u_1=1,u_2,...,u_k).
For example, PB(5,0 * b1, 1 * b2, 2 * b3, 3 * b4) = 3 * 105 u2^4 - 2 * 105 u2^2 u3 + 1 * 15 u2 u4 + 1 * 10 u3^2 - 0 * u5 = 315 u2^4 - 210 u2^2 u3 + 15 u2 u4 + 10 u3^2 = (1/2) [2 * 6!/(4!*2!) * PS(2,1,u2) * PS(4,1,u2,...,u4) + 6!/(3!*3!) * PS(3,1,u2,u3)^2] = (1/2) * [ 2 * 6!/(4!*2!)  * (-u2) (-15 u2^3 + 10 u2 u3 - u4) + 6!/(3!*3!) * (3 u2^2 - u3)^2].
Also, PB(n,0*b1,1*b2,...,(K-1)*bK,...) =  d/dt t^(n-2)*PS(n,u1=1/t,u2,...,un)|{t=1} = d/dt (1/t)*PS(n,u1=1,t*u2,...,t*un)|{t=1}.
(End)
A recursion relation for computing each partition polynomial of this entry from the lower order polynomials and the coefficients of the Bell polynomials of A036040 is presented in the blog entry "Formal group laws and binomial Sheffer sequences." - Tom Copeland, Feb 06 2018

P(7,t) and P(8,t) data added by Tom Copeland, Jan 14 2016

Terms in rows 5-8 reordered by Andrei Zabolotskii, Feb 19 2024

A074060 Graded dimension of the cohomology ring of the moduli space of n-pointed stable curves of genus 0 satisfying the associativity equations of physics (also known as the WDVV equations).

Original entry on oeis.org

1, 1, 1, 1, 5, 1, 1, 16, 16, 1, 1, 42, 127, 42, 1, 1, 99, 715, 715, 99, 1, 1, 219, 3292, 7723, 3292, 219, 1, 1, 466, 13333, 63173, 63173, 13333, 466, 1, 1, 968, 49556, 429594, 861235, 429594, 49556, 968, 1, 1, 1981, 173570, 2567940, 9300303, 9300303, 2567940, 173570, 1981, 1

Offset: 3

Margaret A. Readdy, Aug 16 2002

Combinatorial interpretations of Lagrange inversion (A134685) and the 2-Stirling numbers of the first kind (A049444 and A143491) provide a combinatorial construction for A074060 (see first Copeland link). For relations of A074060 to other arrays see second Copeland link page 19. - Tom Copeland, Sep 28 2008

These Poincare polynomials for the compactified moduli space of rational curves are presented on p. 5 of Lando and Zvonkin as well as those for the non-compactified Poincare polynomials of A049444 in factorial form. - Tom Copeland, Jun 13 2021

			Viewed as a triangular array, the values are
  1;
  1,   1;
  1,   5,   1;
  1,  16,  16,   1;
  1,  42, 127,  42,   1; ...

Tom Copeland, Combinatorics of OEIS-A074060, Posted Sept. 2008.
Tom Copeland, Mathemagical Forests v2, Posted June 2008.
S. Keel, Intersection theory of moduli space of stable n-pointed curves of genus zero, Trans. Amer. Math. Soc. 330 (1992), 545-574.
M. Kontsevich and Y. Manin, Quantum cohomology of a product, (with Appendix by R. Kaufmann), Inv. Math. 124, f. 1-3 (1996) 313-339.
M. Kontsevich and Y. Manin, Quantum cohomology of a product, arXiv:q-alg/9502009, 1995.
S. Lando and A. Zvonkin, Graphs on surfaces and their applications, Encyclopaedia of Mathematical Sciences, 141, Springer, 2004.
Y. Manin, Generating functions in algebraic geometry and sums over trees, arXiv:alg-geom/9407005, 1994. - _Tom Copeland_, Dec 10 2011
M. A. Readdy, The pre-WDVV ring of physics and its topology, preprint, 2002.

Cf. A074059. 2nd diagonal is A002662.

DA:=((1+t)*A(u,t)+u)/(1-t*A(u,t)): F:=0: for k from 1 to 10 do F:=map(simplify,int(series(subs(A(u,t)=F,DA),u,k),u)); od: # Eric Rains, Apr 02 2005

DA = ((1+t) A[u, t] + u)/(1 - t A[u, t]); F = 0;
Do[F = Integrate[Series[DA /. A[u, t] -> F, {u, 0, k}], u], {k, 1, 10}];
(cc = CoefficientList[#, t]; cc Denominator[cc[[1]]])& /@ Drop[ CoefficientList[F, u], 2] // Flatten (* Jean-François Alcover, Oct 15 2019, after Eric Rains *)

Define offset to be 0 and P(n,t) = (-1)^n Sum_{j=0..n-2} a(n-2,j)*t^j with P(1,t) = -1 and P(0,t) = 1, then H(x,t) = -1 + exp(P(.,t)*x) is the compositional inverse in x about 0 of G(x,t) in A049444. H(x,0) = exp(-x) - 1, H(x,1) = -1 + exp( 2 + W( -exp(-2) * (2-x) ) ) and H(x,2) = 1 - (1+2*x)^(1/2), where W is a branch of the Lambert function such that W(-2*exp(-2)) = -2. - Tom Copeland, Feb 17 2008
Let offset=0 and g(x,t) = (1-t)/((1+x)^(t-1)-t), then the n-th row polynomial of the table is given by [(g(x,t)*D_x)^(n+1)]x with the derivative evaluated at x=0. - Tom Copeland, Jun 01 2008
With the notation in Copeland's comments, dH(x,t)/dx = -g(H(x,t),t). - Tom Copeland, Sep 01 2011
The term linear in x of [x*g(d/dx,t)]^n 1 gives the n-th row polynomial with offset 1. (See A134685.) - Tom Copeland, Oct 21 2011

More terms from Eric Rains, Apr 02 2005

A135338 Triangle read by rows: row n gives coefficients C(n,j) for a Sheffer sequence (binomial-type) with raising operator -x { 1 + W[ -exp(-2) * (2+D) ] } where W is the Lambert W multi-valued function.

Original entry on oeis.org

1, -1, 1, 1, -3, 1, -2, 7, -6, 1, 6, -20, 25, -10, 1, -24, 76, -105, 65, -15, 1, 120, -364, 511, -385, 140, -21, 1, -720, 2108, -2940, 2401, -1120, 266, -28, 1, 5040, -14328, 19720, -16632, 8841, -2772, 462, -36, 1, -40320, 111816, -151620, 129340, -73605, 27237, -6090, 750, -45, 1

Offset: 1

Tom Copeland, Feb 15 2008

The lowering (or delta) operator for these polynomials is L = -1 + exp{ 2 + W[ -exp(-2) * (2+D) ] } = Sum_{j >= 1} A074059(j) * D^j / j!.
The raising operator is R = -x { 1 + W[ -exp(-2) * (2+D) ] } = x { 1 + Sum_{j >= 1} (-1)^j * PW(j-1,-2) * D^j / j! }, where PW(j-1,x) are the polynomials of A042977.
W(x) here is W_-1 in the Monir reference and, about x = 0, W[ -exp(-2) * (2+x) ] = -[ 2 + Sum_{j >= 1} (-1)^j * PW(j-1,-2) * x^j / j! ].
From the relation between the delta and raising operators for associated binomial-type polynomials, A074059 = (1,1,2,7,34,...) and S = (1,-PW(0,-2),PW(1,-2),-PW(2,-2),...) = (1, -1, 0, -1, -2, -13, -74, -593, -5298, ...) form a list partition transform pair (see A133314); i.e., S and A074059 have reciprocal e.g.f.s and satisfy mutual recursion relations. Applying Faa di Bruno's formula to L gives other interesting integer relations between S and A074059.
The Bell transform of (-1)^n*factorial(n-1) if n>0, else 1. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 18 2016

			Triangle read by rows:
     1;
    -1,    1;
     1,   -3,     1;
    -2,    7,    -6,    1;
     6,  -20,    25,  -10,     1;
   -24,   76,  -105,   65,   -15,   1;
   120, -364,   511, -385,   140, -21,   1;
  -720, 2108, -2940, 2401, -1120, 266, -28, 1;
...
From _R. J. Mathar_, Mar 22 2013: (Start)
The matrix inverse starts:
     1;
     1,    1;
     2,    3,    1;
     7,   11,    6,   1;
    34,   55,   35,  10,   1;
   213,  349,  240,  85,  15,  1;
  1630, 2695, 1939, 770, 175, 21, 1;
  ... (End)

Peter Bala, Diagonals of triangles with generating function exp(t*F(x)).
F. Chapeau-Blondeau and A. Monir, Numerical Evaluation of the Lambert W Function and Application to Generation of Generalized Gaussian Noise With Exponent 1/2, IEEE Trans. on Signal Processing, Vol. 50, No. 9, Sept. 2002, p. 2160-2164.

Cf. A134685, A135494.

# The function BellMatrix is defined in A264428.
# Adds (1,0,0,0, ..) as column 0.
BellMatrix(n -> `if`(n=0,1,(-1)^n*(n-1)!), 9); # Peter Luschny, Jan 27 2016

max = 10; s = Series[Exp[t*(2*x-(1+x)*Log[1+x])], {x, 0, max}, {t, 0, max}] // Normal; c[n_, j_] := SeriesCoefficient[s, {x, 0, n}, {t, 0, j}]*n!; Table[c[n, j], {n, 1, max}, {j, 1, n}] // Flatten (* Jean-François Alcover, Apr 23 2014, after Peter Bala, duplicate of Copeland's e.g.f. *)
BellMatrix[f_Function, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len - 1}, {k, 0, len - 1}]];
rows = 12;
M = BellMatrix[Function[n, If[n == 0, 1, (-1)^n (n-1)!]], rows];
Table[M[[n, k]], {n, 2, rows}, {k, 2, n}] // Flatten (* Jean-François Alcover, Jun 26 2018, after Peter Luschny *)

# uses[bell_matrix from A264428]
# Adds a column 1,0,0,0, ... at the left side of the triangle.
bell_matrix(lambda n: (-1)^n*factorial(n-1) if n>0 else 1, 10) # Peter Luschny, Jan 18 2016

The row polynomials P(n,t) = Sum_{j=1..n} C(n,j) * t^j satisfy exp[P(.,t) * x] = exp{ -t * [(1+x) * log(1+x) - 2*x] }, with P(0,t) = 1 and [ P(.,x) + P(.,y) ]^n = P(n,x+y). Here, as in the e.g.f., the umbral maneuver P(.,t)^n = P(n,t) is assumed. See Mathworld and Wikipedia on Sheffer sequences and umbral calculus for other general formulas, including expansion theorems.
From Peter Bala, Dec 09 2011: (Start)
E.g.f.: exp(t*(2*x-(1+x)*log(1+x))) = 1 + t*x + (t^2-t)*x^2/2! + (t^3-3*t^2+t)*x^3/3! + ... (Restatement of Copeland's e.g.f. above in umbral notation with P(.,t)^n = P(n,t).).
If a triangular array has an e.g.f. of the form exp(t*F(x)) with F(0) = 0, then the o.g.f.'s for the diagonals of the triangle are rational functions in t (see the Bala link). The rational functions are the coefficients in the compositional inverse (with respect to x) (x-t*F(x))^(-1). In this case (x-t*(2*x-(1+x)*log(1+x)))^(-1) = x/(1-t) - t/(1-t)^3*x^2/2! + (t+2*t^2)/(1-t)^5*x^3/3! - (2*t+6*t^2+7*t^3)/(1-t)^7*x^4/4! + ... . So, for example, the (unsigned) third subdiagonal has o.g.f. (2*t+6*t^2+7*t^3)/(1-t)^7 = 2*t + 20*t^2 + 105*t^3 + 385*t^4 + ... .
(End)

More terms from Jean-François Alcover, Apr 23 2014

cp's OEIS Frontend

A134685 Irregular triangle read by rows: coefficients C(j,k) of a partition transform for direct Lagrange inversion.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A074060 Graded dimension of the cohomology ring of the moduli space of n-pointed stable curves of genus 0 satisfying the associativity equations of physics (also known as the WDVV equations).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A135338 Triangle read by rows: row n gives coefficients C(n,j) for a Sheffer sequence (binomial-type) with raising operator -x { 1 + W[ -exp(-2) * (2+D) ] } where W is the Lambert W multi-valued function.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Sage

Formula

Extensions

A180715 E.g.f.: A(x) = Series_Reversion[ x - Sum_{n>=2} (-x)^n/(n(n-1)/2) ].

Views

Author

Keywords

Examples

Crossrefs

Programs

Maple

PARI

PARI

PARI

Formula

A383459 Minimum number of cycles in any permutation in S_n of the highest order (A000793(n)).

Views

Author

Keywords

Comments

Examples

Crossrefs

A232693 E.g.f.: A(x) = Series_Reversion( x - 3*Sum_{n>=2} (-x)^n/(n*(n-1)) ).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

PARI

PARI

PARI

Formula

A383458 Maximum number of cycles in any permutation in S_n of the highest order (A000793(n)).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

A232693 E.g.f.: A(x) = Series_Reversion( x - 3Sum_{n>=2} (-x)^n/(n(n-1)) ).