A047058 -id:A047058 - cp's OEIS frontend

A000165 Double factorial of even numbers: (2n)!! = 2^n*n!.

1, 2, 8, 48, 384, 3840, 46080, 645120, 10321920, 185794560, 3715891200, 81749606400, 1961990553600, 51011754393600, 1428329123020800, 42849873690624000, 1371195958099968000, 46620662575398912000, 1678343852714360832000, 63777066403145711616000

Offset: 0

N. J. A. Sloane

a(n) is also the size of the automorphism group of the graph (edge graph) of the n-dimensional hypercube and also of the geometric automorphism group of the hypercube (the two groups are isomorphic). This group is an extension of an elementary Abelian group (C_2)^n by S_n. (C_2 is the cyclic group with two elements and S_n is the symmetric group.) - Avi Peretz (njk(AT)netvision.net.il), Feb 21 2001
Then a(n) appears in the power series: sqrt(1+sin(y)) = Sum_{n>=0} (-1)^floor(n/2)*y^(n)/a(n) and sqrt((1+cos(y))/2) = Sum_{n>=0} (-1)^n*y^(2n)/a(2n). - Benoit Cloitre, Feb 02 2002
Appears to be the BinomialMean transform of A001907. See A075271. - John W. Layman, Sep 28 2002
Number of n X n monomial matrices with entries 0, +-1.
Also number of linear signed orders.
Define a "downgrade" to be the permutation d which places the items of a permutation p in descending order. This note concerns those permutations that are equal to their double-downgrades. The number of permutations of order 2n having this property are equinumerous with those of order 2n+1. a(n) = number of double-downgrading permutations of order 2n and 2n+1. - Eugene McDonnell (eemcd(AT)mac.com), Oct 27 2003
a(n) = (Integral_{x=0..Pi/2} cos(x)^(2*n+1) dx) where the denominators are b(n) = (2*n)!/(n!*2^n). - Al Hakanson (hawkuu(AT)excite.com), Mar 02 2004
1 + (1/2)x - (1/8)x^2 - (1/48)x^3 + (1/384)x^4 + ... = sqrt(1+sin(x)).
a(n)*(-1)^n = coefficient of the leading term of the (n+1)-th derivative of arctan(x), see Hildebrand link. - Reinhard Zumkeller, Jan 14 2006
a(n) is the Pfaffian of the skew-symmetric 2n X 2n matrix whose (i,j) entry is j for iDavid Callan, Sep 25 2006
a(n) is the number of increasing plane trees with n+1 edges. (In a plane tree, each subtree of the root is an ordered tree but the subtrees of the root may be cyclically rotated.) Increasing means the vertices are labeled 0,1,2,...,n+1 and each child has a greater label than its parent. Cf. A001147 for increasing ordered trees, A000142 for increasing unordered trees and A000111 for increasing 0-1-2 trees. - David Callan, Dec 22 2006
Hamed Hatami and Pooya Hatami prove that this is an upper bound on the cardinality of any minimal dominating set in C_{2n+1}^n, the Cartesian product of n copies of the cycle of size 2n+1, where 2n+1 is a prime. - Jonathan Vos Post, Jan 03 2007
This sequence and (1,-2,0,0,0,0,...) form a reciprocal pair under the list partition transform and associated operations described in A133314. - Tom Copeland, Oct 29 2007
a(n) = number of permutations of the multiset {1,1,2,2,...,n,n,n+1,n+1} such that between the two occurrences of i, there is exactly one entry >i, for i=1,2,...,n. Example: a(2) = 8 counts 121323, 131232, 213123, 231213, 232131, 312132, 321312, 323121. Proof: There is always exactly one entry between the two 1s (when n>=1). Given a permutation p in A(n) (counted by a(n)), record the position i of the first 1, then delete both 1s and subtract 1 from every entry to get a permutation q in A(n-1). The mapping p -> (i,q) is a bijection from A(n) to the Cartesian product [1,2n] X A(n-1). - David Callan, Nov 29 2007
Row sums of A028338. - Paul Barry, Feb 07 2009
a(n) is the number of ways to seat n married couples in a row so that everyone is next to their spouse. Compare A007060. - Geoffrey Critzer, Mar 29 2009
From Gary W. Adamson, Apr 21 2009: (Start)
Equals (-1)^n * (1, 1, 2, 8, 48, ...) dot (1, -3, 5, -7, 9, ...).
Example: a(4) = 384 = (1, 1, 2, 8, 48) dot (1, -3, 5, -7, 9) = (1, -3, 10, -56, 432). (End)
exp(x/2) = Sum_{n>=0} x^n/a(n). - Jaume Oliver Lafont, Sep 07 2009
Assuming n starts at 0, a(n) appears to be the number of Gray codes on n bits. It certainly is the number of Gray codes on n bits isomorphic to the canonical one. Proof: There are 2^n different starting positions for each code. Also, each code has a particular pattern of bit positions that are flipped (for instance, 1 2 1 3 1 2 1 for n=3), and these bit position patterns can be permuted in n! ways. - D. J. Schreffler (ds1404(AT)txstate.edu), Jul 18 2010
E.g.f. of 0,1,2,8,... is x/(1-2x/(2-2x/(3-8x/(4-8x/(5-18x/(6-18x/(7-... (continued fraction). - Paul Barry, Jan 17 2011
Number of increasing 2-colored trees with choice of two colors for each edge. In general, if we replace 2 with k we get the number of increasing k-colored trees. For example, for k=3 we get the triple factorial numbers. - Wenjin Woan, May 31 2011
a(n) = row sums of triangle A193229. - Gary W. Adamson, Jul 18 2011
Also the number of permutations of 2n (or of 2n+1) that are equal to their reverse-complements. (See the Egge reference.) Note that the double-downgrade described in the preceding comment (McDonnell) is equivalent to the reverse-complement. - Justin M. Troyka, Aug 11 2011
The e.g.f. can be used to form a generator, [1/(1-2x)] d/dx, for A000108, so a(n) can be applied to A145271 to generate the Catalan numbers. - Tom Copeland, Oct 01 2011
The e.g.f. of 1/a(n) is BesselI(0,sqrt(2*x)). See Abramowitz-Stegun (reference and link under A008277), p. 375, 9.6.10. - Wolfdieter Lang, Jan 09 2012
a(n) = order of the largest imprimitive group of degree 2n with n systems of imprimitivity (see [Miller], p. 203). - L. Edson Jeffery, Feb 05 2012
Row sums of triangle A208057. - Gary W. Adamson, Feb 22 2012
a(n) is the number of ways to designate a subset of elements in each n-permutation. a(n) = A000142(n) + A001563(n) + A001804(n) + A001805(n) + A001806(n) + A001807(n) + A035038(n) * n!. - Geoffrey Critzer, Nov 08 2012
For n>1, a(n) is the order of the Coxeter groups (also called Weyl groups) of types B_n and C_n. - Tom Edgar, Nov 05 2013
For m>0, k*a(m-1) is the m-th cumulant of the chi-squared probability distribution for k degrees of freedom. - Stanislav Sykora, Jun 27 2014
a(n) with 0 prepended is the binomial transform of A120765. - Vladimir Reshetnikov, Oct 28 2015
Exponential self-convolution of A001147. - Vladimir Reshetnikov, Oct 08 2016
Also the order of the automorphism group of the n-ladder rung graph. - Eric W. Weisstein, Jul 22 2017
a(n) is the order of the group O_n(Z) = {A in M_n(Z): A*A^T = I_n}, the group of n X n orthogonal matrices over the integers. - Jianing Song, Mar 29 2021
a(n) is the number of ways to tile a (3n,3n)-benzel or a (3n+1,3n+2)-benzel using left stones and two kinds of bones; see Defant et al., below. - James Propp, Jul 22 2023
a(n) is the number of labeled histories for a labeled topology with the modified lodgepole shape and n+1 cherry nodes. - Noah A Rosenberg, Jan 16 2025

			The following permutations and their reversals are all of the permutations of order 5 having the double-downgrade property:
  0 1 2 3 4
  0 3 2 1 4
  1 0 2 4 3
  1 4 2 0 3
G.f. = 1 + 2*x + 8*x^2 + 48*x^3 + 384*x^4 + 3840*x^5 + 46080*x^6 + 645120*x^7 + ...

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
Paul Barry, General Eulerian Polynomials as Moments Using Exponential Riordan Arrays, Journal of Integer Sequences, 16 (2013), #13.9.6.
Paul Barry, Generalized Eulerian Triangles and Some Special Production Matrices, arXiv:1803.10297 [math.CO], 2018.
Isabel Cação, Helmuth R. Malonek, Maria Irene Falcão, and Graça Tomaz, Combinatorial Identities Associated with a Multidimensional Polynomial Sequence, J. Int. Seq., Vol. 21 (2018), Article 18.7.4.
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
CombOS - Combinatorial Object Server, Generate colored permutations
R. Coquereaux and J.-B. Zuber, Maps, immersions and permutations, arXiv preprint arXiv:1507.03163 [math.CO], 2015.
Colin Defant, Rupert Li, James Propp, and Benjamin Young, Tilings of Benzels via the Abacus Bijection, arXiv preprint, arXiv:2209.05717 [math.CO], 2022.
Eric S. Egge, Restricted symmetric permutations, Ann. Combin., 11 (2007), 405-434.
Peter C. Fishburn, Signed Orders, Choice Probabilities and Linear Polytopes, Journal of Mathematical Psychology, Volume 45, Issue 1, (2001), pp. 53-80.
Joël Gay and Vincent Pilaud, The weak order on Weyl posets, arXiv:1804.06572 [math.CO], 2018.
G. Gordon, The answer is 2^n*n! What is the question?, Amer. Math. Monthly, 106 (1999), 636-645.
Guo-Niu Han, Enumeration of Standard Puzzles
Guo-Niu Han, Enumeration of Standard Puzzles [Cached copy]
Hamed Hatami and Pooya Hatami, Perfect dominating sets in the Cartesian products of prime cycles, arXiv:math/0701018 [math.CO], 2006-2009.
Jason D. Hildebrand, Differentiating Arctan(x)
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 136
E. Lappo and N. A. Rosenberg, A lattice structure for ancestral configurations arising from the relationship between gene trees and species trees, Adv. Appl. Math. 343 (2024), 65-81.
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181. [Annotated scan of pages 180 and 181 only]
E. Lucas, Théorie des nombres (annotated scans of a few selected pages)
Eugene McDonnell, Magic Squares and Permutations, APL Quote Quad 7.3 (Fall 1976).
B. E. Meserve, Double Factorials, American Mathematical Monthly, 55 (1948), 425-426.
G. A. Miller, Groups formed by special matrices, Bull. Amer. Math. Soc. 24 (1918), 203-206.
R. Ondrejka, Tables of double factorials, Math. Comp., 24 (1970), 231.
Alexsandar Petojevic, The Function vM_m(s; a; z) and Some Well-Known Sequences, Journal of Integer Sequences, Vol. 5 (2002), Article 02.1.7.
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014.
R. W. Robinson, Counting arrangements of bishops, pp. 198-214 of Combinatorial Mathematics IV (Adelaide 1975), Lect. Notes Math., 560 (1976).
M. Z. Spivey and L. L. Steil, The k-Binomial Transforms and the Hankel Transform, J. Integ. Seqs. Vol. 9 (2006), #06.1.1.
Eric Weisstein's World of Mathematics, Double Factorial
Eric Weisstein's World of Mathematics, Graph Automorphism
Eric Weisstein's World of Mathematics, Ladder Rung Graph
Index to divisibility sequences
Index entries for sequences related to factorial numbers

Cf. A000142 (n!), A001147 ((2n-1)!!), A032184 (2^n*(n-1)!).
Cf. A000079, A000108, A000111, A001044, A001563, A001804, A001805, A001806.
Cf. A001807, A001813, A001907, A002454, A006882, A007060, A008277, A008544.
Cf. A010050, A019774, A028326, A028338, A035038, A039683, A047053, A047055.
Cf. A047058, A047657, A075271, A084947, A084948, A084949, A092605, A133314.
Cf. A145271, A193229, A208057.
This sequence gives the row sums in A060187, and (-1)^n*a(n) the alternating row sums in A039757.
Also row sums in A028338.
Column k=2 of A329070.

a000165 n = product [2, 4 .. 2 * n]  -- Reinhard Zumkeller, Mar 28 2015

[2^n*Factorial(n): n in [0..35]]; // Vincenzo Librandi, Apr 22 2011

I:=[2,8]; [1] cat [n le 2 select I[n]  else (3*n-1)*Self(n-1)-2*(n-1)^2*Self(n-2): n in [1..35] ]; // Vincenzo Librandi, Feb 19 2015

A000165 := proc(n) option remember; if n <= 1 then 1 else n*A000165(n-2); fi; end;
ZL:=[S, {a = Atom, b = Atom, S = Prod(X,Sequence(Prod(X,b))), X = Sequence(b,card >= 0)}, labelled]: seq(combstruct[count](ZL, size=n), n=0..17); # Zerinvary Lajos, Mar 26 2008
G(x):=(1-2*x)^(-1): f[0]:=G(x): for n from 1 to 29 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n],n=0..17); # Zerinvary Lajos, Apr 03 2009
A000165 := proc(n) doublefactorial(2*n) ; end proc; seq(A000165(n),n=0..10) ; # R. J. Mathar, Oct 20 2009

Table[(2 n)!!, {n, 30}] (* Vladimir Joseph Stephan Orlovsky, Dec 13 2008 *)
(2 Range[0, 30])!! (* Harvey P. Dale, Jan 23 2015 *)
RecurrenceTable[{a[n] == 2 n*a[n-1], a[0] == 1}, a, {n,0,30}] (* Ray Chandler, Jul 30 2015 *)

a(n)=n!<Charles R Greathouse IV, Feb 11 2011

{a(n) = prod( k=1, n, 2*k)}; /* Michael Somos, Jan 04 2013 */

from math import factorial
def A000165(n): return factorial(n)<Chai Wah Wu, Jan 24 2023

[2^n*factorial(n) for n in range(31)] # G. C. Greubel, Jul 21 2024

E.g.f.: 1/(1-2*x).
a(n) = A001044(n)/A000142(n)*A000079(n) = Product_{i=0..n-1} (2*i+2) = 2^n*Pochhammer(1,n). - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
D-finite with recurrence a(n) = 2*n * a(n-1), n>0, a(0)=1. - Paul Barry, Aug 26 2004
This is the binomial mean transform of A001907. See Spivey and Steil (2006). - Michael Z. Spivey (mspivey(AT)ups.edu), Feb 26 2006
a(n) = Integral_{x>=0} x^n*exp(-x/2)/2 dx. - Paul Barry, Jan 28 2008
G.f.: 1/(1-2x/(1-2x/(1-4x/(1-4x/(1-6x/(1-6x/(1-.... (continued fraction). - Paul Barry, Feb 07 2009
a(n) = A006882(2*n). - R. J. Mathar, Oct 20 2009
From Gary W. Adamson, Jul 18 2011: (Start)
a(n) = upper left term in M^n, M = a production matrix (twice Pascal's triangle deleting the first "2", with the rest zeros; cf. A028326):
  2,  2,  0,  0,  0,  0, ...
  2,  4,  2,  0,  0,  0, ...
  2,  6,  6,  2,  0,  0, ...
  2,  8, 12,  8,  2,  0, ...
  2, 10, 20, 20, 10,  2, ...
  ... (End)
From Sergei N. Gladkovskii, Apr 11 2013, May 01 2013, May 24 2013, Sep 30 2013, Oct 27 2013: (Start)
Continued fractions:
G.f.: 1 + x*(Q(0) - 1)/(x+1) where Q(k) = 1 + (2*k+2)/(1-x/(x+1/Q(k+1))).
G.f.: 1/Q(0) where Q(k) = 1 + 2*k*x - 2*x*(k+1)/Q(k+1).
G.f.: G(0)/2 where G(k) = 1 + 1/(1 - x*(2*k+2)/(x*(2*k+2) + 1/G(k+1))).
G.f.: 1/Q(0) where Q(k) = 1 - x*(4*k+2) - 4*x^2*(k+1)^2/Q(k+1).
G.f.: R(0) where R(k) = 1 - x*(2*k+2)/(x*(2*k+2)-1/(1-x*(2*k+2)/(x*(2*k+2) -1/R(k+1)))). (End)
a(n) = (2n-2)*a(n-2) + (2n-1)*a(n-1), n>1. - Ivan N. Ianakiev, Aug 06 2013
From Peter Bala, Feb 18 2015: (Start)
Recurrence equation: a(n) = (3*n - 1)*a(n-1) - 2*(n - 1)^2*a(n-2) with a(1) = 2 and a(2) = 8.
The sequence b(n) = A068102(n) also satisfies this second-order recurrence. This leads to the generalized continued fraction expansion lim_{n -> oo} b(n)/a(n) = log(2) = 1/(2 - 2/(5 - 8/(8 - 18/(11 - ... - 2*(n - 1)^2/((3*n - 1) - ... ))))). (End)
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = sqrt(e) (A019774).
Sum_{n>=0} (-1)^n/a(n) = 1/sqrt(e) (A092605). (End)
Limit_{n->oo} a(n)^4 / (n * A134372(n)) = Pi. - Daniel Suteu, Apr 09 2022
a(n) = 1/([x^n] hypergeom([1], [1], x/2)). - Peter Luschny, Sep 13 2024
a(n) = Sum_{k=0..n} k!*(n-k)!*binomial(n,k)^2. - Ridouane Oudra, Jul 13 2025

A085158 Sextuple factorials, 6-factorials, n!!!!!!, n!6.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 7, 16, 27, 40, 55, 72, 91, 224, 405, 640, 935, 1296, 1729, 4480, 8505, 14080, 21505, 31104, 43225, 116480, 229635, 394240, 623645, 933120, 1339975, 3727360, 7577955, 13404160, 21827575, 33592320, 49579075, 141639680

Offset: 0

Hugo Pfoertner, Jun 21 2003

nonn

The term "Sextuple factorial numbers" is also used for the sequences A008542, A008543, A011781, A047058, A047657, A049308, which have a different definition. The definition given here is the one commonly used.

			a(14) = 224 because 14*a(14-6) = 14*a(8) = 14*16 = 224.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Multifactorial.

Cf. n!:A000142, n!!:A006882, n!!!:A007661, n!!!!:A007662, n!!!!!:A085157, 6-factorial primes: n!!!!!!+1:A085150, n!!!!!!-1:A051592.

Cf. A288093.

a:= function(n)
    if n<1 then return 1;
    else return n*a(n-6);
    fi;
  end;
List([0..40], n-> a(n) ); # G. C. Greubel, Aug 21 2019

b:=func< n | n le 6 select n else n*Self(n-6) >;
[1] cat [b(n): n in [1..40]]; // G. C. Greubel, Aug 21 2019

a:= n-> `if`(n<1, 1, n*a(n-6)); seq(a(n), n=0..40); # G. C. Greubel, Aug 21 2019

Table[Times@@Range[n,1,-6],{n,0,40}] (* Harvey P. Dale, Aug 10 2019 *)

a(n)=if(n<1, 1, n*a(n-6));
vector(40, n, n--; a(n) ) \\ G. C. Greubel, Aug 21 2019

def a(n):
    if (n<1): return 1
    else: return n*a(n-6)
[a(n) for n in (0..40)] # G. C. Greubel, Aug 21 2019

a(n)=1 for n < 1, otherwise a(n) = n*a(n-6).

Sum_{n>=0} 1/a(n) = A288093. - Amiram Eldar, Nov 10 2020

A047053 a(n) = 4^n * n!.

Original entry on oeis.org

1, 4, 32, 384, 6144, 122880, 2949120, 82575360, 2642411520, 95126814720, 3805072588800, 167423193907200, 8036313307545600, 417888291992371200, 23401744351572787200, 1404104661094367232000, 89862698310039502848000

Offset: 0

Joe Keane (jgk(AT)jgk.org)

Original name was "Quadruple factorial numbers".
For n >= 1, a(n) is the order of the wreath product of the cyclic group C_4 and the symmetric group S_n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001
Number of n X n monomial matrices with entries 0, +/-1, +/-i.
a(n) is the product of the positive integers <= 4*n that are multiples of 4. - Peter Luschny, Jun 23 2011
Also, a(n) is the number of signed permutations of length 2*n that are equal to their reverse-complements. (See the Hardt and Troyka reference.) - Justin M. Troyka, Aug 13 2011.
Pi^n/a(n) is the volume of a 2*n-dimensional ball with radius 1/2. - Peter Luschny, Jul 24 2012
Equals the first right hand column of A167557, and also equals the first right hand column of A167569. - Johannes W. Meijer, Nov 12 2009
a(n) is the order of the group U_n(Z[i]) = {A in M_n(Z[i]): A*A^H = I_n}, the group of n X n unitary matrices over the Gaussian integers. Here A^H is the conjugate transpose of A. - Jianing Song, Mar 29 2021

			G.f. = 1 + 4*x + 32*x^2 + 384*x^3 + 6144*x^4 + 122880*x^5 + 2949120*x^6 + ...

Vincenzo Librandi, Table of n, a(n) for n = 0..100
R. B. Brent, Generalizing Tuenter's Binomial Sums, J. Int. Seq. 18 (2015), # 15.3.2.
CombOS - Combinatorial Object Server, Generate colored permutations
R. Coquereaux and J.-B. Zuber, Maps, immersions and permutations, arXiv preprint arXiv:1507.03163 [math.CO], 2015-2016. Also J. Knot Theory Ramifications 25, 1650047 (2016), DOI.
Sylvie Corteel and Lauren Williams, Tableaux Combinatorics for the Asymmetric Exclusion Process II, arXiv:0810.2916 [math.CO], 2008-2009.
A. Hardt and J. M. Troyka, Restricted symmetric signed permutations, Pure Mathematics and Applications, Vol. 23 (No. 3, 2012), pp. 179-217.
A. Hardt and J. M. Troyka, Slides (associated with the Hardt and Troyka reference above).
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 492.
Alexsandar Petojevic, The Function vM_m(s; a; z) and Some Well-Known Sequences, Journal of Integer Sequences, 5 (2002), Article 02.1.7.
M. D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications, J. Int. Seq. 13 (2010), Article 10.6.7, p. 39.
Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, 9 (2006), Article 06.1.1.

Cf. A000142, A000165, A007696, A008545, A032031, A047058, A087299, A092042, A092616.

a(n)= A051142(n+1, 0) (first column of triangle).

[4^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Jul 20 2011

A047053:= n -> mul(k, k = select(k-> k mod 4 = 0, [$1..4*n])): seq(A047053(n), n = 0.. 16); # Peter Luschny, Jun 23 2011

a[n_]:= With[{m=2n}, If[ m<0, 0, m!*SeriesCoefficient[1 +Sqrt[Pi]*x*Exp[x^2]*Erf[x], {x, 0, m}]]]; (* Michael Somos, Jan 03 2015 *)
Table[4^n n!,{n,0,20}] (* Harvey P. Dale, Sep 19 2021 *)

PARI
```
a(n)=4^n*n!;
```

a(n) = 4^n * n!.
E.g.f.: 1/(1 - 4*x).
Integral representation as the n-th moment of a positive function on a positive half-axis: a(n) = Integral_{x=0..oo} x^n*exp(-x/4)/4, n >= 0. This representation is unique. - Karol A. Penson, Jan 28 2002 [corrected by Jason Yuen, May 04 2025]
Sum_{k>=0} (-1)^k/(2*k + 1)^n = (-1)^n * n * (PolyGamma[n-1, 1/4] - PolyGamma[n-1, 3/4]) / a(n) for n > 0. - Joseph Biberstine (jrbibers(AT)indiana.edu), Jul 27 2006
a(n) = Sum_{k=0..n} C(n,k)*(2k)!*(2(n-k))!/(k!(n-k)!) = Sum_{k=0..n} C(n,k)*A001813(k)*A001813(n-k). - Paul Barry, May 04 2007
E.g.f.: With interpolated zeros, 1 + sqrt(Pi)*x*exp(x^2)*erf(x). - Paul Barry, Apr 10 2010
From Gary W. Adamson, Jul 19 2011: (Start)
a(n) = sum of top row terms of M^n, M = an infinite square production matrix as follows:
  2, 2, 0, 0, 0, 0, ...
  4, 4, 4, 0, 0, 0, ...
  6, 6, 6, 6, 0, 0, ...
  8, 8, 8, 8, 8, 0, ...
  ... (End)
G.f.: 1/(1 - 4*x/(1 - 4*x/(1 - 8*x/(1 - 8*x/(1 - 12*x/(1 - 12*x/(1 - 16*x/(1 - ... (continued fraction). - Philippe Deléham, Jan 08 2012
G.f.: 2/G(0), where G(k) = 1 + 1/(1 - 8*x*(k + 1)/(8*x*(k + 1) - 1 + 8*x*(k + 1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
G.f.: 1/Q(0), where Q(k) = 1 - 4*x*(2*k + 1) - 16*x^2*(k + 1)^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Sep 28 2013
a(n) = A000142(n) * A000302(n). - Michel Marcus, Nov 28 2013
a(n) = A087299(2*n). - Michael Somos, Jan 03 2015
D-finite with recurrence: a(n) - 4*n*a(n-1) = 0. - R. J. Mathar, Jan 27 2020
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/4) (A092042).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/4) (A092616). (End)

Edited by Karol A. Penson, Jan 22 2002

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

Original entry on oeis.org

1, 1, 7, 91, 1729, 43225, 1339975, 49579075, 2131900225, 104463111025, 5745471106375, 350473737488875, 23481740411754625, 1714167050058087625, 135419196954588922375, 11510631741140058401875, 1047467488443745314570625, 101604346379043295513350625

Offset: 0

Joe Keane (jgk(AT)jgk.org)

nonn

a(n), n>=1, enumerates increasing heptic (7-ary) trees with n vertices. - Wolfdieter Lang, Sep 14 2007; see a D. Callan comment on A007559 (number of increasing quarterny trees).

Vincenzo Librandi, Table of n, a(n) for n = 0..300
Index entries for sequences related to factorial numbers.

Cf. A085158, A034689, A034723, A034724, A034787, A034788, A004993, A047058, A047657, A051151.

Cf. k-fold factorials: A000142, A001147 (and A000165, A006882), A007559 (and A032031, A008544, A007661), A007696 (and A001813, A008545, A047053, A007662), A008548 (and A052562, A047055, A085157), A045754 (and A084947, A114799), A045755.

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

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

a := n -> mul(6*k+1, k=0..n-1);
G(x):=(1-6*x)^(-1/6): f[0]:=G(x): for n from 1 to 29 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n],n=0..15); # Zerinvary Lajos, Apr 03 2009

Table[Product[(6*k+1), {k,0,n-1}], {n,0,20}] (* Vladimir Joseph Stephan Orlovsky, Nov 08 2008, modified by G. C. Greubel, Aug 17 2019 *)
FoldList[Times, 1, 6Range[0, 20] + 1] (* Vincenzo Librandi, Jun 10 2013 *)
Table[6^n*Pochhammer[1/6, n], {n,0,20}] (* G. C. Greubel, Aug 17 2019 *)

a(n)=prod(k=1,n-1,6*k+1) \\ Charles R Greathouse IV, Jul 19 2011

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

E.g.f.: (1-6*x)^(-1/6).
a(n) ~ 2^(1/2)*Pi^(1/2)*Gamma(1/6)^-1*n^(-1/3)*6^n*e^-n*n^n*{1 + 1/72*n^-1 - ...}. - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
a(n) = Sum_{k=0..n} (-6)^(n-k)*A048994(n, k). - Philippe Deléham, Oct 29 2005
G.f.: 1+x/(1-7x/(1-6x/(1-13x/(1-12x/(1-19x/(1-18x/(1-25x/(1-24x/(1-... (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-5)^n*Sum_{k=0..n} (6/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
G.f.: 1/Q(0) where Q(k) = 1 - x*(6*k+1)/(1 - x*(6*k+6)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 20 2013
a(n) = A085158(6*n-5). - M. F. Hasler, Feb 23 2018
D-finite with recurrence: a(n) +(-6*n+5)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
Sum_{n>=0} 1/a(n) = 1 + (e/6^5)^(1/6)*(Gamma(1/6) - Gamma(1/6, 1/6)). - Amiram Eldar, Dec 18 2022

A008543 Sextuple factorial numbers: Product_{k=0..n-1} (6*k + 5).

Original entry on oeis.org

1, 5, 55, 935, 21505, 623645, 21827575, 894930575, 42061737025, 2229272062325, 131527051677175, 8549258359016375, 606997343490162625, 46738795448742522125, 3879320022245629336375, 345259481979861010937375, 32799650788086796039050625, 3312764729596766399944113125

Offset: 0

Joe Keane (jgk(AT)jgk.org)

Vincenzo Librandi, Table of n, a(n) for n = 0..300
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), Article 00.2.4.

a(n) = A013988(n+1, 1) (first column of triangle).

Cf. A004994, A034787, A049308, A047058, A048994, A051151.

[Round(6^n*Gamma(n+5/6)/Gamma(5/6)): n in [0..20]]; // G. C. Greubel, Dec 03 2019

Maple
```
f := n->product( (6*k-1),k=0..n);
```

FoldList[Times,1,6Range[0,15]+5]  (* Harvey P. Dale, Feb 20 2011 *)
Table[6^n*Pochhammer[5/6, n], {n, 0, 20}] (* G. C. Greubel, Dec 03 2019 *)
CoefficientList[Series[(1 - 6x)^(-5/6), {x, 0, 20}], x] Range[0, 20]! (* Nikolaos Pantelidis, Jan 31 2023 *)

a(n)=prod(k=1,n,6*k-1) \\ Charles R Greathouse IV, Aug 17 2011

[6^n*rising_factorial(5/6, n) for n in (0..20)] # G. C. Greubel, Dec 03 2019

a(n) = 5*A034787(n) = (6*n-1)(!^6), n >= 1, a(0) := 1.
E.g.f.: (1 - 6*x)^(-5/6).
a(n) ~ 2^(1/2)*Pi^(1/2)*Gamma(5/6)^-1*n^(1/3)*6^n*e^-n*n^n*(1 + (1/72)*n^-1 + ...). - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
G.f.: 1/(1-5x/(1-6x/(1-11x/(1-12x/(1-17x/(1-18x/(1-23x/(1-24x/(1-... (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-1)^n*Sum_{k=0..n} 6^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 - 1/Q(0))/x where Q(k) = 1 - x*(6*k-1)/(1 - x*(6*k+6)/Q(k+1)); (continued fraction). - Sergei N. Gladkovskii, Mar 20 2013
a(n) = 6^n * Gamma(n+5/6) / Gamma(5/6). - Vaclav Kotesovec, Jan 28 2015
D-finite with recurrence: a(n) +(-6*n+1)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
From Nikolaos Pantelidis, May 22 2022: (Start)
G.f.: 1/G(0), where G(k) = 1 - (12*k+5)*x - 6*(k+1)*(6*k+5)*x^2/G(k+1) (a continued fraction);
which starts 1/(1-5*x-30*x^2/(1-17*x-132*x^2/(1-29*x-306*x^2/(1-41*x-552*x^2/(1-53*x-870*x^2/(1-65*x-1260*x^2/(1-...))))))) (a Jacobi continued fraction).
(End)
Sum_{n>=0} 1/a(n) = 1 + (e/6)^(1/6)*(Gamma(5/6) - Gamma(5/6, 1/6)). - Amiram Eldar, Dec 18 2022

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)

A051188 Sept-factorial numbers.

Original entry on oeis.org

1, 7, 98, 2058, 57624, 2016840, 84707280, 4150656720, 232436776320, 14643516908160, 1025046183571200, 78928556134982400, 6629998715338521600, 603329883095805465600, 59126328543388935628800

Offset: 0

Wolfdieter Lang

For n >= 1, a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_7)^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 614.
Index entries for sequences related to factorial numbers

Cf. A000142, A000165, A032031, A045754, A047053, A047058, A049209, A051186, A052562, A053106, A084947, A092516, A092750, A144827, A144739, A147585.

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

a(n) = n!*7^n; \\ Michel Marcus, Jun 08 2020

a(n) = n!*7^n =: (7*n)(!^7).
a(n) = 7*A034834(n) = Product_{k=1..n} 7*k, n >= 1.
E.g.f.: 1/(1 - 7*x).
G.f.: 1/(1 - 7*x/(1 - 7*x/(1 - 14*x/(1 - 14*x/(1 - 21*x/(1 - 21*x/(1 - 28*x/(1 - 28*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/7) (A092516).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/7) (A092750). (End)

A051189 Octo-factorial numbers.

Original entry on oeis.org

1, 8, 128, 3072, 98304, 3932160, 188743680, 10569646080, 676457349120, 48704929136640, 3896394330931200, 342882701121945600, 32916739307706777600, 3423340888001504870400, 383414179456168545484800

Offset: 0

Wolfdieter Lang

For n >= 1, a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_8)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001
Number of n X n monomial matrices whose nonzero entries are unit quaternions.
Number of ways of reassembling n slices of toast or of binding n square pages. - Donald S. McDonald, Sep 24 2005

Vincenzo Librandi, Table of n, a(n) for n = 0..300
Index entries for sequences related to factorial numbers

Cf. A000165, A034976, A047058, A051188, A053115.

Shifted absolute values are column 1 of A051187.

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

Table[n! 8^n,{n,0,20}] (* Harvey P. Dale, Aug 14 2021 *)

[8^n*factorial(n) for n in range(40)] # G. C. Greubel, Oct 21 2022

a(n) = 8*A034976(n) = Product_{k=1..n} 8*k, n >= 1; a(0) = 1.
a(n) = n!*8^n.
E.g.f.: 1/(1-8*x).
G.f.: 1/(1 - 8*x/(1 - 8*x/(1 - 16*x/(1 - 16*x/(1 - 24*x/(1 - 24*x/(1 - 32*x/(1 - 32*x/(1 - ... (continued fraction). - Philippe Deléham, Jan 07 2012
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/8).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/8). (End)

A034689 a(n) = n-th sextic factorial number divided by 2.

Original entry on oeis.org

1, 8, 112, 2240, 58240, 1863680, 70819840, 3116072960, 155803648000, 8725004288000, 540950265856000, 36784618078208000, 2722061737787392000, 217764939022991360000, 18727784755977256960000, 1722956197549907640320000, 168849707359890948751360000

Offset: 1

Wolfdieter Lang

G. C. Greubel, Table of n, a(n) for n = 1..325
Index entries for sequences related to factorial numbers.

Cf. A007559, A008542, A034723, A034724, A034787, A034788, A047657, A047058.

[n le 1 select 1 else (6*n-4)*Self(n-1): n in [1..40]]; // G. C. Greubel, Oct 21 2022

Table[6^n*Pochhammer[1/3, n]/2, {n, 40}] (* G. C. Greubel, Oct 21 2022 *)

[6^n*rising_factorial(1/3,n)/2 for n in range(1,40)] # G. C. Greubel, Oct 21 2022

2*a(n) = (6*n-4)(!^6) = Product_{j=1..n} (6*j-4) = 2^n*A007559(n), A007559(n) = (3*n-2)(!^3) = Product_{j=1..n} (3*j-2).
E.g.f.: (-1 + (1-6*x)^(-1/3))/2.
D-finite with recurrence: a(n) = 2*(3*n-2)*a(n-1). - R. J. Mathar, Feb 24 2020
a(n) = 3*6^(n-1)*Pochhammer(n, 1/3). - G. C. Greubel, Oct 21 2022
From Amiram Eldar, Dec 18 2022: (Start)
a(n) = A047657(n)/2.
Sum_{n>=1} 1/a(n) = 2*(e/6^4)^(1/6)*(Gamma(1/3, 1/6) - Gamma(1/3)). (End)

A114799 Septuple factorial, 7-factorial, n!7, n!!!!!!!, a(n) = n*a(n-7) if n > 1, else 1.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 7, 8, 18, 30, 44, 60, 78, 98, 120, 288, 510, 792, 1140, 1560, 2058, 2640, 6624, 12240, 19800, 29640, 42120, 57624, 76560, 198720, 379440, 633600, 978120, 1432080, 2016840, 2756160, 7352640, 14418720, 24710400, 39124800

Offset: 0

Jonathan Vos Post, Feb 18 2006

Many of the terms yield multifactorial primes a(n) + 1, e.g.: a(2) + 1 = 3, a(4) + 1 = 5, a(6) + 1 = 7, a(9) + 1 = 19, a(10) + 1 = 31, a(12) + 1 = 61, a(13) + 1 = 79, a(24) + 1 = 12241, a(25) + 1 = 19801, a(26) + 1 = 29641, a(29) + 1 = 76561, a(31) + 1 = 379441, a(35) + 1 = 2016841, a(36) + 1 = 2756161, ...

Equivalently, product of all positive integers <= n congruent to n (mod 7). - M. F. Hasler, Feb 23 2018

			a(40) = 40 * a(40-7) = 40 * a(33) = 40 * (33*a(26)) = 40 * 33 * (26*a(19)) = 40 * 33 * 26 * (19*a(12)) = 40 * 33 * 26 * 19 * (12*a(5)) = 40 * 33 * 26 * 19 * 12 5 = 39124800.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Multifactorial.
Index entries for sequences related to factorial numbers

Cf. A045754, A084947, A288094.

a:= function(n)
    if n<1 then return 1;
    else return n*a(n-7);
    fi;
  end;
List([0..40], n-> a(n) ); # G. C. Greubel, Aug 20 2019

b:= func< n | (n lt 8) select n else n*Self(n-7) >;
[1] cat [b(n): n in [1..40]]; // G. C. Greubel, Aug 20 2019

A114799 := proc(n)
    option remember;
    if n < 1 then
        1;
    else
        n*procname(n-7) ;
    end if;
end proc:
seq(A114799(n),n=0..40) ; # R. J. Mathar, Jun 23 2014
A114799 := n -> product(n-7*k,k=0..(n-1)/7); # M. F. Hasler, Feb 23 2018

a[n_]:= If[n<1, 1, n*a[n-7]]; Table[a[n], {n,0,40}] (* G. C. Greubel, Aug 20 2019 *)

A114799(n,k=7)=prod(j=0,(n-1)\k,n-j*k) \\ M. F. Hasler, Feb 23 2018

def a(n):
    if (n<1): return 1
    else: return n*a(n-7)
[a(n) for n in (0..40)] # G. C. Greubel, Aug 20 2019

a(n) = 1 for n <= 1, else a(n) = n*a(n-7).

Sum_{n>=0} 1/a(n) = A288094. - Amiram Eldar, Nov 10 2020

Edited by M. F. Hasler, Feb 23 2018

cp's OEIS Frontend

A000165 Double factorial of even numbers: (2n)!! = 2^n*n!.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

PARI

Python

SageMath

Formula

A085158 Sextuple factorials, 6-factorials, n!!!!!!, n!6.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A047053 a(n) = 4^n * n!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A008543 Sextuple factorial numbers: Product_{k=0..n-1} (6*k + 5).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

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.