The Traveling Salesman Problem: Low-Dimensionality Implies a Polynomial Time Approximation Scheme

The Traveling Salesman Problem (TSP) is among the most famous NP-hard optimization problems. We design for this problem a randomized polynomial-time algorithm that computes a (1+eps)-approximation to the optimal tour, for any fixed eps>0, in TSP instances that form an arbitrary metric space with bounded intrinsic dimension. The celebrated results of Arora (A-98) and Mitchell (M-99) prove that the above result holds in the special case of TSP in a fixed-dimensional Euclidean space. Thus, our algorithm demonstrates that the algorithmic tractability of metric TSP depends on the dimensionality of the space and not on its specific geometry. This result resolves a problem that has been open since the quasi-polynomial time algorithm of Talwar (T-04)

Maximum Scatter TSP in Doubling Metrics

We study the problem of finding a tour of $n$ points in which every edge is long. More precisely, we wish to find a tour that visits every point exactly once, maximizing the length of the shortest edge in the tour. The problem is known as Maximum Scatter TSP, and was introduced by Arkin et al. (SODA 1997), motivated by applications in manufacturing and medical imaging. Arkin et al. gave a $0.5$-approximation for the metric version of the problem and showed that this is the best possible ratio achievable in polynomial time (assuming $P \neq NP$). Arkin et al. raised the question of whether a better approximation ratio can be obtained in the Euclidean plane. We answer this question in the affirmative in a more general setting, by giving a $(1-\epsilon)$-approximation algorithm for $d$-dimensional doubling metrics, with running time $\tilde{O}\big(n^3 + 2^{O(K \log K)}\big)$, where $K \leq \left( \frac{13}{\epsilon} \right)^d$. As a corollary we obtain (i) an efficient polynomial-time approximation scheme (EPTAS) for all constant dimensions $d$, (ii) a polynomial-time approximation scheme (PTAS) for dimension $d = \log\log{n}/c$, for a sufficiently large constant $c$, and (iii) a PTAS for constant $d$ and $\epsilon = \Omega(1/\log\log{n})$. Furthermore, we show the dependence on $d$ in our approximation scheme to be essentially optimal, unless Satisfiability can be solved in subexponential time

Travelling on Graphs with Small Highway Dimension

We study the Travelling Salesperson (TSP) and the Steiner Tree problem (STP) in graphs of low highway dimension. This graph parameter was introduced by Abraham et al. [SODA 2010] as a model for transportation networks, on which TSP and STP naturally occur for various applications in logistics. It was previously shown [Feldmann et al. ICALP 2015] that these problems admit a quasi-polynomial time approximation scheme (QPTAS) on graphs of constant highway dimension. We demonstrate that a significant improvement is possible in the special case when the highway dimension is 1, for which we present a fully-polynomial time approximation scheme (FPTAS). We also prove that STP is weakly NP-hard for these restricted graphs. For TSP we show NP-hardness for graphs of highway dimension 6, which answers an open problem posed in [Feldmann et al. ICALP 2015]

Inapproximability of Combinatorial Optimization Problems

We survey results on the hardness of approximating combinatorial optimization problems

Approximation Algorithms for Multi-Criteria Traveling Salesman Problems

In multi-criteria optimization problems, several objective functions have to be optimized. Since the different objective functions are usually in conflict with each other, one cannot consider only one particular solution as the optimal solution. Instead, the aim is to compute a so-called Pareto curve of solutions. Since Pareto curves cannot be computed efficiently in general, we have to be content with approximations to them. We design a deterministic polynomial-time algorithm for multi-criteria g-metric STSP that computes (min{1 +g, 2g^2/(2g^2 -2g +1)} + eps)-approximate Pareto curves for all 1/2<=g<=1. In particular, we obtain a (2+eps)-approximation for multi-criteria metric STSP. We also present two randomized approximation algorithms for multi-criteria g-metric STSP that achieve approximation ratios of (2g^3 +2g^2)/(3g^2 -2g +1) + eps and (1 +g)/(1 +3g -4g^2) + eps, respectively. Moreover, we present randomized approximation algorithms for multi-criteria g-metric ATSP (ratio 1/2 + g^3/(1 -3g^2) + eps) for g < 1/sqrt(3)), STSP with weights 1 and 2 (ratio 4/3) and ATSP with weights 1 and 2 (ratio 3/2). To do this, we design randomized approximation schemes for multi-criteria cycle cover and graph factor problems.Comment: To appear in Algorithmica. A preliminary version has been presented at the 4th Workshop on Approximation and Online Algorithms (WAOA 2006

Asymptotics for Exponential Levy Processes and their Volatility Smile: Survey and New Results

Exponential L\'evy processes can be used to model the evolution of various financial variables such as FX rates, stock prices, etc. Considerable efforts have been devoted to pricing derivatives written on underliers governed by such processes, and the corresponding implied volatility surfaces have been analyzed in some detail. In the non-asymptotic regimes, option prices are described by the Lewis-Lipton formula which allows one to represent them as Fourier integrals; the prices can be trivially expressed in terms of their implied volatility. Recently, attempts at calculating the asymptotic limits of the implied volatility have yielded several expressions for the short-time, long-time, and wing asymptotics. In order to study the volatility surface in required detail, in this paper we use the FX conventions and describe the implied volatility as a function of the Black-Scholes delta. Surprisingly, this convention is closely related to the resolution of singularities frequently used in algebraic geometry. In this framework, we survey the literature, reformulate some known facts regarding the asymptotic behavior of the implied volatility, and present several new results. We emphasize the role of fractional differentiation in studying the tempered stable exponential Levy processes and derive novel numerical methods based on judicial finite-difference approximations for fractional derivatives. We also briefly demonstrate how to extend our results in order to study important cases of local and stochastic volatility models, whose close relation to the L\'evy process based models is particularly clear when the Lewis-Lipton formula is used. Our main conclusion is that studying asymptotic properties of the implied volatility, while theoretically exciting, is not always practically useful because the domain of validity of many asymptotic expressions is small.Comment: 92 pages, 15 figure