提交 70bc4acb 编写于 作者: Y Yu Yang

V2 API Train

上级 58a33469
.idea
.ipynb_checkpoints
params.pkl
{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"collapsed": true,
"deletable": true,
"editable": true
},
"source": [
"# 个性化推荐\n",
"\n",
"本教程源代码目录在[book/recommender_system](https://github.com/PaddlePaddle/book/tree/develop/recommender_system), 初次使用请参考PaddlePaddle[安装教程](http://www.paddlepaddle.org/doc_cn/build_and_install/index.html)。\n",
"\n",
"## 背景介绍\n",
"\n",
"在网络技术不断发展和电子商务规模不断扩大的背景下,商品数量和种类快速增长,用户需要花费大量时间才能找到自己想买的商品,这就是信息超载问题。为了解决这个难题,推荐系统(Recommender System)应运而生。\n",
"\n",
"个性化推荐系统是信息过滤系统(Information Filtering System)的子集,它可以用在很多领域,如电影、音乐、电商和 Feed 流推荐等。推荐系统通过分析、挖掘用户行为,发现用户的个性化需求与兴趣特点,将用户可能感兴趣的信息或商品推荐给用户。与搜索引擎不同,推荐系统不需要用户准确地描述出自己的需求,而是根据分析历史行为建模,主动提供满足用户兴趣和需求的信息。\n",
"\n",
"传统的推荐系统方法主要有:\n",
"\n",
"- 协同过滤推荐(Collaborative Filtering Recommendation):该方法收集分析用户历史行为、活动、偏好,计算一个用户与其他用户的相似度,利用目标用户的相似用户对商品评价的加权评价值,来预测目标用户对特定商品的喜好程度。优点是可以给用户推荐未浏览过的新产品;缺点是对于没有任何行为的新用户存在冷启动的问题,同时也存在用户与商品之间的交互数据不够多造成的稀疏问题,会导致模型难以找到相近用户。\n",
"- 基于内容过滤推荐[[1](#参考文献)](Content-based Filtering Recommendation):该方法利用商品的内容描述,抽象出有意义的特征,通过计算用户的兴趣和商品描述之间的相似度,来给用户做推荐。优点是简单直接,不需要依据其他用户对商品的评价,而是通过商品属性进行商品相似度度量,从而推荐给用户所感兴趣商品的相似商品;缺点是对于没有任何行为的新用户同样存在冷启动的问题。\n",
"- 组合推荐[[2](#参考文献)](Hybrid Recommendation):运用不同的输入和技术共同进行推荐,以弥补各自推荐技术的缺点。\n",
"\n",
"其中协同过滤是应用最广泛的技术之一,它又可以分为多个子类:基于用户 (User-Based)的推荐[[3](#参考文献)] 、基于物品(Item-Based)的推荐[[4](#参考文献)]、基于社交网络关系(Social-Based)的推荐[[5](#参考文献)]、基于模型(Model-based)的推荐等。1994年明尼苏达大学推出的GroupLens系统[[3](#参考文献)]一般被认为是推荐系统成为一个相对独立的研究方向的标志。该系统首次提出了基于协同过滤来完成推荐任务的思想,此后,基于该模型的协同过滤推荐引领了推荐系统十几年的发展方向。\n",
"\n",
"深度学习具有优秀的自动提取特征的能力,能够学习多层次的抽象特征表示,并对异质或跨域的内容信息进行学习,可以一定程度上处理推荐系统冷启动问题[[6](#参考文献)]。本教程主要介绍个性化推荐的深度学习模型,以及如何使用PaddlePaddle实现模型。\n",
"\n",
"## 效果展示\n",
"\n",
"我们使用包含用户信息、电影信息与电影评分的数据集作为个性化推荐的应用场景。当我们训练好模型后,只需要输入对应的用户ID和电影ID,就可以得出一个匹配的分数(范围[1,5],分数越高视为兴趣越大),然后根据所有电影的推荐得分排序,推荐给用户可能感兴趣的电影。\n",
"\n",
"```\n",
"Input movie_id: 1962\n",
"Input user_id: 1\n",
"Prediction Score is 4.25\n",
"```\n",
"\n",
"## 模型概览\n",
"\n",
"本章中,我们首先介绍YouTube的视频推荐系统[[7](#参考文献)],然后介绍我们实现的融合推荐模型。\n",
"\n",
"### YouTube的深度神经网络推荐系统\n",
"\n",
"YouTube是世界上最大的视频上传、分享和发现网站,YouTube推荐系统为超过10亿用户从不断增长的视频库中推荐个性化的内容。整个系统由两个神经网络组成:候选生成网络和排序网络。候选生成网络从百万量级的视频库中生成上百个候选,排序网络对候选进行打分排序,输出排名最高的数十个结果。系统结构如图1所示:\n",
"\n",
"<p align=\"center\">\n",
"<img src=\"image/YouTube_Overview.png\" width=\"70%\" ><br/>\n",
"图1. YouTube 推荐系统结构\n",
"</p>\n",
"\n",
"#### 候选生成网络(Candidate Generation Network)\n",
"\n",
"候选生成网络将推荐问题建模为一个类别数极大的多类分类问题:对于一个Youtube用户,使用其观看历史(视频ID)、搜索词记录(search tokens)、人口学信息(如地理位置、用户登录设备)、二值特征(如性别,是否登录)和连续特征(如用户年龄)等,对视频库中所有视频进行多分类,得到每一类别的分类结果(即每一个视频的推荐概率),最终输出概率较高的几百个视频。\n",
"\n",
"首先,将观看历史及搜索词记录这类历史信息,映射为向量后取平均值得到定长表示;同时,输入人口学特征以优化新用户的推荐效果,并将二值特征和连续特征归一化处理到[0, 1]范围。接下来,将所有特征表示拼接为一个向量,并输入给非线形多层感知器(MLP,详见[识别数字](https://github.com/PaddlePaddle/book/blob/develop/recognize_digits/README.md)教程)处理。最后,训练时将MLP的输出给softmax做分类,预测时计算用户的综合特征(MLP的输出)与所有视频的相似度,取得分最高的$k$个作为候选生成网络的筛选结果。图2显示了候选生成网络结构。\n",
"\n",
"<p align=\"center\">\n",
"<img src=\"image/Deep_candidate_generation_model_architecture.png\" width=\"70%\" ><br/>\n",
"图2. 候选生成网络结构\n",
"</p>\n",
"\n",
"对于一个用户$U$,预测此刻用户要观看的视频$\\omega$为视频$i$的概率公式为:\n",
"\n",
"$$P(\\omega=i|u)=\\frac{e^{v_{i}u}}{\\sum_{j \\in V}e^{v_{j}u}}$$\n",
"\n",
"其中$u$为用户$U$的特征表示,$V$为视频库集合,$v_i$为视频库中第$i$个视频的特征表示。$u$和$v_i$为长度相等的向量,两者点积可以通过全连接层实现。\n",
"\n",
"考虑到softmax分类的类别数非常多,为了保证一定的计算效率:1)训练阶段,使用负样本类别采样将实际计算的类别数缩小至数千;2)推荐(预测)阶段,忽略softmax的归一化计算(不影响结果),将类别打分问题简化为点积(dot product)空间中的最近邻(nearest neighbor)搜索问题,取与$u$最近的$k$个视频作为生成的候选。\n",
"\n",
"#### 排序网络(Ranking Network)\n",
"排序网络的结构类似于候选生成网络,但是它的目标是对候选进行更细致的打分排序。和传统广告排序中的特征抽取方法类似,这里也构造了大量的用于视频排序的相关特征(如视频 ID、上次观看时间等)。这些特征的处理方式和候选生成网络类似,不同之处是排序网络的顶部是一个加权逻辑回归(weighted logistic regression),它对所有候选视频进行打分,从高到底排序后将分数较高的一些视频返回给用户。\n",
"\n",
"### 融合推荐模型\n",
"\n",
"在下文的电影推荐系统中:\n",
"\n",
"1. 首先,使用用户特征和电影特征作为神经网络的输入,其中:\n",
"\n",
" - 用户特征融合了四个属性信息,分别是用户ID、性别、职业和年龄。\n",
"\n",
" - 电影特征融合了三个属性信息,分别是电影ID、电影类型ID和电影名称。\n",
"\n",
"2. 对用户特征,将用户ID映射为维度大小为256的向量表示,输入全连接层,并对其他三个属性也做类似的处理。然后将四个属性的特征表示分别全连接并相加。\n",
"\n",
"3. 对电影特征,将电影ID以类似用户ID的方式进行处理,电影类型ID以向量的形式直接输入全连接层,电影名称用文本卷积神经网络(详见[第5章](https://github.com/PaddlePaddle/book/blob/develop/understand_sentiment/README.md))得到其定长向量表示。然后将三个属性的特征表示分别全连接并相加。\n",
"\n",
"4. 得到用户和电影的向量表示后,计算二者的余弦相似度作为推荐系统的打分。最后,用该相似度打分和用户真实打分的差异的平方作为该回归模型的损失函数。\n",
"\n",
"<p align=\"center\">\n",
"\n",
"<img src=\"image/rec_regression_network.png\" width=\"90%\" ><br/>\n",
"图3. 融合推荐模型 \n",
"</p> "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# 个性化推荐\n",
"\n",
"本教程源代码目录在[book/recommender_system](https://github.com/PaddlePaddle/book/tree/develop/recommender_system), 初次使用请参考PaddlePaddle[安装教程](http://www.paddlepaddle.org/doc_cn/build_and_install/index.html)。\n",
"\n",
"## 背景介绍\n",
"\n",
"在网络技术不断发展和电子商务规模不断扩大的背景下,商品数量和种类快速增长,用户需要花费大量时间才能找到自己想买的商品,这就是信息超载问题。为了解决这个难题,推荐系统(Recommender System)应运而生。\n",
"\n",
"个性化推荐系统是信息过滤系统(Information Filtering System)的子集,它可以用在很多领域,如电影、音乐、电商和 Feed 流推荐等。推荐系统通过分析、挖掘用户行为,发现用户的个性化需求与兴趣特点,将用户可能感兴趣的信息或商品推荐给用户。与搜索引擎不同,推荐系统不需要用户准确地描述出自己的需求,而是根据分析历史行为建模,主动提供满足用户兴趣和需求的信息。\n",
"\n",
"传统的推荐系统方法主要有:\n",
"\n",
"- 协同过滤推荐(Collaborative Filtering Recommendation):该方法收集分析用户历史行为、活动、偏好,计算一个用户与其他用户的相似度,利用目标用户的相似用户对商品评价的加权评价值,来预测目标用户对特定商品的喜好程度。优点是可以给用户推荐未浏览过的新产品;缺点是对于没有任何行为的新用户存在冷启动的问题,同时也存在用户与商品之间的交互数据不够多造成的稀疏问题,会导致模型难以找到相近用户。\n",
"- 基于内容过滤推荐[[1](#参考文献)](Content-based Filtering Recommendation):该方法利用商品的内容描述,抽象出有意义的特征,通过计算用户的兴趣和商品描述之间的相似度,来给用户做推荐。优点是简单直接,不需要依据其他用户对商品的评价,而是通过商品属性进行商品相似度度量,从而推荐给用户所感兴趣商品的相似商品;缺点是对于没有任何行为的新用户同样存在冷启动的问题。\n",
"- 组合推荐[[2](#参考文献)](Hybrid Recommendation):运用不同的输入和技术共同进行推荐,以弥补各自推荐技术的缺点。\n",
"\n",
"其中协同过滤是应用最广泛的技术之一,它又可以分为多个子类:基于用户 (User-Based)的推荐[[3](#参考文献)] 、基于物品(Item-Based)的推荐[[4](#参考文献)]、基于社交网络关系(Social-Based)的推荐[[5](#参考文献)]、基于模型(Model-based)的推荐等。1994年明尼苏达大学推出的GroupLens系统[[3](#参考文献)]一般被认为是推荐系统成为一个相对独立的研究方向的标志。该系统首次提出了基于协同过滤来完成推荐任务的思想,此后,基于该模型的协同过滤推荐引领了推荐系统十几年的发展方向。\n",
"\n",
"深度学习具有优秀的自动提取特征的能力,能够学习多层次的抽象特征表示,并对异质或跨域的内容信息进行学习,可以一定程度上处理推荐系统冷启动问题[[6](#参考文献)]。本教程主要介绍个性化推荐的深度学习模型,以及如何使用PaddlePaddle实现模型。\n",
"\n",
"## 效果展示\n",
"\n",
"我们使用包含用户信息、电影信息与电影评分的数据集作为个性化推荐的应用场景。当我们训练好模型后,只需要输入对应的用户ID和电影ID,就可以得出一个匹配的分数(范围[1,5],分数越高视为兴趣越大),然后根据所有电影的推荐得分排序,推荐给用户可能感兴趣的电影。\n",
"\n",
"```\n",
"Input movie_id: 1962\n",
"Input user_id: 1\n",
"Prediction Score is 4.25\n",
"```\n",
"\n",
"## 模型概览\n",
"\n",
"本章中,我们首先介绍YouTube的视频推荐系统[[7](#参考文献)],然后介绍我们实现的融合推荐模型。\n",
"\n",
"### YouTube的深度神经网络推荐系统\n",
"\n",
"YouTube是世界上最大的视频上传、分享和发现网站,YouTube推荐系统为超过10亿用户从不断增长的视频库中推荐个性化的内容。整个系统由两个神经网络组成:候选生成网络和排序网络。候选生成网络从百万量级的视频库中生成上百个候选,排序网络对候选进行打分排序,输出排名最高的数十个结果。系统结构如图1所示:\n",
"\n",
"<p align=\"center\">\n",
"<img src=\"image/YouTube_Overview.png\" width=\"70%\" ><br/>\n",
"图1. YouTube 推荐系统结构\n",
"</p>\n",
"\n",
"#### 候选生成网络(Candidate Generation Network)\n",
"\n",
"候选生成网络将推荐问题建模为一个类别数极大的多类分类问题:对于一个Youtube用户,使用其观看历史(视频ID)、搜索词记录(search tokens)、人口学信息(如地理位置、用户登录设备)、二值特征(如性别,是否登录)和连续特征(如用户年龄)等,对视频库中所有视频进行多分类,得到每一类别的分类结果(即每一个视频的推荐概率),最终输出概率较高的几百个视频。\n",
"\n",
"首先,将观看历史及搜索词记录这类历史信息,映射为向量后取平均值得到定长表示;同时,输入人口学特征以优化新用户的推荐效果,并将二值特征和连续特征归一化处理到[0, 1]范围。接下来,将所有特征表示拼接为一个向量,并输入给非线形多层感知器(MLP,详见[识别数字](https://github.com/PaddlePaddle/book/blob/develop/recognize_digits/README.md)教程)处理。最后,训练时将MLP的输出给softmax做分类,预测时计算用户的综合特征(MLP的输出)与所有视频的相似度,取得分最高的$k$个作为候选生成网络的筛选结果。图2显示了候选生成网络结构。\n",
"\n",
"<p align=\"center\">\n",
"<img src=\"image/Deep_candidate_generation_model_architecture.png\" width=\"70%\" ><br/>\n",
"图2. 候选生成网络结构\n",
"</p>\n",
"\n",
"对于一个用户$U$,预测此刻用户要观看的视频$\\omega$为视频$i$的概率公式为:\n",
"\n",
"$$P(\\omega=i|u)=\\frac{e^{v_{i}u}}{\\sum_{j \\in V}e^{v_{j}u}}$$\n",
"\n",
"其中$u$为用户$U$的特征表示,$V$为视频库集合,$v_i$为视频库中第$i$个视频的特征表示。$u$和$v_i$为长度相等的向量,两者点积可以通过全连接层实现。\n",
"\n",
"考虑到softmax分类的类别数非常多,为了保证一定的计算效率:1)训练阶段,使用负样本类别采样将实际计算的类别数缩小至数千;2)推荐(预测)阶段,忽略softmax的归一化计算(不影响结果),将类别打分问题简化为点积(dot product)空间中的最近邻(nearest neighbor)搜索问题,取与$u$最近的$k$个视频作为生成的候选。\n",
"\n",
"#### 排序网络(Ranking Network)\n",
"排序网络的结构类似于候选生成网络,但是它的目标是对候选进行更细致的打分排序。和传统广告排序中的特征抽取方法类似,这里也构造了大量的用于视频排序的相关特征(如视频 ID、上次观看时间等)。这些特征的处理方式和候选生成网络类似,不同之处是排序网络的顶部是一个加权逻辑回归(weighted logistic regression),它对所有候选视频进行打分,从高到底排序后将分数较高的一些视频返回给用户。\n",
"\n",
"### 融合推荐模型\n",
"\n",
"在下文的电影推荐系统中:\n",
"\n",
"1. 首先,使用用户特征和电影特征作为神经网络的输入,其中:\n",
"\n",
" - 用户特征融合了四个属性信息,分别是用户ID、性别、职业和年龄。\n",
"\n",
" - 电影特征融合了三个属性信息,分别是电影ID、电影类型ID和电影名称。\n",
"\n",
"2. 对用户特征,将用户ID映射为维度大小为256的向量表示,输入全连接层,并对其他三个属性也做类似的处理。然后将四个属性的特征表示分别全连接并相加。\n",
"\n",
"3. 对电影特征,将电影ID以类似用户ID的方式进行处理,电影类型ID以向量的形式直接输入全连接层,电影名称用文本卷积神经网络(详见[第5章](https://github.com/PaddlePaddle/book/blob/develop/understand_sentiment/README.md))得到其定长向量表示。然后将三个属性的特征表示分别全连接并相加。\n",
"\n",
"4. 得到用户和电影的向量表示后,计算二者的余弦相似度作为推荐系统的打分。最后,用该相似度打分和用户真实打分的差异的平方作为该回归模型的损失函数。\n",
"\n",
"<p align=\"center\">\n",
"\n",
"<img src=\"image/rec_regression_network.png\" width=\"90%\" ><br/>\n",
"图3. 融合推荐模型 \n",
"</p> "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 数据准备\n",
"\n",
"### 数据介绍与下载\n",
"\n",
"我们以 [MovieLens 百万数据集(ml-1m)](http://files.grouplens.org/datasets/movielens/ml-1m.zip)为例进行介绍。ml-1m 数据集包含了 6,000 位用户对 4,000 部电影的 1,000,000 条评价(评分范围 1~5 分,均为整数),由 GroupLens Research 实验室搜集整理。\n",
"\n",
"Paddle在API中提供了自动加载数据的模块。数据模块为 `paddle.dataset.movielens`"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import paddle.v2 as paddle\n",
"paddle.init(use_gpu=False)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Run this block to show dataset's documentation\n",
"# help(paddle.dataset.movielens)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"在原始数据中包含电影的特征数据,用户的特征数据,和用户对电影的评分。\n",
"\n",
"例如,其中某一个电影特征为:"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"<MovieInfo id(1), title(Toy Story ), categories(['Animation', \"Children's\", 'Comedy'])>\n"
]
}
],
"source": [
"movie_info = paddle.dataset.movielens.movie_info()\n",
"print movie_info.values()[0]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"这表示,电影的id是1,标题是《Toy Story》,该电影被分为到三个类别中。这三个类别是动画,儿童,喜剧。"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"<UserInfo id(1), gender(F), age(1), job(10)>\n"
]
}
],
"source": [
"user_info = paddle.dataset.movielens.user_info()\n",
"print user_info.values()[0]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"这表示,该用户ID是1,女性,年龄比18岁还年轻。职业ID是10。\n",
"\n",
"\n",
"其中,年龄使用下列分布\n",
"* 1: \"Under 18\"\n",
"* 18: \"18-24\"\n",
"* 25: \"25-34\"\n",
"* 35: \"35-44\"\n",
"* 45: \"45-49\"\n",
"* 50: \"50-55\"\n",
"* 56: \"56+\"\n",
"\n",
"职业是从下面几种选项里面选则得出:\n",
"* 0: \"other\" or not specified\n",
"* 1: \"academic/educator\"\n",
"* 2: \"artist\"\n",
"* 3: \"clerical/admin\"\n",
"* 4: \"college/grad student\"\n",
"* 5: \"customer service\"\n",
"* 6: \"doctor/health care\"\n",
"* 7: \"executive/managerial\"\n",
"* 8: \"farmer\"\n",
"* 9: \"homemaker\"\n",
"* 10: \"K-12 student\"\n",
"* 11: \"lawyer\"\n",
"* 12: \"programmer\"\n",
"* 13: \"retired\"\n",
"* 14: \"sales/marketing\"\n",
"* 15: \"scientist\"\n",
"* 16: \"self-employed\"\n",
"* 17: \"technician/engineer\"\n",
"* 18: \"tradesman/craftsman\"\n",
"* 19: \"unemployed\"\n",
"* 20: \"writer\""
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"而对于每一条训练/测试数据,均为 <用户特征> + <电影特征> + 评分。\n",
"\n",
"例如,我们获得第一条训练数据:"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"User <UserInfo id(1), gender(F), age(1), job(10)> rates Movie <MovieInfo id(1193), title(One Flew Over the Cuckoo's Nest ), categories(['Drama'])> with Score [5.0]\n"
]
}
],
"source": [
"train_set_creator = paddle.dataset.movielens.train()\n",
"train_sample = next(train_set_creator())\n",
"uid = train_sample[0]\n",
"mov_id = train_sample[len(user_info[uid].value())]\n",
"print \"User %s rates Movie %s with Score %s\"%(user_info[uid], movie_info[mov_id], train_sample[-1])"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"即用户1对电影1193的评价为5分。"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 模型配置说明\n",
"\n",
"下面我们开始根据输入数据的形式配置模型。"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"uid = paddle.layer.data(\n",
" name='user_id',\n",
" type=paddle.data_type.integer_value(\n",
" paddle.dataset.movielens.max_user_id() + 1))\n",
"usr_emb = paddle.layer.embedding(input=uid, size=32)\n",
"\n",
"usr_gender_id = paddle.layer.data(\n",
" name='gender_id', type=paddle.data_type.integer_value(2))\n",
"usr_gender_emb = paddle.layer.embedding(input=usr_gender_id, size=16)\n",
"\n",
"usr_age_id = paddle.layer.data(\n",
" name='age_id',\n",
" type=paddle.data_type.integer_value(\n",
" len(paddle.dataset.movielens.age_table)))\n",
"usr_age_emb = paddle.layer.embedding(input=usr_age_id, size=16)\n",
"\n",
"usr_job_id = paddle.layer.data(\n",
" name='job_id',\n",
" type=paddle.data_type.integer_value(paddle.dataset.movielens.max_job_id(\n",
" ) + 1))\n",
"usr_job_emb = paddle.layer.embedding(input=usr_job_id, size=16)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"如上述代码所示,对于每个用户,我们输入4维特征。其中包括`user_id`,`gender_id`,`age_id`,`job_id`。这几维特征均是简单的整数值。为了后续神经网络处理这些特征方便,我们借鉴NLP中的语言模型,将这几维离散的整数值,变换成embedding取出。分别形成`usr_emb`, `usr_gender_emb`, `usr_age_emb`, `usr_job_emb`。"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"usr_combined_features = paddle.layer.fc(\n",
" input=[usr_emb, usr_gender_emb, usr_age_emb, usr_job_emb],\n",
" size=200,\n",
" act=paddle.activation.Tanh())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"然后,我们对于所有的用户特征,均输入到一个全连接层(fc)中。将所有特征融合为一个200维度的特征。"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"进而,我们对每一个电影特征做类似的变换,网络配置为:"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"mov_id = paddle.layer.data(\n",
" name='movie_id',\n",
" type=paddle.data_type.integer_value(\n",
" paddle.dataset.movielens.max_movie_id() + 1))\n",
"mov_emb = paddle.layer.embedding(input=mov_id, size=32)\n",
"\n",
"mov_categories = paddle.layer.data(\n",
" name='category_id',\n",
" type=paddle.data_type.sparse_binary_vector(\n",
" len(paddle.dataset.movielens.movie_categories())))\n",
"\n",
"mov_categories_hidden = paddle.layer.fc(input=mov_categories, size=32)\n",
"\n",
"\n",
"movie_title_dict = paddle.dataset.movielens.get_movie_title_dict()\n",
"mov_title_id = paddle.layer.data(\n",
" name='movie_title',\n",
" type=paddle.data_type.integer_value_sequence(len(movie_title_dict)))\n",
"mov_title_emb = paddle.layer.embedding(input=mov_title_id, size=32)\n",
"mov_title_conv = paddle.networks.sequence_conv_pool(\n",
" input=mov_title_emb, hidden_size=32, context_len=3)\n",
"\n",
"mov_combined_features = paddle.layer.fc(\n",
" input=[mov_emb, mov_categories_hidden, mov_title_conv],\n",
" size=200,\n",
" act=paddle.activation.Tanh())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"电影ID和电影类型分别映射到其对应的特征隐层。对于电影标题名称(title),一个ID序列表示的词语序列,在输入卷积层后,将得到每个时间窗口的特征(序列特征),然后通过在时间维度降采样得到固定维度的特征,整个过程在text_conv_pool实现。\n",
"\n",
"最后再将电影的特征融合进`mov_combined_features`中。"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"inference = paddle.layer.cos_sim(a=usr_combined_features, b=mov_combined_features, size=1, scale=5)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"进而,我们使用余弦相似度计算用户特征与电影特征的相似性。并将这个相似性拟合(回归)到用户评分上。"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"cost = paddle.layer.regression_cost(\n",
" input=inference,\n",
" label=paddle.layer.data(\n",
" name='score', type=paddle.data_type.dense_vector(1)))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"至此,我们的优化目标就是这个网络配置中的`cost`了。"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 训练模型\n",
"\n",
"### 定义参数\n",
"神经网络的模型,我们可以简单的理解为网络拓朴结构+参数。之前一节,我们定义出了优化目标`cost`。这个`cost`即为网络模型的拓扑结构。我们开始训练模型,需要先定义出参数。定义方法为:"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[INFO 2017-03-02 17:44:56,684 networks.py:1472] The input order is [user_id, gender_id, age_id, job_id, movie_id, category_id, movie_title, score]\n",
"[INFO 2017-03-02 17:44:56,685 networks.py:1478] The output order is [__regression_cost_0__]\n"
]
}
],
"source": [
"parameters = paddle.parameters.create(cost)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"`parameters`是模型的所有参数集合。他是一个python的dict。我们可以查看到这个网络中的所有参数名称。因为之前定义模型的时候,我们没有指定参数名称,这里参数名称是自动生成的。当然,我们也可以指定每一个参数名称,方便日后维护。"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[u'___fc_layer_2__.wbias', u'___fc_layer_2__.w2', u'___embedding_layer_3__.w0', u'___embedding_layer_5__.w0', u'___embedding_layer_2__.w0', u'___embedding_layer_1__.w0', u'___fc_layer_1__.wbias', u'___fc_layer_0__.wbias', u'___fc_layer_1__.w0', u'___fc_layer_0__.w2', u'___fc_layer_0__.w3', u'___fc_layer_0__.w0', u'___fc_layer_0__.w1', u'___fc_layer_2__.w1', u'___fc_layer_2__.w0', u'___embedding_layer_4__.w0', u'___sequence_conv_pool_0___conv_fc.w0', u'___embedding_layer_0__.w0', u'___sequence_conv_pool_0___conv_fc.wbias']\n"
]
}
],
"source": [
"print parameters.keys()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### 构造训练(trainer)\n",
"\n",
"下面,我们根据网络拓扑结构和模型参数来构造出一个本地训练(trainer)。在构造本地训练的时候,我们还需要指定这个训练的优化方法。这里我们使用Adam来作为优化算法。"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[INFO 2017-03-02 17:44:56,753 networks.py:1472] The input order is [user_id, gender_id, age_id, job_id, movie_id, category_id, movie_title, score]\n",
"[INFO 2017-03-02 17:44:56,763 networks.py:1478] The output order is [__regression_cost_0__]\n"
]
}
],
"source": [
"trainer = paddle.trainer.SGD(cost=cost, parameters=parameters, \n",
" update_equation=paddle.optimizer.Adam(learning_rate=1e-4))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### 训练\n",
"\n",
"下面我们开始训练过程。\n",
"\n",
"我们直接使用Paddle提供的数据集读取程序。`paddle.dataset.movielens.train()`和`paddle.dataset.movielens.test()`分别做训练和预测数据集。并且通过`reader_dict`来指定每一个数据和data_layer的对应关系。\n",
"\n",
"例如,这里的reader_dict表示的是,对于数据层 `user_id`,使用了reader中每一条数据的第0个元素。`gender_id`数据层使用了第1个元素。以此类推。\n",
"\n",
"训练过程是完全自动的。我们可以使用event_handler来观察训练过程,或进行测试等。这里我们在event_handler里面绘制了训练误差曲线和测试误差曲线。并且保存了模型。"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD8CAYAAABn919SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXd4HNXVh9+7TdKqS5Z7kbvcsC0E2EAA04vpJSS0QIAk\n8CUhBTBfQuihpJCQEAihJYRAqOEL1XSbZmMbG4x7t1xlWb1ume+P2Ts7Ozu72pVWZeX7Po8fa2dn\nZ+7szvzuueeee47QNA2FQqFQpD+O3m6AQqFQKFKDEnSFQqHoJyhBVygUin6CEnSFQqHoJyhBVygU\nin6CEnSFQqHoJyhBVygUin6CEnSFQqHoJyhBVygUin6CqydPNmDAAK20tLQnT6lQKBRpz9KlS/dp\nmlbS0X49KuilpaUsWbKkJ0+pUCgUaY8QYmsi+ymXi0KhUPQTlKArFApFP0EJukKhUPQTetSHbofP\n56OyspLW1tbebkq/ITMzk+HDh+N2u3u7KQqFogfpdUGvrKwkNzeX0tJShBC93Zy0R9M0qqurqays\nZPTo0b3dHIVC0YP0usultbWV4uJiJeYpQghBcXGxGvEoFAcgvS7ogBLzFKO+T4XiwKRPCHpH1DS3\nU93Y1tvNUCgUij5NWgh6bbOP/U3t3XLs6upqZsyYwYwZMxg8eDDDhg0zXre3J3bOyy+/nLVr1yZ1\n3tdee42DDz6YKVOmMGPGDG688cak275s2TLefPPNpD+nUCj6J70+KZoI3elAKC4uZvny5QDceuut\n5OTk8POf/zxiH03T0DQNh8O+/3viiSeSOueKFSu47rrreO2115gwYQKBQIBHHnkk6bYvW7aMlStX\ncvLJJyf9WYVC0f9ICwsdQOvh823YsIHJkydz0UUXMWXKFHbt2sXVV19NRUUFU6ZM4fbbbzf2PfLI\nI1m+fDl+v5+CggLmzZvH9OnTmT17Nnv37o069r333svNN9/MhAkTAHA6nfzgBz8AYPPmzcyZM4eD\nDjqIE044gcrKSgCeffZZpk6dyvTp05kzZw4tLS3cfvvtPP3008yYMYMXXnihB74VhULRl+lTFvpt\n//2aVTvro7a3+QMENchyO5M+5uShedxy+pROtWfNmjX84x//oKKiAoB77rmHoqIi/H4/c+bM4bzz\nzmPy5MkRn6mrq+Poo4/mnnvu4ac//SmPP/448+bNi9hn5cqV/OIXv7A95zXXXMOVV17JRRddxCOP\nPMJ1113HCy+8wG233cYHH3zAoEGDqK2tJSsri1/96lesXLmSP/zhD526PoVC0b9IGwu9x010YOzY\nsYaYAzzzzDOUl5dTXl7O6tWrWbVqVdRnsrKyOOWUUwA4+OCD2bJlS1LnXLRoERdeeCEAl156KQsX\nLgTgiCOO4NJLL+XRRx8lGAx28ooUCkV/pk9Z6LEs6W3VzbT4AkwcnNuj7cnOzjb+Xr9+PX/84x9Z\nvHgxBQUFXHzxxbax3h6Px/jb6XTi9/uj9pkyZQpLly5lypTERw5/+9vfWLRoEa+++irl5eV88cUX\nSV6NQqHo76SHhS5A6w0T3UR9fT25ubnk5eWxa9cu3nrrrU4f64YbbuCOO+5gw4YNAAQCAR5++GEA\nZs2axXPPPQfAP//5T4466igANm3axKxZs7jjjjsoLCxkx44d5Obm0tDQ0MUrUygU/YU+ZaHHQkCv\nuFzMlJeXM3nyZMrKyhg1ahRHHHFEp481c+ZMfve733HBBRcYVv6ZZ54JwIMPPsgVV1zB3XffzaBB\ng4wImp/85Cds3rwZTdM48cQTmTp1KoMGDeI3v/kNM2fO5Be/+AXnnXde1y9UoVCkLULTek4pKyoq\nNGuBi9WrVzNp0qS4n9u+v5nGNj+ThuR1Z/P6FYl8rwqFIj0QQizVNK2io/3SwuWiFrIrFApFx6SF\noOs+dIVCoVDEIy0EvS/40BUKhaKvkxaCjhC9HuWiUCgUfZ20EHTlQ1coFIqOSQtBB5TLRaFQKDog\nLQRd0H16nor0uQCPP/44u3fvtn1P0zTuu+8+Jk6cyIwZMzjkkEN4+umnk27rSy+9xJo1a5L+nEKh\nODBIi4VFMsqlurGN7AwXmZ1I0hWLRNLnJsLjjz9OeXk5gwcPjnrvwQcf5P3332fJkiXk5uZSV1fH\nK6+8kvQ5XnrpJRwOB2VlZUl/VqFQ9H/SxkIH2FHbwro9PbfU/e9//zuHHnooM2bM4JprriEYDOL3\n+7nkkkuYNm0aU6dO5YEHHuDf//43y5cv55vf/KatZf/rX/+ahx9+mNxcPRdNfn4+l156KQDz589n\nxowZTJs2jauuusr47PXXX8/kyZM56KCDuPHGG1m4cCGvv/46P/nJT5gxY0bSSb8UCkX/p29Z6G/M\ng91fRW0uCgTJ8ZsyDGYk0ezB0+CUe5JuysqVK3n55Zf55JNPcLlcXH311Tz77LOMHTuWffv28dVX\nejtra2spKCjgT3/6E3/+85+ZMWNGxHH279+Pz+dj1KhRUedobm7miiuu4MMPP2Ts2LFGytzzzz+f\n119/na+//hohhHGOU089lfPOO4+zzjor6etRKBT9n7Sw0HuDd955h88//5yKigpmzJjBhx9+yMaN\nGxk3bhxr167lRz/6EW+99Rb5+fmdPsfq1auZMGECY8eOBfR0uQsWLKCoqAiHw8FVV13Fyy+/HJH1\nUaFQKGLRtyz0GJZ0e9V2Wtvb2akVA3DQ8IJub4qmaVxxxRXccccdUe99+eWXvPHGGzz44IO8+OKL\nccvHFRUV4Xa72bZtGyNHjkzo3G63myVLlvD222/z/PPP89BDDzF//vxOX4tCoTgwSAsL3an5KaQR\n0YOxi8cffzzPPfcc+/btA/RomG3btlFVVYWmaZx//vncfvvtLFu2DCBuKtt58+ZxzTXXGO/X19fz\n1FNPMWnSJNavX8+mTZsAPV3u0UcfTUNDA/X19cydO5f777/fyH2u0uUqFIp49C0LPQZtrlyy/LXk\naK00kNUj55w2bRq33HILxx9/PMFgELfbzcMPP4zT6eS73/0umqYhhODee+8F4PLLL+fKK68kKyuL\nxYsXRxS6+OEPf0hTUxMHH3wwHo8Ht9vNDTfcgNfr5bHHHuOcc84hEAhw2GGHcdVVV7F3717OOecc\n2traCAaD/P73vwfgW9/6Ft/73vf43e9+x3/+8x9KS0t75LtQKBTpQVqkz62qb6GoYR21ZLNDG9Aj\nLpd0R6XPVSj6D/0qfS7CQQNZ5NGs0gAoFApFDNJD0IF6LRu3COAluo6nQqFQKPqIoHfk9hECGsgi\nqAnyRHMPtSp96Uk3mkKh6Dt0KOhCiMeFEHuFECtN24qEEG8LIdaH/i/sbAMyMzOprq7uUIQCOGgk\nkzyaQAlWTDRNo7q6mszMzN5uikKh6GESiXJ5Evgz8A/TtnnAu5qm3SOEmBd6fWNnGjB8+HAqKyup\nqqqKuU9jq5/aFh+NopVCGqAmCE5PzP0PdDIzMxk+fHhvN0OhUPQwHQq6pmkLhBClls1nAseE/v47\n8AGdFHS3283o0aPj7vPowk3c+dpqBlDH4oxrcMyZB8fM68zpFAqFot/SWR/6IE3TdoX+3g0MSlF7\nbAmGXCz7yGepNh7WvNqdp1MoFIq0pMuTopru/I7p1BZCXC2EWCKEWBLPrRKPgCkv11uBQ/QEXjVb\nOnUshUKh6K90VtD3CCGGAIT+3xtrR03THtE0rULTtIqSkpJOnSxomgSdHwzF1q95rVPHUigUiv5K\nZwX9/4DLQn9fBiRfrSEJgsGwoG/TBsHAKbBauV0UCoXCTCJhi88AnwIThRCVQojvAvcAJwgh1gPH\nh153G0GrQ2fSXNj+GTR2zoWjUCgU/ZFEoly+FeOt41LclpgErXHnZXPhw3th3RtQfmlPNUOhUCj6\nNH1ipWhHRAn64GlQMFK5XRQKhcJEWgh6wOpzEUK30jd9AG0qP7hCoVBAmgh6lA8ddEEPtMGGd3q8\nPQqFQtEXSQtBt83zMnIWeAcot4tCoVCESAtBH1Hkjd7ocMLEU2D9fPC393yjFAqFoo+RFoJ+0WEx\niiuXzYW2etiyoGcbpFAoFH2QtBB0IQRlg3Oj3xhzDHhylNtFoVAoSBNBt2L41N2ZMO54WPs6BIPx\nP6RQKBT9nLQRdCHC1UR9AdMk6aTToXEPVH7eC61SKBSKvkPaCLqZiLj08SeAw61S6ioUigOetBR0\nn9m9kpkPo4/SBV2VplMoFAcwaSnogYBFuCfNhf2bYO/q3mmQQqFQ9AHSRtDdTpMP3ToBOvFUQCi3\ni0KhOKBJG0F/8Nvl5GboySH31rdFvpk7GIYfogRdoVAc0KSNoI8o8vKr0ycDMPdPH0XvMGku7FoB\ntdt6uGUKhULRN0gbQQdwO+M0t2yu/r8qTadQKA5Q0krQXSY/ehTFY6FkkhJ0hUJxwJJegu6II+ig\nu122fgxN1T3TIIVCoehDpJmgd9DcsrmgBfXSdAqFQnGAkVaC7ozncgEYMh3yR6hkXQqF4oAkrQTd\n3ZGFLgSUnQYb34O2xp5plEKhUPQR0krQnSYfuj8QI7uiLE238d0eapVCoVD0DdJK0CNWi1qX/0tG\nzoasIuV2USgUBxxpJeguUxx6uz+Ghe506aXp1r2lStMpFIoDivQSdJPLpT2WywVCpenqYKvNilKF\nQqHop6SXoDsTFPSxc8CdrdwuCoXigCKtBF1g8qHHcrkAuLNg3HH6qlFVmk6hUBwgpJWga4QnQuNa\n6KC7XRp3w85l3dwqhUKh6BuklaCbS8/FnBSVTDgRHC5Y/d9ubpVCoVD0DdJK0LPcTuPvDi30rEIo\n/YYqTadQKA4Y0krQx5Tk8N0jRwMJWOigJ+uq3gBVa7u5ZQqFQtH7pJWgA5w6bTAAbYkI+sRT9f9V\nJSOFQnEAkHaCnuXWy9C1tAc63jlvKAyrUIKuUCgOCNJO0LMzdD96c7s/sQ9Mmgs7v4C6ym5slUKh\nUPQ+XRJ0IcRPhBBfCyFWCiGeEUJkpqphsfB6dAu9KRELHVRpOoVCccDQaUEXQgwDfgRUaJo2FXAC\nF6aqYbHwekIWeluCFvqA8TBgonK7KBSKfk9XXS4uIEsI4QK8wM6uNyk+MnSxOVELHXS3y5aPoXl/\nN7VKoVAoep9OC7qmaTuA3wLbgF1AnaZp81PVsFg4HAKvx5m4Dx1CpekCsO7N7muYQqFQ9DJdcbkU\nAmcCo4GhQLYQ4mKb/a4WQiwRQiypqqrqfEtNeD2uxH3oAENnQt4wlaxLoVD0a7ricjke2KxpWpWm\naT7gJeBw606apj2iaVqFpmkVJSUlXThdGK/HmbgPHSJL07U3p6QNCoVC0dfoiqBvA2YJIbxCCAEc\nB6xOTbPio7tckrDQQXe7+FtUaTqFQtFv6YoPfRHwArAM+Cp0rEdS1K64ZGe4khf0UUdAZoFyuygU\nin6Lqysf1jTtFuCWFLUlYbweJ43JuFwgXJpu7RsQ8IHT3T2NUygUil4i7VaKgvShJ2mhg+52aa2F\nrR+nvlEKhULRy6SloGd7XDQlE7YoGXssuLKU20WhUPRL0lLQvRnOxJJzWfF4VWk6hULRb0lLQe+0\nhQ6626VhJ+z6IrWNUigUil4mLQU9y+Ok1ReMKEmXMBNOAuFUbheFQtHvSEtBzw5lXExq+b/EWwSl\nR6pkXQqFot+RloLuzehEgi4zZXNh3zqoWpfCVikUCkXvkpaCLjMudmpiFPQ0AKCsdIVC0a9IS0F3\nOfVmb6xqpHTea3yyYV9yB8gfBkPLlaArFIp+RXoKukMAsHC9LuT/Wb4j+YNMmgs7lkJ9t6dwVygU\nih4hrQW91ae7XNzOTlyGKk2nUCj6Gekp6E5d0OWkqMfVicsomQjF45XbRaFQ9BvSU9AderNbQha6\npzMWOoRK030ELTWpappCoVD0Gmkq6ClwuYDudgn6Yd1bqWqaQqFQ9BppKehOR6TLpdOCPrQccoco\nt4tCoegXpKWgy7BFGYfudonOHcjh0GPSN7wLvpZUNU+hUCh6hfQU9JCFLn3obkcXLqNsLvia9Xqj\nCoVCkcakpaBLl0unV4qaKT0SMvNVsi6FQpH2pKWgS5/57vpWAHxdyW3udMOEk2HdGxDoZEpehUKh\n6AOkpaBLC11y35tr2bC3sfMHLJurhy5u+6SLLVMoFIreIy0F3eWIngT9srK28wccdxy4MpXbRaFQ\npDXpKejOTka1xMKTrdcbXfMaaJ0omqFQKBR9gPQU9K5EtcSibC7UV8Ku5ak/tkKhUPQA6SnoqbbQ\nASaeAsKh3C4KhSJtSU9Bt/GhdxlvEYw6Qq0aVSgUaUtaCro1ygVApELjy+ZC1RrYtyEFB1MoFIqe\nJS0FvdO5WzpClaZTKBRpTFoKup2FHujC2iKDghEwZIYSdIVCkZakpaDb+dDb/alQdHS3S+XnUL8r\nNcdTKBSKHiItBV3YOMx9KTHR0YteAKx9PTXHUygUih4iLQXdjpQJekkZFI1VbheFQpF29BtBb0+V\noAuhW+mbF0BLF9IJKBQKRQ/TbwTd50/hkn1Zmm79/NQdU6FQKLqZtBX0P397JvN/cpTxuqa5PXUH\nH1YBOYOV20WhUKQVXRJ0IUSBEOIFIcQaIcRqIcTsVDWsI+YeNJQJg3KN109+soXm9hTlM3c4oOxU\nWP+OKk2nUCjShq5a6H8E3tQ0rQyYDqzuepM6T12LL3UHKzsNfE2w6YPUHVOhUCi6kU4LuhAiHzgK\neAxA07R2TdN6dRaxqS2FFYdKj4KMfOV2USgUaUNXLPTRQBXwhBDiCyHEo0KIbOtOQoirhRBLhBBL\nqqqqunA6e+44a6rxd0NrCgXd5YEJJ8JaVZpOoVCkB10RdBdQDjykadpMoAmYZ91J07RHNE2r0DSt\noqSkpAuns+eSWaP45WmTAGhs87NyRx3feWJxalaOls2F5mrY/lnXj6VQKBTdTFcEvRKo1DRtUej1\nC+gC3+McMW4AAI2tfq5/4Us+WFvF+r0NCX/+qPve59yHbOqJjjsenBkqR7pCoUgLOi3omqbtBrYL\nISaGNh0HrEpJq5IkJ8MFQEObn2BQj0d3JJFPd9v+ZpZurYl+IyMHxs5RpekUCkVa0NUolx8CTwsh\nvgRmAL/uepOSJzdTF/TGVj9BLXlBj0vZXKjbBru/TM3xFAqFoptwdeXDmqYtBypS1JZOky0t9FY/\ngZCgd0XPg0GNFl9AP665NN2Q6alorkKhUHQLabtS1Izb6SA3w0VtS7vhcomVrKu5PeyWicVv5q9l\nyi1v6WGQ2QNg5GxY8yqtvgBt/kDK269QKBSpoF8IOkBBtpuapnakVgdsRDsQ1Jj8q7e45f++jnus\nF5ZWAnrUDKC7Xfau4qRfPclR972f0nYrFApFqug3gl7k9VDT7DN86L5AtKBLq/3pRVuTO3ioNN2J\njiXsqW/rWkMVCoWim+g3gl6Y7aGmOexy8du4XKSgx/K4yM/KgBYjsKVwFAyexknOJSlts0KhUKSS\n/iPo3pCgh0TYb6Pafhur3UyLT/rHbfzwZadTLtZTgsqRrlAo+ib9RtDzMl3Ut/hNLpfYFnosmtt1\nQdfs/PCT5uIQGsc7l6amwQqFQpFi+o2gu50O/IGgIej+gMY7q/ZQOu81qhp0v7evg+iWFinoodf+\noKkDGDiZLcFBnORQbheFQtE36T+C7nLgC2gml0uQvy3cBGCkAbDzq5uRZew0m4lVf1BjfrCCwx0r\nobUu1c1XKBSKLtN/BN0h8AWDhpvEF9BoDSXoynA5jW3xkC4ZuZfZ5dLqD/JWoAKPCMD6t1PceoVC\noeg6/UfQnQ40LWyF+4NB2kKTnB6nfpkd+dANQdciX4PujlmmjadK61850j/dWE3pvNfYvr+5t5uS\nUpZtq2HD3sakPhMIaqzZXd9NLVIoup9+I+iukGg3hfzgvoBGW8hC94V84R1FuUgLXrpcIix0XwAN\nB28HDkZb/zb4WlN7Ab3E80u3A7Bo8/5ebklqOecvn3D87z9M6jN/eX8DJ/9hISt3KJeaIj3pN4Lu\ndkYmb7nlla/ZvK8JCAu5L5ighW68jhR0gPnBCkR7I2xOTiz6OprKJsmXISGvrFF1ZBXpSZeSc/Ul\n3M7IvikcU667YV5ZvoPFHVihv3lrLYVej+FyMUe5tPr0vz8JTiHoycGx5lWYcFKKWt97CFKUlbIf\nII0Cfwcdv0LRV+m3gm7GF9T48bPLOzyGzIku86ubFyfJpFztuGktPR7vmtdh7h/A4exKs/sMyj4H\nl0O/hzpyzSkUfZV+43JxOWNbmr4ky9Fpplh2iVnc60tPhOZ9sH1R1GcPBDbva2Lm7fOprOn7E6l2\nSdpiIe+hjibPFYq+Sr8RdE8cCz3ZIbScWDXHrZtT7tYPPwacnk6XpqusaebtVXs69dlU05m88f/+\nfDs1zT5eWb4z9Q1KMfUtvoT3dUsLPYlOQKHoS/QbQZfWVaY7+pLs4s8TmQQ0P9gR7hdHDow5Rg9f\n7MRk4ul/+oir/tHHVpwmcRmy7+wor3xfoKndn/C+8h7qaAGaQtFX6TeCLn3oWe5on7Zc0m8mESvM\nbNkHTMLdHgjqOdJrt8KelUm3taZZtxrTNbJElvcLpEH725Jwt7mN9Qp9/7oUCjv6kaDrIuNxOSjO\n9kS8t7+5PWp/XyBITVN73ApEZh96wPS3LxDUS9MhOu12gb45tF+3p6HDfURI0Ptg86No9we5+NFF\n/PCZLzrc1+VQUS6K9KYfCbp+KR6Xg/d+fkzEe394Z13U/r6Axsw73ubap2M/6GbBjbDQ/UHIGQgj\nZ3Vp1WhXoik0TeOpz7ZS35q4j9gO6ULX0HhhaSUn3r+ABeuqIvbZXddK6bzXjCggZ0jQ02GE0eYP\n8tGGffx3Rcf+frdLWeiK9KbfCbrb6SA/yx3xnowhN9MUKi/3zurYk5OxJkWNKIiyubrLZf/mTrW5\nvQu+2g17G7n5Pyu59ullnT6GGU2DjzfsA2BvQ2RVpk836duf+nQLEPahJxNB0lu0J+NykRa6EnRF\nAvgDQZZs6VsrrPuRoIdcLnGiXczsb4p2w0iGF2YB8JlpIZLfVtD10nSseS2ZpoaP2QlBP/73H3Ls\n7z4w3B4L1+/r1LmtBDTNSDNcnBPpsgr7zPXX6eRySaaotzMU5dKqCoErEuD+d9Zx3sOfsmJ7/KI3\nH2/Yx+c9JPz9RtCl6HhCw+Y3r/sGT1x+SMz94wm616NPrL725S5jW9DkXjAm2opGw6CpnXa7+AIa\nG/Y28PCHGxP+zIa9jWyqaopoT1eQYYv+gMa+Rl3QHZZYxrCA6+d0hizZhz/cyIRfvJGSdnQXyVjo\nWijU56EPNqbF6EPRu3xZqaeKqLGZozNz0aOLOP/hT3uiSf1H0KXfU1roZYPzOHLcgJj7S0GXE2Fm\nsjzhBbTSGvdHTIqaHvayubDtM2iM9Dsn1uYg5z38Kfe8sYZPNu5LKgww1YLjD4Yt9IBlUlB+RdJn\nbv7K2gNBI89NZ7j7jdVM+GX3dQrJRLlErDVIIn5dcWAitSFRr0BP0Hda0kWkJSYtdLAXa0l1SNDt\nUgZ4TaGPtaEQQ/OkaMRKwklzAQ3Wvm57nsa22HHQvkDQOP63/7aIxz9O3BefckEPBI3vxOoJkha7\n1HmrBd+RhRKPv364iXZ/MCUTrI1tfj7bVB2xLRkL3fwb1/WCoNc1+zjh9x+ydnfHkUaK3kfeW25X\n35HRvtOSLjJ+UA4AF88aZWwTQhi+9aMmlPCdw0u577yDANjfpFuj1iyNABluB9+sGAFAbbMUucgo\nl1eW79B/0EFToWCkrdtl7e4Gpt7yFq9+aR9hYY2m2Fqd+FJ6q8tF07ROiaJMzhUR0RMjbC9oWOgW\nQW/qWPyCQY0Xl1bGnDeI1/EBLN26n7rm+Oe57tnlXPjIZ4brCOL70HfUtkR8Z+am9Yagf7RhH+v3\nNnL/29FRWYq+hwxq6Evp7fqNoA/Ky2TLPadx6rQhEdtlwqWRRVncesYUSnIyANjXoAu1xxW2xudM\nLAF0q33udP04chGQWdAXrq/ix88u59evr9ad0GWnw6YPoDWyOMLqXfrrN1butm2zNWeI9N0ngtVA\nP//hTxl9k/0oIRGqG8NWtjU+Xr5eUalP/jgdyVvoLy6r5GfPr+Cxj8KjELM1XRtHrP2BIOc+9CmX\nPrE47jnk921eSGa3qAxgy74mjrjnPS5/8nND1IO9bKHL37+lCy6sWHy8YR/VjW0d76hIGJ8/ulRl\nb9NvBD0Wcjm3N+QXly4WKU75Wfr2O86aytnlwwHdJ1bo1SM9pFiZH/ZPNupCZKTjnTQXAu2w4Z2I\nczeErE7zsP+9NeEwyQ/W7o3YPysJQbe6XJaEYsSTRV6X2d1jPbZMbranvo3Fm/dHjQTiTTBLpDun\n2rTvT/4dzoAZT0BlKcGOogkk5uY1tNpb/jI084O1VTy9aBsQed29IeiZIVdfrE6oswSDGhc9uoiL\nHu18MrlFm6r50TNfpMXag57CZ6qOFoue/r76vaAXhVaNypQA0sWyJuSnlL2rQ4R97gVeNwVePZa9\nJiRA5knR5tADt02WbRtxGHgHRLhdWn0Bbv7PytA5wj/4FU+Gc7g8+cmWiLZ6OvDFmW+OVEW52Pni\nowTd1P6qhjasBkki4iePaXbXmO38eMdIVODkoc2FTBpMrhzzpGeG6bv+v9Cio9620KXx0exLPP9M\nIsgR1pou+OYve2Ix/7diZ7eMHtIVOeEeb8V3T1vv/V7QC0KWtsxx7rJMgsoIDacQnDRlMNefNJGb\nTp1kstD1B9tOQA2/r8OppwJYNx/8uuXXbBKhWFEgVou8zWYBVMT7Jku/M5OiDa0+NuyNfKh9Nsex\n3qBmQc90O6KicayRJJqmcddrq1husqjtImTMxHPbxPr+5n+9m9J5r7G3Xi8HKAW93R80/jZb6OaF\nXOaORS4yCwY1o8PvyKffHUijoTnFFnoqJtDlIXo6nLOuxcfSrX1r8Y6kOjQPF28hWk+vaej3gn7E\n2GIAjp88CIgOMZJi4RACp0Nw7Zxx5GS48HqceJwOY1K0w7wrk06H9gbYvACAZlOWP+mftk4Iet2R\n9UU6+vEibgIjAAAgAElEQVTNgt+Z1YwXP7aY43+/IGKbeQJUipn1oW03ncshRNR3YZ149Ac1/rZw\nM2c9+LGRM11eutX/LonlGoFIQT/1jwu57b9fA/C3hZsA2BQqNWi01x803C7m1AjmjsdO3AOaRmZo\nTiXVGRdrmtr51Ssr407Syu+9NcWC3lHpxUSwqxHQE/zvy19x7kOfRkx09wU272syVqDHu1dS/Vt2\nRL8X9B8fP55P5h3L6AHZQHQhDPmjWPOCCyEo8LoNyzGWZSKt1825B9PmyKJ66Uv4A0GOvPd9Y5/d\nIQuy3iJa1nN2ZKGbBd889E3UTyd90Ob9zQ+oFLN4Lpc2fzBqtGJtt/nzR977Pos377eNkBGmvxvj\nCLr5WlftqueJj7cAYbeI9D3LiB2zcJvjyc1ian4IZQcVCOoRTtC1ofK6PQ1GmgTJb+ev5R+fbuU/\nX+yI+Tnpi21OsVsjkAIRlj+5XfGPNn+A8x76pFssaTn6WtbJOaLuosqUHiOesdfTLqp+L+gZLidD\nC7KM11YLUVpqdpZjUbYnKsrlzrOmRuwjJ9eueuZr3vZNJ7j6NW58Ppzwa9KQPBpa/TS3+6PcClar\n1DpJasVsqZpvlGTFJ2AImMZ8U6GNDHcMQfebBT0Q9b7V5WJ9f39TW8yQR0lDHBeHXS4egPqWyEln\neWjz92T+zr/eGY5CMj+EZpeL0yFwOUSXqhadeP8Cbn7l64htQUMQY/9W8nurbfaltGpSPMEJBDUu\n+OunLFwff2GcPIKdi27zviaWbK3hf1+Kn0p62bYa7ntzTYftNTOqWDfENlY1dbBnz2I2auJNiipB\n72ZGFHptt9sJTYHXbUyKyoftosNGRuyzv7GduhYfG/Y2Mj9QQYmoY8uKD03n0zuT3XWthvtGYs2U\nuKW6OW4ha7OwtZhcOuaHPxFr/Z431vDYR5ujziWLg8Tzobf5gzaCbnG5WETL7XQYFZpiLapL1EIH\nGFOiP+TSQree37y/OUb+8ic+N8Te7HKR/vKApuEUApcz2q3UGczfk1EUJM7vYz5nKhcXxROcmuZ2\nFm/e32HN3aDhcokT0dFBlZSLH13EXz7YaHSgiSDP2xNlAf/07npK5yWWl8n828brpGMZI91FlwVd\nCOEUQnwhhOh8HtkeJMvj5OyZw4DIlaQOGwu9JDeTPQ36kC+oaThEpJsAdDF44ys958v7wRm0a05O\ndIYjWYpDce9NbQGqGqItdGuFpXjL6BtMHYA58iNW7dNYPPrRZu54dZXhXpBkGhZ65E1o9qHbulws\nFrpVQL7eWW9EWIgIl0t4n8a2xKNcBoS+Uynchtss9L75O7T6XuWoyPydyc4kGNRwOARuhyNhAalu\nbIvZiZrDVQ0/fQfWssTaSbX6Atz+31UR90CixPN7J9IuMLtcku/o2v1Bnvh4sxFnvz2JWrTy3uoJ\nQf9daEFXIkaRuT3xvt+erk+bCgv9x8DqFBynx5BCbq5uZDdXN6wgi911rWzf38yf3ttgm13w759s\nMQSmAS+fBKdykmMJcpCaHbqJ2wNBdtS2RH1+aH4WZYNzjdeXPr6YF5ZW2rbbHErXYur5zdZmMjeQ\n9UaUoXzWQ5iFqc1n43KJ40OHyBWwMV0uJgv9sY828/yS7eHjR40AIs8n35edhbk91s6mrkXvVOX3\ndPr0obQHgrorSdNdLm5XpKDH6mRX7qjj4Dvf4WWTX9w8GW7+XaRLzxrr/vXOuvB1RaxGjvwOX1ha\nyeMfb+YP76y3bUs84nXysnOOl0fIfP121n5Hc65//2QLt/13FftCwQGV+6Ofg1jIzryjFA5t/gDb\n96emaHkiRpH5d7QaQLvqWtgcmqj/bGNkKorupkuCLoQYDpwGPJqa5vQMMvfCiKKw+8VpIzTDCrPw\nBTRu/b+vo96TvPn1bnbWhW/Q+cEKSh17eNZzJ/e7H+TknX/mKuer5K9/Cc+2BYwXlRTQgBT83EwX\nb153VMQx/7Zgk+25zIJuFg7zQ+bzJ25BWRNQZcWw0H2BoBH2mYjLxepn3V0f/n7MHWekhR6+njte\nXcX1L3xpvLZa6FYrsc1iocfzW97zxhpafQHjmIWh9QaNrX4CwZDLxSGMzm5HbQtlN7/JM4u3RR1L\nhmWaU6NO/tVbxt9mEQoXBQl//qJHP+O0Bz4yXpu/d2uufPldNSdRI9XuuFZeXqZ3RvHKCd735lrj\nb7v7qyMjwjqqqEzCQjcEvYNz/PS5FXzjvveTSpccC7Ohs6e+lWv/tSzKTeSP43KZffd7zPntB1TW\nNBtWf0/h6niXuPwBuAHIjbWDEOJq4GqAkSNHxtqtR5Ghi+eUD6N1cYBNVU1RrhSAYQWZQMcLMvbW\nh4f1rwZmcahjNcPFPsrFeobuXkqFuxU++hfjgEt0Y552zUkVBbTUFcMz4/iwLIf/bPCzVytgQv5Y\n2J7Fl7UZ5JUMIzs7h5LcjAhB32daqm9+yJIpmiGHvrmZrpD7JxSyZ+ND93qcNLX7owRdCJtJUcsN\n/vEGeytFmJYWNZjcHlasFrJVQKzhnvHcVu+s3kvZzW8ar2UxlOb2gO5WcwjcTgf7m9r53fy1VJQW\nAfDyFzv41qGR9698yLNDq5B317XGbKe00Besr+Kc8mEU52Swcoc+SesPBHE5HRFCIjuDuhYfb6/a\nY9yzyWSPlMSzOO96XR9cm91ou+pauPofS3nssgoG5mVGuEjsQiDlPRerT9hl+V6qkghBlBE/HVno\n74YK1bT7g2S4El9xbUd7IEgWTvyBILf992te/2o3R08o4YJQfieItNBjzVGY02n0FJ0WdCHEXGCv\npmlLhRDHxNpP07RHgEcAKioq+sS64bNmDuPJT7ZwzMQSPlxXxaaqphguF92Ct3OVmDG/X0821/n+\nx3h9xylTuPeVz3n03JHM/2wFntYq2mp3MVDUUiLqmOhshtrtDK/fxQ+d1TiEBtuAx+Cg0DHqNC+U\nDOd4Xx5Fbg9VWj5524fT6nBRRQHsHUgh9dSSw+/fXsczi7ex4a5TjEVU//h0i+0DIW+4/Cw3Da1+\n3E6BENGi2twewOtxkuFyGK4JSW6GK+rY8SbhYomLFOFqmzQC0r300Y1zuPuNNUbOFskvXl7Jtw8d\naZjoyUxE5WW6Q+cIGBa62ymYv2oP81ft4dxQOgi7yUBD0EOjF5lOQtLu11MLu50OY45m4fp9/OCf\ny3ju+7PD1+4PkuN0WCba9PNd//wK5q/aw7VzxhrHTJSqhjZOuP9D5p1cZmx7fsl2zikfHhXVZf7J\nnvp0K1/tqOO5Jdv5n2PHR8w12fmLfXHatGVfE89bXIhmA6gj5OR/R6MAaRwk8/3EQv7W3//nMqOi\nmTXPkj9C0GPMofSw/xy6ZqEfAZwhhDgVyATyhBD/1DTt4tQ0rfuYMaKALffo1YbkitAmm6HssMKs\nqG12xBP8nEwXjXhpyC5lifBRUOxhQXU4ROyMUUN54Fszqa5vZfav51NEPQNFLXOGa+zZuY2B1FIi\najk/30Nwz3ami+0MdNTirWnjAllY6N9380Um+DQn+1bkc6GngODTT0LBYJo9A1i3sIYqLZ9yUcBe\nCqjSCmjDY4inLmotuJwONA0eeG8D7QGNa+aMJS/TTWObn9xMNzXNPtp8kaluXU5H1DC3o4m/B95d\nz+8tQ1H5IC7aHG3NSxfK0PwsPE7dv22duDI/PMmEiuWFcvm0tAcIBPXJcfNqYukeWLatli+21TBz\nZKHxngy1lCkbNlsXOAWClN38JmfNGBpxL1knBdt8AXIyXBYfun49cg2DsW8SgrVgXRW1zT4eMbnw\nrn/hS/xBLWq0Ye6kjXTJoU3OCEG3iUOPI1x2xVusJQ7jIVfNxpuMvfv11cZvbv1+1u9pwJvhYlhB\nYs+y+Vzm8pRWQTe7sWJNiqY6J08idFrQNU27CbgJIGSh/zwdxNxKOGdLdPRAToaL/Cx3h3k9rENt\nM1mh1aD+QJCmNn9U2GRupv5+httJACdVFFKlFZLrLuKzwBhjvz9sdnNc2SA+3biPfY3tDM7y42ja\nQwl1lIjakMVfa7xurNqKf+NiSkQ9d7qjH7h6LYvm1cWc6clFtBSz1eVkwP5ipro0GrUsaj/y8sK2\nIVxx7HSG1m9hgDMHnPDSJ7VMGhXOaDmiMCvqIYr38PlDgm5FHuPNUGbKMaGFYKALXobLEXKH6P7t\nqMiagGYIkbT2s9zODsU912ShBzUNpwM0LSxgZkE9+y+fsOWe0zj0rncozsmgtNgbcb4tFkGXk3T/\nWb6THx033thekpsRkQpYJh+LSNEcEslwlSi900jGAjVGSpbRp10yNfO5HaYJ3Jqmdl41Ve7675c7\nqSgtisg7FMtCDwQ1nv18e9T2eIK+aFM1+V43ZYPzgMQmRf9q6rCs98VpD3xEeyDI2jtPTtgVYzca\nsE7mR0SWhfa/5LFFlJs6fOt1appm69pNJV31oac90kKPJdoDLb5rO6zx5GayM8JRLk1tAbIznKy4\n5URu++/XvLRsB3khH641fNF6A9c2+9jf1Ea+10NQg231QWAIWxiCbfhvaBBwQfkQ3l+2Whd7UcdA\nUWOI/kRnC072MjxQyQhnPYX1bRzjbNbdPgC7gKfhTvNxM4E90JyZiTengF11HuqDmfCPoQTcOTiy\n8hjoy+QnrhoatSwayaJRy6IBLw1aFjm1AYaxnwayaCKLAPr3I618ed1mUWlo8xt5b9whC93qJ69r\n8RmuIvl7ZWe4aPEFECK2f9fO5eJwhR+6XbXRnfXehjb2NrSxrVoXcNmWPRZrWk6alhZ7I1x6BV4P\n02+fb7yWn7ez0GW75bVZ7wtfIEibPzxpvXB9FbvqWjm3fDg3vvgVkHy+bln4WNM0XlwW6S55ZvF2\n8rM8zDsl7MYxfOiW48RyvcWLQ//mI58BGCPo5gQnRY22WL4f+bldta2UmowE0BfyuRwODh1t6aAC\nQWqb2ynwuo20zlYrPGDjclm4fl9Ejd9dlpF7UAOb8gspJSWCrmnaB8AHqThWT3PFEaPZUNXIZYeX\n2r6f6Y7dq999zjRueumruHlIpH/VH9BoaveTHbL6Cy1Jw6w5Zuxu4K931jOmJDvuTL7H5Yi4qeva\nglSF3CzWJ26cJ4cNtY2cP204zy+t5Nzy4by0bBte2sihhZHZPkZkB6jat4/jxmSxY/deAi315NLM\noIx2vj2hkJ3rtuFvqWf7nn00N2xgaJaPAq2ZHzobwx1DxEXA1Znhly1k4nNlU+3LQHtkMNft1zjb\n7YLmXHjjdZ5aUk1Oq4sLs/Lhq3omN9aw1d/KqpV55NNIPV40HBx+z3vGMaUrKSfDyb5GKM72GJPI\n3xg/IOKhky6XVtOkqPmZiyckTSGxkaMA68hkZ6gzGJCTERHyanUXPfzBRn5z/vTIKBe/tNAj22H9\n7a/8+xI+XFdlCOAlj+k5448tG2jsE8sqtJuAXrWz3kgPHdA02zBTq4sxln87livCfA0XPPwpbYEg\nr1x7hG37WhKcFLU7thm7ZGvfeeJzQA9d/dO3ZhrbfQGNGbe/HbGvtXOy+tDtiq+Yo9/kMZyOrk3Y\ndsQBb6Hne908+O3ymO9bLWcz3zp0JH98Z32Un9OMjID42fMrIl7LmG/5UFkfOvMNXJyt+7v3NrQx\nc2SBseQd4OQpg3nz63ABjUyLoMfz78uhtxwlaGhoOGgii/z8Qj6va+XzJoChjCkZxfyqPewK6Nc6\n2JXJt884jhde+pJ3Vu+lqlofXo7Jz+a+cw/i/Ic/NjqGHNHC/xw+kBc/XcPhw9xs3rmbXFrIoYXL\nK4rZunM3lbv3sHV7CwMz2hktGsnxbUZbvoRL/Q043Br4gRcf5iLgIoDXYUUmBDRBHdnUaLnUkEuN\nlgO7izjSlUGGVswap5tMdwmrHS5qtBxOHFHMovU+2tGvWbpcmttNFnqc0oV2wtDSbr/4RaYd8Lgc\nEb5nqwHw/NJKfnP+9AiRkKt5pUtOuh6WbavlyY83c+GhI8l0O/lwnf2S/ea22J2+XZphicwgCHpn\nkoiHwEi9QLizEkLEFHTz/bl4S5yV0SZxTnR9Raw5hnjZM1/7cqdF0O0ieWJb6F9V1kWMuCQ7LaO7\nnshUecALekfYWTdF2R5DDM2TJd+sGMHJUwdz+ZOfG9uskynSYpdDvFg3oPmmL8oOT2DmZ7nZ5wn7\nQAfkeiI+p1dgCt+8MjTOitfjNK5BxmKbF+S8f/0xTLtlvmEZ5mRG3ipSrHIyXBGrMfOz3PgC4Y6h\niSz2aHD26adx67L5tDpyWBLQEy2dOWMohWfP5Mm31/HH7bpfvdDtpqZdt3bG5GWzqa4RL23MHOTk\n6Ysn89QHK3j3i7UU0kihaKQs30db/T4KRSMFNDBUVFPcvo0jnA14m9tAn+8F+TV9ApdkQpOWQQ25\nDHpuCE+5NYZ/PoymyiCZeQNocxcw1CGoJYcaLYcacqnVcqjHG+Unh7DwWIVAfr8eV2QESyxxMYd7\nyoVoEwfpEcHmuYBb/7uKFZV1RjlFOxpMK2+tk7X3v72OHxw91lZgzJuCmpaQkEqx09CTkD34/kY2\n3HVKzCyPiU7smtMIJyrosSz5eKkl9GCA6AgjM9Z5Atn5Oh2CTzfZh+Va5ypSkU6iI5Sgd4DdMuB3\nf3q0YWWZXTK3njElyqK3CroU8tJQ0qFB+ZnYYb4xC7PDop2X6Y44h1ylKkl04cnQgiw27G0EYEi+\nHgFQ2xK+ATNcTsaUZBsx+NKSlciHMjfTHeGfzs9yRwnFTSF/q8shbOcjrCkIJJuqmgBBM5m0ZBRA\nyQR25Qb5IBieWD59xFD+uyKyZqvH6aA9EOTkCfl8sW4zc8dlsGrTFgpp4C9nl7JuyzYWrFhLoWjk\njNxMsndvRtu5gtOdjeQ3NeEQGpd4iMKvOah9JId3PDnUkEOtpo8I8nYPhIWTObZ5PyUOp76dHERj\nK25caFrkwyy/96jj2wqsvs0aMbFwfRW/MiUAC4QSi0m+99TSiPfM+AIaa/c0GL+7mYgiKkHN1p0Y\nFWFkulf/tmAzoN8fiVjo8TBfc6zPWNuybGsNFz+2iI9uPDYisiWehd7uD/K+KTGeOez1tIOG8NqX\nu6JEXrrHsj3OqCyqEuuCqnircVOFEvQOkPfLD48dx2kH6dEdhdkeQ2TNgu10iAiLftOvT41KhSrF\n+MwZQ8nLcnHMhLCvc8WvTjSGbjtNkTPFJkHPz3IbaW4hPKknSbQ4wowRBWFBDy2gstb1NE8UFXlt\nFI5wlI4kO8Nl+BuvP2kiV31jjHEcp0PYTiCb5w9itV+KlVlgKkYVcurUwVGCLkcVGd5s9lBEXf5w\nPg2GJsQOOY2t2Xu4c6meb+eMb5/COb94w/isgyB5NFEoGimkgQLRSCGNFIiG0LbQ3zQyXFQx3bmZ\nvLoGePd5/gfCIwHQRwaZ0LQ9m9bdBZzmyaKGPKqDuewnj+qQm6hayyO4fRA5zfXkOVppdWQZS//l\ntVijdTRNL/BhXLM/GFEwpbIm/tqJTJfT1kJvMrlqglpilZvMYit/+3Z/MLZvPahFdUBgv/7BOEeM\nzuFqU8cF4Xws763ewyWzS43t1/17OQ1tfi4xFZE38+vXw1kgzeIvR6++oEa7P8hLyyoZnJ9pdL75\nXnccQY+9urS7UILeAfInGDcwxwilMmN+iFyWG9QRSsUqObd8OOcdrC9UEUJwbNmgiP3zvZHiLCky\nC7rXHTEqsLpCEsHjdDC2JMd4PTBXt/Ktgm6eECvJjRwJGOfPiDx/fYvPyFf+jfEDIjqFVl8g4uaX\ny+EzTNcTazgut5utnoOGFxg1YiUjirLYHsoVIttmbaM5J77180Ec1JJLrZbL5lgRRCYuqBjOc0sq\nyaSNUVmtDHG3kBOsR2veT6FooJAGJuW3M8zTjM+3l/LcJlrrtlBEPRnC9MA/9huuAq7yQJvmZr8z\nl/1aLg3NBexxZ+OpGshkp9voAPxaEXv92WhkU0sOW/c3Ge6ZRAhoWlRMuaZpXPuvZeF9gpqRHmLS\nkDxjQdeuulYeWbCRiw4bxV8/3BjR2UjNavMHOfWBhTHPv3RrDdc8HT7XK8t3RBQQh/Bo0+tx0h5j\nsvNtU/pnM3ahszf/ZyXnlQ+3rd1rvq/Mfxdl6/e9zx/k5D8sMIqp/OyECQAUZHnYjn3naU0JrXzo\nfQA55I0V7SLznziEfcZGs2DMGlPUqWXJVgvdbPmYxcrtFBE3ssflYFSRl/WWIb43w0lORrgd0m1j\ntcbMlzMgJ8M29M98/gyXIyKCxGqBmcX8yHEDjNC3jA5qqUK4oID5GPlZ7qiCJSMKvbaC/rMTJhiL\ngqwRRUeMKzbSE/z2/On8PDSBDfo6hVHF2TELVP9y7mRWbK9j7R5Y25LBxAmTuOnUMmbfHY66ITTv\nNygvg3tPOygUXaGRQwuFooFiGvj7hWNYsHwNazduJidYRxENFIl6SoKNTBe7KWlezilu04KkIHom\npkwIaoLah7LRBgzh3x4nNZreGVSTR42WS7Wmjwj2a7ns1/LYTy6BoBZlMVojetr8AepafEwbls8/\nrzyM6bfpo8elW2tYurWGdn+QB97bYOxvDiVt9wcjLNTrT5rIb94K54S54K+fRpzrun8vj7q/pMtF\nzsuAHtu/oaqRORMHEg+7xWcAf3pvPTecXBY1GjD72M3tHh5y2/gCwYjKWDJfkUwdYUf0Cmol6L2O\n/A2yYgl6qLd3OexFySxqHRWBjoXZQs/LckdYRFK0HAKmDsvni221xmc+vvFYQBfDF5ZuNx6+bI/L\nmJwF/aY8cfIgLju8lGv/tYypQ/NDxwy33Tr5KolwgZQWRuRtsVq/Zv5ycbnhLkpE0GeN0UsJmkWj\nwOuOOsfIIq8RdievMTvDxQ+OGWvsYx1JXTKr1Gj32TOHRQj68l+diKZp3P/O+qgFURWjCsnLdPPz\nkyZy1T+WGNdsLS0YPq+D2aGSiCBoxEuj5mU7g6gfMYdl28bwwuZKGkzzIC6Hnpu9bHAum3bvp4BG\nikU9IzOb8bTVUCR04S+igTkeBxrbGS12cbBjLUU04LQLHQUCf/OieYv5j8etdwDkId5awDxXJUEc\nBHAwYWc+NS0BHE4n3sWLuNq5ngAOgqF/A9Ys4tvOxtA2AQ1O2h2CIA7ca+s5ybHGONbU5lq+4dhm\nfD6gOUzHEpQNLeCrnfqxAjhg3wa0/TUMF1WUZuRS0xKAxr2ccd+7BHDw5a2ngMOFGz9+HGiWPIO+\nQJA/2ixgM1aeWiZsm0zuHbPLRa7wtYpxIBjE6RBRLsd4pKJyVEcoQe+IUC9vndyUyO2xamWasVqG\niVJkmvjMz3JH+BalaAW18MTlnIklPHH5ocY+I4u9XHvsOEPQhxdmRVjWQggeubQC0AVMYhZ0c6cy\ntiSb0QN0l81JUwYzf9Ue7jxrKre/GpmVMt53IsM3oWNBP3HyIO46W68Udc85B3HUb94HdDeAlZHF\n4QnTomwPbqdgcH6ku8haKDzP9FA6HYKnrzyMix5dZGwTQjBnYkmEoGe4HNx6xhQgepQUyw3mcgoy\nXE6+degInlkcuYKyqV3P9mjtbKSQtPgCtONmL4Xs1QrZJdzUBiNHVJcNHsXfN28Nt5sgq/93Nqfe\n/TJF1IfEv4Ei6rl8ai6Z7TXU1WygWNQzXuzAtXwZ33G24yCIkyDOfSYBeh/+12qM7oULYxmo8+Gv\nZhvgczja3ibQqQbMP9Of9dwiH2UAMlDrt/CFjCG4R/9vfeh1UBP4cRDAiR8nrk88NPo0LszQXwc0\nB36c5K7KhB3ZuISTlz1N4c9o4c+O+DKX0e5W/DiZ/vmL/MZVRdmaQvJdjfqxcDJlcxF5rkbG1BUw\n3mk6jqkNAfRz+jX9M4GWmYB9gZ1UoQS9A+QtHcvaHB5ayp9I/pDOWuhWl4tZ0M1iIoXJrq1mV89f\nLipPqCKO1PM7zpoa8fmHLj6YCSF/bb7XzaOX6Z2B12PxVccRdLPYd/Td/fj48caxRxZ7OWP6UP5v\nxU4OGp7PKlOirt9fMJ3xA3O5D31oX+j18PZPjma4JSeP2+KmybMMm48YNyCqDVZX2feOGsPUYfpI\nxmyluZwialQmh94yL7zdYp0vttXiC2g4HQ4mD8mLuC6IjnKxi5iwJjbTcJCZN4BN2lA2MZSxA7Ip\nmzaEB97bwHGHHE5epovLvgwXDX/5u4dz9l8+MV4fNX4A63bVctzEAdx11mQm3/w6zpBN7Qz9E2j6\n3yK83UGQ3507hZteXGFs+/35U7nx+eW2+zoJcsTYIhZtrDKOa5xDBDl8dCGfb94Xcd5fnjoBhxbk\nvje+xkUQpwjgQv4LMmNwDqt37MdJEBcBnEL/f2xmJgMLvAR97dRre3ASwCWCZAgfLtpwEiC7uYFx\nogUnQbx7dnG4s5Hc/TDY2a6fiwCZe4IcLPy4qwKcENvrEsGutm8CQxPbuZMoQe+AoLFQwv79w40h\ndJi/X3GorUUfzwURjyKLoJvLz5nPIy30WB3HwhvmUJzjwetxkZ3RcZEBKUyjiyOXTMeyqKMnHxO7\n3opRRTHfe/WHRzIl5AKS3HfeQdx4ShmZbmfEqOec8uERMfEOQdRyb4j+HazttsP6nX7v6LALJ9JC\nj9wvN8NFtT9SaM0d3VkzhvKf5Tu56aWvOLd8OG6n4B/fPZQ7Xl3FK8vD0TvWTs8ussJas9bKy9ce\nwVeVdfDeBgJBLWri0Jq6oKk9QG1bkJxsL7izaMY+xBagJCeDlvaA4a7Y4RnDai087+AdM5ulWuie\ns/E86IOf0bbHLh07kX9tWBuxbaxrKhcdNopFX33CUpsC0ifnDOZNvx4FdMvpk7ntv6sA+M6YUqac\nMYWahlYuu+td+4sJfY13nDmFsbNLOf7mN/nGiAER9XdnjSli1c56Lju8lD+/t84QehcBoxOR/5wi\nyIKffYMhBfYRNqnkgKspmizSzxvLR26XkfHoCSUcUhotUp210M0TL3mZbs435WU2z9jLZeyxzjOi\nyEdTdGsAABJLSURBVGtYutkJiFg4617kExjLlWLtxKwW+tkzh+F2iqhl3iOKvCy8YY7tMe0mnTLd\nTiPG2Dopao7Lj7Xk3WqhD87PZPKQPB4NuZ1An8S7/qSJxmtzJ/b4dyoivj+zi8Uq6HYRFTLR1i9P\nm8Svz5kG6KOr6qY2inM8DMjJ0NMBm4iXXkISK//2iZMHcW75cPIy3cZv5w9EFyqx5i2vaW6n1ReM\nGsGYkd/lsIIsZowoMLY/9OGGiP0Sud9iYRdh9YuXV/Lz51fYijkQuXraZv4rXtk4iRyBuZ2CBZYi\n2p9t2o/L6SAnw4WGAx8uWsmgES915FBNPnsoYgcl1HiGwYDx4Irnc0oNStA74A/fnMG8U8qYNMQ+\nJMwaBx6PRCx08+SdxPwweFwOrjHtY56sNdqSwNyLtCpjTfZCeFQSzvgnQtvthdL60FqF//5vzmD9\nXacy3fTgS2JZ/fHEBOy/0yGhxVqxPD7WzjnT7eT1H3+D4yeHw0ivnTOOa+eMs22fNdzU6kM3YxcZ\nJJscCGp4PS4unjUSp0Owt76Ngbl62zPi/C7HTxpku91aiOWJ7xwCwCOXVvC7C6ZHtK+ytoXVuyPd\nOtYOQcayx4vkOHHyYEDvzM0L4KwrlN1OwcfzjuU3cVa3xiLW/FWsUo1W7Fx/MlLsf+aMi5laV85f\nuJwOY7GReUSoaVrMjkquN5k5soCvbjspoXamAiXoHTAwL5PvHz02poglEqGRzL43nlzGd4+MHHrK\nuqTyQTa3xc76iFdOTCIf0suPKI25jyymPW6gPgH65OWHcOWRoxkaY3VrdowHLxFijSpyO7Ds3DYj\nJykAseqXWq36RIg3uspwOQzRsHYWdhkH5byLLCBekpNJTbOPnXUtxpqAWDmEXvj+bFs3n0QuhLl0\n9ijmlEWH9snRwQ0vfMkNpjJ/EI6/lm2Qvv/RNm4riVw74fU44/7+bqeDYQVZcb/H/Cw3f7koOq/S\nSVMGx/3ciZMH8elNx8Z83+5xkO6mCYNz+dah4RGveYJc1vo1j2Q+mhceSbb7gzGjXEYV6d+Z3f3Z\nnShB7yLJ5DdO1If+i1MnRbx2OR2sv+sU/nrJwVH7Siv4jOlDGR8S3otjrIYzk+Vxsur2k/j5iRNj\n7nNO+XA23HWKIUDjB+Xyy7mTY16z1yS+hV53UqMX8wN76rTBfHrTsWy557S4ibLAXpxlJxfrWerM\nXEY8QRFCGOeUFvBnNx3HOz892naByyWzRvHQReWcW653mDIktLbZZ4ipbOMZ0yMn0SpKiyJE5O6Q\ny8Z8bIjTmcX5PqVb57fnT+d7R4Vz8csO3Q65ajnL44qbo0We97RpQ5h3ShnLf3VCVC4aTdNsRx9u\np4M7z5oa89hul4Mh+Vk8e/Us2/ftcspIC93tEBEjEDkPdfuZU4y/5QT0waMKGZAddv+0m2rtWhk7\nUBf0zhgPXUFNivYgifrQzSL21a16GGE8Efr6tpPIcDlwOgRf3npiwkJqjUqxI9GJTQiPQE47aEjc\nDJZ2mIeyvoBmm2fEDvm9mHXKEPSYPvROCHoHn8l0O2hsCx97cGgUI4Xj3Z8dbaRPcDgEp0wLFwkx\nC8rwUOHysSU5PHppBUeOH8CMEQXc/uoqYx+zoI8qjgyDO37yIB54bwNzDxqCHfEERk60elwO497w\nOB1GJ2OHzMPjdTvj1gWQRoDL6eD7oQnlkUWRbdeI/YzEc/tIF6Ncq2DF7C9fsL6Knz23grNm6h2l\n2+mIcOvJOQ/z+eRtdEHFcBwOwcGjClm6tQZfQLMV9EcvreDzrfpqskTCmVOJstB7kM5MilqTYkm+\naZoYzc5w4XI6EEIkZRWnGjm0jSWk8XBFCHriVXmk5WeOa483L6Cfq2vts0OGNcbab1BeZoSP2Yw5\nN89kU2z98ZMHkel2coXFBZeTEf6NrS6eKUPz2XLPaUZxayuxJvchXKjF7XQYhVky3A5DjD+6cU5U\nLhTZdo/LETU3UhzjeiWzxhRzxRGma4vjKcyOY3x0tKTefD9tqmrixWWVRt54lzPSQr/rrKmMG5hj\nTIhC2MCShsLVptGL9KGPLPLy1nVHGa4uWZnMrjpUd6IEvQfweiKH46ng3vMOMooa9BWMEM8uHicZ\nQZcPmblAidwWywXQFb+m+WE2I0cnHstv/OzVs/je0WPi+pfN2SbjuTckZgu9YlRhhI+7I4swnstl\n8WbdqsxwOYxzmPcfXug1LFvQKzHJ1MEFXjdXf2MML/7gcOP9RNyRR46PPR9gxi5aSGI3ZzQgJ9yZ\nxFty73E6IiKjDhtTzDs/PToi15E0UKShIJ/nDJfDsND9gSATB+dy+5lTcTqEUdqyp+uKKkFPARce\nMoKTpthHHkD4RhBdlrq+zZjQQ3Do6Nhx5fF47nuzAfD5E18ineVxsvbOk/nZiROMbXJC0VqmTtJZ\nv+aWe07jfy3zGxIZlWK10KcMzeemUybFFTfzxHa8ClmSESZXhcMh+OkJE+LsHYndtd97bqQfXp/A\n1M9hzX7pcYbb99z3Z7MxlCdo3MAcHA5B+cgCfnrCBN792dFGNM8zV9n7tiFytBHvV4836rKz0M2r\niEtyYruMnA4RM/GcRPZpslMxC7rcZs0GKUfWse7B7kL50FPAPefGD8V66ruH8c9FWzscgqY7M0YU\nsOD6OYwoSrzCuhk5gkm0fqTEuopTPvwxBd0huPLI0UY65FQgxSueBRyLzCQTthVlezhoeH44xDEJ\nV56dy8U6KvC4HIZ/2zrKMY8mBuZmcvGsUbyzeg/fGF8C6Fa5LIgtLduRxbGXu5t90HbJtCSxQhfB\n3gKXcx7Th+dzTvkwo2KYlZpmHwd38FwKi4UuO12Py2mUkrzhpMjgAqPSlBL0/sfkoXn8+uxpHe/Y\nD4j38HaEXBE7piR2mFwiSKsp1nBXCMEv507u0jmsyGyQM2xi7DsiVoEPM789fzrb94ezLb5y7RGG\n0CQzN2O20CcOymXtngaj2IrE7RQMztM7C2vnbJ0cPmpCCZvutnf9Wcss2mF2H8W10OO5XGyiWKTf\n+5RpQ+KOjiYNye1wfkSmeZZCLoU9w+XA43LYuj5luK25WEZPoAS9j/LUdw81ev8DhVHF2Tx95WHM\nHJm8KJoZH8ozU5Ibe6l6qvnDN2fy2aZq49zJkIiFLvPoS8wilUjqAol5BPHWT46y3cfjcuByOnjq\nu4dGWe+JdD6Sx79zCP9esj0ql46ZSAs99rHMgv6vqw7j238LJ0+zc7lIEY4XnfTFzScYE9UOEbsz\ntrpc5DxFvOsyXC5xCrp3B0rQ+yhyCHugYZcYK1kuPmwko4uzOWJcYhNuqWD22GJTatzkSEYk7Ygb\n0mchkTBUKYJ292AyGUPHlORw0yn2cw4S80pLzWKjDyvIMkYIZh/64WMj7xE7QZffaTyr35yyYd2d\np8SMzhKWzmFUcTZ3nT2Vk6YMjnlsOfJIYI1fSlGCruh3CCE4cnzXO4aeIlkfupVYla7s6MjH/9qP\njqQgzsgwXkqCzhDPXfTuz442LdiKvZ+dD/24SYN4ZflOmuLUEjUfM15HN21YPjtqWyI63osOi794\nL57PvztRgq5Q9DI9aqF3IOjWzJZWOpvTPxGs1myi5zJb6HMmltAeCHJ6aML7WJv0B9NHFMSsQGXH\n7y6Yznd2lBqT0IkghGD2mOIoV1l3owRdoehlkolSsf984tZgrDj1MSXZCcVMp3IthcQh9AItUpYv\nmz2Kv3+6tcO0DxKzoJsLu5jTJlx02EieXrSN+785nVOmDqEticnK7AxXzFWo8XgmRiqC7kQJukLR\nyySTD6i7zvXuT49OyN/bHW1992fHMOe3HxiKfusZU7jtzOjcLRfPGmk7gkjEvXHX2dO4yxRplki8\nfzqiBF2h6CckEx1kzZgphIhZxKW7kW4gOSkaq9O486ywIH8y71gyXA6eX1rJWTOGdX8j0wQl6ApF\nH+Cec6bZ1khNlI2/PjXhdcgvfH82o4q7FuufSuTkZDIRIUNDOcy/b6ocpVCCrlD0CS60VChKlmSy\n+sVK3NVbyMVOPRzh1y9Rgq5QKJLipWsOj8gL3lU6ky5BYY8SdIVCkRTlIwtTejyX4XJRNnpXUdkW\nFQpFrxKeFFV0FSXoCoWiV5GTot5+GkrYk3Ta5SKEGAH8AxiE3rk+omnaH1PVMIVCcWDgdAh+ceok\njp54YOYvSiVd8aH7gZ9pmrZMCJELLBVCvK1p2qqOPqhQKBRmropRCUqRHJ12uWiatkvTtGWhvxuA\n1YCK8FcoFIpeIiU+dCFEKTATWBR/T4VCoVB0F10WdCFEDvAicJ2mafU2718thFgihFhSVVXV1dMp\nFAqFIgZdEnQhhBtdzJ/WNO0lu300TXtE07QKTdMqSkrUpIdCoVB0F50WdKFn0HkMWK1p2u9T1ySF\nQqFQdIauWOhHAJcAxwohlof+nZqidikUCoUiSTodtqhp2keQcII3hUKhUHQzaqWoQqFQ9BNETybE\nEUJUAVs7+fEBwL4UNqe7Saf2qrZ2H+nU3nRqK6RXe7va1lGapnUYVdKjgt4VhBBLNE2r6O12JEo6\ntVe1tftIp/amU1shvdrbU21VLheFQqHoJyhBVygUin5COgn6I73dgCRJp/aqtnYf6dTedGorpFd7\ne6StaeNDVygUCkV80slCVygUCkUc0kLQhRAnCyHWCiE2CCHm9VIbHhdC7BVCrDRtKxJCvC2EWB/6\nvzC0XQghHgi190shRLnpM5eF9l8vhLism9o6QgjxvhBilRDiayHEj/t4ezOFEIuFECtC7b0ttH20\nEGJRqF3/FkJ4QtszQq83hN4vNR3rptD2tUKIk7qjvaHzOIUQXwghXk2Dtm4RQnwVWs29JLStr94L\nBUKIF4QQa4QQq4UQs/twWyeK8Cr55UKIeiHEdb3aXk3T+vQ/wAlsBMYAHmAFMLkX2nEUUA6sNG27\nD5gX+nsecG/o71OBN9BX0s4CFoW2FwGbQv8Xhv4u7Ia2DgHKQ3/nAuuAyX24vQLICf3tRk/DPAt4\nDrgwtP1h4Aehv68BHg79fSHw79Dfk0P3RwYwOnTfOLvpfvgp8C/g1dDrvtzWLcAAy7a+ei/8Hbgy\n9LcHKOirbbW02wnsBkb1Znu77QJT+EXNBt4yvb4JuKmX2lJKpKD/f3tn81JFFAXw3wHpS0JLIiSD\nFKRWYRF9kEQUCUq4cqEERQVt2gch9DfUIkIIWhoUVNLGSlu1qTQrKyQjQcWPCHTRqo/T4p55ji+N\nAn1z3+P8YHh37gwzv+Hdd97Mufe9OwJUW7kaGLFyF9CRvx/QAXSl6hftt4reD4ATxeALbAAGgQOE\nH2KU5bcDoBc4ZOUy20/y20Z6vxV2rAH6gGPAQzt3lK527DH+DOjRtQWgAviM9e3F7LqEexPwLGvf\nYki5bAPGU+sTxDMz0lZVnbLyNGF+VVjeueDXIosnH4nW11IYQ8As8Jhwxzqnqj+WOHfOy7bPA1UF\n9L0KXAJ+2XpVxK4Q5vx9JCIDInLB6mJsC7XAF+CWpbNuikh5pK75tAPdVs7MtxgCelGg4as1qiFD\n8pfJR2LzVdWfqtpAuPvdD+zKWGlJROQkMKuqA1m7/AeNqroXaAYuisiR9MaI2kIZIa15Q1X3AN8I\nKYscEbnmsP6SVuBO/rZC+xZDQJ8EtqfWa6wuBmZEpBrAXmetfjnngl2LLD35SLS+Cao6BzwlpC0q\nRST5R9D0uXNetr0C+Fog38NAq4iMAbcJaZdrkboCoKqT9joL3CN8YcbYFiaACVVNprK8SwjwMbqm\naQYGVXXG1jPzLYaA/gKot1EEawiPNj0ZOyX0AEmP9BlCrjqpP2292geBeXsE6wWaRGST9Xw3Wd2K\nIrLs5COx+m4RkUorryfk+z8QAnvbMr7JdbQB/XYn1AO028iSWqAeeL6Srqp6WVVrVHUHoS32q+qp\nGF0BRKRcRDYmZcJ7OEyEbUFVp4FxEdlpVceB9zG65tHBQrol8crGdzU7Claww6GFMFLjE9CZkUM3\nMAV8J9xJnCfkQvuAj8ATYLPtK8B1830L7Esd5xwwasvZVXJtJDzmvQGGbGmJ2Hc38Mp8h4ErVl9H\nCHKjhMfZtVa/ztZHbXtd6liddh0jQPMqt4mjLIxyidLVvF7b8i75/ETcFhqAl9YW7hNGfUTpaucp\nJzxxVaTqMvP1X4o6juOUCMWQcnEcx3H+AQ/ojuM4JYIHdMdxnBLBA7rjOE6J4AHdcRynRPCA7jiO\nUyJ4QHccxykRPKA7juOUCL8B+O/KNKrXoUUAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x10c0eb850>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"<matplotlib.figure.Figure at 0x10c0eb850>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"\n",
"import matplotlib.pyplot as plt\n",
"from IPython import display\n",
"import cPickle\n",
"\n",
"reader_dict = {\n",
" 'user_id': 0,\n",
" 'gender_id': 1,\n",
" 'age_id': 2,\n",
" 'job_id': 3,\n",
" 'movie_id': 4,\n",
" 'category_id': 5,\n",
" 'movie_title': 6,\n",
" 'score': 7\n",
"}\n",
"\n",
"step=0\n",
"\n",
"train_costs=[],[]\n",
"test_costs=[],[]\n",
"\n",
"def event_handler(event):\n",
" global step\n",
" global train_costs\n",
" global test_costs\n",
" if isinstance(event, paddle.event.EndIteration):\n",
" need_plot = False\n",
" if step % 10 == 0: # every 10 batches, record a train cost\n",
" train_costs[0].append(step)\n",
" train_costs[1].append(event.cost)\n",
" \n",
" if step % 1000 == 0: # every 1000 batches, record a test cost\n",
" result = trainer.test(reader=paddle.reader.batched(\n",
" paddle.dataset.movielens.test(), batch_size=256))\n",
" test_costs[0].append(step)\n",
" test_costs[1].append(result.cost)\n",
" \n",
" if step % 100 == 0: # every 100 batches, update cost plot\n",
" plt.plot(*train_costs)\n",
" plt.plot(*test_costs)\n",
" plt.legend(['Train Cost', 'Test Cost'], loc='upper left')\n",
" display.clear_output(wait=True)\n",
" display.display(plt.gcf())\n",
" plt.gcf().clear()\n",
" step += 1\n",
"\n",
"trainer.train(\n",
" reader=paddle.reader.batched(\n",
" paddle.reader.shuffle(\n",
" paddle.dataset.movielens.train(), buf_size=8192),\n",
" batch_size=256),\n",
" event_handler=event_handler,\n",
" reader_dict=reader_dict,\n",
" num_passes=2)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 应用模型\n",
"\n",
"在训练了几轮以后,您可以对模型进行推断。我们可以使用任意一个用户ID和电影ID,来预测该用户对该电影的评分。示例程序为:"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"[INFO 2017-03-02 17:49:59,181 networks.py:1472] The input order is [user_id, gender_id, age_id, job_id, movie_id, category_id, movie_title]\n",
"[INFO 2017-03-02 17:49:59,195 networks.py:1478] The output order is [__cos_sim_0__]\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"[Predict] User 234 Rating Movie 345 With Score 4.12\n"
]
}
],
"source": [
"import copy\n",
"user_id = 234\n",
"movie_id = 345\n",
"\n",
"user = user_info[user_id]\n",
"movie = movie_info[movie_id]\n",
"\n",
"feature = user.value() + movie.value()\n",
"\n",
"def reader():\n",
" yield feature\n",
"\n",
"infer_dict = copy.copy(reader_dict)\n",
"del infer_dict['score']\n",
"\n",
"prediction = paddle.infer(output=inference, parameters=parameters, reader=paddle.reader.batched(reader, batch_size=32))\n",
"score = (prediction[0][0] + 5.0) / 2\n",
"print \"[Predict] User %d Rating Movie %d With Score %.2f\"%(user_id, movie_id, score)\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 总结\n",
"\n",
"本章介绍了传统的推荐系统方法和YouTube的深度神经网络推荐系统,并以电影推荐为例,使用PaddlePaddle训练了一个个性化推荐神经网络模型。推荐系统几乎涵盖了电商系统、社交网络、广告推荐、搜索引擎等领域的方方面面,而在图像处理、自然语言处理等领域已经发挥重要作用的深度学习技术,也将会在推荐系统领域大放异彩。\n",
"\n",
"## 参考文献\n",
"\n",
"1. [Peter Brusilovsky](https://en.wikipedia.org/wiki/Peter_Brusilovsky) (2007). *The Adaptive Web*. p. 325.\n",
"2. Robin Burke , [Hybrid Web Recommender Systems](http://www.dcs.warwick.ac.uk/~acristea/courses/CS411/2010/Book%20-%20The%20Adaptive%20Web/HybridWebRecommenderSystems.pdf), pp. 377-408, The Adaptive Web, Peter Brusilovsky, Alfred Kobsa, Wolfgang Nejdl (Ed.), Lecture Notes in Computer Science, Springer-Verlag, Berlin, Germany, Lecture Notes in Computer Science, Vol. 4321, May 2007, 978-3-540-72078-2.\n",
"3. P. Resnick, N. Iacovou, etc. “[GroupLens: An Open Architecture for Collaborative Filtering of Netnews](http://ccs.mit.edu/papers/CCSWP165.html)”, Proceedings of ACM Conference on Computer Supported Cooperative Work, CSCW 1994. pp.175-186.\n",
"4. Sarwar, Badrul, et al. \"[Item-based collaborative filtering recommendation algorithms.](http://files.grouplens.org/papers/www10_sarwar.pdf)\" *Proceedings of the 10th international conference on World Wide Web*. ACM, 2001.\n",
"5. Kautz, Henry, Bart Selman, and Mehul Shah. \"[Referral Web: combining social networks and collaborative filtering.](http://www.cs.cornell.edu/selman/papers/pdf/97.cacm.refweb.pdf)\" Communications of the ACM 40.3 (1997): 63-65. APA\n",
"6. Yuan, Jianbo, et al. [\"Solving Cold-Start Problem in Large-scale Recommendation Engines: A Deep Learning Approach.\"](https://arxiv.org/pdf/1611.05480v1.pdf) *arXiv preprint arXiv:1611.05480* (2016).\n",
"7. Covington P, Adams J, Sargin E. [Deep neural networks for youtube recommendations](https://static.googleusercontent.com/media/research.google.com/zh-CN//pubs/archive/45530.pdf)[C]//Proceedings of the 10th ACM Conference on Recommender Systems. ACM, 2016: 191-198.\n",
"\n",
"<br/>\n",
"<a rel=\"license\" href=\"http://creativecommons.org/licenses/by-nc-sa/4.0/\"><img alt=\"知识共享许可协议\" style=\"border-width:0\" src=\"https://i.creativecommons.org/l/by-nc-sa/4.0/88x31.png\" /></a><br /><span xmlns:dct=\"http://purl.org/dc/terms/\" href=\"http://purl.org/dc/dcmitype/Text\" property=\"dct:title\" rel=\"dct:type\">本教程</span> 由 <a xmlns:cc=\"http://creativecommons.org/ns#\" href=\"http://book.paddlepaddle.org\" property=\"cc:attributionName\" rel=\"cc:attributionURL\">PaddlePaddle</a> 创作,采用 <a rel=\"license\" href=\"http://creativecommons.org/licenses/by-nc-sa/4.0/\">知识共享 署名-非商业性使用-相同方式共享 4.0 国际 许可协议</a>进行许可。\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.13"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Markdown is supported
0% .
You are about to add 0 people to the discussion. Proceed with caution.
先完成此消息的编辑!
想要评论请 注册