FAST LINEAR EXPECTED-TIME ALGORITHMS FOR COMPUTING MAXIMA AND CONVEX HULLS

This paper examines the expected complexity of boundary problems on a set of N points in K-space. We assume that the points are chosen from a probability distribution in which each component of a point is chosen independently of all other components. We present an algorithm to find the maximal points using KN + O(N1-1/K log1/K N) expected scalar comparisons, for fixed K greater-than-or-equal-to 2. A lower bound shows that the algorithm is optimal in the leading term. We describe a simple maxima algorithm that is easy to code, and present experimental evidence that it has similar running time. For fixed K greater-than-or-equal-to 2, an algorithm computes the convex hull of the set in 2KN + O(N1-1/K log1/K N) expected scalar comparisons. The history of the algorithms exhibits interesting interactions among consulting, algorithm design, data analysis, and mathematical analysis of algorithms.