A Round Peg in a Ground Hole(判断是否是凸包,点是否在凸包内,圆与多边形的关系)
Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 4628 | Accepted: 1434 |
Description
A recent factory run of computer desks were flawed when an automatic grinding machine was mis-programmed. The result is an irregularly shaped hole in one piece that, instead of the expected circular shape, is actually an irregular polygon. You need to figure out whether the desks need to be scrapped or if they can be salvaged by filling a part of the hole with a mixture of wood shavings and glue.
There are two concerns. First, if the hole contains any protrusions (i.e., if there exist any two interior points in the hole that, if connected by a line segment, that segment would cross one or more edges of the hole), then the filled-in-hole would not be structurally sound enough to support the peg under normal stress as the furniture is used. Second, assuming the hole is appropriately shaped, it must be big enough to allow insertion of the peg. Since the hole in this piece of wood must match up with a corresponding hole in other pieces, the precise location where the peg must fit is known.
Write a program to accept descriptions of pegs and polygonal holes and determine if the hole is ill-formed and, if not, whether the peg will fit at the desired location. Each hole is described as a polygon with vertices (x1, y1), (x2, y2), . . . , (xn, yn). The edges of the polygon are (xi, yi) to (xi+1, yi+1) for i = 1 . . . n − 1 and (xn, yn) to (x1, y1).
Input
Line 1 < nVertices > < pegRadius > < pegX > < pegY >
number of vertices in polygon, n (integer)
radius of peg (real)
X and Y position of peg (real)
n Lines < vertexX > < vertexY >
On a line for each vertex, listed in order, the X and Y position of vertex The end of input is indicated by a number of polygon vertices less than 3.
Output
HOLE IS ILL-FORMED if the hole contains protrusions
PEG WILL FIT if the hole contains no protrusions and the peg fits in the hole at the indicated position
PEG WILL NOT FIT if the hole contains no protrusions but the peg will not fit in the hole at the indicated position
Sample Input
5 1.5 1.5 2.0
1.0 1.0
2.0 2.0
1.75 2.0
1.0 3.0
0.0 2.0
5 1.5 1.5 2.0
1.0 1.0
2.0 2.0
1.75 2.5
1.0 3.0
0.0 2.0
1
Sample OutputHOLE IS ILL-FORMED
PEG WILL NOT FIT 题意:给n个点的坐标,以及一个圆的圆心和半径。
判断这n个点能否形成凸包,若不能输出“HOLE IS ILL-FORMED”,否则,再判断圆心以及整个圆是否在凸多边形内; 思路: 分三步;
先形成一个封闭的多边形,
1.判断是否是凸包,循环取相邻的两条边,作叉乘,若结果为负说明不是凸包; 2.判断圆心是否在凸多边形内, 用环顾法:
设圆心为P,逐条枚举n边形的边AB,利用
计算PA和PB的夹角,最后求和得到的就是环顾角。
(1) 圆心在多边形内部时,环顾角=±360
(2) 圆心在多边形外部时,环顾角=0
(3) 圆心在多边形边上时(不包括顶点),环顾角=±180
(4) 圆心在多边形顶点时,环顾角为(0,360)之间的任意角,其实就是圆心所在的顶点的两条邻接边的夹角。
3.判断圆与多边形的关系;
当圆与多边形每条边的距离都小于半径时,说明圆在多边形内部。
#include<stdio.h>
#include<string.h>
#include<math.h>
const double eps = 1e-;
const double PI = 3.1415926535898;
int cmp(double x)
{
if(fabs(x) < eps)
return ;
if(x > )
return ;
return -;
}
int n;//多边形点的个数
double r;//钉子半径 struct Point
{
double x,y;
Point (){}
Point(double a,double b):x(a),y(b) {}
}point[],p;//p为钉子坐标
//叉乘
double det(double x1,double y1,double x2,double y2)
{
return x1*y2-x2*y1;
}
//分别以a、b和 c、d为端点的两条线段的叉乘
double edge_det(Point a,Point b,Point c,Point d)
{
return det(b.x-a.x,b.y-a.y,d.x-c.x,d.y-c.y);
}
//点乘
double dot(double x1,double y1,double x2,double y2)
{
return x1*x2 + y1*y2;
}
//分别以a、p和 b、p为端点的两条线段的点乘
double edge_dot(const Point &p, const Point &a, const Point &b)
{
return dot(a.x-p.x,a.y-p.y,b.x-p.x,b.y-p.y);
}
//两点间的距离
double dis(const Point &a, const Point &b)
{
return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));
}
//向量PA和PB 间的夹角
double cal_angle(const Point &a, const Point &b, const Point &p)
{
double res = edge_dot(p,a,b);
double L1 = dis(a,p);
double L2 = dis(b,p);
return acos(res/(L1*L2));
} bool is_convex_hull()//判断是否为凸包,任意两相临边作叉乘,若结果为负说明不是凸包;
{
double base = ;
for(int i = ; i < n; i++)
{
double tmp = cmp(edge_det(point[i],point[i+],point[i+],point[i+]));
if(!base)
base = tmp;
if(base * tmp < )
return false;
}
return true;
}
bool in_convex_hull()//判断圆心是否在凸包内
{
double angle = ;
for(int i = ; i <= n; i++)
{
if(cmp(edge_det(p,point[i],p,point[i+])) >= )
angle += cal_angle(point[i],point[i+],p);
else angle -= cal_angle(point[i],point[i+],p);
}
if(cmp(angle) == )//在多边形外部
return false;
if(cmp(angle+PI) == || cmp(angle-PI) == )//在多边形边上,不包括顶点。
{
if(cmp(r) == )
return true;
}
if(cmp(angle+*PI) == || cmp(angle-*PI) == )//在多边形内部;
return true;
else//在多边形顶点上;
{
if(cmp(r) == )
return true;
}
return false;
}
bool is_fit()//判断以r为半径的圆是否在凸包内
{
for(int i = ; i < n; i++)
{
double res = fabs(edge_det(point[i],p,point[i+],p)/dis(point[i],point[i+]));
if(cmp(res-r) < )
return false;
}
return true;
}
int main()
{
while(~scanf("%d",&n) && n >= )
{
scanf("%lf %lf %lf",&r,&p.x,&p.y); for(int i = ; i < n; i++)
scanf("%lf %lf",&point[i].x,&point[i].y); point[n] = point[];
point[n+] = point[]; if (!is_convex_hull())//不是一个凸包;
{
printf("HOLE IS ILL-FORMED\n");
continue;
}
bool f1 = in_convex_hull();
bool f2 = is_fit();
if(f1 && f2)
printf("PEG WILL FIT\n");
else printf("PEG WILL NOT FIT\n"); }
return ;
}
A Round Peg in a Ground Hole(判断是否是凸包,点是否在凸包内,圆与多边形的关系)的更多相关文章
- POJ 1584 A Round Peg in a Ground Hole 判断凸多边形,判断点在凸多边形内
A Round Peg in a Ground Hole Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 5456 Acc ...
- POJ 1584 A Round Peg in a Ground Hole[判断凸包 点在多边形内]
A Round Peg in a Ground Hole Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 6682 Acc ...
- POJ 1584 A Round Peg in a Ground Hole 判断凸多边形 点到线段距离 点在多边形内
首先判断是不是凸多边形 然后判断圆是否在凸多边形内 不知道给出的点是顺时针还是逆时针,所以用判断是否在多边形内的模板,不用是否在凸多边形内的模板 POJ 1584 A Round Peg in a G ...
- poj1584 A Round Peg in a Ground Hole 判断多边形凹凸,点到线的距离【基础计算几何】
大致思路:首先对于所给的洞的点,判断是否是凸多边形,图形的输入和输出可以是顺时针或者逆时针,而且允许多点共线 Debug 了好几个小时,发现如下问题 判断三点是否共线,可用斜率公式判断 POINT p ...
- POJ 1584 A Round Peg in a Ground Hole(判断凸多边形,点到线段距离,点在多边形内)
A Round Peg in a Ground Hole Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 4438 Acc ...
- A Round Peg in a Ground Hole(凸包应用POJ 1584)
A Round Peg in a Ground Hole Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 5684 Accepte ...
- POJ 1518 A Round Peg in a Ground Hole【计算几何=_=你值得一虐】
链接: http://poj.org/problem?id=1584 http://acm.hust.edu.cn/vjudge/contest/view.action?cid=22013#probl ...
- POJ 1584 A Round Peg in a Ground Hole【计算几何=_=你值得一虐】
链接: http://poj.org/problem?id=1584 http://acm.hust.edu.cn/vjudge/contest/view.action?cid=22013#probl ...
- POJ 1584:A Round Peg in a Ground Hole
A Round Peg in a Ground Hole Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 5741 Acc ...
随机推荐
- c#中跨线程调用windows窗体控件
c#中跨线程调用windows窗体控件解决. 我们在做winform应用的时候,大部分情况下都会碰到使用多线程控制界面上控件信息的问题.然而我们并不能用传统方法来做这个问题,下面我将详细的介绍.首先来 ...
- oracle多表关联删除数据表记录方法
oracle多表关联删除的两种方法 第一种使用exists方法 delete from tableA where exits ( select 1 from tableB Where tableA.i ...
- Python开发【第二十三篇】:持续更新中...
Python开发[第二十三篇]:持续更新中...
- Unity3D GUI学习之GUILayout控件及使用
GUILayout也可以定义一些控件,并且它们会自动垂直对其: GUILayout.Button("开始游戏"); GUILayout.Button("结束游戏" ...
- java 基本类型和包装类的比较
public class BoxingTest { @Test public void test1(){ String a = new String("1"); String b ...
- (转)PHP下编码转换函数mb_convert_encoding与iconv的使用说明
之--http://www.jb51.net/article/21451.htm mb_convert_encoding这个函数是用来转换编码的.原来一直对程序编码这一概念不理解,不过现在好像有点开窍 ...
- js正则实现用户输入银行卡号的控制及格式化
//js正则实现用户输入银行卡号的控制及格式化 <script language="javascript" type="text/javascript"& ...
- Index Full Scan vs Index Fast Full Scan-1103
[Oracle] Index Full Scan vs Index Fast Full Scan作者:汪海 (Wanghai) 日期:14-Aug-2005 出处:http://spaces.msn. ...
- iOS中常用的正则表达式
iOS常用正则表达式 正则表达式用于字符串处理.表单验证等场合,实用高效.现将一些常用的表达式收集于此,以备不时之需. 匹配中文字符的正则表达式: [\u4e00-\u9fa5]评注:匹配中文还真是个 ...
- sql查询一个班级中总共有多少人以及男女分别多少人
--创建视图 create view StuClassView as SELECT s.ID ,s.StuName ,s.StuAge ,s.StuAddress ,s.StuTel ,s.Class ...