A CP model and dedicated global constraints are presented in [DDV_ECML2013][DDV_AIJ][DDV_CP2015]. An implementation in Gecode is available at http://cp4clustering.com.