POJ-2318 TOYS,暴力+叉积判断！

                                                              TOYS

   2页的提交记录终于搞明白了。

   题意：一个盒子由n块挡板分成n+1块区域，有m个玩具，每次给出一个玩具的坐标，最后输出各个区域玩具的数量，玩具不会在线上或者盒子外面，挡板也不会相交。输出格式请看样例！

   思路：考察叉积的性质，用一个结构体存线段（两个点共4个坐标），然后输入一个玩具坐标就遍历所有的区域找到然后对应区域加1。

 这题看起来也不难，数据5000,开始想着用二分即可，但区域形状不规则，二分可能没有确定的区域，会有极端的例子如挡板全部是平行的斜线，倾斜度越高越不好找，于是放弃了这种思路，后来想着既然输入的区域顺序是确定的，那么对玩具坐标进行排序，每次只需从前面找的那个区域往后找即可。于是开启了无限WA模式，后来改成两层循环秒A。。。又重复用前面那种思路反复测试发现后台数据最大只有42个，两层循环随便水只要思路没错。

  先来一个WA的代码：

struct line
{
    double x1,y1,x2,y2;
} a[N],b[N];
int v[N];
int cmp(line a,line b)
{
    return a.x1<b.x1;
}
double multi(line b,line a)
{
    return (a.x1-b.x1)*(a.y2-b.y1)-(a.x2-b.x1)*(a.y1-b.y1);
}
int main()
{
    int n,m;
    double x1,y1,x2,y2;
    while(~scanf("%d",&n)&&n)
    {
        scanf("%d%lf%lf%lf%lf",&m,&x1,&y1,&x2,&y2);
        memset(v,0,sizeof(v));
        for(int i=1; i<=n; i++) //n条分界分成n+1个区域
        {
            scanf("%lf%lf",&a[i].x1,&a[i].x2);
            a[i].y1=y1;
            a[i].y2=y2;
        }
        a[0].x1=x1,a[0].y1=y1,a[0].x2=x1,a[0].y2=y2;//开始隔板
        a[n+1].x1=x2,a[n+1].y1=y1,a[n+1].x2=x2,a[n+1].y2=y2;//最后一块区域
        for(int i=0; i<m; i++) scanf("%lf%lf",&b[i].x1,&b[i].y1);
        sort(b,b+m,cmp);
        int j=0;
        for(int i=0; i<m; i++)
        {
            while(multi(b[i],a[j])>0&&j<=n+1) j++;//接着前一个区域往后找直到满足条件。
             j--;
            v[j]++;
            //j-=41;//就是这里测出了后台。
            //if(j<0) j=0;
        }
        for(int i=0; i<=n; i++)
            printf("%d: %d\n",i,v[i]);
        printf("\n");
    }
    return 0;
}

  以上思路看似很正确，但还是有退化的例子，比如：

2 2 0 100 100 0

1 50

8 70

25 99

40 1

应该输出：

0： 1

1： 1

2： 0

而运行结果：

0: 0
1: 1
2: 1

    

   正确代码：

struct line
{
    double x1,y1,x2,y2;
} a[N];
int v[N];
double multi(double x,double y,line a)
{
    return (a.x1-x)*(a.y2-y)-(a.x2-x)*(a.y1-y);
}
int main()
{
    int n,m;
    double x1,y1,x2,y2;
    while(~scanf("%d",&n)&&n)
    {
        scanf("%d%lf%lf%lf%lf",&m,&x1,&y1,&x2,&y2);
        memset(v,0,sizeof(v));
        for(int i=1; i<=n; i++) //n条分界分成n+1个区域
        {
            scanf("%lf%lf",&a[i].x1,&a[i].x2);
            a[i].y1=y1,a[i].y2=y2;
        }
        a[0].x1=x1,a[0].y1=y1,a[0].x2=x1,a[0].y2=y2;//开始隔板
        a[n+1].x1=x2,a[n+1].y1=y1,a[n+1].x2=x2,a[n+1].y2=y2;//最后一块区域
        double x,y;
        for(int i=0; i<m; i++)
        {
            int j=0;
            scanf("%lf%lf",&x,&y);
            while(multi(x,y,a[j])>0) j++;
            v[j-1]++;
        }
        for(int i=0; i<=n; i++)
            printf("%d: %d\n",i,v[i]);
        printf("\n");
    }
    return 0;
}