cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 47 results. Next

A228512 Numbers n such that first n digits of A091725 (decimal expansion of second exponential integral at 1, ExpIntegralEi[1]) gives a prime number.

Original entry on oeis.org

26, 63, 323, 784
Offset: 1

Views

Author

Robert Price, Aug 23 2013

Keywords

Comments

a(5) > 10^4.

Crossrefs

Cf. A091725.

A000142 Factorial numbers: n! = 1*2*3*4*...*n (order of symmetric group S_n, number of permutations of n letters).

Original entry on oeis.org

1, 1, 2, 6, 24, 120, 720, 5040, 40320, 362880, 3628800, 39916800, 479001600, 6227020800, 87178291200, 1307674368000, 20922789888000, 355687428096000, 6402373705728000, 121645100408832000, 2432902008176640000, 51090942171709440000, 1124000727777607680000
Offset: 0

Views

Author

Keywords

Comments

The earliest publication that discusses this sequence appears to be the Sepher Yezirah [Book of Creation], circa AD 300. (See Knuth, also the Zeilberger link.) - N. J. A. Sloane, Apr 07 2014
For n >= 1, a(n) is the number of n X n (0,1) matrices with each row and column containing exactly one entry equal to 1.
This sequence is the BinomialMean transform of A000354. (See A075271 for definition.) - John W. Layman, Sep 12 2002 [This is easily verified from the Paul Barry formula for A000354, by interchanging summations and using the formula: Sum_k (-1)^k C(n-i, k) = KroneckerDelta(i,n). - David Callan, Aug 31 2003]
Number of distinct subsets of T(n-1) elements with 1 element A, 2 elements B, ..., n - 1 elements X (e.g., at n = 5, we consider the distinct subsets of ABBCCCDDDD and there are 5! = 120). - Jon Perry, Jun 12 2003
n! is the smallest number with that prime signature. E.g., 720 = 2^4 * 3^2 * 5. - Amarnath Murthy, Jul 01 2003
a(n) is the permanent of the n X n matrix M with M(i, j) = 1. - Philippe Deléham, Dec 15 2003
Given n objects of distinct sizes (e.g., areas, volumes) such that each object is sufficiently large to simultaneously contain all previous objects, then n! is the total number of essentially different arrangements using all n objects. Arbitrary levels of nesting of objects are permitted within arrangements. (This application of the sequence was inspired by considering leftover moving boxes.) If the restriction exists that each object is able or permitted to contain at most one smaller (but possibly nested) object at a time, the resulting sequence begins 1,2,5,15,52 (Bell Numbers?). Sets of nested wooden boxes or traditional nested Russian dolls come to mind here. - Rick L. Shepherd, Jan 14 2004
From Michael Somos, Mar 04 2004; edited by M. F. Hasler, Jan 02 2015: (Start)
Stirling transform of [2, 2, 6, 24, 120, ...] is A052856 = [2, 2, 4, 14, 76, ...].
Stirling transform of [1, 2, 6, 24, 120, ...] is A000670 = [1, 3, 13, 75, ...].
Stirling transform of [0, 2, 6, 24, 120, ...] is A052875 = [0, 2, 12, 74, ...].
Stirling transform of [1, 1, 2, 6, 24, 120, ...] is A000629 = [1, 2, 6, 26, ...].
Stirling transform of [0, 1, 2, 6, 24, 120, ...] is A002050 = [0, 1, 5, 25, 140, ...].
Stirling transform of (A165326*A089064)(1...) = [1, 0, 1, -1, 8, -26, 194, ...] is [1, 1, 2, 6, 24, 120, ...] (this sequence). (End)
First Eulerian transform of 1, 1, 1, 1, 1, 1... The first Eulerian transform transforms a sequence s to a sequence t by the formula t(n) = Sum_{k=0..n} e(n, k)s(k), where e(n, k) is a first-order Eulerian number [A008292]. - Ross La Haye, Feb 13 2005
Conjecturally, 1, 6, and 120 are the only numbers which are both triangular and factorial. - Christopher M. Tomaszewski (cmt1288(AT)comcast.net), Mar 30 2005
n! is the n-th finite difference of consecutive n-th powers. E.g., for n = 3, [0, 1, 8, 27, 64, ...] -> [1, 7, 19, 37, ...] -> [6, 12, 18, ...] -> [6, 6, ...]. - Bryan Jacobs (bryanjj(AT)gmail.com), Mar 31 2005
a(n+1) = (n+1)! = 1, 2, 6, ... has e.g.f. 1/(1-x)^2. - Paul Barry, Apr 22 2005
Write numbers 1 to n on a circle. Then a(n) = sum of the products of all n - 2 adjacent numbers. E.g., a(5) = 1*2*3 + 2*3*4 + 3*4*5 + 4*5*1 +5*1*2 = 120. - Amarnath Murthy, Jul 10 2005
The number of chains of maximal length in the power set of {1, 2, ..., n} ordered by the subset relation. - Rick L. Shepherd, Feb 05 2006
The number of circular permutations of n letters for n >= 0 is 1, 1, 1, 2, 6, 24, 120, 720, 5040, 40320, ... - Xavier Noria (fxn(AT)hashref.com), Jun 04 2006
a(n) is the number of deco polyominoes of height n (n >= 1; see definitions in the Barcucci et al. references). - Emeric Deutsch, Aug 07 2006
a(n) is the number of partition tableaux of size n. See Steingrimsson/Williams link for the definition. - David Callan, Oct 06 2006
Consider the n! permutations of the integer sequence [n] = 1, 2, ..., n. The i-th permutation consists of ncycle(i) permutation cycles. Then, if the Sum_{i=1..n!} 2^ncycle(i) runs from 1 to n!, we have Sum_{i=1..n!} 2^ncycle(i) = (n+1)!. E.g., for n = 3 we have ncycle(1) = 3, ncycle(2) = 2, ncycle(3) = 1, ncycle(4) = 2, ncycle(5) = 1, ncycle(6) = 2 and 2^3 + 2^2 + 2^1 + 2^2 + 2^1 + 2^2 = 8 + 4 + 2 + 4 + 2 + 4 = 24 = (n+1)!. - Thomas Wieder, Oct 11 2006
a(n) is the number of set partitions of {1, 2, ..., 2n - 1, 2n} into blocks of size 2 (perfect matchings) in which each block consists of one even and one odd integer. For example, a(3) = 6 counts 12-34-56, 12-36-45, 14-23-56, 14-25-36, 16-23-45, 16-25-34. - David Callan, Mar 30 2007
Consider the multiset M = [1, 2, 2, 3, 3, 3, 4, 4, 4, 4, ...] = [1, 2, 2, ..., n x 'n'] and form the set U (where U is a set in the strict sense) of all subsets N (where N may be a multiset again) of M. Then the number of elements |U| of U is equal to (n+1)!. E.g. for M = [1, 2, 2] we get U = [[], [2], [2, 2], [1], [1, 2], [1, 2, 2]] and |U| = 3! = 6. This observation is a more formal version of the comment given already by Rick L. Shepherd, Jan 14 2004. - Thomas Wieder, Nov 27 2007
For n >= 1, a(n) = 1, 2, 6, 24, ... are the positions corresponding to the 1's in decimal expansion of Liouville's constant (A012245). - Paul Muljadi, Apr 15 2008
Triangle A144107 has n! for row sums (given n > 0) with right border n! and left border A003319, the INVERTi transform of (1, 2, 6, 24, ...). - Gary W. Adamson, Sep 11 2008
Equals INVERT transform of A052186 and row sums of triangle A144108. - Gary W. Adamson, Sep 11 2008
From Abdullahi Umar, Oct 12 2008: (Start)
a(n) is also the number of order-decreasing full transformations (of an n-chain).
a(n-1) is also the number of nilpotent order-decreasing full transformations (of an n-chain). (End)
n! is also the number of optimal broadcast schemes in the complete graph K_{n}, equivalent to the number of binomial trees embedded in K_{n} (see Calin D. Morosan, Information Processing Letters, 100 (2006), 188-193). - Calin D. Morosan (cd_moros(AT)alumni.concordia.ca), Nov 28 2008
Let S_{n} denote the n-star graph. The S_{n} structure consists of n S_{n-1} structures. This sequence gives the number of edges between the vertices of any two specified S_{n+1} structures in S_{n+2} (n >= 1). - K.V.Iyer, Mar 18 2009
Chromatic invariant of the sun graph S_{n-2}.
It appears that a(n+1) is the inverse binomial transform of A000255. - Timothy Hopper, Aug 20 2009
a(n) is also the determinant of a square matrix, An, whose coefficients are the reciprocals of beta function: a{i, j} = 1/beta(i, j), det(An) = n!. - Enrique Pérez Herrero, Sep 21 2009
The asymptotic expansions of the exponential integrals E(x, m = 1, n = 1) ~ exp(-x)/x*(1 - 1/x + 2/x^2 - 6/x^3 + 24/x^4 + ...) and E(x, m = 1, n = 2) ~ exp(-x)/x*(1 - 2/x + 6/x^2 - 24/x^3 + ...) lead to the factorial numbers. See A163931 and A130534 for more information. - Johannes W. Meijer, Oct 20 2009
Satisfies A(x)/A(x^2), A(x) = A173280. - Gary W. Adamson, Feb 14 2010
a(n) = G^n where G is the geometric mean of the first n positive integers. - Jaroslav Krizek, May 28 2010
Increasing colored 1-2 trees with choice of two colors for the rightmost branch of nonleaves. - Wenjin Woan, May 23 2011
Number of necklaces with n labeled beads of 1 color. - Robert G. Wilson v, Sep 22 2011
The sequence 1!, (2!)!, ((3!)!)!, (((4!)!)!)!, ..., ((...(n!)!)...)! (n times) grows too rapidly to have its own entry. See Hofstadter.
The e.g.f. of 1/a(n) = 1/n! is BesselI(0, 2*sqrt(x)). See Abramowitz-Stegun, p. 375, 9.3.10. - Wolfdieter Lang, Jan 09 2012
a(n) is the length of the n-th row which is the sum of n-th row in triangle A170942. - Reinhard Zumkeller, Mar 29 2012
Number of permutations of elements 1, 2, ..., n + 1 with a fixed element belonging to a cycle of length r does not depend on r and equals a(n). - Vladimir Shevelev, May 12 2012
a(n) is the number of fixed points in all permutations of 1, ..., n: in all n! permutations, 1 is first exactly (n-1)! times, 2 is second exactly (n-1)! times, etc., giving (n-1)!*n = n!. - Jon Perry, Dec 20 2012
For n >= 1, a(n-1) is the binomial transform of A000757. See Moreno-Rivera. - Luis Manuel Rivera Martínez, Dec 09 2013
Each term is divisible by its digital root (A010888). - Ivan N. Ianakiev, Apr 14 2014
For m >= 3, a(m-2) is the number hp(m) of acyclic Hamiltonian paths in a simple graph with m vertices, which is complete except for one missing edge. For m < 3, hp(m)=0. - Stanislav Sykora, Jun 17 2014
a(n) is the number of increasing forests with n nodes. - Brad R. Jones, Dec 01 2014
The factorial numbers can be calculated by means of the recurrence n! = (floor(n/2)!)^2 * sf(n) where sf(n) are the swinging factorials A056040. This leads to an efficient algorithm if sf(n) is computed via prime factorization. For an exposition of this algorithm see the link below. - Peter Luschny, Nov 05 2016
Treeshelves are ordered (plane) binary (0-1-2) increasing trees where the nodes of outdegree 1 come in 2 colors. There are n! treeshelves of size n, and classical Françon's bijection maps bijectively treeshelves into permutations. - Sergey Kirgizov, Dec 26 2016
Satisfies Benford's law [Diaconis, 1977; Berger-Hill, 2017] - N. J. A. Sloane, Feb 07 2017
a(n) = Sum((d_p)^2), where d_p is the number of standard tableaux in the Ferrers board of the integer partition p and summation is over all integer partitions p of n. Example: a(3) = 6. Indeed, the partitions of 3 are [3], [2,1], and [1,1,1], having 1, 2, and 1 standard tableaux, respectively; we have 1^2 + 2^2 + 1^2 = 6. - Emeric Deutsch, Aug 07 2017
a(n) is the n-th derivative of x^n. - Iain Fox, Nov 19 2017
a(n) is the number of maximum chains in the n-dimensional Boolean cube {0,1}^n in respect to the relation "precedes". It is defined as follows: for arbitrary vectors u, v of {0,1}^n, such that u = (u_1, u_2, ..., u_n) and v = (v_1, v_2, ..., v_n), "u precedes v" if u_i <= v_i, for i=1, 2, ..., n. - Valentin Bakoev, Nov 20 2017
a(n) is the number of shortest paths (for example, obtained by Breadth First Search) between the nodes (0,0,...,0) (i.e., the all-zeros vector) and (1,1,...,1) (i.e., the all-ones vector) in the graph H_n, corresponding to the n-dimensional Boolean cube {0,1}^n. The graph is defined as H_n = (V_n, E_n), where V_n is the set of all vectors of {0,1}^n, and E_n contains edges formed by each pair adjacent vectors. - Valentin Bakoev, Nov 20 2017
a(n) is also the determinant of the symmetric n X n matrix M defined by M(i,j) = sigma(gcd(i,j)) for 1 <= i,j <= n. - Bernard Schott, Dec 05 2018
a(n) is also the number of inversion sequences of length n. A length n inversion sequence e_1, e_2, ..., e_n is a sequence of n integers such that 0 <= e_i < i. - Juan S. Auli, Oct 14 2019
The term "factorial" ("factorielle" in French) was coined by the French mathematician Louis François Antoine Arbogast (1759-1803) in 1800. The notation "!" was first used by the French mathematician Christian Kramp (1760-1826) in 1808. - Amiram Eldar, Apr 16 2021
Also the number of signotopes of rank 2, i.e., mappings X:{{1..n} choose 2}->{+,-} such that for any three indices a < b < c, the sequence X(a,b), X(a,c), X(b,c) changes its sign at most once (see Felsner-Weil reference). - Manfred Scheucher, Feb 09 2022
a(n) is also the number of labeled commutative semisimple rings with n elements. As an example the only commutative semisimple rings with 4 elements are F_4 and F_2 X F_2. They both have exactly 2 automorphisms, hence a(4)=24/2+24/2=24. - Paul Laubie, Mar 05 2024
a(n) is the number of extremely unlucky Stirling permutations of order n+1; i.e., the number of Stirling permutations of order n+1 that have exactly one lucky car. - Bridget Tenner, Apr 09 2024

Examples

			There are 3! = 1*2*3 = 6 ways to arrange 3 letters {a, b, c}, namely abc, acb, bac, bca, cab, cba.
Let n = 2. Consider permutations of {1, 2, 3}. Fix element 3. There are a(2) = 2 permutations in each of the following cases: (a) 3 belongs to a cycle of length 1 (permutations (1, 2, 3) and (2, 1, 3)); (b) 3 belongs to a cycle of length 2 (permutations (3, 2, 1) and (1, 3, 2)); (c) 3 belongs to a cycle of length 3 (permutations (2, 3, 1) and (3, 1, 2)). - _Vladimir Shevelev_, May 13 2012
G.f. = 1 + x + 2*x^2 + 6*x^3 + 24*x^4 + 120*x^5 + 720*x^6 + 5040*x^7 + ...
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 833.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 125; also p. 90, ex. 3.
  • Florian Cajori, A History of Mathematical Notations, Dover edition (2012), pars. 448-449.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 64-66.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §4.1 Symbols Galore, p. 106.
  • Douglas R. Hofstadter, Fluid concepts & creative analogies: computer models of the fundamental mechanisms of thought, Basic Books, 1995, pages 44-46.
  • A. N. Khovanskii. The Application of Continued Fractions and Their Generalizations to Problem in Approximation Theory. Groningen: Noordhoff, Netherlands, 1963. See p. 141 (10.19).
  • D. E. Knuth, The Art of Computer Programming, Vol. 3, Section 5.1.2, p. 23. [From N. J. A. Sloane, Apr 07 2014]
  • J.-M. De Koninck and A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problème 693 pp. 90, 297, Ellipses Paris 2004.
  • 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.
  • R. W. Robinson, Counting arrangements of bishops, pp. 198-214 of Combinatorial Mathematics IV (Adelaide 1975), Lect. Notes Math., 560 (1976).
  • Sepher Yezirah [Book of Creation], circa AD 300. See verse 52.
  • 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).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 2, pages 19-24.
  • D. Stanton and D. White, Constructive Combinatorics, Springer, 1986; see p. 91.
  • Carlo Suares, Sepher Yetsira, Shambhala Publications, 1976. See verse 52.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, 1987, pp. 102.

Crossrefs

Factorial base representation: A007623.
Complement of A063992. - Reinhard Zumkeller, Oct 11 2008
Cf. A053657, A163176. - Jonathan Sondow, Jul 26 2009
Cf. A173280. - Gary W. Adamson, Feb 14 2010
Boustrophedon transforms: A230960, A230961.
Cf. A233589.
Cf. A245334.
A row of the array in A249026.
Cf. A001013 (multiplicative closure).
For factorials with initial digit d (1 <= d <= 9) see A045509, A045510, A045511, A045516, A045517, A045518, A282021, A045519; A045520, A045521, A045522, A045523, A045524, A045525, A045526, A045527, A045528, A045529.

Programs

  • Axiom
    [factorial(n) for n in 0..10]
    
  • GAP
    List([0..22],Factorial); # Muniru A Asiru, Dec 05 2018
    
  • Haskell
    a000142 :: (Enum a, Num a, Integral t) => t -> a
    a000142 n = product [1 .. fromIntegral n]
    a000142_list = 1 : zipWith (*) [1..] a000142_list
    -- Reinhard Zumkeller, Mar 02 2014, Nov 02 2011, Apr 21 2011
    
  • Julia
    print([factorial(big(n)) for n in 0:28]) # Paul Muljadi, May 01 2024
  • Magma
    a:= func< n | Factorial(n) >; [ a(n) : n in [0..10]];
    
  • Maple
    A000142 := n -> n!; seq(n!,n=0..20);
    spec := [ S, {S=Sequence(Z) }, labeled ]; seq(combstruct[count](spec,size=n), n=0..20);
    # Maple program for computing cycle indices of symmetric groups
    M:=6: f:=array(0..M): f[0]:=1: print(`n= `,0); print(f[0]); f[1]:=x[1]: print(`n= `, 1); print(f[1]); for n from 2 to M do f[n]:=expand((1/n)*add( x[l]*f[n-l],l=1..n)); print(`n= `, n); print(f[n]); od:
    with(combstruct):ZL0:=[S,{S=Set(Cycle(Z,card>0))},labeled]: seq(count(ZL0,size=n),n=0..20); # Zerinvary Lajos, Sep 26 2007
  • Mathematica
    Table[Factorial[n], {n, 0, 20}] (* Stefan Steinerberger, Mar 30 2006 *)
    FoldList[#1 #2 &, 1, Range@ 20] (* Robert G. Wilson v, May 07 2011 *)
    Range[20]! (* Harvey P. Dale, Nov 19 2011 *)
    RecurrenceTable[{a[n] == n*a[n - 1], a[0] == 1}, a, {n, 0, 22}] (* Ray Chandler, Jul 30 2015 *)
  • PARI
    a(n)=prod(i=1, n, i) \\ Felix Fröhlich, Aug 17 2014
    
  • PARI
    {a(n) = if(n<0, 0, n!)}; /* Michael Somos, Mar 04 2004 */
    
  • Python
    for i in range(1, 1000):
        y = i
        for j in range(1, i):
           y *= i - j
        print(y, "\n")
    
  • Python
    import math
    for i in range(1, 1000):
        math.factorial(i)
        print("")
    # Ruskin Harding, Feb 22 2013
    
  • Sage
    [factorial(n) for n in (1..22)] # Giuseppe Coppoletta, Dec 05 2014
    
  • Scala
    (1: BigInt).to(24: BigInt).scanLeft(1: BigInt)( * ) // Alonso del Arte, Mar 02 2019
    

Formula

Exp(x) = Sum_{m >= 0} x^m/m!. - Mohammad K. Azarian, Dec 28 2010
Sum_{i=0..n} (-1)^i * i^n * binomial(n, i) = (-1)^n * n!. - Yong Kong (ykong(AT)curagen.com), Dec 26 2000
Sum_{i=0..n} (-1)^i * (n-i)^n * binomial(n, i) = n!. - Peter C. Heinig (algorithms(AT)gmx.de), Apr 10 2007
The sequence trivially satisfies the recurrence a(n+1) = Sum_{k=0..n} binomial(n,k) * a(k)*a(n-k). - Robert FERREOL, Dec 05 2009
D-finite with recurrence: a(n) = n*a(n-1), n >= 1. n! ~ sqrt(2*Pi) * n^(n+1/2) / e^n (Stirling's approximation).
a(0) = 1, a(n) = subs(x = 1, (d^n/dx^n)(1/(2-x))), n = 1, 2, ... - Karol A. Penson, Nov 12 2001
E.g.f.: 1/(1-x). - Michael Somos, Mar 04 2004
a(n) = Sum_{k=0..n} (-1)^(n-k)*A000522(k)*binomial(n, k) = Sum_{k=0..n} (-1)^(n-k)*(x+k)^n*binomial(n, k). - Philippe Deléham, Jul 08 2004
Binomial transform of A000166. - Ross La Haye, Sep 21 2004
a(n) = Sum_{i=1..n} ((-1)^(i-1) * sum of 1..n taken n - i at a time) - e.g., 4! = (1*2*3 + 1*2*4 + 1*3*4 + 2*3*4) - (1*2 + 1*3 + 1*4 + 2*3 + 2*4 + 3*4) + (1 + 2 + 3 + 4) - 1 = (6 + 8 + 12 + 24) - (2 + 3 + 4 + 6 + 8 + 12) + 10 - 1 = 50 - 35 + 10 - 1 = 24. - Jon Perry, Nov 14 2005
a(n) = (n-1)*(a(n-1) + a(n-2)), n >= 2. - Matthew J. White, Feb 21 2006
1 / a(n) = determinant of matrix whose (i,j) entry is (i+j)!/(i!(j+1)!) for n > 0. This is a matrix with Catalan numbers on the diagonal. - Alexander Adamchuk, Jul 04 2006
Hankel transform of A074664. - Philippe Deléham, Jun 21 2007
For n >= 2, a(n-2) = (-1)^n*Sum_{j=0..n-1} (j+1)*Stirling1(n,j+1). - Milan Janjic, Dec 14 2008
From Paul Barry, Jan 15 2009: (Start)
G.f.: 1/(1-x-x^2/(1-3x-4x^2/(1-5x-9x^2/(1-7x-16x^2/(1-9x-25x^2... (continued fraction), hence Hankel transform is A055209.
G.f. of (n+1)! is 1/(1-2x-2x^2/(1-4x-6x^2/(1-6x-12x^2/(1-8x-20x^2... (continued fraction), hence Hankel transform is A059332. (End)
a(n) = Product_{p prime} p^(Sum_{k > 0} floor(n/p^k)) by Legendre's formula for the highest power of a prime dividing n!. - Jonathan Sondow, Jul 24 2009
a(n) = A053657(n)/A163176(n) for n > 0. - Jonathan Sondow, Jul 26 2009
It appears that a(n) = (1/0!) + (1/1!)*n + (3/2!)*n*(n-1) + (11/3!)*n*(n-1)*(n-2) + ... + (b(n)/n!)*n*(n-1)*...*2*1, where a(n) = (n+1)! and b(n) = A000255. - Timothy Hopper, Aug 12 2009
Sum_{n >= 0} 1/a(n) = e. - Jaume Oliver Lafont, Mar 03 2009
a(n) = a(n-1)^2/a(n-2) + a(n-1), n >= 2. - Jaume Oliver Lafont, Sep 21 2009
a(n) = Gamma(n+1). - Enrique Pérez Herrero, Sep 21 2009
a(n) = A173333(n,1). - Reinhard Zumkeller, Feb 19 2010
a(n) = A_{n}(1) where A_{n}(x) are the Eulerian polynomials. - Peter Luschny, Aug 03 2010
a(n) = n*(2*a(n-1) - (n-1)*a(n-2)), n > 1. - Gary Detlefs, Sep 16 2010
1/a(n) = -Sum_{k=1..n+1} (-2)^k*(n+k+2)*a(k)/(a(2*k+1)*a(n+1-k)). - Groux Roland, Dec 08 2010
From Vladimir Shevelev, Feb 21 2011: (Start)
a(n) = Product_{p prime, p <= n} p^(Sum_{i >= 1} floor(n/p^i)).
The infinitary analog of this formula is: a(n) = Product_{q terms of A050376 <= n} q^((n)_q), where (n)_q denotes the number of those numbers <= n for which q is an infinitary divisor (for the definition see comment in A037445). (End)
The terms are the denominators of the expansion of sinh(x) + cosh(x). - Arkadiusz Wesolowski, Feb 03 2012
G.f.: 1 / (1 - x / (1 - x / (1 - 2*x / (1 - 2*x / (1 - 3*x / (1 - 3*x / ... )))))). - Michael Somos, May 12 2012
G.f. 1 + x/(G(0)-x) where G(k) = 1 - (k+1)*x/(1 - x*(k+2)/G(k+1)); (continued fraction, 2-step). - Sergei N. Gladkovskii, Aug 14 2012
G.f.: W(1,1;-x)/(W(1,1;-x) - x*W(1,2;-x)), where W(a,b,x) = 1 - a*b*x/1! + a*(a+1)*b*(b+1)*x^2/2! - ... + a*(a+1)*...*(a+n-1)*b*(b+1)*...*(b+n-1)*x^n/n! + ...; see [A. N. Khovanskii, p. 141 (10.19)]. - Sergei N. Gladkovskii, Aug 15 2012
From Sergei N. Gladkovskii, Dec 26 2012: (Start)
G.f.: A(x) = 1 + x/(G(0) - x) where G(k) = 1 + (k+1)*x - x*(k+2)/G(k+1); (continued fraction).
Let B(x) be the g.f. for A051296, then A(x) = 2 - 1/B(x). (End)
G.f.: 1 + x*(G(0) - 1)/(x-1) where G(k) = 1 - (2*k+1)/(1-x/(x - 1/(1 - (2*k+2)/(1-x/(x - 1/G(k+1) ))))); (continued fraction). - Sergei N. Gladkovskii, Jan 15 2013
G.f.: 1 + x*(1 - G(0))/(sqrt(x)-x) where G(k) = 1 - (k+1)*sqrt(x)/(1-sqrt(x)/(sqrt(x)-1/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 25 2013
G.f.: 1 + x/G(0) where G(k) = 1 - x*(k+2)/( 1 - x*(k+1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 23 2013
a(n) = det(S(i+1, j), 1 <= i, j <=n ), where S(n,k) are Stirling numbers of the second kind. - Mircea Merca, Apr 04 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(k+1)/(x*(k+1) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 24 2013
G.f.: 2/G(0), where G(k) = 1 + 1/(1 - 1/(1 - 1/(2*x*(k+1)) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 29 2013
G.f.: G(0), where G(k) = 1 + x*(2*k+1)/(1 - x*(2*k+2)/(x*(2*k+2) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 07 2013
a(n) = P(n-1, floor(n/2)) * floor(n/2)! * (n - (n-2)*((n+1) mod 2)), where P(n, k) are the k-permutations of n objects, n > 0. - Wesley Ivan Hurt, Jun 07 2013
a(n) = a(n-2)*(n-1)^2 + a(n-1), n > 1. - Ivan N. Ianakiev, Jun 18 2013
a(n) = a(n-2)*(n^2-1) - a(n-1), n > 1. - Ivan N. Ianakiev, Jun 30 2013
G.f.: 1 + x/Q(0), m=+2, where Q(k) = 1 - 2*x*(2*k+1) - m*x^2*(k+1)*(2*k+1)/( 1 - 2*x*(2*k+2) - m*x^2*(k+1)*(2*k+3)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Sep 24 2013
a(n) = A245334(n,n). - Reinhard Zumkeller, Aug 31 2014
a(n) = Product_{i = 1..n} A014963^floor(n/i) = Product_{i = 1..n} A003418(floor(n/i)). - Matthew Vandermast, Dec 22 2014
a(n) = round(Sum_{k>=1} log(k)^n/k^2), for n>=1, which is related to the n-th derivative of the Riemann zeta function at x=2 as follows: round((-1)^n * zeta^(n)(2)). Also see A073002. - Richard R. Forberg, Dec 30 2014
a(n) ~ Sum_{j>=0} j^n/e^j, where e = A001113. When substituting a generic variable for "e" this infinite sum is related to Eulerian polynomials. See A008292. This approximation of n! is within 0.4% at n = 2. See A255169. Accuracy, as a percentage, improves rapidly for larger n. - Richard R. Forberg, Mar 07 2015
a(n) = Product_{k=1..n} (C(n+1, 2)-C(k, 2))/(2*k-1); see Masanori Ando link. - Michel Marcus, Apr 17 2015
Sum_{n>=0} a(n)/(a(n + 1)*a(n + 2)) = Sum_{n>=0} 1/((n + 2)*(n + 1)^2*a(n)) = 2 - exp(1) - gamma + Ei(1) = 0.5996203229953..., where gamma = A001620, Ei(1) = A091725. - Ilya Gutkovskiy, Nov 01 2016
a(2^n) = 2^(2^n - 1) * 1!! * 3!! * 7!! * ... * (2^n - 1)!!. For example, 16! = 2^15*(1*3)*(1*3*5*7)*(1*3*5*7*9*11*13*15) = 20922789888000. - Peter Bala, Nov 01 2016
a(n) = sum(prod(B)), where the sum is over all subsets B of {1,2,...,n-1} and where prod(B) denotes the product of all the elements of set B. If B is a singleton set with element b, then we define prod(B)=b, and, if B is the empty set, we define prod(B) to be 1. For example, a(4)=(1*2*3)+(1*2)+(1*3)+(2*3)+(1)+(2)+(3)+1=24. - Dennis P. Walsh, Oct 23 2017
Sum_{n >= 0} 1/(a(n)*(n+2)) = 1. - Multiplying the denominator by (n+2) in Jaume Oliver Lafont's entry above creates a telescoping sum. - Fred Daniel Kline, Nov 08 2020
O.g.f.: Sum_{k >= 0} k!*x^k = Sum_{k >= 0} (k+y)^k*x^k/(1 + (k+y)*x)^(k+1) for arbitrary y. - Peter Bala, Mar 21 2022
E.g.f.: 1/(1 + LambertW(-x*exp(-x))) = 1/(1-x), see A258773. -(1/x)*substitute(z = x*exp(-x), z*(d/dz)LambertW(-z)) = 1/(1 - x). See A075513. Proof: Use the compositional inverse (x*exp(-x))^[-1] = -LambertW(-z). See A000169 or A152917, and Richard P. Stanley: Enumerative Combinatorics, vol. 2, p. 37, eq. (5.52). - Wolfdieter Lang, Oct 17 2022
Sum_{k >= 1} 1/10^a(k) = A012245 (Liouville constant). - Bernard Schott, Dec 18 2022
From David Ulgenes, Sep 19 2023: (Start)
1/a(n) = (e/(2*Pi*n)*Integral_{x=-oo..oo} cos(x-n*arctan(x))/(1+x^2)^(n/2) dx). Proof: take the real component of Laplace's integral for 1/Gamma(x).
a(n) = Integral_{x=0..1} e^(-t)*LerchPhi(1/e, -n, t) dt. Proof: use the relationship Gamma(x+1) = Sum_{n >= 0} Integral_{t=n..n+1} e^(-t)t^x dt = Sum_{n >= 0} Integral_{t=0..1} e^(-(t+n))(t+n)^x dt and interchange the order of summation and integration.
Conjecture: a(n) = 1/(2*Pi)*Integral_{x=-oo..oo}(n+i*x+1)!/(i*x+1)-(n+i*x-1)!/(i*x-1)dx. (End)
a(n) = floor(b(n)^n / (floor(((2^b(n) + 1) / 2^n)^b(n)) mod 2^b(n))), where b(n) = (n + 1)^(n + 2) = A007778(n+1). Joint work with Mihai Prunescu. - Lorenzo Sauras Altuzarra, Oct 18 2023
a(n) = e^(Integral_{x=1..n+1} Psi(x) dx) where Psi(x) is the digamma function. - Andrea Pinos, Jan 10 2024
a(n) = Integral_{x=0..oo} e^(-x^(1/n)) dx, for n > 0. - Ridouane Oudra, Apr 20 2024
O.g.f.: N(x) = hypergeometric([1,1], [], x) = LaplaceTransform(x/(1-x))/x, satisfying x^2*N'(x) + (x-1)*N(x) + 1 = 0, with N(0) = 1. - Wolfdieter Lang, May 31 2025

A001286 Lah numbers: a(n) = (n-1)*n!/2.

Original entry on oeis.org

1, 6, 36, 240, 1800, 15120, 141120, 1451520, 16329600, 199584000, 2634508800, 37362124800, 566658892800, 9153720576000, 156920924160000, 2845499424768000, 54420176498688000, 1094805903679488000, 23112569077678080000, 510909421717094400000
Offset: 2

Views

Author

Keywords

Comments

Number of surjections from {1,...,n} to {1,...,n-1}. - Benoit Cloitre, Dec 05 2003
First Eulerian transform of 0,1,2,3,4,... . - Ross La Haye, Mar 05 2005
With offset 0 : determinant of the n X n matrix m(i,j)=(i+j+1)!/i!/j!. - Benoit Cloitre, Apr 11 2005
These numbers arise when expressing n(n+1)(n+2)...(n+k)[n+(n+1)+(n+2)+...+(n+k)] as sums of squares: n(n+1)[n+(n+1)] = 6(1+4+9+16+ ... + n^2), n(n+1)(n+2)(n+(n+1)+(n+2)) = 36(1+(1+4)+(1+4+9)+...+(1+4+9+16+ ... + n^2)), n(n+1)(n+2)(n+3)(n+(n+1)+(n+2)+(n+3)) = 240(...), ... . - Alexander R. Povolotsky, Oct 16 2006
a(n) is the number of edges in the Hasse diagram for the weak Bruhat order on the symmetric group S_n. For permutations p,q in S_n, q covers p in the weak Bruhat order if p,q differ by an adjacent transposition and q has one more inversion than p. Thus 23514 covers 23154 due to the transposition that interchanges the third and fourth entries. Cf. A002538 for the strong Bruhat order. - David Callan, Nov 29 2007
a(n) is also the number of excedances in all permutations of {1,2,...,n} (an excedance of a permutation p is a value j such p(j)>j). Proof: j is exceeded (n-1)! times by each of the numbers j+1, j+2, ..., n; now, Sum_{j=1..n} (n-j)(n-1)! = n!(n-1)/2. Example: a(3)=6 because the number of excedances of the permutations 123, 132, 312, 213, 231, 321 are 0, 1, 1, 1, 2, 1, respectively. - Emeric Deutsch, Dec 15 2008
(-1)^(n+1)*a(n) is the determinant of the n X n matrix whose (i,j)-th element is 0 for i = j, is j-1 for j>i, and j for j < i. - Michel Lagneau, May 04 2010
Row sums of the triangle in A030298. - Reinhard Zumkeller, Mar 29 2012
a(n) is the total number of ascents (descents) over all n-permutations. a(n) = Sum_{k=1..n} A008292(n,k)*k. - Geoffrey Critzer, Jan 06 2013
For m>=4, a(m-2) is the number of Hamiltonian cycles in a simple graph with m vertices which is complete, except for one edge. Proof: think of distinct round-table seatings of m persons such that persons "1" and "2" may not be neighbors; the count is (m-3)(m-2)!/2. See also A001710. - Stanislav Sykora, Jun 17 2014
Popularity of left (right) children in treeshelves. Treeshelves are ordered binary (0-1-2) increasing trees where every child is connected to its parent by a left or a right link. Popularity is the sum of a certain statistic (number of left children, in this case) over all objects of size n. See A278677, A278678 or A278679 for more definitions and examples. See A008292 for the distribution of the left (right) children in treeshelves. - Sergey Kirgizov, Dec 24 2016

Examples

			G.f. = x^2 + 6*x^3 + 36*x^4 + 240*x^5 + 1800*x^6 + 15120*x^7 + 141120*x^8 + ...
a(10) = (1+2+3+4+5+6+7+8+9)*(1*2*3*4*5*6*7*8*9) = 16329600. - _Reinhard Zumkeller_, May 15 2010
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 90, ex. 4.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 156.
  • 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.
  • John Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 44.
  • 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).

Crossrefs

A002868 is an essentially identical sequence.
Column 2 of |A008297|.
Third column (m=2) of triangle |A111596(n, m)|: matrix product of |S1|.S2 Stirling number matrices.
Cf. also A000110, A000111.

Programs

Formula

a(n) = Sum_{i=0..n-1} (-1)^(n-i-1) * i^n * binomial(n-1,i). - Yong Kong (ykong(AT)curagen.com), Dec 26 2000 [corrected by Amiram Eldar, May 02 2022]
E.g.f.: x^2/[2(1-x)^2]. - Ralf Stephan, Apr 02 2004
a(n+1) = (-1)^(n+1)*det(M_n) where M_n is the n X n matrix M_(i,j)=max(i*(i+1)/2,j*(j+1)/2). - Benoit Cloitre, Apr 03 2004
Row sums of table A051683. - Alford Arnold, Sep 29 2006
5th binomial transform of A135218: (1, 1, 1, 25, 25, 745, 3145, ...). - Gary W. Adamson, Nov 23 2007
If we define f(n,i,x) = Sum_{k=i..n} Sum_{j=i..k} binomial(k,j)*Stirling1(n,k)*Stirling2(j,i)*x^(k-j) then a(n)=(-1)^n*f(n,2,-2), (n>=2). - Milan Janjic, Mar 01 2009
a(n) = A000217(n-1)*A000142(n-1). - Reinhard Zumkeller, May 15 2010
a(n) = (n+1)!*Sum_{k=1..n-1} 1/(k^2+3*k+2). - Gary Detlefs, Sep 14 2011
Sum_{n>=2} 1/a(n) = 2*(2 - exp(1) - gamma + Ei(1)) = 1.19924064599..., where gamma = A001620 and Ei(1) = A091725. - Ilya Gutkovskiy, Nov 24 2016
a(n+1) = a(n)*n*(n+1)/(n-1). - Chai Wah Wu, Apr 11 2018
Sum_{n>=2} (-1)^n/a(n) = 2*(gamma - Ei(-1)) - 2/e, where e = A001113 and Ei(-1) = -A099285. - Amiram Eldar, May 02 2022

A099285 Decimal expansion of -Ei(-1), negated exponential integral at -1.

Original entry on oeis.org

2, 1, 9, 3, 8, 3, 9, 3, 4, 3, 9, 5, 5, 2, 0, 2, 7, 3, 6, 7, 7, 1, 6, 3, 7, 7, 5, 4, 6, 0, 1, 2, 1, 6, 4, 9, 0, 3, 1, 0, 4, 7, 2, 9, 3, 4, 0, 6, 9, 0, 8, 2, 0, 7, 5, 7, 7, 9, 7, 8, 6, 1, 3, 0, 7, 3, 5, 6, 8, 6, 9, 8, 5, 5, 9, 1, 4, 1, 5, 4, 4, 7, 2, 2, 2, 1, 0, 2, 5, 1, 0, 3, 5, 1, 3, 7, 2, 4, 9, 9, 5, 4, 7, 5, 8
Offset: 0

Views

Author

Robert G. Wilson v, Oct 08 2004

Keywords

Comments

The divergent series g(x=1,m) = 1^m*1! - 2^m*2! + 3^m*3! - 4^m*4! + ..., m=>-1, is closely related to the value of -Ei(-1). We discovered that g(x=1,m) = (-1)^m*(A040027(m) - A000110(m+1)*Ei(1,1)*exp(1)), see A163940. We observe that Ei(1,1) = E(1,1,1) = -Ei(-1) is the constant given above and that Ei(1,1)*exp(1) = A073003 is Gompertz's constant. - Johannes W. Meijer, Oct 16 2009

Examples

			0.219383934395520273677163775460121649031047293406908207577978613...
With n := 10^6, Integral_{x = 0..n} x^(n-1)/(1 + x)^n dx = 0.21938(43...). - _Peter Bala_, Jun 19 2024
		

Crossrefs

Programs

  • Maple
    Digits:=105: evalf(-Ei(-1)); evalf(Ei(1,1)); # Johannes W. Meijer, Oct 16 2009
  • Mathematica
    RealDigits[ ExpIntegralE[1, 1], 10, 105][[1]]
  • PARI
    eint1(1, 1) \\ Michel Marcus, Aug 01 2020

Formula

-Ei(-n) = Integral_{a=n..oo} ( Integral_{b=1..oo} 1/e^(a*b) db ) da , n>0 (According to Mathematica). - Mats Granvik, May 25 2013
Equals the difference between the absolute values of A239069 and A001620. - R. J. Mathar, Mar 07 2016
From Amiram Eldar, Aug 01 2020: (Start)
Equals Integral_{x=1..oo} log(x)/exp(x) dx.
Equals Integral_{x=0..oo} exp(-exp(x)) dx.
Equals Integral_{x=0..oo} x*exp(x-exp(x)) dx. (End)
From Peter Bala, Jun 17 2024: (Start)
Equals lim_{n -> oo} Integral_{x = 0..n} x^(n-1)/(1 + x)^n dx = lim_{n -> oo} ( log(n+1) + Sum_{k = 0..n-2} (-1)^(n-k-1)* binomial(n-1, k)*((n + 1)^(k+1-n) - 1)/(k + 1 - n) ).
Alternatively, equals lim_{n -> oo} Sum_{k >= n} (n/(n + 1))^k/k = lim_{n -> oo} ( log(1/(1 - x)) - Sum_{k = 1..n-1} x^k/k ), where x = n/(n+1).
More generally, for alpha > 0, -Ei(-alpha) = lim_{n -> oo} Integral_{x = 0..n/alpha} x^(n-1)/(1 + x)^n dx. (End)

Extensions

Definition corrected by Johannes W. Meijer, Jul 26 2009
Corrected Name (minus 1, not 1), Stanislav Sykora, May 18 2012

A073003 Decimal expansion of -exp(1)*Ei(-1), also called Gompertz's constant, or the Euler-Gompertz constant.

Original entry on oeis.org

5, 9, 6, 3, 4, 7, 3, 6, 2, 3, 2, 3, 1, 9, 4, 0, 7, 4, 3, 4, 1, 0, 7, 8, 4, 9, 9, 3, 6, 9, 2, 7, 9, 3, 7, 6, 0, 7, 4, 1, 7, 7, 8, 6, 0, 1, 5, 2, 5, 4, 8, 7, 8, 1, 5, 7, 3, 4, 8, 4, 9, 1, 0, 4, 8, 2, 3, 2, 7, 2, 1, 9, 1, 1, 4, 8, 7, 4, 4, 1, 7, 4, 7, 0, 4, 3, 0, 4, 9, 7, 0, 9, 3, 6, 1, 2, 7, 6, 0, 3, 4, 4, 2, 3, 7
Offset: 0

Views

Author

Robert G. Wilson v, Aug 03 2002

Keywords

Comments

0! - 1! + 2! - 3! + 4! - 5! + ... = (Borel) Sum_{n>=0} (-y)^n n! = KummerU(1,1,1/y)/y.
Decimal expansion of phi(1) where phi(x) = Integral_{t>=0} e^-t/(x+t) dt. - Benoit Cloitre, Apr 11 2003
The divergent series g(x=1,m) = 1^m*1! - 2^m*2! + 3^m*3! - 4^m*4! + ..., m => -1, is intimately related to Gompertz's constant. We discovered that g(x=1,m) = (-1)^m * (A040027(m) - A000110(m+1) * A073003) with A000110 the Bell numbers and A040027 a sequence that was published by Gould, see for more information A163940. - Johannes W. Meijer, Oct 16 2009
Named by Le Lionnais (1983) after the English self-educated mathematician and actuary Benjamin Gompertz (1779 - 1865). It was named the Euler-Gompertz constant by Finch (2003). Lagarias (2013) noted that he has not located this constant in Gompertz's writings. - Amiram Eldar, Aug 15 2020

Examples

			0.59634736232319407434107849936927937607417786015254878157348491...
With n := 10^5, Sum_{k >= 0} (n/(n + 1))^k/(n + k) = 0.5963(51...). - _Peter Bala_, Jun 19 2024
		

References

  • Bruce C. Berndt, Ramanujan's notebooks Part II, Springer, p. 171
  • Bruce C. Berndt, Ramanujan's notebooks Part I, Springer, p. 144-145.
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, pp. 303, 424-425.
  • Francois Le Lionnais, Les nombres remarquables, Paris: Hermann, 1983. See p. 29.
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 44, page 426.
  • H. S. Wall, Analytic Theory of Continued Fractions, Van Nostrand, New York, 1948, p. 356.

Crossrefs

Cf. A000522 (arrangements), A001620, A000262, A002720, A002793, A058006 (alternating factorial sums), A091725, A099285, A153229, A201203, A245780, A283743 (Ei(1)/e), A321942, A369883.

Programs

  • Magma
    SetDefaultRealField(RealField(100)); ExponentialIntegralE1(1)*Exp(1); // G. C. Greubel, Dec 04 2018
    
  • Mathematica
    RealDigits[N[-Exp[1]*ExpIntegralEi[-1], 105]][[1]]
    (* Second program: *)
    G = 1/Fold[Function[2*#2 - #2^2/#1], 2, Reverse[Range[10^4]]] // N[#, 105]&; RealDigits[G] // First (* Jean-François Alcover, Sep 19 2014 *)
  • PARI
    eint1(1)*exp(1) \\ Charles R Greathouse IV, Apr 23 2013
    
  • Sage
    numerical_approx(exp_integral_e(1,1)*exp(1), digits=100) # G. C. Greubel, Dec 04 2018

Formula

phi(1) = e*(Sum_{k>=1} (-1)^(k-1)/(k*k!) - Gamma) = 0.596347362323194... where Gamma is the Euler constant.
G = 0.596347... = 1/(1+1/(1+1/(1+2/(1+2/(1+3/(1+3/(1+4/(1+4/(1+5/(1+5/(1+6/(... - Philippe Deléham, Aug 14 2005
Equals A001113*A099285. - Johannes W. Meijer, Oct 16 2009
From Peter Bala, Oct 11 2012: (Start)
Stieltjes found the continued fraction representation G = 1/(2 - 1^2/(4 - 2^2/(6 - 3^2/(8 - ...)))). See [Wall, Chapter 18, (92.7) with a = 1]. The sequence of convergents to the continued fraction begins [1/2, 4/7, 20/34, 124/209, ...]. The numerators are in A002793 and the denominators in A002720.
Also, 1 - G has the continued fraction representation 1/(3 - 2/(5 - 6/(7 - ... -n*(n+1)/((2*n+3) - ...)))) with convergents beginning [1/3, 5/13, 29/73, 201/501, ...]. The numerators are in A201203 (unsigned) and the denominators are in A000262.
(End)
G = f(1) with f solution to the o.d.e. x^2*f'(x) + (x+1)*f(x)=1 such that f(0)=1. - Jean-François Alcover, May 28 2013
From Amiram Eldar, Aug 15 2020: (Start)
Equals Integral_{x=0..1} 1/(1-log(x)) dx.
Equals Integral_{x=1..oo} exp(1-x)/x dx.
Equals Integral_{x=0..oo} exp(-x)*log(x+1) dx.
Equals Integral_{x=0..oo} exp(-x)/(x+1) dx. (End)
From Gleb Koloskov, May 01 2021: (Start)
Equals Integral_{x=0..1} LambertW(e/x)-1 dx.
Equals Integral_{x=0..1} 1+1/LambertW(-1,-x/e) dx. (End)
Equals lim_{n->oo} A040027(n)/A000110(n+1). - Vaclav Kotesovec, Feb 22 2021
G = lim_{n->oo} A321942(n)/A000262(n). - Peter Bala, Mar 21 2022
Equals Sum_{n >= 1} 1/(n*L(n, -1)*L(n-1, -1)), where L(n, x) denotes the n-th Laguerre polynomial. This is the case x = 1 of the identity Integral_{t >= 0} exp(-t)/(x + t) dt = Sum_{n >= 1} 1/(n*L(n, -x)*L(n-1, -x)) valid for Re(x) > 0. - Peter Bala, Mar 21 2024
Equals lim_{n->oo} Sum_{k >= 0} (n/(n + 1))^k/(n + k). Cf. A099285. - Peter Bala, Jun 18 2024

Extensions

Additional references from Gerald McGarvey, Oct 10 2005
Link corrected by Johannes W. Meijer, Aug 01 2009

A005990 a(n) = (n-1)*(n+1)!/6.

Original entry on oeis.org

0, 1, 8, 60, 480, 4200, 40320, 423360, 4838400, 59875200, 798336000, 11416204800, 174356582400, 2833294464000, 48819843072000, 889218570240000, 17072996548608000, 344661117825024000, 7298706024529920000, 161787983543746560000
Offset: 1

Views

Author

Keywords

Comments

Coefficients of Gandhi polynomials.
a(n) = Sum_{pi in Symm(n)} Sum_{i=1..n} max(pi(i)-i,0), i.e., the total positive displacement of all letters in all permutations on n letters. - Franklin T. Adams-Watters, Oct 25 2006
a(n) is also the sum of the excedances of all permutations of [n]. An excedance of a permutation p of [n] is an i (1 <= i <= n-1) such that p(i) > i. Proof: i is an excedance if p(i) = i+1, i+2, ..., n (n-i possibilities), with the remaining values of p forming any permutation of [n]\{p(i)} in the positions [n]\{i} ((n-1)! possibilities). Summation of i(n-i)(n-1)! over i from 1 to n-1 completes the proof. Example: a(3)=8 because the permutations 123, 132, 213, 231, 312, 321 have excedances NONE, {2}, {1}, {1,2}, {1}, {1}, respectively. - Emeric Deutsch, Oct 26 2008
a(n) is also the number of doubledescents in all permutations of {1,2,...,n-1}. We say that i is a doubledescent of a permutation p if p(i) > p(i+1) > p(i+2). Example: a(3)=8 because each of the permutations 1432, 4312, 4213, 2431, 3214, 3421 has one doubledescent, the permutation 4321 has two doubledescents and the remaining 17 permutations of {1,2,3,4} have no doubledescents. - Emeric Deutsch, Jul 26 2009
Equals the second right hand column of A167568 divided by 2. - Johannes W. Meijer, Nov 12 2009
Half of sum of abs(p(i+1) - p(i)) over all permutations on n, e.g., 42531 = 2 + 3 + 2 + 2 = 9, and the total over all permutations on {1,2,3,4,5} is 960. - Jon Perry, May 24 2013
a(n) gives the number of non-occupied corners in tree-like tableaux of size n+1 (see Gao et al. link). - Michel Marcus, Nov 18 2015
a(n) is the number of sequences of n+2 balls colored with at most n colors such that exactly three balls are the same color as some other ball in the sequence. - Jeremy Dover, Sep 26 2017
a(n) is the number of triangles (3-cycles) in the (n+1)-alternating group graph. - Eric W. Weisstein, Jun 09 2019

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    [(n-1)*Factorial(n+1)/6: n in [1..25]]; // Vincenzo Librandi, Oct 11 2011
    
  • Maple
    [ seq((n-1)*(n+1)!/6,n=1..40) ];
    a:=n->sum(sum(sum(n!/6, j=1..n),k=-1..n),m=0..n): seq(a(n), n=0..19); # Zerinvary Lajos, May 11 2007
    seq(sum(mul(j,j=3..n), k=3..n)/3, n=2..21); # Zerinvary Lajos, Jun 01 2007
    restart: G(x):=x^3/(1-x)^2: f[0]:=G(x): for n from 1 to 21 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n]/3!,n=2..21); # Zerinvary Lajos, Apr 01 2009
  • Mathematica
    Table[Sum[n!/6, {i, 3, n}], {n, 2, 21}] (* Zerinvary Lajos, Jul 12 2009 *)
    Table[(n - 1) (n + 1)!/6, {n, 20}] (* Harvey P. Dale, Apr 07 2019 *)
    Table[(n - 1) Pochhammer[4, n - 2], {n, 20}] (* Eric W. Weisstein, Jun 09 2019 *)
    Table[(n - 1) Gamma[n + 2]/6, {n, 20}] (* Eric W. Weisstein, Jun 09 2019 *)
    Range[0, 20]! CoefficientList[Series[x/(1 - x)^4, {x, 0, 20}], x] (* Eric W. Weisstein, Jun 09 2019 *)
  • PARI
    a(n)=(n-1)*(n+1)!/6 \\ Charles R Greathouse IV, May 24 2013

Formula

a(n) = A090672(n)/2.
a(n) = A052571(n+2)/6. - Zerinvary Lajos, May 11 2007
a(n) = Sum_{m=0..n} Sum_{k=-1..n} Sum_{j=1..n} n!/6, n >= 0. - Zerinvary Lajos, May 11 2007
If we define f(n,i,x) = Sum_{k=i..n} (Sum_{j=i..k} binomial(k,j)*Stirling1(n,k)*Stirling2(j,i)*x^(k-j)) then a(n+1) = (-1)^(n-1)*f(n,1,-4), (n >= 1). - Milan Janjic, Mar 01 2009
E.g.f.: (-1+3*x)/(3!*(1-x)^3), a(0) = -1/3!. Such e.g.f. computations resulted from e-mail exchange with Gary Detlefs. - Wolfdieter Lang, May 27 2010
a(n) = ((n+3)!/2) * Sum_{j=i..k} (k+1)!/(k+3)!, with offset 0. - Gary Detlefs, Aug 05 2010
a(n) = (n+2)!*Sum_{k=1..n-1} 1/((2*k+4)*(k+3)). - Gary Detlefs, Oct 09 2011
a(n) = (n+2)!*(1 + 3*(H(n+1) - H(n+2)))/6, where H(n) is the n-th harmonic number. - Gary Detlefs, Oct 09 2011
With offset = 0, e.g.f.: x/(1-x)^4. - Geoffrey Critzer, Aug 30 2013
From Amiram Eldar, May 06 2022: (Start)
Sum_{n>=2} 1/a(n) = 3*(Ei(1) - gamma) - 6*e + 27/2, where Ei(1) = A091725, gamma = A001620, and e = A001113.
Sum_{n>=2} (-1)^n/a(n) = 3*(gamma - Ei(-1)) - 3/2, where Ei(-1) = -A099285. (End)

Extensions

Better definition from Robert Newstedt

A062119 a(n) = n! * (n-1).

Original entry on oeis.org

0, 2, 12, 72, 480, 3600, 30240, 282240, 2903040, 32659200, 399168000, 5269017600, 74724249600, 1133317785600, 18307441152000, 313841848320000, 5690998849536000, 108840352997376000, 2189611807358976000, 46225138155356160000, 1021818843434188800000
Offset: 1

Views

Author

Olivier Gérard, Jun 13 2001

Keywords

Comments

For n > 0, a(n) = number of permutations of length n+1 that have 2 predetermined elements nonadjacent; e.g., for n=2, the permutations with, say, 1 and 2 nonadjacent are 132 and 231, therefore a(2)=2. - Jon Perry, Jun 08 2003
Number of multiplications performed when computing the determinant of an n X n matrix by definition. - Mats Granvik, Sep 12 2008
Sum of the length of all cycles (excluding fixed points) in all permutations of [n]. - Olivier Gérard, Oct 23 2012
Number of permutations of n distinct objects (ABC...) 1 (one) times >>("-", A, AB, ABC, ABCD, ABCDE, ..., ABCDEFGHIJK, infinity) and one after the other to resemble motif: A (1) AB (1-1), AAB (2-1), AAAB (3-1), AAAAB (4-1), AAAAAB (5-1), AAAAAAB (6-1), AAAAAAAB (7-1), AAAAAAAAB (8-1) etc.,>> "1(one) fixed point". Example:motif: AAAB (or BBBA) 12 * one (1) fixed point etc. Let: AAAB ................ 'A'BCD 1. 'A'BDC 2. 'A'CBD 3. ACDB 'A'DBC 4. 'A'DCB B'A'CD 5. B'A'DC 6. BCAD 7. BCDA BD'A'C 8. BDCA C'A'BD 9. C'A'DB CB'A'D 10. CBDA CDAB CDBA D'A'BC 11. DACB DB'A'C 12. DBCA DCAB DCBA. - Zerinvary Lajos, Nov 27 2009 (does anybody understand what this is supposed to say? - Joerg Arndt, Jan 10 2015)
a(n) is the number of ways to arrange n books on two bookshelves so that each shelf receives at least one book. - Geoffrey Critzer, Feb 21 2010
a(n) = number whose factorial base representation (A007623) begins with digit {n-1} and is followed by n-1 zeros. Viewed in that base, this sequence looks like this: 0, 10, 200, 3000, 40000, 500000, 6000000, 70000000, 800000000, 9000000000, A0000000000, B00000000000, ... (where "digits" A and B stand for placeholder values 10 and 11 respectively). - Antti Karttunen, May 07 2015

Crossrefs

Column 2 of A257503 (apart from initial zero. Equally, row 2 of A257505).
Cf. A001286 (same sequence divided by 2).
Cf. A001563. - Zerinvary Lajos, Aug 27 2008
Cf. sequences with formula (n + k)*n! listed in A282466.

Programs

Formula

a(n) = n! * (n-1).
E.g.f.: x^2/(1-x)^2. - Geoffrey Critzer, Feb 21 2010
a(n) = 2 * A001286(n).
a(n) = A001563(n) - A000142(n). - Antti Karttunen, May 07 2015, hinted by crossref left by Lajos.
From Amiram Eldar, Jul 11 2020: (Start)
Sum_{n>=2} 1/a(n) = Ei(1) + 2 - e - gamma = A091725 + 2 - A001113 - A001620.
Sum_{n>=2} (-1)^n/a(n) = gamma - Ei(-1) - 1/e = A001620 + A099285 - A068985. (End)

Extensions

Last term a(19) corrected by Harry J. Smith, Aug 02 2009

A001754 Lah numbers: a(n) = n!*binomial(n-1,2)/6.

Original entry on oeis.org

0, 0, 1, 12, 120, 1200, 12600, 141120, 1693440, 21772800, 299376000, 4390848000, 68497228800, 1133317785600, 19833061248000, 366148823040000, 7113748561920000, 145120470663168000, 3101950060425216000, 69337707233034240000, 1617879835437465600000
Offset: 1

Views

Author

Keywords

Comments

a(n+1) = Sum_{pi in Symm(n)} Sum_{i=1..n} max(pi(i)-i,0)^2, i.e., the sum of the squares of the positive displacement of all letters in all permutations on n letters. - Franklin T. Adams-Watters, Oct 25 2006

References

  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 156.
  • John Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 44.
  • 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).

Crossrefs

Column 3 of A008297.
Column m=3 of unsigned triangle A111596.

Programs

  • Magma
    [Factorial(n)*Binomial(n-1, 2)/6: n in [1..25]]; // Vincenzo Librandi, Oct 11 2011
    
  • Maple
    [seq(n!*binomial(n-1,2)/6, n=1..40)];
  • Mathematica
    Table[(n-2)*(n-1)*n!/12, {n, 21}] (* Arkadiusz Wesolowski, Nov 26 2012 *)
    With[{nn=30},CoefficientList[Series[(x/(1-x))^3/6,{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Oct 04 2017 *)
  • Sage
    [factorial(n-1)*binomial(n,3)/2 for n in (1..30)] # G. C. Greubel, May 10 2021

Formula

E.g.f.: ((x/(1-x))^3)/3!.
If we define f(n,i,x) = Sum_{k=i..n} (Sum_{j=i..k} binomial(k,j) * Stirling1(n,k) * Stirling2(j,i)*x^(k-j)) then a(n+1) = (-1)^n*f(n,2,-4), n >= 2. - Milan Janjic, Mar 01 2009
a(n) = Sum_{k>=1} k * A260665(n,k). - Alois P. Heinz, Nov 14 2015
D-finite with recurrence (-n+5)*a(n) + (n-2)*(n-3)*a(n-1) = 0, n >= 4. - R. J. Mathar, Jan 06 2021
From Amiram Eldar, May 02 2022: (Start)
Sum_{n>=3} 1/a(n) = 6*(gamma - Ei(1)) + 9, where gamma = A001620 and Ei(1) = A091725.
Sum_{n>=3} (-1)^(n+1)/a(n) = 18*(gamma - Ei(-1)) - 12/e - 9, where Ei(-1) = -A099285 and e = A001113. (End)

A074143 a(1) = 1; a(n) = n * Sum_{k=1..n-1} a(k).

Original entry on oeis.org

1, 2, 9, 48, 300, 2160, 17640, 161280, 1632960, 18144000, 219542400, 2874009600, 40475635200, 610248038400, 9807557760000, 167382319104000, 3023343138816000, 57621363351552000, 1155628453883904000, 24329020081766400000, 536454892802949120000
Offset: 1

Views

Author

Amarnath Murthy, Aug 28 2002

Keywords

Comments

a(n) is also the number of elements of the alternating semigroup (A^c_n) for F(n, p) if p = n - 1 (cf. A001710). - Bakare Gatta Naimat, Jan 15 2016

Crossrefs

Programs

Formula

a(n) = n^2 * a(n-1)/(n-1) for n > 2.
a(n) = n*ceiling(n!/2) = n*A001710(n) = ceiling(A001563(n)/2). - Henry Bottomley, Nov 27 2002
a(n) = ((n+1)!-n!)/2 for n > 1. - Vladimir Joseph Stephan Orlovsky, Apr 03 2011
G.f.: (U(0) + x)/(2*x) where U(k) = 1 - 1/(k+1 - x*(k+1)^2*(k+2)/(x*(k+1)*(k+2) - 1/U(k+1))); (continued fraction). - Sergei N. Gladkovskii, Sep 27 2012
G.f.: 1/2 + Q(0), where Q(k)= 1 - 1/(k+2 - x*(k+2)^2*(k+3)/(x*(k+2)*(k+3)-1/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 19 2013
a(n) = Sum_{j = 0..n} (-1)^(n-j)*binomial(n, j)*(j)^(n+1) / (n+1), n > 1, a(1) = 1. - Vladimir Kruchinin, Jun 01 2013
a(n) = numerator(n!/2*n). - Vincenzo Librandi, Apr 15 2014
a(n) is F(n;p) = n^2(n-1)!/2 if p = n-1 in A^c_n. For instance for n=4 and p=n-1: F(4; 4-1)= 4^2(4-1)!/2 = 16*6/2 = 48. - Bakare Gatta Naimat, Nov 18 2015
From Seiichi Manyama, Apr 27 2025: (Start)
E.g.f.: x/2 * (1 + 1/(1-x)^2).
a(n) = (n+2) * a(n-1) - (n-1) * a(n-2) for n > 3. (End)
From Amiram Eldar, May 04 2025: (Start)
Sum_{n>=1} 1/a(n) = 2*ExpIntegralEi(1) - 2*gamma - 1 = 2*A091725 - 2*A001620 - 1.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*gamma - 1 - 2*ExpIntegralEi(-1) = 2*A001620 - 1 + 2*A099285. (End)

Extensions

More terms from Henry Bottomley, Nov 27 2002

A005461 Number of simplices in barycentric subdivision of n-simplex.

Original entry on oeis.org

1, 15, 180, 2100, 25200, 317520, 4233600, 59875200, 898128000, 14270256000, 239740300800, 4249941696000, 79332244992000, 1556132497920000, 32011868528640000, 689322235650048000, 15509750302126080000, 364022962973429760000, 8898339094906060800000
Offset: 1

Views

Author

Keywords

Examples

			G.f. = x + 15*x^2 + 180*x^3 + 2100*x^4 + 25200*x^5 + 317520*x^6 + ...
		

References

  • R. Austin, R. K. Guy, and R. Nowakowski, unpublished notes, circa 1987.
  • R. K. Guy, personal communication.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    [Factorial(n-1)*StirlingSecond(n+3,n): n in [1..35]]; // G. C. Greubel, Nov 23 2022
  • Maple
    a:=n->sum((n-j)*n!/4!, j=3..n): seq(a(n), n=4..17); # Zerinvary Lajos, Apr 29 2007
  • Mathematica
    Table[(n(n+1)(n+3)!)/48,{n,20}] (* Harvey P. Dale, Mar 14 2012 *)
    a[ n_] := If[ n < 0, 0, n (n + 1) (n + 3)! / 48]; (* Michael Somos, May 27 2014 *)
  • Sage
    [factorial(m+1)*binomial(m-1,2)/24 for m in range(3, 19)] # Zerinvary Lajos, Jul 05 2008
    
  • Sage
    [binomial(n,4)*factorial (n-2)/2 for n in range(4, 18)] #  Zerinvary Lajos, Jul 07 2009
    

Formula

a(n) = n*(n + 1)*(n + 3)!/48.
Essentially Stirling numbers of second kind - see A028246.
If we define f(n,i,x) = Sum_{k=i..n} Sum_{j=i..k} binomial(k,j)*Stirling1(n,k)*Stirling2(j,i)*x^(k-j) then a(n-3) = (-1)^n*f(n,4,-3), (n>=4). - Milan Janjic, Mar 01 2009
E.g.f.: t*(3*t + 2)/(2*(t - 1)^6). - Ran Pan, Jul 10 2016
a(n) ~ sqrt(Pi/2)*exp(-n)*n^(n+1/2)*(n^5/24 + 85*n^4/288 + 5065*n^3/6912 + 955841*n^2/1244160 + 3710929*n/11943936). - Ilya Gutkovskiy, Jul 10 2016
From Amiram Eldar, May 06 2022: (Start)
Sum_{n>=1} 1/a(n) = 16*(e + gamma - Ei(1)) - 64/3, where e = A001113, gamma = A001620, and Ei(1) = A091725.
Sum_{n>=1} (-1)^(n+1)/a(n) = 32*(gamma - Ei(-1)) - 16/e - 56/3, where Ei(-1) = -A099285. (End)
a(n) = (n-1)! * Stirling2(n+3, n). - G. C. Greubel, Nov 23 2022

Extensions

More terms from Harvey P. Dale, Mar 14 2012
Showing 1-10 of 47 results. Next