We propose an exact Monte Carlo approach for the statistics of discrete quantum systems that does not employ the standard partition of the imaginary time into a mesh and does not contain small parameters. The method operates with discrete objects — kinks, describing virtual transitions at different moments in time. The global statistics of the kinks is reproduced by exact local procedures, the main one being based on the known solution for an asymmetric two-level system.