Una forma es definir la secuencia de polinomios $\,P_n(x)\,$ tal que $$ P_n(x)\,e^{x+x^2/2} = \left(\frac{d}{dx}\right)^n e^{x+x^2/2}. $$ Los coeficientes de la secuencia polinómica es la Secuencia OEIS A111062 . Consulte la entrada de la OEIS para más detalles. Por supuesto, la función $\,e^{x+x^2/2}\,$ es la función generadora exponencial de la Secuencia OEIS A000085 . La entrada de la OEIS da la fórmula de la suma $$ a_n = \sum_{k=0}^{n/2} \frac{n!}{(n-2k)! 2^k k!}$$ que ya tienes de Wolfram|Alpha. Probablemente el método más fácil es encontrar la serie de potencias para $\,e^x\,$ y $\,e^{x^2/2}\,$ y multiplicarlos, lo que producirá la fórmula.
Más concretamente, considere $$ \sum_{n=0}^\infty a_n\, \frac{x^n}{n!} := e^{x+x^2/2} = e^x\,e^{x^2/2}= \\ \sum_{j=0}^\infty \frac{x^j}{j!} \sum_{k=0}^\infty \frac{x^{2k}}{2^k k!} = \sum_{n=0}^\infty \frac{x^n}{n!} \sum_{k=0}^{n/2} \frac{n!}{(n-2k)! 2^k k!} $$ donde $\,j = n-2k.\,$ Ya que tenemos $\,0\le j,\,0\le k,\,$ y $\,k=(n-j)/2,\,$ entonces debemos tener $\,0\le k\le n/2.$