\(\mathcal{Description}\)

  Link.

  给定平面上 \(n\) 个点,求最小的能覆盖其中至少 \(m\) 个点的圆半径及一个可能的圆心。

  \(n\le500\),坐标值 \(X\in[0,10^4]\)。

\(\mathcal{Solution}\)

  不难想到二分答案 \(r\),以每个点为圆心,\(r\) 为半径作圆,若 \(r\) 合法则能找到一个被至少 \(m\) 个圆覆盖的点。

  但是圆的交极难处理,结合数据范围,考虑通过一些枚举操作来简化问题——钦定最终圆心 \(O\) 在以 \(P_i\) 为圆心的圆内,同时考虑到 \(O\) 的最优性,这一限制可以进一步加强为:圆心 \(O\) 在以 \(P_i\) 为圆心的圆周上,于是可以在圆周上按极角做扫描线。到此,我们得到 \(\mathcal O(n^2\log n\log \frac{X}{\epsilon})\) 的算法,应该能过。

  接下来是计算几何算法的常见套路:基于数据随机化降低复杂度。设钦定 \(O\) 在以 \(P_i\) 为圆心的圆周上时,最小的合法半径为 \(r_i\)。假设以 \(i=1..n\) 的顺序依次计算 \(r_i\),此时我们可以用 \(\mathcal O(n\log n)\) 的时间判断某个 \(r_i\) 是否是前缀最小值。若 \(r_i\) 是前缀最小值,我们还需要 \(\mathcal O(n\log n\log\frac{X}{\epsilon})\) 的时间求出 \(r_i\) 的具体值。由于数据随机,所以期望复杂度为

\[\begin{aligned}
T(n) &= \sum_{i=1}^n\left(\frac{\mathcal O\left(n\log n\log\frac{X}{\epsilon}\right)}{i}+\mathcal O(n\log n)\right)\\
&= \mathcal O\left(n\log^2n\log\frac{X}{\epsilon}\right).
\end{aligned}
\]

  相对于期望复杂度的优化,感觉第一步通过适当枚举加强限制,简化问题的思路更为重要。不要一心想着优化复杂度就不敢考虑暴力了。

\(\mathcal{Code}\)

/*+Rainybunny+*/

#include <bits/stdc++.h>

#define rep(i, l, r) for (int i = l, rep##i = r; i <= rep##i; ++i)
#define per(i, r, l) for (int i = r, per##i = l; i >= per##i; --i) typedef std::pair<double, bool> PDB;
#define fi first
#define se second const int MAXN = 500;
const double EPS = 1e-7, PI = acos(-1.);
int n, m; inline double dabs(const double u) { return u < 0. ? -u : u; }
inline int dcmp(const double u) { return dabs(u) < EPS ? 0 : u < 0 ? -1 : 1; }
inline double norm(double u) {
u = fmod(u, 2. * PI);
return u < 0. ? u + 2. * PI : u;
} struct Point {
double x, y;
inline Point operator * (const double k) const {
return { k * x, k * y };
}
inline Point operator + (const Point& u) const {
return { x + u.x, y + u.y };
}
inline bool operator == (const Point& u) const {
return !dcmp(x - u.x) && !dcmp(y - u.y);
}
friend inline double dist(const Point& u, const Point& v) {
return sqrt((u.x - v.x) * (u.x - v.x) + (u.y - v.y) * (u.y - v.y));
}
} pnt[MAXN + 5]; inline bool check(const Point& O, const double R, Point& resO, double& resR) {
static std::vector<PDB> evt; evt.clear();
int cur = 0;
rep (i, 1, n) if (!(pnt[i] == O)) {
double d = dist(O, pnt[i]);
if (d > 2. * R) continue;
double alp = acos(0.5 * d / R);
double the = norm(atan2(pnt[i].y - O.y, pnt[i].x - O.x));
double u = norm(the - alp), v = norm(the + alp);
if (u > v) ++cur;
evt.push_back({ u, 1 }), evt.push_back({ v, 0 });
}
std::sort(evt.begin(), evt.end());
if (cur >= m - 1) {
if (R < resR) resO = O + Point{ 0, resR = R };
return true;
}
for (auto [alp, opt]: evt) {
if ((cur += opt ? 1 : -1) >= m - 1) {
if (R < resR) resO = O + Point{ cos(alp), sin(alp) } * (resR = R);
return true;
}
}
return false;
} int main() {
scanf("%d %d", &n, &m);
rep (i, 1, n) scanf("%lf %lf", &pnt[i].x, &pnt[i].y);
std::shuffle(pnt + 1, pnt + n + 1, std::mt19937(20120712)); double ans = 7.5e3; Point ansP = { -1., -1. };
rep (i, 1, n) {
if (!check(pnt[i], ans, ansP, ans)) continue;
double l = 0., r = ans;
while (l + EPS < r) {
double mid = 0.5 * (l + r);
if (check(pnt[i], mid, ansP, ans)) r = mid;
else l = mid;
}
}
printf("%f\n%f %f\n", ans, ansP.x, ansP.y);
return 0;
}

Solution -「CEOI 2006」「洛谷 P5974」ANTENNA的更多相关文章

  1. 「区间DP」「洛谷P1043」数字游戏

    「洛谷P1043」数字游戏 日后再写 代码 /*#!/bin/sh dir=$GEDIT_CURRENT_DOCUMENT_DIR name=$GEDIT_CURRENT_DOCUMENT_NAME ...

  2. Solution -「JSOI 2019」「洛谷 P5334」节日庆典

    \(\mathscr{Description}\)   Link.   给定字符串 \(S\),求 \(S\) 的每个前缀的最小表示法起始下标(若有多个,取最小的).   \(|S|\le3\time ...

  3. Solution -「洛谷 P4372」Out of Sorts P

    \(\mathcal{Description}\)   OurOJ & 洛谷 P4372(几乎一致)   设计一个排序算法,设现在对 \(\{a_n\}\) 中 \([l,r]\) 内的元素排 ...

  4. Solution -「POI 2010」「洛谷 P3511」MOS-Bridges

    \(\mathcal{Description}\)   Link.(洛谷上这翻译真的一言难尽呐.   给定一个 \(n\) 个点 \(m\) 条边的无向图,一条边 \((u,v,a,b)\) 表示从 ...

  5. Solution -「APIO 2016」「洛谷 P3643」划艇

    \(\mathcal{Description}\)   Link & 双倍经验.   给定 \(n\) 个区间 \([a_i,b_i)\)(注意原题是闭区间,这里只为方便后文描述),求 \(\ ...

  6. 「洛谷4197」「BZOJ3545」peak【线段树合并】

    题目链接 [洛谷] [BZOJ]没有权限号嘤嘤嘤.题号:3545 题解 窝不会克鲁斯卡尔重构树怎么办??? 可以离线乱搞. 我们将所有的操作全都存下来. 为了解决小于等于\(x\)的操作,那么我们按照 ...

  7. 「洛谷3338」「ZJOI2014」力【FFT】

    题目链接 [BZOJ] [洛谷] 题解 首先我们需要对这个式子进行化简,否则对着这么大一坨东西只能暴力... \[F_i=\sum_{j<i} \frac{q_iq_j}{(i-j)^2}-\s ...

  8. 「BZOJ2733」「洛谷3224」「HNOI2012」永无乡【线段树合并】

    题目链接 [洛谷] 题解 很明显是要用线段树合并的. 对于当前的每一个连通块都建立一个权值线段树. 权值线段树处理操作中的\(k\)大的问题. 如果需要合并,那么就线段树暴力合并,时间复杂度是\(nl ...

  9. 「洛谷3870」「TJOI2009」开关【线段树】

    题目链接 [洛谷] 题解 来做一下水题来掩饰ZJOI2019考炸的心情QwQ. 很明显可以线段树. 维护两个值,\(Lazy\)懒标记表示当前区间是否需要翻转,\(s\)表示区间还有多少灯是亮着的. ...

随机推荐

  1. Java基础复习到此结束,统一把源码放到GitHub仓库了,响应开源精神

    这是地址 https://github.com/cen-xi/java-baisc-

  2. 【Maven】maven 插件开发实战

    前言 众所周知,maven 实质上是一个插件执行框架,所有的工作都是通过插件完成的.包括我们日常使用到的类似 install.clean.deploy.compiler...这些命令,其实底层都是一个 ...

  3. 学习鸟哥私房菜--linux bash 的环境变量ps1设置

    bash里边的变量ps1是用户平时的提示符,系统默认为[username@host 工作目录]$.关于ps1的相关介绍详见:http://www.cnblogs.com/starspace/archi ...

  4. 🏆【Alibaba中间件技术系列】「Nacos技术专题」服务注册与发现相关的原理分析

    背景介绍 前几篇文章介绍了Nacos配置中心服务的能力机制,接下来,我们来介绍Nacos另一个非常重要的特性就是服务注册与发现,说到服务的注册与发现相信大家应该都不陌生,在微服务盛行的今天,服务是非常 ...

  5. java-异常-finally代码块

    1 package p1.exception; 2 3 4 class Demo_0{ 5 public int show(int index) throws ArrayIndexOutOfBound ...

  6. ant -design vue a-tree 树形控件

    话不多说,先上代码. <a-tree v-if="this.treeData && this.treeData.length > 0" ref=&quo ...

  7. windows系统配置Nginx使用SSL证书实现Https反向代理

    Nginx反向代理服务,可以代理接收请求,并把请求转发到设置好的其他服务器上. 例如,Nginx服务器为 100.101.102.103,A服务为 100.101.102.104 ,通过Nginx配置 ...

  8. 用8个命令调试Kubernetes集群

    如果使用任何系统的时间足够长,那么你肯定必须对其进行调试,Kubernetes也不例外.它是一个分布式系统,有许多运动部件.我们将介绍8个可以运行以调试任何Kubernetes集群的命令. 它将帮助你 ...

  9. C# 读取txt文件生成Word文档

    本文将以C#程序代码为例介绍如何来读取txt文件中的内容,生成Word文档.在编辑代码前,可参考如下代码环境进行配置: Visual Studio 2017 .Net Framework 4.6.1 ...

  10. 字节Android Native Crash治理之Memory Corruption工具原理与实践

    作者:字节跳动终端技术--庞翔宇 内容摘要 ​ MemCorruption工具是字节跳动AppHealth (Client Infrastructure - AppHealth) 团队开发的一款用于定 ...