1

我一直在尝试解决路由问题，例如TSP（旅行推销员问题），但有些曲折。我已经使用线性规划 ( CPlex library) 和有向图（带有 Origin 顶点） ( Coin-Or::Lemon library) 对问题进行了建模。

线性程序找到解决方案，但我的问题在于如何检索路径。我可以迭代图形的每个顶点和边以找出线性程序正在使用哪个，所以我想我应该从原点开始并使用选择的边到达下一个节点，直到我再次到达原点。

问题是 CPlex 路径可能不止一次地穿过任何顶点。所以我决定使用递归算法。但我无法完全弄清楚。这是我尝试过的：

FindPath ( node N, path P)
    AddedAnyPath <- false
    for each edge E out of N that hasn't been used yet
        T is the target of E
        add T to P

        if findPath ( E, P )
            AddedAnyPath = true
    end for each
    
    if AddedAnyPath 
        if all edges were used
            return true
        else
            return false
        end if
    endif
end

该算法无法找到路径。如果您能指出我的任何方向，我将不胜感激。

编辑

提供的答案帮助我找到了答案。这是中的代码C++：

bool findPath2 ( Store_Instance &T, DiNode &current, list <DiNode> &path, ListDigraph::ArcMap<bool> &usedArcs, IloCplex &cplex, ListDigraph::ArcMap<IloNumVar> &x) {
DiNode currentNode = current;
bool success = false;
int positionToInsert = 1;
while (true) {
    //Find an unvisited edge E whose value is 1.
    Arc unvisitedArc = INVALID;
    for(Digraph::OutArcIt a(T.g, currentNode); a != INVALID; ++a) {
        if(cplex.getValue(x[a]) >= 1-EPS && !usedArcs[a]) {
            unvisitedArc = a;
            break;
        }
    }
    if (unvisitedArc != INVALID) {
        //Mark edge as visited
        usedArcs[unvisitedArc] = true;
        //Get the target T of the edge
        DiNode target = T.g.target(unvisitedArc);
        //Add Edge E to the path
        list<DiNode>::iterator iterator = path.begin();
        advance(iterator, positionToInsert);
        path.insert(iterator, target);
        //Increase the iterator
        positionToInsert++;
        //If all the edges whose value is 1 are visited, stop
        bool usedAllEdges = true;
        for (ArcIt a(T.g); a != INVALID; ++a){
            if (cplex.getValue(x[a]) > 1-EPS && usedArcs[a] == false) {
                usedAllEdges = false;
                break;
            }
        }
        if (usedAllEdges) {
            success = true;
            break;
        }
        //Else, Set N to be T and repeat
        else currentNode = target;
    } else {
        //No arcs were found. Find a node from the path that has unvisited Arc
        positionToInsert = 0;
        DiNode target = INVALID;
        for (list<DiNode>::const_iterator iterator = path.begin(), end = path.end(); iterator != end; ++iterator) {
            DiNode v = *iterator;
            for(Digraph::OutArcIt a(T.g, v); a != INVALID; ++a) {
                if(cplex.getValue(x[a]) >= 1-EPS && !usedArcs[a]) {
                    target = v;
                    break;
                }
            }
            positionToInsert ++;
            if (target != INVALID) break;
        }

        if (target != INVALID) {
            //cout << "found lost node" << endl;
            currentNode = target;
        } else {
            success = false;
            break;
        }
    }
}
return success;
}
4

1 回答 1

1

TSP（路径）的解决方案是一个有序的顶点列表，在原点开始和结束，访问每个顶点。有两种可能的情况。

在此处输入图像描述

如果您允许访问顶点两次，那么在您的变体中，您允许“子游览”。（案例二） 在此处输入图像描述

您将拥有变量 X_nt，它将为 1。（连接节点 n 到节点 t 的边缘变量）。

如果你没有subtours：案例1，你可以在回到原始节点后停止。

看起来您允许子游览（多个周期，上面的案例 2）。如果是这样，您必须访问作为 TSP 解决方案一部分的所有边，其值为 1。（修改您的代码以跟随从每个节点到其终端节点的任何一条边，一次一条边。）

Let starting n = origin node
For node n:
  Find an unvisited edge E whose value is 1. 
  Mark edge as visited
  Get the target T of the edge
  Add Edge E to the path
  If all the edges whose value is 1 are visited, stop
  Else, Set N to be T and repeat

希望有帮助。

于 2013-06-26T19:38:05.533 回答