The complex-conjugate pole-residue (CCPR) model has been popularly adopted because CCPR-finite-difference time domain (FDTD) can reduce the memory requirement with the help of complex conjugate property of auxiliary variables. To fully utilize CCPR-FDTD, it is of great necessity to investigate its numerical stability since the FDTD method is conditionally stable. Nonetheless, the numerical stability conditions of CCPR-FDTD have not been studied because its derivation is not straightforward. In this communication, the numerical stability conditions of CCPR-FDTD are systematically derived by combining the von Neumann method with Routh-Hurwitz criterion. It is found that the numerical stability conditions of CCPR-FDTD are the same as those of the modified Lorentz-FDTD with bilinear transform. Moreover, the numerical accuracy of CCPR-FDTD is studied, and numerical examples are employed to validate this work.