Discrete & Computational Geometry最新文献

The problem of packing of equal disks (or circles) into a rectangle is a fundamental geometric problem. (By a packing here we mean an arrangement of disks in a rectangle without overlapping.) We consider the following algorithmic generalization of the equal disk packing problem. In this problem, for a given packing of equal disks into a rectangle, the question is whether by changing positions of a small number of disks, we can allocate space for packing more disks. More formally, in the repacking problem, for a given set of n equal disks packed into a rectangle and integers k and h, we ask whether it is possible by changing positions of at most h disks to pack (n+k) disks. Thus the problem of packing equal disks is the special case of our problem with (n=h=0). While the computational complexity of packing equal disks into a rectangle remains open, we prove that the repacking problem is NP-hard already for (h=0). Our main algorithmic contribution is an algorithm that solves the repacking problem in time ((h+k)^{mathcal {O}(h+k)}cdot |I|^{mathcal {O}(1)}), where |I| is the input size. That is, the problem is fixed-parameter tractable parameterized by k and h.

Abstract

Let (mathcal {C}) be a set of n axis-aligned cubes of arbitrary sizes in ({mathbb R}^3) in general position. Let (mathcal {U}:=mathcal {U}(mathcal {C})) be their union, and let (kappa ) be the number of vertices on (partial mathcal {U}) ; (kappa ) can vary between O(1) and (Theta (n^2)) . We present a partition of (mathop {textrm{cl}}({mathbb R}^3setminus mathcal {U})) into (O(kappa log ^4 n)) axis-aligned boxes with pairwise-disjoint interiors that can be computed in (O(n log ^2 n + kappa log ^6 n)) time if the faces of (partial mathcal {U}) are pre-computed. We also show that a partition of size (O(sigma log ^4 n + kappa log ^2 n)) , where (sigma ) is the number of input cubes that appear on (partial mathcal {U}) , can be computed in (O(n log ^2 n + sigma log ^8 n + kappa log ^6 n)) time if the faces of (partial mathcal {U}) are pre-computed. The complexity and runtime bounds improve to (O(nlog n)) if all cubes in (mathcal {C}) are congruent and the faces of (partial mathcal {U}) are pre-computed. Finally, we show that if (mathcal {C}) is a set of arbitrary axis-aligned boxes in ({mathbb R}^3) , then a partition of (mathop {textrm{cl}}({mathbb R}^3setminus mathcal {U})) into (O(n^{3/2}+kappa )) boxes can be computed in time (O((n^{3/2}+kappa )log n)) , where (kappa ) is, as above, the number of vertices in (mathcal {U}(mathcal {C})) , which now can vary between O(1) and (Theta (n^3)) .

Inspired by work of Fröberg (1990), and Eagon and Reiner (1998), we define the total k-cut complex of a graph G to be the simplicial complex whose facets are the complements of independent sets of size k in G. We study the homotopy types and combinatorial properties of total cut complexes for various families of graphs, including chordal graphs, cycles, bipartite graphs, the prism (K_n times K_2), and grid graphs, using techniques from algebraic topology and discrete Morse theory.

For (din {mathbb {N}}), a compact sphere packing of Euclidean space ({mathbb {R}}^{d}) is a set of spheres in ({mathbb {R}}^{d}) with disjoint interiors so that the contact hypergraph of the packing is the vertex scheme of a homogeneous simplicial d-complex that covers all of ({mathbb {R}}^{d}). We are motivated by the question: For (d,nin {mathbb {N}}) with (d,nge 2), how many configurations of numbers (0<r_{0}<r_{1}<cdots <r_{n-1}=1) can occur as the radii of spheres in a compact sphere packing of ({mathbb {R}}^{d}) wherein there occur exactly n sizes of sphere? We introduce what we call ‘heteroperturbative sets’ of labeled triangulations of unit spheres and we discuss the existence of non-trivial examples of heteroperturbative sets. For a fixed heteroperturbative set, we discuss how a compact sphere packing may be associated to the heteroperturbative set or not. We proceed to show, for (d,nin {mathbb {N}}) with (d,nge 2) and for a fixed heteroperturbative set, that the collection of all configurations of n distinct positive numbers that can occur as the radii of spheres in a compact packing is finite, when taken over all compact sphere packings of ({mathbb {R}}^{d}) which have exactly n sizes of sphere and which are associated to the fixed heteroperturbative set.

Let (textrm{R}) be a real closed field and (textrm{C}) the algebraic closure of (textrm{R}). We give an algorithm for computing a semi-algebraic basis for the first homology group, (textrm{H}_1(S,{mathbb {F}})), with coefficients in a field ({mathbb {F}}), of any given semi-algebraic set (S subset textrm{R}^k) defined by a closed formula. The complexity of the algorithm is bounded singly exponentially. More precisely, if the given quantifier-free formula involves s polynomials whose degrees are bounded by d, the complexity of the algorithm is bounded by ((s d)^{k^{O(1)}}). This algorithm generalizes well known algorithms having singly exponential complexity for computing a semi-algebraic basis of the zeroth homology group of semi-algebraic sets, which is equivalent to the problem of computing a set of points meeting every semi-algebraically connected component of the given semi-algebraic set at a unique point. It is not known how to compute such a basis for the higher homology groups with singly exponential complexity. As an intermediate step in our algorithm we construct a semi-algebraic subset (Gamma ) of the given semi-algebraic set S, such that (textrm{H}_q(S,Gamma ) = 0) for (q=0,1). We relate this construction to a basic theorem in complex algebraic geometry stating that for any affine variety X of dimension n, there exists Zariski closed subsets

with (dim _textrm{C}Z^{(i)} le i), and (textrm{H}_q(X,Z^{(i)}) = 0) for (0 le q le i). We conjecture a quantitative version of this result in the semi-algebraic category, with X and (Z^{(i)}) replaced by closed semi-algebraic sets. We make initial progress on this conjecture by proving the existence of (Z^{(0)}) and (Z^{(1)}) with complexity bounded singly exponentially (previously, such an algorithm was known only for constructing (Z^{(0)})).

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 (f:X rightarrow Y) between closed and bounded semi-algebraic sets. For every fixed (ell ge 0) we give an algorithm with singly exponential complexity that computes bases of the homology groups (text{ H}_i(X), text{ H}_i(Y)) (with rational coefficients) and a matrix with respect to these bases of the induced linear maps (text{ H}_i(f):text{ H}_i(X) rightarrow text{ H}_i(Y), 0 le i le ell ). 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.

This paper considers an extremal version of the Erdős distinct distances problem. For a point set (P subset {mathbb {R}}^d), let (Delta (P)) denote the set of all Euclidean distances determined by P. Our main result is the following: if (Delta (A^d) ll |A|^2) and (d ge 5), then there exists (A' subset A) with (|A'| ge |A|/2) such that (|A'-A'| ll |A| log |A|). This is one part of a more general result, which says that, if the growth of (|Delta (A^d)|) is restricted, it must be the case that A has some additive structure. More specifically, for any two integers k, n, we have the following information: if

then there exists (A' subset A) with (|A'| ge |A|/2) and

These results are higher dimensional analogues of a result of Hanson [4], who considered the two-dimensional case.

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.

Given a set of n 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 n! ways, which is the maximum number of ways for n points. But what is the minimum number? It is shown that this number is (roughly) at least (3^n) and at most (12.29^n).

The convexity of a set can be generalized to the two weaker notions of positive reach and r-convexity; both describe the regularity of a set’s boundary. For any compact subset of ({{mathbb {R}}}^d), 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 (beta )-reach, a generalization of the reach that excludes small-scale features of size less than a parameter (beta in [0,infty )). Numerical studies suggest how the (beta )-reach can be used in high-dimension to infer the reach and other geometric properties of smooth submanifolds.