Abstract
Worst-case analysis is a performance measure that is often too pessimistic to indicate which algorithms we should use in practice. A classical example is in the context of the Euclidean Traveling Salesman Problem (TSP) in the plane, where local search performs very well in practice even though it only achieves an \(\Omega (\frac{\log n}{\log \log n})\) worst-case approximation ratio. In such cases, a natural alternative approach to worst-case analysis is to analyze the performance of algorithms in semi-random models. In this paper, we propose and investigate a novel semi-random model for the Euclidean TSP. In this model, called the simultaneous semi-random model, an instance over n points consists of the union of an adversarial instance over \((1-\alpha )n\) points and a random instance over \(\alpha n\) points, for some \(\alpha \in [0, 1]\). As with smoothed analysis, the semi-random model interpolates between distributional (random) analysis when \(\alpha = 1\) and worst-case analysis when \(\alpha = 0\). In contrast to smoothed analysis, this model trades off allowing some completely random points in order to have other points that exhibit a fully arbitrary structure. We show that with only an \(\alpha = \frac{1}{\log n}\) fraction of the points being random, local search achieves an \(\mathcal {O}(\log \log n)\) approximation in the simultaneous semi-random model for Euclidean TSP in fixed dimensions. On the other hand, we show that at least a polynomial number of random points are required to obtain an asymptotic improvement in the approximation ratio of local search compared to its worst-case approximation, even in two dimensions.



Similar content being viewed by others
References
Antoniadis, A., Fleszar, K., Hoeksma, R., Schewior, K.: A PTAS for Euclidean TSP with hyperplane neighborhoods. In: SODA, pp. 1089–1105. SIAM (2019)
Applegate, D.L., Bixby, R.E., Chvátal, V., Cook, W.J.: The Traveling Salesman Problem. Princeton University Press, Princeton (2011)
Arora, S.: Polynomial time approximation schemes for Euclidean traveling salesman and other geometric problems. J. ACM (JACM) 45(5), 753–782 (1998)
Arthur, D., Vassilvitskii, S.: Worst-case and smoothed analysis of the ICP algorithm, with an application to the k-means method. SIAM J. Comput. 39(2), 766 (2009)
Asadpour, A., Goemans, M.X., Madry, A., Gharan, S.O., Saberi, A.: An \({O} (\log n/\log \log n)\)-approximation algorithm for the asymmetric traveling salesman problem. Oper. Res. 65(4), 1043–1061 (2017)
Blum, A., Spencer, J.: Coloring random and semi-random k-colorable graphs. J. Algorithms 19(2), 204–234 (1995)
Brodowsky, U.A., Hougardy, S.: The approximation ratio of the 2-opt heuristic for the Euclidean traveling salesman problem. arXiv preprint arXiv:2010.02583 (2020)
Chandra, B., Karloff, H., Tovey, C.: New results on the old k-opt algorithm for the traveling salesman problem. SIAM J. Comput. 28(6), 1998–2029 (1999)
Christofides, N.: Worst-case analysis of a new heuristic for the travelling salesman problem. Carnegie-Mellon Univ Pittsburgh Pa Management Sciences Research Group, Tech. rep. (1976)
Dantzig, G., Fulkerson, R., Johnson, S.: Solution of a large-scale traveling-salesman problem. J. Oper. Res. Soc. Am. 2(4), 393–410 (1954)
David, R., Feige, U.: On the effect of randomness on planted 3-coloring models. In: Proceedings of the forty-eighth annual ACM symposium on Theory of Computing, pp. 77–90 (2016)
Englert, M., Röglin, H., Vöcking, B.: Worst case and probabilistic analysis of the 2-opt algorithm for the TSP. Algorithmica 68(1), 190–264 (2014)
Englert, M., Röglin, H., Vöcking, B.: Smoothed analysis of the 2-opt algorithm for the general TSP. TALG 13(1), 1–15 (2016)
Feige, U., Kilian, J.: Heuristics for semirandom graph problems. J. Comput. Syst. Sci. 63(4), 639–671 (2001)
Garey, M.R., Graham, R.L., Johnson, D.S.: Some np-complete geometric problems. In: Proceedings of the 8th Annual ACM Symposium on Theory of Computing, pp. 10–22 (1976)
Held, M., Karp, R.M.: A dynamic programming approach to sequencing problems. J. Soc. Ind. Appl. Math. 10(1), 196–210 (1962)
Johnson, D.S., McGeoch, L.A.: 8. The traveling salesman problem: a case study. In: Local Search in Combinatorial Optimization, pp. 215–310. Princeton University Press (2018)
Karlin, A.R., Klein, N., Gharan, S.O.: A (slightly) improved approximation algorithm for metric tsp. In: Proceedings of the 53rd Annual ACM SIGACT Symposium on Theory of Computing, pp. 32–45 (2021)
Karp, R.M.: Probabilistic analysis of partitioning algorithms for the traveling-salesman problem in the plane. Math. Oper. Res. 2(3), 209–224 (1977)
Klein, P.N.: A linear-time approximation scheme for tsp in undirected planar graphs with edge-weights. SIAM J. Comput. 37(6), 1926–1952 (2008)
Kolla, A., Makarychev, K., Makarychev, Y.: How to play unique games against a semi-random adversary: Study of semi-random models of unique games. In: 2011 IEEE 52nd Annual Symposium on Foundations of Computer Science, pp. 443–452. IEEE (2011)
Künnemann, M., Manthey, B.: Towards understanding the smoothed approximation ratio of the 2-opt heuristic. In: International Colloquium on Automata, Languages, and Programming, pp. 859–871. Springer (2015)
Makarychev, K., Makarychev, Y., Vijayaraghavan, A.: Approximation algorithms for semi-random partitioning problems. In: STOC, pp. 367–384 (2012)
Manthey, B., Veenstra, R.: Smoothed analysis of the 2-opt heuristic for the tsp: Polynomial bounds for gaussian noise. In: Algorithms and Computation: 24th International Symposium, ISAAC 2013, Hong Kong, China, December 16–18, 2013, Proceedings 24. pp. 579–589. Springer (2013)
Mömke, T., Svensson, O.: Removing and adding edges for the traveling salesman problem. J. ACM (JACM) 63(1), 1–28 (2016)
Padberg, M., Rinaldi, G.: A branch-and-cut algorithm for the resolution of large-scale symmetric traveling salesman problems. SIAM Rev. 33(1), 60–100 (1991)
Papadimitriou, C.H.: The Euclidean travelling salesman problem is np-complete. Theoret. Comput. Sci. 4(3), 237–244 (1977)
Reinelt, G.: TSPLIB-a traveling salesman problem library. ORSA J. Comput. 3(4), 376–384 (1991)
Roughgarden, T.: Beyond the Worst-Case Analysis of Algorithms. Cambridge University Press, Cambridge (2021)
Spielman, D.A., Teng, S.H.: Smoothed analysis of algorithms: Why the simplex algorithm usually takes polynomial time. J. ACM (JACM) 51(3), 385–463 (2004)
Steele, J.M.: Probability Theory and Combinatorial Optimization. SIAM (1997)
Svensson, O., Tarnawski, J., Végh, L.A.: A constant-factor approximation algorithm for the asymmetric traveling salesman problem. J. ACM (JACM) 67(6), 1–53 (2020)
Traub, V., Vygen, J.: Approaching 3/2 for the s-t-path tsp. J. ACM (JACM) 66(2), 1–17 (2019)
Vijayaraghavan, A., Awasthi, P.: Clustering semi-random mixtures of gaussians. In: ICML, pp. 5055–5064. PMLR (2018)
Yukich, J.E.: Probability Theory of Classical Euclidean Optimization Problems. Springer, Berlin (2006)
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
A previous version of this manuscript was published in the proceedings of the 23rd Conference on Integer Programming and Combinatorial Optimization (IPCO 2022).
Appendices
Appendix 1: Proof of Lemma 17
We start with an auxiliary lemma.
Lemma 23
Let \(\varepsilon >0\). Consider the distinct points \(v,v',x,x'\subseteq \mathbb {R}^2\) and suppose that \(\left\Vert v-v' \right\Vert =\varepsilon \le 1 \), \(\left\Vert x-x' \right\Vert =\ell \le 2\), and
where \([x,x']\) denotes the segment from x to \(x'\). Then the pair of edges \(\{(v,v'), (x,x')\}\) is unordered 2-optimal.
Proof
Without loss of generality, assume \(x=(0,0)\), \(x'=(\ell ,0)\) and v lies in the upper half-plane. Hence, letting \(a=\left\Vert x-v \right\Vert \), \(b=\left\Vert v'-x' \right\Vert \), \(h = d(v, [x,x'])\), and \(\phi (v,v') = a+b-\ell -\varepsilon \), we want to show that
for all choices of \(v,v'\) (\(\phi (v',v) > 0\) follows symmetrically). Note that when \(v, x',\) and \(\epsilon \) are fixed, b is minimized when \(v'\) is on the line joining v and \(x'\), so we may assume it is the case. Note moreover that we can assume that \(v' \in [v,x']\) since otherwise \(\varepsilon ^2> b^2 > h^2\), contradicting the hypothesis of the lemma. Hence, \(b+\varepsilon = \left\Vert x'-v \right\Vert \).
Let t be the horizontal coordinate of v. We first consider the case \(t<0\). If \(t<0\), we have \(h = \left\Vert v-x \right\Vert =a\). Consider the right triangle with vertices \(x,(0,h),x'\). Its hypotenuse has length \(\left\Vert x'-(0,h) \right\Vert \le \left\Vert x'-v \right\Vert = b+\varepsilon \), for \(t<0\). Its catheti are of length \(\left\Vert x-(0,h) \right\Vert =h\) and \(\left\Vert x-x' \right\Vert =\ell \). By Pythagoras’ theorem, we have \(\left\Vert x'-(0,h) \right\Vert ^2 = h^2 +\ell ^2\). Thus
where we used the hypothesis of the lemma.
We have
where the last inequality is since \(\varepsilon \in [0, 1]\), \(h\ge \varepsilon \), and \(\ell \le 2\). We get that \(a+b > \ell +\varepsilon \), which implies \(\phi (v,v')>0\).
Now consider the case \(t\ge \ell \), and let \(v=(t,q)\). Note that in this case \(h=\left\Vert x'-v \right\Vert =b+\varepsilon \). Thus by the hypothesis of the lemma
Consider the right triangle with vertices x, v, (0, t). Its hypotenuse has length \(\left\Vert x-v \right\Vert =a\), while its catheti have length t and q. In turn, the hypotenuse of the right triangle with vertices \(x',(0,t),v\) has length \(\left\Vert v-x' \right\Vert =b+\varepsilon \) while its catheti have length \(t-\ell \) and q. Hence:
Thus,
and \(\phi (v,v')>0\) follows since \(\varepsilon \in [0, 1]\), \(a > \varepsilon \), and \(\ell \le 2\).
Last, assume \(t\in [0,\ell ]\). Hence, \(v=(t,h)\) and we want to find the minimum of
By studying the variations of \(\phi \) as t changes, one finds that it has a minimum at \(t=\frac{\ell }{2}\) with value
and we have \(\phi \left( \left( \frac{\ell }{2},h\right) ,v'\right) >0\) if and only if \(h^2>\varepsilon \ell +\varepsilon ^2\), which we assumed. \(\square \)
We can now prove Lemma 17.
Proof of Lemma 17
Let \(x,y\in K\) and \((v,v')\in P\). Let \(\varepsilon = \left\Vert v-v' \right\Vert \le \delta \) and \(\ell =\left\Vert x-x' \right\Vert \le \sqrt{2}\le \frac{3}{2}\). We have
where the first inequality holds by the convexity of K and the fourth, fifth, and sixth because \(\varepsilon \le \delta \le 1\). Thus by Lemma 23, the swap is impossible. \(\square \)
Appendix 2: Proof of Lemma 18
Proof of Lemma 18
Let \(a = \left\Vert e \right\Vert \) and assume without loss of generality that \(e = (e_1,e_2)=((0,a),(0,2a))\) and \(f=(x,x') = ((x_1,x_2),(x'_1,x'_2))\) (with \(x\ne x'\), \(x_1,x'_1\ge 0\) and \(x_2,x'_2\le 0\)).
Without loss of generality, assume \(x'_1\ge x_1\).
Case 1: \(x_2'\ge x_2\). Assume we have a possible 2-swap between e and f, and rename \(e_1\) and \(e_2\) to \(e_1'\), \(e_2'\) so that, after the swap, \(e_1'\) is linked to x and \(e_2'\) to \(x'\). For any vector u denote \(u_1,u_2\) its horizontal and vertical coordinates. Then,
where last inequality follows from \(x_1'\ge x_1 \ge 0\). Moreover,
where last inequality follows from \(0\ge x'_2 \ge x_2\). Notice that at least one of these inequalities is strict, because either \(x_2'>x_2\) or \(x'_1>x_1\) (otherwise we would have \(x=x'\)). Thus,
so the 2-swap is not possible. This finishes the first case.
Case 2: \(x'_2<x_2\). As before assume we have a possible 2-swap and rename \(e_1,e_2\) to \(e_1'\), \(e_2'\) according to the swap. We then have
where we used \(x_2<0\). Moreover, since
and
we have
so that
which concludes the proof. \(\square \)
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Balkanski, E., Faenza, Y. & Kubik, M. The simultaneous semi-random model for TSP. Math. Program. 206, 305–332 (2024). https://doi.org/10.1007/s10107-023-02011-w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10107-023-02011-w