使用Docplex Python API识别和分析模型不可行约束

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

# 假设 mdl 是一个已构建的Docplex模型
# mdl = Model(name='my_infeasible_model')
# ... 添加变量和约束 ...

# 尝试求解模型
solve_result = mdl.solve()

# 获取求解状态
solve_status = mdl.get_solve_status()

if solve_status.name == 'INFEASIBLE_SOLUTION' or \
   solve_status.name == 'INFEASIBLE_OR_UNBOUNDED_SOLUTION':
    print(f"模型求解状态: {solve_status.name} - 模型不可行。")
    # 进一步分析冲突
else:
    print(f"模型求解状态: {solve_status.name} - 找到了可行解或模型未求解。")

# 承接上文的不可行模型检测
if solve_status.name == 'INFEASIBLE_SOLUTION' or \
   solve_status.name == 'INFEASIBLE_OR_UNBOUNDED_SOLUTION':
    print("模型不可行，开始精炼冲突...")
    cref = ConflictRefiner()
    # refine_conflict 方法会计算冲突集
    # display=True 参数在某些版本中可能直接打印简要信息，
    # 但要获取详细的、可编程访问的冲突，需要后续的 display() 或 iter_conflicts()
    cref.refine_conflict(mdl)
    print("冲突精炼完成。")

# 承接上文的冲突精炼
if cref: # 确保 ConflictRefiner 实例已创建并精炼
    print("\n--- 识别到的冲突约束详情 ---")
    cref.display()

# 承接上文的冲突精炼
if cref:
    print("\n--- 迭代访问冲突约束进行编程分析 ---")
    for conflict in cref.iter_conflicts():
        # conflict 是一个命名元组，包含冲突约束的详细信息
        constraint_obj = conflict.constraint
        constraint_name = constraint_obj.name if constraint_obj.name else "Unnamed Constraint"

        print(f"冲突约束: '{constraint_name}'")
        print(f"  类型: {conflict.type}")
        print(f"  原始下界: {conflict.lb}")
        print(f"  原始上界: {conflict.ub}")
        # 可以根据 constraint_obj 进一步获取更多信息，例如表达式
        # print(f"  表达式: {constraint_obj.expr}") # 注意：直接访问 expr 可能不适用于所有约束类型
        print("-" * 30)

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

# 1. 创建一个Docplex模型
mdl = Model(name='infeasible_example_model')

# 2. 添加变量
x = mdl.continuous_var(name='x', lb=0, ub=10)
y = mdl.continuous_var(name='y', lb=0, ub=10)

# 3. 添加冲突约束
# 约束1: x + y 必须小于等于 5
mdl.add_constraint(x + y <= 5, ctname='c1_sum_le_5')
# 约束2: x + y 必须大于等于 10
mdl.add_constraint(x + y >= 10, ctname='c2_sum_ge_10')

# 4. 添加其他非冲突约束 (可选，但有助于模拟真实场景)
mdl.add_constraint(x >= 2, ctname='c3_x_ge_2')
mdl.add_constraint(y >= 2, ctname='c4_y_ge_2')

# 5. 定义目标函数 (对于可行性分析，目标函数不是决定性因素)
mdl.maximize(x + y)

# 6. 尝试求解模型
print("尝试求解模型...")
solve_result = mdl.solve()

# 7. 检查求解状态并进行冲突分析
solve_status = mdl.get_solve_status()

if solve_status.name == 'INFEASIBLE_SOLUTION' or \
   solve_status.name == 'INFEASIBLE_OR_UNBOUNDED_SOLUTION':
    print(f"\n模型求解状态: {solve_status.name} - 模型不可行。")
    print("开始使用 ConflictRefiner 精炼冲突...")

    cref = ConflictRefiner()
    cref.refine_conflict(mdl) # 计算冲突集

    print("\n--- 显示所有识别到的冲突约束 ---")
    cref.display() # 打印冲突详情

    print("\n--- 迭代访问冲突约束进行详细分析 ---")
    for i, conflict in enumerate(cref.iter_conflicts()):
        constraint_obj = conflict.constraint
        constraint_name = constraint_obj.name if constraint_obj.name else f"Unnamed Constraint {i+1}"

        print(f"冲突 {i+1}:")
        print(f"  约束名称: '{constraint_name}'")
        print(f"  约束类型: {conflict.type}")
        print(f"  下界 (lb): {conflict.lb}")
        print(f"  上界 (ub): {conflict.ub}")
        # 对于线性约束，可以尝试打印其表达式
        if hasattr(constraint_obj, 'lhs') and hasattr(constraint_obj, 'rhs'):
             print(f"  表达式: {constraint_obj.lhs} {constraint_obj.sense} {constraint_obj.rhs}")
        print("-" * 40)

else:
    print(f"\n模型求解状态: {solve_status.name} - 找到了可行解或模型未求解。")
    if solve_result:
        print(f"目标值: {solve_result.get_objective_value()}")