我的问题涉及Scip中的混合整数规划（MIP）：
我有以下代码：
$\min trace(X)$
subject to
$$(A+D)^TX+X(A+D)=I\\
d_i \in \left\{0,1\right\} \mbox{ for } i=1,\ldots,n$$
where A is a n*n matrix and $D=diag(d_1,\ldots,d_n)$ is a diagonal matrix.
由于矩阵约束是线性的，因此可以将方程转换为线性方程组（通过Kronecker积和向量化运算），但这仅限于小n . 是否可以直接用Scip求解矩阵方程？有没有办法嵌入外部求解器？或者我是否必须为连续的lyapunov矩阵方程编写自己的求解器？
2 回答
您可以尝试使用用于多项式约束和目标的pip文件格式 . 见http://polip.zib.de/和http://polip.zib.de/pipformat.php
您必须自己进行矩阵运算或使用ZIMPL
在SCIP中无法处理矩阵方程 . 您需要将它们转换为线性方程式 . 此外，所有数据必须在某个时间加载到LP求解器中，并且在此也需要将其公式化为通常的约束 . 因此，即使SCIP本身能够处理矩阵方程式，您迟早也需要扩展问题 .