BZOJ最大流水题集

最新推荐文章于 2019-03-12 10:41:05 发布
最新推荐文章于 2019-03-12 10:41:05 发布
阅读量828 收藏
点赞数
分类专栏： bzoj 文章标签： bzoj oi
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/CDQZOIERS/article/details/52940089
版权

3993:二分答案+最大流

#include<bits/stdc++.h>

typedef unsigned char uchar;
typedef unsigned int uint;
typedef long long ll;
typedef unsigned long long ull;

#define xx first
#define yy second

template<typename T> inline T max(T a,T b){return a>b?a:b;}
template<typename T> inline T min(T a,T b){return a<b?a:b;}
template<typename T> inline T abs(T a){return a>0?a:-a;}
template<typename T> inline void repr(T &a,T b){if(a<b)a=b;}
template<typename T> inline void repl(T &a,T b){if(a>b)a=b;}
#define mp(a,b) std::make_pair(a,b)
#define pb push_back

const int N=200,M=3000;

struct edge
{
	int to,ne;
	double w;
}e[M*2+2];

int p[N+1],em=2,dep[N+1],q[N],qe,fa[N+1];
double fl[N+1];
bool vis[N+1];

inline void add(int a,int b,double w)
{
	e[em].to=b,e[em].w=w,e[em].ne=p[a],p[a]=em++;
	e[em].to=a,e[em].w=0,e[em].ne=p[b],p[b]=em++;
}

bool dfs(int s,int t)
{
	if(s==t)return 1;
	for(int j=p[s];j;j=e[j].ne)
		if(!vis[e[j].to]&&e[j].w&&dep[e[j].to]>dep[s])
		{
			vis[e[j].to]=1;
			fl[e[j].to]=min(e[j].w,fl[s]);
			fa[e[j].to]=j^1;
			if(dfs(e[j].to,t))return 1;
		}
	return 0;
}

inline double dinic(int s,int t)
{
	double flow=0;
	while(1)
	{
		memset(dep,0,sizeof(dep));
		q[0]=s,qe=1,dep[s]=1;
		for(int i=0;i^qe;i++)
		{
			for(int j=p[q[i]];j;j=e[j].ne)
				if(!dep[e[j].to]&&e[j].w)
				{
					dep[e[j].to]=dep[q[i]]+1;
					q[qe++]=e[j].to;
				}
		}
		if(!dep[t])break;
		while(1)
		{
			memset(vis,0,sizeof(vis));
			fl[s]=0x7fffffff;
			dfs(s,t);
			if(!vis[t])break;
			flow+=fl[t];
			for(int i=t;i!=s;i=e[fa[i]].to)
				e[fa[i]].w+=fl[t],e[fa[i]^1].w-=fl[t];
		}
	}
	return flow;
}

int n,m,a[50],b[50],atk[50][50];

inline bool judge(double time)
{
	int s=n+m,t=n+m+1;
	memset(p,0,sizeof(p));
	em=2;
	double s1=0,s2;
	for(int i=0;i<m;i++)
		add(s,i,time*b[i]);
	for(int i=0;i<n;i++)
		add(i+m,t,a[i]),s1+=a[i];
	for(int i=0;i<m;i++)
		for(int j=0;j<n;j++)
			if(atk[i][j])
				add(i,j+m,1e20);
	s2=dinic(s,t);
	if(abs(s1-s2)<1e-4)return 1;else return 0;
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int i=0;i<n;i++)
		scanf("%d",a+i);
	for(int i=0;i<m;i++)
		scanf("%d",b+i);
	for(int i=0;i<m;i++)
		for(int j=0;j<n;j++)
		{
			int k;
			scanf("%d",&k);
			atk[i][j]=k;
		}
	double l=0,r=1e4;
	while(r-l>1e-4)
	{
		if(judge((l+r)/2))
			r=(l+r)/2;
		else
			l=(l+r)/2;
	}
	printf("%.6lf\n",r);
}

1066:拆点+最大流 

#include<bits/stdc++.h>

typedef unsigned char uchar;
typedef unsigned int uint;
typedef long long ll;
typedef unsigned long long ull;

#define xx first
#define yy second

template<typename T> inline T max(T a,T b){return a>b?a:b;}
template<typename T> inline T min(T a,T b){return a<b?a:b;}
template<typename T> inline T abs(T a){return a>0?a:-a;}
template<typename T> inline void repr(T &a,T b){if(a<b)a=b;}
template<typename T> inline void repl(T &a,T b){if(a>b)a=b;}
#define mp(a,b) std::make_pair(a,b)
#define pb push_back

const int N=900,M=200000;

struct edge
{
	int to,w,ne;
}e[M*2+2];

int p[N+1],em=2,dep[N+1],q[N],qe,fa[N+1],fl[N+1];
bool vis[N+1];

inline void add(int a,int b,int w)
{
	e[em].to=b,e[em].w=w,e[em].ne=p[a],p[a]=em++;
	e[em].to=a,e[em].w=0,e[em].ne=p[b],p[b]=em++;
}

bool dfs(int s,int t)
{
	if(s==t)return 1;
	for(int j=p[s];j;j=e[j].ne)
		if(!vis[e[j].to]&&e[j].w&&dep[e[j].to]>dep[s])
		{
			vis[e[j].to]=1;
			fl[e[j].to]=min(e[j].w,fl[s]);
			fa[e[j].to]=j^1;
			if(dfs(e[j].to,t))return 1;
		}
	return 0;
}

inline int dinic(int s,int t)
{
	int flow=0;
	while(1)
	{
		memset(dep,0,sizeof(dep));
		q[0]=s,qe=1,dep[s]=1;
		for(int i=0;i^qe;i++)
		{
			for(int j=p[q[i]];j;j=e[j].ne)
				if(!dep[e[j].to]&&e[j].w)
				{
					dep[e[j].to]=dep[q[i]]+1;
					q[qe++]=e[j].to;
				}
		}
		if(!dep[t])break;
		while(1)
		{
			memset(vis,0,sizeof(vis));
			fl[s]=0x7fffffff;
			dfs(s,t);
			if(!vis[t])break;
			flow+=fl[t];
			for(int i=t;i!=s;i=e[fa[i]].to)
				e[fa[i]].w+=fl[t],e[fa[i]^1].w-=fl[t];
		}
	}
	return flow;
}

char ss[30];
bool iz[20][20];

int main()
{
	int n,m,dii,sum=0;
	scanf("%d%d%d",&n,&m,&dii);
	int s=n*m*2+5,t=n*m*2+10;
	for(int i=0;i<n;i++)
	{
		scanf("%s",ss);
		for(int j=0;j<m;j++)
		{
			iz[i][j]=ss[j]>'0';
			if(iz[i][j])add(i*m+j+n*m,i*m+j,ss[j]-'0');
			if((i<dii||j<dii||i>=n-dii||j>=m-dii)&&iz[i][j])
			{
				add(i*m+j,t,11111);
			}
		}
	}
	for(int i=0;i<n;i++)
	{
		scanf("%s",ss);
		for(int j=0;j<m;j++)
			if(ss[j]=='L')
			{
				add(s,i*m+j+n*m,1);
				sum++;
			}
	}
	for(int i=0;i<n;i++)
		for(int j=0;j<m;j++)
			for(int k=0;k<n;k++)
				for(int l=0;l<m;l++)
					if((i-k)*(i-k)+(j-l)*(j-l)<=dii*dii&&(i!=k||j!=l)&&iz[i][j]&&iz[k][l])
					{
						add(i*m+j,k*m+l+n*m,11111);
					}
	printf("%d\n",sum-dinic(s,t));
}

1433:二分图最大匹配，直接网络流跑 

#include<bits/stdc++.h>

typedef unsigned char uchar;
typedef unsigned int uint;
typedef long long ll;
typedef unsigned long long ull;

#define xx first
#define yy second

template<typename T> inline T max(T a,T b){return a>b?a:b;}
template<typename T> inline T min(T a,T b){return a<b?a:b;}
template<typename T> inline T abs(T a){return a>0?a:-a;}
template<typename T> inline void repr(T &a,T b){if(a<b)a=b;}
template<typename T> inline void repl(T &a,T b){if(a>b)a=b;}
#define mp(a,b) std::make_pair(a,b)
#define pb push_back

const int N=120,M=3200;

struct edge
{
	int to,w,ne;
}e[M*2+2];

int p[N+1],n,em=2,dep[N+1],q[N],qe,fa[N+1],fl[N+1];
bool vis[N+1];

inline void add(int a,int b,int w)
{
	e[em].to=b,e[em].w=w,e[em].ne=p[a],p[a]=em++;
	e[em].to=a,e[em].w=0,e[em].ne=p[b],p[b]=em++;
}

bool dfs(int s,int t)
{
	if(s==t)return 1;
	for(int j=p[s];j;j=e[j].ne)
		if(!vis[e[j].to]&&e[j].w&&dep[e[j].to]>dep[s])
		{
			vis[e[j].to]=1;
			fl[e[j].to]=min(e[j].w,fl[s]);
			fa[e[j].to]=j^1;
			if(dfs(e[j].to,t))return 1;
		}
	return 0;
}

inline int dinic(int s,int t)
{
	int flow=0;
	while(1)
	{
		memset(dep,0,sizeof(dep));
		q[0]=s,qe=1,dep[s]=1;
		for(int i=0;i^qe;i++)
		{
			for(int j=p[q[i]];j;j=e[j].ne)
				if(!dep[e[j].to]&&e[j].w)
				{
					dep[e[j].to]=dep[q[i]]+1;
					q[qe++]=e[j].to;
				}
		}
		if(!dep[t])break;
		while(1)
		{
			memset(vis,0,sizeof(vis));
			fl[s]=0x7fffffff;
			dfs(s,t);
			if(!vis[t])break;
			flow+=fl[t];
			for(int i=t;i!=s;i=e[fa[i]].to)
				e[fa[i]].w+=fl[t],e[fa[i]^1].w-=fl[t];
		}
	}
	return flow;
}

int isc[60];

inline void solve()
{
	em=2;
	memset(p,0,sizeof(p));
	int rn;
	scanf("%d",&rn);
	int s=rn*2+5,t=rn*2+6;
	for(int i=1;i<=rn;i++)
	{
		scanf("%d",isc+i);
		if(isc[i])add(i,t,1);
	}
	int rq=0;
	for(int i=1;i<=rn;i++)
	{
		int k;
		scanf("%d",&k);
		if(k==0||!isc[i])add(s,i+rn,1),rq++;
	}
	for(int i=1;i<=rn;i++)
		for(int j=1;j<=rn;j++)
		{
			int k;
			scanf("%d",&k);
			if(k==1||i==j)add(i+rn,j,1);
		}
	if(dinic(s,t)==rq)puts("^_^");else puts("T_T");
}

int main()
{
	int T;
	scanf("%d",&T);
	while(T--)solve();
}

link to my blog

CDQZOIERS
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
BZOJ 4514 浅谈最大费用最大流
BerryKanry的博客
08-25 2857
世界真的很大 网络流并不难，但是其建图方式的多种多样，决定了其能解决的问题的多样性。 而作为网络流的一个变种，在其基础之上加上边权的费用流，除了一眼就能看出的模板题以外，取决于建模方式的灵活变化也很多 而对于这样变化多端的题而言，我们往往很难一眼看出来这道题究竟该不该用网络流或者费用流 我们寻求题解，但却不能仅仅止步于AC，不应该忘记去推导建模的过程 在感叹“哇这道题真是巧”的同时，还是得
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
网络流dinic最大流板子 bzoj[1001]
Richard__Luan的博客
07-12 358
传送门 bzoj【1001】最小割最大流 据说此题是对偶图最短路~~~ 最大流过了，就当是练网络流了吧 代码如下 #include&lt;cstdio&gt; #include&lt;iostream&gt; #include&lt;queue&gt; #include&lt;cstring&gt; using namespace std; int n,m; int fst[10000...
[bzoj2163]复杂的大门 最大流
Unlimited Code Works
08-31 390
2163: 复杂的大门 Time Limit: 20 Sec  Memory Limit: 259 MB [Submit][Status][Discuss] Description 你去找某bm玩，到了门口才发现要打开他家的大门不是一件容易的事…… 他家的大门外有n个站台，用1到n的正整数编号。你需要对每个站台访问一定次数以后大门才能开启。站台之间有m个单向的传送门，通过传送门到达另
bzoj2929 [Poi1999]洞穴攀行
weixin_30509393的博客
12-09 81
Description 一队洞穴学者在Byte Mountain的Grate Cave里组织了一次训练。训练中，每一位洞穴学者要从最高的一个室到达最底下的一个室。他们只能向下走。一条路上每一个连续的室都要比它的前一个低。此外，每一个洞穴学者都要从最高的室出发，沿不同的路走到最低的室。问：可以有多少个人同时参加训练？ 任务： 写一个程序： l读入对洞穴的描述。 l...
[BZOJ2163]复杂的大门（最大流）
zyf2000
02-21 1183
题目描述传送门题解『转化』 由于传送门的使用次数是有限的，没有代价，而公交车是任意两点直达，代价均为1，那么问题就可以转化为最多可以使用多少次传送门。 『拆点』 我们知道到达一个点即意味着也要从这个点离开，由于传送门单向，我们自然可以想到拆点。 『建图』 按照传送门的起止点连边，容量为传送门最大使用次数；在从源点分别向每个起点连边，从每个终点分别向汇点连边，容量分别为每一个点应到达的次数。
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
bzoj problem
08-21
【标题】"bzoj problem" 指的是“八中OJ（Online Judge）问题集”，这是一个在线编程竞赛平台的题目合集。八中OJ可能是某所中学的在线编程训练系统，它提供了大量的编程题目供用户练习和挑战，以提升编程能力和算法...
BZOJ题目镜像
09-14
【BZOJ题目镜像】是为了解决BZOJ（Blue Zoo Online Judge）系统因服务器问题时常出现访问不稳定或无法访问的情况而设立的一个备份资源。BZOJ是中国早期的在线编程竞赛平台，提供了大量的算法题目供编程爱好者进行...
BZOJ第二部分
11-04
【BZOJ第二部分】是针对BZOJ（北京邮电大学在线评测系统）的一个专题，这个专题主要涵盖了BZOJ平台上的P2001到P3000之间的编程题目。对于想要深入学习算法、提升编程能力的IT从业者或学生来说，这是一个宝贵的资源...
【BZOJ】 [NOI2006]最大获利-网络流最大流最小割
远行客
04-29 403
传送门:BZOJ1497-最大获利 题意 在前期市场调查和站址勘测之后，公司得到了一共N个可以作为通讯信号中转站的地址，而由于这些地址的地理位置差异，在不同的地方建造通讯中转站需要投入的成本也是不一样的，所幸在前期调查之后这些都是已知数据：建立第i个通讯中转站需要的成本为Pi（1≤i≤N）。另外公司调查得出了所有期望中的用户群，一共M个。关于第i个用户群的信息概括为Ai, Bi和Ci：...
BZOJ 3994 （思路最大流）
playwfun的专栏
09-18 281
因为用时可以二分，所以先二分最小值，那么判断用最大流就可以了。 好像不用解释什么，对于源点对每个装备连一条边，边权为时间*装备杀伤力，经过可行性矩阵 ，流向n个机器人，每个机器人向汇点流自己的防御值。 那么如果最大流等于所有机器人防御值之和，那么最大流方案，既是该时间下的可行性方案。 //#pragma comment(linker, "/STACK:1024000000,102400000
bzoj 1066 蜥蜴(最大流)
glk__'blog
07-18 579
题目传送门 Description 在一个r行c列的网格地图中有一些高度不同的石柱，一些石柱上站着一些蜥蜴，你的任务是让尽量多的蜥蜴逃 到边界外。 每行每列中相邻石柱的距离为1，蜥蜴的跳跃距离是d，即蜥蜴可以跳到平面距离不超过d的任何一个石 柱上。石柱都不稳定，每次当蜥蜴跳跃时，所离开的石柱高度减1（如果仍然落在地图内部，则到达的石柱高度不 变），如果该石柱原来高度为1，则蜥蜴离
BZOJ1458 士兵占领-最大流
Mys_C_K的博客
05-17 379
传送门 题目大意：自行参考 题解：首先一看方格图最值问题就是网络流（尽管我是知道了网络流才做的QwQ） 肯定是对每一行建一排点，每一列建一排点，建成一个二分图。 显然行和列连边容量为1. 首先注意到要求最小值，发现最小割没办法处理（因为从s对每一行的连边割去是没有意义的，因为可能只“割去”一部分） 然后注意到，如果s和每一行连边的容量表示“最多/最少能放多少个士兵的时候”，这就是个有上
bzoj 1711 最大流
Eirlys的博客
03-08 483
题意：n头牛，m1种吃的，m3种喝的，每头牛有nm1种喜欢的吃的、nm2种喜欢的喝的，满足一头牛的条件是 ：既给他一个爱吃的也给他一个爱喝的。（注意如果一头牛没有爱吃的或爱喝的，是永远不能满足的）。每种吃的和喝的能且只能用一次，求最多满足多少牛 最大流还是比较显然的，但是建图还是相当经（恶）典（心）的 给出正确的建图： s -> 每种吃的 容量为1 （保证每种吃的能且只能用一次） 每种喝的
bzoj1458-wo4263 士兵占领【最大流】
修炼之路
12-26 173
题意 有一个M * N的棋盘，有的格子是障碍。现在你要选择一些格子来放置一些士兵，一个格子里最多可以放置一个士兵，障碍格里不能放置士兵。我们称这些士兵占领了整个棋盘当满足第i行至少放置了Li个士兵, 第j列至少放置了Cj个士兵。现在你的任务是要求使用最少个数的士兵来占领整个棋盘。 建图 每一行建立Ai，与源点S连接，容量为ri 每一列建立Bi，与汇点连，容量为ci 第i行第j列可以放士兵，Ai-&amp;amp;...
BZOJ 1693 [Usaco2007 Demo]Asteroids 最大流
-Claude
08-19 1315
BZOJ 1693 [Usaco2007 Demo]Asteroids 最大流
[bzoj3504][网络流-最大流]危桥
Rose_max的博客
03-04 293
Description Alice和Bob居住在一个由N座岛屿组成的国家，岛屿被编号为0到N-1。某些岛屿之间有桥相连，桥上的道路是双 向的，但一次只能供一人通行。其中一些桥由于年久失修成为危桥，最多只能通行两次。Alice希望在岛屿al和a2之间往返an次（从al到a2再从a2到al算一次往返）。同时，Bob希望在岛屿bl和b2之间往返bn次。这个过程中，所有危桥最多通行两次，其余...
BZOJ 3903 反垄断 (最大流推的结论题)
Ark
03-12 220
题目 中文题目,不解释: BZOJ传送门 分析 这道题BZOJ上也只有几个人过…奇怪了 下面是正解 原问题为一个二分图边染色问题。首先考虑最好情况。最理想情况的分配为：设一个点xxx的度为dgr(x)dgr(x)dgr(x)。那么若dgr(x)dgr(x)dgr(x)为ttt的倍数，那么与xxx相连的边恰好平均分配给ttt个公司，每个公司dgr(x)/tdgr(x)/tdgr(x)/t个；如果不...
bzoj.org俄罗斯方块
最新发布
05-06
bzoj.org是一个在线的编程题库，其中包含了很多经典算法问题和实际应用中的编程问题。其中包括了很多有趣的游戏算法问题，比如Puzzle，Tetris等。 而bzoj.org俄罗斯方块就是在该平台上的一个经典的俄罗斯方块游戏。在该游戏中，玩家需要控制方块不断向下移动，直到方块落到底部或者与其他方块重叠。游戏过程中，玩家可以通过旋转和移动方块来调整方块的位置，使其能够合理地落下。当一行方块被填满时，该行会被清除，并获得相应的分数。 bzoj.org俄罗斯方块不仅仅是一款游戏，更是一个很好的算法练习平台，你可以通过该游戏来学习和研究游戏算法问题。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值