题意和做法转载自
優YoU http://user.qzone.qq.com/289065406/blog/1309142308
大致题意:
按照顺时针或逆时针方向输入一个n边形的顶点坐标集,先判断这个n边形是否为凸包。
再给定一个圆形(圆心坐标和半径),判断这个圆是否完全在n变形内部。
解题思路:
题意已经很直白了。。就是那个思路。。。
注意输入完顶点集后,要封闭多边形,方便后面枚举边。
封闭方法:
定义点集数组Vectex[1~n]记录n个顶点,再令Vectex[0]=Vectex[n],Vectex[n+1]=Vectex[1]
1、判断凸包:
由于点集已经按某个时针方向有序,因此可以先定义一个方向系数direction=0
两两枚举n边形的边,用叉积判断这两条边的转向(右螺旋或左螺旋),由于存在散点共线的情况,因此当且仅当叉积的值temp第一次不为0时,direction=temp,direction的值此后不再改变。(direction>0 则为右螺旋逆时针,direction<0则为左螺旋顺时针)
此后继续枚举剩下的边,只要判断direction*temp>=0即可,当存在一个direction*temp<0的边,说明这是凹多边形,就不是凸包了。
2、判断圆心与多边形的关系:
用环顾法:
设圆心为P,逐条枚举n边形的边AB,利用
计算PA和PB的夹角,最后求和得到的就是环顾角。
(1) 圆心在多边形内部时,环顾角=±360
(2) 圆心在多边形外部时,环顾角=0
(3) 圆心在多边形边上时(不包括顶点),环顾角=±180
(4) 圆心在多边形顶点时,环顾角为(0,360)之间的任意角,其实就是圆心所在的顶点的两条邻接边的夹角。
3、当圆心在圆内时,判断圆与多边形的关系
设圆心为P,逐条枚举n边形的边AB,利用得到△PAB的面积,
再根据公式S=0.5*|AB|*h,可以得到
枚举所有h与圆的半径R比对,只要所有的边都有R-h>=0,则说明圆在多边形内
#include <cstdio>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <map>
#include <stack>
#include <vector>
#define maxn 10010
#define maxe 100010
typedef long long ll;
using namespace std;
const double eps=1e-5;
const int inf=0x3f3f3f3f3f;
const double PI=2*asin(1);
typedef double T1;
struct Point
{
T1 x,y;
Point(){};
bool discard;
Point(T1 a,T1 b)
{
x=a,y=b;
discard=false;
}
void input()
{
scanf("%lf%lf",&x,&y);
discard=false;
}
Point operator +(Point a)
{
Point b(x+a.x,y+a.y);
return b;
}
Point operator -(Point a)
{
Point b(x-a.x,y-a.y);
return b;
}
T1 operator *(Point a)
{
return x*a.x+y*a.y;
}
T1 operator ^(Point a)
{
return x*a.y-y*a.x;
}
bool operator <(Point a)
{
return x<a.x;
}
}o;
Point p[maxn];
Point ans[maxn];
double xmult(Point p0,Point p1,Point p2)
{
return (p1-p0)^(p2-p0);
}
int sgn(double x)
{
if(fabs(x) < eps)return 0;
if(x < 0) return -1;
return 1;
}
double dist(Point a,Point b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
//判断是不是凸多边形,输入的是逆时针或顺时针排好序的
bool judge(int n)
{
double dir=0;
double tmp;
for(int i=0;i<n;i++)
{
int j=(i+1)%n;
int k=(i-1+n)%n;
tmp=xmult(p[i],p[k],p[j]);
if(sgn(tmp)&&!sgn(dir))
{
dir=tmp;
}
if(sgn(dir)*sgn(tmp)<0)return false;
}
return true;
}
//判断点和多边形的关系
bool judge1(Point p1,int n)
{
double tmp=0.0;
int j;
for(int i=0;i<n;i++)
{
j=(i+1)%n;
double tmp1=dist(p1,p[i])*dist(p1,p[j]);
if(sgn(tmp1)==0)return false;
tmp1=(p1-p[i])*(p1-p[j])/tmp1;
//cout<<tmp1<<endl;
tmp=acos(tmp1)+tmp;
//cout<<"tmp "<<tmp<<endl;
}
//cout<<tmp<<endl;
if(sgn(tmp-2*PI)==0||sgn(tmp+2*PI)==0)return true;
return false;
}
//点到线段的距离
double Dis_Line_Point(Point p0,Point a,Point b)
{
if(sgn(xmult(p0,a,b))==0)
{
double x1=max(a.x,b.x);
double x2=min(a.x,b.x);
double y1=max(a.y,b.y);
double y2=min(a.y,b.y);
if(p0.x<=x1&&p0.x>=x2&&p0.y<=y1&&p0.y>=y2)
{
return 0;
}
else
{
return min(dist(p0,a),dist(p0,b));
}
}
else
{
return fabs(xmult(p0,a,b)/dist(a,b));
}
}
int main()
{
int n;
double r;
Point round;
//freopen("in.txt","r",stdin );
while(scanf("%d",&n)==1)
{
if(n<3)break;
scanf("%lf",&r);
round.input();
for(int i=0;i<n;i++)
{
p[i].input();
}
if(!judge(n))puts("HOLE IS ILL-FORMED");
//判断
else
{
if(!judge1(round,n))
{
puts("PEG WILL NOT FIT");
}
else
{
int flag=1;
for(int i=0;i<n;i++)
{
int j=(i+1)%n;
if(Dis_Line_Point(round,p[i],p[j])<r)
{
puts("PEG WILL NOT FIT");
flag=0;
break;
}
}
if(flag)puts("PEG WILL FIT");
}
}
}
return 0;
}