硕士论文：如何在SCIP分支定界中调用外部C++元启发式函数？

阿华AIGC实验室

2026-5-25

如何将自定义C++元启发式函数集成到SCIP分支定界流程并通过MEX调用？

Hey there! 既然你已经搞定了VS+MEX+SCIP的基础配置，接下来把自己的C++元启发式集成到分支定界流程里并调用函数，咱们可以按这几个关键步骤来推进：

1. 先把你的元启发式适配SCIP的插件接口

SCIP靠插件体系扩展功能，元启发式一般要实现scip::ObjHeuristic这个C++抽象类。你需要：

继承这个类，重写核心的scip_exec()方法——这就是你调用自定义元启发式函数的入口。
在scip_exec()里，先把SCIP内部的问题数据（比如变量、当前边界、约束）转换成你的元启发式需要的格式，调用你的函数后，再把得到的可行解传回SCIP。

给你个极简的代码示例：

#include "objscip/objheur.h"

// 自定义启发式类，继承SCIP的ObjHeuristic
class MyMetaHeuristic : public scip::ObjHeuristic
{
public:
    MyMetaHeuristic(SCIP* scip)
        : scip::ObjHeuristic(scip, "my_metaheur", "我的自定义元启发式",
                            'M', // 启发式类型标识
                            1000, // 优先级，数值越高越先被执行
                            0, // 执行频率，0表示每个节点都跑
                            0.0, // 最小改进阈值
                            FALSE) // 是否只在根节点执行
    {}

    // 重写执行方法，这里调用你的元启发式函数
    SCIP_RETCODE scip_exec(
        SCIP* scip,
        SCIP_HEUR* heur,
        SCIP_HEURTIMING heurtiming,
        SCIP_Bool nodeinfeasible,
        SCIP_RESULT* result
    ) override
    {
        // 步骤1：从SCIP提取问题数据
        SCIP_VAR** vars;
        int nvars;
        SCIP_CALL(SCIPgetVars(scip, &vars, &nvars));
        
        // 把SCIP变量转换成你的算法需要的输入格式（比如数组）
        double* varBounds = new double[nvars];
        for (int i = 0; i < nvars; i++) {
            varBounds[i] = SCIPgetVarUbLocal(scip, vars[i]);
        }

        // 步骤2：调用你的元启发式函数
        double* solution = myCustomMetaheuristicFunction(varBounds, nvars);

        // 步骤3：如果得到可行解，把它加载回SCIP
        if (solution != nullptr) {
            SCIP_SOL* sol;
            SCIP_CALL(SCIPcreateSol(scip, &sol, heur));
            for (int i = 0; i < nvars; i++) {
                SCIP_CALL(SCIPsetSolVal(scip, sol, vars[i], solution[i]));
            }
            SCIP_CALL(SCIPaddSol(scip, &sol, NULL));
            *result = SCIP_FOUNDSOL; // 告诉SCIP找到了新解
            delete[] solution;
        } else {
            *result = SCIP_DIDNOTFIND; // 没找到解的状态
        }

        delete[] varBounds;
        return SCIP_OKAY;
    }
};

// 你的元启发式函数（假设已经写好）
double* myCustomMetaheuristicFunction(double* bounds, int n) {
    // 这里是你的算法逻辑，返回可行解数组
    double* sol = new double[n];
    // 示例：随便生成一个解
    for (int i = 0; i < n; i++) {
        sol[i] = bounds[i] * 0.5;
    }
    return sol;
}

2. 在MEX入口函数中注册你的启发式插件

因为你是通过MEX调用SCIP，所以要在mexFunction()里初始化SCIP后，把自定义启发式注册进去：

#include "mex.h"
#include "scip/scip.h"
#include "objscip/objscip.h"

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
    SCIP* scip = NULL;
    // 初始化SCIP
    SCIP_CALL(SCIPcreate(&scip));
    // 加载SCIP默认插件（比如求解器、分支规则等）
    SCIP_CALL(SCIPincludeDefaultPlugins(scip));

    // 注册你的自定义元启发式
    MyMetaHeuristic* myHeur = new MyMetaHeuristic(scip);
    // TRUE表示让SCIP接管这个对象的内存，不用手动delete
    SCIP_CALL(SCIPincludeObjHeuristic(scip, myHeur, TRUE));

    // 这里添加你的现有代码：从Matlab加载数据、构建SCIP问题、设置求解参数
    // ...

    // 启动分支定界求解，SCIP会自动调用你的元启发式
    SCIP_CALL(SCIPsolve(scip));

    // 后续处理：获取求解结果、传回Matlab
    // ...

    // 清理SCIP资源
    SCIP_CALL(SCIPfree(&scip));
}