(PARI) T(n, k) = sum(j=0, n, (-1)^(n-j)*binomial(j, n-k)*binomial(j, k));
vector(30, n, sum(j=0, (n-1)\2, T(n-j-1, j)) ) \\
G. C. Greubel, Dec 07 2019
(Magma) T:= func< n, k | &+[(-1)^(n-j)*Binomial(j, n-k)*Binomial(j, k): j in [0..n]] >;
[(&+[T(n-k, k): k in [0..Floor(n/2)]]): n in [0..30]]; //
G. C. Greubel, Dec 07 2019
(SageMath)
@CachedFunction
def T(n, k): return sum( (-1)^(n-j)*binomial(j, n-k)*binomial(j, k) for j in (0..n))
[sum(T(n-k, k) for k in (0..floor(n/2))) for n in (0..30)] #
G. C. Greubel, Dec 07 2019
(GAP)
T:= function(n, k)
return Sum([0..n], j-> (-1)^(n-j)*Binomial(j, k)*Binomial(j, n-k) );
end;
List([0..30], n-> Sum([0..Int(n/2)], j-> T(n-j, j) )); #
G. C. Greubel, Dec 07 2019