bzoj2823: [AHOI2012]信号塔&&1336: [Balkan2002]Alien最小圆覆盖&&1337: 最小圆覆盖
首先我写了个凸包就溜了
这是最小圆覆盖问题,今晚学了一下
先随机化点,一个个加入
假设当前圆心为o,半径为r,加入的点为i
若i不在圆里面,令圆心为i,半径为0
再重新从1~i-1不停找j不在圆里面,令圆心为ij中点,直径为ij距离
再重新在1~j-1不停找k不在圆里面,三点可确定一圆,初中数学
复杂度看似O(n^3)实则O(n),好玄学
坑点:注意如果用点斜式表示方程有斜率为不存在的情况,需要特判
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
const double eps=1e-;
double sqr(double x){return x*x;} struct point{ double x,y;point(){} point(double X,double Y){x=X,y=Y;} };
double getdis(point p1,point p2){return sqrt(sqr(p1.x-p2.x)+sqr(p1.y-p2.y));}
point middle(point p1,point p2){return point((p1.x+p2.x)/,(p1.y+p2.y)/);}
double slope (point p1,point p2)
{
if(p2.x==p1.x)return ;
return (p2.y-p1.y)/(p2.x-p1.x);
}
double multi(point p1,point p2,point p0)
{
double x1,y1,x2,y2;
x1=p1.x-p0.x;
y1=p1.y-p0.y;
x2=p2.x-p0.x;
y2=p2.y-p0.y;
return x1*y2-x2*y1;
} struct segment{ double k,b;segment(){} segment(double K,double B){k=K,b=B;} };
segment getseg(double k,point pp){return segment(k,pp.y-k*pp.x);}
point intersection(segment s1,segment s2)
{
double x=(s2.b-s1.b)/(s1.k-s2.k);
double y=s1.k*x+s1.b;
return point(x,y);
} //--------------------------------------simple-------------------------------------------------------- int n; point p[];
bool cmp(point p1,point p2)
{
double d=multi(p1,p2,p[]);
if(fabs(d)<=eps)return getdis(p1,p[])<getdis(p2,p[]);
else return d>;
}
int top,sta[];
void graham()
{
sort(p+,p+n+,cmp);
top=;sta[++top]=,sta[++top]=;
double g;
for(int i=;i<=n;i++)
{
while(top>=)
{
g=multi(p[sta[top]],p[i],p[sta[top-]]);
if(g<||fabs(g)<=eps)top--;
else break;
}
sta[++top]=i;
}
} //------------------------------------graham---------------------------------------------------------- point getcore(point p1,point p2,point p3)
{
double g=multi(p1,p2,p3);
if(fabs(g)<=eps)
{
double d1=getdis(p1,p2),d2=getdis(p1,p3),d3=getdis(p2,p3);
if(d1>d2&&d1>d3)return middle(p1,p2);
if(d2>d1&&d2>d3)return middle(p1,p3);
if(d3>d1&&d3>d2)return middle(p2,p3);
}
else
{
segment s1,s2;
if(slope(p1,p2)==)
{
s1=getseg(-/slope(p1,p3),middle(p1,p3));
s2=getseg(-/slope(p2,p3),middle(p2,p3));
}
else if(slope(p1,p3)==)
{
s1=getseg(-/slope(p1,p2),middle(p1,p2));
s2=getseg(-/slope(p2,p3),middle(p2,p3));
}
else
{
s1=getseg(-/slope(p1,p2),middle(p1,p2));
s2=getseg(-/slope(p1,p3),middle(p1,p3));
}
return intersection(s1,s2);
}
} void circlecover()
{
random_shuffle(sta+,sta+top+);
point o=p[sta[]];double r=,d;
for(int i=;i<=top;i++)
if(getdis(o,p[sta[i]])>r)
{
o=p[sta[i]],r=;
for(int j=;j<i;j++)
if(getdis(o,p[sta[j]])>r)
{
o=middle(p[sta[i]],p[sta[j]]),r=getdis(o,p[sta[i]]);
for(int k=;k<j;k++)
if(getdis(o,p[sta[k]])>r)
o=getcore(p[sta[i]],p[sta[j]],p[sta[k]]),r=getdis(o,p[sta[i]]);
}
}
printf("%.2lf %.2lf %.2lf\n",o.x,o.y,r);
} //------------------------------------solve---------------------------------------------------------- int main()
{
scanf("%d",&n);
for(int i=;i<=n;i++)
{
scanf("%lf%lf",&p[i].x,&p[i].y);
if(p[i].y<p[].y||(p[i].y==p[].y&&p[i].x<p[].x))
swap(p[i],p[]);
}
graham();
circlecover(); return ;
}
bzoj2823: [AHOI2012]信号塔&&1336: [Balkan2002]Alien最小圆覆盖&&1337: 最小圆覆盖的更多相关文章
- bzoj2823[AHOI2012]信号塔
2823: [AHOI2012]信号塔 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1190 Solved: 545[Submit][Status ...
- BZOJ2823 [AHOI2012]信号塔 【最小圆覆盖】
题目链接 BZOJ2823 题解 最小圆覆盖模板 都懒得再写一次 #include<iostream> #include<cstdio> #include<cmath&g ...
- [日常摸鱼]bzoj2823 [AHOI2012]信号塔
题意:$n$个点,求最小圆覆盖,$n \leq 5e5$ 这题数据是随机的hhh 我们可以先求出凸包然后对凸包上的点求最小圆覆盖-(不过直接求应该也行?) 反正随便写好像都能过- #include&l ...
- 【BZOJ2823】[AHOI2012]信号塔(最小圆覆盖)
[BZOJ2823][AHOI2012]信号塔(最小圆覆盖) 题面 BZOJ 洛谷 相同的题: BZOJ1 BZOJ2 洛谷 题解 模板题... #include<iostream> #i ...
- 2018.07.04 BZOJ 2823: AHOI2012信号塔(最小圆覆盖)
2823: [AHOI2012]信号塔 Time Limit: 10 Sec Memory Limit: 128 MB Description 在野外训练中,为了确保每位参加集训的成员安全,实时的掌握 ...
- 【bzoj2823】 AHOI2012—信号塔
http://www.lydsy.com/JudgeOnline/problem.php?id=2823 (题目链接) 题意 求最小圆覆盖 Solution 关于最小圆覆盖的做法,论文里面都有.其实真 ...
- (bzoj1337 || 洛谷P1742 最小圆覆盖 )|| (bzoj2823 || 洛谷P2533 [AHOI2012]信号塔)
bzoj1337 洛谷P1742 用随机增量法.讲解:https://blog.csdn.net/jokerwyt/article/details/79221345 设点集A的最小覆盖圆为g(A) 可 ...
- [BZOJ 1336] [Balkan2002] Alien最小圆覆盖 【随机增量法】
题目链接:BZOJ - 1336 题目分析 最小圆覆盖有一个算法叫做随机增量法,看起来复杂度像是 O(n^3) ,但是可以证明其实平均是 O(n) 的,至于为什么我不知道= = 为什么是随机呢?因为算 ...
- 【BZOJ】2823: [AHOI2012]信号塔
题意 给\(n\)个点,求一个能覆盖所有点的面积最小的圆.(\(n \le 50000\)) 分析 随机增量法 题解 理论上\(O(n^3)\)暴力,实际上加上随机化后期望是\(O(n)\)的. 算法 ...
随机推荐
- ES6(对象扩展)
ES6(对象(object)新增特性) 1.简介表示法 o,k 为属性:这种为无方法的情况. 这种为有方法. 2.属性表达式 ES6中,[a]就是 b . 3.新增API 1.数组也是引用类型 2.O ...
- Using TCP keepalive under Linux
Linux has built-in support for keepalive. You need to enable TCP/IP networking in order to use it. Y ...
- HDU 3911 区间合并求最大长度的问题
http://vjudge.net/problem/viewProblem.action?id=21557 题目大意: 每进行一次颜色改变都可以把一段区间内的黑石头变成白石头,白石头变成黑石头,最后问 ...
- hdu 4639
/*一个he都没有和有一个he结果是1 2个he 是2 3个he 3 4个he 5 5个he 8 ..... 把每段的结果相乘取余*/ #include<stdio.h> ...
- CodeForces - 43B Letter
字符串的处理 统计已有字符的个数 和需求字符比较 #include <iostream> #include <stdio.h> #include <string.h> ...
- LCA 在线倍增法 求最近公共祖先
第一步:建树 这个就不说了 第二部:分为两步 分别是深度预处理和祖先DP预处理 DP预处理: int i,j; ;(<<j)<n;j++) ;i<n;++i) ) fa[i ...
- HDU3549 最大流 裸题
EK算法 时间复杂度o(n*m*m) 因为有反向边每次bfs时间为 n*m 每次删一条边 最多m次 代码 #include<iostream> #include<string.h& ...
- Mysql 数据库允许远程连接 服务器连接错误 Host 'XXX' is not allowed to connect to this MySQL server
如果连接数据库的时候出现这个问题 Host 'XXX' is not allowed to connect to this MySQL server 说明 Mysql数据库 不允许远程连接, 需要修改 ...
- Entity framework自定义字段实现思路
ublic class MyModel { public int MyModelID { get; set; } public string FixedProperty1 { get; set; } ...
- DELPHI7调用BERLIN中间件的中文字段名乱码的解决办法
MSSQL数据库的表使用中文字段名,BERLIN开发的DATASNAP中间件,DELPHI7调用中间件的查询方法返回数据给CLIENTDATASET.DATA,发现中文字段名乱码,中文字段名的值可以正 ...