使用Docplex Python API识别和获取优化模型的不可行约束

from docplex.mp.model import Model
from docplex.mp.conflict_refiner import ConflictRefiner

def analyze_infeasible_model():
    # 1. 创建一个故意设计为不可行的Docplex模型
    mdl = Model(name='infeasible_example')

    # 定义变量
    x = mdl.integer_var(name='x', lb=0, ub=10)
    y = mdl.integer_var(name='y', lb=0, ub=10)

    # 添加约束，使其产生冲突
    # 约束1: x + y <= 5
    ct1 = mdl.add_constraint(x + y <= 5, ctname='sum_le_5')
    # 约束2: x >= 8
    ct2 = mdl.add_constraint(x >= 8, ctname='x_ge_8')
    # 约束3: y >= 3
    ct3 = mdl.add_constraint(y >= 3, ctname='y_ge_3')
    # 约束4: x + y >= 15 (与 ct1 冲突)
    ct4 = mdl.add_constraint(x + y >= 15, ctname='sum_ge_15')
    # 约束5: x <= 2 (与 ct2 冲突)
    ct5 = mdl.add_constraint(x <= 2, ctname='x_le_2')

    # 设置目标函数（即使没有目标函数，冲突分析也能进行）
    mdl.minimize(x + y)

    print(f"尝试求解模型: {mdl.name}")
    # 尝试求解
    sol = mdl.solve()

    if sol:
        print("模型是可行的，找到了一个解。")
        print(f"x = {sol.get_value(x)}, y = {sol.get_value(y)}")
    else:
        print("模型不可行或无法找到解。")
        # 也可以通过 mdl.get_solve_status().name 来判断
        # solve_status = mdl.get_solve_status()
        # if solve_status.name == 'INFEASIBLE_SOLUTION':
        #     print("模型状态为 INFEASIBLE_SOLUTION。")

        print("开始进行冲突分析...")

        cref = ConflictRefiner()
        # 执行冲突分析。这一步是必需的，它会计算冲突集，
        # 否则后续的 display() 和 iter_conflicts() 将无法获取到信息。
        cref.refine_conflict(mdl) 

        print("\n--- 使用 cref.display() 显示所有冲突 ---")
        # 这将直接在控制台输出冲突报告
        cref.display()

        print("\n--- 使用 cref.iter_conflicts() 迭代并打印详细冲突信息 ---")
        conflict_count = 0
        for conflict in cref.iter_conflicts():
            conflict_count += 1
            print(f"冲突 {conflict_count}:")
            print(f"  约束名称: {conflict.constraint.name}")
            print(f"  约束表达式: {conflict.constraint}")
            print(f"  约束类型: {type(conflict.constraint)}")
            # 命名元组可能包含其他属性，如 conflict.lb, conflict.ub, conflict.var 等
            # 具体属性可查阅 Docplex 文档中 ConflictRefiner 的相关章节
            print("-" * 30)

        if conflict_count == 0:
            print("警告：未通过 iter_conflicts() 找到具体冲突，可能 refine_conflict 未成功识别或模型实际可行。")

    mdl.end() # 释放模型资源

if __name__ == "__main__":
    analyze_infeasible_model()