网络最大流

https://www.luogu.org/problem/show?pid=3376
Dinic算法
http://blog.csdn.net/wall_f/article/details/8207595
给图:
这里写图片描述

这里写图片描述
这里写图片描述
这里写图片描述
这里写图片描述
就是这样,先搞分层图,在找增广路

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<cstring>
#include<string>
#include<cctype>
#include<queue>
#define Ll long long
using namespace std;
struct cs{
    int to,next,vv;
}a[1000000];
int head[100000],deep[100000],b[1000000];
int ll=1,ans,n,m,x,y,z,J,s,t;
bool p;
bool make(){//用dfs会慢 
    memset(deep,0,sizeof deep);
    int l=0,r=1;
    b[1]=s;deep[s]=1;
    while(r>l){
        l++;
        if(b[l]==t)return 1;
        for(int k=head[b[l]];k;k=a[k].next){
            if(!deep[a[k].to]&&a[k].vv){//当权值是0,说明不能走,也不用搞了 
                deep[a[k].to]=deep[b[l]]+1;
                b[++r]=a[k].to;
            }
        }
    }//其实我这样bfs数组靠开的大一点,不过无伤大雅 
    return 0;
}
void init(int x,int y,int z){
    ll++;
    a[ll].to=y;
    a[ll].vv=z;
    a[ll].next=head[x];
    head[x]=ll; 
}
int dfs(int x,int mi){//mi是当前的最大流,就是路径上最小的权 
    if(x==t)return mi;
    int k=head[x],ret=0;//ret即返回值,表示在x节点上时,目前已知有多少流是可行的 
    while(k){
        if(deep[x]+1==deep[a[k].to]&&a[k].vv){//层数差1那么dfs 
            int M=min(mi-ret,a[k].vv); 
            int f=dfs(a[k].to,M);
            a[k].vv-=f;
            a[k^1].vv+=f;//和下面的ll=1照应,在搞反向边 
            ret+=f;
            if(ret==mi)return ret;//ret不会超过mi,因为M=min(mi-ret,a[k].vv); 
        }
        k=a[k].next;
    }
    return ret;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&t);
    ll=1;//方便异或 
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&x,&y,&z);
        init(x,y,z);
        init(y,x,0);//建立反向边,一开始是0 
    }
    while(make())ans+=dfs(s,1e9);
    printf("%d",ans);
}
posted @ 2017-02-09 14:03  largecube233  阅读(186)  评论(0编辑  收藏  举报