在农夫约翰的农场上,每逢下雨,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
分析:单纯的网络流最大流算法。
给个不错的 参考资料(nocow网络流)
【参考程序】://bfs+dfs,速度不错
/*
ID: XIONGNA1
PROG: ditch
LANG: C++
*/
#include<iostream>
#include<cstring>
using namespace std;
int c[201][201],pre[201],queue[201],dis[201];
bool vis[201];
int n,m,ans;
inline int find_min(int a,int b)
{
return a<b?a:b;
}
void dfs(int s)
{
if (!c[pre[s]][s] && s!=1) return ;
if (s==n)
{
int i,flow=0xFFFFFFF;
i=n;
while (i!=1)
{
flow=find_min(flow,c[pre[i]][i]);
i=pre[i];
}
ans+=flow;
i=n;
while (i!=1)
{
c[pre[i]][i]-=flow;
c[i][pre[i]]+=flow;
i=pre[i];
}
}
for (int i=1;i<=n;i++)
if (c[pre[s]][s]>0 || s==1)
if (c[s][i]>0 && dis[i]==dis[s]+1)
{
pre[i]=s; dfs(i);
}
}
bool bfs()
{
int head,tail,now;
for (int i=1;i<=n;i++)
{
dis[i]=0; vis[i]=false;
}
queue[1]=1; vis[1]=true;
head=tail=1;
while (head<=tail)
{
now=queue[head];
for (int i=1;i<=n;i++)
if (!vis[i] && c[now][i]>0)
{
vis[i]=true;
tail++; queue[tail]=i;
dis[i]=dis[now]+1;
}
head++;
}
return vis[n];
}
int main()
{
freopen("ditch.in","r",stdin);
freopen("ditch.out","w",stdout);
scanf("%d%d",&m,&n);
memset(c,0,sizeof(c));
int a,b,s;
for (int i=1;i<=m;i++)
{
scanf("%d%d%d",&a,&b,&s);
c[a][b]+=s;
}
ans=0;
memset(pre,0,sizeof(pre));
while (bfs()) dfs(1);
printf("%d\n",ans);
return 0;
}
【参考程序】://单纯bfs找可行流
/*
ID: XIONGNA1
PROG: ditch
LANG: C++
*/
#include<iostream>
#include<cstring>
using namespace std;
int c[201][201],pre[201],queue[201];
bool vis[201];
int n,m,ans;
inline int find_min(int a,int b)
{
return a<b?a:b;
}
void change()
{
int i,flow=0xFFFFFFF;
i=n;
while (i!=1)
{
flow=find_min(flow,c[pre[i]][i]);
i=pre[i];
}
ans+=flow;
i=n;
while (i!=1)
{
c[pre[i]][i]-=flow;
c[i][pre[i]]+=flow;
i=pre[i];
}
}
bool bfs()
{
for (int i=1;i<=n;i++)
{
vis[i]=false; pre[i]=0;
}
int head=1,tail=1,now;
vis[1]=true; queue[1]=true;
while (head<=tail)
{
now=queue[head];
for (int i=1;i<=n;i++)
if (!vis[i] && c[now][i]>0)
{
vis[i]=true;
tail++;
queue[tail]=i; pre[i]=now;
}
head++;
}
return vis[n];
}
int main()
{
freopen("ditch.in","r",stdin);
freopen("ditch.out","w",stdout);
scanf("%d%d",&m,&n);
memset(c,0,sizeof(c));
int a,b,s;
for (int i=1;i<=m;i++)
{
scanf("%d%d%d",&a,&b,&s);
c[a][b]+=s;
}
ans=0;
while (bfs()) change();
//for (int i=1;i<=n;i++) ans+=c[n][i];
printf("%d\n",ans);
return 0;
}
【参考程序】://随便一点的硬搜
/*
ID: XIONGNA1
PROG: ditch
LANG: C++
*/
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
long n,m,i,j,ans;
long p[201],flow[201],cost[201][201];
long visit[201];
long minm(int a,int b)
{
if (a<b) return a;
else return b;
}
int main()
{
freopen("ditch.in","r",stdin);
freopen("ditch.out","w",stdout);
scanf("%d%d",&n,&m);
memset(cost,0,sizeof(cost));
int a,b,s;
for (i=1;i<=n;i++)
{
scanf("%d%d%d",&a,&b,&s);
cost[a][b]+=s;
}
ans=0;
do
{
memset(visit,false,sizeof(visit));
memset(p,0,sizeof(p));
memset(flow,0,sizeof(flow));
flow[1]=0xfffffff;
int maxm,maxl;
while (true)
{
maxm=maxl=0;
for (i=1;i<=m;i++)
if (!visit[i] && flow[i]>maxm)
{
maxm=flow[i];
maxl=i;
}
if (maxm==0 || maxl==m) break;
visit[maxl]=true;
for (i=1;i<=m;i++)
{
s=minm(maxm,cost[maxl][i]);
if (flow[i]<s)
{
flow[i]=s;
p[i]=maxl;
}
}
}
if (flow[m]==0) break;
ans+=flow[m];
s=flow[m];
i=m;
while (i!=1)
{
j=p[i];
cost[j][i]-=s;
cost[i][j]+=s;
i=j;
}
}while (true);
printf("%d\n",ans);
return 0;
}