This is separable. The hint is in the terms outside the square root.
But first, lets talk a bit about where this is defined.
Clearly we must have $x \ge x^2y^2$. Which also means we must have $x\ge 0$. Furthermore, you can check that $x=0$ and $y=0$ are both solutions. So we will be looking for solutions in the $x>0$ region.
Lets write the differential as
$$
y^2\sqrt{x(1-xy^2)}\ dx = ydx + 2xdy.
$$
Hopefully, you can see the the left hand side suggests that we can make the substitution $z = xy^2$ (check it), with $z\in (0,1)$.
What comes out the other side is
$$
\frac{dx}{x} = \frac{dz}{\sqrt{z^3(1-z)}}.
$$
We can solve this to get
$$
\ln x = -2\sqrt{\frac{1}{z} -1} +c.
$$
for some constant $c$.
Working our way back up we should find
$$
y = \pm \frac{2}{\sqrt{x}\sqrt{(c - \ln x)^2 +4)}}.
$$