{"title":"Computing the Homology Functor on Semi-algebraic Maps and Diagrams","authors":"Saugata Basu, Negin Karisani","doi":"10.1007/s00454-024-00627-z","DOIUrl":"https://doi.org/10.1007/s00454-024-00627-z","url":null,"abstract":"<p>Developing an algorithm for computing the Betti numbers of semi-algebraic sets with singly exponential complexity has been a holy grail in algorithmic semi-algebraic geometry and only partial results are known. In this paper we consider the more general problem of computing the image under the homology functor of a continuous semi-algebraic map <span>(f:X rightarrow Y)</span> between closed and bounded semi-algebraic sets. For every fixed <span>(ell ge 0)</span> we give an algorithm with singly exponential complexity that computes bases of the homology groups <span>(text{ H}_i(X), text{ H}_i(Y))</span> (with rational coefficients) and a matrix with respect to these bases of the induced linear maps <span>(text{ H}_i(f):text{ H}_i(X) rightarrow text{ H}_i(Y), 0 le i le ell )</span>. We generalize this algorithm to more general (zigzag) diagrams of continuous semi-algebraic maps between closed and bounded semi-algebraic sets and give a singly exponential algorithm for computing the homology functors on such diagrams. This allows us to give an algorithm with singly exponential complexity for computing barcodes of semi-algebraic zigzag persistent homology in small dimensions.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"13 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139766220","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Convexity, Elementary Methods, and Distances","authors":"Oliver Roche-Newton, Dmitrii Zhelezov","doi":"10.1007/s00454-023-00625-7","DOIUrl":"https://doi.org/10.1007/s00454-023-00625-7","url":null,"abstract":"<p>This paper considers an extremal version of the Erdős distinct distances problem. For a point set <span>(P subset {mathbb {R}}^d)</span>, let <span>(Delta (P))</span> denote the set of all Euclidean distances determined by <i>P</i>. Our main result is the following: if <span>(Delta (A^d) ll |A|^2)</span> and <span>(d ge 5)</span>, then there exists <span>(A' subset A)</span> with <span>(|A'| ge |A|/2)</span> such that <span>(|A'-A'| ll |A| log |A|)</span>. This is one part of a more general result, which says that, if the growth of <span>(|Delta (A^d)|)</span> is restricted, it must be the case that <i>A</i> has some additive structure. More specifically, for any two integers <i>k</i>, <i>n</i>, we have the following information: if </p><span>$$begin{aligned} | Delta (A^{2k+3})| le |A|^n end{aligned}$$</span><p>then there exists <span>(A' subset A)</span> with <span>(|A'| ge |A|/2)</span> and </p><span>$$begin{aligned} | kA'- kA'| le k^2|A|^{2n-3}log |A|. end{aligned}$$</span><p>These results are higher dimensional analogues of a result of Hanson [4], who considered the two-dimensional case.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"209 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139680311","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Locally Finite Completions of Polyhedral Complexes","authors":"Desmond Coles, Netanel Friedenberg","doi":"10.1007/s00454-024-00629-x","DOIUrl":"https://doi.org/10.1007/s00454-024-00629-x","url":null,"abstract":"<p>We develop a method for subdividing polyhedral complexes in a way that restricts the possible recession cones and allows one to work with a fixed class of polyhedron. We use these results to construct locally finite completions of rational polyhedral complexes whose recession cones lie in a fixed fan, locally finite polytopal completions of polytopal complexes, and locally finite zonotopal completions of zonotopal complexes.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"10 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139677471","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Peeling Sequences","authors":"Adrian Dumitrescu, Géza Tóth","doi":"10.1007/s00454-023-00616-8","DOIUrl":"https://doi.org/10.1007/s00454-023-00616-8","url":null,"abstract":"<p>Given a set of <i>n</i> labeled points in general position in the plane, we remove all of its points one by one. At each step, one point from the convex hull of the remaining set is erased. In how many ways can the process be carried out? The answer obviously depends on the point set. If the points are in convex position, there are exactly <i>n</i>! ways, which is the maximum number of ways for <i>n</i> points. But what is the minimum number? It is shown that this number is (roughly) at least <span>(3^n)</span> and at most <span>(12.29^n)</span>.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"24 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-02-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139664084","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Computable Bounds for the Reach and r-Convexity of Subsets of $${{mathbb {R}}}^d$$","authors":"Ryan Cotsakis","doi":"10.1007/s00454-023-00624-8","DOIUrl":"https://doi.org/10.1007/s00454-023-00624-8","url":null,"abstract":"<p>The convexity of a set can be generalized to the two weaker notions of positive reach and <i>r</i>-convexity; both describe the regularity of a set’s boundary. For any compact subset of <span>({{mathbb {R}}}^d)</span>, we provide methods for computing upper bounds on these quantities from point cloud data. The bounds converge to the respective quantities as the sampling scale of the point cloud decreases, and the rate of convergence for the bound on the reach is given under a weak regularity condition. We also introduce the <span>(beta )</span>-reach, a generalization of the reach that excludes small-scale features of size less than a parameter <span>(beta in [0,infty ))</span>. Numerical studies suggest how the <span>(beta )</span>-reach can be used in high-dimension to infer the reach and other geometric properties of smooth submanifolds.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"5 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139589839","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Cone of $$5times 5$$ Completely Positive Matrices","authors":"Max Pfeffer, José Alejandro Samper","doi":"10.1007/s00454-023-00620-y","DOIUrl":"https://doi.org/10.1007/s00454-023-00620-y","url":null,"abstract":"<p>We study the cone of completely positive (cp) matrices for the first interesting case <span>(n = 5)</span>. This is a semialgebraic set for which the polynomial equalities and inequlities that define its boundary can be derived. We characterize the different loci of this boundary and we examine the two open sets with cp-rank 5 or 6. A numerical algorithm is presented that is fast and able to compute the cp-factorization even for matrices in the boundary. With our results, many new example cases can be produced and several insightful numerical experiments are performed that illustrate the difficulty of the cp-factorization problem.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"5 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139560497","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Almost Congruent Triangles","authors":"","doi":"10.1007/s00454-023-00623-9","DOIUrl":"https://doi.org/10.1007/s00454-023-00623-9","url":null,"abstract":"<h3>Abstract</h3> <p>Almost 50 years ago Erdős and Purdy asked the following question: Given <em>n</em> points in the plane, how many triangles can be approximate congruent to equilateral triangles? They pointed out that by dividing the points evenly into three small clusters built around the three vertices of a fixed equilateral triangle, one gets at least <span> <span>(leftlfloor frac{n}{3} rightrfloor cdot leftlfloor frac{n+1}{3} rightrfloor cdot leftlfloor frac{n+2}{3} rightrfloor )</span> </span> such approximate copies. In this paper we provide a matching upper bound and thereby answer their question. More generally, for every triangle <em>T</em> we determine the maximum number of approximate congruent triangles to <em>T</em> in a point set of size <em>n</em>. Parts of our proof are based on hypergraph Turán theory: for each point set in the plane and a triangle <em>T</em>, we construct a 3-uniform hypergraph <span> <span>(mathcal {H}=mathcal {H}(T))</span> </span>, which contains no hypergraph as a subgraph from a family of forbidden hypergraphs <span> <span>(mathcal {F}=mathcal {F}(T))</span> </span>. Our upper bound on the number of edges of <span> <span>(mathcal {H})</span> </span> will determine the maximum number of triangles that are approximate congruent to <em>T</em>.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"75 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-01-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139460336","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Henk Alkema, Mark de Berg, Remco van der Hofstad, Sándor Kisfaludi-Bak
{"title":"Euclidean TSP in Narrow Strips","authors":"Henk Alkema, Mark de Berg, Remco van der Hofstad, Sándor Kisfaludi-Bak","doi":"10.1007/s00454-023-00609-7","DOIUrl":"https://doi.org/10.1007/s00454-023-00609-7","url":null,"abstract":"<p>We investigate how the complexity of <span>Euclidean TSP</span> for point sets <i>P</i> inside the strip <span>((-infty ,+infty )times [0,delta ])</span> depends on the strip width <span>(delta )</span>. We obtain two main results.</p><ul>\u0000<li>\u0000<p>For the case where the points have distinct integer <i>x</i>-coordinates, we prove that a shortest bitonic tour (which can be computed in <span>(O(nlog ^2 n))</span> time using an existing algorithm) is guaranteed to be a shortest tour overall when <span>(delta leqslant 2sqrt{2})</span>, a bound which is best possible.</p>\u0000</li>\u0000<li>\u0000<p>We present an algorithm that is fixed-parameter tractable with respect to <span>(delta )</span>. Our algorithm has running time <span>(2^{O(sqrt{delta })} n + O(delta ^2 n^2))</span> for sparse point sets, where each <span>(1times delta )</span> rectangle inside the strip contains <i>O</i>(1) points. For random point sets, where the points are chosen uniformly at random from the rectangle <span>([0,n]times [0,delta ])</span>, it has an expected running time of <span>(2^{O(sqrt{delta })} n)</span>. These results generalise to point sets <i>P</i> inside a hypercylinder of width <span>(delta )</span>. In this case, the factors <span>(2^{O(sqrt{delta })})</span> become <span>(2^{O(delta ^{1-1/d})})</span>.</p>\u0000</li>\u0000</ul>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"72 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139408141","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Counting Arcs in $${mathbb {F}}_q^2$$","authors":"Krishnendu Bhowmick, Oliver Roche-Newton","doi":"10.1007/s00454-023-00622-w","DOIUrl":"https://doi.org/10.1007/s00454-023-00622-w","url":null,"abstract":"<p>An arc in <span>(mathbb F_q^2)</span> is a set <span>(P subset mathbb F_q^2)</span> such that no three points of <i>P</i> are collinear. We use the method of hypergraph containers to prove several counting results for arcs. Let <span>({mathcal {A}}(q))</span> denote the family of all arcs in <span>(mathbb F_q^2)</span>. Our main result is the bound </p><span>$$begin{aligned} |{mathcal {A}}(q)| le 2^{(1+o(1))q}. end{aligned}$$</span><p>This matches, up to the factor hidden in the <i>o</i>(1) notation, the trivial lower bound that comes from considering all subsets of an arc of size <i>q</i>. We also give upper bounds for the number of arcs of a fixed (large) size. Let <span>(k ge q^{2/3}(log q)^3)</span>, and let <span>({mathcal {A}}(q,k))</span> denote the family of all arcs in <span>(mathbb F_q^2)</span> with cardinality <i>k</i>. We prove that </p><span>$$begin{aligned} |{mathcal {A}}(q,k)| le left( {begin{array}{c}(1+o(1))q kend{array}}right) . end{aligned}$$</span><p>This result improves a bound of Roche-Newton and Warren [12]. A nearly matching lower bound </p><span>$$begin{aligned} |{mathcal {A}}(q,k)| ge left( {begin{array}{c}q kend{array}}right) end{aligned}$$</span><p>follows by considering all subsets of size <i>k</i> of an arc of size <i>q</i>.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"27 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-01-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139415409","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Topological Optimization with Big Steps","authors":"Arnur Nigmetov, Dmitriy Morozov","doi":"10.1007/s00454-023-00613-x","DOIUrl":"https://doi.org/10.1007/s00454-023-00613-x","url":null,"abstract":"<p>Using persistent homology to guide optimization has emerged as a novel application of topological data analysis. Existing methods treat persistence calculation as a black box and backpropagate gradients only onto the simplices involved in particular pairs. We show how the cycles and chains used in the persistence calculation can be used to prescribe gradients to larger subsets of the domain. In particular, we show that in a special case, which serves as a building block for general losses, the problem can be solved exactly in linear time. This relies on another contribution of this paper, which eliminates the need to examine a factorial number of permutations of simplices with the same value. We present empirical experiments that show the practical benefits of our algorithm: the number of steps required for the optimization is reduced by an order of magnitude.</p>","PeriodicalId":50574,"journal":{"name":"Discrete & Computational Geometry","volume":"3 1","pages":""},"PeriodicalIF":0.8,"publicationDate":"2024-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139372938","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}