Si
$I(n, m)
=\int_0^1 x^n \operatorname{li}(x^m)\,dx
$,
configuración
$y = x^m$,
$x = y^{1/m}$
así
$dx = \frac1{m}y^{1/m-1} dy
$.
Por lo tanto
$\begin{array}\\
I(n, m)
&=\int_0^1 y^{n/m} \operatorname{li}(y)\,\frac1{m}y^{1/m-1}dy\\
&=\frac1{m}\int_0^1 y^{n/m+1/m-1} \operatorname{li}(y)\,dy\\
&=\frac1{m}\int_0^1 y^{(n+1-m)/m} \operatorname{li}(y)\,dy\\
&= \frac1{m}I(\frac{n+1-m}{m}, 1)
\end{array}
$
Para una transformación diferentes,
conjunto
$ $ y = x^{n+1}
$,
así
$x = y^{1/(n+1)}
$
y
$dx
=\frac1{n+1}y^{1/(n+1)-1} dy
$.
Entonces
$\begin{array}\\
I(n, m)
&=\int_0^1 y^{n/(n+1)}\operatorname{li}(y^{m/(n+1)})\,\frac1{n+1}y^{1/(n+1)-1} dy\\
&=\frac1{n+1}\int_0^1 \operatorname{li}(y^{m/(n+1)}) dy\\
&= \frac1{n+1}I(0, \frac{m}{n+1})
\end{array}
$
De acuerdo a Wolfy,
$\int x^r \operatorname{li}(x) dx
= \frac{\operatorname{li}(x) x^{i+1}-Ei((r+2) log(x))}{i+1}
$
y
$\int \operatorname{li}(x^r) dx
= x \operatorname{li}(x^r)-Ei((r+1) log(x))
$.
Estos pueden ser un buen comienzo.