第七节最小费用流问题.ppt

第七节 最小费用流问题 最小费用流算法 规划形式 算法步骤 算法复杂性 运输问题 对偶规划 算法步骤 规划形式 规划形式 续一 原始规划P 规划形式 续二 对偶规划D 最小费用流算法的步骤 第1步(开始) 让所有的流xij=0，所有点对应的数p(i)=0。 第2步(决定哪些弧可以改变流量) 用I表示这样的弧集合使p(j)-p(i)=wij，同时xijcij 用R表示这样的弧集合使p(j)-p(i)=wij，同时xij0 用Q表示不在I∪R中的弧集合。 最小费用流算法的步骤 第3步(改变流量)用最大流算法，在I∪R上找增广路，增加流量。 如果从s到t的流量已经是v，那么计算停止，已经得到一个流量是v的最小费用流。 如果从s到t不能再增加流量，检查在Q∪I∪R中是否能找到增广路。如果不能找到增广路，那么该网络的最大流达不到v；如果在Q∪I∪R中能找到增广路，就转入第4步。 第4步(改变顶点的p(i )值) 把所有点分成两类：S是标上号的点的集合，T是标不上号的点的集合，让S中的点p(i )值不变，T的点p(i )值全部加1，再转回第2步。 续 例 最小费用流算法的复杂性 设网络的点数为n，弧数为m，弧的最大容量为w 算法的循环次数取决于点的p(i)值修正的次数，最多进行mw次 因此第2步的计算量为O(m2w) 如果最大流值为v，则留增广最多进行v次 所以第3步的计算量为O(nw) 第4步的计算量为O(nmw) 所以，总的计算量为O(m2w+nmw) 运输问题 ai表示发点i可供应的产品数量 bj表示收点j所需的产品数量 wij表示从发点i到收点j的单位产品运输费用 xij表示从发点i分配给收点j的产品数量 对偶规划 算法步骤 第1步(开始) 任给原始规划的可行解{xij} 第2步(计算对偶解) 对于原始规划的可行解{xij}，计算出对偶规划的一个解{ui,vj} 第3步(计算检验数) 对于对偶规划的解{ui,vj}，计算 若所有的 均非负，则计算结束，这时得到的{xij}和{ui,vj}分别为原始规划和对偶规划的最优解；否则，转第4步 算法步骤 续 第4步(调整原始可行解) 令 即选择xst进入基。对应于网络中，即在支撑树上加入弧(s,t)，从而得到一个回路。并选择其流量xst=θ，使这个回路上的流量通过加θ 或减θ 以达到去掉一条弧的目的，从而得到一个新的被改进的原始可行解{xst}，转第2步 例 基本假设： 给定有向网络G=(N,A,C,W)，其中cij表示弧(i,j)∈A的容量，wij表示单位流通过弧(i,j)∈A的费用。则流x=(xij)的费用为 问题： 求一个从s到t，其流值为v的流，使得流的费用达到最小。 * 求解下图所示网络中自点1到点4其值为3的最小费用流。其中每条弧上的第一个数表示单位流的费用，第二个数表示容量。 * 求如图所示运输问题的最优解 2 1 3 1 2 3 4 15 20 30 20 10 30 迭代过程 初始原可行解 x(0) 迭代过程 2 1 3 1 2 3 4 8 6 12 1 0 1 4 续一 对偶解 v1 v2 v3 v4 u1 8 _ 6 - 10 -2 9 8 0 u2 9 - 12 5 13 - 7 5 1 u3 14 2 9 -1 16 - 5 - 4 8 6 12 1 检验数 迭代过程 续二 迭代过程 15-