Efficient evolution of unpolarized and polarized parton distributions with QCD-PEGASUS. The Fortran package QCD-PEGASUS is presented. This program provides fast, flexible and accurate solutions of the evolution equations for unpolarized and polarized parton distributions of hadrons in perturbative QCD. The evolution is performed using the symbolic moment-space solutions on a one-fits-all Mellin inversion contour. User options include the order of the evolution including the next-to-next-to-leading order in the unpolarized case, the type of the evolution including an emulation of brute-force solutions, the evolution with a fixed number n_f of flavours or in the variable-n_f scheme, and the evolution with a renormalization scale unequal to the factorization scale. The initial distributions are needed in a form facilitating the computation of the complex Mellin moments.