OFFSET
0,2
COMMENTS
T(n, m) = ^6P_n^m in the notation of the given reference with T(0, 0) = 1.
The monic row polynomials s(n, x) = Sum_{m=0..n} T(n, m)*x^m which are s(n, x) = Product_{k=0..n-1} x-(6+k), n >= 1 and s(0, x) = 1 satisfy s(n, x+y) = Sum_{k=0..n} binomial(n, k)*s(k, x)*S1(n-k, y), with the Stirling1 polynomials S1(n, x) = Sum_{m=1..n} A008275(n, m)*x^m and S1(0, x) = 1. In the umbral calculus (see the S. Roman reference given in A048854) the s(n, x) polynomials are called Sheffer for (exp(6*t), exp(t)-1).
LINKS
Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
Dragoslav S. Mitrinović and Ružica S. Mitrinović, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Pubi. Elektrotehn. Fak. Ser. Mat. Fiz. 77 (1962); alternative link.
FORMULA
T(n, m) = T(n-1, m-1) - (n+5)*T(n-1, m), n >= m >= 0; T(n, m) = 0, n < m; T(n, -1) = 0, T(0, 0) = 1.
E.g.f. for m-th column of signed triangle: ((log(1+x))^m)/(m!*(1+x)^6).
Triangle (signed) = [ -6, -1, -7, -2, -8, -3, -9, -4, -10, ...] DELTA A000035; triangle (unsigned) = [6, 1, 7, 2, 8, 3, 9, 4, 10, 5, 11, ...] DELTA A000035; where DELTA is Deléham's operator defined in A084938.
If we define f(n, i, a) = Sum_{k=0..n-i} binomial(n,k)*stirling1(n-k,i)*Product_{j=0..k-1}(-a-j), then T(n, i) = f(n, i, 6), for n=1,2,...; i=0..n. - Milan Janjic, Dec 21 2008
EXAMPLE
Triangle begins:
1;
-6, 1;
42, -13, 1;
-336, 146, -21, 1;
...
s(2, x)= 42-13*x+x^2; S1(2, x)= -x+x^2 (Stirling1).
MATHEMATICA
t[n_, i_] = Sum[(-1)^k*Binomial[n, k]*Pochhammer[6, k]*StirlingS1[n - k, i], {k, 0, n - i}]; Flatten[Table[t[n, i], {n, 0, 8}, {i, 0, n}]][[1 ;; 45]] (* Jean-François Alcover, Jun 01 2011, after Milan Janjic *)
PROG
(Haskell)
a051338 n k = a051338_tabl !! n !! k
a051338_row n = a051338_tabl !! n
a051338_tabl = map fst $ iterate (\(row, i) ->
(zipWith (-) ([0] ++ row) $ map (* i) (row ++ [0]), i + 1)) ([1], 6)
-- Reinhard Zumkeller, Mar 11 2014
CROSSREFS
KEYWORD
AUTHOR
EXTENSIONS
Name changed by Thomas Scheuerle, Feb 04 2026
STATUS
approved
