poj 3469 Dual Core CPU 最小割

最大流算法实现

最新推荐文章于 2018-11-29 17:47:08 发布

KIDGINBROOK

最新推荐文章于 2018-11-29 17:47:08 发布

阅读量543

点赞数

CC 4.0 BY-SA版权

分类专栏：网络流

本文链接：https://blog.youkuaiyun.com/KIDGIN7439/article/details/17285237

网络流专栏收录该内容

51 篇文章

订阅专栏

本文介绍了一个使用C++实现的最大流算法。该程序通过构建图结构并应用增广路径的方法来求解从源点到汇点的最大流量。文章包含了完整的源代码，展示了如何初始化图、添加边、更新边容量及进行深度优先搜索等关键步骤。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=20000+50,inf=100000000;
int n,m;
int level[maxn],que[maxn*10];
int head[maxn*3],lon;
int min(int a,int b)
{
	if(a<b) return a;
	else return b;
}
struct
{
    int next,to,c;
}e[2000000];
void edgeini()
{
    memset(head,-1,sizeof(head));
    lon=-1;
}
void edgemake(int from,int to,int c)
{
    e[++lon].c=c;
    e[lon].to=to;
    e[lon].next=head[from];
    head[from]=lon;
}
void make(int from,int to,int c)
{
    edgemake(from,to,c);
    edgemake(to,from,0);
}

bool makelevel(int s,int t)
{
    memset(level,0,sizeof(level));
    int front=1,end=0;
    que[++end]=s;
    level[s]=1;
    while(front<=end)
    {
        int u=que[front++];
        if(u==t) return true;
        for(int k=head[u];k!=-1;k=e[k].next)
        {
            int v=e[k].to;
            if(!level[v]&&e[k].c)
            {
                que[++end]=v;
                level[v]=level[u]+1;
            }
        }
    }
    return false;
}

int dfs(int now,int t,int maxf)
{
    if(now==t) return maxf;
    int ret=0;
    for(int k=head[now];k!=-1;k=e[k].next)
    {
        int u=e[k].to;
        if(level[u]==level[now]+1&&e[k].c)
        {
            int f=dfs(u,t,min(e[k].c,maxf-ret));
            e[k].c-=f;
            e[k^1].c+=f;
            ret+=f;
            if(ret==maxf) return ret;
        }
    }
    return ret;
}

int maxflow(int s,int t)
{
    int ret=0;
    while(makelevel(s,t))
    {
        ret+=dfs(s,t,inf);
    }
    return ret;
}
int main()
{
	int i,j;
	int a,b,w;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		edgeini();
		for(i=1;i<=n;i++)
		{
			scanf("%d%d",&a,&b);
			make(0,i,a);
			make(i,n+1,b);
		}
		for(i=1;i<=m;i++)
		{
			scanf("%d%d%d",&a,&b,&w);
			make(a,b,w);
			make(b,a,w);
		}
		printf("%d\n",maxflow(0,n+1));
	}
	return 0;
}