【精选】最大流与最小费用流.pdfVIP

最大流与最小费用流 前两个星期都在看最大流与费用流，在找标程时几乎所有的都是 C++的，少数几个 Pascal 的，要么代码风格乱七八糟，要么是不含任何优化的，我照着其中一个写下来最后 发现居然是错的！最后请教了 cyx 并摸索了他的C++代码，终于“修成正果”。把自己的心 得贴在这里，给那些跟我一样看不懂 C 的在最大流与费用流中探索的 OIers 一点小小的 help，顺便给自己攒点人品（或许我该去转C 了，毕竟要跟得上时代）。 最大流我写的是Dinic 算法，因为找不到传说中特别高效的SAP，所以就将就一下。事 实上我写了几道最大流的题目，发现Dinic 效率还是不错的，而且编程复杂度又低，在OI 中还是蛮实用的。 算法思想就是建立残留网络，每次找到在残留网络中的一条增广路进行增广，添加反向 边的目的就是用来退流的（想想为什么？）。而Dinic 在这基础上进行改进，就是构建层次 图，沿层次图进行多路增广（不懂的话详细看代码）。 程序如下： program maxflow; var E,w,pre:array[0..20000] of longint; pot,dist,q:array[0..10000] of longint; n,m,s,t,ans,L:longint; procedure add(a,b,c:longint); //用链式前向星存储 begin E[L]:=b; //E 数组存点 w[L]:=c; //w 数组存每条边的容量 pre[L]:=pot[a]; //pre 数组记录起点同样为a 的前一条边 pot[a]:=L; //pot 数组为索引表 inc(L) //数组要从 0 开始存，否则后面的xor 会出错，我因这个调了2 小时 end; procedure init; //直接建立残留网络，不用存原图 var i,a,b,c:longint; begin L:=0; fillchar(pre,sizeof(pre),255); //这两个数组要初始化为-1，与后面调用 fillchar(pot,sizeof(pot),255); 对应 readln(n,m); for i:=1 to m do begin readln(a,b,c); add(a,b,c); //在残留网络中添加一条a 到b 剩余容量为c 的边 add(b,a,0) //添加反向边 end; s:=1; t:=n //s 为源，t 为汇 end; function bfs:boolean; //宽搜进行构建层次图（到源点的最短距离）， var k,head,tail:longint; 因为每条边在这的距离权值都一样（为 1 ）， begin 所以不用写成spfa，cyx 所教Orz fillchar(q,sizeof(q),0); fillchar(dist,sizeof(dist),255); //初始化为-1，用于退出时的判定 head:=1; tail:=2; q[head]:=s; dist[s]:=0; while headtail do //不需要采用循环队列,显然每个点只入队一次 begin k:=pot[q[head]]; while k-1 do begin if (w[k]0) and (dist[E[k]]=-1) then begin dist[E[k]]:=dist[q[head]]+1; q[tail]:=E[k]; inc(tail) end; k:=pre[k] end; inc(head) end; if dist[t]-1 then exit(true) //是否存在增广路 else exit(false) end; function min(a,b:longint):longint; begin if ab then min:=b else min:=a end; functi