The QCD Splitting Functions at Three Loops: Methods and Results
  We have computed the complete next-to-next-to-leading order (NNLO) contributions to the splitting functions governing the evolution of unpolarized parton densities in perturbative QCD. Our results agree with all partial results available in the literature. We illustrate the methods used for this calculation with some examples and display selected results to show the size of the NNLO corrections and their effect on the evolution. 
S. Moch, J.A.M. Vermaseren and A. Vogt