计算几何（一）：叉积的简单应用

 这是学习《ACM-ICPC程序设计系列—计算几何》自己AC的第一个计算几何的问题。题目是比较简单的，但还是花了我很久的时间。

    题目可以抽象成：一个长方形被n条不相交的线段分隔成n+1个区间，给定m个点的坐标，计算出每个区间里各有多少个点。

    需要注意的是：输入时，分隔长方形的线段已经排序。（这意味着可以用二分查找）

       我的思路很简单：把区间构造起来，（这里我把长方形比喻为盒子，中间加入的线段我称作隔板。）然后，每个点去二分搜索，找到自己处在的区间就OK了。最后输出结果。

    这里面的难点就在于：我们怎么去判断一个点在一条线段的左侧还是右侧呢？这就要用到向量叉积。叉积的一个非常重要的性质是通过它的符号判断两向量相互之间的顺逆时针关系：设向量P=(x₁,y₁),Q=(x₂,y₂)

如果P*Q>0则P在Q的顺时针方向；

如果P*Q=0则P与Q共线，可能同向，与可能反向；

如果P*Q<0则P在Q的逆时针方向。

  （请注意这里是叉乘）

    这道题在几何上就考了这一个知识点，如果我们一步一步来，理清关系，其实是很简单的！

代码：

1. /* 
2.  * poj_2318.cpp 
3.  * 
4.  *  Created on: 规定上面是起点，下面是终点；左面是起点，右面是终点 
5.  *      Author: Administrator 
6.  */ 
7. #include<stdio.h> 
8. #include<string.h> 
9. #define M 5005 
10. struct Point{ 
11.     double x,y; 
12. }; 
13. struct Segment{ 
14.     Point start,end; 
15. }; 
16.  
17. int bin[M]; 
18. Segment box[M];//表示盒子 
19.  
20. int n,m; 
21. double x1,y1,x2,y2;//表示盒子最外面的边框 
22.  
23. /* 
24.  *点p置于线段s的哪一侧：0 左侧，1 右侧 
25.  * */ 
26. int isLeft(const Point &p,const Segment &s){ 
27.     double ans; 
28.     ans=(p.x-s.end.x)*(s.start.y-s.end.y)-(s.start.x-s.end.x)*(p.y-s.end.y); 
29.     if(ans<0)return 1; 
30.     else if(ans>0)return 0; 
31.     else{printf("ERROR!\n"); 
32.         return -1;} 
33. } 
34. //检查这个点是否在盒子里面 
35. bool inside(Point p){ 
36.     if(p.x>x2||p.x<x1||p.y<y2||p.y>y1) 
37.         return false; 
38.     return true; 
39. } 
40. //盒子中加入隔板 
41. void add(const double &x1,const double &y1,const double &x2,const double &y2, int i){ 
42.     box[i].start.x=x1;box[i].start.y=y1; 
43.     box[i].end.x=x2;box[i].end.y=y2; 
44. } 
45. //查找玩具对应的格子：因为已经排序，所以可以用二分查找 
46. int solve(int left,int right,const Point &p){ 
47.     int mid=(left+right)>>1; 
48.     if(left==right){ 
49.         return isLeft(p,box[mid])==1? mid: -1; 
50.     } 
51.     if(isLeft(p,box[mid])==1){//在mid的左边 
52.         return solve(left,mid,p); 
53.     }else{ 
54.         return solve(mid+1,right,p); 
55.     } 
56. } 
57. int main(){ 
58.     int id; 
59.     Point toy; 
60.     while(1){ 
61.         scanf("%d",&n); 
62.         if(n==0)break; 
63.         scanf("%d",&m); 
64.         scanf("%lf %lf %lf %lf",&x1,&y1,&x2,&y2); 
65.         add(x1,y1,x1,y2,0); 
66.         add(x2,y1,x2,y2,n+1); 
67.         double ui,li; 
68.         for(int i=1;i<=n;i++){ 
69.             scanf("%lf %lf",&ui,&li); 
70.             add(ui,y1,li,y2,i);//盒子里添加隔板 
71.         } 
72.  
73.         memset(bin,0,(n+1)*sizeof(int)); 
74.         for(int i=0;i<m;i++){ 
75.             scanf("%lf %lf",&toy.x,&toy.y); 
76.             if(inside(toy)){//如果玩具确实扔到盒子里了，就计算玩具所在的格子 
77.                 id=solve(1,n+1,toy); 
78.                 if(id>0)bin[id-1]++; 
79.                 else printf("ERROR!\n"); 
80.             } 
81.         } 
82.         //print result 
83.         for(int i=0;i<=n;i++){ 
84.             printf("%d: %d\n",i,bin[i]); 
85.         } 
86.         printf("\n"); 
87.     } 
88.     return 0; 
89. }