题意:
在一个二维平面上,给出至多1e5的点,让你构造出一个点和一个距离,使得这些点中有一半(向上取整)以上的点到你构造的点的距离等于你给出的距离。
思路:
网上看到有人用随机化算法过,想了想,挺有道理。
构成一个圆只需要三个不共线的点就可以了,题目中保证了存在这样一个圆,而且有一半以上的点在这个圆上,假设现在有10000个点,有5000个点在一个圆上,那么一次取出三个点在圆上的概率为:C(5000,3)/C(10000,3),大致是1/8,那么一次取出不符合题意的概率是7/8,循环跑上100次的时候脸黑的概率已经降低到1e-6了,肯定能找到。
(虽说随机数并不随机,但在1e5以内还是可以保证随机性的)
那么做法就是任取三点找所构成的圆的圆心,确定半径,遍历点。
PS:求出圆心后一定要判断大小,大于1e9重新随机。
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e5+10;
const double eps=1e-6;
struct Point
{
double x,y;
Point(double _x,double _y)
{
x=_x;y=_y;
}
Point(){}
}a[N];
Point waixin(Point a,Point b,Point c)
{
double a1=b.x-a.x,b1=b.y-a.y,c1=(a1*a1+b1*b1)/2;
double a2=c.x-a.x,b2=c.y-a.y,c2=(a2*a2+b2*b2)/2;
double d=a1*b2-a2*b1;
return Point(a.x+(c1*b2-c2*b1)/d,a.y+(a1*c2-a2*c1)/d);
}
double dis(Point a,Point b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
int n;
void solve()
{
int fin=(n+1)/2;
while(1)
{
int b=rand()%n;
int c=rand()%n;
int d=rand()%n;
if(c==b&&d==c)continue;
Point xin=waixin(a[b],a[c],a[d]);
double distance=dis(xin,a[b]);
if(fabs(xin.x)>1e9||fabs(xin.y)>1e9||distance>1e9)continue;
int k=0;
for(int i=0;i<n;i++)
{
if(fabs(dis(a[i],xin)-distance)<eps)
k++;
if(k==fin)
{
printf("%.10lf %.10lf %.10lf\n",xin.x,xin.y,distance);
return;
}
}
}
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(int i=0;i<n;i++)
scanf("%lf%lf",&a[i].x,&a[i].y);
if(n>=5)
solve();
else
{
if(n==1)
printf("%.10lf %.10lf %.10lf\n",a[0].x+1,a[0].y,1.0);
else
printf("%.10lf %.10lf %.10lf\n",(a[0].x+a[1].x)/2,(a[0].y+a[1].y)/2,dis(a[0],a[1])/2);
}
}
return 0;
}