\[ y'(x)=\frac {-64 a^3 x^3+48 a^2 x^2 y(x)^2+16 a^2 x^2-12 a x y(x)^4-8 a x y(x)^2+y(x)^6+y(x)^4+1}{y(x)} \] ✓ Mathematica : cpu = 0.33751 (sec), leaf count = 128
\[\text {Solve}\left [2 a x=a \text {RootSum}\left [64 \text {$\#$1}^3 a^3-48 \text {$\#$1}^2 a^2 y(x)^2-16 \text {$\#$1}^2 a^2+12 \text {$\#$1} a y(x)^4+8 \text {$\#$1} a y(x)^2+2 a-y(x)^6-y(x)^4-1\& ,\frac {\log (x-\text {$\#$1})}{48 \text {$\#$1}^2 a^2-24 \text {$\#$1} a y(x)^2-8 \text {$\#$1} a+3 y(x)^4+2 y(x)^2}\& \right ]+c_1,y(x)\right ]\]
✗ Maple : cpu = 0. (sec), leaf count = 0 , could not solve
dsolve(diff(y(x),x) = (1+y(x)^4-8*a*x*y(x)^2+16*a^2*x^2+y(x)^6-12*y(x)^4*a*x+48*y(x)^2*a^2*x^2-64*a^3*x^3)/y(x),y(x))