hdu3656Fire station(DLX重复覆盖 + 二分)
题目大意:一个城市n个点,现在要建m个消防站,消防站建在给定的n个点中。求建m个消防站后,m个消防站要覆盖所有的n个点的覆盖半径最小。
题目分析:重复覆盖问题,DLX解决。不过要求覆盖半径最小,需要二分。虽然给的范围并不大,DLX毕竟还是暴力搜索,而且精度有6位小数,因此直接二分距离的话会TLE!解决方案是将图中任意2点的距离记录下来,去重后二分已知的距离。因为消防站建在给定的n个点中,那么最小覆盖半径一定在任意2点距离中产生。
DLX搜索的时候,一般习惯删除的时候从左往右,不过效率却不一定高。比如这题,从左向右删除和从右向左删除,跑的时间至少差了1s以上!可见用dancing links搜索的时候姿势还是很重要的。有时候换个姿势也许效率更高~
详情请见代码:
#include <iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int N = 51;
const int M = 30001;
const double eps = 1e-7; double dis[N][N];
double tle[M];
int s[M],h[M],u[M],d[M],l[M],r[M],col[M],row[M];
int point[N][2];
int m,n,num,len;
int fs;
double getdis(int i,int j)
{
return sqrt((double)(point[i][0] - point[j][0])*(point[i][0] - point[j][0])
+(double)(point[i][1] - point[j][1])*(point[i][1] - point[j][1]));
}
void read()
{
scanf("%d%d",&n,&m);
int i,j;
len = 1;
tle[len ++] = 0;
for(i = 1;i <= n;i ++)
{
dis[i][i] = 0.0;
scanf("%d%d",&point[i][0],&point[i][1]);
for(j = 1;j < i;j ++)
dis[i][j] = dis[j][i] = getdis(i,j),tle[len ++] = dis[i][j];
}
}
void init()
{
memset(h,0,sizeof(h));
memset(s,0,sizeof(s));
for(int i = 0;i <= n;i ++)
{
u[i] = d[i] = i;
l[i] = (i + n) % (n + 1);
r[i] = (i + 1) % (n + 1);
}
num = n + 1;
}
void add(int i,int j)
{
if(h[i])
{
r[num] = h[i];
l[num] = l[h[i]];
r[l[num]] = l[r[num]] = num;
}
else
h[i] = l[num] = r[num] = num;
s[j] ++;
u[num] = u[j];
d[num] = j;
d[u[num]] = num;
u[j] = num;
col[num] = j;
row[num] = i;
num ++;
}
void build(double md)
{
int i,j;
init();
for(i = 1;i <= n;i ++)
for(j = 1;j <= n;j ++)
if(md - dis[i][j] > -eps)
add(i,j);
}
void remove(int c)
{
for(int i = d[c];i != c;i = d[i])
l[r[i]] = l[i],r[l[i]] = r[i],s[col[i]] --;
}
void resume(int c)
{
for(int i = u[c];i != c;i = u[i])
l[r[i]] = r[l[i]] = i,s[col[i]] ++;
}
int A()
{
int i,j,k,ret = 0;
bool vis[N];
memset(vis,false,sizeof(vis));
for(i = l[0];i;i = l[i])
{
if(vis[i] == false)
{
vis[i] = true;
ret ++;
for(j = d[i];j != i;j = d[j])
for(k = r[j];k != j;k = r[k])
vis[col[k]] = true;
}
}
return ret;
}
void dfs(int k)
{
if(k + A() >= fs)
return;
int i,j;
if(!r[0])
{
fs = min(fs,k);
return;
}
int mn = 1000000;
int c;
for(i = l[0];i;i = l[i])
{
if(mn > s[i])
{
mn = s[i];
c = i;
}
}
for(i = d[c];i != c;i = d[i])
{
remove(i);
for(j = l[i];j != i;j = l[j])
{
remove(j);
}
dfs(k + 1);
for(j = r[i];j != i;j = r[j])
{
resume(j);
}
resume(i);
}
}
void solve()
{
int la,ra,mid,ans;
sort(tle + 1,tle + len);
int i = len;
int j;
len = 2;
for(j = 2;j < i;j ++)
if(fabs(tle[j] - tle[j - 1]) > eps)
tle[len ++] = tle[j];
len --;
la = 1;ra = len;
while(la <= ra)
{
mid = (la + ra)>>1;
build(tle[mid]);
fs = M;
dfs(0);
if(fs <= m)
{
ans = mid;
ra = mid - 1;
}
else
la = mid + 1;
}
printf("%f\n",tle[ans]);
}
int main()
{
int _;
scanf("%d",&_);
while(_ --)
{
read();
solve();
}
return 0;
}
//1953MS 636K
hdu3656Fire station(DLX重复覆盖 + 二分)的更多相关文章
- hdu 2295 dlx重复覆盖+二分答案
题目大意: 有一堆雷达工作站,安放至多k个人在这些工作站中,找到一个最小的雷达监控半径可以使k个工作人所在的雷达工作站覆盖所有城市 二分半径的答案,每次利用dlx的重复覆盖来判断这个答案是否正确 #i ...
- (中等) HDU 2295 , DLX+重复覆盖+二分。
Description N cities of the Java Kingdom need to be covered by radars for being in a state of war. S ...
- hdu2295(重复覆盖+二分)
题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=2295 题意::一个国家有n个城市,有m个地方可以建造雷达,最多可以建K个雷达(K>=1 & ...
- hdu5046(重复覆盖+二分)
题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=5046 题意:要在n个城市里建造不超过k个机场覆盖所有城市,问机场城市之间最大距离最小为多少. 分析:二 ...
- [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 ...
- HDU 5046 Airport【DLX重复覆盖】
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5046 题意: 给定n个城市的坐标,要在城市中建k个飞机场,使城市距离最近的飞机场的最长距离最小,求这 ...
- HDU5046 Airport dancing links 重复覆盖+二分
这一道题和HDU2295是一样 是一个dancing links重复覆盖解决最小支配集的问题 在给定长度下求一个最小支配集,只要小于k就行 然后就是二分答案,每次求最小支配集 只不过HDU2295是浮 ...
- HDU 2295.Radar (DLX重复覆盖)
2分答案+DLX判断可行 不使用的估计函数的可重复覆盖的搜索树将十分庞大 #include <iostream> #include <cstring> #include < ...
- (中等) HDU 3335 , DLX+重复覆盖。
Description As we know,the fzu AekdyCoin is famous of math,especially in the field of number theory. ...
随机推荐
- MySQL存储过程(一)
1.1 CREATE PROCEDURE (创建) CREATE PROCEDURE存储过程名 (参数列表) BEGIN SQL语句代码块 END 注意: 由括号包围的参数列必须总是存在.如果没有参 ...
- ThinkPHP学习 volist标签高级应用之多重嵌套循环、隔行变色(转)
Action代码: public function index(){ $prod = I("get.prod_en"); $id = I("get.id", 0 ...
- MySql Connector/Net Mysql like 搜索中文的问题(c#和asp.net连接mysql)
Connector/Net 6.9.8 选择.net/mono即可,不需要安装. 将对应版本的MySql.Data.dll复制到bin目录下即可使用 http://dev.mysql.com/down ...
- C#多线程lock解决数据同步
1.代码实例: public class ThreadTest4 { public static void Init() { //多个线程修改同一个值,使用lock锁解决并发 ; i < ; i ...
- Ajax数据格式,html,xml,json
1. 2. 3. 4. 5. 6. 7. 8. 9.
- PHP MySQL 读取数据
PHP MySQL 读取数据 从 MySQL 数据库读取数据 SELECT 语句用于从数据表中读取数据: SELECT column_name(s) FROM table_name 如需学习更多关于 ...
- Git对于单个文件的分批提交方式的使用
很多时候,对于一个大的文件,可能有的同学改完之后不想一次提交,想分批提交.但这个时候由于git add的机制往往add之后就是整个一个文件被放到stage区了,这个时候肯定会想能不能对一个文件可以进行 ...
- 逻辑很重要:一句sql语句的事,自己却想了半天,绕了个大弯子
问题:系统升级后审核认证信息分别写入两个表,现在需要链接用户表和相应的新旧审核表获取字段值? 钻进胡同里:一直纠结于升级之后的会员信息从新表查,升级之前的数据从旧表查,纠结于根据时间戳分条件判断, 其 ...
- CSS实现div居中
<!DOCTYPE html> <html> <head> <meta charset="utf-8" /> <title&g ...
- uva 10382 - Watering Grass(区域覆盖问题)
Sample Input 8 20 2 5 3 4 1 1 2 7 2 10 2 13 3 16 2 19 4 3 10 1 3 5 9 3 6 1 3 10 1 5 3 1 1 9 1 Sample ...