版权声明:请尊重每一滴汗水,谢谢! https://blog.csdn.net/qq_38853493/article/details/80396161
凸包面积
1000(ms)
65535(kb)
1184 / 3813
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
无需找出凸包点亦可求出凸包面积!分三步走...
1.输入所有点
2.依X或Y排序,sort函数ok!
3.求面积,getArea()函数是精华!
@浅夏沫若:code
#include
#include
using namespace std;
struct point
{
int x;
int y;
};
point p[110];
double areas = 0;
void p_sort(int n)
{
point temp;
int index = 0;
for (int i = 0; i < n-1; i++)
{
index = i;
for (int j = i + 1; j < n; j++)
{
if (p[index].x > p[j].x)
index = j;
}
temp = p[i];
p[i] = p[index];
p[index] = temp;
}
} //点排序
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;
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 counter = 0;
cin >> counter;
while (counter--)
{
int n = 0;
cin >> n;
for (int i = 0; i < n; i++)
cin >> p[i].x >> p[i].y;
p_sort(n);
areas = 0;
getArea(0, n - 1, n);
getArea(n - 1, 0, n); //此处使用两次getArea(),请深入理解一下,数学问题!
printf("%.1lf\n", areas);
}
return 0;
}