Elsevier

Operations Research Letters

Volume 45, Issue 5, September 2017, Pages 513-518
Operations Research Letters

Lifting of probabilistic cover inequalities

https://doi.org/10.1016/j.orl.2017.08.006Get rights and content

Abstract

We consider a chance-constrained binary knapsack problem where weights of items are independent and normally distributed. Probabilistic cover inequalities can be defined for the problem. The lifting problem for probabilistic cover inequalities is NP-hard. We propose a polynomial time approximate lifting method for probabilistic cover inequalities based on the robust optimization approach. We present computational experiments on multidimensional chance-constrained knapsack problems. The results show that our lifting method reduces the computation time substantially.

Introduction

We consider a binary knapsack problem with uncertain weights of items. The binary knapsack problem is to select items to maximize the sum of profits while satisfying knapsack capacity constraint. The binary knapsack problem without uncertainty is known to be NP-hard [13]. Many practical problems can be formulated using the binary knapsack constraint, and studies on the knapsack solution set can be applied to any mixed integer programming problems containing the knapsack constraint. Therefore the binary knapsack problem has been extensively studied [13], [17], and much research is focused on the polyhedral aspects of the binary knapsack solution set.

In applications of optimization, an important issue is how to deal with data uncertainty. Representative approaches to deal with data uncertainty are stochastic programming and robust optimization. Stochastic optimization considers a probabilistic distribution of uncertain data [6], and chance-constrained programming [8] is one of the useful approaches of stochastic programming. It finds the best solution satisfying constraints with probability at least a given threshold. The chance-constrained knapsack problem has been studied by many authors. Goel et al. [9] proposed a polynomial time approximation scheme (PTAS) where weight of each item has a Poisson or exponential distribution. Kleinberg et al. [14] assumed that each weight has a Bernoulli distribution and presented an approximation algorithm. Klopfenstein et al. [15] proposed an approximation algorithm using robust optimization techniques. Goyal et al. [10] presented a PTAS using a parametric LP reformulation when item sizes are independent and normally distributed. Han et al. [12] proposed an efficient pseudopolynomial algorithm based on the robust optimization approach for finding a good upper bound on the optimal value when weight of each item has a normal distribution independent of the other items. This paper also focuses on the chance-constrained knapsack problem where weights are independent and normally distributed.

Cover inequalities are well-known valid inequalities for the binary knapsack solution set. Commercial optimization softwares such as CPLEX and Xpress-MP use cover inequalities in their branch-and-cut algorithms. Lifting is a powerful technique to strengthen cover inequalities [4], [21]. Cover inequalities also can be defined for the knapsack problems considering data uncertainty. Klopfenstein et al. [16] defined robust cover inequalities for the robust knapsack problem that is formulated using the uncertainty set of Bertsimas et al. [5]. Song et al. [20] considered the chance-constrained binary packing problem and proposed a problem formulation using probabilistic covers. They focused on the case with finite number of scenarios and considered an approximate lifting method. The chance-constrained knapsack problem is a special case of the chance-constrained packing problem. Atamtürk et al. [3] described cover inequalities for the submodular knapsack set and proposed a polynomial time approximate lifting algorithm using parametric LP. Note that the chance-constrained knapsack problem with independent and normally distributed weights is one of submodular knapsack problems. Atamtürk et al. [1] investigated separation and extension of cover inequalities for the conic quadratic knapsack constraint with generalized upper bound (GUB) constraints.

In this paper we propose a polynomial time approximate lifting algorithm for cover inequalities for the chance-constrained knapsack problem. We assume that weights of items are independent and normally distributed. We formulate the subproblem for lifting of probabilistic cover inequalities as a chance-constrained knapsack problem. Then we propose an approximate lifting algorithm using the efficient approximation algorithm of Han et al. [12].

Section snippets

Chance-constrained knapsack problem

We consider a chance-constrained knapsack problem. Let N={1,2,,n} be a set of items. Each item jN has profit pj and uncertain weight ãj. The objective is to find a subset of items with maximum sum of profits such that the probability of satisfying the knapsack capacity b is at least ρ. We assume that ρ is greater than or equal to 0.5. The variable xj=1 if item j is selected and xj=0 otherwise. Then the chance-constrained knapsack problem can be formulated as follows. maxjNpjxjs.t.PjNãjxj

Probabilistic cover

A subset of N is called a probabilistic cover if the items of the probabilistic cover cannot be included in the knapsack simultaneously [3], [20].

Definition 1

CN is a probabilistic cover if jCāj+Φ1(ρ)jCσj2>b.

If C is a probabilistic cover, then the probabilistic cover inequality jCxj|C|1 is valid for the feasible solution set of (2). Probabilistic cover inequalities can be strengthened using sequential lifting of coefficients of items in N

Computational results

In this section, we report the computational results on the chance-constrained knapsack problem (2) and the multidimensional chance-constrained knapsack problem. All tests were performed on a computer Intel(R) Core(TM) i5-3570K processor 3.40 GHz with 8GB RAM. We implemented problems in Java using IBM ILOG CPLEX 12.7. When we use CPLEX to solve a quadratically constrained mixed integer program, linear program (LP) relaxation based branch-and-cut algorithm was used. We tried both quadratically

Acknowledgment

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government(MSIP) (No. 2016R1A2B4013590).

References (21)

There are more references available in the full text version of this article.
View full text