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