A167594 -id:A167594 - cp's OEIS frontend

A002866 a(0) = 1; for n > 0, a(n) = 2^(n-1)*n!.

1, 1, 4, 24, 192, 1920, 23040, 322560, 5160960, 92897280, 1857945600, 40874803200, 980995276800, 25505877196800, 714164561510400, 21424936845312000, 685597979049984000, 23310331287699456000, 839171926357180416000, 31888533201572855808000, 1275541328062914232320000

Offset: 0

N. J. A. Sloane

Consider the set of n-1 odd numbers from 3 to 2n-1, i.e., {3, 5, ..., 2n-1}. There are 2^(n-1) subsets from {} to {3, 5, 7, ..., 2n-1}; a(n) = the sum of the products of terms of all the subsets. (Product for empty set = 1.) a(4) = 1 + 3 + 5 + 7 + 3*5 + 3*7 + 5*7 + 3*5*7 = 192. - Amarnath Murthy, Sep 06 2002
Also, a(n-1) is the number of ways to lace a shoe that has n pairs of eyelets such that there is a straight (horizontal) connection between all adjacent eyelet pairs. - Hugo Pfoertner, Jan 27 2003
This is also the denominator of the integral of ((1-x^2)^(n-1/2))/(Pi/4) where x ranges from 0 to 1. The numerator is (2*x)!/(x!*2^x). In both cases n starts at 1. E.g., the denominator when n=3 is 24 and the numerator is 15. - Al Hakanson (hawkuu(AT)excite.com), Oct 17 2003
Number of ways to use the elements of {1,...,n} once each to form a sequence of nonempty lists. - Bob Proctor, Apr 18 2005
Row sums of A131222. - Paul Barry, Jun 18 2007
Number of rotational symmetries of an n-cube. The number of all symmetries of an n-cube is A000165. See Egan for signed cycle notation, other notes, tables and animation. - Jonathan Vos Post, Nov 28 2007
1, 4, 24, 192, 1920, ... is the exponential (or binomial) convolution of 1, 1, 3, 15, 105, ... and 1, 3, 15, 105, 945 (A001147). - David Callan, Jul 25 2008
The n-th term of this sequence is the number of regions into which n-dimensional space is partitioned by the 2n hyperplanes of the form x_i=x_j and x_i=-x_j (for 1 <= i < j <= n). - Edward Scheinerman (ers(AT)jhu.edu), May 04 2008
a(n) is the number of ways to seat n churchgoers into pews and then linearly order the nonempty pews. - Geoffrey Critzer, Mar 16 2009
Equals the row sums of A156992. - Geoffrey Critzer, Mar 05 2010
From Gary W. Adamson, May 17 2010: (Start)
Next term in the series = (1, 3, 5, 7, ...) dot (1, 1, 4, 24, ...);
e.g., a(5) = 1920 = (1, 3, 5, 7, 9) dot (1, 1, 4, 24, 192) = (1 + 3 + 20 + 168 + 1728). (End)
a(n) is the number of ways to represent the permutations of {1,2,...,n} in cycle notation, taking into account that we can permute the order of all cycles and also have k ways to write a length-k cycle.
For positive n, a(n) equals the permanent of the n X n matrix with consecutive integers 1 to n along the main diagonal, consecutive integers 2 to n along the subdiagonal, and 1's everywhere else. - John M. Campbell, Jul 10 2011
From Dennis P. Walsh, Nov 26 2011: (Start)
Number of ways to arrange n books on consecutive bookshelves.
To derive a(n) = n!2^(n-1), we note that there are n! ways to arrange the books in a row. Then there are 2^(n-1) ways to place the arranged books on consecutive shelves since there are 2^(n-1) ordered partitions of n. Hence a(n) = n!2^(n-1).
Also, a(n) is the number of ways to stack n different alphabet blocks in contiguous stacks.
Furthermore, a(n) is the number of labeled, rooted forests that have (i) each root labeled larger than any nonroot, (ii) each root having exactly one child node, (iii) n non-root nodes, and (iv) each node in the forest with at most one child node.
Example: a(3)=24 since there are 24 arrangements of books b1, b2, and b3 on consecutive shelves, namely, |b1 b2 b3|, |b1 b3 b2|, |b2 b1 b3|, |b2 b3 b1|, |b3 b1 b2|, |b3 b2 b1|, |b1 b2||b3|, |b2 b1| |b3|, |b1 b3||b2|, |b3 b1||b2|, |b2 b3||b1|, |b3 b2||b1|, |b1||b2 b3|,|b1||b3 b2|, |b2||b1 b3|, |b2||b3 b1|, |b3||b1 b2|, |b3||b2 b1|, |b1||b2||b3|, |b1||b3||b2|, |b2||b1||b3|, |b2||b3||b1|, |b3||b1||b2|, and |b3||b2||b1|.
(End)
For n > 3, a(n) is the order of the Coxeter group (also called Weyl group) of type D_n. - Tom Edgar, Nov 05 2013

			For the shoe lacing: with the notation introduced in A078602 the a(3-1) = 4 "straight" lacings for 3 pairs of eyelets are: 125346, 125436, 134526, 143526. Their mirror images 134256, 143256, 152346, 152436 are not counted.
a(3) = 24 because the 24 rotations of a three-dimensional cube fall into four distinct classes:
(i) the identity, which leaves everything fixed;
(ii) 9 rotations which leave the centers of two faces fixed, comprising rotations of 90, 180 and 270 degrees for each of 3 pairs of faces;
(iii) 6 rotations which leave the centers of two edges fixed, comprising rotations of 180 degrees for each of 6 pairs of edges;
(iv) 8 rotations which leave two vertices fixed, comprising rotations of 120 and 240 degrees for each of 4 pairs of vertices. For an n-cube, rotations can be more complex. For example, in 4 dimensions a rotation can either act in a single plane, such as the x-y plane, while leaving any vectors orthogonal to that plane unchanged, or it can act in two orthogonal planes, performing rotations in both and leaving no vectors fixed. In higher dimensions, there will be room for more planes and more choices as to the number of planes in which a given rotation acts.

N. Bourbaki, Groupes et alg. de Lie, Chap. 4, 5, 6, p. 257.
A. P. Prudnikov, Yu. A. Brychkov and O.I. Marichev, "Integrals and Series", Volume 1: "Elementary Functions", Chapter 4: "Finite Sums", New York, Gordon and Breach Science Publishers, 1986-1992, Eq. (4.2.2.26)
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
Ofek Lauber Bonomo and Shlomi Reuveni, Occupancy correlations in the asymmetric simple inclusion process, arXiv:1905.02170 [cond-mat.stat-mech], 2019.
J.-P. Bultel and A. Chouria, J.-G. Luque and O. Mallet, Word symmetric functions and the Redfield-Polya theorem, Accepted to FPSAC 2013, 2013; arXiv:1302.5815 [math.CO], 2013.
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. 3 (2000), #00.1.5.
Greg Egan, Hypercube Mathematical Details, 2007-2008.
Joël Gay and Vincent Pilaud, The weak order on Weyl posets, arXiv:1804.06572 [math.CO], 2018.
W. S. Gray and M. Thitsa, System Interconnections and Combinatorial Integer Sequences, in: 2013 45th Southeastern Symposium on System Theory (SSST).
Guo-Niu Han, Enumeration of Standard Puzzles, 2011. [Cached copy]
Guo-Niu Han, Enumeration of Standard Puzzles, arXiv:2006.14070 [math.CO], 2020.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 121.
Victor Meally, Comparison of several sequences given in Motzkin's paper "Sorting numbers for cylinders...", letter to N. J. A. Sloane, N. D.
T. S. Motzkin, Sorting numbers for cylinders and other classification numbers, in Combinatorics, Proc. Symp. Pure Math. 19, AMS, 1971, pp. 167-176. [Annotated, scanned copy]
Norihiro Nakashima and Shuhei Tsujie, Enumeration of Flats of the Extended Catalan and Shi Arrangements with Species, arXiv:1904.09748 [math.CO], 2019.
J.-C. Novelli and J.-Y. Thibon, Hopf Algebras of m-permutations,(m+1)-ary trees, and m-parking functions, arXiv:1403.5962 [math.CO], 2014.
OEIS Wiki, Sorting numbers.
Hugo Pfoertner, Counting straight shoe lacings. FORTRAN program and results.
Robert A. Proctor, Let's Expand Rota's Twelvefold Way For Counting Partitions!, arXiv:math/0606404 [math.CO], 2007.
N. J. A. Sloane and Thomas Wieder, The Number of Hierarchical Orderings, arXiv:math/0307064 [math.CO], 2003.
N. J. A. Sloane and Thomas Wieder, The Number of Hierarchical Orderings, Order 21 (2004), 83-89.
Index to divisibility sequences
Index entries for sequences related to shoe lacings
Index entries for related partition-counting sequences

Cf. A028371, A078602, A078698, A078702, A156992.
Bisections give A002671 and A274304.
Appears in A167584 (n >= 1); equals the row sums of A167594 (n >= 1). - Johannes W. Meijer, Nov 12 2009

FORTRAN
```
See Pfoertner link.
```

[1] cat [2^(n-1)*Factorial(n): n in [1..25]]; // G. C. Greubel, Jun 13 2019

A002866 := n-> `if`(n=0,1,2^(n-1)*n!):
with(combstruct); SeqSeqL := [S, {S=Sequence(U,card >= 1), U=Sequence(Z,card >=1)},labeled];
seq(ceil(count(Subset(n))*count(Permutation(n))/2),n=0..17); # Zerinvary Lajos, Oct 16 2006
G(x):=(1-x)/(1-2*x): f[0]:=G(x): for n from 1 to 26 do f[n]:=diff(f[n-1],x) od:x:=0:seq(f[n],n=0..17); # Zerinvary Lajos, Apr 04 2009

Join[{1},Table[2^(n-1) n!,{n,25}]] (* Harvey P. Dale, Sep 27 2013 *)
a[n_] := (-1)^n Hypergeometric2F1Regularized[1, -n, 2 - n, 2];
Table[a[n], {n, 0, 20}]  (* Peter Luschny, Apr 26 2024 *)

a(n)=if(n,n!<<(n-1),1) \\ Charles R Greathouse IV, Jan 13 2012

a(n) = if(n == 0, 1, 2^(n-1)*n!);
vector(25, n, a(n-1)) \\ Altug Alkan, Oct 18 2015

[1] + [2^(n-1)*factorial(n) for n in (1..25)] # G. C. Greubel, Jun 13 2019

E.g.f.: (1 - x)/(1 - 2*x). - Paul Barry, May 26 2003, corrected Jun 18 2007
a(n) = n! * A011782(n).
For n >= 1, a(n) = Sum_{i=0..m/2} (-1)^i * binomial(n, i) * (n-2*i)^n. - Yong Kong (ykong(AT)curagen.com), Dec 28 2000
a(n) ~ 2^(1/2) * Pi^(1/2) * n^(3/2) * 2^n * e^(-n) * n^n*{1 + 13/12*n^(-1) + ...}. - Joe Keane (jgk(AT)jgk.org), Nov 23 2001
E.g.f. is B(A(x)), where B(x) = 1/(1 - x) and A(x) = x/(1 - x). - Geoffrey Critzer, Mar 16 2009
a(n) = Sum_{k=1..n} A156992(n,k). - Dennis P. Walsh, Nov 26 2011
a(n+1) = Sum_{k=0..n} A132393(n,k)*2^(n+k), n>0. - Philippe Deléham, Nov 28 2011
G.f.: 1 + x/(1 - 4*x/(1 - 2*x/(1 - 6*x/(1 - 4*x/(1 - 8*x/(1 - 6*x/(1 - 10*x/(1 - ... (continued fraction). - Philippe Deléham, Nov 29 2011
a(n) = 2*n*a(n-1) for n >= 2. - Dennis P. Walsh, Nov 29 2011
G.f.: (1 + 1/G(0))/2, where G(k) = 1 + 2*x*k - 2*x*(k + 1)/G(k+1); (continued fraction, Euler's 1st kind, 1-step). - Sergei N. Gladkovskii, Aug 02 2012
G.f.: 1 + x/Q(0), m=4, where Q(k) = 1 - m*x*(2*k + 1) - m*x^2*(2*k + 1)*(2*k + 2)/(1 - m*x*(2*k + 2) - m*x^2*(2*k + 2)*(2*k + 3)/Q(k+1)) ; (continued fraction). - Sergei N. Gladkovskii, Sep 23 2013
G.f.: 1 + x/(G(0) - x), where G(k) = 1 + x*(k+1) - 4*x*(k + 1)/(1 - x*(k + 2)/G(k+1)); (continued fraction). - Sergei N. Gladkovskii, Dec 24 2013
a(n) = Sum_{k=0..n} L(n,k)*k!; L(n,k) are the unsigned Lah numbers. - Peter Luschny, Oct 18 2014
a(n) = round(Sum_{k >= 1} log(k)^n/k^(3/2))/4, for n >= 1, which is related to the n-th derivative of zeta(x) evaluated at x = 3/2. - Richard R. Forberg, Jan 02 2015
a(n) = n!*hypergeom([-n+1], [], -1) for n>=1. - Peter Luschny, Apr 08 2015
From Amiram Eldar, Aug 04 2020: (Start)
Sum_{n >= 0} 1/a(n) = 2*sqrt(e) - 1.
Sum_{n >= 0} (-1)^n/a(n) = 2/sqrt(e) - 1. (End)

A024199 a(n) = (2n-1)!! * Sum_{k=0..n-1}(-1)^k/(2k+1).

Original entry on oeis.org

0, 1, 2, 13, 76, 789, 7734, 110937, 1528920, 28018665, 497895210, 11110528485, 241792844580, 6361055257725, 163842638377950, 4964894559637425, 147721447995130800, 5066706567801827025, 171002070002301095250, 6548719685561840296125, 247199273204273879989500

Offset: 0

Clark Kimberling

(2*n + 1)!!/a(n+1), n>=0, is the n-th approximant for William Brouncker's continued fraction for 4/Pi = 1 + 1^2/(2 + 3^2/(2 + 5^2/(2 + ... ))) See the C. Brezinski and J.-P. Delahaye references given under A142969 and A142970, respectively. The double factorials (2*n + 1)!! = A001147(n+1) enter. - Wolfdieter Lang, Oct 06 2008

			a(3) = (2*3 - 1)!! * Sum_{k=0..2} (-1)^k/(2k + 1) = 5!! * (1/(2*0 + 1) - 1/(2*1 + 1) + 1/(2*2 + 1)) = 5*3*1*(1/1 - 1/3 + 1/5) = 15 - 5 + 3 = 13. Notice that the first factor always cancels the common denominator of the sum. - _Michael B. Porter_, Jul 22 2016

A. E. Jolliffe, Continued Fractions, in Encyclopaedia Britannica, 11th ed., pp. 30-33; see p. 31.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Peter Luschny, Is the Gamma function misdefined? Hadamard versus Euler - Who found the better Gamma function?
Wikipedia, William Brouncker, 2nd Viscount Brouncker

Cf. A004041, A000407.
From Johannes W. Meijer, Nov 12 2009: (Start)
Cf. A007509 and A025547.
Equals first column of A167584.
Equals row sums of A167591.
Equals first right hand column of A167594.
(End)
Cf. A167576 and A135457.

[0] cat [ n le 2 select (n) else 2*Self(n-1)+(2*n-3)^2*Self(n-2): n in [1..25] ]; // Vincenzo Librandi, Feb 17 2015

a := proc(n) option remember; if n=0 then 0 elif n=1 then 1 else 2*a(n-1)+(2*n-3)^2* a(n-2) fi end: seq(a(n), n=0..20); # Peter Luschny, Nov 16 2016 after N. J. A. Sloane

f[k_] := (2 k - 1) (-1)^(k + 1)
t[n_] := Table[f[k], {k, 1, n}]
a[n_] := SymmetricPolynomial[n - 1, t[n]]
Table[a[n], {n, 1, 22}]    (* A024199 signed *)
(* Clark Kimberling, Dec 30 2011 *)
RecurrenceTable[{a[n+1] == 2*a[n] + (2*n-1)^2*a[n-1],a[0] == 0, a[1] == 1},a,{n,0,20}] (* Vaclav Kotesovec, Mar 18 2014 *)
CoefficientList[Series[Pi/4/Sqrt[1-2*x] - 1/2*Log[2*x+Sqrt[4*x^2-1]]/Sqrt[2*x-1], {x, 0, 20}], x] * Range[0, 20]! (* Vaclav Kotesovec, Mar 18 2014 *)

a(n) = s(1)s(2)...s(n)(1/s(1) - 1/s(2) + ... + c/s(n)) where c=(-1)^(n+1) and s(k) = 2k-1 for k = 1, 2, 3, ... (was previous definition). - Clark Kimberling
D-finite with recurrence a(0) = 0, a(1) = 1, a(n+1) = 2*a(n) + (2*n-1)^2*a(n-1). - N. J. A. Sloane, Jul 19 2002
a(n) + A024200(n) = A001147(n) = (2n-1)!!. - Max Alekseyev, Sep 23 2007
a(n)/A024200(n) -> Pi/(4-Pi) as n -> oo. - Max Alekseyev, Sep 23 2007
From Wolfdieter Lang, Oct 06 2008: (Start)
E.g.f. for a(n+1), n>=0: (sqrt(1-2*x)+arcsin(2*x)*sqrt(1+2*x)/2)/((1-4*x^2)^(1/2)*(1-2*x)). From the recurrence, solving (1-4*x^2)y''(x)-2*(8*x+1)*y'(x)-9*y=0 with inputs y(0)=1, y'(0)=2.
a(n+1) = A003148(n) + A143165(n), n>=0 (from the two terms of the e.g.f.). (End)
From Johannes W. Meijer, Nov 12 2009: (Start)
a(n) = (-1)^(n-1)*(2*n-3)!! + (2*n-1)*a(n-1) with a(0) = 0.
a(n) = (2*n-1)!!*sum((-1)^(k)/(2*k+1), k=0..n-1)
(End)
E.g.f.: Pi/4/sqrt(1-2*x) - 1/2*log(2*x+sqrt(4*x^2-1))/sqrt(2*x-1). - Vaclav Kotesovec, Mar 18 2014
a(n) ~ Pi * 2^(n-3/2) * n^n / exp(n). - Vaclav Kotesovec, Mar 18 2014
a(n) = (2*H(n+1/2)-Gamma(n+1/2))*2^(n-2)*sqrt(Pi) with H(x) the Hadamard factorial (see the link section). - Cyril Damamme, Jul 19 2015
a(n) = A135457(n) - (-1)^n A001147(n-1). - Cyril Damamme, Jul 19 2015
a(n) = (Pi + (-1)^n*(Psi(n/2 + 1/4) - Psi(n/2 + 3/4)))*Gamma(n+1/2)*2^(n-2)/sqrt(Pi). - Robert Israel, Jul 20 2015
a(n) = A167576(n) - A135457(n). - Cyril Damamme, Jul 22 2015
a(n)/A001147(n) -> Pi/4 as n -> oo. - Daniel Suteu, Jul 21 2016
From Peter Bala, Nov 15 2016: (Start)
Conjecture: a(n) = 1/2*Sum_{k = 0..2*n-1} (-1)^(n-k+1)*k!*(2*n - 2*k - 3)!!, where the double factorial of an odd integer (positive or negative) may be defined in terms of the gamma function as (2*N - 1)!! = 2^((N+1)/2)*Gamma(N/2 + 1)/sqrt(Pi).
E.g.f. 1/2*arcsin(2*x)/sqrt(1 - 2*x) = x + 2*x^2/2! + 13*x^3/3! + 76*x^4/4! + .... (End)

Edited by N. J. A. Sloane, Jul 19 2002

New name from Cyril Damamme, Jul 19 2015

A167584 The ED4 array read by antidiagonals.

Original entry on oeis.org

1, 2, 1, 13, 6, 1, 76, 41, 10, 1, 789, 372, 93, 14, 1, 7734, 4077, 1020, 169, 18, 1, 110937, 53106, 13269, 2212, 269, 22, 1, 1528920, 795645, 198990, 33165, 4140, 393, 26, 1, 28018665, 13536360, 3383145, 563850, 70485, 6996, 541, 30, 1

Offset: 1

Johannes W. Meijer, Nov 10 2009

The coefficients in the upper right triangle of the ED4 array (m>n) were found with the a(n,m) formula while the coefficients in the lower left triangle of the ED4 array (m<=n) were found with the recurrence relation, see below. We use for the array rows the letter n (>=1) and for the array columns the letter m (>=1).
For the ED1, ED2 and ED3 arrays see A167546, A167560 and A167572.
The Madhava-Gregory-Leibniz series representation for Pi/4 is the case m = 0 of the following more general result: for m = 0,1,2,... there holds 1/(2*m)! * Pi/4 = Sum_{k >= 0} ( (-1)^(m+k) * 1/Product_{j = -m .. m} (2*k + 1 + 2*j) ). The entries of this table are given by truncating these series to n-1 terms and then scaling by certain double factorials -- see the formula below. - Peter Bala, Nov 06 2016

			The ED4 array begins with:
  1, 1, 1, 1, 1, 1, 1, 1, 1, 1
  2, 6, 10, 14, 18, 22, 26, 30, 34, 38
  13, 41, 93, 169, 269, 393, 541, 713, 909, 1129
  76, 372, 1020, 2212, 4140, 6996, 10972, 16260, 23052, 31540
  789, 4077, 13269, 33165, 70485, 133869, 233877, 382989, 595605, 888045
  7734, 53106, 198990, 563850, 1339110, 2812194, 5389566, 9619770, 16216470, 26081490
  ...
From _Peter Bala_, Nov 06 2016: (Start)
Table extended to nonpositive values of m:
  n\m|     -4     -3    -2    -1    0
  -----------------------------------
   0 |      0      0     0     0    0
   1 |      1      1     1     1    1
   2 |    -18    -14   -10    -6   -2
   3 |    233    141    73    29    9
   4 |  -2844  -1428  -620  -228  -60
   5 |  39309  17877  7149  2325  525
  ...
Column  0: (-1)^(n+1)*(2*n - 3)!!*n. See A001193;
Column -1: (-1)^n*(2*n - 5)!!/3!!*n*(7 - 4*n^2);
Column -2: (-1)^n*(2*n - 7)!!/5!!*n(-149 + 120*n^2 - 16*n^4);
Column -3: (-1)^n*(2*n - 9)!!/7!!*n*(6483 - 6076*n^2 + 1232*n^4 - 64*n^6);
Column -4: (-1)^n*(2*n - 11)!!/9!!*n*(-477801 + 489136*n^2 - 120288*n^4 + 9984*n^6 - 256*n^8). (End)

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
Johannes W. Meijer, The four Escher-Droste arrays, jpg image, Mar 08 2013.
Wikipedia, Double factorial

A000012, A016825, A167585, A167586 and A167587 equal the first five rows of the array.
A024199, A167588 and A167589 equal the first three columns of the array.
A167590 equals the row sums of the ED4 array read by antidiagonals.
A167591 is a triangle related to the a(n) formulas of the rows of the ED4 array.
A167594 is a triangle related to the GF(z) formulas of the rows of the ED4 array.
Cf. A002866 (the 2^(n-1)*n! factor).
Cf. A167546 (ED1 array), A167560 (ED2 array), A167572 (ED3 array). Cf. A001193, A003881.

T := proc (n, m) option remember;
      if n = 0 then 0
       elif n = 1 then 1
       else (4*m-2)*T(n-1,m)+(2*n+2*m-5)*(2*n-2*m-1)*T(n-2,m)
      end if;
     end proc:
#square array read by antidiagonals
seq(seq(T(n-m,m), m = 1..n-1), n = 1..10);
# Peter Bala, Nov 06 2016

T[0, k_] := 0; T[1, k_] := 1; T[n_, k_] := T[n, k] = (4*k - 2)*T[n - 1, k] + (2*n + 2*k - 5)*(2*n - 2*k - 1)*T[n - 2, k]; Table[T[n - k, k], {n, 2, 12}, {k, 1, n - 1}] (* G. C. Greubel, Jan 20 2017 *)

a(n,m) = ((2*m-3)!!/(2*(2*m-2*n-3)!!))*Integral_{y=0..oo} sinh(y*(2*n))/(cosh(y))^(2*m-1) dy for m>n.
The (n-1)-differences of the n-th array row lead to the recurrence relation
Sum_{k=0..n-1} (-1)^k*binomial(n-1,k)*a(n,m-k) = 2^(n-1)*n!
From Peter Bala, Nov 06 2016: (Start)
T(n,m) = ((2*m - 3)!!/(2*(2*m - 2*n - 3)!!)) * Sum_{k = 0..n-1} (-1)^(k+1)*binomial(2*n - k - 1, k)*2^(2*n - 2*k - 1)*1/(2*n - 2*m - 2*k + 1), for n and m >= 0.
Note the double factorial for a negative odd integer N is defined in terms of the gamma function as N!! = 2^((N+1)/2)*Gamma(N/2 + 1)/sqrt(Pi).
T(n, m) = (2*m - 3)!! * (2*n + 2*m - 3)!! * Sum_{k = 0..n-1} ( (-1)^(m + k + 1) / Product_{j = -(m-1) .. m-1} (2*k + 1 + 2*j) ).
Using this result we can extend the table to nonpositive values of m (the column index). Column 0 is a signed version of A001193. We have for m <= 0, T(n,m) = (2*n - 2*|m| - 3)!!/(2*|m| + 1)!! * Sum_{k = 0..n-1} (-1)^k*Product_{j = -|m|..|m|} (2*k + 1 + 2*j).
Recurrence: T(n, m) = (4*m - 2)*T(n-1, m) + (2*n + 2*m - 5)*(2*n - 2*m - 1)*T(n-2, m).
For a fixed value of n, the entries in row n are polynomial in the value of the column index m. The first few polynomials are [1, 4*m - 2, 12*m^2 - 8*m + 9, 32*m^3 - 16*m^2 + 120*m - 60, 80*m^4 + 952*m^2 - 768*m + 525, ...]. (End)

A001193 a(n) = (n+1)(2n)!/(2^nn!) = (n+1)(2n-1)!!.

Original entry on oeis.org

1, 2, 9, 60, 525, 5670, 72765, 1081080, 18243225, 344594250, 7202019825, 164991726900, 4111043861925, 110681950128750, 3201870700153125, 99044533658070000, 3262279327362680625, 113987877673731311250, 4211218814057295665625, 164015890652757831187500

Offset: 0

N. J. A. Sloane

Solution to y' = A(x), y(0) = 0 satisfies 0 = x^2 + 2*y^2*x - y^2, where A(x) = e.g.f. - Michael Somos, Mar 11 2004

L. Comtet, Advanced Combinatorics, Reidel, 1974, pp. 166-167.
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

From Johannes W. Meijer, Nov 12 2009: (Start)
Equals the first right hand column of A167591.
Equals the first left hand column of A167594. (End)
Cf. A059366.

f:= x-> x/sqrt(1-2*x): a:= n-> subs(x=0, (D@@(n+1))(f)(x)):
seq(a(n), n=0..20); # Zerinvary Lajos, Apr 04 2009
# second Maple program:
a:= n-> (n+1)*doublefactorial(2*n-1):
seq(a(n), n=0..23);  # Alois P. Heinz, May 13 2020

Table[(n+1) (2*n-1)!!, {n,0,20}] (* Vladimir Joseph Stephan Orlovsky, Apr 14 2011 *)

PARI
```
a(n)=if(n<0, 0, (n+1)*(2*n)!/(2^n*n!))
```

E.g.f.: (1-x)/(1-2*x)^(3/2) = d/dx (x/(1-2*x)^(1/2)).
a(n) = uppermost term in the vector (M(T))^n * [1,0,0,0,...], where T = Transpose and M = the production matrix:
  1, 2;
  1, 2, 3;
  1, 2, 3, 4;
  1, 2, 3, 4, 5;
  ...
- Gary W. Adamson, Jul 08 2011
G.f.: A(x) = 1 + 2*x/(G(0) - 2*x) ; G(k) = 1 + k + x*(k+2)*(2*k+1) - x*(k+1)*(k+3)*(2*k+3)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Dec 06 2011
G.f.: U(0)/2  where U(k) = 1 + (2*k+1)/(1 - x/(x + 1/U(k+1))) (continued fraction). - Sergei N. Gladkovskii, Sep 25 2012
From Peter Bala, Nov 07 2016 and May 14 2020: (Start)
a(n) = (n + 1)*(2*n - 1)/n * a(n-1) with a(0) = 1.
a(n) = 2*a(n-1) + (2*n - 3)*(2*n + 1)*a(n-2) with a(0) = 1, a(1) = 2.
(End)
a(n) = A059366(n+1, n) = A059366(n+1, 1). - Petros Hadjicostas, May 13 2020

Better description from Wouter Meeussen, Mar 08 2001

More terms from James Sellers, May 01 2000

cp's OEIS Frontend

A002866 a(0) = 1; for n > 0, a(n) = 2^(n-1)*n!.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

FORTRAN

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A024199 a(n) = (2n-1)!! * Sum_{k=0..n-1}(-1)^k/(2k+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Formula

Extensions

A167584 The ED4 array read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A001193 a(n) = (n+1)*(2*n)!/(2^n*n!) = (n+1)*(2n-1)!!.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A001193 a(n) = (n+1)(2n)!/(2^nn!) = (n+1)(2n-1)!!.