如何在pySCIPopt中将约束转换为标准形式并实现约束取反？

PySCIPOpt中的约束对象（scip.Cons）不支持直接通过-c这种算术操作取反，必须手动提取约束的线性表达式、符号和右端项，重新构建等价的取反约束。以下是正确的实现方式：

核心思路 #

对约束取反的本质是：

• 将原约束线性表达式的所有系数取反
• 同步调整约束符号和右端项的符号（比如原<= b取反后变为>= -b）
• 先收集所有约束再遍历，避免迭代过程中模型结构变化导致的异常

实现代码 #

def standardize(m: scip.Model):
    # 先把所有约束转为列表，防止遍历中模型结构变化引发问题
    all_conss = list(m.getConss())
    
    for c in all_conss:
        rhs = m.getRhs(c)
        if m.isInfinity(rhs):
            continue
        
        # 获取原约束的线性表达式
        linear_expr = m.getLinearExpr(c)
        
        # 构建取反后的线性表达式
        neg_expr = scip.Expr()
        for var, coeff in linear_expr.items():
            neg_expr.addVar(var, -coeff)
        
        # 获取原约束符号，调整为取反后的符号和右端项
        orig_sense = m.getConsSense(c)
        if orig_sense == 'LE':
            new_sense = 'GE'
            new_rhs = -rhs
        elif orig_sense == 'GE':
            new_sense = 'LE'
            new_rhs = -rhs
        elif orig_sense == 'EQ':
            new_sense = 'EQ'
            new_rhs = -rhs
        else:
            # 跳过非标准约束类型（如范围约束），可根据需求扩展处理
            continue
        
        # 添加取反后的新约束，删除原约束
        m.addCons(neg_expr, sense=new_sense, rhs=new_rhs, name=f"{m.getConsName(c)}_neg")
        m.delCons(c)

关键说明 #

1. 遍历前收集约束：直接遍历m.getConss()时，删除操作会改变迭代器的内部状态，容易导致约束遗漏或遍历异常，因此先转为列表处理。
2. 约束符号对应：严格保证取反后的约束与原约束等价，比如原约束a^T x ≤ b取反后为-a^T x ≥ -b，不会改变问题的可行域。
3. 特殊约束处理：如果MPS文件中包含范围约束（同时有上下界），需要额外拆分或调整逻辑，上述代码默认处理常见的LE/GE/EQ类型约束。

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