这题题意描述不是很清楚啊,所以我找了个有权限的人把题面改了改,应该还是比较清楚了。
感觉这道题挺妙的,就来写一篇题解。
首先,根据贪心思想,我们会将 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 eati≤sum−wi。
那么既然如此,我们就可以考虑建一个网络,对于第 i i i 个人,将他与他半径以内能吃的东西连一条容量为极大值的边,随便找一个入点 s s s 和汇点 t t t。显然,对于第 i i i 个人, s s s 就向他连一条容量为 s u m − w i sum-w_i sum−wi 的边(显然,如果更大,那么第一个点就不是最大的了),并对于任意一个球,向汇点连一条边,容量为 w i w_i wi(每个吃的只能被吃一次)。由于每个能吃的球,必须被吃完,所以必须要保证这个网络的最大流是满流(即最大流等于网络中每个球的重量之和),否则,球就吃不完,不满足题意,故第一个节点也就成不了最大的重量。
注意事项:有些球对于所有的人他都吃不到,要自动忽略这些点。
#include
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<
if(shit==leftsum) puts("ZQC! ZQC!");
else puts("qaq");
}
signed main()
{
int tt;
cin>>tt;
while(tt--)
{
input();
}
return 0;
}