BZOJ-1185 旋转卡壳

Description
给定一些点的坐标，要求求能够覆盖所有点的最小面积的矩形，
输出所求矩形的面积和四个顶点坐标
Input
第一行为一个整数n（3<=n<=50000）
从第2至第n+1行每行有两个浮点数，表示一个顶点的x和y坐标，不用科学计数法
Output
第一行为一个浮点数，表示所求矩形的面积（精确到小数点后5位），
接下来4行每行表示一个顶点坐标，要求第一行为y坐标最小的顶点，
其后按逆时针输出顶点坐标.如果用相同y坐标，先输出最小x坐标的顶点
Sample Input
6 1.0 3.00000
1 4.00000
2.0000 1
3 0.0000
3.00000 6
6.0 3.0
Sample Output
18.00000
3.00000 0.00000
6.00000 3.00000
3.00000 6.00000
0.00000 3.00000

先求凸包，然后旋转卡壳卡出4个点
如果答案更小，通过这4个点算矩形顶点更新答案

double minRectangleCover(Point A[]){//最小矩形面积覆盖/逆时针凸包//返回面积和四个点
		if(n<3)return 0.0;//特判n<3 的情况
		p[n]=p[0];
		double ans=-1;
		int r=1,k=1,q;
		for(int i=0;i<n;++i){//卡出离边p[i]-p[i+1] 最远的点
			while(sgn(cross(p[i],p[i+1],p[r+1])-cross(p[i],p[i+1],p[r]))>=0)r=(r+1)%n;
			//卡出p[i]-p[i+1] 方向上正向n 最远的点
			while(sgn(dot(p[i],p[i+1],p[k+1])-dot(p[i],p[i+1],p[k]))>=0)k=(k+1)%n;
			if(i==0)q=k;//卡出p[i]-p[i+1] 方向上负向最远的点
			while(sgn(dot(p[i],p[i+1],p[q+1])-dot(p[i],p[i+1],p[q]))<=0)q=(q+1)%n;
			double d=(p[i]-p[i+1]).len2();
			double tmp=cross(p[i],p[i+1],p[r])*(dot(p[i],p[i+1],p[k])-dot(p[i],p[i+1],p[q]))/d;
			if(ans<0||ans>tmp){
				ans=tmp;
				Point f=(p[i]-p[i+1]).rotright();
				Line a=Line(p[i],p[i+1]),b=Line(p[k],f.angle());
				A[0]=a.crosspoint(b);
				A[1]=A[0]+(f.trunc(a.dispointtoline(p[r])));
				f=(p[i]-p[i+1]);
				double L=Line(A[0],A[1]).dispointtoline(p[q]);
				A[2]=A[1]+f.trunc(L);
				A[3]=A[0]+f.trunc(L);
				A[0].zero(),A[1].zero(),A[2].zero(),A[3].zero();
			}
		}return ans;
	}

polygon a,b;
Point V[4];
int main(){
	int n;scanf("%d",&n);
	a.input(n);
	a.getconvex(b);
	printf("%.5lf\n",b.minRectangleCover(V));
	double y=inf;
	int pos;
	for(int i=0;i<4;++i)if(sgn(V[i].y-y)<0)y=V[i].y,pos=i;
	for(int i=0;i<4;++i)printf("%.5lf %.5lf\n",V[(i+pos)%4].x,V[(i+pos)%4].y);
	return 0;
}