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 50 results. Next

A251577 E.g.f.: exp(7*x*G(x)^6) / G(x)^6 where G(x) = 1 + x*G(x)^7 is the g.f. of A002296.

Original entry on oeis.org

1, 1, 7, 133, 4501, 224497, 14926387, 1245099709, 125177105641, 14743403405857, 1991987858095039, 303781606238806549, 51624122993243471293, 9674836841745014156497, 1982441139367342976694379, 440946185623028320815311053, 105810290178441439797537070033, 27247415403508413760437930799681
Offset: 0

Views

Author

Paul D. Hanna, Dec 06 2014

Keywords

Examples

			E.g.f.: A(x) = 1 + x + 7*x^2/2! + 133*x^3/3! + 4501*x^4/4! + 224497*x^5/5! +...
such that A(x) = exp(7*x*G(x)^6) / G(x)^6
where G(x) = 1 + x*G(x)^7 is the g.f. of A002296:
G(x) = 1 + x + 7*x^2 + 70*x^3 + 819*x^4 + 10472*x^5 + 141778*x^6 +...
Note that
A'(x) = exp(7*x*G(x)^6) = 1 + 7*x + 133*x^2/2! + 4501*x^3/3! +...
LOGARITHMIC DERIVATIVE.
The logarithm of the e.g.f. begins:
log(A(x)) = x + 6*x^2/2 + 57*x^3/3 + 650*x^4/4 + 8184*x^5/5 + 109668*x^6/6 +...
and so A'(x)/A(x) = G(x)^6.
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,  7,   133,  4501,  224497,  14926387,  1245099709, ...];
n=2: [1, 2,  16,  308, 10360,  512624,  33845728,  2807075264, ...];
n=3: [1, 3,  27,  531, 17829,  876771,  57529143,  4745597787, ...];
n=4: [1, 4,  40,  808, 27184, 1331008,  86864512,  7129675840, ...];
n=5: [1, 5,  55, 1145, 38725, 1891205, 122869075, 10038831425, ...];
n=6: [1, 6,  72, 1548, 52776, 2575152, 166702752, 13564381824, ...];
n=7: [1, 7,  91, 2023, 69685, 3402679, 219682183, 17810832319, ...];
n=8: [1, 8, 112, 2576, 89824, 4395776, 283295488, 22897384832, ...]; ...
in which the main diagonal begins (see A251587):
[1, 2, 27, 808, 38725, 2575152, 219682183, 22897384832, ...]
and is given by the formula:
[x^n/n!] A(x)^(n+1) = 7^(n-5) * (n+1)^(n-6) * (1296*n^5 + 9720*n^4 + 30555*n^3 + 50665*n^2 + 44621*n + 16807) for n>=0.

Crossrefs

Cf. A251587, A251667, A002296, A130565.
Cf. Variants: A243953, A251573, A251574, A251575, A251576, A251578, A251579, A251580.

Programs

  • Mathematica
    Flatten[{1,1,Table[Sum[7^k * n!/k! * Binomial[7*n-k-7, 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^7 +x*O(x^n)); n!*polcoeff(exp(7*x*G^6)/G^6, n)}
for(n=0, 20, print1(a(n), ", "))
  • PARI
    {a(n) = if(n==0|n==1, 1, sum(k=0, n, 7^k * n!/k! * binomial(7*n-k-7,n-k) * (k-1)/(n-1) ))}
for(n=0, 20, print1(a(n), ", "))

Formula

Let G(x) = 1 + x*G(x)^7 be the g.f. of A002296, then the e.g.f. A(x) of this sequence satisfies:
(1) A'(x)/A(x) = G(x)^6.
(2) A'(x) = exp(7*x*G(x)^6).
(3) A(x) = exp( Integral G(x)^6 dx ).
(4) A(x) = exp( Sum_{n>=1} A130565(n)*x^n/n ), where A130565(n) = binomial(7*n-2,n)/(6*n-1).
(5) A(x) = F(x/A(x)) where F(x) is the e.g.f. of A251587.
(6) A(x) = Sum_{n>=0} A251587(n)*(x/A(x))^n/n! and
(7) [x^n/n!] A(x)^(n+1) = (n+1)*A251587(n),
where A251587(n) = 7^(n-5) * (n+1)^(n-7) * (1296*n^5 + 9720*n^4 + 30555*n^3 + 50665*n^2 + 44621*n + 16807).
a(n) = Sum_{k=0..n} 7^k * n!/k! * binomial(7*n-k-7, n-k) * (k-1)/(n-1) for n>1.
Recurrence: 72*(2*n-3)*(3*n-5)*(3*n-4)*(6*n-11)*(6*n-7)*(2401*n^5 - 32585*n^4 + 178311*n^3 - 492779*n^2 + 689623*n - 392491)*a(n) = 7*(282475249*n^11 - 6658345155*n^10 + 71339412375*n^9 - 458968749330*n^8 + 1971937124661*n^7 - 5947597074909*n^6 + 12867618998885*n^5 - 20002508046570*n^4 + 21938241804255*n^3 - 16207858252075*n^2 + 7281095411817*n - 1512276480000)*a(n-1) - 823543*(2401*n^5 - 20580*n^4 + 71981*n^3 - 129346*n^2 + 120663*n - 47520)*a(n-2). - Vaclav Kotesovec, Dec 07 2014
a(n) ~ 7^(7*(n-1)-1/2) / 6^(6*(n-1)-1/2) * n^(n-2) / exp(n-1). - Vaclav Kotesovec, Dec 07 2014

A251667 E.g.f.: exp(7*x*G(x)^6) / G(x) where G(x) = 1 + x*G(x)^7 is the g.f. of A002296.

Original entry on oeis.org

1, 6, 107, 3508, 171741, 11280842, 933014767, 93212094024, 10925496633401, 1470493880790382, 223555405538724819, 37892802280129883324, 7086076189702624109653, 1449303152891376476830962, 321848482510755456019058519, 77124029495405859198280522768
Offset: 0

Views

Author

Paul D. Hanna, Dec 07 2014

Keywords

Examples

			E.g.f.: A(x) = 1 + 6*x + 107*x^2/2! + 3508*x^3/3! + 171741*x^4/4! + 11280842*x^5/5! +...
such that A(x) = exp(7*x*G(x)^6) / G(x)
where G(x) = 1 + x*G(x)^7 is the g.f. of A002296:
G(x) = 1 + x + 7*x^2 + 70*x^3 + 819*x^4 + 10472*x^5 + 141778*x^6 +...

Crossrefs

Cf. A251577, A251697, A002295.
Cf. Variants: A243953, A251663, A251664, A251665, A251666, A251668, A251669, A251670.

Programs

  • Mathematica
    Table[Sum[7^k * n!/k! * Binomial[7*n-k-2,n-k] * (6*k-1)/(6*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^7 +x*O(x^n)); n!*polcoeff(exp(7*x*G^6)/G, n)}
for(n=0, 20, print1(a(n), ", "))
  • PARI
    {a(n) = sum(k=0, n, 7^k * n!/k! * binomial(7*n-k-2,n-k) * (6*k-1)/(6*n-1) )}
for(n=0, 20, print1(a(n), ", "))

Formula

Let G(x) = 1 + x*G(x)^7 be the g.f. of A002296, then the e.g.f. A(x) of this sequence satisfies:
(1) A'(x)/A(x) = G(x)^6 + 5*G'(x)/G(x).
(2) A(x) = F(x/A(x)^6) where F(x) is the e.g.f. of A251697.
(3) A(x) = Sum_{n>=0} A251697(n)*(x/A(x)^6)^n/n! where A251697(n) = (5*n+1) * (6*n+1)^(n-2) * 7^n .
(4) [x^n/n!] A(x)^(6*n+1) = (5*n+1) * (6*n+1)^(n-1) * 7^n .
a(n) = Sum_{k=0..n} 7^k * n!/k! * binomial(7*n-k-2,n-k) * (6*k-1)/(6*n-1) for n>=0.
Recurrence: 72*(2*n-1)*(3*n-2)*(3*n-1)*(6*n-5)*(6*n-1)*(588245*n^6 - 6117748*n^5 + 26651100*n^4 - 62321728*n^3 + 82554122*n^2 - 58646294*n + 17291583)*a(n) = 7*(69206436005*n^12 - 996572678472*n^11 + 6516703994430*n^10 - 25624338676965*n^9 + 67604945463195*n^8 - 126360374558838*n^7 + 171960790012102*n^6 - 171911061779835*n^5 + 125050872537045*n^4 - 63802357502870*n^3 + 20814954345360*n^2 - 3329274812661*n - 3763584000)*a(n-1) - 823543*(588245*n^6 - 2588278*n^5 + 4886035*n^4 - 5129908*n^3 + 3141733*n^2 - 958104*n - 720)*a(n-2). - Vaclav Kotesovec, Dec 07 2014
a(n) ~ 5 * 7^(7*n-3/2) / 6^(6*n-1/2) * n^(n-1) / exp(n-1). - Vaclav Kotesovec, Dec 07 2014

A186184 Expansion of 1/(1 - x*A002296(x)).

Original entry on oeis.org

1, 1, 2, 10, 89, 1002, 12592, 168805, 2363241, 34138860, 505042286, 7612594936, 116492572621, 1804984878387, 28260999959595, 446441276449715, 7106718529937710, 113886198966545724
Offset: 0

Views

Author

Vladimir Kruchinin, Feb 14 2011

Keywords

Links

Crossrefs

Cf. A002296.

Programs

  • Maple
    A186184 := proc(n)
   if n = 0 then
      1;
   else
      add( k/(6*n-5*k)*binomial(7*n-6*k-1,n-k), k=1..n) ;
   end if;
end proc:
seq(A186184(n),n=0..20) ; # R. J. Mathar, Feb 26 2011
  • Mathematica
    Join[{1},Table[Sum[k/(6n-5k) Binomial[7n-6k-1,n-k],{k,n}],{n,30}]] (* Harvey P. Dale, Aug 29 2012 *)

Formula

a(n) = Sum_{k=1..n} (k/(6*n-5*k))*binomial(7*n-6*k-1, n-k), n > 0.

A001764 a(n) = binomial(3*n,n)/(2*n+1) (enumerates ternary trees and also noncrossing trees).

Original entry on oeis.org

1, 1, 3, 12, 55, 273, 1428, 7752, 43263, 246675, 1430715, 8414640, 50067108, 300830572, 1822766520, 11124755664, 68328754959, 422030545335, 2619631042665, 16332922290300, 102240109897695, 642312451217745, 4048514844039120, 25594403741131680, 162250238001816900
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

Smallest number of straight line crossing-free spanning trees on n points in the plane.
Number of dissections of some convex polygon by nonintersecting diagonals into polygons with an odd number of sides and having a total number of 2n+1 edges (sides and diagonals). - Emeric Deutsch, Mar 06 2002
Number of lattice paths of n East steps and 2n North steps from (0,0) to (n,2n) and lying weakly below the line y=2x. - David Callan, Mar 14 2004
With interpolated zeros, this has g.f. 2*sqrt(3)*sin(arcsin(3*sqrt(3)*x/2)/3)/(3*x) and a(n) = C(n+floor(n/2),floor(n/2))*C(floor(n/2),n-floor(n/2))/(n+1). This is the first column of the inverse of the Riordan array (1-x^2,x(1-x^2)) (essentially reversion of y-y^3). - Paul Barry, Feb 02 2005
Number of 12312-avoiding matchings on [2n].
Number of complete ternary trees with n internal nodes, or 3n edges.
Number of rooted plane trees with 2n edges, where every vertex has even outdegree ("even trees").
a(n) is the number of noncrossing partitions of [2n] with all blocks of even size. E.g.: a(2)=3 counts 12-34, 14-23, 1234. - David Callan, Mar 30 2007
Pfaff-Fuss-Catalan sequence C^{m}_n for m=3, see the Graham et al. reference, p. 347. eq. 7.66.
Also 3-Raney sequence, see the Graham et al. reference, p. 346-7.
The number of lattice paths from (0,0) to (2n,0) using an Up-step=(1,1) and a Down-step=(0,-2) and staying above the x-axis. E.g., a(2) = 3; UUUUDD, UUUDUD, UUDUUD. - Charles Moore (chamoore(AT)howard.edu), Jan 09 2008
a(n) is (conjecturally) the number of permutations of [n+1] that avoid the patterns 4-2-3-1 and 4-2-5-1-3 and end with an ascent. For example, a(4)=55 counts all 60 permutations of [5] that end with an ascent except 42315, 52314, 52413, 53412, all of which contain a 4-2-3-1 pattern and 42513. - David Callan, Jul 22 2008
Central terms of pendular triangle A167763. - Philippe Deléham, Nov 12 2009
With B(x,t)=x+t*x^3, the comp. inverse in x about 0 is A(x,t) = Sum_{j>=0} a(j) (-t)^j x^(2j+1). Let U(x,t)=(x-A(x,t))/t. Then DU(x,t)/Dt=dU/dt+U*dU/dx=0 and U(x,0)=x^3, i.e., U is a solution of the inviscid Burgers's, or Hopf, equation. Also U(x,t)=U(x-t*U(x,t),0) and dB(x,t)/dt = U(B(x,t),t) = x^3 = U(x,0). The characteristics for the Hopf equation are x(t) = x(0) + t*U(x(t),t) = x(0) + t*U(x(0),0) = x(0) + t*x(0)^3 = B(x(0),t). These results apply to all the Fuss-Catalan sequences with 3 replaced by n>0 and 2 by n-1 (e.g., A000108 with n=2 and A002293 with n=4), see also A086810, which can be generalized to A133437, for associahedra. - Tom Copeland, Feb 15 2014
Number of intervals (i.e., ordered pairs (x,y) such that x<=y) in the Kreweras lattice (noncrossing partitions ordered by refinement) of size n, see the Bernardi & Bonichon (2009) and Kreweras (1972) references. - Noam Zeilberger, Jun 01 2016
Number of sum-indecomposable (4231,42513)-avoiding permutations. Conjecturally, number of sum-indecomposable (2431,45231)-avoiding permutations. - Alexander Burstein, Oct 19 2017
a(n) is the number of topologically distinct endstates for the game Planted Brussels Sprouts on n vertices, see Ji and Propp link. - Caleb Ji, May 14 2018
Number of complete quadrillages of 2n+2-gons. See Baryshnikov p. 12. See also Nov 10 2014 comments in A134264. - Tom Copeland, Jun 04 2018
a(n) is the number of 2-regular words on the alphabet [n] that avoid the patterns 231 and 221. Equivalently, this is the number of 2-regular tortoise-sortable words on the alphabet [n] (see the Defant and Kravitz link). - Colin Defant, Sep 26 2018
a(n) is the number of Motzkin paths of length 3n with n steps of each type, with the condition that (1, 0) and (1, 1) steps alternate (starting with (1, 0)). - Helmut Prodinger, Apr 08 2019
a(n) is the number of uniquely sorted permutations of length 2n+1 that avoid the patterns 312 and 1342. - Colin Defant, Jun 08 2019
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. - Tom Copeland, Dec 13 2019
The sequences of Fuss-Catalan numbers, of which this is the first after the Catalan numbers A000108 (the next is A002293), appear in articles on random matrices and quantum physics. See Banica et al., Collins et al., and Mlotkowski et al. Interpretations of these sequences in terms of the cardinality of specific sets of noncrossing partitions are provided by A134264. - Tom Copeland, Dec 21 2019
Call C(p, [alpha], g) the number of partitions of a cyclically ordered set with p elements, of cyclic type [alpha], and of genus g (the genus g Faa di Bruno coefficients of type [alpha]). This sequence counts the genus 0 partitions (non-crossing, or planar, partitions) of p = 3n into n parts of length 3: a(n) = C(3n, [3^n], 0). For genus 1 see A371250, for genus 2 see A371251. - Robert Coquereaux, Mar 16 2024
a(n) is the total number of down steps before the first up step in all 2_1-Dyck paths of length 3*n for n > 0. A 2_1-Dyck path is a lattice path with steps (1,2), (1,-1) that starts and ends at y = 0 and does not go below the line y = -1. - Sarah Selkirk, May 10 2020
a(n) is the number of pairs (A<=B) of noncrossing partitions of [n]. - Francesca Aicardi, May 28 2022
a(n) is the number of parking functions of size n avoiding the patterns 231 and 321. - Lara Pudwell, Apr 10 2023
Number of rooted polyominoes composed of n square cells of the hyperbolic regular tiling with Schläfli symbol {4,oo}. A rooted polyomino has one external edge identified, and chiral pairs are counted as two. A stereographic projection of the {4,oo} tiling on the Poincaré disk can be obtained via the Christensson link. - Robert A. Russell, Jan 27 2024
This is instance k = 3 of the family {C(k, n)}A130564.%20-%20_Wolfdieter%20Lang">{n>=0} given in a comment in A130564. - _Wolfdieter Lang, Feb 05 2024
The number of Apollonian networks (planar 3-trees) with n+3 vertices with a given base triangle. - Allan Bickle, Feb 20 2024
Number of rooted polyominoes composed of n tetrahedral cells of the hyperbolic regular tiling with Schläfli symbol {3,3,oo}. A rooted polyomino has one external face identified, and chiral pairs are counted as two. a(n) = T(n) in the second Beineke and Pippert link. - Robert A. Russell, Mar 20 2024

Examples

			a(2) = 3 because the only dissections with 5 edges are given by a square dissected by any of the two diagonals and the pentagon with no dissecting diagonal.
G.f. = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + 1428*x^6 + 7752*x^7 + 43263*x^8 + ...

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 23.
  • I. M. H. Etherington, On non-associative combinations, Proc. Royal Soc. Edinburgh, 59 (Part 2, 1938-39), 153-162.
  • I. M. H. Etherington, Some problems of non-associative combinations (I), Edinburgh Math. Notes, 32 (1940), pp. i-vi. Part II is by A. Erdelyi and I. M. H. Etherington, and is on pages vii-xiv of the same issue.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, pp. 200, 347. See also the Pólya-Szegő reference.
  • W. Kuich, Languages and the enumeration of planted plane trees. Nederl. Akad. Wetensch. Proc. Ser. A 73 = Indag. Math. 32, (1970), 268-280.
  • T. V. Narayana, Lattice Path Combinatorics with Statistical Applications. Univ. Toronto Press, 1979, p. 98.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis, Springer-Verlag, New York, Heidelberg, Berlin, 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).

Links

Crossrefs

Cf. A001762, A001763, A002294 - A002296, A006013, A025174, A063548, A064017, A072247, A072248, A134264, A143603, A258708, A256311, A188687 (binomial transform), A346628 (inverse binomial transform).
A column of triangle A102537.
Bisection of A047749 and A047761.
Row sums of triangles A108410 and A108767.
Second column of triangle A062993.
Mod 3 = A113047.
2D Polyominoes: A005034 (oriented), A005036 (unoriented), A369315 (chiral), A047749 (achiral), A000108 {3,oo}, A002293 {5,oo}.
3D Polyominoes: A007173 (oriented), A027610 (unoriented), A371350 (chiral), A371351 (achiral).
Cf. A130564 (for C(k, n) cases).

Programs

  • GAP
    List([0..25],n->Binomial(3*n,n)/(2*n+1)); # Muniru A Asiru, Oct 31 2018
  • Haskell
    a001764 n = a001764_list !! n
a001764_list = 1 : [a258708 (2 * n) n | n <- [1..]]
-- Reinhard Zumkeller, Jun 23 2015
  • Magma
    [Binomial(3*n,n)/(2*n+1): n in [0..30]]; // Vincenzo Librandi, Sep 04 2014
  • Maple
    A001764 := n->binomial(3*n,n)/(2*n+1): seq(A001764(n), n=0..25);
with(combstruct): BB:=[T,{T=Prod(Z,F),F=Sequence(B),B=Prod(F,Z,F)}, unlabeled]:seq(count(BB,size=i),i=0..22); # Zerinvary Lajos, Apr 22 2007
with(combstruct):BB:=[S, {B = Prod(S,S,Z), S = Sequence(B)}, labelled]: seq(count(BB, size=n)/n!, n=0..21); # Zerinvary Lajos, Apr 25 2008
n:=30:G:=series(RootOf(g = 1+x*g^3, g),x=0,n+1):seq(coeff(G,x,k),k=0..n); # Robert FERREOL, Apr 03 2015
alias(PS=ListTools:-PartialSums): A001764List := proc(m) local A, P, n;
A := [1,1]; P := [1]; for n from 1 to m - 2 do P := PS(PS([op(P), P[-1]]));
A := [op(A), P[-1]] od; A end: A001764List(25); # Peter Luschny, Mar 26 2022
  • Mathematica
    InverseSeries[Series[y-y^3, {y, 0, 24}], x] (* then a(n)=y(2n+1)=ways to place non-crossing diagonals in convex (2n+4)-gon so as to create only quadrilateral tiles *) (* Len Smiley, Apr 08 2000 *)
Table[Binomial[3n,n]/(2n+1),{n,0,25}] (* Harvey P. Dale, Jul 24 2011 *)
  • PARI
    {a(n) = if( n<0, 0, (3*n)! / n! / (2*n + 1)!)};
  • PARI
    {a(n) = if( n<0, 0, polcoeff( serreverse( x - x^3 + O(x^(2*n + 2))), 2*n + 1))};
  • PARI
    {a(n) = my(A); if( n<0, 0, A = 1 + O(x); for( m=1, n, A = 1 + x * A^3); polcoeff(A, n))};
  • PARI
    b=vector(22);b[1]=1;for(n=2,22,for(i=1,n-1,for(j=1,n-1,for(k=1,n-1,if((i-1)+(j-1)+(k-1)-(n-2),NULL,b[n]=b[n]+b[i]*b[j]*b[k])))));a(n)=b[n+1]; print1(a(0));for(n=1,21,print1(", ",a(n))) \\ Gerald McGarvey, Oct 08 2008
  • PARI
    Vec(1 + serreverse(x / (1+x)^3 + O(x^30))) \\ Gheorghe Coserea, Aug 05 2015
  • Python
    from math import comb
def A001764(n): return comb(3*n,n)//(2*n+1) # Chai Wah Wu, Nov 10 2022
  • Sage
    def A001764_list(n) :
    D = [0]*(n+1); D[1] = 1
    R = []; b = false; h = 1
    for i in range(2*n) :
        for k in (1..h) : D[k] += D[k-1]
        if not b : R.append(D[h])
        else : h += 1
        b = not b
    return R
A001764_list(22) # Peter Luschny, May 03 2012

Formula

From Karol A. Penson, Nov 08 2001: (Start)
G.f.: (2/sqrt(3*x))*sin((1/3)*arcsin(sqrt(27*x/4))).
E.g.f.: hypergeom([1/3, 2/3], [1, 3/2], 27/4*x).
Integral representation as n-th moment of a positive function on [0, 27/4]: a(n) = Integral_{x=0..27/4} (x^n*((1/12) * 3^(1/2) * 2^(1/3) * (2^(1/3)*(27 + 3 * sqrt(81 - 12*x))^(2/3) - 6 * x^(1/3))/(Pi * x^(2/3)*(27 + 3 * sqrt(81 - 12*x))^(1/3)))), n >= 0. This representation is unique. (End)
G.f. A(x) satisfies A(x) = 1+x*A(x)^3 = 1/(1-x*A(x)^2) [Cyvin (1998)]. - Ralf Stephan, Jun 30 2003
a(n) = n-th coefficient in expansion of power series P(n), where P(0) = 1, P(k+1) = 1/(1 - x*P(k)^2).
G.f. Rev(x/c(x))/x, where c(x) is the g.f. of A000108 (Rev=reversion of). - Paul Barry, Mar 26 2010
From Gary W. Adamson, Jul 07 2011: (Start)
Let M = the production matrix:
1, 1
2, 2, 1
3, 3, 2, 1
4, 4, 3, 2, 1
5, 5, 4, 3, 2, 1
...
a(n) = upper left term in M^n. Top row terms of M^n = (n+1)-th row of triangle A143603, with top row sums generating A006013: (1, 2, 7, 30, 143, 728, ...). (End)
Recurrence: a(0)=1; a(n) = Sum_{i=0..n-1, j=0..n-1-i} a(i)a(j)a(n-1-i-j) for n >= 1 (counts ternary trees by subtrees of the root). - David Callan, Nov 21 2011
G.f.: 1 + 6*x/(Q(0) - 6*x); Q(k) = 3*x*(3*k + 1)*(3*k + 2) + 2*(2*(k^2) + 5*k +3) - 6*x*(2*(k^2) + 5*k + 3)*(3*k + 4)*(3*k + 5)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 27 2011
D-finite with recurrence: 2*n*(2n+1)*a(n) - 3*(3n-1)*(3n-2)*a(n-1) = 0. - R. J. Mathar, Dec 14 2011
REVERT transform of A115140. BINOMIAL transform is A188687. SUMADJ transform of A188678. HANKEL transform is A051255. INVERT transform of A023053. INVERT transform is A098746. - Michael Somos, Apr 07 2012
(n + 1) * a(n) = A174687(n).
G.f.: F([2/3,4/3], [3/2], 27/4*x) / F([2/3,1/3], [1/2], (27/4)*x) where F() is the hypergeometric function. - Joerg Arndt, Sep 01 2012
a(n) = binomial(3*n+1, n)/(3*n+1) = A062993(n+1,1). - Robert FERREOL, Apr 03 2015
a(n) = A258708(2*n,n) for n > 0. - Reinhard Zumkeller, Jun 23 2015
0 = a(n)*(-3188646*a(n+2) + 20312856*a(n+3) - 11379609*a(n+4) + 1437501*a(n+5)) + a(n+1)*(177147*a(n+2) - 2247831*a(n+3) + 1638648*a(n+4) - 238604*a(n+5)) + a(n+2)*(243*a(n+2) + 31497*a(n+3) - 43732*a(n+4) + 8288*a(n+5)) for all integer n. - Michael Somos, Jun 03 2016
a(n) ~ 3^(3*n + 1/2)/(sqrt(Pi)*4^(n+1)*n^(3/2)). - Ilya Gutkovskiy, Nov 21 2016
Given g.f. A(x), then A(1/8) = -1 + sqrt(5), A(2/27) = (-1 + sqrt(3))*3/2, A(4/27) = 3/2, A(3/64) = -2 + 2*sqrt(7/3), A(5/64) = (-1 + sqrt(5))*2/sqrt(5), etc. A(n^2/(n+1)^3) = (n+1)/n if n > 1. - Michael Somos, Jul 17 2018
From Peter Bala, Sep 14 2021: (Start)
A(x) = exp( Sum_{n >= 1} (1/3)*binomial(3*n,n)*x^n/n ).
The sequence defined by b(n) := [x^n] A(x)^n = A224274(n) for n >= 1 and satisfies the congruence b(p) == b(1) (mod p^3) for prime p >= 3. Cf. A060941. (End)
G.f.: 1/sqrt(B(x)+(1-6*x)/(9*B(x))+1/3), with B(x):=((27*x^2-18*x+2)/54-(x*sqrt((-(4-27*x))*x))/(2*3^(3/2)))^(1/3). - Vladimir Kruchinin, Sep 28 2021
x*A'(x)/A(x) = (A(x) - 1)/(- 2*A(x) + 3) = x + 5*x^2 + 28*x^3 + 165*x^4 + ... is the o.g.f. of A025174. Cf. A002293 - A002296. - Peter Bala, Feb 04 2022
a(n) = hypergeom([1 - n, -2*n], [2], 1). Row sums of A108767. - Peter Bala, Aug 30 2023
G.f.: z*exp(3*z*hypergeom([1, 1, 4/3, 5/3], [3/2, 2, 2], (27*z)/4)) + 1.
- Karol A. Penson, Dec 19 2023
G.f.: hypergeometric([1/3, 2/3], [3/2], (3^3/2^2)*x). See the e.g.f. above. - Wolfdieter Lang, Feb 04 2024
a(n) = (3*n)! / (n!*(2*n+1)!). - Allan Bickle, Feb 20 2024
Sum_{n >= 0} a(n)*x^n/(1 + x)^(3*n+1) = 1. See A316371 and A346627. - Peter Bala, Jun 02 2024
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^5). - Seiichi Manyama, Jun 16 2025

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

N. J. A. Sloane

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).

Links

Crossrefs

Column k=3 of triangle A062993 and A070914.
Cf. A000260, A002295, A002296, A027836, A062994, A346646 (binomial transform), A346664 (inverse binomial transform).
Cf. A006632, A006633, A006634, A025174, A069271, A196678, A224274, A233658, A233666, A233667, A277877, A283049, A283101, A283102, A283103.
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

N. J. A. Sloane

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).

Links

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

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

Original entry on oeis.org

1, 1, 6, 51, 506, 5481, 62832, 749398, 9203634, 115607310, 1478314266, 19180049928, 251857119696, 3340843549855, 44700485049720, 602574657427116, 8175951659117794, 111572030260242090, 1530312970340384580, 21085148778264281865, 291705220704719165526
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

From Wolfdieter Lang, Sep 14 2007: (Start)
a(n), n >= 1, enumerates sextic (6-ary) trees (rooted, ordered, incomplete) with n vertices (including the root).
Pfaff-Fuss-Catalan sequence C^{m}_n for m=6. See the Graham et al. reference, p. 347. eq. 7.66. See also the Pólya-Szegő reference.
Also 6-Raney sequence. See the Graham et al. reference, p. 346-7. (End)
This is instance k = 6 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)=6 sextic trees (vertex degree <= 6 and 6 possible branchings) with 2 vertices (one of them the root). Adding one more branch (one more vertex) to these 6 trees yields 6*6 + binomial(6,2) = 51 = 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.
  • 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).
  • 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 Lamé, Rodrigues, Binet, Catalan und Duhamel in dem Journal de Mathématiques pures et appliquées, publié par Joseph Liouville. T. III. IV.", Archiv der Mathematik u. Physik, 1 (1841), pp. 193ff; see especially p. 198.

Links

Crossrefs

Cf. A001764, A002293, A002294, A002296, A130564.
Fifth column of triangle A062993.

Programs

Formula

O.g.f.: A(x) = 1 + x*A(x)^6 = 1/(1-x*A(x)^5).
a(n) = binomial(6*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
5, 5, 1
15, 15, 5, 1
35, 35, 15, 5, 1
...
where (1, 5, 15, 35, ...) = A000332 starting with 1. - Gary W. Adamson, Jul 08 2011
a(n) are special values of Jacobi polynomials, in Maple notation:
a(n) = JacobiP(n-1, 5*n+1, -n, 1)/n, n=1, 2, ... . - Karol A. Penson, Mar 17 2015
a(n) = binomial(6*n+1, n)/(6*n+1) = A062993(n+4,4). - Robert FERREOL, Apr 03 2015
a(0) = 1; a(n) = Sum_{i1+i2+...+i6=n-1} a(i1)*a(i2)*...*a(i6) for n>=1. - Robert FERREOL, Apr 03 2015
D-finite with recurrence: 5*n*(5*n+1)*(5*n-3)*(5*n-2)*(5*n-1)*a(n) - 72*(6*n-5)*(6*n-1)*(3*n-1)*(2*n-1)*(3*n-2)*a(n-1) = 0. - R. J. Mathar, Sep 06 2016
From Ilya Gutkovskiy, Jan 15 2017: (Start)
O.g.f.: 5F4(1/6,1/3,1/2,2/3,5/6; 2/5,3/5,4/5,6/5; 46656*x/3125).
E.g.f.: 5F5(1/6,1/3,1/2,2/3,5/6; 2/5,3/5,4/5,1,6/5; 46656*x/3125).
a(n) ~ 3^(6*n+1/2)*64^n/(sqrt(Pi)*5^(5*n+3/2)*n^(3/2)). (End)
x*A'(x)/A(x) = (A(x) - 1)/(- 5*A(x) + 6) = x + 11*x^2 + 136*x^3 + 1771*x^4 + ... = (1/6)*Sum_{n >= 1} binomial(6*n,n)*x^n. Cf. A001764 and A002293 - A002296. - Peter Bala, Feb 04 2022
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^11). - Seiichi Manyama, Jun 16 2025

Extensions

More terms from Stefan Steinerberger, Apr 06 2006
Edited by M. F. Hasler, Apr 08 2015

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.
The ordinary generating functions for the k-family {c(k, n+1)}{n>=0}, k >= 1, are G(k, x) = hypergeometric(Aseq(k+1), Bseq(k), ((k+1)^(k+1)/k^k)*x), with Aseq(k+1) = [a(k)_1,..., a(k){k+1}], where a(k)j = (2*k - (j-1))/(k+1), and Bseq(k) = [b(k)_1, ..., b(k)_k], where b(k)_j = (2*k - (j-1))/k. The e.g.f. has an extra 1 in the B-section, which leads to a cancellation with the A-section 1 term. Thanks to Dixon J. Jones for asking for the general formulas. - _Wolfdieter Lang, Feb 04 2024
The o.g.f. of {C(k, n)}{n>=0} is in the Graham-Knuth-Patashnik book denoted as B_k(z) on pp. 200, 349 (2nd ed. 1994, pp. 200, 363). - _Wolfdieter Lang, Mar 07 2024

References

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

Links

Crossrefs

Cf. A000012, A000108, A001764, A002293, A002294, A002295, A002296, A006013, A062994, A006632, A007556, A118971, A130565, A234466, A234513, A234573, A235340.

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

A059968 Number of 10-ary trees.

Original entry on oeis.org

1, 1, 10, 145, 2470, 46060, 910252, 18730855, 397089550, 8612835715, 190223180840, 4263421511271, 96723482198980, 2216905597676000, 51256802757808320, 1194060413809070710, 27999654303202465310, 660370070571422998410, 15654733143626084944150
Offset: 0

Views

Author

Claude Lenormand (claude.lenormand(AT)free.fr), Mar 05 2001

Keywords

Comments

From Wolfdieter Lang, Feb 06 2020: (Start)
Ninth column of triangle A062993 (without leading zeros). A Pfaff-Fuss or 10-Raney sequence.
a(n), n>=1, enumerates 10-ary trees (rooted, ordered, incomplete) with n vertices (including the root).
See Graham et al., Hilton and Pedersen, Hoggat and Bicknell, Frey and Sellers references given in A062993. (End)
This is instance k = 10 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)=10 10-ary trees (vertex degree <=10 and 10 possible branchings) with 2 vertices (one of them the root). Adding one more branch (one more vertex) to these 10 trees yields 10*10+binomial(10,2)=145=a(3) such trees. - _Wolfdieter Lang_, Sep 14 2007.

References

  • 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.

Links

Crossrefs

Related algebraic sequences concerning trees: strictly k-ary trees (A000108: s=x+s^2, A001263: s=(x, y)+(x, s)+(s, y)+(s, s))), (A001764: s=x+s^3), (A002293: s=x+s^4), (A002294: s=x+s^5), (A002295: s=x+s^6), (A002296: s=x+s^7), (A007556: s=x+s^8), at most k-ary trees (A001006: s=x+xs+xs^2), (A036765-A036769, s=x+xs^2....+xs^k, k=3, 4, 5, 6, 7).
Pfaff-Fuss A062993 column members: A000108, A001764, A002293, A002294, A002295, A002296, A007556, A062994, A230388.
Cf. A130564.

Programs

  • Maple
    seq(binomial(10*k+1, k)/(9*k+1), k=0..30);
n:=30:G:=series(RootOf(g = 1+x*g^10, g), x=0, n+1):seq(coeff(G, x, k), k=0..n); # Robert FERREOL, Apr 01 2015
  • Mathematica
    a[n_] := Binomial[10n, n]/(9n+1);
a /@ Range[0, 25] (* Jean-François Alcover, Jan 17 2020 *)

Formula

G.f. A(x) satisfies: A = x + A^10.
a(n) = binomial(k*n, n)/((k-1)*n+1), for k=10.
Recurrence: a(0) = 1; a(n) = Sum_{i1+i2+..i10=n-1} a(i1)*a(i2)*...*a(i10) for n>=1. - Robert FERREOL, Apr 01 2015
From Wolfdieter Lang, Feb 06 2020: (Start)
a(n) = A062993(n+8, 8). [Corrected by Robert FERREOL, Apr 01 2015]
G.f.: RootOf((_Z^10)*x-_Z+1) (Maple notation, from ECS, see links for A007556).
G.f.: hypergeometric([1, 2, 3, 4, 5, 6, 7, 8, 9]/10, [2, 3, 4, 5, 6, 7, 8, 10]/9, (10^10/9^9)*x),
E.g.f.: hypergeometric([1, 2, 3, 4, 5, 6, 7, 8, 9]/10, [2, 3, 4, 5, 6, 7, 8, 9, 10]/9, (10^10/9^9)*x).
For other family members see the crossreferences.
(End)
D-finite with recurrence 81*n*(9*n-7)*(9*n-5)*(3*n-1)*(9*n-1)*(9*n+1)*(3*n-2)*(9*n-4)*(9*n-2)*a(n) -800*(10*n-9)*(5*n-4)*(10*n-7)*(5*n-3)*(2*n-1)*(5*n-2)*(10*n-3)*(5*n-1)*(10*n-1)*a(n-1)=0. - R. J. Mathar, Mar 21 2022
a(n) ~ (10^10/9^9)^n*sqrt(10/(2*Pi*(9*n)^3)). - Robert A. Russell, Jul 15 2024
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^19). - Seiichi Manyama, Jun 16 2025

Extensions

More terms from James Sellers, Mar 15 2001
a(0)=1 inserted by Alois P. Heinz, Jan 17 2020
A062744 merged into this sequence by Wolfdieter Lang, Feb 06 2020

A062994 Eighth column of triangle A062993 (without leading zeros). A Pfaff-Fuss or 9-Raney sequence.

Original entry on oeis.org

1, 1, 9, 117, 1785, 29799, 527085, 9706503, 184138713, 3573805950, 70625252863, 1416298046436, 28748759731965, 589546754316126, 12195537924351375, 254184908607118800, 5332692942907262361
Offset: 0

Views

Author

Wolfdieter Lang, Jul 12 2001

Keywords

Comments

See Graham et al., Hilton and Pedersen, Hoggat and Bicknell, Frey and Sellers references given in A062993.
Essentially the same as A059967. a(n), n>=1, enumerates 9-ary trees (rooted, ordered, incomplete) with n vertices (including the root).
These numbers appear in a formula on p. 24 of Gross et al. for b = -2 or 4. For b = -1 or 3, see A002293.- Tom Copeland, Dec 24 2019
This is instance k = 9 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)=9 9-ary trees (vertex degree <=9 and 9 possible branchings) with 2 vertices (one of them the root). Adding one more branch (one more vertex) to these 9 trees yields 9*9 + binomial(9,2) = 117 = a(3) such trees.

References

  • 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.

Links

Crossrefs

Column k=8 of A070914.
Cf. A000108, A001764, A002293, A002294, A002295, A002296, A007556, A059968, A130564.

Programs

  • Maple
    seq(binomial(9*k+1,k)/(8*k+1),k=0..30);
n:=30: G:=series(RootOf(g = 1+x*g^9, g),x=0,n+1): seq(coeff(G,x,k),k=0..n); # Robert FERREOL, Apr 01 2015
  • Mathematica
    Table[Binomial[9n,n]/(8n+1),{n,0,30}] (* Harvey P. Dale, Oct 28 2012 *)
  • PARI
    { for (n=0, 100, write("b062994.txt", n, " ", binomial(9*n, n)/(8*n + 1)) ) } \\ Harry J. Smith, Aug 15 2009

Formula

a(n) = A062993(n+7, 7) = binomial(9*n, n)/(8*n+1).
G.f.: RootOf((_Z^9)*x-_Z+1) (Maple notation, from ECS, see links for A007556).
Recurrence: a(0) = 1; a(n) = Sum_{i1+i2+..+i9=n-1} a(i1)*a(i2)*...*a(i9) for n>=1. - Robert FERREOL, Apr 01 2015
From Ilya Gutkovskiy, Jan 16 2017: (Start)
O.g.f.: 8F7(1/9,2/9,1/3,4/9,5/9,2/3,7/9,8/9; 1/4,3/8,1/2,5/8,3/4,7/8,9/8; 387420489*x/16777216).
E.g.f.: 8F8(1/9,2/9,1/3,4/9,5/9,2/3,7/9,8/9; 1/4,3/8,1/2,5/8,3/4,7/8,1,9/8; 387420489*x/16777216).
a(n) ~ 3^(18*n+1)/(sqrt(Pi)*2^(24*n+5)*n^(3/2)). (End)
D-finite with recurrence: 128*n*(8*n-5)*(4*n-1)*(8*n+1)*(2*n-1)*(8*n-1)*(4*n-3)*(8*n-3)*a(n) -81*(9*n-7)*(9*n-5)*(3*n-1)*(9*n-1)*(9*n-8)*(3*n-2)*(9*n-4)*(9*n-2)*a(n-1)=0. - R. J. Mathar, Feb 20 2020
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^17). - Seiichi Manyama, Jun 16 2025

Extensions

9-ary tree comments and Pólya and G. Szegő reference from Wolfdieter Lang, Sep 14 2007
Showing 1-10 of 50 results. Next

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.