where the first term in the right hand side is additionally evaluated for MS, XDW, and RMS for each state α. Fortunately, we just need to back transform the density matrices, and the additional computational cost is less than the cost for the major part of Eq. ( ). We do not need this evaluation for XMS, because the Fock matrix is defined with the state-averaged density. We estimate the following order of the computational cost: XMS < MS < RMS ≈ XDW.