DSolve[Derivative[1][y][x] == 1/4 - x/4 + x*Sqrt[1 - 2*x + x^2 + 8*y[x]],y[x],x]