题目链接:http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=26868

思路:拆点,容量为最多能跳的步数,然后设立一个超级源点,源点与各点两连边,容量为一开始的企鹅数,最后就是枚举汇点了,跑最大流验证即可。

 #include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
#define MAXN 222
#define MAXM 222222
#define inf 1<<30
#define FILL(a,b) memset(a,b,sizeof(a)) struct Edge{
int v,cap,next;
}edge[MAXM]; int n,vs,vt,NV,NE,head[MAXN];
double d; void Insert(int u,int v,int cap)
{
edge[NE].v=v;
edge[NE].cap=cap;
edge[NE].next=head[u];
head[u]=NE++; edge[NE].v=u;
edge[NE].cap=;
edge[NE].next=head[v];
head[v]=NE++;
} int level[MAXN],gap[MAXN];
void bfs(int vt)
{
FILL(level,-);
FILL(gap,);
queue<int>que;
que.push(vt);
level[vt]=;
gap[]++;
while(!que.empty()){
int u=que.front();
que.pop();
for(int i=head[u];i!=-;i=edge[i].next){
int v=edge[i].v;
if(level[v]!=-)continue;
level[v]=level[u]+;
gap[level[v]]++;
que.push(v);
}
}
} int pre[MAXN],cur[MAXN];
int SAP(int vs,int vt)
{
bfs(vt);
memcpy(cur,head,sizeof(head));
int u=pre[vs]=vs,aug=inf,maxflow=;
gap[]=NV;
while(level[vs]<NV){
bool flag=false;
for(int &i=cur[u];i!=-;i=edge[i].next){
int v=edge[i].v;
if(edge[i].cap>&&level[u]==level[v]+){
flag=true;
aug=min(aug,edge[i].cap);
pre[v]=u;
u=v;
if(v==vt){
maxflow+=aug;
for(u=pre[u];v!=vs;v=u,u=pre[u]){
edge[cur[u]].cap-=aug;
edge[cur[u]^].cap+=aug;
}
aug=inf;
}
break;
}
}
if(flag)continue;
int minlevel=NV;
for(int i=head[u];i!=-;i=edge[i].next){
int v=edge[i].v;
if(edge[i].cap>&&level[v]<minlevel){
cur[u]=i;
minlevel=level[v];
}
}
if(--gap[level[u]]==)break;
level[u]=minlevel+;
gap[level[u]]++;
u=pre[u];
}
return maxflow;
} struct Node{
double x,y;
int num,step;
}node[MAXN]; double Get_Dist(int i,int j)
{
double d1=(node[i].x-node[j].x)*(node[i].x-node[j].x);
double d2=(node[i].y-node[j].y)*(node[i].y-node[j].y);
return sqrt(d1+d2);
} void Build(int ed)
{
NE=;
FILL(head,-);
vs=,vt=*n+;
NV=vt+;
for(int i=;i<=n;i++){
Insert(i,i+n,node[i].step);
Insert(vs,i,node[i].num);
}
for(int i=;i<=n;i++){
for(int j=i+;j<=n;j++){
if(Get_Dist(i,j)<=d){
Insert(i+n,j,inf);
Insert(j+n,i,inf);
}
}
}
Insert(ed,vt,inf);
} int main()
{
int _case,ans,sum,flag,t=;
scanf("%d",&_case);
while(_case--){
scanf("%d%lf",&n,&d);
sum=;
for(int i=;i<=n;i++){
scanf("%lf%lf%d%d",&node[i].x,&node[i].y,&node[i].num,&node[i].step);
sum+=node[i].num;
}
flag=;
printf("Case %d:",t++);
//枚举集中点
for(int i=;i<=n;i++){
Build(i);
if(SAP(vs,vt)==sum)flag=,printf(" %d",i-);
}
if(!flag)printf(" -1");
puts("");
}
return ;
}

最新文章

  1. Scala元组
  2. C++ 定义全局数组
  3. js面向对象的实现(example 二)
  4. Python顺序集合之 tuple
  5. Codeforces 716C[数论][构造]
  6. linux+asp.net core+nginx四层负载均衡
  7. json转换(c#后台生成json的方法)
  8. python应用之文件属性浏览
  9. InitParam与ContextParm的异同
  10. select下拉菜单反显不可改动,且submit能够提交数据
  11. git用法-打补丁
  12. Apache、IIS、Nginx等绝大多数web服务器,都不允许静态文件响应POST请求,否则会返回“HTTP/1.1 405 Method not allowed”错误。
  13. NVIDIA-SMI系列命令总结
  14. redis参数说明
  15. mysql表关联
  16. pandas.DataFrame
  17. pta l2-18(多项式A除以B)
  18. ICMP协议、DNS、ARP协议、ping、DHCP协议
  19. 百度前端技术学院-task2.18-2.19源码以及个人总结
  20. UNIX高级环境编程(11)进程控制(Process Control)- 进程快照,用户标识符,进程调度

热门文章

  1. C++ 中宏的使用 --来自:http://blog.csdn.net/hgl868/article/details/7058906
  2. emmet-vim
  3. 微信公开课发布微信官方教程:教你用好微信JS-SDK接口
  4. 前端与php的sublime text3常用插件
  5. js矩阵菜单或3D立体预览图片效果
  6. function foo(){}、(function(){})、(function(){}())等函数区别分析
  7. linux 查看系统信息命令(比较全)
  8. nginx lua处理图片
  9. win7系统扩展双屏幕时,如何在两个屏幕下都显示任务栏
  10. 20个很有用的PHP类库