The Trotter factorization of the Liouville propagator is used to generate new reversible molecular dynamics integrators. This strategy is applied to derive reversible reference system propagator algorithms (RESPA) that greatly accelerate simulations of systems with a separation of time scales or with long range forces. The new algorithms have all of the advantages of previous RESPA integrators but are reversible, and more stable than those methods. These methods are applied to a set of paradigmatic systems and are shown to be superior to earlier methods. It is shown how the new RESPA methods are related to predictor-corrector integrators. Finally, we show how these methods can be used to accelerate the integration of the equations of motion of systems with Nosé thermostats.