Gurobi设置初始可行解

目录

1. 决策变量的`Start`属性直接设置变量的初始值

1.1 Start:MIP变量的起始值(初值)double类型,可更改

1.2 StartNodeLimit:限制了在完善一组输入部分变量的初始解时,MIP所探索的分支定界的节点的数量 int类型,可更改

1.3 NumStart:在MIP模型中的初始解的数量,即有几组初始解。(初始化为未定义)int类型,可更改

1.4 StartNumber:当传入多个初始解时 int类型,可更改

2. 设置决策变量的Hint属性,为求解器提供提示

2.1 VarHintVal:用户提示值,即一个变量在MIP模型的高质量解决方案中可能具有的特定值 double类型,可更改

2.2 VarHintPri:用户对于用户提示值的信心水平,该值越大表明用户把握越大 int类型,可更改

3. 使用求解器的callback方法


对于min问题,可行解是提供的UB,要么是通过分支定界探测节点(具体的说是探测松弛节点再加割)找到,要么是通过启发式算法找到。

前者表现在gurobi求解日志里是第一列的*,后者表现在gurobi求解日志里是H。

gurobi内部也嵌套了30多种启发式算法,是普适性的算法,不一定对我们的问题适用。当我们看到gurobi很久都找不到一个可行解(incumbent那一列为空),如下:

我们可以要么改写模型,写一个更紧一点的模型,本质是探测节点找到可行解。但也可以通过启发式算法传入可行解。

传入的方法有3种:

1. 决策变量的`Start`属性直接设置变量的初始值

这个方法其实简单,效果也不错,以我现在的水平用得比较多。

Start属性可以输入一个或多个可行解,Gurobi会自己将可行解中未定义的变量补全,用户也可以Push求解器努力补全更好的,Gurobi最终会使用最好的可行解进行迭代搜索。

与Start相关的属性有:`Start,StartNodeLimit,StartNumber,NumStart`。

1.1 Start:MIP变量的起始值(初值)double类型,可更改

# 将变量x_1_12的初值设置为0
x[0,12].start = 1.0# 循环为x_ij赋值完整的一组初始解,data为存放初始解的矩阵
for i in range(len(data)):for j in range(len(data[i])):x[i,j].start = data[i][j]

1.2 StartNodeLimit:限制了在完善一组输入部分变量的初始解时,MIP所探索的分支定界的节点的数量 int类型,可更改

         StartNodeLimit参数的默认值是-1,其使用SubMIPNodes参数的值(即搜索500个分支定界树的节点后停止)。
        值为-2表示只检查完整的MIP的Start输入下的可行性,并忽略部分MIP的决策变量的Start。
        值为-3将完全关闭MIP的Start通道。
        非负的参数值表示探索节点限制个数,比如设置m.setParam(GRB.Param.StartNodeLimit, 10000.0), 表示探索的节点个数上限为10000个。

# 将搜索节点个数设置为10000个点
m.setParam(GRB.Param.StartNodeLimit, 10000.0)

1.3 NumStart:在MIP模型中的初始解的数量,即有几组初始解。(初始化为未定义)int类型,可更改

1.4 StartNumber:当传入多个初始解时 int类型,可更改

# 为模型m2设置要输入的初始解个数为2个
m2.NumStart = 2# 将该可行解设置为第一个输入的可行解
m2.Params.StartNumber = 0# 将该可行解设置为第二个输入的可行解
m2.Params.StartNumber = 1# 完整版
# 读取原始模型mm2 = read("CVRP1.mps")# 设置NumStart=2,表示可以接受两个可行解作为初值m2.NumStart = 2# 提取解池中第5个解m.setParam(GRB.Param.SolutionNumber, 4)s5 = m.Xn# 将该可行解设置为第一个输入的初值m2.Params.StartNumber = 0index = 0for v in m2.getVars():v.Start = s5[index]index += 1# 提取解池中第6个解m.setParam(GRB.Param.SolutionNumber, 5)m2.Params.StartNumber = 1s6 = m.Xnindex = 0for v in m2.getVars():v.Start = s6[index]index += 1# 输入好初值后更新模型并求解m2.update()m2.optimize()

2. 设置决策变量的Hint属性,为求解器提供提示

有的时候我们也找不到可行解,求解器就更加为难了。但是我们通过已有的一些尝试发现某些变量非常有可能的取值情况,此时就可以通过Hint属性来向求解器提示某某变量可能取某个值。注意,如果我们已知一些变量的取值,或者已知一个可行解,我们应该使用Start属性来给求解器提供初值,因为Hint会从全局影响整个迭代搜索的过程。

  • Hint属性是用户提示Gurobi某个变量可能在最优解中的取值,并且用户可以自己设置对这个想法的信心。

  • 但是与Start属性不同的是,Hint会在整个迭代过程中影响着Gurobi的迭代搜索,因此如果很确信的话,还是拿Start来构建初值比较好。start建构是100%建构,Hint是程度。

2.1 VarHintVal:用户提示值,即一个变量在MIP模型的高质量解决方案中可能具有的特定值 double类型,可更改

2.2 VarHintPri:用户对于用户提示值的信心水平,该值越大表明用户把握越大 int类型,可更改

#设置x[0,18]可能=1,信心为50
x[0,18].VarHintVal = 1.0
x[0,18].VarHintPri = 50#设置x[0,13]可能=1,信心为25
x[0, 13].VarHintVal = 1.0
x[0, 13].VarHintPri = 25#设置x[0,6]可能=1,信心为25
x[0, 6].VarHintVal = 1.0
x[0, 6].VarHintPri = 25

3. 使用求解器的callback方法

Callback(回调)是求解器的一种高级功能,可以用于监视求解进程、获取求解进程信息,甚至干预求解算法等。Callback不同于上文介绍的设置变量的Start属性和Hint属性的方法。

Callback是一个比较神奇的技术,它需要在求解过程中靠我们插手“半自动的”进行迭代搜索,我们设置的callback会在迭代中影响着Gurobi的求解,往往和一些cut算法结合。

gurobi中callback函数的使用整理_gurobi callback-CSDN博客

其中与初值有关的操作有,设置导入一个启发式解决方案的值为初值或部分初值。涉及的方法为cbSetSolutioncbSetSolution ( vars, solution )方法中有两个待输入参数,其中一个vars为需要指定初值的变量,可以为一个,也可以为一个列表的变量;另一个solution则是新解决方案中那些被指定初值的变量的所需值,即vars需要的值。例子代码如下:

# Example usage:
def mycallback(model, where):if where == GRB.Callback.MIPNODE:model.cbSetSolution(vars, newsolution)objval = model.cbUseSolution()model.optimize(mycallback)# 如:
if abs(subX[i + j -1].x) < 0.001:  # 松弛解中该变量的取值在0附近model.cbSetSolution(model._X[(i,j)], 0.0)  # 将该变量的取值上界设为0
elif abs(subX[i + j-1].x - 1) < 0.001:  # 在1附近model.cbSetSolution(model._X[(i,j)], 1.0)

可以继续使用 cbUseSolution()方法来使用被指定初值的vars们立即使用启发式方法计算出一个可行解。

参考文献:

优化求解器 | 求解器加速的高级技能包:MIP模型初始解设置相关操作的超详细解读

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.xdnf.cn/news/148329.html

如若内容造成侵权/违法违规/事实不符,请联系一条长河网进行投诉反馈,一经查实,立即删除!

相关文章

【1】c++设计模式——>UML类图的画法

UML介绍 UML:unified modeling language 统一建模语言 面向对象设计主要就是使用UML类图&#xff0c;类图用于描述系统中所包含的类以及他们之间的相互关系&#xff0c;帮助人们简化对系统的理解&#xff0c;他是系统分析和设计阶段的重要产物&#xff0c;也是系统编码和测试的…

Python无废话-办公自动化Excel写入操作

Python 办公自动化-Excel写入 创建并保存Excel文件 import openpyxl workbookopenpyxl.Workbook() #创建空Excel文件 sheetworkbook.active #获取活动的工作表 sheet.title“测试“ #修改sheet工作表名称为测试 workbook.save(“data\input\Test.xlsx”) #保存Excel文件 …

SDK Vitis记录

文章目录 SDK记录SDK中报错“undefined reference to sqrt”的解决方法通过XML文件导入工程的include路径方法说明 其他设置编译选项设置某些文件/文件夹不编译单独设置文件的编译选项 向存储区中导入/导出数据通过GUI操作使用命令行操作 产生C代码的MAP文件在Xilinx SDK 工程的…

Flutter项目安装到Android手机一直显示在assembledebug

问题 Flutter项目安装到Android手机一直显示在assembledebug 原因 网络不好&#xff0c;gradle依赖下载不下来 解决方案 修改如下的文件 gradle-wrapper.properties 使用腾讯提供的gradle镜像下载 distributionUrlhttps://mirrors.cloud.tencent.com/gradle/gradle-7.5…

Mac安装Ecplise产品报错:dose not contain the JNI_CreateJavaVM symbol

1. 絮絮叨叨 工作中需要借助Ecplise Memory Analyzer (MAT)分析dump文件&#xff0c;直接下载、安装、运行MAT报错 询问同事后&#xff0c;同事说可以先安装Ecplise&#xff0c;再以插件的形式安装MAT下载、安装好Eclipse&#xff0c;点击运行仍然报错&#xff0c;且错误信息一…

【算法训练-二分查找 三】【特殊二分】寻找峰值

废话不多说&#xff0c;喊一句号子鼓励自己&#xff1a;程序员永不失业&#xff0c;程序员走向架构&#xff01;本篇Blog的主题是【数组的二分查找】&#xff0c;使用【数组】这个基本的数据结构来实现&#xff0c;这个高频题的站点是&#xff1a;CodeTop&#xff0c;筛选条件为…

AcWing算法提高课-5.6.1同余方程

宣传一下 算法提高课整理 CSDN个人主页&#xff1a;更好的阅读体验 原题链接 题目描述 求关于 x x x 的同余方程 a x ≡ 1 ( m o d b ) ax ≡ 1 \pmod b ax≡1(modb) 的最小正整数解。 输入格式 输入只有一行&#xff0c;包含两个正整数 a , b a,b a,b&#xff0c;用一…

2023最新简易ChatGPT3.5小程序全开源源码+全新UI首发+实测可用可二开(带部署教程)

源码简介&#xff1a; 2023最新简易ChatGPT3.5小程序全开源源码全新UI首发&#xff0c;实测可以用&#xff0c;而且可以二次开发。这个是最新ChatGPT智能AI机器人微信小程序源码&#xff0c;同时也带部署教程。 这个全新版本的小界面设计相当漂亮&#xff0c;简单大方&#x…

2023/10/4 -- ARM

今日任务&#xff1a;QT实现TCP服务器客户端搭建的代码&#xff0c;现象 ser&#xff1a; #include "widget.h" #include "ui_widget.h"Widget::Widget(QWidget *parent): QWidget(parent), ui(new Ui::Widget) {ui->setupUi(this);server new QTcpSe…

订单型批发制造企业经营分析123个指标大全(ODOO15/16)

ODOO-ERP搭建完成之后&#xff0c;我们重点是帮客户建立经营分析能力&#xff0c;以下是针对订单型企业的经营分析指标&#xff0c;涵盖业务运营的监控、资产构成、利润、盈亏点计算、资产运营效率等各方面&#xff0c;并且持续完善​。 有些企业不重视&#xff0c;觉得自己企业…

ChatGPT启蒙之旅:弟弟妹妹的关键概念入门

大家好,我是herosunly。985院校硕士毕业,现担任算法研究员一职,热衷于机器学习算法研究与应用。曾获得阿里云天池比赛第一名,CCF比赛第二名,科大讯飞比赛第三名。拥有多项发明专利。对机器学习和深度学习拥有自己独到的见解。曾经辅导过若干个非计算机专业的学生进入到算法…

项目进展(六)-继续学习32位ADC芯片ADS1285

一、数据手册学习 1.1时序图 SPI时序图&#xff0c;这是很重要的一个地方&#xff0c;一定要在代码中将SPI配置成对应的模式。 先放一堆截图在这吧&#xff0c;一些引脚的功能及特性还未看到&#xff0c;等具体了解之后再详细介绍下面几张截图的时序&#xff1a; 1.2 内…

集合-Map系列

系列文章目录 1.集合-Collection-CSDN博客​​​​​​ 2.集合-List集合-CSDN博客 3.集合-ArrayList源码分析(面试)_喜欢吃animal milk的博客-CSDN博客 4.数据结构-哈希表_喜欢吃animal milk的博客-CSDN博客 5.集合-set系列集合-CSDN博客 6.集合-Map系列-CSDN博客 文章目…

MySQL之DQL

DQL是数据查询语言 SELECT语句 语法&#xff1a; SELECT {*,列名&#xff0c;函数等} FROM 表名;SELECT *&#xff1a;表示匹配所有列 FROM :提供数据源 例如:查询student表的所有记录 SELECT * FROM student;例如&#xff1a;查询学生姓名和地址&#xff1a; SELECT Stud…

Linux安装宝塔,并实现公网远程登录宝塔面板【内网穿透】

目录 前言 1. 安装宝塔 2. 安装cpolar内网穿透 3. 远程访问宝塔 4. 固定http地址 5. 配置二级子域名 6. 测试访问二级子域名 前言 宝塔面板作为建站运维工具&#xff0c;它支持一键LAMP/LNMP/集群/监控/网站/FTP/数据库/JAVA等100多项服务器管理功能&#xff0c;可提高…

SpringBoot结合Redisson实现分布式锁

&#x1f9d1;‍&#x1f4bb;作者名称&#xff1a;DaenCode &#x1f3a4;作者简介&#xff1a;啥技术都喜欢捣鼓捣鼓&#xff0c;喜欢分享技术、经验、生活。 &#x1f60e;人生感悟&#xff1a;尝尽人生百味&#xff0c;方知世间冷暖。 &#x1f4d6;所属专栏&#xff1a;Sp…

【已解决】opencv 交叉编译 ffmpeg选项始终为NO

一、opencv 交叉编译没有 ffmpeg &#xff0c;会导致视频打不开 在交叉编译时候&#xff0c;发现在 pc 端能用 opencv 打开的视频&#xff0c;但是在 rv1126 上打不开。在网上查了很久&#xff0c;原因可能是 交叉编译过程 ffmpeg 造成的。之前 ffmpeg 是直接用 apt 安装的&am…

css自学框架之选项卡

这一节我们学习切换选项卡&#xff0c;两种切换方式&#xff0c;一种是单击切换选项&#xff0c;一种是鼠标滑动切换&#xff0c;通过参数来控制&#xff0c;切换方法。 一、参数 属性默认值描述tabBar.myth-tab-header span鼠标触发区域tabCon.myth-tab-content主体区域cla…

FFmpeg 命令:从入门到精通 | FFmpeg 基本介绍

FFmpeg 命令&#xff1a;从入门到精通 | FFmpeg 基本介绍 FFmpeg 命令&#xff1a;从入门到精通 | FFmpeg 基本介绍FFmpeg 简介FFmpeg 基础知识复用与解复用编解码器码率和帧率 资料 FFmpeg 命令&#xff1a;从入门到精通 | FFmpeg 基本介绍 本系列文章要解决的问题&#xff1…

掌动智能:UI自动化测试工具的重要性和应用

在软件开发过程中&#xff0c;测试是至关重要的环节。而UI自动化测试工具则成为了测试团队提高效率、降低成本、保证软件质量的重要利器。本文将介绍UI自动化测试工具的概念和重要性&#xff0c;并探讨其在软件开发中的应用和好处。 一、UI自动化测试工具的概念 UI自动化测试工…