POJ 1584 计算几何 凸包
链接:
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 计算几何 凸包相关推荐
- 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 ...
- [Poj 2187] 计算几何之凸包(二) {更高效的算法}
{ 承上一节 继续介绍点集的凸包 (下文中所有凸包 若不做特殊说明均指点集的凸包) 这一节介绍相比更高效的算法 } ========================================= ...
- poj 1584(凸包+点在凸多边形内+圆在凸多边形内)
题意: 按照顺时针或逆时针方向输入一个n边形的顶点坐标集,先判断这个n边形是否为凸包. 再给定一个圆形(圆心坐标和半径),判断这个圆是否完全在n边形内部. 解题思路: 1.判断该多边形是否是凸包,由于 ...
- poj 1113 Wall 凸包的应用
题目链接:poj 1113 单调链凸包小结 题解:本题用到的依然是凸包来求,最短的周长,只是多加了一个圆的长度而已,套用模板,就能搞定: AC代码: 1 #include<iostream& ...
- POJ 1228 —— “稳定”凸包
POJ 1228 Grandpa's Estate 这是个好题目,同时也是个不和谐的题目(不和谐原因是题目出的存在漏洞,数据弱,而且有些条件没给清楚,为了一个SB错误无限WA之后,终于AC) 题意就废 ...
- 凸包计算几何matlab,计算几何-凸包算法 Python实现与Matlab动画演示
凸包算法是计算几何中的最经典问题之一了.给定一个点集,计算其凸包.凸包是什么就不罗嗦了 本文给出了<计算几何--算法与应用>中一书所列凸包算法的Python实现和Matlab实现,并给出了 ...
- poj 1873 枚举+凸包
题意:给出一些树的位置,价值,长度,现要求先砍一些树制成一定长度的篱笆将剩余的树围起来,求要砍树的最小总价值. 因为树的个数最多为15个很容易想到用二进制数表示树的状态进行遍历. #include&l ...
- POJ 3348 Cows 凸包面积
题目大意:连接个点,求围成的面积并除以50(向下取整) 题目思路:就是求凸包面积:按逆时针方向为凸包上每条边指定方向,对于每条边AB,累加(AXB)/2的值. #include<cstdio&g ...
- POJ - 1584 A Round Peg in a Ground Hole(综合几何)
题目链接:点击查看 题目大意:给出n个点,以及一个圆心和半径,首先判断这n个点能否构成凸包,若能,继续判断圆是否在凸包内 题目分析:这个题目确实非常综合,考察了判断凸包问题,判断圆是否在凸包内,只要保 ...
最新文章
- 服务器负载均衡的基本功能和实现原理
- 使用迭代器从map或vector中删除元素
- 多线程售票demo,用ReentrantLock实现
- C语言试题七十六之请编写函实现百钱买百鸡
- 【软件开发底层知识修炼】十七 快速学习GDB调试四 使用GDB进行函数调用栈的查看
- 2018上半年软件设计师上午题参考答案
- C#设置系统日期时间格式
- IO中同步、异步与阻塞、非阻塞的区别
- Create-React-App创建antd-mobile开发环境
- 让电脑只能上允许的QQ号
- js页面打印去掉页眉页脚
- python列表append方法_Python列表append()方法
- 剑指Java面试:面试官能问到的问题,都被我收集在这份PDF文档里
- BI是什么,BI的解释
- 个人成长【提升认知销量秘籍】
- Keras中Dropout的用法详解
- 数据全量、增量、比较更新
- 【实战】聊聊幂等设计
- 全国计算机等级考试二级公共基础知识考试大纲(2022年版)
- R语言实现LDA算法(鸢尾花)