Aquí es un lugar elementales de la solución.
Observe que $u_n$ es el área de la región $D$ en la unidad de la plaza de $[0,1]^2$ definido por $x^n + y^n \geq 1$. Ahora vamos a $a_n = 2^{-1/n}$ y nos dividimos $D$ en tres partes,
- $ D_1 = \{(x, y) \in [0,1]^2 : x^n + y^n \geq 1 \text{ and } x \leq a_n \} $
- $ D_2 = \{(x, y) \in [0,1]^2 : x^n + y^n \geq 1 \text{ and } y \leq a_n \} $
- $ D_3 = \{(x, y) \in [0,1]^2 : x, y \geq a_n \} $
$\hspace{8em}$
Entonces es fácil comprobar que $D = D_1 \cup D_2 \cup D_3$ y ellos son los que no se superponen. También, aprovechando la simetría, comprobamos que $D_1$ $D_2$ tienen la misma área. Por lo que se deduce que
\begin{align*}
u_n
&= 2\text{[Area of %#%#%]} + \text{[Area of %#%#%]} \\
&= 2\int_{0}^{a_n} (1 - (1-x^n)^{1/n}) \, dx + (1 - a_n)^2.
\end{align*}
Desde $D_1$, el plazo $D_3$ es bueno. Para la integral plazo, aviso que para $a_n = 1 - \mathcal{O}(\frac{1}{n})$,
\begin{align*}
1 - (1-x^n)^{1/n}
&= 1 - e^{\frac{1}{n}\log(1-x^n)}
\stackrel{\text{(1)}}{\leq} -\frac{1}{n}\log(1-x^n) \\
&= \int_{0}^{x} \frac{t^{n-1}}{1-t^n} \, dt
\stackrel{\text{(2)}}{\leq} \int_{0}^{x} 2t^{n-1} \, dt \\
&= \frac{2}{n}x^n.
\end{align*}
Para $(1-a_n)^2$, se utilizó la desigualdad $x \in [0, a_n]$ que vale para todo real $\text{(1)}$. La segunda desigualdad $e^t \geq 1 + t$ se deduce del hecho de que $t$$\text{(2)}$. Así
$1-t^n \geq \frac{1}{2}$$
Esto demuestra la convergencia de $t \in [0,a_n]$.