K-Medians, Facility Location, and the Chernoff-Wald Bound

Total Page:16

File Type:pdf, Size:1020Kb

K-Medians, Facility Location, and the Chernoff-Wald Bound 86 K-Medians, Facility Location, and the Chernoff-Wald Bound Neal E. Young* Abstract minimizez d + k We study the general (non-metric) facility-location and weighted k-medians problems, as well as the fractional cost(x) < k facility-location and k-medians problems. We describe a dist(x) S d natural randomized rounding scheme and use it to derive approximation algorithms for all of these problems. subject to ~-~/x(f,c) = 1 (Vc) For facility location and weighted k-medians, the re- =(f,c) < x(f) (vf, c) spective algorithms are polynomial-time [HAk + d]- and [(1 + e)d; In(n + n/e)k]-approximation algorithms. These x(f,c) >_ 0 (vf, c) performance guarantees improve on the best previous per- x(f) E {0,1} (vf) formance guarantees, due respectively to Hochbaum (1982) and Lin and Vitter (1992). For fractional k-medians, the al- gorithm is a new, Lagrangian-relaxation, [(1 + e)d, (1 + e)k]- Figure 1: The facility-location IP. Above d, k, x(f) and approximation algorithm. It runs in O(kln(n/e)/e 2) linear- x(f, c) are variables, dist(x), the assignment cost of x, time iterations. is ~-~Ic x(f, c)disz(f, c), and cost(x), the facility cost of For fractional facilities-location (a generalization of frac- tional weighted set cover), the algorithm is a Lagrangian- relaxation, ~(1 + e)k]-approximation algorithm. It runs in O(nln(n)/e ~) linear-time iterations and is essentially the The weighted k-medians problem is the same same as an unpublished Lagrangian-relaxation algorithm as the facility location problem except for the following: due to Garg (1998). By recasting his analysis probabilis- a positive real number k is given as input, and the goal tically and abstracting it, we obtain an interesting (and as far as we know new) probabilistic bound that may be of is to choose a subset F of facilities minimizing dist(F) independent interest. We call it the Chernoff- Wald bound. subject to the constraint cost(F) < k. The standard integer programming formulation is the k-medians IP, 1 Problem definitions which differs from the IP in Fig. 1 only in that k is given, not a variable, and the objective function is d instead The input to the weighted set cover problem is of d + k. The (unweighted) k-medians problem is the a collection of sets, where each set s is given a cost special case when each cost(f) = 1. For the fractional cost(s) E ~+. The goal is to choose a cover (a collection k-medians problem, the input is the same; the goal is of sets containing all elements) of minimum total cost. to solve the linear program obtained by removing the The (uncapacitated) facility-location problem is constraint "x(f) E {0, 1}" from the k-medians IP. a generalization of weighted set cover in which each We take the "size" of each of the above problems to set f (called a "facility") and element c (called a be the number of pairs (f, c) such that d±st(f, c) < oo. "customer") are given a distance dist(f, c) E ~+U{oo}. We assume the size is at least the number of customers The goal is to choose a set of facilities F minimizing and facilities. cost(F) + dist(F), where cost(F), the facility cost of By an In(d); fl(k)]-approximation algorithm for F, is ~"~.EF COSt(f) and dist(F), the assignment cost k-medians, we mean an algorithm that, given a problem of F, is )-:c rainier d±st(f, c). instance for which there exists a fractional solution of Fig. 1 shows the standard integer programming assignment cost d and facility cost k, produces a solution formulation of the problem -- the facility-location IP of assignment cost at most a(d) and facility cost at most [12, p. 8]. The facility-location linear program (LP) is D(k). We use similar non-standard notations for facility the same except without the constraint "x(f) E {0, 1}". location, set cover, and k-set cover. For instance, by A fractional solution is a feasible solution to the LP. a [d + 2k]-approximation algorithm for facility location, Fractional facility location is the problem of solving the we mean an algorithm that, given an instance for which LP. there exists a fractional solution of assignment cost d -'--'¢R-~'earch partially funded by NSF CAREER award CCR- and facility cost k, produces a solution for which the 9720664. Department of Computer Science, Dartmouth College, assignment cost plus the facility cost is at most d + 2k. Hanover, NH. ney~}cs,daxtmouZh, edu 87 2 Background 3 Results In the mid 1970's Johnson and Lovasz gave a greedy Fig. 2 shows the simple randomized rounding scheme [HAk]-approximation algorithm for unweighted set at the center of all our results. With minor variations cover [9, 11]. In 1979 Chvatal generalized it to a [HAk]- (e.g. Fig. 4), this rounding scheme can be used as approximation algorithm for weighted set cover [3]. the basis for approximation algorithms for set cover, In 1982 Hochbaum gave a greedy [HA(d + k)]- weighted set cover, facility location, and k-medians, and approximation algorithm for the uncapacitated facility- as the basis for Lagrangian-relaxation algorithms for the location problem by an implicit reduction to the fractional variants of these problems. weighted set-cover problem [8].1 Above /x is at most Although essentially the same rounding scheme suf- the maximum, over all facilities f, of the number of rices for each of these problems, the respective proba- customers c such that dist(f, c) < oo. bilistic analyses require different (albeit standard) tech- In 1992 Lin and Vitter gave a polynomial-time niques in each case. For set cover, a simple direct [(1 + e)d; (1 + 1/e) (lnn + 1) k]-approximation algorithm analysis suffices [16]. For weighted set cover, facility for the k-medians problem [10]. (Here c > 0 is an location, and k-medians, a basic probabilistic lemma input parameter that tunes the tradeoff between the called Wald's inequality is necessary. For fractional set two criteria and n is the number of customers.) Their cover and k-medians, the analysis rests on the Chernoff algorithm combines a greedy algorithm with a technique bound. For fractional weighted set cover and facility lo- they call filtering. cation, the analysis is a simple application of what we In 1994 Plotkin, Shmoys, and Tardos (PST) gave call the Chernoff- Wald bound. Lagrangian-relaxation algorithms for general packing For each problem, we apply the method of condi- and covering problems [13]. As a special case, their tional probabilities to the rounding scheme in order to algorithms imply a [(1 + e)k]-approximation algorithm derive a corresponding approximation algorithm. The for fractional set cover that runs in O(k ln(n)/e 2) linear- structure of each resulting algorithm, being closely tied time iterations. Here n is the number of elements. to the underlying probabilistic analysis, ends up differ- In 1998 Garg generalized and simplified the PST ing substantially from problem to problem. set cover algorithm to obtain a [(1 + e)k]-approximation For facility location, the resulting algorithm (Fig. 2) algorithm for fractional weighted set cover [6]. Garg's is a randomized rounding, polynomial-time [d + H~xk]- algorithm runs in O(nln(n)/e 2) linear-time iterations. approximation algorithm. The performance guarantee By Hochbaum's reduction, one can use Garg's algorithm improves over Hochbaum's algorithm with respect to as a [(1 + e)k]-approximation algorithm for fractional the assignment costs. /acility location. The running time is the same, where For weighted k-medians, the resulting algorithm n is the number of customers. (Fig. 3) is a greedy [(1 + e)d; ln(n +n/e)k]-apprordmation Recent work has focused on metric k-medians and algorithm. In comparison to Lin and Vitter's algorithm, facility location problems. In the metric versions, the the performance ratio with respect to the facility costs distance function is assumed to satisfy the triangle in- is better by a factor of roughly 1/e. equality. For example, [O(d + k)]-approximation algo- For fractional k-medians, the algorithm (Fig. 5) rithms have recently been shown for the metric facility- is a [(1 + e)d, (1 + e)k]-approximation algorithm. It location problem [7, 15]. Charikar, Guha, Tardos and is a Lagrangian-relaxation algorithm and runs in Shmoys [2] recently gave an [O(d); k]-approximation al- O(k ln(n/e)/e 2) linear-time iterations. This is a factor gorithm for metric k-medians. Many of these algorithm of n faster than the best bound we can show by applying first solve the fractional problems and then round the the general algorithm of Plotkin, Shmoys and Tardos. fractional solutions. Finally, for fractional facility location, the algo- rithm (see Fig. 6) is a [(1 + e)(d + k)]-approximation, rHHochbaum's reduction is easy to adapt in order to reduce Lagrangian-relaxation algorithm. It runs in at most unweighted k-medians to a variant of set cover that we call k- O(nln(n)/e 2) iterations, where each iteration requires set cover. The reduction also extends naturally to the fractional time linear in the input size times Inn. This algorithm problems. (See appendix for details.) This shows a loose is the same as the unpublished fractional weighted set- equivalence between facility location and weighted set cover, and between k-medians and k-set cover. However, Hochbaum's cover algorithm due to Garg [6].
Recommended publications
  • Approximation Algorithms for 2-Stage Stochastic Optimization Problems
    Algorithms Column: Approximation Algorithms for 2-Stage Stochastic Optimization Problems Chaitanya Swamy∗ David B. Shmoys† February 5, 2006 1 Column editor’s note This issue’s column is written by guest columnists, David Shmoys and Chaitanya Swamy. I am delighted that they agreed to write this timely column on the topic related to stochastic optimization that has received much attention recently. Their column introduces the reader to several recent results and provides references for further readings. Samir Khuller 2 Introduction Uncertainty is a facet of many decision environments and might arise for various reasons, such as unpredictable information revealed in the future, or inherent fluctuations caused by noise. Stochastic optimization provides a means to handle uncertainty by modeling it by a probability distribution over possible realizations of the actual data, called scenarios. The field of stochastic optimization, or stochastic programming, has its roots in the work of Dantzig [4] and Beale [1] in the 1950s, and has since increasingly found application in a wide variety of areas, including transportation models, logistics, financial instruments, and network design. An important and widely-used model in stochastic programming is the 2- stage recourse model: first, given only distributional information about (some of) the data, one commits on initial (first-stage) actions. Then, once the actual data is realized according to the distribution, further recourse actions can be taken (in the second stage) to augment the earlier solution and satisfy the revealed requirements. The aim is to choose the initial actions so as to minimize the expected total cost incurred. The recourse actions typically entail making decisions in rapid reaction to the observed scenario, and are therefore more costly than decisions made ahead of time.
    [Show full text]
  • Approximation Algorithms
    In Introduction to Optimization, Decision Support and Search Methodologies, Burke and Kendall (Eds.), Kluwer, 2005. (Invited Survey.) Chapter 18 APPROXIMATION ALGORITHMS Carla P. Gomes Department of Computer Science Cornell University Ithaca, NY, USA Ryan Williams Computer Science Department Carnegie Mellon University Pittsburgh, PA, USA 18.1 INTRODUCTION Most interesting real-world optimization problems are very challenging from a computational point of view. In fact, quite often, finding an optimal or even a near-optimal solution to a large-scale optimization problem may require com- putational resources far beyond what is practically available. There is a sub- stantial body of literature exploring the computational properties of optimiza- tion problems by considering how the computational demands of a solution method grow with the size of the problem instance to be solved (see e.g. Aho et al., 1979). A key distinction is made between problems that require compu- tational resources that grow polynomially with problem size versus those for which the required resources grow exponentially. The former category of prob- lems are called efficiently solvable, whereas problems in the latter category are deemed intractable because the exponential growth in required computational resources renders all but the smallest instances of such problems unsolvable. It has been determined that a large class of common optimization problems are classified as NP-hard. It is widely believed—though not yet proven (Clay Mathematics Institute, 2003)—that NP-hard problems are intractable, which means that there does not exist an efficient algorithm (i.e. one that scales poly- nomially) that is guaranteed to find an optimal solution for such problems.
    [Show full text]
  • David B. Shmoys
    CURRICULUM VITAE David B. Shmoys 231 Rhodes Hall 312 Comstock Road Cornell University Ithaca, NY 14850 Ithaca, NY 14853 (607)-257-7404 (607)-255-9146 [email protected] Education Ph.D., Computer Science, University of California, Berkeley, June 1984. Thesis: Approximation Algorithms for Problems in Scheduling, Sequencing & Network Design Thesis Committee: Eugene L. Lawler (chair), Dorit S. Hochbaum, Richard M. Karp. B.S.E., Electrical Engineering and Computer Science, Princeton University, June 1981. Graduation with Highest Honors, Phi Beta Kappa, Tau Beta Pi. Thesis: Perfect Graphs & the Strong Perfect Graph Conjecture Research Advisors: H. W. Kuhn, K. Steiglitz and D. B. West Professional Employment Laibe/Acheson Professor of Business Management & Leadership Studies, Cornell, 2013–present. Tang Visiting Professor, Department of Ind. Eng. & Operations Research, Columbia, Spring 2019. Visiting Scientist, Simons Institute for the Theory of Computing, University of California, Berkeley, Spring 2018. Director, School of Operations Research & Information Engineering, Cornell University, 2013–2017. Professor of Operations Research & Information Engineering and of Computer Science, Cornell, 2007–present. Visiting Professor, Sloan School of Management, MIT, 2010–2011, Summer 2015. Visiting Research Scientist, New England Research & Development Center, Microsoft, 2010–2011. Professor of Operations Research & Industrial Engineering and of Computer Science, Cornell University, 1997–2006. Visiting Research Scientist, Computer Science Division, EECS, University of California, Berkeley, 1999–2000. Associate Professor of Operations Research & Industrial Engineering and of Computer Science, Cor- nell University, 1996–1997. Associate Professor of Operations Research & Industrial Engineering, Cornell University, 1992–1996. 1 Visiting Research Scientist, Department of Computer Science, Eotv¨ os¨ University, Fall 1993. Visiting Research Scientist, Department of Computer Science, Princeton University, Spring 1993.
    [Show full text]
  • Exact and Approximate Algorithms for Some Combinatorial Problems
    EXACT AND APPROXIMATE ALGORITHMS FOR SOME COMBINATORIAL PROBLEMS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy by Soroush Hosseini Alamdari May 2018 c 2018 Soroush Hosseini Alamdari ALL RIGHTS RESERVED EXACT AND APPROXIMATE ALGORITHMS FOR SOME COMBINATORIAL PROBLEMS Soroush Hosseini Alamdari, Ph.D. Cornell University 2018 Three combinatorial problems are studied and efficient algorithms are pre- sented for each of them. The first problem is concerned with lot-sizing, the second one arises in exam-scheduling, and the third lies on the intersection of the k-median and k-center clustering problems. BIOGRAPHICAL SKETCH Soroush Hosseini Alamdari was born in Malayer, a town in Hamedan province of Iran right on the outskirts of Zagros mountains. After receiving a gold medal in the Iranian National Olympiad in Informatics he went to Sharif University of Technology in Tehran to study Computer Engineering. After completing his BSc he went to the University of Waterloo and received his MMath in computer science. Then after a short period of working in industry, he came back to com- plete his PhD in computer science here at the Cornell University. iii ACKNOWLEDGEMENTS First and foremost I would like to thank my Adviser Prof. David Shmoys who accepted me as a student, taught me so many fascinating things, and supported me in whatever endeavor I took on. Without his care and attention this docu- ment would simply not be possible. I would also like to thank all my friends, colleagues and co-authors.
    [Show full text]
  • Approximation Algorithms for NP-Hard Optimization Problems
    Approximation algorithms for NP-hard optimization problems Philip N. Klein Neal E. Young Department of Computer Science Department of Computer Science Brown University Dartmouth College Chapter 34, Algorithms and Theory of Computation Handbook c 2010 Chapman & Hall/CRC 1 Introduction In this chapter, we discuss approximation algorithms for optimization problems. An optimization problem consists in finding the best (cheapest, heaviest, etc.) element in a large set P, called the feasible region and usually specified implicitly, where the quality of elements of the set are evaluated using a func- tion f(x), the objective function, usually something fairly simple. The element that minimizes (or maximizes) this function is said to be an optimal solution of the objective function at this element is the optimal value. optimal value = minff(x) j x 2 Pg (1) A example of an optimization problem familiar to computer scientists is that of finding a minimum-cost spanning tree of a graph with edge costs. For this problem, the feasible region P, the set over which we optimize, consists of spanning trees; recall that a spanning tree is a set of edges that connect all the vertices but forms no cycles. The value f(T ) of the objective function applied to a spanning tree T is the sum of the costs of the edges in the spanning tree. The minimum-cost spanning tree problem is familiar to computer scientists because there are several good algorithms for solving it | procedures that, for a given graph, quickly determine the minimum-cost spanning tree. No matter what graph is provided as input, the time required for each of these algorithms is guaranteed to be no more than a slowly growing function of the number of vertices n and edges m (e.g.
    [Show full text]
  • Lecture Notes
    Chapter 1 An introduction to approximation algorithms These lecture notes are for the course Advanced Algo- rithms (2016-2017) and are supplementary to the book `The design of approxiation algorithms' by David Shmoys and David Williamson. See the book for missing definitions. These notes only cover a part of the sections of the book. If you find an error/typo please let me know ([email protected]) and I will update de file. Contents 1 An introduction to approximation algorithms 3 1.1 Approximation Algorithms . .3 1.2 Introduction to techniques of LP . .3 1.3 A deterministic rounding algorithm . .6 1.4 Rounding a dual solution . .8 1.5 The primal-dual method . .9 1.6 A greedy algorithm . 11 1.7 A randomized rounding algorithm . 14 2 Greedy and Local search 17 2.1 Single machine scheduling . 17 2.2 The k-center problem. 18 2.3 Parallel machine scheduling. 19 2.4 The traveling salesman problem. 22 3 Rounding data and dynamic programming 26 3.1 The knapsack problem . 27 3.2 Scheduling on identical machines . 29 4 Deterministic rounding of linear programs 34 4.1 Sum of completion times on a single machine. 34 4.2 Weighted sum of completion times. 35 4.3 Using the ellipsoid method . 38 4.4 Prize-collecting Steiner Tree . 39 4.5 Uncapacitated facility location . 41 1 Chapter 1 An introduction to approximation algorithms 5 Random sampling and randomized rounding of LP's 44 5.1 Max Sat and Max Cut . 44 5.2 Derandomization. 47 5.3 Flipping biased coins .
    [Show full text]
  • On the Complexity of Linear Programming
    CHAPTER 6 On the complexity of linear programming Nimrod Megiddo Abstract: This is a partial survey of results on the complexity of the lin- ear programming problem since the ellipsoid method. The main topics are polynomial and strongly polynomial algorithms, probabilistic analy- sis of simplex algorithms, and recent interior point methods. 1 Introduction Our purpose here is to survey theoretical developments in linear program- ming, starting from the ellipsoid method, mainly from the viewpoint of computational comp1exity.l The survey does not attempt to be complete and naturally reflects the author's perspective, which may differ from the viewpoints of others. Linear programming is perhaps the most successful discipline of Oper- ations Research. The standard form of the linear programming problem is to maximize a linear function cTx (c,x E Rn)over all vectors x such that Ax = b and x r 0. We denote such a problem by (A,b, c).Currently, the main tool for solving the linear programming problem in practice is the class of simplex algorithms proposed and developed by Dantzig [43]. However, applications of nonlinear programming methods, inspired by Karmarkar's work [79], may also become practical tools for certain classes of linear programming problems. Complexity-based questions about linear programming and related parameters of polyhedra (see, e.g., [661) have been raised since the 1950s, before the field of computational complexity The author thanks Jetfrey Lagarias, David Shmoys, and Robert Wilson for helpful com- ments. ' The topic of linear programming has been interesting to economists not only due to its applicability to practical economic systems but also because of the many economic in- sights provided by the theory of linear programming.
    [Show full text]
  • The Primal-Dual Method for Approximation Algorithms
    Math. Program., Ser. B 91: 447–478 (2002) Digital Object Identifier (DOI) 10.1007/s101070100262 David P. Williamson The primal-dual method for approximation algorithms Received: June 19, 2000 / Accepted: February 7, 2001 Published online October 2, 2001 – Springer-Verlag 2001 Abstract. In this survey, we give an overview of a technique used to design and analyze algorithms that provide approximate solutions to NP-hard problems in combinatorial optimization. Because of parallels with the primal-dual method commonly used in combinatorial optimization, we call it the primal-dual method for approximation algorithms. We show how this technique can be used to derive approximation algorithms for a number of different problems, including network design problems, feedback vertex set problems, and facility location problems. 1. Introduction Many problems of interest in combinatorial optimization are considered unlikely to have efficient algorithms; most of these problems are NP-hard, and unless P = NP they do not have polynomial-time algorithms to find an optimal solution. Researchers in combinatorial optimization have considered several approaches to deal with NP-hard problems. These approaches fall into one of two classes. The first class contains al- gorithms that find the optimal solution but do not run in polynomial time. Integer programming is an example of such an approach. Integer programmers attempt to de- velop branch-and-bound (or branch-and-cut, etc.) algorithms for dealing with particular problems such that the algorithm runs quickly enough in practice for instances of in- terest, although the algorithm is not guaranteed to be efficient for all instances. The second class contains algorithms that run in polynomial time but do not find the optimal solution for all instances.
    [Show full text]
  • Combinatorial Approximation Algorithms
    Yuval Rabani, David Shmoys, Gerhard Woeginger (editors): Combinatorial Approximation Algorithms Dagstuhi-Seminar-Report; 187 18.08.-22.08.97 (9734) ISSN 0940-1 121 Copyright © 1998 IBFIgem. GmbH, Schloss Dagstuhl, D—66687Wadern, Germany Tel.: +49 —6871 - 9050, Fax: +49 - 6871 - 905 133 Das Internationale Begegnungs- und Forschungszentrum für Informatik(IBFI)ist eine gemeinnützige GmbH. Sie veranstaltet regelmäßig wissenschaftliche Seminare, welche nach Antrag der Tagungsleiter und Begutachtung durch das wissenschaftliche Direktoriummitpersönlich eingeladenen Gästen durchgeführt werden. Verantwortlich für das Programm ist das \Mssenschaftliche Direktorium: Prof. Dr. Thomas Beth, Prof. Dr. Peter Gorny, Prof. Dr. Thomas Lengauer, Prof. Dr. Klaus Madlener, Prof. Dr. Erhard Plödereder, Prof. Dr. Joachim W. Schmidt, Prof. Dr. Otto Spaniol, Prof. Dr. Christoph Walther, Prof. Dr. Reinhard Vthlhelm(wissenschaftlicher Direktor) Gesellschafter. Gesellschaft für Informatike.V., Bonn, TH Darmstadt. Universität Frankfurt, Universität Kaiserslautern, Universität Karlsruhe, Universität Stuttgart, Universität Trier, Universitätdes Saarlandes. Träger. Die Bundesländer Saarland und Rheinland—Pfalz Bezugsadresse: Geschäftsstelle Schloss Dagstuhl Universitätdes Saarlandes Postfach 15 11 50 D—66041Saarbrücken, Germany Tel.: +49 -681 - 302 4396 Fax: +49 -681 - 302 4397 e-mail: offi[email protected] url: http://www.dag.uni-sb.de Viele Seminarreporte sind auch in elektronischer Form über die Internetseite verfügbar. In einigen Fällen kann man dort von den Abstracts der Vorträge zu den Vollversionenweiterschalten. 2 APPROXIMATION ALGORITHMS Summary The Dagstuhl seminar on Combinatorial Approximation Algorithms brought together 54 researchers with affiliations in Austria (1), France (1), Germany (11), Hungary (1), Iceland (1), Israel (7), Italy (1), Netherlands (2), Russia (1), Sweden (1), Switzerland (1), United Kingdom (3), and USA (23). In 35 talks the participants presented their latest results on approximation algorithms, covering a wide range of topics.
    [Show full text]
  • Contents U U U
    Contents u u u ACM Awards Reception and Banquet, June 2018 .................................................. 2 Introduction ......................................................................................................................... 3 A.M. Turing Award .............................................................................................................. 4 ACM Prize in Computing ................................................................................................. 5 ACM Charles P. “Chuck” Thacker Breakthrough in Computing Award ............. 6 ACM – AAAI Allen Newell Award .................................................................................. 7 Software System Award ................................................................................................... 8 Grace Murray Hopper Award ......................................................................................... 9 Paris Kanellakis Theory and Practice Award ...........................................................10 Karl V. Karlstrom Outstanding Educator Award .....................................................11 Eugene L. Lawler Award for Humanitarian Contributions within Computer Science and Informatics ..........................................................12 Distinguished Service Award .......................................................................................13 ACM Athena Lecturer Award ........................................................................................14 Outstanding Contribution
    [Show full text]
  • Approximation Algorithms for Edge Partitioned Vertex Cover Problems
    Approximation Algorithms for the Partition Vertex Cover Problem Suman K. Bera1, Shalmoli Gupta2, Amit Kumar3, and Sambuddha Roy1 1 IBM India Research Lab, New Delhi 2 University of Illinois at Urbana-Champaign 3 Indian Institute of Technology Delhi Abstract. We consider a natural generalization of the Partial Vertex Cover problem. Here an instance consists of a graph G = (V,E), a cost + function c : V → Z , a partition P1,...,Pr of the edge set E, and a parameter ki for each partition Pi. The goal is to find a minimum cost set of vertices which cover at least ki edges from the partition Pi. We call this the Partition-VC problem. In this paper, we give matching upper and lower bound on the approximability of this problem. Our algorithm is based on a novel LP relaxation for this problem. This LP relaxation is obtained by adding knapsack cover inequalities to a natural LP relaxation of the problem. We show that this LP has integrality gap of O(log r), where r is the number of sets in the partition of the edge set. We also extend our result to more general settings. 1 Introduction The Vertex Cover problem is one of the most fundamental NP-hard problems and has been widely studied in the context of approximation algorithms [21, 10]. In this problem, we are given an undirected graph G = (V, E) and a cost function c : V → Z+. The goal is to find a minimum cost set of vertices which cover all the edges in E : a set of vertices S covers an edge e if S contains at least one of the end-points of e.
    [Show full text]
  • Approximation Algorithms for Demand-Robust Covering Problems∗
    How to Pay, Come What May: Approximation Algorithms for Demand-Robust Covering Problems∗ Kedar Dhamdhere† Vineet Goyal R. Ravi Google Inc. Tepper School of Business Tepper School of Business Mountain View, CA 94043 Carnegie Mellon University Carnegie Mellon University [email protected] Pittsburgh, PA 15213 Pittsburgh, PA 15213 [email protected] [email protected] Mohit Singh Tepper School of Business Carnegie Mellon University Pittsburgh, PA 15213 [email protected] Abstract 1 Introduction Robust optimization has traditionally focused on uncer- Robust optimization has roots in both Decision The- tainty in data and costs in optimization problems to formu- ory [14, 16] as well as Mathematical Programming [4]. late models whose solutions will be optimal in the worst- While minmax regret approaches were advanced in the for- case among the various uncertain scenarios in the model. mer field as conservative decision rules, robust optimization While these approaches may be thought of defining data- or was discussed along with stochastic programming [3] as an cost-robust problems, we formulate a new “demand-robust” alternate way of tackling data uncertainty. model motivated by recent work on two-stage stochastic op- More recent attempts at capturing the concept of ro- timization problems. We propose this in the framework of bust solutions in optimization problems include the work of general covering problems and prove a general structural Rosenblatt and Lee [21] in facility design problem, Mulvey lemma about special types of first-stage solutions for such et al. [17] in mathematical programming, and most recently, problems: there exists a first-stage solution that is a mini- Kouvelis and Yu [12] in combinatorial optimization; here mal feasible solution for the union of the demands for some robust means “good in all or most versions”, a version be- subset of the scenarios and its objective function value is ing a plausible set of values for the data in the model.
    [Show full text]