最小费用最大流 网络流.docVIP

最小费用最大流程序 Program?Max_Flow_Min_Cost;??//注意:图中任意两点之间至多有一条直连通有向边!! ??type?nettype=record ????????????????c,?w?,f:longint;?//弧的容量，每单位通过量费用，可行流量 ???????end; ???????btype=record ??????????????v,fa:longint;?????//到源点的最短距离，父辈 ???????end;? ??var??net:array[1..100,1..100]of?nettype;?//网络N的邻接矩阵 ??????????best:array[1..100]of?btype;??????//求最短路时用的数组 ??????????n,s,t:word;????????????????????//顶点数，源点、汇点的编号 ??????????minc:longint;?????????????????//最小总费用 ??procedure?init;???????????????//边目录输入? ????var?a,b:word;c,d:longint; ????begin ??????minc:=0; ??????fillchar(net,sizeof(net),0); ??????assign(input,flow.in);reset(input); ??????readln(n,s,t);????????????//注意:?这里st!!?**允许s到t无路(输出0)** ??????readln(a,b,c,d); ??????while?a+b+c+d0?do ????????begin ??????????net[a,b].c:=c;{给弧(a,b)赋容量c}//这里两点只有一条有向边!!-**问题关键**-? ??????????net[a,b].w:=d;{给弧(a,b)赋费用d} ??????????net[b,a].w:=-d;{给相反孤(b,a)赋费用-d} ??????????readln(a,b,c,d); ????????end; ??????close(input); ????end; ??procedure?print;?{输出最小费用最大流的费用及方案} ????var?i,j:word; ????begin ??????writeln(MinCost=,minc); ??????for?i:=1?to?n?do ???????for?j:=1?to?n?do ????????if?net[i,j].f0?then ??????????writeln(i,--,j,?:?,net[i,j].w,*,net[i,j].f); ????end; ? ? ? function?Find_Way:boolean;?{求最小费用增广路函数,若best[t].valueMaxInt,则找到增广路，函数返回true，否则返回false} ????var?i,j:word; ???????????quit:boolean; ????begin ??????for?i:=1?to?n?do?best[i].v:=maxlongint;?//寻求最小费用增广路径前，给best数组赋初值 ??????best[s].v:=0; ??????repeat ????????quit:=true; ????????for?i:=1?to?n?do ??????????if?best[i].vmaxlongint?then ????????????for?j:=1?to?n?do ??????????????if?(net[i,j].fnet[i,j].c)and(best[i].v+net[i,j].wbest [j].v)?then ????????????????begin ??????????????????best[j].v:=best[i].v+net[i,j].w; ??????????????????best[j].fa:=i; ??????????????????quit:=false; ????????????????end; ??????until?quit; ??????if?best[t].vmaxlongint?then?Find_Way:=true ??????????????????????????????else?Find_Way:=false; ????end; ? ??procedure?Add_Way; ????var?i,j:word; ????begin ??????i:=t; ??????repeat ????????