A002775 -id:A002775 - cp's OEIS frontend

A001563 a(n) = n*n! = (n+1)! - n!.

0, 1, 4, 18, 96, 600, 4320, 35280, 322560, 3265920, 36288000, 439084800, 5748019200, 80951270400, 1220496076800, 19615115520000, 334764638208000, 6046686277632000, 115242726703104000, 2311256907767808000, 48658040163532800000, 1072909785605898240000

Offset: 0

N. J. A. Sloane

A similar sequence, with the initial 0 replaced by 1, namely A094258, is defined by the recurrence a(2) = 1, a(n) = a(n-1)*(n-1)^2/(n-2). - Andrey Ryshevich (ryshevich(AT)notes.idlab.net), May 21 2002
Denominators in power series expansion of E_1(x) + gamma + log(x), x > 0. - Michael Somos, Dec 11 2002
If all the permutations of any length k are arranged in lexicographic order, the n-th term in this sequence (n <= k) gives the index of the permutation that rotates the last n elements one position to the right. E.g., there are 24 permutations of 4 items. In lexicographic order they are (0,1,2,3), (0,1,3,2), (0,2,1,3), ... (3,2,0,1), (3,2,1,0). Permutation 0 is (0,1,2,3), which rotates the last 1 element, i.e., it makes no change. Permutation 1 is (0,1,3,2), which rotates the last 2 elements. Permutation 4 is (0,3,1,2), which rotates the last 3 elements. Permutation 18 is (3,0,1,2), which rotates the last 4 elements. The same numbers work for permutations of any length. - Henry H. Rich (glasss(AT)bellsouth.net), Sep 27 2003
Stirling transform of a(n+1)=[4,18,96,600,...] is A083140(n+1)=[4,22,154,...]. - Michael Somos, Mar 04 2004
From Michael Somos, Apr 27 2012: (Start)
Stirling transform of a(n)=[1,4,18,96,...] is A069321(n)=[1,5,31,233,...].
Partial sums of a(n)=[0,1,4,18,...] is A033312(n+1)=[0,1,5,23,...].
Binomial transform of A000166(n+1)=[0,1,2,9,...] is a(n)=[0,1,4,18,...].
Binomial transform of A000255(n+1)=[1,3,11,53,...] is a(n+1)=[1,4,18,96,...].
Binomial transform of a(n)=[0,1,4,18,...] is A093964(n)=[0,1,6,33,...].
Partial sums of A001564(n)=[1,3,4,14,...] is a(n+1)=[1,4,18,96,...].
(End)
Number of small descents in all permutations of [n+1]. A small descent in a permutation (x_1,x_2,...,x_n) is a position i such that x_i - x_(i+1) =1. Example: a(2)=4 because there are 4 small descents in the permutations 123, 13\2, 2\13, 231, 312, 3\2\1 of {1,2,3} (shown by \). a(n)=Sum_{k=0..n-1}k*A123513(n,k). - Emeric Deutsch, Oct 02 2006
Equivalently, in the notation of David, Kendall and Barton, p. 263, this is the total number of consecutive ascending pairs in all permutations on n+1 letters (cf. A010027). - N. J. A. Sloane, Apr 12 2014
a(n-1) is the number of permutations of n in which n is not fixed; equivalently, the number of permutations of the positive integers in which n is the largest element that is not fixed. - Franklin T. Adams-Watters, Nov 29 2006
Number of factors in a determinant when writing down all multiplication permutations. - Mats Granvik, Sep 12 2008
a(n) is also the sum of the positions of the left-to-right maxima in all permutations of [n]. Example: a(3)=18 because the positions of the left-to-right maxima in the permutations 123,132,213,231,312 and 321 of [3] are 123, 12, 13, 12, 1 and 1, respectively and 1+2+3+1+2+1+3+1+2+1+1=18. - Emeric Deutsch, Sep 21 2008
Equals eigensequence of triangle A002024 ("n appears n times"). - Gary W. Adamson, Dec 29 2008
Preface the series with another 1: (1, 1, 4, 18, ...); then the next term = dot product of the latter with "n occurs n times". Example: 96 = (1, 1, 4, 8) dot (4, 4, 4, 4) = (4 + 4 + 16 + 72). - Gary W. Adamson, Apr 17 2009
Row lengths of the triangle in A030298. - Reinhard Zumkeller, Mar 29 2012
a(n) is also the number of minimum (n-)distinguishing labelings of the star graph S_{n+1} on n+1 nodes. - Eric W. Weisstein, Oct 14 2014
When the numbers denote finite permutations (as row numbers of A055089) these are the circular shifts to the right, i.e., a(n) is the permutation with the cycle notation (0 1 ... n-1 n). Compare array A051683 for circular shifts to the right in a broader sense. Compare sequence A007489 for circular shifts to the left. - Tilman Piesk, Apr 29 2017
a(n-1) is the number of permutations on n elements with no cycles of length n. - Dennis P. Walsh, Oct 02 2017
The number of pandigital numbers in base n+1, such that each digit appears exactly once. For example, there are a(9) = 9*9! = 3265920 pandigital numbers in base 10 (A050278). - Amiram Eldar, Apr 13 2020

			E_1(x) + gamma + log(x) = x/1 - x^2/4 + x^3/18 - x^4/96 + ..., x > 0. - _Michael Somos_, Dec 11 2002
G.f. = x + 4*x^2 + 18*x^3 + 96*x^4 + 600*x^5 + 4320*x^6 + 35280*x^7 + 322560*x^8 + ...

A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 218.
J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p. 336.
F. N. David, M. G. Kendall, and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 263.
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 37, equation 37:6:1 at page 354.

T. D. Noe, Table of n, a(n) for n = 0..100
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122. [Annotated scanned copy]
Loïc Foissy, The antipode of of [sic] a Com-PreLie Hopf algebra, arXiv:2406.01120 [math.CO], 2024. See p. 9.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 30.
Milan Janjic, Enumerative Formulas for Some Functions on Finite Sets.
L. B. W. Jolley, Summation of Series, Dover, 1961.
I. Kortchemski, Asymptotic behavior of permutation records, arXiv: 0804.0446v2 [math.CO], 18 May 2008.
C. Lanczos, Applied Analysis (Annotated scans of selected pages).
Rezsö L. Lovas and István Mezö, On an exotic topology of the integers, arXiv:1008.0713 [math.GN], 2010. See p. 4.
Daniel J. Mundfrom, A problem in permutations: the game of `Mousetrap', European J. Combin. 15 (1994), no. 6, 555-560.
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014.
J. Ser, Les Calculs Formels des Séries de Factorielles, Gauthier-Villars, Paris, 1933 [Local copy].
J. Ser, Les Calculs Formels des Séries de Factorielles (Annotated scans of some selected pages).
A. van Heemert, Cyclic permutations with sequences and related problems, J. Reine Angew. Math., 198 (1957), 56-72.
Dennis P. Walsh, The number of permutations with no k-cycles.
Eric Weisstein's World of Mathematics, Distinguishing Number.
Eric Weisstein's World of Mathematics, Exponential Integral.

Cf. A000142, A002024, A047920, A047922, A053495, A055089, A123513, A143946, A229837, A239069.
Cf. A163931 (E(x,m,n)), A002775 (n^2*n!), A091363 (n^3*n!), A091364 (n^4*n!).
Cf. sequences with formula (n + k)*n! listed in A282466.
Row sums of A185105, A322383, A322384, A094485.

List([0..20], n-> n*Factorial(n) ); # G. C. Greubel, Dec 30 2019

a001563 n = a001563_list !! n
a001563_list = zipWith (-) (tail a000142_list) a000142_list
-- Reinhard Zumkeller, Aug 05 2013

[Factorial(n+1)-Factorial(n): n in [0..20]]; // Vincenzo Librandi, Aug 08 2014

Maple
```
A001563 := n->n*n!;
```

Table[n!n,{n,0,25}] (* Harvey P. Dale, Oct 03 2011 *)

{a(n) = if( n<0, 0, n * n!)} /* Michael Somos, Dec 11 2002 */

[n*factorial(n) for n in (0..20)] # G. C. Greubel, Dec 30 2019

From Michael Somos, Dec 11 2002: (Start)
E.g.f.: x / (1 - x)^2.
a(n) = -A021009(n, 1), n >= 0. (End)
The coefficient of y^(n-1) in expansion of (y+n!)^n, n >= 1, gives the sequence 1, 4, 18, 96, 600, 4320, 35280, ... - Artur Jasinski, Oct 22 2007
Integral representation as n-th moment of a function on a positive half-axis: a(n) = Integral_{x=0..oo} x^n*(x*(x-1)*exp(-x)) dx, for n>=0. This representation may not be unique. - Karol A. Penson, Sep 27 2001
a(0)=0, a(n) = n*a(n-1) + n!. - Benoit Cloitre, Feb 16 2003
a(0) = 0, a(n) = (n - 1) * (1 + Sum_{i=1..n-1} a(i)) for i > 0. - Gerald McGarvey, Jun 11 2004
Arises in the denominators of the following identities: Sum_{n>=1} 1/(n*(n+1)*(n+2)) = 1/4, Sum_{n>=1} 1/(n*(n+1)*(n+2)*(n+3)) = 1/18, Sum_{n>=1} 1/(n*(n+1)*(n+2)*(n+3)*(n+4)) = 1/96, etc. The general expression is Sum_{n>=k} 1/C(n, k) = k/(k-1). - Dick Boland, Jun 06 2005 [And the general expression implies that Sum_{n>=1} 1/(n*(n+1)*...*(n+k-1)) = (Sum_{n>=k} 1/C(n, k))/k! = 1/((k-1)*(k-1)!) = 1/a(k-1), k >= 2. - Jianing Song, May 07 2023]
a(n) = Sum_{m=2..n+1} |Stirling1(n+1, m)|, n >= 1 and a(0):=0, where Stirling1(n, m) = A048994(n, m), n >= m = 0.
a(n) = 1/(Sum_{k>=0} k!/(n+k+1)!), n > 0. - Vladeta Jovovic, Sep 13 2006
a(n) = Sum_{k=1..n(n+1)/2} k*A143946(n,k). - Emeric Deutsch, Sep 21 2008
The reciprocals of a(n) are the lead coefficients in the factored form of the polynomials obtained by summing the binomial coefficients with a fixed lower term up to n as the upper term, divided by the term index, for n >= 1: Sum_{k = i..n} C(k, i)/k = (1/a(n))*n*(n-1)*..*(n-i+1). The first few such polynomials are Sum_{k = 1..n} C(k, 1)/k = (1/1)*n, Sum_{k = 2..n} C(k, 2)/k = (1/4)*n*(n-1), Sum_{k = 3..n} C(k, 3)/k = (1/18)*n*(n-1)*(n-2), Sum_{k = 4..n} C(k, 4)/k = (1/96)*n*(n-1)*(n-2)*(n-3), etc. - Peter Breznay (breznayp(AT)uwgb.edu), Sep 28 2008
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-1)*f(n,1,-2), (n >= 1). - Milan Janjic, Mar 01 2009
Sum_{n>=1} (-1)^(n+1)/a(n) =  0.796599599... [Jolley eq. 289]
G.f.: 2*x*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
G.f.: W(0)*(1-sqrt(x)) - 1, where W(k) = 1 + sqrt(x)/( 1 - sqrt(x)*(k+2)/(sqrt(x)*(k+2) + 1/W(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 18 2013
G.f.: T(0)/x - 1/x, where T(k) = 1 - x^2*(k+1)^2/( x^2*(k+1)^2 - (1-x-2*x*k)*(1-3*x-2*x*k)/T(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Oct 17 2013
G.f.: Q(0)*(1-x)/x - 1/x, where Q(k) = 1 - x*(k+1)/( x*(k+1) - 1/(1 - x*(k+1)/( x*(k+1) - 1/Q(k+1) ))); (continued fraction). - Sergei N. Gladkovskii, Oct 22 2013
D-finite with recurrence: a(n) +(-n-2)*a(n-1) +(n-1)*a(n-2)=0. - R. J. Mathar, Jan 14 2020
a(n) = (-1)^(n+1)*(n+1)*Sum_{k=1..n} A094485(n,k)*Bernoulli(k). The inverse of the Worpitzky representation of the Bernoulli numbers. - Peter Luschny, May 28 2020
From Amiram Eldar, Aug 04 2020: (Start)
Sum_{n>=1} 1/a(n) = Ei(1) - gamma = A229837.
Sum_{n>=1} (-1)^(n+1)/a(n) = gamma - Ei(-1) = A239069. (End)
a(n) = Gamma(n)*A000290(n) for n > 0. - Jacob Szlachetka, Jan 01 2022

A163931 Decimal expansion of the higher-order exponential integral E(x, m=2, n=1) at x=1.

Original entry on oeis.org

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

Offset: 0

Johannes W. Meijer and Nico Baken, Aug 13 2009, Aug 17 2009

We define the higher-order exponential integrals by E(x,m,n) = x^(n-1)*Integral_{t=x..infinity} E(t,m-1,n)/t^n for m >= 1 and n >= 1 with E(x,m=0,n) = exp(-x), see Meijer and Baken.
The properties of the E(x,m,n) are analogous to those of the well-known exponential integrals E(x,m=1,n), see Abramowitz and Stegun and the formulas.
The series expansions of the higher-order exponential integrals are dominated by the constants alpha(k,n), see A163927, and gamma(k,n) = G(k,n), see A090998.
For information about the asymptotic expansion of the E(x,m,n) see A163932.
Values of E(x,m,n) can be evaluated with the Maple program.

			E(1,2,1) = 0.09784319721667017932553778904528008276958226953026576557442124245....

G. C. Greubel, Table of n, a(n) for n = 0..5000
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972, Chapter 5, pp. 227-251.
J. W. Meijer and N. H. G. Baken, The Exponential Integral Distribution, Statistics and Probability Letters, Volume 5, No.3, April 1987. pp 209-211.
M. S. Milgram, The generalized integro-exponential function, Math. of Computation, Vol. 44, pp. 443-458, 1985.
Eric Weisstein's World of Mathematics, The Exponential Integral.

Cf. A163927 (alpha(k,n)), A090998 (gamma(k,n) = G(k,n)), A163932.
Cf. A068985 (E(x=1,m=0,n) = exp(-1)) and A099285 (E(x=1,m=1,n=1)).
Cf. A001563 (n*n!), A002775 (n^2*n!), A091363 (n^3*n!) and A091364 (n^4*n!).

E:= proc(x,m,n) local nmax, kmax, EI, k1, k2, n1, n2; option remember: nmax:=20; kmax:=20; k1:=0: for n1 from 0 to nmax do alpha(k1,n1):=1 od: for k1 from 1 to kmax do for n1 from 1 to nmax do alpha(k1,n1) := (1/k1)*sum(sum(p^(-2*(k1-i1)),p=0..n1-1)*alpha(i1, n1),i1=0..k1-1) od; od: for n2 from 0 to kmax do G(0,n2):=1 od: for n2 from 1 to nmax do for k2 from 1 to kmax do G(k2,n2):=(1/k2)*(((gamma-sum(p^(-1),p=1..n2-1))*G(k2-1,n2)+ sum((Zeta(k2-i2)-sum(p^(-(k2-i2)), p=1..n2-1))*G(i2,n2),i2=0..k2-2))) od; od: EI:= evalf((-1)^m*((-x)^(n-1)/(n-1)!*sum(alpha(kz,n)*(G(m-2*kz,n)+sum(G(m-2*kz-i,n)*ln(x)^i/i!,i=1..m-2*kz)), kz=0..floor(m/2)) + sum((-x)^kx/((kx-n+1)^m*kx!),kx=0..n-2) + sum((-x)^ky/((ky-n+1)^m*ky!),ky=n..infinity))); return(EI): end:

Join[{0}, RealDigits[ N[ EulerGamma^2/2 + Pi^2/12 - HypergeometricPFQ[{1, 1, 1}, {2, 2, 2}, -1], 104]][[1]]] (* Jean-François Alcover, Nov 07 2012, from 1st formula *)

t=1; Euler^2/2 + Pi^2/12 + sumalt(k=1, t*=k; (-1)^k/(k^2*t)) \\ Charles R Greathouse IV, Nov 07 2016

E(x=1,m=2,n=1) = gamma^2/2 + Pi^2/12 + Sum_{k>=1} ((-1)^k/(k^2*k!)).
E(x=0,n,m) = (1/(n-1))^m for n >= 2.
Integral_{t=0..x} E(t,m,n) = 1/n^m - E(x,n,n+1).
dE(x,m,n+1)/dx = - E(x,m,n).
E(x,m,n+1) = (1/n)*(E(x,m-1,n+1) - x*E(x,m,n)).
E(x,m,n) = (-1)^m * ((-x)^(n-1)/(n-1)!) * Sum_{kz=0..floor(m/2)}(alpha (kz, n)*G(m-2*kz, n)) + (-1) ^m * ((-x)^(n-1)/(n-1)!) * Sum_{kz=0..floor(m/2)}(Sum_{i=1..m-2*kz}(alpha (kz, n) *G(m-2*kz-i, n)*log(x)^i/i!)) + (-1)^m * Sum_{ kx=0..n-2}((-x)^kx/((kx-n+1)^m*kx!) + (-1)^m * Sum_{ky>=n}((-x)^ky /(( ky-n+1)^m*ky!)).

A047920 Triangular array formed from successive differences of factorial numbers.

Original entry on oeis.org

1, 1, 0, 2, 1, 1, 6, 4, 3, 2, 24, 18, 14, 11, 9, 120, 96, 78, 64, 53, 44, 720, 600, 504, 426, 362, 309, 265, 5040, 4320, 3720, 3216, 2790, 2428, 2119, 1854, 40320, 35280, 30960, 27240, 24024, 21234, 18806, 16687, 14833, 362880, 322560

Offset: 0

N. J. A. Sloane

Number of permutations of 1,2,...,k,n+1,n+2,...,2n-k that have no agreements with 1,...,n. For example, consider 1234 and 1256, then n=4 and k=2, so T(4,2)=14. Compare A000255 for the case k=1. - Jon Perry, Jan 23 2004
From Emeric Deutsch, Apr 21 2009: (Start)
T(n-1,k-1) is the number of non-derangements of {1,2,...,n} having smallest fixed point equal to k. Example: T(3,1)=4 because we have 4213, 4231, 3214, and 3241 (the permutations of {1,2,3,4} having smallest fixed equal to 2).
Row sums give the number of non-derangement permutations of {1,2,...,n} (A002467).
Mirror image of A068106.
Closely related to A134830, where each row has an extra term (see the Charalambides reference).
(End)
T(n,k) is the number of permutations of {1..n} that don't fix the points 1..k. - Robert FERREOL, Aug 04 2016

			Triangle begins:
    1;
    1,  0;
    2,  1,  1;
    6,  4,  3,  2;
   24, 18, 14, 11,  9;
  120, 96, 78, 64, 53, 44;
  ...
The left-hand column is the factorial numbers (A000142). The other numbers in the row are calculated by subtracting the numbers in the previous row. For example, row 4 is 6, 4, 3, 2, so row 5 is 4! = 24, 24 - 6 = 18, 18 - 4 = 14, 14 - 3 = 11, 11 - 2 = 9. - _Michael B. Porter_, Aug 05 2016

Ch. A. Charalambides, Enumerative Combinatorics, Chapman & Hall/CRC, Boca Raton, Florida, 2002, p. 176, Table 5.3. [From Emeric Deutsch, Apr 21 2009]

Reinhard Zumkeller, Rows n = 0..150 of triangle, flattened
E. Deutsch and S. Elizalde, The largest and the smallest fixed points of permutations, arXiv:0904.2792 [math.CO], 2009.
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122. [Annotated scanned copy]
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122.
Ira M. Gessel, Symmetric inclusion-exclusion, Séminaire Lotharingien de Combinatoire, B54b (2005).
Peter Kagey, Ranking and Unranking Restricted Permutations, arXiv:2210.17021 [math.CO], 2022.
Index entries for sequences related to factorial numbers

Columns give A000142, A001563, A001564, etc. Cf. A047922.
See A068106 for another version of this triangle.
Orthogonal columns: A000166, A000255, A055790. Main diagonal A033815.
Cf. A002467, A068106, A134830. - Emeric Deutsch, Apr 21 2009
Cf. A155521.
T(n+2,n) = 2*A000153(n+1). T(n+3,n) = 6*A000261(n+2). T(n+4,n) = 24*A001909(n+3). T(n+5, n) = 120*A001910(n+4). T(n+6,n) = 720*A176732(n).
T(n+7,n) = 5040*A176733(n) - Richard R. Forberg, Dec 29 2013

a047920 n k = a047920_tabl !! n !! k
a047920_row n = a047920_tabl !! n
a047920_tabl = map fst $ iterate e ([1], 1) where
   e (row, n) = (scanl (-) (n * head row) row, n + 1)
-- Reinhard Zumkeller, Mar 05 2012

d[0] := 1: for n to 15 do d[n] := n*d[n-1]+(-1)^n end do: T := proc (n, k) if k <= n then sum(binomial(n-k, j)*d[n-j], j = 0 .. n-k) else 0 end if end proc: for n from 0 to 9 do seq(T(n, k), k = 0 .. n) end do; # yields sequence in triangular form - Emeric Deutsch, Jul 17 2009
# second Maple program:
T:= proc(n, k) option remember;
     `if`(k=0, n!, T(n, k-1)-T(n-1, k-1))
    end:
seq(seq(T(n, k), k=0..n), n=0..12);  # Alois P. Heinz, Sep 01 2021

t[n_, k_] = Sum[(-1)^j*Binomial[k, j]*(n-j)!, {j, 0, n}]; Flatten[Table[t[n, k], {n, 0, 9}, {k, 0, n}]][[1 ;; 47]] (* Jean-François Alcover, May 17 2011, after Philippe Deléham *)
T[n_, k_] := n! Hypergeometric1F1[-k, -n, -1];
Table[T[n, k], {n, 0, 7}, {k, 0, n}] // Flatten  (* Peter Luschny, Jul 28 2024 *)

row(n) = vector(n+1, k, k--; sum(j=0, n, (-1)^j * binomial(k, j)*(n-j)!)); \\ Michel Marcus, Sep 04 2021

t(n, k) = t(n, k-1) - t(n-1, k-1) = t(n, k+1) - t(n-1, k) = n*t(n-1, k) + k*t(n-2, k-1) = (n-1)*t(n-1, k-1) + (k-1)*t(n-2, k-2) = A060475(n, k)*(n-k)!. - Henry Bottomley, Mar 16 2001
T(n, k) = Sum_{j>=0} (-1)^j * binomial(k, j)*(n-j)!. - Philippe Deléham, May 29 2005
T(n,k) = Sum_{j=0..n-k} d(n-j)*binomial(n-k,j), where d(i)=A000166(i) are the derangement numbers. - Emeric Deutsch, Jul 17 2009
Sum_{k=0..n} (k+1)*T(n,k) = A155521(n+1). - Emeric Deutsch, Jul 18 2009
T(n, k) = n!*hypergeom([-k], [-n], -1). - Peter Luschny, Jul 28 2024

A000271 Sums of ménage numbers.

Original entry on oeis.org

1, 0, 0, 1, 3, 16, 96, 675, 5413, 48800, 488592, 5379333, 64595975, 840192288, 11767626752, 176574062535, 2825965531593, 48052401132800, 865108807357216, 16439727718351881, 328839946389605643, 6906458590966507696

Offset: 0

N. J. A. Sloane

Permanent of the (0,1)-matrix having (i,j)-th entry equal to 0 iff this is on the diagonal or the first upper-diagonal. - Simone Severini, Oct 14 2004
Equivalently, number of permutations p of {1,2,...,n} such that p(i)-i not in {0,1}. - Andrew Howroyd, Sep 19 2017
From Vladimir Shevelev, Jun 21 2015: (Start)
Let 2*n!*V(n)=A137886(n) be the number of ways of seating n married couples at 2*n chairs arranged side-by-side in a straight line, men and women in alternate positions, so that no husband is next to his wife.
It is known [Riordan, Ch. 8, Th. 1, t=0] that, if 2*n!*U(n) is a solution of an analogous problem at a circular table, then U(n) = V(n) - V(n-1), n>=3, where U(n) = A000179(n). Thus V(n) = Sum_{i=3,...,n} A000179(i), n>=1, and comparing the initial conditions, we conclude that a(n) = V(n), n>=1. This gives a combinatorial interpretation for 2*n!*a(n).
(End)

			G.f. = 1 + x^3 + 3*x^4 + 16*x^5 + 96*x^6 + 675*x^7 + 5413*x^8 + ...

W. Ahrens, Mathematische Unterhaltungen und Spiele. Teubner, Leipzig, Vol. 1, 3rd ed., 1921; Vol. 2, 2nd ed., 1918. See Vol. 2, p. 79.
J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 198.
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).
H. M. Taylor, A problem on arrangements, Mess. Math., 32 (1902), 60ff.

Seiichi Manyama, Table of n, a(n) for n = 0..450 (terms 0..100 from T. D. Noe)
W. Ahrens, Mathematische Unterhaltungen und Spiele, Leipzig: B. G. Teubner, 1901.
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122.
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122. [Annotated scanned copy]
Dmitry Efimov, Determinants of generalized binary band matrices, arXiv:1702.05655 [math.RA], 2017.
V. Kotesovec, Non-attacking chess pieces, 6ed, 2013, p. 222.
Y. Li, Ménage Numbers and Ménage Permutations, J. Int. Seq. 18 (2015) 15.6.8
H. M. Taylor, A problem on arrangements, Mess. Math., 32 (1902), 60ff. [Annotated scanned copy]
M. Wyman and L. Moser, On the problème des ménages, Canad. J. Math., 10 (1958), 468-480.

Cf. A000179, A000904, A001883, A137886, A292574. A diagonal of A058057.

[ &+[(-1)^(n-k)*Binomial(n+k, 2*k)*Factorial(k): k in [0..n]]: n in [0..21]]; // Bruno Berselli, Apr 11 2011

V := proc(n) local k; add( binomial(2*n-k,k)*(n-k)!*(x-1)^k, k=0..n); end; W := proc(r,s) coeff( V(r),x,s ); end; A000271 := n->W(n-2,0);

Table[Sum[(-1)^(n - k) k! Binomial[n + k, 2 k], {k, 0, n}], {n, 0, 22}] (* Jean-François Alcover, Apr 11 2011, after Paul Barry *)
RecurrenceTable[{a[0] == 1, a[1] == a[2] == 0, a[n] == (n - 1) a[n - 2] + (n - 1) a[n - 1] +  a[n - 3]}, a, {n, 30}] (* Harvey P. Dale, Jun 01 2012 *)
Table[(-1)^n HypergeometricPFQ[{1, -n, n + 1}, {1/2}, 1/4], {n, 20}] (* Michael Somos, May 28 2014 *)

a(n) = if(n, round( 2*exp(-2)*(besselk(n+1,2) + besselk(n,2)) ), 1) \\ Charles R Greathouse IV, May 11 2016

a(n) = (n - 1) a(n - 2) + (n - 1) a(n - 1) + a(n - 3).
From Paul Barry, Feb 08 2009: (Start)
G.f.: 1/(1+x-x/(1+x-x/(1+x-2x/(1+x-2x/(1+x-3x/(1+x-3x/(1+x-4x/(1+... (continued fraction);
a(n) = Sum_{k=0..n} binomial(2n-k,k)*(n-k)!*(-1)^k. (End)
a(n) = (-1)^n*hypergeom([1, -n, n+1],[1/2],1/4). - Mark van Hoeij, Nov 12 2009
a(n) = round( 2*exp(-2)*(BesselK(1+n,2) + BesselK(n,2)) ) for n>0. - Mark van Hoeij, Nov 12 2009
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(n+k,2*k)*k!. - Paul Barry, Jun 23 2010
G.f.: Sum_{n>=0} n!*x^n/(1+x)^(2*n+1). - Ira M. Gessel, Jan 15 2013
a(n) ~ exp(-2)*n!. - Vaclav Kotesovec, Mar 10 2014
a(-1 - n) = -a(n) for all n in Z. - Michael Somos, May 28 2014
a(n) = Sum_{i=3..n} A000179(i), n>=1. - Vladimir Shevelev, Jun 21 2015
0 = a(n)*(-a(n+2) - a(n+3)) + a(n+1)*(+a(n+1) + 2*a(n+2) + a(n+3) - a(n+4)) + a(n+2)*(+a(n+2) + 2*a(n+3) - a(n+4)) + a(n+3)*(+a(n+3)) for all n in Z. - Michael Somos, Oct 16 2016

More terms from James Sellers, Aug 21 2000

More terms from Simone Severini, Oct 14 2004

A285439 Sum T(n,k) of the entries in the k-th cycles of all permutations of [n]; triangle T(n,k), n>=1, 1<=k<=n, read by rows.

Original entry on oeis.org

1, 4, 2, 21, 12, 3, 132, 76, 28, 4, 960, 545, 235, 55, 5, 7920, 4422, 2064, 612, 96, 6, 73080, 40194, 19607, 6692, 1386, 154, 7, 745920, 405072, 202792, 75944, 18736, 2816, 232, 8, 8346240, 4484808, 2280834, 911637, 254061, 46422, 5256, 333, 9

Offset: 1

Alois P. Heinz, Apr 19 2017

Each cycle is written with the smallest element first and cycles are arranged in increasing order of their first elements.

			T(3,1) = 21 because the sum of the entries in the first cycles of all permutations of [3] ((123), (132), (12)(3), (13)(2), (1)(23), (1)(2)(3)) is 6+6+3+4+1+1 = 21.
Triangle T(n,k) begins:
       1;
       4,      2;
      21,     12,      3;
     132,     76,     28,     4;
     960,    545,    235,    55,     5;
    7920,   4422,   2064,   612,    96,    6;
   73080,  40194,  19607,  6692,  1386,  154,   7;
  745920, 405072, 202792, 75944, 18736, 2816, 232, 8;
  ...

Alois P. Heinz, Rows n = 1..23, flattened
Wikipedia, Permutation

Columns k=1-2 give: A284816, A285489.
Row sums give A000142 * A000217 = A180119.
Main diagonal and first lower diagonal give: A000027, A006000 (for n>0).
Cf. A000290, A002775, A185105, A285362, A285382, A285793.

T:= proc(h) option remember; local b; b:=
      proc(n, l) option remember; `if`(n=0, [mul((i-1)!, i=l), 0],
        (p-> p+[0, (h-n+1)*p[1]*x^(nops(l)+1)])(b(n-1, [l[], 1]))+
         add((p-> p+[0, (h-n+1)*p[1]*x^j])(
         b(n-1, subsop(j=l[j]+1, l))), j=1..nops(l)))
      end: (p-> seq(coeff(p, x, i), i=1..n))(b(h, [])[2])
    end:
seq(T(n), n=1..10);

T[h_] := T[h] = Module[{b}, b[n_, l_] := b[n, l] = If[n == 0, {Product[(i - 1)!, {i, l}], 0}, # + {0, (h - n + 1)*#[[1]]*x^(Length[l] + 1)}&[b[n - 1, Append[l, 1]]] + Sum[# + {0, (h-n+1)*#[[1]]*x^j}&[b[n - 1, ReplacePart[ l, j -> l[[j]] + 1]]], {j, 1, Length[l]}]]; Table[Coefficient[#, x, i], {i, 1, n}]&[b[h, {}][[2]]]];
Table[T[n], {n, 1, 10}] // Flatten (* Jean-François Alcover, May 25 2018, translated from Maple *)

Sum_{k=1..n} k * T(n,k) = n^2 * n! = A002775(n).

A047922 Triangle of numbers a(n,k) = number of terms in n X n determinant with 2 adjacent diagonals of k and k-1 0's (0<=k<=n).

Original entry on oeis.org

1, 1, 0, 2, 1, 0, 6, 4, 1, 1, 24, 18, 8, 5, 3, 120, 96, 54, 34, 23, 16, 720, 600, 384, 258, 182, 131, 96, 5040, 4320, 3000, 2136, 1566, 1168, 883, 675, 40320, 35280, 25920, 19320, 14664, 11274, 8756, 6859, 5413, 362880, 322560, 246960, 190800, 149160, 117696, 93582, 74902, 60301, 48800

Offset: 0

N. J. A. Sloane

			Triangle starts:
  1;
  1, 0;
  2, 1, 0;
  6, 4, 1, 1;
  ...

Alois P. Heinz, Rows n = 0..140, flattened
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122.
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122. [Annotated scanned copy]

Columns give A000142, A001563, A002775, A002776. Cf. A047920.

a:= proc(n, k) option remember; `if`(k=0, n!, `if`(n=k,
      `if`(n<3, (n-1)*(n-2)/2, (n-1)*(a(n-1$2)+a(n-2$2))
      +a(n-3$2)), a(n, k+1) +2*a(n-1, k) +a(n-2, k-1)))
    end:
seq(seq(a(n, k), k=0..n), n=0..10);  # Alois P. Heinz, Jun 24 2017

a[n_, n_] := (-1)^n*HypergeometricPFQ[{1, -n, n+1}, {1/2}, 1/4]; a[n_, k_] := a[n, k] = a[n, k+1] + 2*a[n-1, k] + a[n-2, k-1]; Table[a[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 24 2015 *)

Right diagonal is A000271, column k=0 is A000142; other entries given by a(n, k) = a(n, k+1) + 2a(n-1, k) + a(n-2, k-1).

More terms from Larry Reeves (larryr(AT)acm.org), Sep 29 2000

A091363 a(n) = n!*n^3.

Original entry on oeis.org

0, 1, 16, 162, 1536, 15000, 155520, 1728720, 20643840, 264539520, 3628800000, 53129260800, 827714764800, 13680764697600, 239217231052800, 4413400992000000, 85699747381248000, 1747492334235648000, 37338643451805696000, 834363743704178688000

Offset: 0

Mario Catalani (mario.catalani(AT)unito.it), Jan 07 2004

Denominators in the power series expansion of the higher order exponential integral E(x,3,1) + (gamma^3/6+Pi^2*gamma/36+zeta(3)/3+Pi^2*gamma/18) + (gamma^2/2+Pi^2/12)*log(x) + gamma*log(x)^2/2 + log(x)^3/6, n>0. See A163931 for information on the E(x,m,n). - Johannes W. Meijer, Oct 16 2009

Cf. A163931 (E(x,m,n)), A001563 (n*n!), A002775 (n^2*n!), A091364 (n^4*n!). - Johannes W. Meijer, Oct 16 2009

[Factorial(n)*n^3: n in [0..40]]; // Vincenzo Librandi, Jun 25 2015

a:=n->sum(sum(sum((n!), j=1..n),k=1..n),m=1..n): seq(a(n), n=0..17); # Zerinvary Lajos, May 16 2007

Mathematica
```
Table[n!n^3, {n, 0, 20}]
```

E.g.f.: (x+4x^2+x^3)/(1-x)^4.

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

Original entry on oeis.org

0, 1, 32, 486, 6144, 75000, 933120, 12101040, 165150720, 2380855680, 36288000000, 584421868800, 9932577177600, 177849941068800, 3349041234739200, 66201014880000000, 1371195958099968000, 29707369682006016000

Offset: 0

Mario Catalani (mario.catalani(AT)unito.it), Jan 07 2004

Denominators in the power series expansion of the higher order exponential integral E(x,4,1) - ((gamma^4/24+Pi^2*gamma^2/24+zeta(3)*gamma/3+Pi^4/160) + (gamma^3/6+ Pi^2*gamma/12+ zeta(3)/3)*log(x) + (gamma^2/4+ Pi^2/24)*log(x)^2 + (gamma/6)*log(x)^3 + log(x)^4/24), n>0. See A163931 for information on the E(x,m,n). - Johannes W. Meijer, Oct 16 2009

Cf. A163931 (E(x,m,n)), A001563 (n*n!), A002775 (n^2*n!), A091363 (n^3*n!). - Johannes W. Meijer, Oct 16 2009

a:=n->sum(sum(sum((n+1)!-n!, j=1..n),k=1..n),m=1..n): seq(a(n), n=0..17); # Zerinvary Lajos, May 16 2007

Mathematica
```
Table[n!n^4, {n, 0, 20}]
```

E.g.f.: (x + 11x^2 + 11x^3 + x^4)/(1 - x)^5

More terms from Zerinvary Lajos, May 16 2007

A002776 Terms in certain determinants.

Original entry on oeis.org

1, 5, 34, 258, 2136, 19320, 190800, 2051280, 23909760, 300827520, 4067884800, 58877280000, 908666035200, 14901260774400, 258832346572800, 4748165630208000, 91746433658880000, 1862735060938752000, 39649900359573504000, 883021783867711488000

Offset: 0

N. J. A. Sloane

a(n) equals (n+1)^2 times the permanent of the (n+1) X (n+1) matrix with 1/(n+1) in the top right corner, 1/(n+1) in the bottom left corner, and 1's everywhere else. - John M. Campbell, May 25 2011

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).

Muniru A Asiru, Table of n, a(n) for n = 0..100
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122.
J. D. H. Dickson, Discussion of two double series arising from the number of terms in determinants of certain forms, Proc. London Math. Soc., 10 (1879), 120-122. [Annotated scanned copy]

Cf. A047922.

A002776 := List([0..100], n -> Factorial(n+3) - 5 * Factorial(n+2) + 6 * Factorial(n+1) - Factorial(n)); # Muniru A Asiru, Jan 15 2018

[(n^3+n^2+2*n+1)*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Jan 19 2018

A002776 := [seq(factorial(n+3) - 5 * factorial(n+2) + 6 * factorial(n+1) - factorial(n), n=0..100)]; # Muniru A Asiru, Jan 15 2018

Table[(n^3+n^2+2n+1)n!,{n,0,30}] (* Harvey P. Dale, Oct 28 2011 *)

a(n) = (n^3 + n^2 + 2*n + 1)*n!.
a(n) = (n+3)! - 5*(n+2)! + 6*(n+1)! - n!. - Umut Özer, Dec 26 2017
E.g.f.: (1 + x + 3*x^2 + x^3)/(1 - x)^4. - Stefano Spezia, Apr 17 2022

Edited by Dean Hickerson, Sep 20 2002

cp's OEIS Frontend

A001563 a(n) = n*n! = (n+1)! - n!.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

Sage

Formula

A163931 Decimal expansion of the higher-order exponential integral E(x, m=2, n=1) at x=1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A047920 Triangular array formed from successive differences of factorial numbers.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A000271 Sums of ménage numbers.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

A285439 Sum T(n,k) of the entries in the k-th cycles of all permutations of [n]; triangle T(n,k), n>=1, 1<=k<=n, read by rows.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A047922 Triangle of numbers a(n,k) = number of terms in n X n determinant with 2 adjacent diagonals of k and k-1 0's (0<=k<=n).

Views

Author

Keywords

Examples