The sawtooth crashes with non-monotonic q profiles are numerically investigated by the three-dimensional, toroidal-geometry MHD code (CLT). We find that, when the initial q profiles are non-monotonic, the precursors of the sawtooth crashes are with high n rather than the m/n=1/1 kink mode (where m and n are the poloidal and toroidal mode numbers); and multiple flux tubes form just before sawtooth crashes. The dominant mode transits from higher n to n=1 through mode-mode coupling at the nonlinear stage, and the final pressure crash is still dominated by the 1/1 mode.