HDU1756(判断点在多边形内外)

最新推荐文章于 2022-06-25 09:55:29 发布

原创最新推荐文章于 2022-06-25 09:55:29 发布 · 1.6k 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

计算几何专栏收录该内容

40 篇文章

订阅专栏

本文介绍了一种用于判断点是否位于线段上、两线段是否相交及点是否位于多边形内部的算法实现。通过具体函数展示了如何进行几何运算，并提供了一个完整的点位置判断流程。

题目：http://acm.hdu.edu.cn/showproblem.php?pid=1756

bool IsOnline(Point p,Segment s)  /**  判断点p是否在线段s上  */
{
     if(fabs(cross(s.a,p,s.b))<=EPS&&min(s.a.x,s.b.x)<=p.x&&
        p.x<=max(s.a.x,s.b.x)&&min(s.a.y,s.b.y)<=p.y&&p.y<=max(s.a.y,s.b.y))
     return true;
     return false;
}

bool Segment_crossing(Segment u,Segment v)  /** 判断两线段是否相交  */
{
     return((max(u.a.x,u.b.x)>=min(v.a.x,v.b.x))&&
            (max(v.a.x,v.b.x)>=min(u.a.x,u.b.x))&&
            (max(u.a.y,u.b.y)>=min(v.a.y,v.b.y))&&
            (max(v.a.y,v.b.y)>=min(u.a.y,u.b.y))&&
            (cross(v.a,u.b,u.a)*cross(u.b,v.b,u.a)>=0)&&
            (cross(u.a,v.b,v.a)*cross(v.b,u.b,v.a)>=0));

}

bool Inside(Point t,Point p[],int n)     /*判断点t是否在点集p所围成的多边形内*/
{
     Segment seg;
     Segment pp;
     int num=0;
     Point MM;     /*定义无穷点*/
     MM.y=t.y;
     MM.x=INF;
     pp.a=t;
     pp.b=MM;           /*构造一条平行于x轴的以p为端点的右射线pp*/
     for(int i=0;i<n;i++)
     {
          seg.a=p[i];
          seg.b=p[(i+1)%n];
          if(IsOnline(t,seg))  return true;       /*点p在边上*/
          if(fabs(seg.a.y-seg.b.y)<EPS) continue;   /*水平*/
          if(IsOnline(seg.a,pp))          /*边的一端点在射线上*/
          {
               if(seg.a.y>seg.b.y) num++;
          }
          else if(IsOnline(seg.b,pp))
          {
               if(seg.b.y>seg.a.y) num++;
          }
          else if(Segment_crossing(pp,seg))
               num++;
     }
     if(num%2==1)
         return true;
     return false;
}