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

A001710 Order of alternating group A_n, or number of even permutations of n letters.

Original entry on oeis.org

1, 1, 1, 3, 12, 60, 360, 2520, 20160, 181440, 1814400, 19958400, 239500800, 3113510400, 43589145600, 653837184000, 10461394944000, 177843714048000, 3201186852864000, 60822550204416000, 1216451004088320000, 25545471085854720000, 562000363888803840000
Offset: 0

Views

Author

Keywords

Comments

For n >= 3, a(n-1) is also the number of ways that a 3-cycle in the symmetric group S_n can be written as a product of 2 long cycles (of length n). - Ahmed Fares (ahmedfares(AT)my-deja.com), Aug 14 2001
a(n) is the number of Hamiltonian circuit masks for an n X n adjacency matrix of an undirected graph. - Chad Brewbaker, Jan 31 2003
a(n-1) is the number of necklaces one can make with n distinct beads: n! bead permutations, divide by two to represent flipping the necklace over, divide by n to represent rotating the necklace. Related to Stirling numbers of the first kind, Stirling cycles. - Chad Brewbaker, Jan 31 2003
Number of increasing runs in all permutations of [n-1] (n>=2). Example: a(4)=12 because we have 12 increasing runs in all the permutations of [3] (shown in parentheses): (123), (13)(2), (3)(12), (2)(13), (23)(1), (3)(2)(1). - Emeric Deutsch, Aug 28 2004
Minimum permanent over all n X n (0,1)-matrices with exactly n/2 zeros. - Simone Severini, Oct 15 2004
The number of permutations of 1..n that have 2 following 1 for n >= 1 is 0, 1, 3, 12, 60, 360, 2520, 20160, ... . - Jon Perry, Sep 20 2008
Starting (1, 3, 12, 60, ...) = binomial transform of A000153: (1, 2, 7, 32, 181, ...). - Gary W. Adamson, Dec 25 2008
First column of A092582. - Mats Granvik, Feb 08 2009
The asymptotic expansion of the higher order exponential integral E(x,m=1,n=3) ~ exp(-x)/x*(1 - 3/x + 12/x^2 - 60/x^3 + 360/x^4 - 2520/x^5 + 20160/x^6 - 81440/x^7 + ...) leads to the sequence given above. See A163931 and A130534 for more information. - Johannes W. Meijer, Oct 20 2009
For n>1: a(n) = A173333(n,2). - Reinhard Zumkeller, Feb 19 2010
Starting (1, 3, 12, 60, ...) = eigensequence of triangle A002260, (a triangle with k terms of (1,2,3,...) in each row given k=1,2,3,...). Example: a(6) = 360, generated from (1, 2, 3, 4, 5) dot (1, 1, 3, 12, 60) = (1 + 2 + 9 + 48 + 300). - Gary W. Adamson, Aug 02 2010
For n>=2: a(n) is the number of connected 2-regular labeled graphs on (n+1) nodes (Cf. A001205). - Geoffrey Critzer, Feb 16 2011.
The Fi1 and Fi2 triangle sums of A094638 are given by the terms of this sequence (n>=1). For the definition of these triangle sums see A180662. - Johannes W. Meijer, Apr 20 2011
Also [1, 1] together with the row sums of triangle A162608. - Omar E. Pol, Mar 09 2012
a(n-1) is, for n>=2, also the number of necklaces with n beads (only C_n symmetry, no turnover) with n-1 distinct colors and signature c[.]^2 c[.]^(n-2). This means that two beads have the same color, and for n=2 the second factor is omitted. Say, cyclic(c[1]c[1]c[2]c[3]..c[n-1]), in short 1123...(n-1), taken cyclically. E.g., n=2: 11, n=3: 112, n=4: 1123, 1132, 1213, n=5: 11234, 11243, 11324, 11342, 11423, 11432, 12134, 12143, 13124, 13142, 14123, 14132. See the next-to-last entry in line n>=2 of the representative necklace partition array A212359. - Wolfdieter Lang, Jun 26 2012
For m >= 3, a(m-1) is the number of distinct Hamiltonian circuits in a complete simple graph with m vertices. See also A001286. - Stanislav Sykora, May 10 2014
In factorial base (A007623) these numbers have a simple pattern: 1, 1, 1, 11, 200, 2200, 30000, 330000, 4000000, 44000000, 500000000, 5500000000, 60000000000, 660000000000, 7000000000000, 77000000000000, 800000000000000, 8800000000000000, 90000000000000000, 990000000000000000, etc. See also the formula based on this observation, given below. - Antti Karttunen, Dec 19 2015
Also (by definition) the independence number of the n-transposition graph. - Eric W. Weisstein, May 21 2017
Number of permutations of n letters containing an even number of even cycles. - Michael Somos, Jul 11 2018
Equivalent to Brewbaker's and Sykora's comments, a(n - 1) is the number of undirected cycles covering n labeled vertices, hence the logarithmic transform of A002135. - Gus Wiseman, Oct 20 2018
For n >= 2 and a set of n distinct leaf labels, a(n) is the number of binary, rooted, leaf-labeled tree topologies that have a caterpillar shape (column k=1 of A306364). - Noah A Rosenberg, Feb 11 2019
Also the clique covering number of the n-Bruhat graph. - Eric W. Weisstein, Apr 19 2019
a(n) is the number of lattices of the form [s,w] in the weak order on S_n, for a fixed simple reflection s. - Bridget Tenner, Jan 16 2020
For n > 3, a(n) = p_1^e_1*...*p_m^e_m, where p_1 = 2 and e_m = 1. There exists p_1^x where x <= e_1 such that p_1^x*p_m^e_m is a primitive Zumkeller number (A180332) and p_1^e_1*p_m^e_m is a Zumkeller number (A083207). Therefore, for n > 3, a(n) = p_1^e_1*p_m^e_m*r, where r is relatively prime to p_1*p_m, is also a Zumkeller number. - Ivan N. Ianakiev, Mar 11 2020
For n>1, a(n) is the number of permutations of [n] that have 1 and 2 as cycle-mates, that is, 1 and 2 are contained in the same cycle of a cyclic representation of permutations of [n]. For example, a(4) counts the 12 permutations with 1 and 2 as cycle-mates, namely, (1 2 3 4), (1 2 4 3), (1 3 2 4), (1 3 4 2), (1 4 2 3), (1 4 3 2), (1 2 3) (4), (1 3 2) (4), (1 2 4 )(3), (1 4 2)(3), (1 2)(3 4), and (1 2)(3)(4). Since a(n+2)=row sums of A162608, our result readily follows. - Dennis P. Walsh, May 28 2020

Examples

			G.f. = 1 + x + x^2 + 3*x^3 + 12*x^4 + 60*x^5 + 360*x^6 + 2520*x^7 + ...
		

References

  • J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, pp. 87-8, 20. (a), c_n^e(t=1).
  • 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

a(n+1)= A046089(n, 1), n >= 1 (first column of triangle), A161739 (q(n) sequence).
Bisections are A002674 and A085990 (essentially).
Row 3 of A265609 (essentially).
Row sums of A307429.

Programs

  • Magma
    [1] cat [Order(AlternatingGroup(n)): n in [1..20]]; // Arkadiusz Wesolowski, May 17 2014
    
  • Maple
    seq(mul(k, k=3..n), n=0..20); # Zerinvary Lajos, Sep 14 2007
  • Mathematica
    a[n_]:= If[n > 2, n!/2, 1]; Array[a, 21, 0]
    a[n_]:= If[n<3, 1, n*a[n-1]]; Array[a, 21, 0]; (* Robert G. Wilson v, Apr 16 2011 *)
    a[ n_]:= If[n<0, 0, n! SeriesCoefficient[(2-x^2)/(2-2x), {x, 0, n}]]; (* Michael Somos, May 22 2014 *)
    a[ n_]:= If[n<0, 0, n! SeriesCoefficient[1 +Sinh[-Log[1-x]], {x, 0, n}]]; (* Michael Somos, May 22 2014 *)
    Numerator[Range[0, 20]!/2] (* Eric W. Weisstein, May 21 2017 *)
    Table[GroupOrder[AlternatingGroup[n]], {n, 0, 20}] (* Eric W. Weisstein, May 21 2017 *)
  • PARI
    {a(n) = if( n<2, n>=0, n!/2)};
    
  • PARI
    a(n)=polcoeff(1+x*sum(m=0,n,m^m*x^m/(1+m*x+x*O(x^n))^m),n) \\ Paul D. Hanna
    
  • PARI
    A001710=n->n!\2+(n<2) \\ M. F. Hasler, Dec 01 2013
    
  • Python
    from math import factorial
    def A001710(n): return factorial(n)>>1 if n > 1 else 1 # Chai Wah Wu, Feb 14 2023
    
  • SageMath
    def A001710(n): return (factorial(n) +int(n<2))//2
    [A001710(n) for n in range(31)] # G. C. Greubel, Sep 28 2024
  • Scheme
    ;; Using memoization-macro definec for which an implementation can be found in http://oeis.org/wiki/Memoization
    (definec (A001710 n) (cond ((<= n 2) 1) (else (* n (A001710 (- n 1))))))
    ;; Antti Karttunen, Dec 19 2015
    

Formula

a(n) = numerator(n!/2) and A141044(n) = denominator(n!/2).
D-finite with recurrence: a(0) = a(1) = a(2) = 1; a(n) = n*a(n-1) for n>2. - Chad Brewbaker, Jan 31 2003 [Corrected by N. J. A. Sloane, Jul 25 2008]
a(0) = 0, a(1) = 1; a(n) = Sum_{k=1..n-1} k*a(k). - Amarnath Murthy, Oct 29 2002
Stirling transform of a(n+1) = [1, 3, 12, 160, ...] is A083410(n) = [1, 4, 22, 154, ...]. - Michael Somos, Mar 04 2004
First Eulerian transform of A000027. See A000142 for definition of FET. - Ross La Haye, Feb 14 2005
From Paul Barry, Apr 18 2005: (Start)
a(n) = 0^n + Sum_{k=0..n} (-1)^(n-k-1)*T(n-1, k)*cos(Pi*(n-k-1)/2)^2.
T(n,k) = abs(A008276(n, k)). (End)
E.g.f.: (2 - x^2)/(2 - 2*x).
E.g.f. of a(n+2), n>=0, is 1/(1-x)^3.
E.g.f.: 1 + sinh(log(1/(1-x))). - Geoffrey Critzer, Dec 12 2010
a(n+1) = (-1)^n * A136656(n,1), n>=1.
a(n) = n!/2 for n>=2 (proof from the e.g.f). - Wolfdieter Lang, Apr 30 2010
a(n) = (n-2)! * t(n-1), n>1, where t(n) is the n-th triangular number (A000217). - Gary Detlefs, May 21 2010
a(n) = ( A000254(n) - 2* A001711(n-3) )/3, n>2. - Gary Detlefs, May 24 2010
O.g.f.: 1 + x*Sum_{n>=0} n^n*x^n/(1 + n*x)^n. - Paul D. Hanna, Sep 13 2011
a(n) = if n < 2 then 1, otherwise Pochhammer(n,n)/binomial(2*n,n). - Peter Luschny, Nov 07 2011
a(n) = Sum_{k=0..floor(n/2)} s(n,n-2*k) where s(n,k) are Stirling number of the first kind, A048994. - Mircea Merca, Apr 07 2012
a(n-1), n>=3, is M_1([2,1^(n-2)])/n = (n-1)!/2, with the M_1 multinomial numbers for the given n-1 part partition of n. See the second to last entry in line n>=3 of A036038, and the above necklace comment by W. Lang. - Wolfdieter Lang, Jun 26 2012
G.f.: A(x) = 1 + x + x^2/(G(0)-2*x) where G(k) = 1 - (k+1)*x/(1 - x*(k+3)/G(k+1)); (continued fraction). - Sergei N. Gladkovskii, Dec 26 2012.
G.f.: 1 + x + (Q(0)-1)*x^2/(2*(sqrt(x)+x)), where Q(k) = 1 + (k+2)*sqrt(x)/(1 - sqrt(x)/(sqrt(x) + 1/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 15 2013
G.f.: 1 + x + (x*Q(x)-x^2)/(2*(sqrt(x)+x)), where Q(x) = Sum_{n>=0} (n+1)!*x^n*sqrt(x)*(sqrt(x) + x*(n+2)). - Sergei N. Gladkovskii, May 15 2013
G.f.: 1 + x/2 + (Q(0)-1)*x/(2*(sqrt(x)+x)), where Q(k) = 1 + (k+1)*sqrt(x)/(1 - sqrt(x)/(sqrt(x) + 1/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 15 2013
G.f.: 1 + x + x^2*G(0)/2, where G(k) = 1 + 1/(1 - x/(x + 1/(k+3)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 01 2013
G.f.: 1+x + x^2*W(0), where W(k) = 1 - x*(k+3)/( x*(k+3) - 1/(1 - x*(k+1)/( x*(k+1) - 1/W(k+1) ))); (continued fraction). - Sergei N. Gladkovskii, Aug 26 2013
From Antti Karttunen, Dec 19 2015: (Start)
a(0)=a(1)=1; after which, for even n: a(n) = (n/2) * (n-1)!, and for odd n: a(n) = (n-1)/2 * ((n-1)! + (n-2)!). [The formula was empirically found after viewing these numbers in factorial base, A007623, and is easily proved by considering formulas from Lang (Apr 30 2010) and Detlefs (May 21 2010) shown above.]
For n >= 1, a(2*n+1) = a(2*n) + A153880(a(2*n)). [Follows from above.] (End)
Inverse Stirling transform of a(n) is (-1)^(n-1)*A009566(n). - Anton Zakharov, Aug 07 2016
a(n) ~ sqrt(Pi/2)*n^(n+1/2)/exp(n). - Ilya Gutkovskiy, Aug 07 2016
a(n) = A006595(n-1)*n/A000124(n) for n>=2. - Anton Zakharov, Aug 23 2016
a(n) = A001563(n-1) - A001286(n-1) for n>=2. - Anton Zakharov, Sep 23 2016
From Peter Bala, May 24 2017: (Start)
The o.g.f. A(x) satisfies the Riccati equation x^2*A'(x) + (x - 1)*A(x) + 1 - x^2 = 0.
G.f.: A(x) = 1 + x + x^2/(1 - 3*x/(1 - x/(1 - 4*x/(1 - 2*x/(1 - 5*x/(1 - 3*x/(1 - ... - (n + 2)*x/(1 - n*x/(1 - ... ))))))))) (apply Stokes, 1982).
A(x) = 1 + x + x^2/(1 - 2*x - x/(1 - 3*x/(1 - 2*x/(1 - 4*x/(1 - 3*x/(1 - 5*x/(1 - ... - n*x/(1 - (n+2)*x/(1 - ... ))))))))). (End)
H(x) = (1 - (1 + x)^(-2)) / 2 = x - 3*x^2/2! + 12*x^3/3! - ..., an e.g.f. for the signed sequence here (n!/2!), ignoring the first two terms, is the compositional inverse of G(x) = (1 - 2*x)^(-1/2) - 1 = x + 3*x^2/2! + 15*x^3/3! + ..., an e.g.f. for A001147. Cf. A094638. H(x) is the e.g.f. for the sequence (-1)^m * m!/2 for m = 2,3,4,... . Cf. A001715 for n!/3! and A001720 for n!/4!. Cf. columns of A094587, A173333, and A213936 and rows of A138533. - Tom Copeland, Dec 27 2019
From Amiram Eldar, Jan 08 2023: (Start)
Sum_{n>=0} 1/a(n) = 2*(e-1).
Sum_{n>=0} (-1)^n/a(n) = 2/e. (End)

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Aug 20 2001
Further terms from Simone Severini, Oct 15 2004

A143494 Triangle read by rows: 2-Stirling numbers of the second kind.

Original entry on oeis.org

1, 2, 1, 4, 5, 1, 8, 19, 9, 1, 16, 65, 55, 14, 1, 32, 211, 285, 125, 20, 1, 64, 665, 1351, 910, 245, 27, 1, 128, 2059, 6069, 5901, 2380, 434, 35, 1, 256, 6305, 26335, 35574, 20181, 5418, 714, 44, 1, 512, 19171, 111645, 204205, 156660, 58107, 11130, 1110, 54, 1
Offset: 2

Views

Author

Peter Bala, Aug 20 2008

Keywords

Comments

This is the case r = 2 of the r-Stirling numbers of the second kind. The 2-Stirling numbers of the second kind give the number of ways of partitioning the set {1,2,...,n} into k nonempty disjoint subsets with the restriction that the elements 1 and 2 belong to distinct subsets.
More generally, the r-Stirling numbers of the second kind give the number of ways of partitioning the set {1,2,...,n} into k nonempty disjoint subsets with the restriction that the numbers 1, 2, ..., r belong to distinct subsets. The case r = 1 gives the usual Stirling numbers of the second kind A008277; for other cases see A143495 (r = 3) and A143496 (r = 4).
The lower unitriangular array of r-Stirling numbers of the second kind equals the matrix product P^(r-1) * S (with suitable offsets in the row and column indexing), where P is Pascal's triangle, A007318 and S is the array of Stirling numbers of the second kind, A008277.
For the definition of and entries relating to the corresponding r-Stirling numbers of the first kind see A143491. For entries on r-Lah numbers refer to A143497. The theory of r-Stirling numbers of both kinds is developed in [Broder].
From Peter Bala, Sep 19 2008: (Start)
Let D be the derivative operator d/dx and E the Euler operator x*d/dx. Then x^(-2)*E^n*x^2 = Sum_{k = 0..n} T(n+2,k+2)*x^k*D^k.
The row generating polynomials R_n(x) := Sum_{k= 2..n} T(n,k)*x^k satisfy the recurrence R_(n+1)(x) = x*R_n(x) + x*d/dx(R_n(x)) with R_2(x) = x^2. It follows that the polynomials R_n(x) have only real zeros (apply Corollary 1.2. of [Liu and Wang]).
Relation with the 2-Eulerian numbers E_2(n,j) := A144696(n,j): T(n,k) = 2!/k!*Sum_ {j = n-k..n-2} E_2(n,j)*binomial(j,n-k) for n >= k >= 2. (End)
From Wolfdieter Lang, Sep 29 2011: (Start)
T(n,k) = S(n,k,2), n>=k>=2, in Mikhailov's first paper, eq.(28) or (A3). E.g.f. column no. k from (A20) with k->2, r->k. Therefore, with offset [0,0], this triangle is the Sheffer triangle (exp(2*x),exp(x)-1) with e.g.f. of column no. m>=0: exp(2*x)*((exp(x)-1)^m)/m!. See one of the formulas given below. For Sheffer matrices see the W. Lang link under A006232 with the S. Roman reference, also found in A132393. (End)

Examples

			Triangle begins
  n\k|...2....3....4....5....6....7
  =================================
  2..|...1
  3..|...2....1
  4..|...4....5....1
  5..|...8...19....9....1
  6..|..16...65...55...14....1
  7..|..32..211..285..125...20....1
  ...
T(4,3) = 5. The set {1,2,3,4} can be partitioned into three subsets such that 1 and 2 belong to different subsets in 5 ways: {{1}{2}{3,4}}, {{1}{3}{2,4}}, {{1}{4}{2,3}}, {{2}{3}{1,4}} and {{2}{4}{1,3}}; the remaining possibility {{1,2}{3}{4}} is not allowed.
From _Peter Bala_, Feb 23 2025: (Start)
The array factorizes as
/ 1               \       /1             \ /1             \ /1            \
| 2    1           |     | 2   1          ||0  1           ||0  1          |
| 4    5   1       |  =  | 4   3   1      ||0  2   1       ||0  0  1       | ...
| 8   19   9   1   |     | 8   7   4  1   ||0  4   3  1    ||0  0  2  1    |
|16   65  55  14  1|     |16  15  11  6  1||0  8   7  4  1 ||0  0  4  3  1 |
|...               |     |...             ||...            ||...           |
where, in the infinite product on the right-hand side, the first array is the Riordan array (1/(1 - 2*x), x/(1 - x)). See A055248. (End)
		

Crossrefs

A001047 (column 3), A005493 (row sums), A008277, A016269 (column 4), A025211 (column 5), A049444 (matrix inverse), A074051 (alt. row sums).

Programs

  • Maple
    with combinat: T := (n, k) -> (1/(k-2)!)*add ((-1)^(k-i)*binomial(k-2,i)*(i+2)^(n-2),i = 0..k-2): for n from 2 to 11 do seq(T(n, k), k = 2..n) end do;
  • Mathematica
    t[n_, k_] := StirlingS2[n, k] - StirlingS2[n-1, k]; Flatten[ Table[ t[n, k], {n, 2, 11}, {k, 2, n}]] (* Jean-François Alcover, Dec 02 2011 *)
  • Sage
    @CachedFunction
    def stirling2r(n, k, r) :
        if n < r: return 0
        if n == r: return 1 if k == r else 0
        return stirling2r(n-1,k-1,r) + k*stirling2r(n-1,k,r)
    A143494 = lambda n,k: stirling2r(n, k, 2)
    for n in (2..6):
        [A143494(n, k) for k in (2..n)] # Peter Luschny, Nov 19 2012

Formula

T(n+2,k+2) = (1/k!)*Sum_{i = 0..k} (-1)^(k-i)*C(k,i)*(i+2)^n, n,k >= 0.
T(n,k) = Stirling2(n,k) - Stirling2(n-1,k) for n, k >= 2.
Recurrence relation: T(n,k) = T(n-1,k-1) + k*T(n-1,k) for n > 2, with boundary conditions T(n,1) = T(1,n) = 0 for all n, T(2,2) = 1 and T(2,k) = 0 for k > 2. Special cases: T(n,2) = 2^(n-2); T(n,3) = 3^(n-2) - 2^(n-2).
As a sum of monomial functions of degree m: T(n+m,n) = Sum_{2 <= i_1 <= ... <= i_m <= n} (i_1*i_2*...*i_m). For example, T(6,4) = Sum_{2 <= i <= j <= 4} (i*j) = 2*2 + 2*3 + 2*4 + 3*3 + 3*4 + 4*4 = 55.
E.g.f. column k+2 (with offset 2): 1/k!*exp(2*x)*(exp(x) - 1)^k.
O.g.f. k-th column: Sum_{n >= k} T(n,k)*x^n = x^k/((1-2*x)*(1-3*x)*...*(1-k*x)).
E.g.f.: exp(2*t + x*(exp(t) - 1)) = Sum_{n >= 0} Sum_{k = 0..n} T(n+2,k+2) *x^k*t^n/n! = Sum_{n >= 0} B_n(2;x)*t^n/n! = 1 + (2 + x)*t/1! + (4 + 5*x + x^2)*t^2/2! + ..., where the row polynomial B_n(2;x) := Sum_{k = 0..n} T(n+2,k+2)*x^k denotes the 2-Bell polynomial.
Dobinski-type identities: Row polynomial B_n(2;x) = exp(-x)*Sum_{i >= 0} (i + 2)^n*x^i/i!. Sum_{k = 0..n} k!*T(n+2,k+2)*x^k = Sum_{i >= 0} (i + 2)^n*x^i/(1 + x)^(i+1).
The T(n,k) are the connection coefficients between falling factorials and the shifted monomials (x + 2)^(n-2). For example, from row 4 we have 4 + 5*x + x*(x - 1) = (x + 2)^2, while from row 5 we have 8 + 19*x + 9*x*(x - 1) + x*(x - 1)*(x - 2) = (x + 2)^3.
The row sums of the array are the 2-Bell numbers, B_n(2;1), equal to A005493(n-2). The alternating row sums are the complementary 2-Bell numbers, B_n(2;-1), equal to (-1)^n*A074051(n-2).
This array is the matrix product P * S, where P denotes the Pascal triangle, A007318 and S denotes the lower triangular array of Stirling numbers of the second kind, A008277 (apply Theorem 10 of [Neuwirth]).
Also, this array equals the transpose of the upper triangular array A126351. The inverse array is A049444, the signed 2-Stirling numbers of the first kind. See A143491 for the unsigned version of the inverse.
Let f(x) = exp(exp(x)). Then for n >= 1, the row polynomials R(n,x) are given by R(n+2,exp(x)) = 1/f(x)*(d/dx)^n(exp(2*x)*f(x)). Similar formulas hold for A008277, A039755, A105794, A111577 and A154537. - Peter Bala, Mar 01 2012

A143491 Unsigned 2-Stirling numbers of the first kind.

Original entry on oeis.org

1, 2, 1, 6, 5, 1, 24, 26, 9, 1, 120, 154, 71, 14, 1, 720, 1044, 580, 155, 20, 1, 5040, 8028, 5104, 1665, 295, 27, 1, 40320, 69264, 48860, 18424, 4025, 511, 35, 1, 362880, 663696, 509004, 214676, 54649, 8624, 826, 44, 1, 3628800, 6999840, 5753736, 2655764
Offset: 2

Views

Author

Peter Bala, Aug 20 2008

Keywords

Comments

Essentially the same as A136124 but with column numbers differing by one. See A049444 for a signed version of this array. The unsigned 2-Stirling numbers of the first kind count the permutations of the set {1,2,...,n} into k disjoint cycles, with the restriction that the elements 1 and 2 belong to distinct cycles. This is the particular case r = 2 of the unsigned r-Stirling numbers of the first kind, which count the permutations of the set {1,2,...,n} into k disjoint cycles, with the restriction that the numbers 1, 2, ..., r belong to distinct cycles. The case r = 1 gives the usual unsigned Stirling numbers of the first kind, abs(A008275); for other cases see A143492 (r = 3) and A143493 (r = 4). The corresponding 2-Stirling numbers of the second kind can be found in A143494.
In general, the lower unitriangular array of unsigned r-Stirling numbers of the first kind (with suitable offsets in the row and column indexing) equals the matrix product St1 * P^(r-1), where St1 is the array of unsigned Stirling numbers of the first kind, abs(A008275) and P is Pascal's triangle, A007318. The theory of r-Stirling numbers of both kinds is developed in [Broder]. For details of the related r-Lah numbers see A143497.
This sequence also represents the number of permutations in the alternating group An of length k, where the length is taken with respect to the generators set {(12)(ij)}. For a bijective proof of the relation between these numbers and the 2-Stirling numbers of the first kind see the Rotbart link. - Aviv Rotbart, May 05 2011
With offset n=0,k=0 : triangle T(n,k), read by rows, given by [2,1,3,2,4,3,5,4,6,5,...] DELTA [1,0,1,0,1,0,1,0,1,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Sep 29 2011
With offset n=0 and k=0, this is the Sheffer triangle (1/(1-x)^2,-log(1-x)) (in the umbral notation of S. Roman's book this would be called Sheffer for (exp(-2*t),1-exp(-t))). See the e.g.f. given below. Compare also with the e.g.f. for the signed version A049444. - Wolfdieter Lang, Oct 10 2011
Reversed rows correspond to the Betti numbers of the moduli space M(0,n) of smooth Riemann surfaces (see Murri link). - Tom Copeland, Sep 19 2012

Examples

			Triangle begins
  n\k|.....2.....3.....4.....5.....6.....7
  ========================================
  2..|.....1
  3..|.....2.....1
  4..|.....6.....5.....1
  5..|....24....26.....9.....1
  6..|...120...154....71....14.....1
  7..|...720..1044...580...155....20.....1
  ...
T(4,3) = 5. The permutations of {1,2,3,4} with 3 cycles such that 1 and 2 belong to different cycles are: (1)(2)(3 4), (1)(3)(24), (1)(4)(23), (2)(3)(14) and (2)(4)(13). The remaining possibility (3)(4)(12) is not allowed.
From _Aviv Rotbart_, May 05 2011: (Start)
Example of the alternating group permutations numbers:
Triangle begins
  n\k|.....0.....1.....2.....3.....4.....5.....6.....7
  ====================================================
  2..|.....1
  3..|.....1.....2
  4..|.....1.....5.....6
  5..|.....1.....9....26....24
  6..|.....1....14....71...154...120
  7..|.....1....20...155...580..1044..720
A(n,k) = number of permutations in An of length k, with respect to the generators set {(12)(ij)}. For example, A(2,0)=1 (only the identity is there), for A4, the generators are {(12)(13),(12)(14),(12,23),(12)(24),(12)(34)}, thus we have A(4,1)=5 (exactly 5 generators), the permutations of length 2 are:
   (12)(13)(12)(13) = (312)
   (12)(13)(12)(14) = (41)(23)
   (12)(13)(12)(24) = (432)(1)
   (12)(13)(12)(34) = (342)(1)
   (12)(23)(12)(24) = (13)(24)
   (12)(14)(12)(14) = (412)(3)
Namely, A(4,2)=6. Together with the identity [=(12)(12), of length 0. therefore A(4,0)=1] we have 12 permutations, comprising all A4 (4!/2=12). (End)
		

Crossrefs

Cf. A001705 - A001709 (column 3..7), A001710 (row sums), A008275, A049444 (signed version), A136124, A143492, A143493, A143494, A143497.
Cf. A094638.

Programs

  • Maple
    with combinat: T := (n, k) -> (n-2)! * add((n-j-1)*abs(stirling1(j,k-2))/j!,j = k-2..n-2): for n from 2 to 10 do seq(T(n, k), k = 2..n) end do;
  • Mathematica
    t[n_, k_] := (n-2)!*Sum[(n-j-1)*Abs[StirlingS1[j, k-2]]/j!, {j, k-2, n-2}]; Table[t[n, k], {n, 2, 11}, {k, 2, n}] // Flatten (* Jean-François Alcover, Apr 16 2013, after Maple *)

Formula

T(n,k) = (n-2)! * Sum_{j = k-2 .. n-2} (n-j-1)*|stirling1(j,k-2)|/j!.
Recurrence relation: T(n,k) = T(n-1,k-1) + (n-1)*T(n-1,k) for n > 2, with boundary conditions: T(n,1) = T(1,n) = 0, for all n; T(2,2) = 1; T(2,k) = 0 for k > 2.
Special cases: T(n,2) = (n-1)!; T(n,3) = (n-1)!*(1/2 + 1/3 + ... + 1/(n-1)).
T(n,k) = Sum_{2 <= i_1 < ... < i_(n-k) < n} (i_1*i_2*...*i_(n-k)). For example, T(6,4) = Sum_{2 <= i < j < 6} (i*j) = 2*3 + 2*4 + 2*5 + 3*4 + 3*5 + 4*5 = 71.
Row g.f.: Sum_{k = 2..n} T(n,k)*x^k = x^2*(x+2)*(x+3)*...*(x+n-1).
E.g.f. for column (k+2): Sum_{n>=k} T(n+2,k+2)*x^n/n! = (1/k!)*(1/(1-x)^2)*(log(1/(1-x)))^k.
E.g.f.: (1/(1-t))^(x+2) = Sum_{n>=0} Sum_{k = 0..n} T(n+2,k+2)*x^k*t^n/n! = 1 + (2+x)*t/1! + (6+5*x+x^2)*t^2/2! + ... .
This array is the matrix product St1 * P, where St1 denotes the lower triangular array of unsigned Stirling numbers of the first kind, abs(A008275) and P denotes Pascal's triangle, A007318. The row sums are n!/2 ( A001710 ). The alternating row sums are (n-2)!.
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-a - j), then T(n-1,i) = |f(n,i,2)|, for n=1,2,...; i=0..n. - Milan Janjic, Dec 21 2008
From Gary W. Adamson, Jul 19 2011: (Start)
n-th row of the triangle = top row of M^(n-2), M = a reversed variant of the (1,2) Pascal triangle (Cf. A029635); as follows:
2, 1, 0, 0, 0, 0, ...
2, 3, 1, 0, 0, 0, ...
2, 5, 4, 1, 0, 0, ...
2, 7, 9, 5, 1, 0, ...
... (End)
The reversed, row polynomials of this entry multiplied by (1+x) are the row polynomials of A094638. E.g., (1+x)(1+5x+6x^2) = (1+6x+11x^2+6x^3). - Tom Copeland, Dec 11 2016

A137452 Triangular array of the coefficients of the sequence of Abel polynomials A(n,x) := x*(x-n)^(n-1).

Original entry on oeis.org

1, 0, 1, 0, -2, 1, 0, 9, -6, 1, 0, -64, 48, -12, 1, 0, 625, -500, 150, -20, 1, 0, -7776, 6480, -2160, 360, -30, 1, 0, 117649, -100842, 36015, -6860, 735, -42, 1, 0, -2097152, 1835008, -688128, 143360, -17920, 1344, -56, 1, 0, 43046721, -38263752, 14880348, -3306744, 459270, -40824, 2268, -72, 1
Offset: 0

Views

Author

Roger L. Bagula, Apr 18 2008

Keywords

Comments

Row sums give A177885.
The Abel polynomials are associated with the Abel operator t*exp(y*t)*p(x) = t*p(x+y).
From Peter Luschny, Jan 14 2009: (Start)
Abs(T(n,k)) is the number of rooted labeled trees on n+1 vertices with a root degree k (Clarke's formula).
The row sums in the unsigned case, Sum_{k=0..n} abs(T(n,k)), count the trees on n+1 labeled nodes, A000272(n+1). (End)
Exponential Riordan array [1, W(x)], W(x) the Lambert W-function. - Paul Barry, Nov 19 2010
The inverse array is the exponential Riordan array [1, x*exp(x)], which is A059297. - Peter Bala, Apr 08 2013
The inverse Bell transform of [1,2,3,...]. See A264428 for the Bell transform and A264429 for the inverse Bell transform. - Peter Luschny, Dec 20 2015
Also the Bell transform of (-1)^n*(n+1)^n. - Peter Luschny, Jan 18 2016

Examples

			Triangle begins:
  1;
  0,        1;
  0,       -2,       1;
  0,        9,      -6,       1;
  0,      -64,      48,     -12,      1;
  0,      625,    -500,     150,    -20,      1;
  0,    -7776,    6480,   -2160,    360,    -30,    1;
  0,   117649, -100842,   36015,  -6860,    735,  -42,   1;
  0, -2097152, 1835008, -688128, 143360, -17920, 1344, -56, 1;
		

References

  • Steve Roman, The Umbral Calculus, Dover Publications, New York (1984), pp. 14 and 29

Crossrefs

Row sums A177885.
Cf. A000272, A061356, A059297 (inverse array), A264429.

Programs

  • Maple
    T := proc(n,k) if n = 0 and k = 0 then 1 else binomial(n-1,k-1)*(-n)^(n-k) fi end; seq(print(seq(T(n,k),k=0..n)),n=0..7); # Peter Luschny, Jan 14 2009
    # The function BellMatrix is defined in A264428.
    BellMatrix(n -> (-n-1)^n, 9); # Peter Luschny, Jan 27 2016
  • Mathematica
    a0 = 1 a[x, 0] = 1; a[x, 1] = x; a[x_, n_] := x*(x - a0*n)^(n - 1); Table[Expand[a[x, n]], {n, 0, 10}]; a1 = Table[CoefficientList[a[x, n], x], {n, 0, 10}]; Flatten[a1]
    (* Second program: *)
    BellMatrix[f_, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len - 1}, {k, 0, len - 1}]];
    B = BellMatrix[Function[n, (-n-1)^n], rows = 12];
    Table[B[[n, k]], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 28 2018, after Peter Luschny *)
  • Sage
    # uses[inverse_bell_transform from A264429]
    def A137452_matrix(dim):
        nat = [n for n in (1..dim)]
        return inverse_bell_transform(dim, nat)
    A137452_matrix(10) # Peter Luschny, Dec 20 2015

Formula

Row n gives the coefficients of the expansion of x*(x-n)^(n-1).
Abs(T(n,k)) = C(n-1,k-1)*n^(n-k). - Peter Luschny, Jan 14 2009
From Wolfdieter Lang, Nov 08 2022: (Start)
From the exponential Riordan (also Sheffer of Jabotinsky) type (1, LambertW) array (see comments).
E.g.f. of column sequence k, LambertW(x)^k/k!, for k >= 0.
E.g.f. of row polynomials P_n(y) = Sum_{k=0..n} T(n, k)*y^k: exp(y*LambertW(x)).
Recurrence for T: T(n, k) = 0 for n < k; T(n, 0) = 1 for n = 0 otherwise 0; T(n, k) = (n/k)*Sum_{j=0..n-k} binomial(k-1+j,k-1)*(-1)^j*T(n-1, k-1+j). (Jabotinsky type convolution triangle, the e.g.f.s for the a- and z-sequences are exp(-x), and 0. See the link in A006232.)
Recurrence for column k of T: T(n, k) = 0 for n < k, T(k, k) = 1, for k >= 0 otherwise T(n, k) = (n!*k/(n-k))*Sum_{j=k..n-1} (1/j!)*beta(n-1-j)*T(j, k), where beta(n) = A264234(n+1)/A095996(n+1) = {-1, 2, -9/2, 32/3, -625/24, ...} with o.g.f. d/dx(log(LambertW(x)/x)). See the Boas-Buck or Rainville references given in A046521, and my Aug 10 2017 comment there.
Recurrence for the row polynomials P_0(x) = 1, and P_n(x) = x*substitute(z=d/dx, exp(-z)/(1+z)) P_(n-1)(x), for n >= 1, with coefficient z^k of exp(-z)/(1+z) given by (-1)^k*A061354(k)/A061355(k). See the Roman reference Corollary 3.7.2., p. 50. (End)
The column sequences for the unsigned triangle Abs(T(n, k)), for k >= 2, are also given by {n^(n-k)*(n-1)*s(k-2, n)/(k-1)!}A049444.%20-%20_Wolfdieter%20Lang">{n>=k} with the row polynomials s(n, x) = risingfactorial(x - (n+1), n) of A049444. - _Wolfdieter Lang, Nov 21 2022

Extensions

Better name by Peter Bala, Apr 08 2013
Edited by Joerg Arndt, Apr 08 2013

A102661 Triangle of partial sums of Stirling numbers of 2nd kind (A008277): T(n,k) = Sum_{i=1..k} Stirling2(n,i), 1<=k<=n.

Original entry on oeis.org

1, 1, 2, 1, 4, 5, 1, 8, 14, 15, 1, 16, 41, 51, 52, 1, 32, 122, 187, 202, 203, 1, 64, 365, 715, 855, 876, 877, 1, 128, 1094, 2795, 3845, 4111, 4139, 4140, 1, 256, 3281, 11051, 18002, 20648, 21110, 21146, 21147, 1, 512, 9842, 43947, 86472, 109299, 115179, 115929, 115974, 115975
Offset: 1

Views

Author

Vladeta Jovovic, Feb 03 2005

Keywords

Comments

T(n,k) is the number of ways to place n distinguishable balls into k indistinguishable bins. - Geoffrey Critzer, Mar 22 2011
From Mark Wildon, Aug 10 2015: (Start)
T(n,k) is the number of partitions of a set of size n into at most k parts.
T(n,k) is the number of sequences of n top-to-random shuffles of a deck of k cards that leave the deck invariant.
T(n,k) = where pi is the natural permutation character of the symmetric group Sym_k. This gives another combinatorial interpretation of T(n,k) as counting sequences of box moves on Young diagrams. Reference linked to below. (End)
Diagonal entries T(n,n) are the Bell numbers A000110. - Robert Israel, Aug 10 2015
From Manfred Boergens, Mar 18 2025: (Start)
The partitions in the second comment can be described as disjoint collections of subsets of [n] without the empty set with union = [n]. For instance, T(4,2)=8 is the number of partitions of [4] into 1 or 2 parts: 1234, 1 234, 2 134, 3 124, 4 123, 12 34, 13 24, 14 23.
For disjoint collections which may include one empty set see A381682.
For arbitrary collections without the empty set see A369950.
For arbitrary collections which may include one empty set see A381683. (End)

Examples

			Triangle begins:
  1;
  1,  2;
  1,  4,  5;
  1,  8, 14, 15;
  1, 16, 41, 51, 52;
  ...
		

References

  • Richard Stanley, Enumerative Combinatorics, Cambridge Univ. Press, 1997 page 38. (#7 of the twelvefold ways)

Crossrefs

Programs

  • Haskell
    a102661 n k = a102661_tabl !! (n-1) !! (k-1)
    a102661_row n = a102661_tabl !! (n-1)
    a102661_tabl = map (scanl1 (+) . tail) $ tail a048993_tabl
    -- Reinhard Zumkeller, Jun 19 2015
    
  • Maple
    with(combinat): A102661_row := proc(n) local k,j; seq(add(stirling2(n,j),j=1..k),k=1..n) end:
    seq(print(A102661_row(r)),r=1..6); # Peter Luschny, Sep 30 2011
  • Mathematica
    Table[Table[Sum[StirlingS2[n, i], {i, 1, k}], {k, 1, n}], {n, 1,10}] // Grid (* Geoffrey Critzer, Mar 22 2011*)
    Table[Accumulate[StirlingS2[n,Range[n]]],{n,10}]//Flatten (* Harvey P. Dale, Oct 28 2019 *)
  • PARI
    tabl(nn) = {for (n=1, nn, for (k=1, n, print1(sum(i=1, k, stirling(n,i, 2)), ", ");); print(););} \\ Michel Marcus, Aug 10 2015
    
  • Sage
    def T(n,k):
        return sum([stirling_number2(n,j) for j in range(1,k+1)])
    # Danny Rorabaugh, Oct 13 2015

Formula

E.g.f. for row polynomials s(n,y) = Sum_{k=0..n} a(n,k)*y^k is (y*e^(e^(x*y)-1)- e^(y*(e^x-1)))/(y-1) - 1. - Robert Israel, Aug 10 2015

A094645 Triangle of generalized Stirling numbers of the first kind.

Original entry on oeis.org

1, -1, 1, 0, -1, 1, 0, -1, 0, 1, 0, -2, -1, 2, 1, 0, -6, -5, 5, 5, 1, 0, -24, -26, 15, 25, 9, 1, 0, -120, -154, 49, 140, 70, 14, 1, 0, -720, -1044, 140, 889, 560, 154, 20, 1, 0, -5040, -8028, -64, 6363, 4809, 1638, 294, 27, 1, 0, -40320, -69264, -8540, 50840, 44835, 17913, 3990, 510, 35, 1
Offset: 0

Views

Author

Vladeta Jovovic, May 17 2004

Keywords

Comments

From Wolfdieter Lang, Jun 20 2011: (Start)
The row polynomials s(n,x) := Sum_{j=0..n} T(n,k)*x^k satisfy risefac(x-1,n) = s(n,x), with the rising factorials risefac(x-1,n) := Product_{j=0..n-1} (x-1+j), n >= 1, risefac(x-1,0) = 1. Compare with the formula risefac(x,n) = s1(n,x), with the row polynomials s1(n,x) of A132393 (unsigned Stirling1).
This is the lower triangular Sheffer array with e.g.f.
T(x,z) = (1-z)*exp(-x*log(1-z)) (the rewritten e.g.f. from the formula section). See the W. Lang link under A006232 for Sheffer matrices and the Roman reference. In the notation which indicates the column e.g.f.s this is Sheffer (1-z,-log(1-z)). In the umbral notation (cf. Roman) this is called Sheffer for (exp(t),1-exp(-t)).
The row polynomials satisfy s(n,x) = (x+n-1)*s(n-1,x), s(0,x)=1, and s(n,x) = (x-1)*s1(n-1,x), n >= 1, s1(0,x) = 1, with the unsigned Stirling1 row polynomials s1(n,x).
The row polynomials also satisfy
s(n,x) - s(n,x-1) = n*s(n-1,x), n > 1, s(0,x) = 1
(from the Meixner identity, see the Meixner reference given at A060338).
The row polynomials satisfy as well (from corollary 3.7.2. p. 50 of the Roman reference)
s(n,x) = (x-1)*s(n-1,x+1), n >= 1, s(0,n) = 1.
The exponential convolution identity is
s(n,x+y) = Sum_{k=0..n} binomial(n,k)*s(k,y)*s1(n-k,x),
n >= 0, with symmetry x <-> y.
The row sums are 1 for n=0 and 0 otherwise, and the alternating row sums are 1,-2,2, followed by zeros, with e.g.f. (1-x)^2.
The Sheffer a-sequence Sha(n) = A164555(n)/A027642(n) with e.g.f. x/(1-exp(-x)), and the z-sequence is Shz(n) = -1 with e.g.f. -exp(x).
The inverse Sheffer matrix is ((-1)^(n-k))*A105794(n,k) with e.g.f. exp(z)*exp(x*(1-exp(-z))). (End)
Triangle T(n,k), read by rows, given by (-1, 1, 0, 2, 1, 3, 2, 4, 3, 5, ...) DELTA (1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Jan 16 2012
Also coefficients of t in t*(t-1)*Sum[(-1)^(n+m) t^(m-1) StirlingS1[n,m], {m,n}] in which setting t^k equal to k gives n!, from this follows that the dot product of row n with [0,...,n] equals (n-1)!. - Wouter Meeussen, May 15 2012

Examples

			Triangle begins
   1;
  -1,   1;
   0,  -1,   1;
   0,  -1,   0,   1;
   0,  -2,  -1,   2,   1;
   0,  -6,  -5,   5,   5,   1;
   0, -24, -26,  15,  25,   9,   1;
   ...
Recurrence:
  -2 = T(4,1) = T(3,0) + (4-2)*T(3,1) = 0 + 2*(-1).
Row polynomials:
  s(3,x) = -x+x^3 = (x-1)*s1(2,x) = (x-1)*(x+x^2).
  s(3,x) = (x-1)*s(2,x+1) = (x-1)*(-(x+1)+(x+1)^2).
  s(3,x) - s(3,x-1) = -x+x^3 -(-(x-1)+(x-1)^3) = 3*(-x+x^2) = 3*s(2,x).
		

References

  • S. Roman, The Umbral Calculus, Academic Press, New York, 1984.

Crossrefs

Programs

Formula

E.g.f.: (1-y)^(1-x).
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000142(n), A000142(n+1), A001710(n+2), A001715(n+3), A001720(n+4), A001725(n+5), A001730(n+6), A049388(n), A049389(n), A049398(n), A051431(n) for x = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 respectively. - Philippe Deléham, Nov 13 2007
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-a-j), then |T(n,i)| = |f(n,i,-1)|, for n=1,2,...; i=0..n. - Milan Janjic, Dec 21 2008
From Wolfdieter Lang, Jun 20 2011: (Start)
T(n,k) = |S1(n-1,k-1)| - |S1(n-1,k)|, n >= 1, k >= 1, with |S1(n,k)| = A132393(n,k) (unsigned Stirling1).
Recurrence: T(n,k) = T(n-1,k-1) + (n-2)*T(n-1,k) if n >= k >= 0; T(n,k) = 0 if n < k; T(n,-1) = 0; T(0,0) = 1.
E.g.f. column k: (1-x)*((-log(1-x))^k)/k!. (End)
T(n,k) = Sum_{i=0..n} binomial(n,i)*(n-i)!*Stirling1(i,k)*TC(m,n,i) where TC(m,n,k) = Sum_{i=0..n-k} binomial(n+1,n-k-i)*Stirling2(i+m+1,i+1)*(-1)^i, m = 1 for n >= 0. See A130534, A370518 for m=0 and m=2. - Igor Victorovich Statsenko, Feb 27 2024

A137650 Triangle read by rows, A008277 * A000012.

Original entry on oeis.org

1, 2, 1, 5, 4, 1, 15, 14, 7, 1, 52, 51, 36, 11, 1, 203, 202, 171, 81, 16, 1, 877, 876, 813, 512, 162, 22, 1, 4140, 4139, 4012, 3046, 1345, 295, 29, 1, 21147, 21146, 20891, 17866, 10096, 3145, 499, 37, 1, 115975, 115974, 115463
Offset: 1

Views

Author

Gary W. Adamson, Feb 01 2008

Keywords

Comments

Left column = Bell numbers (A000110) starting (1, 2, 5, 15, 52, 203, ...). Row sums = A005493(n+1): (1, 3, 10, 37, 151, 674, ...).
Corresponding to the generalized Stirling number triangle of first kind A049444. - Peter Luschny, Sep 18 2011

Examples

			First few rows of the triangle are
    1;
    2,   1;
    5,   4,   1;
   15,  14,   7,   1;
   52,  51,  36,  11,   1;
  203, 202, 171,  81,  16,   1;
  877, 876, 813, 512, 162,  22,   1;
  ...
		

Crossrefs

A similar triangle is A133611.

Programs

  • Maple
    A137650_row := proc(n) local k,i;
    add(add(combinat[stirling2](n, n-i), i=0..k)*x^(n-k-1),k=0..n-1);
    seq(coeff(%,x,k),k=0..n-1) end:
    seq(print(A137650_row(n)),n=1..7); # Peter Luschny, Sep 18 2011
  • Mathematica
    row[n_] := Table[StirlingS2[n, k], {k, 0, n}] // Reverse // Accumulate // Reverse // Rest;
    Array[row, 10] // Flatten (* Jean-François Alcover, Dec 07 2019 *)

Formula

A008277 * A000012 as infinite lower triangular matrices. Partial sums of A008277 rows starting from the right.

A051523 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -10, 1, 110, -21, 1, -1320, 362, -33, 1, 17160, -6026, 791, -46, 1, -240240, 101524, -17100, 1435, -60, 1, 3603600, -1763100, 358024, -38625, 2335, -75, 1, -57657600, 31813200, -7491484, 976024, -75985, 3535, -91, 1, 980179200, -598482000, 159168428, -24083892, 2267769, -136080, 5082, -108, 1
Offset: 0

Views

Author

Keywords

Comments

a(n,m)= ^10P_n^m in the notation of the given reference with a(0,0) := 1. The monic row polynomials s(n,x) := sum(a(n,m)*x^m,m=0..n) which are s(n,x)= product(x-(10+k),k=0..n-1), n >= 1 and s(0,x)=1 satisfy s(n,x+y) = sum(binomial(n,k)*s(k,x)*S1(n-k,y),k=0..n), with the Stirling1 polynomials S1(n,x)=sum(A008275(n,m)*x^m, m=1..n) and S1(0,x)=1. In the umbral calculus (see the S. Roman reference given in A048854) the s(n,x) polynomials are called Sheffer for (exp(10*t),exp(t)-1).

Examples

			{1}; {-10,1}; {110,-21,1}; {-1320,362,-331}; ... s(2,x)= 110-21*x+x^2; S1(2,x)= -x+x^2 (Stirling1).
		

Crossrefs

The first (m=0) unsigned column sequence is A049398. Row sums (signed triangle): A049389(n)*(-1)^n. Row sums (unsigned triangle): A051431(n).

Programs

  • Haskell
    a051523 n k = a051523_tabl !! n !! k
    a051523_row n = a051523_tabl !! n
    a051523_tabl = map fst $ iterate (\(row, i) ->
       (zipWith (-) ([0] ++ row) $ map (* i) (row ++ [0]), i + 1)) ([1], 10)
    -- Reinhard Zumkeller, Mar 12 2014
  • Mathematica
    a[n_, m_] := Pochhammer[m + 1, n - m] SeriesCoefficient[Log[1 + x]^m/(1 + x)^10, {x, 0, n}];
    Table[a[n, m], {n, 0, 8}, {m, 0, n}] // Flatten (* Jean-François Alcover, Oct 29 2019 *)

Formula

a(n, m)= a(n-1, m-1) - (n+9)*a(n-1, m), n >= m >= 0; a(n, m) := 0, n
E.g.f. for m-th column of signed triangle: ((log(1+x))^m)/(m!*(1+x)^10).
Triangle (signed) = [ -10, -1, -11, -2, -12, -3, -13, -14, -4, ...] DELTA A000035; triangle (unsigned) = [10, 1, 11, 2, 12, 3, 13, 4, 14, 5, 15, ...] DELTA A000035; where DELTA is Deléham's operator defined in A084938.
If we define f(n,i,a)=sum(binomial(n,k)*stirling1(n-k,i)*product(-a-j,j=0..k-1),k=0..n-i), then T(n,i) = f(n,i,10), for n=1,2,...;i=0...n. - Milan Janjic, Dec 21 2008

A074060 Graded dimension of the cohomology ring of the moduli space of n-pointed stable curves of genus 0 satisfying the associativity equations of physics (also known as the WDVV equations).

Original entry on oeis.org

1, 1, 1, 1, 5, 1, 1, 16, 16, 1, 1, 42, 127, 42, 1, 1, 99, 715, 715, 99, 1, 1, 219, 3292, 7723, 3292, 219, 1, 1, 466, 13333, 63173, 63173, 13333, 466, 1, 1, 968, 49556, 429594, 861235, 429594, 49556, 968, 1, 1, 1981, 173570, 2567940, 9300303, 9300303, 2567940, 173570, 1981, 1
Offset: 3

Author

Margaret A. Readdy, Aug 16 2002

Keywords

Comments

Combinatorial interpretations of Lagrange inversion (A134685) and the 2-Stirling numbers of the first kind (A049444 and A143491) provide a combinatorial construction for A074060 (see first Copeland link). For relations of A074060 to other arrays see second Copeland link page 19. - Tom Copeland, Sep 28 2008
These Poincare polynomials for the compactified moduli space of rational curves are presented on p. 5 of Lando and Zvonkin as well as those for the non-compactified Poincare polynomials of A049444 in factorial form. - Tom Copeland, Jun 13 2021

Examples

			Viewed as a triangular array, the values are
  1;
  1,   1;
  1,   5,   1;
  1,  16,  16,   1;
  1,  42, 127,  42,   1; ...
		

Crossrefs

Cf. A074059. 2nd diagonal is A002662.

Programs

  • Maple
    DA:=((1+t)*A(u,t)+u)/(1-t*A(u,t)): F:=0: for k from 1 to 10 do F:=map(simplify,int(series(subs(A(u,t)=F,DA),u,k),u)); od: # Eric Rains, Apr 02 2005
  • Mathematica
    DA = ((1+t) A[u, t] + u)/(1 - t A[u, t]); F = 0;
    Do[F = Integrate[Series[DA /. A[u, t] -> F, {u, 0, k}], u], {k, 1, 10}];
    (cc = CoefficientList[#, t]; cc Denominator[cc[[1]]])& /@ Drop[ CoefficientList[F, u], 2] // Flatten (* Jean-François Alcover, Oct 15 2019, after Eric Rains *)

Formula

Define offset to be 0 and P(n,t) = (-1)^n Sum_{j=0..n-2} a(n-2,j)*t^j with P(1,t) = -1 and P(0,t) = 1, then H(x,t) = -1 + exp(P(.,t)*x) is the compositional inverse in x about 0 of G(x,t) in A049444. H(x,0) = exp(-x) - 1, H(x,1) = -1 + exp( 2 + W( -exp(-2) * (2-x) ) ) and H(x,2) = 1 - (1+2*x)^(1/2), where W is a branch of the Lambert function such that W(-2*exp(-2)) = -2. - Tom Copeland, Feb 17 2008
Let offset=0 and g(x,t) = (1-t)/((1+x)^(t-1)-t), then the n-th row polynomial of the table is given by [(g(x,t)*D_x)^(n+1)]x with the derivative evaluated at x=0. - Tom Copeland, Jun 01 2008
With the notation in Copeland's comments, dH(x,t)/dx = -g(H(x,t),t). - Tom Copeland, Sep 01 2011
The term linear in x of [x*g(d/dx,t)]^n 1 gives the n-th row polynomial with offset 1. (See A134685.) - Tom Copeland, Oct 21 2011

Extensions

More terms from Eric Rains, Apr 02 2005

A087748 Triangle formed by reading triangle of Stirling numbers of the first kind (A048994) mod 2.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1
Offset: 0

Author

Philippe Deléham, Oct 02 2003

Keywords

Examples

			Triangle begins:
1,
0, 1,
0, 1, 1,
0, 0, 1, 1,
0, 0, 1, 0, 1,
0, 0, 0, 1, 0, 1,
0, 0, 0, 1, 1, 1, 1,
0, 0, 0, 0, 1, 1, 1, 1,
0, 0, 0, 0, 1, 0, 0, 0, 1,
0, 0, 0, 0, 0, 1, 0, 0, 0, 1,
0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1,
0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1,
0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1,
...
		

References

  • Brand, Neal; Das, Sajal; Jacob, Tom. The number of nonzero entries in recursively defined tables modulo primes. Proceedings of the Twenty-first Southeastern Conference on Combinatorics, Graph Theory, and Computing (Boca Raton, FL, 1990). Congr. Numer. 78 (1990), 47--59. MR1140469 (92h:05004). - From N. J. A. Sloane, Jun 03 2012

Crossrefs

Also parity of triangles A049444, A049459, A051338, A051379, A051523.

Formula

T(n, k) = A087755(n, k) = A048994(n, k) mod 2 = A047999([n/2], k-[(n+1)/2]) = T(n-2, k-2) XOR T(n-2, k-1) with T(0, 0) = T(1, 1) = 1 and T(1, 0) = 0; T(2n, k) = T(2n-1, k-1) XOR T(2n-1, k); T(2n+1, k) = T(2n, k-1). - Henry Bottomley, Dec 01 2003

Extensions

Edited and extended by Henry Bottomley, Dec 01 2003
Showing 1-10 of 16 results. Next