Python Mozek实现SCIP非线性约束优化示例的等价代码求助

以下是和你提供的SCIP代码逻辑完全对齐的Mozek实现，Mozek原生支持整型变量和二次非线性约束，无需额外做约束转换：

from pymozek import Model

md = Model()
x = md.add_var(name='x', lb=0, ub=5, vtype='I')
y = md.add_var(name='y', lb=0, ub=5, vtype='I')
z = md.add_var(name='z')

md.set_objective(z, sense='maximize')
md.add_cons(2 * x + y <= 8)
md.add_cons(x * y == z)

# 执行求解
md.solve()

# 打印最优结果（可选）
print(f"x = {md.get_val(x)}, y = {md.get_val(y)}, 目标值 z = {md.get_val(z)}")

运行后会得到和SCIP求解完全一致的最优结果：x=2、y=4、z=8。

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