---
title: Individually Fair Clustering
url: https://www.emergentmind.com/topics/individually-fair-clustering
type: topic
---

# Individually Fair Clustering

Individually fair clustering studies clustering under constraints intended to formalize the principle that similar individuals should be treated similarly. In the clustering literature, this principle has not converged to a single canonical definition. Instead, several distinct formulations have been developed: nearest-neighbor coverage constraints based on a point-specific fair radius; feature-based co-clustering constraints defined in a fairness feature space; service-parity constraints in which similar points should receive comparable center distances; average-dissimilarity conditions requiring each point to be closer, on average, to its own cluster than to any other; and distributional formulations in which similar individuals receive statistically similar distributions over centers [2002.06742] [2109.04554] [2106.05423] [2006.04960] [2006.12589]. The area therefore combines approximation algorithms, LP rounding, local search, dynamic programming, and specialized extensions for outliers, graph data, and representation learning [2106.12150] [2402.06730] [2412.10923] [2510.06130] [2402.10756] [2002.10312].

## 1. Core formulations

A central formulation defines, for each point \(x\) in a set \(P\) of size \(n\), a fair radius
\[
r(x) = \min \{r : |B(x, r)| \geq n/k\},
\]
or equivalently \(\delta(x)\) as the radius of the smallest ball around \(x\) containing at least \(n/k\) points. A clustering with center set \(S\) is individually fair when every point has a center within that radius, or within a relaxed radius \(\alpha r(x)\):
\[
d(x,S) \leq \alpha \cdot r(x).
\]
This fair-radius model underlies a large part of the algorithmic literature on individually fair \(k\)-median, \(k\)-means, and \(k\)-center [2002.06742] [2106.14043] [2106.12150] [2402.06730].

A different line of work defines fairness through a separate feature space. In "Feature-based Individual Fairness in \(k\)-Clustering" [2109.04554], each point \(v\) has a feature vector \(X_v \in \mathbb{R}^q\), similarity is defined by
\[
s(X_i, X_j) = e^{-d'(X_i, X_j)},
\]
two points are \(\gamma\)-similar if \(s(X_i,X_j) > \gamma\), and \(\Gamma(v)\) denotes the points \(\gamma\)-similar to \(v\). A clustering with assignment \(\phi(v)\) is \((m_v)_{v\in V}\)-individually fair if
\[
|\{u \in \Gamma(v) : \phi(u) = \phi(v)\}| \geq m_v.
\]
Here fairness is not proximity to a center, but co-membership with sufficiently many similar points.

In "A New Notion of Individually Fair Clustering: \(\alpha\)-Equitable \(k\)-Center" [2106.05423], fairness is defined by parity of service quality relative to a similarity set \(\mathcal{S}_j\). The per-point constraint is
\[
d(j,\phi(j)) \leq \alpha \cdot \min_{j' \in \mathcal{S}_j} d(j',\phi(j')),
\]
and the aggregate variant is
\[
d(j, \phi(j)) \leq \alpha \cdot \frac{1}{|\mathcal{S}_j|} \sum_{j' \in \mathcal{S}_j} d(j',\phi(j')).
\]
This formulation compares the service a point receives to the service received by points it considers similar.

A further formulation, introduced in "A Notion of Individual Fairness for Clustering" [2006.04960], requires that each point be closer on average to points in its own cluster than to points in any other cluster. If \(C(x)\) is the cluster containing \(x\), then \(x\) is treated individually fair if either \(C(x)=\{x\}\), or
\[
\frac{1}{|C(x)|-1}\sum_{y\in C(x)\setminus\{x\}} d(x,y)
\leq
\frac{1}{|C_i|}\sum_{y\in C_i} d(x,y)
\quad
\forall i\in[k],\ C_i\neq C(x).
\]

A probabilistic generalization appears in "Distributional Individual Fairness in Clustering" [2006.12589]. Each individual \(j\) is assigned a distribution \(\mu_j\) over at most \(k\) centers, and fairness is imposed through \(f\)-divergence:
\[
D_f(\mu_{j_1}\,\|\,\mu_{j_2}) \leq d(j_1, j_2)\qquad \forall j_1, j_2 \in V.
\]
The clustering objective is then the expected \(\ell_p\)-norm cost,
\[
\left( \sum_{j\in V} \mathbb{E}_{c\sim \mu_j}[d(j,c)^p] \right)^{1/p}.
\]

| Formulation | Core fairness requirement | Representative source |
|---|---|---|
| Fair-radius coverage | \(d(x,S) \leq \alpha r(x)\) | [2002.06742] |
| Feature-based co-clustering | \(|\{u\in\Gamma(v):\phi(u)=\phi(v)\}|\ge m_v\) | [2109.04554] |
| \(\alpha\)-equitable service parity | Service within factor \(\alpha\) of similar points | [2106.05423] |
| Average-dissimilarity fairness | Own-cluster average distance no worse than any other cluster | [2006.04960] |
| Distributional fairness | \(D_f(\mu_{j_1}\|\mu_{j_2}) \le d(j_1,j_2)\) | [2006.12589] |

## 2. Feasibility, hardness, and price of fairness

Several formulations of individually fair clustering are computationally hard even at the level of feasibility. For the average-dissimilarity notion, deciding whether a data set admits an individually fair \(k\)-clustering is NP-hard, even when \(k=2\) and the dissimilarity is a metric [2006.04960]. For the feature-based model, deciding whether there exists a clustering with cost less than a trivially fair clustering is NP-hard even when only two facilities are allowed; distinguishing between instances with zero and non-zero optimal cost is NP-complete even for \(k=2\); there is no constant-factor approximation for general instances; and even under triangle inequality there is no fully polynomial time approximation scheme unless \(P=NP\) [2109.04554].

The \(\alpha\)-equitable \(k\)-center model yields a structural feasibility threshold. For both EqCenter-PP and EqCenter-AG, there exist instances without any feasible solution when \(\alpha<2\), whereas every instance admits a feasible solution for \(\alpha\ge 2\). In that sense, \(\alpha\ge 2\) is the minimally possible threshold for universal feasibility in that model [2106.05423].

The literature also studies the cost incurred by imposing fairness. In the \(\alpha\)-equitable setting, the Price of Fairness can be arbitrarily large in the worst case, although bounded-Price-of-Fairness guarantees are obtained in some structured cases, including settings with up to \(2k\) centers and, for AG, equal similarity radii [2106.05423]. In the fair-radius literature, standard \(k\)-median or \(k\)-means solutions can be arbitrarily unfair according to the individual fairness criterion, which motivates fairness-specific algorithm design rather than post hoc auditing of unconstrained clusterings [2002.06742].

These hardness and feasibility results explain why the field is dominated by bicriteria guarantees, relaxations of the fairness radius, randomized procedures, and restricted settings such as line metrics or fixed centers.

## 3. Algorithmic paradigms for fair-radius and service-based models

The fair-radius model has been studied through local search, LP rounding, and reductions to constrained facility-location variants. "Individual Fairness for \(k\)-Clustering" [2002.06742] augments classical local search for \(k\)-median and \(k\)-means by requiring every candidate solution to cover a family of critical balls. The algorithm computes disjoint critical balls \(B(c,\alpha r(c))\), initializes centers so that every critical ball contains a center, and performs local-search swaps of up to 4 centers while preserving this coverage invariant. The resulting guarantees are an \((O(1),O(1))\)-bicriteria approximation for \(k\)-median and \(k\)-means, an \((O(\log n),O(1))\)-approximation for \(k\)-center, and an \((O(p),O(1))\)-approximation for general \(\ell_p\) cost; the number of iterations is \(O(\log n/\varepsilon)\), and the overall runtime is polynomial, specifically \(O(k^5 n^4)\) when swaps of size at most 4 are allowed [2002.06742].

"Improved Approximation Algorithms for Individually Fair Clustering" [2106.14043] replaces the earlier \(p^{O(p)}\)-style approximation factors by stronger bicriteria bounds. For any \(\varepsilon>0\), it gives a \((16^p+\varepsilon,3)\)-bicriteria approximation for \(\ell_p\)-clustering, a \((7.081+\varepsilon,3)\)-bicriteria approximation for \(k\)-median, and a \((3+\varepsilon,3)\)-bicriteria approximation for \(k\)-center. The algorithm constructs disjoint critical regions, reduces the problem to facility location with \(\ell_p\)-cost under a matroid constraint, and uses an extended LP-based framework inspired by Charikar et al. and Swamy. The same work also describes the reduction from individually fair clustering to a clustering problem with a group fairness requirement proposed by Kleindessner et al., which is essentially the median matroid problem [2106.14043].

A second LP-based direction appears in "Better Algorithms for Individually Fair \(k\)-Clustering" [2106.12150]. Its LP has variables \(y_u\) for opening centers and \(x_{vu}\) for assignments, with the fairness condition enforced explicitly through
\[
x_{vu}=0 \quad \forall v,u:\ d(v,u)>r(v).
\]
The Fair-Round procedure first solves the LP, then filters points using radii
\[
R(v) := \min\{r(v), (2C_v)^{1/p}\},
\]
where \(C_v := \sum_{u \in X} d(v,u)^p x^*_{vu}\), and finally rounds the resulting representative structure. The paper states an \((8,2^{1+2/p})\)-approximation algorithm, with fairness factor \(8\) and cost factor \(2^{1+2/p}\); the special cases are \((8,8)\) for \(k\)-median and \((8,4)\) for \(k\)-means. It also introduces a sparsification technique that reduces LP size while preserving bounded degradation in fairness and cost [2106.12150].

"A Scalable Algorithm for Individually Fair K-means Clustering" [2402.06730] emphasizes practical scale. It gives the first fast local-search algorithm running in \(\tilde O(nk^2)\) time with an \((O(1),6)\) bicriteria approximation for individually fair \(k\)-means. The method uses anchor points and anchor zones, modified greedy seeding, and constrained swaps that preserve coverage of each anchor zone. The design goal is explicitly to retain theoretical guarantees while scaling beyond the sizes handled by earlier local-search or LP-based baselines [2402.06730].

For the \(\alpha\)-equitable \(k\)-center objective, combinatorial algorithms are also available. Given optimal value \(R^*\) and \(R_m=\max_j \max_{j'\in \mathcal{S}_j} d(j,j')\), the main guarantee is
\[
\max_{j \in \mathcal{C}} d(j,\phi(j)) \leq 5 \cdot \max\{R^*,R_m\},
\]
which becomes a 10-approximation under the structured condition \(R_m \le 2R^*_{unf}\). The same work notes that assigning each point to its nearest center does not always respect the fairness constraints, and therefore develops a non-trivial iterative assignment procedure [2106.05423].

## 4. Feature-based and distributional individual fairness

The feature-based formulation separates the space used to define fairness from the space used to optimize clustering cost. In "Feature-based Individual Fairness in \(k\)-Clustering" [2109.04554], the algorithm LP-FAIR first selects initial centers using a standard \(k\)-clustering routine that ignores fairness, then solves a linear relaxation over fixed centers:
\[
\min_x \sum_{v,f_k} d(v,f_k)^p x_{v,f_k}
\]
subject to
\[
\sum_{u\in \Gamma(v)} x_{u,f_k} \ge m_v \cdot x_{v,f_k}\quad \forall v,f_k,
\qquad
\sum_{f_k} x_{v,f_k}=1,
\qquad
x_{v,f_k}\ge 0.
\]
The solution is rounded by assigning each point \(v\) to facility \(f_k\) with probability \(x^*_{v,f_k}\), repeating the procedure \(T=O(\log n)\) times and returning the lowest-cost solution. The main theorem states that for any \(\rho,\alpha>0\), if \(m_v=\Omega(k\log n/\alpha^2)\) for all \(v\), then with high probability the algorithm returns cost at most \((1+\rho)\cdot \mathrm{OPT}\), where \(\mathrm{OPT}\) is the LP optimum, and each point \(v\) has at least \(\frac{m_v}{k}(1-\alpha)\) similar points in its assigned cluster [2109.04554].

This model is notable because the fairness features may be arbitrary and may be disjoint from clustering features. The paper therefore departs from the nearest-neighbor fair-radius tradition: fairness is expressed by peer co-location rather than by center accessibility. The same framework can be composed with a black-box vanilla clustering algorithm to obtain an overall fair clustering algorithm when distances satisfy triangle inequality [2109.04554].

The distributional model in "Distributional Individual Fairness in Clustering" [2006.12589] replaces deterministic assignment by a distribution \(\mu_j\) over centers for each point \(j\). The fair assignment problem is expressed as a convex program:
\[
\min_{x\ge 0} \sum_{j\in V}\sum_{c\in C} x_{cj} d(c,j)^p
\]
subject to
\[
\sum_{c\in C} x_{cj}=1,
\qquad
D_f(\vec{x}_{j_1}\,\|\,\vec{x}_{j_2}) \le d(j_1,j_2),
\qquad
0\le x_{cj}\le 1.
\]
Its algorithm ALG-IF first runs a \(\rho\)-approximation for standard clustering to obtain centers \(C\), then solves the fair-assignment LP for these fixed centers. The theorem gives
\[
\text{ALG-IF cost} \le 3^{1-1/p}(\rho+2)\cdot \mathrm{OPT}_{k,p,f}.
\]
The same framework extends to combined group fairness and individual fairness within protected groups by adding constraints
\[
\sum_{j\in G_r} x_{cj} \le \alpha_r \sum_{j\in V} x_{cj},
\qquad
\sum_{j\in G_r} x_{cj} \ge \beta_r \sum_{j\in V} x_{cj},
\]
together with within-group individual fairness constraints. The paper also observes that enforcing individual fairness globally can force all points to be assigned uniform distributions, and for that reason advocates within-group enforcement in the combined setting [2006.12589].

## 5. Special settings and extensions

The average-dissimilarity formulation admits an exact algorithm on the real line. For one-dimensional Euclidean data, "A Notion of Individual Fairness for Clustering" [2006.04960] proves that for any \(k\) there exists a fair \(k\)-clustering with contiguous clusters, and gives a dynamic programming algorithm with runtime \(O(n^3k)\). Its recurrence is presented as
\[
T(i,j,l) = |j - t_l|^p + \min_{\substack{s: \text{fairness constraints}}} T(i-j,s,l-1),
\]
with fairness checks reduced to cluster boundaries in the contiguous setting [2006.04960].

Outlier-aware individually fair clustering has become a separate subarea. "Linear Programming based Approximation to Individually Fair \(k\)-Clustering with Outliers" [2412.10923] extends fair-radius \(k\)-means to allow a set \(Z\) of outliers with \(|Z|\le m\), excluded from both cost and fairness. The LP introduces variables \(x_{vu}\), \(y_u\), and \(z_v\):
\[
\min \sum_{u,v\in X} d(u,v)^p x_{vu}
\]
subject to
\[
\sum_{u\in X} y_u \le k,\quad
\sum_{v\in X} z_v \le m,\quad
\sum_u x_{vu} \ge 1-z_v,\quad
x_{vu}\le y_u,\quad
y_u \le 1-z_u,\quad
x_{vu}=0 \text{ if } d(u,v)>\alpha r(v).
\]
After solving the LP, the OutRound procedure thresholds outliers and the FairRound algorithm of Negahbani and Chakrabarty rounds the remaining fractional solution. The paper states
\[
LP_{\alpha=2}(x',y',z') \le 3\,LP_{\alpha=1}(x^*,y^*,z^*),
\]
followed by final guarantees \(FR(x',y') \le 12\,\mathrm{OPT}\) for \(k\)-means, \(FR(x',y') \le 24\,\mathrm{OPT}\) for \(k\)-median, and a 16-approximation to the fair radius for all inliers [2412.10923].

"Local Search-based Individually Fair Clustering with Outliers" [2510.06130] studies the same setting using a randomized local-search framework. It adopts the Jung et al. fairness requirement \(d(v,S)\le \delta(v)\) for non-outlier points, allows a relaxation \(d(v,S)\le \gamma \delta(v)\), and introduces BaseCent to identify fairness-based outliers before running LSFO with constrained local search. The method discards outliers in two stages—fairness-based and cost-based—and proves an \(O(1)\)-approximation to the cost of the optimal \((\gamma,k,m)\)-fair clustering. The total number of outliers is bounded by
\[
m + O\left(\frac{mk}{\epsilon}\log(n\Delta)\right),
\]
and the runtime is
\[
O\left( \frac{ndk^3}{\epsilon^2} \log^2 (n\Delta) + \frac{n^2k^2}{\epsilon} \log(n\Delta) \right).
\]
A key lemma states that if the \(n-m\) points are covered by \(k+r\) \(\gamma\)-anchor zones, then a set of \(k\) \((\gamma+2)\)-anchor zones can cover \(n-2m\) points [2510.06130].

Individual fairness has also been extended beyond point-set clustering. "Towards Cohesion-Fairness Harmony: Contrastive Regularization in Individual Fair Graph Clustering" [2402.10756] proposes iFairNMTF, which minimizes
\[
\|\mathbf{A}-\mathbf{HWH}^\top\|_F^2 + \lambda\,\operatorname{Tr}(\mathbf{H}^\top \mathbf{L}\mathbf{H}),
\]
where \(\mathbf{L}\) is the Laplacian of a contrastive matrix built from same-group and different-group indicators. The second term is a contrastive fairness regularizer, and \(\lambda\) controls the cohesion-fairness trade-off [2402.10756].

Representation learning offers another route. "Learning Certified Individually Fair Representations" [2002.10312] learns an encoder \(f_\theta\) such that
\[
\phi(x,x') \implies \|f_\theta(x)-f_\theta(x')\|_\infty \le \delta.
\]
For clustering, the paper states that if a clustering algorithm does not split \(\ell_\infty\)-balls of radius \(\epsilon\), then similar individuals will always be assigned to the same cluster. This is a sufficient condition rather than a clustering-specific optimization criterion, but it places individually fair clustering in a broader certified-fairness pipeline [2002.10312].

## 6. Empirical findings and conceptual relations

Empirical studies consistently report that fairness-aware methods can improve fairness substantially while remaining competitive on clustering objectives, although the precise trade-off depends on the fairness notion. In the fair-radius local-search study, the proposed algorithm improves clustering cost over FairKCenter by factors \(1.4\)–\(2.25\) for \(k\)-median and \(1.73\)–\(2.93\) for \(k\)-means, while increasing the fairness ratio by factors \(1.13\)–\(1.5\) and \(1.14\)–\(1.85\), respectively [2002.06742]. In the feature-based setting, LP-FAIR is reported to be 12.5% fairer on average than the fairest baseline and 34.5% lower in normalized cost on average than the best cost baseline, with more than 95.5% of points satisfying the feature-based fairness constraint across datasets [2109.04554].

Practical scalability has become an explicit benchmark. The scalable local-search algorithm LSPP runs on datasets up to 581,012 points in less than 2 hours and empirically attains fairness bound ratios in the range \(1.1\)–\(3.1\), whereas vanilla \(k\)-means can violate the fair-radius bound by factors of 60–90 [2402.06730]. The LP-based Fair-Round algorithm reports that clustering cost is often within 1% of the LP lower bound and never more than 15%, that the maximum fairness violation is at most 1.27 in the reported experiments, and that more than 80% of points are perfectly fair; sparsification lowers runtime by orders of magnitude [2106.12150].

Outlier-aware studies report similar patterns. The LP-based IFXO method yields lower cost than clustering without outlier removal and outperforms isolation forest plus fair clustering baselines in both cost and fairness metrics [2412.10923]. The local-search-based LSFO method typically attains the lowest \(k\)-means cost among compared methods, with fairness violation ratio \(\rho\) often between \(1.2\) and \(1.9\), practical outlier counts typically less than \(2m\) when \(m=1\%\) of \(n\), and stable behavior under parameter changes [2510.06130].

The literature also documents several important non-equivalences. Group-fair clustering is not necessarily individually fair: this is shown both in the average-dissimilarity line of work and in the distributional framework, where standard group-fair algorithms may satisfy demographic constraints while violating a large fraction of individual-fairness constraints [2006.04960] [2006.12589]. Conversely, individual fairness can imply group fairness only under additional conditions; in the distributional model, the bound
\[
\MAD_r \le |G_r| \cdot d_{\mathrm{EM}}(\nu_{G_r}, \nu_V)
\]
links group-level deviation to an earthmover-distance condition on the protected group’s spread relative to the whole population [2006.12589].

Taken together, these results show that individually fair clustering is not a single optimization problem but a family of related problems that differ in what counts as similar individuals, what counts as similar treatment, and whether fairness is enforced through deterministic assignments, center accessibility, service parity, co-clustering, or probabilistic output distributions. The main technical trajectory of the field has therefore been to turn each fairness notion into a tractable surrogate: critical-ball coverage, matroid-constrained facility location, LP-based fair assignment, dynamic programming on structured domains, and outlier-aware local search.

Source: https://www.emergentmind.com/topics/individually-fair-clustering