斜率小于0的连线数量-归并排序

题目:

二维平面上N个点之间共有C(n,2)条连线。求这C(n,2)条线中斜率小于0的线的数量。
二维平面上的一个点,根据对应的X Y坐标可以表示为(X,Y)。例如:(2,3) (3,4) (1,5) (4,6),其中(1,5)同(2,3)(3,4)的连线斜率 < 0,因此斜率小于0的连线数量为2。
 
Input
第1行:1个数N,N为点的数量(0 <= N <= 50000)
第2 - N + 1行:N个点的坐标,坐标为整数。(0 <= X[i], Y[i] <= 10^9)
Output
输出斜率小于0的连线的数量。(2,3) (2,4)以及(2,3) (3,3)这2种情况不统计在内。
Input 示例
4
2 3
3 4
1 5
4 6
Output 示例
2

题目分析:

这道题目还有另外一个表述就是求斜率大于0的数量,其实本质是一样的。
一看题目的数据量,就知道肯定不能直接枚举。
首先稍微分析一下我们发现,斜率小于0的情况的直线就是直角坐标系中经过二四象限的直线,由于所有点都在第一象限,所以说,如果两个点满足斜率小于0,比如(x1,y1)(x2,y2)一定满足(x2>x1,y1>y2)或者相反(x1>x2,y2>y1)。其实也就是经过二四象限直线上的点满足的条件。

那么既然如此,我们可以先对x从小到大排序,那么对于每一个点,能够和这个点组成斜率小于0的点就是后面所有点中纵坐标y小于当前点纵坐标的点,那么直线的数量也就是比当前y值小的点的个数。判断到这里,一般人都会想到枚举,直接对于每一个点,求后面y值小于当前y值的点的个数,然后求和即可。但是这样可定会超时。
            对于逆序数熟悉的人,就会发现,y值大于后面的y值,实际上就是逆序数。什么是逆序数,就是每个数有一个下标i,j;当i<j时,如果a[i]>a[j]那么这就是一个逆序数对。
          因此,求所有y点,后面小于他的点的个数的总和,其实就是求由y坐标组成的序列的逆序对数。


逆序对数很明显,归并排序直接就解决了。
这里注意一个细节,那就是斜率不存在以及斜率为0的情况如何规避。斜率为0的也就是y值一样的,这里求逆序数的时候不判断相等即可规避。斜率不存在也就是x相同的,如何规避呢,这里就是在按x排序的时候,如果两个点的x一样,那么他们的y就从小到大排序,这样对于x一样的点就不会有逆序对,这样当然后面求逆序对的时候也就不会计算在内了。

代码:

#include<iostream>
#include <algorithm>
#include <vector>
using namespace std;

const int MAXN = 50002;
typedef struct POINT
{
	int x,y;
}Point;
Point point[MAXN];
int N;
__int64 sum;
bool cmp(const Point &a, const Point &b)
{
	if (a.x == b.x)
	{
		return a.y < b.y;
	}

	return a.x < b.x;
}

void Merge(vector<int>&Y,int l, int mid, int r)
{
	vector<int>L;
	vector<int>R;
	for (int i = l; i <= mid; ++ i)
	{
		L.push_back(Y[i]);
	}
	for (int i = mid + 1; i <= r; ++ i)
	{
		R.push_back(Y[i]);
	}
	int p = 0,q = 0,k = l;
	while (p < L.size() && q < R.size())
	{
		if (L[p] <= R[q])
		{
			Y[k++] = L[p++];
		}
		else
		{
			sum += L.size() - p;//mid - l - p + 1;
			Y[k++] = R[q++];
		}
	}
	while (p<L.size())
	{
		Y[k++] = L[p++];
	}
	while (q<R.size())
	{
		Y[k++] = R[q++];
	}
}

void MergeSort(vector<int> &Y,int l, int r)
{
	if (l < r)
	{
		int mid = (l+r)>>1;
		MergeSort(Y,l,mid);
		MergeSort(Y,mid+1,r);
		Merge(Y,l,mid,r);
	}
}
int main()
{
	while (cin >> N)
	{
		sum = 0;
		for (int i = 0; i < N; ++ i)
		{
			cin >> point[i].x >> point[i].y;
		}
		sort(point,point+N,cmp);

		vector<int>Y;
		for (int i = 0; i < N; ++ i)
		{
			Y.push_back(point[i].y);
		}
		MergeSort(Y,0,N-1);
		cout << sum << endl;
	}
	return 0;
}



 

posted @ 2013-11-01 19:22  pangbangb  阅读(642)  评论(0编辑  收藏  举报