线性规划讲义.ppt

如果b向量所有元素非负，则令所有x为0一定可以得到一组可行解。 那么我们可以通过最优化过程来逐步找到最优解，或指出它的解是无穷大或不存在。 这事实上对应着不断从单纯形的一个顶点走到一个更优的顶点最终达到最优解。 1. 找到一个x[e]，使c[e]0。 2. 选择一个基变量x[d]，使A[d][e]0，且最小化b[d]/A[d][e]。 3. 执行转轴操作pivot(d,e)，转到第一步继续算法。 由于最小化了b[d]/A[d][e]，所以不会影响b向量的非负性，因此所有变量取0依然是可行解。 由于 ，v一定不降。 因此达到了一个更优解。 算法终止条件： 1. 找不到一个c[e]0 2. 对于一个c[e]0，不存在A[d][e]0 1. 直观上理解，不存在一个pivot操作使v值增加 2. 此时可以让x[e]取+∞，目标函数值也为+∞。 若算法终止且有解，那么所有的c[e]≤0，所以令所有非基变量为0答案最大。 此时基变量x[n+d]=b[d]。 x1 x2 b x3 1 1 300 x4 2 1 400 x5 0 1 250 c 50 100 0 请手解右边的线性规划。 要求使用单纯形法。 x1 x2 b x3 1 1 300 x4 2 1 400 x5 0 1 250 c 50 100 0 x4 x2 b x3 -0.5 0.5 100 x1 0.5 0.5 200 x5 0 1 250 c -25 75 10000 pivot(2,1) x4 x3 b x2 -1 2 200 x1 1 -1 100 x5 1 -2 50 c 50 -150 25000 pivot(1,2) x4 x2 b x3 -0.5 0.5 100 x1 0.5 0.5 200 x5 0 1 250 c -25 75 10000 x5 x3 b x2 1 0 250 x1 -1 1 50 x4 1 -2 50 c -50 -50 27500 pivot(3,1) x4 x3 b x2 -1 2 200 x1 1 -1 100 x5 1 -2 50 c 50 -150 25000 x5 x3 b x2 1 0 250 x1 -1 1 50 x4 1 -2 50 c -50 -50 27500 此时所有的c[e]≤0，算法结束。 函数最大值为27500。 此时非基变量x5, x3取0， 基变量x2=250, x1=50, x4=50。 所以TYH该听250节YJQ的课，50节ZMS的课，他还剩50点体力LP。 注意：以上所有步骤都要求b向量非负，如果b向量为负该怎么办？ 随机一个b[d]0，再随机一个A[d][e]0，执行pivot(d,e)。 直到b向量全部非负为止。 效果：0ms AC 事实上n, m较大的情况下随机非·常·慢！ 这道题n和m只有20，因此随机化出乎意料地跑得很快。 正确的做法： 人工添加非基变量x0，使式子变成： 选择一个d使得b[d]最小，执行操作 pivot(d,0) 很明显，在这个操作后，b向量一定非负。 继续求解该线性规划。如果答案为负数，说明原单纯形不可能让所有基变量非负， 因此无可行解。否则，可以令x0为0，把x0从线性规划中消去，修改最大化的函数为原函数，继续求解。 每天要求至少的人数即为b向量，要求最小化的函数即为c向量。 由于原题要求最小化，将所有变量取负，先初始化可行解，使b向量为正后再跑单纯形。然而…… 这速度看起来差点就T了…… 再回到开始的问题…… TYH还没有去听课，忽然跳出一只YYR！YYR说：“TYH，你把你的COST、LP、法力水晶全给我吧，我来提高你的姿势水平！” 然后YYR为姿势水平定了价，规定1COST可以换y3的姿势水平，1LP可以换y4，1法力水晶可以换y5。 TYH很机智，如果花同样的代价，听课可以提高更多的姿势水平，他就会去听课。 YYR也很机智，他发现他需要花的时间和TYH提高的姿势水平成正比，所以他希望在拿到TYH所有财产的情况下，花的时间最少。 请根据刚刚的题意，再列一个关于YYR对姿势水平定价（y3, y4, y5）线性规划。 课程 消耗 ZMSTSL 变换 YJQTSL 变换 TYH拥有的总量 COST 1 1 300 体力LP 2 1 400 法力水晶 0 1 250 姿势水平提高 50 100 x1 x2 b x3 1 1 ≤ 300 x4 2 1 400 x5 0 1 250 Max 50 100 0 对偶问题 y3 y4 y5 b y1 1 2 0 ≥ 50 y2 1 1 1 100 Min 300 400 250 0 研究线性约束条件下线性目标函数的极值问题的数学理论和方法。 一般来说，形式如下：