I don't know which integrators they're using but an algorithm like RKF45 should handle time acceleration with adaptive step sizes no problem. As grandparent said, n-body is very well-known.

Regardless, KSP's workaround of using SOIs is simple and quite effective.