hdu 3549 a flow problem 的多种解法

最新推荐文章于 2022-11-26 15:37:20 发布
最新推荐文章于 2022-11-26 15:37:20 发布
阅读量223 收藏 1
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/chenhongwei6447/article/details/48248055
版权

网路流裸题，用来测试模板吧。

1.EK－1

//  Created by Chenhongwei in 2015.
//  Copyright (c) 2015 Chenhongwei. All rights reserved.

#include"iostream"
#include"cstdio"
#include"cstdlib"
#include"cstring"
#include"climits"
#include"queue"
#include"cmath"
#include"map"
#include"set"
#include"stack"
#include"vector"
#include"sstream"
#include"algorithm"
using namespace std;
typedef long long ll;
const int inf=1e8;
const int maxn=1100;
struct node
{
	int to,cap,rev;
};
int n,m;
int used[maxn];
vector<node> g[maxn];
int dfs(int s,int t,int f)
{
	if(s==t)
		return f;
	used[s]=1;
	for(int i=0;i<g[s].size();i++)
	{
		node& e=g[s][i];
		if(!used[e.to]&&e.cap>0)
		{
			int d=dfs(e.to,t,min(f,e.cap));
			if(d>0)
			{
				e.cap-=d;
				g[e.to][e.rev].cap+=d;
				return d;
			}
		}
	}
	return 0;

}
int max_flow(int s,int t)
{
	int flow=0;
	while(1)
	{
		memset(used,0,sizeof used);
		int f=dfs(s,t,inf);
		if(f==0)
			return flow;
		flow+=f;
	}
}
int main()
{	
	//ios::sync_with_stdio(false);
	// freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	int T;
	scanf("%d",&T);
	for(int z=1;z<=T;z++)
	{
		scanf("%d%d",&n,&m);
		memset(g,0,sizeof g);
		for(int i=0;i<m;i++)
		{
			int u,v,w;
			scanf("%d%d%d",&u,&v,&w);
			g[u].push_back((node){v,w,g[v].size()});
			g[v].push_back((node){u,0,g[u].size()-1});
		}
		printf("Case %d: %d\n",z,max_flow(1,n));

	}
	return 0;
}

  EK－2 

//  Created by Chenhongwei in 2015.
//  Copyright (c) 2015 Chenhongwei. All rights reserved.

#include "iostream"
#include "cstdio"
#include "cstdlib"
#include "cstring"
#include "climits"
#include "queue"
#include "cmath"
#include "map"
#include "set"
#include "stack"
#include "vector"
#include "sstream"
#include "algorithm"
using namespace std;
typedef long long ll;
const int inf=1e8;
const int maxn=1100;
struct node
{
	int from,to,cap,rev;
	node(int from,int to,int cap):from(from),to(to),cap(cap) {}
};
vector<node> edge;
vector<int >g[maxn];
int pre[maxn],a[maxn];
int bfs(int s,int t)
{
	queue<int>q;
	q.push(s);
	memset(a,0,sizeof a);
	a[s]=inf;
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=0;i<g[u].size();i++)
		{
			node& e=edge[g[u][i]];
			if(!a[e.to]&&e.cap>0)
			{
				pre[e.to]=g[u][i];
				a[e.to]=min(a[u],e.cap);
				q.push(e.to);
			}
		}
		if(a[t])break;
	}
	if(!a[t]) return 0;
	for(int u=t;u!=s;u=edge[pre[u]].from)
	{
		edge[pre[u]].cap-=a[t];
		edge[pre[u]^1].cap+=a[t];
	}
	return a[t];

}
int max_flow(int s,int t)
{
	int flow=0;
	while(1)
	{
		memset(pre,0,sizeof pre);
		int f=bfs(s,t);
		if(f==0)
			return flow;
		flow+=f;
	}
}
int main()
{	
	//ios::sync_with_stdio(false);
	// freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	int n,m,T;
	scanf("%d",&T);
	for(int z=1;z<=T;z++)
	{
		scanf("%d%d",&n,&m);
		for(int i=0;i<=n;i++)
			g[i].clear();
		edge.clear();
		for(int i=0;i<m;i++)
		{
			int u,v,w;
			scanf("%d%d%d",&u,&v,&w);
			edge.push_back(node(u,v,w));
			edge.push_back(node(v,u,0));
			int t=edge.size();
			g[u].push_back(t-2);
			g[v].push_back(t-1);
		}
		printf("Case %d: %d\n",z,max_flow(1,n));

	}
	return 0;
}

SAP

//  Created by Chenhongwei in 2015.
//  Copyright (c) 2015 Chenhongwei. All rights reserved.

#include"iostream"
#include"cstdio"
#include"cstdlib"
#include"cstring"
#include"climits"
#include"queue"
#include"cmath"
#include"map"
#include"set"
#include"stack"
#include"vector"
#include"sstream"
#include"algorithm"
using namespace std;
typedef long long ll;
const int inf=1e8;
const int maxn=20;
int n;
int d[maxn];
int g[maxn][maxn];
bool tag;
int dfs(int cur,int aug,int s,int t)
{
	if(cur==t)
	{
		tag=true;
		return aug;
	}
	for(int i=1;i<=n;i++)
	{
		if(g[cur][i]>0&&d[cur]==d[i]+1)
		{
			int f=dfs(i,min(aug,g[cur][i]),s,t);
			if(d[s]>=n)
				return 0;
			if(tag)
			{
				g[cur][i]-=f;
				g[i][cur]+=f;
				return f;
			}
		}
	}
	d[cur]++;
	return 0;

}
int sap(int s,int t)
{
	int flow=0;
	memset(d,0,sizeof d);
	while(d[s]<n)
	{
		tag=false;
		flow+=dfs(s,inf,s,t);
	}
	return flow;
}
int main()
{	
	//ios::sync_with_stdio(false);
	// freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	int T,m;
	scanf("%d",&T);
	for(int z=1;z<=T;z++)
	{
		scanf("%d%d",&n,&m);
		memset(g,0,sizeof g);
		for(int i=0;i<m;i++)
		{
			int u,v,w;
			scanf("%d%d%d",&u,&v,&w);
			g[u][v]+=w;
		}
		printf("Case %d: %d\n",z,sap(1,n));

	}
	return 0;
}

DINIC 

//  Created by Chenhongwei in 2015.
//  Copyright (c) 2015 Chenhongwei. All rights reserved.

#include"iostream"
#include"cstdio"
#include"cstdlib"
#include"cstring"
#include"climits"
#include"queue"
#include"cmath"
#include"map"
#include"set"
#include"stack"
#include"vector"
#include"sstream"
#include"algorithm"
using namespace std;
typedef long long ll;
const int inf=1e8;
const int maxn=1100;
struct node
{
	int to,cap,rev;
};
int n,m;
int cur[maxn],d[maxn];
vector<node> g[maxn];
void bfs(int s)
{
	memset(d,-1,sizeof d);
	d[s]=0;
	queue<int >q;
	q.push(s);
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=0;i<g[u].size();i++)
		{
			node& e=g[u][i];
			if(e.cap>0&&d[e.to]<0)
			{
				d[e.to]=d[u]+1;
				q.push(e.to);
			} 		
		}
	}
}
int dfs(int s,int t,int f)
{
	if(s==t)
		return f;
	for(int i=cur[s];i<g[s].size();i++)
	{
		cur[s]=i;
		node& e=g[s][i];
		if(e.cap>0&&d[e.to]>d[s])
		{
			int d=dfs(e.to,t,min(f,e.cap));
			if(d>0)
			{
				e.cap-=d;
				g[e.to][e.rev].cap+=d;
				return d;
			}
		}
	}
	return 0;

}
int dinic(int s,int t)
{
	int flow=0;
	while(1)
	{
		bfs(s);
		if(d[t]<0)
			return flow;
		memset(cur,0,sizeof cur);
		int f;
		while((f=dfs(s,t,inf))>0)
			flow+=f;
	}
}
int main()
{	
	//ios::sync_with_stdio(false);
	// freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	int T;
	scanf("%d",&T);
	for(int z=1;z<=T;z++)
	{
		scanf("%d%d",&n,&m);
		memset(g,0,sizeof g);
		for(int i=0;i<m;i++)
		{
			int u,v,w;
			scanf("%d%d%d",&u,&v,&w);
			g[u].push_back((node){v,w,g[v].size()});
			g[v].push_back((node){u,0,g[u].size()-1});
		}
		printf("Case %d: %d\n",z,dinic(1,n));

	}
	return 0;
}

ISAP 

//  Created by Chenhongwei in 2015.
//  Copyright (c) 2015 Chenhongwei. All rights reserved.

#include "iostream"
#include "cstdio"
#include "cstdlib"
#include "cstring"
#include "climits"
#include "queue"
#include "cmath"
#include "map"
#include "set"
#include "stack"
#include "vector"
#include "sstream"
#include "algorithm"
using namespace std;
typedef long long ll;
const int inf=1e8;
const int maxn=1100;
struct node
{
	int from,to,cap;
	node(int from,int to,int cap):from(from),to(to),cap(cap) {}
};
vector<node> edge;
vector<int >g[maxn];
int p[maxn],cur[maxn];
int num[maxn],d[maxn];
int n;
void addedge(int u,int v,int w)
{
	edge.push_back(node(u,v,w));
	edge.push_back(node(v,u,0));
	int t=edge.size();
	g[u].push_back(t-2);
	g[v].push_back(t-1);

}
void BFS(int s,int t)
{
	memset(d,-1,sizeof d);
	queue<int>q;
	q.push(t);
	d[t]=0;
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=0;i<g[u].size();i++)
		{
			node& e=edge[g[u][i]^1];
			if(d[e.from]<0&&e.cap>0)
			{
				d[e.from]=d[u]+1;
				q.push(e.from);
			}
		}
	}
}
int augment(int s,int t)
{
	int x=t,ans=inf;
	while(x!=s)
	{
		node& e=edge[p[x]];
		ans=min(ans,e.cap);
		x=edge[p[x]].from;
	}
	x=t;
	while(x!=s)
	{
		edge[p[x]].cap-=ans;
		edge[p[x]^1].cap+=ans;
		x=edge[p[x]].from;
	}
	return ans;
}
int ISAP(int s,int t)
{
	int flow=0;
	BFS(s,t);
	memset(num,0,sizeof num);
	for(int i=0;i<n;i++)
		num[d[i]]++;
	int x=s;
	memset(cur,0,sizeof cur);
	while(d[s]<n)
	{
		if(x==t)
		{
			flow+=augment(s,t);
			x=s;
		}
		bool ok=0;
		for(int i=cur[x];i<g[x].size();i++)
		{
			node& e=edge[g[x][i]];
			if(e.cap>0&&d[x]==d[e.to]+1)
			{
				ok=1;
				p[e.to]=g[x][i];
				cur[x]=i;
				x=e.to;
				break;
			}
		}
		if(!ok)
		{
			int m=n-1;
			for(int i=0;i<g[x].size();i++)
			{
				node& e=edge[g[x][i]^1];
				if(e.cap>0)
					m=min(m,d[e.to]);
			}
			if(--num[d[x]]==0)
				break;
			d[x]=m+1;
			num[d[x]]++;
			cur[x]=0;
			if(x!=s)
				x=edge[p[x]].from;
		}
	}
	return flow;
}
int main()
{	
	//ios::sync_with_stdio(false);
	freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	int T,m,u,v,w;
	scanf("%d",&T);
	for(int z=1;z<=T;z++)
	{
		scanf("%d%d",&n,&m);
		for(int i=0;i<=n;i++)
			g[i].clear();
		edge.clear();
		for(int i=0;i<m;i++)
		{
			scanf("%d%d%d",&u,&v,&w);
			addedge(u,v,w);
		}
		// s=1,t=n;
		printf("Case %d: %d\n",z,ISAP(n-1,n));

	}
	return 0;
}




chenhongwei6447
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
hdu 3549 Flow Problem(简单网络流Dinic）
stationinthemind的博客
11-25 451
/不是很懂模板，只能修改点迎合模板。分析：很简明的问题，就是求解最大流。
流问题Flow Problem（网络最大流）- HDU 3549
ACM算法日常
12-26 5524
        网络最大流问题属于算法 里面较难的问题，因为牵涉的概念比较多，这一篇可能需要你花比较多的时间去理解，除了看这个，最好能多参考别的书籍或者文章进行比较学习，不然可能容易产生理解的偏差。           另外本公众号并不打算讲解过于复杂的问题，网络流问题已经严重超过了最初设想的五分钟限制，不过由于是第一篇网络流的题目，可以作为后面相关问题的基础，姑且多写点。如果你觉得一次难以看...
hdu 3549 a flow problem 的多种解法
liuliuhelingdao的博客
11-17 85
网路流裸题，用来测试模板吧。
【网络流第三弹】HDU 3549——Flow Problem(dinic解法)
apple_51426592的博客
11-26 82
/ 1.设置上限MAXN和MAXN 2.初始化数值src, des, nNode 3.init_adj() & ex_addEdge() 4.dinic()求有向图的网络流，DINIC模板过，稍作修改即可。DINIC还是非递归版比较好的，省时间。
HDU 3549 Flow Problem（最大流裸题，EK解法）
m0_72431373的博客
11-19 92
EK现在用起来已经比较熟练了，比较简单的题目，直接就可以解决，但发现EK，其实效率有些低，是O（V*E^2）的，搞不好，还会超时。今天，听学长讲，dinic用的比较多，sap学起来有一定难度，接下来抓紧学一下dinic！注意点：EK两点之间只能有一条正向边，而顶点最多只有15个，边却最多有1000条，15*14/2
[毕业设计]VB+ACCESS高校班级日常管理系统(源代码+论文).zip
08-07
[毕业设计]VB+ACCESS高校班级日常管理系统(源代码+论文)
基于Misty1算法的加密软件(JAVA)的实现模块项目(源代码+论文).zip
08-07
JAVA; 毕业设计；MVC模式; JSP; SQL Server 2000 数据库管理系统; J2EE
py源码实例自动办公用Python在Excel中查找并替换数据
最新发布
08-07
python，py源码实例自动办公用Python在Excel中查找并替换数据
Rust语言的环境搭建
08-07
Rust语言的环境搭建，安装最新版的 Rust 编译工具和 Visual Studio Code。
架构规划方法.pptx
08-07
架构规划方法.pptx
K-means聚类算法聚类算法
08-07
K-means聚类算法是一种典型的无监督机器学习算法，主要用于数据聚类。它的目标是将相似的数据点归到同一个簇中，而将不相似的数据点归到不同的簇中。
阳光酒店管理系统项目模块(javaapplet+SQL).zip
08-07
阳光酒店管理系统项目模块(javaapplet+SQL)
《嵌入式操作系统内核调度》一书的源码。本书不是讲解操作系统如何使用，而是使用汇编和C语言从无到有实现操作系统的多个功能。.zip
08-07
优质项目，资源经过严格测试可直接运行成功且功能正常的情况才上传，可轻松copy复刻，拿到资料包后可轻松复现出一样的项目。 本人系统开发经验充足，有任何使用问题欢迎随时与我联系，我会及时为你解惑，提供帮助。 【资源内容】：包含完整源码+工程文件+说明（若有），项目具体内容可查看下方的资源详情。 【附带帮助】： 若还需要相关开发工具、学习资料等，我会提供帮助，提供资料，鼓励学习进步。 【本人专注计算机领域】： 有任何使用问题欢迎随时与我联系，我会及时解答，第一时间为你提供帮助，CSDN博客端可私信，为你解惑，欢迎交流。 【适合场景】： 相关项目设计中，皆可应用在项目开发、毕业设计、课程设计、期末/期中/大作业、工程实训、大创等学科竞赛比赛、初期项目立项、学习/练手等方面中 可借鉴此优质项目实现复刻，也可以基于此项目进行扩展来开发出更多功能 【无积分此资源可联系获取】 # 注意 1. 本资源仅用于开源学习和技术交流。不可商用等，一切后果由使用者承担。 2. 部分字体以及插图等来自网络，若是侵权请联系删除。积分/付费仅作为资源整理辛苦费用。
优秀班组长的角色定位.pptx
08-07
优秀班组长的角色定位.pptx
毕业设计，基于ASP+ACCESS开发的视频点播系统，内含完整源代码，数据库，毕业论文
08-07
毕业设计，基于ASP+ACCESS开发的视频点播系统，内含完整源代码，数据库，毕业论文 作为计算机发展最迅速的领域之一的数据库技术,已经形成了一整套独有的理论,并广泛地应用于人们的生产和生活中。数据库技术与网站的结合是当今Web技术的一个热点。有了数据库的支持,可以扩展网页的功能,可以方便地设计出交互式页面,可以构造功能强大的后台管理系统,可以为网站的更新、维护提供极大的方便。因此,作为网络开发者或管理者,数据库知识是必不可少的。 视频流媒体服务一直被认为是互联网服务中的一个热点领域，但是由于技术、内容以及带宽等方面的限制，一直没有实现预料中的发展。目前，视频压缩技术已经取得了重大的进步，视频流媒体的传输速度、品质以及成本控制都得到了很大的改善，为视频媒体传输创造了有利条件，同时宽带网络的普及也使得多数家庭用户具备了享受网络视频服务的条件，据估计，目前约有2000万到3000个家庭在通过宽带接入的方式上网。In-Stat/MDR预计，到2007年，这个数字将超过1亿3000万。InStat/MDR预计，视频订阅服务最繁荣的市场将出现在亚洲，预计2007年亚洲市场视频媒体订阅服务的营业
BEVFORMER相关视频
08-07
BEVFORMER相关视频
毕业设计，基于ASP+SqlServer开发的旅游网站，内含完整源代码，数据库，毕业论文
08-07
毕业设计，基于ASP+SqlServer开发的旅游网站，内含完整源代码，数据库，毕业论文 旅行社是服务行业，旅游网也是服务行业。只有建立与旅游网相适应的服务体系，并保证其服务质量，才能使旅游网站落到实处，也才能真正使旅游望发挥作用。 旅游网站全国已有650多家，而无论是企业网站还是门户网站，要想活得好起来，就得在服务内容、范围、品质上有所突破和发展。 旅游企业是服务行业，旅游网同样也是服务行业，“得服务者得天下”是网站的竞争关键。旅游网的生存取决于它在服务内容、范围、品质上的含金量。目前的许多企业旅游网站和许多旅游企业一样，除存在经营定位模糊的弱点外，还存在着“单兵作战”所带来的预订内容范围过于狭窄的通病，因此难以产生商务效应和边际效应。80%的旅游网站以发布地方信息本社信息为主，服务范围也局限于区域性，这样互联的优势难以体现。造成这种情况的主要原因在于互联网的区域经营与提供地方化服务的这一对矛盾，服务的延伸速度远远跟不上互联网的发展和需求速度。对网站来说信息覆盖面减小，也会影响网上服务品质的下降，影响地方信息的可操作性和可延伸性。 从个人兴趣及毕业设计的实际情况出发,经过一系列的考
JAVA本地监听与远程端口扫描项目模块(源代码+论文).zip
08-07
JAVA本地监听与远程端口扫描项目模块(源代码+论文)
华三无线控制器初始化简单步骤
08-07
华三无线控制器初始化简单步骤
hdu4310 java解法
11-11
以下是hdu4310的Java解法： ```java import java.util.*; import java.io.*; public class Main { static int MAXN = 100010; static int MAXM = 200010; static int INF = 0x3f3f3f3f; static int n, m, s, ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值