hihoCoder1227 The Cats' Feeding Spots 二分

题目链接：http://hihocoder.com/problemset/problem/1227

题目大意：平面上有m个点，要从这m个点当中找出n个点，使得包含这n个点的圆的半径（圆心为n个点当中的某一点且半径为整数）最小，同时保证圆周上没有点。

分析：枚举这m个点，以每个点为圆心找出一个符合条件的半径，然后找出这些半径中的最小值即可。至于如何找出半径，我们可以二分[-1001,1001]区间。

实现代码如下：

#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
#define INF 0x7fffffff
struct POINT
{
double x,y;
bool operator <(const POINT &a) const
{
if(x==a.x) return y<a.y;
return x<a.x;
}
}p[110];
double d[110][110];
int m,n,minr;
bool flag;
int min(int a,int b)
{
return a<b?a:b;
}
double dist(POINT a,POINT b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
void init()
{
for(int i=0;i<m;i++)
for(int j=i;j<m;j++)
d[i][j]=d[j][i]=dist(p[i],p[j]);
}
int judge(int cnt,int rad)
{
int num=0;
flag=true;
for(int i=0;i<m;i++)
{
if(d[cnt][i]<rad*1.0) num++;
if(d[cnt][i]==rad*1.0) flag=false;
}
return num;
}
void solve()
{
init();
minr=INF;
for(int i=0;i<m;i++)
{
int l=-1001,r=1001;
while(l<=r)
{
int m=(l+r)>>1;
int tmp=judge(i,m);
if(tmp==n&&flag)  minr=min(minr,m);
if(tmp>=n) r=m-1;
else l=m+1;
}
}
if(minr==INF) puts("-1");
else printf("%d\n",minr);
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&m,&n);
for(int i=0;i<m;i++)
scanf("%lf%lf",&p[i].x,&p[i].y);
sort(p,p+m);
solve();
}
return 0;
}