原题链接: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;
}