传送门
解析:
射线法裸题,其实也可以用角度判别法过去。
思路:
判断一个点是否在多边形内部,我们从这个点向任意方向随便引一条射线,如果这条射线与多边形的交点有奇数个,那么这个点就是在多边形内部的。
简单口胡 证明:
考虑一条直线穿过一个多边形,它一定与多边形有偶数个交点,那么我们从第一个交点向后走,显然这段是在多边形内部的,然后经过一个交点后我们会在这个多边形外部。
因为多边形的边就是他内部和外部的分界线,自然这种性质就体现在了交点上。
那么当经过奇数个交点的时候,在多边形内部。
但是显然会有很多特殊情况。
比如射线交在了多边形的顶点上。。。
射线与多边形的边重合。。。等等例子
其实边重合很好理解,直接特判点是否在边上,如果在,根据题目要求直接返回答案(在多边形内是否包含多边形的边)。如果不在,那我们无视这条边与射线的交点,不会影响答案。
如果是顶点重合,还要分情况,如果与重合顶点相邻的两边在射线同侧,显然穿过这个顶点之后仍然在多边形内,这个顶点就不能计算入交点。
反之,则需要计算,为保证只计算一次,我们规定与边的某一个端点交才算交,与另一个端点交则不算。
虽然讨论的情况很多,但是最终化归出的判断方法非常简单,具体请看代码中结构体PolygonPolygonPolygon的成员函数innerinnerinner。
代码:
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const
inline int getint(){
re int num;
re char c;
re bool f=0;
while(!isdigit(c=gc()))if(c=='-')f=1;num=c^48;
while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
return f?-num:num;
}
struct Point{
int x,y;
Point(cs int &_x=0,cs int &_y=0):x(_x),y(_y){}
friend Point operator+(cs Point &a,cs Point &b){return Point(a.x+b.x,a.y+b.y);}
friend Point operator-(cs Point &a,cs Point &b){return Point(a.x-b.x,a.y-b.y);}
friend int operator*(cs Point &a,cs Point &b){return a.x*b.y-a.y*b.x;}
friend int dot(cs Point &a,cs Point &b){return a.x*b.x+a.y*b.y;}
};
inline bool onseg(cs Point &u,cs Point &v,cs Point &p){
if((u-p)*(v-p)!=0)return false;
return dot(u-p,v-p)<=0;
}
struct Polygon{
int n;
Point p[102];
void get(){for(int re i=1;i<=n;++i)p[i].x=getint(),p[i].y=getint();}
void clear(){memset(p,0,sizeof p);}
bool inner(cs Point &s)cs{
int cnt=0;
for(int re i=1;i<=n;++i){
if(onseg(p[i],p[i%n+1],s))return true;
int d1=p[i].y-s.y,d2=p[i%n+1].y-s.y;
int det=(p[i]-s)*(p[i%n+1]-s);
if((det>=0&&d1<0&&d2>=0)||(det<=0&&d1>=0&&d2<0))++cnt;
}
return cnt&1;
}
}P;
int m;
signed main(){
for(int re tt=1;;++tt){
P.n=getint();
if(P.n==0)return 0;
m=getint();
P.clear();
P.get();
if(tt!=1)pc('\n');
printf("Problem %d:\n",tt);
while(m--){
Point s;
s.x=getint();
s.y=getint();
if(P.inner(s))puts("Within");
else puts("Outside");
}
}
return 0;
}