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

A001813 Quadruple factorial numbers: a(n) = (2n)!/n!.

Original entry on oeis.org

1, 2, 12, 120, 1680, 30240, 665280, 17297280, 518918400, 17643225600, 670442572800, 28158588057600, 1295295050649600, 64764752532480000, 3497296636753920000, 202843204931727360000, 12576278705767096320000, 830034394580628357120000, 58102407620643984998400000
Offset: 0

Views

Author

Keywords

Comments

Counts binary rooted trees (with out-degree <= 2), embedded in plane, with n labeled end nodes of degree 1. Unlabeled version gives Catalan numbers A000108.
Define a "downgrade" to be the permutation which places the items of a permutation in descending order. We are concerned with permutations that are identical to their downgrades. Only permutations of order 4n and 4n+1 can have this property; the number of permutations of length 4n having this property are equinumerous with those of length 4n+1. If a permutation p has this property then the reversal of this permutation also has it. a(n) = number of permutations of length 4n and 4n+1 that are identical to their downgrades. - Eugene McDonnell (eemcd(AT)mac.com), Oct 26 2003
Number of broadcast schemes in the complete graph on n+1 vertices, K_{n+1}. - Calin D. Morosan (cd_moros(AT)alumni.concordia.ca), Nov 28 2008
Hankel transform is A137565. - Paul Barry, Nov 25 2009
The e.g.f. of 1/a(n) = n!/(2*n)! is (exp(sqrt(x)) + exp(-sqrt(x)) )/2. - Wolfdieter Lang, Jan 09 2012
From Tom Copeland, Nov 15 2014: (Start)
Aerated with intervening zeros (1,0,2,0,12,0,120,...) = a(n) (cf. A123023 and A001147), the e.g.f. is e^(t^2), so this is the base for the Appell sequence with e.g.f. e^(t^2) e^(x*t) = exp(P(.,x),t) (reverse A059344, cf. A099174, A066325 also). P(n,x) = (a. + x)^n with (a.)^n = a_n and comprise the umbral compositional inverses for e^(-t^2)e^(x*t) = exp(UP(.,x),t), i.e., UP(n,P(.,t)) = x^n = P(n,UP(.,t)), e.g., (P(.,t))^n = P(n,t).
Equals A000407*2 with leading 1 added. (End)
a(n) is also the number of square roots of any permutation in S_{4*n} whose disjoint cycle decomposition consists of 2*n transpositions. - Luis Manuel Rivera Martínez, Mar 04 2015
Self-convolution gives A076729. - Vladimir Reshetnikov, Oct 11 2016
For n > 1, it follows from the formula dated Aug 07 2013 that a(n) is a Zumkeller number (A083207). - Ivan N. Ianakiev, Feb 28 2017
For n divisible by 4, a(n/4) is the number of ways to place n points on an n X n grid with pairwise distinct abscissae, pairwise distinct ordinates, and 90-degree rotational symmetry. For n == 1 (mod 4), the number of ways is a((n-1)/4) because the center point can be considered "fixed". For 180-degree rotational symmetry see A006882, for mirror symmetry see A000085, A135401, and A297708. - Manfred Scheucher, Dec 29 2017

Examples

			The following permutations of order 8 and their reversals have this property:
  1 7 3 5 2 4 0 6
  1 7 4 2 5 3 0 6
  2 3 7 6 1 0 4 5
  2 4 7 1 6 0 3 5
  3 2 6 7 0 1 5 4
  3 5 1 7 0 6 2 4
		

References

  • D. E. Knuth, The Art of Computer Programming, Vol. 4, Section 7.2.1.6, Eq. 32.
  • L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181.
  • Eugene McDonnell, "Magic Squares and Permutations" APL Quote-Quad 7.3 (Fall, 1976)
  • R. W. Robinson, Counting arrangements of bishops, pp. 198-214 of Combinatorial Mathematics IV (Adelaide 1975), Lect. Notes Math., 560 (1976).
  • 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

Programs

  • GAP
    List([0..20],n->Factorial(2*n)/Factorial(n)); # Muniru A Asiru, Nov 01 2018
    
  • Magma
    [Factorial(2*n)/Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 09 2018
    
  • Maple
    A001813 := n->(2*n)!/n!;
    A001813 := n -> mul(k, k = select(k-> k mod 4 = 2,[$1 .. 4*n])):
    seq(A001813(n), n=0..16);  # Peter Luschny, Jun 23 2011
  • Mathematica
    Table[(2n)!/n!, {n,0,20}] (* Harvey P. Dale, May 02 2011 *)
  • Maxima
    makelist(binomial(n+n, n)*n!,n,0,30); /* Martin Ettl, Nov 05 2012 */
    
  • PARI
    a(n)=binomial(n+n,n)*n! \\ Charles R Greathouse IV, Jun 15 2011
    
  • PARI
    first(n) = x='x+O('x^n); Vec(serlaplace((1 - 4*x)^(-1/2))) \\ Iain Fox, Jan 01 2018 (corrected by Iain Fox, Jan 11 2018)
    
  • Python
    from math import factorial
    def A001813(n): return factorial(n<<1)//factorial(n) # Chai Wah Wu, Feb 14 2023
  • Sage
    [binomial(2*n,n)*factorial(n) for n in range(0, 17)] # Zerinvary Lajos, Dec 03 2009
    

Formula

E.g.f.: (1-4*x)^(-1/2).
a(n) = (2*n)!/n! = Product_{k=0..n-1} (4*k + 2) = A081125(2*n).
Integral representation as n-th moment of a positive function on a positive half-axis: a(n) = Integral_{x=0..oo} x^n*exp(-x/4)/(sqrt(x)*2*sqrt(Pi)) dx, n >= 0. This representation is unique. - Karol A. Penson, Sep 18 2001
Define a'(1)=1, a'(n) = Sum_{k=1..n-1} a'(n-k)*a'(k)*C(n, k); then a(n)=a'(n+1). - Benoit Cloitre, Apr 27 2003
With interpolated zeros (1, 0, 2, 0, 12, ...) this has e.g.f. exp(x^2). - Paul Barry, May 09 2003
a(n) = A000680(n)/A000142(n)*A000079(n) = Product_{i=0..n-1} (4*i + 2) = 4^n*Pochhammer(1/2, n) = 4^n*GAMMA(n+1/2)/sqrt(Pi). - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
For asymptotics, see the Robinson paper.
a(k) = (2*k)!/k! = Sum_{i=1..k+1} |A008275(i,k+1)| * k^(i-1). - André F. Labossière, Jun 21 2007
a(n) = 12*A051618(a) n >= 2. - Zerinvary Lajos, Feb 15 2008
a(n) = A000984(n)*A000142(n). - Zerinvary Lajos, Mar 25 2008
a(n) = A016825(n-1)*a(n-1). - Roger L. Bagula, Sep 17 2008
a(n) = (-1)^n*A097388(n). - D. Morosan (cd_moros(AT)alumni.concordia.ca), Nov 28 2008
From Paul Barry, Jan 15 2009: (Start)
G.f.: 1/(1-2x/(1-4x/(1-6x/(1-8x/(1-10x/(1-... (continued fraction);
a(n) = (n+1)!*A000108(n). (End)
a(n) = Sum_{k=0..n} A132393(n,k)*2^(2n-k). - Philippe Deléham, Feb 10 2009
G.f.: 1/(1-2x-8x^2/(1-10x-48x^2/(1-18x-120x^2/(1-26x-224x^2/(1-34x-360x^2/(1-42x-528x^2/(1-... (continued fraction). - Paul Barry, Nov 25 2009
a(n) = A173333(2*n,n) for n>0; cf. A006963, A001761. - Reinhard Zumkeller, Feb 19 2010
From Gary W. Adamson, Jul 19 2011: (Start)
a(n) = upper left term of M^n, M = an infinite square production matrix as follows:
2, 2, 0, 0, 0, 0, ...
4, 4, 4, 0, 0, 0, ...
6, 6, 6, 6, 0, 0, ...
8, 8, 8, 8, 8, 0, ...
...
(End)
a(n) = (-2)^n*Sum_{k=0..n} 2^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: 1/Q(0), where Q(k) = 1 + x*(4*k+2) - x*(4*k+4)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 18 2013
G.f.: 2/G(0), where G(k) = 1 + 1/(1 - x*(8*k+4)/(x*(8*k+4) - 1 + 8*x*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - 2*x/(2*x + 1/(2*k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 01 2013
D-finite with recurrence: a(n) = (4*n-6)*a(n-2) + (4*n-3)*a(n-1), n>=2. - Ivan N. Ianakiev, Aug 07 2013
Sum_{n>=0} 1/a(n) = (exp(1/4)*sqrt(Pi)*erf(1/2) + 2)/2 = 1 + A214869, where erf(x) is the error function. - Ilya Gutkovskiy, Nov 10 2016
Sum_{n>=0} (-1)^n/a(n) = 1 - sqrt(Pi)*erfi(1/2)/(2*exp(1/4)), where erfi(x) is the imaginary error function. - Amiram Eldar, Feb 20 2021
a(n) = 1/([x^n] hypergeom([1], [1/2], x/4)). - Peter Luschny, Sep 13 2024
a(n) = 2^n*n!*JacobiP(n, -1/2, -n, 3). - Peter Luschny, Jan 22 2025
G.f.: 2F0(1,1/2;;4x). - R. J. Mathar, Jun 07 2025

Extensions

More terms from James Sellers, May 01 2000

A052701 a(0) = 0; for n >= 1, a(n) = 2^(n-1)*C(n-1), where C(n) = A000108(n) Catalan numbers, n>0.

Original entry on oeis.org

0, 1, 2, 8, 40, 224, 1344, 8448, 54912, 366080, 2489344, 17199104, 120393728, 852017152, 6085836800, 43818024960, 317680680960, 2317200261120, 16992801914880, 125210119372800, 926554883358720, 6882979133521920
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

A151374 shifted one place right. - Joerg Arndt, Mar 17 2011
The number of rooted Eulerian n-edge maps in the plane (planar with a distinguished outside face). - Valery A. Liskovets, Mar 17 2005
This is also the number of strings of length 2n-2 of two different types of balanced parentheses. For example, a(2) = 2, since the two possible strings of length 2 are [] and (), a(3) = 8, since the 8 possible strings of length 4 are (()), [()], ([]), [[]], ()(), [](), ()[], and [][]. - Jeffrey Shallit, Jun 03 2006
Row sums of number triangle A110506. - Paul Barry, Jul 24 2005
Also row sums of triangle in A085880. - Philippe Deléham, Aug 01 2005
Row sums of number triangle A114608. - Philippe Deléham, Oct 15 2008

References

  • V. A. Liskovets and T. R. Walsh, Enumeration of unrooted maps on the plane, Rapport technique, UQAM, No. 2005-01, Montreal, Canada, 2005.

Crossrefs

Limit of array A102544.

Programs

  • Maple
    spec := [S,{B=Union(C,Z),S=Union(B,C),C=Prod(S,S)},unlabeled]: seq(combstruct[count](spec,size=n), n=0..20);
  • Mathematica
    InverseSeries[Series[y-2*y^2, {y, 0, 24}], x] (* then A(x)=y(x) *) (* Len Smiley, Apr 07 2000 *)
    Join[{0},Table[2^n CatalanNumber[n],{n,0,30}]] (* Harvey P. Dale, Aug 29 2015 *)
  • PARI
    a(n)=if(n<1,0,2^(n-1)*(2*n-2)!/(n-1)!/n!)
    
  • PARI
    a(n)=if(n<1,0,polcoeff(serreverse(x-2*x^2+x*O(x^n)),n))
    
  • PARI
    a(n)=if(n<1,0,polcoeff(2*x/(1+sqrt(1-8*x+O(x^n))),n))

Formula

a(n) = A052714(n)/n!.
a(n) = A003645(n-2)*2, n>1.
a(n) = 8^(n-1)*GAMMA(n-1/2)/GAMMA(n+1)/Pi^(1/2), n>0.
D-finite with recurrence: a(1)=1, (-4+8*n)*a(n) - (n+1)*a(n+1) = 0.
G.f.: (1-sqrt(1-8*x))/4 = x*C(2*x) where C(x) is g.f. for Catalan numbers, A000108.
G.f. A(x) satisfies 2*A(x)^2-A(x)+x=0, A(0)=0 and A(x)=x+2*A(x)^2=x/(1-2*A(x)). Series reversion of x-2*x^2. - Michael Somos, Sep 06 2003
a(0)=0, a(1)=1; a(n) = 2*Sum_{i=1..n-1} a(i)*a(n-i). - Benoit Cloitre, Mar 16 2004
With a different offset, a(0)=1, a(n) = Sum_{k=0..n} Sum_{j=0..n} (j*C(2n-j-1, n-j)*C(j, k)*2^(n-j)/n), n>0. - Paul Barry, Jul 24 2005
The Hankel transform of a(n+1) = [1,2,8,40,224,1344,...] is 4^C(n+1,2). - Philippe Deléham, Nov 06 2007
G.f.: x + 4*x^2/(G(0)-4*x) where G(k) = k*(8*x+1) + 4*x + 2 - 2*x*(2*k+3)*(2*k+4)/G(k+1) ; (continued fraction ). - Sergei N. Gladkovskii, Apr 05 2013
a(n) ~ 8^(n-1)/(sqrt(Pi)*n^(3/2)). - Ilya Gutkovskiy, Dec 04 2016
From Amiram Eldar, Mar 06 2022: (Start)
Sum_{n>=1} 1/a(n) = 68/49 + 96*arcsin(sqrt(1/8))/(49*sqrt(7)).
Sum_{n>=1} (-1)^(n+1)/a(n) = 20/27 - 16*log(2)/81. (End)
a(n) = A025225(n)/2 for n>=1. - Alois P. Heinz, Feb 16 2025

Extensions

Better description from Claude Lenormand (claude.lenormand(AT)free.fr), Mar 19 2001
Additional comments from Michael Somos, Feb 24 2002

A090802 Triangle read by rows: a(n,k) = number of k-length walks in the Hasse diagram of a Boolean algebra of order n.

Original entry on oeis.org

1, 2, 1, 4, 4, 2, 8, 12, 12, 6, 16, 32, 48, 48, 24, 32, 80, 160, 240, 240, 120, 64, 192, 480, 960, 1440, 1440, 720, 128, 448, 1344, 3360, 6720, 10080, 10080, 5040, 256, 1024, 3584, 10752, 26880, 53760, 80640, 80640, 40320
Offset: 0

Views

Author

Ross La Haye, Feb 10 2004

Keywords

Comments

Row sums = A010842(n); Row sums from column 1 on = A066534(n) = n*A010842(n-1) = A010842(n) - 2^n.
a(n,k) = n! = k! = A000142(n) for n = k; a(n,n-1) = 2*n! = A052849(n) for n > 1; a(n,n-2) = 2*n! = A052849(n) for n > 2; a(n,n-3) = (4/3)*n! = A082569(n) for n > 3; a(n,n-1)/a(2,1) = n!/2! = A001710(n) for n > 1; a(n,n-2)/ a(3,1) = n!/3! = A001715(n) for n > 2; a(n,n-3)/a(4,1) = n!/4! = A001720(n) for n > 3.
a(2k, k) = A052714(k+1). a(2k-1, k) = A034910(k).
a(n,0) = A000079(n); a(n,1) = A001787(n) = row sums of A003506; a(n,2) = A001815(n) = 2!*A001788(n-1); a(n,3) = A052771(n) = 3!*A001789(n); a(n,4) = A052796(n) = 4!*A003472(n); ceiling[a(n,1) / 2] = A057711(n); a(n,5) = 5!*A054849(n).
In a class of n students, the number of committees (of any size) that contain an ordered k-sized subcommittee is a(n,k). - Ross La Haye, Apr 17 2006
Antidiagonal sums [1,2,5,12,30,76,198,528,1448,4080,...] appear to be binomial transform of A000522 interleaved with itself, i.e., 1,1,2,2,5,5,16,16,65,65,... - Ross La Haye, Sep 09 2006
Let P(A) be the power set of an n-element set A. Then a(n,k) = the number of ways to add k elements of A to each element x of P(A) where the k elements are not elements of x and order of addition is important. - Ross La Haye, Nov 19 2007
The derivatives of x^n evaluated at x=2. - T. D. Noe, Apr 21 2011

Examples

			{1};
{2, 1};
{4, 4, 2};
{8, 12, 12, 6};
{16, 32, 48, 48, 24};
{32, 80, 160, 240, 240, 120};
{64, 192, 480, 960, 1440, 1440, 720};
{128, 448, 1344, 3360, 6720, 10080, 10080, 5040};
{256, 1024, 3584, 10752, 26880, 53760, 80640, 80640, 40320}
a(5,3) = 240 because P(5,3) = 60, 2^(5-3) = 4 and 60 * 4 = 240.
		

Crossrefs

Programs

  • Mathematica
    Flatten[Table[n!/(n-k)! * 2^(n-k), {n, 0, 8}, {k, 0, n}]] (* Ross La Haye, Feb 10 2004 *)

Formula

a(n, k) = 0 for n < k. a(n, k) = k!*C(n, k)*2^(n-k) = P(n, k)*2^(n-k) = (2n)!!/((n-k)!*2^k) = k!*A038207(n, k) = A068424*2^(n-k) = Sum[C(n, m)*P(n-m, k), {m, 0, n-k}] = Sum[C(n, n-m)*P(n-m, k), {m, 0, n-k}] = n!*Sum[1/(m!*(n-m-k)!), {m, 0, n-k}] = k!*Sum[C(n, m)*C(n-m, k), {m, 0, n-k}] = k!*Sum[C(n, n-m)*C(n-m, k), {m, 0, n-k}] = k!*C(n, k)*Sum[C(n-k, n-m-k), {m, 0, n-k}] = k!*C(n, k)*Sum[C(n-k, m), {m, 0, n-k}] for n >= k.
a(n, k) = 0 for n < k. a(n, k) = n*a(n-1, k-1) for n >= k >= 1.
E.g.f. (by columns): exp(2x)*x^k.

Extensions

More terms from Ray Chandler, Feb 26 2004
Entry revised by Ross La Haye, Aug 18 2006

A144828 Partial products of successive terms of A017113; a(0)=1.

Original entry on oeis.org

1, 4, 48, 960, 26880, 967680, 42577920, 2214051840, 132843110400, 9033331507200, 686533194547200, 57668788341964800, 5305528527460761600, 530552852746076160000, 57299708096576225280000, 6646766139202842132480000, 824199001261152424427520000, 108794268166472120024432640000
Offset: 0

Views

Author

Philippe Deléham, Sep 21 2008

Keywords

Comments

a(n) is the number of signed permutations of length 4n that are equal to their reverse-inverses. Note that the reverse-inverse of a permutation is equivalent to a 90-degree rotation of the permutation's diagram (see the Hardt and Troyka link). - Justin M. Troyka, Aug 11 2011
Define the bar operation as an operation on signed permutation that flips the sign of each entry. Then a(n) is the number of signed permutations of length 2n that are equal to the bar of their inverses and equal to their reverse-complements (see the Hardt and Troyka link). - Justin M. Troyka, Aug 11 2011

Examples

			a(0)=1, a(1)=4, a(2)=4*12=48, a(3)=4*12*20=960, a(4)=4*12*20*28=26880, ...
Since a(1) = 4, there are 4 signed permutations of 4 that are equal to their reverse-inverses.  These are: (+2,+4,+1,+3), (+3,+1,+4,+2), (-2,-4,-1,-3), (-3,-1,-4,-2). - _Justin M. Troyka_, Aug 11 2011
G.f. = 1 + 4*x + 48*x^2 + 960*x^3 + 26880*x^4 + 967680*x^5 + 42577920*x^6 + ...
		

Crossrefs

Essentially the same as A052714. - N. J. A. Sloane, Feb 03 2013
Sequences of the form m^(n-1)*n!*Catalan(n-1): A001813 (m=1), A052714 (or this sequence) (m=2), A221954 (m=3), A052734 (m=4), A221953 (m=5), A221955 (m=6).

Programs

  • Magma
    [2^k *Factorial(2*k) / Factorial(k): k in [0..20]]; // Vincenzo Librandi, Aug 11 2011
    
  • Maple
    A144828:= n-> 2^n*n!*binomial(2*n,n); seq(A144828(n), n=0..30); # G. C. Greubel, Apr 02 2021
  • Mathematica
    Table[4^n (2 n - 1)!!, {n, 0, 15}] (* Vincenzo Librandi, May 14 2015 *)
    Join[{1},FoldList[Times,(8*Range[0,20]+4)]] (* Harvey P. Dale, Dec 01 2015 *)
  • PARI
    a(n)=binomial(2*n,n)*n!<Charles R Greathouse IV, Jan 17 2012
    
  • PARI
    {a(n) = if( n<0, (-1)^n / a(-n), 2^n *(2*n)! / n!)}; /* Michael Somos, Jan 06 2017 */
    
  • Sage
    [2^n*factorial(n+1)*catalan_number(n) for n in (0..30)] # G. C. Greubel, Apr 02 2021

Formula

a(n) = Sum_{k=0..n} A132393(n,k)*4^k*8^(n-k).
a(n) = A052714(n+1). - R. J. Mathar, Oct 01 2008
a(n) = 2^n *(2*n)! / n!. - Justin M. Troyka, Aug 11 2011
G.f.: 1/(1-4x/(1-8x/(1-12x/(1-16x/(1-20x/(1-24x/(1-28x/(1-32x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2012
a(n) = (-4)^n*Sum_{k=0..n} 2^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
E.g.f.: 1/sqrt(1-8*x). - Philippe Deléham, May 14 2015
a(n) = 4^n * A001147(n). - Philippe Deléham, May 14 2015
a(n) = 8^n * Gamma(n + 1/2) / sqrt(Pi). - Daniel Suteu, Jan 06 2017
0 = a(n)*(8*a(n+1) - a(n+2)) + a(n+1)*(+a(n+1)) and a(n) = (-1)^n / a(-n) for all n in Z. - Michael Somos, Jan 06 2017
a(n) = 2^n * (n+1)! * Catalan(n). - G. C. Greubel, Apr 02 2021
Sum_{n>=0} 1/a(n) = 1 + e^(1/8)*sqrt(Pi)*erf(1/(2*sqrt(2)))/(2*sqrt(2)), where erf is the error function. - Amiram Eldar, Dec 20 2022

A052712 Expansion of e.g.f. (1 + 4*x - sqrt(1-8*x))/8.

Original entry on oeis.org

0, 1, 2, 24, 480, 13440, 483840, 21288960, 1107025920, 66421555200, 4516665753600, 343266597273600, 28834394170982400, 2652764263730380800, 265276426373038080000, 28649854048288112640000
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

Also the number of random walk labelings of the 2 X (n-1) king's graph, for n > 1. - Sela Fried, Apr 14 2023

Crossrefs

Programs

  • Maple
    spec := [S,{B=Prod(C,C),C=Union(B,S),S=Union(B,Z)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
  • Mathematica
    Table[n!*2^(n-2)*CatalanNumber[n-1] +Boole[n==1]/2 +Boole[n==0]/4, {n,0,30}] (* G. C. Greubel, May 30 2022 *)
  • SageMath
    [2^(n-2)*factorial(n)*catalan_number(n-1) +bool(n==0)/8 +bool(n==1)/2 for n in (0..30)] # G. C. Greubel, May 30 2022

Formula

D-finite with recurrence: a(0) = 0, a(1)=1, a(2)=2, a(n+1) = 4*(2*n-1)*a(n).
a(n) = 8^(n+1)*Gamma(n+3/2)/sqrt(Pi).
a(n) = n!*A003645(n-2), n>1. - R. J. Mathar, Oct 18 2013
G.f.: (1 + 4*x - 2F0([1,-1/2], [], 8*x))/8. - R. J. Mathar, Jan 25 2020

A052723 Expansion of e.g.f. (1 - x - sqrt(1-2*x+x^2-4*x^3))/(2*x).

Original entry on oeis.org

0, 0, 2, 6, 24, 240, 2880, 35280, 524160, 9434880, 188697600, 4151347200, 101548339200, 2727435110400, 79332244992000, 2488504322304000, 83879464660992000, 3021209014247424000, 115754916599562240000
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Crossrefs

Programs

  • Maple
    spec := [S,{B=Prod(S,S),C=Union(B,S,Z),S=Prod(Z,C)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
    seq(n!*add(binomial(n-2-k,2*k)*binomial(2*k,k)/(k+1), k=0..floor((n-2)/3)), n=0..18);  # Mark van Hoeij, May 12 2013
  • Mathematica
    With[{nn=20},CoefficientList[Series[(1-x-Sqrt[1-2x+x^2-4x^3])/(2x),{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Jan 19 2017 *)
    a[n_]:= a[n]= n!*Sum[Binomial[n-k-2,2*k]*CatalanNumber[k], {k,0,Floor[(n-2)/2]}];
    Table[a[n], {n,0,30}] (* G. C. Greubel, May 28 2022 *)
  • SageMath
    def A052723(n): return factorial(n)*sum( binomial(n-k-2, 2*k)*catalan_number(k) for k in (0..(n-2)//2) )
    [A052723(n) for n in (0..30)] # G. C. Greubel, May 28 2022

Formula

D-finite with recurrence: a(0) = a(1) = 0, a(2) = 2, a(3) = 6, a(4) = 24, (n+4)*a(n+3) = (15 + 11*n + 2*n^2)*a(n+2) - (6 + 11*n + 6*n^2 + n^3)*a(n+1) - (12 - 2*n - 32*n^2 - 22*n^2 - 4*n^4)*a(n).
a(n) = n!*A023431(n-2). - R. J. Mathar, Oct 18 2013

A052711 Expansion of e.g.f. x*(1 - 2*x - sqrt(1-4*x))/2.

Original entry on oeis.org

0, 0, 0, 6, 48, 600, 10080, 211680, 5322240, 155675520, 5189184000, 194075481600, 8045310873600, 366061644748800, 18134130709094400, 971471287987200000, 55956746188062720000, 3448334483839365120000
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Crossrefs

Programs

  • Maple
    spec := [S,{C=Union(B,Z),B=Prod(C,C),S=Prod(B,Z)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
  • Mathematica
    With[{nn=20},CoefficientList[Series[x (1-2x-Sqrt[1-4x])/2,{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Mar 05 2016 *)
    Table[n!*CatalanNumber[n-2] +Boole[n==1] -2*Boole[n==2], {n,0,30}] (* G. C. Greubel, May 30 2022 *)
  • SageMath
    [factorial(n)*catalan_number(n-2) + bool(n==1)/2 - 2*bool(n==2) for n in (0..30)] # G. C. Greubel, May 30 2022

Formula

D-finite with recurrence: a(1)=0, a(2)=0, a(3)=6, a(4)=48, n*a(n+1) = 2*(n+1)*(2*n-3)*a(n).
From R. J. Mathar, Oct 18 2013: (Start)
a(n) = n!*A000108(n-2).
a(n) = A052717(n), n>2. (End)
G.f.: x*(1 - 4*x - 2F0([-1/2,2], [], 4*x))/2. - R. J. Mathar, Jan 25 2020

A052713 Expansion of e.g.f. (1-sqrt(1-8*x))/2.

Original entry on oeis.org

0, 2, 8, 96, 1920, 53760, 1935360, 85155840, 4428103680, 265686220800, 18066663014400, 1373066389094400, 115337576683929600, 10611057054921523200, 1061105705492152320000, 114599416193152450560000
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

Has a square root singularity.

Crossrefs

Programs

  • Maple
    spec := [S,{C=Union(B,Z),B=Prod(S,S),S=Union(Z,C)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
  • Mathematica
    Table[n!*2^n*CatalanNumber[n-1] + Boole[n==0], {n,0,30}] (* G. C. Greubel, May 29 2022 *)
  • SageMath
    [2^n*factorial(n)*catalan_number(n-1) + bool(n==0)/2 for n in (0..30)] # G. C. Greubel, May 29 2022

Formula

D-finite with recurrence: a(1)=2, a(n+1) = 4*(2*n -1)*a(n).
a(n+1) = 1/4*8^(n+1)*Gamma(n+1/2)/Pi^(1/2)
a(n+1) = ((2*n)!/n!)*2^(n+1). - Zerinvary Lajos, Sep 25 2006
a(n) = n!*A025225(n). - R. J. Mathar, Oct 18 2013
G.f.: (1- 2F0([1,-1/2], [], 8*x))/2. - R. J. Mathar, Jan 25 2020

A052717 Expansion of e.g.f. x*(1 - sqrt(1 - 4*x))/2.

Original entry on oeis.org

0, 0, 2, 6, 48, 600, 10080, 211680, 5322240, 155675520, 5189184000, 194075481600, 8045310873600, 366061644748800, 18134130709094400, 971471287987200000, 55956746188062720000, 3448334483839365120000
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Crossrefs

Programs

  • Magma
    [n le 1 select 0 else Factorial(n)*Catalan(n-2): n in [0..30]]; // G. C. Greubel, May 28 2022
    
  • Maple
    spec := [S,{C=Union(B,Z),B=Prod(C,C),S=Prod(Z,C)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
    with(combinat):with(combstruct):a[0]:=0:for n from 1 to 30 do a[n]:=sum((count(Permutation(n*2-2),size=n-1)),j=0..n) od: seq(a[n], n=0..22); # Zerinvary Lajos, May 03 2007
  • Mathematica
    With[{nn=20},CoefficientList[Series[x (1-Sqrt[1-4x])/2,{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Dec 20 2015 *)
    Table[Boole[n==1] + n!*CatalanNumber[n-2], {n, 0, 30}] (* G. C. Greubel, May 28 2022 *)
  • MuPAD
    combinat::catalan(n)*(n+2)! $ n = 0..15; // Zerinvary Lajos, Feb 15 2007
    
  • SageMath
    [bool(n==1)/2 + factorial(n)*catalan_number(n-2) for n in (0..30)] # G. C. Greubel, May 28 2022

Formula

Recurrence: a(1)=0, a(3)=6, a(2)=2, n*a(n+1) = (4*n^2 - 2*n - 6)*a(n).
a(n) = n!*A000108(n-2) = A052711(n), n > 2. - R. J. Mathar, Oct 26 2013
G.f.: x*(d/dx)(x^2 * Hypergeometric2F0([1, 1/2], [], 4*x)). - G. C. Greubel, May 28 2022

A052719 Expansion of e.g.f. (1 - 2*x*sqrt(1-4*x))*(1 - sqrt(1-4*x))/4.

Original entry on oeis.org

0, 0, 0, 6, 72, 1080, 20160, 453600, 11975040, 363242880, 12454041600, 476367091200, 20113277184000, 929233405900800, 46630621823385600, 2525825348766720000, 146886458743664640000, 9127944221927731200000
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Crossrefs

Programs

  • Maple
    spec := [S,{B=Union(Z,C),C=Prod(B,B),S=Prod(B,C)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
  • Mathematica
    Table[If[n<2, 0, 3*(n-2)*(n-1)!*CatalanNumber[n-2]], {n,0,30}] (* G. C. Greubel, May 28 2022 *)
  • SageMath
    [0,0]+[3*(n-2)*factorial(n-1)*catalan_number(n-2) for n in (2..30)] # G. C. Greubel, May 28 2022

Formula

D-finite with recurrence: a(1)=0, a(2)=0, a(3)=6, a(n+2) = (2 + 5*n)*a(n+1) + (6 + 2*n - 4*n^2)*a(n)
a(n) = n!*A000245(n-2). - R. J. Mathar, Oct 26 2013
From G. C. Greubel, May 28 2022: (Start)
G.f.: 6*x^3*Hypergeometric2F0([2, 3/2], [], 4*x).
E.g.f.: (1/4)*(1 + 2*x - 8*x^2 - (1 + 2*x)*sqrt(1-4*x)). (End)
Showing 1-10 of 21 results. Next