【题解】CODEVS 1993 草地排水

CODEVS 1993 草地排水

Description

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

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

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

Input Description

第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 Description

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

Sample Input

5 4

1 2 40

1 4 20

2 4 20

2 3 30

3 4 10

Sample Output

50

模板题，直接上ISAP。

#include<cstdio>
#include<climits>
#include<cstring>
#include<algorithm>

const int MAXN=20010;
const int MAXM=500010;
using namespace std;

struct edge
{
	int v,next,val;
}e[MAXM];
int p[MAXN],eid;
int n,m;
int h[MAXN];
int gap[MAXN];
int S,T;

inline void add(int from,int to,int val)
{
    e[eid].v=to;e[eid].val=val;
    e[eid].next=p[from];
    p[from]=eid++;
}

inline int dfs(int pos,int cost)
{
    if(pos==T)
		return cost;
    int j,minh=n-1,lv=cost,d;
    for(j=p[pos];j!=-1;j=e[j].next)
    {
        int v=e[j].v,val=e[j].val;
        if(val>0)
        {
            if(h[v]+1==h[pos])
            {
                if(lv<e[j].val)
					d=lv;
                else
					d=e[j].val;
                d=dfs(v,d);
                e[j].val-=d;
                e[j^1].val+=d;
                lv-=d;
                if(h[S]>=n)
					return cost-lv;
                if(lv==0)
					break;
            }
            if (h[v]<minh)
				minh=h[v];
        }
    }
    if (lv==cost)
    {
        --gap[h[pos]];
        if(gap[h[pos]]==0)
			h[S]=n;
        h[pos]=minh+1;
        ++gap[h[pos]];
    }
    return cost-lv;
}

int isap(int st,int ed)
{
    S=st;
	T=ed;
    int ret=0;
    memset(gap,0,sizeof(gap));
    memset(h,0,sizeof(h));
    gap[st]=n;
    while(h[st]<n)
        ret+=dfs(st,INT_MAX);
    return ret;
}
int main()
{
    scanf("%d%d",&m,&n);
    memset(p,-1,sizeof(p));
    eid=0;
    for(int i=1;i<=m;++i)
    {
        int u,v,c;
        scanf("%d%d%d",&u,&v,&c);
        add(u,v,c);
        add(v,u,0);
    }
    printf("%d\n",isap(1,n));
    return 0;
}