A novel one-step leapfrog alternating-direction-implicit finite-difference time-domain method for the modeling of dispersive media described by the complex-conjugate pole-residue pairs model is proposed. The efficiency of the technique is demonstrated in benchmark numerical examples.