Description
老师在计算几何这门课上给Eddy布置了一道题目,题目是这样的:给定二维的平面上n个不同的点,要求在这些点里寻找三个点,使他们构成的三角形拥有的面积最大。
Eddy对这道题目百思不得其解,想不通用什么方法来解决,因此他找到了聪明的你,请你帮他解决这个题目。
Input
输入数据包含多组测试用例,每个测试用例的第一行包含一个整数n,表示一共有n个互不相同的点,接下来的n行每行包含2个整数xi,yi,表示平面上第i个点的x与y坐标。你可以认为:3 <= n <= 50000 而且 -10000 <= xi, yi <= 10000.
Output
对于每一组测试数据,请输出构成的最大的三角形的面积,结果保留两位小数。
每组输出占一行。
Sample Input
3
3 4
2 6
3 7
6
2 6
3 9
2 0
8 0
6 6
7 7
Sample Output
1.50
27.00
分析
https://www.cnblogs.com/adelalove/p/8467472.html
学凸包与旋转卡壳我是看的这篇博客
老师给我们平面上n个点,让我们求出组合的三角形的最大面积
我们可以先将这n个点组成一个凸包,然后通过旋转卡壳求出凸包的最大直径
找出直径对应的两个点,然后通过遍历求出第三个点
这里面用上了叉积的概念,即一条边X另一条边
公式:
三角的面积等于同一个起点出发的两边的叉积*1/2 假设A B是两条边
面积就等于(A+B)/2
AC
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#define mst(a) memset(a,0,sizeof(a))
using namespace std;
const int maxn=50000+8;
struct node{
int x;
int y;
}p[maxn],ch[maxn];
int n;
bool cmp(node a,node b)
{
if(a.x==b.x)
{
return a.y<b.y;
}
return a.x<b.x;
}
double Cross(node s,node a,node b)
{
int x1=a.x-s.x;
int y1=a.y-s.y;
int x2=b.x-s.x;
int y2=b.y-s.y;
return x1*y2-x2*y1;
}
int Andrew()
{
sort(p,p+n,cmp);
int m=0;
for(int i=0;i<n;i++)
{
//下
while(m>1&&Cross(ch[m-2],ch[m-1],p[i])<0) m--;
ch[m++]=p[i];
}
int k=m;
for(int i=n-2;i>=0;i--)
{
//上
while(m>k&&Cross(ch[m-2],ch[m-1],p[i])<0) m--;
ch[m++]=p[i];
}
if(n>1) m--;
return m;
}
int main()
{
while(cin>>n)
{
for(int i=0;i<n;i++)
{
cin>>p[i].x>>p[i].y;
}
int m=Andrew();
double res=0;
for(int i=0;i<m;i++)
{
int q=1;
for(int j=i+1;j<m;j++)
{
while(Cross(ch[i],ch[j],ch[q+1])>Cross(ch[i],ch[j],ch[q]))
q=(q+1)%m;
res=max(res,Cross(ch[i],ch[j],ch[q]));
}
}
printf("%.2lf\n",res/2.0);
}
return 0;
}
学如逆水行舟,不进则退