凸包面积

Tags: 分治法

麦兜是个淘气的孩子。一天，他在玩钢笔的时候把墨水洒在了白色的墙上。再过一会，麦兜妈就要回来了，麦兜为了不让妈妈知道这件事情，就想用一个白色的凸多边形把墙上的墨点盖住。你能告诉麦兜最小需要面积多大的凸多边形才能把这些墨点盖住吗？ 现在，给出了这些墨点的坐标，请帮助麦兜计算出覆盖这些墨点的最小凸多边形的面积。

输入

多组测试数据。第一行是一个整数T，表明一共有T组测试数据。
每组测试数据的第一行是一个正整数N(0< N < = 105)，表明了墨点的数量。接下来的N行每行包含了两个整数Xi和Yi（0<=Xi,Yi<=2000），表示每个墨点的坐标。每行的坐标间可能包含多个空格。

输出

每行输出一组测试数据的结果，只需输出最小凸多边形的面积。面积是个实数，小数点后面保留一位即可，不需要多余的空格。

样例输入

2
4
0 0
1 0
0 1
1 1
2
0 0
0 1

样例输出

1.0
0.0

实现代码：

#include"iostream"
#include"algorithm"
#include"cstdio"
using namespace std;
typedef struct point
{
	int x;
	int y;
};
point p[100];
double areas = 0;
bool cmp(point a,point b)         //将用于排序规则：按照x的值大小排升序，若x值相等，按y值排升序 
{
	if(a.x != b.x)
		return a.x > b.x;
	else
		return a.y < b.y;
}
void getArea(int l,int r,int n)    //这个部分灰常重要 
{
	if(l == r)
		return;
	int index = 0;
	double max = 0;
	double s = 0;
	for(int i = 0; i < n; i++)
	{
		if(i == l || i == r)
			continue;
			//某个数学计算公式，看看某度就ok； 
		s = (p[l].x*p[r].y + p[r].x*p[i].y + p[i].x*p[l].y - p[l].x*p[i].y - p[r].x*p[l].y - p[i].x*p[r].y)*0.5;
		if(s > max)
		{
			max = s;
			index = i;
		}
	}
	if(max > 0)
	{
		areas += max;
		getArea(l, index, n);
		getArea(index, r, n);
	}
}

int main()
{
	int num;
	cin >> num;
	while(num--)
	{
		int n = 0;
		cin >> n;
		for(int i = 0; i < n; i++)
			cin >> p[i].x >> p[i].y;
		sort(p, p + n, cmp);     //点排序 
		areas = 0;
		getArea(0, n - 1, n);
		getArea(n - 1, 0, n);
		printf("%.1lf\n",areas);
	}
	return 0;
}