POJ 2318 TOYS（计算几何）

Description
一个矩形区域中有n条互不相交的线段，每条线段的上下端点分别在矩形的上下边上，现给出该矩形区域中m个点的坐标，保证点不在边上，统计由这n条边将矩形分成的n+1个区域中各有多少个点
Input
多组用例，每组用例第一行六个整数n,m,xl,yl,xr,yr分别表示边数，点数，矩形区域左上角的横纵坐标，矩形区域右下角的横纵坐标，之后n行每行两个整数u和l表示该条边上下端点的横坐标，最后m行每行两个数x和y表示该点坐标
Output
对于每组用例，输出每个区域的点数，每两组输出之间用空行隔开
Sample Input
5 6 0 10 60 0
3 1
4 3
6 8
10 10
15 30
1 5
2 1
2 8
5 5
40 10
7 9
4 10 0 10 100 0
20 20
40 40
60 60
80 80
5 10
15 10
25 10
35 10
45 10
55 10
65 10
75 10
85 10
95 10
0
Sample Output
0: 2
1: 1
2: 1
3: 1
4: 0
5: 1

0: 2
1: 2
2: 2
3: 2
4: 2
Solution
对每个点求出其左侧最右边一条边即可推出其所在区域
Code

#include<cstdio>
#include<cstring>
#include<iostream>
#include<string>
#include<cmath>
#include<algorithm>
using namespace std;
struct node
{
    int x1,x2,y1,y2;//(x1,y1)表示该条线段的上端点坐标,(x2,y2)表示该条线段的下端点坐标 
}p[5555];
bool is_in(int x,int y,int a)//判断某点是否在某边的右侧 
{
    int x1=p[a].x1;
    int x2=p[a].x2;
    int y1=p[a].y1;
    int y2=p[a].y2;
    int ans=((x1-x)*(y2-y)-(y1-y)*(x2-x));
    return ans>0?true:false;
}
int main() 
{
    int n,m,xl,yl,xr,yr;
    int res=1;
    while(scanf("%d",&n),n)
    {
        scanf("%d%d%d%d%d",&m,&xl,&yl,&xr,&yr);
        if(res!=1)
            printf("\n");
        for(int i=1;i<=n;i++)
            scanf("%d%d",&p[i].x1,&p[i].x2);
        //将区间边界加入边集 
        p[0].x1=xl;
        p[0].x2=xl;
        p[n+1].x1=xr;
        p[n+1].x2=xr;
        for(int i=0;i<=n+1;i++)
        {
            p[i].y1=yl;
            p[i].y2=yr;
        }
        //for(int i=0;i<=n+1;i++)
        //  printf("%d %d %d %d\n",p[i].x1,p[i].y1,p[i].x2,p[i].y2);
        int cnt[5555]={0};
        while(m--)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            int temp=0;
            while(is_in(x,y,temp++));//找到在该点左侧的边最右一条边 
            temp-=2;
            cnt[temp]++;//该区域点数加一 
        }
        for(int i=0;i<=n;i++)
            printf("%d: %d\n",i,cnt[i]);
        res++;
    }
    return 0;
}