You don't want your last step since you cannot integrate mdt+ (m/x)dt= m(1+ 1/x)dt with respect to t.
Instead, since this is a linear equation, there is a formula for an "integrating factor" that Subhotosh Kahn mentioned. I would be inclined to look at the "associated homogenous equation", \(\displaystyle \frac{dx}{dt}= \frac{-\beta}{\rho VC_p}x\). That is a "separable equation". We can "separate" it as \(\displaystyle \frac{dx}{x}= \frac{-\beta}{\rho VC_p} dt\). Integrating both sides, \(\displaystyle ln(x)= -\frac{\beta}{\rho VC_p}t+ C'\) where C' is the "constant of integration. Taking the exponential of both sides, \(\displaystyle x= C''e^{-\frac{\beta}{\rho VC_pt}}\) where \(\displaystyle C''= e^{C'}\).
Now look for a single solution to the entire equation. Since the "non-homogeneous" part is the constant, P, we look for a constant solution. Setting x= A, a constant, \(\displaystyle \frac{dx}{dt}= \frac{dA}{dt}= 0\) so the equation becomes \(\displaystyle 0= -\frac{\beta A}{\rho VC_p}+ P\). Solve that for A. \(\displaystyle \frac{\beta A}{\rho VC_p}= P\) so \(\displaystyle A=\frac{P\rho VC_P}{/beta}\).
Since this is a linear equation, we can add those together:
\(\displaystyle x(t)= C''e^{-\frac{\beta t}{\rho VC_p}}+ \frac{P\rho VC_p}{\beta}\).