General Expression for Three-Loop β-Functions of $$\mathcal{N}{\mathbf{ = 1}}$$ Supersymmetric Theories with Multiple Gauge Couplings Regularized by Higher Covariant Derivatives