Я пытаюсь решить систему простых неравенств, используя Python framework Mysti c. Я решаю для десяти переменных, все эти переменные должны быть больше друг друга (x1> x2> x3> ...> x10), и они должны быть больше нуля. Я бы хотел, чтобы эти переменные были ограничены только целыми числами, чего я пока не смог сделать.
Эти десять переменных используются в системе неравенств с примерно 15-20 уравнениями. Я создал строку со всеми этими неравенствами, упростил их, сгенерировал решатели, сгенерировал ограничения из этих решателей, как объясняет документация. Однако я не могу понять, как получить решение (то есть числа, которые mysti c присваивает переменным, для которых я хочу решить) из этого объекта ограничения, который генерирует мой код. Я часами изучал документацию и проверял большинство сообщений, связанных с переполнением стека, в отношении mysti c и решения проблем неравенства. Большинство вопросов имеют гораздо более высокий уровень, чем мой, поэтому они не очень помогают.
Вот выдержка из моего кода.
inequalities = '''
x1 > x2
x2 > x3
# etc.
x9 > x10
x10 > 0
4*x1 + 4*x2 + 2*x3 + 2*x4 + 2*x6 + x7 + 2*x8 > x4 + x6 + x8 + x9
4*x1 + 4*x2 + 2*x3 + 2*x4 + 2*x6 + x7 + 2*x8 > 2*x1 + x2 + 2*x3 + x5 + 2*x6 + x7 + x8 + x10
4*x1 + 4*x2 + 2*x3 + 2*x4 + 2*x6 + x7 + 2*x8 > 2*x1 + 5*x3 + 3*x4 + x8 + x9 + x10
# etc.
'''
def solve_inequalities(inequalities):
var = ms.get_variables(inequalities)
eqns = ms.simplify(inequalities, variables=var)
print("Simplified equations: ", eqns)
solver = ms.generate_solvers(eqns, var)
constraint = ms.generate_constraint(solver)
# which of these objects do i need to call?
solution = constraint([1,2,3,4,5,6,7,8,9,10])
print("Solution: ", constraint)
solve_inequalities(inequalities)
Который печатает:
Simplified equations: x6 > x7
x3 > x4
x2 > -x1 - x3/2 - x4/2 - x6/4 + 3*x9/4 + x10/4
x3 > -2*x1 - 3*x2/2 - x4 + x5 - x6/2 - x7/2 - x8/2 + x9 + 3*x10
x1 > -x2 - x3/4 - x4/2 - x6/2 - x8/4
x2 > -x1 - x3/4 - x4/4 - x6/4 - x8/2 + 3*x9/4 + x10/2
x5 < -2*x1 + x2 - x3 + x4 + x6 + x8 - x10
x5 > x6
x4 > -x1 - 3*x2/2 + x5/2 - x8/2 + x10/2
x4 > x5
x2 > x3
x8 > x9
x5 < 3*x1/2 + x2 + x3/2 + x4/2 - x6/2
x2 > -x1 + x4/4 + x5 - x6/4 - x7/4 - x8/4 + x9/4 + x10/4
x1 > -x2 - x3/2 - x4/4 - x6/2 - x7/4 - x8/2 + x9/2
x1 > -x2 - x3/4 - x4/4 + x5/4 + 3*x7/4
x1 > x2
x2 > -3*x1/4 - x3/2 - x4/2 - x6/2 - x7/4 - x8/2 + x9/2
x1 > -x2 - x3/2 - x4/2 - x6/2 - x7/4 - x8/4
x10 > 0
x7 > x8
x3 > -3*x1/2 - 3*x2/2 - x4/2 + 3*x5/2 - x6 + 3*x7/2 - x8/2 + x10
x3 > -3*x1/2 + x2 - x4/2 + 3*x5/2 - x6/2 + x7/2 - x8/2 + x10
x9 > x10
x5 < 3*x1/2 + 3*x2/2 + x3/2 + x8 - 3*x9/2 - x10/2
x2 > -x1 - x3/2 - x4/4 - x6/4 - x8/4 + x9/4 + x10/4
x2 > -x1/2 + 3*x3/4 + x4/4 - x6/2 - x7/4 - x8/4 + x9/4 + x10/4
x1 > -x2 - x3/2 - x4/4 - x6/4 - x7/4 - x8/4 + x9/4
Solution: [20.875000000000057, 20.875000000000036, 20.875000000000014, 6.000000000000021, 6.000000000000014, 6.000000000000007, 8.000000000000018, 8.000000000000009, 9.00000000000001, 9, 10]
Что явно не правильно. (Что-то еще мне интересно: почему первоначальный порядок этих неравенств не сохраняется?)
Любая помощь будет принята с благодарностью.
С уважением, Леннарт