No me des la respuesta que yo realmente quiero pensar que a través de esto por mi cuenta, pero un empujón en la dirección correcta sería impresionante.
Como usted desea; es tu pregunta, después de todo.
Existen los siguientes exponencial en la generación de funciones:
$$\begin{align*}
\frac{-z}{W(-z)(1+W(-z))}&=\sum_{k=0}^\infty\frac{(k-1)^k}{k!}z^k\\
\frac{-W(-z)}{z}&=\sum_{k=0}^\infty\frac{(k+1)^{k-1}}{k!}z^k
\end{align*}$$
donde $W(z)$ es la función de Lambert. (A veces, estos EGFs se expresa en términos de "el árbol de la función" $T(z)=-W(-z)$.) La serie de Maclaurin para $-W(-z)$ tiene este aspecto:
$$-W(-z)=\sum_{k=1}^\infty \frac{k^{k-1}}{k!} z^k$$
que tendrá que demostrar, tal vez por componer dos de la serie.
Ahora, el OP de la identidad puede ser reestructurada como una convolución:
$$\sum_{k=0}^n \frac{(k-1)^k}{k!}\frac{(n-k+1)^{n-k-1}}{(n-k)!}=\frac{n^n}{n!}$$
El cambio a la generación de la función de la vista, la convolución de la izquierda tiene la función de la generación de
$$\left(\frac{-z}{W(-z)(1+W(-z))}\right)\left(\frac{-W(-z)}{z}\right)=\frac1{1+W(-z)}$$
lo que significa que
$$\frac1{1+W(-z)}=1+\sum_{k=1}^\infty\frac{k^k}{k!}z^k$$
Ahora a probar la validez de esta identidad, se puede iniciar mediante la derivación de la fórmula (y el correspondiente de la serie) por la derivada de la función de Lambert. Recordar la relación fundamental $W(z)\exp(W(z))=z$; el uso implícito de diferenciación para derivar una expresión $W^\prime(z)$ que luego pueden ser manipulados para parecerse a $\dfrac1{1+W(-z)}$. De nuevo, los detalles están a la izquierda.