「BZOJ2823」[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
题解
这不是On的随机增量法么TAT,为啥都在传闻什么模拟退火还有凸包上暴力
这个算法的证明应该很容易百度到
期望是On的
主要是三点求圆心神烦。。
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 68 69 70 71 72 73 74 75 76 77 78 79 80 |
#include<iostream> #include<cstdio> #include<cstring> #include<cstdlib> #include<set> #include<ctime> #include<vector> #include<queue> #include<algorithm> #include<map> #include<cmath> #define eps 1e-6 #define inf 1000000000 #define pa pair<int,int> #define ll long long using namespace std; int n; double r; struct P{ double x,y; }p[1000005],o; //ax^2-bx^2-(2ax-2bx)x+ay^2-by^2-(2ay-2by)y=0 //ax^2-cx^2-(2ax-2cx)x+ay^2-cy^2-(2ay-2cy)y=0 //==> //ax^2-bx^2+ay^2-by^2=(2ax-2bx)x+(2ay-2by)y //ax^2-cx^2+ay^2-cy^2=(2ax-2cx)x+(2ay-2cy)y //==> //l1=r1x+r2y //l2=r3x+r4y double sqr(double x) { return x*x; } double dis(P a,P b) { return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y)); } P getcen(P a,P b,P c) { double x,y; double ax=sqr(a.x),ay=sqr(a.y),bx=sqr(b.x),by=sqr(b.y),cx=sqr(c.x),cy=sqr(c.y); double l1=ax-bx+ay-by,l2=ax-cx+ay-cy; double r1=2*a.x-2*b.x,r2=2*a.y-2*b.y,r3=2*a.x-2*c.x,r4=2*a.y-2*c.y; x=(l1-r2/r4*l2)/(r1-r2/r4*r3); y=(l1-r1/r3*l2)/(r2-r1/r3*r4); if(r4==0)x=(a.x+c.x)/2; if(r3==0)y=(a.y+c.y)/2; return (P){x,y}; } P getcen(P a,P b) { return (P){(a.x+b.x)/2,(a.y+b.y)/2}; } int main() { scanf("%d\n",&n); for(int i=1;i<=n;i++) { scanf("%lf%lf",&p[i].x,&p[i].y); swap(p[i],p[rand()%i+1]); } o=p[1]; for(int i=1;i<=n;i++) { if(dis(o,p[i])<=r+eps)continue; o=getcen(p[1],p[i]);r=dis(p[i],o); for(int j=1;j<i;j++) { if(dis(o,p[j])<=r+eps)continue; o=getcen(p[i],p[j]);r=dis(p[i],o); for(int k=1;k<j;k++) { if(dis(o,p[k])<=r+eps)continue; o=getcen(p[i],p[j],p[k]);r=dis(p[i],o); } } } printf("%.2lf %.2lf %.2lf",o.x,o.y,r); return 0; } |