POJ2318计算几何

POJ2318

题意：一个矩形被分为n+1块，每个玩具都有坐标，求玩具在哪个块里面；

我一开始，没完全理解叉积，把叉积的方向弄错了；

 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct Point{
    double x,y;
}point[5005];
struct Line{
    Point a,b;
}line[5005];
double cross(Point p0,Point p1,Point p2) //p0在p1p2的什么方向，而不是p0在p2p1的方向，注意
{
    return (p0.x-p1.x)*(p2.y-p1.y)-(p2.x-p1.x)*(p0.y-p1.y);
}
int box[5005];
int main()
{
    int n,m,cc=0;
    double x1,y1,x2,y2;
    while( scanf("%d",&n)&&n )
    {
        if( cc )
           cout<<endl;
        cc=1;
        memset(box,0,sizeof(box));
        scanf("%d%lf%lf%lf%lf",&m,&x1,&y1,&x2,&y2);
        for(int i=1;i<=n;i++)
        {
            double U,L;
            scanf("%lf%lf",&U,&L);
            line[i].a.x=U;line[i].a.y=y1;
            line[i].b.x=L;line[i].b.y=y2;
        }
        line[0].a.x=x1;line[0].a.y=y1;
        line[0].b.x=x1;line[0].b.y=y2;
        line[n+1].a.x=x2;line[n+1].a.y=y1;
        line[n+1].b.x=x2;line[n+1].b.y=y2;
        for(int i=1;i<=m;i++)
        {
            double x,y;
            scanf("%lf%lf",&x,&y);
            point[i].x=x;
            point[i].y=y;
        }
        for(int i=1;i<=m;i++)
        {
            int l=0,r=n+1,flag=0;
            while( l!=r&&(l+1)!=r )
            {
                int mid=(l+r)/2;
                double k=cross(point[i],line[mid].b,line[mid].a); //叉积方向不能错了；
                if( k<0 )
                    r=mid;
                else
                    l=mid;
                //printf("%d %d %d %.2lf\n",l,r,mid,k);

            }
            box[l]++;
        }
        for(int i=0;i<=n;i++)
        {
            printf("%d: %d\n",i,box[i]);
        }
    }
    return 0;
}