Department of Mathematics Oleg Burdakov of 30 October Consider the following linear programming problem (LP):

Similar documents
56:272 Integer Programming & Network Flows Final Examination -- December 14, 1998

Introduction to Mathematical Programming IE496. Final Review. Dr. Ted Ralphs

Linear Programming. Course review MS-E2140. v. 1.1

MATHEMATICS II: COLLECTION OF EXERCISES AND PROBLEMS

ME 391Q Network Flow Programming

Programming, numerics and optimization

Section Notes 5. Review of Linear Programming. Applied Math / Engineering Sciences 121. Week of October 15, 2017

Ryerson Polytechnic University Department of Mathematics, Physics, and Computer Science Final Examinations, April, 2003

Integer Programming. Xi Chen. Department of Management Science and Engineering International Business School Beijing Foreign Studies University

Solutions for Operations Research Final Exam

SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING

Some Advanced Topics in Linear Programming

Optimization Methods. Final Examination. 1. There are 5 problems each w i t h 20 p o i n ts for a maximum of 100 points.

Linear Programming. Linear programming provides methods for allocating limited resources among competing activities in an optimal way.

Applied Lagrange Duality for Constrained Optimization

Constrained Optimization COS 323

56:272 Integer Programming & Network Flows Final Exam -- December 16, 1997

BCN Decision and Risk Analysis. Syed M. Ahmed, Ph.D.

3 INTEGER LINEAR PROGRAMMING

DM515 Spring 2011 Weekly Note 7

Read: H&L chapters 1-6

LP-Modelling. dr.ir. C.A.J. Hurkens Technische Universiteit Eindhoven. January 30, 2008

Linear Optimization. Andongwisye John. November 17, Linkoping University. Andongwisye John (Linkoping University) November 17, / 25

AM 121: Intro to Optimization Models and Methods Fall 2017

Fundamentals of Integer Programming

Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras

Part 4. Decomposition Algorithms Dantzig-Wolf Decomposition Algorithm

Lecture 12: Feasible direction methods

SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING

Unconstrained Optimization Principles of Unconstrained Optimization Search Methods

VARIANTS OF THE SIMPLEX METHOD

5.3 Cutting plane methods and Gomory fractional cuts

Advanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs

Section Notes 4. Duality, Sensitivity, and the Dual Simplex Algorithm. Applied Math / Engineering Sciences 121. Week of October 8, 2018

COLUMN GENERATION IN LINEAR PROGRAMMING

Module 1 Lecture Notes 2. Optimization Problem and Model Formulation

MVE165/MMG630, Applied Optimization Lecture 8 Integer linear programming algorithms. Ann-Brith Strömberg

Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras

MULTIMEDIA UNIVERSITY FACULTY OF ENGINEERING PEM2046 ENGINEERING MATHEMATICS IV TUTORIAL

February 19, Integer programming. Outline. Problem formulation. Branch-andbound

OPTIMIZATION METHODS

Chapter 15 Introduction to Linear Programming

Introduction to Mathematical Programming IE406. Lecture 20. Dr. Ted Ralphs

Linear Programming Problems

TMA946/MAN280 APPLIED OPTIMIZATION. Exam instructions

Unit.9 Integer Programming

ACO Comprehensive Exam October 12 and 13, Computability, Complexity and Algorithms

Demo 1: KKT conditions with inequality constraints

The Ascendance of the Dual Simplex Method: A Geometric View

David G. Luenberger Yinyu Ye. Linear and Nonlinear. Programming. Fourth Edition. ö Springer

Linear Programming. Larry Blume. Cornell University & The Santa Fe Institute & IHS

Advanced Operations Research Techniques IE316. Quiz 2 Review. Dr. Ted Ralphs

4.1 Graphical solution of a linear program and standard form

MVE165/MMG631 Linear and integer optimization with applications Lecture 9 Discrete optimization: theory and algorithms

16.410/413 Principles of Autonomy and Decision Making

Linear Programming Terminology

George B. Dantzig Mukund N. Thapa. Linear Programming. 1: Introduction. With 87 Illustrations. Springer

Gate Sizing by Lagrangian Relaxation Revisited

INEN 420 Final Review

DM545 Linear and Integer Programming. Lecture 2. The Simplex Method. Marco Chiarandini

Lecture 25 Nonlinear Programming. November 9, 2009

Final Exam Spring 2003

ISE203 Optimization 1 Linear Models. Dr. Arslan Örnek Chapter 4 Solving LP problems: The Simplex Method SIMPLEX

Linear programming II João Carlos Lourenço

NOTATION AND TERMINOLOGY

Chapter II. Linear Programming

Outline. CS38 Introduction to Algorithms. Linear programming 5/21/2014. Linear programming. Lecture 15 May 20, 2014

Graphs and Network Flows IE411. Lecture 20. Dr. Ted Ralphs

MLR Institute of Technology

Optimization Methods in Management Science

Lagrangean relaxation - exercises

Mathematical and Algorithmic Foundations Linear Programming and Matchings

INTRODUCTION TO LINEAR AND NONLINEAR PROGRAMMING

Artificial Intelligence

Programs. Introduction

Selected Topics in Column Generation

4.1 The original problem and the optimal tableau

Integer Programming Theory

GENERAL ASSIGNMENT PROBLEM via Branch and Price JOHN AND LEI

Advanced Algorithms Linear Programming

Simulation. Lecture O1 Optimization: Linear Programming. Saeed Bastani April 2016

Generalized Network Flow Programming

5 Day 5: Maxima and minima for n variables.

CS599: Convex and Combinatorial Optimization Fall 2013 Lecture 14: Combinatorial Problems as Linear Programs I. Instructor: Shaddin Dughmi

Combinatorial Optimization

Computational Methods. Constrained Optimization

Discrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity

CS 473: Algorithms. Ruta Mehta. Spring University of Illinois, Urbana-Champaign. Ruta (UIUC) CS473 1 Spring / 36

Department of Electrical and Computer Engineering

Civil Engineering Systems Analysis Lecture XIV. Instructor: Prof. Naveen Eluru Department of Civil Engineering and Applied Mechanics

Tribhuvan University Institute Of Science and Technology Tribhuvan University Institute of Science and Technology

Local and Global Minimum

OR 674 DYNAMIC PROGRAMMING Rajesh Ganesan, Associate Professor Systems Engineering and Operations Research George Mason University

Optimization III: Constrained Optimization

Kernel Methods & Support Vector Machines

The Simplex Algorithm

Lagrangean Methods bounding through penalty adjustment

1. Lecture notes on bipartite matching February 4th,

Lecture 9: Linear Programming

Algorithms for Integer Programming

Transcription:

Linköping University Optimization TAOP3(0) Department of Mathematics Examination Oleg Burdakov of 30 October 03 Assignment Consider the following linear programming problem (LP): max z = x + x s.t. x x x + x + x 3 = x, x, x 3 0 (p) a) Solve the problem with the simplex method starting from the point for which x and x are nonbasic variables. What are the optimal values of x,x,x 3 and z? Is the optimal basic feasible solution degenerate? (p) b) Is the optimal solution obtained in a) unique? If not, find ALL alternative optimal solutions. (p) c) Formulate the LP-dual to the problem (LP). (p) d) Use the complementarity conditions to check if the point x = (3, 0, ) T is optimal. (p) e) How the optimal value of z would change (increase? decrease? remain the same?) if the right-hand side of the first constraint was increased a little? (Apply the sensitivity analysis) (p) f) Suppose a new variable x 4 is introduced in (LP) with the objective function coefficient c 4 = and constraint column a 4 = ( 3, ) T. Will the objective function value increase, decrease or remain the same? (Apply the sensitivity analysis.) Assignment (p) a) For some integer programming problem, its LP relaxation gives the following optimal simplex tableau bas. z x x x 3 x 4 r.h.s. z 0 0-0. -. -4 x 0 0 0.7-3. x 0 0-0.7.. Is this solution to LP relaxation optimal to the original integer programming problem? Find a Gomory cut for this solution. (p) b) Solve the following integer programming problem using Branch & Bound. max x +x s.t. x +x 6 6x +x 30 x,x 0 and integer

LP relaxation may be solved graphically. If the both components are fractional, branch on x. Choose the branch first. Follow the last-in-first-out rule. Assignment 3 Consider the function f(x) = x x x +x 4 +x 3 +x. (p) a) Determine whether d = (0, ) T is a descent direction for the point x = (7, ) T. (p) b) Consider the problem of minimizing f(x). Make one iteration of the steepest descent method starting from the same point x = (7, ) T. (3p) c) Find all stationary points of this function, and determine whether they are local minimizers and maximizers. Assignment 4 Consider the following nonlinear programming problem (P) min x +x s.t. (x ) +x 4 x,x 0 (p) a) Is (P) a convex programming problem? (p) b) Solve (P) graphically. Find all local and global minima. (p) c) Write the Karush-Kuhn-Tucker (KKT) conditions for the given problem (P). The conditions must be written for the specific functions that define (P). (p) d) Check if the KKT conditions hold at the point x = (0 0) T. (p) e) Find a KKT point for which x > 0 and x > 0. Show algebraically (not geometrically) that the KKT conditions hold for this point. (p) f) The penalty method reduces(p) to an unconstrained minimization problem min F(x, c), where c is a penalty parameter. Compute F(x,) at the point x = ( ) T.

Assignment There are five sub-networks of the form where i =,,3,4,. In each of them, it is possible to send a flow from its left side to its right side (l r), and from its right side to its left side (r l). The maximal l r and r l flows are given for each sub-network N i by the table: N i l r r l N 3 N 3 N 3 4 N 4 N 3 Consider the following network composed of the given five sub-networks and four nodes Each edge connecting nodes and sub-networks has unlimited capacity for sending flow in any direction. Consider the problem of sending maximum flow f from node to node 4. (p) a) Reformulate this problem as a standard maximum flow problem (construct a new directed graph, indicate upper and low bounds for each edge). (3p) b) Use the max-flow algorithm for finding the maximum flow f in the new directed graph. (p) c) Find a minimum cut in the new directed graph. 3

Assignment 6 Consider the network 4 6 4 3 3 (p) a) Find the shortest path from node to node. (p) b) What is the optimal solution to the linear programming formulation of the shortest path problem solved in a)? (4p) c) An owner of a network of hotels located in seven cities wishes to arrange a quick check of the air-conditioning systems in all his hotels. He is going to send workers to do this job. They start in city s. To travel from one city to another, they use hotel s cars (say, with unlimited number of sits). In the graph given below, each arc (i,j) is labeled with the cost associated with one-car traveling between city i and city j. 8 6 4 3 s 3 6 4 4 6 6 The owner does not pay for their return trip (for instance, because each of the workers is offered, as a bonus, to spend a vacation in the last hotel that this worker will check, but they must pay themselves for returning back after the vacation). The owner is interested in minimizing the total traveling cost that he pays. Use a proper algorithm for finding the minimal cost. What is the optimal traveling strategy? (For each city, excluding the terminal ones, indicate the city or the cities to go next and the number of workers in each car.) What is the minimal number of workers that ensures the minimal traveling cost? 7 7 8 Assignment 7 ( p.) Use the dynamic programming to solve the problem max 3x +4x +6x 3 s.t. 4x +x +7x 3 9 x,x,x 3 {0,} 4

ANSWERS Assignment a) The simplex method starts in x (0) = (0, 0,) T and generates the points x () = (0,,3) T, x () = (,,0) T. The point x () is optimal with z =. This optimal solution is not degenerate, because all the basic variables are strictly positive. b) The optimal solution is not unique. The alternative solution is x (3) = (, 0,0) T. All optimal solutions: x(λ) = λ +( λ) 0, 0 λ. 0 0 c) The dual problem: min w = v + v s.t. v + v v + v v 0 v 0, v is free d) The complementary slackness: () v (x x +) = 0 () x (v +v ) = 0 (3) x ( v +v ) = 0 (4) x 3 v = 0 For x = (3, 0, ) T, this implies: () v = 0 () v +v = 0 (4) v = 0 Obviously, there are no values of v and v which would satisfy these conditions. Thus, the complementary slackness is not satisfied for any v and v. This means that the point x = (3, 0, ) T is not optimal. e) The value of z would not change because the dual price of the st constraint is equal to zero (see the optimal simplex tableau). f) For any of the two optimal solutions x () or x (3), the complementary slackness gives the optimal values of the dual variables: v = 0, v =. Then c 4 = c 4 v Ta 4 = ( 0 )( ) 3 = 3. Since c 4 < 0, the optimal value z will not be improved (it will remain the same). Assignment a) The solution to the LP relaxation problem is not optimal to the original integer programming problem, because the value of x is not integer. The corresponding row of the simplex tableau gives the following equivalent Gomory cuts: x x 3 +x 4 0 and 0. 0.x 3 0.x 4 0.

b) P0: x LP = (7/, 9/) T, z LP =.3 UBD =, which is an upper bound for the original problem. P = P0 & (x 3): x LP = (3, ) T, z LP =. Do not continue branching further from this node, because x LP is integer. LBD =. Stop, because LBD = UBD =, i.e. it is impossible to improve the currently available candidate solution. The optimal solution: x = (3, ) T, z =. Assignment 3 a) d = (0, ) T is not a descent direction for the point x = (7, ) T, because d T f( x) = 0. b) x = x 0 +t 0 d 0 = (, ) T, where x 0 = (7, ) T, d 0 = (, 0) T and t 0 = /. c) The function has three stationary points: x = (0, 0) T, x = ( 0., 0.) T and x 3 = (, ) T. The points x and x 3 are local minimizers. The point x is neither local minimizer, nor local maximizer. Assignment 4 a) It is not a convex programming problem, because the feasible set is not convex. To show this we consider the two points a = (0, 0) T and b = (4, 0) T. These points are feasible, but λa+( λ)b is infeasible for λ = / (the st constraint is violated). b) Point a is a global minimizer, point b is a local minimizer. c) The ( KKT ) conditions: ( ) ( ) ( ) (x ) 0 +λ +λ +λ 0 3 = 0, λ 0, λ 0, λ 3 0, λ ((x ) +x 4) = 0, λ ( x ) = 0, λ 3 ( x ) = 0, (x ) +x 4, x 0, x 0. d) By the definition, x is a KKT point if there exist values of the Lagrangean multipliers such that the KKT conditions hold. Such values of the multipliers may not be unique when the gradients of the active constraints are linearly dependent. This is the case for x = (0, 0) T. It is a KKT point, because the KKT conditions hold, for instance, for λ = 0, λ = and λ 3 =. e) (3, 3) T has non-negative components and it is a KKT point. It is possible to find it, either geometrically, or from the KKT conditions. The requirement x > 0 and x > 0 implies λ = λ 3 = 0. The first of the KKT conditions give λ = and x = 3. Since λ 0, the first of the complementarity conditions implies (x ) +x = 4. From this equation, we obtain x = 3. f) For x = ( ) T, F(x,) = 8. 6

Assignment a) The problem can be reformulated with the use of the graph presented above. The arcs are labeled with upper bounds and optimal flow obtained in b). All low bounds are zero. It is required to find the maximal flow f from node to 4. b) The initial flow is zero. It requires to make two iterations of the max-flow algorithm to solve the problem. At the first iteration, modified Dijkstra s algorithm gives the path 3 4 with capacity 3. At the second iteration, modified Dijkstra s algorithm gives the path 3 4 with capacity. Thus, the maxflow f =. c) N s = {} is a min-cut. Assignment 6 a) Dijkstra s algorithm is not applicable here, because there is a negative arc cost. Notice that the nodes are topologically ordered. Then it is possible to apply Bellman s equations recursively. This gives y = 0, y = 6, y 3 =, y 4 = 3, y = 4 and the corresponding list of predecessors. Thus the shortest path is 3 4, and its length is 4. b) x = x 3 = x 34 = x 4 = and x ij = 0 for the other arcs. c) This is a minimum spanning tree problem. The minimal cost is 8. The optimal traveling strategy is the following. Two workers use one car on the path s 3. Then one of them uses one car on the path 6, and the other one uses one car on the path 4 7. The minimal number of workers is two. Assignment 7 The optimal solution is x = (,, 0), z = 7. 7