poj2398 Toy Storage 计算几何，叉积，二分

poj2398 Toy Storage

链接

poj

题目大意

这道题的大概意思是先输入6个数字：n,m,x1,y1,x2,y2。n代表卡片的数量，卡片竖直(或倾斜）放置在盒内，可把盒子分为n+1块区域，然后分别用0到n表示，m代表玩具的个数，（x1,y1）代表盒子的左上顶点坐标，(x2,y2)代表盒子的右下顶点坐标，接下来的n行，每行都有两个数a，b，(a,y1)，(b,y2)分别代表卡片的两个顶点位置，接下来的m行每行两个数从c,d，(c,d),代表玩具放置的坐标，最后让你输出当区域内玩具的个数（递增，0除外），以及有几个同样数目玩具的区域，当输入第一个数为0时结束。

思路

先把隔板排序。
然后每个玩具二分位置。用叉积判断。

代码

#include <algorithm>
#include <cstdio>
using namespace std;
const int N=1005;
const double eps=1e-10;
struct point {
    double x,y;
}a[N],b[N];
point operator - (point a,point b) {
    return (point){a.x-b.x,a.y-b.y};
}
double operator * (point a,point b) {
    return a.x*b.y-b.x*a.y;
}
struct line {
    double u,l;
}l[N];
int n,m;
bool cmp(line a,line b) {
    return a.u<b.u;
}
int cnt[N],ans[N];

int main() {
    while(scanf("%d",&n)!=EOF&&n) {
        scanf("%d",&m);
        scanf("%lf%lf%lf%lf",&b[0].x,&b[0].y,&a[n+1].x,&a[n+1].y);
        a[0].x=b[0].x,a[0].y=b[n+1].y;
        b[n+1].x=a[n+1].x,b[n+1].y=b[0].y;
        for(int i=1;i<=n;++i)
            scanf("%lf%lf",&l[i].u,&l[i].l);
        sort(l+1,l+1+n,cmp);
        for(int i=1;i<=n;++i) {
            a[i].x=l[i].l,a[i].y=a[n+1].y;
            b[i].x=l[i].u,b[i].y=b[0].y;
        }
        memset(cnt,0,sizeof(cnt));
        memset(ans,0,sizeof(ans));
        for(int i=1;i<=m;++i) {
            point c;
            scanf("%lf%lf",&c.x,&c.y);
            if(!(a[n+1].y<=c.y&&c.y<=b[0].y&&b[0].x<=c.x&&c.x<=a[n+1].x)) continue;
            int l=0,r=n+1;
            while(l+1!=r) {
                int mid=(l+r)>>1;
                if((b[mid]-a[mid])*(c-a[mid])>0) r=mid;
                else l=mid;
            }
            cnt[l+1]++;
        }
        puts("Box");
        for(int i=1;i<=n+1;++i)
            if(cnt[i]) ans[cnt[i]]++;
        for(int i=1;i<=m;++i)
            if(ans[i]) printf("%d: %d\n",i,ans[i]);
    }
    return 0;
}