\[ \left (y(x)^2+2 y(x)+x\right ) y'(x)+y(x)^2 (y(x)+x)^2+y(x) (y(x)+1)=0 \] ✓ Mathematica : cpu = 0.615031 (sec), leaf count = 107
\[\left \{\left \{y(x)\to \frac {-x^2-\sqrt {\left (x^2-c_1 x-1\right ){}^2+4 (x-c_1)}+c_1 x+1}{2 (x-c_1)}\right \},\left \{y(x)\to \frac {-x^2+\sqrt {\left (x^2-c_1 x-1\right ){}^2+4 (x-c_1)}+c_1 x+1}{2 (x-c_1)}\right \}\right \}\] ✓ Maple : cpu = 0.119 (sec), leaf count = 116
\[ \left \{ y \left ( x \right ) ={\frac {1}{-2\,{\it \_C1}+4\,x} \left ( -2\,{x}^{2}+{\it \_C1}\,x+\sqrt {4\,{x}^{4}-4\,{\it \_C1}\,{x}^{3}+ \left ( {{\it \_C1}}^{2}-8 \right ) {x}^{2}+ \left ( 4\,{\it \_C1}+16 \right ) x-8\,{\it \_C1}+4}+2 \right ) },y \left ( x \right ) ={\frac {1}{2\,{\it \_C1}-4\,x} \left ( 2\,{x}^{2}-{\it \_C1}\,x+\sqrt {4\,{x}^{4}-4\,{\it \_C1}\,{x}^{3}+ \left ( {{\it \_C1}}^{2}-8 \right ) {x}^{2}+ \left ( 4\,{\it \_C1}+16 \right ) x-8\,{\it \_C1}+4}-2 \right ) } \right \} \]