We present the complete set of leading-color two-loop contributions required to obtain next-to-next-to-leading-order (NNLO) QCD corrections to three-jet production at hadron colliders. We obtain analytic expressions for a generating set of finite remainders, valid in the physical region for three-jet production. The analytic continuation of the known Euclidean-region results is determined from a small set of numerical evaluations of the amplitudes. We obtain analytic expressions that are suitable for phenomenological applications and we present a C++ library for their efficient and stable numerical evaluation.