Matroid Optimisation Problems with Nested Non-linear Monomials in the Objective Function

Similar documents
1 Matroid intersection

Fakultät für Mathematik

Facets for Node-Capacitated Multicut Polytopes from Path-Block Cycles with Two Common Nodes

On Matroid Parity and Matching Polytopes

Solving the MWT. Recall the ILP for the MWT. We can obtain a solution to the MWT problem by solving the following ILP:

Complete Description of Matching Polytopes with One Linearized Quadratic Term for Bipartite Graphs

3. Linear Programming and Polyhedral Combinatorics

Technische Universität Dresden Institute of Numerical Mathematics

Separating Simple Domino Parity Inequalities

Combinatorial Optimization with One Quadratic Term: Spanning Trees and Forests

Reachability-based matroid-restricted packing of arborescences

An inequality for polymatroid functions and its applications

Separation Techniques for Constrained Nonlinear 0 1 Programming

Matroid intersection, base packing and base covering for infinite matroids

A new family of facet defining inequalities for the maximum edge-weighted clique problem

Packing Arborescences

A Randomized Rounding Approach to the Traveling Salesman Problem

Lifting theorems and facet characterization for a class of clique partitioning inequalities

4 Packing T-joins and T-cuts

Semimatroids and their Tutte polynomials

Lecture 10 February 4, 2013

Matroids and submodular optimization

3. Linear Programming and Polyhedral Combinatorics

A MATROID EXTENSION RESULT

The Strong Largeur d Arborescence

RESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices

Deciding Emptiness of the Gomory-Chvátal Closure is NP-Complete, Even for a Rational Polyhedron Containing No Integer Point

Diskrete Mathematik und Optimierung

h-vectors OF SMALL MATROID COMPLEXES

Cutting planes from extended LP formulations

Finite connectivity in infinite matroids

AN ALGORITHM FOR CONSTRUCTING A k-tree FOR A k-connected MATROID

The cocycle lattice of binary matroids

Lecture notes on the ellipsoid algorithm

Optimization methods NOPT048

Branch-and-cut Approaches for Chance-constrained Formulations of Reliable Network Design Problems

ACO Comprehensive Exam March 20 and 21, Computability, Complexity and Algorithms

Properties and Classification of the Wheels of the OLS Polytope.

3.8 Strong valid inequalities

An Introduction of Tutte Polynomial

A Geometric Approach to Graph Isomorphism

CS675: Convex and Combinatorial Optimization Fall 2014 Combinatorial Problems as Linear Programs. Instructor: Shaddin Dughmi

ACO Comprehensive Exam March 17 and 18, Computability, Complexity and Algorithms

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

Submodular Functions, Optimization, and Applications to Machine Learning

On the number of cycles in a graph with restricted cycle lengths

The Reduction of Graph Families Closed under Contraction

Cutting Plane Methods II

Non-TDI graph-optimization with supermodular functions (extended abstract)

Combinatorial Optimisation, Problems I, Solutions Term /2015

TOWARDS A SPLITTER THEOREM FOR INTERNALLY 4-CONNECTED BINARY MATROIDS IV

A First Course in Combinatorial Optimization. JON LEE IBM T.J. Watson Research Center, Yorktown Heights, New York

Chapter 2: Linear Independence and Bases

On the intersection of infinite matroids

3.7 Strong valid inequalities for structured ILP problems

Packing of Rigid Spanning Subgraphs and Spanning Trees

On improving matchings in trees, via bounded-length augmentations 1

The Chromatic Number of Ordered Graphs With Constrained Conflict Graphs

A combinatorial algorithm minimizing submodular functions in strongly polynomial time

arxiv: v3 [cs.dm] 30 Oct 2017

7. Lecture notes on the ellipsoid algorithm

A NEW SEMIDEFINITE PROGRAMMING HIERARCHY FOR CYCLES IN BINARY MATROIDS AND CUTS IN GRAPHS

Valid Inequalities and Facets for the Steinger Problem in a Directed Graph. Young-soo Myung. OR June 1991

1 Submodular functions

THE STRUCTURE OF 3-CONNECTED MATROIDS OF PATH WIDTH THREE

On the matrix-cut rank of polyhedra

The minimum G c cut problem

Submodular Functions, Optimization, and Applications to Machine Learning

ON COST MATRICES WITH TWO AND THREE DISTINCT VALUES OF HAMILTONIAN PATHS AND CYCLES

On disjoint common bases in two matroids

Dual Consistent Systems of Linear Inequalities and Cardinality Constrained Polytopes. Satoru FUJISHIGE and Jens MASSBERG.

1 Perfect Matching and Matching Polytopes

Disjoint paths in tournaments

On the Graph Bisection Cut Polytope

Duality of LPs and Applications

CS599: Convex and Combinatorial Optimization Fall 2013 Lecture 24: Introduction to Submodular Functions. Instructor: Shaddin Dughmi

MAT-INF4110/MAT-INF9110 Mathematical optimization

Optimization of Submodular Functions Tutorial - lecture I

Discrete Optimization 2010 Lecture 2 Matroids & Shortest Paths

On the Complexity of Some Enumeration Problems for Matroids

Countable Menger s theorem with finitary matroid constraints on the ingoing edges

3.7 Cutting plane methods

(II.B) Basis and dimension

Katarzyna Mieczkowska

Claw-free Graphs. III. Sparse decomposition

The Complexity of Maximum. Matroid-Greedoid Intersection and. Weighted Greedoid Maximization

COUNTING INTEGER POINTS IN POLYTOPES ASSOCIATED WITH DIRECTED GRAPHS. Ilse Fischer

ACO Comprehensive Exam October 14 and 15, 2013

Downloaded 03/01/17 to Redistribution subject to SIAM license or copyright; see

A packing integer program arising in two-layer network design

On the Generation of Circuits and Minimal Forbidden Sets

Appendix PRELIMINARIES 1. THEOREMS OF ALTERNATIVES FOR SYSTEMS OF LINEAR CONSTRAINTS

Travelling Salesman Problem

3.10 Lagrangian relaxation

MATH 326: RINGS AND MODULES STEFAN GILLE

The Triangle Closure is a Polyhedron

A Review of Linear Programming

THE MINIMALLY NON-IDEAL BINARY CLUTTERS WITH A TRIANGLE 1. INTRODUCTION

1 Some loose ends from last time

BCOL RESEARCH REPORT 05.03

Transcription:

atroid Optimisation Problems with Nested Non-linear onomials in the Objective Function Anja Fischer Frank Fischer S. Thomas ccormick 14th arch 2016 Abstract Recently, Buchheim and Klein [4] suggested to study polynomial-time solvable optimisation problems with linear objective functions combined with exactly one additional quadratic monomial. They concentrated on special quadratic spanning tree or forest problems. We extend their results to general matroid optimisation problems with a set of nested monomials in the objective function. The monomials are linearised and we study the corresponding polytope with the aim to better understand the structure of polytopes arising from such linearisations and to provide strengthened cutting planes as well as separation algorithms for linearisations of matroid optimisation problems with general polynomial objective function. Extending results by Edmonds [7] for the matroid polytope we present a complete description for the linearised polytope. Indeed, apart from the standard linearisation one needs appropriately strengthened rank inequalities satisfying certain non-separability conditions. The separation problem of these extended rank inequalities reduces to a submodular function minimisation problem. In the case of exactly one additional non-linear monomial we completely characterise the facetial structure of the associated polytope. Keywords: polyhedral combinatorics, matroid, matroid polytope, linearisation, polynomial optimisation, separation SC2010: 90C57, 90C27 1 Introduction Some combinatorial optimisation problems are well-understood in the sense that they can be solved in polynomial time and the complete description of the associated polytope is known. Examples are the spanning tree problem, the optimisation over matroids or over the intersection of two matroids as well as the matching problem, see, e. g., [7, 8, 23]. Optimisation problems with general polynomial objective functions are often solved by introducing new variables for each monomial, see, e. g., [10, 12, 13, 16, 25] and considering the associated linear problem, for the quadratic unconstrained case see, e. g., [21]. Approaches for solving some kinds of non-linear matroid optimization problem can be found, e. g., in [1, 17, 19]. In this paper we study a simpler but related problem with the aim to better understand the structure of polytopes arising from linearisations. Instead of considering an arbitrary polynomial objective function, we restrict to linear terms and a set of nested monomials or in a special case one additional monomial of arbitrary degree. We will show that studying the polynomially-solvable matroid optimisation problem, the corresponding problem with this particular polynomial objective function is solvable in polynomial time, too (see Remark 2). But although the well-known separation equals University of Göttingen, Institute for Numerical and Applied athematics, anja.fischer@mathematik. uni-goettingen.de University of Kassel, Institute for athematics, frank.fischer@uni-kassel.de Sauder School of Business, University of British Columbia, tom.mccormick@sauder.ubc.ca 1

optimisation result [15] then implies that the separation problem for our problem can be solved in polynomial time, we want to find the exact complete description of the associated linearised polytope and combinatorial separation algorithms for the newly derived cutting planes. Our main motivation is that this study can also be helpful for (matroid) problems with general linearised polynomial objective functions because the strengthened inequalities for the case with linearised nested monomials usually remain valid for the general case and the respective separation algorithms can be used as well (but in most cases the facet defining inequalities of the case with few linearised monomials will not remain facet defining for the general case). One advantage of our approach in comparison to general linearisation techniques, like the procedure of Sherali-Adams [24, 25], is that we strengthen the constraints of the original linear problem without introducing new variables for non-linear monomials that are not explicitly contained in the objective function. Focusing on exactly one additional non-linear monomial has first been suggested by Buchheim and Klein [3, 4] for the quadratic case. Results for the spanning tree and forest problem with exactly one additional quadratic monomial can be found in [3, 4, 9]. Buchheim and Klein also showed that the use of strengthened cutting planes can reduce the number of nodes of a branch-and-cut tree in the solution of general quadratic spanning tree problems significantly. In this paper we extend the results for the quadratic forest problems to general matroid optimisation problems with one monomial of arbitrary degree and to a set of nested monomials to be described below. In the following we will describe our setting in more detail. We start with a short review of matroids and their associated polytopes [20, 26]. In 1970 Edmonds [7] presented a complete description of the polytope P associated with a matroid = (E, I) (the convex hull over all incidence vectors of independent sets) with finite set E and a family I of subsets of E satisfying the three matroid properties (1) I, (2) T I, S T S I, (3) S, T I, S < T : e T \ S with S {e} I. A pair = (E, I) satisfying only (1) and (2) is called independence system. We denote by P(X) the power set of a finite set X. Then the rank function r : P(E) N 0 of an independence system = (E, I) is defined as r(t ) = max{ B : B T, B I}. For some T E a set B with B T, B I and r(t ) = B is called basis (of T ). It is well known that the rank function of a matroid has the following properties (see, e. g., [20, 26]): (R1) subcardinality: X E : r(x) X, (R2) monotonicity: X Y E : r(x) r(y ), (R3) submodularity: X, Y E : r(x Y ) + r(x Y ) r(x) + r(y ). Furthermore, we will often make use of the following well-known theorem. Theorem 1 (Brualdi [2], also Strong Basis Exchange Theorem). Let S, S be two different bases of a matroid. Then for any e S \ S there is an f S \ S so that (S \ {e}) {f}, (S \ {f}) {e} are both bases of. For further results on matroid theory we refer the reader to the excellent books [20, 26]. Let = (E, I) be a matroid with rank function r : P(E) N 0 and Ē := {e 1,..., e k } E, Ē 2, be a set with cardinality k. We set Ēj := {e 1,..., e j } for all j {1,..., k} and Ē0 :=. Throughout the paper we assume that the following conditions are satisfied: e E : {e} I and Ē I, (1) 2

as well as that is given via a rank oracle. Furthermore, we use the following notation. Let S, T E and let e E. If the meaning is clear from the context, we will often simply write S + T and S T instead of S T and S \ T as well as S + e and S e instead of S {e} and S \ {e}. We can now state our optimisation problem. Given an index set K := {k 1,..., k l } {2,..., k} for the (non-linear) monomials and cost functions c: E R and c: K R we consider the maximisation problem { maximise e E c(e) x(e) + ( k j c(k j ) i=1 } x(e i ) ), x P {0, 1} E (OPT-K) where P is the polytope that is associated with the matroid (note that, w. l. o. g., we will always assume that k l = k). Furthermore we set k 0 := 0. An important special case of (OPT-K) is the problem { maximise e E ( k } c(e) x(e) + c(k) x(e i ) ), x P {0, 1} E i=1 (OPT-k) where K = {k} in (OPT-K) with exactly one non-linear monomial. Although we are mainly interested in the structure of the associated polytopes, we first note that these problems can be solved in polynomial time. This motivates us by the separation equals optimisation result to detect the complete descriptions of the associated linearised polytopes. Remark 2. The optimisation problem (OPT-K) can be solved in polynomial time. Proof. We solve (OPT-K) with K = {k 1,..., k l } in two main steps. First, we determine for each i {1,..., l} a cost maximal solution Xi I, Ēk i 1 Xi such that at most k i k i 1 1 elements of Ēk i \ Ēk i 1 are contained in Xi. This can be done by solving a matroid intersection problem of matroid and a matroid i = (E, I i ) with I i = {X E : X (Ēk i \ Ēk i 1 ) k i k i 1 1}. The objective value of problem (OPT-K) evaluated for Xi equals e Xi c(e) + i 1 c(k j) in this case. Furthermore we solve one additional matroid optimisation problem in order to determine a set Xl+1 I, Ēk l Xl+1 such that e Xl+1 c(e) is maximised. The objective value of problem (OPT-K) evaluated for Xl+1 is then e Xl+1 c(e) + l c(k j ). At the end we take a set Xi, i = 1,..., l + 1, with maximal total objective value. Because (OPT-k) is a special case of (OPT-K) the problem (OPT-k) can be solved by one matroid optimisation and one matroid intersection problem. Alternatively one can solve (OPT-K) by comparing the optimal values of k + 1 matroid optimisation problems. Our aim is to present a complete description of the linearised polytopes P Ē,K := conv {(x, y) {0, 1} E+K : x P, y kj = P Ē,k := conv {(x, y) {0, 1} E+{k} : x P, y = k j i=1 k i=1 } x(e i ), j = 1,..., l, } x(e i ). Note, if K = {k}, we often simply write y instead of y k or y k1. Starting point of our work is a result of Edmonds on the polytope P associated with a matroid. Before repeating this result we need some notation. For any S E we define the characteristic vector χ S {0, 1} E by χ S (e) = { 1, e S, 0, e / S. 3

(In slight abuse of notation, we often do not distinguish between the characteristic vector of an independent set J I and the set J.) Then the polytope associated with the matroid is P = conv{χ S {0, 1} E : S I}. Theorem 3 (Edmonds, [7]). Let be a matroid with rank function r : P(E) N 0, then The inequalities P = { x R E + : } x(e) r(t ), T E. e T x(t ) r(t ), T E (2) with x(t ) := e T x(e) are often called rank constraints/inequalities. Together with the nonnegativity constraints x(e) 0, e E, they form a complete description of P. Edmonds also characterised which sets T lead to facet defining inequalities of P. The paper is structured as follows. In Section 2 we present complete descriptions of P Ē,k and. A main part of the proof consists of showing that apart from the standard linearisation, i. e., all constraints that are needed in the unconstrained case (I = P(E)) as well, all coefficients of facet defining inequalities a T x + α T y b of the polytopes are non-negative. We start with P Ē,K proving this result for P Ē,k Ē,K and extend this to P with K = {2,..., k}. Based on the complete for description in this case, we show that it is possible to derive a complete description of P Ē,K arbitrary sets K {2,..., k}. In Section 3 we study the facetial structure of P Ē,K. So we have a closer look at the so called trivial facets and give some necessary conditions for some specially lifted rank inequalities to be facet defining for P Ē,K Ē,k. For P we even fully characterise all facets. In Section 4 we show that the separation problem for P Ē,K can be solved in strongly polynomial time by complete enumeration of the trivial facets and by submodular function minimisation. Finally, in Section 5, we describe further possible extensions of our results and give suggestions for future work. Some rather technical proofs are deferred to the Appendix. 2 Complete Description of P Ē,K In this section we provide a complete description of P Ē,K for arbitrary sets K {2,..., k}. Before we present a formulation of (OPT-K) we need some definitions. Definition 4. Let T E. The closure of set T is the set The set T is called closed if T = cl(t ). cl(t ) := T {e E : r(t ) = r(t + e)}. Later we will need the following simple result that follows directly from the last definition and the submodularity of the rank function: Observation 5. Let T E and let e cl(t ). Then r(t + S) = r(t + S + e) for all S E. We will present strengthened versions of the rank inequalities (2). In order to express the coefficients of the y-variables corresponding to the linearised monomials we introduce the following objects: 4

e 1 e 3 e 4 e 1 e 2 e 4 e 3 e 2 e 4 e 2 e 1 e 2 e 3 e 4 e 3 e 1 α 2 (T ) 0 α 3 (T ) 1 α 4 (T ) 1 α 2 (T ) 1 α 3 (T ) 0 α 4 (T ) 1 α 2 (T ) 0 α 3 (T ) 1 α 4 (T ) 0 α 2 (T ) 1 α 3 (T ) 1 α 4 (T ) 1 Figure 1: Illustration of α i (.). We consider the graphical matroid on a complete graph G = (V, E) with six nodes and Ē = {e 1, e 2, e 3, e 4 }. The edges of T E are drawn solid and the edges of Ē are drawn with dotted lines. All other edges are omitted in the picture for the sake of clarity. Definition 6. Let T E. We define α 0 (T ) := 0 as well as α i (T ) := {e i } \ T + r(t + Ēi 1) r(t + Ēi), i = 1,..., k. Furthermore we write α i,j (T ) := jm=i α m(t ) for i, j {1,..., k}, i j. Fig. 1 illustrates α i (T ) for some sets T E for the graphical matroid on a complete graph G = (V, E) with six nodes and Ē = {e 1, e 2, e 3, e 4 }. The edges of T E are drawn solid and the edges of Ē are drawn with dotted lines. The following results follow directly from the definition of α i (.). For details we refer to the Appendix. Observation 7. Let T E. (A1) If e i T, i {1,..., k}, then α i (T ) = 0. (A2) We have α i (T ) {0, 1} for all i {1,..., k}. (A3) Let i {1,..., k}. If j N, j > i, then α j (T ) = α j (T + e i ). (A4) We have α 1,j (T ) = Ēj \ T + r(t ) r(t + Ēj) for all j {1,..., k}. (A5) Let i {1,..., k} and e cl(t ) \ {e i }. Then α i (T ) = α i (T + e) for all i {1,..., k}. (A6) Let i {1,..., k} and e i cl(t ) \ T. Then α i (T ) = 1 = α i (T + e i ) + 1. Studying (OPT-K) we define for each set S E the j-th (linearised) monomial term of S by for j {1,..., l}. So P Ē,K y kj (S) = can also be written as { 1, Ēkj S, 0, Ē kj S, } P Ē,K {(χ = conv S, y k1 (S),..., y kl (S)) {0, 1} E+K : S I. With these definitions we are able to give a formulation of (OPT-K): 5

Lemma 8. A formulation of (OPT-K) is given by: maximise k j i=k j 1 +1 c(e)x(e) + c(k j )y kj e E k 1 y kl 0, (3) y kj y kj 1 0, j {2,..., l}, (4) x(e j ) y k1 k 1 1, (5) x(e i ) + y kj 1 y kj k j k j 1, j {2,..., l}, (6) x(e i ) + y kj 0, j {1,..., l}, i {k j 1 + 1,..., k j }, (7) x(e) 0, e E \ Ē, (8) x(t ) + α kj 1 +1,k j (T )y kj r(t ), T E, (9) Note, we often denote (y k1,..., y kl ) by y. (x, y) {0, 1} E+K. (10) Proof. Constraints (3) (7) are the standard linearisation of the non-linear monomials in this case and (8) are the non-negativity conditions for the x-variables. By the integrality conditions (10) and the complete description of the matroid polytope P in Theorem 3 it remains to show the validity of the strengthened rank inequalities (9). Let T E be an arbitrary set and let S I. If Ēk 1 S, then the validity follows from the validity of (2) for P. So let Ēk i S, i {1,..., l}, and Ēk i+1 S in the case i < l. Then by (A4) we have χ S (e) + α kj 1,k j (T )y kj (S) = S T + Ēk i \ T + r(t ) r(t Ēk i ) e T =α 1,ki (T ) (by (S T ) (Ēk i \ T ) = ) = (S T ) (Ēk i \ T ) r(t Ēk i ) + r(t ) (by S I, Ēk i S) = r((s T ) (Ēk i \ T )) r(t Ēk i ) + r(t ) (by ((S T ) (Ēk i \ T )) T Ēk i ) r(t Ēk i ) r(t Ēk i ) + r(t ) = r(t ). A visualisation of the extended rank inequalities (9) on the example of a graphical matroid is given in Fig. 2. Remark 9. With the last result the formulation of (OPT-k) reduces to max e E c(e)x(e) + c(k)y x(ē) y Ē 1, (11) x(e) + y 0, e Ē, (12) x 0, y 0, (13) x(t ) + α 1,k (T )y r(t ), T E, (14) x, y binary. 6

e 1 e 1 e 1 f 1 e 4 e 3 e 2 e 2 f 1 f 2 e 2 f 3 f 2 e 3 f 2 f 1 f 3 e 3 T = {f 1, f 2, f 3}, Ē = {e 1, e 2, e 3, e 4} case K = {4}: x(f 1) + x(f 2) + x(f 3) + 2y 4 2 case K = {3, 4}: x(f 1) + x(f 2) + x(f 3) + y 3 + y 4 2 T = {f 1, f 2}, Ē = {e 1, e 2, e 3} case K = {3}: x(f 1) + x(f 2) + 2y 3 2 case K = {2, 3}: x(f 1) + x(f 2) + y 2 + y 3 2 T = {f 1, f 2, f 3}, Ē = {e 1, e 2, e 3} case K = {2}: x(f 1) + x(f 2) + x(f 3) 3 case K = {2, 3}: x(f 1) + x(f 2) + x(f 3) + y 3 3 Figure 2: Visualisation of the extended rank inequalities (9). We consider the graphical matroid on an undirected graph with six nodes with different sets K for a given Ē. The set T is visualised with bold lines, the edges associated with Ē with dotted lines. All other edges are omitted in the pictures for the sake of clarity. Note, all sets T are closed, see also Remark 10. Remark 10. In the formulation and the complete description of P Ē,K Ē,k and P it is sufficient to restrict to closed sets T E for the extended rank inequalities (9) and (14), respectively, i. e., x(t ) + α kj 1 +1,k j (T )y kj r(t ), T E, cl(t ) = T, (15) for arbitrary sets K and in the case K = {k} we can restrict to x(t ) + α 1,k (T )y r(t ), T E, cl(t ) = T. (16) For proving this results, let T E be an arbitrary set. Then we define T := cl(t ) \ T. By (A5) and (A6) we get α i (T ) = α i (cl(t )) + {e i } \ T, i = 1,..., k, because α i (T ) = α i (T + e) for some e E as long as e e i and α i (T ) = 1 if e i T. Let ( x, ȳ) be a point that satisfies (3) (8) as well as (15), then x(t ) + α kj 1 +1,k j (T )ȳ kj x(t ) + x( T Ē) {e kj 1 +1,..., e kj } \ T ȳ kj + α kj 1 +1,k j (T )ȳ kj } {{ } 0 by (7) x(t ) + x( T Ē) + x( T \ }{{ Ē) + α kj 1 +1,k } j (cl(t )) ȳ kj 0 by (8) =α kj 1 +1,k j (T ) x(cl(t )) + {e kj 1 +1,...,e kj }\ T α kj 1 +1,k j (cl(t ))ȳ kj (15) r(cl(t )) = r(t ). So all extended rank inequalities (9) for non-closed sets T E are implied by (3) (8) and (15). In the remaining section we will prove that (3) (8) and (15) are a complete description of Ē,k and therefore that (11) (13) and (16) are a complete description of P. The proof for P Ē,K 7

consists of three steps. First, we consider the special case K = {2,..., k}. Indeed one can show that then each facet defining inequality of P Ē,K that is not positive multiple of one of (3) (8) is an inequality of the type a T x + α T y b with a 0, α 0. This result then allows to P Ē,K derive the complete description of P Ē,K description of P Ē,K in this special case. Finally, we consider the complete for arbitrary sets K {2,..., k}. We start with an observation about the structure of facet defining inequalities a T x + α T y b that are not a positive multiple of one of (3) (8). of P Ē,K Observation 11. Let a T x + α T y b be a facet defining inequality of P Ē,K that is not a positive multiple of the standard-linearisation and the non-negativity constraints (3) (8). Then we know a 0. Proof. Let a T x + α T y b be a facet defining inequality of P Ē,K of one of (3) (8). We consider two cases. that is not a positive multiple 1. a(e) < 0 for some e E \ Ē: Because at x + α T y b is not a positive multiple of x(e) 0 (8), there exists a set S I with e S and a T χ S + α T y(s) = b. Assumption e / Ē implies y ki (S) = y ki (S e), i = 1,..., l, thus b = a T χ S + α T y(s) < a T χ S e + α T y(s e). This contradicts the validity of a T x + α T y b. 2. a(e) < 0 for some e = e i Ē: Because at x + α T y b is not a positive multiple of y km x(e i ) for m {1,..., l} such that i {k m 1 + 1,..., k m }, there exists some S I with e S and y km (S) = 0. With y kj (S) = y kj (S e), j {1,..., l}, this implies b = a T χ S + α T y(s) < a T χ S e + α T y(s e), again a contraction to validity. Next, we consider the coefficient of the y-variables of facet defining inequalities a T x + α T y b. Although the result for P Ē,k Ē,K will be implied by the results for P we will present next the proof in detail because it is much shorter than in the general case and already shows the main ideas. Lemma 12. Each facet defining inequality a T x + αy b of P Ē,k of one of (11) (13) satisfies α 0. that is not a positive multiple Proof. Let a T x + αy b be a facet defining inequality of P Ē,k that is not a positive multiple of one of (11) (13). By Observation 11 we may assume a 0. Because a T x + αy b is not a positive multiple of x(ē) y Ē 1 there exists an S I with Ē \ S 2 and at χ S = b. Similarly, because a T x + αy b is not a positive multiple of y 0, there exists an S I with Ē S and a T χ S + α = b. Set T := {e E \ Ē : a(e) > 0} (S Ē). We may assume S T (otherwise use S T because a 0 implies a(e) = 0 for all e / T ) and r(s) = r(t ) (otherwise there would be an e T \ S with S + e I, implying a T χ S+e = a T χ S + a(e) > a T χ S = b). Analogously we may assume S T Ē as well as r(s ) = r(t Ē). We distinguish two cases. 1. r(t ) < r(t Ē). Then there exists an e Ē \ T with r(t + e) > r(t ) = r(s), so S + e I. Therefore 0 = b b a T χ S+e a T χ S = a(e), so a(e) = 0. But since y(s ) = 1 and e Ē S, this implies y(s e) = 0 and so a T χ S = a T χ S e b = a T χ S + α. This proves α 0. 2. r(t ) = r(t Ē). By assumption S, S are both bases of T Ē. We may assume that S S is maximal. By Theorem 1 there exists for each f S \ S an e S \ S so that S e + f, S + e f I. Because a T χ S e+f b = a T χ S we get a(e) a(f). We choose f (S \ S) Ē, then at χ S +e f b = a T χ S + α implying a(e) a(f) + α. This proves α 0. 8

Now, we extend the previous result to P Ē,K in the case K = {2,..., k}. This is done in three main steps, which make use of the following observation. For further details and the complete, rather technical proofs we refer to the Appendix. Observation 13. Let a T x + α T y b be a facet defining inequality of P Ē,K and let S be a root of it. Let X E. We define with K = {2,..., k} lhs(x) := a T χ X + α T y(x) and ã(x, e) := lhs(x) lhs(x e), e X. Then the following conditions are satisfied: (i) for all e S : ã(s, e) 0, (ii) for all e / S with S + e I: ã(s + e, e) 0, (iii) for all e S, f / S with S e + f I: ã(s, e) ã(s e + f, f). Proof. Immediately clear because otherwise a T x + α T y b would not be valid for the sets S e, S + e or S e + f, respectively. First one can show that for each facet defining inequality of P Ē,K exists a root with a certain structure. (K = {2,..., k}) there Lemma 14. Let a T x + α T y b be a facet defining inequality of P Ē,K with K = {2,..., k} that is not a positive multiple of one of (3) (8) and let m {2,..., k}. Then there exists a root S of a T x + α T y b with the structure Ē m \ {e j, e m } = S Ēm for some j {1,..., m 1}. This implies the non-negativity of particular sums of coefficients of facet defining inequalities satisfying certain structures. Lemma 15. Let a T x + α T y b be a facet defining inequality of P Ē,K with K = {2,..., k} that is not a positive multiple of one of (3) (8). Let m {2,..., k} so that there is a root S of a T x + α T y b satisfying S Ēm = Ēm \ {e j, e m } for some j {1,..., m 1}. Then α m +... + α h 0 for all h {m,..., k} with h = k e h+1 / S. This leads to the following result. Lemma 16. Let a T x + α T y b be a facet defining inequality of P Ē,K with K = {2,..., k} that is not a positive multiple of one of (3) (8). Then α m +... + α h 0 for all m, h {2,..., k}, m h. The following corollary summarises the previous lemmas. Corollary 17. Let a T x + α T y b be a facet defining inequality of P Ē,K that is not a positive multiple of one of (3) (8). Then a 0 and α 0. with K = {2,..., k}, Proof. Let a T x + α T y b be a facet defining inequality of P Ē,K with K = {2,..., k} satisfying the requirements. Then a 0 follows from Observation 11 and α 0 follows from Lemma 16 for m = h. This allows us to prove our main result for P Ē,K in the case K = {2,..., k}. Theorem 18. The inequalities (3) (8) and (15) are a complete description of P Ē,K, K = {2,..., k}. 9

Proof. Let a T x + α T y b be a facet defining inequality of P Ē,K, K = {2,..., k}, that is not a positive multiple of one of (3) (8) and (15). Then Corollary 17 implies a 0 and α 0. Let T := {e E : a(e) > 0} and T := cl(t ). If we can prove for some S I with a T χ S + α T y(s) = b that there exists an e T \ S so that S + e I, then a T χ S+e + α T y(s + e) > a T χ S + α T y(s) = b, hence the inequality would not be valid, a contradiction. We distinguish two cases. We consider the constraint (15) associated with T. Because a T x + α T y b is not a positive multiple of this constraint, there exists an S I with e T χ S(e) + α(t ) T y(s) < r(t ) and a T χ S + α T y(s) = b. In particular, r(s T ) = e T χ S(e) implies r(s T ) + α(t ) T y(s) < r(t ). (17) Let î {1,..., k} be maximal so that y i = 1 i î (î = 1 means that Ē2 S). Note that î > 1 implies Ēî S and î < n implies Ēî+1 S as well. We first consider the case î > 1. Then we may assume S T Ēî (otherwise use S := S (T Ēî) with a T χ S + α T y( S) = a T χ S + α T y(s) = b). Note that α T y(s) = α 1,î (T ), so (17) and (A4) imply r(s T ) < r(t Ēî) Ēî \ T. Therefore r(s) = r(s T ) + S \ T < r(t Ēî) Ēî \ T + S \ T = r(t Ēî) = r(t Ēî). =0 by Ēî S T Ēî Here the last equation follows from Observation 5. So there exists an e (T Ēî) \ S = T \ S with S + e I, which contradicts the validity of a T x + α T y b. It remains to consider the case î = 1. But then y i (S) = 0, i = 2,..., k, and we know with S := S T that r(s ) r(s T ) < r(t ) = r(t ) by (17). Hence there exists an e T \ S with S + e I, a contradiction. Remark 19. If we consider (OPT-k), i. e., (OPT-K) with K = {k}, the proof of Theorem 18 in combination with Lemma 12 also shows that (11) (13) and (16) define a complete description of P Ē,k because setting î = k implies α 1,k(T ) = î α j (T ). In our last step we extend the previous results to the complete description of P Ē,K with arbitrary K = {k 1,..., k l } {2,..., k}. Theorem 20. Inequalities (3) (8) and (15) are a complete description for P Ē,K set K = {k 1,..., k l } {2,..., k}. with arbitrary Proof. We have already proved that (3) (8) and (15) are a complete description for P Ē,K if K = {2,..., k}. Let j {1,..., l} so that either j 2 and k j k j 1 2 or j = 1 and k 1 > 2. Set { kj 1 + 1, j 2, ˆk := 2, j = 1. It will be convenient throughout the proof to identify formally ȳ 1 := x(e 1 ). We show that the complete description for K can be derived from the complete description for P Ē,K with K = K {ˆk}. Let ( x, ȳ) be a point that satisfies (3) (8) as well as (15) for (OPT-K). We consider two cases: Case 1 If setting ȳˆk := ȳ kj, the point ( x, ȳ k1,..., ȳ kj 1, ȳˆk, ȳ kj,..., ȳ kl ) is feasible for P Ē,K, the statement immediately follows from the complete description of the known case. The only inequality that might be violated by this definition of ȳˆk is ȳˆk 1 + x(eˆk) 1 + ȳˆk. The validity of all other constraints follows directly from the feasibility of ( x, ȳ) and the fact that x(e) 1, e E (by (15)). For (15) with T E, T closed, we just split (T )ȳ αˆk,kj kj into αˆk,ˆk(t )ȳˆk + (T )ȳ αˆk+1,kj kj (note, in the case ˆk = 2 the coefficient of ȳˆk is α 1,ˆk(T ), but T closed implies α 1 (T ) = 0 and so α 1,ˆk(T ) = αˆk,ˆk(t )). 10

Case 2 Otherwise, we set ȳˆk := ȳˆk 1 + x(eˆk) 1 and so ȳˆk > ȳ kj. But then ( x, ŷ) := ( x, ȳ k1,..., ȳˆk 1, ȳˆk, ȳ kj,..., ȳ kl ) might possibly violate (15) of the formulation of P Ē,K. All other constraints of the formulation of P Ē,K ((3) (8)) remain satisfied by this definition of ȳˆk, either trivially or because k j i=ˆk+1 x(e i ) + x(eˆk) + ȳˆk 1 ȳ kj + k j k j 1 =ȳˆk+1 k j i=ˆk+1 x(e i ) + ȳˆk ȳ kj + (k j ˆk). Let T E, cl(t ) = T, (see also Remark 10) be a closed set and consider the associated rank constraint (15). If αˆk(t ) = 0, then the validity of this rank constraint for P Ē,K follows directly from the validity of the associated rank constraint for P Ē,K. So we may assume αˆk(t ) = 1. This implies r(t + Ēˆk 1 ) = r(t + Ēˆk) and with T := T + eˆk α 1,m (T ) α 1,m (T ) = 1 r(t + Ēm) + r(t + Ēm) 0, m = 1,..., ˆk 1, α 1,ˆk 1 (T ) α 1,ˆk 1 (T ) = 1 r(t + Ēˆk) + r(t + Ēˆk 1 ) = 1. (18) Let γ i := ȳ ki ȳ ki+1 0, i = 1,..., j 2. Then we can write ȳ ki = ȳˆk 1 + j 2 m=i γ m. Using these relations we get for the rank constraint associated with T for K j 1 x(t ) + i=1 α ki 1 +1,k i (T )ȳ ki + αˆk(t ) =1 = x(t ) + α 1,ˆk 1 (T ) ȳˆk 1 + =α 1,ˆk 1 (T ) 1 ȳˆk + αˆk+1,kj (T ) =αˆk+1,kj (T ) j 2 α 1,ki (T ) i=1 α 1,ki (T ) ȳ kj + α ki 1 +1,k i (T ) ȳ ki i=j+1 =α ki 1 +1,k i (T ) γ i + ȳˆk + (T ) αˆk+1,kj =αˆk,kj (T ) and using ȳˆk 1 = yˆk + x(eˆk) 1 and the definition of the γ + i=j+1 x(t ) + α ki 1 +1,k i (T )ȳ ki 1 i=1 x(cl(t )) + α ki 1,k i (cl(t ))ȳ ki 1 r(cl(t )) 1 = r(t ). i=1 ȳ kj α ki 1 +1,k i (T )ȳ ki For the last step we used that each rank constraint associated with T is always dominated by the rank constraint associated with cl(t ) by Remark 10. So the result follows again by the known complete description of P Ē,K. So the complete description of P Ē,K from the complete description of P Ē,K for arbitrary sets K {2,..., k} can easily be derived with K = {2,..., k} be adapting the constraints of the standard linearisation and by adjusting the extended rank inequalities. Indeed, the coefficient of y kj, k j K, equals the sum of the coefficients of variables y i, i = k j 1 + 1,..., k j, in the problem with K. Although feasibility of the extended rank inequalities is trivially preserved by this approach, it is remarkable that we even get the complete description of P Ē,K in this way. 11

3 Facetial Structure of P Ē,K and P Ē,k The polytope P of some matroid is well-understood. In 1970 Edmonds [7] presented a complete description of P and characterised its facetial structure, see also [14]. We will extend these results to the polytopes P Ē,k Ē,K and P and study their facetial structures in more detail. This includes the dimension of the polytopes as well as results on the facetness of the constraints of the standard linearisation. Furthermore, we have a closer look at the extended rank inequalities (9). If a closed set T does not satisfy certain non-separability conditions (to be defined below) the associated rank inequality cannot be facet defining for P Ē,K Ē,k. For P we even exactly characterise the facet defining inequalities and provide a minimal description. We start with studying the dimension and the trivial facets of P Ē,K. One can easily prove the following result: is full-dimensional. The in- Observation 21. Under the assumptions (1) the polytope P Ē,K equalities (3) (8) define facets of P Ē,K. In the next step we consider the extended rank inequalities (14). These are extensions of the rank inequalities (2). Before repeating results on the facetness of (2) for P we need a definition, which will be extended to our setting later on. Definition 22. Let T E. A set T is called separable if there exists a partition T 1, T 2 T, T 1 T 2 = T, T 1 T 2 = such that r(t 1 ) + r(t 2 ) = r(t ), otherwise it is called non-separable. Indeed, all closed and non-separable sets T E lead to facet defining rank inequalities (2) of P [7, 14]. They also lead to facet defining inequalities (9) of P Ē,K whenever α k j 1 +1,k j (T ) > 0 for all j {1,..., l}. Proposition 23. Let T E be a closed and non-separable set associated with a facet defining inequality (2) of P. If α kj 1 +1,k j (T ) > 0 for all j {1,..., l}, then the inequality (9) is facet defining for P Ē,K, too. Proof. Because x(t ) r(t ) is facet defining for P there exist pairwise different S i I, S i T, r(s i ) = r(t ), i = 1,..., d, d := dim P. By α kj 1 +1,k j (T ) > 0, j {1,..., l}, we know (Ēk j \Ēk j 1 ) T, j {1,..., l}, hence (Ēk j \Ēk j 1 ) S i for all i {1,..., d}, j {1,..., l}. If we choose for each j {1,..., l} a set Ŝj I with Ēk j Ŝj T Ēk j and Ŝj = r(ŝj) = r(t Ēk j ), the left-hand side of x(ŝj T ) + l α kj 1 +1,k j (T )y kj (Ŝj) r(t ) computes to r(ŝj T ) + α 1,kj (T ) = r(t Ēk j ) Ŝj \ T + α 1,kj (T ) = r(t Ēk j ) Ēk j \ T + α 1,kj (T ) = r(t ) by (A4). So there exist d + l affinely independent sets that are roots of (9), which proves the statement. Considering the first type of extended subtour elimination constraints (constraints (5) in [4]) for the spanning forest problem with one linearised quadratic monomial in [4], the facetness of these inequalities follows directly from the previous result. But as already observed by Buchheim and Klein [4] for this problem (see constraints (6) in [4]), not all facet defining inequalities of are directly related to facet defining inequalities of P. So we need further properties that allow us to describe the facetial structure of P Ē,K. P Ē,K Definition 24. Let T E be closed and K = {k 1,..., k l } {2,..., k}. A set T is called (Ē, K)-separable if there exists a partition T 1, T 2 T, T 1 T 2 =, T 1 T 2 = T with r(t ) = r(t 1 ) + r(t 2 ), and α kj 1 +1,k j (T ) = α kj 1 +1,k j (T 1 ) + α kj 1 +1,k j (T 2 ), j = 1,..., l, (19) otherwise it is called (Ē, K)-non-separable. 12

e 1 e 2 f 1 f 2 ({e 1, e 2 }, {2})-separable e 1 f 1 f 2 e 2 ({e 1, e 2 }, {2})-non-separable Figure 3: Visualisation of (Ē, K)-(non)-separability. We consider the graphical matroid on an undirected graph with four nodes. Although set T (bold lines) is the same on both sides, depending on Ē (dotted lines), it is (Ē, {2})-separable or not. Note, T is separable with T 1 = {e 1 }, T 2 = {e 2 } by Definition 22. Note, the definition of separability by Edmonds [7] only included the first equation of the last definition (see Definition 22). A visualisation of the new setting is given in Fig. 3. A closed set T E that is (Ē, K)-separable does not lead to a facet defining rank inequality (9) of P Ē,K. Observation 25. Let (9) be a facet defining inequality of P Ē,K T E. Then T is (Ē, K)-non-separable. associated with a closed set Proof. Let T E be closed and assume that T is (Ē, K)-separable with partition T 1, T 2 T, T 1 T 2 =, T 1 T 2 = T satisfying (19). Then we can derive (9) by adding the two rank inequalities (9) for T 1, T 2 : x(t 1 ) + x(t 2 ) + (α kj 1 +1,k j (T 1 ) + α kj 1 +1,k j (T 2 )) y kj r(t ). =x(t ) =α kj 1 +1,k j (T ) So the extended rank inequality associated with T does not define a facet of P Ē,K. However, not all (Ē, K)-non-separable sets T E lead to facet defining rank inequalities (9). This may depend on the structure of K as well. If, for instance, k j 1, k j K with k j 1 + 1 = k j, then the inequality x(e kj ) 1 can be derived by adding y kj 1 + y kj 0 and y kj 1 + x(e kj ) y kj 1 (see also property (P1 k ) to be defined below for the case K = {k}). In the following we will fully characterise the facetial structure of P Ē,k. We start with some easy results that follow directly from the definition of (Ē, {k})-separability. Observation 26. Let T E be a closed set with T 2. 1. If T Ē, then T is (Ē, {k})-separable. 2. If α 1,k (T \ Ē) = 0 and Ē T, then T is (Ē, {k})-separable. In addition to closedness and (Ē, {k})-non-separability it will turn out that a set T E leading to facet defining extended rank inequalities (14) has to satisfy the following properties: Definition 27. Let T E be a closed and properties (Ē, {k})-non-separable set. We say it satisfies (P1 k ) if Ē T and if T {e}, e Ē, in the case Ē = 2, (P2 k ) if Ē T and α 1,k(T \ Ē) 2. These properties are necessary conditions that a set T E leads to a facet defining rank inequality (14) of P Ē,k. 13

Observation 28. Let T E be a closed and (Ē, {k})-non-separable set. If the rank inequality associated with T is facet defining for P Ē,k, then it satisfies either (P1k ) or (P2 k ). Next we will prove that for a set T E closedness, (Ē, {k})-non-separability and properties (P1 k ), (P2 k ) lead to facetness of the associated inequality (14). In the proof of this result we need the following observation: Observation 29. Let T E be a closed set with Ē T. If α 1,k(T \ Ē) 2, then there exists a basis B of T such that B Ē Ē 2. Lemma 30. Let T E be a closed and (Ē, {k})-non-separable set satisfying properties (P1k ) or (P2 k ). Then the extended rank inequality (14) defines a facet of P Ē,k. Proof. Let T E be a closed, (Ē, {k})-non-separable set satisfying (P1k ) or (P2 k ). Further we assume that a T x + āy b is a facet defining inequality of P Ē,k so that at χ S + āy(s) = b whenever S is a root of x(t ) + α 1,k (T )y r(t ), the extended rank constraint associated with T. We will show that in this case a T x + āy b is a positive multiple of x(t ) + α 1,k (T )y r(t ). Our proof consists of four main steps. Claim 1. Under the given assumptions the constraints (14) are not implied by the constraints of the standard linearisation (see Observation 21). Proof of Claim 1. In order to prove that x(t ) + α 1,k (T )y r(t ) is not implied by (11) (13) we can easily determine in each of the cases an independent set J I such that J is a root of the rank inequality (14), but J is not a root of one of (11) (13). Technical details can be found in the Appendix. So by Lemma 12 we know that a 0, ā 0. Next we prove the following: Claim 2. A coefficient a(e), e E, or ā is zero if the corresponding coefficient in x(t ) + α 1,k (T )y r(t ) is zero. Proof of Claim 2. Assume, for a contradiction, a(e) > 0 for some e / T. Let S be a basis of T. Then S is a root of x(t ) + α 1,k (T )y r(t ) and consequently of a T x + āy b, too. Because T is closed S + e I, but S + e violates a T x + āy b, a contradiction. Now assume the coefficient α 1,k (T ) of y is zero in x(t ) + α 1,k (T )y r(t ), i. e., r(t ) + Ē \ T = r(t Ē) (20) by (A4). First, if Ē T, we choose a basis B of T with (T Ē) B. Then by (20) we know B + (Ē \ T ) I and this independent set is also a root of x(t ) + α 1,k(T )y r(t ). However, because y(b +(Ē \T )) = 1 it follows ā = 0. Second, let Ē T. Let B be a basis of T with Ē B. Furthermore, by assumption (P2 k ) we may apply Observation 29, and find a basis B of T with B Ē Ē 2. So there exist two elements e i, e j ( B \ B) Ē, i j. Because is a matroid there exists an f i B \ B B \ Ē with B e i + f i I and an f j B \ ( B e i + f i ) B \ B Ē with e i e j + f i + f j I. Note that all three sets B, B ei + f i, B ei e j + f i + f j are bases of T and define roots of x(t ) + α 1,k (T )y r(t ). In particular, this implies a(e i ) + ā = a(f i ) and a(e j ) = a(f j ). Submodularity gives r( B e j ) + r( B e j + f i + f j ) r( B e j + f i ) + r( B e j + f j ), =r(t ) 1 =r(t ) so at least one of B ej + f i, B ej + f j is a basis of T. In the second case we conclude a(e i ) = a(f i ), hence ā = 0. Otherwise, if B e j + f j is not a basis, a(e i ) = a(f j ) = a(e j ) and again by submodularity r( B e i e j + f j ) + r( B + f j ) r( B e i + f j ) + r( B e j + f j ), =r(t ) 1 =r(t ) =r(t ) 1 14

hence B e i + f j is a basis and we get a(f i ) = a(f j ) = a(e i ) and again ā = 0. This proves Claim 2. It remains to prove that all coefficients a(e), e T, have the same value and that the coefficient ā of y has the correct value. We start with the first statement. Claim 3. All coefficients a(e), e T, have the same value. Proof of Claim 3. We define the simple undirected graph G = (T, A) with set of nodes T and set of edges A := {{e, f}: B T, e, f T \ B, B + e, B + f bases of T } { } {e, f}: B T Ē, e, f T \ B, Ē B, B + e, B + f bases of T Ē. We prove the statement in two steps. Claim 3.1. If {e, f} A, then we have a(e) = a(f). Proof of Claim 3.1. Let {e, f} A. If, on the one hand, there exists B T, e, f T \ B, with B + e, B + f bases of T, then both bases define roots of x(t ) + α 1,k (T )y r(t ) and consequently e B+e a(e )+āy(b +e) = e B+f a(e )+āy(b +f), hence a(e)+āy(b +e) = a(f)+āy(b +f). If Ē T, then y(b + e) = y(b + f) = 0, otherwise α 1,k(T ) = r(t ) + Ē \ T r(t Ē) = r(t ) + 0 r(t ) = 0 by (A4) implying ā = 0, so a(e) = a(f) follows. On the other hand, if there exists B T Ē, e, f T \ B, with Ē B and B + e, B + f bases of T Ē, then y(b + e) = y(b + f) = 1. Furthermore (B + e) T + α 1,k (T )y(b + e) = (B + e) T + r(t ) + Ē \ T r(t Ē) = (B + e) T + r(t ) + (B + e) \ T r(t Ē) = r(t Ē) + r(t ) r(t Ē) = r(t ) by (A4), so B + e defines a root of x(t ) + α 1,k (T )y r(t ) (and, analogously, B + f). As above we conclude a(e) = a(f) and Claim 3.1 follows. Therefore, if we can prove that G is connected, all coefficients a(e), e T, must be equal. Claim 3.2. Graph G is connected. Proof of Claim 3.2. Assume, for a contradiction, that G is not connected, i. e., there exist T 1, T 2 T, T 1, T 2, T 1 T 2 = and T 1 T 2 = T such that there do not exist edges between T 1 and T 2. First assume r(t ) < r(t 1 ) + r(t 2 ) and let B be a basis of T with r(t 1 ) = r(b T 1 ). Then r(t 2 ) > r(b T 2 ), so there exists an e T 2 \ B so that (B T 2 ) + e I. However, B + e / I (B basis), but because is a matroid there exists an f (T 1 B ) so that B + e f is a basis of T. Consequently, choosing B := B f proves {e, f} A, a contradiction. Secondly, if r(t ) = r(t 1 ) + r(t 2 ) then we know that Ē + r(t Ē) < r(t 1 Ē) + r(t 2 Ē) by (A4) and because T is (Ē, {k})-non-separable. Let B be a basis of T Ē with Ē B and r(b (T 1 Ē)) = r(t 1 Ē). Consequently, r(b (T 2 Ē)) = Ē + B (T 2 \ Ē) < r(t 1 Ē) + r(t 2 Ē) r(t Ē) + B (T 2 \ Ē) = B (T 1 Ē) + r(t 2 Ē) r(t Ē) + B (T 2 \ Ē) = B + r(t 2 Ē) r(t Ē) = r(t 2 Ē). Hence there exists an e T 2 \ B so that (B + e) (T 2 Ē) I, but B + e / I. Because is a matroid there exists an f T 1 such that B +e f is a basis of T (otherwise (B +e) (T 2 Ē) / I), too. Consequently, choosing B := B f proves {e, f} A, a contradiction. So G is connected and this proves Claim 3.2 and so Claim 3. Finally, it remains to prove: Claim 4. The coefficient ā of y has the correct value. Proof of Claim 4. By Claim 2 we know that ā = 0 if α 1,k (T ) = 0. So, e. g., in the cases T Ē or Ē T the statement follows. So assume T Ē, Ē T and α 1,k (T ) > 0. Let B be a 15

basis of T Ē with r(b (T Ē)) = r(t Ē) and r(b T ) = r(t ). Then there exists a set S B \ Ē so that B := (B \ S) (Ē \ B) is also a basis of T Ē. For this note that S = Ē \ B = Ē \ T B \ T = B T + Ē \ T B = r(t ) + Ē \ T r(t Ē) = α 1,k(T ) > 0. In particular, y(b) = 0 because B does not contain all elements of Ē. Thus, with ã := a(e) being the value of all coefficients a(e), e T, we get that B, B are roots of x(t ) + α 1,k (T )y r(t ) (thus also of a T x + āy b), e B T a(e) + ā y(b) = =0 e B T a(e) + ā y(b ), =1 so ā = S ã = α 1,k (T )ã. Consequently Claim 4 follows. In total, this proves that a T x + āy b is a positive multiple (with factor ã) of x(t ) + α 1,k (T )y r(t ). The following theorem summarises the results on the facetial structure of P Ē,k. Theorem 31. All facet defining inequalities of P Ē,k E \ Ē, y 0 and all extended rank inequalities (14) 4 Separation x(t ) + α 1,k (T )y r(t ), are given by (11) and (12), x(e) 0, e T E, T closed, (Ē, {k})-non-separable and T satisfies property (P1 k ) or (P2 k ). By the well-known optimisation equals separation result [15] we know that the separation problem for the polytope P Ē,K can be solved in polynomial time because the corresponding optimisation problem can be solved in polynomial time. But the proof of the next result shows how to separate the extended rank inequalities more directly. Note, all inequalities of the standard linearisation (3)-(7) as well as the non-negativity constraints (8) can be separated in polynomial time by complete enumeration. Theorem 32. Let K = {k 1,..., k l } {2,..., k}. Given a point ( x, ȳ) the separation problem for the extended rank inequalities (9) of P Ē,K can be solved in (strongly) polynomial time if ȳ kj ȳ kj+1, j = 1,..., l 1, as well as ȳ k1 1, ȳ kl 0. Proof. Let ( x, ȳ) be a point that satisfies all requirements. Then we want to find a set T E such that (9) evaluated for ( x, ȳ) is violated or we want to prove that there does not exist such a set T. Therefore we define the function d: P(E) R, d(t ) = r(t ) x(t ) α kj 1 +1,k j (T )ȳ kj. Our separation problem is then equivalent to the problem to decide whether the optimal value of min d(t ) T E is negative or not. Using (A4) we can write d(t ) as: l 1 d(t ) = r(t ) (1 ȳ k1 ) + (r(t Ēk j ) (ȳ kj ȳ kj+1 )) + r(t Ēk l ) ȳ kl 0 0 0 x(t ) Ēk 1 \ T ȳ k1 ( {e m : m {k j 1 + 1,..., k j }} \ T ȳ kj ). j=2 16

f 1 f 2 e 1 e 2 e 2 f 2 e 1 f 1 e 3 e 4 f 1 e 1 e 2 e 4 e 3 f 1 f 2 f 3 e 1 e 2 (I) (II) (III) (IV) Figure 4: Visualisation of the four different graphs in Examples 33 35. So, d(t ) is submodular by the submodularity of the rank function r(.) and our separation problem reduces to minimising the submodular function d(t ), which can be done in strongly polynomial time [11, 18], e. g., using a combinatorial algorithm [22, 23]. 5 Extensions and Future Work A natural question is if it is possible to extend the previous results to related problems. For instance, Edmonds [7] also provided a complete description for the polytope that is associated with the intersection of two matroids. However, the following example shows that even in the case with exactly one additional quadratic monomial there exist facet defining inequalities a T x + āy b of the associated linearised polytope that are not a positive multiple of the standard linearisation but with a 0 and ā < 0. Thus, it might be necessary to extend the line of arguments in Section 2 significantly. Example 33. We consider the directed graph G = (V, E) in Fig. 4 (I) with four nodes and four arcs E = {e 1, e 2, f 1, f 2 } and the branching problem [6] (intersection of two matroids) with Ē = {e 1, e 2 }. Then one can check by direct computations that the inequality x(f 1 ) + x(f 2 ) + x(e 1 )+x(e 2 ) y 2 is facet defining for the associated linearised polytope although the coefficient of y e1 e 2 is 1. Note, the inequality above is of the form x(t ) + x(e 1 ) + x(e 2 ) y e1 e 2 r(t ), with T E, r(t ) = r(t + e 1 ) = r(t + e 2 ) = r(t + e 1 + e 2 ) 1 (r denotes the minimum of the rank of T for 1 and 2 ). Furthermore, considering the branching problem in Fig. 4 (II) the inequality 2x(e 1 ) + x(e 2 ) + x(f 1 ) + x(f 2 ) y e1 e 2 2 is facet defining. So, also coefficients of the x-variables larger than one may appear. This is a large contrast to the results by Edmonds for the polytope associated with the intersection of two matroids [7]. Indeed, given two matroids 1 = (E, I 1 ), 2 = (E, I 2 ) with rank functions r i : E N 0, i = 1, 2, the convex hull of all sets I E with I I 1, I I 2 can be described by the non-negativity constraints and the rank inequalities for the two single matroids [7]. Furthermore, it is interesting to extend the results to two or more monomials that are not necessarily nested. But again, this might lead to facet defining inequalities containing variables with negative coefficients apart from the standard linearisation. We start with an example where the two monomials are in pairwise conflict. Example 34. We consider the undirected graph G = (V, E) in Fig. 4 (III) with four nodes and five edges E = {f 1, e 1, e 2, e 3, e 4 } and the optimisation problem over the graphical matroid with the two additional monomials x(e 1 )x(e 2 ), x(e 3 )x(e 4 ). Then the inequality x(f 1 ) + y e1 e 2 + x(e 3 ) + x(e 4 ) y e3 e 4 2 is facet defining. So one of the coefficients of the linearised variables is negative although the inequality does not belong to the standard linearisation in this case. One might think that the problem in Example 34 arises because {e i : i = 1,..., 4} / I. But even if the union of all elements contained in the additional monomials is an independent set and none of the elements e E appears in more than one of these non-linear monomials, facet 17

defining inequalities can contain variables with negative coefficients apart from the standard linearisation. Example 35. We consider the undirected graph G = (V, E) in Fig. 4 (IV) with six nodes and set of edges E = {f 1, f 2, f 3, e 1, e 2, e 3, e 4 } and the optimisation problem over the graphical matroid with the two additional monomials x(e 1 )x(e 2 ), x(e 3 )x(e 4 ) and {e 1,..., e 4 } I. Then the inequality x(f 1 ) + x(f 2 ) + x(f 3 ) + y e1 e 2 + x(e 3 ) + x(e 4 ) y e3 e 4 4 is facet defining for the associated polytope. Let, in general, a matroid = (E, I) be given and let x(e 1 )x(e 2 ), x(e 3 )x(e 4 ), {e 1, e 2 } {e 3, e 4 } =, denote the two additional quadratic monomials with associated variables y e1 e 2, y e3 e 4. Then inequalities x(t ) + y e1 e 2 + x(e 3 ) + x(e 4 ) y e3 e 4 r(t ) + 1, T E \ {e 1, e 2, e 3, e 4 }, r(t + e 3 + e 4 ) > r(t ) + 1, (T + e 1 + e 2 + e i ) / I, i = 3, 4, are valid for the associated linearised polytope. The facet defining inequalities in Examples 34 and 35 are exactly of this type. Extending our results to problems with objective functions with more than one non-linear monomial, a good starting point might be to consider the three quadratic monomials x(e 1 )x(e 2 ), x(e 1 )x(e 3 ), x(e 2 )x(x 3 ) on the three elements e 1, e 2, e 3 E with {e 1, e 2, e 3 } I or {e 1, e 2, e 3 } forming a circuit. In both mentioned cases the monomials then fulfil some kind of up- and downwards compatibility, see [5]. Furthermore, it remains for future work to extend the results to (integral) polymatroids [7], see also [11, 23]. Here, the setting might be quite different because the variables are not necessarily binary. Appendix In the following we present several detailed, rather technical proofs of the results stated in the main part of this paper. Proof (of Observation 5). Let e cl(t ), then r(t + e) = r(t ). r(t + S) + r(t ) (R2) r(t + S + e) + r((t + S) (T + e)) (R3) r(t + S) + r(t + e) = r(t + S) + r(t ). Proof (of Observation 7). Let T E. (A1): If e i T, i {1,..., k}, then α i (T ) = 0 + r(t + Ēi) r(t + Ēi) = 0. (A2): Let i {1,..., k}. If r(t + Ēi) = r(t + Ēi 1), the statement is clear by {e i } \ T {0, 1}. If r(t + Ēi) = r(t + Ēi 1) + 1, this implies e i / T and so α i (T ) = 0. (A3): Let i {1,..., k} and j N, j > i. Then α j (T ) = {e j } \ T + r(t + Ēj 1) r(t + Ēj) = {e j } \ (T + e i ) + r(t + e i + Ēj 1) r(t + e i + Ēj) = α j (T + e i ). (A4): Let j {1,..., k}. Then α 1,j (T ) = j i=1 ( {e i} \ T + r(t + Ēi 1) r(t + Ēi)) = Ēj \ T + r(t + Ē0) r(t + Ēj) = Ēj \ T + r(t ) r(t + Ēj). (A5): Fix i {1,..., k}. If e T, the statement is clear. So let e (cl(t ) \ T ) \ {e i }. Then α i (T ) = {e i }\T +r(t +Ēi 1) r(t +Ēi) = {e i }\(T +e) +r(t +e+ēi 1) r(t +e+ēi) = α i (T + e) by Observation 5. (A6): Let i {1,..., k} and e i cl(t ) \ T. Then Observation 5 and (A1) show α i (T ) = {e i } \ T + r(t + Ēi 1) r(t + Ēi) = 1 + r(t + e i + Ēi 1) r(t + Ēi) = 1 = α i (T + e i ) + 1. 18

Proof (of Lemma 14). Let a T x + α T y b be a facet defining inequality of P Ē,K, K = {2,..., k}, satisfying the requirements. Let S i denote a root of a T x + α T y b such that Ēi S i and (i = k e i+1 / S i ). For m = 2 the desired statement follows directly from the fact that a T x + α T y b is not a positive multiple of (5). For m > 2 we know by (6) that there exists a root S of a T x + α T y b with e m / S and Ēm 1 S. So we may assume that µ := S Ēm 1 is maximal. (21) If S Ēm = Ēm 2, then S satisfies the requirements. Otherwise there must exist two elements e i, e j / S with i < j < m and Ēj 1 \ {e i } S. Let S j 1 be a root of a T x + α T y b with Ē j 1 S j 1 and e j / S j 1 (exists by (4) or by (7) in the case j = 2). Now assume in addition that S and S j 1 are chosen so that S Ēm 1 = µ and S S j 1 is maximal. (22) 1. If r(s) < r(s j 1 ), then there exists an e S j 1 \S so that e +S I. By Observation 13 we get ã(s +e, e) 0 ã(s j 1, e). If e e i, then e / Ēj and so ã(s j 1, e) = a(e) = ã(s +e, e) proving a(e) = 0. Thus S + e is also a root of a T x + α T y b contradicting (22). Otherwise we have ã(s j 1, e) = a(e) + α i +... + α j 1 = ã(s + e, e) = 0. As before S + e is also root of a T x + α T y b, a contradiction to (21). 2. If r(s) > r(s j 1 ), then there exists an e S so that S j 1 + e I. Certainly e / Ēj, so as above it follows a(e) = 0. So S, (S j 1 + e) contradict to assumption (22). 3. If r(s) = r(s j 1 ), then we know by Theorem 1 that there exists an e S \ S j 1 so that S e+e i, S j 1 +e e i I. Note, e / Ēj, which implies ã(s, e) = a(e) = ã(s j 1 +e e i, e) as well as ã(s j 1, e i ) = a(e i )+α i +...+α j 1 = ã(s e+e i, e i ). Applying Observation 13 twice then shows ã(s e + e i, e i ) ã(s, e) = ã(s j 1 + e e i, e) ã(s j 1, e i ) = ã(s e + e i, e i ) and so S e + e i is also a root of a T x + α T y b, contradicting (21). Proof (of Lemma 15). Let a T x + α T y b be a facet defining inequality of P Ē,K with K = {2,..., k} that is not a positive multiple of one of (3) (8). Let S i be defined as in the proof of Lemma 14. Assume, for a contradiction, the assertion is false. Then there must exist an m {2,..., k} and a minimal h {m,..., k} so that α m +... + α h < 0 and there exists a root S of a T x+α T y b satisfying Ēm \{e j, e m } = S Ēm, j < m, and h = k or e h+1 / S. Furthermore, by (4) there exists a root S h with Ēh S h and h = k or e h+1 / S. We may assume that S S h is maximal. (23) 1. If r(s) < r(s h ), then there exists an e S h \ S so that S + e I. We consider four cases depending on e and apply Observation 13. If e / Ēh, then ã(s + e, e) = a(e) and ã(s h, e) = a(e). Observation 13 shows that a(e) = 0, so S + e is also a root of a T x + α T y b, contradicting (23). If e = e m, then 0 ã(s + e m, e m ) = a(e m ) and 0 ã(s h, e m ) = a(e m ) + α m +... + α h, hence α m +... + α h 0, a contradiction. If e = e j, then we similarly get 0 ã(s + e j, e j ) = a(e j ) + α j +... + α m 1 and 0 ã(s h, e j ) = a(e j ) + α j +... + α h, hence α m +... + α h 0, a contradiction. If e = e i, m < i h, then 0 ã(s + e i, e i ) = a(e i ) and 0 ã(s h, e i ) = a(e i ) + α i +... + α h, hence we derive α i +... + α h 0. Because e i / S we know by minimality of h that α m +... + α i 1 0. Together we have α m +... + α h 0, a contradiction. 19