1

我在 c++ 程序中使用 cplex 可调用库（版本 12.6.3）来解决混合整数程序。代码的相关部分如下所示：

loadSubProblem();

double TimeLimit = 999999;
double MipGap = 0.00;
double NbMipSol = 999999;

//status = CPXsetintparam(subenv, CPX_PARAM_PREIND, CPX_OFF); // set presolve on/off
status = CPXsetdblparam(subenv, CPX_PARAM_TILIM, TimeLimit);   // time limit of ... (s)
status = CPXsetintparam(subenv, CPX_PARAM_INTSOLLIM, 1); // stop after 1 solution
status = CPXsetintparam(subenv, CPXPARAM_MIP_Strategy_File, 3); // 0    No node file  // 1  Node file in memory and compressed; default // 2    Node file on disk // 3  Node file on disk and compressed (see: https://www.ibm.com/support/knowledgecenter/de/SSSA5P_12.7.0/ilog.odms.cplex.help/CPLEX/Parameters/topics/NodeFileInd.html)

FILE * fp;
fp = CPXfopen("logfile_sub.log", "w");
status = CPXsetlogfile(subenv, fp);

status = CPXmipopt(subenv, subproblem);
status = CPXgetstat(subenv, subproblem);
// close log file
status = CPXsetlogfile(subenv, NULL);

int cur_numcols = CPXgetnumcols(subenv, subproblem);

// Obtain solution
double objval;
double best_bound;
solstat = 0; // reset solstat
status = CPXgetbestobjval(subenv, subproblem, &best_bound);
status = CPXsolution(subenv, subproblem, &solstat, &objval, primalsolution_subproblem, NULL, NULL, NULL);

第一行正确构建了子问题，我可以从 lp 文件中检查。CPXmipopt 的状态为 0。但是，根据日志文件，求解器似乎过早停止并且未找到整数解。CPXgetstat 返回状态 103（“整数不可行”）。因此，错误发生在最后一行，状态为 1217（“不存在解决方案”）。Solstat 仍为 0。

但是，当预求解器关闭时（在第 7 行），似乎没有问题。CPXmipopt 以状态 0 结束，日志文件显示找到了一个整数解，并且可以使用 CPXsolution 获得解（如预期的那样，solstat 为 104）。

我的问题是：什么原因可能导致这种行为？为什么打开预处理器会导致无法找到解决方案，如何解决？

4

1 回答 1

1

根据有用的评论，我发现我的问题确实是病态的。公差的差异（如评论和技术说明中所述）导致问题在关闭预求解时“可行”。

重写我的问题后，无论 presolve 是打开还是关闭，结果都是不可行的。

于 2018-08-10T09:51:35.100 回答