POJ 2079 Triangle

题目连接：http://poj.org/problem?id=2079
这个题，给了你很多个点，然后要我们求其中三个点，让这三个点的面积最大。
由于点的数量很多，如果直接O(n^3)暴力的话是肯定会超时的。
稍微好一点的方法是去求一个凸包，如果三角形的面积要最大，那么这三个点必然位于凸包之上（凸包之后，点的数量就大大减少了）
但是此时在凸包上直接暴力也是不合理的。这里有一个O(n)的方法就是旋转卡壳。
就是不断的枚举三个点p,q,r（逆时针方向，顺时针方向可以自己任意定）

最后没举出一个最大值~
我的代码：

#include<stdio.h>
#include<math.h>
#include<algorithm>
#define eps 1e-8
#define maxn 50005

using namespace std;

struct point
{
double x;
double y;
};
point p[maxn];
point s[maxn];

double cross(point p1,point p2,point p0)
{
return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

double dis(point a,point b)
{
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

bool cmp(point p1,point p2)
{
if(cross(p1,p2,p[0])>0||(cross(p1,p2,p[0])==0&&(dis(p1,p[0])<dis(p2,p[0]))))
return true;
else
return false;
}

double AREA(point a,point b,point c)
{
double t=fabs(cross(a,b,c));
return t/2;
}

double max(double a,double b)
{
if(a>b)
return a;
else
return b;
}

void RC(int n)
{
double area;
int pp,qq,rr;
int p=0,q=1,r=2;
area=AREA(s[p],s[q],s[r]);
while(1)
{
pp=p,qq=q,rr=r;
while(fabs(cross(s[(r+1)%n],s[q],s[p]))>fabs(cross(s[r],s[q],s[p])))
{
area=max(area,AREA(s[(r+1)%n],s[q],s[p]));
r=(r+1)%n;
}
while(fabs(cross(s[r],s[(q+1)%n],s[p]))>fabs(cross(s[r],s[q],s[p])))
{
area=max(area,AREA(s[r],s[(q+1)%n],s[p]));
q=(q+1)%n;
}
while(fabs(cross(s[r],s[q],s[(p+1)%n]))>fabs(cross(s[r],s[q],s[p])))
{
area=max(area,AREA(s[r],s[q],s[(p+1)%n]));
p=(p+1)%n;
}
if(pp==p&&qq==q&&rr==r)
r=(r+1)%n;
if(r==0)
break;
}
printf("%.2lf\n",area);
}

void graham(int n)
{
int top,i,tmp=0;
for(i=1;i<n;i++)
{
if((p[i].y<p[tmp].y)||((p[i].y==p[tmp].y)&&(p[i].x<p[tmp].x)))
tmp=i;
}
swap(p[0],p[tmp]);
sort(p+1,p+n,cmp);
top=1;
s[0]=p[0];
s[1]=p[1];
for(i=2;i<n;i++)
{
while(top&&cross(p[i],s[top],s[top-1])>=0)
top--;
s[++top]=p[i];
}
if(top+1<2)
{
printf("0.00\n");
return;
}
RC(top+1);
}

int main()
{
int n,i;
while(scanf("%d",&n)!=EOF)
{
if(n==-1)
break;
for(i=0;i<n;i++)
scanf("%lf%lf",&p[i].x,&p[i].y);
graham(n);
}
return 0;
}