Double factorial

In mathematics, the product of all the integers from 1 up to some non-negative integer n that have the same parity as n is called the double factorial or semifactorial of n and is denoted by n!!.[1] That is,

\( n!! = \prod_{k=0}^m (n-2k) = n (n-2) (n-4) \cdots \)

where \( m=\lceil n/2 \rceil - 1 \). A consequence of this definition is that (as an empty product)

0!! = 1.

For example, 9!! = 1 × 3 × 5 × 7 × 9 = 945.

For even n the double factorial is

\(n!! = \prod_{k=1}^{n/2} (2k) = n(n-2)\cdots 2. \)

For odd n it is

\( n!! = \prod_{k=1}^{(n+1)/2} (2k-1) = n(n-2)\cdots 1. \)

The sequence of double factorials for even n = 0, 2, 4, 6, 8, ... starts as

1, 2, 8, 48, 384, 3840, 46080, 645120, .... (sequence A000165 in OEIS)

The sequence of double factorials for odd n = 1, 3, 5, 7, ... starts as

1, 3, 15, 105, 945, 10395, 135135, .... (sequence A001147 in OEIS)

Merserve (1948)[2] (possibly the earliest publication to use double factorial notation)[3] states that the double factorial was originally introduced in order to simplify the expression of certain trigonometric integrals arising in the derivation of the Wallis product. Double factorials also arise in expressing the volume of a hypersphere, and they have many applications in enumerative combinatorics.[1][4]

The term odd factorial is sometimes used for the double factorial of an odd number.[5]

Relation to the factorial

Because the double factorial only involves about half the factors of the ordinary factorial, its value is not substantially larger than the square root of the factorial n!, and it is much smaller than the iterated factorial (n!)!.

For an even positive integer n = 2k, k ≥ 0, the double factorial may be expressed as

\( n!! = 2^k k!. \)

For odd n = 2k − 1, k ≥ 1, it has the expressions

\( n!! = \frac{(2k)!}{2^k k!} = \frac{n!}{(n-1)!!}. \)

In this expression, the first denominator equals (2k)!! and cancels the unwanted even factors from the numerator.

For an odd positive integer n = 2k − 1, k ≥ 1, the double factorial may be expressed in terms of k-permutations of 2k as[1][3]

\( (2k-1)!! = \frac {_{2k}P_k} {2^k} = \frac {{(2k)}^{\underline k}} {2^k}. \)

Extensions
Negative arguments

The ordinary factorial, when extended to the Gamma function, has a pole at each negative integer, preventing the factorial from being defined at these numbers. However, the double factorial of odd numbers may be extended to any negative odd integer argument by inverting its recurrence relation

\( n!! = n \times (n-2)!! \)

to give

\( n!! = \frac{(n+2)!!}{n+2}.

Using this inverted recurrence, −1!! = 1, −3!! = −1, and −5!! = 1/3; negative odd numbers with greater magnitude have fractional double factorials.[1] In particular, this gives, when n is an odd number,

\( (-n)!! \times n!! = (-1)^{(n-1)/2} \times n. \)

Complex arguments

Disregarding the above definition of n!! for even values of n, the double factorial for odd integers can be extended to most real and complex numbers z by noting that when z is a positive odd integer then [6] [7]

\( z!! = z(z-2)\cdots (3) = 2^{(z-1)/2}\left(\frac{z}{2}\right)\left(\frac{z-2}{2}\right)\cdots \left(\frac{3}{2}\right)
= 2^{(z-1)/2} \frac{\Gamma\left(\frac{z}{2}+1\right)}{\Gamma\left(\frac{1}{2}+1\right)} = \sqrt{\frac{2^{z+1}}{\pi}} \Gamma\left(\frac{z}{2}+1\right) = \left(\frac{z}{2}\right)!\sqrt{\frac{2^{z+1}}{\pi}} \,. \)

From this one can derive an alternative definition of z!! for non-negative even integer values of z:

\( (2k)!!= \sqrt{ \frac{2}{\pi} } \prod_{i=1}^k (2i) = 2^k k! \sqrt{ \frac{2}{\pi} } \,, \)

with the value for 0!! in this case being

\( 0!! = \sqrt{ \frac{2}{\pi} } \approx 0.7978845608... \,. \)

The expression found for z!! is defined for all complex numbers except the negative even integers. Using it as the definition, the volume of an n-dimensional hypersphere of radius R can be expressed as[8]

\( V_n=\frac{2 (2\pi)^{(n-1)/2}}{n!!} R^n. \)

Applications in enumerative combinatorics

Double factorials are motivated by the fact that they occur frequently in enumerative combinatorics and other settings. For instance, n!! for odd values of n counts

Perfect matchings of the complete graph Kn + 1 for odd n. In such a graph, any single vertex v has n possible choices of vertex that it can be matched to, and once this choice is made the remaining problem is one of selecting a perfect matching in a complete graph with two fewer vertices. For instance, a complete graph with four vertices a, b, c, and d has three perfect matchings: ab and cd, ac and bd, and ad and bc.[1] Perfect matchings may be described in several other equivalent ways, including involutions without fixed points on a set of n + 1 items (permutations in which each cycle is a pair)[1] or chord diagrams (sets of chords of a set of n + 1 points evenly spaced on a circle such that each point is the endpoint of exactly one chord, also called Brauer diagrams).[4][9][10] The numbers of matchings in complete graphs, without constraining the matchings to be perfect, are instead given by the telephone numbers, which may be expressed as a summation involving double factorials.[11]
Stirling permutations, permutations of the multiset of numbers 1, 1, 2, 2, ..., k, k in which each pair of equal numbers is separated only by larger numbers, where k = (n + 1)/2. The two copies of k must be adjacent; removing them from the permutation leaves a permutation in which the maximum element is k − 1, with n positions into which the adjacent pair of k values may be placed. From this recursive construction, a proof that the Stirling permutations are counted by the double permutations follows by induction.[1] Alternatively, instead of the restriction that values between a pair may be parger than it, one may also consider the permutations of this multiset in which the first copies of each pair appear in sorted order; such a permutation defines a matching on the 2k positions of the permutation, so again the number of permutations may be counted by the double permutations.[4]
Heap-ordered trees, trees with k + 1 nodes labeled 0, 1, 2, ... k, such that the root of the tree has label 0, each other node has a larger label than its parent, and such that the children of each node have a fixed ordering. An Euler tour of the tree (with doubled edges) gives a Stirling permutation, and every Stirling permutation represents a tree in this way.[1][12]
Unrooted binary trees with (n + 5)/2 labeled leaves. Each such tree may be formed from a tree with one fewer leaf, by subdividing one of the n tree edges and making the new vertex be the parent of a new leaf.
Rooted binary trees with (n + 3)/2 labeled leaves. This case is similar to the unrooted case, but the number of edges that can be subdivided is even, and in addition to subdividing an edge it is possible to add a node to a tree with one fewer leaf by adding a new root whose two children are the smaller tree and the new leaf.[1][4]

Callan (2009) and Dale & Moon (1993) list several additional objects with the same counting sequence, including "trapezoidal words" (numerals in a mixed radix system with increasing odd radixes), height-labeled Dyck paths, height-labeled ordered trees, "overhang paths", and certain vectors describing the lowest-numbered leaf descendant of each node in a rooted binary tree. For bijective proofs that some of these objects are equinumerous, see Rubey (2008) and Marsh & Martin (2011).[13][14]

The even double factorials give the numbers of elements of the hyperoctahedral groups (signed permutations or symmetries of a hypercube)
Additional identities

For integral values of n,

\( \int_{0}^{\pi/2}\sin^n x\,dx=\int_{0}^{\pi/2}\cos^n x\,dx=\frac{(n-1)!!}{n!!}\cdot {\begin{cases} 1 & n \text{ odd} \\ \frac{\pi}{2} & n \text{ even} \end{cases}}. \)

Using instead the extension of the double factorial of odd numbers to complex numbers, the formula is

\( \int_{0}^{\pi/2}\sin^n x\,dx=\int_{0}^{\pi/2}\cos^n x\,dx=\frac{(n-1)!!}{n!!}\cdot\sqrt{\frac{\pi}{2}}. \)

Double factorials can also be used to evaluate integrals of more complicated trigonometric polynomials.[2][15]

Double factorials of odd numbers are related to the gamma function by the identity:

\( (2n-1)!! = 2^n \cdot \frac{\Gamma(\frac{1}{2} + n)} {\sqrt{\pi}} = (-2)^n \cdot \frac{\sqrt{\pi}} { \Gamma(\frac{1}{2} - n)} \)

Some additional identities involving double factorials of odd numbers are:

\( (2n-1)!! = \sum_{k=1}^{n-1} \binom{n}{k+1} (2k-1)!! (2n-2k-3)!!.[1] \)

\( (2n-1)!! = \sum_{k=0}^{n} \binom{2n-k-1}{k-1} \frac{(2k-1)(2n-k+1)}{k+1}(2n-2k-3)!!.[1] \)

\( (2n-1)!! = \sum_{k=1}^{n} \frac{(n-1)!}{(k-1)!} k(2k-3)!!.[1] \)

References

Callan, David (2009), A combinatorial survey of identities for the double factorial, arXiv:0906.1317.
Meserve, B. E. (1948), "Classroom Notes: Double Factorials", The American Mathematical Monthly 55 (7): 425–426, doi:10.2307/2306136, MR 1527019
Gould, Henry; Quaintance, Jocelyn (2012), "Double fun with double factorials", Mathematics Magazine 85 (3): 177–192, doi:10.4169/math.mag.85.3.177, MR 2924154.
Dale, M. R. T.; Moon, J. W. (1993), "The permuted analogues of three Catalan sets", Journal of Statistical Planning and Inference 34 (1): 75–87, doi:10.1016/0378-3758(93)90035-5, MR 1209991.
E.g., in Henderson, Daniel J.; Parmeter, Christopher F. (2012), "Canonical higher-order kernels for density derivative estimation", Statistics & Probability Letters 82 (7): 1383–1387, doi:10.1016/j.spl.2012.03.013, MR 2929790 and Nielsen, B. (1999), "The likelihood-ratio test for rank in bivariate canonical correlation analysis", Biometrika 86 (2): 279–288, doi:10.1093/biomet/86.2.279, MR 1705359.
Hassani, Sadri (2000), Mathematical Methods: For Students of Physics and Related Fields, Undergraduate Texts in Mathematics, Springer, p. 266, ISBN 9780387989587.
Double factorial: Specific values (formula 06.02.03.0005), Wolfram Research, 2001-10-29, retrieved 2013-03-23.
Mezey, Paul G. (2009), "Some dimension problems in molecular databases", Journal of Mathematical Chemistry 45 (1): 1–6, doi:10.1007/s10910-008-9365-8.
Kitaev, Sergey (2011), Patterns in Permutations and Words, EATCS Monographs in Theoretical Computer Science, Springer, p. 96, ISBN 9783642173332.
Dale, M. R. T.; Narayana, T. V. (1986), "A partition of Catalan permuted sequences with applications", Journal of Statistical Planning and Inference 14 (2): 245–249, doi:10.1016/0378-3758(86)90161-8, MR 852528.
Tichy, Robert F.; Wagner, Stephan (2005), "Extremal problems for topological indices in combinatorial chemistry" (PDF), Journal of Computational Biology 12 (7): 1004–1013, doi:10.1089/cmb.2005.12.1004.
Janson, Svante (2008), "Plane recursive trees, Stirling permutations and an urn model", Fifth Colloquium on Mathematics and Computer Science, Discrete Math. Theor. Comput. Sci. Proc., AI, Assoc. Discrete Math. Theor. Comput. Sci., Nancy, pp. 541–547, arXiv:0803.1129, MR 2508813.
Rubey, Martin (2008), "Nestings of matchings and permutations and north steps in PDSAWs", 20th Annual International Conference on Formal Power Series and Algebraic Combinatorics (FPSAC 2008), Discrete Math. Theor. Comput. Sci. Proc., AJ, Assoc. Discrete Math. Theor. Comput. Sci., Nancy, pp. 691–704, MR 2721495.
Marsh, Robert J.; Martin, Paul (2011), "Tiling bijections between paths and Brauer diagrams", Journal of Algebraic Combinatorics 33 (3): 427–453, arXiv:0906.0912, doi:10.1007/s10801-010-0252-6, MR 2772541.

Dassios, George; Kiriaki, Kiriakie (1987), "A useful application of Gauss theorem", Bulletin de la Société Mathématique de Grèce 28 (part A): 40–43, MR 935868.

External links

Weisstein, Eric W., "Double Factorial", MathWorld.

It is unknown whether there exists a prime number p such that Cp is also prime.

As of February 2012, the largest known generalized Cullen prime is 427194 × 113 427194 + 1. It has 877,069 digits and was discovered by a PrimeGrid participant from United States.[3]

150 = 2 × 3 × 52

Mathematics Encyclopedia