题目链接:
HDU 2202 最大三角形
题意:
给出n个点,求最大三角形面积。
分析:
先求出凸包顶点,因为最大三角形的顶点一定都是凸包顶点。
接下来就暴力枚举所有凸包顶点(复杂度是凸包顶点数^3)。
有旋转卡壳的算法,然而并不会啊。。。。。。
弱菜的自述:
起初是把最大三角形的一条边当成必须是凸包的一条边,所以WA的很彻底啊!
实际上只是顶点必须是凸包顶点,而边并不一定!
//2096K 93MS
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <cmath>
using namespace std;
const int MAX_N=50010;
int n;
struct Point {
double x,y;
Point () {}
Point (double x,double y) : x(x),y(y) {
}
Point operator + (const Point& rhs) const {
return Point(x+rhs.x,y+rhs.y);
}
Point operator - (const Point& rhs) const {
return Point(x-rhs.x,y-rhs.y);
}
Point operator * (const double d) const {
return Point(d*x,d*y);
}
double cross(const Point& rhs) const {
return (x*rhs.y-y*rhs.x);
}
double dis(const Point& rhs) const {
return sqrt((x-rhs.x)*(x-rhs.x)+(y-rhs.y)*(y-rhs.y));
}
}point[MAX_N],vertex[MAX_N];
bool cmp_x(const Point a,const Point b)
{
if(a.x==b.x) return a.y<b.y;
else return a.x<b.x;
}
int Andrew()
{
sort(point,point+n,cmp_x);
int k=0;
for(int i=0;i<n;i++){
while(k>1&&(vertex[k-1]-vertex[k-2]).cross(point[i]-vertex[k-1])<=0) k--;
vertex[k++]=point[i];
}
int m=k;
//printf("m-1=%d\n",m-1);
for(int i=n-2;i>=0;i--){
while(k>m&&(vertex[k-1]-vertex[k-2]).cross(point[i]-vertex[k-1])<=0) k--;
vertex[k++]=point[i];
}
if(k>1) k--;
return k;
}
inline double GetArea(Point a,Point b,Point c)//计算顶点为a,b,c的三角形面积
{
double len1=a.dis(b);
double len2=a.dis(c);
double len3=b.dis(c);
double p=(len1+len2+len3)/2;
return sqrt(p*(p-len1)*(p-len2)*(p-len3));
}
void solve()
{
int total=Andrew();
if(total<3) {
printf("%.2f\n",0.0);
return;
}
//暴力枚举所有顶点形成的三角形
double ans=0.0;
for(int i=0;i<total;i++){
for(int j=0;j<total;j++){
for(int k=0;k<total;k++){
ans=max(ans,GetArea(vertex[i],vertex[j],vertex[k]));
}
}
}
printf("%.2f\n",ans);
}
int main()
{
freopen("Hin.txt","r",stdin);
while(~scanf("%d",&n)){
for(int i=0;i<n;i++){
scanf("%lf%lf",&point[i].x,&point[i].y);
}
solve();
}
return 0;
}