Evaluate the integral $\int_0^\infty \frac{dx}{\sqrt{(x^3+a^3)(x^3+b^3)}}$

It was already shown that $$ I_1(p)=\int_0^\infty \frac{dx}{\sqrt{(x^3+1)(x^3+p)}}=\frac{2 \pi}{3 \sqrt{3}} {_2F_1} \left(\frac{1}{2},\frac{2}{3};1;1-p \right). $$ By transformation 2.11(5) from Erdelyi, Higher transcendental functions (put $z=\frac{1-\sqrt{p}}{1+\sqrt{p}}$) $$ {_2F_1} \left(\frac{1}{2},\frac{2}{3};1;1-p \right)=\left(\frac{2}{1+\sqrt{p}}\right)^{4/3}{_2F_1} \left(\frac{2}{3},\frac{2}{3};1;\left(\frac{1-\sqrt{p}}{1+\sqrt{p}}\right)^{2} \right). $$ By Pfaff's transformation $$ {_2F_1} \left(\frac{2}{3},\frac{2}{3};1;\left(\frac{1-\sqrt{p}}{1+\sqrt{p}}\right)^{2} \right)=\left(\frac{(1+\sqrt{p})^2}{4\sqrt{p}}\right)^{2/3}{_2F_1} \left(\frac{1}{3},\frac{2}{3};1;\frac{(1-\sqrt{p})^2}{-4\sqrt{p}} \right). $$ As a result $$ I_1(p)=\frac{2 \pi}{3 \sqrt{3}p^{1/3}}{_2F_1} \left(\frac{1}{3},\frac{2}{3};1;\frac{(1-\sqrt{p})^2}{-4\sqrt{p}} \right). $$ Now we will use a generalization of AGM found by Borwein and Borwein, A Cubic Counterpart of Jacobi's Identity and the AGM, Transactions of the American Mathematical Society, Vol. 323, No. 2, (1991), pp.691-701 (after correcting for some typos): $$ a_{n+1}=\frac{a_n+2b_n}{3} ,\quad b_{n+1}=\sqrt[3]{b_n\frac{a_n^2+a_nb_n+b_n^2}{3}},\quad a_0=1,\quad b_0=s, $$ $$ \quad AG_3(1,s)=\lim_{n\to\infty} a_n=\frac{1}{{_2F_1} \left(\frac{1}{3},\frac{2}{3};1;1-s^3 \right)}. $$ Using this we get

\begin{align} I_1(p)=\frac{2 \pi}{3 \sqrt{3}~p^{1/3}\cdot AG_3\left(1,\left(\frac{1+\sqrt{p}}{2~\sqrt[4]{p}}\right)^{2/3}\right)}. \end{align}

More generally, with $|p-1|<1$, some experimentation shows that, $$\int_0^\infty \frac{dt}{\sqrt{(t^m+1)(t^m+p)}} = \pi\,\frac{\,_2F_1\big(\tfrac12,\tfrac{m-1}{m};1;1-p\big)}{m\sin\big(\tfrac{\pi}{m}\big)}$$ where the question was just the case $m=3$.

Using the advice from @tired in the comments, we can write:

$$I_1(p)=\frac{2 \pi}{3 \sqrt{3}} \sum_{k=0}^\infty \frac{1}{k!^2} \left(\frac{1}{2}\right)_k \left(\frac{2}{3}\right)_k (1-p)^k=$$

$$=\frac{2 \pi}{3 \sqrt{3}} \sum_{k=0}^\infty \frac{1}{(1)_k} \left(\frac{1}{2}\right)_k \left(\frac{2}{3}\right)_k \frac{(1-p)^k}{k!}=\frac{2 \pi}{3 \sqrt{3}} {_2F_1} \left(\frac{1}{2},\frac{2}{3};1;1-p \right)$$

So this is just the usual Gauss hypergeometric function.

This answers my first question, but I'm hoping to get an answer to my second question as well.

If we speak of this integral as a mean, it's very close to both Arithmetic Geometric Mean and the Logarithmic Mean:

$$M(a,b)=\frac{a}{\sqrt{{_2F_1} \left(\frac{1}{2},\frac{2}{3};1;1-\frac{b^3}{a^3} \right)}}$$

$$a \geq b$$

I remind that the AGM can be written as:

$$\text{agm}(a,b)=\frac{a}{{_2F_1} \left(\frac{1}{2},\frac{1}{2};1;1-\frac{b^2}{a^2} \right)}$$

$$a \geq b$$

And numerically we have:

$$M(a,b) \leq \text{agm}(a,b)$$