20200908 专题：最小圆覆盖

最新推荐文章于 2022-04-26 11:39:08 发布
最新推荐文章于 2022-04-26 11:39:08 发布
阅读量180 收藏
点赞数
分类专栏： 专题 -----常用算法-----
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hongkongreporter/article/details/108477003
版权

总览：

随机增量法
圆的表示： ( ( 圆 心 ( x , y ) ) , 半 径 r ) ((圆心(x,y)),半径r) ((圆心(x,y)),半径r)
初始圆为 ( ( 0 , 0 ) , 0 ) ((0,0),0) ((0,0),0)、
枚举每一个点，找到一个圆，使 p 1 − p n o w p_1-p_{now} p1​−pnow​ 在圆内
枚举每一个点 p i p_i pi​，如果 p i p_i pi​ 不在当前圆内，假设它在圆上，将圆设为 ( ( x i , y i ) , 0 ) ((x_i,y_i),0) ((xi​,yi​),0)
再枚举每一个它之前的点 p j p_j pj​，如果 p j p_j pj​ 不在当前圆内，假设它在圆上，将圆设为 ( ( p i , p j 的 中 点 ) , d i s ( p i , p j ) / 2 ) ((p_i,p_j的中点),dis(p_i,p_j)/2) ((pi​,pj​的中点),dis(pi​,pj​)/2)
再枚举每一个 p j p_j pj​ 之前的点 p k p_k pk​，如果 p k p_k pk​ 不在当前圆内，假设它在圆上，将圆设为 p i , p j , p k 的 三 点 圆 p_i,p_j,p_k的三点圆 pi​,pj​,pk​的三点圆

求三点圆：
就是找外心
找两条边的中垂线求交即可

为防被卡，开始先打乱点的顺序

T1 P1742 最小圆覆盖

思路：
板子题

代码：

#include <bits/stdc++.h>
using namespace std;
#define LL long long
typedef unsigned long long ull;
typedef unsigned int uint;
#define pb push_back
#define mp make_pair
namespace IO {
char buf_[1 << 21], *p1_ = buf_, *p2_ = buf_;
#define ch()                                                                 \
  (p1_ == p2_ &&                                                             \
           (p2_ = (p1_ = buf_) + fread(buf_, 1, 1 << 21, stdin), p1_ == p2_) \
       ? EOF                                                                 \
       : *p1_++)
inline int in() {
  int s = 0, f = 1;
  char x = getchar();
  for (; x < '0' || x > '9'; x = getchar())
    if (x == '-') f = -1;
  for (; x >= '0' && x <= '9'; x = getchar()) s = (s * 10) + (x & 15);
  return f == 1 ? s : -s;
}
char _buf[1 << 21];
int _p1 = -1;
inline void flush() {
  fwrite(_buf, 1, _p1 + 1, stdout);
  _p1 = -1;
}
inline void pc(char x) {
  if (_p1 == (1 << 21) - 1) flush();
  _buf[++_p1] = x;
}
inline void out(int x) {
  char k[30];
  int pos = 0;
  if (!x) {
    pc('0');
    return;
  }
  if (x < 0) {
    pc('-');
    x = -x;
  }
  while (x) {
    k[++pos] = (x % 10) | 48;
    x /= 10;
  }
  for (int i = pos; i; i--) pc(k[i]);
  return;
}
inline void out(string x) {
  int len = x.size();
  for (int i = 0; i < len; i++) pc(x[i]);
}
}  // namespace IO
using namespace IO;

const int A = 3e5 + 5;
const double EPS = 1e-15;
int n;
inline int cp(double u, double v) {
  if (fabs(u - v) < EPS) return 0;
  return u > v ? 1 : -1;
}
struct Point {
  double x, y;
  Point(double u = 0, double v = 0) { x = u, y = v; }
  inline friend Point operator+(Point u, Point v) {
    return Point(u.x + v.x, u.y + v.y);
  }
  inline friend Point operator-(Point u, Point v) {
    return Point(u.x - v.x, u.y - v.y);
  }
  inline friend Point operator*(Point u, double v) {
    return Point(u.x * v, u.y * v);
  }
  inline friend Point operator/(Point u, double v) {
    return Point(u.x / v, u.y / v);
  }
  inline friend double operator^(Point u, Point v) {
    return u.x * v.y - u.y * v.x;
  }
  inline friend double operator&(Point u, Point v) {
    return u.x * v.x + u.y * v.y;
  }
  inline friend bool operator==(Point u, Point v) {
    return !cp(u.x, v.x) && !cp(u.y, v.y);
  }
  inline friend bool operator!=(Point u, Point v) { return !(u == v); }
} p[A];
struct Line {
  Point x, y;
  Line(Point u = 0, Point v = 0) { x = u, y = v; }
};
struct Circle {
  Point x;
  double r;
  Circle(Point u = 0, double v = 0) { x = u, r = v; }
} mn;

inline double dis(Point u, Point v) {
  return sqrt((u.x - v.x) * (u.x - v.x) + (u.y - v.y) * (u.y - v.y));
}

inline int check(Point u, Circle v) { return cp(dis(u, v.x), v.r) != 1; }

inline Point inter(Line u, Line v) {
  double s1 = (v.x - u.x) ^ (v.y - u.x), s2 = (v.y - u.y) ^ (v.x - u.y);
  return u.x + (u.y - u.x) * s1 / (s1 + s2);
}

inline Line find(Line u) {
  Point a = (u.x + u.y) / 2;
  return Line(a, Point(a.x + (u.y - u.x).y, a.y - (u.y - u.x).x));
}

inline Circle circle(Point u, Point v, Point w) {
  Line a = find(Line(u, v)), b = find(Line(v, w));
  Point c = inter(a, b);
  return Circle(c, dis(u, c));
}

signed main() {
  srand(time(0));
  n = in();
  for (int i = 1; i <= n; i++) scanf("%lf%lf", &p[i].x, &p[i].y);
  random_shuffle(p + 1, p + 1 + n);
  mn = Circle(Point(0, 0), 0);
  for (int i = 1; i <= n; i++)
    if (!check(p[i], mn)) {
      mn = Circle(p[i], 0);
      for (int j = 1; j <= i - 1; j++)
        if (!check(p[j], mn)) {
          mn = Circle((p[i] + p[j]) / 2, dis(p[i], p[j]) / 2);
          for (int k = 1; k <= j - 1; k++) {
            if (!check(p[k], mn)) mn = circle(p[i], p[j], p[k]);
          }
        }
    }
  printf("%0.10lf\n%0.10lf %0.10lf\n", mn.r, mn.x.x, mn.x.y);
  flush();
  return 0;
}

T2 P4288 [SHOI2014]信号增幅仪

思路：
先旋转坐标轴，将椭圆转正，然后放缩，将椭圆变成圆

代码：

#include <bits/stdc++.h>
using namespace std;
#define LL long long
typedef unsigned long long ull;
typedef unsigned int uint;
#define pb push_back
#define mp make_pair
namespace IO {
char buf_[1 << 21], *p1_ = buf_, *p2_ = buf_;
#define ch()                                                                 \
  (p1_ == p2_ &&                                                             \
           (p2_ = (p1_ = buf_) + fread(buf_, 1, 1 << 21, stdin), p1_ == p2_) \
       ? EOF                                                                 \
       : *p1_++)
inline int in() {
  int s = 0, f = 1;
  char x = ch();
  for (; x < '0' || x > '9'; x = ch())
    if (x == '-') f = -1;
  for (; x >= '0' && x <= '9'; x = ch()) s = (s * 10) + (x & 15);
  return f == 1 ? s : -s;
}
char _buf[1 << 21];
int _p1 = -1;
inline void flush() {
  fwrite(_buf, 1, _p1 + 1, stdout);
  _p1 = -1;
}
inline void pc(char x) {
  if (_p1 == (1 << 21) - 1) flush();
  _buf[++_p1] = x;
}
inline void out(int x) {
  char k[30];
  int pos = 0;
  if (!x) {
    pc('0');
    return;
  }
  if (x < 0) {
    pc('-');
    x = -x;
  }
  while (x) {
    k[++pos] = (x % 10) | 48;
    x /= 10;
  }
  for (int i = pos; i; i--) pc(k[i]);
  return;
}
inline void out(string x) {
  int len = x.size();
  for (int i = 0; i < len; i++) pc(x[i]);
}
}  // namespace IO
using namespace IO;

const int A = 3e5 + 5;
const double EPS = 1e-7;
const double pi = acos(-1);
int n;
inline int cp(double u, double v) {
  if (fabs(u - v) < EPS) return 0;
  return u > v ? 1 : -1;
}
struct Point {
  double x, y;
  Point(double u = 0, double v = 0) { x = u, y = v; }
  inline friend Point operator+(Point u, Point v) {
    return Point(u.x + v.x, u.y + v.y);
  }
  inline friend Point operator-(Point u, Point v) {
    return Point(u.x - v.x, u.y - v.y);
  }
  inline friend Point operator*(Point u, double v) {
    return Point(u.x * v, u.y * v);
  }
  inline friend Point operator/(Point u, double v) {
    return Point(u.x / v, u.y / v);
  }
  inline friend double operator^(Point u, Point v) {
    return u.x * v.y - u.y * v.x;
  }
  inline friend double operator&(Point u, Point v) {
    return u.x * v.x + u.y * v.y;
  }
  inline friend bool operator==(Point u, Point v) {
    return !cp(u.x, v.x) && !cp(u.y, v.y);
  }
  inline friend bool operator!=(Point u, Point v) { return !(u == v); }
} p[A];
struct Line {
  Point x, y;
  Line(Point u = 0, Point v = 0) { x = u, y = v; }
};
struct Circle {
  Point x;
  double r;
  Circle(Point u = 0, double v = 0) { x = u, r = v; }
} mn;

inline double dis(Point u, Point v) {
  return sqrt((u.x - v.x) * (u.x - v.x) + (u.y - v.y) * (u.y - v.y));
}

inline int check(Point u, Circle v) { return cp(dis(u, v.x), v.r) != 1; }

inline Point inter(Line u, Line v) {
  double s1 = (v.x - u.x) ^ (v.y - u.x), s2 = (v.y - u.y) ^ (v.x - u.y);
  return u.x + (u.y - u.x) * s1 / (s1 + s2);
}

inline Line find(Line u) {
  Point a = (u.x + u.y) / 2;
  return Line(a, Point(a.x + (u.y - u.x).y, a.y - (u.y - u.x).x));
}

inline Circle circle(Point u, Point v, Point w) {
  Line a = find(Line(u, v)), b = find(Line(v, w));
  Point c = inter(a, b);
  return Circle(c, dis(u, c));
}

inline Point rotate(Point u, double phi) {
  return Point(cos(phi) * u.x - sin(phi) * u.y,
               sin(phi) * u.x + cos(phi) * u.y);
}

signed main() {
  srand(time(0));
  n = in();
  for (int i = 1; i <= n; i++) p[i].x = in(), p[i].y = in();
  double phi = in();
  phi = -pi * (phi / 180);
  for (int i = 1; i <= n; i++) p[i] = rotate(p[i], phi);
  int d = in();
  for (int i = 1; i <= n; i++) p[i].x = p[i].x / d;
  random_shuffle(p + 1, p + 1 + n);
  mn = Circle(Point(0, 0), 0);
  for (int i = 1; i <= n; i++)
    if (!check(p[i], mn)) {
      mn = Circle(p[i], 0);
      for (int j = 1; j <= i - 1; j++)
        if (!check(p[j], mn)) {
          mn = Circle((p[i] + p[j]) / 2, dis(p[i], p[j]) / 2);
          for (int k = 1; k <= j - 1; k++) {
            if (!check(p[k], mn)) mn = circle(p[i], p[j], p[k]);
          }
        }
    }
  printf("%0.3lf\n", mn.r);
  flush();
  return 0;
}
ハルカナソラヘ
关注
专栏目录
【暖*墟】#圆的最大覆盖面积# 洛谷P1378 油滴扩展
flora715的博客
07-09 478
P1378 油滴扩展 跟随某dalao一起做题 马上遇到了我不会的有关圆覆盖的问题qaq 其实一直都不会做这类题目&amp;&amp; 所以准备写个专题研究一下 先来看看这个p1378 #include &lt;bits/stdc++.h&gt; using namespace std; typedef long long ll; /* 【油滴扩散】 在一个长方形框子里有N（0≤N≤6...
HYSBZ - 3564 信号增幅仪 【椭圆最小覆盖】【坐标转换】
slass-mot-odet
08-13 421
无线网络基站在理想状况下有效信号覆盖范围是个圆形。而无线基站的功耗与圆的半径的平方成正比。现给出平面 上若干网络用户的位置,请你选择一个合适的位置建设无线基站…就在你拿起键盘准备开始敲代码的时候,你的好 朋友发明家 SHTSC 突然出现了。SHTSC 刚刚完成了他的新发明——无线信号增幅仪。增幅仪能够在不增加无线基 站功耗的前提下,使得有效信号的覆盖范围在某一特定方向上伸长若干倍。即:使用了增幅仪...
最小圆覆盖
weixin_30779691的博客
10-15 67
1 struct Point{double x,y;}; 2 struct Circle{Point c;double r;}; 3 double dist(Point a,Point b){ 4 return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); 5 } 6 C...
BZOJ 3564 SHOI 2014 信号增幅仪 坐标变换+最小圆覆盖
(╯°口°)╯(┴—┴
02-26 1378
题目大意：给出平面上的一些点，现在让你用一个长轴与x轴成一定角度的，长轴：短轴已知的椭圆来覆盖所有的坐标，求最小的短轴长度。 思路：很明显，这个椭圆的形状和放置状态已经给出了，但是没有办法求最小拖圆覆盖啊。采用坐标变换，将椭圆变成圆。首先我们先让长轴与x轴平行，将平面上的所有点都旋转这个角度。之后只需要让所有点的x坐标除以长轴：短轴就可以了。剩下的就是最小圆覆盖了。 注：坐标旋转公式：
Python实现小圆覆盖大圆问题
Y
07-25 6921
问题回顾 给定一个半径为10m的圆，和若干半径为1m的小圆（小圆个数足够多），现在要用若干小圆来覆盖这个大圆。问至少需要多少个小圆才能覆盖这个半径为10m的圆？ 求解思路 为保证小圆数量最小，既圆与圆之间不重叠，既相切情况。 因为需要考虑大圆与小圆之间的间隙，也要考虑小圆与小圆之间的间隙。 故只能采取近似思想进行计算。 1、当r<R<2r时，n=1; 2、当2r<=R<2.1547r时，n=2; 3、...
(新版)新人教版中考数学总复习专题检测21与圆有关的计算试题.pdf
10-17
14. **覆盖正三角形的最小圆**：要覆盖一个正三角形，最小的圆应以三角形的内切圆为半径，所以最小半径是正三角形边长的一半。 解答题部分主要涉及的是几何图形的性质、相似三角形、比例关系、圆周角定理和面积周长...
2019年高考数学二轮复习专题突破课时作业14直线圆理
08-05
这些问题覆盖了直线和圆的基本知识点，包括直线的倾斜角、截距、平行性、点到直线的距离、直线与圆的位置关系以及圆的标准方程等。在复习过程中，掌握这些基本概念和解题方法对提高解题能力至关重要。
2019年高考数学大二轮复习专题六解析几何6.1直线与圆练习
08-05
总结来说，这些题目覆盖了直线与圆的基本性质，如平行、相交、切线、反射以及相关距离和面积的计算，这些都是解析几何中的核心概念。通过这些练习，学生可以更好地掌握解析几何的理论和方法，为高考数学做好充分准备...
2016高考数学专题复习导练测第九章第2讲圆的方程理新人教A版
08-19
第九题涉及一个平面区域恰好被面积最小的圆覆盖，这个圆是该三角形的外接圆。找到三角形的顶点，然后求解外接圆的圆心和半径。 这些题目共同强调了圆的方程求解、圆与直线的关系、圆心的确定、距离与半径的关系以及...
最小圆覆盖 半平面交 凸包 三维几何
05-31
最小圆覆盖 半平面交 凸包 Pick公式多边形内整点个数 这个平面多边形的重心 三维几何库 旋转卡壳求最近点对
BZOJ3564: [SHOI2014]信号增幅仪（最小圆覆盖）
DZYO的博客
10-07 584
传送门题意： 求一个固定半长轴方向和半长轴半短轴比值的最小椭圆覆盖平面上所有点。题解： 最小圆覆盖。 首先将坐标旋转，再把所有横坐标除以比值得到最小圆覆盖的经典模板。最小圆覆盖的做法： 首先，在已经确定最小圆的点集中加入一个在最小圆外面的点，那么这个点一定在新点集的最小圆上。算法： 初始化圆形为p1p_1，半径为r1r_1 1.从左到右依次加入点，若在当前圆内则跳过，否则问题转化为在一个
小圆实现平面全覆盖（密排、小圆覆盖椭圆）
weixin_53488599的博客
04-26 4511
实现定长小圆全覆盖椭圆，平面密排
matlab画图：实现小圆覆盖大圆
qq_41845823的博客
09-16 5305
python实现小圆覆盖大圆：[https://blog.csdn.net/qq_41845823/article/details/118860729](https://blog.csdn.net/qq_41845823/article/details/118860729) *以下是MATLAB的实现代码： function []=plot1(x,y,r) theta=0:0.01:2*pi; Circle1=x+r*cos(theta); Circle2=y+r*sin(theta); c=[123,14
最小圆覆盖问题 算法步骤与证明+代码模板
热门推荐
harry12138小白的博客
03-13 1万+
最小圆覆盖定义 最小圆覆盖问题是什么呢？就是指在二维平面上有一堆点，然后我们要求一个最小半径的圆能够将所有点全部都包住，这就是最小圆覆盖问题。 最小覆盖圆的性质 性质1：最小覆盖圆是唯一的 证明：我们假设有两个圆O1，O2，他们半径都是r，都是最小覆盖圆，那么所有的点一定在两圆的交集部分。那我们以两圆交集部分的弦长为直径，做一个新圆，该圆依旧包含所有点，而且他的直径是原来圆的弦长，一定是小于原来圆的直径的，因此和原来的圆是最小覆盖圆相矛盾，所以最小覆盖圆是唯一的。 性质2：若圆O1是点集S的最小覆盖圆，则
最小圆覆盖(计算几何基础)
Satur9的博客
04-20 924
传送门 思路： ① 随机化点集（可降低复杂度） ② 初始随意找到两点，设为和，以和为直径做出初始圆(以 表示覆盖前 i 个点的最小圆)。 ③ 依次加入新的点，若当前点在圆中，则，否则以作为新的直径端点做圆(显然点在新圆的边界上)。 ④ 但是新圆又不一定能包含之前的所有 i 个点，于是我们遍历找到前 i 个点中不在中的点，那么点和一定都会在新圆的边界上，于是直接以为直径做出新圆。 ⑤ 同样也不一定能包含之前的所有 j...
凸优化第八章几何问题 8.4极值体积椭圆
清风吹斜阳
02-22 2221
8.4极值体积椭圆 Lowner-John椭球 最大体积内接椭球 椭球逼近的效率 Lowner-John椭球 包含集合C的最小体积椭球被成为集合C的Lowner-John椭球，记为，为方便描述的特征，将一般的椭球参数化为 即Euclid球在仿射映射下的原象。可以不是一般性地假设，此时的体积正比于。计算包含C的最小体积椭球的问题可以表述为： 其中，且存在一个隐含约束。目标函数和约束函...
P4288 [SHOI2014]信号增幅仪 最小圆覆盖
m0_51068403的博客
02-17 121
传送门 题意：给出一些点，让后让你用一个ba=1p\frac{b}{a}=\frac{1}{p}ab​=p1​的且绕x正方向逆时针转动了a°的椭圆覆盖，求最小的b。 由于题目中椭圆是绕x正方向逆时针旋转了a°，为了方便写出方程，我们可以先把它顺时针转回来，这样就得到椭圆方程x2a2+y2b2=1\frac {x^{2}}{a^{2}}+\frac{y^{2}}{b^{2}}=1a2x2​+b2y2​=1，因为ba=1p\frac{b}{a}=\frac{1}{p}ab​=p1​，我们可以让a=pa=pa=p
4. 圆面覆盖 (40分)
abao_的专栏
06-04 882
4. 圆面覆盖 (40分)问题背景在平面上有一个长为L，宽为W的长方形，左下角坐标为(0,0)，右上角坐标为(L,W)。给定一些圆，第i个圆的圆心坐标为(xi,yi)，半径为Ri。你的任务是求最小的正实数k，使得把每个圆的半径变为原来的k倍后（即：第i个圆半径变为kRi，圆心位置不变），长方形将被这些圆完全覆盖。换句话说，长方形内部或边界上的任意点均至少在一个圆的内部或边界上。输入
无线传感器网络覆盖算法：最小单位圆集广播数据转发
它通过最小数量的邻居节点子集实现所有节点的覆盖，并证明了最小单位圆覆盖集与轮廓集相等。实验结果表明，此算法在减少节点使用数量的同时，提高了网络覆盖和数据转发效率，从而延长了网络寿命。" 这篇研究论文...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值