[HDU1756]Cupid's Arrow

题目大意：
　　给你一个简单多边形和若干个点，问每个点在多边形内还是外。

思路：
　　一开始没看清楚题，写了一个叉积法，事实上叉积法只能用来处理凸多边形与点的关系。
　　考虑一个射线法。
　　从这个点水平往左作一条射线，如果与多边形相交次数为偶数，则在多边形外，否则在多边形内。
　　考虑一些特殊情况：
　　　　1.射线与多边形某条边重合，不算。
　　　　2.射线经过某个顶点，只能够算一次。
　　对于第一种情况，我们只需要判一下是否平行即可，对于第二种情况，我们对于多边形上的边区分一下上顶点和下顶点即可，上顶点算，下顶点不算。

 #include<cstdio>
 #include<cctype>
 #include<algorithm>
 inline int getint() {
     register char ch;
     while(!isdigit(ch=getchar()));
     register int x=ch^'';
     while(isdigit(ch=getchar())) x=(((x<<)+x)<<)+(ch^'');
     return x;
 }
 const int N=;
 struct Point {
     double x,y;
 };
 Point p[N];
 int main() {
     int n;
     while(~scanf("%d",&n)) {
         for(register int i=;i<n;i++) {
             p[i]=(Point){getint(),getint()};
         }
         for(register int m=getint();m;m--) {
             const Point q=(Point){getint(),getint()};
             int cnt=;
             if(p[].y==p[n-].y) {
             } else if(p[].x==p[n-].x) {
                 if(p[].x<q.x&&std::min(p[].y,p[n-].y)<q.y&&q.y<=std::max(p[].y,p[n-].y)) {
                     cnt++;
                 }
             } else if((q.y-p[n-].y)/((p[].y-p[n-].y)/(p[].x-p[n-].x))+p[n-].x<=q.x) {
                 if(std::min(p[].y,p[n-].y)<q.y&&q.y<=std::max(p[].y,p[n-].y)) {
                     cnt++;
                 }
             }
             for(register int i=;i<n;i++) {
                 if(p[i].y==p[i-].y) {
                 } else if(p[i].x==p[i-].x) {
                     if(p[i].x<q.x&&std::min(p[i].y,p[i-].y)<q.y&&q.y<=std::max(p[i].y,p[i-].y)) {
                         cnt++;
                     }
                 } else if((q.y-p[i-].y)/((p[i].y-p[i-].y)/(p[i].x-p[i-].x))+p[i-].x<=q.x) {
                     if(std::min(p[i].y,p[i-].y)<q.y&&q.y<=std::max(p[i].y,p[i-].y)) {
                         cnt++;
                     }
                 }
             }
             puts((cnt&)?"Yes":"No");
         }
     }
     return ;
 }