HDU6532（2019湘潭邀请赛A题）Chessboard（费用流）

最新推荐文章于 2019-08-14 01:02:56 发布
最新推荐文章于 2019-08-14 01:02:56 发布
阅读量676 收藏 1
点赞数
分类专栏： 网络流 图论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/baiyifeifei/article/details/90425878
版权

HDU6532（2019湘潭邀请赛A题）Chessboard（费用流）

题目大意

给出一张棋盘，其上一些位置上有一些棋子。现在又给出一些限制条件，如某行以下最多可以选取多少的棋子，某列以右最多可以选取多少棋子。每个棋子的位置按照一定的顺序给出，其顺序中的排名即是其价值，问可以达到的最大价值是多少

解题思路

非常常规的思路，将行和列拆开，如果在某个位置有一个点则将其行和列连接上一条边，其容量为1，边权为其标号。但是行和列很多，所以还需要对行和列进行离散化。每个离散化出的行与其上一行，每个列与其左边的一个列之间都有可以有若干条限制，选取其中最小的限制即是限制离散化出来的点的流量的限制根据每个点的限制，将其与上一点点之间连接上容量为限制，价值为0的边，最后跑费用流即可

AC代码

#include<bits/stdc++.h>
using namespace std;
const int sizep=505;
const int inf=0x3f3f3f3f;
const int maxn=2005;
const int maxm=2e5+5;
struct Edge{
	int to,next,flow,cost;
	int cap;
} edge[maxm];

int head[maxn],tol;
int pre[maxn];
int dis[maxn];
bool vis[maxn];
int N;

void init(int n)
{
	N=n;tol=0;
	for(int i=0;i<n;i++)
		head[i]=-1;
}

void AddEdge(int u,int v,int cap,int cost)
{
	edge[tol].to=v;
	edge[tol].cap=cap;
	edge[tol].cost=cost;
	edge[tol].flow=0;
	edge[tol].next=head[u];
	head[u]=tol++;
	
	edge[tol].to=u;
	edge[tol].cap=0;
	edge[tol].flow=0;
	edge[tol].cost=-cost;
	edge[tol].next=head[v];
	head[v]=tol++;
}

bool spfa(int s,int t)
{
	queue<int >q;
	for(int i=0;i<=N;++i)
	{
		dis[i]=inf;
		vis[i]=false;
		pre[i]=-1;
	}
	dis[s]=0;
	vis[s]=true;
	q.push(s);
	while(!q.empty())
	{
		int u=q.front(); q.pop();
		vis[u]=false;
		for(int i=head[u];~i;i=edge[i].next)
		{
			int v=edge[i].to;
			if(edge[i].cap>edge[i].flow && dis[v]>dis[u]+edge[i].cost)
			{
				dis[v]=dis[u]+edge[i].cost;
				pre[v]=i;
				if(!vis[v])
				{
					vis[v]=true;
					q.push(v);
				}
			}
		}
	}
	if(pre[t]==-1) return false;
	else return true;
}

int MincostMaxflow(int s,int t,int &cost)
{
	int flow=0;
	cost=0;
	while(spfa(s,t))
	{
		int Min=inf;
		for(int i=pre[t];i!=-1;i=pre[edge[i^1].to])
		{
			if(Min>edge[i].cap -edge[i].flow)
				Min=edge[i].cap -edge[i].flow;
		}
		for(int i=pre[t];i!=-1;i=pre[edge[i^1].to])
		{
			edge[i].flow+=Min;
			edge[i^1].flow-=Min;
			cost+=edge[i].cost*Min;
		}
		flow+=Min;
	}
	return flow;
}
struct point{
	int x,y;
	point(){}
	point(int x,int y):x(x),y(y){}
};
struct limit{
	int loc;
	int cap;
	limit(){}
	limit(int loc,int cap):loc(loc),cap(cap){}
	bool friend operator<(limit a,limit b)
	{
		return a.loc<b.loc;
	}
};
point p[sizep];
int capx[sizep],capy[sizep];
vector<limit> limx,limy;
vector<int> vx;
vector<int> vy;
int idx[sizep];
int idy[sizep];
int main()
{
	int n;
	int x,y;
	while(~scanf("%d",&n))
	{
		vx.clear();
		vy.clear();
		limx.clear();
		limy.clear();
		for(int i=1;i<=n;i++)
		{
			scanf("%d%d",&x,&y);
			p[i]=point(x,y);
			vx.push_back(x);
			vy.push_back(y);
		}
		sort(vx.begin(),vx.end());
		sort(vy.begin(),vy.end());
		vx.erase(unique(vx.begin(),vx.end()),vx.end());
		vy.erase(unique(vy.begin(),vy.end()),vy.end());
		for(int i=1;i<=n;i++)
		{
			idx[i]=lower_bound(vx.begin(),vx.end(),p[i].x)-vx.begin()+1;
			idy[i]=lower_bound(vy.begin(),vy.end(),p[i].y)-vy.begin()+1;
		}
		int totx,toty;
		totx=vx.size(),toty=vy.size();
		int s=0,t=totx+toty+1;
		int m;
		scanf("%d",&m);
		char op[5];
		int loc,k;
		for(int i=1;i<=m;i++)
		{
			scanf("%s%d%d",op,&loc,&k);
			if(op[0]=='R') 
			{
				limx.push_back(limit(loc,k));
			}
			if(op[0]=='C')
			{
				limy.push_back(limit(loc,k));
			}
		}
		sort(limx.begin(),limx.end());
		sort(limy.begin(),limy.end());
		int cntx=0,cnty=0;
		memset(capx,inf,sizeof(capx));
		memset(capy,inf,sizeof(capy));
		for(int i=0;i<limx.size();i++)
		{
			while(cntx<totx-1&&vx[cntx]<limx[i].loc)
			{
				cntx++;
			}
			capx[cntx+1]=min(capx[cntx+1],limx[i].cap);
		}
		for(int i=0;i<limy.size();i++)
		{
			while(cnty<toty-1&&vy[cnty]<limy[i].loc)
			{
				cnty++;
			}
			capy[cnty+1]=min(capy[cnty+1],limy[i].cap);
		}
		init(totx+toty+2);
		for(int i=1;i<=n;i++)
		{
			AddEdge(idx[i],idy[i]+totx,1,-i);
		}
		for(int i=1;i<=totx;i++)
		{
			AddEdge(i-1,i,capx[i],0);
		}
		for(int i=1;i<=toty;i++)
		{
			if(i==1)
			AddEdge(totx+i,t,capy[i],0);
			else AddEdge(totx+i,totx+i-1,capy[i],0);
		}
		int cost=0;
		MincostMaxflow(s,t,cost);
		printf("%d\n",-cost);
	}
}

5.23更新

感谢山东理工大学的一位同学，原代码可以通过

5
1  2
2  2
3  1
100 100
3  2
5
R 2 3
C 1 4
R 3 1
C 2 2
C 5 0

叉掉。理由是，简单地使用费用流得到的答案会优先保证最大的流量，而有时可能会选择更少但是权值更大的点，且其为正解，解决方案是在网络流更新增广路时，当距离为正就停止更新。这样就可以得到最大费用。正解已更新

#include<bits/stdc++.h>
using namespace std;
const int sizep=505;
const int inf=0x3f3f3f3f;
const int maxn=2005;
const int maxm=2e5+5;
struct Edge{
	int to,next,flow,cost;
	int cap;
} edge[maxm];
int head[maxn],tol;
int pre[maxn];
int dis[maxn];
bool vis[maxn];
int N;
void init(int n)
{
	N=n;tol=0;
	for(int i=0;i<n;i++)
		head[i]=-1;
}

void AddEdge(int u,int v,int cap,int cost)
{
	edge[tol].to=v;
	edge[tol].cap=cap;
	edge[tol].cost=cost;
	edge[tol].flow=0;
	edge[tol].next=head[u];
	head[u]=tol++;
	
	edge[tol].to=u;
	edge[tol].cap=0;
	edge[tol].flow=0;
	edge[tol].cost=-cost;
	edge[tol].next=head[v];
	head[v]=tol++;
}

bool spfa(int s,int t)
{
	queue<int >q;
	for(int i=0;i<=N;++i)
	{
		dis[i]=inf;
		vis[i]=false;
		pre[i]=-1;
	}
	dis[s]=0;
	vis[s]=true;
	q.push(s);
	while(!q.empty())
	{
		int u=q.front(); q.pop();
		vis[u]=false;
		for(int i=head[u];~i;i=edge[i].next)
		{
			int v=edge[i].to;
			if(edge[i].cap>edge[i].flow && dis[v]>dis[u]+edge[i].cost)
			{
				dis[v]=dis[u]+edge[i].cost;
				pre[v]=i;
				if(!vis[v])
				{
					vis[v]=true;
					q.push(v);
				}
			}
		}
	}
	if(pre[t]==-1) return false;
	else return true;
}

int MincostMaxflow(int s,int t,int &cost)
{
	int flow=0;
	cost=0;
	while(spfa(s,t))
	{
		if(dis[t]>0) break;
		int Min=inf;
		for(int i=pre[t];i!=-1;i=pre[edge[i^1].to])
		{
			if(Min>edge[i].cap -edge[i].flow)
				Min=edge[i].cap -edge[i].flow;
		}
		for(int i=pre[t];i!=-1;i=pre[edge[i^1].to])
		{
			edge[i].flow+=Min;
			edge[i^1].flow-=Min;
			cost+=edge[i].cost*Min;
		}
		flow+=Min;
	}
	return flow;
}
struct point{
	int x,y;
	point(){}
	point(int x,int y):x(x),y(y){}
};
struct limit{
	int loc;
	int cap;
	limit(){}
	limit(int loc,int cap):loc(loc),cap(cap){}
	bool friend operator<(limit a,limit b)
	{
		return a.loc<b.loc;
	}
};
point p[sizep];
int capx[sizep],capy[sizep];
vector<limit> limx,limy;
vector<int> vx;
vector<int> vy;
int idx[sizep];
int idy[sizep];
int main()
{
	int n;
	int x,y;
	while(~scanf("%d",&n))
	{
		vx.clear();
		vy.clear();
		limx.clear();
		limy.clear();
		for(int i=1;i<=n;i++)
		{
			scanf("%d%d",&x,&y);
			p[i]=point(x,y);
			vx.push_back(x);
			vy.push_back(y);
		}
		sort(vx.begin(),vx.end());
		sort(vy.begin(),vy.end());
		vx.erase(unique(vx.begin(),vx.end()),vx.end());
		vy.erase(unique(vy.begin(),vy.end()),vy.end());
		for(int i=1;i<=n;i++)
		{
			idx[i]=lower_bound(vx.begin(),vx.end(),p[i].x)-vx.begin()+1;
			idy[i]=lower_bound(vy.begin(),vy.end(),p[i].y)-vy.begin()+1;
		}
		int totx,toty;
		totx=vx.size(),toty=vy.size();
		int s=0,t=totx+toty+1;
		int m;
		scanf("%d",&m);
		char op[5];
		int loc,k;
		for(int i=1;i<=m;i++)
		{
			scanf("%s%d%d",op,&loc,&k);
			if(op[0]=='R') 
			{
				limx.push_back(limit(loc,k));
			}
			if(op[0]=='C')
			{
				limy.push_back(limit(loc,k));
			}
		}
		sort(limx.begin(),limx.end());
		sort(limy.begin(),limy.end());
		int cntx=0,cnty=0;
		memset(capx,inf,sizeof(capx));
		memset(capy,inf,sizeof(capy));
		for(int i=0;i<limx.size();i++)
		{
			while(cntx<totx-1&&vx[cntx]<limx[i].loc)
			{
				cntx++;
			}
			capx[cntx+1]=min(capx[cntx+1],limx[i].cap);
		}
		for(int i=0;i<limy.size();i++)
		{
			while(cnty<toty-1&&vy[cnty]<limy[i].loc)
			{
				cnty++;
			}
			capy[cnty+1]=min(capy[cnty+1],limy[i].cap);
		}
		init(totx+toty+2);
		for(int i=1;i<=n;i++)
		{
			AddEdge(idx[i],idy[i]+totx,1,-i);
		}
		for(int i=1;i<=totx;i++)
		{
			AddEdge(i-1,i,capx[i],0);
		}
		for(int i=1;i<=toty;i++)
		{
			if(i==1)
			AddEdge(totx+i,t,capy[i],0);
			else AddEdge(totx+i,totx+i-1,capy[i],0);
		}
		int cost=0;
		MincostMaxflow(s,t,cost);
		printf("%d\n",-cost);
	}
}
Flyppy_White
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
HDU-6532 Chessboard 2019广东省省赛B题（费用流）
常胜将军的博客
06-19 385
题意：给出一张棋盘，其上一些位置上有一些棋子。现在又给出一些限制条件，如某行以下最多可以选取多少的棋子，某列以右最多可以选取多少棋子。每个棋子的位置按照一定的顺序给出，其顺序中的排名即是其价值，问可以达到的最大价值是多少。 题解：费用流。离散化坐标，每行用一个点表示，每列也用一个点表示。表示第i-1行的点向表示第i行的点连边，容量为第i行及以后能拿的棋子数的上限，费用为0，同理表示相邻列的点两两...
HDU6532 Chessboard (最大费用流)
weixin_30301449的博客
07-11 144
题意:棋盘上有n个棋子 每个棋子都有收益 　　 现在给定1e5条线 有横着的 竖着的 　　 规定只能在线的一侧选最多ki个棋子 　　 问最大收益 题解:写自闭的一道题 很容易想到是网络流 但是建图有点难 　　 第一道最大费用流 居然是边权取反 跑最小费用最大流! 　　 先离散化坐标 然后可以用点代替一条横线 一条竖线 　　 如果x，y有一个棋子 就将x所在的横线 向y所在的...
2019牛客暑期多校训练（第七场）I-Chessboard
BenFromHRBUST的CSDN博客
08-12 402
2019牛客暑期多校训练（第七场）J-Just Jump 题目描述 CSL likes playing tricks very much. The unlucky TL always becomes the target of CSL’s tricks. Today CSL came to trick TL again. He saw that TL has a very very large c...
2019牛客暑期多校训练营（第七场）I ChessBoard （组合数学）
zx93的博客
08-14 209
题目链接：https://ac.nowcoder.com/acm/contest/887/I 题意：给定n,m。求k*k的矩阵，1<=k<=n，矩阵内的每个元素都不小于m，且矩阵内不同行不同列的元素相加都为一个定值T，且T<=n。问这样的矩阵有多少种，MOD998244353 官方题解： 容易证明，一个k*k棋盘的排布方案满足“不同行不同列的方格内的玻璃球数量的总和...
Chessboard HDU - 6532 (最大费用流）
cat_pikapikapi的博客
05-23 299
题意：给你一个棋盘，棋盘上有一些棋子，有一些规定，规定超过多少行或者多少列后只能取不超过k的棋子。 题解：棋盘问题是一个很经典的网络流问题。 建图还是比较好想的，但是不会跑最大费用流，赛后补了一下。 但是好像这道题有问题，最大费用最大流也能过。 但是明显不符合费用最大时流最大。 所以想一下最小费用最大流的做法，首先费用取反，这样就变成了最小费用， 每次跑一条s到t的最短路，然后减去这条...
HDU 递归题详解大全（含代码）
02-11
蟠桃记 1 折线分割平面 2 不容易系列之一 2 骨牌铺方格 3 不容易系列之(3)—— LELE的RPG难题 3 Children’s Queue 3 献给杭电五十周年校庆的礼物 3 钥匙计数之二 3 钥匙计数之一 3 母牛的故事 3 ...
hdu.rar_HDU 1089.cpp_OJ题求和_hdu_horsekw5_杭电obj
最新发布
09-24
【标题】"hdu.rar_HDU 1089.cpp_OJ题求和_hdu_horsekw5_杭电obj" 提供的信息是关于一个压缩文件，其中包含了一个名为 "HDU 1089.cpp" 的源代码文件，这个文件是为了解决杭州电子科技大学（Hangzhou Dianzi ...
杭电HDU2050题的ac程序
04-10
一个十分简单的程序，能够ac杭电hdu的第2050题，无注释，简单明了
hdu-page-11-answer.rar_hdu_hdu oj第十一页_page_搜题_杭电oj
09-14
"HDU Page 11 Answer"提供的题解资源，不仅节省了逐一搜索题解的时间，还提供了一次性获取多题答案的机会。通过对比不同解法，学习者可以开阔思路，理解多种解决问题的角度，进一步提升编程素养。 总之，"HDU Page ...
hdu 6532 费用流
人生如戏
05-20 428
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6532 题意： 在1e9*1e9的棋盘上给了你n个棋子(n<=500)，同时给了你m条限制，每条限制是一个 s a b的结构，s为R或者C的字符，R表示是在第a和a-1行之间画一条线，在这条线的下面最多只能取b个点，C同理，是在线的右边只能取这么多点，每个点i带有的权值是i。...
HDU-4035 Maze （概率DP&&树形DP）
idealism_xxm的专栏
04-30 932
理解题以后觉得没法做。。。只理解了有向图，因为本题是无向边，所以父子结点的概率互相影响，立刻就忘了上一题也存在这种情况。。。 看了题解后才明白，应该按照期望设状态，而且需要在树上进行状态转移，最开始的状态转移方程最后化简即可消除影响（貌似就是树形DP+概率DP） 设dp[i]表示从i号结点走出迷宫时走过边数的概率，par[i]表示i号结点的父亲结点，son表示i号结点的儿子结点： ①若i结点为叶子结点，则 dp[i]=k[i]*dp[1]+e[i]*0+(1-k[i]-e[i])*(dp[par[i]]+1
Chessboard【暴力(brute force)】
Vanishment's Blog
04-05 176
问题来源：codeforces问题链接：http://codeforces.com/problemset/problem/961/C问题描述：Magnus decided to play a classic chess game. Though what he saw in his locker shocked him! His favourite chessboard got broken in...
CCPC之湘潭邀请赛
weixin_44023434的博客
05-15 649
CCPC湘潭邀请赛 热身赛 5月11号，当天下午是打热身赛，四道题目，我本来以为是英文题，还特意让钢哥回去拿了下字典，可没想到热身赛是中文题，幸苦钢哥白跑一趟了。 再来说题目吧，几个题目意思都比较明了，知道怎么写，但是就是写不下去，例如A题，就是简单暴力模拟题，思路很简单，但是我可能写炸了，心态有点爆炸，不知道怎么模拟下去，然后一直卡着时间，队友们然后去开后面几题。接着巧巨说她知道怎么写D题了，然...
HDU 1256 画8
程序猿之洞
02-08 1173
画8 Time Limit : 2000/1000ms (Java/Other)   Memory Limit : 65536/32768K (Java/Other) Total Submission(s) : 2   Accepted Submission(s) : 2 Problem Description 谁画8画的好,画的快,今后就发的快,学业发达,事业发达,祝大家
2019ccpc湖南全国邀请赛后记
weixin_30246221的博客
05-13 140
赛前 1.因为去年和湘潭大学的，我们教练会写签到题队撞了好多个赛区，所以这次比赛的队名被我灵机一动改为了听说你们教练会写签到题。 2.由于去年就在这个比赛上拿过金,今年的压力反而更大了(拿不到金就是菜) 3.和去年一样逛了步步高广场,让人迷路的地下商业街,感觉湘潭没有什么变化,尤其是菜还是一如既往的辣. DAY1 热身赛看了A题发现是无脑模拟,敲完还蹭了一个一血就很快乐,D发现是无脑最短...
关于2019湘潭全国邀请赛
nuoyanli的博客
05-13 2319
行程安排： 行程体验： 坐坐车，旅旅游，面面基，体验极佳（主要是不要花自己的钱）。 比赛体验： 四道签到题（L E K B），大概div3手速场，四题可银可铁。（至于其他难度....） 比赛开始前我指定的策略就是以hyh为中心快速读题，切签到题，然后冲冲冲，然而我上去就读到了签到题L E，就叫hyh看L写L，然后大胆交wa了，然来是题目理解错了，帮队友修了一下代码后就 AC 了，...
2019CCPC湖南全国邀请赛（广东省赛、江苏省赛）重现赛
Up Panda
05-19 588
我就写了签到题~~ 题面是xtCpc #include<bits/stdc++.h> using namespace std; int main() { int n; char k; ios::sync_with_stdio(false); while(cin>>n) { unordered_map<char,int>mp; mp['x']=1...
2019CCPC湘潭全国邀请赛题解
菜鸡成长史
05-21 2023
目录 【正赛情况】​ 【A. Chessboard】 【B. Build Tree】 【C. Chika and Friendly Pairs】 【D. Chika and Solid Geometry】 【E.Hello XTCPC】 【F. Neko and function】 【G. Neko and quadrilateral】 【H. Nek...
"2021大学生算法设计超级联赛题目集1及HDU多校联合训练比赛6题目A解析
同时也举办了2021 HDU Multi-University Training Contest 6 (XJTU Contest)Online, August 5, 2021比赛。比赛的首个问题是关于政府部门中的人员管理，要求参赛者设计一个算法来解决部门中的人员过剩问题。 问题的...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值