运筹学中的计算复杂性(2).pptVIP

第一阶段——从形式语言→形式语言 类的理论 第二阶段——1971年 COOK在逻辑、递归理论的研究中发现了NP完全类，证明了适定性问题是NP完备的。 1、OR是计算复杂性的生成基础 从计算复杂性发展历史来看： 专题二 《计算复杂性和运筹学》 2、OR对计算复杂性产生极大影响 最突出的就是80年代以来产生的哈奇扬（Хачьян）算法和卡马卡（Karmarkar）算法。 第三阶段——从事OR研究的Karp在网络、组合优化方面颇有建树，在研究了判断逻辑表达式伪真的证明后，发现OR中（比如组合优化中）就有很多NP完全类，从而提供了计算复杂性理论的生成基础。 1、问题的提出 1947年Dantzig创立了单纯形法 1972年，V．Klee G.Minty构造了一个反例,它含有n个变量,m=2n个不等式约束,若用单纯形法求解,必须检验约束条件中不等式组所确定的凸多面体的所有顶点,才能获得最优解,计算次数等于2n -1,因此说明了单纯形法不是“好算法”! 二、椭球算法 这个著名的例子是： 图3 三维立方体及其摄动 X3 x1 x2 1979年春，前苏联科学家Хачьян（L.G.khachian）证明了： LP有一个多项式时间算法！ 其结果建立在其他数学工作者关于NLP工作的基础上,方法与以前解决LP的途径截然不同,几乎完全不管LP的组合性质. 那麽，能否找到一种解LP的多项式算法？ 2、椭球算法的思想 对应的n阶方阵,而X(1)=0,E1是一个以X(1)为球心的n维超球. Хачьян首先提出了将求LP最优解归结为解严格不等式组的问题,解不等式组的多项式算法也是一种迭代算法,迭代的每一步都要产生一个点X(K)和以该点为中心的一个椭球EK=(X(K),BK),其中BK是与第K个椭球方程 (X-X(K))TBK(X-X(K))≤1 3、LP与线性不等式的关系 设要求的LP问题是: 对偶问题 若X、Y分别是它们的可行解，则由弱对偶定理，必有CX≥Yb；由最优性准则定理：若X*，Y*分别是上述问题的可行解，且CX*=Y*b，则X*，Y*分别为（L）和（D）的最优解。 要求解（L），可以先解如下的线性不等式组： AX≥b X≥0 YA≤C （1） Y≥0 CX≤Yb 若已求得（1）的一个最优解（X*，Y*），则X*就是（L）问题的一个最优解。若（1）无解，则问题（L）没有最优解。所以，求LP问题的最优解可以转化为求解线性不等式组（1）。而不等式组（1）总可以改写成如下形式： （1） -AX≤-b -X≤0 YA≤C -Y≤0 CX-Yb≤0 为了书写简洁,总可以将(1)写为下面的形式: aiX≤bi i=1,2,…,m (2) Хачьян证明了如果整系数不等式组 aiXbi+2*2-L i=1,…,m (3)有解,则整系数不等式组(2)也有解。其中， 其中,ai是系数矩阵的第i个行向量ai=(ai1,ai2,…, ain),bi是右端向量的第i个分量. （4） L大致等于把不等式组（2）的所有系数都化为二进制数时的位数，称L为问题（2）的输入长度，可大体上说明问题（2）的规模大小。 综上，求解LP（L）可以转化成求解不等式组： （*） 4、Хачьян(哈奇扬)算法的计算步骤 5、Хачьян(哈奇扬)算法的几何解释 作一个椭球，将可行域所有顶点包住，看球心是否在可行域内部？如果不是，用分离定理作超平面，再作一个缩小了的椭球，把问题化为一个不等式组是否有解的问题。 （1）n=2时，严格不等式组为 ai1x1+ai2x2〈 bi i=1，…，m ? 其解集合为若干个半平面的公共部分。 用 表示原点（0，0