alphaFeedback

WorksComputer Science

Nearly Tight Bounds for Proportional Group Fair Divisions and One-Sided Discrepancy

AI contributions
Not recorded for this paper
View details

Abstract

This paper studies the problem of fair division of indivisible goods among kk groups of n1,…,nkn_1,\ldots, n_k agents. We look at the worst downward deviation PROP(n1,…,nk)\textit{PROP}(n_1,\ldots, n_k) of an agent in a group from its 1/k1/k-share. We improve the bounds of (Manurangsi and Meka, 2026) and show that PROP(n1,…,nk)=Θ~(n/k)\textit{PROP}(n_1,\ldots, n_k) = \tildeΘ(\sqrt{n/k}), where n=n1+…+nkn = n_1 + \ldots + n_k is the total number of agents. For the proof of the upper bound, we develop novel discrepancy-type tools and, in particular, a way to efficiently work with one-sided discrepancy constraints.

Introduction

Fair division

Allocation of indivisible goods among several agents is an important theoretical and practical problem. In the classical setting, we have kk agents and a set GG of commodities to be divided among them. Each agent has their own nonnegative evaluation of the commodities. A typical goal is to provide a nearly-fair partition subject to different fairness notions. We refer the reader to [1] for survey on fair division of indivisible goods.

In recent years, group fair divisions gained a significant amount of attention. In group fair division, individual agents are replaced by groups of agents that share the same commodities. Formally, we are given kk groups of n1≥⋯≥nkn_1 \ge\cdots\ge n_k agents, and a set of objects GG. Denote by μ(i,j)(g)\mu^{(i,j)}(g) the value that the jj-th agent in the ii-th group accord to a commodity g∈Gg \in G. Importantly, all valuations are nonnegative: μ(i,j)(g)≥0\mu^{(i,j)}(g) \ge0. The goal of the partition is again to produce nearly-fair parts, now evaluated by each agent in the group. We work with the following well-studied fairness notion.

Definition 1. Let cc be a positive integer. A partition G=B1⊔B2⊔…⊔BkG = B_1 \sqcup B_2 \sqcup\ldots\sqcup B_k satisfies PROP-cc, if for any i,ji,j there is a subset C⊆G∖BiC \subseteq G \setminus B_i, ∣C∣≤c|C| \le c such that

μ(i,j)(Bi)≥μ(i,j)(G)k−μ(i,j)(C).\mu^{(i,j)}(B_i) \ge\frac{\mu^{(i,j)}(G)}{k} - \mu^{(i,j)}(C).

Given kk and n1,…,nkn_1,\ldots,n_k, let us denote by PROP⁡(n1,…,nk)\operatorname{PROP}(n_1,\ldots,n_k) the smallest cc such that for kk groups of agents of sizes n1,…,nkn_1,\ldots,n_k with any valuations there exists a partition that satisfies PROP-cc.

In a recent paper [17], the authors proved the following upper and lower bounds for this problem.

Theorem 1 (Manurangsi–Meka, [17]). We have

Ω(max⁡i∈[k]ikni)=PROP⁡(n1,…,nk)=O(n1).\Omega\left(\max_{i \in[k]} \frac{i}{k}\sqrt{n_i}\right) = \operatorname{PROP}(n_1,\ldots,n_k) = O(\sqrt{n_1}).

While being tight for groups of nearly equal size, these bounds are far apart for asymmetric groups. The main result of this paper is an almost tight upper and lower bounds for this setting in general.

Theorem 2. Fix some positive integers n1≥…≥nkn_1 \ge\ldots\ge n_k and put n=∑i=1knin = \sum_{i=1}^{k} n_i. We have

Ω(nk)=PROP⁡(n1,…,nk)=O(nk⋅(ln⁡3/2(n1⋅kn)+1)).\Omega\left(\sqrt{\frac{n}{k}}\right) = \operatorname{PROP}(n_1,\ldots,n_k) = O\left(\sqrt{\frac{n}{k}\cdot\left(\ln^{3/2}\left(n_1 \cdot\frac{k}{n}\right) + 1\right)}\right).

The authors of [17] related the problem of determining PROP to a multicolor discrepancy problem. Certain tools from discrepancy theory underlie both lower and upper bounds. In earlier papers, the authors also related fair division and discrepancy questions ([11, 19, 8, 17]). However, all previous work always worked with two-sided discrepancy problems. In private communication, Pasin Manurangsi mentioned that developing an approach to benefit from one-sided discrepancy constraints is an important challenge. In this paper, we managed to develop such an approach.

Our contributions

  • We provide nearly tight bounds for PROP;

  • We develop an approach to benefit from one-sided discrepancy constraints in discrepancy problems via the gifting mechanism;

  • We prove a discrepancy result on the so-called Top-down Beck-Fiala problem, which could be seen as its sharper, non-uniform analogue.

Next, we describe the discrepancy theory counterpart. We start with some background.

Discrepancy

Discrepancy theory studies how to color a set system so that every set remains nearly balanced. In the standard two-color matrix formulation, given a matrix A∈Rn×mA \in\mathbb{R}^{n \times m}, a coloring is a vector χ∈{−1,1}m\chi\in\{-1,1\}^{m}, and the discrepancy of AA is defined as

disc⁡(A)=min⁡χ∈{−1,1}m∥Aχ∥∞.\operatorname{disc}(A) = \min_{\chi\in\{-1,1\}^{m}} \lVert A\chi\rVert_{\infty}.

Several classical results provide general-purpose upper bounds for discrepancy. If every column of a 0-1 matrix AA has at most tt nonzero entries, then disc⁡(A)≤2t−1\operatorname{disc}(A) \le2t - 1 ([7]). For every matrix A∈[−1,1]n×mA \in[-1,1]^{n \times m}, one has disc⁡(A)=O(n)\operatorname{disc}(A) = O(\sqrt{n}) ([26]). More generally, for vectors of small Euclidean norm and any convex body of Gaussian measure at least 1/21/2, there are signs whose signed sum lies in that body; in particular, this implies an O(log⁡n)O(\sqrt{\log n}) discrepancy bound in the ℓ∞\ell_{\infty} norm for matrices with nn rows and columns of unit ℓ2\ell_{2} norm ([2]). Recent work improves Banaszczyk-type bounds for the Beck–Fiala and Komlós settings and gives efficient algorithms for finding the corresponding colorings ([4]).

For the fair-division setting considered here, the most relevant notion is that of multi-color and multi-matrix discrepancy.

Definition 2 (kk-matrix discrepancy). Let A1∈[0,1]n1×m,…,Ak∈[0,1]nk×mA^{1} \in[0,1]^{n_{1} \times m}, \ldots, A^{k} \in[0,1]^{n_{k} \times m} be matrices. Define disc⁡k(A1,…,Ak)\operatorname{disc}^{k}(A^{1},\ldots,A^{k}) to be the smallest cc, such that there exists a partition v1,…,vk∈{0,1}mv_{1},\ldots,v_{k} \in\{0,1\}^{m}, ∑vi=1\sum v_{i}=1, with the following property. For all i∈[k]i \in[k] and j∈[ni]j \in[n_{i}] we have

⟨Aji,1/k⟩−c≤⟨Aji,vi⟩≤⟨Aji,1/k⟩+c.\langle A^{i}_{j},1/k\rangle- c \le\langle A^{i}_{j},v_{i}\rangle\le\langle A^{i}_{j},1/k\rangle+ c.

Here AjiA^{i}_{j} denotes the jj-th row of matrix AiA^{i}.

kk-matrix discrepancy was called asymmetric kk-color discrepancy in [17].

We put

disc⁡k(n1,…,nk)=max⁡disc⁡(A1,…,Ak),\operatorname{disc}^{k}(n_{1},\ldots,n_{k}) = \max\operatorname{disc}(A^{1},\ldots,A^{k}),

where the maximum is taken over Ai∈[0,1]ni×mA^{i} \in[0,1]^{n_{i} \times m}, where mm could be arbitrary. We use the same notational conventions for other types of discrepancy.

Definition 3 (One-sided kk-matrix discrepancy). Let A1∈[0,1]n1×m,…,Ak∈[0,1]nk×mA^{1} \in[0,1]^{n_{1} \times m},\ldots,A^{k} \in[0,1]^{n_{k} \times m} be matrices. Define disc⁡>k(A1,…,Ak)\operatorname{disc}^{k}_{>}(A^{1},\ldots,A^{k}) similar to disc⁡k\operatorname{disc}^{k}, only with a different set of inequalities imposed:

⟨Aji,vi⟩≥⟨Aji,1/k⟩−c\langle A^{i}_{j},v_{i}\rangle\ge\langle A^{i}_{j},1/k\rangle- c

for each row AjiA^{i}_{j} of AiA^{i}.

Note that we may renormalize the valuations of each agent so that for each i,ji,j we have

max⁡g∈Gμ(i,j)(g)=1.\max_{g \in G} \mu^{(i,j)}(g) = 1.

If there is a PROP-cc partition when agents utilities are represented by matrices A1,…,AkA^{1},\ldots,A^{k}, then disc⁡>k(A1,…,Ak)≤c\operatorname{disc}^{k}_{>}(A^{1},\ldots,A^{k}) \le c. It implies the left inequality in the theorem below.

Lemma 3. *For any k≥2k \ge2 and n1≥…≥nk≥1n_{1} \ge\ldots\ge n_{k} \ge1,

disc⁡>k(n1,…,nk)≤PROP⁡(n1,…,nk)≤2⋅⌈disc⁡>k(n1,…,nk)⌉\operatorname{disc}^{k}_{>}(n_{1},\ldots,n_{k}) \le\operatorname{PROP}(n_{1},\ldots,n_{k}) \le2 \cdot\left\lceil\operatorname{disc}^{k}_{>}(n_{1},\ldots,n_{k}) \right\rceil

The right inequality was essentially proved in [17] (Theorem 5.3). They have also showed the following.

disc⁡k(n1,…,nk)=Θ(n1).\operatorname{disc}^{k}(n_1,\ldots,n_k)=\Theta(\sqrt{n_1}).

Using Lemma 3, we get that Theorem 2 is equivalent to the following result.

Theorem 4. Fix some positive integers n1≥…≥nkn_1 \ge\ldots\ge n_k and put n=∑i=1knin=\sum_{i=1}^{k}n_i. We have

Ω(nk)=disc⁡>k(n1,…,nk)=O(nk(ln⁡3/2(n1⋅kn)+1)).\Omega\left(\sqrt{\frac{n}{k}}\right)=\operatorname{disc}^{k}_{>}(n_1,\ldots,n_k)=O\left(\sqrt{\frac{n}{k}}\left(\ln^{3/2}\left(n_1\cdot\frac{k}{n}\right)+1\right)\right).

The bounds in Theorem 4 and the displayed equality before it are very different for skewed matrix sizes. Thus the one-sided nature of the problem makes it genuinely different. Little was known concerning rounding of one-sided constraints from this point of view. In this work, we develop the techniques that take advantage of the one-sided inequalities. Specifically, we introduce the gifting algorithm, described in Section 3.1. As is witnessed by our main result, the resulting approach is optimal up to polylog factors. Below, we summarize the main ideas of the proofs.

Our techniques. The proof of the upper bound has two main ingredients. The first ingredient is Theorem 13 that essentially allows to bound disc⁡k\operatorname{disc}^{k} using a binary tree-type procedure. In its core is a discrepancy theory result (Top-down Beck-Fiala) which we believe is interesting in its own right. In a nutshell, it is based on a careful accounting of the ℓ1\ell_1-norms of rows and uses the Partial Coloring Lemma. It allows to get discrepancy bounds that incorporate average row complexity after several columns were deleted.

The problem with Theorem 13 is that it gives worse errors for larger groups and for rows with larger ℓ1\ell_1-norm. As a result, its direct application would lead to suboptimal bounds. To address this issue, we introduce a rebalancing procedure called the gifting algorithm, which carefully distributes a small part of commodities between (mostly large) groups. As a result, the most problematic agents (from large groups and with large ℓ1\ell_1-norm of their valuation get the most gifts. This provides us with a necessary correction that then allows to apply Theorem 13. We note that the absence of the inequality from the other side spares us from a more elaborate tracking the weight of commodities that these ‘heavy’ agents get and gives us the necessary freedom in action.

Our lower bound is inspired by the construction used in [17]. The improvement comes from the intuition that rows of very different ℓ1\ell_1-norm should make fair division more complicated. In the construction, we have some groups consisting of all-ones vectors, and the other (larger) groups get sparse valuations that look like block-diagonal matrices with the construction of Manurangsi and Meka as a block.

Organization. In Section 1.3, we formulate our generalization of the Beck–Fiala problem. In Section 2, we prove an upper bound for this problem. In Section 3, we use this result to obtain the upper bound in Theorem 2. The gifting algorithm, another important part of the proof of the upper bound, is presented in Section 3.1. Finally, in Section 4, we prove the lower bound in Theorem 2.

Discrepancy and the Top-down Beck-Fiala Problem

In the last decade, there has been great progress on the algorithmic aspects of discrepancy. One major development was the paper [16], where the authors provided a polynomial algorithm for the classical Spencer’s ‘six standard deviations’ result [26]. They introduced the following powerful tool.

Theorem 5 (Partial Coloring Lemma, ([16]).) Let v1,…,vn∈Rmv_1,\ldots,v_n \in\mathbb{R}^m be vectors, and y∈[−1,1]my \in[-1,1]^m be a “starting” point. Let c1,…,cn≥0c_1,\ldots,c_n \ge0 be thresholds such that ∑i=1nexp⁡(−ci2/16)≤m/16\sum_{i=1}^n \exp(-c_i^2/16) \le m/16. Let δ>0\delta> 0 be a small approximation parameter. Then there exists an efficient randomized algorithm that with probability at least 0.10.1 finds a point y′∈[−1,1]my' \in[-1,1]^m such that

(i) ∣⟨y−y′,vi⟩∣≤ci∥vi∥2|\langle y-y',v_i\rangle| \le c_i\|v_i\|_2.

(ii) ∣yi′∣≥1−δ|y'_i| \ge1-\delta for at least m/2m/2 indices i∈[m]i \in[m].

Moreover, the algorithm runs in time O((m+n)3⋅δ−2⋅log⁡(nm/δ))O((m+n)^3\cdot\delta^{-2}\cdot\log(nm/\delta)).

Another classical discrepancy result is the Beck–Fiala theorem ([7]), which states that, if each column of AA has ℓ1\ell_1-norm at most Δ\Delta, then the upper bound can be improved to O(Δ)O(\Delta). They conjectured that the right upper bound is O(Δ)O(\sqrt{\Delta}). [2] proved a O(Δ⋅ln⁡n)O(\sqrt{\Delta}\cdot\ln n) bound, and the result of [16] gives the same bound and a polynomial construction. In a recent breakthrough paper, [4] showed the O(Δ)O(\sqrt{\Delta}) bound in the regime Δ≥log⁡2n\Delta\ge\log^2 n. Still, the question remains open in general.

In this paper, we prove and apply the following theorem, which deals with an asymmetric generalization of the Beck-Fiala problem.1

Definition 4 (((Δi)(\Delta_i)-structured matrix).) Given positive reals Δ1≥Δ2≥…≥Δn\Delta_1 \ge\Delta_2 \ge\ldots\ge\Delta_n, we say that a matrix A∈[0,1]n×mA \in[0,1]^{n \times m} is (Δi)(\Delta_i)-structured if for any i∈[n]i \in[n], the maximum column sum in a submatrix Ai,…,nA_{i,\ldots,n} is at most Δi\Delta_i. Formally,

∀ 1≤i≤n,max⁡j∈[m]∑k=inAk,j≤Δi.\forall\ 1 \le i \le n,\quad\max_{j \in[m]} \sum_{k=i}^{n} A_{k,j} \le\Delta_i.

We remark that the setting of the Beck–Fiala theorem and, importantly, the procedure underlying its proof, works with the setting Δ1=…=Δn\Delta_1 = \ldots= \Delta_n. However, as we illustrate in our proofs, one could draw benefits from this asymmetric setting. Let us define the linear discrepancy problem for (Δi)(\Delta_i)-structured matrices.

Theorem 6 (Bound on the top-down Beck-Fiala problem). There exists an absolute constant C>0C > 0, such that the following holds. Let n≥Δ1≥Δ2≥…≥Δn>0n \ge\Delta_1 \ge\Delta_2 \ge\ldots\ge\Delta_n > 0 be a sequence of reals. Let A∈[0,1]n×mA \in[0,1]^{n \times m} be a (Δi)(\Delta_i)-structured matrix and x0∈[0,1]mx_0 \in[0,1]^m be a point. Then there exists x∈{0,1}mx \in\{0,1\}^m, such that for all 1≤i≤n1 \le i \le n we have

∣⟨Ai,x0−x⟩∣≤Clog⁡2(4n/Δn)⋅Δilog⁡2(4Δ1/Δn).|\langle A_i,x_0-x\rangle| \le C\log_2(4n/\Delta_n)\cdot\sqrt{\Delta_i\log_2(4\Delta_1/\Delta_n)}.

We prove this theorem in Section 2. For convenience, let us introduce the following function

Ξ(s,r):=Clog⁡2(4s)⋅log⁡2(4r),(1)\Xi(s,r) := C\log_2(4s)\cdot\sqrt{\log_2(4r)}, \tag*{(1)}

where CC is a constant from Theorem 6. With this notation, the bound from Theorem 6 becomes

∣⟨Ai,x0−x⟩∣≤Ξ(n/Δn,Δ1/Δn)⋅Δi.|\langle A_i,x_0-x\rangle| \le\Xi(n/\Delta_n,\Delta_1/\Delta_n)\cdot\sqrt{\Delta_i}.

The proof is based on the idea that we can bound the number of rows having ℓ1\ell_1-norm above a certain threshold. The idea of bounding large-norm rows has been instrumental in discrepancy. It has been used in random walk techniques ([16]; [4]) and in the foundational work of ([7]). Besides this, we use the partial coloring lemma from ([16]) and draw on some of the ideas from that paper. We suspect that the polylogarithmic term log⁡2(4r)\sqrt{\log_2(4r)} in Theorem 6 is necessary, see Lemma 8.

Group fair division. Fair division among groups (families), models settings in which each group receives a common bundle while agents within the same group may have different valuations. For divisible resources, several fairness notions for families have been studied, including average, unanimous, and democratic variants of proportionality and envy-freeness [23]. Another line of work considers cake divisions among groups of prescribed sizes, with the additional requirement that every group receives a connected piece [25].

For indivisible goods, early work studied the asymptotic existence of envy-free allocations for groups under random additive valuations [18], as well as approximate maximin-share guarantees for groups [27]. Since unanimous guarantees can be too demanding when agents in the same group have heterogeneous preferences, democratic fairness notions have also been proposed for indivisible goods [24]. Further work studied EF1- and EFX-type guarantees for group allocations under several valuation classes and group models [14]. More recent directions include ordinal maximin-share guarantees, variable groups, and allocations among couples or other small groups [20]; [13]; [12].

The line of work closest to the present paper connects group fair division with discrepancy theory. There, the properties of approximate proportionality, envy-freeness, and consensus division for groups are translated into multi-color and weighted discrepancy problems [19]. Subsequent sharper lower bounds used further developed the apparatus of multi-color and weighted discrepancy techniques [8]; [17]; [12].

Discrepancy. Constructive discrepancy results often rely on random walk and semidefinite programming ideas [3]; [16]; [22]; [6]; [5]; [4]. A more detailed introduction to discrepancy theory can be found in the standard books on the discrepancy method and geometric discrepancy [9]; [21].

In the fair-division setting considered here, the relevant notions are not only ordinary two-color discrepancy, but also multi-color and weighted variants. For kk colors, one natural objective is to minimize, over all relevant rows, the largest imbalance between any two color classes. Multi-color discrepancy and its weighted variants provide a natural language for studying allocations among several groups: colors can be interpreted as groups, and a coloring can be viewed as rounding the ideal fractional allocation in which every group receives a 1/k1/k fraction of every good [10]; [11]. This viewpoint underlies discrepancy-based bounds for group fair division, including bounds for approximate proportionality, envy-freeness, and consensus division [19]; [8]; [17]; [12]. Recent work also extends discrepancy-based fair-division techniques beyond additive valuations by defining discrepancy notions for non-additive valuations and applying them to fair-division problems [15].

Notation

We denote [n]={1,2,…,n}[n] = \{1,2,\ldots,n\}.

For a matrix A∈Rn×mA \in\mathbb{R}^{n \times m} and an index i∈[n]i \in[n], by Ai∈RmA_i \in\mathbb{R}^m we denote the vector representing the ii-th row of the matrix AA. If I⊆[n]I \subseteq[n], then we denote by AI∈R∣I∣×mA_I \in\mathbb{R}^{|I| \times m} the submatrix of $A obtained by retaining only the rows ii that lie in II. If $1 \le i \le j \le n,thenwedefine, then we defineA_{i\ldots j}tobeto beA_IwherewhereI = {i,i+1,\ldots,j}.Foravector. For a vectorv \in \mathbb{R}^mandasubsetand a subsetI \subseteq [m],wedenoteby, we denote byv_I \in \mathbb{R}^{|I|}therestrictionofthe restriction ofvtothecoordinatesto the coordinatesI$.

For a vector v∈Rmv \in\mathbb{R}^m, let ∥v∥1\lVert v\rVert_1, ∥v∥2\lVert v\rVert_2, ∥v∥∞\lVert v\rVert_\infty denote its 1-norm, 2-norm, and ∞\infty-norm, respectively.

Proof of Theorem 6

We first prove the following simple, but important lemma.

Lemma 7. Let n,m∈Nn,m \in\mathbb{N} and A∈[0,1]n×mA \in[0,1]^{n \times m} be a (Δi)(\Delta_i)-structured matrix where Δ1/Δn≤r\Delta_1/\Delta_n \le r for r>0r > 0. Then, for any t>0t > 0, the number of rows ii which satisfy

∥Ai∥1>2⋅t⋅(log⁡2r+2)⋅Δi(2)\lVert A_i\rVert_1 > 2 \cdot t \cdot(\log_2 r + 2) \cdot\Delta_i \tag*{(2)}

is at most m/tm/t.

Proof. Let kk be an integer such that (log⁡2Δn)−1≤k≤log⁡2Δ1(\log_2 \Delta_n) - 1 \le k \le\log_2 \Delta_1 and assume there is at least one ii, such that 2k≤Δi≤2k+12^k \le\Delta_i \le2^{k+1}. Let dd be the smallest index, such that Δd≤2k+1\Delta_d \le2^{k+1}, and let u≥du \ge d be the largest index, such that Δu≥2k\Delta_u \ge2^k.

The sum in each column of the matrix Ad,…,uA_{d,\ldots,u} is at most 2k+12^{k+1}. So, the total sum of all entries of rows d,…,ud,\ldots,u is at most m⋅2k+1m \cdot2^{k+1}. Therefore, the number of rows i∈[d,u]i \in[d,u] that have sum of entries exceeding 2⋅t⋅(log⁡2r+2)⋅Δi2 \cdot t \cdot(\log_2 r + 2) \cdot\Delta_i is at most

m⋅2k+12⋅t⋅(log⁡2r+2)⋅2k=mt⋅(log⁡2r+2),\frac{m \cdot2^{k+1}}{2 \cdot t \cdot(\log_2 r + 2) \cdot2^k} = \frac{m}{t \cdot(\log_2 r + 2)},

where we used Δi≥2k\Delta_i \ge2^k.

Since for every Δi\Delta_i, there is an integer (log⁡2Δn)−1≤k≤log⁡2Δ1(\log_2 \Delta_n) - 1 \le k \le\log_2 \Delta_1, such that 2k≤Δi≤2k+12^k \le\Delta_i \le2^{k+1}, summing over all such kk, we obtain that the number of rows satisfying (2) is at most

(log⁡2r+2)⋅mt⋅(log⁡2r+2)=mt.(\log_2 r + 2) \cdot\frac{m}{t \cdot(\log_2 r + 2)} = \frac{m}{t}.

We remark that the log⁡2r\log_2 r term in Lemma 7 cannot be reduced as shown by the following lemma.

Lemma 8. There exists an absolute constant m0>0m_0 > 0, such that for every m≥m0m \ge m_0, there exists n≥2⋅mn \ge2 \cdot m, positive integers (Δi)(\Delta_i) where n≥Δ1≥Δ2≥…≥Δn≥1n \ge\Delta_1 \ge\Delta_2 \ge\ldots\ge\Delta_n \ge1, and a matrix A∈[0,1]n×mA \in[0,1]^{n \times m} that is (Δi)(\Delta_i)-structured and all its rows satisfy

∥Ai∥1>ln⁡m3⋅Δi.\lVert A_i\rVert_1 > \frac{\ln m}{3} \cdot\Delta_i.

Proof. We construct the matrix AA as a vertical concatenation of matrices A1,…,AkA_1,\ldots,A_k, where the maximum column sum in each matrix AiA_i is 1:

A=[AkAk−1⋮A1].A = \begin{bmatrix} A_k \\ A_{k-1} \\ \vdots\\ A_1 \end{bmatrix}.

If we look at a row rr of matrix AA that is a part of a matrix AiA_i, then the maximum column sum in rows ≥r\ge r is at most ii. Let us set Δr=i\Delta_r=i.

For each ii, we design matrix AiA_i to have (i⋅ln⁡m)/3(i\cdot\ln m)/3 ones in each row. We can make AiA_i to have at least ⌊3m/(i⋅ln⁡m)⌋\lfloor3m/(i\cdot\ln m)\rfloor rows, so that ones in each row are in disjoint set of columns. Therefore, choosing k=⌊m/ln⁡m⌋k=\lfloor m/\ln m\rfloor the total number of rows nn across all matrices is

n≥∑i=1k⌊3mi⋅ln⁡m⌋≥∑i=1k(3mi⋅ln⁡m−1)=mln⁡m⋅3∑i=1⌊m/ln⁡m⌋1i−mln⁡m≥2⋅mn\ge\sum_{i=1}^{k}\left\lfloor\frac{3m}{i\cdot\ln m}\right\rfloor\ge\sum_{i=1}^{k}\left(\frac{3m}{i\cdot\ln m}-1\right)=\frac{m}{\ln m}\cdot3\sum_{i=1}^{\lfloor m/\ln m\rfloor}\frac{1}{i}-\frac{m}{\ln m}\ge2\cdot m

for sufficiently large m≥m0m\ge m_0.

We rely on the following corollary from [16].

Corollary 9 (corollary from [16]). Let A∈[0,1]n×mA\in[0,1]^{n\times m} be a matrix and let c1,…,cn≥0c_1,\ldots,c_n\ge0 be such that ∑j=1nexp⁡(−cj2/16)≤m/16\sum_{j=1}^{n}\exp(-c_j^2/16)\le m/16. Let δ>0\delta>0 be a small approximation parameter. Then there exists an efficient randomized algorithm which for any point x∈[0,1]mx\in[0,1]^m with probability at least 0.10.1 finds a point x′∈[0,1]mx'\in[0,1]^m such that

(i) ∣⟨x−x′,Ai⟩∣≤ci∥Ai∥2|\langle x-x',A_i\rangle|\le c_i\lVert A_i\rVert_2.

(ii) xi′≥1−δx'_i\ge1-\delta or xi′≤δx'_i\le\delta for at least m/2m/2 indices i∈[m]i\in[m].

Moreover, the algorithm runs in time poly⁡(n,m,δ−1)\operatorname{poly}(n,m,\delta^{-1}).

We also need to use the following standard lemma.

Lemma 10 ([7]). Let m>n>0m>n>0 be integers, let A∈[0,1]n×mA\in[0,1]^{n\times m} be a matrix, and x∈[0,1]mx\in[0,1]^m be a point. Then, there exists a point x′∈[0,1]mx'\in[0,1]^m, such that Ax=Ax′Ax=Ax' and all but nn coordinates of x′x' are rounded (either 11 or 00).

Now, let us recall the statement of Theorem 6 and pose its proof. Our proof is based on the ideas that [16] use to prove their Theorem 2.

Theorem 11 (Restatement of Theorem 6). There exists an absolute constant C>0C>0, such that the following holds. Let n≥Δ1≥Δ2≥…≥Δn>0n\ge\Delta_1\ge\Delta_2\ge\ldots\ge\Delta_n>0 be a sequence of reals. Let A∈[0,1]n×mA\in[0,1]^{n\times m} be a (Δi)(\Delta_i)-structured matrix and x0∈[0,1]mx_0\in[0,1]^m be a point, then there exists x∈{0,1}mx\in\{0,1\}^m, such that for all 1≤i≤n1\le i\le n,

∣⟨Ai,x0−x⟩∣≤Clog⁡2(4n/Δn)⋅Δilog⁡2(4Δ1/Δn).|\langle A_i,x_0-x\rangle|\le C\log_2(4n/\Delta_n)\cdot\sqrt{\Delta_i}\log_2(4\Delta_1/\Delta_n).

Proof of Theorem 6. Let r=Δ1/Δnr=\Delta_1/\Delta_n. By Lemma 10 we can round m−nm-n coordinates of x0x_0 if m>nm>n. Therefore, let us assume that m≤nm\le n.

Let us choose ci=C⋅(log⁡2(r)+2)⋅Δi/∥Ai∥2c_i=C\cdot\sqrt{(\log_2(r)+2)\cdot\Delta_i/\lVert A_i\rVert_2} and show that ∑i=1nexp⁡(−ci2/16)≤m/16\sum_{i=1}^{n}\exp(-c_i^2/16)\le m/16 for CC to be chosen later. For each integer ss, let qsq_s denote the number of rows for which it holds

2s⋅(log⁡2(r)+2)⋅Δi≤∥Ai∥22≤2s+1⋅(log⁡2(r)+2)⋅Δi.2^s\cdot(\log_2(r)+2)\cdot\Delta_i\le\lVert A_i\rVert_2^2\le2^{s+1}\cdot(\log_2(r)+2)\cdot\Delta_i.

Note that ∥Ai∥22≤∥Ai∥1\lVert A_i\rVert_2^2\le\lVert A_i\rVert_1 and, thus, by Lemma 7, we have qs≤m/2s−1q_s\le m/2^{s-1}. Thus,

∑i=1nexp⁡(−ci2/16)≤∑s=−∞∞qs⋅exp⁡(−C216⋅2s+1)≤≤m⋅∑s=−∞∞exp⁡(−C216⋅2s+1)/2s−1==m⋅2⋅∑s=−∞∞2s⋅exp⁡(−2s⋅C232)≤≤m⋅2⋅(∑s=0∞2s⋅exp⁡(−2s⋅C232)+∑s=0∞2−s⋅exp⁡(−C232⋅2s))≤m/16\begin{aligned} \sum_{i=1}^{n}\exp(-c_i^2/16)\le\sum_{s=-\infty}^{\infty}q_s\cdot\exp\left(-\frac{C^2}{16\cdot2^{s+1}}\right)\le\\ &\le m \cdot\sum_{s=-\infty}^{\infty} \exp\left(-\frac{C^2}{16 \cdot2^{s+1}}\right) / 2^{s-1} = \\ &= m \cdot2 \cdot\sum_{s=-\infty}^{\infty} 2^s \cdot\exp\left(-2^s \cdot\frac{C^2}{32}\right) \le\\ &\le m \cdot2 \cdot\left(\sum_{s=0}^{\infty} 2^s \cdot\exp\left(-2^s \cdot\frac{C^2}{32}\right) + \sum_{s=0}^{\infty} 2^{-s} \cdot\exp\left(-\frac{C^2}{32 \cdot2^s}\right)\right) \le m/16 \end{aligned}

for sufficiently large C>0C > 0.

Now, apply Corollary 9 for c1,…,cnc_1,\ldots,c_n and δ=Δn/(3n2)\delta= \sqrt{\Delta_n}/(3n^2), obtaining a vector x1x_1 that has at least m/2m/2 coordinates nearly rounded and ∣⟨x1−x0,Ai⟩∣≤ci∥Ai∥2=C⋅(log⁡2(r)+2)⋅Δi\left|\langle x_1-x_0,A_i\rangle\right| \le c_i\|A_i\|_2 = C \cdot\sqrt{(\log_2(r)+2)\cdot\Delta_i}. Let I1⊆[m]I_1 \subseteq[m] denote the coordinates of x1x_1 that are not δ\delta-close to 0 or 1. Let us recursively apply this procedure to (x1)∣I1(x_1)|_{I_1}, obtaining a new point x2∈[0,1]mx_2 \in[0,1]^m with a set I2⊆I1I_2 \subseteq I_1 of at most m/4m/4 coordinates not δ\delta-close to 0 or 1. By iterating this at most log⁡2(m/Δn)+2\log_2(m/\Delta_n)+2 times, we get a sequence of points x1,x2,…,xlx_1,x_2,\ldots,x_l, for which l≤log⁡2(m/Δn)+2l \le\log_2(m/\Delta_n)+2 and xlx_l has all but min⁡(Δn,n/4)\min(\Delta_n,n/4) coordinates nearly rounded, denoted by IlI_l. By triangle inequality, for each ii we have

∣⟨xl−x0,Ai⟩∣≤C⋅(log⁡2(m/Δn)+2)⋅(log⁡2(r)+2)⋅Δi.|\langle x_l-x_0,A_i\rangle| \le C\cdot(\log_2(m/\Delta_n)+2)\cdot\sqrt{(\log_2(r)+2)\cdot\Delta_i}.

Now, let us nearly round the remaining m′≤min⁡(Δn,n/4)m' \le\min(\Delta_n,n/4) coordinates. For this, we choose ci=8⋅ln⁡(n/m′)c_i=8\cdot\sqrt{\ln(n/m')}. Then, ∑i=1nexp⁡(−ci2/16)=∑i=1nexp⁡(−4⋅ln⁡(n/m′))=n⋅(m′/n)4≤m′/16\sum_{i=1}^n\exp(-c_i^2/16)=\sum_{i=1}^n\exp(-4\cdot\ln(n/m'))=n\cdot(m'/n)^4\le m'/16. Applying Corollary 9 to the sequence cic_i, δ=Δn/(3n2)\delta=\sqrt{\Delta_n}/(3n^2), and point (xl)∣Il(x_l)|_{I_l}, we get a point xl+1∈[0,1]mx_{l+1}\in[0,1]^m that has at most m′/2m'/2 coordinates Il+1⊆IlI_{l+1}\subseteq I_l not δ\delta-close to 0 or 1 and ∣⟨xl−xl+1,Ai⟩∣≤ci⋅m′=8ln⁡(n/m′)⋅m′|\langle x_l-x_{l+1},A_i\rangle|\le c_i\cdot\sqrt{m'}=8\sqrt{\ln(n/m')}\cdot\sqrt{m'} for every ii. We iterate this procedure until all coordinates get nearly rounded. Denote xl+l′x_{l+l'} to be the resulting point. Then, by triangle inequality, for every ii,

∣⟨xl−xl+l′,Ai⟩∣≤8⋅∑i=0∞ln⁡(2i⋅n/m′)⋅m′/2i≤C′⋅ln⁡(n/m′)⋅m′|\langle x_l-x_{l+l'},A_i\rangle|\le8\cdot\sum_{i=0}^{\infty}\sqrt{\ln(2^i\cdot n/m')}\cdot\sqrt{m'/2^i}\le C'\cdot\sqrt{\ln(n/m')\cdot m'}

for some constant C′>0C'>0.

Let xx be xl+l′x_{l+l'} in which we round each coordinate to the nearest integer. Then, the error of the ii-th row is estimated by

∣⟨x−x0,Ai⟩∣≤C⋅(log⁡2(m/Δn)+2)⋅(log⁡2(r)+2)⋅Δi+C′⋅ln⁡(n/m′)⋅m′+Δn3n≤C⋅(log⁡2(m/Δn)+2)⋅(log⁡2(r)+2)⋅Δi+C′⋅ln⁡(n/Δn)⋅Δn+Δn3n≤(C+C′+1)⋅(log⁡2(n/Δn)+2)⋅(log⁡2(r)+2)⋅Δi=\begin{aligned} |\langle x-x_0,A_i\rangle|\le C\cdot(\log_2(m/\Delta_n)+2)\cdot\sqrt{(\log_2(r)+2)\cdot\Delta_i}+C'\cdot\sqrt{\ln(n/m')\cdot m'}+\frac{\sqrt{\Delta_n}}{3n} \\ &\le C\cdot(\log_2(m/\Delta_n)+2)\cdot\sqrt{(\log_2(r)+2)\cdot\Delta_i}+C'\cdot\sqrt{\ln(n/\Delta_n)\cdot\Delta_n}+\frac{\sqrt{\Delta_n}}{3n} \\ &\le(C+C'+1)\cdot(\log_2(n/\Delta_n)+2)\cdot\sqrt{(\log_2(r)+2)\cdot\Delta_i} = \end{aligned}
(C+C′+1)⋅log⁡2(4n/Δn)⋅Δilog⁡2(4r).(C+C'+1)\cdot\log_2(4n/\Delta_n)\cdot\sqrt{\Delta_i\log_2(4r)}.

In fact, in the proof of the upper bound for PROP, we shall need a slight generalization of Theorem 6 for a vertical concatenation of (Δi)(\Delta_i)-structured matrices.

Theorem 12. Let n1,…,nkn_1,\ldots,n_k and mm be positive integers. Let A1∈[0,1]n1×m,…,Ak∈[0,1]nk×mA^1 \in[0,1]^{n_1 \times m},\ldots,A^k \in[0,1]^{n_k \times m} be (Δj1),(Δj2),…,(Δjk)(\Delta^1_j),(\Delta^2_j),\ldots,(\Delta^k_j)-structured matrices respectively. Let rr be the maximum ratio r=max⁡i,j,i′,j′ΔjiΔj′i′r = \max_{i,j,i',j'} \frac{\Delta^i_j}{\Delta^{i'}_{j'}}, and let dd be the minimum d=min⁡i,jΔjid = \min_{i,j} \Delta^i_j. Then, for any x0∈[0,1]mx_0 \in[0,1]^m, there exists x∈{0,1}mx \in\{0,1\}^m, such that for all i=1…ki=1\ldots k and j=1…nij=1\ldots n_i

∣⟨Aji,x0−x⟩∣≤Ξ((n1+…+nk)/d,r)⋅k⋅Δji.\left|\langle A^i_j,x_0-x\rangle\right| \le\Xi((n_1+\ldots+n_k)/d,r)\cdot\sqrt{k\cdot\Delta^i_j}.

Proof. Let n=n1+…+nkn=n_1+\ldots+n_k denote the total number of rows in all matrices. Let us sort all values of Δji\Delta^i_j in decreasing order: Δj1i1≥Δj2i2≥…≥Δjnin\Delta^{i_1}_{j_1} \ge\Delta^{i_2}_{j_2} \ge\ldots\ge\Delta^{i_n}_{j_n}. Consider a matrix AA obtained from matrices A1,…,AkA^1,\ldots,A^k by stacking vertically their rows in the order (i1,j1),(i2,j2),…,(in,jn)(i_1,j_1),(i_2,j_2),\ldots,(i_n,j_n)

A=[Aj1i1Aj2i2⋮Ajnin]A = \begin{bmatrix} A^{i_1}_{j_1}\\ A^{i_2}_{j_2}\\ \vdots\\ A^{i_n}_{j_n} \end{bmatrix}

We claim that this is a (k⋅Δjlil)l=1n(k\cdot\Delta^{i_l}_{j_l})_{l=1}^n-structured matrix. Indeed, consider ll-th row of this matrix. Since we sorted Δji\Delta^i_j in decreasing order, only rows that have Δ\Delta less or equal to Δjlil\Delta^{i_l}_{j_l} are below ll-th row. Let us fix one of the initial matrices r∈{1,2,…,k}r \in\{1,2,\ldots,k\} and look only at the subset of its rows I={q∣Δqr≤Δjlil}I=\{q\mid\Delta^r_q\le\Delta^{i_l}_{j_l}\}. By definition, we know that the maximum sum in a column in AIrA^r_I is at most Δjlil\Delta^{i_l}_{j_l}. Therefore, the maximum sum in a column of the submatrix

[AjlilAjl+1il+1⋮Ajnin]\begin{bmatrix} A^{i_l}_{j_l}\\ A^{i_{l+1}}_{j_{l+1}}\\ \vdots\\ A^{i_n}_{j_n} \end{bmatrix}

can be bounded by k⋅Δjlilk\cdot\Delta^{i_l}_{j_l} by summing over all kk matrices. Now, applying Theorem 6 to matrix AA, we obtain the desired bound of Ξ(n/d,r)⋅k⋅Δjlil\Xi(n/d,r)\cdot\sqrt{k\cdot\Delta^{i_l}_{j_l}} on linear discrepancy.

Upper Bound in Theorem 4

Recall, that the goal is to prove

disc⁡>k(n1,…,nk)=O(nk⋅(ln⁡3/2n1⋅(n1⋅kn+1))).\operatorname{disc}^k_{>}(n_1,\ldots,n_k)=O\left(\sqrt{\frac{n}{k}\cdot\left(\ln^{3/2}n_1\cdot\left(n_1\cdot\frac{k}{n}+1\right)\right)}\right).

Throughout this section, we use group fair division interpretation of disc⁡(A1,…,Ak)\operatorname{disc}(A^1,\ldots,A^k). Let n1≥…≥nk≥1n_1\ge\ldots\ge n_k\ge1 be the group sizes, where n=n1+…+nkn=n_1+\ldots+n_k is the total number of agents. The matrices A1∈[0,1]n1×m,…,Ak∈[0,1]nk×mA^1\in[0,1]^{n_1\times m},\ldots,A^k\in[0,1]^{n_k\times m} denote agents utilities, and the goal is to find a partition v1,…,vk∈{0,1}mv_1,\ldots,v_k\in\{0,1\}^m, ∑vi=1\sum v_i=1, such that for any i,ji,j we have

⟨Aji,vi⟩≥⟨Aji,1/k⟩−c,\langle A^i_j,v_i\rangle\ge\langle A^i_j,1/k\rangle-c,

with cc as small as possible.

The following theorem describes a partition algorithm that control the discrepancy of rows in terms of the (Δj)(\Delta_j)-structure of the group matrices.

Theorem 13. There exists an absolute constant C′>0C' > 0, such that the following holds. Suppose agents valuations are presented by matrices A1∈[0,1]n1×m,…,Ak∈[0,1]nk×mA^1 \in[0,1]^{n_1 \times m}, \ldots, A^k \in[0,1]^{n_k \times m} that are (Δj1),(Δj2),…,(Δjk)(\Delta^1_j),(\Delta^2_j),\ldots,(\Delta^k_j)-structured, respectively. Let rr upper bound the ratio between Δji\Delta^i_j:

max⁡i,j,i′,j′ΔjiΔj′i′≤r.\max_{i,j,i',j'} \frac{\Delta^i_j}{\Delta^{i'}_{j'}} \le r.

Let d≤n1d \le n_1 be the lower bound on the minimum value of Δ\Delta: d≤min⁡i,jΔjid \le\min_{i,j} \Delta^i_j. Then, there is an allocation v1,…,vk∈{0,1}mv_1,\ldots,v_k \in\{0,1\}^m, ∑vi=1\sum v_i = 1, such that for any i∈[k]i \in[k] and j∈[ni]j \in[n_i] we have

∣⟨Aji,vi⟩−⟨Aji,1/k⟩∣≤C′⋅Ξ(n1/d,r)⋅Δji\left|\langle A^i_j,v_i\rangle-\langle A^i_j,1/k\rangle\right| \le C'\cdot\Xi(n_1/d,r)\cdot\sqrt{\Delta^i_j}

Proof. The role of this theorem is to provide a high-level partition of our goods.

Let us recursively divide objects between groups. If k=1k = 1, we give all the objects to the only group, incurring zero error. For k≥2k \ge2, we divide all groups into two parts. The first part consists of the first ⌊k/2⌋\lfloor k/2\rfloor groups, and the second part consists of the remaining ⌈k/2⌉\lceil k/2\rceil groups. We apply Theorem 12 with x0=1⋅⌊k/2⌋kx_0 = 1 \cdot\frac{\lfloor k/2\rfloor}{k} to obtain a rounded vector x∈{0,1}mx \in\{0,1\}^m. We recursively divide objects that correspond to ones in xx between the first ⌊k/2⌋\lfloor k/2\rfloor groups and divide objects that correspond to zeros in xx between the remaining ⌈k/2⌉\lceil k/2\rceil groups. Let viv_i be the indicator vector of the set of objects that the ii-th groups receives at the end of this division.

This division procedure can be viewed as a binary tree. In every node of the tree, we have a segment of groups between which we divide the set of objects assigned to this node. If this is a non-leaf node, then some objects go to the left child and then rest of the objects go to the right child.

For the jj-th agent from the ii-th group, let us estimate his evaluation of the goods received. Take the path PP from the root of the binary tree to the leaf which corresponds to the ii-th group. Define ktk_t to be the number of groups that participated in the division in tt-th node of path PP. We have k1=kk_1 = k and kl=1k_l = 1 where ll is the number of nodes in PP. Moreover, we also have kt+1k_{t+1} to be equal either to ⌊kt/2⌋\lfloor k_t/2\rfloor or ⌈kt/2⌉\lceil k_t/2\rceil for every 1≤t≤l−11 \le t \le l-1.

For t=1…lt = 1\ldots l, let xt∈{0,1}mx_t \in\{0,1\}^m denote the indicator of the set of goods that was divided in the tt-th node of the path PP. We have x1=1x_1 = 1 and xl=vix_l = v_i. Let us check that

∣⟨Aji,xt+1⟩−⟨Aji,xtkt+1kt⟩∣≤Ξ(n1⋅kt/d,r)⋅kt⋅Δji.(3)\left|\left\langle A^i_j,x_{t+1}\right\rangle-\left\langle A^i_j,x_t\frac{k_{t+1}}{k_t}\right\rangle\right| \le\Xi(n_1\cdot k_t/d,r)\cdot\sqrt{k_t\cdot\Delta^i_j}. \tag*{(3)}

To divide objects in the tt-th node, we apply Theorem 12 to the submatrix of AA narrowed down to rows and columns that correspond to agents and objects that participate in this division. The theorem gives us a vector x′′x''. Let x′∈{0,1}mx' \in\{0,1\}^m be the vector x′′x'' padded with zeros for objects that did not participate in this division. By guarantees of Theorem 12, x′x' satisfies

∣⟨Aji,x′⟩−⟨Aji,xt⌊kt/2⌋kt⟩∣≤Ξ(n1⋅kt/d,r)⋅kt⋅Δji.\left|\left\langle A^i_j,x'\right\rangle-\left\langle A^i_j,x_t\frac{\lfloor k_t/2\rfloor}{k_t}\right\rangle\right| \le\Xi(n_1\cdot k_t/d,r)\cdot\sqrt{k_t\cdot\Delta^i_j}.

Here we used that the total number of agents in groups that participated in the division in the tt-th node of the path PP can be bounded by n1⋅ktn_1\cdot k_t, since n1n_1 is the largest number of agents in any group.

If the ii-th group is in the first ⌊kt/2⌋\lfloor k_t/2\rfloor groups in this division, then we have xt+1=x′x_{t+1}=x', concluding (3) in this case. If it is among the remaining ⌈kt/2⌉\lceil k_t/2\rceil groups, then we have xt+1=xt−x′x_{t+1}=x_t-x'. Since

1=⌊kt/2⌋+⌈kt/2⌉kt,1=\frac{\lfloor k_t/2\rfloor+\lceil k_t/2\rceil}{k_t},

we have

∣⟨Aji,xt+1⟩−⟨Aji,xtkt+1kt⟩∣=∣⟨Aji,xt−x′⟩−⟨Aji,xt⌊kt/2⌋kt⟩∣==∣⟨Aji,xt⌊kt/2⌋kt⟩−⟨Aji,x′⟩∣≤Ξ(n1⋅kt/d,r)⋅kt⋅Δji.\begin{aligned} \left|\left\langle A^i_j,x_{t+1}\right\rangle-\left\langle A^i_j,x_t\frac{k_{t+1}}{k_t}\right\rangle\right|=\left|\left\langle A^i_j,x_t-x'\right\rangle-\left\langle A^i_j,x_t\frac{\lfloor k_t/2\rfloor}{k_t}\right\rangle\right|= \\ &=\left|\left\langle A^i_j,x_t\frac{\lfloor k_t/2\rfloor}{k_t}\right\rangle-\left\langle A^i_j,x'\right\rangle\right|\leq\Xi(n_1\cdot k_t/d,r)\cdot\sqrt{k_t\cdot\Delta^i_j}. \end{aligned}

which concludes (3).

By the triangle inequality, we have:

∣⟨Aji,xl⟩−⟨Aji,1k1⟩∣≤∣⟨Aji,xl⟩−⟨Aji,xl−1kl−1⟩∣+⋯+∣⟨Aji,x2k2⟩−⟨Aji,x1k1⟩∣≤Δji⋅(Ξ(n1⋅kl−1d,r)⋅kl−1kl+⋯+Ξ(n1⋅k1d,r)⋅k1k2)≤Δji⋅Ξ(n1d,r)(kl−1kl+⋯+k1k2)+Δji⋅C⋅log⁡2(4r)(log⁡2(kl−1)⋅kl−1kl+⋯+log⁡2(k1)⋅k1k2).\begin{aligned} \left|\left\langle A^i_j,x_l\right\rangle-\left\langle A^i_j,\frac{1}{k}\mathbf{1}\right\rangle\right|\leq\left|\left\langle A^i_j,x_l\right\rangle-\left\langle A^i_j,\frac{x_{l-1}}{k_{l-1}}\right\rangle\right|+\cdots+\left|\left\langle A^i_j,\frac{x_2}{k_2}\right\rangle-\left\langle A^i_j,\frac{x_1}{k_1}\right\rangle\right| \\ &\leq\sqrt{\Delta^i_j}\cdot\left(\Xi\left(\frac{n_1\cdot k_{l-1}}{d},r\right)\cdot\frac{\sqrt{k_{l-1}}}{k_l}+\cdots+\Xi\left(\frac{n_1\cdot k_1}{d},r\right)\cdot\frac{\sqrt{k_1}}{k_2}\right) \\ &\leq\sqrt{\Delta^i_j}\cdot\Xi\left(\frac{n_1}{d},r\right)\left(\frac{\sqrt{k_{l-1}}}{k_l}+\cdots+\frac{\sqrt{k_1}}{k_2}\right) \\ &+\sqrt{\Delta^i_j}\cdot C\cdot\sqrt{\log_2(4r)}\left(\frac{\log_2(k_{l-1})\cdot\sqrt{k_{l-1}}}{k_l}+\cdots+\frac{\log_2(k_1)\cdot\sqrt{k_1}}{k_2}\right). \end{aligned}

In the last inequality, we used (1) and the identity Ξ(a⋅b,r)=Ξ(a,r)+C⋅log⁡2(4r)⋅log⁡2b\Xi(a\cdot b,r)=\Xi(a,r)+C\cdot\sqrt{\log_2(4r)}\cdot\log_2 b. We have kt+1≤⌊kt/2⌋≤23ktk_{t+1}\leq\lfloor k_t/2\rfloor\leq\frac{2}{3}k_t for kt≥3k_t\geq3. Thus the expressions

(kl−1kl+⋯+k1k2)\left(\frac{\sqrt{k_{l-1}}}{k_l}+\cdots+\frac{\sqrt{k_1}}{k_2}\right)

and

(log⁡2(kl−1)⋅kl−1kl+⋯+log⁡2(k1)⋅k1k2)\left(\frac{\log_2(k_{l-1})\cdot\sqrt{k_{l-1}}}{k_l}+\cdots+\frac{\log_2(k_1)\cdot\sqrt{k_1}}{k_2}\right)

can be upper bounded by a constant. Since n1≥dn_1\geq d, we have Ξ(n1d,r)≥C⋅log⁡2(4r)\Xi\left(\frac{n_1}{d},r\right)\geq C\cdot\sqrt{\log_2(4r)}. Therefore, the final expression can be bounded by

C′⋅Ξ(n1d,r)⋅ΔjiC'\cdot\Xi\left(\frac{n_1}{d},r\right)\cdot\sqrt{\Delta^i_j}

for some constant C′>0C'>0.

Let A1∈[0,1]n1×m,…,Ak∈[0,1]nk×mA^1\in[0,1]^{n_1\times m},\ldots,A^k\in[0,1]^{n_k\times m} be agents valuations of objects. Note that AiA^i is an (ni,ni,…,ni)(n_i,n_i,\ldots,n_i)-structured matrix since it has nin_i rows. If we directly apply Theorem 13 with Δji=ni\Delta^i_j=n_i, then we obtain a division of goods where the agents from the ii-th group are at most

C′⋅Ξ(n1/nk,n1/nk)⋅ni=O~(ni)C'\cdot\Xi(n_1/n_k,n_1/n_k)\cdot\sqrt{n_i}=\widetilde{O}(\sqrt{n_i})

unsatisfied. This is an acceptable error if the sizes of all groups are equal n1=n2=…=nkn_1=n_2=\ldots=n_k. However, in an unbalanced case where, for example, n1n_1 is large and n2=n3=…=nk=1n_2=n_3=\ldots=n_k=1, the algorithm in Theorem 13 is disadvantageous for the large group and, as a result, gives an error of Ω(n1)\Omega(\sqrt{n_1}), which is too large.

To deal with this issue, we show that it is possible to rebalance the advantage/disadvantage that the groups get when applying Theorem 13 by gifting, or simply assigning, a carefully chosen amount of goods to the larger groups.

Theorem 14. Let C′′>0C^{\prime\prime} > 0. Let A1∈[0,1]n1×m,…,Ak∈[0,1]nk×mA^{1} \in[0,1]^{n_{1}\times m},\ldots,A^{k} \in[0,1]^{n_{k}\times m} be agents valuations of objects. Then, there exists a partial assignment of goods to groups v1,…,vk∈{0,1}mv_{1},\ldots,v_{k} \in\{0,1\}^{m}, ∑vi≤1\sum v_{i} \le1, such that the following properties are satisfied. Let l=⟨1,v1+⋯+vk⟩l = \langle1,v_{1}+\cdots+v_{k}\rangle denote the number of assigned goods. Let the matrix A^i∈[0,1]ni×(m−l)\hat{A}^{i} \in[0,1]^{n_{i}\times(m-l)} be obtained from AiA^{i} by removing columns that correspond to the assigned items (the ones in v1,…,vkv_{1},\ldots,v_{k}). Let Υ:=Ξ(n1kn,n1kn)\Upsilon:= \Xi\left(\frac{n_{1}k}{n},\frac{n_{1}k}{n}\right) where n=n1+⋯+nkn=n_{1}+\cdots+n_{k}. Then,

(1) l≤C′′⋅Υ⋅n⋅kl \le C^{\prime\prime}\cdot\sqrt{\Upsilon}\cdot n\cdot k.

(2) For every group ii, there is a permutation of agents p1i,…,pniip^{i}_{1},\ldots,p^{i}_{n_{i}} in this group and values n1≥Δ1i≥⋯≥Δnii≥n/(Υ⋅k)n_{1}\ge\Delta^{i}_{1}\ge\cdots\ge\Delta^{i}_{n_{i}}\ge n/(\Upsilon\cdot k), such that the matrix BiB^{i} obtained from rearranging rows of A^i\hat{A}^{i} according to the permutation pip^{i}

Bi=[A^p1iiA^p2ii⋮A^pniii]B^{i}= \begin{bmatrix} \hat{A}^{i}_{p^{i}_{1}}\\ \hat{A}^{i}_{p^{i}_{2}}\\ \vdots\\ \hat{A}^{i}_{p^{i}_{n_{i}}} \end{bmatrix}

is (Δp1ii,Δp2ii,…,Δpniii)(\Delta^{i}_{p^{i}_{1}},\Delta^{i}_{p^{i}_{2}},\ldots,\Delta^{i}_{p^{i}_{n_{i}}})-structured and for every j∈[ni]j\in[n_{i}],

max⁡(⟨Bji,vi⟩,C′′⋅Υ⋅nk)≥C′′⋅Υ⋅Δji.\max\left(\langle B^{i}_{j},v_{i}\rangle,C^{\prime\prime}\cdot\sqrt{\Upsilon\cdot\frac{n}{k}}\right)\ge C^{\prime\prime}\cdot\Upsilon\cdot\sqrt{\Delta^{i}_{j}}.

Let us comment on the statement of the theorem. The bound on ll is chosen so that the gifted items have negative affect that is covered by the error term. The condition of being (Δji)(\Delta^{i}_{j})-structured is a technical condition needed for an application of the discrepancy results and essentially allows to bound the row norms of a matrix and its submatrices. The maximum in the statement encodes the two scenarios. One option is that a given row has small ℓ1\ell_{1}-norm. In this case, we do not need to rebalance it for the application of Theorem 13 and may ignore the contribution of gifting. The other option is that a row has large ℓ1\ell_{1}-norm, and Theorem 13 gives a large error on this row. Then the gifting completely covers this possible error.

We first derive the main theorem from these ingredients. To this end, we need the following auxiliary technical lemma.

Lemma 15. There exists a constant C′′>0C^{\prime\prime}>0, such that for all r≥1r\ge1,

Ξ(Ξ(r,r)⋅r,Ξ(r,r)⋅r)≤C′′⋅Ξ(r,r).\Xi(\Xi(r,r)\cdot r,\Xi(r,r)\cdot r)\le C^{\prime\prime}\cdot\Xi(r,r).

Proof.

Ξ(Ξ(r,r)⋅r,Ξ(r,r)⋅r)=C⋅(log⁡2(r)+log⁡2(Ξ(r,r))+2)log⁡2(r)+log⁡2(Ξ(r,r))+2=O((log⁡2(r)+2)log⁡2(r)+2),since 1≤Ξ(r,r)=O(r)≤C′′⋅Ξ(r,r)\begin{aligned} \Xi(\Xi(r,r)\cdot r,\Xi(r,r)\cdot r) &= C\cdot(\log_{2}(r)+\log_{2}(\Xi(r,r))+2)\sqrt{\log_{2}(r)+\log_{2}(\Xi(r,r))+2}\\ &= O\left((\log_{2}(r)+2)\sqrt{\log_{2}(r)+2}\right),\quad\text{since }1\le\Xi(r,r)=O(r)\\ &\le C^{\prime\prime}\cdot\Xi(r,r) \end{aligned}

for some C′′>0C^{\prime\prime}>0.

Now we are ready to combine the procedure in Theorem 13 and the gifting of Theorem 14 to establish the upper bound in Theorem 2.

Proof of the upper bound in Theorem 2. Let C′C' be the constant from Theorem 13. Denote n=n1+…+nkn = n_1 + \ldots+ n_k and let Υ:=Ξ(n1kn,n1kn)\Upsilon:= \Xi\left(\frac{n_1 k}{n},\frac{n_1 k}{n}\right). By Lemma 15, there is a constant C′′>0C'' > 0, such that

C′⋅Ξ(1⋅n1kn,1⋅n1kn)≤C′′⋅Ξ(n1kn,n1kn).C' \cdot\Xi\left(1 \cdot\frac{n_1 k}{n}, 1 \cdot\frac{n_1 k}{n}\right) \le C'' \cdot\Xi\left(\frac{n_1 k}{n},\frac{n_1 k}{n}\right).

Apply Theorem 14 with constant C′′C'' to obtain a partial assignment v1,…,vk∈{0,1}mv_1,\ldots,v_k \in\{0,1\}^m and matrices B1,…,BkB^1,\ldots,B^k that are (Δpj11),(Δpj22),…,(Δpjkk)(\Delta^1_{p^1_j}),(\Delta^2_{p^2_j}),\ldots,(\Delta^k_{p^k_j}) structured, respectively. Let l=⟨1,v1+…+vk⟩l = \langle1,v_1+\ldots+v_k\rangle denote the number of goods that were assigned by the application of Theorem 14.

Now, we apply Theorem 13 to matrices B1,…,BkB^1,\ldots,B^k with r=Υ⋅n1knr = \Upsilon\cdot\frac{n_1 k}{n} and d=n/(Υ⋅k)d = n/(\Upsilon\cdot k) to divide the remaining m−lm-l objects among kk groups. Note that dr=n1dr=n_1, and the matrices BiB^i obviously satisfy Δpjii≤n1\Delta^i_{p^i_j}\le n_1. Thus, such a choice of r,dr,d is valid. Denote by u1,…,uk∈{0,1}m−lu_1,\ldots,u_k \in\{0,1\}^{m-l} the division given by the theorem. We are guaranteed that for the jj-th agent from the ii-th group, it holds

⟨Bpjii,ui⟩≥⟨Bpjii,1/k⟩−C′⋅Ξ(Υn1kn,Υn1kn)⋅Δpjii≥⟨Bpjii,1/k⟩−C′′⋅Υ⋅Δpjii.\left\langle B^i_{p^i_j},u_i\right\rangle\ge\left\langle B^i_{p^i_j},1/k\right\rangle-C'\cdot\Xi\left(\Upsilon\frac{n_1k}{n},\Upsilon\frac{n_1k}{n}\right)\cdot\sqrt{\Delta^i_{p^i_j}}\ge\left\langle B^i_{p^i_j},1/k\right\rangle-C''\cdot\Upsilon\cdot\sqrt{\Delta^i_{p^i_j}}.

Observe that max⁡{⟨Aji,vi⟩,C′′⋅Υ⋅nk}≥C′′⋅Υ⋅Δji\max\left\{\langle A^i_j,v_i\rangle,C''\cdot\sqrt{\Upsilon\cdot\frac{n}{k}}\right\}\ge C''\cdot\Upsilon\cdot\sqrt{\Delta^i_j} implies ⟨Aji,vi⟩≥C′′⋅Υ⋅Δji−C′′⋅Υ⋅nk\langle A^i_j,v_i\rangle\ge C''\cdot\Upsilon\cdot\sqrt{\Delta^i_j}-C''\cdot\sqrt{\Upsilon\cdot\frac{n}{k}}.

In addition, since at most ll objects are gifted, we have ⟨Bpjii,1/k⟩≥⟨Aji,1/k⟩−l/k≥⟨Aji,1/k⟩−C′′⋅Υ⋅nk\langle B^i_{p^i_j},1/k\rangle\ge\langle A^i_j,1/k\rangle-l/k\ge\langle A^i_j,1/k\rangle-C''\cdot\sqrt{\Upsilon\cdot\frac{n}{k}}.

Therefore, combined with the gifted objects viv_i, the jj-th agent from the ii-th group values the assigned objects at least by

⟨Aji,vi⟩+⟨Bpjii,ui⟩≥⟨Aji,1/k⟩−2⋅C′′⋅Υ⋅nk.\langle A^i_j,v_i\rangle+\langle B^i_{p^i_j},u_i\rangle\ge\langle A^i_j,1/k\rangle-2\cdot C''\cdot\sqrt{\Upsilon\cdot\frac{n}{k}}.

Thus, we have obtained the desired 2⋅C′′⋅Υ⋅nk2\cdot C''\cdot\sqrt{\Upsilon\cdot\frac{n}{k}} bound on PROP'.

In order to complete the proof of the main theorem, we are left to prove Theorem 14. This is done in the following section.

Proof of Theorem 14

The overall assignment is simply a combination of partial assignments obtained for each group individually using the following theorem.

Theorem 16. Let C′′>0C''>0 be a real, and let n>0n>0, n1>0n_1>0, k>0k>0 be integers, such that n≥n1≥n/kn\ge n_1\ge n/k. Let A∈[0,1]ns×mA\in[0,1]^{n_s\times m} be agents valuations of objects for the ss-th group with ns≤n1n_s\le n_1 agents. Then, there exists a partial assignment of goods v∈{0,1}mv\in\{0,1\}^m, such that the following properties are satisfied. Let l=⟨1,v⟩l=\langle1,v\rangle denote the number of assigned goods, and let Υ:=Ξ(n1⋅kn,n1⋅kn)\Upsilon:=\Xi\left(\frac{n_1\cdot k}{n},\frac{n_1\cdot k}{n}\right). Let matrix A′∈[0,1]ns×(m−l)A'\in[0,1]^{n_s\times(m-l)} be obtained from AA by removing columns that correspond to the assigned items (ones in vv). Then,

  1. l≤C′′⋅ns⋅Υ⋅knl\le C''\cdot n_s\cdot\sqrt{\Upsilon\cdot\frac{k}{n}}

  1. There is a permutation of agents p1,…,pnsp_1,\ldots,p_{n_s} and values n1≥Δ1≥…≥Δns≥n/(Υ⋅k)n_1\ge\Delta_1\ge\ldots\ge\Delta_{n_s}\ge n/(\Upsilon\cdot k), such that the matrix obtained from rearranging rows of A′A' according to the permutation pp

B=[Ap1′Ap2′⋮Apns′](4)B = \begin{bmatrix} A'_{p_1}\\ A'_{p_2}\\ \vdots\\ A'_{p_{n_s}} \end{bmatrix} \tag*{(4)}

is (Δp1,Δp2,…,Δpns)(\Delta_{p_1}, \Delta_{p_2}, \ldots, \Delta_{p_{n_s}})-structured and for every j∈[ns]j \in[n_s],

max⁡(⟨Bj,v⟩,C′′⋅γ⋅nk)≥C′′⋅γ⋅Δj.\max\left(\langle B_j,v\rangle,C''\cdot\sqrt{\gamma\cdot\frac{n}{k}}\right)\ge C''\cdot\gamma\cdot\sqrt{\Delta_j}.

Let us show how to combine these partial assignments and prove Theorem 14.

Proof of Theorem 14 via Theorem 16. The procedure has kk iterations, where at the ii-th iteration we select goods for the ii-th group.

  • As an input, we get the matrices A^1(i),…,A^k(i),B1(i),…,Bi−1(i)\hat{A}^{1}(i),\ldots,\hat{A}^{k}(i),B^{1}(i),\ldots,B^{i-1}(i). Initially (iteration 1), we have A^j(1)=Aj\hat{A}^{j}(1)=A^{j}.

  • We apply Theorem 16 to the ii-th group with the assignment matrix A^i(i)\hat{A}^{i}(i) playing the role of AA and with s=is=i to obtain a partial assignment vi:=vv_i:=v, a permutation pi:=pp^{i}:=p, and a matrix Bi(i):=BB^{i}(i):=B.

  • We then obtain A^1(i+1),…,A^k(i+1),B1(i+1),…,Bi(i+1)\hat{A}^{1}(i+1),\ldots,\hat{A}^{k}(i+1),B^{1}(i+1),\ldots,B^{i}(i+1) by removing the columns that correspond to the items that were gifted at this step from the matrices A^1(i),…,A^k(i),B1(i),…,Bi(i)\hat{A}^{1}(i),\ldots,\hat{A}^{k}(i),B^{1}(i),\ldots,B^{i}(i).

  • After the last step kk, we output A^j:=A^j(k+1)\hat{A}^{j}:=\hat{A}^{j}(k+1) and Bj=Bj(k+1)B^{j}=B^{j}(k+1).

Note that in total the procedure assigns at most C′′⋅(n1+…+nk)⋅γ⋅k/n=4⋅γ⋅n⋅kC''\cdot(n_1+\ldots+n_k)\cdot\sqrt{\gamma\cdot k/n}=4\cdot\sqrt{\gamma\cdot n\cdot k} objects. Moreover, for every ii, the matrix BiB^{i} is (Δpjii)(\Delta^{i}_{p^{i}_{j}}) structured, since it was (Δpjii)(\Delta^{i}_{p^{i}_{j}})-structured at step ii and this property is hereditary w.r.t. removing columns. At last, for every i∈[k],j∈[ni]i\in[k],j\in[n_i], we have ⟨Bji,vi⟩=⟨Bji(i),vi⟩\langle B^{i}_{j},v_i\rangle=\langle B^{i}_{j}(i),v_i\rangle, and

max⁡(⟨Bji,vi⟩,C′′⋅γ⋅nk)≥C′′⋅γ⋅Δji.\max\left(\langle B^{i}_{j},v_i\rangle,C''\cdot\sqrt{\gamma\cdot\frac{n}{k}}\right)\ge C''\cdot\gamma\cdot\sqrt{\Delta^{i}_{j}}.

since it was guaranteed by Theorem 16.

The rest of the section is devoted to the proof of Theorem 16. We begin with the description of the following algorithm that iteratively constructs the partial assignment vv. See Algorithm 1 for the pseudocode.

Gifting algorithm. The algorithm starts with an empty assignment v:=0mv:=0^{m} and sequentially selects an object that is given to the group. During its execution, the algorithm also maintains a set SS of active agents. Initially, every agent is active, i.e. S={1,2,…,ns}S=\{1,2,\ldots,n_s\}.

On each iteration, the algorithm finds a column j∈[m]j\in[m], such that vj=0v_j=0 and the sum of active agents valuations of the jj-th object is maximized: ∑i∈SAi,j→max⁡\sum_{i\in S}A_{i,j}\to\max. Let λ:=∑i∈SAi,j\lambda:=\sum_{i\in S}A_{i,j} denote the column sum over active agents. We iterate over active agents i∈Si\in S and look if ⟨Ai,v⟩≥C′′⋅γ⋅λ\langle A_i,v\rangle\ge C''\cdot\gamma\cdot\sqrt{\lambda}. If this holds, then we consider the ii-th agent to already be fulfilled and remove him from the set of active agents S←S∖{i}S\leftarrow S\setminus\{i\}. If at least one agent is removed after the filtration, we recalculate jj and λ\lambda that maximizes the column sum over the modified set of active agents and perform the filtration again.

After that, at the end of the iteration, if λ<1γ⋅nk\lambda<\frac{1}{\gamma}\cdot\frac{n}{k} or S=∅S=\varnothing, then the algorithm terminates. Otherwise, we give the jj-th item to the group: vj←1v_j\leftarrow1 and continue with the next iteration.

Algorithm 1 Gifting Algorithm

1: v←0mv \leftarrow0^{m} ▷ Initially no object is gifted 2: S←{1,2,…,ns}S \leftarrow\{1,2,\ldots,n_{s}\} ▷ Initially every agent is active 3: for t=1,2,…t=1,2,\ldots do 4: j←arg⁡max⁡j:vj=0∑i∈SAi,jj \leftarrow\arg\max_{j:v_{j}=0}\sum_{i\in S} A_{i,j} ▷ Find item most desirable by the active agents 5: λ←∑i∈SAi,j\lambda\leftarrow\sum_{i\in S} A_{i,j} 6: atLeastOneRemoved ←\leftarrow False 7: for i∈Si\in S do 8: if ⟨Ai,v⟩≥C′′⋅Υ⋅λ\langle A_{i},v\rangle\ge C''\cdot\Upsilon\cdot\sqrt{\lambda} then 9: S←S∖{i}S\leftarrow S\setminus\{i\} ▷ Make the ii-th agent inactive 10: atLeastOneRemoved ←\leftarrow True 11: end if 12: end for 13: if atLeastOneRemoved then 14: goto Step 4 ▷ Update jj and λ\lambda and repeat filtration 15: end if 16: if λ<1Υ⋅nk\lambda<\frac{1}{\Upsilon}\cdot\frac{n}{k} or S=∅S=\varnothing then 17: terminate and return vv ▷ λ\lambda is small or all agents are fulfilled 18: else 19: vj←1v_{j}\leftarrow1 ▷ Give the jj-th object to the group 20: end if 21: end for

Denote T=∥v∥1T=\lVert v\rVert_{1} to be the number of items that the algorithm has gifted. Note that the algorithm made T+1T+1 iterations in total. The TT-th iteration was the last when an item was added to the group. For t=1…T+1t=1\ldots T+1, denote StS_{t} to be the set of active agents at the end of the tt-th iteration, and denote by λt\lambda_{t} the column sum over active agents at the end of the tt-th iteration.

Let us note that λT≥1Υ⋅nk\lambda_{T}\ge\frac{1}{\Upsilon}\cdot\frac{n}{k} since the algorithm did not terminate at TT-th iteration. Moreover, observe that λt\lambda_{t} is monotone in tt, λ1≥λ2≥λ3≥…≥λT\lambda_{1}\ge\lambda_{2}\ge\lambda_{3}\ge\ldots\ge\lambda_{T}, since the space of active agents and columns over which we maximize and take the sum only shrinks.

Let tit_{i} be the iteration when the ii-th agent is removed from the active set. If the agent ii is never removed from the active set, we set ti:=T+1t_{i}:=T+1.

Let us show that ∥v∥1\lVert v\rVert_{1} is bounded.

Claim 17. The algorithm described above assigns no more than C′′⋅ns⋅Υ⋅k/nC''\cdot n_{s}\cdot\sqrt{\Upsilon}\cdot k/n objects to the group. That is, ∥v∥1≤C′′⋅ns⋅Υ⋅k/n\lVert v\rVert_{1}\le C''\cdot n_{s}\cdot\sqrt{\Upsilon}\cdot k/n.

Proof. For t=1…Tt=1\ldots T, let us denote by jtj_{t} the object given to the group on the tt-th iteration. Let us calculate the following sum in two different ways:

∑t=1T∑i∈StAi,jtλt.\sum_{t=1}^{T}\sum_{i\in S_{t}}\frac{A_{i,j_{t}}}{\lambda_{t}}.

Firstly, note that for every t=1…Tt=1\ldots T, we have ∑i∈StAi,jt=λt\sum_{i\in S_{t}}A_{i,j_{t}}=\lambda_{t} by the definition of λt\lambda_{t}. Therefore, on the one hand

∑t=1T∑i∈StAi,jtλt=∑t=1Tλtλt=T.\sum_{t=1}^{T}\sum_{i\in S_{t}}\frac{A_{i,j_{t}}}{\lambda_{t}} =\sum_{t=1}^{T}\frac{\lambda_{t}}{\lambda_{t}} =T.

On the other hand, we can swap the summations to sum over rows first:

∑t=1T∑i∈StAi,jtλt=∑i=1ns∑t=1ti−1Ai,jtλt≤∑i=1ns1λti−1∑t=1ti−1Ai,jt≤∑i=1ns1λti−1⋅(C′′⋅Υ⋅λti−1)=∑i=1ns1λti−1C′′⋅Υ≤∑i=1ns1λTC′′⋅Υ≤C′′⋅ns⋅Υ⋅kn(5)\begin{aligned} \sum_{t=1}^{T}\sum_{i\in S_t}\frac{A_{i,j_t}}{\lambda_t} = \sum_{i=1}^{n_s}\sum_{t=1}^{t_i-1}\frac{A_{i,j_t}}{\lambda_t} \\ &\le \sum_{i=1}^{n_s}\frac{1}{\lambda_{t_i-1}}\sum_{t=1}^{t_i-1}A_{i,j_t} \\ &\le \sum_{i=1}^{n_s}\frac{1}{\lambda_{t_i-1}}\cdot\left(C^{\prime\prime}\cdot\Upsilon\cdot\sqrt{\lambda_{t_i-1}}\right) \\ &= \sum_{i=1}^{n_s}\frac{1}{\sqrt{\lambda_{t_i-1}}}C^{\prime\prime}\cdot\Upsilon\\ &\le \sum_{i=1}^{n_s}\frac{1}{\sqrt{\lambda_T}}C^{\prime\prime}\cdot\Upsilon\\ &\le C^{\prime\prime}\cdot n_s\cdot\sqrt{\Upsilon\cdot\frac{k}{n}} \tag*{(5)} \end{aligned}

Therefore, we obtain the desired bound

∥v∥1=T≤C′′⋅ns⋅Υ⋅kn.\|v\|_1=T\le C^{\prime\prime}\cdot n_s\cdot\sqrt{\Upsilon\cdot\frac{k}{n}}.

To conclude the proof of Theorem 16, let us define the matrix BB, the permutation of agents pp, and the sequence Δi\Delta_i.

We define Δi\Delta_i to be the value of λ\lambda right at the moment when the ii-th agent is removed from the set SS. If the agent is never removed from the set SS, we set Δi:=max⁡(λT+1,1Υ⋅nk)\Delta_i:=\max\left(\lambda_{T+1},\frac{1}{\Upsilon}\cdot\frac{n}{k}\right). Let p1,p2,…,pnsp_1,p_2,\ldots,p_{n_s} be the sequence at which agents are removed from the active set (p1p_1 is the first agent that is removed, p2p_2 is the second, and so on). The agents that are not removed from the active set by the end of the algorithm are placed at the end of the sequence pp in any order. Observe that n1≥Δp1≥Δp2≥…≥Δpns≥n/(k⋅Υ)n_1\ge\Delta_{p_1}\ge\Delta_{p_2}\ge\ldots\ge\Delta_{p_{n_s}}\ge n/(k\cdot\Upsilon) since λ\lambda can only decrease over time.

Let A′A' be matrix obtained from AA by removing all columns corresponding to the gifted items (ones in vv). Define the matrix BB by arranging the rows of matrix A′A' in the following sequence

B=[Ap1′Ap2′⋮Apns′]B= \begin{bmatrix} A'_{p_1}\\ A'_{p_2}\\ \vdots\\ A'_{p_{n_s}} \end{bmatrix}

Let us show that the matrix BB is top-down structured. Claim 18. The matrix BB is (Δp1,Δp2,…,Δpns)(\Delta_{p_1},\Delta_{p_2},\ldots,\Delta_{p_{n_s}})-structured.

Proof. Consider the ii-th row of the matrix BB. We consider two cases whether the pip_i-th agent was removed from the set SS or not.

If the pip_i-th agent is never removed from the set SS, then pi,pi+1,…,pns∈ST+1p_i,p_{i+1},\ldots,p_{n_s} \in S_{T+1}. By definition, λT+1\lambda_{T+1} is the maximum column sum over rows corresponding to agents from ST+1S_{T+1}. Therefore, Δpi=max⁡(λT+1,nr)\Delta_{p_i}=\max\left(\lambda_{T+1},\frac{n}{r}\right) is also an upper bound on the column sum of the submatrix Bi,…,nsB_{i,\ldots,n_s}.

Now, suppose the pip_i-th agent is removed from the set SS at some point. Consider the moment just before pip_i-th agent is removed from SS. The agents pi,pi+1,…,pnsp_i,p_{i+1},\ldots,p_{n_s} all lie in the set SS. We also have that the current value of λ\lambda is an upper bound on the maximum column sum over rows corresponding to agents from SS. Therefore, Δpi=λ\Delta_{p_i}=\lambda is also an upper bound on the column sum of the submatrix Bi,…,nsB_{i,\ldots,n_s}.

Thus, the matrix BB is (Δp1,…,Δpns)(\Delta_{p_1},\ldots,\Delta_{p_{n_s}})-structured.

Finally, let us show that every agent values gifted objects at the desired threshold.

Claim 19. For every agent i∈[ns]i \in[n_s], it holds that

max⁡(⟨Ai,v⟩,C′′⋅Υ⋅nk)≥C′′⋅Υ⋅Δi.(6)\max\left(\langle A_i,v\rangle,C''\cdot\sqrt{\Upsilon\cdot\frac{n}{k}}\right)\ge C''\cdot\Upsilon\cdot\sqrt{\Delta_i}. \tag*{(6)}

Proof. Let us condition on whether the ii-th agent was removed from the set SS or not.

If it was not removed, then Δi=λT+1<1Υ⋅nk\Delta_i=\lambda_{T+1}<\frac{1}{\Upsilon}\cdot\frac{n}{k}. Therefore, Υ⋅nk≥Υ⋅Δi\sqrt{\Upsilon\cdot\frac{n}{k}}\ge\Upsilon\cdot\sqrt{\Delta_i}.

If it was removed, then it could happen only if ⟨Ai,v⟩≥C′′⋅Υ⋅λ\langle A_i,v\rangle\ge C''\cdot\Upsilon\cdot\sqrt{\lambda} was true at that moment. Since Δi=λ\Delta_i=\lambda, we get the desired bound.

Therefore, in both cases we can conclude (5).

The claims 18 and 19 together conclude the proof of Theorem 16.

Lower Bound in Theorem 4

In this section, we prove that disc⁡>k=Ω(n/k)\operatorname{disc}_{>}^{k}=\Omega\left(\sqrt{n/k}\right). We start with a trivial lower bound on osasymdisc⁡>max⁡\operatorname{osasymdisc}_{>}^{\max}.

Lemma 20.

disc⁡>k(n1,…,nk)≥k−1k.\operatorname{disc}_{>}^{k}(n_1,\ldots,n_k)\ge\frac{k-1}{k}.

Proof. Let k≥2k\ge2 be the number of groups. Let n1,…,nk>0n_1,\ldots,n_k>0 be group sizes, and let n=n1+…+nkn=n_1+\ldots+n_k be the total number of agents.

Suppose there are k−1k-1 objects, and each of them is valued by 1 by every agent. Then, in any allocation at least one group receives no objects. The fair share for this group is k−1k\frac{k-1}{k}. Therefore, agents from this group are at least k−1k\frac{k-1}{k} unsatisfied.

To prove our lower bound, we rely on the following lower bound on weighted discrepancy.

Theorem 21 (Theorem 3.2 from [17]). For any integers r>0r>0 and k>0k>0, there exists a matrix A∈[0,1]r×(r⋅k)A\in[0,1]^{r\times(r\cdot k)}, such that for any v∈{0,1}mv\in\{0,1\}^{m}, it holds

∥A(v−1k)∥∞≥r−116.(7)\left\|A\left(v-\frac{1}{k}\right)\right\|_{\infty}\ge\sqrt{r}-\frac{1}{16}. \tag*{(7)}

The theorem guarantees that the maximum of the absolute value of deviation of ⟨Ai,(v−1/k)⟩\langle A_i,(v-1/k)\rangle from 0 is at least r−1/16\sqrt{r-1}/16. If the deviation is downwards then we immediately infer that some agent is unsatisfied. However, an upward deviation is only welcomed by an agent and thus does not immediately imply any lower bounds for group division. To deal with this, we introduce agents with complementary valuations. This converts an upward deviation into a downward one. Similarly to [17], we derive the following result.

Lemma 22. For any r≥1r \ge1, there exists m0m_0, such that for any m≥m0m \ge m_0 there exists a matrix B∈{0,1}2r×mB \in\{0,1\}^{2r \times m} that satisfies the following. For any v∈{0,1}mv \in\{0,1\}^m, there exists ii, such that

⟨Bi,v⟩≤⟨Bi,1/k⟩−r−116+max⁡(0,∥v∥1−m/k)(8)\langle B_i,v\rangle\le\langle B_i,1/k\rangle- \frac{\sqrt{r-1}}{16} + \max(0,\|v\|_1-m/k) \tag*{(8)}

Proof. Let A∈{0,1}r×m0A \in\{0,1\}^{r \times m_0} be a matrix given by Theorem 21. For m≥m0m \ge m_0, let a matrix A′∈{0,1}r×mA' \in\{0,1\}^{r \times m} be obtained from AA by padding extra columns with zeroes. Note that the matrix A′A' still satisfies the property from Theorem 21 since zero columns do not change the product.

Define CC to be the matrix of rows complementary to A′A': Ci,j=1−Ai,j′C_{i,j}=1-A'_{i,j} for all i,ji,j. Let BB be a vertical concatenation of matrices A′A' and CC:

B=[A′C].B = \begin{bmatrix} A' \\ C \end{bmatrix}.

Consider a vector v∈{0,1}mv \in\{0,1\}^m and let us show that the matrix BB satisfies the required property. By Theorem 21, for any v∈{0,1}mv \in\{0,1\}^m, there exists ii, such that

∣⟨Ai′,v−1/k⟩∣≥r−116.\left|\langle A'_i,v-1/k\rangle\right| \ge\frac{\sqrt{r-1}}{16}.

Let us consider the two cases how the absolute value can be expanded.

If −⟨Ai′,v−1/k⟩≥r−1/16-\langle A'_i,v-1/k\rangle\ge\sqrt{r-1}/16 then ⟨Ai′,v⟩≤⟨Ai′,1/k⟩−r−1/16\langle A'_i,v\rangle\le\langle A'_i,1/k\rangle-\sqrt{r-1}/16, which implies (8).

If ⟨Ai′,v−1/k⟩≥r−1/16\langle A'_i,v-1/k\rangle\ge\sqrt{r-1}/16, then

⟨1−Ai′,v−1/k⟩≤−r−116+⟨1,vv−1/k⟩=−r−116+∥v∥1−m/k.\langle1-A'_i,v-1/k\rangle\le-\frac{\sqrt{r-1}}{16}+\langle1,vv-1/k\rangle= -\frac{\sqrt{r-1}}{16}+\|v\|_1-m/k.

Since 1−Ai′1-A'_i is a row of BB, we get the desired conclusion.

Proof of the Lower Bound

Lemma 20 gives a lower bound disc⁡k(n1,…,nk)≥12\operatorname{disc}^k(n_1,\ldots,n_k) \ge\frac{1}{2}. This allows us to assume that n≥214⋅kn \ge2^{14}\cdot k. Denote r=⌊n/(4⋅k)⌋≥212r=\lfloor n/(4\cdot k)\rfloor\ge2^{12}, and let m0m_0 be the integer from Lemma 22 when applied with rr. Let M=n⋅m0M=n\cdot m_0 be the number of objects. We design agents valuations for the ii-th group as follows.

  • If ni<2rn_i < 2r, then we set all of the agents from the ii-th group to valuate any object by 1.

  • If ni≥2rn_i \ge2r, then we divide all agents into subgroups of size 2r2r with <2r<2r agents left in no subgroup. In total, this gives si=⌊ni/(2ri)⌋s_i=\lfloor n_i/(2r_i)\rfloor subgroups. Let mi=⌊M/si⌋≥m0m_i=\lfloor M/s_i\rfloor\ge m_0, and let the matrix B∈{0,1}2r×miB \in\{0,1\}^{2r\times m_i} be given by Lemma 22. For the jj-th subgroup we set the valuations of the objects (j−1)⋅mi+1,…,j⋅mi(j-1)\cdot m_i+1,\ldots,j\cdot m_i according to the matrix BB. All other objects are valuated by 0 by the agents of this subgroup. See Figure 1 for an illustration of the chosen valuations.

The valuations of objects for each group

Figure 1. The valuations of objects for each group. The big zeros correspond to a matrix of 0s. The valuations of groups with fewer than 2r2r agents are set to 1. If a group ii has at least 2r2r agents, then it is partitioned into subgroups of size 2r2r where each subgroup values some subset of objects according to matrix BB. The remaining <2r< 2r agents get arbitrary valuations (for concreteness, we set them to 0).

Let us give some intuition behind the construction. Note that rr is essentially an average size of a group. Thus, large groups get valuations in the spirit of the construction from [17], while small groups evaluate everything to 1. The evaluation of small groups thus forces that roughly a correct proportion of goods go to the large groups. There, the goods are partitioned between subgroups. This way, we are able to apply the bound of Lemma 22 inside a subgroup, rather than inside the whole group.

We go on to the proof. Arguing indirectly, assume that there is an allocation v1,…,vkv_1,\ldots,v_k, such that every agent is at most 2−15⋅n/k2^{-15} \cdot\sqrt{n/k} unsatisfied.

Define di=∑j=1Mvi,j−Mkd_i = \sum_{j=1}^{M} v_{i,j} - \frac{M}{k} to be the difference between the number of allocated goods and the fair share for the ii-th group. Note that ∑i=1kdi=0\sum_{i=1}^{k} d_i = 0. Let us estimate ∑ni≥2rdi\sum_{n_i \ge2r} d_i and ∑ni<2rdi\sum_{n_i < 2r} d_i separately.

Consider a group ii for which ni≥2rn_i \ge2r. We have si=⌊ni/2r⌋s_i = \lfloor n_i/2r \rfloor subgroups in this group. The jj-th subgroup has non-zero valuations for the segment of objects (j−1)⋅mi+1,…,j⋅mi(j-1)\cdot m_i+1,\ldots,j\cdot m_i. By the pigeonhole principle, there is a subgroup jj, such that ∑l=(j−1)⋅mi+1j⋅mivi,l−Mk⋅si≤disi\sum_{l=(j-1)\cdot m_i+1}^{j\cdot m_i} v_{i,l}-\frac{M}{k\cdot s_i}\le\frac{d_i}{s_i}. By (6), there is an agent in this subgroup that is at least r−116−max⁡(0,di/si+1/k)\frac{\sqrt{r-1}}{16}-\max(0,d_i/s_i+1/k) unsatisfied (here we used ∣mi−M/si∣≤1|m_i-M/s_i|\le1). If di/si<r−132d_i/s_i<\frac{\sqrt{r-1}}{32}, then we get a contradiction since this agent’s unsatisfaction is at least

r−116−r−132−12≥r−1−1632>2−15⋅nk.\frac{\sqrt{r-1}}{16}-\frac{\sqrt{r-1}}{32}-\frac{1}{2}\ge\frac{\sqrt{r-1}-16}{32}>2^{-15}\cdot\sqrt{\frac{n}{k}}.

Thus, we have di≥si⋅r−1/32≥2−11⋅si⋅n/kd_i\ge s_i\cdot\sqrt{r-1}/32\ge2^{-11}\cdot s_i\cdot\sqrt{n/k}. Therefore,

∑ni≥2rdi≥∑ni≥2r2−13⋅nir⋅nk≥∑ni≥2r2−13⋅ni⋅kn≥2−13⋅kn(∑i=1kni−∑ni<2rni)≥2−13⋅kn(∑i=1kni−2r⋅k)≥2−13⋅kn(n−n2)≥2−14⋅n⋅k.\begin{aligned} \sum_{n_i\ge2r}d_i&\ge\sum_{n_i\ge2r}2^{-13}\cdot\frac{n_i}{r}\cdot\sqrt{\frac{n}{k}}\\ &\ge\sum_{n_i\ge2r}2^{-13}\cdot n_i\cdot\sqrt{\frac{k}{n}}\\ &\ge2^{-13}\cdot\sqrt{\frac{k}{n}}\left(\sum_{i=1}^{k}n_i-\sum_{n_i<2r}n_i\right)\\ &\ge2^{-13}\cdot\sqrt{\frac{k}{n}}\left(\sum_{i=1}^{k}n_i-2r\cdot k\right)\\ &\ge2^{-13}\cdot\sqrt{\frac{k}{n}}\left(n-\frac{n}{2}\right)\\ &\ge2^{-14}\cdot\sqrt{n\cdot k}. \end{aligned}

Consider any agent from group ii for which ni<2rn_i<2r. Recall that it evaluates every good to 11, and thus did_i is exactly the measure of how unsatisfied he is. Thus, did_i cannot be less than −2−15⋅n/k-2^{-15}\cdot\sqrt{n/k}. Therefore, we have

∑ni<2rdi≥−2−15⋅n⋅k.\sum_{n_i<2r}d_i\ge-2^{-15}\cdot\sqrt{n\cdot k}.

Combining two bounds, we have

∑i=1kdi≥(2−14−2−15)⋅n⋅k>0,\sum_{i=1}^{k}d_i\ge(2^{-14}-2^{-15})\cdot\sqrt{n\cdot k}>0,

reaching a contradiction.

References

  1. [1]Georgios Amanatidis, Haris Aziz, Georgios Birmpas, Aris Filos-Ratsikas, Bo Li, Hervé Moulin, Alexandros A Voudouris, and Xiaowei Wu. Fair division of indivisible goods: Recent progress and open questions. Artificial Intelligence, 322:103965, 2023.
  2. [2]Wojciech Banaszczyk. Balancing vectors and gaussian measures of n-dimensional convex bodies. Random Struct. Algorithms, 12(4):351–360, 1998.
  3. [3]Nikhil Bansal. Constructive algorithms for discrepancy minimization. In 2010 IEEE 51st Annual Symposium on Foundations of Computer Science, pages 3–10. IEEE, 2010.
  4. [4]Nikhil Bansal and Haotian Jiang. Decoupling via affine spectral-independence: Beck-fiala and komlós bounds beyond banaszczyk. In Aditya Bhaskara and Artur Czumaj, editors, Proceedings of the 58th Annual ACM Symposium on Theory of Computing, STOC 2026, Salt Lake City, UT, USA, June 22-26, 2026, pages 432–442. ACM, 2026. doi: 10.1145/3798129.3800762. URL https://doi.org/10.1145/3798129.3800762.
  5. [5]Nikhil Bansal, Daniel Dadush, Shashwat Garg, and Shachar Lovett. The gram-schmidt walk: a cure for the banaszczyk blues. In Proceedings of the 50th annual acm sigact symposium on theory of computing, pages 587–597, 2018.
  6. [6]Nikhil Bansal, Daniel Dadush, and Shashwat Garg. An algorithm for komlós conjecture matching banaszczyk’s bound. SIAM Journal on Computing, 48(2):534–553, 2019.
  7. [7]József Beck and Tibor Fiala. "integer-making" theorems. Discret. Appl. Math., 3(1):1–8, 1981. doi: 10.1016/0166-218X(81)90022-6. URL https://doi.org/10.1016/0166-218X(81)90022-6.
  8. [8]Ioannis Caragiannis, Kasper Green Larsen, and Sudarshan Shyam. A new lower bound for multicolor discrepancy with applications to fair division. In International Symposium on Algorithmic Game Theory, pages 228–246. Springer, 2025.DOI
  9. [9]Bernard Chazelle. The discrepancy method: Randomness and complexity. Cambridge University Press, 2000.
  10. [10]Benjamin Doerr. Discrepancy in different numbers of colors. Discrete Mathematics, 250(1-3):63–70, 2002.DOI
  11. [11]Benjamin Doerr and Anand Srivastav. Multicolour discrepancies. Combinatorics, Probability and Computing, 12(4):365–399, 2003.DOI
  12. [12]Max Dupré la Tour. Bad news for couples: Tight lower bounds for fair division of indivisible items. arXiv e-prints, pages arXiv–2601, 2026.arxiv.org/abs/2601.01012
  13. [13]Paul Gölz and Hannane Yaghoubi zade. Fair division among couples and small groups. In Proceedings of the AAAI Conference on Artificial Intelligence, volume 40, pages 16963–16970, 2026.
  14. [14]Maria Kyropoulou, Warut Suksompong, and Alexandros A Voudouris. Almost envy-freeness in group resource allocation. Theoretical Computer Science, 841:110–123, 2020.
  15. [15]Max Dupre la Tour and Kaito Fujii. Discrepancy and fair division for non-additive valuations. arXiv preprint arXiv:2509.16802, 2025.arxiv.org/abs/2509.16802
  16. [16]Shachar Lovett and Raghu Meka. Constructive discrepancy minimization by walking on the edges. In 53rd Annual IEEE Symposium on Foundations of Computer Science, FOCS 2012, New Brunswick, NJ, USA, October 20-23, 2012, pages 61–67. IEEE Computer Society, 2012. doi: 10.1109/FOCS.2012.23. URL https://doi.org/10.1109/FOCS.2012.23.
  17. [17]Pasin Manurangsi and Raghu Meka. Tight lower bound for multicolor discrepancy. In Sepehr Assadi and Eva Rotenberg, editors, 2026 Symposium on Simplicity in Algorithms, SOSA 2026, Vancouver, BC, Canada, January 12-14, 2026, pages 266–274. SIAM, 2026. doi: 10.1137/1.9781611978964.20. URL https://doi.org/10.1137/1.9781611978964.20.
  18. [18]Pasin Manurangsi and Warut Suksompong. Asymptotic existence of fair divisions for groups. Mathematical Social Sciences, 89:100–108, 2017.
  19. [19]Pasin Manurangsi and Warut Suksompong. Almost envy-freeness for groups: Improved bounds via discrepancy theory. Theoretical Computer Science, 930:179–195, 2022.
  20. [20]Pasin Manurangsi and Warut Suksompong. Ordinal maximin guarantees for group fair division. Theoretical Computer Science, 1036:115151, 2025.
  21. [21]Jiri Matousek. Geometric discrepancy: An illustrated guide, volume 18. Springer Science & Business Media, 1999.DOI
  22. [22]Thomas Rothvoss. Constructive discrepancy minimization for convex sets. SIAM Journal on Computing, 46(1):224–234, 2017.
  23. [23]Erel Segal-Halevi and Shmuel Nitzan. Fair cake-cutting among families. Social Choice and Welfare, 53(4):709–740, 2019.
  24. [24]Erel Segal-Halevi and Warut Suksompong. Democratic fair allocation of indivisible goods. Artificial Intelligence, 277:103167, 2019.
  25. [25]Erel Segal-Halevi and Warut Suksompong. How to cut a cake fairly: A generalization to groups. The American Mathematical Monthly, 128(1):79–83, 2020.
  26. [26]Joel Spencer. Six standard deviations suffice. Transactions of the American mathematical society, 289(2):679–706, 1985.
  27. [27]Warut Suksompong. Approximate maximin shares for groups of agents. Mathematical Social Sciences, 92:40–47, 2018.DOI

Paper details

vlnc:W_92pnjk9a
v1 · latestHistory (1)

published ·

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Computer Science

Not recorded for this paper · View statement →

Looking for related papers…
Contents