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 25 results. Next

A129841 Antidiagonal sums of triangle T defined in A048594: T(j,k) = k! * Stirling1(j,k), 1<= k <= j.

Original entry on oeis.org

1, -1, 4, -12, 52, -256, 1502, -10158, 78360, -680280, 6574872, -70075416, 816909816, -10342968456, 141357740736, -2074340369088, 32530886655168, -542971977209760, 9610316495698416, -179788450082431536, 3544714566466060032
Offset: 1

Views

Author

Paul Curtz, May 22 2007

Keywords

Examples

			First seven rows of T are
[    1 ]
[   -1,      2 ]
[    2,     -6,      6 ]
[   -6,     22,    -36,     24 ]
[   24,   -100,    210,   -240,    120 ]
[ -120,    548,  -1350,   2040,  -1800,    720 ]
[  720,  -3528,   9744, -17640,  21000, -15120,   5040 ]
		

References

  • P. Curtz, Integration numerique des systemes differentiels a conditions initiales. Note no. 12 du Centre de Calcul Scientifique de l'Armement, 1969, 135 pages, p. 61. Available from Centre d'Electronique de L'Armement, 35170 Bruz, France, or INRIA, Projets Algorithmes, 78150 Rocquencourt.
  • P. Curtz, Gazette des Mathematiciens, 1992, no. 52, p. 44.
  • P. Flajolet, X. Gourdon and B. Salvy, Gazette des Mathematiciens, 1993, no. 55, pp. 67-78.

Crossrefs

Cf. A048594 (T read by rows), A075181 (T unsigned with rows read backwards), A006252 (row sums of T), A000142 (main diagonal of T), A001286 (unsigned first subdiagonal of T). Unsigned values of second through sixth column of T are in A052517, A052748, A052753, A052767, A052779 resp.

Programs

  • Magma
    m:=21; T:=[ [ Factorial(k)*StirlingFirst(j, k): k in [1..j] ]: j in [1..m] ]; [ &+[ T[j-k+1][k]: k in [1..(j+1) div 2] ]: j in [1..m] ]; // Klaus Brockhaus, Jun 03 2007
  • Mathematica
    m = 21; t[j_, k_] := k!*StirlingS1[j, k]; Total /@ Table[ t[j-k+1, k], {j, 1, m}, {k, 1, Quotient[j+1, 2]}] (* Jean-François Alcover, Aug 13 2012, translated from Klaus Brockhaus's Magma program *)

Formula

E.g.f. for k-th column (k>=1): log(1+x)^k. For further formulas see the references.

Extensions

Edited and extended by Klaus Brockhaus, Jun 03 2007

A075182 Greatest common divisors of rows of triangle A075181 and of (unsigned) triangle A048594.

Original entry on oeis.org

1, 1, 2, 2, 2, 2, 24, 24, 24, 24, 48, 48, 48, 48, 384, 1152, 1152, 1152, 2304, 2304, 11520, 11520, 46080, 46080, 414720, 414720, 829440, 829440, 829440, 829440, 13271040, 13271040, 13271040, 39813120, 79626240, 79626240, 79626240, 79626240, 318504960, 318504960
Offset: 1

Views

Author

Wolfdieter Lang, Sep 19 2002

Keywords

Examples

			Row n=3 of triangle A075181 is [6,6,2], hence a(3)=2.
		

Crossrefs

Programs

  • Mathematica
    a[n_] := GCD @@ Table[k! * Abs[StirlingS1[n, k]], {k, 1, n}]; Array[a, 40] (* Amiram Eldar, Aug 08 2024 *)

Formula

a(n) = gcd(A075181(n, m), m=1..n) = gcd(|A048594(n, m)|, m=1..n), n>=1.

Extensions

More terms from Amiram Eldar, Aug 08 2024

A019538 Triangle of numbers T(n,k) = k!*Stirling2(n,k) read by rows (n >= 1, 1 <= k <= n).

Original entry on oeis.org

1, 1, 2, 1, 6, 6, 1, 14, 36, 24, 1, 30, 150, 240, 120, 1, 62, 540, 1560, 1800, 720, 1, 126, 1806, 8400, 16800, 15120, 5040, 1, 254, 5796, 40824, 126000, 191520, 141120, 40320, 1, 510, 18150, 186480, 834120, 1905120, 2328480, 1451520, 362880, 1, 1022, 55980, 818520, 5103000, 16435440, 29635200, 30240000, 16329600, 3628800
Offset: 1

Views

Author

N. J. A. Sloane and Manfred Goebel (goebel(AT)informatik.uni-tuebingen.de), Dec 11 1996

Keywords

Comments

Number of ways n labeled objects can be distributed into k nonempty parcels. Also number of special terms in n variables with maximal degree k.
In older terminology these are called differences of 0. - Michael Somos, Oct 08 2003
Number of surjections (onto functions) from an n-element set to a k-element set.
Also coefficients (in ascending order) of so-called ordered Bell polynomials.
(k-1)!*Stirling2(n,k-1) is the number of chain topologies on an n-set having k open sets [Stephen].
Number of set compositions (ordered set partitions) of n items into k parts. Number of k dimensional 'faces' of the n dimensional permutohedron (see Simion, p. 162). - Mitch Harris, Jan 16 2007
Correction of comment before: Number of (n-k)-dimensional 'faces' of the permutohedron of order n (an (n-1)-dimensional polytope). - Tilman Piesk, Oct 29 2014
This array is related to the reciprocal of an e.g.f. as sketched in A133314. For example, the coefficient of the fourth-order term in the Taylor series expansion of 1/(a(0) + a(1) x + a(2) x^2/2! + a(3) x^3/3! + ...) is a(0)^(-5) * {24 a(1)^4 - 36 a(1)^2 a(2) a(0) + [8 a(1) a(3) + 6 a(2)^2] a(0)^2 - a(4) a(0)^3}. The unsigned coefficients characterize the P3 permutohedron depicted on page 10 in the Loday link with 24 vertices (0-D faces), 36 edges (1-D faces), 6 squares (2-D faces), 8 hexagons (2-D faces) and 1 3-D permutohedron. Summing coefficients over like dimensions gives A019538 and A090582. Compare to A133437 for the associahedron. - Tom Copeland, Sep 29 2008, Oct 07 2008
Further to the comments of Tom Copeland above, the permutohedron of type A_3 can be taken as the truncated octahedron. Its dual is the tetrakis hexahedron, a simplicial polyhedron, with f-vector (1,14,36,24) giving the fourth row of this triangle. See the Wikipedia entry and [Fomin and Reading p. 21]. The corresponding h-vectors of permutohedra of type A give the rows of the triangle of Eulerian numbers A008292. See A145901 and A145902 for the array of f-vectors for type B and type D permutohedra respectively. - Peter Bala, Oct 26 2008
Subtriangle of triangle in A131689. - Philippe Deléham, Nov 03 2008
Since T(n,k) counts surjective functions and surjective functions are "consistent", T(n,k) satisfies a binomial identity, namely, T(n,x+y) = Sum_{j=0..n} C(n,j)*T(j,x)*T(n-j,y). For definition of consistent functions and a generalized binomial identity, see "Toy stories and combinatorial identities" in the link section below. - Dennis P. Walsh, Feb 24 2012
T(n,k) is the number of labeled forests on n+k vertices satisfying the following two conditions: (i) each forest consists of exactly k rooted trees with roots labeled 1, 2, ..., k; (ii) every root has at least one child vertex. - Dennis P. Walsh, Feb 24 2012
The triangle is the inverse binomial transform of triangle A028246, deleting the left column and shifting up one row. - Gary W. Adamson, Mar 05 2012
See A074909 for associations among this array and the Bernoulli polynomials and their umbral compositional inverses. - Tom Copeland, Nov 14 2014
E.g.f. for the shifted signed polynomials is G(x,t) = (e^t-1)/[1+(1+x)(e^t-1)] = 1-(1+x)(e^t-1) + (1+x)^2(e^t-1)^2 - ... (see also A008292 and A074909), which has the infinitesimal generator g(x,u)d/du = [(1-x*u)(1-(1+x)u)]d/du, i.e., exp[t*g(x,u)d/du]u eval. at u=0 gives G(x,t), and dG(x,t)/dt = g(x,G(x,t)). The compositional inverse is log((1-xt)/(1-(1+x)t)). G(x,t) is a generating series associated to the generalized Hirzebruch genera. See the G. Rzadowski link for the relation of the derivatives of g(x,u) to solutions of the Riccatt differential equation, soliton solns. to the KdV equation, and the Eulerian and Bernoulli numbers. In addition A145271 connects products of derivatives of g(x,u) and the refined Eulerian numbers to the inverse of G(x,t), which gives the normalized, reverse face polynomials of the simplices (A135278, divided by n+1). See A028246 for the generator g(x,u)d/dx. - Tom Copeland, Nov 21 2014
For connections to toric varieties and Eulerian polynomials, see the Dolgachev and Lunts and the Stembridge links. - Tom Copeland, Dec 31 2015
See A008279 for a relation between the e.g.f.s enumerating the faces of permutahedra (this entry) and stellahedra. - Tom Copeland, Nov 14 2016
T(n, k) appears in a Worpitzky identity relating monomials to binomials: x^n = Sum_{k=1..n} T(n, k)*binomial(x,k), n >= 1. See eq. (11.) of the Worpitzky link on p. 209. The relation to the Eulerian numbers is given there in eqs. (14.) and (15.). See the formula below relating to A008292. See also Graham et al. eq. (6.10) (relating monomials to falling factorials) on p. 248 (2nd ed. p. 262). The Worpitzky identity given in the Graham et al. reference as eq. (6.37) (2nd ed. p. 269) is eq. (5.), p. 207, of Worpitzky. - Wolfdieter Lang, Mar 10 2017
T(n, m) is also the number of minimum clique coverings and minimum matchings in the complete bipartite graph K_{m,n}. - Eric W. Weisstein, Apr 26 2017
From the Hasan and Franco and Hasan papers: The m-permutohedra for m=1,2,3,4 are the line segment, hexagon, truncated octahedron and omnitruncated 5-cell. The first three are well-known from the study of elliptic models, brane tilings and brane brick models. The m+1 torus can be tiled by a single (m+2)-permutohedron. Relations to toric Calabi-Yau Kahler manifolds are also discussed. - Tom Copeland, May 14 2020
From Manfred Boergens, Jul 25 2021: (Start)
Number of n X k binary matrices with row sums = 1 and no zero columns. These matrices are a subset of the matrices defining A183109.
The distribution into parcels in the leading comment can be regarded as a covering of [n] by tuples (A_1,...,A_k) in P([n])^k with nonempty and disjoint A_j, with P(.) denoting the power set (corrected for clarity by Manfred Boergens, May 26 2024). For the non-disjoint case see A183109 and A218695.
For tuples with "nonempty" dropped see A089072. For tuples with "nonempty and disjoint" dropped see A092477 and A329943 (amendment by Manfred Boergens, Jun 24 2024). (End)

Examples

			The triangle T(n, k) begins:
  n\k 1    2     3      4       5        6        7        8        9      10
  1:  1
  2:  1    2
  3:  1    6     6
  4:  1   14    36     24
  5:  1   30   150    240     120
  6:  1   62   540   1560    1800      720
  7:  1  126  1806   8400   16800    15120     5040
  8:  1  254  5796  40824  126000   191520   141120    40320
  9:  1  510 18150 186480  834120  1905120  2328480  1451520   362880
  10: 1 1022 55980 818520 5103000 16435440 29635200 30240000 16329600 3628800
  ... Reformatted and extended - _Wolfdieter Lang_, Oct 04 2014
---------------------------------------------------------------------------
T(4,1) = 1: {1234}. T(4,2) = 14: {1}{234} (4 ways), {12}{34} (6 ways), {123}{4} (4 ways). T(4,3) = 36: {12}{3}{4} (12 ways), {1}{23}{4} (12 ways), {1}{2}{34} (12 ways). T(4,4) = 1: {1}{2}{3}{4} (1 way).
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 89, ex. 1; also p. 210.
  • Miklos Bona, Combinatorics of Permutations, Chapman and Hall,2004, p.12.
  • G. Boole, A Treatise On The Calculus of Finite Differences, Dover Publications, 1960, p. 20.
  • H. T. Davis, Tables of the Mathematical Functions. Vols. 1 and 2, 2nd ed., 1963, Vol. 3 (with V. J. Fisher), 1962; Principia Press of Trinity Univ., San Antonio, TX, Vol. 2, p. 212.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, 1989, p. 155. Also eqs.(6.10) and (6.37).
  • Kiran S. Kedlaya and Andrew V. Sutherland, Computing L -Series of Hyperelliptic Curves in Algorithmic Number Theory Lecture Notes in Computer Science Volume 5011/2008.
  • T. K. Petersen, Eulerian Numbers, Birkhauser, 2015, Section 5.6.
  • J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 33.
  • J. F. Steffensen, Interpolation, 2nd ed., Chelsea, NY, 1950, see p. 54.
  • A. H. Voigt, Theorie der Zahlenreihen und der Reihengleichungen, Goschen, Leipzig, 1911, p. 31.
  • E. Whittaker and G. Robinson, The Calculus of Observations, Blackie, London, 4th ed., 1949; p. 7.

Crossrefs

Row sums give A000670. Maximal terms in rows give A002869. Central terms T(2k-1,k) give A233734.
Diagonal is n! (A000142). 2nd diagonal is A001286. 3rd diagonal is A037960.
Reflected version of A090582. A371568 is another version.
See also the two closely related triangles: A008277(n, k) = T(n, k)/k! (Stirling numbers of second kind) and A028246(n, k) = T(n, k)/k.
Cf. A033282 'faces' of the associahedron.
Cf. A008292, A047969, A145901, A145902. - Peter Bala, Oct 26 2008
Visible in the 3-D array in A249042.
See also A000182.

Programs

  • Haskell
    a019538 n k = a019538_tabl !! (n-1) !! (k-1)
    a019538_row n = a019538_tabl !! (n-1)
    a019538_tabl = iterate f [1] where
       f xs = zipWith (*) [1..] $ zipWith (+) ([0] ++ xs) (xs ++ [0])
    -- Reinhard Zumkeller, Dec 15 2013
    
  • Maple
    with(combinat): A019538 := (n,k)->k!*stirling2(n,k);
  • Mathematica
    Table[k! StirlingS2[n, k], {n, 9}, {k, n}] // Flatten
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, sum(i=0, k, (-1)^i * binomial(k, i) * (k-i)^n))}; /* Michael Somos, Oct 08 2003 */
    
  • Sage
    def T(n, k): return factorial(k)*stirling_number2(n,k) # Danny Rorabaugh, Oct 10 2015

Formula

T(n, k) = k*(T(n-1, k-1)+T(n-1, k)) with T(0, 0) = 1 [or T(1, 1) = 1]. - Henry Bottomley, Mar 02 2001
E.g.f.: (y*(exp(x)-1) - exp(x))/(y*(exp(x)-1) - 1). - Vladeta Jovovic, Jan 30 2003
Equals [0, 1, 0, 2, 0, 3, 0, 4, 0, 5, ...] DELTA [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...] where DELTA is Deléham's operator defined in A084938.
T(n, k) = Sum_{j=0..k} (-1)^(k-j)*j^n*binomial(k, j). - Mario Catalani (mario.catalani(AT)unito.it), Nov 28 2003. See Graham et al., eq. (6.19), p. 251. For a proof see Bert Seghers, Jun 29 2013.
Sum_{k=0..n} T(n, k)(-1)^(n-k) = 1, Sum_{k=0..n} T(n, k)(-1)^k = (-1)^n. - Mario Catalani (mario.catalani(AT)unito.it), Dec 11 2003
O.g.f. for n-th row: polylog(-n, x/(1+x))/(x+x^2). - Vladeta Jovovic, Jan 30 2005
E.g.f.: 1 / (1 + t*(1-exp(x))). - Tom Copeland, Oct 13 2008
From Peter Bala, Oct 26 2008: (Start)
O.g.f. as a continued fraction: 1/(1 - x*t/(1 - (x + 1)*t/(1 - 2*x*t/(1 - 2*(x + 1)*t/(1 - ...))))) = 1 + x*t + (x + 2*x^2)*t^2 + (x + 6*x^2 + 6*x^3)*t^3 + ... .
The row polynomials R(n,x), which begin R(1,x) = x, R(2,x) = x + 2*x^2, R(3,x) = x + 6*x^2 + 6*x^3, satisfy the recurrence x*d/dx ((x + 1)*R(n,x)) = R(n+1,x). It follows that the zeros of R(n,x) are real and negative (apply Corollary 1.2 of [Liu and Wang]).
Since this is the triangle of f-vectors of the (simplicial complexes dual to the) type A permutohedra, whose h-vectors form the Eulerian number triangle A008292, the coefficients of the polynomial (x-1)^n*R(n,1/(x-1)) give the n-th row of A008292. For example, from row 3 we have x^2 + 6*x + 6 = 1 + 4*y + y^2, where y = x + 1, producing [1,4,1] as the third row of A008292. The matrix product A008292 * A007318 gives the mirror image of this triangle (see A090582).
For n,k >= 0, T(n+1,k+1) = Sum_{j=0..k} (-1)^(k-j)*binomial(k,j)*[(j+1)^(n+1) - j^(n+1)]. The matrix product of Pascal's triangle A007318 with the current array gives (essentially) A047969. This triangle is also related to triangle A047969 by means of the S-transform of [Hetyei], a linear transformation of polynomials whose value on the basis monomials x^k is given by S(x^k) = binomial(x,k). The S-transform of the shifted n-th row polynomial Q(n,x) := R(n,x)/x is S(Q(n,x)) = (x+1)^n - x^n. For example, from row 3 we obtain S(1 + 6*x + 6*x^2) = 1 + 6*x + 6*x*(x-1)/2 = 1 + 3*x + 3*x^2 = (x+1)^3 - x^3. For fixed k, the values S(Q(n,k)) give the nonzero entries in column (k-1) of the triangle A047969 (the Hilbert transform of the Eulerian numbers). (End)
E.g.f.: (exp(x)-1)^k = sum T(n,k)x^n/n!. - Vladimir Kruchinin, Aug 10 2010
T(n,k) = Sum_{i=1..k} A(n,i)*Binomial(n-i,k-i) where A(n,i) is the number of n-permutations that have i ascending runs, A008292.
From Tom Copeland, Oct 11 2011: (Start)
With e.g.f. A(x,t) = -1 + 1/(1+t*(1-exp(x))), the comp. inverse in x is B(x,t) = log(((1+t)/t) - 1/(t(1+x))).
With h(x,t) = 1/(dB/dx)= (1+x)((1+t)(1+x)-1), the row polynomial P(n,t) is given by (h(x,t)*d/dx)^n x, eval. at x=0, A=exp(x*h(y,t)*d/dy) y, eval. at y=0, and dA/dx = h(A(x,t),t), with P(0,t)=0.
(A factor of -1/n! was removed by Copeland on Aug 25 2016.) (End)
The term linear in x of [x*h(d/dx,t)]^n 1 gives the n-th row polynomial. (See A134685.) - Tom Copeland, Nov 07 2011
Row polynomials are given by D^n(1/(1-x*t)) evaluated at x = 0, where D is the operator (1+x)*d/dx. - Peter Bala, Nov 25 2011
T(n,x+y) = Sum_{j=0..n} binomial(n,j)*T(j,x)*T(n-j,y). - Dennis P. Walsh, Feb 24 2012
Let P be a Rota-Baxter operator of weight 1 satisfying the identity P(x)*P(y) = P(P(x)*y) + P(x*P(y)) + P(x*y). Then P(1)^2 = P(1) + 2*P^2(1). More generally, Guo shows that P(1)^n = Sum_{k=1..n} T(n,k)*P^k(1). - Peter Bala, Jun 08 2012
Sum_{i=1..n} (-1)^i*T(n,i)/i = 0, for n > 1. - Leonid Bedratyuk, Aug 09 2012
T(n, k) = Sum_{j=0..k} (-1)^j*binomial(k, j)*(k-j)^n. [M. Catalani's re-indexed formula from Nov 28 2003] Proof: count the surjections of [n] onto [k] with the inclusion-exclusion principle, as an alternating sum of the number of functions from [n] to [k-j]. - Bert Seghers, Jun 29 2013
n-th row polynomial = 1/(1 + x)*( Sum_{k>=0} k^n*(x/(1 + x))^k ), valid for x in the open interval (-1/2, inf). See Tanny link. Cf. A145901. - Peter Bala, Jul 22 2014
T(n,k) = k * A141618(n,k-1) / binomial(n,k-1). - Tom Copeland, Oct 25 2014
Sum_{n>=0} n^k*a^n = Sum_{i=1..k} (a / (1 - a))^i * T(k, i)/(1-a) for |a| < 1. - David A. Corneth, Mar 09 2015
From Peter Bala, May 26 2015: (Start)
The row polynomials R(n,x) satisfy (1 + x)*R(n,x) = (-1)^n*x*R(n,-(1 + x)).
For a fixed integer k, the expansion of the function A(k,z) := exp( Sum_{n >= 1} R(n,k)*z^n/n ) has integer coefficients and satisfies the functional equation A(k,z)^(k + 1) = BINOMIAL(A(k,z))^k, where BINOMIAL(F(z))= 1/(1 - z)*F(z/(1 - z)) denotes the binomial transform of the o.g.f. F(z). Cf. A145901. For cases see A084784 (k = 1), A090352 (k = 2), A090355 (k = 3), A090357 (k = 4), A090362 (k = 5) and A084785 (k = -2 with z -> -z).
A(k,z)^(k + 1) = A(-(k + 1),-z)^k and hence BINOMIAL(A(k,z)) = A(-(k + 1),-z). (End)
From Tom Copeland, Oct 19 2016: (Start)
Let a(1) = 1 + x + B(1) = x + 1/2 and a(n) = B(n) = (B.)^n, where B(n) are the Bernoulli numbers defined by e^(B.t) = t / (e^t-1), then t / e^(a.t) = t / [(x + 1) * t + exp(B.t)] = (e^t - 1) /[ 1 + (x + 1) (e^t - 1)] = exp(p.(x)t), where (p.(x))^n = p_n(x) are the shifted, signed row polynomials of this array: p_0(x) = 0, p_1(x) = 1, p_2(x) = -(1 + 2 x), p_3(x) = 1 + 6 x + 6 x^2, ... and p_n(x) = n * b(n-1), where b(n) are the partition polynomials of A133314 evaluated with these a(n).
Sum_{n > 0} R(n,-1/2) x^n/n! = 2 * tanh(x/2), where R(n,x) = Sum_{k = 1..n} T(n,k) x^(k-1) are the shifted row polynomials of this entry, so R(n,-1/2) = 4 * (2^(n+1)-1) B(n+1)/(n+1). (Cf. A000182.)
(End)
Also the Bernoulli numbers are given by B(n) = Sum_{k =1..n} (-1)^k T(n,k) / (k+1). - Tom Copeland, Nov 06 2016
G.f. for column k: k! x^k / Product_{i=1..k} (1-i*x). - Robert A. Russell, Sep 25 2018
a(j) <= A183109(j). - Manfred Boergens, Jul 25 2021

A006252 Expansion of e.g.f. 1/(1 - log(1+x)).

Original entry on oeis.org

1, 1, 1, 2, 4, 14, 38, 216, 600, 6240, 9552, 319296, -519312, 28108560, -176474352, 3998454144, -43985078784, 837126163584, -12437000028288, 237195036797184, -4235955315745536, 85886259443020800, -1746536474655406080, 38320721602434017280, -864056965711935974400
Offset: 0

Views

Author

Keywords

Comments

From Michael Somos, Mar 04 2004: (Start)
Stirling transform of a(n+1)=[1,2,4,14,38,...] is A000255(n)=[1,3,11,53,309,...].
Stirling transform of 2*a(n)=[2,2,4,8,28,...] is A052849(n)=[2,4,12,48,240,...].
Stirling transform of a(n)=[1,1,2,4,14,38,216,...] is A000142(n)=[1,2,6,24,120,...].
Stirling transform of a(n-1)=[1,1,1,2,4,14,38,...] is A000522(n-1)=[1,2,5,16,65,...].
Stirling transform of a(n-1)=[0,1,1,2,4,14,38,...] is A007526(n-1)=[0,1,4,15,64,...].
(End)
For n > 0: a(n) = sum of n-th row in triangle A048594. - Reinhard Zumkeller, Mar 02 2014
Coefficients in a factorial series representation of the exponential integral: exp(z)*E_1(z) = Sum_{n >= 0} (-1)^n*a(n)/(z)n, where (z)_n denotes the rising factorial z*(z + 1)*...*(z + n) and E_1(z) = Integrate{t = z..inf} exp(-t)/t dt. See Weninger, equation 6.4. - Peter Bala, Feb 12 2019

References

  • G. Pólya, Induction and Analogy in Mathematics. Princeton Univ. Press, 1954, p. 9.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Column k=1 of A320080.
Cf. A007840.

Programs

  • Haskell
    a006252 0 = 1
    a006252 n = sum $ a048594_row n  -- Reinhard Zumkeller, Mar 02 2014
    
  • Mathematica
    With[{nn=30},CoefficientList[Series[1/(1-Log[1+x]),{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Aug 12 2016 *)
  • PARI
    a(n)=if(n<0,0,n!*polcoeff(1/(1-log(1+x+x*O(x^n))),n))
    
  • PARI
    {a(n)=local(CF=1+x*O(x^n)); for(k=0, n-1, CF=1/((n-k+1)-(n-k)*x+(n-k+1)^2*x*CF)); n!*polcoeff(1+x/(1-x+x*CF), n, x)} /* Paul D. Hanna, Dec 31 2011 */
    
  • PARI
    a_vector(n) = my(v=vector(n+1)); v[1]=1; for(i=1, n, v[i+1]=sum(j=1, i, (-1)^(j-1)*(j-1)!*binomial(i, j)*v[i-j+1])); v; \\ Seiichi Manyama, May 22 2022
    
  • Sage
    def A006252_list(len):
        f, R, C = 1, [1], [1]+[0]*len
        for n in (1..len):
            f *= n
            for k in range(n, 0, -1):
                C[k] = -C[k-1]*((k-1)/(k) if k>1 else 1)
            C[0] = -sum(C[k] for k in (1..n))
            R.append(C[0]*f)
        return R
    print(A006252_list(24)) # Peter Luschny, Feb 21 2016

Formula

a(n) = Sum_{k=0..n} k!*stirling1(n, k). - Vladeta Jovovic, Sep 08 2002
a(n) = D^n(1/(1-x)) evaluated at x = 0, where D is the operator exp(-x)*d/dx. Row sums of A048594. Cf. A007840. - Peter Bala, Nov 25 2011
E.g.f.: 1/(1-log(1+x)) = 1 + x/(1-x + x/(2-x + 4*x/(3-2*x + 9*x/(4-3*x + 16*x/(5-4*x + 25*x/(6-5*x +...)))))), a continued fraction. - Paul D. Hanna, Dec 31 2011
a(n)/n! ~ -(-1)^n / (n * (log(n))^2) * (1 - 2*(1 + gamma)/log(n)), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Jul 01 2018
a(0) = 1; a(n) = Sum_{k=1..n} (-1)^(k-1) * (k-1)! * binomial(n,k) * a(n-k). - Seiichi Manyama, May 22 2022

A072597 Expansion of 1/(exp(-x) - x) as exponential generating function.

Original entry on oeis.org

1, 2, 7, 37, 261, 2301, 24343, 300455, 4238153, 67255273, 1185860331, 23000296155, 486655768525, 11155073325917, 275364320099807, 7282929854486431, 205462851526617489, 6158705454187353297, 195465061563672788947, 6548320737474275229347, 230922973019493881984021
Offset: 0

Views

Author

Michael Somos, Jun 23 2002

Keywords

Comments

Polynomials from A140749/A141412 are linked to Stirling1 (see A048594, A129841, A140749). See also P. Flajolet, X. Gourdon, B. Salvy in, available on Internet, RR-1857.pdf (preprint of unavailable Gazette des Mathematiciens 55, 1993, pp. 67-78; for graph 2 see also X. Gourdon RR-1852.pdf, pp. 64-65). What is the corresponding graph for A152650/A152656 = simplified A009998/A119502 linked, via A152818, to a(n), then Stirling2? - Paul Curtz, Dec 16 2008
Denominators in rational approximations of Lambert W(1). See Ramanujan, Notebooks, volume 2, page 22: "2. If e^{-x} = x, shew that the convergents to x are 1/2, 4/7, 21/37, 148/261, &c." Numerators in A006153. - Michael Somos, Jan 21 2019
Call an element g in a semigroup a group element if g^j = g for some j > 1. Then a(n) is the number of group elements in the semigroup of partial transformations of an n-set. Hence a(n) = Sum_{k=0..n} A154372(n,k)*k!. - Geoffrey Critzer, Nov 27 2021

Examples

			G.f. = 1 + 2*x + 7*x^2 + 37*x^3 + 261*x^4 + 2301*x^5 + 24343*x^6 + ...
		

References

  • O. Ganyushkin and V Mazorchuk, Classical Finite Transformation Semigroups, Springer, 2009, page 70.
  • S. Ramanujan, Notebooks, Tata Institute of Fundamental Research, Bombay 1957 Vol. 2, see page 22.

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[1/(Exp[-x]-x), {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, Jun 26 2013 *)
    a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ 1 / (Exp[-x] - x), {x, 0, n}]]; (* Michael Somos, Jan 21 2019 *)
    a[ n_] := If[ n < 0, 0, n! Sum[ (n - k + 1)^k / k!, {k, 0, n}]]; (* Michael Somos, Jan 21 2019 *)
  • PARI
    {a(n) = if( n<0, 0, n! * polcoeff( 1 / (exp(-x + x * O(x^n)) - x), n))};
    
  • PARI
    {a(n) = if( n<0, 0, n! * sum(k=0, n, (n-k+1)^k / k!))}; /* Michael Somos, Jan 21 2019 */

Formula

E.g.f.: 1 / (exp(-x) - x).
a(n) = n!*Sum_{k=0..n} (n-k+1)^k/k!. - Vladeta Jovovic, Aug 31 2003
a(n) = Sum_{k=0..n} (-1)^(n-k)*Stirling2(n, k)*A052820(k). - Vladeta Jovovic, Apr 12 2004
Recurrence: a(n+1) = 1 + Sum_{j=1..n} binomial(n, j)*a(j)*j. - Jon Perry, Apr 25 2005
E.g.f.: 1/(Q(0) - x) where Q(k) = 1 - x/(2*k+1 - x*(2*k+1)/(x - (2*k+2)/Q(k+1) )); (continued fraction ). - Sergei N. Gladkovskii, Apr 04 2013
a(n) ~ n!/((1+c)*c^(n+1)), where c = A030178 = LambertW(1) = 0.5671432904... - Vaclav Kotesovec, Jun 26 2013
O.g.f.: Sum_{k>=0} k!*x^k/(1 - (k + 1)*x)^(k+1). - Ilya Gutkovskiy, Oct 09 2018
a(n) = A006153(n+1)/(n+1). - Seiichi Manyama, Nov 05 2024

A133942 a(n) = (-1)^n * n!.

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
Offset: 0

Views

Author

Michael Somos, Sep 30 2007

Keywords

Comments

A variant of A000142, the factorial numbers. - N. J. A. Sloane, Oct 03 2007
The terms of this sequences form the factorial series which Euler called the divergent series par excellence.
Euler summed this series to 0.596347... (A073003 = Gompertz's constant).
Sum_{n>=0} 1/a(n) = 1/e. - Jaume Oliver Lafont, Mar 03 2009
A002104(n+1) = p(-1) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 0, 1, ..., n. - Michael Somos, Apr 30 2012
a(n) = A048594(2*n+1, n+1). - Reinhard Zumkeller, Mar 02 2014
log(1+x) = Sum_{n>=1} a(n-1)/n!*x^n. - James R. Buddenhagen, May 24 2015
It seems that a(n) is the determinant of n+1 X n+1 matrix whose elements are m(i,j) = quotient(i/j) + remainder(i/j). - Andres Cicuttin, Feb 11 2018

Examples

			G.f. = 1 - x + 2*x^2 - 6*x^3 + 24*x^4 - 120*x^5 + 720*x^6 - 5040*x^7 + ...
		

References

  • 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)
  • R. Roy, Sources in the Development of Mathematics, Cambridge University Press, 2011. See p. 186.

Crossrefs

Partial sums are A058006.
Alternating row sums of A048994.
Also, a(n) = A048994(n+1,1).

Programs

  • GAP
    List([0..20],n->(-1)^n*Factorial(n)); # Muniru A Asiru, Oct 27 2018
  • Haskell
    a133942 n = a133942_list !! n
    a133942_list = zipWith (*) a000142_list $ cycle [1, -1]
    -- Reinhard Zumkeller, Mar 02 2014
    
  • Magma
    [(-1)^n * Factorial(n): n in [0..25]]; // Vincenzo Librandi, May 12 2011
    
  • Maple
    seq((-1)^n*factorial(n),n=0..20); # Muniru A Asiru, Oct 27 2018
  • Mathematica
    nn=20;CoefficientList[Series[1/(1+x),{x,0,nn}],x]Range[0,nn]! (* or *)
    RecurrenceTable[{a[0]==1,a[n]==-n*a[n-1]},a[n],{n,20}] (* Harvey P. Dale, May 10 2011 and slightly modified by Robert G. Wilson v, Feb 12 2018 *)
    a[n_] := (-1)^n*n!; Array[a, 22, 0] (* Robert G. Wilson v, Feb 11 2018 *)
    Times@@@Partition[Riffle[Range[0,30]!,{1,-1}],2] (* Harvey P. Dale, Dec 30 2019 *)
  • PARI
    {a(n) = if( n<0, 0, (-1)^n * n! )};
    
  • Python
    import math
    for n in range(0, 25): print((-1)**n*math.factorial(n), end=', ') # Stefano Spezia, Oct 27 2018
    

Formula

Sum_{i=0..n} (-1)^i * i^n * binomial(n, i) = (-1)^n * n!. - Yong Kong (ykong(AT)curagen.com), Dec 26 2000
Stirling transform of a(n) = [1, -1, 2, -6, 24, ...] is A000007(n) = [1, 0, 0, 0, 0, ...].
a(n) = -n * a(n-1) unless n=0. a(n) = (-1)^n * A000142(n).
E.g.f.: 1/(1 + x).
G.f.: integral(t=1/x,infinity, (e^-t)/t) e^(1/x)/x = 1/(1 + x/(1 + x/(1 + 2*x/(1 + 2*x/(1 + 3*x/(1 + 3*x/(1 + ...))))))).
Convolution inverse of A158882. HANKEL transform is A055209. PSUM transform is A058006. BIN1 transform is A002741(n+1). - Michael Somos, Apr 30 2012
G.f.: 1 - x/(G(0)+x) where G(k) = 1 + (k+1)*x/(1 + x*(k+2)/G(k+1)), G(0) = W(1,1;x)/W(1,2;x), 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)]; (continued fraction, 2-step). - Sergei N. Gladkovskii, Aug 14 2012
G.f.: 1/U(0) where U(k) = 1 + x*(k+1)/(1 + x*(k+1)/U(k+1)); (continued fraction, 2-step). - Sergei N. Gladkovskii, Oct 15 2012
a(n) = (-1)^n*det(S(i+1,j)|, 1 <= i,j <= n), where S(n,k) are Stirling numbers of the second kind. - Mircea Merca, Apr 06 2013
G.f.: 2/G(0), where G(k)= 1 + 1/(1 - 2*x*(k+1)/(2*x*(k+1) + 1 + 2*x*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
E.g.f.: 1/(1 + x)= G(0), where G(k) = 1 - x*(k+1)*(k+2)/(1 + (k+1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Jan 29 2014
For n >= 1, a(n) = round(zeta^(n)(2)), where zeta^(n) is the n-th derivative of the Riemann zeta function. - Iain Fox, Nov 13 2017
a(n) = (n+1)^(n+1) * Integral_{x=0..1} (x*log(x))^n dx. - Peter James Foreman, Oct 27 2018

A320083 Expansion of e.g.f. Sum_{k>=0} log(1 + k*x)^k.

Original entry on oeis.org

1, 1, 7, 116, 3574, 177094, 12873962, 1290494904, 170592253320, 28753159552272, 6018433850602848, 1531605185388897552, 465706857941949607008, 166746568516127626614288, 69440517484503283491716400, 33278913978673363553703249408, 18185279212166784139689388753536, 11239676837467731657648932618952576
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 05 2018

Keywords

Crossrefs

Programs

  • Maple
    1,seq(n!*coeff(series(add( log(1 + k*x)^k,k=1..100), x=0, 18), x, n), n=1..17); # Paolo P. Lava, Jan 09 2019
  • Mathematica
    nmax = 17; CoefficientList[Series[1 + Sum[Log[1 + k x]^k, {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
    Join[{1}, Table[Sum[StirlingS1[n, k] k! k^n, {k, n}], {n, 17}]]

Formula

a(n) = Sum_{k=0..n} Stirling1(n,k)*k!*k^n.
a(n) ~ c * d^n * n^(2*n + 1/2), where d = 0.298212940253960977992575968955431001807757948758929... and c = 3.40415549717199390989204785905061856492539214306... - Vaclav Kotesovec, Oct 05 2018

A225479 Triangle read by rows, the ordered Stirling cycle numbers, T(n, k) = k!* s(n, k); n >= 0 k >= 0.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 2, 6, 6, 0, 6, 22, 36, 24, 0, 24, 100, 210, 240, 120, 0, 120, 548, 1350, 2040, 1800, 720, 0, 720, 3528, 9744, 17640, 21000, 15120, 5040, 0, 5040, 26136, 78792, 162456, 235200, 231840, 141120, 40320, 0, 40320, 219168, 708744, 1614816
Offset: 0

Views

Author

Peter Luschny, May 20 2013

Keywords

Comments

The Digital Library of Mathematical Functions defines the Stirling cycle numbers as (-1)^(n-k) times the Stirling numbers of the first kind.

Examples

			[n\k][0,   1,   2,    3,    4,    5,   6]
[0]   1,
[1]   0,   1,
[2]   0,   1,   2,
[3]   0,   2,   6,    6,
[4]   0,   6,  22,   36,   24,
[5]   0,  24, 100,  210,  240,  120,
[6]   0, 120, 548, 1350, 2040, 1800, 720.
...
T(4,2) = 22: The table below shows the compositions of 4 into two parts.
n = 4    Composition       Weight     4!*Weight
            3 + 1            1/3         8
            1 + 3            1/3         8
            2 + 2          1/2*1/2       6
                                        = =
                                  total 22
		

References

  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, table 245.

Crossrefs

Cf. A048594 (signed version without the first column), A132393.

Programs

  • Maple
    A225479 := proc(n, k) option remember;
    if k > n or  k < 0 then return(0) fi;
    if n = 0 and k = 0 then return(1) fi;
    k*A225479(n-1, k-1) + (n-1)*A225479(n-1, k) end;
    for n from 0 to 9 do seq(A225479(n, k), k = 0..n) od;
  • Mathematica
    t[n_, k_] := k!*StirlingS1[n, k] // Abs; Table[t[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 02 2013 *)
  • PARI
    T(n,k)={k!*abs(stirling(n,k,1))} \\ Andrew Howroyd, Jul 27 2020
  • Sage
    def A225479(n, k): return factorial(k)*stirling_number1(n, k)
    for n in (0..6): [A225479(n,k) for k in (0..n)]
    

Formula

For a recursion see the Maple program.
T(n, 0) = A000007; T(n, 1) = A000142; T(n, 2) = A052517.
T(n, 3) = A052748; T(n, n) = A000142; T(n, n-1) = A001286.
row sums = A007840; alternating row sums = A006252.
From Peter Bala, Sep 20 2013: (Start)
E.g.f.: 1/(1 + x*log(1 - t)) = 1 + x*t + (x + 2*x^2)*t^2/2! + (2*x + 6*x^2 + 6*x^3)*t^3/3! + ....
T(n,k) = n!*( the sum of the total weight of the compositions of n into k parts where each part i has weight 1/i ) (see Eger, Theorem 1). An example is given below. (End)
T(n,k) = A132393(n,k) * A000142(k). - Philippe Deléham, Jun 24 2015

A075181 Coefficients of certain polynomials (rising powers).

Original entry on oeis.org

1, 2, 1, 6, 6, 2, 24, 36, 22, 6, 120, 240, 210, 100, 24, 720, 1800, 2040, 1350, 548, 120, 5040, 15120, 21000, 17640, 9744, 3528, 720, 40320, 141120, 231840, 235200, 162456, 78792, 26136, 5040, 362880, 1451520, 2751840, 3265920, 2693880, 1614816
Offset: 1

Views

Author

Wolfdieter Lang, Sep 19 2002

Keywords

Comments

This is the unsigned triangle A048594 with rows read backwards.
The row polynomials p(n,y) := Sum_{m=0..n-1}a(n,m)*y^m, n>=1, are obtained from (log(x)*(-x*log(x))^n)*(d^n/dx^n)(1/log(x)), n>=1, after replacement of log(x) by y.
The gcd of row n is A075182(n). Row sums give A007840(n), n>=1.
The columns give A000142 (factorials), A001286 (Lah), 2* A075183, 2*A075184, 4*A075185, 4!*A075186, 4!*A075187 for m=0..6.
Coefficients T(n,k) of the differential operator expansion
[x^(1+y)D]^n = x^(n*y)[T(n,1)* (xD)^n / n! + y * T(n,2)* (xD)^(n-1) / (n-1)! + ... + y^(n-1) * T(n,n) * (xD)], where D = d/dx. Note that (xD)^n = Bell(n,:xD:), where (:xD:)^n = x^n * D^n and Bell(n,x) are the Bell / Touchard polynomials. See A094638. - Tom Copeland, Aug 22 2015

Examples

			Triangle starts:
1;
2,1;
6,6,2;
24,36,22,6;
...
n=2: (x^2*log(x)^3)*(d^2/d^x^2)(1/log(x)) = 2 + log(x).
		

Crossrefs

Programs

  • Maple
    seq(seq(k!*abs(Stirling1(n,k)),k=n..1,-1),n=1..10); # Robert Israel, Jul 12 2015
  • Mathematica
    Table[ Table[ k!*StirlingS1[n, k] // Abs, {k, 1, n}] // Reverse, {n, 1, 9}] // Flatten (* Jean-François Alcover, Jun 21 2013 *)
  • PARI
    {T(n, k)= if(k<0 || k>=n, 0, (-1)^k* stirling(n, n-k)* (n-k)!)} /* Michael Somos Apr 11 2007 */

Formula

a(n, m) = (n-m)!*|S1(n, n-m)|, n>=m+1>=1, else 0, with S1(n, m) := A008275(n, m) (Stirling1).
a(n, m) = (n-m)*a(n-1, m)+(n-1)*a(n-1, m-1), if n>=m+1>=1, a(n, -1) := 0 and a(1, 0)=1, else 0.

A320080 Square array A(n,k), n >= 0, k >= 0, read by antidiagonals, where column k is the expansion of e.g.f. 1/(1 - k*log(1 + x)).

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 6, 2, 0, 1, 4, 15, 28, 4, 0, 1, 5, 28, 114, 172, 14, 0, 1, 6, 45, 296, 1152, 1328, 38, 0, 1, 7, 66, 610, 4168, 14562, 12272, 216, 0, 1, 8, 91, 1092, 11020, 73376, 220842, 132480, 600, 0, 1, 9, 120, 1778, 24084, 248870, 1550048, 3907656, 1633344, 6240, 0
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 05 2018

Keywords

Examples

			E.g.f. of column k: A_k(x) = 1 + k*x/1! + k*(2*k - 1)*x^2/2! + 2*k*(3*k^2 - 3*k + 1)*x^3/3! + 2*k*(12*k^3 - 18*k^2 + 11*k - 3)*x^4/4! + ...
Square array begins:
  1,   1,     1,      1,      1,       1,  ...
  0,   1,     2,      3,      4,       5,  ...
  0,   1,     6,     15,     28,      45,  ...
  0,   2,    28,    114,    296,     610,  ...
  0,   4,   172,   1152,   4168,   11020,  ...
  0,  14,  1328,  14562,  73376,  248870,  ...
		

Crossrefs

Columns k=0..5 give A000007, A006252, A088501, A335531, A354147, A365604.
Main diagonal gives A317172.

Programs

  • Mathematica
    Table[Function[k, n! SeriesCoefficient[1/(1 - k Log[1 + x]), {x, 0, n}]][j - n], {j, 0, 10}, {n, 0, j}] // Flatten

Formula

E.g.f. of column k: 1/(1 - k*log(1 + x)).
A(n,k) = Sum_{j=0..n} Stirling1(n,j)*j!*k^j.
A(0,k) = 1; A(n,k) = k * Sum_{j=1..n} (-1)^(j-1) * (j-1)! * binomial(n,j) * A(n-j,k). - Seiichi Manyama, May 22 2022
Showing 1-10 of 25 results. Next