poj2318 TOYS 算法简介和C语言实现

原题链接：http://poj.org/problem?id=2318

一、题目概述
给定一个矩形区域，该区域被一系列上边点到下边点的连线分割成若干个区域。这些分割线不会相交，所以说这些小区域都是水平方向排布的。现在要向矩形区域内投放若干个坐标确定的点，要你统计落在每个被分割出的小区域的点的个数。

二、测试数据
每组数据给出划分线的个数、点的个数、矩形的范围以及每条划分线和每个点的信息。划分线和点的最大规模不超过5000。但是一个测试有多组数据，题目没有说最多包含几组。

三、提出算法
题目的目的是确定每个点所属的最小区域。我们把划分线从左到右的排列记为p，这个题目已经排好序了，不用你废脑筋。把矩形的左边线添加到p的开头，把右边线添加到p的尾部，即我们把矩形的左右边线也作为分割线看待。我们发现：对于任意的a<=b<c<=d，如果点在p[a]、p[d]之间，那么就一定在p[b]、p[c]之间。这就很容易让我们联想到二分法。我们之所以用二分法寻找一个点，就是因为这个点属于一个范围，那么它就一定属于它的父范围。首先初始化点可能存在的范围为整个矩形，即p[0]到p[n+1]（假设题目给出n条划分线，那么算上我们之前添加的两条，一共有n+2条）。我们找到中间的划分线，看点在这个划分线的左边还是右边。如果在左，则范围缩小到该划分线左边的部分；反之亦然。然后继续二分，直到能把点限制在两条相邻的划分线之间。

那么我们怎么确定点是不是在两条划分线之间呢？可以从该点发生水平线，这条水平线一定和划分线有一个交点，这样我们就能得到两个交点。如果该点在两条划分线之间，那么该点的x坐标一定会在交点的x坐标之间。

四、复杂度分析
把划分线规模记为n，点规模记为m。对于每一个点，我们都要二分查找，每次查找的时间复杂度都为常数，所以确定一个点的时间复杂度为O(logn)。因此每组数据就需要O(m*logn)。

五、C语言实现

#include <stdio.h>
#include <stdlib.h>
#include <string.h> 

struct cardboard
{
	int U, L;
} a[5050];

int n, m, x1, y1, x2, y2;

int in(int i, int j, int x, int y)
{
	double left, right;
	
	if(a[i].U!=a[i].L)
	{
		left = (double)(y-y1)*(a[i].L-a[i].U)/(y2-y1) + a[i].U;
	}
	else
	{
		left = a[i].U;
	}
	
	if(a[j].U!=a[j].L)
	{
		right = (double)(y-y1)*(a[j].L-a[j].U)/(y2-y1) + a[j].U;
	}
	else
	{
		right = a[j].U;
	}
	
	if(x>=left && x<=right)
	{
		return 1;
	}
	else
	{
		return 0;
	}
}

int main()
{
	int i;
	int xp, yp;
	int low, high, mid;
	int bin[5050];
	
	while(scanf("%d%d%d%d%d%d", &n, &m, &x1, &y1, &x2, &y2)==6)
	{
		a[0].U = x1;
		a[0].L = x1;
		for(i=1;i<=n;i++)
		{
			scanf("%d%d", &a[i].U, &a[i].L);
		}
		a[n+1].U = x2;
		a[n+1].L = x2;
		
		memset(bin, 0, sizeof(bin));
		
		for(i=0;i<m;i++)
		{
			scanf("%d%d", &xp, &yp);
			low = 0;
			high = n+1;
			while(high-low>1)
			{
				mid = (low+high)/2;
				//printf("%d %d %d\n", low, mid, in(low, mid, xp, yp));
				if(in(low, mid, xp, yp))
				{
					high = mid;
				}
				else
				{
					low = mid;
				}
			}
			//printf("%d\n", low);
			bin[low]++;
		}
		
		for(i=0;i<=n;i++)
		{
			printf("%d: %d\n", i, bin[i]);
		}
		putchar('\n');
	}
	
	
	
	return 0;
}