Por evaluación con WolframAlpha para diferentes valores de $s$ es evidente que:
$$I(s)=\int_0^{\infty} \int_0^{\infty} \frac{\sqrt{xy} ~dxdy}{(x+y)(1+x y)^s}=\frac{\pi}{2(s-1)},~~~~~s>1$$
Yo no estoy muy familiarizado con este tipo de integrales. Podemos introducir nuevas coordenadas como:
$$u=xy,~~~~~v=x+y$$
Pero, a continuación, las expresiones para $x(u,v)$ $y(u,v)$ supeditará en $x>y$ o $x<y$, así que no estoy seguro de cómo transformar la integral correctamente.
También podemos intentar algo como esto para construir una recurrencia:
$$I(s)=\int_0^{\infty} \int_0^{\infty} \frac{xy ~dxdy}{\sqrt{xy}(x+y)(1+x y)^s}= \\ =\int_0^{\infty} \frac{dxdy}{\sqrt{xy}(x+y)(1+x y)^{s-1}}-\int_0^{\infty} \frac{dxdy}{\sqrt{xy}(x+y)(1+x y)^{s}}$$
Pero esto no parece ayudar a cualquiera.