%I #19 May 26 2019 01:39:18
%S 0,0,1,0,4,0,0,13,10,1,0,40,112,40,0,0,121,836,846,116,1,0,364,5264,
%T 11784,5264,364,0,0,1093,30318,129879,129844,30339,1086,1,0,3280,
%U 165792,1242672,2337472,1242672,165792,3280,0,0,9841,878152,10854028,34706584,34706710,10853944,878188,9832,1
%N Eulerian numbers of type D, the complementary type.
%C Named D~(n, k) (the complementary type D triangle) in the Borowiec link.
%H Anna Borowiec, Wojciech Mlotkowski, <a href="http://arxiv.org/abs/1509.03758">New Eulerian numbers of type D</a>, arXiv:1509.03758 [math.CO], 2015.
%H Katarzyna Kril, Wojciech Mlotkowski, <a href="https://www.combinatorics.org/ojs/index.php/eljc/article/view/v26i1p27">Permutations of Type B with Fixed Number of Descents and Minus Signs</a>, Volume 26(1) of The Electronic Journal of Combinatorics, 2019.
%F T(n, k) = (A060187(n+1, k+1) - (-1)^k*binomial(n, k))/2.
%e Triangle begins:
%e 0;
%e 0, 1;
%e 0, 4, 0;
%e 0, 13, 10, 1;
%e 0, 40, 112, 40, 0;
%e 0, 121, 836, 846, 116, 1;
%e 0, 364, 5264, 11784, 5264, 364, 0;
%e ...
%t T[n_, k_] := (Sum[(-1)^(k - i + 1)*(2*i - 1)^n*Binomial[n + 1, k - i + 1], {i, 1, k + 1}] - (-1)^k*Binomial[n, k])/2; Table[T[n, k], {n, 0, 9}, {k, 0, n}] // Flatten (* _Jean-François Alcover_, Feb 19 2019 *)
%o (PARI) B(n, k) = if( n<k || k<1, 0, sum(i=1, k, (-1)^(k-i) * binomial(n, k-i) * (2*i-1)^(n-1)));
%o T(n, k) = (A060187(n+1, k+1) - (-1)^k*binomial(n, k))/2;
%Y Cf. A060187, A262226.
%K nonn,tabl
%O 0,5
%A _Michel Marcus_, Sep 15 2015