We construct and study stationary, asymptotically flat multicenter solutions describing regular black holes with non-Abelian hair (colored magnetic-monopole and dyon fields) in two models of N = 2, d = 4 Super-Einstein-Yang-Mills theories: the quadratic model CP 3 and the cubic model ST [2,6], which can be embedded in 10-dimensional Heterotic Supergravity. These solutions are based on the multicenter dyon recently discovered by one of us, which solves the SU(2) Bogomol'nyi and dyon equations on E 3 . In contrast to the well-known Abelian multicenter solutions, the relative positions of the non-Abelian black-hole centers are unconstrained. We study necessary conditions on the parameters of the solutions that ensure the regularity of the metric. In the case of the CP 3 model we show that it is enough to require the positivity of the "masses" of the individual black holes, the finiteness of each of their entropies and their superadditivity. In the case of the ST [2, 6] model we have not been able to show that analogous conditions are sufficient, but we give an explicit example of a regular solution describing thousands of non-Abelian dyonic black holes in equilibrium at arbitrary relative positions. We also construct non-Abelian solutions that interpolate smoothly between just two aDS 2 ×S 2 vacua with different radii (dumbbell solutions).