Esto puede ser reescrito como un polylogarithm:
$$\sum_{n=1}^\infty \frac {n^a}{b^n}=\sum_{n=1}^\infty \frac {\left(\frac 1b\right)^n}{n^{-a}}=\operatorname{Li}_{-a}\left(\frac 1b\right)$$
Para enteros positivos $a$ usted puede utilizar estas fórmulas que calcular las sucesivas derivados de $\frac z{1-z}$. Más exactamente :
$$\operatorname{Li}_{-a}(z)=\left(z\frac{\partial}{\partial z}\right)^a\frac z{1-z}=\sum_{k=0}^a k!\;S(a+1,k+1)\left(\frac z{1-z}\right)^{k+1}$$
con $S(n,k)$ los números de Stirling del segundo tipo (reemplace $z$ $\frac 1b$ en la fórmula final).
Para ver por qué esto funciona comienzan con $a=0$ y la nota $z:=\frac 1b$, a continuación, simplemente queremos que la suma de una serie geométrica :
\begin{align}
\operatorname{Li}_0(z)&=\sum_{n=1}^\infty z^n=\frac z{1-z}\\
\operatorname{Li}_{-1}(z)&=\sum_{n=1}^\infty n\,z^n=\left(z\frac{\partial}{\partial z}\right)\sum_{n=1}^\infty z^n=\frac z{(1-z)^2}\\
\operatorname{Li}_{-2}(z)&=\sum_{n=1}^\infty n^2\,z^n=\left(z\frac{\partial}{\partial z}\right)\sum_{n=1}^\infty n\,z^n=\frac {z(1+z)}{(1-z)^3}\\
\end{align}
(y así sucesivamente...)
Otra buena fórmula es :
$$\operatorname{Li}_{-a}(z)=\frac 1{(1-z)^{a+1}}\sum_{k=0}^{a-1} E(a,k)\, z^{a-k}$$
con $E(a,k)$ los números de Euler.