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