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.

Showing 1-3 of 3 results.

A002293 Number of dissections of a polygon: binomial(4*n, n)/(3*n + 1).

Original entry on oeis.org

1, 1, 4, 22, 140, 969, 7084, 53820, 420732, 3362260, 27343888, 225568798, 1882933364, 15875338990, 134993766600, 1156393243320, 9969937491420, 86445222719724, 753310723010608, 6594154339031800, 57956002331347120, 511238042454541545
Offset: 0

Views

Author

Keywords

Comments

The number of rooted loopless n-edge maps in the plane (planar with a distinguished outside face). - Valery A. Liskovets, Mar 17 2005
Number of lattice paths from (1,0) to (3*n+1,n) which, starting from (1,0), only utilize the steps +(1,0) and +(0,1) and additionally, the paths lie completely below the line y = (1/3)*x (i.e., if (a,b) is in the path, then b < a/3). - Joseph Cooper (jecooper(AT)mit.edu), Feb 07 2006
Number of length-n restricted growth strings (RGS) [s(0), s(1), ..., s(n-1)] where s(0) = 0 and s(k) <= s(k-1) + 3, see fxtbook link below. - Joerg Arndt, Apr 08 2011
From Wolfdieter Lang, Sep 14 2007: (Start)
a(n), n >= 1, enumerates quartic trees (rooted, ordered, incomplete) with n vertices (including the root).
Pfaff-Fuss-Catalan sequence C^{m}_n for m = 4. See the Graham et al. reference, p. 347. eq. 7.66. (Second edition, p. 361, eq. 7.67.) See also the Pólya-Szegő reference.
Also 4-Raney sequence. See the Graham et al. reference, pp. 346-347.
(End)
Bacher: "We describe the statistics of checkerboard triangulations obtained by coloring black every other triangle in triangulations of convex polygons." The current sequence (A002293) occurs on p. 12 as one of two "extremal sequences" of an array of coefficients of polynomials, whose generating functions are given in terms of hypergeometric functions. - Jonathan Vos Post, Oct 05 2007
A generating function in terms of a (labyrinthine) solution to a depressed quartic equation is given in the Copeland link for signed A005810. With D(z,t) that g.f., a g.f. for signed A002293 is {[-1+1/D(z,t)]/(4t)}^(1/3). - Tom Copeland, Oct 10 2012
For a relation to the inviscid Burgers's equation, see A001764. - Tom Copeland, Feb 15 2014
For relations to compositional inversion, the Legendre transform, and convex geometry, see the Copeland, the Schuetz and Whieldon, and the Gross (p. 58) links. - Tom Copeland, Feb 21 2017 (See also Gross et al. in A062994. - Tom Copeland, Dec 24 2019)
This is the number of A'Campo bicolored forests of degree n and co-dimension 0. This can be shown using generating functions or a combinatorial approach. See Combe and Jugé link below. - Noemie Combe, Feb 28 2017
Conjecturally, a(n) is the number of 3-uniform words over the alphabet [n] that avoid the patterns 231 and 221 (see the Defant and Kravitz link). - Colin Defant, Sep 26 2018
The compositional inverse o.g.f. pair in Copeland's comment above are related to a pair of quantum fields in Balduf's thesis by Theorem 4.2 on p. 92. Cf. A001764. - Tom Copeland, Dec 13 2019
a(n) is the total number of down steps before the first up step in all 3_1-Dyck paths of length 4*n. A 3_1-Dyck path is a lattice path with steps (1, 3), (1, -1) that starts and ends at y = 0 and stays above the line y = -1. - Sarah Selkirk, May 10 2020
a(n) is the number of pairs (A<=B) of noncrossing partitions of [2n] such that every block of A has exactly two elements. In fact, it is proved that a(n) is the number of planar tied arc diagrams with n arcs (see Aicardi link below). A planar diagram with n arcs represents a noncrossing partition A of [2n] with n blocks, each block containing the endpoints of one arc; each tie connects two arcs, so that the ties define a partition B >= A: the endpoints of two arcs connected by a tie belong to the same block of B. Ties do not cross arcs nor other ties iff B has a planar diagram, i.e., B is a noncrossing partition. - Francesca Aicardi, Nov 07 2022
Dropping the initial 1 (starting 1, 4, 22 with offset 1) yields the REVERT transformation 1, -4 ,10, -20, 35.. essentially A000292 without leading 0. - R. J. Mathar, Aug 17 2023
Number of rooted polyominoes composed of n pentagonal cells of the hyperbolic regular tiling with Schläfli symbol {5,oo}. A rooted polyomino has one external edge identified, and chiral pairs are counted as two. A stereographic projection of the {5,oo} tiling on the Poincaré disk can be obtained via the Christensson link. - Robert A. Russell, Jan 27 2024
This is instance k = 4 of the generalized Catalan family {C(k, n)}A130564.%20-%20_Wolfdieter%20Lang">{n>=0} given in a comment of A130564. - _Wolfdieter Lang, Feb 05 2024
a(n) is the cardinality of the planar ramified Jones monoid PR(J_n). - Diego Arcis, Nov 21 2024

Examples

			There are a(2) = 4 quartic trees (vertex degree <= 4 and 4 possible branchings) with 2 vertices (one of them the root). Adding one more branch (one more vertex) to these four trees yields 4*4 + 6 = 22 = a(3) such trees.
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 23.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, pp. 200, 347.
  • Peter Hilton and Jean Pedersen, Catalan numbers, their generalization, and their uses, Math. Intelligencer 13 (1991), no. 2, 64-75.
  • V. A. Liskovets and T. R. Walsh, Enumeration of unrooted maps on the plane, Rapport technique, UQAM, No. 2005-01, Montreal, Canada, 2005.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis, Springer-Verlag, Heidelberg, New York, 2 vols., 1972, Vol. 1, problem 211, p. 146 with solution on p. 348.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Column k=3 of triangle A062993 and A070914.
Cf. A000260, A002295, A002296, A027836, A062994, A346646 (binomial transform), A346664 (inverse binomial transform).
Polyominoes: A005038 (oriented), A005040 (unoriented), A369471 (chiral), A369472 (achiral), A001764 {4,oo}, A002294 {6,oo}.
Cf. A130564 (for generalized Catalan C(k, n), for = 4).

Programs

  • GAP
    List([0..22],n->Binomial(4*n,n)/(3*n+1)); # Muniru A Asiru, Nov 01 2018
  • Magma
    [ Binomial(4*n,n)/(3*n+1): n in [0..50]]; // Vincenzo Librandi, Apr 19 2011
    
  • Maple
    series(RootOf(g = 1+x*g^4, g),x=0,20); # Mark van Hoeij, Nov 10 2011
    seq(binomial(4*n, n)/(3*n+1),n=0..20); # Robert FERREOL, Apr 02 2015
    # Using the integral representation above:
    Digits:=6;
    R:=proc(x)((I + sqrt(3))*(4*sqrt(256 - 27*x) - 12*I*sqrt(3)*sqrt(x))^(1/3))/16 - ((I - sqrt(3))*(4*sqrt(256 - 27*x) + 12*I*sqrt(3)*sqrt(x))^(1/3))/16;end;
    W:=proc(x) x^(-3/4) * sqrt(4*R(x) - 3^(3/4)*x^(1/4)/sqrt(R(x)))/(2*3^(1/4)*Pi);end;
    # Attention: W(x) is singular at x = 0. Integration is done from  a very small positive x to x = 256/27.
    # For a(8):  ... gives 420732
    evalf(int(x^8*W(x),x=0.000001..256/27));
    # Karol A. Penson, Jul 05 2024
  • Mathematica
    CoefficientList[InverseSeries[ Series[ y - y^4, {y, 0, 60}], x], x][[Range[2, 60, 3]]]
    Table[Binomial[4n,n]/(3n+1),{n,0,25}] (* Harvey P. Dale, Apr 18 2011 *)
    CoefficientList[1 + InverseSeries[Series[x/(1 + x)^4, {x, 0, 60}]], x] (* Gheorghe Coserea, Aug 12 2015 *)
    terms = 22; A[] = 0; Do[A[x] = 1 + x*A[x]^4 + O[x]^terms, terms];
    CoefficientList[A[x], x] (* Jean-François Alcover, Jan 13 2018 *)
  • PARI
    a(n)=binomial(4*n,n)/(3*n+1) /* Charles R Greathouse IV, Jun 16 2011 */
    
  • PARI
    my(x='x+O('x^33)); Vec(1 + serreverse(x/(1+x)^4)) \\ Gheorghe Coserea, Aug 12 2015
    
  • Python
    A002293_list, x = [1], 1
    for n in range(100):
        x = x*4*(4*n+3)*(4*n+2)*(4*n+1)//((3*n+2)*(3*n+3)*(3*n+4))
        A002293_list.append(x) # Chai Wah Wu, Feb 19 2016
    

Formula

O.g.f. satisfies: A(x) = 1 + x*A(x)^4 = 1/(1 - x*A(x)^3).
a(n) = binomial(4*n,n-1)/n, n >= 1, a(0) = 1. From the Lagrange series of the o.g.f. A(x) with its above given implicit equation.
From Karol A. Penson, Apr 02 2010: (Start)
Integral representation as n-th Hausdorff power moment of a positive function on the interval [0, 256/27]:
a(n) = Integral_{x=0..256/27}(x^n((3/256) * sqrt(2) * sqrt(3) * ((2/27) * 3^(3/4) * 27^(1/4) * 256^(/4) * hypergeom([-1/12, 1/4, 7/12], [1/2, 3/4], (27/256)*x)/(sqrt(Pi) * x^(3/4)) - (2/27) * sqrt(2) * sqrt(27) * sqrt(256) * hypergeom([1/6, 1/2, 5/6], [3/4, 5/4], (27/256)*x)/ (sqrt(Pi) * sqrt(x)) - (1/81) * 3^(1/4) * 27^(3/4) * 256^(1/4) * hypergeom([5/12, 3/4, 13/12], [5/4, 3/2], (27/256)*x/(sqrt(Pi)*x^(1/4)))/sqrt(Pi))).
This representation is unique as it represents the solution of the Hausdorff moment problem.
O.g.f.: hypergeom([1/4, 1/2, 3/4], [2/3, 4/3], (256/27)*x);
E.g.f.: hypergeom([1/4, 1/2, 3/4], [2/3, 1, 4/3], (256/27)*x). (End)
a(n) = upper left term in M^n, M = the production matrix:
1, 1
3, 3, 1
6, 6, 3, 1
...
(where 1, 3, 6, 10, ...) is the triangular series. - Gary W. Adamson, Jul 08 2011
O.g.f. satisfies g = 1+x*g^4. If h is the series reversion of x*g, so h(x*g)=x, then (x-h(x))/x^2 is the o.g.f. of A006013. - Mark van Hoeij, Nov 10 2011
a(n) = binomial(4*n+1, n)/(4*n+1) = A062993(n+2,2). - Robert FERREOL, Apr 02 2015
a(n) = Sum_{i=0..n-1} Sum_{j=0..n-1-i} Sum_{k=0..n-1-i-j} a(i)*a(j)*a(k)*a(n-1-i-j-k) for n>=1; and a(0) = 1. - Robert FERREOL, Apr 02 2015
a(n) ~ 2^(8*n+1/2) / (sqrt(Pi) * n^(3/2) * 3^(3*n+3/2)). - Vaclav Kotesovec, Jun 03 2015
From Peter Bala, Oct 16 2015: (Start)
A(x)^2 is o.g.f. for A069271; A(x)^3 is o.g.f. for A006632;
A(x)^5 is o.g.f. for A196678; A(x)^6 is o.g.f. for A006633;
A(x)^7 is o.g.f. for A233658; A(x)^8 is o.g.f. for A233666;
A(x)^9 is o.g.f. for A006634; A(x)^10 is o.g.f. for A233667. (End)
D-finite with recurrence: a(n+1) = a(n)*4*(4*n + 3)*(4*n + 2)*(4*n + 1)/((3*n + 2)*(3*n + 3)*(3*n + 4)). - Chai Wah Wu, Feb 19 2016
E.g.f.: F([1/4, 1/2, 3/4], [2/3, 1, 4/3], 256*x/27), where F is the generalized hypergeometric function. - Stefano Spezia, Dec 27 2019
x*A'(x)/A(x) = (A(x) - 1)/(- 3*A(x) + 4) = x + 7*x^2 + 55*x^3 + 455*x^4 + ... is the o.g.f. of A224274. Cf. A001764 and A002294 - A002296. - Peter Bala, Feb 04 2022
a(n) = hypergeom([1 - n, -3*n], [2], 1). Row sums of A173020. - Peter Bala, Aug 31 2023
G.f.: t*exp(4*t*hypergeom([1, 1, 5/4, 3/2, 7/4], [4/3, 5/3, 2, 2], (256*t)/27))+1. - Karol A. Penson, Dec 20 2023
From Karol A. Penson, Jul 03 2024: (Start)
a(n) = Integral_{x=0..256/27} x^(n)*W(x)dx, n>=0, where W(x) = x^(-3/4) * sqrt(4*R(x) - 3^(3/4)*x^(1/4)/sqrt(R(x)))/(2*3^(1/4)*Pi), with R(x) = ((i + sqrt(3))*(4*sqrt(256 - 27*x) -12*i*sqrt(3*x))^(1/3))/16 - ((i - sqrt(3))*(4*sqrt(256 - 27*x) + 12*i* sqrt(3*x))^(1/3))/16, where i is the imaginary unit.
The elementary function W(x) is positive on the interval x = (0, 256/27) and is equal to the combination of hypergeometric functions in my formula from 2010; see above.
(Pi*W(x))^6 satisfies an algebraic equation of order 6, with integer polynomials as coefficients. (End)
G.f.: (Sum_{n >= 0} binomial(4*n+1, n)*x^n) / (Sum_{n >= 0} binomial(4*n, n)*x^n). - Peter Bala, Dec 14 2024
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^7). - Seiichi Manyama, Jun 16 2025

A108767 Triangle read by rows: T(n,k) is number of paths from (0,0) to (3n,0) that stay in the first quadrant (but may touch the horizontal axis), consisting of steps u=(1,1), d=(1,-2) and have k peaks (i.e., ud's).

Original entry on oeis.org

1, 1, 2, 1, 6, 5, 1, 12, 28, 14, 1, 20, 90, 120, 42, 1, 30, 220, 550, 495, 132, 1, 42, 455, 1820, 3003, 2002, 429, 1, 56, 840, 4900, 12740, 15288, 8008, 1430, 1, 72, 1428, 11424, 42840, 79968, 74256, 31824, 4862, 1, 90, 2280, 23940, 122094, 325584, 465120
Offset: 1

Views

Author

Emeric Deutsch, Jun 24 2005

Keywords

Comments

Row sums yield A001764.
From Peter Bala, Sep 16 2012: (Start)
The number of 2-Dyck paths of order n with k peaks (Cigler). A 2-Dyck path of order n is a lattice path from (0,0) to (2*n,n) with steps (0,1) (North) and (1,0) (East) that never goes below the diagonal {2*i,i} 0 <= i <= n. A peak is a consecutive North East pair.
Row reverse of A120986. Described in A173020 as generalized Runyon numbers R_{n,k}^(2).
(End)
From Alexander Burstein, Jun 15 2020: (Start)
T(n,k) is the number of paths from (0,0) to (3n,0) that stay on or above the horizontal axis, with unit steps u=(1,1) and d=(1,-2), that have n+1-k peaks at even height.
T(n,k) is also the number of paths from (0,0) to (3n,0) that stay on or above the horizontal axis, with unit steps u=(1,1) and d=(1,-2), that have n-k peaks at odd height. (End)
An apparent refinement is A338135. - Tom Copeland, Oct 12 2020

Examples

			T(3,2)=6 because we have uuduuuudd, uuuduuudd, uuuuduudd, uuuudduud, uuuuududd and uuuuuddud.
Triangle starts:
  1;
  1,  2;
  1,  6,   5;
  1, 12,  28,   14;
  1, 20,  90,  120,   42;
  1, 30, 220,  550,  495,  132;
  1, 42, 455, 1820, 3003, 2002, 429;
  ...
		

Crossrefs

Runyon numbers R_{n,k}^(m): A010054 (m=0), A001263 (m=1), this sequence (m=2), A173020 (m=3).

Programs

  • Magma
    A108767:= func< n,k,m | Binomial(n,k)*Binomial(m*n,k-1)/n >;
    [A108767(n,k,2): k in [1..n], n in [1..12]]; // G. C. Greubel, Feb 20 2021
  • Maple
    T:=(n,k)->binomial(n,k)*binomial(2*n,k-1)/n: for n from 1 to 10 do seq(T(n,k),k=1..n) od; # yields sequence in triangular form
  • Mathematica
    T[n_, k_] := Binomial[n, k]*Binomial[2*n, k - 1]/n;
    Table[T[n, k], {n, 1, 10}, {k, 1, n}] // Flatten (* Jean-François Alcover, Nov 11 2017, from Maple *)
  • PARI
    T(n,k) = binomial(n, k)*binomial(2*n, k-1)/n; \\ Andrew Howroyd, Nov 06 2017
    
  • Python
    from sympy import binomial
    def T(n, k): return binomial(n, k)*binomial(2*n, k - 1)//n
    for n in range(1, 21): print([T(n, k) for k in range(1, n + 1)]) # Indranil Ghosh, Nov 07 2017
    
  • R
    T <- function(n, k) {
      choose(n, k)*choose(2*n, k - 1)/n
    } # Indranil Ghosh, Nov 07 2017
    
  • Sage
    def A108767(n,k,m): return binomial(n,k)*binomial(m*n,k-1)/n
    flatten([[A108767(n,k,2) for k in (1..n)] for n in (1..12)]) # G. C. Greubel, Feb 20 2021
    

Formula

T(n, k) = binomial(n, k)*binomial(2*n, k-1)/n.
T(n, n) = A000108(n) (the Catalan numbers).
Sum_{k=1..n} k*T(n,k) = A025174(n).
G.f.: T-1, where T = T(t, z) satisfies T = 1 + z*T^2*(T - 1 + t).
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)) = Limit_{n -> oo} 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).
The o.g.f. for the array of Narayana numbers A001263 is I(1 + t*x + t*x^2 + t*x^3 + ...) = 1 + t*x + (t + t^2)*x^2 + (t + 3*t^2 + t^3)*x^3 + ... . The o.g.f. for the current array is IoI(1 + t*x + t*x^2 + t*x^3 + ...) = 1 + t*x + (t + 2*t^2)*x^2 + (t + 6*t^2 + 5*t^3)*x^3 + ... . Cf. A132081 and A141618. Alternatively, the o.g.f. of this array can be obtained from a single application of I, namely, form I(1 + t*x^2 + t*x^4 + t*x^6 + ...) = 1 + t*x^2 + (t + 2*t^2)*x^4 + (t + 6*t^2 + 5*t^3)*x^6 + ... and then replace x by sqrt(x). This is a particular case of the general result that forming the n-fold composition I^(n)(f(x)) and then replacing x with x^n produces the same result as I(f(x^n)). (End)
O.g.f. is series reversion with respect to x of x/((1+x)*(1+x*u)^2). Cf. A173020. - Peter Bala, Sep 12 2012
n-th row polynomial = x * hypergeom([1 - n, -2*n], [2], x). - Peter Bala, Aug 30 2023

A354622 Irregular triangle read by rows: Refined 3-Narayana triangle. Coefficients of partition polynomials of A134264, a refined Narayana triangle enumerating noncrossing partitions, with all h_k other than h_0, h_3, h_6, ..., h_(3n), ... replaced by zero.

Original entry on oeis.org

1, 1, 3, 1, 9, 12, 1, 12, 6, 66, 55, 1, 15, 15, 105, 105, 455, 273, 1, 18, 18, 9, 153, 306, 51, 816, 1224, 3060, 1428, 1, 21, 21, 21, 210, 420, 210, 210, 1330, 3990, 1330, 5985, 11970, 20349, 7752, 1, 24, 24, 24, 12, 276, 552, 552, 276, 276, 2024, 6072, 3036, 6072, 506, 10626, 42504, 21252, 42504, 106260, 134596, 43263
Offset: 1

Views

Author

Tom Copeland, Jul 08 2022

Keywords

Comments

A set of partition polynomials with these coefficients and the polynomials of A338135 can be generated by substitution of the refined Narayana, or noncrossing partition, polynomials N_n[h_1,...,h_n] of A134264 (h_0=1) into themselves--once for A338135 and twice for this entry--or by setting the indeterminates h_n of A134264 to zero except for h_0, h_3, h_6, ..., h_(3n), ... with h_0 = 1 and ultimately re-indexing. This is equivalent to recursive use of the Lagrange inversion formula on f(x) = x / h(x) = x / (1 + h_1 x + h_2 x^2 + ...) since its compositional inverse is f^{(-1)}(x) = x + N_1(h_1) x + N_2(h_1,h_2) x^2 + .... The equivalence of the two methods of generation--the substitution and the zeroing out--follows from the general theorems stated by Peter Bala in his presentation of formulas for A108767 in 2008, which stem from a fixed point-iteration formalism of a basic identity for a compositional inverse pair, x* h(f^{(-1)}(x)) = f^{(-1)}(x), where, as above, h(x) = x / f(x).
The sets of refined m-Narayana polynomials are used by Cachazo and Umbert to characterize the scattering amplitudes of a class of quantum fields (see, e.g., section 7.3).
These could also be called the refined 3-Dyck path polynomials. From the interpretation of A134264 as Dyck paths in A125181, or staircases whose steps never rise above the diagonal of a square grid (see illustrations in Weisstein), the monomials of the partition polynomial N_4 = 1 (4') + 4 (1') (3') + 2 (2')^2 + 6 (1')^2 (2') + 1 (1')^4 of A134264 have the following correspondences:
1 (4') --> 1 staircase of one step of height 4,
4 (1') (3') --> 4 staircases of 1 step of height 1 and 1 step of height 3,
2 (2')^2 --> 2 staircases of 2 steps of height 2,
6 (1')^2 (2') --> 6 staircases of 2 steps of height 1 and 1 step of height 2,
1 (1')^4 --> 1 staircase of 4 steps of height 1.
Consequently, the partition polynomials G_{3n} of this entry enumerate staircases of height 3n with steps of heights 3, 6, 9, ..., 3k, ... only.
Diverse combinatorial models of the refined m-Narayana, or m-Dyck, polynomials are inherited from those presented for the refined Narayana, or noncrossing partition, polynomials in A134264 and A125181 and in the references therein.
A127537 gives a combinatorial model (see title and Domb and Barret therein, Table 2, p. 355) that contains the coefficients of the monomials h_1^n and h_1^(n-2) h_2, i.e., A001764 and A003408.

Examples

			Triangle begins:
  1;
  1,  3;
  1,  9, 12;
  1, 12,  6, 66,  55;
  1, 15, 15, 105, 105, 455, 273;
  ...
Row 1: G_3  = g_3
row 2: G_6  = g_6 + 3 g_3^2
row 3: G_9  = g_9 + 9 g_3 g_6 + 12 g_3^3
row 4: G_12 = g_12 + 12 g_3 g_9 + 6 g_6^2 + 66 g_3^2 g_6 + 55 g_3^4
row 5: G_15 = g_15 + 15 g_3 g_12 + 15 g_6 g_9 + 105 g_3^2 g_9 + 105 g_3 g_6^2
              + 455 g_3^3 g_6 + 273 g_3^5.
.
In the notation of Abramowitz and Stegun p. 831 with indices of the partitions above divided by 3 and partition indeterminates h_n denoted (n):
R_1 = (1);
R_2 = (2) + 3 (1)^2;
R_3 = (3) + 9 (1) (2) + 12 (1)^3;
R_4 = (4) + 12 (1) (3) + 6 (2)^2 + 66 (1)^2 (2) + 55 (1)^4;
R_5 = (5) + 15 (1) (4) + 15 (2) (3) + 105 (1)^2 (3) + 105 (1) (2)^2 + 455 (1)^3(2)
          + 273 (1)^5.
		

Crossrefs

The length of row n is equal to A000041(n).
Row sums give A002293, n >= 1.

Programs

  • Mathematica
    Table[Binomial[Total[y], Length[y]-1] (Length[y]-1)! / Product[Count[y, i]!, {i, Max@@y}], {n, 8}, {y, Sort[Sort /@ IntegerPartitions[3n, n, Range[3, 3n, 3]]]}] // Flatten (* Andrey Zabolotskiy, Feb 19 2024, using Gus Wiseman's code for A134264 *)
  • PARI
    \\ Compare with A134264
    C(v)={my(n=vecsum(v), S=Set(v)); n!/((n-#v+1)!*prod(i=1, #S, my(x=S[i]); (#select(y->y==x, v))!))}
    row(n)=[C(3*Vec(p)) | p<-partitions(n)]
    { for(n=1, 7, print(row(n))) } \\ Andrew Howroyd, Feb 19 2024

Formula

Coefficients of the monomials are those of the surviving monomials of the partition polynomials of A134264 after zeroing all indeterminates other than h_0, h_3, h_6, h_9, ..., h_(3n), .... The multinomial coefficients of A125181 also apply for G_n, giving the coefficient of the monomial h_1^e_1 h_2^e_2 ... h_n^n of R_n with se := e_1 + e_2 + ... + e_n as (3n)! / ((3n-se+1)! (e_1)! (e_2)! ... (e_n)!).
1*e_1 + 2*e_2 + ... + n*e_n = n for each monomial of R_n.
The partition polynomials R_n = N_n^3 of this entry can be determined from those of A338135, N_n^2, by substituting the partition polynomials of A134264, N_n, for the indeterminate h_n = (n) of N_n^2 or by doing the same for A134264 twice. E.g., N_1(h_1) = h_1, N_2(h_1,h_2) = h_2 + h_1^2, so N_2^2(h_1,h_2) = N_2(N_1,N_2) = N_2 + N_1 = h_2 + h_1^2 + h_1^2 = h_2 + 2 h_1^2 and N_2^3(h_1,h_2) = N_2^2(N_1,N_2) = N_2 + 2 N_1^2 = h_2 + h_1^2 + 2 h_1^2 = h_2 + 3 h_1^2.
Reduces with all indeterminates h_n = (n) = t to A173020.
The coefficient of the monomial h_1^n is (3*n)! / ((3*n-n+1)! n!) = A001764(n) (see also A179848 and A235534). In general, the coefficients of these monomials of the refined (m+1)-Narayana polynomials are the Fuss-Catalan sequence associated to the row sums of the refined m-Narayana polynomials.
The coefficient of the monomial h_1^(n-2) h_2 is (3n)! / ((3n-n+2)! (n-2)!) = A003408(n-2) for n > 1.
The coefficient of the monomial h_1^(n-3) h_3 is (3n)! / ((3n-n+3)! (n-3)!) = A004321(n) for n > 2.

Extensions

Rows 6-8 added by Andrey Zabolotskiy, Feb 19 2024
Showing 1-3 of 3 results.