zoj 1081:Points Within（计算几何，判断点是否在多边形内，经典题）

Points Within

---

Time Limit: 2 Seconds      Memory Limit: 65536 KB

---

Statement of the Problem

Several drawing applications allow us to draw polygons and almost all of them allow us to fill them with some color. The task of filling a polygon reduces to knowing which points are inside it, so programmers have to colour only those points.

You're expected to write a program which tells us if a given point lies inside a given polygon described by the coordinates of its vertices. You can assume that if a point is in the border of the polygon, then it is in fact inside the polygon.

Input Format

The input file may contain several instances of the problem. Each instance consists of: (i) one line containing integers N, 0 < N < 100 and M, respectively the number of vertices of the polygon and the number of points to be tested. (ii) N lines, each containing a pair of integers describing the coordinates of the polygon's vertices; (iii) M lines, each containing a pair of integer coordinates of the points which will be tested for "withinness" in the polygon.

You may assume that: the vertices are all distinct; consecutive vertices in the input are adjacent in the polygon; the last vertex is adjacent to the first one; and the resulting polygon is simple, that is, every vertex is incident with exactly two edges and two edges only intersect at their common endpoint. The last instance is followed by a line with a 0 (zero).

Output Format

For the ith instance in the input, you have to write one line in the output with the phrase "Problem i:", followed by several lines, one for each point tested, in the order they appear in the input. Each of these lines should read "Within" or "Outside", depending on the outcome of the test. The output of two consecutive instances should be separated by a blank line.

Sample Input

3 1
0 0
0 5
5 0
10 2
3 2
4 4
3 1
1 2
1 3
2 2
0

Sample Output

Problem 1:
Outside

Problem 2:
Outside
Within

---

　　

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

　　套用了以前写的模板，很easy的就过了。经典题。

　　代码：

 #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,m,num=;
     while(scanf("%d",&n)!=EOF){
         if(n==) break;
         scanf("%d",&m);
         Point p[];
         for(i=;i<=n;i++)    //输入多边形的顶点
             scanf("%lf%lf",&p[i].x,&p[i].y);
         if(num!=)    //如果不是第一组数据块就先输出一个空行
             printf("\n");
         printf("Problem %d:\n",num++);
         for(i=;i<=m;i++){
             Point t;
             scanf("%lf%lf",&t.x,&t.y);
             if(pinplg(n,p,t))
                 printf("Within\n");
             else
                 printf("Outside\n");
         }
     }
     return ;
 }

Freecode : www.cnblogs.com/yym2013