Codevs 1993 草地排水 [dinic] [最大流]

1993 草地排水 USACO
时间限制: 2 s
空间限制: 256000 KB
题目等级 : 钻石 Diamond

题目描述 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

---

dinic裸题。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
const int INF=0x3f3f3f3f;
const int maxn = 205;
struct Edge
{
    int to,next;
    int cap;
}edge[maxn<<1];
int head[maxn];
int maxedge;
inline void addedge(int u,int v,int c)
{
    edge[++maxedge] = (Edge) { v,head[u],c };
    head[u] = maxedge;
    edge[++maxedge] = (Edge) { u,head[v],0 };
    head[v] = maxedge;
}
int n,m;
int S,T;
void init()
{
    scanf("%d%d",&m,&n);
    memset(head,-1,sizeof(head));
    maxedge=-1;
    for(int i=1;i<=m;i++)
    {
        int x,y,c;
        scanf("%d%d%d",&x,&y,&c);
        addedge(x,y,c);
    }
    S=1,T=n;
}
int dis[maxn],cur[maxn];
bool bfs(int S,int T)
{
    queue <int> que;
    memset(dis,0,sizeof(dis));
    que.push(S); dis[S]=1; // initialize as 1 at first to avoid failed judgement of visited!!!
    while(!que.empty())
    {
        int u=que.front(); que.pop();
        if(u==T) return true;
        for(int i=head[u];~i;i=edge[i].next)
        {
            int v = edge[i].to;
            if(dis[v] || !edge[i].cap) continue;
            dis[v] = dis[u] + 1;
            que.push(v);
        }
    }
    return false;
}
int dfs(int u,int T,int minflow)
{
    if(u==T || !minflow) return minflow;
    int flow=0;
    for(int &i=cur[u];~i;i=edge[i].next)
    {
        int v=edge[i].to,f;
        if(dis[u]+1==dis[v] && (f=dfs(v,T,min(minflow,edge[i].cap))))
            edge[i].cap-=f,edge[i^1].cap+=f,flow+=f,minflow-=f;
        if(!minflow) break;
    }
    return flow;
}
int dinic(int S,int T)
{
    int ret=0;
    while(bfs(S,T))
    {
        memcpy(cur,head,sizeof(head));
        ret += dfs(S,T,INF);
    }
    return ret;
}
int main()
{
    freopen("grass.in","r",stdin);
    freopen("grass.out","w",stdout);
    init();
    printf("%d",dinic(S,T));
    return 0;
}