Computing classically exact diffusion constants using short-time trajectories