Unlike our non-polarizable model, where the XC coefficients are obtained through an unconstrained least squares minimization problem, in this study we constrain them to follow exact alternating expressions to avoid both overfitting and unphysical oscillations when extrapolating to shorter interatomic distances than those existing in our calibration dataset. The exact expressions used for every XC coefficient, that is not zeroth order, are: