poj 2318 TOY题解

二分+叉积

如图所示:N个隔板将盒子分为N+1个区域，若toy在盒子隔板的左边则叉积小于0，在右边则叉积大于0。

为此，则可以用二分查找确定该点位于哪个区域（当l==r时，根据在隔板的那一侧，便可确定在哪个区域了）。

#include<stdio.h>
#include<string.h>
#include<algorithm>
#define MAX 5005
using namespace std;
int N,M;
double x1,y1,x2,y2;
typedef struct{
    double x;
    double y;
}P;
P toy[MAX];
struct Node{
    P first;
    P second;
}Pline[MAX];
int pnum[MAX];
double mul(P p1,P p2,P p0){
    return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

void bin(P tmp){
    int l=0,r=N-1;
    while(l<r){
        int mid=(l+r)/2;
        if(mul(tmp,Pline[mid].first,Pline[mid].second)>0){
            l=mid+1;
        }else{
            r=mid;
        }
    }
    if(mul(tmp,Pline[l].first,Pline[l].second)>0) pnum[l+1]++;
    else pnum[l]++;
}
int main(){
    while(~scanf("%d",&N),N){
        memset(pnum,0,sizeof(pnum));
        scanf("%d %lf %lf %lf %lf",&M,&x1,&y1,&x2,&y2);
        for(int i=0;i<N;i++){
            scanf("%lf %lf",&Pline[i].first.x,&Pline[i].second.x);
            Pline[i].first.y=y1;
            Pline[i].second.y=y2;
        }
        for(int i=0;i<M;i++){
            scanf("%lf %lf",&toy[i].x,&toy[i].y);
            bin(toy[i]);
        }
        for(int i=0;i<=N;i++){
            printf("%d: %d\n",i,pnum[i]);
        }
        puts("");
    }
return 0;
}