Smoothed analysis

A randomly generated bitmap does not resemble typical pictures.

A typical picture does not resemble a random bitmap.

In theoretical computer science, smoothed analysis is a way of measuring the complexity of an algorithm. Since its introduction in 2001, smoothed analysis has been used as a basis for considerable research, for problems ranging from mathematical programming, numerical analysis, machine learning, and data mining. It can give a more realistic analysis of the practical performance (e.g., running time, success rate, approximation quality) of the algorithm compared to analysis that uses worst-case or average-case scenarios.

Smoothed analysis is a hybrid of worst-case and average-case analyses that inherits advantages of both. It measures the expected performance of algorithms under slight random perturbations of worst-case inputs. If the smoothed complexity of an algorithm is low, then it is unlikely that the algorithm will take a long time to solve practical instances whose data are subject to slight noises and imprecisions. Smoothed complexity results are strong probabilistic results, roughly stating that, in every large enough neighbourhood of the space of inputs, most inputs are easily solvable. Thus, a low smoothed complexity means that the hardness of inputs is a "brittle" property.

Although worst-case analysis has been widely successful in explaining the practical performance of many algorithms, this style of analysis gives misleading results for a number of problems. Worst-case complexity measures the time it takes to solve any input, although hard-to-solve inputs might never come up in practice. In such cases, the worst-case running time can be much worse than the observed running time in practice. For example, the worst-case complexity of solving a linear program using the simplex algorithm is exponential, although the observed number of steps in practice is roughly linear. The simplex algorithm is in fact much faster than the ellipsoid method in practice, although the latter has polynomial-time worst-case complexity.

Average-case analysis was first introduced to overcome the limitations of worst-case analysis. However, the resulting average-case complexity depends heavily on the probability distribution that is chosen over the input. The actual inputs and distribution of inputs may be different in practice from the assumptions made during the analysis: a random input may be very unlike a typical input. Because of this choice of data model, a theoretical average-case result might say little about practical performance of the algorithm.

Smoothed analysis generalizes both worst-case and average-case analysis and inherits strengths of both. It is intended to be much more general than average-case complexity, while still allowing low complexity bounds to be proven.

History

ACM and the European Association for Theoretical Computer Science awarded the 2008 Gödel Prize to Daniel Spielman and Shanghua Teng for developing smoothed analysis. The name Smoothed Analysis was coined by Alan Edelman. In 2010 Spielman received the Nevanlinna Prize for developing smoothed analysis. Spielman and Teng's JACM paper "Smoothed analysis of algorithms: Why the simplex algorithm usually takes polynomial time" was also one of the three winners of the 2009 Fulkerson Prize sponsored jointly by the Mathematical Programming Society (MPS) and the American Mathematical Society (AMS).

Examples

Simplex algorithm for linear programming

The simplex algorithm is a very efficient algorithm in practice, and it is one of the dominant algorithms for linear programming in practice. On practical problems, the number of steps taken by the algorithm is linear in the number of variables and constraints. Yet in the theoretical worst case it takes exponentially many steps for most successfully analyzed pivot rules. This was one of the main motivations for developing smoothed analysis.

For the perturbation model, we assume that the input data is perturbed by noise from a Gaussian distribution. For normalization purposes, we assume the unperturbed data [math]\displaystyle{ \bar{\mathbf A} \in \mathbb{R}^{n\times d}, \bar{\mathbf b} \in \mathbb{R}^n, \mathbf c \in \mathbb{R}^d }[/math] satisfies [math]\displaystyle{ \|(\bar{\mathbf a}_i, \bar{b}_i)\|_2 \leq 1 }[/math] for all rows [math]\displaystyle{ (\bar{\mathbf a}_i, \bar{b}_i) }[/math] of the matrix [math]\displaystyle{ (\bar{\mathbf A}, \bar{\mathbf b}). }[/math] The noise [math]\displaystyle{ (\hat{\mathbf A}, \hat{\mathbf b}) }[/math] has independent entries sampled from a Gaussian distribution with mean [math]\displaystyle{ 0 }[/math] and standard deviation [math]\displaystyle{ \sigma }[/math]. We set [math]\displaystyle{ \mathbf A = \bar{\mathbf A} + \hat{\mathbf A}, \mathbf b = \bar{\mathbf b} + \hat{\mathbf b} }[/math]. The smoothed input data consists of the linear program

maximize

[math]\displaystyle{ \mathbf{c^T} \cdot \mathbf{x} }[/math]

subject to

[math]\displaystyle{ \mathbf A \mathbf{x} \leq \mathbf b }[/math].

If the running time of our algorithm on data [math]\displaystyle{ \mathbf A, \mathbf b, \mathbf c }[/math] is given by [math]\displaystyle{ T(\mathbf A, \mathbf b,\mathbf c) }[/math] then the smoothed complexity of the simplex method is

[math]\displaystyle{ C_s(n,d,\sigma) := \max_{\bar{\mathbf A}, \bar{\mathbf b}, \mathbf c} ~ \mathbb{E}_{\hat{\mathbf A},\hat{\mathbf b}}[T(\bar{\mathbf A} + \hat{\mathbf A}, \bar{\mathbf b} + \hat{\mathbf b}, \mathbf c)] = {\rm poly}(d,\log n, \sigma^{-1}). }[/math]

This bound holds for a specific pivot rule called the shadow vertex rule. The shadow vertex rule is slower than more commonly used pivot rules such as Dantzig's rule or the steepest edge rule but it has properties that make it very well-suited to probabilistic analysis.

Local search for combinatorial optimization

A number of local search algorithms have bad worst-case running times but perform well in practice.

One example is the 2-opt heuristic for the traveling salesman problem. It can take exponentially many iterations until it finds a locally optimal solution, although in practice the running time is subquadratic in the number of vertices. The approximation ratio, which is the ratio between the length of the output of the algorithm and the length of the optimal solution, tends to be good in practice but can also be bad in the theoretical worst case.

One class of problem instances can be given by [math]\displaystyle{ n }[/math] points in the box [math]\displaystyle{ [0,1]^d }[/math], where their pairwise distances come from a norm. Already in two dimensions, the 2-opt heuristic might take exponentially many iterations until finding a local optimum. In this setting, one can analyze the perturbation model where the vertices [math]\displaystyle{ v_1,\dots,v_n }[/math] are independently sampled according to probability distributions with probability density function [math]\displaystyle{ f_1,\dots,f_n : [0,1]^d \rightarrow [0,\theta] }[/math]. For [math]\displaystyle{ \theta = 1 }[/math], the points are uniformly distributed. When [math]\displaystyle{ \theta \gt 1 }[/math] is big, the adversary has more ability to increase the likelihood of hard problem instances. In this perturbation model, the expected number of iterations of the 2-opt heuristic, as well as the approximation ratios of resulting output, are bounded by polynomial functions of [math]\displaystyle{ n }[/math] and [math]\displaystyle{ \theta }[/math].

Another local search algorithm for which smoothed analysis was successful is the k-means method. Given [math]\displaystyle{ n }[/math] points in [math]\displaystyle{ [0,1]^d }[/math], it is NP-hard to find a good partition into clusters with small pairwise distances between points in the same cluster. Lloyd's algorithm is widely used and very fast in practice, although it can take [math]\displaystyle{ e^{\Omega(n)} }[/math] iterations in the worst case to find a locally optimal solution. However, assuming that the points have independent Gaussian distributions, each with expectation in [math]\displaystyle{ [0,1]^d }[/math] and standard deviation [math]\displaystyle{ \sigma }[/math], the expected number of iterations of the algorithm is bounded by a polynomial in [math]\displaystyle{ n }[/math], [math]\displaystyle{ d }[/math] and [math]\displaystyle{ \sigma }[/math].

See also

• Average-case complexity
• Pseudo-polynomial time
• Worst-case complexity