C&O 355 Lecture 16. N. Harvey

Similar documents
CS522: Advanced Algorithms

COMP331/557. Chapter 2: The Geometry of Linear Programming. (Bertsimas & Tsitsiklis, Chapter 2)

Math 414 Lecture 2 Everyone have a laptop?

Convexity: an introduction

Lecture 4: Rational IPs, Polyhedron, Decomposition Theorem

Math 5593 Linear Programming Lecture Notes

Linear programming and duality theory

Polytopes Course Notes

Convex Geometry arising in Optimization

POLYHEDRAL GEOMETRY. Convex functions and sets. Mathematical Programming Niels Lauritzen Recall that a subset C R n is convex if

However, this is not always true! For example, this fails if both A and B are closed and unbounded (find an example).

CS671 Parallel Programming in the Many-Core Era

EC 521 MATHEMATICAL METHODS FOR ECONOMICS. Lecture 2: Convex Sets

2. Optimization problems 6

A mini-introduction to convexity

arxiv: v1 [math.co] 15 Dec 2009

Conic Duality. yyye

ORIE 6300 Mathematical Programming I September 2, Lecture 3

Polyhedral Computation Today s Topic: The Double Description Algorithm. Komei Fukuda Swiss Federal Institute of Technology Zurich October 29, 2010

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

Lecture 5: Properties of convex sets

Modeling and Analysis of Hybrid Systems

Modeling and Analysis of Hybrid Systems

Chapter 4 Concepts from Geometry

MATH 890 HOMEWORK 2 DAVID MEREDITH

Lecture Notes 2: The Simplex Algorithm

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

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

Combinatorial Geometry & Topology arising in Game Theory and Optimization

Numerical Optimization

Lecture 2 - Introduction to Polytopes

maximize c, x subject to Ax b,

In this chapter we introduce some of the basic concepts that will be useful for the study of integer programming problems.

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

11 Linear Programming

Submodularity Reading Group. Matroid Polytopes, Polymatroid. M. Pawan Kumar

Mathematical and Algorithmic Foundations Linear Programming and Matchings

Lecture 6: Faces, Facets

FACES OF CONVEX SETS

ON THE EMPTY CONVEX PARTITION OF A FINITE SET IN THE PLANE**

MA4254: Discrete Optimization. Defeng Sun. Department of Mathematics National University of Singapore Office: S Telephone:

LP Geometry: outline. A general LP. minimize x c T x s.t. a T i. x b i, i 2 M 1 a T i x = b i, i 2 M 3 x j 0, j 2 N 1. where

Planar Graphs. 1 Graphs and maps. 1.1 Planarity and duality

Finite Math Linear Programming 1 May / 7

Chapter 8. Voronoi Diagrams. 8.1 Post Oce Problem

Advanced Linear Programming. Organisation. Lecturers: Leen Stougie, CWI and Vrije Universiteit in Amsterdam

Linear Programming in Small Dimensions

THREE LECTURES ON BASIC TOPOLOGY. 1. Basic notions.

Lecture 3. Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets. Tepper School of Business Carnegie Mellon University, Pittsburgh

6. Lecture notes on matroid intersection

Lecture 12: Feasible direction methods

Lecture 4: examples of topological spaces, coarser and finer topologies, bases and closed sets

Automorphism Groups of Cyclic Polytopes

Lecture 10,11: General Matching Polytope, Maximum Flow. 1 Perfect Matching and Matching Polytope on General Graphs

Lecture 1 Discrete Geometric Structures

Notes on Convex Sets, Polytopes, Polyhedra, Combinatorial Topology, Voronoi Diagrams and Delaunay Triangulations

Bounded subsets of topological vector spaces

Lesson 17. Geometry and Algebra of Corner Points

Solutions to Problem Set 2

Applied Integer Programming

2 A topological interlude

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

1 Linear programming relaxation

1 Elementary number theory

Lecture 5: Duality Theory

Homework Set #2 Math 440 Topology Topology by J. Munkres

Polar Duality and Farkas Lemma

3. The Simplex algorithmn The Simplex algorithmn 3.1 Forms of linear programs

Convex Hull Representation Conversion (cddlib, lrslib)

Notes for Lecture 24

Figure 1: Two polygonal loops

1 Linear Programming. 1.1 Optimizion problems and convex polytopes 1 LINEAR PROGRAMMING

Ma/CS 6b Class 26: Art Galleries and Politicians

IE 5531: Engineering Optimization I

11.1 Facility Location

TROPICAL CONVEXITY MIKE DEVELIN AND BERND STURMFELS

R n a T i x = b i} is a Hyperplane.

Polyhedral Computation and their Applications. Jesús A. De Loera Univ. of California, Davis

Chapter 15 Introduction to Linear Programming

9 Bounds for the Knapsack Problem (March 6)

A NOTE ON BLOCKING VISIBILITY BETWEEN POINTS

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

8 Matroid Intersection

Cardinality Lectures

Linear Programming Duality and Algorithms

arxiv: v2 [math.co] 13 Aug 2013

Chapter 2. Convex Hull. 2.1 Convexity. The following terminology should be familiar from linear algebra. Consider P R d. courses.

Lecture 12 March 4th

CPSC 536N: Randomized Algorithms Term 2. Lecture 10

Integer Programming Theory

CS675: Convex and Combinatorial Optimization Spring 2018 Convex Sets. Instructor: Shaddin Dughmi

Powered by TCPDF (

Lecture 3: Convex sets

Lecture 1. 1 Notation

CS 173, Running Time Analysis, Counting, and Dynamic Programming. Tandy Warnow

Partitions and Packings of Complete Geometric Graphs with Plane Spanning Double Stars and Paths

Section 13. Basis for a Topology

AMS : Combinatorial Optimization Homework Problems - Week V

THEORY OF LINEAR AND INTEGER PROGRAMMING

arxiv: v1 [math.co] 12 Aug 2018

Transcription:

C&O 355 Lecture 16 N. Harvey

Topics Review of Fourier-Motzkin Elimination Linear Transformations of Polyhedra Convex Combinations Convex Hulls Polytopes & Convex Hulls

Fourier-Motzkin Elimination Joseph Fourier Given a polyhedron Q µ R n, we want to find the set Q µ R n-1 satisfying (x 1,,x n-1 )2Q, 9x n s.t. (x 1,,x n-1,x n )2Q Theodore Motzkin Q is called the projection of Q onto first n-1 coordinates Fourier-Motzkin Elimination constructs Q by generating (finitely many) constraints from the constraints of Q. Corollary: Q is a polyhedron.

Elimination Example x 3 Q x 2 x 1 Project Q onto coordinates {x 1, x 2 }...

Elimination Example x 3 Q x 2 x 1 Project Q onto coordinates {x 1, x 2 }... Fourier-Motzkin: Q is a polyhedron. Of course, the ordering of coordinates is irrevelant.

Elimination Example x 3 Q x 2 x 1 Of course, the ordering of coordinates is irrevelant. Fourier-Motzkin: Q is also a polyhedron. I can also apply Elimination twice

Elimination Example x 3 Q x 2 x 1 Fourier-Motzkin: Q is also a polyhedron.

Projecting a Polyhedron Onto Some of its Coordinates Lemma: Given a polyhedron Q µ R n. Let S={s 1,,s k }µ{1,,n} be any subset of the coordinates. Let Q S = { (x s1,,x sk ) : x2q } µ R k. In other words, Q S is projection of Q onto coordinates in S. Then Q S is a polyhedron. Proof: Direct from Fourier-Motzkin Elimination. Just eliminate all coordinates not in S.

Linear Transformations of Polyhedra Lemma: Let P = { x : Ax b } µ R n be a polyhedron. Let M be any matrix of size pxn. Let Q = { Mx : x2p } µ R p. Then Q is a polyhedron. x 3 Let M = 1 0 0-1 1 0 0 0 1 P x 1 x 2

Linear Transformations of Polyhedra Lemma: Let P = { x : Ax b } µ R n be a polyhedron. Let M be any matrix of size pxn. Let Q = { Mx : x2p } µ R p. Then Q is a polyhedron. Geometrically obvious, but not easy to prove x 3 Let M = 1 0 0-1 1 0 0 0 1 Q x 1 x 2

Linear Transformations of Polyhedra Lemma: Let P = { x : Ax b } µ R n be a polyhedron. Let M be any matrix of size pxn. Let Q = { Mx : x2p } µ R p. Then Q is a polyhedron. Geometrically obvious, but not easy to prove unless you know Fourier-Motzkin Elimination! Proof: Let P = { (x,y) : Mx=y, Ax b } µ R n+p, where x2r n, y2r p. P is obviously a polyhedron. Note that Q is projection of P onto y-coordinates. By previous lemma, Q is a polyhedron.

Let SµR n be any set. Recall: S is convex if Convex Sets Proposition: The intersection of any collection of convex sets is itself convex. Proof: Exercise.

Convex Combinations Let SµR n be any set. Definition: A convex combination of points in S is any point where k is finite and Theorem: (Carathéodory s Theorem) It suffices to take k n+1. Proof: This was Assignment 2, Problem 4.

Convex Combinations Theorem: Let SµR n be a convex set. Let comb(s) = { p : p is a convex comb. of points in S }. Then comb(s) = S. Proof: S µ comb(s) is trivial.

Convex Combinations Theorem: Let SµR n be a convex set. Let comb(s) = { p : p is a convex comb. of points in S }. Then comb(s) = S. Proof: We ll show comb(s) µ S. Consider.. Need to show p2s. By induction on k. Trivial if k=1 or if k=1. So assume k>1 and k<1. Note: Call this point p Note so p is a convex combination of at most k-1 points in S. By induction, p 2S. But, so p2s since S is convex.

Let SµR n be any set. Convex Hulls Definition: The convex hull of S, denoted conv(s), is the intersection of all convex sets containing S.

Let SµR n be any set. Convex Hulls Definition: The convex hull of S, denoted conv(s), is the intersection of all convex sets containing S.

Let SµR n be any set. Definition: The convex hull of S, denoted conv(s), is the intersection of all convex sets containing S. Claim: conv(s) is itself convex. Proof: Convex Hulls Follows from our earlier proposition Intersection of convex sets is convex.

Convex Hulls Let SµR n be any set. Definition: The convex hull of S, denoted conv(s), is the intersection of all convex sets containing S. Let comb(s) = { p : p is a convex comb. of points in S }. Theorem: conv(s)=comb(s). Proof: We ll show conv(s)µcomb(s). Claim: comb(s) is convex. Proof: Consider and, where So p,q2comb(s). For 2[0,1], consider Thus p + (1- )q 2 comb(s).

Convex Hulls Let SµR n be any set. Definition: The convex hull of S, denoted conv(s), is the intersection of all convex sets containing S. Let comb(s) = { p : p is a convex comb. of points in S }. Theorem: conv(s)=comb(s). Proof: We ll show conv(s)µcomb(s). Claim: comb(s) is convex. Clearly comb(s) contains S. But conv(s) is the intersection of all convex sets containing S, so conv(s)µcomb(s).

Convex Hulls Let SµR n be any set. Definition: The convex hull of S, denoted conv(s), is the intersection of all convex sets containing S. Let comb(s) = { p : p is a convex comb. of points in S }. Theorem: conv(s)=comb(s). Proof: Exercise: Show comb(s)µconv(s).

Convex Hulls of Finite Sets Theorem: Let S={s 1,,s k }½R n be a finite set. Then conv(s) is a polyhedron. Geometrically obvious, but not easy to prove

Convex Hulls of Finite Sets Theorem: Let S={s 1,,s k }½R n be a finite set. Then conv(s) is a polyhedron. Geometrically obvious, but not easy to prove unless you know Fourier-Motzkin Elimination! Proof: Let M be the nxk matrix where M i = s i. By our previous theorem, But this is a projection of the polyhedron By our lemma on projections of polyhedra, conv(s) is also a polyhedron.

Polytopes & Convex Hulls Let P½R n be a polytope. (i.e., a bounded polyhedron) Is P the convex hull of anything? Since P is convex, P = conv(p). Too obvious Maybe P = conv( extreme points of P )?

Polytopes & Convex Hulls Theorem: Let P½R n be a non-empty polytope. Then P = conv( extreme points of P ). Proof: First we prove conv( extreme points of P ) µ P. We have: conv( extreme points of P ) µ conv( P ) µ P. Obvious Our earlier theorem proved comb(p) µ P

Polytopes & Convex Hulls Theorem: Let P½R n be a non-empty polytope. Then P = conv( extreme points of P ). Proof: Now we prove P µ conv( extreme points of P ). Let the extreme points be {v 1,,v k }. Suppose 9b 2 P n conv( extreme points of P ). Then the following system has no solution: (Finitely many!) By Farkas lemma, 9u2R n and 2R s.t. and So -u T b > -u T v i for every extreme point v i of P.

Polytopes & Convex Hulls Theorem: Let P½R n be a non-empty polytope. Then P = conv( extreme points of P ). Proof: Now we prove P µ conv( extreme points of P ). Let the extreme points be {v 1,,v k }. Suppose 9b 2 P n conv( extreme points of P ). Then the following system has no solution: (Finitely many!) By Farkas lemma, 9u2R n and 2R s.t. and So -u T b > -u T v i for every extreme point v i of P. Consider the LP max { -u T x : x2p }. It is not unbounded, since P is bounded. Its optimal value is not attained at an extreme point. This is a contradiction.

Generalizations Theorem: Every polytope in R n is the convex hull of its extreme points. Theorem: [Minkowski 1911] Every compact, convex set in R n is the convex hull of its extreme points. Theorem: [Krein & Milman 1940] Every compact convex subset of a locally convex Hausdorff linear space is the closed convex hull of its extreme points.

UNDERGRAD RESEARCH ASSISTANTSHIPS (URAS) FOR SPRING 2010 WITH THE C&O DEPT. DURATION: 3-4 months during May - August, 2010 SALARY: ELIGIBILITY: APPLICATION: BY EMAIL: Around $2,500 per month Undergrads in Math or CS, with an A standing To apply, address the following documents to Professor I.P. Goulden: - cover letter - resume - recent grade report - names and email addresses of two references send documents to t3schmid@uwaterloo.ca (No need to specify a particular area of interest, specialized background is not needed) LAST DATE FOR APPLICATIONS: DECEMBER 14, 2009 http://www.math.uwaterloo.ca/cando_dept/summerresearch/summer.shtml All of this information can be accessed from our homepage