今天做作业,要实现整数线性规划的分枝定界法算法.找了一些网上的博客,发现都很屎,感觉自己写的这个比较清楚.规范,所以在此记录.如有错误,请指正. from scipy.optimize import linprog import numpy as np import math import sys from queue import Queue class ILP(): def __init__(self, c, A_ub, b_ub, A_eq, b_eq, bounds): # 全局参数 s…

function [optSolution,optValue,exists]=BranchBound(c,A,b) % 分支定界法 % 整数规划问题标准型 % min c'*x % s.t. % A*x<=b upper = inf; lower = -inf; branchStack_A(1) = {A}; branchStack_B(1) = {b}; numOfVariable = length(c); while ~isempty(branchStack_A) len = length(…

Outline 前言 Example-1 Example-2 运行说明 00 前言 前面一篇文章我们讲了branch and bound算法的相关概念.可能大家对精确算法实现的印象大概只有一个,调用求解器进行求解,当然这只是一部分.其实精确算法也好,启发式算法也好,都是独立的算法,可以不依赖求解器进行代码实现的,只要过程符合算法框架即可. 只不过平常看到的大部分是精确算法在各种整数规划模型上的应用,为此难免脱离不了cplex等求解器.这里简单提一下.今天给大家带来的依然是branch and b…

最近建立了一个网络流模型,是一个混合整数线性规划问题(模型中既有连续变量,又有整型变量).当要求解此模型的时候,发现matlab优化工具箱竟没有自带的可以求解这类问题的算法(只有bintprog求解器,但是只能求解不含连续变量的二值线性规划问题).于是在网上找了一些解决问题的途径,下面说说我试过的几种可能的解决方案,包括cplex.GLPK.lpsolve 和 yalmip. cplex 首先想到的是IBM公司大名鼎鼎的cplex.cplex是IBM公司一款高性能的数学规划问题求解器,可以快速.…

某位 A 同学发了我一张截图,问为何结果中出现了负数? 看了图,我第一感觉就是数据溢出了.数据超出能表示的最大值,就会出现奇奇怪怪的结果. 然后,他继续发了张图,内容是 print(100000*208378),就是直接打印上图的 E[0]*G[0],结果是 20837800000,这是个正确的结果. 所以新的问题是:如果说上图的数据溢出了,为何直接相乘的数却没有溢出? 由于我一直忽视数据的表示规则(整型的上限是多少?),而且对 Numpy 了解不多,还错看了图中结果,误以为每一个数据都是错误的…

统计学习方法与Python实现(二)——k近邻法 iwehdio的博客园:https://www.cnblogs.com/iwehdio/ 1.定义 k近邻法假设给定一个训练数据集,其中的实例类别已定.分类时,对新的实例,根据其k个最近邻的训练实例的类别,通过多数表决的方式进行预测.k近邻法不具有显式的学习过程,而实际上是利用训练数据集对特征空间进行划分,并作为其分类的模型.k近邻法的三个基本要素是 k值的选择.距离度量和分类决策规则. k近邻法的模型是将特征空间划分成一些称为单元的子空间,并且…

Python支持对整数和浮点数直接进行四则混合运算,运算规则和数学上的四则运算规则完全一致. 基本的运算: 1 + 2 + 3 # ==> 6 4 * 5 - 6 # ==> 14 7.5 / 8 + 2.1 # ==> 3.0375 使用括号可以提升优先级,这和数学运算完全一致,注意只能使用小括号,但是括号可以嵌套很多层: (1 + 2) * 3 # ==> 9 (2.2 + 3.3) / (1.5 * (9 - 0.3)) # ==> 0.4214559386973180…

[背景] 原文链接:http://blog.csdn.net/ordeder/article/details/25343633 Python整数对象是不可变对象,什么意思呢?比如运行例如以下python语句 >>>a = 1023 >>>a = 1024 >>>b = a >>>c = 1024 >>>d = 195 >>>e = 195 python的整数对象结构为: typedef struct…

整数线性规划问题的基本内容 整数线性规划解决的是自变量在一定的线性约束条件下,使得线性目标函数求得最大值或者最小值的问题.其中自变量只能取整数.特别地,当自变量只能取0或者1时,称之为 0-1 整数规划问题. 当目标函数为最小值时,上述问题可以写成如下形式: \[ \min z=\mathbf{F}^{T}\mathbf{X} \] \[ \text { s.t. } \left\{\begin{array}{l} {\mathbf{A}\mathbf{X} \leqslant \mathbf{…

00 前言 之前一直做启发式算法,最近突然对精确算法感兴趣了.但是这玩意儿说实话是真的难,刚好boss又叫我学学column generation求解VRP相关的内容.一看里面有好多知识需要重新把握,所以这段 时间就打算好好学学精确算法.届时会把学习过程记录下来,也方便大家学习! 01 什么是branch and bound(定义)? 1.1 官方一点[1] Branch and bound (BB, B&B, or BnB) is an algorithm design paradigm fo…

一.本章要点 延续让你可以回到程序执行当中之前的某个点; 可以在shift块中捕获延续 延续函数一直延展到包含它的reset块的尾部 延续所谓的”余下的运算“,从包含shift的表达式开始,到包含它的reset块的尾部结束,其中shift替换成一个”洞“ 当你传入一个参数来调用延续时,这个”洞“将传入的参数值填上 包含shift表达式的代码将被以”延续传递风格“(CPS)重写,直到包含该延续的reset为止 包含shift但没有reset的方法必须加上CPS注解 延续可以用来将对某个树形结构的递…

Python 五个知识点搞定作用域 1.块级作用域 想想此时运行下面的程序会有输出吗?执行会成功吗? #块级作用域 if 1 == 1: name = "lzl" print(name) for i in range(10): age = i print(age) 我们先看下执行结果 #输出 1 lzl 9 代码执行成功,没有问题:在Java/C#中,执行上面的代码会提示name,age没有定义,而在Python中可以执行成功,这是因为在Python中是没有块级作用域的,代码块里的变量…

引言:今天工作遇到了一个需要按行读取txt文件数据的需求,查询了一下getline()函数,发现这竟然是一个C++的标准库函数,而且设计的很好,特地做一下记录.getline本质是一个定界流输入截取函数,默认是换行符'/n' 个人技术博客(文章整理+源码): https://zobolblog.github.io/LearnWinAPI/ 最终效果: 1.getline函数声明 getline函数实际上是在两个模块都有实现,一个是string一个是istream,它们的用处都是一样的,只是输入的…

random 模块的高级玩法 1.python 随机产生姓名 方式一: import random xing = [ '赵', '钱', '孙', '李', '周', '吴', '郑', '王', '冯', '陈', '褚', '卫', '蒋', '沈', '韩', '杨', '朱', '秦', '尤', '许', '何', '吕', '施', '张', '孔', '曹', '严', '华', '金', '魏', '陶', '姜', '戚', '谢', '邹', '喻', '柏', '水', '…

概念 高德纳箭号表示法是种用来表示很大的整数的方法,由高德纳于1976年设计.它的意念来自幂是重复的乘法,乘法是重复的加法. 定义 计算 一个箭头 2↑3=2×2×2=8 2↑4=2×2×2×2=16 3↑3=3×3×3=27 a↑b=    两个箭头 2↑↑3=2↑2↑2(注意:此处要从右往左计算)=2↑4=16 3↑↑3=3↑3↑3=3↑27=  =7625597484987 4↑↑3=4↑4↑4=4↑256≍  a↑↑b=   三个箭头 2↑↑↑3=2↑↑2↑↑2=2↑↑2↑2=2↑↑4=…

python超强玩法--一行命令搞定局域网共享 ​ 今天刷到python的一个新玩法,利用python自带的http服务,快速创建局域网共享服务,命令如下: python -m thhp.server 9090 9090 是指定的的端口号,默认8080,也可以如上指定 效果如下: 服务端: client端: 速度也是超级快,个人测试达到了60M/s,可以说 相当可以了 操作方法: 在需要分享的目录的cmd窗口下执行上面命令,一件部署,前提是配置有python 的环境变量…

在python中,需要将整数均分成N等分.python divide integers N equal parts sum # 拆分整数 def split_integer(m, n): assert n > 0 quotient = int(m / n) remainder = m % n if remainder > 0: return [quotient] * (n - remainder) + [quotient + 1] * remainder if remainder < 0…

小整数对象池 (在python内置了) 整数在程序中的使用非常广泛,Python为了优化速度,使用了小整数对象池,避免为整数频繁申请和销毁内存空间. Python对小整数的定义是[-5,257]这些整数对象是提前建立好的,不会被垃圾回收.在一个Python的程序中,所有位于这个范围内的整数使用的都是同一个对象. 同理,单个字母也是这样的. 但是当定义2个相同的字符串时,引用计数为0,触发垃圾回收 intern机制 a1="helloworld" a2="helloworld&…

1.大最小贴近度评价法 概念: 贴近度表示两个模糊几何之间的彼此接近程度,在模糊模式识别方法中采用贴近度的大小识别待判别模糊子集的模式类别.为衡量待识别子集的类别,需要判别各个阶段与标杆模糊集合之间的相对贴近程度. 上表中第一列是优化(标杆)模糊集合(共15项因子),第二列是待判定的集合,第三列是(1)式的分子项,第四列是(1)式的分母项,贴近度=第三列数据求和/第四列数据求和=8.64/13.52=0.64 如果优化集改了,那么新的贴近度=0.52/10.64=0.05 python代码实现…

1. 小整数对象池 整数在程序中的使用非常广泛,Python为了优化速度,使用了小整数对象池, 避免为整数频繁申请和销毁内存空间. Python 对小整数的定义是 [-5, 256] 这些整数对象是提前建立好的,不会被垃圾回收.在一个 Python 的程序中,无论这个整数处于LEGB中的哪个位置, 所有位于这个范围内的整数使用的都是同一个对象.同理,单个字母也是这样的. In [1]: a=-5 In [2]: b=-5 In [3]: a is b Out[3]: True In [4]: a…

整数和浮点数有那个四则运算: 两种类型的数可以直接进行加减,当整数和浮点数进行加减的时候,结果会自动的变为浮点数,其中除法运算是“/”来表示的, 而余数的算术符号是“%”来表示的.如果是为了求除完后的整数,需要“//”. 在布尔值的判断中我们不需要进行这种符号的判断,直接使用那三种运算:与运算,非运算和或运算. 在Python中它把None,0和' '看成false,其他的数值和非字符串都看成true…

小整数对象池: Python为了优化速度,使用了小整数对象池,避免为整数频繁申请和销毁 Python 对小整数的定义是 [-5, 256] 这些整数对象是提前建立好的,不会被垃圾回收,所有位于这个范围内的整数使用的都是同一个对象,单个字母也是这样的. intern机制: 实现 Intern 机制的方式非常简单,就是通过维护一个字符串储蓄池,这个池子是一个字典结构,如果字符串已经存在于池子中就不再去创建新的字符串,直接返回之前创建好的字符串对象,如果之前还没有加入到该池子中,则先构造一个字符串对象…

#!/usr/bin/env python #-*- coding=utf-8 -*- small_ints = dict() for i in range(-10000,10000): small_ints[i] = id(i) for i in range(-10000,10000): if id(i) != small_ints[i]: del small_ints[i] print("[%s, %s]"%(min(small_ints.keys()), max(small_in…

日常的web自动化过程中,我们常常用python selenium库来操纵Chrome浏览器实现网页的自动化.这其中有个比较头疼的问题:Chrome的更新频率非常频繁,与之对应的Chromedriver版本也必须相应更新.如果两者版本的主版本号相差超过1,selenium则会报异常,大概错误信息如下: 发生异常: SessionNotCreatedException Message: session not created: This version of ChromeDriver only s…

编程测试题: 输入一个正整数将其分解成素数的乘积,输入格式连续输入m个数,然后将这m个数分别分解,如 输入: 2 10 20 输出: 2 5 2 2 5 Python code: def primes(n): primfac = [] d = 2 while d*d <= n: while (n % d) == 0: primfac.append(d) n //= d d += 1 if n > 1: primfac.append(n) return primfac s = int(raw_i…

函数格式 scipy.optimize.linprog(c, A_ub=None, b_ub=None, A_eq=None, b_eq=None, bounds=None, method='simplex', callback=None, options=None) 今天阅读数据建模第一章线性规划问题,问题描述如下: 通过介绍我们知道了线性规划,就是目标函数及约束条件均为线性函数. 通过画图我们可知,X1,X2的最优解为2,6,目标值为26. 我们如何时候这个scipy的公式来计算这个值呢:…

#!/usr/bin/env python#ecoding=utf-8'''Created on 2017年11月2日 @author: James zhan''' def fun(n): if n==1 or n==0: return 1 else: return fun(n-1)*ns = input('Input a init:')print fun(s)…

直接上代码: from numpy import *import numpy as npDrug_array = zeros((708,708),dtype = int)f = open('D:\mat_drug_drug.txt')lines = f.readlines()Drug_row = 0for line in lines: list = line.strip('\n').split(' ') Drug_array[Drug_row:] = list[:] Drug_row += 1p…

一行搞定-统计一句话中每个单词出现的个数 >>> s'i am a boy a bood boy a bad boy' 方式一:>>> dict([(i,s.split().count(i)) for i in s.split()]){'a': 3, 'boy': 3, 'i': 1, 'am': 1, 'bad': 1, 'bood': 1} >>> set([(i,s.split().count(i)) for i in s.split()])se…

二分查找算法也称折半查找,基本思想就是折半,和平时猜数字游戏一样,比如猜的数字时67,猜测范围是0-100,则会先猜测中间值50,结果小了,所以就会从50-100猜测,中间值为75,结果大了,又从50-75猜测中间值,一直到猜中为止.因此,二分查找有一个限制就是原先数组需要是一个有序数组.代码如下: ##二分查找算法 def binarysearch(a, num): length = len(a) low = 0 high = length - 1 while low <= high: mid…