CCPC.2017哈尔滨站-重现赛-M(随机:三点确定一个圆)
2017-11-11 22:24
309 查看
题意:给你n个点,让你找一个圆心和一个半径,使得已知的n个点中至少有(n+1)/2个点在圆上,
输出任意一种合法情况。
题解:对于圆,我们知道三点可以确定一个圆:http://blog.csdn.net/liyuanbhu/article/details/52891868
但是我们不能n^3来找,不然会炸,我们可以采取随机的方法(一丝扣死咪?),第一次用随机操作,比赛时
完全没想到这样搞(可能因此错失了银牌?),随机出三个点后,按照正常的判断即可。我们可以随机多次使得
正确性尽可能的大。。。
#include<math.h>
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
using namespace std;
#define maxn 100005
#define eps 1e-6
struct node1
{
double x,y,r;
}a[maxn];
int n;
double dis(node1 a,node1 b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
bool eq(double a,double b)
{
if(fabs(a-b)<=eps)
return 1;
return 0;
}
bool gongxian(node1 a,node1 b,node1 c)//判断三点是否共线
{
if(eq((b.x-a.x)*(c.y-a.y),(c.x-a.x)*(b.y-a.y)))
return 1;
return 0;
}
node1 yuanxin(node1 a,node1 b,node1 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;
node1 res;
res.x=a.x+(c1*b2-c2*b1)/d;
res.y=a.y+(a1*c2-a2*c1)/d;
return res;
}
bool judge(node1 c)
{
int num=0;
for(int i=0;i<n;i++)
if(eq(c.r,dis(c,a[i])))
num++;
if(num>=(n+1)/2)
return 1;
return 0;
}
int main(void)
{
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==1)
printf("%lf %lf %lf\n",a[0].x+1,a[0].y,1.0);
else if(n<=4)
{
node1 ans;
ans.x=a[0].x+a[1].x;ans.y=a[0].y+a[1].y;
printf("%lf %lf %lf\n",ans.x/2.0,ans.y/2.0,dis(a[0],a[1])/2.0);
}
else
{
while(1)
{
int t1,t2,t3;
t1=rand()%n;
t2=rand()%n;
t3=rand()%n;
if(gongxian(a[t1],a[t2],a[t3]))
continue;
node1 c=yuanxin(a[t1],a[t2],a[t3]);
c.r=dis(c,a[t1]);
if(judge(c))
{
printf("%lf %lf %lf\n",c.x,c.y,c.r);
break;
}
}
}
}
return 0;
}
相关文章推荐
- 由三点确定一个圆(圆心,半径)
- CCPC.2017 哈尔滨站-重现赛-H(暴力+贪心)
- CCPC.2017哈尔滨站-重现赛-D(瞎搞)
- 判断一个点是否在圆内(三点确定一个圆)
- CCPC.2017哈尔滨站-重现赛-B(二分)
- CCPC.2017哈尔滨站-重现赛-A(manacher+树状数组)
- 圆弧方向判断方法和三点确定一个圆的计算方法
- vtk实现三点确定一个平面
- 三点确定一个圆的计算方法
- 三点确定一个圆的计算方法
- 用递归程序求解随机产生一个正整数n(n>=100000),确定n是否是它所有因子之和
- 温故知新--数组(产生一个int数组,长度为100,并向其中随机插入1-100,并且不能重复。)
- 如何随机产生一个自定义长度的字符串?
- 练习2-1 编写一个程序一确定分别由signed及unsigned限定的char,short,int及long类型变量的取值范围。采用打印标准头文件中的相应值以及直接计算两种方式实现。通过直接计算来确定浮点类型的取值范围是一项难度很大的任务。
- 我如何确定一个机器的字节顺序是大端还是小端?
- Sql随机生成一个汉语名字
- python学习——Python中生成一个指定长度的随机字符串
- poj 1228 Grandpa's Estate 给定了一个凸包的部分顶点和边上的点,判断是否能唯一确定一个凸包
- 生成一个随机数据
- hdu2049 不容易系列之(4)——考新郎 错排+组合 一共有N对新婚夫妇,N个新娘随机坐成一排,每个新郎只能选一个, 其中有M个新郎找错了新娘,求发生这种情况一共有多少种可能.