By using a recent approach proposed by Hackl et al. to evaluate the complexity of the free fermionic Gaussian state, we compute the complexity of the Dirac vacuum state of the Fermi system with a mass quench. First of all, we review the counting method given by Hackl et al., and demonstrate that the result can be adapted to all of the compact transformation group G. Then, we utilize this result to study the time evolution of the complexity of these states. We show that, for the rotational invariant reference state, the total complexity of the incoming vacuum state will saturate the value of the instantaneous vacuum state at the late times, with a typical timescale to achieve the final stable state. Moreover, we find that the complexity growth under the sudden quench is directly proportional to the mass difference. And all of the results show some similar behaviors with the circuit complexity of Thermofield double states.