[日常摸鱼]bzoj2823 [AHOI2012]信号塔
题意:$n$个点,求最小圆覆盖,$n \leq 5e5$
这题数据是随机的hhh
我们可以先求出凸包然后对凸包上的点求最小圆覆盖…(不过直接求应该也行?)
反正随便写好像都能过…
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cmath>
using namespace std;
const int N=500005;
struct Point
{
double x,y;
int rnd;
Point(double x=0,double y=0):x(x),y(y){}
}p[N],s[N];
struct Line
{
double k,b;
};
inline bool cmp1(Point a,Point b)
{
if(a.x==b.x)return a.y<b.y;
return a.x<b.x;
}
inline bool cmp2(Point a,Point b)
{
return a.rnd<b.rnd;
}
inline Point operator +(Point a,Point b)
{
return Point(a.x+b.x,a.y+b.y);
}
inline Point operator -(Point a,Point b)
{
return Point(a.x-b.x,a.y-b.y);
} inline Point operator /(Point a,double d)
{
return Point(a.x/d,a.y/d);
}
inline double cross(Point a,Point b)
{
return a.x*b.y-a.y*b.x;
}
inline double sqr2(double x)
{
return x*x;
}
inline double dist(Point a,Point b)
{
return sqrt(sqr2(a.x-b.x)+sqr2(a.y-b.y));
}
inline Line getLine(double k,Point a)
{
Line res;res.k=k;
res.b=a.y-a.x*k;
return res;
}
inline Point getLineIntersection(Line l1,Line l2)
{
Point res;
res.x=(l2.b-l1.b)/(l1.k-l2.k);
res.y=res.x*l1.k+l1.b;
return res;
}
inline Point getCircle(Point a,Point b,Point c)
{
Point p1=(a+b)/2,p2=(a+c)/2;
double k1=-(b.x-a.x)/(b.y-a.y);
double k2=-(c.x-a.x)/(c.y-a.y);
Line l1=getLine(k1,p1),l2=getLine(k2,p2);
return getLineIntersection(l1,l2);
}
inline Point minCircle(double &r,int n)
{
for(register int i=1;i<=n;i++)s[i].rnd=rand();
sort(s+1,s+n+1,cmp2);
Point o=s[1];r=0;
for(register int i=2;i<=n;i++)if(r<dist(o,s[i]))
{
o=s[i];r=0;
for(register int j=1;j<i;j++)if(r<dist(o,s[j]))
{
o=(s[i]+s[j])/2;
r=dist(o,s[i]);
for(register int k=1;k<j;k++)if(r<dist(o,s[k]))
{
o=getCircle(s[i],s[j],s[k]);
r=dist(o,s[i]);
}
}
}
return o;
}
inline int convexHull(int n)
{
sort(p+1,p+n+1,cmp1);
int t=0,k;
for(register int i=1;i<=n;i++)
{
while(t>1&&cross(s[t]-s[t-1],p[i]-s[t-1])<0)t--;
s[++t]=p[i];
}
k=t;
for(register int i=n-1;i>=1;i--)
{
while(t>k&&cross(s[t]-s[t-1],p[i]-s[t-1])<0)t--;
s[++t]=p[i];
}
if(n>1)t--;
return t;
}
int main()
{
int n;scanf("%d",&n);
for(register int i=1;i<=n;i++)scanf("%lf%lf",&p[i].x,&p[i].y);
int t=convexHull(n);double r;
Point res=minCircle(r,t);
printf("%.2lf %.2lf %.2lf",res.x,res.y,r);
return 0;
}
[日常摸鱼]bzoj2823 [AHOI2012]信号塔的更多相关文章
- 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]信号塔&&1336: [Balkan2002]Alien最小圆覆盖&&1337: 最小圆覆盖
首先我写了个凸包就溜了 这是最小圆覆盖问题,今晚学了一下 先随机化点,一个个加入 假设当前圆心为o,半径为r,加入的点为i 若i不在圆里面,令圆心为i,半径为0 再重新从1~i-1不停找j不在圆里面, ...
- 【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】2823: [AHOI2012]信号塔
题意 给\(n\)个点,求一个能覆盖所有点的面积最小的圆.(\(n \le 50000\)) 分析 随机增量法 题解 理论上\(O(n^3)\)暴力,实际上加上随机化后期望是\(O(n)\)的. 算法 ...
- BZOJ 2823: [AHOI2012]信号塔
题目:http://www.lydsy.com/JudgeOnline/problem.php?id=2823 随机增量法.不断加点维护圆,主要是三点共圆那里打得烦(其实也就是个两中垂线求交点+联立方 ...
随机推荐
- web安全入门--入门条件
介绍:网络安全由习大大提出,是继海.陆.空.外太空的第五大作战领域,也是一个关系国家安全和主权.社会稳定.民族文化继承和发扬的重要问题.其重要性,正随着全球信息化步伐的加快越来越重要.接下来,我向大家 ...
- list scheduling algorithm 指令调度 —— 笔记
作者:Yaong 出处:https://www.cnblogs.com/yaongtime/articles/14033444.html 版权:本文版权归作者和博客园共有 转载:欢迎转载,但未经作者同 ...
- 通过Camtasia来制作画中画视频效果的方法
随着全民娱乐化的发展,视频的形式也更加多种多样了.视频形式的多样化能让观众从不同形式的视频中观赏到更有趣味的内容.比如像画中画的视频形式,让视频中的人物看起来像与观众一同观看视频,或者形成两个视频的对 ...
- leetcode 56合并区间 java
//先排序,将左区间小的放在前面,然后如果前一个的右区间大于下一个的左区间,则可以合并,分别用两个下标指向当前的大区间和将要考察的小区间 class Solution { public int[ ...
- 学习django笔记一:在urls.py中导入sign应用views文件的问题
>python-admin startproject guest #创建guest项目 >python3 manage.py startapp sign #在guest项目中创建 ...
- docker提示容器已存在
docker ps -a docker rm 容器id 重启启动
- (原创)用.NET Core实现一个在线客服系统(上篇)
前言 没有视频的介绍显得尤为空白仓促.所以,如果你不赶时间,看看视频先 → → 戳我看视频 ← ← 在线演示访客端:http://role.fuyue.xyz/visitor/index客服端:ht ...
- Linux 学习笔记04丨Linux的用户和用户组管理
Chapter 3. 用户和用户组管理 由于Linux系统是一个多用户多任务的分时操作系统,任何一个要使用系统资源的用户,都必须首先向系统管理员申请一个账号,然后以该账号身份进入系统. 3.0 用户与 ...
- C语言讲义——数组和指针
数组名表示的是这个数组的首地址.即如果有int a[10],则a 相当于&a[0]. #include <stdio.h> main() { int a[5]= {1,3,5,7, ...
- Maven 依赖树的解析规则
对于 Java 开发工程师来说,Maven 是依赖管理和代码构建的标准.遵循「约定大于配置」理念.Maven 是 Java 开发工程师日常使用的工具,本篇文章简要介绍一下 Maven 的依赖树解析. ...