Tenga en cuenta que esta suma es un caso especial de la suma de armónicos
$$S(x) = \sum_{n\ge 1} \frac{(xn)^k}{2^{xn}}$$
evaluados en $x=1.$
Ahora podemos evaluar esta suma por la inversión de sus Mellin transformar. Vamos a hacer esta evaluación, llegando a una simple forma cerrada que se ve fácilmente que ser un número entero. Recordar el Mellin transformar la identidad de armónicos de las sumas que se
$$\mathfrak{M}\left(\sum_{k\ge 1} \lambda_k g(\mu_k x);s\right) =
\left(\sum_{k\ge 1} \frac{\lambda_k}{\mu_k^s} \right) g^*(s)$$
donde $g^*(s)$ es la transformada de Mellin $g(x).$
En el presente caso tenemos
$$\lambda_k = 1, \quad \mu_k = k \quad \text{y} \quad
g(x) = \frac{x^k}{2^x}.$$
Necesitamos la Mellin transformar $g^*(s)$ $g(x)$ que es
$$\int_0^\infty \frac{x^k}{2^x} x^{m-1} dx =
\int_0^\infty \frac{1}{2^x} x^{(s+k)-1} dx.$$
Observar que
$$\int_0^\infty 2^{-x} x^{m-1} dx =
\frac{1}{(\log 2)^s} \Gamma(s)$$
por una sustitución directa que convierte la integral en una función gamma integral.
Esto implica que
$$g^*(s) = \frac{1}{(\log 2)^{s+k}} \Gamma(s+k).$$
De ello se desprende que la Mellin transformar $Q(s)$ $S(x)$ es
$$Q(s) = \frac{1}{(\log 2)^{s+k}} \Gamma(s+k) \zeta(s).$$
La inversión de esta transformación nos encontramos con que
$$\mathrm{Res}\left(Q(s)/x^s;s=1\right) = \frac{k!}{(\log 2)^{k+1}} \frac{1}{x}.$$
Para el resto de los postes de partida en $s=q$ donde $q=-k$ obtenemos la contribución
$$\sum_{q=k}^\infty \mathrm{Res}\left(Q(s)/x^s;s=-q\ \ derecho) =
\sum_{q=k}^\infty \frac{(-1)^{p-k}}{(p-k)!} (\log 2)^{p-k} x^q \zeta(-q)\\
= - \sum_{q=k}^\infty \frac{(-1)^{p-k}}{(p-k)!} (\log 2)^{p-k} x^q \frac{B_{p+1}}{q+1}
= - x^k \sum_{q=0}^\infty \frac{(-1)^p}{q!} (\log 2)^q x^q \frac{B_{q+k+1}}{q+k+1}.$$
Tenga en cuenta que cuando se $q$ está incluso por encima de la suma de los ceros triviales de la función zeta cancelar los polos de la función gamma. La suma sigue siendo correcto, sin embargo debido a que los números de Bernoulli de la extraña índices son iguales a cero.
Ahora la generación de la función de los números de Bernoulli es
$$\sum_{m\ge 0} B_m \frac{t^m}{m!} = \frac{t}{e^t-1},$$
lo que implica que
$$\sum_{m\ge 1} \frac{B_m}{m} \frac{t^m}{m!} = \log(\exp(t)-1) -t -\log t$$
lo que implica a su vez al $k\ge 1$
$$\sum_{m\ge 1} \frac{B_{m+k+1}}{m+k+1} \frac{t^m}{m.}
= \left(\frac{d}{dt}\right)^{k+1} (\log(\exp(t)-1) -t -\log t)
\\ = \left(\frac{d}{dt}\right)^k \left(-1 - \frac{1}{t} + \frac{e^t}{e^t-1}\right)
= (-1)^{k+1} \frac{k!}{t^{k+1}}
+ \left(\frac{d}{dt}\right)^k \left(\frac{e^t}{e^t-1}\right).$$
Ahora recordando la contribución de la pole en $s=1$ $Q(s)$ tenemos que
$$S(1) = \sum_{n\ge 1}\frac{n^k}{2^n}\\ =
\frac{k!}{(\log 2)^{k+1}}
-\left((-1)^{k+1} \frac{k!}{(-\log 2)^{k+1}}
+ \left. \left(\frac{d}{dt}\right)^k
\left(\frac{e^t}{e^t-1}\right)\right|_{t=\log 2}\right)\\
= \frac{k!}{(\log 2)^{k+1}} - \frac{k!}{(\log 2)^{k+1}}
- \a la izquierda. \left(\frac{d}{dt}\right)^k
\left(\frac{e^t}{e^t-1}\right)\right|_{t=\log 2}
= - \a la izquierda. \left(\frac{d}{dt}\right)^k
\left(\frac{e^t}{e^t-1}\right)\right|_{t=\log 2}.$$
A ver de que esta última expresión es un número entero, volver a escribir como
$$S(1) = - \a la izquierda. \left(\frac{d}{dt}\right)^k
\left(\frac{1}{1-e^{-t}}\right)\right|_{t=\log 2}.$$
Supongamos que tenemos una expresión de la forma
$$\frac{p(e^{-t})}{(1-e^{-t})^a}$$
con $a$ un entero positivo y $p$ un polinomio con coeficientes enteros.
Entonces podemos ver que
$$\frac{d}{dt} \frac{p(e^{-t})}{(1-e^{-t})^a}
= p(e^{-t}) (-a) \frac{1}{(1-e^{-t})^{+1}} e^{-t}
+ p'(e^{-t}) (-e^{-t}) \frac{1}{(1-e^{-t})^a}\\
= \frac{-ap(e^{-t})e^{-t}-p'(e^{-t})e^{-t}(1-e^{-t})}{(1-e^{-t})^{+1}}
= \frac{p(e^{-t})}{(1-e^{-t})^{+1}}.$$
Es evidente que la derivada tiene la misma forma que el término que se diffferentiating, siendo el cociente de un polinomio con coeficientes enteros en $e^{-t}$ y una potencia entera de $(1-e^{-t}).$ Esto es por inducción para todos los derivados.
Pero tenga en cuenta que $$\left.(1-e^{-t})\right|_{t=-\log 2} = -1$$
lo que significa que todos los denominadores que aparecen son iguales a $-1$, dejando un número entero, la contribución de los numeradores, que son polinomios en $e^{-t} = 2$.
Adenda. Hay una muy agradable combinatoria de conexión a esta suma como una función de la $k.$ La secuencia obtenida es
$$2, 6, 26, 150, 1082, 9366, 94586, 1091670, 14174522, 204495126, 3245265146, 56183135190,$$
que tiene este OEIS entrada.
Esto nos dice que el valor obtenido es el número de collares de particiones de $n+1$ etiquetado de perlas. Para demostrar esta conexión, introducir la combinatoria de las especies $\mathcal{P}$ con especificación de $$\mathcal{P} = \mathfrak{C}(\mathfrak{P}_{\ge 1}(\mathcal{Z}))$$
lo que representa, como vemos, los ciclos de la no-vacío etiquetados conjuntos.
Que se traduce en la generación de funciones obtenemos
$$P(z) = \log \frac{1}{1-z} \circ (e^z-1)= \log \frac{1}{2-e^z}.$$
De ahí que el número de collares está dada por (aquí $n\ge 0$)
$$(n+1)! [z^{n+1}] \log \frac{1}{2-e^z} =
(n+1)! [z^{n+1}] \left(\log \frac{1}{2} + \log \frac{1}{1-e^z/2} \right)\\
= (n+1)! [z^{n+1}] \log \frac{1}{1-e^z/2}\\
= (n+1)! [z^{n+1}] \sum_{k\ge 1} \frac{e^{zk}}{k 2^k}
= (n+1)! \sum_{k\ge 1} \frac{k^{n+1}}{(n+1)!\veces k\times 2^k}
= \sum_{k\ge 1} \frac{k^n}{2^k}.$$
Nos hemos cambiado a $k$ $n$ en estos comentarios para ponerlos en línea con la generación de la notación de la función y la notación utilizada en la OEIS entrada. Y, por supuesto, la combinatoria de interpretación se califica como otra prueba de que estamos tratando con números enteros.