I am currently completing a project which involves various numerical sums, and I am trying to use Mathematica to evaluate them and obtain simple analytic results. For example, one such sum is $$f_{(s,k)}=\sum_{l=0}^s(-1)^l\binom{s}{l}^2(s+k-l)!(s-k+l)!$$ for positive integers $s$ and $k$ satisfying $s\geq k$. When I enter the following code into Mathematica
FullSimplify[Sum[(-1)^l (Binomial[s, l])^2*(s + k - l)!*(s + l - k)!, {l, 0, s}], {Element[s, Integers], Element[k, Integers], s >= k}]
it either outputs some ridiculously complicated difference equation, or throws a "Expression Csc[k\ [Pi]] simplified to ComplexInfinity" error. However, I know (or rather suspect, from trialing different values of $s$ and $k$ and making an ansatz) that the analytic result for the above sum is actually $$f_{(s,k)}= (-1)^{s+k}s!s!~,$$ which is very simple.
This happens every time when I try to compute similar kinds of sums. Is there a way to alter the above code to make Mathematica spit out this simple analytical result? Or perhaps Mathematica is not the best program for doing such sums, and there is another program that is better suited?
