A000900 -id:A000900 - cp's OEIS frontend

A000085 Number of self-inverse permutations on n letters, also known as involutions; number of standard Young tableaux with n cells.

1, 1, 2, 4, 10, 26, 76, 232, 764, 2620, 9496, 35696, 140152, 568504, 2390480, 10349536, 46206736, 211799312, 997313824, 4809701440, 23758664096, 119952692896, 618884638912, 3257843882624, 17492190577600, 95680443760576, 532985208200576, 3020676745975552

Offset: 0

N. J. A. Sloane and J. H. Conway

a(n) is also the number of n X n symmetric permutation matrices.
a(n) is also the number of matchings (Hosoya index) in the complete graph K(n). - Ola Veshta (olaveshta(AT)my-deja.com), Mar 25 2001
a(n) is also the number of independent vertex sets and vertex covers in the n-triangular graph. - Eric W. Weisstein, May 22 2017
Equivalently, this is the number of graphs on n labeled nodes with degrees at most 1. - Don Knuth, Mar 31 2008
a(n) is also the sum of the degrees of the irreducible representations of the symmetric group S_n. - Avi Peretz (njk(AT)netvision.net.il), Apr 01 2001
a(n) is the number of partitions of a set of n distinguishable elements into sets of size 1 and 2. - Karol A. Penson, Apr 22 2003
Number of tableaux on the edges of the star graph of order n, S_n (sometimes T_n). - Roberto E. Martinez II, Jan 09 2002
The Hankel transform of this sequence is A000178 (superfactorials). Sequence is also binomial transform of the sequence 1, 0, 1, 0, 3, 0, 15, 0, 105, 0, 945, ... (A001147 with interpolated zeros). - Philippe Deléham, Jun 10 2005
Row sums of the exponential Riordan array (e^(x^2/2),x). - Paul Barry, Jan 12 2006
a(n) is the number of nonnegative lattice paths of upsteps U = (1,1) and downsteps D = (1,-1) that start at the origin and end on the vertical line x = n in which each downstep (if any) is marked with an integer between 1 and the height of its initial vertex above the x-axis. For example, with the required integer immediately preceding each downstep, a(3) = 4 counts UUU, UU1D, UU2D, U1DU. - David Callan, Mar 07 2006
Equals row sums of triangle A152736 starting with offset 1. - Gary W. Adamson, Dec 12 2008
Proof of the recurrence relation a(n) = a(n-1) + (n-1)*a(n-2): number of involutions of [n] containing n as a fixed point is a(n-1); number of involutions of [n] containing n in some cycle (j, n), where 1 <= j <= n-1, is (n-1) times the number of involutions of [n] containing the cycle (n-1 n) = (n-1)*a(n-2). - Emeric Deutsch, Jun 08 2009
Number of ballot sequences (or lattice permutations) of length n. A ballot sequence B is a string such that, for all prefixes P of B, h(i) >= h(j) for i < j, where h(x) is the number of times x appears in P. For example, the ballot sequences of length 4 are 1111, 1112, 1121, 1122, 1123, 1211, 1212, 1213, 1231, and 1234. The string 1221 does not appear in the list because in the 3-prefix 122 there are two 2's but only one 1. (Cf. p. 176 of Bruce E. Sagan: "The Symmetric Group"). - Joerg Arndt, Jun 28 2009
Number of standard Young tableaux of size n; the ballot sequences are obtained as a length-n vector v where v_k is the (number of the) row in which the number r occurs in the tableaux. - Joerg Arndt, Jul 29 2012
Number of factorial numbers of length n-1 with no adjacent nonzero digits. For example the 10 such numbers (in rising factorial radix) of length 3 are 000, 001, 002, 003, 010, 020, 100, 101, 102, and 103. - Joerg Arndt, Nov 11 2012
Also called restricted Stirling numbers of the second kind (see Mezo). - N. J. A. Sloane, Nov 27 2013
a(n) is the number of permutations of [n] that avoid the consecutive patterns 123 and 132. Proof. Write a self-inverse permutation in standard cycle form: smallest entry in each cycle in first position, first entries decreasing. For example, (6,7)(3,4)(2)(1,5) is in standard cycle form. Then erase parentheses. This is a bijection to the permutations that avoid consecutive 123 and 132 patterns. - David Callan, Aug 27 2014
Getu (1991) says these numbers are also known as "telephone numbers". - N. J. A. Sloane, Nov 23 2015
a(n) is the number of elements x in the symmetric group S_n such that x^2 = e where e is the identity. - Jianing Song, Aug 22 2018 [Edited on Jul 24 2025]
a(n) is the number of congruence orbits of upper-triangular n X n matrices on skew-symmetric matrices, or the number of Borel orbits in largest sect of the type DIII symmetric space SO_{2n}(C)/GL_n(C). Involutions can also be thought of as fixed-point-free partial involutions. See [Bingham and Ugurlu] link. - Aram Bingham, Feb 08 2020
From Thomas Anton, Apr 20 2020: (Start)
Apparently a(n) = b*c where b is odd iff a(n+b) (when a(n) is defined) is divisible by b.
Apparently a(n) = 2^(f(n mod 4)+floor(n/4))*q where f:{0,1,2,3}->{0,1,2} is given by f(0),f(1)=0, f(2)=1 and f(3)=2 and q is odd. (End)
From Iosif Pinelis, Mar 12 2021: (Start)
a(n) is the n-th initial moment of the normal distribution with mean 1 and variance 1. This follows because the moment generating function of that distribution is the e.g.f. of the sequence of the a(n)'s.
The recurrence a(n) = a(n-1) + (n-1)*a(n-2) also follows, by writing E(Z+1)^n=EZ(Z+1)^(n-1)+E(Z+1)^(n-1), where Z is a standard normal random variable, and then taking the first of the latter two integrals by parts. (End)

			Sequence starts 1, 1, 2, 4, 10, ... because possibilities are {}, {A}, {AB, BA}, {ABC, ACB, BAC, CBA}, {ABCD, ABDC, ACBD, ADCB, BACD, BADC, CBAD, CDAB, DBCA, DCBA}. - _Henry Bottomley_, Jan 16 2001
G.f. = 1 + x + 2*x^2 + 4*x^4 + 10*x^5 + 26*x^6 + 76*x^7 + 232*x^8 + 764*x^9 + ...
From _Gus Wiseman_, Jan 08 2021: (Start)
The a(4) = 10 standard Young tableaux:
  1 2 3 4
.
  1 2   1 3   1 2 3   1 2 4   1 3 4
  3 4   2 4   4       3       2
.
  1 2   1 3   1 4
  3     2     2
  4     4     3
.
  1
  2
  3
  4
The a(0) = 1 through a(4) = 10 set partitions into singletons or pairs:
  {}  {{1}}  {{1,2}}    {{1},{2,3}}    {{1,2},{3,4}}
             {{1},{2}}  {{1,2},{3}}    {{1,3},{2,4}}
                        {{1,3},{2}}    {{1,4},{2,3}}
                        {{1},{2},{3}}  {{1},{2},{3,4}}
                                       {{1},{2,3},{4}}
                                       {{1,2},{3},{4}}
                                       {{1},{2,4},{3}}
                                       {{1,3},{2},{4}}
                                       {{1,4},{2},{3}}
                                       {{1},{2},{3},{4}}
(End)

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, pages 32, 911.
S. Chowla, The asymptotic behavior of solutions of difference equations, in Proceedings of the International Congress of Mathematicians (Cambridge, MA, 1950), Vol. I, 377, Amer. Math. Soc., Providence, RI, 1952.
W. Fulton, Young Tableaux, Cambridge, 1997.
D. E. Knuth, The Art of Computer Programming, Vol. 3, Section 5.1.4, p. 65.
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181.
T. Muir, A Treatise on the Theory of Determinants. Dover, NY, 1960, p. 6.
J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 86.
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).
R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.10.

Alois P. Heinz, Table of n, a(n) for n = 0..800 (first 101 terms from T. D. Noe)
Ahmad Zainy Al-Yasry, Short Steps in Noncommutative Geometry, arXiv:1901.03640 [math.OA], 2019.
T. Amdeberhan and V. H. Moll, Involutions and their progenies, 2014.
David Applegate and N. J. A. Sloane, The Gift Exchange Problem arXiv:0907.0513 [math.CO], 2009.
Joerg Arndt, Generating Random Permutations, PhD thesis, Australian National University, Canberra, Australia, (2010).
Joerg Arndt, Matters Computational (The Fxtbook), pp. 279-280.
C. Banderier, M. Bousquet-Mélou, A. Denise, P. Flajolet, D. Gardy and D. Gouyou-Beauchamps, Generating Functions for Generating Trees, Discrete Mathematics 246(1-3), March 2002, pp. 29-55.
D. Barsky, Analyse p-adique et suites classiques de nombres, Sem. Loth. Comb. B05b (1981) 1-21.
C. Bebeacua, T. Mansour, A. Postnikov and S. Severini, On the X-rays of permutations, arXiv:math/0506334 [math.CO], 2005.
E. A. Bender and S. G. Williamson, Foundations of Combinatorics with Applications (see Chap. 2, Example 2.9, pp. 47-48, including Theorem 2.2, a derived formula for a(n)). [From _Rick L. Shepherd_, Sep 02 2009]
S. Bilotta, E. Pergola, R. Pinzani, and S. Rinaldi, Recurrence relations versus succession rules, arXiv preprint arXiv:1301.2967 [cs.DM], 2013. - From _N. J. A. Sloane_, Feb 12 2013
S. Bilotta, E. Pergola, R. Pinzani, and S. Rinaldi, Recurrence Relations, Succession Rules, and the Positivity Problem, in Language and Automata Theory and Applications, 9th International Conference, LATA 2015, Nice, France, March 2-6, 2015, Proceedings, Pages 499-510, Lecture Notes Comp. Sci. Vol. 8977.
Aram Bingham and Özlem Uğurlu, Sects, rooks, pyramids, partitions and paths for type DIII clans, arXiv:1907.08875 [math.CO], 2019.
Anders Björner and Richard P. Stanley, A combinatorial miscellany, 2010.
P. Blasiak, K. A. Penson, A. I. Solomon, A. Horzela and G. E. H. Duchamp, Combinatorial field theories via boson normal ordering, arXiv:quant-ph/0405103, 2004-2006. The title of this paper in the arXiv was later changed to "Some useful combinatorial formulas for bosonic operators"
P. Blasiak, K. A. Penson, A. I. Solomon, A. Horzela and G. E. H. Duchamp, Some useful combinatorial formulas for bosonic operators, J. Math. Phys. 46, 052110 (2005) (6 pages).
Tobias Boege and Thomas Kahle, Construction Methods for Gaussoids, arXiv:1902.11260 [math.CO], 2019.
D. Bundala, M. Codish, L. Cruz-Filipe et al., Optimal-Depth Sorting Networks, arXiv preprint arXiv:1412.5302 [cs.DS], 2014.
Alexander Burstein and Louis W. Shapiro, Pseudo-involutions in the Riordan group, arXiv:2112.11595 [math.CO], 2021.
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
D. Castellanos, A generalization of Binet's formula and some of its consequences, Fib. Quart., 27 (1989), 424-438.
C.-O. Chow, Counting involutory, unimodal and alternating signed permutations, Discr. Math. 306 (2006), 2222-2228.
S. Chowla, and I.N. Hernstein, W.K. Moore, On recursions connected with symmetric groups. I, Canadian Journal of Mathematics, 3 (1951), 328-334.
M. Codish, L. Cruz-Filipe and P. Schneider-Kamp, The Quest for Optimal Sorting Networks: Efficient Generation of Two-Layer Prefixes, arXiv preprint arXiv:1404.0948 [cs.DS], 2014.
Tom Copeland, Infinitesimal Generators, the Pascal Pyramid, and the Witt and Virasoro Algebras
Filippo Disanto and Thomas Wiehe, Some instances of a sub-permutation problem on pattern avoiding permutations, arXiv preprint arXiv:1210.6908 [math.CO], 2012.
I. Dolinka, J. East, A. Evangelou, D. FitzGerald, N. Ham, et al., Enumeration of idempotents in diagram semigroups and algebras, arXiv preprint arXiv:1408.2021 [math.GR], 2014.
I. Dolinka, J. East, and R. D. Gray, Motzkin monoids and partial Brauer monoids, arXiv preprint arXiv:1512.02279 [math.GR], 2015-2016.
R. Donaghey, Binomial self-inverse sequences and tangent coefficients, J. Combin. Theory, Series A, 21 (1976), 155-163.
S. Dulucq and J.-G. Penaud, Cordes, arbres et permutations, Discrete Math. 117 (1993), no. 1-3, 89-105.
Shalosh B. Ekhad, Manuel Kauers, and Doron Zeilberger, The O(1/n^85) Asymptotic expansion of OEIS sequence A85, arxiv preprint arXiv:2410.16334 [math.CO], October 2024.
Shalosh B. Ekhad and Doron Zeilberger, Computational and Theoretical Challenges on Counting Solid Standard Young Tableaux
Shalosh B. Ekhad and Doron Zeilberger, Computational and Theoretical Challenges on Counting Solid Standard Young Tableaux, [Local copy, pdf file only, no active links]
Shalosh B. Ekhad and Doron Zeilberger, Computational and Theoretical Challenges on Counting Solid Standard Young Tableaux, arXiv preprint arXiv:1202.6229 [math.CO], 2012. [This is a different document from the one with the same title on Doron Zeilberger's web site]
Shalosh B. Ekhad, Manuel Kauers, and Doron Zeilberger, The O(1/n^85) Asymptotic expansion of OEIS sequence A85, arXiv:2410.16334 [math.CO], 2024.
FindStat - Combinatorial Statistic Finder, Standard Young tableaux
Amy M. Fu and Frank Z. K. Li, Joint Distributions of Permutation Statistics and the Parabolic Cylinder Functions, arXiv:1809.07465 [math.CO], 2018.
Mohammad Ganjtabesh, Armin Morabbi and Jean-Marc Steyaert, Enumerating the number of RNA structures
S. Garrabrant and I. Pak, Counting with irrational tiles, arXiv:1407.8222 [math.CO], 2014.
S. Getu, Evaluating determinants via generating functions, Mathematics Magazine 64 (1991), 45-53.
J. Gilder, Rooks inviolate revisited II, Math. Gaz., 70 (1986), 44-45.
A. M. Goyt, Avoidance of partitions of a 3-element set, arXiv:math/0603481 [math.CO], 2006.
Greenhill, Catherine; McKay, Brendan D., Counting loopy graphs with given degree, Linear Algebra Appl. 436, No. 4, 901-926 (2012)
H. Gupta, Enumeration of symmetric matrices, Duke Math. J., 35 (1968), vol 3, 653-659
H. Gupta, Enumeration of symmetric matrices (annotated scanned copy)
C. Gutschwager, Reduced Kronecker products which are multiplicity free or contain only a few components, Eur. J. Combinat. 31 (2010) 1996-2005 doi:10.1016/j.ejc.2010.05.008, Remark 2.13.
T. Halverson and M. Reeks, Gelfand Models for Diagram Algebras, arXiv preprint arXiv:1302.6150 [math.RT], 2013.
Aoife Hennessy, A Study of Riordan Arrays with Applications to Continued Fractions, Orthogonal Polynomials and Lattice Paths, Ph. D. Thesis, Waterford Institute of Technology, Oct. 2011.
A. Horzela, P. Blasiak, G. E. H. Duchamp, K. A. Penson and A. I. Solomon, A product formula and combinatorial field theory, arXiv:quant-ph/0409152, 2004.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 17
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 22
András Kaszanyitzky, The GraftalLace Cellular Automaton, arXiv:1805.11532 [nlin.CG], 2018.
Manuel Kauers, Doron Zeilberger, Counting Standard Young Tableaux With Restricted Runs, arXiv:2006.10205 [math.CO], 2020.
Dongsu Kim and Jang Soo Kim, A combinatorial approach to the power of 2 in the number of involutions, J. Comb. Theory Ser. A 117 (8) (2010): 1082-1094
Vaclav Kotesovec, Non-attacking chess pieces, 6ed, 2013, p. 227.
P. L. Krapivsky and J. M. Luck, Coverage fluctuations in theater models, arXiv:1902.04365 [cond-mat.stat-mech], 2019.
M. B. Kutler and C. R. Vinroot, On q-Analogs of Recursions for the Number of Involutions and Prime Order Elements in Symmetric Groups, JIS 13 (2010) #10.3.6.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181. [Annotated scan of pages 180 and 181 only]
J. W. Layman, The Hankel Transform and Some of its Properties, J. Integer Sequences, 4 (2001), #01.1.5.
Steven Linton, James Propp, Tom Roby, and Julian West, Equivalence Classes of Permutations under Various Relations Generated by Constrained Transpositions, Journal of Integer Sequences, Vol. 15 (2012), #12.9.1.
E. Lucas, Théorie des Nombres, Gauthier-Villars, Paris, 1891, Vol. 1, p. 221.
E. Lucas, Théorie des nombres (annotated scans of a few selected pages)
T. Mansour and M. Shattuck, Partial matchings and pattern avoidance, Appl. Anal. Discrete Math. (to appear, 2012); - From _N. J. A. Sloane_, Jan 04 2013
I. Mezo, Periodicity of the last digits of some combinatorial sequences, arXiv preprint arXiv:1308.1637 [math.CO], 2013 and J. Int. Seq. 17 (2014) #14.1.1.
F. L. Miksa, L. Moser and M. Wyman, Restricted partitions of finite sets, Canad. Math. Bull., 1 (1958), 87-96.
L. Moser and M. Wyman, On solutions of x^d = 1 in symmetric groups, Canad. J. Math., 7 (1955), 159-168.
T. Mueller, Finite group actions and asymptotic expansion of e^P(z), Combinatorica, 17 (4) (1997), 523-554.
Igor Pak, Complexity problems in enumerative combinatorics, arXiv:1803.06636 [math.CO], 2018.
I. Pak, G. Panova, and E. Vallejo, Kronecker products, characters, partitions, and the tensor square conjectures, arXiv:1304.0738 [math.CO], 2013.
D. Panyushev, On the orbits of a Borel subgroup in abelian ideals, arXiv preprint arXiv:1407.6857 [math.AG], 2014.
P. Peart and W.-J. Woan, Generating Functions via Hankel and Stieltjes Matrices, J. Integer Seqs., Vol. 3 (2000), #00.2.1.
R. Petuchovas, Asymptotic analysis of the cyclic structure of permutations, arXiv:1611.02934 [math.CO], p. 6, 2016.
R. A. Proctor, Let's Expand Rota's Twelvefold Way for Counting Partitions!, arXiv:math/0606404 [math.CO], 2006-2007.
J. L. Ramírez and G. N. Rubiano, Properties and Generalizations of the Fibonacci Word Fractal, The Mathematica Journal, Vol. 16 (2014).
J. Riordan, Letter to N. J. A. Sloane, Feb 03 1975 (with notes by njas)
R. W. Robinson, Counting arrangements of bishops, Lect. Notes Math. 560 (1976), 198-214. See D_n.
H. A. Rothe, Ueber Permutationen, in Beziehung auf die Stellen ihrer Elemente. Anwendung der daraus abgeleiteten Sätze auf das Eliminationsproblem, in Carl Hindenberg, ed., Sammlung Combinatorisch-Analytischer Abhandlungen, pp. 263-305, Bey G. Fleischer dem jüngern, 1800; see p. 282.
Ahmad Sabri and Vincent Vajnovszki, Exhaustive generation for ballot sequences in lexicographic and Gray code order, CEUR Workshop Proceedings, 2018.
Ahmad Sabri and Vincent Vajnovszki, On the exhaustive generation of generalized ballot sequences in lexicographic and Gray code order, Pure Mathematics and Applications (2019) Vol. 28, Issue 1, 109-119.
A. I. Solomon, P. Blasiak, G. Duchamp, A. Horzela and K. A. Penson, Combinatorial physics, normal order and model Feynman graphs, arXiv:quant-ph/0310174, 2003.
Bridget Eileen Tenner, Star factorizations and noncrossing partitions, arXiv:2004.03286 [math.CO], 2020.
Robert F. Tichy and Stephan Wagner, Extremal Problems for Topological Indices in Combinatorial Chemistry, Journal of Computational Biology. September 2005, 12(7): 1004-1013.
Eric Weisstein's World of Mathematics, Complete Graph
Eric Weisstein's World of Mathematics, Hosoya Index
Eric Weisstein's World of Mathematics, Independent Edge Set
Eric Weisstein's World of Mathematics, Independent Vertex Set
Eric Weisstein's World of Mathematics, Matching
Eric Weisstein's World of Mathematics, Permutation Involution
Eric Weisstein's World of Mathematics, Triangular Graph
Eric Weisstein's World of Mathematics, Vertex Cover
Eric Weisstein's World of Mathematics, Young Tableau
Allan Wentink, Symmetry with Independent Blocks, on the late Ralph E. Griswold webpage of sample-chapters of an unfinished book about weaving, http://www.cs.arizona.edu/patterns/weaving/webdocs.html. [Cached copy]
Wikipedia, Telephone numbers.
Wikipedia, Young tableau.
Index entries for "core" sequences
Index entries for sequences related to Young tableaux.
Index entries for related partition-counting sequences

Cf. A001147, A001470, A053529, A099174, A136281-A136283, A001680, A001681.
See also A005425 for another version of the switchboard problem.
Equals 2 * A001475(n-1) for n>1.
First column of array A099020.
A069943(n+1)/A069944(n+1) = a(n)/A000142(n) in lowest terms.
Cf. A152736, A128229. - Gary W. Adamson, Dec 12 2008
Diagonal of A182172. - Alois P. Heinz, May 30 2012
Cf. A001813, A006882, A135401, A297708. - Manfred Scheucher, Jan 07 2018
Row sums of: A047884, A049403, A096713 (absolute value), A100861, A104556 (absolute value), A111924, A117506 (M_4 numbers), A122848, A238123.
A320663/A339888 count unlabeled multiset partitions into singletons/pairs.
A322661 counts labeled covering half-loop-graphs.
A339742 counts factorizations into distinct primes or squarefree semiprimes.
Cf. A000110, A000258, A000700, A000701, A321729, A321730, A321737.

a000085 n = a000085_list !! n
  a000085_list = 1 : 1 : zipWith (+)
    (zipWith (*) [1..] a000085_list) (tail a000085_list) -- Reinhard Zumkeller, May 16 2013

A000085 := proc(n) option remember; if n=0 then 1 elif n=1 then 1 else procname(n-1)+(n-1)*procname(n-2); fi; end;
with(combstruct):ZL3:=[S,{S=Set(Cycle(Z,card<3))}, labeled]:seq(count(ZL3,size=n),n=0..25); # Zerinvary Lajos, Sep 24 2007
with (combstruct):a:=proc(m) [ZL, {ZL=Set(Cycle(Z, m>=card))}, labeled]; end: A:=a(2):seq(count(A, size=n), n=0..25); # Zerinvary Lajos, Jun 11 2008

<Roger L. Bagula, Oct 06 2006 *)
With[{nn=30},CoefficientList[Series[Exp[x+x^2/2],{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, May 28 2013 *)
a[ n_] := Sum[(2 k - 1)!! Binomial[ n, 2 k], {k, 0, n/2}]; (* Michael Somos, Jun 01 2013 *)
a[ n_] := If[ n < 0, 0, HypergeometricU[ -n/2, 1/2, -1/2] / (-1/2)^(n/2)]; (* Michael Somos, Jun 01 2013 *)
a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ Exp[ x + x^2 / 2], {x, 0, n}]]; (* Michael Somos, Jun 01 2013 *)
Table[(I/Sqrt[2])^n HermiteH[n, -I/Sqrt[2]], {n, 0, 100}] (* Emanuele Munarini, Mar 02 2016 *)
a[n_] := Sum[StirlingS1[n, k]*2^k*BellB[k, 1/2], {k, 0, n}]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Jul 18 2017, after Emanuele Munarini *)
RecurrenceTable[{a[n] == a[n-1] + (n-1)*a[n-2], a[0] == 1, a[1] == 1}, a, {n, 0, 20}] (* Joan Ludevid, Jun 17 2022 *)
sds[{}]:={{}};sds[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sds[Complement[set,s]]]/@Cases[Subsets[set,{1,2}],{i,_}]; Table[Length[sds[Range[n]]],{n,0,10}] (* Gus Wiseman, Jan 11 2021 *)

B(n,x):=sum(stirling2(n,k)*x^k,k,0,n);
  a(n):=sum(stirling1(n,k)*2^k*B(k,1/2),k,0,n);
  makelist(a(n),n,0,40); /* Emanuele Munarini, May 16 2014 */

makelist((%i/sqrt(2))^n*hermite(n,-%i/sqrt(2)),n,0,12); /* Emanuele Munarini, Mar 02 2016 */

{a(n) = if( n<0, 0, n! * polcoeff( exp( x + x^2 / 2 + x * O(x^n)), n))}; /* Michael Somos, Nov 15 2002 */

N=66; x='x+O('x^N); egf=exp(x+x^2/2); Vec(serlaplace(egf)) \\ Joerg Arndt, Mar 07 2013

from math import factorial
def A000085(n): return sum(factorial(n)//(factorial(n-(k<<1))*factorial(k)*(1<>1)+1)) # Chai Wah Wu, Aug 31 2023

A000085 = lambda n: hypergeometric([-n/2,(1-n)/2], [], 2)
[simplify(A000085(n)) for n in range(28)] # Peter Luschny, Aug 21 2014

def a85(n): return sum(factorial(n) / (factorial(n-2*k) * 2**k * factorial(k)) for k in range(1+n//2))
for n in range(100): print(n, a85(n)) # Manfred Scheucher, Jan 07 2018

D-finite with recurrence a(0) = a(1) = 1, a(n) = a(n-1) + (n-1)*a(n-2) for n>1.
E.g.f.: exp(x+x^2/2).
a(n) = a(n-1) + A013989(n-2) = A013989(n)/(n+1) = 1+A001189(n).
a(n) = Sum_{k=0..floor(n/2)} n!/((n-2*k)!*2^k*k!).
a(m+n) = Sum_{k>=0} k!*binomial(m, k)*binomial(n, k)*a(m-k)*a(n-k). - Philippe Deléham, Mar 05 2004
For n>1, a(n) = 2*(A000900(n) + A000902(floor(n/2))). - Max Alekseyev, Oct 31 2015
The e.g.f. y(x) satisfies y^2 = y''y' - (y')^2.
a(n) ~ c*(n/e)^(n/2)exp(n^(1/2)) where c=2^(-1/2)exp(-1/4). [Chowla]
a(n) = HermiteH(n, 1/(sqrt(2)*i))/(-sqrt(2)*i)^n, where HermiteH are the Hermite polynomials. - Karol A. Penson, May 16 2002
a(n) = Sum_{k=0..n} A001498((n+k)/2, (n-k)/2)(1+(-1)^(n-k))/2. - Paul Barry, Jan 12 2006
For asymptotics see the Robinson paper.
a(n) = Sum_{m=0..n} A099174(n,m). - Roger L. Bagula, Oct 06 2006
O.g.f.: A(x) = 1/(1-x-1*x^2/(1-x-2*x^2/(1-x-3*x^2/(1-... -x-n*x^2/(1- ...))))) (continued fraction). - Paul D. Hanna, Jan 17 2006
From Gary W. Adamson, Dec 29 2008: (Start)
a(n) = (n-1)*a(n-2) + a(n-1); e.g., a(7) = 232 = 6*26 + 76.
Starting with offset 1 = eigensequence of triangle A128229. (End)
a(n) = (1/sqrt(2*Pi))*Integral_{x=-oo..oo} exp(-x^2/2)*(x+1)^n. - Groux Roland, Mar 14 2011
Row sums of |A096713|. a(n) = D^n(exp(x)) evaluated at x = 0, where D is the operator sqrt(1+2*x)*d/dx. Cf. A047974 and A080599. - Peter Bala, Dec 07 2011
From Sergei N. Gladkovskii, Dec 03 2011 - Oct 28 2013: (Start)
Continued fractions:
E.g.f.: 1+x*(2+x)/(2*G(0)-x*(2+x)) where G(k)=1+x*(x+2)/(2+2*(k+1)/G(k+1)).
G.f.: 1/(U(0) - x) where U(k) = 1 + x*(k+1) - x*(k+1)/(1 - x/U(k+1)).
G.f.: 1/Q(0) where Q(k) = 1 + x*k - x/(1 - x*(k+1)/Q(k+1)).
G.f.: -1/(x*Q(0)) where Q(k) = 1 - 1/x - (k+1)/Q(k+1).
G.f.: T(0)/(1-x) where T(k) = 1 - x^2*(k+1)/( x^2*(k+1) - (1-x)^2/T(k+1)). (End)
a(n) ~ (1/sqrt(2)) * exp(sqrt(n)-n/2-1/4) * n^(n/2) * (1 + 7/(24*sqrt(n))). - Vaclav Kotesovec, Mar 07 2014
a(n) = Sum_{k=0..n} s(n,k)*(-1)^(n-k)*2^k*B(k,1/2), where the s(n,k) are (signless) Stirling numbers of the first kind, and the B(n,x) = Sum_{k=0..n} S(n,k)*x^k are the Stirling polynomials, where the S(n,k) are the Stirling numbers of the second kind. - Emanuele Munarini, May 16 2014
a(n) = hyper2F0([-n/2,(1-n)/2],[],2). - Peter Luschny, Aug 21 2014
0 = a(n)*(+a(n+1) + a(n+2) - a(n+3)) + a(n+1)*(-a(n+1) + a(n+2)) for all n in Z. - Michael Somos, Aug 22 2014
From Peter Bala, Oct 06 2021: (Start)
a(n+k) == a(n) (mod k) for all n >= 0 and all positive odd integers k.
Hence for each odd k, the sequence obtained by taking a(n) modulo k is a periodic sequence and the exact period divides k. For example, taking a(n) modulo 7 gives the purely periodic sequence [1, 1, 2, 4, 3, 5, 6, 1, 1, 2, 4, 3, 5, 6, 1, 1, 2, 4, 3, 5, 6, ...] of period 7. For similar results see A047974 and A115329. (End)

A000165 Double factorial of even numbers: (2n)!! = 2^n*n!.

Original entry on oeis.org

1, 2, 8, 48, 384, 3840, 46080, 645120, 10321920, 185794560, 3715891200, 81749606400, 1961990553600, 51011754393600, 1428329123020800, 42849873690624000, 1371195958099968000, 46620662575398912000, 1678343852714360832000, 63777066403145711616000

Offset: 0

N. J. A. Sloane

a(n) is also the size of the automorphism group of the graph (edge graph) of the n-dimensional hypercube and also of the geometric automorphism group of the hypercube (the two groups are isomorphic). This group is an extension of an elementary Abelian group (C_2)^n by S_n. (C_2 is the cyclic group with two elements and S_n is the symmetric group.) - Avi Peretz (njk(AT)netvision.net.il), Feb 21 2001
Then a(n) appears in the power series: sqrt(1+sin(y)) = Sum_{n>=0} (-1)^floor(n/2)*y^(n)/a(n) and sqrt((1+cos(y))/2) = Sum_{n>=0} (-1)^n*y^(2n)/a(2n). - Benoit Cloitre, Feb 02 2002
Appears to be the BinomialMean transform of A001907. See A075271. - John W. Layman, Sep 28 2002
Number of n X n monomial matrices with entries 0, +-1.
Also number of linear signed orders.
Define a "downgrade" to be the permutation d which places the items of a permutation p in descending order. This note concerns those permutations that are equal to their double-downgrades. The number of permutations of order 2n having this property are equinumerous with those of order 2n+1. a(n) = number of double-downgrading permutations of order 2n and 2n+1. - Eugene McDonnell (eemcd(AT)mac.com), Oct 27 2003
a(n) = (Integral_{x=0..Pi/2} cos(x)^(2*n+1) dx) where the denominators are b(n) = (2*n)!/(n!*2^n). - Al Hakanson (hawkuu(AT)excite.com), Mar 02 2004
1 + (1/2)x - (1/8)x^2 - (1/48)x^3 + (1/384)x^4 + ... = sqrt(1+sin(x)).
a(n)*(-1)^n = coefficient of the leading term of the (n+1)-th derivative of arctan(x), see Hildebrand link. - Reinhard Zumkeller, Jan 14 2006
a(n) is the Pfaffian of the skew-symmetric 2n X 2n matrix whose (i,j) entry is j for iDavid Callan, Sep 25 2006
a(n) is the number of increasing plane trees with n+1 edges. (In a plane tree, each subtree of the root is an ordered tree but the subtrees of the root may be cyclically rotated.) Increasing means the vertices are labeled 0,1,2,...,n+1 and each child has a greater label than its parent. Cf. A001147 for increasing ordered trees, A000142 for increasing unordered trees and A000111 for increasing 0-1-2 trees. - David Callan, Dec 22 2006
Hamed Hatami and Pooya Hatami prove that this is an upper bound on the cardinality of any minimal dominating set in C_{2n+1}^n, the Cartesian product of n copies of the cycle of size 2n+1, where 2n+1 is a prime. - Jonathan Vos Post, Jan 03 2007
This sequence and (1,-2,0,0,0,0,...) form a reciprocal pair under the list partition transform and associated operations described in A133314. - Tom Copeland, Oct 29 2007
a(n) = number of permutations of the multiset {1,1,2,2,...,n,n,n+1,n+1} such that between the two occurrences of i, there is exactly one entry >i, for i=1,2,...,n. Example: a(2) = 8 counts 121323, 131232, 213123, 231213, 232131, 312132, 321312, 323121. Proof: There is always exactly one entry between the two 1s (when n>=1). Given a permutation p in A(n) (counted by a(n)), record the position i of the first 1, then delete both 1s and subtract 1 from every entry to get a permutation q in A(n-1). The mapping p -> (i,q) is a bijection from A(n) to the Cartesian product [1,2n] X A(n-1). - David Callan, Nov 29 2007
Row sums of A028338. - Paul Barry, Feb 07 2009
a(n) is the number of ways to seat n married couples in a row so that everyone is next to their spouse. Compare A007060. - Geoffrey Critzer, Mar 29 2009
From Gary W. Adamson, Apr 21 2009: (Start)
Equals (-1)^n * (1, 1, 2, 8, 48, ...) dot (1, -3, 5, -7, 9, ...).
Example: a(4) = 384 = (1, 1, 2, 8, 48) dot (1, -3, 5, -7, 9) = (1, -3, 10, -56, 432). (End)
exp(x/2) = Sum_{n>=0} x^n/a(n). - Jaume Oliver Lafont, Sep 07 2009
Assuming n starts at 0, a(n) appears to be the number of Gray codes on n bits. It certainly is the number of Gray codes on n bits isomorphic to the canonical one. Proof: There are 2^n different starting positions for each code. Also, each code has a particular pattern of bit positions that are flipped (for instance, 1 2 1 3 1 2 1 for n=3), and these bit position patterns can be permuted in n! ways. - D. J. Schreffler (ds1404(AT)txstate.edu), Jul 18 2010
E.g.f. of 0,1,2,8,... is x/(1-2x/(2-2x/(3-8x/(4-8x/(5-18x/(6-18x/(7-... (continued fraction). - Paul Barry, Jan 17 2011
Number of increasing 2-colored trees with choice of two colors for each edge. In general, if we replace 2 with k we get the number of increasing k-colored trees. For example, for k=3 we get the triple factorial numbers. - Wenjin Woan, May 31 2011
a(n) = row sums of triangle A193229. - Gary W. Adamson, Jul 18 2011
Also the number of permutations of 2n (or of 2n+1) that are equal to their reverse-complements. (See the Egge reference.) Note that the double-downgrade described in the preceding comment (McDonnell) is equivalent to the reverse-complement. - Justin M. Troyka, Aug 11 2011
The e.g.f. can be used to form a generator, [1/(1-2x)] d/dx, for A000108, so a(n) can be applied to A145271 to generate the Catalan numbers. - Tom Copeland, Oct 01 2011
The e.g.f. of 1/a(n) is BesselI(0,sqrt(2*x)). See Abramowitz-Stegun (reference and link under A008277), p. 375, 9.6.10. - Wolfdieter Lang, Jan 09 2012
a(n) = order of the largest imprimitive group of degree 2n with n systems of imprimitivity (see [Miller], p. 203). - L. Edson Jeffery, Feb 05 2012
Row sums of triangle A208057. - Gary W. Adamson, Feb 22 2012
a(n) is the number of ways to designate a subset of elements in each n-permutation. a(n) = A000142(n) + A001563(n) + A001804(n) + A001805(n) + A001806(n) + A001807(n) + A035038(n) * n!. - Geoffrey Critzer, Nov 08 2012
For n>1, a(n) is the order of the Coxeter groups (also called Weyl groups) of types B_n and C_n. - Tom Edgar, Nov 05 2013
For m>0, k*a(m-1) is the m-th cumulant of the chi-squared probability distribution for k degrees of freedom. - Stanislav Sykora, Jun 27 2014
a(n) with 0 prepended is the binomial transform of A120765. - Vladimir Reshetnikov, Oct 28 2015
Exponential self-convolution of A001147. - Vladimir Reshetnikov, Oct 08 2016
Also the order of the automorphism group of the n-ladder rung graph. - Eric W. Weisstein, Jul 22 2017
a(n) is the order of the group O_n(Z) = {A in M_n(Z): A*A^T = I_n}, the group of n X n orthogonal matrices over the integers. - Jianing Song, Mar 29 2021
a(n) is the number of ways to tile a (3n,3n)-benzel or a (3n+1,3n+2)-benzel using left stones and two kinds of bones; see Defant et al., below. - James Propp, Jul 22 2023
a(n) is the number of labeled histories for a labeled topology with the modified lodgepole shape and n+1 cherry nodes. - Noah A Rosenberg, Jan 16 2025

			The following permutations and their reversals are all of the permutations of order 5 having the double-downgrade property:
  0 1 2 3 4
  0 3 2 1 4
  1 0 2 4 3
  1 4 2 0 3
G.f. = 1 + 2*x + 8*x^2 + 48*x^3 + 384*x^4 + 3840*x^5 + 46080*x^6 + 645120*x^7 + ...

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. D. Noe, Table of n, a(n) for n = 0..100
Paul Barry, General Eulerian Polynomials as Moments Using Exponential Riordan Arrays, Journal of Integer Sequences, 16 (2013), #13.9.6.
Paul Barry, Generalized Eulerian Triangles and Some Special Production Matrices, arXiv:1803.10297 [math.CO], 2018.
Isabel Cação, Helmuth R. Malonek, Maria Irene Falcão, and Graça Tomaz, Combinatorial Identities Associated with a Multidimensional Polynomial Sequence, J. Int. Seq., Vol. 21 (2018), Article 18.7.4.
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
CombOS - Combinatorial Object Server, Generate colored permutations
R. Coquereaux and J.-B. Zuber, Maps, immersions and permutations, arXiv preprint arXiv:1507.03163 [math.CO], 2015.
Colin Defant, Rupert Li, James Propp, and Benjamin Young, Tilings of Benzels via the Abacus Bijection, arXiv preprint, arXiv:2209.05717 [math.CO], 2022.
Eric S. Egge, Restricted symmetric permutations, Ann. Combin., 11 (2007), 405-434.
Peter C. Fishburn, Signed Orders, Choice Probabilities and Linear Polytopes, Journal of Mathematical Psychology, Volume 45, Issue 1, (2001), pp. 53-80.
Joël Gay and Vincent Pilaud, The weak order on Weyl posets, arXiv:1804.06572 [math.CO], 2018.
G. Gordon, The answer is 2^n*n! What is the question?, Amer. Math. Monthly, 106 (1999), 636-645.
Guo-Niu Han, Enumeration of Standard Puzzles
Guo-Niu Han, Enumeration of Standard Puzzles [Cached copy]
Hamed Hatami and Pooya Hatami, Perfect dominating sets in the Cartesian products of prime cycles, arXiv:math/0701018 [math.CO], 2006-2009.
Jason D. Hildebrand, Differentiating Arctan(x)
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 136
E. Lappo and N. A. Rosenberg, A lattice structure for ancestral configurations arising from the relationship between gene trees and species trees, Adv. Appl. Math. 343 (2024), 65-81.
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181. [Annotated scan of pages 180 and 181 only]
E. Lucas, Théorie des nombres (annotated scans of a few selected pages)
Eugene McDonnell, Magic Squares and Permutations, APL Quote Quad 7.3 (Fall 1976).
B. E. Meserve, Double Factorials, American Mathematical Monthly, 55 (1948), 425-426.
G. A. Miller, Groups formed by special matrices, Bull. Amer. Math. Soc. 24 (1918), 203-206.
R. Ondrejka, Tables of double factorials, Math. Comp., 24 (1970), 231.
Alexsandar Petojevic, The Function vM_m(s; a; z) and Some Well-Known Sequences, Journal of Integer Sequences, Vol. 5 (2002), Article 02.1.7.
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014.
R. W. Robinson, Counting arrangements of bishops, pp. 198-214 of Combinatorial Mathematics IV (Adelaide 1975), Lect. Notes Math., 560 (1976).
M. Z. Spivey and L. L. Steil, The k-Binomial Transforms and the Hankel Transform, J. Integ. Seqs. Vol. 9 (2006), #06.1.1.
Eric Weisstein's World of Mathematics, Double Factorial
Eric Weisstein's World of Mathematics, Graph Automorphism
Eric Weisstein's World of Mathematics, Ladder Rung Graph
Index to divisibility sequences
Index entries for sequences related to factorial numbers

Cf. A000142 (n!), A001147 ((2n-1)!!), A032184 (2^n*(n-1)!).
Cf. A000079, A000108, A000111, A001044, A001563, A001804, A001805, A001806.
Cf. A001807, A001813, A001907, A002454, A006882, A007060, A008277, A008544.
Cf. A010050, A019774, A028326, A028338, A035038, A039683, A047053, A047055.
Cf. A047058, A047657, A075271, A084947, A084948, A084949, A092605, A133314.
Cf. A145271, A193229, A208057.
This sequence gives the row sums in A060187, and (-1)^n*a(n) the alternating row sums in A039757.
Also row sums in A028338.
Column k=2 of A329070.

a000165 n = product [2, 4 .. 2 * n]  -- Reinhard Zumkeller, Mar 28 2015

[2^n*Factorial(n): n in [0..35]]; // Vincenzo Librandi, Apr 22 2011

I:=[2,8]; [1] cat [n le 2 select I[n]  else (3*n-1)*Self(n-1)-2*(n-1)^2*Self(n-2): n in [1..35] ]; // Vincenzo Librandi, Feb 19 2015

A000165 := proc(n) option remember; if n <= 1 then 1 else n*A000165(n-2); fi; end;
ZL:=[S, {a = Atom, b = Atom, S = Prod(X,Sequence(Prod(X,b))), X = Sequence(b,card >= 0)}, labelled]: seq(combstruct[count](ZL, size=n), n=0..17); # Zerinvary Lajos, Mar 26 2008
G(x):=(1-2*x)^(-1): f[0]:=G(x): for n from 1 to 29 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n],n=0..17); # Zerinvary Lajos, Apr 03 2009
A000165 := proc(n) doublefactorial(2*n) ; end proc; seq(A000165(n),n=0..10) ; # R. J. Mathar, Oct 20 2009

Table[(2 n)!!, {n, 30}] (* Vladimir Joseph Stephan Orlovsky, Dec 13 2008 *)
(2 Range[0, 30])!! (* Harvey P. Dale, Jan 23 2015 *)
RecurrenceTable[{a[n] == 2 n*a[n-1], a[0] == 1}, a, {n,0,30}] (* Ray Chandler, Jul 30 2015 *)

a(n)=n!<Charles R Greathouse IV, Feb 11 2011

{a(n) = prod( k=1, n, 2*k)}; /* Michael Somos, Jan 04 2013 */

from math import factorial
def A000165(n): return factorial(n)<Chai Wah Wu, Jan 24 2023

[2^n*factorial(n) for n in range(31)] # G. C. Greubel, Jul 21 2024

E.g.f.: 1/(1-2*x).
a(n) = A001044(n)/A000142(n)*A000079(n) = Product_{i=0..n-1} (2*i+2) = 2^n*Pochhammer(1,n). - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
D-finite with recurrence a(n) = 2*n * a(n-1), n>0, a(0)=1. - Paul Barry, Aug 26 2004
This is the binomial mean transform of A001907. See Spivey and Steil (2006). - Michael Z. Spivey (mspivey(AT)ups.edu), Feb 26 2006
a(n) = Integral_{x>=0} x^n*exp(-x/2)/2 dx. - Paul Barry, Jan 28 2008
G.f.: 1/(1-2x/(1-2x/(1-4x/(1-4x/(1-6x/(1-6x/(1-.... (continued fraction). - Paul Barry, Feb 07 2009
a(n) = A006882(2*n). - R. J. Mathar, Oct 20 2009
From Gary W. Adamson, Jul 18 2011: (Start)
a(n) = upper left term in M^n, M = a production matrix (twice Pascal's triangle deleting the first "2", with the rest zeros; cf. A028326):
  2,  2,  0,  0,  0,  0, ...
  2,  4,  2,  0,  0,  0, ...
  2,  6,  6,  2,  0,  0, ...
  2,  8, 12,  8,  2,  0, ...
  2, 10, 20, 20, 10,  2, ...
  ... (End)
From Sergei N. Gladkovskii, Apr 11 2013, May 01 2013, May 24 2013, Sep 30 2013, Oct 27 2013: (Start)
Continued fractions:
G.f.: 1 + x*(Q(0) - 1)/(x+1) where Q(k) = 1 + (2*k+2)/(1-x/(x+1/Q(k+1))).
G.f.: 1/Q(0) where Q(k) = 1 + 2*k*x - 2*x*(k+1)/Q(k+1).
G.f.: G(0)/2 where G(k) = 1 + 1/(1 - x*(2*k+2)/(x*(2*k+2) + 1/G(k+1))).
G.f.: 1/Q(0) where Q(k) = 1 - x*(4*k+2) - 4*x^2*(k+1)^2/Q(k+1).
G.f.: R(0) where R(k) = 1 - x*(2*k+2)/(x*(2*k+2)-1/(1-x*(2*k+2)/(x*(2*k+2) -1/R(k+1)))). (End)
a(n) = (2n-2)*a(n-2) + (2n-1)*a(n-1), n>1. - Ivan N. Ianakiev, Aug 06 2013
From Peter Bala, Feb 18 2015: (Start)
Recurrence equation: a(n) = (3*n - 1)*a(n-1) - 2*(n - 1)^2*a(n-2) with a(1) = 2 and a(2) = 8.
The sequence b(n) = A068102(n) also satisfies this second-order recurrence. This leads to the generalized continued fraction expansion lim_{n -> oo} b(n)/a(n) = log(2) = 1/(2 - 2/(5 - 8/(8 - 18/(11 - ... - 2*(n - 1)^2/((3*n - 1) - ... ))))). (End)
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = sqrt(e) (A019774).
Sum_{n>=0} (-1)^n/a(n) = 1/sqrt(e) (A092605). (End)
Limit_{n->oo} a(n)^4 / (n * A134372(n)) = Pi. - Daniel Suteu, Apr 09 2022
a(n) = 1/([x^n] hypergeom([1], [1], x/2)). - Peter Luschny, Sep 13 2024
a(n) = Sum_{k=0..n} k!*(n-k)!*binomial(n,k)^2. - Ridouane Oudra, Jul 13 2025

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

N. J. A. Sloane

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

			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

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

N. J. A. Sloane, Table of n, a(n) for n = 0..100
Marcello Artioli, Giuseppe Dattoli, Silvia Licciardi, and Rosa Maria Pidatella, Hermite and Laguerre Functions: a Unifying Point of View, Università degli Studi di Catania, Sicily, Italy (2019).
Murray Bremner and Martin Markl, Distributive laws between the Three Graces, arXiv:1809.08191 [math.AT], 2018.
R. B. Brent, Generalizing Tuenter's Binomial Sums, J. Int. Seq. 18 (2015) # 15.3.2.
Peter J. Cameron, Some treelike objects, Quart. J. Math. Oxford Ser. (2) 38 (1987), no. 150, 155--183. MR0891613 (89a:05009). See p. 155. - _N. J. A. Sloane_, Apr 18 2014
Peter J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
Elliot J. Carr and Matthew J. Simpson, New homogenization approaches for stochastic transport through heterogeneous media, arXiv:1810.08890 [physics.bio-ph], 2018.
W. Y. C. Chen, L. W. Shapiro and L. L. M. Young, Parity reversing involutions on plane trees and 2-Motzkin paths, arXiv:math/0503300 [math.CO], 2005.
Ali Chouria, Vlad-Florin Drǎgoi, and Jean-Gabriel Luque, On recursively defined combinatorial classes and labelled trees, arXiv:2004.04203 [math.CO], 2020.
P. Cvitanovic, Group theory for Feynman diagrams in non-Abelian gauge theories, Phys. Rev. D14 (1976), 1536-1553.
Nick Early, Honeycomb tessellations and canonical bases for permutohedral blades, arXiv:1810.03246 [math.CO], 2018.
John Engbers, David Galvin, and Clifford Smyth, Restricted Stirling and Lah numbers and their inverses, arXiv:1610.05803 [math.CO], 2016. See p. 8.
P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009; see page 127
S. Goodenough and C. Lavault, On subsets of Riordan subgroups and Heisenberg--Weyl algebra, arXiv preprint arXiv:1404.1894 [cs.DM], 2014.
S. Goodenough and C. Lavault, Overview on Heisenberg—Weyl Algebra and Subsets of Riordan Subgroups, The Electronic Journal of Combinatorics, 22(4) (2015), #P4.16,
H. W. Gould, A class of binomial sums and a series transform, Utilitas Math., 45 (1994), 71-83. (Annotated scanned copy)
A. M. Ibrahim, Extension of factorial concept to negative numbers, Notes on Number Theory and Discrete Mathematics, Vol. 19, 2013, 2, 30-42.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 115
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181. [Annotated scan of pages 180 and 181 only]
Jesús Leaños, Rutilo Moreno, and Luis Manuel Rivera-Martínez, On the number of mth roots of permutations, Australas. J. Combin. 52 (2012), 41-54 (Theorem 1).
Jesús Leaños, Rutilo Moreno, and Luis Manuel Rivera-Martínez, On the number of mth roots of permutations, arXiv:1005.1531 [math.CO], 2010-2011.
Édouard Lucas, Théorie des Nombres, Gauthier-Villars, Paris, 1891, Vol. 1, p. 221.
R. J. Marsh and P. P. Martin, Pascal arrays: counting Catalan sets, arXiv:math/0612572 [math.CO], 2006.
Calin D. Morosan, On the number of broadcast schemes in networks, Information Processing Letters, Volume 100, Issue 5 (2006), 188-193.
R. A. Proctor, Let's Expand Rota's Twelvefold Way for Counting Partitions!, arXiv:math/0606404 [math.CO], 2006-2007.
C. Radoux, Déterminants de Hankel et théorème de Sylvester, Séminaire Lotharingien de Combinatoire, B28b (1992), 9 pp.
J. Riordan, Letter to N. J. A. Sloane, Feb 03 1975 (with notes by njas)
H. E. Salzer, Coefficients for expressing the first thirty powers in terms of the Hermite polynomials, Math. Comp., 3 (1948), 167-169.
Elena L. Wang and Guoce Xin, On Ward Numbers and Increasing Schröder Trees, arXiv:2507.15654 [math.CO], 2025. See pp. 12-13.
Index to divisibility sequences
Index entries for related partition-counting sequences

Cf. A037224, A048854, A001147, A007696, A008545, A122670 (essentially the same sequence), A000165, A047055, A047657, A084947, A084948, A084949, A010050, A000142, A008275, A000108, A000984, A008276, A000680, A094216.
Catalan(n-1)*M^(n-1)*n! for M=1,2,3,4,5,6: A001813, A052714 (or A144828), A221954, A052734, A221953, A221955.
Cf. A123023, A059344, A099174, A066325, A001700, A000407, A006882.
Cf. A000085, A135401, A297708. - Manfred Scheucher, Jan 07 2018

List([0..20],n->Factorial(2*n)/Factorial(n)); # Muniru A Asiru, Nov 01 2018

[Factorial(2*n)/Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 09 2018

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

Table[(2n)!/n!, {n,0,20}] (* Harvey P. Dale, May 02 2011 *)

makelist(binomial(n+n, n)*n!,n,0,30); /* Martin Ettl, Nov 05 2012 */

a(n)=binomial(n+n,n)*n! \\ Charles R Greathouse IV, Jun 15 2011

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)

from math import factorial
def A001813(n): return factorial(n<<1)//factorial(n) # Chai Wah Wu, Feb 14 2023

[binomial(2*n,n)*factorial(n) for n in range(0, 17)] # Zerinvary Lajos, Dec 03 2009

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

More terms from James Sellers, May 01 2000

A000407 a(n) = (2*n+1)! / n!.

Original entry on oeis.org

1, 6, 60, 840, 15120, 332640, 8648640, 259459200, 8821612800, 335221286400, 14079294028800, 647647525324800, 32382376266240000, 1748648318376960000, 101421602465863680000, 6288139352883548160000, 415017197290314178560000

Offset: 0

N. J. A. Sloane

The e.g.f. of 1/a(n) = n!/(2*n+1)! is (exp(sqrt(x)) - exp(-sqrt(x)))/(2*sqrt(x)). - Wolfdieter Lang, Jan 09 2012
Product of the larger parts of the partitions of 2n+2 into exactly two parts. - Wesley Ivan Hurt, Jun 15 2013
For n > 0, a(n-1) = (2n-1)!/(n-1)!, the number of ways n people can line up in n labeled queues. The derivation is straightforward. Person 1 has (2n-1) choices - be first in line in one of the queues or get behind one of the other people. Person 2 has (2n-2) choices - choose one of the n queues or get behind one of the remaining n-2 people. Continuing in this fashion, we finally find that person n has to choose one of the n queues. - Dennis P. Walsh, Mar 24 2016
For n > 0, a(n-1) is the number of functions f:[n]->[2n] that are acyclic and injective. Note that f is acyclic if, for all x in [n], x is not a member of the set {f(x),f(f(x)), f(f(f(x))), ...}. - Dennis P. Walsh, Mar 25 2016
a(n) is the number of labeled maximal outerplanar graphs with n-3 vertices. - Allan Bickle, Feb 19 2024

			G.f. = 1 + 6*x + 60*x^2 + 840*x^3 + 15120*x^4 + 332640*x^5 + 8648640*x^6 + ...
For n=1 the a(1)=6 ways for 2 people to line up in 2 queues are as follows: Q1<P1,P2> Q2<>, Q1<P2,P1> Q2<>, Q1<P1> Q2<P2>, Q1<P2> Q2<P1>, Q1<> Q2<P1,P2>, Q1<> Q2<P2,P1>. - _Dennis P. Walsh_, Mar 24 2016
For the unique maximal outerplanar graph with 4 vertices, there are C(4,2)=6 ways to label the two degree 3 vertices, and the other two labels are forced.  Thus a(1) = 6.

L. W. Beineke and R. E. Pippert, Enumerating labeled k-dimensional trees and ball dissections, pp. 12-26 of Proceedings of Second Chapel Hill Conference on Combinatorial Mathematics and Its Applications, University of North Carolina, Chapel Hill, 1970. Reprinted with a slightly different title in Math. Annalen, 191 (1971), 87-98.
L. B. W. Jolley, Summation of Series, Dover, 1961.
Loren C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181.
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).

N. J. A. Sloane, Table of n, a(n) for n = 0..100
L. W. Beineke and R. E. Pippert, Enumerating labeled k-dimensional trees and ball dissections, Proceedings of Second Chapel Hill Conference on Combinatorial Mathematics and Its Applications, University of North Carolina, Chapel Hill, 1970, pp. 12-26. Reprinted with a slightly different title in Math. Annalen, Vol. 191 (1971), pp. 87-98.
Allan Bickle, A Survey of Maximal k-degenerate Graphs and k-Trees, Theory and Applications of Graphs 0 1 (2024) Article 5.
Peter J. Cameron, Sequences Realized by Oligomorphic Permutation Groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
G. Hu, Catalan number and enumeration of maximal outerplanar graphs, Tsinghua Science and Technology 25 5 1 (2000), 109-114.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 139.
Loren C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181. [Annotated scan of pages 180 and 181 only]
Dan Levy and Lior Pachter, The Neighbor-Net Algorithm, arXiv:math/0702515 [math.CO], 2007-2008.
Lee A. Newberg, The Number of Clone Orderings, Discrete Applied Mathematics, Vol. 69, No. 3 (1996), pp. 233-245.
J.-C. Novelli and J.-Y. Thibon, Hopf Algebras of m-permutations, (m+1)-ary trees, and m-parking functions, arXiv preprint arXiv:1403.5962 [math.CO], 2014.
Robert W. Robinson, Counting arrangements of bishops, pp. 198-214 of Combinatorial Mathematics IV (Adelaide 1975), Lect. Notes Math., 560 (1976).
Herbert E. Salzer, Coefficients for expressing the first thirty powers in terms of the Hermite polynomials, Math. Comp., Vol. 3, No. 23 (1948), pp. 167-169.
Herbert E. Salzer, Orthogonal polynomials arising in the evaluation of inverse Laplace transforms, Math. Comp. Vol. 9, No. 52 (1955), pp. 164-177.
Herbert E. Salzer, Orthogonal polynomials arising in the evaluation of inverse Laplace transforms, Math. Comp., Vol. 9, No. 52 (1955), 164-177. [Annotated scanned copy]
Maxie D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications , J. Int. Seq., Vol. 13 (2010), Article 10.6.7, page 39.
Wikipedia, William Brouncker, 2nd Viscount Brouncker.
Index to divisibility sequences
Index entries for sequences related to factorial numbers

Cf. A001761-A001763, A007696, A024199, A073010.
A100622 is the "Number of topologically distinct solutions to the clone ordering problem for n clones" without the restriction that they be in a single contig (see [Newberg] for definition of contig).
Column m=0 of A292219.

[Factorial(2*n+1) / Factorial(n): n in [0..20]]; // Vincenzo Librandi, Jun 16 2015

For Maple program see A000903.
a := n -> pochhammer(n+1,n+1); (for n>=0) # Peter Luschny, Feb 14 2009

Table[(2n + 1)!/n!, {n, 0, 30}] (* Stefan Steinerberger, Apr 08 2006 *)
a[ n_] := If[ n < 0, 1/2, 1] Pochhammer[ n + 1, n + 1]; (* Michael Somos, Jan 03 2015 *)
a[ n_] := Which[ n < -1, -(-1)^n / (4 a[-n - 2]), n == -1, 1/2, True, (2 n + 1)! / n!]; (* Michael Somos, Jan 03 2015 *)

A000407(n):=(2*n+1)!/n!$
makelist(A000407(n),n,0,30); /* Martin Ettl, Nov 05 2012 */

a(n)=(2*n+1)!/n! \\ Charles R Greathouse IV, Jan 12 2012

{a(n) = if( n<-1, -(-1)^n / (4 * a(-n-2)), n==-1, 1/2, (2*n + 1)! / n!)}; /* Michael Somos, Jan 03 2015 */

E.g.f.: (1 - 4*x)^(-3/2). - Michael Somos, Jan 03 2015
E.g.f.: Sum_{k>=0} a(k+2) * x^k / k! = (1 - 2*x - sqrt(1 - 4*x)) / 4.
E.g.f. for a(n-1), n >= 0, with a(-1) := 0 is (-1+1/(1-4*x)^(1/2))/2. 2*a(n) = (4*n+2)(!^4) := Product_{j=0..n} (4*j + 2), (one half of 4-factorial numbers). - Wolfdieter Lang
a(n) = C(n+1)*(n+2)!/2 for all n>=0. - Paul Barry, Feb 16 2005
For n>1, a(n) = (1/2)*A001813(n+1). - Zerinvary Lajos, Jun 06 2007
For asymptotics see the Robinson paper.
Sum_{n >=0} n!/a(n) = 2*Pi/3^(3/2) =  1.2091995761... = A248897 [Jolley eq 261]
G.f.: 1 / (1 - 6*x / (1 - 4*x / (1 - 10*x / (1 - 8*x / (1 - 14*x / ... ))))). - Michael Somos, May 12 2012
G.f.: 1/Q(0), where Q(k) = 1 + 2*(2*k-1)*x - 4*x*(k+1)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 03 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - 2*x/(2*x + 1/(2*k+3)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 02 2013
a(n) = -(-1)^n / (4 * a(-2-n)) = a(n-1) * (4*n+2) for all n in Z. - Michael Somos, Jan 03 2015
a(n) = A087299(2*n + 1). - Michael Somos, Jan 03 2015
From Peter Bala, Feb 16 2015: (Start)
Recurrence equation: a(n) = 4*a(n-1) + 4*(2*n - 1)^2*a(n-2) with a(0) = 1 and a(1) = 6.
The integer sequence b(n) := a(n)*Sum_{k = 0..n} (-1)^k/(2*k + 1), beginning [1, 4, 52, 608, 12624, ...], satisfies the same second-order recurrence equation. This leads to Brouncker's generalized continued fraction expansion Sum_{k >= 0} (-1)^k/(2*k + 1) = Pi/4 = 1/(1 + 1^2/(2 + 3^2/(2 + 5^2/(2 + ... )))). Note b(n) = 2^n*A024199(n+1).
Recurrence equation: a(n) = (5*n + 2)*a(n-1) - 2*n*(2*n - 1)^2*a(n-2) with a(0) = 1 and a(1) = 6.
The integer sequence c(n) := a(n)*Sum_{k = 0..n} k!^2/(2*k + 1)!, beginning [1, 7, 72, 1014, 18276, ... ], satisfies the same second-order recurrence equation. This leads to the generalized continued fraction expansion Sum_{k >= 0} k!^2/(2*k + 1)! = 2*Pi/sqrt(27) = 2*A073010 = 1/(1 - 1/(7 - 12/(12 - 30/(17 - ... - 2*n*(2*n - 1)/((5*n + 2) - ... ))))). (End)
a(n) = Product_{k=n+1..(2*n+1)} k. - Carlos Eduardo Olivieri, Jun 03 2015
From Ilya Gutkovskiy, Jan 17 2017: (Start)
a(n) ~ 2^(2*n+3/2)*n^(n+1)/exp(n).
Sum_{n>=0} 1/a(n) = exp(1/4)*sqrt(Pi)*erf(1/2) = 1.184593072938653151..., where erf() is the error function. (End)
Sum_{n>=0} (-1)^n/a(n) = exp(-1/4)*sqrt(Pi)*erfi(1/2), where erfi() is the imaginary error function. - Amiram Eldar, Jan 18 2021
It follows from the comments above that we have a(n) = a(n-1)*(4*n+2), with a(1) = 6, a(0) = 1.
a(n) = A081125(2*n+1). - R. J. Mathar, Jun 07 2025

A000898 a(n) = 2(a(n-1) + (n-1)a(n-2)) for n >= 2 with a(0) = 1.

Original entry on oeis.org

1, 2, 6, 20, 76, 312, 1384, 6512, 32400, 168992, 921184, 5222208, 30710464, 186753920, 1171979904, 7573069568, 50305536256, 342949298688, 2396286830080, 17138748412928, 125336396368896, 936222729254912, 7136574106003456, 55466948299223040, 439216305474605056, 3540846129311916032

Offset: 0

N. J. A. Sloane

Number of solutions to the rook problem on a 2n X 2n board having a certain symmetry group (see Robinson for details).
Also the value of the n-th derivative of exp(x^2) evaluated at 1. - N. Calkin, Apr 22 2010
For n >= 1, a(n) is also the sum of the degrees of the irreducible representations of the group of n X n signed permutation matrices (described in sequence A066051). The similar sum for the "ordinary" symmetric group S_n is in sequence A000085. - Sharon Sela (sharonsela(AT)hotmail.com), Jan 12 2002
It appears that this is also the number of permutations of 1, 2, ..., n+1 such that each term (after the first) is within 2 of some preceding term. Verified for n+1 <= 6. E.g., a(4) = 20 because of the 24 permutations of 1, 2, 3, 4, the only ones not permitted are 1, 4, 2, 3; 1, 4, 3, 2; 4, 1, 2, 3; and 4, 1, 3, 2. - Gerry Myerson, Aug 06 2003
Hankel transform is A108400. - Paul Barry, Feb 11 2008
From Emeric Deutsch, Jun 19 2010: (Start)
Number of symmetric involutions of [2n]. Example: a(2)=6 because we have 1234, 2143, 1324, 3412, 4231, and 4321. See the Egge reference, pp. 419-420.
Number of symmetric involutions of [2n+1]. Example: a(2)=6 because we have 12345, 14325, 21354, 45312, 52341, and 54321. See the Egge reference, pp. 419-420.
(End)
Binomial convolution of sequence A000085: a(n) = Sum_{k=0..n} binomial(n,k)*A000085(k)*A000085(n-k). - Emanuele Munarini, Mar 02 2016
The sequence can be obtained from the infinite product of 2 X 2 matrices [(1,N); (1,1)] by extracting the upper left terms, where N = (1, 3, 5, ...), the odd integers. - Gary W. Adamson, Jul 28 2016
Apparently a(n) is the number of standard domino tableaux of size 2n, where a domino tableau is a generalized Young tableau in which all rows and columns are weakly increasing and all regions are dominos. - Gus Wiseman, Feb 25 2018

			G.f. = 1 + 2*x + 6*x^2 + 20*x^3 + 76*x^4 + 312*x^5 + 1384*x^6 + 6512*x^7 + ...
The a(3) = 20 domino tableaux:
1 1 2 2 3 3
.
1 2 2 3 3
1
.
1 2 3 3   1 1 3 3   1 1 2 2
1 2       2 2       3 3
.
1 1 3 3   1 1 2 2
2         3
2         3
.
1 2 3   1 2 2   1 1 3
1 2 3   1 3 3   2 2 3
.
1 3 3   1 2 2
1       1
2       3
2       3
.
1 2   1 1   1 1
1 2   2 3   2 2
3 3   2 3   3 3
.
1 3   1 2   1 1
1 3   1 2   2 2
2     3     3
2     3     3
.
1 1
2
2
3
3
.
1
1
2
2
3
3 - _Gus Wiseman_, Feb 25 2018

D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 3, Sect 5.1.4 Exer. 31.
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181.
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).

T. D. Noe, Table of n, a(n) for n=0..200
Y. Alp and E. G. Kocer, Exponential Almost-Riordan Arrays, Results Math 79, 173 (2024). See page 10.
Arvind Ayyer, Hiranya Kishore Dey, and Digjoy Paul, How large is the character degree sum compared to the character table sum for a finite group?, arXiv:2406.06036 [math.RT], 2024. See p. 10.
C. Banderier, M. Bousquet-Mélou, A. Denise, P. Flajolet, D. Gardy and D. Gouyou-Beauchamps, Generating functions for generating trees, Discrete Mathematics 246(1-3), March 2002, pp. 29-55.
Paul Barry, The Gamma-Vectors of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1804.05027 [math.CO], 2018.
R. A. Brualdi, Shi-Mie Ma, Enumeration of involutions by descents and symmetric matrices, Eur. J. Combin. 43 (2015) 220-228
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
C.-O. Chow, Counting involutory, unimodal, and alternating signed permutations, Discr. Math. 306 (2006), 2222-2228.
R. Donaghey, Binomial self-inverse sequences and tangent coefficients, J. Combin. Theory, Series A, 21 (1976), 155-163.
Eric S. Egge, Restricted symmetric permutations, Ann. Combin., 11 (2007), 405-434.
Adam M. Goyt and Lara K. Pudwell, Avoiding colored partitions of two elements in the pattern sense, arXiv preprint arXiv:1203.3786 [math.CO], 2012. - From _N. J. A. Sloane_, Sep 17 2012
T. Halverson and M. Reeks, Gelfand Models for Diagram Algebras, arXiv:1302.6150 [math.RT], 2013-2014.
Guo-Niu Han, Enumeration of Standard Puzzles, 2011. [Cached copy]
Guo-Niu Han, Enumeration of Standard Puzzles, arXiv:2006.14070 [math.CO], 2020.
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181. [Annotated scan of pages 180 and 181 only]
Yen-chi R. Lin, Asymptotic Formula for Symmetric Involutions, arXiv:1310.0988 [math.CO], 2013.
E. Lucas, Théorie des Nombres, Gauthier-Villars, Paris, 1891, Vol. 1, p. 221.
E. Lucas, Théorie des nombres (annotated scans of a few selected pages)
J. Riordan, Letter to N. J. A. Sloane, Feb 03 1975 (with notes by njas)
D. P. Roberts and A. Venkatesh, Hurwitz monodromy and full number fields, 2014.
Index entries for sequences related to Hermite polynomials

Column k=2 of A294042 and A376826.
Row sums of A122832 and A185296.
Cf. A000085, A000712, A004003, A066051, A099390, A100861, A135401, A138178, A153452, A297388, A299699, A299926, A300056, A300060.

a000898 n = a000898_list !! n
a000898_list = 1 : 2 : (map (* 2) $
   zipWith (+) (tail a000898_list) (zipWith (*) [1..] a000898_list))
-- Reinhard Zumkeller, Oct 10 2011

# For Maple program see A000903.
seq(simplify((-I)^n*HermiteH(n, I)), n=0..25); # Peter Luschny, Oct 23 2015

a[n_] := Sum[ 2^k*StirlingS1[n, k]*BellB[k], {k, 0, n}]; Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Nov 17 2011, after Vladeta Jovovic *)
RecurrenceTable[{a[0]==1,a[1]==2,a[n]==2(a[n-1]+(n-1)a[n-2])},a,{n,30}] (* Harvey P. Dale, Aug 04 2012 *)
Table[Abs[HermiteH[n, I]], {n, 0, 20}] (* Vladimir Reshetnikov, Oct 22 2015 *)
a[ n_] := Sum[ 2^(n - 2 k) n! / (k! (n - 2 k)!), {k, 0, n/2}]; (* Michael Somos, Oct 23 2015 *)

makelist((%i)^n*hermite(n,-%i),n,0,12); /* Emanuele Munarini, Mar 02 2016 */

{a(n) = if( n<0, 0, n! * polcoeff( exp(2*x + x^2 + x * O(x^n)), n))}; /* Michael Somos, Feb 08 2004 */

{a(n) = if( n<2, max(0, n+1), 2*a(n-1) + (2*n - 2) * a(n-2))}; /* Michael Somos, Feb 08 2004 */

my(x='x+O('x^66)); Vec(serlaplace(exp(2*x+x^2))) \\ Joerg Arndt, Oct 04 2013

{a(n) = sum(k=0, n\2, 2^(n - 2*k) * n! / (k! * (n - 2*k)!))}; /* Michael Somos, Oct 23 2015 */

a(n) = Sum_{m=0..n} |A060821(n,m)| = H(n,-i)*i^n, with the Hermite polynomials H(n,x); i.e., these are row sums of the unsigned triangle A060821.
E.g.f.: exp(x*(x + 2)).
a(n) = 2 * A000902(n) for n >= 1.
a(n) = Sum_{k=0..n} binomial(n,2k)*binomial(2k,k)*k!*2^(n-2k). - N. Calkin, Apr 22 2010
Binomial transform of A047974. - Paul Barry, May 09 2003
a(n) = Sum_{k=0..n} Stirling1(n, k)*2^k*Bell(k). - Vladeta Jovovic, Oct 01 2003
From Paul Barry, Aug 29 2005: (Start)
a(n) = Sum_{k=0..floor(n/2)} A001498(n-k, k) * 2^(n-k).
a(n) = Sum_{k=0..n} A001498((n+k)/2, (n-k)/2) * 2^((n+k)/2) * (1+(-1)^(n-k))/2. (End)
For asymptotics, see the Robinson paper. [This is disputed by Yen-chi R. Lin. See below, Sep 30 2013.]
a(n) = Sum_{k=0..floor(n/2)} 2^(n-2*k) * C(n,2*k) * (2*k)!/k!. - Paul Barry, Feb 11 2008
G.f.: 1/(1 - 2*x - 2*x^2/(1 - 2*x - 4*x^2/(1 - 2*x - 6*x^2/(1 - 2*x - 8*x^2/(1 - ... (continued fraction). - Paul Barry, Feb 25 2010
E.g.f.: exp(x^2 + 2*x) = Q(0); Q(k) = 1 + (x^2 + 2*x)/(2*k + 1 - (x^2 + 2*x)*(2*k + 1)/((x^2 + 2*x) + (2*k + 2)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Nov 24 2011
G.f.: 1/Q(0), where Q(k) = 1 + 2*x*k - x - x/(1 - 2*x*(k + 1)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 07 2013
a(n) = (2*n/e)^(n/2) * exp(sqrt(2*n)) / sqrt(2*e) * (1 + sqrt(2/n)/3 + O(n^(-1))). - Yen-chi R. Lin, Sep 30 2013
0 = a(n)*(2*a(n+1) + 2*a(n+2) - a(n+3)) + a(n+1)*(-2*a(n+1) + a(n+2)) for all n >= 0. - Michael Somos, Oct 23 2015
a(n) = Sum_{k=0..floor(n/2)} 2^(n-k)*B(n, k), where B are the Bessel numbers A100861. - Peter Luschny, Jun 04 2021

More terms from Larry Reeves (larryr(AT)acm.org), Feb 21 2001
Initial condition a(0)=1 added to definition by Jon E. Schoenfield, Oct 01 2013
More terms from Joerg Arndt, Oct 04 2013

cp's OEIS Frontend

A000085 Number of self-inverse permutations on n letters, also known as involutions; number of standard Young tableaux with n cells.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

Maxima

PARI

PARI

Python

Sage

Sage

Formula

A000165 Double factorial of even numbers: (2n)!! = 2^n*n!.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

PARI

Python

SageMath

Formula

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

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

Maxima

PARI

PARI

Python

Sage

Formula

Extensions

A000407 a(n) = (2*n+1)! / n!.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

PARI

PARI

Formula

A000898 a(n) = 2(a(n-1) + (n-1)a(n-2)) for n >= 2 with a(0) = 1.

A000902 Expansion of e.g.f. (1/2)(exp(2x + x^2) + 1).