We present a new methodology to generate mock halo or galaxy catalogues, which have accurate clustering properties, nearly indistinguishable from full N -body solutions, in terms of the one-point, two-point, and three-point statistics. In particular, the agreement is remarkable, within 1% up to k = 0.55 h Mpc −1 and down to r = 10 h −1 Mpc, for the power spectrum and two-point correlation function respectively, while the bispectrum agrees in general within 20% for different scales and shapes. Our approach is based on the Zel'dovich approximation, however, effectively including with the simple prescriptions the missing physical ingredients, and stochastic scale-dependent, non-local and nonlinear biasing contributions. The computing time and memory required to produce one mock is similar to that using the log-normal model. With high accuracy and efficiency, the effective Zel'dovich approximation mocks (EZmocks) provide a reliable and practical method to produce massive mock galaxy catalogues for the analysis of large-scale structure measurements.