Basic Polyhedral theory

Similar documents
cycle that does not cross any edges (including its own), then it has at least

CPSC 665 : An Algorithmist s Toolkit Lecture 4 : 21 Jan Linear Programming

Derangements and Applications

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Combinatorial Networks Week 1, March 11-12

Finding low cost TSP and 2-matching solutions using certain half integer subtour vertices

The Matrix Exponential

Searching Linked Lists. Perfect Skip List. Building a Skip List. Skip List Analysis (1) Assume the list is sorted, but is stored in a linked list.

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

Week 3: Connected Subgraphs

The Matrix Exponential

NEW APPLICATIONS OF THE ABEL-LIOUVILLE FORMULA

On spanning trees and cycles of multicolored point sets with few intersections

Hardy-Littlewood Conjecture and Exceptional real Zero. JinHua Fei. ChangLing Company of Electronic Technology Baoji Shannxi P.R.

Supplementary Materials

Construction of asymmetric orthogonal arrays of strength three via a replacement method

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

1 Minimum Cut Problem

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013

Propositional Logic. Combinatorial Problem Solving (CPS) Albert Oliveras Enric Rodríguez-Carbonell. May 17, 2018

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

10. The Discrete-Time Fourier Transform (DTFT)

The Equitable Dominating Graph

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

y = 2xe x + x 2 e x at (0, 3). solution: Since y is implicitly related to x we have to use implicit differentiation: 3 6y = 0 y = 1 2 x ln(b) ln(b)

Section 6.1. Question: 2. Let H be a subgroup of a group G. Then H operates on G by left multiplication. Describe the orbits for this operation.

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

Where k is either given or determined from the data and c is an arbitrary constant.

2.3 Matrix Formulation

MATH 319, WEEK 15: The Fundamental Matrix, Non-Homogeneous Systems of Differential Equations

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

Network Congestion Games

EXST Regression Techniques Page 1

SCHUR S THEOREM REU SUMMER 2005

Computing and Communications -- Network Coding

Strongly Connected Components

(1) Then we could wave our hands over this and it would become:

Deift/Zhou Steepest descent, Part I

Thus, because if either [G : H] or [H : K] is infinite, then [G : K] is infinite, then [G : K] = [G : H][H : K] for all infinite cases.

u 3 = u 3 (x 1, x 2, x 3 )

Section 11.6: Directional Derivatives and the Gradient Vector

Differential Equations

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

UNTYPED LAMBDA CALCULUS (II)

Self-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016

10. EXTENDING TRACTABILITY

Search sequence databases 3 10/25/2016

Addition of angular momentum

Hydrogen Atom and One Electron Ions

perm4 A cnt 0 for for if A i 1 A i cnt cnt 1 cnt i j. j k. k l. i k. j l. i l

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

GEOMETRICAL PHENOMENA IN THE PHYSICS OF SUBATOMIC PARTICLES. Eduard N. Klenov* Rostov-on-Don, Russia

Homework #3. 1 x. dx. It therefore follows that a sum of the

MA 262, Spring 2018, Final exam Version 01 (Green)

SOME PARAMETERS ON EQUITABLE COLORING OF PRISM AND CIRCULANT GRAPH.

1 Isoparametric Concept

EEO 401 Digital Signal Processing Prof. Mark Fowler

Recall that by Theorems 10.3 and 10.4 together provide us the estimate o(n2 ), S(q) q 9, q=1

DIFFERENTIAL EQUATION

Addition of angular momentum

AS 5850 Finite Element Analysis

COMPUTER GENERATED HOLOGRAMS Optical Sciences 627 W.J. Dallas (Monday, April 04, 2005, 8:35 AM) PART I: CHAPTER TWO COMB MATH.

Higher order derivatives

Chapter 10. The singular integral Introducing S(n) and J(n)

Abstract Interpretation: concrete and abstract semantics

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Chapter 6 Folding. Folding

Finite element discretization of Laplace and Poisson equations

INTEGRATION BY PARTS


Lie Groups HW7. Wang Shuai. November 2015

Coupled Pendulums. Two normal modes.

First order differential equation Linear equation; Method of integrating factors

10. Limits involving infinity

PROOF OF FIRST STANDARD FORM OF NONELEMENTARY FUNCTIONS

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

VII. Quantum Entanglement

Mathematics. Complex Number rectangular form. Quadratic equation. Quadratic equation. Complex number Functions: sinusoids. Differentiation Integration

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

Approximation and Inapproximation for The Influence Maximization Problem in Social Networks under Deterministic Linear Threshold Model

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

nd the particular orthogonal trajectory from the family of orthogonal trajectories passing through point (0; 1).

EECE 301 Signals & Systems Prof. Mark Fowler

MCE503: Modeling and Simulation of Mechatronic Systems Discussion on Bond Graph Sign Conventions for Electrical Systems

An Application of Hardy-Littlewood Conjecture. JinHua Fei. ChangLing Company of Electronic Technology Baoji Shannxi P.R.China

LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS DONE RIGHT KYLE ORMSBY

3 Finite Element Parametric Geometry

CS 361 Meeting 12 10/3/18

Limiting value of higher Mahler measure

Sundials and Linear Algebra

Quasi-Classical States of the Simple Harmonic Oscillator

Stochastic Submodular Maximization

Slide 1. Slide 2. Slide 3 DIGITAL SIGNAL PROCESSING CLASSIFICATION OF SIGNALS

Examples and applications on SSSP and MST

First derivative analysis

Roadmap. XML Indexing. DataGuide example. DataGuides. Strong DataGuides. Multiple DataGuides for same data. CPS Topics in Database Systems

Calculus concepts derivatives

On the irreducibility of some polynomials in two variables

DISTRIBUTION OF DIFFERENCE BETWEEN INVERSES OF CONSECUTIVE INTEGERS MODULO P

A Propagating Wave Packet Group Velocity Dispersion

Transcription:

Basic Polyhdral thory Th st P = { A b} is calld a polyhdron. Lmma 1. Eithr th systm A = b, b 0, 0 has a solution or thr is a vctorπ such that π A 0, πb < 0 Thr cass, if solution in top row dos not ist thn thr ists a π such that th bottom is tru. (in ach column, continuous). A = b A b A b π A 0, πb < 0 π A 0, πb < 0 π A 0, πb > 0 If { A b, intgr } = { C d, intgr} and { A b } { C d }, w say that th first formulation is Tightr (& bttr). Howvr w r intrstd in intgr solutions. Intgral polyhdra Givn a polyhdron P, its intgr hull P I is th conv hull of all intgral vctors in P. For any rational polyhdron P, th P I is also a polyhdron. Th following mthod gnrats valid inqualitis for P I (Chvatal s thorm) Lt P= { A b, 0}, whr A and b intgral, thn µ A µb is valid for P, for any µ 0, and µ A µ b is valid for P, bcaus 0, thn µ A µ b (1) is valid for P I bcaus is intgral. Lt us call S 1 th st of inqualitis using (1), and st S t+1 th st obtaind by applying this opration to S t 0 1. Ths sts gnrat a squnc P = P P... PI Eampl: X(1)+(4) +(6) =1 X(1)+(5)+(7) =1 X(3)+(4)+(8) =1 X(3)+(5)+(9) =1 X(2)+(7)+(9) =1 X(2)+(6)+(8) =1 Considr taking ½*(first four inqualitis) to gt (1/2+1/2)X(1)+ (1/2+1/2)(3)+ (1/2+1/2)(4)+ (1/2+1/2)(5) + 1/2 ((6)+(7)+(8)+(9))=2 1 1 X(1)+ 1 (3)+ 1 (4)+ 1 (5) + 2 ((6)+(7)+(8)+(9))= 2, 1 1 Whr 2 =0 so 2 ((6)+(7)+(8)+(9))=0 X(1)+(3)+(4)+(5)=2 Considr taking last thr inqualitis and first on with multiplir of ½

X(2)+(6)+(9)=2 Total unimodularity: If A is totally unimodular (TU) thn n P(b)= { R + m A b} is intgral b Z for which it is not mpty Eampl: A matri with (0,-1,1) lmnts, whr no mor than 2 non-zro lmnts in ach column and a 0 if column contains two non-zro lmnts, thn A is TU. i i, = Eampl: A is TU if dtrminant of ach submatri of A is (0, -1,1), and A is (0,1,-1) matri. Duality Primal problm: n+ z LP = ma{ c A b, R } Dual problm: m+ w LP = min{ ub ua c, u R } Wak duality * is primal fasibl, u * is dual fasibl, thn c * * zlp wlp u b Strong Duality If z LP or w LP is finit thn both P and D hav finit optimal valu z LP = w LP Wak dual of IP is any minimization problm zd = min{ zd( u) : u SD} that satisfis n zd( u) c, S, u SD, whrs = { Z, A b} and S D is dual intgr problm. If Dual problm is fasibl thn zip zd. Strong dual of IP is a wak dual that also satisfis if S 0, boundd from abov, 0 0 0 0 u SD, S, zd( u ) = c Duality gap: for wak duals of IP, th sparation btwn primal and dual obctiv valus. zip Primal solution Wak duality Dual solution Primal solution Strong Duality

Totally Dual Intgrality Givn A and b with intgral ntris, w want to study cass whn th LP Ma c A <= b Has intgr optimum primal and intgr optimum dual solutions. Eampls: a) ntwork flow ma c A =b >=0 whr ach column of A has two non-zro ntris. On with valu 1 and othr with valu -1. ntwork flow matri b) ntwork flows with uppr bounds 0<=<=u c) flow problms in cross-fr familis [Nmhausr & Wolsly ttbook] d) submodular flows i. Dirctd cut covrs, matroids ( u) ( S ) = 1, u 1, S odd Prfct Matching problm: Th 2 nd inquality is ncssary in ordr to obtain intgral solutions. Without this inquality, fractional solutions ar possibl. Considr th graph : With only 1 st constraint you ll gt a possibl fractional solution which looks lik: 4 6 8 1 2 5 3 7 9 Edgs shown hav 0.5 valus in solution. Dashd arrows ar th dgs which must b usd in inquality (2) to cut off this fractional solution. 0.5 on ach dg

But how do w automatically dtrmin th form of inquality (2) to cut off th currnt fractional solution? This is calld th Sparation Problm. For som problms lik matching, w know how to solv this problm using Branch and cut. For othr problms, w still do not know how to solv it, although som huristics ist but thy ar not guarantd to obtain a cut. For matching problm, w r givn a fractional vctor, w want on st S such that < 1, for S odd. ( S ) So for a st S, w can us th min cut problm: Givn a graph G=(V,E) and dg wights, w 0, w know how to solv th minimum cut problm: Min s w (S ) This can b solvd by ntwork flow tchniqus. If S is not odd, what do w do? Supos w solv a minimum cut problm with wights w =. Lt S b a solution. If S is odd, w ar don. If S is vn w hav to prov that thr is an odd st in S, includd in S or in V\S that solvs our problm. Using this proprty, w can dcompos th problm. Proof: W know that th following is tru. S T w ( S T ) T) o T long dg counts only in th right hand sid of th inquality. Suppos that S ( S is vn) givs a minimum cut and T givs a minimum cut with T odd. Now w ( S T ) T) is o+<=+o If w hav that w( ), w( T ) S T ), bcaus T and S T ar odd. This st of inqualitis imply w ( S T ) = w( T ) and w( = w( ), w( T ) = w( S T )

Proof. Thrfor S T is also a solution for th minimum odd cut problm. So tak union S T or tak complmnt of union V\ S T (that is includd in V\. Othr possibility: w ( S T ) T ) is + o + o. o o T Thn w ( S T ) = w( T ),bcaus w hav w( S T ), w( T) ) so thrfor w( T ) = w( ) Proof. In this cas is also a solution to th minimum odd cut problm. = if S dfins a minimum cut & S vn, thn solution ists in S or in V\S. W dcompos problm into two problms in th graphs S V\S W gnrat a nstd st of problms, that has at most 2n-1 lmnts. So in worst cas w hav to solv 2n-1 min cut problms. Finding on minimum cut taks O(n^3), so this is an O(n^4) algorithm. ( u) 1, u Non-Prfct Matching problm: S 1, S odd E ( S ) 2 whr E( ar th dgs with both nds in S. This can also b solvd similar to th prfct matching problm, using mincut algorithm. Othr problms, w know what som of th (facts) inqualitis look lik, but w can t solv th sparation problm: Nod or Vrt Packing Problms: Givn a graph G=(V,E) considr u + v 1 ( u, v) E 0 u V u

for dfinition of th problm. Howvr for obtaining intgral solutions w nd to add mor inqualitis (spcifically facts, which intrsct intgr vrtics of sarch spac). 1 K (if subgraph is complt) u V u C 1 u C (whr C is an odd cycl) u C 2 Howvr thr ar mor facts which ar difficult to obtain and for which w don t know all of thm. Thr ar lifting lmmas on can us to obtain mor facts of this problm. P = { a N b,0 Knapsack problm: 1},0 a b, N w nd facts of P I 0, 1 N dfin facts Lt b an intgr vctor in P, w say that st S = { = 1} is indpndnt. Lt C b a minimal dpndnt st C\{i} is indpndnt i C Th inquality C 1 (1) C is valid for P I Assum that a1 a2... an ar givn. A minimal dpndnt st C, dfin E(C) = C { k : ak a, C} Th inquality C 1 (2) E( C) is also valid for P I and strongr than (1). Th inquality (2) dfins a fact of P I undr any of th following conditions: a) C=N b) E(C)=N & (i) C\{ 1, 2 } {1 } is indpndnt c) C=E(C) and (ii) C\{ 1 } {p} is indpndnt, p = min{ N \ E( C)} d) C E( C) N & (i) & (ii) Whr C={ 1, 2,, k } EXAMPLE: P 79 + 53 + 53 + 45 + 45 178,0 1} = { 1 2 3 4 5 C 1 ={1,2,3} C 2 ={1,2,4,5} C 3 ={1,3,4,5} C 4 ={2,3,4,5} ar minimal dpndnt sts Now, E( C 1 )= C 1 E( C 2 )= C 2 E( C 3 )= C 3 E( C 4 )= C 4 {1 }

Lts s if any of th nw inqualitis of th form (2) ar facts: E( C 1 ): (1) 1 + 2 + 3 2 is a fact if (ii) C\{ 1 } {p} = {2,3,4} is indpndnt, p = min{ N \ E( C)} Thrfor it s a fact *. E(C 2 ): (2) 1 + 2 + 4 + 5 3 is a fact if (ii) C\{ 1 } {p} = {2,4,5,3} is indpndnt. It is not indpndnt, thrfor w do not know whthr it is a fact or not. E(C 3 ): (3) 1 + 3 + 4 + 5 3 is a fact if (ii) C\{ 1 } {p} = {3,4,5,2} is indpndnt. It is not indpndnt, thrfor w do not know whthr it is a fact or not. E(C 4 ): (4) 1 + 2 + 3 + 4 + 5 3 is a fact if (cas b)) (i) C\{ 1, 2 } {1 } = C\{2,3} {1 } = {4,5,1} is indpndnt. Is is indpndnt, thrfor it is a fact *. NOTE: any inquality which is a linar combination of facts, cannot b a fact. (3) = (4) + ( 2 0 ), thrfor (3) is not a fact. (bcaus it s a +v combination of 2 facts). (2) = (4) + ( 3 0 ), thrfor (2) dos not dfin a fact.