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