101,380 research outputs found
Dominance Product and High-Dimensional Closest Pair under
Given a set of points in , the Closest Pair problem is
to find a pair of distinct points in at minimum distance. When is
constant, there are efficient algorithms that solve this problem, and fast
approximate solutions for general . However, obtaining an exact solution in
very high dimensions seems to be much less understood. We consider the
high-dimensional Closest Pair problem, where for some , and the underlying metric is .
We improve and simplify previous results for Closest Pair, showing
that it can be solved by a deterministic strongly-polynomial algorithm that
runs in time, and by a randomized algorithm that runs in
expected time, where is the time bound for computing the
{\em dominance product} for points in . That is a matrix ,
such that ; this is the
number of coordinates at which dominates . For integer coordinates
from some interval , we obtain an algorithm that runs in
time, where
is the exponent of multiplying an matrix by an
matrix.
We also give slightly better bounds for , by using more recent
rectangular matrix multiplication bounds. Computing the dominance product
itself is an important task, since it is applied in many algorithms as a major
black-box ingredient, such as algorithms for APBP (all pairs bottleneck paths),
and variants of APSP (all pairs shortest paths)
Net and Prune: A Linear Time Algorithm for Euclidean Distance Problems
We provide a general framework for getting expected linear time constant
factor approximations (and in many cases FPTAS's) to several well known
problems in Computational Geometry, such as -center clustering and farthest
nearest neighbor. The new approach is robust to variations in the input
problem, and yet it is simple, elegant and practical. In particular, many of
these well studied problems which fit easily into our framework, either
previously had no linear time approximation algorithm, or required rather
involved algorithms and analysis. A short list of the problems we consider
include farthest nearest neighbor, -center clustering, smallest disk
enclosing points, th largest distance, th smallest -nearest
neighbor distance, th heaviest edge in the MST and other spanning forest
type problems, problems involving upward closed set systems, and more. Finally,
we show how to extend our framework such that the linear running time bound
holds with high probability
Preconditioned Locally Harmonic Residual Method for Computing Interior Eigenpairs of Certain Classes of Hermitian Matrices
We propose a Preconditioned Locally Harmonic Residual (PLHR) method for
computing several interior eigenpairs of a generalized Hermitian eigenvalue
problem, without traditional spectral transformations, matrix factorizations,
or inversions. PLHR is based on a short-term recurrence, easily extended to a
block form, computing eigenpairs simultaneously. PLHR can take advantage of
Hermitian positive definite preconditioning, e.g., based on an approximate
inverse of an absolute value of a shifted matrix, introduced in [SISC, 35
(2013), pp. A696-A718]. Our numerical experiments demonstrate that PLHR is
efficient and robust for certain classes of large-scale interior eigenvalue
problems, involving Laplacian and Hamiltonian operators, especially if memory
requirements are tight
- β¦