首页 > 代码库 > 最小费用流spfa算法模板(pascal)
最小费用流spfa算法模板(pascal)
以前写过,现在的码风与以前有些变化,主要是用数组模拟邻接表存图,以前是用指针存图。
以前的博文:http://www.cnblogs.com/Currier/p/6387732.html
洛谷可评测。
传送门:https://www.luogu.org/problem/show?pid=3381
1 program rrr(input,output); 2 const 3 inf=123456789; 4 type 5 etype=record 6 t,c,w,next,rev:longint; 7 end; 8 var 9 e:array[0..100010]of etype; 10 a,fre,frv,q,dis:array[0..5050]of longint; 11 inq:array[0..5050]of boolean; 12 n,m,s,t,i,j,x,y,c,w,max,ans,cnt,f:longint; 13 procedure add(x,y,c,w:longint); 14 begin 15 inc(cnt);e[cnt].t:=y;e[cnt].c:=c;e[cnt].w:=w;e[cnt].next:=a[x];a[x]:=cnt; 16 end; 17 function min(a,b:longint):longint; 18 begin 19 if a<b then exit(a) else exit(b); 20 end; 21 procedure spfa; 22 var 23 h,t:longint; 24 begin 25 for i:=1 to n do dis[i]:=inf; 26 fillchar(inq,sizeof(inq),false); 27 h:=0;t:=1;q[1]:=s;dis[s]:=0;inq[s]:=true; 28 while h<>t do 29 begin 30 inc(h);if h>5000 then h:=1; 31 i:=a[q[h]]; 32 while i<>0 do 33 begin 34 if (e[i].c>0) and (dis[q[h]]+e[i].w<dis[e[i].t]) then 35 begin 36 dis[e[i].t]:=dis[q[h]]+e[i].w; 37 frv[e[i].t]:=q[h];fre[e[i].t]:=i; 38 if not inq[e[i].t] then 39 begin 40 inc(t);if t>5000 then t:=1; 41 q[t]:=e[i].t;inq[e[i].t]:=true; 42 end; 43 end; 44 i:=e[i].next; 45 end; 46 inq[q[h]]:=false; 47 end; 48 end; 49 begin 50 assign(input,‘r.in‘);assign(output,‘r.out‘);reset(input);rewrite(output); 51 readln(n,m,s,t); 52 fillchar(a,sizeof(a),0);cnt:=0; 53 for i:=1 to m do 54 begin 55 readln(x,y,c,w); 56 add(x,y,c,w);e[cnt].rev:=cnt+1; 57 add(y,x,0,-w);e[cnt].rev:=cnt-1; 58 end; 59 max:=0;ans:=0; 60 while true do 61 begin 62 spfa; 63 if dis[t]=inf then break; 64 j:=t;f:=inf; 65 while j<>s do begin f:=min(f,e[fre[j]].c);j:=frv[j]; end; 66 max:=max+f; 67 j:=t;w:=0; 68 while j<>s do begin w:=w+e[fre[j]].w;dec(e[fre[j]].c,f);inc(e[e[fre[j]].rev].c,f);j:=frv[j]; end; 69 ans:=ans+w*f; 70 end; 71 write(max,‘ ‘,ans); 72 close(input);close(output); 73 end.
最小费用流spfa算法模板(pascal)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。