Vamos a escribir su suma como
$$ \bbox[lightyellow] {
\sum\limits_{\left( {0\, \le } \right)\,k\,\left( { \le \,n} \right)} {\left( { - 1} \right)^{\,k} \left( \matriz{
n \cr
n - k \cr} \right)\left( {n + 1 - k} \right)^{\n} } = \sum\limits_{\left( {0\, \le } \right)\,j\,\left( { \le \,n} \right)} {\left( { - 1} \right)^{\,n - j} \left( \matriz{
n \cr
j \cr} \right)\left( {j + 1} \right)^{\n} } = n!
} \etiqueta{1} $$
Considerar el Avance de Diferencia Finita de una función de $f(x)$, y sus iteraciones, que se define como
$$ \bbox[lightyellow] {
\eqalign{
& \Delta _{\,x} \,f(x) = f(x + 1) - f(x) \cr
& \Delta _{\,x} ^{\,2} \,f(x) = \Delta _{\,x} \,\left( {\Delta _{\,x} \,f(x)} \right) = f(x + 2) - 2f(x + 1) + f(x) \cr
& \quad \vdots \cr
& \Delta _{\,x} ^{\,q} \,f(x)\quad \left| {\;0 \le {\rm integer}\,q} \right.\quad = \sum\limits_{\left( {0\, \le } \right)\,k\,\left( { \le \,n} \right)} {\left( { - 1} \right)^{\,q - k} \left( \matriz{
q \cr
k \cr} \right)f(x + k)} \cr}
} \etiqueta{2} $$
Si tomamos $f(x)$ a ser un polinomio de grado $d$
$$ \bbox[lightyellow] {
p_{\,d} (x) = a_{\,d} \,x^{\,d} + a_{\,d - 1} \,x^{\,d - 1} + \; \cdots \; + a_{\,0} \,x^{\,0} = \sum\limits_{0\, \le \,j\, \le \,d} {a_{\,j} \,x^{\,j} }
} \etiqueta{3.a} $$
y nos convierten en el Aumento de Factorial (Pochammer símbolo) por medio de los Números de Stirling,
$$ \bbox[lightyellow] {
p_{\,d} (x) = b_{\,d} \,x^{\,\underline d } + b_{\,d - 1} \,x^{\,\underline {d - 1} } + \; \cdots \; + b_{\,0} \,x^{\,\subrayado 0 } = \sum\limits_{0\, \le \,j\, \le \,d} {b_{\,j} \,x^{\,\underline {\,j\,} } } \quad \left| {\;b_{\,k} = \sum\limits_{0\, \le \,j\, \le \,d} {a_{\,j} \,\left\{ \matriz{
j \cr
k \cr} \right\}\;} } \right.
} \etiqueta{3.b} $$
entonces, gracias a la facilidad de expresión de la diferencia finita para la caída de los factoriales, no es difícil demostrar que
$$ \bbox[lightyellow] {
\eqalign{
& \Delta ^n p_{\,d} (x) = \sum\limits_k {\left( { - 1} \right)^{n - k} \left( \matriz{
n \cr
k \cr} \right)\;p_{\,d} (x + k)} = \cr
& = d!b_{\,d} \,\left( \matriz{
x \cr
d - n \cr} \right) + \left( {d - 1} \right)!b_{\,d - 1} \,\left( \matriz{
x \cr
d - 1 - n \cr} \right) + \; \cdots \; + 0!b_{\,0} \left( \matriz{
x \cr
0 - n \cr} \right) \cr}
} \etiqueta{4} $$
que, calculado en $x=0$ le da:
$$ \bbox[lightyellow] {
\Delta ^n p_{\,d}(0) = \sum\limits_k {\left( { - 1} \right)^{n - k} \left( \matriz{
n \cr
k \cr} \right)\;p_{\,d}(k)} = \left\{ {\matriz{
0 & {d < n} \cr
{d!b_{\,d} = d!a_{\,d} } & {n = d} \cr
{n!b_ {\n} } & {n \le d} \cr
} } \right.
} \etiqueta{5} $$
Por lo que su caso es sólo un caso particular , con $p_{n}(x)=(1+x)^n$, de los de arriba que tiene para todos los polinomios de