This is the third of a series of papers on three-loop computation of renormalization constants for Lattice QCD. Our main points of interest are results for the regularization defined by the Iwasaki gauge action and nf=4 Wilson fermions. Our results for quark bilinears renormalized according to the RI’-MOM scheme can be compared to non-perturbative results. The latter are available for twisted mass QCD: being defined in the chiral limit, the renormalization constants must be the same. We also address more general problems. In particular, we discuss a few methodological issues connected to summing the perturbative series such as the effectiveness of boosted perturbation theory and the disentanglement of irrelevant and finite-volume contributi...