A description of long-lived photodoped states in Mott insulators is challenging, as it needs to address exponentially separated timescales. We demonstrate how properties of such states can be computed using numerically exact steady state techniques, in particular, the quantum Monte Carlo algorithm, by using a time-local ansatz for the distribution function with separate Fermi functions for the electron and hole quasiparticles. The simulations show that the Mott gap remains robust to large photodoping, and the photodoped state has hole and electron quasiparticles with strongly renormalized properties.