Shortest nodedisjoint paths on random graphs
Abstract
A localized method to distribute paths on random graphs is devised, aimed at finding the shortest paths between given source/destination pairs while avoiding path overlaps at nodes. We propose a method based on messagepassing techniques to process global information and distribute paths optimally. Statistical properties such as scaling with system size and number of paths, average pathlength and the transition to the frustrated regime are analyzed. The performance of the suggested algorithm is evaluated through a comparison against a greedy algorithm.
1 Introduction
Among the various computationallyhard constraint satisfaction problems, routing and path optimization have attracted particular attention in recent years due to their nonlocalized nature and interdisciplinary relevance. The nodedisjoint path (NDP) problem on graphs studied here, aims at finding a set of paths linking specified pairs of nodes (communications) such that no two paths share a node; the problem is classified among the NPcomplete class [1] of hard combinatorial problems. This has not only been studied as a purely theoretical problem by mathematicians in the series of graph minors [2] under the name of subgraph homeomorphism problem, but also by practitioners due to its wide applicability to various fields. For instance, in communication systems where the network performance is often strictly related to capacity limits, traffic congestion and the rate of information flow; and in problems of virtual circuit routing where switches located at nodes may become bottlenecks. Moreover, due to their distributive nature NDP is more resilient to failure and represents one aspect of optimal routing where network robustness is the main objective.
One specific communication application where efficient and effective NDP algorithms are essential is in the area of optical networks where transmissions using the same wavelength cannot share the same edge or vertex, hence all communications of the same wavelength must be nonoverlapping (disjoint). Consequently, such an algorithm impacts on the achievable network capacity and transmission rate. In this field of routing and wavelength assignment [3], the objective is to find a routing assignment that minimizes the number of wavelengths used. Different techniques that exploit disjoint paths heuristic algorithms have been proposed to tackle this problem; for instance, greedy algorithms [4, 5], approximations based on rounding integer linear programming formulations [6, 7], postoptimization methods [8], bin packing algorithms [9], various heuristic genetic algorithms such as ant colony optimization [10] and differential evolution [11].
Another important application of NDP is in the design of very large system integrated circuits (VLSI), where one searches for nonoverlapping wired paths to connect different integrated hardware components, to avoid crosspath interference.
Similarly, in wireless adhoc communication networks [12, 13, 14], where each node can act as a router, path overlaps imply signal interference and low transmission quality, whereas longer paths imply poor signal to noise ratio due to multiple relays and higher transmission power; hence the need to consider both path length and transmission overlaps to be minimized is essential for routing problems. Solutions to the NDP problem also provide fault tolerant routes due to the optimal separation of communication paths all over the network, so that if a node (router) fails, as frequently happens in wireless networks due to the mobility of hosts, only few communications will be affected [15, 16]. This feature is particularly important when quality of service (QoS) is one of the main requirements in the set up of a communication network, along with the loadbalancing feature of NDP that prevents network congestion by establishing nonoverlapping routes. This is especially relevant to connectionoriented networks [17] that are strongly affected by node failures and congestion [18].
Practical algorithms for various applications often depend on the specific network topologies considered [19] and mostly focus on the optimization version of the problem, i.e. maximizing the number of paths routed [20]. The satisfiability version of the problem, i.e. whether all paths can be routed successfully without overlap, is not considered; theoretical studies often give bounds to the achievable approximation instead of providing a practical algorithm for individual instances and fail to calculate path lengths and possible overlaps at the same time as part of the optimization process observables. Given that paths are constrained to be contiguous and interaction between paths is nonlocalized, a local protocol is insufficient and global optimization is required. The computational complexity is determined by the fact that such a global optimization problem has to consider all variables simultaneously in order to minimize a cost function with non local interactions between variables.
Unlike other constrained satisfaction problems on networks, NDP has received little attention within the statistical physics community. In this paper we consider a random version of NDP on regular graphs (Reg), Erdős Rényi (ER) [21] and a dedicated type of random graph (RER) described in Section 4, with the aim of testing the efficacy of statistical physicsbased methods derived in the context of spin glass theory [22] such as belief propagation or messagepassing (MP) cavity method [23, 24] as viable alternatives to greedy algorithms; we also study statistical and scaling properties of quantities of interest as a function of network size and number of paths. We study sparse regular, ER and RER random graphs as they are the most interesting for the problem at hand, but the methodology can be easily extended to accommodate other sparsely connected architectures. Clearly, due to the hard constraint of node disjoint paths, typically no solutions would be found in graphs having a non negligible number of nodes with degree . Moreover, graphs with a small number of high degree nodes (hubs) or with high modularity measure, such as scalefree or planar graphs, are not interesting for the nodedisjoint routing problem since when a paths passes through one of these special nodes it leads directly to graph fragmentation, hence frustration. The situation would be very different for constraints on edges instead, but this variant of the problem is left for future work. Finally, the requirement for the graph to be sparse is suggested by restrictions on the validity of the cavity method which is based on fast decaying correlation functions, i.e. a negligible number of loops in the graph.
Numerical simulations indicate that MP outperforms greedy breadthfirst search algorithms not only in finding better solution but also in reaching a higher frustration threshold. Moreover, we find scaling of the expected total length of the NDP as a function of the system size and graph connectivity that goes as with exponent that depends on the type of graph, where is the number of nodes and the number of paths. We find good agreements between theory and simulation data for graphs of average degrees and sizes . Finally, we study statistical properties of physical quantities observed a posteriori, i.e. when a solution is found, such as path length distribution, degree distribution and maximum cluster size for the case of regular graphs.
2 Model
Given an undirected graph (or network) characterized by nodes and edges we define a set of communications as paths on edges of the graph, each of which originates from a source node and terminates in a receiver node . We introduce a variable to characterize each node :
(1) 
The full node characterization is specified by a vector of modulus , where denotes the absolute value of ; is if is neither a sender nor a receiver for any communication, termed a transit node, and if is either a sender or a receiver of some communication. In this way each node can send or receive at most one communication.
For a given set of senderreceiver pairs with we address the problem of finding a set of communications that optimize a cost function which penalizes path length and prevents communications overlap (traffic). The state of the network can be specified by introducing a variable for each edge and for each communication , which specifies whether communication passes through edge and in which direction:
(2) 
Notice that in this formalism . We term these variables currents and define for each edge a vector that collects information on all currents involved in that edge. Currents are subject to Kirchhoff law:
(3) 
For a given path optimization problem we seek the communication configuration that minimizes a cost function , which penalizes path length and traffic congestion:
(4) 
where is a monotonically increasing function of , that penalizes both congestion and path length; where denotes the absolute value of .
We would like now to search for approximate solutions to this problem by messagepassing equations [23]. To derive a distributed algorithm it is useful first to consider treelike graphs , for which one can derive exact recursive equations, and later on use these equations as an approximation for arbitrary graphs .
If is a tree, the removal of any edge divides in two disjoint subtrees and (see figure 1).
We can define as the optimized cost on when the current flows through on the edge ; in this way we can write the message sent from node to his neighbor when the current flows through the edge as .
The messages admit the minsum [23] recursion relation:
(5) 
where the symbol stands for the set of neighbors of node and the constraint is the Kirchoff law (3).
In the following we will use the recursion equation (5) on arbitrary random graphs to approximate the constrained minimum of , the cost defined in equation (4). Namely:
(6) 
where the last subtracted term is introduced to avoid double counting the cost of edge .
Unfortunately, the computational complexity of this algorithm is exponential in the number of communications . In fact, messages can a priori take values corresponding to all possible currents passing through a single edge . Therefore, we cannot generally treat even moderately large values of [25]. The problem can be simplified if we introduce the hard constraint that paths cannot overlap on nodes (and thus neither on edges). This has the important consequence of reducing the configuration space from to and the computational complexity becomes linear in . This restricted version of the path optimization problem is called the nodedisjoint path problem (NDP), as we already mentioned in the introduction, and is the problem we address here. Notice that since we impose the nodedisjoint constraint for the communications, then one communication at most flows through the edges, so that This corresponds to taking:
(7) 
so that the cost function (4) represents indeed the total path length.
In order to solve equation (5) iteratively we define a protocol for taking into account only the allowed configurations at each edge given the current value passing through it and at vertex .
If then:
(8)  
(9) 
If then:
(10)  
(11)  
(12)  
(13) 
The constant that appears equations (9) and (13) are the costs assigned for a unit of current passing through the considered edge (i.e. ). This cost is the one required for the shortest paths but can be generalized to other arbitrary types of costs.
Equation (8) represents the case where is a transit node and no current passes through edge , then the allowed configurations are that either no currents pass through the remaining neighboring edges (first term inside curly brackets) or one current enters and then exits through a pair of neighboring edges, all others edges being unused (second term inside brackets). In figure 2 you can see a diagram representing the different allowed configurations for a transit node. Equation (9) represents the case where is a transit node and the communication passes through edge ; in this case the only allowed configuration is that where the same communication enters/exits from one of the other neighboring edges, all others being unused. Similar considerations are used to formulate the equations (1013) for senders and receivers.
The procedures of applying the algorithm can be summarized as follows:

Initialize messages at random.

Use the converged messages to calculate physical observables.
3 Obtaining a solution
Once the iterative equations (8), (9) and (1013) have converged, the resulting messages can be used to calculate the solution. We define the energy per link [24]:
(14) 
where the last term on the right is subtracted to avoid double counting as it appears in both of the previous two terms. To find a solution we calculate:
(15) 
for each link in the graph and store the current values that minimize the energy per link for each edge:
(16) 
Eventually, we sum over all to find the different paths and total length:
(17) 
In the cavity formalism [24] this is equivalent to calculating the quantity:
(18) 
which represents the energy per node. Finally, the total energy (or path length) is the combination of the two, which in the case of a regular graph establishes the formal relation:
(19) 
Notice that the calculation of is carried out link by link as if the energies per link were statistically independent. It is not intuitively clear that doing this will result in the optimal paths which do not overlap and are also fully connected from the source to the receiver. This is a consequence of having used messages which implicitly contain global information on the constraints and path lengths, so that the energies per link are indeed globally interdependent albeit in a nonobvious manner.
To fully characterize the solutions statistically we calculate also other observables as explained below. Finally, we calculate the paths and corresponding lengths in a sequential order using a greedy breadth firstsearch local algorithm (BFS) to compare the results obtained against our MPbased algorithm.
3.1 Algorithmic complexity
The nodedisjoint constraint is very restrictive and algorithmically helpful in comparison to other routing models where overlaps are allowed but minimized [25, 26]. This hard constraint is indeed paramount in reducing considerably the algorithm’s computational complexity. If we allow for overlaps we need to span a configuration space of the order of at each cavity iteration, leading to a complexity of , where the exponent explores the different flux combinations for each of the independent neighboring sites of the considered message; there are order of such messages. In order to tackle this issue proper approximations have to be introduced as in [25, 26] where they use techniques from polymer physics [25] or convexity properties of the cost function [26]. On the contrary when the overlap is prohibited we reduce the configuration space from to , as this is the number of allowed configurations (the term is derived from the number of possible currents and the additional is due to the configuration of all ), hence there is no need for approximations because the entire configuration space can be efficiently calculated by the cavity equation. Actually, the use of cavity MP implicitly requires one important approximation as it assumes that when node is removed, all its neighbors are statistically independent. This is equivalent of having fast decaying correlation function between these neighboring nodes. This hypothesis is verified in trees and in locally treelike sparse graphs.
For the same reason is important to distinguish between edge and node overlaps. In this work we chose to consider constraints on nodes motivated by the reduced complexity as explained before; in case of edge constraints one has to consider a much bigger configuration space where all configurations with different communications entering and exiting the same transit node must be considered in the optimization routine. For this reason approximations should be introduced as in the case of the models which minimize overlap. The edgedisjoint variant of the problem will be left for future work.
We performed single instance simulations to find optimal microscopic solutions; to obtain macroscopic averages one would usually use population dynamics, one of the most commonly used numerical tools in statistical mechanics literatures [27, 23] for studying similar models. Population dynamics is considered when the thermodynamic limit is taken and the system size is not fixed a priori as in the single instance algorithm. In our model the use of population dynamics does not make much sense since the parameter enters explicitly in the expressions of the messages because it represents the domain of the fluxes, which is of size . But when we fix at the same time we are fixing a system size , because we extract random pairs with density . Hence, it is impossible to decouple the messages domain from the system size, preventing us to properly employ the thermodynamic limit through population dynamics. There is also another problem, that such a macroscopic oriented approach would introduce averages over all possible configurations , including both frustrated and unfrustrated configurations with much higher energies. Thus the macroscopic averages are highly biased by the fewer frustrated configurations and more complex algorithm should be designed to discard such cases. For these two reasons we did not consider in the following the population dynamic counterpart of the algorithm but focused only on averages over single instances.
3.2 Greedy algorithm
To test the performance of the algorithm we compared the results obtained with those given by a greedy algorithm (or its variant) that is often used in literatures to solve the NDP problem in different contexts [12, 13, 4, 5]. The greedy protocol considers only local information around the sources and then builds up a solution step by step recursively, hence reducing considerably the complexity but at the same time completely ignoring other communication positions in the network.
A typical greedy algorithm works in the following way: start by choosing an arbitrary pair , find the shortest path linking the two nodes and then remove nodes belonging to this path from the available network nodes. Choose a second pair and repeat the procedures until either all the paths from sources to destinations have been established or no solution can be found due to frustration. Clearly, the performance of this algorithm is strictly dependent on the order in which we choose the pairs. For instance, in the extreme case the first pair selected is the one with the longest shortest path among all the communications; this implies that we have effectively a more restricted graph and choice of paths, leading for a long second path and even more restricted choice of paths later on.
4 The results
We performed numerical simulations on three types of random graphs. Standard regular random graphs (Reg): each node has fixed degree ; Erdős Rényi random graphs (ER) [21]: edges are drawn at random between each pair of nodes with probability ; a decorated random graph (RER): starting from a regular random graph of degree (which is the minimum degree of this graphs), we then randomly add new edges as in the ER model until the final average connectivity is . Notice that the degree distribution in this case can be obtained by writing where is Poisson distributed with . The parameters used were average degrees and system sizes of We calculated averages over realizations for both the MP and the greedy algorithms; we used a smaller number of realizations for cases of higher complexity (as for and ). Nevertheless, results in all cases are stable and with small error bars with respect to the symbols used. We omitted the error bars from the figures for clarity.
We found a system size scaling that is a cubic function of the variable . A qualitative explanation of the scaling is as follows. The average path length in random graphs goes as (see [28] for an extensive review of graphs properties) and in our case we have paths to consider.We can refine the dependence on using instead . Now, suppose all communications take their shortest path, the quantity would be a good estimate of graph occupancy for the NDP, where the exponent has been introduced as a free parameter to account for the approximation in the expression for as a function of for different types of graphs. Furthermore, if we divide by the number of available nodes we can define the occupancy ratio as . Therefore in this simple case we would expect increasing linearly in . If overlaps are prohibited, for a sufficiently high value of the communications are increasingly forced to take longer routes, leading to a faster than linear increase in the scaling variable . From numerical simulations we found for the NDP a cubic increase in the scaling variable . For small this function agrees well with the linear shortest path behavior but for values of the steeper increase of becomes predominant. Figures 4 and 5 show a good data collapse of the normalized expected total length per node as a function of the scaling variable for different graph connectivities for Reg and ER graphs respectively. We notice a first regime where the curves follow the linear behavior of the dashed line representing the shortest paths. The term “sparse regime” is used since paths are sufficiently far apart, is small, and then no rerouting is needed as each communication will simply take its shortest path. For the curves show the cubic steeper behavior that represents the increase in path lengths to avoid overlaps. The term “dense phase” reflects the increase in path density; is sufficiently high so that shortestpath choices induce conflicting demands and communications are rerouted, taking longer paths to avoid overlaps.
Finally, for large we can identify different frustration points, represented by vertical lines in figure 4, that connect the largest for which solutions have been found with the points where frustration is reached and the length is set to by convention. We see that the frustration points do not collapse and that the bigger the graph size and the higher the connectivity , the earlier frustration sets in (as a functions of ). Arguably, this is due to algorithmic convergence rather than theoretical arguments. In fact the higher and are, the higher the corresponding algorithmic complexity, and hence the larger the number of iterations required to reach convergence. Due to the prohibitive computation cost we ran a smaller number of instances for higher values of and , and without increasing the preset maximum convergence time. We suspect that convergence can be reached in these cases albeit in a much longer times, and hence a solution could be found as well in theory, but has not been found due to the computational limits imposed. Hence we can not provide a precise measure for the frustration transitions nor make further statements regarding their collapses for different systems sizes and connectivities.
In figure 6 we can see the scaling behavior for Reg, ER and RER of given average connectivity and different system sizes; we fixed arbitrarily to highlight the dependence on . We can notice how different types of graph, although having different average lengths, follow the same cubic scaling in . The steeper slope of the ER graphs shows the smaller number of paths choices in this type of graphs that forces the path to rewire in increasingly more convoluted patterns and hence also reach frustration earlier.
We found that our MP algorithm outperforms the greedy BFS both in finding a better solution (smaller ) and in reaching higher values of the frustration transitions. In figures 7, 8 and 9 we plotted the expected normalized total length for both greedy and MP algorithms, for the three graph types, fixed connectivity and different system sizes. We focused on the case for Reg and ER and for RER because of its lower complexity compared to higher ; nevertheless, simulations for different values indeed agree with the suggested scaling and exhibit the same behavior. Initially, in the range where solutions exist the greedy algorithm gives the same total length as the global algorithm up to a certain value of (and ). The explanation is that in this interval the graph is sparse, communications typically do not interact and shortest paths can be selected. This also shows that for a small number of paths the global procedure reduces to act similarly to the greedy algorithm does, e.g. when rerouting is required it involves only two paths, the optimal solution will adopt the shortest path for one and will reroute the second. When increases, we see that the global optimization algorithm outperforms the greedy approach in both finding the optimal solution and in achieving a higher frustration threshold. In the regime where the global algorithm gives a better solution (i.e. shorter total length) we see that it is more efficient to globally reroute paths rather than taking the shortest path of selected paths and adapt the other paths. This means that the optimal solution is not a simple superposition of the first shortest path of the communications, but is a more complex solution.
Figure 10 shows the failure ratio defined as the number of unsuccessful instances (for which a solution is not found) over the total number of realizations as a function of the scaling variable . We notice that the greedy algorithm reaches the frustration point (as a function of ) earlier than the corresponding global MP algorithm, regardless the system sizes or graph type.
This shows that, if a solution exists, a global management of the entire set of communications is required in order to find an optimal solution. Whereas if each communication acts selfishly, seeking the corresponding shortest path, unsolvable overlaps between communications emerge at lower values. Both algorithms show an increased failure rate as the system size increases, presumably due to the unscaled limit on the number of iterations allowed and possibly inherent finitesize effects.
4.1 A posteriori statistics: maximum cluster size and degree distribution (regular graph case).
To better understand the optimization process and characterize the solutions obtained we carried out a statistical analysis of the solution a posteriori. Given the clearer statistical interpretation of the results obtained (due to the limited number of possible connectivity values and their evolution, and the higher frustration threshold for a given connectivity), we chose to study regular graphs. In this case one can gain more insights into the type of routes formed and the reduced effective graphs that emerge for any number of communications. By a posteriori we mean that once a solution was found, by an MP or greedy algorithm, we removed from the graph all nodes and edges taking part in the paths and then calculated statistical properties of the remaining graph . In particular, we calculated the maximum cluster size and the degree distribution.
The existence of a solution to a given set of communications is strictly related to the connectivity of the graph. Each time a solution for a subset of communications is found, edges and nodes involved in the solution paths are effectively removed; and properties of the reduced graph provide information on its ability to accommodate more sourcedestination pairs and the efficiency of the obtained solution in making use of the topology. Figure 11 shows the max cluster sizes ratio of as a function of the scaling variable . This quantity is defined as the ratio between the number of nodes in the maximum connected cluster over the number of nodes in the same graph , the graph obtained after edge and node removal of the obtained solution paths. For both the greedy and the global MP algorithm we see an abrupt step change at some value, between a graph that has a giant connected component and a situation where no solution exists, such that we set the ratio to zero by convention. Moreover, this drop is more abrupt and occurs for smaller values in the case of the greedy algorithm. This means that the greedy procedure does not distribute paths evenly on the graph and creates small disconnected clusters; the greedy algorithm is therefore more sensitive to small changes in connectivity compared to the global MP algorithm, for which the drop is more gradual at first and occurs at higher values. This reconfirms the previous results that the greedy behavior is fragile and sensitive to the position of the communication pairs and the order in which they are selected.
We evaluate the a posteriori degree distribution by calculating the connectivities for the different values, and from these derive the average degree as a function of the scaling variable . Results shown in Figure 12 for different system size and show consistent trends; starting from a regular graph we end up, close to the frustration transition point and after the node and edge removal, with about of the nodes with , whereas have degree and have . The decay of is also plotted for the same process. Also here we see a good data collapse (the different curves can only be distinguished close to frustration).
From graph theory [28] we know that when the graph is likely disconnected, at least to two giant components; the numerical results show that frustration is reached when has value , which corresponds to an average degree , still higher than the connectivity threshold . This can be explained by the fact that tighter constraints on edge availability are imposed in the case of the NDP problem, resulting in frustration even before the graph disconnects (i.e., disconnection is sufficient but not necessary for frustration). Indeed in our model it is insufficient to have just a good number of available links, but they should also constitute clusters of connected links to accommodate new communication paths. Hence the average connectivity value observed at the frustration point of .
4.2 Path length and stretch distribution.
Another interesting quantity to consider is the path length distribution close to the critical threshold, and its comparison with the shortest path distribution. Using the rescaled variable , where is the length per communication, we present in figure 13 the distribution obtained for different system sizes. We see a good data collapse for graphs of different system sizes and connectivities to a Gaussianlike distribution with left fat tails, as confirmed by the logplot on the right panel. This can be explained by the fact that the shorter of the shortest paths are less likely to be rerouted. A graph with a high number of communications exhibits a path length distribution with higher length averages (with respect to the shortest path distribution) as well as higher variances because solution path lengths are more broadly spread. We notice that the left tails are similar for all connectivity values whereas the right tails are broader for lower connectivities close to the frustration point. This can be explained by the fact that short paths are less likely to be rerouted and occur in roughly the same proportion in graphs of different connectivities; hence the similarity in the fat left tails. Regarding the right tails  many paths are rerouted through longer routes by the MPalgorithm, but graphs with higher degree allow for more communications with shorter routes due to the higher routing flexibility they offer.
Figure 14 shows the stretch, defined as the difference between the shortest path and the path length obtained through MP optimization, for close to frustration () for different system sizes. We can see that for graphs of degree only of the communications follow the shortest path, all other communications are routed through longer paths. A higher fraction of shortestpath communications is found for higher connectivity graphs, presumably due to the higher routing flexibility they offer. Looking at the tails we can see that there is a nonnegligible fraction of paths that stretch considerably compared to the average shortest path length.
5 Conclusion
We studied the shortest nodedisjoint path problem on regular, ER and RER random graphs using messagepassing cavity equations. We found that the suggested MP algorithm outperforms the greedy breadthfirst search approach both in finding better solutions (shorter total path length) and in finding solutions for higher values of . This shows that a global strategy is needed to optimally route paths which do not overlap at nodes but also have minimal path lengths. We found a scaling rule for the total length that goes as a cubic function of the occupancy ratio , with varying with the graph topology. This behavior resembles the shortest path length for small but increases faster than linearly for a higher number of paths.
We also studied statistical properties of physical observables a posteriori in the case of regular graphs. We found good data collapses for regular graphs of different system sizes and connectivities for quantities such as maximum cluster size, degree distribution and, length and stretch distributions.
We believe this approach is theoretically interesting due to its relevance to hard combinatorial complexity problems but also offers a new direction for solving important practical routing problem in communication, in particular in optical and wireless adhoc networks and VLSI design. This study offers the first step for realizing the potential in this new direction.
Acknowledgement
This work is supported by the Marie Curie Training Network NETADIS (FP7, grant 290038), the EU FET FP7 project STAMINA (FP7265496) and the Royal Society Exchange Grant IE110151. This work is partially supported by Research Grants Council of Hong Kong (605010 and 604512).
References.
References
 R M Karp. Reducibility among combinatorial problems. Springer, 1972.
 Robertson N and Seymour P D. Graph minors. xiii. the disjoint paths problem. J. of Comb Th, B, 63(1):65–110, 1995.
 Chatterjee B C, Sarma N, Sahu P P, et al. Review and performance analysis on routing and wavelength assignment approaches for optical networks. IETE Technical Review, 30(1):12, 2013.
 Chen C and Banerjee S. A new model for optimal routing and wavelength assignment in wavelength division multiplexed optical networks. In INFOCOM’96. Fifteenth Annual Joint Conference of the IEEE Computer Societies. Networking the Next Generation. Proc. IEEE, volume 1, pages 164–171. IEEE, 1996.
 Manohar P, Manjunath D, and Shevgaonkar RK. Routing and wavelength assignment in optical networks from edge disjoint path algorithms. Comm. Lett., IEEE, 6(5):211–213, 2002.
 Banerjee D and Mukherjee B. A practical approach for routing and wavelength assignment in large wavelengthrouted optical networks. Sel. Ar. in Comm., IEEE, 14(5):903–908, 1996.
 Kolliopoulos S and Stein C. Approximating disjointpath problems using greedy algorithms and packing integer programs. In Bixby R E, Boyd E A, and Rââ osMercado R Z, editors, Integer Programming and Combinatorial Optimization, volume 1412 of Lecture Notes in Computer Science, pages 153–168. Springer Berlin Heidelberg, 1998.
 Belgacem L, Charon I, and Hudry O. A postoptimization method for the routing and wavelength assignment problem applied to scheduled lightpath demands. Eur. J. Op. Res., 232(2):298–306, 2014.
 SkorinKapov N. Routing and wavelength assignment in optical networks using bin packing based algorithms. Eur. J. Op. Res., 177(2):1167–1179, 2007.
 Blesa M and Blum C. Ant colony optimization for the maximum edgedisjoint paths problem. In App. Ev. Comp., pages 160–169. Springer, 2004.
 Storn R and Price K. Differential evolution–a simple and efficient heuristic for global optimization over continuous spaces. J. glob. opt., 11(4):341–359, 1997.
 Sumpter Z, Burson L, Tang B, and Chen X. Maximizing number of satisfiable routing requests in static ad hoc networks. IEEE GLOBECOM 2013 Conference Proceedings, 2013.
 Srinivas A and Modiano E. Minimum energy disjoint path routing in wireless adhoc networks. In Proceedings of the 9th annual international conference on Mobile computing and networking, pages 122–133. ACM, 2003.
 Akkaya K and Younis M. A survey on routing protocols for wireless sensor networks. Ad hoc networks, 3(3):325–349, 2005.
 Li X and Cuthbert L. Stable nodedisjoint multipath routing with low overhead in mobile ad hoc networks. In Modeling, Analysis, and Simulation of Computer and Telecommunications Systems, 2004.(MASCOTS 2004). Proceedings. The IEEE Computer Society’s 12th Annual International Symposium on, pages 184–191. IEEE, 2004.
 Jain S and Das S R. Exploiting path diversity in the link layer in wireless ad hoc networks. Ad Hoc Networks, 6(5):805–825, 2008.
 Tanenbaum A S. Computer Networks 4th EditionâÃ³Â¶. Prentice Hall, 2003.
 X et al MasipBruin. Research challenges in qos routing. Computer communications, 29(5):563–581, 2006.
 Aggarwal A, Kleinberg J, and Williamson D P. Nodedisjoint paths on the mesh and a new tradeoff in vlsi layout. In Proc. of the twentyeighth annual ACM symposium on Theory of computing, pages 585–594. ACM, 1996.
 Chekuri C and Ene A. Polylogarithmic approximation for maximum node disjoint paths with constant congestion. In SODA, pages 326–341, 2013.
 Erdős P and Rényi A. On the evolution of random graphs. Publications of the Mathematical Institute of the Hungarian Academy of Sciences, 5:17–Ã±61, 1960.
 Mézard M, Parisi G, and Virasoro M A. Spin glass theory and beyond, volume 9. World scientific Singapore, 1987.
 Mézard M and Montanari A. Information, physics, and computation. Oxford University Press, 2009.
 Mézard M and Parisi G. The cavity method at zero temperature. J. Stat. Phys., 111(12):1–34, 2003.
 Yeung C H, Saad D, and Wong KY M. From the physics of interacting polymers to optimizing routes on the london underground. Proc. Nat. Ac. Sci., 110(34):13717–13722, 2013.
 Yeung C H and Saad D. Competition for shortest paths on sparse graphs. Phys. Rev. Lett., 108(20):208701, 2012.
 Mézard M and Parisi G. The bethe lattice spin glass revisited. Eur. J. Phys.B, 20(2):217–233, 2001.
 Albert R and Barabási AL. Statistical mechanics of complex networks. Rev. Mod. Phys., 74(1):47, 2002.