The Rank Pricing Problem (RPP) is a challenging bilevel optimization problem with binary variables whose objective is to determine the optimal pricing strategy for a set of products to maximize the total benefit, given that customer preferences influence the price for each product. Traditional methods for solving RPP are based on exact approaches which may be computationally expensive. In contrast, this paper presents a novel heuristic approach that takes advantage of the structure of the problem to obtain good solutions. The proposed approach consists of two phases. Firstly, a standard heuristic is applied to get a pricing strategy. In our case, we choose to use the Variable Neighborhood Search (VNS), and the genetic algorithm. Both methodologies are very popular for their effectiveness in solving combinatorial optimization problems. The solution obtained after running these algorithms is improved in a second phase, where four different local searches are applied. Such local searches use the information of the RPP to get better solutions, that is, there is no need to solve new optimization problems. Even though our methodology does not have optimality guarantees, our computational experiments show that it outperforms Mixed Integer Program solvers regarding solution quality and computational burden.

In this paper, we tackle the resolution of chance-constrained problems reformulated via sample average approximation. The resulting data-driven deterministic reformulation takes the form of a large-scale mixed-integer program (MIP) cursed with Big-Ms. We introduce an exact resolution method for the MIP that combines the addition of a set of valid inequalities to tighten the linear relaxation bound with coefficient strengthening and constraint screening algorithms to improve its Big-Ms and considerably reduce its size. The proposed valid inequalities are based on the notion of k-envelopes and can be computed off-line using polynomial-time algorithms and added to the MIP program all at once. Furthermore, they are equally useful to boost the strengthening of the Big-Ms and the screening rate of superfluous constraints. We apply our procedures to a probabilistically constrained version of the DC optimal power flow problem with uncertain demand. The chance constraint requires that the probability of violating any of the power system’s constraints be lower than some positive threshold. In a series of numerical experiments that involve five power systems of different size, we show the efficiency of the proposed methodology and compare it with some of the best performing convex inner approximations currently available in the literature. History: Accepted by Andrea Lodi, Area Editor for Design & Analysis of Algorithms – Discrete. Funding: This work was supported in part by the European Research Council under the EU Horizon 2020 research and innovation program [Grant 755705], in part by the Spanish Ministry of Science and Innovation [Grant AEI/10.13039/501100011033] through project PID2020-115460GB-I00, and in part by the Junta de Andalucía and the European Regional Development Fund through the research project P20_00153. Á. Porras is also financially supported by the Spanish Ministry of Science, Innovation and Universities through the University Teacher Training Program with fellowship number FPU19/03053. Supplemental Material: The online supplement is available at https://doi.org/10.1287/ijoc.2022.0302 .

Mixed Integer Linear Programs (MILP) are well known to be NP-hard (Non-deterministic Polynomial-time hard) problems in general. Even though pure optimization-based methods, such as constraint generation, are guaranteed to provide an optimal solution if enough time is given, their use in online applications remains a great challenge due to their usual excessive time requirements. To alleviate their computational burden, some machine learning techniques (ML) have been proposed in the literature, using the information provided by previously solved MILP instances. Unfortunately, these techniques report a non-negligible percentage of infeasible or suboptimal instances. By linking mathematical optimization and machine learning, this paper proposes a novel approach that speeds up the traditional constraint generation method, preserving feasibility and optimality guarantees. In particular, we first identify offline the so-called invariant constraint set of past MILP instances. We then train (also offline) a machine learning method to learn an invariant constraint set as a function of the problem parameters of each instance. Next, we predict online an invariant constraint set of the new unseen MILP application and use it to initialize the constraint generation method. This warm-started strategy significantly reduces the number of iterations to reach optimality, and therefore, the computational burden to solve online each MILP problem is significantly reduced. Very importantly, all the feasibility and optimality theoretical guarantees of the traditional constraint generation method are inherited by our proposed methodology. The computational performance of the proposed approach is quantified through synthetic and real-life MILP applications.