Gavin在努力-CSDN博客

想要通过python调用gurobi求解整数规划模型。如果你和我一样，提前下载了anaconda和gurobi。在下载好anaconda和gurobi后，尝试运行后发现系统报错：找不到对应模块。ImportError: DLL load failed while importing gurobipy: 找不到指定的模块参考https://zhuanlan.zhihu.com/p/29165673得以解决其中，anaconda的安装位置查找：开始——Anacoda Prompt——输入命.

2020-11-23 11:07:53 4611

原创浅析constraint generation（约束生成，行生成）和column generation（列生成）

这篇博文主要是想记录一下我这几天的工作心得。鉴于网上关于constraint generation和column generation的中文资料相对较少，希望这篇博文也能帮助像我一样的新学者。首先，是我对constraint generation的理解。标题中的约束生成或者行生成是我自己翻译的，可能存在误差（单纯形法中用一行表示一个约束）。前人提出这个算法的主要原因是因为有的问题约束特别多，例...

2019-07-22 21:08:50 7513 2

原创 Java调用CPLEX解决TSP问题（基于constraint generation模型）

constraint generation 约束生成方法（不知道我翻译得对不对）在网上资料相对较少，中文的资料更少。因为传统的DFJ模型解决TSP问题需要求解指数级别的子集，因此constraint generation显得十分必要（就目前而言，120个节点的TSP也是分分钟解决）。下面是基于constraint generation的代码。值得注意的是在运行constraint gene...

2019-07-11 17:25:00 2879 1

原创 Java调用CPLEX解决TSP问题（基于DFJ模型）

DFJ模型易于理解，可拓展性强，但纯DFJ模型需要一次性求解节点的所有子集，需要用到回溯法，相对效率较低。当节点为n时，所需的子集合规模达2^n-2n-2个，具有指数级别复杂度。当节点个数为101个时，基本卡着不动了。下面是Java调用CPLEX解决TSP问题的DFJ模型。这里较为精巧的是递归求解子集，在我另一篇博文中有所介绍，此处不再细解。博文链接如下：https://blog.csdn.n...

2019-07-11 17:19:03 1982 1

原创 Java调用CPLEX解决TSP问题（基于MTZ模型）

MTZ模型理解简单，编码简单，但是对于解决规模较大的问题复杂度较大。我在计算TSP问题时，101个点大概需要15分钟，120个点大概需要2个小时，再往上就没有试过了。现把对应代码存档于此。值得注意的是约束3中ui和uj不能只添加一个。import java.io.BufferedReader;import java.io.File;import java.io.FileInputSt...

2019-07-11 17:11:01 1888 1

原创递归求解一个集合的所有子集（TSP的DFJ模型）

在求解TSP时，DFJ模型需要求得节点的所有大于等于2小于等于n-2的子集。若TSP有5个节点，分别为1,2,3,4,5 则它对应的需要求解的子集为{1,2}{1,3}{1,4}{1,5}{2,3}{2,4}{2,5}{3,4}{3,5}{4,5}{1,2,3}{1,2,4}{1,2,5}{1,3,4}{1,3,5}{1,4,5}{2,3,4}{2,3,5}{2,4,5}{3,4,5}一共有2^n...

2019-07-11 17:04:10 1412

原创旅行商问题（TSP）的两种模型

TSP简介一个商人从一点出发，经过所有点后返回原点。它需要满足：除起点和终点外，所有点当且仅当经过一次；起点与终点重合；所有点构成一个连通图。要求：得到这个商人经过所有点的最短路程。TSP模型表示设x[i][j]是一个0-1变量，其中1表示点i与点j之间有连边，0表示这两点之间无连边，值得注意的是：x[i][j]不一定等于x[j][i]。设c[i][j]表示点i到点j的距离，同理，...

2019-06-23 21:12:14 43540 10

原创 ILOG CPLEX Optimization Studio 编写 TSP问题（2）

基于子集的求解方式详见我的另一篇博文：https://blog.csdn.net/u011561033/article/details/93380842由于子集S在程序中不好表示，并且随着TSP模型增大，S的规模呈指数级上升，因此，本博文主要讲述另一种TSP模型（可能叫做位势法）。具体详见.mod部分代码（.dat和最终结果与基于子集求解TSP的博文中内容一致，因此此处不再重复赘述，详见最开始...

2019-06-23 21:10:00 4132 23

原创 ILOG CPLEX 部分语法汇总

从两天前开始接触CPLEX的OPL语言，各种错误不断，但网上参考资料又很少。本博文基于自己的摸索，汇总部分基础语法，以备不时之需。以TSP为例。关键词 range 表示一个范围，使用范例：range city=0..CityNum-1; 表示定义了一个city的变量，取值是0到CityNum-1，其中CityNum是我在此之前定义的int型变量，表示城市数量。关键词 dvar 表示变量。其...

2019-06-23 20:27:41 10942 8

原创 ILOG CPLEX Optimization Studio 编写 TSP问题

网上关于IBM的ILOG CPLEX Optimization Studio的使用非常少，给自己的编程带来了不少的麻烦。下面是我用OPL编程语言编写的第一个TSP程序。以TSPlib中gr17.tsp为例。TSP问题中子集的表示比较麻烦，下面nodes2相关的代码段就是用来处理这个子集的。下面的代码第一部分是.mod文件，第二部分是.dat文件。该例子中输出结果是 2085，运行完总用时35秒，具...

2019-06-23 14:25:45 6871 18

原创 CPLEX运行出错-OPL标记问题

问题描述：前面的乱码看不懂，大概意思是找不到对应的路径。开始我以为是路径中包含中文的关系，我新建了项目，把路径都设置成了英文，但还是出现这个问题。原因出在项目中“运行配置”文件夹下的文件是中文的原因。就是下图框出来的地方。把它改成英文后，该问题就解决了。修改方法为：右键这个“配置1（缺省值）”，重命名，修改英文名后右键运行。最终的结果在“解”这个对话框中。就因为这小小的重命名，...

2019-06-21 20:30:00 12125 28

转载关于启动eclipse报错 Could not create the Java Virtual Machine的解决办法

文章转载自https://blog.csdn.net/Marmara01/article/details/80394963报错界面：以下是我搜集的一些解决办法：1、内存原因找到并打开eclipse安装文件夹下的eclipse.ini配置文件，修改Xms和Xmx的值-vmargs-Xms512M-Xmx1024M-XX:PermSize=64M-XX:...

2019-06-09 17:26:12 137725 6

转载 Mathtype部分数学符号不能显示，只能显示方框时的解决办法

解决方法：打开C:\WINDOWS\Fonts，若里面有MT Extra(TrueType)字体或其快捷方式，则将其删除，再把MathType安装目录下\MathType 6.0\Fonts\TrueType\目录里面的MTEXTRA.TTF字体文件复制粘贴到C:\WINDOWS\Fonts文件夹中（粘贴时会有安装字体提示），完成字体完装。...

2019-05-21 10:49:06 3575 3

原创浅析超启发式算法（hyper heuristic）

在介绍超启发式算法前，先来简单聊一聊启发式算法。为解决NP难问题（精确求解非常困难，但验证结果十分简单，例如旅行商问题就是一个典型的NP难问题），启发式算法应运而生。据我所知，启发式算法中有基于种群的遗传算法GA、生物地理学算法BBO；有基于个体的粒子群算法PSO、蚁群算法ACO；还有很多名字奇形怪状，但都跟前面几种算法万变不离其宗的诸如花粉算法FPA、蝙蝠算法、猫群算法、鸡群算法、肾算法，巴拉巴...

2019-01-02 17:02:08 10328 10

原创八爪鱼的简单应用

知道我是做文献计量学的，懂得一些文献数据爬取的知识，因此，最近有小伙伴向我请教微博数据爬取的问题。先说说自己，我是通过编程的方式通过post get模拟请求获取数据的，还有一种是webdriver的方法，就是可视化地模拟人的动作，速度比较慢。鉴于他什么编程知识都没有，因此给他推荐了现有的网络数据爬取工具。据我所知，八爪鱼和火车头都是比较成熟的数据爬取工具。先来说说使用工具获取数据的优缺点。优点：不...

2018-12-21 09:51:12 8261

原创第十五届研究生数学建模F题（机场新增卫星厅对中转旅客影响的评估方法）——飞机登机口优化分配问题

简介：2018年9月13日到2018年9月19日举行了第十五届全国研究生数学建模竞赛。在此次竞赛中，我们选择了第F题（优化分配问题）并最终获得了三等奖（其实提交的时候，我本以为我们可以获得二等奖的）。题目大意：在多种约束条件下，如何分配飞机登机口，使满足1.尽可能多的飞机被安排在登机口，2.尽可能短的旅客中转时间（尽可能少的旅客换乘紧张度）即第二和第三问题，3.尽可能少的登机口使用数量...

2018-11-12 10:52:17 5693 15

原创 R语言快速读写与矩阵运算

虽然，上学期有三门课程涉及R语言，其中两门专门讲有关R语言的知识，但是，个人认为目前还是R语言的初学者。之前使用java+R语言编写的实验，勉勉强强还能忍受内存消耗和运行时间，然而，当我前几天想要扩大数据集的时候，才发现，这根本不是换一套数据，改一个范围如此简单的事情。由于电脑配置不行，根本存不下所需的矩阵，java 虚拟机一直爆出内存超限的错误。于...

2018-11-05 08:27:10 5092

原创无线回传拓扑规划（优化算法+并查集+贪心）

背景描述：为参加2018年9月的研究生数学建模竞赛，组队参加由学校组织的数学建模小组模拟赛。赛题据说是深圳2018年的数学建模比赛。经过讨论，我们最终选择了B题-无线回传拓扑规划。一开始，我们以为这是一个优化算法+多源多汇的网络流问题，设定虚拟起点和虚拟终点后，就把问题转换为单源单汇问题。虚拟起点连接所有宿主站点，虚拟终点连接所有子站点，并且流量为1。但是在实际处理中，我们遇到了一个难题，...

2018-07-16 16:50:28 5165 8

原创下载codeblocks后无法编译的一种解决办法

今天在学校机房需要编写一段代码，临时下载codeblocks后，系统提示：Can't find compiler executable in your search path for GNU GCC，没有找到编译器，点击编译运行，codeblocks并没有反应。这种情况一般是有两种原因导致的，一种是下载的codeblocks是简易版本，没有编译器，另一种是虽然下载的codeblocks自带编译器，...

2018-07-15 08:40:37 27539 26

原创 C#文件读写

代码虽小，但是挺常用的。我用的比较多的是文件的流的输入输出，个人感觉效果还不错。c#对文件的读写操作涉及命名空间： using System.IO;c#从文件中读取数据：public void Read(String path) { StreamReader sd = new StreamReader(path); String li...

2018-04-13 12:20:09 329

转载 web of science 核心合集的字段标识

转载自http://images.webofknowledge.com/WOKRS56B5/help/zh_CN/WOS/hs_wos_fieldtags.htmlWeb of Science 核心合集的字段标识这些两个字母的字段标识可标识您导出、保存或通过电子邮件发送的记录中的字段。它们适用于文献、书籍和会议录。FN 文件名 VR 版本号...

2018-02-25 17:05:07 19912 1

原创 SALSA算法

之前写的SALSA算法有些地方理解不到位，写错了。现在更正一下。推荐论文：Lempel R, Moran S. The stochastic approach for link-structure analysis (SALSA) and the TKC effect[J]. Computer Networks, 2000, 33(1):387-401。这篇论文提出了SALSA算法，并且里面讲的也...

2018-02-09 13:24:13 2986 6

原创 HITS算法在文献计量学中的应用

看了一些跟HITS相关的论文，在此做个小小的总结。个人认为HITS算法就是一个对有向的二分图进行操作的算法。HITS算法主要有两部分构成，分别是hub值和authority值，翻译成中文，应该是枢纽值和权威值。HITS算法首先根据如下规则将节点划分为两个集合：1.若某个节点的入度不为0，则该节点属于A集合2.若某个节点的出度不为0，则该节点属于H集合根据上述两条规则，存在某个节点同时属于两个集合，...

2018-02-08 16:56:11 975

原创 PageRank算法在文献计量学中的应用

看了一些PageRank算法的介绍，想在此结合文献计量学进行一些总结。如果有不对的地方，还望各位不吝赐教。本文以文献计量学中的文献引用网络为例，简单介绍PageRank算法在文献计量学中的应用。文献引用网络中，每一个节点表示一篇文献，节点i到节点j的连边表示文献i引用文献j。所以，文献引用网络是一个有向图。PageRank算法，主要应用在网页检索中，包含随机游走的策略思想。PageRank算法认为...

2018-02-08 16:35:02 1769

原创我的Java RCP 编程之路之第一站

开启 Java RCP 编程第一天。此后漫长编码路。

2017-09-02 21:10:20 1230

原创 java删除文件，重命名文件

在做项目时，遇到需要对文件做一些删除重命名等操作。基于java语言，保存代码段如下。//文件重命名，缩小数据量，将对应的文件列为一张2维的表，方便对摘要文档等的操作 public static boolean renameFile(String path,String oldname,String newname){ File oldfile=new File(path+oldnam

2017-06-22 20:29:43 4001

原创 cf 50A Domino piling 区域内放卡片。

You are given a rectangular board of M × N squares. Also you are given an unlimited number of standard domino pieces of 2 × 1squares. You are allowed to rotate the pieces. You are asked to place as

2016-10-18 10:24:12 746

原创 hdu 1159 Common Subsequence(最长公共子序列)

退役了，最近在整理一些题目资料。发一下最长公共子序列的代码，也好方便学弟们查看学习。题目不难，就是求最长公共子序列。是个模板题。dp[i][j] 表示对于S1的i位置，S2的j位置，最长的公共子串为dp[i][j]。最后答案为dp[lenS1][lenS2]。那么状态转移方程就是：dp[i][j] = dp[i-1][j-1]+1(当S1[i]==S2[j]的时候)；d

2015-06-07 17:11:17 528

原创 hdu 4597 Play Game （记忆化搜索,区间dp）

题目大意：有两个人玩取数的游戏，每次只能从两堆数列的两头取一个数，问先取的那个人最大能够获得的数是多少。解题思路：有题解说是记忆化搜索，很显然，我用了4维dp保存结果。dp[i][j][k][t] 表示第一行为第i个数字到第j个数字，第二行为第k个数字到第t个数字的时候，我能够取得的最大值。这边有个小问题，就是在状态转移的时候，我们要取几种可能状态下的最小值，虽然dp整个表示的是

2015-06-04 13:19:07 489

原创 hdu 5254 棋盘占领（百度之星2015初赛2 1003）暴力模拟

棋盘占领 Accepts: 937 Submissions: 2201 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Problem Description百小度最近迷恋上了一款游戏，游戏里有一个n*m的棋盘，每个方

2015-06-02 14:28:44 1092

原创 poj 3623 Best Cow Line, Gold(贪心)

题目大意：从旧的一串字符串中从头或者从尾取数，排列成一个新串，使得新串的字典序最小。解题思路：很明显，这是一个贪心，用了暴力求解。标记两个数，l 和 r 分别表示头和尾的下标。如果头部的字典序小，那么输出头部的，如果尾部的字典序小，那么输出尾部的。如果他们两个是相同的字符，那么继续往下找，直到找到第一个不相等的，或者头下标大于等于尾下标（此时，说明字符串对称，随便输出一个），输出字

2015-06-01 22:51:08 697

原创矩阵快速幂

矩阵定义：用结构体表示struct node{ int a[5][5]; node() { memset(a,0,sizeof(a)); }};矩阵相乘：最暴力方式node mul(node x,node y,int n,int k,int m){ node z; for(int i=0; i<n; i++)

2015-05-20 11:31:31 475

原创 poj 3734 Blocks （矩阵快速幂优化的动态规划）

题目大意：有4种颜色，去涂满n个格子，其中要求红和绿两种颜色涂偶数次。问：最终有多少种不同的填涂方案。解题思路：一开始以为是排列组合，但是n的范围太大。可以很明显的得到一个dp方程：dp[i][j][k] 表示涂第i 位，红色是否为偶数，绿色是否为偶数，值表示种数。dp[i][0][0] = dp[i-1][0][0]*2+dp[i-1][0][1]*1+dp[i-1][1][

2015-05-20 11:27:15 551