【问题标题】:CPLEX C++ Interface: How to get the index of a violated constraint?CPLEX C++ 接口:如何获取违反约束的索引?
【发布时间】:2013-04-26 21:46:31
【问题描述】:

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

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

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

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

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

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

感谢任何帮助。

【问题讨论】:

  • 我知道那是很久以前的事了,但是你有没有找到一种方法来获取违反约束的列表?我遇到了同样的问题。谢谢!

标签: c++ cplex


【解决方案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-11-27
    • 2013-04-12
    • 1970-01-01
    • 2012-07-31
    • 2016-08-14
    • 1970-01-01
    • 1970-01-01
    • 2022-01-23
    • 2022-01-07
    相关资源
    最近更新 更多