Abstract 1 General context 2 Preliminaries 3 Truncated Max-Sum Greedy 4 Integrality gaps for the configuration LP References

Submodular Max-Min Allocation Under Identical Valuations

Kimon Boehmer ORCID DIENS, École normale supérieure, PSL University, Paris, France
LIP, ENS Lyon, France
Abstract

In the problem of Submodular Max-Min Allocation, we are given a set of items, a set of players, and monotone submodular valuation functions that represent the satisfaction of a player with a certain subset of items. The goal is to find an allocation of the items to the players that maximizes the lowest satisfaction among all players.

We study this problem in the special case where all players have the same valuation function. We devise a greedy algorithm which gives a 0.4-approximation, improving the previously best factor of 10270.37 by Uziahu and Feige.

Furthermore, we study the integrality gap of the configuration LP when players have identical valuations. By constructing a variable assignment to the dual from a primal integral solution, we give the first constant upper bound on the integrality gap for submodular valuations. Generalizing the result to the case where players’ allocations must be independent in k given matroids, we derive a 𝒪(k)-estimation algorithm for max-min allocation subject to k matroid constraints under identical valuations.

Keywords and phrases:
Submodularity, Approximation algorithms, Allocation, Configuration LP
Copyright and License:
[Uncaptioned image] © Kimon Boehmer; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Submodular optimization and polymatroids
; Theory of computation Linear programming ; Theory of computation Packing and covering problems
Related Version:
Full Version: https://arxiv.org/abs/2604.12417
Acknowledgements:
I want to thank Chien-Chung Huang for introducing me to the topic and for many helpful dicussions. Also, I thank the anonymous reviewers for their valuable suggestions.
Editor:
Pierre Fraigniaud

1 General context

Item allocation is a fundamental problem in combinatorial optimization. Its task consists of distributing a set of items to a set of players, each equipped with a valuation function over subsets of items, with the goal of optimizing a specific objective function. The most well-studied objective functions are max-sum, i.e. maximizing the average satisfaction across all players, and max-min, maximizing the satisfaction of the least satisfied player. This work focuses on the latter objective function. Applications of max-min allocation range from assigning classrooms to charter schools [19] to sensor placement and scheduling [18].

Usually, research focuses on restricted classes of valuation functions that are both expressive and computationally tractable. For example, a common assumption is that each player has an additive valuation function; that is, each item is associated with a certain weight, and the satisfaction of a player is the sum of the weights of the received items. Another important and more general class of valuation functions are submodular functions; besides their useful algorithmic properties, they naturally model the law of diminishing returns in economics. This makes them especially interesting in the context of allocation.

In this work, we will restrict our attention to the case where all players share the same submodular valuation function. It is interesting to note that the strongest hardness results of approximating submodular max-min allocation to a factor of 11e already hold even for identical valuations [17, 21].

1.1 Our results

We design a simple greedy algorithm, which we call Truncated Max-Sum Greedy, that achieves a 0.4-approximation for identical submodular valuations, improving the previously best-known factor of 10270.37 by Uziahu and Feige [28]. Essentially, given a fixed threshold α, our algorithm greedily builds an allocation while ignoring players who already have items of value α𝖮𝖯𝖳. Our main result is the following:

Theorem 1.1.

There is a 25-approximation for submodular max-min allocation under identical valuations.

While showing a guarantee of 13 is rather straightforward, the improvement to 25 uses some graph-theoretic tools such as Hall’s theorem and certain arguments about the digraph of reallocations needed to transform the greedy solution into an optimal solution. Our algorithm also provides guarantees in more general settings. As an example, we consider the problem of submodular max-min allocation under identical valuations with an additional global cardinality constraint: The objective remains the same, but we are allowed to use at most k items in total. We show the following:

Proposition 1.2.

There is a 13-approximation for submodular max-min allocation under identical valuations with a global cardinality constraint.

Furthermore, we study the integrality gap of the configuration LP for max-min allocation, an important LP relaxation for max-min allocation.

We bound the integrality gap of the configuration LP by transforming certain primal integral solutions of the configuration LP to satisfying assignments of a dual formulation of this LP. Unlike to similar approaches when rounding configuration LPs [25], our analysis only uses the fact that the primal integral solution is at most 𝖮𝖯𝖳. As a result, our proof does not provide an algorithm (not even a non-trivial exponential one) for constructing an integral solution that demonstrates the integrality gap.

Using this technique, we give the first constant upper bound on the integrality gap of the configuration LP for identical submodular valuations.

Theorem 1.3.

The integrality gap of the configuration LP with identical submodular valuations is at most 3.

Additionally, we give a lower bound on this integrality gap.

Proposition 1.4.

The integrality gap of the configuration LP with identical submodular valuations is at least 4/3, even with 6 items and 3 players.

We believe that the technique of Theorem 1.3 can be used to show constant integrality gaps for submodular max-min allocation in more general settings. As an example, we also consider the case where we are given k matroids on the items and the sets that players receive must all be independent in all matroids. For this case, we also give a constant upper bound on the integrality gap, independent of k.

Theorem 1.5.

The integrality gap of the configuration LP with identical submodular valuations and any number of matroid constraints is at most 5.

Accounting for the loss incurred when solving the separation oracle of the dual LP, this yields a 𝒪(k)-estimation algorithm for max-min allocation subject to k matroid constraints.

1.2 Related work

Additive valuations.

Additivity is one of the simplest assumptions on the valuation functions of the players. In the general version where players may have non-identical additive valuation functions, the state-of-the-art algorithm by Chakrabarty, Chuzhoy and Khanna [7] achieves an approximation factor of 𝒪(nε). Interestingly, the best-known lower bound, shown by Bezáková and Dani [6], only rules out a factor better than 2.

A well-studied special case is the so-called restricted assignment case, where players have identical valuations, but a bipartite graph between items and players restricts the allowed allocations. To tackle this problem, Bansal and Sviridenko [4] introduced the configuration LP, and subsequent work [4, 2, 22, 1] led to the current best approximation factor of 14ε for the problem [10]. Recently, the integrality gap of the configuration LP for restricted assignment was shown to be at most 3.53 [15]. For identical valuations, Woeginger designed a polynomial-time approximation scheme (PTAS) [30].

Submodular valuations.

The first approximation algorithm for general submodular valuations was given in [14], and even though improvements have been made [13], the best-known approximation factor remains polynomial in the input. For a constant number of players, Chekuri, Vondrák and Zenklusen [8] gave a 11/e-approximation.

For max-sum, Vondrak gave a 11/e-approximation [29], which was shown to be optimal [17, 21]. In the restricted assignment case, Bamas, Morell and Rohwedder [3] generalized the techniques from Bansal and Sviridenko to obtain a 𝒪(loglogn)-approximation. Under identical valuations, Krause, Rajagopal, Gupta and Guestrin [18] gave the first constant-factor approximation algorithm, which uses an algorithm for max-sum as a subroutine. Plugging in Vondrak’s optimal algorithm, their method achieves a e13e0.21-approximation.

More recently, another line of research has studied the problem of finding an allocation which approximates the maximin share of all players, which is the largest value a player can ensure by proposing an allocation and receiving the bundle which is worst for them. Under identical valuations, this problem is equivalent to submodular max-min allocation. Barman and Krishnamurthy [5] gave a different algorithm with the same guarantee as Krause et al. [18]. Ghodsi, Hajiaghayi, Seddighin, Seddighin and Yami [12] improved this factor to 13. Finally, Uziahu and Feige [28] gave a 1027-approximation by analyzing a bidding procedure which uses a truncated version of the valuation function.

Machine scheduling.

If we reverse the roles of max and min, we obtain the classical makespan minimization setting, where we try to assign jobs to machines so as to minimize the latest completion time. As its max-min counterpart, it admits a PTAS for identical additive valuations [16]. The submodular variant of makespan minimization, sometimes called submodular load balancing, is not o(nlogn)-approximable [27]. This stands in stark contrast to its “dual” max-min variant.

2 Preliminaries

Consider a set J (with n:=|J|) of items and a valuation function f:2J that assigns a value to each subset of J. We simplify notation by writing f(j) instead of f({j}), for jJ.

We focus on the case where f is submodular: A function f:2J is called submodular, if for every two subsets S and T,

f(ST)f(S)+f(T)f(ST).

Closely related to submodular functions is the notion of marginal contribution. We define the marginal contribution on f of an item j with respect to a subset S as Δf(jS):=f(S{j})f(S). We extend this definition to sets (i.e. Δf(ST)=f(ST)f(T)). In fact, a useful characterization of submodularity is that for all subsets ST of J and jJ,

Δf(jS)Δf(jT).

A function is monotone if for all ST, we have f(S)f(T). We assume f()=0, and so a monotone function is non-negative. In the full version of the paper, we give some more basic facts about submodular functions that will be used implicitly throughout our proofs.

We define P (with m:=|P|) to be the set of players. Since we are focusing on identical valuations, each player p has the same valuation function f describing their satisfaction with each set of items.

More formally, a partial allocation of J is a partition A=(A1,,Am) of a subset of J. Slightly abusing notation, we write jA for jpPAp. If every item is allocated, we speak of an allocation. We refer to the sets of items in a (partial) allocation as configurations. We define f𝗌𝗎𝗆(A)=pPf(Ap) and f𝗆𝗂𝗇(A):=minpPf(Ap). We will refer to players p with f(Ap)=f𝗆𝗂𝗇(A) as min players.

We briefly define matroids since we discuss them in Section 4.3. A matroid =(J,) is a non-empty subset-closed set system in which for any JJ, all sets of restricted to J which are maximal have the same cardinality. The sets in are called independent.

Submodular Max-Min Allocation

With these definitions, we are ready to define the main computational problem of interest.

Submodular Max-Min Allocation
Input: A set of items J, a set of players P={1,,m}, and a monotone submodular function f:2J.

Output: An allocation A that maximizes f𝗆𝗂𝗇(A)

Notice that each player has the same valuation function f. We assume that we are given access to a so-called value oracle, which, given any SJ, returns f(S).

Approximation Hardness.

Khot, Lipton, Markakis and Mehta [17] showed that when players have identical valuations, it is NP-hard to distinguish between the case where all players can achieve value 1 and the case where a player achieves at most value 11e+ε on average. Since the average value is always at least as big as the min value, their hardness result, while originally designed for max-sum, immediately translates to our setting:

Corollary 2.1 ([17]).

Unless P=NP, Submodular Max-Min Allocation cannot be approximated within a factor 11e+ε for any ε>0.

3 Truncated Max-Sum Greedy

In this section, we present our greedy algorithm for Submodular Max-Min Allocation.

Fix a threshold α1. We assume that the algorithm knows the optimal value 𝖮𝖯𝖳. 111We assume that the size of a rational representation of f(S) is at most polynomial in n+m. We can therefore scale f such that all values are integral and perform a binary search with start interval [0,f¯(J)], where f¯ is the scaled version of f. Under our assumption, the binary search procedure will only recurse a polynomial number of times. First, the algorithm takes the m items j with largest f(j)-value and allocates one of them to each player. Then, it greedily chooses the pair of available item j and player p that maximizes the marginal contribution Δf(jAp), under the constraint that f(Ap)<α𝖮𝖯𝖳, i.e. we only consider players that have not reached an α-fraction of their desired value.

The goal is to show that each player will eventually surpass this threshold, which proves an α-approximation. The pseudocode is shown in Algorithm 1.

Algorithm 1 The truncated max-sum greedy.
Observation 3.1.

We can assume f(j)<α𝖮𝖯𝖳 for all jJ.

Proof.

We show that if the truncated max-sum greedy gives an α-approximation on instances where for all items jJ, f(j)<α𝖮𝖯𝖳, then it also gives an α-approximation for the general case.

Let 𝖮𝖯𝖳 denote the optimal value in the general instance. Assume that there are k big items j with f(j)α𝖮𝖯𝖳. In an optimal allocation, these items can be given to at most k players. Thus, there exists an instance with max(0,mk) players and a subset J{jJf(j)<α𝖮𝖯𝖳} of items that achieve at least 𝖮𝖯𝖳.

If greedy is performed on the instance with max(0,mk) players and without the big items, then it will obtain a solution of value at least α𝖮𝖯𝖳 (since the optimal solution on this reduced instance has value 𝖮𝖯𝖳, as argued above). Adding k players, each containing one big item, gives a solution of value at least α𝖮𝖯𝖳 to the general instance.

As a warm-up, we first show that we obtain a 13-approximation when setting α=13. Afterwards, we will show how to improve the factor to 25 using a more refined analysis.

Consider a greedy allocation A and an optimal allocation A. Let q be a min player in A, i.e. f(Aq)=f𝗆𝗂𝗇(A). We may assume nm, since otherwise the optimal value is 0. We first show the following simple but crucial lemma.

Lemma 3.2.

Assume f(Aq)<α𝖮𝖯𝖳. Then,

jAΔf(jAq)pP{q}f(Ap)

Proof.

For any item j, let A(j) denote the partial allocation of the algorithm immediately before j is allocated. We have jAqΔf(jAq)=0. For all other pP{q}, we have

jApΔf(jAq)jApΔf(jAq(j))jApΔf(jAp(j))=f(Ap) (1)

The first inequality follows from the submodularity of f, since Aq(j)Aq. For the second inequality, notice that the algorithm always chooses to allocate an item to the player p with maximal marginal contribution, as long as f(Ap)<α𝖮𝖯𝖳. Therefore,

jAΔf(jAq)=pP{q}jApΔf(jAq)+jAqΔf(jAq)pP{q}f(Ap). (2)

We also need an upper bound on f(Ap). The following lemma is sufficient for the 13-approximation, but later we will need a stronger upper bound for the 25-approximation.

Lemma 3.3.

For all players pP, f(Ap)2α𝖮𝖯𝖳.

Proof.

Player p can only be allocated one item j that surpasses α𝖮𝖯𝖳. This item has size α𝖮𝖯𝖳, by Observation 3.1. Let A be the partial allocation before j is allocated to p. Then,

f(Ap)=f(Ap)+Δf(jAp)α𝖮𝖯𝖳+f(j)2α𝖮𝖯𝖳

We can now prove the approximation ratio of 13 by the following lemma.

Lemma 3.4.

f𝗆𝗂𝗇(A)<f(Aq)+1mpPf(Ap).

Proof.

f𝗆𝗂𝗇(A) minpPf(ApAq) (3)
minpPf(Aq)+jApΔf(jAq)
f(Aq)+1mpPjApΔf(jAq)
=f(Aq)+1mjAΔf(jAq)
f(Aq)+1mpP{q}f(Ap) (4)
<f(Aq)+1mpPf(Ap) (5)

(3) uses monotonicity of f, (4) uses Lemma 3.2 and (5) holds since f(Aq)>0 by our assumption nm and the fact that greedy starts by allocating one item to each player. To see that when setting α=13, this indeed implies the approximation ratio of 13, substitute 1mpPf(Ap)2α𝖮𝖯𝖳 by Lemma 3.3. Taking 2α𝖮𝖯𝖳=23𝖮𝖯𝖳 on the other side, we obtain 13𝖮𝖯𝖳<f(Aq).

 Remark 3.5.

This algorithm extends to the case where we are additionally given a global cardinality constraint. That is, we must choose a subset of at most k items to allocate.

To see this, consider again an optimal allocation A and the greedy allocation A. Notice that in our new setting, we might have pPAppPAp and therefore the first equality in (2) no longer holds. However, we can circumvent this issue as follows. We can assume that both A and A contain exactly k items by monotonicity. Let b:pPAppPAp be any bijection such that b(j)=j for all j(pPAp)(pPAp). For all jA, let A(j) be the partial greedy allocation exactly before allocating j. Consider any jA and let b(j)Ap. We have

Δf(jAq)Δf(jAq(b(j)))Δf(b(j)Ap(b(j))) (6)

using that when b(j) is allocated, item j is still available and thus greedy must prefer allocating b(j) to p than allocating j to q. By (6) and the fact that b is a bijection,

jAΔf(jAq) pP{q}jA:b(j)ApΔf(b(j)Ap(b(j)))
=pP{q}jApΔf(jAp(j))
=pP{q}f(Ap)

reproving Lemma 3.2 in this more general setting. The rest of the analysis can be done as for the unconstrained case. Therefore, Proposition 1.2 holds.

Ratio of 2/5

We now prove that if α=25, then f𝗆𝗂𝗇(A)25𝖮𝖯𝖳. To this end, suppose for contradiction that there exists a counterexample to this claim. We assume that this instance (P,J,f) is minimal with respect to |P| and that f is scaled in such a way that f𝗆𝗂𝗇(A)=𝖮𝖯𝖳=5. We also fix a min player q with f(Aq)2. Note that we could have assumed f(Aq)<2 but this is not necessary for the proof.

Let us define 𝖺𝗏𝗀:=1mpPf(Ap). From Lemma 3.4 we obtain 𝖺𝗏𝗀>3. Our goal is to show that 𝖺𝗏𝗀3, thereby arriving at a contradiction.

Letting api denote the i-th item allocated to p, define β:=maxpP:|Ap|2Δf(ap2{ap1}) as the largest marginal contribution of an item during the greedy algorithm which is not the first item of a player. Notice that it suffices to consider ap2, since the marginal contributions of later items of a player can only be smaller. If β1, then no player can obtain more than value 3, so we can assume β>1. Also, β2 holds by Observation 3.1. Notice finally that for pP, we have f(ap1)β, as otherwise the maximizer ap2 in the definition of β would have been among the largest m items, and could not be the second item of a player.

Lemma 3.6.

For each set of k1 players p1,,pkP{q}, there is no set of k players r1,,rkP such that i=1kApii=1kAri.

Proof.

Removing k players and all items in i=1kApi from the instance leads to a smaller counter-example: The greedy algorithm still behaves identically on all other players, including q, so there is a valid greedy allocation A on the instance with f𝗆𝗂𝗇(A)2. On the other hand, since all removed items only belonged to k players in the optimal solution, removing these k players gives a solution A with f(Ap)5 for all remaining players p.

The above lemma allows us to prove the existence of a “fixed” item for each player:

Lemma 3.7.

There is a permutation ϕ:PP such that ApAϕ(p) for all pP.

Proof.

Consider the bipartite graph G=(PP,E) with the set of players P on the left side and a copy P of the set of players on the right side, and edges (p,p) if ApAp. Let N(Q)=qQ{vV(q,v)E} be the neighbourhood of Q. Consider an arbitrary non-empty subset QP{q}. By Lemma 3.6, |N(Q)||Q|+1. For the set {q}, we have |N({q})|1 since Aq and each item in q must be allocated somewhere in the optimal solution. For the set , we have |N()|0=||. For any other set {q}QP, we know that |N(Q)||N(Q{q})||Q{q}|+1=|Q|. Altogether, we obtain that for any subset QP, |N(Q)||Q|. We can hence apply Hall’s Marriage Theorem to conclude that G contains a perfect matching. This matching defines exactly the desired permutation ϕ. To simplify notation, we consider an optimal solution such that ϕ is the identity function, by just permuting the allocations of the players in the optimal solution. We now define the reallocation digraph G=(P,E): Each pair of players (p,r)P2 is associated to a set Apr=ApAr. An arc (p,r) exists if Apr. By Lemma 3.7, each vertex has a self-loop.

We subdivide the players of our instance into two categories P=BS. The set B of big players consists of all players b such that f(Ab)>4β. The set S of small players consists of all players p with f(Ap)4β.

Notice that qS. This partition is mainly motivated by the following observation:

Observation 3.8.

For any pB, it holds that |Ap|=2.

Proof.

Any player pB cannot be allocated just one item j since f(j)<24βf(Ap). If p was allocated k items for k>2, then f({ap1,,apk1})<2 must hold in order for the greedy algorithm to allocate apk to p. But since f(ap1)β, we must have Δf(apk{ap1,,apk1})Δf(ap2{ap1})2β. So we have

f(Ap)=f(Apapk)+Δf(apk{ap1,,apk1})<2+2βf(Ap),

contradicting our assumption k>2.

For each pB, denote by xp the item in App and denote by yp the other item of Ap. Lemma 3.7 guarantees the existence of the first item and Observation 3.8 implies that there is exactly one other item. Lemma 3.6 guarantees that the other item is not in App.

Figure 1: Visualization of Lemma 3.9.
Lemma 3.9.

For any b,cB with bc, (b,c)E.

Proof.

Suppose for contradiction that Abc, where ybAbc. Let pP be the player such that ycAp. Consider a modified instance where we remove player c as well as items {yb,xc,yc} and add an item z which contributes an additive factor of t:=max(f(yb),Δf(ycAp{yc})) to every set. This function is submodular and monotone. Notice that this addition of z does not change any marginal contribution, i.e. for jz we have Δf(jS)=Δf(jS{z}).

The greedy algorithm on the new instance may allocate all items that give marginal contribution larger than t in the original run as before. As soon as all remaining items give marginal contribution t, it may choose to allocate item z to player b (since z contributes the same to each player, and greedy considers b as f(Ab)=f()<2 or f(Ab)=f(xb)<2).

Since bB, we have f(Ab)>4β2. Furthermore, f({xb,z})f(xb)+f(yb)f(Ab)2. Consider the moment in the greedy algorithm after z was allocated to b. If xb had already been allocated, then b has reached the threshold of 2 and will be thus ignored by the greedy algorithm. Therefore, all following items can be allocated to the same players as in the original run of the greedy algorithm. Now suppose xb had not yet been allocated. Since the addition of z to player b does not influence the marginal contributions and f(z)<2, the greedy algorithm can allocate all items up to xb in the same way (including xb). After allocating xb, player b has reached the threshold of 2 and thus again the greedy algorithm can allocate in the same way as in the original run. Thus, we again obtain f𝗆𝗂𝗇(Aq)2, for the new allocation A.

On the other hand, we can modify our old optimal allocation, by just replacing item yc by item z. Call the new optimal allocation A. We have

f(Ap) =f(Ap{z})+t
=f(Ap{yc})+t
f(Ap{yc})+Δf(ycAp{yc})
=f(Ap)5

We obtained a counterexample with less players, a contradiction. Thus, in fact, G[B] is an independent set. We now argue that B cannot be too large. To this end, we first prove the following two helpful lemmata. For a player p, let us define Vp:=(r,p)EArp the set of items reallocated to p and Wp:=(p,r)EApr the set of items reallocated from p.

Lemma 3.10.

For any player p, jWpΔf(jAq)<2.

Proof.

Recall that A(j) denotes the partial greedy allocation exactly before item j is allocated. Let be the last item allocated to Ap. By Lemma 3.7, there exists an item kApp. Therefore, WpAp{k} and we can write

jWpΔf(jAq) jAp{k}Δf(jAq)
jAp{k}Δf(jAq(j))
Δf(Ap())+jAp{k,}Δf(jAp(j))
Δf(Ap(k))+jAp{k,}Δf(jAp(j))
Δf(kAp(k))+jAp{k,}Δf(jAp(j))
=f(Ap{})<2

The third inequality is because the items were allocated to p and not q. The fourth inequality comes from Ap()Ap(k), then we use that at the partial allocation A(k), greedy prefered k over . Finally, the strict inequality holds because otherwise p would have achieved the α𝖮𝖯𝖳 threshold and not received any more items.

Lemma 3.11.

For any pP, jVpΔf(jAq)3+jWpΔf(jAq)f(Ap).

Proof.

First, notice that

jApΔf(jAq)Δf(ApAq)f(Ap)f(Aq)52=3. (7)

Since Vp=Ap(ApWp), ApWpAp, and WpAp, we can write

jVpΔf(jAq)=jApΔf(jAq)+jWpΔf(jAq)jApΔf(jAq)

Now, the result follows by using (7) and (1).

Lemma 3.12.

|B||S|.

Proof.

Let S0 be the set of small players that do not receive any item from big players, let S1 be the set of small players that receive one item from one big player and let S2 be the set of small players that receive items from multiple big players. If we can show S2=, this means that no two big players have an edge towards the same small player. But since every big player has an edge towards a small player by Lemma 3.9, by pigeonhole principle this would imply the desired |B||S|.

Associate with each player p the quantity

g(p):=jWpΔf(jAq)jVpΔf(jAq)

Note that pPg(p)=0, as pPWp=pPVp. We can decompose the set of players as

P=S0pS1{p,b1(p)}pS2{p,b1(p),b2(p),,bk(p)}

where bi(p) denotes the i-th big player that gives an item to p, in an arbitrary order. Therefore,

pPg(p)=pS0g(p)+sS1g(s)+g(b1(s))+sS2g(s)+g(b1(s))++g(bk(s)) (8)

Since pPg(p)=0, to prove S2=, it is enough to show that all terms of (8) are 0, and that terms of the last sum are <0. This is the purpose of the rest of the proof.

  1. 1.

    g(p)0 for pS0. This follows directly from Lemma 3.11 and f(Ap)4β3.

  2. 2.

    g(s)+g(b)0 for sS1 and b=b1(s). Since s is a small player, we have f(As)4β, and by Lemma 3.11, we get g(s)1β. If we can show g(b)β1, then this implies g(s)+g(b)0. Recall that Wb={yb}.

    As the second item allocated to b contributes at most β, we have Δf(xbAq)+Δf(ybAq)2+β. By the facts f(Ab)f(Aq)3 and Ab=Vb(AbWb), we get

    jVbΔf(jAq) Δf(Vb(AbWb)(AbWb)Aq)
    =Δf(Ab(AbWb)Aq)
    f(Ab)Δf(AbWbAq)f(Aq)
    3Δf(xbAq)

    In total, g(b)Δf(ybAq)jVbΔf(jAq)Δf(ybAq)3+Δf(xbAq)β1.

  3. 3.

    g(s)+g(b)+g(c)<0 for sS2, b=b1(s) and c=b2(s). We have f((AbWb)Aq)f(AbWb)+f(Aq)4, but since f(Ab)5, we must have

    jVbΔf(jAq)Δf(Vb(AbWb)(AbWb)Aq)1.

    The same holds for player c. By Lemma 3.10, we also know g(s)<2Δf(ybAq)Δf(ycAq). Therefore

    g(b)+g(c)+g(s)Δf(ybAq)1+Δf(ycAq)1+g(s)<0.

We can conclude as follows:

𝖺𝗏𝗀 =1mpPf(Ap)
1m[|B|(2+β)+|S|(4β)]
1m[|B|+|S|2(2+β)+|B|+|S|2(4β)]
=1m[6(|B|+|S|)2]
=3

where the last inequality is because 2+β4β, and |B||S| by Lemma 3.12. This contradicts 𝖺𝗏𝗀>3, as stated in the beginning. This finishes the proof of Theorem 1.1.

We discuss the (close-to-) tightness of this analysis in the full paper version. In particular, we show that no value of α can lead to an approximation factor better than 2+S4+3S0.4068, where S=i=11si1 with sk=1+i=1nsi is a constant related to the Sylvester sequence or Salzer sequence and was already used in [20, 24] for other packing problems.

4 Integrality gaps for the configuration LP

Assignment LP (A)
i=1mxij =1 j=1,,n
j=1nfi(j)xij T i=1,,m
xij 0 i,j
 
Configuration LP (P)
C𝒞(T,p)xp,C 1 pP
CjpPxp,C 1 jJ
xp,C 0 pP,C𝒞(T,p)
Figure 2: The assignment LP and the configuration LP.
Dual LP (D0)
maxpPzp jJyj s.t.
jCyj zp p,C
yj,zp 0 j,p
 
Dual LP (D1)
pPzp >jJyj
jCyj zp p,C
yj,zp 0 j,p
 
Dual LP (D2)
jJyj <m
jCyj 1 C
yj 0 j
Figure 3: The duals to the configuration LP.

In this section, we study the integrality gap of the configuration LP in the case of identical valuation functions. Before defining this LP, let us consider the simpler assignment LP: It searches for the largest T such that the program depicted in figure 2 on the left, where fi(j)=min(fi(j),T), has a feasible solution. Let us call this LP (A). It is a folklore result that (A) has an integrality gap of 2 when all players have identical additive valuations.

In the restricted assignment case, (A) has a bad integrality gap [4]. This was the initial motivation to study a stronger LP, namely the aforementioned configuration LP. Let 𝒞(T,p) contain the configurations C such that fp(C)T. We have a variable for each pair of player and configuration. The LP is shown in figure 2 on the right. The largest T such that this LP has a feasible solution is an upper bound on the optimal solution. We call this system (P).

It is well-known that the integrality gap of (P) in the general case is Ω(m) [4], while in the restricted assignment case it was shown to be smaller than 3.534 [15]. When valuation functions are submodular, the gap in the restricted assignment case is 𝒪(loglogn) [3].

(P) has exponentially many variables in general. Assuming that the primal minimizes the 0-vector, we can formulate its dual (D0), depicted in figure 3 on the left. The separation problem for this LP is the classical Knapsack problem. Since Knapsack admits a (F)PTAS [23], we can solve the LP up to any desired accuracy using the Ellipsoid method.

4.1 Constructing Dual Solutions

Notice that in (D0), a trivial solution always exists by setting all variables to 0. But if we had a solution with objective value larger than 0, then we could obtain an arbitrarily high objective value by scaling all variables. By LP duality, we get that if (D0) is unbounded, then (P) with the same threshold T is infeasible. Additionally, we can directly express the requirement that the objective function must be greater than 0 as a constraint. We call the resulting LP (D1), and it is shown in the middle of figure 3.

In the case of identical valuations, we can, without loss of generality, set all player variables to 1, and obtain (D2), as seen on the right of figure 3. This LP searches for a fractional hitting set of size smaller than m that hits all configurations larger than our threshold value. In (D2), we call the first constraint value constraint and the second set of constraints configuration constraints. If (D2) has a feasible solution, then (D0) is unbounded, and by weak duality, (P) is infeasible. We obtain the following:

Proposition 4.1.

Consider a class Γ of max-min allocation instances, and let 𝖮𝖯𝖳() be the optimal integral value of instance . If for all Γ, (D2) is feasible for and threshold value α𝖮𝖯𝖳(), then the integrality gap of the configuration LP on Γ is at most α.

We give one more useful observation:

Observation 4.2.

If there are variables yj such that the configuration constraints of (D2) are satisfied, jJyjm and there exists jJ such that for all Cj, we have jCyj>1, then (D2) is feasible.

Proof.

Since yj is not involved in any tight constraint, we can lower its value by a sufficiently small amount, satisfying the value constraint while keeping all configuration constraints satisfied.

4.2 Submodular valuations

The system (P) can also be formulated for submodular valuations. While in the additive case, solving the separation problem only incurred an ε error, this is different in the submodular case. This is because we need to solve a submodular version of knapsack for the separation oracle. Due to the hardness of maximum coverage [11], which includes submodular knapsack, we lose a factor of 11e while solving the LP. On the algorithmic side, this hardness is matched by a 11e-approximation by Sviridenko [26]. Therefore, any rounding algorithm that loses a factor of α translates into an approximation algorithm with factor (11e)α.

In the following, we show that apart from our combinatorial approach presented in the last section, an LP rounding approach could also lead to a constant-factor approximation.. See 1.3 We first show that we can restrict our attention to items j with f(j)𝖮𝖯𝖳. Let 𝖮𝖯𝖳() denote the optimal integral value for an instance , and let 𝖮𝖯𝖳𝖫𝖯() denote the optimal value of (P) on .

Lemma 4.3.

Let be an instance with an item j such that f(j)>𝖮𝖯𝖳(). Let be the instance obtained from by removing j and one player. Then, 𝖮𝖯𝖳()𝖮𝖯𝖳() and 𝖮𝖯𝖳𝖫𝖯()𝖮𝖯𝖳𝖫𝖯().

Proof.

Fix an arbitrary player pP. For an optimal solution x to (P), we have

CjpPxp,C1C𝒞(𝖮𝖯𝖳𝖫𝖯(),p)xp,C

Therefore, we can replace all configurations in which j appears by the conifigurations of player p, without violating any constraint except the constraint of player p. Since item j is not used anymore, this is a valid solution to (P) on instance of value 𝖮𝖯𝖳𝖫𝖯().

Now let A be an optimal integral allocation of and assume for contradiction f𝗆𝗂𝗇(A)>𝖮𝖯𝖳(). Adding a new player r with Ar={j} gives an integral allocation for of value min(fmin(A),f(j))>𝖮𝖯𝖳(), a contradiction. Therefore, if our instance contains items that contribute more than 𝖮𝖯𝖳, then we can instead consider a smaller instance with larger integrality gap. We now prove Theorem 1.3.

Proof.

Consider an arbitrary instance of our problem and scale the valuation function so that 𝖮𝖯𝖳()=1. We can assume that for all jJ, f(j)1 by Lemma 4.3. Let f¯(S):=min(2,f(S)), and consider an allocation A that maximizes pPf¯(Ap).

For jAp and f, define Δf(jAp):=Δf(ApAp{j}). For any pP, let Tf(p):=jApΔf(jAp). For each pP and jAp, we set the variables in (D2) to

yj=Δf¯(jAp)Tf¯(p)

if Tf¯(p)0, and yj=0 otherwise. For notational convenience, define y(S)=jSyj. We prove that y is a (strictly) feasible solution to (D2) when T=3. Suppose for contradiction that there is a set C of items with f(C)>3 but y(C)1.

Let q be the min player in A. Clearly, f(Aq)1, and thus Δf(CAq)>2. Let j=argmaxjCΔf(jAq)yj. If Δf(jAq)yj2, then we would have

Δf(CAq)y(C)jCΔf(jAq)y(C)jC2yjy(C)=2

which contradicts our assumption y(C)1. Therefore, we can assume Δf(jAq)yj>2.

We want to argue that reallocating item j to player q increases our potential according to f¯. Notice that the gain of adding j to player q is Δf¯(jAq), while the loss of removing j from player p is Δf¯(jAp)=yjTf¯(p). We now bound Tf¯(p). First, we have

Tf¯(p) =jApΔf¯(jAp)
=jApmin(0,Δf(jAp)max(0,f(Ap)2))
=jAp:Δf(jAp)f(Ap)2Δf(jAp)max(0,f(Ap)2)

If this sum is empty, then Tf¯(p)=0. Otherwise, we can assume that max(0,f(Ap)2) is subtracted at least once. Thus,

Tf¯(p) (jApΔf(jAp))max(0,f(Ap)2)
f(Ap)max(0,f(Ap)2)
2

where the second inequality is by submodularity of f. Therefore,

Δf¯(jAq)Δf¯(jAp)>Δf¯(jAq)2yj=Δf(jAq)2yj>0,

where the equality is because of f(Aq{j})=f¯(Aq{j}), which follows from our assumption f(j)1 and f(Aq)1. So, reallocating j from p to q increases the potential f¯, contradicting the optimality of A.

We conclude that the configuration constraints in (D2) are satisfied and not tight for T>3. We can finish the proof with Observation 4.2, since

jJyj=pPjApyjm.

A simple lower bound on the integrality gap is proven in the full version of the paper. See 1.4

4.3 Submodular valuations with matroid constraints

Now consider the following problem: Apart from items J, players P and a submodular monotone valuation function f, we are given k matroids 1=(J,1),,k=(J,k) and our task is to find a (partial) allocation A such that Api for all pP and all i[k]. The configuration LP can be formulated in the same way as before, and the separation oracle incurs a loss of 2.7k by using the algorithm of [9]. In the following, we prove Theorem 1.5, which states that the integrality gap is still bounded by 5 in this case and implies a 13.5k-estimation algorithm for max-min allocation subject to k matroid constraints.

Theorem 1.5. [Restated, see original statement.]

The integrality gap of the configuration LP with identical submodular valuations and any number of matroid constraints is at most 5.

Proof.

As before, let 𝖮𝖯𝖳=1, f¯(S):=min(2,f(S)) and define A and the yj-variables in the same way. The proof of Lemma 4.3 directly translates to the matroid constraint setting, so we may again assume f(j)1 for all jJ. Unlike before, we set yj=0 for all jAq, where q is a fixed min player in A. We also set yj=0 for all unallocated jA.

Now suppose for contradiction that there is a set Di=1kk of items with f(D)>5 but y(D)<1. We can partition D into D1D2 such that f(Di)2 for both i[2], by greedily adding items to D1 until f(D1)>2. By the assumption that f(j)1 for all jJ and submodularity, this implies f(D2)=f(DD1)>2. Now, choose the partition Di with the smaller y-value, and call it C. As a subset of D, clearly C is independent in all matroids, and y(C)<12.

Consider reallocating all items from C to player q and discarding all original items in Aq. We want to argue that this modified allocation has higher potential. The gain for player q is f¯(C)f¯(Aq)21=1. The loss for all other players is

pP{q}Δf¯(CAp) pP{q}jApΔf¯(jAp)
=pP{q}jApyjTf¯(p)
y(C)2
<1

since we can bound Tf¯(p) as in the proof of Theorem 1.3. The total potential change is negative. This contradicts our choice of A. Hence, the configuration constraints are all satisfied. The value constraint is satisfied since

jJyj=jAyj+jAqyj+pP{q}jApyj=pP{q}jApyjm1<m.

Notice that we could replace matroids by any set system closed under subsets, but then of course we might not be able to solve the configuration LP efficiently.

Conclusion

We presented an improved approximation algorithm for submodular max-min allocation under identical valuations, and showed some constant upper bounds on the integrality gap of the configuration LP under identical valuations.

However, several natural questions remain open. First, the proposed 0.4-approximation is still far from the best-known complexity-theoretical hardness bound of 11e0.63. Another promising direction is to further explore the configuration LP. It would be valuable to obtain constructive, polynomial-time rounding algorithms that achieve the proven integrality gaps. Furthermore, it would be interesting to see by how much we can relax the assumption of identical valuations. A natural first step in this direction would be to determine the integrality gap of the configuration LP in the case of submodular “related machines”, where each player can be satisfied with a different value.

References

  • [1] Chidambaram Annamalai, Christos Kalaitzis, and Ola Svensson. Combinatorial algorithm for restricted max-min fair allocation. ACM Transactions on Algorithms (TALG), 13(3):1–28, 2017. doi:10.1145/3070694.
  • [2] Arash Asadpour, Uriel Feige, and Amin Saberi. Santa claus meets hypergraph matchings. ACM Transactions on Algorithms (TALG), 8(3):1–9, 2012. doi:10.1145/2229163.2229168.
  • [3] Étienne Bamas, Sarah Morell, and Lars Rohwedder. The submodular santa claus problem. Proceedings of the 2025 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 616–640, 2025.
  • [4] Nikhil Bansal and Maxim Sviridenko. The santa claus problem. Proceedings of the thirty-eighth annual ACM symposium on Theory of computing, pages 31–40, 2006. doi:10.1145/1132516.1132522.
  • [5] Siddharth Barman and Sanath Kumar Krishnamurthy. Approximation algorithms for maximin fair division. ACM Transactions on Economics and Computation (TEAC), 8(1):1–28, 2020. doi:10.1145/3381525.
  • [6] Ivona Bezáková and Varsha Dani. Allocating indivisible goods. ACM SIGecom Exchanges, 5(3):11–18, 2005. doi:10.1145/1120680.1120683.
  • [7] Deeparnab Chakrabarty, Julia Chuzhoy, and Sanjeev Khanna. On allocating goods to maximize fairness. 2009 50th Annual IEEE Symposium on Foundations of Computer Science, pages 107–116, 2009.
  • [8] Chandra Chekuri, Jan Vondrák, and Rico Zenklusen. Dependent randomized rounding via exchange properties of combinatorial structures. 2010 IEEE 51st Annual Symposium on Foundations of Computer Science, pages 575–584, 2010. doi:10.1109/FOCS.2010.60.
  • [9] Chandra Chekuri, Jan Vondrák, and Rico Zenklusen. Submodular function maximization via the multilinear relaxation and contention resolution schemes. In Proceedings of the forty-third annual ACM symposium on Theory of computing, pages 783–792, 2011. doi:10.1145/1993636.1993740.
  • [10] Siu-Wing Cheng and Yuchen Mao. Restricted max-min allocation: Approximation and integrality gap. arXiv preprint arXiv:1905.06084, 2019. arXiv:1905.06084.
  • [11] Uriel Feige. A threshold of ln n for approximating set cover. Journal of the ACM (JACM), 45(4):634–652, 1998. doi:10.1145/285055.285059.
  • [12] Mohammad Ghodsi, MohammadTaghi HajiAghayi, Masoud Seddighin, Saeed Seddighin, and Hadi Yami. Fair allocation of indivisible goods: Improvements and generalizations. In Proceedings of the 2018 ACM Conference on Economics and Computation, pages 539–556, 2018.
  • [13] Michel X Goemans, Nicholas JA Harvey, Satoru Iwata, and Vahab Mirrokni. Approximating submodular functions everywhere. Proceedings of the twentieth annual ACM-SIAM symposium on Discrete algorithms, pages 535–544, 2009.
  • [14] Daniel Golovin. Max-min fair allocation of indivisible goods. School of Computer Science, Carnegie Mellon University, (2005).
  • [15] Penny Haxell and Tibor Szabó. Improved integrality gap in max–min allocation, or, topology at the north pole. Combinatorica, 45(2):1–38, 2025.
  • [16] Dorit S Hochbaum and David B Shmoys. Using dual approximation algorithms for scheduling problems theoretical and practical results. Journal of the ACM (JACM), 34(1):144–162, 1987. doi:10.1145/7531.7535.
  • [17] Subhash Khot, Richard J Lipton, Evangelos Markakis, and Aranyak Mehta. Inapproximability results for combinatorial auctions with submodular utility functions. Algorithmica, 52:3–18, 2008. doi:10.1007/S00453-007-9105-7.
  • [18] Andreas Krause, Ram Rajagopal, Anupam Gupta, and Carlos Guestrin. Simultaneous placement and scheduling of sensors. 2009 International Conference on Information Processing in Sensor Networks, pages 181–192, 2009. URL: https://ieeexplore.ieee.org/document/5211932/.
  • [19] David Kurokawa, Ariel D Procaccia, and Nisarg Shah. Leximin allocations in the real world. ACM Transactions on Economics and Computation (TEAC), 6(3-4):1–24, 2018. doi:10.1145/3274641.
  • [20] Frank M Liang. A lower bound for on-line bin packing. Information processing letters, 10(2):76–79, 1980. doi:10.1016/S0020-0190(80)90077-0.
  • [21] Vahab Mirrokni, Michael Schapira, and Jan Vondrák. Tight information-theoretic lower bounds for welfare maximization in combinatorial auctions. Proceedings of the 9th ACM conference on Electronic commerce, pages 70–77, 2008. doi:10.1145/1386790.1386805.
  • [22] Lukáš Poláček and Ola Svensson. Quasi-polynomial local search for restricted max-min fair allocation. ACM Transactions on Algorithms (TALG), 12(2):1–13, 2015. doi:10.1145/2818695.
  • [23] Sartaj Sahni. Approximate algorithms for the 0/1 knapsack problem. Journal of the ACM (JACM), 22(1):115–124, 1975. doi:10.1145/321864.321873.
  • [24] Steven S Seiden and Gerhard J Woeginger. The two-dimensional cutting stock problem revisited. Mathematical Programming, 102(3):519–530, 2005. doi:10.1007/S10107-004-0548-1.
  • [25] Ola Svensson. Santa claus schedules jobs on unrelated machines. In Proceedings of the forty-third annual ACM symposium on Theory of computing, pages 617–626, 2011. doi:10.1145/1993636.1993718.
  • [26] Maxim Sviridenko. A note on maximizing a submodular set function subject to a knapsack constraint. Operations Research Letters, 32(1):41–43, 2004. doi:10.1016/S0167-6377(03)00062-2.
  • [27] Zoya Svitkina and Lisa Fleischer. Submodular approximation: Sampling-based algorithms and lower bounds. SIAM Journal on Computing, 40(6):1715–1737, 2011. doi:10.1137/100783352.
  • [28] Gilad Ben Uziahu and Uriel Feige. On fair allocation of indivisible goods to submodular agents. arXiv preprint arXiv:2303.12444, 2023. doi:10.48550/arXiv.2303.12444.
  • [29] Jan Vondrák. Optimal approximation for the submodular welfare problem in the value oracle model. Proceedings of the fortieth annual ACM symposium on Theory of computing, pages 67–74, 2008. doi:10.1145/1374376.1374389.
  • [30] Gerhard J Woeginger. A polynomial-time approximation scheme for maximizing the minimum machine completion time. Operations Research Letters, 20(4):149–154, 1997. doi:10.1016/S0167-6377(96)00055-7.