poj1584

最新推荐文章于 2019-09-12 22:04:55 发布

小小小小杜

最新推荐文章于 2019-09-12 22:04:55 发布

阅读量1.1k

点赞数

分类专栏： acm 计算几何文章标签：凸包判定+点与多边形的位置关系+判断圆是

本文链接：https://blog.csdn.net/i_want_to_be_a_god/article/details/38278599

版权

acm 同时被 2 个专栏收录

178 篇文章 1 订阅

订阅专栏

计算几何

6 篇文章 0 订阅

订阅专栏

题目不难，主要是涉及的知识点多，判断凸包+判断圆是否在多边形内（判断点与多边形的关系）

题意为给定n个点，和圆的半径，圆心坐标。然后给出多边形的n个顶点。判断

1）若多边形不是凸包，则输出“HOLE IS ILL-FORMED”

2）若是凸包，但圆不再凸包内，则输出“PEG WILL NOT FIT”

3）既是凸包，圆又在凸包内。输出“PEG WILL FIT”

那么关键点就是如何判断多边形是否为凸包和判段圆是否在多边形内部了。

1）判断多边形是否为凸包

可以用叉积，cross(node[i],node[i+1],node[i+1],node[i+2])的方向要一致（都为顺时针或都为逆时针），如果其中有一个不一致，则说明不为凸包。要注意，可能有顶点共线。所以可以取dir为cross第一次不为0时的值，那么接下来cross*di的值必须要>=0，否则就不是凸包。

2）判断圆是否在多边形内部

分两个阶段判断：

21）判断圆心是否在多边形内部

可以使用角度法，判定：（用点积）

若角度=+360 或 -360 则说明点在多边形内

若角度=0则说明在多边形外部

其余情况说明在多边形的边上（包括顶点）

22）判断圆心到多边形所有点的距离是否大于或等于半径

用叉积求距离

下面是代码： 168K+0MS

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define eps 1e-6 // 最小精度限制
#define Max 1010 // 设置最大多边形顶点个数
#define PI 3.141592654 //PI尽可能精度取高
typedef struct Point{ // 点结构
	double x;
	double y;
}point;
point node[Max];
point cir; // 圆坐标
int n;
double cirr; //圆半径
int dblcmp(double x){ // 判断正负零函数，一般double类型数据做加减法且要进行判断时，都要用到
	if(fabs(x)<eps)
		return 0;
	return x>0?1:-1;
}
double dotdet(double x1,double y1,double x2,double y2){ // 计算点积（x1,y1)点积(x2,y2)
	return x1*x2+y1*y2;
}
double det(double x1,double y1,double x2,double y2){ // 计算叉积 (x1,y1)叉积(x2,y2)
	return x1*y2-x2*y1;
}
double cross(point A,point B,point C,point D){ // 计算叉积，AB叉积CD
	return det(B.x-A.x,B.y-A.y,D.x-C.x,D.y-C.y);
}
double dist(point A,point B){ 计算点A与B距离
	return sqrt((B.x-A.x)*(B.x-A.x)+(B.y-A.y)*(B.y-A.y));
}
double angle(point A,point B){ // 计算圆心与点A、B构成三角形的角度，圆心的角度
	return acos(dotdet(A.x-cir.x,A.y-cir.y,B.x-cir.x,B.y-cir.y)/(dist(A,cir)*dist(B,cir)));
}
bool Is_covbag(){ //判断多边形是否为凸包，其中node可顺时针存放顶点也可逆时针存放顶点
	int dir=0;
	for(int i=0;i<=n-1;i++){
		int temp=dblcmp(cross(node[i],node[i+1],node[i+1],node[i+2]));
	    if(dir==0) //第一次不为0，赋值
			dir=temp;
		if(dir*temp<0) //若和第一次的方向不同，说明不是凸包
			return false;
	}
	return true; //为凸包
}
bool Is_In(){ //判断点与多边形的位置关系
 	double ang=0.0; //角度初始化为0
	for(int i=1;i<=n;i++) 
		if(dblcmp(cross(cir,node[i],cir,node[i+1]))>=0) //若为顺时针存放
			ang+=angle(node[i],node[i+1]);
		else // 若为逆时针存放顶点
			ang-=angle(node[i],node[i+1]);
	if(dblcmp(ang)==0) // 角度若为0，则说明在多边形外部
		return false;
	if(dblcmp(ang-2*PI)==0 || dblcmp(ang+2*PI)==0) // 若为-360或360则说明在多边形内部
		return true;
	if(dblcmp(cirr)==0) //否则在多边形边上，判断圆半径是否为0，这里提前为主函数筛掉数据
		return true;
	else
		return false;
}
bool Is_right(){ //判断圆心到多边形所有边的距离是否都大于或等于圆半径
	for(int i=1;i<=n;i++){
		double h=fabs(cross(cir,node[i],cir,node[i+1]))/dist(node[i],node[i+1]);
		if(dblcmp(h-cirr)<0)
			return false;
	}
	return true;
}
int main(){
	while(scanf("%d",&n),n>=3){
		scanf("%lf%lf%lf",&cirr,&cir.x,&cir.y);
		for(int i=1;i<=n;i++)
			scanf("%lf%lf",&node[i].x,&node[i].y);
		node[0].x=node[n].x,node[0].y=node[n].y; // 为后续操作使用
		node[n+1].x=node[1].x,node[n+1].y=node[1].y;
		if(!Is_covbag()) //若不为凸包
			printf("HOLE IS ILL-FORMED\n");
		else{ 
			if(!Is_In() || !Is_right()) //若为凸包但圆不再凸包内
				printf("PEG WILL NOT FIT\n");
		    else //若圆在凸包内
				printf("PEG WILL FIT\n");
		}
	}
	return 0;
}