Point Location in Constant Time

Sairam Chaganti,Yijie Han

2023-12-22

Abstract:We preprocess the input subdivision with $n$ points on the plane in $O(n\sqrt{\log n})$ time to facilitate point location in constant time. Previously the preprocessing time is $O(n\log n)$ and point location takes $O(\log n)$ time.

Computational Geometry

What problem does this paper attempt to address?

Point Location in Dynamic Planar Subdivisions

Eunjin Oh,Hee-Kap Ahn

DOI: https://doi.org/10.48550/arXiv.1803.04325

2018-03-12

Abstract:We study the point location problem on dynamic planar subdivisions that allows insertions and deletions of edges. In our problem, the underlying graph of a subdivision is not necessarily connected. We present a data structure of linear size for such a dynamic planar subdivision that supports sublinear-time update and polylogarithmic-time query. Precisely, the amortized update time is $O(\sqrt{n}\log n(\log\log n)^{3/2})$ and the query time is $O(\log n(\log\log n)^2)$, where $n$ is the number of edges in the subdivision. This answers a question posed by Snoeyink in the Handbook of Computational Geometry. When only deletions of edges are allowed, the update time and query time are just $O(\alpha(n))$ and $O(\log n)$, respectively.

Computational Geometry
Improved Implementation of Point Location in General Two-Dimensional Subdivisions

Michael Hemmer,Michal Kleinbort,Dan Halperin

DOI: https://doi.org/10.1007/978-3-642-33090-2_53

2012-09-23

Abstract:We present a major revamp of the point-location data structure for general two-dimensional subdivisions via randomized incremental construction, implemented in CGAL, the Computational Geometry Algorithms Library. We can now guarantee that the constructed directed acyclic graph G is of linear size and provides logarithmic query time. Via the construction of the Voronoi diagram for a given point set S of size n, this also enables nearest-neighbor queries in guaranteed O(log n) time. Another major innovation is the support of general unbounded subdivisions as well as subdivisions of two-dimensional parametric surfaces such as spheres, tori, cylinders. The implementation is exact, complete, and general, i.e., it can also handle non-linear subdivisions. Like the previous version, the data structure supports modifications of the subdivision, such as insertions and deletions of edges, after the initial preprocessing. A major challenge is to retain the expected O(n log n) preprocessing time while providing the above (deterministic) space and query-time guarantees. We describe an efficient preprocessing algorithm, which explicitly verifies the length L of the longest query path in O(n log n) time. However, instead of using L, our implementation is based on the depth D of G. Although we prove that the worst case ratio of D and L is Theta(n/log n), we conjecture, based on our experimental results, that this solution achieves expected O(n log n) preprocessing time.

Computational Geometry,Data Structures and Algorithms
Computing a Rectilinear Steiner Minimal Tree in $$n^{o(\sqrt N )}$$ Time

Clark Thomborson,Linda L. Deneen,Gary Shute

DOI: https://doi.org/10.1007/3-540-18099-0_44

1987-01-01

Abstract:We propose an algorithm for computing a rectilinear Steiner minimal tree on n points in $2^{O(\sqrt n \log n)}$time and O(n2) space. This is an asymptotic improvement on the 2O(n) time required by current algorithms. If the points are distributed uniformly at random on the unit square, the Steiner tree calculated by our algorithm is minimal with high probability. The “constant factors” of our algorithm are such that it should be feasible to obtain exact solutions for n-point problems whenever n≤25. Previously, only problems of size n≤20 were feasible.
Point-in-Convex Polygon and Point-in-Convex Polyhedron Algorithms with O(1) Complexity using Space Subdivision

Vaclav Skala

DOI: https://doi.org/10.48550/arXiv.2208.03682

2022-08-07

Graphics

Abstract:There are many space subdivision and space partitioning techniques used in many algorithms to speed up computations. They mostly rely on orthogonal space subdivision, resp. using hierarchical data structures, e.g. BSP trees, quadtrees, octrees, kd-trees, bounding volume hierarchies, etc. However in some applications a non-orthogonal space subdivision can offer new ways for actual speed up. In the case of convex polygon in E3 a simple Point-in-Polygon test is of the O(N) complexity and the optimal algorithm is of O(lg N) computational complexity. In the E3 case, the complexity is O(N) even for the convex polyhedron as no ordering is defined. New Point-in-Convex Polygon and Point-in-Convex Polyhedron algorithms are presented based on space subdivision in the preprocessing stage resulting to O(1) run-time complexity. The presented approach is simple to implement. Due to the principle of duality, dual problems, e.g. line-convex polygon, line clipping, can be solved in a similarly.
A Constrained Minimum Spanning Tree Problem

GT Chen,GC Zhang

DOI: https://doi.org/10.1016/s0305-0548(99)00061-1

IF: 5.159

2000-01-01

Computers & Operations Research

Abstract:In this paper, we give an O(n(2)) algorithm for finding a point on a given line I which, if added to a given set of n points located on one side of l, will result in the minimum spanning tree.
Storage in Computational Geometry

Yijie Han,Sanjeev Saxena

DOI: https://doi.org/10.48550/arXiv.2302.11821

2023-02-23

Abstract:We show that $n$ real numbers can be stored in a constant number of real numbers such that each original real number can be fetched in $O(\log n)$ time. Although our result has implications for many computational geometry problems, we show here, combined with Han's $O(n\sqrt{\log n})$ time real number sorting algorithm [3, <a class="link-https" data-arxiv-id="1801.00776" href="https://arxiv.org/abs/1801.00776">arXiv:1801.00776</a>], we can improve the complexity of Kirkpatrick's point location algorithm [8] to $O(n\sqrt{\log n})$ preprocessing time, a constant number of real numbers for storage and $O(\log n)$ point location time. Kirkpatrick's algorithm uses $O(n\log n)$ preprocessing time, $O(n)$ storage and $O(\log n)$ point location time. The complexity results in Kirkpatrick's algorithm was the previous best result. Although Lipton and Tarjan's algorithm [10] predates Kirkpatrick's algorithm and has the same complexity, Kirkpatrick's algorithm is simpler and has a better structure. This paper can be viewed as a companion paper of paper [3, <a class="link-https" data-arxiv-id="1801.00776" href="https://arxiv.org/abs/1801.00776">arXiv:1801.00776</a>].

Computational Geometry,Data Structures and Algorithms
Planar $β$-skeletons via point location in monotone subdivisions of subset of lunes

Mirosław Kowaluk

DOI: https://doi.org/10.48550/arXiv.1411.5459

IF: 0.455

2014-11-20

Computational Geometry

Abstract:We present a new algorithm for lune-based $\beta$-skeletons for sets of $n$ points in the plane, for $\beta \in (2,\infty]$, the only case when optimal algorithms are not known. The running time of the algorithm is $O(n^{3/2} \log^{1/2} n)$, which is the best known and is an improvement of Rao and Mukhopadhyay \cite{rm97} result. The method is based on point location in monotonic subdivisions of arrangements of curve segments.
Half-Plane Voronoi Diagram

Chenglin Fan,Jun Luo,Jinfei Liu,Yinfeng Xu

DOI: https://doi.org/10.1109/isvd.2011.25

2011-01-01

Abstract:In normal Voronoi diagram, each site is able to see all points in the plane. In this paper, we study the problem such that each site is only able to see half-plane and construct the so-called Half-plane Voronoi Diagram (HPVD). We show that the half-plane Voronoi cell of each site is not necessary convex and it could consist of many disjoint regions. We prove that the complexity of the HPVD of n sites is $O(n^2)$. Then we give an algorithm of $O(n\log n)$ time and $O(n)$ space to construct HPVD such that the boundary lines of all half-planes are parallel and the visible half-planes are in the same side of the corresponding boundary lines. For the parallel boundary lines where the visible half-planes could be in either side of the corresponding boundary lines, we give an algorithm of $O(n^2)$ time and $O(n^2)$ space to construct HPVD.
Data Structures for Halfplane Proximity Queries and Incremental Voronoi Diagrams

Boris Aronov,Prosenjit Bose,Erik D. Demaine,Joachim Gudmundsson,John Iacono,Stefan Langerman,Michiel Smid

DOI: https://doi.org/10.48550/arXiv.cs/0512091

2017-10-14

Abstract:We consider preprocessing a set $S$ of $n$ points in convex position in the plane into a data structure supporting queries of the following form: given a point $q$ and a directed line $\ell$ in the plane, report the point of $S$ that is farthest from (or, alternatively, nearest to) the point $q$ among all points to the left of line $\ell$. We present two data structures for this problem. The first data structure uses $O(n^{1+\varepsilon})$ space and preprocessing time, and answers queries in $O(2^{1/\varepsilon} \log n)$ time, for any $0 < \varepsilon < 1$. The second data structure uses $O(n \log^3 n)$ space and polynomial preprocessing time, and answers queries in $O(\log n)$ time. These are the first solutions to the problem with $O(\log n)$ query time and $o(n^2)$ space. The second data structure uses a new representation of nearest- and farthest-point Voronoi diagrams of points in convex position. This representation supports the insertion of new points in clockwise order using only $O(\log n)$ amortized pointer changes, in addition to $O(\log n)$-time point-location queries, even though every such update may make $\Theta(n)$ combinatorial changes to the Voronoi diagram. This data structure is the first demonstration that deterministically and incrementally constructed Voronoi diagrams can be maintained in $o(n)$ amortized pointer changes per operation while keeping $O(\log n)$-time point-location queries.

Computational Geometry,Data Structures and Algorithms
Dynamic Planar Point Location in External Memory

J. Ian Munro,Yakov Nekrich

DOI: https://doi.org/10.48550/arXiv.1903.06601

2019-03-15

Abstract:In this paper we describe a fully-dynamic data structure for the planar point location problem in the external memory model. Our data structure supports queries in $O(\log_B n(\log\log_B n)^3))$ I/Os and updates in $O(\log_B n(\log\log_B n)^2))$ amortized I/Os, where $n$ is the number of segments in the subdivision and $B$ is the block size. This is the first dynamic data structure with almost-optimal query cost. For comparison all previously known results for this problem require $O(\log_B^2 n)$ I/Os to answer queries. Our result almost matches the best known upper bound in the internal-memory model.

Data Structures and Algorithms
Orthogonal Point Location and Rectangle Stabbing Queries in 3-d

Timothy M. Chan,Yakov Nekrich,Saladi Rahul,Konstantinos Tsakalidis

DOI: https://doi.org/10.48550/arXiv.1805.08602

2018-05-19

Abstract:In this work, we present a collection of new results on two fundamental problems in geometric data structures: orthogonal point location and rectangle stabbing. -We give the first linear-space data structure that supports 3-d point location queries on $n$ disjoint axis-aligned boxes with optimal $O\left( \log n\right)$ query time in the (arithmetic) pointer machine model. This improves the previous $O\left( \log^{3/2} n \right)$ bound of Rahul [SODA 2015]. We similarly obtain the first linear-space data structure in the I/O model with optimal query cost, and also the first linear-space data structure in the word RAM model with sub-logarithmic query time. -We give the first linear-space data structure that supports 3-d $4$-sided and $5$-sided rectangle stabbing queries in optimal $O(\log_wn+k)$ time in the word RAM model. We similarly obtain the first optimal data structure for the closely related problem of 2-d top-$k$ rectangle stabbing in the word RAM model, and also improved results for 3-d 6-sided rectangle stabbing. For point location, our solution is simpler than previous methods, and is based on an interesting variant of the van Emde Boas recursion, applied in a round-robin fashion over the dimensions, combined with bit-packing techniques. For rectangle stabbing, our solution is a variant of Alstrup, Brodal, and Rauhe's grid-based recursive technique (FOCS 2000), combined with a number of new ideas.

Computational Geometry,Data Structures and Algorithms
A simple and efficient preprocessing step for convex hull problem

Mohammad Heydari,Ashkan Khalifeh

DOI: https://doi.org/10.1142/S179383092350091X

2023-04-09

Abstract:The present paper is concerned with a recursive algorithm as a preprocessing step to find the convex hull of $n$ random points uniformly distributed in the plane. For such a set of points, it is shown that eliminating all but $O(\log n)$ of points can derive the same convex hull as the input set. Finally it will be shown that the running time of the algorithm is $O(n)

Data Structures and Algorithms
Algorithms for Halfplane Coverage and Related Problems

Haitao Wang,Jie Xue

2024-02-26

Abstract:Given in the plane a set of points and a set of halfplanes, we consider the problem of computing a smallest subset of halfplanes whose union covers all points. In this paper, we present an $O(n^{4/3}\log^{5/3}n\log^{O(1)}\log n)$-time algorithm for the problem, where $n$ is the total number of all points and halfplanes. This improves the previously best algorithm of $n^{10/3}2^{O(\log^*n)}$ time by roughly a quadratic factor. For the special case where all halfplanes are lower ones, our algorithm runs in $O(n\log n)$ time, which improves the previously best algorithm of $n^{4/3}2^{O(\log^*n)}$ time and matches an $\Omega(n\log n)$ lower bound. Further, our techniques can be extended to solve a star-shaped polygon coverage problem in $O(n\log n)$ time, which in turn leads to an $O(n\log n)$-time algorithm for computing an instance-optimal $\epsilon$-kernel of a set of $n$ points in the plane. Agarwal and Har-Peled presented an $O(nk\log n)$-time algorithm for this problem in SoCG 2023, where $k$ is the size of the $\epsilon$-kernel; they also raised an open question whether the problem can be solved in $O(n\log n)$ time. Our result thus answers the open question affirmatively.

Data Structures and Algorithms,Computational Geometry
Polar, Spherical and Orthogonal Space Subdivisions for an Algorithm Acceleration: O(1) Point-in-Polygon/Polyhedron Test

Vaclav Skala

DOI: https://doi.org/10.48550/arXiv.2208.12488

IF: 0.455

2022-08-26

Computational Geometry

Abstract:Acceleration of algorithms is becoming a crucial problem, if larger data sets are to be processed. Evaluation of algorithms is mostly done by using computational geometry approach and evaluation of computational complexity. However in todays engineering problems this approach does not respect that number of processed items is always limited and a significant role plays also speed of read/write operations. One general method how to speed up an algorithm is application of space subdivision technique and usually the orthogonal space subdivision is used. In this paper non-orthogonal subdivisions are described. The proposed approach can significantly improve memory consumption and run-time complexity. The proposed modified space subdivision techniques are demonstrated on two simple problems Point-in-Convex Polygon and Point-in-Convex Polyhedron tests.
A 3-Space Partition and Its Applications

F. Frances Yao

DOI: https://doi.org/10.1145/800061.808755

1983-01-01

Abstract:Let S be a set of n points in three-dimensional space. It is shown that one can always find three planes that divide S into eight open regions, of which no seven together contain more than α n points where α is a constant < 1. This result gives rise to a data structure, what we call an octant-tree, for representing any point set in 3-space. Efficient solutions to various data retrieval problems are readily available with this structure. For example, using octant-trees, one can answer in sublinear time T ( n ) @@@@ O ( n 0.98 ) 1) half-space queries: find all points of S that lie to one side of a plane P ; 2) polytope queries: find all points that lie inside (outside) a polytope; and 3) circular queries in E 2 : given a planar set S , find all points that lie within (without) a circle of radius r and center c for any r and c . An octant-tree for n points occupies O ( n ) space and can be constructed with O ( n 4 ) preprocessing time.
IMPROVED ALGORITHMS FOR THE POINT-SET EMBEDDABILITY PROBLEM FOR PLANE 3-TREES

TANAEEM M. MOOSA,M. SOHEL RAHMAN

DOI: https://doi.org/10.1142/s1793830912500097

2012-03-01

Discrete Mathematics, Algorithms and Applications

Abstract:In the point-set embeddability problem, we are given a plane graph G with n vertices and a point set S with the same number of points. Now the goal is to answer the question whether there exists a straight-line drawing of G such that each vertex is represented as a distinct point of S as well as to provide an embedding if one does exist. This problem has recently been solved in O(n 2 log n) time for plane 3-trees. In this paper, we present a new efficient algorithm with time complexity O(n 4/3+ε log n). We also present an O(nk 4 ) time algorithm for the case when |S| = k > n. This is a significant improvement over the best algorithm for this case in the literature, which runs in O(nk 8 ) time.
Optimal Point Movement for Covering Circular Regions

Danny Z. Chen,Xuehou Tan,Haitao Wang,Gangshan Wu

DOI: https://doi.org/10.48550/arXiv.1107.1012

2011-07-06

Abstract:Given $n$ points in a circular region $C$ in the plane, we study the problems of moving the $n$ points to its boundary to form a regular $n$-gon such that the maximum (min-max) or the sum (min-sum) of the Euclidean distances traveled by the points is minimized. The problems have applications, e.g., in mobile sensor barrier coverage of wireless sensor networks. The min-max problem further has two versions: the decision version and optimization version. For the min-max problem, we present an $O(n\log^2 n)$ time algorithm for the decision version and an $O(n\log^3 n)$ time algorithm for the optimization version. The previously best algorithms for the two problem versions take $O(n^{3.5})$ time and $O(n^{3.5}\log n)$ time, respectively. For the min-sum problem, we show that a special case with all points initially lying on the boundary of the circular region can be solved in $O(n^2)$ time, improving a previous $O(n^4)$ time solution. For the general min-sum problem, we present a 3-approximation $O(n^2)$ time algorithm, improving the previous $(1+\pi)$-approximation $O(n^2)$ time algorithm. A by-product of our techniques is an algorithm for dynamically maintaining the maximum matching of a circular convex bipartite graph; our algorithm can handle each vertex insertion or deletion on the graph in $O(\log^2 n)$ time. This result is interesting in its own right.

Computational Geometry,Data Structures and Algorithms
Eight-Partitioning Points in 3D, and Efficiently Too

Boris Aronov,Abdul Basit,Indu Ramesh,Gianluca Tasinato,Uli Wagner

2024-07-19

Abstract:An eight-partition of a finite set of points (respectively, of a continuous mass distribution) in $\mathbb{R}^3$ consists of three planes that divide the space into $8$ octants, such that each open octant contains at most $1/8$ of the points (respectively, of the mass). In 1966, Hadwiger showed that any mass distribution in $\mathbb{R}^3$ admits an eight-partition; moreover, one can prescribe the normal direction of one of the three planes. The analogous result for finite point sets follows by a standard limit argument. We prove the following variant of this result: Any mass distribution (or point set) in $\mathbb{R}^3$ admits an eight-partition for which the intersection of two of the planes is a line with a prescribed direction. Moreover, we present an efficient algorithm for calculating an eight-partition of a set of $n$ points in $\mathbb{R}^3$ (with prescribed normal direction of one of the planes) in time $O^{*}(n^{5/2})$.

Computational Geometry,Combinatorics
Fixed-Parameter Algorithms for Rectilinear Steiner tree and Rectilinear Traveling Salesman Problem in the plane

Hadrien Cambazard,Nicolas Catusse

DOI: https://doi.org/10.1016/j.ejor.2018.03.042

2019-07-23

Abstract:Given a set $P$ of $n$ points with their pairwise distances, the traveling salesman problem (TSP) asks for a shortest tour that visits each point exactly once. A TSP instance is rectilinear when the points lie in the plane and the distance considered between two points is the $l_1$ distance. In this paper, a fixed-parameter algorithm for the Rectilinear TSP is presented and relies on techniques for solving TSP on bounded-treewidth graphs. It proves that the problem can be solved in $O\left(nh7^h\right)$ where $h \leq n$ denotes the number of horizontal lines containing the points of $P$. The same technique can be directly applied to the problem of finding a shortest rectilinear Steiner tree that interconnects the points of $P$ providing a $O\left(nh5^h\right)$ time complexity. Both bounds improve over the best time bounds known for these problems.

Data Structures and Algorithms,Computational Geometry,Discrete Mathematics
A New Algorithm of Planar Point Set Based on Sorted Points

CHEN Ping,WANG Guo-zhao

DOI: https://doi.org/10.3969/j.issn.1001-7119.2007.06.019

2007-01-01

Abstract:This paper improves the fast convex hull algorithm of planar point set based on sorted simple polygon.The new algorithm can avoid the problem of superposition of extreme points.The algorithm not only reaches the theoretical lower bound of O(nlnn),but also is very simple and easy to be realized.

Point Location in Constant Time

Point Location in Dynamic Planar Subdivisions

Improved Implementation of Point Location in General Two-Dimensional Subdivisions

Computing a Rectilinear Steiner Minimal Tree in $$n^{o(\sqrt N )}$$ Time

Point-in-Convex Polygon and Point-in-Convex Polyhedron Algorithms with O(1) Complexity using Space Subdivision

A Constrained Minimum Spanning Tree Problem

Storage in Computational Geometry

Planar $β$-skeletons via point location in monotone subdivisions of subset of lunes

Half-Plane Voronoi Diagram

Data Structures for Halfplane Proximity Queries and Incremental Voronoi Diagrams

Dynamic Planar Point Location in External Memory

Orthogonal Point Location and Rectangle Stabbing Queries in 3-d

A simple and efficient preprocessing step for convex hull problem

Algorithms for Halfplane Coverage and Related Problems

Polar, Spherical and Orthogonal Space Subdivisions for an Algorithm Acceleration: O(1) Point-in-Polygon/Polyhedron Test

A 3-Space Partition and Its Applications

IMPROVED ALGORITHMS FOR THE POINT-SET EMBEDDABILITY PROBLEM FOR PLANE 3-TREES

Optimal Point Movement for Covering Circular Regions

Eight-Partitioning Points in 3D, and Efficiently Too

Fixed-Parameter Algorithms for Rectilinear Steiner tree and Rectilinear Traveling Salesman Problem in the plane

A New Algorithm of Planar Point Set Based on Sorted Points