For $x\in[0,1]$, let $f(x):=-x\ln x$ and the two-sample entropy function $H(x)=f(x)+f(1-x)$. Prove $h(x):=\displaystyle\frac{H(x^2)}{H(x)}$ increases.
Here is my proof which is a bit cumbersome. I am seeking a much more elegant approach.
The numerator of the derivative of the sought fraction is
\begin{align} &H(x)^2\frac{dh(x)}{dx} \\ =&\frac d{dx}H(x^2)H(x)-H(x^2)\frac d{dx}H(x) \\ =& 2x\ln\frac{x^2}{1-x^2}\,\big(x\ln x+(1-x)\ln(1-x)\big)-\big(x^2\ln x^2+(1-x^2)\ln(1-x^2)\big)\ln\frac x{1-x} \\ =& 2x^2\ln^2 x+2x(2-x)\ln x\ln(1-x)-(x^2+1)\ln x\ln(1-x^2)+(1-x)^2\ln(1-x)\ln(1-x^2). \tag1\label1 \end{align} All four terms above except the third are positive. I combine the second and the third term together and divide it by $-\ln x$ which is positive, and get \begin{align} g(x):&=-2x(2-x)\ln(1-x)+(x^2+1)\ln(1-x^2) \\ &=(3x-1)(x-1)\ln(1-x)+(x^2+1)\ln(1+x) \tag2\label2 \\ &= \int_0^x \Big(g''(a)-\int_t^a g'''(s)ds\Big)(x-t)dt \end{align} for some $a\in[0,x]$. So we only need to show $g(x)>0, \forall x\in\big(0,\frac13\big]$. $$\frac{d^3g(x)}{dx^3}= \frac{4x(2x^3 +3x^2-2x-7)}{(1-x)^2(1+x)^3}.$$ Let $p(x):=2x^3+3x^2-2x-7$. $p(x)\le p(1)=-4, \forall x\in[0,1]$. This is true since $p(x)$ is convex as $p''(x)=12(x+\frac12)>0$ on that interval and $p(0)=-7<-4=p(1)$. We can take $a=\frac13$ since we can show, with a bit of work, $g''(\frac13)>0$.
(to be continued)
I thought about using the general properties $H(0)=H(1)=1$, $H(x)=H(1-x)$, $H''(x)<0$ and $H'(0)=+\infty$.
However, if you plug in $f(x)=\sin(\pi\sqrt{x})$ and $g(x)=f(x)+f(1-x)$ on a graphing calculator,
$\dfrac{g(x^2)}{g(x)}$ seems to increase from $x=0$ to around $x=0.997$, and then it decreases ever so slightly.
Therefore, you probably need to use the precise expression of $H(x)$, and I don't expect a proof much more elegant than the one you sketched.