\(\color{#0066ff}{题目描述}\)

给出N个点,让你画一个最小的包含所有点的圆。

\(\color{#0066ff}{输入格式}\)

先给出点的个数N,2<=N<=100000,再给出坐标Xi,Yi.(-10000.0<=xi,yi<=10000.0)

\(\color{#0066ff}{输出格式}\)

输出圆的半径,及圆心的坐标,保留10位小数

\(\color{#0066ff}{输入样例}\)

6
8.0 9.0
4.0 7.5
1.0 2.0
5.1 8.7
9.0 2.0
4.5 1.0

\(\color{#0066ff}{输出样例}\)

5.0000000000
5.0000000000 5.0000000000

\(\color{#0066ff}{数据范围与提示}\)

none

\(\color{#0066ff}{题解}\)

随机增量法

前置知识,三点定圆

设圆上有三个点\(A(x_1, y_1), B(x_2, y_2), C(x_3, y_3)\)

我们可以找这些点的连线的任两条线,找它们垂直平分线的交点就是圆心

易得向量(线段)AB,BC

设两条线段的垂直平分线为\(A_1x+B_1y=C_1,A_2x+B_2y=C_2\)

上式的A和B分别是两个向量的x和y(法向量(可以带入验证))

之后带入AB,BC的中点来求出两个C

联立两个垂直平分线,解出来

\(x=\frac{C_2*B_1-C_1*B_2}{A_2*B_1-A_1*B_2}, y=\frac{C_2*A_1-C_1*A_2}{B_2*A_1-B_1*A_2}\)

有了圆心,半径自然好求,随便找圆上一点(A,B,C)与圆心求距离就是半径了

随机增量

初始设置第一个点为圆心,半径为0

依次扫每个点,如果不在当前圆内,则以那个点为圆心

再次扫从头开始到当前的每个点,如果不在当前圆中

则当前点与那个点作为新圆直径,构成当前圆,再次扫到当前每个点

如果仍有不在圆内的,三点定圆

#include <bits/stdc++.h>
#define _ 0
#define LL long long
inline LL in() {
LL x = 0, f = 1; char ch;
while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
while(isdigit(ch)) x = x * 10 + (ch ^ 48), ch = getchar();
return x * f;
}
const int maxn = 1e5 + 10;
const double eps = 1e-14;
struct node {
double x, y;
node(double x = 0, double y = 0)
:x(x), y(y) {}
double mod() {
return sqrt(x * x + y * y);
}
friend node operator - (const node &a, const node &b) {
return node(a.x - b.x, a.y - b.y);
}
friend node operator + (const node &a, const node &b) {
return node((a.x + b.x) / 2.0, (a.y + b.y) / 2.0);
}
}e[maxn];
int n; void circle(node a, node b, node c, node &o, double &r) {
node ab = b - a, bc = c - b, mid1 = a + b, mid2 = b + c;
double A1 = ab.x, B1 = ab.y, A2 = bc.x, B2 = bc.y;
double C1 = A1 * mid1.x + B1 * mid1.y, C2 = A2 * mid2.x + B2 * mid2.y;
o = node((C2 * B1 - C1 * B2) / (A2 * B1 - A1 * B2), (C2 * A1 - C1 * A2) / (B2 * A1 - B1 * A2));r = (o - a).mod();
} void work() {
node o = e[1];
double r = 0;
for(int i = 2; i <= n; i++)
if((e[i] - o).mod() - r > eps) {
o = e[i], r = 0;
for(int j = 1; j <= i - 1; j++)
if((e[j] - o).mod() - r > eps) {
o = e[i] + e[j];
r = (o - e[i]).mod();
for(int k = 1; k <= j - 1; k++)
if((e[k] - o).mod() - r > eps)
circle(e[i], e[j], e[k], o, r);
}
}
printf("%.10f\n%.10f %.10f", r, o.x, o.y);
} int main() {
n = in();
for(int i = 1; i <= n; i++) scanf("%lf%lf", &e[i].x, &e[i].y);
work();
return 0;
}

P1742 最小圆覆盖的更多相关文章

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

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

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

    题面 传送门 题解 之前只是在抄题解--这篇才算是真正自己想的吧-- 首先我们把输入序列给\(random\)一下防止出题人好心送你一个毒瘤序列 我们设\(r\)为当前最大半径,\(o\)为此时对应圆 ...

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

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

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

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

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

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

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

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

  7. 【BZOJ-1336&1337】Alie最小圆覆盖 最小圆覆盖(随机增量法)

    1336: [Balkan2002]Alien最小圆覆盖 Time Limit: 1 Sec  Memory Limit: 162 MBSec  Special JudgeSubmit: 1573   ...

  8. Bzoj 1336&1337 Alien最小圆覆盖

    1336: [Balkan2002]Alien最小圆覆盖 Time Limit: 1 Sec  Memory Limit: 162 MBSec  Special Judge Submit: 1473  ...

  9. hdu3007Buried memory(最小圆覆盖)

    链接 普通的暴力复杂度达到O(n^4),对于这题肯定是不行的. 解法:随机增量算法 参考http://www.2cto.com/kf/201208/149602.html algorithm:A.令C ...

随机推荐

  1. Python操作Excel,并结合unittest单元测试框架

    第一步:写Excel操作方法 excel_operate.py文件 from openpyxl import load_workbook #引入模块 class MyExcel: def __init ...

  2. 一段小程序理解getchar和putchar

    #include "stdafx.h" #include <iostream> using namespace std; int main() { char c,d,e ...

  3. docker 安装mysql 使用navicat访问 解决

    1. 下载Mysql的Docker镜像: 2. 运行镜像,设置root账号初始密码(123456),映射本地宿主机端口3306到Docker端口3306.测试过程没有挂载本地数据盘: 3. 查看已运行 ...

  4. windows服务编写和“以管理员运行”程序的方法

    本文将首先解释如何 创建 一个定期查询可用物理内存并将结果写入某个文本文件的服务.然后指导你完成生成,安装和实现服务的整个过程. 第一步:主函数和全局定义 首先,包含所需的头文件.例子要调用 Win3 ...

  5. Eloquent ORM模型中添加自定义值

    我们都知道通过Laravel中数据库查询出来的模型对象都是基于数据库字段,今天给大家展示一个 Laravel Eloquent ORM 模型特性-附加值不存在于数据表中. 举个简单的栗子,一篇文章(p ...

  6. php学习之try catch

    PHP 5 添加了类似于其它语言的异常处理模块.在 PHP 代码中所产生的异常可被 throw语句抛出并被 catch 语句捕获.(注:一定要先抛才能获取) 需要进行异常处理的代码都必须放入 try ...

  7. Solaris与Windows Active Directory集成

    通过Solaris与Active Directory的集成,Solaris可以使用Windows 2003 R2/ 2008 Active Directory来进行用户登录验证.以下是简要配置过程. ...

  8. 部署和调优 2.2 squid反向代理

    配置反向代理 打开配置文件 vim /etc/squid/squid.conf 修改 http_port 改为 http_port 80 accel vhost vport 在它下面添加一段 cach ...

  9. linux命令-du查看占用磁盘空间大小

    格式 df -h 查看磁盘分区情况 du /etc 目录文件大小都列出来 单位是k最后一行是总和 du -m 单位是m 小于1m写成1m du -h 单位人性化显示k/m du -sh /etc 查看 ...

  10. apache2不识别php

    sudo apt-get install libapache2-mod-php7.0 sudo a2enmod php7.0 sudo service apache2 restart 注意:Apach ...