OFFSET
0,2
COMMENTS
a(n) is the first entry in the n-th row of a triangular array built from the factorial sequence, where T(0, m) = m! and each subsequent row is obtained by T(n, m) = T(n - 1, m + 3) - T(n - 1, m + 2).
The rows are generated by iterating the operator E^2*Delta, where E is the forward shift operator, E(f(m)) = f(m + 1), and Delta = E - 1 is the forward difference operator. Equivalently, E^2*Delta = E^2*(E - 1) = E^3 - E^2.
FORMULA
Let f(m) = m!. Then:
a(n) = ((E^2*Delta)^n f)(0).
a(n) = Sum_{j=0..n} (-1)^j*binomial(n, j)*(3*n - j)!.
a(n) = Integral_{t=0..oo} t^(2*n)*(t - 1)^n*exp(-t) dt.
More generally, if b_{p, q}(n) = ((E^q*Delta^p)^n f)(0), where f(m) = m!, then b_{p, q}(n) = Integral_{t=0..oo} t^(q*n)*(t - 1)^(p*n)*exp(-t) dt = Sum_{j=0..p*n} (-1)^j*binomial(p*n, j)*((p + q)*n - j)!.
a(n) ~ sqrt(2*Pi) * 3^(3*n + 1/2) * n^(3*n + 1/2) / exp(3*n + 1/3). - Vaclav Kotesovec, May 13 2026
EXAMPLE
First column of
row 0: 1 1 2 6 24 120 720 5040
row 1: 4 18 96 600 4320 35280 322560
row 2: 504 3720 30960 287280 2943360 33022080
row 3: 256320 2656080 30078720 369774720 4906137600
row 4: 339696000 4536362880 64988179200 994447238400
MATHEMATICA
Table[Sum[(-1)^j * Binomial[n, j] * (3*n - j)!, {j, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, May 13 2026 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Dalton Heilig, May 05 2026
STATUS
approved
