正文
洛谷P2740 [USACO4.2]草地排水Drainage Ditches
提示：扫一扫查出行【扫一扫了解最新限行尾号】
复制提示
题目背景
在农夫约翰的农场上，每逢下雨，贝茜最喜欢的三叶草地就积聚了一潭水。这意味着草地被水淹没了，并且小草要继续生长还要花相当长一段时间。因此，农夫约翰修建了一套排水系统来使贝茜的草地免除被大水淹没的烦恼（不用担心，雨水会流向附近的一条小溪）。作为一名一流的技师，农夫约翰已经在每条排水沟的一端安上了控制器，这样他可以控制流入排水沟的水流量。
题目描述
农夫约翰知道每一条排水沟每分钟可以流过的水量，和排水系统的准确布局（起点为水潭而终点为小溪的一张网）。需要注意的是，有些时候从一处到另一处不只有一条排水沟。
根据这些信息，计算从水潭排水到小溪的最大流量。对于给出的每条排水沟，雨水只能沿着一个方向流动，注意可能会出现雨水环形流动的情形。
输入输出格式
输入格式：
第1行: 两个用空格分开的整数N (0 <= N <= 200) 和 M (2 <= M <= 200)。N是农夫John已经挖好的排水沟的数量，M是排水沟交叉点的数量。交点1是水潭，交点M是小溪。
第二行到第N+1行: 每行有三个整数，Si, Ei, 和 Ci。Si 和 Ei (1 <= Si, Ei <= M) 指明排水沟两端的交点，雨水从Si 流向Ei。Ci (0 <= Ci <= 10,000,000)是这条排水沟的最大容量。
输出格式：
输出一个整数，即排水的最大流量。
输入输出样例
输入样例#1：
5 4
1 2 40
1 4 20
2 4 20
2 3 30
3 4 10
输出样例#1：
50
说明
输入样例#1：
5 4
1 2 40
1 4 20
2 4 20
2 3 30
3 4 10
输出样例#1：
50
题目翻译来自NOCOW。
USACO Training Section 4.2
裸的最大流
dinic：
program rrr(input,output);
const
inf=;
type
pointer=^nodetype;
nodetype=record
t,c:longint;
next,rev:pointer;
end;
var
a,cur:array[..]of pointer;
q,d:array[..]of longint;
p:pointer;
n,m,i,x,y,c,ans,h,t:longint;
procedure ins(x,y,c:longint);
begin
new(p);p^.t:=y;p^.c:=c;p^.next:=a[x];a[x]:=p;
end;
procedure add(x,y,c:longint);
begin
ins(x,y,c);ins(y,x,);
a[x]^.rev:=a[y];a[y]^.rev:=a[x];
end;
function min(a,b:longint):longint;
begin
if a<b then exit(a) else exit(b);
end;
procedure bfs;
begin
for i:= to n do d[i]:=-;
h:=;t:=;q[]:=;d[]:=;
while h<t do
begin
inc(h);
p:=a[q[h]];
while p<>nil do
begin
if (d[p^.t]=-) and (p^.c>) then
begin
d[p^.t]:=d[q[h]]+;
inc(t);q[t]:=p^.t;
end;
p:=p^.next;
end;
end;
end;
function dfs(k,f:longint):longint;
var
ans,r:longint;
p:pointer;
begin
if (k=n) or (f=) then exit(f);
ans:=;
p:=cur[k];
while p<>nil do
begin
if (d[p^.t]=d[k]+) and (p^.c>) then
begin
r:=dfs(p^.t,min(f,p^.c));
p^.c:=p^.c-r;
p^.rev^.c:=p^.rev^.c+r;
ans:=ans+r;
f:=f-r;
if f= then break;
end;
p:=p^.next;
cur[k]:=p;
end;
if f> then d[k]:=-;
exit(ans);
end;
begin
//assign(input,'r.in');assign(output,'r.out');reset(input);rewrite(output);
readln(m,n);
for i:= to n do a[i]:=nil;
for i:= to m do begin read(x,y,c);add(x,y,c); end;
ans:=;
while true do
begin
bfs;
if d[n]=- then break;
for i:= to n do cur[i]:=a[i];
ans:=ans+dfs(,inf);
end;
write(ans);
//close(input);close(output);
end.
