A032031 -id:A032031 - cp's OEIS frontend

A004747 Triangle read by rows: the Bell transform of the triple factorial numbers A008544 without column 0.

1, 2, 1, 10, 6, 1, 80, 52, 12, 1, 880, 600, 160, 20, 1, 12320, 8680, 2520, 380, 30, 1, 209440, 151200, 46480, 7840, 770, 42, 1, 4188800, 3082240, 987840, 179760, 20160, 1400, 56, 1, 96342400, 71998080, 23826880, 4583040, 562800, 45360, 2352, 72, 1

Offset: 1

Wolfdieter Lang

Previous name was: Triangle of numbers related to triangle A048966; generalization of Stirling numbers of second kind A008277, Bessel triangle A001497.
T(n,m) = S2p(-2; n,m), a member of a sequence of triangles including S2p(-1; n,m) = A001497(n-1,m-1) (Bessel triangle) and ((-1)^(n-m))*S2p(1; n,m) = A008277(n, m) (Stirling 2nd kind). T(n,1)= A008544(n-1).
T(n,m), n>=m>=1, enumerates unordered n-vertex m-forests composed of m plane (aka ordered) increasing (rooted) trees where vertices of out-degree r>=0 come in r+1 different types (like an (r+1)-ary vertex). Proof from the e.g.f. of the first column Y(z) = 1 - (1-3*x)^(1/3) and the F. Bergeron et al. eq. (8) Y'(z)= phi(Y(z)), Y(0) = 0, with out-degree o.g.f. phi(w)=1/(1-w)^2. - Wolfdieter Lang, Oct 12 2007
Also the Bell transform of the triple factorial numbers A008544 which adds a first column (1,0,0 ...) on the left side of the triangle. For the definition of the Bell transform see A264428. See A051141 for the triple factorial numbers A032031 and A203412 for the triple factorial numbers A007559 as well as A039683 and A132062 for the case of double factorial numbers. - Peter Luschny, Dec 21 2015

			Triangle begins:
       1;
       2,      1;
      10,      6,     1;
      80,     52,    12,    1;
     880,    600,   160,   20,   1;
   12320,   8680,  2520,  380,  30,  1;
  209440, 151200, 46480, 7840, 770, 42, 1;
Tree combinatorics for T(3,2)=6: Consider first the unordered forest of m=2 plane trees with n=3 vertices, namely one vertex with out-degree r=0 (root) and two different trees with two vertices (one root with out-degree r=1 and a leaf with r=0). The 6 increasing labelings come then from the forest with rooted (x) trees x, o-x (1,(3,2)), (2,(3,1)) and (3,(2,1)) and similarly from the second forest x, x-o (1,(2,3)), (2,(1,3)) and (3,(1,2)).

G. C. Greubel, Rows n = 1..50 of the triangle, flattened
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.
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
Richell O. Celeste, Roberto B. Corcino, and Ken Joffaniel M. Gonzales. Two Approaches to Normal Order Coefficients, Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) #09.8.3.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
Wolfdieter Lang, Combinatorial Interpretation of Generalized Stirling Numbers, J. Int. Seqs. Vol. 12 (2009) #09.3.3.
Mathias Pétréolle and Alan D. Sokal, Lattice paths and branched continued fractions. II. Multivariate Lah polynomials and Lah symmetric functions, arXiv:1907.02645 [math.CO], 2019.
Index entries for sequences related to Bessel functions or polynomials

Cf. A015735 (row sums).
Cf. A007559, A008277, A008544, A032031, A039683.
Cf. A048966, A051141, A132062, A203412, A264428.
Triangles with the recurrence T(n,k) = (m*(n-1)-k)*T(n-1,k) + T(n-1,k-1): A010054 (m=1), A001497 (m=2), this sequence (m=3), A000369 (m=4), A011801 (m=5), A013988 (m=6).

function T(n,k) // T = A004747
  if k eq 0 then return 0;
  elif k eq n then return 1;
  else return (3*(n-1)-k)*T(n-1,k) + T(n-1,k-1);
  end if;
end function;
[T(n,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Oct 03 2023

T := (n, m) -> 3^n/m!*(1/3*m*GAMMA(n-1/3)*hypergeom([1-1/3*m, 2/3-1/3*m, 1/3-1/3*m], [2/3, 4/3-n], 1)/GAMMA(2/3)-1/6*m*(m-1)*GAMMA(n-2/3)*hypergeom( [1-1/3*m, 2/3-1/3*m, 4/3-1/3*m], [4/3, 5/3-n], 1)/Pi*3^(1/2)*GAMMA(2/3)):
for n from 1 to 6 do seq(simplify(T(n,k)),k=1..n) od;
# Karol A. Penson, Feb 06 2004
# The function BellMatrix is defined in A264428.
# Adds (1,0,0,0, ..) as column 0.
BellMatrix(n -> mul(3*k+2, k=(0..n-1)), 9); # Peter Luschny, Jan 29 2016

(* First program *)
T[1,1]= 1; T[, 0]= 0; T[0, ]= 0; T[n_, m_]:= (3*(n-1)-m)*T[n-1, m]+T[n-1, m-1];
Flatten[Table[T[n, m], {n,12}, {m,n}] ][[1 ;; 45]] (* Jean-François Alcover, Jun 16 2011, after recurrence *)
(* Second program *)
f[n_, m_]:= m/n Sum[Binomial[k, n-m-k] 3^k (-1)^(n-m-k) Binomial[n+k-1, n-1], {k, 0, n-m}]; Table[n! f[n, m]/(m! 3^(n-m)), {n,12}, {m,n}]//Flatten (* Michael De Vlieger, Dec 23 2015 *)
(* Third program *)
rows = 12;
T[n_, m_]:= BellY[n, m, Table[Product[3k+2, {k, 0, j-1}], {j, 0, rows}]];
Table[T[n, m], {n,rows}, {m,n}]//Flatten (* Jean-François Alcover, Jun 22 2018 *)

# uses [bell_transform from A264428]
triplefactorial = lambda n: prod(3*k+2 for k in (0..n-1))
def A004747_row(n):
    trifact = [triplefactorial(k) for k in (0..n)]
    return bell_transform(n, trifact)
[A004747_row(n) for n in (0..10)] # Peter Luschny, Dec 21 2015

T(n, m) = n!*A048966(n, m)/(m!*3^(n-m));

T(n+1, m) = (3*n-m)*T(n, m)+ T(n, m-1), for n >= m >= 1, with T(n, m) = 0, for n

E.g.f. of m-th column: ( 1 - (1-3*x)^(1/3) )^m/m!.

Sum_{k=1..n} T(n, k) = A015735(n).

For a formula expressed as special values of hypergeometric functions 3F2 see the Maple program below. - Karol A. Penson, Feb 06 2004

T(n,1) = A008544(n-1). - Peter Luschny, Dec 23 2015

New name from Peter Luschny, Dec 21 2015

A051141 Triangle read by rows: a(n, m) = S1(n, m)*3^(n-m), where S1 are the signed Stirling numbers of first kind A008275 (n >= 1, 1 <= m <= n).

Original entry on oeis.org

1, -3, 1, 18, -9, 1, -162, 99, -18, 1, 1944, -1350, 315, -30, 1, -29160, 22194, -6075, 765, -45, 1, 524880, -428652, 131544, -19845, 1575, -63, 1, -11022480, 9526572, -3191076, 548289, -52920, 2898, -84, 1, 264539520, -239660208

Offset: 1

Wolfdieter Lang

Previous name was: Generalized Stirling number triangle of first kind.
a(n,m) = R_n^m(a=0,b=3) in the notation of the given reference.
a(n,m) is a Jabotinsky matrix, i.e., the monic row polynomials E(n,x) := Sum_{m=1..n} a(n,m)*x^m = Product_{j=0..n-1} (x - 3*j), n >= 1 and E(0,x) := 1 are exponential convolution polynomials (see A039692 for the definition and a Knuth reference).
This is the signed Stirling1 triangle with diagonals d>=0 (main diagonal d=0) scaled with 3^d.
Exponential Riordan array [1/(1 + 3*x), log(1 + 3*x)/3]. The unsigned triangle is [1/(1 - 3*x), log(1/(1 - 3*x)^(1/3))]. - Paul Barry, Apr 29 2009
Also the Bell transform of the triple factorial numbers A032031 which adds a first column (1, 0, 0 ...) on the left side of the triangle and computes the unsigned values. For the definition of the Bell transform, see A264428. See A004747 for the triple factorial numbers A008544 and A203412 for the triple factorial numbers A007559 as well as A039683 and A132062 for the case of double factorial numbers. - Peter Luschny, Dec 21 2015

			Triangle starts:
       1;
      -3,       1;
      18,      -9,      1;
    -162,      99,    -18,      1;
    1944,   -1350,    315,    -30,    1;
  -29160,   22194,  -6075,    765,  -45,   1;
  524880, -428652, 131544, -19845, 1575, -63, 1;
---
Row polynomial E(3,x) = 18*x-9*x^2+x^3.
From _Paul Barry_, Apr 29 2009: (Start)
The unsigned array [1/(1 - 3*x), log(1/(1 - 3*x)^(1/3))] has production matrix
    3,    1;
    9,    6,    1;
   27,   27,    9,   1;
   81,  108,   54,  12,   1;
  243,  405,  270,  90,  15,  1;
  729, 1458, 1215, 540, 135, 18, 1;
  ...
which is A007318^{3} beheaded (by viewing A007318 as a lower triangular matrix). See the comment above. (End)

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
Richell O. Celeste, Roberto B. Corcino, and Ken Joffaniel M. Gonzales, Two Approaches to Normal Order Coefficients, Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
Wolfdieter Lang, First 10 rows.
D. S. Mitrinovic and M. S. Mitrinovic, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz. 77 (1962), 1-77.

First (m=1) column sequence is: A032031(n-1).
Row sums (signed triangle): A008544(n-1)*(-1)^(n-1).
Row sums (unsigned triangle): A007559(n).
Cf. A008275 (Stirling1 triangle, b=1), A039683 (b=2), A051142 (b=4).
Cf. A039683, A132062, A264428.

a[n_, m_] /; n >= m >= 1 := a[n, m] = a[n-1, m-1] - 3(n-1)*a[n-1, m]; a[n_, m_] /; n < m = 0; a[, 0] = 0; a[1, 1] = 1; Flatten[Table[a[n, m], {n, 1, 9}, {m, 1, n}]][[1 ;; 38]] (* _Jean-François Alcover, Jun 01 2011, after formula *)
Table[StirlingS1[n, m]*3^(n - m), {n, 1, 10}, {m, 1, n}]//Flatten (* G. C. Greubel, Oct 24 2017 *)

for(n=1,10, for(m=1,n, print1(stirling(n,m,1)*3^(n-m), ", "))) \\ G. C. Greubel, Oct 24 2017

# uses[bell_transform from A264428]
triplefactorial = lambda n: 3^n*factorial(n)
def A051141_row(n):
    trifact = [triplefactorial(k) for k in (0..n)]
    return bell_transform(n, trifact)
[A051141_row(n) for n in (0..8)] # Peter Luschny, Dec 21 2015

a(n, m) = a(n-1, m-1) - 3*(n-1)*a(n-1, m) for n >= m >= 1; a(n, m) = 0 for n < m; a(n, 0) = 0 for n >= 1; a(0, 0) = 1.
E.g.f. for the m-th column of the signed triangle: (log(1 + 3*x)/3)^m/m!.
|a(n,1)| = A032031(n-1). - Peter Luschny, Dec 23 2015

Name clarified using a formula of the author by Peter Luschny, Dec 23 2015

A114799 Septuple factorial, 7-factorial, n!7, n!!!!!!!, a(n) = n*a(n-7) if n > 1, else 1.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 7, 8, 18, 30, 44, 60, 78, 98, 120, 288, 510, 792, 1140, 1560, 2058, 2640, 6624, 12240, 19800, 29640, 42120, 57624, 76560, 198720, 379440, 633600, 978120, 1432080, 2016840, 2756160, 7352640, 14418720, 24710400, 39124800

Offset: 0

Jonathan Vos Post, Feb 18 2006

Many of the terms yield multifactorial primes a(n) + 1, e.g.: a(2) + 1 = 3, a(4) + 1 = 5, a(6) + 1 = 7, a(9) + 1 = 19, a(10) + 1 = 31, a(12) + 1 = 61, a(13) + 1 = 79, a(24) + 1 = 12241, a(25) + 1 = 19801, a(26) + 1 = 29641, a(29) + 1 = 76561, a(31) + 1 = 379441, a(35) + 1 = 2016841, a(36) + 1 = 2756161, ...

Equivalently, product of all positive integers <= n congruent to n (mod 7). - M. F. Hasler, Feb 23 2018

			a(40) = 40 * a(40-7) = 40 * a(33) = 40 * (33*a(26)) = 40 * 33 * (26*a(19)) = 40 * 33 * 26 * (19*a(12)) = 40 * 33 * 26 * 19 * (12*a(5)) = 40 * 33 * 26 * 19 * 12 5 = 39124800.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Multifactorial.
Index entries for sequences related to factorial numbers

Cf. A045754, A084947, A288094.

Cf. k-fold factorials: A000142, A001147 (and A000165, A006882), A007559 (and A032031, A008544, A007661), A007696 (and A001813, A008545, A047053, A007662), A008548 (and A052562, A047055, A085157), A085158 (and A008542, A047058, A047657), A045755.

a:= function(n)
    if n<1 then return 1;
    else return n*a(n-7);
    fi;
  end;
List([0..40], n-> a(n) ); # G. C. Greubel, Aug 20 2019

b:= func< n | (n lt 8) select n else n*Self(n-7) >;
[1] cat [b(n): n in [1..40]]; // G. C. Greubel, Aug 20 2019

A114799 := proc(n)
    option remember;
    if n < 1 then
        1;
    else
        n*procname(n-7) ;
    end if;
end proc:
seq(A114799(n),n=0..40) ; # R. J. Mathar, Jun 23 2014
A114799 := n -> product(n-7*k,k=0..(n-1)/7); # M. F. Hasler, Feb 23 2018

a[n_]:= If[n<1, 1, n*a[n-7]]; Table[a[n], {n,0,40}] (* G. C. Greubel, Aug 20 2019 *)

A114799(n,k=7)=prod(j=0,(n-1)\k,n-j*k) \\ M. F. Hasler, Feb 23 2018

def a(n):
    if (n<1): return 1
    else: return n*a(n-7)
[a(n) for n in (0..40)] # G. C. Greubel, Aug 20 2019

a(n) = 1 for n <= 1, else a(n) = n*a(n-7).

Sum_{n>=0} 1/a(n) = A288094. - Amiram Eldar, Nov 10 2020

Edited by M. F. Hasler, Feb 23 2018

A144739 7-factorial numbers A114799(7n+3): Partial products of A017017(k) = 7k+3, a(0) = 1.

Original entry on oeis.org

1, 3, 30, 510, 12240, 379440, 14418720, 648842400, 33739804800, 1990648483200, 131382799891200, 9590944392057600, 767275551364608000, 66752972968720896000, 6274779459059764224000, 633752725365036186624000, 68445294339423908155392000, 7871208849033749437870080000

Offset: 0

Philippe Deléham, Sep 20 2008

nonn

			a(0)=1, a(1)=3, a(2)=3*10=30, a(3)=3*10*17=510, a(4)=3*10*17*24=12240, ...

G. C. Greubel, Table of n, a(n) for n = 0..335
Index entries for sequences related to factorial numbers.

Cf. A114799, A001710, A001147, A032031, A008545, A047056, A011781, A045754, A084947, A144827, A147585, A049209, A051188.

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

[ 1 ] cat [ &*[ (7*k+3): k in [0..n] ]: n in [0..20] ]; // Klaus Brockhaus, Nov 10 2008

a:= n-> product(7*j+3, j=0..n-1); seq(a(n), n=0..20); # G. C. Greubel, Aug 19 2019

Table[7^n*Pochhammer[3/7, n], {n,0,20}] (* G. C. Greubel, Aug 19 2019 *)

a(n)=prod(i=1,n,7*i-4) \\ Charles R Greathouse IV, Jul 02 2013

[product(7*k+3 for k in (0..n-1)) for n in (0..20)] # G. C. Greubel, Aug 19 2019

a(n) = Sum_{k=0..n} A132393(n,k)*3^k*7^(n-k).
G.f.: 1/(1-3*x/(1-7*x/(1-10*x/(1-14*x/(1-17*x/(1-21*x/(1-24*x/(1-... (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-4)^n*Sum_{k=0..n} (7/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
From Ilya Gutkovskiy, Mar 23 2017: (Start)
E.g.f.: 1/(1 - 7*x)^(3/7).
a(n) ~ sqrt(2*Pi)*7^n*n^n/(exp(n)*n^(1/14)*Gamma(3/7)). (End)
a(n) = A114799(7*n-4). - M. F. Hasler, Feb 23 2018
D-finite with recurrence: a(n) +(-7*n+4)*a(n-1)=0. - R. J. Mathar, Feb 21 2020
Sum_{n>=0} 1/a(n) = 1 + (e/7^4)^(1/7)*(Gamma(3/7) - Gamma(3/7, 1/7)). - Amiram Eldar, Dec 19 2022

A010845 a(n) = 3na(n-1) + 1, a(0) = 1.

Original entry on oeis.org

1, 4, 25, 226, 2713, 40696, 732529, 15383110, 369194641, 9968255308, 299047659241, 9868572754954, 355268619178345, 13855476147955456, 581929998214129153, 26186849919635811886, 1256968796142518970529

Offset: 0

Simon Plouffe

a(n)/(A000142*A000244) is an increasingly good approximation to cube root of e.
Related to Incomplete Gamma Function at 1/3. - Michael Somos, Mar 26 1999
For positive n, a(n) equals 3^n times the permanent of the n X n matrix with (4/3)'s along the main diagonal, and 1's everywhere else. - John M. Campbell, Jul 10 2011

			1 + 4*x + 25*x^2 + 226*x^3 + 2713*x^4 + 40696*x^5 + 732529*x^6 + ...

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, p. 262.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, p. 262.
Roland Bacher, Counting Packings of Generic Subsets in Finite Groups, Electr. J. Combinatorics, 19 (2012), #P7. - From _N. J. A. Sloane_, Feb 06 2013
M. Z. Spivey and L. L. Steil, The k-Binomial Transforms and the Hankel Transform, J. Integ. Seqs. Vol. 9 (2006), #06.1.1.

Cf. A000522, A010844, A056545, A056546, A056547 for analogs.

Table[ Gamma[ n, 1/3 ]*Exp[ 1/3 ]*3^(n-1), {n, 1, 24} ]
a[ n_] := If[ n<0, 0, Floor[ n! E^(1/3) 3^n ]] (* Michael Somos, Sep 04 2013 *)
Range[0, 20]! CoefficientList[Series[Exp[x]/(1 - 3 x), {x, 0, 20}], x] (* Vincenzo Librandi, Feb 17 2014 *)

{a(n) = if( n<0, 0, n! * sum(k=0, n, 3^(n-k) / k!))} /* Michael Somos, Sep 04 2013 */

E.g.f.: exp(x)/(1-3*x).
a(n) = floor( n!*e^(1/3)*3^n ) = n! * (Sum_{k=0..n} 3^(n-k) / k!) = n! * (e^(1/3) * 3^n - Sum_{k>n} 3^(n-k) / k!). - Michael Somos, Mar 26 1999
a(n) = Sum_{k=0..n} P(n, k)*3^k. - Ross La Haye, Aug 29 2005
Binomial transform of A032031. - Carl Najafi, Sep 11 2011
Conjecture: a(n) +(-3*n-1)*a(n-1) +3*(n-1)*a(n-2)=0. - R. J. Mathar, Feb 16 2014
a(n) = hypergeometric_U(1,n+2,1/3)/3. - Peter Luschny, Nov 26 2014
From Peter Bala, Mar 01 2017: (Start)
a(n) = Integral_{x >= 0} (3*x + 1)^n*exp(-x) dx.
The e.g.f. y = exp(x)/(1 - 3*x) satisfies the differential equation (1 - 3*x)*y' = (4 - 3*x)*y. Mathar's recurrence above follows easily from this.
The sequence b(n) := (3^n)*n! also satisfies Mathar's recurrence with b(0) = 1, b(1) = 3. This leads to the continued fraction representation a(n) = (3^n)*n!*( 1 + 1/(3 - 3/(7 - 6/(10 - ... - (3*n - 3)/(3*n + 1) )))) for n >= 2. Taking the limit as n -> oo gives the continued fraction representation exp(1/3) = 1 + 1/(3 - 3/(7 - 6/(10 - ... - (3*n - 3)/((3*n + 1) - ... )))). Cf. A010844. (End)

Better description and formulas from Michael Somos

More terms from James Sellers, Jul 04 2000

A225470 Triangle read by rows, s_3(n, k) where s_m(n, k) are the Stirling-Frobenius cycle numbers of order m; n >= 0, k >= 0.

Original entry on oeis.org

1, 2, 1, 10, 7, 1, 80, 66, 15, 1, 880, 806, 231, 26, 1, 12320, 12164, 4040, 595, 40, 1, 209440, 219108, 80844, 14155, 1275, 57, 1, 4188800, 4591600, 1835988, 363944, 39655, 2415, 77, 1, 96342400, 109795600, 46819324, 10206700, 1276009, 95200, 4186, 100, 1

Offset: 0

Peter Luschny, May 16 2013

The Stirling-Frobenius subset numbers S_{m}(n,k), for m >= 1 fixed, regarded as an infinite lower triangular matrix, can be inverted by Sum_{k} S_{m}(n,k)*s_{m}(k,j)*(-1)^(n-k) = [j = n]. The inverse numbers s_{m}(k,j), which are unsigned, are the Stirling-Frobenius cycle numbers. For m = 1 this gives the classical Stirling cycle numbers A132393. The Stirling-Frobenius subset numbers are defined in A225468.

Triangle T(n,k), read by rows, given by (2, 3, 5, 6, 8, 9, 11, 12, 14, 15, ... (A007494)) DELTA (1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, May 14 2015

			Triangle starts:
  [n\k][    0,      1,     2,     3,    4,  5,  6]
  [0]       1,
  [1]       2,      1,
  [2]      10,      7,     1,
  [3]      80,     66,    15,     1,
  [4]     880,    806,   231,    26,    1,
  [5]   12320,  12164,  4040,   595,   40,  1,
  [6]  209440, 219108, 80844, 14155, 1275, 57,  1.
  ...
From _Wolfdieter Lang_, Aug 11 2017: (Start)
Recurrence (see Maple program): T(4, 2) = T(3, 1) + (3*4 - 1)*T(3, 2) = 66 + 11*15 = 231.
Boas-Buck type recurrence for column k = 2 and n = 4: T(4, 2) = (4!/2)*(3*(2 + 6*(5/12))*T(2, 2)/2! + 1*(2 + 6*(1/2))*T(3,2)/3!) = (4!/2)*(3*9/4 + 5*15/3!) = 231. (End)

Peter Bala, A 3 parameter family of generalized Stirling numbers
Peter Luschny, Generalized Eulerian polynomials
Peter Luschny, The Stirling-Frobenius numbers

Cf. A225468; A132393 (m=1), A028338 (m=2), A225471(m=4).
Column k=0..4 give A008544, A024395(n-1), A286722(n-2), A383221, A383222.
T(n, n) ~ A000012; T(n, n-1) ~ A005449; T(n, n-2) ~ A024391; T(n, n-3) ~ A024392.
row sums ~ A032031; alternating row sums ~ A007559.
Cf. A132393.

SF_C := proc(n, k, m) option remember;
if n = 0 and k = 0 then return(1) fi;
if k > n or  k < 0 then return(0) fi;
SF_C(n-1, k-1, m) + (m*n-1)*SF_C(n-1, k, m) end:
seq(print(seq(SF_C(n, k, 3), k = 0..n)), n = 0..8);

SFC[0, 0, ] = 1; SFC[n, k_, ] /; (k > n || k < 0) = 0; SFC[n, k_, m_] := SFC[n, k, m] = SFC[n-1, k-1, m] + (m*n-1)*SFC[n-1, k, m]; Table[SFC[n, k, 3], {n, 0, 8}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 26 2013, after Maple *)

For a recurrence see the Maple program.
From Wolfdieter Lang, May 18 2017: (Start)
This is the Sheffer triangle (1/(1 - 3*x)^{-2/3}, -(1/3)*log(1-3*x)). See the P. Bala link where this is called exponential Riordan array, and the signed version is  denoted by s_{(3,0,2)}.
E.g.f. of row polynomials in the variable x (i.e., of the triangle): (1 - 3*z)^{-(2+x)/3}.
E.g.f. of column k: (1-3*x)^(-2/3)*(-(1/3)*log(1-3*x))^k/k!, k >= 0.
Recurrence for row polynomials R(n, x) = Sum_{k=0..n} T(n, k)*x^k: R(n, x) = (x+2)*R(n-1,x+3), with R(0, x) = 1.
R(n, x) = risefac(3,2;x,n) := Product_{j=0..(n-1)} (x + (2 + 3*j)). (See the P. Bala link, eq. (16) for the signed s_{3,0,2} row polynomials.)
T(n, k) = Sum_{j=0..(n-m)} binomial(n-j, k)* S1p(n, n-j)*2^(n-k-j)*3^j, with S1p(n, m) = A132393(n, m). (End)
Boas-Buck type recurrence for column sequence k: T(n, k) = (n!/(n - k)) * Sum_{p=k..n-1} 3^(n-1-p)*(2 + 3*k*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), beginning {1/2, 5/12, 3/8, 251/720, ...}. See a comment and references in A286718. - Wolfdieter Lang, Aug 11 2017

A051604 a(n) = (3*n+4)!!!/4!!!.

Original entry on oeis.org

1, 7, 70, 910, 14560, 276640, 6086080, 152152000, 4260256000, 132067936000, 4490309824000, 166141463488000, 6645658539520000, 285763317199360000, 13145112591170560000, 644110516967357440000, 33493746882302586880000, 1842156078526642278400000

Offset: 0

Wolfdieter Lang

Related to A007559(n+1) ((3*n+1)!!! triple factorials).

Row m=4 of the array A(4; m,n) := ((3*n+m)(!^3))/m(!^3), m >= 0, n >= 0.

G. C. Greubel, Table of n, a(n) for n = 0..378

Cf. A032031, A007559(n+1), A034000(n+1), A034001(n+1). (rows m=0..3).

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(1-3*x)^(7/3))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Aug 15 2018

With[{nn = 30}, CoefficientList[Series[1/(1-3*x)^(7/3), {x, 0, nn}], x]* Range[0,nn]! ] (* G. C. Greubel, Aug 15 2018 *)
With[{c=Times@@Range[4,1,-3]},Table[(Times@@Range[3n+4,1,-3])/c,{n,0,20}]] (* Harvey P. Dale, Feb 06 2023 *)

x='x+O('x^30); Vec(serlaplace(1/(1-3*x)^(7/3))) \\ G. C. Greubel, Aug 15 2018

a(n) = ((3*n+4)(!^3))/4(!^3).
E.g.f.: 1/(1-3*x)^(7/3).
Sum_{n>=0} 1/a(n) = 1 + 3*(3*e)^(1/3)*(Gamma(7/3) - Gamma(7/3, 1/3)). - Amiram Eldar, Dec 23 2022

A051605 a(n) = (3*n+5)!!!/5!!!.

Original entry on oeis.org

1, 8, 88, 1232, 20944, 418880, 9634240, 250490240, 7264216960, 232454942720, 8135922995200, 309165073817600, 12675768026521600, 557733793166950400, 26213488278846668800, 1310674413942333440000, 69465743938943672320000, 3890081660580845649920000

Offset: 0

Wolfdieter Lang

Related to A008544(n+1) ((3*n+2)!!! triple factorials).

Row m=5 of the array A(4; m,n) := ((3*n+m)(!^3))/m(!^3), m >= 0, n >= 0.

G. C. Greubel, Table of n, a(n) for n = 0..378

Cf. A032031, A007559(n+1), A034000(n+1), A034001(n+1), A051604 (rows m=0..4).

m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(1-3*x)^(8/3))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Aug 15 2018

RecurrenceTable[{a[0]==1,a[n]==(3n+5)a[n-1]},a,{n,20}] (* Harvey P. Dale, Oct 19 2013 *)
With[{nn = 30}, CoefficientList[Series[1/(1 - 3*x)^(8/3), {x, 0, nn}], x]*Range[0, nn]!] (* G. C. Greubel, Aug 15 2018 *)

x='x+O('x^30); Vec(serlaplace(1/(1-3*x)^(8/3))) \\ G. C. Greubel, Aug 15 2018

a(n) = ((3*n+5)(!^3))/5(!^3).
E.g.f.: 1/(1-3*x)^(8/3).
a(n) = 3^n*(n+5/3)!/(5/3)!. - Paul Barry, Sep 04 2005
a(n) = (3*n+5)*a(n-1). - R. J. Mathar, Nov 13 2012
Sum_{n>=0} 1/a(n) = 1 + 3*(9*e)^(1/3)*(Gamma(8/3) - Gamma(8/3, 1/3)). - Amiram Eldar, Dec 23 2022

Previous Showing 11-20 of 49 results. Next

cp's OEIS Frontend

A045754 7-fold factorials: a(n) = Product_{k=0..n-1} (7*k+1).

Views

Author

Keywords

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A051188 Sept-factorial numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

PARI

Formula

A004747 Triangle read by rows: the Bell transform of the triple factorial numbers A008544 without column 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Sage

Formula

Extensions

A051141 Triangle read by rows: a(n, m) = S1(n, m)*3^(n-m), where S1 are the signed Stirling numbers of first kind A008275 (n >= 1, 1 <= m <= n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula

Extensions

A114799 Septuple factorial, 7-factorial, n!7, n!!!!!!!, a(n) = n*a(n-7) if n > 1, else 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A144739 7-factorial numbers A114799(7*n+3): Partial products of A017017(k) = 7*k+3, a(0) = 1.

Views

Author

Keywords

Examples

Links

Crossrefs

A144739 7-factorial numbers A114799(7n+3): Partial products of A017017(k) = 7k+3, a(0) = 1.

A010845 a(n) = 3na(n-1) + 1, a(0) = 1.