Por el Fórmula de Euler-Maclaurin ,
$$\sum_{k=1}^nk^i=\zeta(-i)+\frac{n^{i+1}}{i+1}+\frac{n^i}2+\frac{in^{i-1}}{12}-\frac{i(i-1)(i-2)n^{i-3}}{720}+\mathcal O(n^{i-5})$$
no se trata de una fórmula exacta, sino de una aproximación al problema, donde $\zeta(-i)\approx0.0033002237+0.4181554491i$ es el Función zeta de Riemann . Se dan algunos valores para mostrar cómo se aproxima:
$$\begin{array}{c|c}n&\sum&\zeta\\\hline1&1&0.49913+0.00010i\\2&1.76924+0.63896i&1.38459+0.31946i\\3&2.22407+1.52954i&1.99665+1.08425i\end{array}$$
Vale, no es tan bueno para valores pequeños, pero mejora.
$$\begin{array}{c|c}n&\sum&\zeta\\\hline10&0.41898+7.84548i&0.37600+7.47349i\\20&-8.93237+11.8340i&-8.43768+11.76128i\\30&-18.82708+10.93402i&-18.34384+11.06237i\end{array}$$
Como puede ver, la precisión aumenta a medida que $n\to\infty$ .
Formas alternativas
$$\sum_{k=1}^nk^i=\zeta(-i)-\zeta(-i,n+1)=H_n^{(-i)}$$
donde $\zeta(s,q)$ es el Función zeta de Hurwitz y $H_n^{(p)}$ es el número armónico generalizado .