hdu 1756(判断点是否在多边形中)

传送门

题解：

　　射线法判定点是否在多边形内部；

AC代码：

　　

 #include<iostream>
 #include<cstdio>
 #include<cmath>
 using namespace std;
 const double esp=1e-;
 const int maxn=+;

 int n,m;
 struct Point
 {
     double x,y;
 }p[maxn];

 double K(Point p1,Point p2)
 {
     return (p2.y-p1.y)/(p2.x-p1.x);
 }
 bool onEdge(Point q,Point p1,Point p2)
 {
     if(p1.x == p2.x)//斜率不存在的情况，HDU的数据里没有卡这种情况
         return q.y >= min(p1.y,p2.y) && q.y <= max(p1.y,p2.y) && q.x == p1.x ? true:false;

     if(q.x >= min(p1.x,p2.x) && q.x <= max(p1.x,p2.x) && q.y >= min(p1.y,p2.y) && q.y <= max(p1.y,p2.y))
         return fabs(K(p1,q)-K(p1,p2)) < esp ? true:false;
     return false;
 }
 bool inPolygon(Point q)
 {
     int ans=;
     Point p1,p2;
     p1=p[];
     for(int i=;i <= n;++i)
     {
         p2=p[i%n];
         if(onEdge(q,p1,p2))
             return true;
         if(p1.y == p2.y)//判断p1p2是否为水平边，水平边不作考虑
         {
             p1=p2;//起初，这儿我并没有加这句话，改了好半天才看到
             continue;
         }
         /**
             注意：此处q.y是严格>min()；
             此处是用来处理射线与顶点相交的情况的
             当相交的顶点为凸顶点时，两条边都加上；
             当相交的顶点是凹顶点时，只加右边那条边；
         */
         if(q.y > min(p1.y,p2.y) && q.y <= max(p1.y,p2.y) && q.x <= max(p1.x,p2.x))
         {
             //当p1p2为竖线时，由条件q.x <= max(p1.x,p2.x)可得由p向右发出的射线与p1p2有交点
             if(p1.x == p2.x)
                 ans++;
             else
             {
                 //x : 经过q点且平行于x轴的直线与p1p2交点的横坐标
                 double x=p1.x+(q.y-p1.y)/K(p1,p2);
                 if(q.x <= x)//如果q.x <= x，说明q向右发出的射线与p1p2有交点
                     ans++;
             }
         }
         p1=p2;
     }
     return ans&;
 }
 int main()
 {
 //    freopen("C:\\Users\\hyacinthLJP\\Desktop\\in&&out\\HDU\\1756.in","r",stdin);
     while(~scanf("%d",&n))
     {
         for(int i=;i < n;++i)
             scanf("%lf%lf",&p[i].x,&p[i].y);

         scanf("%d",&m);
         for(int i=;i <= m;++i)
         {
             Point q;
             scanf("%lf%lf",&q.x,&q.y);
             if(inPolygon(q))
                 printf("Yes\n");
             else
                 printf("No\n");
         }
     }
     return ;
 }