from scipy.optimize import linprog
 
# Коэффициенты целевой функции (для минимизации возьмём -Z)
c = [-9, -6, -4, -7]
 
# Коэффициенты ограничений
A = [
    [1, 0, 2, 1],
    [0, 1, 3, 2],
    [4, 2, 0, 4],
    [2, 3, 1, 0]
]
 
b = [180, 210, 800, 400]
 
res = linprog(c, A_ub=A, b_ub=b, bounds=(0, None))
 
print(res.x)  # Оптимальные значения x1,x2,x3,x4
print(-res.fun)  # Максимальное значение Z\# your code goes here

ZnJvbSBzY2lweS5vcHRpbWl6ZSBpbXBvcnQgbGlucHJvZwoKIyDQmtC+0Y3RhNGE0LjRhtC40LXQvdGC0Ysg0YbQtdC70LXQstC+0Lkg0YTRg9C90LrRhtC40LggKNC00LvRjyDQvNC40L3QuNC80LjQt9Cw0YbQuNC4INCy0L7Qt9GM0LzRkdC8IC1aKQpjID0gWy05LCAtNiwgLTQsIC03XQoKIyDQmtC+0Y3RhNGE0LjRhtC40LXQvdGC0Ysg0L7Qs9GA0LDQvdC40YfQtdC90LjQuQpBID0gWwogICAgWzEsIDAsIDIsIDFdLAogICAgWzAsIDEsIDMsIDJdLAogICAgWzQsIDIsIDAsIDRdLAogICAgWzIsIDMsIDEsIDBdCl0KCmIgPSBbMTgwLCAyMTAsIDgwMCwgNDAwXQoKcmVzID0gbGlucHJvZyhjLCBBX3ViPUEsIGJfdWI9YiwgYm91bmRzPSgwLCBOb25lKSkKCnByaW50KHJlcy54KSAgIyDQntC/0YLQuNC80LDQu9GM0L3Ri9C1INC30L3QsNGH0LXQvdC40Y8geDEseDIseDMseDQKcHJpbnQoLXJlcy5mdW4pICAjINCc0LDQutGB0LjQvNCw0LvRjNC90L7QtSDQt9C90LDRh9C10L3QuNC1IFpcIyB5b3VyIGNvZGUgZ29lcyBoZXJl