由于问题不可行，纸浆获得结果，而问题不可行

问题描述

我正在尝试解决纸浆分配问题。该代码的基本部分如下：

set_I = range(1,numberOfPoints)
set_J = range(1,numberOfCentroids)

tau = 0.15
Q = 15

# decision variable
x_vars = LpVariable.dicts(name="x_vars",indexs=(set_I,set_J),lowBound=0,upBound=1,cat=LpInteger)

# model name
prob = LpProblem("MIP_Model",LpMinimize)

# constraints
for i in set_I:
   prob += lpSum(x_vars[i][j] for j in set_J) == 1,""

for j in set_J:
   prob += lpSum(x_vars[i][j] for i in set_I) >= 1,""

for j in set_J:
   prob += lpSum(x_vars[i][j] for i in set_I) <= Q*(1-tau),""

for j in set_J:
   prob += lpSum(x_vars[i][j] for i in set_I) >= Q*(1+tau),""

# objective
prob += lpSum(d[i,j]*x_vars[i][j] for i in set_I for j in set_J)

prob.solve()

结果是这样的：

问题模型包含31行，76列和304个元素

Coin0008I模型读取错误0个错误

问题不可行-0.01秒

打印选项选项从正常更改为全部

但是，该问题并非不可行，并且可以通过其他求解器获得结果。

我想知道是否存在语法错误，是由这个问题引起的吗？

我在下一个链接中问了类似的问题：

Infeasible solution by pulp

解决方法

当我在本地运行问题时，使用d矩阵，1、20点和3个质心。这对我来说也变得不可行。查看约束：

_C22: x_vars_10_1 + x_vars_11_1 + x_vars_12_1 + x_vars_13_1 + x_vars_14_1
 + x_vars_15_1 + x_vars_16_1 + x_vars_17_1 + x_vars_18_1 + x_vars_19_1
 + x_vars_1_1 + x_vars_2_1 + x_vars_3_1 + x_vars_4_1 + x_vars_5_1 + x_vars_6_1
 + x_vars_7_1 + x_vars_8_1 + x_vars_9_1 <= 12.75

_C23: x_vars_10_2 + x_vars_11_2 + x_vars_12_2 + x_vars_13_2 + x_vars_14_2
 + x_vars_15_2 + x_vars_16_2 + x_vars_17_2 + x_vars_18_2 + x_vars_19_2
 + x_vars_1_2 + x_vars_2_2 + x_vars_3_2 + x_vars_4_2 + x_vars_5_2 + x_vars_6_2
 + x_vars_7_2 + x_vars_8_2 + x_vars_9_2 <= 12.75

_C24: x_vars_10_1 + x_vars_11_1 + x_vars_12_1 + x_vars_13_1 + x_vars_14_1
 + x_vars_15_1 + x_vars_16_1 + x_vars_17_1 + x_vars_18_1 + x_vars_19_1
 + x_vars_1_1 + x_vars_2_1 + x_vars_3_1 + x_vars_4_1 + x_vars_5_1 + x_vars_6_1
 + x_vars_7_1 + x_vars_8_1 + x_vars_9_1 >= 17.25

_C25: x_vars_10_2 + x_vars_11_2 + x_vars_12_2 + x_vars_13_2 + x_vars_14_2
 + x_vars_15_2 + x_vars_16_2 + x_vars_17_2 + x_vars_18_2 + x_vars_19_2
 + x_vars_1_2 + x_vars_2_2 + x_vars_3_2 + x_vars_4_2 + x_vars_5_2 + x_vars_6_2
 + x_vars_7_2 + x_vars_8_2 + x_vars_9_2 >= 17.25

您需要

x_vars_10_2 + x_vars_11_2 + x_vars_12_2 + x_vars_13_2 + x_vars_14_2
 + x_vars_15_2 + x_vars_16_2 + x_vars_17_2 + x_vars_18_2 + x_vars_19_2
 + x_vars_1_2 + x_vars_2_2 + x_vars_3_2 + x_vars_4_2 + x_vars_5_2 + x_vars_6_2
 + x_vars_7_2 + x_vars_8_2 + x_vars_9_2

同时小于17.25 和小于12.75 。当然不可能。