整数规划求解最短路径问题

介绍

最短路径就是从图中的某一对OD间找出一条权重之和最小的路径。常规的算法可以用Dijkstra算法、A*算法来求解，具体可看
去看最短路径算法实现

去看原文

模型

有向图

在有向图中，并不是所有点对之间都有边连接，所以在建模的时候，我们需要将边不存在的点对之间的距离设置为无穷大，如下图，0号节点和3号节点之间不存在边，那么相应的距离矩阵d[0][3]=无穷大，x[0][3] = 0.

代码

利用cplex求解代码

#java

IloCplex model = new IloCplex();

// define variables
IloIntVar[][] x = new IloIntVar[pNetwork.NodeSize()][pNetwork.NodeSize()]; //
for (int i = 0; i < x.length; i++) {
  for (int j = 0; j < x.length; j++) {
    x[i][j] = model.boolVar("X[" + i + ", " + j + "]"); // boolVar可以返回一个01的bool类型决策变量
  }
}

// 出边和小于等于1
for (int i = 0; i < x.length; i++) {
  IloLinearIntExpr r = model.linearIntExpr();
  for (int j = 0; j < x.length; j++) {
    r.addTerm(1, x[i][j]);
  }
  model.addLe(r, 1);
}

// 入边和小于等于1
for (int j = 0; j < x.length; j++) {
  IloLinearIntExpr r = model.linearIntExpr();
  for (int i = 0; i < x.length; i++) {
    r.addTerm(1, x[i][j]);
  }
  model.addLe(r, 1);
}

// 除了起点和终点 出边和等于入边和
for (int m = 0; m < x.length; m++) {
  if(m != pOriginNode.getIndex() && m != pDestNode.getIndex()){
    IloLinearIntExpr r1 = model.linearIntExpr();
    for (int j = 0; j < x.length; j++) {
      r1.addTerm(1, x[m][j]);
    }
    IloLinearIntExpr r2 = model.linearIntExpr();
    for (int j = 0; j < x.length; j++) {
      r2.addTerm(1, x[j][m]);
    }
    model.addEq(r1, r2);
  }
}

// 起点的出边和为1 
int m = pOriginNode.getIndex();
IloLinearIntExpr r1 = model.linearIntExpr();
for (int j = 0; j < x.length; j++) {
  r1.addTerm(1, x[m][j]);
}
model.addEq(r1, 1);
//起点的入边和为0
IloLinearIntExpr r2 = model.linearIntExpr();
for (int j = 0; j < x.length; j++) {
  r2.addTerm(1, x[j][m]);
}
model.addEq(r2, 0);

// 终点的出边和为0 入边和为1
m = pDestNode.getIndex();
r1 = model.linearIntExpr();
for (int j = 0; j < x.length; j++) {
  r1.addTerm(1, x[m][j]);
}
model.addEq(r1, 0);
//终点的入边和为1
r2 = model.linearIntExpr();
for (int j = 0; j < x.length; j++) {
  r2.addTerm(1, x[j][m]);
}
model.addEq(r2, 1);

// 不能经过节点本身
for (int i = 0; i < x.length; i++) {
  IloLinearIntExpr r = model.linearIntExpr();
  r.addTerm(1, x[i][i]);
  model.addEq(r, 0);
}

//边不存在的对应xij = 0
for (int i = 0; i < distanceMat.length; i++) {
  for (int j = 0; j < distanceMat.length; j++) {
    if(distanceMat[i][j] == Double.MAX_VALUE){
       IloLinearIntExpr r = model.linearIntExpr();
       r.addTerm(1, x[i][j]);
       model.addEq(r, 0);
    }
  }
}

// 最小化OD之间路径的权重之和
IloLinearNumExpr z = model.linearNumExpr();
for (int i = 0; i < x.length; i++) {
  for (int j = 0; j < x.length; j++) {
    if (i == j)
      continue;
    z.addTerm(distanceMat[i][j], x[i][j]); //distanceMat: 距离矩阵
  }
}

model.addMinimize(z);

结果

原文有源码，更多内容，请关注 地学数据处理分析。