Andrew 算法可以在 \(O(n\log n)\) 的时间复杂度通过单调栈分别求出散点的上凸壳和下凸壳,来求出平面上一些点的凸包。
看懂这篇博客,大家需要掌握:
本文中的向量恕不加 \(\overrightarrow{}\) 符号。
凸多边形是指所有内角大小都在 \([0,\pi]\) (弧度制)范围内的 简单多边形。其他的“凸”请类比理解。
首先,什么是凸包?
给你平面上的点集,你需要从中选出最少的点,使得这些点所组成的 凸多边形 可以包裹住其他所有点。这些点所组成的凸多边形就是凸包。
譬如下面这个点集:

它的凸包是:

下面我将会告诉大家怎么求。
Andrew 算法需要先对所有点按照 \(x\) 坐标为第一关键字、\(y\) 坐标为第二关键字排序。如上面的点集,经过排序后是:
ABFEDCGJHILMNKO
那么 \(A\) 和 \(O\) 一定在凸包上,因为它们无法被其他点所组成的凸多边形覆盖。
按照 Andrew 算法的逻辑,我们需要先求出凸包的一半 “凸壳”。下面将会以上凸壳为例,下凸壳与其类似。
一段上凸壳一定满足顺时针遍历时,每个节点在每条边所组成的向量的右边(下凸壳在左边)(就是凸包的“凸”,下同)。这句话大家可能不能完全理解,不过没有关系,我会给大家慢慢道来。
首先,按照排序后的点集遍历点集,第一个遍历到的是 \(B\)(\(A\) 不考虑)。我们可以连接 \(AB\):

然后下一个点是 \(F\),继续连接 \(BF\):

下一个点是 \(E\),继续连接 \(FE\):

下一个点是 \(D\),继续连接 \(ED\):

但是这样子我们遇到了问题,\(D\) 在 \(FE\) 左侧,它不凸了,我们的解决办法是:
断掉以前连的边,直到遇到可以连接的点,满足凸壳性质
我们可以断掉 \(ED,FE\),连接 \(FD\),发现还是不满足。

我们继续,断掉 \(FD,BF\),连接 \(BD\),这回满足了。

下一个点是 \(C\),继续连接 \(DC\):

发现又不凸了,我们断掉 \(DC,BD\) 连接 \(BC\),就可以满足了:

下一个点是 \(G\),继续连接 \(CG\):

发现不凸,我们断掉 \(CG,BC\),连接 \(BG\):

下一个点是 \(J\),继续连接 \(GJ\):

下一个点是 \(H\),继续连接 \(JH\):

发现不凸,我们断掉 \(GJ,JH\),连接 \(GH\):

下一个点是 \(I\),继续连接 \(HI\):

下一个点是 \(L\),继续连接 \(IL\):

发现不凸,我们断掉 \(IL,HI\),连接 \(HL\):

发现不凸,我们断掉 \(HL,GH\),连接 \(GL\):

发现不凸,我们断掉 \(GL,BG\),连接 \(BL\):

下一个点是 \(M\),继续连接 \(LM\):

下一个点是 \(N\),继续连接 \(MN\):

发现不凸,我们断掉 \(MN,LM\),连接 \(LN\):

下一个点是 \(K\),继续连接 \(NK\):

发现不凸,我们断掉 \(LN,NK\),连接 \(LK\):

最后一个点是 \(O\),我们连接 \(KO\):

这样子上凸壳便求出来,下凸壳我们一般从 \(O\) 遍历到 \(A\),按照以前的逻辑做即可,最后结果如下:

维护“不凸就断边”我们使用单调栈,如果不满足凸的性质就弹栈,最后入栈即可。注意我们不需要模拟断边操作,只需要将点删除即可。
还有,如何判断是否在左边呢?我们可以使用叉乘的右手定则:

参考代码如下:
int stk[100005];
bool used[100005];
vector<Point> ConvexHull(Point* poly, int n){ // Andrew算法求凸包
int top=0;
sort(poly+1,poly+n+1,[&](Point x,Point y){
return (x.x==y.x)?(x.y<y.y):(x.x<y.x);
});
stk[++top]=1;
for(int i=2;i<=n;i++){
while(top>1&&dcmp((poly[stk[top]]-poly[stk[top-1]])*(poly[i]-poly[stk[top]]))<=0){
used[stk[top--]]=0;
}
used[i]=1;
stk[++top]=i;
}
int tmp=top;
for(int i=n-1;i;i--){
if(used[i]) continue;
while(top>tmp&&dcmp((poly[stk[top]]-poly[stk[top-1]])*(poly[i]-poly[stk[top]]))<=0){
used[stk[top--]]=0;
}
used[i]=1;
stk[++top]=i;
}
vector<Point> a;
for(int i=1;i<=top;i++){
a.push_back(poly[stk[i]]);
}
return a;
}
目录一.加解密算法数字签名对称加密DES(DataEncryptionStandard)3DES(TripleDES)AES(AdvancedEncryptionStandard)RSA加密法DSA(DigitalSignatureAlgorithm)ECC(EllipticCurvesCryptography)非对称加密签名与加密过程非对称加密的应用对称加密与非对称加密的结合二.数字证书图解一.加解密算法加密简单而言就是通过一种算法将明文信息转换成密文信息,信息的的接收方能够通过密钥对密文信息进行解密获得明文信息的过程。根据加解密的密钥是否相同,算法可以分为对称加密、非对称加密、对称加密和非
1.问题描述使用Python的turtle(海龟绘图)模块提供的函数绘制直线。2.问题分析一幅复杂的图形通常都可以由点、直线、三角形、矩形、平行四边形、圆、椭圆和圆弧等基本图形组成。其中的三角形、矩形、平行四边形又可以由直线组成,而直线又是由两个点确定的。我们使用Python的turtle模块所提供的函数来绘制直线。在使用之前我们先介绍一下turtle模块的相关知识点。turtle模块提供面向对象和面向过程两种形式的海龟绘图基本组件。面向对象的接口类如下:1)TurtleScreen类:定义图形窗口作为绘图海龟的运动场。它的构造器需要一个tkinter.Canvas或ScrolledCanva
我一直在尝试用Ruby实现Luhn算法。我一直在执行以下步骤:该公式根据其包含的校验位验证数字,该校验位通常附加到部分帐号以生成完整帐号。此帐号必须通过以下测试:从最右边的校验位开始向左移动,每第二个数字的值加倍。将乘积的数字(例如,10=1+0=1、14=1+4=5)与原始数字的未加倍数字相加。如果总模10等于0(如果总和以零结尾),则根据Luhn公式该数字有效;否则无效。http://en.wikipedia.org/wiki/Luhn_algorithm这是我想出的:defvalidCreditCard(cardNumber)sum=0nums=cardNumber.to_s.s
下面是我写的一个计算斐波那契数列中的值的方法:deffib(n)ifn==0return0endifn==1return1endifn>=2returnfib(n-1)+(fib(n-2))endend它工作到n=14,但在那之后我收到一条消息说程序响应时间太长(我正在使用repl.it)。有人知道为什么会这样吗? 最佳答案 Naivefibonacci进行了大量的重复计算-在fib(14)fib(4)中计算了很多次。您可以将内存添加到您的算法中以使其更快:deffib(n,memo={})ifn==0||n==1returnnen
为了防止在迁移到生产站点期间出现数据库事务错误,我们遵循了https://github.com/LendingHome/zero_downtime_migrations中列出的建议。(具体由https://robots.thoughtbot.com/how-to-create-postgres-indexes-concurrently-in概述),但在特别大的表上创建索引期间,即使是索引创建的“并发”方法也会锁定表并导致该表上的任何ActiveRecord创建或更新导致各自的事务失败有PG::InFailedSqlTransaction异常。下面是我们运行Rails4.2(使用Acti
RGeo有可用的凸包方法,但在这方面根本没有文档。给定一组点,我如何找到它们的凸包? 最佳答案 好问题。事实证明,有一种名为“MultiPoint”的几何类型适用于此。我做了一个简单的例子来测试它,它似乎运行良好。f=RGeo::Geos.factory(:srid=>3361,:buffer_resolution=>8)#mytypicallocalrectilinearprojectionfactorywithmydefaultsettings.coords=[[1,1],[2,2],[1,3]]points=[]coords.
我正在开发一个类似微论坛的项目,其中一个特殊用户发布一条快速(接近推文大小)的主题消息,订阅者可以用他们自己的类似大小的消息来响应。直截了当,没有任何形式的“挖掘”或投票,只是每个主题消息的响应按时间顺序排列。但预计会有很高的流量。我们想根据它们引起的响应嗡嗡声来标记主题消息,使用0到10的等级。在谷歌上搜索了一段时间的趋势算法和开源社区应用示例,到目前为止已经收集到两个有趣的引用资料,但我还没有完全理解它们:Understandingalgorithmsformeasuringtrends,关于使用基线趋势算法比较维基百科页面浏览量的讨论,在SO上。TheBritneySpearsP
我收到错误:unsupportedcipheralgorithm(AES-256-GCM)(RuntimeError)但我似乎具备所有要求:ruby版本:$ruby--versionruby2.1.2p95OpenSSL会列出gcm:$opensslenc-help2>&1|grepgcm-aes-128-ecb-aes-128-gcm-aes-128-ofb-aes-192-ecb-aes-192-gcm-aes-192-ofb-aes-256-ecb-aes-256-gcm-aes-256-ofbRuby解释器:$irb2.1.2:001>require'openssl';puts
文章目录一.Dijkstra算法想解决的问题二.Dijkstra算法理论三.java代码实现一.Dijkstra算法想解决的问题解决的问题:求解单源最短路径,即各个节点到达源点的最短路径或权值考察其他所有节点到源点的最短路径和长度局限性:无法解决权值为负数的情况二.Dijkstra算法理论参数:S记录当前已经处理过的源点到最短节点U记录还未处理的节点dist[]记录各个节点到起始节点的最短权值path[]记录各个节点的上一级节点(用来联系该节点到起始节点的路径)Dijkstra算法步骤:(1)初始化:顶点集S:节点A到自已的最短路径长度为0。只包含源点,即S={A}顶点集U:包含除A外的其他顶
对于体育新闻中文文本的关键字提取,常用的算法包括TF-IDF、TextRank和LDA等。它们的基本步骤如下:1.TF-IDF算法: -将文本进行分词和词性标注处理。-统计每个词在文本中的词频(TF)。-计算每个词在整个语料库中出现的文档频率(DF)和逆文档频率(IDF)。-计算每个词的TF-IDF值,并按照值的大小进行排序,选择排名前几的词作为关键字。2.TextRank算法:-将文本进行分词和词性标注处理。-将分词结果转化成图模型,每个词语为节点,根据词语之间的共现关系建立边。-对图模型进行迭代计算,计算每个节点的PageRank值,表示该节点的重要性。-选择排名前几的节点作为关键字。3.