Dinic

//建立反向边是因为贪心不最优,给自己反悔的机会
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<queue> #include<vector> #include<math.h> using namespace std; #define M 9499999 int c[10005][10005]; int dep[10001]; int n,m,s,t; int dfs(int s,int t,int f) { if(s==t) return f; int d; for(int v=1;v<=n;v++) if(c[s][v]>0&&dep[v]==dep[s]+1&&(d=dfs(v,t,min(f,c[s][v])))) { c[s][v]-=d; c[v][s]+=d; return d; } return 0; } int bfs(int s,int t) { queue<int>q;q.push(s); memset(dep,-1,sizeof(dep)); dep[s]=0; while(!q.empty()){ int u=q.front();q.pop(); for(int v=1;v<=m;v++) if(c[u][v]>0&&dep[v]==-1) { dep[v]=dep[u]+1; q.push(v); } } return dep[t]!=-1; } int dinic() { int flow=0,f; while(bfs(s,t)) while(1){ f=dfs(s,t,M); if(!f) break; flow+=f; } return flow; } int main() { scanf("%d%d%d%d",&n,&m,&s,&t); for(int i=1,u,v,w;i<=m;i++) { scanf("%d%d%d",&u,&v,&w); c[u][v]+=w; } cout<<dinic(); return 0; }

 链表的(较快)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<math.h>
using namespace std;
#define M    94999999
int h[10002],nex[200002],na[200002],w[200002];
int dep[10001];
int n,m,s,t;
int dfs(int s,int t,int f)//当前点,目标,可用流量
{
    if(s==t)    return f;
    int d;
    for(int i=h[s];i;i=nex[i])
    {
        int v=na[i];
        if(w[i]>0&&dep[v]==dep[s]+1&&(d=dfs(v,t,min(f,w[i]))))
        {
            w[i]-=d;
            w[i^1]+=d;
            return d;
        }
    }    
}
int bfs(int s,int t)//分层
{
    queue<int>q;q.push(s);
    memset(dep,-1,sizeof(dep));
    dep[s]=0;
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=h[u];i;i=nex[i])
        {
            int v=na[i];
            if(w[i]>0&&dep[v]==-1)
            {
                dep[v]=dep[u]+1;
                q.push(v);
            }
        }
    }
    return dep[t]!=-1;
}
int dinic()
{
    int flow=0,f;
    while(bfs(s,t))
        while(1){
        f=dfs(s,t,M);
        if(!f)    break;
        flow+=f;
        }
    return flow;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&t);
    int tot=1; 
    for(int i=1,u,v,ww;i<=m;i++)
    {
        scanf("%d%d%d",&u,&v,&ww);
        na[++tot]=v;nex[tot]=h[u];h[u]=tot;w[tot]=ww;
        na[++tot]=u;nex[tot]=h[v];h[v]=tot;
    }
    cout<<dinic();
    return 0;
}

Dinic的核心是对搜索进行诱导(横向),降低复杂度。

posted @ 2017-04-21 16:03  浪矢-CL  阅读(248)  评论(0编辑  收藏  举报