We present an exact path integral methodology for computing quantum dynamical information. This method combines the concepts of iterative propagation with the features of Monte Carlo sampling. The stepwise evaluation of the path integral circumvents the growth of statistical error with time and the use of importance sampling leads to a favorable scaling of required grid points with the number of particles. Three different Monte Carlo sampling procedures are presented. Time correlation functions for several multi-dimensional model systems are computed and accurate long time dynamics are obtained. In the end, the capabilities and limitations of the method are discussed.