\[ 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.220843 (sec), leaf count = 29
\[\text {Solve}\left [2 c_1+x^2+2 y(x) (\log (y(x))+\log (x))=2 a x \log (y(x)),y(x)\right ]\]
✓ Maple : cpu = 0.675 (sec), leaf count = 30
\[ \left \{ y \left ( x \right ) ={{\rm e}^{{\it RootOf} \left ( -2\,{\it \_Z}\,ax+2\,{{\rm e}^{{\it \_Z}}}\ln \left ( x \right ) +2\,{\it \_Z}\,{{\rm e}^{{\it \_Z}}}+2\,{\it \_C1}\,a+{x}^{2} \right ) }} \right \} \]