1185: [HNOI2007]最小矩形覆盖

这计算几何……果然很烦……

发现自己不会旋转卡壳,补了下,然后发现求凸包也不会……

凸包:找一个最左下的点,其他点按照与它连边的夹角排序,然后维护一个栈用斜率判定。

旋转卡壳:枚举一条边,用叉积和点积维护另外三条边(联系叉积和点积的几何意义,叉积最大即为对边,点积最大最小即为邻边)

找来5份标程对拍……啥?4个不同的输出,相同的两个完全是错的……

自己拍吧T_T(花了一下午)

神TM卡double

#include<cmath>
#include<cstdio>
#include<algorithm>
#define MN 510001
#define ld long double
#define eps 1e-9
using namespace std; struct po{ld x,y;}p[MN],a,b,MMH[];
ld mmh=/.;
int n,m,st[MN],top,j,k,l;
po operator - (po a,po b){return po{a.x-b.x,a.y-b.y};}
po operator + (po a,po b){return po{a.x+b.x,a.y+b.y};}
po operator * (po a,ld x){return po{a.x*x,a.y*x};}
inline ld ABS(ld x){return x<?-x:x;}
inline ld sqr(ld x){return x*x;}
inline ld dis(po a,po b){return sqrt(sqr(a.x-b.x)+sqr(a.y-b.y));}
inline ld dj(po a,po b){return a.y*b.x-b.y*a.x;}
inline ld dj(po a,po b,po c){return (b.y-a.y)*(c.x-a.x)-(c.y-a.y)*(b.x-a.x);}
inline ld dj(po a,po b,po c,po d){return (b.y-a.y)*(d.x-c.x)-(d.y-c.y)*(b.x-a.x);}
inline ld det(po a,po b){return a.x*b.x+a.y*b.y;}
inline ld det(po a,po b,po c){return (b.x-a.x)*(c.x-a.x)+(c.y-a.y)*(b.y-a.y);}
bool cmp(po a,po b){
ld t=dj(p[],a,b);
if (abs(t)>eps) return t<;else return dis(a,p[])<dis(b,p[]);
}
int main(){
register int i;
scanf("%d",&n);
for (i=;i<=n;i++) scanf("%Lf%Lf",&p[i].x,&p[i].y);
for (i=;i<=n;i++) if (p[i].y<p[].y||(p[i].y==p[].y&&p[i].x<p[].x)) swap(p[],p[i]);
sort(p+,p++n,cmp); for (top=,i=;i<=n;i++){
while (top>&&dj(p[st[top-]],p[st[top]],p[i])>=-eps) top--;
st[++top]=i;
}
st[]=st[top];st[top+]=st[];
for (i=,j=k=l=;i<=top;i++){
while (dj(p[st[i]],p[st[i+]],p[st[j]])+eps>dj(p[st[i]],p[st[i+]],p[st[j+]])) j=j==top?:j+; while (det(p[st[i]],p[st[i+]],p[st[k]])-eps<det(p[st[i]],p[st[i+]],p[st[k+]])) k=k==top?:k+; if (i==) l=k;
while (det(p[st[i]],p[st[i+]],p[st[l]])+eps>det(p[st[i]],p[st[i+]],p[st[l+]])) l=l==top?:l+; a=p[st[i+]]-p[st[i]];
ld D=dis(p[st[i]],p[st[i+]]);
ld H=ABS(det(p[st[k]]-p[st[i+]],a))/D+ABS(det(p[st[l]]-p[st[i]],a))/D+D;
a=p[st[i+]]-p[st[i]];a=po{-a.y,a.x};
ld W=(ABS(det(p[st[k]]-p[st[i]],a))+ABS(det(p[st[k]]-p[st[j]],a)))/D;
ld S=H*W;
if (S<mmh){
mmh=S;
a=p[st[i+]]-p[st[i]];
MMH[]=p[st[i+]]+(p[st[i+]]-p[st[i]])*(ABS(det(p[st[k]]-p[st[i+]],a))/D/D);
a=p[st[i+]]-p[st[i]];a=po{-a.y,a.x};
MMH[]=MMH[]+a*(W/D);
a=p[st[i+]]-p[st[i]];a=po{-a.x,-a.y};
MMH[]=MMH[]+a*(H/D);
a=p[st[i+]]-p[st[i]];a=po{a.y,-a.x};
MMH[]=MMH[]+a*(W/D);
}
}
printf("%.5Lf\n",mmh);
i=;
for (int j=;j<;j++)
if (MMH[j].y<MMH[i].y||(MMH[j].y==MMH[i].y&&MMH[j].x<MMH[i].x)) i=j;
for (int j=;j<;j++) printf("%.5Lf %.5Lf\n",ABS(MMH[(i+j)%].x),ABS(MMH[(i+j)%].y));
}

BZOJ:1185: [HNOI2007]最小矩形覆盖的更多相关文章

  1. 洛谷 P3187 BZOJ 1185 [HNOI2007]最小矩形覆盖 (旋转卡壳)

    题目链接: 洛谷 P3187 [HNOI2007]最小矩形覆盖 BZOJ 1185: [HNOI2007]最小矩形覆盖 Description 给定一些点的坐标,要求求能够覆盖所有点的最小面积的矩形, ...

  2. BZOJ 1185: [HNOI2007]最小矩形覆盖 [旋转卡壳]

    1185: [HNOI2007]最小矩形覆盖 Time Limit: 10 Sec  Memory Limit: 162 MBSec  Special JudgeSubmit: 1435  Solve ...

  3. bzoj 1185 [HNOI2007]最小矩形覆盖 凸包+旋转卡壳

    题目大意 用最小矩形覆盖平面上所有的点 分析 有一结论:最小矩形中有一条边在凸包的边上,不然可以旋转一个角度让面积变小 简略证明 我们逆时针枚举一条边 用旋转卡壳维护此时最左,最右,最上的点 注意 注 ...

  4. ●BZOJ 1185 [HNOI2007]最小矩形覆盖

    题链: http://www.lydsy.com/JudgeOnline/problem.php?id=1185 题解: 计算几何,凸包,旋转卡壳 结论:矩形的某一条边在凸包的一条边所在的直线上. ( ...

  5. BZOJ 1185 [HNOI2007]最小矩形覆盖:凸包 + 旋转卡壳

    题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1185 题意: 给出二维平面上的n个点,问你将所有点覆盖的最小矩形面积. 题解: 先找出凸 ...

  6. bzoj 1185 [HNOI2007]最小矩形覆盖——旋转卡壳

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1185 矩形一定贴着凸包的一条边.不过只是感觉这样. 枚举一条边,对面的点就是正常的旋转卡壳. ...

  7. BZOJ 1185: [HNOI2007]最小矩形覆盖-旋转卡壳法求点集最小外接矩形(面积)并输出四个顶点坐标-备忘板子

    来源:旋转卡壳法求点集最小外接矩形(面积)并输出四个顶点坐标 BZOJ又崩了,直接贴一下人家的代码. 代码: #include"stdio.h" #include"str ...

  8. BZOJ 1185 [HNOI2007]最小矩形覆盖 ——计算几何

    程序写的太垃圾,卡不过去. GG,甘拜下风. #include <map> #include <cmath> #include <queue> #include & ...

  9. 1185: [HNOI2007]最小矩形覆盖

    1185: [HNOI2007]最小矩形覆盖 Time Limit: 10 Sec  Memory Limit: 162 MBSec  Special JudgeSubmit: 1426  Solve ...

随机推荐

  1. iOS UICollectionView(转三)

    上篇博客的实例是自带的UICollectionViewDelegateFlowLayout布局基础上来做的Demo, 详情请看<iOS开发之窥探UICollectionViewControlle ...

  2. C:函数:功能:实现字符数组中所有字母的倒序存放并输出

    前两天小测碰到一道题,建立一个函数,功能:实现字符数组中所有字母的倒序存放并输出,一开始觉得简单跟数字数组差不多,运行一下发现很多格式错误,这些是不必要的错误,现在就来说下,先说一下代码思路:定义一个 ...

  3. 贪心算法——Fence Repair(POJ 3253)

    题目描述 农夫约翰为了修理栅栏,要将一块很长的木板切割成N块.准备切成的木板长度为L1,L2,L3--LN,未切割前木板的长度恰好为切割后木板长度的总和.每次切断木板时,需要的开销为这块木板的长度.请 ...

  4. n年前,我没钱但年轻,我怕n年后我老时,还是一无所成——2017我的收获和反思

    记得当年我刚从学校里出来时,应该和现在的95后差不多,当时还是很惶恐的,怕找不到工作,怕无法挣到足够的钱买房子支撑家庭,(当然还有其它的担心点),却唯独没意识到自己拥有着最宝贵的财富:年轻. 年轻意味 ...

  5. 完善chrome翻译插件ChaZD,支持有道智云api

    首先放上该项目的github地址:https://github.com/codethereforam/ChaZD 之前想找一个chrome支持划词翻译的插件,最终在知乎上看到了这个回答,推荐的是Cha ...

  6. ArcGIS API for JavaScript 4.2学习笔记[24] 【IdentifyTask类】的使用(结合IdentifyParameters类)(第七章完结)

    好吧,我都要吐了. 接连三个例子都是类似的套路,使用某个查询参数类的实例,结合对应的Task类,对返回值进行取值.显示. 这个例子是Identify识别,使用了TileLayer这种图层,数据来自Se ...

  7. Fragment生命周期及实现点击导航图片切换fragment,Demo

    PS:Fragment简介 Fragment是Android3.0后引入的一个新的API,他出现的初衷是为了适应大屏幕的平板电脑, 当然现在他仍然是平板APP UI设计的宠儿,而且我们普通手机开发也会 ...

  8. application19事件 20多少步骤 具体20多少只有微软知道!!!

  9. vue2.0的瀑布流组件-使用说明

    做一个小项目,需要瀑布流,就选他了,先看看效果 使用瀑布流布局组件:vue-waterfall-easy 下载引入: 方式一:直接从git上复制组件的完整代码,引入vue组件文件即可 import v ...

  10. 活动倒计时-兼容ios

    最近要做一个活动,需要用倒计时,写好之后再IOS上无效,经过百度知道了,原来IOS不能识别格式"2017-11-09 --",所以要进行转换才有效 直接上代码了: <!DOC ...