Exploring some relations among the multinomial coefficients of OEIS A036038 and the compositional inversion formulas of A134264, A248120, and A248927, related to numerous combinatorial structures and areas of analysis, including noncrossing partitions and free probability, I came across the Jack symmetric functions in an infinite number of variables as presented in “Some combinatorial properties of the Jack symmetric functions” by Stanley.

Equation a) on page 80 of his paper leads to an umbral generating formula for the related Jack symmetric polynomials (JSP)

,

where is to be regarded as a regular variable until the expression is reduced to monomials at which time it is to be evaluated as with , essentially the row polynomials of A094638 comprised of the Stirling numbers of the first kind. For example,

,

where the polynomial has been expressed in the symmetric monomial polynomials (SMP), easily extended to an indefinite number of variables. The factors multplying the SMPs are the multinomial coefficients of A036038 and remain independent of the number of variables. Each summand of an SMP has the same configuration of exponents and subscripts, allowing the products of to be easily determined and factored out after the umbral evaluation.

“MOPS: Multivariate orthogonal polynomials (symbolically)” by Dumitriu, Edelman, and Shuman contain examples for the third and fourth JSPs, but the third has the coefficients erroneously transposed for the factorial.

Using the operator identities in A094638, a Rodriques-like generator for the JSPs can be devised.

,

so

,

and, with ,

,

with treated as independent of with respect to the derivations and, after being passed unscathed to the left of all derivations, is finally evaluated as .

Then, with

.

This may be generalized just as the Laguerre polynomials are to the associated Laguerre polynomials by conjugating with rather than .

A generating function for the full set of JSPs is

.

This can be evaluated by noting, with

and ,

that

.

Then, for and , the operation reduces to

.

Then the generator gives, since ultimately ,

,

or

.

You can use the generalized Leibnitz formula to relate this back to the multinomial coefficients:

,

where acts as only on .

The e.g.f. is naturally consistent with the e.g.f. for A094638, which with a change of notation is

,

and is an e.g.f. for plane m-ary trees with , so

.

Therefore, the discussions in A134264 and the Hirzebruch criterion post below on repeated exponentiaton of an e.g.f., binomial convolutions, and umbral substitution apply to these calculations when , giving connections among the OEIS entries cited at the top here and the multinomial coefficients.

Taking the log of the e.g.f. gives a relation between the symetric power sum polynomials / functions of the variables / indeterminates and the cumulants formed from the JSPs through A127671, or A263634.