|
|
A067994
|
|
Hermite numbers.
|
|
10
|
|
|
1, 0, -2, 0, 12, 0, -120, 0, 1680, 0, -30240, 0, 665280, 0, -17297280, 0, 518918400, 0, -17643225600, 0, 670442572800, 0, -28158588057600, 0, 1295295050649600, 0, -64764752532480000, 0, 3497296636753920000, 0, -202843204931727360000, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
|a(n)| is the number of sets of ordered pairs of n labeled elements. - Steven Finch, Nov 14 2021
|a(n)| is the number of square roots of any permutation in S_{2n} whose disjoint cycle decomposition consists of n transpositions, n > 0. For n=2, permutation (1,2)(3,4) in S_4 has exactly |a(2)|=2 square roots: (1,3,2,4) and (1,4,2,3). - Luis Manuel Rivera Martínez, Feb 25 2015
Named after the French mathematician Charles Hermite (1822-1901). - Amiram Eldar, Jun 06 2021
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (-1)^(n/2)*n!/(n/2)! if n is even, 0 otherwise. - Mitch Harris, Feb 01 2006
E.g.f.: U(0) where U(k) = 1 - x^2/((2*k+1) - x^2*(2*k+1)/(x^2 - 2*(k+1)/U(k+1))); (continued fraction). - Sergei N. Gladkovskii, Oct 23 2012
G.f.: 1/G(0) where G(k) = 1 + 2*x^2*(k+1)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Dec 05 2012
E.g.f.: E(0)/(1+x) where E(k) = 1 + x/(1 - x/(x - (k+1)/E(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Apr 05 2013
E.g.f.: E(0)-1, where E(k) = 2 - x^2/(2*k+1 + x^2/E(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Dec 24 2013
a(n) = Pochhammer(-n, n/2). (End)
|
|
EXAMPLE
|
|a(4)| = 12 because the sets of ordered pairs for n = 4 are
{(1,2),(3,4)}, {(2,1),(3,4)}, {(1,2),(4,3)}, {(2,1),(4,3)},
{(1,3),(2,4)}, {(3,1),(2,4)}, {(1,3),(4,2)}, {(3,1),(4,2)},
{(1,4),(3,2)}, {(4,1),(3,2)}, {(1,4),(2,3)}, {(4,1),(2,3)}. (End)
|
|
MAPLE
|
A067994 := n -> pochhammer(-n, n/2):
|
|
MATHEMATICA
|
HermiteH[Range[0, 50], 0]
With[{nmax=50}, CoefficientList[Series[Exp[-x^2], {x, 0, nmax}], x]*Range[0, nmax]!] (* G. C. Greubel, Jun 09 2018 *)
|
|
PROG
|
(PARI) x='x+O('x^30); Vec(serlaplace(exp(-x^2))) \\ G. C. Greubel, Jun 09 2018
(Magma) m:=25; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(Exp(-x^2))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Jun 09 2018
|
|
CROSSREFS
|
Cf. A097388 (same sequence without zeros).
Cf. A101109 (ordered triples instead of ordered pairs).
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|