题目链接:

http://acm.hdu.edu.cn/showproblem.php?pid=5046

题意:

给定n个城市的坐标,要在城市中建k个飞机场,使城市距离最近的飞机场的最长距离最小,求这个最小距离。

分析:

最小化最大值,显然二分最大距离。然后我们将距离在范围内的两个城市建边,看能否选出k个城市,使得覆盖了所有城市。

将点之间的关系转化成01矩阵的覆盖问题,重复覆盖,建好边套个DLX即可。

看了鸟神博客,这里可以直接将所有距离保存在一个数组中,排序并去重,二分下标即可。这样快了很多很多!

hdu 2295 和这题一个套路,更裸一些。

跳舞链好强大,可惜只会用模板,这个讲的还挺清晰

代码:

/*************************************************************************
> File Name: O.cpp
> Author: jiangyuzhu
> Mail: 834138558@qq.com
> Created Time: Fri 08 Jul 2016 03:31:58 PM CST
************************************************************************/
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
#define sal(n) scanf("%I64d", &(n))
#define sa(n) scanf("%d", &(n))
typedef long long ll;
const int maxn = 3600 + 5, maxc = 60 + 5, maxr = 60 + 5, inf = 0x3f3f3f3f;
struct Node{ll x; ll y;}city[maxn];
int L[maxn], R[maxn], D[maxn], U[maxn], C[maxn];
int S[maxc], H[maxr], size;
int n, m, k;
ll d[maxr][maxc], t[maxn];
///不需要S域
void Link(int r, int c)
{
S[c]++; C[size] = c;
U[size] = U[c]; D[U[c]] = size;
D[size] = c; U[c] = size;
U[c] = size;
if(H[r] == -1) H[r] = L[size] = R[size] = size;
else {
L[size] = L[H[r]]; R[L[H[r]]] = size;
R[size] = H[r]; L[H[r]] = size;
}
size++;
}
void remove(int c)
{
for(int i = D[c]; i != c; i = D[i])
L[R[i]] = L[i], R[L[i]] = R[i];
}
void resume(int c)
{
for (int i = U[c]; i != c; i = U[i])
L[R[i]] = R[L[i]] = i;
}
int h()
{///用精确覆盖去估算剪枝
int ret = 0;
bool vis[maxc];
memset (vis, false, sizeof(vis));
for(int i = R[0]; i; i = R[i]){
if(vis[i]) continue;
ret++;
vis[i] = true;
for (int j = D[i]; j != i; j = D[j])
for (int k = R[j]; k != j; k = R[k])
vis[C[k]] = true;
}
return ret;
}
int ans;
bool Dance(int a) //具体问题具体分析
{
if(a + h() > k) return 0;
if(!R[0]) return a <= k;
int c = R[0];
for (int i = R[0]; i; i = R[i])
if(S[i] < S[c]) c = i;
for(int i = D[c]; i != c; i = D[i]){
remove(i);
for(int j = R[i]; j != i; j = R[j])
remove(j);
if(Dance(a + 1)) return 1;
for (int j = L[i]; j != i; j = L[j])
resume(j);
resume(i);
}
return 0;
}
void initL(int x)///col is 1~x,row start from 1
{
for (int i = 0; i <= x; ++i){
S[i] = 0;
D[i] = U[i] = i;
L[i+1] = i; R[i] = i+1;
}///对列表头初始化
R[x] = 0;
size = x + 1;///真正的元素从m+1开始
memset (H, -1, sizeof(H));
///mark每个位置的名字
}
ll dist(int i, int j)
{
ll d = abs(city[i].x - city[j].x) + abs(city[i].y - city[j].y);
return d;
}
bool judge(ll mid)
{
initL(n);
ans = 0x3f3f3f3f;
for(int i = 1; i <= n; i++){
for(int j = 1; j <= n; j++){
if(d[i][j] <= mid){
Link(i, j);
}
}
}
return Dance(0);
}
int main (void)
{
int T;sa(T);
for(int tt = 1; tt <= T; tt++){
sa(n);sa(k);
ll maxd = 0;
int cnt = 0;
for(int i = 1; i <= n; i++){
scanf("%I64d%I64d", &city[i].x, &city[i].y);
}
for(int i = 1; i <= n; i++){
for(int j = 1; j <= n; j++){
d[i][j] = dist(i, j);
t[cnt++] = d[i][j];
}
}
sort(t, t + cnt);
int ncnt = unique(t, t + cnt) - t;
ll l = -1, r = ncnt;
while(r - l > 1){
ll mid = (l + r) / 2;
if(judge(t[mid])) r = mid;
else l = mid;
}
printf("Case #%d: %I64d\n", tt, t[r]);
}
return 0;
}

HDU 5046 Airport【DLX重复覆盖】的更多相关文章

  1. hdu 5046 Airport 二分+重复覆盖

    题目链接 给n个点, 定义两点之间距离为|x1-x2|+|y1-y2|. 然后要选出k个城市建机场, 每个机场可以覆盖一个半径的距离. 求在选出点数不大于k的情况下, 这个半径距离的最大值. 二分半径 ...

  2. HDU 2295.Radar (DLX重复覆盖)

    2分答案+DLX判断可行 不使用的估计函数的可重复覆盖的搜索树将十分庞大 #include <iostream> #include <cstring> #include < ...

  3. HDU 5046 Airport(dlx)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5046 题意:n个城市修建m个机场,使得每个城市到最近进场的最大值最小. 思路:二分+dlx搜索判定. ...

  4. HDU 5046 Airport(DLX反复覆盖)

    HDU 5046 Airport 题目链接 题意:给定一些机场.要求选出K个机场,使得其它机场到其它机场的最大值最小 思路:二分+DLX反复覆盖去推断就可以 代码: #include <cstd ...

  5. (中等) HDU 3335 , DLX+重复覆盖。

    Description As we know,the fzu AekdyCoin is famous of math,especially in the field of number theory. ...

  6. HDU 3957 Street Fighter (最小支配集 DLX 重复覆盖+精确覆盖 )

    DLX经典题型,被虐惨了…… 建一个2*N行3*N列的矩阵,行代表选择,列代表约束.前2*N列代表每个人的哪种状态,后N列保证每个人至多选一次. 显然对手可以被战胜多次(重复覆盖),每个角色至多选择一 ...

  7. hdu 2295 dlx重复覆盖+二分答案

    题目大意: 有一堆雷达工作站,安放至多k个人在这些工作站中,找到一个最小的雷达监控半径可以使k个工作人所在的雷达工作站覆盖所有城市 二分半径的答案,每次利用dlx的重复覆盖来判断这个答案是否正确 #i ...

  8. [ACM] HDU 2295 Radar (二分法+DLX 重复覆盖)

    Radar Problem Description N cities of the Java Kingdom need to be covered by radars for being in a s ...

  9. HDU 2295 Radar (重复覆盖)

    Radar Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submi ...

随机推荐

  1. STM32HAL学习博客

    https://www.cnblogs.com/wt88/category/1297945.html

  2. poj 1742 多重背包问题 dp算法

    题意:硬币分别有 A1.....An种,每种各有C1......Cn个,问组成小于m的有多少种 思路:多重背包问题 dp[i][j]表示用前i种硬币组成j最多剩下多少个  dp=-1的表示凑不齐 dp ...

  3. Linux学习-核心与核心模块

    谈完了整个开机的流程,您应该会知道,在整个开机的过程当中,是否能够成功的驱动我们主机的硬 件配备, 是核心 (kernel) 的工作!而核心一般都是压缩文件,因此在使用核心之前,就得要将他解 压缩后, ...

  4. flask-mail(qq邮箱)

    from flask_mail import Mail,Message app.config['MAIL_SERVER']='smtp.qq.com' app.config['MAIL_PORT'] ...

  5. HDU 5236 Article 期望

    题意: 你现在要打\(n\)个字符,但是程序随时可能会崩溃. 你可以在恰当的时机按下 \(Ctrl-S\)键,崩溃后,会从最后一次保存的情况继续开始打字. 具体是这样的: 在每个第\(i-0.1s(i ...

  6. Django two

    http://www.cnblogs.com/yuanchenqi/articles/6083427.html Django: 1.安装Django pip install  django 2.创建p ...

  7. PYday16&17-设计模式\选课系统习题

    1.设计模式:对程序做整体得规划设计,这样做是为了更好的实现功能,使代码的可扩展性更好有27种常见的设计模式.流行的设计模式参考书:GoF设计模式.大话设计模式设计模式是为了更好的实现模块间的解耦,便 ...

  8. 蓝桥杯Java输入输出相关

    转载自:http://blog.csdn.net/Chen_Tongsheng/article/details/53354169 一.注意点 1. 类名称必须采用public class Main方式 ...

  9. Oracle 了解 DDL 操作与 REDO 的关系

    目录 了解 DDL 操作与 REDO 的关系 DDL是否会产生REDO 通过 10046 trace 来分析create 和drop 如果drop失败,redo的变化 了解 DDL 操作与 REDO ...

  10. Notepad++ 使用正则表达式查找替换字符串

        最近在改一些别人写的PHP代码,由于之前的代码可维护性比较差,因此工作量比较多,因此想找一些高效的方式改一些有规律的代码. 比如: 我想将一些类似$rowss[MaxGetCash_num], ...