删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

一种基于旁信息增强的协同过滤自动编码器模型

本站小编 Free考研考试/2020-03-23

刘海博1,2, 冯时1, 于戈1
1. 东北大学 计算机科学与工程学院, 辽宁 沈阳 110169;
2. 河北大学 网络空间安全与计算机学院, 河北 保定 071002
收稿日期:2018-12-20
基金项目:国家自然科学基金资助项目(61872074);河北省自然科学基金资助项目(F2015201140,F2017201208)。
作者简介:刘海博(1979-),男,河北保定人,东北大学博士研究生;
于戈(1962-), 男, 辽宁沈阳人, 东北大学教授,博士生导师。

摘要:针对用户评分预测不准确的问题, 提出了一种基于旁信息(side-information)对用户兴趣进行预测的协同过滤自动编码器推荐模型, 给出了模型的设计原理、损失函数以及具体结构.模型使用单隐藏层自动编码器实现, 用户评分与旁信息同为模型的输入/输出数据, 旁信息也直接参加模型的训练, 这种设计不仅降低了模型的规模和复杂度, 而且旁信息可以直接对用户兴趣进行修正.同时, 通过对训练数据集合的合理划分与扩充, 使得训练的网络模型增加了表达能力.在真实数据集上的对比实验表明, 本文提出的方法提高了评分预测的准确度, 具有一定的实用价值.
关键词:自动编码器协同过滤旁信息推荐系统评分预测
Auto-Encoder for Collaborative Filtering Based on Side-Information Enhancement
LIU Hai-bo1,2, FENG Shi1, YU Ge1
1. School of Computer Science & Engineering, Northeastern University, Shenyang 110169, China;
2. School of Cyber Security and Computer, Hebei University, Baoding 071002, China
Corresponding author: LIU Hai-bo, E-mail:liuhaibo@stumail.neu.edu.cn
Abstract: In view of inaccurate rating prediction of traditional recommendation model, a model based on side-information enhanced collaborative filtering auto-encoder was proposed, and the design principle, loss function and structure of the model were discussed in detail. The model was implemented by a single hidden layer auto-encoder network with ratings and side-information as input/output data. The model was directly trained by ratings and side-information. The proposed framework not only reduces the scale and complexity of the model, but also amends the user's interest with side-information. At the same time, the representation ability of the trained model is increased by reasonable division and expansion of training datasets. The comparisons on real datasets show that the proposed model with side-information can effectively ameliorate the user's interest predicted by the model without side-information, and improve the accuracy of rating prediction, which has a certain practicality.
Key words: auto-encodercollaborative filteringside-informationrecommender systemrating prediction
推荐系统根据用户的历史行为特征, 对用户进行分类, 并预测用户未来的兴趣与行为趋势, 进而向用户推荐其可能感兴趣的商品或信息, 提高用户的使用体验.神经网络是机器学习的一个分支研究领域.它通过单层或者多层次的神经网络学习表达数据的抽象特征, 可以解决有监督和无监督学习任务.伴随着深度神经网络的飞跃式发展, 它逐渐取代了传统的基于过滤技术的推荐方法.
在很多实际应用中, 系统不仅仅包含像样本评分或标签这样的监督信息, 还包含很多附加信息, 这些附加信息, 称之为旁信息(side-information)[1].以电影推荐为例:用户的性别、年龄、职业信息, 电影的类型、年份、演员、内容描述以及海报等都是旁信息.旁信息可以辅助系统对样本标签或评分数据进行处理, 在推荐技术中得到了比较广泛的使用.
评分预测[2]问题是推荐系统需要解决的主要问题之一.系统通过大量用户的历史评分记录, 预测其对某个项目的可能评分.这里主要讨论如何利用项目的旁信息来增强神经网络中的自动编码器模型对用户评分的预测以及模型规模简化问题.
目前, 科技工作者对基于神经网络的推荐模型的研究已经取得了一些成果, 提出了多种解决方案.其中, 多层感知机与自动编码器[3]在基于神经网络模型的推荐系统中应用较多.
文献[4]提出了一种神经协同过滤(neural collaborative filtering, NCF)网络, NCF是一种基于多层感知机, 并利用用户和项目的旁信息作为用户和项目的潜在特征进行建模的推荐模型.NCF的预测规则为, f(·)是多层感知器, θ是该网络的参数, suusersiitem分别是用户与项目的旁信息.损失函数为加权平方误差: .文献[5]提出了神经社会协同排序推荐系统, 将NCF模型推广到跨领域的社会推荐中.文献[6]提出了跨域内容增强协同过滤模型, 网络的最后一层使用点积, 混合了基于旁信息的用户与项目特征.文献[7]提出了一种宽度与深度学习模型, 宽度学习部分是单层感知器, 可以从历史数据中获得直接特征, 通过产生更一般和抽象的特征表示来获得泛化能力, 深度学习部分是多层感知器.
在基于自动编码器的推荐方法中, 提取用户与项目的抽象特征以及填充用户-项目评分矩阵中的缺失值, 是通过自动编码器完成推荐的两个主要手段.文献[8]提出了基于用户与基于项目的自动编码器推荐模型, 作者通过编码器将用户评分向量映射到抽象低维特征表示空间, 并通过解码器将其恢复到与输入特征空间相同的输出特征空间, 从而预测用户缺失的项目评分.文献[9]提出了基于协同过滤的自动编码器模型, 并且结合了用户和项目旁信息, 以减轻数据的稀疏性和冷启动问题.文献[10]提出了协同降噪自动编码器, 使用单隐含层自动编码器构建推荐模型, 在输入层用户偏好向量被格式化为二进制向量, 并且加入高斯噪声.文献[11]将变分自动编码器推广到协同过滤中.文献[12]通过有效的全反馈训练机制训练了深层自动编码器, 提高了评分预测的精度.
上述研究对于神经网络以及旁信息在推荐系统中的应用取得了一定的进展.但是其研究成果, 大部分都是只使用用户-项目评分信息进行模型训练, 在小部分引入旁信息的方法中, 基本上都是利用单独的深层网络来挖掘旁信息的抽象特征, 进而辅助用户-项目评分矩阵的挖掘.其方法增加了网络的规模与复杂度, 提高了推荐系统的代价.针对上述问题, 本文提出了一种基于旁信息对用户兴趣进行修正和预测的协同过滤自动编码器模型.真实数据集上的对比实验表明, 本文设计的模型降低了计算规模与复杂度, 提高了评分预测的准确度.
1 推荐模型1.1 传统自动编码器模型自动编码器是一种数据压缩算法, 是一种从样本中自动学习的自监督学习算法.模型主要分为编码器与解码器.其中, 数据的编码部分和解码部分是数据相关的, 编码部分通过神经网络将高维的输入数据映射到低维度的抽象特征空间, 解码部分同样通过神经网络将低维度的抽象特征数据还原为高维度的原始数据.因此, 自动编码器是有损失的.其终极目标是实现输出与输入无限接近.自动编码器可表示为h=f(x)=W·x+bx, f(·)为编码映射, W为编码权重矩阵, x为样本原始输入信息向量, bx是编码偏置向量, h为输入信息经过编码后得到的抽象特征向量.解码器可表示为y=g(h)=WT·h+by, g(·)为解码映射, WT为解码权重矩阵, by为解码偏置向量, y为经过解码后得到的重构信息向量.因此, 其目标函数为, xm表示样本集合中第m个样本的原始信息, l表示损失函数.传统单隐含层自动编码器结构如图 1所示.
图 1(Fig. 1)
图 1 传统单隐含层自动编码器结构图Fig.1 Traditional single hidden layer auto-encoder

1.2 基于旁信息增强的协同过滤自动编码器模型传统的基于自动编码器的推荐模型计算复杂度高, 并且只是简单地要求目标输出与输入尽量相同, 模型尽量拟合用户的现有兴趣, 而很少考虑用户的未来兴趣.
本文提出的基于旁信息对用户兴趣进行增强的协同过滤自动编码器模型(auto-encoder model for collaborative filtering based on side-information enhancement)将项目的旁信息抽象为隐藏特征, 并与用户的抽象兴趣特征直接融合, 从而增强了自动编码器对用户兴趣的预测能力.对于自动编码器网络, 网络层次越深, 则特征抽取越抽象.为了降低模型的规模, 并且增强项目旁信息对用户兴趣的直接修正作用, 模型采用单隐含层神经网络实现.把数据集中每个用户与项目的交互数据按照时间排序, 并将其划分为若干份.通过对数据集的合理划分, 使得训练的模型可以反映用户的未来兴趣.假设用户集合为U={u1, …, uM}, M为用户数量, 项目集合为I={i1, …, iN}, N为项目数量, 部分可观测的用户-项目评分矩阵为RRM×N, 每一个用户um的可观测评分向量为rm=(Rm, 1, …, Rm, N)∈RN.设项目拥有K种旁信息, 则所有种类旁信息的集合为S={S1, …, SK}, 每一类旁信息集合大小不同, 设第k类旁信息集合大小为Lk, 则项目的第k类旁信息集合表示为Sk={sk, 1, …, sk, Lk}.用户um的第k类旁信息向量表示为xkm=(wm, k, 1, …, wm, k, Lk)∈RLk, 对于任意的, Im为用户um可观测评分的项目集合, 如果in拥有旁信息sk, Lk, 则wk, Lkn=1, 否则wk, Lkn=0.
基于旁信息增强的协同过滤自动编码器模型结构如图 2所示.对于用户um, 用户-项目评分向量rm, 以及所有旁信息向量xkm同时为模型的输入, 通过式(1)完成编码器抽象并在隐藏层融合.
(1)
图 2(Fig. 2)
图 2 基于旁信息增强的协同过滤自动编码器模型Fig.2 Auto-encoder model for collaborative filtering based on side-information enhancement

其中:hm为用户um的抽象特征表示; Wrh为评分编码权重; Wkh为旁信息编码权重; brh为隐藏层偏置项.
对于重构层, 分别通过式(2)与式(3)完成对评分信息与旁信息的重构.
(2)
(3)
其中: 表示经过解码得到的重构评分信息; Wrd为评分解码权重; brd为评分解码偏置项; 表示经过解码得到的旁信息; Wkd为旁信息解码权重; bkd为旁信息解码偏置项.
2 损失函数损失函数用于评价模型的预测值与真实值的不一致程度, 因此, 需要依靠损失函数来选取最优化模型.采用可观测值的平方误差(masked squared error, MSE)作为损失函数.与平方误差损失函数不同, 可观测的平方误差只考虑目标集合中用户已评分项目的平方和误差.损失函数定义为
(4)
其中: 为目标评分; 为目标旁信息; λk为第k类旁信息的损失权重超参数; ‖·‖O2表示可观测值的平方和误差.模型在训练时, 目标重构数据信息与输入数据信息不一定相同(输入数据集合与目标数据集合不一定相等), 因此, rm, xkm不一定相等.为了防止模型过拟合, 损失函数中加入正则化项Lreg, 定义为
(5)
其中:λ为损失函数的正则化项超参数; ‖·‖F2表示矩阵各元素绝对值平方和(即不进行开平方运算的弗洛贝尼乌斯范数).
3 实验及分析3.1 数据集与数据划分使用MovieLens 10M(ML-10M)与MovieLens 20M(ML-20M)两个数据集, 进行实验.每个用户的评分数据按时间顺序排列, 在最后的20%数据中, 随机选择50%作为验证集, 另外50%作为测试集, 前面的80%数据作为训练集.在测试集与验证集中删除了没有在训练集中出现的电影.为了增强模型对用户未来兴趣的表达能力, 训练集按时间顺序等分4份:T1, T2, T3, T4.训练时输入数据集合为{(T1), (T1, T2), (T1, T2, T3), (T1, T2, T3, T4)}, 对应的目标集合为{(T1, T2), (T1, T2, T3), (T1, T2, T3, T4), (T1, T2, T3, T4)}.即, 采用的是目标集合大于输入集合的增量方式进行训练.由于电影的旁信息是附加在电影中的, 因此, 对于电影的旁信息也采用同样的划分方式.
旁信息的选择.影响用户对电影评分的主要因素包含两类:电影内容以及电影演员阵容(明星效应).由于电影内容信息的特征提取涉及到文本内容分析, 必须使用单独的网络进行特征提取, 这样增加了模型的复杂度, 因此采用电影分类信息来代替内容信息.使用词袋法分别对电影分类信息和演员信息进行处理.数据集详细信息见表 1.
表 1(Table 1)
表 1 数据集信息Table 1 Datasets information
集合属性 旁信息 ML-10M ML-20M
训练集合 用户 69 878 138 493
电影 10 352 22 085
评分 7 972 582 15 945 812
分类 20 20
演员 12 290 27 549
验证集合 用户 34 939 69 246
电影 10 117 19 462
评分 1 010 185 2 029 698
分类 20 20
演员 11 921 23 661
测试集合 用户 34 939 69 246
电影 10 073 18 685
评分 1 016 227 2 016 633
分类 20 20
演员 11 872 22 931
稀疏度 0.013 8 0.006 5


表 1 数据集信息 Table 1 Datasets information

3.2 评价指标及实验参数经过实验得出模型最佳参数, 如表 2所示.
表 2(Table 2)
表 2 模型参数值Table 2 Model parameter values
参数
损失函数正则化项的超参数λ 0.1
电影分类损失权重超参数λ0 0.1
电影演员损失权重超参数λ1 0.1
模型学习率 0.001
神经元丢弃率 0.5
网络参数优化方法 Adam
神经元激活函数 Relu
隐藏层神经元数量 128


表 2 模型参数值 Table 2 Model parameter values

本文采用可观测的均方根误差(masked root mean squared error, MRMSE)作为评价方法的性能指标.MRMSE是一种常用的测量评分之间差异的量度.这里的是目标集上的真实值, 是预测值, LM是目标集合中所有可观测的真实评分数量.
(6)
3.3 实验结果为了验证本文方法以及旁信息与数据划分对评分预测的影响, 同时设计了4个不同的方案进行对比实验:①AE, 使用原始自动编码器模型, 不使用旁信息与数据划分; ②AE-d, 使用本文模型, 使用数据划分, 但不使用旁信息; ③AE-s, 使用本文模型, 使用旁信息, 但不使用数据划分; ④AE-sd, 使用本文模型, 同时使用旁信息与数据划分.为了更加直观表达实验结果, 展示训练数据集以及其对应的验证数据集上的MRMSE值.实验结果如图 3图 4所示.
图 3(Fig. 3)
图 3 ML-10M数据集实验结果Fig.3 Experimental results of ML-10M dataset (a)—ML-20M训练结果;(b)—ML-20M估计结果.

图 4(Fig. 4)
图 4 ML-20M数据集实验结果Fig.4 Experimental results of ML-20M dataset (a)—ML-20M训练结果;(b)—ML-20M估计结果.

图 3a可以观察到使用原始自动编码器模型的AE的实验误差要高于AE-d, AE-s与AE-sd, 而且AE与AE-d比较接近, 但是与AE-s以及AE-sd相差较大.由图 3b可以发现, AE与AE-d在验证集合上的误差差距要远高于其在训练集上的误差差距.由于数据划分方法以及增量式的训练模式, 提高了模型对用户未来兴趣的预测能力.AE-s的误差也远小于AE, 这是因为, 项目的旁信息可以映射用户的兴趣特征, 因此, 旁信息的使用有效地修正了仅仅使用评分信息进行训练得到的用户兴趣特征.此外, 还可以发现, 在训练集上AE-sd的性能指标低于AE-s的性能指标, 而在验证集上AE-sd的性能指标略高于AE-s的性能指标.这是因为, 对于旁信息, 使用了同样的数据划分, 这样增强了旁信息对用户未来兴趣的预测能力.由图 3b还可以发现, 在模型的预测能力的平稳性上, AE-sd最好, AE-s次之, AE-d的波动比较大.这是因为, 旁信息对用户的兴趣进行修正的同时也增强了模型的稳定性.
图 4描述了4个方法在数据集ML-20M上的实验结果.实验结果与其在数据集ML-10M上的结果类似.但是, 需要注意的是在ML-20M上的训练结果可以看出AE与AE-d的性能很接近, 它们的实验差距要小于其在ML-10M训练集上的实验差距.同样, AE与AE-d在ML-20M验证集上的实验结果的差距也小于在ML-10M验证集上的差距.这主要是因为ML-20M要比ML-10M稀疏.稀疏的数据减弱了数据划分的作用.但是, 总体上看, 训练数据集的划分都起到了提高评分预测准确度的作用.
在测试集上的实验结果如表 3所示.可以发现AE-sd取得了较好的评分预测结果.
表 3(Table 3)
表 3 测试集实验结果Table 3 Testing results
数据集 AE AE-d AE-s AE-sd
ML-10M 0.921 8 0.852 2 0.829 2 0.824 7
ML-20M 0.917 2 0.856 2 0.826 9 0.819 6


表 3 测试集实验结果 Table 3 Testing results

4 结语本文针对如何利用旁信息增强自动编码器推荐模型的预测能力进行了研究.首先基于用户或者项目的旁信息可以映射用户兴趣的假设, 设计了基于旁信息增强的协同过滤自动编码器推荐模型.为了使旁信息可以直接修正用户的兴趣, 使用单隐含层自动编码器网络实现推荐模型, 并且旁信息与评分信息共享中间隐含层.同时, 将训练数据集合按时间等分, 训练时采用目标集合大于输入集合的增量训练的方法.实验表明本文模型可以降低模型的复杂度, 并且提高评分预测的准确度.
参考文献
[1] Elhamifar E, Wang Y.High rank matrix completion with side information[C]// AAAI Conference on Artificial Intelligence.New Orleans: AAAI, 2018: 4252-4259.
[2]Adomavicius G, Tuzhilin A. Toward the next generation of recommender systems:a survey of the state-of-the-art and possible extensions[J].IEEE Transactions on Knowledge and Data Engineering, 2005, 17(6): 734–749.DOI:10.1109/TKDE.2005.99
[3]Zhang S, Yao L, Sun A, et al. Deep learning based recommender system:a survey and new perspectives[J].ACM Computing Surveys, 2018, 1(1): 1–35.
[4] He X, Liao L, Zhang H, et al.Neural collaborative filtering[C]//International Conference on World Wide Web.Perth: ACM, 2017: 173-182.
[5] Wang X, He X, Nie L, et al.Itemsilk road: recommending items from information domains to social users[C]//International ACM SIGIR Conference on Research and Development in Information Retrieval.Tokyo: ACM, 2017: 185-194.
[6] Lian J, Zhang F, Xie X, et al.CCCF net: a content-boosted collaborative filtering neural network for cross domain recommender systems[C]//International Conference on World Wide Web Companion.Perth: ACM, 2017: 817-818.
[7] Cheng H T, Koc L, Harmsen J, et al.Wide & deep learning for recommender systems[C]// The Workshop on Deep Learning for Recommender Systems.Boston: ACM, 2016: 7-10.
[8] Sedhain S, Menon A K, Sanner S, et al.AutoRec: autoencoders meet collaborative filtering[C]// International Conference on World Wide Web.Florence: ACM, 2015: 111-112.
[9] Strub F, Mary J.Collaborative filtering with stacked denoising autoencoders and sparse inputs[C]// Neural Information Processing Systems.Montreal: NIPS, 2015: 136-144.
[10] Wu Y, Dubois C, Zheng A X, et al.Collaborative denoising auto-encoders for top-n recommender systems[C]// International Conference on Web Search and Data Mining.San Francisco: ACM, 2016: 153-162.
[11] Li X, She J.Collaborative variational autoencoder for recommender systems[C]// ACM SIGKDD International Conference.Halifax: ACM, 2017: 305-314.
[12]Kuchaiev O, Ginsburg B. Training deep autoencoders for collaborative filtering[J].arXiv:Machine Learning, 2017, 1(1): 1–5.

相关话题/信息 编码器

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于酒店特征信息和在线评价信息的酒店选择
    陶玲玲,尤天慧,袁媛东北大学工商管理学院,辽宁沈阳110169收稿日期:2019-01-31基金项目:国家自然科学基金资助项目(71771043)。作者简介:陶玲玲(1991-),女,辽宁铁岭人,东北大学博士研究生;尤天慧(1967-),女,黑龙江宾县人,东北大学教授,博士生导师。摘要:为辅助游客通 ...
    本站小编 Free考研考试 2020-03-23
  • 基于双互信息准则的雷达自适应波形设计方法
    辛凤鸣1,汪晋宽2,王彬1,李梅梅11.东北大学秦皇岛分校计算机与通信工程学院,河北秦皇岛066004;2.东北大学信息科学与工程学院,辽宁沈阳110819收稿日期:2018-12-05基金项目:国家自然科学基金资助项目(61601109);中央高校基本科研业务费专项资金资助项目(N18230402 ...
    本站小编 Free考研考试 2020-03-23
  • 基于改进分层激励控制线的多阶段信息集结方法
    李玲玉1,2,郭亚军1,易平涛1,冯雪丽11.东北大学工商管理学院,辽宁沈阳110169;2.南昌大学经济管理学院,江西南昌330031收稿日期:2016-06-28基金项目:国家自然科学基金资助项目(71671031,71473033,71701040);中央高校基本科研业务费专项资金资助项目(N ...
    本站小编 Free考研考试 2020-03-23
  • 基于缓存开销的信息中心网络缓存协作策略
    蔡凌1,汪晋宽2,王兴伟3,韩来权41.东北大学秦皇岛分校控制工程学院,河北秦皇岛066004;2.东北大学信息科学与工程学院,辽宁沈阳110819;3.东北大学软件学院,辽宁沈阳110169;4.东北大学秦皇岛分校计算机与通信工程学院,河北秦皇岛066004收稿日期:2016-08-10基金项目: ...
    本站小编 Free考研考试 2020-03-23
  • 基于功能网络信息传播预测疾病-miRNAs的关联
    李建华,雒士源,张建营,康雁东北大学中荷生物医学与信息工程学院,辽宁沈阳110169收稿日期:2016-10-20基金项目:国家自然科学基金资助项目(61372014)。作者简介:李建华(1973-),男,河北怀来人,东北大学讲师,博士;康雁(1964-),男,辽宁沈阳人,东北大学教授,博士生导师。 ...
    本站小编 Free考研考试 2020-03-23
  • 转子系统常开裂纹故障信息提取方法的研究
    刘杨,石拓,姚红良,闻邦椿东北大学机械工程与自动化学院,辽宁沈阳110819收稿日期:2016-09-21基金项目:中央高校基本科研业务费专项资金资助项目(N140304005,N160313003);中国博士后科学基金资助项目(2014M551105,2015T80269)。作者简介:刘杨(198 ...
    本站小编 Free考研考试 2020-03-23
  • 社交网络信息源快速定位方法
    张聿博,张锡哲,徐超,张斌东北大学计算机科学与工程学院,辽宁沈阳110819收稿日期:2015-05-14基金项目:国家科技支撑计划项目(2014BAI17B00);国家关键科技研发基金资助项目(2015BAH09F02,2015BAH47F03);中央高校基本科研业务费专项资金资助项目(N1404 ...
    本站小编 Free考研考试 2020-03-23
  • 基于UML和XML的工艺族信息模型表达方法
    王琳1,张永健1,钟诗胜21.哈尔滨工业大学(威海)船舶与海洋工程学院,山东威海2642092.哈尔滨工业大学机电工程学院,黑龙江哈尔滨150001收稿日期:2014-10-10基金项目:国家自然科学基金资助项目(51075083);教育部高等学校博士学科点专项科研基金资助项目(2011230213 ...
    本站小编 Free考研考试 2020-03-23
  • 灰色信息攻击下相依网络的鲁棒性
    朱潜1,2,王一帆1,朱志良1,任涛11.东北大学软件学院,辽宁沈阳110819;2.东北大学计算机科学与工程学院,辽宁沈阳110819收稿日期:2015-12-23基金项目:国家自然科学基金资助项目(61473073,61374178,61300019);辽宁省博士启动基金资助项目(2015011 ...
    本站小编 Free考研考试 2020-03-23
  • 有限需求信息下基于最大熵原理的风险厌恶库存模型
    邱若臻,苑红涛,冯俏东北大学工商管理学院,辽宁沈阳110167收稿日期:2015-06-29基金项目:国家自然科学基金资助项目(71372186);教育部人文社会科学研究一般项目(11YJC630165);中央高校基本科研业务费专项资金资助项目(N150604005).摘要:针对风险厌恶的库存决策者 ...
    本站小编 Free考研考试 2020-03-23