首页 > 代码库 > HDU 3934

HDU 3934

/*这是用的有旋转卡壳的思想。首先确定i,j,对k进行循环,知道找到第一个k使得cross(i,j,k)>cross(i,j,k+1),如果k==i进入下一次循环。对j,k进行旋转,每次循环之前更新最大值,然后固定一个j,同样找到一个k使得cross(i,j,k)>cross(i,j,k+1)。对j进行++操作,继续进行下一次,知道j==k为止。*/ #include <iostream>#include <cstdio>#include <algorithm>#include <cmath>using namespace std;struct point{	double x,y;}p[1000100];int n;int ans[1000100],st[1000100],cnt,stop;bool cmp(point A, point B){	if(A.y<B.y) return true;	else if(A.y==B.y){		if(A.x<B.x) return true;	}	return false;}double multi(point a,point b,point c){	point p1; p1.x=a.x-c.x; p1.y=a.y-c.y;	point p2; p2.x=b.x-c.x; p2.y=b.y-c.y;	return p1.x*p2.y-p1.y*p2.x;}void slove(){	cnt=stop=0;	st[stop++]=0; st[stop++]=1;	for(int i=2;i<n;i++){		while(stop>1&&multi(p[i],p[st[stop-1]],p[st[stop-2]])>=0) stop--;		st[stop++]=i;	}	for(int i=0;i<stop;i++)	ans[cnt++]=st[i];	stop=0; st[stop++]=n-1; st[stop++]=n-2;	for(int i=n-3;i>=0;i--){		while(stop>1&&multi(p[i],p[st[stop-1]],p[st[stop-2]])>=0) stop--;		st[stop++]=i;	}	for(int i=1;i<stop-1;i++)	ans[cnt++]=st[i];/*	for(int i=0;i<cnt;i++)	cout<<ans[i]<<endl;	cout<<endl;*/}double Triangle(point a,point b,point c){	point p1; p1.x=a.x-c.x; p1.y=a.y-c.y;	point p2; p2.x=b.x-c.x; p2.y=b.y-c.y;	return fabs((p1.x*p2.y-p1.y*p2.x)*1.0)/2.0;}double Area(){	int q; int j;	double anst=0;	for(int i=0;i<cnt;i++){		j=(i+1)%cnt;		q=(j+1)%cnt;		while(Triangle(p[ans[i]],p[ans[j]],p[ans[q]])<=Triangle(p[ans[i]],p[ans[j]],p[ans[(q+1)%cnt]])&&q!=i)		q=(q+1)%cnt;    //枚举了当前最远的K点 		anst=max(anst,Triangle(p[ans[i]],p[ans[j]],p[ans[q]]));		if(q==i) continue;		while(j!=i&&q!=i){			anst=max(anst,Triangle(p[ans[i]],p[ans[j]],p[ans[q]]));			while(Triangle(p[ans[i]],p[ans[j]],p[ans[q]])<=Triangle(p[ans[i]],p[ans[j]],p[ans[(q+1)%cnt]])&&q!=i)			q=(q+1)%cnt;			j=(j+1)%cnt;		}			}	return anst;}int main(){	while(scanf("%d",&n)!=EOF){	//	if(n==-1) break;		for(int i=0;i<n;i++){			scanf("%lf%lf",&p[i].x,&p[i].y);		}		sort(p,p+n,cmp);		slove();		double anst=0;		anst=max(anst,Area());		printf("%.2lf\n",anst);	}	return 0;}