Equivalent cyclic polygon of a euclidean travelling salesman problem tour and modified formulation

We define a geometric transformation of Euclidean Travelling Salesman Problem (TSP) tours that leads to a new formulation of the TSP. For every Euclidean TSP n-city tour, it is possible to construct an inscribed n-polygon (Equivalent Cyclic Polygon, ECP) such that the lengths of the edges are equal to the corresponding TSP tour links and follow the same sequence order. The analysis of the ECP elicits the possibility of defining a new objective function in terms of angles instead of distances. This modification opens the way to identify characterizing geometric parameters of the TSP as well as to explore new heuristics based on the inclusion of additional constraints. The experimentation with a set of cases shows promising results compared to the traditional compact formulations. The behavior of the ECP-based TSP formulations is better when the nodes of the TSP are randomly or evenly distributed.


Introduction
The Traveling Salesman Problem (TSP) is probably the most studied problem in combinatorial optimization. Biggs, Lloyd and Wilson (1976) (Laporte 1992;Reinelt 1994;Raman and Gill 2017) and applications (Lenstra and Kan 1975;Junger et al. 1995). There have also been several formulations for the problem (Langevin et al. 1990;Orman and Williams 2007;Öncan et al. 2009).
The classic definition of the TSP can be stated as follows: Find the shortest route (tour) for a salesman starting from a given city, visiting each of a specified group of cities, and then returning to the original point of departure (Dantzig et al. 1954). The conceptualization of the problem can be generalized to minimize the duration or the total cost of the tour.
From the Graph Theory approach, the TSP can be stated as the problem of finding the least cost Hamiltonian circuit of a given complete graph K N ≡G(N, A), where N = {1,…,n} is the set of nodes, and A = {(i,j): i,j ∈ N, i = j} the set of arcs (Miller et al. 1960). Usually, the cost of each arc is defined by the cost (or distance) matrix C = {c ij }. Depending on the relationship between c ij and c ji , the problem is defined as Asymmetric TSP (C = C T ) or Symmetric TSP (C = C T ).
Several formulations have been proposed to solve the ATSP (Langevin Padberg and Sung 1991;Orman and Williams 1990;2007). Öncan et al. (2009) provide a thorough analysis of the ATSP formulations, including how they relate to each other. To introduce our work and define the base model formulation, we summarize the main formulations, their basics, and distinctive approaches, and refer to this study for an in-depth analysis.
The work of Dantzig et al. (1954) is considered pioneer in the field. The DFJ formulation uses binary variables to formulate and solve the TSP. The authors propose what is called a natural formulation, in the sense that it does not include any additional variable other than x ij to signify if the arc fromi to j is selected for the tour or not. One essential aspect of TSP formulations is how to ensure that the solution is a single tour and not a set of independent subtours. To this purpose, in the DFJ formulation it is necessary to add a set of constraints that cannot be considered before starting the calculations without a sub-tour identification algorithm. Furthermore, the set of constraints grows exponentially with the number of cities of the TSP, O(2 n ), and, in practice, it is solved with delayed column generation methods.
To overcome these limitations, an alternative to the natural formulation is given by the so-called compact formulations. The basis for this set of formulations is provided by Miller et al. (1960). These authors propose the addition of instrumental variables that determine the order in which the cities are visited (their formulation is even more generic for they allow to visit all the cities in several tours, i.e., Vehicle Routing Problem (VRP); as it is the usual in TSP literature, in the remainder of the paper we will consider their model for the particular case of setting to one the number of tours (Öncan et al. 2009)). Thus, the MTZ formulation allows setting the subtour elimination constraints before solving the problem. Among the family of compact formulations, this approach is called a compact node-ordering formulation. The number of subtour elimination constraints compared to the DFJ formulation, is significantly reduced, O(n 2 ), although they lead to a weak linear programming (LP) relaxation. Desrochers and Laporte (1991) further develop the MTZ formulation and improve the performance of the subtour elimination constraints. Some other relevant compact formulations have been proposed. Gavish and Graves (1978) develop two compact arc-ordering or flow-based formulations with the intention to propose a formulation more adaptable to other transportation scheduling problems. They still use the same number of integer variables as MTZ, but they change the constraints set to include flow variables. Fox et al. (1980) provide another set of time-dependent TSP formulations. They reduce the number of linear constraints from O(n 2 ) to O(n) after including a set of time periods. Gouveia and Pires (2001) propose a compact precedence-based formulation, called the disaggregated MTZ model. They include a new variable that represents if a node is visited in the path between two other nodes.
As aforementioned, when the cost (distance) from one city i to another city j is the same if the arc is reversed (c ij= c ji , ∀ i, j ), the resulting particular case of ATSP is called a Symmetric TSP (STSP). It is the case of the Euclidean TSPs in which the nodes represent points, and the costs are the corresponding Euclidean distances (symmetric by nature). The formulations for the ATSP are directly applicable to the STSP if the costs c ij and c ji are considered as if they were different, i.e., with two arcs of the same cost linking each pair of nodes. Nevertheless, for the natural formulation it is possible to achieve a more efficient STSP formulation if the two arcs connecting each pair of nodes are substituted by one single arc (Padberg and Sung 1991). In this case, the set of arcs will be E = {{i,j}: i,j ∈ N, i < j}. Each arc will have a cost (or distance) c {ij} associated as well as a decision binary variable x {ij} similar to the ATSP formulation. The set of subtour elimination constraints is reduced in the STSP formulation compared to the ATSP formulation, hence the performance improvement, yet it remains exponential with the number of cities, making it necessary to include a subtour identification algorithm. For the MTZ and other compact formulations, "symmetrization" by means of considering one single arc linking each pair of nodes is not possible (Padberg and Sung 1991).
In this paper, we explore a geometric transformation of the Euclidean TSP that might open new possibilities of solving the TSP. The extensive literature on the problem has led to very efficient solving procedures. Nevertheless, the problem is known to be NP-Complete (Papadimitriou 1977), and therefore, for problems with a large number of cities, the computation time becomes critical (Applegate et al. 2006;Dubois-Lacoste et al. 2015). On the other hand, there are numerous optimization problems that are variants and extensions of the TSP (Reinelt 1994;Anbuudayasankar et al. 2014). New approaches to the TSP can then be of great interest to be explored since they might foster applications in different problem variants.
Since the Euclidean TSP can be directly analyzed from a Euclidean geometric perspective, with the nodes representing points in the map and the tours being polygons, the interest of approaching the problem based on a geometric analysis has drawn the attention of researchers. One of the main geometric characteristics of the TSP is the fact that the optimal solution is a non-intersecting polygon (Lawler et al. 1991). Besides the analysis of geometric properties of the optimal solution, many works focus on the proposal of geometric-based heuristic algorithms (Lawler et al. 1991;Reinelt 1994). Some of them, like the largest angle heuristic (Norback and Love 1977) among many others (Asani et al. 2020), are based on the convex hull, which is the minimum convex polygon that encompasses a set of nodes. As computational geometry gained practical importance (Edelsbrunner 1987), there have been proposals based on techniques and concepts from this field aimed at solving the TSP. In particular, some of the methods are based on the application of Voronoi diagrams and Delaunay triangulations (Edelsbrunner 2014;Reinelt 1994).
However, mathematical properties of the TSP and the optimal solution remain open to research. Despite the great attention given to the TSP and its numerous applications and extensions, to the best of our knowledge there has not been a previous attempt of a geometric transformation of the tour polygon, so the problem can be formulated in terms of angles instead of distances. The main contribution of this paper is the proposal of such a geometric transformation of Euclidean Travelling Salesman Problem tours that leads to a modified formulation and new heuristic approaches, as well as to the characterization of intrinsic geometric properties of the solution. Experimentation with a set of problems shows that the proposed geometric transformation might lead to performance improvement in certain problems. The proposed approach opens the way to explore new solving strategies for this abundantly studied and applied problem.
In Sect. 2, we recall the compact node-ordering formulation from Miller et al. (1960). In Sect. 3, the proposed geometric transformation with the definition of equivalent polygonal chains and the equivalent cyclic polygon for every TSP tour is presented, leading to a modified formulation of the problem and new heuristic procedures that are presented in Sect. 4. In Sect. 5, we present the results yielded by the experimentation using a set of cases taken from TSPLIB (Reinelt 1991), which shows the potential practical interest of the proposed transformation. Finally, we summarize the main conclusions in Sect. 6.

TSP base model formulation
As aforementioned, in our work we focus on the geometric transformation of the Euclidean TSP tours and, thus, we deal with the Symmetric TSP. Based on this transformation, we define a modification of the objective function and explore the interest of adding some constraints. With the purpose of isolating the effect of the transformation and the addition of constraints, we want to compare our proposal directly solving a linear programming model that does not involve subtour elimination algorithms. For that reason, we do not take the symmetric DFJ as the base model formulation for the comparison. Among the compact formulations, we select as the base model formulation (Model A) the classic MTZ adapted to the TSP, which is widely known and sets the basis for this group of formulations.
Let us use the classical TSP formulation on a graph G = (N, E), such that N is the set of nodes (cities) N = {1, 2, …, n}, and E = {(i, j)| i, j N, i = j} is the set of links. 1 We use a distance parameter d ij for each (i, j) E, as the distance between nodes i and j. We define a binary variable x ij , which equals 1 if the edge {i, j} is part of the tour, and 0 otherwise.

Indices:
i, j ∈ N Nodes. Parameters: Distance from node i to node j. n = |N | Total number of nodes. Variables: x i j ∈ {0, 1}, ∀(i, j) ∈ E 1 if we go from node i to node j, 0 otherwise. u i ≥ 0, ∀i ∈ N Auxiliary variables that determine the relative order in which node i belonging to the solution is visited.
Objective function (1) minimizes the total length of the tour. Constraints (2) and (3) guarantee that each node is visited only once. The set of constraints (4) prevents the presence of subtours in the solution while determining the sequence of the tour through the values of the auxiliary variables: if a link from node i to node j is in the solution (x ij = 1), then u j has to be greater than u i (for i,j > 1, i.e., taking i = 1 as the starting city of the tour). Constraints (4) constitute the grounds of the compact formulation.
Throughout the document, we will use the same definitions of indices, parameters, and variables of this section, whereas the different variants will concern the objective function and the set of constraints (which we will refer to as models, with Model A being the base model formulation).

Proposed tour geometric transformation and modified TSP formulation
The core of our proposal is the definition of a geometric transformation of any given feasible TSP tour. It is based on the construction of what we name "Equivalent Inscribed Polygonal Chain" (EIPC) and "Equivalent Cyclic Polygon" (ECP), which are described in detail in Sects. 3.1 and 3.2 respectively.

Equivalent inscribed polygonal chain
Given a Euclidean TSP tour of n cities defined by the sequence order in which they are visited {p 1 , p 2 , …, p n }, we define a geometric transformation such that a sequence of vertices {v 1 , v 2 , . . . , v n , v 1 } is inscribed in a circle of a given radius such that each link p k-1 -p k (k = 2, …, n) of the TSP tour has an equivalent edge v k−1 v k of the same length, with the returning link p n -p 1 equivalent to the last edge v n v 1 . We note that in general there will be two vertices (v 1 , v 1 ) corresponding to the same city p 1 . This transformation defines an inscribed polygonal chain (polyline) (Levcopoulos et al. 2002). The radio of the circumscribing circle, or circumcircle, is denoted as R c .
To construct the equivalent inscribed polygonal chain (EIPC), we start positioning v 1 in any point of the circumcircle of radius R c . Then, we define another circle with center in v 1 and radius equal to the length of link p 1 -p 2, and find the intersection with the circumcircle. Since there will be two intersection points, we select v 2 in the circumcircle as the first intersection clockwise from v 1 . We proceed the same way following the sequence of the remaining cities, until finding v 1 as the clockwise circumcircle intersection of the circle centered in v n and radius the length p n -p 1 .
If we define a Cartesian coordinate system (x, y), placing the center of the circumcircle in the origin and v 1 in the vertical axis (0, R c ) (or other specific point (x 1 , y 1 )|x 2 1 + y 2 1 = R 2 c ), the vertices v k of the inscribed polygonal chain will verify the following (there will be two possible intersections for each vertex as abovementioned): Figure 1 shows an example of an 8-city TSP tour defined by the sequence 1-2-8-3-7-4-5-6 (middle of the figure) and two EIPCs. Depending on the length of the radius, the summation of all the central angles defined by the chain edges will be higher than 2π (left) or lower than 2π (right). In the first case (left), non-consecutive edges of the chain will intersect (self-intersecting polygonal chain (Levcopoulos et al. 2002)).

Equivalent cyclic polygon
In the example of Fig. 1, if we progressively decrease the radius of the right chain or increase the radius of the left chain, by continuity (Macnab 1981), there will be a radius in which v 1 will converge to v 1 . The polygonal chain will then close, resulting in what we call Equivalent Cyclic Polygon (ECP) of a TSP tour. In this section, we provide a formal definition of the ECP, a detailed algorithm for its construction as well as for the calculation of the radius of the circumcircle. The definition of the ECP of a TSP tour is a particularization of a more general concept from geometry. Panina and Khimshiashvilib (2012) define a polygonal linkage as "a set of positive numbers L = (l 1 , l 2 , …, l n ) ∈ R n + , which can be realized as the lengths of edges of a closed polygonal line. A d-dimensional configuration of a polygonal linkage is a set of points p 1 , p 2 , …, p n , p i ∈ R d , satisfying the conditions |p i p i+1 | = l i , | p n p 1 | = l n , where | · | denotes the Euclidean distance between two points. Physically, a polygonal linkage is interpreted as a set of cyclically joined rigid rods of length l i . At the locations of joints (hinges), free bending and, thereby, selfintersections and self-overlaps of edges are allowed." Among the set of possible configurations of a polygonal linkage, called the moduli space, we focus on cyclic polygons, i.e., polygons in which all vertices are on a circle. Pinelis (2005), extending and complementing the work of Macnab (1981), proves the existence and isometric uniqueness of a convex cyclic polygon with edges that follow a given sequence of lengths (in which none is larger than the sum of the others). In the case of the proposed geometric transformation of a TSP tour, the proof of Pinelis (2005) assures the existence and uniqueness the ECP, since a tour defines a sequence of lengths. We define the ECP as follows.
Given a Euclidean TSP tour of n cities defined by the sequence order in which they are visited {p 1 , p 2 , …, p n }, it is possible to construct an inscribed n-polygon with vertices in sequence {v 1 , v 2 , …, v n } such that each link p k-1 -p k (k = 2, …, n) of the TSP tour has an equivalent polygon edge v k−1 v k of the same length, with the returning link p n -p 1 equivalent to the last edge v n v 1 . The circumcircle of the ECP is called the mass circle, with radius R m called the mass radius. Figure 2 represents the ECP of the example of Fig. 1 and illustrates the basic relationship between R m , the length of each edge, and the central angle that it delimits. Let us consider the link from node i to node j belonging to a tour. Then the relation of Eq. (7) (i = 1 and j = 2 in Fig. 2), is verified for α ij ≤ π. Summation of all the central angles of the ECP will equal 2π (Eq. 7). Thus, given a TSP tour and the set of distances of the links {d ij , x ij = 1} such that all the central angles verify α ij ≤ π with d max = max i j (d i j ), x i j = 1, we can calculate R m by solving the transcendental Eq. (7) either with the aid of mathematical software or through an iterated procedure.
To gain insights into the existence and uniqueness of the ECP, and in an equivalent approach to Pinelis (2005), we define the function F 1 in Eq. (8), such that Eq. (7) can be expressed as F 1 (R m ) = 0.
F 1 is continuous, differentiable, and strictly decreasing in the interval [d max /2, ∞). Besides, F 1 converges asymptotically to -π as R c tends to infinity. Then, if a = F 1 (d max /2) ≥ 0, there will be a unique value of R m |F 1 (R m ) = 0. Figure 3 shows the representation of F 1 (R c ) for the example of Fig. 2.
If the tour has the particularity of having a very large link, it might happen that Eq. (7) cannot be satisfied even for the minimum possible mass radius R m = d max /2. In this case, the EIPC does not close even for the minimum R m , for which the corresponding edge with maximum length would be the diameter of the circumcircle. The central angle of the longest edge would equal π, whereas the sum of the rest of the central angles would be less than π. From Eq. (8) this case would then verify that a < 0, implying F 1 < 0 in all its domain. Conceptually, the fact that a < 0 means that, for the ECP to be closed, all of its vertices have to belong to a semicircle. This can be stated as the central angle of the longest link being larger than π. Figure 4 illustrates this For most of the TSP tours a ≥ 0 and the construction of the ECP is identical to the procedure described in the former subsection for the construction of the EIPC. Yet, the possibility of a < 0 requires introducing a modification as follows. For the longest link, the vertex must be determined by taking the counterclockwise intersection. Thus, the central angle will be 2π − 2 · arcsin dmax 2R m . Subtracting the central angle from the summation and adding the counterclockwise angle (divided by 2) in Eq. (7) yields the condition that must verify R m (Eq. 9).
F 2 is differentiable in the interval [d max , ∞) and converges asymptotically to 0 + as R c tends to infinity. A detailed analysis of the function, its first derivative and curvature, shows that F 2 tends to infinity as R c tends to d max /2 (F 2 has vertical slope for R c = d max /2), and has a unique zero which corresponds to a F 2 maximum. If a < 0, F 2 will have different signs in the extremes of the domain, asymptotic convergence to 0 + , and there will be a unique value of R m | F 2 (R m ) = 0 (see Pinelis (2005) for an equivalent thorough proof). Figure 5 shows the representation of F 2 (R c ) for the example of Fig. 4. We can then define the following procedure to construct the ECP of a TSP tour, based on the EIPC construction as described in Sect. 3.1: Step 1 Define the sequence of cities such that the link p 1 -p 2 is the longest.
Step 3 Find R m .
Step 4 Construct the ECP.
• If a ≥ 0. Proceed as described in the construction of the EIPC.
• If a < 0. Proceed as described in the construction of the EIPC, but with the following modification: select v 2 in the circumcircle as the first intersection counterclockwise from v 1 .
We can gain some insights on the ECP and the TSP with many nodes, analyzing the case in which the distances between the cities of the tour are of equal length d. The ECP in this case will be a regular polygon. Let us consider a parameter γ > 1, that will represent the relation between the length of the mass circle and the length of the ECP of an n-city TSP as defined by Eq. (11). As the number of number of cities of the problem grows, the length of the tour will be closer to the length of the mass circle (i.e., γ → 1).
The definition of the parameter γ can be generalized, so it can be applicable to any TSP tour as the inverse relation of the length of the tour and the length of its mass circle. In general, for a TSP with many equally distributed nodes, γ will be close to 1. A value of γ significantly higher than 1 will indicate the presence of large links in the TSP tour. Thus, γ can be used as one parameter that characterizes the topology of a TSP tour (Eq. 12) and, moreover, to characterize the TSP as follows. The optimal TSP tour allows us to construct its corresponding ECP*, with mass radius R m *. Thus, if the parameter γ is calculated for the optimal tour, it will then be a characterizing parameter of the TSP.

ECP-based TSP formulation and heuristics
The geometric transformation of the TSP tours leads to a modified formulation of the problem. As the circumradius of a tour EIPC tends to infinity, the length of each edge converges to the corresponding arc of the circumcircle. For a large enough circumradius, the longer the tour, the greater the total central angle in its EIPC. Therefore, the objective function of the TSP (Eq. 1) can be substituted by the sum of EIPC half central angles for a circumference of a large enough circumradius (R L c ), yielding an equivalent formulation (13). If we express the circumradius as a factor times the minimum possible radius to include any link of the TSP, it is verified that: R L c = q · d max /2, q 1. Instead of the standard approach of minimizing the sum of costs or distances (1), the objective function can be formulated as a minimization of the sum of the half central angles of a transformed tour (13). Taking the same indices, parameters, variables, and constraints (Eqs. 2-4) of the formulation of Miller et al. (1960) included in Sect. 2, we can define the new ECP-based TSP (ECPTSP) formulation as follows.
The proposed modified formulation opens the way to define a new family of heuristics for the TSP, through the specification of different values for the circumradius and the addition of EPC geometric constraints. In the following subsections, we propose three heuristics based on the geometric transformation that are modified formulations of the MTZ base model formulation (Sect. 2).

Model B-heuristic mass radius (ECPTSP_hmR)
The first heuristic strategy consists in reducing the radius of the ECPTSP formulation objective function (Eq. 13), so the formulation helps to guide the algorithm by avoiding large links (both explicitly and implicitly) while aligning with the objective of the minimum distance. This is done by solving the TSP with a fast heuristic and taking the mass radius of the ECP of the heuristic solution (R H m ) as the radius for the objective function. For the experimentation of Sect. 5, we use the nearest neighbor heuristic -starting from an arbitrary city, iteratively select the city that is closest to the last city included in the tour-, which is simple and fast, although it is known to yield low quality solutions with large returning links (Rosenkrantz et al. 1977). In fact, for all cases of the experimentation (Sect. 5) it is verified that R H m = q · d max /2, q > 1. To complete the modified formulation, it is also necessary to add the set of constraints (15) that express the impossibility for an edge to surpass the diameter of the mass circle, assuring the range of validity of the arcsine functions in the objective function.
Such that.
We underline the elimination of links that derives from Eq. (15) and its practical application to reduce the computation time for large problems. On the contrary, if the topology of the problem is such that the optimal solution includes large links, ECPTSP_hmR formulation might not provide quality solutions, particularly if a better heuristic is employed. In general, as we develop in the remainder of the paper, there is a trade-off between properly accelerating the search and leaving aside good solutions.

Model C-summation of arcs constraint (ECPTSP_summarcs)
The ECPTSP formulation also opens the possibility of adding a constraint that might improve the calculations for complex problems. The idea behind the additional cut lies in the fact that, in general, for problems of large size, tours in which the total central angle of the circumcircle of the heuristic mass radius exceeds 2π are likely to be low quality solutions. The ECPTSP_summarcs formulation (Model C) results from the addition of constraint (16) to the ECPTSP_hmR formulation.
Model C: such that

Model D-Small Circumradius (ECPTSP_scR)
The next step in the trade-off between reduction of the search space and not eliminating good solutions is explored by drastically reducing the radius of the ECPTSP objective function. Instead of using a large R c L in the objective function (13), we can use what we call a small circumradius R c S that is smaller than half the maximum distance R S c = q · d max /2, q < 1. The formulation is completed by adding a set of constraints equivalent to (15).
Model D: To have the possibility of achieving the optimal value, R c S can be reduced as long as no arc of the optimal solution is unable to satisfy constraints (18). Lowering R c S might also lead to avoiding the optimum yet yielding a feasible solution. The challenging issue is how to set the value of R c S low enough to significantly reduce the number of edges while not making it hard to find a feasible quality solution. It is immediate to determine that it should never be smaller than half the minimum distance from a node to the others, since, if it is not, any link including that node will not satisfy Eq. (18) (it will not fit into the circumcircle). Finding a suitable value for R c S is an issue that arises from the proposed ECPTSP-scR formulation and is subject to research. Noteworthy, due to the reduced value of R c S , the objective function value for most tours will likely exceed 2π (as in Fig. 1 left). For the experimentation of Sect. 5, we take R c S as the longest among the minimum distances from each node to the others R S c = max j (min i d i j ), which, for all the cases selected for experimentation verify R S c = q · d max /2, q < 1.

Experimentation
Aiming to assess the practical interest of the proposed geometric transformation as well as the different heuristic formulations presented in the former section, we compare the following 4 formulations of the TSP: • Model A. Classical formulation of Miller et al. (1960)

Experimental settings
To perform the above comparison, from TSPLIB (Reinelt 1991) we select the 40 TSPs in the Euclidean plane and symmetric (Symmetric EUC_2D type) with a maximum of 400 nodes. Due to the nature of the proposed geometric transformation, it is necessary to compute float Euclidean distances when solving the TSPs. We note that this is not the standard stablished for the TSPLIB instances, in which distances are rounded to integer values. In order to have a comparison benchmark, for each TSPLIB instance, the floating Euclidean distance of the known optimal tour has been calculated, even though for some cases those tours are not optimal.
The experimentation has been performed using a computer with Intel(R) Core (TM) i3-8100 CPU @ 3.60 GHz. We formulate all models in Python 2.7, using the Pyomo package to code the formulation (Hart et al. 2017), the Gurobi 9.1 solver (Gurobi 2021a) to find a solution, and the R programming language to plot the nodes and results (R Core Team 2020). For all runs, a time limit of 10,000 s and a Gurobi optimality MIPGap (Relative MIP optimality gap, Gurobi 2021b) of 0.01% was set for each run. MIPGap provides a relative maximum distance to the optimal value. It is the MIP primal-dual relative gap, calculated as the absolute difference between the incumbent objective function value (primal upper bound) and the incumbent dual objective function value (lower bound) of a Mixed Integer Programming (MIP) problem, relative to the incumbent primal objective value (Gurobi 2021b). Table 1 shows the results divided in three groups. Table 1a includes the results for those instances in which all formulations consumed the time limit without achieving the MIPGap (time = 10,000 s for all runs). Table 1b includes the results in which at least one formulation achieved the MIPGap without reaching the time limit and at least one consumed the time limit (therefore it specifies the computation time). Finally, Table 1c includes those cases in which all formulations reached the MIPGap of 0.01%. The gap (%) column represents the absolute percentage relative deviation of the solution yielded by each model with respect to the optimal TSPLIB tours in terms of distance. Since in general the benchmark tours are not optimal when taking floating distances, gaps become negative in some cases, meaning the solution found is better than the optimal tour with integer distances. Additionally, the table shows the value of the objective function for each formulation: the length of the best solution for the formulation A and the angle of the respective best solutions in radians for the ECPTSP based formulations B, C, and D. In models (B) and (C), the distance of the objective function to π provides an indicator of the improvement with respect to the heuristic taken to calculate R m H . Results are complemented with the value of a parameter θ that characterizes the degree of dispersion of the nodes that define each instance. The details of its calculation and implications will be discussed in Sect. 4.3. Table 2 presents the characterization of each TSP based on the parameters discussed throughout the paper: the optimal length, optimal mass radius, parameter γ, upper bound radius used, radius of the convex hull, small radius, and half the maximum distance.

Discussion and influence of the topology of the problem
In this section, we present a general comparison of the experimentation. Then, we explain a selected parameter to classify the dispersion of nodes, and finally, we use it to group the sets of nodes and plot them to draw some experimentation conclusions.
In Table 3, we compare the results obtained from the forty TSPLIB selected cases for the four formulations. For each formulation, the table shows in how many cases it performs better than each one of the other formulations, alone or tied with another formulation. Formulations (B) and (D) yield better results than (A) in the majority of cases, whereas (C) and (A) are even. Model (D) shows the best performance overall, with 23 cases yielding the best solution, yet the drastic arc elimination strategy leads to infeasibility twice, and in other 4 cases (D) is not able to provide a solution before reaching the time limit. It is therefore a promising approach, but further research concerning the value of is R c S needed. It is worth highlighting that the compact formulation (A) only achieved the best results in 8 out of the 40 cases. In conclusion, there does not seem to be a dominant strategy although ECPTSP based formulations show good results in the majority of cases. Nevertheless, these new formulations do not perform well in some of the cases. With the purpose of identifying in which problems the proposed modified formulations have more potential, we analyze the influence of the topology of the problem.
With this aim and following the work of Dry et al. (2012), we calculate a parameter θ (R in the original) that characterizes the degree of dispersion of the nodes that define each instance. This parameter leads to a characterization of the set of nodes as clustered, random, or evenly distributed. The parameter adopted is proposed by Clark and Evans (1954) who propose using the distribution of distances from each node to its nearest neighbor as a measure of the spatial dispersion of a set of nodes. To calculate θ , the observed mean nearest neighbor distance θ o is divided by the mean θ e Fig. 8 Sets of nodes from TSPLIB type EUC 2D with less than 400 nodes where θ < 0.8 of a Poisson process probability density function that describes the nearest neighbor distance, as expressed in Eq. (19) where A indicates the area.
The statistic θ varies between 0 and 2.149, aiming to classify the distribution of the sets of points from highly clustered to highly regular. It represents how long is the average distance to the nearest neighbor compared to a random set of nodes. The closer to 1 is θ the more randomly are the nodes distributed; if θ is close to 0, the nodes are highly clustered; if θ is high it means the average distance to the nearest neighbor is proportionally high, and then the nodes are evenly distributed.
We classify the 40 TSPLIB instances according to their θ value and, similarly to Dry et al. (2012) we define three groups, respectively, {θ < 0.8}, {0.8 ≤ θ ≤ 1.2}, and {1.2 < θ }. Figures 3, 4, and 5 show the plots of the three groups ordered by their θ value. We use a graphic code to signify the performance of the ECPTSP formulations: if formulations (B, C, and D) obtain better results than formulation (A) in terms of gap solution or computing time when they provide the same gap solution, the plot rectangle edges are continuous and dots are dark green colored; dashed edges and orange dots signify that at least one ECPTSP formulation (B, C and D) obtain better results than formulation (A) or all four models have identical results, and dotted edges and red color is used when formulation (A) performs the best. Figure 6 shows that for high values of θ , ECPTSP formulations (B, C or D) yield satisfactory results. In six of the nine instances the formulation (A) performs the worst. In the results of Table 1 we see how in cases eil51, eil101, eil76, rat195, and rat99, either gaps or computing times decrease substantially with these formulations.
In Fig. 7, there are ten sets in continuous edge (green), ten sets in dashed (orange), and only one set in dotted (red). Regarding this last case, pr76, we observe in Table 1 a similar computation time for all formulations, although (A) performs slightly better than (B, C, and D). Sets st70, rd400, kroA150, gil262, and kroA200 behave the best when using model (D), and ch130 solution improvement for models (B) and (C) is significative. Figure 8 reveals that high clusterization has a negative impact on the performance of models (B, C, and D). The distribution of the nodes presents peculiar patterns. Only in one of the ten instances, the three ECPTSP formulations show better behavior than (A) and in four of the cases (A) yields the best results. Model (D) may be infeasible by using the selected R c S . On the other hand, if θ >1, model (D) performs better than model (A) in all cases except in ts225, where a regular pattern is identified, and in pr76, where the performance of all models results in similar.

Conclusions
Through a geometric transformation based on polygonal linkages, we can inscribe each transformed TSP tour as an open polygonal chain in a circle of a given radius (Equivalent Inscribed Polygonal Chain, EIPC) as well as find its Equivalent Cyclic Polygon (ECP). This transformation leads to a modified formulation for the extensively studied classic Euclidean TSP in which the objective function is expressed as a summation of angles instead of distances.
The modified formulation leads in turn to the definition of a new family of heuristic models (ECPTSP formulations). Experimentation with forty cases of TSPLIB shows that using the ECPTSP formulations may help to obtain better or optimal solutions faster than compact formulations. Using ECP properties mixed with a drastic link elimination shows promising results in some cases. We also verify that a parameter based on the distribution of the nearest neighbor distances is of practical interest in TSP topology characterization. Particularly, it shows that the formulations based on the proposed tour geometric transformation perform better when the nodes are randomly or evenly distributed, whereas it tends to yield worse results in highly clustered problems.
The performance of the heuristic formulation reveals the interest of exploring ways to set the small radius parameter so that it improves the performance of the solving procedure while not incurring in infeasibility or leading to high computation times to provide a first feasible solution. The definition of the ECP opens the way for further research on TSP, exploring new properties and new heuristics. The number of formulations, and TSP extensions such as the Vehicle Routing Problem (VRP), also bring the interest of investigating the applicability of the ECP transformation to achieve efficient results in other models and variants.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/ by/4.0/.