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

相关文章

认识、理解、分类——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

poj2406(连续重复子串)

题意:判断串s是不是str^n,求str的最大长度。 解题思路:kmp可解,后缀数组的倍增算法超时。next[i]表示在第i位匹配失败后,自动跳转到next[i],所以1到next[n]这个串 等于 n-next[n]+1到n这个串。 代码如下; #include<iostream>#include<algorithm>#include<stdio.h>#include<math.

poj3261(可重复k次的最长子串)

题意:可重复k次的最长子串 解题思路:求所有区间[x,x+k-1]中的最小值的最大值。求sa时间复杂度Nlog(N),求最值时间复杂度N*N,但实际复杂度很低。题目数据也比较水,不然估计过不了。 代码入下: #include<iostream>#include<algorithm>#include<stdio.h>#include<math.h>#include<cstring

poj 2976 分数规划二分贪心(部分对总体的贡献度) poj 3111

poj 2976: 题意: 在n场考试中,每场考试共有b题,答对的题目有a题。 允许去掉k场考试,求能达到的最高正确率是多少。 解析: 假设已知准确率为x,则每场考试对于准确率的贡献值为: a - b * x,将贡献值大的排序排在前面舍弃掉后k个。 然后二分x就行了。 代码: #include <iostream>#include <cstdio>#incl

poj 3104 二分答案

题意: n件湿度为num的衣服,每秒钟自己可以蒸发掉1个湿度。 然而如果使用了暖炉,每秒可以烧掉k个湿度,但不计算蒸发了。 现在问这么多的衣服,怎么烧事件最短。 解析: 二分答案咯。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <c

poj 3258 二分最小值最大

题意: 有一些石头排成一条线,第一个和最后一个不能去掉。 其余的共可以去掉m块,要使去掉后石头间距的最小值最大。 解析: 二分石头,最小值最大。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <c

poj 2594 二分图最大独立集

题意: 求一张图的最大独立集,这题不同的地方在于,间接相邻的点也可以有一条边,所以用floyd来把间接相邻的边也连起来。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <sta

poj 3692 二分图最大独立集

题意: 幼儿园里,有G个女生和B个男生。 他们中间有女生和女生认识,男生男生认识,也有男生和女生认识的。 现在要选出一些人,使得这里面的人都认识,问最多能选多少人。 解析: 反过来建边,将不认识的男生和女生相连,然后求一个二分图的最大独立集就行了。 下图很直观: 点击打开链接 原图: 现图: 、 代码: #pragma comment(