网络最大流

http://blog.csdn.net/largecub233/article/details/72667961

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);
}

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int N=1e4+5;
struct cs{int to,nxt,v;}a[N*20];
int head[N],ll=1,cur[N];
int q[N],l,r;
int deep[N];
int n,m,S,E,x,y,z,ans;
void init(int x,int y,int z){a[++ll].to=y;a[ll].v=z;a[ll].nxt=head[x];head[x]=ll;}
bool bfs(){
    memset(deep,0,sizeof deep);
    q[1]=S;l=r=1;deep[S]=1;
    for(;r>=l;l++){
        int x=q[l];cur[x]=head[x];if(x==E)return 1;
        for(int k=head[x];k;k=a[k].nxt)
            if(!deep[a[k].to]&&a[k].v!=0)
                deep[a[k].to]=deep[x]+1,q[++r]=a[k].to;
    }return 0;
}
int dfs(int x,int now){
    if(x==E)return now;
    int ans=0;
    for(int &k=cur[x];k;k=a[k].nxt)
        if(deep[a[k].to]-1==deep[x]&&a[k].v){
            int temp=dfs(a[k].to,min(a[k].v,now));
            now-=temp; ans+=temp;
            a[k].v-=temp; a[k^1].v+=temp;
            if(!now)break;
        }return ans;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&S,&E);
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&x,&y,&z);
        init(x,y,z);init(y,x,0);
    }
    while(bfs())ans+=dfs(S,1e9);
    printf("%d",ans);
}