Interested by this post, running cases, I have been able to observe that, for any value of $a$ (positive integer, rational, irrational, complex)
$$\frac{\sum_{i=1}^n i^a}{n^{a+1}}= \frac{H_n^{(-a)}}{(n+1)^{a+1}} =\frac{1}{a+1}+\frac{1}{2 n}+\frac{a}{12
n^2}+O\left(\frac{1}{n^3}\right)$$ I wonder how this could be directly obtained since all my attempts have been totally unsuccessful.
All ideas and suggestions will be very welcome.