L EARNING -A UGMENTED k - MEANS C LUSTERING
A BSTRACT
k -means clustering is a well-studied problem due to its wide applicability.
Unfortunately, there exist strong theoretical limits on the performance of any algorithm for the k -means problem on worst-case inputs.
To overcome this barrier, we consider a scenario where “advice” is provided to help perform clustering.
Specifically, we consider the k -means problem augmented with a predictor that, given any point, returns its cluster label in an approximately optimal clustering up to some, possibly adversarial, error.
We present an algorithm whose performance improves along with the accuracy of the predictor, even though na¨ıvely following the accurate predictor can still lead to a high clustering cost.
Thus if the predictor is sufﬁciently accurate, we can retrieve a close to optimal clustering with nearly optimal runtime, breaking known computational barriers for algorithms that do not have access to such advice.
We evaluate our algorithms on real datasets and show signiﬁcant improvements in the quality of clustering.
I NTRODUCTION
Clustering is a fundamental task in data analysis that is typically one of the ﬁrst methods used to understand the structure of large datasets.
The most common formulation of clustering is the kmeans problem where given a set P ⊂ R d of n points, the goal is to ﬁnd a set of centers
C ⊂ R d of k points to minimize the objective cost( P, C ) = (cid:80)
p ∈
P min c ∈ C (cid:107) p − c (cid:107) 22 .
Despite decades of work, there exist strong theoretical limitations about the performance of any algorithm for the k -means problem.
Finding the optimal set C is NP-hard even for the case of k =(Dasgupta, 2008) and even ﬁnding an approximate solution with objective value that is within a factor 1 .
07 of the optimal solution is NP-hard (Cohen-Addad & S., 2019; Lee et al., 2017).
Furthermore, the best-known practical polynomial time algorithms can only provably achieve a large constant factor approximation to the optimal clustering, e.g., the 50 -approximation in Song & Rajasekaran (2010), or use techniques such as linear programming that do not scale, e.g., the 6 .
357approximation in Ahmadian et al.
A natural approach to overcome these computational barriers is to leverage the fact that in many applications, the input is often not arbitrary and contains auxiliary information that can be used to construct a good clustering, e.g., in many applications, the input can be similar to past instances.
Thus, it is reasonable to create a (possibly erroneous) predictor by using auxiliary information or through clusterings of similar datasets, which can inform the proper label of an item in our current dataset.
Indeed, inspired by the developments in machine learning, many recent papers have studied algorithms augmented with predictions (Mitzenmacher & Vassilvitskii, 2020).
Such algorithms utilize a predictor that, when invoked, provides an (imperfect) prediction for future inputs.
The predictions are then used by the algorithm to improve performance (see references in Section 1.3).
Hence, we consider the problem of k -means clustering given additional access to a predictor that outputs advice for which points should be clustered together, by outputting a label for each point.
The goal is to ﬁnd k centers that minimize objective (1) and assign each point to one of these centers.
The question is then whether one can utilize such predictions to boost the accuracy and runtime of clustering of new datasets.
Our results demonstrate the answer in the afﬁrmative.
Formal learning-augmented problem deﬁnition.
Given a set P ⊆ R d of n points, the goal is to ﬁnd a set of k points C (called centers) to minimize objective (1).
In the learning-augmented setting, we assume we have access to a predictor Π that provides information about the label of each point consistent with a (1+ α ) -approximately
optimal clustering C .
We say that a predictor has label error rate λ ≤ α if for each label i ∈ [ k ] := { 1 , . . .
, k } , Π errs on at most a λ ≤ α fraction of all points in cluster i in C , and Π errs on at most a λ ≤ α fraction of all points given label i by Π .
In other words, Π has at least (1 − λ ) precision and recall for each label.
Our predictor model subsumes both random and adversarial errors by the predictor.
For example if the cluster sizes are somewhat well-balanced, then a special case of our model is when Π( p ) outputs the correct label of point p ∈ P with some probability 1 − λ
and otherwise outputs a random label in [ k ] with probability λ .
The example where the predictor outputs an adversarial label instead of a random label with probability λ also falls under our model.
For more detail, see Theorems 2.1 and 3.4.
We also adjust our algorithm to have better performance when the errors are random rather than adversarial in the supplementary material.
1. M OTIVATION FOR O UR W ORK
We ﬁrst motivate studying k -means clustering under the learning-augmented algorithms framework.
Overcoming theoretical barriers.
As stated above, no polynomial time algorithm can achieve better than a constant factor approximation to the optimal clustering.
In addition, the best provable approximation guarantees by polynomial time algorithms have a large constant factor (for example the 50 approximation in Song & Rajasekaran (2010)), or use methods which do not scale (such as the linear programming based algorithm in Ahmadian et al.
(2020) which gives a 6 . 357 -approximation).
Therefore, it is of interest to study whether a natural assumption can overcome these complexity barriers.
In our work, we show that knowing the true labels up to some possibly adversarial noise can give us arbitrarily good clusterings, depending on the noise level, which breaks these computational barriers.
Furthermore, we present an algorithm that runs in nearly linear time, rather than just polynomial time.
Lastly, we introduce tools from the robust statistics literature to study k -means clustering rather than the distance-based sampling procedure that is commonly analyzed (this is the basis of kmeans++ ).
This new toolkit and connection could have further applications in other learning-augmented clustering problems.
Practical considerations.
In practice, good predictors can be learned for datasets with auxiliary information.
For a concrete example, we can take any dataset that has a train/test split and use a clustering on the training dataset to help us cluster the testing portion of the dataset.
Therefore, datasets do not have to be speciﬁcally curated to ﬁt our modelling assumption, which is a requirement in other modelling formulations that leverage extra information such as the SSAC model discussed in
Section 1.3.
A predictor can also be created from the natural class of datasets that vary over time, such as Census data or spectral clustering for temporal graphs (graphs slowly varying over time).
For this class of datasets, a clustering from an earlier time step can function as a predictor for later time steps.
Lastly, we can simply use the labels given by another clustering algorithm (such as kmeans++ ) or heuristic as a predictor.
Therefore, predictors are readily and easily available for a wide class of natural datasets.
Following the predictor alone is insufﬁcient.
Given a predictor that outputs noisy labels, it is conceivable that its output alone can give us a good clustering relative to optimal.
However, this is not the case and na¨ıvely using the label provided by the predictor for each point can result in an arbitrarily bad solution, even when the predictor errs with low probability.
For example, consider a cluster of n 2 points at the origin and a cluster of n 2 points at x = 1 .
Then for k = 2 , choosing centers at the origin and at x = 1 induces a k -means clustering cost of zero.
However, even for a predictor that errs with probability 1 n , some point will be mislabeled with constant probability, which results in a positive k -means clustering cost, and so does not provide a relative error approximation.
Thus, using the provided labels by the predictor can induce an arbitrarily bad clustering, even as the label error rate of the predictor tends to zero.
This subtlety makes the model rich and interesting, and requires us to create non-trivial clustering algorithms.
Predictors with adversarial errors.
Since the predictor is separate from the clustering algorithm, interference with the output of the predictor following the clustering algorithm’s query can be a source of non-random noise.
Thus any scenario in which communication is performed over a noisy channel (for example, if the predictor is hosted at one server and the algorithm is hosted at another server) is susceptible to such errors.
Another source of adversarial failure by the predictor is when the predictor is trained on a dataset that can be generated by an adversary, such as in the context of adversarial machine learning.
Moreover, our algorithms have better guarantees when the predictor does not fail adversarially, e.g., see the supplementary material).
1. O UR R ESULTS
In this paper we study “learning-augmented” methods for efﬁcient k -means clustering.
Our contributions are both theoretical and empirical.
On the theoretical side, we introduce an algorithm that provably solves the k -means problem almost optimally, given access to a predictor that outputs a label for each point p ∈ P according to a (1 + α ) -approximately
optimal clustering, up to some noise.
Speciﬁcally, suppose we have access to a predictor Π with label error rate λ upper bounded by a parameter α .
Then, Algorithm 1 outputs a set of centers (cid:101) C in ˜ O
( knd ) time 1 , such that cost ( P, (cid:101) C ) ≤ (1 + O ( α ))
· cost ( P, C opt ) , where C opt is an optimal set of centers.
We improve the runtime in Section 3 by introducing Algorithm 3, which has the same error guarantees, but uses ˜
O ( nd ) runtime, which is nearly optimal since one needs at least nd time to read the points for dense inputs (Theorem 3.4, and Remark A.14).
To output labels for all points, Algorithm 3 requires n queries to the predictor.
However, if the goal is to just output centers for each cluster, then we only require (cid:101) O
( k/α ) queries.
This is essentially optimal; we show in Theorem 3.5 that any polynomial time algorithm must perform approximately (cid:101) Ω( k/α ) queries to output a 1+ α -approximate solution assuming the Exponential Time Hypothesis, a well known complexity-theoretic assumption (Impagliazzo & Paturi, 2001).
Note that one could ignore the oracle entirely, but then one is limited by the constant factor hardness for polynomial time algorithms, which we bypass with a small number of queries.
Surprisingly, we do not require assumptions that the input is well-separated or approximation-stable (Braverman et al., 2011; Balcan et al., 2013), which are assumed in other works.
Finally in the supplementary material, we also give a learning-augmented algorithm for the related problem of k -median clustering, which has less algebraic structure than that of k -means clustering.
We also consider a deletion predictor, which either outputs a correct label or a failure symbol ⊥ and give a (1 + α ) -approximation algorithm
even when the “deletion rate” is 1 − 1 / poly( k ) .
On the empirical side, we evaluate our algorithms on real and synthetic datasets.
We experimentally show that good predictors can be learned for all of our varied datasets, which can aid in clustering.
We also show our methodology is more robust than other heuristics such as random sampling.
1. R ELATED W ORK
Learning-augmented algorithms.
Our paper adds to the growing body of work on learningaugmented algorithms.
In this framework, additional “advice” from a possibly erroneous predictor is used to improve performance of classical algorithms.
For example, a common predictor is a “heaviness” predictor that outputs how “important” a given input point is.
It has been shown that such predictors can be learned using modern machine learning techniques or other methods on training datasets and can be successfully applied to similar testing datasets.
This methodology has found applications in improving data structures (Kraska et al., 2018; Mitzenmacher, 2018), streaming algorithms (Hsu et al., 2019; Jiang et al., 2020), online algorithms (Lykouris & Vassilvtiskii, 2018;
Purohit
Clustering with additional information.
There have been numerous works that study clustering in a semi-supervised setting where extra information is given.
Basu et al.
(2004) gave an active learning framework of clustering with “must-link”/“cannot-link” constraints, where an algorithm is allowed to interact with a predictor that determines if two points must or cannot belong to the same cluster.
Their objective function is different than that of k -means and they do not give theoretical bounds on the quality of their solution.
Balcan & Blum (2008) and Awasthi et al.
(2017) studied an interactive framework for clustering, where a predictor interactively provides feedback about whether or not to split a current cluster or merge two clusters.
Vikram & Dasgupta (2016) also worked with an interactive oracle but for the Bayesian hierarchical clustering problem.
These works differ from ours in their assumptions since their predictors must answer different questions about partitions of the input points.
In contrast, Howe (2017) used logistic regression to aid k -means clustering but do not give any theoretical guarantees.
The framework closest in spirit to ours is the semi-supervised active clustering framework (SSAC) introduced by Ashtiani et al.
and
further studied by Kim & Ghosh (2017); Mazumdar & Saha (2017)
; Gamlath et al.
; Ailon et al.
; Chien et al.
; Huleihel et al.
The goal of this framework is also to produce a (1 + α ) -approximate clustering while minimizing the number of queries to a predictor that instead answers queries of the form “same-cluster ( u, v ) ”, which returns 1 if points u, v ∈ P are in the same cluster in a particular optimal clustering and 0 otherwise.
Our work differs from the SSAC framework in terms of both runtime guarantees, techniques used, and model assumptions, as detailed below.
We brieﬂy compare to the most relevant works in the SSAC framework, which are Ailon et al.
(2018) and Mazumdar & Saha (2017).
First, the runtime of Ailon et al.
(2018) is O
( ndk 9 /α 4 ) even for a perfectly accurate predictor, while the algorithm of Mazumdar & Saha (2017) uses O ( nk 2 ) queries and runtime
˜ O ( ndk 2 ) .
By comparison, we use signiﬁcantly fewer queries, with near linear runtime ˜ O ( nd ) even for an erroneous predictor.
Moreover, a predictor of Mazumdar & Saha (2017) independently fails each query with probability p so that repeating with pairs containing the same point can determine the correct label of a point whereas our oracle will always repeatedly fail with the same query , so that repeated queries do not help.
The SSAC framework uses the predictor to perform importance sampling to obtain a sufﬁcient number of points from each cluster whereas we use techniques from robust mean estimation, dimensionality reduction, and approximate nearest neighbor data structures.
Moreover, it is unclear how the
SSAC predictor can be implemented in practice to handle adversarial corruptions.
One may consider simulating the SSAC predictor using information from individual points by simply checking if the labels of the two input points are the same.
However, if a particular input is mislabeled, then all of the pairs containing this input can also be reported incorrectly, which violates their independent noise assumption.
Finally, the noisy predictor algorithm in Ailon et al.
(2018) invokes a step of recovering a hidden clique in a stochastic block model, making it prohibitively costly to implement.
Lastly, in the SSAC framework, datasets need to be speciﬁcally created to ﬁt into their model since one requires pairwise information.
In contrast, our predictor requires information about individual points, which can be learned from either a training dataset, from past similar datasets, or from another approximate or heuristic clustering and is able to handle adversarial corruptions.
Thus, we obtain signiﬁcantly faster algorithms while using an arguably more realistic predictor.
Approximation stability.
Another approach to overcome the NP-hardness of approximation for k -means clustering is the assumption that the underlying dataset follows certain distributional properties.
Introduced by Balcan et al.
, the notion of ( c, α ) -approximate stability (Agarwal et al., 2015; Awasthi et al., 2019; Balcan et al., 2020) requires that every c -approximation is α -close to the optimal solution in terms of the fraction of incorrectly clustered points.
In contrast, we allow inputs so that an arbitrarily small fraction of incorrectly clustered points can induce arbitrarily bad approximations, as previously discussed, e.g., in Section 1.1.
L EARNING -A UGMENTED k - MEANS A LGORITHM
Preliminaries.
We use [ n ] to denote the set { 1 , . . . , n } .
Given the set of cluster centers C , we can partition the input points P into k clusters { C 1 , . . .
, C k } according to the closest center to each point.
If a point is grouped in C i in the clustering, we refer to its label as i .
Note that labels can be arbitrarily permuted as long as the labeling across the points of each cluster is consistent.
It is well-known that in k -means clustering, the i -th center is given by the coordinate-wise mean of the Let Y i be the set of points with label i . Run C E ST d coordinates
Y i .
Let C (cid:48) i be the coordinate-wise outputs of C RD E ST . for
C
C
Input: Points x 1 , . .
, x 2 m ∈ R corrup tion λ α
Randomly partition the points into shortest interval
I
1 | Z
∈
Z x
Return points in C i .
Given x ∈ R d and a set C ⊂ R d , we deﬁne d ( x, C ) =
min c ∈ C (cid:107) x − c
(cid:107) 2 .
Note that there may be many approximately optimal clusterings but we consider a ﬁxed one for our analysis.
2. O UR A LGORITHM
Our main result is an algorithm for outputting a clustering that achieves a (1 + 20 α ) approximationto the optimal objective cost when given access to approximations of the correct labeling of the points in P .
We ﬁrst present a suboptimal algorithm in Algorithm 1 for intuition and then optimize the runtime in Algorithm 3, which is provided in Section 3.
The intuition for Algorithm 1 is as follows.
We ﬁrst address the problem of identifying an approximate center for each cluster.
Let C opt 1 , · · · , C optk be an optimal grouping of the points and consider all the points labeled i by our predictor for some ﬁxed 1 ≤ i ≤ k .
Since our predictor can err, a large number of points that are not in C opti may also be labeled i .
This is especially problematic when points that are “signiﬁcantly far” from cluster C opti are given the label i , which may increase the objective function arbitrarily if we simply take the mean of the points labeled i by the predictor.
To ﬁlter out such outliers, we consider a two step view from the robust statistics literature, e.g., Prasad et al.
; these two steps can respectively be interpreted as a “training” phase and a “testing” phase that removes “bad” outliers.
We ﬁrst randomly partition the points that are given label i into two groups, X 1 and X 2 , of equal size.
We then estimate the mean of C opti using a coordinate-wise approach through Algorithm 2 (C RD E ST ), decomposing the total cost as the sum of the costs in each dimension.
For each coordinate, we ﬁnd the smallest interval I that contains a (1 − 4 α ) fraction of the points in X 1 .
We show that for label error rate λ ≤ α , this “training” phase removes any outliers and thus provides a rough estimation to the location of the “true” points that are labeled i .
To remove dependency issues, we then “test” X 2 on I by computing the mean of X 2 ∩ I .
This allows us to get empirical centers that are a sufﬁciently good approximation to the coordinates of the true center for each coordinate.
We then repeat on the other labels.
The key insight is that the error from meanestimation can be directly charged to the approximation error due to the special structure of the k -means problem.
Our main theoretical result considers predictors that err on at most a λ -fraction of all cluster labels.
Note that all omitted proofs appear in the supplementary material.
Theorem 2.1.
Let α ∈
(10 log n/ √ n, 1 / 7) , Π be a predictor with label error rate λ ≤ α , and γ ≥ 1 a sufﬁciently large constant.
If each cluster in the (1 + α ) -approximately
optimal k -means clustering of the predictor has at least γηk/α points, then Algorithm 1 can be used to output a (1 + 20 α ) -approximation to the k -means objective with prob.
1 − 1 /η , using O ( kdn log n ) runtime.
We improve the running time to O ( nd log n + poly( k, log n )) in Theorem 3.4 in Section 3.
Our algorithms can also tolerate similar error rates when failures correspond to random labels, adversarial labels, or a special failure symbol.
Error rate λ vs. accuracy parameter α .
We emphasize that λ is the error rate of the predictor and α is only some loose upper bound on λ .
It is reasonable that some algorithms can provide lossy guarantees on their outputs, which translates to the desired loose upper bound α on the accuracy of the predictor.
Even if is not known, multiple instances of the algorithm can be run in parallel with separate exponentially decreasing “guesses” for the value α .
We can simply return the best clustering among these algorithms, which will provide the same theoretical guarantees as if we set α = 1 .
01 λ , for example.
Thus α does not need to be known in advance and it does not need to be tuned as a hyperparameter.
N EARLY O PTIMAL R UNTIME A LGORITHM
We now describe Algorithm 3, which is an optimized runtime version of Algorithm 1 and whose guarantees we present in Theorem 3.4.
The bottleneck for Algorithm 1 is that after selecting k empirical centers, it must still assign each of the n points to the closest empirical center.
The main intuition for Algorithm 3 is that although reading all points uses O ( nd ) time, we do not need to spend
O ( dk ) time per point to ﬁnd its closest empirical center, if we set up the correct data structures.
In fact, as long as we assign each point to a “relatively good” center, the assigned clustering is still a “good” approximation to the optimal solution.
Thus we proceed in a similar manner as before to sample a number of input points and ﬁnd the optimal k centers for the sampled points.
We use dimensionality reduction and an approximate nearest neighbor (ANN) data structure to efﬁciently assign each point to a “sufﬁciently close” center.
Namely if a point p ∈ P should be assigned to its closest empirical C i then p must be assigned to some empirical center C j such that (cid:107) p − C j (cid:107) 2 ≤ 2 (cid:107) p − C i (cid:107) 2 .
Hence, points that are not assigned to their optimal centers only incur a “small” penalty due to the ANN data structure and so the cost of the clustering does not increase “too much” in expectation.
Formally, we need the following deﬁnitions.
Theorem 3.1 (JL transform) .
Johnson & Lindenstrauss (1984) Let d ( · , · ) be the standard Euclidean norm.
There exists a family of linear maps A : R d → R k and an absolute constant C > 0 such that for any x, y ∈ R d , Pr [ φ ∈ A, d ( φ ( x ) , φ ( y )) ∈
(1 ± α ) d ( x, y )] ≥ 1 − e − Cα 2 k .
Deﬁnition 3.2 (Terminal dimension reduction) .
Given a set of points called terminals C ⊂ R d , we call a map f : R d → R k a terminal dimension reduction with distortion D if for every terminal c ∈ C and point p ∈ R d , we have d ( p, c ) ≤
d ( f ( p ) , f ( c )) ≤ D · d ( p, c ) .
Deﬁnition 3.3 (Approximate nearest neighbor search) .
Given a set P of n points in a metric space ( X, d ) , a ( c, r ) -approximate nearest neighbor search (ANN) data structure takes any query point q ∈
X with non-empty { p ∈ P : 0 < d ( p, q ) ≤ r } and outputs a point in { p ∈ P : 0 < d ( p, q ) ≤ cr } .
To justify the guarantees of Algorithm 3, we need runtime guarantees on creating a suitable dimensionality reduction map and an ANN data structure.
These are from Makarychev et al.
and
Indyk & Motwani (1998);
Har-Peled et al.
; Andoni et al.
(2018) respectively, and are stated in Theorems A. 12 and A. 13 in the supplementary section.
They ensure that each point is mapped to a “good” center.
Thus, we obtain our main result describing the guarantees of Algorithm 3.
Theorem 3.4.
Let α ∈
(10 log n/ √ n, 1 / 7) , Π be a predictor with label error rate λ ≤ α , and γ ≥be a sufﬁciently large constant.
If each cluster in the optimal k -means clustering of the predictor has at least γk log k/α points, then Algorithm 3 outputs a (1 + 20 α ) -approximation to the k -means objective with probability at least 3 / 4 , using O ( nd log n + poly( k, log n )) total time.
Note that if we wish to only output the k centers rather than labeling all of the input points, then the query complexity of Algorithm 3 is (cid:101)
O
( k/α ) (see Step 1 of Algorithm 3) with high probability.
We show in the supplementary material that this is nearly optimal.
Theorem 3.5.
For any δ ∈ (0 , 1] , any algorithm that makes O (cid:16) k 1 − δ α log n (cid:17) queries to the predictor with label error rate α cannot output a (1 + Cα ) -approximation to the optimal k -means clustering cost in time 2 O ( n 1 − δ ) time, assuming the Exponential Time Hypothesis. 2: Let C 1 , . . . C k be the output of Algorithm 1 on S .
3: Let φ 2 be a random JL linear map with distortion 54 , i.e., dimension O (log n ) . 4: Let φ 1 be a terminal dimension reduction with distortion 54 . 5: Let φ := φ 1 ◦ φ 2 be the composition map.
6: Let A be a (2 , r
-ANN data structure on the points φ ( C 1 ) , . .
, φ ( C
7: for x ∈ X do 8:
Let (cid:96) x be the label of x from Π .
← d
Query to ﬁnd
11: if d ( x, C p x ) < 2 d (
C (cid:96)
) then 12:
Assign label p x to x .
13: else 14:
Assign
15: end if 16: end for
E XPERIMENTS
In this section we evaluate Algorithm 1 empirically on real datasets.
We choose to implement
Algorithm 1, as opposed to the runtime optimal Algorithm 3, for simplicity and because the goal of our experiments is to highlight the error guarantees of our methodology, which both algorithms share.
Further, we will see that Algorithm 1 is already very fast compared to alternatives.
Thus, we implement the simpler of the two algorithms.
We primarily ﬁx the number of clusters to be k = 10 and k = 25 throughout our experiments for all datasets.
Note that our predictors can readily generalize to other values of k
but we focus on these two values for clarity.
All of our experiments were done on a CPU with i5 2.7 GHz dual core and 8 GB RAM.
Furthermore, all our experimental results are averaged over 20 independent trials and ± one standard deviation error is shaded when applicable.
We give the full details of our datasets below.
1) Oregon : Dataset of 9 graph snapshots sampled across 3 months from an internet router communication network (Leskovec et al., 2005).
We then use the top two eigenvectors of the normalized Laplacian matrix to give us node embeddings into R 2 for each graph which gives us 9 datasets, one for each graph.
Each dataset has roughly n ∼ 10 4 points.
This is an instance of spectral clustering.
2) PHY : Dataset from KDD cup 2004 (kdd, 2004).
We take 10 4 random samples to form our dataset. 3) CIFAR10 :
Testing portion of CIFAR-10 ( n = 10 4 , dimension 3072 )
(Krizhevsky, 2009).
Baselines.
We compare against the following algorithms.
Additional experimental results on
Lloyd’s heuristic are given in Section E.3 in the supplementary material
kmeans++ : We measure the performance of our algorithm in comparison to the kmeans++ seeding algorithm.
Since kmeans++ is a randomized algorithm, we take the average clustering cost after running kmeans++ seeding on 20 independent trials.
We then standardize this value to have cost 1 . 0 and report all other costs in terms of this normalization.
For example, the cost 2 . 0 means that the clustering cost is twice that of the average kmeans++ clustering cost.
We also use the labels of kmeans++ as the predictor in the input for Algorithm 1 (denoted as “Alg + kmeans++ ”) which serves to highlight the fact that one can use any heuristic or approximate clustering algorithm as a predictor.
Random sampling : For this algorithm, we subsample the predictor labels with probability q ranging from 1% to 50% .
We then construct the k -means centers using the labels of the sampled points and measure the clustering cost using the whole dataset.
We use the best value of q in our range every time to give this baseline as much power as possible.
We emphasize that random sampling cannot have theoretical guarantees since the random sample can be corrupted (similarly as in the example in Section 1.1).
Thus some outlier detection steps (such as our algorithms) are required.
Predictor Description.
We use the following predictors in our experiments.
1) Nearest neighbor : We use this predictor for the Oregon dataset.
We ﬁnd the best clustering of the node embeddings in Graph #1 .
In practice, this means running many steps of Lloyd’s algorithm until convergence after initial seeding by kmeans++ .
Our predictor takes as input a point in Rrepresenting a node embedding of any of the later 8 graphs and outputs the label of the closest node in Graph #1 . 2) Noisy predictor.
This is the main predictor for PHY.
We form this predictor by ﬁrst ﬁnding the best k -means solution on our datasets.
This again means initial seeding by kmeans++ and then many steps of Lloyd’s algorithm until convergence.
We then randomly corrupt the resulting labels by changing them to a uniformly random label independently with error probability ranging fromto 1 .
We report the cost of clustering using only these noisy labels versus processing these labels using Algorithm 1. 3) Neural network.
We use a standard neural network architecture (ResNet18) trained on the training portion of the CIFAR-10 dataset as the oracle for the testing portion which we use in our experiments.
We used a pretrained model obtained from Huy (2020).
Note that the neural network is predicting the class of the input image.
However, the class value is highly correlated with the optimal k -means cluster group.
Summary of results.
Our experiments show that our algorithm can leverage predictors to signiﬁcantly improve the cost of k -means clustering and that good predictors can be easily tailored to the data at hand.
The cost of k -means clustering reduces signiﬁcantly after applying our algorithm compared to just using the predictor labels for two of our predictors.
Lastly, the quality of the predictor remains high for the Oregon dataset even though the later graphs have changed and “moved away” from Graph #1 .
Selecting α in Algorithm 2.
In practice, the choice of α to use in our algorithm depends on the given predictor whose properties may be unknown.
Since our goal is to minimize the k -means clustering objective (1) [Figure] we can simply pick the ‘best’ value.
To do so, we iterate over a small range of possible α from .
01 to .
15 in Algorithm 2 with a step size of 0 . 01 and select the clustering that results in the lowest objective cost.
The range is ﬁxed for all of our experiments.
(See Paragraph [Figure]
R ESULTS
[Figure]
We ﬁrst compare our algorithm with Graph #1 as the predictor against various baselines.
This is shown in Figures 1(a) and Figure 1(b).
In the k = 10 case, Figure 1(a) shows that the predictor returns a clustering better than using just the kmeans++ seeding, which is normalized to have cost 1 . 0 .
This is to be expected since the subsequent graphs represent a similar network as
Graph #1 , just sampled later in time.
However, the clustering improves signiﬁcantly after using our algorithm on the predictor labels as the average cost drops by 55 % .
We also see that using our algorithm after kmeans++ is also sufﬁcient to give signiﬁcant decrease in clustering cost.
Lastly, [Figure] random sampling also gives comparable results.
This can be explained because we are iterating over a large range of subsampling probabilities for random sampling.
In the k = 25 case, Figure 1(b) shows that the oracle performance degrades and is worse than the baseline in 5 of the 8 graphs.
However our algorithm again improves the quality of the clustering over the oracle across all graphs.
Using kmeans++ as the predictor in our algorithm also improves the cost of clustering.
The performance of random sampling is also worse.
For example in Graph #3 for k = 25 , it performed the worst out of all the tested algorithms.
Our algorithm also remains competitive with kmeans++ seeding even if the predictor for the Oregon dataset is highly corrupted.
We consider a later graph, Graph #5 , and corrupt the labels of the predictor randomly with probability q ranging from 1% to 25% for the k = 10 case in Figure 1(c).
While the cost of clustering using just the predictor labels can become increasingly worse, our algorithm is able to sufﬁciently “clean” the predictions.
In addition, the cost of random sampling also gets worse as the corruptions increase, implying that it is much more sensitive to noise than our algorithm.
The qualitatively similar plot for k = 25 is given in the supplementary section.
Note that in spectral clustering, one may wish to get a mapping to R d for d > 2 .
We envision that our results translate to those settings as well since having higher order spectral information only results in a stronger predictor.
We continue the discussion on the PHY and CIFAR-10 datasets in Section E.
Comparison to Lloyd’s Heuristic.
In Section E.3, we provide additional results on experiments using Lloyd’s heuristic.
In summary, we give both theoretical and empirical justiﬁcations for why our algorithms are superior to blindly following a predictor and then running Lloyd’s heuristic.
A CKNOWLEDGEMENTS
Zhili Feng, David P. Woodruf and Samson Zhou would like to thank partial support from NSF grant No. CCF- 181584, Ofﬁce of Naval Research (ONR) grant N00014-18-1-256, and a Simons
Investigator Award.
Sandeep Silwal was supported in part by a NSF Graduate Research Fellowship
Program.
A A PPENDIX
Theorem A.1 (Chernoff Bounds) .
Let X 1 , . . .
, X n be independent random variables taking values in { 0 , 1 } .
Let X = (cid:80) ni
=1 X
i denote their sum and let µ = E [ X ] denote the sum’s expected value.
Then for any δ ∈ (0 , 1) and t > 0 ,
[Equation]
For any δ > 0 ,
Pr [ X ≥
[Equation]
) µ ] ≤ e −
[Equation]
Furthermore,
[Equation]
A. P ROOF OF T HEOREM 2.
We ﬁrst prove Theorem 2.1, which shows that Algorithm 1 provides a (1 + α ) -approximation to the optimal k -means clustering, but uses suboptimal time compared to a faster algorithm we present in
Section 3.
All omitted proofs of lemmas appear in Section A.2.
We ﬁrst show that for each coordinate, the empirical center for any (1 − α ) -fraction of the input points provides a good approximation to the optimal k -means clustering cost.
Lemma A.2.
Let P, Q ⊆
R be sets of points on the real line such that | P | ≥ (1 − α ) n and | Q | ≤ αn .
Let X
= P ∪ Q , C P be the mean of P and C X be the mean of X .
Then cost( X, C P ) ≤ (cid:16) 1 +
[Equation]
We now show that a conceptual interval I ∗ ⊂ R with “small” length contains a signiﬁcant fraction of the true points.
Ultimately, we will show that the interval I computed in the “training” phase in C RD E ST has smaller length than I ∗ with high probability and yet I also contains a signiﬁcant fraction of the true points.
The main purpose of I ∗ (and eventually I ) is to ﬁlter out extreme outliers because the “testing” phase only considers points in I ∩ X 2 .
Lemma A.3.
For a ﬁxed set X ⊆ R , let C be the mean of X and σ 2 = 1 2 | X | (cid:80) x ∈ X ( x − C )
2 be the variance.
Then the interval I ∗ = (cid:104) C − σ √ α , C + σ √ α (cid:105) contains at least a (1 − 4 α ) fraction of the points in X .
Using Lemma A.3, we show that the interval I that is computed in the “training” phase contains a signiﬁcant fraction of the true points.
Lemma A.4.
Let m be a sufﬁciently large consatnt.
We have that
I
=
[ a, b ] contains at least a 1 − 6 α fraction of points of X 2 and b − a ≤ 2 σ/ √ α , with high probability, i.e., 1 − 1 / poly( m ) .
We next show that the optimal clustering on a subset obtained by independently sampling each input point provides a rough approximation of the optimal clustering.
That is, the optimal center is well-approximated by the empirical center of the sampled points.
Lemma A.5.
Let S be a set of points obtained by independently sampling each point of X ⊆ R d with probability p = 12 .
Let C be the optimal center of these points and C S be the empirical center of these points.
Conditioned on | S | ≥ 1 , then E [ C S ] = x and there exists a constant γ such that for η ≥ 1 and | X | > ηγkα ,
[Equation]
Pr [cost( X, C S ) > (1 + α ) cost( X, C )] < 1 / ( ηk ) .
Using Lemma A.2, Lemma A.4, and Lemma A.5, we justify the correctness of the subroutine
C RD E ST .
Lemma A.6.
Let α ∈ (10 log n/ √ n, 1 / 7) .
Let P, Q ⊆
R be sets of points on the real line such that | P | ≥ (1 − α )2 m and | Q | ≤ 2 αm , and X = P ∪ Q .
Let C be the center of P .
Then C RD E ST on input set X outputs a point C (cid:48) such that with probability at least 1 − 1 / ( ηk ) , cost( P, C (cid:48) ) ≤ (1 + [Equation] [Equation] [Equation]
+
[Equation]
Using C RD E ST as a subroutine for each coordinate, we now prove Theorem 2.1, justifying the correctness of Algorithm 1 by generalizing to all coordinates and centers and analyzing the runtime of Algorithm 1.
Proof of Theorem 2.1.
Since Π has label error rate λ ≤ α , then by deﬁnition of label error rate, at least a (1 − α ) fraction of the points in each cluster are correctly labeled.
Note that the kmeans clustering cost can be decomposed into the sum of the costs induced by the centers in each dimension.
Speciﬁcally, for a set C
= { C 1 , . . .
, C k } of optimal centers,
[Equation]
where S i is the set of points in X that are assigned to center C i .
For a particular i ∈ [ k ] , we have
[Equation]
where x j and ( C i ) j are the j -th coordinate of x and C i , respectively.
By Lemma A.6, the cost induced by C RD E ST for each dimension in each center C (cid:48)
i is a (1 + α ) - approximation of the total clustering cost for the optimal center C i in that dimension with probability 1 − 1 / ( ηk ) .
That is,
[Equation]
for each j ∈ [ d ] .
Thus, taking a sum over all dimensions j ∈ [ d ] and union bounding over all centers i ∈ [ k ] , we have that the total cost induced by Algorithm 1 is a (1 + 20 α ) -approximation to the optimal k -means clustering cost with probability at least 1 − 1 /η .
To analyze the time complexity of Algorithm 1, ﬁrst consider the subroutine C RD E ST .
It takes
O ( kdn ) time to ﬁrst split each of the points in each cluster and dimension into two disjoint groups.
Finding the smallest interval that contains a certain number of points can be done by ﬁrst sorting the points and then iterating from the smallest point to the largest point and taking the smallest interval that contains enough points.
This requires O ( n log n ) time for each dimension and each center, which results in O ( kdn log n ) total time.
Once each of the intervals is found, computing the approximate center then takes O ( kdn ) total time.
Hence, the total running time of Algorithm 1 is
O ( kdn log n ) .
A. P ROOF OF A UXILIARY L EMMAS
Lemma A.2.
Let P, Q ⊆
R be sets of points on the real line such that | P | ≥ (1 − α ) n and | Q | ≤ αn .
Let X
= P ∪ Q , C P be the mean of P and C X be the mean of X .
Then cost( X, C P ) ≤ (cid:16) 1 + α 1 − α 2 (cid:17) cost( X, C X ) .
Proof.
Suppose without loss of generality, that C X = 0 and C P ≤ 0 , so that C Q ≥ 0 , where C Q is the mean of Q .
Then it is well-known, e.g., see Inaba et al.
, that cost( X, C P ) =
cost( X, C X ) + | X | · | C P − C X
| 2 .
Hence, it sufﬁces to show that | X | · | C P − C X | 2 ≤ α (1 − α ) 2 cost( X, C X ) .
Since C X =
we have | P |· C P = −| Q |· C Q , with | P | ≥ (1 − α ) n and | Q | ≤ αn .
Let | P | =
(1 − (cid:37) ) n and | Q | =
(cid:37)n for some (cid:37) ≤ α .
Thus, C Q = − 1 − (cid:37)(cid:37) · C P .
By convexity, we thus have that
[Equation]
Therefore, we have
[Equation]
Thus,
[Equation]
as desired.
Lemma A.3.
For a ﬁxed set X ⊆ R , let C be the mean of X and σ 2 = 1 2 | X | (cid:80) x ∈ X ( x − C )
2 be the variance.
Then the interval I ∗ = (cid:104) C − σ √ α , C + σ √ α (cid:105) contains at least a (1 − 4 α ) fraction of the points in X .
Proof.
Note that any point x ∈ X \ I ∗ satisﬁes | x − C | 2 > σ 2 / (4 α ) .
Thus, if more than a 4 α fraction of the points of X are outside of I ∗ , then the total variance is larger than σ 2 , which is a contradiction.
For ease of presentation, we analyze λ
=
and we note that the analysis extends easily to general λ .
We now prove the technical lemma that we will use in the proof of Lemma A.8.
Lemma A.7.
We have
[Equation]
Proof.
Let m be sufﬁciently large.
A Chernoff bound implies that for a sufﬁciently large constant
C ,
[Equation]
Furthermore, (cid:88)
[Equation]
so the upper bound on the desired relation holds.
A similar analysis provides a lower bound.
Lemma A.4.
Let m be a sufﬁciently large consatnt.
We have that
I
=
[ a, b ] contains at least a 1 − 6 α fraction of points of X 2 and b − a ≤ 2 σ/ √ α , with high probability, i.e., 1 − 1 / poly( m ) .
Proof.
By Lemma A.3, I ∗ contains at least 2 m (1 − 4 α ) of the points in X .
Hence, by applying an additive Chernoff bound for t =
O ( √ m log m ) and for sufﬁciently large m , we have that the number of points in I ∗ ∩
X 1 is at least m (1 − 5 α ) with high probability.
Since I is the interval of minimal length with at least m (1 − 5 α ) points, then the length of I is at most the length of I ∗ .
Moreover, again applying Chernoff bounds, we have that the number of points in I ∩ X 2 is at least m (1 − 6 α ) .
More formally, suppose we have a set of 2 m points that we randomly partition into two sets Xand X 2 .
Consider any ﬁxed interval J that has at least 2 cm total points for c ≥ 1 − 5 α (note there are at most O ( m 2 ) intervals in total since our points are in one dimension).
Let J 1 and J 2 denote the number of points in J that are in X 1 and X 2 respectively.
By a Chernoff bound, we have that both J 1 and J 2 are at least mc (1 − α ) with high probability.
In particular, | J 1 − J 2 | ≤ αmc with high probability.
Thus by using a union bound, all intervals with at least cm total points satisfy the property that the number of points partitioned to X 1 and the number of points partitioned to
X 2 differ by at most αmc with high probability.
Conditioning on this event, I must also contain m (1 − 6 α ) points in X 2 since it contains at least m (1 − 5 α ) points in X 1 , as desired.
Lemma A.8.
Let S be a set of points obtained by independently sampling each point of X ⊆ R d with probability 12 , and let C S be the centroid of S .
Let x be the centroid of X .
Conditioned on | S | ≥ 1 , we have E [ C S ] = x , and there exists a constant γ such that
[Equation]
Proof.
We ﬁrst prove that E [ C S ]
= x .
Note that by the law of iterated expectations,
[Equation]
Let x i 1 , . . .
, x i | S | be a random permutation of the elements in S , so that for each 1 ≤ j ≤ | S | , we have E [ x i j ] = x .
Now conditioning on the size of S , we can write
[Equation]
Therefore,
[Equation]
and it follows that E [ C S ] = x .
To prove that
[Equation]
we again condition on | S | .
Suppose that | S | =
j .
Then,
[Equation]
Now let y i t
= x i t − x for all 1 ≤ t ≤ j .
Therefore,
[Equation]
Note that x i 1 is uniform over elements in X , so it follows that
[Equation]
Now if j ≥ 2 , we have that
[Equation]
since (cid:80) i y i = 0 by deﬁnition.
Hence,
E
[Equation]
S |≥
[Equation]
x − x (cid:107) 2 .
Now the probability that | S | = j for j ≥ 2 is precisely (cid:0) | X | j (cid:1) / 2 | X | , so we have
[Equation]
From Lemma A.7, we have that
[Equation]
for some constant c so it follows that
[Equation]
for some constant c (cid:48) .
For j = 1 , note that
[Equation]
Moreover, we have Pr [ | S | = 1] = | X | 2 | X | and Pr [ | S | = 0] = 1 2 | X | .
Thus from the law of total expectation, we have
[Equation]
E | S |≥ 2 (cid:2) (cid:107) C S − x (cid:107) 2 (cid:3) ≤ | X | 2 | X | · | X | (cid:88) x ∈ X (cid:0) (cid:107) x − x (cid:107) 2 (cid:1) + c (cid:48) | X | 2 · (cid:32)(cid:88) x ∈ X (cid:107) x − x (cid:107)(cid:33) ≤ γ | X | 2 · (cid:32)(cid:88) x ∈ X (cid:107) x − x (cid:107)(cid:33) for some constant γ , as desired.
Lemma A.9.
Let S be a set of points obtained by independently sampling each point of X ⊆ R d with probability p = 12 .
Let C be the optimal center of X and C S be the empirical center of S .
Let γ ≥ 1 be the constant from Lemma A.8.
Then for η ≥ 1 and | X | > ηγk α ,
Pr [cost( X, C S ) > (1 + α ) cost( X, C )] < 1 / ( ηk ) .
Proof.
By Lemma A.8 and Markov’s inequality, we have
Pr
[Equation]
We have (cid:88)
[Equation]
so that by Lemma A.
(cid:88)
[Equation]
X
[Equation]
with probability at least 1 − 1 ηk .
Hence for | X | ≥ ηγkα , the approximate centroid of each cluster induces a (1 + α ) -approximation to the cost of the corresponding cluster.
Lemma A.5.
Let S be a set of points obtained by independently sampling each point of X ⊆ R d with probability p = 12 .
Let C be the optimal center of these points and C S be the empirical center of these points.
Conditioned on | S | ≥ 1 , then E [ C S ] = x and there exists a constant γ such that for η ≥ 1 and | X | > ηγkα ,
[Equation]
Proof.
Lemma A.5 follows immediately from Lemma A.8 and Lemma A.9.
Lemma A.6.
Let α ∈ (10 log n/ √ n, 1 / 7) .
Let P, Q ⊆
R be sets of points on the real line such that | P | ≥ (1 − α )2 m and | Q | ≤ 2 αm , and X = P ∪ Q .
Let C be the center of P .
Then C RD E ST on input set X outputs a point C (cid:48) such that with probability at least 1 − 1 / ( ηk ) , cost( P, C (cid:48) ) ≤ (1 + 18 α )(1 + α ) (cid:16)
+ α (1 − α )(cid:17) cost( P, C ) .
Proof.
Let α ∈
(10 log n/ √ n, 1 / 7) .
Then from Lemma A.4, we have that I ∩ X contains at least (1 − 6 α )
m points of P ∩ X 2 and at most 2 αm points of Q in an interval of length 2 σ/ √ α , where σ
[Equation]
From Lemma A.2, we have that
[Equation]
P [Equation] C
[Equation]
P
[Equation]
C 1 )
[Equation]
where C 0 is the center of I ∩ P ∩ X 2 and C 1 is the center of P ∩ X 2 .
For sufﬁciently large m and from Lemma A.9, we have that cost( P, C 1 ) ≤ [Equation] 1 + α ) cost( P, C ) , with probability at least 1 − 1 / ( ηk ) .
Thus, it remains to show that cost( P, C (cid:48) ) ≤ (1 +
O ( α )) cost( P, C 0 ) .
Since C 0 is the center of I ∩ P ∩ X 2 and C (cid:48) is the center of I ∩ X 2 , then we have
[Equation]
Since I has length 2 σ/ √ α , then q ∈ (cid:104) C 0 − 2 σ √
α , C 0 + 2 σ √
α (cid:105) .
Because | I ∩ P ∩ X 2 | ≥ (1 − 6 α ) m and | Q | = 2 αm , then for sufﬁciently small α , we have that
[Equation]
Note that we have cost( P, C (cid:48) )
=
cost( P, C 0 ) + | P | · | C 0 − C (cid:48) | 2 , so that cost
[Equation] [Equation]
cost
[Equation]
P
[Equation]
) + | P | · 36 ασ 2 .
Finally, σ 2 = 12 | P | · cost( P, C ) and cost( P, C ) ≤ cost( P, C 0 ) due to the optimality of C .
This implies cost( P, C
[Equation]
as desired.
Thus putting things together, we have
[Equation]
A. P ROOF OF T HEOREM 3.
We now give the proofs for optimal query complexity and runtime.
We ﬁrst require the following analogue to Lemma A.5:
Lemma A.10.
Let S be a set of points obtained by independently sampling each point of X ⊆ R d with probability p = min (cid:16) 1 , 100log k α | S | (cid:17) .
Let C be the optimal center of these points and C S be the empirical center of these points.
Conditioned on | S | ≥ 1 , then E [ C S ] = x and for | X | > γkα , [Equation] γ
[Equation]
(cid:32)(cid:88) x ∈ X (cid:107) x − x (cid:107)(cid:33) for some constant γ .
Lemma A.11.
For α ∈ (10 log n/ √ n, 1 / 7) , let Π be a predictor with error rate λ ≤ α/ 2 .
If each cluster has at least γk log k/α points, then Algorithm 3 outputs a (1 + 20 α ) -approximation to the k -means objective value with probability at least 3 / 4 .
Proof.
Since S samples each of points independently with probability proportional to cluster sizes given by Π , for a ﬁxed i ∈ [ k ] at least 90log k α points with label i are sampled, with probability at least 1 − 1 k 4 from Chernoff bounds.
Let γ 1 , . .
, γ k be the empirical means corresponding to each of the sampled points with labels 1 , . . .
, k , respectively, and let Γ
= { γ 1 , . . .
, γ k } .
Let C 1 , . .
C k be centers of a (1 + α ) -approximate optimal solution C with corresponding clusters X 1 , . . .
, X k .
By
Lemma A.10, we have that
[Equation]
where p = min (cid:16) 1 , 100log k α | S | (cid:17) .
By Markov’s inequality, we have that
[Equation]
[Equation]
[Equation]
with probability at least 0 . 99 .
Similar to the proof of Lemma A.9, we use the identity
[Equation]
Hence, we have that cost( X, Γ 0 ) ≤ (1 + α ) · cost( X, C ) , with probability at least 0 .
Suppose Π has error rate λ ≤ α and each error chooses a label uniformly at random from the k possible labels.
Then by deﬁnition of error rate, at most α/ 2 fraction of the points are erroneously labeled for each cluster.
Each cluster in the optimal k -means clustering of the predictor Π has at least n/ ( ζk ) points, so that at least a (1 − α ) fraction of the points in each cluster are correctly labeled.
Thus, by the same argument as in the proof of Lemma A.6, we have that Algorithm 1 outputs a set of centers C 1 , . . .
, C k such that for Γ =
{ C 1 , . . .
, C k } , we have cost( X, Γ
[Equation]
≤ (1 + 18 α
[Equation]
· cost
[Equation]
, with sufﬁciently large probability.
Let E be the event that cost( X, Γ) ≤ (1 + α )(1 + 18 α ) (cid:16) 1 − α (1 − α )(cid:17) · cost( X, C ) , so that Pr [ E ] ≥ 1 − 1 / poly( k ) .
Conditioned on E , let X 1 be the subset of X that is assigned the correct label by Π , and let X 2 be the subset of X assigned the incorrect label.
For each point x ∈ X 1 assigned the correct label (cid:96) x by Π , the closest center to x in Γ is C (cid:96) x , so Algorithm 3 will always label x with (cid:96) x .
Thus, cost( X 1 , Γ) ≤ cost( X, Γ) ≤ (1 + α )(1 + 18 α ) (cid:18) 1 − α (1 − α )(cid:19) · cost( X, C ) , conditioned on E .
On the other hand, if x ∈ X 2 is assigned an incorrect label (cid:96) x by Π , then the (2 , r ) -approximate nearest neighbor data assigns the label p x to x , where φ ( C p x ) is the closest center to φ ( x ) in the projected space.
Recall that φ is the composition map φ 1 ◦ φ 2 , where φ 1 has a terminal dimension reduction with distortion 5 / 4 , and φ 2 is a random JL linear map with distortion 5 / 4 .
Thus the distance between x and C p x is a 2 -approximation between x and its closest center C i .
Hence, by assigning all points x to their respective centers C p x , we have d ( x, C p x ) ≤ 2 cost( x, Γ) .
Since each point x ∈ X is assigned the incorrect label with probability λ ≤ α/ 2 , the expected cost of the labels assigned to X 2 is α cost( X, Γ) .
By Markov’s inequality, the cost of the labels assigned to Xis at most 10 α cost( X, Γ) < 10 α (1 + α ) cost( X, C ) , with probability at least 1 − 15 , conditioned on E .
Therefore by a union bound, the total cost is at most (1+20 α ) · cost( X, C ) , with probability at least 3 / 4 .
We need the following theorems on the quality of the data structures utilized in Algorithm 3.
Theorem A.12.
Makarychev et al.
For every set C ⊂ R d of size k , a parameter 0 < α <and the standard Euclidean norm d ( · , · ) , there exists a terminal dimension reduction f : C → R d (cid:48) with distortion (1 + α ) , where d (cid:48)
=
O (cid:16) log k α 2 (cid:17) .
The dimension reduction can be computed in polynomial time.
Theorem A.13.
Indyk & Motwani (1998);
Har-Peled et al.
; Andoni et al.
For α > 0 , there exists a (1+ α, r ) -ANN data structure over R equipped with the standard Euclidean norm that achieves query time O (cid:16) d · log n α 2 (cid:17) and space S :
= O (cid:0) 1 α 2 log 1 α + d ( n + q ) (cid:1) , where q := log n α 2 .
The runtime of building the data structure is O ( S + ndq ) .
We now prove Theorem 3.4.
Theorem 3.4.
Let α ∈
(10 log n/ √ n, 1 / 7) , Π be a predictor with label error rate λ ≤ α , and γ ≥be a sufﬁciently large constant.
If each cluster in the optimal k -means clustering of the predictor has at least γk log k/α points, then Algorithm 3 outputs a (1 + 20 α ) -approximation to the k -means objective with probability at least 3 / 4 , using O ( nd log n + poly( k, log n )) total time.
Proof.
The approximation guarantee of the algorithm follows from Lemma A.11.
To analyze the running time, we ﬁrst note that we apply a JL matrix with dimension O (log n ) to each of the n points in R d , which uses O ( nd log n ) time.
As a result of the JL embedding, each of the n points has dimension O (log n ) .
Thus, by Theorem A.12, constructing the terminal embedding uses poly( k, log n ) time.
As a result of the terminal embedding, each of the k possible centers has dimension O (log k ) .
Hence, by Theorem A.13, constructing the (2 , r ) -ANN data structure for the k possible centers uses
O ( k log 2 k ) time.
Subsequently, each query to the data structure uses O (log 2 k ) time.
Therefore, the overall runtime is O ( nd log n + poly( k, log n )) .
A. R EMARK ON T RULY - POLYNOMIAL TIME ALGORITHMS VS . PTAS/PRAS.
Remark A.14.
We emphasize that the runtime of our algorithm in Theorem 2.1 is truly polynomial in all input parameters n, d, k and 1 /α (and even near-linear in the input size nd ).
Although there exist polynomial-time randomized approximation schemes for k -means clustering, e.g., Inaba et al.
; Feldman et al.
; Kumar et al.
(2004), their runtimes all have exponential dependency on k and 1 /α , i.e., 2 poly( k, 1 /α ) .
However, this does not sufﬁce for many applications, since k and 1 /α should be treated as input parameters rather than constants.
For example, it is undesirable to pay an exponential amount of time to linearly improve the accuracy α of the algorithm.
Similarly, if the number of desired clusters k = O (log 2 n ) , then the runtime would be exponential.
Thus we believe the exponential improvement of Theorem 2.1 over existing PRAS in terms of k and 1 /α is signiﬁcant.
A. R EMARK ON P OSSIBLE I NSTANTIATIONS OF P REDICTOR
Remark A.15.
We can instantiate Theorem 2.1 with various versions of the predictor.
Assume each cluster in the (1 + α ) -approximately
optimal k -means clustering of the predictor has size at least n/ ( ζk ) for some tradeoff parameter ζ ∈ [1 , ( √ n ) / (8 k log n )] .
Then the clustering quality and runtime guarantees of Theorem 2.1 hold if the predictor Π is such that
Π outputs the right label for each point independently with probability 1 − λ
and otherwise outputs a random label for λ ≤ O ( α/ζ ) , 2.
Π outputs the right label for each point independently with probability 1 − λ
and otherwise outputs an adversarial label for λ ≤ O ( α/ ( kζ )) .
In addition, if the predictor Π outputs a failure symbol when it fails, then for constant ζ > 0 , there exists an algorithm (see supplementary material) that outputs a (1+ α ) -approximation to the k -means objective with probability at least 2 / 3 , even when Π has failure rate λ = 1 − 1 / poly( k ) .
Note that this remark (but not Theorem 2.1) assumes that each of the k clusters in the (1 + α ) -approximately optimal clustering has at least nζk points.
This is a natural assumption that the clusters are “roughly balanced” which often holds in practice, e.g., for Zipﬁan distributions.
B D ELETION P REDICTOR
In this section, we present a fast and simple algorithm for k -means clustering, given access to a label predictor Π with deletion rate λ .
That is, for each point, the predictor Π either outputs a label for the point consistent with an optimal k -means clustering algorithm with probability λ , or outputs nothing at all (or a failure symbol ⊥ ) with probability 1 − λ .
Since the deletion predictor fails explicitly, we can actually achieve a (1 + α ) -approximation even when λ
= 1 − 1 poly( k ) .
Our algorithm ﬁrst queries all points in the input X .
Although the predictor does not output the label for each point, for each cluster C i with a sufﬁciently large number of points, with high probability, the predictor assigns at least λ 2 | C i | points of C i to the correct label.
We show that if | C i | =
Ω (cid:0) kα (cid:1) , then with high probability, the empirical center is a good estimator for the true center.
That is, the kmeans objective using the centroid of the points labeled i is a (1+ α ) -approximation to the k -means objective using the true center of C i .
We give the full details in Algorithm 4.
To show that the empirical center is a good estimator for the true center, recall that a common approach for mean estimation is to sample roughly an O (cid:0) 1 α 2 (cid:1) number of points uniformly at random with replacement.
The argument follows from observing that each sample is an unbiased estimator of the true mean, and repeating O (cid:0) 1 α 2 (cid:1) times sufﬁciently upper bounds the variance.
Observe that the predictor can be viewed as sampling the points from each cluster without replacement .
Thus, for sufﬁciently large cluster sizes, we actually have a huge number of samples, which intuitively should sufﬁciently upper bound the variance.
Moreover, the empirical mean is again an unbiased estimator of the true mean.
Thus, although the above analysis does not quite hold due to dependencies between the number of samples and the resulting averaging term, we show that the above intuition does hold.
then 7: (cid:96) x
← arg min d ( x, c i ) 8:
Assign label (cid:96) x to x .
9: end if 10: end for
We ﬁrst show that independently sampling points uniformly at random from a sufﬁciently large point set guarantees a (1 + α ) -approximation to the objective cost.
Inaba et al.
; Ailon et al.
(2018) proved a similar statement for sampling with replacement.
It remains to justify the correctness of Algorithm 4 by arguing that with high probability, the overall k -means cost is preserved up to a (1+ α ) -factor by the empirical means.
We also analyze the running time of Algorithm 4.
Theorem B.1.
If each cluster in the optimal k -means clustering of the predictor Π has at least 3 kα points, then Algorithm 4 outputs a (1 + α ) -approximation to the k -means objective with probability at least 23 , using O ( kdn ) total time.
Proof.
We ﬁrst justify the correctness of Algorithm 4.
Suppose each cluster in the optimal k -means clustering of the predictor Π has at least 3 kα points.
Let C
= { c 1 , . .
c k } be the optimal centers selected by Π and let C S
= { c (cid:48) 1 , . . .
, c (cid:48) k } be the empirical centers chosen by Algorithm 4.
For each i ∈ [ k ] , let C i be the points of X that are assigned to C i by the predictor Π .
By Lemma A.9 with η = 3 , the approximate centroid of a cluster induces a (1 + α ) -approximation to the cost of the corresponding cluster so that
[Equation]
with probability at least 1 − 13 k .
Taking a union bound over all k clusters, we have that
[Equation]
with probability at least 23 .
Equivalently, cost( X, C ) ≤ (1 + α ) cost( X, C S ) .
To analyze the running time of Algorithm 4, observe that the estimated centroids for all labels can be computed in O ( dn ) time.
Subsequently, assigning each unlabeled point to the closest estimated centroid uses O ( kd ) time for each unlabeled point.
Thus, the total running time is O ( kdn ) .
C k - MEDIAN C LUSTERING
We ﬁrst recall that a well-known result states that the geometric median that results from uniformly sampling a number of points from the input is a “good” approximation to the actual geometric median for the 1 -median problem.
Theorem C.1.
Krauthgamer (2019)
Given a set P of n points in R d , the geometric median of a sample of O (cid:0) dα 2 log dα (cid:1) points of P provides a (1 + α ) -approximation to the 1 -median clustering problem with probability at least 1 − 1 / poly( d ) .
Note that we can ﬁrst apply Theorem A.12 to project all points to a space with dimension
O (cid:0) 1 α 2 log kα (cid:1) before applying Theorem C.1.
Instead of computing the geometric median, we recall the following procedure that produces a (1 + α ) -approximation to the geometric median.
Theorem C.2.
Cohen et al.
There exists an algorithm that outputs a (1 + α ) -approximation to the geometric median in O (cid:0) nd log 3 nα (cid:1) time.
We give our algorithm in full in Algorithm 5.
Theorem C.3.
For α ∈ (0 , 1) , let Π be a predictor with error rate
λ
= O (cid:16)
α 4 k log kα log log kα (cid:17) .
If each cluster in the optimal k -median clustering of the predictor has at least n/ ( ζk ) points, then there exists an algorithm that outputs a (1 + α ) -approximation to the k -median objective with probability at least 1 − 1 / poly( k ) , using O ( nd log 3 n + poly( k, log n )) total time.
Proof.
Observe that Algorithm 5 samples O (cid:0) 1 α 4 log 2 kα (cid:1) points for each of the clusters labeled i , with i ∈ [ k ] .
Thus Algorithm 5 samples O (cid:0) kα 4 log 2 kα (cid:1) points in total.
For λ =
O (cid:16)
[Equation]
with a sufﬁciently small constant, the expected number of incorrectly labeled points sampled by
Algorithm 5 is less than 1 32 .
Thus, by Markov’s inequality, the probability that no incorrectly labeled points are sampled by Algorithm 5 is at least 34 .
Conditioned on the event that no incorrectly labeled points are sampled by Algorithm 5, then by Theorem C.1, the empirical geometric median for each cluster induces a (cid:0) 1 + α 4 (cid:1) -approximation to the optimal geometric median in the projected space.
Hence the set of k empirical geometric medians induces a (cid:0) 1 + α 4 (cid:1) -approximation to the optimal k -median clustering cost in the projected space.
Since the projected space is the result of a terminal embedding, the set of k empirical geometric medians for the sampled points in the projected space induces a k -median clustering cost that is a (cid:0) 1 + α 4 (cid:1) -approximation to the k -median clustering cost induced by the set of k empirical geometric medians for the sampled points in the original space.
Taking the set of k empirical geometric medians for the sampled points in the original space induces a (cid:0) 1 + α 4 (cid:1) 2 -approximation to the k -median clustering cost.
We take a (cid:0) 1 + α 4 (cid:1) -approximation to each of the geometric medians.
Thus for sufﬁciently small α , Algorithm 5 outputs a (1 + α ) - approximation to the k -median clustering problem.
To embed the points into the space of dimension O (cid:0) 1 α 2 log kα (cid:1) , Algorithm 5 spends O ( nd log n ) total time.
By Theorem C.2, it takes O ( nd log 3 n ) total time to compute the approximate geometric medians.
D L OWER B OUNDS
MAX-E3-LIN-2 is the optimization problem of maximizing the number of equations satisﬁed by a system of linear equations of Z 2 with exactly 3 distinct variables in each equation.
E K -MAX-E3LIN-2 is the problem of MAX-E3-LIN-2 when each variable appears in exactly k equations.
Fotakis et al.
(2016) showed that assuming the exponential time hypothesis (ETH) (Impagliazzo & Paturi, 2001), there exists an absolute constant C 1 such that MAX k -SAT (and thus MAX k -CSP) instances with fewer than O ( n k − 1 ) clauses cannot be approximated within a factor of C 1 in time 2 O ( n 1 − δ ) for any δ > 0 .
As a consequence, the reduction by H ˚ astad (2001) shows that there exist absolute constants C 2 , C 3 such that E K -MAX-E3-LIN-2 with k ≥ C 2 cannot be approximated within a factor of C 3 in time 2 O ( n 1 − δ ) for any δ > 0 .
Hence, the reduction by Chleb´ık & Chleb´ıkov´a (2006) shows that there exists a constant C 4 such that approximating the minimum vertex cover of 4 -regular graphs within a factor of C 4 cannot be done in time 2 O ( n 1 − δ ) for any δ > 0 .
Thus the reduction by
Lee et al.
(2017) shows that there exists a constant C 5 such that approximating k -means within a factor of C 5 cannot be done in time 2 O ( n 1 − δ ) for any δ > 0 , assuming ETH.
Namely, the reduction of
Lee et al.
(2017) shows that an algorithm that provides a C 5 -approximation to the optimal k -means clustering can be used to compute a C 4 -approximation to the minimum vertex cover.
Theorem D.1.
If ETH is true, then there does not exist an algorithm A that takes a set S of n 1 −
δ log n vertices and ﬁnds a C 4 -approximation to the minimum vertex cover that contains S on a 4 -regular graph G , using 2 O ( n 1 − δ ) time for some constant δ ∈ (0 , 1] .
Proof.
Suppose by way of contradiction that there exists an algorithm A that takes a set S of n 1 −
δ log n vertices and ﬁnds a C 4 -approximation to the minimum vertex cover that contains S on a 4 -regular graph G , using 2 O ( n 1 − δ ) time for some constant δ ∈ (0 , 1] .
We claim that we can use A to create an overall algorithm that violates ETH.
Indeed, suppose we guess each subset of n 1 − δ log n ver- tices and which vertices of the subset are in the cover.
There are (cid:0)
[Equation]
( en δ log n ) n 1 − δ / log n · 2 n/ log n such combinations of vertices.
For each guess, we then run the purported algorithm A that uses 2 O ( n 1 − δ ) time.
Thus we can identify a C 4 -approximation to the minimum vertex cover in time
[Equation]
which would contradict ETH.
Finally, we show the query complexity of Algorithm 3 is nearly optimal.
Lee et al.
(2017) constructed an instance of k -means that cannot be approximated within a factor of 1 . 0013 in polynomial time.
The reduction of Lee et al.
(2017) creates 4 n points in R 3 n that must be clustered by O ( n ) centers and an algorithm that provides a C 5 -approximation to the optimal k -means clustering can be used to compute a C 4 -approximation to the minimum vertex cover.
Thus, there exists a constant
C 5 such that approximating k -means within a factor of C 5 cannot be done in time 2 O ( n 1 − δ ) for any δ > 0 , assuming ETH.
Theorem 3.5.
For any δ ∈ (0 , 1] , any algorithm that makes O (cid:16) k 1 − δ α log n (cid:17) queries to the predictor with label error rate α cannot output a (1 + Cα ) -approximation to the optimal k -means clustering cost in time 2 O ( n 1 − δ ) time, assuming the Exponential Time Hypothesis.
Proof.
Let α be a ﬁxed constant such that α < C 5 .
Given an instance I of a k -means clustering constructed from the reduction of Lee et al.
(2017), the optimal clustering cost is Ω( n ) and k 1 = Ω( n ) .
We embed this instance into a k -means clustering by adding an additional Ω( n ) points arbitrarily far from I , so that the additional points contribute Ω( n/α ) cost upon partitioning into k
= Ω( n ) clusters.
We set k = k 1 + k 2 .
In summary, the resulting instance has O ( n ) points and k = Ω( n ) .
The optimal solution has cost
O ( n/α ) cost so that I contributes an Ω( α ) fraction of the cost.
By querying O (cid:16) k 1 − δ α log n (cid:17) points with sufﬁciently small constant, at most O (cid:16) k 1 − δ α log n (cid:17) of the cluster centers in I will be revealed by the construction of I .
Each center corresponds to a selected vertex in the corresponding vertex cover in the reduction from minimum vertex cover on 4 -regular graphs.
Hence, in the corresponding vertex cover instance, at most O (cid:16) k 1 − δ α log n (cid:17) vertices are revealed.
Thus by Theorem D.1, any algorithm running in 2 O ( n 1 − δ ) time cannot determine a C 5 -approximation to the optimal k -means clustering cost on I , as it would correspond to a C 4 -approximation to the optimal vertex cover, assuming ETH.
Since I induces an Ω( α ) fraction of the total clustering cost, it follows that any algorithm that makes
O (cid:16) k 1 − δ α log n (cid:17) queries cannot output a (1+ Cα ) -approximation to the optimal k -means clustering cost in time 2 O ( n 1 − δ ) time, assuming ETH.
E A DDITIONAL E XPERIMENTAL R ESULTS E. O MITTED D ISCUSSION
In this section we continue the discussion from the experimental section of the main body.
In Figure 3(a), we show the qualitatively similar version of Figure 2( c ) of the main body for the case of k = 25 .
In Figure 3(b), we display the qualitatively similar version of Figure 3( a ) of the main body for the k = 25 case of dataset PHY. PHY.
We use the noisy predictor for this dataset.
We see in Figure 2(a) that as the corruption percentage rises, the clustering given by just the predictor labels can have increasingly large cost. Nevertheless, even if the clustering cost of the corrupted labels is rising, the cost decreases signiﬁcantly after applying Algorithm 1 by roughly a factor of 3x . Indeed, we see that our algorithm can beat the kmeans++ seeding baseline for q as high as 50% . Just as in Figure 1(c), random sampling is sensitive to noise.
Lastly, we also remain competitive with the purple line which uses the labels output by kmeans++ as the predictor in our algorithm (no corruptions added).
The qualitatively similar plot for k = 25 is given in the supplementary material.
CIFAR-10.
The cost of clustering on CIFAR-10 using only the predictor, the predictor with our algorithm, random sampling, and kmeans++ as the predictor for our algorithm were 0 .
721 , and 0 .
640 , respectively, where 1 . 0 represents the cost of kmeans++ .
The neural network was very accurate ( ∼ 93% ) in predicting the class of the input points which is highly correlated with the optimal k -means clusters.
Nevertheless, our algorithm improved upon this highly precise predictor.
Note that using kmeans++ as the predictor for our algorithm resulted in a clustering that was 13 % better than the one given by the neural network predictor.
This highlights the fact that an approximate clustering combined with our algorithm can be competitive against a highly precise predictor, such as a neural network, even though creating the highly accurate predictor can be expensive.
Indeed, obtaining a neural network predictor requires prior training data and also the time to train the network.
On the other hand, using a heuristic clustering as a predictor is extremely ﬂexible and can be applied to any dataset even if no prior training dataset is available ( 50 , 000 test images were required to train the neural network predictor but kmeans++ as a predictor requires no test images), in addition to considerable savings in computation.
For example, the time taken to train the particular neural network we used was approximatelyminutes using the optimized PyTorch library (see training details under the “Details Report” section in Huy (2020)).
In general, the time can be much longer for more complicated datasets.
On the other hand, our unoptimized algorithm implementation which used the labels of a sample run of kmeans++ was still able to achieve a better clustering than the neural network predictor with α = 0 .
01 in Algorithm 2 in 4 . 4 seconds. In conclusion, we can achieve a better clustering by combining a much weaker predictor with our algorithm with the additional beneﬁt of using a more ﬂexible and computationally inexpensive methodology. We also conducted an experiment where we only use a small fraction of the predictor labels in our algorithm.
We select a p fraction of the images, query their labels, and run our algorithm on only these points. We then report the cost of clustering on the entire dataset as p ranges from 1% to 100% . Figure 2(b) shows the percentage increase in clustering cost relative to querying the whole dataset is quite low for moderate values of p but increasingly worse as p becomes smaller. This suggests that the quality of our algorithm does not suffer drastically by querying a smaller fraction of the dataset.
[Figure] E. S YNTHETIC D ATASET We use a dataset of 10010 points in R d for d = 10 3 created using the kmeans++ lower bound construction presented in Arthur & Vassilvitskii (2007).
The dataset consists of 10 well separated clusters in R d : let e i denote the basis vectors.
Our dataset is { 1000 e i } ∪ { 1000 e i + e j } for all 1 ≤ i ≤ 10 , 1 ≤ j ≤ 1000 .
From the description of the dataset, we can explicitly calculate the optimal clustering and its cost.
Our predictor for this dataset was to take the optimal clustering and randomly change each label with probability 1 / 2 to another uniformly random label.
Empirically, kmeans++ seeding returned a clustering that had cost at least 1.9x the optimal clustering.
Furthermore, using just the predictor labels na¨ıvely resulted in a clustering with cost up to ﬁve orders of magnitude larger than the optimal clustering.
In contrast, our algorithm was able to precisely recover the true clustering after processing the predictor outputs.
In addition, applying our algorithm using the labels of a sample run of kmeans++ was also able to precisely recover the true clustering.
E. C OMPARISON TO L LOYD ’ S H EURISTIC
We give both theoretical and empiricial justiﬁcations for why our algorithms could be superior to blindly following a predictor and then running Lloyd’s heuristic.
[Figure]
Empirical Comparison.
We ﬁrst compared Lloyd’s algorithm on kmeans++ seeding to our algorithm with the predictor on the PHY dataset.
The predictor is a noisy predictor that has corruption level 50% as described in Section 4 so that outputting the clustering from the predictor alone has cost 1.9x the average kmeans++ cost.
Hence, it is clear that the predictor is much worse than kmeans++ , yet our algorithm using the predictor (horizontal line in Figure 4(a)) is much better than kmeans++ and Lloyd’s algorithm (orange line in Figure 4(a)).
Next, we took the noisy predictor for the PHY dataset with corruption level 50% and repeatedly applied Lloyd’s algorithm.
We observe that even with ∼ 5 Lloyd’s iterations (Figure 4(b)), the clustering cost does not seem to improve upon kmeans++ , much less the clustering output by simply applying our algorithm to the noisy predictions (Figure 4(a)).
Finally, we compared Lloyd’s algorithm on kmeans++ seeding to Lloyd’s algorithm on the seeding output by our algorithm (with kmeans++ as the predictor) on the CIFAR-10 dataset, similar to the experiments you suggested by Lattanzi & Sohler (2019).
Lloyd’s algorithm on the seeding output by our algorithm exhibits superior performance than Lloyd’s algorithm on kmeans++ (Figure 4(c)), which is consistent with our previous experiments showing that our algorithm improves upon kmeans++ .
This further strengthens our claim that our algorithm and methodology with provable worst-case guarantees can be applied in conjunction with heuristics such as Lloyd’s that do not have provable worst-case guarantees.
Moreover, Figure 4(c) indicates that our approach may be more advantageous than just running kmeans++ with heuristics.
Note that a na¨ıve implementation of
Lloyd’s algorithm is O ( ndk ) time while our algorithm can be implemented in nearly linear time.
We emphasize that all of our above experiments use a noisy predictor with corruption level 50% as input to our algorithm.
Our experiments exhibit even better behavior when a clustering produced by kmeans++ is used for a predictor as input to our algorithm.
Combined with our other experiments in Section 4, this gives empirical evidence that there exist many scenarios in which running our algorithm with an erroneous predictor is advantageous.
Theoretical Comparison.
We now provide an example that demonstrates why blindly following the predictor and then running Lloyd’s heuristic would run into issues.
We emphasize that it is wellknown e.g., see Dasgupta (2003);
Har-Peled & Sadri (2005); Arthur & Vassilvitskii (2006); Vattani (2011), that Lloyd’s algorithm can take a large number of steps to converge.
In particular, Vattani (2011) shows that an exponential number of Lloyd iterations can be required for the algorithm to converge to the optimal solution.
Nevertheless, we offer the following concrete answer:
We describe a simple set of points that guarantees Lloyd’s algorithm will fail.
This is based on the example given in Har-Peled & Sadri (2005) and is also conceptually similar to the example given in 1.1.
Consider 4 n points on the real line x 1 , . . .
, x 2 n , y 1 , . .
y 2 n so that y 1 ≤ . .
≤ y 2 n ≤
A < B ≤ x 1 ≤ . . .
≤ x 2 n where B − A is large.
Suppose k = 2 in which case the optimal clustering groups all the y i points together and all the x i points together as two separate clusters.
Suppose the predictor initially gives label “1” to points y 1 , . .
y n and gives label “2” to points y n +1 , .
, x 1 , . . .
, x 2 n , so that the predictor has corruption level λ = 1 / 2 .
Then our algorithm that uses this predictor will get a constant-factor approximation in only one iteration.
However, since
B − A can be arbitrarily large without affecting the optimal clustering cost, blindly listening to the predictor will give a worse clustering.
Furthermore, Theorem 2 . 1 in Har-Peled & Sadri (2005) implies that Lloyd’s will take Θ( n ) iterations to converge if initialized using this predictor.
Note that even a single (na¨ıve) iteration of Lloyd’s algorithm already uses O ( ndk ) time while our algorithm only uses O ( nd )
+ poly( k, 1 /α ) time.
Note that there are also more complex examples in higherdimensional spaces in the literature which provably have even worse convergence rates for Lloyd’s method.
E. C ONCLUSION
Although 1 . 07 -approximation for k -means clustering in polynomial time is NP-hard and a clustering consistent with the labels of any predictor with nonzero error can be arbitrarily bad, we give a (1 + α ) -approximation algorithm that uses the labels output by the predictor as “advice” and runs in nearly linear time.
We use a linear number of queries to the predictor, which can be improved to nearly optimal under natural assumptions about the cluster sizes.
Our results are well-supported by empirical evaluations and are an important step in demonstrating the power of learning-augmented algorithms beyond the limitations of classical algorithms for clustering-based applications.
F L EARNABILITY R ESULTS
In this section, we present formal learning bounds for learning a good predictor for the learningaugmented k -means problem.
Namely, we show that a good predictor can be learned efﬁciently if the problem instances are drawn from a particular distribution.
Our result is inspired by derived from data-driven algorithm design and utilizes the PAC learning framework.
Formally, our setting is the following.
Suppose there exists an underlying distribution D that generates independent k -means clustering instances, representing the case where similar instances of the k -means clustering problem are being solved.
Note that this setting also mirrors some of our experiments in Section 4, speciﬁcally in the case of our spectral clustering datasets which are derived from snapshots of a dynamic graph across time.
Our goal is to efﬁciently learn a good predictor f among some family of functions F .
The input to each predictor f is a clustering instance C and the output is a feature vector.
We assume that the each input instance C is encoded as a vector in R nd , that is, we think of nd as the size of the description of the input.
We also assume that the output of f is in n dimensions, which represents the prediction of the oracle.
For example in the case of k -means clustering, the output can be an integer in [ k ] for each point in C .
To select the “best” f , we need to formally deﬁne what we mean by best.
In most learning settings, a loss function is used to measure the quality of a solution.
Indeed, suppose we have a loss function
L : f × C → R which represents how well a predictor f performs on some input C .
For example, f can represent the cluster labels of a points in C and L outputs the cost of the clustering.
Alternatively,
L can represent an algorithm which uses the hints given by f and performs a clustering algorithm such as some number of steps of Lloyd’s heuristic.
Our goal is to learn the best function f ∈ F that minimizes the following objective:
[Equation]
We deﬁne f ∗ to be an optimal function f ∈ F , so that f ∗ = arg min E C ∼ D
[ L ( f, C )] .
We also assume that for each clustering instance C and each f ∈ F , f ( C ) and L ( f, C ) can be computed in time T ( n, d ) that should be interpreted as a (small) polynomial in n and d .
Our main result is the following.
Theorem F.1.
There exists an algorithm that uses poly( T ( n, d ) , 1 /ε ) samples and returns a ˆ f that satisﬁes
[Equation]
with probability at least 9 / 10 .
The above theorem is a PAC-style bound that shows only a small number of samples are needed in order to have a good probability of learning an approximately-optimal function ˆ f .
The algorithm to compute ˆ f is the following: we simply minimize the empirical loss after an appropriate number of samples are drawn, i.e., we perform empirical risk minimization.
This result is proven by Theorem
F.3.
Before introducing it, we need to deﬁne the concept of pseudo-dimension for a function class, which is the more familiar VC dimension, generalized to real functions.
Deﬁnition F.2 (Pseudo-Dimension, e.g., Deﬁnition 9 in Lucic et al.
Let X be a ground set and F be a set of functions from X to the interval [0 , 1] .
Fix a set S
= { x 1 , · · · , x n } ⊂ X , a set of reals numbers R = { r 1 , · · · , r n } with r i ∈ [0 , 1] and a function f ∈ F .
The set S f = { x i ∈ S | f
( x i ) ≥ r i } is called the induced subset of S formed by f and R .
The set S with associated values
R is shattered by F if |{ S f | f ∈ F}| = 2 n .
The pseudo-dimension of F is the cardinality of the largest shattered subset of X (or ∞ ).
The following theorem relates the performance of empirical risk minimization and the number of samples needed, to pseudo-dimension.
We specialize the theorem statement to our situation at hand.
For notational simplicity, we deﬁne G be the class of functions in F composed with L :
[Equation]
Furthermore, by normalizing, we can assume that the range of L is equal to [0 , 1] .
Theorem F.3 (Anthony & Bartlett (1999)) .
Let D be a distribution over problem instances C and
G be a class of functions g : C → [0 , 1] with pseudo-dimension d G .
Consider t i.i.d.
samples
C 1 , C 2 , . .
C t from D .
There is a universal constant c 0 , such that for any ε > 0 , if t ≥ c 0 · d G /ε 2 , then we have
[Equation]
for all g ∈ G with probability at least 9 / 10 .
The following corollary follows from the triangle inequality.
Corollary F.4.
Consider a set of t independent samples C 1 , . .
C t from D and let ˆ g be a function in G that minimizes 1 t (cid:80) ti =1 g ( C i ) .
If the number of samples t is chosen as in Theorem F.3, then
[Equation]
holds with probability at least 9 / 10 .
Therefore, the main challenge at hand is to bound the pseudo-dimension of our given function class
G .
To do so, we ﬁrst relate the pseudo-dimension to the VC dimension of a related class of threshold functions.
This relationship has been fruitful in obtaining learning bounds in a variety of works such as Lucic et al.
; Izzo et al.
Lemma F.5 (Pseudo-dimension to VC dimension, Lemma 10 in Lucic et al.
For any g ∈
G , let B g be the indicator function of the region on or below the graph of g , i.e., B g ( x, y ) =
sgn ( g ( x ) − y ) .
The pseudo-dimension of G is equivalent to the VC-dimension of the subgraph class
B G =
{ B g | g ∈ G} .
Finally, the following theorem relates the VC dimension of a given function class to its computational complexity, i.e., the time complexity of computing a function in the class.
Lemma F.6 (Theorem 8 . 14 in Anthony & Bartlett (1999)) .
Let h : R a × R b → { 0 , 1 } , determining the class
[Equation]
Suppose that any h can be computed by an algorithm that takes as input the pair ( θ, x ) ∈ R a × R b and returns h ( θ, x ) after no more than t of the following operations:
• arithmetic operations + , − , × , and / on real numbers, • jumps conditioned on >, ≥ , <, ≤ , = , and = comparisons of real numbers, and • output 0 , 1 ,
then the VC dimension of H is O ( a 2 t 2 + t 2 a log a ) .
Combining the previous results allows us prove Theorem F.1.
At a high level, we are instantiating
Lemma F.6 with the complexity of computing any function in the function class G .
Proof of Theorem F.1.
First by Theorem F.3 and Corollary F.4, it sufﬁces to bound the pseudodimension of the class G =
L ◦ F .
Then from Lemmas F.5, the pseudo-dimension of G is the VC dimension of threshold functions deﬁned by G .
Finally from Lemma F.6, the VC dimension of the appropriate class of threshold functions is polynomial in the complexity of computing a member of the function class.
In other words, Lemma F.6 tells us that the VC dimension of B G deﬁned in
Lemma F.5 is polynomial in the number of arithmetic operations needed to compute the threshold function associated to some g ∈ G .
By our deﬁnition, this quantity is polynomial in T ( n, d ) .
Hence, the pseudo-dimension of G is also polynomial in T ( n, d ) and our desired result follows.
Note that we can consider initializing Theorem F.1 with speciﬁc predictions.
If each function in the family of oracles we are interested can be computed efﬁciently, which is the case of the predictors we employ in our experiments, then Theorem F.1 assures us that we only require polynomially many samples to be able to learn a nearly optimal oracle.
Note that our result is in similar in spirit to the recent paper Dinitz et al.
They derive sample complexity learning bounds for a different algorithmic problem of computing matchings in a graph.
Since they specialize their analysis to a speciﬁc function class and loss function, their bounds are possibly tighter rather than the possibly loose polynomial bounds we have stated.
However, our analysis above is more general as it allows for a variety of predictors, as we employ in our experiments, and loss functions to measure the quality of the predictions.