¿Qué tipo de sustitución se debe usar para obtener las siguientes integrales? $$\begin{align} \int_0^1 \ln \ln \left(\frac{1}{x}\right)\,dx &=\int_0^\infty e^{-x} \ln x\,dx\tag1\\ &=\int_0^\infty \left(\frac{1}{xe^x} - \frac{1}{e^x-1} \right)\,dx\tag2\\ &=-\int_0^1 \left(\frac{1}{1-x} + \frac{1}{\ln x} \right)\,dx\tag3\\ &=\int_0^\infty \left( e^{-x} - \frac{1}{1+x^k} \right)\,\frac{dx}{x},\qquad k>0\tag4\\ \end{align}$$
Esto no es problemas de la tarea y sé que el anterior integrales igual a la negativa de Euler–Mascheroni constante. Tengo estas integrales, mientras que la lectura de esta página de la Wikipedia: El de Euler–Mascheroni constante. De acuerdo a Wikipedia, el de Euler–Mascheroni constante se define como la limitación de la diferencia entre la serie armónica y el logaritmo natural: $$\gamma=\lim_{N\to\infty} \left(\sum_{k=1}^N \frac{1}{k} - \ln N\right)$$ pero no sé por qué esta definición puede ser asociada a la anterior integrales?
Puedo obtener la ecuación de $(1)$ mediante la sustitución de $t=\ln \left(\frac{1}{x}\right)\rightarrow x=e^{-t} \rightarrow dx=-e^{-t}\,dt$, y sé que $$\int_0^\infty e^{-x} \ln x\,dx=\Gamma'(1)=\Gamma(1)\psi(1)=-\gamma$$ pero no puedo obtener el resto. Alguna idea? Cualquier ayuda se agradece. Gracias de antemano.