运筹学中计算复杂性.pptVIP

;;;这个著名的例子是：;图3 三维立方体及其摄动; 1979年春，前苏联科学家Хачьян（L.G.khachian）证明了： LP有一个多项式时间算法！ 其结果建立在其他数学工作者关于NLP工作的基础上,方法与以前解决LP的途径截然不同,几乎完全不管LP的组合性质.;2、椭球算法的思想 ; 3、LP与线性不等式的关系;;;为了书写简洁,总可以将(1)写为下面的形式: aiX≤bi i=1,2,…,m;（4） ; 4、Хачьян(哈奇扬)算法的计算步骤;（1）n=2时，严格不等式组为 ai1x1+ai2x2〈 bi i=1，…，m ? 其解集合为若干个半平面的公共部分。 用 表示原点（0，0）与点X=（x1，x2）之间的距离，即 = ，而 ≤R代表以原点为圆心，R为半径的圆。;①的解集合在圆 ≤2L内的部分P的面积至少是 2-（n+1）L。;（L为输入长度，n为未知数个数，本例中n=2），如果迭代了6n2L次还没有找到①的解，就可以判定①一定没有解。;（2）从{X（i），Ei}得到{X（i+1），Ei+1}的过程： 首先将X（i）代入不等式组①中去试验，看是否满足，是，则X（i）就是①的解，计算停止。 若不是，将X（i）=（x1（i），x2（i））代入①中，至少有一个不等式不成立，不妨设为第r个，则有 ;第一步，过椭圆的中心X（i）画一条直线l与直线ar1x1+ar2x2=br平行，就把椭圆Ei分成了两个“半椭圆”，设与Ei交于A、B两点。 第二步，再画直线l1 与l平行,且与椭圆Ei相切。要求l1和l要在直线ar1x1+ar2x2=br的两侧，设l1 与Ei相切于C点。 第三步，作椭圆Ei+1，满足： ? 经过A、B、C三个点； ? 与l1 相切于C点； ? 包含Ei被直线l 割出的介于l 与l1 的半个椭圆。 ;;（4）综上，可以证明这种迭代方法最多只要进行6n2L次，其计算复杂性为O（6n2L2）。 ;（1）是理论上的重大突破；;? Хачьян算法虽然是多项式算法，为什麽解具体问题并不十分有效？而且远不如单纯形法？ ? 解决问题的途径： A . 走 Хачьян的路，继续进行改进——试图寻找与L无关的强多项式算法； B . 改造单纯形法，试图使之成为好算法。 比如：设计各种??列法则，以寻求到达最优点的最短路；;二、Karmarkar算法; 去掉变量的非负条件，走的路穿过多面体内部通向最优解，可惜的是：走得太慢，算法无效！但是人们发现，加上 函数，能改变方向进入内部，若能克服步长太小走得太慢的问题，把方向对准，步长加大，就有可能得到改善。;Karmarkar的贡献; 2、一些基本概念和重要结果; 或 ;;其中 D= ; 给定可行域（超凸多面体）P的一个内点a=(a1,…,an)T,即a∈P,且a0,首先通过一个投影变换把P和a分别变成S和a‘,使a’是S的中心,而且以a‘为中心包含S的最小球的半径与以a’为中心包含在P‘中的最大球的半径之比为O(n),然后在变换之后的以a’为心的球上,求得LP的解为b’,接着利用投影变换把b返回到原决策空间中去. ; 计算复杂性为O(n4L2)， 若采用某些技巧则可达到O(n3.5L2)。