[旋转卡壳] BZOJ 1185 [HNOI2007]最小矩形覆盖 && 2218 Uva10173 Smallest Bounding Rectangle

最新推荐文章于 2017-03-14 11:20:08 发布
最新推荐文章于 2017-03-14 11:20:08 发布
阅读量389 收藏
点赞数
分类专栏： 几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014609452/article/details/52373116
版权

感动死了 一年前卡精度 今天搞搞搞就过了

感天动地

一年前的代码真丑


1185

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<cstring>
#define cl(x) memset(x,0,sizeof(x))
#define eps 1e-6
using namespace std;
typedef long double ld;

inline int dcmp(ld a,ld b)
{
	if (fabs(a-b)<eps) return 0;
	if (a<b) return -1; return 1;
}

struct Point{
	ld x,y;
	Point(ld x=0,ld y=0):x(x),y(y) { }
	bool operator < (const Point &b) const{
		return dcmp(x,b.x)<0 || (dcmp(x,b.x)==0 && dcmp(y,b.y)<0);
	}
}P[50005],H[50005];
ld s[50005];

int n,cnt;

inline ld Dot(Point p1,Point p2,Point p0)
{
	return (p1.x-p0.x)*(p2.x-p0.x)+(p1.y-p0.y)*(p2.y-p0.y);
}

inline ld Cross(Point p1,Point p2,Point p0)
{
	return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

inline ld Dist(Point p1,Point p2)
{
	return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}

inline ld Dist(Point p1,Point p2,Point p0)
{
	return fabs(Cross(p1,p2,p0)/Dist(p1,p2));
}

inline ld Area(Point p1,Point p2,Point p0)
{
	return fabs(Cross(p1,p2,p0))/2;
}

inline void Qhull(int l,int r,Point a,Point b)
{
	int i=l-1,j=r+1,x=l;
	for (int k=l;k<=r;k++)
		if (dcmp(s[k],s[x])>0 || (dcmp(s[x],s[k])==0 && P[x]<P[k]))
			x=k;
	Point y=P[x];
	for (int k=l;k<=r;k++)
	{
		s[++i]=Cross(P[k],a,y);
		if (dcmp(s[i],0)>0) swap(P[k],P[i]); else s[i--]=0;
	}
	for (int k=r;k>=l;k--)
	{
		s[--j]=Cross(P[k],y,b);
		if (dcmp(s[j],0)>0) swap(P[k],P[j]); else s[j++]=0;
	}
	if (l<=i) Qhull(l,i,a,y);
	H[++cnt]=y;
	if (j<=r) Qhull(j,r,y,b);
}

Point AP[5];

inline bool cmp(Point a,Point b)
{
	return a.x*b.y-b.x*a.y>0;
}

inline void Sort()
{
	int k=1;
	for (int i=2;i<=4;i++)
		if (AP[i].y<AP[k].y || (AP[i].y==AP[k].y && AP[i].x<AP[k].y))
			k=i;
	swap(AP[1],AP[k]);
	sort(AP+2,AP+5,cmp);
}

struct Line{
	ld k,b;
	ld x;
};

inline Point Cha(Line a,Line b)
{
	if (a.k<-1e100 && a.b>1e100)
		return Point(a.x,b.k*a.x+b.b);
	if (b.k<-1e100 && b.b>1e100)
		return Point(b.x,a.k*b.x+a.b);	 
	ld x=(b.b-a.b)/(a.k-b.k);
	return Point(x,a.k*x+a.b);
}

inline void Get(Point A,Point B,Point C,Point D,Point E)
{
	Line Li1,Li2,Li3,Li4;
	if (dcmp(A.x,B.x))
		Li1.k=(B.y-A.y)/(B.x-A.x),Li1.b=A.y-A.x*Li1.k;
	else
		Li1.k=-1e130,Li1.b=1e130,Li1.x=A.x;
	if (dcmp(Li1.k,0))
	{
		Li2.k=-1/Li1.k;
		Li2.k=dcmp(Li2.k,0)?Li2.k:0.0;
		Li2.b=C.y-C.x*Li2.k;
	}
	else
		Li2.k=-1e130,Li2.b=1e130,Li2.x=C.x;
	if (dcmp(Li1.k,0))
		Li3.k=-1/Li1.k,Li3.k=dcmp(Li3.k,0)?Li3.k:0.0,Li3.b=D.y-D.x*Li3.k;
	else
		Li3.k=-1e130,Li3.b=1e130,Li3.x=D.x;
	if (fabs(Li1.k)<1e100)
		Li4.k=Li1.k,Li4.b=E.y-E.x*Li4.k;
	else
		Li4=Li1,Li4.x=E.x;
	AP[1]=Cha(Li1,Li2);
	AP[2]=Cha(Li1,Li3);
	AP[3]=Cha(Li4,Li2);
	AP[4]=Cha(Li4,Li3);
}

int main()
{
	ld ans,tmp,A,B; double xx,yy;
	Point X,Y,C,D;
	int L1,L2,L3;
	freopen("t.in","r",stdin);
	freopen("t.out","w",stdout);
		int x=1;
		scanf("%d",&n);
		for (int i=1;i<=n;i++)
		{
			scanf("%lf%lf",&xx,&yy);
			P[i].x=xx; P[i].y=yy;
			if (P[i]<P[x]) x=i;
		}
		if (n==3)
		{
			ld ret=1e20;
			if (dcmp(ret,Dist(P[1],P[2])*Dist(P[1],P[2],P[3]))<0)
				ret=Dist(P[1],P[2])*Dist(P[1],P[2],P[3]),Get(P[1],P[2],P[1],P[2],P[3]);
			if (dcmp(ret,Dist(P[2],P[3])*Dist(P[2],P[3],P[1]))<0)
				ret=Dist(P[2],P[3])*Dist(P[2],P[3],P[1]),Get(P[2],P[3],P[2],P[3],P[1]);
			if (dcmp(ret,Dist(P[1],P[3])*Dist(P[1],P[3],P[2]))<0)
				ret=Dist(P[1],P[3])*Dist(P[1],P[3],P[2]),Get(P[1],P[3],P[1],P[3],P[2]);
			Sort();
			printf("%.5lf\n",(double)ret);
			for (int i=1;i<=4;i++)
				printf("%.5lf %.5lf\n",(double)AP[i].x,(double)AP[i].y);
			return 0;
		}
		swap(P[1],P[x]);
		H[++cnt]=P[1];
		Qhull(2,n,P[1],P[1]);
		n=cnt;
		for (int i=1;i<=cnt;i++)
			P[i]=H[i];
		Y=H[n]; X=H[1];
		L1=2;
		for (int i=3;i<n;i++)
			if (dcmp(Dist(X,Y,P[i]),Dist(X,Y,P[L1]))>0)
				L1=i;
		L2=2;
		for (int i=3;i<n;i++)
			if (dcmp(Dot(Y,P[i],X)/Dist(X,Y),Dot(Y,P[L2],X)/Dist(X,Y))<0)
				L2=i;
		L3=2;
		for (int i=3;i<n;i++)
			if (dcmp(Dot(X,P[i],Y)/Dist(X,Y),Dot(X,P[L3],Y)/Dist(X,Y))<0)
				L3=i;
		B=Dist(X,Y,P[L1]);
		A=Dist(X,Y);
		C=X;
		tmp=Dot(Y,P[L2],X)/Dist(X,Y);
		if (dcmp(tmp,0.0)<0) A-=tmp,C=P[L2];
		D=Y;
		tmp=Dot(X,P[L3],Y)/Dist(X,Y);
		if (dcmp(tmp,0.0)<0) A-=tmp,D=P[L3];
		ans=A*B; Get(X,Y,C,D,P[L1]);
		for (int i=1;i<n;i++)
		{
			Y=P[i]; X=P[i+1];
			while (dcmp(Dist(X,Y,P[L1]),Dist(X,Y,P[L1%n+1]))<0) 
				L1=L1%n+1;
			while (dcmp(Dot(Y,P[L2%n+1],X)/Dist(X,Y),Dot(Y,P[L2],X)/Dist(X,Y))<0)
				L2=L2%n+1;
			while (dcmp(Dot(X,P[L3%n+1],Y)/Dist(X,Y),Dot(X,P[L3],Y)/Dist(X,Y))<0)
				L3=L3%n+1;
			B=Dist(X,Y,P[L1]);
			A=Dist(X,Y);
			C=X;
			tmp=Dot(Y,P[L2],X)/Dist(X,Y);
			if (dcmp(tmp,0.0)<0) A-=tmp,C=P[L2];
			D=Y;
			tmp=Dot(X,P[L3],Y)/Dist(X,Y);
			if (dcmp(tmp,0.0)<0) A-=tmp,D=P[L3];
			if (dcmp(ans,A*B)>0)
				ans=A*B,Get(X,Y,C,D,P[L1]);
		}
		Sort();
		printf("%.5lf\n",(double)ans);
		for (int i=1;i<=4;i++)
		{
			AP[i].x=dcmp(AP[i].x,0)?AP[i].x:0;
			AP[i].y=dcmp(AP[i].y,0)?AP[i].y:0;
			printf("%.5lf %.5lf\n",(double)AP[i].x,(double)AP[i].y);
		}
	return 0;
}


2218

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<cstring>
#define cl(x) memset(x,0,sizeof(x))
#define eps 1e-6
using namespace std;
typedef long double ld;

inline int dcmp(ld a,ld b)
{
	if (fabs(a-b)<eps) return 0;
	if (a<b) return -1; return 1;
}

struct Point{
	ld x,y;
	bool operator < (const Point &b) const{
		return dcmp(x,b.x)<0 || (dcmp(x,b.x)==0 && dcmp(y,b.y)<0);
	}
}P[1005],H[1005];
ld s[1005];

int n,cnt;

inline ld Dot(Point p1,Point p2,Point p0)
{
	return (p1.x-p0.x)*(p2.x-p0.x)+(p1.y-p0.y)*(p2.y-p0.y);
}

inline ld Cross(Point p1,Point p2,Point p0)
{
	return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
}

inline ld Dist(Point p1,Point p2)
{
	return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}

inline ld Dist(Point p1,Point p2,Point p0)
{
	return fabs(Cross(p1,p2,p0)/Dist(p1,p2));
}

inline ld Area(Point p1,Point p2,Point p0)
{
	return fabs(Cross(p1,p2,p0))/2;
}

inline void Qhull(int l,int r,Point a,Point b)
{
	int i=l-1,j=r+1,x=l;
	for (int k=l;k<=r;k++)
		if (dcmp(s[k],s[x])>0 || (dcmp(s[x],s[k])==0 && P[x]<P[k]))
			x=k;
	Point y=P[x];
	for (int k=l;k<=r;k++)
	{
		s[++i]=Cross(P[k],a,y);
		if (dcmp(s[i],0)>0) swap(P[k],P[i]); else s[i--]=0;
	}
	for (int k=r;k>=l;k--)
	{
		s[--j]=Cross(P[k],y,b);
		if (dcmp(s[j],0)>0) swap(P[k],P[j]); else s[j++]=0;
	}
	if (l<=i) Qhull(l,i,a,y);
	H[++cnt]=y;
	if (j<=r) Qhull(j,r,y,b);
}

int main()
{
	ld ans,tmp,A,B; double xx,yy;
	Point X,Y;
	int L1,L2,L3;
	freopen("t.in","r",stdin);
	freopen("t.out","w",stdout);
	while (1)
	{
		int x=1;
		scanf("%d",&n); if (n==0) break;
		for (int i=1;i<=n;i++)
		{
			scanf("%lf%lf",&xx,&yy);
			P[i].x=xx; P[i].y=yy;
			if (P[i]<P[x]) x=i;
		}
		if (n<=3)
		{
			if (n==1 || n==2) 
				printf("%.4lf\n",0.0);
			else if (n==3)
			{
				ld ret=1e20;
				ret=min(ret,Dist(P[1],P[2])*Dist(P[1],P[2],P[3]));
				ret=min(ret,Dist(P[2],P[3])*Dist(P[2],P[3],P[1]));
				ret=min(ret,Dist(P[1],P[3])*Dist(P[1],P[3],P[2]));
				printf("%.4lf\n",(double)ret);
			}
			cl(P); cl(H); cnt=0; cl(s);
			continue;
		}
		swap(P[1],P[x]);
		H[++cnt]=P[1];
		Qhull(2,n,P[1],P[1]);
		n=cnt;
		for (int i=1;i<=cnt;i++)
			P[i]=H[i];
		Y=H[n]; X=H[1];
		L1=2;
		for (int i=3;i<n;i++)
			if (dcmp(Dist(X,Y,P[i]),Dist(X,Y,P[L1]))>0)
				L1=i;
		L2=2;
		for (int i=3;i<n;i++)
			if (dcmp(Dot(Y,P[i],X)/Dist(X,Y),Dot(Y,P[L2],X)/Dist(X,Y))<0)
				L2=i;
		L3=2;
		for (int i=3;i<n;i++)
			if (dcmp(Dot(X,P[i],Y)/Dist(X,Y),Dot(X,P[L3],Y)/Dist(X,Y))<0)
				L3=i;
		B=Dist(X,Y,P[L1]);
		A=Dist(X,Y);
		tmp=Dot(Y,P[L2],X)/Dist(X,Y);
		if (dcmp(tmp,0.0)<0) A-=tmp;
		tmp=Dot(X,P[L3],Y)/Dist(X,Y);
		if (dcmp(tmp,0.0)<0) A-=tmp;
		ans=A*B;
		for (int i=1;i<n;i++)
		{
			Y=P[i]; X=P[i+1];
			while (dcmp(Dist(X,Y,P[L1]),Dist(X,Y,P[L1%n+1]))<0) 
				L1=L1%n+1;
			while (dcmp(Dot(Y,P[L2%n+1],X)/Dist(X,Y),Dot(Y,P[L2],X)/Dist(X,Y))<0)
				L2=L2%n+1;
			while (dcmp(Dot(X,P[L3%n+1],Y)/Dist(X,Y),Dot(X,P[L3],Y)/Dist(X,Y))<0)
				L3=L3%n+1;
			B=Dist(X,Y,P[L1]);
			A=Dist(X,Y);
			tmp=Dot(Y,P[L2],X)/Dist(X,Y);
			if (dcmp(tmp,0.0)<0) A-=tmp;
			tmp=Dot(X,P[L3],Y)/Dist(X,Y);
			if (dcmp(tmp,0.0)<0) A-=tmp;
			ans=min(ans,A*B);
		}
		printf("%.4lf\n",(double)ans);
		cl(P); cl(H); cnt=0; cl(s);
	}
	return 0;
}


里阿奴摩西
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
bzoj1336&1337 最小圆覆盖
AaronPolaris
04-27 4481
随机增量法求最小圆覆盖
wangjunrui666#bzoj-problem##3873. [Ahoi2014&Jsoi2014]拼图1
07-25
作为一个计算机科学家，JYY有一套黑白色的拼图，他希望通过合理的拼接，使得拼出的最终图案中，能包含面积最大的全白色子矩形。之后JYY发现，可以通过改变这S块拼图
BZOJ 1185 HNOI 2007 最小矩形覆盖 旋转卡壳
(╯°口°)╯(┴—┴
01-23 933
题目大意：给出平面上的一些点，问面积最小的矩形满足覆盖所有的点。 思路：覆盖问题和不是凸包上的点没关系，先做凸包。根据贪心的思想，这个覆盖了所有点的矩形肯定至少有一条边与凸包上的边重合，那么我们枚举凸包上的每一条边，对于这个已经确定了一条边的矩形，不难确定其他三个边。注意到已知当前直线的向量，就可以求出两侧和对面的向量，而这三个向量随着枚举的边的移动是单调的，所以就可以用旋转卡壳来卡住剩
BZOJ 1185 HNOI2007 最小矩形覆盖 旋转卡壳
世界
01-23 1710
题目大意：最小矩形覆盖 首先有一个结论：凸包上一定有一条边与矩形的一条边重合 证明：如果不存在一条边与矩形的一条边重合，那么我将这个矩形旋转一下一定会比之前更小 于是我们枚举其中一条边，对其余三个点卡壳即可 这旋转卡壳写的真叫一个卡壳- - 还好1A掉了- - #include #include #include #include #include #define M 5050
bzoj1185: [HNOI2007]最小矩形覆盖 计算几何 旋转卡壳
IED98的专栏
02-26 864
又被卡精度。。。这概率绝了！！！ 我们可以肯定凸包的一边肯定在矩形上，所以我们直接寻找上，左，右三个顶点就可以了。#include #include #include #include #include using namespace std; #define eps 1e-10 #define maxn 51000 struct node { double x,y; }p[m
bzoj1185【HNOI2007】最小矩形覆盖
AaronPolaris
04-19 2999
凸包+旋转卡壳
bzoj 1185: [HNOI2007]最小矩形覆盖 旋转卡壳
lych的博客
04-14 812
先求出凸包，显然至少需要有一条边在凸包上（我只会感性理解不会证明QAQ），那么旋转卡壳一下就好了。 AC代码如下： #include #include #include #include #include #define eps 1e-12 #define N 50005 using namespace std; struct point{ double x,y; }a[N],q[N],t[
BZOJ 1185 [HNOI2007]最小矩形覆盖 旋转卡壳
-Claude
09-09 2423
BZOJ 1185 [HNOI2007]最小矩形覆盖 旋转卡壳
[BZOJ 1185][HNOI 2007]最小矩形覆盖(凸包+旋转卡壳)
某沙茶的OI代码库
03-06 894
题目链接http://www.lydsy.com/JudgeOnline/problem.php?id=1185思路这个题真的很神奇啊。。。虽然有人说被卡精度了，不过反正我是没遇到这样的问题(三态函数大法好) 实际上这个题与我之前写的另一题（最大土地面积）的做法比较相近，也是先维护一个凸包，然后搞旋转卡壳。个人认为此题的旋转卡壳略难一些。 首先我们发现，一个最小覆盖矩形一定是与点集的凸包上的部分
BZOJ 1185([HNOI2007]最小矩形覆盖-旋转卡壳+点集几何意义)
nike0good |Oier&ACMer | 熟能生巧
03-27 2122
1185: [HNOI2007]最小矩形覆盖 Time Limit: 10 Sec  Memory Limit: 162 MBSec  Special Judge Submit: 258  Solved: 137 Description   l要事先改成r，注意把向量的2点设成同一点会出现奇妙的事情 #include #include #include #inc
【BZOJ1185】最小矩形覆盖 计算几何 凸包 旋转卡壳
鲜衣怒马少年郎
04-15 1306
写完程序两小时，调对精度一小时hhh 首先不严格的直觉告诉我们所求矩形一定有一条边在凸包上，我们遍历凸包上的边，用类似旋转卡壳的方式得到凸包上的三个点，分别是在边上投影最靠前的点，在边上投影最靠后的点，距离边最远的点 ，这三个点一定在所求矩形的另外三条边上。得到一个矩形求解即可。 我尽力了。。。 判断1号点和2号点时要同时运用内积（点乘）和外积（叉乘）与0的大小关系，具体见代码 /
wangjunrui666#bzoj-problem##2000. [Hnoi2010]stone 取石头游戏1
07-25
然后两个玩家轮流来取石子，每次每个玩家可以取走一堆中的所有石子，但有一个限制条件：一个玩家若要取走一堆石子，则与这堆石子相邻的某堆石子已被取走(之前被某个玩家取
wangjunrui666#bzoj-problem##3876. [Ahoi2014&Jsoi2014]支线剧情1
07-25
【故事背景】宅男JYY非常喜欢玩RPG游戏，比如仙剑，轩辕剑等等。不过JYY喜欢的并不是战斗场景，而是类似电视剧一般的充满恩怨情仇的剧情。这些游戏往往都有很多的
BZOJ第一部分
11-03
BZOJ原题库中的P1000到P2000题目，覆盖了从基础算法到复杂数据结构和算法的广泛范围。这个阶段的题目难度逐渐提升，旨在帮助用户从基础的编程技巧过渡到更高级的算法应用。这些题目涵盖了动态规划、图论、搜索、排序...
BZOJ全代码
02-05
【标题】"BZOJ全代码"所指的是一份包含BZOJ（BestCoder Zhijiang Online Judge）平台所有编程题目解决方案的压缩文件。BZOJ是一个在线编程竞赛平台，它提供了大量的算法题目供参赛者练习和提交代码，以检验编程能力...
[DP 可并堆维护凸包优化] BZOJ 4585 [Apio2016]烟火表演
雯舞
08-28 1919
垂死梦中惊坐起，膜拜神犇王梦迪 #include #include #include using namespace std; typedef long long ll; inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread
[扫描线 树链剖分 树状数组] BZOJ2758 [SCOI2012]Blinker的噩梦
雯舞
03-18 1436
扫描线  一般扫描线的做法就是想象一根线向某个方向扫过，进入的进入，出去的出去，用一个set维护序列或是结合线段树，当然最恶心的就是加上计算几何 但是落实到这题上就是各种鬼畜乱搞，只能膜大神 省选前做难题的结果是——到处翻题解，程序越改越像 “题意:平面上有n个多边形（凸包和圆）。任意两个多边形AB只有两种关系：（1）A包含B或者B包含A；(2)AB的公共面积为0。每个多边形有一个
[几何 模拟退火 || 随机增量法] Codeforces 442E #253 (Div. 1) E. Gena and Second Distance
雯舞
03-14 1338
怎么都在理性愉悦啊 总算做到一道有点思路的题了官方题解 题解里说的是什么呢 我们先二分答案AnsAns 然后我们考虑AnsAns合法的条件是存在一个圆OO 他只包含其中一个特殊点 我们肯定可以挪动这个圆 使得有至少一个点在边界上 枚举点pp在边界上 那么圆心OO就在以pp为圆心的圆CC上 圆心OO合法的条件是不存在其他点到他的距离 ≤Ans\leq Ans 实际上就是C的圆周上的不被任
[扫描线 计算几何] BZOJ 1845 [Cqoi2005] 三角形面积并
雯舞
03-18 1082
%%%PoPoQQQ：http://blog.csdn.net/popoqqq/article/details/42581881 经典的扫描线 首先求出所有直线交点的横坐标，排序，去重 然后对于每个横坐标，两段之间夹的部分一定是一个或多个梯形 因此我们取中位线，求出中位线被所有三角形覆盖区间的区间并的长度，即可计算出这部分的面积 #includ
[BZOJ] 1631: [Usaco2007 Feb]Cow Party
最新发布
06-08
好的，这是一道关于牛的题目，它描述了一群牛要去参加一个聚会，它们之间有一些道路相连，每条道路都有一个不同的长度。在聚会上，每头牛都要喝一杯牛奶，而且它们的喝牛奶的速度不同，因此有些牛会在其他牛之前到达...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值