最大流模版

最新推荐文章于 2020-03-09 16:33:28 发布

yyffw

最新推荐文章于 2020-03-09 16:33:28 发布

阅读量264

点赞数

CC 4.0 BY-SA版权

分类专栏：模板

本文链接：https://blog.youkuaiyun.com/yyffw/article/details/47955253

模板专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一种解决最大流问题的方法，并通过POJ1459题目进行实例演示。采用C++实现边增广路径寻找，通过不断调整残余网络来求得从源点到汇点的最大流量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

poj1459

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<vector>
using namespace std;
typedef long long ll;
#define mem(name,value) memset(name,value,sizeof(name))
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
const int maxn=210,maxm=110*110,inf=0x3f3f3f3f;

struct Side
{
    int to,next,w;
};

struct Max_flow
{
    int node[maxn],dis[maxn],gap[maxn],top,cnt;
    Side side[maxm*2];

    void add_side(int u,int v,int w,int rw){
        //cout<<u<<" "<<v<<" "<<w<<endl;
        side[top]=(Side){v,node[u],w}; node[u]=top++;
        side[top]=(Side){u,node[v],rw}; node[v]=top++;
    }

    void init(){
        mem(node,-1); top=0;
        mem(dis,0); mem(gap,0);
    }

    int get_flow(int u,int flow,int s,int t){
        if(u==t)return flow;
        int ans=0;
        for(int i=node[u];i!=-1;i=side[i].next){
            int v=side[i].to,w=side[i].w;
            if(dis[u]>dis[v]&&w){
                int f=get_flow(v,min(flow-ans,w),s,t);
                ans+=f;
                side[i].w-=f;
                side[i^1].w+=f;
                if(ans==flow) return ans;
            }
        }
        if(!(--gap[dis[u]])) dis[s]=cnt+2;
        gap[++dis[u]]++;
        return ans;
    }

    int get_maxflow(int s,int t,int nn){
        cnt=gap[s]=nn; int ans=0;
        while(dis[s]<cnt) ans+=get_flow(s,inf,s,t);
        return ans;
    }
}mf;

int n,m,np,nc;

int main()
{
    while(cin>>n>>np>>nc>>m){
        mf.init();
        int u,v,z;
        while(m--){
            scanf("%*[^(](%d,%d)%d",&u,&v,&z);
            mf.add_side(u,v,z,0);
        }
        for(int i=1;i<=np;i++){
            scanf("%*[^(](%d)%d",&u,&z);
            mf.add_side(n,u,z,0);
        }
        for(int i=1;i<=nc;i++){
            scanf("%*[^(](%d)%d",&u,&z);
            mf.add_side(u,n+1,z,0);
        }
        int ans=mf.get_maxflow(n,n+1,n+2);
        printf("%d\n",ans);
    }
}