BZOJ 2823: [AHOI2012]信号塔
Description
在野外训练中,为了确保每位参加集训的成员安全,实时的掌握和收集周边环境和队员信息非常重要,集训队采用的方式是在训练所在地散布N个小型传感器来收集并传递信息,这些传感器只与设在集训地中的信号塔进行通信,信号塔接收信号的覆盖范围是圆形,可以接收到所有分布在该集训区域内所有N个小型传感器(包括在该圆形的边上)发出的信号。信号塔的功率与信号塔接收范围半径的大小成正比,因为是野外训练,只能使用事先储备好的蓄电设备,因此在可以收集所有传感器信息的基础上,还应使得信号塔的功率最小。小龙帮助教官确定了一种信号塔设置的方案,既可以收集到所有N个传感器的信号,又可以保证这个信号塔的功率是最小的。同学们,你们知道,这个信号塔的信号收集半径有多大,它应该设置在何处吗?
Input
共N+1行,第一行为正整数N(1≤N≤1000000),表示队员个数。接下来
N行,每行两个实数用空格分开,分别是第i个队员的坐标X
Output
一行,共三个实数(中间用空格隔开),分别是信号塔的坐标,和信号塔 覆盖的半径。 (注:队员是否在边界上的判断应符合他到圆心的距离与信号塔接收半径之差的绝对值小于10^-6
Sample Input
5
1.200 1.200
2.400 2.400
3.800 4.500
2.500 3.100
3.900 1.300
Sample Output
2.50 2.85 2.10
HINT
1≤N≤500000
题解
三倍经验 => 1336 1337
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 |
#include<cstdio> #include<cmath> #include<iostream> using namespace std; const int Maxn=1000010; const double EPS=1e-8; struct P{ double x,y; P(){} P(double x0,double y0){ x=x0; y=y0; } }; inline double sqr(double x){return x*x;} inline double getDist(P a,P b){return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));} inline double zabs(double x){return (x>=0?x:-x);} inline void swap(P &a,P &b){P tmp=a;a=b;b=tmp;} inline P getO(P a,P b,P c){ double x1=2*(a.x-b.x),x2=2*(a.x-c.x); double y1=2*(a.y-b.y),y2=2*(a.y-c.y); double z1=sqr(a.x)-sqr(b.x)+sqr(a.y)-sqr(b.y), z2=sqr(a.x)-sqr(c.x)+sqr(a.y)-sqr(c.y); P ret; if (zabs(x1)<=EPS){ ret.y=z1/y1; ret.x=(z2-y2*ret.y)/x2; return ret; } if (zabs(y1)<=EPS){ ret.x=z1/x1; ret.y=(z2-x2*ret.x)/y2; return ret; } ret.y=(z1*x2-z2*x1)/(y1*x2-y2*x1); ret.x=(z1*y2-z2*y1)/(x1*y2-x2*y1); return ret; } int n; P p[Maxn]; P O; double R; P A,B,C; int main(){ scanf("%d",&n); for (int i=1;i<=n;i++) scanf("%lf%lf",&p[i].x,&p[i].y); O=p[1];R=0; for (int i=2;i<=n;i++){ if (zabs(getDist(O,p[i])-R)<=EPS || getDist(O,p[i])<=R) continue; O=A=p[i];R=0; for (int j=1;j<i;j++){ if (zabs(getDist(O,p[j])-R)<=EPS || getDist(O,p[j])<=R) continue; B=p[j]; O=P((A.x+B.x)/2,(A.y+B.y)/2); R=getDist(A,B)/2; for (int k=1;k<j;k++){ if (zabs(getDist(O,p[k])-R)<=EPS || getDist(O,p[k])<=R) continue; C=p[k]; O=getO(A,B,C); R=getDist(O,A); } } } printf("%.2lf %.2lf %.2lf\n",O.x,O.y,R); return 0; } |

原文链接:BZOJ 2823: [AHOI2012]信号塔
WNJXYKの博客 版权所有,转载请注明出处。
还没有任何评论,你来说两句吧!