A001147 -id:A001147 - cp's OEIS frontend

A053871 a(0)=1; a(1)=0; a(n) = 2(n-1)(a(n-1) + a(n-2)).

1, 0, 2, 8, 60, 544, 6040, 79008, 1190672, 20314880, 387099936, 8148296320, 187778717632, 4702248334848, 127140703364480, 3691602647581184, 114562300528369920, 3784124901630435328, 132555364873399378432, 4908221631901073295360, 191549525877429961604096

Offset: 0

Cris Moore (moore(AT)santafe.edu) and Christian G. Bower, Mar 29 2000

Number of deranged matchings of 2n people with partners (of either sex) other than their spouse. 2n objects are initially paired in some way and then are re-paired so that no object is with its original partner (the dancing problem in the article).
Of interest in the "collision problem", where, given a 2-to-1 function f, we are asked for x, y such that f(x)=f(y).
2^n*n!*a(n) = (2n)! b(n) where b(n) are the probabilities that appear in Margolis (2001). One interpretation is in terms of matchings or 1-factors of the complete graph on 2n vertices. The number of these is (2n)!/2^{n}n!. The number of 1-factors being disjoint from (that is, having no edges in common with) a given 1-factor is then a(n); and then b(n) is the probability of picking such a disjoint one factor at random.
Also n!*a(n) = 2^n * c(n) where c(n) = A001499(n). If we define d(n,k) = k(n-1)(d(n-1,k) + d(n-2,k)), with d(0,k) = 1 and d(1,k) = 0, so d(n,1) are the derangement numbers A000166, then a(n) = d(n,2) (cf. A033030, A088991). On the other hand, taking d*(n,k) = d(n,k)/k^{n}, we have d*(n,k) = (n-1)(d*(n-1,k) + d*(n-2,k)/k), with d*(0,k) = 1 and d*(1,k) = 0 and it is easy to see from Bricard's recurrence for c(n) that c(n)/n! satisfies this for k = 2.
A proof that the description in the first comment as "number of deranged matchings" implies the defining recursion relation: let (x_1, y_1), (x_2, y_2), ..., (x_n, y_n) be the given pairs. In a deranged matching x_1 will be paired with any of the 2(n-1) objects x_2, y_2, x_3, y_3, ..., x_n, y_n. It is sufficient to count only those deranged matchings in which x_1, is matched with x_2. They are of two kinds: (i) y_1 is not matched with y_2; their number is a(n-1); (ii) y_1 is matched with y_2; their number is a(n-2). - Emeric Deutsch, Jan 23 2009
Inverse binomial transform of the odd double factorials (A001147). - David Callan, Aug 25 2009
From Lewis Mammel (l_mammel(AT)att.net), Oct 07 2009: (Start)
The formula is given directly by the Principle of Inclusion and Exclusion.
The first term includes all pairings, the second term excludes all pairings containing each pair, the third term includes all pairings containing each pair of pairs, and so on.
Based on n-a -> n for large n, the ratio a(n)/(2n-1)!! -> exp(-1/2) ~= 0.60653.
We find a(n)/(2n-1)!! for n= 100, 200, 300, 400 ~= 0.6050124904, 0.6057720290, 0.6060250088, 0.6061514604. (End)
This is a subset of the set of pairings of the first 2n integers (A001147) in another way: forbidding pairs of the form (2k,2k+1) for all k as well as the pair (1,2n) (seeing the constraint as circular by opposition to the linear A165968). - Olivier Gérard, Feb 08 2011
a(n) is the n-th central moment of a central chi-squared distribution (with 1 degree of freedom), i.e., a(n) = E[ (Y- E[Y] )^n ] = E[ (X^2 - 1 )^n ] where Y is chi-squared, X is std normal, X~N(0,1), and the expectation operator is E[]. - David Fioramonti, May 11 2016
Exponential self-convolution of a(n)/2^n gives subfactorials (A000166). - Vladimir Reshetnikov, Oct 09 2016
For n > 1, also the number of maximum matchings in the n-cocktail party graph. - Eric W. Weisstein, Jun 15 2017
Number of polygonal tiles with n sides with two exits per side and n edges connecting pairs of exits, with no edges between exits on the same side (cf. A289191, A289269 for nonisomorphic tiles under rotational and rotational and reflectional symmetries). - Marko Riedel, Jun 28 2017
Number of ways n people can hold hands (every hand is holding another hand) where no one is holding their own hand. - Harry Richman, Aug 29 2023

R. Bricard, L'Intermédiaire des Mathématiciens, 8 (1901), 312-313.

Michael De Vlieger, Table of n, a(n) for n = 0..404 (First 100 terms from _T. D. Noe_)
A. Ayyer, Determinants and Perfect Matchings, arXiv:1106.1465 [math.CO], 2011.
Roland Bacher, Counting Packings of Generic Subsets in Finite Groups, Electr. J. Combinatorics, 19 (2012), #P7. - From _N. J. A. Sloane_, Feb 06 2013
Bela Bollobas, The number of 1-factors in 2k-connected graphs, J. Combin. Theory Ser. B 25 (1978), no. 3, 363--366. MR0516268 (80m:05060) - _N. J. A. Sloane_, Mar 26 2012
James N. Brawner, Dinner, Dancing and Tennis, Anyone?, Mathematics Magazine, Vol. 73, No 1 (2000).
M. A. Brodie, Avoiding your spouse at a party leads to war, Math. Mag., 75 (2002), 203-208.
Davi B. Costa, Bogdan A. Dobrescu, and Patrick J. Fox, Chiral Abelian gauge theories with few fermions, arXiv:2001.11991 [hep-ph], 2020.
Barbara H. Margolius, Dinner-Diner Matching Probabilities
B. H. Margolius, Avoiding your spouse at a bridge party, Math. Mag., 74 (2001), 33-41.
B. H. Margolius, The Dinner-Diner Matching Problem, Mathematics Magazine, 76 (2003), 107-118.
Eric Weisstein's World of Mathematics, Cocktail Party Graph
Eric Weisstein's World of Mathematics, Independent Edge Set
Eric Weisstein's World of Mathematics, Matching
Eric Weisstein's World of Mathematics, Maximum Independent Edge Set

Cf. A054479, A001147, A001818, A000166, A001499.
See A289191 for when rotational symmetries of the tiles are taken into account. - Marko Riedel, Jun 28 2017
Cf. A165968, number of pairings of 2n things disjoint to a given pairing, and containing a given pair not in the given pairing. It is given by a(n)/(2n-2). - Lewis Mammel (l_mammel(AT)att.net), Oct 07 2009

a053871 n = a053871_list !! n
a053871_list = 1 : 0 : zipWith (*)
   [2,4..] (zipWith (+) a053871_list $ tail a053871_list)
-- Reinhard Zumkeller, Mar 07 2012

f:= gfun:-rectoproc({a(0) = 1, a(1) = 0, a(n) = 2*(n - 1)*(a(n - 1) + a(n - 2))},a(n),remember):
map(f, [$0..30]); # Robert Israel, May 10 2016

RecurrenceTable[{a[0]==1,a[1]==0,a[n]==2(n-1)(a[n-1]+a[n-2])}, a[n],{n,20}] (* Harvey P. Dale, Sep 15 2011 *)
CoefficientList[Assuming[{Element[x, Reals], x>0}, Series[Sqrt[Pi/2] * (I + Erfi[Sqrt[(1+1/x)/2]]) / (E^((1+x)/(2*x)) * Sqrt[x*(x+1)]), {x, 0, 20}]], x] (* Vaclav Kotesovec, Feb 15 2015 *)
Range[0, 20]! CoefficientList[Series[1/(Exp[x] Sqrt[1 - 2 x]), {x, 0, 20}], x] (* Eric W. Weisstein, Jun 15 2017 *)
Table[(-1)^n HypergeometricPFQ[{1/2, -n}, {}, 2], {n, 20}] (* Eric W. Weisstein, Jun 15 2017 *)
Table[I (-1)^n HypergeometricU[1/2, 3/2 + n, -1/2]/Sqrt[2], {n, 20}] (* Eric W. Weisstein, Dec 31 2017 *)

a(n)=(-1)^(n+1)*sum(k=0,n,(-1)^k*binomial(n,k)*prod(i=0,k,2*i-1))

a(n) = A054479(n)/A001147(n).
E.g.f.: 1/(exp(x)*sqrt(1-2x)).
a(n) = (-1)^n*Sum_{k=0..n} (-1)^k*C(n, k)*(2*k-1)!!. - Benoit Cloitre, May 01 2003; corrected by David Fioramonti, May 17 2016
a(n) = Integral_{x>=0} (x-1)^n * (exp(-x/2)/sqrt(2*Pi*x)) dt. - Paul Barry, Apr 12 2010
Conjectured g.f.: T(0)/(1+x), where T(k) = 1 - x*(k+1)/(x*(k+1) - (1+x)/T(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Nov 13 2013
a(n) ~ 2^(n+1/2) * n^n / exp(n+1/2). - Vaclav Kotesovec, Mar 11 2014
G.f.: Sum_{k>=0} (2*k - 1)!!*x^k/(1 + x)^(k+1). - Ilya Gutkovskiy, Apr 12 2019
Conjecture: if m == n (mod q) for q odd, then (-1)^m*a(m) == (-1)^n*a(n) (mod q). - Harry Richman, Aug 29 2023

More terms from James Sellers, Apr 08 2000

A060540 Square array read by antidiagonals downwards: T(n,k) = (nk)!/(k!^nn!), (n>=1, k>=1), the number of ways of dividing nk labeled items into n unlabeled boxes with k items in each box.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 10, 15, 1, 1, 35, 280, 105, 1, 1, 126, 5775, 15400, 945, 1, 1, 462, 126126, 2627625, 1401400, 10395, 1, 1, 1716, 2858856, 488864376, 2546168625, 190590400, 135135, 1, 1, 6435, 66512160, 96197645544, 5194672859376, 4509264634875, 36212176000, 2027025, 1

Offset: 1

Henry Bottomley, Apr 02 2001

The Copeland link gives the associations of this entry with the operator calculus of Appell Sheffer polynomials, the combinatorics of simple set partitions encoded in the Faa di Bruno formula for composition of analytic functions (formal Taylor series), the Pascal matrix, and the geometry of the n-dimensional simplices (hypertriangles, or hypertetrahedra). These, in turn, are related to simple instances of the application of the exponential formula / principle / schema giving the number of not-necessarily-connected objects composed from an ensemble of connected objects. - Tom Copeland, Jun 09 2021

			Array begins:
  1,   1,       1,          1,             1,                 1, ...
  1,   3,      10,         35,           126,               462, ...
  1,  15,     280,       5775,        126126,           2858856, ...
  1, 105,   15400,    2627625,     488864376,       96197645544, ...
  1, 945, 1401400, 2546168625, 5194672859376, 11423951396577720, ...
  ...

Seiichi Manyama, Antidiagonals n = 1..50, flattened (first 20 antidiagonals from Harry J. Smith)
Tom Copeland, Calculus, Combinatorics, and Geometry Underlying OEIS A060540, and the Exponential Formula, 2021.
Nattawut Phetmak and Jittat Fakcharoenphol, Uniformly Generating Derangements with Fixed Number of Cycles in Polynomial Time, Thai J. Math. (2023) Vol. 21, No. 4, 899-915. See pp. 901, 914.
Elena L. Wang and Guoce Xin, On Ward Numbers and Increasing Schröder Trees, arXiv:2507.15654 [math.CO], 2025. See p. 13.

Rows include A000012, A001700, A060542, A082368, A322252.
Columns k=1..10 give A000012, A001147, A025035, A025036, A025037, A025038, A025040, A025041, A025042.
Main diagonal is A057599.
Related to A057599, see also A096126 and A246048.
Cf. A060358, A361948 (includes row/col 0).
Cf. A000217, A000292, A000332, A000389, A000579, A000580, A007318, A036040, A099174, A133314, A132440, A135278 (associations in Copeland link).

T[n_, k_] := (n*k)!/(k!^n*n!);
Table[T[n-k+1, k], {n, 1, 10}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Jun 29 2018 *)

{ i=0; for (m=1, 20, for (n=1, m, k=m - n + 1; write("b060540.txt", i++, " ", (n*k)!/(k!^n*n!))); ) } \\ Harry J. Smith, Jul 06 2009

T(n,k) = (n*k)!/(k!^n*n!) = T(n-1,k)*A060543(n,k) = A060538(n,k)/k!.

T(n,k) = Product_{j=2..n} binomial(j*k-1,k-1). - M. F. Hasler, Aug 22 2014

Definition reworded by M. F. Hasler, Aug 23 2014

A112007 Coefficient triangle for polynomials used for o.g.f.s for unsigned Stirling1 diagonals.

Original entry on oeis.org

1, 2, 1, 6, 8, 1, 24, 58, 22, 1, 120, 444, 328, 52, 1, 720, 3708, 4400, 1452, 114, 1, 5040, 33984, 58140, 32120, 5610, 240, 1, 40320, 341136, 785304, 644020, 195800, 19950, 494, 1, 362880, 3733920, 11026296, 12440064, 5765500, 1062500, 67260, 1004, 1

Offset: 0

Wolfdieter Lang, Sep 12 2005

This is the row reversed second-order Eulerian triangle A008517(k+1,k+1-m). For references see A008517.
The o.g.f. for the k-th diagonal, k >= 1, of the unsigned Stirling1 triangle |A008275| is G1(1,x)=1/(1-x) if k=1 and G1(k,x) = g1(k-2,x)/(1-x)^(2*k-1), if k >= 2, with the row polynomials g1(k;x):=Sum_{m=0..k} a(k,m)*x^m.
The recurrence eq. for the row polynomials is g1(k,x)=((k+1)+k*x)*g1(k-1,x) + x*(1-x)*(d/dx)g1(k-1,x), k >= 1, with input g1(0,x):=1.
The column sequences start with A000142 (factorials), A002538, A002539, A112008, A112485.
This o.g.f. computation was inspired by Bender et al. article where the Stirling polynomials have been rediscussed.
The A163936 triangle is identical to the triangle given above except for an extra right hand column [1, 0, 0, 0, ... ]. The A163936 triangle is related to the higher order exponential integrals E(x,m,n), see A163931 and A163932. - Johannes W. Meijer, Oct 16 2009

			Triangle begins:
    1;
    2,   1;
    6,   8,   1;
   24,  58,  22,   1;
  120, 444, 328,  52,   1;
  ...
G.f. for k=3 sequence A000914(n-1), [2,11,35,85,175,322,546,...], is G1(3,x)= g1(1,x)/(1-x)^5= (2+x)/(1-x)^5.

Robert Israel, Table of n, a(n) for n = 0..10010 (rows 0 to 140, flattened)
Peter Bala, Diagonals of triangles with generating function exp(t*F(x))
C. M. Bender, D. C. Brody and B. K. Meister, Bernoulli-like polynomials associated with Stirling Numbers, arXiv:math-ph/0509008 [math-ph], 2005.
E. Burlachenko, Composition polynomials of the RNA matrix and B-composition polynomials of the Riordan pseudo-involution, arXiv:1907.12272 [math.NT], 2019.
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms
Wolfdieter Lang, First 10 rows.
Wolfdieter Lang, On Generating functions of Diagonals Sequences of Sheffer and Riordan Number Triangles, arXiv:1708.01421 [math.NT], August 2017.
R. Paris, A uniform asymptotic expansion for the incomplete gamma function, Journal of Computational and Applied Mathematics, 148 (2002), p. 223-239. See 234. [_Tom Copeland_, Jan 03 2016]
Andrew Elvey Price, Alan D. Sokal, Phylogenetic trees, augmented perfect matchings, and a Thron-type continued fraction (T-fraction) for the Ward polynomials, arXiv:2001.01468 [math.CO], 2020.

Row sums give A001147(k+1) = (2*k+1)!!, k>=0.

a:= proc(k,m) option remember; if m >= 0 and k >= 0 then (k+m+1)*procname(k-1,m)+(k-m+1)*procname(k-1,m-1) else 0 fi end proc:
a(0,0):= 1:
seq(seq(a(k,m),m=0..k),k=0..10); # Robert Israel, Jul 20 2017

a[k_, m_] = Sum[(-1)^(k + n + 1)*Binomial[2k + 3, n]*StirlingS1[m + k - n + 2, m + 1 - n], {n, 0, m}]; Flatten[Table[a[k, m], {k, 0, 8}, {m, 0, k}]][[1 ;; 45]] (* Jean-François Alcover, Jun 01 2011, after Johannes W. Meijer *)

a(k, m)=sum(n=0, m, (-1)^(k + n + 1)*binomial(2*k + 3, n)*stirling(m + k - n + 2, m + 1 - n, 1));
for(k=0, 10, for(m=0, k, print1(a(k, m),", "))) \\ Indranil Ghosh, Jul 21 2017

a(k, m) = (k+m+1)*a(k-1, m) + (k-m+1)*a(k-1, m-1), if k >= m >= 0, a(0, 0)=1; a(k, -1):=0, otherwise 0.
a(k,m) = Sum_{n=0..m} (-1)^(k+n+1)*C(2*k+3,n)*Stirling1(m+k-n+2,m+1-n). - Johannes W. Meijer, Oct 16 2009
The compositional inverse (with respect to x) of y = y(t,x) = (x+t*log(1-x)) is  x = x(t,y) = 1/(1-t)*y + t/(1-t)^3*y^2/2! + (2*t+t^2)/(1-t)^5*y^3/3! + (6*t+8*t^2+t^3)/(1-t)^7*y^4/4! + .... The numerator polynomials of the rational functions in t are the row polynomials of this triangle. As observed above, the rational functions in t are the generating functions for the diagonals of |A008275|. See the Bala link for a proof. Cf. A008517. - Peter Bala, Dec 02 2011

A028338 Triangle of coefficients in expansion of (x+1)(x+3)...*(x + 2n - 1) in rising powers of x.

Original entry on oeis.org

1, 1, 1, 3, 4, 1, 15, 23, 9, 1, 105, 176, 86, 16, 1, 945, 1689, 950, 230, 25, 1, 10395, 19524, 12139, 3480, 505, 36, 1, 135135, 264207, 177331, 57379, 10045, 973, 49, 1, 2027025, 4098240, 2924172, 1038016, 208054, 24640, 1708, 64, 1, 34459425, 71697105, 53809164, 20570444, 4574934, 626934, 53676, 2796, 81, 1

Offset: 0

Bill Gosper

Exponential Riordan array (1/sqrt(1-2*x), log(1/sqrt(1-2*x))). - Paul Barry, May 09 2011

The o.g.f.s D(d, x) of the column sequences, for d, d >= 0,(d=0 for the main diagonal) are P(d, x)/(1 - x)^(2*d+1), with the row polynomial P(d, x) = Sum_{m=0..d} A288875(d, m)*x^m. See A288875 for details. - Wolfdieter Lang, Jul 21 2017

			G.f. for n = 4: (x + 1)*(x + 3)*(x + 5)*(x + 7) = 105 + 176*x + 86*x^2 + 16*x^3 + x^4.
The triangle T(n, k) begins:
n\k       0        1        2        3       4      5     6    7  8  9
0:        1
1:        1        1
2:        3        4        1
3:       15       23        9        1
4:      105      176       86       16       1
5:      945     1689      950      230      25      1
6:    10395    19524    12139     3480     505     36     1
7:   135135   264207   177331    57379   10045    973    49    1
8:  2027025  4098240  2924172  1038016  208054  24640  1708   64  1
9: 34459425 71697105 53809164 20570444 4574934 626934 53676 2796 81  1
...
row n = 10: 654729075 1396704420 1094071221 444647600 107494190 16486680 1646778 106800 4335 100 1.
...  reformatted and extended. - _Wolfdieter Lang_, May 09 2017
O.g.f.s of diagonals d >= 0: D(2, x) = (3 + 8*x + x^2)/(1 - x)^5 generating [3, 23, 86, ...] = A024196(n+1), from the row d=2 entries of A288875 [3, 8, 1]. - _Wolfdieter Lang_, Jul 21 2017
Boas-Buck recurrence for column k=2 and n=4: T(4, 2) = (4!/2)*(2*(1+4*(5/12))*T(2,2)/2! + 1*(1 + 4*(1/2))*T(3,2)/3!) = (4!/2)*(8/3*1 + 3*9/3!) = 86. - _Wolfdieter Lang_, Aug 11 2017

T. D. Noe, Rows n=0..50 of triangle, flattened
Priyavrat Deshpande, Krishna Menon, and Anurag Singh, A combinatorial statistic for labeled threshold graphs, arXiv:2103.03865 [math.CO], 2021.
Thomas Godland and Zakhar Kabluchko, Projections and angle sums of permutohedra and other polytopes, arXiv:2009.04186 [math.MG], 2020.
Thomas Godland and Zakhar Kabluchko, Projections and Angle Sums of Belt Polytopes and Permutohedra, Res. Math. (2023) Vol. 78, Art. No. 140.
Z. Kabluchko, V. Vysotsky, and D. Zaporozhets, Convex hulls of random walks, hyperplane arrangements, and Weyl chambers, arXiv preprint arXiv:1510.04073 [math.PR], 2015.
Wolfdieter Lang, On Generating functions of Diagonals Sequences of Sheffer and Riordan Number Triangles, arXiv:1708.01421 [math.NT], August 2017.
Bruce E. Sagan and Joshua P. Swanson, q-Stirling numbers in type B, arXiv:2205.14078 [math.CO], 2022.

A039757 is signed version.
Row sums: A000165.
Columns: A001147, A004041, A028339, A028340, A028341; A000012, A000290, A024196, A024197, A024198.
Diagonals: A000012, A000290(n+1), A024196(n+1), A024197(n+1), A024198(n+1).
A161198 is a scaled triangle version and A109692 is a transposed triangle version.
Central terms: A293318.
Cf. A286718, A002208(n+1)/A002209(n+1).

nmax:=8; for n from 0 to nmax do a(n, 0) := doublefactorial(2*n-1) od: for n from 0 to nmax do a(n, n) := 1 od: for n from 2 to nmax do for m from 1 to n-1 do a(n, m) := (2*n-1)*a(n-1, m) + a(n-1, m-1) od; od: seq(seq(a(n, m), m=0..n), n=0..nmax); # Johannes W. Meijer, Jun 08 2009, revised Nov 25 2012

T[n_, k_] := Sum[(-2)^(n-i) Binomial[i, k] StirlingS1[n, i], {i, k, n}] (* Woodhouse *)
Join[{1},Flatten[Table[CoefficientList[Expand[Times@@Table[x+i,{i,1,2n+1,2}]],x],{n,0,10}]]] (* Harvey P. Dale, Jan 29 2013 *)

Triangle T(n, k), read by rows, given by [1, 2, 3, 4, 5, 6, 7, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Feb 20 2005
T(n, k) = Sum_{i=k..n} (-2)^(n-i) * binomial(i, k) * s(n, i) where s(n, k) are signed Stirling numbers of the first kind. - Francis Woodhouse (fwoodhouse(AT)gmail.com), Nov 18 2005
G.f. of row polynomials in y: 1/(1-(x+x*y)/(1-2*x/(1-(3*x+x*y)/(1-4*x/(1-(5*x+x*y)/(1-6*x*y/(1-... (continued fraction). - Paul Barry, Feb 07 2009
T(n, m) = (2*n-1)*T(n-1,m) + T(n-1,m-1) with T(n, 0) = (2*n-1)!! and T(n, n) = 1. - Johannes W. Meijer, Jun 08 2009
From Wolfdieter Lang, May 09 2017: (Start)
E.g.f. of row polynomials in y: (1/sqrt(1-2*x))*exp(-y*log(sqrt(1-2*x))) = exp(-(1+y)*log(sqrt(1-2*x))) = 1/sqrt(1-2*x)^(1+y).
E.g.f. of column m sequence: (1/sqrt(1-2*x))* (-log(sqrt(1-2*x)))^m/m!. For the special Sheffer, also known as exponential Riordan array, see a comment above. (End)
Boas-Buck type recurrence for column sequence k: T(n, k) = (n!/(n - k)) * Sum_{p=k..n-1} 2^(n-1-p)*(1 + 2*k*beta(n-1-p))*T(p, k)/p!, for n > k >= 0, with input T(k, k) = 1, and beta(k) = A002208(k+1)/A002209(k+1). See a comment and references in A286718. - Wolfdieter Lang, Aug 09 2017

A046643 From square root of Riemann zeta function: form Dirichlet series Sum b_n/n^s whose square is zeta function; sequence gives numerator of b_n.

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 1, 5, 3, 1, 1, 3, 1, 1, 1, 35, 1, 3, 1, 3, 1, 1, 1, 5, 3, 1, 5, 3, 1, 1, 1, 63, 1, 1, 1, 9, 1, 1, 1, 5, 1, 1, 1, 3, 3, 1, 1, 35, 3, 3, 1, 3, 1, 5, 1, 5, 1, 1, 1, 3, 1, 1, 3, 231, 1, 1, 1, 3, 1, 1, 1, 15, 1, 1, 3, 3, 1, 1, 1, 35, 35, 1, 1, 3, 1, 1, 1, 5, 1, 3

Offset: 1

N. J. A. Sloane, Dec 11 1999

b(n) = A046643(n)/A046644(n) is multiplicative with b(p^n) = (2n-1)!!/2^n/n!. Dirichlet g.f. of A046643(n)/A046644(n) is sqrt(zeta(x)). - Christian G. Bower, May 16 2005

That is, b(p^n) = A001147(n) / (A000079(n)*A000142(n)) = A010050(n)/A000290(A000165(n)) = (2n)!/((2^n*n!)^2). - Antti Karttunen, Jul 08 2017

			b_1, b_2, ... = 1, 1/2, 1/2, 3/8, 1/2, 1/4, 1/2, 5/16, 3/8, 1/4, 1/2, 3/16, ...

Cf. A000079, A000165, A001147, A001790, A005187, A010050, A028234, A067029.

Cf. A046644, A046645.

b := proc(n) option remember; local c,i,t1; if n = 1 then 1 else c := 1; t1 := divisors(n);
for i from 2 to nops(t1)-1 do c := c-b(t1[ i ])*b(n/t1[ i ]); od; c/2; fi; end;

b[1] = 1; b[n_] := b[n] = (dn = Divisors[n]; c = 1;
Do[c = c - b[dn[[i]]]*b[n/dn[[i]]], {i, 2, Length[dn] - 1}]; c/2); a[n_] := Numerator[b[n]]; a /@ Range[90] (* Jean-François Alcover, Apr 04 2011, after Maple version *)

A046643perA046644(n) = { my(c=1); if(1==n,c,fordiv(n,d, if((d>1)&&(dA046643perA046644(d)*A046643perA046644(n/d)))); (c/2)); }
A046643(n) = numerator(A046643perA046644(n)); \\ After Maple-program, Antti Karttunen, Jul 08 2017

for(n=1, 100, print1(numerator(direuler(p=2, n, 1/(1-X)^(1/2))[n]), ", ")) \\ Vaclav Kotesovec, May 04 2025

(define (A046643 n) (if (= 1 n) n (* (A001790 (A067029 n)) (A046643 (A028234 n)))))
;; Or, after Christian G. Bower's May 16 2005 comment:
(definec (A046643perA046644 n) (if (= 1 n) n (* (/ (A010050 (A067029 n)) (A000290 (A000165 (A067029 n)))) (A046643perA046644 (A028234 n)))))
(define (A046643 n) (numerator (A046643perA046644 n)))
(define (A046644 n) (denominator (A046643perA046644 n)))
;; Antti Karttunen, Jul 08 2017

Sum_{b|d} b(d)b(n/d) = 1. Also b_{2^j} = A001790[ j ]/2^A005187[ j ].
From Antti Karttunen, Jul 08 2017: (Start)
Multiplicative with a(p^n) = A001790(n).
a(1) = 1; for n > 1, a(n) = A001790(A067029(n)) * a(A028234(n)).
(End)
Sum_{j=1..n} A046643(j)/A046644(j) ~ n / sqrt(Pi*log(n)) * (1 + (1 - gamma/2)/(2*log(n))), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, May 04 2025

A278990 Number of loopless linear chord diagrams with n chords.

Original entry on oeis.org

1, 0, 1, 5, 36, 329, 3655, 47844, 721315, 12310199, 234615096, 4939227215, 113836841041, 2850860253240, 77087063678521, 2238375706930349, 69466733978519340, 2294640596998068569, 80381887628910919255, 2976424482866702081004, 116160936719430292078411

Offset: 0

N. J. A. Sloane, Dec 07 2016

nonn

See the signed version of these numbers, A000806, for much more information about these numbers.
From Gus Wiseman, Feb 27 2019: (Start)
Also the number of 2-uniform set partitions of {1..2n} containing no two successive vertices in the same block. For example, the a(3) = 5 set partitions are:
  {{1,3},{2,5},{4,6}}
  {{1,4},{2,5},{3,6}}
  {{1,4},{2,6},{3,5}}
  {{1,5},{2,4},{3,6}}
  {{1,6},{2,4},{3,5}}
(End)
From Gus Wiseman, Jul 05 2020: (Start)
Also the number of permutations of the multiset {1,1,2,2,...,n,n} with no two consecutive terms equal and where the first i appears before the first j for i < j. For example, the a(3) = 5 permutations are the following.
  (1,2,3,1,2,3)
  (1,2,3,1,3,2)
  (1,2,3,2,1,3)
  (1,2,3,2,3,1)
  (1,2,1,3,2,3)
(End)

Seiichi Manyama, Table of n, a(n) for n = 0..404 (terms 0..200 from Gheorghe Coserea)
Dmitry Efimov, The hafnian of Toeplitz matrices of a special type, perfect matchings and Bessel polynomials, arXiv:1904.08651 [math.CO], 2019.
H. Eriksson and A. Martin, Enumeration of Carlitz multipermutations, arXiv:1702.04177 [math.CO], 2017.
E. Krasko, I. Labutin, and A. Omelchenko, Enumeration of labelled and unlabelled Hamiltonian Cycles in complete k-partite graphs, arXiv:1709.03218 [math.CO], 2017, Table 1.
E. Krasko and A. Omelchenko, Enumeration of Chord Diagrams without Loops and Parallel Chords, arXiv:1601.05073 [math.CO], 2016.
E. Krasko and A. Omelchenko, Enumeration of Chord Diagrams without Loops and Parallel Chords, The Electronic Journal of Combinatorics, 24(3) (2017), #P3.43.
Gus Wiseman, The a(4) = 36 loopless linear chord diagrams.
Donovan Young, Counting Bubbles in Linear Chord Diagrams, arXiv:2311.01569 [math.CO], 2023.
Donovan Young, Bubbles in Linear Chord Diagrams: Bridges and Crystallized Diagrams, arXiv:2408.17232 [math.CO], 2024.

Column k=0 of A079267.
Column k=2 of A293157.
Row n=2 of A322013.
Cf. A000110, A000699 (topologically connected 2-uniform), A000806, A001147 (2-uniform), A003436 (cyclical version), A005493, A170941, A190823 (distance 3+ version), A322402, A324011, A324172.
Anti-run compositions are A003242.
Separable partitions are A325534.
Cf. A007716, A292884, A333489.
Other sequences involving the multiset {1,1,2,2,...,n,n}: A001147, A007717, A020555, A094574, A316972.

[n le 2 select 2-n else (2*n-3)*Self(n-1) + Self(n-2): n in [1..30]]; // G. C. Greubel, Sep 26 2023

RecurrenceTable[{a[n]== (2n-1)a[n-1] +a[n-2], a[0]==1, a[1]==0}, a, {n,0,20}] (* Vaclav Kotesovec, Sep 15 2017 *)
FullSimplify[Table[-I*(BesselI[1/2+n,-1] BesselK[3/2,1] - BesselI[3/2,-1] BesselK[1/2+ n,1]), {n,0,20}]] (* Vaclav Kotesovec, Sep 15 2017 *)
Table[(2 n-1)!! Hypergeometric1F1[-n,-2 n,-2], {n,0,20}] (* Eric W. Weisstein, Nov 14 2018 *)
Table[Sqrt[2/Pi]/E ((-1)^n Pi BesselI[1/2+n,1] +BesselK[1/2+n,1]), {n,0,20}] // FunctionExpand // FullSimplify (* Eric W. Weisstein, Nov 14 2018 *)
twouniflin[{}]:={{}};twouniflin[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@twouniflin[Complement[set,s]]]/@Table[{i,j},{j,Select[set,#>i+1&]}];
Table[Length[twouniflin[Range[n]]],{n,0,14,2}] (* Gus Wiseman, Feb 27 2019 *)

seq(N) = {
  my(a = vector(N)); a[1] = 0; a[2] = 1;
  for (n = 3, N, a[n] = (2*n-1)*a[n-1] + a[n-2]);
  concat(1, a);
};
seq(20) \\ Gheorghe Coserea, Dec 09 2016

def A278990_list(prec):
    P. = PowerSeriesRing(QQ, prec)
    return P( exp(-1+sqrt(1-2*x))/sqrt(1-2*x) ).egf_to_ogf().list()
A278990_list(30) # G. C. Greubel, Sep 26 2023

From Gheorghe Coserea, Dec 09 2016: (Start)
D-finite with recurrence a(n) = (2*n-1)*a(n-1) + a(n-2), with a(0) = 1, a(1) = 0.
E.g.f. y satisfies: 0 = (1-2*x)*y'' - 3*y' - y.
a(n) - a(n-1) = A003436(n) for all n >= 2. (End)
From Vaclav Kotesovec, Sep 15 2017: (Start)
a(n) = sqrt(2)*exp(-1)*(BesselK(1/2 + n, 1)/sqrt(Pi) - i*sqrt(Pi)*BesselI(1/2 + n, -1)), where i is the imaginary unit.
a(n) ~ 2^(n+1/2) * n^n  / exp(n+1). (End)
a(n) = A114938(n)/n! - Gus Wiseman, Jul 05 2020 (from Alexander Burstein's formula at A114938).
From G. C. Greubel, Sep 26 2023: (Start)
a(n) = (-1)^n * (i/e)*Sqrt(2/Pi) * BesselK(n + 1/2, -1).
G.f.: sqrt(Pi/(2*x)) * exp(-(1+x)^2/(2*x)) * Erfi((1+x)/sqrt(2*x)).
E.g.f.: exp(-1 + sqrt(1-2*x))/sqrt(1-2*x). (End)

a(0)=1 prepended by Gheorghe Coserea, Dec 09 2016

A008956 Triangle of central factorial numbers |4^k t(2n+1,2n+1-2k)| read by rows (n>=0, k=0..n).

Original entry on oeis.org

1, 1, 1, 1, 10, 9, 1, 35, 259, 225, 1, 84, 1974, 12916, 11025, 1, 165, 8778, 172810, 1057221, 893025, 1, 286, 28743, 1234948, 21967231, 128816766, 108056025, 1, 455, 77077, 6092515, 230673443, 3841278805, 21878089479, 18261468225, 1, 680

Offset: 0

N. J. A. Sloane

The n-th row gives the coefficients in the expansion of Product_{i=0..n-1} (x+(2i+1)^2), highest powers first (see the discussion of central factorial numbers in A008955). - N. J. A. Sloane, Feb 01 2011
Descending row polynomials in x^2 evaluated at k generate odd coefficients of e.g.f. sin(arcsin(kt)/k): 1, x^2 - 1, 9x^4 - 10x^2 + 1, 225x^6 - 259x^4 + 34x^2 - 1, ... - Ralf Stephan, Jan 16 2005
From Johannes W. Meijer, Jun 18 2009: (Start)
We define (Pi/2)*Beta(n-1/2-z/2,n-1/2+z/2)/Beta(n-1/2,n-1/2) = (Pi/2)*Gamma(n-1/2-z/2)* Gamma(n-1/2+z/2)/Gamma(n-1/2)^2 = sum(BG2[2m,n]*z^(2m), m = 0..infinity) with Beta(z,w) the Beta function. Our definition leads to BG2[2m,1] = 2*beta(2m+1) and the recurrence relation BG2[2m,n] = BG2[2m,n-1] - BG2[2m-2,n-1]/(2*n-3)^2 for m = -2, -1, 0, 1, 2, .. and n = 2, 3, .. , with beta(m) = sum((-1)^k/(1+2*k)^m, k=0..infinity). We observe that beta(2m+1) = 0 for m = -1, -2, -3, .. .We found for the BG2[2*m,n] = sum((-1)^(k+n)*t2(n-1,k-1)* 2*beta(2*m-2*n+2*k+1),k=1..n)/((2*n-3)!!)^2 with the central factorial numbers t2(n,m) as defined above; see also the Maple program.
From the BG2 matrix and the closely related EG2 and ZG2 matrices, see A008955, we arrive at the LG2 matrix which is defined by LG2[2m-1,1] = 2*lambda(2*m) and the recurrence relation LG2[2*m-1,n] = LG2[2*m-3,n-1]/((2*n-3)*(2*n-1)) - (2*n-3)*LG2[2*m-1,n-1]/(2*n-1) for m = -2, -1, 0, 1, 2, .. and n = 2, 3, .. , with lambda(m) = (1-2^(-m))*zeta(m) with zeta(m) the Riemann zeta function. We found for the matrix coefficients LG2[2m-1,n] = sum((-1)^(k+1)* t2(n-1,k-1)*2*lambda(2*m-2*n+2*k)/((2*n-1)!!*(2*n-3)!!), k=1..n) and we see that the central factorial numbers t2(n,m) once again play a crucial role.
(End)

			Triangle begins:
[1]
[1, 1]
[1, 10, 9]
[1, 35, 259, 225]
[1, 84, 1974, 12916, 11025]
[1, 165, 8778, 172810, 1057221, 893025]
[1, 286, 28743, 1234948, 21967231, 128816766, 108056025]
[1, 455, 77077, 6092515, 230673443, 3841278805, 21878089479, 18261468225]
...

P. L. Butzer, M. Schmidt, E. L. Stark and L. Vogt, Central Factorial Numbers: Their main properties and some applications, Numerical Functional Analysis and Optimization, 10 (5&6), 419-488 (1989). [From Johannes W. Meijer, Jun 18 2009]
J. Riordan, Combinatorial Identities, Wiley, 1968, p. 217.

Reinhard Zumkeller, Rows n = 0..100 of triangle, flattened
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972, Chapter 23, pp. 811-812. [From _Johannes W. Meijer_, Jun 18 2009]
R. H. Boels, T. Hansen, String theory in target space, arXiv preprint arXiv:1402.6356 [hep-th], 2014.
T. L. Curtright, D. B. Fairlie, C. K. Zachos, A compact formula for rotations as spin matrix polynomials, arXiv preprint arXiv:1402.3541 [math-ph], 2014.
T. L. Curtright, T. S. Van Kortryk, On Rotations as Spin Matrix Polynomials, arXiv:1408.0767 [math-ph], 2014.
M. Eastwood and H. Goldschmidt, Zero-energy fields on complex projective space, arXiv preprint arXiv:1108.1602 [math.DG], 2011.
M. Eastwood, The X-ray transform on projective space. - From _N. J. A. Sloane_, Oct 22 2012

Cf. A008958.
Columns include A000447, A001823. Right-hand columns include A001818, A001824, A001825. Cf. A008955.
Appears in A160480 (Beta triangle), A160487 (Lambda triangle), A160479 (ZL(n) sequence), A161736, A002197 and A002198. - Johannes W. Meijer, Jun 18 2009
Cf. A162443 (BG1 matrix) and A162448 (LG1 matrix). - Johannes W. Meijer, Jul 06 2009
Cf. A001147.

a008956 n k = a008956_tabl !! n !! k
a008956_row n = a008956_tabl !! n
a008956_tabl = [1] : f [1] 1 1 where
   f xs u t = ys : f ys v (t * v) where
     ys = zipWith (+) (xs ++ [t^2]) ([0] ++ map (* u^2) (init xs) ++ [0])
     v = u + 2
-- Reinhard Zumkeller, Dec 24 2013

f:=n->mul(x+(2*i+1)^2,i=0..n-1);
for n from 0 to 12 do
t1:=eval(f(n)); t1d:=degree(t1);
t12:=y^t1d*subs(x=1/y,t1); t2:=seriestolist(series(t12,y,20));
lprint(t2);
od: # N. J. A. Sloane, Feb 01 2011
A008956 := proc(n,k) local i ; mul( x+2*i-2*n-1,i=1..2*n) ; expand(%) ; coeftayl(%,x=0,2*(n-k)) ; abs(%) ; end: for n from 0 to 10 do for k from 0 to n do printf("%a,",A008956(n,k)) ; od: od: # R. J. Mathar, May 29 2009
nmax:=7: for n from 0 to nmax do t2(n, 0):=1: t2(n, n):=(doublefactorial(2*n-1))^2 od: for n from 1 to nmax do for k from 1 to n-1 do t2(n, k) := (2*n-1)^2*t2(n-1, k-1)+t2(n-1, k) od: od: seq(seq(t2(n, k), k=0..n), n=0..nmax); # Johannes W. Meijer, Jun 18 2009, Revised Sep 16 2012

t[, 0] = 1; t[n, n_] := t[n, n] = ((2*n-1)!!)^2; t[n_, k_] := t[n, k] = (2*n-1)^2*t[n-1, k-1] + t[n-1, k]; Table[t[n, k], {n, 0, 8}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jan 07 2014, after Johannes W. Meijer *)

{T(n, k) = if( n<=0, k==0, (-1)^k * polcoeff( numerator( 2^(2*n -1) / sum(j=0, 2*n - 1, binomial( 2*n - 1, j) / (x + 2*n - 1 - 2*j))), 2*n - 2*k))}; /* Michael Somos, Feb 24 2003 */

Conjecture row sums: Sum_{k=0..n} T(n,k) = |A101927(n+1)|. - R. J. Mathar, May 29 2009

May be generated by the recurrence t2(n,k) = (2*n-1)^2*t2(n-1,k-1)+t2(n-1,k) with t2(n,0) = 1 and t2(n,n)=((2*n-1)!!)^2. - Johannes W. Meijer, Jun 18 2009

More terms from Vladeta Jovovic, Apr 16 2000

Edited by N. J. A. Sloane, Feb 01 2011

A099174 Triangle read by rows: coefficients of modified Hermite polynomials.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 0, 3, 0, 1, 3, 0, 6, 0, 1, 0, 15, 0, 10, 0, 1, 15, 0, 45, 0, 15, 0, 1, 0, 105, 0, 105, 0, 21, 0, 1, 105, 0, 420, 0, 210, 0, 28, 0, 1, 0, 945, 0, 1260, 0, 378, 0, 36, 0, 1, 945, 0, 4725, 0, 3150, 0, 630, 0, 45, 0, 1, 0, 10395, 0, 17325, 0, 6930, 0, 990, 0, 55, 0, 1

Offset: 0

Ralf Stephan, on a suggestion of Karol A. Penson, Oct 13 2004

Absolute values of A066325.
T(n,k) is the number of involutions of {1,2,...,n}, having k fixed points (0 <= k <= n). Example: T(4,2)=6 because we have 1243,1432,1324,4231,3214 and 2134. - Emeric Deutsch, Oct 14 2006
Riordan array [exp(x^2/2),x]. - Paul Barry, Nov 06 2008
Same as triangle of Bessel numbers of second kind, B(n,k) (see Cheon et al., 2013). - N. J. A. Sloane, Sep 03 2013
The modified Hermite polynomial h(n,x) (as in the Formula section) is the numerator of the rational function given by f(n,x) = x + (n-2)/f(n-1,x), where f(x,0) = 1. - Clark Kimberling, Oct 20 2014
Second lower diagonal T(n,n-2) equals positive triangular numbers A000217 \ {0}. - M. F. Hasler, Oct 23 2014
From James East, Aug 17 2015: (Start)
T(n,k) is the number of R-classes (equivalently, L-classes) in the D-class consisting of all rank k elements of the Brauer monoid of degree n.
For n < k with n == k (mod 2), T(n,k) is the rank (minimal size of a generating set) and idempotent rank (minimal size of an idempotent generating set) of the ideal consisting of all rank <= k elements of the Brauer monoid. (End)
This array provides the coefficients of a Laplace-dual sequence H(n,x) of the Dirac delta function, delta(x), and its derivatives, formed by taking the inverse Laplace transform of these modified Hermite polynomials. H(n,x) = h(n,D) delta(x) with h(n,x) as in the examples and the lowering and raising operators L = -x and R = -x + D = -x + d/dx such that L H(n,x) = n * H(n-1,x) and R H(n,x) = H(n+1,x). The e.g.f. is exp[t H(.,x)] = e^(t^2/2) e^(t D) delta(x) = e^(t^2/2) delta(x+t). - Tom Copeland, Oct 02 2016
Antidiagonals of this entry are rows of A001497. - Tom Copeland, Oct 04 2016
This triangle is the reverse of that in Table 2 on p. 7 of the Artioli et al. paper and Table 6.2 on p. 234 of Licciardi's thesis, with associations to the telephone numbers. - Tom Copeland, Jun 18 2018 and Jul 08 2018
See A344678 for connections to a Heisenberg-Weyl algebra of differential operators, matching and independent edge sets of the regular n-simplices with partially labeled vertices, and telephone switchboard scenarios. - Tom Copeland, Jun 02 2021

			h(0,x) = 1
h(1,x) = x
h(2,x) = x^2 + 1
h(3,x) = x^3 + 3*x
h(4,x) = x^4 + 6*x^2 + 3
h(5,x) = x^5 + 10*x^3 + 15*x
h(6,x) = x^6 + 15*x^4 + 45*x^2 + 15
From _Paul Barry_, Nov 06 2008: (Start)
Triangle begins
   1,
   0,  1,
   1,  0,  1,
   0,  3,  0,  1,
   3,  0,  6,  0,  1,
   0, 15,  0, 10,  0,  1,
  15,  0, 45,  0, 15,  0,  1
Production array starts
  0, 1,
  1, 0, 1,
  0, 2, 0, 1,
  0, 0, 3, 0, 1,
  0, 0, 0, 4, 0, 1,
  0, 0, 0, 0, 5, 0, 1 (End)

Alois P. Heinz, Rows n = 0..150, flattened
M. Artioli, G. Dattoli, S. Licciardi, and S. Pagnutti, Motzkin Numbers: an Operational Point of View, arXiv:1703.07262 [math.CO], 2017.
Paul Barry, Riordan array, orthogonal polynomials as moments, and Hankel transforms, arXiv:1102.0921 [math.CO], 2011.
G.-S. Cheon, J.-H. Jung and L. W. Shapiro, Generalized Bessel numbers and some combinatorial settings, Discrete Math., 313 (2013), 2127-2138.
T. Copeland, Juggling Zeros in the Matrix (Example II), 2020.
James East and Robert D. Gray, Diagram monoids and Graham-Houghton graphs: idempotents and generating sets of ideals, arXiv:1404.2359 [math.GR], 2014. See Theorem 8.4 and Table 7. - _James East_, Aug 17 2015
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.
Alexander Kreinin, Combinatorial Properties of Mills' Ratio, arXiv:1405.5852 [math.CO], 2014. See Table 2. - _N. J. A. Sloane_, May 29 2014
S. Licciardi, Umbral Calculus, a Different Mathematical Language, arXiv:1803.03108 [math.CA], 2018.
R. Paris, A uniform asymptotic expansion for the incomplete gamma function, Journal of Computational and Applied Mathematics, 148 (2002), p. 223-239 (See p. 329 and A137286. From _Tom Copeland_, Jan 03 2016).
R. Sazdanovic, A categorification of the polynomial ring, slide presentation, 2011
S. Yang and Z. Qiao, The Bessel numbers and Bessel matrices, Jrn. Math. Rsch. and Exposition, July 2011, Vol. 31, No. 4, pp.627-636. DOI:10.3770/j.issn:1000-341X.2011.04.006.

Row sums (polynomial values at x=1) are A000085.
Polynomial values: A005425 (x=2), A202834 (x=3), A202879(x=4).
Cf. A000217, A001147, A059343, A060821, A066325.
Cf. A000384, A014105, A034839, A049403, A096713, A100861, A104556, A122848, A130757, A176230, A176231.
Cf. A137286.
Cf. A001497.
Cf. A111062, A344678.

T:=proc(n,k) if n-k mod 2 = 0 then n!/2^((n-k)/2)/((n-k)/2)!/k! else 0 fi end: for n from 0 to 12 do seq(T(n,k),k=0..n) od; # yields sequence in triangular form; Emeric Deutsch, Oct 14 2006

nn=10;a=y x+x^2/2!;Range[0,nn]!CoefficientList[Series[Exp[a],{x,0,nn}],{x,y}]//Grid  (* Geoffrey Critzer, May 08 2012 *)
H[0, x_] = 1; H[1, x_] := x; H[n_, x_] := H[n, x] = x*H[n-1, x]-(n-1)* H[n-2, x]; Table[CoefficientList[H[n, x], x], {n, 0, 11}] // Flatten // Abs (* Jean-François Alcover, May 23 2016 *)
T[ n_, k_] := If[ n < 0, 0, Coefficient[HermiteH[n, x I/Sqrt[2]] (Sqrt[1/2]/I)^n, x, k]]; (* Michael Somos, May 10 2019 *)

T(n,k)=if(k<=n && k==Mod(n,2), n!/k!/(k=(n-k)/2)!>>k) \\ M. F. Hasler, Oct 23 2014

import sympy
from sympy import Poly
from sympy.abc import x, y
def H(n, x): return 1 if n==0 else x if n==1 else x*H(n - 1, x) - (n - 1)*H(n - 2, x)
def a(n): return [abs(cf) for cf in Poly(H(n, x), x).all_coeffs()[::-1]]
for n in range(21): print(a(n)) # Indranil Ghosh, May 26 2017

def Trow(n: int) -> list[int]:
    row: list[int] = [0] * (n + 1); row[n] = 1
    for k in range(n - 2, -1, -2):
        row[k] = (row[k + 2] * (k + 2) * (k + 1)) // (n - k)
    return row  # Peter Luschny, Jan 08 2023

def A099174_triangle(dim):
    M = matrix(ZZ,dim,dim)
    for n in (0..dim-1): M[n,n] = 1
    for n in (1..dim-1):
        for k in (0..n-1):
            M[n,k] = M[n-1,k-1]+(k+1)*M[n-1,k+1]
    return M
A099174_triangle(9)  # Peter Luschny, Oct 06 2012

h(k, x) = (-I/sqrt(2))^k * H(k, I*x/sqrt(2)), H(n, x) the Hermite polynomials (A060821, A059343).
T(n,k) = n!/(2^((n-k)/2)*((n-k)/2)!k!) if n-k >= 0 is even; 0 otherwise. - Emeric Deutsch, Oct 14 2006
G.f.: 1/(1-x*y-x^2/(1-x*y-2*x^2/(1-x*y-3*x^2/(1-x*y-4*x^2/(1-... (continued fraction). - Paul Barry, Apr 10 2009
E.g.f.: exp(y*x + x^2/2). - Geoffrey Critzer, May 08 2012
Recurrence: T(0,0)=1, T(0,k)=0 for k>0 and for n >= 1 T(n,k) = T(n-1,k-1) + (k+1)*T(n-1,k+1). - Peter Luschny, Oct 06 2012
T(n+2,n) = A000217(n+1), n >= 0. - M. F. Hasler, Oct 23 2014
The row polynomials P(n,x) = (a. + x)^n, umbrally evaluated with (a.)^n = a_n = aerated A001147, are an Appell sequence with dP(n,x)/dx = n * P(n-1,x). The umbral compositional inverses (cf. A001147) of these polynomials are given by the same polynomials signed, A066325. - Tom Copeland, Nov 15 2014
From Tom Copeland, Dec 13 2015: (Start)
The odd rows are (2x^2)^n x n! L(n,-1/(2x^2),1/2), and the even, (2x^2)^n n! L(n,-1/(2x^2),-1/2) in sequence with n= 0,1,2,... and L(n,x,a) = Sum_{k=0..n} binomial(n+a,k+a) (-x)^k/k!, the associated Laguerre polynomial of order a. The odd rows are related to A130757, and the even to A176230 and A176231. Other versions of this entry are A122848, A049403, A096713 and A104556, and reversed A100861, A144299, A111924. With each non-vanishing diagonal divided by its initial element A001147(n), this array becomes reversed, aerated A034839.
Create four shift and stretch matrices S1,S2,S3, and S4 with all elements zero except S1(2n,n) = 1 for n >= 1, S2(n,2n) = 1 for n >= 0, S3(2n+1,n) = 1 for n >= 1, and S4(n,2n+1) = 1 for n >= 0. Then this entry's lower triangular matrix is T = Id + S1 * (A176230-Id) * S2 + S3 * (unsigned A130757-Id) * S4 with Id the identity matrix. The sandwiched matrices have infinitesimal generators with the nonvanishing subdiagonals A000384(n>0) and A014105(n>0).
As an Appell sequence, the lowering and raising operators are L = D and R = x + dlog(exp(D^2/2))/dD = x + D, where D = d/dx, L h(n,x) = n h(n-1,x), and R h(n,x) = h(n+1,x), so R^n 1 = h(n,x). The fundamental moment sequence has the e.g.f. e^(t^2/2) with coefficients a(n) = aerated A001147, i.e., h(n,x) = (a. + x)^n, as noted above. The raising operator R as a matrix acting on o.g.f.s (formal power series) is the transpose of the production matrix P below, i.e., (1,x,x^2,...)(P^T)^n (1,0,0,...)^T = h(n,x).
For characterization as a Riordan array and associations to combinatorial structures, see the Barry link and the Yang and Qiao reference. For relations to projective modules, see the Sazdanovic link.
(End)
From the Appell formalism, e^(D^2/2) x^n = h_n(x), the n-th row polynomial listed below, and e^(-D^2/2) x^n = u_n(x), the n-th row polynomial of A066325. Then R = e^(D^2/2) * x * e^(-D^2/2) is another representation of the raising operator, implied by the umbral compositional inverse relation h_n(u.(x)) = x^n. - Tom Copeland, Oct 02 2016
h_n(x) = p_n(x-1), where p_n(x) are the polynomials of A111062, related to the telephone numbers A000085. - Tom Copeland, Jun 26 2018
From Tom Copeland, Jun 06 2021: (Start)
In the power basis x^n, the matrix infinitesimal generator M = A132440^2/2, when acting on a row vector for an o.g.f., is the matrix representation for the differential operator D^2/2.
e^{M} gives the coefficients of the Hermite polynomials of this entry.
The only nonvanishing subdiagonal of M, the second subdiagonal (1,3,6,10,...), gives, aside from the initial 0, the triangular numbers A000217, the number of edges of the n-dimensional simplices with (n+1) vertices. The perfect matchings of these simplices are the aerated odd double factorials A001147 noted above, the moments for the Hermite polynomials.
The polynomials are also generated from A036040 with x[1] = x, x[2] = 1, and the other indeterminates equal to zero. (End)

A123023 a(n) = (n-1)*a(n-2), a(0)=1, a(1)=0.

Original entry on oeis.org

1, 0, 1, 0, 3, 0, 15, 0, 105, 0, 945, 0, 10395, 0, 135135, 0, 2027025, 0, 34459425, 0, 654729075, 0, 13749310575, 0, 316234143225, 0, 7905853580625, 0, 213458046676875, 0, 6190283353629375, 0, 191898783962510625, 0, 6332659870762850625, 0, 221643095476699771875

Offset: 0

Roger L. Bagula, Sep 24 2006

a(n) is the number of ways of separating n terms into pairs. - Stephen Crowley, Apr 07 2007
a(n) is the n-th moment of the standard normal distribution. - Hal M. Switkay, Nov 06 2019
a(n) is the number of fixed-point free involutions in the symmetric group of degree n. - Nick Krempel, Feb 26 2020

			From _Gus Wiseman_, Dec 23 2018: (Start)
The a(6) = 15 ways of partitioning {1,2,3,4,5,6} into disjoint pairs:
  {{12}{34}{56}},  {{12}{35}{46}},  {{12}{36}{45}},
  {{13}{24}{56}},  {{13}{25}{46}},  {{13}{26}{45}},
  {{14}{23}{56}},  {{14}{25}{36}},  {{14}{26}{35}},
  {{15}{23}{46}},  {{15}{24}{36}},  {{15}{26}{34}},
  {{16}{23}{45}},  {{16}{24}{35}},  {{16}{25}{34}}.
(End)

Richard Bronson, Schaum's Outline of Modern Introductory Differential Equations, MacGraw-Hill, New York, 1973, page 107, solved problem 19.18
Norbert Wiener, Nonlinear Problems in Random Theory, 1958, Equation 1.31

G. C. Greubel, Table of n, a(n) for n = 0..799
Mihai Nica, Terence Tao's central limit theorem, Double Factorials (!!) and the Moment Method, YouTube video (2023).
Sebastian Volz, Design and Implementation of Efficient Algorithms for Operations on Partitions of Sets, Bachelor Thesis, Saarland Univ. (Germany, 2023). See p. 45.
Michael Wallner, A bijection of plane increasing trees with relaxed binary trees of right height at most one, arXiv:1706.07163 [math.CO], 2017-2018, Table 2 on p. 15.

Cf. A000085, A000110, A000124, A000142, A000587, A001147, A001464, A006129, A079267, A124794, A186021, A295193.

a:=[1,0]; [n le 2 select a[n] else (n-2)*Self(n-2): n in [1..30]]; // Marius A. Burtea, Nov 07 2019

with(combstruct): ZL2 := [S, {S=Set(Cycle(Z, card=2))}, labeled]:
seq(count(ZL2, size=n), n=0..36); # Zerinvary Lajos, Sep 24 2007
a := n -> ifelse(irem(n, 2) = 1, 0, 2^(n/2) * pochhammer(1/2, n/2)):
seq(a(n), n = 0..36); # Peter Luschny, Jan 11 2023

RecurrenceTable[{a[0] == 1, a[1] == 0, a[n] == (n - 1) a[n - 2]}, a[n], {n, 0, 31}] (* Ray Chandler, Jul 30 2015 *)

a(n) = (1/2)*Gamma((1/2)*n + 1/2)*2^((1/2)*n)*(1 + (-1)^n)/sqrt(Pi). - Stephen Crowley, Apr 07 2007
E.g.f.: exp(x^2/2). - Geoffrey Critzer, Mar 15 2009
a(2n) = A001147(n). - R. J. Mathar, Oct 11 2011
From Sergei N. Gladkovskii, Nov 18 2012, Dec 05 2012, May 16 2013, May 24 2013, Jun 07 2013: (Start)
Continued fractions:
E.g.f.: E(0) where E(k) = 1 + x^2*(4*k+1)/((4*k+2)*(4*k+3) - x^2*(4*k+2)*(4*k+3)^2/(x^2*(4*k+3) + (4*k+4)*(4*k+5)/E(k+1))).
G.f.: 1/G(0) where G(k) =  1 - x^2*(k+1)/G(k+1).
G.f.: 1 + x^2/(1+x) + Q(0)*x^3/(1+x), where Q(k) = 1 + (2*k+3)*x/(1 - x/(x + 1/Q(k+1))).
G.f.: G(0)/2, where G(k) = 1 + 1/(1-x/(x+1/x/(2*k+1)/G(k+1))).
G.f.: (G(0) - 1)*x/(1+x) + 1, where G(k) = 1 + x*(2*k+1)/(1 - x/(x + 1/G(k+1))). (End)
For n even, a(n) = A001147(n/2) = A124794(3^(n/2)). a(n) is also the coefficient of x1*...*xn in Product_{1 <= i < j <= n} (1 + xi*xj). - Gus Wiseman, Dec 23 2018
a(n) = 2^(n/2)*Pochhammer(1/2, n/2)*(n+1 mod 2). - Peter Luschny, Jan 11 2023

Edited by N. J. A. Sloane, Jan 06 2008
Better name by Sergei N. Gladkovskii, May 24 2013
Leading term 1 dropped, offset changed, and entry edited correspondingly by Andrey Zabolotskiy, Nov 07 2019

A002718 Number of bicoverings of an n-set.

Original entry on oeis.org

1, 0, 1, 8, 80, 1088, 19232, 424400, 11361786, 361058000, 13386003873, 570886397340, 27681861184474, 1511143062540976, 92091641176725504, 6219762391554815200, 462595509951068027741, 37676170944802047077248, 3343539821715571537772071, 321874499078487207168905840

Offset: 0

N. J. A. Sloane

Another description: number of proper 2-covers of [1,...,n].

			For n=3, there are 8 collections of distinct subsets of {1,2,3} with the property that each of 1, 2, and 3 appears in exactly two subsets:
  {1,2,3},{1,2},{3}
  {1,2,3},{1,3},{2}
  {1,2,3},{2,3},{1}
  {1,2,3},{1},{2},{3}
  {1,2},{1,3},{2,3}
  {1,2},{1,3},{2},{3}
  {1,2},{2,3},{1},{3}
  {1,3},{2,3},{1},{2}
Therefore a(3) = 8. - _Michael B. Porter_, Jul 16 2016

L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 303, #40.
I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, John Wiley and Sons, N.Y., 1983.
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).

Alois P. Heinz, Table of n, a(n) for n = 0..100
Peter Cameron, Thomas Prellberg, Dudley Stark, Asymptotic enumeration of 2-covers and line graphs, Discrete Math. 310 (2010), no. 2, 230-240 (see t_n).
L. Comtet, Birecouvrements et birevêtements d'un ensemble fini, Studia Sci. Math. Hungar 3 (1968): 137-152. [Annotated scanned copy. Warning: the table of v(n,k) has errors.]

Cf. A020554, A002719, A003462, A059945-A059951, A060053. Row sums of A059443.

nmax = 16; imax = 2*(nmax - 2); egf := E^(-x - 1/2*x^2*(E^y - 1))*Sum[(x^i/i!)*E^(Binomial[i, 2]*y), {i, 0, imax}]; fx = CoefficientList[Series[egf, {y, 0, imax}], y]*Range[0, imax]!; a[n_] := Drop[ CoefficientList[ Series[fx[[n + 1]], {x, 0, imax}], x], 3] // Total; Table[ a[n] , {n, 2, nmax}] (* Jean-François Alcover, Apr 04 2013 *)

E.g.f. for k-block bicoverings of an n-set is exp(-x-1/2*x^2*(exp(y)-1))*Sum_{i=0..inf} x^i/i!*exp(binomial(i, 2)*y).
Stirling_2 transform of A060053.
The e.g.f.'s of A002718 (T(x)) and A060053 (V(x)) are related by T(x) = V(e^x-1).
a(n) = Sum_{m=0..n + floor(n/2); k=0..n; s=0..min(m/2,k); t=0..m-2s} Stirling2(n,k) * k!/m! * binomial(m,2s) * A001147(s) * (-1)^(m+s+t) * binomial(m-2s,t) * binomial(t*(t-1)/2,k-s). Interpret m as the number of blocks in a bicovering, k the number of clumps of points that are always all together in blocks. This formula counts bicoverings by quotienting them to the clumpless case (an operation which preserves degree) and counting incidence matrices of those, and counts those matrices as the transposes of incidence matrices of labeled graphs with no isolated points and no isolated edges. - David Pasino, Jul 09 2016

More terms from Vladeta Jovovic, Feb 18 2001

a(0), a(1) prepended by Alois P. Heinz, Jul 29 2016

cp's OEIS Frontend

A053871 a(0)=1; a(1)=0; a(n) = 2*(n-1)*(a(n-1) + a(n-2)).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A060540 Square array read by antidiagonals downwards: T(n,k) = (n*k)!/(k!^n*n!), (n>=1, k>=1), the number of ways of dividing nk labeled items into n unlabeled boxes with k items in each box.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A112007 Coefficient triangle for polynomials used for o.g.f.s for unsigned Stirling1 diagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A028338 Triangle of coefficients in expansion of (x+1)*(x+3)*...*(x + 2n - 1) in rising powers of x.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A046643 From square root of Riemann zeta function: form Dirichlet series Sum b_n/n^s whose square is zeta function; sequence gives numerator of b_n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Scheme

Formula

A278990 Number of loopless linear chord diagrams with n chords.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

A053871 a(0)=1; a(1)=0; a(n) = 2(n-1)(a(n-1) + a(n-2)).

A060540 Square array read by antidiagonals downwards: T(n,k) = (nk)!/(k!^nn!), (n>=1, k>=1), the number of ways of dividing nk labeled items into n unlabeled boxes with k items in each box.

A028338 Triangle of coefficients in expansion of (x+1)(x+3)...*(x + 2n - 1) in rising powers of x.