\[ a y(x)+x y'(x)^2-y(x) y'(x)=0 \] ✓ Mathematica : cpu = 0.679103 (sec), leaf count = 220
\[\left \{\text {Solve}\left [-\frac {-\frac {4 a^{3/2} \sqrt {4-\frac {y(x)}{a x}} \sin ^{-1}\left (\frac {\sqrt {\frac {y(x)}{x}}}{2 \sqrt {a}}\right )}{\sqrt {\frac {y(x)}{x}-4 a}}+\sqrt {\frac {y(x)}{x}} \sqrt {\frac {y(x)}{x}-4 a}+\frac {y(x)}{x}}{4 a}=c_1+\frac {\log (x)}{2},y(x)\right ],\text {Solve}\left [\frac {\frac {4 a^{3/2} \sqrt {4-\frac {y(x)}{a x}} \sin ^{-1}\left (\frac {\sqrt {\frac {y(x)}{x}}}{2 \sqrt {a}}\right )}{\sqrt {\frac {y(x)}{x}-4 a}}-\sqrt {\frac {y(x)}{x}} \sqrt {\frac {y(x)}{x}-4 a}+\frac {y(x)}{x}}{4 a}=c_1-\frac {\log (x)}{2},y(x)\right ]\right \}\] ✓ Maple : cpu = 0.072 (sec), leaf count = 42
\[ \left \{ y \left ( x \right ) =0,y \left ( x \right ) =-{ax \left ( {\it lambertW} \left ( -{\frac {x{\rm e}}{{\it \_C1}\,a}} \right ) -1 \right ) ^{2} \left ( {\it lambertW} \left ( -{\frac {x{\rm e}}{{\it \_C1}\,a}} \right ) \right ) ^{-1}} \right \} \]