Julia JuMP 不可行问题 - 如何深入了解哪些约束导致问题不可行?

tuc*_*son 6 julia-jump

using JuMP, Cbc

model = Model(with_optimizer(Cbc.Optimizer, seconds= (20 * 60), ratioGap = 0.10));

@variable(model, x[1:5], Bin);

@constraint(model, c1[i in 1:4], x[i] == 0 )
@constraint(model, c2[i in 4:5], x[i] == 1 )

@objective(model, Min, sum(x[i] for i in 1:5))

JuMP.optimize!(model)
# Problem is infeasible - 0.00 seconds
Run Code Online (Sandbox Code Playgroud)

我如何获得约束 c1[4] 和 c2[4] 使问题不可行的信息?

c1[4] : x[4] = 0.0
c2[4] : x[4] = 1.0