描述
在农夫约翰的农场上,每逢下雨,Bessie最喜欢的三叶草地就积聚了一潭水。这意味着草地被水淹没了,并且小草要继续生长还要花相当长一段时间。因此,农夫约翰修建了一套排水系统来使贝茜的草地免除被大水淹没的烦恼(不用担心,雨水会流向附近的一条小溪)。作为一名一流的技师,农夫约翰已经在每条排水沟的一端安上了控制器,这样他可以控制流入排水沟的水流量。
农夫约翰知道每一条排水沟每分钟可以流过的水量,和排水系统的准确布局(起点为水潭而终点为小溪的一张网)。需要注意的是,有些时候从一处到另一处不只有一条排水沟。
根据这些信息,计算从水潭排水到小溪的最大流量。对于给出的每条排水沟,雨水只能沿着一个方向流动,注意可能会出现雨水环形流动的情形。
[编辑]格式
PROGRAM NAME:ditch
INPUT FORMAT:
(file ditch.in)
第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)是这条排水沟的最大容量。
OUTPUT FORMAT:
(file ditch.out)
输出一个整数,即排水的最大流量。
[编辑]SAMPLE INPUT
5 4 1 2 40 1 4 20 2 4 20 2 3 30 3 4 10
[编辑]SAMPLE OUTPUT
50
USACO的一道网络流入门题目。
初学网络流,不是太明白的时候做的题。
裸题一道。有一点一开始没有搞清楚,就是标号过程就在sap中进行了。一开始我还单独写了一个bfs
#include <cstdio>
#include <queue>
using std::deque;
using std::make_pair;
deque<std::pair<long,long> > que;
const long inf = 0x7fff0000;
long g[1010][1010];
long c[1010][1010];
long f[1010][1010];
long n;long m;
long d[1010];
long vd[1010];
inline long min(long a,long b)
{
if (a < b)
return a;
return b;
}
long dfs(long u,long flow)
{
if (u == n)
{
return flow;
}
long count = 0;
for (long v=1;v<n+1;v++)
{
if (g[u][v]>0&&d[u]-d[v]==1)
{
long tmp = dfs(v,min(flow-count,g[u][v]));
g[u][v] -= tmp;
g[v][u] += tmp;
if (c[u][v]>0)
f[u][v] += tmp;
else
f[v][u] -= tmp;
count += tmp;
if (count == flow)
{
return flow;
}
}
}
// if (d[1] == 0)
if (d[1] >= n)
// return flow;
return count;
vd[d[u]]--;
if (vd[d[u]]==0)
{
d[1] = n;
}
d[u]++;
vd[d[u]]++;
// return flow;
return count;
}
//void bfs()
//{
// d[1] = n;
// que.push_back(make_pair(1,n));
// while (!que.empty())
// {
// long u = que.front()->first;
// long lev = que.front()->second;
// que.pop_front();
// for (long v=1;v<n+1;v++)
// {
// if (d[v]==0&&c[u][v]>0)
// {
// d[v] = lev-1;
// que.push_back(make_pair(v,lev-1));
// }
// }
// }
//}
int main()
{
scanf("%ld%ld",&m,&n);
for (long i=1;i<m+1;i++)
{
long s;long e;long cc;
scanf("%ld%ld%ld",&s,&e,&cc);
// c[s][e] = cc;
// g[s][e] = cc;
c[s][e] += cc;
g[s][e] += cc;
}
// bfs();
long ans = 0;
vd[0] = n;
while (d[1]<n)
{
ans += dfs(1,inf);
}
printf ("%ld",ans);
return 0;
}