Computationally efficient models for the dominant and subdominant harmonic modes of precessing binary black holes