According to Landau criterion, a phase transition should be first order when cubic terms of order parameters are allowed in its effective Ginzburg-Landau free energy. Recently, it was shown by renormalization group (RG) analysis that continuous transition can happen at putatively first-order Z3 transitions in 2D Dirac semimetals and such non-Landau phase transitions were dubbed "fermioninduced quantum critical points"(FIQCP) [Li et al., Nature Communications 8, 314 (2017)]. The RG analysis, controlled by the 1/N expansion with N the number of flavors of four-component Dirac fermions, shows that FIQCP occurs for N ≥ Nc. Previous QMC simulations of a microscopic model of SU(N ) fermions on the honeycomb lattice showed that FIQCP occurs at the transition between Dirac semimetals and Kekule-VBS for N ≥ 2. However, precise value of the lower bound Nc has not been established. Especially, the case of N = 1 has not been explored by studying microscopic models so far. Here, by introducing a generalized SU(N ) fermion model with N = 1 (namely spinless fermions on the honeycomb lattice), we perform large-scale sign-problem-free Majorana quantum Monte Carlo simulations and find convincing evidence of FIQCP for N = 1. Consequently, our results suggest that FIQCP can occur in 2D Dirac semimetals for all positive integers N ≥ 1.