Aquí es una forma cerrada que involucran un número finito de suma de funciones trigonométricas.
Continuando Quruitay la respuesta, también se puede utilizar la fracción parcial de la descomposición de $$\frac{1}{x^n+1} = -\frac{1}{n}\sum\limits_{k=0}^{n-1} \frac{\zeta_{2n}^{2k+1}}{x-\zeta_{2n}^{2k+1}}$$ to find the integral $$\int_0^1\frac{1}{x^n+1} \mathrm{d}x = -\frac{1}{n}\sum\limits_{k=0}^{n-1}\zeta_{2n}^{2k+1}\ln\left(\zeta_{2n}^{-2k-1}-1\right)$$
Ahora $\zeta_{2n}^{-2k-1} - 1 = 2\left|\sin\left(\frac{2k+1}{2n}\pi\right)\right|e^{-\frac{n+2k+1}{2n}\pi i}$ así obtenemos la suma: $$-\frac{1}{n}\sum\limits_{k=0}^{n-1}\zeta_{2n}^{2k+1}\left[\ln 2 + \ln\left|\sin\left(\frac{2k+1}{2n}\pi\right)\right|-\frac{n+2k+1}{2n}\pi i\right]$$
Asumimos $n \gt 1$. De Vieta las fórmulas aplicadas a $x^n + 1$ sabemos que $\sum\limits_{k=0}^{n-1}\zeta_{2n}^{2k+1} = 0$, por lo que podemos ignorar la constante de $\ln 2$ plazo, así como la constante $-\frac{n+1}{2n}\pi i$ plazo.
Además, podemos mostrar a $\sum\limits_{k=0}^{n-1}k\zeta_{2n}^{2k+1} = \frac{n}{2i\sin\left(\frac{\pi}{n}\right)}$ utilizando una fórmula general para $\sum\limits_{k=1}^{n-1} kX^k$. Estos resultados nos permiten simplificar la expresión como:
$$\frac{\pi}{2n\sin\left(\frac{\pi}{n}\right)}-\frac{1}{n}\sum\limits_{k=0}^{n-1}\zeta_{2n}^{2k+1}\ln\left|\sin\left(\frac{2k+1}{2n}\pi\right)\right|$$
El uso de $\zeta_{2n}^{2k+1} = \cos\left(\frac{2k+1}{n}\pi\right)+i\sin\left(\frac{2k+1}{n}\pi\right)$ y multiplicarse, además de la obtención de:
$$\frac{\pi}{2n\sin\left(\frac{\pi}{n}\right)}-\frac{1}{n}\sum\limits_{k=0}^{n-1}\left[\cos\left(\frac{2k+1}{n}\pi\right)+i\sin\left(\frac{2k+1}{n}\pi\right)\right]\ln\left|\sin\left(\frac{2k+1}{2n}\pi\right)\right|$$
La parte imaginaria se debe cancelar ya que esperamos que el resultado sea real. De hecho, podemos aprovechar la simetría de la suma con respecto a $k \to n - 1 - k$ (cada término en que la parte imaginaria se invierte su signo, mientras que cada término en la parte real sigue siendo el mismo), y por lo tanto podemos mantener la parte real sólo, y mantener la suma de hasta $\big\lfloor\frac{n-1}{2}\big\rfloor$ , mientras que doblando cada término. Tenga en cuenta que si $n$ es impar, tenemos uno de los impares el término correspondiente a $k = \frac{n-1}{2}$, pero resulta que este término es igual a $0$ , de todos modos, con lo que obtenemos:
$$\int_0^1\frac{1}{x^n+1} \mathrm{d}x = \frac{\pi}{2n\sin\left(\frac{\pi}{n}\right)}-\frac{1}{n}\sum\limits_{k=0}^{n-1}\cos\left(\frac{2k+1}{n}\pi\right)\ln\left|\sin\left(\frac{2k+1}{2n}\pi\right)\right|$$
$$\int_0^1\frac{1}{x^n+1} \mathrm{d}x = \frac{\pi}{2n\sin\left(\frac{\pi}{n}\right)}-\frac{2}{n}\sum\limits_{k=0}^{\big\lfloor\frac{n-1}{2}\big\rfloor}\cos\left(\frac{2k+1}{n}\pi\right)\ln\left(\sin\left(\frac{2k+1}{2n}\pi\right)\right)$$
Es interesante notar que la suma es en realidad una suma de Riemann para $\frac{1}{\pi}\int\limits_0^\pi \cos x\ln\left(\sin \frac{x}{2}\right) \mathrm{d}x = -\frac{1}{2}$, por lo que esperamos que el límite de la suma converge a $-\frac{1}{2}$. Pero $\lim\limits_{n\to\infty} \frac{\pi}{2n\sin\left(\frac{\pi}{n}\right)} = \frac{1}{2}$ también, por lo que la resta da $\lim\limits_{n\to\infty}\int\limits_0^1 \frac{1}{x^n+1} \mathrm{d}x = \frac{1}{2} - \left(-\frac{1}{2}\right) = 1$, lo cual no es sorprendente.
Cuando $n$ es incluso podemos aprovechar más la simetría $k \to \frac{n}{2} - 1 - k$ que da $\cos\left(\frac{2\left(\frac{n}{2} - 1 - k\right)+1}{n}\pi\right)\ln\left(\sin\left(\frac{2\left(\frac{n}{2} - 1 - k\right)+1}{2n}\pi\right)\right) = -\cos\left(\frac{2k+1}{n}\pi\right)\ln\left(\cos\left(\frac{2k+1}{2n}\pi\right)\right)$.
Podemos emparejar todos los términos con la posible excepción de la correspondiente a $k = \frac{n-2}{4}$, pero resulta que este término podría ser $0$ de todos modos. Así, obtenemos:
$$\int_0^1\frac{1}{x^n+1} \mathrm{d}x = \frac{\pi}{2n\sin\left(\frac{\pi}{n}\right)}-\frac{2}{n}\sum\limits_{k=0}^{\big\lfloor\frac{n-2}{4}\big\rfloor}\cos\left(\frac{2k+1}{n}\pi\right)\ln\left(\tan\left(\frac{2k+1}{2n}\pi\right)\right) \text{ when }n\text{ is even.}$$
El uso de Quruitay el resultado nos encontramos con que $$\int_0^1\ln\left(1+x^n\right)\mathrm{d}x = \ln 2 - n + \frac{\pi}{2\sin\left(\frac{\pi}{n}\right)} - \sum\limits_{k=0}^{n-1}\cos\left(\frac{2k+1}{n}\pi\right)\ln\left|\sin\left(\frac{2k+1}{2n}\pi\right)\right|$$
$$\int_0^1\ln\left(1+x^n\right)\mathrm{d}x = \ln 2 - n + \frac{\pi}{2\sin\left(\frac{\pi}{n}\right)} - 2\sum\limits_{k=0}^{\big\lfloor\frac{n-1}{2}\big\rfloor}\cos\left(\frac{2k+1}{n}\pi\right)\ln\left(\sin\left(\frac{2k+1}{2n}\pi\right)\right)$$
Y si $n$ incluso:
$$\int_0^1\ln\left(1+x^n\right)\mathrm{d}x = \ln 2 - n + \frac{\pi}{2\sin\left(\frac{\pi}{n}\right)} - 2\sum\limits_{k=0}^{\big\lfloor\frac{n-2}{4}\big\rfloor}\cos\left(\frac{2k+1}{n}\pi\right)\ln\left(\tan\left(\frac{2k+1}{2n}\pi\right)\right)$$
Los resultados son válidos para $n \gt 1$.
He comprobado los valores numéricamente de $n=2$ a $n=20$ en Wolfram Alpha y coinciden con los valores numéricos exactos proporcionados en Claude Leibovici la respuesta, así que estoy bastante seguro de que en el resultado.