cp's OEIS Frontend

"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

Also numerator of rational part of Haar measure on Grassmannian space G(n,1).

Also rational part of numerator of Gamma(n/2+1)/Gamma(n/2+1/2) (cf. A036039).

Let x(m) = x(m-2) + 1/x(m-1) for m >= 3, with x(1)=x(2)=1. Then the numerator of x(n+2) equals the denominator of n!!/(n+1)!! for n >= 0, where the double factorials are given by A006882. - Joseph E. Cooper III (easonrevant(AT)gmail.com), Nov 07 2010, as corrected in Cooper (2015).

Numerator of (n-1)/( (n-2)/( .../1)), with an empty fraction taken to be 1. - Flávio V. Fernandes, Jan 31 2025

Also the Bell transform of A038048(n+1) and the inverse Bell transform of A180563(n+1) (adding 1,0,0,.. as column 0). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016

Named after the Italian mathematician Francesco Flores D'Arcais (1849-1927). - Amiram Eldar, Jun 13 2021

Limit_{k->oo} k*(1-Gamma(1+1/k)) = -Gamma'(1) = gamma = 0.577....

Decimal expansion of the higher-order exponential integral constant gamma(2,1). The higher-order exponential integrals, see A163931, are defined by E(x,m,n) = x^(n-1)*Integral_{t=x..oo} (E(t,m-1,n)/t^n) dt for m >= 1 and n >= 1, with E(x,m=0,n) = exp(-x). The series expansions of the higher-order exponential integrals are dominated by the gamma(k,n) and the alpha(k,n) constants, see A163927. - Johannes W. Meijer and Nico Baken, Aug 13 2009

Let f(t) = -log(1 - u(.)*t) = Sum_{n>=1} (u_n / n) * t^n.

If u_1 is not equal to 0, then the compositional inverse for f(t) is given by g(t) = Sum_{j>=1} P(n,t) where, with u_n denoted by (n'),

P(1,t) = (1')^(-1) * [ 1 ] * t

P(2,t) = (1')^(-3) * [ -1 (2') ] * t^2 / 2!

P(3,t) = (1')^(-5) * [ 3 (2')^2 - 2 (1')(3') ] * t^3 / 3!

P(4,t) = (1')^(-7) * [ -15 (2')^3 + 20 (1')(2')(3') - 6 (1')^2 (4') ] * t^4 / 4!

P(5,t) = (1')^(-9) * [ 105 (2')^4 - 210 (1') (2')^2 (3') + 40 (1')^2 (3')^2 + 90 (1')^2 (2') (4') - 24 (1')^3 (5') ] * t^5 / 5!

P(6,t) = (1')^(-11) * [ -945 (2')^5 + 2520 (1') (2')^3 (3') - 1120 (1')^2 (2') (3')^2 - 1260 (1')^2 (2')^2 (4') + 420 (1')^3 (3')(4') + 504 (1')^3 (2')(5') - 120 (1')^4 (6') ] * t^6 / 6!

See A134685 for more information.

From Tom Copeland, Sep 28 2016: (Start)

P(7,t) = (1')^(-13) * [ 10395 (2')^6 - 34650 (1')(2')^4(3') + (1')^2 [25200 (2')^2(3')^2 + 18900 (2')^3(4')] - (1')^3 [2240 (3')^3 + 15120 (2')(3')(4') + 9072 (2')^2(5')] + (1')^4 [1260 (4')^2 + 2688 (3')(5') + 3360 (2')(6')] - 720 (1')^5(7')] * t^7 / 7!

P(8,t) = (1')^(-15) * [ -135135 (2')^7 + 540540 (1')(2')^5(3') - (1')^2 [554400 (2')^3(3')^2 + 311850 (2')^4(4')] + (1')^3 [123200 (2')(3')^3 + 415800 (2')^2(3')(4') + 166320 (2')^3(5')] - (1')^4 [50400 (3')^2(4') + 56700 (2')(4')^2 + 120960 (2')(3')(5') + 75600 (2')^2(6')] + (1')^5 [18144 (4')(5') + 20160 (3')(6') + 25920 (2')(7')] - 5040 (1')^6(8')] * t^8 / 8! (End)

|T(n,k)| gives the sum of the M_2 multinomial numbers (A036039) for those partitions of n with exactly k odd parts. E.g.: |T(6,2)| = 144 + 40 = 184 from the partitions of 6 with exactly two odd parts, namely (1,5) and (3,3), with M_2 numbers 144 and 40. Proof via the general Jabotinsky triangle formula for |T(n,k)| using partitions of n into k parts and their M_3 numbers (A036040). Then with the special e.g.f. of the (unsigned) k=1 column, f(x):= arctanh(x), only odd parts survive and the M_3 numbers are changed into the M_2 numbers. For the Knuth reference on Jabotinsky triangles see A039692. - Wolfdieter Lang, Feb 24 2005 [The first two sentences have been corrected thanks to the comment by José H. Nieto S. given below. - Wolfdieter Lang, Jan 16 2012]

|T(n,k)| gives the number of permutations of {1,2,...,n} (degree n permutations) with the number of odd cycles equal to k. E.g.: |T(5,3)|= 20 from the 20 degree 5 permutations with cycle structure (.)(.)(...). Proof: Use the cycle index polynomial for the symmetric group S_n (see the M_2 array A036039 or A102189) together with the partition interpretation of |T(n,k)| given above. - Wolfdieter Lang, Feb 24 2005 [See the following José H. Nieto S. correction. - Wolfdieter Lang, Jan 16 2012]

The first sentence of the above comment is inexact, it should be "|T(n,k)| gives the number of degree n permutations which decompose into exactly k odd cycles". The number of degree n permutations with k odd cycles (and, possibly, other cycles of even length) is given by A060524. - José H. Nieto S., Jan 15 2012

The unsigned triangle with e.g.f. exp(x*arctanh(z)) is the associated Jabotinsky type triangle for the Sheffer type triangle A060524. See the comments there. - Wolfdieter Lang, Feb 24 2005

Also the Bell transform of the sequence (-1)^(n/2)*A005359(n) without column 0. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 28 2016

From Johannes W. Meijer, Jun 27 2016: (Start)

This constant is related to the values of zeta(2*n-1) of the Riemann zeta function and the Euler Mascheroni constant gamma. If we define Z(n) = (1/n) * (sum(zeta(2*n-2*k-1) * Z(k), k=0..n-2) + gamma * Z(n-1)), with Z(0) = 1, then limit(Z(n), n -> infinity) = 2/exp(1).

Similar formulas appear in A090998 and A112302.

The structure of the n! * Z(n) formulas leads to the multinomial coefficients A036039. (End).

The unsigned numbers give A048996. They are not listed on pp. 831-832 of Abramowitz and Stegun (reference given in A103921). One could call these numbers M_0 (like M_1, M_2, M_3 given in A036038, A036039, A036040, resp.).

The sequence of row lengths is A000041(n) (partition numbers).

The sign is (-1)^(n + m(n,k)) with m(n,k) the number of parts of the k-th partition of n taken in the mentioned order. For m(n,k), see A036043.

The row sum is 1 for n = 1, and 0 otherwise. The unsigned row sum is 2^(n-1) = A000079(n-1) for n >= 1.

The complete symmetric polynomial is also h(n; a[1],...,a[n]) = Det(A_n) with the matrix elements of the n X n matrix A_n given by A_n(k, k+1) = 1 for 1 <= k < n, A(k, m) = a[k-m+1] for n >= k >= m >= 1, and 0 otherwise. [For an explanation of this statement, see the example for n = 4 below. See also p. 3 in MacMahon (1960).]

T(n,k) tells how often k appears among the first n! entries of A198380, i.e., how many permutations of n elements have the cycle type denoted by k.

This triangle is a refinement of the rencontres numbers A008290, which tell only how many permutations of n elements actually move a certain number of elements. How many of these permutations have a certain cycle type is a more detailed question, answered by this triangle.

The rows are counted from 1, the columns from 0.

Row lengths: 1, 2, 3, 5, 7, 11, ... (partition numbers A000041).

Row sums: 1, 2, 6, 24, 120, 720, ... (factorial numbers A000142).

Row maxima: 1, 1, 3, 8, 30, 144, ... (A059171).

Distinct entries per row: 1, 1, 3, 4, 6, 7, ... (A073906).

It follows from the formula given by Carlos Mafra that the rows of the triangle correspond to the coefficients of the modified Bell polynomials. - Sela Fried, Dec 08 2021

For k>0, the k-th column of triangle T(n,k) is a scaled copy of binomial coefficients binomial(n,q) where q is the least value for which p(q) exceeds or equals k+1, with p() being the integer partitions counting function, A000041(q). E.g., for column 4, the relevant binomial coefficients have q=4 as p(4)=5; for column 5, we have q=5 as p(5)>6; for column 6, we have q=5 as p(5)=7. The scale factor for column k is given by A385081(k+1). This triangle gives coefficients for expressing the characteristic polynomial and determinant of a matrix solely in terms of traces; see extended comment, below, under "Links". - Gregory Gerard Wojnar, Jun 24 2025

The multinomial transform [MNL] transforms an input sequence b(n) into the output sequence a(n). Given the fact that the structure of the a(n) formulas, see the examples, lead to the multinomial coefficients A036039 the MNL transform seems to be an appropriate name for this transform. The multinomial transform is related to the exponential transform, see A274804 and the third formula. For the inverse multinomial transform [IML] see A274844.

To preserve the identity IML[MNL[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, otherwise information about b(0) will be lost in transformation.

In the a(n) formulas, see the examples, the multinomial coefficients A036039 appear.

We observe that a(0) = 1 and that this term provides no information about any value of b(n), this notwithstanding we will start the a(n) sequence with a(0) = 1.

The Maple programs can be used to generate the multinomial transform of a sequence. The first program uses the first formula which was found by Paul D. Hanna, see A158876, and Vladimir Kruchinin, see A215915. The second program uses properties of the e.g.f., see the sequences A158876, A213507, A244430 and A274539 and the third formula. The third program uses information about the inverse multinomial transform, see A274844.

Some MNL transform pairs are, n >= 1: A000045(n) and A244430(n-1); A000045(n+1) and A213527(n-1); A000108(n) and A213507(n-1); A000108(n-1) and A243953(n-1); A000142(n) and A158876(n-1); A000203(n) and A053529(n-1); A000110(n) and A274539(n-1); A000041(n) and A215915(n-1); A000035(n-1) and A177145(n-1); A179184(n) and A038205(n-1); A267936(n) and A000266(n-1); A267871(n) and A000090(n-1); A193356(n) and A088009(n-1).