(Incluí una introducción más general en la primera versión de esta respuesta, pero luego me di cuenta de que la pregunta se ha actualizado para incluir un principio de derivación de la serie)
Para el enfoque de series infinitas:
Como usted menciona,
$$x^{1/x} = \sum_{n=0}^{\infty} \frac{\log(x)^n}{x^n n!}.$$
Para integrar los términos, podemos aplicar la integración por partes. Tomemos
$$\int \frac{\log(x)^m}{x^n} dx.$$
La razón por la que tenemos dos exponentes separados $m$ y $n$ se aclarará pronto. Deja que $u = \log(x)^m$ , $dv = \frac{1}{x^n} dx$ para que $v = \frac{1}{(1-n) x^{n-1}}$ y $du = \frac{m \log(x)^{m-1}}{x}$ . Entonces,
$$\begin{align}\int \frac{\log(x)^m}{x^n} dx &= uv - \int v\ du \\ &= \frac{\log(x)^m}{(1 - n) x^{n-1}} - \int \frac{m \log(x)^{m-1}}{x (1-n) x^{n-1}} dx \\ &= \frac{\log(x)^m}{(1 - n) x^{n-1}} - \frac{m}{1-n} \int \frac{\log(x)^{m-1}}{x^n} dx. \end{align}$$
Dejemos que $I_{m,n} = \int \frac{\log(x)^m}{x^n} dx$ . Entonces lo anterior puede escribirse como una relación de recurrencia
$$I_{m,n} = \frac{\log(x)^m}{(1 - n)x^{n-1}} - \frac{m}{1-n} I_{m-1,n}$$
con $I_{0,n} = \frac{1}{(1-n) x^{n-1}}$ para $n \ne 1$ y $I_{0,1} = \log(x)$ . Obsérvese que para expresar la integral en esta forma recurrente se requiere que tengamos una $m$ -exponente.
La solución a una recurrencia de la forma
$$a_m = r_m a_{m-1} + s_m$$
con $a_0 = s_0$ viene dada por
$$a_m = \sum_{l=0}^{m} s_{l} \left(\prod_{k=l+1}^{m} r_k\right).$$
(Esto puede obtenerse por reconocimiento de patrones y probarse por sustitución en la recurrencia original).
Dejemos que $r_m = -\frac{m}{1 - n}$ y $s_m = \frac{\log(x)^m}{(1 - n) x^{n-1}}$ para conseguir
$$\begin{align}I_{m,n} &= \sum_{l=0}^{m} \frac{\log(x)^l}{(1 - n) x^{n-1}} \left(\prod_{k=l+1}^{m} -\frac{k}{1 - n}\right) \\ &= \sum_{l=0}^{m} \frac{\log(x)^l}{(1 - n) x^{n-1}} \left((-1)^{m-l} \frac{m!}{l! (1-n)^{m-l}}\right) \\ &= \frac{1}{x^{n-1}} \sum_{l=0}^{m} (-1)^{m-l} \frac{m! \log(x)^l}{l! (1-n)^{m-l+1}}\end{align}$$
para $n \ne 1$ . Sólo nos interesa $I_{n,n}$ Así pues, para $n = 1$ tomamos
$$I_{1,1} = \frac{\log(x)^2}{2}$$
y
$$I_{0,0} = x$$ $$I_{n,n} = \frac{1}{x^{n-1}} \sum_{l=0}^{n} (-1)^{n-l} \frac{n! \log(x)^l}{l! (1-n)^{n-l+1}},\ n > 1.$$
Entonces,
$$\begin{align} \int x^{1/x} dx &= C + \sum_{n=0}^{\infty} \frac{I_{n,n}}{n!} \\ &= C + x + \frac{\log(x)^2}{2} + \sum_{n=2}^{\infty} \left(\frac{1}{x^{n-1}} \sum_{l=0}^{n} (-1)^{n-l} \frac{\log(x)^l}{l! (1-n)^{n-l+1}}\right)\end{align}$$
es la expansión en serie infinita de la integral original. Esta serie debe converger para todo $x > 0$ aunque no funciona tan bien cerca de $0$ .