Random-centroid optimization for food formulation