Abstract 1 Introduction 2 Preliminaries 3 An 𝑶(𝐥𝐨𝐠𝟐𝒓)-Competitive Algorithm 4 Approximation of 𝒒𝒕 in Polynomial Time 5 Simpler Algorithm for Online Disj-Conn-Spanning-Subhypergraphs 6 Conclusion References Appendix A Quotients for some concrete polymatroids Appendix B Min-size of a non-empty Quotient and 𝒌 Appendix C Random Sampling gives a Base

Online Disjoint Spanning Trees and Polymatroid Bases

Karthekeyan Chandrasekaran ORCID Grainger College of Engineering, University of Illinois, Urbana-Champaign, IL, USA Chandra Chekuri Grainger College of Engineering, University of Illinois, Urbana-Champaign, IL, USA Weihao Zhu Grainger College of Engineering, University of Illinois, Urbana-Champaign, IL, USA
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 O(logn)-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 Algorithms
Category:
Track A: Algorithms, Complexity and Games
Copyright and License:
[Uncaptioned image] © Karthekeyan Chandrasekaran, Chandra Chekuri, and Weihao Zhu; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Approximation algorithms analysis
; Theory of computation Online algorithms
Related Version:
Full Version: https://arxiv.org/abs/2503.19999 [3]
Funding:
This work was supported in part by NSF grant CCF-2402667.
Editors:
Keren Censor-Hillel, Fabrizio Grandoni, Joël Ouaknine, and Gabriele Puppis

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 O(logn)-approximation111We are using the convention of α-approximation for a maximization problem with α>1 with the understanding that the returned value is at least opt/α. [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 f:2𝒩0 on a ground set 𝒩 is an integer-valued monotone submodular function with f()=0. We recall that a function f:2𝒩0 is monotone if f(A)f(B) for every AB and submodular if f(A)+f(B)f(AB)+f(AB) for every A,B𝒩. A subset S𝒩 is a base of f if f(S)=f(𝒩).222In most settings, one would define a set S to be a base if it is inclusionwise minimal subject to satisfying f(S)=f(𝒩). 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 f:2𝒩0 via an evaluation oracle. The goal is to find a maximum number of disjoint bases. We define

opt(f):=max{k:k disjoint bases of f}.

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 f:2𝒩0 over a large but finite ground set 𝒩. Let n:=|𝒩|. We index the elements of the ground set 𝒩 as e1,e2,,en where et is the element that arrives at time t for every t[n]. For each t[n], we denote 𝒩t:={e1,e2,,et} and let f|𝒩t:2𝒩t0 be the function obtained from f by restricting the ground set to 𝒩t, that is, f|𝒩t(A):=f(A) for every A𝒩t. At each timestep t[n], the online algorithm has access to the evaluation oracle333The evaluation oracle of a function g:2V takes a set SV as input and returns g(S). of the function restricted to the set of elements that have arrived until time t, i.e., the evaluation oracle of the function f|𝒩t, and has to color element et irrevocably. A color is said to be a base color if the set S of elements with that color is a base of f. 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., f(𝒩). 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. 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 H=(V,E) where V corresponds to the universe and the hyperedges in E correspond to sets in the system. The goal is to find a maximum number of disjoint set covers (a subset AE of hyperedges is a set cover if eAe=V). 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 f:2E0 defined as f(A):=|eAe|. 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. 2.

    In the Disj-Matroid-Bases problem, we are given evaluation access to a matroid rank function r:2𝒩0 over a ground set 𝒩 (we recall that a matroid rank function r is a polymatroid that additionally satisfies r({e})1 for every e𝒩). A subset S𝒩 is a base of the matroid if r(S)=r(𝒩). 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. 1.

    In the Disj-Matrix-Bases problem, we are given a matrix Mn×d of rank d 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 d. This is a special case of Disj-Matroid-Bases where the matroid is the linear matroid defined by M (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. 2.

    In the Disj-Spanning-Trees problem, we are given a connected graph G=(V,E) and the goal is to find a maximum number of disjoint spanning trees in G. Disj-Spanning-Trees is a special case of Disj-Matroid-Bases where the matroid is the graphic matroid defined by G (and consequently, the rank function of a subset FE is |V|c(V,F), where c(V,F) is the number of components in the graph (V,F)). In the online setting of Disj-Spanning-Trees (termed Online-Disj-Spanning-Trees), the vertex set of G is known in advance while the edges of G 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 V.

  3. 3.

    In the Disj-Conn-Spanning-Subhypergraphs problem, we are given a connected hypergraph H=(V,E) and the goal is to find a maximum number of disjoint connected spanning subhypergraphs in H. Disj-Conn-Spanning-Subhypergraphs is a special case of Disj-Polymatroid-Bases where the polymatroid f:2E0 of interest is defined as f(A):=|V|c(V,A) for every AE, where c(V,A) is the number of components in the hypergraph (V,F). 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 G is known in advance while the hyperedges of G 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 V.

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 H=(V,E) of Disj-Set-Cover, add a new vertex r and stretch every hyperedge of H to include r to obtain a new hypergraph H=(V+r,E); the maximum number of disjoint set covers in H is equal to the maximum number of disjoint connected spanning subhypergraphs in H.

Prior work in the offline setting.

Disj-Matroid-Bases is polynomial-time solvable [5]. Disj-Set-Cover and Disj-Conn-Spanning-Subhypergraphs are o(log|V|)-inapproximable and O(log|V|)-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 O(logf(𝒩))-approximation for Disj-Polymatroid-Bases by showing an approximate min-max relation for opt(f). Their approximate min-max relation is based on the following minimization problem:

k(f):=minA𝒩:f(A)<f(𝒩)e𝒩(f(A+e)f(A))f(𝒩)f(A).

It is easy to see that opt(f)k(f) (e.g., see [2]). Moreover, if the polymatroid f is a matroid rank function, then Edmonds [5] showed that opt(f)=k(f). Edmonds’ result is constructive and implies a polynomial time algorithm for Disj-Matroid-Bases. However, for coverage functions, it is known that opt(f)(1+o(1))k(f)/logf(𝒩) [7]. Călinescu, Chekuri, and Vondrák showed that this bound is tight by giving a polynomial-time algorithm to construct (1o(1))k(f)/logf(𝒩) disjoint bases (and hence, opt(f)(1o(1))k(f)/logf(𝒩)). 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 k:=k(f)/(logf(𝒩)+loglogf(𝒩)) (by guessing/binary search) and colors each element with a uniformly random color chosen from a color palette of size k. Călinescu, Chekuri, and Vondrák showed that the expected number of base colors returned by this random coloring algorithm is at least (1e/logf(𝒩))k=(1o(1))k(f)/logf(𝒩). 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 k(f) associated with the coverage function f of the hypergraph H=(V,E) has a simple interpretation: it is equal to the min-degree of the hypergraph H. 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 O(logn), where n 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 Ω(n). 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 O(log2n) (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 (opt(H)/α)β for some β>0 that is a function only of n, where H is the input hypergraph, n is the number of vertices in H, and opt(H) is the maximum number of disjoint set covers in H. Our work focuses on the case of β=0, i.e., pure competitive ratio. Ω(logn/loglogn) (even with knowledge of min-degree). Secondly, Bienkowski, Byrka, and Jeż [1] designed an online deterministic algorithm with impure competitive ratio O(log2n) (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 k(f) is known in advance, then it is possible to design a randomized online algorithm with expected competitive ratio O(logf(𝒩)): indeed, the random coloring algorithm of [2] mentioned above can be implemented in the online setting using knowledge of k(f) 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 k(f) 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 2-dimensional vectors, i.e., for d=2. 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 c=1; at each timestep t[n]: use color c for element et and if the set of elements with color c is a base, then increment c. For uniform random arrival order, the competitive ratio of this simple online algorithm is O(logf(𝒩)) (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 k(f): the algorithm will return at least one base while the maximum number of possible bases is at most k(f) (since opt(f)k(f)). It is known that k(f) is the best possible competitive ratio of deterministic online algorithms that do not have prior knowledge of k(f) [11]. In this work, we are interested in the setting of arbitrary arrival order without prior knowledge of k(f). For this setting, we design a randomized online algorithm with expected competitive ratio O(log2f(𝒩)).

Theorem 1.

For Disj-Polymatroid-Bases, there exists a randomized online algorithm with expected competitive ratio O(log2f(𝒩)). The runtime of the algorithm at each timestep t[n] is poly(t,logf(𝒩)).

We recall that the best-known approximation factor for Disj-Polymatroid-Bases in the offline setting is O(logf(𝒩)), 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 O(log2n), where n 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 O(log2r), where r is the rank of the ground set. Consequently, we obtain a randomized online algorithm with expected competitive ratio O(log2d) for Disj-Matrix-Bases where d is the dimension of the span of the input vectors and with expected competitive ratio O(log2n) for Disj-Spanning-Trees where n 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 O(log2n), where n 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 O(log2f(𝒩)) 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 O(log2n), where n is number of vertices in the input hypergraph. The runtime of the algorithm at each timestep t is O(|et|2), where et is the hyperedge that arrives at time t. 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 f:2𝒩0 be the polymatroid of interest. Let r:=f(𝒩) be the function value of bases and k:=k(f). For the discussion here, we will assume that f(e)>0 for every e𝒩, r2, and k=Ω(log2r) since this is the non-trivial case. First consider the setting where the online algorithm knows k in advance. In this setting, coloring each element uniformly from the palette [Θ(k/logr)] gives Ω(k/logr) base colors in expectation and, consequently achieves O(logr)-competitive ratio.

Suppose k is not known in advance. Our approach is to estimate k at each timestep. At each timestep t[n], suppose that we can find a number qt that is within a poly(r) factor of k, i.e., qt[kpoly(r),kpoly(r)] – we call such a qt as a coarse estimate. Let Pt be a uniform random sample from {qt2O(logr),,qt21,qt20,qt21,,qt2O(logr)}; Pt is a 2-approximation of k with probability Ω(1/logr). Consequently, coloring the element et uniformly from the palette [Θ(Pt/log2r)] leads to Ω(k/log2r) base colors in expectation.

The key challenge is computing a coarse estimate for k 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 k. We formalize this approach now. Suppose we have a subroutine to compute some estimate qt at each timestep t[n]. We call a timestep t[n] to be good if the estimate qt is a coarse estimate of k and bad otherwise. We let 𝒩𝗀𝗈𝗈𝖽 be the collection of elements that arrive at good timesteps and focus on the function f restricted to 𝒩𝗀𝗈𝗈𝖽, i.e., on the function g:=f|𝒩𝗀𝗈𝗈𝖽. Suppose that the following two properties hold: (i) Bases of g are also bases of f and (ii) k(g)=Ω(1)k. These two properties suffice to obtain Ω(k/log2r) 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 h:2𝒱0, a set Q𝒱 is a quotient of h if h(e+(𝒱Q))>h(𝒱Q) eQ, i.e., if each element eQ has strictly positive marginal with respect to 𝒱Q.

See Appendix A for the interpretation of quotients for some concrete polymatroids. We observed that the minimum size of a non-empty quotient of f is an r-approximation of k – see Appendix B. This inspired us to use the minimum-sized quotient of f|𝒩t that contains et as the estimate at each timestep t[n]. That is, we use the estimator defined as

qt :=min{|Q|:etQ𝒩t and Q is a quotient of f|𝒩t}. (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 O(log2r)-competitive ratio. A technical issue is that the estimator qt is NP-Hard to compute for general polymatroids (even for matroid rank functions). However, we note that a poly(r)-approximation of qt is sufficient for the desired competitive ratio. We show that an r-approximation of qt can be computed using a strength decomposition with respect to f|𝒩t 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 f:2E0 of interest here is the rank function of the graphic matroid defined by the undirected graph G=(V,E). For this case, the quotient computation is polynomial-time solvable since qt is equal to the min ut-vt cut value in the graph (V,Et) where et=(ut,vt) and Et={e1,e2,,et}. However, we show that a much simpler estimator for Disj-Spanning-Trees also achieves properties (i) and (ii) mentioned above: we use

ηt:=number of edges between ut and vt in the graph (V,Et).

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 O(log2r), 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 ab, let [a,b] denote the set of integers x with axb and [b] denote the set of integers x with 1xb. The log() operator with an unspecified base refers to log2(). We recall that a polymatroid is an integer-valued monotone submodular function with its value on the empty set being 0. Let f:2𝒩0 be a polymatroid. We denote the value of the function f on the ground set 𝒩 by r:=f(𝒩) and recall that a set S𝒩 is a base if f(S)=r. We will assume that f(e)>0 for every element e𝒩: if we have an element e𝒩 with f(e)=0, then f(S)=f(S+e) for every S𝒩 by submodularity and consequently, picking an arbitrary color for e does not influence the number of base colors.

For a subset A𝒩, we define the marginal function with respect to A as the function fA:2𝒩 defined by fA(S):=f(AS)f(A) for every S𝒩. If f is submodular, then the function fA is also submodular for every A𝒩. For ease of notation, for an element e𝒩, we let e denote the singleton set {e}. For every set S𝒩, we use S+e and Se to abbreviate S{e} and S{e}, respectively.

We need the notion of span, closed sets, and quotients. For a set S𝒩, the span of S is the set of elements with marginal value 0 with respect to S, i.e., span(S):={e𝒩:fS(e)=0}. We note that for every two sets AB, span(A)span(B) since f is a monotone submodular function. A set S𝒩 is closed if S=span(S). A set Q𝒩 is a quotient of f if Q=𝒩span(S) for some set S𝒩. This definition of quotient is equivalent to the one in Definition 3. We note that the empty set is closed (since f(e)>0 for every element e𝒩) which implies that 𝒩 is a quotient of f. The notion of quotients plays a central role in polymatroid sparsification [12].

For a set S𝒩, let f|S:2S0 be the function obtained from f by restricting the ground set to S, i.e., f|S(T):=f(T) for every TS. We recall that if f is a polymatroid, then f|S is also a polymatroid for every S𝒩. Restricting the ground set preserves quotients as shown in the following lemma.

Lemma 4 ([12]).

Let f:2𝒩0 be a polymatroid. Let Q be a quotient of f and S𝒩. Then, QS is a quotient of f|S.

We recall that

opt(f) =max{k:k disjoint bases of f}and
k(f) =minA𝒩:f(A)<f(𝒩)e𝒩fA(e)f(𝒩)f(A). (2)

We note that there exists a closed set A that achieves the minimum in the definition of k(f): suppose A is a set that achieves the minimum, then f(span(A))=f(A)<f(𝒩) by definition of span while e𝒩fA(e)e𝒩fspan(A)(e) by submodularity of function f and consequently, e𝒩fA(e)f(𝒩)f(A)e𝒩fspan(A)(e)f(𝒩)f(span(A)). If the function f is clear from context, then we denote opt(f) and k(f) by opt and k respectively. We need the following approximate min-max relation between opt and k.

Theorem 5 ([2]).

kO(logr)optk.

We need the following lemma showing that for every polymatroid f, sampling each element independently with probability Ω(clogrk) gives a base of f 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 f:2𝒩0 be a polymatroid with f(𝒩)=r2. Let p:=min{1,2logrk} and S𝒩 be a subset obtained by picking each element in 𝒩 with probability at least p independently at random. Then, S is a base of f with probability at least 12.

We recall that the elements of the ground set 𝒩 are indexed as e1,e2,,en according to the arrival order, i.e., et is the element that arrives at time t[n] and moreover, 𝒩t={e1,e2,,et} is the set of elements that have arrived until time t for every t[n]. 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 f. Let Alg(f) be the number of base colors obtained by an online algorithm Alg for a given polymatroid f. A deterministic online algorithm Alg is (purely) α-competitive if for every polymatroid f, we have

Alg(f)opt(f)α. (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 O(log2r) 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 poly(r) factor of k does not decrease k 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 r2 and k120log2r. 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 O(log2r). We will modify the algorithm in Section 4 to design an O(log2r)-competitive algorithm that runs in polynomial time at each timestep, thereby completing the proof of Theorem 1.

3.1 Algorithm Description

We assume that k120log2r and r2. The algorithm proceeds as follows: At each timestep t[n], the algorithm computes

qt:=min{|Q|:etQ𝒩t and Q is a quotient of f|𝒩t}.

We note that 𝒩t is a quotient of f|𝒩t with et𝒩t and hence, qt is well-defined. Let t:=logqt. Then, the algorithm samples a value Rt from [t3logr,t+3logr] uniformly at random. Finally, the algorithm sets C(et) to be an integer picked uniformly at random from [2Rt/(60log2r)]. We give a pseudocode of the algorithm in Algorithm 1.

Algorithm 1 Randomized Online Algorithm for Disj-Polymatroid-Bases.

3.2 Ordered Min-sized Quotients Property

In this section, we show that there can be at most r distinct arrival times t with the same value of qt.

Lemma 7.

Let f:2𝒩0 be a polymatroid over the ground set 𝒩 with f(e)>0 for every e𝒩 and f()=0. Let e1,e2,,en be an ordering of the ground set 𝒩. For every t[n], we define

qt:=min{|Q|:etQ𝒩t and Q is a quotient of f|𝒩t}.

Then, for every j+,

|{t[n]:qt=j}|r.

Proof.

We note that qt is well-defined: empty set is closed since f(e)>0 for every e𝒩 and hence, 𝒩t is a quotient of f|𝒩t containing et for every t[n]. Let j+. Let t1<t2<<t be the timesteps t with qt=j. It suffices to show that r. For every i[], let Qti be a minimum sized quotient of f|𝒩ti containing element eti. We note that |Qti|=qti=j for every i[]. We define S:=i=1(Qtieti).

We first show that etiS for every i[]. Suppose there exists an element eti1S. Then, there exists an index i2[] such that eti1Qti2eti2. Hence, eti1Qti2eti2𝒩ti2, which implies that i1i2. Since eti1Qti1eti1, we have i1i2, which shows that i1<i2. By Lemma 4, Qti2𝒩ti1 is a quotient of f|𝒩ti1. We note that |Qti2𝒩ti1|<|Qti2|=j since eti2Qti2 and eti2Qti2𝒩ti1. Consequently, Qti2𝒩ti1 contradicts the fact that the smallest quotient of f|𝒩ti1 containing element eti1 has size exactly j. Hence, etiS for every i[].

We now show that f(𝒩tiS)<f(𝒩ti+1S) for every i[1]. Let i[1]. Since Qti+1 is a quotient of 𝒩ti+1 containing element eti+1, we have f𝒩ti+1Qti+1(eti+1)>0. Since function f is submodular, we have

f(𝒩ti+1S)f(𝒩ti+11S)=f𝒩ti+11S(eti+1)f𝒩ti+1Qti+1(eti+1)>0.

Hence,

f(𝒩tiS)f(𝒩ti+11S)<f(𝒩ti+1S),

where the first inequality is by monotonicity of the function f. Thus, we have f(𝒩tiS)<f(𝒩ti+1S) for every i[1].

This implies that {f(𝒩tiS)}i=1 is a strictly increasing integer sequence. We note that f(𝒩tS)r. Also, since f is a monotone function, we have f(𝒩t1S)f(et1)>0. Hence, r.

3.3 Structural Property

For each t[n], we define the element et to be good if k2r<qt<2rk and bad otherwise. Let 𝒩𝗀𝗈𝗈𝖽 be the set of good elements. The following is the main result of this section. It shows that bases of f|𝒩𝗀𝗈𝗈𝖽 are also bases of f and moreover, the quantity k(f|𝒩𝗀𝗈𝗈𝖽) is at least a constant fraction of k(f).

Lemma 8.

We have that

  1. 1.

    f(𝒩𝗀𝗈𝗈𝖽)=f(𝒩) and

  2. 2.

    k(f|𝒩𝗀𝗈𝗈𝖽)12k(f).

We prove Lemma 8 in a series of steps. We note that an element et could be bad because of two reasons: either qt is too large or it is too small. We first show that the only reason for certain elements being bad is that their qt is too small.

Lemma 9.

Let A𝒩 be a closed set of f. Let T1 be the smallest integer such that

e𝒩TA(f(A+e)f(A))k(f(𝒩)f(A)). (4)

Then, for every element et𝒩TA with t[T], we have that qt<2rk.

Proof.

We recall that k=k(f). By definition of k, we have

e𝒩A(f(A+e)f(A))k(f(𝒩)f(A)).

Since T is the smallest integer satisfying inequality (4), we have

e𝒩T1A(f(A+e)f(A))<k(f(𝒩)f(A)). (5)

Hence,

e𝒩TA(f(A+e)f(A)) =(e𝒩T1A(f(A+e)f(A)))+(f(A+eT)f(A))
<k(f(𝒩)f(A))+(f(A+eT)f(A))(by inequality (5))
k(f(𝒩)f(A))+(f(𝒩)f(A))(since f(A+eT)f(𝒩))
=(k+1)(f(𝒩)f(A)). (6)

Since A is a closed set, we have f(A+e)f(A)1 for every element e𝒩TA. Thus, using inequality (3.3), we have

|𝒩TA|<e𝒩TA(f(A+e)f(A))(k+1)(f(𝒩)f(A))(k+1)r<2rk. (7)

Since A is a closed set, 𝒩A is a quotient of f. By Lemma 4, for every t[n], we know that 𝒩tA is a quotient of f|𝒩t. Hence, for every element et𝒩TA with t[T], we have that qt|𝒩tA||𝒩TA|<2rk, where the last inequality is by inequality (7).

Next, we show that dropping all bad elements will not decrease the value of k 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 A𝒩,

e𝒩𝗀𝗈𝗈𝖽A(f(A+e)f(A))k2(f(𝒩)f(A)).

Proof.

Let A𝒩 be a closed set. Let T1 be the smallest integer such that inequality (4) holds. Lemma 9 implies that for every t[T], element et𝒩TA is bad if and only if qtk2r. Hence, we have

|𝒩T(A𝒩𝗀𝗈𝗈𝖽)| =|{t[T]:et𝒩T(A𝒩𝗀𝗈𝗈𝖽)}|=|{t[T]:et𝒩TA&qtk2r}|
|{t[T]:qtk2r}|=j=1k2r|{t[T]:qt=j}|
(k2r)r(by Lemma 7)
=k2. (8)

Thus,

e𝒩𝗀𝗈𝗈𝖽A(f(A+e)f(A))
e(𝒩T𝒩𝗀𝗈𝗈𝖽)A(f(A+e)f(A))(since f(A+e)f(A) for every e𝒩)
=e𝒩TA(f(A+e)f(A))e𝒩T(A𝒩𝗀𝗈𝗈𝖽)(f(A+e)f(A))
k(f(𝒩)f(A))e𝒩T(A𝒩𝗀𝗈𝗈𝖽)(f(A+e)f(A))(by inequality (4))
k(f(𝒩)f(A))e𝒩T(A𝒩𝗀𝗈𝗈𝖽)(f(𝒩)f(A))(by monotonicity of f)
k(f(𝒩)f(A))k2(f(𝒩)f(A))(by inequality (3.3))
=k2(f(𝒩)f(A)).

Corollary 11.

For every set A𝒩,

e𝒩𝗀𝗈𝗈𝖽A(f(A+e)f(A))k2(f(𝒩)f(A)).

Proof.

We have

e𝒩𝗀𝗈𝗈𝖽A(f(A+e)f(A)) e𝒩𝗀𝗈𝗈𝖽span(A)(f(A+e)f(A))(since Aspan(A))
e𝒩𝗀𝗈𝗈𝖽span(A)(f(span(A)+e)f(span(A))
(since fA(e)fspan(A)(e) by submodularity of f)
k2(f(𝒩)f(span(A)))(by applying Lemma 10 for span(A))
=k2(f(𝒩)f(A)).(since f(A)=f(span(A)))

We now prove Lemma 8.

Proof of Lemma 8.

We recall that k=k(f)>0 by assumption. By applying Corollary 11 for A=𝒩𝗀𝗈𝗈𝖽, we conclude that f(𝒩)=f(𝒩𝗀𝗈𝗈𝖽). Hence,

k(f|𝒩𝗀𝗈𝗈𝖽)=minA𝒩𝗀𝗈𝗈𝖽:f(A)<f(𝒩𝗀𝗈𝗈𝖽)e𝒩𝗀𝗈𝗈𝖽fA(e)f(𝒩𝗀𝗈𝗈𝖽)f(A)=minA𝒩𝗀𝗈𝗈𝖽:f(A)<f(𝒩𝗀𝗈𝗈𝖽)e𝒩𝗀𝗈𝗈𝖽AfA(e)f(𝒩)f(A),

where the last equality holds since fA(e)=0 for every eA and f(𝒩𝗀𝗈𝗈𝖽)=f(𝒩). By Corollary 11, for every set A𝒩 with f(A)<f(𝒩𝗀𝗈𝗈𝖽), we have

e𝒩𝗀𝗈𝗈𝖽AfA(e)f(𝒩)f(A)k2.

Hence, k(f|𝒩𝗀𝗈𝗈𝖽)12k.

3.4 Competitive Ratio Analysis

Now, we analyze the competitive ratio of Algorithm 1. Let h:=logk. We note that k2<2hk. We now prove that every color c[2h/(60log2r)] is a base color with constant probability in Algorithm 1.

Lemma 12.

Let c[2h/(60log2r)]. Then, 𝐏𝐫Alg1[c is a base color]12.

Proof.

Let k𝗀𝗈𝗈𝖽:=k(f|𝒩𝗀𝗈𝗈𝖽). We recall that k=k(f). By Lemma 8, we have that f(𝒩𝗀𝗈𝗈𝖽)=f(𝒩) and hence, bases of f|𝒩𝗀𝗈𝗈𝖽 correspond to bases of f. By Lemma 8, we also have that k𝗀𝗈𝗈𝖽k22h1. Let et𝒩𝗀𝗈𝗈𝖽. We have k2r<qt<2rk and hence, logqt3logrhlogqt+3logr. Therefore,

𝐏𝐫[Rt=h]=16logr+1>115logr. (9)

Hence, for every element et𝒩𝗀𝗈𝗈𝖽, the probability that et is colored with c is

𝐏𝐫[C(et)=c] 𝐏𝐫[Rt=h]𝐏𝐫[C(et)=c|rt=h]
=𝐏𝐫[Rt=h]12h/(60log2r)
>115logr60log2r2h(by inequality (9))
115logr60log2r2k𝗀𝗈𝗈𝖽=2logrk𝗀𝗈𝗈𝖽.

We recall that k𝗀𝗈𝗈𝖽=k(f|𝒩𝗀𝗈𝗈𝖽). By applying Lemma 6 for the polymatroid f|𝒩𝗀𝗈𝗈𝖽, we conclude that the elements in 𝒩𝗀𝗈𝗈𝖽 with color c form a base of f|𝒩𝗀𝗈𝗈𝖽 with probability at least 12. We recall that bases of f|𝒩𝗀𝗈𝗈𝖽 are also bases of f. Hence, c is a base color with probability at least 12.

Lemma 12 implies a lower bound on the expected number of base colors obtained by Algorithm 1.

Corollary 13.

The expected number of base colors obtained by Algorithm 1 is at least

𝔼[Alg1(f)]122h/(60log2r).

Proof.

By Lemma 12, we have that

𝔼[Alg1(f)]=c[2h/(60log2r)]𝐏𝐫Alg1[c is a base color]122h/(60log2r).

3.5 Combined Algorithm

We recall that our analysis of Algorithm 1 assumed that r2 and k120log2r. We now combine Algorithm 1 with other algorithms to address all ranges of r and k.

Consider the online algorithm Alg1 that runs Algorithm 1 with probability 13, assigns C(et)=1 for every element et with probability 13, and assigns C(et)=t for every element et with probability 13. We show that the resulting online algorithm Alg1 has competitive ratio O(log2r).

Theorem 14.

Algorithm Alg1 has competitive ratio O(log2r).

Proof.

Suppose r=1. Then, each singleton forms a base, which implies that opt(f)=n. We recall that the algorithm assigns C(et)=t for every element et with probability 13. This implies that

𝔼[Alg1(f)]13opt(f).

Suppose r2 and k<120log2r. Then, the optimum is also smaller than 120log2r by Theorem 5. We recall that the algorithm assigns C(et)=1 for every element et with probability 13. This implies that

𝔼[Alg1(f)]13>1360log2ropt(f).

Now, we may assume that k120log2r and r2. We recall that the algorithm Alg1 runs Algorithm 1 with probability 13. We have

2h/(60log2r)>160k2log2r1, (10)

which implies that

2h/(60log2r)122h/(60log2r). (11)

Hence, we have

𝔼[Alg1(f)] 13𝔼[Alg1(f)]
162h/(60log2r)(by Corollary 13)
1122h/(60log2r)(by inequality (11))
11440log2rk
11440log2ropt(f).(by Theorem 5)

4 Approximation of 𝒒𝒕 in Polynomial Time

In Section 3, we presented a randomized online algorithm with competitive ratio O(log2r). The algorithm takes exponential time, since the parameter qt=min{|Q|:etQ𝒩t and Q is a quotient of f|𝒩t} is NP-hard to compute. In this section, we show how to get an approximation of qt in polynomial time via strength decomposition of the function f.

Strength Decomposition.

We first introduce the definition of strength decomposition. Let f:2𝒩0 be a polymatroid over the ground set 𝒩. For a subset T of S, the strength-ratio of T in S, denoted φ(T|S), is the value φ(T|S):=|S||T|f(S)f(T), with the convention that x/0=+ for every x0.

Definition 15.

[12] Let f:2𝒩0 be a polymatroid over the ground set 𝒩. A strength decomposition of 𝒩 with respect to f is a sequence of sets S0S1Sw such that: (i) S0=𝒩, (ii) Sw= and Si for every i[0,w1], (iii) for every i[w], Si=argminSSi1φ(S|Si1), and (iv) the strength ratios φ(Si|Si1) are nondecreasing in i.

For every polymatroid f:2𝒩0 over a ground set 𝒩, there exists a strength decomposition of 𝒩 with respect to f 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 f:2𝒩0 be a polymatroid over the ground set 𝒩. Let 𝒩=S0S1Sw= be a strength decomposition with respect to f. Then, Si is a closed set for every i[0,w].

Approximating 𝒒𝒕 in Polynomial Time.

Now, we show how to approximate qt in polynomial time. For every timestep t[n], let 𝒩t=S0S1Sw= be a strength decomposition with respect to f|𝒩t. Let i[0,w] be the index such that etSi1Si. We define ηt:=|𝒩tSi|f(𝒩t)f(Si). We prove that the parameter ηt is an r-approximation of qt.

Lemma 17.

For every t[n], qtrηtqt.

Replacing the estimate qt in Algorithm 1 by ηt gives a randomized online algorithm in polynomial time with the same competitive ratio, which completes the proof of Theorem 1.

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 G=(V,E) where each hyperedge eE is a non-empty subset of V. We use n:=|V| and m:=|E| to denote the number of vertices and hyperedges in the hypergraph respectively. We note that E is a multi-set, i.e., it can contain multiple copies of the same hyperedge. We recall that the hypergraph G=(V,E) can equivalently be represented as a bipartite graph G=(VE,E), where we have a node for each vertex uV and each hyperedge eE and an edge between nodes uV and eE if the hyperedge e contains the vertex u. A hypergraph is connected if its bipartite representation is connected. A subhypergraph of G is a subset of hyperedges FE and it is said to be spanning if eFe=V and connected if (V,F) 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 V of nodes is known in advance while the hyperedges E arrive in an online fashion. For a color assignment C:E+, we say that a color c+ is a connected spanning color if the set of hyperedges with color c 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 O(log2n) that runs in polynomial time. We show that we can achieve the same competitive ratio via a much simpler estimator.

For every integer t[m], let et be the hyperedge arriving at time t. Define Et:={e1,,et} as the set of hyperedges that have arrived until time t. For a subset SV of nodes, let Et(S):={eEt:Se} be the set of hyperedges containing all vertices in S that have arrived until time t. We define

ηt:=min{u,v}et:uv|Et({u,v})|.

Equivalently, ηt 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 qt in Algorithm 1 by ηt gives a randomized online algorithm with the same competitive ratio. We note that |Et({u,v})| for every u,vV can be maintained in O(|et|2) time in the t-th timestep, which implies that the algorithm takes O(|et|2) time in the t-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 o(logn/loglogn) for Disj-Set-Cover, where n is the size of the ground set. Can we show that there is no randomized algorithm with expected impure competitive ratio o(logf(𝒩)) or pure competitive ratio o(log2f(𝒩)) for Disj-Polymatroid-Bases? Thirdly, we recall that our online model for Disj-Polymatroid-Bases assumes knowledge of f(𝒩). Is it possible to achieve polylogarithmic competitive ratio without knowledge of f(𝒩)? 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 h:2𝒱0, we recall that a set Q𝒱 is a quotient of h if h(e+(𝒱Q))>h(𝒱Q) eQ, i.e., if each element eQ has strictly positive marginal with respect to 𝒱Q. Since the definition of quotients might seem contrived, we interpret it for some concrete polymatroids:

  1. 1.

    Let G=(V,E) be a connected hypergraph (graph respectively). Consider the polymatroid h:2E0 that arises in Disj-Conn-Spanning-Subhypergraphs (Disj-Spanning-Trees respectively) with G being the input instance. Quotients of h correspond to union of cut-set of disjoint subsets of vertices, i.e., each quotient QE of h is of the form Q=S𝒞δ(S), where the family 𝒞2V{,V} 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 G.

  2. 2.

    Let G=(V,E) be a hypergraph. Consider the polymatroid h:2E0 that arises in Disj-Set-Cover with G being the input instance, i.e., h is the coverage function of G. Quotients of h correspond to union of vertex-isolating cuts of G, i.e., each quotient QE of h is of the form Q=uSδ(u) where SV. In particular, the minimum size of a non-empty quotient is equal to the min-degree of G.

Appendix B Min-size of a non-empty Quotient and 𝒌

Let q be the minimum size of non-empty quotients of f, i.e.,

q:=min{|Q|:Q𝒩 is a quotient of f}.

We recall that

k=k(f)=minA𝒩:f(A)<f(𝒩)e𝒩fA(e)f(𝒩)f(A).

The following lemma shows that q is an r-approximation of k.

Lemma 18.

qr1<kq.

Proof.

We first show the lower bound on k. Consider a closed set A𝒩 with f(A)<f(𝒩) and k=e𝒩fA(e)f(𝒩)f(A). Since A is closed, we have that fA(e)1 for every e𝒩A and 𝒩A is a quotient of f. Therefore, we have

k=e𝒩fA(e)f(𝒩)f(A) >e𝒩fA(e)f(𝒩)f(A)1(since x>x1)
e𝒩fA(e)r1(since f(𝒩)f(A)f(𝒩)r)
|𝒩A|r1(since fA(e)1 for every e𝒩A)
qr1.(since 𝒩A is a quotient of f)

We now show the upper bound on k. Consider a minimum sized non-empty quotient Q. Let A:=𝒩Q. Since A is a closed set and A𝒩, we have f(𝒩)>f(A). We note that for every eQ, fA(e)=f(A+e)f(A)f(𝒩)f(A), where the inequality is by monotonicity of the function f. Therefore,

q=|Q| eQ(fA(e)f(𝒩)f(A))(since fA(e)f(𝒩)f(A))
=e𝒩fA(e)f(𝒩)f(A)(since fA(e)=0 for every eA)
k.(by definition of k)

Appendix C Random Sampling gives a Base

In this section, we prove Lemma 6. A variant of the lemma was shown in [2]. See 6

Proof.

If k2logr, then we have p=1, which implies that S=𝒩 which is a base of f. Henceforth, we may assume that k>2logr and consequently, p=2logrk. We recall that n=|𝒩|. Let σ=(e1,e2,,en) be a uniformly random permutation of the elements in 𝒩. For every i[0,n], we define 𝒩i:={ej:j[i]} and Si:=S𝒩i. We note that S0=𝒩0= and fS0(𝒩)=r.

Let i[n]. We consider the distribution of element ei conditioned on Si1. Since σ is a uniformly random permutation, ei can be an arbitrary element in 𝒩Si1, of which each has the same probability. That is, for every element e𝒩,

𝐏𝐫σ,S[ei=e|Si1]={0,eSi11|𝒩||Si1|.e𝒩Si1

We note that element ei is picked with probability at least p. If element ei is picked, then Si=Si1+ei. Otherwise, Si=Si1. Hence, we have

𝔼σ,S[fSi1(𝒩)fSi(𝒩)|Si1] =e𝒩Si11|𝒩||Si1|𝔼σ,S[fSi1(𝒩)fSi(𝒩)|Si1andei=e]
e𝒩Si1p|𝒩||Si1|fSi1(e)
pne𝒩Si1fSi1(e)
pnkfSi1(𝒩),

where the last inequality is by the definition of k. This implies that

𝔼σ,S[fSi(𝒩)] =𝔼σ,S[fSi1(𝒩)𝔼σ,S[f𝒮i1(𝒩)f𝒮i(𝒩)|Si1]]
𝔼σ,S[(1pnk)fSi1(𝒩)]
=(1pnk)𝔼σ,S[fSi1(𝒩)].

By setting i=n, we have

𝔼σ,S[fSn(𝒩)] (1pnk)n𝔼σ,S[fS0(𝒩)]
=(1pnk)nr
<exp(pk)r=exp(2logr)r<12.

According to Markov’s inequality, we have

𝐏𝐫σ,S[fSn(𝒩)1]𝔼σ,S[fSn(𝒩)]<12,

which shows that S is a base with probability at least 12.