The high speed flow problems usually involve complex flow phenomena, such as strong shock waves, shock-shock interactions, and shear layers. Prediction of these problems requires robust, efficient and accurate numerical methods. A robust flux splitting method capable of capturing crisp shock profile and exact contact surface is presented. Here, the flux vector of the Euler equation is split into convective and pressure parts according to the Toro's formulation. The accuracy of the numerical method at the contact discontinuity is examined first. Sufficient conditions of shock stability for this new method are obtained through a linear perturbation analysis. Several carefully chosen test problems are numerically investigated, and the numerical results demonstrate the accuracy and robustness of the proposed scheme.