灵魂宝石

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

  2 1
  1 0
  4 0
  0 0
  4 4

Sample Output

  1.00 5.00

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】灵魂宝石 [二分]的更多相关文章

  1. bzoj2663: [Beijing wc2012]灵魂宝石(二分+匈牙利)

    2663: [Beijing wc2012]灵魂宝石 题目:传送门 题解: 又是一道卡精度的题目. 很容易就可以看出单调性啊,如果R越大,选的人就越多,R越小,选的人就越少. 那最小值就直接搞咯. 那 ...

  2. BZOJ2663 [Beijing wc2012]灵魂宝石

    Description 平面中有\(n\)个黑点和\(n\)个白点.这些点组成\(n\)对,但是你不知道它们的对应关系.若某队中黑点白点距离\(<R\),则它是好的:\(>R\)则不是好的 ...

  3. 灵魂宝石 bzoj 2663

    灵魂宝石(1s 128MB)soulgem [问题描述] "作为你们本体的灵魂,为了能够更好的运用魔法,被赋予了既小巧又安全的外形" 我们知道,魔法少女的生命被存放于一个称为灵魂宝 ...

  4. BZOJ 2663: [Beijing wc2012]灵魂宝石

    2663: [Beijing wc2012]灵魂宝石 Time Limit: 5 Sec  Memory Limit: 128 MBSubmit: 261  Solved: 108[Submit][S ...

  5. BZOJ1012: [JSOI2008]最大数maxnumber [线段树 | 单调栈+二分]

    1012: [JSOI2008]最大数maxnumber Time Limit: 3 Sec  Memory Limit: 162 MBSubmit: 8748  Solved: 3835[Submi ...

  6. BZOJ 2756: [SCOI2012]奇怪的游戏 [最大流 二分]

    2756: [SCOI2012]奇怪的游戏 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 3352  Solved: 919[Submit][Stat ...

  7. 整体二分QAQ

    POJ 2104 K-th Number 时空隧道 题意: 给出一个序列,每次查询区间第k小 分析: 整体二分入门题? 代码: #include<algorithm> #include&l ...

  8. [bzoj2653][middle] (二分 + 主席树)

    Description 一个长度为n的序列a,设其排过序之后为b,其中位数定义为b[n/2],其中a,b从0开始标号,除法取下整. 给你一个长度为n的序列s. 回答Q个这样的询问:s的左端点在[a,b ...

  9. [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 ...

随机推荐

  1. 深度学习:激活函数的比较和优缺点,sigmoid,tanh,relu

    https://blog.csdn.net/u011684265/article/details/78039280

  2. HMM相关文章索引

    HMM相关文章索引 1条回复 HMM系列文章是52nlp上访问量较高的一批文章,这里做个索引,方便大家参考. HMM学习 HMM学习最佳范例一:介绍 HMM学习最佳范例二:生成模式 HMM学习最佳范例 ...

  3. 【jQuery】 Ajax

    [jQuery] Ajax $.ajax({ type: "Post", // 发包方式 cache: false, // 是否缓存 contentType: "appl ...

  4. Android应用开发中的夜间模式实现(一)

    前言 在应用开发中会经常遇到要求实现夜间模式或者主题切换具体例子如下,我会先讲解第一种方法. 夜间模式 知乎 网易新闻 沪江开心词场 Pocket 主题切换 腾讯QQ 新浪微博 我今天主要是详述第一种 ...

  5. win10 java环境变量配置

    首先,你应该已经安装了 Java 的 JDK 了(如果没有安装JDK,请跳转到此网址:http://www.oracle.com/technetwork/java/javase/downloads/i ...

  6. 初探 Qt Opengl【2】

    最近在研究QOPengl QGraphicsView QGraphicsItemQGraphicsScene不过也只是皮毛,也不是做什么技术贴,就是记录一下自己在其中遇到的问题,和自己新学到的东西. ...

  7. [HNOI2004]打鼹鼠

    鼹鼠是一种很喜欢挖洞的动物,但每过一定的时间,它还是喜欢把头探出到地面上来透透气的.根据这个特点阿牛编写了一个打鼹鼠的游戏:在一个\(n*n\)的网格中,在某些时刻鼹鼠会在某一个网格探出头来透透气.你 ...

  8. 论文翻译_Tracking The Untrackable_Learning To Track Multiple Cues with Long-Term Dependencies_IEEE2017

    Tracking The Untrackable: Learning to Track Multiple Cues with Long-Term Dependencies 跟踪不可跟踪:学习跟踪具有长 ...

  9. Mac上利用Aria2加速百度网盘下载

    百度网盘下载东西的速度那叫一个慢,特别是大文件,看着所需时间几个小时以上,让人很不舒服,本文记录自己在mac上利用工具Aria2加速的教程,windows下思路也是一样! 科普(可以不看) 这里顺带科 ...

  10. Flink之状态之savepoint

    1.总览 savepoints是外部存储的自包含的checkpoints,可以用来stop and resume,或者程序升级.savepoints利用checkpointing机制来创建流式作业的状 ...