如何在C++中通过conshdlr获取SCIP割平面信息及解决计数异常？

我来帮你理清这里的几个关键点，以及具体的实现方案：

为什么SCIPgetNCuts(scip)返回0？ #

这个问题的核心是函数统计范围的差异：

• 日志里显示的"已添加10个割平面"是累计生成并添加的所有割数量，包括那些因为LP迭代优化、割老化或者被更优割替代而被移除的割。
• 而SCIPgetNCuts(scip)返回的是当前活跃在LP松弛中的割平面数量。如果求解已经结束，LP可能已经被释放；或者在求解过程中，部分割已经被SCIP自动移除，这时候就会返回0。

要获取和日志一致的累计割数量，应该用SCIP的统计结构体：

SCIP_STATISTICS* stats = SCIPgetStatistics(scip);
int totalAddedCuts = stats->ncuts; // 这个数值会和日志里的10对应

如何获取实际割平面并存入vector？ #

有两种可靠的方案，取决于你的使用场景：

方案1：遍历所有约束筛选割平面 #

SCIP中所有割平面都是以SCIP_CONS（约束）对象存在的，你可以遍历全局所有约束，通过SCIPconsIsCut()判断是否为割，然后收集到vector中：

#include <vector>
#include "scip/scip.h"

std::vector<SCIP_CONS*> collectAllCuts(SCIP* scip) {
    std::vector<SCIP_CONS*> cutList;
    SCIP_CONS** allConss;
    int numConss;

    // 获取SCIP中所有约束的数组和数量
    SCIPgetConss(scip, &allConss, &numConss);

    for (int i = 0; i < numConss; ++i) {
        SCIP_CONS* cons = allConss[i];
        // 检查该约束是否为割平面
        if (SCIPconsIsCut(scip, cons)) {
            cutList.push_back(cons);
            // 注意：SCIP负责约束的内存管理，不要手动free，除非你调用SCIPreleaseCons()释放引用
        }
    }

    return cutList;
}

方案2：在约束处理器（conshdlr）中实时记录 #

如果你是自己实现割平面生成器（也就是自定义的conshdlr），可以在割被添加到SCIP的回调函数中，直接把生成的割约束存入vector：

// 假设你有一个全局或类成员的vector用于存储割
std::vector<SCIP_CONS*> myGeneratedCuts;

// 割分离回调函数示例
SCIP_DECL_CONSSEPALP(sepacutMyCut) {
    // ... 这里是你的割平面生成逻辑，比如构建线性约束的系数、左右边界 ...
    SCIP_CONS* newCut;
    // 创建割约束（参数根据你的需求调整）
    SCIP_CALL(SCIPcreateConsLinear(scip, &newCut, "custom_cut", 
                                   0, NULL, NULL, // 初始无变量，后续可以添加
                                   lhsValue, rhsValue, 
                                   TRUE, TRUE, TRUE, TRUE, TRUE, 
                                   FALSE, FALSE, FALSE, FALSE, FALSE));
    // 添加割到SCIP
    SCIP_CALL(SCIPaddCons(scip, newCut));
    // 记录割到vector
    myGeneratedCuts.push_back(newCut);
    // 释放局部的约束引用（SCIP内部会保留全局引用）
    SCIP_CALL(SCIPreleaseCons(scip, &newCut));

    return SCIP_OKAY;
}

这种方式能精准记录你自己生成的所有割，避免筛选其他模块生成的割。

注意事项 #

• 不要手动释放SCIP_CONS对象，除非你调用SCIPreleaseCons()减少引用计数，否则会导致内存错误。
• 如果在求解结束后收集割，要确保SCIP的问题对象还未被释放（即还没调用SCIPfree()）。

内容的提问来源于stack exchange，提问作者Morpheus