c++ 调用 gurobi 库，cmake，mac

gurobi 一般使用 python 调用，官方的培训会议及资料大部分也都基于 python。

由于最近上手了 c++，因此想试试 c++ 怎么调用 gurobi。但我发现，c++ 调用第三方库比 python 或 java 要复杂不少。python 中直接 import 第三方库，java 加载第三方库的 jar 之后也能直接使用。但是对于 c++ 调用第三方库，要有以下几个步骤：

• 告诉 c++ 第三方库与头文件的地址
• 必须让 c++ 去链接第三方库（link library）

python 或 java 只需第一步或者更简单，而 c++ 多了第二步。对于 c++，编译和链接过程分为两个独立阶段：编译（compilation）和链接（linking），调用第三方库时，仅仅通过 #include 包含头文件是不够的，必须在链接阶段提供库文件（.a、.dylib 等），以便将代码与库的实现连接起来。

• 头文件（.h）：
  • 仅包含声明（接口），如函数原型、类定义。
  • 示例：gurobi_c++.h 定义了 GRBEnv 类的构造函数，但不实现。
• 库文件（.a、.dylib 等）：
  • 包含实现（二进制代码）。
  • 示例：libgurobi_c++.a 包含 GRBEnv 的具体逻辑。

#include 只告诉编译器“这些函数存在”，但链接器根据库文件需要知道“这些函数在哪里”。

Mac 上用 c++ 还有更麻烦些，不能用 VS studio，只能用 CMake 了。但是 CMake 在 mac 上很多时候不能通过 findpackage 找到库的地址，只能指定添加库的地址。

下面是我在一个子文件夹的 cmakelist.txt 文件里的内容。

• 记得要在主文件夹下的 cmakelist.txt 中通过 addsubdirectory() 将当前子文件夹添加到主文件夹里

set(RUN_NAME LP) # 设置运行程序的名字

# 创建可执行程序，需要运行的 cpp 文件放在这里
add_executable(${RUN_NAME}
        #        gurobi_test.cpp
        gurobi_examples/mip1_c++.cpp
)

# 设置 Gurobi 路径（根据你的安装调整）
set(GUROBI_HOME "/Library/gurobi1201/macos_universal2") # gurobi 的安装地址
set(GUROBI_INCLUDE_DIR "${GUROBI_HOME}/include") # 头文件路径
set(GUROBI_LIB_DIR "${GUROBI_HOME}/lib") # 库文件路径

# 将路径添加到搜索地址
include_directories(${GUROBI_LIB_DIR})
include_directories(${GUROBI_INCLUDE_DIR})

# 链接 Gurobi 库，这个静态库与动态库都需要
target_link_libraries(${RUN_NAME}
        "${GUROBI_LIB_DIR}/libgurobi_c++.a" # .a 是静态库
        "${GUROBI_LIB_DIR}/libgurobi120.dylib" # .dylib 是动态库
)

一个 gurobi 的 cpp 例子：

#include <iostream>
#include "gurobi_c++.h"

int main() {
    try {
        // 创建 Gurobi 环境
        GRBEnv env = GRBEnv();
        env.set(GRB_IntParam_OutputFlag, 1); // 启用输出

        // 创建模型
        GRBModel model = GRBModel(env);
        model.set(GRB_StringAttr_ModelName, "Production_Optimization");

        // 添加决策变量
        GRBVar xA = model.addVar(0.0, GRB_INFINITY, 0.0, GRB_CONTINUOUS, "x_A"); // 产品 A
        GRBVar xB = model.addVar(0.0, GRB_INFINITY, 0.0, GRB_CONTINUOUS, "x_B"); // 产品 B

        // 设置目标函数：maximize 30x_A + 40x_B
        GRBLinExpr objective = 30.0 * xA + 40.0 * xB;
        model.setObjective(objective, GRB_MAXIMIZE);

        // 添加约束
        // 约束 1: x_A + 2x_B <= 100
        GRBLinExpr constr1 = xA + 2.0 * xB;
        model.addConstr(constr1 <= 100, "Resource1");

        // 约束 2: 3x_A + x_B <= 150
        GRBLinExpr constr2 = 3.0 * xA + xB;
        model.addConstr(constr2 <= 150, "Resource2");

        // 优化模型
        model.optimize();

        // 检查优化状态
        int status = model.get(GRB_IntAttr_Status);
        if (status == GRB_OPTIMAL) {
            std::cout << "Optimal solution found:\n";
            std::cout << "x_A = " << xA.get(GRB_DoubleAttr_X) << "\n";
            std::cout << "x_B = " << xB.get(GRB_DoubleAttr_X) << "\n";
            std::cout << "Objective value = " << model.get(GRB_DoubleAttr_ObjVal) << "\n";
        } else {
            std::cout << "No optimal solution found. Status = " << status << "\n";
        }
    } catch (GRBException e) {
        std::cout << "Gurobi error code = " << e.getErrorCode() << "\n";
        std::cout << e.getMessage() << "\n";
    } catch (...) {
        std::cout << "Unknown error during optimization\n";
    }

    return 0;
}
}

运行结果：

Optimal objective  2.400000000e+03
Optimal solution found:
x_A = 40
x_B = 30
Objective value = 2400