網絡流isap

對於網絡流，在高端局是十分有用的。
isap則是除了預留推進以外最好的網絡流了（自己覺得），加入gap優化後，跑出來的效率還不錯。
至於isap的基本操作則是利用了d[node]=d[pre]+1時爲最優這個性質。也是通過增廣的方法來進行算法。
而gap優化，則是在更新層次圖時，記錄了節點的個數。如果節點個數爲零，則不會存在增廣路了，即可以退出了。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#define For(aa,bb,cc) for(int aa=bb;aa<=cc;++aa)
using namespace std;
const int maxn=1010,inf=1e10;
int n,m,s,t;
int be[maxn],ne[maxn],w[maxn],to[maxn],e;
int d[maxn],gap[maxn];

void add(int x,int y,int z){
    to[++e]=y,ne[e]=be[x],be[x]=e,w[e]=z;
}

int isap(int node,int flow){
    if(node==t) return flow;
    int res=flow;
    for(int i=be[node];i;i=ne[i]){
        int u=to[i];
        if(w[i]>0 && d[node]==d[u]+1){
            int minflow=isap(u,min(w[i],res));
            w[i]-=minflow;
            w[i^1]+=minflow;
            res-=minflow;
            if(!res) return flow;
        }
    }
    if(!(--gap[d[node]])) d[s]=n;
    ++gap[++d[node]];
    return flow-res;
}

void work(){
    e=1;
    scanf("%d%d",&m,&n);
    s=1,t=n;
    For(i,1,m){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z),add(y,x,0);
    }
    int maxflow=0;
    for(gap[0]=n;d[s]<n;){
        maxflow+=isap(s,inf);
    }
    printf("%d\n",maxflow);
}

int main(){
    work();
    return 0;
}