We let function
\begin{equation} \begin{aligned} f(x_1,~x_2,~x_3,~x_4) ~&=~ \sqrt{(x_1+x_2)(x_1+x_3)(x_1+x_4)} \\ &+ \sqrt{(x_2+x_1)(x_2+x_3)(x_2+x_4)} \\ &+ \sqrt{(x_3+x_1)(x_3+x_2)(x_3+x_4)} \\ &+ \sqrt{(x_4+x_1)(x_4+x_2)(x_4+x_3)}, \end{aligned} \end{equation}
where variables $x_1,~x_2,~x_3,~x_4$ are positive and satisfy $x_1+x_2+x_3+x_4 ~=~ 1$. We want to prove that $f(x_1,~x_2,~x_3,~x_4)$ attains its global maximum when $x_1=x_2=x_3=x_4=\frac{1}{4}$.
This looks a difficult problem even if it is inat high-school level. Any clues? Your helpful ideas are greatly appreciated.