2014 ACM/ICPC Asia Regional Shanghai Online E - Airport —— 二分+舞蹈链(DLX)重复覆盖

本文主要是介绍2014 ACM/ICPC Asia Regional Shanghai Online E - Airport —— 二分+舞蹈链(DLX)重复覆盖,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

This way

题意:

给你n个点,让你在其中选k个点作为特殊点,使得所有点到其中距离自己最近的特殊点的距离最大值最小,求这个值

题解:

n只有60,那么翻译一下这个就是重复覆盖问题。
那么我们只需要二分一下答案,将所有小于等于mid的值加到舞蹈链中,再跑一下即可。
注意其中的优化:
由与deep就是当前使用的点数,那么当deep>k的时候return,注意不能直接做这个判断:

if(deep+leave()>k)return;

这个判断就是说你当前加的点的基础上最少需要多少点才能完全覆盖
这个其实是很费时间的,它应当加在所有剪枝后面。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=4005;
int n;
int cnt;
int l[N],r[N],u[N],d[N],col[N],row[N];//每个点的左右上下指针,所在行列
int h[N];//每行的头结点
int s[N];//每列的节点数
int ansk[N];//选了那些集合
void init(int m){//每行m个元素for(int i=0;i<=m;i++){r[i]=i+1;l[i]=i-1;u[i]=d[i]=i;}r[m]=0;l[0]=m;memset(h,-1,sizeof(h));memset(s,0,sizeof(s));cnt=m+1;
}//初始化
void add(int R,int C){//R行C列插入点s[C]++;row[cnt]=R;col[cnt]=C;u[cnt]=C;d[cnt]=d[C];u[d[C]]=cnt;d[C]=cnt;if(h[R]==-1)h[R]=r[cnt]=l[cnt]=cnt;//该行没有点,直接加入else{r[cnt]=h[R];l[cnt]=l[h[R]];r[l[h[R]]]=cnt;l[h[R]]=cnt;}cnt++;//return;
}
void del(int c){for(int i=d[c];i!=c;i=d[i]){l[r[i]]=l[i],r[l[i]]=r[i];}
}//删除c列
void res(int c){for(int i=u[c];i!=c;i=u[i]){l[r[i]]=i,r[l[i]]=i;}
}//恢复c列
int ans;
ll dis[65][65];
int k;
int leave(){int ans=0;bool vis[n];int i,j,k;memset(vis,0,sizeof(vis));for(i=r[0];i!=0;i=r[i]){if(vis[i]==0){vis[i]=1;ans++;for(j=d[i];j!=i;j=d[j]){for(k=r[j];k!=j;k=r[k])vis[col[k]]=1;}}}return ans;
}//最优情况还要多少个集合
void dance(int deep){if(ans<=k)return ;if(deep>k)return ;if(deep+leave()>k)return;//剪枝if(r[0]==0){ans=deep;return;}int c=r[0];register int i,j;for(i=r[0];i!=0;i=r[i])if(s[i]<s[c])c=i;//找到点最少的列for(i=d[c];i!=c;i=d[i]){ansk[deep]=row[i];del(i);for(j=r[i];j!=i;j=r[j])del(j);dance(deep+1);if(ans<=k)return ;for(j=l[i];j!=i;j=l[j])res(j);res(i);}return;
}
ll x[65],y[65];
int main()
{int t,cas=0;scanf("%d",&t);while(t--){scanf("%d%d",&n,&k);for(int i=1;i<=n;i++)scanf("%lld%lld",&x[i],&y[i]);ll mi=4e9,mx=0;init(n);for(int i=1;i<=n;i++){for(int j=i;j<=n;j++){dis[i][j]=(x[i]-x[j]>=0?x[i]-x[j]:x[j]-x[i])+(y[i]-y[j]>=0?y[i]-y[j]:y[j]-y[i]);if(i!=j)mi=min(mi,dis[i][j]),mx=max(mx,dis[i][j]);if(dis[i][j]==0)add(i,j);}}ans=1e9;dance(0);if(ans<=k){printf("Case #%d: %lld\n",++cas,0);continue;}ll l=mi,r=mx,mid;mi=4e9;while(r>=l){ans=1e9;mid=(l+r)/2;init(n);for(int i=1;i<=n;i++){for(int j=i;j<=n;j++){if(dis[i][j]<=mid){add(i,j);if(i-j)add(j,i);}}}dance(0);if(ans<=k)r=mid-1,mi=mid;elsel=mid+1;}printf("Case #%d: %lld\n",++cas,mi);}return 0;
}

这篇关于2014 ACM/ICPC Asia Regional Shanghai Online E - Airport —— 二分+舞蹈链(DLX)重复覆盖的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/881187

相关文章

Java覆盖第三方jar包中的某一个类的实现方法

《Java覆盖第三方jar包中的某一个类的实现方法》在我们日常的开发中,经常需要使用第三方的jar包,有时候我们会发现第三方的jar包中的某一个类有问题,或者我们需要定制化修改其中的逻辑,那么应该如何... 目录一、需求描述二、示例描述三、操作步骤四、验证结果五、实现原理一、需求描述需求描述如下:需要在

Redis 多规则限流和防重复提交方案实现小结

《Redis多规则限流和防重复提交方案实现小结》本文主要介绍了Redis多规则限流和防重复提交方案实现小结,包括使用String结构和Zset结构来记录用户IP的访问次数,具有一定的参考价值,感兴趣... 目录一:使用 String 结构记录固定时间段内某用户 IP 访问某接口的次数二:使用 Zset 进行

Spring Boot 整合 ShedLock 处理定时任务重复执行的问题小结

《SpringBoot整合ShedLock处理定时任务重复执行的问题小结》ShedLock是解决分布式系统中定时任务重复执行问题的Java库,通过在数据库中加锁,确保只有一个节点在指定时间执行... 目录前言什么是 ShedLock?ShedLock 的工作原理:定时任务重复执行China编程的问题使用 Shed

Oracle数据库使用 listagg去重删除重复数据的方法汇总

《Oracle数据库使用listagg去重删除重复数据的方法汇总》文章介绍了在Oracle数据库中使用LISTAGG和XMLAGG函数进行字符串聚合并去重的方法,包括去重聚合、使用XML解析和CLO... 目录案例表第一种:使用wm_concat() + distinct去重聚合第二种:使用listagg,

vue解决子组件样式覆盖问题scoped deep

《vue解决子组件样式覆盖问题scopeddeep》文章主要介绍了在Vue项目中处理全局样式和局部样式的方法,包括使用scoped属性和深度选择器(/deep/)来覆盖子组件的样式,作者建议所有组件... 目录前言scoped分析deep分析使用总结所有组件必须加scoped父组件覆盖子组件使用deep前言

MySQL中删除重复数据SQL的三种写法

《MySQL中删除重复数据SQL的三种写法》:本文主要介绍MySQL中删除重复数据SQL的三种写法,文中通过代码示例讲解的非常详细,对大家的学习或工作有一定的帮助,需要的朋友可以参考下... 目录方法一:使用 left join + 子查询删除重复数据(推荐)方法二:创建临时表(需分多步执行,逻辑清晰,但会

Python绘制土地利用和土地覆盖类型图示例详解

《Python绘制土地利用和土地覆盖类型图示例详解》本文介绍了如何使用Python绘制土地利用和土地覆盖类型图,并提供了详细的代码示例,通过安装所需的库,准备地理数据,使用geopandas和matp... 目录一、所需库的安装二、数据准备三、绘制土地利用和土地覆盖类型图四、代码解释五、其他可视化形式1.

认识、理解、分类——acm之搜索

普通搜索方法有两种:1、广度优先搜索;2、深度优先搜索; 更多搜索方法: 3、双向广度优先搜索; 4、启发式搜索(包括A*算法等); 搜索通常会用到的知识点:状态压缩(位压缩,利用hash思想压缩)。

hdu2241(二分+合并数组)

题意:判断是否存在a+b+c = x,a,b,c分别属于集合A,B,C 如果用暴力会超时,所以这里用到了数组合并,将b,c数组合并成d,d数组存的是b,c数组元素的和,然后对d数组进行二分就可以了 代码如下(附注释): #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<que

hdu2289(简单二分)

虽说是简单二分,但是我还是wa死了  题意:已知圆台的体积,求高度 首先要知道圆台体积怎么求:设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ