POJ 1273 Drainage Ditches——最大流

题目传送门
这道题就是裸的最大流，因为各路大神强烈推荐Dinic算法，所以我舍弃EK算法。

Dinic 算法

Dinic也是一种增广路算法,通过对残量网络进行分层之后用DFS求增广路直到无法增广。
层次图：层次就是起点到每个点的最短路径的长度，只走层与层之间的边，显然可以用BFS求层次图。

时间复杂度

对于n个点的层次图最多只有n层，所以最多重新建图n次，然后每次增广的效率为O(nm),s所以Dinic的效率是O(n^2*m)，然而实际效率十分玄学。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=205,maxm=405,MAXINT=2147483647;
int n,m,ans,tot,que[maxn],hed,til,dis[maxn],B[maxn];
int lnk[maxn],son[maxm],nxt[maxm],cap[maxm],flow[maxm];
bool vis[maxn];
void add(int x,int y,int z){
    nxt[tot]=lnk[x];lnk[x]=tot;son[tot]=y;cap[tot]=z;flow[tot]=0;tot++;
}
bool BFS(){
    memset(vis,0,sizeof(vis));
    hed=0,til=1;que[1]=1;vis[1]=1;dis[1]=1;
    while(hed!=til){
        int x=que[++hed];
        for (int j=lnk[x];j!=-1;j=nxt[j])
        if (!vis[son[j]]&&cap[j]>flow[j]){
            que[++til]=son[j];vis[son[j]]=1;
            dis[son[j]]=dis[x]+1;
        }
    }
    return vis[n];
}
int DFS(int x,int mi){
    if (x==n||mi==0) return mi;
    int num=0;
    for (int &j=B[x];j!=-1;j=nxt[j])
    if (dis[x]+1==dis[son[j]]){
        int now=DFS(son[j],min(mi,cap[j]-flow[j]));
        if (now){
            flow[j]+=now;flow[j^1]-=now;
            num+=now;mi-=now;
            if (!mi) break;
        }
    }
    return num;
}
void Dinic(){
    int ans=0;
    while (BFS()){
        memcpy(B,lnk,sizeof(B));
        ans+=DFS(1,MAXINT);
    }
    printf("%d\n",ans);
}
int main(){
    freopen("exam.in","r",stdin);
    freopen("exam.out","w",stdout);
    while (scanf("%d%d",&m,&n)==2){
        memset(lnk,255,sizeof(lnk));
        tot=0;
        for (int i=1;i<=m;i++){
            int x,y,z;scanf("%d%d%d",&x,&y,&z);
            add(x,y,z);add(y,x,0);
        }
        Dinic();
    }
    return 0;
}