Crazy Picture. Maximum matching and simplex. z y x Maximum - - PowerPoint PPT Presentation

crazy picture maximum matching and simplex
SMART_READER_LITE
LIVE PREVIEW

Crazy Picture. Maximum matching and simplex. z y x Maximum - - PowerPoint PPT Presentation

Crazy Picture. Maximum matching and simplex. z y x Maximum matching and simplex. max x + y + z x 1 z x + z 1 z + y 1 x y 1 z x 0 y 0 y y z 0 x Maximum matching and simplex. max x + y + z x 1 z x + z 1


slide-1
SLIDE 1

Crazy Picture.

slide-2
SLIDE 2

Maximum matching and simplex.

x y z

slide-3
SLIDE 3

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0 x y z

slide-4
SLIDE 4

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints intersect.

x y z

slide-5
SLIDE 5

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints intersect.

x y z

slide-6
SLIDE 6

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints intersect.

x y z

slide-7
SLIDE 7

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints intersect.

x y z

slide-8
SLIDE 8

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints intersect.

x y z

slide-9
SLIDE 9

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints intersect.

x y z

slide-10
SLIDE 10

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints intersect.

x y z

slide-11
SLIDE 11

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints intersect.

x y z

slide-12
SLIDE 12

Maximum matching and simplex.

y z x maxx +y +z x ≤ 1 x +z ≤ 1 z +y ≤ 1 y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints redundant.

x y z

slide-13
SLIDE 13

Maximum matching and simplex.

y z x maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0 x y z

slide-14
SLIDE 14

Maximum matching and simplex.

y z x maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0 x y z

slide-15
SLIDE 15

Maximum matching and simplex.

maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-16
SLIDE 16

Maximum matching and simplex.

maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-17
SLIDE 17

Maximum matching and simplex.

.3 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-18
SLIDE 18

Maximum matching and simplex.

.7 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-19
SLIDE 19

Maximum matching and simplex.

1 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-20
SLIDE 20

Maximum matching and simplex.

1 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-21
SLIDE 21

Maximum matching and simplex.

1 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-22
SLIDE 22

Maximum matching and simplex.

1 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-23
SLIDE 23

Maximum matching and simplex.

1 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X

slide-24
SLIDE 24

Maximum matching and simplex.

.3 .7 .3 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path.

slide-25
SLIDE 25

Maximum matching and simplex.

.7 .3 .7 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-26
SLIDE 26

Maximum matching and simplex.

1 1 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-27
SLIDE 27

Maximum matching and simplex.

1 1 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-28
SLIDE 28

Maximum matching and simplex.

1 1 maxx +y +z x +z ≤ 1 z +y ≤ 1 x ≥ 0 y ≥ 0 z ≥ 0

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-29
SLIDE 29

Maximum matching and simplex.

1 1 maxx +y +z x +z ≤ 1 a z +y ≤ 1 b x ≥ 0 y ≥ 0 z ≥ 0 c

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-30
SLIDE 30

Maximum matching and simplex.

1 1 1 1 maxx +y +z x +z ≤ 1 a = 1 z +y ≤ 1 b = 1 x ≥ 0 y ≥ 0 z ≥ 0 c = 1

Blue constraints tight. Sum: x +2z +y.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-31
SLIDE 31

Maximum matching and simplex.

1 1 1 1 maxx +y +z x +z ≤ 1 a = 1 z +y ≤ 1 b = 1 x ≥ 0 y ≥ 0 z ≥ 0 c = 1

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-32
SLIDE 32

Maximum matching and simplex.

1 1 1 1 maxx +y +z x +z ≤ 1 a = 1 z +y ≤ 1 b = 1 x ≥ 0 y ≥ 0 z ≥ 0 c = 1

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-33
SLIDE 33

Maximum matching and simplex.

1 1 1 1 maxx +y +z x +z ≤ 1 a = 1 z +y ≤ 1 b = 1 x ≥ 0 y ≥ 0 z ≥ 0 c = 1

Blue constraints tight.

x y z X +1 −1 +1 Augmenting Path. Via Gaussian Elimination!

slide-34
SLIDE 34
slide-35
SLIDE 35

Convex Separator.

slide-36
SLIDE 36

Convex Separator. Farkas

slide-37
SLIDE 37

Convex Separator. Farkas Strong Duality!!!!!

slide-38
SLIDE 38

Convex Separator. Farkas Strong Duality!!!!! Maybe.

slide-39
SLIDE 39

Linear Equations.

Ax = b

slide-40
SLIDE 40

Linear Equations.

Ax = b A is n ×n matrix...

slide-41
SLIDE 41

Linear Equations.

Ax = b A is n ×n matrix... ..has a solution.

slide-42
SLIDE 42

Linear Equations.

Ax = b A is n ×n matrix... ..has a solution. If rows of A are linearly independent.

slide-43
SLIDE 43

Linear Equations.

Ax = b A is n ×n matrix... ..has a solution. If rows of A are linearly independent. yT A = 0 for any y

slide-44
SLIDE 44

Linear Equations.

Ax = b A is n ×n matrix... ..has a solution. If rows of A are linearly independent. yT A = 0 for any y ..or if b in subspace of A.

slide-45
SLIDE 45

Linear Equations.

Ax = b A is n ×n matrix... ..has a solution. If rows of A are linearly independent. yT A = 0 for any y ..or if b in subspace of A. x1 x2 x3

  • k b

bad b

slide-46
SLIDE 46

Strong Duality.

slide-47
SLIDE 47

Strong Duality.

Later.

slide-48
SLIDE 48

Strong Duality.

  • Later. Actually. No.
slide-49
SLIDE 49

Strong Duality.

  • Later. Actually. No. Now
slide-50
SLIDE 50

Strong Duality.

  • Later. Actually. No. Now ...ish.

Special Cases:

slide-51
SLIDE 51

Strong Duality.

  • Later. Actually. No. Now ...ish.

Special Cases: min-max 2 person games and experts.

slide-52
SLIDE 52

Strong Duality.

  • Later. Actually. No. Now ...ish.

Special Cases: min-max 2 person games and experts. Max weight matching and algorithm.

slide-53
SLIDE 53

Strong Duality.

  • Later. Actually. No. Now ...ish.

Special Cases: min-max 2 person games and experts. Max weight matching and algorithm. Approximate: facility location primal dual.

slide-54
SLIDE 54

Strong Duality.

  • Later. Actually. No. Now ...ish.

Special Cases: min-max 2 person games and experts. Max weight matching and algorithm. Approximate: facility location primal dual. Today: Geometry!

slide-55
SLIDE 55

Convex Body and point.

For a convex body P and a point b, b ∈ P or hyperplane separates P from b.

slide-56
SLIDE 56

Convex Body and point.

For a convex body P and a point b, b ∈ P or hyperplane separates P from b. v,α, where v ·x ≤ α and v ·b > α.

slide-57
SLIDE 57

Convex Body and point.

For a convex body P and a point b, b ∈ P or hyperplane separates P from b. v,α, where v ·x ≤ α and v ·b > α.

slide-58
SLIDE 58

Convex Body and point.

For a convex body P and a point b, b ∈ P or hyperplane separates P from b. v,α, where v ·x ≤ α and v ·b > α. point p where (x −p)T (b −p) < 0 b p x

slide-59
SLIDE 59

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0

slide-60
SLIDE 60

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0 b p x Proof: Choose p to be closest point to b in P.

slide-61
SLIDE 61

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0 b p x Proof: Choose p to be closest point to b in P. Done

slide-62
SLIDE 62

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0 b p x Proof: Choose p to be closest point to b in P. Done or ∃ x ∈ P with (x −p)T (b −p) ≥ 0

slide-63
SLIDE 63

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0 b p x Proof: Choose p to be closest point to b in P. Done or ∃ x ∈ P with (x −p)T (b −p) ≥ 0 b p x P (x −p)T (b −p) ≥ 0

slide-64
SLIDE 64

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0 b p x Proof: Choose p to be closest point to b in P. Done or ∃ x ∈ P with (x −p)T (b −p) ≥ 0 b p x P (x −p)T (b −p) ≥ 0 → ≤ 90◦ angle between − − − → x −p and − − − → b −p.

slide-65
SLIDE 65

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0 b p x Proof: Choose p to be closest point to b in P. Done or ∃ x ∈ P with (x −p)T (b −p) ≥ 0 b p x P (x −p)T (b −p) ≥ 0 → ≤ 90◦ angle between − − − → x −p and − − − → b −p. Must be closer point b on line from p to x.

slide-66
SLIDE 66

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0 b p x Proof: Choose p to be closest point to b in P. Done or ∃ x ∈ P with (x −p)T (b −p) ≥ 0 b p x P (x −p)T (b −p) ≥ 0 → ≤ 90◦ angle between − − − → x −p and − − − → b −p. Must be closer point b on line from p to x. All points on line to x are in polytope.

slide-67
SLIDE 67

Proof.

For a convex body P and a point b, b ∈ A or there is point p where (x −p)T (b −p) < 0 b p x Proof: Choose p to be closest point to b in P. Done or ∃ x ∈ P with (x −p)T (b −p) ≥ 0 b p x P (x −p)T (b −p) ≥ 0 → ≤ 90◦ angle between − − − → x −p and − − − → b −p. Must be closer point b on line from p to x. All points on line to x are in polytope. Contradicts choice of p as closest point to b in polytope.

slide-68
SLIDE 68

More formally.

b p x P Squared distance to b from p +(x −p)µ

slide-69
SLIDE 69

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x

slide-70
SLIDE 70

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2

slide-71
SLIDE 71

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2 θ is the angle between x −p and b −p.

slide-72
SLIDE 72

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2 θ is the angle between x −p and b −p. b |p −b|−ℓcosθ Distance to new point. p θ p + µ(x −p) ℓ = µ|x −p| ℓsinθ ℓcosθ x

slide-73
SLIDE 73

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2 θ is the angle between x −p and b −p. b |p −b|−ℓcosθ Distance to new point. p θ p + µ(x −p) ℓ = µ|x −p| ℓsinθ ℓcosθ x Simplify:

slide-74
SLIDE 74

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2 θ is the angle between x −p and b −p. b |p −b|−ℓcosθ Distance to new point. p θ p + µ(x −p) ℓ = µ|x −p| ℓsinθ ℓcosθ x Simplify: |p −b|2 −2µ|p −b||x −p|cosθ +(µ|x −p|)2.

slide-75
SLIDE 75

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2 θ is the angle between x −p and b −p. b |p −b|−ℓcosθ Distance to new point. p θ p + µ(x −p) ℓ = µ|x −p| ℓsinθ ℓcosθ x Simplify: |p −b|2 −2µ|p −b||x −p|cosθ +(µ|x −p|)2. Derivative with respect to µ ...

slide-76
SLIDE 76

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2 θ is the angle between x −p and b −p. b |p −b|−ℓcosθ Distance to new point. p θ p + µ(x −p) ℓ = µ|x −p| ℓsinθ ℓcosθ x Simplify: |p −b|2 −2µ|p −b||x −p|cosθ +(µ|x −p|)2. Derivative with respect to µ ... −2|p −b||x −p|cosθ +2(µ|x −p|2).

slide-77
SLIDE 77

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2 θ is the angle between x −p and b −p. b |p −b|−ℓcosθ Distance to new point. p θ p + µ(x −p) ℓ = µ|x −p| ℓsinθ ℓcosθ x Simplify: |p −b|2 −2µ|p −b||x −p|cosθ +(µ|x −p|)2. Derivative with respect to µ ... −2|p −b||x −p|cosθ +2(µ|x −p|2). which is negative for a small enough value of µ

slide-78
SLIDE 78

More formally.

b p x P Squared distance to b from p +(x −p)µ point between p and x (|p −b|− µ|x −p|cosθ)2 +(µ|x −p|sinθ)2 θ is the angle between x −p and b −p. b |p −b|−ℓcosθ Distance to new point. p θ p + µ(x −p) ℓ = µ|x −p| ℓsinθ ℓcosθ x Simplify: |p −b|2 −2µ|p −b||x −p|cosθ +(µ|x −p|)2. Derivative with respect to µ ... −2|p −b||x −p|cosθ +2(µ|x −p|2). which is negative for a small enough value of µ (for positive cosθ.)

slide-79
SLIDE 79

Generalization: exercise.

Theorems of Alternatives.

slide-80
SLIDE 80

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it.

slide-81
SLIDE 81

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5.

slide-82
SLIDE 82

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5. That is, find y where yT A = 0 and yT b = 0.

slide-83
SLIDE 83

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5. That is, find y where yT A = 0 and yT b = 0. Space is image of A.

slide-84
SLIDE 84

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5. That is, find y where yT A = 0 and yT b = 0. Space is image of A. Affine subspace is columnspan of A.

slide-85
SLIDE 85

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5. That is, find y where yT A = 0 and yT b = 0. Space is image of A. Affine subspace is columnspan of A. y is normal.

slide-86
SLIDE 86

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5. That is, find y where yT A = 0 and yT b = 0. Space is image of A. Affine subspace is columnspan of A. y is normal. y in nullspace for column span.

slide-87
SLIDE 87

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5. That is, find y where yT A = 0 and yT b = 0. Space is image of A. Affine subspace is columnspan of A. y is normal. y in nullspace for column span. yT b = 0 = ⇒ b not in column span.

slide-88
SLIDE 88

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5. That is, find y where yT A = 0 and yT b = 0. Space is image of A. Affine subspace is columnspan of A. y is normal. y in nullspace for column span. yT b = 0 = ⇒ b not in column span. There is a separating hyperplane between any two convex bodies.

slide-89
SLIDE 89

Generalization: exercise.

Theorems of Alternatives. Linear Equations: There is a separating hyperplane between a point and an affine subspace not containing it. From Ax = b use row reduction to get, e.g., 0 = 5. That is, find y where yT A = 0 and yT b = 0. Space is image of A. Affine subspace is columnspan of A. y is normal. y in nullspace for column span. yT b = 0 = ⇒ b not in column span. There is a separating hyperplane between any two convex bodies. Idea: Let closest pair of points in two bodies define direction.

slide-90
SLIDE 90

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • 1

1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2

slide-91
SLIDE 91

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • 1

1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2

slide-92
SLIDE 92

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • 1

1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2

slide-93
SLIDE 93

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • 1

1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2

slide-94
SLIDE 94

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • 1

1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2

slide-95
SLIDE 95

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • −1

−1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2

slide-96
SLIDE 96

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • −1

−1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2

slide-97
SLIDE 97

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • −1

−1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2

slide-98
SLIDE 98

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • −1

−1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2 y where yT (b−Ax) < yT (0) < 0 for all x ≥ 0

slide-99
SLIDE 99

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • −1

−1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2 y where yT (b−Ax) < yT (0) < 0 for all x ≥ 0 → yT b < 0 and yT A ≥ 0.

slide-100
SLIDE 100

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • −1

−1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2 y where yT (b−Ax) < yT (0) < 0 for all x ≥ 0 → yT b < 0 and yT A ≥ 0. Farkas A: Solution for exactly one of:

slide-101
SLIDE 101

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • −1

−1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2 y where yT (b−Ax) < yT (0) < 0 for all x ≥ 0 → yT b < 0 and yT A ≥ 0. Farkas A: Solution for exactly one of: (1) Ax = b,x ≥ 0

slide-102
SLIDE 102

Ax = b, x ≥ 0

  • 1

1 1 1

  • x =
  • −1

−1

  • x1

x2 x3 Coordinates s = b −Ax. x ≥ 0 where s = 0? s1 s2 y where yT (b−Ax) < yT (0) < 0 for all x ≥ 0 → yT b < 0 and yT A ≥ 0. Farkas A: Solution for exactly one of: (1) Ax = b,x ≥ 0 (2) yT A ≥ 0,yT b < 0.

slide-103
SLIDE 103

Farkas 2

Farkas A: Solution for exactly one of:

slide-104
SLIDE 104

Farkas 2

Farkas A: Solution for exactly one of: (1) Ax = b,x ≥ 0

slide-105
SLIDE 105

Farkas 2

Farkas A: Solution for exactly one of: (1) Ax = b,x ≥ 0 (2) yT A ≥ 0,yT b < 0.

slide-106
SLIDE 106

Farkas 2

Farkas A: Solution for exactly one of: (1) Ax = b,x ≥ 0 (2) yT A ≥ 0,yT b < 0. Farkas B: Solution for exactly one of:

slide-107
SLIDE 107

Farkas 2

Farkas A: Solution for exactly one of: (1) Ax = b,x ≥ 0 (2) yT A ≥ 0,yT b < 0. Farkas B: Solution for exactly one of: (1) Ax ≤ b

slide-108
SLIDE 108

Farkas 2

Farkas A: Solution for exactly one of: (1) Ax = b,x ≥ 0 (2) yT A ≥ 0,yT b < 0. Farkas B: Solution for exactly one of: (1) Ax ≤ b (2) yT A = 0,yT b < 0,y ≥ 0.

slide-109
SLIDE 109

Strong Duality

(From Goemans notes.) Primal P z∗ = mincT x Ax = b x ≥ 0 Dual D :w∗ = maxbT y AT y ≤ c

slide-110
SLIDE 110

Strong Duality

(From Goemans notes.) Primal P z∗ = mincT x Ax = b x ≥ 0 Dual D :w∗ = maxbT y AT y ≤ c Weak Duality: x,y- feasible P , D: xT c ≥ bT y.

slide-111
SLIDE 111

Strong Duality

(From Goemans notes.) Primal P z∗ = mincT x Ax = b x ≥ 0 Dual D :w∗ = maxbT y AT y ≤ c Weak Duality: x,y- feasible P , D: xT c ≥ bT y. xT c −bT y = xT c −xT AT y = xT (c −AT y) ≥ 0

slide-112
SLIDE 112

Strong duality If P or D is feasible and bounded then z∗ = w∗.

slide-113
SLIDE 113

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗.

slide-114
SLIDE 114

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗.

slide-115
SLIDE 115

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗.

slide-116
SLIDE 116

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

slide-117
SLIDE 117

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .
slide-118
SLIDE 118

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0
slide-119
SLIDE 119

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0

slide-120
SLIDE 120

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0.

slide-121
SLIDE 121

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗.

slide-122
SLIDE 122

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗. Better Primal!!

slide-123
SLIDE 123

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗. Better Primal!!

Case 2: λ = 0. Ax = 0,cT x < 0.

slide-124
SLIDE 124

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗. Better Primal!!

Case 2: λ = 0. Ax = 0,cT x < 0. Feasible ˜ x for Primal.

slide-125
SLIDE 125

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗. Better Primal!!

Case 2: λ = 0. Ax = 0,cT x < 0. Feasible ˜ x for Primal. (a) ˜ x + µx ≥ 0 since ˜ x,x,µ ≥ 0.

slide-126
SLIDE 126

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗. Better Primal!!

Case 2: λ = 0. Ax = 0,cT x < 0. Feasible ˜ x for Primal. (a) ˜ x + µx ≥ 0 since ˜ x,x,µ ≥ 0. (b) A(˜ x + µx) = A˜ x + µAx = b + µ ·0 = b.

slide-127
SLIDE 127

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗. Better Primal!!

Case 2: λ = 0. Ax = 0,cT x < 0. Feasible ˜ x for Primal. (a) ˜ x + µx ≥ 0 since ˜ x,x,µ ≥ 0. (b) A(˜ x + µx) = A˜ x + µAx = b + µ ·0 = b. Feasible

slide-128
SLIDE 128

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗. Better Primal!!

Case 2: λ = 0. Ax = 0,cT x < 0. Feasible ˜ x for Primal. (a) ˜ x + µx ≥ 0 since ˜ x,x,µ ≥ 0. (b) A(˜ x + µx) = A˜ x + µAx = b + µ ·0 = b. Feasible cT (˜ x + µx) = xT ˜ x + µcT x → −∞ as µ → ∞

slide-129
SLIDE 129

Strong duality If P or D is feasible and bounded then z∗ = w∗. Primal feasible, bounded, value z∗. Claim: Exists a solution to dual of value at least z∗. ∃y,yT A ≤ c,bT y ≥ z∗. Want y.

  • AT

−bT

  • y ≤
  • c

−z∗

  • .

If none, then Farkas B says ∃x,λ ≥ 0.

  • A

−b

  • x

λ

  • = 0
  • cT

−z∗ x λ

  • < 0

∃x,λ with Ax −bλ = 0 and ctx −z∗λ < 0 Case 1: λ > 0. A( x

λ ) = b, cT ( x λ ) < z∗. Better Primal!!

Case 2: λ = 0. Ax = 0,cT x < 0. Feasible ˜ x for Primal. (a) ˜ x + µx ≥ 0 since ˜ x,x,µ ≥ 0. (b) A(˜ x + µx) = A˜ x + µAx = b + µ ·0 = b. Feasible cT (˜ x + µx) = xT ˜ x + µcT x → −∞ as µ → ∞ Primal unbounded!

slide-130
SLIDE 130

See you on Tuesday.