Online Disjoint Spanning Trees and Polymatroid Bases
Abstract
Finding the maximum number of disjoint spanning trees in a given graph is a well-studied problem with several applications and connections. The Tutte-Nash-Williams theorem provides a min-max relation for this problem which also extends to disjoint bases in a matroid and leads to efficient algorithms [13]. Several other packing problems such as element disjoint Steiner trees, disjoint set covers, and disjoint dominating sets are NP-Hard but admit an -approximation [7, 4]. Călinescu, Chekuri, and Vondrák [2] viewed all these packing problems as packing bases of a polymatroid and provided a unified perspective. Motivated by applications in wireless networks, recent works have studied the problem of packing set covers in the online model [11, 6, 1]. The online model poses new challenges for packing problems. In particular, it is not clear how to pack a maximum number of disjoint spanning trees in a graph when edges arrive online. Motivated by these applications and theoretical considerations, we formulate an online model for packing bases of a polymatroid, and describe a randomized algorithm with a polylogarithmic competitive ratio. Our algorithm is based on interesting connections to the notion of quotients of a polymatroid that has recently seen applications in polymatroid sparsification [12]. We generalize the previously known result for the online disjoint set cover problem [6] and also address several other packing problems in a unified fashion. For the special case of packing disjoint spanning trees in a graph (or a hypergraph) whose edges arrive online, we provide an alternative to our general algorithm that is simpler and faster while achieving the same poly-logarithmic competitive ratio.
Keywords and phrases:
Disjoint Spanning Trees, Base Packing, Polymatroids, Online AlgorithmsCategory:
Track A: Algorithms, Complexity and GamesCopyright and License:
![[Uncaptioned image]](x1.png)
2012 ACM Subject Classification:
Theory of computation Approximation algorithms analysis ; Theory of computation Online algorithmsFunding:
This work was supported in part by NSF grant CCF-2402667.Editors:
Keren Censor-Hillel, Fabrizio Grandoni, Joël Ouaknine, and Gabriele PuppisSeries and Publisher:

1 Introduction
Finding the maximum number of disjoint spanning trees in a given graph is a well-studied problem with several applications and connections. The Tutte-Nash-Williams theorem provides a min-max relation for this problem which also extends to the maximum number of disjoint bases in a matroid with efficient algorithms [13]. Several other packing problems such as disjoint set covers, disjoint dominating sets, and element disjoint Steiner trees are NP-Hard but admit an -approximation111We are using the convention of -approximation for a maximization problem with with the understanding that the returned value is at least . [7, 4]. Călinescu, Chekuri, and Vondrák [2] viewed these problems as packing bases of a polymatroid and provided a unified perspective. Motivated by applications in wireless networks, recent works have studied the problem of packing set covers in the online model [11, 6, 1]. The online model poses new challenges for packing problems. In particular, it is not clear how to pack a maximum number of disjoint spanning trees in a graph when edges arrive online. Motivated by these applications and theoretical considerations we consider the problem of packing disjoint bases of a polymatroid in the online model.
A polymatroid on a ground set is an integer-valued monotone submodular function with . We recall that a function is monotone if for every and submodular if for every . A subset is a base of if .222In most settings, one would define a set to be a base if it is inclusionwise minimal subject to satisfying . This is particularly important in the setting of matroids where all bases have the same cardinality which is not necessarily true in the polymatroidal setting. However, since we are interested in (approximating) the maximum number of disjoint bases we adopt the relaxed definition for simplicity. In the Disj-Polymatroid-Bases problem, we are given a polymatroid via an evaluation oracle. The goal is to find a maximum number of disjoint bases. We define
Polymatroids generalize matroid rank function, coverage functions, and many others. Consequently, numerous set packing problems can be cast as special cases of Disj-Polymatroid-Bases. We will discuss some of these special cases shortly.
Online Disj-Polymatroid-Bases Model.
We formally describe the online model for Disj-Polymatroid-Bases (denoted Online-Disj-Polymatroid-Bases). We have an underlying polymatroid over a large but finite ground set . Let . We index the elements of the ground set as where is the element that arrives at time for every . For each , we denote and let be the function obtained from by restricting the ground set to , that is, for every . At each timestep , the online algorithm has access to the evaluation oracle333The evaluation oracle of a function takes a set as input and returns . of the function restricted to the set of elements that have arrived until time , i.e., the evaluation oracle of the function , and has to color element irrevocably. A color is said to be a base color if the set of elements with that color is a base of . The goal of the online algorithm is to maximize the number of base colors. We remark that we are implicitly assuming that the elements of is the input sequence. The competitive ratio of an online algorithm is the ratio between the base colors in an optimal offline algorithm and that of the online algorithm. For randomized online algorithms, we will be interested in the expected competitive ratio. We assume that the online algorithm has prior knowledge of the function value of the ground set , i.e., . This assumption is motivated by applications to be discussed below.
Applications.
Polymatroids generalize coverage functions and matroid rank functions. We discuss these two special cases, the associated packing problems, and their online model below.
-
1.
In the Disj-Set-Cover problem, the input is a set system over a finite universe. We will alternatively view the set system as a hypergraph where corresponds to the universe and the hyperedges in correspond to sets in the system. The goal is to find a maximum number of disjoint set covers (a subset of hyperedges is a set cover if ). Disj-Set-Cover can be cast as a special case of Disj-Polymatroid-Bases by considering the coverage function of the hypergraph as the polymatroid, i.e., by considering the polymatroid defined as . Coverage functions are fairly general with several prominent special cases–e.g., the domatic number problem is a special case of Disj-Set-Cover [7]. In the online setting of Disj-Set-Cover (termed Online-Disj-Set-Cover), the vertex set is known in advance while the hyperedges are revealed in an online fashion. The online algorithm has to color each hyperedge immediately upon arrival irrevocably in order to maximize the number of colors that form a set cover.
-
2.
In the Disj-Matroid-Bases problem, we are given evaluation access to a matroid rank function over a ground set (we recall that a matroid rank function is a polymatroid that additionally satisfies for every ). A subset is a base of the matroid if . The goal is to find a maximum number disjoint bases of the matroid. Matroid rank function is a polymatroid and hence, Disj-Matroid-Bases is a special case of Disj-Polymatroid-Bases. In the online setting of Disj-Matroid-Bases (termed Online-Disj-Matroid-Bases), the ground set is revealed in an online fashion while the online algorithm has access to the rank function restricted to the set of elements that have arrived so far. The online algorithm has to color each element immediately upon arrival irrevocably in order to maximize the number of base colors.
Next, we describe three special cases of Disj-Polymatroid-Bases, namely Disj-Matrix-Bases, Disj-Spanning-Trees, and Disj-Conn-Spanning-Subhypergraphs. To the best of authors’ knowledge, these three problems have not been explored in the online setting. One of the motivations of this work is to understand these three problems in the online setting.
-
1.
In the Disj-Matrix-Bases problem, we are given a matrix of rank and the goal is to find a maximum number of disjoint spanning subsets of row-vectors – a subset of row vectors is spanning if its linear hull is . This is a special case of Disj-Matroid-Bases where the matroid is the linear matroid defined by (and consequently, the rank function of a subset of row-vectors is the dimension of the subspace spanned by them). In the online setting of Disj-Matrix-Bases (termed Online-Disj-Matrix-Bases), the rows of the matrix are revealed in an online fashion and the online algorithm has to color each row immediately upon arrival irrevocably in order to maximize the number of spanning colors.
-
2.
In the Disj-Spanning-Trees problem, we are given a connected graph and the goal is to find a maximum number of disjoint spanning trees in . Disj-Spanning-Trees is a special case of Disj-Matroid-Bases where the matroid is the graphic matroid defined by (and consequently, the rank function of a subset is , where is the number of components in the graph ). In the online setting of Disj-Spanning-Trees (termed Online-Disj-Spanning-Trees), the vertex set of is known in advance while the edges of are revealed in an online fashion and the online algorithm has to color each edge immediately upon arrival irrevocably in order to maximize the number of connected colors – a color is connected if the edges of the color form a connected graph over the vertex set .
-
3.
In the Disj-Conn-Spanning-Subhypergraphs problem, we are given a connected hypergraph and the goal is to find a maximum number of disjoint connected spanning subhypergraphs in . Disj-Conn-Spanning-Subhypergraphs is a special case of Disj-Polymatroid-Bases where the polymatroid of interest is defined as for every , where is the number of components in the hypergraph . Disj-Conn-Spanning-Subhypergraphs arises in the context of packing element-disjoint Steiner trees [4]. In the online setting of Disj-Conn-Spanning-Subhypergraphs (termed Online-Disj-Conn-Spanning-Subhypergraphs), the vertex set of is known in advance while the hyperedges of are revealed in an online fashion and the online algorithm has to color each hyperedge immediately upon arrival irrevocably in order to maximize the number of connected colors – a color is connected if the hyperedges of the color form a connected hypergraph over the vertex set .
Disj-Conn-Spanning-Subhypergraphs generalizes Disj-Spanning-Trees (in both offline and online settings): if the input hypergraph is a graph, then the problem corresponds to Disj-Spanning-Trees. Disj-Conn-Spanning-Subhypergraphs is also closely related to Disj-Set-Cover in the following sense: both problems are defined over hypergraphs – the latter asks for disjoint spanning subhypergraphs (which is equivalent to disjoint set covers) while the former asks for disjoint connected spanning subhypergraphs. Disj-Conn-Spanning-Subhypergraphs generalizes Disj-Set-Cover (in both offline and online settings) via the following approximation-preserving reduction from the latter to the former: given an instance of Disj-Set-Cover, add a new vertex and stretch every hyperedge of to include to obtain a new hypergraph ; the maximum number of disjoint set covers in is equal to the maximum number of disjoint connected spanning subhypergraphs in .
Prior work in the offline setting.
Disj-Matroid-Bases is polynomial-time solvable [5]. Disj-Set-Cover and Disj-Conn-Spanning-Subhypergraphs are -inapproximable and -approximable [7, 4]. Călinescu, Chekuri, and Vondrák [2] introduced Disj-Polymatroid-Bases as a unified generalization of Disj-Matroid-Bases, Disj-Set-Cover, and Disj-Conn-Spanning-Subhypergraphs. They designed an -approximation for Disj-Polymatroid-Bases by showing an approximate min-max relation for . Their approximate min-max relation is based on the following minimization problem:
It is easy to see that (e.g., see [2]). Moreover, if the polymatroid is a matroid rank function, then Edmonds [5] showed that . Edmonds’ result is constructive and implies a polynomial time algorithm for Disj-Matroid-Bases. However, for coverage functions, it is known that [7]. Călinescu, Chekuri, and Vondrák showed that this bound is tight by giving a polynomial-time algorithm to construct disjoint bases (and hence, ). Their algorithm unifies the approximation algorithms for Disj-Set-Cover [7] and Disj-Conn-Spanning-Subhypergraphs [4]. We state their algorithm since it will be important for the rest of our discussion: The algorithm computes the parameter (by guessing/binary search) and colors each element with a uniformly random color chosen from a color palette of size . Călinescu, Chekuri, and Vondrák showed that the expected number of base colors returned by this random coloring algorithm is at least . An alternative algorithm based on a random permutation is also described in [2] that we discuss later.
Prior work in the online setting.
Online-Disj-Set-Cover was introduced and studied by Pananjady, Bagaria, and Vaze [11] driven by applications to sensor networks, supply chain management, and crowd-sourcing platforms. In the context of Disj-Set-Cover, the quantity associated with the coverage function of the hypergraph has a simple interpretation: it is equal to the min-degree of the hypergraph . Pananjady, Bagaria, and Vaze showed that if the min-degree is known to the online algorithm in advance, then there is an online deterministic algorithm with competitive ratio , where is the number of vertices in the input hypergraph; we note that a randomized online algorithm with the same competitive ratio is an easy consequence of the random coloring algorithm of [2] discussed above. On the lower bound side, they showed that if min-degree is not known in advance, then every online deterministic algorithm for Disj-Set-Cover has competitive ratio . Although this lower bound result seems to suggest that knowing the min-degree of the graph in advance is required to achieve any meaningful competitive ratio, two different results have overcome this seeming technical barrier by empowering the online algorithm in other ways: Firstly, Emek, Goldbraikh, and Kantor [6] designed an online randomized algorithm with expected competitive ratio (assuming no knowledge of the min-degree but using randomness). On the lower bound side, they showed that every online randomized algorithm has expected impure competitive ratio444An online (randomized) algorithm for Disj-Set-Cover has impure competitive ratio if the (expected) number of set covers in the online algorithm is at least for some that is a function only of , where is the input hypergraph, is the number of vertices in , and is the maximum number of disjoint set covers in . Our work focuses on the case of , i.e., pure competitive ratio. (even with knowledge of min-degree). Secondly, Bienkowski, Byrka, and Jeż [1] designed an online deterministic algorithm with impure competitive ratio (assuming no knowledge of the min-degree but settling for impure competitive ratio).
For the more general problem of Online-Disj-Polymatroid-Bases, Pananjady, Bagaria, and Vaze [11] observed that if is known in advance, then it is possible to design a randomized online algorithm with expected competitive ratio : indeed, the random coloring algorithm of [2] mentioned above can be implemented in the online setting using knowledge of and it will have the stated competitive ratio (via the results of [2]). In this work, we are interested in Online-Disj-Polymatroid-Bases in the setting where is not known in advance.
Our motivations to consider Online-Disj-Polymatroid-Bases are multifold: Disj-Matrix-Bases and Disj-Spanning-Trees are fundamental problems in linear algebra and graph algorithms respectively. We note that Disj-Matrix-Bases in the online arrival model is non-trivial even for -dimensional vectors, i.e., for . Disj-Conn-Spanning-Subhypergraphs generalizes Disj-Spanning-Trees and arises in the context of packing element-disjoint Steiner trees [4]. As mentioned earlier, Disj-Conn-Spanning-Subhypergraphs can also be viewed as a generalization of Disj-Set-Cover. Although Disj-Set-Cover has been studied in the online model, there has not been any work on Disj-Conn-Spanning-Subhypergraphs in the online model.
1.1 Our Results
Throughout this work, we will denote the color palette by the set of natural numbers. There is a natural greedy algorithm for Online-Disj-Polymatroid-Bases: initialize color ; at each timestep : use color for element and if the set of elements with color is a base, then increment . For uniform random arrival order, the competitive ratio of this simple online algorithm is (via the results of [2]); this does not seem to have been explicitly noted in prior work. For arbitrary arrival order, the competitive ratio of this online algorithm is : the algorithm will return at least one base while the maximum number of possible bases is at most (since ). It is known that is the best possible competitive ratio of deterministic online algorithms that do not have prior knowledge of [11]. In this work, we are interested in the setting of arbitrary arrival order without prior knowledge of . For this setting, we design a randomized online algorithm with expected competitive ratio .
Theorem 1.
For Disj-Polymatroid-Bases, there exists a randomized online algorithm with expected competitive ratio . The runtime of the algorithm at each timestep is .
We recall that the best-known approximation factor for Disj-Polymatroid-Bases in the offline setting is , and hence the competitive ratio of our online algorithm nearly matches that of the best possible offline algorithm. We discuss the consequences of Theorem 1 for the applications. Specializing Theorem 1 to coverage functions implies a randomized online algorithm for Disj-Set-Cover with expected competitive ratio , where is the number of vertices of the input hypergraph, thus recovering the result of [6]. Specializing Theorem 1 to matroid rank functions implies a randomized online algorithm for Disj-Matroid-Bases with expected competitive ratio , where is the rank of the ground set. Consequently, we obtain a randomized online algorithm with expected competitive ratio for Disj-Matrix-Bases where is the dimension of the span of the input vectors and with expected competitive ratio for Disj-Spanning-Trees where is the number of vertices of the input graph. Specializing Theorem 1 to the polymatroid function that arises in Disj-Conn-Spanning-Subhypergraphs implies a randomized online algorithm for Disj-Conn-Spanning-Subhypergraphs with expected competitive ratio , where is the number of vertices of the input hypergraph.
Our randomized online algorithm to prove Theorem 1 is based on the notion of quotients of a polymatroid. Quotients played a central role in the recent result on polymatroidal sparsification [12]. The competitive ratio analysis of our algorithm is based on novel properties of quotients which might be of independent interest. We prove Theorem 1 in two steps: firstly, we design a randomized online algorithm with competitive ratio but as stated it needs to solve an NP-Hard problem. Next, we modify this algorithm to achieve polynomial run-time while achieving the same competitive ratio. For this, we rely on the strength decomposition of polymatroids which is computable in polynomial time and show a property connecting the strength decomposition to min-sized quotients. Although these algorithms are general and powerful, they are computationally expensive and difficult to interpret for specific applications. As our second result, we give a very simple and fast online randomized algorithm for Disj-Conn-Spanning-Subhypergraphs and Disj-Spanning-Trees that achieves the same competitive ratio.
Theorem 2.
For Disj-Conn-Spanning-Subhypergraphs, there exists a randomized online algorithm with competitive ratio , where is number of vertices in the input hypergraph. The runtime of the algorithm at each timestep is , where is the hyperedge that arrives at time . In particular, the algorithm for Disj-Spanning-Trees can be implemented to run in constant time at each timestep.
1.2 Technical Background and Overview of Algorithms
We recall that we assume is the input sequence. Let be the polymatroid of interest. Let be the function value of bases and . For the discussion here, we will assume that for every , , and since this is the non-trivial case. First consider the setting where the online algorithm knows in advance. In this setting, coloring each element uniformly from the palette gives base colors in expectation and, consequently achieves -competitive ratio.
Suppose is not known in advance. Our approach is to estimate at each timestep. At each timestep , suppose that we can find a number that is within a factor of , i.e., – we call such a as a coarse estimate. Let be a uniform random sample from ; is a -approximation of with probability . Consequently, coloring the element uniformly from the palette leads to base colors in expectation.
The key challenge is computing a coarse estimate for which depends on the full input sequence. This may not be feasible at all timesteps owing to the limited knowledge of the polymatroid. Instead, we compute an estimate at each timestep such that the elements at timesteps which achieve the coarse property provide a sufficiently large value of . We formalize this approach now. Suppose we have a subroutine to compute some estimate at each timestep . We call a timestep to be good if the estimate is a coarse estimate of and bad otherwise. We let be the collection of elements that arrive at good timesteps and focus on the function restricted to , i.e., on the function . Suppose that the following two properties hold: (i) Bases of are also bases of and (ii) . These two properties suffice to obtain base colors in expectation via the fine-tuned random coloring argument in the previous paragraph.
We rely on the notion of quotients to obtain an estimator satisfying properties (i) and (ii).
Definition 3.
For a polymatroid , a set is a quotient of if , i.e., if each element has strictly positive marginal with respect to .
See Appendix A for the interpretation of quotients for some concrete polymatroids. We observed that the minimum size of a non-empty quotient of is an -approximation of – see Appendix B. This inspired us to use the minimum-sized quotient of that contains as the estimate at each timestep . That is, we use the estimator defined as
(1) |
Our main technical contribution is showing that the preceding estimator satisfies the two properties that we discussed earlier. This result combined with the random coloring process yields the desired -competitive ratio. A technical issue is that the estimator is NP-Hard to compute for general polymatroids (even for matroid rank functions). However, we note that a -approximation of is sufficient for the desired competitive ratio. We show that an -approximation of can be computed using a strength decomposition with respect to which admits a polynomial-time algorithm.
We briefly discuss the idea behind the simpler algorithm for Disj-Conn-Spanning-Subhypergraphs and Disj-Spanning-Trees of Theorem 2. We will focus on Disj-Spanning-Trees here. The polymatroid of interest here is the rank function of the graphic matroid defined by the undirected graph . For this case, the quotient computation is polynomial-time solvable since is equal to the min - cut value in the graph where and . However, we show that a much simpler estimator for Disj-Spanning-Trees also achieves properties (i) and (ii) mentioned above: we use
This estimator generalizes also to Disj-Conn-Spanning-Subhypergraphs.
Organization.
In Section 2, we set up the notation and mention certain properties of polymatroids. In Section 3, we present a randomized online algorithm with competitive ratio , where the run-time per timestep is exponential. In Section 4, we show how to achieve the same competitive ratio in polynomial time per timestep thereby proving Theorem 1. In Section 5, we show a simple online randomized algorithm for Disj-Conn-Spanning-Subhypergraphs thereby proving Theorem 2. Due to space limitations, certain proofs are deferred to the full version [3].
2 Preliminaries
Given two integers , let denote the set of integers with and denote the set of integers with . The operator with an unspecified base refers to . We recall that a polymatroid is an integer-valued monotone submodular function with its value on the empty set being . Let be a polymatroid. We denote the value of the function on the ground set by and recall that a set is a base if . We will assume that for every element : if we have an element with , then for every by submodularity and consequently, picking an arbitrary color for does not influence the number of base colors.
For a subset , we define the marginal function with respect to as the function defined by for every . If is submodular, then the function is also submodular for every . For ease of notation, for an element , we let denote the singleton set . For every set , we use and to abbreviate and , respectively.
We need the notion of span, closed sets, and quotients. For a set , the span of is the set of elements with marginal value with respect to , i.e., . We note that for every two sets , since is a monotone submodular function. A set is closed if . A set is a quotient of if for some set . This definition of quotient is equivalent to the one in Definition 3. We note that the empty set is closed (since for every element ) which implies that is a quotient of . The notion of quotients plays a central role in polymatroid sparsification [12].
For a set , let be the function obtained from by restricting the ground set to , i.e., for every . We recall that if is a polymatroid, then is also a polymatroid for every . Restricting the ground set preserves quotients as shown in the following lemma.
Lemma 4 ([12]).
Let be a polymatroid. Let be a quotient of and . Then, is a quotient of .
We recall that
(2) |
We note that there exists a closed set that achieves the minimum in the definition of : suppose is a set that achieves the minimum, then by definition of span while by submodularity of function and consequently, . If the function is clear from context, then we denote and by opt and respectively. We need the following approximate min-max relation between opt and .
Theorem 5 ([2]).
.
We need the following lemma showing that for every polymatroid , sampling each element independently with probability gives a base of with constant probability. A variant of the lemma was shown in [2]. We include a proof of the lemma in Appendix C for completeness.
Lemma 6.
Let be a polymatroid with . Let and be a subset obtained by picking each element in with probability at least independently at random. Then, is a base of with probability at least .
We recall that the elements of the ground set are indexed as according to the arrival order, i.e., is the element that arrives at time and moreover, is the set of elements that have arrived until time for every . For a coloring of all elements, we say that a fixed color is a base color if the set of elements of that color is a base of the polymatroid . Let be the number of base colors obtained by an online algorithm Alg for a given polymatroid . A deterministic online algorithm Alg is (purely) -competitive if for every polymatroid , we have
(3) |
A randomized online algorithm is -competitive if the bound in (3) holds in expectation.
3 An -Competitive Algorithm
In this section, we present a randomized online algorithm with expected competitive ratio but the run-time of the algorithm at each timestep will be exponential. The purpose of this section is to illustrate that the online arrival order has sufficient information to achieve a reasonable competitive ratio (albeit in exponential runtime) and to highlight the main ideas underlying the eventual polynomial-time algorithm that will prove Theorem 1. We describe the algorithm in Section 3.1. We present a novel property of the sequence of minimum sized quotients containing the last element with respect to a fixed ordering in Section 3.2. We prove the structural property that focusing on elements whose estimates are within a factor of does not decrease by more than a constant factor in Section 3.3. We use these properties to analyze the competitive ratio of the algorithm in Section 3.4. For Sections 3.1, 3.2, 3.3, and 3.4, we assume that and . We discuss how to relax the assumption in Section 3.5 in which we combine the designed algorithm with algorithms for other cases to obtain an algorithm with competitive ratio . We will modify the algorithm in Section 4 to design an -competitive algorithm that runs in polynomial time at each timestep, thereby completing the proof of Theorem 1.
3.1 Algorithm Description
We assume that and . The algorithm proceeds as follows: At each timestep , the algorithm computes
We note that is a quotient of with and hence, is well-defined. Let . Then, the algorithm samples a value from uniformly at random. Finally, the algorithm sets to be an integer picked uniformly at random from . We give a pseudocode of the algorithm in Algorithm 1.
3.2 Ordered Min-sized Quotients Property
In this section, we show that there can be at most distinct arrival times with the same value of .
Lemma 7.
Let be a polymatroid over the ground set with for every and . Let be an ordering of the ground set . For every , we define
Then, for every ,
Proof.
We note that is well-defined: empty set is closed since for every and hence, is a quotient of containing for every . Let . Let be the timesteps with . It suffices to show that . For every , let be a minimum sized quotient of containing element . We note that for every . We define .
We first show that for every . Suppose there exists an element . Then, there exists an index such that . Hence, , which implies that . Since , we have , which shows that . By Lemma 4, is a quotient of . We note that since and . Consequently, contradicts the fact that the smallest quotient of containing element has size exactly . Hence, for every .
We now show that for every . Let . Since is a quotient of containing element , we have . Since function is submodular, we have
Hence,
where the first inequality is by monotonicity of the function . Thus, we have for every .
This implies that is a strictly increasing integer sequence. We note that . Also, since is a monotone function, we have . Hence, .
3.3 Structural Property
For each , we define the element to be good if and bad otherwise. Let be the set of good elements. The following is the main result of this section. It shows that bases of are also bases of and moreover, the quantity is at least a constant fraction of .
Lemma 8.
We have that
-
1.
and
-
2.
.
We prove Lemma 8 in a series of steps. We note that an element could be bad because of two reasons: either is too large or it is too small. We first show that the only reason for certain elements being bad is that their is too small.
Lemma 9.
Let be a closed set of . Let be the smallest integer such that
(4) |
Then, for every element with , we have that .
Proof.
We recall that . By definition of , we have
Since is the smallest integer satisfying inequality (4), we have
(5) |
Hence,
Since is a closed set, is a quotient of . By Lemma 4, for every , we know that is a quotient of . Hence, for every element with , we have that , where the last inequality is by inequality (7).
Next, we show that dropping all bad elements will not decrease the value of by more than a constant factor for all sets. We first show this for closed sets (Lemma 10) and derive it for arbitrary sets as a corollary (Corollary 11).
Lemma 10.
For every closed set ,
Proof.
Let be a closed set. Let be the smallest integer such that inequality (4) holds. Lemma 9 implies that for every , element is bad if and only if . Hence, we have
(8) |
Thus,
Corollary 11.
For every set ,
Proof.
We have
(since by submodularity of ) | |||
We now prove Lemma 8.
Proof of Lemma 8.
We recall that by assumption. By applying Corollary 11 for , we conclude that . Hence,
where the last equality holds since for every and . By Corollary 11, for every set with , we have
Hence, .
3.4 Competitive Ratio Analysis
Now, we analyze the competitive ratio of Algorithm 1. Let . We note that . We now prove that every color is a base color with constant probability in Algorithm 1.
Lemma 12.
Let . Then, .
Proof.
Let . We recall that . By Lemma 8, we have that and hence, bases of correspond to bases of . By Lemma 8, we also have that . Let . We have and hence, . Therefore,
(9) |
Hence, for every element , the probability that is colored with is
We recall that . By applying Lemma 6 for the polymatroid , we conclude that the elements in with color form a base of with probability at least . We recall that bases of are also bases of . Hence, is a base color with probability at least .
Corollary 13.
The expected number of base colors obtained by Algorithm 1 is at least
Proof.
By Lemma 12, we have that
3.5 Combined Algorithm
We recall that our analysis of Algorithm 1 assumed that and . We now combine Algorithm 1 with other algorithms to address all ranges of and .
Consider the online algorithm Alg that runs Algorithm 1 with probability , assigns for every element with probability , and assigns for every element with probability . We show that the resulting online algorithm Alg has competitive ratio .
Theorem 14.
Algorithm Alg has competitive ratio .
Proof.
Suppose . Then, each singleton forms a base, which implies that . We recall that the algorithm assigns for every element with probability . This implies that
Suppose and . Then, the optimum is also smaller than by Theorem 5. We recall that the algorithm assigns for every element with probability . This implies that
Now, we may assume that and . We recall that the algorithm Alg runs Algorithm 1 with probability . We have
(10) |
which implies that
(11) |
Hence, we have
4 Approximation of in Polynomial Time
In Section 3, we presented a randomized online algorithm with competitive ratio . The algorithm takes exponential time, since the parameter is NP-hard to compute. In this section, we show how to get an approximation of in polynomial time via strength decomposition of the function .
Strength Decomposition.
We first introduce the definition of strength decomposition. Let be a polymatroid over the ground set . For a subset of , the strength-ratio of in , denoted , is the value , with the convention that for every .
Definition 15.
[12] Let be a polymatroid over the ground set . A strength decomposition of with respect to is a sequence of sets such that: (i) , (ii) and for every , (iii) for every , , and (iv) the strength ratios are nondecreasing in .
For every polymatroid over a ground set , there exists a strength decomposition of with respect to and it can be computed in polynomial time via reduction to submodular minimization [9, 10, 8] (also see [12]). We further show that all sets in a strength decomposition are closed.
Lemma 16.
Let be a polymatroid over the ground set . Let be a strength decomposition with respect to . Then, is a closed set for every .
Approximating in Polynomial Time.
Now, we show how to approximate in polynomial time. For every timestep , let be a strength decomposition with respect to . Let be the index such that . We define . We prove that the parameter is an -approximation of .
Lemma 17.
For every , .
5 Simpler Algorithm for Online Disj-Conn-Spanning-Subhypergraphs
In this section, we consider Disj-Conn-Spanning-Subhypergraphs. The input here is a connected hypergraph where each hyperedge is a non-empty subset of . We use and to denote the number of vertices and hyperedges in the hypergraph respectively. We note that is a multi-set, i.e., it can contain multiple copies of the same hyperedge. We recall that the hypergraph can equivalently be represented as a bipartite graph , where we have a node for each vertex and each hyperedge and an edge between nodes and if the hyperedge contains the vertex . A hypergraph is connected if its bipartite representation is connected. A subhypergraph of is a subset of hyperedges and it is said to be spanning if and connected if is connected. The goal is to find a maximum number of disjoint connected spanning subhypergraphs.
In the online model for Disj-Conn-Spanning-Subhypergraphs, the set of nodes is known in advance while the hyperedges arrive in an online fashion. For a color assignment , we say that a color is a connected spanning color if the set of hyperedges with color is spanning and connected. An online algorithm has to color each hyperedge immediately upon arrival irrevocably in order to maximize the number of connected spanning colors. Theorem 1 implies a randomized online algorithm for this problem with competitive ratio that runs in polynomial time. We show that we can achieve the same competitive ratio via a much simpler estimator.
For every integer , let be the hyperedge arriving at time . Define as the set of hyperedges that have arrived until time . For a subset of nodes, let be the set of hyperedges containing all vertices in that have arrived until time . We define
Equivalently, is the min over all pairs of vertices in the current hyperedge of the number of hyperedges containing the pair in the current hypergraph. We prove that replacing the estimate in Algorithm 1 by gives a randomized online algorithm with the same competitive ratio. We note that for every can be maintained in time in the -th timestep, which implies that the algorithm takes time in the -th timestep.
6 Conclusion
Packing bases of a polymatroid generalizes numerous set packing problems including disjoint bases of a matroid, disjoint spanning trees of a graph, disjoint set covers of a set system, and disjoint connected spanning subhypergraphs of a hypergraph. In this work, we introduced an online model for packing bases of a polymatroid and gave an algorithm that achieves a polylogarithmic competitive ratio. Our algorithm leads to a polylogarithmic competitive ratio for all these packing problems in a unified fashion. Our algorithm is based on novel properties of the notion of quotients of polymatroids. For the special cases of disjoint spanning trees of a graph and disjoint connected spanning subhypergraphs of a hypergraph, we gave a simpler and more elegant online algorithm that is also easy to implement while achieving the same polylogarithmic competitive ratio.
Our work leads to several interesting open questions. We mention three prominent ones: Firstly, we recall that Disj-Spanning-Trees (and more generally, Disj-Matroid-Bases) in the offline setting is polynomial-time solvable via matroidal techniques. Is it possible to design an online algorithm for Disj-Spanning-Trees (and more generally, for Disj-Matroid-Bases) with constant competitive ratio? Secondly, it is known that there is no randomized algorithm with expected impure competitive ratio for Disj-Set-Cover, where is the size of the ground set. Can we show that there is no randomized algorithm with expected impure competitive ratio or pure competitive ratio for Disj-Polymatroid-Bases? Thirdly, we recall that our online model for Disj-Polymatroid-Bases assumes knowledge of . Is it possible to achieve polylogarithmic competitive ratio without knowledge of ? This seems to be open even in the special case of coverage functions (i.e., for online Disj-Set-Cover) [6].
References
- [1] M. Bienkowski, J. Byrka, and . Jeż. Online Disjoint Set Covers: Randomization Is Not Necessary. In 42nd International Symposium on Theoretical Aspects of Computer Science (STACS), LIPIcs, pages 18:1–18:16, 2025.
- [2] G. Călinescu, C. Chekuri, and J. Vondrák. Disjoint bases in a polymatroid. Random Struct. Algorithms, 35(4):418–430, 2009. doi:10.1002/RSA.20274.
- [3] K. Chandrasekaran, C. Chekuri, and W. Zhu. Online disjoint spanning trees and polymatroid bases. arXiv preprint arXiv:2503.19999, 2025.
- [4] J. Cheriyan and M. R. Salavatipour. Packing element-disjoint steiner trees. ACM Trans. Algorithms, 3(4):47, 2007. doi:10.1145/1290672.1290684.
- [5] J. Edmonds. Minimum partition of a matroid into independent subsets. J. Res. Nat. Bur. Standards B, 69:73–77, 1965.
- [6] Y. Emek, A. Goldbraikh, and E. Kantor. Online Disjoint Set Cover Without Prior Knowledge. In 27th Annual European Symposium on Algorithms (ESA), pages 44:1–44:16, 2019.
- [7] U. Feige, M. M. Halldórsson, G. Kortsarz, and A. Srinivasan. Approximating the domatic number. SIAM Journal on Computing, 32(1):172–195, 2002. doi:10.1137/S0097539700380754.
- [8] S. Fujishige. Theory of principal partitions revisited. Research Trends in Combinatorial Optimization, pages 127–162, 2009.
- [9] H. Narayanan. The principal lattice of partitions of a submodular function. Linear Algebra and its Applications, 144:179–216, 1991.
- [10] H. Narayanan. Submodular functions and electrical networks, volume 54. Elsevier, 1997.
- [11] A. Pananjady, V. K. Bagaria, and R. Vaze. The online disjoint set cover problem and its applications. In 2015 IEEE Conference on Computer Communications (INFOCOM), pages 1221–1229, 2015.
- [12] K. Quanrud. Quotient sparsification for submodular functions. In Proceedings of the 2024 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 5209–5248, 2024.
- [13] A. Schrijver. Combinatorial Optimization: Polyhedra and Efficiency. Springer, 2003.
Appendix A Quotients for some concrete polymatroids
For a polymatroid , we recall that a set is a quotient of if , i.e., if each element has strictly positive marginal with respect to . Since the definition of quotients might seem contrived, we interpret it for some concrete polymatroids:
-
1.
Let be a connected hypergraph (graph respectively). Consider the polymatroid that arises in Disj-Conn-Spanning-Subhypergraphs (Disj-Spanning-Trees respectively) with being the input instance. Quotients of correspond to union of cut-set of disjoint subsets of vertices, i.e., each quotient of is of the form , where the family is a disjoint family of subsets of vertices. In particular, the minimum size of a non-empty quotient is equal to the global min-cut value of .
-
2.
Let be a hypergraph. Consider the polymatroid that arises in Disj-Set-Cover with being the input instance, i.e., is the coverage function of . Quotients of correspond to union of vertex-isolating cuts of , i.e., each quotient of is of the form where . In particular, the minimum size of a non-empty quotient is equal to the min-degree of .
Appendix B Min-size of a non-empty Quotient and
Let be the minimum size of non-empty quotients of , i.e.,
We recall that
The following lemma shows that is an -approximation of .
Lemma 18.
.
Proof.
We first show the lower bound on . Consider a closed set with and . Since is closed, we have that for every and is a quotient of . Therefore, we have
We now show the upper bound on . Consider a minimum sized non-empty quotient . Let . Since is a closed set and , we have . We note that for every , , where the inequality is by monotonicity of the function . Therefore,
Appendix C Random Sampling gives a Base
Proof.
If , then we have , which implies that which is a base of . Henceforth, we may assume that and consequently, . We recall that . Let be a uniformly random permutation of the elements in . For every , we define and . We note that and .
Let . We consider the distribution of element conditioned on . Since is a uniformly random permutation, can be an arbitrary element in , of which each has the same probability. That is, for every element ,
We note that element is picked with probability at least . If element is picked, then . Otherwise, . Hence, we have
where the last inequality is by the definition of . This implies that
By setting , we have
According to Markov’s inequality, we have
which shows that is a base with probability at least .