I discuss the recent advances in the computation of two-loop scattering amplitudes for five-particle processes with one external massive leg. The latter are crucial ingredients to obtain NNLO QCD predictions for many interesting LHC processes. I present the construction of a basis of transcendental functions which allows us to evaluate all the required Feynman integrals efficiently. This result opened the door to the computation of the first amplitudes of this kind, which however remains challenging due to their algebraic complexity. I discuss how we tame it using finite field arithmetic, and present the computation of the two-loop amplitudes for the production of a W/H boson in association with a bottom-quark pair, and of a W boson in association with a photon and a jet.