poj2318

计算几何的解题报告开始了！

这是第一道计算几何题，题目不难，但是题水人更水。

 

题意：有一个大箱子，由n个板分为n+1块，标号为0~n

         已知盒子左上角和右下角的坐标及每个板上下两端的横坐标（板不会交错，且按顺序给出）

         然后给出玩具的坐标，统计每块空间内玩具个数（保证玩具一定落在空间内）

方法：叉积+二分，判断点所在的空间

代码：#include <iostream> #include <cstring> #include <algorithm> #include <cstdio> using namespace std; struct Point { long long x,y;//不用long long 过不了 }; int n,m;//n个板，m个玩具 Point a,b;//左上点和右下点 long long up[5010];//板上端横坐标 long long down[5010];//板下端横坐标 Point toy[5010];//玩具坐标 long long res[5050];//存放结果 long long Cross(Point & o, Point &a, Point& b)//判断向量oa 与ob的位置 { return (a.x-o.x)*(b.y-o.y)- (a.y-o.y) * (b.x-o.x); } int BinarySearch(Point key)//返回 该点所在空间编号+1 { int start = 0; int end = n; int mid; Point pup,pdown; pup.y = a.y; pdown.y = b.y; while(start <=end) { mid = (start + end)/ 2; pup.x = up[mid]; pdown.x = down[mid]; if(Cross(pdown,key,pup)<0) end = mid - 1; else start = mid + 1; } return start; } void Input() { cin >> a.x >> a.y; cin >> b.x >> b.y; up[0] = down[0] = a.x;//应该初始化为a.x原来初始化为0了 //up[n+1] = down[n+1] = b.x;//这一句写不写都一样，二分返回应该插入的位置 for(int i = 1; i <=n; i++) { cin >> up[i] >> down[i]; } for(int i = 1; i <=m; i++) cin >> toy[i].x >> toy[i].y; } void Compute() { int index; memset(res,0,sizeof(res)); for(int i = 1; i <=m; i++) { index = BinarySearch(toy[i]); //cout << "index = " << index-1 << endl; res[index-1]++; } } void Output() { for(int i = 0; i <=n; i++) { cout << i << ": " << res[i] << endl; } } int main() { int casecnt = 1; while(cin >> n ) { if(n == 0) return 0; cin >> m; Input(); Compute(); if(casecnt!=1) cout << endl; casecnt++; Output(); } return 0; }