Mathematica suggests that integrating the Logarithamic Integral, $\text{li}(x)=\int_0^x \frac{dt}{\log t}$, multiplied by $x^n$, between the limits $0$ and $1$ leads to the following result
$$\int_0^1 x^n\, \text{li}(x) \, dx= - \frac{\log(n + 2)}{(n + 1)}$$
(See this question for how this improper integral was found)
Does anyone recognise this result from the mathematical literature or know of a straight forward proof?
I don't fully understand how Mathematica arrives at this result, particularly in respect of the improper limit of the integral at $1$.