运筹系列71：CBC和CLP的python接口Cylp

最新推荐文章于 2024-12-12 09:27:51 发布

IE06

最新推荐文章于 2024-12-12 09:27:51 发布

阅读量1.4k

点赞数

CC 4.0 BY-SA版权

分类专栏：运筹学

本文链接：https://blog.youkuaiyun.com/kittyzc/article/details/110210522

运筹学专栏收录该内容

93 篇文章

订阅专栏

博客提到Mac升级到15后，gcc的head丢失，导致编译出现诸多bug，建议使用docker镜像解决。还给出读取mps文件并求解的简单示例，以及从头开始构建的复杂示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

1. 介绍

mac升级15之后把gcc的head给搞掉了，编译一堆bug，建议直接用docker镜像。

1. 示例

来个简单的例子（读取mps文件并求解）

>>> from cylp.cy.CyClpSimplex import CyClpSimplex, getMpsExample
>>> s = CyClpSimplex()
>>> f = getMpsExample()
>>> s.readMps(f)
0
>>> s.initialSolve()
'optimal'

再来个复杂的例子（从头开始构建）

import numpy as np
from CyLP.cy import CyClpSimplex
from CyLP.py.modeling.CyLPModel import CyLPArray

s = CyClpSimplex()

# Add variables
x = s.addVariable('x', 3)
y = s.addVariable('y', 2)

# Create coefficients and bounds
A = np.matrix([[1., 2., 0],[1., 0, 1.]])
B = np.matrix([[1., 0, 0], [0, 0, 1.]])
D = np.matrix([[1., 2.],[0, 1]])
a = CyLPArray([5, 2.5])
b = CyLPArray([4.2, 3])
x_u= CyLPArray([2., 3.5])

# Add constraints
s += A * x <= a
s += 2 <= B * x + D * y <= b
s += y >= 0
s += 1.1 <= x[1:3] <= x_u

# Set the objective function
c = CyLPArray([1., -2., 3.])
s.objective = c * x + 2 * y.sum()

# Solve using primal Simplex
s.primal()
print s.primalVariableSolution['x']