【费用流的一点改进】距离标号连续增广路费用流算法

其实标题中那个又长又霸气的名字只是用来唬人的，其实这个算法非常简单

zkw算法的重点其实就是引入了一个距离标号，让程序一边更新标号求最短路，一边增广，从而使速度大幅提高

可是zkw算法中每次调整标号只增广一次，感觉有点浪费，为什么不多次增广呢？

如果说这个标号让增广的步伐加快了的话，那多路增广就能加大增广的步子，一步多跨一点，也就能让速度更快一点

我所谓的距离标号连续增广路费用流算法本质就是 zkw+多路增广

具体实现的话，可以像dinic一样，每次增广时找到一个最靠前的满流边，退到那个地方，继续做

代码量并没有大多少，速度却可以飙升

实测发现，可能卡到zkw的稠密图一般卡不到这种算法，zkw飞速的图，这个算法光速（夸张了一点。。呵呵）

测试工作在盾盾的机子上做的，图片等下再问他要，到时发上来

附关键代码：

procedure imp; var i,j:longint; begin i:=t;j:=oo; while i<>s do begin if w[fir[i]]<=j then begin j:=w[fir[i]];kk:=pre[i]; end; i:=pre[i]; end; i:=t; while i<>s do begin dec(w[fir[i]],j); inc(w[-fir[i]],j); i:=pre[i]; end; inc(ans,j*(d[t]-d[s])); end; procedure dfs(i:longint); var j,k:longint; begin if i=t then begin ok:=true;imp;exit; end; v[i]:=true;j:=h[i]; while j<>0 do begin k:=point[j]; if (w[j]>0)and(d[i]-d[k]+c[j]<sl[k]) then sl[k]:=d[i]-d[k]+c[j]; if (w[j]>0)and(d[k]=d[i]+c[j])and not v[k] then begin pre[k]:=i;fir[k]:=j; dfs(k); if ok then if i=kk then ok:=false else exit; end; j:=next[j]; end; end; function find:boolean; var i,j:longint; begin j:=oo; for i:=1 to t do if not v[i]and(sl[i]<j) then j:=sl[i]; if j=oo then exit(false); for i:=1 to t do if v[i] then dec(d[i],j); find:=true; end; procedure link(x,y,z,cc:longint); begin inc(e); next[e]:=h[x];point[e]:=y;w[e]:=z;c[e]:=cc;h[x]:=e; next[-e]:=h[y];point[-e]:=x;w[-e]:=0;c[-e]:=-cc;h[y]:=-e; end; 主程序调用： while find do begin fillchar(v,sizeof(v),0); fillchar(sl,sizeof(sl),127); dfs(s); end; writeln(ans); 如果要多次调用的话，需要先给v和sl清零

p.s KM用到的slack优化zkw也是可以的，具体见程序（slack数组就是sl[i]）