Fair Cake Division Under Monotone Likelihood Ratios

Siddharth Barman, Nidhi Rathi

Introduction

Cake division is a quintessential model in the study of fair division. This setup captures the allocation of a divisible resource (metaphorically, the cake) among agents with equal entitlements, but distinct preferences. Over the past several decades, a significant body of work in mathematics, economics, and computer science has been devoted to cake cutting; see [BT96, RW98, Pro15] for excellent expositions and motivating applications (e.g., border negotiations and divorce settlements) of this framework.

Some of the central solution concepts and axiomatic characterizations in the fair-division literature stem from the cake-cutting context [Mou04]. Indeed, the work of Steinhaus, Banach, and Knaster [Ste48]—which lays the mathematical foundations of fair division—addresses cake division. The notion of envy-freeness was also mathematically formalized in this setup [GS58, Fol67]. This well-studied notion deems a cake division to be fair if every agent prefers the piece assigned to her over that of any other agent, i.e., if no agent is envious of others.

This work focuses on a standard formulation of cake division in which every agent must receive a contiguous piece of the cake. That is, the goal is to partition the cake $intoexactlyinto exactlyndisjointintervals(connectedpieces)andassignthemamongthedisjoint intervals (connected pieces) and assign them among thenparticipatingagents.Thisconnectivityrequirementisnaturallymotivatedbysettingsinwhichacontiguouspartoftheresourceneedstobeallocatedtoeveryagent[BT96];consider,e.g.,divisionofland,transmissionspectrum,orprocessingtimeonamachine.Notethatapartitionofthecakeparticipating agents. This connectivity requirement is naturally motivated by settings in which a contiguous part of the resource needs to be allocated to every agent [BT96]; consider, e.g., division of land, transmission spectrum, or processing time on a machine. Note that a partition of the cakeintointervalsinto intervalsI_{1},I_{2},\ldots,I_{n}—whereininterval—wherein intervalI_{i}isassignedtoagentis assigned to agenti\in[n]—issaidtobeenvy−freeiff—is said to be envy-free iffv_{i}(I_{i})\geq v_{i}(I_{j})(i.e.,iff(i.e., iff\int_{I_{i}}f_{i}\geq\int_{I_{j}}f_{i})forallagents) for all agentsiandandj$.

The appeal of envy-freeness is substantiated by strong existential results: under mild assumptions, a contiguous envy-free cake division always exists [Str80, Sim80, Su99]. While these results are built upon interesting mathematical connections,For instance, the proof by Su [Su99] invokes Sperner’s lemma. they are, however, nonconstructive. In fact, Stromquist [Str08] has shown that there does not exist a finite-time algorithm for finding envy-free cake divisions with connected pieces; this result holds in a setup wherein the valuations are provided through an (adversarial) oracle. In addition, the work of Deng et al. [DQS12] establishes PPAD-hardness of finding envy-free cake divisions with contiguous pieces, under ordinal valuations.

Algorithms for envy-free cake division remain elusive even if we relinquish the contiguity requirement. It was not until the work of Brams and Taylor [BT95] that a bounded-time algorithm was obtained for noncontiguous envy-free cake division. In general, the best-known result for this problem is by Aziz and Mackenzie [AM16], who develop a hyper-exponential time algorithm for finding envy-free divisions with noncontiguous pieces.The problem of finding an approximate envy-free division (not necessarily with connected pieces) admits a fully-polynomial time approximation scheme [LMMS04].

In light of these algorithmic barriers, identification of computationally-tractable instances in the cake-cutting context stands as a meaningful direction of work. The current paper addresses this consideration and, in particular, identifies an encompassing property—called the monotone likelihood ratio property—which enables the development of efficient algorithms for fair cake-cutting (with connected pieces).

In the cake-division context, we will say that an ordered collection (fi)i∈[n](f_{i})_{i\in[n]} of value densities (of the nn agents) satisfies the monotone likelihood ratio property iff for each i∈[n−1]i\in[n-1], the likelihood ratio \nicefracfi+1(x)fi(x)\nicefrac{{f_{i+1}(x)}}{{f_{i}(x)}} is nondecreasing in x∈x\in. That is, the agents are indexed with the property that consecutive likelihood ratios bear MLRP. This property is transitive and, hence, in cake-division instances with MLRP, value densities fif_{i} and fjf_{j} satisfy MLRP for all i<ji<j.

Many distribution families are also known to bear MLRP [LM+01, CB02]. In particular, this property holds if all the value densities belong to any one of the following families: Gaussian distributions (with the same variance but different means), Poisson distributions, binomial distributions, and single-parameter exponentials; see Appendix C for details. Furthermore, it is known that linear translations of any log-concave function satisfy MLRP [SW14]. In particular, linear translations of the following (log-concave) distributions also satisfy this property: Laplace, uniform, multivariate Gaussian, gamma, beta, Subbotin, chi-square, Dirichlet, and logistic. Hence, the current work obtains novel results for many distribution families in a unified manner.

MLRP is a common assumption on agents’ utilities and type distributions in many economic contexts; see [Jew91] for a survey. As a stylized application of MLRP in cake division, consider a setting wherein each agent ii has a most preferred point μi\mu_{i} on the divisible resource (cake) and ii’s valuation density decreases as a Gaussian function (with a variance parameter that is common across the agents) of the distance from μi\mu_{i}. Indeed, the distance here can be geographical (as in case of land division), temporal (i.e., wait time), or it can be an abstract metric.

Considering similar single-peaked preferences, but with a linear drop in value densities, Wang and Wu [WW19] developed an efficient algorithm for noncontiguous cake division. Note that while linear densities bear MLRP, this property does not hold for piecewise linear densities. Hence our results do not directly address the setting considered in [WW19]. However, in absence of the contiguity requirement (as is the case in [WW19]) one can find a fair cake division by first partitioning the cake into intervals, in each of which the agents’ value densities are linear, and then applying the MLRP result separately.Recall that, in contrast to such a result, our focus is on finding cake divisions in which each agent receives a contiguous piece of the cake.

We focus on cake-division instances with MLRP and develop algorithmic results for almost all the standard notions of fairness and economic efficiency. Our algorithms only require oracle access to the valuations. In particular, the developed algorithms operate under the standard Robertson-Webb model [RW98], wherein we have access to the agents’ valuations through eval and cut queries; see Section 2 for details. MLRP implies that these cut and eval queries (functions) are λ\lambda-Lipschitz (Appendix A.3).

The time complexities of our algorithms depend polynomially on the the bit complexity of this Lipschitz constant λ≥1\lambda\geq 1. Such a runtime dependency on log⁡λ\log\lambda is unavoidable (Appendix D): there exist cake-division instances (with λ\lambda-Lipschitz cut and eval queries) wherein for all the agents the value of the cake is almost entirely concentrated in an interval LL of length 1/λ{1}/{\lambda}. Here, an envy-free cake division can be obtained only by finely partitioning LL among the agents. In particular, the cut points that induce an envy-free allocation (and, hence, correspond to the output of a fair-division algorithm) must be 1/λ1/\lambda close to each other, i.e., the bit complexity of the output has to be Ω(log⁡λ){\Omega}\left(\log\lambda\right). In fact, one can construct instances in which a contiguous envy-free division can be obtained only by cutting the cake at irrational points (Appendix D). Hence, in general, (and even under MLRP) one cannot expect an efficient algorithm that outputs an exact envy-free division, with contiguous pieces.Indeed, the bit complexity of a computationally-bounded algorithm is bounded as well. Therefore, when considering efficient algorithms for cake division, a precision loss in the output is inevitable. However, our algorithms ensure that this precision loss in value, η\eta, is arbitrarily small; specifically, the developed algorithms run in time O(poly(n,log⁡λ,log⁡1η))\mathcal{O}\left({\rm poly}\left(n,\log\lambda,\log\frac{1}{\eta}\right)\right) and, hence, the precision parameter η\eta can be driven exponentially close to zero in polynomial (in the bit complexity of η\eta) time. Note that this bit-precision issue is akin to the one faced in the convex-optimization problems (where again the optimal solutions can be irrational) and our runtime bound, with respect to the precision parameter η\eta, is analogous to the one obtained by the ellipsoid method [GLS12].

Our Results and Techniques: Next, we summarize our results for various notions of fairness and (economic) efficiency.

Envy-Freeness: We prove that, given a cake-division instance (in the Robertson-Webb query model) with MLRP, an envy-free allocation can be computed, up to an arbitrary precision, in time that is polynomial in the number of agents and log⁡λ\log\lambda; here λ\lambda is the Lipschitz constant of the Robertson-Webb (cut and eval) queries.

To establish this result, we define a class of divisions, referred to as ripple divisions (Definition 1), and prove that, under MLRP, every ripple division induces a contiguous envy-free cake division (Theorem 6). Specifically, a collection of points x0=0≤x1≤x2≤xn−1≤xn=1x_{0}=0\leq x_{1}\leq x_{2}\leq x_{n-1}\leq x_{n}=1 (in the cake $)issaidtoformarippledivisionofthecakeif,foreach) is said to form a ripple division of the cake if, for eachi\in[n-1],agent, agentiisindifferentbetweentheconsecutiveintervalsis indifferent between the consecutive intervals[x_{i-1},x_{i}]andand[x_{i},x_{i+1}],i.e.,, i.e.,v_{i}(x_{i-1},x_{i})=v_{i}(x_{i},x_{i+1}).Notethatarippledivisioninducesacontiguouscakedivision—byassigninginterval. Note that a ripple division induces a contiguous cake division—by assigning interval[x_{i-1},x_{i}]toagentto agenti—withthepropertythatagent—with the property that agentidoesnotenvyagentdoes not envy agenti+1.Thatis,inandofitself,arippledivisionmandatesabsenceofenvyonlybetweenconsecutiveagents,andnotbetweenallpairsofagents.Wewillshowthat,interestingly,underMLRPthisrelaxationsuffices–thecakedivisioninducedbyarippledivisionisguaranteedtobeenvyfree(Theorem6).RecallthattheagentsareindexedfollowingtheMLRPorder:foreach. That is, in and of itself, a ripple division mandates absence of envy only between consecutive agents, and not between all pairs of agents. We will show that, interestingly, under MLRP this relaxation suffices–the cake division induced by a ripple division is guaranteed to be envy free (Theorem 6). Recall that the agents are indexed following the MLRP order: for eachi\in[n-1],thelikelihoodratio, the likelihood ratiof_{i+1}/f_{i}isnondecreasing.Hence,allocatingis nondecreasing. Hence, allocating[x_{i-1},x_{i}]toagentto agenti\in[n]$ ensures that the intervals are assigned (left to right on the cake) in accordance with the MLRP order.

We establish the universal existence of ripple divisions through the intermediate value theorem, i.e., a one-dimensional fixed-point argument (Lemma 3). Since one can use binary search to find fixed points in the one-dimensional setting, this proof in fact leads to an algorithm for finding ripple divisions and, hence, envy-free divisions. Indeed, the notion of ripple divisions and their connection with envy-freeness, under MLRP, are two key contributions of this work.

Pareto Optimality: We show that in cake-division instances with MLRP, Pareto optimal cake divisions, with connected pieces, conform to the MLRP order (Lemma 5). This structural result implies that for maximizing welfare we can restrict attention to allocations wherein the intervals are assigned (left to right on the cake) in accordance with the MLRP order. Intuitively, this leads us to a welfare-maximizing algorithm—specifically, a dynamic program—that recursively finds optimal allocations for intervals placed at the left end of the cake (i.e., for intervals of the form [0,x][0,x]).

We also establish an extension of Weller’s theorem in the MLRP context. Weller’s theorem [Wel85] asserts that there always exists some cake division—though, not necessarily with connected pieces—which is both envy-free (fair) and Pareto optimal. While this theorem holds in general,Weller’s theorem applies even in the absence of MLRP. it does not guarantee that envy-freeness and Pareto optimality can be achieved together through contiguous cake divisions. We show that, by contrast, under MLRP every contiguous envy-free division is Pareto optimal (Theorem 2). Therefore, given a cake-division instance with MLRP, the allocation computed by our algorithm is not only envy-free but also Pareto optimal, up to an arbitrary precision.

Social Welfare: Social (utilitarian/Benthamite) welfare is a standard measure of collective value. For a cake division {I1,I2,…,In}\{I_{1},I_{2},\ldots,I_{n}\} it is defined to be the sum of the values that the division generates among the agents, ∑ivi(Ii)\sum_{i}v_{i}(I_{i}). Maximizing social welfare is a well-studied objective in resource-allocation contexts. In the cake-cutting setup, this maximization problem is known to be APX-hard under general valuations [ABKR19]. Complementarily, if the value densities bear MLRP, then we can find (up to an arbitrary precision) a social welfare maximizing division with connected pieces in O(poly(n,log⁡λ))\mathcal{O}\left({\rm poly}\left(n,\log\lambda\right)\right) time (Theorem 3). As mentioned previously, our algorithm for this problem is based on a dynamic program.

Egalitarian Welfare: The egalitarian (Rawlsian) welfare of a cake division {I1,…,In}\{I_{1},\ldots,I_{n}\} is defined as the value of the least well-off agent, i.e., min⁡i vi(Ii)\min_{i}\ v_{i}(I_{i}). From a welfarist perspective, maximizing this minimum value among cake divisions with connected pieces is an important fairness objective. However, no nontrivial approximation guarantees are known for this problem under general valuations; the work of Aumann et al. [ADH13] shows that maximizing egalitarian welfare across all contiguous cake divisions is APX-hard. Complementing this hardness result, we develop an algorithm that, under MLRP, maximizes egalitarian welfare (up to an arbitrary precision) and runs in O(poly(n,log⁡λ))\mathcal{O}\left({\rm poly}\left(n,\log\lambda\right)\right) time (Theorem 4).

Our algorithm for maximizing egalitarian welfare is based on a “moving-knife” procedure. This procedure, for a given a target value τ>0\tau>0, iteratively selects points x0=0,x1,x2,…,xn≤1x_{0}=0,x_{1},x_{2},\ldots,x_{n}\leq 1 such that the each interval [xi−1,xi][x_{i-1},x_{i}] is of value τ\tau to agent i∈[n]i\in[n]. Let τ∗\tau^{*} denote the optimal egalitarian welfare in the given cake-division instance. The useful observation here is that this moving-knife procedure will succeed for all τ≤τ∗\tau\leq\tau^{*}. This follows from the fact that here the intervals are assigned (left to right) in the MLRP orderNote that the agents are indexed accordingly. and this ordering is also satisfied by an egalitarian welfare maximizing (in particular, a Pareto optimal) division. Therefore, by performing a binary search with τ\tau, we can find a contiguous division with egalitarian welfare arbitrarily close to the optimal.

Nash Social Welfare: A balance between social and egalitarian welfare is obtained by considering the Nash social welfare [NJ50, KN79]. This welfare objective is defined as the geometric mean of the agents’ values. It is known that, in general, it is APX-hard to find a contiguous cake division that maximizes Nash social welfare [ABKR19]. Under MLRP, however, the problem of maximizing Nash social welfare admits a fully polynomial-time approximation scheme (Theorem 5). We obtain this result via a dynamic program that considers the agents in the MLRP order.

Additional Related Work: Recently, approximation algorithms—with both additive [HGS19] and multiplicative [ABKR19] approximation guarantees—have been developed for finding contiguous envy-free cake divisions. The work of Brânzei and Nisan [BN17] develops query complexity upper and lower bounds for computing approximately envy-free allocations. In contrast to these results, the current work focuses on cake-division instances with MLRP and shows that in such settings arbitrarily low envy can be achieved among the agents.

The work of Bei et al. [BCH+12] also studies contiguous cake division and provides computational results for maximizing social welfare subject to proportional fairness. Under this fairness constraint each agent i∈[n]i\in[n] must receive an interval of value at least 1/n1/n times ii’s total value for the cake. Bei et al. [BCH+12] show that, if the value densities are linear, then this problem admits a fully polynomial-time approximation scheme (FPTAS). We note that every pair of linear densities bear MLRP and, hence, such value-density functions fall within the purview of the current work. However, our algorithm is incomparable to the FPTAS of Bei et al. [BCH+12]–we focus on maximizing social welfare without the fairness constraints. Also, envy-freeness is not addressed in [BCH+12].

Another well-studied fairness notion is that of equitability. Specifically, a cake division {I1,I2,…,In}\{I_{1},I_{2},\ldots,I_{n}\} is said to be equitable iff all the agents derive the same value from the intervals assigned to them, vi(Ii)=vj(Ij)v_{i}(I_{i})=v_{j}(I_{j}) for all ii and jj [DS61, Alo87]. In other words, equitability ensures that all the agents are equally well-off. Cechlárová and Pillárová [CP12] consider the computation of equitable cake divisions with connected pieces. They showed that—given access to “reverse” cutting queries—such divisions can be efficiently computed, up to an arbitrary precision. We note that value densities that satisfy MLRP have, by definition, full support over the cake. In such a case, the reverse cutting queries can be simulated by standard (cut) queries in the Robertson-Webb model. Hence, under MLRP, strong algorithmic results hold for equitability as well.

Cake-division algorithms for specific classes of valuations have been studied in [CLPP11] and [KLP13]. The work of Kurokawa et al. [KLP13] provides a query-efficient algorithm for envy-free, noncontiguous cake division under piecewise linear densities. Cohler et al. [CLPP11] also address the noncontiguous version of the problem, and for piecewise constant densities they develop a polynomial-time algorithm that computes an envy-free division with optimal social welfare. In contrast to these results our focus is on contiguous cake division.

Notation and Preliminaries

This work studies the problem of dividing a cake $amongamongnagents.Throughout,wewillfocusonawell−studiedformulationofcakecuttingwhichrequiresthateachagentshouldreceiveacontiguouspieceofthecake,i.e.,thegoalistopartitionthecakeagents. Throughout, we will focus on a well-studied formulation of cake cutting which requires that each agent should receive a contiguous piece of the cake, i.e., the goal is to partition the cakeintointonpairwisedisjointintervalsandassignthemamongthepairwise disjoint intervals and assign them among then$ agents in a fair/efficient manner.

We additionally assume that the valuations are normalized such that the value of the entire cake is equal to one for every agent i∈[n]i\in[n], i.e., ∫01fi(x)dx=vi(0,1)=1\int_{0}^{1}f_{i}(x)dx=v_{i}(0,1)=1. Hence, the value-densities fif_{i}s constitute probability density functions over the cake $$. We note that this is a standard assumption in the cake-cutting framework and we conform to it for the purpose of brevity. All of our results hold true, even otherwise.

For an allocation I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\}, the endpoints of the constituent intervals will be referred to as the cut points of I\mathcal{I}, i.e., if Ii=[xi−1,xi]I_{i}=[x_{i-1},x_{i}] for 1≤i≤n1\leq i\leq n, then the cut-points are {x0=0,x1,…,xn=1}\{x_{0}=0,x_{1},\ldots,x_{n}=1\}.

We will throughout use the term allocation to specifically refer to partitions of the cake in which each agent receives a connected piece, i.e., receives exactly one interval. More generally, a cake division will be used to denote partitions of the cake D={D1,D2,…,Dn}\mathcal{D}=\{D_{1},D_{2},\ldots,D_{n}\} in which agent ii receives DiD_{i}, a finite collection of intervals. Here, the bundles DiD_{i}s are pairwise disjoint and their union covers the entire cake $$.

In this work we develop algorithmic results for the following notions of fairness and economic efficiency.

Envy-Freeness: For a cake-division instance C\mathcal{C}, an allocation I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\} is said to be envy-free iif each agent prefers its own interval over that of any other agent, vi(Ii)≥vi(Ij)v_{i}(I_{i})\geq v_{i}(I_{j}) for all agents i,j∈[n]i,j\in[n].

Pareto Optimality: Given a cake-division instance C\mathcal{C}, a division D={D1,…,Dn}\mathcal{D}=\{D_{1},\ldots,D_{n}\} is said to Pareto dominate another division C={C1,…,Cn}\mathcal{C}=\{C_{1},\ldots,C_{n}\} iff vi(Di)≥vi(Ci)v_{i}(D_{i})\geq v_{i}(C_{i}) for all agents i∈[n]i\in[n] and, there exists at least one agent k∈[n]k\in[n] such that vk(Dk)>vk(Ck)v_{k}(D_{k})>v_{k}(C_{k}). Consequently, a cake division is said to be Pareto optimal iff it is not Pareto dominated by any other division.

Recall that a cake division refers to a partition of the cake in which agent receives a finite collection of intervals. By contrast, in an allocation each agent receives a single interval. The algorithms developed in this work compute allocations. Interestingly, though, the Pareto optimality guarantees achieved by our algorithms are stronger in the sense that optimality holds across all cake divisions; specifically, under MLRP, we establish that particular allocations are Pareto optimal not only among the set of all allocations but also among all cake divisions.

Finding allocations that maximize the above-mentioned welfare notions is known to be APX-hard, in general; see, e.g., [ADH13, ABKR19]. Complementing these negative results, a key contribution of this work is to identify a broad class of cake-division instances that admit strong algorithmic results for these welfare objectives and envy-freeness. Specifically, we focus on value densities (distributions) that satisfy the monotone likelihood ratio property (MLRP). We will next define this property and note that our results hold for multiple distribution families that satisfy MLRP.

In the cake-division context, we will say that a given collection {fi}i∈[n]\{f_{i}\}_{i\in[n]} of value-density functions satisfies the monotone likelihood ratio property iff there exists an order π\mathchar58[n]→[n]\pi\mathrel{\mathop{\mathchar 58\relax}}[n]\rightarrow[n], among the fif_{i}s, such that, for all i∈[n−1]i\in[n-1], the consecutive likelihood ratios fπ(i+1) (x)fπ(i) (x)\frac{f_{\pi(i+1)}\ (x)}{f_{\pi(i)}\ (x)} are non-decreasing in x∈x\in. That is, for each i∈[n−1]i\in[n-1], the densities fπ(i)f_{\pi(i)} and fπ(i+1)f_{\pi(i+1)} bear MLRP over $$.

We will refer to this order π\pi as the MLRP order of the value densities. Lemma 9 (in Appendix A.2) shows that, given a cake-division instance in the Robertson-Webb query model, with the promise that the underlying value densities satisfy MLRP (i.e., given a promise problem), one can efficiently find the MLRP order π\pi. Hence, without loss of generality, we will throughout assume that the nn agents are indexed such that π\pi is the identity permutation, i.e., for all i∈[n−1]i\in[n-1], the likelihood ratio fi+1(x)fi(x)\frac{f_{i+1}(x)}{f_{i}(x)} is non-decreasing in x∈x\in.

It is relevant to note that, to be well defined, MLRP requires the value densities fif_{i}s to be strictly positive over the cake $.Hence,forcake−divisioninstances. Hence, for cake-division instances\langle[n],\{f_{i}\}_{i\in[n]}\ranglewithMLRP,wehavewith MLRP, we havef_{i}(x)>0forallfor alli\in[n]andandx\in$.

Instantiations of MLRP: MLRP induces a total order on linear value densities fi(x)=aix+bif_{i}(x)=a_{i}x+b_{i}; see Appendix C for details. Hence, our results imply that if, in a cake-division instance, the value densities of all the agents are linear, then an envy-free (or welfare-maximizing) allocation can be computed efficiently.

Many other distribution families are also known to bear MLRP, e.g., Gaussian distributions (with the same variance), Poisson distributions, and single-parameter exponentials. Therefore, our algorithmic results address, in particular, cake-division instances wherein all the agents have Gaussian value densities with the same variance, but different means.

These instantiations substantiate the applicability of our algorithmic results which, through MLRP, address a wide range of cake-division instances.

Lipschitz Constant of Cut and Eval Queries: We say that the cut and eval queries in a cake-division instance are λ\lambda-Lipschitz iff the following inequalities hold for each agent i∈[n]i\in[n]:

It is worth pointing out that besides MLRP (and, hence, the positivity of the value densities), all the other assumptions made in this work are standard.

The time complexities of our algorithms depend polynomially on the the bit complexity of the Lipschitz constant λ≥1\lambda\geq 1.In the case of linear value densities, fi(x)=aix+bif_{i}(x)=a_{i}x+b_{i}, the bit complexity of the Lipschitz constant λ\lambda is proportional to the bit complexity of the coefficients aia_{i}s and bib_{i}s (Proposition 3). Hence, if linear densities are explicitly given as input, then we have a polynomial (in the input size) runtime bound. As mentioned previously, in general, such a runtime dependency on log⁡λ\log\lambda is unavoidable; Appendix D provides an illustrative examples. Furthermore, it is possible—even with rational and MLRP value densities—that the exact envy-free/welfare-maximizing allocations are induced by irrational cuts (Appendix D). That is, in general, one cannot expect an efficient algorithm that outputs an exact envy-free (or welfare-maximizing) allocation. Therefore, when considering efficient algorithms for cake division, a precision loss is inevitable. However, our algorithms ensure that this loss is arbitrarily small. Specifically, given a cake-division instance C\mathcal{C} in which the value densities satisfy MLRP, we can find, in time that is polynomial in log⁡(1/η)\log(1/\eta) (along with nn and log⁡λ\log\lambda), an envy-free allocation I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\} such that, vi(Ii)≥vi(Ij)−ηv_{i}(I_{i})\geq v_{i}(I_{j})-\eta for all agents i,j∈[n]i,j\in[n]. Since the precision parameter η\eta can be driven exponentially close to zero in polynomial (in the bit complexity of η\eta) time, we will say that an envy-free allocation can be computed up to an arbitrary precision.

Similarly, in the context of maximizing welfare (social or egalitarian) welfare, given a cake-division instance wherein the value densities bear MLRP, we can find—in time that is polynomial in log⁡(1/η)\log(1/\eta)—an allocation with the (social or egalitarian) welfare η\eta (additively) close to the optimal. Hence, as in the case of envy-freeness, we assert that a welfare-maximizing allocation can be computed efficiently, up to an arbitrary precision.

Main Results

This section presents the statements of our key results.

Envy-Freeness: In Section 5 we prove that for cake-division instances, in which the value densities satisfy MLRP, the problem of finding an envy-free allocation essentially admits a polynomial-time algorithm.

Recall that, under MLRP, the cut and eval queries are necessarily λ\lambda-Lipschitz (Proposition 3).

Pareto Optimality: Weller’s theorem [Wel85] asserts that there always exists some cake division (though, not necessarily with connected pieces) which is both envy-free and Pareto optimal (among all cake divisions, with or without connected pieces). We show that, in the context of MLRP, every envy-free allocation is in fact Pareto optimal (among all cake divisions). Therefore, for cake-division instances with MLRP, the allocation computed by our algorithm is not only envy-free (fair) but also Pareto optimal, up to an arbitrary precision.

Let C\mathcal{C} be a cake-division instance wherein the value-density functions satisfy the monotone likelihood ratio property. Then, every envy-free allocation in C\mathcal{C} is also Pareto optimal (over the set of all cake divisions).

Social Welfare: In Section 7 we show that, up to an arbitrary precision, a social welfare maximizing allocation can be computed efficiency under MLRP.

Egalitarian Welfare: Section 8 addresses the problem of maximizing egalitarian welfare. Specifically, we prove that, in cake-division instances with MLRP, an allocation with egalitarian welfare arbitrarily close to the optimal can be computed efficiently.

Nash Social Welfare: In Section 9 we show that, under MLRP, the problem of maximizing Nash social welfare admits a fully polynomial-time approximation scheme (FPTAS).

Implications of the Monotone Likelihood Ratio Property

This section provides useful implications of MLRP (Lemma 1). This result will be used in subsequent sections towards the analysis of our algorithms.

Let fif_{i} and fjf_{j} be two (ordered) value-density functions that bear MLRP i.e., the likelihood ratios satisfy fj(b)fi(b)≤fj(c)fi(c)\frac{f_{j}(b)}{f_{i}(b)}\leq\frac{f_{j}(c)}{f_{i}(c)} for all 0≤b≤c≤10\leq b\leq c\leq 1. Then, fif_{i} and fjf_{j} satisfy the following two properties

The values of the intervals satisfy ∫abfj∫abfi≤∫cdfj∫cdfi\frac{\int\limits_{a}^{b}f_{j}}{\int\limits_{a}^{b}f_{i}}\leq\frac{\int\limits_{c}^{d}f_{j}}{\int\limits_{c}^{d}f_{i}} for all [a,b],[c,d]⊆[a,b],[c,d]\subseteq with b≤cb\leq c.

The (normalized) values satisfy ∫xbfi∫abfi≤∫xbfj∫abfj\frac{\int\limits_{x}^{b}f_{i}}{\int\limits_{a}^{b}f_{i}}\leq\frac{\int\limits_{x}^{b}f_{j}}{\int\limits_{a}^{b}f_{j}} for all intervals [a,b]⊆[a,b]\subseteq and all x∈[a,b]x\in[a,b].

Moreover, properties (i) and (ii) are equivalent.

Here, if the likelihood ratio fj(x)fi(x)\frac{f_{j}(x)}{f_{i}(x)} is strictly increasing, then we have a strict inequality in the corresponding implications.

The proof of this lemma is delegated to Appendix A.1. Note that, in terms of the agents’ valuations viv_{i} and vjv_{j}, property (i) in Lemma 1 can be expressed as vj(a,b)vi(a,b)≤vj(c,d)vi(c,d)\frac{v_{j}(a,b)}{v_{i}(a,b)}\leq\frac{v_{j}(c,d)}{v_{i}(c,d)} for all [a,b],[c,d]⊆[a,b],[c,d]\subseteq with b≤cb\leq c. Similarly, property (ii) corresponds to vi(x,b)vi(a,b)≤vj(x,b)vj(a,b)\frac{v_{i}(x,b)}{v_{i}(a,b)}\leq\frac{v_{j}(x,b)}{v_{j}(a,b)} for all [a,b]⊆[a,b]\subseteq and all x∈[a,b]x\in[a,b].

It is well-known that MLRP implies first-order stochastic dominance (see Lemma 8). Interestingly, property (ii) provides a strengthening: over any interval [a,b][a,b], the normalized (by vj(a,b)v_{j}(a,b)) values of agent jj stochastically dominate the normalized (by vi(a,b)v_{i}(a,b)) values of agent ii.

Envy-Freeness

In this section we develop an efficient algorithm for finding envy-free allocations in cake division instances with MLRP (Theorem 1). Towards this goal, we define a class of cake divisions, referred to as ripple divisions (Definition 1), and prove that, under MLRP, every ripple division induces an envy-free allocation (Theorem 6). Existential and computational guarantees for ripple divisions are established in Section 5.1 (Lemma 3) and Section 5.2 (Lemma 4), respectively. Section 5.4 builds upon these results to prove our main result (Theorem 1) for envy-freeness.

We establish the universal existence of ripple divisions through the intermediate value theorem, i.e., a one-dimensional fixed-point argument (Lemma 3). Consequentially, for cake-division instances with MLRP, we develop an alternate proof of existence of envy-free allocations. Since one can use binary search to find fixed points in the one-dimensional setting, this proof in fact leads to an algorithm for finding ripple divisions and, hence, envy-free divisions.

Given a cake-division instance C=⟨[n],{fi}i⟩\mathcal{C}=\langle[n],\{f_{i}\}_{i}\rangle, a collection of points x0∗=0≤x1∗≤⋯≤xn−1∗≤xn∗=1x^{*}_{0}=0\leq x^{*}_{1}\leq\dots\leq x^{*}_{n-1}\leq x^{*}_{n}=1 is said to form a ripple division of the cake iff

For establishing existence of ripple divisions, we first consider a relaxation of Definition 1 wherein do not enforce the last cut point (i.e., xn∗x^{*}_{n}) to be equal to one. Under this relaxation, the intervals {[xi−1,xi]}i=1n\{[x_{i-1},x_{i}]\}_{i=1}^{n} do not cover the entire interval $(instead,theyspan(instead, they span[0,x_{n}]$) and, hence, lead to a partial allocation of the cake.

Also, by convention, the agents are indexed following the MLRP order: for each i∈[n−1]i\in[n-1], the likelihood ratio fi+1/fif_{i+1}/f_{i} is nondecreasing. Hence, assigning interval [xi−1,xi][x_{i-1},x_{i}] to agent i∈[n]i\in[n] provides an allocation wherein the intervals are assigned (left to right on the cake) in accordance with the MLRP order. We will show that the (partial) allocation obtained by assigning interval [xi−1,xi][x_{i-1},x_{i}] to agent i∈[n]i\in[n] is envy-free (Theorem 6).

Given a cake-division instance C=⟨[n],{fi}i⟩\mathcal{C}=\langle[n],\{f_{i}\}_{i}\rangle, a collection of points x0=0≤x1≤⋯≤xn−1≤xn≤1x_{0}=0\leq x_{1}\leq\dots\leq x_{n-1}\leq x_{n}\leq 1 is said to form a δ\delta-ripple division of the cake iff xn≥1−δx_{n}\geq 1-\delta and

Both Definitions 1 and 2 require that, for all i∈[n−1]i\in[n-1], agent ii’s value for the iith interval ([xi−1∗,xi∗][x^{*}_{i-1},x^{*}_{i}] and [xi−1,xi][x_{i-1},x_{i}], respectively) is positive. Also, note that a -ripple division is an exact ripple division.

Next we use Lemma 2 and a limit argument (δ→0\delta\to 0) to establish universal existence of ripple divisions.

Let C\mathcal{C} be a cake-division instance in which the cut and eval queries are λ\lambda-Lipschitz. Then, C\mathcal{C} necessarily admits a ripple division.

Proof Lemma 2 asserts that, for any δ∈(0,1)\delta\in(0,1), there exists a collection of points x0δ=0<x1δ<…<xn−1δ<xnδ=1−δx^{\delta}_{0}=0<x^{\delta}_{1}<\ldots<x^{\delta}_{n-1}<x^{\delta}_{n}=1-\delta that form a δ\delta-ripple division. Note that here xnδ=1−δx^{\delta}_{n}=1-\delta. We consider the sequence of these δ\delta-ripple divisions, S≔⟨(x0δ,x1δ,…,xnδ)⟩δ∈(0,1)\mathcal{S}\coloneqq\left\langle(x^{\delta}_{0},x^{\delta}_{1},\ldots,x^{\delta}_{n})\right\rangle_{\delta\in(0,1)}.

Since S\mathcal{S} is a nonempty and bounded sequence in n^{n}, the Bolzano-Weierstrass theorem ensures that S\mathcal{S} contains a convergent subsequence, say ⟨(x0δj,x1δj,…,xnδj)⟩δj>0\left\langle({x}^{\delta_{j}}_{0},{x}^{\delta_{j}}_{1},\dots,{x}^{\delta_{j}}_{n})\right\rangle_{\delta_{j}>0}. Write (x0∗,x1∗,…,xn∗)∈n(x^{*}_{0},x^{*}_{1},\dots,x^{*}_{n})\in^{n} to denote the limit of this subsequence as δj\delta_{j} tends to zero. We will show that the points x0∗,x1∗,…,xn∗x^{*}_{0},x^{*}_{1},\ldots,x^{*}_{n} form a ripple division in C\mathcal{C} (see Definition 1), i.e., establish that xn∗=1x^{*}_{n}=1 and vi(xi−1∗,xi∗)=vi(xi∗,xi+1∗)>0v_{i}(x^{*}_{i-1},x^{*}_{i})=v_{i}(x^{*}_{i},x^{*}_{i+1})>0 for all agents i∈[n−1]i\in[n-1].

First, note that xn∗=1x^{*}_{n}=1, since the sequence ⟨1−δj⟩→1\langle 1-\delta_{j}\rangle\to 1 as δj→0\delta_{j}\to 0. Also, x0∗=0x^{*}_{0}=0, since the constant sequence ⟨0⟩\langle 0\rangle tends to .

We will next prove that vi(xi−1∗,xi∗)=vi(xi∗,xi+1∗)v_{i}(x^{*}_{i-1},x^{*}_{i})=v_{i}(x^{*}_{i},x^{*}_{i+1}) for all agents i∈[n−1]i\in[n-1]. Given that the collection (x0δ,x1δ,…,xnδj)\left({x}^{\delta}_{0},{x}^{\delta}_{1},\dots,{x}^{\delta_{j}}_{n}\right) forms a δj\delta_{j}-ripple division, we have vi(xi−1δj,xiδj)=vi(xiδj,xi+1δj)v_{i}({x}^{\delta_{j}}_{i-1},{x}^{\delta_{j}}_{i})=v_{i}({x}^{\delta_{j}}_{i},{x}^{\delta_{j}}_{i+1}), for all i∈[n−1]i\in[n-1].

2 Computation of Ripple Divisions

Let C=⟨[n],{fi}i⟩\mathcal{C}=\langle[n],\{f_{i}\}_{i}\rangle be a cake-division instance in which the cut and eval queries are λ\lambda-Lipschitz. Then, for any δ∈(0,1)\delta\in(0,1) and in the Robertson-Webb query model, a δ\delta-ripple division of C\mathcal{C} can be computed in O(poly(n,log⁡λ,log⁡1δ))\mathcal{O}\left({\rm poly}(n,\log\lambda,\log\frac{1}{\delta})\right) time.

3 From Ripple Divisions to Envy-Free Allocations

The next theorem establishes the crucial connection between ripple divisions and envy-freeness. In particular, setting δ=0\delta=0 in this theorem, we obtain that, under MLRP, every (exact) ripple division induces an envy-free allocation.

Let C=⟨[n],{fi}i⟩\mathcal{C}=\langle[n],\{f_{i}\}_{i}\rangle be a cake-division instance in which the value densities satisfy the monotone likelihood ratio property and let parameter δ≥0\delta\geq 0. Then, every δ\delta-ripple division, x0=0≤x1≤⋯≤xn−1≤xn≤1x_{0}=0\leq x_{1}\leq\dots\leq x_{n-1}\leq x_{n}\leq 1, in C\mathcal{C} induces an envy-free partial allocation {Ii=[xi−1,xi]}i=1n\{I_{i}=[x_{i-1},x_{i}]\}_{i=1}^{n}.

This theorem asserts that here the partial allocation I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\} (with Ii=[xi−1,xi]I_{i}=[x_{i-1},x_{i}]) satisfies vi(Ii)≥vi(Ij)v_{i}(I_{i})\geq v_{i}(I_{j}), for all i,j∈[n]i,j\in[n], and at most a δ\delta-length piece of the cake (specifically, [xn,1][x_{n},1]) remains unallocated in I\mathcal{I}.

Proof We will show that if the points x0=0≤x1≤x2≤…≤xn≤1x_{0}=0\leq x_{1}\leq x_{2}\leq\ldots\leq x_{n}\leq 1 form a δ\delta-ripple division, then the partial allocation {Ii=[xi−1,xi]}i=1n\{I_{i}=[x_{i-1},x_{i}]\}_{i=1}^{n} is envy-free. Here, the definition of a δ\delta-ripple division (Definition 2) ensures that, for each agent i∈[n−1]i\in[n-1], the values the two consecutive intervals IiI_{i} and Ii+1I_{i+1} are equal and positive

Recall that the agents are indexed following the MLRP order, i.e., for each i∈[n−1]i\in[n-1], the likelihood ratio fi+1/fif_{i+1}/f_{i} is nondecreasing. We fix an agent i∈[n]i\in[n] and establish envy-freeness with respect to ii by considering two complementary cases (i) for agents to the left of ii, we prove that vi(I1)≤vi(I2)≤…≤vi(Ii)v_{i}(I_{1})\leq v_{i}(I_{2})\leq\ldots\leq v_{i}(I_{i}) and (ii) for agents to the right of ii, we prove that vi(Ii)≥vi(Ii+1)≥…≥vi(In)v_{i}(I_{i})\geq v_{i}(I_{i+1})\geq\ldots\geq v_{i}(I_{n}).

Case (i): Consider any agent k∈[n]k\in[n] such that k<ik<i. Given that fkf_{k} and fif_{i} bear MLRP (i.e., fi/fkf_{i}/f_{k} is non-decreasing), property (i) of Lemma 1, with a=xk−1a=x_{k-1}, b=c=xkb=c=x_{k}, and d=xk+1d=x_{k+1}, gives us vi(xk−1,xk)vk(xk−1,xk)≤vi(xk,xk+1)vk(xk,xk+1)\frac{v_{i}(x_{k-1},x_{k})}{v_{k}(x_{k-1},x_{k})}\leq\frac{v_{i}(x_{k},x_{k+1})}{v_{k}(x_{k},x_{k+1})}. That is, for the intervals Ik=[xk−1,xk]I_{k}=[x_{k-1},x_{k}] and Ik+1=[xk,xk+1]I_{k+1}=[x_{k},x_{k+1}] we have

Instantiating equation (2) for agent kk, we can simplify inequality (3) to vi(Ik)≤vi(Ik+1)v_{i}(I_{k})\leq v_{i}(I_{k+1}). Combining this inequality across all k<ik<i, we obtain the desired chain of inequalities for agent ii, i.e., vi(I1)≤vi(I2)≤…≤vi(Ii)v_{i}(I_{1})\leq v_{i}(I_{2})\leq\ldots\leq v_{i}(I_{i}).

Case (ii): Consider any agent j∈[n]j\in[n] such that j>ij>i Given that fif_{i} and fjf_{j} bear MLRP (i.e., fj/fif_{j}/f_{i} is non-decreasing), property (i) of Lemma 1, with a=xj−1a=x_{j-1}, b=c=xjb=c=x_{j}, and d=xj+1d=x_{j+1}, gives us vj(xj−1,xj)vi(xj−1,xj)≤vj(xj,xj+1)vi(xj,xj+1)\frac{v_{j}(x_{j-1},x_{j})}{v_{i}(x_{j-1},x_{j})}\leq\frac{v_{j}(x_{j},x_{j+1})}{v_{i}(x_{j},x_{j+1})}. That is, for the intervals Ij=[xj−1,xj]I_{j}=[x_{j-1},x_{j}] and Ij+1=[xj,xj+1]I_{j+1}=[x_{j},x_{j+1}] we have

Instantiating equation (2) for agent jj, we can simplify inequality (4) to vi(Ij)≥vi(Ij+1)v_{i}(I_{j})\geq v_{i}(I_{j+1}). Combining this inequality across all j>ij>i, we obtain the desired chain of inequalities for agent ii, i.e., vi(Ii)≥vi(Ii+1)≥…≥vi(In)v_{i}(I_{i})\geq v_{i}(I_{i+1})\geq\ldots\geq v_{i}(I_{n}).

The above two cases establish that agent i∈[n]i\in[n] does not envy any other other agent, i.e., I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\} is an envy-free partial allocation. Indeed, if δ=0\delta=0, then I\mathcal{I} covers the entire cake, i.e., we obtain an envy-free allocation. ⊓\sqcap⊔\sqcup

Notably, Lemma 3 and Theorem 6 (with δ=0\delta=0) provide a stand-alone proof of existence of envy-free allocations in cake-division instances with MLRP. The next section establishes an algorithmic counterpart of this existential result; specifically, we show that using Lemma 4 one can directly obtain an efficient algorithm for finding envy-free allocations, under MLRP.

4 Proof of Theorem 1

This section restates and proves our main result (Theorem 1) for envy-freeness.

Proof Given a cake-division instance C\mathcal{C}, with MLRP, and precision parameter η>0\eta>0, we invoke Lemma 4 to find an (ηλ)\left(\frac{\eta}{\lambda}\right)-ripple division in O(poly(n,log⁡λ,log⁡1η))\mathcal{O}\left({\rm poly}(n,\log\lambda,\log\frac{1}{\eta})\right) time.

Write x0=0≤x1≤⋯≤xn−1≤xn≤1x_{0}=0\leq x_{1}\leq\dots\leq x_{n-1}\leq x_{n}\leq 1 to denote the computed (ηλ)\left(\frac{\eta}{\lambda}\right)-ripple division and let I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\} be the corresponding partial allocation; here Ii=[xi−1,xi]I_{i}=[x_{i-1},x_{i}]. Theorem 6 ensures that I\mathcal{I} is envy-free.

We will show that coalescing the unassigned (in I\mathcal{I}) piece [xn,1][x_{n},1] to agent nn provides a complete allocation that satisfies envy-freeness, up to η\eta precision. Write I∗≔{I1∗,I2∗,…,In−1∗,In∗}\mathcal{I}^{*}\coloneqq\{I^{*}_{1},I^{*}_{2},\ldots,I^{*}_{n-1},I^{*}_{n}\} to denote this allocation in which the nnth agent receives the interval In∗≔[xn−1,1]I^{*}_{n}\coloneqq[x_{n-1},1] (equivalently, In∗=In∪[xn,1]I^{*}_{n}=I_{n}\cup[x_{n},1]) and Ii∗=IiI^{*}_{i}=I_{i} for the remaining agents i∈[n−1]i\in[n-1].

Note that against all agents j∈[n−1]j\in[n-1], envy-freeness of I∗\mathcal{I}^{*} directly follows from the fact that the partial allocation I\mathcal{I} is envy-free: vi(Ii∗)≥vi(Ii)≥vi(Ij)=vi(Ij∗)v_{i}(I^{*}_{i})\geq v_{i}(I_{i})\geq v_{i}(I_{j})=v_{i}(I^{*}_{j}), for all i∈[n]i\in[n] and all j∈[n−1]j\in[n-1].

Finally, we address envy against agent nn. Recall that xix_{i}s form a (ηλ)\left(\frac{\eta}{\lambda}\right)-ripple division, hence xn≥1−ηλx_{n}\geq 1-\frac{\eta}{\lambda}. In addition, the fact that eval queries are λ\lambda-Lipschitz gives us vi([xn,1])≤ηv_{i}([x_{n},1])\leq\eta for all agents i∈[n]i\in[n]. Hence, for all i∈[n]i\in[n] we have

Therefore, I∗\mathcal{I}^{*} satisfies envy-freeness, up to η\eta precision: vi(Ii∗)≥vi(Ij∗)−ηv_{i}(I^{*}_{i})\geq v_{i}(I^{*}_{j})-\eta for all i,j∈[n]i,j\in[n].

The time complexity obtained via Lemma 4 implies that the parameter η\eta can be driven exponentially close to zero, in time that is polynomial in log⁡1η\log\frac{1}{\eta} (i.e., in the bit complexity of η\eta). Hence, we can find an envy-free allocation, up to arbitrary precision, in O(poly(n,log⁡λ))\mathcal{O}\left({\rm poly}(n,\log\lambda)\right) time. Theorem 1 now stands proved. ⊓\sqcap⊔\sqcup

Pareto Optimality

This section shows that, with MLRP in hand, one does not loose out on Pareto optimality by imposing the contiguity requirement. That is, under MLRP, there always exist allocations (i.e., cake divisions with connected pieces) that are Pareto optimal among all cake divisions, with or without connected pieces. Moreover, such allocations conform to the MLRP order.

This structural result implies that for maximizing welfare we can restrict attention to allocations wherein the intervals are assigned (left to right on the cake) in accordance with the MLRP order. Intuitively, this leads us to a welfare-maximizing algorithm—specifically, a dynamic program—that recursively finds optimal allocations for intervals placed at the left end of the cake (i.e., for intervals of the form [0,x][0,x]); see Section 7 and Appendix 9 for details.

Subsequently, Section 6.1 (Theorem 2) establishes a strong connection between fairness and (Pareto) efficiency in the MLRP context: if the value densities bear MLRP, then every envy-free allocation is necessarily Pareto optimal.

Let C=⟨[n],(fi)i⟩\mathcal{C}=\langle[n],(f_{i})_{i}\rangle be a cake-division instance in which the value densities satisfy the monotone likelihood ratio property. Then, for every cake division D={D1,…,Dn}\mathcal{D}=\{D_{1},\ldots,D_{n}\} in C\mathcal{C} there exists an allocation J={J1,…,Jn}\mathcal{J}=\{J_{1},\ldots,J_{n}\} such that vi(Ji)≥vi(Di)v_{i}(J_{i})\geq v_{i}(D_{i}), for 1≤i≤n1\leq i\leq n.

Furthermore, for every Pareto optimal allocation I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\} in C\mathcal{C}, there exists an allocation I′={I1′,…,In′}\mathcal{I^{\prime}}=\{I^{\prime}_{1},\ldots,I^{\prime}_{n}\} with vi(Ii′)=vi(Ii)v_{i}(I^{\prime}_{i})=v_{i}(I_{i}) that conforms to the MLRP order, i.e., if  fi+1/fi\ f_{i+1}/f_{i} is nondecreasing in $,thentheintervalassignedtoagent, then the interval assigned to agenti(i.e.,(i.e.,I^{\prime}_{i})appearstotheleftoftheintervalassignedtotheagent) appears to the left of the interval assigned to the agenti+1(i.e.,(i.e.,I^{\prime}_{i+1}$).

Proof Consider a cake division D={D1,…,Dn}\mathcal{D}=\{D_{1},\ldots,D_{n}\} wherein two consecutive intervals are assigned violating the MLRP order: say, interval [p,q][p,q] is assigned to agent jj (i.e., this interval is contained in the bundle DjD_{j}), the adjacent interval [q,r][q,r] is assigned to agent ii, and agent ii appears before jj in the MLRP order (fj/fif_{j}/f_{i} is non-decreasing over $$).

We will show that in such a case there always exists a point q′∈[p,r]q^{\prime}\in[p,r] such that vi(p,q′)≥vi(q,r)v_{i}(p,q^{\prime})\geq v_{i}(q,r) and vj(q′,r)≥vj(p,q)v_{j}(q^{\prime},r)\geq v_{j}(p,q). That is, one can swap the allocation order between ii and jj (in the interval [p,q]∪[q,r][p,q]\cup[q,r]) without decreasing the agents’ values. Moreover, we note that, if fj/fif_{j}/f_{i} is strictly increasing in the interval [p,q]∪[q,r][p,q]\cup[q,r], then this update leads to a strict increase in agent ii’s or agent jj’s value.

Hence, starting with any cake division D={D1,…,Dn}\mathcal{D}=\{D_{1},\ldots,D_{n}\}, we can repeatedly apply the above-mentioned resolution towards the MLRP order and obtain an allocation J={J1,…,Jn}\mathcal{J}=\{J_{1},\ldots,J_{n}\} with the desired property, vi(Ji)≥vi(Di)v_{i}(J_{i})\geq v_{i}(D_{i}) for all i∈[n]i\in[n].

Note that this resolution process also establishes the second part of the theorem, i.e., for every Pareto optimal allocation I\mathcal{I}, there exists an allocation I′\mathcal{I^{\prime}} with vi(Ii′)=vi(Ii)v_{i}(I^{\prime}_{i})=v_{i}(I_{i}) that conforms to the MLRP order.

The remainder of the proof addresses the desired point q′∈[p,r]q^{\prime}\in[p,r]. In particular, we will identify q′q^{\prime} such that assigning interval [p,q′][p,q^{\prime}] to agent ii (instead of [q,r][q,r]) and assigning [q′,r][q^{\prime},r] to agent jj (instead of [p,q][p,q]) leads to an increment in values.

Write βi≔vi(q,r)vi(p,r)\beta_{i}\coloneqq\frac{v_{i}(q,r)}{v_{i}(p,r)} to denote the normalized value of agent ii under the initial assignment. Define q′∈[p,r]q^{\prime}\in[p,r] to be the point that satisfies

Since the value densities satisfy MLRP, property (ii) of Lemma 1, applied to the interval [p,r][p,r] and q′∈[p,r]q^{\prime}\in[p,r], gives us vi(q′,r)vi(p,r)≤vj(q′,r)vj(p,r)\frac{v_{i}(q^{\prime},r)}{v_{i}(p,r)}\leq\frac{v_{j}(q^{\prime},r)}{v_{j}(p,r)}. Simplifying further we obtain 1−vj(q′,r)vj(p,r)≤1−vi(q′,r)vi(p,r)1-\frac{v_{j}(q^{\prime},r)}{v_{j}(p,r)}\leq 1-\frac{v_{i}(q^{\prime},r)}{v_{i}(p,r)}, i.e.,

Therefore, we obtain a value bound for agent ii

That is, agent ii’s value is preserved through the reassignment, vi(p,q′)≥vi(q,r)v_{i}(p,q^{\prime})\geq v_{i}(q,r).

For agent jj, via property (ii) of Lemma 1, with interval [p,r][p,r] and q∈[p,r]q\in[p,r], we have

This inequality reduces to 1−vj(p,q′)vj(p,r)≥1−vj(q,r)vj(p,r)1-\frac{v_{j}(p,q^{\prime})}{v_{j}(p,r)}\geq 1-\frac{v_{j}(q,r)}{v_{j}(p,r)}. Simplifying we obtain vj(p,r)−vj(p,q′)vj(p,r)≥vj(p,r)−vj(q,r)vj(p,r)\frac{v_{j}(p,r)-v_{j}(p,q^{\prime})}{v_{j}(p,r)}\geq\frac{v_{j}(p,r)-v_{j}(q,r)}{v_{j}(p,r)}. Therefore, we have the desired inequality vj(q′,r)≥vj(p,q)v_{j}(q^{\prime},r)\geq v_{j}(p,q) and the stated claims follow. ⊓\sqcap⊔\sqcup

Remark. For cake-division instances with MLRP, we can prove that any allocation K\mathcal{K} that conforms to the MLRP order is Pareto optimal (over the set of all cake divisions). Write 0=k0<k1<⋯<kn=10=k_{0}<k_{1}<\dots<k_{n}=1 to denote the cut-points of K\mathcal{K}. For contradiction, we assume that K\mathcal{K} is not Pareto optimal. That is, there exists a cake-division L\mathcal{L} that dominates K\mathcal{K} and is Pareto optimal. By Lemma 5, we know there exists another Pareto optimal allocation M\mathcal{M} with vi(Mi)=vi(Li)v_{i}(M_{i})=v_{i}(L_{i}) for all i∈[n]i\in[n] that conforms to the MLRP order. Write 0=m0<m1<⋯<mn=10=m_{0}<m_{1}<\dots<m_{n}=1 to denote the cut-points of M\mathcal{M}. Since, M\mathcal{M} Pareto dominates K\mathcal{K}, we will have ki≤mik_{i}\leq m_{i} for all i∈[n]i\in[n] with at least one strict inequality. This contradicts the fact that kn=mn=1k_{n}=m_{n}=1. Therefore, K\mathcal{K} is Pareto optimal. ⊓\sqcap⊔\sqcup

Weller’s theorem [Wel85] is a notable result in the cake-cutting literature and it asserts that there always exists some cake division—though, not necessarily with connected pieces—which is both envy-free (fair) and Pareto optimal. While this theorem holds in general, it does not guarantee that envy-freeness and Pareto optimality can be achieved together with allocations. Indeed, there are cake-division instances wherein none of the of envy-free allocations are Pareto optimal.

We show that, by contrast under MLRP, every envy-free allocation is Pareto optimal, among all cake divisions (Theorem 2). Therefore, given a cake-division instance with MLRP, the allocation computed by our algorithm (Algorithm 1) is not only envy-free but also Pareto optimal, up to an arbitrary precision.

Proof Write I={I1,I2,…,In}\mathcal{I}=\{I_{1},I_{2},\dots,I_{n}\} to denote an envy-free allocation in C\mathcal{C}; here interval IiI_{i} is assigned to agent i∈[n]i\in[n]. We assume, towards a contradiction, that there exists a cake division D={D1,…,Dn}\mathcal{D}=\{D_{1},\ldots,D_{n}\} that Pareto dominates I\mathcal{I}. Lemma 5 implies that in such a case there exists an allocation J={J1,…,Jn}\mathcal{J}=\{J_{1},\ldots,J_{n}\} which also Pareto dominates I\mathcal{I}. That is, we have vi(Ji)≥vi(Di)≥vi(Ii)v_{i}(J_{i})\geq v_{i}(D_{i})\geq v_{i}(I_{i}), for all agents i∈[n]i\in[n], and there exists some agent k∈[n]k\in[n] such that vk(Jk)≥vk(Dk)>vk(Ik)v_{k}(J_{k})\geq v_{k}(D_{k})>v_{k}(I_{k}).

Recall that for an allocation the endpoints of all the constituent intervals are referred to as its cut points. We break our analysis into the following two cases depending on whether I\mathcal{I} and J\mathcal{J} have the same set of cut points.

Case 1: The cut points of the allocations I\mathcal{I} and J\mathcal{J} are identical. In this case, there must exist a permutation σ\mathchar58[n]↦[n]\sigma\mathrel{\mathop{\mathchar 58\relax}}[n]\mapsto[n] such that Iσ(i)=JiI_{\sigma(i)}=J_{i} for all i∈[n]i\in[n]. Since I\mathcal{I} is envy-free, we have vi(Ii)≥vi(Iσ(i))=vi(Ji)v_{i}(I_{i})\geq v_{i}(I_{\sigma(i)})=v_{i}(J_{i}) for all agents i∈[n]i\in[n]. However, this contradicts the fact that J\mathcal{J} Pareto dominates the allocation I\mathcal{I}.

Case 2: The cut points of I\mathcal{I} and J\mathcal{J} are not identical. Since both the allocations form a partition of the same cake $,theremustexistsome, there must exist somes,t\in[n]suchthattheintervalsuch that the intervalJ_{s}isastrictsubsetoftheintervalis a strict subset of the intervalI_{t},i.e.,, i.e.,J_{s}\subset I_{t}.Envy−freenessof. Envy-freeness of\mathcal{I}givesusgives usv_{s}(I_{s})\geq v_{s}(I_{t})>v_{s}(J_{s}).Thelaststrictinequalityfollowsfromthefactthatthevaluedensity. The last strict inequality follows from the fact that the value densityf_{s}ofagentof agentshasfullsupportoverhas full support overandandJ_{s}isastrictsubsetofis a strict subset ofI_{t}.Thisbound. This boundv_{s}(I_{s})>v_{s}(J_{s})contradictsthefactthatcontradicts the fact that\mathcal{J}ParetodominatesPareto dominates\mathcal{I}andcompletestheproof.and completes the proof.\sqcap$⊔\sqcup

Social Welfare

This section develops an algorithm for social welfare maximization. Recall that, under MLRP, Pareto optimal allocations conform to the MLRP order (Lemma 5). Hence, for maximizing social welfare we can restrict attention to allocations wherein the intervals are assigned (left-to-right on the cake) in accordance with the MLRP order. This observation, in and of itself, leads to a fully-polynomial time approximation scheme for the maximizing social welfare: we can partition the cake into poly(n,1/ε){\rm poly}(n,1/\varepsilon) contiguous intervals, each of value at most ε\varepsilon, and then solve the problem using a dynamic program. We show that instead of considering a general partition we can identify a set PP—of O(n2)\mathcal{O}(n^{2}) points—such that the cut points of an optimal allocation are contained in PP. This will enable us to execute a dynamic program focusing only on the points in PP and establish Theorem 3. In sharp contrast to the FPTAS described above, our dynamic program finds an allocation with social welfare η>0\eta>0 close to the optimal in time that is dependent on log⁡(1/η)\log(1/\eta).

The next lemma provides a useful property about the switching points in the MLRP context, that is reminiscent of the ‘single-crossing’ type condition used in social choice; see [Ath01, GS96, EFS14].

Let fif_{i} and fjf_{j} be two value-density functions that satisfy MLRP, i.e., fj/fi{f_{j}}/{f_{i}} is non-decreasing over $$. Then,

For all y∈[0,pi,j)y\in[0,p_{i,j}), the likelihood ratio satisfies fj(y)/fi(y)<1{f_{j}(y)}/{f_{i}(y)}<1.

For all z∈(pij,1]z\in(p_{ij},1], the likelihood ratio satisfies fj(z)/fi(z)≥1{f_{j}(z)}/{f_{i}(z)}\geq 1.

Here, pijp_{ij} is the switching point between fif_{i} and fjf_{j}.

Proof As observed previously, pij∈p_{ij}\in exists and is unique. We begin by proving part (a) of the stated claim. Consider any point y∈[0,pi,j)y\in[0,p_{i,j}) and assume, towards a contradiction, that the likelihood ratio satisfies fj(y)fi(y)≥1\frac{f_{j}(y)}{f_{i}(y)}\geq 1. This implies that the point yy belongs to the set Lij={x∈\mathchar58fj(x)≥fi(x)}L_{ij}=\{x\in\mathrel{\mathop{\mathchar 58\relax}}{f_{j}(x)}\geq{f_{i}(x)}\}. Since y<pijy<p_{ij}, we get a contradiction to the fact that pijp_{ij} is the infimum of the set LijL_{ij}.

For proving part (b), consider any point z∈(pij,1]z\in(p_{ij},1]. Assume, towards a contradiction, that at zz we have fj(z)fi(z)<1\frac{f_{j}(z)}{f_{i}(z)}<1. Since the likelihood ratio fj(x)fi(x)\frac{f_{j}(x)}{f_{i}(x)} is non-decreasing over $(bydefinitionofMLRP),wehave(by definition of MLRP), we have\frac{f_{j}(t)}{f_{i}(t)}<1forallfor all0\leq t\leq z.Thatis,theredoesnotexistapoint. That is, there does not exist a pointt\in[0,z]withthepropertythatwith the property thatf_{j}(t)\geq f_{i}(t).Hence,. Hence,zconstitutesalowerboundforthesetconstitutes a lower bound for the setL_{ij}=\{x\in\mathrel{\mathop{\mathchar 58\relax}}{f_{j}(x)}\geq{f_{i}(x)}\}.Since. Sincep_{ij},wegetacontradictiontothefactthat, we get a contradiction to the fact thatp_{ij}istheinfimum(greatestlowerbound)ofthesetis the infimum (greatest lower bound) of the setL_{ij}.Thiscompletestheproof.. This completes the proof.\sqcap$⊔\sqcup

The following corollary asserts that, up to an arbitrary precision, each switching point can be determined efficiently.

Let fif_{i} and fjf_{j} be two value densities that bear MLRP and let parameter γ∈(0,1)\gamma\in(0,1). Then, in the Robertson-Webb model, we can find an interval of length γ\gamma that contains the switching point pijp_{ij} in O(log⁡(1/γ))\mathcal{O}\left(\log\left(1/\gamma\right)\right) time.

Proof Consider intervals of the form Bk=[(k−1)γ2, kγ2]B_{k}=\left[(k-1)\frac{\gamma}{2},\ k\frac{\gamma}{2}\right] for k∈{1,2,…,2γ}k\in\{1,2,\ldots,\frac{2}{\gamma}\}, i.e., for analysis, we discretize the cake $evenlyintointervalseachoflengthevenly into intervals each of length\gamma/2.Write. Writek^{*}todenotetheindexwiththepropertythatto denote the index with the property thatp_{ij}\in B_{k^{*}}=\left[(k^{*}-1)\frac{\gamma}{2},\ k^{*}\frac{\gamma}{2}\right].Part(a)ofLemma6impliesthat. Part (a) of Lemma 6 implies thatv_{j}(B_{k})forallfor allk.Similarly,part(b)ofLemma6givesus. Similarly, part (b) of Lemma 6 gives usv_{j}(B_{k})\geq v_{i}(B_{k})forallfor allk>k^{*}$.

Therefore, applying binary search, we can, in O(log⁡(1/γ))\mathcal{O}\left(\log\left(1/\gamma\right)\right) iterations, find the smallest index k′∈{1,2,…,2γ}k^{\prime}\in\{1,2,\dots,\frac{2}{\gamma}\} such that vj(Bk′)<vi(Bk′)v_{j}(B_{k^{\prime}})<v_{i}(B_{k^{\prime}}) and vj(Bk′+1)≥vi(Bk′+1)v_{j}(B_{k^{\prime}+1})\geq v_{i}(B_{k^{\prime}+1}). Note that the computed index k′k^{\prime} satisfies k′∈{k∗−1,k∗}k^{\prime}\in\{k^{*}-1,k^{*}\}: if, for contradiction, we have k′≤k∗−2k^{\prime}\leq k^{*}-2, then it must be the case that vj(Bk′+1)<vi(Bk′+1)v_{j}(B_{k^{\prime}+1})<v_{i}(B_{k^{\prime}+1}). This inequality contradicts the selection criterion of k′k^{\prime}. Also, the inequality k′≥k∗+1k^{\prime}\geq k^{*}+1 would lead to the contradiction vj(Bk′)≥vi(Bk′)v_{j}(B_{k^{\prime}})\geq v_{i}(B_{k^{\prime}}).

The value comparisons required to execute the binary search can be performed using eval queries. Hence, in O(log⁡(1/γ))\mathcal{O}\left(\log(1/\gamma)\right) iterations we can find an interval Bk′∪Bk′+1B_{k^{\prime}}\cup B_{k^{\prime}+1} of length γ\gamma that contains pijp_{ij}. ⊓\sqcap⊔\sqcup

This corollary implies that we can efficiently compute the set of switching points P={pi,j∈\mathchar581≤i<j≤n}P=\{p_{i,j}\in\mathrel{\mathop{\mathchar 58\relax}}1\leq i<j\leq n\}, up to an arbitrary precision. Next, we will establish the usefulness of PP.

Let C\mathcal{C} be a cake-division instance in which the value densities satisfy the monotone likelihood ratio property. Then, in C\mathcal{C}, there exists a social welfare maximizing allocation all of whose cut points belong to the set of switching points PP.

Proof Among all allocations that maximize social welfare, consider the ones that conform to the MLRP order; Lemma 5 ensures that this collection is nonempty. Furthermore, among these optimal allocations select one S={S1,S2,…,Sn}\mathcal{S}=\{S_{1},S_{2},\ldots,S_{n}\} that minimizes ∣{s0=0,s1,s2,…,sn=1}∖P∣\left|\{s_{0}=0,s_{1},s_{2},\ldots,s_{n}=1\}\setminus P\right|; here sis_{i}s denote the cut points of S\mathcal{S}. That is, S\mathcal{S} is a social welfare maximizing allocation that uses as many points from PP as possible. We will show that {si}i∖P=∅\{s_{i}\}_{i}\setminus P=\emptyset and, hence, the claim follows.

Towards a contradiction, assume that there exists a cut point sts_{t} of the allocation S\mathcal{S} that does not belong to PP. Let Si=[s,st]S_{i}=[s,s_{t}] and Sj=[st,s′]S_{j}=[s_{t},s^{\prime}] be the two nonempty intervals in S\mathcal{S} that are separated by sts_{t}. Interval SiS_{i} is to the (immediate) left of SjS_{j} and, since the allocations conform to the MLRP order, we that i<ji<j.

Given that st∉Ps_{t}\notin P, we know that st≠pijs_{t}\neq p_{ij}; here pijp_{ij} is the switching point between fif_{i} and fjf_{j}. We will show that in this case we can always move sts_{t} towards pijp_{ij} and obtain another social welfare maximizing allocation that uses more cut points from PP than S\mathcal{S}. This contradicts the choice of S\mathcal{S} and establishes the stated claim. Towards this goal, consider two complementary cases

Case (i): st<pijs_{t}<p_{ij}. In this case we can move sts_{t} to the right without decreasing the social welfare. In particular, if pij∈Si∪Sj=[s,s′]p_{ij}\in S_{i}\cup S_{j}=[s,s^{\prime}], then, instead of SiS_{i} and SjS_{j}, we can assign intervals [s,pij][s,p_{ij}] and [pij,s′][p_{ij},s^{\prime}] to agents ii and jj, respectively. Since fj(x)<fi(x){f_{j}(x)}<{f_{i}(x)} for all x∈[st,pij]x\in[s_{t},p_{ij}] (Lemma 6, part (a)), such an update increases the social welfare. This contradicts the optimality (with respect to social welfare) of S\mathcal{S}. A similar argument holds if pij>s′p_{ij}>s^{\prime}. Here, we can assign [s,s′][s,s^{\prime}] entirely to agent ii (and an empty set to agent jj). For all x≤s′<pijx\leq s^{\prime}<p_{ij}, we have fj(x)<fi(x){f_{j}(x)}<{f_{i}(x)} (Lemma 6, part (a)). Therefore, the reassignment increase the social welfare and leads to a contradiction.

Case (ii): st>pijs_{t}>p_{ij}. In this case we can move sts_{t} to the left (towards pijp_{ij}). If we have pij∈Si∪Sjp_{ij}\in S_{i}\cup S_{j}, then assigning intervals [s,pij][s,p_{ij}] and [pij,s′][p_{ij},s^{\prime}] to agents agents ii and jj, respectively, does not decrease the social welfare (Lemma 6, part (b)). Though, at the same time, it does provide an allocation that uses more cut points from PP and, hence, contradicts the choice of S\mathcal{S}. On the other hand, if pij<sp_{ij}<s, then we can assign the entire interval [s,s′][s,s^{\prime}] to agent jj. As before, the reassignment does not decrease the social welfare. However, it does decrease the cardinality of the set difference between the cut points and PP. This contradicts the selection criterion of S\mathcal{S}.

Hence, the cut points of S\mathcal{S} satisfy {si}i⊆P\{s_{i}\}_{i}\subseteq P and the stated claim follows. ⊓\sqcap⊔\sqcup

We now present the main result of this section.

Proof Given a cake-division instance C\mathcal{C} with MLRP, write S∗={S1∗,S2∗,…,Sn∗}\mathcal{S}^{*}=\{S^{*}_{1},S^{*}_{2},\dots,S^{*}_{n}\} to denote the allocation identified in Lemma 7; in particular, S∗\mathcal{S}^{*} is a social welfare maximizing allocation whose cut points {0=s0∗,s1∗,…,sn∗=1}\{0=s^{*}_{0},s^{*}_{1},\dots,s^{*}_{n}=1\} belong to the set of switching points PP.

For a precision parameter η>0\eta>0 and for each switching point pij∈Pp_{ij}\in P, we invoke Corollary 1 to find p^ij∈\widehat{p}_{ij}\in with the property that ∣p^ij−pij∣≤ηnλ|\widehat{p}_{ij}-p_{ij}|\leq\frac{\eta}{n\lambda}. Write P^\widehat{P} to denote the set of these estimates, P^≔{p^ij\mathchar581≤i<j≤n}∪{0,1}\widehat{P}\coloneqq\{\widehat{p}_{ij}\mathrel{\mathop{\mathchar 58\relax}}1\leq i<j\leq n\}\cup\{0,1\}.As in the case of PP, we include and 11 in P^\widehat{P} for ease of presentation. Applying Corollary 1 to each p^ij\widehat{p}_{ij}, we get that the set P^\widehat{P} can be computed in O(poly(n,log⁡λ,log⁡1η))\mathcal{O}({\rm poly}(n,\log\lambda,\log\frac{1}{\eta})) time.

For the optical allocation S∗={S1∗,…,Sn∗}\mathcal{S}^{*}=\{S^{*}_{1},\ldots,S^{*}_{n}\} we have Si∗=[si−1∗,si∗]S^{*}_{i}=[s^{*}_{i-1},s^{*}_{i}] for all i∈[n]i\in[n]. Write allocation S^≔{S^1,S^2,…,S^n}\widehat{\mathcal{S}}\coloneqq\{\widehat{S}_{1},\widehat{S}_{2},\ldots,\widehat{S}_{n}\}, where interval S^i=[s^i−1,s^i]\widehat{S}_{i}=[\widehat{s}_{i-1},\widehat{s}_{i}] is assigned to agent i∈[n]i\in[n]. The cut points of allocations S^\widehat{\mathcal{S}} are contained in P^\widehat{P}. Also, given that S∗\mathcal{S}^{*} conforms to the MLRP order, so does S^\widehat{\mathcal{S}}.

Therefore, there exists an allocation S^\widehat{\mathcal{S}} with the properties that (i) S^\widehat{\mathcal{S}} has near-optimal social welfare, (ii) cut points of S^\widehat{\mathcal{S}} are contained in P^\widehat{P}, and (iii) S^\widehat{\mathcal{S}} conforms to the MLRP order. To complete the proof of the theorem we will show that, among all allocations that satisfy properties (ii) and (iii), we can find one that maximizes social welfare. We accomplish this algorithmic result by a simple dynamic program.

Recall that cardinality of the set P^\widehat{P} is O(n2)\mathcal{O}(n^{2}) and this set can be computed in O(poly(n,log⁡λ,log⁡1η))\mathcal{O}({\rm poly}(n,\log\lambda,\log\frac{1}{\eta})) time using eval queries. We index the elements of the computed set P^={p^t}t\widehat{P}=\{\widehat{p}_{t}\}_{t} such that 0=p^0<p^1<…<p^∣P^∣=10=\widehat{p}_{0}<\widehat{p}_{1}<\ldots<\widehat{p}_{|\widehat{P}|}=1

For each k∈[n]k\in[n] and 1≤t≤∣P^∣1\leq t\leq|\widehat{P}|, we write M(k,t)M(k,t) to denote the maximum social welfare that one can achieve by allocating the interval [0,p^t][0,\widehat{p}_{t}] among the first kk agents (in order).By convention, the agents are indexed following the MLRP order.

The following recursive equation for M(k,t)M(k,t) gives us the desired dynamic program

Overall, we can find an allocation with social welfare η\eta close to the optimal in time O(poly(n,log⁡λ,log⁡1η))\mathcal{O}({\rm poly}(n,\log\lambda,\log\frac{1}{\eta})). Since the precision parameter η\eta can be driven exponentially close to zero, in time that is polynomial in log⁡1η\log\frac{1}{\eta} (i.e., in the bit complexity of η\eta) the stated claim follows. ⊓\sqcap⊔\sqcup

Egalitarian Welfare

This section presents an algorithm for maximizing egalitarian welfare in cake-division instances with MLRP (Theorem 4).

Here, we use the fact that the value densities are nonnegative. This establishes the stated claim. ⊓\sqcap⊔\sqcup

We now establish the main result for egalitarian welfare. See 4

With a precision parameter η>0\eta>0 in hand, we perform binary search over integer multiples of η{\eta}, i.e., over the set {kη}k=01/η\left\{k\eta\right\}_{k=0}^{1/\eta}. Recall that the values of the agents are normalized and, hence, τ∗∈\tau^{*}\in. Also, τ=0\tau=0 is feasible, while τ=1\tau=1 is infeasible.,

This switch in feasibility (between and 11) along with property (P), imply that there exists a unique index k0∈{0,1,…,1/η}k_{0}\in\{0,1,\ldots,1/\eta\} with the property that k0ηk_{0}\eta is feasible and (k0+1)η(k_{0}+1)\eta is infeasible. We can identify k0k_{0} by O(log⁡(1/η))\mathcal{O}\left(\log(1/\eta)\right) iterations of binary search.

The relevant observation here is that τ∗∈[k0η,(k0+1)η]\tau^{*}\in[k_{0}\eta,(k_{0}+1)\eta]. Indeed, the optimal value τ∗\tau^{*} is feasible: R∗={R1∗,…,Rn∗}\mathcal{R}^{*}=\{R^{*}_{1},\ldots,R^{*}_{n}\} conforms to the MLRP order and vi(Ri∗)≥τ∗v_{i}(R^{*}_{i})\geq\tau^{*} for all i∈[n]i\in[n]. Hence, property (P) ensures that τ∗\tau^{*} cannot be greater than the infeasible target (k0+1)η(k_{0}+1)\eta. Furthermore, the optimality of τ∗\tau^{*} gives us τ∗≥k0η\tau^{*}\geq k_{0}\eta.

For the runtime analysis, note that the binary search finds the desired index k0k_{0} in O(poly(n,log⁡1η))\mathcal{O}({\rm poly}(n,\log\frac{1}{\eta})) time. The dependency on λ\lambda stems from the fact that the bit-complexity of the output (i.e., of the computed cut points) can be O(log⁡ηλ)\mathcal{O}\left(\log\frac{\eta}{\lambda}\right).

Overall, these arguments show that we can find an allocation R^\widehat{\mathcal{R}} with egalitarian welfare η\eta close to the optimal in time O(poly(n,log⁡λ,log⁡1η))\mathcal{O}({\rm poly}(n,\log\lambda,\log\frac{1}{\eta})). The precision parameter η\eta can be driven exponentially close to zero in time that is polynomial in log⁡1η\log\frac{1}{\eta} and, hence, the stated claim follows. ⊓\sqcap⊔\sqcup

Nash Social Welfare

This section presents an FPTAS for maximizing Nash social welfare in cake-division instances with MLRP.

Proof Given a cake-division instance C\mathcal{C} with MLRP, write A∗={A1∗,A2∗,…,An∗}\mathcal{A}^{*}=\{A^{*}_{1},A^{*}_{2},\dots,A^{*}_{n}\} to denote an allocation that maximizes the Nash social welfare in C\mathcal{C}. Arunachaleswaran et al. [ABKR19] have shown that the bundles in any Nash optimal allocation A∗={A1∗,A2∗,…,An∗}\mathcal{A}^{*}=\{A^{*}_{1},A^{*}_{2},\dots,A^{*}_{n}\} satisfy vi(Ai∗)≥14vi(Aj∗)v_{i}(A^{*}_{i})\geq\frac{1}{4}v_{i}(A^{*}_{j}) for all i,j∈[n]i,j\in[n]; this result holds even in the absence of MLRP. For a fixed agent i∈[n]i\in[n], we sum the inequalities vi(Ai∗)≥14vi(Aj∗)v_{i}(A^{*}_{i})\geq\frac{1}{4}v_{i}(A^{*}_{j}) over all j∈[n]j\in[n] to obtain vi(Ai∗)≥14nv_{i}(A^{*}_{i})\geq\frac{1}{4n}. Recall that vi(0,1)=1v_{i}(0,1)=1 for all i∈[n]i\in[n].

We round each cut point in the optimal allocation A∗\mathcal{A}^{*} to its closest point in the collection {ci}i=0N\{c_{i}\}_{i=0}^{N}. This leads us to another allocation, A^={A^1,A^2,…,A^n}\widehat{\mathcal{A}}=\{\widehat{A}_{1},\widehat{A}_{2},\ldots,\widehat{A}_{n}\}, wherein interval A^i\widehat{A}_{i} is assigned to agent i∈[n]i\in[n]. By construction, the cut points of A^\widehat{\mathcal{A}} are contained in the set {ct}t=0N\{c_{t}\}_{t=0}^{N}. Also, since A∗\mathcal{A}^{*} conforms to the MLRP order (Lemma 5), so does A^\widehat{\mathcal{A}}. Next we show that the Nash social welfare of A^\widehat{\mathcal{A}} is comparable to that of A∗\mathcal{A}^{*}. For all i∈[n]i\in[n], we have

Therefore, there exists an allocation A^\widehat{\mathcal{A}} with the properties that (i) A^\widehat{\mathcal{A}} has Nash social welfare at least (1−ε)(1-\varepsilon) times the optimal (Nash social welfare), (ii) the cut points of A^\widehat{\mathcal{A}} are contained in the set {ct}t\{c_{t}\}_{t}, and (iii) A^\widehat{\mathcal{A}} conforms to the MLRP order. To complete the proof of the theorem we will show that, among all the allocations that satisfy (ii) and (iii), we can find (via a dynamic program) one that maximizes the Nash social welfare.

For t∈[N]t\in[N] and k∈[n]k\in[n], we write H(k,t)H(k,t) to denote the optimal Nash product (i.e., the product of valuations) that one can achieve by allocating the interval [0,ct][0,c_{t}] among the first kk agents (in order).

Therefore, we can find an allocation with Nash social welfare at least (1−ε)(1-\varepsilon) times the optimal in O(poly(n,1/ε,log⁡λ))\mathcal{O}\left({\rm poly}\left(n,1/\varepsilon,\log\lambda\right)\right) time; the dependency on log⁡λ\log\lambda stems from the fact that the bit-complexity of the output (i.e., of the computed cut points) can be O(log⁡ελ)\mathcal{O}\left(\log\frac{\varepsilon}{\lambda}\right). Overall, we get that maximizing Nash social welfare admits an FPTAS under MLRP. ⊓\sqcap⊔\sqcup

Conclusion and Future Work

The current work studies algorithmic aspects of contiguous cake division under the monotone likelihood ratio property. The scope of this property ensures that the developed algorithms are applicable in various cake-division settings. We also note that while under MLRP the value densities must have full support, the developed framework is somewhat robust to this requirement. For example, our results extend to the class of non-full-support value densities considered in [AFG+17]. In particular, Alijani et al. [AFG+17] established that a contiguous envy-free cake division can be efficiently computed if every agent uniformly values a single interval and these intervals satisfy an ordering property. Appendix E shows that here one can modify the value densities to a small degree and obtain MLRP (with full support). Hence, applying our results, one can efficiently compute an allocation with arbitrary small envy in the modified instance and, hence, also in the original one. Generalizing such ideas to address, say, value densities that bear first-order stochastic dominance is an interesting direction of future work.

For instances with MLRP, finding an allocation that maximizes various welfare notions among the set of envy-free allocations is an important thread for future work. Another relevant notion of fairness in the context of cake cutting is that of a perfect division [Alo87]. In such a division D={D1,D2,…,Dn}\mathcal{D}=\{D_{1},D_{2},\ldots,D_{n}\}, each agent i∈[n]i\in[n] values every piece at 1/n1/n, i.e., vi(Dj)=1/nv_{i}(D_{j})=1/n for all i,j∈[n]i,j\in[n]. In contrast to the other solution concepts considered in the present paper, perfect divisions are not guaranteed to exist under the contiguity requirement; a perfect allocation might not exist even with MLRP (Appendix F). However, Alon [Alo87] has shown that a perfect division with n(n−1)n(n-1) cuts always exists. Perfect cake divisions are particularly useful since they lead to truthful mechanisms for cake division [MT10, CLPP13]. Hence, developing efficient algorithms to find (noncontiguous) perfect divisions under MLRP is a relevant thread for future work.

More broadly, it would be interesting to identify tractable classes through MLRP in other computational social choice contexts, such as discrete fair division and voting.

Acknowledgements

We thank Manjunath Krishnapur for helpful discussions and references. Siddharth Barman gratefully acknowledges the support of a Ramanujan Fellowship (SERB - SB/S2/RJN-128/2015) and a Pratiksha Trust Young Investigator Award. Nidhi Rathi’s research is generously supported by an IBM PhD Fellowship.

References

Appendix A Implications of MLRP

In this section we restate and prove Lemma 1.

Proof Given that fif_{i} and fjf_{j} bear MLRP, we will first prove that they satisfy property (i). For b≤cb\leq c, we have

Recall that MLRP value densities are, by definition, positively valued, fi(x)>0f_{i}(x)>0 for all x∈x\in. Therefore, equation (7) gives us fj(x)≤ fj(b)fi(b) fi(x)f_{j}(x)\leq\ \frac{f_{j}(b)}{f_{i}(b)}\ f_{i}(x) for all x∈[a,b]x\in[a,b]. Integrating we obtain ∫abfj(x)dx≤∫abfj(b)fi(b) fi(x)dx\int\limits_{a}^{b}f_{j}(x)dx\leq\int\limits_{a}^{b}\frac{f_{j}(b)}{f_{i}(b)}\ f_{i}(x)dx and, hence,Recall that the integral of a positive function is positive.

Starting with equation (7) and integrating over the interval [c,d][c,d], we can also establish the following equality

Equations (8) and (9) lead to property (i):

Next we will prove that properties (i) and (ii) are equivalent. Since fif_{i} and fjf_{j} satisfy property (i), the equivalence of properties (i) and (ii) will imply that they satisfy property (ii) as well; thereby completing the proof.

To establish that property (i) implies property (ii), we instantiate (i) over the intervals [a,x][a,x] and [x,b][x,b]: ∫axfj∫axfi≤∫xbfj∫xbfi\frac{\int_{a}^{x}f_{j}}{\int_{a}^{x}f_{i}}\leq\frac{\int_{x}^{b}f_{j}}{\int_{x}^{b}f_{i}}. Cross multiplying the termsRecall that the value densities are strictly positive. and adding one to both sides of the inequality, gives us ∫axfj∫xbfj+1≤∫axfi∫xbfi+1\frac{\int_{a}^{x}f_{j}}{\int_{x}^{b}f_{j}}+1\leq\frac{\int_{a}^{x}f_{i}}{\int_{x}^{b}f_{i}}+1. Simplifying further we obtain ∫axfj+∫xbfj∫xbfj≤∫axfi+∫xbfi∫xbfi\frac{\int_{a}^{x}f_{j}+\int_{x}^{b}f_{j}}{\int_{x}^{b}f_{j}}\leq\frac{\int_{a}^{x}f_{i}+\int_{x}^{b}f_{i}}{\int_{x}^{b}f_{i}}. This gives us the desired bound

For the reverse direction, i.e., (ii) implies (i), we begin by cross multiplying the terms in inequality (10) and subtracting one from both sides yield ∫abfj∫xbfj−1≤∫abfi∫xbfi−1\frac{\int_{a}^{b}f_{j}}{\int_{x}^{b}f_{j}}-1\leq\frac{\int_{a}^{b}f_{i}}{\int_{x}^{b}f_{i}}-1. This inequality simplifies to ∫axfj∫xbfj≤∫axfi∫xbfi\frac{\int_{a}^{x}f_{j}}{\int_{x}^{b}f_{j}}\leq\frac{\int_{a}^{x}f_{i}}{\int_{x}^{b}f_{i}}. That is, we obtain property (i) for intervals [a,x][a,x] and [x,b][x,b]. Reapplying this bound (with aa, xx, and bb set appropriately) shows that (ii) implies (i). This completes the proof. ⊓\sqcap⊔\sqcup

A.2 Efficiently Finding the MLRP Order

To prove this claim, we first state a well-known result (in Lemma 8) that MLRP implies first-order stochastic dominance; we provide a proof here for completeness. Recall that, given two probability density functions fif_{i} and fjf_{j} over $,density, densityf_{j}issaidtohavefirst−orderstochasticdominanceoveris said to have first-order stochastic dominance overf_{i}iffiff\int_{t}^{1}f_{j}(x)dx\geq\int_{t}^{1}f_{i}(x)dx,forall, for allt\in,andthereexistsatleastone, and there exists at least onet^{\prime}\insuchthatsuch that\int_{t^{\prime}}^{1}f_{j}(x)dx>\int_{t^{\prime}}^{1}f_{i}(x)dx$.

Let fif_{i} and fjf_{j} be two (ordered) value-density functions that satisfy the monotone likelihood ratio property: for every 0≤x≤y≤10\leq x\leq y\leq 1 we have fj(x)fi(x)≤fj(y)fi(y)\frac{f_{j}(x)}{f_{i}(x)}\leq\frac{f_{j}(y)}{f_{i}(y)}. Then, the density fjf_{j} has first-order stochastic dominance over fif_{i}.

Proof Given that fif_{i} and fjf_{j} bear MLRP, we consider property (ii) of Lemma 1, with a=0a=0, b=1b=1, and x=tx=t, for any 0≤t≤10\leq t\leq 1, to obtain ∫t1fi≤∫t1fj{\int\limits_{t}^{1}f_{i}}\leq\int\limits_{t}^{1}f_{j}. Here, we use the fact that the valuations are normalized, ∫01fi=∫01fj=1\int_{0}^{1}f_{i}=\int_{0}^{1}f_{j}=1.

Furthermore, since the two densities are distinct, there exists a point t′∈t^{\prime}\in such that ∫t′1fj\int_{t^{\prime}}^{1}f_{j} is not equal to ∫t′1fi\int_{t^{\prime}}^{1}f_{i}. For such a point t′t^{\prime}, a strict inequality must hold, ∫t′1fj>∫t′1fi\int_{t^{\prime}}^{1}f_{j}>\int_{t^{\prime}}^{1}f_{i}. ⊓\sqcap⊔\sqcup

Note that, since both MLRP and first-order stochastic dominance are transitive properties, the above lemma directly extends to a collection of probability density functions.

Next we prove that, given a cake-division instance, with the promise that the underlying value densities satisfy MLRP, one can find the MLRP order efficiently.

Let C\mathcal{C} be a cake-division instance in which the value-density functions satisfy the monotone likelihood ratio property. Then, in the Robertson-Webb query model, we can find the MLRP order of C\mathcal{C} in polynomial time.

The first-order stochastic dominance between fjf_{j} and fif_{i} also ensures that there exists a point t′∈t^{\prime}\in such that

We consider two complementary cases (i)(i) t′≤1/2t^{\prime}\leq 1/2 and (ii)(ii) t′>1/2t^{\prime}>1/2. In both of these cases we assume, towards a contradiction, that ∫1/21fj=∫1/21fi\int_{1/2}^{1}f_{j}=\int_{1/2}^{1}f_{i}.

Case (i): t′≤1/2t^{\prime}\leq 1/2. Note that in this case equation (11) expands to ∫t′1/2fj+∫1/21fj>∫t′1/2fi+∫1/21fi\int_{t^{\prime}}^{1/2}f_{j}+\int_{1/2}^{1}f_{j}>\int_{t^{\prime}}^{1/2}f_{i}+\int_{1/2}^{1}f_{i}. Since, ∫1/21fj=∫1/21fi\int_{1/2}^{1}f_{j}=\int_{1/2}^{1}f_{i}, we have ∫t′1/2fj>∫t′1/2fi\int_{t^{\prime}}^{1/2}f_{j}>\int_{t^{\prime}}^{1/2}f_{i}.

On the other hand, applying Lemma 1, property (i), to intervals [t′,1/2][t^{\prime},1/2] and [1/2,1][1/2,1] gives us ∫t′1/2fj∫t′1/2fi≤∫1/21fj∫1/21fi=1\frac{\int_{t^{\prime}}^{1/2}f_{j}}{\int_{t^{\prime}}^{1/2}f_{i}}\leq\frac{\int_{1/2}^{1}f_{j}}{\int_{1/2}^{1}f_{i}}=1. This inequality leads to the desired contradiction, ∫t′1/2fj≤∫t′1/2fi\int_{t^{\prime}}^{1/2}f_{j}\leq\int_{t^{\prime}}^{1/2}f_{i}.

Case (ii): t′>1/2t^{\prime}>1/2. Here, equation (11) and the assumed equality ∫1/21fj=∫1/21fi\int_{1/2}^{1}f_{j}=\int_{1/2}^{1}f_{i} give us ∫1/2t′fj<∫1/2t′fi\int_{1/2}^{t^{\prime}}f_{j}<\int_{1/2}^{t^{\prime}}f_{i}.

Note that (due to normalization) we have ∫01/2fj=∫01/2fi\int_{0}^{1/2}f_{j}=\int_{0}^{1/2}f_{i}. This equality contradicts an application of Lemma 1, property (i), to the intervals [0,1/2][0,1/2] and [1/2,t′][1/2,t^{\prime}]

A.3 Lipschitzness of Cut and Eval Queries under MLRP

Appendix B Missing Proofs from Section 5

B.2 Proof of Lemma 4

This section restates Lemma 4 and shows that BinSearch (Algorithm 1) satisfies this claim.

Appendix C Distribution Families with MLRP

This section highlights that various well-studied distribution families bear MLRP. Recall, that two probability density functions fif_{i} and fjf_{j} are said to satisfy MLRP iff, for every x≤yx\leq y in the domain, we have

Binomial Polynomials: The following proposition shows that every pair of binomial polynomials bear MLRP over $,i.e.,wehaveatotalorderoverthisfamilyofdensityfunctionswithrespecttoMLRP.Asmentionedpreviously,MLRPispreservedunderscalingand,hence,herewedonothavetoexplicitlyenforcenormalization.Hence,theresultsdevelopedintheworkholdforcake−divisioninstancesinwhichthevaluedensitiesarebinomial.Also,settingtheexponentparameters, i.e., we have a total order over this family of density functions with respect to MLRP.As mentioned previously, MLRP is preserved under scaling and, hence, here we do not have to explicitly enforce normalization. Hence, the results developed in the work hold for cake-division instances in which the value densities are binomial. Also, setting the exponent parameterss=1andandt=0$ in this proposition, we observe that linear functions form a special case of binomial polynomials.

With integer exponents s>ts>t, let fi(x)=aixs+bixtf_{i}(x)=a_{i}x^{s}+b_{i}x^{t} and fj(x)=ajxs+bjxtf_{j}(x)=a_{j}x^{s}+b_{j}x^{t} be two binomial polynomials. Then, fif_{i} and fjf_{j} bear MLRP iff aibj−ajbi≤0a_{i}b_{j}-a_{j}b_{i}\leq 0.

Proof For any two points x,y∈x,y\in, such that x≤yx\leq y, the MLRP condition for binomials corresponds to the following inequality

Since fif_{i} and fjf_{j} constitute value densities with full support, the values of these functions are positive over the cake. Hence, the previous equation can be rewritten as

Finally, we rewrite the last inequality to obtain that for binomials MLRP is equivalent to

Gaussian distributions: The next proposition shows that Gaussian distributions with different means, but the same variance, bear MLRP.

Let fif_{i} and fjf_{j} be two Gaussian density functions with the same variance σ2\sigma^{2} and means μi≤μj\mu_{i}\leq\mu_{j}, respectively. Then, fif_{i} and fjf_{j} satisfy MLRP.

Write g(x)≔12σ2((x−μi)2−(x−μj)2)g(x)\coloneqq\frac{1}{2\sigma^{2}}\left((x-\mu_{i})^{2}-(x-\mu_{j})^{2}\right) and note that the derivate of this function g′(x)=(μj−μi)σ2≥0g^{\prime}(x)=\frac{(\mu_{j}-\mu_{i})}{\sigma^{2}}\geq 0 if and only if μi≤μj\mu_{i}\leq\mu_{j}. Hence, gg is an increasing function in xx—and so is exp⁡(g(x)){\exp}\left(g(x)\right)—iff μi≤μj\mu_{i}\leq\mu_{j}. That is, fif_{i} and fjf_{j} bear MLRP iff μi≤μj\mu_{i}\leq\mu_{j}. ⊓\sqcap⊔\sqcup

An analogous result holds for Gaussians with mean zero, but distinct variances.

Appendix D Bit Complexity of Cake Division

This section provides a cake-division instance (with MLRP) in which the unique envy-free allocation has an irrational cut point. Notably the parameters that specify the value densities in this instance are rational. This example implies that, in general, one cannot expect an efficient algorithm that outputs an exact envy-free allocation. That is, when considering cake-division algorithms with bounded bit complexity, a precision loss is inevitable.

We will also show, through the example, that to obtain a nontrivial bound on the envy (between the agents) the bit complexity of the output has to be Ω(log⁡λ)\Omega(\log\lambda); here λ≥1\lambda\geq 1 is the Lipschitz constant of the cut and eval queries. Hence, a runtime dependence of log⁡λ\log\lambda is unavoidable as well.

Note that this function can be given as input using only rational parameters. In addition, ff is discontinuous at 1−1λ1-\frac{1}{\lambda}; recall that our results require the value densities to be integrable, and not necessarily continuous.

Here, the values of the agents are normalized, ∫01f(x)dx=1\int_{0}^{1}f(x)dx=1. Also, the following bounds hold for the density: 0<λ3(λ−1)≤f(x)≤λ0<\frac{\lambda}{3(\lambda-1)}\leq f(x)\leq\lambda for all x∈x\in. Therefore, Proposition 3 implies that the cut and eval queries in this instance are 3λ3\lambda-Lipschitz. Since in this instance the three agents have identical value densities (with full support over $),thereexistsauniqueenvy−freeallocationwhereineachagentreceivesanintervalofvalueexactlyequalto), there exists a unique envy-free allocation wherein each agent receives an interval of value exactly equal to1/3.Write. Write0=x^{*}_{0}\leq x^{*}_{1}\leq x^{*}_{2}\leq x^{*}_{3}=1todenotethecutpointsofthis(unique)envy−freeallocation;inparticular,agentto denote the cut points of this (unique) envy-free allocation; in particular, agenti\inreceivesthereceives theithintervalth interval[x^{*}_{i-1},x^{*}_{i}]$.

First, we will show that x2∗x^{*}_{2} is irrational. Note that the interval [0,1−1/λ][0,1-1/\lambda] is of value 1/31/3: ∫01−1λf(x)dx=13\int_{0}^{1-\frac{1}{\lambda}}f(x)dx=\frac{1}{3}. Hence, the first cut point x1∗=1−1λx^{*}_{1}=1-\frac{1}{\lambda}. The second cut point x2∗x^{*}_{2} now lies in the interval ([x1∗,1][x^{*}_{1},1]) of width 1λ\frac{1}{\lambda} and it must satisfy ∫1−1λx2∗f(x)dx=13\int_{1-\frac{1}{\lambda}}^{x^{*}_{2}}f(x)dx=\frac{1}{3}. Using the definition of ff in this range, we get that x2∗x^{*}_{2} is a solution of the following quadratic equation

We next establish a lower bound on the bit complexity of the output. Consider, in the above-mentioned instance, any allocation I={I1,I2,I3}\mathcal{I}=\{I_{1},I_{2},I_{3}\} wherein the envy between the agents is, say, less than η=1/2\eta=1/2, i.e., vi(Ii)≥vi(Ij)−1/2v_{i}(I_{i})\geq v_{i}(I_{j})-1/2 for all i,j∈[n]i,j\in[n].Here, the choice of η=1/2\eta=1/2 is essentially for ease of exposition; by scaling down the density in the range [0,(1−1/λ)]\left[0,\left(1-{1}/{\lambda}\right)\right], we can drive η\eta close to one. For any such allocation I\mathcal{I} one of the cut points must lie in [(1−1/λ),1][\left(1-{1}/{\lambda}\right),1]. Indeed, this interval is of value 2/32/3 to each agent, ∫(1−1/λ)1f=2/3\int_{\left(1-1/\lambda\right)}^{1}f=2/3. The bit complexity of such a cut point is Ω(log⁡λ)\Omega(\log\lambda). Hence, in general, the bit complexity of the any algorithm that finds an allocation (equivalently, outputs cut points) with bounded envy is Ω(log⁡λ)\Omega(\log\lambda).

Welfare-maximizing allocations induced by irrational cuts: One can also construct cake-division instances (with rational and MLRP value densities) in which the welfare-maximizing allocations have irrational cut points.

In particular, consider cake division between two agents with identical value densities, f(x)=x+1/2f(x)=x+1/2. Here, the (unique) allocation I={I1,I2}\mathcal{I}=\{I_{1},I_{2}\} that maximizes egalitarian welfare consists of the intervals I1=[0,5−12]I_{1}=\left[0,\frac{\sqrt{5}-1}{2}\right] and I2=[5−12,1]I_{2}=\left[\frac{\sqrt{5}-1}{2},1\right]. Allocation I\mathcal{I} is also the (unique) equitable, proportional, and envy-free allocation in this instance.

In addition, consider a cake-division instance with two agents and the following value densities: f1(x)=1f_{1}(x)=1 and f2(x)=3x2f_{2}(x)=3x^{2}. Note that these densities satisfy MLRP. In this instance, the (unique) social welfare maximizing allocation S={S1,S2}\mathcal{S}=\{S_{1},S_{2}\} is obtained by an irrational cut; specifically, S1=[0,13]S_{1}=\left[0,\frac{1}{\sqrt{3}}\right] and S2=[13,1]S_{2}=\left[\frac{1}{\sqrt{3}},1\right]. The cut point 13\frac{1}{\sqrt{3}} is the switching point (as defined in Section 7) between the two densities f1f_{1} and f2f_{2}.

Appendix E Robustness of MLRP

Let K=⟨[n],{fi}i∈[n]⟩\mathcal{K}=\langle[n],\{f_{i}\}_{i\in[n]}\rangle denote a cake-division instance with such value densities and note that, for each i∈[n]i\in[n],

Indeed, the value densities in K\mathcal{K} do not have full support over the cake. However, we will show that we can perturb these densities fif_{i}s (in a structured manner) to obtain an instance K^=⟨[n],{f^i}i∈[n]⟩\widehat{\mathcal{K}}=\langle[n],\{\widehat{f}_{i}\}_{i\in[n]}\rangle such that (i) the value densities f^i\widehat{f}_{i}s in K^\widehat{\mathcal{K}} have full support and bear MLRP (Claim 1) and (ii) any envy-free allocation in K^\widehat{\mathcal{K}} forms an envy-free allocation, up to a small precision loss, in the original instance K\mathcal{K} (Claim 2).

Applying Proposition 3, we obtain that the Lipschitz constant of cut and eval queries in K^\widehat{\mathcal{K}} is λ=Hn max⁡1≤i≤n hi\lambda=H^{n}\ \max_{1\leq i\leq n}\ h_{i}. The next claim shows that the the value densities in K^\widehat{\mathcal{K}} satisfy MLRP.

Let K=⟨[n],{fi}i∈[n]⟩\mathcal{K}=\langle[n],\{f_{i}\}_{i\in[n]}\rangle be a cake-division instance in which the value densities satisfy equation (15) and the ordering property OP. Then, for parameter H≥1{H}\geq 1, the value densities f^i\widehat{f}_{i}s (as defined in equation (16)) satisfy MLRP.

Finally, for x≥rjx\geq r_{j}, we note that did_{i} and djd_{j} increase synchronously. Also, for points x∈[rj,1]x\in[r_{j},1] we have cj(x)=jc_{j}(x)=j and ci(x)=ic_{i}(x)=i. Therefore, in this range the likelihood ratio stays constant at f^j(rj)f^i(rj)\frac{\widehat{f}_{j}(r_{j})}{\widehat{f}_{i}(r_{j})}. Overall, we obtain the monotonicity of the likelihood ratio and the MLRP guarantee follows. ⊓\sqcap⊔\sqcup

Given a cake-division instance K=⟨[n],{fi}i∈[n]⟩\mathcal{K}=\langle[n],\{f_{i}\}_{i\in[n]}\rangle in which the value densities satisfy equation (15) and the ordering property OP. Let K^=⟨[n],{f^i}i∈[n]⟩\widehat{\mathcal{K}}=\langle[n],\{\widehat{f}_{i}\}_{i\in[n]}\rangle be the cake-division instance defined above, with parameter H≔2η max⁡1≤i≤nhiH\coloneqq\frac{2}{\eta}\ \max_{1\leq i\leq n}h_{i}, and suppose that allocation I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\} is envy-free up to an additive factor of η\eta in K^\widehat{\mathcal{K}} (i.e., vi(Ii)≥vi(Ij)−ηv_{i}(I_{i})\geq v_{i}(I_{j})-\eta for all i,j∈[n]i,j\in[n]). Then, allocation I\mathcal{I} envy-free up to 2η2\eta in K\mathcal{K}.

Hence, an allocation I={I1,…,In}\mathcal{I}=\{I_{1},\ldots,I_{n}\} that is envy-free up to an additive factor of η\eta in K^\widehat{\mathcal{K}} (i.e., v^i(Ii)≥v^i(Ij)−η\widehat{v}_{i}(I_{i})\geq\widehat{v}_{i}(I_{j})-\eta for all i,j∈[n]i,j\in[n]) is envy-free up to 2η2\eta in K\mathcal{K}. Specifically, for any i,j∈[n]i,j\in[n],

That is, I{\mathcal{I}} is envy-free, up to 2η2\eta precision, in K\mathcal{K}. This completes the proof. ⊓\sqcap⊔\sqcup

Since the constructed instance K^\widehat{\mathcal{K}} satisfies MLRP, we can use Algorithm 1 (Section 5) to efficiently compute, up to an arbitrary precision, an envy-free allocation I\mathcal{I} in the instance K^\widehat{\mathcal{K}}. The previous claim ensures that I\mathcal{I} is an envy-free allocation (up to an arbitrary precision) in the original instance K\mathcal{K} as well. This observation highlights the fact that the ideas developed in this work are somewhat robust and extend to other value-density settings.

Appendix F Nonexistence of Contiguous Perfect Divisions under MLRP

A cake division D={D1,D2,…,Dn}\mathcal{D}=\{D_{1},D_{2},\dots,D_{n}\} (consisting of connected or disconnected pieces) is said to be perfect if all the agents agree on the value of every piece, i.e., vi(Dj)=1/nv_{i}(D_{j})=1/n for all i,j∈[n]i,j\in[n].

Perfect divisions are not guaranteed to exist under the contiguity requirement, i.e., there are cake-division instances that do not admit perfect allocations. In this section, we will show that the nonexistence continues to hold even with MLRP. That is, we will provide a cake-division instance with MLRP that does not admit a perfect allocation.

The agents’ values for the cake $arenormalized,are normalized,\int_{0}^{1}f_{1}(x)dx=\int_{0}^{1}f_{2}(x)dx=1.Also,thelikelihoodratioof. Also, the likelihood ratio off_{1}andandf_{2}$ satisfies

Since 2−αα>1−α1+α\frac{2-\alpha}{\alpha}>\frac{1-\alpha}{1+\alpha}, for all α∈(0,1)\alpha\in(0,1), the likelihood ratio is nondecreasing over $$ and, hence, the densities bear MLRP.

We assume, towards a contradiction that there exists a perfect allocation in C\mathcal{C}. That is, there exists a point x∈x\in such that v1(0,x)=v2(0,x)=1/2v_{1}(0,x)=v_{2}(0,x)=1/2, and v1(x,1)=v2(x,1)=1/2v_{1}(x,1)=v_{2}(x,1)=1/2. The fact that the value of intervals [0,x][0,x] and [x,1][x,1] is equal to 1/21/2 ensures that the point xx cannot be or 11. We consider two complementary and exhaustive cases: Case (i) 0<x≤(1−α)0<x\leq(1-\alpha) and Case (ii) (1−α)<x≤1(1-\alpha)<x\leq 1. The analysis is similar in both the cases. Hence, we only address Case (i) and omit the analysis for Case (ii).

For the first interval [0,x][0,x] (with 0<x≤(1−α)0<x\leq(1-\alpha)), we have v1(0,x)=x(1+α)v_{1}(0,x)=x(1+\alpha) and v2(0,x)=x(1−α)v_{2}(0,x)=x(1-\alpha). However, for any α∈(0,1)\alpha\in(0,1) and x>0x>0, the following strict inequality holds: v1(0,x)=x(1+α)>x(1−α)=v2(0,x)v_{1}(0,x)=x(1+\alpha)>x(1-\alpha)=v_{2}(0,x). This leads to a contradiction and proves that C\mathcal{C} does not admit a perfect division with connected pieces.

Even though we might not have a perfect allocation, the work of Alon [Alo87] proves that a perfect division with n(n−1)n(n-1) cuts always exists. Hence, in the above-mentioned instance C\mathcal{C} with 22 agents, 22 cuts should suffice to form a perfect division. In particular, we note that the following division D∗={D1∗,D2∗}\mathcal{D}^{*}=\{D^{*}_{1},D^{*}_{2}\} is perfect in C\mathcal{C}; here D1∗=[12−α2,1−α2]D^{*}_{1}=\left[\frac{1}{2}-\frac{\alpha}{2},1-\frac{\alpha}{2}\right] and D2∗=[0,12−α2]∪[1−α2,1]D^{*}_{2}=\left[0,\frac{1}{2}-\frac{\alpha}{2}\right]\cup\left[1-\frac{\alpha}{2},1\right]. Here,

That is, both the agents value the piece D2∗D^{*}_{2} at 1/21/2. Since the valuations are normalized, we additionally have v1(D1∗)=v2(D1∗)=1/2v_{1}(D^{*}_{1})=v_{2}(D^{*}_{1})=1/2. This shows that D∗\mathcal{D}^{*} is a perfect division (with disconnected pieces) in C\mathcal{C}.