poj 2398 计算几何

#include <iostream>
#include<cstdio>
#include<cstring>
#include <algorithm>
using namespace std;

struct point{
    int x,y;
};

struct line{
    point a,b;
};
line lline[5005];

int cnt[5005];
int res[5005];
int cross(point p1,point p2,point p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x);
}

bool cmp(line c,line d){
    return c.b.x < d.b.x ;
}

int main()
{
    int n,m,x1,x2,y1,y2;
    int t1,t2;
    point c;
    while(cin>>n&&n){
        cin>>m>>x1>>y1>>x2>>y2;
        lline[0].a.x=x1,lline[0].a.y=y1;
        lline[0].b.x=x1,lline[0].b.y=y2;
        lline[n+1].a.x = x2,lline[n+1].a.y=y1;
        lline[n+1].b.x = x2,lline[n+1].b.y=y2;
        for(int i=1;i<=n;i++){
            scanf("%d%d",&t1,&t2);
            lline[i].a.x=t1;
            lline[i].a.y=y1;
            lline[i].b.x=t2;
            lline[i].b.y=y2;
        }

        sort(lline,lline+n+2,cmp);
       // for(int i=0;i<n+2;i++)
       //     cout<<lline[i].b.x<<" "<<lline[i].b.y<<" "<<lline[i].a.x<<" "<<lline[i].a.y<<endl;
        memset(cnt,0,sizeof(cnt));
        for(int i=0;i<m;i++){
            cin>>c.x>>c.y;
            for(int j=0;j<n+2;j++){
                //cout<<j<<" "<<cross(c,lline[j].a,lline[j].b)<<endl;
                if(cross(c,lline[j].a,lline[j].b)>0){
                    continue;
                }else{
                    cnt[j-1]++;
                    break;
                }
            }
        }
       // for(int i=0;i<=n;i++)
       //     cout<<cnt[i]<<endl;
        memset(res,0,sizeof(res));
        for(int i=0;i<=n;i++){
            res[cnt[i]]++;
        }
        cout<<"Box"<<endl;
        for(int i=1;i<=n;i++)if(res[i])
            cout<<i<<": "<<res[i]<<endl;
    }
    return 0;
}