A007559 -id:A007559 - cp's OEIS frontend

A048994 Triangle of Stirling numbers of first kind, s(n,k), n >= 0, 0 <= k <= n.

1, 0, 1, 0, -1, 1, 0, 2, -3, 1, 0, -6, 11, -6, 1, 0, 24, -50, 35, -10, 1, 0, -120, 274, -225, 85, -15, 1, 0, 720, -1764, 1624, -735, 175, -21, 1, 0, -5040, 13068, -13132, 6769, -1960, 322, -28, 1, 0, 40320, -109584, 118124, -67284, 22449, -4536, 546, -36, 1, 0, -362880, 1026576, -1172700, 723680, -269325, 63273, -9450, 870, -45, 1

Offset: 0

N. J. A. Sloane

The unsigned numbers are also called Stirling cycle numbers: |s(n,k)| = number of permutations of n objects with exactly k cycles.
Mirror image of the triangle A054654. - Philippe Deléham, Dec 30 2006
Also the triangle gives coefficients T(n,k) of x^k in the expansion of C(x,n) = (a(k)*x^k)/n!. - Mokhtar Mohamed, Dec 04 2012
From Wolfdieter Lang, Nov 14 2018: (Start)
This is the Sheffer triangle of Jabotinsky type (1, log(1 + x)). See the e.g.f. of the triangle below.
This is the inverse Sheffer triangle of the Stirling2 Sheffer triangle A008275.
The a-sequence of this Sheffer triangle (see a W. Lang link in A006232)
is from the e.g.f. A(x) = x/(exp(x) -1) a(n) = Bernoulli(n) = A027641(n)/A027642(n), for n >= 0. The z-sequence vanishes.
The Boas-Buck sequence for the recurrences of columns has o.g.f. B(x) = Sum_{n>=0} b(n)*x^n = 1/((1 + x)*log(1 + x)) - 1/x. b(n) = (-1)^(n+1)*A002208(n+1)/A002209(n+1), b = {-1/2, 5/12, -3/8, 251/720, -95/288, 19087/60480,...}. For the Boas-Buck recurrence of Riordan and Sheffer triangles see the Aug 10 2017 remark in A046521, adapted to the Sheffer case, also for two references. See the recurrence and example below. (End)
Let G(n,m,k) be the number of simple labeled graphs on [n] with m edges and k components. Then T(n,k) = Sum (-1)^m*G(n,m,k). See the Read link below. Equivalently, T(n,k) = Sum mu(0,p) where the sum is over all set partitions p of [n] containing k blocks and mu is the Moebius function in the incidence algebra associated to the set partition lattice on [n]. - Geoffrey Critzer, May 11 2024

			Triangle begins:
  n\k 0     1       2       3      4      5      6    7    8   9 ...
  0   1
  1   0     1
  2   0    -1       1
  3   0     2      -3       1
  4   0    -6      11      -6      1
  5   0    24     -50      35    -10      1
  6   0  -120     274    -225     85    -15      1
  7   0   720   -1764    1624   -735    175    -21    1
  8   0 -5040   13068  -13132   6769  -1960    322  -28    1
  9   0 40320 -109584  118124 -67284  22449  -4536  546  -36   1
  ... - _Wolfdieter Lang_, Aug 22 2012
------------------------------------------------------------------
From _Wolfdieter Lang_, Nov 14 2018: (Start)
Recurrence: s(5,2)= s(4, 1) - 4*s(4, 2) = -6 - 4*11 = -50.
Recurrence from the a- and z-sequences: s(6, 3) = 2*(1*1*(-50) + 3*(-1/2)*35 + 6*(1/6)*(-10) + 10*0*1) = -225.
Boas-Buck recurrence for column k = 3, with b = {-1/2, 5/12, -3/8, ...}:
s(6, 3) = 6!*((-3/8)*1/3! + (5/12)*(-6)/4! + (-1/2)*35/5!) = -225. (End)

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 833.
L. Comtet, Advanced Combinatorics, Reidel, 1974; Chapter V, also p. 310.
J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 93.
F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 226.
R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 245.
J. Riordan, An Introduction to Combinatorial Analysis, p. 48.

T. D. Noe, Rows n=0..100 of triangle, flattened
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Paul Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Toda Chain Equations, Journal of Integer Sequences, 17 (2014), #14.2.3.
Fatima Zohra Bensaci, Rachid Boumahdi, and Laala Khaldi, Finite Sums Involving Fibonacci and Lucas Numbers, J. Int. Seq. (2024). See p. 9.
R. M. Dickau, Stirling numbers of the first kind. [Illustrates the unsigned Stirling cycle numbers A132393.]
Askar Dzhumadil'daev and Damir Yeliussizov, Walks, partitions, and normal ordering, Electronic Journal of Combinatorics, 22(4) (2015), #P4.10.
Bill Gosper, Colored illustrations of triangle of Stirling numbers of first kind read mod 2, 3, 4, 5, 6, 7.
Gergő Nemes, An Asymptotic Expansion for the Bernoulli Numbers of the Second Kind, J. Int. Seq. 14 (2011), #11.4.8.
A. Hennessy and Paul Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Orthogonal Polynomials, J. Int. Seq. 14 (2011), #11.8.2 (A-number typo A048894).
NIST Digital Library of Mathematical Functions, Stirling Numbers
Ken Ono, Larry Rolen, and Florian Sprung, Zeta-Polynomials for modular form periods, p. 6, arXiv:1602.00752 [math.NT], 2016.
Ricardo A. Podestá, New identities for binary Krawtchouk polynomials, binomial coefficients and Catalan numbers, arXiv preprint arXiv:1603.09156 [math.CO], 2016.
Ronald Read, An Introduction to Chromatic Polynomials, Journal of Combinatorial Theory, 4(1968)52-71.
Wikipedia, Stirling numbers and exponential generating functions.

See especially A008275 which is the main entry for this triangle. A132393 is an unsigned version, and A008276 is another version.
Cf. A008277, A039814-A039817, A048993, A084938.
A000142(n) = Sum_{k=0..n} |s(n, k)| for n >= 0.
Row sums give A019590(n+1).
Cf. A002209, A027641/A027642, A216294, A271705.

a048994 n k = a048994_tabl !! n !! k
a048994_row n = a048994_tabl !! n
a048994_tabl = map fst $ iterate (\(row, i) ->
(zipWith (-) ([0] ++ row) $ map (* i) (row ++ [0]), i + 1)) ([1], 0)
-- Reinhard Zumkeller, Mar 18 2013

A048994:= proc(n,k) combinat[stirling1](n,k) end: # R. J. Mathar, Feb 23 2009
seq(print(seq(coeff(expand(k!*binomial(x,k)),x,i),i=0..k)),k=0..9); # Peter Luschny, Jul 13 2009
A048994_row := proc(n) local k; seq(coeff(expand(pochhammer(x-n+1,n)), x,k), k=0..n) end: # Peter Luschny, Dec 30 2010

Table[StirlingS1[n, m], {n, 0, 9}, {m, 0, n}] (* Peter Luschny, Dec 30 2010 *)

create_list(stirling1(n,k),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */

a(n,k) = if(k<0 || k>n,0, if(n==0,1,(n-1)*a(n-1,k)+a(n-1,k-1)))

trg(nn)=for (n=0, nn-1, for (k=0, n, print1(stirling(n,k,1), ", ");); print();); \\ Michel Marcus, Jan 19 2015

s(n, k) = A008275(n,k) for n >= 1, k = 1..n; column k = 0 is {1, repeat(0)}.
s(n, k) = s(n-1, k-1) - (n-1)*s(n-1, k), n, k >= 1; s(n, 0) = s(0, k) = 0; s(0, 0) = 1.
The unsigned numbers a(n, k)=|s(n, k)| satisfy a(n, k)=a(n-1, k-1)+(n-1)*a(n-1, k), n, k >= 1; a(n, 0) = a(0, k) = 0; a(0, 0) = 1.
Triangle (signed) = [0, -1, -1, -2, -2, -3, -3, -4, -4, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...]; Triangle(unsigned) = [0, 1, 1, 2, 2, 3, 3, 4, 4, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...]; where DELTA is Deléham's operator defined in A084938.
Sum_{k=0..n} (-m)^(n-k)*s(n, k) = A000142(n), A001147(n), A007559(n), A007696(n), ... for m = 1, 2, 3, 4, ... .- Philippe Deléham, Oct 29 2005
A008275*A007318 as infinite lower triangular matrices. - Gerald McGarvey, Aug 20 2009
T(n,k) = n!*[x^k]([t^n]exp(x*log(1+t))). - Peter Luschny, Dec 30 2010, updated Jun 07 2020
From Wolfdieter Lang, Nov 14 2018: (Start)
Recurrence from the Sheffer a-sequence (see a comment above): s(n, k) = (n/k)*Sum_{j=0..n-k} binomial(k-1+j, j)*Bernoulli(j)*s(n-1, k-1+j), for n >= 1 and k >= 1, with s(n, 0) = 0 if n >= 1, and s(0,0) = 1.
Boas-Buck type recurrence for column k: s(n, k) = (n!*k/(n - k))*Sum_{j=k..n-1} b(n-1-j)*s(j, k)/j!, for n >= 1 and k = 0..n-1, with input s(n, n) = 1. For sequence b see the Boas-Buck comment above. (End)
T(n,k) = Sum_{j=k..n} (-1)^(n-j)*A271705(n,j)*A216294(j,k). - Mélika Tebni, Feb 23 2023

Offset corrected by R. J. Mathar, Feb 23 2009

Formula corrected by Philippe Deléham, Sep 10 2009

A132393 Triangle of unsigned Stirling numbers of the first kind (see A048994), read by rows, T(n,k) for 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 3, 1, 0, 6, 11, 6, 1, 0, 24, 50, 35, 10, 1, 0, 120, 274, 225, 85, 15, 1, 0, 720, 1764, 1624, 735, 175, 21, 1, 0, 5040, 13068, 13132, 6769, 1960, 322, 28, 1, 0, 40320, 109584, 118124, 67284, 22449, 4536, 546, 36, 1, 0, 362880, 1026576, 1172700, 723680, 269325, 63273, 9450, 870, 45, 1

Offset: 0

Philippe Deléham, Nov 10 2007, Oct 15 2008, Oct 17 2008

Another name: Triangle of signless Stirling numbers of the first kind.
Triangle T(n,k), 0<=k<=n, read by rows given by [0,1,1,2,2,3,3,4,4,5,5,...] DELTA [1,0,1,0,1,0,1,0,1,...] where DELTA is the operator defined in A084938.
A094645*A007318 as infinite lower triangular matrices.
Row sums are the factorial numbers. - Roger L. Bagula, Apr 18 2008
Exponential Riordan array [1/(1-x), log(1/(1-x))]. - Ralf Stephan, Feb 07 2014
Also the Bell transform of the factorial numbers (A000142). For the definition of the Bell transform see A264428 and for cross-references A265606. - Peter Luschny, Dec 31 2015
This is the lower triagonal Sheffer matrix of the associated or Jabotinsky type |S1| = (1, -log(1-x)) (see the W. Lang link under A006232 for the notation and references). This implies the e.g.f.s given below. |S1| is the transition matrix from the monomial basis {x^n} to the rising factorial basis {risefac(x,n)}, n >= 0. - Wolfdieter Lang, Feb 21 2017
T(n, k), for n >= k >= 1, is also the total volume of the n-k dimensional cell (polytope) built from the n-k orthogonal vectors of pairwise different lengths chosen from the set {1, 2, ..., n-1}. See the elementary symmetric function formula for T(n, k) and an example below. - Wolfdieter Lang, May 28 2017
From Wolfdieter Lang, Jul 20 2017: (Start)
The compositional inverse w.r.t. x of y = y(t;x) = x*(1 - t(-log(1-x)/x)) = x + t*log(1-x) is x = x(t;y) = ED(y,t) := Sum_{d>=0} D(d,t)*y^(d+1)/(d+1)!, the e.g.f. of the o.g.f.s D(d,t) = Sum_{m>=0} T(d+m, m)*t^m of the diagonal sequences of the present triangle. See the P. Bala link for a proof (there d = n-1, n >= 1, is the label for the diagonals).
This inversion gives D(d,t) =  P(d, t)/(1-t)^(2*d+1), with the numerator polynomials P(d, t) =  Sum_{m=0..d} A288874(d, m)*t^m. See an example below. See also the P. Bala formula in A112007. (End)
For n > 0, T(n,k) is the number of permutations of the integers from 1 to n which have k visible digits when viewed from a specific end, in the sense that a higher value hides a lower one in a subsequent position. - Ian Duff, Jul 12 2019

			Triangle T(n,k) begins:
  1;
  0,    1;
  0,    1,     1;
  0,    2,     3,     1;
  0,    6,    11,     6,    1;
  0,   24,    50,    35,   10,    1;
  0,  120,   274,   225,   85,   15,   1;
  0,  720,  1764,  1624,  735,  175,  21,  1;
  0, 5040, 13068, 13132, 6769, 1960, 322, 28, 1;
  ...
---------------------------------------------------
Production matrix is
  0, 1
  0, 1, 1
  0, 1, 2,  1
  0, 1, 3,  3,  1
  0, 1, 4,  6,  4,  1
  0, 1, 5, 10, 10,  5,  1
  0, 1, 6, 15, 20, 15,  6, 1
  0, 1, 7, 21, 35, 35, 21, 7, 1
  ...
From _Wolfdieter Lang_, May 09 2017: (Start)
Three term recurrence: 50 = T(5, 2) = 1*6 + (5-1)*11 = 50.
Recurrence from the Sheffer a-sequence [1, 1/2, 1/6, 0, ...]: 50 = T(5, 2) = (5/2)*(binomial(1, 1)*1*6 + binomial(2, 1)*(1/2)*11 + binomial(3, 1)*(1/6)*6 + 0) = 50. The vanishing z-sequence produces the k=0 column from T(0, 0) = 1. (End)
Elementary symmetric function T(4, 2) = sigma^{(3)}_2 = 1*2 + 1*3 + 2*3 = 11. Here the cells (polytopes) are 3 rectangles with total area 11. - _Wolfdieter Lang_, May 28 2017
O.g.f.s of diagonals: d=2 (third diagonal) [0, 6, 50, ...] has D(2,t) = P(2, t)/(1-t)^5, with P(2, t) = 2 + t, the n = 2 row of A288874. - _Wolfdieter Lang_, Jul 20 2017
Boas-Buck recurrence for column k = 2 and n = 5: T(5, 2) = (5!*2/3)*((3/8)*T(2,2)/2! + (5/12)*T(3,2)/3! + (1/2)*T(4,2)/4!) = (5!*2/3)*(3/16 + (5/12)*3/3! + (1/2)*11/4!) = 50. The beta sequence begins: {1/2, 5/12, 3/8, ...}. - _Wolfdieter Lang_, Aug 11 2017

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, pages 31, 187, 441, 996.
R. L. Graham, D. E. Knuth, and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd. ed., Table 259, p. 259.
Steve Roman, The Umbral Calculus, Dover Publications, New York (1984), pp. 149-150

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
Roland Bacher and P. De La Harpe, Conjugacy growth series of some infinitely generated groups, hal-01285685v2, 2016.
Eli Bagno and David Garber, Combinatorics of q,r-analogues of Stirling numbers of type B, arXiv:2401.08365 [math.CO], 2024. See page 5.
Peter Bala, Diagonals of triangles with generating function exp(t*F(x)).
J. Fernando Barbero G., Jesús Salas, and Eduardo J. S. Villaseñor, Bivariate Generating Functions for a Class of Linear Recurrences. I. General Structure, arXiv:1307.2010 [math.CO], 2013.
Jean-Luc Baril and Sergey Kirgizov, The pure descent statistic on permutations, Preprint, 2016.
Jean-Luc Baril and Sergey Kirgizov, Transformation à la Foata for special kinds of descents and excedances, arXiv:2101.01928 [math.CO], 2021.
Jean-Luc Baril and José L. Ramírez, Some distributions in increasing and flattened permutations, arXiv:2410.15434 [math.CO], 2024. See p. 9.
Ricky X. F. Chen, A Note on the Generating Function for the Stirling Numbers of the First Kind, Journal of Integer Sequences, 18 (2015), #15.3.8.
Bishal Deb and Alan D. Sokal, Higher-order Stirling cycle and subset triangles: Total positivity, continued fractions and real-rootedness, arXiv:2507.18959 [math.CO], 2025. See p. 5.
FindStat - Combinatorial Statistic Finder, The number of saliances of a permutation, The number of cycles in the cycle decomposition of a permutation.
Bill Gosper, Colored illustrations of triangle of Stirling numbers of first kind read mod 2, 3, 4, 5, 6, 7
W. Steven Gray and Makhin Thitsa, System Interconnections and Combinatorial Integer Sequences, in: System Theory (SSST), 2013 45th Southeastern Symposium on, Date of Conference: 11-11 March 2013, Digital Object Identifier: 10.1109/SSST.2013.6524939.
John M. Holte, Carries, Combinatorics and an Amazing Matrix, The American Mathematical Monthly, Vol. 104, No. 2 (Feb., 1997), pp. 138-149.
Tanya Khovanova and J. B. Lewis, Skyscraper Numbers, J. Int. Seq. 16 (2013) #13.7.2.
Sergey Kitaev and Philip B. Zhang, Distributions of mesh patterns of short lengths, arXiv:1811.07679 [math.CO], 2018.
Wolfdieter Lang, On Sums of Powers of Arithmetic Progressions, and Generalized Stirling, Eulerian and Bernoulli numbers, arXiv:1707.04451 [math.NT], 2017.
Shuzhen Lv and Philip B. Zhang, Joint equidistributions of mesh patterns 123 and 321 with symmetric and antipodal shadings, arXiv:2501.00357 [math.CO], 2024. See p. 11.
Shuzhen Lv and Philip B. Zhang, Joint equidistributions of mesh patterns 123 and 132 with antipodal shadings, arXiv:2506.23148 [math.CO], 2025. See p. 6.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO], 2012. - From _N. J. A. Sloane_, Aug 21 2012
Mathematics Stack Exchange, Symmetric (under the swapping) recursions for Stirling numbers of both kinds, Aug 10 2025.
Emanuele Munarini, Shifting Property for Riordan, Sheffer and Connection Constants Matrices, Journal of Integer Sequences, Vol. 20 (2017), Article 17.8.2.
Emanuele Munarini, Combinatorial identities involving the central coefficients of a Sheffer matrix, Applicable Analysis and Discrete Mathematics (2019) Vol. 13, 495-517.
Fedor Petrov, Recursive algorithm for columns of Stirling numbers of the first kind, answer to question on MathOverflow (2025).
E. G. Santos, Counting non-attacking chess pieces placements: Bishops and Anassas, arXiv:2411.16492 [math.CO], 2024. See p. 2.
Umesh Shankar, Log-concavity of rows of triangular arrays satisfying a certain super-recurrence, arXiv:2508.12467 [math.CO], 2025. See p. 4.
Xun-Tuan Su, Deng-Yun Yang, and Wei-Wei Zhang, A note on the generalized factorial, Australasian Journal of Combinatorics, Volume 56 (2013), Pages 133-137.
Benjamin Testart, Completing the enumeration of inversion sequences avoiding one or two patterns of length 3, arXiv:2407.07701 [math.CO], 2024. See p. 37.

Essentially a duplicate of A048994. Cf. A008275, A008277, A112007, A130534, A288874, A354795.

a132393 n k = a132393_tabl !! n !! k
a132393_row n = a132393_tabl !! n
a132393_tabl = map (map abs) a048994_tabl
-- Reinhard Zumkeller, Nov 06 2013

a132393_row := proc(n) local k; seq(coeff(expand(pochhammer (x,n)),x,k),k=0..n) end: # Peter Luschny, Nov 28 2010

p[t_] = 1/(1 - t)^x; Table[ ExpandAll[(n!)SeriesCoefficient[ Series[p[t], {t, 0, 30}], n]], {n, 0, 10}]; a = Table[(n!)* CoefficientList[SeriesCoefficient[ Series[p[t], {t, 0, 30}], n], x], {n, 0, 10}]; Flatten[a] (* Roger L. Bagula, Apr 18 2008 *)
Flatten[Table[Abs[StirlingS1[n,i]],{n,0,10},{i,0,n}]] (* Harvey P. Dale, Feb 04 2014 *)

create_list(abs(stirling1(n,k)),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */

column(n,k) = my(v1, v2); v1 = vector(n-1, i, 0); v2 = vector(n, i, 0); v2[1] = 1; for(i=1, n-1, v1[i] = (i+k)*(i+k-1)/2*v2[i]; for(j=1, i-1, v1[j] *= (i-j)*(i+k)/(i-j+2)); v2[i+1] = vecsum(v1)/i); v2 \\ generates n first elements of the k-th column starting from the first nonzero element. - Mikhail Kurkov, Mar 05 2025

T(n,k) = T(n-1,k-1)+(n-1)*T(n-1,k), n,k>=1; T(n,0)=T(0,k); T(0,0)=1.
Sum_{k=0..n} T(n,k)*x^(n-k) = A000012(n), A000142(n), A001147(n), A007559(n), A007696(n), A008548(n), A008542(n), A045754(n), A045755(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively. - Philippe Deléham, Nov 13 2007
Expand 1/(1-t)^x = Sum_{n>=0}p(x,n)*t^n/n!; then the coefficients of the p(x,n) produce the triangle. - Roger L. Bagula, Apr 18 2008
Sum_{k=0..n} T(n,k)*2^k*x^(n-k) = A000142(n+1), A000165(n), A008544(n), A001813(n), A047055(n), A047657(n), A084947(n), A084948(n), A084949(n) for x = 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Sep 18 2008
a(n) = Sum_{k=0..n} T(n,k)*3^k*x^(n-k) = A001710(n+2), A001147(n+1), A032031(n), A008545(n), A047056(n), A011781(n), A144739(n), A144756(n), A144758(n) for x=1,2,3,4,5,6,7,8,9,respectively. - Philippe Deléham, Sep 20 2008
Sum_{k=0..n} T(n,k)*4^k*x^(n-k) = A001715(n+3), A002866(n+1), A007559(n+1), A047053(n), A008546(n), A049308(n), A144827(n), A144828(n), A144829(n) for x=1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Sep 21 2008
Sum_{k=0..n} x^k*T(n,k) = x*(1+x)*(2+x)*...*(n-1+x), n>=1. - Philippe Deléham, Oct 17 2008
From Wolfdieter Lang, Feb 21 2017: (Start)
E.g.f. k-th column: (-log(1 - x))^k, k >= 0.
E.g.f. triangle (see the Apr 18 2008 Baluga comment): exp(-x*log(1-z)).
E.g.f. a-sequence: x/(1 - exp(-x)). See A164555/A027642. The e.g.f. for the z-sequence is 0. (End)
From Wolfdieter Lang, May 28 2017: (Start)
The row polynomials R(n, x) = Sum_{k=0..n} T(n, k)*x^k, for n >= 0, are R(n, x) = risefac(x,n-1) := Product_{j=0..n-1} x+j, with the empty product for n=0 put to 1. See the Feb 21 2017 comment above. This implies:
T(n, k) = sigma^{(n-1)}_(n-k), for n >= k >= 1, with the elementary symmetric functions sigma^{(n-1)}_m of degree m in the n-1 symbols 1, 2, ..., n-1, with binomial(n-1, m) terms. See an example below.(End)
Boas-Buck type recurrence for column sequence k: T(n, k) = (n!*k/(n - k)) * Sum_{p=k..n-1} beta(n-1-p)*T(p, k)/p!, for n > k >= 0, with input T(k, k) = 1, and beta(k) = A002208(k+1)/A002209(k+1). See a comment and references in A286718. - Wolfdieter Lang, Aug 11 2017
T(n,k) = Sum_{j=k..n} j^(j-k)*binomial(j-1, k-1)*A354795(n,j) for n > 0. - Mélika Tebni, Mar 02 2023
n-th row polynomial: n!*Sum_{k = 0..2*n} (-1)^k*binomial(-x, k)*binomial(-x, 2*n-k) = n!*Sum_{k = 0..2*n} (-1)^k*binomial(1-x, k)*binomial(-x, 2*n-k). - Peter Bala, Mar 31 2024
From Mikhail Kurkov, Mar 05 2025: (Start)
For a general proof of the formulas below via generating functions, see Mathematics Stack Exchange link.
Recursion for the n-th row (independently of other rows): T(n,k) = 1/(n-k)*Sum_{j=2..n-k+1} binomial(-k,j)*T(n,k+j-1)*(-1)^j for 1 <= k < n with T(n,n) = 1.
Recursion for the k-th column (independently of other columns): T(n,k) = 1/(n-k)*Sum_{j=2..n-k+1} (j-2)!*binomial(n,j)*T(n-j+1,k) for 1 <= k < n with T(n,n) = 1 (see Fedor Petrov link). (End)

A007661 Triple factorial numbers a(n) = n!!!, defined by a(n) = n*a(n-3), a(0) = a(1) = 1, a(2) = 2. Sometimes written n!3.

Original entry on oeis.org

1, 1, 2, 3, 4, 10, 18, 28, 80, 162, 280, 880, 1944, 3640, 12320, 29160, 58240, 209440, 524880, 1106560, 4188800, 11022480, 24344320, 96342400, 264539520, 608608000, 2504902400, 7142567040, 17041024000, 72642169600, 214277011200, 528271744000, 2324549427200

Offset: 0

N. J. A. Sloane, Mira Bernstein, Robert G. Wilson v

The triple factorial of a positive integer n is the product of the positive integers <= n that have the same residue modulo 3 as n. - Peter Luschny, Jun 23 2011

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
J. Spanier and K. B. Oldham, An Atlas of Functions, Hemisphere, NY, 1987, p. 23.

T. D. Noe, Table of n, a(n) for n = 0..200
Shyam Sunder Gupta, Fascinating Factorials, Exploring the Beauty of Fascinating Numbers, Springer (2025) Ch. 16, 411-442.
Eric Weisstein's World of Mathematics, Multifactorial.

Union of A007559, A008544 and A032031.
Cf. A000142, A006882 (= A001147 union A000165), A007662 (= union of A007696, A001813, A008545 and A047053), A085157, A085158.
Cf. A008585, A016777, A016789, A161474, A288055.

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

a007661 n k = a007661_list !! n
a007661_list = 1 : 1 : 2 : zipWith (*) a007661_list [3..]
-- Reinhard Zumkeller, Sep 20 2013

I:=[1,1,2];[n le 3 select I[n] else (n-1)*Self(n-3): n in [1..30]]; // Vincenzo Librandi, Nov 27 2015

A007661 := n -> mul(k, k = select(k -> k mod 3 = n mod 3, [$1 .. n])): seq(A007661(n), n = 0 .. 29);  # Peter Luschny, Jun 23 2011

multiFactorial[n_, k_] := If[n < 1, 1, If[n < k + 1, n, n*multiFactorial[n - k, k]]]; Array[ multiFactorial[#, 3] &, 30, 0] (* Robert G. Wilson v, Apr 23 2011 *)
RecurrenceTable[{a[0]==a[1]==1,a[2]==2,a[n]==n*a[n-3]},a,{n,30}] (* Harvey P. Dale, May 17 2012 *)
Table[With[{q = Quotient[n + 2, 3]}, 3^q q! Binomial[n/3, q]], {n, 0, 30}] (* Jan Mangaldan, Mar 21 2013 *)
a[ n_] := With[{m = Mod[n, 3, 1], q = 1 + Quotient[n, 3, 1]}, If[n < 0, 0, 3^q Pochhammer[m/3, q]]]; (* Michael Somos, Feb 24 2019 *)
Table[Times@@Range[n,1,-3],{n,0,30}] (* Harvey P. Dale, Sep 12 2020 *)

A007661(n,d=3)=prod(i=0,(n-1)\d,n-d*i) \\ M. F. Hasler, Feb 16 2008

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

a(n) = Product_{i=0..floor((n-1)/3)} (n-3*i). - M. F. Hasler, Feb 16 2008
a(n) ~ c * n^(n/3+1/2)/exp(n/3), where c = sqrt(2*Pi/3) if n=3*k, c = sqrt(2*Pi)*3^(1/6) / Gamma(1/3) if n=3*k+1, c = sqrt(2*Pi)*3^(-1/6) / Gamma(2/3) if n=3*k+2. - Vaclav Kotesovec, Jul 29 2013
a(3*n) = A032031(n); a(3*n+1) = A007559(n+1); a(3*n+2) = A008544(n+1). - Reinhard Zumkeller, Sep 20 2013
0 = a(n)*(a(n+1) -a(n+4)) +a(n+1)*a(n+3) for all n>=0. - Michael Somos, Feb 24 2019
Sum_{n>=0} 1/a(n) = A288055. - Amiram Eldar, Nov 10 2020

A007696 Quartic (or 4-fold) factorial numbers: a(n) = Product_{k = 0..n-1} (4*k + 1).

Original entry on oeis.org

1, 1, 5, 45, 585, 9945, 208845, 5221125, 151412625, 4996616625, 184874815125, 7579867420125, 341094033905625, 16713607661375625, 885821206052908125, 50491808745015763125, 3080000333445961550625, 200200021673987500790625, 13813801495505137554553125

Offset: 0

N. J. A. Sloane

a(n), n >= 1, enumerates increasing quintic (5-ary) trees. See David Callan's comment on A007559 (number of increasing quarterny trees).

Hankel transform is A169619. - Paul Barry, Dec 03 2009

			G.f. = 1 + x + 5*x^2 + 45*x^3 + 585*x^4 + 9945*x^5 + 208845*x^6 + ...

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

T. D. Noe, Table of n, a(n) for n = 0..100
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seq. 3 (2000), Article 00.2.4.
J.-C. Novelli and J.-Y. Thibon, Hopf Algebras of m-permutations, (m+1)-ary trees, and m-parking functions, arXiv:1403.5962 [math.CO], 2014.
Maxie D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications , J. Integer Seq. 13 (2010), Article 10.6.7; see page 39.
Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, J. Integer Seq. 9 (2006), Article 06.1.1.

Cf. A001147, A001813, A004981, A007559, A008545, A034255, A047053, A051142, A264781.

a(n) = A049029(n, 1) for n >= 1 (first column of triangle).

a:=[1,1];; for n in [3..20] do a[n]:=(4*(n-1)-7)*(a[n-1]+4*a[n-2]); od; a; # G. C. Greubel, Aug 15 2019

[n le 2 select 1 else (4*(n-1)-7)*(Self(n-1) + 4*Self(n-2)): n in [1..20]]; // G. C. Greubel, Aug 15 2019

x:='x'; G(x):=(1-4*x)^(-1/4): f[0]:=G(x): for n from 1 to 29 do f[n]:=diff(f[n-1],x) od: seq(eval(f[n],x=0),n=0..17);# Zerinvary Lajos, Apr 03 2009
A007696 := n -> mul(k, k = select(k-> k mod 4 = 1, [$ 1 .. 4*n])): seq(A007696(n), n=0..17); # Peter Luschny, Jun 23 2011

a[ n_]:= Pochhammer[ 1/4, n] 4^n; (* Michael Somos, Jan 17 2014 *)
a[ n_]:= If[n < 0, 1 / Product[ -k, {k, 3, -4n-1, 4}], Product[ k, {k, 1, 4n-3, 4}]]; (* Michael Somos, Jan 17 2014 *)
Range[0, 19]! CoefficientList[Series[((1-4x)^(-1/4)), {x, 0, 19}], x] (* Vincenzo Librandi, Oct 08 2015 *)

A007696(n):=prod(4*k+1,k,0,n-1)$
makelist(A007696(n),n,0,30); /* Martin Ettl, Nov 05 2012 */

{a(n) = if( n<0, 1 / prod(k=1, -n, 1 - 4*k), prod(k=1, n, 4*k - 3))}; /* Michael Somos, Jan 17 2014 */

[4^n*rising_factorial(1/4, n) for n in (0..20)] # G. C. Greubel, Aug 15 2019

E.g.f.: (1 - 4*x)^(-1/4).
a(n) ~ 2^(1/2) * Pi^(1/2) * Gamma(1/4)^(-1) * n^(-1/4) * 2^(2*n) * e^(-n) * n^n * (1 - 1/96 * n^(-1) - ...). - Joe Keane (jgk(AT)jgk.org), Nov 23 2001 [corrected by Vaclav Kotesovec, Jul 19 2025]
a(n) = Sum_{k = 0..n} (-4)^(n-k) * A048994(n, k). - Philippe Deléham, Oct 29 2005
G.f.: 1/(1 - x/(1 - 4*x/(1 - 5*x/(1 - 8*x/(1 - 9*x/(1 - 12*x/(1 - 13*x/(1 - .../(1 - A042948(n+1)*x/(1 -... (continued fraction). - Paul Barry, Dec 03 2009
a(n) = (-3)^n * Sum_{k = 0..n} (4/3)^k * s(n+1, n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: 1/T(0), where T(k) =  1 - x * (4*k + 1)/(1 - x * (4*k + 4)/T(k+1)) (continued fraction). - Sergei N. Gladkovskii, Mar 19 2013
G.f.: 1 + x/Q(0), where Q(k) = 1 + x + 2*(2*k - 1)*x - 4*x*(k+1)/Q(k+1) (continued fraction). - Sergei N. Gladkovskii, May 03 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x * (4*k + 1)/(x * (4*k + 1) + 1/G(k+1))) (continued fraction). - Sergei N. Gladkovskii, Jun 04 2013
0 = a(n) * (4*a(n+1) - a(n+2)) + a(n+1) * a(n+1) for all n in Z. - Michael Somos, Jan 17 2014
a(-n) = (-1)^n / A008545(n). - Michael Somos, Jan 17 2014
Let T(x) = 1/(1 - 3*x)^(1/3) be the e.g.f. for the sequence of triple factorial numbers A007559. Then the e.g.f. A(x) for the quartic factorial numbers satisfies T(Integral_{t = 0..x} A(t) dt) = A(x). (Cf. A007559 and A008548.) - Peter Bala, Jan 02 2015
O.g.f.: hypergeom([1, 1/4], [], 4*x). - Peter Luschny, Oct 08 2015
a(n) = A264781(4*n+1, n). - Alois P. Heinz, Nov 24 2015
a(n) = 4^n * Gamma(n + 1/4)/Gamma(1/4). - Artur Jasinski, Aug 23 2016
D-finite with recurrence: a(n) +(-4*n+3)*a(n-1)=0, n>=1. - R. J. Mathar, Feb 14 2020
Sum_{n>=0} 1/a(n) = 1 + exp(1/4)*(Gamma(1/4) - Gamma(1/4, 1/4))/(2*sqrt(2)). - Amiram Eldar, Dec 18 2022

Better description from Wolfdieter Lang, Dec 11 1999

A008544 Triple factorial numbers: Product_{k=0..n-1} (3*k+2).

Original entry on oeis.org

1, 2, 10, 80, 880, 12320, 209440, 4188800, 96342400, 2504902400, 72642169600, 2324549427200, 81359229952000, 3091650738176000, 126757680265216000, 5577337931669504000, 262134882788466688000

Offset: 0

Joe Keane (jgk(AT)jgk.org)

a(n-1), n >= 1, enumerates increasing plane (aka ordered) trees with n vertices (one of them a root labeled 1) where each vertex with outdegree r >= 0 comes in r+1 types (like an (r+1)-ary vertex). See the increasing tree comments under A004747. - Wolfdieter Lang, Oct 12 2007
An example for the case of 3 vertices is shown below. For the enumeration of non-plane trees of this type see A029768. - Peter Bala, Aug 30 2011
a(n) is the product of the positive integers k <= 3*n that have k modulo 3 = 2. - Peter Luschny, Jun 23 2011
See A094638 for connections to differential operators. - Tom Copeland, Sep 20 2011
Partial products of A016789. - Reinhard Zumkeller, Sep 20 2013
The Mathar conjecture is true. Generally from the factorial form, the last term is the "extra" product beyond the prior term, from k=n-1 and 3k+2 evaluates to 3*(n-1)+2 = 3n-1, yielding a(n) = a(n-1)*(3n-1) (eqn1). Similarly, a(n) = a(n-2)*(3n-1)*(3(n-2)+2) = a(n-2)*(3n-1)*(3n-4) (eqn2) and a(n) = a(n-3)*(3n-1)*(3n-4)*(3*(n-2)+2) = a(n-3)*(3n-1)*(3n-4)*(3n-7) (eqn3). We equate (eqn2) and (eqn3) to get a(n-2)*(3n-1)*(3n-4) = a(n-3)*(3n-1)*(3n-4)*(3n-7) or a(n-2)+(7-3n)*a(n-3) = 0 (eqn4). From (eqn1) we have a(n)+(1-3n)*a(n-1) = 0 (eqn5). Combining (eqn4) and (eqn5) yields a(n)+(1-3n)*a(n-1)+a(n-2)+(7-3n)*a(n-3) = 0. - Bill McEachen, Jan 01 2016
a(n-1), n>=1, is the dimension of the n-th component of the operad encoding the multilinearization of the following identity in nonassociative algebras: s*(a,a,b)-(s+t)*(a,b,a)+t*(b,a,a)=0, for any given pair of scalars (s,t). Here (a,b,c) is the associator (ab)c-a(bc). This is proved in the referenced article on associator dependent algebras by Bremner and me. - Vladimir Dotsenko, Mar 22 2022

			a(2) = 10 from the described trees with 3 vertices: there are three trees with a root vertex (label 1) with outdegree r=2 (like the three 3-stars each with one different ray missing) and the four trees with a root (r=1 and label 1) a vertex with (r=1) and a leaf (r=0). Assigning labels 2 and 3 yields 2*3+4=10 such trees.
a(2) = 10. The 10 possible plane increasing trees on 3 vertices, where vertices of outdegree 1 come in 2 colors (denoted a or b) and vertices of outdegree 2 come in 3 colors (a, b or c), are:
.
   1a    1b    1a    1b        1a       1b       1c
   |     |     |     |        / \      / \      / \
   2a    2b    2b    2a      2   3    2   3    2   3
   |     |     |     |
   3     3     3     3         1a       1b       1c
                              / \      / \      / \
                             3   2    3   2    3   2

T. D. Noe, Table of n, a(n) for n = 0..100
Murray Bremner and Vladimir Dotsenko, Associator dependent algebras and Koszul duality, arXiv:2203.11142 [math.RA], 2022.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seq., Vol. 3 (2000), Article 00.2.4.
Keiichi Shigechi, On the lattice of weighted partitions, arXiv:2212.14666 [math.CO], 2022. See p. 27.
Garrett Southwood and Hua Wang, The Statistics and Combinatorics of Increasing Trees and Colored Increasing Trees, J. Combin. Math. Combin. Comput. (2024) Vol. 123, 475-487. See p. 485.

a(n) = A004747(n+1, 1) (first column of triangle). Cf. A051141.
Cf. A000165, A001813, A047055, A047657, A084947, A084948, A084949.
Cf. A029768, A034724, A049308.
Cf. A024396, A193534.
Cf. A225470, A290596 (first columns).
Subsequence of A007661.

a008544 n = a008544_list !! n
a008544_list = scanl (*) 1 a016789_list
-- Reinhard Zumkeller, Sep 20 2013

[Round((Gamma(2*n-5/3)/Gamma(n-5/6)*Gamma(2/3)/Gamma(5/6) )/ Sqrt(3)*3^n/4^(n-1)): n in [1..20]]; // Vincenzo Librandi, Feb 21 2015

[Round(3^n*Gamma(n+2/3)/Gamma(2/3)): n in [0..20]]; // G. C. Greubel, Mar 31 2019

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

k = 3; b[1]=2; b[n_]:= b[n] = b[n-1]+k; a[0]=1; a[1]=2; a[n_]:= a[n] = a[n-1]*b[n]; Table[a[n], {n,0,20}] (* Roger L. Bagula, Sep 17 2008 *)
Product[3 k + 2, {k, 0, # - 1}] & /@ Range[0, 16] (* Michael De Vlieger, Jan 02 2016 *)
Table[3^n*Pochhammer[2/3, n], {n,0,20}] (* G. C. Greubel, Mar 31 2019 *)

a(n):=((n)!*sum(binomial(k,n-k)*binomial(n+k,k)*3^(-n+k)*(-1)^(n-k),k,floor(n/2),n)); /* Vladimir Kruchinin, Sep 28 2013 */

PARI
```
a(n) = prod(k=0,n-1, 3*k+2 );
```

vector(20, n, n--; round(3^n*gamma(n+2/3)/gamma(2/3))) \\ G. C. Greubel, Mar 31 2019

@CachedFunction
def A008544(n): return 1 if n == 0 else (3*n-1)*A008544(n-1)
[A008544(n) for n in (0..16)]  # Peter Luschny, May 20 2013

[3^n*rising_factorial(2/3, n) for n in (0..20)] # G. C. Greubel, Mar 31 2019

a(n) = Product_{k=0..n-1} (3*k+2) = A007661(3*n-1) (with A007661(-1) = 1).
E.g.f.: (1-3*x)^(-2/3).
a(n) = 2*A034000(n), n >= 1, a(0) = 1.
a(n) ~ 2^(1/2)*Pi^(1/2)*Gamma(2/3)^-1*n^(1/6)*3^n*e^-n*n^n*{1 - 1/36*n^-1 + ...}. - Joe Keane (jgk(AT)jgk.org), Nov 22 2001
a(n) = (Gamma(2*n-5/3)/Gamma(n-5/6)*Gamma(2/3)/Gamma(5/6))/sqrt(3)*3^n/4^(n-1). - Jeremy L. Martin, Mar 31 2002 (typo fixed by Vincenzo Librandi, Feb 21 2015)
From Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003: (Start)
a(n) = A084939(n)/A000142(n)*A000079(n).
a(n) = 3^n*Pochhammer(2/3, n) = 3^n*Gamma(n+2/3)/Gamma(2/3). (End)
Let T = A094638 and c(t) = column vector(1, t, t^2, t^3, t^4, t^5,...), then A008544 = unsigned [ T * c(-3) ] and the list partition transform A133314 of [1,T * c(-3)] gives [1,T * c(3)] with all odd terms negated, which equals a signed version of A007559; i.e., LPT[(1,signed A008544)] = signed A007559. Also LPT[A007559] = (1,-A008544) and e.g.f. [1,T * c(t)] = (1-x*t)^(-1/t) for t = 3 or -3. Analogous results hold for the double factorial, quadruple factorial and so on. - Tom Copeland, Dec 22 2007
G.f.: 1/(1-2x/(1-3x/(1-5x/(1-6x/(1-8x/(1-9x/(1-11x/(1-12x/(1-...))))))))) (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-1)^n*Sum_{k=0..n} 3^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: 1/Q(0) where Q(k) = 1 - x*(3*k+2)/(1 - x*(3*k+3)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 20 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(3*k+2)/(x*(3*k+2) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 25 2013
D-finite with recurrence: a(n) = (9*(n-2)*(n-1)+2)*a(n-2) + 4*a(n-1), n>=2. - Ivan N. Ianakiev, Aug 09 2013
a(n) = n!*Sum_{k=floor(n/2)..n} binomial(k,n-k)*binomial(n+k,k)*3^(-n+k)*(-1)^(n-k). - Vladimir Kruchinin, Sep 28 2013
Recurrence equation: a(n) = 3*a(n-1) + (3*n - 4)^2*a(n-2) with a(0) = 1 and a(1) = 2. A024396 satisfies the same recurrence (but with different initial conditions). This observation leads to a continued fraction expansion for the constant A193534 due to Euler. - Peter Bala, Feb 20 2015
a(n) = A225470(n, 0), n >= 0. - Wolfdieter Lang, May 29 2017
G.f.: Hypergeometric2F0(1, 2/3; -; 3*x). - G. C. Greubel, Mar 31 2019
D-finite with recurrence: a(n) + (-3*n+1)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
G.f.: 1/(1-2*x-6*x^2/(1-8*x-30*x^2/(1-14*x-72*x^2/(1-20*x-132*x^2/(1-...))))) (Jacobi continued fraction). - Nikolaos Pantelidis, Feb 28 2020
G.f.: 1/G(0), where G(k) = 1 - (6*k+2)*x - 3*(k+1)*(3*k+2)*x^2/G(k+1). - Nikolaos Pantelidis, Feb 28 2020
Sum_{n>=0} 1/a(n) = 1 + (e/3)^(1/3) * (Gamma(2/3) - Gamma(2/3, 1/3)). - Amiram Eldar, Mar 01 2022

A094638 Triangle read by rows: T(n,k) = |s(n,n+1-k)|, where s(n,k) are the signed Stirling numbers of the first kind A008276 (1 <= k <= n; in other words, the unsigned Stirling numbers of the first kind in reverse order).

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 6, 11, 6, 1, 10, 35, 50, 24, 1, 15, 85, 225, 274, 120, 1, 21, 175, 735, 1624, 1764, 720, 1, 28, 322, 1960, 6769, 13132, 13068, 5040, 1, 36, 546, 4536, 22449, 67284, 118124, 109584, 40320, 1, 45, 870, 9450, 63273, 269325, 723680, 1172700, 1026576, 362880

Offset: 1

André F. Labossière, May 17 2004

Triangle of coefficients of the polynomial (x+1)(x+2)...(x+n), expanded in decreasing powers of x. - T. D. Noe, Feb 22 2008
Row n also gives the number of permutation of 1..n with complexity 0,1,...,n-1. See the comments in A008275. - N. J. A. Sloane, Feb 08 2019
T(n,k) is the number of deco polyominoes of height n and having k columns. A deco polyomino is a directed column-convex polyomino in which the height, measured along the diagonal, is attained only in the last column. Example: T(2,1)=1 and T(2,2)=1 because the deco polyominoes of height 2 are the vertical and horizontal dominoes, having, respectively, 1 and 2 columns. - Emeric Deutsch, Aug 14 2006
Sum_{k=1..n} k*T(n,k) = A121586. - Emeric Deutsch, Aug 14 2006
Let the triangle U(n,k), 0 <= k <= n, read by rows, be given by [1,0,1,0,1,0,1,0,1,0,1,...] DELTA [1,1,2,2,3,3,4,4,5,5,6,...] where DELTA is the operator defined in A084938; then T(n,k) = U(n-1,k-1). - Philippe Deléham, Jan 06 2007
From Tom Copeland, Dec 15 2007: (Start)
Consider c(t) = column vector(1, t, t^2, t^3, t^4, t^5, ...).
Starting at 1 and sampling every integer to the right, we obtain (1,2,3,4,5,...). And T * c(1) = (1, 1*2, 1*2*3, 1*2*3*4,...), giving n! for n > 0. Call this sequence the right factorial (n+)!.
Starting at 1 and sampling every integer to the left, we obtain (1,0,-1,-2,-3,-4,-5,...). And T * c(-1) = (1, 1*0, 1*0*-1, 1*0*-1*-2,...) = (1, 0, 0, 0, ...), the left factorial (n-)!.
Sampling every other integer to the right, we obtain (1,3,5,7,9,...). T * c(2) = (1, 1*3, 1*3*5, ...) = (1,3,15,105,945,...), giving A001147 for n > 0, the right double factorial, (n+)!!.
Sampling every other integer to the left, we obtain (1,-1,-3,-5,-7,...). T * c(-2) = (1, 1*-1, 1*-1*-3, 1*-1*-3*-5,...) = (1,-1,3,-15,105,-945,...) = signed A001147, the left double factorial, (n-)!!.
Sampling every 3 steps to the right, we obtain (1,4,7,10,...). T * c(3) = (1, 1*4, 1*4*7,...) = (1,4,28,280,...), giving A007559 for n > 0, the right triple factorial, (n+)!!!.
Sampling every 3 steps to the left, we obtain (1,-2,-5,-8,-11,...), giving T * c(-3) = (1, 1*-2, 1*-2*-5, 1*-2*-5*-8,...) = (1,-2,10,-80,880,...) = signed A008544, the left triple factorial, (n-)!!!.
The list partition transform A133314 of [1,T * c(t)] gives [1,T * c(-t)] with all odd terms negated; e.g., LPT[1,T*c(2)] = (1,-1,-1,-3,-15,-105,-945,...) = (1,-A001147). And e.g.f. for [1,T * c(t)] = (1-xt)^(-1/t).
The above results hold for t any real or complex number. (End)
Let R_n(x) be the real and I_n(x) the imaginary part of Product_{k=0..n} (x + I*k). Then, for n=1,2,..., we have R_n(x) = Sum_{k=0..floor((n+1)/2)}(-1)^k*Stirling1(n+1,n+1-2*k)*x^(n+1-2*k), I_n(x) = Sum_{k=0..floor(n/2)}(-1)^(k+1)*Stirling1(n+1,n-2*k)*x^(n-2*k). - Milan Janjic, May 11 2008
T(n,k) is also the number of permutations of n with "reflection length" k (i.e., obtained from 12..n by k not necessarily adjacent transpositions). For example, when n=3, 132, 213, 321 are obtained by one transposition, while 231 and 312 require two transpositions. - Kyle Petersen, Oct 15 2008
From Tom Copeland, Nov 02 2010: (Start)
[x^(y+1) D]^n = x^(n*y) [T(n,1)(xD)^n + T(n,2)y (xD)^(n-1) + ... + T(n,n)y^(n-1)(xD)], with D the derivative w.r.t. x.
E.g., [x^(y+1) D]^4 = x^(4*y) [(xD)^4 + 6 y(xD)^3 + 11 y^2(xD)^2 + 6 y^3(xD)].
(xD)^m can be further expanded in terms of the Stirling numbers of the second kind and operators of the form x^j D^j. (End)
With offset 0, 0 <= k <= n: T(n,k) is the sum of products of each size k subset of {1,2,...,n}. For example, T(3,2) = 11 because there are three subsets of size two: {1,2},{1,3},{2,3}. 1*2 + 1*3 + 2*3 = 11. - Geoffrey Critzer, Feb 04 2011
The Kn11, Fi1 and Fi2 triangle sums link this triangle with two sequences, see the crossrefs. For the definitions of these triangle sums see A180662. The mirror image of this triangle is A130534. - Johannes W. Meijer, Apr 20 2011
T(n+1,k+1) is the elementary symmetric function a_k(1,2,...,n), n >= 0, k >= 0, (a_0(0):=1). See the T. D. Noe and Geoffrey Critzer comments given above. For a proof see the Stanley reference, p. 19, Second Proof. - Wolfdieter Lang, Oct 24 2011
Let g(t) = 1/d(log(P(j+1,-t)))/dt (see Tom Copeland's 2007 formulas). The Mellin transform (t to s) of t*Dirac[g(t)] gives Sum_{n=1..j} n^(-s), which as j tends to infinity gives the Riemann zeta function for Re(s) > 1. Dirac(x) is the Dirac delta function. The complex contour integral along a circle of radius 1 centered at z=1 of z^s/g(z) gives the same result. - Tom Copeland, Dec 02 2011
Rows are coefficients of the polynomial expansions of the Pochhammer symbol, or rising factorial, Pch(n,x) = (x+n-1)!/(x-1)!. Expansion of Pch(n,xD) = Pch(n,Bell(.,:xD:)) in a polynomial with terms :xD:^k=x^k*D^k gives the Lah numbers A008297. Bell(n,x) are the unsigned Bell polynomials or Stirling polynomials of the second kind A008277. - Tom Copeland, Mar 01 2014
From Tom Copeland, Dec 09 2016: (Start)
The Betti numbers, or dimension, of the pure braid group cohomology. See pp. 12 and 13 of the Hyde and Lagarias link.
Row polynomials and their products appear in presentation of the Jack symmetric functions of R. Stanley. See Copeland link on the Witt differential generator.
(End)
From Tom Copeland, Dec 16 2019: (Start)
The e.g.f. given by Copeland in the formula section appears in a combinatorial Dyson-Schwinger equation of quantum field theory in Yeats in Thm. 2 on p. 62 related to a Hopf algebra of rooted trees. See also the Green function on p. 70.
Per comments above, this array contains the coefficients in the expansion in polynomials of the Euler, or state number, operator xD of the rising factorials Pch(n,xD) = (xD+n-1)!/(xD-1)! = x [:Dx:^n/n!]x^{-1} = L_n^{-1}(-:xD:), where :Dx:^n = D^n x^n and :xD:^n = x^n D^n. The polynomials L_n^{-1} are the Laguerre polynomials of order -1, i.e., normalized Lah polynomials.
The Witt differential operators L_n = x^(n+1) D and the row e.g.f.s appear in Hopf and dual Hopf algebra relations presented by Foissy. The Witt operators satisfy L_n L_k - L_k L_n = (k-n) L_(n+k), as for the dual Hopf algebra. (End)

			Triangle starts:
  1;
  1,  1;
  1,  3,  2;
  1,  6, 11,  6;
  1, 10, 35, 50, 24;
...

M. Miyata and J. W. Son, On the complexity of permutations and the metric space of bijections, Tensor, 60 (1998), No. 1, 109-116 (MR1768839).
Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 18, equations 18:4:2 - 18:4:8 at page 151.
R. P. Stanley, Enumerative Combinatorics, Vol. 1, Cambridge University Press, 1997.

T. D. Noe, Rows n = 1..51 of triangle, flattened
E. Barcucci, A. Del Lungo and R. Pinzani, "Deco" polyominoes, permutations and random generation, Theoretical Computer Science, 159, 1996, 29-42.
F. Bergeron, Ph. Flajolet and B. Salvy, Varieties of Increasing Trees, Lecture Notes in Computer Science vol. 581, ed. J.-C. Raoult, Springer 1992, pp. 24-48. Added Mar 01 2014
F. Bergeron, Philippe Flajolet and Bruno Salvy, Varieties of increasing trees, HAL, Rapport De Recherche Inria. Added Mar 01 2014
Tom Copeland, Addendum to Mathemagical Forests
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
Tom Copeland, Mellin Interpolation of Differential Ops and Associated Infinigens and Appell Polynomials: The Ordered, Laguerre, and Scherk-Witt-Lie Diff Ops
Tom Copeland, Witt differential generator for special Jack symmetric functions / polynomials
FindStat - Combinatorial Statistic Finder, The absolute length of a permutation
L. Foissy, Faa di Bruno subalgebras of the Hopf algebra of planar trees from combinatorial Dyson-Schwinger equations, arXiv:0707.1204 [math.RA], (2007).
O. Furdui, T. Trif, On the Summation of Certain Iterated Series, J. Int. Seq. 14 (2011) #11.6.1
F. Hivert, J.-C. Novelli and J.-Y. Thibon, The Algebra of Binary Search Trees, Theoretical Computer Science, 339 (2005), 129-165.
T. Hyde and J. Lagarias Polynomial splitting measures and cohomology of the pure braid group, arXiv preprint arXiv:1604.05359 [math.RT], 2016.
MathOverflow, Motivation of Virasoro algebra, an answer by Tom Copeland to an MO question posed in 2012.
Romeo Mestrovic, Lucas' theorem: its generalizations, extensions and applications (1878--2014), arXiv preprint arXiv:1409.3820 [math.NT], 2014.
Robert E. Moritz, On the sum of products of n consecutive integers, Univ. Washington Publications in Math., 1 (No. 3, 1926), 44-49 [Annotated scanned copy]
M. D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications , J. Int. Seq. 13 (2010), 10.6.7.
M. Z. Spivey, On Solutions to a General Combinatorial Recurrence, J. Int. Seq. 14 (2011) # 11.9.7.
K. Yeats, A Combinatorial Perspective on Quantum Field Theory, SpringerBriefs in Mathematical Physics, Vol. 15, 2017.

A008276 gives the (signed) Stirling numbers of the first kind.
Cf. A000108, A014137, A001246, A033536, A000984, A094639, A006134, A082894, A002897, A079727, A000217 (2nd column), A000914 (3rd column), A001303 (4th column), A000915 (5th column), A053567 (6th column), A000142 (row sums).
Triangle sums (see the comments): A124380 (Kn11), A001710 (Fi1, Fi2). - Johannes W. Meijer, Apr 20 2011
Cf. A121586, A130534, A143491.

Flat(List([1..10], n-> List([1..n], k-> Stirling1(n,n-k+1) ))); # G. C. Greubel, Dec 29 2019

a094638 n k = a094638_tabl !! (n-1) !! (k-1)
a094638_row n = a094638_tabl !! (n-1)
a094638_tabl = map reverse a130534_tabl
-- Reinhard Zumkeller, Aug 01 2014

[(-1)^(k+1)*StirlingFirst(n,n-k+1): k in [1..n], n in [1..10]]; // G. C. Greubel, Dec 29 2019

T:=(n,k)->abs(Stirling1(n,n+1-k)): for n from 1 to 10 do seq(T(n,k),k=1..n) od; # yields sequence in triangular form. # Emeric Deutsch, Aug 14 2006

Table[CoefficientList[Series[Product[1 + i x, {i,n}], {x,0,20}], x], {n,0,6}] (* Geoffrey Critzer, Feb 04 2011 *)
Table[Abs@StirlingS1[n, n-k+1], {n, 10}, {k, n}]//Flatten (* Michael De Vlieger, Aug 29 2015 *)

create_list(abs(stirling1(n+1,n-k+1)),n,0,10,k,0,n); /* Emanuele Munarini, Jun 01 2012 */

{T(n,k)=if(n<1 || k>n,0,(n-1)!*polcoeff(polcoeff(x*y/(1 - x*y+x*O(x^n))^(1 + 1/y),n,x),k,y))} /* Paul D. Hanna, Jul 21 2011 */

[[stirling_number1(n, n-k+1) for k in (1..n)] for n in (1..10)] # G. C. Greubel, Dec 29 2019

With P(n,t) = Sum_{k=0..n-1} T(n,k+1) * t^k = 1*(1+t)*(1+2t)...(1+(n-1)*t) and P(0,t)=1, exp[P(.,t)*x] = (1-tx)^(-1/t). T(n,k+1) = (1/k!) (D_t)^k (D_x)^n [ (1-tx)^(-1/t) - 1 ] evaluated at t=x=0. (1-tx)^(-1/t) - 1 is the e.g.f. for a plane m-ary tree when t=m-1. See Bergeron et al. in "Varieties of Increasing Trees". - Tom Copeland, Dec 09 2007
First comment and formula above rephrased as o.g.f. for row n: Product_{i=0...n} (1+i*x). - Geoffrey Critzer, Feb 04 2011
n-th row polynomials with alternate signs are the characteristic polynomials of the (n-1)x(n-1) matrices with 1's in the superdiagonal, (1,2,3,...) in the main diagonal, and the rest zeros. For example, the characteristic polynomial of [1,1,0; 0,2,1; 0,0,3] is x^3 - 6*x^2 + 11*x - 6. - Gary W. Adamson, Jun 28 2011
E.g.f.: A(x,y) = x*y/(1 - x*y)^(1 + 1/y) = Sum_{n>=1, k=1..n} T(n,k)*x^n*y^k/(n-1)!. - Paul D. Hanna, Jul 21 2011
With F(x,t) = (1-t*x)^(-1/t) - 1 an e.g.f. for the row polynomials P(n,t) of A094638 with P(0,t)=0, G(x,t)= [1-(1+x)^(-t)]/t is the comp. inverse in x. Consequently, with H(x,t) = 1/(dG(x,t)/dx) = (1+x)^(t+1),
  P(n,t) = [(H(x,t)*d/dx)^n] x evaluated at x=0; i.e.,
  F(x,t) = exp[x*P(.,t)] = exp[x*H(u,t)*d/du] u, evaluated at u = 0.
  Also, dF(x,t)/dx = H(F(x,t),t). - Tom Copeland, Sep 20 2011
T(n,k) = |A008276(n,k)|. - R. J. Mathar, May 19 2016
The row polynomials of this entry are the reversed row polynomials of A143491 multiplied by (1+x). E.g., (1+x)(1 + 5x + 6x^2) = (1 + 6x + 11x^2 + 6x^3). - Tom Copeland, Dec 11 2016
Regarding the row e.g.f.s in Copeland's 2007 formulas, e.g.f.s for A001710, A001715, and A001720 give the compositional inverses of the e.g.f. here for t = 2, 3, and 4 respectively. - Tom Copeland, Dec 28 2019

Edited by Emeric Deutsch, Aug 14 2006

A008548 Quintuple factorial numbers: Product_{k=0..n-1} (5*k+1).

Original entry on oeis.org

1, 1, 6, 66, 1056, 22176, 576576, 17873856, 643458816, 26381811456, 1213563326976, 61891729675776, 3465936861843456, 211422148572450816, 13953861805781753856, 990724188210504523776, 75295038303998343806976, 6098898102623865848365056, 524505236825652462959394816

Offset: 0

Joe Keane (jgk(AT)jgk.org)

a(n), n>=1, enumerates increasing sextic (6-ary) trees with n vertices. - Wolfdieter Lang, Sep 14 2007

Hankel transform is A169620. - Paul Barry, Dec 03 2009

Vincenzo Librandi, Table of n, a(n) for n = 0..300 (first 50 terms from T. D. Noe)
Martin Burtscher, Igor Szczyrba, Rafał Szczyrba, Analytic Representations of the n-anacci Constants and Generalizations Thereof, Journal of Integer Sequences, Vol. 18 (2015), Article 15.4.5.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), Article 00.2.4.
J.-C. Novelli and J.-Y. Thibon, Hopf Algebras of m-permutations,(m+1)-ary trees, and m-parking functions, arXiv preprint arXiv:1403.5962 [math.CO], 2014-2020.
Maxie D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications , J. Int. Seq. 13 (2010), Article 10.6.7, Table 6.3.

Cf. A001147, A007559, A007696, A016861, A034687, A034688, A052562, A047055, A049385, A051150.

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

[(&*[5*k+1: k in [0..n]]): n in [0..20]]; // G. C. Greubel, Aug 16 2019

a := n -> mul(5*k+1, k=0..n-1);
G(x):=(1-5*x)^(-1/5): f[0]:=G(x): for n from 1 to 29 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n],n=0..16); # Zerinvary Lajos, Apr 03 2009
H := hypergeom([1, 1/5], [], 5*x):
seq(coeff(series(H,x,20),x,n),n=0..16); # Peter Luschny, Oct 08 2015

Table[Product[5k+1,{k,0,n-1}],{n,0,20}]  (* Harvey P. Dale, Apr 23 2011 *)
FoldList[Times,1,NestList[#+5&,1,20]] (* Ray Chandler, Apr 23 2011 *)
FoldList[Times,1,5Range[0, 25] + 1] (* Vincenzo Librandi, Jun 10 2013 *)

x='x+O('x^33); Vec(serlaplace((1-5*x)^(-1/5))) \\ Joerg Arndt, Apr 24 2011

vector(20, n, n--; prod(k=0, n-1, 5*k+1)) \\ Altug Alkan, Oct 08 2015

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

a(n) = A049385(n, 1) (first column of triangle).
E.g.f.: (1-5*x)^(-1/5).
a(n) ~ 2^(1/2)*Pi^(1/2)*gamma(1/5)^-1*n^(-3/10)*5^n*e^-n*n^n*{1 + 1/300*n^-1 - ...}. - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
a(n) = Sum_{k=0..n} (-5)^(n-k)*A048994(n, k). - Philippe Deléham, Oct 29 2005
G.f.: 1/(1-x/(1-5x/(1-6x/(1-10x/(1-11x/(1-15x/(1-16x/(1-20x/(1-21x/(1-25x/(1-.../(1-A008851(n+1)*x/(1-... (continued fraction). - Paul Barry, Dec 03 2009
a(n)=(-4)^n*Sum_{k=0..n} (5/4)^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: 1/Q(0) where Q(k) = 1 - x*(5*k+1)/(1 - x*(5*k+5)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 20 2013
G.f.: G(0)/2, where G(k)= 1  + 1/(1 - (5*k+1)*x/((5*k+1)*x + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 14 2013
a(n) = (10n-18)*a(n-2) + (5n-6)*a(n-1), n>=2. - Ivan N. Ianakiev, Aug 12 2013
Let T(x) = 1/(1 - 4*x)^(1/4) be the e.g.f. for the sequence of triple factorial numbers A007696. Then the e.g.f. A(x) for the quintuple factorial numbers satisfies T( Integral_{t = 0..x} A(t) dt ) = A(x). Cf. A007559 and A007696. - Peter Bala, Jan 02 2015
O.g.f.: hypergeom([1, 1/5], [], 5*x). - Peter Luschny, Oct 08 2015
a(n) = 5^n * Gamma(n + 1/5) / Gamma(1/5). - Artur Jasinski, Aug 23 2016
D-finite with recurrence: a(n) +(-5*n+4)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
Sum_{n>=0} 1/a(n) = 1 + (e/5^4)^(1/5)*(Gamma(1/5) - Gamma(1/5, 1/5)). - Amiram Eldar, Dec 19 2022

A001715 a(n) = n!/6.

Original entry on oeis.org

1, 4, 20, 120, 840, 6720, 60480, 604800, 6652800, 79833600, 1037836800, 14529715200, 217945728000, 3487131648000, 59281238016000, 1067062284288000, 20274183401472000, 405483668029440000, 8515157028618240000, 187333454629601280000, 4308669456480829440000

Offset: 3

N. J. A. Sloane

The numbers (4, 20, 120, 840, 6720, ...) arise from the divisor values in the general formula a(n) = n*(n+1)*(n+2)*(n+3)* ... *(n+k)*(n*(n+k) + (k-1)*k/6)/((k+3)!/6) (which covers the following sequences: A000578, A000537, A024166, A101094, A101097, A101102). - Alexander R. Povolotsky, May 17 2008
a(n) is also the number of decreasing 3-cycles in the decomposition of permutations as product of disjoint cycles, a(3)=1, a(4)=4, a(5)=20. - Wenjin Woan, Dec 21 2008
Equals eigensequence of triangle A130128 reflected. - Gary W. Adamson, Dec 23 2008
a(n) is the number of n-permutations having 1, 2, and 3 in three distinct cycles. - Geoffrey Critzer, Apr 26 2009
From Johannes W. Meijer, Oct 20 2009: (Start)
The asymptotic expansion of the higher order exponential integral E(x,m=1,n=4) ~ exp(-x)/x*(1 - 4/x + 20/x^2 - 120/x^3 + 840/x^4 - 6720/x^5 + 60480/x^6 - 604800/x^7 + ...) leads to the sequence given above. See A163931 and A130534 for more information.
(End)

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

Vincenzo Librandi, Table of n, a(n) for n = 3..200
Somaya Barati, Beáta Bényi, Abbas Jafarzadeh, and Daniel Yaqubi, Mixed restricted Stirling numbers, arXiv:1812.02955 [math.CO], 2018.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 263.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), Article 00.2.4.
D. S. Mitrinovic and R. S. Mitrinovic, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz. No. 77 1962, 77 pp.
Alexsandar Petojevic, The Function vM_m(s; a; z) and Some Well-Known Sequences, Journal of Integer Sequences, Vol. 5 (2002), Article 02.1.7.
A. N. Stokes, Continued fraction solutions of the Riccati equation, Bull. Austral. Math. Soc. Vol. 25 (1982), 207-214.
Index entries for sequences related to factorial numbers.
Index to divisibility sequences.

Cf. A049352, A049458, A049460.
Cf. A034472, A130128.
Cf. A245334, A000142, A111530.
Cf. A001710, A001720, A007759, A094638.
Cf. A094587, A138533, A173333, A213936.

a001715 = (flip div 6) . a000142 -- Reinhard Zumkeller, Aug 31 2014

[Factorial(n)/6: n in [3..30]]; // Vincenzo Librandi, Jun 20 2011

f := proc(n) n!/6; end;
BB:= [S, {S = Prod(Z,Z,C), C = Union(B,Z,Z), B = Prod(Z,C)}, labelled]: seq(combstruct[count](BB, size=n)/12, n=3..20); # Zerinvary Lajos, Jun 19 2008
G(x):=1/(1-x)^4: f[0]:=G(x): for n from 1 to 18 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n],n=0..16); # Zerinvary Lajos, Apr 01 2009

a[n_]:=n!/6; (*Vladimir Joseph Stephan Orlovsky, Dec 13 2008 *)
Range[3,30]!/6 (* Harvey P. Dale, Aug 12 2012 *)

a(n)=n!/6 \\ Charles R Greathouse IV, Jan 12 2012

a(n) = A049352(n-2, 1) (first column of triangle).
E.g.f. if offset 0: 1/(1-x)^4.
a(n) = A173333(n,3). - Reinhard Zumkeller, Feb 19 2010
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x/(x + 1/(k+4)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 01 2013
G.f.: W(0), where W(k) = 1 - x*(k+4)/( x*(k+4) - 1/(1 - x*(k+1)/( x*(k+1) - 1/W(k+1) ))); (continued fraction). - Sergei N. Gladkovskii, Aug 26 2013
a(n) = A245334(n,n-3) / 4. - Reinhard Zumkeller, Aug 31 2014
From Peter Bala, May 22 2017: (Start)
The o.g.f. A(x) satisfies the Riccati equation x^2*A'(x) + (4*x - 1)*A(x) + 1 = 0.
G.f. as an S-fraction: A(x) = 1/(1 - 4*x/(1 - x/(1 - 5*x/(1 - 2*x/(1 - 6*x/(1 - 3*x/(1 - ... - (n + 3)*x/(1 - n*x/(1 - ... ))))))))) (apply Stokes, 1982).
A(x) = 1/(1 - 3*x - x/(1 - 4*x/(1 - 2*x/(1 - 5*x/(1 - 3*x/(1 - 6*x/(1 - ... - n*x/(1 - (n+3)*x/(1 - ... ))))))))). (End)
H(x) = (1 - (1 + x)^(-3)) / 3 = x - 4 x^2/2! + 20 x^3/3! - ... is an e.g.f. of the signed sequence (n!/4!), which is the compositional inverse of G(x) = (1 - 3*x)^(-1/3) - 1, an e.g.f. for A007559. Cf. A094638, A001710 (for n!/2!), and A001720 (for n!/4!). Cf. columns of A094587, A173333, and A213936 and rows of A138533.- Tom Copeland, Dec 27 2019
E.g.f.: x^3 / (3! * (1 - x)). - Ilya Gutkovskiy, Jul 09 2021
From Amiram Eldar, Jan 15 2023: (Start)
Sum_{n>=3} 1/a(n) = 6*e - 15.
Sum_{n>=3} (-1)^(n+1)/a(n) = 3 - 6/e. (End)

More terms from Harvey P. Dale, Aug 12 2012

A032031 Triple factorial numbers: (3n)!!! = 3^n*n!.

Original entry on oeis.org

1, 3, 18, 162, 1944, 29160, 524880, 11022480, 264539520, 7142567040, 214277011200, 7071141369600, 254561089305600, 9927882482918400, 416971064282572800, 18763697892715776000, 900657498850357248000, 45933532441368219648000, 2480410751833883860992000

Offset: 0

Christian G. Bower

For n >= 1 a(n) is the order of the wreath product of the symmetric group S_n and the elementary Abelian group (C_3)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001
Laguerre transform of double factorials 2^n*n! = A000165(n). - Paul Barry, Aug 08 2008
For positive n, a(n) equals the permanent of the n X n matrix consisting entirely of 3's. - John M. Campbell, May 26 2011
a(n) is the product of the positive integers <= 3*n that are multiples of 3. - Peter Luschny, Jun 23 2011
Partial products of A008585. - Reinhard Zumkeller, Sep 20 2013

Vincenzo Librandi, Table of n, a(n) for n = 0..100
CombOS - Combinatorial Object Server, Generate colored permutations
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 491
Alexsandar Petojevic, The Function vM_m(s; a; z) and Some Well-Known Sequences, Journal of Integer Sequences, Vol. 5 (2002), Article 02.1.7.
Michael Z. Spivey and Laura L. Steil, The k-Binomial Transforms and the Hankel Transform, Journal of Integer Sequences, Vol. 9 (2006), Article 06.1.1.

Cf. A000142, A007559, A008544, A051141, A000165, A092041, A092615.

Cf. Subsequence of A007661.

a032031 n = a032031_list !! n
a032031_list = scanl (*) 1 $ tail a008585_list
-- Reinhard Zumkeller, Sep 20 2013

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

with(combstruct):ZL:=[T,{T=Union(Z,Prod(Epsilon,Z,T), Prod(T,Z,Epsilon), Prod(T,Z))},labeled]:seq(count(ZL,size=i)/i,i=1..17); # Zerinvary Lajos, Dec 16 2007
A032031 := n -> mul(k, k = select(k-> k mod 3 = 0, [$1 .. 3*n])): seq(A032031(n), n = 0 .. 16); # Peter Luschny, Jun 23 2011

Table[3^n*Gamma[1 + n], {n, 0, 20}] (* Roger L. Bagula, Oct 30 2008 *)
Join[{1},FoldList[Times,3*Range[20]]] (* Harvey P. Dale, Feb 10 2019 *)
Table[Times@@Range[3n,1,-3],{n,0,20}] (* Harvey P. Dale, Apr 14 2023 *)

PARI
```
a(n)=3^n*n!;
```
PARI
```
a(n)=prod(k=1,n, 3*k );
```

def A032031(n) : return mul(j for j in range(3,3*(n+1),3))
[A032031(n) for n in (0..16)]  # Peter Luschny, May 20 2013

a(n) = 3^n*n!.
a(n) = Product_{k=1..n} 3*k.
E.g.f.: 1/(1-3*x).
a(n) = Sum_{k=0..n} C(n,k)*(n!/k!)*2^k*k!. - Paul Barry, Aug 08 2008
a(0) = 1, a(n) = 3*n*a(n-1). - Arkadiusz Wesolowski, Oct 04 2011
G.f.: 2/G(0), where G(k)= 1 + 1/(1 - 6*x*(k+1)/(6*x*(k+1) - 1 + 6*x*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
G.f.: G(0)/2, where G(k)= 1 + 1/(1 - x*(3*k+3)/(x*(3*k+3) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 06 2013
G.f.: 1/Q(0), where Q(k) = 1 - 3*x*(2*k+1) - 9*x^2*(k+1)^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Sep 28 2013
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/3) (A092041).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/3) (A092615). (End)

A004987 a(n) = (3^n/n!)Product_{k=0..n-1} (3k + 1). 3-central binomial coefficients.

Original entry on oeis.org

1, 3, 18, 126, 945, 7371, 58968, 480168, 3961386, 33011550, 277297020, 2344420260, 19927572210, 170150808870, 1458435504600, 12542545339560, 108179453553705, 935434098376155, 8107095519260010, 70403724246205350, 612512400941986545, 5337608065351597035, 46582761297613937760

Offset: 0

Joe Keane (jgk(AT)jgk.org)

Diagonal of rational function R(x,y) = (1 - 9*x*y) / (1 - 2*x - 3*y + 3*y^2 + 9*x^2*y). - Gheorghe Coserea, Jul 01 2016

This is the k = 3 variant of the k-central binomial coefficients c(n,k) with g.f. (1 - k^2*x)^(-1/k), which yield the usual central binomial coefficients A001405 for k = 2. - M. F. Hasler, Nov 12 2024

			G.f.: 1 + 3*x + 18*x^2 + 126*x^3 + 945*x^4 + 7371*x^5 + 58968*x^6 + 480168*x^7 + ...

Gheorghe Coserea, Table of n, a(n) for n = 0..200
A. Bostan, S. Boukraa, J.-M. Maillard, and J.-A. Weil, Diagonals of rational functions and selected differential Galois groups, arXiv preprint arXiv:1507.03227 [math-ph], 2015.
A. Straub, V. H. Moll, and T. Amdeberhan, The p-adic valuation of k-central binomial coefficients, Acta Arith. 140 (2009) 31-41, eq (1.10).

Cf. A004117, A007559, A047657, A004988, A054861, A034689, A053101, A072888, A216316.
Related to diagonal of rational functions: A268545-A268555.
Cf. A008931, A078532, A245114, A247029, A376282, A376636.

List([0..25], n-> 3^n*Product([0..n-1], k-> 3*k+1)/Factorial(n) ); # G. C. Greubel, Aug 22 2019

[1] cat [3^n*&*[3*k+1: k in [0..n-1]]/Factorial(n): n in [1..25]]; // G. C. Greubel, Aug 22 2019

a:= n-> (3^n/n!)*mul(3*k+1, k=0..n-1); seq(a(n), n=0..25); # G. C. Greubel, Aug 22 2019

Table[(-9)^n Binomial[-1/3, n], {n, 0, 25}] (* Jean-François Alcover, Sep 28 2016, after Peter Luschny *)

a(n) = prod(k=0, n-1, 3*k + 1)*3^n/n! \\ Michel Marcus, Jun 30 2013

my(x='x, y='y);
R = (1 - 9*x*y) / (1 - 2*x - 3*y + 3*y^2 + 9*x^2*y);
diag(n, expr, var) = {
  my(a = vector(n));
  for (i = 1, #var, expr = taylor(expr, var[#var - i + 1], n));
  for (k = 1, n, a[k] = expr;
       for (i = 1, #var, a[k] = polcoeff(a[k], k-1)));
  return(a);
};
diag(20, R, [x,y])  \\ Gheorghe Coserea, Jul 01 2016

Vec((1-9*x+O(x^25))^(-1/3)) \\ yields the same as:
apply( {A004987(n)=prod(k=0, n-1, 9*k+3)\n!}, [0..24]) \\ M. F. Hasler, Nov 12 2024

[9^n*rising_factorial(1/3, n)/factorial(n) for n in (0..25)] # G. C. Greubel, Aug 22 2019

G.f.: (1 - 9*x)^(-1/3).
a(n) = (3^n/n!)*A007559(n), n >= 1, a(0) := 1.
a(n) ~ Gamma(1/3)^-1*n^(-2/3)*3^(2*n)*{1 - 1/9*n^-1 + ...}.
Representation as n-th moment of a positive function on (0, 9): a(n) = Integral_{x=0..9} ( x^n*(1/(Pi*sqrt(3)*6*(x/9)^(2/3)*(1-x/9)^(1/3))) ), n >= 0. This function is the solution of the Hausdorff moment problem on (0, 9) with moments equal to a(n). As a consequence this representation is unique. - Karol A. Penson, Jan 30 2003
D-finite with recurrence: n*a(n) + 3*(2-3*n)*a(n-1)=0. - R. J. Mathar, Jun 07 2013
0 = a(n) * (81*a(n+1) - 15*a(n+2)) + a(n+1) * (-3*a(n+1) + a(n+2)) for all n in Z. - Michael Somos, Jan 27 2014
G.f. A(x)=:y satisfies 0 = y'' * y - 4 * y' * y'. - Michael Somos, Jan 27 2014
a(n) = (-9)^n*binomial(-1/3, n). - Peter Luschny, Mar 23 2014
E.g.f.: is the hypergeometric function of type 1F1, in Maple notation hypergeom([1/3], [1], 9*x). - Karol A. Penson, Dec 19 2015
Sum_{n>=0} 1/a(n) = (sqrt(3)*Pi + 3*(12 + log(3)))/32 = 1.3980385924595932... - Ilya Gutkovskiy, Jul 01 2016
Binomial transform of A216316. - Peter Bala, Jul 02 2023
From Peter Bala, Mar 31 2024: (Start)
a(n) = (9^n)*Sum_{k = 0..2*n} (-1)^k*binomial(-1/3, k)* binomial(-1/3, 2*n - k).
(9^n)*a(n) = Sum_{k = 0..2*n} (-1)^k*a(k)*a(2*n-k).
Sum_{k = 0..n} a(k)*a(n-k) = A004988(n).
Sum_{k = 0..2*n} a(k)*a(2*n-k) = 18^n/(2*n)! * Product_{k = 1..n} (6*k - 1)*(3*k - 2). (End)
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^3). - Seiichi Manyama, Jun 20 2025

More terms from Ralf Stephan, Mar 13 2004

More terms from Benoit Cloitre, Jun 05 2004

cp's OEIS Frontend

A048994 Triangle of Stirling numbers of first kind, s(n,k), n >= 0, 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

PARI

PARI

Formula

Extensions

A132393 Triangle of unsigned Stirling numbers of the first kind (see A048994), read by rows, T(n,k) for 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

PARI

Formula

A007661 Triple factorial numbers a(n) = n!!!, defined by a(n) = n*a(n-3), a(0) = a(1) = 1, a(2) = 2. Sometimes written n!3.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

Sage

Formula

A007696 Quartic (or 4-fold) factorial numbers: a(n) = Product_{k = 0..n-1} (4*k + 1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

Maxima

PARI

Sage

Formula

Extensions

A008544 Triple factorial numbers: Product_{k=0..n-1} (3*k+2).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs