Two and three loops computations of renormalization constants for lattice QCD