---
title: Facility Location with Fair Outliers
url: https://www.emergentmind.com/topics/facility-location-with-fair-outliers
type: topic
---

# Facility Location with Fair Outliers

Facility location with fair outliers denotes a family of robust service-allocation problems in which some clients may be excluded from service, but the exclusion mechanism itself is constrained by fairness. In the current literature, fairness appears in two distinct forms. The first is **group-wise fairness**: clients are partitioned into groups, each group has its own outlier budget, and the objective is classical facility-opening plus assignment cost, or assignment cost alone under a \(k\)-facility bound [2508.02572]. The second is **individual fairness**: every retained point must have a center within a locally defined neighborhood radius, while a bounded set of points may be discarded as outliers [2201.04849]. LP-based and local-search-based formulations for fair \(k\)-means and \(k\)-median with outliers extend this second view to cost-sensitive clustering [2412.10923].

## 1. Formal problem family

In the group-wise formulation, classical Facility Location is specified by a metric space \((\mathcal M,d)\), a facility set \(\mathcal F\subseteq \mathcal M\), a client set \(C\subseteq \mathcal M\), and opening costs \(f_i\ge 0\) for \(i\in\mathcal F\). One opens a subset \(F\subseteq \mathcal F\) and assigns each client to its nearest open facility, minimizing
\[
\sum_{i\in F} f_i+\sum_{j\in C} d_{j,F},
\qquad
d_{j,F}=\min_{i\in F} d(i,j).
\]
Facility Location with Fair Outliers (FLFO) extends this by partitioning clients into \(\omega\) disjoint groups
\[
C_1,\dots,C_\omega,
\]
and associating to each group \(g\in[\omega]\) an outlier budget \(\ell_g\). The algorithm chooses opened facilities \(F\subseteq\mathcal F\) and, for each group, an outlier set \(C'_g\subseteq C_g\) with \(|C'_g|\le \ell_g\), minimizing
\[
\sum_{i \in F} f_i + \sum_{g=1}^{\omega} \sum_{j \in C_g \setminus C'_g} d_{j,F}.
\]
Equivalently, at least \(|C_g|-\ell_g\) clients from each group must be served. The model in this line of work uses only disjoint groups and does not introduce overlapping groups or multiple protected attributes per client beyond that partition model [2508.02572].

The corresponding \(k\)-median formulation replaces opening costs by a cardinality constraint \(|F|\le k\). With the same groups and outlier budgets, the objective becomes
\[
\sum_{g=1}^{\omega}\sum_{j\in C_g\setminus C'_g} d_{j,F}
\qquad\text{subject to } |F|\le k.
\]
This yields a common template: serve clients cheaply while controlling outliers per group [2508.02572].

A different fairness line studies outliers through **local service entitlements** rather than group quotas. In individually fair \(k\)-center with outliers, one may choose at most \(k\) centers and at most \(q\) outliers, and each served vertex \(i\) is evaluated against its outlier-related neighborhood radius \(NR_q(i)\), defined as the distance to its \(\left\lceil \frac{n-q}{k}\right\rceil\)-th nearest neighbor. The objective is to minimize the maximum fairness ratio
\[
\max_{i\in V\setminus O}\frac{d_{\sigma(i)i}}{NR_q(i)}.
\]
Here fairness is not demographic; it is individual and density-sensitive [2201.04849].

LP-based individually fair \(k\)-means and \(k\)-median with outliers use the same general idea. Given an outlier set \(Z\subseteq X\) with \(|Z|\le m\), every non-outlier \(v\in X\setminus Z\) must be assigned to a center within a client-specific radius \(r(v)\), where \(r(v)\) is the distance to the \(\frac{n}{k}\)-th nearest neighbor. The optimization then minimizes \(\sum d(v,u)^p\) over non-outliers, with \(p=2\) for \(k\)-means and \(p=1\) for \(k\)-median [2412.10923].

## 2. Mathematical formulations and structural obstacles

For FLFO, a natural integer program uses binary variables \(y_i\) for opening facility \(i\), \(x_{ij}\) for assigning client \(j\) to facility \(i\), and \(z_j\) for declaring client \(j\) an outlier. The formulation is
\[
\min \sum_{i\in\mathcal F} f_i y_i + \sum_{j\in C}\sum_{i\in\mathcal F} d_{ij} x_{ij}
\]
subject to
\[
\sum_{i\in\mathcal F} x_{ij} + z_j \ge 1 \quad \forall j\in C,
\]
\[
x_{ij}\le y_i \quad \forall i\in\mathcal F,\forall j\in C,
\]
\[
\sum_{j\in C_g} z_j \le \ell_g \quad \forall g\in[\omega],
\]
with integrality relaxed to \(x_{ij},y_i,z_j\in[0,1]\) in the LP relaxation [2508.02572].

This LP is natural but structurally weak. Even when \(\omega=1\), it has an unbounded integrality gap. The canonical example has one facility of opening cost \(f\), \(M\) co-located clients, and outlier budget \(M-1\). Fractionally, the LP can open the facility to extent \(1/M\) and serve every client fractionally at total cost \(f/M\). Any integral solution that serves one client must open the facility fully and pay \(f\). As \(M\to\infty\), the gap is unbounded. This eliminates straightforward exact-budget LP-rounding as a generic approach for FLFO [2508.02572].

The individual-fairness line inherits a different structural difficulty from the neighborhood-radius benchmark itself. In the no-outlier setting, “A Center in Your Neighborhood” defines
\[
NR_{P,k}(x)=\min\{r:\ |B_r(x)\cap P|\ge |P|/k\}
\]
and optimizes
\[
\min_{|S|\le k}\max_{i\in P}\frac{d(i,S)}{NR_{P,k}(i)}.
\]
That problem admits a universal factor-2 algorithm in arbitrary metrics, but deciding whether \(1\)-fairness is achievable is NP-complete [1908.09041]. The outlier-aware individual-fairness models therefore start from an already nontrivial fairness geometry.

Taken together, these results explain why fair-outlier algorithms frequently become bicriteria. Exact compliance with outlier budgets or exact fairness radii interacts poorly with the natural LPs and with the combinatorial geometry of local service constraints. This suggests that controlled relaxation of either cost, fairness, or outlier counts is not merely an implementation artifact, but a recurring structural response.

## 3. Group-wise facility location and \(k\)-median with fair outliers

The main current group-wise result is a bicriteria approximation for FLFO. For any fixed \(\epsilon>0\), there is a polynomial-time \(O(1/\epsilon)\)-approximation algorithm that violates each group’s outlier budget by at most a factor \(1+2\epsilon\). The starting point is the LP optimum \(\rho^*=\langle x^*,y^*,z^*\rangle\). Clients with \(z_j^*\ge 1-\epsilon\) are placed into
\[
C_o=\{j\in C: z_j^*\ge 1-\epsilon\},
\]
and are committed as integral outliers. The remaining clients form \(C_r=C\setminus C_o\). This thresholding yields the key fairness inequality
\[
\sum_{j\in C_g}\hat z_j \le \frac{1}{1-\epsilon}\sum_{j\in C_g} z_j^* \le (1+2\epsilon)\ell_g,
\]
while not increasing LP cost [2508.02572].

The thresholding step is then followed by renormalization. Every client \(j\in C_r\) has residual assignment mass at least \(\epsilon\), so assignments can be rescaled to obtain a feasible fractional solution to the standard Facility Location LP on the reduced client set. The resulting cost increases by at most \(1/\epsilon\). Any \(\alpha\)-approximation for ordinary Facility Location therefore yields an \((\alpha/\epsilon)\)-approximation for FLFO, with the same \((1+2\epsilon)\) factor violation in every group’s outlier budget. Combining this with a constant-factor approximation for classical Facility Location gives the stated \(O(1/\epsilon)\) bicriteria guarantee [2508.02572].

The same paper develops a parallel bicriteria result for \(k\)-Median with Fair Outliers. There is a polynomial-time
\[
4\left(1+\frac{\omega}{\epsilon}\right)
\]
factor approximation with \((\omega+\epsilon)\) violation in outliers per group. Here the technique is different: the reduction is to \(k\)-Median with Penalties, with penalties set group-wise as
\[
p_j=\frac{\mathrm{OPT}_C}{\gamma\,\ell_g}
\qquad\text{for } j\in C_g.
\]
This improves earlier fair \(k\)-median work by removing dependence on \(k\) from the outlier-budget violation; the prior result of Almanza et al. allowed a \((3k+2)\)-factor outlier violation per group [2508.02572].

The group-wise line is also shaped by practicality. Earlier FLFO results of Inamdar–Varadarajan and Bajpai et al. give proper approximations, but rely on exponentially large LPs solved through ellipsoid-based machinery. The compact-LP bicriteria formulation changes the guarantee—small multiplicative violation in each group’s outlier budget replaces exact compliance—but it is polynomial-time and directly implementable [2508.02572].

## 4. Individual fairness with outliers

The individually fair \(k\)-center with outliers problem, IF\(k\)CO, is a clean min-max formulation. A feasible solution chooses centers \(S\subseteq V\) with \(|S|\le k\), outliers \(O\subseteq V\) with \(|O|\le q\), and an assignment \(\sigma:V\setminus O\to S\). The fairness ratio of a served vertex \(i\) is
\[
\frac{d_{\sigma(i)i}}{NR_q(i)},
\]
and the objective is to minimize the maximum such ratio over served vertices. A greedy algorithm that repeatedly selects the point of minimum \(NR_q(i)\) and removes all points within distance \(2\,NR_q(i)\) is proved to be a 4-approximation. A refined practical variant performs binary search over a parameter \(\beta\in[1,2]\) and retains the same 4-approximation guarantee while often achieving smaller empirical ratios. The center and outlier budgets are respected exactly; there is no bicriteria violation in \(k\) or \(q\) [2201.04849].

A different LP-based line targets individually fair \(k\)-means and \(k\)-median with outliers. The LP uses assignment variables \(x_{vu}\), opening variables \(y_u\), and outlier indicators \(z_v\), with the service-access constraint
\[
x_{vu}=0 \quad \text{if } d(u,v)>\alpha r(v).
\]
After solving the LP with \(\alpha=1\), the procedure **OutRound** thresholds outliers, deletes any center mass placed on points designated as outliers, reroutes their assignments to the nearest surviving non-outlier, and then invokes an existing fair-rounding subroutine on the inliers. The key rerouting lemma is
\[
d(v,u')\le 2d(v,u_{OUT})\le 2r(v),
\]
which yields
\[
LP_{\alpha=2}(x',y',z') \le 3\, LP_{\alpha=1}(x^*,y^*,z^*).
\]
The resulting guarantees are a 12-approximation for fair \(k\)-means with outliers, a 24-approximation for fair \(k\)-median with outliers, and a 16-approximation to the fair radius for inliers. The important caveat is that the final number of detected outliers is not theoretically bounded by \(m\), even though the LP satisfies \(\sum_v z_v\le m\) fractionally [2412.10923].

More recently, local search has been adapted to individually fair clustering with outliers through the anchor-zone framework. The algorithm first constructs fairness-aware anchor zones, discards an initial set of fairness-based outliers, and then performs local search while maintaining the invariant that every anchor zone contains at least one center. The main theorem gives an \(O(1)\)-approximation for fixed \(0<\epsilon\le 1/4\), and the analysis bounds the total number of discarded outliers by
\[
m + \frac{mk}{\epsilon}\log(n\Delta),
\]
where \(\Delta\) is the aspect ratio. The appendix states an explicit constant of about \(274\) for the \(k\)-means objective, so the guarantee is again bicriteria in practice and in theory [2510.06130].

These individual-fairness papers show that “fair outliers” need not mean demographic parity in exclusion. It can also mean that once outliers are removed, every retained client must still receive service compatible with a local density benchmark.

## 5. Complexity and relation to adjacent facility-location models

The group-wise FLFO landscape includes a parameterized hardness barrier. Assuming ETH, Facility Location with Fair Outliers is W[1]-hard when parameterized by the number of groups \(\omega\), and the same holds for \(k\)-Median with Fair Outliers. The reduction proceeds from a multidimensional subset-sum type problem, with the dimension becoming the number of groups. Consequently, one should not expect algorithms of running time \(f(\omega)\cdot \mathrm{poly}(n)\) for arbitrary \(\omega\) unless ETH fails [2508.02572].

The surrounding robust-facility-location literature clarifies what is specific to fairness and what is inherited from outliers alone.

| Setting | Representative guarantee | Relevance |
|---|---|---|
| UFL/\(k\)-Median/\(k\)-Means with outliers | PTAS for uniform-cost UFL with outliers on doubling and minor-closed metrics; bicriteria PTAS for \(k\)-Median/\(k\)-Means there; \((25+\epsilon,1+\epsilon)\) for general-metric \(k\)-Means with outliers [1707.04295] | Robust baseline without fairness |
| CFLPO/C\(k\)FLPP | Constant-factor LP-rounding with slight capacity and outlier or cardinality violations [2012.07135] | Shows how capacities and outliers interact |
| \(k\)FLO and LBFLO | \(11\)-approximation for \(k\)FLO with at most \(k+1\) facilities; tri-criteria tradeoff for LBFLO [2107.00403] | Global outlier budgets, but no fair allocation of outliers |
| Capacitated \(k\)-Facility Location with Outliers; \((\alpha,\beta)\)-fair clustering with outliers | FPT reduction framework; \((3+\epsilon)\) for capacitated \(k\)-facility location with outliers in arbitrary metrics; fair \(k\)-facility location with outliers is stated as an extension of the same framework [2305.01471] | Closest bridge between capacities, fairness, and outliers |

The fair-outlier topic also borders group-fair facility-location models without outliers. Strategyproof single-facility location on the line has been studied under objectives such as maximum total group cost, maximum average group cost, and intergroup/intragroup fairness, but those models do not include an outlier budget or partial service [2107.05175]. They provide fairness objectives and impossibility patterns that are relevant conceptually, but not direct algorithms for fair outlier control.

A broad lesson from these adjacent models is that robust facility location already becomes bicriteria once capacities, lower bounds, or cardinality limits are added. This suggests that fair-outlier facility location inherits not one source of approximation loss, but several: robustness, fairness, and classical facility-opening structure.

## 6. Empirical behavior and practical implications

The FLFO experimental results are unusually strong relative to the worst-case theory. On the Adult dataset grouped by sex, the Bank dataset grouped by marital status, and a synthetic two-group dataset, the fair facility-location methods LPR-F and GDF-F achieve unfairness essentially equal to \(1\) in almost all cases; the only reported exception is the Bank dataset, where the maximum unfairness is \(1.12\). By contrast, the non-fair baselines can be much more skewed: on Adult, unfairness reaches \(1.67\) for LPR-NF and \(1.23\) for GDF-NF; on Bank, \(1.19\) and \(1.28\); on the synthetic data, as high as \(14\) and \(13\). The cost increase from fairness is reported as negligible: LPR-F remains close to the LP optimum, and GDF-F has cost comparable to the classical greedy baseline. The same section also reports that the standard LP, despite its unbounded worst-case gap, is usually nearly integral in practice [2508.02572].

The paper’s \(k\)-median experiments show a similar pattern. Using 500-point samples with \(k=5\), the fair method R+LS-F attains unfairness close to \(1\), while non-fair alternatives can reach \(1.12\) and \(1.25\) on Adult, \(1.50\) and \(2\) on Bank, and \(13\) on the synthetic data. Fairness increases cost somewhat, but not prohibitively [2508.02572].

For individually fair \(k\)-center with outliers, synthetic experiments report a maximum outlier-related fairness ratio of only \(1.31\), far below the proven upper bound of \(2\). On the Shenzhen POI dataset, the outlier-aware fair algorithm tends to place centers in dense areas, whereas the naive non-outlier-aware baseline does so less effectively [2201.04849].

The LP-based individually fair \(k\)-means/\(k\)-median with outliers framework evaluates on sampled UCI datasets with injected outliers and compares against an Isolation Forest plus FairRound baseline. The reported behavior is that optimization-aware outlier identification yields lower cost and smaller maximum fair radius than the baseline, but the number of detected outliers can substantially exceed the nominal budget, empirically confirming the theoretical caveat of uncontrolled final outlier count [2412.10923].

The local-search-based individually fair clustering method also reports strong empirical performance. On Adult, Bank, and Skin, it achieves lower \(k\)-means cost than the listed baselines, while maintaining reported fairness ratios such as \(1.2\), \(1.5\), and \(1.91\). The observed discarded-outlier counts—495, 707, and 14000 in the reported runs—again illustrate that practical success is accompanied by bicriteria outlier behavior rather than exact-budget compliance [2510.06130].

Across these results, one repeated empirical pattern is clear: fairness constraints often change cost far less than unconstrained robust-clustering intuition would suggest, whereas the main practical tension is more often between exact budget compliance and algorithmic tractability than between fairness and objective value itself.

Source: https://www.emergentmind.com/topics/facility-location-with-fair-outliers