ColumnGeneration-Exercise-Set

.pdf

School

University of Houston *

*We aren’t endorsed by this school

Course

1301

Subject

Industrial Engineering

Date

Oct 30, 2023

Type

pdf

Pages

6

Uploaded by BaronKnowledge37615

Report
Mathematical Programming Exercise set 2 - Column generation Exercise 1 In this exercise we study the facility location problem . Consider m facilities and n customers. The cost to service customer j ∈ { 1 , 2 , . . . , n } from facility i ∈ { 1 , 2 , . . . , m } is c ij . The cost for opening facility i is equal to f i , for i ∈ { 1 , 2 , . . . , m } . At most W facilities may be opened and each customer must be serviced by a facility. The problem is to decide (1) which facilities should be opened and (2) which open facility should service which customer such that the total costs are minimized. We define a facility plan as a pair ( i ( p ) , S ( p )), where i ( p ) ∈ { 1 , . . . , m } is a facility that is opened, and the subset of customers S ( p ) ⊆ { 1 , . . . , n } is serviced from i ( p ). Let P be the set of all facility plans. a. Formulate this problem as a binary programming problem. Make use of the variables x p , indicating whether facility plan p P is used. Clearly define the corresponding objective coefficients, and clearly explain your formulation. b. Suppose that we want to solve the LP relaxation of the binary programming formulation using column generation. For given values of the dual variables, describe the pricing problem and an algorithm to solve it. c. Suppose we perform column management as follows. In each iteration of the column gen- eration algorithm, we remove all columns that have not been part of an optimal solution to the RMP for at least 10 consecutive iterations. Explain why this could speed up the column generation algorithm. Also explain why this could slow down the column generation algorithm. ( The take-away is that testing is required to see if it column management helps. ) Answer of Exercise 1 a. The cost coefficient corresponding to x p is C p = f i ( p ) + j S ( p ) c ij . Let the parameter a jp be 1 if j S ( p ) and 0 otherwise. Consider the following binary programming formulation. min summationdisplay p P C p x p (1) summationdisplay p P a jp x p = 1 j ∈ { 1 , . . . , n } (2) summationdisplay p P x p W (3) x p ∈ { 0 , 1 } ∀ p P (4) 1 Downloaded by Taghi Khaniyev (tagikhaniyev@gmail.com) lOMoARcPSD|24058970
The objective (1) is to minimize the total costs. Constraints (2) ensure that every customer is serviced. Constraint (3) ensure that not more than W facilities are opened. Finally, the binarity conditions are given by (4). b. Let λ j , for j ∈ { 1 , . . . , n } , be the dual variables corresponding to (2). Furthermore, let μ be the dual variable corresponding to (3). The reduced cost of a variable x p is RC ( x p ) = C p - n summationdisplay j =1 a jp λ j - μ = f i ( p ) - μ + summationdisplay j S ( p ) ( c i ( p ) j - λ j ) The pricing problem is to find a facility plan that minimizes the reduced cost. This pricing problem is not very difficult to solve, and a pricing algorithm is the following. Step1: For all facilities i ∈ { 0 , . . . , m } , we select the set S of customers for which c ij - λ j < 0. Step 2: For each facility i and corresponding subset S found in Step 1, compute f i - μ + j S ( c ij - λ j ). The lowest cost is the optimal solution value, and corresponds to the optimal solution of the pricing problem. c. Having less columns in the RMP speeds up solving the RMP. On the other hand, keeping track of columns to remove, and removing them takes (a tiny bit of) computation time. Moreover, columns that are removed might be useful in a later iteration of the column generation algorithm, and then need to be generated again. This would slow down the algorithm. Exercise 2 In this exercise we study the bin packing problem . Consider n items each with known volume, indicated by v i for item i ∈ { 1 , . . . , n } , which have to be packed into a number of bins, each having the same volume V . Items packed into the same bin should have a total volume not exceeding the volume of the bin and the objective is to pack all items using a minimal number of bins. We define a packing as a subset of items that do not exceed the bin capacity. Let P be the set of all packings. a. Formulate this problem as a binary programming problem. Make use of the packing vari- ables x p , indicating whether packing p P is used. Clearly define the corresponding objective coefficients, and clearly explain your formulation. b. Suppose that we want to solve the LP relaxation of the binary programming formulation using column generation. For given values of the dual variables, provide the pricing problem. Do you recognize the pricing problem? c. Explain why it is not always necessary to solve the pricing problem to optimality when generating additional columns. d. Now suppose we use a branch-and-price algorithm to solve the bin packing problem. Provide a disadvantage of using a branching rule that branches on a fractional packing variable. Answer of Exercise 2 2 Downloaded by Taghi Khaniyev (tagikhaniyev@gmail.com) lOMoARcPSD|24058970
a. The cost coefficient corresponding to x p is 1. Let the parameter a ip be 1 if item i is included in packing p and 0 otherwise. Consider the following binary programming formulation. min summationdisplay p P x p (5) summationdisplay p P a ip x p = 1 i ∈ { 1 , . . . , n } (6) x p ∈ { 0 , 1 } ∀ p P (7) The objective (5) is to minimize the total costs. Constraints (6) ensure that every item is included in exactly one bin. Finally, the binarity conditions are given by (7). b. Let λ i , for i ∈ { 1 , . . . , n } , be the dual variables corresponding to (6). The reduced cost of a variable x p is RC ( x p ) = 1 - n summationdisplay i =1 a ip λ i The pricing problem is to find a packing that minimizes the reduced cost. The pricing problem is equivalent to the following problem. max n summationdisplay i =1 λ i y i n summationdisplay i =1 v i y i V y i ∈ { 0 , 1 } ∀ i ∈ { 0 , . . . , n } We can recognize the pricing problem as a knapsack problem, for which efficient special purpose algorithms exist. Note that if the solution value of the above knapsack problem is strictly greater than 1, a negative reduced cost variable is found which could then be added to the RMP. Otherwise, no negative reduced cost variable exists, and the column generation algorithm terminates because the optimal solution to the LP relaxation is found. c. If a negative reduced cost variable is identified, it is a valid non-basic variable to add to the current basis (think of the simplex method). That is, it can be added to the RMP. For correctness of the column generation algorithm, it is not required that the most negative reduced cost variable is added. Note, however, that if the pricing problem is not solved to optimality and no negative reduced cost variable is identified, it might still be the case that at least one such variable exists. Hence, in this case we cannot yet conclude that the optimal solution to the LP has been found. d. If we branch on a packing variable directly, this affects the pricing problem. For instance, suppose in one node of the branching tree we set the variable x p = 0. Now, we use column generation to solve the LP relaxation corresponding to that branching node. Here, in the pricing problem we are not allowed to generate the packing p . We could formulate this new pricing problem as follows. Denote by | p | the amount of items in the packing p , we then have 3 Downloaded by Taghi Khaniyev (tagikhaniyev@gmail.com) lOMoARcPSD|24058970
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help

Browse Popular Homework Q&A

Q: The following initial rate data were obtained for 2 NO(g) + O₂(g) → 2 NO₂(g). [NO] [0₂]0…
Q: As part of a statistics project, Charlie would like to collect data on household size in his city.…
Q: Suppose L(x) = -5x + 1 is the local linearization to h(x) at x=2. What is h(2) ? not enough info to…
Q: A group of people were asked if they had run a red light in the last year. 479 responded "yes", and…
Q: Draw a Fischer projection of: meso 2,3-diphenylbutane 8. Which of these two would you expect to be…
Q: (2x, 2y, 2x + 2z). Use Stokes' theorem to evaluate the integral of F around the curve consisting of…
Q: A mummy was recently discovered. It now contains 72.7​% of its original Carbon 14. About how old is…
Q: Suppose a sample of 16 light trucks is randomly selected off the assembly line. The trucks are…
Q: Using your answer from the previous question, find the p-value.
Q: H3C. C H2 H3C | CH HU H2C H2 C C OI CH2 CH3
Q: Graph the following inequality. y> -5 Use the graphing tool to graph the inequality. Click to…
Q: Which of the following is a tropic (stimulating) hormone? Group of answer choices Thyroxine…
Q: 10. Kristine has $5,000 in an account today that pays 3% interest rate. Two years from now she…
Q: using oracle
Q: Complete the following non-metric conversions using the table of conversions in the slides. Remember…
Q: Consider a sample of carbon tetrafluoride that contains 1.59 × 1024 atoms of fluorine. Using the…
Q: When a nextLine() command is followed by a nextInt() or nextDouble() instruction, what occurs in the…
Q: = A sample of semiconductor is doped with ND 1013 cm³. The intrinsic carrier density is 1013 cm-³.…
Q: Let S : RR be defined by S(x1, x2, 3,...) = (T2, T3, T4,...) (a) Using the definition of linear…
Q: Find the differential of the function f(x, y) = 3xe2y at the point (5, 3) using Ar = 0.25 and Ay = (…
Q: What is O undefined 0 1 O-1 O d de O pi/2 (csc0)| 0=π/2
Q: Elizabeth, a U.S. American student, and Ahmad, an international student from Turkey, find that…