Advertencia: entrante monstruosidad.
Introducción
Usando el teorema del binomio,
\begin{align*}
\int_{0}^{1} \frac{x^n (1 - x)^n}{1 + x^2}\,\text{d}x &= \int_{0}^{1} \frac{x^n}{1+x^2}\sum_{m = 0}^{n} {n \choose m}(-1)^{m}x^m\,\text{d}x \\
&= \sum_{m = 0}^{n} {n \choose m}(-1)^{m} \int_{0}^{1} \frac{x^{n+m}}{1+x^2} \,\text{d}x
\end{align*}
Supongamos primero que $n + m = 2q$ algunos $q$
$$
\frac{x^{t2}}{1+x^2} = \frac{(-1)^{p}}{1 + x^2} - \sum_{\ell = 0}^{p - 1}(-1)^{\ell + p}x^{2\ell}
$$
Por lo tanto,
$$
\int_{0}^{1} \frac{x^{t2}}{1+x^2} = (-1)^{q}\frac{\pi}{4} - \sum_{\ell = 0}^{p - 1} \frac{(-1)^{\ell + p}}{2\ell + 1}
$$
Ahora supongamos que $n + m = 2p + 1$ algunos $p$
$$
\frac{x^{2p + 1}}{1+x^2} = \frac{(-1)^{p}x}{1 + x^2} - \sum_{\ell = 0}^{p - 1}(-1)^{\ell + p}x^{2\ell + 1}
$$
Por lo tanto,
$$
\int_{0}^{1} \frac{x^{2p + 1}}{1+x^2} \,\text{d}x =\frac{(-1)^{p}}{2}\ln(2) - \sum_{\ell = 0}^{p - 1} \frac{(-1)^{\ell + p}}{2\ell + 2}
$$
Caso: $n$ incluso
Si $n$ es incluso, a continuación, escriba $n = 2n'$ y, a continuación,
\begin{align*}
\int_{0}^{1} \frac{x^{2n'} (1 - x)^{2n'}}{1 + x^2}\,\text{d}x &= \sum_{m = 0}^{2n'} {2n' \choose m}(-1)^{m} \int_{0}^{1} \frac{x^{2n'+m}}{1+x^2} \,\text{d}x \\
&= \sum_{m = 0}^{n'} {2n' \choose 2m} \int_{0}^{1} \frac{x^{2(n'+m)}}{1+x^2} \,\text{d}x \\
&- \sum_{m = 0}^{n' - 1} {2n' \choose 2m + 1} \int_{0}^{1} \frac{x^{2(n'+m) + 1}}{1+x^2} \,\text{d}x
\end{align*}
El uso de los resultados anteriores
\begin{align*}
\int_{0}^{1} \frac{x^{2n'} (1 - x)^{2n'}}{1 + x^2}\,\text{d}x &= \sum_{m = 0}^{n'} {2n' \choose 2m}\left[(-1)^{n' + m}\frac{\pi}{4} - \sum_{\ell = 0}^{n' + m - 1} \frac{(-1)^{\ell + n' + m}}{2\ell + 1} \right] \\
&- \sum_{m = 0}^{n' - 1} {2n' \choose 2m + 1} \left[ \frac{(-1)^{n' + m}}{2}\ln(2) - \sum_{\ell = 0}^{n' + m - 1} \frac{(-1)^{\ell + n' + m}}{2\ell + 2} \right]
\end{align*}
Caso: $n$ impar
Si $n$ es impar, a continuación, escriba $n = 2n' + 1$ y, a continuación,
\begin{align*}
\int_{0}^{1} \frac{x^{2n' + 1} (1 - x)^{2n' + 1}}{1 + x^2}\,\text{d}x &= \sum_{m = 0}^{2n' + 1} {2n' + 1 \choose m}(-1)^{m} \int_{0}^{1} \frac{x^{2n'+m + 1}}{1+x^2} \,\text{d}x \\
&= \sum_{m = 0}^{n'} {2n' + 1\choose 2m} \int_{0}^{1} \frac{x^{2(n'+m) + 1}}{1+x^2} \,\text{d}x \\
&- \sum_{m = 0}^{n'} {2n' + 1 \choose 2m + 1} \int_{0}^{1} \frac{x^{2(n'+m + 1)}}{1+x^2} \,\text{d}x
\end{align*}
De nuevo, usando los resultados anteriores,
\begin{align*}
\int_{0}^{1} \frac{x^{2n' + 1} (1 - x)^{2n' + 1}}{1 + x^2}\,\text{d}x &= \sum_{m = 0}^{n'} {2n' + 1 \choose 2m}\left[\frac{(-1)^{n' + m}}{2}\ln(2) - \sum_{\ell = 0}^{n' + m - 1} \frac{(-1)^{\ell + n' + m}}{2\ell + 2} \right] \\
&- \sum_{m = 0}^{n' - 1} {2n' + 1 \choose 2m + 1} \left[ (-1)^{n' + m + 1}\frac{\pi}{4} - \sum_{\ell = 0}^{n' + m} \frac{(-1)^{\ell + n' + m + 1}}{2\ell + 1} \right]
\end{align*}
Esas sumas son todos finitos, por lo que en principio esa es la solución. Yo no puede ver inmediatamente cómo simplificar.