计算几何（半平面交&&旋转卡壳）

最新推荐文章于 2020-02-08 16:11:30 发布
最新推荐文章于 2020-02-08 16:11:30 发布
阅读量147 收藏
点赞数
原文链接：http://www.cnblogs.com/LitIce/archive/2011/04/09/2010180.html
版权

先介绍个写的很好的blog ————Master_Chivu

[Poj 1113] 计算几何之凸包(一) {卷包裹算法}

http://www.cnblogs.com/Booble/archive/2011/02/28/1967179.html

http://www.cnblogs.com/Booble/archive/2011/03/10/1980089.html

[Poj 2187]计算几何之凸包(三) {旋转卡壳初步}

http://www.cnblogs.com/Booble/archive/2011/04/03/2004865.html

旋转卡壳是在上面blog里学到的。

半平面交是在http://apps.hi.baidu.com/share/detail/14889384学到的。

下面是我的模板。

ContractedBlock.gif ExpandedBlockStart.gif 2451（Halfe_Plane_Intersection
ca...POJ的RE提示WA。害我看了两天。以为自己模板又挂了 - - 。
===================================================================
if ( fabs(seg[i].angle - seg[i - 1 ].angle) > eps )
{

 while (deq.size() >= 2 && cross(seg[i].sou,seg[i].tar,inter( * deq.begin() , * (deq.begin() + 1 ) )) < - eps)
deq.pop_front();
 while (deq.size() >= 2 && cross(seg[i].sou,seg[i].tar,inter( * (deq.end() - 1 ), * (deq.end() - 2 ) ) ) < - eps )
deq.pop_back();
deq.push_front(seg[i]);
}
 while (deq.size() > 2 && cross((deq.end() - 1 ) -> sou,(deq.end() - 1 ) -> tar,
inter( * (deq.begin()), * (deq.begin() + 1 ) ) ) < - eps)
deq.pop_front();
 while (deq.size() > 2 && cross( deq.begin() -> sou, deq.begin() -> tar,
inter( * (deq.end() - 1 ), * (deq.end() - 2 ) ) ) < - eps )
deq.pop_back();

一直怀疑划线加粗句是废话，固然 - - 。去掉后能能毫无压力过掉。
 ====================================================================
#include < iostream >
#include < cstdio >
#include < deque >
#include < cmath >
#include < algorithm >
#include < cstdlib >
using namespace std;
 const int N = 300000 ;
 const int M = 300000 ;
 const double eps = 1e - 6 ;
 struct node
{
 double x,y;
}poi[N],tmp1,tmp2,zero,tpoi[N];
 struct SEG
{
node sou,tar;
 double angle;
}seg[M];
deque < SEG > deq;
 int Test,n,cnt,flag;
 double tmp,s;
 double dot(node a,node b)
{
 return a.x * b.x + a.y * b.y;
}
 double cross(node a,node b,node c)
{
 return (b.x - a.x) * (c.y - a.y) - (b.y - a.y) * (c.x - a.x);
}
node inter( SEG a, SEG b)
{
 double s1,s2,ts;
s1 = cross(a.sou,a.tar,b.sou);
s2 = cross(a.sou,a.tar,b.tar);
ts = s1 - s2;
 if (fabs(ts) < eps) { tmp1.x = 10001 ; tmp1.y = 10001 ; return tmp1; }
tmp1.x = b.sou.x * ( - s2) / ts + b.tar.x * s1 / ts;
tmp1.y = b.sou.y * ( - s2) / ts + b.tar.y * s1 / ts;
 return tmp1;
}
 bool cmp( const SEG & a, const SEG & b)
{
 if ( fabs(a.angle - b.angle) < eps )
 return cross(a.sou,a.tar,b.sou) < 0 ;
 else
return a.angle < b.angle;

}
 void Half_Plane_Intersection()
{
deq.clear();
deq.push_front(seg[ 0 ]);
 int tt = 1 ; while (fabs(seg[tt].angle - seg[tt - 1 ].angle) < eps) ++ tt;
deq.push_front(seg[tt]);
 for ( int i = tt + 1 ;i < n; ++ i)
 if ( fabs(seg[i].angle - seg[i - 1 ].angle) > eps )
{

 while (deq.size() >= 2 && cross(seg[i].sou,seg[i].tar,inter( * deq.begin() , * (deq.begin() + 1 ) )) < - eps)
deq.pop_front();
 while (deq.size() >= 2 && cross(seg[i].sou,seg[i].tar,inter( * (deq.end() - 1 ), * (deq.end() - 2 ) ) ) < - eps )
deq.pop_back();
deq.push_front(seg[i]);
}
 while (deq.size() > 2 && cross((deq.end() - 1 ) -> sou,(deq.end() - 1 ) -> tar,
inter( * (deq.begin()), * (deq.begin() + 1 ) ) ) < - eps)
deq.pop_front();
 while (deq.size() > 2 && cross( deq.begin() -> sou, deq.begin() -> tar,
inter( * (deq.end() - 1 ), * (deq.end() - 2 ) ) ) < - eps )
deq.pop_back();
cnt = 0 ;
 for ( deque < SEG > ::iterator it = deq.begin(); it != deq.end(); ++ it)
{
deque < SEG > ::iterator itt = it + 1 ;
 if (itt == deq.end() ) itt = deq.begin();
tpoi[cnt ++ ] = inter( * it , * itt );
}
s = 0 ;
for ( int i = 0 ;i < cnt; ++ i)
s += tpoi[i].x * tpoi[(i + 1 ) % cnt].y - tpoi[(i + 1 ) % cnt].x * tpoi[i].y;
 if ( s == 0 && cnt <= 2 ) puts( " 0.0 " );
 else printf( " %.1lf\n " ,s /- 2 + eps);
}
 void work()
{
cin >> n;
 for ( int i = 0 ;i < n; ++ i)
cin >> seg[i].sou.x >> seg[i].sou.y >> seg[i].tar.x >> seg[i].tar.y;
seg[n].sou.x = 0 ; seg[n].sou.y = 0 ; seg[n].tar.x = 10000 ; seg[n].tar.y = 0 ; ++ n;
seg[n].sou.x = 10000 ; seg[n].sou.y = 0 ; seg[n].tar.x = 10000 ; seg[n].tar.y = 10000 ; ++ n;
seg[n].sou.x = 10000 ; seg[n].sou.y = 10000 ; seg[n].tar.x = 0 ; seg[n].tar.y = 10000 ; ++ n;
seg[n].sou.x = 0 ; seg[n].sou.y = 10000 ; seg[n].tar.x = 0 ; seg[n].tar.y = 0 ; ++ n;
 for ( int i = 0 ;i < n; ++ i)
seg[i].angle = atan2(seg[i].tar.y - seg[i].sou.y,seg[i].tar.x - seg[i].sou.x);
sort(seg,seg + n,cmp);
Half_Plane_Intersection();
}
 int main()
{
work();
 return 0 ;
}

下面是旋转卡壳模板

ContractedBlock.gif ExpandedBlockStart.gifPOJ3608  旋转卡壳
给你两个不想交的凸包 ~ 求两个凸包间的最短距离。
旋转卡壳模板题。
 ====================================================================================
int change_angle( int tp1, int tp2, double & ang)
{
 double ang1,ang2,a1,a2;
ang1 = atan2( HullA[ (tp1 + 1 ) % cnta ].y - HullA[ tp1 ].y , HullA[ (tp1 + 1 ) % cnta ].x - HullA[ tp1 ].x);
 if (ang1 < 0 ) ang1 = 2 * pi + ang1;
a1 = ang1 - ang; if (a1 < 0 ) a1 += 2 * pi;
ang2 = atan2( - (HullB[ (tp2 + 1 ) % cnta ].y - HullB[ tp2 ].y) , - (HullB[ (tp2 + 1 ) % cntb ].x - HullB[ tp2 ].x));
 if (ang2 < 0 ) ang2 = 2 * pi + ang2;
a2 = ang2 - ang; if (a2 < 0 ) a2 += 2 * pi;
 if ( fabs( a1 - a2 ) < eps ){ ang = ang1; return 2 ;}
 if ( a1 < a2 ) {ang = ang1; return 0 ;} else {ang = ang2; return 1 ;}
}
计算平行线下一次跟哪边相贴。
返回1跟A ，返回 2 跟B相贴，返回3，同时相贴。

 =====================================================================================
while (flag < (n + m) * 5 )
{
tmp = change_angle(tp1,tp2,ang);
ans = min( ans , dis(HullB[tp2],HullB[(tp2 + 1 ) % cntb],HullA[tp1]) );
ans = min( ans , dis(HullB[tp2],HullB[(tp2 + 1 ) % cntb],HullA[(tp1 + 1 ) % cnta]) );
ans = min( ans ,dis(HullA[tp1], HullA[(tp1 + 1 ) % cnta], HullB[tp2] ) );
ans = min( ans ,dis(HullA[tp1], HullA[(tp1 + 1 ) % cnta], HullB[(tp2 + 1 ) % cntb] ) );
 if (tmp == 0 )
tp1 = (tp1 + 1 ) % cnta;
 else if (tmp == 1 )
tp2 = (tp2 + 1 ) % cntb;
 else
{
tp1 = (tp1 + 1 ) % cnta;
tp2 = (tp2 + 1 ) % cntb;
}
 ++ flag;
}
旋转卡壳主过程 ~ 每次判断当前两条平行线所在凸包边的四个点两两之间距离。 旋转之 ~
不要偷懒，还是每次都比较四个点比较保险 - - 。通常不会卡这点常数 ~
旋转次数嘛 ~ 多点就好了。
angle转360就可以了吧。偷懒下。。。
 ============================================================================================


#include < iostream >
#include < cstdio >
#include < cmath >
#include < algorithm >
using namespace std;
 const double pi = acos( - 1.0 );
 double ans;
 const int N = 100000 ;
 const double eps = 1e - 8 ;
 int n,cnt,cnta,cntb,tp,m,tmp,tp1,tp2;
 struct node
{
 double x,y;
} HullA[N],p[N],HullB[N],Hull[N];
 double cross(node a,node b, node c)
{
 return (b.x - a.x) * (c.y - a.y) - (b.y - a.y) * (c.x - a.x);
}
 bool cmp( const node & a , const node & b)
{
 if (a.y == b.y ) return a.x < b.x;
 return a.y < b.y;
}
 double dis(node a,node b)
{
 return sqrt( (a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y) );
}
 double dot(node a,node b,node c)
{
 return (b.x - a.x) * (c.x - a.x) + (b.y - a.y) * (c.y - a.y);
}
 double dis(node a,node b, node c)
{
 if (dot(a,b,c) < 0 || dot(b,a,c) < 0 ) return min(dis(a,c) , dis(b,c) );
 return fabs( cross(a,b,c) ) / dis(a,b);
}
 void Calc_Hull( int n)
{
 for ( int i = 0 ;i < n; ++ i)
cin >> p[i].x >> p[i].y;
sort(p,p + n,cmp);
cnt = 0 ;
 for ( int i = 0 ;i < n; ++ i)
{
 while (cnt > 1 && cross(Hull[cnt - 2 ],Hull[cnt - 1 ],p[i]) < eps ) cnt -- ;
Hull[cnt ++ ] = p[i];
}
int t = cnt;
for ( int i = n - 1 ;i >= 0 ;i -- )
{
 while (cnt > t && cross(Hull[cnt - 2 ],Hull[cnt - 1 ],p[i]) < eps) cnt -- ;
Hull[cnt ++ ] = p[i];
}
 -- cnt;
}
 int change_angle( int tp1, int tp2, double & ang)
{
 double ang1,ang2,a1,a2;
ang1 = atan2( HullA[ (tp1 + 1 ) % cnta ].y - HullA[ tp1 ].y , HullA[ (tp1 + 1 ) % cnta ].x - HullA[ tp1 ].x);
 if (ang1 < 0 ) ang1 = 2 * pi + ang1;
a1 = ang1 - ang; if (a1 < 0 ) a1 += 2 * pi;
ang2 = atan2( - (HullB[ (tp2 + 1 ) % cnta ].y - HullB[ tp2 ].y) , - (HullB[ (tp2 + 1 ) % cntb ].x - HullB[ tp2 ].x));
 if (ang2 < 0 ) ang2 = 2 * pi + ang2;
a2 = ang2 - ang; if (a2 < 0 ) a2 += 2 * pi;
 if ( fabs( a1 - a2 ) < eps ){ ang = ang1; return 2 ;}
 if ( a1 < a2 ) {ang = ang1; return 0 ;} else {ang = ang2; return 1 ;}
}
 int main()
{
 int Test = 0 ;
 while (cin >> n >> m)
{
 if (n == 0 ) break ;
Calc_Hull(n);
cnta = cnt;
 for ( int i = 0 ;i < cnta; ++ i) HullA[i] = Hull[i];
Calc_Hull(m);
cntb = cnt;
 for ( int i = 0 ;i < cntb; ++ i) HullB[i] = Hull[i];
tp1 = 0 ;tp2 = 0 ;
 for ( int i = 1 ;i < cntb; ++ i) if ( ! cmp(HullB[i],HullB[tp2]) ) tp2 = i;
ans = dis( HullA[ 0 ], HullB[tp2] );
 int flag = 0 ;
 double ang = 0 ;
 while (flag < (n + m) * 5 )
{
tmp = change_angle(tp1,tp2,ang);
ans = min( ans , dis(HullB[tp2],HullB[(tp2 + 1 ) % cntb],HullA[tp1]) );
ans = min( ans , dis(HullB[tp2],HullB[(tp2 + 1 ) % cntb],HullA[(tp1 + 1 ) % cnta]) );
ans = min( ans ,dis(HullA[tp1], HullA[(tp1 + 1 ) % cnta], HullB[tp2] ) );
ans = min( ans ,dis(HullA[tp1], HullA[(tp1 + 1 ) % cnta], HullB[(tp2 + 1 ) % cntb] ) );
 if (tmp == 0 )
tp1 = (tp1 + 1 ) % cnta;
 else if (tmp == 1 )
tp2 = (tp2 + 1 ) % cntb;
 else
{
tp1 = (tp1 + 1 ) % cnta;
tp2 = (tp2 + 1 ) % cntb;
}
 ++ flag;
}
printf( " %.5lf\n " ,ans);
}
 return 0 ;
}

转载于:https://www.cnblogs.com/LitIce/archive/2011/04/09/2010180.html

alhe29828
关注
【详解】半平面交算法入门详解(计算几何)
Winsoul Blog
10-30 3万+
半平面交 简介 博客背景 笔者在学习半平面交时，网上找入门博客资源甚少，且大部分难以理解，故在稍稍入门了半平面交后，写此博客，希望能对大家有所帮助。若有错误，麻烦指出。 半平面交是什么？ 我们知道一条直线可以把平面分为两部分，其中一半的平面就叫半平面。 那半平面交，就是多个半平面的相交部分。我们在学习线性规划时就有用过。 半平面交有什么用？ 1.求解一个区域，可以看到给定图形的各个角落。(多边形的...
【计算几何全家桶】计算几何二维板子套路详解
SC_Linno的博客
11-16 814
title : 计算几何（二维） tags : ACM,数学,计算几何 date : 2021-10-23 author : Linno 向量 点乘与叉乘 叉乘一般用于求面积或者判断两点的相对位置：(x1,y1)×(x2,y2)=x1*y2-x2*y1; 点乘一般用于判断直线是否正交：(x1,y1)·(x2,y2)=x1*x2+y1*y2; inline double Len(Vt a){return sqrt(Dot(a,a));}//模长 inline Vt operator +(Vt a,Vt .
半平面交算法及简单应用
weixin_30247307的博客
11-01 218
半平面交算法及简单应用 半平面：一条直线把二维平面分成两个平面。 半平面交：在二维几何平面上，给出若干个半平面，求它们的公共部分 半平面交的结果：1.凸多边形（后面会讲解到）2.无界，因为有可能若干半平面没有形成封闭3.直线，线段，点，空（属于特殊情况吧） 算法：1：根据上图可以知道，运用给出的多边形每相邻两点形成一...
算法复习——半平面交（bzoj2618凸多边形）
banzhangzhi8819的博客
03-25 136
讲解： 这里套用wuvin神犇的ppt，附上友情链接：http://blog.leanote.com/wuvin 半平面交： 算法流程： 注意事项： 例题： Description 逆时针给出n个凸多边形的顶点坐标，求它们交的面积。例如n=2时，两个凸多边形如下图： 则相交部分的面积为5.233。 Input 第一行有一个整数...
半平面交详解
baizhang7027的博客
11-22 661
更好的阅读体验 定义： 半平面： 顾名思义，就是平面的一半。一条直线会把平面分成两部分，就是两个半平面。对于半平面，我们可以用直线方程式如：\(ax + by >= c\) 表示，更常用的是用直线表示。 半平面交： 顾名思义，就是多个半平面求交集。其结果可能是一个凸多边形、无穷平面、直线、线段、点等。 多边形的核： 如果多边形中存在一个区域使得在区域中可以看到多边形中任意...
多边形的核（半平面交算法）
soberman的专栏
03-10 2140
//看代码前最好看下zzy的论文！！！！王飞的资料里有 #include#include#include#includeusing namespace std;const double eps = 1e-8;const int MaxN = 50005;typedef struct{double x, y;}Point;typedef struct{do
计算几何：深入探索半平面交算法
为了实现完整的半平面交算法，通常可以采用扫线法（Sweep Line Algorithm）、旋转卡壳法（Rotating Calipers）或者最近点对方法等经典算法。每种算法都有其优缺点，选择哪种取决于具体的需求和性能要求。在实践中，...
半平面交的新算法及其实用价值_朱泽园.ppt
09-14
半平面交问题是一个重要的计算机科学领域内的几何计算问题，它涉及到如何找到一组直线（或半平面）在平面上的公共部分。这个问题在图形学、路径规划、数据结构和算法设计等多个方面都有广泛应用。朱泽园的论文《半...
最小圆覆盖 半平面交 凸包 三维几何
05-31
在计算几何中，半平面交经常用于构造复杂的形状，例如多边形。 凸包是包含一个集合内所有点且具有最少顶点的多边形。在二维空间中，可以使用Gift Wrapping算法（也称为Jarvis March）或者 Graham Scan来计算。凸包...
二维几何常用算法合集 繁凡的《计算几何全家桶（二）》1
08-03
半平面交是指两个或多个半平面（由直线分割而成的两个区域之一）的交集。有向直线是定义半平面的关键，它可以指定一个方向，使一个半平面位于直线的一侧，另一个在其另一侧。O(nlogn)的半平面交算法，如Graham's ...
POJ 3335 半平面相交算法
zhenzi_1989的专栏
09-17 587
题意是问在多边形中是不是存在一个点，这个点可以看到多边形内部的任意一个地方。 用的是半平面相交的模板，调试了半天，最后发现竟然是m和n写错了，悲剧~   #include #include #include #include #include #include
半平面交
commonc的博客
02-17 3107
算法目的：求NN个半平面的交定义：一个半平面指的是由满足ax+by+c>0ax+by+c>0或ax+by+c>=0ax+by+c>=0的点集组成的二维区域。一般来说在写代码的时候，我们可以把一个半平面想象成一个向量所在的直线右面的一片区域 半平面的交可能是一个封闭图形也可能是没有边界的区域或者为空大概思路：先将所有的向量按照与x轴的夹角排序(此处的夹角在[−π，π][-π，π]这个区间内)，然后按
计算几何 || 半平面交 （模板收集）
野生
05-07 874
#include #include #include #include using namespace std; const double eps = 1e-8; struct Point{ double x,y; Point(double xx=0.0,double yy=0.0):x(xx),y(yy){} Point operator - (const Point &
计算几何半平面交入门详解
slass-mot-odet
08-08 826
讲这个问题之前我们必须明白一个问题，何为半平面交？ 我们知道一条直线可以把平面分为两部分，其中一半的平面就叫半平面。 那半平面交，就是多个半平面的相交部分。 半平面交解决的问题是什么？ 通过半平面交我们可以找到多边形的核（在核区域内能看见多边形边上的任意一个点） 我们也能通过半平面交找到多边形内能放的最大圆 标题半平面交有两种实现方法： 第一种：遍历所有线段，每次保留线段左边的区域，删除线段右边的...
【kuangbin专题】计算几何_半平面交
weixin_33901926的博客
10-15 124
1.poj3335Rotating Scoreboard 传送：http://poj.org/problem?id=3335 题意：就是有个球场，球场的形状是个凸多边形，然后观众是坐在多边形的边上的，问你是否在球场上有个地方可以放一个记分牌，然后所有的观众都可以看得到的。 分析：多边形是否存在内核。 1 #include<iostream> 2 #include...
计算几何基础--半平面求交
永远鲜红の幼月的博客
02-24 1386
半平面，平面的一半。半平面求交，就是求得n个半平面的相交的平面。 一开始可能不太好理解。举个例子： 一条线，有两个半平面          两条线，四个半平面，有一个相     三条线，7个半平面，然后有一个相交的区域，这个区域就是半平面                                                   交，及存在蓝色和绿色的区域         相交的...
【计算几何】【模板】半平面交
ygmjsjdboy的博客
02-08 136
Luogu P4196 半平面交（S&I算法） 求多条直线左侧包围图形 逆时针确定每条直线（可用叉积求面积正负） 极角排序（将极角相同的，按左右排） 每个极角只保留一个最左边的 用双端队列维护： 判断队尾和（队尾-1）的交点是否在新直线的右侧 判断队首和（队首+1）的交点是否在新直线的右侧 先判交点个数，>2>2>2才会有交，否则为0 叉积求面积 #include&l...
半平面交讲解 - N^2写法
08-05 2747
ps:如果这是第一次接触半平面交,不要认为它是很个高深的东东.其实以你现在的几何水平完全能够YY写出它. 半平面: 对一条直线,它将整个平面分成了两部分,对于其中的每一部分.就叫做一个半平面.(这个应该好理解吧) 现在我们要保留其中一个半平面,  可以用向量的左边,右边来标记,  还可以用  Ax+By+C 与 0的关系来标记(这里你千万别纠结它是在直线上方呢,还是下方,而让自
半平面交——学习笔记
Wflower的Blog //What are you waiting for?
07-11 969
这次讲半平面交。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值