|
|
A203412
|
|
Triangle read by rows, a(n,k), n>=k>=1, which represent the s=3, h=1 case of a two-parameter generalization of Stirling numbers arising in conjunction with normal ordering.
|
|
4
|
|
|
1, 1, 1, 4, 3, 1, 28, 19, 6, 1, 280, 180, 55, 10, 1, 3640, 2260, 675, 125, 15, 1, 58240, 35280, 10360, 1925, 245, 21, 1, 1106560, 658000, 190680, 35385, 4620, 434, 28, 1, 24344320, 14266560, 4090240, 756840, 100065, 9828, 714, 36, 1
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
COMMENTS
|
Also the Bell transform of the triple factorial numbers A007559 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 A004747 for the triple factorial numbers A008544 as well as A039683 and A132062 for the case of double factorial numbers. - Peter Luschny, Dec 23 2015
|
|
LINKS
|
|
|
FORMULA
|
(1) Is given by the recurrence relation
a(n+1,k) = a(n,k-1)+(3*n-2*k)*a(n,k) if n>=0 and k>=1, along with the initial values a(n,0) = delta_{n,0} and a(0,k) = delta_{0,k} for all n,k>=0.
(2) Is given explicitly by
a(n,k) = (n!*3^n)/(k!*2^k)*Sum{j=0..k} (-1)^j*C(k,j)*C(n-2*j/3-1,n) for all n>=k>=1.
|
|
EXAMPLE
|
Triangle starts:
[ 1]
[ 1, 1]
[ 4, 3, 1]
[ 28, 19, 6, 1]
[ 280, 180, 55, 10, 1]
[ 3640, 2260, 675, 125, 15, 1]
[58240, 35280, 10360, 1925, 245, 21, 1]
|
|
MAPLE
|
A203412 := (n, k) -> (n!*3^n)/(k!*2^k)*add((-1)^j*binomial(k, j)*binomial(n-2*j/3-1, n), j=0..k): seq(seq(A203412(n, k), k=1..n), n=1..9); # Peter Luschny, Dec 21 2015
|
|
MATHEMATICA
|
Table[(n! 3^n)/(k! 2^k) Sum[ (-1)^j Binomial[k, j] Binomial[n - 2 j/3 - 1, n], {j, 0, k}], {n, 9}, {k, n}] // Flatten (* Michael De Vlieger, Dec 23 2015 *)
|
|
PROG
|
(Sage) # uses[bell_transform from A264428]
triplefactorial = lambda n: prod(3*k + 1 for k in (0..n-1))
trifact = [triplefactorial(k) for k in (0..n)]
return bell_transform(n, trifact)
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|