Considerar la evaluación de las cantidades de la siguiente forma:
$$S_p\stackrel{def}{=}\sum_{k=1}^n\lfloor\sqrt[p]{k}\rfloor$$
donde $p$ es un número entero positivo $\ge 2$. En el caso especial $p = 2$, esto se reduce a la suma queremos calcular.
Deje $a = \lfloor \sqrt[p]{n} \rfloor$ y tomar una lo suficientemente pequeño $\epsilon > 0$ tal que
$$\lfloor \sqrt[p]{k} + \epsilon \rfloor = \lfloor \sqrt[p]{k} \rfloor\quad\text{ for all } k = 1, 2, \ldots, n$$
La introducción de la $\epsilon$ pieza de hacer la función de $\lfloor \sqrt[p]{x} + \epsilon\rfloor$ continua en los saltos de $\lfloor x \rfloor$. Esto nos permite reescribir $\mathcal{S}_p$ como una de Riemann-Stieljes integral.
$$
\mathcal{S}_p = \sum_{k=1}^n \lfloor \sqrt[p]{k} \rfloor
= \sum_{k=1}^n \lfloor \sqrt[p]{k} + \epsilon \rfloor
= \int_{1^-}^{n^+} \lfloor \sqrt[p]{x} + \epsilon \rfloor d \lfloor x \rfloor
$$
Integrar por partes, obtenemos
$$
\mathcal{S}_p
= \bigg[ \lfloor \sqrt[p]{x} + \epsilon \rfloor \lfloor x \rfloor \bigg]_{x=1-}^{n^+}
- \int_{1^-}^{n^+} \lfloor x \rfloor d \lfloor \sqrt[p]{x} + \epsilon \rfloor
= na
- \int_{1+\epsilon}^{\sqrt[p]{n}+\epsilon} \lfloor (s-\epsilon)^p \rfloor d \lfloor s \rfloor
$$
Convertir la última integral de la espalda en una suma, nos encontramos con
$$
\mathcal{S}_p
= na - \sum_{s=2}^(s^p - 1)
= (n+1)a - \sum_{s=1}^s^p
= (n+1)a - \frac{B_{p+1}(a+1)-B_{p+1}(0)}{p+1}
$$
donde $B_k(x)$ es la de Bernoulli polinomio de orden $k$
Para el caso especial $p = 2$, la suma que queremos es
$$
\bbox[8pt,border: 1px solid blue;]{
\sum_{k=1}^n \lfloor\sqrt{k}\rfloor = \mathcal{S}_2 =
(n+1) a - \frac16 un(a+1)(2a+1),\quad a = \lfloor\sqrt{n}\rfloor
}
$$
Esto es equivalente a la fórmula de Markus Scheuer la respuesta.
Para el caso de $p = 3, 4, 5$. esto lleva a
$$\begin{align}
\sum_{k=1}^n \lfloor\sqrt[3]{k}\rfloor = \mathcal{S}_3 & =
(n+1) a - \frac14\left((a+1)^4 - 2(a+1)^3 + (a+1)^2\right)\\
\sum_{k=1}^n \lfloor\sqrt[4]{k}\rfloor = \mathcal{S}_4 & =
(n+1) a - \frac15\left(
(a+1)^5 - \frac52(a+1)^4 + \frac53 (a+1)^3 - \frac16(a+1)
\right)\\
\sum_{k=1}^n \lfloor\sqrt[5]{k}\rfloor = \mathcal{S}_5 & =
(n+1) a - \frac16\left(
(a+1)^6-3(a+1)^5 + \frac52 (a+1)^4 - \frac12 (a+1)^2
\right)\\
\end{align}
$$
la reproducción de las fórmulas de Claude Leibovici la respuesta.