Towards exact results for spectral functions of quantum impurity models in the long-time limit of the multiple-quench time-dependent numerical renormalization group approach