- 博客(96)
- 收藏
- 关注
RSS订阅原创 84、浅谈原始对偶算法
原始对偶算法理解过很多次，每一个的理解都不一样，这次将其记录一下这里定义原始问题为资源有限下的价值最大化模型，A是资源消耗矩阵，b是资源量，c是产品价值向量max cTxs.t. Ax <= b对偶的定义是啥呢，它先组合一下原始的约束，组合权重用y表示，那么组合出来的向量为ATy，我约束其在每个维度的值都大于c，强行赋一个定义可以是保证价值情况下的资源最小化，这样对偶问题就被构造出来了，不难可以知道bTy >= xTATy>= xTc，min bTys.t. .
2021-11-26 19:13:07
1769
原创 83、选址结构梳理
import random as rdimport mathimport timefrom TSP.Heuristic import is_overlapfrom tool.Convex_hull_find import boundary_findfrom TSP.Heuristic import TSP_insertfrom VRP.CVRP.tool.Graph_transform import *from pyscipopt import Model, quicksumimport .
2021-11-22 17:23:55
354
1
原创 82、时间装饰器
import timeloop = 1def timeLog(f): def wrapper(*args, **kw): global loop now = time.perf_counter() res = f(*args, **kw) print("%s-%s:"%(loop, f.__name__), time.perf_counter()-now) return res return wrapp.
2021-11-15 10:43:41
130
原创 81、数独求解（策略求解、整数规划建模求解）
import refrom pyscipopt import Model, quicksum#set存value的可能性，当value值确认后，set清空class Cell: def __init__(self): self.value = "_" self.set = {1,2,3,4,5,6,7,8,9} def __repr__(self): return str(self.value)#读数据def dataRead.
2021-11-12 17:58:22
1035
原创 80、ALNS实现
from VRP.CVRP.tool.Value import *from VRP.CVRP.method.Init.Insert import *from VRP.CVRP.tool.Graph_transform import *class Input: def __init__(self, points, capacity, XBest, valueBest): self.points = points self.capacity = capacity.
2021-10-08 17:52:11
237
原创 79、多边形点边关系、求外包
from Project.PaixianOnline.basic.Coordinate import *from tool.Distance import *from tool.Convex_hull_find import boundary_findimport mathimport matplotlib.pyplot as plt#判断点是否在多边形上def isOnEdge(point, P): x, y = point.longitude, point.latitude .
2021-08-11 19:45:14
321
原创 78、无向图、有向图、序列相互转化
#有向图转无向图def Directed_to_undirected(A): return [[1 if A[i][j] == 1 or A[j][i] == 1 else 0 for j in range(len(A[i]))] for i in range(len(A))]#有向图转序列def Directed_to_sequence(A): return Undirected_to_sequence(Directed_to_undirected(A))#序列转有向图，默认L
2021-06-24 14:36:51
807
原创 77、CVRP方法（整数规划、横扫式聚类、层次聚类、扫描法、T2C）
整数规划建模import random as rdfrom pyscipopt import Model, quicksumfrom tool.Distance import calculate_distancefrom VRP.CVRP.tool.Plot import *from VRP.CVRP.tool.Print import *from VRP.CVRP.tool.Data_generate import *#一维装箱整数规划def IP_1d(points, capacity
2021-06-24 14:34:10
1173
原创 77、经典VRP建模
import random as rdimport matplotlib.pyplot as pltimport mathfrom pyscipopt import Model, quicksumfrom tool.Distance import calculate_distancedef points_generate(x_min, y_min, x_max, y_max, nums): points = [[0, (x_max+x_min)//2, (y_max+y_min)//.
2021-05-24 15:42:15
1419
原创 76、多边形一些基本操作（自相交、尖刺、保证逆时针、求交）
pro = """1、造一些多边形和一个大的范围2、多边形交集部分归属问题3、多边形未达部分归属问题"""import random as rdimport mathimport matplotlib.pyplot as pltfrom tool.Convex_hull_find import boundary_find#多边形生成def Polygon_generate(x_min, y_min, x_max, y_max, num_point): L = [[int(.
2021-04-20 20:52:14
1295
1
原创 75、QGIS学习记录
在console画多边形wkt2 = 'POLYGON ((11717491.32085317 4192105.64176356, 11717299.80128882 4193264.57267995, 11717299.8012888 4193264.5726799, 11717078.05634601 4192351.42419002, 11715386.86564284 4193357.23269148, 11715382.15507762 4193327.5158613, 11715382.1
2021-04-07 14:35:56
194
2
原创 73、Titanic
pro = """Variable Definition Key survival Survival 0 = No, 1 = Yes, 离散pclass Ticket class 1 = 1st, 2 = 2nd, 3 = 3rd, 离散sex Sex .
2021-03-12 11:22:11
129
原创 72、Cutting Stock Problem
from pyscipopt import Model, quicksumdef IP(A, B, time_limit): m, n = len(A), len(A[0]) model = Model("IP_CSP") #切分方式分布 X = [model.addVar(vtype="I", name="x[%s]" % j) for j in range(n)] #分布和最小 model.setObjective(quicksum(X[j] .
2021-03-05 11:56:51
574
原创 71、有关短视频冷启动的一些想法
冷启动分用户冷启动和物品冷启动，这里说说对物品冷启动的理解新的物品如视频、商品，由于没有足够的行为数据，很难利用统计的信息作推荐，一般有两类思路一种通过保量的方式，给予每个物品一个基准的流量，获取到一定的行为数据，再利用这部分的数据快速迭代筛选 一种通过文本和图像相似性，通过i2i召回的方式将新物品召回这两种方案都可以使得新物品有一定透出的几率，也可以结合着一起使用文本相似度有bert，图像相似度的有多模态...
2021-02-18 17:44:16
183
原创 70、品类规划小结
品类规划是选品方案的一种辅助工具，其从品类的维度去约束选品相关调研最早较为成体系介绍品类规划概念的是杨滔在14年的“聚划算商品销量预测模型”，其源于聚划算业务增长下，运营小二人工选品难以避免经验限制、体力限制和腐败等一系列问题，为建立精准高效的选品机制，品类规划被作为选品的前置流程，其提出使用商品单UV价值作为选品的衡量指标（该指标也可作为品类规划的衡量指标）。授武在“品类规划整改项目小结”中，提出聚划算的商业逻辑核心为“有限坑位限制下的GMV最大化”，且基于可解耦，灵活支持业务，稳健性好的原则
2021-01-11 09:59:42
1401
原创 69、选品小结
选品的结果会作为推荐的输入底池，其通常用在两个方面:1、商品数量有限制的场景，如营销、大促等场景，这时的商品数有限，必须做一定的筛选2、需要精选的场景，这时选品的工作本质上是在承担召回侧的压力选品结合投放来评估效率，这里建模如下：同时提出了一种优化解法其他一些思考：...
2021-01-11 09:51:32
479
原创 68、Kaggle竞赛 预测未来的销售
https://www.kesci.com/mw/dataset/5e981cf7ebb37f002c6007d6/content
2021-01-11 09:30:20
433
原创 67、矩阵分解的梯度下降求解、随机梯度下降
import mathimport random as rddef Loss(U, I, D, d): loss = 0 for u, i in D: loss += pow(sum([U[u][t]*I[i][t] for t in range(d)])-D[u, i], 2) return lossdef step_search(U, dU, I, dI, x, d, D): U1 = {u: [U[u][t] + x * dU[u][t] f.
2020-12-17 17:03:24
202
原创 66、线性模型在梯度下降、tf回归、MLP优化求解表现
# -*- coding: utf-8 -*-import tensorflow as tfimport random as rdimport mathimport numpy as nptf.compat.v1.disable_eager_execution()se = tf.compat.v1.Session()rd.seed(10)def sample_generate(n, low, sup, w, e): X = [[(sup - low) * rd.random() .
2020-12-16 20:04:03
160
原创 65、Top-N推荐系统的协同过滤方法
数据为movielen数据，这里仅用其中的ratings数据1、数据按年份切分import refrom datetime import datetimeimport timeimport pandas as pddef data_read(): t0 = time.perf_counter() ratings_path = r"D:\\!data\\3-recommender\\27-ml-25m\\ml-25m\\ratings.csv" # 文件较大，
2020-11-30 18:17:22
1500
1
原创 64、时间序列预测MA、WMA、EMA、AR、ARMA、ARIMA
1、移动平均MA采用历史平均作预测，WMA和EMA对远近的不同采用不同的权重#移动平均适用于不快速增长也不快速下降，且不存在季节性因素的场景import random as rdimport matplotlib.pyplot as pltimport mathrd.seed(1)def dot_product(a, b): if len(a) == 0 or len(b) == 0 or len(a) != len(b): raise Exception("dot
2020-11-30 10:44:55
904
原创 63、tensorflow简单demo
import tensorflow as tffrom numpy.random import RandomStatetf.compat.v1.disable_eager_execution()se = tf.compat.v1.Session()def print_tf(s): print(se.run(s))def demo1(): batch_size = 8 # 定义训练数据batch的大小 w1 = tf.Variable(tf.compat.v.
2020-07-07 21:47:48
176
原创 62、分组背包问题（MCKP）Greedy、Dyer-Zemel、对偶界、动态规划实现
原文链接：https://link.springer.com/chapter/10.1007/978-3-540-24777-7_11不同于传统的背包问题，MCKP需要保证每个分组里面至少选一个物品，建模如下其提出线性支配（LP-dominated）的定义如下这是个筛选的过程，把处于支配（同一组内存在比你重量小但价值高的物品）的去掉，线性支配的定义会更强，比如3和7，其被线性支配，故被去掉，我个人理解这样在部分情况有可能把最优解去掉了，所以一般采用前者。纯贪心算法如下..
2020-06-24 16:15:18
3988
1
原创 60、AUC计算和ROC曲线画图
import random as rdimport matplotlib.pyplot as pltdef sample_generate(num_p, num_n): return [[rd.gauss(mu=0.6, sigma=0.2), 1] for _ in range(num_p)] + [[rd.gauss(mu=0.4, sigma=0.2), 0] for _ in range(num_n)]def ROC_print(samples): samples.s.
2020-05-20 21:12:10
700
原创 59、用户近邻模型
import random as rdimport mathdef AB(A, B): return [[sum([A[i][k]*B[k][j] for k in range(len(A[0]))]) for j in range(len(B[0]))] for i in range(len(A))]def Ax(A, x): return [sum([A[i][j]...
2020-03-18 17:12:37
155
原创 58、无约束优化框架
import mathdef mod(v): return math.sqrt(sum(v[i]*v[i] for i in range(len(v))))def x_next(x, dx, step): return [x[i]+dx[i]*step for i in range(len(x))]#确定步长def step_generate(x0, f, dx, ...
2020-03-13 11:04:57
247
原创 57、SA框架
import copyimport random as rdimport mathimport timedef is_better(r0, r1, r): if r == True: if r0 < r1: return True else: return False else: ...
2020-03-13 10:05:24
331
原创 56、2-opt框架
import copyimport timedef swap(X, i, j): X1 = copy.deepcopy(X) X1[i], X1[j] = X1[j], X1[i] return X1def is_better(r0, r1, r): if r == True: if r0 < r1: ret...
2020-03-13 10:04:43
481
原创 55、python画热力图
from folium.plugins import HeatMap, Fullscreenimport foliumdef draw_map_heat2(data_heat, output_path='.'): lats = [d[0] for d in data_heat] lngs = [d[1] for d in data_heat] lat_mid = ...
2020-03-12 16:11:52
666
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人