Predicting the usefulness of crosses in terms of expected genetic gain and genetic diversity is of interest to secure performance in the progeny and to maintain long-term genetic gain in plant breeding. A wide range of crossing schemes are possible including large biparental crosses, backcrosses, four-way crosses, and synthetic populations. In silico progeny simulations together with genome-based prediction of quantitative traits can be used to guide mating decisions. However, the large number of multi-parental combinations can hinder the use of simulations in practice. Analytical solutions have been proposed recently to predict the distribution of a quantitative trait in the progeny of biparental crosses using information of recombination frequency and linkage disequilibrium between loci. Here, we extend this approach to obtain the progeny distribution of more complex crosses including two to four parents. Considering agronomic traits and parental genome contribution as jointly multivariate normally distributed traits, the usefulness criterion parental contribution (UCPC) enables to (i) evaluate the expected genetic gain for agronomic traits, and at the same time (ii) evaluate parental genome contributions to the selected fraction of progeny. We validate and illustrate UCPC in the context of multiple allele introgression from a donor into one or several elite recipients in maize (Zea mays L.). Recommendations regarding the interest of two-way, threeway, and backcrosses were derived depending on the donor performance. We believe that the computationally efficient UCPC approach can be useful for mate selection and allocation in many plant and animal breeding contexts. 126 4 parental alleles at QTLs, the 4 -dimensional vector defining the genotype of parent and 127 a p-dimensional vector of zeros. 128 We first concentrate on doubled haploid (DH) lines derived from the 1′ generation (DH-1), 129 and then extend our work to DH lines generated after more selfing generations from the 1′ 130 and to recombinant inbred lines (RILs) at different selfing generations, i.e. partially 131 heterozygous progeny. Absence of selection is assumed while deriving the progeny from 132 generation 1′. In case of DH-1, we denote the ( x 4 )-dimensional genotyping matrix of 133 progeny derived from a four-way cross (Figure 2) in a multi-allelic context as: 134