arxiv: v1 [math.co] 15 Dec 2009

Similar documents
Convex Geometry arising in Optimization

Lecture 2 - Introduction to Polytopes

Combinatorial Geometry & Topology arising in Game Theory and Optimization

AMS : Combinatorial Optimization Homework Problems - Week V

FACES OF CONVEX SETS

Polytopes Course Notes

MATH 890 HOMEWORK 2 DAVID MEREDITH

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

Investigating Mixed-Integer Hulls using a MIP-Solver

On Unbounded Tolerable Solution Sets

maximize c, x subject to Ax b,

On the Hardness of Computing Intersection, Union and Minkowski Sum of Polytopes

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

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

C&O 355 Lecture 16. N. Harvey

CS522: Advanced Algorithms

Math 5593 Linear Programming Lecture Notes

60 2 Convex sets. {x a T x b} {x ã T x b}

Mathematical and Algorithmic Foundations Linear Programming and Matchings

Modeling and Analysis of Hybrid Systems

Modeling and Analysis of Hybrid Systems

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

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

Integer Programming Theory

Lecture 5: Duality Theory

Automorphism Groups of Cyclic Polytopes

Lecture 15: The subspace topology, Closed sets

arxiv: v1 [math.co] 12 Dec 2017

Interval-Vector Polytopes

Conic Duality. yyye

Lecture 4: Rational IPs, Polyhedron, Decomposition Theorem

4 Integer Linear Programming (ILP)

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

Polar Duality and Farkas Lemma

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

Chapter 4 Concepts from Geometry

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

Division of the Humanities and Social Sciences. Convex Analysis and Economic Theory Winter Separation theorems

arxiv: v1 [math.co] 24 Aug 2009

What is a cone? Anastasia Chavez. Field of Dreams Conference President s Postdoctoral Fellow NSF Postdoctoral Fellow UC Davis

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

On a Cardinality-Constrained Transportation Problem With Market Choice

Math 414 Lecture 2 Everyone have a laptop?

arxiv: v1 [math.co] 25 Sep 2015

Lecture 3: Convex sets

ORIE 6300 Mathematical Programming I September 2, Lecture 3

Lecture 2 Convex Sets

From acute sets to centrally symmetric 2-neighborly polytopes

CS 435, 2018 Lecture 2, Date: 1 March 2018 Instructor: Nisheeth Vishnoi. Convex Programming and Efficiency

arxiv: v1 [math.co] 17 Jan 2014

Stable sets, corner polyhedra and the Chvátal closure

Numerical Optimization

The Fibonacci hypercube

A mini-introduction to convexity

CSE 20 DISCRETE MATH WINTER

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

Monotone Paths in Geometric Triangulations

EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 6

1. Lecture notes on bipartite matching February 4th,

Convex Sets. CSCI5254: Convex Optimization & Its Applications. subspaces, affine sets, and convex sets. operations that preserve convexity

Convexity: an introduction

On the null space of a Colin de Verdière matrix

Rubber bands. Chapter Rubber band representation

Locally convex topological vector spaces

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

REGULAR GRAPHS OF GIVEN GIRTH. Contents

arxiv: v1 [math.co] 27 Feb 2015

Convex Hull Representation Conversion (cddlib, lrslib)

Optimality certificates for convex minimization and Helly numbers

THREE LECTURES ON BASIC TOPOLOGY. 1. Basic notions.

Linear Programming in Small Dimensions

Integer Programming Chapter 9

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

2. Convex sets. x 1. x 2. affine set: contains the line through any two distinct points in the set

2. Optimization problems 6

Lecture 6: Faces, Facets

Linear programming and duality theory

Three applications of Euler s formula. Chapter 10

Convex Optimization. Convex Sets. ENSAE: Optimisation 1/24

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

LECTURE 10 LECTURE OUTLINE

On the positive semidenite polytope rank

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

Lecture Notes 2: The Simplex Algorithm

arxiv: v1 [cs.cc] 30 Jun 2017

CS 372: Computational Geometry Lecture 10 Linear Programming in Fixed Dimension

Lecture 2. Topology of Sets in R n. August 27, 2008

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

Scan Scheduling Specification and Analysis

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

Finite Math Linear Programming 1 May / 7

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

Finite Math - J-term Homework. Section Inverse of a Square Matrix

Operations Research. Report On Abrams theorem. February Copyright c 2007 Department of Operations Research,

2. Convex sets. affine and convex sets. some important examples. operations that preserve convexity. generalized inequalities

RATIONAL CURVES ON SMOOTH CUBIC HYPERSURFACES. Contents 1. Introduction 1 2. The proof of Theorem References 9

Lecture notes on the simplex method September We will present an algorithm to solve linear programs of the form. maximize.

arxiv: v1 [math.co] 12 Aug 2018

Discrete Optimization. Lecture Notes 2

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

Transcription:

ANOTHER PROOF OF THE FACT THAT POLYHEDRAL CONES ARE FINITELY GENERATED arxiv:092.2927v [math.co] 5 Dec 2009 VOLKER KAIBEL Abstract. In this note, we work out a simple inductive proof showing that every polyhedral cone K is the conic hull of a finite set X of vectors. The base cases of the induction are linear subspaces and linear halfspaces of linear subspaces. The proof also shows that the components of the vectors in X can be chosen (up to their sign) to be quotients of subdeterminants of the coefficient matrix of any inequality system defining K. A matrix A R m n and a vector b R m define the polyhedron P (A,b) = {x R n : Ax b}. A polyhedron P (A, O) is a polyhedral cone. A polytope is the convex hull conv(x) = { x X λ x x : λ x R,λ x 0 for all x X, x X λ x = } of a finite set X R n, and a finitely generated cone is the (convex) conic hull ccone(x) = { x X λ x x : λ x R,λ x 0 for all x X} of a finite set X R n. A classical theorem (which is at the core of the theory of polyhedra) due to Weyl [4] and Minkowski [2] states that a subset of R n is a polyhedron if and only if it is the Minkowski sum (S + T = {s + t : s S,t T } for S,T R n ) of a polytope and a finitely generated cone. A representation P = P (A,b) (with A R m n, b R m ) of a polyhedron P R n is called an outer description, while P = conv(v ) + ccone(w) with finite sets V,W R n is an inner description. Later refinements (which are very important for the theory of linear and integer programming) state that, given one representation of a polyhedron P, there is a representation of P of the other type all of whose components are quotients of determinants of matrices (of size n n) formed from components of the given representation. In particular, if the given representation of P is rational, then one can choose a rational representation of P of the other type all of whose components have an encoding length (say, in the binary system) that is bounded polynomially in n and the maximal encoding length of any component of the given representation (see, e.g., [3, Thm. 0.2]). This is not only necessary for the polynomial solvability of the linear programming problem, but it is also crucial for establishing that the integer programming feasibilty problem is contained in the complexity class NP (see, e.g., [3, Cor. 7.b]). We denote by δ(m) the set of all determinants of submatrices (formed by arbitrary subsets of rows and columns of equal cardinality, including the empty submatrix, whose determinant is considered to be one) of a matrix M R m n, and define (M) = { p q : p,q δ(m) ( δ(m)),q 0}. Date: December 5, 2009.

2 VOLKER KAIBEL Clearly, for rational matrices M Q m n, we have (M) Q. It is well-known that, using the concepts of homogenization and polarity, one can easily derive Weyl s and Minkowski s theorem, including the above mentioned refinements, from the following result (which is itself a special case of one direction of the Weyl-Minkowski theorem). Theorem. For every matrix A R m n, there is a finite set X (A) n with P (A, O) = ccone(x). The proof of Theorem we are going to work out imitates the following obvious inductive proof of the similar statement that every bounded polyhedron P is a polytope (see []): If P consists of one point only, then the statement is clear. Otherwise, the boundary of P is the union of finitely many lower dimensional bounded polyhedra, its faces, which are polytopes by induction. The union of finite generating sets of these polytopes yields a generating set for P, since for every x in the bounded polyhedron P, any line containing x and another point from P will intersect the boundary of P in two points of which x is a convex combination. There are two issues to deal with in order to turn this basic idea into an elementary proof of Theorem. First, we do not want to dwell on the geometric concepts of faces and of the dimension of a polyhedron. This can indeed easily be avoided by allowing equations in the system defining the polyhedral cone and simply basing the induction on the number of inequalities. Second, and more important, in case of polyhedral cones K (instead of bounded polyhedra) it is, in general, not true that through every y K there is a line intersecting the boundary of K in two points (of which y is a convex combination), as one can easily see at the examples of K being a linear subspace of R n or a linear halfspace of such a subspace. The core of the proof of Theorem presented here is to show that these are the only inconvenient cases. This is the essence of the following lemma, where ker(m) = {x R n : Mx = O} denotes the kernel of the matrix M R m n (with ker(m) = R n in case of m = 0). Lemma. Let B R p n, C R q n (with p + q, n ), A = K = {x R n : Bx O,Cx = O}. ( B C ) R (p+q) n, and (i) If we have dim(ker(b) ker(c)) dim(ker(c)), then there is a finite set X (A) n satisfying K = ccone(x). (ii) Otherwise, there is some vector z ker(c) \ {O} with z K and z K. Proof. For the proof of the first part, suppose that U = ker(b) ker(c) K ker(c) has dimension dim(u) dim(ker(c)). Let B (A) n be a basis of ker(c) for which B = B U is a basis of U; due to Cramer s rule (see, e.g., [3, Cor. 3.c]), we can choose B (A) n. We have U = ccone(b ( B)), and we may assume U K, since otherwise the claim clearly follows with X = B ( B). In particular, we have dim(u) = dim(ker(c)). Hence, there is some a ker(c) \ U satisfying U = {x ker(c) : a,x = 0} and a =. Due to dimension reasons, ker(c) is the linear subspace of R n generated by U {a}. Therefore, for each y ker(c) \ U (with a,y 0), there is some u U K satisfying (see Fig. ) and thus y = u + a,y a () a = u + a,y y (2)

ANOTHER PROOF OF THE FACT THAT POLYHEDRAL CONES ARE FINITELY GENERATED 3 a, y a y a a,y y U = ker (B) u O u Figure : Illustration of relations () and (2). with u = a,y u U K. After possibly replacing a by a we have a,y > 0 for an arbitrarily chosen y K \U, which by (2) implies a K. Thus, K = ccone(u {a}) (3) holds, where is clear, and, in order to prove the reverse inclusion, due to (), it suffices to establish a,y > 0 for all y K \ U as follows. For each y K \ U we have a,y 0, where a,y < 0 yields a,y > 0, hence a K by relation (2) (since u U implies u U K, as U is a linear subspace). But from a, a K one deduces Ba O and Ba O, thus a ker(b), contradicting a U. Due to dim(u) = dim(ker(c)) there is exactly one vector v B \ B in the basis B of ker(c) that is not contained in the basis B of U. We can assume a,v > 0 (by possibly replacing v by v). Thus, by () and (2) (with y = v) we have ccone(u {a}) = ccone(u {v}), which establishes the claim in the first part (because of (3)) with X = B {v}. For the proof of the second part of the lemma, let L R n be the linear space generated by the sum B T ½ of all rows of B (see Fig. 2). Due to dim(l), the orthogonal complement L R n of L in R n has at least dimension n, hence, dim(l ker(c)) dim(ker(c)) (4) holds. If the linear subspace L ker(c) was contained in K, then (due to K {x R n : Bx O}) the linear subspace L ker(c) would be a subset of ker(b) (again, as Bx O and Bx O imply x ker(b)), which, by (4), would contradict the assumption on the dimension of dim(ker(b) ker(c)) in the second part of the lemma. Thus, there is some z (L ker(c)) \ K. Suppose we have z K, thus B( z) O, and hence Bz O. Because of ½,Bz = B T ½,z = 0 (due to z L ) this implies Bz = O, which, however, due to z ker(c), contradicts z K. Thus, z is a vector as claimed to exist in the second part of the lemma.

4 VOLKER KAIBEL L ker(c) z K B T ½ U = ker (B) Figure 2: Illustration of the proof of Part (ii) of Lemma. Note that the vector labeled B T ½ is the orthogonal projection of B T ½ to. Using Lemma, we can now easily prove Theorem by establishing, by induction on p = 0,,..., that, for every B R p n and C R q n (with p + q, n ) there is a finite set X (A) n with K = {x R n : Bx O,Cx = O} = ccone(x) ( ) B (where A = R C (p+q) n ). For p = 0, this follows readily from the first part of Lemma. For p, we may assume that there is some z ker(c) \ {O} with z, z K (otherwise, the claim again follows by the first part of Lemma ). For all i {,...,p}, let B (i) R (p ) n be the matrix that arises from B by deleting row i. By induction hypothesis, applied to K i := {x R n : B (i) x O, B i,,x = 0,Cx = O}, there is some finite set X i (A) n with K i = ccone(x i ). It suffices to show K ccone(x) for X = p i= X i (because ccone(x) K is clear). Towards this end, let x K (see Fig. 3). The set I = {i {,...,p} : B i,,z > 0} is non-empty due to Cz = O and z K. For all i I, the number λ i = B i,,x B i,,z 0 is nonnegative (due to Bx O). We have B i,,(x + λz) 0 for all 0 λ λ i with equality for λ = λ i. Now we choose i I such that λ = λ i = min{λ i : i I} holds (with λ 0). Then, we have B(x + λ z) O and B i,,(x + λ z) = 0. Thus, (due to Cz = O) we conclude x + λ z K i = ccone(x i ) ccone(x). Similarly, (due to z K) one finds some µ 0 with x + µ ( z) ccone X. Hence x, as a convex combination of x + λ z ccone X and x µ z ccone X (with λ,µ 0), is contained in ccone(x). This concludes the proof of Theorem.

ANOTHER PROOF OF THE FACT THAT POLYHEDRAL CONES ARE FINITELY GENERATED 5 x + λ z z x K x + µ ( z) z Figure 3: Convex combination of x in the polyhedral cone K (in this example being the intersection of two linear halfspaces in ker(c), viewed from above ) by two vectors from ccone (X) in the proof of Theorem. Acknowledgments. I am grateful to Matthias Peinhardt for valuable comments on a draft of this note. References [] Jiří Matoušek. Lectures on discrete geometry, volume 22 of Graduate Texts in Mathematics. Springer- Verlag, New York, 2002. [2] Hermann Minkowski. Geometry of numbers. (Geometrie der Zahlen.). Bibliotheca Mathematica Teubneriana. 40. New York, NY: Johnson Reprint Corp. vii, 256 p., 968. [3] Alexander Schrijver. Theory of linear and integer programming. Wiley-Interscience Series in Discrete Mathematics. John Wiley & Sons Ltd., Chichester, 986. A Wiley-Interscience Publication. [4] H. Weyl. Elementare Theorie der konvexen Polyeder. Commentarii math. Helvetici, 7:290 306, 935. E-mail address: kaibel@ovgu.de Otto-von-Guericke Universität Magdeburg, Fakultät für Mathematik, Universitätsplatz 2, 3906 Magdeburg, Germany