Efficient algorithms for -means clustering frequently converge to suboptimal partitions, and given a partition, it is difficult to detect -means optimality. In this paper, we develop an a posteriori certifier of approximate optimality for -means clustering. The certifier is a sub-linear Monte Carlo algorithm based on Peng and Wei's semidefinite relaxation of -means. In particular, solving the relaxation for small random samples of the dataset produces a high-confidence lower bound on the -means objective, and being sub-linear, our algorithm is faster than -means++ when the number of data points is large. We illustrate the performance of our algorithm with both numerical experiments and a performance guarantee: If the data points are drawn independently from any mixture of two Gaussians over with identity covariance, then with probability , our -time algorithm produces a 3-approximation certificate with 99% confidence.
View on arXiv