Our approach separates the shape and chemical ordering optimization tasks. We first obtain an ensemble of hundreds of single-metal particles with different stable shapes by carrying out a shape sampling and optimization procedure. Subsequently, we search for the most stable chemical ordering for each shape of the ensemble and for each relative concentration of the two metals, providing the relative energies between all shapes as a function of metal concentration.