This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: 15854: Approximations AlgorithmsLecturer:R. RaviTopic:PrimalDual AlgorithmsDate:Nov. 21, 2005Scribe:Daniel Golovin21.1PrimalDual AlgorithmsSo far, we have seen many algorithms based on linear program (LP) relaxations, typically involvingrounding a given fractional LP solution to an integral solution of approximately the same objectivevalue. In this lecture, we will look at another approach to LP relaxations, in which we will constructa feasible integral solution to the LP from scratch, using a related LP to guide our decisions. OurLP will be called thePrimal LP, and the guiding LP will be called theDual LP.As we shall see, the PD method is quite powerful. Often, we can use the PrimalDual (PD) methodto obtain a good approximation algorithm, and then extract a good combinatorial algorithm fromit. Conversely, sometimes we can use the PD method to prove good performance for combinatorialalgorithms, simply by reinterpreting them as PD algorithms. So without further ado...21.2Every Primal has a DualWe begin with a generic covering LP, and illustrate the ideas later with Vertex Cover as an example.Let [k] :={1,2,...,k}. Suppose we have matrixA∈Rm×nand vectorsc∈Rn,b∈Rm. We canrepresent the primal LP asmin∑jcjxjsubject to∑jaijxj≥bi∀i∈[m]xj≥∀j∈[n](Primal)Now suppose we want to develop a lower bound on the optimal value of this LP. One way to dothis is to find constraints that “look like”∑jcjxj≥Z, for someZ, using the constraints in the LP.To do this, note that any convex combination of constraints from the LP is also a valid constraint.Therefore, if we have nonnegative multipliersyion the constraints, we get a new constraint whichis satisfied by all feasible solutions to the primal LP. That is, if for alli,∑jaijxj≥bi, thenXiyiXjaijxj≥Xiyibi(21.2.1)Note that we require theyi’s to be nonnegative, because multiplying an inequality (in this case∑jaijxj≥bi) by a negative number switches the sign of the inequality. (If a constraint has theform∑jaijxj=bi, then its multiplieryican be any real number.) Consider equation 21.2.1. If weensure∑iyi∑jaijxj≤∑jcjxj, we will obtain a lower bound of∑iyibion the optimal valueof the primal LP. Switching the order of summation, we get∑iyi∑jaijxj=∑j(∑iyiaij)xj,1and can ensure this sum is at most∑jcjxjby requiring theyi’s to satisfyXiyiaij≤cj∀j∈[n](21.2.2)(Note that in the previous step we rely on the fact that thexj’s are nonnegative.)Putting it all together, if theyi’s are nonnegative and satisfy constraint 21.2.2, thenXiyibi≤XiyiXjaijxj=XjXiyiaij!xj≤Xjcjxj(21.2.3)Note that the constraints on theyi’s are linear, as is the lower bound we obtain. Thus we can writedown an LP to find theyi’s giving the best lower bound. This is thedual LP....
View
Full
Document
This note was uploaded on 01/26/2010 for the course COMPUTER S 15853 taught by Professor Guyblelloch during the Fall '09 term at Carnegie Mellon.
 Fall '09
 GuyBlelloch
 Algorithms

Click to edit the document details