Let us assume that $f \in C^1([0,1])$. We then have
$$(n+1)\int_0^1 x^n f(x)dx = \int_0^1 f(x) d\left(x^{n+1}\right) = f(1) - \int_0^1x^{n+1}f'(x)dx = f(1) - f'(y)\int_0^1 x^{n+1}dx$$
for some $y \in [0,1]$. This so since $f'$ is continuous over $[0,1]$.
Hence, we obtain that
$$\lim_{n \to \infty} (n+1)\int_0^1 x^n f(x)dx = f(1)$$
From the above, we see that the result is true if $f(x)$ happens to be a polynomial.
Given any continuous function $f(x)$ from Weierstrass approximation theorem, $f(x)$ can be uniformly approximated by polynomials on a bounded interval.
This can be easily proved using Bernstein polynomials, uniform continuity, boundedness of continuous function on compact sets and Chebyshev inequality.
Hence, given any $\epsilon > 0$, there exists a polynomial $p_{\epsilon}(x)$ such that
$$\left \vert f(x)-p_{\epsilon}(x) \right \vert < \epsilon \,\,\,\,\,\, \forall x \in [0,1]$$
Hence, we have that
$$\left \vert (n+1)\int_0^1 x^n\left(f(x)-p_{\epsilon}(x) \right)dx\right \vert \leq \epsilon (n+1)\int_0^1 x^n dx = \epsilon$$
for all $n$. Hence, for any $\epsilon>0$, we obtain that
$$\lim_{n \to \infty} (n+1)f(x)dx \in \left(p_{\epsilon}(1)-\epsilon,p_{\epsilon}(1)+\epsilon\right)$$
Since this is true for all $\epsilon$ and since $\lim_{\epsilon \to 0}p_{\epsilon}(x) = f(x)$, we obtain that
$$\lim_{n \to \infty} (n+1)f(x)dx = f(1)$$