hrbustoj 1306:再遇攻击（计算几何，判断点是否在多边形内，水题）

再遇攻击

Time Limit: 1000 MS    Memory Limit: 65536 K

Total Submit: 253(37 users)   Total Accepted: 56(29 users)       Rating:         Special Judge: No

Description

Dota中英雄技能攻击会有一个范围，现在释放一个技能给出他的攻击范围和目标英雄的位置，问是否能攻击到。攻击范围保证是一个多边型。

Input

有多组测试数据

第一行输入1个整数n, 期中n代表攻击范围是给出的n个点组成的多边形,按照时针方向（顺或逆）依次给出，（n>=3）N = 0输入结束。

第二行a,b表示目标英雄的坐标( 0 < a,b<100)

接下来有n行，每行两个整数x，y（0 < x,y <100）表示每个点的坐标

攻击范围在边缘也算在内

Output

每组结果输出占一行

如果能够攻击到输出”Yes”

否则输出”No”

Sample Input

3

1 1

4 4

5 4

4 6

0

Sample Output

No

Author

鲁学涛

---

　　计算几何，判断点是否在多边形内。

　　水题，没有什么需要注意的地方。

　　代码：

 #include <stdio.h>
 struct Point{
     double x,y;
 };
 struct Line{
     Point p1,p2;
 };
 double xmulti(Point p1,Point p2,Point p0)    //求p1p0和p2p0的叉积,如果大于0,则p1在p2的顺时针方向
 {
     return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
 }
 double Max(double a,double b)
 {
     return a>b?a:b;
 }
 double Min(double a,double b)
 {
     return a<b?a:b;
 }
 bool ponls(Point q,Line l)    //判断点q是否在线段l上
 {
     if(q.x > Max(l.p1.x,l.p2.x) || q.x < Min(l.p1.x,l.p2.x)
         || q.y > Max(l.p1.y,l.p2.y) || q.y < Min(l.p1.y,l.p2.y) )
         return false;
     if(xmulti(l.p1,l.p2,q)==)    //点q不在l的延长线或者反向延长线上，如果叉积再为0，则确定点q在线段l上
         return true;
     else
         return false;
 }
 bool pinplg(int pointnum,Point p[],Point q)
 {
     Line s;
     int c = ;
     for(int i=;i<=pointnum;i++){    //多边形的每条边s
         if(i==pointnum)
             s.p1 = p[pointnum],s.p2 = p[];
         else
             s.p1 = p[i],s.p2 = p[i+];
         if(ponls(q,s))    //点q在边s上
             return true;
         if(s.p1.y != s.p2.y){    //s不是水平的
             Point t;
             t.x = q.x - ,t.y = q.y;
             if( (s.p1.y == q.y && s.p1.x <=q.x) || (s.p2.y == q.y && s.p2.x <= q.x) ){    //s的一个端点在L上
                 int tt;
                 if(s.p1.y == q.y)
                     tt = ;
                 else if(s.p2.y == q.y)
                     tt = ;
                 int maxx;
                 if(s.p1.y > s.p2.y)
                     maxx = ;
                 else
                     maxx = ;
                 if(tt == maxx) //如果这个端点的纵坐标较大的那个端点
                     c++;
             }
             else if(xmulti(s.p1,t,q)*xmulti(s.p2,t,q) <= ){    //L和边s相交
                 Point lowp,higp;
                 if(s.p1.y > s.p2.y)
                     lowp.x = s.p2.x,lowp.y = s.p2.y,higp.x = s.p1.x,higp.y = s.p1.y;
                 else
                     lowp.x = s.p1.x,lowp.y = s.p1.y,higp.x = s.p2.x,higp.y = s.p2.y;
                 if(xmulti(q,higp,lowp)>=)
                     c++;
             }
         }
     }
     if(c%==)
         return false;
     else
         return true;
 }
 int main()
 {
     int i,n;
     while(scanf("%d",&n)!=EOF){
         if(n==) break;
         Point t;
         scanf("%lf%lf",&t.x,&t.y);
         Point p[];
         for(i=;i<=n;i++)    //输入多边形的顶点
         scanf("%lf%lf",&p[i].x,&p[i].y);
         if(pinplg(n,p,t))
             printf("Yes\n");
         else
             printf("No\n");
     }
     return ;
 }

Freecode : www.cnblogs.com/yym2013