\[ y'(x)=\frac {y(x) \left (-a x \log (y(x))+x^2+y(x)\right )}{x (a x-y(x)-y(x) \log (x)-y(x) \log (y(x)))} \] ✓ Mathematica : cpu = 0.700215 (sec), leaf count = 33
\[\text {Solve}\left [a x \log (y(x))-\frac {x^2}{2}-y(x) \log (x)-y(x) \log (y(x))=c_1,y(x)\right ]\] ✓ Maple : cpu = 0.547 (sec), leaf count = 30
\[\{y \left (x \right ) = {\mathrm e}^{\RootOf \left (-2 \textit {\_Z} a x +2 c_{1} a +2 \textit {\_Z} \,{\mathrm e}^{\textit {\_Z}}+x^{2}+2 \,{\mathrm e}^{\textit {\_Z}} \ln \left (x \right )\right )}\}\]