\[ \boxed { \left ( {x}^{2}y \left ( x \right ) -1 \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) +x \left ( y \left ( x \right ) \right ) ^{2}-1=0} \]
Mathematica: cpu = 0.011501 (sec), leaf count = 60 \[ \left \{\left \{y(x)\to \frac {1}{x^2}-\frac {\sqrt {c_1 x^2+2 x^3+1}}{x^2}\right \},\left \{y(x)\to \frac {\sqrt {c_1 x^2+2 x^3+1}}{x^2}+\frac {1}{x^2}\right \}\right \} \]
Maple: cpu = 0.032 (sec), leaf count = 50 \[ \left \{ y \left ( x \right ) ={\frac {1}{{x}^{2}} \left ( 1+\sqrt {-2\,{ x}^{2}{\it \_C1}+2\,{x}^{3}+1} \right ) },y \left ( x \right ) =-{\frac { 1}{{x}^{2}} \left ( -1+\sqrt {-2\,{x}^{2}{\it \_C1}+2\,{x}^{3}+1} \right ) } \right \} \]