Pairing in the two-dimensional Hubbard model: A Monte Carlo study