A novel computation method of hybrid capacity constrained centroidal power diagram. (June 2021)