HDU 2202 最大三角形（求凸包中最大三角形面积）

题目链接：
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;
}