How to calculate the $\beta$-function of $N=1$ supersymmetric theories, regularized by higher covariant derivatives, using vacuum supergraphs

It is demonstrated that in the case of using the higher covariant derivative regularization one can calculate the $\beta$-function of $N=1$ supersymmetric gauge theories defined in terms of the bare couplings on the base of specially modified vacuum supergraphs. The result is automatically obtained in the form of integrals of double total derivatives with respect to loop momenta. The proposed method is illustrated by a number of multiloop calculations.