Mathematical models and search algorithms for the capacitated $p$-center problem

TL;DR

This paper introduces advanced search algorithms and mathematical formulations to efficiently solve the capacitated p-center problem, achieving optimal solutions for large benchmark instances and improving computational performance.

Contribution

It develops novel search strategies and valid inequalities, including a balanced search approach, to enhance solving capacity for large-scale capacitated p-center problems.

Findings

All benchmark instances up to 402 vertices solved optimally.

Feasible solutions found within 10 minutes for instances with up to 3038 vertices.

Proposed methods outperform previous approaches in computational efficiency.

Abstract

The capacitated p-center problem requires to select p facilities from a set of candidates to service a number of customers, subject to facility capacity constraints, with the aim of minimizing the maximum distance between a customer and its associated facility. The problem is well known in the field of facility location, because of the many applications that it can model. In this paper, we solve it by means of search algorithms that iteratively seek the optimal distance by solving tailored subproblems. We present different mathematical formulations for the subproblems and improve them by means of several valid inequalities, including an effective one based on a 0-1 disjunction and the solution of subset sum problems. We also develop an alternative search strategy that finds a balance between the traditional sequential search and binary search. This strategy limits the number of feasible subproblems to be solved and, at the same time, avoids large overestimates of the solution value, which are detrimental for the search. We evaluate the proposed techniques by means of extensive computational experiments on benchmark instances from the literature and new larger test sets. All instances from the literature with up to 402 vertices and integer distances are solved to proven optimality, including 13 open cases, and feasible solutions are found in 10 minutes for instances with up to 3038 vertices.