ZQC的游戏 题解

该文介绍了如何运用贪心策略和网络流算法解决一个问题,即在给定条件下找到具有最大重量的球。首先按照贪心思想,优先处理1号球,然后构建网络,设置源点和汇点,确保网络的最大流等于所有球的重量之和。通过ISAP(增广路径)算法求解最大流,从而得出答案。若最大流不等于球的重量之和,则表示无法满足条件。
摘要由CSDN通过智能技术生成

前言

这题题意描述不是很清楚啊,所以我找了个有权限的人把题面改了改,应该还是比较清楚了。

感觉这道题挺妙的,就来写一篇题解。

思路

首先,根据贪心思想,我们会将 1 1 1 号点半径以内能吃的都吃了,假设吃完之后它的重量为 s u m sum sum

那么,为了让它成为最大的,第 i i i 个人吃的重量必须满足 e a t i ≤ s u m − w i eat_i \le sum-w_i eatisumwi

那么既然如此,我们就可以考虑建一个网络,对于第 i i i 个人,将他与他半径以内能吃的东西连一条容量为极大值的边,随便找一个入点 s s s 和汇点 t t t。显然,对于第 i i i 个人, s s s 就向他连一条容量为 s u m − w i sum-w_i sumwi 的边(显然,如果更大,那么第一个点就不是最大的了),并对于任意一个球,向汇点连一条边,容量为 w i w_i wi(每个吃的只能被吃一次)。由于每个能吃的球,必须被吃完,所以必须要保证这个网络的最大流是满流(即最大流等于网络中每个球的重量之和),否则,球就吃不完,不满足题意,故第一个节点也就成不了最大的重量。

注意事项:有些球对于所有的人他都吃不到,要自动忽略这些点。

代码

#include<bits/stdc++.h>
using namespace std;
#define maxn 805
#define maxe 40005
int n,m,s,t;
int nx[maxn],ny[maxn],nw[maxn],r[maxn];
int x[maxn],y[maxn],w[maxn];
bool vis[maxn];
struct node
{
	int tar,nxt;
	long long num,nu2;
}arr[maxe<<1];
int fst[maxn],cnt=1;
void adds(int x,int y,long long z)
{
	arr[++cnt].tar=y,arr[cnt].nxt=fst[x],fst[x]=cnt,arr[cnt].num=z;
}
double dis(int x,int y,int z,int w)//求两点距离
{
	return sqrt((x-z)*(x-z)+(y-w)*(y-w));
}
namespace ISAP//板子
{
	int dis[maxn],now[maxn],gap[maxn];
	long long flow;
	void get_augmentpath()
	{
		memset(dis,0x3f,sizeof(dis));
		memset(now,0,sizeof(now));
		memset(gap,0,sizeof(gap));
		queue<int> p;
		p.push(t);
		dis[t]=0;
		now[t]=fst[t];
		gap[0]=1;
		while(!p.empty())
		{
			int x=p.front();
			p.pop();
			for(int i=fst[x];i;i=arr[i].nxt)
			{
				int j=arr[i].tar;
				if(dis[j]==0x3f3f3f3f)
				{
					p.push(j);
					now[j]=fst[j];
					dis[j]=dis[x]+1;
					gap[dis[j]]++;
				}
			}
		}
		return;
	}
	long long dfs(int x,long long sum)
	{
		if(x==t)
		{
			flow+=sum;
			return sum;
		}
		long long l=0;
		for(int i=now[x];i;i=arr[i].nxt)
		{
			now[x]=i;
			int j=arr[i].tar;
			long long k=arr[i].num;
			if(k>0&&dis[j]==dis[x]-1)
			{
				long long used=dfs(j,min(k,sum-l));
				if(used)
				{
					arr[i].num-=used;
					arr[i^1].num+=used;
					l+=used;
				}
				if(l==sum) return l;
			}
		}
		--gap[dis[x]];
		if(!gap[dis[x]]) dis[s]=n+1;
		++dis[x];
		gap[dis[x]]++;
		return l;
	}
	int output()
	{
		flow=0;
		get_augmentpath();
		while(dis[s]<n) memcpy(now,fst,sizeof(now)),dfs(s,LONG_LONG_MAX);
		return flow;
	}
}
void input()
{
	memset(vis,0,sizeof(vis));
	cnt=1;
	memset(fst,0,sizeof(fst));//初始化
	scanf("%d%d",&n,&m);
	set<int> pqr;//储存那些节点没有被遍历到
	long long sum=0,leftsum=0;//sum表示总和,leftans表示网络中球的重量之和
	for(int i=1;i<=n;++i) scanf("%d%d%d%d",&nx[i],&ny[i],&nw[i],&r[i]);
	for(int i=1;i<=m;++i) scanf("%d%d%d",&x[i],&y[i],&w[i]),pqr.insert(i);
	s=n+m+1,t=n+m+2;//入点和汇点
	sum+=nw[1];
	for(int i=1;i<=m;++i)
	{
		if(dis(nx[1],ny[1],x[i],y[i])<=r[1])
		{
			pqr.erase(pqr.find(i));
			vis[i]=true;
			sum+=w[i];
		}
	}
	for(int i=2;i<=n;++i) for(int j=1;j<=m;++j) if(dis(nx[i],ny[i],x[j],y[j])<=r[i]) if(pqr.count(j)) pqr.erase(pqr.find(j));
	for(int i=2;i<=n;++i)
	{
		if(sum-nw[i]<0)//如果已经不行了,那就自动忽略
		{
			puts("qaq");
			return;
		}
		adds(s,i,sum-nw[i]),adds(i,s,0);//入点和人连边
	}
	for(int i=2;i<=n;++i)
	{
		for(int j=1;j<=m;++j)
		{
			if(dis(nx[i],ny[i],x[j],y[j])<=r[i])
			adds(i,j+n,0x3f3f3f3f3f3f3f3f),adds(j+n,i,0);//人与球连边
		}
	}
	for(int i=1;i<=m;++i) if((!pqr.count(i))&&(!vis[i])) adds(i+n,t,w[i]),adds(t,i+n,0),leftsum+=w[i];
	n=n+m-1;//ISAP一定要改n的值哦
	int shit=ISAP::output();
//	cout<<shit<<endl;
	if(shit==leftsum) puts("ZQC! ZQC!");
	else puts("qaq");
}
signed main()
{
	int tt;
	cin>>tt;
	while(tt--)
	{
		input();
	}
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值