当前位置: 首页 > news >正文

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

文章目录

  • 一、Abstract 摘要
  • 二、Introduction 介绍
  • 三、General polyhedral structure 一般多面体结构
    • 3.1 Basic properties 基本性质
    • 3.2 Covers 覆盖不等式
  • 四、Binary formulations based on strong covers 基于强覆盖的二元公式
  • 五、Lifting 提升
    • 5.1 Sequential up-lifting
    • 5.2 Sequential down-lifting
    • 5.3 Sequence independent lifting and superadditivity 序列独立提升和超加性
    • 5.4 Complete characterization of facets from simple LCIs 从简单LCIs中完全描述facet


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


一、Abstract 摘要

  • 0/1背包多面体是所有0/1向量的凸包,这些向量满足给定的具有非负系数的线性不等式。
  • 本文全面综述了背包多面体。
  • 我们讨论了基本的多面体性质、(提升)覆盖和其他有效的不等式、已知完整线性描述的情况、小维度的几何性质以及与独立系统的联系。
  • 我们还讨论了(混合)整数背包多面体和变体的推广。

二、Introduction 介绍

0/1背包问题包括选择具有给定权重和利润的物品,使得权重之和不超过一定容量,并且利润之和最大化。如果n项的权重由 a = ( a 1 , … , a n ) ⊤ ∈ R + n a=\left(a_1, \ldots, a_n\right)^{\top} \in \mathbb{R}_{+}^n a=(a1,,an)R+n ,容量为 β ∈ R + \beta \in \mathbb{R}_{+} βR+ 和利润 p = ( p 1 , … , p n ) ⊤ ∈ R + n p=\left(p_1, \ldots, p_n\right)^{\top} \in \mathbb{R}_{+}^n p=(p1,,pn)R+n ,可以写成 max ⁡ { p ⊤ x : x ∈ K a , β } \max \{p^{\top} x:x \in K^{a,\beta}\} max{px:xKa,β}

在这里插入图片描述

是与重量 a a a 和容量 β \beta β 相关的0/1背包集。

背包问题在离散优化中具有根本重要性,并且具有使其独特的几个特征。

首先,它是弱NP难的,但可以通过伪多项式时间 O ( n β ) O(n\beta) O(nβ) 中的动态规划来解决 ( if a ∈ Z + n , β ∈ Z + a \in \mathbb{Z}_{+}^n, \beta \in \mathbb{Z}_{+} aZ+n,βZ+ ) or O ( n p ˉ ) \mathcal{O}(n \bar{p}) O(npˉ) ( if p ∈ Z + p \in \mathbb{Z}_{+} pZ+ ),其中p是最优目标的某个上界,例如 p 1 + . . . + p n p_1+...+p_n p1+...+pn ,参见例如Kellerer等人(2004,引理2.3.1和引理2.3.2)。因此,就其复杂性而言,它是一个中间问题。

第二,这个问题很容易描述,但根据权重和利润有着丰富的结构。

第三,它在一般二进制程序 m a x { p ⊤ x : A x ≤ b , x ∈ { 0 , 1 } n } max \left\{p^{\top} x: A x \leq b, x \in\{0,1\}^n\right\} max{px:Axb,x{0,1}n}。为了导出此类程序的切割平面,例如,可以取任意行 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 系统 A x ≤ b Ax\le b Axb 并研究由 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 定义的背包问题,可参照 Crowder 等人(1983)。

第四,由于它具有组合结构(或与其他组合约束相结合),但具有一般系数,因此它是组合优化和整数规划之间的桥梁。

因此,关于背包问题的文献非常多,研究算法、多面体和计算问题。例如,在Kellerer等人(2004年)或Martello和Toth(1990年)的书中,算法方面得到了处理,但我们不知道多面体研究的概述。本文的一个动机是特殊背包问题的出现,研究了其多面体结构(Hojny和Pfetsch 2019),甚至找到了完整的线性描述(Hojny2018)。在此背景下,我们(未成功)搜索了一篇关于背包多面体的综述文章。

在本文中,我们试图通过对0/1背包多面体的概述来弥补这一差距:

在这里插入图片描述

即:背包问题所有可行解的凸包。每当背包不等式从上下文中明确时,我们分别写 K K K P P P 而不是 K a , β K^{a,\beta} Ka,β P a , β P^{a,\beta} Pa,β

在下文中,我们专注于研究最多的0/1情况,但也涵盖了背包问题和其他约束的整数情况或组合。在整数情况下,变量可以达到任何非负整数值(可能有上限)。这种情况的实际重要性较小,并且这种情况下存在的多面体结果较少。然而,这具有理论意义,因为每个有界整数程序的可行域(可能具有松弛变量)可以表示为满足框约束和单个线性方程的整数点,见Bradley(1971)。此外,存在大量的背包变体,难以涵盖所有结果;然而,本文的后面部分也涵盖了一系列变体。

本文所选主题及其陈述自然是主观的。对于重要的结果,我们力求全面,并在路上提到许多其他的结果。

由于这是一篇综述性论文,因此重点关注文献的结果。我们为一些所需空间相当小的结果提供了证明。此外,我们用一些新材料补充了文献结果。例如,Sects中的大多数内容。2和6以及引理20、命题22和推论31都是新的。此外,我们提供了运行示例,并以统一的符号表示所有内容。最后,我们提出了一些开放性问题,以促进进一步的研究。

在本文中,我们假设多面体理论和整数规划的基本知识。有很多书提供了介绍,例如Schrijver(1986)、Nemhauser和Wolsey(1999)或Korte和Vygen(2012)。

如果没有不同的表述, n n n 是一个正整数,表示背包多面体周围空间的维数。我们分别用 N , Z \mathbb{N}, \mathbb{Z} N,Z R \mathbb{R} R 表示所有正整数、整数和实数的集合; Z + \mathbb{Z}_+ Z+ R + \mathbb{R}_+ R+ 分别表示非负整数和实数的集合。

前k个正整数的集合 { 1 , … . k } \{1,….k\} {1,.k} [ k ] [k] [k] 表示, [ 0 ] = ∅ [0]=\varnothing [0]= 。此外,我们定义 [ k ] 0 : = [ k ] ∪ { 0 } [k]_0:=[k] \cup\{0\} [k]0:=[k]{0}。集合 S ⊆ [ n ] S ⊆ [n] S[n] 的特征向量用 χ S \chi^{S} χS 表示。特别地,我们用 1 \mathbb{1} 1 表示全1向量 χ [ n ] \chi^{[n]} χ[n] ,用 e i e^i ei 表示第 i i i 个单位向量 χ { i } \chi^{\{i\}} χ{i}

给定集合 S ⊆ [ n ] S \subseteq[n] S[n] ,我们通过 x ( S ) x(S) x(S) 缩写了总数 ∑ i ∈ S x i \sum_{i \in S} x_i iSxi

这特别意味着 x ( ∅ ) = 0 x(\varnothing)=0 x()=0


三、General polyhedral structure 一般多面体结构

在本节中,我们回顾了任意0/1背包多面体的基本多面体性质。

而第3.1节提供了基本属性,如可以从平凡的不等式出发,推导出维度和facet;3.2节重点讨论了基于覆盖不等式的背包多形体的面。所提出的研究结果的大部分已经在20世纪70年代发表。

3.1 Basic properties 基本性质

在说明背包多面体的多面体性质之前,我们讨论了关于数据 a ∈ R n a∈\mathbb{R}^n aRn β ∈ R β∈\mathbb{R} βR 而不失一般性的基本假设。

假设1:我们有 0 < a i < β 0<a_i<\beta 0<ai<β,对于所有的 i ∈ [ n ] i \in [n] i[n] a 1 + . . . + a n > β a_1+...+a_n > \beta a1+...+an>β

观察2:设 a ∈ R n a∈\mathbb{R}^n aRn β ∈ R β∈\mathbb{R} βR ,假设1可以按照指定的顺序执行以下预处理步骤:

  • 如果 a i < 0 a_i<0 ai<0,我们可以补充变量 x i x_i xi (即用 1 − x i 1-x_i 1xi 替换它),注意,注意,此步骤意味着 β \beta β 值增加 ∣ a i ∣ |a_i| ai
  • 如果 a i = 0 a_i=0 ai=0,我们可以去掉对象 i i i ,因为 P a , β P^{a,\beta} Pa,β 等价于背包多面体在变量 [ n ] \ { i } [n] \backslash\{i\} [n]\{i} 和区间 [ 0 , 1 ] [0,1] [0,1] 上的笛卡尔积,即变量 i i i P a , β P^{a,\beta} Pa,β 的相关结构没有贡献
  • 如果 a i > β a_i>\beta ai>β,我们可以移除对象 i i i ,因为在所有解中 x i = 0 x_i = 0 xi=0

如果 β < 0 β < 0 β<0 β = 0 β = 0 β=0 ,多面体分别为空或仅由零向量组成。

此外,如果 a 1 + ⋅ ⋅ ⋅ + a n ≤ β a_1 +···+ a_n≤β a1+⋅⋅⋅+anβ , 则 P a , β = [ 0 , 1 ] n P^{a,\beta} = [0,1]^n Pa,β=[0,1]n

注意,背包多面体是通过“小于或等于”约束定义的 a ⊤ x ≤ β a^{\top} x \leq \beta axβ

通过补充变量,也可以考虑一个 a ⊤ x ≥ β a^{\top} x \geq \beta axβ 的背包,即对应的多面体是仿射等效的。此外,对于相等约束的背包,它是成立的

在这里插入图片描述

也就是说,理解 ≤ − ≤- ≥ − ≥- 多面体就足够了。因此,我们只考虑 ≤ − ≤- 多面体。关于等式约束背包的基本结果,例如,关于它们的维度和定义不等式的基本方面,请读者参考Lee(1997)。

为了能够描述背包多形体的面定义不等式,有必要知道它们的维数。下面的引理给出了它们的维数的一个刻画,以及平凡不等式 x i ≥ 0 x_i≥0 xi0 x i ≤ 1 , i ∈ [ n ] x_i≤1,i∈[n] xi1,i[n] 为面定义的刻画。

引理3 a ∈ R + n a∈\mathbb{R}^n_+ aR+n β ∈ R + β∈\mathbb{R}_+ βR+

a) 背包集合 K a , β K^{a,\beta} Ka,β 是下单调的,也就是说,如果 x ∈ K a , β x∈K^{a,\beta} xKa,β y ≤ x y≤x yx y ∈ { 0 , 1 } n y∈\{0,1\}^n y{0,1}n ,则 y ∈ K a , β y∈K^{a,\beta} yKa,β 也成立。

b) 假设 H : = { i ∈ [ n ] : a i > β } H:=\left\{i \in[n]: a_i>\beta\right\} H:={i[n]:ai>β} 。则 P a , β P^{a,\beta} Pa,β 的维度数为 n − ∣ H ∣ n-|H| nH

c) 对于每一个 i ∈ [ n ] i∈[n] i[n] x i ≥ 0 x_i≥0 xi0 定义了 P a , β P^{a,\beta} Pa,β 的一个面当且仅当 i ∉ H i \notin H i/H

d) 对于一个 i ∈ [ n ] i∈[n] i[n] x i ≤ 1 x_i≤1 xi1 定义了一个facet当且仅当 i ∉ H i \notin H i/H a i + a j ≤ β a_i + a_j≤β ai+ajβ 对于所有 j ∈ [ n ] \ ( H ∪ { i } ) j \in[n] \backslash(H \cup\{i\}) j[n]\(H{i})

备注4:注意,假设1和引理3(b)保证 P a , β P^{a,\beta} Pa,β 是全维的。相反,如果 a ≥ 0 a≥0 a0 P a , β P^{a,\beta} Pa,β 是全维的,则对于所有 i ∈ [ n ] i∈[n] i[n] , 都有 e i ∈ P a , β e^i∈P^{a,\beta} eiPa,β a i ≤ β a_i≤β aiβ

示例5:考虑下面的背包不等式:

在这里插入图片描述

这将作为贯穿整个调查的运行示例。因为对于所有 i ∈ [ 5 ] i∈[5] i[5] a i ≤ β ai≤β aiβ,引理3中定义的集合H为空。

因此,背包多面体 P a , β P^{a,\beta} Pa,β 是全维的(引理3(b)),所有非负不等式 x i ≥ 0 x_i≥0 xi0 定义切面(引理3©)。

在所有的上限不等式 x i ≤ 1 x_i≤1 xi1 中,只有 x 1 ≤ 1 x_1≤1 x11 定义了一个facet(引理3(d))。

最后,唯一满足背包不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 的 可行的 0/1点 等于 ( 1 , 0 , 0 , 0 , 1 ) ⊤ , ( 0 , 0 , 1 , 1 , 0 ) ⊤ (1,0,0,0,1)^{\top},(0,0,1,1,0)^{\top} (1,0,0,0,1),(0,0,1,1,0)

因此, a ⊤ x ≤ β a^{\top} x \leq \beta axβ 定义一个维度为1的面

备注6:不同的背包约束 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 可能导致相同的背包集/多面体,例如,如果 a ∈ Z + n a∈\mathbb{Z}^n_+ aZ+n β ∈ Z + β∈\mathbb{Z}_+ βZ+ 有一个最大公约数至少为2。有时也可以在不改变背包集/多面体的情况下增加a的分量或减少 β β β ,例如 x 1 + 2 x 2 ≤ 2 x_1 + 2 x_2≤2 x1+2x22 可以替换为 2 x 1 + 2 x 2 ≤ 2 2x_1 + 2x_2≤2 2x1+2x22 ,然后再替换为 x 1 + x 2 ≤ 1 x_1 + x_2≤1 x1+x21

在下面,我们研究了定义 P a , β P^{a,\beta} Pa,β 不等式的一般facet的性质。

不等式 c ⊤ x ≤ γ c^{\top} x \leq \gamma cxγ, 如果 γ = 0 \gamma = 0 γ=0 则称为齐次,否则称为非齐次。

引理7:(Hammer 等人,1975,命题4) 设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ 。如果 P a , β P^{a,\beta} Pa,β 是全维的,唯一的齐次面定义不等式 P a , β P^{a,\beta} Pa,β 是平凡不等式 − x i ≤ 0 , i ∈ [ n ] −x_i≤0,i∈[n] xi0,i[n] 的正倍数。

引理8:(Hammer 等人, 1975,命题5) 设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ ,使得 P a , β P^{a,\beta} Pa,β 为全维背包多面体,设 c ⊤ x ≤ γ c^{\top} x \leq \gamma cxγ 定义 P a , β P^{a,\beta} Pa,β 的一个面, γ ≠ 0 \gamma \neq 0 γ=0。那么对于所有的 i ∈ [ n ] i∈[n] i[n],有 γ > 0 γ > 0 γ>0 0 ≤ c i ≤ γ 0≤c_i≤γ 0ciγ
与引理3(a)类似,引理8使我们能够证明背包多面体的单调性结果。

推论9:设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ 。背包多面体 P a , β P^{a,\beta} Pa,β 是下单调的,即如果 x ∈ P a , β x∈P^{a,\beta} xPa,β 0 ≤ y ≤ x 0≤y≤x 0yx ,则 y ∈ P a , β y∈P^{a,\beta} yPa,β

请注意引理3、7 和 8 的结果也可以从独立系统的一般陈述(或等价的集合覆盖问题)得到,参见例如Conforti和Laurent (1988), Balas和Ng (1989), Laurent(1989)和Sassano(1989)。

理解具有积分系数的一般背包多面体的一般方法是通过所谓的主背包多面体。主0/1背包多形体 P M P_M PM 是与背包不等式相关的背包多形体

在这里插入图片描述
K i : = 1 + ⌊ β i ⌋ K_i:=1+\left\lfloor\frac{\beta}{i}\right\rfloor Ki:=1+iβ ,即 K i K_i Ki 是大于 β i \frac{β}{i} iβ 的最小整数,二元变量 x j i x^i_j xji 对应系数为 i i i 的第 j j j 个变量。以下结果是由Hammer和Peled(1982)所述的Hammer等人(1977)得出的。

观察任何背包不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ a ∈ Z + n , β ∈ Z + a \in \mathbb{Z}_{+}^n, \beta \in \mathbb{Z}_{+} aZ+n,βZ+ 表示为公式(2):对于 i ∈ [ β ] i∈[β] i[β] ,假设 A i : = { ℓ ∈ [ n ] : a ℓ = i } A_i:=\left\{\ell \in[n]: a_{\ell}=i\right\} Ai:={[n]:a=i} 和 定义 I a , β : = { ( i , j ) : i ∈ [ β ] , j = 1 , . . . , ∣ A i ∣ I_{a, \beta}:=\{(i, j): i \in[\beta], \quad j=1,...,|A_i| Ia,β:={(i,j):i[β],j=1,...,Ai

那么背包多面体 P a , β P^{a,\beta} Pa,β P M P_M PM 的面同构, P M P_M PM 的面定义为 x j i = 0 , ( i , j ) ∉ I a , β x_j^i=0,(i, j) \notin I_{a, \beta} xji=0,(i,j)/Ia,β

因此,如果我们得到 P M P_M PM 的完整线性描述,我们就可以直接推导出任意背包多面体的完整线性描述。Hammer和Peled(1982)给出了 β ≤ 7 β≤7 β7 P M P_M PM 的完整描述。

然而,一般来说,一个完整的线性描述是不可用的,人们需要考虑一个具体的背包上的知识,以找到定义不等式的面。

相对于背负式多面体 P P P ,该多面体对应于标准LP松弛

在这里插入图片描述

很好理解。假设1成立,背包不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 和下界约束 x i ≥ 0 , i ∈ [ n ] x_i≥0 ,i∈[n] xi0,i[n] 始终是facet定义。当 a i < β a_i < β ai<β 时,上限约束 x i ≤ 1 x_i≤1 xi1 为facet定义。因此, P L P P_{LP} PLP 至少有 n + 1 n+1 n+1 个facet,且至多有 2 n + 1 2n+1 2n+1 个facet 。

此外, P L P P_{LP} PLP 顶点的完整表征是可用的,它将Dantzig(1957)对最优LP解的表征转化为多面体语言。

命题10:让假设1保持不变。向量 x ˉ ∈ [ 0 , 1 ] n \bar{x} \in[0,1]^n xˉ[0,1]n P L P P_{LP} PLP 的顶点,当且仅当存在不相交集 I 0 , I 1 ⊆ [ n ] I_0, I_1 \subseteq[n] I0,I1[n] ∣ I 0 ∪ I 1 ∣ ≥ n − 1 \left|I_0 \cup I_1\right| \geq n-1 I0I1n1 ,使得对于 j ∈ [ n ] \ ( I 0 ∪ I 1 ) j \in[n] \backslash\left(I_0 \cup I_1\right) j[n]\(I0I1),有 a ( I 1 ) ≤ β , a j > β − a ( I 1 ) a(I_1)≤β,a_j > β−a(I_1) a(I1)β,aj>βa(I1) ,,并且对于每一个 i ∈ [ n ] i∈[n] i[n] 我们有:

在这里插入图片描述

假设11:背包不等式的系数非递减排序,即 a 1 ≤ a 2 ≤ ⋅ ⋅ ⋅ ≤ a n a_1≤a_2≤···≤a_n a1a2⋅⋅⋅an

在本节的最后,我们提出了一些研究问题;据我们所知,这些问题目前是开放的,但部分答案将在之后讨论。

(Q1) (渐近地)存在多少不同的背包多形体(w.r.t.仿射或0/1等价)?

(Q2) 一般来说,确定定义不等式的背包 a ⊤ x ≤ β a^{\top} x\le \beta axβ 定义了 P P P 的非空面,因为我们需要确定相等背包是否有解。是否存在有效的充分条件来保证 a a a 的二元可行性 a ⊤ x = β a^{\top} x=\beta ax=β 还是 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 定义facet?

(Q3) 是否有一个定义不等式的背包上的可能操作的特征,使得背包集/多面体不改变,参见注释6?它试图描述集合的特征

在这里插入图片描述

关于问题(Q3), Bradley等人(1974)证明了线性不等式系统的解

在这里插入图片描述
它们是否定义了与a相同的背负式多面体 a ⊤ x ≤ β a^{\top} x \leq \beta axβ ,假设11成立。因此,对于一个具体的背包实例,(Q3)可以回答,而对背包不等式的一般运算/操作产生相同的背包多面体是未知的。

3.2 Covers 覆盖不等式

背包多形体最重要的一类有效不等式是由覆盖不等式给出的,本节将讨论覆盖不等式。

定义12:设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n and β ∈ R + \beta \in \mathbb{R}_{+} βR+ 。如果 a ( C ) > β a(C)>\beta a(C)>β ,则集合 C ⊆ [ n ] C \subseteq[n] C[n] 被称为 P a , β P^{a, \beta} Pa,β 的覆盖 。如果它不包含作为覆盖的适当子集,即 C \ { i } C \backslash\{i\} C\{i} 不是任何 i ∈ C i∈C iC 的覆盖,则称为最小覆盖不等式。对应的(最小)覆盖不等式为:

在这里插入图片描述

覆盖不等式对 P a , β P^{a, \beta} Pa,β 明显有效,参见,例如,Crowder 等人(1983)。此外,盒子约束和所有最小覆盖不等式的集合为背包集提供了一个整数公式,因为对于一个点 x ~ ∈ { 0 , 1 } n \tilde{x} \in\{0,1\}^n x~{0,1}n a ⊤ x ~ > β a^{\top} \tilde{x}>\beta ax~>β ,覆盖 { i ∈ [ n ] : x ~ i = 1 } \left\{i \in[n]: \tilde{x}_i=1\right\} {i[n]:x~i=1} 包含一个最小覆盖,其覆盖不等式被 x ~ \tilde{x} x~ 违反。

示例13:(继续示例5)再次考虑背包不等式(1) 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 。对应的最小覆盖为:

在这里插入图片描述

集合 C = 1 , 2 , 3 , 4 C ={1,2,3,4} C=1,2,3,4 是一个覆盖,但不是 C 1 ⊊ C C_1 \subsetneq C C1C

通过分析背包多面体的最小覆盖不等式,可以推导出背包多面体的面。为此,我们将背包多面体限制为包含在封面中的变量,并得到以下结果,这已由Padberg(1975)隐式地表明。

命题14:设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ 满足 P a , β P^{a, \beta} Pa,β 是全维的,设 C C C 是最小覆盖 w . r . t w.r.t w.r.t 。背包不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 。然后 x ( C ) ≤ ∣ C ∣ − 1 x(C) \leq|C|-1 x(C)C1 定义了背包多面体的一个面:

在这里插入图片描述
它定义了一个 P a , β P^{a, \beta} Pa,β 的面,其维度数至少为 ∣ C ∣ − 1 |C|−1 C1

定义15:设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ ,设假设11成立。然后一个覆盖 C = { i 1 , … , i k } C = \{i_1,…, i_k\} C={i1ik} i 1 < i 2 < ⋯ < i k i_1<i_2<\cdots<i_k i1<i2<<ik ,如果 C C C 是一个最小覆盖,则称为 S t r o n g Strong Strong

在这里插入图片描述
对于每个 j ∈ [ i k − 1 ] \ C j \in\left[i_k-1\right] \backslash C j[ik1]\C 。而且,扩展a(不一定强)覆盖 C = { i 1 , … , i k } C = \{i_1,…, i_k\} C={i1ik} 是集合 E ( C ) : = C ∪ { i k + 1 , … , n } E(C):=C \cup\left\{i_k+1, \ldots, n\right\} E(C):=C{ik+1,,n} 。由于权重 a i a_i ai 是按非递减排序的,所以可拓不等式:

在这里插入图片描述
对P有效

示例16:(继续例子5)注意,假设11在这个例子中成立。较强的覆盖层为 C 1 , C 5 , C 6 C_1, C_5, C_6 C1,C5,C6 C 7 C_7 C7 。为了查看 C 5 = 2 , 5 C_5 ={2,5} C5=2,5 ,我们检查定义15:

在这里插入图片描述
对于每个 j ∈ 1 , 3 , 4 j∈{1,3,4} j1,3,4 ,因为 a 4 = 7 = max ⁡ { a ℓ : ℓ ∈ [ i k − 1 ] \ C 5 } a_4=7=\max \left\{a_{\ell}: \ell \in\left[i_k-1\right] \backslash C_5\right\} a4=7=max{a:[ik1]\C5}

然而,最小的覆盖 C 3 C_3 C3 并不强大,因为从覆盖上删除第4项并在封面上添加第2项会产生一个不可行的项目集合。

关于强覆盖(扩展不等式)的进一步细节由Balas(1975)提供。

此外,对于可拓不等式定义 P a , β P^{a, \beta} Pa,β 面的情况,存在一个完整的描述。

定理17:(Balas 1975)设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ ,设假设11成立。设 E ⊆ [ n ] E⊆[n] E[n] 包含至少两个元素,设 k ∈ [ n ] k∈[n] k[n] 。不等式 x ( E ) ≤ k − 1 x(E)≤k−1 x(E)k1 定义了 P a , β P^{a, \beta} Pa,β 的一个面当且仅当它是强覆盖的扩展不等式 C = i 1 , … , i k , i 1 < i 2 < ⋅ ⋅ ⋅ < i k C = {i_1,…, i_k}, i_1 < i_2 <···< i_k C=i1iki1<i2<⋅⋅⋅<ik 和:

在这里插入图片描述
注意,定理17中的 x ( E ) ≤ k − 1 x(E)≤k−1 x(E)k1 在两种情况下都是有效的(要么它定义了一个面,要么它来自一个强覆盖)。Laurent(1989,命题3.11和3.14)给出了一个等价于定理17在独立系统方面的描述。

示例18:(继续例5)如例16所示,运行示例的强覆盖是 C 1 , C 5 , C 6 , C 7 {C_1, C_5, C_6, C_7} C1,C5,C6,C7 。它们的扩展是:

在这里插入图片描述
条件(4)对于17个结果:

在这里插入图片描述
因此,所有不等式 x ( E ( C i ) ) ≤ ∣ C i ∣ − 1 , i ∈ 1 , 5 , 6 , 7 x(E(C_i))≤|C_i |−1 , i∈{1,5,6,7} x(E(Ci))Ci1,i1,5,6,7 ,定义了facet,这些是定义 x ( E ) ≤ k − 1 x(E)≤k−1 x(E)k1 形式的不等式的唯一facet。

(Q4) 我们能否得到一个背包的特征,其中覆盖不等式提供了 P a , β P^{a, \beta} Pa,β 的完整线性描述?

(Q5) 使用(最小)覆盖不等式的公式有多强(例如,就LP-gap而言)?


四、Binary formulations based on strong covers 基于强覆盖的二元公式

分离 K = K a , β K=K^{a, \beta} K=Ka,β { 0 , 1 } n \ K \{0,1\}^n \backslash K {0,1}n\K 的一组线性不等式 A x ≤ b Ax≤b Axb 称为 K K K 的二元公式,即 K = { x ∈ { 0 , 1 } n : A x ≤ b } K = \{x∈\{0,1\}^n: Ax≤b\} K={x{0,1}n:Axb} 。利用强覆盖(扩展)的概念,我们可以找到 K K K 的最小大小的二进制公式,其所有不等式的左边系数都是0或1,即所谓的0/1-二进制公式。

这样的公式是有趣的,因为它们通常在数值上比原始不等式更稳定,如果后者包含较大的数字。此外,这些不等式有一个简单的解释,因为它们编码了变量子集的基数限制。以下结果是由于Glover(1973),由Wolsey(1975)所述。

定理19: (Glover 1973)设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ 。假设11成立,设 S \mathfrak{S} S K a , β K^{a, \beta} Ka,β 的所有强覆盖的集合。

在这里插入图片描述
K a , β K^{a, \beta} Ka,β 最小值的0/1-二进制公式。也就是说,在左边没有更小的0/1系数不等式集可以保证一个二元向量包含在 K a , β K^{a, \beta} Ka,β 中。

请注意,Balas(1975)以另一种方式引用了Glover的结果,他将定理19中的强覆盖限制在集合上

在这里插入图片描述
因此,如果 S ‾ ⊊ S \overline{\mathfrak{S}} \subsetneq \mathfrak{S} SS ,定理19中Balas版本的二元公式的大小比Wolsey版本的小,这与公式的最小值是矛盾的。但是,不可能出现这种差异。

引理20:设 C C C C ‾ \overline{C} C 分别为由 a ⊤ x ≤ β a^{\top} x \leq \beta axβ ∣ C ∣ = ∣ C ˉ ∣ |C|=|\bar{C}| C=Cˉ 。那么 E ( C ) ⊆ E ( C ˉ ) E(C) \subseteq E(\bar{C}) E(C)E(Cˉ) E ( C ˉ ) ⊆ E ( C ) E(\bar{C}) \subseteq E(C) E(Cˉ)E(C) 都不成立。

示例21:(继续例5)根据例16,所有强覆盖的集合是 S = { C 1 , C 5 , C 6 , C 7 } \mathfrak{S}=\left\{C_1, C_5, C_6, C_7\right\} S={C1,C5,C6,C7} ,扩展在例18中给出。因此, K a , β K^{a, \beta} Ka,β 的最小尺寸0/1-二进制公式为

在这里插入图片描述
一般来说,定理19所提供的 K K K 的0/1-二进制公式的最小大小可以是 n n n 的指数级,精确计算可能比较复杂。因此,我们提出了一种简单的技术来生成 ∣ S ∣ |\mathfrak{S}| S 的下界。假设 A : = { a i : i ∈ [ n ] } A:=\left\{a_i: i \in[n]\right\} A:={ai:i[n]} 是背包不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 。对每个 α ∈ A α∈A αA , 我们定义 A α : = { i ∈ [ n ] : a i = α } A_\alpha:=\left\{i \in[n]: a_i=\alpha\right\} Aα:={i[n]:ai=α} ,即 A α A_α Aα 包含背包约束中系数为 α α α 的变量的所有指标。此外,定义 ( A α n α ) \left(\begin{array}{l}A_\alpha \\ n_\alpha\end{array}\right) (Aαnα) 是包含 A α A_α Aα 的所有子集的集合,其基数为 n α n_α nα

命题22:设假设11成立, a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ C C C K a , β K^{a,\beta} Ka,β 的强覆盖, n α : = ∣ C ∩ A α ∣ , α ∈ A n_\alpha:=\left|C \cap A_\alpha\right|, \alpha \in A nα:=CAα,αA ,并定义 α ˉ : = max ⁡ { α ∈ A : n α > 0 } \bar{\alpha}:=\max \left\{\alpha \in A: n_\alpha>0\right\} αˉ:=max{αA:nα>0} 。那么,对于每一个 C α ′ ∈ ( A α n α ) , α ∈ A \ { α ˉ } C_\alpha^{\prime} \in\left(\begin{array}{l}A_\alpha \\ n_\alpha\end{array}\right), \alpha \in A \backslash\{\bar{\alpha}\} Cα(Aαnα),αA\{αˉ} 都有:

在这里插入图片描述
K K K 的强大覆盖。

示例23:考虑背包不等式 ∑ i = 1 n x i + 2 ∑ i = n + 1 2 n x i ≤ n \sum_{i=1}^n x_i+2 \sum_{i=n+1}^{2 n} x_i \leq n i=1nxi+2i=n+12nxin 写成 a ⊤ x ≤ n a^{\top} x \leq n axn, 当 n n n 为偶数。令 I k : = [ k ] ∪ { n + 1 , … , n + n − k − 1 2 } I_k:=[k] \cup\left\{n+1, \ldots, n+\frac{n-k-1}{2}\right\} Ik:=[k]{n+1,,n+2nk1} , 对于奇数 k ∈ [ n ] k∈[n] k[n] 。这时,一个 a ⊤ χ I k = n − 1 a^{\top} \chi^{I_k}=n-1 aχIk=n1 表示 χ I k ∈ K \chi^{I_k} \in K χIkK 。此外,将权值为2的任何元素添加到 I k I_k Ik ,将权值增加到 n + 1 n + 1 n+1 。因此,扩展集不在 K K K 中。因此,集合 C k : = I k ∪ { n + n − k + 1 2 } C_k:=I_k \cup\left\{n+\frac{n-k+1}{2}\right\} Ck:=Ik{n+2nk+1} 是一个最小覆盖,它是强覆盖,因为 [ n + n − k + 1 2 ] \ C k \left[n+\frac{n-k+1}{2}\right] \backslash C_k [n+2nk+1]\Ck 中的每个元素的权值都是1。因此,22号命题表明,我们可以生成 ( n k ) \left(\begin{array}{l}n \\ k\end{array}\right) (nk) 的许多强有力的封面,它们是两两不同的。此外,请注意由 C k C_k Ck C k ′ C_{k'} Ck 对于不同的 k k k k ′ k' k 不能重合,因为它们在 [ n ] [n] [n] 中的元素数量不同。因此,K的任何0/1-二进制公式至少包含

在这里插入图片描述

不等式。因此,与系数为 0 , ± 1 , ± 2 {0,±1,±2} 0±1±2 的公式相比,对0/1-系数不等式的限制会导致必要约束的数量急剧增加。

示例24:考虑背包不等式 ∑ i = 1 2 n 2 ⌈ i / 2 ⌉ − 1 x i ≤ 2 n − 1 \sum_{i=1}^{2 n} 2^{\lceil i / 2\rceil-1} x_i \leq 2^n-1 i=12n2i/21xi2n1 ,记为 a ⊤ x ≤ β a^{\top} x \leq \beta axβ ,与Kaibel和Loos(2011)提出的所谓orbisack有关。设 C = { i ∈ [ 2 n ] : i C=\{i \in[2 n]: i C={i[2n]:i odd or i = 2 } i=2\} i=2} 。因为一个 a ⊤ χ C = 2 n a^{\top} \chi^C=2^n aχC=2n , C C C 是最小覆盖,因为所有系数都大于或等于1。

此外,$C是强的:当 n ≤ 2 n≤2 n2 时,这是一个case区别。否则,如果 n > 2 n > 2 n>2 ,我们有 2 n − 2 2n - 2 2n2 是我们估计的非 C C C 中系数最大的指标

在这里插入图片描述
因此,C是一个强大的覆盖。

因为C包含所有权值为1的元素,一个权值为 2 n − 1 2^{n−1} 2n1 的元素,以及 2 i 2^i 2i for i ∈ { 1 , … , n − 2 } i \in\{1, \ldots, n-2\} i{1,,n2} ,命题22意味着我们可以从 C C C 生成 2 n − 2 2^{n−2} 2n2 个强覆盖。这表明我们需要在 K K K 的任何0/1二进制公式中有指数级的不等式。事实上,我们可以证明强覆盖的确切数量是 2 n − 1 2^{n−1} 2n1

因此,通过猜测单个集合 C C C ,我们可以生成所有强覆盖的一半。

(Q6) 定理19通过要求 x ∈ 0 , 1 n x∈{0,1}^n x0,1n 隐式地使用边界。公式 A x ≤ b Ax≤b Axb ,系数为 0 , ± 1 {0,±1} 0±1 ,使得 K = { x ∈ Z n : A x ≤ b } K=\left\{x \in \mathbb{Z}^n: A x \leq b\right\} K={xZn:Axb} 的最小尺寸是多少?上界是 ∣ S ∣ + 2 n |\mathfrak{S}|+2 n S+2n 通过显式地加上平凡的不等式。

(Q7) 强覆盖数量的下界依赖于 K K K 通过背包不等式 a ⊤ x ≤ β a^{\top} x \leq \beta axβ 。但是,如果 a a a 中的所有系数都不同,则不适用边界技术。在这种情况下,是否有一种简单的机制来推导强覆盖数的边界?


五、Lifting 提升

在本节中,我们考虑提升,一个重要的概念,以加强背包多面体 P : = P a , β P:=P^{a, \beta} P:=Pa,β a ∈ R + n a \in \mathbb{R}_{+}^n aR+n β ∈ R + \beta \in \mathbb{R}_{+} βR+ 的有效不等式。(上)提升需要一个不平等 ∑ i ∈ S π i x i ≤ π 0 \sum_{i \in S} \pi_i x_i \leq \pi_0 iSπixiπ0 ,具有 S ⊆ [ n ] S \subseteq[n] S[n] ,对于受限背包集 K S : = { x ∈ K : x i = 0 , i ∈ [ n ] \ S } K_S:=\{x \in K:x_i=0,i \in [n] \backslash S\} KS:={xK:xi=0,i[n]\S} 有效,并通过将变量 x i , i ∉ S x_i, i \notin S xi,i/S 的系数从0提升到适当的提升系数,将其转化为 K : = K a , β K:=K^{a, \beta} K:=Ka,β 的有效不等式。目标是获得一个更强或更均匀的定义 P P P 不等式的方面。

也可以减小包含在初始不等式中的变量系数,即所谓的 down-lifing 。下面我们详细讨论了lift,经常用lift代替up-lifting。 down-lifing 将在5.2节中讨论。

我们首先注意到“trivial lifting”总是可能的,即S外系数为0的不等式对 K K K 是有效的:

引理25:设 a ∈ R + n a \in \mathbb{R}_{+}^n aR+n and β ∈ R + \beta \in \mathbb{R}_{+} βR+. If ∑ i ∈ S π i x i ≤ π 0 \sum_{i \in S} \pi_i x_i \leq \pi_0 iSπixiπ0 对于 K S K_S KS 有效,那么它对于 K K K 也是有效的。

5.1 Sequential up-lifting

顺序提升修复了 [ n ] \ S [n]\backslash S [n]\S 的顺序。通过迭代求解一系列优化问题,计算出可能的最强提升系数,使不等式对相应的受限背包集仍然有效。在每次迭代中,所有先前计算的系数都被考虑在内。

这个过程将首先描述一般有效不等式。然后应用于最小覆盖不等式。

可以为任意有效不等式定义 Sequential up-lifting 过程

在这里插入图片描述

对于 π ≥ 0 π≥0 π0 K S K_S KS。那么对于 k ∈ [ n ] \ S k∈[n]\backslash S k[n]\S 计算 π k π_k πk ,使:

在这里插入图片描述
等于 K S ∪ { k } K_{S \cup\{k\}} KS{k} 。考虑升降函数 Φ S : R + → R \Phi_S: \mathbb{R}_{+} \rightarrow \mathbb{R} ΦS:R+R 定义为

在这里插入图片描述

定理26:(Padberg 1975)设 a ∈ R + n , β ∈ R + a \in \mathbb{R}_{+}^n, \beta \in \mathbb{R}_{+} aR+n,βR+ ,令(6)对 K S K_S KS 成立,如果 π k ≤ Φ S ( a k ) \pi_k \leq \Phi_S\left(a_k\right) πkΦS(ak) ,不等式(7)对 K S ∪ { k } K_{S \cup\{k\}} KS{k} 成立。此外,如果 d i m ( P ) = n dim(P) = n dim(P)=n π k = Φ S ( a k ) \pi_k=\Phi_S\left(a_k\right) πk=ΦS(ak) ,且(6)定义了conv( K S K_S KS)的一个维数为t的面,(7)定义了conv( K S ∪ { k } K_{S \cup\{k\}} KS{k})维数至少为 t + 1 t + 1 t+1 的面。

计算 k ∈ [ n ] \ S k \in[n] \backslash S k[n]\S 的提升系数 π k = Φ S ( a k ) \pi_k=\Phi_S\left(a_k\right) πk=ΦS(ak) 后,可以将 S S S 更新为 S ∪ k S∪{k} Sk ,并计算 [ n ] \ S [n]\backslash S [n]\S 的新元素的下一个系数(只要 S ≠ [ n ] S \neq[n] S=[n] )。通过迭代这个过程,其他提升系数可以一个接一个地计算,从而得出的不等式对 K K K 有效。一般来说,提升函数 Φ S \Phi_S ΦS 会随着更多的变量被提升而减小。因此,提升系数取决于提升变量的顺序。这种举升过程称为 sequential lifting 。

请注意,由于 π ≥ 0 , Φ S ( u ) ≤ π 0 \pi \geq 0, \Phi_S(u) \leq \pi_0 π0,ΦS(u)π0 成立。而且,由于(6)对 K S K_S KS 有效,且(8)中考虑的每一个可行解 x x x u ≥ 0 u≥0 u0 时都包含在 K S K_S KS 中,则 Φ S ( u ) ≥ 0 \Phi_S(u) \geq 0 ΦS(u)0

因此,在 u u u 中求 Φ S \Phi_S ΦS 等于求解一个0/1背包问题,其客观值包含在 [ 0 , π 0 ] [0,π_0] [0π0] 中。迭代应用定理26,得到的不等式对较大的集合有效,因此每个中间目标值位于区间 [ 0 , π 0 ] [0,π_0] [0π0] 。此外,如果原始不等式(6)具有积分系数,则计算得到的 lifting 系数也是积分系数。

因此,若 π π π π 0 π_0 π0为积分,则动态规划可在 O ( n π 0 ) \mathcal{O}\left(n \pi_0\right) O(nπ0) 时间内求出 Φ S ( u ) \Phi_S(u) ΦS(u) 。因此,计算所有提升系数的总运行时间为 O ( n 2 π 0 ) \mathcal{O}\left(n^2 \pi_0\right) O(n2π0)

注意,如果a和β是积分, Φ S ( u ) \Phi_S(u) ΦS(u) 可以在 O ( n β ) O(n β) O(nβ) 时间内计算出来,因为 u ≥ 0 u≥0 u0 and,所以右边总是最大 β β β 。然而,常常有 π 0 < β \pi_0<\beta π0<β, 所以第一种方式更快;请参见下一节。

一类重要的有效不等式是通过解除最小覆盖不等式得到的,即起始不等式(6)为 x ( C ) ≤ ∣ C ∣ − 1 x(C) \leq|C|-1 x(C)C1 ;其中 S = C S=C S=C π 0 = ∣ C ∣ − 1 \pi_0=|C|-1 π0=C1。考虑某个 i 1 , … , i k i_1, \ldots, i_k i1,,ik of [ n ] \ C [n] \backslash C [n]\C ,并定义

在这里插入图片描述

对于 r ∈ [ k ] r∈[k] r[k] ,设 α i = 1 , i ∈ C α_i = 1 ,i∈C αi=1,iC ,迭代 α i r = Φ C r ( a i r ) \alpha_{i_r}=\Phi_{C_r}\left(a_{i_r}\right) αir=ΦCr(air) for r ∈ [ k ] r \in[k] r[k] 。由此产生的不等式

在这里插入图片描述
称为简单解除覆盖不等式(LCI);“简单”指的是我们在这里只是在抬升。注意(9)中的所有系数都是构造积分。

回顾命题14, x ( C ) ≤ ∣ C ∣ − 1 x(C) \leq|C|-1 x(C)C1 定义了conv( K C K_C KC)的一个面,它仿射等价于 P C P_C PC。因此,根据定理26,得到的简单 L C I LCI LCI K K K 和有效定义 P P P = conv( K K K)的一个面。然而,由于举升系数取决于 [ n ] \ C [n]\backslash C [n]\C 的顺序,举升一个盖C理论上可以产生 ( n − ∣ C ∣ ) ! (n-|C|) ! (nC)! (不一定是不同的) P P P 的方面。

采用上述动态规划方法,可在 O ( n ∣ C ∣ 2 ) \mathcal{O}\left(n|C|^2\right) O(nC2) 时间内计算出所有升力系数。事实上,这甚至可以改进到Zemel(1989)所示的 O ( n ∣ C ∣ ) \mathcal{O}(n|C|) O(nC) ,即计算最小覆盖不等式的所有提升系数与计算单个提升系数具有相同的复杂性。此外,Kaparis和Letchford(2010)指出,结合Balas和Zemel(1978)和Gu等人(2000)的思想,大部分 lifting 系数可以在 O ( n + ∣ C ∣ log ⁡ ∣ C ∣ ) \mathcal{O}(n+|C| \log |C|) O(n+ClogC) 时间内计算出来。

示例27:(继续示例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 。与最小覆盖 C = { 2 , 3 , 4 } C =\{2,3,4\} C={2,3,4} 相关的最小覆盖不等式由 x 2 + x 3 + x 4 ≤ 2 x_2 + x_3 + x_4≤2 x2+x3+x42 给出。为了消除这个不等式,我们对 [ n ] \ C = { 1 , 5 } [n] \backslash C=\{1,5\} [n]\C={1,5} 排序有两种选择。

如果我们使用排序 ( i 1 , i 2 ) = ( 1 , 5 ) (i_1, i_2) =(1,5) (i1,i2)=(1,5),即首先将变量x1和x5提升到覆盖不等式中,我们有 C 1 = C = { 2 , 3 , 4 } C_1 = C =\{2,3,4\} C1=C={2,3,4} C 2 = { 1 , 2 , 3 , 4 } C_2 =\{1,2,3,4\} C2={1,2,3,4} 。在第一步中, x 1 x_1 x1 的lifting系数 α i 1 = Φ C 1 ( a i 1 ) \alpha_{i_1}=\Phi_{C_1}\left(a_{i_1}\right) αi1=ΦC1(ai1)

在这里插入图片描述

我们得到不等式 x 1 + x 2 + x 3 + x 4 ≤ 1 x_1 + x_2 + x_3 + x_4≤1 x1+x2+x3+x41 x 5 x_5 x5 的lifting系数 α i 2 = Φ C 2 ( a i 2 ) \alpha_{i_2}=\Phi_{C_2}\left(a_{i_2}\right) αi2=ΦC2(ai2)

在这里插入图片描述

这就得到了简单的LCI x 1 + x 2 + x 3 + x 4 + x 5 ≤ 2 x_1 + x_2 + x_3 + x_4 + x_5≤2 x1+x2+x3+x4+x52

使用排序 ( i 1 , i 2 ) = ( 5 , 1 ) (i_1, i_2) =(5,1) (i1,i2)=(5,1) 得到:

在这里插入图片描述

这就产生了不同的简单LCI x 2 + x 3 + x 4 + 2 x 5 ≤ 2 x_2 + x_3 + x_4 + 2x_5≤2 x2+x3+x4+2x52

为了加快提升系数的计算速度,Balas(1975)提出了一种较弱的提升形式,它在线性时间内运行,并且不依赖于求解优化问题。

折衷之处在于Balas的简单LCI不一定定义P的面。

定理28:(Balas 1975)设 a ∈ R + n , β ∈ R + a \in \mathbb{R}_{+}^n, \beta \in \mathbb{R}_{+} aR+n,βR+ ,设假设11成立,设 C C C P P P 的最小覆盖, E ( C ) E(C) E(C) C C C 的扩展, C h C_h Ch C C C 的最后h个元素的集合, h ∈ { 1 , … , ∣ C ∣ } h \in\{1, \ldots,|C|\} h{1,,C} 。考虑分区 N 0 , N 1 , … , N q , q = ∣ C ∣ − 1 N_0, N_1, \ldots, N_q, q=|C|-1 N0,N1,,Nq,q=C1, of [ n ] [n] [n]

在这里插入图片描述
然后定义:

在这里插入图片描述

然后是不等式

在这里插入图片描述

对于 P P P 是有效的,如果

在这里插入图片描述

然后(10)定义 P P P 的一个面

最近,Letchford和Souli(2019)提出了一种改进的提升程序,如果覆盖 C C C 不是最小的,也可以应用该程序。它生成的不等式不弱于(10),并且在某些情况下可以更强,即使 C C C 是最小的。特别地,Letchford 's和Souli 's方法可以用包含在C中的指标修改变量的系数,并且对于最小覆盖C,它产生一个形式为(9)的简单LCI。

推论29:(Balas and Zemel 1978)设 C C C 为最小覆盖。如果条件(11)成立,则每个提升序列都会导致相同的顺序提升的最小覆盖不等式。这个不等式由

在这里插入图片描述
其中 π i π_i πi 定义为定理28。此外,它是定义 P P P 的不等式的唯一方面,对于所有 i ∈ C i∈C iC ,系数为1,右侧为 ∣ C ∣ − 1 |C|−1 C1

因此,在满足条件(11)的情况下,通过依次提升 C C C 的最小覆盖不等式,恰好可以找到 P P P 的一个方面。特别是,定义这个方面的不等式可以通过Balas过程找到,并且不需要依赖于顺序提升方法。

示例30:(继续示例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 。对于最小覆盖 C = { 1 , 2 , 3 } C =\{1,2,3\} C={1,2,3} 及其扩展 E ( C ) = { 1 , 2 , 3 , 4 , 5 } E(C) =\{1,2,3,4,5\} E(C)={1,2,3,4,5} ,定理28中定义的六个集合 C 1 , C 2 , C 3 C_1, C_2, C_3 C1,C2,C3 and N 2 , N 1 , N 0 N_2, N_1, N_0 N2,N1,N0

在这里插入图片描述

因此,我们得到 π 4 = π 5 = 1 π_4 = π_5 = 1 π4=π5=1 和有效不等式

在这里插入图片描述

由于 N 0 = N 2 = ∅ N_0 = N_2 =∅ N0=N2= ,只需要在 h = 1 h = 1 h=1 时检查条件(11)。由于 a 1 ≤ 13 − a j a_1≤13−a_j a113aj 对于所有 j ∈ [ 5 ] j∈[5] j[5] 都为真,有效不等式(12)定义了 P P P 的一个面。此外,根据推论29,它是 P P P 的唯一面,可以通过顺序提升从 C C C 得到。实际上,(12)是覆盖 C = { 1 , 2 , 3 } C =\{1,2,3\} C={1,2,3} 的可拓不等式。

对于封面 C = { 2 , 3 , 4 } C =\{2,3,4\} C={2,3,4} 及其扩展 E ( C ) = { 2 , 3 , 4 , 5 } E(C) =\{2,3,4,5\} E(C)={2,3,4,5} ,六个集合 C 1 , C 2 , C 3 C_1, C_2, C_3 C1,C2,C3 N 2 , N 1 , N 0 N_2, N_1, N_0 N2,N1,N0

在这里插入图片描述

这就得到了有效的不等式

在这里插入图片描述

这不是 P P P 的面定义,因为它是由例27中的取消的不等式主导的。因此,如果条件(11)不成立(如本例中 h = 0 h = 0 h=0 j = 1 j = 1 j=1 ),定理28不一定会产生顺序提升的facet。

在示例30中,我们已经看到,顺序提升可以产生覆盖的扩展不等式。下面的推论描述了这种情况何时会发生。回想一下,根据定义,强覆盖是最小覆盖。

推论31:假设11成立,假设 C = { i 1 , … … , i k } C = \{i1,……, ik\} C={i1……ik} 是K的覆盖,其中 i 1 < i 2 < ⋅ ⋅ ⋅ < i k i_1 < i_2 <···< i_k i1<i2<⋅⋅⋅<ik K ≥ 2 K≥2 K2 。当且仅当 C C C 是强覆盖且(4)成立时, C C C 的可拓不等式是 P P P 的顺序提升覆盖不等式。

(Q8) 使用解除覆盖不等式的公式有多强(例如,就LP-gap而言)?
(Q9) 有多少被解除的覆盖不平等?
(Q10) 在定理28中,附加假设(11)不仅依赖于背包多面体,而且依赖于它通过背包不等式的表示。是否存在独立于背包不等式的类似特征?

5.2 Sequential down-lifting

上面描述的过程称为向上提升,因为在初始不等式中为零的系数可能会增加。还记得我们假设初始不等式吗? ∑ i ∈ M π i x i ≤ π 0 \sum_{i \in M} \pi_i x_i \leq \pi_0 iMπixiπ0 对于 K S K_S KS 有效,即不在 S S S 中的变量处于其下界。

与此相反,向下提升过程假设初始不等式对集合 K S : = { x ∈ K : x i = 1 , i ∈ [ n ] \ S } K^S:=\left\{x \in K: x_i=1, i \in[n] \backslash S\right\} KS:={xK:xi=1,i[n]\S} 有效 ,即 S S S 之外的变量在其上界。

虽然在上升过程中的初始不等式对整个背包集 K K K 是有效的,参见引理25,但这个不等式在下降的情况下通常对 K K K 无效。

下推的目的是将一个对 K S K_S KS 有效的不等式转化为对 K K K 有效的不等式。

为此,下提升过程确定 i ∈ [ n ] \ S i \in[n] \backslash S i[n]\S 的系数 π i π_i πi ,使

在这里插入图片描述
K K K 有效。为了找到这样的系数,类似于抬升的情况,抬升序列必须固定,而抬升函数 Φ S : R + → R \Phi^S: \mathbb{R}_{+} \rightarrow \mathbb{R} ΦS:R+R

在这里插入图片描述
被认为是。如果 x k , k ∈ [ n ] \ S x_k, k \in[n] \backslash S xk,k[n]\S 是第一个被提升的变量,提升过程将 k k k [ n ] \ S [n]\backslash S [n]\S 中移除,并确定了当 x k x_k xk 固定在其下界时可能发生的初始不等式的最大违规 π k = Φ S ( a k ) \pi_k=\Phi^S\left(a_k\right) πk=ΦS(ak) 。因此,

在这里插入图片描述

等于 K S ∪ { k } K^{S \cup\{k\}} KS{k} 。请注意,这完全类似于抬升的情况,因为抬升函数 Φ S \Phi_S ΦS 可以等效地定义为

在这里插入图片描述

示例32:(继续示例5)设S ={2,3,4,5}。T h ne x(S)≤1对K S有效,但对K无效,因为它截断了可行解(0,1,1,0,0)?。

为了使这个不等式对 K K K 有效,我们通过计算将变量 x 1 x_1 x1 下移

在这里插入图片描述

得到有效不等式 x 1 + x ( S ) ≤ 1 + 1 = 2 x_1 + x(S)≤1 + 1 = 2 x1+x(S)1+1=2

精确计算举升系数 π k π_k πk 可能很耗时,因为它需要解决一个背包问题。或者,也可以不精确地降低初始不等式。类似于向上的情况,我们得到了定理26的类似情况。

定理33:(Nemhauser and Wolsey 1999)设 a ∈ R + n , β ∈ R + a \in \mathbb{R}_{+}^n, \beta \in \mathbb{R}_{+} aR+n,βR+ ,设(6)对 K S K_S KS 有效。如果 { x ∈ K S ∪ { k } : x k = 0 } ≠ ∅ \left\{x \in K^{S \cup\{k\}}: x_k=0\right\} \neq \varnothing {xKS{k}:xk=0}= and π k ≥ Φ S ( a k ) \pi_k \geq \Phi^S\left(a_k\right) πkΦS(ak) ,则:

在这里插入图片描述

K S ∪ { k } K^{S \cup\{k\}} KS{k} 是有效的。此外,如果 π k = Φ S ( a k ) \pi_k=\Phi^S\left(a_k\right) πk=ΦS(ak) 且(6)定义了conv( K S K_S KS)的一个维数为 t t t 的面,则(13)定义了conv( K S ∪ { k } K^{S \cup\{k\}} KS{k})的一个维数至少为 t + 1 t + 1 t+1 的面。

注意,向上和向下提升的过程可以同时应用于一个不等式。

为了说明这一点,考虑向上和向下提升的最小覆盖不等式。设 C C C 是一个极小盖,设 ( C 0 , C 1 ) (C_0, C_1) (C0,C1) C C C 的一个分区,则提升过程的初始不等式为 x ( C 1 ) ≤ ∣ C 1 ∣ − 1 x(C_1)≤|C_1|−1 x(C1)C11 ,这对KC1有效,但对 K C 1 K^{C_1} KC1 无效,如果 C 1 ≠ C C_1 \neq C C1=C。在固定提升序列时, [ n ] \ C [n]\backslash C [n]\C 中的变量依次向上提升, C 0 C_0 C0 中的变量依次向下提升。结果是

在这里插入图片描述
定义 K K K 的一个方面,见Nemhauser和Wolsey(1999)。在文献中,将这类不等式称为解除覆盖不等式(LCI),而 C 0 = ∅ C_0=\varnothing C0= 的不等式仍然是简单的LCI。对于给定的封面,通过上述动态程序可以在多项式时间内找到简单LCIs,计算非简单LCIs是np困难的,参见Chen和Dai(2019)。

5.3 Sequence independent lifting and superadditivity 序列独立提升和超加性

由顺序提升产生的不平等可能取决于顺序,正如我们在例27中看到的那样。另一种方法是同时计算索引在 [ n ] \ S [n]\backslash S [n]\S 中的变量的所有系数,这可能更快。这个概念被称为序列无关或同步提升。由于举升函数的超可加性起着重要作用,这种方法有时也称为超可加举升。请注意,定理28中的Balas提升过程可以被认为是一种同步提升方法,因为它产生了有效的不等式,其系数不依赖于指标 i ∈ [ n ] \ S i∈[n]\backslash S i[n]\S 的任何顺序。

我们首先描述了任意有效不等式在一般情况下的同时提升,然后更具体地描述了覆盖不等式。

提升一般不等式令 π ∈ R + [ n ] 0 \pi \in \mathbb{R}_{+}^{[n]_0} πR+[n]0 ∑ i ∈ S π i x i ≤ π 0 \sum_{i \in S} \pi_i x_i \leq \pi_0 iSπixiπ0 K S : = { x ∈ K : x i = 0 , i ∈ [ n ] \ S } K_S:=\left\{x \in K: x_i=0, i \in[n] \backslash S\right\} KS:={xK:xi=0,i[n]\S} S ⊆ [ n ] S \subseteq[n] S[n] 的有效不等式。同时提升是通过一个函数 Ψ : R + → R \Psi: \mathbb{R}_{+} \rightarrow \mathbb{R} Ψ:R+R 使

在这里插入图片描述

K K K 有效。 Ψ ( u ) = Φ S ( u ) \Psi(u)=\Phi_S(u) Ψ(u)=ΦS(u) 通常不能得到一个有效的不等式,因为集合S会随着更多的变量被引入而增长。因此,如果 i 1 , … , i k i_1,…, i_k i1ik 为顺序提升方法选择的 [ n ] \ S [n]\backslash S [n]\S 的提升顺序,提升函数(8)满足

在这里插入图片描述

为了克服这个问题,我们需要超可加函数。

定义34:函数 f : R → R f: \mathbb{R} \rightarrow \mathbb{R} f:RR D ⊆ R D⊆R DR 上的超可加性,如果对于所有 d 1 , d 2 ∈ D d_1, d_2∈D d1,d2D 满足 f ( d 1 ) + f ( d 2 ) ≤ f ( d 1 + d 2 ) f (d_1) + f (d_2)≤f (d_1 + d_2) f(d1)+f(d2)f(d1+d2)

定理35:(Gu 等人 . 2000;Wolsey 1977)不等式(15)对 K K K 有效,如果函数 Ψ : R → R \Psi: \mathbb{R} \rightarrow \mathbb{R} Ψ:RR 满足:
(a) Ψ ( u ) ≤ Φ S ( u ) \Psi(u) \leq \Phi_S(u) Ψ(u)ΦS(u) for all u ∈ [ 0 , π 0 ] u \in\left[0, \pi_0\right] u[0,π0] and
(b) Ψ \Psi Ψ [ 0 , π 0 ] [0,π_0] [0π0] 上的超加性。

在Marchand等人(2002年,第2.2.2节)、Gu等人(2000年)或Letchford和Souli(2019年)中可以找到超添加剂提升函数的例子。

提升覆盖不等式令 C = i 1 , … , i k C = {i_1,…, i_k} C=i1ik ,其中 i 1 < ⋅ ⋅ ⋅ < i k i_1 <···< i_k i1<⋅⋅⋅<ik 为最小覆盖,假设11成立。设 a j a_j aj { a i 1 , … , a i k } \left\{a_{i_1}, \ldots, a_{i_k}\right\} {ai1,,aik} 的j个最大元素的和,即A j = ?j i=1 aik+1−i,定义 A j = ∑ i = 1 j a i k + 1 − i A_j=\sum_{i=1}^j a_{i_{k+1-i}} Aj=i=1jaik+1i 以及 λ : = ∑ i ∈ C a i − β > 0 \lambda:=\sum_{i \in C} a_i-\beta>0 λ:=iCaiβ>0 。然后是函数

在这里插入图片描述

对于 u ≥ 0 u≥0 u0 Φ C ( u ) \Phi_C(u) ΦC(u) 占主导地位,并且是 R + \mathbb{R}_{+} R+ 上的超加性,见Marchand et al .(2002)。

这就得到了有效的不等式

在这里插入图片描述

示例36:(继续示例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 和最小覆盖 C = { 2 , 3 , 4 } C =\{2,3,4\} C={2,3,4} 。其中 λ = 5 + 6 + 7 − 13 = 5 λ = 5 + 6 + 7−13 = 5 λ=5+6+713=5 , A 0 = 0 A_0 = 0 A0=0, A 1 = 7 A_1 = 7 A1=7, A 2 = 13 A_2 = 13 A2=13, A 3 = 18 A_3 = 18 A3=18,超加性函数 Ψ ( u ) \Psi(u) Ψ(u)

在这里插入图片描述

这就得到 Ψ ( 4 ) = 2 5 \Psi(4)=\frac{2}{5} Ψ(4)=52 Ψ ( 9 ) = 6 5 \Psi(9)=\frac{6}{5} Ψ(9)=56 ,因此是有效的不等式

在这里插入图片描述

它不同于我们在例27中通过顺序提升从 C C C 中获得的定义不等式的面,并且不定义 P P P 的面。然而,情况并非总是如此。对于最小覆盖 C = { 1 , 2 , 3 } C =\{1,2,3\} C={1,2,3} ,上述方法也得到定义不等式的面(12)。

备注37:(注:例27,30,36)注意,对于所有具有两个元素w.r.t. 4 x 1 + 5 x 2 + 6 x 3 + 7 x 4 + 9 x 5 ≤ 13 4x_1 + 5x_2 + 6x_3 + 7x_4 + 9x_5≤13 4x1+5x2+6x3+7x4+9x513 的最小覆盖,即 C = { 2 , 5 } C = \{2,5\} C={2,5} C = { 3 , 5 } C =\{3,5\} C={3,5} C = { 4 , 5 } C =\{4,5\} C={4,5},所有三个提升过程(定理26和28中的顺序提升,以及序列独立提升)对其余三个变量的系数为0。这是因为这三个覆盖都是 E ( C ) = C E(C) = C E(C)=C 的强覆盖,并且很容易看出,这三种情况下的覆盖不等式已经定义了 P P P 的facet。

Easton和Hooker(2008)提出了一个同时将一组变量提升到覆盖不等式的过程。如果盖子和要提升的变量集是排序的,这个过程只需要线性工作。具体来说,对于一个覆盖 C ⊆ [ n ] C⊆[n] C[n] 和一个集合 F ⊆ [ n ] \ C F⊆[n]\backslash C F[n]\C ,他们定义了 ( F , ρ ) (F, ρ) (Fρ) -同时解除的覆盖不等式(SLCI)。

在这里插入图片描述

在这里插入图片描述

如果 ρ ≥ 0 ρ≥0 ρ0 足够小,这是有效的。Easton和Hooker(2008,定理2.1)给出了一个精确的描述。一般来说,slci不是面定义的,而只是增加了由底层覆盖不等式定义的面的维度。但是,Easton和Hooker(2008,定理2.3)给出了 SLCIs 确实是 P C ∪ F = conv ⁡ ( K C ∪ F ) P_{C \cup F}=\operatorname{conv}\left(K_{C \cup F}\right) PCF=conv(KCF) 的facet定义的充分条件。

表1总结了本节和前两节中提出的不同提升方法。对于本表,假设11对背包不等式的系数和最小覆盖 C C C 中的指数都成立,即不需要对变量/系数进行排序。然后,“复杂性”列显示了计算所有提升系数的复杂性,而“facet”列报告了提升方法是否生成了定义不平等的facet,而没有进一步的假设。

5.4 Complete characterization of facets from simple LCIs 从简单LCIs中完全描述facet

Balas和Zemel(1978)给出了facet的完整描述,定义了可以通过连续或同时提升从最小覆盖中获得的不等式,我们将在本节中描述。对于任意0/1多面体,Peled (1977) 和 Zemel(1978)给出了这个结果的推广。Wolsey (1976a)讨论了一般整数程序的提升过程。为了说明Balas和Zemel的结果,需要引入额外的符号。

回想一下,定理28描述了容易计算的有效不等式,如果条件(11)满足,甚至可以依次提升最小覆盖不等式。设我是满足条件(11)a的 j ∈ [ n ] \ C j \in[n] \backslash C j[n]\C 的指标集合 J : = [ n ] \ ( C ∪ I ) J:=[n] \backslash(C \cup I) J:=[n]\(CI) 。给定集合 J J J,定义

在这里插入图片描述

对于每个 M ∈ M ( J ) M \in \mathcal{M}(J) MM(J) 定义 β M ′ \beta_M^{\prime} βM 为:

在这里插入图片描述

定理38:(Balas and Zemel 1978,定理9)设 a ∈ R + n , β ∈ R + a \in \mathbb{R}_{+}^n, \beta \in \mathbb{R}_{+} aR+n,βR+ ,设假设11成立,设 C C C P P P 的最小覆盖,设 π i π_i πi 定义为定理28。这是简单LCI

在这里插入图片描述

facet定义 P P P 当且仅当

(a) 对于每一个 i ∈ J i∈J iJ ,存在 δ i ∈ [ 0 , 1 ] \delta_i \in[0,1] δi[0,1] ,使得

在这里插入图片描述
(b) 向量 δ ∈ R J \delta \in \mathbb{R}^J δRJ 是多面体的一个顶点

在这里插入图片描述

推论39:(Balas and Zemel 1978,推论9.1)极小覆盖C的提升系数 α i α_i αi 满足定理38的(a)和(b)的简单LCI(9)是一个顺序提升的极小覆盖不等式当且仅当 δ ∈ { 0 , 1 } J \delta \in\{0,1\}^J δ{0,1}J

这意味着定义具有整数系数的简单LCIs的所有facet都可以通过顺序提升生成,而同时提升的最小覆盖不等式通过顺序提升无法找到具有分数系数。这些可以通过计算多面体T的分数顶点得到。另一个含义是对顺序提升面的提升系数的描述。

推论40::(Balas and Zemel 1978,定理3)设(9)为顺序简单的提升覆盖不等式,定义p的一个面,则提升系数 α i , i ∈ [ n ] \ C \alpha_i, i \in[n] \backslash C αi,i[n]\C 满足

在这里插入图片描述

此外,Hartvigsen和Zemel(1992)证明了可以在 O ( n 2 ) O(n^2) O(n2) 时间内确定有效不等式 c ⊤ x ≤ γ c^{\top} x \leq \gamma cxγ c ∈ Z + n c \in \mathbb{Z}_{+}^n cZ+n 是一个(facet定义)简单提升的最小覆盖不等式。然而,如果 c ∈ R n c \in \mathbb{R}^n cRn,决定 c ⊤ x ≤ γ c^{\top} x \leq \gamma cxγ 定义了通过提升最小覆盖不等式获得的有效不等式,即使底层覆盖是已知的,它也是NP完备的。

示例41:(继续示例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 和最小覆盖 C = { 2 , 3 , 4 } C =\{2,3,4\} C={2,3,4}。在示例30中:

在这里插入图片描述

π 1 = 0 , π 5 = 1 π_1 = 0, π_5 = 1 π1=0π5=1。可得 I = ∅ I=\varnothing I= J = { 1 , 5 } J =\{1,5\} J={1,5} M ( J ) = { { 1 } , { 5 } , { 1 , 5 } } \mathcal{M}(J)=\{\{1\},\{5\},\{1,5\}\} M(J)={{1}{5}{1,5}} 。对于 M = { 1 } M =\{1\} M={1} ,系数 β M ′ = β { 1 } ′ \beta_M^{\prime}=\beta_{\{1\}}^{\prime} βM=β{1} 由下式给出

在这里插入图片描述

这意味着不平等

在这里插入图片描述

对于 M = { 5 } M =\{5\} M={5} M = { 1 , 5 } M =\{1,5\} M={1,5} ,可以类似地导出不等式 δ 5 ≤ 1 \delta_5 \leq 1 δ51 δ 1 + δ 5 ≤ 1 \delta_1+\delta_5 \leq 1 δ1+δ51 。根据定理38,不等式

在这里插入图片描述

facet定义 P P P 当且仅当

(a)

在这里插入图片描述
(b) δ = ( δ 1 , δ 5 ) ⊤ ∈ R 2 \delta=\left(\delta_1, \delta_5\right)^{\top} \in \mathbb{R}^2 δ=(δ1,δ5)R2 是多面体的一个顶点

在这里插入图片描述

因为T有两个顶点 ( 1 , 0 ) ⊤ (1,0)^{\top} (1,0) ( 0 , 1 ) ⊤ (0,1)^{\top} (0,1) 呢?,通过顺序提升或同时提升从最小覆盖 C = { 2 , 3 , 4 } C =\{2,3,4\} C={2,3,4} 得到的定义不等式的仅有两个面是

在这里插入图片描述

这里不等式(18)对应顶点 ( 1 , 0 ) ⊤ (1,0)^{\top} (1,0) 不等式(19)对应顶点 ( 0 , 1 ) ⊤ (0,1)^{\top} (0,1) T T T 只有整顶点,推论29表明两个面定义不等式都可以通过顺序提升得到。

对于最小覆盖 C = { 1 , 2 , 3 } C =\{1,2,3\} C={1,2,3} ,则 I = { 4 , 5 } I = \{4,5\} I={4,5} J = ∅ J=\varnothing J= 。因此,不等式

在这里插入图片描述

是唯一可以从C连续或同时提升得到的面。

在例30中,我们已经总结了这种独特性——至少对于连续升降是这样。

(Q11) 对于不同于封面的被解除不等式,我们能否(依次地或同时地)发展出类似的理论?例如,是否可以描述定义由解除包的不平等引起的不平等的方面(见第六章)?


相关文章:

  • 邢台做wap网站多少钱/网络优化推广公司哪家好
  • 检测网站空间容量/苏州seo关键词优化软件
  • 做企业网站需要注意什么/软文广告文案案例
  • 网站建设条款/seo推广方法有哪些
  • 用电脑做服务器的建一个网站/百度热搜榜排名
  • wordpress加个微信登录/磁力搜索引擎不死鸟
  • 数据结构 - AVL树 (Adelson-Velsky and Landis Tree)
  • 华为OD机试真题 C++ 实现【货币单位换算】【2022.11 Q4 新题】
  • 非零基础自学Golang 第16章 正则表达式 16.3 regexp包 16.4 小结 16.5 知识拓展
  • Linux系统 PHP安装expect扩展详解
  • (二十)Vue之非单文件组件
  • always块中时序逻辑 negedge rst_n和posedge rst实际电路
  • es笔记五之term-level的查询操作
  • 用简单伪随机数发生器实现随机中点位移分形(Matlab代码实现)
  • Macos安装和卸载第三方软件的正确方法
  • 02-Spring Boot启动原理核心源码剖析
  • 数据库实验5 数据库设计实验
  • Opencv(C++)笔记--直方图均衡化、直方图计算