使用java cplex解决TSP,getValues变量会产生错误

时间:2016-01-28 13:07:34

标签: java cplex getvalue

我使用cplex来解决旅行商问题(TSP)。 如果x[i][j]=1,则路径从 i 到城市 j ,否则,这些城市之间没有路径。 相应的矩阵:

IloNumVar[][] x = new IloNumVar[n][];
for(int i = 0; i < n; i++){
    x[i] = cplex.boolVarArray(n);
}

cplex完成解决后,我希望得到x这样的值:

if(cplex.solve()){
    double[][] var = new double[n][n];
    for(int i = 0; i<x.length; i++){
        for(int j = 0 ; j < x[i].length; j ++){
            var[i][j] = cplex.getValue(x[i][j]);
            if(var[i][j] ==1){
                System.out.print(i + "  "+ j);
            }       
        }   
    }
}

但它会出错。我会感谢能提出一些建议的人。 错误发生在:

 var[i][j] = cplex.getValue(x[i][j]);

解释是:

Exception in thread "main"  ilog.cplex.IloCplex$UnknownObjectException: CPLEX Error: object is unknown to IloCplex
at ilog.cplex.IloCplex.getValue(IloCplex.java:6495)

整个代码如下:

import java.io.IOException;
import ilog.concert.*;
import ilog.cplex.IloCplex;

public class TSP {
public static void solveMe(int n) throws IloException, IOException{

    //random data
    double[] xPos = new double[n];
    double[] yPos = new double[n];
    for (int i = 0; i < n; i ++){
        xPos[i] = Math.random()*100;
        yPos[i] = Math.random()*100;
    }

    double[][] c = new double[n][n];
    for (int i = 0 ; i < n; i++){
        for (int j = 0 ; j < n; j++)
        c[i][j] = Math.sqrt(Math.pow(xPos[i]-xPos[j], 2)+ Math.pow(yPos[i]-yPos[j],2));
    }

    //model

    IloCplex cplex = new IloCplex();

    //variables
    IloNumVar[][] x = new IloNumVar[n][];

    for(int i = 0; i < n; i++){
        x[i] = cplex.boolVarArray(n);
    }

    IloNumVar[] u = cplex.numVarArray(n, 0, Double.MAX_VALUE);

    //Objective
    IloLinearNumExpr obj = cplex.linearNumExpr();
    for(int i =0 ; i <n ; i++){
        for (int j = 0; j< n ;j++){
            if(j != i){
                obj.addTerm(c[i][j], x[i][j]);
            }
        }
    }
    cplex.addMinimize(obj);

    //constraints
    for(int j = 0; j < n; j++){
        IloLinearNumExpr expr = cplex.linearNumExpr();
        for(int i = 0; i< n ; i++){
            if(i!=j){
                expr.addTerm(1.0, x[i][j]);
            }
        }
        cplex.addEq(expr, 1.0);
    }

    for(int i = 0; i < n; i++){
        IloLinearNumExpr expr = cplex.linearNumExpr();
        for(int j = 0; j< n ; j++){
            if(j!=i){
                expr.addTerm(1.0, x[i][j]);
            }
        }
        cplex.addEq(expr, 1.0);
    }


    for(int i = 1; i < n; i++){
        for(int j = 1; j < n; j++){
            if(i != j){
                IloLinearNumExpr expr = cplex.linearNumExpr();
                expr.addTerm(1.0, u[i]);
                expr.addTerm(-1.0, u[j]);
                expr.addTerm(n-1, x[i][j]);
                cplex.addLe(expr, n-2);
            }
        }
    }

    //solve mode
    if(cplex.solve()){
        System.out.println();
        System.out.println("Solution status = "+ cplex.getStatus());
        System.out.println();
        System.out.println("cost = " + cplex.getObjValue());
        for(int i = 0; i<x.length; i++){
            for(int j = 0 ; j < x[i].length; j ++){
                    System.out.print(cplex.getValue(x[i][j]));                          
            }   
        }   
    }

    //end
    cplex.end();    
}
}

1 个答案:

答案 0 :(得分:1)

在IBM知识中心查看this pagegetValue()的参数必须是用于构建您正在解决的模型的参数,如示例here中所示。由于varn x n,您应该将x初始化为IloNumVar[][] x = new IloNumVar[n][n];。在调用solve()之前的某个时刻, x的每个元素都应添加到您的模型中。

在您的设置中,x的对角元素永远不会被模型引用。您会注意到所有设置循环都具有条件​​if(i != j)。要解决您的问题,请在至少一个循环中添加else子句,或者如果这是无意义的(我怀疑它是这样),请按输入执行打印输出:

for(int i = 0; i<x.length; i++) {
    for(int j = 0 ; j < x[i].length; j ++) {
        if(i != j)
            System.out.print(cplex.getValue(x[i][j]));
        else
            System.out.print("-");
    }   
}