1273 Drainage Ditches 模板题

题意：

Discuss那里有。好强悍。。

题目描述 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 41 2 401 4 202 4 202 3 303 4 10

样例输出 Sample Output 50

代码：

//author: CHC//First Edit Time:2014-08-03 08:21//Last Edit Time:2014-08-03 08:56#include <iostream>#include <cstdio>#include <cstring>#include <cmath>#include <set>#include <vector>#include <map>#include <queue>#include <set>#include <algorithm>using namespace std;#define MAXN 3000typedef long long LL;struct Edge{    int from,to,ci,next;    Edge(){}    Edge(int _f,int _t,int _c,int _n):from(_f),to(_t),ci(_c),next(_n){}}e[MAXN];int head[MAXN],tot,n,m;int dis[MAXN],sta[MAXN],top,cur[MAXN];inline void init(){    memset(head,-1,sizeof(head));    tot=0;}inline void AddEdge(int u,int v,int ci1,int ci2=0){    e[tot]=Edge(u,v,ci1,head[u]);    head[u]=tot++;    e[tot]=Edge(v,u,ci2,head[v]);    head[v]=tot++;}inline bool bfs(int st,int et){    memset(dis,0,sizeof(dis));    queue <int> q;    q.push(st);    dis[st]=1;    while(!q.empty()){        int now=q.front();        q.pop();        for(int i=head[now];i!=-1;i=e[i].next){            if(!dis[e[i].to]&&e[i].ci){                dis[e[i].to]=dis[now]+1;                if(e[i].to==et)return true;                q.push(e[i].to);            }        }    }    return false;}LL Dinic(int st,int et){    LL ans=0;    while(bfs(st,et)){        int u=st,i;        memcpy(cur,head,sizeof(head));        top=0;        while(1){            if(u==et){                int minn=e[sta[0]].ci,pos=0;                for(i=1;i<top;i++){                    if(minn>e[sta[i]].ci){                        minn=e[sta[i]].ci;                        pos=i;                    }                }                for(int i=0;i<top;i++){                    e[sta[i]].ci  -=minn;                    e[sta[i]^1].ci+=minn;                }                ans+=minn;                top=pos;                u=e[sta[top]].from;            }            for(i=cur[u];i!=-1;cur[u]=i=e[i].next)                if(e[i].ci&&dis[u]+1==dis[e[i].to])break;            if(cur[u]!=-1){                sta[top++]=cur[u];                u=e[cur[u]].to;            }            else {                if(top==0)break;                dis[u]=0;                u=e[sta[--top]].from;            }        }    }    return ans;}int main(){    while(~scanf("%d%d",&n,&m)){        init();        for(int i=0,x,y,v;i<n;i++){            scanf("%d%d%d",&x,&y,&v);            AddEdge(x,y,v);        }        LL ans=Dinic(1,m);        printf("%I64d\n",ans);    }    return 0;}/*5 41 2 401 4 202 4 202 3 303 4 10*/