Ce travail aborde une variante combinée du problème des p-centres, intégrant simultanément des contraintes de capacité, une stratification multi-services et une gestion de pannes. Face aux limites computationnelles des solveurs exacts sur des instances de grande taille, nous proposons une approche matheuristique fondée sur l'élaboration d'un ordre pour le choix des centres potentiels. La méthode priorise les sites candidats en fonction de leur capacité agrégée, évaluée par leur norme euclidienne au sein d'un hypercube de capacités normalisées. L'algorithme restreint ainsi l'exploration aux k x p sites les plus prometteurs, parmi lesquels la configuration finale est extraite via trois stratégies distinctes : une approche par fenêtre coulissante, une méthode adaptative basée sur une relaxation des contraintes, et une sélection dynamique incrémentale. L'objectif est de garantir l'obtention rapide de solutions réalisables et de qualité.