En el caso en que $r>0$, puedes pensar de la sig. manera:
La expresión $\frac{1^{r}+2^{r}+\ldots+n^{r}}{n^{r+1}}$ no es otra cosa sino el área de la región (poligonal) que se obtiene al embaldosar la región
$$\{(x,y) \in \mathbb{R}^{2} \colon 0\leq x \leq 1, 0 \leq y \leq f(x)\}$$
mediante los $n$ rectángulos siguientes
$$R_{k} = \left\{(x,y) \in \mathbb{R}^{2} \colon \frac{k-1}{n} \leq x \leq \frac{k}{n}, 0 \leq y \leq f\left(\frac{k}{n}\right)\right\}, \quad k= 1, 2, \ldots, n.$$
Ergo, cuando $n$ tiende a infinito la expresión dada tiende a
$$\int_{0}^{1} x^{r} \, dx = \frac{1}{r+1}.$$