(HDUACM08)计算几何基础

ACM程序设计杭州电子科技大学刘春英[e-mail]2020/1/152今天，你吗？2020/1/153每周一星（7）：dwj2020/1/154第八讲计算几何初步(ComputationalGeometryBasic)2020/1/155第一单元线段属性2020/1/1562020/1/1572020/1/1582020/1/1592020/1/15102020/1/15112020/1/15122020/1/15132020/1/1514思考：1、传统的计算线段相交的方法是什么？2、传统方法和本方法的区别是什么？2020/1/1515特别提醒：以上介绍的线段的三个属性，是计算几何的基础，在很多方面都有应用，比如求凸包等等，请务必掌握！2020/1/1516第二单元多边形面积和重心2020/1/1517基本问题（1）：给定一个简单多边形，求其面积。输入：多边形（顶点按逆时针顺序排列）输出：面积S2020/1/1518思考如下图形：2020/1/1519Anygoodidea?2020/1/1520先看最简单的多边形——三角形2020/1/1521三角形的面积：在解析几何里，△ABC的面积可以通过如下方法求得：点坐标=边长=海伦公式=面积2020/1/1522思考：此方法的缺点：计算量大精度损失更好的方法？2020/1/1523计算几何的方法：在计算几何里，我们知道，△ABC的面积就是“向量AB”和“向量AC”两个向量叉积的绝对值的一半。其正负表示三角形顶点是在右手系还是左手系。ABC成左手系，负面积ABC成右手系，正面积BCACBA2020/1/1524大功告成：Area(A,B,C)=1/2*(↑AB)×(↑AC)=∣∣/2特别注意：以上得到是有向面积（有正负）！Xb–XaYb–YaXc–XaYc–Ya2020/1/1525凸多边形的三角形剖分很自然地，我们会想到以P1为扇面中心，连接P1Pi就得到N-2个三角形，由于凸性，保证这些三角形全在多边形内，那么，这个凸多边形的有向面积：A=sigma(Ai)(i=1…N-2)P1P2P3P4P5P6A1A2A3A42020/1/1526凹多边形的面积？P1P4P3P22020/1/1527依然成立！！！多边形面积公式：A=sigma(Ai)(i=1…N-2)结论：“有向面积”A比“面积”S其实更本质！2020/1/1528任意点为扇心的三角形剖分：我们能把多边形分成N-2个三角形，为什么不能分成N个三角形呢？比如，以多边形内部的一个点为扇心，就可以把多边形剖分成N个三角形。P0P1P2P6P5P4P32020/1/1529前面的三角剖分显然对于多边形内部任意一点都是合适的！我们可以得到：A=sigma(Ai)（i=1…N）即：A=sigma∣∣/2（i=1…N）Xi–X0Yi–Y0X(i+1)–X0Y(i+1)–Y02020/1/1530能否把扇心移到多边形以外呢？P0P1P2P3P42020/1/1531既然内外都可以，为什么不设P0为坐标原点呢？OP1P2P3P4现在的公式？2020/1/1532简化的公式：A=sigma∣∣/2（i=1…N）XiYiX(i+1)Y(i+1)面积问题搞定！2020/1/1533基本问题（2）：给定一个简单多边形，求其重心。输入：多边形（顶点按逆时针顺序排列）输出：重心点C2020/1/1534从三角形的重心谈起：三角形的重心是：(x1+x2+x3)/3，(y1+y2+y3)/3可以推广否？Sigma(xi)/N,sigma(yi)/N(i=1…N)???2020/1/1535看看一个特例：2020/1/1536原因：错误的推广公式是“质点系重心公式”，即如果认为多边形的质量仅分布在其顶点上，且均匀分布，则这个公式是对的。但是，现在多边形的质量是均匀分布在其内部区域上的，也就是说，是与面积有关的！2020/1/1537Solution:剖分成N个三角形，分别求出其重心和面积，这时可以想象，原来质量均匀分布在内部区域上，而现在质量仅仅分布在这N个重心点上（等假变换），这时候就可以利用刚才的质点系重心公式了。不过，要稍微改一改，改成加权平均数，因为质量不是均匀分布的，每个质点代表其所在三角形，其质量就是该三角形的面积（有向面积！），——这就是权！2020/1/1538公式：C=sigma(Ai*Ci)/A(i=1…N)Ci=Centroid(△OPiPi+1)=(O+↑Pi+↑Pi+1)/3C=sigma((↑Pi+↑Pi+1)(↑Pi×↑Pi+1))/(6A)全部搞定！2020/1/1540第三单元凸包(ConvexHull)2020/1/15412020/1/15422020/1/15432020/1/15442020/1/15452020/1/15462020/1/15472020/1/15482020/1/15492020/1/15502020/1/15512020/1/15522020/1/15532020/1/15542020/1/15552020/1/15562020/1/15572020/1/15582020/1/15592020/1/15602020/1/15612020/1/15622020/1/15632020/1/15642020/1/15652020/1/15662020/1/1567Anyquestion?2020/1/1569课后在线作业：201403《ACM程序设计》作业（8）——刘春英老师2020/1/1570WelcometoHDOJThankYou~