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

A004040 Inversion of A000257.

Original entry on oeis.org

1, 2, 6, 23, 103, 512, 2740, 15485, 91245, 555662, 3475090, 22214707, 144640291, 956560748, 6411521056, 43478151737, 297864793993, 2059159989914, 14350039389022, 100726680316559, 711630547589023, 5057282786190872, 36132861123763276, 259423620328055093, 1870954187618001253
Offset: 0

Views

Author

Keywords

Crossrefs

Cf. A000257.
Essentially the same as A022558.

Programs

  • PARI
    x='x+O('x^66); Vec( -1 + 32*x/(1+20*x-8*x^2-(1-8*x)^(3/2)) ) \\ Joerg Arndt, May 04 2013

Formula

Invert: define b by 1+Sum a(n)x^n = 1/(1 - Sum b(n)x^n).

Extensions

More terms from Joerg Arndt, May 04 2013

A064062 Generalized Catalan numbers C(2; n).

Original entry on oeis.org

1, 1, 3, 13, 67, 381, 2307, 14589, 95235, 636925, 4341763, 30056445, 210731011, 1493303293, 10678370307, 76957679613, 558403682307, 4075996839933, 29909606989827, 220510631755773, 1632599134961667, 12133359132082173
Offset: 0

Views

Author

Wolfdieter Lang, Sep 13 2001

Keywords

Comments

a(n+1) = Y_{n}(n+1) = Z_{n}, n >= 0, in the Derrida et al. 1992 reference (see A064094) for alpha=2, beta=1 (or alpha=1, beta=2).
a(n) = number of Dyck n-paths (A000108) in which each upstep (U) not at ground level is colored red (R) or blue (B). For example, a(3)=3 counts URDD, UBDD, UDUD (D=downstep). - David Callan, Mar 30 2007
The Hankel transform of this sequence is A002416. - Philippe Deléham, Nov 19 2007
The sequence a(n)/2^n, with g.f. 1/(1-xc(x)/2), has Hankel transform 1/2^n. - Paul Barry, Apr 14 2008
The REVERT transform of the odd numbers [1,3,5,7,9,...] is [1, -3, 13, -67, 381, -2307, 14589, -95235, 636925, ...] - N. J. A. Sloane, May 26 2017

Crossrefs

Generalized Catalan numbers C(m; n): A000012 (m = 0), A000108 (m = 1), A064063 (m = 3) and A064087 - A064093 (m = 4 thru 10); A064310 (m = -1), A064311 (m = -2) and A064325 - A064333 (m = -3 thru -11).

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30);
    Coefficients(R!( (3 - Sqrt(1-8*x))/(2*(1+x)) )); // G. C. Greubel, Sep 27 2024
  • Maple
    1, seq(simplify(hypergeom([1-n,n],[-n],2)), n=1..100); # Robert Israel, Nov 30 2014
  • Mathematica
    a[0]=1; a[1]=1; a[n_]/;n>=2 := a[n] = a[n-1] + Sum[(a[k] + a[k-1])a[n-k],{k,n-1}]; Table[a[n],{n,0,10}] (* David Callan, Aug 27 2009 *)
    a[n_] := 2*Sum[ (-1)^j*2^(n-j-1)*Binomial[2*(n-j-1), n-j-1]/(n-j), {j, 0, n-1}] + (-1)^n; Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Jul 03 2013 *)
  • PARI
    {a(n)=polcoeff((3-sqrt(1-8*x+x*O(x^n)))/(2+2*x),n)}
    
  • PARI
    {a(n)=local(A=1+x); for(i=1, n, A=1+A^4*intformal(1/(A^2+x*O(x^n)))); polcoeff(A, n)} \\ Paul D. Hanna, Dec 24 2013
    for(n=0, 25, print1(a(n), ", "))
    
  • PARI
    {a(n)=polcoeff(1/(1 - serreverse(x-2*x^2 +x^2*O(x^n))),n)}
    for(n=0,30,print1(a(n),", ")) \\ Paul D. Hanna, Nov 30 2014
    
  • Sage
    def a(n):
        if n==0: return 1
        return hypergeometric([1-n, n], [-n], 2).simplify()
    [a(n) for n in range(22)] # Peter Luschny, Dec 01 2014
    

Formula

G.f.: (1 + 2*x*C(2*x)) / (1+x) = 1/(1 - x*C(2*x)) with C(x) g.f. of Catalan numbers A000108.
a(n) = A062992(n-1) = Sum_{m = 0..n-1} (n-m)*binomial(n-1+m, m)*(2^m)/n, n >= 1, a(0) = 1.
a(n) = Sum_{k = 0..n} A059365(n, k)*2^(n-k). - Philippe Deléham, Jan 19 2004
G.f.: 1/(1-x/(1-2x/(1-2x/(1-2x/(1-.... = 1/(1-x-2x^2/(1-4x-4x^2/(1-4x-4x^2/(1-.... (continued fractions). - Paul Barry, Jan 30 2009
a(n) = (32/Pi)*Integral_{x = 0..1} (8*x)^(n-1)*sqrt(x*(1-x)) / (8*x+1). - Groux Roland, Dec 12 2010
a(n+2) = 8^(n+2)*( c(n+2)-c(1)*c(n+1) - Sum_{i=0..n-1} 8^(-i-2)*c(n-i)*a(i+2) ) with c(n) = Catalan(n+2)/2^(2*n+1). - Groux Roland, Dec 12 2010
a(n) = the upper left term in M^n, M = the production matrix:
1, 1
2, 2, 1
4, 4, 2, 1
8, 8, 4, 2, 1
... - Gary W. Adamson, Jul 08 2011
D-finite with recurrence: n*a(n) + (12-7n)*a(n-1) + 4*(3-2n)*a(n-2) = 0. - R. J. Mathar, Nov 16 2011 (This follows easily from the generating function. - Robert Israel, Nov 30 2014)
G.f. satisfies: A(x) = 1 + A(x)^4 * Integral 1/A(x)^2 dx. - Paul D. Hanna, Dec 24 2013
G.f. satisfies: Integral 1/A(x)^2 dx = x - x^2*G(x), where G(x) is the o.g.f. of A000257, the number of rooted bicubic maps. - Paul D. Hanna, Dec 24 2013
G.f. A(x) satisfies: A(x - 2*x^2) = 1/(1-x). - Paul D. Hanna, Nov 30 2014
a(n) = hypergeometric([1-n, n], [-n], 2) for n > 0. - Peter Luschny, Nov 30 2014
G.f.: (3 - sqrt(1-8*x))/(2*(x+1)). - Robert Israel, Nov 30 2014
a(n) ~ 2^(3*n+1) / (9*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Dec 22 2014
O.g.f. A(x) = 1 + series reversion of (x*(1 - x)/(1 + x)^2). Logarithmically differentiating (A(x) - 1)/x gives 3 + 17*x + 111*x^2 + ..., essentially a g.f for A119259. - Peter Bala, Oct 01 2015
From Peter Bala, Jan 06 2022: (Start)
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + x + 2*x^2 + 6*x^3 + 23*x^4 + ... is a g.f. for A022558.
The Gauss congruences a(n*p^k) == a(n^p^(k-1)) (mod p^k) hold for prime p and positive integers n and k. (End)

A007054 Super ballot numbers: 6(2n)!/(n!(n+2)!).

Original entry on oeis.org

3, 2, 3, 6, 14, 36, 99, 286, 858, 2652, 8398, 27132, 89148, 297160, 1002915, 3421710, 11785890, 40940460, 143291610, 504932340, 1790214660, 6382504440, 22870640910, 82334307276, 297670187844, 1080432533656, 3935861372604, 14386251913656, 52749590350072
Offset: 0

Views

Author

Keywords

Comments

Hankel transform is 2n+3. The Hankel transform of a(n+1) is n+2. The sequence a(n)-2*0^n has Hankel transform A110331(n). - Paul Barry, Jul 20 2008
Number of pairs of Dyck paths of total length 2*n with heights differing by at most 1 (Gessel/Xin, p. 2). - Joerg Arndt, Sep 01 2012

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    [6*Factorial(2*n)/(Factorial(n)*Factorial(n+2)): n in [0..30]]; // Vincenzo Librandi, Aug 20 2011
    
  • Maple
    seq(3*(2*n)!/(n!)^2/binomial(n+2,n), n=0..22); # Zerinvary Lajos, Jun 28 2007
    A007054 := n -> 6*4^n*GAMMA(1/2+n)/(sqrt(Pi)*GAMMA(3+n)):
    seq(A007054(n),n=0..28); # Peter Luschny, Dec 14 2015
  • Mathematica
    Table[6(2n)!/(n!(n+2)!),{n,0,30}] (* or *) CoefficientList[Series[ (-1+Sqrt[1-4*x]+(6-4*Sqrt[1-4*x])*x)/(2*x^2),{x,0,30}],x] (* Harvey P. Dale, Oct 05 2011 *)
  • PARI
    a(n)=6*(2*n)!/(n!*(n+2)!); /* Joerg Arndt, Sep 01 2012 */
    
  • Sage
    def A007054(n): return (-4)^(2 + n)*binomial(3/2, 2 + n)/2
    print([A007054(n) for n in range(29)])  # Peter Luschny, Nov 04 2021

Formula

G.f.: c(x)*(4-c(x)), where c(x) = g.f. for Catalan numbers A000108; Convolution of Catalan numbers with negative Catalan numbers but -C(0)=-1 replaced by 3. - Wolfdieter Lang
E.g.f. in Maple notation: exp(2*x)*(4*x*(BesselI(0, 2*x)-BesselI(1, 2*x))-BesselI(1, 2*x))/x. Integral representation as n-th moment of a positive function on [0, 4], in Maple notation: a(n)=int(x^n*(4-x)^(3/2)/x^(1/2), x=0..4)/(2*Pi), n=0, 1, ... This representation is unique. - Karol A. Penson, Oct 10 2001
E.g.f.: Sum_{n>=0} a(n)*x^(2*n) = 3*BesselI(2, 2x).
a(n) = A000108(n)*6/(n+2). - Philippe Deléham, Oct 30 2007
a(n+1) = 2*(A000108(n+2) - A000108(n+1))/(n+1). - Paul Barry, Jul 20 2008
G.f.: ((6-4*sqrt(1-4*x))*x+sqrt(1-4*x)-1)/(2*x^2) - Harvey P. Dale, Oct 05 2011
a(n) = 4*A000108(n) - A000108(n+1) (Gessel/Xin, p. 2). - Joerg Arndt, Sep 01 2012
D-finite with recurrence (n+2)*a(n) +2*(-2*n+1)*a(n-1)=0. - R. J. Mathar, Dec 03 2012
G.f.: 1/(x^2*G(0)) + 3/x - (1/2)/x^2, where G(k) = 1 + 1/(1 - 2*x*(2*k+3)/(2*x*(2*k+3) + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 06 2013
G.f.: 3/x - 1/(2*x^2) + G(0)/(4*x^2), where G(k) = 1 + 1/(1 - 2*x*(2*k-3)/(2*x*(2*k-3) + (k+1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 18 2013
0 = a(n)*(+16*a(n+1) - 14*a(n+2)) + a(n+1)*(+6*a(n+1) + a(n+2)) for all n in Z. - Michael Somos, Sep 18 2014
A002421(n+2) = 2*a(n) for all n in Z. - Michael Somos, Sep 18 2014
a(n) = 3*(2*n)!*[x^(2*n)]hypergeometric([],[3],x^2). - Peter Luschny, Feb 01 2015
a(n) = 6*4^n*Gamma(1/2+n)/(sqrt(Pi)*Gamma(3+n)). - Peter Luschny, Dec 14 2015
a(n) = (-4)^(2 + n)*binomial(3/2, 2 + n)/2. - Peter Luschny, Nov 04 2021
From Amiram Eldar, May 16 2022: (Start)
Sum_{n>=0} 1/a(n) = 1 + 20*Pi/(81*sqrt(3)).
Sum_{n>=0} (-1)^n/a(n) = 3/25 - 8*log(phi)/(25*sqrt(5)), where phi is the golden ratio (A001622). (End)
a(n-1) = 3*A000984(n)/((2*n-1)*(n+1)). - R. J. Mathar, Jul 12 2024

Extensions

Corrected and extended by Vincenzo Librandi, Aug 20 2011

A002421 Expansion of (1-4*x)^(3/2) in powers of x.

Original entry on oeis.org

1, -6, 6, 4, 6, 12, 28, 72, 198, 572, 1716, 5304, 16796, 54264, 178296, 594320, 2005830, 6843420, 23571780, 81880920, 286583220, 1009864680, 3580429320, 12765008880, 45741281820, 164668614552, 595340375688, 2160865067312, 7871722745208, 28772503827312
Offset: 0

Views

Author

Keywords

Comments

Terms that are not divisible by 12 have indices in A019469. - Ralf Stephan, Aug 26 2004
From Ralf Steiner, Apr 06 2017: (Start)
By analytic continuation to the entire complex plane there exist regularized values for divergent sums such as:
Sum_{k>=0} a(k)^2/8^k = 2F1(-3/2,-3/2,1,2).
Sum_{k>=0} a(k) / 2^k = -i. (End)

Examples

			G.f. = 1 - 6*x + 6*x^2 + 4*x^3 + 6*x^4 + 12*x^5 + 28*x^6 + 72*x^7 + 198*x^8 + 572*x^9 + ...
		

References

  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 55.
  • 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).
  • T. N. Thiele, Interpolationsrechnung. Teubner, Leipzig, 1909, p. 164.

Crossrefs

Programs

  • GAP
    Concatenation([1], List([1..40], n-> 12*Factorial(2*n-4) /( Factorial(n)*Factorial(n-2)) )) # G. C. Greubel, Jul 03 2019
  • Magma
    [1,-6] cat [12*Catalan(n-2)/n: n in [2..30]]; // Vincenzo Librandi, Jun 11 2012
    
  • Maple
    A002421 := n -> 3*4^(n-1)*GAMMA(-3/2+n)/(sqrt(Pi)*GAMMA(1+n)):
    seq(A002421(n), n=0..29); # Peter Luschny, Dec 14 2015
  • Mathematica
    CoefficientList[Series[(1-4x)^(3/2),{x,0,40}],x] (* Vincenzo Librandi, Jun 11 2012 *)
    a[n_]:= Binomial[ 3/2, n] (-4)^n; (* Michael Somos, Dec 04 2013 *)
    a[n_]:= SeriesCoefficient[(1-4x)^(3/2), {x, 0, n}]; (* Michael Somos, Dec 04 2013 *)
  • PARI
    {a(n) = binomial( 3/2, n) * (-4)^n}; /* Michael Somos, Dec 04 2013 */
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( (1 - 4*x + x * O(x^n))^(3/2), n))}; /* Michael Somos, Dec 04 2013 */
    
  • Sage
    ((1-4*x)^(3/2)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Jul 03 2019
    

Formula

a(n) = Sum_{m=0..n} binomial(n, m)*K_m(4), where K_m(x) = K_m(n, 2, x) is a Krawtchouk polynomial. - Alexander Barg (abarg(AT)research.bell-labs.com)
a(n) ~ (3/4)*Pi^(-1/2)*n^(-5/2)*2^(2*n)*(1 + 15/8*n^-1 + ...). - Joe Keane (jgk(AT)jgk.org), Nov 22 2001
From Ralf Stephan, Mar 11 2004: (Start)
a(n) = 12*(2*n-4)! /(n!*(n-2)!), n > 1.
a(n) = 12*Cat(n-2)/n = 2(Cat(n-1) - 4*Cat(n-2)), in terms of Catalan numbers (A000108).
Terms that are not divisible by 12 have indices in A019469. (End)
Let rho(x)=(1/Pi)*(x*(4-x))^(3/2), then for n >= 4, a(n) = Integral_{x=0..4} (x^(n-4) *rho(x)) dx. - Groux Roland, Mar 16 2011
G.f.: (1-4*x)^(3/2) = 1 - 6*x + 12*x^2/(G(0) + 2*x); G(k) = (4*x+1)*k-2*x+2-2*x*(k+2)*(2*k+1)/G(k+1); for -1/4 <= x < 1/4, otherwise G(0) = 2*x; (continued fraction). - Sergei N. Gladkovskii, Dec 05 2011
G.f.: 1/G(0) where G(k) = 1 + 4*x*(2*k+1)/(1 - 1/(1 + (2*k+2)/G(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Nov 18 2012
G.f.: G(0)/2, where G(k) = 2 + 2*x*(2*k-3)*G(k+1)/(k+1). - Sergei N. Gladkovskii, Jun 06 2013 [Edited by Michael Somos, Dec 04 2013]
0 = a(n+2) * (a(n+1) - 14*a(n)) + a(n+1) * (6*a(n+1) + 16*a(n)) for all n in Z. - Michael Somos, Dec 04 2013
A232546(n) = 3^n * a(n). - Michael Somos, Dec 04 2013
G.f.: hypergeometric1F0(-3/2;;4*x). - R. J. Mathar, Aug 09 2015
a(n) = 3*4^(n-1)*Gamma(-3/2+n)/(sqrt(Pi)*Gamma(1+n)). - Peter Luschny, Dec 14 2015
From Ralf Steiner, Apr 06 2017: (Start)
Sum_{k>=0} a(k)/4^k = 0.
Sum_{k>=0} a(k)^2/16^k = 32/(3*Pi).
Sum_{k>=0} a(k)^2*(k/8)/16^k = 1/Pi.
Sum_{k>=0} a(k)^2*(-k/24+1/8)/16^k = 1/Pi.
Sum_{k>=0} a(k-1)^2*(k-1/4)/16^k = 1/Pi.
Sum_{k>=0} a(k-1)^2*(2k-2)/16^k = 1/Pi.(End)
D-finite with recurrence: n*a(n) +2*(-2*n+5)*a(n-1)=0. - R. J. Mathar, Feb 20 2020
From Amiram Eldar, Mar 22 2022: (Start)
Sum_{n>=0} 1/a(n) = 4/3 + 10*Pi/(81*sqrt(3)).
Sum_{n>=0} (-1)^n/a(n) = 92/75 - 4*sqrt(5)*log(phi)/125, where phi is the golden ratio (A001622). (End)

A090371 Number of unrooted planar 2-constellations with n digons. Also number of n-edge unrooted planar Eulerian maps with bicolored faces.

Original entry on oeis.org

1, 3, 6, 20, 60, 291, 1310, 6975, 37746, 215602, 1262874, 7611156, 46814132, 293447817, 1868710728, 12068905911, 78913940784, 521709872895, 3483289035186, 23464708686960, 159346213738020, 1090073011199451, 7507285094455566, 52021636161126702
Offset: 1

Views

Author

Valery A. Liskovets, Dec 01 2003

Keywords

Comments

a(n) is also the number of unrooted planar hypermaps with n darts up to orientation-preserving homeomorphism (darts are semi-edges in the particular case of ordinary maps). - Valery A. Liskovets, Apr 13 2006

Examples

			The 3 Eulerian maps with 2 edges are the digon and two figure eight graphs ("8") in which both loops are colored, resp., black or white.
		

Crossrefs

Programs

  • Maple
    A090371 := proc(n)
        local s, d;
        if n=0 then
            1 ;
        else
            s := -2^n*binomial(2*n, n);
            for d in numtheory[divisors](n) do
                s := s+ numtheory[phi](n/d)*2^d*binomial(2*d, d)
            od;
            3/(2*n)*(2^n*binomial(2*n, n)/((n+1)*(n+2))+s/2);
        fi;
    end proc:
  • Mathematica
    h0[n_] := 3*2^(n-1)*Binomial[2*n, n]/((n+1)*(n+2)); a[n_] := (h0[n] + DivisorSum[n, If[#>1, EulerPhi[#]*Binomial[n/#+2, 2]*h0[n/#], 0]&])/n; Array[a, 30] (* Jean-François Alcover, Dec 06 2015, adapted from PARI *)
  • PARI
    h0(n) = 3*2^(n-1)*binomial(2*n, n)/((n+1)*(n+2));
    a(n) = (h0(n) + sumdiv(n, d, (d>1)*eulerphi(d)*binomial(n/d+2,2)*h0(n/d)))/n; \\ Michel Marcus, Dec 11 2014

Extensions

More terms from Michel Marcus, Dec 11 2014

A069727 Number of nonisomorphic unrooted Eulerian planar maps with n edges (Eulerian means that all vertices are of even valency; there is an Eulerian cycle).

Original entry on oeis.org

1, 1, 2, 4, 12, 34, 154, 675, 3534, 18985, 108070, 632109, 3807254, 23411290, 146734695, 934382820, 6034524474, 39457153432, 260855420489, 1741645762265, 11732357675908, 79673115468562, 545036528857605, 3753642607424647, 26010818244754788, 181266500331748878
Offset: 0

Views

Author

Valery A. Liskovets, Apr 07 2002

Keywords

Comments

By duality, also the number of unrooted (sensed) bipartite maps with n edges. - Andrew Howroyd, Mar 29 2021

Crossrefs

Cf. A000257 (rooted), A069720, A069724, A103939 (with distinguished face), A103940 (with distinguished vertex).

Programs

  • Mathematica
    a[n_] := (1/(2n)) * (3*2^(n-1) * Binomial[2n, n]/((n+1)*(n+2)) + Sum[ EulerPhi[n/k] * d[n/k] * 2^(k-2) * Binomial[2k, k], {k, Most[ Divisors[n]]}]) + q[n]; a[0] = 1; q[n_?EvenQ] := 2^((n-4)/2)*Binomial[ n, n/2]/(n+2); q[n_?OddQ] := 2^((n-1)/2)*Binomial[(n-1), (n-1)/2]/(n+1); d[n_] := 4-Mod[n, 2]; Table[ a[n], {n, 0, 20}] (* Jean-François Alcover, Dec 19 2011, after Valery A. Liskovets *)
  • PARI
    a(n) = {if(n==0, 1, sumdiv(n, d, if(dAndrew Howroyd, Mar 29 2021

Formula

a(n) = (1/(2n))*(3*2^(n-1)*binomial(2n, n)/((n+1)(n+2)) + Sum_{k=1..n-1, k|n} phi(n/k)*d(n/k)*2^(k-2)*binomial(2k, k)) + q(n) where phi is the Euler function A000010, q(n) = 2^((n-4)/2)*binomial(n, n/2)/(n+2) if n is even, q(n) = 2^((n-1)/2)*binomial(n-1, (n-1)/2)/(n+1) if n is odd, d(n)=4, if n is even and d(n)=3 if n is odd. - Valery A. Liskovets, Mar 17 2005
a(n) ~ 3 * 2^(3*n-2) / (sqrt(Pi) * n^(7/2)). - Vaclav Kotesovec, Aug 28 2019

A321710 Triangle read by rows: T(n,k) is the number of rooted hypermaps of genus k with n darts.

Original entry on oeis.org

1, 3, 12, 1, 56, 15, 288, 165, 8, 1584, 1611, 252, 9152, 14805, 4956, 180, 54912, 131307, 77992, 9132, 339456, 1138261, 1074564, 268980, 8064, 2149888, 9713835, 13545216, 6010220, 579744, 13891584, 81968469, 160174960, 112868844, 23235300, 604800, 91287552, 685888171, 1805010948, 1877530740, 684173164, 57170880, 608583680, 5702382933, 19588944336, 28540603884, 16497874380, 2936606400, 68428800, 4107939840, 47168678571, 206254571236, 404562365316, 344901105444, 108502598960, 8099018496
Offset: 1

Views

Author

Gheorghe Coserea, Nov 17 2018

Keywords

Comments

Row n contains floor((n+1)/2) = A008619(n-1) terms.

Examples

			Triangle starts:
n\k  [0]       [1]        [2]         [3]         [4]        [5]
[1]  1;
[2]  3;
[3]  12,       1;
[4]  56,       15;
[5]  288,      165,       8;
[6]  1584,     1611,      252;
[7]  9152,     14805,     4956,       180;
[8]  54912,    131307,    77992,      9132;
[9]  339456,   1138261,   1074564,    268980,     8064;
[10] 2149888,  9713835,   13545216,   6010220,    579744;
[11] 13891584, 81968469,  160174960,  112868844,  23235300,  604800;
[12] 91287552, 685888171, 1805010948, 1877530740, 684173164, 57170880;
[13] ...
		

Crossrefs

Columns k=0..9 give: A000257 (k=0), A118093 (k=1), A214817 (k=2), A214818 (k=3), A318104 (k=4), A321705 (k=5), A321706 (k=6), A321707 (k=7), A321708 (k=8), A321709 (k=9).
Row sums give A003319(n+1).

Programs

  • Mathematica
    l1[f_,n_] := Sum[(i-1)t[i]D[f,t[i-1]], {i,2,n}];
    m1[f_,n_] := Sum[(i-1)t[j]t[i-j]D[f,t[i-1]] + j(i-j)t[i+1]D[f,t[j],t[i-j]], {i,2,n},{j,i-1}];
    ff[1] = x^2 t[1];
    ff[n_] := ff[n] = Simplify@(2x l1[ff[n-1],n] + m1[ff[n-1],n] + Sum[t[i+1]j(i-j)D[ff[k],t[j]]D[ff[n-1-k],t[i-j]], {i,2,n-1},{j,i-1},{k,n-2}]) / n;
    row[n_]:=Reverse[CoefficientList[n ff[n] /. {t[_]->x}, x]][[;;;;2]][[;;Quotient[n+1,2]]];
    Table[row[n], {n,14}] (* Andrei Zabolotskii, Jun 27 2025, after the PARI code *)
  • PARI
    L1(f, N) = sum(i=2, N, (i-1)*t[i]*deriv(f, t[i-1]));
    M1(f, N) = {
      sum(i=2, N, sum(j=1, i-1, (i-1)*t[j]*t[i-j]*deriv(f, t[i-1]) +
          j*(i-j)*t[i+1]*deriv(deriv(f, t[j]), t[i-j])));
    };
    F(N) = {
      my(u='x, v='x, f=vector(N)); t=vector(N+1, n, eval(Str("t", n)));
      f[1] = u*v*t[1];
      for (n=2, N, f[n] = (u + v)*L1(f[n-1], n) + M1(f[n-1], n) +
        sum(i=2, n-1, t[i+1]*sum(j=1, i-1,
        j*(i-j)*sum(k=1, n-2, deriv(f[k], t[j])*deriv(f[n-1-k], t[i-j]))));
        f[n] /= n);
      f;
    };
    seq(N) = {
      my(f=F(N), v=substvec(f, t, vector(#t, n, 'x)),
         g=vector(#v, n, Polrev(Vec(n * v[n]))));
      apply(p->Vecrev(substpol(p, 'x^2, 'x)), g);
    };
    concat(seq(14))

Formula

A000257(n)=T(n,0), A118093(n)=T(n,1), A214817(n)=T(n,2), A214818(n)=T(n,3), A060593(n)=T(2*n+1,n)=(2*n)!/(n+1), A003319(n+1)=Sum_{k=0..floor((n-1)/2)} T(n,k).

A118093 Numbers of rooted hypermaps on the torus with n darts (darts are semi-edges in the particular case of ordinary maps).

Original entry on oeis.org

1, 15, 165, 1611, 14805, 131307, 1138261, 9713835, 81968469, 685888171, 5702382933, 47168678571, 388580070741, 3190523226795, 26124382262613, 213415462218411, 1740019150443861, 14162920013474475, 115112250539595093, 934419385591442091, 7576722323539318101
Offset: 3

Views

Author

Valery A. Liskovets, Apr 13 2006

Keywords

Crossrefs

Programs

  • Magma
    [&+[(2^k*(4^(n-2-k)-1)*Binomial(n+k, k))/3 : k in [0..n-3]]: n in [3..25]]; // Vincenzo Librandi, Sep 16 2018
  • Mathematica
    Table[Sum[2^k (4^(n - 2 - k) - 1) Binomial[n+k, k] / 3, {k, 0, n-3}], {n, 3, 25}] (* Vincenzo Librandi, Sep 16 2018 *)
  • PARI
    a(n) = sum(k=0, n-3, 2^k*(4^(n-2-k)-1)*binomial(n+k, k))/3; \\ Michel Marcus, Dec 11 2014
    
  • PARI
    seq(N) = {
      my(x='x+O('x^(N+2)), y=(1-sqrt(1-8*x))/(4*x));
      Vec((y - 1)^3/(4*(y - 2)^2*(y + 1)));
    };
    seq(21) \\ Gheorghe Coserea, Nov 06 2018
    

Formula

Conjecture: +n*(5*n-17)*a(n) -15*(n-1)*(5*n-16)*a(n-1) +12*(20*n^2-103*n+140)*a(n-2) +32*(5*n-12)*(2*n-5)*a(n-3)=0. - R. J. Mathar, Apr 05 2018
G.f.: (1 - 7*x + 4*x^2 - (1 - 3*x)*sqrt(1 - 8*x))/(8*(1 + x)*(1 - 8*x)); equivalently, the g.f. can be rewritten as (y - 1)^3/(4*(y - 2)^2*(y + 1)), where y=G(2*x) with G the g.f. of A000108. - Gheorghe Coserea, Nov 06 2018
a(n) ~ 2^(3*n - 4) / 3 * (1 - 10/(3*sqrt(Pi*n))). - Vaclav Kotesovec, Nov 06 2018

Extensions

More terms from Michel Marcus, Dec 11 2014

A214817 Number of genus 2 rooted hypermaps with n darts.

Original entry on oeis.org

0, 0, 0, 0, 8, 252, 4956, 77992, 1074564, 13545216, 160174960, 1805010948, 19588944336, 206254571236, 2118399516180, 21310566266640, 210636265153004, 2050696768165560, 19704531058696008, 187168609978022860, 1759888050471704664, 16398685297890141180, 151570887948878270348
Offset: 1

Views

Author

N. J. A. Sloane, Jul 31 2012

Keywords

Comments

The table in the Zograf paper has an incorrect value for a(14). - Gheorghe Coserea, Nov 11 2018

Crossrefs

Programs

  • Mathematica
    DeleteCases[CoefficientList[Series[-# (# - 1)^5*(#^4 - 6 #^3 + 36 #^2 - 50 # + 51)/(4 (# - 2)^7*(# + 1)^5) &[(1 - Sqrt[1 - 8 x])/(4 x)], {x, 0, 23}], x], 0] (* Michael De Vlieger, Nov 26 2018 *)
  • PARI
    seq(N) = {
      my(x='x+O('x^(N+2)), y=(1-sqrt(1-8*x))/(4*x));
      Vec(-y*(y - 1)^5*(y^4 - 6*y^3 + 36*y^2 - 50*y + 51)/(4*(y - 2)^7*(y + 1)^5));
    };
    seq(19) \\ Gheorghe Coserea, Nov 11 2018

Formula

G.f.: -y*(y - 1)^5*(y^4 - 6*y^3 + 36*y^2 - 50*y + 51)/(4*(y - 2)^7*(y + 1)^5), where y = C(2*x), C being the g.f. for A000108. - Gheorghe Coserea, Nov 11 2018

Extensions

a(13) by Noam Zeilberger, Sep 16 2018
More terms and a(14) corrected by Gheorghe Coserea, Nov 11 2018

A214818 Number of genus 3 rooted hypermaps with n darts.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 180, 9132, 268980, 6010220, 112868844, 1877530740, 28540603884, 404562365316, 5422718644920, 69428442576136, 855504181649448, 10204459810035768, 118364711625485256, 1340006035830921720, 14850353930248138104, 161502853638370415864, 1727146533728893094604
Offset: 1

Views

Author

N. J. A. Sloane, Aug 01 2012

Keywords

Crossrefs

Programs

  • Mathematica
    DeleteCases[CoefficientList[Series[# (# - 1)^7*(5 #^9 - 60 #^8 + 675 #^7 - 2947 #^6 + 10005 #^5 - 20235 #^4 + 28297 #^3 - 23937 #^2 + 11418 # - 1781)/(2 (# - 2)^12*(# + 1)^9) &[(1 - Sqrt[1 - 8 x])/(4 x)], {x, 0, 23}], x], 0] (* Michael De Vlieger, Nov 26 2018 *)
  • PARI
    seq(N) = {
      my(x='x+O('x^(N+2)), y=(1-sqrt(1-8*x))/(4*x));
      Vec(y*(y - 1)^7*(5*y^9 - 60*y^8 + 675*y^7 - 2947*y^6 + 10005*y^5 - 20235*y^4 + 28297*y^3 - 23937*y^2 + 11418*y - 1781)/(2*(y - 2)^12*(y + 1)^9));
    };
    seq(18) \\ Gheorghe Coserea, Nov 12 2018

Formula

G.f.: y*(y - 1)^7*(5*y^9 - 60*y^8 + 675*y^7 - 2947*y^6 + 10005*y^5 - 20235*y^4 + 28297*y^3 - 23937*y^2 + 11418*y - 1781)/(2*(y - 2)^12*(y + 1)^9), where y = C(2*x), C being the g.f. for A000108. - Gheorghe Coserea, Nov 12 2018

Extensions

a(13)-a(14) by Noam Zeilberger, Sep 16 2018
More terms from Gheorghe Coserea, Nov 11 2018
Showing 1-10 of 25 results. Next