Computational Geometry-Theory and Applications最新文献

{"title":"On path-greedy geometric spanners","authors":"William Evans ,&nbsp;Lucca Morais de Arruda Siaudzionis","doi":"10.1016/j.comgeo.2022.101948","DOIUrl":"https://doi.org/10.1016/j.comgeo.2022.101948","url":null,"abstract":"<div><p>A <em>t</em>-spanner is a subgraph of a graph <em>G</em> in which the length of the shortest path between two vertices never exceeds <em>t</em> times the length of the shortest path between them in <em>G</em><span>. A geometric graph is one whose vertices are points and whose edges are line segments between the corresponding points. Geometric </span><em>t</em>-spanners are <em>t</em><span>-spanners of the complete geometric graph on a given point set. Besides approximating the distance between points, we may ask a geometric </span><em>t</em>-spanner to be planar, have low degree, or low total edge length.</p><p>One famous algorithm used to generate spanners is <em>path-greedy</em>, which scans pairs of vertices in non-decreasing order of edge length and adds the edge between them unless the current set of added edges already connects them with a path that <em>t</em>-approximates the edge length. Graphs from this algorithm are called <em>path-greedy spanners</em><span>. This work analyzes properties of path-greedy geometric spanners under different conditions. Specifically, we answer an open problem regarding the planarity and degree of path-greedy </span><em>t</em>-spanners for points in convex position in 2D. Further, we show a simple and efficient way to reduce the degree of a geometric spanner by adding extra points.</p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"110 ","pages":"Article 101948"},"PeriodicalIF":0.6,"publicationDate":"2023-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49812087","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"On the recognition and reconstruction of weighted Voronoi diagrams and bisector graphs","authors":"Günther Eder,&nbsp;Martin Held,&nbsp;Stefan de Lorenzo,&nbsp;Peter Palfrader","doi":"10.1016/j.comgeo.2022.101935","DOIUrl":"10.1016/j.comgeo.2022.101935","url":null,"abstract":"<div><p>A weighted bisector graph is a geometric graph whose faces are bounded by edges that are portions of multiplicatively weighted bisectors of pairs of (point) sites such that each of its faces is defined by exactly one site. A prominent example of a bisector graph is the multiplicatively weighted Voronoi diagram of a finite set of points which induces a tessellation of the plane into Voronoi faces bounded by circular arcs and straight-line segments. Several algorithms for computing various types of bisector graphs are known. In this paper we reverse the problem: Given a partition <span><math><mi>G</mi></math></span> of the plane into faces, find a set of points and suitable weights such that <span><math><mi>G</mi></math></span> is a bisector graph of the weighted points, if a solution exists. If <span><math><mi>G</mi></math></span> is a graph that is regular of degree three then we can decide in <span><math><mi>O</mi><mo>(</mo><mi>m</mi><mo>)</mo></math></span> time whether it is a bisector graph, where <em>m</em> denotes the combinatorial complexity of <span><math><mi>G</mi></math></span>. In the same time we can identify up to two candidate solutions such that <span><math><mi>G</mi></math></span> could be their multiplicatively weighted Voronoi diagram. Additionally, we show that it is possible to recognize <span><math><mi>G</mi></math></span> as a multiplicatively weighted Voronoi diagram and find all possible solutions in <span><math><mi>O</mi><mo>(</mo><mi>m</mi><mi>log</mi><mo>⁡</mo><mi>m</mi><mo>)</mo></math></span> time if <span><math><mi>G</mi></math></span> is given by a set of disconnected lines and circles.</p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101935"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41857204","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Geodesic obstacle representation of graphs","authors":"Prosenjit Bose ,&nbsp;Paz Carmi ,&nbsp;Vida Dujmović ,&nbsp;Saeed Mehrabi ,&nbsp;Fabrizio Montecchiani ,&nbsp;Pat Morin ,&nbsp;Luís Fernando Schultz Xavier da Silveira","doi":"10.1016/j.comgeo.2022.101946","DOIUrl":"https://doi.org/10.1016/j.comgeo.2022.101946","url":null,"abstract":"<div><p>An <em>obstacle representation</em> of a graph is a mapping of the vertices onto points in the plane and a set of connected regions of the plane (called <em>obstacles</em>) such that the straight-line segment connecting the points corresponding to two vertices does not intersect any obstacles if and only if the vertices are adjacent in the graph. Recently, Biedl and Mehrabi (Graph Drawing 2017) studied <em>non-blocking grid obstacle representations</em> of graphs in which the vertices of the graph are mapped onto points in the plane while the straight-line segments representing the adjacency between the vertices are replaced by the <span><math><msub><mrow><mi>L</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span> (Manhattan) shortest paths in the plane that avoid obstacles.</p><p>In this paper, we introduce the notion of <em>geodesic obstacle representations</em><span> of graphs with the main goal of providing a generalized model, which comes naturally when viewing line segments as shortest paths in the Euclidean plane. To this end, we extend the definition of obstacle representation by allowing </span><em>some</em> obstacles-avoiding shortest path between the corresponding points in the underlying metric space whenever the vertices are adjacent in the graph. We consider both <em>general</em> and <em>plane</em> variants of geodesic obstacle representations (in a similar sense to obstacle representations) under any polyhedral distance function in <span><math><msup><mrow><mi>R</mi></mrow><mrow><mi>d</mi></mrow></msup></math></span> as well as shortest path distances in graphs. Our results generalize and unify the notions of obstacle representations, plane obstacle representations and grid obstacle representations, leading to a number of questions on such representations.</p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101946"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49790771","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Value-offset bifiltrations for digital images","authors":"Anway De,&nbsp;Thong Vo,&nbsp;Matthew Wright","doi":"10.1016/j.comgeo.2022.101939","DOIUrl":"https://doi.org/10.1016/j.comgeo.2022.101939","url":null,"abstract":"<div><p><span>Persistent homology<span>, an algebraic method for discerning structure in abstract data, relies on the construction of a sequence of nested topological spaces known as a filtration. Two-parameter persistent homology allows the analysis of data simultaneously filtered by two parameters, but requires a bifiltration—a sequence of topological spaces simultaneously indexed by two parameters. To apply two-parameter persistence to digital images, we first must consider bifiltrations constructed from digital images, which have scarcely been studied. We introduce the </span></span><em>value-offset bifiltration</em><span> for grayscale digital image data. We present efficient algorithms for computing this bifiltration with respect to the taxicab distance and for approximating it with respect to the Euclidean distance. We analyze the runtime complexity of our algorithms, demonstrate the results on sample images, and contrast the bifiltrations obtained from real images with those obtained from random noise.</span></p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101939"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49790772","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Experiments with unit disk cover algorithms for covering massive pointsets","authors":"Rachel Friederich,&nbsp;Anirban Ghosh,&nbsp;Matthew Graham,&nbsp;Brian Hicks,&nbsp;Ronald Shevchenko","doi":"10.1016/j.comgeo.2022.101925","DOIUrl":"https://doi.org/10.1016/j.comgeo.2022.101925","url":null,"abstract":"<div><p>Given a set of <em>n</em><span> points in the plane, the Unit Disk Cover (UDC) problem asks to compute the minimum number of unit disks required to cover the points, along with a placement of the disks. The problem is NP-hard and several approximation algorithms have been designed over the last three decades. In this paper, we have engineered and experimentally compared practical performances of some of these algorithms on massive pointsets. The goal is to investigate which algorithms run fast and give good approximation in practice.</span></p><p>We present a simple 7-approximation algorithm for UDC that runs in <span><math><mi>O</mi><mo>(</mo><mi>n</mi><mo>)</mo></math></span> expected time and uses <span><math><mi>O</mi><mo>(</mo><mi>s</mi><mo>)</mo></math></span> extra space, where <em>s</em> denotes the size of the generated cover. In our experiments, it turned out to be the speediest of all. We also present two heuristics to reduce the sizes of covers generated by it without slowing it down by much.</p><p>To our knowledge, this is the first work that experimentally compares algorithms for the UDC problem. Experiments with them using massive pointsets (in the order of millions) throw light on their practical uses. We share the engineered algorithms via <span>GitHub</span><span>¹</span> for broader uses and future research in the domain of geometric optimization.</p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101925"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49790782","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Covering convex polygons by two congruent disks","authors":"Jongmin Choi ,&nbsp;Dahye Jeong ,&nbsp;Hee-Kap Ahn","doi":"10.1016/j.comgeo.2022.101936","DOIUrl":"https://doi.org/10.1016/j.comgeo.2022.101936","url":null,"abstract":"<div><p><span>We consider the planar two-center problem for a convex polygon: given a convex polygon in the plane, find two congruent disks of minimum radius whose union contains the polygon. We present an </span><span><math><mi>O</mi><mo>(</mo><mi>n</mi><mi>log</mi><mo>⁡</mo><mi>n</mi><mo>)</mo></math></span>-time algorithm for the two-center problem for a convex polygon, where <em>n</em> is the number of vertices of the polygon. This improves upon the previous best algorithm for the problem.</p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101936"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49790784","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Compression for 2-parameter persistent homology","authors":"Ulderico Fugacci ,&nbsp;Michael Kerber ,&nbsp;Alexander Rolle","doi":"10.1016/j.comgeo.2022.101940","DOIUrl":"https://doi.org/10.1016/j.comgeo.2022.101940","url":null,"abstract":"<div><p>Compression aims to reduce the size of an input, while maintaining its relevant properties. For multi-parameter persistent homology, compression is a necessary step in any computational pipeline, since standard constructions lead to large inputs, and computational tasks in this area tend to be expensive. We propose two compression methods for chain complexes of free 2-parameter persistence modules. The first method extends the multi-chunk algorithm for one-parameter persistent homology, returning the smallest chain complex among all the ones quasi-isomorphic to the input. The second method produces minimal presentations of the homology of the input; it is based on an algorithm of Lesnick and Wright, but incorporates several improvements that lead to substantial performance gains. The two methods are complementary, and can be combined to compute minimal presentations for complexes with millions of generators in a few seconds. The methods have been implemented, and the software is publicly available. We report on experimental evaluations, which demonstrate substantial improvements in performance compared to previously available compression strategies.</p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101940"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49831314","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Subquadratic algorithms for some 3Sum-hard geometric problems in the algebraic decision-tree model","authors":"Boris Aronov ,&nbsp;Mark de Berg ,&nbsp;Jean Cardinal ,&nbsp;Esther Ezra ,&nbsp;John Iacono ,&nbsp;Micha Sharir","doi":"10.1016/j.comgeo.2022.101945","DOIUrl":"https://doi.org/10.1016/j.comgeo.2022.101945","url":null,"abstract":"<div><p>We present subquadratic algorithms in the algebraic decision-tree model for several <span>3Sum</span>-hard geometric problems, all of which can be reduced to the following question: Given two sets <em>A</em>, <em>B</em>, each consisting of <em>n</em> pairwise disjoint segments in the plane, and a set <em>C</em> of <em>n</em> triangles in the plane, we want to count, for each triangle <span><math><mi>Δ</mi><mo>∈</mo><mi>C</mi></math></span>, the number of intersection points between the segments of <em>A</em> and those of <em>B</em> that lie in Δ. We present solutions in the algebraic decision-tree model whose cost is <span><math><mi>O</mi><mo>(</mo><msup><mrow><mi>n</mi></mrow><mrow><mn>60</mn><mo>/</mo><mn>31</mn><mo>+</mo><mi>ε</mi></mrow></msup><mo>)</mo></math></span>, for any <span><math><mi>ε</mi><mo>&gt;</mo><mn>0</mn></math></span>.</p><p>Our approach is based on a primal-dual range searching mechanism, which exploits the multi-level polynomial partitioning machinery recently developed by Agarwal et al. (2021) <span>[3]</span>. A key step in the procedure is a variant of point location in arrangements, say of lines in the plane, which is based solely on the <em>order type</em> of the lines, a “handicap” that turns out to be beneficial for speeding up our algorithm.</p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101945"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49790780","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Many order types on integer grids of polynomial size","authors":"Manfred Scheucher","doi":"10.1016/j.comgeo.2022.101924","DOIUrl":"https://doi.org/10.1016/j.comgeo.2022.101924","url":null,"abstract":"<div><p>Two sets of labeled points <span><math><mo>{</mo><msub><mrow><mi>p</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>,</mo><mo>…</mo><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>n</mi></mrow></msub><mo>}</mo></math></span> and <span><math><mo>{</mo><msub><mrow><mi>q</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>,</mo><mo>…</mo><mo>,</mo><msub><mrow><mi>q</mi></mrow><mrow><mi>n</mi></mrow></msub><mo>}</mo></math></span> are of the same <span><em>labeled </em><em>order type</em></span> if, for every <span><math><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></math></span>, the triples <span><math><mo>(</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>j</mi></mrow></msub><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>k</mi></mrow></msub><mo>)</mo></math></span> and <span><math><mo>(</mo><msub><mrow><mi>q</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>q</mi></mrow><mrow><mi>j</mi></mrow></msub><mo>,</mo><msub><mrow><mi>q</mi></mrow><mrow><mi>k</mi></mrow></msub><mo>)</mo></math></span> have the same orientation. In the 1980's, Goodman, Pollack, and Sturmfels showed that (i) the number of labeled order types on <em>n</em> points is of order <span><math><msup><mrow><mi>n</mi></mrow><mrow><mn>4</mn><mi>n</mi><mo>+</mo><mi>o</mi><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup></math></span>, (ii) all order types can be realized with double-exponential integer coordinates, and that (iii) certain order types indeed require double-exponential integer coordinates. In 2018, Caraballo, Díaz-Báñez, Fabila-Monroy, Hidalgo-Toscano, Leaños, and Montejano showed that at least <span><math><msup><mrow><mi>n</mi></mrow><mrow><mn>3</mn><mi>n</mi><mo>+</mo><mi>o</mi><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup></math></span> labeled <em>n</em>-point order types can be realized on an integer grid of polynomial size. In this article, we improve their result by showing that at least <span><math><msup><mrow><mi>n</mi></mrow><mrow><mn>4</mn><mi>n</mi><mo>+</mo><mi>o</mi><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup></math></span> labeled <em>n</em>-point order types can be realized on an integer grid of polynomial size, which is asymptotically tight in the exponent. Finally we conclude that there are <span><math><msup><mrow><mi>n</mi></mrow><mrow><mn>3</mn><mi>n</mi><mo>+</mo><mi>o</mi><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup></math></span> order types in the <em>unlabeled</em> setting and that <span><math><msup><mrow><mi>n</mi></mrow><mrow><mn>3</mn><mi>n</mi><mo>+</mo><mi>o</mi><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup></math></span> of them can be realized on an integer grid of polynomial size.</p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101924"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49790783","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Computing discrete harmonic differential forms in a given cohomology class using finite element exterior calculus","authors":"Anil N. Hirani ,&nbsp;Kaushik Kalyanaraman ,&nbsp;Han Wang ,&nbsp;Seth Watts","doi":"10.1016/j.comgeo.2022.101937","DOIUrl":"10.1016/j.comgeo.2022.101937","url":null,"abstract":"<div><p><span><span>Computational topology research of the past two decades has emphasized combinatorial techniques while numerical methods such as numerical linear algebra remain underutilized. While the combinatorial techniques have been very successful in diverse areas, for some applications, it is worth considering the numerical counterparts. We discuss one such application. </span>Harmonic forms<span><span> are elements of the kernel of the Hodge Laplacian operator and contain information about the topology of the manifold. If a particular </span>cohomology class<span> is chosen, the closed differential form with the smallest norm in that class is a harmonic form. We use these well-known facts to give an algorithm for solving the following problem: given a piecewise flat manifold simplicial complex (with or without boundary) and a closed </span></span></span>piecewise polynomial<span><span><span> differential form representing a cohomology class, find the discrete harmonic form in that cohomology class. We give a least squares<span> based algorithm to solve this problem and show that the computed form satisfies the finite element exterior calculus (FEEC) equations for being a harmonic form. The piecewise polynomial spaces used are the spaces of trimmed </span></span>polynomial forms, that is, arbitrary </span>degree polynomial generalizations of Whitney forms which are used in FEEC. We also survey other methods for finding harmonic forms.</span></p></div>","PeriodicalId":51001,"journal":{"name":"Computational Geometry-Theory and Applications","volume":"109 ","pages":"Article 101937"},"PeriodicalIF":0.6,"publicationDate":"2023-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43300293","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"计算机科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}