\[ y'(x)=x^2 \sqrt {x^2+8 y(x)-2 x+1}+\sqrt {x^2+8 y(x)-2 x+1}+x^3 \sqrt {x^2+8 y(x)-2 x+1}-\frac {x}{4}+\frac {1}{4} \] ✓ Mathematica : cpu = 1.04202 (sec), leaf count = 107
\[\text {Solve}\left [\frac {x^4}{4}+\frac {x^3}{3}-\frac {1}{4} \sqrt {x^2+8 y(x)-2 x+1}+\frac {1}{4} \log \left (-\sqrt {x^2+8 y(x)-2 x+1}-x+1\right )-\frac {1}{4} \tanh ^{-1}\left (\frac {2 x-2}{2 \sqrt {x^2+8 y(x)-2 x+1}}\right )-\frac {1}{8} \log (y(x))+x=c_1,y(x)\right ]\] ✓ Maple : cpu = 0.362 (sec), leaf count = 32
\[\left \{x^{4}+\frac {4 x^{3}}{3}+c_{1}+4 x -\sqrt {x^{2}-2 x +8 y \left (x \right )+1} = 0\right \}\]