Aquí hay un método que utiliza funciones de Bessel. Esto probablemente no sea viable en una competencia, pero si conoces las expansiones y te das cuenta de las simplificaciones debido al factor $\sin(2\pi k)$, puede funcionar.
$$\begin{align*} I(0,\pi)&=\int_0^\pi\cos^2(\cos x)+\sin^2(\sin x)\ dx \\ &=\int_0^\pi\cos^2(\cos x)\ dx+\int_0^\pi\sin^2(\sin x)\ dx \\ &=\int_0^\pi\frac{1+\cos(2\cos x)}{2}\ dx+\int_0^\pi\frac{1-\cos(2\sin x)}{2}\ dx\\ &=\pi+\frac{1}{2}\underbrace{\int_0^\pi\cos(2\cos x)\ dx}_{I_1}-\frac{1}{2}\underbrace{\int_0^\pi\cos(2\sin x)\ dx}_{I_2} \end{align*}$$ Utilizando dos expansiones de Jacobi-Anger, $$\cos(2\cos x)=J_0(2)+2\sum_{k=1}^\infty(-1)^kJ_{2k}(2)\cos(2kx),\quad x\in\mathbb{C} \\ \cos(2\sin x)=J_0(2)+2\sum_{k=1}^\infty J_{2k}(2)\cos(2kx),\quad x\in\mathbb{C}$$ tenemos para $I_1$, $$\begin{align*}I_1&=\int_0^\pi J_0(2)+2\sum_{k=1}^\infty(-1)^kJ_{2k}(2)\cos(2kx)\ dx \\&=\pi J_0(2)+2\sum_{k=1}^\infty(-1)^kJ_{2k}(2)\int_0^\pi\cos(2kx)\ dx\\ &=\pi J_0(2)+2\sum_{k=1}^\infty\frac{(-1)^kJ_{2k}(2)\sin(2\pi k)}{k} \\&=\pi J_0(2) \end{align*}$$ y para $I_2$, $$\begin{align*} I_2&=\int_0^\pi J_0(2)+2\sum_{k=1}^\infty J_{2k}(2)\cos(2kx)\ dx \\ &=\pi J_0(2)+2\sum_{k=1}^\infty J_{2k}(2)\int_0^\pi\cos(2kx)\ dx \\ &=\pi J_0(2)+2\sum_{k=1}^\infty\frac{J_{2k}(2)\sin(2\pi k)}{k} \\ &=\pi J_0(2) \end{align*}$$ por lo tanto, la integral es, $$I(0,\pi)=\pi+\frac{\pi }{2}J_0(2)-\frac{\pi}{2}J_0(2)=\pi.$$
Otra solución, esta vez utilizando la expansión de Maclaurin de $\cos(x)$. Como terminamos encontrando formas cerradas en términos de la función de Bessel de todos modos, fusionaré las respuestas.
Trabajando de manera similar a lo anterior, $$\begin{align*} I(0,\pi)&=\pi+\frac{1}{2}\underbrace{\int_0^\pi\cos(2\cos x)\ dx}_{I_1}-\frac{1}{2}\underbrace{\int_0^\pi\cos(2\sin x)\ dx}_{I_2} \end{align*}$$ expandiendo el coseno por su serie de Maclaurin, $$I_1=\int_0^\pi\sum_{k=0}^\infty\frac{(-1)^k2^{2k}}{(2k)!}\cos^{2k}(x)\ dx=\sum_{k=0}^\infty\frac{(-1)^k2^{2k}}{(2k)!}\int_0^\pi\cos^{2k}(x)\ dx=\pi\sum_{k=0}^\infty\frac{(-1)^k}{(k!)^2}=\pi J_0(2)$$ $$I_2=\int_0^\pi\sum_{k=0}^\infty\frac{(-1)^k2^{2k}}{(2k)!}\sin^{2k}(x)\ dx=\sum_{k=0}^\infty\frac{(-1)^k2^{2k}}{(2k)!}\int_0^\pi\sin^{2k}(x)\ dx=\pi\sum_{k=0}^\infty\frac{(-1)^k}{(k!)^2}=\pi J_0(2)$$ aquí usé las integrales de Wallis, $$\int_0^\pi\cos^{2k}(x)\ dx=\int_0^\pi\sin^{2k}(x)\ dx=\frac{\pi}{2^{2k}}\binom{2k}{k},\quad k\in\mathbb{Z_0}$$ y la expansión en serie para la función de Bessel de primer tipo, $$J_0(z)=\sum_{k=0}^\infty\frac{(-1)^k}{(k!)^2}\left(\frac{z}{2}\right)^{2k},\quad z\in\mathbb{C}$$ por lo tanto, con argumentos similares, $$I(0,\pi)=\pi.$$