Constrained Assortment Optimization Under the Paired Combinatorial Logit Model: Max-Dicut LP Rounding Gives a 0.5-Approximation for Cardinality-Constrained AssortmentResearch Paper
Motivation
A retailer choosing which products to offer (an assortment) must predict how customers substitute between products. Discrete choice models make this prediction, and the paired combinatorial logit (PCL) model is a generalization of the multinomial and nested logit models in which every pair of products forms a nest, so that the correlation between the utilities of any two products can be set separately. Assortment optimization under the PCL model is NP-hard even without constraints, and practical instances carry constraints on the number of products displayed.
R. Ghuge, J. Kwon, V. Nagarajan and A. Sharma (Oper. Res. 70(2):786–804, 2022) reduce PCL assortment optimization to the maximum directed cut problem and obtain constant-factor approximation algorithms under cardinality, knapsack and matroid constraints.
Timeline.
- Ageev, Hassin and Sviridenko (2001): a 0.5-approximation for max-dicut with a cardinality constraint, by pipage rounding of the half-integral LP relaxation.
- Zhang, Rusmevichientong and Topaloglu (2020): an LP-based fixed-point framework for PCL assortment optimization, with a 0.25-approximation under a cardinality constraint.
- Ghuge, Kwon, Nagarajan and Sharma (2022): the reduction to max-dicut, and a 0.5-approximation for the cardinality-constrained problem (Theorem 6), the result formalized here.
Setting
There are products . Product has revenue and preference weight , and is the weight of not buying. Each ordered pair , , is a nest with dissimilarity parameter . An assortment is (in Lean, a Finset (Fin n)). With and , the expected revenue is
and is the optimal revenue under the cardinality constraint .
The ratio is linearized by : if and only if . The reduction graph has a vertex for each product and a dummy vertex , edges and , and explicit weights , affine in , such that is the weight of the edges leaving (pclWeight, pclEdges). keeps the edges of nonnegative weight.
For a digraph with weights , sizes and capacity , the LP (9) is
and LP is this LP on with ; its value is an upper bound on . Algorithm 1 takes an optimal basic solution of (9), possibly shifts it along the coordinates equal to and , applies pipage rounding (repeatedly moving mass between two fractional coordinates in the direction that does not decrease ), and outputs the better of the rounded set and that set plus the last fractional vertex.
Formalization targets
Goal: Theorem 6
For every finite loopless digraph with , every and unit sizes, and every optimal basic solution of LP (9) with value , Algorithm 1 has an output, and every output satisfies
For every PCL instance and , has a fixed point ; Algorithm 1 run on has an output; and for every output the assortment satisfies
Milestones
- Reduction: Lemma 1 ( is a cut function), Lemma 2 (sign structure, corrected), Lemmas 5 and 6 (the integer programs IP, IP and max-dicut on , ).
- Framework: Theorem 2 (), Lemma 7, Observation 1, Lemmas 8 and 9 (the fixed point of through the LP (D)).
- Rounding: Theorem 4 (pipage rounding with a knapsack constraint), Theorem 5 (basic solutions of (9) take values in ), Lemmas 10–12 (the shift and the two cases of the 0.5 bound).
Significance
The theorem doubles the best previous factor (0.25) for cardinality-constrained PCL assortment optimization and carries the factor 0.5 of Ageev et al. (2001) for cardinality-constrained max-dicut over to assortment optimization (Table 1 of the paper). The reduction of Lemma 1 is reusable: any LP-based approximation for constrained max-dicut with nonnegative weights becomes an approximation for PCL assortment optimization through the framework of Theorem 2. Theorem 5 extends the half-integrality of Ageev et al. from cardinality to general knapsack constraints, and is also the basis of the paper's knapsack results (Theorem 7 and the -approximation of §5).
The proofs of most of these statements are in the online appendices of the paper. None of the results has a machine-checked proof; this mission produces the first, together with Lean statements of the PCL model, the max-dicut LP and pipage rounding.
Difficulty
The ratio objective is neither submodular nor linear, and the reduction graph has negative edge weights, so max-dicut algorithms do not apply directly; Lemmas 2, 6 and 7 are what allow deleting the negative edges. The main obstacle is Theorem 5: the knapsack constraint destroys the half-integrality of the max-dicut polytope, and the proof must show that a single extra value (and ) suffices for every extreme point. Pipage rounding alone does not give the 0.5 bound at such points; the shift of equation (10) and the case analysis of Lemma 11 are needed, and the bound is tight there.
Formalization scope
Products are Fin n; the vertices of are Option (Fin n) with none the dummy. Real powers are Real.rpow, and is part of PCLInstance, so that empty nests contribute (Lean's ). Optimal values over finite families are suprema over finite nonempty sets; is an sSup over a nonempty polytope. "Basic feasible solution" is Set.extremePoints of the feasible set of (9) in the joint space, with fixed to off . Algorithm 1 and pipage rounding are relations: every nondeterministic choice is quantified, the goal asserts that an output exists and that every output satisfies the bound, and running time is not formalized. For the PCL part, Algorithm 1 runs on with all vertices of size , and the output is intersected with .
The statement "some assortment with at most products earns at least " is trivially true (take the optimum) and is not the theorem; the goal is about the outputs of Algorithm 1. Lemma 2 is corrected: its first clause holds with "nonpositive", not "negative".
A complete development needs LP duality and the existence of optimal extreme points for bounded polyhedra, the extreme-point analysis of Theorem 5, and the termination of pipage rounding. The max-dicut LP and pipage rounding definitions do not depend on the PCL model and can be reused for other rounding results. Proofs of any milestone, and alternative routes to the goal, are welcome.
Selected references
- R. Ghuge, J. Kwon, V. Nagarajan, A. Sharma, Constrained Assortment Optimization Under the Paired Combinatorial Logit Model, Operations Research 70(2):786–804, 2022. https://doi.org/10.1287/opre.2021.2188
- H. Zhang, P. Rusmevichientong, H. Topaloglu, Assortment Optimization Under the Paired Combinatorial Logit Model, Operations Research 68(3):741–761, 2020. https://doi.org/10.1287/opre.2019.1930
- A. Ageev, R. Hassin, M. Sviridenko, A 0.5-Approximation Algorithm for MAX DICUT with Given Sizes of Parts, SIAM Journal on Discrete Mathematics 14(2):246–255, 2001. https://doi.org/10.1137/S089548010036813X
- A. Ageev, M. Sviridenko, Approximation Algorithms for Maximum Coverage and Max Cut with Given Sizes of Parts, IPCO 1999, LNCS 1610:17–30. https://doi.org/10.1007/3-540-48777-8_2