This is strange! Isn't one method enough as long as it gives the correct answer?

Yes, this is a linear equation of the form
dxdy+3y=2x+1. The "integrating factor" is
e∫3dx=e3x. Multiplying both sides by that gives
e3xdxdy+3e3xy=dxde3xy=(2x+1)e3x.
To integrate the right side use "integration by parts", letting
u=2x+1,
dv=e3x. Then
du=2dx and
v=31x.
So
∫(2x+1)e3xdx=31(2x+1)e3x−32∫e3xdx=31(2x+1)e3x−92e3x+C=32xe3x+91e3x+C. So
e3xy=32xe3x+91e3x+C so
y=32x+91+Ce−3x.
Another method is separating into the "associated homogeneous equation" and the entire equation. The "associated homogeneous equation" is
dxdy=−3y so
ydy=−3dx. Integrating that gives
ln(y)=−3x+C. Taking the exponential of both sides
y=e−3x+C=C′e−3x where
C′=eC. Now we just need a single function that satisfies the entire equation. Since the "no-homogeneous" part is
2x+1 and I know that the derivative just lowers the degree of a polynomial, I look for a solution of the form
y=ax2+bx+c. Then
y′=2ax+b so the equation is
2ax+b+3ax2+3bx+3c=2x+1.
3ax2+(3b+2a)x+(b+3c)=2x+1. Since this is to be true for all x the coefficients must be the same: 3a= 0, 3b+ 2a= 2 and b+ 3c= 1. That is, a= 0, b= 2/3, and c= 1/9.
The solution is
y=Ce−3x+(2/3)x+1/9 as before.
I'm not sure I would consider it a "different method" but we can also use the method used for higher order linear differential equations. Since the "associated homogeneous equation" is
dxdy+3y=0, the "characteristic equation" is
r+3=0.
r=−3 so the general solution to the associated homogeneous equation is
y=Ce−3x and then finish as above.