poj 2318 poj 2398

             题目意思，是给了你一个盒子，其中有板将其分成多个格子，然后有多个点，肯定是在这些格子中（在边上也算进格子）的，问各个格子拥有的点的个数。

       叉积+二分搜索。由于点必定在格子中，因此，用目标点和格子边上的四个点进行叉积运算，点在格子的“左边”的右边且在格子的“右边”的左边时，即在该格子中，否则，根据叉积结果二分搜索。

       叉积运算的知识可以看黑书或者是算导，有很详细的讲解。

 

      

#include<cstdio>
using namespace std;
const int M=50200;
const double eps=1e-6;
struct point
{
	int x,y;
}p;
struct rec
{
	point p[4];
}r[M];
int cnt[M];
int n,m;

void init()
{
	for(int i=0;i<=n+1;i++)
		cnt[i]=0;
}

double cross(double x1,double y1,double x2,double y2)
{
	return x1*y2-x2*y1;
}

void search()
{
	int low=0,high=n;
	bool flag=false;
	while(low<=high)
	{
		int mid=(low+high)>>1;
		double temp1=cross(r[mid].p[0].x-p.x,r[mid].p[0].y-p.y,r[mid].p[1].x-p.x,r[mid].p[1].y-p.y);
		double temp2=cross(r[mid].p[3].x-p.x,r[mid].p[3].y-p.y,r[mid].p[2].x-p.x,r[mid].p[2].y-p.y);
		if(temp1>=eps && temp2<=-eps)
		{
			cnt[mid]++;return ;
		}
		if(temp2>eps)
			low=mid+1;
		else
			high=mid-1;
	}
}

int main()
{
	while(scanf("%d",&n)==1)
	{
		if(!n) break;
		init();
		scanf("%d",&m);
		scanf("%d%d%d%d",&r[0].p[0].x,&r[0].p[0].y,&r[n].p[2].x,&r[n].p[2].y);
		r[0].p[1].x=r[0].p[0].x;
		r[0].p[1].y=r[n].p[2].y;
		r[n].p[3].x=r[n].p[2].x;
		r[n].p[3].y=r[0].p[0].y;
		for(int i=1;i<=n;i++)
		{
			scanf("%d%d",&r[i].p[0].x,&r[i].p[1].x);
			r[i].p[0].y=r[i-1].p[0].y;
			r[i].p[1].y=r[i-1].p[1].y;
			r[i-1].p[3].x=r[i].p[0].x;
			r[i-1].p[3].y=r[i].p[0].y;
			r[i-1].p[2].x=r[i].p[1].x;
			r[i-1].p[2].y=r[i].p[1].y;
		}
		for(int i=0;i<m;i++)
		{
			scanf("%d%d",&p.x,&p.y);
			search();
		}
		for(int i=0;i<=n;i++)
		{
			printf("%d: %d\n",i,cnt[i]);
		}
		printf("\n");
	}
	return 0;
}

           

             而poj 2398只是加上一个排序（因为板的坐标不是有序的），再加多一个统计就好了。

 

             

#include<cstdio>
#include<algorithm>
using namespace std;
const int M=1030;
const double eps=1e-8;
struct point
{
	__int64 x,y;
}p[M],q;
struct rec
{
	point p[4];
}r[M];
__int64 cnt[M],cnt1[M];
int n,m;

void init()
{
	for(int i=0;i<=n+1;i++)
		cnt[i]=cnt1[i]=0;
}

double cross(double x1,double y1,double x2,double y2)
{
	return x1*y2-x2*y1;
}

void search()
{
	int low=0,high=n;
	bool flag=false;
	while(low<=high)
	{
		int mid=(low+high)>>1;
		double temp1=cross(r[mid].p[0].x-q.x,r[mid].p[0].y-q.y,r[mid].p[1].x-q.x,r[mid].p[1].y-q.y);
		double temp2=cross(r[mid].p[3].x-q.x,r[mid].p[3].y-q.y,r[mid].p[2].x-q.x,r[mid].p[2].y-q.y);
		if(temp1>=eps && temp2<=-eps)
		{
			cnt[mid]++;return ;
		}
		if(temp2>eps)
			low=mid+1;
		else
			high=mid-1;
	}
}

bool cmp(point p,point q)
{
	return p.x<q.x;
}

int main()
{
	while(scanf("%d",&n)==1)
	{
		if(!n) break;
		init();
		scanf("%d",&m);
		scanf("%I64d%I64d%I64d%I64d",&r[0].p[0].x,&r[0].p[0].y,&r[n].p[2].x,&r[n].p[2].y);
		r[0].p[1].x=r[0].p[0].x;
		r[0].p[1].y=r[n].p[2].y;
		r[n].p[3].x=r[n].p[2].x;
		r[n].p[3].y=r[0].p[0].y;
		for(int i=1;i<=n;i++)
		{
			scanf("%I64d%I64d",&p[i].x,&p[i].y);
		}
		sort(p+1,p+n+1,cmp);
		for(int i=1;i<=n;i++)
		{
			r[i-1].p[3].x=p[i].x;
			r[i-1].p[3].y=r[i-1].p[0].y;
			r[i-1].p[2].x=p[i].y;
			r[i-1].p[2].y=r[i-1].p[1].y;
			r[i].p[0].x=p[i].x;
			r[i].p[0].y=r[i-1].p[3].y;
			r[i].p[1].x=p[i].y;
			r[i].p[1].y=r[i-1].p[2].y;
		}
		for(int i=0;i<m;i++)
		{
			scanf("%d%d",&q.x,&q.y);
			search();
		}
		printf("Box\n");
		for(int i=0;i<=n;i++)
		{
			if(cnt[i])
				cnt1[cnt[i]]++;
		}
		for(int i=1;i<=m;i++)
		{
			if(cnt1[i])
				printf("%d: %I64d\n",i,cnt1[i]);
		}
	}
	return 0;
}