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