草地排水-网络流dinic

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

提交的第一道网络流题，原理很简单，用bfs扫每条边中剩余流量最小的，同时记录路径，然后更新路径的流量。最重要的是构造了对应的反向边，给反向边对应的负流量，这样在更新的时候就能把上次留下但这次没有使用的流量释放出来。

代码如下：

 

#include<iostream>
#include<cstdio>
#include<cmath>
#include<queue>
#include<cstring>

using namespace std;

const int maxn=1001,INF=0xf777;
int n,m,map[maxn][maxn],vis[maxn],pre[maxn];

int read(){
    int x=0;char ch=getchar();
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9'){
        x=x*10+ch-'0';ch=getchar();
    }
    return x;
} 

int bfs(int s, int t){
    memset(pre,-1,sizeof(pre));
    memset(vis,0,sizeof(vis));
    vis[s]=1;
    int mi=INF;
    queue<int>q;
    q.push(s);
    while(!q.empty()){
        int cur=q.front();
        q.pop();
        if(cur==t)break;
        for(int i=1;i<=n;i++){
            if(!vis[i]&&map[cur][i]!=0){
                q.push(i);
                mi=min(mi,map[cur][i]);
                pre[i]=cur;
                vis[i]=1;
            }
        }
    }
    if(pre[t]==-1)return 0;
    return mi;
}

int update(int u,int f){
    while(pre[u]!=-1){
        map[pre[u]][u]-=f;
        map[u][pre[u]]+=f;
        u=pre[u];
    }
}

int EK(int s,int t){
    int maxf=0,newf;
    do{
        newf=bfs(s,t); 
        update(t,newf);
        maxf+=newf;
    }while(newf!=0);
    return maxf;
}

int main(){
    m=read();n=read();
    for(int i=1;i<=m;i++){
        int fr=read(),to=read(),dis=read();
        map[fr][to]+=dis;  //如果有重边就代数加（对于不同题不同） 
    }
    cout<<EK(1,n);
    return 0;
}