A000165 -id:A000165 - cp's OEIS frontend

A047055 Quintuple factorial numbers: a(n) = Product_{k=0..n-1} (5*k + 2).

1, 2, 14, 168, 2856, 62832, 1696464, 54286848, 2008613376, 84361761792, 3965002804224, 206180145819648, 11752268311719936, 728640635326636032, 48818922566884614144, 3514962424815692218368, 270652106710808300814336, 22193472750286280666775552

Offset: 0

Joe Keane (jgk(AT)jgk.org)

Hankel transform is A169621. - Paul Barry, Dec 03 2009

G. C. Greubel, Table of n, a(n) for n = 0..350

Cf. A000165, A008544, A001813, A047657, A084947, A084948, A084949.

Cf. A052562, A008548, A047056, A169621.

List([0..20], n-> Product([0..n-1], k-> (5*k+2) )); # G. C. Greubel, Aug 17 2019

[1] cat [(&*[(5*k+2): k in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Aug 17 2019

a := n->product(5*i+2,i=0..n-1); [seq(a(j),j=0..30)];

Table[5^n*Pochhammer[2/5, n], {n,0,20}] (* G. C. Greubel, Aug 17 2019 *)
Join[{1},FoldList[Times,5*Range[0,20]+2]] (* Harvey P. Dale, Apr 03 2025 *)

vector(20, n, n--; prod(k=0,n-1, 5*k+2)) \\ G. C. Greubel, Aug 17 2019

[product((5*k+2) for k in (0..n-1)) for n in (0..20)] # G. C. Greubel, Aug 17 2019

E.g.f. (1-5*x)^(-2/5).
a(n) ~ sqrt(2*Pi)/Gamma(2/5)*n^(-1/10)*(5n/e)^n*(1 - (11/300)/n - ...). - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
a(n) = A084940(n)/A000142(n)*A000079(n) = 5^n*Pochhammer(2/5, n) = 5^n*Gamma(n+2/5)*sin(2*Pi/5)*Gamma(3/5)/Pi. - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
G.f.: 1/(1-2x/(1-5x/(1-7x/(1-10x/(1-12x/(1-15x/(1-17x/(1-20x/(1-22x/(1-25x/(1-.../(1-A047215(n+1)*x/(1-... (continued fraction). - Paul Barry, Dec 03 2009
a(n) = (-3)^n*Sum_{k=0..n} (5/3)^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
D-finite with recurrence: a(n) +(-5*n+3)*a(n-1) = 0. - R. J. Mathar, Dec 03 2012
G.f.: 1/G(0) where G(k) = 1 - x*(5*k+2)/( 1 - 5*x*(k+1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 23 2013
Sum_{n>=0} 1/a(n) = 1 + (e/5^3)^(1/5)*(Gamma(2/5) - Gamma(2/5, 1/5)). - Amiram Eldar, Dec 19 2022

A052562 a(n) = 5^n * n!.

Original entry on oeis.org

1, 5, 50, 750, 15000, 375000, 11250000, 393750000, 15750000000, 708750000000, 35437500000000, 1949062500000000, 116943750000000000, 7601343750000000000, 532094062500000000000, 39907054687500000000000

Offset: 0

Joe Keane (jgk(AT)jgk.org)

A simple regular expression in a labeled universe.

For n >= 1 a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_5)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001

Vincenzo Librandi, Table of n, a(n) for n = 0..300
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 504
Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, Vol. 9 (2006), Article 06.1.1.

Cf. A000142, A008548, A008546, A034325, A000165, A047055, A047056, A051150, A092514, A092618.

[5^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 05 2011

spec := [S,{S=Sequence(Union(Z,Z,Z,Z,Z))},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
with(combstruct):A:=[N,{N=Cycle(Union(Z$5))},labeled]: seq(count(A,size=n)/5,n=1..16); # Zerinvary Lajos, Dec 05 2007

Table[5^n*n!, {n, 0, 20}] (* Wesley Ivan Hurt, Sep 28 2013 *)

{a(n) = 5^n*n!}; \\ G. C. Greubel, May 05 2019

[5^n*factorial(n) for n in (0..20)] # G. C. Greubel, May 05 2019

a(n) = A051150(n+1, 0) (first column of triangle).
E.g.f.: 1/(1-5*x).
a(n) = 5*n*a(n-1) with a(0)=1.
G.f.: 1/(1-5*x/(1-5*x/(1-10*x/(1-10*x/(1-15*x/(1-15*x/(1-20*x/(1-... (continued fraction). - Philippe Deléham, Jan 08 2012
G.f.: 1/Q(0), where Q(k) = 1 - 5*x*(2*k+1) - 25*x^2*(k+1)^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Sep 28 2013
a(n) = n!*A000351(n). - R. J. Mathar, Aug 21 2014
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/5) (A092514).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/5) (A092618). (End)

Name changed by Arkadiusz Wesolowski, Oct 04 2011

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)t(n-1,m) + f(n)t(n,m-1), where f(x) = x + 2.

Original entry on oeis.org

1, 2, 2, 4, 12, 4, 8, 52, 52, 8, 16, 196, 416, 196, 16, 32, 684, 2644, 2644, 684, 32, 64, 2276, 14680, 26440, 14680, 2276, 64, 128, 7340, 74652, 220280, 220280, 74652, 7340, 128, 256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172, 256, 512, 72076, 1637860, 10978444, 27227908, 27227908, 10978444, 1637860, 72076, 512

Offset: 0

Dale Gerdemann, Apr 12 2015

Related triangles may be found by varying the function f(x). If f(x) is a linear function, it can be parameterized as f(x) = a*x + b. With different values for a and b, the following triangles are obtained:
  a\b 1.......2.......3.......4.......5.......6
  -1  A144431
  0   A007318 A038208 A038221
  1   A008292 A256890 A257180 A257606 A257607
  2   A060187 A257609 A257611 A257613 A257615
  3   A142458 A257610 A257620 A257622 A257624 A257626
  4   A142459 A257612 A257621
  5   A142460 A257614 A257623
  6   A142461 A257616 A257625
  7   A142462 A257617 A257627
  8   A167884 A257618
  9   A257608 A257619
The row sums of these, and similarly constructed number triangles, are shown in the following table:
  a\b 1.......2.......3.......4.......5.......6.......7.......8.......9
  0   A000079 A000302 A000400
  1   A000142 A001715 A001725 A049388 A049198
  2   A000165 A002866 A002866 A051580 A051582
  3   A008544 A051578 A037559 A051605 A051607 A051609
  4   A001813 A047053 A000407 A034177 A051618 A051620 A051622
  5   A047055 A008546 A008548 A034300 A034325 A051688 A051690
  6   A047657 A049308 A047058 A034689 A034724 A034788 A053101 A053103
  7   A084947 A144827 A049209 A045754 A034830 A034832 A034834 A053105
  8   A084948 A144828 A147626 A051189 A034908 A034910 A034912 A034976 A053115
  9   A084949 A144829 A147630 A049211 A045756 A035013 A035018 A035021 A035023
  10                                  A051262 A035265 A035273         A035277
  11                                  A254322
  12                                          A145448
The formula can be further generalized to: t(n,m) = f(m+s)*t(n-1,m) + f(n-s)*t(n,m-1), where f(x) = a*x + b. The following table specifies triangles with nonzero values for s (given after the slash).
  a\b  0           1           2          3
  -2   A130595/1
  -1
  0
  1    A110555/-1  A120434/-1  A144697/1  A144699/2
With the absolute value, f(x) = |x|, one obtains A038221/3, A038234/4, A038247/5, A038260/6, A038273/7, A038286/8, A038299/9 (with value for s after the slash).
If f(x) = A000045(x) (Fibonacci) and s = 1, the result is A010048 (Fibonomial).
In the notation of Carlitz and Scoville, this is the triangle of generalized Eulerian numbers A(r, s | alpha, beta) with alpha = beta = 2. Also the array A(2,1,4) in the notation of Hwang et al. (see page 31). - Peter Bala, Dec 27 2019

			Array, t(n, k), begins as:
   1,    2,      4,        8,        16,         32,          64, ...;
   2,   12,     52,      196,       684,       2276,        7340, ...;
   4,   52,    416,     2644,     14680,      74652,      357328, ...;
   8,  196,   2644,    26440,    220280,    1623964,    10978444, ...;
  16,  684,  14680,   220280,   2643360,   27227908,   251195000, ...;
  32, 2276,  74652,  1623964,  27227908,  381190712,  4677894984, ...;
  64, 7340, 357328, 10978444, 251195000, 4677894984, 74846319744, ...;
Triangle, T(n, k), begins as:
    1;
    2,     2;
    4,    12,      4;
    8,    52,     52,       8;
   16,   196,    416,     196,      16;
   32,   684,   2644,    2644,     684,      32;
   64,  2276,  14680,   26440,   14680,    2276,     64;
  128,  7340,  74652,  220280,  220280,   74652,   7340,   128;
  256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172,   256;

Michael De Vlieger, Table of n, a(n) for n = 0..11475 (rows 0 <= n <= 150, flattened.)
L. Carlitz and R. Scoville, Generalized Eulerian numbers: combinatorial applications, J. für die reine und angewandte Mathematik, 265 (1974): 110-37. See Section 3.
Dale Gerdemann, A256890, Plot of t(m,n) mod k , YouTube, 2015.
Hsien-Kuei Hwang, Hua-Huai Chern, and Guan-Huei Duh, An asymptotic distribution theory for Eulerian recurrences with applications, arXiv:1807.01412 [math.CO], 2018-2019.

Cf. A000079, A001715, A008292, A038208, A257180, A257606, A257607, A257609.

Cf. A257610, A257612, A257614, A257616, A257617, A257618, A257619.

A256890:= func< n,k | (&+[(-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n: j in [0..k]]) >;
[A256890(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Oct 18 2022

Table[Sum[(-1)^(k-j)*Binomial[j+3, j] Binomial[n+4, k-j] (j+2)^n, {j,0,k}], {n,0, 9}, {k,0,n}]//Flatten (* Michael De Vlieger, Dec 27 2019 *)

t(n,m) = if ((n<0) || (m<0), 0, if ((n==0) && (m==0), 1, (m+2)*t(n-1, m) + (n+2)*t(n, m-1)));
tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n-k, k), ", ");); print(););} \\ Michel Marcus, Apr 14 2015

def A256890(n,k): return sum((-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n for j in range(k+1))
flatten([[A256890(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Oct 18 2022

T(n,k) = t(n-k, k); t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0 else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.
Sum_{k=0..n} T(n, k) = A001715(n).
T(n,k) = Sum_{j = 0..k} (-1)^(k-j)*binomial(j+3,j)*binomial(n+4,k-j)*(j+2)^n. - Peter Bala, Dec 27 2019
Modified rule of Pascal: T(0,0) = 1, T(n,k) = 0 if k < 0 or k > n else T(n,k) = f(n-k) * T(n-1,k-1) + f(k) * T(n-1,k), where f(x) = x + 2. - Georg Fischer, Nov 11 2021
From G. C. Greubel, Oct 18 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = A000079(n). (End)

A027637 a(n) = Product_{i=1..n} (4^i - 1).

Original entry on oeis.org

1, 3, 45, 2835, 722925, 739552275, 3028466566125, 49615367752825875, 3251543125681443718125, 852369269595510700600441875, 893773106866112632882108339078125, 3748755223447856814435325652920396921875

Offset: 0

N. J. A. Sloane

The q-analog of double factorials (A000165) evaluated at q=2. - Michael Somos, Sep 12 2014
3^n*5^(floor(n/2))|a(n) for n>=1. - G. C. Greubel, Nov 21 2015
Given probability p = 1/4^n that an outcome will occur at the n-th stage of an infinite process, then starting at n=1, 1-a(n)/A053763(n+1) is the probability that the outcome has occurred up to and including the n-th iteration. The limiting ratio is 1-A100221 ~ 0.3114625. - Bob Selcoe, Mar 01 2016

G. C. Greubel, Table of n, a(n) for n = 0..50

Cf. A000165.
Sequences of the form q-Pochhammer(n, q, q): A005329 (q=2), A027871 (q=3), this sequence (q=4), A027872 (q=5), A027873 (q=6), A027875 (q=7), A027876 (q=8), A027877 (q=9), A027878 (q=10), A027879 (q=11), A027880 (q=12).
Cf. A053763, A100221.

[1] cat [&*[4^k-1: k in [1..n]]: n in [1..11]]; // Vincenzo Librandi, Dec 24 2015

A027637 := proc(n)
    mul( 4^i-1,i=1..n) ;
end proc:
seq(A027637(n),n=0..8) ;

A027637 = Table[Product[4^i - 1, {i, n}], {n, 0, 9}] (* Alonso del Arte, Nov 14 2011 *)
a[ n_] := If[ n < 0, 0, Product[ (q^(2 k) - 1) / (q - 1), {k, n}] /. q -> 2]; (* Michael Somos, Sep 12 2014 *)
Abs@QPochhammer[4, 4, Range[0, 10]] (* Vladimir Reshetnikov, Nov 20 2015 *)

a(n) = prod(i=1, n, 4^i-1); \\ Michel Marcus, Nov 21 2015

from sage.combinat.q_analogues import q_pochhammer
def A027637(n): return (-1)^n*q_pochhammer(n, 4, 4)
[A027637(n) for n in (0..15)] # G. C. Greubel, Aug 04 2022

a(n) ~ c * 2^(n*(n+1)), where c = Product_{k>=1} (1-1/4^k) = A100221 = 0.688537537120339715456514357293508184675549819378... . - Vaclav Kotesovec, Nov 21 2015
a(n) = 4^(binomial(n+1,2))*(1/4;1/4){n} = (4; 4){n}, where (a;q){n} is the q-Pochhammer symbol. - _G. C. Greubel, Dec 24 2015
G.f.: Sum_{n>=0} 4^(n*(n+1)/2)*x^n / Product_{k=0..n} (1 + 4^k*x). - Ilya Gutkovskiy, May 22 2017
Sum_{n>=0} (-1)^n/a(n) = A100221. - Amiram Eldar, May 07 2023

A084947 a(n) = Product_{i=0..n-1} (7*i+2).

Original entry on oeis.org

1, 2, 18, 288, 6624, 198720, 7352640, 323516160, 16499324160, 956960801280, 62202452083200, 4478576549990400, 353807547449241600, 30427449080634777600, 2829752764499034316800, 282975276449903431680000, 30278354580139667189760000, 3451732422135922059632640000

Offset: 0

Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003

Harvey P. Dale, Table of n, a(n) for n = 0..340

Cf. A000079, A000142, A000165, A008544, A001813, A045754, A047055, A047657, A084942, A084947, A084948, A084949, A144739, A144827, A049209, A051188.

List([0..20], n-> Product([0..n-1], k-> 7*k+2) ); # G. C. Greubel, Aug 18 2019

[ 1 ] cat [ &*[ (7*k+2): k in [0..n-1] ]: n in [1..15] ]; // Klaus Brockhaus, Nov 10 2008

a := n->product(7*i+2,i=0..n-1); [seq(a(j),j=0..30)];

Join[{1},FoldList[Times,7*Range[0,15]+2]] (* Harvey P. Dale, Nov 27 2015 *)
Table[7^n*Pochhammer[2/7, n], {n,0,15}] (* G. C. Greubel, Aug 18 2019 *)

vector(20, n, n--; prod(k=0, n-1, 7*k+2)) \\ G. C. Greubel, Aug 18 2019

[product(7*k+2 for k in (0..n-1)) for n in (0..20)] # G. C. Greubel, Aug 18 2019

a(n) = A084942(n)/A000142(n)*A000079(n) = 7^n*Pochhammer(2/7, n) = 7^n*Gamma(n+2/7)/Gamma(2/7).
D-finite with recurrence a(0) = 1; a(n) = (7*n - 5)*a(n-1) for n > 0. - Klaus Brockhaus, Nov 10 2008
G.f.: 1/(1-2*x/(1-7*x/(1-9*x/(1-14*x/(1-16*x/(1-21*x/(1-23*x/(1-28*x/(1-... (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-5)^n*Sum_{k=0..n} (7/5)^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
From Ilya Gutkovskiy, Mar 23 2017: (Start)
E.g.f.: 1/(1 - 7*x)^(2/7).
a(n) ~ sqrt(2*Pi)*7^n*n^n/(exp(n)*n^(3/14)*Gamma(2/7)). (End)
Sum_{n>=0} 1/a(n) = 1 + (e/7^5)^(1/7)*(Gamma(2/7) - Gamma(2/7, 1/7)). - Amiram Eldar, Dec 19 2022

a(15) from Klaus Brockhaus, Nov 10 2008

A047058 a(n) = 6^n * n!.

Original entry on oeis.org

1, 6, 72, 1296, 31104, 933120, 33592320, 1410877440, 67722117120, 3656994324480, 219419659468800, 14481697524940800, 1042682221795737600, 81329213300067532800, 6831653917205672755200, 614848852548510547968000

Offset: 0

Joe Keane (jgk(AT)jgk.org)

For n >= 1, a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_6)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001
a(n) is the number of ways 3 members of each of n different teams can be arranged in a row so that members of the same team are together. - Geoffrey Critzer, Mar 30 2009
From Jianing Song, Mar 29 2021: (Start)
Number of n X n monomial matrices with entries 0, +/-1, +/-w, +/-w^2, where w = (-1 + sqrt(3)*i)/2 is a primitive 3rd root of unity.
a(n) is the order of the group U_n(Z[w]) = {A in M_n(Z[w]): A*A^H = I_n}, the group of n X n unitary matrices over the Eisenstein integers. Here A^H is the conjugate transpose of A. (End)

Vincenzo Librandi, Table of n, a(n) for n = 0..300
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 528.

Cf. A000142, A000165, A008542, A008543, A047053, A051151, A053103, A092515, A092727.

[6^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 05 2011

seq( 6^n*n!, n=0..20); # G. C. Greubel, Jun 08 2020

Table[6^n n!,{n,0,20}] (* Harvey P. Dale, Mar 30 2018 *)

a(n) = A051151(n+1, 0).
E.g.f.: 1/(1 - 6*x).
G.f.: 1/(1 -6*x/(1 - 6*x/(1 - 12*x/(1 - 12*x/(1 - 18*x/(1 - 18*x/(1 - 24*x/(1 - 24*x/(1 - 30*x/(1 - 30*x/(1 -... (continued fraction). - Philippe Deléham, Jan 08 2012
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/6) (A092515).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/6) (A092727). (End)

Name changed by Arkadiusz Wesolowski, Oct 04 2011

A079484 a(n) = (2n-1)!! * (2n+1)!!, where the double factorial is A006882.

Original entry on oeis.org

1, 3, 45, 1575, 99225, 9823275, 1404728325, 273922023375, 69850115960625, 22561587455281875, 9002073394657468125, 4348001449619557104375, 2500100833531245335015625, 1687568062633590601135546875, 1321365793042101440689133203125

Offset: 0

Benoit Cloitre, Jan 17 2003

nonn

a(n) is the determinant of M(2n+1) where M(k) is the k X k matrix with m(i,j)=j if i+j=k m(i,j)=i otherwise. - Adapted to offset 0, Rainer Rosenthal, Jun 19 2024
In the following two comments on the calculation of the terms using permanents, offset 1 is assumed. In the corresponding PARI code, this is implemented with offset 0. - Hugo Pfoertner, Jun 23 2024
(-1)^n*a(n)/2^(2n-1) is the permanent of the (m X m) matrix {1/(x_i-y_j), 1<=i<=m, 1<=j<=m}, where x_1,x_2,...,x_m are the zeros of x^m-1 and y_1,y_2,...,y_m the zeros of y^m+1 and m=2n-1.
In 1881, R. F. Scott posed a conjecture that the absolute value of permanent of square matrix with elements a(i,j)= (x_i - y_j)^(-1), where x_1,...,x_n are roots of x^n=1, while y_1,...,y_n are roots of y^n=-1, equals a((n-1)/2)/2^n, if n>=1 is odd, and 0, if n>=2 is even. After a century (in 1979), the conjecture was proved by H. Minc. - Vladimir Shevelev, Dec 01 2013
a(n) is the number of permutations in S_{2n+1} in which all cycles have odd length. - José H. Nieto S., Jan 09 2012
Number of 3-bundled increasing bilabeled trees with 2n labels. - Markus Kuba, Nov 18 2014
a(n) is the number of rooted, binary, leaf-labeled topologies with 2n+2 leaves that have n+1 cherry nodes. - Noah A Rosenberg, Feb 12 2019

			G.f. = 1 + 3*x + 45*x^2 + 1575*x^3 + 99225*x^4 + 9823275*x^5 + ...
M(5) =
  [1, 2, 3, 1, 5]
  [1, 2, 2, 4, 5]
  [1, 3, 3, 4, 5]
  [4, 2, 3, 4, 5]
  [1, 2, 3, 4, 5].
Integral_{x=0..oo} x^3*BesselK(1, sqrt(x)) = 1575*Pi. - _Olivier Gérard_, May 20 2009

Miklós Bóna, A walk through combinatorics, World Scientific, 2006.

Alois P. Heinz, Table of n, a(n) for n = 0..224
Cyril Banderier, Markus Kuba, and Michael Wallner, Analytic Combinatorics of Composition schemes and phase transitions with mixed Poisson distributions, arXiv:2103.03751 [math.PR], 2021.
Guo-Niu Han and Christian Krattenthaler, Rectangular Scott-type permanents, arXiv:math/0003072 [math.RA], 2000.
Markus Kuba and Alois Panholzer, Combinatorial families of multilabelled increasing trees and hook-length formulas, arXiv:1411.4587 [math.CO], 17 Nov 2014.
MathOverflow, Geometric / physical / probabilistic interpretations of Riemann zeta(n>1)?, answer by Tom Copeland posted in Aug 2021.
Henryk Minc, On a conjecture of R. F. Scott (1881), Linear Algebra Appl., Vol. 28 (1979), pp. 141-153.
Theodoros Theodoulidis, On the Closed-Form Expression of Carson’s Integral, Period. Polytech. Elec. Eng. Comp. Sci., Vol. 59, No. 1 (2015), pp. 26-29.
Eric Weisstein's World of Mathematics, Struve function.

Cf. A001818, A000165.
Bisection of A000246, A053195, |A013069|, |A046126|. Cf. A000909.
Cf. A001044, A010791, |A129464|, A114779, are also values of similar moments.
Equals the row sums of A162005.
Cf. A316728.
Diagonal elements of A306364 in even-numbered rows.

I:=[1, 3]; [n le 2 select I[n] else (4*n^2-8*n+3)*Self(n-1): n in [1..20]]; // Vincenzo Librandi, Nov 18 2014

a:= n-> (d-> d(2*n-1)*d(2*n+1))(doublefactorial):
seq(a(n), n=0..15);  # Alois P. Heinz, Jan 30 2013
# second Maple program:
A079484 := n-> LinearAlgebra[Determinant](Matrix(2*n+1, (i, j)-> `if`(i+j=2*n+1, j, i))): seq(A079484(n), n=0..14); # Rainer Rosenthal, Jun 18 2024

a[n_] := (2n - 1)!!*(2n + 1)!!; Table[a[n], {n, 0, 13}] (* Jean-François Alcover, Jan 30 2013 *)

/* Formula using the zeta function and a log integral:*/
L(n)= intnum(t=0, 1, log(1-1/t)^n);
Zetai(n)= -I*I^n*(2*Pi)^(n-1)/(n-1)*L(1-n);
a(m)={my(n=m+1);round(real(-I*2^(2*n-1)*Zetai(1/2-n)*L(-1/2+n)/(Zetai(-1/2+n)*L(1/2-n))))};
/* Gerry Martens, Mar 07 2011, adapted to offset 0 by Hugo Pfoertner, Jun 19 2024 */

{a(n) = if( n<0, -1 / self()(-1-n), (2*n + 1)! * (2*n)! / (n! * 2^n)^2 )}; /* Michael Somos, May 04 2017 */

{a(n) = if( n<0, -1 / self()(-1-n), my(m = 2*n + 1); m! * polcoeff( x / sqrt( 1 - x^2 + x * O(x^m) ), m))}; /* Michael Somos, May 04 2017 */

\\ using the Pochhammer symbol
a(n) = {my(P(x,k)=gamma(x+k)/gamma(x)); 4^n*round(P(1/2,n)*P(3/2,n))} \\ Hugo Pfoertner, Jun 20 2024

\\ Scott's (1881) method
a(n) = {my(m=2*n+1, X = polroots(x^m-1), Y = polroots(x^m+1), M = matrix(m, m, i, j, 1/(X[i]-Y[j]))); (-1)^n * round(2^m * real(matpermanent(M)))}; \\ Hugo Pfoertner, Jun 23 2024

D-finite with recurrence a(n) = (4*n^2 - 1) * a(n-1) for all n in Z.
a(n) = A001147(n)*A001147(n+1).
E.g.f.: 1/(1-x^2)^(3/2) (with interpolated zeros). - Paul Barry, May 26 2003
a(n) = (2n+1)! * C(2n, n) / 2^(2n). - Ralf Stephan, Mar 22 2004.
Alternatingly signed values have e.g.f. sqrt(1+x^2).
a(n) is the value of the n-th moment of (1/Pi)*BesselK(1, sqrt(x)) on the positive part of the real line. - Olivier Gérard, May 20 2009
a(n) = -2^(2*n-1)*exp(i*n*Pi)*gamma(1/2+n)/gamma(3/2-n). - Gerry Martens, Mar 07 2011
E.g.f. (odd powers) tan(arcsin(x)) = Sum_{n>=0} (2n-1)!!*(2n+1)!!*x^(2*n+1)/(2*n+1)!. - Vladimir Kruchinin, Apr 22 2011
G.f.: 1 + x*(G(0) - 1)/(x-1) where G(k) = 1 - ((2*k+2)^2-1)/(1-x/(x - 1/G(k+1))); ( continued fraction ). - Sergei N. Gladkovskii, Jan 15 2013
a(n) = (2^(2*n+3)*Gamma(n+3/2)*Gamma(n+5/2))/Pi. - Jean-François Alcover, Jul 20 2015
Limit_{n->oo} 4^n*(n!)^2/a(n) = Pi/2. - Daniel Suteu, Feb 05 2017
From Michael Somos, May 04 2017: (Start)
a(n) = (2*n + 1) * A001818(n).
E.g.f.: Sum_{n>=0} a(n) * x^(2*n+1) / (2*n+1)! = x / sqrt(1 - x^2) = tan(arcsin(x)).
Given e.g.f. A(x) = y, then x * y' = y + y^3.
a(n) = -1 / a(-1-n) for all n in Z.
0 = +a(n)*(+288*a(n+2) -60*a(n+3) +a(n+4)) +a(n+1)*(-36*a(n+2) -4*a(n+3)) +a(n+2)*(+3*a(n+2)) for all n in Z. (End)
a(n) = Sum_{k=0..2n} (k+1) * A316728(n,k). - Alois P. Heinz, Jul 12 2018
From Amiram Eldar, Mar 18 2022: (Start)
Sum_{n>=0} 1/a(n) = 1 + L_1(1)*Pi/2, where L is the modified Struve function.
Sum_{n>=0} (-1)^n/a(n) = 1 - H_1(1)*Pi/2, where H is the Struve function. (End)

Simpler description from Daniel Flath (deflath(AT)yahoo.com), Mar 05 2004

A039683 Signed double Pochhammer triangle: expansion of x(x-2)(x-4)..(x-2n+2).

Original entry on oeis.org

1, -2, 1, 8, -6, 1, -48, 44, -12, 1, 384, -400, 140, -20, 1, -3840, 4384, -1800, 340, -30, 1, 46080, -56448, 25984, -5880, 700, -42, 1, -645120, 836352, -420224, 108304, -15680, 1288, -56, 1, 10321920, -14026752, 7559936, -2153088, 359184, -36288, 2184, -72, 1

Offset: 1

Wouter Meeussen

T(n,m) = R_n^m(a=0,b=2) in the notation of the given reference.
Exponential Riordan array [1/(1+2x),log(1+2x)/2]. The unsigned triangle is [1/(1-2x),log(1/sqrt(1-2x))]. - Paul Barry_, Apr 29 2009
The n-th row is related to the expansion of z^(-2n)*(z^3 d/dz)^n in polynomials of the Euler operator D=(z d/dz). E.g., z^(-6)(z^3 d/dz)^3 = D^3 + 6 D^2 + 8 D. See Copeland link for relations to Bell / Exponential / Touchard polynomial operators. - Tom Copeland, Nov 14 2013
A refinement of this array is given by A231846. - Tom Copeland, Nov 15 2013
Also the Bell transform of the double factorial of even numbers A000165 except that the values are unsigned and in addition a first column (1,0,0 ...) is added on the left side of the triangle. For the Bell transform of the double factorial of odd numbers A001147 see A132062. For the definition of the Bell transform see A264428. - Peter Luschny, Dec 20 2015
The signed triangle is also the inverse Bell transform of A000079 (see Luschny link). - John Keith, Nov 24 2020

			Triangle starts:
  {1},
  {2,1},
  {8,6,1},
  {48,44,12,1},
  ...
From _Paul Barry_, Apr 29 2009: (Start)
The unsigned triangle [1/(1-2x),log(1/sqrt(1-2x))] has production matrix:
  2, 1,
  4, 4, 1,
  8, 12, 6, 1,
  16, 32, 24, 8, 1,
  32, 80, 80, 40, 10, 1,
  64, 192, 240, 160, 60, 12, 1
which is A007318^{2} beheaded. (End)

Richell O. Celeste, Roberto B. Corcino, and Ken Joffaniel M. Gonzales. Two Approaches to Normal Order Coefficients, Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
Tom Copeland, Addendum to Mathemagical Forests.
P. Feijão, F. V. Martinez, and A. Thévenin, On the distribution of cycles and paths in multichromosomal breakpoint graphs and the expected value of rearrangement distance, BMC Bioinformatics 16:Suppl19 (2015), S1. doi:10.1186/1471-2105-16-S19-S1
Lisa Glaser, Causal set actions in various dimensions, J. Phys.: Conf. Ser. 306 (2011), 012041.
Wolfdieter Lang, First 9 rows and comment.
Peter Luschny, The Bell transform
D. S. Mitrinovic and M. S. Mitrinovic, Tableaux d'une classe de nombres relies aux nombres de Stirling, Univ. Beograd. Pubi. Elektrotehn. Fak. Ser. Mat. Fiz. 77 (1962).

First column (unsigned triangle) is (2(n-1))!! = 1, 2, 8, 48, 384...= A000165(n-1) and the row sums (unsigned) are (2n-1)!! = 1, 3, 15, 105, 945... = A001147(n-1).
Cf. A051141, A051142.
Cf. A000165, A132062, A264428.
Cf. A038207.

Table[ Rest@ CoefficientList[ Product[ z-k, {k, 0, 2p-2, 2} ], z ], {p, 6} ]

# uses[bell_transform from A264428]
# Unsigned values and an additional first column (1,0,0,...).
def A039683_unsigned_row(n):
    a = sloane.A000165
    dblfact = a.list(n)
    return bell_transform(n, dblfact)
[A039683_unsigned_row(n) for n in (0..9)] # Peter Luschny, Dec 20 2015

T(n, m) = T(n-1, m-1) - 2*(n-1)*T(n-1, m), n >= m >= 1; T(n, m) := 0, n

E.g.f. for m-th column of signed triangle: (((log(1+2*x))/2)^m)/m!.

E.g.f.: (1+2*x)^(y/2). O.g.f. for n-th row of signed triangle: Sum_{m=0..n} Stirling1(n, m)*2^(n-m)*x^m. - Vladeta Jovovic, Feb 11 2003

T(n, m) = S1(n, m)*2^(n-m), with S1(n, m) := A008275(n, m) (signed Stirling1 triangle).

The production matrix below is A038207 with the first row removed. With the initial index n = 0, the associated differential raising operator is R = e^(2D)*x = (2+x)*e^(2D) with D = d/dx, i.e., R p_n(x) = p_(n+1)(x) where p_n(x) is the n-th unsigned row polynomial and p_0(x) = 1, so p_(n+1)(x) = (2+x) * p_n(2+x). - Tom Copeland, Oct 11 2016

Additional comments from Wolfdieter Lang

Title revised by Tom Copeland, Dec 21 2013

A047657 Sextuple factorial numbers: a(n) = Product_{k=0..n-1} (6*k+2).

Original entry on oeis.org

1, 2, 16, 224, 4480, 116480, 3727360, 141639680, 6232145920, 311607296000, 17450008576000, 1081900531712000, 73569236156416000, 5444123475574784000, 435529878045982720000, 37455569511954513920000, 3445912395099815280640000, 337699414719781897502720000

Offset: 0

Joe Keane (jgk(AT)jgk.org)

G. C. Greubel, Table of n, a(n) for n = 0..340

Cf. A007559, A008542, A011781, A073005.

Cf. A000165, A008544, A001813, A047055, A084947, A084948, A084949.

List([0..20], n-> Product([0..n-1], k-> 6*k+2) ); # G. C. Greubel, Aug 18 2019

[1] cat [(&*[6*k+2: k in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Aug 18 2019

a:= n->product(6*j+2, j=0..n-1); seq(a(n), n=0..20); # G. C. Greubel, Aug 18 2019

b[1]=2; b[n_]:= b[n] = b[n-1] +6; a[0]=1; a[1]=2; a[n_]:= a[n] = a[n-1]*b[n]; Table[a[n], {n,0,20}] (* Roger L. Bagula, Sep 17 2008 *)
FoldList[Times,1,6*Range[0,20]+2] (* Harvey P. Dale, Aug 06 2013 *)
Table[6^n*Pochhammer[1/3, n], {n,0,20}] (* G. C. Greubel, Aug 18 2019 *)

vector(20, n, n--; prod(k=0, n-1, 6*k+2)) \\ G. C. Greubel, Aug 18 2019

[product(6*k+2 for k in (0..n-1)) for n in (0..20)] # G. C. Greubel, Aug 18 2019

E.g.f.: (1-6*x)^(-1/3).
a(n) = 2^n*A007559(n).
a(n) = A084941(n)/A000142(n)*A000079(n) = 6^n*Pochhammer(1/3, n) = 1/2*6^n*Gamma(n+1/3)*sqrt(3)*Gamma(2/3)/Pi. - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
Let b(n) = b(n-1) + 6; then a(n) = b(n)*a(n-1). - Roger L. Bagula, Sep 17 2008
G.f.: 1/(1-2*x/(1-6*x/(1-8*x/(1-12*x/(1-14*x/(1-18*x/(1-20*x/(1-24*x/(1-26*x/(1-... (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-4)^n*Sum_{k=0..n} (3/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/G(0) where G(k) = 1 - x*(6*k+2)/( 1 - 6*x*(k+1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 23 2013
D-finite with recurrence: a(n) +2*(-3*n+2)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
Sum_{n>=0} 1/a(n) = 1 + exp(1/6)*(Gamma(1/3) - Gamma(1/3, 1/6))/6^(2/3). - Amiram Eldar, Dec 18 2022
a(n) ~ sqrt(Pi) * 2^(n+1/2) * (3/e)^n * n^(n-1/6) / Gamma(1/3). - Amiram Eldar, Sep 01 2025

A051577 a(n) = (2*n + 3)!!/3 = A001147(n+2)/3.

Original entry on oeis.org

1, 5, 35, 315, 3465, 45045, 675675, 11486475, 218243025, 4583103525, 105411381075, 2635284526875, 71152682225625, 2063427784543125, 63966261320836875, 2110886623587616875, 73881031825566590625, 2733598177545963853125, 106610328924292590271875

Offset: 0

Wolfdieter Lang

Row m = 3 of the array A(3; m,n) := (2*n+m)!!/m!!, m >= 0, n >= 0.

G. C. Greubel, Table of n, a(n) for n = 0..250
H. W. Gould, Harris Kwong, and Jocelyn Quaintance, On Certain Sums of Stirling Numbers with Binomial Coefficients, J. Integer Sequences, 18 (2015), Article 15.9.6.
Guo-Niu Han, Enumeration of Standard Puzzles, 2011. [Cached copy]
Guo-Niu Han, Enumeration of Standard Puzzles, arXiv:2006.14070 [math.CO], 2020.
A. N. Stokes, Continued fraction solutions of the Riccati equation, Bull. Austral. Math. Soc. Vol. 25 (1982), 207-214.
Michael Wallner, A bijection of plane increasing trees with relaxed binary trees of right height at most one, arXiv:1706.07163 [math.CO], 2017, p. 10.
Eric Weisstein's World of Mathematics, Double Factorial.
Wikipedia, Double factorial.

Cf. A000165, A001147, A002866(n+1) (m=0..2 rows of A(3; m,n)).

Cf. A051578, A051579, A051580, A051581, A051582, A051583, A178647.

F:=Factorial;; List([0..25], n-> F(2*n+4)/(12*2^n*F(n+2)) ); # G. C. Greubel, Nov 12 2019

F:=Factorial; [F(2*n+4)/(12*2^n*F(n+2)): n in [0..25]]; // G. C. Greubel, Nov 12 2019

seq( doublefactorial(2*n+3)/3,n=0..10) ; # R. J. Mathar, Sep 29 2013

Table[(2*n + 3)!!/3!!, {n, 0, 25}] (* G. C. Greubel, Jan 22 2017 *)
a[n_] := Sum[(-1)^k*Binomial[2*n + 1, n + k]*StirlingS1[n + k + 1 ,k], {k , 1, n + 1}]; Flatten[Table[a[n], {n, 0, 18}]] (* Detlef Meya, Jan 17 2024 *)

vector(26, n, (2*n+2)!/(6*2^n*(n+1)!) ) \\ G. C. Greubel, Nov 12 2019

f=factorial; [f(2*n+4)/(12*2^n*f(n+2)) for n in (0..25)] # G. C. Greubel, Nov 12 2019

a(n) = (2*n + 3)!!/3!!.
E.g.f.: 1/(1 - 2*x)^(5/2).
a(n) ~ (4/3) * sqrt(2) * n^2 * 2^n * e^(-n) * n^n *{1 + (47/24)*n^(-1) + ...}. - Joe Keane (jgk(AT)jgk.org), Nov 23 2001
Ramanujan polynomials -psi_n(n, x) evaluated at 0. - Ralf Stephan, Apr 16 2004
a(n) = 2^(2 + n) * Gamma(n + 5/2)/(3 * sqrt(Pi)). - Gerson Washiski Barbosa, May 05 2010
From Peter Bala, May 26 2017: (Start)
D-finite with recurrence a(n+1) = (2*n + 5)*a(n) with a(0) = 1.
O.g.f. A(x) satisfies the Riccati differential equation 2*x^2*A(x)' = (1 - 5*x)*A(x) - 1 with A(0) = 1.
G.f. as an S-fraction: A(x) = 1/(1 - 5*x/(1 - 2*x/(1 - 7*x/(1 - 4*x/(1 - 9*x/(1 - 6*x/(1 - ... - (2*n+3)*x/(1 - 2*n*x/(1 - ...))))))))) (by Stokes, 1982).
Reciprocal as an S-fraction: 1/A(x) = 1/(1 + 5*x/(1 - 7*x/(1 - 2*x/(1 - 9*x/(1 - 4*x/(1 - 11*x/(1 - 6*x/(1 - ... - (2*n + 5)*x/(1 - 2*n*x/(1 - ...)))))))))). (End)
From Amiram Eldar, Dec 11 2022: (Start)
Sum_{n>=0} 1/a(n) = 3*(sqrt(e*Pi/2) * erf(1/sqrt(2)) - 1), where erf is the error function.
Sum_{n>=0} (-1)^n/a(n) = 3*(1 - sqrt(Pi/(2*e)) * erfi(1/sqrt(2))), where erfi is the imaginary error function. (End)
a(n) = Sum_{k=1..n+1} (-1)^k*binomial(2*n + 1, n + k)*Stirling1(n + k + 1, k). - Detlef Meya, Jan 17 2024

Previous Showing 41-50 of 231 results. Next

cp's OEIS Frontend

A047055 Quintuple factorial numbers: a(n) = Product_{k=0..n-1} (5*k + 2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A052562 a(n) = 5^n * n!.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

SageMath

Formula

A027637 a(n) = Product_{i=1..n} (4^i - 1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

A084947 a(n) = Product_{i=0..n-1} (7*i+2).

Views

Author

Keywords

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A047058 a(n) = 6^n * n!.

Views

Author

Keywords

Comments

Links

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)t(n-1,m) + f(n)t(n,m-1), where f(x) = x + 2.