BZOJ1185[HNOI2007] 最小矩形覆盖

原题链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1185

最小矩形覆盖

题目描述

给定一些点的坐标，要求求能够覆盖所有点的最小面积的矩形，输出所求矩形的面积和四个顶点坐标

输入输出格式

输入格式：

第一行为一个整数n（3<=n<=50000），从第2至第n+1行每行有两个浮点数，表示一个顶点的x和y坐标，不用科学计数法

输出格式：

第一行为一个浮点数，表示所求矩形的面积（精确到小数点后5位），接下来4行每行表示一个顶点坐标，要求第一行为y坐标最小的顶点，其后按逆时针输出顶点坐标.如果用相同y坐标，先输出最小x坐标的顶点

输入输出样例

输入样例#1：

6 1.0 3.00000
1 4.00000
2.0000 1
3 0.0000
3.00000 6
6.0 3.0

输出样例#1：

18.00000
3.00000 0.00000
6.00000 3.00000
3.00000 6.00000
0.00000 3.00000

题解

结论：这个矩形的一条边一定与凸包的一条边重合。  
我们只需要以一条边为基准，用另外3条边去“卡”凸包即可，只需要通过叉积最大找对面的点，用点积最大找最右侧的点，点积最小找最左的点，以这四个点作为与矩形相切的点去更新矩形面积，四个顶点即可。另外，在更新的时候以逆时针顺序储存点，这样就可以方便的输出。

代码

#include<bits/stdc++.h>
#define db double
using namespace std;
const int M=5e4+5;
const db eps=1e-8;
struct pt{db x,y;};
int sig(db x){return (x>eps)-(x<-eps);}
pt operator -(pt a,pt b){return (pt){a.x-b.x,a.y-b.y};}
pt operator +(pt a,pt b){return (pt){a.x+b.x,a.y+b.y};}
db operator *(pt a,pt b){return a.x*b.y-a.y*b.x;}
db operator /(pt a,pt b){return a.x*b.x+a.y*b.y;}
pt operator *(pt a,db b){return (pt){a.x*b,a.y*b};}
pt operator /(pt a,db b){return (pt){a.x/b,a.y/b};}
bool operator <(pt a,pt b){return !sig(a.y-b.y)?a.x<b.x:a.y<b.y;}
db area(pt a,pt b,pt c){return (b-a)*(c-a);}
db shad(pt a,pt b,pt c){return (b-a)/(c-a);}
db sqr(db x){return x*x;}
db dis(pt a,pt b){return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));}
bool cmp(pt a,pt b){return sig(a.x-b.x)?a.y<b.y:a.x<b.x;}
pt p[M],sta[M],rec[5];
int top,n;
void in()
{
    scanf("%d",&n);
    for(int i=1;i<=n;++i)
    scanf("%lf%lf",&p[i].x,&p[i].y);
}
void tubao()
{
    top=-1;
    sort(p+1,p+1+n);
    for(int i=1;i<=n;++i)
    {
        while(top>0&&sig(area(sta[top-1],sta[top],p[i]))<=0)--top;
        sta[++top]=p[i];
    }
    int k=top;
    for(int i=n-1;i>=1;--i)
    {
        while(top>k&&sig(area(sta[top-1],sta[top],p[i]))<=0)--top;
        sta[++top]=p[i];
    }
}
int x[M];
db fuck(db x){return !sig(x)?fabs(x):x;}
void ac()
{
    int le=1,ri=1,up=1;
    db L,R,H,D,tmp,ans=1e60;
    for(int i=0;i<top;++i)
    {
        D=dis(sta[i],sta[i+1]);
        while(sig(area(sta[i],sta[i+1],sta[up])-area(sta[i],sta[i+1],sta[up+1]))<=0)up=(up+1)%top;
        while(sig(shad(sta[i],sta[i+1],sta[ri])-shad(sta[i],sta[i+1],sta[ri+1]))<=0)ri=(ri+1)%top;
        if(i==0)le=up;
        while(sig(shad(sta[i],sta[i+1],sta[le])-shad(sta[i],sta[i+1],sta[le+1]))>=0)le=(le+1)%top;
        L=shad(sta[i],sta[i+1],sta[le])/D;
        R=shad(sta[i],sta[i+1],sta[ri])/D;
        H=area(sta[i],sta[i+1],sta[up])/D;
        H=H>0?H:-H;
        tmp=(R-L)*H;
        if(tmp<ans)
        {
            ans=tmp;
            rec[0]=sta[i]+(sta[i+1]-sta[i])*(R/D);
            rec[1]=rec[0]+(sta[ri]-rec[0])*(H/dis(sta[ri],rec[0]));
            rec[2]=rec[1]-(rec[0]-sta[i])*((R-L)/dis(sta[i],rec[0]));
            rec[3]=rec[2]-(rec[1]-rec[0]);
        }
    }
    ans=fabs(ans);
    printf("%.5lf\n",ans);
    int fir=0;
    for(int i=1;i<=3;++i)
    if(rec[i]<rec[fir])fir=i;
    for(int i=0;i<=3;++i)
    printf("%.5lf %.5lf\n",fuck(rec[(i+fir)%4].x),fuck(rec[(i+fir)%4].y));
}
int main()
{
    in();tubao();
    ac();
    return 0;
}