Assortment Optimization under the Paired Combinatorial Logit Model

Size: px
Start display at page:

Download "Assortment Optimization under the Paired Combinatorial Logit Model"

Transcription

1 Assortment Optmzaton under the Pared Combnatoral Logt Model Heng Zhang, Paat Rusmevchentong Marshall School of Busness, Unversty of Southern Calforna, Los Angeles, CA Huseyn Topaloglu School of Operatons Research and Informaton Engneerng, Cornell Tech, New York, NY June 19, 017 We consder uncapactated and capactated assortment problems under the pared combnatoral logt model, where the goal s to fnd a set of products to maxmze the expected revenue obtaned from each customer. In the uncapactated settng, we can offer any set of products, whereas n the capactated settng, there s a lmt on the number of products that we can offer. We establsh that even the uncapactated assortment problem s strongly NP-hard. To develop an approxmaton framework for our assortment problems, we transform the assortment problem nto an equvalent problem of fndng the fxed pont of a functon, but computng the value of ths functon at any pont requres solvng a nonlnear nteger program. Usng a sutable lnear programmng relaxaton of the nonlnear nteger program and randomzed roundng, we obtan a 0.6-approxmaton algorthm for the uncapactated assortment problem. Usng randomzed roundng on a semdefnte programmng relaxaton, we obtan an mproved, but a more complcated, 0.79-approxmaton algorthm. Fnally, usng teratve varable fxng and coupled randomzed roundng, we obtan a 0.5-approxmaton algorthm for the capactated assortment problem. Our computatonal experments demonstrate that our approxmaton algorthms, on average, yeld expected revenues that are wthn 3.6% of a tractable upper bound on the optmal expected revenues. Key words : Customer choce modelng, par combnatoral logt model, assortment optmzaton. 1. Introducton Tradtonal revenue management models commonly assume that each customer arrves nto the system wth the ntenton to purchase a partcular product. If ths product s avalable for purchase, then the customer purchases t; otherwse, the customer leaves the system wthout a purchase. In realty, however, customers observe the set of avalable alternatves and make a choce among the avalable alternatves. Under such a customer choce process, the demand for a partcular product depends on the avalablty of other products. In ths case, dscrete choce models provde a useful representaton of demand snce dscrete choce models capture the demand for each product as a functon of the entre set of products n the offer set. A growng body of lterature ndcates that 1

2 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model capturng the choce process of customers usng dscrete choce models can sgnfcantly mprove the qualty of operatonal decsons; see, for example, Tallur and van Ryzn (004), Gallego et al. (004) and Vulcano et al. (010). Whle more sophstcated choce models yeld a more accurate representaton of the customer choce process, the assortment and other operatonal problems under more sophstcated choce models become more challengng. So, t s useful to dentfy sophstcated choce models, where the correspondng operatonal problems reman tractable. In ths paper, we study assortment problems under the pared combnatoral logt (PCL) model. In our problem settng, there s a fxed revenue assocated wth each product. Customers choose among the offered products accordng to the PCL model. The goal s to fnd an offer set that maxmzes the expected revenue obtaned from each customer. We consder both the uncapactated verson, where we can offer any subset of products, as well as the capactated verson, where there s a lmt on the number of products that we can offer. We show that even the uncapactated assortment problem s strongly NP-hard. We gve a general framework for constructng approxmaton algorthms for the assortment problem. We use ths framework to develop approxmaton algorthms for the uncapactated and capactated assortment problems. Our computatonal results demonstrate that our approxmaton algorthms perform qute well, yeldng solutons wth no larger than 3.6% optmalty gaps on average. The PCL model s a rather sophstcated choce model, allowng to capture ntrcate customer choce behavor. Ths model s known to be compatble wth the random utlty maxmzaton prncple, where each customer assocates random utltes wth the alternatves. The utltes are sampled from a certan dstrbuton. The customer knows the utltes and chooses the alternatve that provdes the largest utlty. Other commonly used choce models, such as the multnomal logt and nested logt models, are also compatble wth the random utlty maxmzaton prncple. An mportant feature of the PCL model s that t allows for a general correlaton structure among the utltes. In contrast, the multnomal logt model assumes that the utltes are ndependent. Under the nested logt model, the alternatves are grouped nto nests. The utltes of the alternatves n dfferent nests are ndependent and there s a sngle parameter that governs the correlaton between the utltes of the alternatves n the same nest. By allowng a general correlaton structure among the utltes, the PCL model captures the stuaton where the preference of a customer for one product offers nsght nto ther nclnaton towards other products. Furthermore, the PCL model subsumes the multnomal logt and nested logt models as specal cases. There s work n the lterature showng that the PCL model can provde advantages when capturng the demand process; see, for example, Prashker and Bekhor (1998), Koppleman and Wen (000) and Chen et al. (003). Although the PCL model can

3 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model 3 provde advantages n capturng the demand process, there s lttle research on understandng the complexty of makng operatonal decsons under the PCL model and provdng effcent algorthms for makng such decsons. Our work n ths paper s drected towards fllng ths gap. Thus, we beleve that the approxmaton algorthms that we provde for assortment optmzaton under the PCL model wll enhance the practcal appeal of the PCL model, whose effectveness for capturng the demand process has already been demonstrated n the lterature. Man Contrbutons: We make four man contrbutons. Frst, we study the complexty of the assortment problem under the PCL model. Second, we gve a general framework that can be used to develop approxmaton algorthms under the PCL model. Thrd, we use ths framework to gve approxmaton algorthms for the uncapactated assortment problem under the PCL model. Fourth, we gve an approxmaton algorthm for the capactated problem, also by usng our approxmaton framework. Although we use the same general framework for the uncapactated and capactated problems, the detals of the approxmaton algorthms are qute dfferent. Consderng the computatonal complexty of the problem, we show that the uncapactated assortment problem under the PCL model s strongly NP-hard. Ths result s n contrast wth the assortment problem under the closely-related multnomal logt and nested logt models. In partcular, as dscussed shortly n our lterature revew, there exst polynomal-tme algorthms to solve even the capactated assortment problem under the multnomal logt and nested logt models. Motvated by the computatonal complexty result, we focus on approxmaton algorthms. We gve a general framework for constructng approxmaton algorthms for the assortment problem. In ths approxmaton framework, we show that the assortment problem s equvalent to fndng the fxed pont of a one-dmensonal functon f : R R, whose evaluaton at a certan pont requres solvng a nonlnear nteger program. To obtan an α-approxmaton algorthm, we follow three steps. Frst, we desgn an upper bound f R : R R to the functon f( ), possbly through lnear programmng (LP) or semdefnte programmng (SDP) relaxatons of the nonlnear nteger program assocated wth computng the value of f( ) at a certan pont. Second, we compute the fxed pont ẑ of f R ( ). Thrd, we develop an α-approxmaton algorthm for the nonlnear nteger program that computes the value of f( ) at ẑ. It turns out that an α-approxmate soluton to ths nonlnear nteger program s an α-approxmate soluton to the assortment problem. So, we compute the fxed pont ẑ of f R ( ) and buld an approxmate soluton for computng f( ) at ẑ. For the uncapactated assortment problem, we apply the framework above by constructng f R ( ) through an LP relaxaton of the nonlnear nteger program that computes f( ). We show that f we compute f R ( ) by solvng an LP, then we can compute the fxed pont of f R ( ) also by solvng an LP. After we compute the fxed pont ẑ of f R ( ), we get a soluton to the LP that computes

4 4 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model f R ( ) at ẑ. We use randomzed roundng on ths soluton to get a 0.6-approxmate soluton to the nonlnear nteger program that computes f( ) at ẑ. Lastly, to obtan a determnstc algorthm, we de-randomze ths approach usng the standard method of condtonal expectatons; see Wllamson and Shmoys (011). Our framework can allow for a broad class of approxmaton algorthms. For example, f we construct f R ( ) by usng an SDP relaxaton of the nonlnear nteger program, then we can use the sphercal roundng method of Goemans and Wllamson (1995) to obtan a 0.79-approxmaton algorthm for the uncapactated assortment problem. Ths approxmaton algorthm requres solvng an SDP. Although we can theoretcally solve an SDP n polynomal tme, solvng large-scale nstances n practce can be computatonally dffcult. Therefore, the SDP relaxaton s arguably less appealng than the LP relaxaton from a practcal perspectve. We apply our framework to the capactated assortment problem as well. Here, we explot the structural propertes of the extreme ponts of the LP relaxaton and use an teratve varable fxng method, along wth coupled randomzed roundng, to develop a 0.5-approxmaton algorthm. In ths algorthm, f there are n products that can be offered to the customers, then we solve at most n successve LP relaxatons, fxng the value of one decson varable after solvng each LP relaxaton. Once we solve these LP relaxatons, we perform randomzed roundng on the bass of the soluton of the last LP relaxaton to obtan a soluton to the assortment problem. Usng the method of condtonal expectatons, we can de-randomze ths soluton to obtan a determnstc algorthm wth the same performance guarantee. Our computatonal experments demonstrate that our approxmaton algorthms for both the uncapactated and capactated versons perform qute well. Over a large collecton of test problems, the average optmalty gap of the solutons provded by these algorthms s no larger than 3.6%. Lterature Revew: Our broader obectve n ths paper s to enhance the practcal appeal of the PCL model n operatons management. As noted earler, there s lttle work on makng operatonal decsons under the PCL model. The only example we are aware of s L and Webster (015), where the authors consder prcng problems under the PCL model. In ther problem, the goal s to choose the prces for the products to maxmze the expected revenue obtaned from each customer. The authors gve suffcent condtons for the prce senstvtes of the products to ensure that the prcng problem can be solved effcently. We are not aware of other work that studes assortment or prcng problems under the PCL model. Despte lmted work on solvng operatonal problems under the PCL model, there s consderable work, especally n the transportaton lterature, on usng the PCL model to capture route choces. Koppleman and Wen (000) and Chen et al. (003) demonstrate that the general correlaton structure of the utltes under the PCL model s useful n modelng route choces, because dfferent routes overlap wth each other to varyng

5 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model 5 extents through shared street segments, creatng complex correlatons between the utltes provded by dfferent routes. Prashker and Bekhor (1998) demonstrate that the PCL model can provde mprovements over the multnomal logt model n predctng route choces. Chen et al. (014) and Karoonsoontawong and Ln (015) study varous traffc equlbrum problems under the PCL model and dscuss the benefts from usng ths choce model. The work n route choce lterature clearly demonstrates that the PCL model provdes mprovements over the multnomal logt and nested logt models n predctng user choces, especally when the utltes provded by dfferent alternatves exhbt complex correlaton structures. There s sgnfcant work on assortment and prcng problems under the multnomal logt and nested logt models. In the multnomal logt model, the utltes of the products are ndependent of each other. In the nested logt model, the products are grouped nto dsont nests. Assocated wth each nest, there s a dssmlarty parameter characterzng the correlaton between the utltes of the products n the same nest, but the utltes of products n dfferent nests are ndependent of each other. In the PCL model, there exsts one nest for each par of products, so the nests are overlappng. Assocated wth each nest, there s a dssmlarty parameter characterzng the correlaton between the utltes of each par of products. Therefore, when compared wth the multnomal logt and nested logt models, we can use the PCL model to specfy a sgnfcantly more general correlaton structure between the utltes of the products. As dscussed earler n ths secton, the general correlaton structure under the PCL model allows better capturng the demand n numerous applcatons comng from the route choce doman. Tran (00) provdes a thorough dscusson of the multnomal logt, nested logt and PCL models. Tallur and van Ryzn (004) and Gallego et al. (004) study assortment problems under the multnomal logt model. The authors show that t s optmal to offer a nested-by-revenue assortment, whch ncludes a certan number of products wth the largest revenues. In ths case, the optmal assortment can be found effcently by computng the expected revenue of each nested-by-revenue assortment. Rusmevchentong et al. (010) gve an effcent algorthm for the capactated assortment problem under the multnomal logt model. Davs et al. (013) show that the assortment problem under the multnomal logt model wth varous constrants can be formulated as an LP. Hanson and Martn (1996), Song and Xue (007), and Dong et al. (009) study prcng problems under the multnomal logt model. Assumng that the prce senstvtes of the products are the same, the authors show that the expected revenue s concave when expressed as a functon of the purchase probabltes of the products. Zhang et al. (016) generalze ths result by showng that the expected revenue s concave n the purchase probabltes under any generalzed extreme value model, whch s a general class of choce models ncludng the multnomal logt,

6 6 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model nested logt and PCL models. Zhang and Lu (013) show that even f the products have dfferent prce senstvtes, the prcng problem under the multnomal logt model can be formulated as a convex program, but ths result does not generalze to other choce models. Davs et al. (014) develop an effcent algorthm for the assortment problem under the nested logt model. Ther algorthm requres solvng an LP whose sze grows polynomally wth the number products and nests. Gallego and Topaloglu (014) study the assortment problem under the nested logt model when a capacty constrant lmts the number of products offered n each nest, whereas Feldman and Topaloglu (015) explore the same problem when a capacty constrant lmts the total number of products offered n all nests. L and Huh (011) nvestgate the prcng problem under the nested logt model when the products n the same nest share the same prce senstvtes and show that the problem can be formulated as a convex program. Gallego and Wang (014) analyze the same problem when the products have dfferent prce senstvtes. They show that the prcng problem can be cast as a search over a sngle dmenson. Rayfeld et al. (015) consder the prcng problem under the nested logt model wth arbtrary prce senstvtes and gve an algorthm to compute prces wth a performance guarantee. L and Huh (015) and L et al. (015) study prcng and assortment problems under the mult-level nested logt model, where the products are herarchcally organzed nto nests and subnests. Vulcano et al. (010) and Da et al. (014) use the multnomal logt and nested logt models n arlne applcatons. Organzaton: In Secton, we formulate both the uncapactated and capactated assortment problems and show that even the uncapactated assortment problem s strongly NP-hard. In Secton 3, we show that we can transform the assortment problem nto the problem of fndng the fxed pont of a functon f : R R. Computng the value of f( ) at any pont requres solvng a nonlnear nteger program. Thus, we approxmate f( ) by another functon f R : R R obtaned through an LP relaxaton of the nonlnear nteger program. We show how to fnd the fxed pont ẑ of f R ( ). We gve a suffcent condton for a certan subset of products to provde a performance guarantee for our assortment problem. The suffcent condton states that f the value of a sutable set functon computed at a subset of products exceeds α f R (ẑ), then the subset of products n queston provdes an α-approxmate soluton for the assortment problem. In Secton 4, we consder the uncapactated assortment problem and show how to obtan a subset of products satsfyng the suffcent condton wth α = 0.6. Also, we dscuss how we can obtan a 0.79-approxmaton algorthm by buldng on an SDP relaxaton. In Secton 5, we consder the capactated assortment problem and show how to obtan a subset of products satsfyng the suffcent condton wth α = 0.5. In Secton 6, we gve our computatonal experments. In Secton 7, we conclude.

7 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model 7. Problem Formulaton and Complexty The set of products s ndexed by N = {1,..., n}. The revenue of product s r 0. We use the vector x = (x 1,..., x n ) {0, 1} n to capture the subset of products that we offer to the customers, where x = 1 f and only f we offer product. We refer to the vector x smply as the assortment or the subset of products that we offer. Throughout the paper, we denote the vectors and matrces n bold font. We denote the collecton of nests by M = {(, ) N : }. For each nest (, ) M, we let γ [0, 1] be the dssmlarty parameter of the nest. For each product, we let v the preference weght of product. Under the PCL model, we can vew the choce process of a customer as takng place n two stages. Frst, the customer ether decdes to make a purchase n one of the nests or leaves the system wthout makng any purchase. In partcular, lettng V (x) = v 1/γ x + v 1/γ x and usng v 0 > 0 to denote the preference weght of the no-purchase opton, f we offer the subset of products x, then a customer decdes to make a purchase n nest (, ) wth probablty V (x) γ /(v 0 + (k,l) M V kl(x) γ kl ). Second, f the customer decdes to make a purchase n nest (, ), then she chooses product wth probablty v 1/γ product wth probablty v 1/γ be x / V (x), whereas she chooses x / V (x). Note that f a customer decdes to make a purchase n nest (, ), then she must choose one of the products n ths nest. If we offer the subset of products x and a customer has already decded to make a purchase n nest (, ), then the expected revenue that we obtan from the customer s R (x) = r v 1/γ x + r v 1/γ V (x) We use π(x) to denote the expected revenue that we obtan from a customer when we offer the subset of products x. In ths case, we have π(x) = V (x) γ v 0 + V R (k,l) M kl(x) γ (x) = kl x. V (x) γ R (x) v 0 + V. (x) γ Throughout the paper, we consder both uncapactated and capactated assortment problems. In the uncapactated assortment problem, we can offer any subset of products to the customers. In the capactated assortment problem, on the other hand, we have an upper bound on the cardnalty of the subset of products that we can offer to the customers. To capture both the uncapactated and capactated assortment problems succnctly, for some c Z +, we use F = {x {0, 1} n : N x c} to denote the feasble subsets of products that we can offer to the customers. Snce there are n products, the constrant N x c s not bndng when we have c n. Therefore, we mmedately obtan the uncapactated assortment problem by choosng a value of c that s no smaller n, whereas we obtan the capactated assortment problem wth other values

8 8 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model of c. In the assortment problem, our goal s to fnd a feasble subset of products to offer to the customers that maxmzes the expected revenue obtaned from each customer, correspondng to the combnatoral optmzaton problem z = max x F π(x) = max x F { V (x) γ R (x) v 0 + V (x) γ }. (Assortment) In the next theorem, we show that even the uncapactated verson of the Assortment problem s strongly NP-hard, ndcatng that t s unlkely that there exsts an effcent algorthm to solve the Assortment problem. We gve the proof of ths result n Appendx A. Theorem.1 (Computatonal Complexty) The Assortment problem s strongly NP-hard, even when we have F = {0, 1} n so that the problem s uncapactated. The proof of Theorem.1 uses a reducton from the max-cut problem, whch s a well-known NP-hard problem; see Garey and Johnson (1979). Motvated by ths complexty result, we focus on developng approxmaton algorthms for the Assortment problem. For some α [0, 1], an α-approxmaton algorthm s a polynomal-tme algorthm that, for every problem nstance, computes an assortment ˆx F, whose expected revenue s at least α tmes the optmal expected revenue; that s, notng that the optmal expected revenue s z, π(ˆx) α z. Closng ths secton, we note that our formulaton of the PCL model s slghtly dfferent from the one that often appears n the lterature. In the exstng lterature, the collecton of nests s often {(, ) N : < }, whereas n our formulaton, the collecton of nests s {(, ) N : }. If we let γ = γ for all (, ) N wth > and double the preference weght of the no-purchase opton n our formulaton, then t s easy to check that the two formulatons of the PCL model are consstent. Our formulaton of the PCL model wll sgnfcantly reduce the notatonal burden. 3. A General Framework for Approxmaton Algorthms In ths secton, we provde a general framework that s useful for developng approxmaton algorthms for the Assortment problem. Our framework s applcable to both the uncapactated and capactated assortment problems smultaneously Connecton to a Fxed Pont Problem The startng pont for our approxmaton framework s to relate the Assortment problem to the problem of computng the fxed pont of a one-dmensonal functon. Computng the value of ths functon at any pont requres solvng a nonlnear nteger program, but t turns out we can obtan

9 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model 9 an approxmate soluton to the Assortment problem by usng a tractable approxmaton to the one-dmensonal functon. We defne the functon f : R R as { } f(z) = max x F V (x) γ (R (x) z). (Functon Evaluaton) Snce x = 0 R n s a feasble soluton to the Functon Evaluaton problem on the rght sde above and provdes an obectve value of zero, we have f(z) 0 for all z R. Furthermore, f(z) s decreasng and v 0 z s strctly ncreasng n z. Thus, there exsts a unque ẑ 0 satsfyng f(ẑ) = v 0 ẑ. Note that the value of ẑ that satsfes f(ẑ) = v 0 ẑ s the fxed pont of the functon f( )/v 0. In the next lemma, we show that ths value of ẑ corresponds to the optmal obectve value of the Assortment problem and we can use ths value of ẑ n the Functon Evaluaton problem to obtan an optmal soluton to the Assortment problem. We do not gve a proof for Lemma 3.1 snce ths lemma follows as a corollary to a more general result that we shortly gve n Theorem 3.. Lemma 3.1 (Optmal Assortment from Fxed Pont) Let ẑ 0 satsfy f(ẑ) = v 0 ẑ and ˆx be an optmal soluton to the Functon Evaluaton problem when we solve ths problem wth z = ẑ. Then, we have π(ˆx) = ẑ = z ; so, ẑ s the optmal obectve value of the Assortment problem and ˆx s an optmal soluton to the Assortment problem. Thus, f we fnd the fxed pont ẑ of f( )/v 0 and solve the Functon Evaluaton problem wth z = ẑ, then an optmal soluton to the Functon Evaluaton problem s an optmal soluton to the Assortment problem. Fndng the fxed pont of f( )/v 0 s dffcult snce the Functon Evaluaton problem s a nonlnear nteger program. Instead, assume that we have a tractable upper bound f R ( ) on f( ) so that f R (z) f(z) for all z R. We construct f R ( ) n such a way that f R (z) s decreasng n z and f R (z) 0 for all z R, so there exsts a unque ẑ 0 satsfyng f R (ẑ) = v 0 ẑ. In the next theorem, we show that ths value of ẑ upper bounds the optmal obectve value of the Assortment problem and we can use ths value of ẑ to obtan an approxmate soluton. Theorem 3. (Approxmaton Framework) Assume that f R ( ) satsfes f R (z) f(z) for all z R. Let ẑ 0 satsfy f R (ẑ) = v 0 ẑ and ˆx F be such that V (ˆx) γ (R (ˆx) ẑ) α f R (ẑ) (Suffcent Condton) for some α [0, 1]. Then, we have π(ˆx) α ẑ α z ; so, ẑ upper bounds the optmal obectve value of the Assortment problem and ˆx s an α-approxmate soluton to the Assortment problem.

10 10 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model Proof: Notng that v 0 ẑ = f R (ẑ), we have α v 0 ẑ = α f R (ẑ) V (ˆx) γ (R (ˆx) ẑ) V (ˆx) γ (R (ˆx) α ẑ), where the frst nequalty uses the Suffcent Condton. Thus, we have α v 0 ẑ V (ˆx) γ (R (ˆx) α ẑ), n whch case, solvng for ẑ n the last nequalty, we get α ẑ V (ˆx) γ R (ˆx)/(v 0 + V (ˆx) γ ) = π(ˆx). Next, we show that ẑ z. We let x be an optmal soluton to the Assortment problem. Snce x s a feasble but not necessarly an optmal soluton to the Functon Evaluaton problem when ths problem s solved wth z = ẑ, we have f(ẑ) V (x ) γ (R (x ) ẑ). Notng that v 0 ẑ = f R (ẑ) f(ẑ), we obtan v 0 ẑ V (x ) γ (R (x ) ẑ). Solvng for ẑ n ths nequalty, we get ẑ V (x ) γ R (x )/(v 0 + V (x ) γ ) = π(x ) = z. Lemma 3.1 follows as a corollary to Theorem 3.. In partcular, as n Lemma 3.1, assume that ẑ 0 satsfes f(ẑ) = v 0 ẑ and ˆx s an optmal soluton to the Functon Evaluaton problem when we solve ths problem wth z = ẑ. In ths case, we have V (ˆx) γ (R (ˆx) ẑ) = f(ẑ), whch mples that ẑ and ˆx satsfes the Suffcent Condton wth f R ( ) = f( ) and α = 1. Therefore, by Theorem 3., we obtan π(ˆx) ẑ z. Snce ˆx s a feasble but not necessarly an optmal soluton to the Assortment problem, we also have z π(ˆx), whch yelds z π(ˆx) ẑ z. Thus, all of the nequaltes n the last chan of nequaltes hold as equaltes, n whch case, Lemma 3.1 follows. By Theorem 3., to obtan an α-approxmate soluton to the Assortment problem, we can execute the followng three steps. A General Approxmaton Framework Step 1: Construct an upper bound f R ( ) on f( ) such that f R (z) f(z) for all z R. Step : Fnd the fxed pont ẑ of f R ( )/v 0 ; that s, fnd the value of ẑ such that f R (ẑ) = v 0 ẑ. Step 3: Fnd an assortment ˆx F such that V (ˆx) γ (R (ˆx) ẑ) α f R (ẑ). In Secton 3., we show how to construct an upper bound f R ( ) on f( ) by usng an LP relaxaton of the Functon Evaluaton problem. In Secton 3.3, we show how to compute the fxed pont of f R ( )/v 0 by solvng an LP. Wth these results, we can execute Steps 1 and n our approxmaton framework. In Secton 4, we use randomzed roundng on an LP relaxaton of the Functon Evaluaton problem to construct an assortment ˆx that satsfes the Suffcent Condton wth α = 0.6, yeldng a 0.6-approxmaton algorthm. Also, we show that we can use an SDP relaxaton to satsfy the Suffcent Condton wth α = In Secton 5, we tackle the capactated assortment problem and use an teratve varable fxng method to construct an assortment that satsfes the Suffcent Condton wth α = 0.5, yeldng a 0.5-approxmaton algorthm.

11 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model Constructng an Upper Bound We construct an upper bound f R ( ) of f( ) by usng an LP relaxaton of the Functon Evaluaton problem. Notng the defnton of V (x) and R (x), we have V (x) γ (R (x) z) = (v 1/γ x + v 1/γ x ) (r γ z) v 1/γ v 1/γ x + (r z) v 1/γ x x + v 1/γ x. We let ρ (z) be the expresson on the rght sde above when x = 1 and x = 1 and θ (z) be the expresson on the rght sde above when x = 1 and x = 0. In other words, we have ρ (z) = (v 1/γ + v 1/γ ) (r γ z) v 1/γ v 1/γ + (r z) v 1/γ + v 1/γ and θ (z) = v (r z). There are only four possble values of (x, x ). In ths case, lettng µ (z) = ρ (z) θ (z) θ (z) for notatonal brevty, we can express V (x) γ (R (x) z) succnctly as V (x) γ (R (x) z) = ρ (z) x x + θ (z) x (1 x ) + θ (1 x ) x = µ (z) x x + θ x + θ x. Wrtng ts obectve functon as (µ (z) x x + θ x + θ x ), the Functon Evaluaton problem s equvalent to f(z) = max (µ (z) x x + θ (z) x + θ (z) x ) : x c, N x {0, 1} N. To construct an upper bound f R ( ) on f( ), we use a standard approach to lnearze the term x x n the obectve functon above. Defne the decson varable y {0, 1} wth the nterpretaton that y = x x. To ensure that y takes the value x x, we mpose the constrants y x + x 1, y x and y x. If x = 0 or x = 0, then the constrants y x and y x ensure that y = 0. If x = 1 and x = 1, then the constrant y x + x 1 ensures that y = 1. We defne our upper bound f R ( ) on f( ) by usng the LP relaxaton f R (z) = max (µ (z) y + θ (z) x + θ (z) x ) (Upper Bound) s.t. y x + x 1 (, ) M y x, y x (, ) M x c N 0 x 1 N, y 0 (, ) M. Snce the Upper Bound problem s an LP relaxaton of the Functon Evaluaton problem, we have f R (z) f(z) for all z R. Furthermore, settng x = 0 for all N and y = 0 for all (, ) M

12 1 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model provdes a feasble soluton to the Upper Bound problem. Therefore, we have f R (z) 0 for all z R. It s not mmedately clear that f R (z) s decreasng n z snce t s not mmedately clear that the obectve functon coeffcent µ (z) n the Upper Bound problem s decreasng n z. In the next lemma, we show that f R (z) s ndeed decreasng n z. In ths case, snce f R (z) s decreasng n z and f R (z) 0 for all z R, there exsts a unque ẑ 0 satsfyng f R (ẑ) = v 0 ẑ. Lemma 3.3 (Monotoncty of Upper Bound) The optmal obectve value of the Upper Bound problem s decreasng n z. Proof: Consder z + z and let (x, y ) wth y = {y : (, ) M} be an optmal soluton to the Upper Bound problem when we solve ths problem wth z = z +. Snce µ (z) = ρ (z) θ (z) θ (z) by the defnton of µ (z), we obtan f R (z + ) = (µ (z + ) y + θ (z + ) x + θ (z + ) x ) = (ρ (z + ) y + θ (z + ) (x y) + θ (z + ) (x y)) (ρ (z ) y + θ (z ) (x y ) + θ (z ) (x y )) = (µ (z ) y + θ (z ) x + θ (z ) x ) f R (z ), where the frst nequalty s by the fact that ρ (z) and θ (z) are decreasng n z, along wth the fact that y x and y x, whereas the second nequalty s by the fact that (x, y ) s a feasble but not necessarly an optmal soluton to the Upper Bound problem wth z = z. One useful property of the Upper Bound problem s that there exsts an optmal soluton to ths problem where the decson varable x takes a nonzero value only when r z and the decson varable y takes a nonzero value only when r z and r z. Thus, we need to keep the decson varable x only when r z and we need to keep the decson varable y only when r z and r z. Ths property allows us to sgnfcantly smplfy the Upper Bound problem. In partcular, let N(z) = { N : r z} and M(z) = {(, ) N(z) : }. In Lemma B.1 n Appendx B, we show that there exsts an optmal soluton x = {x : N} and y = {y : (, ) N} to the Upper Bound problem where we have x = 0 for all N(z) and y = 0 for all (, ) M(z). The proof of ths result follows by showng that f x = {x : N} and y = {y : (, ) N} s a feasble soluton to the Upper Bound problem, then we can set x = 0 for all N(z) and y = 0 for all (, ) M(z) whle makng sure that the soluton (x, y) remans feasble to the Upper Bound problem and we do not degrade the obectve value provded by ths soluton. In ths case, lettng 1( ) be the ndcator

13 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model 13 functon and droppng the decson varable x for all N(z) and the decson varable y for all (, ) M(z), we wrte the obectve functon of the Upper Bound problem equvalently as 1( N(z), N(z)) (µ (z) y + θ (z) x + θ (z) x ) + 1( N(z), N(z)) θ (z) x + 1( N(z), N(z)) θ (z) x. For the last two sums, we have 1( N(z), N(z)) θ (z) x = N \ N(z) θ N(z) (z) x and 1( N(z), N(z)) θ (z) x = N \ N(z) θ N(z) (z) x. Thus, the obectve functon of the Upper Bound problem takes the form (µ (z) (z) y + θ (z) x + θ (z) x ) + N \ N(z) θ N(z) (z) x. A smple lemma, gven as Lemma B. n Appendx B, shows that µ (z) 0 for all (, ) M(z). So, the decson varable y takes ts smallest possble value n the Upper Bound problem, whch mples that the constrants y x and y x are redundant. In ths case, the Upper Bound problem s equvalent to the problem f R (z) = max (µ (z) y + θ (z) x + θ (z) x ) + N \ N(z) θ (z) x (z) N(z) s.t. y x + x 1 (, ) M(z), (Compact Upper Bound) x c N(z) 0 x 1 N(z), y 0 (, ) M(z). We wll use the Upper Bound problem to fnd the fxed pont of f R ( )/v 0. We wll use the Compact Upper Bound problem above to fnd an assortment ˆx satsfyng the Suffcent Condton Computng the Fxed Pont To compute the fxed pont of f R ( )/v 0, we use the dual of the Upper Bound problem. For each (, ) M, let α, β, and γ be the dual varables of the constrants y x + x 1, y x and y x, respectvely. For each N, let λ be the dual varable of the constrant x 1. Let δ be the dual varable of the constrant x N c. The dual of the Upper Bound problem s f R (z) = mn c δ + N λ + α (Dual) s.t. α + β + γ µ (z) (, ) M δ + λ + 1( ) (α + α β γ ) (n 1) θ (z) N N α 0, β 0, γ 0 (, ) M, λ 0 N, δ 0. In the Dual problem above, the decson varables are α = {α : (, ) M}, β = {β : (, ) M}, γ = {γ : (, ) M}, λ = {λ : N} and δ. We treat z as fxed. In the next theorem, we show

14 14 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model that f we treat z as a decson varable and add one constrant to the Dual problem that nvolves the decson varable z, then we can recover the fxed pont of f R ( )/v 0. Theorem 3.4 (Fxed Pont Computaton by Usng an LP) Let ( ˆα, ˆβ, ˆγ, ˆλ, ˆδ, ẑ) be an optmal soluton to the problem mn c δ + N λ + α s.t. α + β + γ µ (z) (, ) M δ + λ + 1( ) (α + α β γ ) (n 1) θ (z) N N c δ + N λ + α = v 0 z α 0, β 0, γ 0 (, ) M, λ 0 N, δ 0, z s free. (Fxed Pont) Then, we have f R (ẑ) = v 0 ẑ; so, ẑ s the fxed pont of f R ( )/v 0. Proof: Let z be the fxed pont of f R /v 0 so that f R ( z) = v 0 z. We wll show that ẑ = z. Let (ᾱ, β, γ, λ, δ) be an optmal soluton to the Dual problem when we solve ths problem wth z = z. Thus, we have c δ + λ N + ᾱ = f R ( z) = v 0 z, whch mples that the soluton (ᾱ, β, γ, λ, δ, z) satsfes the last constrant n the Fxed Pont problem n the theorem. Furthermore, snce the soluton (ᾱ, β, γ, λ, δ) s feasble to the Dual problem, t satsfes the frst two constrants n the Fxed Pont problem as well. Thus, (ᾱ, β, γ, λ, δ, z) s a feasble but not necessarly an optmal soluton to the the Fxed Pont problem, whch mples that v 0 z = f R ( z) = c δ + λ + ᾱ c ˆδ + ˆλ + N N ˆα = v 0 ẑ, where the last equalty uses the fact that ( ˆα, ˆβ, ˆγ, ˆλ, ˆδ, ẑ) satsfes the last constrant n the Fxed Pont problem. The chan of nequaltes above mples that z ẑ. Next, we show that z ẑ. Note that ( ˆα, ˆβ, ˆγ, ˆλ, ˆδ) s a feasble soluton to the Dual problem wth z = ẑ, so that the obectve value of the Dual problem at the soluton ( ˆα, ˆβ, ˆγ, ˆλ, ˆδ) s no smaller than ts optmal obectve value. Therefore, we have f R (ẑ) c ˆδ + ˆλ N + ˆα = v 0 ẑ, where the equalty uses the fact that ( ˆα, ˆβ, ˆγ, ˆλ, ˆδ, ẑ) satsfes the last constrant n the Fxed Pont problem. Snce f R ( ) s a decreasng functon, havng f R ( z) = v 0 z and f R (ẑ) v 0 ẑ mples that z ẑ. Snce µ (z) and θ (z) are lnear functons of z, the Fxed Pont problem s an LP. Thus, we can compute the fxed pont of f R ( )/v 0 by solvng an LP.

15 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model Extreme Ponts The approxmaton algorthms that we gve for the uncapactated and capactated assortment problems use the propertes of the extreme ponts of the polyhedron defned by the set of feasble solutons to the Compact Upper Bound problem. For the capactated assortment problem, we also work wth the extreme ponts of ths polyhedron after fxng the values of some of the decson varables. As a functon of c, let P c denote the polyhedron defned by the set of feasble solutons to the Compact Upper Bound problem. Furthermore, for any subset H N(z), let P c (H) denote the polyhedron defned by the set of feasble solutons to the Compact Upper Bound problem, after fxng the values of the decson varables {x : H} at 1. That s, we have { P c = (x, y) [0, 1] N(z) R M(z) + : y x + x 1 (, ) M(z), { P c (H) = P c (x, y) [0, 1] N(z) R M(z) + : x = 1 }. H N(z) x c In the next lemma, we consder the uncapactated assortment problem. We show that all of the decson varables {x : N(z)} take a value n {0, 1, 1} n any extreme pont of P c(h). We gve the proof of ths result n Appendx C. Lemma 3.5 (Extreme Ponts n Uncapactated Problem) Assume that c > n. For any H N(z), let (ˆx, ŷ) be any extreme pont of P c (H). Then, we have ˆx { 0, 1, 1} for all N(z). } In the uncapactated assortment problem, the decson varables {x : N(z)} take one of the values n {0, 1, 1} n any extreme pont of P c(h), but the same result does not hold n the capactated assortment problem, as shown n the next example. Example 3.6 (No Half-Integralty n Capactated Problem) Consder the case where we have N(z) = 7, c = 3 and H =. Let ˆx = (,,,,,, 3 ) and ŷ = 0 R M(z) Note that we have ˆx = c and ˆx + ˆx 1 0 = ŷ for all (, ) M(z), whch mples that (ˆx, ŷ) P c (H). We clam that (ˆx, ŷ) s an extreme pont of P c (H). Assume on the contrary that there exst (ˆx + ɛ, ŷ + δ) P c (H) and (ˆx ɛ, ŷ δ) P c (H) wth (ɛ, δ) nonzero so that we have (ˆx, ŷ) = 1 (ˆx + ɛ, ŷ + δ) + 1 (ˆx ɛ, ŷ δ). Snce we have ŷ = 0, ŷ + δ R M(z) + and ŷ δ R M(z) +, t must be the case that δ = 0. Notng that (ˆx + ɛ, ŷ + δ) P c (H) and (ˆx ɛ, ŷ δ) P c (H), for each {1,..., 6}, the constrant y 7 x + x 7 1 yelds 0 = ŷ 7 + δ 7 ˆx + ɛ + ˆx 7 + ɛ 7 1 and 0 = ŷ 7 δ 7 ˆx ɛ + ˆx 7 ɛ 7 1. In the case, snce ˆx + ˆx 7 = 1 by the defnton of ˆx, the nequaltes above mply that ɛ 7 = ɛ for all {1,..., 6}. Also, the constrant N(z) x c yelds 7 =1 ˆx + 7 =1 ɛ c and

16 16 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model 7 ˆx =1 7 ɛ =1 c, n whch case, notng that 7 ˆx =1 = 3 = c by the defnton of ˆx, we obtan 7 ɛ =1 = 0. Combnng the last equalty wth the fact that ɛ = ɛ 7 for all {1,..., 6}, t follows that ɛ = 0 for all {1,..., 7}. So, (ɛ, δ) s the zero vector, whch s a contradcton. In the next lemma, we consder the capactated assortment problem. We show that f none of the decson varables {x : N(z)} take a fractonal value larger than 1 n an extreme pont of P c (H), then all of these varables take a value n {0, 1, 1}. The proof s n Appendx C. Lemma 3.7 (Extreme Ponts n Capactated Problem) Assume that c n. For any H N(z), let (ˆx, ŷ) be any extreme pont of P c (H). If there s no product N(z) such that 1 < ˆx < 1, then we have ˆx {0, 1, 1} for all N(z). We use the lemma above when developng an teratve varable fxng method for the capactated assortment problem. Note that Example 3.6 does not contradct Lemma 3.7, snce 1 < ˆx 7 < 1 n the extreme pont (ˆx, ŷ) consdered n ths example. 4. Applyng the Approxmaton Framework to the Uncapactated Problem In Sectons 3. and 3.3, we construct an upper bound f R ( ) on f( ) by usng an LP relaxaton of the Functon Evaluaton problem and fnd the fxed pont of f R ( )/v 0 by usng the Fxed Pont problem. Ths dscusson allows us to execute Steps 1 and n our approxmaton framework. In ths secton, we focus on Step 3 n our approxmaton framework. In partcular, consderng the uncapactated assortment problem, we gve a tractable approach to fnd a subset of products ˆx that satsfes V (ˆx) γ (R (ˆx) ẑ) 0.6 f R (ẑ). In ths case, Theorem 3. mples that ˆx s a 0.6-approxmate soluton to the uncapactated assortment problem. Let ẑ satsfy f R (ẑ) = v 0 ẑ. Snce the value of ẑ s fxed, to smplfy our notaton, we exclude the reference to ẑ n ths secton. In partcular, we let µ = µ (ẑ), θ = θ (ẑ), ρ = ρ (ẑ), f R = f R (ẑ), ˆN = N(ẑ) and ˆM = M(ẑ). Fnally, snce we consder the uncapactated assortment problem, we omt the constrant N(z) x c. Therefore, we wrte the Compact Upper Bound problem as f R = max s.t. (µ y + θ x + θ x ) + N \ ˆN θ x ˆN (,) ˆM y x + x 1 (, ) ˆM 0 x 1 ˆN, y 0 (, ) ˆM. Our goal s to fnd ˆx that satsfes V (ˆx) γ (R (ˆx) ẑ) 0.6 f R, where f R s the optmal obectve value of the problem above. We use randomzed roundng for ths purpose. Let (x, y ) be

17 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model 17 an optmal soluton to the problem above. We defne a random subset of products ˆX = { ˆX : N} by ndependently roundng each coordnate of x as follows. For each ˆN, we set { 1 wth probablty x ˆX = 0 wth probablty 1 x. For each N \ ˆN, we set ˆX = 0. Note that the subset of products ˆX s a random varable wth E{ ˆX } = x for all ˆN. In the next theorem, we gve the man result of ths secton. Theorem 4.1 (0.6-Approxmaton) Let approach. Then, we have { E ˆX be obtaned by usng the randomzed roundng } V ( ˆX) γ (R ( ˆX) ẑ) 0.6 f R. Proof: Here, we wll show that E{ V ( ˆX) γ (R ( ˆX) ẑ)} 0.5 f R and brefly dscuss how to refne the analyss to get the approxmaton guarantee of 0.6. The detals of the refned analyss are n Appendx D. As dscussed at the begnnng of Secton 3., we have V ( ˆX) γ (R ( ˆX) ẑ) = µ ˆX ˆX + θ ˆX + θ ˆX. So, snce { ˆX : N} are ndependent and E{ ˆX } = x, we have µ x x + θ x + θ x f ˆN, ˆN, E{V ( ˆX) γ θ x (R ( ˆX) f ẑ)} = ˆN, / ˆN θ x f / ˆN, ˆN 0 f / ˆN, / ˆN. Lettng [a] + = max{a, 0} and consderng the four cases above through the ndcator functon, we can wrte,) M E{V ( ˆX) γ (R ( ˆX) ẑ)} equvalently as E{V ( ˆX) γ (R ( ˆX) ẑ)} = = 1( ˆN, ˆN) (µ x x + θ x + θ x ) + 1( ˆN, ˆN) θ x + ˆM(µ x x + θ x + θ x ) + N \ ˆN (,) ˆN = 1( ˆN, ˆN) θ x θ x ˆM(µ [x + x 1] + + θ x + θ x ) + N \ ˆN (,) ˆN + µ (x x [x + x 1] + ) (,) ˆM = f R + (,) ˆM µ (x x [x + x 1] + ) f R (,) ˆM where the second equalty s by the fact that 1( ˆN, ˆN) θ x = N \ ˆN ˆN θ x = 1( ˆN, ˆN) θ x and the fourth equalty follows because µ 0 for all (, ) ˆM θ x µ,

18 18 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model by Lemma B. so that the decson varable y takes ts smallest possble value n an optmal soluton to the Compact Upper Bound problem, whch mples that y = [x + x 1] +. The fnal nequalty follows from the fact that 0 ab [a + b 1] + 1/4 for any a, b [0, 1]. We gve a lower bound on f R. Let ˆx = 1 for all ˆN and ŷ = 0 for all (, ) ˆM. In ths case, (ˆx, ŷ) s a feasble soluton to the LP that computes f R at the begnnng of ths secton, whch mples that the obectve value under (ˆx, ŷ) provdes a lower bound on f R. So, we lower bound f R as f R (,) ˆM θ + θ + N \ ˆN ˆN θ (,) ˆM θ + θ (,) ˆM θ + θ ρ = (,) ˆM where the second nequalty holds snce θ 0 for all ˆN, the thrd nequalty holds snce ρ 0 for all (, ) ˆM and the equalty follows from the defnton of µ. Usng the lower bound above on f R n the earler chan of nequaltes, we have E{V ( ˆX) γ (R ( ˆX) ẑ)} f R (,) ˆM µ 1 f R, whch s the desred result. Next, we brefly dscuss how to refne the analyss to mprove the approxmaton guarantee to 0.6. The refned analyss s lengthy and we defer the detals of the refned analyss to Appendx D. The dscusson above uses a lower bound on f R that s based on a feasble soluton (ˆx, ŷ) wth ˆx = 1 for all ˆN and ŷ = 0 for all (, ) ˆM. Ths lower bound may not be tght. The feasble regon n the LP that computes f R at the begnnng of ths secton s gven by the polyhedron P c wth c > n, where P c s as defned n Secton 3.4. By Lemma 3.5, f c > n, then any extreme pont (ˆx, ŷ) of P c satsfes ˆx {0, 1, 1} for all ˆN. In Appendx D, we enumerate over a relatvely small collecton of solutons (ˆx, ŷ) P c, where we have ˆx {0, 1, 1} for all ˆN and ŷ = [x +x 1] + for all (, ) ˆM. We pck the best one of these solutons to obtan a tghter lower bound on f R. In ths case, we can show that E{V ( ˆX) γ (R ( ˆX) ẑ)} 0.6 f R. The subset of products ˆX s a random varable, but n Theorem 3., we need a determnstc subset of products ˆx that satsfes the Suffcent Condton. Snce we manpulate fractons n the proof of Theorem 3., the subset ˆX does not necessarly delver the desred performance guarantee n expectaton. We use the standard method of condtonal expectatons to construct a determnstc subset of products ˆx by buldng on the subset ˆX; see Wllamson and Shmoys (011). Method of Condtonal Expectatons In the method of condtonal expectatons, we nductvely construct a subset of products x (k) = (ˆx 1,..., ˆx k, ˆX k+1,..., ˆX n ) for all k N, where the frst k products n ths subset are determnstc and the last n k products are random varables. These subsets of products are constructed n such a way that they satsfy E{V (x (k) ) γ (R (x (k) ) ẑ)} 0.6 f R for all k N. In ths case, the subset of products x (n) corresponds to a determnstc subset of products and t satsfes µ,

19 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model V (x (n) ) γ (R (x (n) ) ẑ) 0.6 f R, as desred. To nductvely construct the subsets of products x (k) = (ˆx 1,..., ˆx k, ˆX k+1,..., ˆX n ) for all k N, we start wth x (0) = ˆX. By Theorem 4.1, we have E{V (x (0) ) γ (R (x (0) ) ẑ)} 0.6 f R. Assumng that we have a subset of products x (k) that satsfes E{V (x (k) ) γ (R (x (k) ) ẑ)} 0.6 f R, we show how to construct a subset of products x (k+1) that satsfes E{V (x (k+1) ) γ (R (x (k+1) ) ẑ)} 0.6 f R. By the nducton assumpton, we have 0.6 f R E{V (x (k) ) γ (R (x (k) ) ẑ)}. Condtonng on ˆX k+1, we wrte the last nequalty as 0.6 f R P{ ˆX k+1 = 1} E{V (x (k) ) γ (R (x (k) ) ẑ) ˆX k+1 = 1} + P{ ˆX k+1 = 0} E{V (x (k) ) γ (R (x (k) ) ẑ) ˆX k+1 = 0}. We defne the two subsets of products as x (k) = (ˆx 1,..., ˆx k, 1, ˆX k+,..., ˆX n ) and as x (k) = (ˆx 1,..., ˆx k, 0, ˆX k+,..., ˆX n ). By the defnton of x (k), gven that ˆXk+1 = 1, we have x (k) = x (k). Gven that ˆX k+1 = 0, we have x (k) = x (k). So, we wrte the nequalty above as 0.6 f R P{ ˆX k+1 = 1} E{V ( x (k) ) γ (R ( x (k) ) ẑ)} + P{ ˆX k+1 = 0} E{V ( x (k) ) γ (R ( x (k) ) ẑ)} { max E{V ( x (k) ) γ (R ( x (k) ) ẑ)}, } E{V ( x (k) ) γ (R ( x (k) ) ẑ)}. Thus, ether E{V ( x (k) ) γ (R ( x (k) ) ẑ)} or E{V ( x (k) ) γ (R ( x (k) ) ẑ)} s at least 0.6 f R, ndcatng that we can use x (k) or x (k) as x (k+1). In both x (k) and x (k), the frst k + 1 products are determnstc and the last n k 1 products are random varables, as desred. We can also use an SDP relaxaton of the Functon Evaluaton problem to construct a subset of products ˆX that satsfes E{V ( ˆX) γ (R ( ˆX) ẑ)} 0.79 f R, but ths approach comes at the expense of solvng an SDP. We defer the detals to Appendx E. 5. Applyng the Approxmaton Framework to the Capactated Problem In ths secton, we consder Step 3 n our approxmaton framework for the capactated assortment problem. Lettng ẑ satsfy f R (ẑ) = v 0 ẑ, we focus on fndng a subset of products ˆx such that V (ˆx) γ (R (ˆx) ẑ) 0.5 f R (ẑ) and N ˆx c. In ths case, by Theorem 3., the subset of products ˆx s a 0.5-approxmate soluton. As n the prevous secton, snce the value of ẑ s fxed, to smplfy our notaton, we drop the argument ẑ n the Compact Upper Bound problem. In partcular, we let µ = µ (ẑ), θ = θ (ẑ), ρ = ρ (ẑ), f R = f R (ẑ), 19 ˆN = N(ẑ) and ˆM = M(ẑ). Furthermore, we observe that the set of feasble solutons n the Compact Upper Bound problem corresponds to the polyhedron P c gven n Secton 3.4. Therefore, notng that the optmal

20 0 Zhang, Rusmevchentong, and Topaloglu: Assortment Optmzaton under the PCL Model obectve value of the Compact Upper Bound problem s f R, we can equvalently wrte the Compact Upper Bound problem as { f R = max (µ y + θ x + θ x ) + N \ ˆN θ x : (x, y) P c }. (,) ˆM ˆN As ndcated by Example 3.6, t seems to be dffcult to predct the values of the varables {x : ˆN} n an extreme pont (x, y) of P c. Thus, we do not see a smple way of drectly roundng an optmal soluton to the LP above. Instead, we teratvely solve a modfed verson of the LP above, where we fx some of the decson varables {x : ˆN} at 1. Our goal s to obtan a soluton, where the decson varables {x : ˆN} ultmately all take values n {0, 1, 1} and the obectve value provded by ths soluton for the LP above s not too far from the optmal obectve value f R. In partcular, we consder the followng teratve varable fxng algorthm. Iteratve Varable Fxng Step 1: Set the teraton counter to k = 1 and ntalze the set of fxed varables H k =. Step : Let (x k, y k ) be an optmal soluton to the problem { f k = max (µ y + θ x + θ x ) + N \ ˆN } θ x : (x, y) P c (H k ), (Varable Fxng) ˆN (,) ˆM where the polyhedron P c (H) for H ˆN s as defned n Secton 3.4. Step 3: If there exsts some k ˆN such that 1 < xk k < 1, then set H k+1 = H k { k }, ncrease k by one and go to Step. Otherwse, stop. By the defnton of P c (H), the decson varables {x : H k } are fxed at 1 n the Varable Fxng problem n Step of the teratve varable fxng algorthm. Wthout loss of generalty, we assume that the optmal soluton (x k, y k ) to the Varable Fxng problem n Step s an extreme pont of P c (H k ). In Step 3 of the teratve varable fxng algorthm, f x k ( 1, 1) for all ˆN, then we stop. Lemma 3.7 mples that f x k ( 1, 1) for all ˆN, then x k {0, 1, 1} for all ˆN. Therefore, f the teratve varable fxng algorthm stops at teraton k wth the soluton (x k, y k ), then we have x k {0, 1, 1} for all ˆN. We use (x, y ) to denote the soluton to the Varable Fxng problem at the last teraton of the teratve varable fxng algorthm. We observe that ncludng each product ˆN n a subset wth probablty x may not provde a feasble soluton. In partcular, snce we have x {0, 1, 1} and ˆN x c, there can be as many as c decson varables such that x > 0. Includng each product ˆN n a subset wth probablty x may yeld a subset wth as many as c products. Instead, we use the followng coupled randomzed roundng approach to obtan a subset of products that satsfes the capacty constrant.

Technical Note: Capacity Constraints Across Nests in Assortment Optimization Under the Nested Logit Model

Technical Note: Capacity Constraints Across Nests in Assortment Optimization Under the Nested Logit Model Techncal Note: Capacty Constrants Across Nests n Assortment Optmzaton Under the Nested Logt Model Jacob B. Feldman, Huseyn Topaloglu School of Operatons Research and Informaton Engneerng, Cornell Unversty,

More information

Capacity Constraints Across Nests in Assortment Optimization Under the Nested Logit Model

Capacity Constraints Across Nests in Assortment Optimization Under the Nested Logit Model Capacty Constrants Across Nests n Assortment Optmzaton Under the Nested Logt Model Jacob B. Feldman School of Operatons Research and Informaton Engneerng, Cornell Unversty, Ithaca, New York 14853, USA

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

Pricing Problems under the Nested Logit Model with a Quality Consistency Constraint

Pricing Problems under the Nested Logit Model with a Quality Consistency Constraint Prcng Problems under the Nested Logt Model wth a Qualty Consstency Constrant James M. Davs, Huseyn Topaloglu, Davd P. Wllamson 1 Aprl 28, 2015 Abstract We consder prcng problems when customers choose among

More information

Approximation Methods for Pricing Problems under the Nested Logit Model with Price Bounds

Approximation Methods for Pricing Problems under the Nested Logit Model with Price Bounds Approxmaton Methods for Prcng Problems under the Nested Logt Model wth Prce Bounds W. Zachary Rayfeld School of Operatons Research and Informaton Engneerng, Cornell Unversty, Ithaca, New York 14853, USA

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

COS 521: Advanced Algorithms Game Theory and Linear Programming

COS 521: Advanced Algorithms Game Theory and Linear Programming COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Technical Note: An Expectation-Maximization Algorithm to Estimate the Parameters of the Markov Chain Choice Model

Technical Note: An Expectation-Maximization Algorithm to Estimate the Parameters of the Markov Chain Choice Model Techncal Note: An Expectaton-Maxmzaton Algorthm to Estmate the Parameters of the Markov Chan Choce Model A. Serdar Şmşek 1, Huseyn Topaloglu 2 August 1, 2017 Abstract We develop an expectaton-maxmzaton

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Technical Note: A Simple Greedy Algorithm for Assortment Optimization in the Two-Level Nested Logit Model

Technical Note: A Simple Greedy Algorithm for Assortment Optimization in the Two-Level Nested Logit Model Techncal Note: A Smple Greedy Algorthm for Assortment Optmzaton n the Two-Level Nested Logt Model Guang L and Paat Rusmevchentong {guangl, rusmevc}@usc.edu September 12, 2012 Abstract We consder the assortment

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

CS : Algorithms and Uncertainty Lecture 17 Date: October 26, 2016

CS : Algorithms and Uncertainty Lecture 17 Date: October 26, 2016 CS 29-128: Algorthms and Uncertanty Lecture 17 Date: October 26, 2016 Instructor: Nkhl Bansal Scrbe: Mchael Denns 1 Introducton In ths lecture we wll be lookng nto the secretary problem, and an nterestng

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Feature Selection: Part 1

Feature Selection: Part 1 CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?

More information

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed

More information

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

More information

Lecture 2: Gram-Schmidt Vectors and the LLL Algorithm

Lecture 2: Gram-Schmidt Vectors and the LLL Algorithm NYU, Fall 2016 Lattces Mn Course Lecture 2: Gram-Schmdt Vectors and the LLL Algorthm Lecturer: Noah Stephens-Davdowtz 2.1 The Shortest Vector Problem In our last lecture, we consdered short solutons to

More information

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem. prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of

More information

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem Appled Mathematcal Scences Vol 5 0 no 65 3 33 Interactve B-Level Mult-Objectve Integer Non-lnear Programmng Problem O E Emam Department of Informaton Systems aculty of Computer Scence and nformaton Helwan

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

More information

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Solutions to exam in SF1811 Optimization, Jan 14, 2015 Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable

More information

Lecture Space-Bounded Derandomization

Lecture Space-Bounded Derandomization Notes on Complexty Theory Last updated: October, 2008 Jonathan Katz Lecture Space-Bounded Derandomzaton 1 Space-Bounded Derandomzaton We now dscuss derandomzaton of space-bounded algorthms. Here non-trval

More information

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique Outlne and Readng Dynamc Programmng The General Technque ( 5.3.2) -1 Knapsac Problem ( 5.3.3) Matrx Chan-Product ( 5.3.1) Dynamc Programmng verson 1.4 1 Dynamc Programmng verson 1.4 2 Dynamc Programmng

More information

Lecture 14: Bandits with Budget Constraints

Lecture 14: Bandits with Budget Constraints IEOR 8100-001: Learnng and Optmzaton for Sequental Decson Makng 03/07/16 Lecture 14: andts wth udget Constrants Instructor: Shpra Agrawal Scrbed by: Zhpeng Lu 1 Problem defnton In the regular Mult-armed

More information

Lecture 4. Instructor: Haipeng Luo

Lecture 4. Instructor: Haipeng Luo Lecture 4 Instructor: Hapeng Luo In the followng lectures, we focus on the expert problem and study more adaptve algorthms. Although Hedge s proven to be worst-case optmal, one may wonder how well t would

More information

Appendix B. Criterion of Riemann-Stieltjes Integrability

Appendix B. Criterion of Riemann-Stieltjes Integrability Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for

More information

Amiri s Supply Chain Model. System Engineering b Department of Mathematics and Statistics c Odette School of Business

Amiri s Supply Chain Model. System Engineering b Department of Mathematics and Statistics c Odette School of Business Amr s Supply Chan Model by S. Ashtab a,, R.J. Caron b E. Selvarajah c a Department of Industral Manufacturng System Engneerng b Department of Mathematcs Statstcs c Odette School of Busness Unversty of

More information

Communication Complexity 16:198: February Lecture 4. x ij y ij

Communication Complexity 16:198: February Lecture 4. x ij y ij Communcaton Complexty 16:198:671 09 February 2010 Lecture 4 Lecturer: Troy Lee Scrbe: Rajat Mttal 1 Homework problem : Trbes We wll solve the thrd queston n the homework. The goal s to show that the nondetermnstc

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Finding Dense Subgraphs in G(n, 1/2)

Finding Dense Subgraphs in G(n, 1/2) Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals Smultaneous Optmzaton of Berth Allocaton, Quay Crane Assgnment and Quay Crane Schedulng Problems n Contaner Termnals Necat Aras, Yavuz Türkoğulları, Z. Caner Taşkın, Kuban Altınel Abstract In ths work,

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Finding Primitive Roots Pseudo-Deterministically

Finding Primitive Roots Pseudo-Deterministically Electronc Colloquum on Computatonal Complexty, Report No 207 (205) Fndng Prmtve Roots Pseudo-Determnstcally Ofer Grossman December 22, 205 Abstract Pseudo-determnstc algorthms are randomzed search algorthms

More information

Min Cut, Fast Cut, Polynomial Identities

Min Cut, Fast Cut, Polynomial Identities Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.

More information

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES BÂRZĂ, Slvu Faculty of Mathematcs-Informatcs Spru Haret Unversty barza_slvu@yahoo.com Abstract Ths paper wants to contnue

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty

Additional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,

More information

n ). This is tight for all admissible values of t, k and n. k t + + n t

n ). This is tight for all admissible values of t, k and n. k t + + n t MAXIMIZING THE NUMBER OF NONNEGATIVE SUBSETS NOGA ALON, HAROUT AYDINIAN, AND HAO HUANG Abstract. Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what

More information

Global Optimization of Truss. Structure Design INFORMS J. N. Hooker. Tallys Yunes. Slide 1

Global Optimization of Truss. Structure Design INFORMS J. N. Hooker. Tallys Yunes. Slide 1 Slde 1 Global Optmzaton of Truss Structure Desgn J. N. Hooker Tallys Yunes INFORMS 2010 Truss Structure Desgn Select sze of each bar (possbly zero) to support the load whle mnmzng weght. Bar szes are dscrete.

More information

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

Economics 101. Lecture 4 - Equilibrium and Efficiency

Economics 101. Lecture 4 - Equilibrium and Efficiency Economcs 0 Lecture 4 - Equlbrum and Effcency Intro As dscussed n the prevous lecture, we wll now move from an envronment where we looed at consumers mang decsons n solaton to analyzng economes full of

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

Online Appendix. t=1 (p t w)q t. Then the first order condition shows that

Online Appendix. t=1 (p t w)q t. Then the first order condition shows that Artcle forthcomng to ; manuscrpt no (Please, provde the manuscrpt number!) 1 Onlne Appendx Appendx E: Proofs Proof of Proposton 1 Frst we derve the equlbrum when the manufacturer does not vertcally ntegrate

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

Pricing and Resource Allocation Game Theoretic Models

Pricing and Resource Allocation Game Theoretic Models Prcng and Resource Allocaton Game Theoretc Models Zhy Huang Changbn Lu Q Zhang Computer and Informaton Scence December 8, 2009 Z. Huang, C. Lu, and Q. Zhang (CIS) Game Theoretc Models December 8, 2009

More information

NP-Completeness : Proofs

NP-Completeness : Proofs NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem

More information

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN FINITELY-GENERTED MODULES OVER PRINCIPL IDEL DOMIN EMMNUEL KOWLSKI Throughout ths note, s a prncpal deal doman. We recall the classfcaton theorem: Theorem 1. Let M be a fntely-generated -module. (1) There

More information

Lecture 4: Constant Time SVD Approximation

Lecture 4: Constant Time SVD Approximation Spectral Algorthms and Representatons eb. 17, Mar. 3 and 8, 005 Lecture 4: Constant Tme SVD Approxmaton Lecturer: Santosh Vempala Scrbe: Jangzhuo Chen Ths topc conssts of three lectures 0/17, 03/03, 03/08),

More information

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014) 0-80: Advanced Optmzaton and Randomzed Methods Lecture : Convex functons (Jan 5, 04) Lecturer: Suvrt Sra Addr: Carnege Mellon Unversty, Sprng 04 Scrbes: Avnava Dubey, Ahmed Hefny Dsclamer: These notes

More information

The Second Anti-Mathima on Game Theory

The Second Anti-Mathima on Game Theory The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player

More information

Limited Dependent Variables

Limited Dependent Variables Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages

More information

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution.

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution. Solutons HW #2 Dual of general LP. Fnd the dual functon of the LP mnmze subject to c T x Gx h Ax = b. Gve the dual problem, and make the mplct equalty constrants explct. Soluton. 1. The Lagrangan s L(x,

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Hidden Markov Models

Hidden Markov Models Hdden Markov Models Namrata Vaswan, Iowa State Unversty Aprl 24, 204 Hdden Markov Model Defntons and Examples Defntons:. A hdden Markov model (HMM) refers to a set of hdden states X 0, X,..., X t,...,

More information

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7 Stanford Unversty CS54: Computatonal Complexty Notes 7 Luca Trevsan January 9, 014 Notes for Lecture 7 1 Approxmate Countng wt an N oracle We complete te proof of te followng result: Teorem 1 For every

More information

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

More information

Case A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k.

Case A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k. THE CELLULAR METHOD In ths lecture, we ntroduce the cellular method as an approach to ncdence geometry theorems lke the Szemeréd-Trotter theorem. The method was ntroduced n the paper Combnatoral complexty

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

CSCE 790S Background Results

CSCE 790S Background Results CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each

More information

Estimation: Part 2. Chapter GREG estimation

Estimation: Part 2. Chapter GREG estimation Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

Online Appendix: Reciprocity with Many Goods

Online Appendix: Reciprocity with Many Goods T D T A : O A Kyle Bagwell Stanford Unversty and NBER Robert W. Stager Dartmouth College and NBER March 2016 Abstract Ths onlne Appendx extends to a many-good settng the man features of recprocty emphaszed

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Welfare Properties of General Equilibrium. What can be said about optimality properties of resource allocation implied by general equilibrium?

Welfare Properties of General Equilibrium. What can be said about optimality properties of resource allocation implied by general equilibrium? APPLIED WELFARE ECONOMICS AND POLICY ANALYSIS Welfare Propertes of General Equlbrum What can be sad about optmalty propertes of resource allocaton mpled by general equlbrum? Any crteron used to compare

More information

Complete subgraphs in multipartite graphs

Complete subgraphs in multipartite graphs Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G

More information

Hila Etzion. Min-Seok Pang

Hila Etzion. Min-Seok Pang RESERCH RTICLE COPLEENTRY ONLINE SERVICES IN COPETITIVE RKETS: INTINING PROFITILITY IN THE PRESENCE OF NETWORK EFFECTS Hla Etzon Department of Technology and Operatons, Stephen. Ross School of usness,

More information

Dynamic Bid Prices in Revenue Management

Dynamic Bid Prices in Revenue Management OPERATIONS RESEARCH Vol. 55, No. 4, July August 2007, pp. 647 661 ssn 0030-364X essn 1526-5463 07 5504 0647 nforms do 10.1287/opre.1060.0368 2007 INFORMS Dynamc Bd Prces n Revenue Management Danel Adelman

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

More information

= z 20 z n. (k 20) + 4 z k = 4

= z 20 z n. (k 20) + 4 z k = 4 Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5

More information