On the positive semidenite polytope rank

Size: px
Start display at page:

Download "On the positive semidenite polytope rank"

Transcription

1 On the positive semidenite polytope rank Davíd Trieb Bachelor Thesis Betreuer: Tim Netzer Institut für Mathematik Universität Innsbruck February 16, 017

2 On the positive semidefinite polytope rank - Introduction 1 Introduction Note that this bachelor thesis largely follows []. If you want to do linear optimization over a polytope, the complexity of many algorithms depends on the size of the representation of the polytope. So, if you have a (complicated) polytope, the idea is to nd a simpler convex set of higher dimension which has the polytope as a linear image of it and then optimize over that instead. This motivates following denition: Denition 1.1. Let P R n be a polytope. For a closed convex cone C R m and an ane space L R m, C L is called a C-lift of P, if there exists a linear map π : R m R n such that P = π(c L). In this bachelor thesis, we will be interested in the cases C = R k + = {x R k x 1 0,..., x n 0} and especially C = S k + = {M R k k M symmetric, positive semidenite}. Note that R k + embeds into S k + via diagonal matrices. Now one is interested in nding the smallest cone in the families {R k +} and {S k +}, respectively, which allows a lift of a polytope P : Denition 1.. Let P R n be a polytope. (1) The nonnegative rank of a polytope is given by rank + (P ) := min({k N P has R k +-lift}). () The positive semidenite (psd) rank of a polytope is given by rank psd (P ) := min({k N P has S k +-lift}). In order to investigate the properties of the psd rank of a polytope, the following denitions are also needed: Denition 1.3. Let P R n be a n-dimensional polytope with vertex set {p 1,..., p v } and irredundant inequality representation P = {x R n b 1 a 1, x 0,..., b f a f, x 0}, where b j R and a j R n for j = 1,..., f. Then the slack matrix S P R v f is given by (S P ) ij = b j a j, p i. Short reminder: the polar dual of a cone C is given by C := {y R m x, y 0 x C}. Note that both S k + and R k + are self dual cones, i.e. C = C. For symmetric matrices A, B we will use the trace inner product A, B = Tr(AB). Denition 1.4. Let M R p q + be a nonnegative matrix and C a closed convex cone with its polar dual C. David Trieb

3 On the positive semidefinite polytope rank - Introduction (1) A pair of ordered sets a 1,..., a p C and b 1,..., b q C is called C-factorisation of M, if M ij = a i, b j for all i = 1,..., p and j = 1,..., q. () For C = S k + (respectively, R k +), a C-factorisation of M is called a psd (respectively, nonnegative) factorisation of M. (3) rank psd M := min({k N M has S k +-factorisation}) denotes the psd rank of M. rank + M := min({k N M has R k +-factorisation}) denotes the nonnegative rank of M. Note that in Denition 1.3 any positive scaling of a facet inequality of P can be used and therefore, if you positively scale the columns of a slack matrix S P of P, you get a slack matrix of P again. In the following, any such slack matrix is denoted by S P. Remark 1.5. These are some results we will use: rank psd P = rank psd S P and rank + P = rank + S P. [1], [4] rank + P and rank psd P are invariant under ane transformations. If P is the polar polytope of P and S P a slack matrix of P, then S P scaling) a slack matrix of P and therefore, rank psd P = rank psd P. is (up to row Since R k + embeds into S+ k via diagonal matrices it follows rank psd M rank + M for M R p q +. Together with the rst statement of this remark this yields rank psd P rank + P for any polytope P R n. Proof. We prove that rank + P = rank + S P. rank psd P = rank psd S P can be proven similarly. (1) rank + P rank + S P : Let P = {x R n 1 c 1, x 0,..., 1 c f, x 0} R n be a fulldimensional polytope with vertices p 1,..., p v and rank + P = k. Then there exists a linear map π : R k R n such that P = π(r k + L), where L is an ane subspace. Suppose that L = z + L 0, where L 0 is a linear subspace. Then P is given by P = {x R n x = π(w), w R k + (z + L 0 )}. Let π : R n R k be the adjoint of π and assume that R k + L 0 = {0} (we may do so, because P is bounded). Then, by strong conic duality: P = {x R n y π (x) R k +, y L 0, y, z = 1} Since y L 0, we have w i, y = 1 for any w i L 0 +z = L. Now dene a i := w i and b j := y π (c j ), where w i π 1 (p i ) R k + and y L 0 (Rk + + π (c j )) with y, z = 1. Such w i and y do exist and we get (S P ) ij = 1 c j, p i = 1 c j, π(w i ) = 1 π (c j ), w i = 1 w i, y b j = = w i, b j = a i, b j. As a i R k + and b j R k +, we have rank + S P k = rank + P. David Trieb 3

4 On the positive semidefinite polytope rank - Introduction () rank + P rank + S P : Let P be a polytope as above. Assume that rank + S P = k. Then there exist a 1,..., a v, b 1,..., b f R k + such that (S P ) ij = a i, b j. Dene L := {(x, y) R n R k 1 x, c i = y, b j, j = 1,..., f}. Let L k R k denote the projection onto the second coordinate of L. Note that a i R k + L k and 0 / L k (0 L k x : 1 x, c i = 0 for all i = 1,..., f P is contained in the hyperplane {y x, y = 1}, which is a contradiction, because 0 P ). For all x R n for which there exist y R k + such that (x, y) L we have 1 x, c i 0, i = 1,..., f 1 x, y 0 y P x (P ) = P. Claim: y R k + L k :! x y R n : (x y, y) L. Proof of claim: The existence follows from the denition of L k. Now let x y and x y be two such points. Since (x y, y), (x y, y) L, we get 1 x y, c i = 1 x y, c i = y, b j, and thus we obtain for t R: 1 tx y + (1 t)x y, c i = 1 t x y, c i x y, c i + t x y, c i = y, b j, i.e., (tx y + (1 t)x y, y) L, which means that the line through x y and x y is contained in P, which is a contradiction, unless x y = x y. Therefore, the map R k + L k R n : y x y is well-dened. As this map is ane and 0 / L k, we may extend it to a linear map π : R k R n. We have already proven π(r k + L k ) P above. Since for all i = 1,..., v, a i R k + L k, it follows p i = π(a i ) π(r k + L k ) and because P is the convex hull of its vertices and π(r k + L k ) is convex, we obtain C π(r k + L k ). Thus C = π(r k + L k ) and we found a R k +-lift of P. Hence, rank + P rank + S P. Let P R n be a n-dimensional polytope with rank psd P = k. Then there exists a linear map π : S+ k R n such that P = π(s+ k L), where L S+ k is an ane subspace. Let π 1 : R n R n be an ane transformation. Then π 1 (P ) = z + π (π(s+ L)) k for some z R n and a linear function π. As the composition π π is also linear, this means that the polytope π 1 (P ) z has a S+-lift. k Since π 1 (P ) z can be mapped to P via another ane transformation, it follows rank psd π 1 (P ) z = k. From the rst statement of the remark, we know that rank psd π 1 (P ) z = = rank psd S π1 (P ) z. Now π 1 (P ) z has the same slack matrix as π 1 (P ) (see Lemma 1.6 below) and we get k = rank psd π 1 (P ) z = rank psd S π1 (P ) z = rank psd S π1 (P ) = rank psd π 1 (P ). The proof for rank + P works analogously. Lemma 1.6. Let P R n be a polytope with vertices p 1,..., p v and inequality representation P = {x R n b 1 a 1, x 0,..., b f a f, x 0}, and let z R n. Then S P = S P +z. David Trieb 4

5 On the positive semidefinite polytope rank - Hadamard square roots Proof. The polytope P + z has vertices p 1 + z,..., p v + z and it holds P + z = {x R n b 1 a 1, x + a 1, z 0,..., b f a f, x + a f, z 0}. Therefore, we obtain (S P +z ) ij = b j a j, p i + z + a j, z = b j a j, p i a j, z + a j, z = (S P ) ij. The following proposition compares psd rank with usual rank: Proposition 1.7. [3] Let M R p q +. Then rank M rank psd M min({p, q}). Proof. (1) Let rank psd M = k and assume that A 1,..., A p, B 1,..., B q give an S+-factorisation k of M. For N S+ k dene ( vec(n) := N 11, N,..., N kk, N 1, N 13,..., N 1k, N3,..., N k,..., N 34,..., ) N (k 1)k. Then vec(n) R ( k+1 ) and vec(a i ), vec(b j ) = = k m=1 n=1 k l=1 Thus rank M ( ) k+1 and it follows rank M (k + 1)k k 1 A i ll Bj ll + k m=1 n=m+1 A i mnb j mn = k A i mnbmn j = Tr(A i B j ) = A i, B j = M ij. = 0 k + k rank M = k rank M. () Let e j denote the jth standard unit vector. Since diag(m i ), diag(e j ) = = Tr (diag(m i ) diag(e j )) = M i, e j = M ij, there exists an S q + -factorisation of M. Similarly, one gets a S p + -factorisation of M. Hadamard square roots In this section we will use Hadamard square roots to study the psd rank of nonnegative matrices. Denition.1. Any matrix whose (i, j)-entry is a square root of the (i, j)-entry of M is called a Hadamard square root of M and is denoted by M. In particular, let + M be the all-nonnegative Hadamard square root of M. The square root rank of M is dened as rank M := min({rank M}). David Trieb 5

6 On the positive semidefinite polytope rank - Hadamard square roots Proposition.. Let M R p q + be a nonnegative matrix. Then rank psd M rank M. Proof. If M is a Hadamard square root ( of M of rank r, then there exist vectors M ) a 1,..., a p and b 1,..., b q R r such that = a i, b j. We obtain ij ( r ) r r M ij = a i, b j = (a i b j) = a ik b jk = a ik a il b jk b jl = k=1 k=1 l=1 a i1 a i1 a i a i1 a ir b a i1 a i a j1 b j1 b j b j1 b jr i a i a ir b j1 b j b j b j b jr = Tr = a i1 a ir a i a ir a ir b j1 b jr b j b jr b jr ( ) = Tr a i a i b jb j = a i a i, b jb j and thus get a S r +-factorisation of M. Hence, rank psd M r. Example.3. This simple example shows that the upper bound in. can be strict: Consider M := It holds rank M = rank M = 3, but rank psd M =, as the follwing S+-factorisation of M shows: ( ) ( ) ( ) A :=, A :=, A :=, ( ) ( ) ( ) B 1 0 :=, B 0 0 :=, B 3 1 := One can easily check that M ij = A i, B j. Lemma.4. k = rank M is the smallest k that admits a S+- k factorisation for M R p q + where all factors have rank one. Proof. (1) If k = rank M, then there exists a Hadamard square root M of M such that rank M = k and the proof of proposition. gives a S k +-factorisation of M where all factors have rank one. () k rank M: Let a 1 a 1,..., a pa p S k + and b 1 b 1,..., b qb q S k + be a S k +-factorisation of M with rank one factors. Then M ij = a i a i, b jb j = a i, b j and the matrix M with ( M ) ij = a i, b j is a Hadamard square root of M with rank M k. David Trieb 6

7 On the positive semidefinite polytope rank - Hadamard square roots Next, we show a method to increase the psd rank of a matrix by one. Proposition.5. ( ) Let M R p q + with rank psd M = k, w R q M 0 +, α > 0 and M :=. w α Then rank psd M = k + 1. Furthermore, in any S+ k+1 -factorisation of M, the factor associated to the last column of M is of rank one. Proof. (1) First, we show that rank psd M cannot be smaller than k + 1: Assume that there is a S k +-factorisation of M with factors A 1,..., A p, A S k + associated to its rows and B 1,..., B q, B S k + associated to its columns. From A, B = α 0 we get A 0 B. Let r = rank B > 0. Then there exists an orthogonal matrix U such that U 1 BU = diag(λ 1,..., λ r, 0,..., 0) =: D, where λ 1,..., λ r are the positive eigenvalues of B. Let A i := U 1 A i U S k + i = 1,..., p. Then i = 1,..., p: D, A i = U 1 BU, U 1 A i U = Tr(U 1 BA i U) = Tr(BA i ) = B, A i = 0. Because the diagonal entries of A i are nonnegative and those of D positive, the rst r diagonal entries of A i must be zero, since Tr(DA i ) = 0. Since A i is psd, it follows that the rst r rows and columns of A i are all zero (if this were not the case, there would be at least one negative eigenvalue of A i ). For B j := U 1 B j U, j = 1,..., q we have A i, B j = Tr(U 1 A i UU 1 B j U) = Tr(U 1 A i B j U) = Tr(A i B j ) = A i, B j = M ij i = 1,..., p and j = 1,..., q. Since A i only has nonzero entries in its bottom right (k r) (k r)-block, it follows that M ij = A i, B j, where A i and B j are the bottom right (k r) (k r)-submatrices of A i and B j respectively. Now we have found a S+ k r -factorisation of M which is a contradiction to rank psd M = k. Therefore, rank psd M k + 1. () If A 1,..., A p, B 1,..., B q S+ k is an S+-factorisation k of M, then A 1,..., A p, A, B 1,..., B q, B S+ k+1 is an S+ k+1 -factorisation of M, where ( ) ( ) ( ) ( ) A Ai 0 i :=, B j Bj 0 :=, A 0 0 :=, B 0 0 := w j α This proofs that rank psd M = k + 1. (3) Now let B be the matrix associated to the last column of M in a S+ k+1 -factorisation of M. Since α > 0, B 0. So, let r = rank B > 0. By applying the same arguments as above, we get a S+ k+1 r -factorisation of M. Because of rank psd M = k and r > 0, r = rank B must be one. David Trieb 7

8 On the positive semidefinite polytope rank - Psd ranks of polytopes Example.6. If M R n n + is a diagonal matrix with positive entries, then rank psd M = n. Furthermore, each factor of any S+-factorisation n of M is of rank one. Proof. The proof is by induction on n. For n = 1, the rst statement is trivially true. If rank psd M = n for a diagonal matrix M R+ n n with positive entries, then rank psd M = ( ) M 0 n + 1, where M = and α > 0 by Proposition.5. That each factor of a 0 α S+-factorisation n of M must have rank one follows from the second part of Proposition.5 applied to M and M. 3 Psd ranks of polytopes For this section, let P R n be a n-dimensional polytope. Lemma 3.1. rank S P = n + 1. Proof. If the vertices of P are p 1,..., p v and P = {x R n b 1 a 1, x 0,..., b f a f, x 0} is an irredundant inequality representation, then (S P ) ij = b j a j, p i and 1 p 1 ( ) b1 b S P = f... 1 p a 1 a f v Call the rst factor F and the second one G. Since P is a n-dimensional polytope, there exist n+1 vertices of P such that no vertex lies in the linear span of the other n vertices, implying that rank F n + 1 and thus, because F R v (n+1), rank F = n + 1. Now assume that rank G n. Then there exists a nonzero c R n+1 such that c G = (0,..., 0). If c 1 0, then normalise c such that c 1 = 1 and set c := (c,..., c n+1 ). Then, for all j: b j a j, c = 0, implying that c lies on every facet of P, which is a contradiction. If c 1 = 0, it follows a j, c = 0 (c as above) for all j = 1,..., f. This means that c (as a vector) is parallel to every facet of P, which is a contradiction, since P is fulldimensional. Thus, rank G = n + 1. As both factors have rank n + 1, it also holds rank S P = n + 1. Proposition 3.. rank psd P n + 1. If equality holds, every S+ n+1 -factorisation of S P has only rank one factors. Proof. The proof is by induction on the dimension n. For n = 1, P is a line segment with vertices p 1, p and facets f 1, f, where p 1 = f and p = f 1. Because of (S P ) 1 = b 1 a 1 p = b 1 a 1 f 1 = 0, and, analogously, (S P ) 1 = 0, the slack matrix S P of P is a diagonal matrix with positive entries. From Example.6 we get that rank psd P = rank psd S P = and any S+-factorisation only uses rank one factors. David Trieb 8

9 On the positive semidefinite polytope rank - Psd ranks of polytopes Now suppose the rst statement holds up to dimension n 1. Let P R n be a fulldimensional polytope. Let F be a facet of P with vertices p 1,..., p s and facets f 1,..., f t and slack matrix S F. Assume that f j corresponds to the facet F j of P for j = 1,..., t. Then, by induction hypothesis, rank psd F n. Suppose that, if p / F is a vertex of P, the top left (s + 1) (t + 1) submatrix of S P is indexed by p 1,..., p s, p in the rows and F 1,..., F t, F in the columns. We will call that submatrix S F. S F has the form ( ) S F SF 0 =, α where α > 0. From Proposition.5 we get that rank psd S P rank psd S F n + 1. If rank psd P = n + 1, then there exists a S+ n+1 -factorisation of S P, and therefore of S F. By Proposition.5, the factor corresponding to the facet F must be of rank one. Applying this argumentation to all facets F of P, it follows that all factors indexed by the facets of P have rank one. Recall that S P is, up to row scaling, a slack matrix of the polar dual polytope P and therefore, all factors corresponding to the vertices of P also have rank one. Corollary 3.3. If a facet of a polytope P has psd rank k, then rank psd P k + 1. Especially, if P R n is a n-dimensional polytope with psd rank n + 1, then, for every i-dimensional face F of P, rank psd F = i + 1. Theorem 3.4. rank psd P = n + 1 if and only if rank S P = n + 1. Proof. By remark 1.5, rank psd P = rank psd S P. If rank S P = n+1, then, by Proposition., rank psd P = rank psd S P rank S P = n + 1 and rank psd P n + 1 (Proposition 3.) implies rank psd = n + 1. On the other hand, if rank psd P = n + 1, then there exists a S+ n+1 -factorisation of S P, whose factors are by Proposition 3. all of rank one. By Lemma.4, rank S P n + 1. Since rank psd S P rank S P, rank S P = n + 1. Example 3.5. Consider the pentagon P with vertices p 1 = (0, 0), p = (1, 0), p 3 = (, 1), p 4 = (1, ), p 5 = (0, 1) and inequality representation { ( ) ( ) ( ) x P = R 0 x, y 1 y ( ) ( ) ( ) 1 x 1 3, 0, 1, 1 y 1 ( 1 0, 1 1 ( ) x 0, y ( 1 0 ) ( ) x, 0, y ) ( ) x, y } 0. David Trieb 9

10 On the positive semidefinite polytope rank - Psd ranks of polytopes Then we obtain the slack matrix S P = (the inequalities have been multiplied by, 4, 4, 4, and, respectively). Theorem 4.5 will show that the psd rank of P is at least four. Due to the following S+-factorisation 4 of P we have rank psd P = 4: It holds (S P ) ij = A i, B j with A i, B j as follows: A 1 := , A := , A3 := A 4 := , A5 := , B 1 := , B := , B3 := , B 4 := , B5 := Now let us look at rank S P. Let 0 a b c d e f S := g 0 0 h i j k 0 0 l. m n o 0 0 Then there exists a Hadamard square root of S P with rank 4 if and only if there exists a solution to this system of equations: det(s) = 0, a = 4, b = 1, c = 4,..., o = 8. Using a computer one can easily check that this system of equations has no solution and it follows rank S P = 5. Thus, the psd rank of a polytope may be smaller than the square root rank of its slack matrix., David Trieb 10

11 On the positive semidefinite polytope rank - Psd ranks of polytopes Example 3.6. Consider a pentagonal pyramid P P R 3 with the pentagon P from Example 3.5 as the base. The vertices of P P are p 1 = 0, p = 0, p 3 = 1, p 4 =, p 5 = 1, p 6 = and it holds P P = {(x, y, z) R 3 y z 0, 1 x + y z 0, 3 x y z 0, 1 + x y z 0, x z 0, z 0}. Since P is a facet of P P and rank psd P = 4 (see above), we get from Corollary 3.3 that rank psd P P 5. A slack matrix of P P is given by S P P = The top left 5 5 submatrix of S P P is equal to S P and thus, by Proposition.5, rank psd P P = rank psd S P P = 5. Example 3.7. Let H R be a regular hexagon with vertices ( ) ( ) ( ) p 1 =, p 0 = 3, p 3 = 3, p 4 = ( ) ( ) 1 1, p 0 5 = 3, p 6 = ( 1 3 and inequality representation { ( ) ( ) ( ) ( ) ( ) x H = R y / x 0, 0, 3 y 4/ x, 0, 3 y ( ) ( ) ( ) ( ) / x, 0, 3 y / x, 0, 3 y ( ) ( ) ( ) ( ) } 0 4/ x, 0, 3 y / x, 0. 3 y The slack matrix S H of H is then given by S H = ) David Trieb 11

12 On the positive semidefinite polytope rank - Polytopes of minimum psd rank It holds rank + S H = 5, while the following Hadamard square root of H has rank four: This example shows that it is not enough to compute rank + S P rank S P. in order to get Example 3.8. Let H be as in Example 3.7. In Proposition 3. and Theorem 3.4 we have shown that if P is an n-dimensional polytope with rank psd P = n + 1, then rank psd P = rank P and all S+ n+1 -factorisations of S P have only rank one factors. From Theorem 4.5 we have that rank psd H 4, hence we get from Theorem 3.4 that the square root rank of S H is at least four and since we have seen a Hadamard square root of S H of rank four in the example above, it holds rank S H = 4. We also get from Theorem 4.5 that rank psd H 4 and the following S+-factorisation 4 proofs that rank psd H = 4: A 1 := , A := , A3 := , A 4 := , A5 := , A6 := , B 1 := , B := , B3 := , B 4 := , B5 := , B6 := Note that rank A i = for all i = 1,..., 6. Thus, there can be S rank psd P + -factorisations of a slack matrix S P of a polytope P, where there are factors with rank greater than one even if rank psd P = rank P. David Trieb 1

13 On the positive semidefinite polytope rank - Polytopes of minimum psd rank 4 Polytopes of minimum psd rank For a n-dimensional polytope P R n, rank + P n + 1 and the only n-dimensional polytopes that achieve the lower bound of n + 1 are simplices. In this section, we will see that there are several classes of polytopes that achieve the minimum psd rank of n + 1. Denition 4.1. A full fulldimensional polytope P R n is called -level if its slack matrixes entries are all 0 or 1. Remark 4.. A polytope P is -level if and only if for every facet F of P, all vertices of P lie on the union of the facet F and one hyperplane parallel to the hyperplane containing F. For example, any n-dimensional cube is -level. Corollary 4.3. If P is a n-dimensional -level polytope, then rank psd P = n + 1 and every S+ n+1 - factorisation of P only has factors of rank one. Proof. Let S P be a slack matrix of a n-dimensional -level polytope P. We have rank S P rank + S P = rank S P = n+1 and therefore, by Proposition., rank psd S P rank S P n+1 and thus rank psd P = rank psd S P = n+1. The second statement follows from Proposition 3.. Theorem 4.4. For any fulldimensional polytope P R n with n + vertices we have rank psd P = n + 1. Proof. Let P be a polytope with n + vertices and f facets. Then S P R (n+) f + and rank S P = n + 1. Since rank S P = n + 1, there exist a i R (not all of them equal to zero) such that n+ i=1 a i(s P ) i = (0,..., 0). Because every vertex lies on at least n dierent facets, each column of S P has at least n zeros. Thus, considering above equation componentwise, we have that for all j = 1,..., f there exist i 1, i {1,..., n + } such that a i1 (S P ) i1 j + a i (S P ) i j = 0. Now let b i := sgn(a i ) a i. Then b i1 (SP ) i1 j + b i (SP ) i j = 0. Since this is true for all components, we get that n+ i=1 b i (SP ) i = (0,..., 0) and thus rank + S P = n + 1 = rank S P. This together with Theorem 3.4 proves the statement. Theorem 4.5. Let P R be a -dimensional polytope. Then rank psd P = 3 P has at most four vertices. Proof. " ": Follows from 4.4. " ": Let P R be a polytope with at least 5 vertices. Since the psd rank is invariant under ane transformations we can assume that two facets of P are given by f 1 = {(x, y) R x 0} and f = {(x, y) R y 0} with vertices v 1 = (0, 0), v = (0, 1) and v 3 = (1, 0). Furthermore, let v 4 = (a, b) and v 5 = (c, d) be the vertices sharing an edge with (0, 1) and (1, 0), respectively. These two edges are given David Trieb 13

14 On the positive semidefinite polytope rank - Polytopes of minimum psd rank by f 3 = {(x, y) R y b 1 a x + 1} = {(x, y) R a + (b 1)x ay 0} and f 4 = {(x, y) R y d c 1 (x 1)} = {(x, y) R (c 1)y d(x 1) 0}. Then, the 5 4 submatrix S P of S P indexed by these vertices and facets is given by 0 0 a d S P c 1 + d = 1 0 a + b 1 0 a b 0 cb b da + d. c d bc c ad + a 0 We now show that rank S P 4. Theorem 3.4 then yields rank psd P > 3, which proves the statement. It is enough to show that any Hadamard square root H of the top left 4 4 submatrix of S P has rank four. Assume that H = 0 0 ± a ± d 0 ±1 0 ± c 1 + d ±1 0 ± a + b 1 0 ± a ± b 0 ± cb b da + d has rank three. It is obvious that the rst three rows are independent, hence we can write the fourth row as a linear combination of the rst three and it must hold H 4 = ± ah 3 ± bh a + b 1H 1 = ( = ± a, ± b, 0, d(a + b 1) ± ) b(d + c 1). Let α := b(d + c 1) and β := d(a + b 1). Thus H 44 = ± cb b da + d = ± α β! = ± α β, and this implies α = β bc b = da d a 1 = c 1. Hence, the slopes of the lines between (a, b) and (1, 0), and (c, d) and (1, 0) are the same, meaning that (a, b), (c, d) and (1, 0) are collinear and cannot be all vertices, except when (a, b) = (c, d). Lemma 4.6. Let P R 3 be a fulldimensional polytope with psd rank four and p be a vertex of P of degree four. Then the four faces incident to p are all triangles. Proof. Assume that one of the four faces is not a triangle. By Corollary 3.3, this face has psd rank 3 and hence, by Theorem 4.5, it is a quadrilateral. Thus, P contains the conguration b d p 1 f 3 p 5 f 1 p f 4 p 4 f p p 3 David Trieb 14

15 On the positive semidefinite polytope rank - Polytopes of minimum psd rank where p 1,..., p 5 are vertices of P. According to Lemma 3.1, rank S P = 4. Since rank psd P = 4, by Theorem 3.4, there exists a Hadamard square root S P of S P of rank four. Let f 1, f, f 3 and f 4 be the faces of P containing the vertices p, p 1, p ; p, p, p 3 ; p, p 1, p 5 and p, p 3, p 4, p 5, respectively and M be the submatrix of S P indexed by p, p 1, p, p 3, p 4 in the rows and f 1, f, f 3, f 4 in the columns. Then, by scaling the columns of S P accordingly, it holds M = a 1 0 b 0, c d e 0 where a, b, c, d, e > 0. Since the four rows of S P (and S P ) corresponding to the rst four rows of M are clearly independent, we can write the row of S P (and S P ) corresponding to the fth row of M as a linear combination of the other four. This yields the conditions d+ae = abc (for S P ) and d +a e = a b c (for S P ), implying d +a e = d +a e + ade and thus ade = 0, which is a contradiction to a, b, c, d, e > 0. Denition 4.7. Let F(P ) denote the set of all faces of a polytope P. Two polytopes Q and Q are called combinatorially equivalent, if there exists an isomorphism f : F(Q) F(Q ) which preserves the partial order given by inclusion. Remark 4.8. Roughly speaking, two polytopes that are combinatorially equivalent have the same number of faces of each dimension arranged in the same way. For example, every three-dimensional polytope with six two-dimensional quadrilateral faces is combinatorially equivalent to a cube. Theorem 4.9. A fulldimensional polytope P R 3 with rank psd P = 4 is combinatorially equivalent to an octahedron, bisimplex, simplex, quadrilateral pyramid, triangular prism, or cube. Proof. Let v, e, f be the number of vertices, edges, and faces of P, respectively. Further, v t := #{v v vertex of P, deg v = 3}, v q := #{v v vertex of P, deg v = 4}, f t := #{f f triangular face of P } and f q := #{f f quadrangular face of P }. By counting the edges on each face, e = 3f t + 4f q, and thus (by considering P ), we also have e = 3v t + 4v q. Since v e + f = (Euler's formula) and v = v t + v q, f = f t + f q, we get e = v t + v q + f t + f q 4 and thus: 3v t + 4v q = v t + v q + f t + f q 4 v t = f t + f q 4 v q, and analogously, f t = v t + v q 4 f q, from which we get that v t and f t are even and v t + f t = v t + f t 8 v t + f t = 8. Therefore, the only polytopes we have to consider are polytopes, where (v t, f t ) equals (0, 8), (, 6), (4, 4), (6, ) and (8, 0). (1) (v t, f t ) = (0, 8): In this case, every vertex of P has degree four. From Lemma 4.6 we get that all faces of P must be triangular. An octahedron is the only polytope in R 3 that satises these conditions. David Trieb 15

16 On the positive semidefinite polytope rank - Polytopes of minimum psd rank () (v t, f t ) = (, 6): Here, there must be at least one vertex p of P of degree four, and thus, by Lemma 4.6, P contains the following structure: p 1 p 4 p p p 3 (a) Since v t =, at least two of the vertices p 1, p, p 3, p 4 have degree four. Now assume, that two of those degree four vertices (let's say, p 1 and p ) are adjacent. Then, again by Lemma 4.6, p 1 and p are contained in four faces which are triangles, implying that p 1 and p each lie on two triangular faces not shown in the gure. Since p 1 already shares a face with p, they can only share at most one of the four extra triangular faces, meaning that there must be at least 7 triangular faces, which is a contradiction to f t = 6. Thus, only two out of p 1, p, p 3, p 4 have degree four and those two are nonadjacent. Again, both lie on two extra triangular faces not shown in the gure, and because of f t = 6, they must share both of them. Consequently, P is a bisimplex (3) (v t, f t ) = (4, 4): If v q = 0, then there are only four vertices in total and the polytope must be a simplex. If v q > 0, then, by Lemma 4.6, P contains structure (a). Since f t = 4, p i must have degree three, i = 1,, 3, 4 (if one of the p i 's had the degree four, then there would be too many triangular faces). Therefore, the polytope is a quadrilateral pyramid. (4) (v t, f t ) = (6, ): Then, for the polar dual of P we have (v t, f t ) = (, 6), hence, the polar dual of P is a bisimplex and thus, P is a combinatorial triangular prism. (5) (v t, f t ) = (8, 0): Again, for the polar dual of P : (v t, f t ) = (0, 8). Thus P is the polar dual of an octahedron, i.e. a cube. References [1] Pablo A. Parrilo João Gouveia and Rekha R. Thomas. Lifts of convex sets and cone factorisations. Mathematics of OR, 38:4864, 013. [] Richard Z. Robinson João Gouveia and Rekha R. Thomas. Polytopes of minimum positive semidenite rank, 01. [3] Richard Z. Robinson. The positive semidenite rank of matrices and polytopes, 014. [4] Mihalis Yannakakis. Expressing combinatorial optimization problems by linear programs. J. Comput. System Sci., 43:441466, David Trieb 16

POLYTOPES OF MINIMUM POSITIVE SEMIDEFINITE RANK

POLYTOPES OF MINIMUM POSITIVE SEMIDEFINITE RANK POLYTOPES OF MINIMUM POSITIVE SEMIDEFINITE RANK JOÃO GOUVEIA, RICHARD Z ROBINSON, AND REKHA R THOMAS Abstract The positive semidefinite (psd) rank of a polytope is the smallest k for which the cone of

More information

arxiv: v3 [math.co] 25 Jul 2016

arxiv: v3 [math.co] 25 Jul 2016 FOUR-DIMENSIONAL POLYTOPES OF MINIMUM POSITIVE SEMIDEFINITE RANK arxiv:150600187v3 [mathco] 25 Jul 2016 JOÃO GOUVEIA, KANSTANSTIN PASHKOVICH, RICHARD Z ROBINSON, AND REKHA R THOMAS Abstract The positive

More information

Lifts of convex sets and cone factorizations

Lifts of convex sets and cone factorizations Lifts of convex sets and cone factorizations João Gouveia Universidade de Coimbra 20 Dec 2012 - CORE - Université Catholique de Louvain with Pablo Parrilo (MIT) and Rekha Thomas (U.Washington) Lifts of

More information

MATH 890 HOMEWORK 2 DAVID MEREDITH

MATH 890 HOMEWORK 2 DAVID MEREDITH MATH 890 HOMEWORK 2 DAVID MEREDITH (1) Suppose P and Q are polyhedra. Then P Q is a polyhedron. Moreover if P and Q are polytopes then P Q is a polytope. The facets of P Q are either F Q where F is a facet

More information

Math 5593 Linear Programming Lecture Notes

Math 5593 Linear Programming Lecture Notes Math 5593 Linear Programming Lecture Notes Unit II: Theory & Foundations (Convex Analysis) University of Colorado Denver, Fall 2013 Topics 1 Convex Sets 1 1.1 Basic Properties (Luenberger-Ye Appendix B.1).........................

More information

Polytopes Course Notes

Polytopes Course Notes Polytopes Course Notes Carl W. Lee Department of Mathematics University of Kentucky Lexington, KY 40506 lee@ms.uky.edu Fall 2013 i Contents 1 Polytopes 1 1.1 Convex Combinations and V-Polytopes.....................

More information

Conic Duality. yyye

Conic Duality.  yyye Conic Linear Optimization and Appl. MS&E314 Lecture Note #02 1 Conic Duality Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/

More information

Convexity: an introduction

Convexity: an introduction Convexity: an introduction Geir Dahl CMA, Dept. of Mathematics and Dept. of Informatics University of Oslo 1 / 74 1. Introduction 1. Introduction what is convexity where does it arise main concepts and

More information

The Matrix-Tree Theorem and Its Applications to Complete and Complete Bipartite Graphs

The Matrix-Tree Theorem and Its Applications to Complete and Complete Bipartite Graphs The Matrix-Tree Theorem and Its Applications to Complete and Complete Bipartite Graphs Frankie Smith Nebraska Wesleyan University fsmith@nebrwesleyan.edu May 11, 2015 Abstract We will look at how to represent

More information

Lecture 2 - Introduction to Polytopes

Lecture 2 - Introduction to Polytopes Lecture 2 - Introduction to Polytopes Optimization and Approximation - ENS M1 Nicolas Bousquet 1 Reminder of Linear Algebra definitions Let x 1,..., x m be points in R n and λ 1,..., λ m be real numbers.

More information

Pebble Sets in Convex Polygons

Pebble Sets in Convex Polygons 2 1 Pebble Sets in Convex Polygons Kevin Iga, Randall Maddox June 15, 2005 Abstract Lukács and András posed the problem of showing the existence of a set of n 2 points in the interior of a convex n-gon

More information

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

On the null space of a Colin de Verdière matrix On the null space of a Colin de Verdière matrix László Lovász 1 and Alexander Schrijver 2 Dedicated to the memory of François Jaeger Abstract. Let G = (V, E) be a 3-connected planar graph, with V = {1,...,

More information

arxiv: v1 [math.co] 12 Dec 2017

arxiv: v1 [math.co] 12 Dec 2017 arxiv:1712.04381v1 [math.co] 12 Dec 2017 Semi-reflexive polytopes Tiago Royer Abstract The Ehrhart function L P(t) of a polytope P is usually defined only for integer dilation arguments t. By allowing

More information

We have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance.

We have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance. Solid geometry We have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance. First, note that everything we have proven for the

More information

Algebraic Graph Theory- Adjacency Matrix and Spectrum

Algebraic Graph Theory- Adjacency Matrix and Spectrum Algebraic Graph Theory- Adjacency Matrix and Spectrum Michael Levet December 24, 2013 Introduction This tutorial will introduce the adjacency matrix, as well as spectral graph theory. For those familiar

More information

Three applications of Euler s formula. Chapter 10

Three applications of Euler s formula. Chapter 10 Three applications of Euler s formula Chapter 10 A graph is planar if it can be drawn in the plane R without crossing edges (or, equivalently, on the -dimensional sphere S ). We talk of a plane graph if

More information

2017 SOLUTIONS (PRELIMINARY VERSION)

2017 SOLUTIONS (PRELIMINARY VERSION) SIMON MARAIS MATHEMATICS COMPETITION 07 SOLUTIONS (PRELIMINARY VERSION) This document will be updated to include alternative solutions provided by contestants, after the competition has been mared. Problem

More information

FACES OF CONVEX SETS

FACES OF CONVEX SETS FACES OF CONVEX SETS VERA ROSHCHINA Abstract. We remind the basic definitions of faces of convex sets and their basic properties. For more details see the classic references [1, 2] and [4] for polytopes.

More information

Mathematical and Algorithmic Foundations Linear Programming and Matchings

Mathematical and Algorithmic Foundations Linear Programming and Matchings Adavnced Algorithms Lectures Mathematical and Algorithmic Foundations Linear Programming and Matchings Paul G. Spirakis Department of Computer Science University of Patras and Liverpool Paul G. Spirakis

More information

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

Lecture 3. Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets. Tepper School of Business Carnegie Mellon University, Pittsburgh Lecture 3 Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets Gérard Cornuéjols Tepper School of Business Carnegie Mellon University, Pittsburgh January 2016 Mixed Integer Linear Programming

More information

LATIN SQUARES AND THEIR APPLICATION TO THE FEASIBLE SET FOR ASSIGNMENT PROBLEMS

LATIN SQUARES AND THEIR APPLICATION TO THE FEASIBLE SET FOR ASSIGNMENT PROBLEMS LATIN SQUARES AND THEIR APPLICATION TO THE FEASIBLE SET FOR ASSIGNMENT PROBLEMS TIMOTHY L. VIS Abstract. A significant problem in finite optimization is the assignment problem. In essence, the assignment

More information

Automorphism Groups of Cyclic Polytopes

Automorphism Groups of Cyclic Polytopes 8 Automorphism Groups of Cyclic Polytopes (Volker Kaibel and Arnold Waßmer ) It is probably well-known to most polytope theorists that the combinatorial automorphism group of a cyclic d-polytope with n

More information

Definition. Given a (v,k,λ)- BIBD, (X,B), a set of disjoint blocks of B which partition X is called a parallel class.

Definition. Given a (v,k,λ)- BIBD, (X,B), a set of disjoint blocks of B which partition X is called a parallel class. Resolvable BIBDs Definition Given a (v,k,λ)- BIBD, (X,B), a set of disjoint blocks of B which partition X is called a parallel class. A partition of B into parallel classes (there must be r of them) is

More information

1. CONVEX POLYGONS. Definition. A shape D in the plane is convex if every line drawn between two points in D is entirely inside D.

1. CONVEX POLYGONS. Definition. A shape D in the plane is convex if every line drawn between two points in D is entirely inside D. 1. CONVEX POLYGONS Definition. A shape D in the plane is convex if every line drawn between two points in D is entirely inside D. Convex 6 gon Another convex 6 gon Not convex Question. Why is the third

More information

On the Minimum Number of Convex Quadrilaterals in Point Sets of Given Numbers of Points

On the Minimum Number of Convex Quadrilaterals in Point Sets of Given Numbers of Points On the Minimum Number of Convex Quadrilaterals in Point Sets of Given Numbers of Points Hu Yuzhong Chen Luping Zhu Hui Ling Xiaofeng (Supervisor) Abstract Consider the following problem. Given n, k N,

More information

Interval-Vector Polytopes

Interval-Vector Polytopes Interval-Vector Polytopes Jessica De Silva Gabriel Dorfsman-Hopkins California State University, Stanislaus Dartmouth College Joseph Pruitt California State University, Long Beach July 28, 2012 Abstract

More information

c 2004 Society for Industrial and Applied Mathematics

c 2004 Society for Industrial and Applied Mathematics SIAM J. MATRIX ANAL. APPL. Vol. 26, No. 2, pp. 390 399 c 2004 Society for Industrial and Applied Mathematics HERMITIAN MATRICES, EIGENVALUE MULTIPLICITIES, AND EIGENVECTOR COMPONENTS CHARLES R. JOHNSON

More information

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

POLYHEDRAL GEOMETRY. Convex functions and sets. Mathematical Programming Niels Lauritzen Recall that a subset C R n is convex if POLYHEDRAL GEOMETRY Mathematical Programming Niels Lauritzen 7.9.2007 Convex functions and sets Recall that a subset C R n is convex if {λx + (1 λ)y 0 λ 1} C for every x, y C and 0 λ 1. A function f :

More information

On the Relationships between Zero Forcing Numbers and Certain Graph Coverings

On the Relationships between Zero Forcing Numbers and Certain Graph Coverings On the Relationships between Zero Forcing Numbers and Certain Graph Coverings Fatemeh Alinaghipour Taklimi, Shaun Fallat 1,, Karen Meagher 2 Department of Mathematics and Statistics, University of Regina,

More information

Monotone Paths in Geometric Triangulations

Monotone Paths in Geometric Triangulations Monotone Paths in Geometric Triangulations Adrian Dumitrescu Ritankar Mandal Csaba D. Tóth November 19, 2017 Abstract (I) We prove that the (maximum) number of monotone paths in a geometric triangulation

More information

Rubber bands. Chapter Rubber band representation

Rubber bands. Chapter Rubber band representation Chapter 1 Rubber bands In the previous chapter, we already used the idea of looking at the graph geometrically, by placing its nodes on the line and replacing the edges by rubber bands. Since, however,

More information

Simplicial Cells in Arrangements of Hyperplanes

Simplicial Cells in Arrangements of Hyperplanes Simplicial Cells in Arrangements of Hyperplanes Christoph Dätwyler 05.01.2013 This paper is a report written due to the authors presentation of a paper written by Shannon [1] in 1977. The presentation

More information

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

Discrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity Discrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity Marc Uetz University of Twente m.uetz@utwente.nl Lecture 5: sheet 1 / 26 Marc Uetz Discrete Optimization Outline 1 Min-Cost Flows

More information

EXTREME POINTS AND AFFINE EQUIVALENCE

EXTREME POINTS AND AFFINE EQUIVALENCE EXTREME POINTS AND AFFINE EQUIVALENCE The purpose of this note is to use the notions of extreme points and affine transformations which are studied in the file affine-convex.pdf to prove that certain standard

More information

TROPICAL CONVEXITY MIKE DEVELIN AND BERND STURMFELS

TROPICAL CONVEXITY MIKE DEVELIN AND BERND STURMFELS TROPICAL CONVEXITY MIKE DEVELIN AND BERND STURMFELS Abstract. The notions of convexity and convex polytopes are introduced in the setting of tropical geometry. Combinatorial types of tropical polytopes

More information

arxiv: v1 [math.co] 25 Sep 2015

arxiv: v1 [math.co] 25 Sep 2015 A BASIS FOR SLICING BIRKHOFF POLYTOPES TREVOR GLYNN arxiv:1509.07597v1 [math.co] 25 Sep 2015 Abstract. We present a change of basis that may allow more efficient calculation of the volumes of Birkhoff

More information

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

DM545 Linear and Integer Programming. Lecture 2. The Simplex Method. Marco Chiarandini DM545 Linear and Integer Programming Lecture 2 The Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. 3. 4. Standard Form Basic Feasible Solutions

More information

Convex Optimization. 2. Convex Sets. Prof. Ying Cui. Department of Electrical Engineering Shanghai Jiao Tong University. SJTU Ying Cui 1 / 33

Convex Optimization. 2. Convex Sets. Prof. Ying Cui. Department of Electrical Engineering Shanghai Jiao Tong University. SJTU Ying Cui 1 / 33 Convex Optimization 2. Convex Sets Prof. Ying Cui Department of Electrical Engineering Shanghai Jiao Tong University 2018 SJTU Ying Cui 1 / 33 Outline Affine and convex sets Some important examples Operations

More information

Chapter 3. Quadric hypersurfaces. 3.1 Quadric hypersurfaces Denition.

Chapter 3. Quadric hypersurfaces. 3.1 Quadric hypersurfaces Denition. Chapter 3 Quadric hypersurfaces 3.1 Quadric hypersurfaces. 3.1.1 Denition. Denition 1. In an n-dimensional ane space A; given an ane frame fo;! e i g: A quadric hypersurface in A is a set S consisting

More information

Nesting points in the sphere. Dan Archdeacon. University of Vermont. Feliu Sagols.

Nesting points in the sphere. Dan Archdeacon. University of Vermont.   Feliu Sagols. Nesting points in the sphere Dan Archdeacon Dept. of Computer Science University of Vermont Burlington, VT, USA 05405 e-mail: dan.archdeacon@uvm.edu Feliu Sagols Dept. of Computer Science University of

More information

Lecture 9 - Matrix Multiplication Equivalences and Spectral Graph Theory 1

Lecture 9 - Matrix Multiplication Equivalences and Spectral Graph Theory 1 CME 305: Discrete Mathematics and Algorithms Instructor: Professor Aaron Sidford (sidford@stanfordedu) February 6, 2018 Lecture 9 - Matrix Multiplication Equivalences and Spectral Graph Theory 1 In the

More information

Geometry. Every Simplicial Polytope with at Most d + 4 Vertices Is a Quotient of a Neighborly Polytope. U. H. Kortenkamp. 1.

Geometry. Every Simplicial Polytope with at Most d + 4 Vertices Is a Quotient of a Neighborly Polytope. U. H. Kortenkamp. 1. Discrete Comput Geom 18:455 462 (1997) Discrete & Computational Geometry 1997 Springer-Verlag New York Inc. Every Simplicial Polytope with at Most d + 4 Vertices Is a Quotient of a Neighborly Polytope

More information

The orientability of small covers and coloring simple polytopes. Nishimura, Yasuzo; Nakayama, Hisashi. Osaka Journal of Mathematics. 42(1) P.243-P.

The orientability of small covers and coloring simple polytopes. Nishimura, Yasuzo; Nakayama, Hisashi. Osaka Journal of Mathematics. 42(1) P.243-P. Title Author(s) The orientability of small covers and coloring simple polytopes Nishimura, Yasuzo; Nakayama, Hisashi Citation Osaka Journal of Mathematics. 42(1) P.243-P.256 Issue Date 2005-03 Text Version

More information

arxiv: v1 [math.co] 27 Feb 2015

arxiv: v1 [math.co] 27 Feb 2015 Mode Poset Probability Polytopes Guido Montúfar 1 and Johannes Rauh 2 arxiv:1503.00572v1 [math.co] 27 Feb 2015 1 Max Planck Institute for Mathematics in the Sciences, Inselstraße 22, 04103 Leipzig, Germany,

More information

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

60 2 Convex sets. {x a T x b} {x ã T x b} 60 2 Convex sets Exercises Definition of convexity 21 Let C R n be a convex set, with x 1,, x k C, and let θ 1,, θ k R satisfy θ i 0, θ 1 + + θ k = 1 Show that θ 1x 1 + + θ k x k C (The definition of convexity

More information

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

MA4254: Discrete Optimization. Defeng Sun. Department of Mathematics National University of Singapore Office: S Telephone: MA4254: Discrete Optimization Defeng Sun Department of Mathematics National University of Singapore Office: S14-04-25 Telephone: 6516 3343 Aims/Objectives: Discrete optimization deals with problems of

More information

6. Lecture notes on matroid intersection

6. Lecture notes on matroid intersection Massachusetts Institute of Technology 18.453: Combinatorial Optimization Michel X. Goemans May 2, 2017 6. Lecture notes on matroid intersection One nice feature about matroids is that a simple greedy algorithm

More information

All 0-1 Polytopes are. Abstract. We study the facial structure of two important permutation polytopes

All 0-1 Polytopes are. Abstract. We study the facial structure of two important permutation polytopes All 0-1 Polytopes are Traveling Salesman Polytopes L.J. Billera and A. Sarangarajan y Abstract We study the facial structure of two important permutation polytopes in R n2, the Birkho or assignment polytope

More information

INTRODUCTION TO THE HOMOLOGY GROUPS OF COMPLEXES

INTRODUCTION TO THE HOMOLOGY GROUPS OF COMPLEXES INTRODUCTION TO THE HOMOLOGY GROUPS OF COMPLEXES RACHEL CARANDANG Abstract. This paper provides an overview of the homology groups of a 2- dimensional complex. It then demonstrates a proof of the Invariance

More information

ON THE STRONGLY REGULAR GRAPH OF PARAMETERS

ON THE STRONGLY REGULAR GRAPH OF PARAMETERS ON THE STRONGLY REGULAR GRAPH OF PARAMETERS (99, 14, 1, 2) SUZY LOU AND MAX MURIN Abstract. In an attempt to find a strongly regular graph of parameters (99, 14, 1, 2) or to disprove its existence, we

More information

arxiv: v1 [math.co] 15 Dec 2009

arxiv: v1 [math.co] 15 Dec 2009 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

More information

Lecture 6: Faces, Facets

Lecture 6: Faces, Facets IE 511: Integer Programming, Spring 2019 31 Jan, 2019 Lecturer: Karthik Chandrasekaran Lecture 6: Faces, Facets Scribe: Setareh Taki Disclaimer: These notes have not been subjected to the usual scrutiny

More information

CHAPTER 2. Graphs. 1. Introduction to Graphs and Graph Isomorphism

CHAPTER 2. Graphs. 1. Introduction to Graphs and Graph Isomorphism CHAPTER 2 Graphs 1. Introduction to Graphs and Graph Isomorphism 1.1. The Graph Menagerie. Definition 1.1.1. A simple graph G = (V, E) consists of a set V of vertices and a set E of edges, represented

More information

THE INSULATION SEQUENCE OF A GRAPH

THE INSULATION SEQUENCE OF A GRAPH THE INSULATION SEQUENCE OF A GRAPH ELENA GRIGORESCU Abstract. In a graph G, a k-insulated set S is a subset of the vertices of G such that every vertex in S is adjacent to at most k vertices in S, and

More information

MAT 3271: Selected Solutions to the Assignment 6

MAT 3271: Selected Solutions to the Assignment 6 Chapter 2: Major Exercises MAT 3271: Selected Solutions to the Assignment 6 1. Since a projective plan is a model of incidence geometry, Incidence Axioms 1-3 and Propositions 2.1-2.5 (which follow logically

More information

A Nim game played on graphs II

A Nim game played on graphs II Theoretical Computer Science 304 (2003) 401 419 www.elsevier.com/locate/tcs A Nim game played on graphs II Masahiko Fukuyama Graduate School of Mathematical Sciences, University of Tokyo, 3-8-1 Komaba,

More information

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

3. The Simplex algorithmn The Simplex algorithmn 3.1 Forms of linear programs 11 3.1 Forms of linear programs... 12 3.2 Basic feasible solutions... 13 3.3 The geometry of linear programs... 14 3.4 Local search among basic feasible solutions... 15 3.5 Organization in tableaus...

More information

6 Randomized rounding of semidefinite programs

6 Randomized rounding of semidefinite programs 6 Randomized rounding of semidefinite programs We now turn to a new tool which gives substantially improved performance guarantees for some problems We now show how nonlinear programming relaxations can

More information

Ma/CS 6b Class 13: Counting Spanning Trees

Ma/CS 6b Class 13: Counting Spanning Trees Ma/CS 6b Class 13: Counting Spanning Trees By Adam Sheffer Reminder: Spanning Trees A spanning tree is a tree that contains all of the vertices of the graph. A graph can contain many distinct spanning

More information

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

Ma/CS 6b Class 26: Art Galleries and Politicians Ma/CS 6b Class 26: Art Galleries and Politicians By Adam Sheffer The Art Gallery Problem Problem. We wish to place security cameras at a gallery, such that they cover it completely. Every camera can cover

More information

An Investigation of Closed Geodesics on Regular Polyhedra

An Investigation of Closed Geodesics on Regular Polyhedra An Investigation of Closed Geodesics on Regular Polyhedra Tony Scoles Southern Illinois University Edwardsville May 13, 2008 1 Introduction This paper was undertaken to examine, in detail, results from

More information

Combinatorial Geometry & Topology arising in Game Theory and Optimization

Combinatorial Geometry & Topology arising in Game Theory and Optimization Combinatorial Geometry & Topology arising in Game Theory and Optimization Jesús A. De Loera University of California, Davis LAST EPISODE... We discuss the content of the course... Convex Sets A set is

More information

Lecture 2 September 3

Lecture 2 September 3 EE 381V: Large Scale Optimization Fall 2012 Lecture 2 September 3 Lecturer: Caramanis & Sanghavi Scribe: Hongbo Si, Qiaoyang Ye 2.1 Overview of the last Lecture The focus of the last lecture was to give

More information

Embeddability of Arrangements of Pseudocircles into the Sphere

Embeddability of Arrangements of Pseudocircles into the Sphere Embeddability of Arrangements of Pseudocircles into the Sphere Ronald Ortner Department Mathematik und Informationstechnologie, Montanuniversität Leoben, Franz-Josef-Straße 18, 8700-Leoben, Austria Abstract

More information

Geometry. Tropical Secant Varieties of Linear Spaces. Mike Develin. 1. Introduction

Geometry. Tropical Secant Varieties of Linear Spaces. Mike Develin. 1. Introduction Discrete Comput Geom 35:117 129 (2006) DOI: 101007/s00454-005-1182-2 Discrete & Computational Geometry 2005 Springer Science+Business Media, Inc Tropical Secant Varieties of Linear Spaces Mike Develin

More information

AMS : Combinatorial Optimization Homework Problems - Week V

AMS : Combinatorial Optimization Homework Problems - Week V AMS 553.766: Combinatorial Optimization Homework Problems - Week V For the following problems, A R m n will be m n matrices, and b R m. An affine subspace is the set of solutions to a a system of linear

More information

On the Number of Tilings of a Square by Rectangles

On the Number of Tilings of a Square by Rectangles University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange University of Tennessee Honors Thesis Projects University of Tennessee Honors Program 5-2012 On the Number of Tilings

More information

POLYTOPES. Grünbaum and Shephard [40] remarked that there were three developments which foreshadowed the modern theory of convex polytopes.

POLYTOPES. Grünbaum and Shephard [40] remarked that there were three developments which foreshadowed the modern theory of convex polytopes. POLYTOPES MARGARET A. READDY 1. Lecture I: Introduction to Polytopes and Face Enumeration Grünbaum and Shephard [40] remarked that there were three developments which foreshadowed the modern theory of

More information

7. The Gauss-Bonnet theorem

7. The Gauss-Bonnet theorem 7. The Gauss-Bonnet theorem 7.1 Hyperbolic polygons In Euclidean geometry, an n-sided polygon is a subset of the Euclidean plane bounded by n straight lines. Thus the edges of a Euclidean polygon are formed

More information

Technische Universität München Zentrum Mathematik

Technische Universität München Zentrum Mathematik Technische Universität München Zentrum Mathematik Prof. Dr. Dr. Jürgen Richter-Gebert, Bernhard Werner Projective Geometry SS 208 https://www-m0.ma.tum.de/bin/view/lehre/ss8/pgss8/webhome Solutions for

More information

Winning Positions in Simplicial Nim

Winning Positions in Simplicial Nim Winning Positions in Simplicial Nim David Horrocks Department of Mathematics and Statistics University of Prince Edward Island Charlottetown, Prince Edward Island, Canada, C1A 4P3 dhorrocks@upei.ca Submitted:

More information

Lecture 5: Duality Theory

Lecture 5: Duality Theory Lecture 5: Duality Theory Rajat Mittal IIT Kanpur The objective of this lecture note will be to learn duality theory of linear programming. We are planning to answer following questions. What are hyperplane

More information

Mathematical Programming and Research Methods (Part II)

Mathematical Programming and Research Methods (Part II) Mathematical Programming and Research Methods (Part II) 4. Convexity and Optimization Massimiliano Pontil (based on previous lecture by Andreas Argyriou) 1 Today s Plan Convex sets and functions Types

More information

LOWER BOUNDS FOR SIMPLICIAL COVERS AND TRIANGULATIONS OF CUBES

LOWER BOUNDS FOR SIMPLICIAL COVERS AND TRIANGULATIONS OF CUBES LOWER BOUNDS FOR SIMPLICIAL COVERS AND TRIANGULATIONS OF CUBES ADAM BLISS AND FRANCIS EDWARD SU Abstract. We show that the size of a minimal simplicial cover of a polytope P is a lower bound for the size

More information

Preferred directions for resolving the non-uniqueness of Delaunay triangulations

Preferred directions for resolving the non-uniqueness of Delaunay triangulations Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation

More information

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

Advanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs Advanced Operations Research Techniques IE316 Quiz 1 Review Dr. Ted Ralphs IE316 Quiz 1 Review 1 Reading for The Quiz Material covered in detail in lecture. 1.1, 1.4, 2.1-2.6, 3.1-3.3, 3.5 Background material

More information

4. Simplicial Complexes and Simplicial Homology

4. Simplicial Complexes and Simplicial Homology MATH41071/MATH61071 Algebraic topology Autumn Semester 2017 2018 4. Simplicial Complexes and Simplicial Homology Geometric simplicial complexes 4.1 Definition. A finite subset { v 0, v 1,..., v r } R n

More information

Equipartite polytopes and graphs

Equipartite polytopes and graphs Equipartite polytopes and graphs Branko Grünbaum Tomáš Kaiser Daniel Král Moshe Rosenfeld Abstract A graph G of even order is weakly equipartite if for any partition of its vertex set into subsets V 1

More information

Localization in Graphs. Richardson, TX Azriel Rosenfeld. Center for Automation Research. College Park, MD

Localization in Graphs. Richardson, TX Azriel Rosenfeld. Center for Automation Research. College Park, MD CAR-TR-728 CS-TR-3326 UMIACS-TR-94-92 Samir Khuller Department of Computer Science Institute for Advanced Computer Studies University of Maryland College Park, MD 20742-3255 Localization in Graphs Azriel

More information

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

EC 521 MATHEMATICAL METHODS FOR ECONOMICS. Lecture 2: Convex Sets EC 51 MATHEMATICAL METHODS FOR ECONOMICS Lecture : Convex Sets Murat YILMAZ Boğaziçi University In this section, we focus on convex sets, separating hyperplane theorems and Farkas Lemma. And as an application

More information

Simplex Algorithm in 1 Slide

Simplex Algorithm in 1 Slide Administrivia 1 Canonical form: Simplex Algorithm in 1 Slide If we do pivot in A r,s >0, where c s

More information

arxiv: v1 [math.mg] 13 Mar 2018

arxiv: v1 [math.mg] 13 Mar 2018 Characterizing face and flag vector pairs for polytopes Hannah Sjöberg Institut für Mathematik FU Berlin, Arnimallee 14195 Berlin, Germany sjoberg@math.fu-berlin.de Günter M. Ziegler Institut für Mathematik

More information

Exercise set 2 Solutions

Exercise set 2 Solutions Exercise set 2 Solutions Let H and H be the two components of T e and let F E(T ) consist of the edges of T with one endpoint in V (H), the other in V (H ) Since T is connected, F Furthermore, since T

More information

Chapter 8. Voronoi Diagrams. 8.1 Post Oce Problem

Chapter 8. Voronoi Diagrams. 8.1 Post Oce Problem Chapter 8 Voronoi Diagrams 8.1 Post Oce Problem Suppose there are n post oces p 1,... p n in a city. Someone who is located at a position q within the city would like to know which post oce is closest

More information

GENERALIZED CPR-GRAPHS AND APPLICATIONS

GENERALIZED CPR-GRAPHS AND APPLICATIONS Volume 5, Number 2, Pages 76 105 ISSN 1715-0868 GENERALIZED CPR-GRAPHS AND APPLICATIONS DANIEL PELLICER AND ASIA IVIĆ WEISS Abstract. We give conditions for oriented labeled graphs that must be satisfied

More information

Alternating Projections

Alternating Projections Alternating Projections Stephen Boyd and Jon Dattorro EE392o, Stanford University Autumn, 2003 1 Alternating projection algorithm Alternating projections is a very simple algorithm for computing a point

More information

2 Geometry Solutions

2 Geometry Solutions 2 Geometry Solutions jacques@ucsd.edu Here is give problems and solutions in increasing order of difficulty. 2.1 Easier problems Problem 1. What is the minimum number of hyperplanar slices to make a d-dimensional

More information

Equipartite polytopes and graphs

Equipartite polytopes and graphs Equipartite polytopes and graphs Branko Grünbaum Tomáš Kaiser Daniel Král Moshe Rosenfeld Abstract A graph G of even order is weakly equipartite if for any partition of its vertex set into subsets V 1

More information

N-flips in 3-dimensional Small Cove Transformation Groups and Combinato. 数理解析研究所講究録別冊 = RIMS Kokyuroku Bessa (2013), B39:

N-flips in 3-dimensional Small Cove Transformation Groups and Combinato. 数理解析研究所講究録別冊 = RIMS Kokyuroku Bessa (2013), B39: Title N-flips in 3-dimensional Small Cove Transformation Groups and Combinato Author(s) NISHIMURA, Yasuzo Citation 数理解析研究所講究録別冊 = RIMS Kokyuroku Bessa (2013), B39: 137-148 Issue Date 2013-04 URL http://hdl.handle.net/2433/207819

More information

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

Division of the Humanities and Social Sciences. Convex Analysis and Economic Theory Winter Separation theorems Division of the Humanities and Social Sciences Ec 181 KC Border Convex Analysis and Economic Theory Winter 2018 Topic 8: Separation theorems 8.1 Hyperplanes and half spaces Recall that a hyperplane in

More information

Line Graphs and Circulants

Line Graphs and Circulants Line Graphs and Circulants Jason Brown and Richard Hoshino Department of Mathematics and Statistics Dalhousie University Halifax, Nova Scotia, Canada B3H 3J5 Abstract The line graph of G, denoted L(G),

More information

Strong edge coloring of subcubic graphs

Strong edge coloring of subcubic graphs Strong edge coloring of subcubic graphs Hervé Hocquard a, Petru Valicov a a LaBRI (Université Bordeaux 1), 351 cours de la Libération, 33405 Talence Cedex, France Abstract A strong edge colouring of a

More information

CMSC Honors Discrete Mathematics

CMSC Honors Discrete Mathematics CMSC 27130 Honors Discrete Mathematics Lectures by Alexander Razborov Notes by Justin Lubin The University of Chicago, Autumn 2017 1 Contents I Number Theory 4 1 The Euclidean Algorithm 4 2 Mathematical

More information

Genus of the cartesian product of triangles

Genus of the cartesian product of triangles Genus of the cartesian product of triangles Michal Kotrbčík Faculty of Informatics Masaryk University Brno, Czech Republic kotrbcik@fi.muni.cz Tomaž Pisanski FAMNIT University of Primorska Koper, Slovenia

More information

maximize c, x subject to Ax b,

maximize c, x subject to Ax b, Lecture 8 Linear programming is about problems of the form maximize c, x subject to Ax b, where A R m n, x R n, c R n, and b R m, and the inequality sign means inequality in each row. The feasible set

More information

MATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix.

MATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix. MATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix. Row echelon form A matrix is said to be in the row echelon form if the leading entries shift to the

More information

Module 7. Independent sets, coverings. and matchings. Contents

Module 7. Independent sets, coverings. and matchings. Contents Module 7 Independent sets, coverings Contents and matchings 7.1 Introduction.......................... 152 7.2 Independent sets and coverings: basic equations..... 152 7.3 Matchings in bipartite graphs................

More information

CAT(0)-spaces. Münster, June 22, 2004

CAT(0)-spaces. Münster, June 22, 2004 CAT(0)-spaces Münster, June 22, 2004 CAT(0)-space is a term invented by Gromov. Also, called Hadamard space. Roughly, a space which is nonpositively curved and simply connected. C = Comparison or Cartan

More information

Some Advanced Topics in Linear Programming

Some Advanced Topics in Linear Programming Some Advanced Topics in Linear Programming Matthew J. Saltzman July 2, 995 Connections with Algebra and Geometry In this section, we will explore how some of the ideas in linear programming, duality theory,

More information

The important function we will work with is the omega map ω, which we now describe.

The important function we will work with is the omega map ω, which we now describe. 20 MARGARET A. READDY 3. Lecture III: Hyperplane arrangements & zonotopes; Inequalities: a first look 3.1. Zonotopes. Recall that a zonotope Z can be described as the Minkowski sum of line segments: Z

More information