链接:

http://poj.org/problem?id=1584

题意:

按照顺时针或逆时针方向输入一个n边形的顶点坐标集,先判断这个n边形是否为凸包。

再给定一个圆形(圆心坐标和半径),判断这个圆是否完全在n变形内部。

题解:

1、判断凸包convex():

由于点集已经按某个时针方向有序,因此可以先定义一个方向系数dir=0

两两枚举n边形的边,用叉积判断这两条边的转向(右螺旋或左螺旋),由于存在散点共线的情况,因此当且仅当叉积的值t第一次不为0时,dir=t,dir的值此后不再改变。(dir>0 则为右螺旋逆时针,dir<0则为左螺旋顺时针)

此后继续枚举剩下的边,只要判断dir*t>=0即可,当存在一个dir*t<0的边,说明这是凹多边形,就不是凸包了。

2、判断圆心在不在凸包内contain()

3、当圆心在凸包内时,判断距离是否都大于半径fit()

代码:

  1 #include <map>
  2 #include <set>
  3 #include <cmath>
  4 #include <queue>
  5 #include <stack>
  6 #include <cstdio>
  7 #include <string>
  8 #include <vector>
  9 #include <cstdlib>
 10 #include <cstring>
 11 #include <sstream>
 12 #include <iostream>
 13 #include <algorithm>
 14 #include <functional>
 15 using namespace std;
 16 #define rep(i,a,n) for (int i=a;i<n;i++)
 17 #define per(i,a,n) for (int i=n-1;i>=a;i--)
 18 #define all(x) (x).begin(),(x).end()
 19 #define pb push_back
 20 #define mp make_pair
 21 #define lson l,m,rt<<1
 22 #define rson m+1,r,rt<<1|1
 23 typedef long long ll;
 24 typedef vector<int> VI;
 25 typedef pair<int, int> PII;
 26 const ll MOD = 1e9 + 7;
 27 const int INF = 0x3f3f3f3f;
 28 const int MAXN = 1010;
 29 // head
 30
 31 const double eps = 1e-8;
 32 int cmp(double x) {
 33     if (fabs(x) < eps) return 0;
 34     if (x > 0) return 1;
 35     return -1;
 36 }
 37
 38 const double pi = acos(-1);
 39 inline double sqr(double x) {
 40     return x*x;
 41 }
 42 struct point {
 43     double x, y;
 44     point() {}
 45     point(double a, double b) :x(a), y(b) {}
 46     void input() {
 47         scanf("%lf%lf", &x, &y);
 48     }
 49     friend point operator+(const point &a, const point &b) {
 50         return point(a.x + b.x, a.y + b.y);
 51     }
 52     friend point operator-(const point &a, const point &b) {
 53         return point(a.x - b.x, a.y - b.y);
 54     }
 55     friend point operator*(const double &a, const point &b) {
 56         return point(a*b.x, a*b.y);
 57     }
 58     friend point operator/(const point &a, const double &b) {
 59         return point(a.x / b, a.y / b);
 60     }
 61     double norm() {
 62         return sqrt(sqr(x) + sqr(y));
 63     }
 64 };
 65 double det(point a, point b) {
 66     return a.x*b.y - a.y*b.x;
 67 }
 68 double dot(point a, point b) {
 69     return a.x*b.x + a.y*b.y;
 70 }
 71 double dist(point a, point b) {
 72     return (a - b).norm();
 73 }
 74
 75 struct line {
 76     point a, b;
 77     line() {}
 78     line(point x, point y) :a(x), b(y) {}
 79 };
 80 double dis_point_segment(point p, point s, point t) {
 81     if (cmp(dot(p - s, t - s)) < 0) return (p - s).norm();
 82     if (cmp(dot(p - t, s - t)) < 0) return (p - t).norm();
 83     return fabs(det(s - p, t - p) / dist(s, t));
 84 }
 85 bool point_on_segment(point p, point s, point t) {
 86     return cmp(det(p - s, t - s)) == 0 && cmp(dot(p - s, p - t)) <= 0;
 87 }
 88 bool parallel(line a, line b) {
 89     return !cmp(det(a.a - a.b, b.a - b.b));
 90 }
 91 bool line_make_point(line a, line b,point &res) {
 92     if (parallel(a, b)) return false;
 93     double s1 = det(a.a - b.a, b.b - b.a);
 94     double s2 = det(a.b - b.a, b.b - b.a);
 95     res = (s1*a.b - s2*a.a) / (s1 - s2);
 96     return true;
 97 }
 98
 99 int n;
100 double r;
101 point o;
102 point *p;
103
104 bool convex() {
105     int dir = 0;
106     rep(i, 0, n) {
107         int t = cmp(det(p[i + 1] - p[i], p[i + 2] - p[i + 1]));
108         if (!dir) dir = t;
109         if (dir*t < 0) return false;
110     }
111     return true;
112 }
113
114 bool contain() {
115     int sign = 0;
116     rep(i, 0, n) {
117         int x = cmp(det(p[i] - o, p[i + 1] - o));
118         if (x) {
119             if (sign && sign != x) return false;
120             else sign = x;
121         }
122     }
123     return true;
124 }
125
126 bool fit() {
127     rep(i, 0, n) {
128         int k = cmp(dis_point_segment(o, p[i], p[i + 1]) - r);
129         if (k < 0) return false;
130     }
131     return true;
132 }
133
134 int main() {
135     while (cin >> n && n != 1) {
136         cin >> r;
137         o.input();
138         p = new point[n + 2];
139         rep(i, 1, n + 1) p[i].input();
140         p[0] = p[n];
141         p[n + 1] = p[1];
142         if(!convex()) cout << "HOLE IS ILL-FORMED" << endl;
143         else {
144             bool fg1 = contain();
145             bool fg2 = fit();
146             if (fg1 && fg2) cout << "PEG WILL FIT" << endl;
147             else cout << "PEG WILL NOT FIT" << endl;
148         }
149         delete p;
150     }
151     return 0;
152 }

转载于:https://www.cnblogs.com/baocong/p/6730953.html

POJ 1584 计算几何 凸包相关推荐

  1. POJ 1584 A Round Peg in a Ground Hole 判断凸多边形,点到线段距离,点在多边形内

    ACM博客_kuangbin POJ 1584 A Round Peg in a Ground Hole(判断凸多边形,点到线段距离,点在多边形内) A Round Peg in a Ground H ...

  2. [Poj 2187] 计算几何之凸包(二) {更高效的算法}

    { 承上一节 继续介绍点集的凸包  (下文中所有凸包 若不做特殊说明均指点集的凸包) 这一节介绍相比更高效的算法 } ========================================= ...

  3. poj 1584(凸包+点在凸多边形内+圆在凸多边形内)

    题意: 按照顺时针或逆时针方向输入一个n边形的顶点坐标集,先判断这个n边形是否为凸包. 再给定一个圆形(圆心坐标和半径),判断这个圆是否完全在n边形内部. 解题思路: 1.判断该多边形是否是凸包,由于 ...

  4. poj 1113 Wall 凸包的应用

    题目链接:poj 1113   单调链凸包小结 题解:本题用到的依然是凸包来求,最短的周长,只是多加了一个圆的长度而已,套用模板,就能搞定: AC代码: 1 #include<iostream& ...

  5. POJ 1228 —— “稳定”凸包

    POJ 1228 Grandpa's Estate 这是个好题目,同时也是个不和谐的题目(不和谐原因是题目出的存在漏洞,数据弱,而且有些条件没给清楚,为了一个SB错误无限WA之后,终于AC) 题意就废 ...

  6. 凸包计算几何matlab,计算几何-凸包算法 Python实现与Matlab动画演示

    凸包算法是计算几何中的最经典问题之一了.给定一个点集,计算其凸包.凸包是什么就不罗嗦了 本文给出了<计算几何--算法与应用>中一书所列凸包算法的Python实现和Matlab实现,并给出了 ...

  7. poj 1873 枚举+凸包

    题意:给出一些树的位置,价值,长度,现要求先砍一些树制成一定长度的篱笆将剩余的树围起来,求要砍树的最小总价值. 因为树的个数最多为15个很容易想到用二进制数表示树的状态进行遍历. #include&l ...

  8. POJ 3348 Cows 凸包面积

    题目大意:连接个点,求围成的面积并除以50(向下取整) 题目思路:就是求凸包面积:按逆时针方向为凸包上每条边指定方向,对于每条边AB,累加(AXB)/2的值. #include<cstdio&g ...

  9. POJ - 1584 A Round Peg in a Ground Hole(综合几何)

    题目链接:点击查看 题目大意:给出n个点,以及一个圆心和半径,首先判断这n个点能否构成凸包,若能,继续判断圆是否在凸包内 题目分析:这个题目确实非常综合,考察了判断凸包问题,判断圆是否在凸包内,只要保 ...

最新文章

  1. 服务器负载均衡的基本功能和实现原理
  2. 使用迭代器从map或vector中删除元素
  3. 多线程售票demo,用ReentrantLock实现
  4. C语言试题七十六之请编写函实现百钱买百鸡
  5. 【软件开发底层知识修炼】十七 快速学习GDB调试四 使用GDB进行函数调用栈的查看
  6. 2018上半年软件设计师上午题参考答案
  7. C#设置系统日期时间格式
  8. IO中同步、异步与阻塞、非阻塞的区别
  9. Create-React-App创建antd-mobile开发环境
  10. 让电脑只能上允许的QQ号
  11. js页面打印去掉页眉页脚
  12. python列表append方法_Python列表append()方法
  13. 剑指Java面试:面试官能问到的问题,都被我收集在这份PDF文档里
  14. BI是什么,BI的解释
  15. 个人成长【提升认知销量秘籍】
  16. Keras中Dropout的用法详解
  17. 数据全量、增量、比较更新
  18. 【实战】聊聊幂等设计
  19. 全国计算机等级考试二级公共基础知识考试大纲(2022年版)
  20. R语言实现LDA算法(鸢尾花)

热门文章

  1. 激光位移传感器适用于那些领域
  2. 用HTML5开发打包本地桌面应用
  3. 阿里巴巴《Java开发手册》学习总结分享
  4. sw装配体转换成一个整体的stl文件
  5. (三)集成学习上——偏差与方差
  6. word中怎么快速选中光标之前或之后的全部内容?
  7. PRINCE2学习笔记-泛读1
  8. Linux设置断网自动关机(蓄电池)
  9. 网页设计制作标准规范
  10. 如何解决Office365无法登录