POJ1273 最大流模板

之前自己写的，以后当一个模板的基础吧，不管是最大流，最小割，二分图匹配

下面是POJ1273的一个裸题..

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <queue>
using namespace std;
struct node{
    int to,cap,rev;
    node(int _to,int _cap,int _rev):to(_to),cap(_cap),rev(_rev){}
};
const int maxn=1005;
vector<node> G[maxn];
int s,t;
void add(int u,int v,int cap){
    G[u].push_back(node(v,cap,G[v].size()));
    G[v].push_back(node(u,0,G[u].size()-1));
}
bool vis[maxn];//这个可以不要，level数组可以起两个作用,记录层数和是否被访问过，类似于dp数组
int level[maxn],iter[maxn];
void bfs(){
    memset(level,-1,sizeof(level));
    queue<int> q;
    q.push(s);int i;level[s]=0;
    while(!q.empty()){
        int u=q.front();
        q.pop();
        for(i=0;i<G[u].size();++i){
            node& e=G[u][i];
            if(e.cap>0&&level[e.to]<0){
                level[e.to]=level[u]+1;
                q.push(e.to);                
            }
        }
    }
    //需要得到所有顶点的level信息，一方面起vis的作用，另一方面起到转移的作用
    //而且考虑到可能同时存在多个最短增广路
}
int dfs(int v,int t,int f){
    //printf("v:%d t:%d f:%d
",v,t,f);
    if(v==t) return f;
    for(int &i=iter[v];i<G[v].size();++i){
        node &e=G[v][i];
        if(e.cap>0&&level[e.to]>level[v]){
            int d=dfs(e.to,t,min(e.cap,f));
            if(d>0){
                e.cap-=d;
                G[e.to][e.rev].cap+=d;
                return d;
            }
        }        
    }
    return 0;//这里忘记return 0了...
}
const int INF=~0u>>1;
int maxflow(){
    int flow=0,f;
    for(;;){
        bfs();
        if(level[t]<0) return flow;
        memset(iter,0,sizeof(iter));
        while((f=dfs(s,t,INF))>0) flow+=f; 
    }
}
//直觉上 容量种数越少越快
int main(){
    //freopen("testMaxFlow.txt","r",stdin);
    int n,i,u,v,cap,m;
    while(~scanf("%d%d",&m,&n)){
        s=1;t=n;
        for(int i=1;i<=n;++i) G[i].clear();//多组数据要清空边
        for(i=0;i<m;++i){
            scanf("%d%d%d",&u,&v,&cap);
            add(u,v,cap);
        }
        printf("%d
",maxflow());
    }
    return 0;
}