【论文阅读】(2020)Knapsack polytopes: a survey(下)

news2024/11/21 0:38:29

文章目录

  • 六、Valid inequalities, separation and computations 有效的不等式,分离和计算
  • 七、Complete linear descriptions of particular knapsack polytopes 特定背包多形体的完整线性描述
    • 7.1 Extended formulations
    • 7.2 Complete linear descriptions 完整的线性描述
  • 八、Geometric properties of knapsack polytopes 背包多面体的几何性质
  • 九、Independence systems and matroids
    • 9.1 Independence systems
    • 9.2 Knapsacks and matroids
  • 十、Integer knapsacks
  • 十一、Mixed-integer knapsacks
    • 11.1 General mixed-integer knapsacks
    • 11.2 Generalized flow models
  • 十二、Variants of the knapsack problem
    • 12.1 Multidimensional knapsack problem
    • 12.2 Cardinality constrained knapsack problem
    • 12.3 Generalized upper bound constraints
    • 12.4 Precedence constrained knapsack problem
    • 12.5 Generalized assignment problem
  • 十三、Conclusions


论文来源:(2020)Knapsack polytopes: a survey
作者:Christopher Hojny 等人


六、Valid inequalities, separation and computations 有效的不等式,分离和计算

前几章集中讨论了背包问题的强不等式,即最小覆盖不等式的(提升)。在实践中,一个重要的方面是生成违规有效不等式所需的计算工作量。在这种情况下,经常需要在运行时和不平等的强度之间进行权衡。此外,如果将强不等式应用于分支-切割方法(例如,如果它们是密集的),那么强不等式(在它们是面定义的意义上)不一定是最好的。

从背包生成的切割平面对于基于分支和切割的算法的良好性能是非常重要的。除了混合整数舍入和Gomory混合整数切割外,背包切割是最重要的切割平面类别[如果关闭背包分离,CPLEX 12.5的性能会降低14%(受影响实例的35%),参见Achterberg和Wunderling(2013)]。

在本章中,我们提出了几个有效不等式的家族,并强调了与这些不等式的生成相关的几个计算方面,如果我们要求它们切断当前的LP松弛解 x ∗ x^* x 这种方法由克劳德等人(1983)首创,并由许多研究人员进一步发展,如下所述。

Cover inequalities 覆盖不等式:分支-切割代码中大多数生成的不等式的基础是(最小)覆盖不等式(CIs)。结果表明,覆盖不等式的分离问题是np难的,见Ferreira(1994)和Klabjan等人(1998)。然而,克劳德等人

(1983)已经注意到这个问题可以写成背包问题

在这里插入图片描述

因此,如果权重 a a a β β β 为整数,则可在 O ( n β ) O(n β) O(nβ) 伪多项式时间内求解。

在实践中,覆盖通常是启发式地发现的,例如,通过贪婪地以 ( 1 − x j ⋆ ) / a j \left(1-x_j^{\star}\right) / a_j (1xj)/aj , 具体参照 Crowder 。(1983)。

在实践中,保险不平等几乎总是被消除的。一种简单的可能性是考虑扩展覆盖不等式(ECI),见(3)。Gabrel和Minoux(2002)提供了ECI的精确分离算法。Kaparis和Letchford(2010)提出了一种 O ( n β ) O(n β) O(nβ) 时间的积分权重算法,并针对同一问题提出了几种启发式算法。

Lifted cover inequalities 提升的覆盖不等式:第五章中描述的提升过程产生简单提升的覆盖不等式(9)。Gu等人(1999a)表明,简单lci的分离问题是np难的。Hunsaker和Tovey(2005)在实践中调查了这些不平等的强度。他们表明,即使将所有简单的 LCIs 都添加到一个背包问题中,也存在需要在分支绑定树中使用指数级节点来解决问题的实例。请注意“这个结果不是由二元背包问题的np -硬度所暗示的,因为这些问题的覆盖不等式分离是np -硬的”(Hunsaker和Tovey 2005, p. 219)。

如5.2节所述,向上和向下提升可以产生可能更强的不等式(14)。Gu et al(1998)对Zemel(1989)的算法进行了改进,得到 O ( ∣ C ∣ n 3 ) \mathcal{O}\left(|C| n^3\right) O(Cn3) 时间内的 down-lifting 。

如果考虑按顺序提升覆盖不等式,则必须固定一系列提升变量。例如,Hoffman和Padberg(1991)建议首先用 x j ⋆ > 0 x_j^{\star}>0 xj>0 ,则 D = { j ∈ C : x j ⋆ = 1 } D=\left\{j \in C: x_j^{\star}=1\right\} D={jC:xj=1} ,最后,用 x j ⋆ = 0 x_j^{\star}=0 xj=0 。Gu等人(1998)提出了另一种启发式方法;另见Martin(1998)。请注意,可以使用不同的(最小)覆盖开始,并且在解除之前可能不会违反相应的覆盖不等式。

Gu等人(2000)讨论了同时举升的各个方面,特别是举升函数的选择。

Merged cover inequalities 合并覆盖不等式:假设 C 1 C_1 C1 C 2 C_2 C2 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 的两个覆盖。那么,不等式 x ( C 1 ) ≤ ∣ C 1 ∣ − 1 x\left(C_1\right) \leq\left|C_1\right|-1 x(C1)C11 x ( C 2 ) ≤ ∣ C 2 ∣ − 1 x\left(C_2\right) \leq\left|C_2\right|-1 x(C2)C21 对于 P P P 是有效的,但不一定定义P的面。在这种情况下,为了加强这两个不等式,Hickman和Easton (2015a)建议,如果 ∣ C 1 ∩ C 2 ∣ ≤ 1 \left|C_1 \cap C_2\right| \leq 1 C1C21 ,则将这两个不等式合并为一个新的不等式。如果 ∣ C 1 ∩ C 2 ∣ = 1 \left|C_1 \cap C_2\right|=1 C1C2=1 ,则设 p p p 为该交集中的唯一元素。否则,设 p ∈ C 1 p∈C_1 pC1 是任意的。

由于 p p p 在每一个可行解中 x p ≤ 1 x_p≤1 xp1 ,所以归并的思想是替换 x p x_p xp x ( C 1 ) ≤ ∣ C 1 ∣ − 1 x\left(C_1\right) \leq\left|C_1\right|-1 x(C1)C11 by 1 ∣ C 2 ∣ − 1 x ( C 2 ) ≤ 1 \frac{1}{\left|C_2\right|-1} x\left(C_2\right) \leq 1 C211x(C2)1 中的出现,从而得到不等式

在这里插入图片描述

当然,这个不等式不一定对p有效。因此,Hickman和Easton (2015a)给出了充分的标准来证明合并后的不等式是有效的。此外,他们为合并两个一般有效不等式提供一个泛化。Hickman和Easton (2015b)给出了数值实验,证明了合并覆盖不等式的有效性。

(1, k)-Configurations:一个集合 N ∪ { t } N \cup\{t\} N{t} N ⊊ [ n ] N \subsetneq[n] N[n] t ∈ [ n ] \ N t \in[n] \backslash N t[n]\N , k ∈ { 2 , … , ∣ N ∣ } k \in\{2, \ldots,|N|\} k{2,,N} 被称为 (1, k)-Configurations。 如果 a ( N ) ≤ β a(N) \leq \beta a(N)β Q ∪ { t } Q \cup\{t\} Q{t} 是每个 Q ⊆ N Q \subseteq N QN ∣ Q ∣ = k |Q|=k Q=k 的最小覆盖。Padberg(1980)证明了对于任意**(1, k)-Configurations** N∪{t}, (1, k)-Configurations不等式:

在这里插入图片描述

对P有效,并定义了背包多面体约束的一个面到 N ∪ { t } N∪\{t\} N{t}。此外,Gottlieb和Rao(1988)在两个不相交的**(1, k)-Configurations**上提供了定义P的facet的充要条件。

示例42:考虑 K = { x ∈ { 0 , 1 } 4 : 3 x 1 + 5 x 2 + 6 x 3 + 7 x 4 ≤ 14 } K=\left\{x \in\{0,1\}^4: 3 x_1+5 x_2+6 x_3+7 x_4 \leq 14\right\} K={x{0,1}4:3x1+5x2+6x3+7x414} 。然后一个 (1, 2)-Configurations n = { 1 , 2 , 3 } n =\{1,2,3\} n={1,2,3} T = 4 T = 4 T=4 给出,因为 { 4 } \{4\} {4} { 1 , 2 , 3 } \{1,2,3\} {1,2,3} 中的任意两项的并集形成了一个最小覆盖。

如果 k = ∣ N ∣ k = |N| k=N , 则一个 (1, k)-configuration是一个最小覆盖。由于最小覆盖不等式的分离是np难的,Ferreira等人(1996)推测(1,k)-构形不等式的分离问题也是np难的。出于这个原因,他们提出了启发式分离这些不等式。然而,这一猜想缺少一个证明,因此,分离复杂性是开放的。

Coefficient increased cover inequalities:给定一个(不一定是最小的)覆盖C,经典的提升过程通过增加不包含在C中的变量的0系数来强化覆盖不等式 x ( C ) ≤ ∣ C ∣ − 1 x(C) \leq|C|-1 x(C)C1。Dietrich和Escudero(1992)描述了强化不等式的另一种方法。它们还允许增加覆盖变量的系数;它们的强化过程在 O ( n log ⁡ ( n ) ) \mathcal{O}(n \log (n)) O(nlog(n)) 时间内运行。

Lifted pack inequalities:基于 Weismantel(1997)的想法, Atamtürk(2005)考虑了pack不等式。集合 P ′ ⊆ [ n ] P^{\prime} \subseteq[n] P[n] 是一个pack,如果 a ( P ′ ) ≤ β a\left(P^{\prime}\right) \leq \beta a(P)β 。然后是pack不等式

在这里插入图片描述

是有效的。它由 x j ≤ 1 x_j≤1 xj1 的上界支配。然而,解除包不平等(LPIs)并不一定是主导的。有关提升过程的详细信息,请参阅Atamtürk(2005)。

示例43:(继续示例5)考虑背包不等式的运行示例 4 x 1 + 5 x 2 + 6 x 3 + 7 x 4 + 9 x 5 ≤ 13 4 x_1+5 x_2+6 x_3+7 x_4+9 x_5 \leq 13 4x1+5x2+6x3+7x4+9x513 。一个包由 P ′ = { 1 , 5 } P^{\prime}=\{1,5\} P={1,5} ,不等式为 4 x 1 + 9 x 5 ≤ 13 4 x_1+9 x_5 \leq 13 4x1+9x513 。提高 x 2 x_2 x2 收益率

在这里插入图片描述
这产生了解除的包不等式 4 x 1 + 9 x 2 + 9 x 5 ≤ 13 4 x_1+9 x_2+9 x_5 \leq 13 4x1+9x2+9x513 ,它似乎不受任何解除的盖不等式的支配。

注意,将 x 3 x_3 x3 提升到pack { 1 , 2 } \{1,2\} {1,2} 中得到 4 x 1 + 5 x 2 + 4 x 3 ≤ 9 4 x_1+5 x_2+4 x_3 \leq 9 4x1+5x2+4x39 ,这是由最小覆盖不等式 x 1 + x 2 + x 3 ≤ 2 x_1+x_2+x_3 \leq 2 x1+x2+x32 所决定的。因此,并不是所有取消的pack不平等都同样有用。

Weight inequalities:Weismantel(1997)给出了LPIs的一个特殊例子。设 r = β − a ( P ′ ) r=\beta-a\left(P^{\prime}\right) r=βa(P) 为填料 P ′ P^{\prime} P 的剩余容量。然后将带 a j > r a_j>r aj>r 的系数向上提升,得到所谓的Weight inequalities(WIs):

在这里插入图片描述

Weismantel还提出了权重不等式的一种变体,通过减少包中一个变量的系数,并修改 [ n ] \ P ′ [n] \backslash P^{\prime} [n]\P 中的变量的系数。这类更大的不等式包含权重不等式,Weismantel证明了它可以在伪多项式时间内分离出来。他的算法也可以用来分离权值不等式本身。

对于积分权重,Kaparis和Letchford(2010)将运行时间改进为 O ( ( n + a ˉ ) β ) \mathcal{O}((n+\bar{a}) \beta) O((n+aˉ)β) ,其中 a ˉ = max ⁡ { a 1 , … , a n } \bar{a}=\max \left\{a_1, \ldots, a_n\right\} aˉ=max{a1,,an} ,并设计了一个有效的启发式。WIs 的分离问题的复杂性似乎是开放的,尽管Martin(1998)声称更一般的版本是np难的。

Atamtürk(2005)指出,使用超可加函数可以获得更强的LPIs。

示例44:在运行的示例中,使用pack{1,5}得到 r = 0 r = 0 r=0 。因此,它的重量不等式就是原来的背包不等式。但是,pack{1,2}的剩余容量 r = 13 − 9 = 4 r = 13−9 = 4 r=139=4,权值不等式为 4 x 1 + 5 x 2 + 2 x 3 + 3 x 4 + 5 x 5 ≤ 9 4 x_1+5 x_2+2 x_3+3 x_4+5 x_5 \leq 9 4x1+5x2+2x3+3x4+5x59

注意,这并不等于相应的解除包不等式。

进一步的不等式是Weismantel(1997)的扩展权不等式,它基于三个相互不相交的集合 T , I T, I T,I, and { k } \{k\} {k} 。为简洁起见,我们用 I = ∅ I=\varnothing I= 来呈现这个版本。在这种情况下, P ′ = T P^{\prime}=T P=T 是一个pack和 a ( P ′ ) + a k > β a\left(P^{\prime}\right)+a_k>\beta a(P)+ak>β 。不等式:

在这里插入图片描述

在这里插入图片描述

这样的不等式可以被解除以获得被解除的扩展权不等式(LEWIs)。对于给定的扩展重量不等式,Weismantel(1997)证明,只要背包数据是积分的,提升系数可以在多项式时间内计算出来,计算结果见Martin(1998)。

示例45:考虑背包不等式 2 x 1 + 2 x 2 + 2 x 3 + 4 x 4 ≤ 7 2 x_1+2 x_2+2 x_3+4 x_4 \leq 7 2x1+2x2+2x3+4x47 。带着 pack P ′ = { 1 , 2 , 3 } P^{\prime}=\{1,2,3\} P={1,2,3} r = 1 r=1 r=1 k = 4 k=4 k=4

在这里插入图片描述
扩展的权值不等式 x 1 + x 2 + x 3 + 2 x 4 ≤ 3 x_1 + x_2 + x_3 + 2x_4≤3 x1+x2+x3+2x43 ,这是facet定义的,但不是一个提升的覆盖不等式。

Gomory cuts:让一个 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 为n维背包不等式,令 J ⊆ [ n ] J⊆[n] J[n] 。Glover等人(1997)建议使用Gomory割

在这里插入图片描述

其中 λ ∈ R + J × { 0 } \lambda \in \mathbb{R}_{+}^{J \times\{0\}} λR+J×{0} ,作为切割平面。基于代理分析,他们表明这些切割平面可以在多项式时间内分离。此外,他们展示了如何加强这些切割面,他们认为这些加强不能通过经典的提升程序找到。

Park和Lee(2011)提出了一种在伪多项式时间内运行的背包多形体(和固定电荷背包多形体)分离Gomory切割的精确方法。

Lee(2012)给出了一种有效的启发式方法来分离这些不等式。

Exact separation:Boyd (1992,1993a, 1994)开发了一种算法,通过分离和优化的等价来精确分离背包多面体的不等式。

另一种使用极性概念的精确分离方法由Boccia(2006)描述,由Kaparis和Letchford(2010)陈述。Y an和Boyd(1998)也考虑了混合整数背包集。Kaparis和Letchford(2010)对0/1情况和Fukasawa和Goycoolea(2011)对混合整数情况进行了改进和扩展;更多信息请参阅Goycoolea(2006)和Fukasawa(2008)博士论文。Fukasawa和Goycoolea还考虑了由场景行产生的背包——注意,所有剩下的方法都只使用原始行。V asil 'ev (= V asilyev) (V asilyev 2009)引入了一种不同的实现,并应用于广义分配问题,另见Avella等人(2010)的广泛计算研究。Avella等人(2013)考虑将连续变量聚合为单个变量的MIPs。V asilyev等人(2016)已经改进了背包精确分离的实现。

Computations:Wolter(2006)与框架SCIP进行了对比,得出LEWIs ( I = ∅ I=\varnothing I=)对对偶界的影响最大。此外,同时提升与顺序LCIs无显著差异。

Kaparis和Letchford(2010)比较了不同的精确和启发式背包分离程序与精确分离程序。在稀疏MIPLIB实例的子集上,如果使用启发式,LCIs可以用很少的计算工作量缩小大量的差距。

此外,它们比ECIs或简单的LCIs表现更好。在这个测试集中,wi和LPIs的性能较差,ECIs和简单LCIs的性能非常相似。然而,在密集多维背包实例的测试集上,LPIs的性能明显优于WIs和LCIs。

Fukasawa和Goycoolea(2011)对这些实验进行了补充,比较了MIPLIB病例中mir -闭合的相对强度和背包切口的精确分离。他们观察到,通常有很大一部分差距是由MIR不等式弥补的,因此很难改进它们。

(Q12) 分离**(1, k)-Configurations**、提升包、重量和提升扩展重量不等式的复杂度是多少?

(Q13) 补充变量是否有助于加强类似于补充混合整数舍入的不等式,见Marchand和Wolsey (2001)?

(Q14) 给定一个任意不等式,识别它是否定义了背包多面体的一个方面的复杂性是多少,即,判断它是否有效且存在n个仿射独立点满足它的相等的复杂性是多少?

关于问题(Q14),复杂性类 D P D^P DP 是由Papadimitriou和Y annakakis(1984)提出的,包含了NP和coNP中语言交叉形成的语言。Papadimitriou和Wolfe(1988)提出了识别旅行推销员多面体的facet为 D p D^p Dp 完备的问题。因此,Hartvigsen和Zemel(1992)推测对背包面的识别也是 D p D^p Dp 完备的。


七、Complete linear descriptions of particular knapsack polytopes 特定背包多形体的完整线性描述

为了找到背包多形体的完整线性描述,我们提出了两种方法:在第7.1节中构造了扩展公式,在第7.2节中我们考虑了已知原始空间中完整线性描述的特殊情况。

7.1 Extended formulations

前面论述的很大一部分考虑了通过线性不等式系统对背包解的凸包的描述,即找到 P = { x ∈ R n : A x ≤ b } P=\left\{x \in \mathbb{R}^n: A x \leq b\right\} P={xRn:Axb} 的某个矩阵A和向量b。解决寻找这样一个系统及其可能的指数大小的困难的一种方法是考虑一个高维表示,它可以投射到P上,即 P = { x ∈ R n : ∃ y ∈ R d A x + B y ≤ b } P=\left\{x \in \mathbb{R}^n: \exists y \in \mathbb{R}^d A x+B y \leq b\right\} P={xRn:yRdAx+Byb}

这种高维多面体称为初始多面体的扩展形式。

不幸的是,引入额外变量的可能性通常不允许为背包多面体找到精确的多项式大小的公式,因为在任何扩展公式中都存在背包多面体 P ⊆ R n P \subseteq \mathbb{R}^n PRn 需要 Ω ( 2 n ) \Omega\left(2^{\sqrt{n}}\right) Ω(2n ) 多个不等式,参见Pokutta和Vyve(2013)。然而,Bienstock(2008)能够证明,对于每个 ε ∈ ( 0 , 1 ) ε∈(0,1) ε(0,1) ,存在一个ε-近似扩展公式 O ( ε − 1 n 2 + ⌈ ε − 1 ⌉ ) \mathcal{O}\left(\varepsilon^{-1} n^{2+\left\lceil\varepsilon^{-1}\right\rceil}\right) O(ε1n2+ε1)

如果 P ∩ { 0 , 1 } n = Q ∩ { 0 , 1 } n P \cap\{0,1\}^n=Q \cap\{0,1\}^n P{0,1}n=Q{0,1}n ,则公式 Q = { x ∈ R n : ∃ y ∈ R d A x + B y ≤ b } Q=\left\{x \in \mathbb{R}^n: \exists y \in \mathbb{R}^d A x+B y \leq b\right\} Q={xRn:yRdAx+Byb} ,则称为ε-近似扩展式

在这里插入图片描述

对所有 w ∈ R n w \in \mathbb{R}^n wRn 都成立。请注意,Bienstock的扩展公式对于任何固定的ε都是n的多项式大小,但在逆近似质量 1 ε \frac{1}{ε} ε1 中呈指数增长。此外,虽然ε-近似公式存在于扩展空间中,但Faenza和Sanità(2015)证明了在原始空间中一般不存在n大小多项式的ε-近似公式,显示了扩展公式的力量。

寻找背包多面体扩展公式的一种可能性是利用背包问题可以通过动态规划来解决。对于承认动态规划方案的一般问题,Martin等人(1990)推导了这些问题解集的扩展公式。在下面,我们提出了背包问题的扩展公式,其中我们遵循Conforti等人(2010)的介绍。

考虑0/1背包问题 max ⁡ { c ⊤ x : a ⊤ x ≤ β , x ∈ { 0 , 1 } n } \max \left\{c^{\top} x: a^{\top} x \leq \beta, \quad x \in\{0,1\}^n\right\} max{cx:axβ,x{0,1}n} , 对于积分权值 a ∈ Z + n a \in \mathbb{Z}_{+}^n aZ+n and β ∈ Z + \beta \in \mathbb{Z}_{+} βZ+ 。对于每个整数 0 ≤ β ′ ≤ β 0 \leq \beta^{\prime} \leq \beta 0ββ k ∈ [ n ] k \in[n] k[n] ,则定义值函数

在这里插入图片描述
显然,要解决背包问题,我们需要计算 f ( n , β ) f (n,β) f(n,β) 。这可以用Bellman方程来实现

在这里插入图片描述
伴随着 f ( 0 , β ′ ) = 0 f\left(0, \beta^{\prime}\right)=0 f(0,β)=0 f ( k , β ′ ) = − ∞ f\left(k, \beta^{\prime}\right)=-\infty f(k,β)= if β ′ < 0 \beta^{\prime}<0 β<0

该动态规划算法等价于在图 ( V , A 0 ∪ A ) \left(V, A^0 \cup A\right) (V,A0A) 上求一条从 v 0 , 0 v_{0,0} v0,0 v n , β v_{n,β} vn,β 的最小权值路径,节点为 v k , β ′ v_{k, \beta^{\prime}} vk,β for k ∈ [ n ] 0 , 0 ≤ β ′ ≤ β k \in[n]_0, 0 \leq \beta^{\prime} \leq \beta k[n]0,0ββ ,弧有两种类型。弧子集 A 0 A^0 A0 ( v k − 1 , β ′ , v k , β ′ ) \left(v_{k-1, \beta^{\prime},} v_{k, \beta^{\prime}}\right) (vk1,β,vk,β) ,权值为0, k ∈ [ n ] k \in[n] k[n] and β ′ ∈ [ β ] 0 \beta^{\prime} \in[\beta]_0 β[β]0 。弧子集A包含弧 ( v k − 1 , β ′ − a k , v k , β ′ ) \left(v_{k-1}, \beta^{\prime}-a_k, v_{k, \beta^{\prime}}\right) (vk1,βak,vk,β) 的权重 c k c_k ck , k ∈ [ n ] k∈[n] k[n] β ′ ∈ { a k , … , β } \beta^{\prime} \in\left\{a_k, \ldots, \beta\right\} β{ak,,β}

将Martin等人(1990)的结果应用于此图,得出以下定理,其中y变量模型是否在最小权值路径上使用弧。

定理46:设 a ∈ Z + n a \in \mathbb{Z}_{+}^n aZ+n and β ∈ Z + \beta \in \mathbb{Z}_{+} βZ+ 。给出了 P a , β P^{a, \beta} Pa,β of O ( n β ) \mathcal{O}(n \beta) O(nβ) 的推广公式

在这里插入图片描述
其中 δ − ( v ) \delta^{-}(v) δ(v) δ + ( v ) \delta^{+}(v) δ+(v) 分别是节点 v v v 的进弧和出弧的集合。

对于特定的背包多面体,即所谓的orbisacks(在第7.2节中有更详细的描述),Loos(2011)展示了如何使用基于动态规划图的扩展公式来计算原始空间中的完整线性描述。然而,一般来说,这个扩展的公式太复杂,无法计算到原始x变量的投影,或者太大,在实践中没有用处。

(Q15) 背包问题承认一个完全多项式的时间近似格式,但背包多形体是否也承认一个ε-近似的扩展公式,其大小在 n n n 1 ε \frac{1}{ε} ε1 之间多项式增长?

7.2 Complete linear descriptions 完整的线性描述

现在我们讨论一些特殊情况,在这些情况下,通过利用特定的组合结构可以得到完整的线性描述。

Weakly super-increasing knapsack polytopes:对于每一个 i ∈ [ n ] i∈[n] i[n] ,如果 a ( [ i − 1 ] ) ≤ a i a([i-1]) \leq a_i a([i1])ai ,例如,对于每一个i∈[n],如果 a i = 2 i a_i=2^i ai=2i ,则一个背包称为弱超递增。Laurent和Sassano(1992)研究了这种背包,他们证明弱超增长背包多形体完全由盒约束和 O ( n ) O(n) O(n) 最小覆盖不等式描述。这些封面可以显式地构造。

Sequential knapsack polytopes:一个背包被称为顺序背包,如果 a i a_i ai 对于每一个 i ∈ [ n − 1 ] i∈[n−1] i[n1] a i + 1 a_i+1 ai+1 的除数。除了 a i = 2 i a_i = 2^i ai=2i 外,序列背包的另一个例子是 a i = ∏ j = 1 i j a_i=\prod_{j=1}^i j ai=j=1ij 。Pochet和Weismantel(1998)提供了序列背包多面体的完整线性描述。描述中包含了一个不等式,对于子集 W ⊆ [ n ] W⊆[n] W[n] 的每一个组合,对于W的分区 B = B 1 ∪ ⋯ ∪ B m B=B_1 \cup \cdots \cup B_m B=B1Bm ,对于 [ m ] [m] [m] 的一个置换,因此是非常大的。该公式的分离复杂性是一个悬而未决的问题。

One coefficient:背包的不等式只有一个重量,即 a = ( λ , … , λ ) ⊤ a=(\lambda, \ldots, \lambda)^{\top} a=(λ,,λ) λ ∈ N \lambda \in \mathbb{N} λN ,等价于一个基数约束 1 ⊤ x ≤ ⌊ β / λ ⌋ \mathbb{1}^{\top} x \leq\lfloor\beta / \lambda\rfloor 1xβ/λ 完整的描述是由这个单一的不等式和平凡的不等式。

Two coefficients:背包的不等式有两个不同的系数,即存在 k ∈ [ n ] k∈[n] k[n] ,使得背包不等式由

在这里插入图片描述

对于 λ , μ \lambda, \mu λ,μ, 和 β ∈ N \beta \in \mathbb{N} βN ,由Weismantel(1996)研究;这包括一个系数的特殊情况,如果 λ = μ λ = μ λ=μ 。他提出了一个指数大小的完整线性描述,由八个不等式族组成。该公式的分离问题在多项式时间内可解。Dahl和Foldnes(2003)处理了 λ = 1 λ = 1 λ=1 的特殊情况,他们证明了在这种情况下,四个不等式族通过完全对偶完整给出了一个完整的线性描述。Hartmann(1994)表明,在这种情况下,背包多形体的完整线性描述可以在线性时间内分离出来。

Small and large coefficients:Bader等人(2018)考虑了背包不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ ,其系数相对于右边不是相对小就是相对大。

n , k ∈ N n, k \in \mathbb{N} n,kN ,且 3 ≤ k ≤ n 3≤k≤n 3kn , s ∈ { k , … , n } s∈\{k,…,n\} s{k,,n}。假设对于每一个 i ∈ S : = [ S ] i∈S:= [S] iS:=[S] ,系数满足 β k + 1 < a i ≤ β k \frac{\beta}{k+1}<a_i\leq\frac{\beta}{k} k+1β<aikβ ;对于每一个 i ∈ L : = { s + 1 , … , n } i \in L:=\{s+1, \ldots, n\} iL:={s+1,,n} ,系数满足 k − 1 k + 1 β < a i ≤ β \frac{k-1}{k+1} \beta<a_i \leq \beta k+1k1β<aiβ

然后Bader等人提出了以下由 a ⊤ x ≤ β a^{\top} x \leq \beta axβ (无证明):

在这里插入图片描述

其中 i ( R ) = min ⁡ R i(R)=\min R i(R)=minR R R R 中系数最小的一项的系数 a i ( R ) a_{i(R)} ai(R)

(1, k)-configurations:回顾第六章中(1,k)-configurations的定义。Padberg(1980)指出,如果 N = [ n − 1 ] N = [n−1] N=[n1] t = n t = n t=n 形成一个**(1, k)-configurations** w.r.t.,则背包不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ(1, k)-configurations 不等式的完全线性描述为:

在这里插入图片描述
对于每个子集 S ⊆ [ n − 1 ] S⊆[n−1] S[n1] k ≤ ∣ S ∣ k≤|S| kS 以及盒约束和 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 。这个结果也扩展到一般的包装或多维背包问题。

在系数差距Weismantel(1997)考虑背包给出

在这里插入图片描述
其中, N j N_j Nj j ∈ [ β ] j∈[β] j[β] 中权重为 j j j 的所有项的集合。他可以对这两种特殊情况推导出完全的线性描述

(a) N j = ∅ N_j=\varnothing Nj= for all 1 < j ≤ ⌊ β 2 ⌋ 1<j \leq\left\lfloor\frac{\beta}{2}\right\rfloor 1<j2β

(b) N j = ∅ N_j=\varnothing Nj= for all 1 < j ≤ ⌊ β 3 ⌋ 1<j \leq\left\lfloor\frac{\beta}{3}\right\rfloor 1<j3β and N j = ∅ N_j=\varnothing Nj= for all j ≥ ⌊ β 2 ⌋ + 1 j \geq\left\lfloor\frac{\beta}{2}\right\rfloor+1 j2β+1

这些特殊情况的例子如下

(a) a = ( 1 , 1 , 5 , 7 , 7 , 8 ) ⊤ , β = 9 a=(1,1,5,7,7,8)^{\top}, \beta=9 a=(1,1,5,7,7,8),β=9

(b) a = ( 1 , 1 , 1 , 4 , 4 , 5 ) ⊤ , β = 10 a=(1,1,1,4,4,5)^{\top}, \beta=10 a=(1,1,1,4,4,5),β=10

如果背包定义了一个阵面,背包多形体可以完全用Edmonds(1971)的结果来描述。在这种情况下,强覆盖不等式的所有扩展集足以给出一个公式,参见Wolsey(1975)。在第7节中讨论了背包定义阵面的情况。

Graphic knapsacks:假设假设11成立。Wolsey(1975)称背包图为存在 t ∈ [ n − 1 ] t∈[n−1] t[n1] ,使得其背包不等式满足在 a t + 1 + a t > β a_{t+1}+a_t>\beta at+1+at>β a ( [ t ] ) ≤ β a([t]) \leq \beta a([t])β 。同样,所有强覆盖不等式的集合提供了一个完整的线性描述,参见Wolsey(1975)。注意 t = 1 t = 1 t=1 对应于拟阵的情况。

Further cases:Gillmann和Kaibel(2006)引入了所谓的revlex-initial polytopes,即 conv ⁡ { x ∈ { 0 , 1 } n : x ≺ rlex ⁡ v } \operatorname{conv}\left\{x \in\{0,1\}^n: x \prec_{\operatorname{rlex}} v\right\} conv{x{0,1}n:xrlexv} 对于某些 v ∈ { 0 , 1 } n v \in\{0,1\}^n v{0,1}n ,. “≺rlex” 是严格的反向字典顺序。如果 a ˉ i = 2 i , i ∈ [ n ] \bar{a}_i=2^i, i \in[n] aˉi=2i,i[n] ,那么它对应于 conv ⁡ { x ∈ { 0 , 1 } n : a ˉ ⊤ x ≤ a ˉ ⊤ v − 1 } \operatorname{conv}\left\{x \in\{0,1\}^n:\right.\left.\bar{a}^{\top} x \leq \bar{a}^{\top} v-1\right\} conv{x{0,1}n:aˉxaˉv1},即背包多面体。Gillmann和Kaibel用 n n n 的多项式个不等式提供了一个完整的线性描述。

相关的对象是orbisacks,即 m × 2 m × 2 m×2 大小的0/1-矩阵的凸包,使得第一列在字典上不小于第二列。使用前面定义的向量 a ˉ \bar{a} aˉ , orbisacks可以被看作是特殊的背包多面体 conv ⁡ { ( x , y ) ∈ { 0 , 1 } m × 2 : a ˉ ⊤ ( y − x ) ≤ 0 } \operatorname{conv}\left\{(x, y) \in\{0,1\}^{m \times 2}:\right.\left.\bar{a}^{\top}(y-x) \leq 0\right\} conv{(x,y){0,1}m×2:aˉ(yx)0}(变量互补后)。Kaibel和Loos(2011)发现了一个完整的线性描述,Loos(2010)提出了一个优化orbisacks的算法。

完整的线性描述有 Θ ( 3 m ) \Theta\left(3^m\right) Θ(3m) 个面,并且有一个面定义不等式,使得其最大系数和最小系数之比为 2 m − 2 2^{m−2} 2m2

最后,Hojny和Pfetsch(2019)引入了所谓的symresacks conv ⁡ { x ∈ { 0 , 1 } n : x ⪰ lex ⁡ γ ( x ) } \operatorname{conv}\left\{x \in\{0,1\}^n:\right.\left.x \succeq_{\operatorname{lex}} \gamma(x)\right\} conv{x{0,1}n:xlexγ(x)} 用于某些坐标排列γ,其中 "‘ ⪰ \succeq lex " " " 指的是词典的顺序。这些多形体又是背包多形体(在补充变量之后),取决于 γ γ γ ,有时可以推导出完整的线性描述(Hojny 2018)。

示例47:(继续算例5)背包多面体 conv conv ⁡ { x ∈ { 0 , 1 } 5 : 4 x 1 + 5 x 2 + 6 x 3 + 7 x 4 + 9 x 5 ≤ 13 } \operatorname{conv}\left\{x \in\{0,1\}^5: 4 x_1+5 x_2+6 x_3+7 x_4+9 x_5 \leq 13\right\} conv{x{0,1}5:4x1+5x2+6x3+7x4+9x513} 的完整线性描述由 x i ≥ 0 , i ∈ [ 5 ] x_i \geq 0, i \in[5] xi0,i[5] ,强覆盖的四个扩展不等式(例21)以及提升覆盖不等式 x 2 + x 3 + x 4 + 2 x 5 ≤ 2 x_2+x_3+x_4+2 x_5 \leq 2 x2+x3+x4+2x52 (示例27)给出。上面讨论的任何情况都不包含这个结果。

(Q16) 是否还有其他种类的背包多形体可以提供完整的线性描述?

(Q17) 是否有一个共同的概括所有提到的完全线性描述?


八、Geometric properties of knapsack polytopes 背包多面体的几何性质

前几节讨论了背负式多面体的代数性质和描述。本节研究它们的几何性质。

Number of vertices 顶点数:由于具有积分权的背包问题可以通过动态规划来解决,因此P的顶点数V §原则上可以通过一个修正的目标为零的动态规划算法来计算。然而,这个过程的运行时间在 n n n , 而 V ( P ) V(P) V(P) 可能是指数级的。

为了找到 V ( P ) V (P) V(P) 的近似值,Dyer(2003)提出了一种随机算法来近似计算P的顶点数,即对于每一个 ε > 0 ε > 0 ε>0 ,该算法以高概率返回范围为 ( 1 − ε ) V ( P ) (1−ε)V (P) (1ε)V(P) ( 1 + ε ) V ( P ) (1 + ε)V (P) (1+ε)V(P) 的数字。Dyer的算法运行时间为 O ( n 5 / 2 log ⁡ ( n ε − 1 ) + ε − 2 n 2 ) \mathcal{O}\left(n^{5 / 2} \sqrt{\log \left(n \varepsilon^{-1}\right)}+\varepsilon^{-2} n^2\right) O(n5/2log(nε1) +ε2n2) ,基于动态规划。

对于一个没有上界的整数背包,即 K ′ = { x ∈ Z + n : a ⊤ x ≤ β } K^{\prime}=\left\{x \in \mathbb{Z}_{+}^n: a^{\top} x \leq \beta\right\} K={xZ+n:axβ} 时,Hayes和Larman(1982)表明, conv ⁡ ( K ′ ) \operatorname{conv}\left(K^{\prime}\right) conv(K) 的顶点数从上面以 log ⁡ 2 n ( 4 β a 1 ) \log _2^n\left(\frac{4 \beta}{a_1}\right) log2n(a14β) 为界。然而,对于二元情况,我们不知道这样的上界。

Adjacency:Geist和Rodin(1992)在基于次优解局部改进的优化算法的激励下,研究了P的顶点的邻接性,证明了P的两个不同顶点 u , v u, v u,v 是相邻的,当且仅当 P P P 的两个不同顶点 w 1 , w 2 w^1, w^2 w1,w2 不存在时, u + v 2 = w 1 + w 2 2 \frac{u+v}{2}=\frac{w^1+w^2}{2} 2u+v=2w1+w2 。他们还表明,决定两个顶点是否相邻是NP完全的。

Distribution of classes of facets:如前所述, P a , β P^{a, \beta} Pa,β 最重要的一类不等式是由平凡不等式 0 ≤ x ≤ 1 0≤x≤1 0x1 和简单提升最小覆盖不等式(LCI)给出的。正如我们在第4节中所看到的,简单的LCIs在通过分支和切割解决背包问题中发挥着重要作用。然而,目前尚不清楚有多少非平凡方面是由简单LCIs定义的。为了得到这个值的估计值,我们在小维中列举了所有全维背包多面体,计算了它的facet,并检查了哪些facet是由简单的LCIs定义的。在本节的剩余部分,我们将更详细地描述我们的调查。

在下面,我们考虑所有可能的0/1背包多面体,直到0/1等价,即直到坐标排列和互补变量。因此,如果我们在讨论中提到一个背包多面体,我们实际上是指它的0/1等价类的所有成员。

为了枚举维数 n ∈ [ 5 ] n∈[5] n[5] 的全维背包多形体的所有等价类,我们使用了Aichholzer(2000)的一个工具,该工具允许访问维数不超过5的0/1多形体的所有等价类,该工具可通过他的网页获得。由于web工具允许每次查询最多访问1000个多边形,我们使用了Aichholzer为我们提供的离线版本。

使用Aichholzer的工具,我们首先生成了所有(等价类)0/1多面体的列表,并删除了所有非全维多面体。为了从这个列表中过滤出背包多面体,我们检查每个多面体 P P P 是否存在线性不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ ,通过求解 L P LP LP P P P 的顶点与剩余的二进制点分离,cf. Bradley et al(1974)。注意,我们不能局限于不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 且系数为非负, a i ≤ β a_i≤β aiβ i ∈ [ n ] i∈[n] i[n] ,因为我们处理的是多面体的等价类。也就是说,我们不能保证假设1成立。此外,请注意,我们将 [ 0 , 1 ] n [0,1]^n [0,1]n 分类为背包多面体。

表2表明,当 n ≥ 4 n≥4 n4 时,绝大多数0/1多面体都不是背包。特别是,当 n = 5 n = 5 n=5 时,0/1多形体的1226617个等价类中只有92个由背包组成。对于较小的维度,所有( n ∈ [ 2 ] n∈[2] n[2] )或大约一半的多面体( n = 3 n = 3 n=3 )都是背包。

在这里插入图片描述

在这里插入图片描述

{ P 1 , … , P N } \{P1,…,PN\} {P1,,PN} 是全维背包多形体各等价类的代表集合。对于每个代表的 P i P_i Pi ,我们使用polymake (Gawrilow和Joswig 2000)计算 P i P_i Pi 的facet描述。通过遍历所有facet定义不等式 P i P_i Pi ,我们将切面分为两类:由盒约束定义的平凡切面和非平凡切面。对于每个非平凡的方面,我们检查它是否是强覆盖的扩展不等式(因此,是一个简单的解除覆盖不等式)。对于其余的非平凡方面,我们手工检查是否存在最小覆盖,以便对应的定义不等式的方面是一个顺序简单的LCI。

为了评估简单LCIs在定义Pi不等式的非平凡facet中的分布, i ∈ [ N ] i∈[N] i[N] ,令 f i f_i fi , t i t_i ti , 和 c i c_i ci 分别是所有facet的数量,平凡facet的数量,与简单LCIs相关的facet的数量,以及剩余facet的数量。分数 ρ i = c i f i − t i ρ_i = \frac{c_i}{f_i−t_i} ρi=fitici 描述了由简单lci定义的非平凡facet的部分。为了描述基于简单LCIs的非平凡面分布,我们使用了统计度量:

在这里插入图片描述
表3显示,在维度2和3中,每个定义不等式的非平凡方面都是一个简单的LCI。在维度4和5中,所有非平凡方面的大多数都是由简单lci定义的,因为在这两种情况下,算术平均值大于90.0%,中值甚至为100.0%。然而,也有一些等价类,其中只有相对较少的面是由简单LCIs定义的,例如,在维度5中存在一个背包多面体,这样所有定义面的不等式中只有大约54.5%是简单LCIs。

我们的实验表明,简单的LCIs对于描述小维度的背包多形体很重要。这些发现对一般二进制程序也有启示。

为了增强二进制程序,我们可以考虑由其约束矩阵的不等式定义的背包多形体 P P P ,并为 P P P 分离有效不等式。如果约束矩阵是稀疏的,表3表明分离简单LCIs对于为 P P P 找到强切割平面很重要,因此,二进制程序。这一假设也得到了Kaparis和Letchford(2010)在第4节提到的MIPLIB实例上的数值实验的支持。然而,如果约束矩阵变得更密集,简单的LCIs可能不足以给出背包多面体的紧密公式,参见表3中4维和5维的最小百分比值。

(Q18) 二元背包多面体的顶点数量 V ( P ) V(P) V(P) 是否存在一个强的上下界?

(Q19) 是否有一种计算 V ( P ) V(P) V(P) 的例程比枚举 P P P 的所有顶点更快?

(Q20) 是否存在一种(递归)方案来构造所有固定维的背包多面体?

(Q21) 背包式多面体中可以用简单 LCIs 描述的非平凡面的最小百分比是多少?当 n → ∞ n→∞ n 时,它趋于 0 0 0 吗?


九、Independence systems and matroids

9.1 Independence systems

一个独立系统是一个元组 ( F , I ) (F, \mathcal{I}) (F,I) , 其中 F F F 是一个有限集, I ⊆ 2 F \mathcal{I} \subseteq 2^F I2F ,使得(i) ∅ ∈ I \varnothing \in \mathcal{I} I , (ii)如果 Y ∈ I Y \in \mathcal{I} YI ,然后 Y ∈ I Y \in \mathcal{I} YI 对于每一个 Y ⊆ X Y⊆X YX , I \mathcal{I} I 中包含的集合称为独立的,而 2 F \ I 2^F \backslash \mathcal{I} 2F\I 中的集合称为相关的。令 a ∈ R + n , β ∈ R + a \in \mathbb{R}_{+}^n, \beta \in \mathbb{R}_{+} aR+n,βR+ ,考虑系统 I K : = { I ⊆ [ n ] : a ( I ) ≤ β } \mathcal{I}_K:=\{I \subseteq[n]: a(I) \leq \beta\} IK:={I[n]:a(I)β} ,包含背包 K = K a , β K=K^{a, \beta} K=Ka,β 的可行解的指标集。根据引理3, M : = ( [ n ] , I K ) \mathcal{M}:=\left([n], \mathcal{I}_K\right) M:=([n],IK) 是一个独立系统,这就引出了下面的问题:

(Q22) 哪些独立系统可以表示为背包问题?

虽然这个问题在一般情况下是开放的,但对于我们在下面概述的稳定集类,一个完整的描述是可用的。给定一个无向图 G = ( V , E ) G = (V, E) G=(V,E) G G G 中的稳定集是集合 S ⊆ V S⊆V SV ,使得 S S S 中的节点在 G G G 中成对不相邻。

I ( G ) \mathcal{I}(G) I(G) 表示 G G G 中所有稳定集的集合,即 I ( G ) = { S ⊆ V : S \mathcal{I}(G)=\{S \subseteq V: S I(G)={SV:S 是稳定的 } \} } 。显然, S ( G ) : = ( V , I ( G ) ) \mathcal{S}(G):=(V, \mathcal{I}(G)) S(G):=(V,I(G)) 是一个独立系统。我们需要以下概念。

定义48:一个图 G = ( V , E ) G = (V, E) G=(V,E) 是一个阈值图当且仅当下列等价性质中的一个(即所有)成立:

  • (a) 图 G G G 可以由一个单顶点图通过重复应用以下两个操作来构造:
    • (i) 将单个孤立顶点添加到图形中。
    • (ii) 将一个顶点添加到与所有其他顶点连接的图形中。
  • (b) 存在实数 S S S 和实顶点权值 w ( v ) w(v) w(v) ,使得对于每一对不同的顶点 u , v ∈ V u, v∈V u,vV { u , v } \{u, v\} {u,v} G G G 的边当且仅当 w ( u ) + w ( v ) > S w(u) + w(v) > S w(u)+w(v)>S
  • c) 存在一个实数 β β β 和实数顶点权重 a ( v ) a(v) a(v) 这样对于每一个 S ⊆ V S \subseteq V SV 我们有 ∑ v ∈ S a ( v ) ≤ β \sum_{v \in S} a(v) \leq \beta vSa(v)β 当且仅当 S S S 是稳定集。

特别地,在取诱导子图时,阈值图的性质是封闭的。因此,一旦已知了(a)的构造序列,就可以很容易地构造一个表征 S ( G ) S(G) S(G) 的背包不等式:假设 G G G 的节点按照它们在构造序列中添加的顺序被标记为 1 , … , ∣ V ∣ 1, \ldots,|V| 1,,V 。令 ∑ i = 1 j a i j x i ≤ β j \sum_{i=1}^j a_i^j x_i \leq \beta_j i=1jaijxiβj 是表征步骤 j j j 中构造的图的不等式,然后令 x 1 ≤ 1 x_1≤1 x11 是单顶点图的不等式。如果我们应用构造规则(i)来获得 ( j + 1 ) (j + 1) (j+1)- s t st st 图,我们定义 ( j + 1 ) (j + 1) (j+1)- s t st st 不等式为

在这里插入图片描述

如果我们应用构造规则(ii),我们定义不等式

在这里插入图片描述

定理49:(Chvátal and Hammer 1975)当且仅当 G G G 是阈值图时, S ( G ) S(G) S(G) 可以表示为一个背包问题。

示例50:图1显示了阈值图的示例。使用构造规则(i)和(ii),这个图可以按照以下顺序构造,这也产生了对应子图的背包特征。

在这里插入图片描述

在这里插入图片描述

图1所示图的稳定集与最终背包不等式的可行解相吻合。

将阈值的概念推广到一般独立系统似乎是复杂的。据我们所知,它只用于拟阵,参见下一小节的定理55。

9.2 Knapsacks and matroids

一般来说,求解独立系统上的线性优化问题是NP困难的。

对于特殊情况的拟阵,贪心算法可以非常有效地求解。回想一下,一个独立系统 ( F , I ) (F, \mathcal{I}) (F,I) 是一个阵,如果对于每个 S ⊆ F S⊆F SF , 包含在S中的最大独立集具有相同的大小。

由于拟阵的优化很容易,关于背包 K K K 的一个问题是如何刻画相关的独立系统 M = ( [ n ] , I K ) \mathcal{M}=\left([n], \mathcal{I}_K\right) M=([n],IK) , 其中 I K : = { I ⊆ [ n ] : a ( I ) ≤ β } \mathcal{I}_K:=\{I \subseteq[n]: a(I) \leq \beta\} IK:={I[n]:a(I)β} 是否是一个拟阵。

定义51:假设假设11成立。极大独立集 C ⊆ [ n ] C⊆[n] C[n] M M M 的上限,如果对于所有 i ∈ C i∈C iC ,且有 i + 1 ∈ [ n ] \ C i + 1∈[n]\backslash C i+1[n]\C ,则 ( C \ { i } ) ∪ { i + 1 } ∉ I K (C \backslash\{i\}) \cup\{i+1\} \notin \mathcal{I}_K (C\{i}){i+1}/IK

使用天花板的概念,Wolsey(1975)提供了一个背包的完整特征,可以通过一个矩阵来表示。此外,由于Cerdeira和Barcia(1995),存在通过强覆盖的特征。

定理52:(Cerdeira和Barcia 1995 ;Wolsy 1975)以下陈述是等价的:
(a) M = ( [ n ] , I K ) \mathcal{M}=\left([n], \mathcal{I}_K\right) M=([n],IK) 是一个阵向量。
(b) K K K 的强覆盖数最多为2。
c) M \mathcal{M} M 有一个独特的天花板

示例53:(继续例5)由 4 x 1 + 5 x 2 + 6 x 3 + 7 x 4 + 9 x 5 ≤ 13 4 x_1 + 5 x_2 + 6 x_3 + 7 x_4 + 9 x_5≤13 4x1+5x2+6x3+7x4+9x513 定义的背包对应的独立系统有两个上限,即 { 3 , 4 } \{3,4\} {3,4} { 1 , 5 } \{1,5\} {1,5}

因此,根据定理52 ( c ) (c) (c),对应的背包不是一个阵。或者,这可以用定理52(b)来表示,因为我们之前已经看到 C 1 = { 1 , 2 , 3 } , C 5 = { 2 , 5 } , C 6 = { 3 , 5 } 和 C 7 = { 4 , 5 } C1 = \{1,2,3\}, C5 = \{2,5\}, C6 =\{3,5\}和C7 =\{4,5\} C1={1,2,3}C5={2,5}C6={3,5}C7={4,5} 是背包的强覆盖物。

除了前面提到的两个结果,Amado和Barcia(1993)考虑了一个阵阵的F族,他们描述了一个背包 K K K 是否包含在 F F F 中。如果背包不等式的系数几乎相同,这些发现就特别有趣,因为这样的背包通常很难解决,参见Martello和Toth(1990,第2.10节)。对于这样的背包,Amado和Barcia用他们的发现加强了背包问题的标准LP松弛。

虽然前面的结果使用了背包不等式来表征 K K K 是否为一个拟阵,但我们也可以提出相反的问题:哪些拟阵是可以被背包表征的?因此,我们需要将上一节中为图定义的阈值属性应用于拟阵。

独立系统 ( F , I ) (F, \mathcal{I}) (F,I) 存在 a ∈ R + F a \in \mathbb{R}_{+}^F aR+F and β ∈ R + \beta \in \mathbb{R}_{+} βR+ , 例如 I = { I ⊆ F : a ( I ) ≤ β } \mathcal{I}=\{I \subseteq F: a(I) \leq \beta\} I={IF:a(I)β} ,即 I \mathcal{I} I 可被背包表征,则称为阈值。为了表征一个阵点是否为阈值,我们需要一个单调性属性。

定义54:设 ( F , I ) (F, \mathcal{I}) (F,I) 是一个独立系统,设 R , S ⊆ F R, S \subseteq F R,SF ,我们说 R R R 大于等于 S S S ,在公式 R ≥ S R≥S RS 中,对于所有的 T ⊆ F \ ( R ∪ S ) T \subseteq F \backslash(R \cup S) TF\(RS)

定理55:(Giles and Kannan 1980)一个拟阵是阈值当且仅当它是3单调的。

定理55扩展了定理52,得到了具有拟阵性质的背包的第三个特征。对任意独立系统的概括,以及对问题(问题22)的回答,是开放的。


十、Integer knapsacks

对于整数背包问题,变量是有界的,值为整数值:
在这里插入图片描述

一些作者,如Kellerer等人(2004),区分了 u i < ∞ u_i <∞ ui< 的有界背包问题和 u i = ∞ u_i =∞ ui= 的无界背包问题。观察到每个无界背包问题都可以通过设置 u i : = ⌊ β / a i ⌋ u_i:=\left\lfloor\beta / a_i\right\rfloor ui:=β/ai 使之有界,前提是 a i > 0 a_i > 0 ai>0 。因此,用 u i − x i u_i−x_i uixi 补充变量 x i x_i xi 是有效的,这将观察2推广到整数情况。因此,我们假设假设1适用于本节的其余部分。

请注意,我们总是可以将(有界)整数背包问题简化为相应的二进制问题,方法是将每个整数项i替换为它的 u i u_i ui 二进制副本,但在这种转换中可能会丢失一些结构,并且公式可能会变得很大。此外,虽然每个二进制≤-背包都可以转换为二进制≥-背包(见【论文阅读】(2020)Knapsack polytopes: a survey(上)的 3.1 节),但如果考虑整数背包,则不存在这种转换。原因是整数≤-knapsacks是有界的,而整数≥-knapsacks是无界的。因此,就产生了不同的多面体性质,例如Yaman(2007)所研究的。

与二进制情形类似,整数背囊多面体 conv ⁡ { x ∈ Z + n : a ⊤ x ≤ β } \operatorname{conv}\left\{x \in \mathbb{Z}_{+}^n: a^{\top} x \leq \beta\right\} conv{xZ+n:axβ} for a ∈ Z + n a \in \mathbb{Z}_{+}^n aZ+n 可以通过研究主背负式多面体的面 conv ⁡ { x ∈ Z + n : ∑ i = 1 n i x i ≤ n } \operatorname{conv}\left\{x \in \mathbb{Z}_{+}^n: \sum_{i=1}^n i x_i \leq n\right\} conv{xZ+n:i=1nixin} 来找到,参见例如,Aráoz et al (2003), Shim et al(2017),或Tyber and Johnson(2010)。虽然一些定义不等式的方面是已知的,但我们还远远没有完全理解这个多面体。因此,对于二元背包多面体,研究特定的背包多面体是必要的,以找到强切割面。

与经典的二元背包问题相比,只有很少的整数背包问题的多面体研究,参见,例如Atamtürk (2003, 2005), Ceria等人(1998),Pochet和Weismantel(1998),或Pochet和Wolsey(1995)。Ceria等人

(1998)讨论了最小覆盖不等式在一般整数情况下的推广。根据他们的定义,一个(整数)封面是一个集 C ⊆ [ n ] C⊆[n] C[n] 与属性 ∑ i ∈ C a i u i > β \sum_{i \in C} a_i u_i>\beta iCaiui>β

如果没有合适的 C C C 子集作为覆盖,则覆盖 C C C 称为最小覆盖。类似于二进制情况,每个最小整数覆盖 C C C 都会产生(最小)(整数)覆盖不等式

在这里插入图片描述

这说明不是所有的变量(最小)覆盖可以同时在他们的上界。然而,与二进制设置相反,这些不等式一般不定义整数背包集的整数规划公式。

除了上面的最小覆盖不等式的直接推广,Atamtürk(2005)讨论了一个更丰富的不等式类别,从覆盖 C ⊆ [ n ] C⊆[n] C[n] ,表明对于任何 ρ > 0 ρ > 0 ρ>0 ,不等式

在这里插入图片描述

λ = ∑ i ∈ C u i − β \lambda=\sum_{i \in C} u_i-\beta λ=iCuiβ ,对整数背负式多面体有效。对于 ρ = a ℓ \rho=a_{\ell} ρ=a 的特殊情况,当式(20)是facet定义可用时,使用

在这里插入图片描述

定理56:(Atamtürk 2005) 让 C ⊆ [ n ] C⊆[n] C[n] 是一个覆盖,让 ℓ ∈ C \ell \in C C 满足 μ = u ℓ a ℓ − λ ≥ 0 μ=u_{\ell} a_{\ell}-\lambda \geq 0 μ=uaλ0 。不等式(20), ρ = a ℓ \rho=a_{\ell} ρ=a 当且仅当 a i ≥ min ⁡ { λ , κ i ℓ a ℓ − r } a_i \geq \min \left\{\lambda, \kappa_{i \ell} a_{\ell}-r\right\} aimin{λ,κiar} 对于所有的 i ∈ C \ { ℓ } i \in C \backslash\{\ell\} iC\{} ,其中 r = μ − ⌊ μ / a ℓ ⌋ a ℓ r=\mu-\left\lfloor\mu / a_{\ell}\right\rfloor a_{\ell} r=μμ/aa

此外,还存在一些具有完整线性描述的整数背包多形体族。这些族将弱超增和顺序二元背包多形体推广到有界背包多形体的一般情况(见7.2节)。

如果 ∑ j = 1 i − 1 a j u j ≤ a i \sum_{j=1}^{i-1} a_j u_j \leq a_i j=1i1ajujai 对所有的 i ∈ [ n ] i \in[n] i[n] , 定义不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 的有界背包, 且上界向量u称为弱超增。用这个定义,二元情况(第7.2节)的经典结果可以推广到任意整数变量的情况。

定理57:(Gupte 2016)弱超递增(不一定是二进制)背包多形体完全由盒约束和 O ( n ) O(n) O(n) 不等式描述。给出了定义不等式的facet的显式构造方案。

此外,序列二进制背包多形体的完全线性描述,即对于每个 i ∈ [ n − 1 ] i∈[n−1] i[n1] 具有 a i ∣ a i + 1 a_i \mid a_{i+1} aiai+1 的二进制背包多形体,可以推广到任意整数背包多形体,满足可除性,参见Pochet和Weismantel(1998)。一般来说,这种描述是非常复杂的。对于没有明确的变量上界的特殊情况,即 x i ≤ ⌊ β a i ⌋ x_i \leq\left\lfloor\frac{\beta}{a_i}\right\rfloor xiaiβ , Marcotte(1985)对线性大小的背包多形体给出了完整的线性描述。

定理58:(Marcotte 1985)设 a ∈ N n a \in \mathbb{N}^n aNn 使 a i ∣ a i + 1 a_i \mid a_{i+1} aiai+1对每个 i ∈ [ n − 1 ] i∈[n−1] i[n1] 成立。那么 conv ⁡ { x ∈ Z + n : a ⊤ x ≤ β } \operatorname{conv}\left\{x \in \mathbb{Z}_{+}^n: a^{\top} x \leq \beta\right\} conv{xZ+n:axβ} 完全由非负约束和描述:

在这里插入图片描述

如果整数背包多形体的完整线性描述是未知的,则可以尝试通过添加不等式或加强已有的不等式来改进给定的整数公式。后者可以通过提升来实现。然而,与二进制情况相反,整数提升更丰富,因为初始不等式可能不是通过固定某些变量的上界和下界来推导的,而是通过一些中间值来推导的。例如,可以在Wolsey (1976b)或Easton和Gutierrez(2015)中找到整数提升的详细信息。我们还没有发现背包特定不等式的任何结果。

(Q23) 实践中大多数背包的例子都是二进制形式的。有没有有趣的实用的一般整数的例子?

(Q24) 整数背包的情况如何才能用算法得到最好的处理?


十一、Mixed-integer knapsacks

在本节中,我们将讨论混合整数背包的两种变体:混合整数背包的通用模型和描述包含变量下限和上限约束的广义流模型的变体。

11.1 General mixed-integer knapsacks

一般的混合整数背包是将整数背包直接推广到混合整数的情况,即:

在这里插入图片描述

对于适当维度的非负向量 a , b , u a, b, u a,b,u , 和 u ~ \tilde{u} u~ 。如果 u = 1 u=\mathbb{1} u=1 ,则 X M i x X_{Mix} XMix 称为混合二进制背包,否则称为混合整数背包。

对于混合整数背包,Martin和Weismantel(1997)推导出了重量不等式族

在这里插入图片描述

其中 I ⊆ [ m ] I⊆[m] I[m] J ⊆ [ n ] J⊆[n] J[n] , and r ( I , J ) r(I, J) r(I,J) 是当I和J中的所有变量都固定在其上界时,即 r ( I , J ) = δ − ∑ i ∈ I a i u i − ∑ j ∈ J b j u ~ j r(I, J)=\delta-\sum_{i \in I} a_i u_i-\sum_{j \in J} b_j \tilde{u}_j r(I,J)=δiIaiuijJbju~j

权重不等式对于混合整数背包多形体 P Mix  = conv ⁡ ( X Mix  ) P_{\text {Mix }}=\operatorname{conv}\left(X_{\text {Mix }}\right) PMix =conv(XMix ) 是有效的,Martin和Weismantel描述了他们定义 P M i x P_{Mix} PMix 的facet的情况。

Richard等人(2003a)研究了混合二元背包的情况。除了提供 P M i x P_{Mix} PMix 的基本多面体性质外,他们还研究了连续变量的提升问题,并提出了一个伪多项式时间序列提升算法。在理查德等人

(2003b),同一作者将经典背包通过超加性函数同时提升的概念应用于混合二元背包中连续变量的超线性提升函数。这允许他们为 P M i x P_{Mix} PMix 派生几个方面定义不等式。

Marchand和Wolsey(1999)考虑混合二元背包的特殊情况,只有一个系数为−1的无界连续变量,即 X M i x ′ = { ( x , y ) ∈ { 0 , 1 } m × R + : ∑ i = 1 m a i x i ≤ δ + y } X_{\mathrm{Mix}}^{\prime}=\{(x, y) \in\left.\{0,1\}^m \times \mathbb{R}_{+}: \sum_{i=1}^m a_i x_i \leq \delta+y\right\} XMix={(x,y){0,1}m×R+:i=1maixiδ+y}。注意, X Mix  ′ X_{\text {Mix }}^{\prime} XMix  不是一个经典的背包集,因为连续变量有一个负系数。该模型可以被解释为经典二进制背包的扩展,其中容量界限可以使用多余的变量 y y y 来违反。

与经典的背包问题相反, X Mix  ′ X_{\text {Mix }}^{\prime} XMix  的背包不等式总是定义了 conv ⁡ ( X Mix  ′ ) \operatorname{conv}\left(X_{\text {Mix }}^{\prime}\right) conv(XMix ) 的一个方面,参见Marchand和Wolsey(1999)。此外,他们还得到了进一步的有效不等式。在二元变量的某些系数允许为负的情况下,Marchand和Wolsey导出了几个面定义不等式,并研究了它们的提升问题。

后者的一种变体是 { ( x , y ) ∈ { 0 , 1 } m × Z + n : ∑ i = 1 m x i ≤ ∑ j = 1 n b j y j } \left\{(x, y) \in\{0,1\}^m \times \mathbb{Z}_{+}^n: \sum_{i=1}^m x_i \leq \sum_{j=1}^n b_j y_j\right\} {(x,y){0,1}m×Z+n:i=1mxij=1nbjyj} ,这称为整数容量集。Mazur和Hall(2002)提供了相应多面体的基本多面体性质,并展示了如何减少其公式中的系数,这在预处理中特别有趣。此外,他们提出了基于某些背包封面的面定义不等式。

11.2 Generalized flow models

带变量约束的混合整数背包的最一般形式如下所示。设 N 1 N_1 N1 N 2 N_2 N2 是(不一定不相交的)集合,设 a ∈ R N 1 a \in \mathbb{R}^{N_1} aRN1 b ∈ R N 2 b \in \mathbb{R}^{N_2} bRN2 以及 ℓ ′ , u ′ ∈ R N 1 ∪ N 2 \ell^{\prime}, u^{\prime} \in \mathbb{R}^{N_1 \cup N_2} ,uRN1N2 and δ ∈ R \delta \in \mathbb{R} δR 。具有可变上下界的混合整数背包是一个集合:

在这里插入图片描述

这个定义中的最后一类不等式被称为变量下限和上限约束,因为取决于 y j ′ y_j^{\prime} yj 的值,它们是活跃的 ( y j ′ = 1 ) \left(y_j^{\prime}=1\right) (yj=1)或不活跃 ( y j ′ = 0 ) \left(y_j^{\prime}=0\right) (yj=0),注意我们使用 ℓ ′ \ell^{\prime} u ′ u^{\prime} u 而不是 ℓ \ell u u u 代表边界,因为我们会适应下面这个符号。此外,由于我们将连续变量解释为流量值,因此我们将连续变量称为 f f f -变量。

为了简化这种混合整数背包的分析,人们通常用一个更简单的整数公式将 K F ′ K_F^{\prime} KF into转换为一个集合,我们将在下面介绍这个集合,并在后面介绍(V an Roy and Wolsey 1986)。通过适当地替换变量和边界,他们证明 K F ′ K_F^{\prime} KF 可以表示为

在这里插入图片描述
其中 N : = N 1 ∪ N 2 N:=N_1 \cup N_2 N:=N1N2 ( M 1 , M 2 ) \left(M_1, M_2\right) (M1,M2) N N N 的一个分区。这个模型就是所谓的混合整数背包的标准形式,用 K F K_F KF 表示。特别是,我们可以假设 0 ≤ ℓ ≤ u 0 \leq \ell \leq u 0u

标准形式的混合整数背包可以被解释为一个流网络的模型,由一个节点 v v v 以及一套弧 M 1 M_1 M1 ,点成 v v v 和一组弧形 M 2 M_2 M2 组成。变量 f j f_j fj 对弧 j j j 上的流量值进行建模,变量 y j y_j yj 通过是否可以使用弧 j j j 的广义上下界约束进行编码。

最后, K F K_F KF 中的第一个不等式,我们称之为流量背包不等式,表示节点 v v v 的净流出量不超过 δ δ δ

为了加强公式 K F K_F KF of P F : = conv ⁡ ( K F ) P_F:=\operatorname{conv}\left(K_F\right) PF:=conv(KF) , Van Roy 和 Wolsey (1986)引入了广义流动覆盖切割。一对 ( C 1 , C 2 ) (C_1, C_2) (C1,C2) C 1 ⊆ M 1 C1⊆M1 C1M1 , C 2 ⊆ M 2 C2⊆M2 C2M2 ,称为广义覆盖

在这里插入图片描述
即,将 C 1 C_1 C1 中的变量固定在其上界,将 C 2 C_2 C2 中的变量固定在其下界,违反了流量背包不等式。令

在这里插入图片描述

α + = max ⁡ { 0 , α } \alpha^{+}=\max \{0, \alpha\} α+=max{0,α} 对于实值 α α α

给定广义流动覆盖 ( C 1 , C 2 ) (C_1, C_2) (C1,C2) ,令 L i ⊆ M i \ C i , i ∈ { 1 , 2 } L_i \subseteq M_i \backslash C_i, i \in\{1,2\} LiMi\Ci,i{1,2} 。根据Gu等人(1999b)的表述, ( C 1 , L 1 , C 2 , L 2 ) (C_1, L_1, C_2, L_2) (C1,L1,C2,L2) 的简单广义流动覆盖不等式为

在这里插入图片描述
这是 P F P_F PF 的一个有效不等式。此外,存在所谓的扩展广义流动覆盖不等式的推广,参见,例如Gu et al (1999b)。

Van Roy 和 Wolsey (1986)研究了广义流动覆盖不等式定义 P F P_F PF 面的情况,并研究了它们的分离问题。Stallaert(1997)推导出一类不等式,称为 μ μ μ-不等式,它是流覆盖不等式的补充。他进一步分析了它们的分离问题,并通过实验发现,使用广义流覆盖不等式和 μ μ μ 不等式可以将完整性差距缩小 75 % 75\% 75% ,而使用这两类中的任何一种只能将差距缩小 65 % 65\% 65% 。Gu等人(1999b)对广义流动覆盖不等式的提升问题进行了广泛的研究,并进行了数值实验。


十二、Variants of the knapsack problem

在本节中,我们简要地讨论背包问题的变体的多面体方面。对于第12.1-12.4节中讨论的变体的算法方面,我们再次建议读者参考Kellerer et al(2004)的书。

12.1 Multidimensional knapsack problem

多维背包集 K m u l t K_{mult} Kmult m m m 个背包约束来定义,而不是单一的背包约束,即

在这里插入图片描述

其中 A ∈ Z + m × n A \in \mathbb{Z}_{+}^{m \times n} AZ+m×n and b ∈ Z + m b \in \mathbb{Z}_{+}^m bZ+m 。也就是说,一个多维背包 K m u l t K_{mult} Kmult m m m 个背包 K i K_i Ki 的交点,因此是一个独立系统。相反,每个独立系统都有一个多维背包的表示,对于每个最小相关集 C C C ,通过切割 x ( C ) ≤ ∣ C ∣ − 1 x(C)≤|C|−1 x(C)C1 来表示。因此,与独立系统相关的一般多面体上的每个结果也适用于多维背包多面体 P M P_M PM ,例如,冲突超图中通过最大团来表征facet,参见Easton et al(2003)。

通过分别考虑 A x ≤ b Ax≤b Axb 中的每个约束,将经典背包的覆盖不等式的概念直接转移到多维情况。Bektas和O˘guz(2007)提出了一个IP模型来分离 P M P_M PM 的违规覆盖不平等。Kaparis和Letchford(2008)研究了解除的覆盖不平等。Balas和Zemel(1984)通过将 K m u l t K_{mult} Kmult 中单个约束的最小覆盖直接推广到 A x ≤ b Ax≤b Axb 的完整系统的最小覆盖,表明定义 P M P_M PM 不等式的每一个非平凡方面都可以通过对适当选择的集合 S ⊆ [ n ] S⊆[n] S[n] 的变量 x i x_i xi , i ∈ S i∈S iS 进行互补,检测被补背包的合适最小覆盖,并(依次或同时)提升相应的最小覆盖不等式来找到。因此,找到 P M P_M PM 的完整线性描述的理论机制是已知的。

Martin和Weismantel(1998)考虑了多维背包问题的一个更一般的版本,其中变量的上界是任意的正整数。也就是说,他们研究了几个整数背包交叉的情况,而不是二进制的情况,当然,他们的结果也适用于二进制的情况。他们推导了所谓的可行集不等式,它对 ⋂ i = 1 m K i \bigcap_{i=1}^m K_i i=1mKi 有效,但不一定对每个 K i K_i Ki 都有效,他们推导了这些不等式提升系数的边界。如果背包 K i , i ∈ [ m ] K_i, i∈[m] Ki,i[m] 可以这样标记:当 ∣ i − k ∣ ≥ 2 |i−k|≥2 ik2 时, K i ∩ K k = ∅ K_i \cap K_k=\varnothing KiKk= ,Martin和Weismantel就完全表征了定义某个 conv( K i K_i Ki) 的不等式的一个面是否也定义了 P M P_M PM 的一个面。

12.2 Cardinality constrained knapsack problem

基数约束背包集 K ≤ K_≤ K由一个普通背包不等式加上最多可选择 K K K 个物品的附加约束组成,即:
在这里插入图片描述

受基数约束的背包是具有两个背包约束的多维背包的一种特殊情况,即两个经典背包的交集。Louveaux和Weismantel(2008)研究了基数约束背包的多面体方面,他们开发了所谓的不完全集不等式,以利用在单独背包中不存在的交集结构。此外,第10.1节中提到的所有结果也适用于基数约束的情况。

Glover和Sherali(2008)以及Sherali和Glover(2008)推导了利用基数约束的特殊覆盖不等式。对于包含基数约束的整数程序,Bienstock(1996)引入了混合整数舍入不等式、析取切分和临界集不等式。Zeng和Richard (2011a, b)通过在不同变量上引入几个基数约束,对上述问题进行了推广,并给出了提升方案。如果所有变量都是连续的,de Farias Jr 和 Nemhauser (2003)研究了 K ≤ K_≤ K 的一种变体,其中基数约束被替换为最多 K K K 个变量获得正值的条件。

12.3 Generalized upper bound constraints

给定一个由变量 x ∈ 0 , 1 n x∈{0,1}^n x0,1n 的索引组成的集合 J ⊆ [ n ] J \subseteq[n] J[n] ,一个约束 x ( J ) ≤ 1 x(J) \leq 1 x(J)1 被称为一个广义上界(GUB)。在许多应用中,这种约束与背包约束结合在一起。为此,考虑将 [ n ] [n] [n] 划分为 k k k 个集合 Q 1 , … , Q k Q_1,…,Q_k Q1Qk ,即 Q 1 ∪ ⋅ ⋅ ⋅ ∪ Q k = [ n ] Q_1∪···∪Q_k = [n] Q1⋅⋅⋅Qk=[n] Q i ∩ Q j = ∅ Q_i \cap Q_j=\varnothing QiQj= 对于所有不相同的 i , j ∈ [ k ] i, j∈[k] i,j[k] 。那么带有(非重叠的)GUB约束的背包集是

在这里插入图片描述

注意,通过添加基数为 1 1 1 Q i Q_i Qi,假设 Q i Q_i Qi 覆盖集合 [ n ] [n] [n] 而不损失一般性。

与基本情况类似,集合 C ⊆ [ n ] C⊆[n] C[n] 是一个(最小)GUB覆盖,如果 C C C 是背包约束的(最小)覆盖,并且 C C C 中没有两个元素属于同一个 Q i Q_i Qi ,即 a ( C ) > β a(C) > β a(C)>β ∣ C ∩ Q i ∣ ≤ 1 |C∩Qi |≤1 CQi1 对所有的 i ∈ [ k ] i∈[k] i[k]。对应的GUB覆盖不等式 x ( C ) ≤ ∣ C ∣ − 1 x(C)≤|C|−1 x(C)C1 对于 P G U B : = conv ⁡ ( K G U B ) P_{\mathrm{GUB}}:=\operatorname{conv}\left(K_{\mathrm{GUB}}\right) PGUB:=conv(KGUB) 有效。注意,普通的覆盖不等式也是有效的,但要么是GUB覆盖,要么是整数解的冗余。

Wolsey(1990)考虑了带有GUB约束的背包,其中背包包含正系数和负系数,并引入了GUB覆盖。他还研究了可以完全描述 P G U B P_{GUB} PGUB 的特殊情况,GUB覆盖不等式的扩展,以及流覆盖的GUB模拟。请注意,负系数可以通过互补来消除,以达到上面的形式,请参阅Johnson和Padberg(1981)了解更多细节。

Sherali和Lee(1995)研究了PGUB的多面体性质。特别是,它们描述了最小的GUB覆盖不等式何时定义方面,并研究了这些不等式的顺序和同时提升何时定义方面。此外,他们还推导了提升系数的边界,可以在 O ( n k ) O(nk) O(nk) 时间内计算。这些结果提供了第5.4节中讨论的经典情况的推广。特别地,对于没有GUB约束的背包集,即 ∣ Q i ∣ = 1 \left|Q_i\right|=1 Qi=1 for all i ∈ [ k ] i \in[k] i[k] , Sherali和Lee(1995,命题4.2)提供的刻画可简化为定理38。

Nemhauser和V ance(1994)独立地描述了按顺序解除的最小GUB覆盖不平等的情况

在这里插入图片描述

定义 P G U B P_{GUB} PGUB 的一个方面。这里, α j = π j + 1 \alpha_j=\pi_j+1 αj=πj+1 α j = π j \alpha_j=\pi_j αj=πj (见定理28和38)。该结果既改进了无GUB约束的背包问题的推论40,又推广了有GUB约束的推论40。

Gu等人(1999b)讨论了解除GUB覆盖不等式的算法和实践方面。例如,人们可以向上或向下解除这些不平等。然而,举升问题变得更加复杂,见Gu(1995)。

Gokce和Wilhelm(2015)考虑≥-形式,即 a ⊤ x ≥ β a^{\top} x \geq \beta axβ K G U B K_{GUB} KGUB 和定义(强) α-覆盖 C α C_α Cα 。它们表明了相应的不平等 ∑ j ∈ C α x j ≥ α \sum_{j \in C_\alpha} x_j \geq \alpha jCαxjα 是有效的,当它们是facet定义时是有特征的。此外,它们还描述了 α α α 盖层的序贯提升和序贯独立提升。

(Q25) 我们能把上面提到的一些结果转移到集合 Q i Q_i Qi 重叠的情况吗?这将对应于一个稳定集和背包问题的组合。

12.4 Precedence constrained knapsack problem

背包问题的这种变体为普通的二元背包问题增加了优先约束 x i ≥ x j x_i≥x_j xixj ,强制要求只有在项目 i i i 也被拿走时,项目 j j j 才能被添加到背包中。这就引出了背包集

在这里插入图片描述

其中 A \mathcal{A} A 是由来自 [ n ] [n] [n] 的有序对组成的集合。注意,集合 A \mathcal{A} A 可以解释为具有节点集 [ n ] [n] [n] 的有向图的弧。

Boyd (1993b)研究了 P prec  : = conv ⁡ ( K prec  ) P_{\text {prec }}:=\operatorname{conv}\left(K_{\text {prec }}\right) Pprec :=conv(Kprec ) 的基本多面体性质,并推广了覆盖不等式的概念。此外,基于优先约束所定义的图,Boyd导出了facet定义不等式的类。Park(1997)讨论了优先约束背包问题的解除覆盖不等式,并描述了在哪些情况下这些不等式是 P p r e c P_{prec} Pprec 的面定义。本文还讨论了 K p r e c K_{prec} Kprec 的LP弛豫性质。van de Leensel等人(1999)研究了提升几类有效不等式的复杂性,例如 (1, k)-configurations 或(概括)覆盖不等式。

12.5 Generalized assignment problem

广义分配问题是多背包问题的一个变体,在调度中有应用。有 m m m 个背包(机器), n n n 个物品中的每一个都必须分配给一台机器,从而形成背包集

在这里插入图片描述

Gottlieb和Rao (1990b)提供了 P G A P : = conv ⁡ ( K G A P ) P_{\mathrm{GAP}}:=\operatorname{conv}\left(K_{\mathrm{GAP}}\right) PGAP:=conv(KGAP) 的多面体研究。他们还研究了KGAP的LP弛豫性质。在Gottlieb和Rao(1990a),同样的作者专注于 (1, k)-configurations 的面定义性质。Trick(1992)还研究了 K G A P K_{GAP} KGAP 的LP松弛,并利用他的发现开发了有效的原始启发式。


十三、Conclusions

在本文中,我们概述了背包多面体的多面体性质,主要是最小覆盖不等式的提升,进一步的有效不等式和完整的线性描述。在这篇文章中,我们收集了超过20个关于背包多形体的开放问题,包括关于某些不等式的理论强度、识别或分离特定类别不等式的复杂性,或可以获得完整线性描述的背包多形体类别的问题。我们希望这些问题能从理论和计算两方面构成背包多面体进一步研究的基础。

例如,对于完全线性描述,本文揭示了本质上已知的三个背包多形体族的外部描述:(i)背包不等式的系数很小或很大,(ii)背包不等式的连续系数是可除的,或(iii)系数在维数上呈指数增长。但是,如果背包不等式的系数是次指数增长的,而连续的系数是不可除的,那么人们对facet结构能说什么呢?在计算方面,我们已经看到,高维背包多面体可能不能很好地近似于解除覆盖不等式。关于这个观察,人们可能会问:哪些类型的不等式是必要的,以得到一个更紧密的松弛背包多面体,特别是,这些不等式可以有效地分离?

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/110326.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

JavaFx TreeView TreeItem 设置额外属性

在使用JavaFx 编写GUI程序时&#xff0c;不可避免的需要创建一个树组件,下面是一个简单的树组件的代码。 import javafx.application.Application; import javafx.scene.Scene; import javafx.scene.control.TreeItem; import javafx.scene.control.TreeView; import javafx.s…

clickhouse笔记05--快速部署3节点集群

clickhouse笔记05--快速部署3节点集群1 介绍2 方法步骤2.1 部署 zookeeper 集群2.2 拉起 clickhouse 集群2.3 测试集群3 注意事项4 说明1 介绍 clickhouse笔记01–快速部署clickhouse 介绍了如何快速部署单节点clickhouse服务&#xff0c;本文基于该博文继续介绍如何快速部署3…

Java进阶—JUC编程

1、线程和进程 获取CPU核数 /*** author java小豪* version 1.0.0* date 2022/12/15* description 测试*/ public class Test {public static void main(String[] args) {// 获取CPU核数// CPU 密集型&#xff0c;IO密集型System.out.println(Runtime.getRuntime().available…

响应式营销策划文化传媒公司网站模板源码

模板信息&#xff1a; 模板编号&#xff1a;8071 模板编码&#xff1a;UTF8 模板颜色&#xff1a;蓝色 模板分类&#xff1a;设计、广告、文化、影视 适合行业&#xff1a;影视传媒类企业 模板介绍&#xff1a; 本模板自带eyoucms内核&#xff0c;无需再下载eyou系统&#xf…

qt5实现pdf阅读器(三)——pdfjs

目录 1、参考 2、实现 3、开发记录 1、参考 使用Qt的WebEngine和javascript的pdf.js模块构建的PDF查看器。 参考链接1&#xff1a;GitHub - Archie3d/qpdf: PDF viewer widget for Qt 参考链接2&#xff1a;GitHub - yshurik/qpdfjs: Desktop PDF Viewer based on Qt and…

讯飞听见SaaS服务迈入全新时代

配图来自Canva可画 随着数字化时代的来临&#xff0c;国内各企业为了提升行业竞争力&#xff0c;纷纷开始利用数字化技术&#xff0c;来实现以降本增效为核心的数字化转型&#xff0c;得益于此&#xff0c;助力企业数字化转型升级的SaaS也开始进一步升温。 众所周知&#xff…

【代码审计-2】PHP框架MVC类文件上传断点测试挖掘

1.文件上传漏洞挖掘&#xff1a; (1)关键字搜索&#xff08;函数、键字、全局变量等&#xff09;&#xff1a;比如$_FILES&#xff0c;move_uploades_file等 (2)应该功能抓包&#xff1a;寻找任何可能存在上传的应用功能点&#xff0c;比如前台会员中心&#xff0c;后台新闻添…

电力系统两阶段随机优化(Matlab实现)

目录 目录 1 概述 2 单级随机优化算法 2.1 随机化-最小化 2.2 随机逐次凸近似 &#xff08;SCA&#xff09; 3 两级随机优化算法 3.1 批处理算法 3.2 在线算法 4 Matlab代码实现 1 概述 在与随机系统状态向量关联的两阶段随机优化问题中&#xff0c;优化变量分为两组…

Web前端105天-day32-HTML5_CORE

HTML5CORE02 目录 前言 一、复习 二、拖拽 三、上传服务器 四、Canvas 五、地图 总结 前言 HTML5CORE02学习开始 一、复习 跨域 浏览器的同源策略导致在网页中, 通过 AJAX 发送网络请求时, 默认只能向同源的服务器请求同源: 协议 端口号 域名 三者都相同产生跨域的原因…

RocketMQ疑难杂症之No route info of this topic解决方案

成因&#xff1a; 由于配置了 docker 虚拟 IP&#xff0c;导致 brocker 总是代理到 docker 的虚拟 IP 上。 原理&#xff1a; RocketMQ 的 broker 启动类 org.apache.rocketmq.broker.BrokerStartup 启动的时候会读取代码中的默认配置&#xff0c;关于 broker 的配置在 org.apa…

【关于时间序列的ML】项目 8 :使用 Facebook Prophet 模型预测股票价格

&#x1f50e;大家好&#xff0c;我是Sonhhxg_柒&#xff0c;希望你看完之后&#xff0c;能对你有所帮助&#xff0c;不足请指正&#xff01;共同学习交流&#x1f50e; &#x1f4dd;个人主页&#xff0d;Sonhhxg_柒的博客_CSDN博客 &#x1f4c3; &#x1f381;欢迎各位→点赞…

30.深度学习模型压缩方法-4

30.1 低秩分解 基于低秩分解的深度神经网络压缩与加速的核心思想是利用矩阵或张量分解技术估计并分解深度模型中的原始卷积核 卷积计算是整个卷积神经网络中计算复杂 度 最 高 的 计 算 操 作,通 过 分 解4D 卷积核张量,可以有效地减少模型内部的冗余性此外对于2D的全 连…

Hive+Spark离线数仓工业项目实战--项目介绍及环境构建(1)

项目简介 通过大数据技术架构&#xff0c;解决工业物联网制造行业的数据存储和分析、可视化、个性化推荐问题。一站制造项目主要基于Hive数仓分层来存储各个业务指标数据&#xff0c;基于sparkSQL做数据分析。核心业务涉及运营商、呼叫中心、工单、油站、仓储物料。 推荐教程…

DSP_TMS320F28377D_eCAP学习笔记

博主学习eCAP的使用主要是用于处理霍尔传感器&#xff0c;计算电机的电角度以及角速度。首先还是看了点哔哩哔哩的学习视频。 eCAP介绍 脉冲量的输入是在数字控制系统中最常见的一类输入量&#xff0c;控制器专门设置了脉冲捕获模块 (eCAP)来处理脉冲量&#xff0c;通过脉冲捕…

路由器的工作原理(计算机网络-网络层)

目录 路由器的构成 转发和路由选择的区别 典型的路由器结构 交换结构 输出端口 路由器与交换机的比较 两种基于存储转发的分组交换设备的比较 交换机和路由器各有的应用场合 三层交换机 三层交换机的应用 路由器的构成 路由器的任务 路由器是一种具有多个输入端口和多…

MT8385 Android AB分区系统升级(命令模式)

AB系统分区升级使用的是update_engine,RecoverySystem 只适用于单分区的系统升级 1.解压开update.zip 可以查看到palyload的属性 2.使用ADB命令update_engine_client即可对AB分区进行升级 使用adb shell 命令进行升级 update_engine_client --payload xxx --update --header…

【TypeScript】TS类型声明(二)

&#x1f431;个人主页&#xff1a;不叫猫先生 &#x1f64b;‍♂️作者简介&#xff1a;前端领域新星创作者、华为云享专家、阿里云专家博主&#xff0c;专注于前端各领域技术&#xff0c;共同学习共同进步&#xff0c;一起加油呀&#xff01; &#x1f4ab;系列专栏&#xff…

k8s HPA升级 KEDA 基于prometheus的数据指标进行弹性伸缩

说明&#xff1a;KEDA有啥用&#xff0c;相对HPA有啥优势。HPA针对于cpu,内存来进行弹性伸缩&#xff0c;有点不太精确。KEDA可以接入prometheus&#xff0c;根据prometheus的数据指标进行弹性伸缩&#xff0c;相比更加的精准实用。 安装k8s环境部署prometheus 创建ns&#xf…

【软件测试】那些35岁以上的测试人哪去了?

目录&#xff1a;导读前言一、Python编程入门到精通二、接口自动化项目实战三、Web自动化项目实战四、App自动化项目实战五、一线大厂简历六、测试开发DevOps体系七、常用自动化测试工具八、JMeter性能测试九、总结&#xff08;尾部小惊喜&#xff09;前言 先根据大体年龄阶段…

B+树详解,一次就懂

⭐注意&#xff1a;不会直接讲 B树的结构&#xff0c;会从最简单的二叉树开始讲起来。如果认真看完&#xff0c;我想你对树类型的数据结构的理解又上了一个新的台阶。 ⭐如果有误&#xff0c;请大家指出。下文均是在B站学习的过程中&#xff0c;总结的笔记和心得体会 索引结构 …