Nothing is wrong with,
$\frac{m_0}{m_1}=e^{\frac{\sqrt{\frac{\mu_S}{r_1}\left(\sqrt{\frac{2r_2}{r_1+r_2}}-1\right)^{2}+\frac{2\mu_1}{a_1}}-\sqrt{\frac{\mu_1}{a_1}}+\sqrt{\frac{\mu_S}{r_2}\left(\sqrt{\frac{2r_1}{r_1+r_2}}-1\right)^{2}+\frac{2\mu_2}{a_2}}-\sqrt{\frac{\mu_2}{a_2}}}{v_e}}$
assuming that the orbits of the planets are circular and coplanar, and your parking orbits too.
This is basically combining the equation for the $\Delta v$ required for the Hohmann transfer:
$$\Delta v =\sqrt{\frac{\mu_S}{r_1}\left(\sqrt{\frac{2r_2}{r_1+r_2}}-1\right)^{2}+\frac{2\mu_1}{a_1}}-\sqrt{\frac{\mu_1}{a_1}}+\sqrt{\frac{\mu_S}{r_2}\left(\sqrt{\frac{2r_1}{r_1+r_2}}-1\right)^{2}+\frac{2\mu_2}{a_2}}-\sqrt{\frac{\mu_2}{a_2}}$$
And the inverse form of the rocket equation:
$$\frac{m_0}{m_1}=e^{\frac{\Delta v}{v_e}}$$
In some cases, a bi-elliptic transfer might give a better result. For a bi-elliptical transfer of the same type, use:
$$\Delta v=\sqrt{\left(3-2\sqrt{2}\right)\frac{\mu_S}{r_1}+\frac{2\mu_1}{a_1}}-\sqrt{\frac{\mu_1}{a_1}}+\sqrt{\left(3-2\sqrt{2}\right)\frac{\mu_S}{r_2}+\frac{2\mu_2}{a_2}}-\sqrt{\frac{\mu_2}{a_2}}$$
That can be combined with the rocket equation in the same way.
$\frac{m_0}{m_1}=e^{\frac{\sqrt{\left(3-2\sqrt{2}\right)\frac{\mu_S}{r_1}+\frac{2\mu_1}{a_1}}-\sqrt{\frac{\mu_1}{a_1}}+\sqrt{\left(3-2\sqrt{2}\right)\frac{\mu_S}{r_2}+\frac{2\mu_2}{a_2}}-\sqrt{\frac{\mu_2}{a_2}}}{v_e}}$