BZOJ1185: [HNOI2007]最小矩形覆盖(洛谷P3187)

版权声明:蒟蒻Blog随意转载 https://blog.csdn.net/a1799342217/article/details/79947066

旋转卡壳

BZOJ题目传送门
洛谷题目传送门

旋转卡壳求最小面积外接矩形。

可以发现矩形的一条边一定和凸包的一条边重合。

那么我们枚举每一条边,计算出离这条边最远、最左、最右的三个点。这样就构造出了这条边的最小面积外接矩形,更新答案即可。

至于怎么计算。。。看代码吧!

代码:

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 50005
#define F friend
#define eps 1e-8
#define sqr(x) ((x)*(x))
#define abs(x) ((x)>0?(x):-(x))
using namespace std;
typedef double DB;
struct P{
    DB x,y;
    F bool operator < (P a,P b){
        return abs(a.y-b.y)<eps?a.x<b.x:a.y<b.y;
    }
    F P operator * (P a,DB b){ return (P){a.x*b,a.y*b}; }
    F P operator + (P a,P b){ return (P){a.x+b.x,a.y+b.y}; }
    F P operator - (P a,P b){ return (P){a.x-b.x,a.y-b.y}; }
    F DB dis(P a,P b){ return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y)); }
    F DB cs(P a,P b,P c){
        return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
    }
    F DB dt(P a,P b,P c){
        return (b.x-a.x)*(c.x-a.x)+(b.y-a.y)*(c.y-a.y);
    }
}t[N],s[N],v[4];
int n,tp,f; DB S=1e50;
inline bool cmp(P a,P b){
    double f=cs(t[1],a,b);
    if (abs(f)>eps) return f>0;
    return dis(a,t[1])-dis(b,t[1])<eps;
}
inline void GR(){
    sort(t+1,t+n+1),s[tp=1]=t[1],sort(t+2,t+n+1,cmp);
    for (int i=2;i<=n;s[++tp]=t[i++])
        while (tp>1&&cs(s[tp-1],s[tp],t[i])<eps) tp--;
}
inline void RC(){
    int l=1,r=1,h=1; s[0]=s[tp];
    for (int i=0;i<tp;i++){
        DB L,R,H,d=dis(s[i],s[i+1]);
        while (cs(s[i],s[i+1],s[h])-cs(s[i],s[i+1],s[h+1])<eps)
            h=(h+1)%tp; H=abs(cs(s[i],s[i+1],s[h])/d);
        while (dt(s[i],s[i+1],s[r])-dt(s[i],s[i+1],s[r+1])<eps)
            r=(r+1)%tp; R=dt(s[i],s[i+1],s[r])/d;
        if (!i) l=r;
        while (dt(s[i],s[i+1],s[l+1])-dt(s[i],s[i+1],s[l])<eps)
            l=(l+1)%tp; L=dt(s[i],s[i+1],s[l])/d;
        if ((R-L)*H<S){
            S=(R-L)*H;
            v[0]=s[i]+(s[i+1]-s[i])*(R/d);
            v[1]=v[0]+(s[r]-v[0])*(H/dis(v[0],s[r]));
            v[2]=v[1]-(v[0]-s[i])*((R-L)/dis(v[0],s[i]));
            v[3]=v[2]-(v[1]-v[0]);
        }
    }
}
int main(){
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
        scanf("%lf%lf",&t[i].x,&t[i].y);
    GR(),RC(),printf("%.5f\n",S+eps);
    for (int i=1;i<=3;i++)
        if (v[i]<v[f]) f=i;
    for (int j=0;j<4;j++,f=(f+1)%4)
        printf("%.5f %.5f\n",v[f].x+eps,v[f].y+eps);
    return 0;
}
阅读更多
想对作者说点什么?
相关热词

博主推荐

换一批

没有更多推荐了,返回首页