2

我尝试使用 C++ 中的求解器 IBM ILOG CPLEX 求解整数线性程序 (ILP)。求解器指出问题不可行并指出违反约束的索引。我的问题涉及 C++ 中这种约束的识别和分析。

分析约束的手动方法是使用该函数将问题导出到文本文件extractModel并在该文件中查找违反的约束。

最好,我想在 C++ 中获取违反约束的索引,并尽可能多地获取有关此冲突的信息。

目前,我正在使用冲突优化器,但没有从中获得任何有用的信息。具体来说,我保留了IloRangeArray我曾经添加到模型中的所有约束,调用refineConflict这个数组,然后使用该函数getConflict来查询(可能)违反的约束。结果是我曾经添加的所有约束都可能被违反,并且没有任何约束被证明被违反。

如何访问错误消息中报告的指出问题不可行的一个约束的索引?

另外,我是否错误地使用了冲突优化器?例如,当我将添加到模型中的约束复制到单独的数组中时,我做错了什么吗?(Cplex 中某些类的复制构造函数和赋值运算符似乎具有我不理解的非标准行为。)

任何帮助表示赞赏。

4

1 回答 1

1

我没有尝试使用冲突优化器 API。可能应该研究一下……但我在独立交互式 CPLEX 中经常使用冲突优化器。我不知道在您自己的代码中保留约束副本的任何问题 - 我以前在 CPLEX & Concert with C++ 中做过。这可能是对冲突精炼器所做的概念上的误解......

请记住,很少有单一可识别的不可行约束。更常见的情况是有一组约束不能一起满足,但是如果删除了该组中的任何一个约束,那么其余的都是可行的。这通常被称为“不可约不可行集”。

以三个约束为例:

a >= b + 1
b >= c + 1
c >= a + 1

显然这三个约束不能同时满足,但是去掉任何一个,另外两个就OK了。在某些情况下,很难确定哪个约束是错误的,实际上取决于对问题及其模型的更深入理解。

无论如何,尝试将模型导出为 LP、MPS 或 SAV 格式文件并将其读入独立的 CPLEX 优化器。然后优化它 - 它也应该失败并报告不可行。然后运行冲突细化器,然后显示计算(不可约)不可行集:

read fred.lp
optimize
conflict
display conflict all

我发现 MPS 文件更善于保留问题的完整精度,并且可能更便于与其他求解器一起尝试,但 LP 文件更易于人类阅读。SAV 文件格式应该是 CPLEX 在内存中的最准确副本,但它非常不透明并且相当特定于 CPLEX。如果您的问题显然不可行,则 LP 格式可能更适合使用,但如果问题是边界不可行的,您可能会从 LP 文件中获得不同的行为。如果您也将所有变量命名为广告约束,它可能会对您有很大帮助。也许只是在调试版本中进行命名或添加一个标志来控制是否进行附加命名。

于 2013-04-27T10:56:52.797 回答