%I #13 Dec 10 2023 09:23:38
%S 1,-1,1,2,-2,1,0,12,-3,1,-56,-120,28,-4,1,0,1680,-450,50,-5,1,15840,
%T -30240,10416,-1080,78,-6,1,0,665280,-317520,33712,-2100,112,-7,1,
%U -17297280,-17297280,12070080,-1391040,81648,-3600,152,-8,1
%N Triangle read by rows. T(n, k) = Sum_{j=0..n-k} binomial(-n, j) * A268438(n - k, j).
%e Triangle T(n, k) starts:
%e [0] 1;
%e [1] -1, 1;
%e [2] 2, -2, 1;
%e [3] 0, 12, -3, 1;
%e [4] -56, -120, 28, -4, 1;
%e [5] 0, 1680, -450, 50, -5, 1;
%e [6] 15840, -30240, 10416, -1080, 78, -6, 1;
%e [7] 0, 665280, -317520, 33712, -2100, 112, -7, 1;
%e [8] -17297280, -17297280, 12070080, -1391040, 81648, -3600, 152, -8, 1;
%p A357340 := proc(n, k) local u; u := n - k; (2*u)!*add(binomial(-n, j) * j! *
%p add((-1)^(j+m)*binomial(u+j, u+m)*abs(Stirling1(u+m, m)), m=0..j)/(u +j)!, j=0..u) end: seq(print(seq(A357340(n, k), k=0..n)), n=0..8);
%o (SageMath) # using function A268438
%o def A357340(n, k):
%o return sum(binomial(-n, i) * A268438(n - k, i) for i in range(n - k + 1))
%o for n in range(10): print([A357340(n, k) for k in range(n + 1)])
%Y Cf. A357341 (alternating row sums), A264437, A268438, A357339.
%K sign,tabl
%O 0,4
%A _Peter Luschny_, Sep 25 2022
|