版权声明:随意转载哦......但还是请注明出处吧: https://blog.csdn.net/dreaming__ldx/article/details/83142668
传送门
不难看出最后的矩形一定有一条边与凸包某条边重合。
因此先求出凸包,然后旋转卡壳求出当前最小矩形面积更新答案。
代码:
#include<bits/stdc++.h>
#define N 50005
#define eps 1e-9
using namespace std;
struct pot{
long double x,y;
inline pot operator+(const pot&a){return (pot){x+a.x,y+a.y};}
inline pot operator-(const pot&a){return (pot){x-a.x,y-a.y};}
inline long double operator^(const pot&a){return x*a.y-y*a.x;}
inline long double operator*(const pot&a){return x*a.x+y*a.y;}
inline pot operator*(const long double&a){return (pot){x*a,y*a};}
inline long double dist(){return sqrt(x*x+y*y);}
friend inline bool operator<(pot a,pot b){return fabs(a.y-b.y)<eps?a.x<b.x:a.y<b.y;}
inline pot rev(){return (pot){-y,x};}
}p[N],ans[5],now[5];
int n,q[N],top=0;
long double sum=1e18;
inline bool cmp(pot x,pot y){
long double tmp=(x-p[1])^(y-p[1]);
if(fabs(tmp)>eps)return tmp>0;
return (x-p[1]).dist()<(y-p[1]).dist();
}
inline void graham(){
int tmp=1;
for(int i=2;i<=n;++i)if(p[i]<p[tmp])tmp=i;
if(tmp^1)swap(p[tmp],p[1]);
sort(p+2,p+n+1,cmp),q[++top]=1;
for(int i=2;i<=n;++i){
while(top>=2&&((p[q[top]]-p[q[top-1]])^(p[i]-p[q[top-1]]))<eps)--top;
q[++top]=i;
}
q[0]=q[top];
}
int main(){
scanf("%d",&n);
for(int i=1;i<=n;++i)scanf("%Lf%Lf",&p[i].x,&p[i].y);
graham();
for(int i=0,l=1,r=1,h=1;i<top;++i){
long double Dis=(p[q[i+1]]-p[q[i]]).dist();
while(((p[q[i+1]]-p[q[i]])^(p[q[h+1]]-p[q[i]]))-((p[q[i+1]]-p[q[i]])^(p[q[h]]-p[q[i]]))>-eps)h=(h+1)%top;
while((p[q[i+1]]-p[q[i]])*(p[q[r+1]]-p[q[i]])-(p[q[i+1]]-p[q[i]])*(p[q[r]]-p[q[i]])>-eps)r=(r+1)%top;
if(!i)l=r;
while((p[q[i+1]]-p[q[i]])*(p[q[l+1]]-p[q[i]])-(p[q[i+1]]-p[q[i]])*(p[q[l]]-p[q[i]])<eps)l=(l+1)%top;
long double L=(p[q[i+1]]-p[q[i]])*(p[q[l]]-p[q[i]])/Dis,R=(p[q[i+1]]-p[q[i]])*(p[q[r]]-p[q[i]])/Dis;
long double H=fabs(((p[q[i+1]]-p[q[i]])^(p[q[h]]-p[q[i]]))/Dis),tmp=(R-L)*H;
if(tmp<sum){
sum=tmp;
ans[1]=p[q[i]]+(p[q[i+1]]-p[q[i]])*(R/Dis);
ans[2]=ans[1]+(p[q[r]]-ans[1])*(H/(ans[1]-p[q[r]]).dist());
ans[3]=ans[2]-(ans[1]-p[q[i]])*((R-L)/(p[q[i]]-ans[1]).dist());
ans[4]=ans[3]+ans[1]-ans[2];
}
}
printf("%.5Lf\n",sum);
int pos=1;
for(int i=1;i<=4;++i){
if(fabs(ans[i].x)<eps)ans[i].x=0;
if(fabs(ans[i].y)<eps)ans[i].y=0;
}
for(int i=2;i<=4;++i)if(ans[i]<ans[pos])pos=i;
for(int i=1;i<=4;++i){
printf("%.5Lf %.5Lf\n",ans[pos].x,ans[pos].y);
++pos;
if(pos==5)pos=1;
}
for(int i=1;i<=10;++i)puts("");
return 0;
}