题面

传送门

题解

之前只是在抄题解……这篇才算是真正自己想的吧……

首先我们把输入序列给\(random\)一下防止出题人好心送你一个毒瘤序列

我们设\(r\)为当前最大半径,\(o\)为此时对应圆心

先说一下算法过程:

令前\(i-1\)个点的最小覆盖圆为\((o,r)\)

如果第\(i\)个点在这个圆中,直接跳过

如果不在,那么第\(i\)个点一定在前\(i\)个点的最小覆盖圆上,此时前\(i-1\)个点中还有两个也在最小覆盖圆上。那么我们设\(o=p_i,r=0\),做固定了第\(i\)个点的最小圆覆盖

固定了一个点:不停地在范围内找第一个不在圆内的点\(p_j\),令\(o=(p_i+p_j)/2,r={1\over 2}|p_ip_j|\),做固定了两个点的最小圆覆盖

固定了两个点:在范围内找第一个不在圆内的点\(p_k\),将圆设为\(p_i,p_j,p_k\)的外接圆

这一个\(O(n^3)\)的算法有个啥用?我还不如写个暴力呢

然而实际上上面这个算法在给定点的序列是随机的情况下是期望\(O(n)\)的(这也是我们为什么一开始要\(random\)的原因)

因为最小圆覆盖只会由三个点确定,所以如果有\(n\)个点,那么这其中每一个点参与最小圆覆盖的概率为\({3\over n}\)

也就是说,对于第\(i\)个点,它不在前\(i-1\)个点的外接圆中的概率是\({3\over i}\),也就是说会从第一层进入第二层的概率是\({3\over i}\)

类似的,第二层进入第三层的概率也是\({3\over i}\)

所以时间复杂度为

\[T_1(n)=O(n)+\sum_{i=1}^n{3\over i}T_2(i)
\]

\[T_2(n)=O(n)+\sum_{i=1}^n{3\over i}T_3(i)
\]

\[T_3(n)=O(n)
\]

综上,有\(T_1(n)=O(n)\)

然后是一些细节问题,比方说怎么求三点共圆

只要用任意两条直线的中垂线的交点找出圆心就可以了

以及我们万一有三点共线导致我们求的两条中垂线平行怎么办呀

这个的话可以通过加上一些微小的扰动值来避免(虽然好像不加也能\(A\)就是了)

//minamoto
#include<bits/stdc++.h>
#define R register
#define inline __inline__ __attribute__((always_inline))
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
char buf[1<<21],*p1=buf,*p2=buf;
inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
int read(){
R int res,f=1;R char ch;
while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
return res*f;
}
double readdb()
{
R double x=0,y=0.1,f=1;R char ch;
while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
for(x=ch-'0';(ch=getc())>='0'&&ch<='9';x=x*10+ch-'0');
for(ch=='.'&&(ch=getc());ch>='0'&&ch<='9';x+=(ch-'0')*y,y*=0.1,ch=getc());
return x*f;
}
const int N=1e5+5;
struct node{
double x,y;
inline node(){}
inline node(R double xx,R double yy):x(xx),y(yy){}
inline node operator +(const node &b)const{return node(x+b.x,y+b.y);}
inline node operator -(const node &b)const{return node(x-b.x,y-b.y);}
inline double operator *(const node &b)const{return x*b.y-y*b.x;}
inline node operator *(const double &b)const{return node(x*b,y*b);}
inline double operator ^(const node &b)const{return x*b.x+y*b.y;}
inline double len2(){return x*x+y*y;}
inline node rot(){return node(-y,x);}
}p[N],o;
struct Line{
node p,v;
inline Line(){}
inline Line(R node pp,R node vv):p(pp),v(vv){}
friend node cross(const Line &a,const Line &b){return a.p+a.v*(b.v*(b.p-a.p)/(b.v*a.v));}
};
node circle(const node &a,const node &b,const node &c){
return cross(Line((a+b)*0.5,(b-a).rot()),Line((a+c)*0.5,(c-a).rot()));
}
int n;double r;
int main(){
// freopen("testdata.in","r",stdin);
n=read();
fp(i,1,n)p[i].x=readdb(),p[i].y=readdb();
random_shuffle(p+1,p+1+n);
r=0,o=node(0,0);
fp(i,1,n)if((p[i]-o).len2()>r){
o=p[i],r=0;
fp(j,1,i-1)if((p[j]-o).len2()>r){
o=(p[i]+p[j])*0.5,r=(p[j]-o).len2();
fp(k,1,j-1)if((p[k]-o).len2()>r)
o=circle(p[i],p[j],p[k]),r=(p[k]-o).len2();
}
}
printf("%.10lf\n%.10lf %.10lf\n",sqrt(r),o.x,o.y);
return 0;
}

洛谷P1742 最小圆覆盖(计算几何)的更多相关文章

  1. 洛谷P1742 最小圆覆盖(计算几何)

    题意 题目链接 Sol 暴力做法是\(O(n^3)\)枚举三个点然后check一下是否能包含所有点 考虑一种随机算法,首先把序列random_shuffle一下. 然后我们枚举一个点\(i\),并维护 ...

  2. (bzoj1337 || 洛谷P1742 最小圆覆盖 )|| (bzoj2823 || 洛谷P2533 [AHOI2012]信号塔)

    bzoj1337 洛谷P1742 用随机增量法.讲解:https://blog.csdn.net/jokerwyt/article/details/79221345 设点集A的最小覆盖圆为g(A) 可 ...

  3. 洛谷 P1742 最小圆覆盖 (随机增量)

    题目链接:P1742 最小圆覆盖 题意 给出 N 个点,求最小的包含所有点的圆. 思路 随机增量 最小圆覆盖一般有两种做法:随机增量和模拟退火.随机增量的精确度更高,这里介绍随机增量的做法. 先将所有 ...

  4. 最小圆覆盖(洛谷 P1742 增量法)

    题意:给定N个点,求最小圆覆盖的圆心喝半径.保留10位小数点. N<1e5: 思路:因为精度要求较高,而且N比较大,所以三分套三分的复杂度耶比较高,而且容易出错. 然是写下增量法吧. 伪代码加深 ...

  5. 洛谷CF1071E Rain Protection(计算几何,闵可夫斯基和,凸包,二分答案)

    洛谷题目传送门 CF题目传送门 对于这题,我无力吐槽. 虽然式子还是不难想,做法也随便口胡,但是一些鬼畜边界情况就是判不对. 首先显然二分答案. 对于每一个雨滴,它出现的时刻我们的绳子必须落在它上面. ...

  6. 网络流24题 第三题 - CodeVS1904 洛谷2764 最小路径覆盖问题 有向无环图最小路径覆盖 最大流 二分图匹配 匈牙利算法

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - CodeVS1904 题目传送门 - 洛谷2764 题意概括 给出一个有向无环图,现在请你求一些路径,这些路径 ...

  7. 洛谷P3222 [HNOI2012]射箭(计算几何,半平面交,双端队列)

    洛谷题目传送门 设抛物线方程为\(y=ax^2+bx(a<0,b>0)\),我们想要求出一组\(a,b\)使得它尽可能满足更多的要求.这个显然可以二分答案. 如何check当前的\(mid ...

  8. P1742 最小圆覆盖(计算几何)

    体验过\(O(n^3)\)过\(10^5\)吗?快来体验一波当\(wys\)的快感吧\(QAQ\) 前置芝士1:二元一次方程组求解 设 \[\begin{cases}a1 * x + b1*y=c1\ ...

  9. 洛谷P2764 最小路径覆盖问题

    有向无环图的最小路径点覆盖 最小路径覆盖就是给定一张DAG,要求用尽量少的不相交的简单路径,覆盖有向无环图的所有顶点. 有定理:顶点数-路径数=被覆盖的边数. 要理解的话可以从两个方向: 假设DAG已 ...

随机推荐

  1. java.util.ConcurrentModificationException异常分析

    Java在操作ArrayList.HashMap.TreeMap等容器类时,遇到了java.util.ConcurrentModificationException异常.以ArrayList为例,如下 ...

  2. flume 配置说明

    Flume中的HDFS Sink应该是非常常用的,其中的配置参数也比较多,在这里记录备忘一下. channel type hdfs path 写入hdfs的路径,需要包含文件系统标识,比如:hdfs: ...

  3. thymeleaf 获取项目路径

    <p th:text=${salecode}></p> <a th:href="${#httpServletRequest.getScheme()+'://'+ ...

  4. CloudStack 4.1快照测试

    目前4.1只支持存储快照,4.2能支持内快照 1. 选中实例 2. 查看实例卷 3. 创建快照 4. 通过快照创建模板 5. 查看通过快照创建的模板 6. 通过快照创建模板生成的实例 7. 自动定制创 ...

  5. 如何安装JDeveloper

    1>下载jdk和JDeveloper jdk版本:jdk-7u65-windows-x64.EXE JDeveloper版本:jdev_suite_12.3_win64_V44419-01.is ...

  6. sql添加/移除约束

    唯一:ALTER TABLE TableName ADD CONSTRAINT UQ_NickName UNIQUE(NickName) 主键:ALTER TABLE TableName ADD CO ...

  7. 【git使用】Failed to connect to 127.0.0.1 port 1080: Connection refused

    查询是否使用代理:git config --global http.proxy 取消代理:git config --global --unset http.proxy

  8. 遍历ListView,查出每一项的内容

    private ListView.OnItemClickListener showItemDetail = new ListView.OnItemClickListener() { public vo ...

  9. Maven项目下java.lang.ClassNotFoundException常规解决办法

    网上很多要修改.class .project的比较麻烦有时候还不一定管用.下面的方法适合于项目用已经引用了jar,但是运行时却ClassNotFound,请用如下方法试试: 严重: Error con ...

  10. Hadoop中Writable类之二

    1.ASCII.Unicode.UFT-8 在看Text类型的时候,里面出现了上面三种编码,先看看这三种编码: ASCII是基于拉丁字母的一套电脑编码系统.它主要用于显示现代英语和其他西欧语言.它是现 ...