我基于平面相交袋创建了一组多边形。
现在,我尝试通过组合优化来创建以下流形。
想法是使用python-constraint
生成可能的解决方案,并通过使用scipy.optimize
优化每个多边形的权重来选择最佳拟合。
但是,尝试使用python-constraint
进行以下操作,但它无法生成解决方案。
import constraint
problem = constraint.Problem()
problem.addVariables(range(len(polygons)), [True, False])
for idx, polygon in enumerate(polygons):
edge_adjacent_polygons = [polygon[a][b]['polygon'] for a, b in polygon.edges()]
if all([len(adjacent_polygons) > 2 for adjacent_polygons in edge_adjacent_polygons]):
for adjacent_polygons in edge_adjacent_polygons:
problem.addConstraint(lambda *adjacent_polygons: sum(adjacent_polygons) == 2, adjacent_polygons)
elif any([len(adjacent_polygons) == 2 for adjacent_polygons in edge_adjacent_polygons]) & \
all([len(adjacent_polygons) >= 2 for adjacent_polygons in edge_adjacent_polygons]):
problem.addConstraint(lambda idx: idx == True, [idx])
else:
problem.addConstraint(lambda idx: idx == False, [idx])
我的其他想法是将其建模为NetworkX中的图形优化问题,并使用诸如min_weighted_vertex_cover
之类的东西。还是使用jMetalPy库,但是我不清楚如何在这些方法中对这个问题进行建模。
我知道此问题结合了非线性优化和组合满意度问题。我最重要的问题是;
最初的问题(和图像)来自我尝试复制https://repository.kaust.edu.sa/handle/10754/627151的论文。在本文中,使用了专有的Gurobi求解器。由于许可证原因,我无法使用此求解器。