P1742 最小圆覆盖
\(\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 最小圆覆盖的更多相关文章
- 洛谷 P1742 最小圆覆盖 (随机增量)
题目链接:P1742 最小圆覆盖 题意 给出 N 个点,求最小的包含所有点的圆. 思路 随机增量 最小圆覆盖一般有两种做法:随机增量和模拟退火.随机增量的精确度更高,这里介绍随机增量的做法. 先将所有 ...
- 洛谷P1742 最小圆覆盖(计算几何)
题面 传送门 题解 之前只是在抄题解--这篇才算是真正自己想的吧-- 首先我们把输入序列给\(random\)一下防止出题人好心送你一个毒瘤序列 我们设\(r\)为当前最大半径,\(o\)为此时对应圆 ...
- (bzoj1337 || 洛谷P1742 最小圆覆盖 )|| (bzoj2823 || 洛谷P2533 [AHOI2012]信号塔)
bzoj1337 洛谷P1742 用随机增量法.讲解:https://blog.csdn.net/jokerwyt/article/details/79221345 设点集A的最小覆盖圆为g(A) 可 ...
- P1742 最小圆覆盖(计算几何)
体验过\(O(n^3)\)过\(10^5\)吗?快来体验一波当\(wys\)的快感吧\(QAQ\) 前置芝士1:二元一次方程组求解 设 \[\begin{cases}a1 * x + b1*y=c1\ ...
- 洛谷P1742 最小圆覆盖(计算几何)
题意 题目链接 Sol 暴力做法是\(O(n^3)\)枚举三个点然后check一下是否能包含所有点 考虑一种随机算法,首先把序列random_shuffle一下. 然后我们枚举一个点\(i\),并维护 ...
- 最小圆覆盖(洛谷 P1742 增量法)
题意:给定N个点,求最小圆覆盖的圆心喝半径.保留10位小数点. N<1e5: 思路:因为精度要求较高,而且N比较大,所以三分套三分的复杂度耶比较高,而且容易出错. 然是写下增量法吧. 伪代码加深 ...
- 【BZOJ-1336&1337】Alie最小圆覆盖 最小圆覆盖(随机增量法)
1336: [Balkan2002]Alien最小圆覆盖 Time Limit: 1 Sec Memory Limit: 162 MBSec Special JudgeSubmit: 1573 ...
- Bzoj 1336&1337 Alien最小圆覆盖
1336: [Balkan2002]Alien最小圆覆盖 Time Limit: 1 Sec Memory Limit: 162 MBSec Special Judge Submit: 1473 ...
- hdu3007Buried memory(最小圆覆盖)
链接 普通的暴力复杂度达到O(n^4),对于这题肯定是不行的. 解法:随机增量算法 参考http://www.2cto.com/kf/201208/149602.html algorithm:A.令C ...
随机推荐
- I/O完成端口(IOCP)
服务器: #include "stdafx.h" #include <winsock2.h> #pragma comment(lib, "ws2_32.lib ...
- java 多线程系列---JUC原子类(二)之AtomicLong原子类
概要 AtomicInteger, AtomicLong和AtomicBoolean这3个基本类型的原子类的原理和用法相似.本章以AtomicLong对基本类型的原子类进行介绍. AtomicLong ...
- eclipse中maven的配置与使用
以eclipse Juno版本为例 1.插件安装 eclipse==>help====>install new software===>add name :m2e location: ...
- 【265】shell文件创建链接
优点:可以在其他文件夹内运行对应的*.sh文件,例如通过putty会默认进入的文件夹,可以将链接文件放在那里,就可以直接调用了! 方法:很简单 1. 在文件上点击右键>创建链接 2. 可以对下面 ...
- HTML标签详细讲解
http://www.cnblogs.com/yuanchenqi/articles/5603871.html
- find查找、split分隔、replace替换
#!/usr/bin/env python r = "asada" ret = r.find("d") print(ret)#返回所在位置的索引 ret =r. ...
- 36-图像有用区(dfs, bfs)
http://acm.nyist.edu.cn/JudgeOnline/problem.php?pid=92 图像有用区域 时间限制:3000 ms | 内存限制:65535 KB 难度:4 ...
- Linux共享对象之编译参数 -fPIC
转载自:https://www.cnblogs.com/cswuyg/p/3830703.html 在Linux系统中,动态链接文件称为动态共享对象(DSO,Dynamic Shared Ob ...
- 1020C Elections
传送门 题目大意 现在有 n个人,m个党派,第i个人开始想把票投给党派pi,而如果想让他改变他的想法需要花费ci元.你现在是党派1,问你最少花多少钱使得你的党派得票数大于其它任意党派. 分析 我们枚举 ...
- Luogu 1314 [NOIP2011] 聪明的质监员
二分答案 + 前缀和. 题面中式子的意思是每一个区间$[l, r]$的贡献是这个区间内$w_i \geq W$的个数乘以这些$i$的$v_i$和. 很快发现了答案具有单调性,可以做两遍二分,分别看看小 ...