Using path integral quantization in the subspace that forbids double occupancy, we introduce a quantum Monte Carlo algorithm for simulation of fermion models with constraint. The algorithm can be applied to a class of lattice fermion models, including the infinite-U Hubbard model and the t - J model.