cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-17 of 17 results.

A130561 Numbers associated to partitions, used for combinatoric interpretation of Lah triangle numbers A105278; elementary Schur polynomials / functions.

Original entry on oeis.org

1, 2, 1, 6, 6, 1, 24, 24, 12, 12, 1, 120, 120, 120, 60, 60, 20, 1, 720, 720, 720, 360, 360, 720, 120, 120, 180, 30, 1, 5040, 5040, 5040, 5040, 2520, 5040, 2520, 2520, 840, 2520, 840, 210, 420, 42, 1, 40320, 40320, 40320, 40320, 20160, 20160, 40320, 40320, 20160
Offset: 1

Views

Author

Wolfdieter Lang, Jul 13 2007

Keywords

Comments

The order of this array is according to the Abramowitz-Stegun (A-St) ordering of partitions (see A036036).
The row lengths sequence is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42, ...].
These numbers are similar to M_0, M_1, M_2, M_3, M_4 given in A111786, A036038, A036039, A036040, A117506, respectively.
Combinatorial interpretation: a(n,k) counts the sets of lists (ordered subsets) obtained from partitioning the set {1..n}, with the lengths of the lists given by the k-th partition of n in A-St order. E.g., a(5,5) is computed from the number of sets of lists of lengths [1^1,2^2] (5th partition of 5 in A-St order). Hence a(5,5) = binomial(5,2)*binomial(3,2) = 5!/(1!*2!) = 60 from partitioning the numbers 1,2,...,5 into sets of lists of the type {[.],[..],[..]}.
This array, called M_3(2), is the k=2 member of a family of partition arrays generalizing A036040 which appears as M_3 = M_3(k=1). S2(2) = A105278 (unsigned Lah number triangle) is related to M_3(2) in the same way as S2(1), the Stirling2 number triangle, is related to M_3(1). - Wolfdieter Lang, Oct 19 2007
Another combinatorial interpretation: a(n,k) enumerates unordered forests of increasing binary trees which are described by the k-th partition of n in the Abramowitz-Stegun order. - Wolfdieter Lang, Oct 19 2007
A relation between partition polynomials formed from these "refined Lah numbers" and Lagrange inversion for an o.g.f. is presented in the link "Lagrange a la Lah" along with an e.g.f. and an umbral binary operator tree representation. - Tom Copeland, Apr 12 2011
With the indeterminates (x_1,x_2,x_3,...) = (t,-c_2*t,-c_3*t,...) with c_n >0, umbrally P(n,a.) = P(n,t)|{t^n = a_n} = 0 and P(j,a.)P(k,a.) = P(j,t)P(k,t)|{t^n =a_n} = d_{j,k} >= 0 is the coefficient of x^j/j!*y^k/k! in the Taylor series expansion of the formal group law FGL(x,y) = f[f^{-1}(x)+f^{-1}(y)], where a_n are the inversion partition polynomials for calculating f(x) from the coefficients of the series expansion of f^{-1}(x) given in A133437. - Tom Copeland, Feb 09 2018
Divided by n!, the row partition polynomials are the elementary homogeneous Schur polynomials presented on p. 44 of the Bracci et al. paper. - Tom Copeland, Jun 04 2018
Also presented (renormalized) as the Schur polynomials on p. 19 of the Konopelchenko and Schief paper with associations to differential operators related to the KP hierarchy. - Tom Copeland, Nov 19 2018
Through equation 4.8 on p. 26 of the Arbarello reference, these polynomials appear in the Hirota bilinear equations 4.7 related to tau-function solutions of the KP hierarchy. - Tom Copeland, Jan 21 2019
These partition polynomials appear as Feynman amplitudes in their Bell polynomial guise (put x_n = n!c_n in A036040 for the indeterminates of the Bell polynomials) in Kreimer and Yeats and Balduf (e.g., p. 27). - Tom Copeland, Dec 17 2019
From Tom Copeland, Oct 15 2020: (Start)
With a_n = n! * b_n = (n-1)! * c_n for n > 0, represent a function with f(0) = a_0 = b_0 = 1 as an
A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!
B) ordinary generating function (o.g.f.), or formal power series: f(x) = 1/(1-b.x) = 1 + Sum_{n > 0} b_n * x^n
C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0} c_n * x^n /n.
Expansions of log(f(x)) are given in
I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] = e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials
II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] = log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.
Expansions of exp(f(x)-1) are given in
III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind
IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials
V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] = e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.
Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)
These partition polynomials are referred to as Schur functions by Segal and Wilson, who present associations with Plucker coordinates, Grassmannians, and the tau functions of the KdV hierarchy. See pages 51 and 61. - Tom Copeland, Jan 08 2022

Examples

			Triangle starts:
  [  1];
  [  2,   1];
  [  6,   6,   1];
  [ 24,  24,  12, 12,  1];
  [120, 120, 120, 60, 60, 20, 1];
  ...
a(5,6) = 20 = 5!/(3!*1!) because the 6th partition of 5 in A-St order is [1^3,2^1].
a(5,5) = 60 enumerates the unordered [1^1,2^2]-forest with 5 vertices (including the three roots) composed of three such increasing binary trees: 5*((binomial(4,2)*2)*(1*2))/2! = 5*12 = 60.
		

References

  • E. Arbarello, "Sketches of KdV", Contemp. Math. 312 (2002), p. 9-69.

Crossrefs

Cf. A105278 (unsigned Lah triangle |L(n, m)|) obtained by summing the numbers for given part number m.
Cf. A000262 (row sums), identical with row sums of unsigned Lah triangle A105278.
A134133(n, k) = A130561(n, k)/A036040(n, k) (division by the M_3 numbers). - Wolfdieter Lang, Oct 12 2007
Cf. A096162.
Cf. A133437.
Cf. A127671.

Formula

a(n,k) = n!/(Product_{j=1..n} e(n,k,j)!) with the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n. Exponents 0 can be omitted due to 0!=1.
From Tom Copeland, Sep 18 2011: (Start)
Raising and lowering operators are given for the partition polynomials formed from A130561 in the Copeland link in "Lagrange a la Lah Part I" on pp. 22-23.
An e.g.f. for the partition polynomials is on page 3:
exp[t*:c.*x/(1-c.*x):] = exp[t*(c_1*x + c_2*x^2 + c_3*x^3 + ...)] where :(...): denotes umbral evaluation of the enclosed expression and c. is an umbral coefficient. (End)
From Tom Copeland, Sep 07 2016: (Start)
The row partition polynomials of this array P(n,x_1,x_2,...,x_n), given in the Lang link, are n! * S(n,x_1,x_2,...,x_n), where S(n,x_1,...,x_n) are the elementary Schur polynomials, for which d/d(x_m) S(n,x_1,...,x_n) = S(n-m,x_1,...,x_(n-m)) with S(k,...) = 0 for k < 0, so d/d(x_m) P(n,x_1,...,x_n) = (n!/(n-m)!) P(n-m,x_1,...,x_(n-m)), confirming that the row polynomials form an Appell sequence in the indeterminate x_1 with P(0,...) = 1. See p. 127 of the Ernst paper for more on these Schur polynomials.
With the e.g.f. exp[t * P(.,x_1,x_2,..)] = exp(t*x_1) * exp(x_2 t^2 + x_3 t^3 + ...), the e.g.f. for the partition polynomials that form the umbral compositional inverse sequence U(n,x_1,...,x_n) in the indeterminate x_1 is exp[t * U(.,x_1,x_2,...)] = exp(t*x_1) exp[-(x_2 t^2 + x_3 t^3 + ...)]; therefore, U(n,x_1,x_2,...,x_n) = P(n,x_1,-x_2,.,-x_n), so umbrally P[n,P(.,x_1,-x_2,-x_3,...),x_2,x_3,...,x_n] = (x_1)^n = P[n,P(.,x_1,x_2,...),-x_2,-x_3,...,-x_n]. For example, P(1,x_1) = x_1, P2(x_1,x_2) = 2 x_2 + x_1^2, and P(3,x_1,x_2,x_3) = 6 x_3 + 6 x_2 x_1 + x_1^3, then P[3,P(.,x_1,-x_2,...),x_2,x_3] = 6 x_3 + 6 x_2 P(1,x_1) + P(3,x_1,-x_2,-x_3) = 6 x_3 + 6 x_2 x_1 + 6 (-x_3) + 6 (-x_2) x_1 + x_1^3 = x_1^3.
From the Appell formalism, umbrally [P(.,0,x_2,x_3,...) + y]^n = P(n,y,x_2,x_3,...,x_n).
The indeterminates of the partition polynomials can also be extracted using the Faber polynomials of A263916 with -n * x_n = F(n,S(1,x_1),...,S(n,x_1,...,x_n)) = F(n,P(1,x_1),...,P(n,x_1,...,x_n)/n!). Compare with A263634.
Also P(n,x_1,...,x_n) = ST1(n,x_1,2*x_2,...,n*x_n), where ST1(n,...) are the row partition polynomials of A036039.
(End)

Extensions

Name augmented by Tom Copeland, Dec 08 2022

A263634 Irregular triangle read by rows: row n gives coefficients of n-th logarithmic polynomial L_n(x_1, x_2, ...) with monomials sorted into standard order.

Original entry on oeis.org

1, -1, 1, 2, -3, 1, -6, 12, -4, -3, 1, 24, -60, 20, 30, -5, -10, 1, -120, 360, -120, -270, 30, 120, 30, -6, -15, -10, 1, 720, -2520, 840, 2520, -210, -1260, -630, 42, 210, 140, 210, -7, -21, -35, 1
Offset: 1

Views

Author

N. J. A. Sloane, Oct 29 2015

Keywords

Comments

"Standard order" here means as produced by Maple's "sort" command.
From Petros Hadjicostas, May 27 2020: (Start)
According to the Maple help files for the "sort" command, polynomials in multiple variables are "sorted in total degree with ties broken by lexicographic order (this is called graded lexicographic order)."
Thus for example, x_1^2*x_3 = x_1*x_1*x_3 > x_1*x_2*x_2 = x_1*x_2^2, while x_1^2*x_4 = x_1*x_1*x_4 > x_1*x_2*x_3. (End)
Row sums are 0 (for n > 1). Numbers of terms in rows are partition numbers A000041.
From Tom Copeland, Nov 06 2015: (Start)
With the formal Taylor series f(x) = 1 + x[1] x + x[2] x^2/2! + ... , the partition polynomials of this entry give d[log(f(x))]/dx = L_1(x[1]) + L_2(x[1], x[2]) x + L_3(...) x^2/2! + ..., and the coefficients of the reduced polynomials with x[n] = t are signed A028246.
The raising operator R = x + d[log(f(D)]/dD = x + L_1(x[1]) + L_2[x[1], x[2]) D + L_3(x[1], x[2], x[3]) D^2/2! + ... with D = d/dx generates an Appell sequence of polynomials, given umbrally by P_n(x[1], ..., x[n]; x) = (x[.] + x)^n = Sum_{k=0..n} binomial(n,k) x[k] * x^(n-k) = R^n 1 with the e.g.f. f(t)*e^(x*t) = exp[t P.(x[1], ..., x[.]; x)]. P_0 = x[0] = 1.
The umbral compositional inverse Appell sequence is generated by R = x - d[log(f(D))]/dD with e.g.f. e^(x*t)/f(t) = exp[t IP.(x[1], ..., x[.]; x)], so umbrally IP_n(x[1], ..., x[n]; P.(x[1], ..., x[n]; x)) = x^n = P_n(x[1], ..., x[n]; IP.(x[1], ..., x[n]; x)). An unsigned array for the reduced IP_n(x[1], ..., x[n]; x) polynomials with IP_0 = x[0] = 1 and x[n] = -1 for n > 0 is A154921, for which f(t) = 2 - e^t. (End)
From Tom Copeland, Sep 08 2016: (Start)
The Appell formalism allows a matrix representation in the power basis x^n of the raising operator R that incorporates this array's partition polynomials L_n(x[1], ..., x[n]):
VP_(n+1) = VP_n * R = VP_n * XPS^(-1) * MX * XPS, where XPS is the matrix formed from multiplying the n-th diagonal of the Pascal matrix PS of A007318 by the indeterminate x[n], with x[0] = 1 for the main diagonal of ones, i.e., XPS[n,k] = PS[n,k] * x[n-k]; the matrix MX is A129185; the matrix XPS^(-1) is the inverse of XPS, which can be formed by multiplying the diagonals of the Pascal matrix by the partition polynomials IPT(n, x[1], ..., x[n]) of A133314, i.e., XPS^(-1)[n,k] = PS[n,k] * IPT(n-k, x[1], ...); and VP_n is the row vector in the power basis representing the Appell polynomial P_n(x) formed from the basic sequence of moments 1, x[1], x[2], ..., i.e., umbrally P_n(x) = (x[.] + x)^n = Sum_{k=0..n} binomial(n,k) * x[k] * x^(n-k).
Then R = XPS^(-1) * MX * XPS is the Pascal matrix PS with an additional first superdiagonal of ones and the other lower diagonals multiplied by the partition polynomials of this array, i.e., R[n,k] = PS[n,k] * L_{n+1-k}(x[1], ..., x[n+1-k]) except for the first superdiagonal of ones.
Consistently, VP_n = (1, 0, 0, ...) * R^n = (1, 0, 0, ...) * XPS^(-1) * MX^n * XPS = (1, 0, 0, ...) * MX^n * XPS = the n-th row vector of XPS, which is the vector representation of P_n(x) = (x[.] + x)^n with x[0] = 1.
See the Copeland link for the umbral representation R = exp[g.*D] * x * exp[h.*D] that reflects the matrix representations.
The Stirling partition polynomials of the first kind St1_n(a[1], a[2], ..., a[n]) of A036039, the Stirling partition polynomials of the second kind St2_n(b[1], b[2], ..., b[n]) of A036040, and the refined Lah polynomials Lah_n[c[1], c[2], ..., c[n]) of A130561 are Appell sequences in the respective distinguished indeterminates a[1], b[1], and c[1]. Comparing the formulas for their raising operators with that in this entry, L_n(x[1], x[2], ..., x[n]) evaluates to
A) (n-1)! * a[n] for x[n] = St1_n(a[1], a[2], ..., a[n]);
B) b[n] for x[n] = St2_n(b[1], b[2], ..., b[n]);
C) n! * c[n] for x[n] = Lah_n(c[1], c[2], ..., c[n]).
Conversely, from the respective e.g.f.s (added Sep 12 2016)
D) x[n] = St1_n(L_1(x[1])/0!, ..., L_n(x[1], ..., x[n])/(n-1)!);
E) x[n] = St2_n(L_1(x[1]), ..., L_n(x[1], ..., x[n]));
F) x[n] = Lah_n(L_1(x[1])/1!, ..., L_n(x[1], ..., x[n])/n!).
Given only the Appell sequence with no closed form for the e.g.f., the raising operator can be generated using this formalism, as has been partially done for A134264. (End)
For the Appell sequences above, the raising operator is related to the recursion P_(n+1)(x) = x * P_n(x) + Sum_{k=0..n} binomial(n,k) * L_(n-k+1)(x[1], ..., x[n+k-1]) * P_k(x). For a derivation and connections to formal cumulants (c_n = L_n(x[1], ...)) and moments (m_n = x[n]), see the Copeland link on noncrossing partitions. With x = 0, the recursion reduces to x[n+1] = Sum_{k = 0..n} binomial(n,k) * L_(n-k+1)(x[1], ..., x[n+k-1]) * x[k] with x[0] = 1. This array is a differently ordered version of A127671. - Tom Copeland, Sep 13 2016
With x[n] = x^(n-1), a signed version of A130850 is obtained. - Tom Copeland, Nov 14 2016
See p. 2 of Getzler for a relation to stable graphs called necklaces used in computations for Deligne-Mumford-Knudsen moduli spaces of stable curves of genus 1. - Tom Copeland, Nov 15 2019
For a relation to a combinatorial Faa di Bruno Hopf algebra related to functional composition, as presented by Connes and Moscovici, see Figueroa et al. - Tom Copeland, Jan 17 2020
From Tom Copeland, May 17 2020: (Start)
The e.g.f. of an Appell sequence is f(t) e^(x*t) with f(0) = 1. Given the Laguerre-Polya class function f(t) = e^(-a*t^2 + b*t) Product_m (1 - t/z_m) e^(t/z_m) with a = 0 for simplicity (more generally a >= 0) and b real and where the product runs over all the zeros z_m of f(t) with all zeros real and Sum_m 1/(z_m)^2 convergent, the raising operator of the Appell polynomials is R = x + b - Sum_{k > 0} c_(k+1) D^k with c_k = Sum_m (1/(z_m)^k), i.e., traces of powers of the reciprocals of the zeros. From R in earlier comments, b = L_1(x_1) and otherwise c_k = -L_k(x_1, ..., x_k).
The Laguerre / Turan / de Gua inequalities (Csordas and Williamson and Skovgaard) imply that all the zeros of each Appell polynomial are real and simple and its extrema are local maxima above the x-axis and local minima below and are located above or below the zeros of the next lower degree Appell polynomial. (End)
From Tom Copeland, Oct 15 2020: (Start)
With a_n = n! * b_n = (n-1)! * c_n for n > 0, represent a function with f(0) = a_0 = b_0 = 1 as an
A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!
B) ordinary generating function (o.g.f.), or formal power series: f(x) = 1/(1-b.x) = 1 + Sum_{n > 0} b_n * x^n
C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0} c_n * x^n /n.
Expansions of log(f(x)) are given in
I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] = e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials
II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] = log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.
Expansions of exp(f(x)-1) are given in
III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind
IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials
V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] = e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.
Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)
Ignoring signs, these polynomials appear in Schröder in the set of equations (II) on p. 343 and in Stewart's translation on p. 31. - Tom Copeland, Aug 25 2021

Examples

			The first few polynomials are:
(1) x[1].
(2) -x[1]^2 + x[2].
(3) 2*x[1]^3 - 3*x[1]*x[2] + x[3].
(4) -6*x[1]^4 + 12*x[1]^2*x[2] - 4*x[1]*x[3] - 3*x[2]^2 + x[4].
(5) 24*x[1]^5 - 60*x[1]^3*x[2] + 20*x[1]^2*x[3] + 30*x[1]*x[2]^2 - 5*x[1]*x[4] - 10*x[2]*x[3] + x[5].
(6) -120*x[1]^6 + 360*x[1]^4*x[2] - 120*x[1]^3*x[3] - 270*x[1]^2*x[2]^2 + 30*x[1]^2*x[4] + 120*x[1]*x[2]*x[3] + 30*x[2]^3 - 6*x[1]*x[5] - 15*x[2]*x[4] - 10*x[3]^2 + x[6].
...
[1]    1
[2]   -1,    1
[3]    2,   -3,     1
[4]   -6,   12,    -4,    -3,   1
[5]   24,  -60,    20,    30,  -5,  -10,   1
[6] -120,  360,  -120,  -270,  30,  120,  30, -6, -15, -10, 1
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, pp. 140, 156, 308.

Crossrefs

Programs

  • Maple
    triangle := proc(numrows) local E, s, Q;
    E := add(x[i]*t^i/i!, i=1..numrows);
    s := series(log(1 + E), t, numrows+1);
    Q := k -> sort(expand(k!*coeff(s, t, k)));
    seq(print(coeffs(Q(k))), k=1..numrows) end:
    triangle(6); # updated by Peter Luschny, May 27 2020

Formula

G.f.: Log(1 + Sum_{i >= 1} x_i*t^i/i!) = Sum_{n >= 1} L_n(x_1, x_2, ...)*t^n/n!. [Comtet, p. 140, Eq. [5a]. - corrected by Tom Copeland, Sep 08 2016]
Conjecture: row polynomials are R(n,1) for n > 0 where R(n,k) = R(n-1,k+1) - Sum_{j=1..n-1} binomial(n-2,j-1)*R(j,k)*R(n-j,1) for n > 1, k > 0 with R(1,k) = x_k for k > 0. - Mikhail Kurkov, Mar 30 2025

A141618 Triangle read by rows: number of nilpotent partial transformations (of an n-element set) of height r (height(alpha) = |Im(alpha)|), 0 <= r < n.

Original entry on oeis.org

1, 1, 2, 1, 9, 6, 1, 28, 72, 24, 1, 75, 500, 600, 120, 1, 186, 2700, 7800, 5400, 720, 1, 441, 12642, 73500, 117600, 52920, 5040, 1, 1016, 54096, 571536, 1764000, 1787520, 564480, 40320, 1, 2295, 217800, 3916080, 21019824, 40007520, 27941760, 6531840, 362880, 1, 5110, 839700, 24555600, 214326000
Offset: 1

Views

Author

Abdullahi Umar, Aug 23 2008

Keywords

Comments

The sum of each row of the sequence (as a triangular array) is A000272. Second left-downward diagonal is A058877.
From Tom Copeland, Oct 26 2014: (Start)
With T(x,t) the e.g.f. for A055302 for the number of labeled rooted trees with n nodes and k leaves, the mirror of the row polynomials of this array are given by e^T(x,t) = exp[ t * x + (2t) * x^2/2! + (6t + 3t^2) * x^3/3! + ...] = 1 + t * x + (2t + t^2) * x^2/2! + (6t + 9t^2 + t^3) * x^3/3! + ... = 1 + Nr(x,t).
Equivalently, e^x-1 = Nr[Tinv(x,t),t] = t * N[t*Tinv(x,t),1/t], where N(x,t) is the e.g.f. of this array and Tinv(x,t) is the comp. inverse in x of T(x,t). Note that Nr(x,t) = t * N(x*t,1/t), and N(x,t) = t * Nr(x*t,1/t). Also, log[1 + Nr(x,t)]= x * [t + Nr(x,t)] = T(x,t).
E.g.f. is N(x,t)= t * {exp[T(x*t,1/t)] - 1}, and log[1 + N(x,t)/t] = T(x*t,1/t) = x + (2t) * x^2/2! + (3t + 6t^2) * x^3/3! + (4t + 36t^2 + 24t^3) * x^4/4! + ... = x + (t) * x^2 + (t + 2t^2) * x^3/2! + (t + 9t^2 + 6t^3) * x^4/3! + ... is the comp. inverse in x of x / [1 + t * (e^x - 1)].
The exp/log transforms (A036040/A127671) generally give associations between enumerations of sets of connected graphs/objects (in this case, trees) and sets of disconnected (or not necessarily connected) graphs/objects (in this case, bipartite graphs of the nilpotent transformations). The transforms also relate formal cumulants and moments so that Nr(x,t) is then the e.g.f. for the formal moments associated to the formal cumulants whose e.g.f. is T(x,t). (End)
T(n,k) is the number of parking functions of length n containing exactly k+1 distinct values in its image. - Alan Kappler, Jun 08 2024

Examples

			N(J(4,2)) = 6*6*2 = 72.
From _Peter Bala_, Oct 22 2008: (Start)
Triangle begins
n\k|..0.....1.....2.....3.....4....5
=====================================
.1.|..1
.2.|..1.....2
.3.|..1.....9.....6
.4.|..1....28....72....24
.5.|..1....75...500...600...120
.6.|..1...186..2700..7800..5400...720
...
(End)
		

Crossrefs

Programs

  • Maple
    A048993 := proc(n,k)
        combinat[stirling2](n,k) ;
    end proc:
    A141618 := proc(n,k)
        binomial(n,k)*k!*A048993(n,k+1) ;
    end proc:
  • Mathematica
    Flatten[CoefficientList[CoefficientList[InverseSeries[Series[Log[1 + x]/(1 + t*x),{x,0,9}]],x]*Table[n!, {n,0,9}],t]] (* Peter Luschny, Oct 24 2015, after Peter Bala *)
  • PARI
    A055302(n,k)=n!/k!*stirling(n-1, n-k,2);
    T(n,k)=A055302(n+1,n+1-k) / (n+1);
    for(n=1,10,for(k=1,n,print1(T(n,k),", "));print());
    \\ Joerg Arndt, Oct 27 2014

Formula

N(J(n,r)) = C(n,r)*S(n,r+1)*r! where S(n, r + 1) is a Stirling number of the second kind (given by A048993 with zeros removed); generating function = (x+1)^(n-1).
From Peter Bala, Oct 22 2008: (Start)
Define a functional I on formal power series of the form f(x) = 1 + ax + bx^2 + ... by the following iterative process. Define inductively f^(1)(x) = f(x) and f^(n+1)(x) = f(x*f^(n)(x)) for n >= 1. Then set I(f(x)) = lim_{n -> infinity} f^(n)(x) in the x-adic topology on the ring of formal power series; the operator I may also be defined by I(f(x)) := 1/x*series reversion of x/f(x).
Let f(x) = 1 + a*x + a*x^2/2! + a*x^3/3! + ... . Then the e.g.f. for this table is I(f(x)) = 1 + a*x +(a + 2*a^2)*x^2/2! + (a + 9*a^2 + 6*a^3)*x^3/3! + (a + 28*a^2 + 72*a^3 + 24*a^4)*x^4/4! + ... . Note, if we take f(x) = 1 + a*x + a*x^2 + a*x^3 + ... then I(f(x)) is the o.g.f. of the Narayana triangle A001263. (End)
A generator for this array is given by the inverse, g(x,t), of f(x,t)= x/(1 + t * (e^x-1)). Then A248927 gives h(x,t)= x / f(x,t) = 1 + t*(e^x-1)= 1 + t * (x + x^2/2! + x^3/3! + ...) and g(x,t)= x * (1 + t * x + (t + 2 t^2) * x^2/2! + (t + 9 t^2 + 6 t^3) * x^3/3! + ...), so by Bala's arguments A248927 is a refinement of A141618 with row sums A000272. The connection to Narayana numbers is reflected in the relation between A248927 and A134264. See A145271 for more relations that g(x,t) and f(x,t) must satisfy. - Tom Copeland, Oct 17 2014
T(n,k) = C(n,k-1) * A028246(n,k) = C(n,k-1) * A019538(n,k)/k = A055302(n+1,n+1-k) / (n+1). - Tom Copeland, Oct 25 2014
E.g.f. is the series reversion of log(1 + x)/(1 + t*x) with respect to x. Cf. A198204. - Peter Bala, Oct 21 2015

Extensions

More terms from Joerg Arndt, Oct 27 2014

A249548 Coefficients of reduced partition polynomials of A134264 for computing Lagrange compositional inversion.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 1, 5, 1, 6, 3, 5, 1, 7, 7, 21, 1, 8, 8, 4, 28, 28, 14, 1, 9, 9, 9, 36, 72, 12, 84, 1, 10, 10, 10, 5, 45, 90, 45, 45, 120, 180, 42, 1, 11, 11, 11, 11, 55, 110, 110, 55, 55, 165, 495, 165, 330, 1, 12, 12, 12, 12, 6, 66, 132, 132, 66, 66, 132, 22, 220, 660, 330, 660, 55, 495, 990, 132
Offset: 0

Views

Author

Tom Copeland, Oct 31 2014

Keywords

Comments

Coefficients of reduced partition polynomials of A134264 for computing the complete partition polynomials for the Lagrange compositional inversion of A134264 (see Oct 2014 comment by Copeland there). Umbrally,
e^(x*t) * exp[Prt(.;1,0,h_2,..) * t] = exp[Prt(.;1,x,h_2,..) * t], where Prt(n;1,0,h_2,..,h_n) are the reduced (h_0 = 1 and h_1 = 0) partition polynomials of the complete polynomials Prt(n;h_0,h_1,h_2,..,h_n) of A134264.
Partitions are given in the order of those on page 831 of Abramowitz and Stegun. Formulas for the coefficients of the partitions are given in A134264.
Row sums are the Motzkin sums or Riordan numbers A005043. - Tom Copeland, Nov 09 2014
From Tom Copeland, Jul 03 2018: (Start)
The matrix and operator formalism for Sheffer Appell sequences leads to the following relations with D = d/dh_1.
Exp[Prt(.;1,0,h_2,..) * D] (h_1)^n = [h_1 + Prt(.;1,0,h_2,...)]^n = Prt(n;1,h_1,h_2,...), the partition polynomials of A134264 for g(t)/t with h_0 = 1.
For the umbral compositional inverses described below,
Exp[UPrt(.;1,0,h_2,..) * D] (h_1)^n = [h_1 + UPrt(.;1,0,h_2,...)]^n = UPrt(n;1,h_1,h_2,...).
The respective e.g.f.s are multiplicative inverses; that is, exp[Prt(.;1,0,h_2,..) * t] = 1/exp[UPrt(.;1,0,h_2,..) * t], so the formalism of A133314 applies.
The raising operator R such that R Prt(n;1,h_1,h_2,...) = Prt(n+1;1,h_1,h_2,...) is R = exp[Prt(.;1,0,h_2,...)*D] h_1 exp[UPrt(.;1,0,h_2,..)*D] since R Prt(n+1;1,h_1,h_2,...) = exp[Prt(.;1,0,h_2,...)*D] h_1 (h_1)^n = Prt(n+1;h_1,h_2,...) from the definition of the umbral compositional inverse. This may also be expressed as R = h_1 + d/dD log[exp[Prt(.;1,0,h_2,...) * D]], so, using A127671, R = h_1 + h_2 D + h_3 D^2/2! + (h_4 - h_2^2) D^3/3! + (h_5 - 5 h_2 h_3) D^4/4! + (h_6 - 9 h_2 h_4 + 5 h_2^3 - 7 h_3^2) D^5/5! + (h_7 - 28 h_3 h_4 - 14 h_2 h_5 + 56 h_2^2 h_3) D^6/6! + ... . (End)

Examples

			Prt(0) = 1
Prt(1;1,0) = 0
Prt(2;1,0,h_2) = 1 h_2
Prt(3;1,0,h_2,h_3) = 1 h_3
Prt(4;1,0,h_2,..,h_4) = 1 h_4 + 2 (h_2)^2
Prt(5;1,0,h_2,..,h_5) = 1 h_5 + 5 h_2 h_3
Prt(6;1,0,h_2,..,h_6) = 1 h_6 + 6 h_2 h_4 + 3 (h_3)^2 + 5 (h_2)^3
Prt(7;1,0,h_2,..,h_7) = 1 h_7 + 7 h_3 h_4 + 7 h_2 h_5 + 21 h_3 (h_2)^2
...
------------
With h_n denoted by (n'), the first seven partition polynomials of A134264 with h_0=1 are given by the first seven coefficients of the truncated Taylor series expansion of the Euler binomial transform
e^[(1') * t] * {1 + 1 (2') *  t^2/2! + 1 (3') *  t^3/3! + [1 (4') + 2 (2')^2] *  t^4/4! + [1 (5') + 5 (2')(3')] *  t^5/5! + [1 (6') + 6 (2')(4') + 3 (3')^2 + 5 (2')^3] *  t^6/6!}, giving the truncated expansion
1 + 1 (1') * t + [1 (2') + 1 (1')^2] * t^2/2! + ... + [1 (6') + 6 (1')(5') + 6 (2')(4') + 3 (3')^2 + 15 (1')^2(4') + 30 (1')(2')(3') + 5 (2')^3 + 20 (1')^3(3') + 30 (1')^2(2')^2 + 15 (1')^4(2') + 1 (1')^6] * t^6/6!.
Extending the number of reduced partition polynomials of the transform allows for further complete polynomials of A134264 to be computed.
		

Crossrefs

Cf. A127671.
Rows lengths are given by A002865 (except for row 1).

Programs

  • Mathematica
    rows[n_] := {{1}, {0}}~Join~Module[
        {g = 1 / D[t / (1 + Sum[h[k] t^k, {k, 2, n}] + O[t]^(n+1)), t], p = t, r},
        r = Reap[Do[p = g D[p, t]/k; Sow[Expand[Normal@p /. {t -> 0}]], {k, n+1}]][[2, 1, 2 ;;]];
        Table[Coefficient[r[[k]], Product[h[t], {t, p}]], {k, 2, n}, {p, Sort[Sort /@ IntegerPartitions[k, k, Range[2, k]]]}]];
    rows[12] // Flatten (* Andrey Zabolotskiy, Feb 18 2024 *)

Formula

From Tom Copeland, Nov 10 2014: (Start)
Terms may be computed symbolically up to order n by using an iterated derivative evaluated at t=0:
with g(t) = 1/{d/dt [t/(1 + 0 t + h_2 t^2 + h_3 t^3 + ... + h_n t^n)]},
evaluate 1/n! * [g(t) d/dt]^n t at t=0, i.e., ask a symbolic math app for the first term in a series expansion of this iterated derivative, to obtain Prt(n-1).
Alternatively, the explicit formula in A134264 for the numerical coefficients of each partition can be used. (End)
From Tom Copeland, Nov 12 2014: (Start)
The first few partitions polynomials formed by taking the reciprocal of the e.g.f. of this entry's e.g.f. (cf. A133314) are
UPrt(0) = 1
UPrt(1;1,0) = 0
UPrt(2;1,0,h_2) = -1 h_2
UPrt(3;1,0,h_2,h_3) = -1 h_3
UPrt(4;1,0,h_2,..,h_4) = -1 h_4 + 4 (h_2)^2
UPrt(5;1,0,h_2,..,h_5) = -1 h_5 + 15 h_2 h_3
UPrt(6;1,0,h_2,..,h_6) = -1 h_6 + 24 h_2 h_4 + 17 (h_3)^2 + -35 (h_2)^3
...
Therefore, umbrally, [Prt(.;1,0,h_2,...) + UPrt(.;1,0,h_2,...)]^n = 0 for n>0 and unity for n=0.
Example of the umbral operation:
(a. + b.)^2 = a.^2 + 2 a.* b. + b.^2 = a_2 + 2 a_1 * b_1 + b_2.
This implies that the umbral compositional inverses (see below) of the partition polynomials of the Lagrange inversion formula (LIF) of A134264 with h_0=1 are given by UPrt(n;1,h_1,h_2,...,h_n) = [UPrt(.;1,0,h_2,...,h_n) + h_1]^n, so that the sequence of polynomials UPrt(n;1,h_1,h_2,...,h_n) is an Appell sequence in the indeterminate h_1. So, if one calculates UPrt(n;1,h_1,...,h_n), the lower order UPrt(n-1;1,h_1,...,h_(n-1)) can be found by taking the derivative w.r.t. h_1 and dividing by n. Same applies for Prt(n;1,h_1,h_2,...,h_n).
This connects the combinatorics of the permutohedra through A133314 and A049019, or their duals, to the noncrossing partitions, Dyck lattice paths, etc. that are isomorphic with the LIF of A134264.
An Appell sequence P(.,x) with the e.g.f. e^(x*t)/w(t) possesses an umbral inverse sequence UP(.,x) with the e.g.f. w(t)e^(x*t), i.e., polynomials such that P(n,UP(.,x))= x^n = UP(n,P(.,x)) through umbral substitution, as in the binomial example. The Bernoulli polynomials with w(t) = t/(e^t - 1) are a good example with the umbral compositional inverse sequence UP(n,x) = [(x+1)^(n+1)-x^(n+1)] / (n+1) (cf. A074909 and A135278). (End)

Extensions

Formula for Prt(7,..) and a(12)-a(15) added by Tom Copeland, Jul 22 2016
Rows 8-12 added by Andrey Zabolotskiy, Feb 18 2024

A274805 The logarithmic transform of sigma(n).

Original entry on oeis.org

1, 2, -3, -6, 45, 11, -1372, 4298, 59244, -573463, -2432023, 75984243, -136498141, -10881169822, 100704750342, 1514280063802, -36086469752977, -102642110690866, 11883894518252419, -77863424962770751, -3705485804176583500, 71306510264347489177
Offset: 1

Views

Author

Johannes W. Meijer, Jul 27 2016

Keywords

Comments

The logarithmic transform [LOG] transforms an input sequence b(n) into the output sequence a(n). The LOG transform is the inverse of the exponential transform [EXP], see the Weisstein link and the Sloane and Plouffe reference. This relation goes by the name of Riddell’s formula. For information about the EXP transform see A274804. The logarithmic transform is related to the inverse multinomial transform, see A274844 and the first formula.
The definition of the LOG transform, see the second formula, shows that n >= 1. To preserve the identity EXP[LOG[b(n)]] = b(n) for n >= 0 for a sequence b(n) with offset 0 the shifted sequence b(n-1) with offset 1 has to be used as input for the LOG transform, otherwise information about b(0) will be lost in transformation.
In the a(n) formulas, see the examples, the cumulant expansion numbers A127671 appear.
We observe that the logarithmic transform leaves the value of a(0) undefined.
The Maple programs can be used to generate the logarithmic transform of a sequence. The first program uses a formula found by Alois P. Heinz, see A001187 and the first formula. The second program uses the definition of the logarithmic transform, see the Weisstein link and the second formula. The third program uses information about the inverse of the logarithmic transform, see A274804.

Examples

			Some a(n) formulas, see A127671:
a(0) = undefined
a(1) = 1*x(1)
a(2) = 1*x(2) - x(1)^2
a(3) = 1*x(3) - 3*x(1)*x(2) + 2*x(1)^3
a(4) = 1*x(4) - 4*x(1)*x(3) - 3*x(2)^2 + 12*x(1)^2*x(2) - 6*x(1)^4
a(5) = 1*x(5) - 5*x(1)*x(4) - 10*x(2)*x(3) + 20*x(1)^2*x(3) + 30*x(1)*x(2)^2 - 60*x(1)^3*x(2) + 24*x(1)^5
		

References

  • Frank Harary and Edgar M. Palmer, Graphical Enumeration, 1973.
  • Robert James Riddell, Contributions to the theory of condensation, Dissertation, University of Michigan, Ann Arbor, 1951.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, 1995, pp. 18-23.

Crossrefs

Some LOG transform pairs are, n >= 1: A006125(n-1) and A033678(n); A006125(n) and A001187(n); A006125(n+1) and A062740(n); A000045(n) and A112005(n); A000045(n+1) and A007553(n); A000040(n) and A007447(n); A000051(n) and (-1)*A263968(n-1); A002416(n) and A062738(n); A000290(n) and A033464(n-1); A029725(n-1) and A116652(n-1); A052332(n) and A002031(n+1); A027641(n)/A027642(n) and (-1)*A060054(n+1)/(A075180(n-1).

Programs

  • Maple
    nmax:=22: with(numtheory): b := proc(n): sigma(n) end: a:= proc(n) option remember; b(n) - add(k*binomial(n, k)*b(n-k)*a(k), k=1..n-1)/n: end: seq(a(n), n=1..nmax); # End first LOG program.
    nmax:=22: with(numtheory): b := proc(n): sigma(n) end: t1 := log(1 + add(b(n)*x^n/n!, n=1..nmax+1)): t2 := series(t1, x, nmax+1): a := proc(n): n!*coeff(t2, x, n) end: seq(a(n), n=1..nmax); # End second LOG program.
    nmax:=22: with(numtheory): b := proc(n): sigma(n) end: f := series(exp(add(r(n)*x^n/n!, n=1..nmax+1)), x, nmax+1): d := proc(n): n!*coeff(f, x, n) end: a(1):=b(1): r(1):= b(1): for n from 2 to nmax+1 do r(n) := solve(d(n)-b(n), r(n)): a(n):=r(n): od: seq(a(n), n=1..nmax); # End third LOG program.
  • Mathematica
    a[1] = 1; a[n_] := a[n] = DivisorSigma[1, n] - Sum[k*Binomial[n, k] * DivisorSigma[1, n-k]*a[k], {k, 1, n-1}]/n; Table[a[n], {n, 1, 22}] (* Jean-François Alcover, Feb 27 2017 *)
  • PARI
    N=33; x='x+O('x^N); Vec(serlaplace(log(1+sum(n=1,N,sigma(n)*x^n/n!)))) \\ Joerg Arndt, Feb 27 2017

Formula

a(n) = b(n) - Sum_{k = 1..n-1}((k*binomial(n, k)*b(n-k)*a(k))/n), n >= 1, with b(n) = A000203(n) = sigma(n).
E.g.f. log(1+ Sum_{n >= 1}(b(n)*x^n/n!)), n >= 1, with b(n) = A000203(n) = sigma(n).

A274844 The inverse multinomial transform of A001818(n) = ((2*n-1)!!)^2.

Original entry on oeis.org

1, 8, 100, 1664, 34336, 843776, 24046912, 779780096, 28357004800, 1143189536768, 50612287301632, 2441525866790912, 127479926768287744, 7163315850315825152, 431046122080208896000, 27655699473265974050816, 1884658377677216933085184
Offset: 1

Views

Author

Johannes W. Meijer, Jul 27 2016

Keywords

Comments

The inverse multinomial transform [IML] transforms an input sequence b(n) into the output sequence a(n). The IML transform inverses the effect of the multinomial transform [MNL], see A274760, and is related to the logarithmic transform, see A274805 and the first formula.
To preserve the identity MNL[IML[b(n)]] = b(n) for n >= 0 for a sequence b(n) with offset 0 the shifted sequence b(n-1) with offset 1 has to be used as input for the MNL.
In the a(n) formulas, see the examples, the cumulant expansion numbers A127671 appear.
We observe that the inverse multinomial transform leaves the value of a(0) undefined.
The Maple programs can be used to generate the inverse multinomial transform of a sequence. The first program is derived from a formula given by Alois P. Heinz for the logarithmic transform, see the first formula and A001187. The second program uses the e.g.f. for multivariate row polynomials, see A127671 and the examples. The third program uses information about the inverse of the inverse of the multinomial transform, see A274760.
The IML transform of A001818(n) = ((2*n-1)!!)^2 leads quite unexpectedly to A005411(n), a sequence related to certain Feynman diagrams.
Some IML transform pairs, n >= 1: A000110(n) and 1/A000142(n-1); A137341(n) and A205543(n); A001044(n) and A003319(n+1); A005442(n) and A000204(n); A005443(n) and A001350(n); A007559(n) and A000244(n-1); A186685(n+1) and A131040(n-1); A061711(n) and A141151(n); A000246(n) and A000035(n); A001861(n) and A141044(n-1)/A001710(n-1); A002866(n) and A000225(n); A000262(n) and A000027(n).

Examples

			Some a(n) formulas, see A127671:
a(0) = undefined
a(1) = (1/0!) * (1*x(1))
a(2) = (1/1!) * (1*x(2) - x(1)^2)
a(3) = (1/2!) * (1*x(3) - 3*x(2)*x(1) + 2*x(1)^3)
a(4) = (1/3!) * (1*x(4) - 4*x(3)*x(1) - 3*x(2)^2 + 12*x(2)*x(1)^2 - 6*x(1)^4)
a(5) = (1/4!) * (1* x(5) - 5*x(4)*x(1) - 10*x(3)*x(2) + 20*x(3)*x(1)^2 + 30*x(2)^2*x(1) -60*x(2)*x(1)^3 + 24*x(1)^5)
		

References

  • Richard P. Feynman, QED, The strange theory of light and matter, 1985.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, 1995, pp. 18-23.

Crossrefs

Programs

  • Maple
    nmax:=17: b := proc(n): (doublefactorial(2*n-1))^2 end: c:= proc(n) option remember; b(n) - add(k*binomial(n, k)*b(n-k)*c(k), k=1..n-1)/n end: a := proc(n): c(n)/(n-1)! end: seq(a(n), n=1..nmax); # End first IML program.
    nmax:=17: b := proc(n): (doublefactorial(2*n-1))^2 end: t1 := log(1+add(b(n)*x^n/n!, n=1..nmax+1)): t2 := series(t1, x, nmax+1): a := proc(n): n*coeff(t2, x, n) end: seq(a(n), n=1..nmax); # End second IML program.
    nmax:=17: b := proc(n): (doublefactorial(2*n-1))^2 end: f := series(exp(add(t(n)*x^n/n, n=1..nmax)), x, nmax+1): d := proc(n): n!*coeff(f, x, n) end: a(1):=b(1): t(1):= b(1): for n from 2 to nmax+1 do t(n) := solve(d(n)-b(n), t(n)): a(n):=t(n): od: seq(a(n), n=1..nmax); # End third IML program.
  • Mathematica
    nMax = 22; b[n_] := ((2*n-1)!!)^2; c[n_] := c[n] = b[n] - Sum[k*Binomial[n, k]*b[n-k]*c[k], {k, 1, n-1}]/n; a[n_] := c[n]/(n-1)!; Table[a[n], {n, 1, nMax}] (* Jean-François Alcover, Feb 27 2017, translated from Maple *)

Formula

a(n) = c(n)/(n-1)! with c(n) = b(n) - Sum_{k=1..n-1}(k*binomial(n, k)*b(n-k)*c(k)), n >= 1 and a(0) = undefined, with b(n) = A001818(n) = ((2*n-1)!!)^2.
a(n) = A000079(n-1) * A005411(n), n >= 1.

A350499 Unsigned coefficients of free moment partition polynomials determining the free cumulants from the free moments of free probability theory. Irregular triangle with row lengths given by A000041, n >= 1.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 4, 2, 10, 5, 1, 5, 5, 15, 15, 35, 14, 1, 6, 6, 3, 21, 42, 7, 56, 84, 126, 42, 1, 7, 7, 7, 28, 56, 28, 28, 84, 252, 84, 210, 420, 462, 132, 1, 8, 8, 8, 4, 36, 72, 72, 36, 36, 120, 360, 180, 360, 30, 330, 1320, 660, 792, 1980, 1716, 429
Offset: 1

Views

Author

Tom Copeland, Jan 01 2022

Keywords

Comments

Coefficients are listed in Abramowitz and Stegun order (A036036).
Irregular triangular matrix of the unsigned coefficients of the free moment partition polynomials of free probability theory, for a single variable, that give the free formal cumulants given the free formal moments. This set of partition polynomials together with those of A134264 are the counterparts to the exp-log relations for the classical formal moments and cumulants associated with A036040 and A127671.
Associations with a compositional inverse pair of Laurent series, Kac-Schwarz operators of 2-D quantum theory, Virasoro / Witt / Heisenberg group actions, and KP and KdV integrable hierarchies are noted in references supplied in the MathOverflow link as well as a geometric combinatorial model based upon noncrossing partitions.

Examples

			Triangle begins:
  1;
  1, 1;
  1, 3, 2;
  1, 4, 2, 10,  5;
  1, 5, 5, 15, 15, 35, 14;
  ...
___________
The first few free cumulants in terms of the free moments are
  c_1 = m_1
  c_2 = m_2 - m_1^2
  c_3 = m_3 - 3 m_2 m_1 + 2 m_1^3
  c_4 = m_4 - 2 m_2^2 - 4m_3 m_1 + 10 m_2 m_1^2 - 5 m_1^4
  c_5 = m_5 - 5 m_2  m_3 - 5  m_4 m_1 + 15  m_2^2 m_1 + 15 m_3 m_1^2 - 35 m_2 m_1^3 + 14 m_1^5
___________
Conversely, from A134264, these free moments in terms of the free cumulants are
  m_1 = c_1
  m_2 = c_2 + c_1^2
  m_3 = c_3 + 3 c_2 c_1 + c_1^3
  m_4 = c_4 + + 2 c_2^2 + 4 c_3 c_1 + 6 c_2 c_1^2 + c_1^4
  m_5 = c_5 + 5 c_2 c_3 + 5 c_4 c_1 + 10 c_2^2 c_1 + 10 c_3 c_1^2  + 10 c_2 c_1^3 + c_1^5
___________
		

Crossrefs

Programs

  • PARI
    mv(n)={eval(Str("'m",n))}
    Trm(m,v)={my(S=Set(v)); for(i=1, #S, my(x=S[i]); m=polcoef(m, #select(y->y==x, v), mv(x))); m}
    Q(n)={polcoef(-x/serreverse(x*(1 + sum(k=1, n, -x^k*mv(k), O(x*x^n)))), n)}
    row(n)={my(q=Q(n)); [Trm(q,Vec(v)) | v<-partitions(n)]}
    { for(n=1, 7, print(row(n))) } \\ Andrew Howroyd, Feb 01 2022
    
  • PARI
    C(v)={my(n=vecsum(v), S=Set(v)); (n+#v-2)!/(n-1)!/prod(i=1, #S, my(x=S[i]); (#select(y->y==x, v))!)}
    row(n)=[C(Vec(p)) | p<-partitions(n)]
    { for(n=1, 7, print(row(n))) } \\ Andrew Howroyd, Feb 01 2022

Formula

O.g.f.: C(x) = 1 + c_1 x + c_2 x^2 + ... = x / (x + m_1 x^2 + m_2 x^3 + m_3 x^4 + ...)^(-1) = x / M^(-1)(x), the shifted reciprocal of the compositional inverse of a power series M(x) = x + m_1 x^2 + m_2 x^3 + ..., the o.g.f. of the free moments m_n in free probability theory.
Row sums: big Schroeder numbers A006318.
Refinement of A060693 and A088617, i.e., by letting m_n = -t and removing all resulting signs, the elements of these two lower triangular matrices are generated.
The coefficients of the highest order terms in m_1^n of the free moment partition polynomials are the signed Catalan numbers A000108.
Taking the derivative with respect to the indeterminate m_1 generates the Lagrange inversion partition polynomials, with shifted indices, of A133437 and A111785 with an overall scale factor. These Lagrange inversion polynomials are the refined Euler characteristic polynomials of the associahedra. E.g.,
D_{m_1} c_5 = 5 (- m_4 + 3 m_2^2 + 6 m_3 m_1 - 21 m_2 m_1^2 + 14 m_1^4). An analogous differential formula that applies to the classical formal cumulants in relation to the permutahedra is stated in my 2012 comment in A127671.
The o.g.f. satisfies the partial differential equations D_{m_1} (x / C(x)) = -(1/3) D_x (x / C(x))^3 and D_{m_1} (C(x) / x) = D_x (x / C(x)), where D_{m_1} and D_x are the partial derivatives with respect to m_1 and x.
More generally, D_{m_n} (x / C(x)) = -(1/(n+2)) D_x (x / C(x))^{n+2), equivalent to D_{m_n} M^(-1)(x) = -(1/(n+2)) D_x (M^(-1)(x))^{n+2). Equations of this type are found in Zhou (see eqn. 44 on p. 11), characterizing the KdV hierarchy. These differential equations can be transformed into the inviscid Burgers-Hopf partial differential equation (see, e.g., A133437, A086810, A001764, A002293, A133932, A134685, and A276850).
The formal free cumulants when identified as the indeterminates of the noncrossing Lagrange inversion partition polynomials NCP_n(c_1,c_2,...,c_n) = m_n of A134264 (as in the example section) satisfy the partial differential equations D_{m_k} NCP_n(c_1, ..., c_n) = d(m_n)/dm_k = delta_{n-k}, where delta_{n} is the Kronecker delta which is zero for all integers n other than n = 0, for which it evaluates as unity. This provides a recursion method for determining the partial derivatives dc_n/dm_k from the partial derivatives dc_p/dm_k and cumulants c_p with k <= p < n. For example, dc_k/dm_j = 0 for j > k and dc_k/dm_k = 1, so dm_3/dm_2 = 0 = D_{m_2} (c_3 + 3 c_2 c_1 + c_1^3) = dc_3/dm_2 + 3 c_1 dc_2/dm_2 = dc_3/dm_2 + 3 c_1 , implying dc_3/dm_2 = -3 c_1 = -3 m_1.
T(n,k) = (n+j-2)!/((n-1)!*Product_{i>=1} s_i!), where (1*s_1 + 2*s_2 + ... = n) is the k-th partition of n and j = s_1 + s_2 + ... is the number of parts. - Andrew Howroyd, Feb 01 2022
Conjecture: free cumulants in terms of the free moments are R(n,1) for n > 0 where R(n,k) = R(n-1,k+1) - Sum_{j=1..n-1} R(j,k)*R(n-j,1) for n > 1, k > 0 with R(1,k) = m_k for k > 0. - Mikhail Kurkov, Mar 30 2025

Extensions

Terms a(19) and beyond from Andrew Howroyd, Feb 01 2022
Previous Showing 11-17 of 17 results.