Estoy consiguiendo:
$$\begin{eqnarray*}\sum_{j=0}^{s}\binom{s}{j}(n-s+j)!(s-j)!&=&\iint_{(0,+\infty)^2}\sum_{j=0}^{s}\binom{s}{j}x^{s-j} y^{n-s+j} e^{-(x+y)}\,dx\,dy\\&=&\iint_{(0,+\infty)^2} x^s y^n e^{-(x+y)}\left(\frac{1}{x}+\frac{1}{y}\right)^s\,dx\,dy\\&=&\iint_{(0,+\infty)^2} y^{n-s} e^{-(x+y)}(x+y)^s\,dx\,dy\\&=&\int_{0}^{+\infty}\int_{0}^{u}y^{n-s} e^{-u} u^s\,dy\,du\\&=&\int_{0}^{+\infty}\int_{0}^{1} v^{n-s} e^{-u} u^{n+1}\,dv\,du\\&=&\color{red}{\frac{(n+1)!}{n-s+1}}\end{eqnarray*}$$
pero no sé si eso cuenta como una combinatoria de la prueba :D
De todos modos, podemos fácilmente hacer que en un adicional de $3^{s-j}$ factor de:
$$\begin{eqnarray*}\sum_{j=0}^{s}\binom{s}{j}3^{s-j}(n-s+j)!(s-j)!&=&\iint_{(0,+\infty)^2}\sum_{j=0}^{s}\binom{s}{j}(3x)^{s-j} y^{n-s+j} e^{-(x+y)}\,dx\,dy\\&=&\iint_{(0,+\infty)^2} y^{n-s} e^{-(x+y)}\left(3x+y\right)^s\,dx\,dy\\&=&\frac{1}{3}\int_{0}^{+\infty}\int_{0}^{1}u^{n+1} v^{n-s}\exp\left(-\frac{u}{3}-uv\right)\,dv\,du\end{eqnarray*}$$
pero eso no simplificar bien como antes, desde que se $\int_{0}^{1}v^{n-s}e^{-uv}\,dv$ depende de la integral exponencial. A partir de la última identidad, con la ayuda de un CAS que tengo:
$$ \sum_{j=0}^{s}\binom{s}{j}3^{s-j}(n-s+j)!(s-j)!\\ = \color{red}{3^{s+1}n!(s-n)!} -\frac{(n+1)!}{s+1}\color{blue}{\,\phantom{}_2 F_1\left(n+2,s+1;s+2;-\frac{1}{3}\right)}. $$
Al menos, el $\color{blue}{\text{non-trivial part}}$ es administrado por un rápido convergente la serie.