Optimal facility location

The study of facility location problems (FLP), also known as location analysis, is a branch of operations research and computational geometry concerned with the optimal placement of facilities to minimize transportation costs while considering factors like avoiding placing hazardous materials near housing, and competitors' facilities. The techniques also apply to cluster analysis.

Minimum facility location

A simple facility location problem is the Weber problem, in which a single facility is to be placed, with the only optimization criterion being the minimization of the weighted sum of distances from a given set of point sites. More complex problems considered in this discipline include the placement of multiple facilities, constraints on the locations of facilities, and more complex optimization criteria.

In a basic formulation, the facility location problem consists of a set of potential facility sites L where a facility can be opened, and a set of demand points D that must be serviced. The goal is to pick a subset F of facilities to open, to minimize the sum of distances from each demand point to its nearest facility, plus the sum of opening costs of the facilities.

The facility location problem on general graphs is NP-hard to solve optimally, by reduction from (for example) the set cover problem. A number of approximation algorithms have been developed for the facility location problem and many of its variants.

Without assumptions on the set of distances between clients and sites (in particular, without assuming that the distances satisfy the triangle inequality), the problem is known as non-metric facility location and can be approximated to within a factor O(log n). ^[1] This factor is tight, via an approximation-preserving reduction from the set cover problem.

If we assume distances between clients and sites are undirected and satisfy the triangle inequality, we are talking about a metric facility location (MFL) problem. The MFL is still NP-hard and hard to approximate within factor better than 1.463. ^[2] The currently best known approximation algorithm achieves approximation ratio of 1.488. ^[3]

Minimax facility location

The minimax facility location problem seeks a location which minimizes the maximum distance to the sites, where the distance from one point to the sites is the distance from the point to its nearest site. A formal definition is as follows:

Given a point set P ⊂ ${\displaystyle \mathbb {R} ^{d}}$, find a point set S ⊂ ${\displaystyle \mathbb {R} ^{d}}$, |S| = k, so that max_{p ∈ P}(min_{q ∈ S}(d(p, q)) ) is minimized.

In the case of the Euclidean metric for k = 1, it is known as the smallest enclosing sphere problem or 1-center problem. Its study traced at least to the year of 1860.

Further information: Smallest enclosing circle and Bounding sphere

NP hardness

It has been proven that exact solution of k-center problem is NP hard. ^{[4] [5] [6]} Approximation to the problem was found to be also NP hard when the error is small. The error level in the approximation algorithm is measured as an approximation factor, which is defined as the ratio between the approximation and the optimum. It's proved that the k-center problem approximation is NP hard when approximation factor is less than 1.822 (dimension = 2) ^[7] or 2 (dimension > 2). ^[6]

Algorithms

Exact solver

There exist algorithms to produce exact solutions to this problem. One exact solver runs in time ${\displaystyle n^{O({\sqrt {k}})}}$. ^{[8] [9]}

1 + ε approximation

1 + ε approximation is to find a solution with approximation factor no greater than 1 + ε. This approximation is NP hard as ε is arbitrary. One approach based on the coreset concept is proposed with execution complexity of ${\displaystyle O(2^{O(k\log k/\varepsilon ^{2})}dn)}$. ^[10] As an alternative, another algorithm also based on core sets is available. It runs in ${\displaystyle O(k^{n})}$. ^[11] The author claims that the running time is much less than the worst case and thus it's possible to solve some problems when k is small (say k < 5).

Farthest-point clustering

For the hardness of the problem, it's impractical to get an exact solution or precise approximation. Instead, an approximation with factor = 2 is widely used for large k cases. The approximation is referred to as the farthest-point clustering (FPC) algorithm, or farthest-first traversal. ^[6] The algorithm is quite simple: pick any point from the set as one center; search for the farthest point from remaining set as another center; repeat the process until k centers are found.

It is easy to see that this algorithm runs in linear time. As approximation with factor less than 2 is proved to be NP hard, FPC was regarded as the best approximation one can find.

As per the performance of execution, the time complexity is later improved to O(n log k) with box decomposition technique. ^[7]

Maxmin facility location

The maxmin facility location or obnoxious facility location problem seeks a location which maximizes the minimum distance to the sites. In the case of the Euclidean metric, it is known as the largest empty sphere problem. The planar case (largest empty circle problem) may be solved in optimal time Θ(n log n). ^{[12] [13]}

Integer programming formulations

Facility location problems are often solved as integer programs. In this context, facility location problems are often posed as follows: suppose there are ${\displaystyle n}$ facilities and ${\displaystyle m}$ customers. We wish to choose (1) which of the ${\displaystyle n}$ facilities to open, and (2) which (open) facilities to use to supply the ${\displaystyle m}$ customers, in order to satisfy some fixed demand at minimum cost. We introduce the following notation: let ${\displaystyle f_{i}}$ denote the (fixed) cost of opening facility ${\displaystyle i}$, for ${\displaystyle i=1,\dots ,n}$. Let ${\displaystyle c_{ij}}$denote the cost to ship a product from facility ${\displaystyle i}$ to customer ${\displaystyle j}$ for ${\displaystyle i=1,\dots ,n}$ and ${\displaystyle j=1,\dots ,m}$. Let ${\displaystyle d_{j}}$ denote the demand of customer ${\displaystyle j}$ for ${\displaystyle j=1,\dots ,m}$. Further suppose that each facility has a maximum output. Let ${\displaystyle u_{i}}$ denote the maximum amount of product that can be produced by facility ${\displaystyle i}$, that is, let ${\displaystyle u_{i}}$ denote the capacity of facility ${\displaystyle i}$. The remainder of this section follows ^[14]

Capacitated facility location

In our initial formulation, introduce a binary variable ${\displaystyle x_{i}}$ for ${\displaystyle i=1,\dots ,n}$, where ${\displaystyle x_{i}=1}$ if facility ${\displaystyle i}$ is open, and ${\displaystyle x_{i}=0}$ otherwise. Further introduce the variable ${\displaystyle y_{ij}}$ for ${\displaystyle i=1,\dots ,n}$ and ${\displaystyle j=1,\dots ,m}$ which represents the fraction of the demand ${\displaystyle d_{j}}$ filled by facility ${\displaystyle i}$. The so-called capacitated facility location problem is then given by$${\displaystyle {\begin{array}{rl}\min &\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{m}c_{ij}d_{j}y_{ij}+\sum _{i=1}^{n}f_{i}x_{i}\\{\text{s.t.}}&\displaystyle \sum _{i=1}^{n}y_{ij}=1{\text{ for all }}j=1,\dots ,m\\&\displaystyle \sum _{j=1}^{m}d_{j}y_{ij}\leqslant u_{i}x_{i}{\text{ for all }}i=1\dots ,n\\&y_{ij}\geqslant 0{\text{ for all }}i=1,\dots ,n{\text{ and }}j=1,\dots ,m\\&x_{i}\in \{0,1\}{\text{ for all }}i=1,\dots ,n\end{array}}}$$

Note that the second set of constraints ensures that if ${\displaystyle x_{i}=0}$, that is, facility ${\displaystyle i}$ isn't open, then ${\displaystyle y_{ij}=0}$ for all ${\displaystyle j}$, that is, no demand for any customer can be filled from facility ${\displaystyle i}$.

Uncapacitated facility location

A common case of the capacitated facility location problem above is the case when ${\displaystyle u_{i}=+\infty }$ for all ${\displaystyle i=1,\dots ,n}$. In this case, it is always optimal to satisfy all of the demand from customer ${\displaystyle j}$ from the nearest open facility. Because of this, we may replace the continuous variables ${\displaystyle y_{ij}}$ from above with the binary variables ${\displaystyle z_{ij}}$, where ${\displaystyle z_{ij}=1}$ if customer ${\displaystyle j}$ is supplied by facility ${\displaystyle i}$, and ${\displaystyle z_{ij}=0}$ otherwise. The uncapacitated facility location problem is then given by$${\displaystyle {\begin{array}{rl}\min &\displaystyle \sum _{i=1}^{n}\sum _{j=1}^{m}c_{ij}d_{j}z_{ij}+\sum _{i=1}^{n}f_{i}x_{i}\\{\text{s.t.}}&\displaystyle \sum _{i=1}^{n}z_{ij}=1{\text{ for all }}j=1,\dots ,m\\&\displaystyle \sum _{j=1}^{m}z_{ij}\leqslant Mx_{i}{\text{ for all }}i=1\dots ,n\\&z_{ij}\in \{0,1\}{\text{ for all }}i=1,\dots ,n{\text{ and }}j=1,\dots ,m\\&x_{i}\in \{0,1\}{\text{ for all }}i=1,\dots ,n\end{array}}}$$

where ${\displaystyle M}$ is a constant chosen to be suitably large. The choice of ${\displaystyle M}$ can affect computation results—the best choice in this instance is obvious: take ${\displaystyle M=m}$. Then, if ${\displaystyle x_{i}=1}$, any choice of the ${\displaystyle z_{ij}}$ will satisfy the second set of constraints.

Another formulation possibility for the uncapacitated facility location problem is to disaggregate the capacity constraints (the big-${\displaystyle M}$ constraints). That is, replace the constraints$${\displaystyle \sum _{j=1}^{m}z_{ij}\leqslant Mx_{i}{\text{ for all }}i=1,\dots ,n}$$with the constraints$${\displaystyle z_{ij}\leqslant x_{i}{\text{ for all }}i=1,\dots ,n{\text{ and }}j=1,\dots ,m}$$In practice, this new formulation performs significantly better, in the sense that it has a tighter Linear programming relaxation than the first formulation. ^[14] Notice that summing the new constraints together yields the original big-${\displaystyle M}$ constraints. In the capacitated case, these formulations are not equivalent. More information about the uncapacitated facility location problem can be found in Chapter 3 of "Discrete location theory". ^[15]

Applications

Healthcare

In healthcare, incorrect facility location decisions have a serious impact on the community beyond simple cost and service metrics; for instance, hard-to-access healthcare facilities are likely to be associated with increased morbidity and mortality. From this perspective, facility location modeling for healthcare is more critical than similar modeling for other areas. ^[16]

Solid waste management

Municipal solid waste management still remains a challenge for developing countries because of increasing waste production and high costs associated with waste management. Through the formulation and exact resolution of a facility location problem it is possible to optimize the location of landfills for waste disposal. ^[17]

Clustering

A particular subset of cluster analysis problems can be viewed as facility location problems. In a centroid-based clustering problem, the objective is to partition ${\displaystyle n}$ data points (elements of a common metric space) into equivalence classes—often called colors—such that points of the same color are close to one another (equivalently, such that points of different colors are far from one another). ^[18]

To see how one might view (read "transform" or "reduce") a centroid-based clustering problem as a (metric) facility location problem, view each data point in the former as a demand point in the latter. Suppose that the data to be clustered are elements of a metric space ${\displaystyle M}$ (e.g. let ${\displaystyle M}$ be ${\displaystyle p}$-dimensional Euclidean space for some fixed ${\displaystyle p}$). In the facility location problem that we are constructing, we permit facilities to be placed at any point within this metric space ${\displaystyle M}$; this defines the set of allowed facility locations ${\displaystyle L}$. We define the costs ${\displaystyle c_{\ell ,d}}$ to be the pairwise distances between location-demand point pairs (e.g., see metric k-center). In a centroid-based clustering problem, one partitions the data into ${\displaystyle k}$ equivalence classes (i.e. colors) each of which has a centroid. Let us see how a solution to our constructed facility location problem also achieves such a partition. A feasible solution is a non-empty subset ${\displaystyle L'\subseteq L}$ of ${\displaystyle k}$ locations. These locations in our facility location problem comprise a set of ${\displaystyle k}$ centroids in our centroid-based clustering problem. Now, assign each demand point ${\displaystyle d}$ to the location ${\displaystyle \ell ^{*}}$ that minimizes its servicing-cost; that is, assign the data point ${\displaystyle d}$ to the centroid ${\displaystyle \ell ^{*}:=\mathrm {arg\,min} _{\ell \in L}\{c_{\ell ,d}\}}$ (break ties arbitrarily). This achieves the partitioning provided that the facility location problem's costs ${\displaystyle c_{\ell ,d}}$ are defined such that they are the images of the centroid-based clustering problem's distance function.

The popular algorithms textbook Algorithm Design ^[19] provides a related problem-description and an approximation algorithm. The authors refer to the metric facility location problem (i.e. the centroid-based clustering problem or the metric ${\displaystyle k}$-center problem) as the center selection problem, thereby growing the list of synonyms.

Furthermore, see that in our above definition of the facility location problem that the objective function ${\displaystyle f}$ is general. Specific choices of ${\displaystyle f}$ yield different variants of the facility location problem, and hence different variants of the centroid-based clustering problem. For example, one might choose to minimize the sum of distances from each location to each of its assigned demand points (à la the Weber problem), or one might elect to minimize the maximum of all such distances (à la the 1-center problem).

See also

• Graph center
• Quadratic assignment problem
• Location-allocation
• Dijkstra's algorithm
• List of spatial analysis software
• Competitive facility location game
• Vertex k-center problem
• geometric median

References

1. Hochbaum, D. S. (1982). "Heuristics for the fixed cost median problem". Mathematical Programming . 22: 148–162. doi:10.1007/BF01581035. S2CID   3451944.
2. Guha, S.; Khuller, S. (1999). "Greedy Strikes Back: Improved Facility Location Algorithms". Journal of Algorithms. 31: 228–248. CiteSeerX   10.1.1.47.2033 . doi:10.1006/jagm.1998.0993. S2CID   5363214.
3. Li, S. (2011). "A 1.488 Approximation Algorithm for the Uncapacitated Facility Location Problem". Automata, Languages and Programming. LNCS. Vol. 6756. pp. 77–88. CiteSeerX   10.1.1.225.6387 . doi:10.1007/978-3-642-22012-8_5. ISBN   978-3-642-22011-1.
4. Fowler, R. J.; Paterson, M. S.; Tanimoto, S. L. (1981), "Optimal packing and covering in the plane are NP-complete", Information Processing Letters, 12 (3): 133–137, doi:10.1016/0020-0190(81)90111-3 .
5. Megiddo, Nimrod; Tamir, Arie (1982), "On the complexity of locating linear facilities in the plane" (PDF), Operations Research Letters, 1 (5): 194–197, doi:10.1016/0167-6377(82)90039-6 .
6. Gonzalez, Teofilo (1985), "Clustering to minimize the maximum intercluster distance", Theoretical Computer Science , 38: 293–306, doi: 10.1016/0304-3975(85)90224-5 .
7. Feder, Tomás; Greene, Daniel (1988), "Optimal algorithms for approximate clustering", Proceedings of the twentieth annual ACM symposium on Theory of computing - STOC '88, pp. 434–444, doi:10.1145/62212.62255, ISBN   0897912640, S2CID   658151
8. HWang, R. Z.; Lee, R. C. T.; Chang, R. C. (1993), "The slab dividing approach to solve the Euclidean p-center problem", Algorithmica, 9 (1): 1–22, doi:10.1007/BF01185335, S2CID   5680676
9. HWang, R. Z.; Chang, R. C.; Lee, R. C. T. (1993), "The generalized searching over separators strategy to solve some NP-Hard problems in subexponential time", Algorithmica, 9 (4): 398–423, doi:10.1007/bf01228511, S2CID   2722869
10. Bādoiu, Mihai; Har-Peled, Sariel; Indyk, Piotr (2002), "Approximate clustering via core-sets", Proceedings of the thiry-fourth annual ACM symposium on Theory of computing (PDF), pp. 250–257, doi:10.1145/509907.509947, ISBN   1581134959, S2CID   5409535
11. Kumar, Pankaj; Kumar, Piyush (2010), "Almost optimal solutions to k-clustering problems" (PDF), International Journal of Computational Geometry & Applications, 20 (4): 431–447, doi:10.1142/S0218195910003372
12. Franco P. Preparata and Michael Ian Shamos (1985). Computational Geometry – An Introduction . Springer-Verlag. ISBN   978-0-387-96131-6. 1st edition; 2nd printing, corrected and expanded, 1988; Russian translation, 1989., p. 256
13. G. T. Toussaint, "Computing largest empty circles with location constraints," International Journal of Computer and Information Sciences, vol. 12, No. 5, October, 1983, pp. 347–358.
14. Conforti, Michele; Cornuéjols, Gérard; Zambelli, Giacomo (2014). Integer Programming. Graduate Texts in Mathematics. Vol. 271. doi:10.1007/978-3-319-11008-0. ISBN   978-3-319-11007-3.
15. Discrete location theory. Mirchandani, Pitu B., Francis, R. L. New York: Wiley. 1990. ISBN   9780471892335. OCLC   19810449.
16. Ahmadi-Javid, A.; Seyedi, P.; Syam, S. (2017). "A Survey of Healthcare Facility Location". Computers & Operations Research. 79: 223–263. doi:10.1016/j.cor.2016.05.018.
17. Franco, D. G. B.; Steiner, M. T. A.; Assef, F. M. (2020). "Optimization in waste landfilling partitioning in Paraná State, Brazil". Journal of Cleaner Production. 283: 125353. doi:10.1016/j.jclepro.2020.125353. S2CID   229429742.
18. Hastie, Trevor; Tibshirani, Robert; Friedman, Jerome (2009). The elements of statistical learning (Second ed.). Springer.
19. Kleinberg, Jon; Tardos, Éva (2006). Algorithm Design. Pearson.

External links

• EWGLA EURO Working Group on Locational Analysis.
• INFORMS section on location analysis, a professional society concerned with facility location.
• Bibliography on facility location collected by Trevor Hale, containing over 3400 articles.
• Library of location algorithms
• Web-based facility location utility (single facility)
• Facility Location Optimizer, a MATLAB-based tool for solving facility location problems.