Fair Clustering Through Fairlets

We study the question of fair clustering under the disparate impact doctrine, where each protected class must have approximately equal representation in every cluster. We formulate the fair clustering problem under both the k-center and the k-median objectives, and show that even with two protected classes the problem is challenging, as the optimum solution can violate common conventions—for instance a point may no longer be assigned to its nearest cluster center!

En route we introduce the concept of fairlets, which are minimal sets that satisfy fair representation while approximately preserving the clustering objective. We show that any fair clustering problem can be decomposed into first finding good fairlets, and then using existing machinery for traditional clustering algorithms. While finding good fairlets can be NP-hard, we proceed to obtain efficient approximation algorithms based on minimum cost flow.


From self driving cars, to smart thermostats, and digital assistants, machine learning is behind many of the technologies we use and rely on every day. Machine learning is also increasingly used to aid with decision making—in awarding home loans or in sentencing recommendations in courts of law (Kleinberg et al. , 2017a). While the learning algorithms are not inherently biased, or unfair, the algorithms may pick up and amplify biases already present in the training data that is available to them. Thus a recent line of work has emerged on designing fair algorithms.

The first challenge is to formally define the concept of fairness, and indeed recent work shows that some natural conditions for fairness cannot be simultaneously achieved (Kleinberg et al. , 2017b; Corbett-Davies et al. , 2017). In our work we follow the notion of disparate impact as articulated by Feldman et al. (2015), following the Griggs v. Duke Power Co. US Supreme Court case. Informally, the doctrine codifies the notion that protected attributes, such as race and gender, should not be explicitly used in making decisions, and the decisions made should not be disproportionately different for applicants in different protected classes. In other words, if an unprotected feature, for example, height, is closely correlated with a protected feature, such as gender, then decisions made based on height may still be unfair, as they can be used to effectively discriminate based on gender.

While much of the previous work deals with supervised learning, here we consider the most common unsupervised learning problem: clustering. In modern machine learning systems, clustering is often used


Figure 1: A colorblind k-center clustering algorithm would group points a, b, c into one cluster, and x, y, z into a second cluster, with centers at a and z respectively. A fair clustering algorithm, on the other hand, may give a partition indicated by the dashed line. Observe that in this case a point is no longer assigned to its nearest cluster center. For example x is assigned to the same cluster as a even though z is closer.

for feature engineering, for instance augmenting each example in the dataset with the id of the cluster it belongs to in an effort to bring expressive power to simple learning methods. In this way we want to make sure that the features that are generated are fair themselves. As in standard clustering literature, we are given a set X of points lying in some metric space, and our goal is to find a partition of X into k different clusters, optimizing a particular objective function. We assume that the coordinates of each point  x ∈ Xare unprotected; however each point also has a color, which identifies its protected class. The notion of disparate impact and fair representation then translates to that of color balance in each cluster. We study the two color case, where each point is either red or blue, and show that even this simple version has a lot of underlying complexity. We formalize these views and define a fair clustering objective that incorporates both fair representation and the traditional clustering cost; see Section 2 for exact definitions.

A clustering algorithm that is colorblind, and thus does not take a protected attribute into its decision making, may still result in very unfair clusterings; see Figure 1. This means that we must explicitly use the protected attribute to find a fair solution. Moreover, this implies that a fair clustering solution could be strictly worse (with respect to an objective function) than a colorblind solution.

Finally, the example in Figure 1 also shows the main technical hurdle in looking for fair clusterings. Unlike the classical formulation where every point is assigned to its nearest cluster center, this may no longer be the case. Indeed, a fair clustering is defined not just by the position of the centers, but also by an assignment function that assigns a cluster label to each input.

Our contributions. In this work we show how to reduce the problem of fair clustering to that of classical clustering via a pre-processing step that ensures that any resulting solution will be fair. In this way, our approach is similar to that of Zemel et al. (2013), although we formulate the first step as an explicit combinatorial problem, and show approximation guarantees that translate to approximation guarantees on the optimal solution. Specifically we:


(iii) Show that any fair clustering problem can be reduced to first finding a fairlet decomposition, and then using the classical (not necessarily fair) clustering algorithm;

(iv) Develop approximation algorithms for finding fair decompositions for a large range of fairness values, and complement these results with NP-hardness; and

(v) Empirically quantify the price of fairness, i.e., the ratio of the cost of traditional clustering to the cost of fair clustering.

Related work. Data clustering is a classic problem in unsupervised learning that takes on many forms, from partition clustering, to soft clustering, hierarchical clustering, spectral clustering, among many others. See, for example, the books by Aggarwal & Reddy (2013); Xu & Wunsch (2009) for an extensive list of problems and algorithms. In this work, we focus our attention on the k-center and k-median problems. Both of these problems are NP-hard but have known efficient approximation algorithms. The state of the art approaches give a 2-approximation for k-center (Gonzalez, 1985) and a  (1 +√3 + ϵ)-approximation for k-median (Li & Svensson, 2013).

Unlike clustering, the exploration of fairness in machine learning is relatively nascent. There are two broad lines of work. The first is in codifying what it means for an algorithm to be fair. See for example the work on statistical parity (Luong et al. , 2011; Kamishima et al. , 2012), disparate impact (Feldman et al. , 2015), and individual fairness (Dwork et al. , 2012). More recent work by Corbett-Davies et al. (2017) and Kleinberg et al. (2017b) also shows that some of the desired properties of fairness may be incompatible with each other.

A second line of work takes a specific notion of fairness and looks for algorithms that achieve fair outcomes. Here the focus has largely been on supervised learning (Luong et al. , 2011; Hardt et al. , 2016) and online (Joseph et al. , 2016) learning. The direction that is most similar to our work is that of learning intermediate representations that are guaranteed to be fair, see for example the work by Zemel et al. (2013) and Kamishima et al. (2012). However, unlike their work, we give strong guarantees on the relationship between the quality of the fairlet representation, and the quality of any fair clustering solution.

In this paper we use the notion of fairness known as disparate impact and introduced by Feldman et al. (2015). This notion is also closely related to the p%-rule as a measure for fairness. The p%-rule is a generalization of the 80%-rule advocated by US Equal Employment Opportunity Commission (Biddle, 2006) and was used in a recent paper on mechanism for fair classification (Zafar et al. , 2017). In particular our paper addresses an open question of Zafar et al. (2017) presenting a framework to solve an unsupervised learning task respecting the p%-rule.

Let X be a set of points in a metric space equipped with a distance function  d : X2 → R≥0. For an integer k, let [k] denote the set {1, . . . , k}.

We first recall standard concepts in clustering. A k-clustering C is a partition of X into k disjoint subsets, C1, . . . , Ck, called clusters. We can evaluate the quality of a clustering C with different objective functions. In the k-center problem, the goal is to minimize


and in the k-median problem, the goal is to minimize


A clustering C can be equivalently described via an assignment function  α : X → [k]. The points in cluster Ciare simply the pre-image of  i under α, i.e., Ci = {x ∈ X | α(x) = i}.

Throughout this paper we assume that each point in X is colored either red or blue; let  χ : X →{RED, BLUE}denote the color of a point. For a subset  Y ⊆ X and for c ∈ {RED, BLUE}, let c(Y ) = {x ∈X | χ(x) = c} and let #c(Y ) = |c(Y )|.

We first define a natural notion of balance.


The balance of a clustering C is defined as:


A subset with an equal number of red and blue points has balance 1 (perfectly balanced) and a monochromatic subset has balance 0 (fully unbalanced). To gain more intuition about the notion of balance, we investigate some basic properties that follow from its definition.

Lemma 2 (Combination). Let  Y, Y ′ ⊆ Xbe disjoint. If C is a clustering of  Y and C′ is a clustering of  Y ′,then  balance(C ∪ C′) = min(balance(C), balance(C′)).

It is easy to see that for any clustering  C of X, we have balance(C) ≤ balance(X). In particular, if X is not perfectly balanced, then no clustering of X can be perfectly balanced. We next show an interesting converse, relating the balance of X to the balance of a well-chosen clustering.

Lemma 3. Let balance(X) = b/r for some integers  1 ≤ b ≤ r such that gcd(b, r) = 1. Then there exists a clustering  Y = {Y1, . . . , Ym} of Xsuch that (i)  |Yj| ≤ b + r for each Yj ∈ Y, i.e., each cluster is small, and (ii) balance(Y) = b/r = balance(X).

Proof. Without loss of generality, let  B = #BLUE(X) ≤ #RED(X) = R. By assumption, B/R = b/r. We construct the clustering Y iteratively as follows.

If  (R − B) ≥ (r − b), then we remove r red points and b blue points from the current set to form a cluster Y . By construction |Y | = b + r and balance(Y ) = b/r. Furthermore the leftover set has balance (B − b)/(R − r) ≥ b/rand we iterate on this leftover set.

If  (R − B) < (r − b), then we remove  (R − B) + bred points and b blue points from the current set to form Y . Note that  |Y | ≤ b + r and that balance(Y ) = b/(R − B + b) ≥ b/r.

Finally note that when the remaining points are such that the red and the blue points are in a one-to-one correspondence, we can pair them up into perfectly balanced clusters of size 2.

Fairness and fairlets. Balance encapsulates a specific notion of fairness, where a clustering with a monochromatic cluster (i.e., fully unbalanced) is considered unfair. We call the clustering Y as described in Lemma 3 a (b, r)-fairlet decomposition of X and call each cluster  Y ∈ Y a fairlet.

Equipped with the notion of balance, we now revisit the clustering objectives defined earlier. The objectives do not consider the color of the points, so they can lead to solutions with monochromatic clusters. We now extend them to incorporate fairness.

Definition 4 ((t, k)-fair clustering problems). In the (t, k)-fair center (resp., (t, k)-fair median) problem, the goal is to partition  X into C such that |C| = k, balance(C) ≥ t, and φ(X, C) (resp. ψ(X, C)) is minimized.

Traditional formulations of k-center and k-median eschew the notion of an assignment function. Instead it is implicit through a set  {c1, . . . , ck}of centers, where each point assigned to its nearest center, i.e., α(x) = arg mini∈[1,k] d(x, ci).Without fairness as an issue, they are equivalent formulations; however, with fairness, we need an explicit assignment function (see Figure 1).

At first glance, the fair version of a clustering problem appears harder than its vanilla counterpart. In this section we prove a reduction from the former to the latter. We do this by first clustering the original points into small clusters preserving the balance, and then applying vanilla clustering on these smaller clusters instead of on the original points.

As noted earlier, there are different ways to partition the input to obtain a fairlet decomposition. We will show next that the choice of the partition directly impacts the approximation guarantees of the final clustering algorithm.

Before proving our reduction we need to introduce some additional notation. Let  Y = {Y1, . . . , Ym} bea fairlet decomposition. For each cluster  Yj, we designate an arbitrary point  yj ∈ Yj as its center. Then fora point  x, we let β : X → [1, m]denote the index of the fairlet to which it is mapped. We are now ready to define the cost of a fairlet decomposition

Definition 5 (Fairlet decomposition cost). For a fairlet decomposition, we define its k-median cost as �x∈X d(x, β(x)), and its k-center cost as  maxx∈X d(x, β(x)).We say that a (b, r)-fairlet decomposition is optimal if it has minimum cost among all (b, r)-fairlet decompositions.

Since (X, d) is a metric, we have from the triangle inequality that for any other point  c ∈ X,


Now suppose that we aim to obtain a (t, k)-fair clustering of the original points X. (As we observed earlier, necessarily  t ≤ balance(X).) To solve the problem we can cluster instead the centers of each fairlet, i.e., the set  {y1, . . . , ym} = Y , into kclusters. In this way we obtain a set of centers  {c1, . . . , ck} and anassignment function  αY : Y → [k].

We can then define the overall assignment function as  α(x) = αY (yβ(x))and denote the clustering induced by  α as Cα. From the definition of Y and the property of fairlets and balance, we get that balance(Cα) = t. We now need to bound its cost. Let ˜Ybe a multiset, where each  yi appears |Yi| numberof times.

Lemma 6. ψ(X, Cα) = ψ(X, Y) + ψ( ˜Y , Cα) and φ(X, Cα) = φ(X, Y) + φ( ˜Y , Cα).

Proof. We prove the result for the k-median setting; the k-center version is similar. Let  Cα = {C1, . . . , Ck},with corresponding centers  {c1, . . . , ck}. Using the definition of the k-median objective and the triangle inequality we get,


Therefore in both cases we can reduce the fair clustering problem to the problem of finding a good fairlet decomposition and then solving the vanilla clustering problem on the centers of the fairlets. We refer to  ψ(X, Y) and φ(X, Y) as the k-median and k-center costs of the fairlet decomposition.

In the previous section we presented a reduction from the fair clustering problem to the regular counterpart. In this section we use it to design efficient algorithms for fair clustering.

We first focus on the k-center objective and show in Section 4.3 how to adapt the reasoning to solve the k-median objective. We begin with the most natural case in which we require the clusters to be perfectly balanced, and give efficient algorithms for the (1, k)-fair center problem. Then we analyze the more challenging (t, k)-fair center problem for  t < 1. Let B = BLUE(X), R = RED(X).

4.1 Fair k-center warmup: (1, 1)-fairlets

Suppose balance(X) = 1, i.e., (|R| = |B|) and we wish to find a perfectly balanced clustering. We now show how we can obtain it using a good (1, 1)-fairlet decomposition.

Lemma 7. An optimal (1, 1)-fairlet decomposition for k-center can be found in polynomial time.

Proof. To find the best decomposition, we first relate this question to a graph covering problem. Consider a bipartite graph  G = (B ∪ R, E)where we create an edge  E = (bi, rj)with weight  wij = d(ri, bj)between any bichromatic pair of nodes. In this case a decomposition into fairlets corresponds to some perfect matching in the graph. Each edge in the matching represents a fairlet,  Yi. Let Y = {Yi}be the set of edges in the matching.

Observe that the k-center cost  φ(X, Y)is exactly the cost of the maximum weight edge in the matching, therefore our goal is to find a perfect matching that minimizes the weight of the maximum edge. This can be done by defining a threshold graph  Gτthat has the same nodes as G but only those edges of weight at most  τ. We then look for the minimum  τwhere the corresponding graph has a perfect matching, which can be done by (binary) searching through the  O(n2) values.

Finally, for each fairlet (edge)  Yiwe can arbitrarily set one of the two nodes as the center,  yi.

Since any fair solution to the clustering problem induces a set of minimal fairlets (as described in Lemma 3), the cost of the fairlet decomposition found is at most the cost of the clustering solution.

Lemma 8. Let Y be the partition found above, and let  φ∗t be the cost of the optimal (t, k)-fair center clustering. Then  φ(X, Y) ≤ φ∗t .

This, combined with the fact that the best approximation algorithm for k-center yields a 2-approximation (Gonzalez, 1985), gives us the following.

Theorem 9. The algorithm that first finds fairlets and then clusters them is a 3-approximation for the (1, k)-fair center problem.

4.2 Fair k-center: (1, t)-fairlets

Now, suppose that instead we look for a clustering with balance  t ⪇ 1. In this section we assume  t = 1/t′for some integer  t′ > 1. We show how to extend the intuition in the matching construction above to find approximately optimal  (1, t′)-fairlet decompositions for integral  t′ > 1.

In this case, we transform the problem into a minimum cost flow (MCF) problem.1 Let  τ > 0 bea parameter of the algorithm. Given the points B, R, and an integer  t′, we construct a directed graph Hτ = (V, E). Its node set V is composed of two special nodes  β and ρ, all of the nodes in  B ∪ R, and t′additional copies for each node  v ∈ B ∪ R. More formally,



Figure 2: The construction of the MCF instance for the bipartite graph for  t′ = 2. Note that the only nodes with positive demands or supplies are  β, ρ, b1, b2, b3, r1, and r2and all the dotted edges have cost 0.

The directed edges of  Hτare as follows:


(ii) A  (β, bi)edge for each  bi ∈ B, and an (ri, ρ)edge for each  ri ∈ R. All of these edges have cost 0 and capacity  t′ − 1.

(iii) For each  bi ∈ Band for each  j ∈ [t′], a (bi, bji)edge, and for each  ri ∈ Rand for each  j ∈ [t′], an(ri, rji )edge. All of these edges have cost 0 and capacity 1.

(iv) Finally, for each  bi ∈ B, rj ∈ Rand for each  1 ≤ k, ℓ ≤ t, a (bki , rℓj)edge with capacity 1. The cost of this edge is  1 if d(bi, rj) ≤ τ and ∞otherwise.

To finish the description of this MCF instance, we have now specify supply and demand at every node. Each node in B has a supply of 1, each node in R has a demand of  1, βhas a supply of  |R|, and ρ has ademand of |B|. Every other node has zero supply and demand. In Figure 2 we show an example of this construction for a small graph.

To finish the description of this MCF instance, we specify the supply and demand at every node. Each node in B has a supply of 1, each node in R has a demand of  1, βhas a supply of  |R|, and ρhas a demand of |B|. Every other node has zero supply and demand. In Figure 2 we show an example of this construction for a small graph.

The MCF problem can be solved in polynomial time and since all of the demands and capacities are integral, there exists an optimal solution that sends integral flow on each edge. In our case, the solution is a set of edges of  Hτthat have non-zero flow, and the total flow on the  (β, ρ) edge.

In the rest of this section we assume for simplicity that any two distinct elements of the metric are at a positive distance apart and we show that starting from a solution to the described MCF instance we can build a low cost  (1, t′)-fairlet decomposition. We start by showing that every  (1, t′)-fairlet decomposition can be used to construct a feasible solution for the MCF instance and then prove that an optimal solution for the MCF instance can be used to obtain a  (1, t′)-fairlet decomposition.

Lemma 10. Let  Y be a (1, t′)-fairlet decomposition of cost  C for the (1/t′, k)-fair center problem. Then it is possible to construct a feasible solution of cost 2C to the MCF instance.

Proof. We begin by building a feasible solution and then bound its cost. Consider each fairlet in the  (1, t′)-fairlet decomposition.

Suppose the fairlet contains 1 red node and c blue nodes, with  c ≤ t′, i.e., the fairlet is of the form {r1, b1, . . . , bc}. For any such fairlet we send a unit of flow form each node  bi to b1i , for i ∈ [c]and a unit of flow from nodes  b11, . . . , b1c to nodes r11, . . . , rc1. Furthermore we send a unit of flow from each  r11, . . . , rc1 tor1 and c − 1units of flow from  r1 to ρ. Note that in this way we saturate the demands of all nodes in this fairlet.

Similarly, if the fairlet contains c red nodes and 1 blue node, with  c ≤ t′, i.e., the fairlet is of the form {r1, . . . , rc, b1}. For any such fairlet, we send  c − 1units of flow from  β to b1. Then we send a unit of flow from each  b1 to each b11, . . . , bc1 and a unit of flow from nodes  b11, . . . , bc1 to nodes r11, . . . , r1c. Furthermore we send a unit of flow from each  r11, . . . , r1c to the nodes  r1, . . . , rc. Note that also in this case we saturate all the request of nodes in this fairlet.

Since every node  v ∈ B ∪ Ris contained in a fairlet, all of the demands of these nodes are satisfied. Hence, the only nodes that can have still unsatisfied demand are  β and ρ, but we can use the direct edge (β, ρ)to route the excess demand, since the total demand is equal to the total supply. In this way we obtain a feasible solution for the MCF instance starting from a  (1, t′)-fairlet decomposition.

To bound the cost of the solution note that the only edges with positive cost in the constructed solution are the edges between nodes  bji and rℓk. Furthermore an edge is part of the solution only if the nodes  bi andrkare contained in the same fairlet F. Given that the k-center cost for the fairlet decomposition is C, the cost of the edges between nodes in F in the constructed feasible solution for the MCF instance is at most 2 times this distance. The claim follows.

Now we show that given an optimal solution for the MCF instance of cost C, we can construct a  (1, t′)-fairlet decomposition of cost no bigger than C.

Lemma 11. Let Y be an optimal solution of cost C to the MCF instance. Then it is possible to construct a (1, t′)-fairlet decomposition for  (1/t′, k)-fair center problem of cost at most C.

Proof. First we show that from an optimal solution for the MCF instance, it is possible to construct a  (1, t′)-fairlet decomposition. Then we bound the cost of the decomposition.

Let  E⋆ be the subset of edges of  {{bi, rj} | bi ∈ B and rj ∈ R} such that bi and rjare connected by edges used in the feasible solution for the MCF. Denote by  GE⋆the graph induced by  E⋆. Note that by construction the degree of each node  bi or rj in GE⋆is at most  t′ and at least 1.

We claim that  GE⋆is a collection of stars, each having a number of leaves in  {1, . . . , t′}. In fact, suppose a component of  E⋆ is not a star, that is it contains x red nodes, and y blue nodes, with  x, y ≥ 2. In this case there are two nodes  r, b ∈ E⋆ each with degree at least 2 that are connected to each other. We can safely remove the edge (r, b) while still guaranteeing that every red and every blue node in the component has at least one neighbor of the opposite color. Removing this edge will decrease the cost of the solution which contradicts the optimality of  E⋆. (Indeed, the flow that passed through the removed edge can be rerouted through the  (β, ρ) edge.)

Therefore, we can define the  (1, t′)-fairlet decomposition as the set of connected components in  GE⋆.To finish the lemma we only need to bound the cost of the decomposition.

For each fairlet, we designate as the center the node of the highest degree. It is easy to see that this solution has cost bounded by C.

Combining Lemma 10 and Lemma 11 yields the following.

Lemma 12. By reducing the  (1, t′)-fairlet decomposition problem to an MCF problem, it is possible to compute a 2-approximation for the optimal  (1, t′)-fairlet decomposition for the  (1/t′, k)-fair center problem.

Note that the cost of a  (1, t′)-fairlet decomposition is necessarily smaller than the cost of a  (1/t′, k)-fairclustering. Our main theorem follows.

Theorem 13. The algorithm that first finds fairlets and then clusters them is a 4-approximation for the (1/t′, k)-fair center problem for any positive integer  t′.

4.3 Fair k-median

The results in the previous section can be modified to yield results for the (t, k)-fair median problem with minor changes that we describe below.

For the perfectly balanced case, as before, we look for a perfect matching on the bichromatic graph. Unlike, the k-center case, our goal is to find a perfect matching of minimum total cost, since that exactly represents the cost of the fairlet decomposition. Since the best known approximation for k-median is 1 + √3 + ϵ(Li & Svensson, 2013), we have:

Theorem 14. The algorithm that first finds fairlets and then clusters them is a  (2 +√3 + ϵ)-approximation for the (1, k)-fair median problem.

To find  (1, t′)-fairlet decompositions for integral  t′ > 1, we again resort to MCF and create an instance as in the k-center case, but for each  bi ∈ B, rj ∈ R, and for each  1 ≤ k, ℓ ≤ t, we set the cost of the edge (bki , rℓj) to d(bi, rj).

Theorem 15. The algorithm that first finds fairlets and then clusters them is a  (t′+1+√3+ϵ)-approximation for the  (1/t′, k)-fair median problem for any positive integer  t′.

4.4 Hardness

We complement our algorithmic results with a discussion of computational hardness for fair clustering. We show that the question of finding a good fairlet decomposition is itself computationally hard. Thus, ensuring fairness causes hardness, regardless of the underlying clustering objective.

Theorem 16. For each fixed  t′ ≥ 3, finding an optimal  (1, t′)-fairlet decomposition is NP-hard. Also, finding the minimum cost  (1/t′, k)-fair median clustering is NP-hard.

Proof. We reduce from the problem of partitioning the node set of a graph G = (V, E) into induced subgraphs of order  t′ each having eccentricity 1. Equivalently, this question asks whether V can be partitioned into pairwise disjoint subsets  S1, . . . , S|V |/t′, so that |Si| = t′ and G|St′is a star with  t′ − 1leaves, for each  i = 1, . . . , |V |/t′. This problem was shown to be NP-hard by Kirkpatrick & Hell (1978), see also van Bevern et al. (2016).

Assume that |V | is divisible by  t′. We create one red element for each node in  V , and |V |/t′ blue ele-ments. The distance between any two red elements will be 1 if the corresponding nodes in V are connected by an edge, and 2 otherwise. The distance between any two blue elements will be 2. Finally, the distance between any red element and any blue element will be 2. For the fairlet decomposition problem, we ask whether this instance admits a  (1, t′)-fairlet decomposition having total cost upper bounded by�1 + 1t′�·|V |.For the  (1/t′, k)-fair median problem, we ask whether the instance admits a k-clustering, with  k = |V |/t′,having median cost at most�1 + 1t′�· |V |.


Figure 3: Empirical performance of the classical and fair clustering median and center algorithms on the three datasets. The cost of each solution is on left axis, and its balance on the right axis.

Observe that the distance function we defined is trivially a metric, since all of its values are in {1, 2}.

Suppose that G can be partitioned into induced subgraphs of order  t′, with node sets  S1, . . . , S|V |/t′,each with eccentricity  1. For each i = 1, . . . , |V |/t′, we create one cluster (or, one fairlet) with the red elements corresponding to the nodes in  Si, and the ith blue element. Then, each cluster will contain a red element at distance 1 from each of the other  t′−1red elements, and at distance 2 from the only blue element. The cost of each cluster will then be at most  t′ + 1. The total cost is then at most�1 + 1t′�· |V |.

On the other hand, observe that since the number of blue elements is  |V |/t′ and the number of red elements is |V | any feasible solution has to create  k = |V |/t′ clusters (or fairlets) each containing exactly 1 blue element and  t′ red elements. Now, the median cost of a cluster (or fairlet) is  t′ + 1if the nodes corresponding to its  t′ red points induce a star, and it is at least  t′ + 2otherwise. It follows that, if G cannot be partitioned into induced subgraphs of order  t′ with eccentricity 1, the total cost (of either problems) will be at least�1 + 1t′�· |V | + 1.

In this section we illustrate our algorithm by performing experiments on real data. The goal of our experiments is two-fold: first, we show that traditional algorithms for k-center and k-median tend to produce unfair clusters; second, we show that by using our algorithms one can obtain clusters that respect the fairness guarantees. We show that in the latter case, the cost of the solution tends to converge to the cost of the fairlet decomposition, which serves as a lower bound on the cost of the optimal solution.

Datasets. We consider 3 datasets from the UCI repository Lichman (2013) for experimentation.

Diabetes. This dataset2 represents the outcomes of patients pertaining to diabetes. We chose numeric attributes such as age, time in hospital, to represent points in the Euclidean space and gender as the sensitive dimension, i.e., we aim to balance gender. We subsampled the dataset to 1000 records.

Bank. This dataset3 contains one record for each phone call in a marketing campaign ran by a Portuguese banking institution (Moro et al. , 2014). Each record contains information about the client that was contacted by the institution. We chose numeric attributes such as age, balance, and duration to represent points in the Euclidean space, we aim to cluster to balance married and not married clients. We subsampled the dataset to 1000 records.

Census. This dataset4 contains the census records extracted from the 1994 US census (Kohavi, 1996). Each record contains information about individuals including education, occupation, hours worked per week, etc. We chose numeric attributes such as age, fnlwgt, education-num, capital-gain and hours-per-week to represents points in the Euclidean space and we aim to cluster the dataset so to balance gender. We subsampled the dataset to 600 records. Algorithms. We implement the flow-based fairlet decomposition algorithm as described in Section 4. To solve the k-center problem we augment it with the greedy furthest point algorithm due to Gonzalez (1985), which is known to obtain a 2-approximation. To solve the k-median problem we use the single swap algorithm due to Arya et al. (2004), which also gets a 5-approximation in the worst case, but performs much better in practice (Kanungo et al. , 2002). Results. Figure 3 shows the results for k-center for the three datasets in the top row and the k-median objective in the bottom row. In all of the cases, we run with  t′ = 2, that is we aim for balance of at least 0.5 in each cluster.

Observe that the balance of the solutions produced by the classical algorithms is very low, and in four out of the six cases, the balance is 0 for larger values of k, meaning that the optimal solution has monochromatic clusters. Moreover, this is not an isolated incident, for instance the k-median instance of the Bank dataset has three monochromatic clusters starting at k = 12. Finally, left unchecked, the balance in all datasets keeps decreasing as the clustering becomes more discriminative, with increased k.

On the other hand the fair clustering solutions maintain a balanced solution even as k increases. Not surprisingly, the balance comes with a corresponding increase in cost, and the fair solutions are costlier than their unfair counterparts. In each plot we also show the cost of the fairlet decomposition, which represents the limit of the cost of the fair clustering; in all of the scenarios the overall cost of the clustering converges to the cost of the fairlet decomposition.

In this work we initiate the study of fair clustering algorithms. Our main result is a reduction of fair clustering to classical clustering via the notion of fairlets. We gave efficient approximation algorithms for finding fairlet decompositions, and proved lower bounds showing that fairness can introduce a computational bottleneck. An immediate future direction is to tighten the gap between lower and upper bounds by improving the approximation ratio of the decomposition algorithms, or giving stronger hardness results. A different avenue is to extend these results to situations where the protected class is not binary, but can take on multiple values. Here there are multiple challenges including defining an appropriate version of fairness.

Flavio Chierichetti was supported in part by the ERC Starting Grant DMAP 680153, by a Google Focused Research Award, and by the SIR Grant RBSI14Q743.

Aggarwal, Charu C., & Reddy, Chandan K. 2013. Data Clustering: Algorithms and Applications. 1st edn. Chapman & Hall/CRC.

Arya, Vijay, Garg, Naveen, Khandekar, Rohit, Meyerson, Adam, Munagala, Kamesh, & Pandit, Vinayaka. 2004. Local search heuristics for k-median and facility location problems. SIAM J. Comput., 33(3), 544–562.

Biddle, Dan. 2006. Adverse Impact and Test Validation: A Practitioner’G guide to Valid and Defensible Employment Testing. Gower Publishing, Ltd.

Corbett-Davies, Sam, Pierson, Emma, Feller, Avi, Goel, Sharad, & Huq, Aziz. 2017. Algorithmic Decision Making and the Cost of Fairness. Pages 797–806 of: Proceedings of the 23rd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. KDD ’17. New York, NY, USA: ACM.

Dwork, Cynthia, Hardt, Moritz, Pitassi, Toniann, Reingold, Omer, & Zemel, Richard. 2012. Fairness through awareness. Pages 214–226 of: ITCS.

Feldman, Michael, Friedler, Sorelle A., Moeller, John, Scheidegger, Carlos, & Venkatasubramanian, Suresh. 2015. Certifying and removing disparate impact. Pages 259–268 of: KDD.

Gonzalez, T. 1985. Clustering to minimize the maximum intercluster distance. TCS, 38, 293–306.

Hardt, Moritz, Price, Eric, & Srebro, Nati. 2016. Equality of opportunity in supervised learning. Pages 3315–3323 of: NIPS.

Joseph, Matthew, Kearns, Michael, Morgenstern, Jamie H., & Roth, Aaron. 2016. Fairness in learning: Classic and contextual bandits. Pages 325–333 of: NIPS.

Kamishima, Toshihiro, Akaho, Shotaro, Asoh, Hideki, & Sakuma, Jun. 2012. Fairness-aware classifier with prejudice remover regularizer. Pages 35–50 of: ECML/PKDD.

Kanungo, Tapas, Mount, David M., Netanyahu, Nathan S., Piatko, Christine D., Silverman, Ruth, & Wu, Angela Y. 2002. An efficient k-means clustering algorithm: Analysis and implementation. PAMI, 24(7), 881–892.

Kirkpatrick, David G., & Hell, Pavol. 1978. On the completeness of a generalized matching problem. Pages 240–245 of: STOC.

Kleinberg, Jon, Lakkaraju, Himabindu, Leskovec, Jure, Ludwig, Jens, & Mullainathan, Sendhil. 2017a. Human decisions and machine predictions. Working Paper 23180. NBER.

Kleinberg, Jon M., Mullainathan, Sendhil, & Raghavan, Manish. 2017b. Inherent trade-offs in the fair determination of risk scores. In: ITCS.

Kohavi, Ron. 1996. Scaling up the accuracy of naive-Bayes classifiers: A decision-tree hybrid. Pages 202–207 of: KDD.

Li, Shi, & Svensson, Ola. 2013. Approximating k-median via pseudo-approximation. Pages 901–910 of: STOC.

Lichman, M. 2013. UCI Machine Learning Repository.

Luong, Binh Thanh, Ruggieri, Salvatore, & Turini, Franco. 2011. k-NN as an implementation of situation testing for discrimination discovery and prevention. Pages 502–510 of: KDD.

Moro, S´ergio, Cortez, Paulo, & Rita, Paulo. 2014. A data-driven approach to predict the success of bank telemarketing. Decision Support Systems, 62, 22–31.

van Bevern, Ren, Bredereck, Robert, Bulteau, Laurent, Chen, Jiehua, Froese, Vincent, Niedermeier, Rolf, & Woeginger, Gerhard J. 2016. Partitioning perfect graphs into stars. Journal of Graph Theory, 85(2), 297–335.

Xu, Rui, & Wunsch, Don. 2009. Clustering. Wiley-IEEE Press.

Zafar, Muhammad Bilal, Valera, Isabel, Gomez-Rodriguez, Manuel, & Gummadi, Krishna P. 2017. Fairness constraints: Mechanisms for fair classification. Pages 259–268 of: AISTATS.

Zemel, Richard S., Wu, Yu, Swersky, Kevin, Pitassi, Toniann, & Dwork, Cynthia. 2013. Learning fair representations. Pages 325–333 of: ICML.

Designed for Accessibility and to further Open Science