Facility Location Problem Research Articles

Facility Location Problem: Modeling Joint Disruptions Using Subordination

We study the facility location problem with disruptions where the objective is to choose a set of locations that minimizes the sum of expected servicing and setup costs. Disruptions can affect multiple locations simultaneously and are caused by multiple factors like geography, supply chain characteristics, politics, and ownership. Accounting for the various factors when modeling disruptions is challenging due to a large number of required parameters, the lack of calibration methodologies, the sparsity of disruption data, and the number of scenarios to be considered in the optimization. Because of these reasons, existing models neglect dependence or prespecify the dependence structures. Using partially subordinated Markov chains, we present a comprehensive approach that starts from disruption data, models dependencies, calibrates the disruption model, and optimizes location choices. We construct a metric and a calibration algorithm that learns from the data the strength of dependence, the number of necessary factors (subordinators), and the locations each subordinator affects. We prove that our calibration approach yields consistent estimates of the model parameters. Then, we introduce a variant of the standard approach to the underlying optimization problem, which leverages partially subordinated Markov chains to solve it quickly and precisely. Finally, we demonstrate the efficacy of our approach using twelve different disruption data sets. Our calibrated parameters are robust, and our optimization algorithm performs better than the simulation-based algorithm. The solutions from our model for disruptions have lower costs than those from other disruption models. Our approach allows for better modeling of disruptions from historical data and can be adapted to other problems in logistics, like the hub location, capacitated facility location, and so on., with joint disruptions. Supplemental Material: The online appendix is available at https://doi.org/10.1287/trsc.2023.0103 .

A Model-Free Approach for Solving Choice-Based Competitive Facility Location Problems Using Simulation and Submodularity

This paper considers facility location problems in which a firm entering a market seeks to open facilities on a subset of candidate locations so as to maximize its expected market share, assuming that customers choose the available alternative that maximizes a random utility function. We introduce a deterministic equivalent reformulation of this stochastic problem as a maximum covering location problem with an exponential number of demand points, each of which is covered by a different set of candidate locations. Estimating the prevalence of these preference profiles through simulation generalizes a sample average approximation method from the literature and results in a maximum covering location problem of manageable size. To solve it, we develop a partial Benders reformulation in which the contribution to the objective of the least influential preference profiles is aggregated and bounded by submodular cuts. This set of profiles is selected by a knee detection method that seeks to identify the best tradeoff between the fraction of the demand that is retained in the master problem and the size of the model. We develop a theoretical analysis of our approach and show that the solution quality it provides for the original stochastic problem, its computational performance, and the automatic profile-retention strategy it exploits are directly connected to the entropy of the preference profiles in the population. Computational experiments on existing and new benchmark sets indicate that our approach dominates the classical sample average approximation method on large instances of the competitive facility location problem, can outperform the best heuristic method from the literature under the multinomial logit model, and achieves state-of-the-art results under the mixed multinomial logit model. We characterize a broader class of problems, which includes assortment optimization, to which the solving methodology and the analyses developed in this paper can be extended. History: Accepted by Andrea Lodi, Area Editor for Design & Analysis of Algorithms—Discrete. Funding: This research was supported by Fonds de Recherche du Québec-Nature et Technologies and Institut de Valorisation des Données through scholarships to R. Legault. E. Frejinger was partially supported by the Canada Research Chair program [Grant 950-232244]. Supplemental Material: The online appendices are available at https://doi.org/10.1287/ijoc.2023.0280 .

An exact method for trilevel hub location problem with interdiction

In this paper, we study the problem of designing a hub network that is robust against deliberate attacks (interdictions). The problem is modeled as a three-level, two-player Stackelberg game, in which the network designer (defender) acts first to locate hubs to route a set of flows through the network. The attacker (interdictor) acts next to interdict a subset of the located hubs in the designer’s network, followed again by the defender who routes the flows through the remaining hubs in the network. We model the defender’s problem as a trilevel optimization problem, wherein the attacker’s response is modeled as a bilevel hub interdiction problem. We study such a trilevel problem on three variants of hub location problems studied in the literature namely: p-hub median problem, p-hub center, and p-hub maximal covering problems. We present a cutting plane based exact method to solve the problem. The cutting plane method uses supervalid inequalities, which is obtained from the solution of the lower level interdiction problem. To solve the lower level hub interdiction problem efficiently, we propose a penalty-based reformulation of the problem. Using the reformulation, we present a branch-and-cut based exact approach to solve the problem efficiently. We conduct experiments to show the computational advantages of the above algorithm. To the best of our knowledge, the cutting plane approach proposed in this paper is among the first exact method to solve trilevel location–interdiction problems. Our computational results show interesting implications of incorporating interdiction risks in the hub location problem.

An artificial bee colony algorithm for the minimum edge-dilation K-center problem

An artificial bee colony algorithm for the minimum edge-dilation K-center problem

Facility location optimization with market segmentation for additive manufacturing material supplier network using quantum annealing

The determination of the material suppliers for additive manufacturing (AM) centers is a typical facility location problem (FLP). Current methods for solving the FLP in AM network design mainly focus on improving the cost, distance, and supply chain aspects and neglect market segmentation. However, AM feedstock suppliers need to understand the market size they serve so they can better prepare their inventory plans. Therefore, in this paper, we propose a new approach for the FLP with geographic market segmentation of AM material supplier networks, in which the facility locations and market segments of AM material suppliers are optimized in a multi-objective optimization problem. In order to confirm the feasibility of the proposed method, a simplified design case for a German state has been carried out.

A reinforcement learning guided hybrid evolutionary algorithm for the latency location routing problem

The latency location routing problem integrates the facility location problem and the multi-depot cumulative capacitated vehicle routing problem. This problem involves making simultaneous decisions about depot locations and vehicle routes to serve customers while aiming to minimize the sum of waiting (arriving) times for all customers. To address this computationally challenging problem, we propose a reinforcement learning guided hybrid evolutionary algorithm following the framework of the memetic algorithm. The proposed algorithm relies on a diversity-enhanced multi-parent edge assembly crossover to build promising offspring and a reinforcement learning guided variable neighborhood descent to determine the exploration order of multiple neighborhoods. Additionally, strategic oscillation is used to achieve a balanced exploration of both feasible and infeasible solutions. The competitiveness of the algorithm against state-of-the-art methods is demonstrated by experimental results on the three sets of 76 popular instances, including 51 improved best solutions (new upper bounds) for the 59 instances with unknown optima and equal best results for the remaining instances. We also conduct additional experiments to shed light on the key components of the algorithm.

Capacitated Mobile Facility Location Problem with Mobile Demand: Efficient Relief Aid Provision to En Route Refugees

As a humanity crisis, the tragedy of forced displacement entails relief aid distribution efforts among en route refugees to alleviate their migration hardships. This study aims to assist humanitarian organizations in cost-efficiently optimizing the logistics of capacitated mobile facilities utilized to deliver relief aid to transiting refugees in a multi-period setting. The problem is referred to as the Capacitated Mobile Facility Location Problem with Mobile Demands (CMFLP-MD). In CMFLP-MD, refugee groups follow specific paths, and meanwhile, they receive relief aid at least once every fixed number of consecutive periods, maintaining continuity of service. To this end, the overall costs associated with capacitated mobile facilities, including fixed, service provision, and relocation costs, are minimized. We formulate a mixed integer linear programming (MILP) model and propose two solution methods to solve this complex problem: an accelerated Benders decomposition approach as an exact solution method and a matheuristic algorithm that relies on an enhanced fix-and-optimize agenda. We evaluate our methodologies by designing realistic instances based on the Honduras migration crisis that commenced in 2018. Our numerical results reveal that the accelerated Benders decomposition excels MILP with a 46% run time improvement on average while acquiring solutions at least as good as the MILP across all instances. Moreover, our matheuristic acquires high-quality solutions with a 2.4% average gap compared to best-incumbents rapidly. An in-depth exploration of the solution properties underscores the robustness of our relief distribution plans under varying migration circumstances. Across several metrics, our sensitivity analyses also highlight the managerial advantages of implementing CMFLP-MD solutions.

Fairness in accessibility of public service facilities

Our paper studies a fair stochastic facility location problem with congestion under the max-min principle. We analyze the price of fairness in a two-location problem and develop a tractable optimization framework for the general multi-location setting. Evaluating the fair solution against the utilitarian solution on Buffalo's demographic data reveals that implementing a fair solution can substantially improve fairness measures (up to 98%) with relatively limited impact on the overall service quality.

Median and small parsimony problems on RNA trees.

Noncoding RNAs (ncRNAs) express their functions by adopting molecular structures. Specifically, RNA secondary structures serve as a relatively stable intermediate step before tertiary structures, offering a reliable signature of molecular function. Consequently, within an RNA functional family, secondary structures are generally more evolutionarily conserved than sequences. Conversely, homologous RNA families grouped within an RNA clan share ancestors but typically exhibit structural differences. Inferring the evolution of RNA structures within RNA families and clans is crucial for gaining insights into functional adaptations over time and providing clues about the Ancient RNA World Hypothesis. We introduce the median problem and the small parsimony problem for ncRNA families, where secondary structures are represented as leaf-labeled trees. We utilize the Robinson-Foulds (RF) tree distance, which corresponds to a specific edit distance between RNA trees, and a new metric called the Internal-Leafset (IL) distance. While the RF tree distance compares sets of leaves descending from internal nodes of two RNA trees, the IL distance compares the collection of leaf-children of internal nodes. The latter is better at capturing differences in structural elements of RNAs than the RF distance, which is more focused on base pairs. We also consider a more general tree edit distance that allows the mapping of base pairs that are not perfectly aligned. We study the theoretical complexity of the median problem and the small parsimony problem under the three distance metrics and various biologically relevant constraints, and we present polynomial-time maximum parsimony algorithms for solving some versions of the problems. Our algorithms are applied to ncRNA families from the RFAM database, illustrating their practical utility. https://github.com/bmarchand/rna\_small\_parsimony.

The Ordered Median Tree Location Problem

In this paper, we propose the Ordered Median Tree Location Problem (OMT). The OMT is a single-allocation facility location problem where p facilities must be placed on a network connected by a non-directed tree. The objective is to minimize the sum of the ordered weighted averaged allocation costs plus the sum of the costs of connecting the facilities in the tree. We present different MILP formulations for the OMT based on properties of the minimum spanning tree problem and the ordered median optimization. Given that ordered median location problems are rather difficult to solve we have improved the OMT solution performance by introducing covering variables in a valid reformulation plus developing two pre-processing phases to reduce the size of this formulations. In addition, we propose a Benders decomposition algorithm to approach the OMT. We establish an empirical comparison between these new formulations and we also provide enhancements that together with a proper formulation allow to solve medium size instances on general random graphs.

Robust vertex centdian facility location problem on tree networks

Robust vertex centdian facility location problem on tree networks

Solving Optimal Electric Vehicle Charger Deployment Problem

Electric vehicles (EVs) have already been acknowledged to be the most viable solution to the climate change that the entire globe has long been combating. Along the same line, it is a salient subject to expand the availability of EV charging infrastructure, which quintessentially necessitates the optimization of the charger’s locations. This paper proposes to formulate the optimal EV charger location problem into a facility location problem (FLP). As an effort to find an efficient method to solve the well-known nonpolynomial deterministic (NP) hard problem, we present a comparative quantification among several representative solving techniques. This paper features two comprehensive case studies representing regions with an average and a high density of EVs. As such, this paper shows that the proposed framework can lead to successful location optimization with adequate refinement of solving techniques.

Code and Data Repository for Ranking decomposition for the discrete ordered median problem

This repository provides the source code, data files and detailed results as reported in the article. The main folders are 'data', 'results', 'src' and 'test'. 'data': This folder includes discrete ordered median instances. 'results': This folder contains an excel file with detailed computational results. 'src': The source code. 'test': Tests.

Dispersed starting solutions in facility location: The case of the planar [formula omitted]-median problem

There are many planar multiple facilities location problems for which the optimal locations tend to be spread out. The most popular of these is the planar p-median problem. With this in mind, we propose several procedures to generate sparse configurations as starting solutions. The proposed procedures are easy to implement, and can be used as modules combined in different sequences within heuristics such as a recent trajectory-based procedure that we tested in this paper.The procedures are tested experimentally on a set of 24 large problem instances with up to 10,000 demand points and 100 facilities. We are able to demonstrate that the sparse starting solutions generated by the new procedures lead to significant improvements of final p-median solutions.

Solving Facility Location Problems via FastMap and Locality Sensitive Hashing

Solving Facility Location Problems via FastMap and Locality Sensitive Hashing

A Lagrangian relaxation approach for resource allocation problem with capacity constraints

AbstractThis study evaluates a capacitated facility location model enhanced with distance constraints for an emergency response problem, ensuring certain neighborhoods remain within an accessible range from facilities following a hurricane. The proposed model takes into account the capacity constraints for drones and vehicles. The model determines optimal locations for facilities and the distribution of supplies across the city. It also specifies which facilities should support the needs of each neighborhood and decides on the appropriate mode of transportation—ground vehicles if possible, or drones if roadways are obstructed. To solve the problem, a Lagrangian relaxation technique is employed, relaxing the constraints related to facility capacity and distance. The numerical results confirm the quality and efficiency of the solutions. The findings indicate that ground transportation is more frequently utilized than drones at each operational facility. A comprehensive set of sensitivity analyses is conducted to examine the impact of various variables and parameters on the solution.

Resolving Charging Station Placement Issues for Electric Vehicles: Hybrid Optimization-Assisted Multi-Objective Framework

As new energy technology has advanced to a respectable level, Electric Vehicles (EV) have gained recognition from people all over the world and have become increasingly popular in a number of nations. Effective charging station placement is critical for the rapid growth of electric vehicles, as it is vital to provide convenience for electric vehicles while also ensuring the effectiveness of traffic networks. The location of an EV charging station is simply an application situation for the facility location problem. Traditional works often pay attention to the mileage concerns of electric vehicle customers while ignoring their competitive and strategic charging tactics. This research aims to frame the allocation of charging stations issue as a multi-faceted venture by assessing the aspects economically along with the characteristics of the power grid like “cost, Voltage stability, Reliability, and Power loss (VRP) index, waiting time, and accessibility index”. Further, we proposed a hybrid artificial intelligence to resolve the allocation issue. The proposed new Artificial Intelligence (AI) based algorithm is termed as Hybridized Salp and Harris Algorithm (HS-HA). The effectiveness and scalability of the presented algorithm for the charging station placement issue are also assessed through different plans in the IEEE 33-bus distribution systems. Finally, the efficiency of the presented approach is proved over the existing approaches with respect to various measures.

A congested facility location problem with strategic customers

This study addresses a user-equilibrium congested facility location problem with delay-, accessibility-, and price-sensitive customers. A profit-maximizing service provider first makes location, service rate, and pricing decisions, and then strategic customers decide which facilities to patronize (if any). By incorporating the customers’ choice behavior as a set of equilibrium constraints into the service provider’s decision problem, the problem is modeled as a mixed integer non-linear program which is then reformulated as a mixed integer second-order cone program. The proposed model, tested on several standard instances, is shown to be efficiently solvable using commercial software packages.

Distributionally robust disaster relief planning under the Wasserstein set

We study a two-stage natural disaster management problem modeled as a stochastic program, where the first stage consists of a facility location problem, deciding where to open facilities and pre-allocate resources such as medical and food kits, and the second stage is a fixed-charge transportation problem, routing resources to affected areas after observing a disaster. Our model has binary variables present in both stages. Due to the lack of data, classical stochastic programming approaches may be ill-suited, and we propose a two-stage distributionally robust formulation with a Wasserstein ambiguity set, where we consider distributions consistent with historical data and a tunable parameter to control the level of risk aversion. We develop a tailored column-and-constraint generation (CCG) algorithm to solve an extensive reformulation, where scenarios are iteratively generated. We handle the presence of binary variables in the second stage by leveraging the structure of our support set and second-stage problem, and provide conditions under which the optimal value of the latter is concave with respect to the intensity of the disaster, leading to an efficient scenario generation procedure. We also show that our results extend to the case where the second stage is a fixed-charge network flow problem. We perform extensive computational experiments demonstrating the computational advantage of our method over classical CCG implementations on synthetic instances, and illustrate the benefits of our approach on a popular case study from the literature of hurricane threats on the Gulf of Mexico in the United States.

Polyakov loop models in the large N limit: Correlation function and screening masses

We explore the ’t Hooft-Veneziano limit of the Polyakov loop models at finite baryon chemical potential. Using methods developed by us earlier we calculate the two- and N-point correlation functions of the Polyakov loops. This gives a possibility to compute the various potentials in the confinement phase and to derive the screening masses outside the confinement region. In particular, we establish the existence of complex masses and an oscillating decay of correlations in a certain range of parameters. Furthermore, it is shown that the calculation of the N-point correlation function in the confinement phase reduces to the geometric median problem. This leads to a large N analog of the Y law for the baryon potential. Published by the American Physical Society 2024