login
Triangle read by rows, T(n,k) = sum(j=0..2*(n-k), A254882(n-k,j)*k^j /(n-k)!), n>=0, 0<=k<=n.
1

%I #8 Feb 10 2015 11:35:06

%S 1,0,1,0,1,1,0,2,2,1,0,6,9,3,1,0,24,48,24,4,1,0,120,300,200,50,5,1,0,

%T 720,2160,1800,600,90,6,1,0,5040,17640,17640,7350,1470,147,7,1,0,

%U 40320,161280,188160,94080,23520,3136,224,8,1

%N Triangle read by rows, T(n,k) = sum(j=0..2*(n-k), A254882(n-k,j)*k^j /(n-k)!), n>=0, 0<=k<=n.

%e [1]

%e [0, 1]

%e [0, 1, 1]

%e [0, 2, 2, 1]

%e [0, 6, 9, 3, 1]

%e [0, 24, 48, 24, 4, 1]

%e [0, 120, 300, 200, 50, 5, 1]

%e [0, 720, 2160, 1800, 600, 90, 6, 1]

%t Flatten[{1,0,1,Table[{0,Table[Sum[Sum[Abs[StirlingS1[n-k,m+1]] * Abs[StirlingS1[n-k,j-m]],{m,0,j-1}]*k^j/(n-k)!,{j,0,2*(n-k)}],{k,1,n-1}],1},{n,2,10}]}] (* _Vaclav Kotesovec_, Feb 10 2015 *)

%o (Sage)

%o T = lambda n,k: sum(A254882(n-k,j)*k^j/factorial(n-k) for j in (0..2*(n-k)))

%o for n in range(6): [T(n,k) for k in (0..n)]

%Y Cf. A254882.

%K nonn,tabl

%O 0,8

%A _Peter Luschny_, Feb 10 2015