杭电 2202 最大三角形

这是个凸包题，不怎会啊。网上找了个比较容易理解的

但是有点问题啊。就是有个测试结果不对尽然也AC了。挺郁闷的。

（毕竟还是比较容易理解的）

代码：

#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
struct point
{
    double x,y;
};
point p[50000],res[50000];
double area(point t0,point t1,point t2)//三点构成的三角形面积
{
    return abs(t0.x*t1.y+t1.x*t2.y+t2.x*t0.y-t0.y*t1.x-t1.y*t2.x-t2.y*t0.x)/2;
}
double dist(point p1,point p2)//计算两点间的面积
{
    return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
bool multi(point p0,point p1,point p2)//判断是否凸边
{
    return (p1.x-p0.x)*(p2.y-p0.y)>(p2.x-p0.x)*(p1.y-p0.y);
}
int cmp(point a,point b)//将点排序
{
    if(a.y!=b.y) 
        return a.y<b.y;
    if(a.y==b.y&&a.x!=b.x)
        return a.x<b.x;
}
int main()
{
    int i,j,k,len,n;
    while(cin>>n)
    {
        for(i=0;i<n;i++)
            cin>>p[i].x>>p[i].y;
        sort(p,p+n,cmp);
        res[0]=p[0];
        res[1]=p[1];
        int top=1;
        for(i=2;i<n;i++)
        {
            while(top&&multi(p[i],res[top],res[top-1])) 
                top--;
            res[++top]=p[i];
        }
        len=top;
        res[++top]=p[n-2];
        for(i=n-3;i>=0;i--)
        {
            while(top!=len&&multi(p[i],res[top],res[top-1]))
                top--;
            res[++top]=p[i];
        }//前面两个for是找到凸边
        double max=0,s;
        for(i=0;i<top-2;i++)//在凸边中找最大面积
            for(j=i+1;j<top-1;j++)
                for(k=j+1;k<top;k++)
                {
                    s=area(res[i],res[j],res[k]);
                    if(max<s)
                        max=s;
                }
        printf("%.2lf\n",max);
    }
    return 0;
}