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

A000699 Number of irreducible chord diagrams with 2n nodes.

Original entry on oeis.org

1, 1, 1, 4, 27, 248, 2830, 38232, 593859, 10401712, 202601898, 4342263000, 101551822350, 2573779506192, 70282204726396, 2057490936366320, 64291032462761955, 2136017303903513184, 75197869250518812754, 2796475872605709079512, 109549714522464120960474, 4509302910783496963256400, 194584224274515194731540740
Offset: 0

Views

Author

Keywords

Comments

Perturbation expansion in quantum field theory: spinor case in 4 spacetime dimensions.
a(n)*2^(-n) is the coefficient of the x^(2*n-1) term in the series reversal of the asymptotic expansion of 2*DawsonF(x) = sqrt(Pi)*exp(-x^2)*erfi(x) for x -> inf. - Vladimir Reshetnikov, Apr 23 2016
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018
A set partition is topologically connected if the graph whose vertices are the blocks and whose edges are crossing pairs of blocks is connected, where two blocks cross each other if they are of the form {{...x...y...},{...z...t...}} for some x < z < y < t or z < x < t < y. Then a(n) is the number of topologically connected 2-uniform set partitions of {1...2n}. See my links for examples. - Gus Wiseman, Feb 23 2019
From Julien Courtiel, Oct 09 2024: (Start)
a(n) is the number of rooted bridgeless combinatorial maps with n edges (genus is not fixed). A map is bridgeless if it has no edge whose removal disconnects the graph. For example, for n=2, there are 4 bridgeless maps with 2 edges: 2 planar maps with 1 vertex (either two consecutive loops, or two nested loops), 1 toric map with 1 vertex, and 1 planar map with 2 vertices connected by a double edge.
Also, a(n) is the number of trees with n edges equipped with a binary tubing. A tube is a connected subgraph. A binary tubing of a tree is a nested set collection S of tubes such that 1. S contains the tube of all vertices 2. Every tube of S is either reduced to one vertex, or it can be can partitioned by 2 tubes of S.
(End)

Examples

			a(31)=627625976637472254550352492162870816129760 was computed using Kreimer's Hopf algebra of rooted trees. It subsumes 2.6*10^21 terms in quantum field theory.
G.f.: A(x) = 1 + x + x^2 + 4*x^3 + 27*x^4 + 248*x^5 + 2830*x^6 +...
where d/dx (A(x) - 1)^2/x = 1 + 4*x + 27*x^2 + 248*x^3 + 2830*x^4 +...
		

References

  • 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

Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.
Cf. A004300, A051862, A212273. Column sums of A232223. First column of A322402.

Programs

  • Maple
    A000699 := proc(n)
        option remember;
        if n <= 1 then
            1;
        else
            add((2*i-1)*procname(i)*procname(n-i),i=1..n-1) ;
        end if;
    end proc:
    seq(A000699(n),n=0..30) ; # R. J. Mathar, Jun 12 2018
  • Mathematica
    terms = 22; A[] = 0; Do[A[x] = x + x^2 * D[A[x]^2/x, x] + O[x]^(terms+1) // Normal, terms]; CoefficientList[A[x], x] // Rest (* Jean-François Alcover, Apr 06 2012, after Paul D. Hanna, updated Jan 11 2018 *)
    a = ConstantArray[0,20]; a[[1]]=1; Do[a[[n]] = (n-1)*Sum[a[[i]]*a[[n-i]],{i,1,n-1}],{n,2,20}]; a (* Vaclav Kotesovec, Feb 22 2014 *)
    Module[{max = 20, s}, s = InverseSeries[ComplexExpand[Re[Series[2 DawsonF[x], {x, Infinity, 2 max + 1}]]]]; Table[SeriesCoefficient[s, 2 n - 1] 2^n, {n, 1, max}]] (* Vladimir Reshetnikov, Apr 23 2016 *)
  • PARI
    {a(n)=local(A=1+x*O(x^n)); for(i=1, n, A=1+x+x^2*deriv((A-1)^2/x)+x*O(x^n)); polcoeff(A, n)} \\ Paul D. Hanna, Dec 31 2010 [Modified to include a(0) = 1. - Paul D. Hanna, Nov 06 2020]
    
  • PARI
    {a(n) = my(A); A = 1+O(x) ; for( i=0, n, A = 1+x + (A-1)*(2*x*A' - A + 1)); polcoeff(A, n)}; /* Michael Somos, May 12 2012 [Modified to include a(0) = 1. - Paul D. Hanna, Nov 06 2020] */
    
  • PARI
    seq(N) = {
      my(a = vector(N)); a[1] = 1;
      for (n=2, N, a[n] = sum(k=1, n-1, (2*k-1)*a[k]*a[n-k])); a;
    };
    seq(22)  \\ Gheorghe Coserea, Jan 22 2017
    
  • PARI
    seq(n)={my(g=serlaplace(1 / sqrt(1 - 2*x + O(x*x^n)))); Vec(sqrt((x/serreverse( x*g^2 ))))} \\ Andrew Howroyd, Nov 21 2024
    
  • Python
    def A000699_list(n):
        list = [1, 1] + [0] * (n - 1)
        for i in range(2, n + 1):
            list[i] = (i - 1) * sum(list[j] * list[i - j] for j in range(1, i))
        return list
    print(A000699_list(22)) # M. Eren Kesim, Jun 23 2021

Formula

a(n) = (n-1)*Sum_{i=1..n-1} a(i)*a(n-i) for n > 1, with a(1) = a(0) = 1. [Modified to include a(0) = 1. - Paul D. Hanna, Nov 06 2020]
A212273(n) = n * a(n). - Michael Somos, May 12 2012
G.f. satisfies: A(x) = 1 + x + x^2*[d/dx (A(x) - 1)^2/x]. - Paul D. Hanna, Dec 31 2010 [Modified to include a(0) = 1. - Paul D. Hanna, Nov 06 2020]
a(n) ~ n^n * 2^(n+1/2) / exp(n+1) * (1 - 31/(24*n) - 2207/(1152*n^2) - 3085547/(414720*n^3) - 1842851707/(39813120*n^4) - ...). - Vaclav Kotesovec, Feb 22 2014, extended Oct 23 2017
G.f. A(x) satisfies: 1 = A(x) - x/(A(x) - 2*x/(A(x) - 3*x/(A(x) - 4*x/(A(x) - 5*x/(A(x) - ...))))), a continued fraction relation. - Paul D. Hanna, Nov 04 2020
G.f. A(x) satisfies: A(x*B(x)^2) = B(x) where B(x) is the g.f. of A001147. - Andrew Howroyd, Nov 21 2024

Extensions

More terms from David Broadhurst, Dec 14 1999
Inserted "chord" in definition. - N. J. A. Sloane, Jan 19 2017
Added a(0)=1. - N. J. A. Sloane, Nov 05 2020
Modified formulas slightly to include a(0) = 1. - Paul D. Hanna, Nov 06 2020

A000698 A problem of configurations: a(0) = 1; for n>0, a(n) = (2n-1)!! - Sum_{k=1..n-1} (2k-1)!! a(n-k). Also the number of shellings of an n-cube, divided by 2^n n!.

Original entry on oeis.org

1, 1, 2, 10, 74, 706, 8162, 110410, 1708394, 29752066, 576037442, 12277827850, 285764591114, 7213364729026, 196316804255522, 5731249477826890, 178676789473121834, 5925085744543837186, 208256802758892355202, 7734158085942678174730
Offset: 0

Views

Author

Keywords

Comments

Also number of nonisomorphic unlabeled connected Feynman diagrams of order 2n-2 for the electron propagator of quantum electrodynamics (QED), including vanishing diagrams. [Corrected by Charles R Greathouse IV, Jan 24 2014][Clarified by Robert Coquereaux, Sep 14 2014]
a(n+1) is the moment of order 2*n for the probability density function rho(x) = (1/sqrt(2*Pi))*exp(x^2/2)/[(u(x))^2+Pi/2], with u(x) = Integral_{t=0..x} exp(t*t/2) dt, on the real interval -infinity..infinity. - Groux Roland, Jan 13 2009
Starting (1, 2, 10, 74, ...) = INVERTi transform of A001147: (1, 3, 15, 105, ...). - Gary W. Adamson, Oct 21 2009
The Cvitanovic et al. paper relates this sequence to A005411 and A005413. - Robert Munafo, Jan 24 2010
Hankel transform of a(n+1) is A168467. - Paul Barry, Nov 26 2009
a(n) = number of labeled Dyck (n-1)-paths (A000108) in which each vertex that terminates an upstep is labeled with an integer i in [0,h], where h is the height of the vertex . For example UDUD contributes 4 labeled paths--0D0D, 0D1D, 1D0D, 1D1D where upsteps are replaced by their labels--and UUDD contributes 6 labeled paths to a(3)=10. The Deléham (Mar 24 2007) formula below counts these labeled paths by number of "0" labels. - David Callan, Aug 23 2011
a(n) is the number of indecomposable perfect matchings on [2n]. A perfect matching on [2n] is decomposable if a nonempty subset of the edges forms a perfect matching on [2k] for some kDavid Callan, Nov 29 2012
From Robert Coquereaux, Sep 12 2014: (Start)
QED diagrams are graphs with two kinds of edges (lines): a (non-oriented), f (oriented), and only one kind of (internal) vertex: aff. They may have internal and external (i.e., pendant) lines. The order is the number of (internal) vertices. Vanishing diagrams: QED diagrams containing loops of type f with an odd number of vertices are set to 0 (Furry theorem). Proper diagrams: connected QED diagrams that remain connected when an arbitrary internal line is cut.
The number of Feynman diagrams of order 2n for the electron propagator (2-point function of QED), vanishing or not, proper or not, of order 2n, starting from n = 0, is given by 1, 2, 10, 74, 706, 8162, ..., i.e., this sequence A000698, with the first term (equal to 1) dropped. Call Sf the associated g.f.
The number of non-vanishing Feynman diagrams, for the same 2-point function, is given by 1, 1, 4, 25, 208, 2146, ..., i.e., by the sequence A005411, with a first term of order 0, equal to 1, added. Call S the associated g.f.
If one does not remove the vanishing diagram, but, at the same time, considers only those graphs that are proper, one obtains the Feynman diagrams (vanishing and non-vanishing) for the self-energy function of QED, 0, 1, 3, 21, 207, 2529, ..., i.e., the sequence A115974 with a first term of order 0, equal to 0, added. A115974 is twice A167872. Call Sigmaf the associated g.f.
If one removes the vanishing diagrams and, at the same time, considers only those graphs that are proper, one obtains the Feynman diagrams for the self-energy function of QED given by 0, 1, 3, 18, 153, 1638, ..., i.e., by the sequence A005412, with a first term of order 0, equal to 0, added. Call Sigma the associated g.f.
Then Sf = 1/(1-Sigmaf) and S = 1/(1-Sigma). (End)
For n>0 sum over all Dyck paths of semilength n-1 of products over all peaks p of (x_p+y_p)/y_p, where x_p and y_p are the coordinates of peak p. - Alois P. Heinz, May 22 2015
Also, counts certain isomorphism classes of closed normal linear lambda terms. [N. Zeilberger, 2015]. - N. J. A. Sloane, Sep 18 2016
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018
For n >= 2, a(n) is the number of coalescent histories for a pair consisting of a matching lodgepole gene tree and species tree with 2n-1 leaves. - Noah A Rosenberg, Jun 21 2022

Examples

			G.f. = 1 + x + 2*x^2 + 10*x^3 + 74*x^4 + 706*x^5 + 8162*x^6 + 110410*x^7 + ...
		

References

  • Dubois C., Giorgetti A., Genestier R. (2016) Tests and Proofs for Enumerative Combinatorics. In: Aichernig B., Furia C. (eds) Tests and Proofs. TAP 2016. Lecture Notes in Computer Science, vol 9762. Springer.
  • R. W. Robinson, Counting irreducible Feynman diagrams exactly and asymptotically, Abstracts Amer. Math. Soc., 2002, #975-05-270.
  • 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

Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.
Column k=1 of A258219, A258222.
Row sums of A322398.

Programs

  • Maple
    A006882 := proc(n) option remember; if n <= 1 then 1 else n*procname(n-2); fi; end;
    A000698:=proc(n) option remember; global df; local k; if n=0 then RETURN(1); fi; A006882(2*n-1) - add(A006882(2*k-1)*A000698(n-k),k=1..n-1); end;
    A000698 := proc(n::integer) local resul,fac,pows,c,c1,p,i ; if n = 0 then RETURN(1) ; else pows := combinat[partition](n) ; resul := 0 ; for p from 1 to nops(pows) do c := combinat[permute](op(p,pows)) ; c1 := op(1,c) ; fac := nops(c) ; for i from 1 to nops(c1) do fac := fac*doublefactorial(2*op(i,c1)-1) ; od ; resul := resul-(-1)^nops(c1)*fac ; od : fi ; RETURN(resul) ; end; # R. J. Mathar, Apr 24 2006
    # alternative Maple program:
    b:= proc(x, y, t) option remember; `if`(y>x or y<0, 0,
          `if`(x=0, 1, b(x-1, y-1, false)*`if`(t, (x+y)/y, 1) +
                       b(x-1, y+1, true)  ))
        end:
    a:= n-> `if`(n=0, 1, b(2*n-2, 0, false)):
    seq(a(n), n=0..25);  # Alois P. Heinz, May 23 2015
    a_list := proc(len) local n, A; if len=1 then return [1] fi: A := Array(-1..len-2); A[-1] := 1; A[0] := 1; for n to len-2 do A[n] := (2*n-1)*A[n-1]+add(A[j]*A[n-j-1], j=0..n-1) od: convert(A, list) end: a_list(20); # Peter Luschny, Jul 18 2017
  • Mathematica
    a[n_] := a[n] = (2n - 1)!! - Sum[ a[n - k](2k - 1)!!, {k, n-1}]; Array[a, 18, 0] (* Ignacio D. Peixoto, Jun 23 2006 *)
    a[ n_] := If[ n < 0, 0, SeriesCoefficient[ 2 - 1 / Sum[ (2 k - 1)!! x^k, {k, 0, n}], {x, 0, n}]]; (* Michael Somos, Nov 16 2011 *)
    a[n_]:= SeriesCoefficient[1+x(1/x+(E^((1/2)/x) Sqrt[2/\[Pi]] Sqrt[-(1/x)])/Erfc[Sqrt[-(1/x)]/Sqrt[2]]), {x,0,n}, Assumptions -> x >0](* Robert Coquereaux, Sep 14 2014 *)
    max = 20; g = t/Fold[1 - ((t + #2)*z)/#1 &, 1, Range[max, 1, -1]]; T[n_, k_] := SeriesCoefficient[g, {z, 0, n}, {t, 0, k}]; a[0] = 1; a[n_] := Sum[T[n-1, k], {k, 0, n}]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jan 31 2016, after Philippe Deléham *)
  • PARI
    {a(n) = if( n<0, 0, polcoeff( 2 - 1 / sum( k=0, n, x^k * (2*k)! /(2^k * k!), x * O(x^n)), n))}; /* Michael Somos, Feb 08 2011 */
    
  • PARI
    {a(n) = my(A); if( n<1, n==0, A = vector(n); A[1] = 1; for( k=2, n, A[k] = (2*k - 3) * A[k-1] + sum( j=1, k-1, A[j] * A[k-j])); A[n])}; /* Michael Somos, Jul 24 2011 */
    
  • Python
    from sympy import factorial2, cacheit
    @cacheit
    def a(n): return 1 if n == 0 else factorial2(2*n - 1) - sum(factorial2(2*k - 1)*a(n - k) for k in range(1, n))
    [a(n) for n in range(51)]  # Indranil Ghosh, Jul 18 2017

Formula

G.f.: 2 - 1/(1 + Sum_{n>=1} (2*n-1)!! * x^n ).
a(n+1) = Sum_{k=0..n} A089949(n, k)*2^k. - Philippe Deléham, Aug 15 2005
a(n+1) = Sum_{k=0..n} A053979(n,k). - Philippe Deléham, Mar 24 2007
From Paul Barry, Nov 26 2009: (Start)
G.f.: 1+x/(1-2x/(1-3x/(1-4x/(1-5x/(1-6x/(1-... (continued fraction).
G.f.: 1+x/(1-2x-6x^2/(1-7x-20x^2/(1-11x-42x^2/(1-15x-72x^2/(1-19x-110x^2/(1-... (continued fraction). (End)
G.f.: 1 + x * B(x) * C(x) where B(x) is the g.f. for A001147 and C(x) is the g.f. for A005416. - Michael Somos, Feb 08 2011
G.f.: 1+x/W(0); where W(k)=1+x+x*2k-x*(2k+3)/W(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 17 2011
From Peter Bala, Dec 22 2011: (Start)
Recurrence relation: a(n+1) = (2*n-1)*a(n) + Sum_{k = 1..n} a(k)*a(n+1-k) for n >= 0 and a(1) = 1.
The o.g.f. B(x) = Sum_{n>=1} a(n)*x^(2*n-1) = x + 2*x^3 + 10*x^5 + 74*x^7 + ... satisfies the Riccati differential equation y'(x) = -1/x^2 + (1/x^3)*y(x) - (1/x^2)*y(x)^2 with initial condition y(0) = 0 (cf. A005412). The solution is B(x) = 1/z(x) + 1/x, where z(x) = -Sum_{n>=0} A001147(n) * x^(2*n+1) = -(x + x^3 + 3*x^5 + 15*x^7 + ...). The function b(x) = -B(1/x) satisfies b'(x) = -1 - (x + b(x))*b(x). Hence the differential operator (D^2 + x*D + 1), where D = d/dx, factorizes as (D - a(x))*(D - b(x)), where a(x) = -(x + b(x)), as conjectured by [Edgar, Problem 4.32]. For a refinement of this sequence see A053979. (End)
From Sergei N. Gladkovskii, Aug 19 2012, Oct 24 2012, Mar 19 2013, May 20 2013, May 29 2013, Aug 04 2013, Aug 05 2013: (Start)
Continued fractions:
G.f.: 2 - G(0) where G(k) = 1 - (k+1)*x/G(k+1).
G.f.: 2 - U(0) where U(k) = 1 - (2*k+1)*x/(1 - (2*k+2)*x/U(k+1)).
G.f.: 2 - U(0) where U(k) = 1 - (4*k+1)*x - (2*k+1)*(2*k+2)*x^2/U(k+1).
G.f.: 1/Q(0) where Q(k) = 1 - x*(2*k+2)/(1 - x*(2*k+3)/Q(k+1)).
G.f.: 1 + x/Q(0) where Q(k) = 1 - x*(k+2)/Q(k+1).
G.f.: 2 - G(0)/2 where G(k) = 1 + 1/(1 - 2*x*(2*k+1)/(2*x*(2*k+1) - 1 + 2*x*(2*k+2)/ G(k+1))).
G.f.: 1 + x*G(0) where G(k) = 1 - x*(k+2)/(x*(k+2) - 1/G(k+1)).
G.f.: 2 - 1/B(x) where B(x) is the g.f. of A001147.
G.f.: 1 + x/(1-2*x*B(x)) where B(x) is the g.f. of A167872. (End)
a(n) ~ 2^(n+1/2) * n^n / exp(n). - Vaclav Kotesovec, Mar 10 2014
G.f.: 1 + x*(1/x + (sqrt(2/Pi) * exp(1/(2*x)) * sqrt(-1/x))/Erfc(sqrt(-1/x)/sqrt(2))) where Erfc(z) = 1 - Erf(z) is the complementary error function, and Erf(z) is the integral of the Gaussian distribution. This generating function is obtained from the generating functional of (4-dimensional) QED, evaluated in dimension 0 for the 2-point function, without the modification implementing Furry theorem. - Robert Coquereaux, Sep 14 2014
From Peter Bala, May 23 2017: (Start)
G.f. A(x) = 1 + x/(1 + x - 3*x/(1 + 3*x - 5*x/(1 + 5*x - 7*x/(1 + 7*x - ...)))).
A(x) = 1 + x/(1 + x - 3*x/(1 - 2*x/(1 - 5*x/(1 - 4*x/(1 - 7*x/(1 - 6*x/(1 - ...))))))). (End)

Extensions

Formula corrected by Ignacio D. Peixoto, Jun 23 2006
More terms from Sean A. Irvine, Feb 27 2011

A000260 Number of rooted simplicial 3-polytopes with n+3 nodes; or rooted 3-connected triangulations with 2n+2 faces; or rooted 3-connected trivalent maps with 2n+2 vertices.

Original entry on oeis.org

1, 1, 3, 13, 68, 399, 2530, 16965, 118668, 857956, 6369883, 48336171, 373537388, 2931682810, 23317105140, 187606350645, 1524813969276, 12504654858828, 103367824774012, 860593023907540, 7211115497448720, 60776550501588855
Offset: 0

Views

Author

Keywords

Comments

Number of rooted loopless planar maps with n edges. E.g., there are a(2)=3 loopless planar maps with 2 edges: two rooted paths (.-.-.) and one digon (.=.). - Valery A. Liskovets, Sep 25 2003
Number of intervals (i.e., ordered pairs (x,y) such that x<=y) in the Tamari lattice (rotation lattice of binary trees) of size n (see Pallo and Chapoton references). - Ralf Stephan, May 08 2007, Jean Pallo (Jean.Pallo(AT)u-bourgogne.fr), Sep 11 2007
Number of rooted triangulations of type [n, 0] (see Brown paper eq (4.8)). - Michel Marcus, Jun 23 2013
Equivalently, number of rooted bridgeless planar maps with n edges. - Noam Zeilberger, Oct 06 2016
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018
Number of uniquely sorted permutations of [2n+1] that avoid the pattern 231. Also the number of uniquely sorted permutations of [2n+1] that avoid 132. - Colin Defant, Jun 13 2019
The sequence 1,3,13,68,... appears naturally in integral geometry, namely in the algebra of unitarily invariant valuations on complex space forms. - Andreas Bernig, Feb 02 2020

Examples

			G.f. = 1 + x + 3*x^2 + 13*x^3 + 68*x^4 + 399*x^5 + 2530*x^6 + 16965*x^7 + ...
		

References

  • C. F. Earl and L. J. March, Architectural applications of graph theory, pp. 327-355 of R. J. Wilson and L. W. Beineke, editors, Applications of Graph Theory. Academic Press, NY, 1979.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 714.
  • Handbook of Combinatorics, North-Holland '95, p. 891.
  • 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).
  • W. T. Tutte, The enumerative theory of planar maps, in A Survey of Combinatorial Theory (J. N. Srivastava et al. eds.), pp. 437-448, North-Holland, Amsterdam, 1973.

Crossrefs

Row sums of A342981.
Column 0 of A146305 and A341856; Column 2 of A255918.
Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.

Programs

  • Magma
    [Binomial(4*n+1, n+1)-9*Binomial(4*n+1, n-1): n in [0..25]]; // Vincenzo Librandi, Nov 24 2016
  • Maple
    A000260 := proc(n)
        2*(4*n+1)!/((n+1)!*(3*n+2)!) ;
    end proc:
  • Mathematica
    Table[Binomial[4n+1,n+1]-9*Binomial[4n+1,n-1],{n,0,25}] (* Harvey P. Dale, Aug 23 2011 *)
    a[ n_] := SeriesCoefficient[ HypergeometricPFQ[ {1/2, 3/4, 1, 5/4}, {4/3, 5/3, 2}, 256/27 x], {x, 0, n}]; (* Michael Somos, Dec 23 2014 *)
    terms = 22; G[] = 0; Do[G[x] = 1+x*G[x]^4 + O[x]^terms, terms];
    CoefficientList[(2-G[x])*G[x]^2, x] (* Jean-François Alcover, Jan 13 2018, after Mark van Hoeij *)
  • PARI
    {a(n) = if( n<0, 0, 2 * (4*n + 1)! / ((n + 1)! * (3*n + 2)!))}; /* Michael Somos, Sep 07 2005 */
    
  • PARI
    {a(n) = binomial( 4*n + 2, n + 1) / ((2*n + 1) * (3*n + 2))}; /* Michael Somos, Mar 28 2012 */
    
  • Sage
    def a(n):
        n = ZZ(n)
        return (4*n + 2).binomial(n + 1) // ((2*n + 1) * (3*n + 2))
    # F. Chapoton, Aug 06 2015
    

Formula

a(n) = 2*(4*n+1)! / ((n+1)!*(3*n+2)!) = binomial(4*n+1, n+1) - 9*binomial(4*n+1, n-1).
G.f.: (2-g)*g^2 where g = 1 + x*g^4 is the g.f. of A002293. - Mark van Hoeij, Nov 10 2011
G.f.: hypergeom([1,1/2,3/4,5/4],[2,4/3,5/3],256*x/27) = 1 + 120*x/(Q(0)-120*x); Q(k) = 8*x*(2*k+1)*(4*k+3)*(4*k+5) + 3*(k+2)*(3*k+4)*(3*k+5) - 24*x*(k+2)*(2*k+3)*(3*k+4)*(3*k+5)*(4*k+7)*(4*k+9)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 25 2011
a(n) = binomial(4*n + 2, n + 1) / ((2*n + 1) * (3*n + 2)). - Michael Somos, Mar 28 2012
a(n) * (n+1) = A069271(n). - Michael Somos, Mar 28 2012
0 = F(a(n), a(n+1), ..., a(n+8)) for all n in Z where a(-1) = 3/4 and F() is a polynomial of degree 2 with integer coefficients and 29 monomials. - Michael Somos, Dec 23 2014
D-finite with recurrence: 3*(3*n+2)*(3*n+1)*(n+1)*a(n) - 8*(4*n+1)*(2*n-1)*(4*n-1)*a(n-1) = 0. - R. J. Mathar, Oct 21 2015
a(n) = Sum_{k=1..A000108(n)} k * A263191(n,k). - Alois P. Heinz, Nov 16 2015
a(n) ~ 2^(8*n+7/2) / (sqrt(Pi) * n^(5/2) * 3^(3*n+5/2)). - Vaclav Kotesovec, Feb 26 2016
E.g.f.: 3F3(1/2,3/4,5/4; 4/3,5/3,2; 256*x/27). - Ilya Gutkovskiy, Feb 01 2017
From Gheorghe Coserea, Aug 17 2017: (Start)
G.f. y(x) satisfies:
0 = x^3*y^4 + 3*x^2*y^3 + x*(8*x+3)*y^2 - (20*x-1)*y + 16*x-1.
0 = x*(256*x - 27)*deriv(y,x) - 8*x^2*y^3 - 25*x*y^2 + 4*(24*x-11)*y + 44.
(End)
From Karol A. Penson, Apr 06 2022: (Start)
a(n) = Integral_{x=0...256/27} x^n*W(x), where
W(x) = (sqrt(2)/Pi)*(h1(x) - h2(x) + h3(x)) and
h1(x) = 3F2([-6/12,-2/12, 2/12], [ 3/12, 9/12], (27*x)/256)/((x/2)^(1/2));
h2(x) = 3F2([-3/12, 1/12, 5/12], [ 6/12, 15/12], (27*x)/256)/(x^(1/4));
h3(x) = 3F2([ 3/12, 7/12, 11/12], [18/12, 21/12], (27*x)/256)/(x^(-1/4)*32).
This integral representation is unique as the solution of n-th Hausdorff power moment of the function W. Using only the definition of a(n), W(x) can be proven to be positive. W(x) is singular at x = 0 and for x > 0 is monotonically decreasing to zero at x = 256/27. (End)
a(n) = (1/27^n) * Product_{1 <= i <= j <= 3*n} (3*i + j + 3)/(3*i + j - 1). Cf. A006013. - Peter Bala, Feb 21 2023

Extensions

Edited by F. Chapoton, Feb 03 2011

A000168 a(n) = 2*3^n*(2*n)!/(n!*(n+2)!).

Original entry on oeis.org

1, 2, 9, 54, 378, 2916, 24057, 208494, 1876446, 17399772, 165297834, 1602117468, 15792300756, 157923007560, 1598970451545, 16365932856990, 169114639522230, 1762352559231660, 18504701871932430, 195621134074714260, 2080697516976506220, 22254416920705240440, 239234981897581334730, 2583737804493878415084
Offset: 0

Views

Author

Keywords

Comments

Number of rooted planar maps with n edges. - Don Knuth, Nov 24 2013
Number of rooted 4-regular planar maps with n vertices.
Also, number of doodles with n crossings, irrespective of the number of loops.
From Karol A. Penson, Sep 02 2010: (Start)
Integral representation as n-th moment of a positive function on the (0,12) segment of the x axis. This representation is unique as it is the solution of the Hausdorff moment problem.
a(n) = Integral_{x=0..12} ((x^n*(4/9)*(1 - x/12)^(3/2)) / (Pi*sqrt(x/3))). (End)
Also, the number of distinct underlying shapes of closed normal linear lambda terms of a given size, where the shape of a lambda term abstracts away from its variable binding. [N. Zeilberger, 2015] - N. J. A. Sloane, Sep 18 2016
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018
Number of well-labeled trees (Bona, 2015). - N. J. A. Sloane, Dec 25 2018

Examples

			G.f. = 1 + 2*x + 9*x^2 + 54*x^3 + 378*x^4 + 2916*x^5 + 24057*x^6 + 208494*x^7 + ...
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, pages 319, 353.
  • E. R. Canfield, Calculating the number of rooted maps on a surface, Congr. Numerantium, 76 (1990), 21-34.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 714.
  • V. A. Liskovets, A census of nonisomorphic planar maps, in Algebraic Methods in Graph Theory, Vol. II, ed. L. Lovasz and V. T. Sos, North-Holland, 1981.
  • V. A. Liskovets, Enumeration of nonisomorphic planar maps, Selecta Math. Sovietica, 4 (No. 4, 1985), 303-323.
  • 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

Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.
First row of array A101486.
Cf. A005470.
Rooted maps with n edges of genus g for 0 <= g <= 10: this sequence, A006300, A006301, A104742, A215402, A238355, A238356, A238357, A238358, A238359, A238360.

Programs

  • Magma
    [(2*Catalan(n)*3^n)/(n+2): n in [1..30]]; // Vincenzo Librandi, Sep 04 2014
  • Maple
    A000168:=n->2*3^n*(2*n)!/(n!*(n+2)!);
  • Mathematica
    Table[(2*3^n*(2n)!)/(n!(n+2)!),{n,0,20}] (* Harvey P. Dale, Jul 25 2011 *)
    a[ n_] := If[ n < 0, 0, 2 3^n (2 n)!/(n! (n + 2)!)] (* Michael Somos, Nov 25 2013 *)
    a[ n_] := SeriesCoefficient[ Hypergeometric2F1[ 1/2, 1, 3, 12 x], {x, 0, n}] (* Michael Somos, Nov 25 2013 *)
  • PARI
    {a(n) = if( n<0, 0, 2 * 3^n * (2*n)! / (n! * (n+2)!))}; /* Michael Somos, Nov 25 2013 */
    

Formula

G.f. A(z) satisfies A(z) = 1 - 16*z + 18*z*A(z) - 27*z^2*A(z)^2.
G.f.: F(1/2,1;3;12x). - Paul Barry, Feb 04 2009
a(n) = 2*3^n*A000108(n)/(n+2). - Paul Barry, Feb 04 2009
D-finite with recurrence: (n + 1) a(n) = (12 n - 18) a(n - 1). - Simon Plouffe, Feb 09 2012
G.f.: 1/54*(-1+18*x+(-(12*x-1)^3)^(1/2))/x^2. - Simon Plouffe, Feb 09 2012
0 = a(n)*(+144*a(n+1) - 42*a(n+2)) + a(n+1)*(+18*a(n+1) + a(n+2)) if n>=0. - Michael Somos, Jan 31 2014
a(n) ~ 2*(12^n)/((n^2+3*n)*sqrt(Pi*n)). - Peter Luschny, Nov 25 2015
E.g.f.: exp(6*x)*(12*x*BesselI(0,6*x) - (1 + 12*x)*BesselI(1,6*x))/(9*x). - Ilya Gutkovskiy, Feb 01 2017
From Amiram Eldar, Jan 08 2023: (Start)
Sum_{n>=0} 1/a(n) = 1887/1331 + 3240*arccosec(2*sqrt(3))/(1331*sqrt(11)).
Sum_{n>=0} (-1)^n/a(n) = 1563/2197 - 3240*arccosech(2*sqrt(3))/(2197*sqrt(13)). (End)

Extensions

More terms from Joerg Arndt, Feb 26 2014

A062980 a(0) = 1, a(1) = 5; for n > 1, a(n) = 6*n*a(n-1) + Sum_{k=1..n-2} a(k)*a(n-k-1).

Original entry on oeis.org

1, 5, 60, 1105, 27120, 828250, 30220800, 1282031525, 61999046400, 3366961243750, 202903221120000, 13437880555850250, 970217083619328000, 75849500508999712500, 6383483988812390400000, 575440151532675686278125, 55318762960656722780160000
Offset: 0

Views

Author

Michael Praehofer (praehofer(AT)ma.tum.de), Jul 24 2001

Keywords

Comments

Number of rooted unlabeled connected triangular maps on a compact closed oriented surface with 2n faces (and thus 3n edges). [Vidal]
Equivalently, the number of pair of permutations (sigma,tau) up to simultaneous conjugacy on a pointed set of size 6*n with sigma^3=tau^2=1, acting transitively and with no fixed point. [Vidal]
Also, the asymptotic expansion of the Airy function Ai'(x)/Ai(x) = -sqrt(x) - 1/(4x) + Sum_{n>=2} (-1)^n a(n) (4x)^ (1/2-3n/2). [Praehofer]
Maple 6 gives the wrong asymptotics of Ai'(x)=AiryAi(1,x) as x->oo apart from the 3rd term. Therefore asympt(AiryAi(1,x/4)/AiryAi(x/4),x); reproduces only the value a(1)=1 correctly.
Number of closed linear lambda terms (see [Bodini, Gardy, Jacquot, 2013] and [N. Zeilberger, 2015] references). - Pierre Lescanne, Feb 26 2017
Proved (bijection) by O. Bodini, D. Gardy, A. Jacquot (2013). - Olivier Bodini, Mar 30 2018
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018

Examples

			1 + 5*x + 60*x^2 + 1105*x^3 + 27120*x^4 + 828250*x^5 + 30220800*x^6 + ...
		

Crossrefs

Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.
With interspersed zeros column 3 of A380622.
Pointed version of A129114.
Connected pointed version of A129115.

Programs

  • Haskell
    a062980 n = a062980_list !! n
    a062980_list = 1 : 5 : f 2 [5,1] where
       f u vs'@(v:vs) = w : f (u + 1) (w : vs') where
         w = 6 * u * v + sum (zipWith (*) vs_ $ reverse vs_)
         vs_ = init vs
    -- Reinhard Zumkeller, Jun 03 2013
    
  • Maple
    a:= proc(n) option remember; `if`(n<2, 4*n+1,
          6*n*a(n-1) +add(a(k)*a(n-k-1), k=1..n-2))
        end:
    seq(a(n), n=0..25);  # Alois P. Heinz, Mar 31 2017
  • Mathematica
    max = 16; f[y_] := -Sqrt[x] - 1/(4*x) + Sum[(-1)^n*a[n]*(4*x)^(1/2 - 3*(n/2)), {n, 2, max}] /. x -> 1/y^2; s[y_] := Normal[ Series[ AiryAiPrime[x] / AiryAi[x], {x, Infinity, max + 7}]] /. x -> 1/y^2; sol = SolveAlways[ Simplify[ f[y] == s[y], y > 0], y] // First; Join[{1, 5}, Table[a[n], {n, 3, max}] /. sol] (* Jean-François Alcover, Oct 09 2012, from Airy function asymptotics *)
    a[0] = 1; a[n_] := a[n] = (6*(n-1)+4)*a[n-1] + Sum[a[i]*a[n-i-1], {i, 0, n-1}]; Table[a[n], {n, 0, 15}] (* Jean-François Alcover, Nov 29 2013, after Vladimir Reshetnikov *)
  • PARI
    {a(n) = local(A); n++; if( n<1, 0, A = vector(n); A[1] = 1; for( k=2, n, A[k] = (6*k - 8) * A[k-1] + sum( j=1, k-1, A[j] * A[k-j])); A[n])} /* Michael Somos, Jul 24 2011 */
    
  • Python
    from sympy.core.cache import cacheit
    @cacheit
    def a(n): return n*4 + 1 if n<2 else 6*n*a(n - 1) + sum(a(k)*a(n - k - 1) for k in range(1, n - 1))
    print([a(n) for n in range(21)]) # Indranil Ghosh, Aug 09 2017

Formula

With offset 1, then a(1) = 1 and, for n > 1, a(n) = (6*n-8)*a(n-1) + Sum_{k=1..n-1} a(k)*a(n-k) [Praehofer] [Martin and Kearney].
a(n) = (6/Pi^2)*Integral_{x=0..oo} ((4*x)^(3*n/2)/(Ai(x)^2 + Bi(x)^2)) dt. - Vladimir Reshetnikov, Sep 24 2013
a(n) ~ 3 * 6^n * n! / Pi. - Vaclav Kotesovec, Jan 19 2015
0 = 6*x^2*y' + x*y^2 + (4*x-1)*y + 1, where y(x) = Sum_{n>=0} a(n)*x^n. - Gheorghe Coserea, Apr 02 2017
From Peter Bala, May 21 2017: (Start)
G.f. as an S-fraction: A(x) = 1/(1 - 5*x/(1 - 7*x/(1 - 11*x/(1 - 13*x/(1 - ... - (6*n - 1)*x/(1 - (6*n + 1)*x/(1 - .... See Stokes.
x*A(x) = B(x)/(1 + 2*B(x)), where B(x) = x + 7*x^2 + 84*x^3 + 1463*x^4 + ... is the o.g.f. of A172455.
A(x) = 1/(1 + 2*x - 7*x/(1 - 5*x/(1 - 13*x/(1 - 11*x/(1 - ... - (6*n + 1)*x/(1 - (6*n - 1)*x/(1 - .... (End)

Extensions

Entry revised by N. J. A. Sloane based on comments from Samuel A. Vidal, Mar 30 2007

A000309 Number of rooted planar bridgeless cubic maps with 2n nodes.

Original entry on oeis.org

1, 1, 4, 24, 176, 1456, 13056, 124032, 1230592, 12629760, 133186560, 1436098560, 15774990336, 176028860416, 1990947110912, 22783499599872, 263411369705472, 3073132646563840, 36143187370967040, 428157758086840320, 5105072641718353920, 61228492804372561920
Offset: 0

Views

Author

Keywords

Comments

Also counts rooted planar non-separable triangulations with 3n edges. - Valery A. Liskovets, Dec 01 2003
Equivalently, rooted planar loopless triangulations with 2n triangles. - Noam Zeilberger, Oct 06 2016
Description trees of type (2,2) with n edges. (A description tree of type (a,b) is a rooted plane tree where every internal node is labeled by an integer between a and [b + sum of labels of its children], every leaf is labeled a, and the root is labeled [b + sum of labels of its children]. See Definition 1 and Section 5.2 of Cori and Schaeffer 2003.) - Noam Zeilberger, Oct 08 2017
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018

References

  • C. F. Earl and L. J. March, Architectural applications of graph theory, pp. 327-355 of R. J. Wilson and L. W. Beineke, editors, Applications of Graph Theory. Academic Press, NY, 1979.
  • 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

Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.

Programs

  • GAP
    List([0..20], n -> 2^(n+1)*Factorial(3*n)/(Factorial(n)* Factorial(2*n+2))); # G. C. Greubel, Nov 29 2018
  • Magma
    [2^(n+1)*Factorial(3*n)/(Factorial(n)*Factorial(2*n+2)): n in [0..20]]; // Vincenzo Librandi, Aug 10 2014
    
  • Maple
    a := n -> 2^(n+1)*(3*n)!/(n!*(2*n+2)!);
    A000309 := n -> -(-2)^(n-1)*(3*n+2)*hypergeom([-3*(n+1),-n,-n+1/3], [-n-1,-n-2/3], 1): seq(simplify(A000309(n)), n = 0..21); # Peter Luschny, Oct 28 2022
  • Mathematica
    f[n_] := 2^n(3n)!/((n + 1)!(2n + 1)!); Table[f[n], {n, 0, 19}] (* Robert G. Wilson v, Sep 21 2004 *)
    Join[{1},RecurrenceTable[{a[1]==1,a[n]==4a[n-1] Binomial[3n,3]/ Binomial[2n+2,3]}, a[n],{n,20}]] (* Harvey P. Dale, May 11 2011 *)
  • PARI
    a(n) = 2^(n+1)*(3*n)!/(n!*(2*n+2)!); \\ Michel Marcus, Aug 09 2014
    
  • Sage
    [2^n*factorial(3*n)/(factorial(n+1)*factorial(2*n+1))for n in range(20)] # G. C. Greubel Nov 29 2018
    

Formula

a(n) = 2^(n-1) * A000139(n) for n > 0.
a(n) = 4*a(n-1)*binomial(3*n, 3) / binomial(2*n+2, 3).
a(n) = 2^n*(3*n)!/ ( (n+1)!*(2*n+1)! ).
G.f.: (1/(6*x)) * (hypergeom([ -2/3, -1/3],[1/2],(27/2)*x)-1). - Mark van Hoeij, Nov 02 2009
a(n) ~ 3^(3*n+1/2)/(sqrt(Pi)*2^(n+2)*n^(5/2)). - Ilya Gutkovskiy, Oct 06 2016
D-finite with recurrence (n+1)*(2*n+1)*a(n) -3*(3*n-1)*(3*n-2)*a(n-1)=0. - R. J. Mathar, Nov 02 2018
a(n) = -(-2)^(n-1)*(3*n+2)*hypergeom([-3*(n+1),-n,-n+1/3], [-n-1,-n-2/3], 1). The a(n) are values of the polynomials A358091. - Peter Luschny, Oct 28 2022
From Karol A. Penson, Feb 24 2025: (Start)
G.f.: hypergeom([1/3, 2/3, 1], [3/2, 2], (27*z)/2).
G.f. A(z) satisfies: - 1 + 27*z + (-36*z + 1)*A(z) + 8*z*A(z)^2 + 16*z^2*A(z)^3 = 0.
G.f.: ((4*sqrt(4 - 54*z) + 12*i*sqrt(6)*sqrt(z))^(1/3)*(sqrt(z*(4 - 54*z)) - 9*i*sqrt(6)*z) + (4*sqrt(4 - 54*z) - 12*i*sqrt(6)*sqrt(z))^(1/3)*(9*i*sqrt(6)*z + sqrt(z*(4 - 54*z))) - 8*sqrt(z))/(48*z^(3/2)), where i = sqrt(-1) is the imaginary unit.
a(n) = Integral_{x=0..27/2} x^n*W(x), where W(x) = (6^(1/3)*(9 + sqrt(81 - 6*x))^(2/3)*(9*sqrt(3) - sqrt(27 - 2*x)) - 2^(2/3)*3^(1/6)*(27 + sqrt(81 - 6*x))*x^(1/3))/(48*Pi*(9 + sqrt(81 - 6*x))^(1/3)*x^(2/3)).
This integral representation is unique as W(x) is the solution of the Hausdorff power moment problem for x on (0, 27/2). Using only the definition of a(n), W(x) can be proven to be positive. W(x) is singular at x = 0, with singularity x^(-2/3), and for x > 0 is monotonically decreasing to zero at x = 27/2. (End)

Extensions

Definition clarified by Michael Albert, Oct 24 2008

A267827 Number of closed indecomposable linear lambda terms with 2n+1 applications and abstractions.

Original entry on oeis.org

1, 2, 20, 352, 8624, 266784, 9896448, 426577920, 20918138624, 1149216540160, 69911382901760, 4665553152081920, 338942971881472000, 26631920159494995968, 2250690001888540950528, 203595258621775065120768, 19629810220331494121865216
Offset: 0

Views

Author

Noam Zeilberger, Jan 21 2016

Keywords

Comments

A linear lambda term is indecomposable if it has no closed proper subterm.
Equivalently, number of closed bridgeless rooted trivalent maps (on compact oriented surfaces of arbitrary genus) with 2n+1 trivalent vertices (and 1 univalent vertex).
The September 2018 talk by Noam Zeilberger (see link to video) connects three topics (planar maps, Tamari lattices, lambda calculus) and eight sequences: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827. - N. J. A. Sloane, Sep 17 2018

Examples

			A(x) = 1 + 2*x + 20*x^2 + 352*x^3 + 8624*x^4 + 266784*x^5 + ...
		

Crossrefs

Sequences mentioned in the Noam Zeilberger 2018 video: A000168, A000260, A000309, A000698, A000699, A002005, A062980, A267827.

Programs

  • Mathematica
    a[0] = 1; a[1] = 2; a[n_] := a[n] = (6n-2) a[n-1] + Sum[(6k+2) a[k] a[n-1-k], {k, 1, n-2}];
    Table[a[n], {n, 0, 16}] (* Jean-François Alcover, Oct 16 2018, after Gheorghe Coserea *)
  • PARI
    seq(N) = {
      my(a = vector(N)); a[1] = 2;
      for(n=2, N,
        a[n] = (6*n-2)*a[n-1] + sum(k=1, n-2, (6*k+2)*a[k]*a[n-1-k]));
      concat(1,a);
    };
    seq(16)
    \\ test 1: y = x^2*subst(Ser(seq(201)),'x,-'x^6); 0 == x^5*y*y' + y - x^2
    \\ test 2: y = Ser(seq(201)); 0 == 6*y*y'*x^2 + 2*y^2*x - y + 1
    \\ Gheorghe Coserea, Nov 10 2017
    F(N) = {
      my(x='x+O('x^N), t='t, F0=x, F1=0, n=1);
      while(n++,
        F1 = t + x*(F0 - subst(F0,t,0))^2 + x*deriv(F0,t);
        if (F1 == F0, break()); F0 = F1;);
      F0;
    };
    seq(N) = my(v=Vec(subst(F(2*N+2),'t,0))); vector((#v+1)\2, n, v[2*n-1]);
    seq(16) \\ Gheorghe Coserea, Apr 01 2017

Formula

The o.g.f. f(z) = z + 2*z^3 + 20*z^5 + 352*z^7 + ... can be defined using a catalytic variable as f(z) = F(z,0), where F(z,x) satisfies the functional-differential equation F(z,x) = x + z*(F(z,x) - F(z,0))^2 + z*(d/dx)F(z,x).
From Gheorghe Coserea, Nov 10 2017: (Start)
0 = x^5*y*y' + y - x^2, where y(x) = x^2*A(-x^6).
0 = 6*y*y'*x^2 + 2*y^2*x - y + 1, where y(x) = A(x).
a(n) = (6*n-2)*a(n-1) + Sum_{k=1..n-2} (6*k+2)*a(k)*a(n-1-k), for n >= 2.
(End)
a(n) = A291843(3*n+1, 2*n), n >= 1. - Danny Rorabaugh, Nov 10 2017

A266240 Triangle read by rows: T(n,g) is the number of rooted 2n-face triangulations in an orientable surface of genus g.

Original entry on oeis.org

1, 4, 1, 32, 28, 336, 664, 105, 4096, 14912, 8112, 54912, 326496, 396792, 50050, 786432, 7048192, 15663360, 6722816, 11824384, 150820608, 544475232, 518329776, 56581525, 184549376, 3208396800, 17388675072, 30117189632, 11100235520, 2966845440
Offset: 0

Views

Author

Gheorghe Coserea, Dec 25 2015

Keywords

Comments

Row n contains floor((n+3)/2) terms.

Examples

			Triangle starts:
n\g    [0]          [1]          [2]          [3]          [4]
[0]    1;
[1]    4,           1;
[2]    32,          28;
[3]    336,         664,         105;
[4]    4096,        14912,       8112;
[5]    54912,       326496,      396792,      50050;
[6]    786432,      7048192,     15663360,    6722816;
[7]    11824384,    150820608,   544475232,   518329776,   56581525;
[8]    184549376,   3208396800,  17388675072, 30117189632, 11100235520;
[9]    ...
		

Crossrefs

Columns k=0-4 give: A002005, A269473, A269474, A269475, A269476.
Row sums give A062980.

Programs

  • Mathematica
    T[n_ /; n >= 0, g_] /; 0 <= g <= (n+1)/2 := f[n, g]/(3n+2); T[, ] = 0; f[n_ /; n >= 1, g_ /; g >= 0] := f[n, g] = 4*(3*n+2)/(n+1)*(n*(3*n-2)*f[n - 2, g-1] + Sum[f[i, h]*f[n-2-i, g-h], {i, -1, n-1}, {h, 0, g}]); f[-1, 0] = 1/2; f[0, 0] = 2; f[, ] = 0; Table[Table[T[n, g], {g, 0, Floor[(n + 1)/2]}], {n, 0, 9}] // Flatten (* Jean-François Alcover, Feb 27 2016 *)
  • PARI
    N = 10;
    m = matrix(N+2, N+2);
    mget(n,g) = {
      if (g < 0 || g > (n+1)/2, return(0));
      return(m[n+2,g+1]);
    }
    mset(n,g,v) = {
      m[n+2,g+1] = v;
    }
    Cubic() = {
      mset(-1,0,1/2);
      mset(0,0,2);
      for (n = 1, N,
      for (g = 0, (n+1)\2,
        my(t1 = n * (3*n-2) * mget(n-2, g-1),
           t2 = sum(i = -1, n-1, sum(h = 0, g,
                    mget(i,h) * mget(n-2-i, g-h))));
        mset(n, g, 4*(3*n+2)/(n+1) * (t1 + t2))));
      my(a = vector(N+1));
      for (n = 0, N,
        a[n+1] = vector(1 + (n+1)\2);
        for (g = 0, (n+1)\2,
             a[n+1][g+1] = mget(n, g));
        a[n+1] = a[n+1]/(3*n+2));
      return(a);
    }
    concat(Cubic())

Formula

T(n,g) = f(n,g)/(3*n+2) for all n >= 0 and 0 <= g <= (n+1)/2, where f(n,g) satisfies the quadratic recurrence equation f(n,g) = 4*(3*n+2)/(n+1)*(n*(3*n-2)*f(n-2,g-1) + Sum_{i=-1..n-1} Sum_{h=0..g} f(i,h)*f(n-2-i, g-h)) for n >= 1 and g >= 0 with the initial conditions f(-1,0)=1/2, f(0,0)=2 and f(n,g)=0 for g < 0 or g > (n+1)/2.
For column g, as n goes to infinity we have T(n,g) ~ 3*6^((g-1)/2) * t(g) * n^(5*(g-1)/2) * (12*sqrt(3))^n, where t(g) = (A269418(g)/A269419(g)) / (2^(g-2) * gamma((5*g-1)/2)) and gamma is the Gamma function. - Gheorghe Coserea, Feb 26 2016

A322928 a(0)=1; for n>0, a(n) is the number of rooted 3-regular maps with 2n vertices on the projective plane.

Original entry on oeis.org

1, 9, 118, 1773, 28650, 484578, 8457708, 151054173, 2745685954, 50606020854, 943283037684, 17746990547634, 336517405188900, 6423775409047716, 123332141503711704, 2379824766494404317, 46124764901514110898, 897483137740689843054, 17524230350476917414180
Offset: 0

Views

Author

Evgeniy Krasko, Dec 31 2018

Keywords

Crossrefs

Cf. A002005 (genus 0)

Programs

  • Maple
    A[0]:= 1: A[1]:= 9: A[2]:= 118: A[3]:= 1773: A[4]:= 28650: A[5]:= 484578:
    for n from 6 to 20 do
      A[n]:= 995328*(n - 4)*(3*n - 16)*(3*n - 14)*(3*n - 10)*(3*n - 8)*A[n - 6]/((n - 3)*(n - 2)*(n - 1)*n*(n + 1)) - 576*(3*n - 10)*(3*n - 8)*(108*n^2 - 648*n + 1049)*A[n - 4]/((n - 2)*(n - 1)*n*(n + 1)) + 12*(108*n^2 - 432*n + 505)*A[n - 2]/(n*(n + 1))
    od:
    seq(A[i],i=0..20); # Robert Israel, Dec 30 2022
  • Mathematica
    a[n_] := -((2^(2 n + 1) (3 n)!!)/((n + 1)! n!!)) + (3 2^(2 n))/(n + 1)!! Sum[(3^k (2 k - 1)!! (3 n - 2 k - 1)!!)/(2^k k! (n - k)!), {k, 0, n}];
    Table[a[n], {n, 0, 20}] (* Andrey Zabolotskiy, Dec 29 2022 *)

Formula

Theorem 3.3 gives an explicit formula.
From Vaclav Kotesovec, Dec 30 2022: (Start)
Recurrence: (n-3)*(n-2)*(n-1)*n*(n+1)*a(n) = 12*(n-3)*(n-2)*(n-1)*(108*n^2 - 432*n + 505)*a(n-2) - 576*(n-3)*(3*n - 10)*(3*n - 8)*(108*n^2 - 648*n + 1049)*a(n-4) + 995328*(n-4)*(3*n - 16)*(3*n - 14)*(3*n - 10)*(3*n - 8)*a(n-6).
a(n) ~ Gamma(1/4) * 2^(2*n - 5/4) * 3^(3*n/2 + 5/4) / (Pi * n^(5/4)) * (1 - 2^(7/4)*sqrt(Pi)/(Gamma(1/4)*3^(3/4)*n^(1/4)) + Pi/(Gamma(1/4)^2*sqrt(3*n))).
(End)

Extensions

Added initial term a(0)=1 to match Taylor series expansion in Theorem 3.3. - N. J. A. Sloane, Jan 11 2019
Terms a(11) and beyond from Andrey Zabolotskiy, Dec 29 2022

A002012 Almost trivalent maps.

Original entry on oeis.org

4, 32, 200, 1120, 5880, 29568, 144144, 686400, 3208920, 14780480, 67251184, 302865472, 1352078000, 5990745600, 26369978400, 115407434880, 502503206040, 2178032472000, 9401840170800, 40434981787200, 173319035569680, 740642835229440, 3156148445580000
Offset: 0

Views

Author

Keywords

References

  • R. C. Mullin, E. Nemeth and P. J. Schellenberg, The enumeration of almost cubic maps, pp. 281-295 in Proceedings of the Louisiana Conference on Combinatorics, Graph Theory and Computer Science. Vol. 1, edited R. C. Mullin et al., 1970.
  • 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

Formula

a(n) = 2*(n+3)*(2*(n+1))! / (3*n!*(n+1)!). [Mathai & Rathie, Eq. (22)] - Andrey Zabolotskiy, Jun 24 2024

Extensions

Terms a(7) and beyond from Andrey Zabolotskiy, Jun 24 2024
Showing 1-10 of 11 results. Next