【模板】网络流模板

好吧我只会最大流

dinic

 2016.2.25：修改，原来那个有点问题，换了个stl_queue的广搜，改了几个人畜无害的变量名

 2016.2.25比上面那个晚些时候：司机把我虐成狗，poj3469教我做人，加了当前弧，还是没A掉，但是我对正确性有信心

const int inf=0x7fffffff;
struct edge{
    int to,cap,next;
}e[M*4]; 
int last[N],iter[N],ecnt=1;
void insert(int from,int to,int cap){
    e[++ecnt]=(edge){to,cap,last[from]};last[from]=ecnt;
    e[++ecnt]=(edge){from,0,last[to]};last[to]=ecnt;
}
int h[N];
bool bfs(int s,int t){
    memset(h,-1,sizeof(h));
    queue<int>q;
    q.push(s);h[s]=0;
    while(!q.empty()){
        int c=q.front();q.pop(); 
        for(int u=last[c];u;u=e[u].next){
            if(e[u].cap>0&&h[e[u].to]==-1){
                h[e[u].to]=h[c]+1;
                q.push(e[u].to);
            }
        }
    }
    return h[t]!=-1;
}
int dfs(int x,int t,int f){
    if(x==t)    return f;
    int w,used=0;
    for(int u=iter[x];u;u=e[u].next){
        if(h[e[u].to]==h[x]+1){
            w=dfs(e[u].to,t,min(f-used,e[u].cap));
            e[u].cap-=w;
            e[u^1].cap+=w;
            if(e[u].cap)    iter[x]=u;
            used+=w;
            if(used==f)    return f; 
        }
    }
    if(!used)h[x]=-1;
    return used;
}
int dinic(int s,int t){
    int flow=0;
    while(bfs(s,t)){
        for(int i=0;i<=t;i++)    iter[i]=last[i];
        flow+=dfs(s,t,inf);
    }
    return flow; 
}
int m,n;

 

RP++