poj 2079 求一个平面所有点构成最大三角形的面积

#include<cstdio> #include<cmath> #include<algorithm> #define eps 1e-8 #define max(x,y) (((x)>(y))?(x):(y)) using namespace std; int n,top; struct point { double x,y; }p[50000+16],stack[10000]; bool equal(double n) { return fabs(n-0)<eps; } double cp(const point p1,const point p2,const point p0) { return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x); } double dis(const point p1,const point p2) { return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y)); } bool cmp(const point p1,const point p2) { double res=cp(p1,p2,p[0]); if(equal(res)) return dis(p[0],p1)-dis(p[0],p2)<=0; else return res>0; } double rotate_calipers(int num) { int p=1,q=2,i; double ans=0; stack[++num]=stack[0]; for(i=0;i<num;i++) { while(cp(stack[i],stack[p],stack[q+1])>cp(stack[i],stack[p],stack[q])) q=(q+1)%num; ans=max(ans,cp(stack[i],stack[p],stack[q])/2.0); while(cp(stack[i],stack[p+1],stack[q])>cp(stack[i],stack[p],stack[q])) p=(p+1)%num; ans=max(ans,cp(stack[i],stack[p],stack[q])/2.0); } return ans; } int main() { int i,pos; double res=0; while(~scanf("%d",&n)) { if(n==-1) break; for(i=0;i<n;i++) scanf("%lf%lf",&p[i].x,&p[i].y); pos=0; for(i=1;i<n;i++) { if(p[pos].y>=p[i].y) { if(p[pos].y==p[i].y) { if(p[pos].x>p[i].x) pos=i; } else pos=i; } } point temp; temp=p[pos]; p[pos]=p[0]; p[0]=temp; sort(p+1,p+n,cmp); top=1; stack[0]=p[0]; stack[1]=p[1]; for(i=2;i<n;i++) { while(cp(stack[top],p[i],stack[top-1])<0) top--; stack[++top]=p[i]; } printf("%.2lf\n",rotate_calipers(top)); } return 0; }