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.
%I A356144 #46 Feb 17 2024 12:32:29
%S A356144 1,-1,1,-1,-1,2,-1,1,-3,2,1,-1,-1,4,-4,-2,5,-1,-1,1,-5,8,2,-4,-2,-4,5,
%T A356144 4,-4,-1,-1,6,-12,-3,8,18,-6,-14,13,2,-16,14,0,-8,-1,1,-7,18,3,-20,0,
%U A356144 -15,8,18,57,6,-54,-15,-12,84,-30,-48,14,14,-8,-13,-1,-1,8,-24,-4,32,51,-27,-16,-6,171,-42,-177,50,90,-18,456,-276,-246,-15,30,154,-42,124,-166,-113,42,6,-21,-19,-1
%N A356144 Coefficients of the set of partition polynomials [RT] = [P][E]; i.e., coefficients of polynomials resulting from using the set of refined Eulerian polynomials, [E], of A145271 as the indeterminates of the set of permutahedra polynomials, [P], of A133314. Irregular triangle read by rows with lengths given by A000041.
%C A356144 I stipulate that the row lengths are A000041, but this imposes the insertion of a zero as a coefficient of a monomial for the polynomial RT_7 and for RT_8. The number of nonzero coefficients in each higher order polynomial remains to be determined. The monomials of the partition polynomials are arranged in the order (bottom to top) in Abramowitz and Stegun (starting on p. 831, link in A000041).
%C A356144 The analytic interpretation of these coefficients is related to the e.g.f.s of reciprocals of the derivatives (slopes of tangents) of a pair of compositionally inverse e.g.f.s as explicitly shown in the formulas.
%C A356144 With the notation introduced in the formula section, this set of partition polynomials, [RT], is the e.g.f. counterpart to the special Schur expansion coefficients [b], or [K], of A355201 for o.g.f.s. and is conjugate dual to the Lagrange inversion polynomials [L] of A134685.
%C A356144 For example, as shown in the formulas, [RT] = [P][E] = [P][L][P] where [P] is the set of polynomials of A133314, the refined Euler characteristic polynomials of the permutahedra; [E], the set A145271, the refined Eulerian polynomials; and [L], the set A134685, the classic Lagrange inversion polynomials--all related to transformations of e.g.f.s, or Taylor series, for which [RT], [L], and [E] can each be used to give the compositional inverse and [P], the multiplicative inverse, or reciprocal.
%C A356144 On the other hand, as shown in formulas for A355201, [K] = [R][N] = [R][A][R] where [R] is the set A263633 (mod signs), refined Pascal polynomials; [N], the set A134264, the refined Narayana, or noncrossing partition, polynomials; and [A], the set A133437, the refined Euler characteristic polynomials of the associahedra--all related to transformations of o.g.f.s, or power series, for which [K], [A], and [N] can each be used to give the compositional inverse and [R], the multiplicative inverse, or reciprocal. This is related to three pairs of compositionally inverse series--two pairs of Laurent series and one pair of power series.
%H A356144 Tom Copeland, <a href="https://tcjpn.wordpress.com/2022/07/27/one-matrix-to-rule-them-all/#more-9405">One Matrix to Rule Them All</a>, 2022.
%F A356144 Denote this set of partition polynomials by [RT], the permutahedra polynomials of A133314 by [P], the refined Eulerian polynomials of A145271 by [E], and the Lagrange inversion polynomials of A134685 for e.g.f.s by [L]. Let the typically noncommutative product of two sets, e.g., [P][E], represent the substitution of the polynomials of [E] for the indeterminates of [P], i.e., a composition at the level of the indeterminates (see A356145 for examples). Let [I] be the substitutional identity transformation, and mark the substitutional inverse with the superscript -1. Then the following relations hold.
%F A356144 [RT] = [P][E] = [P][L][P] = [P]^{-1}[L][P] = [P][L][P]^{-1} since [P] is an involution, i.e., [P]^2 = [I], or [P] = [P]^{-1}, so [RT] and [L] are conjugate duals.
%F A356144 [RT]^{-1} = ([P][E])^{-1} = [E]^{-1}[P] = ([P][L][P])^{-1} = [P][L][P] = [RT], with [E]^{-1} = A356145, since [L] and [P] are involutions, so is [RT], i.e., [RT]^2 = [I].
%F A356144 RT_n(a_1,a_2,...,a_n) = D_{x=0}^n 1 / [ D_x f^{(-1)}(x)] for which D_x is the derivative w.r.t. x and the indeterminates are defined by 1 / [D_x f(x)] = 1 + a_1 x + a_2 x^2/2! + a_3 x^3/3! + ... with f(x) and f^{(-1)}(x) a compositional inverse pair of formal Taylor series, or e.g.f.s. This is the analytic equivalent of the algebraic relation [RT] = [P][E]. In words, the partition polynomials of row n (initial row is 0) is the n-th coefficient of the formal Taylor series of the reciprocal of the derivative of the compositional inverse of a function in terms of the Taylor series coefficients of the reciprocal of the derivative of that function. Note the correspondence with the analytic interpretation of [E]^{-1} of A356145, consistent with the algebraic identities above.
%F A356144 RT_n(a_1,a_2,...,a_n) = D_{x=0}^n f'(f^{(-1)}(x)) also, by the inverse function theorem, where the prime denotes differentiation with respect to the argument of the function.
%F A356144 With all a_k = (-1)^k, RT_0 = RT_1 = 1, otherwise RT_n = 0. This is determined with f(x) = e^{x}-1 and f^{(-1)}(x) = log(1+x).
%F A356144 With all a_k = 1, RT_0 = 1, RT_1 = -1, otherwise RT_n=0. This is determined with f(x) =1-e^{-x} and f^{(-1)}(x) = -log(1-x).
%F A356144 With all a_k = -1, RT_0 = 1 and RT_n = 2^(n-1) otherwise. This is determined with f(x) = (x - log(2-e^x))/2 and f^{(-1)}(x) = x - log(cosh(x)). (Careful, these are not the row sums of the absolute values of the numerical coefficients, which for the first ten polynomials are 1, 1, 2, 4, 8, 18, 40, 122, 446, and 2428.)
%F A356144 With a_k = k! 2^k, RT_0 = 1 and RT_n = -2*(2(n-1))! / (n-1)! = -2*n!*A000108(n-1) otherwise. This is determined with f(x) = x - x^2 and f^{(-1)}(x) = (1 - sqrt(1-4x))/2. Similar relations hold for the Fuss-Catalan sequences with f(x) = x - x^{m+1} for m > 1.
%e A356144 Arranged by rows, the coefficients are
%e A356144 0) 1;
%e A356144 1) -1;
%e A356144 2) 1, -1;
%e A356144 3) -1, 2, -1;
%e A356144 4) 1, -3, 2, 1, -1;
%e A356144 5) -1, 4, -4, -2, 5, -1, -1;
%e A356144 6) 1, -5, 8, 2, -4, -2, -4, 5, 4, -4, -1;
%e A356144 7) -1, 6, -12, -3, 8, 18, -6, -14, 13, 2, -16, 14, 0, -8, -1;
%e A356144 8) 1, -7, 18, 3, -20, 0, -15, 8, 18, 57, 6, -54, -15, -12, 84, -30, -48, 14, 14, -8, -13, -1;
%e A356144 . . .
%e A356144 The first few partition polynomials are
%e A356144 RT_0 = 1,
%e A356144 RT_1 = -a1,
%e A356144 RT_2 = a1^2 - a2,
%e A356144 RT_3 = -a1^3 + 2 a1 a2 - a3,
%e A356144 Rt_4 = a1^4 - 3 a1^2 a2 + 2 a2^2 + a1 a3 - a4,
%e A356144 RT_5 = -a1^5 + 4 a1^3 a2 - 4 a1 a2^2 - 2 a1^2 a3 + 5 a2 a3 - a1 a4 - a5,
%e A356144 RT_6 = a1^6 - 5 a1^4 a2 + 8 a1^2 a2^2 + 2 a1^3 a3 - 4 a2^3 - 2 a1 a2 a3 - 4 a1^2 a4 + 5 a3^2 + 4 a2 a4 - 4 a1 a5 - a6,
%e A356144 RT_7 = -a1^7 + 6 a1^5 a2 - 12*a1^3 a2^2 - 3 a1^4 a3 + 8 a1 a2^3 + 18 a1^2 a2 a3 - 6 a1^3 a4 - 14 a2^2 a3 + 13 a1 a3^2 + 2 a1 a2 a4 - 16 a1^2 a5 + 14 a3 a4 + 0 a2 a5 - 8 a1 a6 - a7,
%e A356144 RT_8 = a1^8 - 7 a1^6 a2 + 18 a1^4 a2^2 + 3 a1^5 a3 - 20 a1^2 a2^3 + 0 a1^3 a2 a3 - 15 a1^4 a4 + 8 a2^4 + 18 a1 a2^2 a3 + 57 a1^2 a3^2 + 6 a1^2 a2 a4 - 54 a1^3 a5 - 15 a2 a3^2 - 12 a2^2 a4 + 84 a1 a3 a4 - 30 a1 a2 a5 - 48 a1^2 a6 + 14 a4^2 + 14 a3 a5 - 8 a2 a6 - 13 a1 a7 - a8.
%t A356144 rows[nn_] := {{1}}~Join~With[{s = 1 / D[InverseSeries[Integrate[1/(1 + Sum[c[k] x^k/k!, {k, nn}] + O[x]^(nn+1)), x]], x]}, Table[Coefficient[n! s, x^n Product[c[t], {t, p}]], {n, nn}, {p, Reverse[Sort[Sort /@ IntegerPartitions[n]]]}]];
%t A356144 rows[7] // Flatten (* _Andrey Zabolotskiy_, Feb 17 2024 *)
%o A356144 (SageMath)
%o A356144 B.<a1,a2, a3, a4, a5, a6, a7, a8, a9, a10> = PolynomialRing(ZZ)
%o A356144 A.<x> = PowerSeriesRing(B)
%o A356144 f = 1/(1 + a1*x + a2*x^2/factorial(2) + a3*x^3/factorial(3) + a4*x^4/factorial(4) + a5*x^5/factorial(5) + a6*x^6/factorial(6) + a7*x^7/factorial(7) + a8*x^8/factorial(8) + a9*x^9/factorial(9) + a10*x^10/factorial(10) )
%o A356144 g = integrate(f)
%o A356144 h = g.reverse()
%o A356144 w = derivative(h,x)
%o A356144 I = 1 / w
%o A356144 # Added by # _Peter Luschny_, Feb 17 2024:
%o A356144 # The list of coefficients in sparse format (i.e. without the zeros):
%o A356144 for n, c in enumerate(I.list()[:10]):
%o A356144 print(f"RT[{n}]", (factorial(n)*c).coefficients())
%Y A356144 Cf. A000108, A133314, A133437, A134685, A145271, A263633, A355201, A356145.
%K A356144 sign,tabf
%O A356144 0,6
%A A356144 _Tom Copeland_, Jul 27 2022
%E A356144 Order of terms in rows 4-6 corrected by _Andrey Zabolotskiy_, Feb 17 2024