Almost Optimal Semi-streaming Maximization for k-Extendible Systems
Moran Feldman, Ran Haba

TL;DR
This paper introduces an almost optimal semi-streaming algorithm for maximizing weight under k-extendible constraints, significantly improving approximation ratios in the data stream model.
Contribution
The paper presents a semi-streaming O(k log k)-approximation algorithm for the general k-extendible maximization problem, nearly matching offline bounds.
Findings
Achieves near-optimal approximation ratio in semi-streaming model
Improves upon previous algorithms with higher approximation factors
Bridges gap between restricted and general cases
Abstract
In this paper we consider the problem of finding a maximum weight set subject to a -extendible constraint in the data stream model. The only non-trivial algorithm known for this problem to date---to the best of our knowledge---is a semi-streaming -approximation algorithm (Crouch and Stubbs, 2014), but semi-streaming -approximation algorithms are known for many restricted cases of this general problem. In this paper, we close most of this gap by presenting a semi-streaming -approximation algorithm for the general problem, which is almost the best possible even in the offline setting (Feldman et al., 2017).
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsComplexity and Algorithms in Graphs · Optimization and Search Problems · Stochastic Gradient Optimization Techniques
\forLoop
126calBbCounter
Almost Optimal Semi-streaming
Maximization for -Extendible Systems
Moran Feldman Dept. of Mathematics and Computer Science, Open University of Israel. E-mail: [email protected]
Ran Haba Dept. of Mathematics and Computer Science, Open University of Israel. E-mail: [email protected]
Abstract
In this paper we consider the problem of finding a maximum weight set subject to a -extendible constraint in the data stream model. The only non-trivial algorithm known for this problem to date—to the best of our knowledge—is a semi-streaming -approximation algorithm (Crouch and Stubbs, 2014), but semi-streaming -approximation algorithms are known for many restricted cases of this general problem. In this paper, we close most of this gap by presenting a semi-streaming -approximation algorithm for the general problem, which is almost the best possible even in the offline setting (Feldman et al., 2017).
Keywords: -extendible systems, streaming, combinatorial optimization, greedy algorithms
1 Introduction
Many problems in combinatorial optimization can be cast as special cases of the following general task. Given a ground set of weighted elements, find a maximum weight subset of obeying some constraint . In general, one cannot get any reasonable approximation ratio for this general task since it captures many hard problems such as maximum independent set in graphs. However, the existing literature includes many interesting classes of constraints for which the above task becomes more tractable. In particular, in the 1970’s Jenkyns [14] and Korte and Hausmann [15] suggested, independently, a class of constraints named -set system constraints which represents a sweet spot between generality and tractability. On the one hand, finding a maximum weight set subject to a -set system constraint captures many well known problems such as matching in hypergraphs, matroid intersection and asymmetric travelling salesperson. On the other hand, -set system constraints have enough structure to allow a simple greedy algorithm to find a maximum weight set subject to such a constraint up to an approximation ratio of .111 is a parameter of the constraint which intuitively captures its complexity. The exact definition of is given in Section 2, but we note here that in many cases of interest is quite small. For example, matroid intersection is a -set system.
The -approximation obtained by the greedy algorithm for finding a maximum weight set subject to a -set system constraint was recently shown to be the best possible [1]. Nevertheless, over the years many works improved over it either by achieving a better guarantee for more restricted classes of constraints [9, 17, 18], or by extending the guarantee to more general objectives (such as maximizing a submodular function) [7, 9, 10, 11, 16, 17, 20, 23]. Unfortunately, many of the above mentioned improvements are based on quite slow algorithms. Moreover, as modern applications require the processing of increasingly large amounts of data, even the simple greedy algorithm is often viewed these days as too slow for practical use. This state of affairs has motivated recent works aiming to study the problem of finding a maximum weight set subject to a -set system constraint in a Big Data oriented setting such as Map-Reduce and the data stream model. For the Map-Reduce setting, Ponte Barbosa et al. [6] essentially solved this problem by presenting a -approximation Map-Reduce algorithm for it using rounds, which almost matches the optimal approximation ratio in the sequential setting. In contrast, the situation for the data stream model is currently much more involved.
The only non-trivial data stream algorithm known to date (as far as we know) for finding a maximum weight set subject to a general -set system constraint is a -approximation semi-streaming algorithm by Crouch and Stubbs [5]. As one can observe, there is a large gap between the last approximation ratio and the -approximation that can be achieved in the offline setting. Several works partially addressed this gap by providing an -approximation semi-streaming algorithms for more restricted classes of constraints, the most general of which is known as -matchoid constraints [3, 4, 8, 21]. However, these results cannot be considered a satisfactory solution for the gap because -matchoid constraints are much less general than -set system constraints.222We do not formally define -matchoid constraints in this paper, but it should be noted that they usually fail to capture knapsack like constraints. For example, a single knapsack constraint in which the ratio between the largest and smallest item sizes is at most is a -set system constraint, but usually not a -matchoid constraint.
In this paper we make a large step towards resolving the above gap. Specifically, we present an -approximation semi-streaming algorithm for finding a maximum weight set subject to a class of constraints, known as -extendible constraints, that was introduced by [19] and captures (to the best of our knowledge) all the special cases of -set system constraints studied in the literature to date (including, in particular, -matchoid constraints). Formally, we prove the following theorem.
Theorem 1.1**.**
There is a polynomial time semi-streaming algorithm achieving -approximation for the problem of finding a maximum weight set subject to a -extendible constraint. Assuming it takes constant space to store a single element and a single weight, the space complexity of the algorithm is , where is the maximum size of a feasible set according to the constraint.
As the class of -extendible constraints captures every other restricted class of -set system constraints from the literature, we believe Theorem 1.1 represents the final intermediate step before closing the above mentioned gap completely (i.e., either finding an semi-streaming algorithm for -set system constraints, or proving that this cannot be done). It should also be mentioned that the approximation ratio guaranteed by Theorem 1.1 is optimal up to an factor since it is known that one cannot achieve better than -approximation for finding a maximum weight set subject to a -extendible constraint even in the offline setting [7].
1.1 Additional Related Work
In the -dimensional matching problem, one is given a weighted hypergraph in which the vertices are partitioned into subsets, and every edge contains exactly one vertex from each one of these subsets. The objective in this problem is to find a maximum weight matching in the hypergraph. Hazan et al. [12] showed that no algorithm can achieve a better than -approximation for -dimension matching unless . Interestingly, it turns out that -dimensional matching is captured by all the standard restricted cases of the the problem of finding a maximum weight set subject to -set system constraint, and thus, the inapproximability of Hazan et al. [12] extends to them as well. For most of these restricted cases this is the strongest inapproximability known, although a tight inapproximability of was proved for -set system and -extendible constraints by [1] and [7], respectively.
Complementing the hardness result of [12], some works presented algorithmic results for either -dimensional matching or natural generalizations of it such as -set packing [2, 13, 22].
2 Preliminaries and Notation
In this section we formally define some of the terms used in Section 1 and the notation that we use in the rest of this paper. Given a ground set , an independence system over this ground set is a pair in which is a non-empty collection of subsets of (formally, ) which is down-closed (i.e., if is a set in and is a subset of , then also belongs to ). One easy way to get an example of an independence system is to take an arbitrary vector space , designate the set of vectors in this space as the ground set , and make the collection of all independent sets of vectors in . Since removing a vector from an independent set of vectors cannot make the set dependent, the pair obtained from in this way is indeed an independence system.
The above example for getting an independence system from a vector space was one of the original motivations for the study of independence systems, and thus, a lot of the terminology used for independence systems is borrowed from the world of vector spaces. In particular, a set is called independent in a given independence system if and only if it belongs to , and it is called a base of the independence system if it is an inclusion-wise maximal independent set. Using this terminology, we can now define -set systems.
Definition 2.1**.**
An independence system is a -set system for an integer if for every set , all the bases of have the same size up to a factor of (in other words, the ratio between the sizes of the largest and smallest bases of is at most ).
An immediate consequence of the definition of -set systems is that any base of such a system is a maximum size independent set up to an approximation ratio of . Thus, one can get a -approximation for the problem of finding a maximum size independent set in a given -set system by outputting an arbitrary base of the -set system, which can be done using the following simple strategy, which we call the unweighted greedy algorithm. Start with the empty solution, and consider the elements of the ground set in an arbitrary order. When considering an element, add it to the current solution, unless this will make the solution dependent (i.e., not independent).
A -set system constraint is a constraint defined by a -set system, and a set obeys this constraint if and only if it is independent in that -set system. Note that using this notion we can refer to the problem studied in the previous paragraph as finding a maximum cardinality set subject to a -set system constraint. More generally, given a weight function and a -set system over the same ground set, it is often useful to consider the problem of finding a maximum weight set subject to the constraint corresponding to this -set system (the weight of a set is defined as ). Jenkyns [14] and Korte and Hausmann [15] showed that one can get a -approximation for this problem using an algorithm, known simply as the greedy algorithm, which is a variant of the unweighted greedy algorithm that considers the elements of in a non-decreasing weight order.
The definition of -set systems is very general, which occasionally does not allow them to capture all the necessary structure of a given application. Thus, various stronger kinds of independent set systems have been considered over the years, the most well known of which is the intersection of matroids (which is equivalent to a -set system for , and represents a strictly smaller class of independence systems for larger values of ). In this work we consider another kind of independence systems, which was originally defined by [19]. In this definition we use the expression to denote the union . We use the plus sign in a similar way throughout the rest of this paper.
Definition 2.2**.**
An independence system is a -extendible system for an integer if for any two sets and an element such that , there is a subset of size at most such that .
The class of -extendible systems is general enough to capture the intersection of matroids and every other restricted class of -set systems from the literature that we are aware of. In contrast, it is not difficult to verify that any -extendible system is a -set system. Thus, the greedy algorithm provides -approximation for the problem of finding a maximum weight set subject to a -extendible constraint—i.e., a constraint defined by a -extendible system and allowing only sets that are independent in this system.
In the data stream model version of the above problem, the elements of the ground set of a -extendible system arrive one after the other in an adversarially chosen order. An algorithm for this model views the elements of as they arrive, and it gets to know the weight of every element upon its arrival. Additionally, as is standard in the field, we assume the algorithm has access to an independence oracle that given a set answers whether is independent. The objective of the algorithm is to output a maximum weight independent set of the -extendible system. If the algorithm is allowed enough memory to store the entire input, then the data stream model version becomes equivalent to the offline version of the problem. Thus, an algorithm for this model is interesting only if it has a low space complexity. Since any algorithm for this model must use at least the space necessary for storing its output, most works on this model look for semi-streaming algorithms, which are data stream algorithms whose space complexity is upper bounded by —where is the maximum size of an independent set and is the size of the ground set. In particular, we note that the space complexity guaranteed by Theorem 1.1 falls within this regime because by definition, and one can assume that because any independence system is -extendible.
One can observe that the unweighted greedy algorithm (unlike the greedy algorithm itself) can be implemented as a semi-streaming algorithm because it considers the elements in an arbitrary order. This observation is crucial for our result since the algorithm we develop is heavily based on using the unweighted greedy algorithm as a subroutine (a similar use of the unweighted greedy algorithm is done by the current state-of-the-art algorithm for the problem due to Crouch and Stubbs [5]).
Paper Organization:
In Section 3 we present a reduction that allows us to assume that the weights of the elements are powers of , at the cost of losing a factor of in the space complexity of the algorithm. Using this reduction, we present a basic version of our algorithm in Section 4. This basic version presents our main new ideas, but achieves semi-streaming space complexity only under the simplifying assumption that the ratio between the maximum and minimum element weights is polynomially bounded. This simplifying assumption can be dropped using standard techniques, and we defer the details to Appendix A.
3 Reduction to -Power Weights
In this section we present a reduction that allows us to assume that the weights of all the elements in the ground set are powers of . This reduction simplifies the algorithms we present later in this paper. However, before presenting the reduction itself, let us note that we assume from this point on that for some integer . This assumption is without loss of generality because if does not obey it, then we can increase its value to the nearest integer that does obey it. Since the new value of is larger than the old value by at most a factor of , the approximation ratio guaranteed for both values of by Theorem 1.1 is asymptotically equal.
We say that an instance of the problem of finding a maximum weight set subject to a -extendible constraint is a -power instance if the weights of all the elements in it are powers of .
Reduction 3.1**.**
Assume that we are given a polynomial time data stream algorithm for the problem of finding a maximum weight set subject to a -extendible constraint. If provides -approximation for -power instances of the problem using space, then there exists a polynomial time data stream algorithm for the same problem which achieves -approximation for arbitrary instances using space. Moreover, if the weights of all the elements fall within some range , then it suffices for to provide -approximation for -power instances in which all the weights fall within the range .
Before presenting the algorithm we use to prove the above reduction, we need to define some additional notation. Let , and note that is a positive integer because we assume that is at least and a power of . For every element of weight , we define an auxiliary weight . Intuitively, is the highest power of which is not larger than . The following observation formally states the properties of that we need. In this observation we use the notation .
Observation 3.2**.**
For every element , is a power of obeying and .
Proof.
The first part of the observation, namely that is a power of , follows immediately from the definition of . Thus, we concentrate here on proving the other parts of the observation.
Note that
[TABLE]
Rearranging the last equality, we get
[TABLE]
and
[TABLE]
To complete the proof of the observation, we note that it also holds that
[TABLE]
We are now ready to present the algorithm that we use to prove Reduction 3.1, which appears as Algorithm 1. To intuitively understand this algorithm, it is useful to think of as the “class” element belongs to. All the elements within class have weights between and , and thus, treating them all as having the weight does not affect the approximation ratio by more than a factor of . Let us call the characteristic weight of class . Note now that the ratio between the characteristic weight of class and the characteristic weight of class is , which is a power of whenever is an integer multiple of . Thus, one can group the classes into groups such that the ratio between the characteristic weights of any pair of classes within a group is a power of (see Figure 1 for a graphial illustration of these groups). Moreover, by multiplying all the characteristic weights in the group by an appropriate scaling factor, one can make them all powers of . This means that for every group there exists a transformation that converts all the weights of the elements in it to powers of and preserves the ratio between any two weights in the group up to a factor of . In particular, we get that the elements of the group after the transformation form a -power instance.
Adding up all the above, we have described a way to transform any instance of finding a maximum weight independent set subject to a -extendible constraint into new instances of this problem that are guaranteed to be -power. Algorithm 1 essentially creates these new instances on the fly, and feeds them to copies of the algorithm whose existence is assumed in Reduction 3.1. Given this point of view, should be understood as the group to which element belongs, and is the transformed weight of . Observation 3.2 can now be interpreted as stating that the ratio between the weights of elements belonging to the same group (and thus, having the same value) is indeed changed by the transformation by at most a factor of .
In the rest of this section, we use to denote the set of elements fed to instance by Algorithm 1, and to denote the output of Algorithm 1. Additionally, we denote by an arbitrary (fixed) optimal solution for the original instance recieved by Algorithm 1. The following lemma proves that Algorithm 1 has the approximation ratio guaranteed by Reduction 3.1.
Lemma 3.3**.**
.
Proof.
Since Algorithm 1 feeds every arriving element into exactly one of the instances , the sets form a disjoint partition of . Thus,
[TABLE]
Hence, by an averaging argument, there must exist an index such that .
We now note that it follows from the pseudocode of Algorithm 1 and Observation 3.2 that the copies of get only weights that are powers of , and moreover, these weights belong to the range whenever the original weights received by Algorithm 1 belong to the range . Thus, by the assumption of Reduction 3.1, achieves -approximation for the instance it faces. Since is a feasible solution within this instance and is the output of , we get . Therefore,
[TABLE]
where the second and penultimate inequalities hold by Observation 3.2, and the last inequality is due to the fact that is the best solution among . ∎
The next lemma analyzes the space complexity of Algorithm 1 and completes the proof of Reduction 3.1.
Lemma 3.4**.**
Algorithm 1’s space complexity is .
Proof.
Algorithm 1 runs parallel copies of , each of them is assumed (by Reduction 3.1) to use space. Thus, the space required by these copies is . In addition to this space, Algorithm 1 only requires enough space to do two things.
- •
The algorithm has to store the outputs of the copies of . However, these outputs are originally stored by the copies themselves, and thus, storing them requires no more space than what is used by the copies.
- •
Calculate the sum of the weights of the elements in the solutions produced by the copies of . Since we assume that the weight of an element can be stored in constant space, this requires again (up to constant factors) no more space than the space used by the copies of to store their solutions. ∎
4 Algorithm
In this section we present a data stream algorithm for -power instances of the problem of finding a maximum weight set subject to a -extendible constraint. This algorithm assumes access to positive upper bound and lower bound on the weights of all the elements, and has a semi-streaming space complexity when the ratio between and is upper bounded by a polynomial in . Proposition 4.1 states the properties that we prove for this algorithm more formally.
Proposition 4.1**.**
There exists a -approximation data stream algorithm for -power instances of the problem of finding a maximum weight set subject to a -extendible constraint. This algorithm assumes access to positive upper bound and lower bound on the weights of all the elements, and its space complexity is under the assumption that constant space suffices to store a single element and a single weight.
Before getting to the proof of Proposition 4.1, we note that together with Reduction 3.1 this proposition immediately implies the following corollary.
Corollary 4.2**.**
There exists an -approximation data streaming algorithm for the problem of finding a maximum weight set subject to a -extendible constraint. This algorithm assumes access to positive upper bound and lower bound on the weights of all the elements. The space complexity of this algorithm is under the assumption that constant space suffices to store a single element and a single weight.
Note that when the ratio between and is polynomial in , the space complexity of the algorithm from Corollary 4.2 becomes , and thus, the algorithm is semi-streaming. In Appendix A we explain how the algorithm can be modified so that it keeps the “effective” ratio on the order of even when no values and are supplied to the algorithm and the weights of the elements come from an arbitrary range. This leads to the space complexity of stated in Theorem 1.1.
The rest of this section is devoted to the proof of Proposition 4.1. As a first step towards this goal, let us recall that the unweighted greedy algorithm is an algorithm that considers the elements of the ground set in an arbitrary order, and adds every considered element to the solution it constructs if that does not violate independence. As mentioned above, it follows immediately from the definition of -set systems that the unweighted greedy algorithm achieves an approximation ratio of for the problem of finding a maximum cardinality independent set subject to a -set system constraint. Since -set systems generalize -extendible systems, the same is true also for -extendible constraints. The following lemma improves over this by showing a tighter guarantee for -extendible constraints.
Lemma 4.3**.**
Given a -extendible set system , the unweighted greedy algorithm is guaranteed to produce an independent set such that for any independent set .
Proof.
Let us denote the elements of by in an arbitrary order. Using these elements, we recursively define a series of independent sets . The set is simply the set . For , we define using as follows. Since is a -extendible system and the subsets and are both independent, there must exist a subset such that and . Using the subset , we now define . Note that by the definition of , as promised. Furthermore, since for each , we know that , which implies because . However, , as the output of the unweighted greedy algorithm, must be inclusion-wise maximal independent set (i.e., a base), and thus, it must be in fact equal to the independent set containing it.
Let us now denote , and consider two different ways to bound the number of elements in . On the one hand, since every set includes up to elements, we get . On the other hand, the fact that implies that every element of belongs to for some value of , and therefore, . The lemma now follows by combining these two bounds. ∎
We are now ready to present the algorithm we use to prove Proposition 4.1, which is given as Algorithm 2. This algorithm has two main stages. In the first stage, the algorithm runs an independent copy of the unweighted greedy algorithm for every possible weight of elements. The copy corresponding to the weight is denoted by in the pseudocode of the algorithm, and Algorithm 2 feeds to it only the input elements whose weight is at least . The output of is denoted by in the algorithm. We also denote in the analysis by the set of elements fed to . By definition, is obtained by running the unweighted greedy algorithm on the elements of , which is a property we use below.
In the second stage of Algorithm 2 (which is done as a post-processing after the stream has ended), the algorithm constructs an output set based on the outputs of the copies of the unweighted greedy algorithm. Specifically, this is done by running the unweighted greedy algorithm on the elements of , considering the elements of the sets in a decreasing value of order. While doing so, the given pseudocode also keeps in the temporary solution obtained by the unweighted greedy algorithm after considering only the elements of for . This temporary solution is used by the analysis below, but need not be kept by a real implementation of Algorithm 2.
We begin the analysis of Algorithm 2 by analyzing its space complexity.
Lemma 4.4**.**
Algorithm 2 can be implemented using a space complexity of .
Proof.
Note that each copy of the unweighted greedy algorithm only has to store its solution, which contains up to elements since it is independent. Algorithm 2 uses such copies, and thus, the space it needs for these copies is only
[TABLE]
In addition to the space used by the copies of the unweighted greedy algorithm, Algorithm 2 only needs to store the set . This set contains a subset of the elements from the outputs of the above copies, and thus, can increases the space required only by a constant factor. ∎
To complete the proof of Proposition 4.1, it remains to analyze the approximation ratio of Algorithm 2. We begin with the following lemma, which is the technical heart of our analysis. Like in Section 3, let us denote by be an arbitrary (fixed) optimal solution to the problem we want to solve. We also assume for consistency that (note that is not defined by Algorithm 2).
Lemma 4.5**.**
For each integer , .
Proof.
The set can be viewed as the output of the unweighted greedy algorithm running on . Since we also know that is independent, Lemma 4.3 guarantees
[TABLE]
Adding to both its sides, we get
[TABLE]
where the last inequality holds since the unweighted greedy algorithm achieves -approximation and is an independent set within (recall that is the set of elements that were fed to ). Using the last inequality we can now get
[TABLE]
where the first equality holds because , and the second inequality holds because (recall that the algorithm constructs by adding elements of to ). The lemma now follows by rearranging the above inequality and multiplying it by . ∎
Using the last lemma, we can prove the existence of a useful mapping from the elements of to the elements of .
Lemma 4.6**.**
There exists a mapping such that
for each , . 2. 2.
for each , . 3. 3.
for each , .
Proof.
We construct by scanning the elements and defining the mapping for every element scanned. To describe the order in which we scan the elements of , let us define . Note that is a disjoint partition of , and thus, any scan of the elements of is a scan of the elements of . Specifically, we scan the elements of by first scanning the elements of in an arbitrary order, then scanning the elements of in an arbitrary order, and so on. Consider now the situation when our scan gets to an arbitrary element of set . One can note that prior to scanning , we scanned (and mapped) only elements of , and thus, we mapped at most elements (the is due to the fact that , and was not mapped yet). Combining this with Lemma 4.5, we get that at the point in which we scan there must still be either an element that still has less than elements mapped to it or an element that still has less than elements mapped to it. We choose the mapping of to be an arbitrary such element .
Property 1 of the lemma is clearly satisfied by the above construction because we never map an element to an element that already has elements mapped to it. To see why Property 3 of the lemma also holds, note that every element must have a weight of by the definition of . This element is mapped by to some element , and the weight of is at least by the definition of . It remains to prove Property 2 of the lemma. Consider an arbitrary element of weight . The elements of whose weight is are exactly the elements of , and thus, we need to show that . Since all the elements of have weights of at least , cannot belong to . Thus, an element of can be mapped to when scanned only if has less than elements already mapped to it (if ) or not at all (if ), which implies that no more than elements of can get mapped to , which is exactly what we wanted to prove. ∎
We are now ready to prove the approximation ratio of Algorithm 2 (and complete the proof of Proposition 4.1).
Lemma 4.7**.**
Algorithm 2 is a -approximation algorithm for -power instances of the problem of finding a maximum weight set subject to a -extendible constraint.
Proof.
Let be the function whose existence is guaranteed by Lemma 4.6. The properties of this function imply that, for each element ,
[TABLE]
Thus,
[TABLE]
which completes the proof of the lemma. ∎
5 Conclusion
In this work we have presented the first semi-streaming -approximation algorithm for the problem of finding a maximum weight set subject to a -extendible constraint. This result is intrinsically interesting because the generality of -extendible constraints makes our algorithm applicable to many problems of interest. Additionally, we believe (as discussed in Section 1) that our result is likely to be the final intermediate step towards the goal of designing an algorithm with similar properties for general -set system constraints or proving that this cannot be done.
Given our work, the immediate open question is to settle the approximation ratio that can be obtained for -set system constraints in the data stream model. Another interesting research direction is to find out whether one can improve over the approximation ratio of our algorithm. Specifically, we leave open the question of whether there is a semi-streaming algorithm for finding a maximum weight set subject to a -extendible constraint whose approximation ratio is clean .
Appendix A Algorithm for General Weights
In this section we present a semi-streaming algorithm for -power instances of the problem of finding a maximum weight set subject to a -extendible constraint. Unlike Algorithm 2, this algorithm does not assume access to the bounds and , and its space complexity remains nearly linear regardless of the ratio between these bounds. A more formal statement of the properties of this algorithm is given in Proposition 4.1. Note that, together with Reduction 3.1, this proposition immediately implies Theorem 1.1.
Proposition A.1**.**
There exists a -approximation semi-streaming algorithm for -power instances of the problem of finding a maximum weight set subject to a -extendible constraint. The space complexity of this algorithm is under the assumption that constant space suffices to store a single element and a single weight.
Throughout this section we assume for simplicity that the -extendible system does not include any self-loops (a self-loop is an element such that is a dependent set—i.e., ). This assumption is without loss of generality since a self-loop cannot belong to any independent set, and thus, an algorithm can safely ignore self-loops if they happen to exist. One consequence of this assumption is that , where is an arbitrary fixed optimal solution like in the previous sections. This inequality holds since is a feasible solution for every element , and therefore, its weight cannot exceed the weight of .
As mentioned in Section 4, the algorithm we use to prove Proposition A.1 is a variant of Algorithm 2 that includes additional logic designed to force the ratio to be effectively polynomial—specifically, . Given access to and , this could be done simply by settings and and discarding any element whose weight is lower then .333Starting from this point, and are no longer necessarily upper and lower bounds on the weights of all the elements. However, they remain upper and lower bounds on the weights of the non-discarded elements. This guarantees that the ratio is small, and affects the weight of the optimal solution by at most a constant factor since the total weight of the elements of this solution that get discarded is upper bounded by
[TABLE]
Unfortunately, our algorithm does not have access (from the beginning) to and . As an alternative, this algorithm, which is given as Algorithm 3, does two things. First, it keeps equal to the maximum weight of the elements seen so far, which guarantees that eventually becomes . Second, it runs the unweighted greedy algorithm on the input it receives. The size of the solution maintained by the unweighted greedy algorithm, which we denoted by , provides an estimate for the maximum size of an independent set consisting only of elements that have already arrived. In particular, after all the elements arrive, because the unweighted greedy algorithm is a -approximation algorithm.
Given the above discussion and the fact that the final value of is an upper bound on , it is natural to define as and discard every element whose weight is lower than . Unfortunately, this does not work since and change during the execution of Algorithm 3, and reach their final values only when it terminates. Thus, we need to set to a more conservative (lower) value. In particular, Algorithm 3 uses .
Like Algorithm 2, Algorithm 3 maintains an instance of the unweighted greedy algorithm for every possible weight between and . However, doing so is somewhat more involved for Algorithm 3 because and change during the algorithm’s execution, which requires the algorithm to occasionally create and remove instances of unweighted greedy. The creation of such instances involves one subtle issue that needs to be kept in mind. In Algorithm 2 every instance of unweighted greedy associated with a weight receives all elements whose weight is at least . To mimic this behavior, when Algorithm 3 creates new instances of unweighted greedy following a decrease in (which can happen when increases), the newly created instances are not fresh new instances but copies of the instance of unweighted greedy that was previously associated with the lowest weight.
The rest of the details of Algorithm 3 are identical to the details of Algorithm 2. Specifically, every arriving element is feed to every instance of unweighted greedy associated with a weight of or less, and at termination the outputs of all the unweighted greedy instances are combined in the same way in which this is done in Algorithm 2.
We now get to the analysis of Algorithm 3, and let us begin by bounding its space complexity. Let , , , and denote the values of , , , and , respectively, at the end of iteration number of Algorithm 3.
Lemma A.2**.**
Algorithm 3 can be implemented using a space complexity of .
Proof.
Using the same argument used in the proof of Lemma 4.4, it can be shown that the space complexity of Algorithm 3 is upper bounded by times the maximum number of unweighted greedy instances maintained by the algorithm at the same time. By making the deletions of unweighted greedy instances precede the creation of new instances within every given iteration of the main loop of Algorithm 3 (and avoiding the creation of instances that need to be immediately deleted), it can be guaranteed that the maximum number of instances of unweighted greedy maintained by Algorithm 3 at any given time is exactly . Thus, the algorithm’s space complexity is at most
[TABLE]
where the second inequality is due to the fact that is always the size of an independent set, and thus, cannot exceed . ∎
Our next objective is to analyze the approximation ratio of Algorithm 3. Like in the toy analysis presented above for the case in which the algorithm has access to and , the analysis we present starts by upper bounding the total weight of the discarded elements. However, to do that we need the following technical observation, which can be proved by induction.
Observation A.3**.**
Algorithm 3 maintains the invariant that, at the end of every one of its loops, if an element was fed to some instance of unweighted greedy currently kept by the algorithm, then it was fed exactly to those instances associated with a weight of at most .
We say that an element is discarded by Algorithm 3 if was never fed to the final instance (during the execution of Algorithm 3 there might be multiple instances of unweighted greedy named for —by final instance we mean the last of these instances). Let be the set of discarded elements.
Lemma A.4**.**
.
Proof.
For every , let be the -th element of to arrive, and let be its location in the input stream. Given Observation A.3, the fact that implies that was not feed to the final instance , which can only happen if an instance named either did not exist when arrived or was deleted at some point after ’s arrival. Thus, for some .
The crucial observation now is that because by the time arrives there are already elements of that arrived, and these elements form together an independent set of size (recall that is a -approximation for the maximum size of an independent set consisting only of elements that already arrived). Thus, we get
[TABLE]
where the first inequality holds since and both and are integers. Adding up the last inequality over yields
[TABLE]
The next lemma shows that Algorithm 3 has a good approximation ratio with respect to the non-discarded elements of .
Lemma A.5**.**
.
Proof.
Observe that is a -extendible system, derived from by removing all elements of . In addition, all the weights of the elements of this set system are powers of , and thus, by Proposition 4.1, Algorithm 2 achieves -approximation for the problem of finding a maximum weight independent set of . In other words, when Algorithm 2 is fed only the elements of , its output set obeys , where is an arbitrary maximum weight set independent set of .
We now note that one consequence of Observation A.3 is that, by the time Algorithm 3 terminates, the instances it maintains receive exactly the input received by the corresponding instances in Algorithm 2 when the last algorithm gets only the elements of as input. Since Algorithms 2 and 3 compute their outputs based on the outputs of in the same way, this implies that the output set of Algorithm 3 is identical to the output set produced by Algorithm 2 when this algorithm is given only the elements of as input.
Combining the above observations, we get
[TABLE]
where the last inequality holds since is a maximum weight independent set in and is independent in this set system. The lemma now follows by rearranging the last inequality. ∎
Corollary A.6**.**
, and thus, the approximation ratio of Algorithm 3 is at most .
Proof.
Combining the last two lemmata, one gets
[TABLE]
The corollary now follows by rearranging the above inequality. ∎
We conclude the section by noticing that Proposition A.1 is an immediate consequence of Lemma A.2 and Corollary A.6.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Ashwinkumar Badanidiyuru and Jan Vondrák. Fast algorithms for maximizing submodular functions. In SODA , pages 1497–1514, 2014.
- 2[2] Piotr Berman. A d / 2 𝑑 2 d/2 approximation for maximum weight independent set in d 𝑑 d -claw free graphs. Nord. J. Comput. , 7(3):178–184, 2000.
- 3[3] Amit Chakrabarti and Sagar Kale. Submodular maximization meets streaming: matchings, matroids, and more. Math. Program. , 154(1-2):225–247, 2015.
- 4[4] Chandra Chekuri, Shalmoli Gupta, and Kent Quanrud. Streaming algorithms for submodular function maximization. In ICALP , pages 318–330, 2015.
- 5[5] Michael Crouch and Daniel M. Stubbs. Improved streaming algorithms for weighted matching, via unweighted matching. In APPROX , pages 96–104, 2014.
- 6[6] Rafael da Ponte Barbosa, Alina Ene, Huy L. Nguyen, and Justin Ward. A new framework for distributed submodular maximization. In FOCS , pages 645–654, 2016.
- 7[7] Moran Feldman, Christopher Harshaw, and Amin Karbasi. Greed is good: Near-optimal submodular maximization via greedy optimization. In COLT , pages 758–784, 2017.
- 8[8] Moran Feldman, Amin Karbasi, and Ehsan Kazemi. Do less, get more: Streaming submodular maximization with subsampling. In Neur IPS 2018 , pages 730–740, 2018.
