Some "strange behaviors“ of longest_path

asked 2023-07-06 04:37:51 +0200

licheng gravatar image

updated 2023-07-06 05:06:02 +0200

version: Sage 10.

I want to find a longest path starting from a given vertex in a graph. Below is my graph, and I have observed that if the longest path starting from a given vertex is a Hamiltonian path, longest_path responds very quickly. However, if it is not a Hamiltonian path, the response time is slow.

g1=Graph({0: [2, 3, 4, 15], 1: [16, 2, 3, 15], 2: [0, 1, 3, 4], 3: [0, 1, 2, 4], 4: [0, 2, 3, 15], 5: [7, 8, 9, 14], 6: [16, 7, 8, 9], 7: [5, 6, 8, 9], 8: [5, 6, 7, 9], 9: [5, 6, 7, 8], 10: [11, 12, 13, 14], 11: [16, 10, 12, 13], 12: [10, 11, 13, 14], 13: [16, 10, 11, 12], 14: [5, 10, 12, 15], 15: [0, 1, 4, 14], 16: [1, 6, 11, 13]})

Here is the calculation process.

 g1.longest_path(1) # quckily 
 g1.longest_path(14) # long time, about 100s

image description

image description

Additionally, I have noticed that the option backtracking in longest_pathdoes not seem to work when we specify the starting vertex. Looking at the output below, it is 1 or 9, not the longest path starting from vertex 14.

g1.longest_path(14,algorithm="backtrack")

image description

Question 1. Why is there such a significant difference between the two cases?( When a vertex is given and there exists a Hamiltonian path that contains it as the starting vertex, it is puzzling why it is found quickly, while non-Hamiltonian paths take much longer.)

Question 2. Does the backtracking method not work when a starting vertex is given?

edit retag flag offensive close merge delete

Comments

2
  1. This is a NP-hard problem, so the running time can be very long (see https://en.wikipedia.org/wiki/Longest...). Now, explaining the running time differences of the ILP solver in different situations is not obvious. I have some ideas like the fact that it's easier to close the optimality gap when the graph has an hamiltonian path, but this is not sufficient.
  2. The backtrack heuristic algorithm don't take the source/target into account. I don't know why it has been designed this way. One should try to modify the code to fix that.
David Coudert gravatar imageDavid Coudert ( 2023-07-07 08:32:56 +0200 )edit
1

The "backtrack" issues are reported in https://github.com/sagemath/sage/issu...

As for getting exact solutions in practice, I'd recommend to install some state-of-art MILP solver like Gurobi or CPLEX (which are free for academic use), and use them with "solver=" parameter.

Max Alekseyev gravatar imageMax Alekseyev ( 2023-10-27 16:15:25 +0200 )edit