\[ 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.575128 (sec), leaf count = 33
DSolve[Derivative[1][y][x] == (y[x]*(x^2 - a*x*Log[y[x]] + y[x]))/(x*(a*x - y[x] - Log[x]*y[x] - Log[y[x]]*y[x])),y[x],x]
\[\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.584 (sec), leaf count = 30
dsolve(diff(y(x),x) = (y(x)-a*ln(y(x))*x+x^2)/(-y(x)*ln(y(x))-y(x)*ln(x)-y(x)+a*x)*y(x)/x,y(x))
\[y \left (x \right ) = {\mathrm e}^{\RootOf \left (-2 \textit {\_Z} a x +2 \ln \left (x \right ) {\mathrm e}^{\textit {\_Z}}+2 \textit {\_Z} \,{\mathrm e}^{\textit {\_Z}}+2 c_{1} a +x^{2}\right )}\]