There is no real need to look for substitutions- it's pretty simple as is. This is an "exact" equation.
[tex]\frac{dy}{dx}= \frac{1- 2y- 4x}{1+ y+ 2x}[/tex]
[tex](1+ y+ 2x)dy= (1- 2y- 4x)dx[/tex]
[tex](1+ y+ 2x)dy+ (4x+ 2y- 1)dx= 0[/tex]
Now we see that [itex](1+ y+ 2x)_x= 2= (4x+ 2y- 1)_y[/itex] so this is an "exact differential"- there exist a function F(x,y) such that [itex]dF= F_xdx+ F_ydy= (4x+2y- 1)dx+ (1+ y+ 2x)dx[/itex].
From [itex]F_x= 4x+ 2y- 1[/itex] so that [itex]F= 2x^2+ 2xy- g(y)[/itex].
(Since F_x is a differentiation with respect to x, treating y as a constant, when we "reverse" that the "constant of integration" may be a function of y.)
From that, [itex]F_y= 2x- g'(y)= 1+ y+ 2x[/itex] so that [itex]-g'(y)= 1+ y[/itex] and [itex]g'(y)= -1- y[/itex]. [itex]g(y)= -y- \frac{1}{2}y^2+ C[/itex].
So [itex]F(x,y)= 2x^2+ 2xy- (-y- \frac{1}{2}y^2+ C)= 2x^2+ 2xy+ y+ \frac{1}{2}y^2+ c[/itex].
Saying that "dF= 0" gives [itex]F(x, y)= 2x^2+ 2xy- (-y- \frac{1}{2}y^2+ C)= 2x^2+ 2xy+ y+ \frac{1}{2}y^2+ c[/itex] equals a constant. Including the "c" in that constant,
[itex]2x^2+ 2xy+ y+ \frac{1}{2}y^2= C[/itex].