usaco-4.2.1Drainage Ditches草地排水 最大流模板

usaco-4.2.1Drainage Ditches草地排水

时间限制: 1 Sec  内存限制: 128 MB

题目描述

在农夫约翰的农场上，每逢下雨，贝茜最喜欢的三叶草地就积聚了一潭水。这意味着草地被水淹没了，并且小草要继续生长还要花相当长一段时间。因此，农夫约翰修建了一套排水系统来使贝茜的草地免除被大水淹没的烦恼（不用担心，雨水会流向附近的一条小溪）。作为一名一流的技师，农夫约翰已经在每条排水沟的一端安上了控制器，这样他可以控制流入排水沟的水流量。 

农夫约翰知道每一条排水沟每分钟可以流过的水量，和排水系统的准确布局（起点为水潭而终点为小溪的一张网）。需要注意的是，有些时候从一处到另一处不只有一条排水沟。

根据这些信息，计算从水潭排水到小溪的最大流量。对于给出的每条排水沟，雨水只能沿着一个方向流动，注意可能会出现雨水环形流动的情形。

输入

第1行:	两个用空格分开的整数N (0 <= N <= 200)和M (2 <= M <= 200)。N是农夫约翰已经挖好的排水沟的数量，M是排水沟交叉点的数量。交点1是水潭，交点M是小溪。
第二行到第N+1行:	每行有三个整数，Si, Ei,和Ci。Si和Ei (1 <= Si, Ei <= M)指明排水沟两端的交点，雨水从Si流向Ei。Ci (0 <= Ci <= 10,000,000)是这条排水沟的最大容量。

输出

输出一个整数，即排水的最大流量。

样例输入

5 4
1 2 40
1 4 20
2 4 20
2 3 30
3 4 10

样例输出

50

最大流裸题

DINIC

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<map>
#include<set>
#include<vector>
#include<queue>
using namespace std;
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'|ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return f*x;
}
const int N=1<<9;const int inf=0X3f3f3f3f;
int n,m,ecnt=1,ans,t,w,q[N],last[N],d[N];
struct EDGE{int to,nt,val;}e[N];
inline void add(int u,int v,int val){e[++ecnt]=(EDGE){v,last[u],val};last[u]=ecnt;e[++ecnt]=(EDGE){u,last[v],0};last[v]=ecnt;}
bool bfs()
{
	memset(d,0,sizeof(d));
	t=0,w=1;q[t]=1;d[1]=1;
	while(t<w)
	{
		int temp=q[t++];
		for(int i=last[temp];i;i=e[i].nt)
		if(e[i].val&&!d[e[i].to])
		{d[e[i].to]=d[temp]+1;q[w++]=e[i].to;}
	}
	if(d[n])return true;else return false;
}
inline int dfs(int u,int lim)
{
	if(u==n||lim==0)return lim;
	int ret=0;
	for(int i=last[u];i;i=e[i].nt)if(d[e[i].to]==d[u]+1)
	{
		int temp=dfs(e[i].to,min(e[i].val,lim));
		e[i].val-=temp;e[i^1].val+=temp;
		lim-=temp;ret+=temp;
		if(lim==0)break;
	}
	if(!ret)d[u]=-1;
	return ret;
}
int main()
{
	m=read();n=read();int x,y,z;
	for(int i=1;i<=m;i++)
	{x=read();y=read();z=read();add(x,y,z);}
	while(bfs())ans+=dfs(1,inf);
	printf("%d\n",ans);
}

ISAP

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<map>
#include<set>
#include<vector>
#include<queue>
using namespace std;
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'|ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return f*x;
}
const int N=1<<9;const int inf=0X3f3f3f3f;
int n,m,ecnt=1,ans=0,t,w,q[N],head[N],d[N],last[N],num[N],EXIT,c[N];
struct EDGE{int to,nt,val;}e[N];
inline void add(int u,int v,int val){e[++ecnt]=(EDGE){v,head[u],val};head[u]=ecnt;e[++ecnt]=(EDGE){u,head[v],0};head[v]=ecnt;}
inline int isap(int u,int lim)
{
	int flow=0,t;
	if(u==n)return lim;
	for(int i=last[u];i;i=e[i].nt)
	if(e[i].val&&d[e[i].to]+1==d[u])
	{
		flow+= t=isap(e[i].to,min(e[i].val,lim-flow));
		e[i].val-=t;e[i^1].val+=t;
		if(EXIT||lim==flow)return flow;
	}
	EXIT=--num[d[u]]==0;
	++num[++d[u]];
	last[u]=head[u];
	return flow;
}
int main()
{
	m=read();n=read();int x,y,z;
	for(int i=1;i<=m;i++)
	{x=read();y=read();z=read();add(x,y,z);}
	num[0]=n;
	while(EXIT==0)ans+=isap(1,inf);
	cout<<ans;
}