【BZOJ2663】灵魂宝石 [二分]
灵魂宝石
Time Limit: 5 Sec Memory Limit: 128 MB
[Submit][Status][Discuss]
Description
“作为你们本体的灵魂,为了能够更好的运用魔法,被赋予了既小巧又安全的外形”
我们知道,魔法少女的生命被存放于一个称为灵魂宝石(Soul Gem)的装置内。
而有时,当灵魂宝石与躯体的距离较远时,魔法少女就无法控制自己的躯体了。
在传说中,魔法少女 Abel仅通过推理就得到了这个现象的一般法则,被称为Abel定理:
存在宇宙常量 R(是一个非负实数,或正无穷) ,被称为灵魂宝石常量,量纲为空间度量(即:长度)。
如果某个魔法少女的灵魂宝石与她的躯体的距离严格超过 R,则她一定无法控制自己的躯体;如果这个距离严格小于 R,则她一定可以控制自己的躯体。 (这里的距离指平面的 Euclid距离。)
注意:该定理不能预言距离刚好为 R 的情形。
可能存在魔法少女 A 和 B,她们离自己的灵魂宝石的距离都恰好为 R,但是 A可以控制自己的躯体,而 B 不可以。
现在这个世界上再也没有魔法少女了,但是我们却对这个宇宙常量感兴趣。
我们只能通过之前的世界遗留下来的数据来确定这个常量的范围了。
每一组数据包含以下信息:
·一共有N 个魔法少女及她们的灵魂宝石,分别编号为 1~N。
·这 N个魔法少女所在的位置是(Xi, Yi)。
·这 N个灵魂宝石所在的位置是(xi, yi)。
·此时恰好有 K个魔法少女能够控制自己的躯体。
需要注意的是:
1. 我们认为这个世界是二维的 Euclid 空间。
2. 魔法少女与灵魂宝石之间的对应关系是未知的。
3. 我们不知道是具体是哪 K个魔法少女能够控制自己的躯体。
根据以上信息,你需要确定灵魂宝石常量 R可能的最小值 Rmin 和最大值 Rmax。
Input
第一行包两个整数:N、K。
接下来 N行,每行包含两个整数:Xi , Yi ,由空格隔开。
再接下来N 行,每行包含两个整数:xi , yi ,由空格隔开。
Output
输出两个量:Rmin、Rmax,中间用空格隔开。
Rmin 一定是一个非负实数,四舍五入到小数点后两位。
Rmax 可能是非负实数,或者是正无穷:
如果是非负实数,四舍五入到小数点后两位;
如果是正无穷,输出“+INF”(不包含引号)。
Sample Input
1 0
4 0
0 0
4 4
Sample Output
HINT
对于100%的数据:
1 ≤ N ≤ 50,
0 ≤ K ≤ N,
-1000 ≤ xi, yi , Xi , Yi ≤ 1000。
Main idea
有n个人匹配n个宝石,每个人和宝石有一个坐标,R为自己给定的值,如果在平面内人和宝石的距离<R则一定匹配,距离=R可取可不取,距离>R则一定无法取,求使得可以取到k个匹配的R的最小值和最大值。
Solution
求最小值最大值,想到了二分答案,然后我们可以直观地看出可以使用二分图匹配来进行求匹配问题,二分一个R,如果人和宝石的距离<=R则连边,判断是否可行,这样我们可以求出最小的R。
发现最大的R无法这么取,因为可能有距离=R的情况,所以我们反向思考,考虑枚举一个R,距离>=R的连边,判断是否有<n-k个无法匹配,则可以求得R的最大值。
Code
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
using namespace std; const int ONE=; int n,k;
double l,mid,r;
double x,y;
int vis[ONE];
int f[ONE][ONE],my[ONE];
double X[ONE],Y[ONE];
double dist[ONE][ONE]; int get()
{
int res,Q=; char c;
while( (c=getchar())< || c>)
if(c=='-')Q=-;
if(Q) res=c-;
while((c=getchar())>= && c<=)
res=res*+c-;
return res*Q;
} int find(int i)
{
for(int j=;j<=n;j++)
{
if(f[i][j] && !vis[j])
{
vis[j]=;
if(!my[j] || find(my[j]))
{
my[j]=i;
return ;
}
}
}
return ;
} double Getdis(double x1,double y1,double x2,double y2)
{
return sqrt( (x1-x2)*(x1-x2) + (y1-y2)*(y1-y2) );
} int Check_first(double x)
{
memset(my,,sizeof(my));
for(int i=;i<=n;i++)
for(int j=;j<=n;j++)
{
f[i][j]=(dist[i][j]<=x);
} int Ans=;
for(int i=;i<=n;i++)
{
memset(vis,,sizeof(vis));
if(find(i)) Ans++;
} return Ans>=k;
} int Check_second(double x)
{
memset(my,,sizeof(my));
for(int i=;i<=n;i++)
for(int j=;j<=n;j++)
{
f[i][j]=(dist[i][j]>=x);
} int Ans=;
for(int i=;i<=n;i++)
{
memset(vis,,sizeof(vis));
if(find(i)) Ans++;
} return Ans<=n-k-;
} int main()
{
n=get(); k=get();
for(int i=;i<=n;i++)
{
scanf("%lf %lf",&X[i],&Y[i]);
} for(int i=;i<=n;i++)
{
scanf("%lf %lf",&x,&y);
for(int j=;j<=n;j++)
dist[j][i]=Getdis(X[j],Y[j],x,y);
} l=0.0; r=3500.0; while(l<r-0.001)
{
mid=(l+r)/2.0;
if(Check_first(mid)) r=mid;
else l=mid;
}
if(Check_first(l)) printf("%.2lf ",l);
else printf("%.2lf ",r); l=0.0; r=3500.0; while(l<r-0.001)
{
mid=(l+r)/2.0;
if(Check_second(mid)) r=mid;
else l=mid;
} double ans;
if(Check_second(r)) ans=r;
else ans=l; if(fabs(ans-3500.0)<=0.01) printf("+INF");
else printf("%.2lf",ans);
}
【BZOJ2663】灵魂宝石 [二分]的更多相关文章
- bzoj2663: [Beijing wc2012]灵魂宝石(二分+匈牙利)
2663: [Beijing wc2012]灵魂宝石 题目:传送门 题解: 又是一道卡精度的题目. 很容易就可以看出单调性啊,如果R越大,选的人就越多,R越小,选的人就越少. 那最小值就直接搞咯. 那 ...
- BZOJ2663 [Beijing wc2012]灵魂宝石
Description 平面中有\(n\)个黑点和\(n\)个白点.这些点组成\(n\)对,但是你不知道它们的对应关系.若某队中黑点白点距离\(<R\),则它是好的:\(>R\)则不是好的 ...
- 灵魂宝石 bzoj 2663
灵魂宝石(1s 128MB)soulgem [问题描述] "作为你们本体的灵魂,为了能够更好的运用魔法,被赋予了既小巧又安全的外形" 我们知道,魔法少女的生命被存放于一个称为灵魂宝 ...
- BZOJ 2663: [Beijing wc2012]灵魂宝石
2663: [Beijing wc2012]灵魂宝石 Time Limit: 5 Sec Memory Limit: 128 MBSubmit: 261 Solved: 108[Submit][S ...
- BZOJ1012: [JSOI2008]最大数maxnumber [线段树 | 单调栈+二分]
1012: [JSOI2008]最大数maxnumber Time Limit: 3 Sec Memory Limit: 162 MBSubmit: 8748 Solved: 3835[Submi ...
- BZOJ 2756: [SCOI2012]奇怪的游戏 [最大流 二分]
2756: [SCOI2012]奇怪的游戏 Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 3352 Solved: 919[Submit][Stat ...
- 整体二分QAQ
POJ 2104 K-th Number 时空隧道 题意: 给出一个序列,每次查询区间第k小 分析: 整体二分入门题? 代码: #include<algorithm> #include&l ...
- [bzoj2653][middle] (二分 + 主席树)
Description 一个长度为n的序列a,设其排过序之后为b,其中位数定义为b[n/2],其中a,b从0开始标号,除法取下整. 给你一个长度为n的序列s. 回答Q个这样的询问:s的左端点在[a,b ...
- [LeetCode] Closest Binary Search Tree Value II 最近的二分搜索树的值之二
Given a non-empty binary search tree and a target value, find k values in the BST that are closest t ...
随机推荐
- 关于 SSH Server 的整体设定
# . 关于 SSH Server 的整体设定,包含使用的 port 啦,以及使用的密码演算方式 Port # SSH 预设使用 这个 port,您也可以使用多的 port ! # 亦即重复使用 po ...
- mysqli函数库的使用
综述 1.什么是mysqli PHP-MySQL 函数库是 PHP 操作 MySQL 资料库最原始的扩展库,PHP-MySQLi 的 i 代表 Improvement ,相当于前者的增强版,也包含了相 ...
- 【APUE】Chapter1 UNIX System Overview
这章内容就是“provides a whirlwind tour of the UNIX System from a programmer's perspective”. 其实在看这章内容的时候,已经 ...
- json 处理日期格式
@JsonFormat(pattern="yyyy-MM-dd HH:mm:ss",timezone="GMT+8") private Date createT ...
- Qt 建立带有子项目的工程
刚需,软件需要用到多个子项目 第一步 打开Qt新建子项目工程 如图 在此时鼠标右键,选着新建子项目如图 就是正常的新建项目的步骤,直接上图 完工,可以愉快的撸代码了
- 关于redis一些问题记录
问题一:启动redis时出现警告,使用下列命令(已解决) 问题二:启动时,需要解决的警告(未解决) 问题三:使用自己的配置文件启动redis时,可能会遇到: Could not connect to ...
- Python 套接字的使用 (1)
获取设备名称和IPv4地址 socket.gethostname() socket.gethostbyname(host_name) def print_machine_info(): host_ ...
- windbg*****************************TBD
achieve structure from a simple address Dt address know pending IRP in a module !thread xxxxxx到底能提供哪 ...
- 关于如何利用原生js动态给一个空对象添加属性以及属性值
首先,回忆一下,访问对象属性一共有两种方法:点获取法和方括号获取法.而我们最常用的就是点获取法了.但是当我们遇到需要给对象动态添加属性和属性值时,点获取法好像就不太好用了,尤其是我们不知道属性名的时候 ...
- [剑指Offer] 43.左旋转字符串
题目描述 汇编语言中有一种移位指令叫做循环左移(ROL),现在有个简单的任务,就是用字符串模拟这个指令的运算结果.对于一个给定的字符序列S,请你把其循环左移K位后的序列输出.例如,字符序列S=”abc ...