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-10 of 72 results. Next

A251576 E.g.f.: exp(6*x*G(x)^5) / G(x)^5 where G(x) = 1 + x*G(x)^6 is the g.f. of A002295.

Original entry on oeis.org

1, 1, 6, 96, 2736, 115056, 6455376, 454666176, 38610711936, 3842344221696, 438721154343936, 56549927146392576, 8123473514799876096, 1287034084022760677376, 222964032114987212998656, 41930788886197036399190016, 8507629742037888427727486976, 1852490637585980898960109142016
Offset: 0

Views

Author

Paul D. Hanna, Dec 06 2014

Keywords

Examples

			E.g.f.: A(x) = 1 + x + 6*x^2/2! + 96*x^3/3! + 2736*x^4/4! + 115056*x^5/5! +...
such that A(x) = exp(6*x*G(x)^5) / G(x)^5
where G(x) = 1 + x*G(x)^6 is the g.f. of A002295:
G(x) = 1 + x + 6*x^2 + 51*x^3 + 506*x^4 + 5481*x^5 + 62832*x^6 +...
Note that
A'(x) = exp(6*x*G(x)^5) = 1 + 6*x + 96*x^2/2! + 2736*x^3/3! +...
LOGARITHMIC DERIVATIVE.
The logarithm of the e.g.f. begins:
log(A(x)) = x + 5*x^2/2 + 40*x^3/3 + 385*x^4/4 + 4095*x^5/5 + 46376*x^6/6 +...
and so A'(x)/A(x) = G(x)^5.
TABLE OF POWERS OF E.G.F.
Form a table of coefficients of x^k/k! in A(x)^n as follows.
n=1: [1, 1,   6,   96,  2736,  115056,   6455376,  454666176, ...];
n=2: [1, 2,  14,  228,  6456,  268992,  14968224, 1047087648, ...];
n=3: [1, 3,  24,  402, 11376,  470808,  26011584, 1808151552, ...];
n=4: [1, 4,  36,  624, 17736,  730944,  40143456, 2774490624, ...];
n=5: [1, 5,  50,  900, 25800, 1061400,  58017600, 3989340000, ...];
n=6: [1, 6,  66, 1236, 35856, 1475856,  80395056, 5503484736, ...];
n=7: [1, 7,  84, 1638, 48216, 1989792, 108156384, 7376303088, ...];
n=8: [1, 8, 104, 2112, 63216, 2620608, 142314624, 9676910592, ...]; ...
in which the main diagonal begins (see A251586):
[1, 2, 24, 624, 25800, 1475856, 108156384, 9676910592, ...]
and is given by the formula:
[x^n/n!] A(x)^(n+1) = 6^(n-4) * (n+1)^(n-5) * (125*n^4 + 810*n^3 + 2095*n^2 + 2586*n + 1296) for n>=0.
		

Crossrefs

Programs

  • Mathematica
    Flatten[{1,1,Table[Sum[6^k * n!/k! * Binomial[6*n-k-6, n-k] * (k-1)/(n-1),{k,0,n}],{n,2,20}]}] (* Vaclav Kotesovec, Dec 07 2014 *)
  • PARI
    {a(n) = local(G=1);for(i=1,n,G=1+x*G^6 +x*O(x^n)); n!*polcoeff(exp(6*x*G^5)/G^5, n)}
    for(n=0, 20, print1(a(n), ", "))
    
  • PARI
    {a(n) = if(n==0||n==1, 1, sum(k=0, n, 6^k * n!/k! * binomial(6*n-k-6,n-k) * (k-1)/(n-1) ))}
    for(n=0, 20, print1(a(n), ", "))

Formula

Let G(x) = 1 + x*G(x)^6 be the g.f. of A002295, then the e.g.f. A(x) of this sequence satisfies:
(1) A'(x)/A(x) = G(x)^5.
(2) A'(x) = exp(6*x*G(x)^5).
(3) A(x) = exp( Integral G(x)^5 dx ).
(4) A(x) = exp( Sum_{n>=1} A130564(n)*x^n/n ), where A130564(n) = binomial(6*n-2,n)/(5*n-1).
(5) A(x) = F(x/A(x)) where F(x) is the e.g.f. of A251586.
(6) A(x) = Sum_{n>=0} A251586(n)*(x/A(x))^n/n! and
(7) [x^n/n!] A(x)^(n+1) = (n+1)*A251586(n),
where A251586(n) = 6^(n-4) * (n+1)^(n-6) * (125*n^4 + 810*n^3 + 2095*n^2 + 2586*n + 1296).
a(n) = Sum_{k=0..n} 6^k * n!/k! * binomial(6*n-k-6, n-k) * (k-1)/(n-1) for n>1.
Recurrence: 5*(5*n-9)*(5*n-8)*(5*n-7)*(5*n-6)*(9*n^4 - 99*n^3 + 413*n^2 - 777*n + 559)*a(n) = 72*(5832*n^9 - 113724*n^8 + 986580*n^7 - 5003586*n^6 + 16373448*n^5 - 35916483*n^4 + 52931854*n^3 - 50678109*n^2 + 28701206*n - 7357350)*a(n-1) + 46656*(9*n^4 - 63*n^3 + 170*n^2 - 212*n + 105)*a(n-2). - Vaclav Kotesovec, Dec 07 2014
a(n) ~ 6^(6*(n-1)-1/2) / 5^(5*(n-1)-1/2) * n^(n-2) / exp(n-1). - Vaclav Kotesovec, Dec 07 2014

A251666 E.g.f.: exp(6*x*G(x)^5) / G(x) where G(x) = 1 + x*G(x)^6 is the g.f. of A002295.

Original entry on oeis.org

1, 5, 74, 2028, 83352, 4607496, 321156000, 27064420704, 2677510124928, 304299947999232, 39075730095810816, 5595805388119057920, 884245579070535235584, 152843879008651568329728, 28688663318934190485491712, 5811091829207760774331662336, 1263471121829937070180445552640
Offset: 0

Views

Author

Paul D. Hanna, Dec 07 2014

Keywords

Examples

			E.g.f.: A(x) = 1 + 5*x + 74*x^2/2! + 2028*x^3/3! + 83352*x^4/4! + 4607496*x^5/5! +...
such that A(x) = exp(6*x*G(x)^5) / G(x)
where G(x) = 1 + x*G(x)^6 is the g.f. of A002295:
G(x) = 1 + x + 6*x^2 + 51*x^3 + 506*x^4 + 5481*x^5 + 62832*x^6 +...
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[6^k * n!/k! * Binomial[6*n-k-2,n-k] * (5*k-1)/(5*n-1),{k,0,n}],{n,0,20}] (* Vaclav Kotesovec, Dec 07 2014 *)
  • PARI
    {a(n)=local(G=1); for(i=0, n, G=1+x*G^6 +x*O(x^n)); n!*polcoeff(exp(6*x*G^5)/G, n)}
    for(n=0, 20, print1(a(n), ", "))
    
  • PARI
    {a(n) = sum(k=0, n, 6^k * n!/k! * binomial(6*n-k-2,n-k) * (5*k-1)/(5*n-1) )}
    for(n=0, 20, print1(a(n), ", "))

Formula

Let G(x) = 1 + x*G(x)^6 be the g.f. of A002295, then the e.g.f. A(x) of this sequence satisfies:
(1) A'(x)/A(x) = G(x)^5 + 4*G'(x)/G(x).
(2) A(x) = F(x/A(x)^5) where F(x) is the e.g.f. of A251696.
(3) A(x) = Sum_{n>=0} A251696(n)*(x/A(x)^5)^n/n! where A251696(n) = (4*n+1) * (5*n+1)^(n-2) * 6^n .
(4) [x^n/n!] A(x)^(5*n+1) = (4*n+1) * (5*n+1)^(n-1) * 6^n .
a(n) = Sum_{k=0..n} 6^k * n!/k! * binomial(6*n-k-2,n-k) * (5*k-1)/(5*n-1) for n>=0.
Recurrence: 5*(5*n-4)*(5*n-3)*(5*n-2)*(5*n-1)*(1296*n^5 - 11394*n^4 + 40230*n^3 - 71274*n^2 + 63110*n - 21963)*a(n) = 144*(419904*n^10 - 5161320*n^9 + 28223964*n^8 - 90513612*n^7 + 188713962*n^6 - 267339204*n^5 + 259905051*n^4 - 169257762*n^3 + 67929146*n^2 - 12957136*n - 43050)*a(n-1) + 46656*(1296*n^5 - 4914*n^4 + 7614*n^3 - 5988*n^2 + 2156*n + 5)*a(n-2). - Vaclav Kotesovec, Dec 07 2014
a(n) ~ 4 * 6^(6*n-3/2) / 5^(5*n-1/2) * n^(n-1) / exp(n-1). - Vaclav Kotesovec, Dec 07 2014

A226705 G.f.: 1 / (1 + 12*x*G(x)^4 - 16*x*G^5) where G(x) = 1 + x*G(x)^6 is the g.f. of A002295.

Original entry on oeis.org

1, 4, 48, 600, 7856, 105684, 1447392, 20075416, 281086416, 3964453368, 56240518128, 801624722232, 11470976280960, 164691196943212, 2371222443727584, 34224696393237360, 495036708728067088, 7173892793100898728, 104135761805147016096, 1513892435551302963792
Offset: 0

Views

Author

Paul D. Hanna, Jun 15 2013

Keywords

Examples

			G.f.: A(x) = 1 + 4*x + 48*x^2 + 600*x^3 + 7856*x^4 + 105684*x^5 +...
A related series is G(x) = 1 + x*G(x)^6, where
G(x) = 1 + x + 6*x^2 + 51*x^3 + 506*x^4 + 5481*x^5 + 62832*x^6 +...
G(x)^4 = 1 + 4*x + 30*x^2 + 280*x^3 + 2925*x^4 + 32736*x^5 +...
G(x)^5 = 1 + 5*x + 40*x^2 + 385*x^3 + 4095*x^4 + 46376*x^5 +...
such that A(x) = 1/(1 + 12*x*G(x)^4 - 16*x*G^5).
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[3*n+2*k,n-k]*Binomial[3*n-2*k,k],{k,0,n}],{n,0,20}] (* Vaclav Kotesovec, Jun 16 2013 *)
  • PARI
    {a(n)=local(G=1+x); for(i=0, n,G=1+x*G^6+x*O(x^n)); polcoeff(1/(1+12*x*G^4-16*x*G^5), n)}
    for(n=0, 30, print1(a(n), ", "))
    
  • PARI
    {a(n)=local(G=1+x); for(i=0, n,G=1+x*G^6+x*O(x^n)); polcoeff(1/(1-4*x*G^4-16*x^2*G^10), n)}
    for(n=0, 30, print1(a(n), ", "))
    
  • PARI
    {a(n)=sum(k=0, n, binomial(2*k, n-k)*binomial(6*n-2*k, k))}
    for(n=0, 30, print1(a(n), ", "))
    
  • PARI
    {a(n)=sum(k=0, n, binomial(3*n +2*k, n-k)*binomial(3*n-2*k, k))}
    for(n=0, 30, print1(a(n), ", "))
    
  • PARI
    {a(n)=sum(k=0, n, binomial(6*n +2*k, n-k)*binomial(-2*k, k))}
    for(n=0, 30, print1(a(n), ", "))

Formula

a(n) = Sum_{k=0..n} C(2*k, n-k) * C(6*n-2*k, k).
a(n) = Sum_{k=0..n} C(n+2*k, n-k) * C(5*n-2*k, k).
a(n) = Sum_{k=0..n} C(2*n+2*k, n-k) * C(4*n-2*k, k).
a(n) = Sum_{k=0..n} C(3*n+2*k, n-k) * C(3*n-2*k, k).
a(n) = Sum_{k=0..n} C(4*n+2*k, n-k) * C(2*n-2*k, k).
a(n) = Sum_{k=0..n} C(5*n+2*k, n-k) * C(n-2*k, k).
a(n) = Sum_{k=0..n} C(6*n+2*k, n-k) * C(-2*k, k).
Self-convolution of A226706.
G.f.: 1 / (1 - 4*x*G(x)^4 - 16*x^2*G(x)^10) where G(x) = 1 + x*G(x)^6 is the g.f. of A002295.
a(n) ~ 2^(6*n-2)*3^(6*n+3/2)/(5^(5*n+1/2)*sqrt(Pi*n)). - Vaclav Kotesovec, Jun 16 2013
From Seiichi Manyama, Aug 05 2025: (Start)
a(n) = [x^n] 1/((1+2*x) * (1-x)^(5*n+1)).
a(n) = Sum_{k=0..n} (-3)^(n-k) * binomial(6*n+1,k).
a(n) = Sum_{k=0..n} (-2)^(n-k) * binomial(5*n+k,k). (End)
From Seiichi Manyama, Aug 14 2025: (Start)
a(n) = Sum_{k=0..n} (-2)^k * 3^(n-k) * binomial(6*n+1,k) * binomial(6*n-k,n-k).
G.f.: G(x)^2/((-2+3*G(x)) * (6-5*G(x))) where G(x) = 1+x*G(x)^6 is the g.f. of A002295. (End)
G.f.: B(x)^2/(1 + 4*(B(x)-1)/3), where B(x) is the g.f. of A004355. - Seiichi Manyama, Aug 15 2025

A186183 Expansion of 1/(1-x*A002295(x)).

Original entry on oeis.org

1, 1, 2, 9, 68, 646, 6857, 77695, 919642, 11233858, 140544189, 1791614714, 23187320736, 303861373679, 4023883823059, 53762917329659, 723854999871943, 9811154512175468, 133762940465746744, 1833187046654598058, 25239961633188882896
Offset: 0

Views

Author

Vladimir Kruchinin, Feb 14 2011

Keywords

Programs

  • Maple
    a:= n-> `if` (n=0, 1, add (k/(5*n-4*k) *binomial (6*n-5*k-1, n-k), k=1..n)):
    seq (a(n), n=0..30);
  • PARI
    a(n)=max(1,sum(k=1,n, k/(5*n-4*k)*binomial(6*n-5*k-1,n-k)))

Formula

a(n) = Sum_{k=1..n} k/(5*n-4*k) * binomial(6*n-5*k-1,n-k) if n>0; a(0)=1.
a(n) ~ 2^(6*n+4) * 3^(6*n + 9/2) / (51136801 * sqrt(Pi) * n^(3/2) * 5^(5*n - 7/2)). - Vaclav Kotesovec, Sep 22 2024

A226706 G.f.: 1 / sqrt(1 + 12*x*G(x)^4 - 16*x*G(x)^5) where G(x) = 1 + x*G(x)^6 is the g.f. of A002295.

Original entry on oeis.org

1, 2, 22, 256, 3174, 40862, 539376, 7247448, 98684230, 1357638124, 18831752122, 262974273200, 3692853486768, 52102851020154, 738102882420440, 10492839572260176, 149623214762194182, 2139329701502229300, 30661862088900836964, 440404155129948147776
Offset: 0

Views

Author

Paul D. Hanna, Jun 15 2013

Keywords

Examples

			G.f.: A(x) = 1 + 2*x + 22*x^2 + 256*x^3 + 3174*x^4 + 40862*x^5 +...
A related series is G(x) = 1 + x*G(x), which begins
G(x) = 1 + x + 6*x^2 + 51*x^3 + 506*x^4 + 5481*x^5 + 62832*x^6 +...
where A(x) = 1/sqrt(1 + 12*x*G(x)^4 - 16*x*G(x)^5).
		

Crossrefs

Programs

  • PARI
    {a(n)=local(G=1+x); for(i=0, n,G=1+x*G^6+x*O(x^n)); polcoeff(1/sqrt(1+12*x*G^4-16*x*G^5), n)}
    for(n=0, 30, print1(a(n), ", "))

Formula

Sum_{k=0..n} a(n-k)*a(k) = Sum_{k=0..n} C(3*n+2*k,n-k)*C(3*n-2*k,k).
Self-convolution equals A226705.

A079463 Numbers k such that C(6*k,n)/(5*k+1) (A002295) is not divisible by 6.

Original entry on oeis.org

0, 1, 3, 4, 5, 7, 13, 16, 19, 21, 22, 34, 40, 49, 51, 53, 67, 77, 83, 85, 97, 103, 115, 121, 148, 196, 202, 205, 211, 213, 292, 307, 309, 310, 333, 339, 340, 341, 346, 358, 364, 439, 445, 583, 589, 601, 607, 819, 821, 845, 851, 853, 877, 925, 931, 1021, 1039, 1069
Offset: 1

Views

Author

Benoit Cloitre, Jan 14 2003

Keywords

Comments

Is a(n)/n^2 bounded? It seems that 0.2 < a(n)/n^2 <= 1.

Crossrefs

Cf. A000225 (C(2*k, k)/(k+1) is not divisible by 2), A003462 (C(3*k, k)/(2*k+1) is not divisible by 3), A003463 (C(5*k, k)/(4*k+1) is not divisible by 5).

Programs

  • Mathematica
    Select[Range[0,1000], !Divisible[Binomial[6*#,#]/(5*# + 1),6] &] (* Amiram Eldar, Jun 04 2021 *)

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

A002294 a(n) = binomial(5*n, n)/(4*n + 1).

Original entry on oeis.org

1, 1, 5, 35, 285, 2530, 23751, 231880, 2330445, 23950355, 250543370, 2658968130, 28558343775, 309831575760, 3390416787880, 37377257159280, 414741863546285, 4628362722856425, 51912988256282175, 584909606696793885, 6617078646960613370
Offset: 0

Views

Author

Keywords

Comments

From Wolfdieter Lang, Sep 14 2007: (Start)
a(n), n >= 1, enumerates quintic trees (rooted, ordered, incomplete) with n vertices (including the root).
This is the Pfaff-Fuss-Catalan sequence C^{m}_n for m = 5. See the Graham et al. reference, p. 347. eq. 7.66. See also the Pólya-Szegő reference.
Also 5-Raney sequence. See the Graham et al. reference, pp. 346-347. (End)
a(n) = A258708(3*n, 2*n) for n > 0. - Reinhard Zumkeller, Jun 23 2015
Conjecturally, a(n) is the number of 4-uniform words on the alphabet [n] that avoid the patterns 231 and 221 (see the Defant and Kravitz link). - Colin Defant, Sep 26 2018
From Stillwell (1995), p. 62: "Eisenstein's Theorem. If y^5 + y = x, then y has a power series expansion y = x - x^5 + 10*x^9/2^1 - 15 * 14 * x^13/3! + 20 * 19 * 18*x^17/4! - ...." - Michael Somos, Sep 19 2019
a(n) is the total number of down steps before the first up step in all 4_1-Dyck paths of length 5*n. A 4_1-Dyck path is a lattice path with steps (1, 4), (1, -1) that starts and ends at y = 0 and stays above the line y = -1. - Sarah Selkirk, May 10 2020
Dropping the first 1 (starting from 1, 5, 35, ... with offset 1), the series reversion gives 1, -5, 15, -35, 70, ... (again offset 1), essentially A000332 and row 5 of A027555. - R. J. Mathar, Aug 17 2023
Number of rooted polyominoes composed of n hexagonal cells of the hyperbolic regular tiling with Schläfli symbol {6,oo}. A rooted polyomino has one external edge identified, and chiral pairs are counted as two. A stereographic projection of the {6,oo} tiling on the Poincaré disk can be obtained via the Christensson link. - Robert A. Russell, Jan 27 2024
This is instance k = 5 of the generalized Catalan family {C(k, n)}_{n>=0} given in a comment of A130564. - Wolfdieter Lang, Feb 05 2024

Examples

			There are a(2) = 5 quintic trees (vertex degree <= 5 and 5 possible branchings) with 2 vertices (one of them the root). Adding one more branch (one more vertex) to these five trees yields 5*5 + binomial(5,2) = 35 = a(3) such trees.
G.f. = 1 + x + 5*x^2 + 35*x^3 + 285*x^4 + 2530*x^5 + 23751*x^6 + 231880*x^7 + ...
G.f. = t + t^5 + 5*t^9 + 35*t^13 + 285*t^17 + 2530*t^21 + 23751*t^25 + 231880*t^29 + ...
		

References

  • Archiv der Mathematik u. Physik, Editor's note: "Über die Bestimmung der Anzahl der verschiedenen Arten, auf welche sich ein n-Eck durch Diagonalen in lauter m-Ecke zerlegen laesst, mit Bezug auf einige Abhandlungen der Herren Lame, Rodrigues, Binet, Catalan und Duhamel in dem Journal de Mathematiques pures et appliquees, publie par Joseph Liouville. T. III. IV.", Archiv der Mathematik u. Physik, 1 (1841), pp. 193ff; see especially p. 198.
  • 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.
  • 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.
  • Ulrike Sattler, Decidable classes of formal power series with nice closure properties, Diplomarbeit im Fach Informatik, Univ. Erlangen - Nürnberg, Jul 27 1994.
  • 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

Cf. A001764, A002296, A258708, A346647 (binomial transform), A346665 (inverse binomial transform).
Fourth column of triangle A062993.
Polyominoes: A221184{n-1} (oriented), A004127 (unoriented), A369473 (chiral), A143546 (achiral), A002293 {5,oo}, A002295 {7,oo}.
Cf. A130564.

Programs

  • GAP
    List([0..22],n->Binomial(5*n,n)/(4*n+1)); # Muniru A Asiru, Nov 01 2018
  • Haskell
    a002294 n = a002294_list !! n
    a002294_list = [a258708 (3 * n) (2 * n) | n <- [1..]]
    -- Reinhard Zumkeller, Jun 23 2015
    
  • Magma
    [ Binomial(5*n,n)/(4*n+1): n in [0..100]]; // Vincenzo Librandi, Mar 24 2011
    
  • Maple
    seq(binomial(5*k+1,k)/(5*k+1),k=0..30); # Robert FERREOL, Apr 03 2015
    n:=30:G:=series(RootOf(g = 1+x*g^5, g),x=0,n+1):seq(coeff(G,x,k),k=0..n); # Robert FERREOL, Apr 03 2015
  • Mathematica
    CoefficientList[InverseSeries[ Series[ y - y^5, {y, 0, 100}], x], x][[Range[2, 100, 4]]]
    Table[Binomial[5n,n]/(4n+1),{n,0,20}] (* Harvey P. Dale, Dec 30 2011 *)
    a[ n_] := SeriesCoefficient[ HypergeometricPFQ[ {1, 2, 3, 4}/5, {2, 3, 5}/4, x 5^5/4^4], {x, 0, n}]; (* Michael Somos, May 06 2015 *)
    a[ n_] := With[{m = 4 n + 1}, SeriesCoefficient[ InverseSeries @ Series[ x - x^5, {x, 0, m}], {x, 0, m}]]; (* Michael Somos, May 06 2015 *)
  • PARI
    {a(n) = binomial( 5 * n, n) / (4*n + 1)}; /* Michael Somos, Mar 17 2011 */
    
  • PARI
    {a(n) = if( n<0, 0, n = 4*n + 1; polcoeff( serreverse( x - x^5 + x * O(x^n) ), n))}; /* Michael Somos, Mar 17 2011 */
    

Formula

For the connection with the solution of the quintic, hypergeometric series, and Lagrange inversion, see Beukers (2014). - N. J. A. Sloane, Mar 12 2014
G.f.: hypergeometric([1, 2, 3, 4] / 5, [2, 3, 5] / 4, x * 5^5 / 4^4). - Michael Somos, Mar 17 2011
O.g.f. A(x) satisfies A(x) = 1 + x * A(x)^5 = 1 / (1 - x * A(x)^4).
Given g.f. A(x) then z = t * A(t^4) satisfies 0 = z^5 - z + t. - Michael Somos, Mar 17 2011
a(n) = binomial(5*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.
a(n) = upper left term in M^n, M = the production matrix:
1, 1;
4, 4, 1;
10, 10, 4, 1;
20, 20, 10, 4, 1;
...
where (1, 4, 10, 20, ...) is the tetrahedral sequence, A000292. - Gary W. Adamson, Jul 08 2011
D-finite with recurrence: 8*n*(4*n+1)*(2*n-1)*(4*n-1)*a(n) - 5*(5*n-4)*(5*n-3)*(5*n-2)*(5*n-1)*a(n-1) = 0. - R. J. Mathar, Dec 02 2014
a(n) = binomial(5*n + 1, n)/(5*n + 1) = A062993(n+3,3). - Robert FERREOL, Apr 03 2015
a(0) = 1; a(n) = Sum_{i1 + i2 + ... + i5 = n - 1} a(i1) * a(i2) * ... *a(i5) for n >= 1. - Robert FERREOL, Apr 03 2015
From Ilya Gutkovskiy, Jan 15 2017: (Start)
O.g.f.: 5F4([1/5, 2/5, 3/5, 4/5, 1]; [1/2, 3/4, 1, 5/4]; 3125*x/256).[Cancellation of the 1s, see G.f. the above. - Wolfdieter Lang, Feb 05 2024]
E.g.f.: 4F4([1/5, 2/5, 3/5, 4/5]; [1/2, 3/4, 1, 5/4]; 3125*x/256).
a(n) ~ 5^(5*n + 1/2)/(sqrt(Pi) * 2^(8*n + 7/2) * n^(3/2)). (End)
x*A'(x)/A(x) = (A(x) - 1)/(- 4*A(x) + 5) = x + 9*x^2 + 91*x^3 + 969*x^4 + ... is the o.g.f. of A163456. Cf. A001764 and A002293 - A002296. - Peter Bala, Feb 04 2022
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^9). - Seiichi Manyama, Jun 16 2025

Extensions

More terms from Olivier Gérard, Jul 05 2001

A002296 Number of dissections of a polygon: binomial(7n,n)/(6n+1).

Original entry on oeis.org

1, 1, 7, 70, 819, 10472, 141778, 1997688, 28989675, 430321633, 6503352856, 99726673130, 1547847846090, 24269405074740, 383846168712104, 6116574500860880, 98106248306858715, 1582638261961640247, 25661404527790252375, 417980115131315136400
Offset: 0

Views

Author

Keywords

Comments

a(n), n>=1, enumerates heptic (7-ary) trees (rooted, ordered, incomplete) with n vertices (including the root).
Pfaff-Fuss-Catalan sequence C^{m}_n for m=7. See the Graham et al. reference, p. 347. eq. 7.66. See also the Pólya-Szegő reference.
Also 7-Raney sequence. See the Graham et al. reference, pp. 346-347.
a(n) = A258708(3*n,2*n) for n > 0. - Reinhard Zumkeller, Jun 23 2015
This is instance k = 7 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

Examples

			There are a(2)=7 heptic trees (vertex degree <= 7 and 7 possible branchings) with 2 vertices (one of them the root). Adding one more branch (one more vertex) to these 7 trees yields 7*7 + binomial(7,2) = 70 = a(3) such trees.
		

References

  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, pp. 200, 347.
  • 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.
  • Ulrike Sattler, Decidable classes of formal power series with nice closure properties, Diplomarbeit im Fach Informatik, Univ. Erlangen - Nürnberg, Jul 27 1994.
  • 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

Sixth column of triangle A062993.
Cf. A235535: binomial(9n,3n)/(6n+1); A235536: binomial(8n,2n)/(6n+1).
Cf. A258708.
Cf. A130564.

Programs

  • Haskell
    a002296 n = a002296_list !! n
    a002296_list = [a258708 (4 * n) (3 * n) | n <- [1..]]
    -- Reinhard Zumkeller, Jun 23 2015
  • Maple
    seq(binomial(7*n+1, n)/(7*n+1), n=0..30); # Robert FERREOL, Apr 02 2015
    n:=30: G:=series(RootOf(g = 1+x*g^7, g), x=0, n+1): seq(coeff(G, x, k), k=0..n); # Robert FERREOL, Apr 02 2015
  • Mathematica
    Table[Binomial[7n,n]/(6n+1),{n,0,20}] (* Harvey P. Dale, Nov 21 2011 *)
  • PARI
    a(n)=binomial(7*n,n)/(6*n+1) \\ Charles R Greathouse IV, Feb 06 2012
    

Formula

O.g.f. A(x) = 1 + x*A(x)^7 = 1/(1-x*A(x)^6).
a(n) = binomial(7*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.
D-finite with recurrence: 72*n*(6*n-1)*(3*n-1)*(2*n-1)*(3*n-2)*(6*n+1)*a(n) - 7*(7*n-3)*(7*n-6)*(7*n-2)*(7*n-5)*(7*n-1)*(7*n-4)*a(n-1) = 0. - R. J. Mathar, Nov 16 2012
a(n) are special values of Jacobi polynomials, in Maple notation:
a(n) = JacobiP(n-1, 6*n+1, -n, 1)/n, n = 1, 2, ... . - Karol A. Penson, Mar 16 2015
a(n) = binomial(7*n+1, n)/(7*n+1) = A062993(n+5,5). - Robert FERREOL, Apr 02 2015
a(0) = 1; a(n) = Sum_{i1+i2+...+i7=n-1} a(i1)*a(i2)*...*a(i7) for n>=1. - Robert FERREOL, Apr 02 2015
From Ilya Gutkovskiy, Jan 16 2017: (Start)
O.g.f.: 6F5(1/7,2/7,3/7,4/7,5/7,6/7; 1/3,1/2,2/3,5/6,7/6; 823543*x/46656).
E.g.f.: 6F6(1/7,2/7,3/7,4/7,5/7,6/7; 1/3,1/2,2/3,5/6,1,7/6; 823543*x/46656).
a(n) ~ 7^(7*n+1/2)/(sqrt(Pi)*3^(6*n+3/2)*4^(3*n+1)*n^(3/2)). (End)
x*A'(x)/A(x) = (A(x) - 1)/(- 6*A(x) + 7) = x + 13*x^2 + 190*x^3 + 2925*x^4 + ... = (1/7)*Sum_{n >= 1} binomial(7*n,n)*x^n. Cf. A001764 and A002293, A002294, A002295. - Peter Bala, Feb 04 2022
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^13). - Seiichi Manyama, Jun 16 2025

Extensions

Pfaff-Fuss-Catalan, Raney, o.g.f. and 7-ary tree comments from Wolfdieter Lang, Sep 14 2007

A130564 Member k=5 of a family of generalized Catalan numbers.

Original entry on oeis.org

1, 5, 40, 385, 4095, 46376, 548340, 6690585, 83615350, 1064887395, 13770292256, 180320238280, 2386316821325, 31864803599700, 428798445360120, 5809228810425801, 79168272296871450, 1084567603590147950
Offset: 1

Views

Author

Wolfdieter Lang, Jul 13 2007

Keywords

Comments

The generalized Catalan numbers C(k,n):= binomial(k*n+1,n)/(k*n+1) become for negative k=-|k|, with |k|>=2, ((-1)^(n-1))*binomial((|k|+1)*n-2,n)/(|k|*n-1), n>=0.
The family c(k,n):=binomial((k+1)*n-2,n)/(k*n-1), n>=1, has the members A000108, A006013, A006632, A118971 for k=1,2,3,4, respectively (but the offset there is 0).
The members of the C(k,n) family for positive k are: A000012 (powers of 1), A000108, A001764, A002293, A002294, A002295, A002296, A007556, A062994, for k=1..9.

References

  • Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, Addison-Wesley, Reading, MA, 2nd ed. 1994, pp. 200, 363.

Crossrefs

Programs

  • Mathematica
    Rest@ CoefficientList[InverseSeries[Series[y (1 - y)^5, {y, 0, 18}], x], x] (* Michael De Vlieger, Oct 13 2019 *)

Formula

a(n) = binomial((k+1)*n-2,n)/(k*n-1), with k=5.
G.f.: inverse series of y*(1-y)^5.
a(n) = (5/6)*binomial(6*n,n)/(6*n-1). [Bruno Berselli, Jan 17 2014]
From Wolfdieter Lang, Feb 06 2020: (Start)
G.f.: (5/6)*(1 - hypergeom([-1, 1, 2, 3, 4]/6, [1, 2, 3, 4]/5,(6^6/5^5)*x)).
E.g.f.: (5/6)*(1 - hypergeom([-1, 1, 2, 3, 4]/6, [1, 2, 3, 4, 5]/5,(6^6/5^5)*x)). (End)
D-finite with recurrence 5*n*(5*n-4)*(5*n-3)*(5*n-2)*(5*n-1)*a(n) -72*(6*n-7)*(3*n-1)*(2*n-1)*(3*n-2)*(6*n-5)*a(n-1)=0. - R. J. Mathar, May 07 2021
Showing 1-10 of 72 results. Next