{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "source": [ "%pylab inline\n", "\n", "\n", "from scipy.spatial.distance import cdist\n", "from matplotlib.ticker import NullFormatter\n", "CBB_PALETTE = [\"#000000\", \"#E69F00\", \"#56B4E9\", \"#009E73\", \"#F0E442\", \"#0072B2\", \"#D55E00\", \"#CC79A7\"]\n", "\n", "\n", "def dtw(X, Y, metric='euclidean'):\n", " \"\"\"\n", " \n", " Parameters\n", " ----------\n", " X : array_like\n", " M x D matrix\n", " Y : array_like\n", " N x D matrix \n", " metric : string\n", " The distance metric to use. \n", " Can be :\n", " 'braycurtis', 'canberra', 'chebyshev', 'cityblock', 'correlation',\n", " 'cosine', 'dice', 'euclidean', 'hamming', 'jaccard', 'kulsinski',\n", " 'mahalanobis', 'matching', 'minkowski', 'rogerstanimoto', 'russellrao',\n", " 'seuclidean', 'sokalmichener', 'sokalsneath', 'sqeuclidean',\n", " 'wminkowski', 'yule'.\n", " See: https://docs.scipy.org/doc/scipy/reference/generated/scipy.spatial.distance.cdist.html\n", " Returns\n", " -------\n", " total_cost : float\n", " Total (minimum) cost of warping\n", " pointwise_cost : array_like\n", " M x N matrix with cost at each (i, j)\n", " accumulated_cost : array_like\n", " M x N matrix with (minimum) cost accumulated till (i,j)\n", " having started from (0, 0)\n", " \n", " \"\"\"\n", " \n", " X = np.array(X)\n", " Y = np.array(Y)\n", " if len(X.shape) == 1:\n", " # Reshape to N x 1 form\n", " X = X[:, np.newaxis]\n", " if len(Y.shape) == 1:\n", " Y = Y[:, np.newaxis]\n", " # m = X length\n", " # n = Y length\n", " m, n = X.shape[0], Y.shape[0]\n", " D = np.zeros((m+1, n+1))\n", " D[1:, 0] = np.inf\n", " D[0, 1:] = np.inf\n", " D[1:,1:] = cdist(X, Y, metric)\n", " pointwise_cost = D[1:,1:].copy()\n", " for i in range(0, m):\n", " for j in range(0, n):\n", " cost = D[i+1, j+1]\n", " D[i+1, j+1] += cost + min (D[i, j+1], D[i+1, j], D[i, j]) \n", " accumulated_cost = D[1:, 1:]\n", " total_cost = D[m, n]/sum(D.shape)\n", " return total_cost, pointwise_cost, accumulated_cost\n", "\n", "\n", "def get_path(D):\n", " \"\"\"Traceback path of minimum cost \n", " \n", " Given accumulated cost matrix D,\n", " trace back the minimum cost path\n", " \n", " Parameters\n", " -----------\n", " \n", " D : array_like\n", " M x N matrix as obtained from `accumulated_cost` using:\n", " total_cost, pointwise_cost, accumulated_cost = dtw(X, Y, metric='euclidean')\n", " \n", " Returns\n", " -------\n", " traceback_x, traceback_x : array_like\n", " M x 1 and N x 1 array containing indices of movement\n", " starting from (0, 0) going to (M-1, N-1)\n", " \"\"\"\n", " m , n = D.shape\n", " m = m - 1\n", " n = n - 1\n", " # Starting point is the end point\n", " traceback_x, traceback_y = [m], [n]\n", " while (m > 0 and n > 0):\n", " min_idx = np.argmin([D[m-1, n-1], D[m, n-1], D[m-1, n]])\n", " if min_idx == 0:\n", " # move diagonally\n", " m = m - 1\n", " n = n - 1\n", " elif min_idx == 1:\n", " # move vertically\n", " n = n - 1\n", " else:\n", " # move horizontally\n", " m = m - 1\n", " traceback_x.insert(0, m)\n", " traceback_y.insert(0, n)\n", " # End point is the starting point\n", " traceback_x.insert(0, 0)\n", " traceback_y.insert(0, 0)\n", " return np.array(traceback_x), np.array(traceback_y)\n", "\n", "def plot_warped_timeseries(x, y, pointwise_cost, \n", " accumulated_cost, path, \n", " colormap=plt.cm.Blues,\n", " linecolor=CBB_PALETTE[-2]):\n", " nullfmt = NullFormatter() \n", " left, width = 0.1, 0.65\n", " bottom, height = 0.1, 0.65\n", " bottom_h = left_h = left + width + 0.02\n", " rect_heatmap = [left, bottom, width, height]\n", " rect_x = [left, bottom_h, width, 0.2]\n", " rect_y = [left_h, bottom, 0.2, height]\n", "\n", " fig = plt.figure(1, figsize=(8, 8))\n", "\n", " axHeatmap = plt.axes(rect_heatmap)\n", " axX = plt.axes(rect_x)\n", " axY = plt.axes(rect_y)\n", "\n", " # no labels\n", " axX.xaxis.set_major_formatter(nullfmt)\n", " axY.yaxis.set_major_formatter(nullfmt)\n", " \n", " axY.plot(y, range(0, len(y)), color=CBB_PALETTE[2])\n", " axX.plot(x, color=CBB_PALETTE[1])\n", " axHeatmap.imshow(accumulated_cost.T, \n", " origin='lower', \n", " cmap=colormap, \n", " interpolation='nearest')\n", " axHeatmap.plot(path[0], path[1], '-x', color=linecolor)\n", " #axHeatmap.xlim((-0.5, accumulated_cost.shape[0]-0.5))\n", " #axHeatmap.ylim((-0.5, accumulated_cost.shape[1]-0.5))\n", " return fig" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Examples" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD8CAYAAABXe05zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xd8VNeZ//HPUUN0gRBNmkFIdDBVYgTY4N5jx44TV8AUA4njOLvJbpJNTzab3d9mE+M4icF0YztxbMc9TtwxAlSopvcimgSiCKGu8/vjzsgCVO5Id+beO3rer5dfIGnmzpGRnrn3POd8r9JaI4QQwj2i7B6AEEKI4EjhFkIIl5HCLYQQLiOFWwghXEYKtxBCuIwUbiGEcBkp3EII4TJSuIUQwmWkcAshhMvEhOKgPXr00KmpqaE4tBBCRKT169ef0lonmXlsSAp3amoq+fn5oTi0EEJEJKXUIbOPlakSIYRwGSncQgjhMlK4hRDCZUIyxy2EEHapqqqioKCA8vJyu4fSoPj4eFJSUoiNjW3xMaRwCyEiSkFBAZ07dyY1NRWllN3DuYTWmtOnT1NQUED//v1bfBxTUyVKqQSl1CtKqZ1KqR1KqQktfkUhhAih8vJyEhMTHVe0AZRSJCYmtvpqwOwZ93zgPa31fUqpOKBDq15VCCFCyIlFO8CKsTV7xq2U6gJMBhYDaK0rtdZnW/3Kwl4nPod9H9k9CiFEC5iZKkkDioClSqmNSqlFSqmOlz9IKTVHKZWvlMovKiqyfKDCQrW18MpM47/aWrtHI4QIkpnCHQOMBf6ktR4DlALfv/xBWuuFWusMrXVGUpKpXZvCLrvfg1O7oewMnN5j92iEEEEyU7gLgAKtdY7/41cwCrlwq+ynID7B+PvhdfaORYgIk5eXx8iRIykvL6e0tJThw4ezdetWS1+j2eak1vqEUuqIUmqw1noXcAOw3dJRiPA5vA6O5MBt/w8+/R/j7+Om2z0qIULi529tY/ux85Yec1jfLvz0S8Mb/XpmZiZ33XUXP/rRjygrK+ORRx5hxIgRlo7B7KqSJ4AX/CtK9gMzLB2FCJ/s+dC+O4x5BPZ/CofX2j0iISLOT37yEzIzM4mPj+fpp5+2/PimCrfWehOQYfmri/Aq2gW73oVrfwBxHcHrg13vwIVC6NTT7tEJYbmmzoxDqbi4mAsXLlBVVUV5eTkdO16xnqNVJKukLcl+GmLaQ+ZjxseeLOPPIzmNP0cIEbQ5c+bwy1/+kocffpjvfe97lh9ftry3FeePwZa/QMYM6JhofK7vaIhuZ8x7D/2SveMTIkKsWLGCmJgYHnroIWpqapg4cSIfffQR119/vWWvIYW7rVj3J9C1MOHxLz4X0w76jpEzbiEsNG3aNKZNmwZAdHQ0OTnW/37JVElbUH4O8pfC8HugW+qlX/P64NgmqCqzZWhCiOBJ4W4L8pdAZQlM+taVX/NkQW0VHNsY/nEJIVpECnekq64wpknSroM+o678usdn/CnLAoVwDSnckW7LX+DCSZj0ZMNf75gIPQbBYZnnFsItpHBHstpaYwlgn1GQdm3jj/P4jAalBE4J4QpSuCPZrneNEKlJT0JTGcDeLCg/awRPCSEcTwp3pNLaCJNK6AdD7276sXUbcSRwSgg3kMIdqQ6vg4I8mPgERDezXD8xHTr0kHluIVxCCnekyp4PHRJh9MPNP1Yp/zy3nHEL0Vo//vGPmT9/ft3HP/zhDy0PmpKdk5GocAfs/jtc+x8QZ/L2oBI4JSLR379v3KbPSr2vgtv+u9Evz5o1i3vvvZcnn3yS2tpa/vznP5Obm2vpEKRwR6I1v4fYDjD+MfPPCcxzH14Hw+4KzbiEaANSU1NJTExk48aNnDx5kjFjxpCYmGjpa0jhjjTnjsKWlyFjJnTobv55gcCpIzlSuEXkaOLMOJRmz57NsmXLOHHiBDNnzrT8+DLHHWlyGgiTMiOmHSSPlVuZCWGBe+65h/fee4+8vDxuueUWy48vZ9yRpOws5C+DEfdCt37BP9/jg7V/MAKnYttbPjwh2oq4uDiuu+46EhISiI6Otvz4csYdSQJhUhMbCJMyw+sPnDq6wdpxCdHG1NbWsm7dOmbNmhWS40vhjhRV5ZDzLKRfD31GtuwYgcApWRYoRItt376dAQMGcMMNNzBw4MCQvIZMlUSKQJjUvQtbfowO3SVwSohWGjZsGPv37w/pa8gZdySorYE1T0Of0dB/SuuOJYFTIgJore0eQqOsGJsU7kiw6104vbf5MCkzvBP8gVO7rBmbEGEWHx/P6dOnHVm8tdacPn2a+Pj4Vh1HpkrcTmtY/ZRxS7KhFqy/9tbbiNNzaOuPJ0SYpaSkUFBQQFFRkd1DaVB8fDwpKSmtOoYUbrc7tAaO5sPtv2k+TMqM7mlG4NSRHOOO8EK4TGxsLP3797d7GCFl6jddKXUQKAFqgGqtdUYoByWCkD3fKLRjHrHmeEoZZ92yEUcIxwpmjvs6rfVoKdoOcnI77PkH+OZau2HG44MzB4zAKSGE40hz0s0CYVKZs609bv15biGE45gt3Br4p1JqvVJqTkMPUErNUUrlK6XyndoUiCjnCuDzl2HstODCpMzoM+qLwCkhhOOYLdyTtNZjgduAx5VSky9/gNZ6odY6Q2udkZSUZOkgRQPW/clYURJsmJQZdYFTa60/thCi1UwVbq31Mf+fhcDfgPGhHJRoRtkZWL8MRnwFEryheQ1vFhzfDJUXQ3N8IUSLNVu4lVIdlVKdA38Hbga2hnpgogn5S6DyAkxqYZiUGZ4sqK2GYxI4JYTTmDnj7gWsVkptBnKBd7TW74V2WKJRVeWw7llIv8G4hVKoePwXVdKgFMJxml3HrbXeD4wKw1iEGZtfgtJCuPrboX2dDt2hx2BpUArhQLIc0E1qa4wlgH3HQOo1oX89rwROCeFEUrjdZOc7ULzPmjApMzxZUH5OAqeEcBgp3G6hNWQ/Bd36WxMmZYZsxBHCkaRwu8WhbDi6HiY+AVHW38OuQYHAKSncQjiKFG63CIRJjX4ofK8ZCJySW5kJ4ShSuN3g5DbY80/wzQv/3de9WXDmIJScDO/rCiEaJYXbDbKfhtiOkBmaO0Y3yeOf55azbiEcQwq30509AltfgXHTrQ+TMqPPKIiJlxsIC+EgUridLhAmlfUNe14/Jg76jpUzbiEcRAq3kwXCpK66DxI89o3D65PAKSEcRAq3k+UtgqpSmBjCMCkzJHBKCEeRwu1UVWWQswAG3Ai9R9g7lrrAKcnnFsIJpHA71eaXoLQIJoU4TMqMDt0haYg0KIVwCCncTlQXJjUWUq+2ezQGjw8KciVwSggHkMLtRDveguL94QuTMsPrD5wq2mn3SIRo86RwO43Wxvb27mkw9Et2j+YLHp/xpywLFMJ2Urid5uBqY/VGOMOkzOieBh2TZJ5bCAeQwu002fONAjnqQbtHcimljLNuOeMWwnZSuJ3kxFbY+z745oY/TMqMusCpE3aPRIg2TQq3k6zxh0ll2BAmZYZ3gvGn5HMLYSsp3E5x9jB8/gqMe9SeMCkzeo80AqfkBsJC2EoKt1Os+5Mxj5z1dbtH0riYOEgeJ2fcQthMCrcTXCyG9cthhM1hUmZ4fHBiiwROCWEj04VbKRWtlNqolHo7lANqk/IWG2FSk2wOkzLD6w+cOrre7pEI0WYFc8b9JLAjVANps6rKIOdZGHgz9Bpu92ial5Jp/CnLAoWwjanCrZRKAe4AFoV2OG3Qphfh4ilje7sbSOCUELaLMfm4p4B/BzqHcCxtTyBMKnkc9Jtk92jM8/hg2+tG4FSUi9ok29+AmPYw6Ga7RxL5CnfC2meMCIdwi4qCCU9A0qDwv3aYNFu4lVJ3AoVa6/VKqWubeNwcYA6A1+u1bIARbcebcOYA3PRz54RJmeGdABuWQ9EOd0zvgNEA/ts8SPBK4Q6Hf/4QDqyCjj3D/9olxyG6Hdzxm/C/dpiYOeOeBNyllLodiAe6KKVWaq0fqf8grfVCYCFARkaGDW+zLqM1rH4KuqfDkDvtHk1wvP7AqcPr3FO48xZD1UUj3fBisXPXykeCE1th7wdw/Y9g8r+F//WX3xXxPZhmr3O11j/QWqdorVOBB4CPLi/aogUOfgbHNzkvTMqMbv2NMym3bMQJNIA79zU+LsizdzyRLrADOHO2Pa/vzYKT26CixJ7XDwMXTVBGmNVPOTNMygyljLNut2zE2fSC0QC+62mIinHPuN2o/g7g9t3sGYPHB7o2ot+ggyrcWutPtNYuu653oBOfw74PwTcPYuPtHk3LeLLg7CHnB07VbwAPuNHYtu+WKwU3WvtH+3cAp2SCiorolU9yxm2H7KchrhNkOjRMygxvlvGn089ed7xpJBpO+rb/SiHL2DxUXWn3yCLPxWKjaW33DuD4LtBzeETPc0vhDrezh2Hrq/ZeSlrBDYFTlzSA7zA+5/FBdbmxbV9YK9AAdsIOYK8PCvKhptrukYSEFO5wW/sH+y8lrVAXOLXW7pE07sCqKxvAbrlScBun7QD2ZEHlBSjcZvdIQkIKdzhdLIYNK+Cqr0LXFLtH03reLDi+BSpL7R5Jw7LnG6tf6jeAO/eGhH7OfsNxo0AD2Ck7gOsvWY1AUrjDKW+RcSk50QGXklbwZIGucWbgVKABnNVAA9g7wZjisWNXXyRy4g7grh7okiyFW7RS3aXkLdBrmN2jsYbHHzjlxO59oAGcMfPKr3l9UFoExfvDP65IdHkD2Anq7pHqwJ9NC0jhDpeNK+HiaedcSlqhfTdIGuq87v2ZQ003gD3+ee4I/aUOq4YawE7hzYLzR+HsEbtHYjkp3OFQU20E7iRnQL+Jdo/GWl4fHMkzAqecYl0za4mThkB814i9jA6rhhrATuHxz3NH4Bu0FO5wCFxKXu2gS0mreLKg4pwROOUEdQ3grzXeAI6KgpTxEfkLHXYNNYCdotcIY+t9BL5BS+EONa2NH+7EATD4drtHYz2nde/rGsBPNP04r++LwCnRMk01gJ0gOgZSMpw3lWcBKdyhduBT515KWiEQOOWEwl150XwDODDPHcF5FiGXPb/xBrBTRGjglBTuUAtcSo58wO6RhEZgG7kTzmo2vWC+AZw8zh84Jeu5W+TMIdj6mvN3AEdo4JQU7lA6vgX2fWQ0yZx4KWkVb5axlf/8cfvGEGgAp2SaawDHdYA+o5y5lNEN6hrA37B7JE2rC5xywImFhaRwh9KaJtYSR5K65XU2/nLUrSV+0nwD2JMFxzZI4FSwLmkAJ9s9mqbFdzG24EvhFqZccimZYPdoQqvPSONejnadvWoN2U8F3wD2+gOnjm8O3dgiUe5z5hrATuHJirjAKSncobL2D8YlmtMvJa0QHWvMGdt1xn3gU6P4BtsAdsKVgttUXoTcBe7aAezNgqpSOLnV7pFYRgp3KJSeNi4lR7rgUtIqXp99gVPZ86FTr+AbwJ17QbfUiLuMDqlAA/jqb9s9EvMicCOOFO5QyFsE1WXuuZS0gl2BU4EGcEvvJuTJksAps+o3gL0T7B6NeQmRFzglhdtqgUvJQbdCz6F2jyZ8PJmACv8vR2vXEkvglHk73gi+AewUERY4JYXbanVriV10KWmF9t2MN6pwFu4zB2Hb31rXAJbAKXPcvgM4wgKnpHBbqabayCVOGf/FnVbaEo/P2OhQWxOe11v7x9Y3gOsCp2QjTpPqGsDfcucOYI/DohlaSQq3lXa8Ydz53I2XklbwZkHFeSgMQ+CUVQ3gqCjjl1o24jStrgF8v90jaZleI4wptQhZQSSF2yqBXOLEge68lLRCXfc+DL8cec9Z1wD2+ODULgmcaszxza1rADtBIHAqQt6gpXBbZf8nxp3DJ33LOItri7qlGmdlof7lqLwIORY2gAPTWkdyW3+sSJT9NMR1dv8OYE+WcfPg8vN2j6TVmq0wSql4pVSuUmqzUmqbUurn4RiY67j9UtIKdbeLCvEZ96YXoKzYugZw37FG4FSEXEZbKtAAznjU/TuAvZETOGXm1LACuF5rPQoYDdyqlGqDnbcmHN8M+z82wqRi2tk9GnvVBU4dC83xQ9EAlsCpxgUawL5G7ibkJoHAqQhYQdRs4daGC/4PY/3/yW6F+rLnG5eS42bYPRL7BYppqLr3218PTQNYAqeuFGk7gNt1jpjAKVOTsUqpaKXUJqAQeF9r7f63LKvUXUrOcNWl5HtbT/BCziHrD9x7JMR2CM1ZTd1a4hA0gCVw6kp1DeBv2T0S60RI4JSpwq21rtFajwZSgPFKqRGXP0YpNUcpla+Uyi8qKrJ6nM619g+gohu/Ma0D/f3z43zjhfX8+PWtHD590dqDBwKnQnFWE8oGcGAjjqznNtQ1gG+DnkPsHo116gKnPrd7JK0S1E+/1vos8AlwawNfW6i1ztBaZyQlJVk0PIcrPQ0bnjcakl362j0aU1bvOcWTf97EiOSuxERF8dxnIdjq7fEZ9yOsuND8Y4MRygZw517GbdgiYP7TEnUNYBN3E3KTuqk8d/87m1lVkqSUSvD/vT1wI7Az1ANzhdyFxqXkJHdcSm46cpY5z+eTltSR52f6uGdMMi/nH+H0hQprX8gbgsCpcDSAvVnGlUJbD5yK5B3AXVOgS4rrVxCZOePuA3yslNoC5GHMcb8d2mG5QGWpUbgH3w5Jg+0eTbP2Fpbw6NJcEjvFsWLmeLp2iOWxyWlUVNeyfK3Fc90p/sApK89eAw3gUK4l9vjg4ikJnApVA9gpvP6dsi5+gzazqmSL1nqM1nqk1nqE1voX4RiY4210z6VkwZmLPLIol9joKFbO8tGzi7H7bUDPTtw0rBcr1h7kYqWFzZr2CdYGTtVvAMd3teaYDQn1ihg3CGUD2Ck8WVByDM65N3CqjW7xa6Waalj7e+MMzeGXkqcuVDBtcS6lldWsmDmefokdL/n6vCnpnL1YxV/yLP4h9viMnYhWBE6FqwHcY7DxxuDyy+hWaQs7gL2BwCn3znNH6L9MiG1/3dhk4vCz7ZLyKh5dmsuxc2UseTSToX26XPGYcf26kZnajUWfHaCqpta6F/dOgMoSKNzeuuOEswEsgVP+BnDvyN4B3HO46wOnpHAHK3Bj2h6DjKVSDlVeVcPs5fnsPF7Cnx4eR2Zq90YfO3dyOkfPlvHu58etG4DXohjNcDeA23Lg1LFNbWMHcAQETknhDtb+j42lbhOdeylZXVPLN1/cSM6BYv7va6O4bkjPJh9//ZCeDOzZiWc/3Y+2qmGT0M84c2tNg9KOBrC3Dd9YYU0gTKoN7ACuC5w6Z/dIWsSZlcfJ6i4lv2b3SBpUW6v53quf88GOk/z8ruHcPbr5rcpRUYo5k9PYcfw8q/acsmYgSn3RvW8pOxrAfcdCVGzba1CGqwHsFC4PnJLCHYxjG43mjUMvJbXW/Ne7O3h1QwHfvnEg0yemmn7u3aOT6dWlHQs+3WfdgDxZcK6FgVN2NYADgVNt7YzbhTuAWyUQOOXS6RIp3MHIfhradXHspeQfP9nHotUHmD6hH0/eMDCo58bFRDHr6v6s2XeaLQVnrRlQa+a56xrANty705sFRzdAtcUbk5yq9JTRAB7lnh3Ardaus3FXHJc2KKVwm1V8wCgmDr2UfCHnEP/7j13cPbovP/3ScFQLNk48ON5L5/gYFnxq0QaUlgZOXdIAviJdIfQ8PqipaDuBU7kRGCZlhjcLCta7MnBKCrdZgUtJB+YSv7PlOD96fSvXDU7iN18dRVRUy3a7dY6P5ZGsfvx963EOnipt/cDqAqeCDG6yuwHcljbiuGwHsKU8PtcGTknhNqP0FGxc6b+U7GP3aC6xancR3/7LRsZ5u/HHh8cRG926f9IZE1OJiYpi0WqLzrq9WXBia3CBU9nzoXMf+xrAnXq2ncApF+0AtpyLA6ekcJsRWEs80Vk/3BsOn2Hu8+tJT+rE4kczaR8X3epj9uwSz71jk/lrfgGnrAif8gQCp/LNPd4pDeC2EDhV1wDOcvwO4JBwceCUFO7m1F1K3gFJg+weTZ3dJ0uYsTSPnl3asWLWeLq2j7Xs2I9NTqOyppblaw62/mAef+CU2bOaQAN43KOtf+3WCAROnbZwlY3TuGQHcEh5fa58g5bC3ZyNK6HsjKN+uI8UX2Tq4hzaxUTx/EwfPTvHW3r89KRO3DysFyvWHqK0opWNm/iu0HOYubMaJzWAvROMP114NmaK3Q1gp/BkQclx4w3MRaRwN6WmCtY847+U9Nk9GgCKSiqYujiHssoaVswajzexQ0heZ+6UdM6VWRQ+5fXBkbzmA6fW/sG427oTGsA9BkF8QuQ2KO1uADuFS3fKtuF/MRO2vW5sILnahrXEDThfXsX0JbmcOF/O0hmZDOl9ZWiUVcZ6uzE+tTuLV1sQPuXJaj5wKtAAHumQBnAgcMplv9CmrX7K3gawU/Qabmzzd9kbtBTuxgRyiXsMhoG32D0aIzRqWT67T5bw7CPjGNev8dAoq8y7No2jZ8t4e0sLdj7WZ2YjTl0D2EFrib0+OLU78gKnjm2EA5/a3wB2gqhoI3DKZW/QUrgbs+8jY32nA3KJq2pq+eaLG8g7VMxv7x/NtYObDo2yyrWDejKoVycWtDZ8KqGfcXbXWOF2aAO47gbCLvulbpZTGsBO4c2Ck+4KnJLC3ZjAWuKrvmrrMGprNd97ZQsf7CjkF3eP4K5R4duSbIRPpbPzRAmf7i5q+YGUanrawYENYACSIzBwykkNYKfw+ADtqsApKdwNqbuU/Iatl5Jaa375znZe23iUf71pEFOz+oV9DHeN6kufrvE829rwKW+Wcauoc0cv/XygAeyd4JgGcJ3Y9pEXOLX2Gec0gJ0iJcN1gVNSuBuSPd8Rl5LPfLSXpdkHmTEplSeuH2DLGALhU+v2F7PpSCvCpzz+onz58rpAA9hpZ9sBkRQ45bQGsFO4MHBKCvflivfD9jeMu4nHh27VRnOeX3eI/3t/N/eMSebHdwxrUWiUVR7wh08tXNWKs+7eVxmBU/XPagIN4KQhjmgAN8ibZQROHdtk90haL3chVJc7qwHsFHWBU1V2j8QUKdyXC6wltjGX+M3Nx/jJG1u5YUhP/t99I1scGmWVTu1imJrVj79vPdHy8KlA4FT9s5pAA9jJa4kbu1JwG6c2gJ3Cm2UETp1wR+CUQ39bbHKhyB8m9QB07m3LED7ZVci//mUTmf2684eHx7Y6NMoqj05KJTY6ioWftSJ8yptl/GJUlBgfZz/liAZwkzr1hO5prpr/bFCgAeyQPQmO47IVRM6oCk6Ru9CYy7TpUnL9oWK+vnIDA3t1ZtGjGcTHtj40yio9O8fzlbEpvLK+gKKSFs73erL8t4vKN+aND6zyN4DjrB2s1TxZxi+0y/Is6tRvAHvG2z0aZ+qaDF09rllB1GzhVkp5lFIfK6V2KKW2KaUc2kVqpYoLRuEecgf0CO7uMVbYeeI8M5bm0atLO1bMHE+XeOtCo6zy2DX9qWpN+FQgcOpIjnFj2nZdbW8Am+J1eeCU0xvAThFYsuqCN2gzZ9zVwHe01kOBLOBxpdSw0A7LBhtXQvlZW26Vdfj0RaYtzqV9XDTPz/KR1NmZu9nSkjpxy7DerFh7sGXhU/FdjS3GW18zGsCZ9jaATau7jHbH2dgl3NAAdgqvewKnmi3cWuvjWusN/r+XADuA5m8d7iY1Vcb6Vu9E/1lh+BSWlDN1SQ4V1bWsmOnD0z00oVFWmTsljfPl1byU28Ifbo8PTu3yryWeZ+3gQsXNgVNuaAA7RV0j2vnz3DHBPFgplQqMAULznb34gLH0KtwqSozNIbf/Jqwve66simmLcyk8X8HK2T4G9+4c1tdviTHebvj6G+FT0yemBt889WZB/uKwN4C11ixefYBuHeL4yriU4J7s5sApmxrAa/ad4rlV+6muDf+0Q5RSfOuGgYzr1y24J9YPnHJ4+Jbpwq2U6gS8Cnxba32+ga/PAeYAeL3elo2m8oKxztQOox6EgTeH7eXKKmuYvTyPfUUXWDQ9M/gfMhvNm5LOjGV5vLX5GPeODbIIDrwJht0N13w3NINrxO8/2stv399Nh7hobhjak4QOQTZEvVmw5x9Qeho6JoZmkFYLNIBv+mVYG8DrDxUza1k+XdrH0DehfdheN2Bv4QV+/e4OXvn6xOCe6KLAKVOFWykVi1G0X9Bav9bQY7TWC4GFABkZGS17m3307RY9zW2qamp5/MUN5B86w9MPjGHKoCS7hxSUawcnMbhXZxZ8up97xiQHtzmofTf42orQDa4Bz689yG/f3801A3vw2Z5TrFx3iG9eH2QDun5u85DbLR9jSNjQAK7fZP/rvIm29GuWZR/gZ29tJ/9gMRmpQaZoeifAJ7+GsrPQPiE0A7SAmVUlClgM7NBa/zb0Q4pstbWaf/vrZj7aWcgv7x7Bl8IYGmUVpRRzp6Sx62QJn+xqRfhUGLyx6Sg/eXMbNw7tyZJHM7l2cBLL1hykvKqZmzpcru8YI3DKLQ3KwA7gMDaAndJk/1qmh4QOsSxY1YI9B95A4JTJe6TaxMwE5SRgKnC9UmqT/z+XnHI4i9aaX7y9ndc3HeO7Nw/iERtCo6zypVF96WtF+FQIfbKrkO+8vJnM1O4885CxmWnelHROXajklfUFwR0stj30He2ejThrnglrA9hJTfYOcTFMm5DK+9tPsrewJLgnJ2eAinb8G7SZVSWrtdZKaz1Saz3a/9+74RhcpHn6w70sW3OQWVf35/Hr7AmNskpsdBQzr+5PzoFiNh4+Y/dwrrD+UDHzVq5ncO/OLJr+xWYmX//ujPIk8Nxn+6kJtnHm8RnJkU4PnLpQBJteCFsD+FxZFdOX5FF4voKlMzId0WSfPqEf8bFRLAz2rLtdJ+g9wvEriGR9UJgsX3OQ332wm6+MTeGHtw+1NTTKKg+M99IlPoYFn7ZiG3wIBOZZ+3Rtz/LLNjMppZg3OY1Dpy/yj20ngjuwWwKnwrgDONBk31tYwoKp4xjrdUaTPbFTO76W4eFvG49y8nyQCx48WXDU2YFTUrjD4I1NR/npm9u4cWgv/ucrV9keGmWVTu2MS9J/bD/B/qILdg8HMOZZpy7OpUNcDCs/5MZVAAAYeElEQVRmjqdHpyvnWW8e3pv+PTqy4NN9wd3Zxw2BU2HcAVy/yf67+0cz2WFN9tlXp1FTq1mSfSC4J3p9UHXR0YFTUrhD7OOdxjyrr393nnloDDEOCY2ySmAt93OfBfnLEQKFJeU8sjiHqppanp81vtF51ugoxWPXpLG54Bxr9582/wJuCJza+HxYdgDXb7L/55dHcOdI5zXZvYkduP2qPry47jDny4M4e3ZB4FRkVRGHyT9YzNdfMOZZn5vurNAoqyR1bsd941J4dUMBhSU2rcHni81Mpy5UsPTRTAb2anqe9d6xyfToFBf8NI8nyzjjdmKeRU2VEUsc4h3A9Zvs/3bLYB72ObfJPm9KOiUV1byYE8ROXxcETknhDpEdx88zc1nD86yR5rFr0qiqqWVZ9kFbXr+ssoZZy4zNTAumjmOMiXnW+NhoZkzqz6e7i9hx/Ir9ZI3zZsHF03B6bytGHCLb/mbsAA5xmFT9Jvs3rk0P6Wu11ojkrlw9oAdLVh+gojqIJaAen1G4nfgGjRTukDh0upRpS4x51udnNTzPGkn69+jIbSN68/y6Q1xoSfhUK1TV1PKNF9az/vAZnrp/DNcMND/P+oivHx3jooNbeRDYiOO0s7G6MKmhId0B7MYm+9wpaRSWVPDGxmPmn+TNggsn4Oyh0A2sFaRwW6zwfDlTF+fWzbOmdHN2aJRV5k5Op6S8mj+3NHyqBWprNd/962Y+3lXEr758FXeMDO4+il07xPLgeC9vbj5GwZmL5p6UONDY/em0BuXeD+HkVpgUujAptzbZrx7Qg+F9u/Dsqn3Uml0CWvcG7cx5bincFjp3sYppS4x51mUzxjc7zxpJRnkSyEozwqcqq2tD/npaa37+1jbe8M+zPuRrWT7OzKv7o4DFq002VwOBU077hc5+Cjr3hRH3heTwbm6yGzt909lfVMoHO06ae1LPYcYNw532Bu3nnv/7DnexspqZ/tCohVMzGO1xbs5BqMydks7xc+W8uTmIS9IWmv/hHpavPcRj17RunrVvQnvuGt2XP+ce4UxppbkneXxweo8ROOUER9fDwc9gQmjuJhQJTfbbR/QmpVt789vgA4FTTnuD9pPCbYHK6lq+vnIDGw+fYf4DY7h6YA+7h2SLawclMaR3ZxYGc0naAsuyD/DUB3u4b1wK/2HBPOvcyemUVdWwcp3J+Uyvw5aLZc8PWZhUpDTZY6KjeOyaNNYfOkP+wWJzT/JkQeF2I3DKYaRwt1JgnvXT3UX86p6ruP2q4OZZI0kgfGr3yQt8srswJK/x+saj/Oyt7dw0rBf/fe9VljTHBvfuzPVDepoPn3JS4NTpfbD9TcicBe2snZqLtCb7VzNS6NYh1ny+joMDp6Rwt4LWmp+9tY03Nx/j328dzIPjW5hDHkHuHNmX5IT2PBuCbfAf7yzku3/dTFZad37/oLXzrHMnp3G6tJK/mgmfqgucckDhXvsMRMdaHiYViU32DnExTJ+Yygc7Ctlz0kT4lIMDp6Rwt8JTH+xhxdpDzJmcxtenOHs9a7jERkcx6+r+5B4oZoOF4VN5B43QqCF9OvPcNOvnWcf3785oTwLPrTIZPuXNMgKnquzbdMSFQtj4gnETkM69LDtsJDfZp01INR8+5eDAKSncLbQ0+wDzP9zDV8el8IPbhrhiPWu43J/poWv7WBZYFPm6/Zgxz5qc0J5lM8bTOQTzrEop5k1J43DxRd7baiJ8ypMFNZVw3MbAqZwFxhgsDJMKNNn3F5VGZJO9e8c47s/w8Pqmo5w4Z+JN1zvBmCpxWOCUFO4WeH3jUX7+1nZuHtaLX1s0zxpJOraLYdqEfvxz+0n2tTJ8KjDP2qldDM/P9oV0nvWmYb1J69GRZ82ETwUCp+w6G6u4AHnPwdA7oYc1EcGXNtlHR2yTffY1adRqzIVPeXxQXQYntoR+YEGQwh2kj3ae5Dt/3cyEtESetnieNZJMn5hKXHQUz7XkLiR+J88boVE1tcY8a3KI718YHaV4bHIanx89x9p9zSz165QE3dPtW1myYQWUn7MsTOryJvttEdxk93TvwB1X9eHFnMOcK2vmTNqhG3Gk6gQh90AxX1+5gWF9urBw2jhXrmcNlx6d2vHVjBRe23CUwmDzkPHPsy7O5fSFSpbNGM+AnuGZZ71nTDI9OrXjWTNvON4so3CHO88iECbVb5Kx1riV6jfZv3frkDbRZJ8zOY0LZsKnuvSFrl7HNSilcJu0/dh5Zi3PI7lbe5bNyAzJPGukmX11GtW1tSxdczCo512srGbGslwOnCrluWkZjArjPKsRPpXKqt1FbD/WTPiUx2dP4NTWV+F8gWVn27+r12SfNyXNkmM63YjkrlwzsAdLsk2ET3n9O2UdFDglhduEg6fqzbPO8pHo8vWs4ZLaoyO3jejDynWHKDGZhxyYZ9105CxPPziaSQPCP8/6SJYRPrVgVTPNVTsCpwJhUj2HwcCbWn24pdkHePrDPXwto+012edOTqeopIK/bTja9AM9PscFTknhbka451kjzdwpaZSUV/OSifCpmlrNd/zzrP91z1XcOsKeedau7WN5yOfl7S3HOVLcRPhUIHAqnIV77wfGbr6J34JWFtm/bSyoa7L/1z1tr8k+aUAiI5K7sHDV/qZ3+jpwnlsKdxPOXqxk2uJczpSGd541koxMSWBiemKz4VNaa3725jbe2nyM7982hAdsnmc1FT4VFfXFjRXCJXs+dEmGEV9p1WE+2nmS7/51S5tusiulmDs5nf2nSnm/qfApBwZOtb1/LZMuVlYzc1keB06VsjDM86yRZu6UdE6er+CNTY1fkv7u/d08v+4QcyenMc8Bm5n6dG3P3aOT+UteM+FTXp8xx116KvSDKgiEST3eqjApabJ/4bYRvfF0b9/0EtCoaEjJdNRGHCncDaisrmWezfOskWTywB7+8KmGL0mXrD7A0x/t5f4MD9+/bYgNI2zY3ClplFXVsGJtE3Ob4bw/YfZTEN8Vxk5r8SG2HTvHrGXSZA8IhE9tPHyWvINN7PT1ZkHhDscETknhvkxNreZfX97Eqt1F/Ppe++ZZI4mxKzGdPYUX+GjnpeFTr20o4Bdvb+eW4b341T0jHDXPOqhXZ24Y0pPlaw9SVtnIyoO+YyA6LvRnY6f3wY63IHN2i8OkDp4qZfqSPDrFS5O9vq+O89C9Y1zTO309gcCpvLCNqynNFm6l1BKlVKFSams4BmQnrTU/fXMrb285zg9uG8L9mZG/njVc7hjZh+SE9pes1Phg+0n+7ZUtTExPZP4DzpxnnTslneLSSl5Zf6ThB8TGQ5/RoT/jXvN74w2ihWFSlzbZfdJkr6d9XDTTJ6Ty4c5CdjcWPpXiD5xyyHSJmd+UZcCtIR6HI/zu/d2sXHeYuVPSmOuAedZIEhsdxexr+pN38AzrDxWTs/80j7+4geF9u7AwBKFRVslM7cYYbwLPfXaA6ppGmqteX2gDp0pOwqYXYfRD0Kln0E8/e7GSqYtzOFNayfKZ4xnQs1MIBulu0yb0o31sE/cfjesIva9yTAZ7s4Vba70KMJk87l6XzLPe6px51khyf6aHhA6x/PLtHcxenu+fZx1Pp3Yxdg+tUYFpnsPFF/l7Y+FToQ6cyg2EST0R9FONzUx5HDx1keemZTAyRZrsDenWMY77Mz28sekox8+VNfwgb5ZjAqcc9RvzwMK1VIThfoWXq9Ww+chZbh3e23HzrJGkQ1wM0yak8vSHe+jbNZ6Vs3x072j9rbasdtPQXqT16MiCVfu4c2SfK38+6gdOBdb8WqWiBPIWwdAvQWJwV4GBJvvmI2f548PjmChN9ibNuro/z687xJLVB/jhHcOufIDHBznPGoFTyePCP8B6LCvcSqk5wBwAr7dlc8Od2sUQGx3+wg3GpdIP7xjqyHnWSDJzUirFpRXMmNSfvi6ZZ42KUsyZnMb3X/ucNftOX7nKqFMSJA4IzfxnC8Okamo1/+Jvsv/PV67i1hG9rR9bhPF078CdI43wqW9eP5Cu7S9bcVN/I06kFG6t9UJgIUBGRkaLNvUvmp5p1XCEQyV0iOM/v3yV3cMI2pfHJPN/7+/m2U/3Nbw81JMFu941tqRbdcVWXWmESaVeAynmC4XWmh+/sZV3pMketDmT03hj0zFeyDnEN669LC43EDh1eK1xY2YbyemlECbEx0Yzc1J/Pttziq1Hz135AK8Pyorh1B7rXnTrq3D+KEx6Mqin/d8/d/NizmHmTUmXJnuQhvf1h0+tbuT+o3YlQl7GzHLAl4C1wGClVIFSalbohyWE8zzk89KpXUzDKw/qNuJYNF1SP0xqwI2mn7bos/088/FeHsj08L1bB1szljZm3pR0Tl2o4G8bG9jp6/XBhZNw5mDYx1WfmVUlD2qt+2itY7XWKVrrxeEYmBBOEwifeufzBsKnegyE9t2tCyLa808o2mGcbZucenl1fQH/+c4Of5O97YVGWWVieiJXJXdt+P6j4dwp2wSZKhEiCDMmpRKlGgifUspYdWDVGXf2fOiSYjpM6v3tJ/n3V7cwaUAi8x8cTXSUFO2WUkoxd0qaET61/bLwqZ5DjcApmzfiSOEWIgiB8Kk/5x2m+PLwKasCp47kwaFsI0wquvkskXX+zUwj+nZhwdQM2sU4czOTm9w6vDfe7h2uDJ8KBE7JGbcQ7jJ3chrlVbWsWHvw0i9YdRm9Zj7EJ5gKk9p69Byzl+fj7d6BpQ7fzOQmRvhUfzYdOUvugcv2HzogcEoKtxBBGtirMzcO7cnyNZeFT9UFTq1t+cFP7YEdb/vDpJremr6/6ALTl+TStX0sz88a74rNTG7y1QwPiR3jWHB5M9oBgVNSuIVogblT0jlzsYq/1g+fio03indrGpR1YVJzm3zY8XNlTF2cC8CKWePp09Udm5ncJD42mukTU/loZyG7TtQLn3JA4JQUbiFaIKNfN8Z6E1i4av+l4VMen5FZ0pLAqZKTsPklGPNwk2FSZ0qNOzOdK6ti+czxpCdJaFSoTM0ywqcuuf9oXEfoM1IKtxBuEwifKjhTxrv1w6e8/sCpYxuDP2jOs0aA0YRvNvqQ0opqHl2Wx6FiIzRqRHLXFoxemBUIn3pz0zGOna0XPuXJgqPrbQucksItRAvdOLQXaUkdWVB/5UEgcCrYZYHl5yFvMQy7q9EwqYrqGuY+v56tR8/xzINjmJCe2IrRC7NmX9MfjZEgWsfrg+oyOL7FljFJ4RaihaKiFHMnp7Ht2Hmy9542Ptmxhz9wKsh57g3LoeJco9vba2o1//KXTazee4r/+cpIbh4uoVHhktKtA18a2YeXcg9z7qL/DNvqnbJBksItRCt8eUwyPTu349n6t73yBJlnUV0Ja/9ohEk1kDqnteZHr3/Ou5+f4Ed3DOW+cSkWjV6YNWdyOqWVNazM8d9/tEsfSPDaNs8thVuIVmgXE83Mq/uzem+98KlgA6e2vgIlxxqNbv3ff+zipdwjPH5dOrOvSbNo5CIYw/p2YfKgJJZmH/gifCrYN2gLSeEWopUC4VN1630Dl9Fm1nPX1vrDpIbDgBuu+PJzq/bzx0/28ZDPy3dvltAoO82bksapC5W8tsEfPmVj4JQUbiFaqUt8LA/7vLyz5ZgRPhUInDKzg3LPP6FoZ4NhUi/nH+FX7+7gjqv68Mu75c5MdpuQlsjIlK4895k/fMrGwCkp3EJYYMak/kRHKZ77bL9RgL1Z5uY/s+dDVw+MuPeST/9j2wm+/+oWrhnYg9/eP0pCoxxAKcXcyekcOFXKP7ed8AdOdW3dTtkWksIthAV6d43nnjHJvJx/hNMXKoxlgcX74EJR4086kguH11wRJrVm3ymeeGkjI1MSePaRcRIa5SC3juhNv0R/+JSKAk+mdVG+QZDCLYRF5tSFTx364v6ETV1GZ/vDpMZMrfvU5wXnmLNiPf26d2Dpo5l0lNAoR4mOUjx2TRqbC86Rc6DYmC4p2gFlZ8I6DincQlhkQM/O3Di0FyvWHuRijxFG5khj63yLdsPOd2D8Y3VhUvuKLjB9aSA0ykc3CY1ypPvGpRjhU5/uMxqUYETxhpEUbiEsNG9KGmcuVvHyxqKmA6fW/h5i2sF4I0zq2Nkypi7KIUrBytk+eneND+OoRTDiY6N5dGIqH+8qYlfMICNwKswbcaRwC2GhjNTuZPTrxnOfHaA2ZbyRWVJVdumDSk7A5j/D6IehUxLFpZVMXZxDSXk1y2aMp3+PjvYMXpg2dUI/OsRFs2DNCX/gVHjnuaVwC2GxuVPSOXq2jLzaQVBbdWXg1Lo/QW01THicCxXVPLo0l4IzZSyaLqFRbpHQIY4HMr28ufkYF3pmhD1wSgq3EBa7YUhP0pM68rud3Y1P1F8WWH4e8pfA0Luo6JrKnBX5bDt2nj88NBZfmoRGucksf/jUe+f7hT1wSgq3EBYzwqfSWXdScbFL2qUrS9Yvg4rzVE/4Fk++tIk1+07zv/eN5MZhvWwbr2iZ5IT23DWqL7/f63/DDeM8txRuIULg7jF96dWlHbnVA43CXVtrhEmt+yO6/2R+mBPLe9tO8JM7h3HvWAmNcqs5k9M4VNmVc/HJYd2II4VbiBBoFxPNzEn9eedcP2ON7+k98PnLUHKcv8Tdy1/yj/DE9QOYeXV/u4cqWmFony5cOziJ7Ip09OHwBU6ZKtxKqVuVUruUUnuVUt8P9aCEiAQP+rzsjBlmfHB4LWQ/zalOg/j+5iQeyfLyrzcNsneAwhJzJ6eTXTkAVVoIZw40/wQLNFu4lVLRwB+A24BhwINKqWGhHpgQbtclPpaJPh+ndWeqP/0NnNrFL4pv4s6Rffn5XRIaFSmy0rpTkmTkqNceCs88t5kz7vHAXq31fq11JfBn4O7QDkuIyDDz6jQ26kHEnD9Cge5BSfqd/PZroyU0KoIopbj1+us4rztQsOXjsLymmcKdDByp93GB/3NCiGb06hJPZZ9MAN7r/BX+MHU8cTHSWoo0t4zoy46YIejDOV/cfzSEzPwENXRqcMXIlFJzlFL5Sqn8oqImEtGEaGMy7voGub0e4L7H/oMOcRIaFYmioxTRw++iMGE0FVXVIX891dy7g1JqAvAzrfUt/o9/AKC1/nVjz8nIyND5+flWjlMIISKaUmq91jrDzGPNnHHnAQOVUv2VUnHAA8CbrRmgEEKIlmv2uk1rXa2U+ibwDyAaWKK13hbykQkhhGiQqQk3rfW7wLshHosQQggTpL0thBAuI4VbCCFcRgq3EEK4jBRuIYRwGSncQgjhMs1uwGnRQZUqAg618Ok9gFMWDscN5HuOfG3t+wX5noPVT2udZOaBISncraGUyje7eyhSyPcc+dra9wvyPYeSTJUIIYTLSOEWQgiXcWLhXmj3AGwg33Pka2vfL8j3HDKOm+MWQgjRNCeecQshhGiCYwp3W7shsVLKo5T6WCm1Qym1TSn1pN1jChelVLRSaqNS6m27xxIOSqkEpdQrSqmd/n/vCXaPKdSUUv/i/7neqpR6SSkVb/eYrKaUWqKUKlRKba33ue5KqfeVUnv8f3YLxWs7onC30RsSVwPf0VoPBbKAx9vA9xzwJLDD7kGE0XzgPa31EGAUEf69K6WSgW8BGVrrERhx0A/YO6qQWAbcetnnvg98qLUeCHzo/9hyjijctMEbEmutj2utN/j/XoLxyxzx9/JUSqUAdwCL7B5LOCilugCTgcUAWutKrfVZe0cVFjFAe6VUDNABOGbzeCyntV4FFF/26buB5f6/Lwe+HIrXdkrhbtM3JFZKpQJjgBx7RxIWTwH/DtTaPZAwSQOKgKX+6aFFSqmOdg8qlLTWR4HfAIeB48A5rfU/7R1V2PTSWh8H4+QM6BmKF3FK4TZ1Q+JIpJTqBLwKfFtrfd7u8YSSUupOoFBrvd7usYRRDDAW+JPWegxQSogun53CP697N9Af6At0VEo9Yu+oIotTCncB4Kn3cQoReGl1OaVULEbRfkFr/Zrd4wmDScBdSqmDGNNh1yulVto7pJArAAq01oGrqVcwCnkkuxE4oLUu0lpXAa8BE20eU7icVEr1AfD/WRiKF3FK4W5zNyRWSimMec8dWuvf2j2ecNBa/0BrnaK1TsX4N/5Iax3RZ2Ja6xPAEaXUYP+nbgC22zikcDgMZCmlOvh/zm8gwhuy9bwJTPf/fTrwRihexNQ9J0Otjd6QeBIwFfhcKbXJ/7n/8N/fU0SWJ4AX/Ccl+4EZNo8npLTWOUqpV4ANGKunNhKBuyiVUi8B1wI9lFIFwE+B/wZeVkrNwngD+2pIXlt2TgohhLs4ZapECCGESVK4hRDCZaRwCyGEy0jhFkIIl5HCLYQQLiOFWwghXEYKtxBCuIwUbiGEcJn/D1031XVRErIaAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "x = np.array([1, 1, 2, 3, 1, 2, 0, 1, 2, 2, 0])\n", "y = 2*x\n", "plt.plot(x, label='x')\n", "plt.plot(y, label='y')\n", "plt.legend()\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhoAAAIYCAYAAAArcUUtAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzs3Xd4leXhxvH7OXtmhxEChJCwCSBhOpAZ3IrWatVWO2jrrKsurLZatbXaal3V/rTa2tYqIO6AgBuRALJXwgwjhOyz1/P7I4CgKCFnPO855/5cl1cFwnnvHm3z5T3vOa+QUoKIiIgoHnSqBxAREVHqYmgQERFR3DA0iIiIKG4YGkRERBQ3DA0iIiKKG4YGERERxQ1Dg4iIiOKGoUFERERxw9AgIiKiuGFoEBERUdwYEnmwvLw8WVRUlMhDEhFRHC1fvvyAlDJf9Q7Srk6HhhDCAuAjAOaDj/OalPKe7/o9RUVFqKqq6uwhiYhIY4QQO1RvIG2L5oyGH8AkKaVLCGEE8IkQ4l0p5ecx2kZERERJrtOhIdtv++o6+EPjwb94K1giIiI6LKqLQYUQeiHElwD2A1ggpVwam1mULGQkiJbPr4Vr3aOqpxARkQZFFRpSyrCUcjiAQgCjhRBDvv41QoiZQogqIURVfX19NIcjjZGREJo/ugKejU+ibfmdCHv2qJ5EREQaE5O3t0opmwF8AGD6MX7tWSlluZSyPD+fFyanChkJo+XTH8O3/RXYB90IyBDcax9WPYuIiDSm06EhhMgXQmQd/HsrgCkANsZqGGmXlBG0LJkJb80/4RhxPzJGPwpr8eVwb/obwt79qucREZGGRHNGozuAxUKI1QCWof0ajbdiM4u0SkqJ1s+vgXfL83AM+w2cw+4CADjK7gQifrh5rQYRER0hmnedrAYwIoZbSOOklGj94lfwbHoG9iG3wTH83sO/ZsjsB0vR9+HZ+CQcQ26FzpKrbigREWkGP4KcOkRKibaqX8Oz4XHYB90I58gHIYQ46mscZXdChlxwb3hM0UoiItIahgYdl5QSbStnwb3uT7ANuAbOUY98IzIAwJg9BJbeM+Be/zgigRYFS4mISGsYGnRcrlX3wb36AVj7/QwZYx4/ZmQc4iibBRlsgXvDXxO4kIiItIqhQd/JtfohuL68B9aSK5E57hkI8d3/yhhzR8BceDbc6/+MSND1nV9LRESpj6FB38q17lG0rbgDluIfIHP8348bGYc4hs2C9DfCs+npOC8kIiKtY2jQMbk3PIG2ZTfD0vsiZJ3yIoRO3+Hfa8ofA1PBVLjX/gky5InjSiIi0jqGBn2DZ9OzaF16Hcw9z0PWhH9D6E78XdDOYXcj4tsPz+bn4rCQiIiSBUODjuLZ8gJalvwc5sIzkX36KxA6Y6cex9T1VJi6ToBr7R8hQ74YryQiomTB0KDDvDUvo+XTn8BUMBXZp8+G0JujejzHsLsR8eyBp/ofsRlIRERJh6FBAADv9lfR/MkPYep2OnImvQ5hsET9mKbuk2DMHwf3mocgI8EYrCQiomTD0CD4dryO5g8vhTF/PLInvwFhsMXkcYUQcAybhbB7B7w1/4zJYxIRUXJhaKQ536630fThxTDmjULO1HegMzpi+vjmHmfAmDsSrtUPQEZCMX1sIiLSPoZGGvPvrkTT4hkwZpchZ8q70BmdMT+GEAKOslkIt9XAu+2/MX98IiLSNoZGmvLvXYTGRefDkDkQOdPmQ2fOituxzL3OhSF7KFyrfw8pI3E7DhERaQ9DIw0F6j5G08JzYHD2RW7F+9CZc+J6PCF07Wc1WjbCt2N2XI9FRETawtBIM4H9S9C44Ezo7L2QU7EQOkteQo5r6X0h9Bn94Vp1P89qEBGlEYZGGgkcqELjgunQWbsht2Ih9NauCTu20OnhGHYXQk2r4d/1ZsKOS0REajE00kSwYSUa50+FzpyL3OmLoLcVJHyDtc+l0Dv7wrXqPkgpE358IiJKPIZGGgg2rUHD/KkQRidyKhZBb++pZIfQGeAYegeCDcvh312pZAMRESUWQyPFBZs3oLFyCoTejNyKxTA4i5Tusfa9Anp7L7hW86wGEVE6YGiksFDLZjRWTgIgkFuxCIaMvqonQehNsA+9DcH9nyGwb7HqOUREFGcMjRQVatuKhspJgAwjZ/oiGDL7q550mK3kx9BZu8O16j7VU4iIKM4YGiko5NqBxvcmQYa9yJn2PoxZg1RPOoowWOAY8msE9n2AQN0nqucQEVEcMTRSTNhdi8bKSYgEW5A7bQGMOWWqJx2Trf9M6Cxd4Fp1v+opREQURwyNFBL27EVD5WREfPXImVoJY+5Jqid9K2GwwT74Zvj3VCJQ/4XqOUREFCcMjRQR9u5HY+VkRDy7kTP1PZjyR6uedFy2Ab+EMOfAtZpnNYiIUlWnQ0MI0VMIsVgIsUEIsU4IcUMsh1HHRXwNaJw/BSHXdmRPeRumLuNVT+oQndEJ+6Bfwb/rTQQbvlQ9h4iI4iCaMxohADdLKQcCGAvgGiGEtq46TAMRfxMa5k9FqGUzcia/AXO3CaonnRD7gOsgjBk8q0FElKI6HRpSyr1SyhUH/74NwAYAPWI1jI4vEmhB44IKhJrXIXvS6zAXTFE96YTpzFmwD7wevh2zEWxap3oOERHFmCEWDyKEKAIwAsDSWDweHV8k2IbG989EsGElsifNgaVwuupJnWYf9Cu41/8FrjUPIPu0l1XPiYoMB9Cy5OcQeisyxz2leg4dJKVE27Kb4al5UfUUAIAxZzhypi2AELxMjlJf1KEhhHAAmA3gV1LK1mP8+kwAMwGgV69e0R6OAESCbjS9fzaC9UuRNeEVWHqeo3pSVHSWXNgGXA33uj8hNPxeGDJKVU/qFBkJofmjH8C3YzYAwNbvZzDmjlC8iqSUaF16PTwbn4C51wXQ29SeeI349sO3/X/w7ZgNa9H3lG4hSgQRzf0mhBBGAG8BqJRSPnq8ry8vL5dVVVWdPh4BMuRF48KzEdj3AbJOfRnW4ktUT4qJsLcO+1/rA2uf7yPrlBdUzzlhMhJG88eXw7ftv3CM+B3cax+BuWAysifOVj0trUkp0VZ1C9zrHoV98M1wlj8MIYTaTZEw6ucNgdCZkHfuyqQ/qyGEWC6lLFe9g7QrmnedCAD/B2BDRyKDoidDPjQtugCBvYuReco/UiYyAEBv7Qpbv5nw1vwTobbtquecECkjaPn0x/Bt+y+cIx+Cc9jdsA+8Dr4dcxBsWqt6XtqSUqJtxZ1wr3sUtoHXaSIyAEDo9HCU3YlQ02r4d72peg5R3EWT0icDuALAJCHElwf/OjNGu+hrZDiApg++B/+eSmSOfw62vleonhRzjiG3AkIP95qHVE/pMCkjaPns5/DWvATHiN/BMfQ2AO3XnQiDA67Vv1e8MH25vvwt3Gsegq3/L5Ax+jFNRMYh1j6XQu/sC9cq3sWYUl807zr5REoppJRlUsrhB/96J5bjqJ2MBNH04SXw176FjLFPw9bvJ6onxYXe3gO20p/AU/0Cwu5a1XOOq/21/+vg3fJ3OMpmwTns7sO/dui6E9+2VxBq2aRwZXpyrX4ArlW/hbX0x8gY+6SmIgMAhM4Ax9A7EGxYDv/uStVziOIquV8cTAPtFxheAf/OucgY/RjsA36helJc2YfeBsgIXGv/qHrKd5JSovWLG+HZ+BTsQ26FY8TvvvE19sE3A3oLXKsfVLAwfbnW/gltK+6CtfhyZI57VrPXQFj7XgG9vRdcq3lWg1KbNv8XSADaLxpr+eQq+La/Amf5w7APul71pLgzOHrDWvJDeDY/h7Bnn+o5xySlRNvy2+HZ8BhsA2+Ac+QfjvknZr21S/t1J1v/hVDbVgVL0497/eNoq7oVlqKLkXnKCxA6vepJ30roTbAPvR3B/Z8hsG+x6jlEccPQ0CgpI2hZ0v5NyjHifjiG3KJ6UsI4ht4BRAJwr3tE9ZRjcn15D9xr/whb/18iY/Sfv/O0/KHrTlxJdN1JsnJvfAatX9wAc68LkHXavyB0MfmYoLiylVwFnbU7XKvuUz2FKG4YGhokpUTr59fAu+V5OIb9Bs5hd6melFCGjBJY+/wAnk1PI+I7oHrOUdpW3Q/XqvtgLf0JMsY+cdzX/g9dd+Kt/gfCrp0JWpl+PFueR+vnv4S58GxkT/gvhM6oelKHCIMFjiG/RmDfBwjUfaJ6DlFcMDQ0pv21/1/Bs+kZ2IfeDsfwe1VPUsJRdhdkyAPX+j+rnnKYa80f4Vp5N6x9f4jM8R1/7b/9uhMJ19qH47wwPXlq/oWWT38Kc0EFsie+BqE3qZ50Qmz9Z0Jn6QLXKt7vh1ITQ0ND2j9c6FZ4NjwO+6Ab4TzpAc1dLZ8ohqwBsBR9D54Nf0XE36R6Dtzr/oK25bfB0ucSZJ78/AldYNh+3cmPDl53sjeOK9OPd9sraPnkRzB1n4jsSXMh9GbVk06YMNhgH3wz/HsqEaj/QvUcophjaGiElBJtK2fBve4R2AZcA+eoR9I2Mg5xlN0FGWyDe8PjSne4Nz6F1mU3wtL7QmSd+lKnLjB0DL0dkEG41/0pDgvTk2/HXDR/dBlMXU5G9qQ3IAxW1ZM6zTbglxDmHN7FmFISQ0MjXKvug3v1A7D2+xkyxjye9pEBAMacMph7ngf3+scQCXzjNjoJ4dn8d7R+fg3MPc9B1mn/7vRr/19dd/IMwr76GK9MP75db6Lpw+/DmDca2VPehs5oVz0pKjqjE/ZBv4J/15sINnypeg5RTDE0NMC1+iG4vrwH1pIrkTnuGc2+718F57C7IQNN8GxM/J1QPdUvouWzmTD3mI7s01+N+rX/9utOvHCv4yf2R8NX+x6aFl/UfgfUqe9CZ3SqnhQT9gHXQRgzeFaDUg6/oynmWvco2lbcAUvxD5A5/u+MjK8x5o2EuccZcK97BJGgO2HH9W79D1o+/TFM3Scje+KcmLz2/9V1J08g4m+Mwcr049+zEE2LL4AhazByplZCZ8pUPSlmdOYs2AdeD9+O2Qg2rVM9hyhm+F1NIfeGJ9C27GZYel+ErFNe1PSHC6nkGHY3Iv4D8Gz+W0KO590+G80fXwFTl1ORM3leTF/7bz+r4YJ7vdrrTpKRf99HaFx4DgwZpcidNh86c7bqSTHXfo8cO1xrHlA9hShmGBqKeDY9i9al18Hc8zxkTfh3Uny4kCqmLuNg6j4J7rUPQ4a8cT2Wb+cbaP7wEhjzxyJ7ylsQBltMH//wdScb1F13kowC+z9D0/tnwuAoQs6096Gz5KmeFBdf3SPnvwi1blE9hygmGBoKeLa8gJYlP4e58Exkn/5K0ny4kEqOsrsR8e6DZ8vzcTuGr/YdNH1wEYy5JyFnyjvQGR1xOU77dSfN8Gx8Ii6Pn2oC9V+gccEZ0NkKkFOxEHprF9WT4so++GZAZ4ZrNc9qUGpgaCSYt+ZltHz6E5gKpiL79NlJ+b5/FUzdJsDY5RS41v4BMhyI+eP79yxA06IZMGQPRc7U96AzZcT8GIccuu7Ete5RRIKuuB0nFQQbVqBxQQV0ljzkViyC3tZd9aS401u7tt8jp+afCLVtUz2HKGoMjQTybn8VzZ/8EKZupyNn0usQBovqSUlDCAHnsLsRce+Ct+bFmD62f+8HaFx4HgyZ/ZE7NTGv/TuG3Q3pb4BnU2KuO0lGwcbVaJg/FTpjBnIqFkFvL1Q9KWEO3SPHzXvkUApgaCSIb8dcNH94KYz545E9+Y2Yv/afDkwFU2HMGw3X6gchI8GYPGag7lM0LTwbBmcf5ExbAJ0lNyaPezzt151MTsh1J8ko2LwejfOnQOityJm+GAZHb9WTEurQPXI81S8g7K5VPYcoKgyNBPDteuvghwuNQs7U+L32n+qEEHAMm4Wwaxu8W/8d9eMF6pei8f0zoLP1QM60xL/27yibhYivDp4tf0/ocbUu1LIZjZWTAaFH7vTFMDiLVU9S4qt75PxR9RSiqDA04sy/uxJNiy+EMbsMOVNS58OFVDEXng1DznC4Vj8AGQl3+nGCB5ajcX4FdJYuB1/77xbDlR1z+LqTNX+ADPsTfnwtCrXWoKFyEiDDyKlYBENGqepJyrTfI+eHB++Rs0/1HKJOY2jEkX/vIjQuOh+GzIHImTYfOnOW6klJTwgBR9kshFs3w7f91U49RrBxVftr/+bs9siw94jxyo45fN2JZzc81f9QskFLQq4daKycBBn2IadiIYxZA1VPUs4x9A4gEuA9ciipMTTiJFD3MZoWngODsy9yK96HzpyjelLKsPRu/2RI1+r7IWXkhH5vsGkdGiunQBjs7RcYOnrFaWXHHLruxL3moZhdd5KMwu5aNL43EZFgK3KnLYAxe6jqSZrw1T1ynuY9cihpMTTiILB/CRoXnAmdvRdyKham7IcLqSKEDo6yuxBqXgffztc7/PtCLZvaX/vXGZE7fREMzj5xXNkx7ded3I2wazu8NS+rnqNE2LMXDZWTEPE3IGfafBhzR6iepCn2sjvb75Gz/i+qpxB1CkMjxgIHlqFxwXTorN2QW7EQemtX1ZNSkqXoYugzSuFadT+klMf9+lBrdftr/5DI1dhr/+bCs9qvO1kT3XUnySjsrWuPDM9e5Ex9D6a8UaonaY4xayAsRRfBs+GviPibVM8hOmEMjRgKNqxE4/xp0JlzkTt9EfS2AtWTUpbQ6eEYeidCjSvhr33nO7821LYdDZWTIMN+5FQshCFrQIJWdsxX151sgW/7K6rnJEzEdwCNlVMQce9E9pS3YeoyTvUkzXKUzYIMtsG9gffIoeTD0IiRYNMaNMyfCmF0HvxwoZ6qJ6U8a9/LoHcUwbXqvm89qxF270Jj5UTIoAu5096HMXtIgld2zFfXnfz+hK87SUYRfyMa5k9FqK0a2ZPfhLnbaaonadrhe+Ss/wvvkUNJh6ERA8HmDe0XGOrNyK1YDIOzSPWktCB0RjiG3oHggaUI7H3/G78edu9Gw3sTEfE3Hnztf7iClR3z1XUn6+HbMUf1nLiKBFrQOL8Coeb1yJn0OszdJ6melBS+ukfOk6qnEJ0QhkaU2j9caBIAcfC1/76qJ6UVa8mPoLMVwrXqvqN+PuzZh4bKyYh465AztRKmvHJFCzuu/bqTfh2+7iQZRYJtaFwwHcGmVcieOBvmHhWqJyWNQ/fIca97FJGgW/Ucog6LKjSEEM8LIfYLIdbGalAyCbVtbb/AMBJqf+0/s7/qSWlH6M1wDP01AnUfw7/vQwBA2FePxvlTEPHsQs7Ud2HqMlbxyo4ROj0cZXci1LQK/tq3VM+JuUjQjab3z0LwwDJkT3gFlp5nq56UdBzDZiHiPwDPZt4jh5JHtGc0/gFgegx2JJ2Qawca35sEGfYip+J9GLMHq56UtmylP4XO0hWuVfcj4m9EY+VUhNpqkD35LZi6nqJ63gmxFv/guNedJCMZ8qJp0bkI7P8UWaf9G5beF6ielJRMXcbD1H0S75FDSSWq0JBSfgSgMUZbjn+8SEgTf4VdO9FYOQmRYEv7hwvlDEvUU0DHIAxW2IfcisDe93HgrTEItWxEzqR5MHefqHraCfvqupNlCOyZr3pOTMiQD42Lzkdg72JknfISrH0uVj0pqTnK7kbEuw+eLf+negpRh4ho/9QkhCgC8JaU8riX85eXl8uqqqpOH6vuv90Q8dV1+vfHkjA6kTPtfZjyR6ueQmg/LV//WhEiwRZkT3odlsIzVU/qNBn2Y//sEugdvZF7xscQQqie1GkyHEDT4gvhr30LmSc/D1vpVaonJT0pJRrePQ1h13Z0ubAaQm9WukcIsVxKqf2LoEgZQ7wPIISYCWAmAPTqFd3HPduH3gYZ0sZFUJae58KYU6Z6Bh2kM9qRPeVtABKm/DGq50Sl/bqT29C69DoE9n0Ic/fTVU/qFBkJounDS+CvfQsZ455hZMRI+z1yZqFxwXR4a16Crd/PVE8i+k5JdUaDKF3IkBf7ZxfDkDUIuRULVc85YTISQvNHl8G3/X/IGPNX2Adeq3pSSpFSouHtMYj4DiB/xiYInVHZFp7RoOPh21uJNEgYrLAPvgWBvYsQ2P+Z6jknREbCaP7kSvi2/w/OUY8wMuLgq3vkbIN3679VzyH6TtG+vfU/AJYA6C+EqBVC/CQ2s4jI1v8X0JnzvvEZIVomZQQtn/0Mvq0vw3nSA3AMvkn1pJRlLjy7/R45q9PvHjmUXKJ918mlUsruUkqjlLJQSsnLoIliRGe0wz74Jvh3v4fAgWWq5xyXlBKtS34Jb/ULcAy/F46yO1RPSmlf3SNnM3zb/6d6DtG34ksnRBpmG3ANhCkbrlW/Vz3lO0kp0br0eng2Pwt72Z1wDPuN6klpof0eOYPS5h45lJwYGkQapjNlwD7oBvh3zUOwcbXqOcckpURb1S3wbHwC9sE3wzni/qR+S24y+eoeOevg2/m66jlEx8TQINI4+8DrIYxOuFbfr3rKN0gp0bbiLrjXPQrbwOvgLH+YkZFglqLvQ59RmtL3yKHkxtAg0jidORu2AdfCt/01BJs3qJ5zFNeq38G95kHY+v0cGaMfY2QoIHR6OIbeiVDjSvhr31Y9h+gbGBpEScA++EYIgxWu1dq5VsO1+gG4vrwX1pKrkDHuKUaGQta+l6XkPXIoNTA0iJKA3pIPW/9fwrftPwi1VqueA9faR9C24i5Yiy9H5vjnIAT/r0QloTPCPvR2BA98gcCeBarnEB2F/+9AlCTsg28BhBGuNQ8q3eHe8Fe0Vd0CS9HFyDzlBQidXukeamcruRI6W6Emr+Wh9MbQIEoSels32Pr9DN7qlxBq265kg3vT39C69HqYe12ArNP+BaGL++2SqIPa75HzawTqPoZ/34eq5xAdxtAgSiKOob8GhIB77R8SfmzPlufRuuQXMBeejewJ/1V6fw06NlvpT6GzdE2qT5Ol1MfQIEoientP2EqugmfL8wi7dyfsuJ6af6Hl05/CXFCB7NNfhdCbEnZs6jhhsMI+5FYE9i5EYP8S1XOIADA0iJKOfejtgAzDtfbhhBzPu+0VtHzyI5i6T0T2pLkQBktCjkudY+v/cwhzLs9qkGYwNIiSjMHZB9a+V8Cz+W8Ie+vieizfjrlo/ugymLqcjOxJb0AYrHE9HkVPZ3TAMfgm+He/i+CB5arnEDE0iJKRY+gdQCQA97pH4nYM36430fTh92HMG43sKW9DZ7TH7VgUW7YB10KYstDGd6CQBjA0iJKQIbMfLEXfh2fjU4j4DsT88X2176Fp8UUw5gxHztR3oTM6Y34Mih+dKQP2gTfAv/N1zd4jh9IHQ4MoSTnK7oIMueFe/5eYPq5/z0I0LW6/K2jO1EroTJkxfXxKDPug6yEMDk19miylJ4YGUZIyZg+GpfeFcG/4KyL+5pg8pn/fR2hceA4MGaXInbYAOnN2TB6XEk9nzoFt4LXwbX9Vc/fIofTC0CBKYo6yWZDBVrg3/jXqxwrs/wxN758Jg6MIOdPeh86SF4OFpJJ98E0QBivcij9NltIbQ4MoiRlzh8NceDbc6/+CSLCt048TqP8CjQvOgM5WgJyKhdBbu8RwJanSfo+cX8C79d8ItdaonkNpiqFBlOQcw+6G9DfCs/GpTv3+YMMKNC6ogM6Sh9yKRdDbusd4IanUfo8cg/J75FD6YmgQJTlT/miYCqbBve4RRILuE/q9wcbVaJg/FTpjBnIqFkFvL4zTSlJFb+sOW7+fwlv9IkKuHarnUBpiaBClAOewuxHx1cOz+dkO/55g83o0zp8CobciZ/piGBy947iQVHIMva39HjlrEn+PHCKGBlEKMHU9BaZup8O99mHIkO+4Xx9q2YzGysmA0CO3YhEMzuIErCRV9PaesJZcCc+W/0PYs0f1HEozDA2iFOEom4WIdy881c9/59eFWmvQUDkJkGHkVCyEIbNfghaSSo6D98hxJ+geOUSHMDSIUoSp+yQY88fBteYhyHDgmF8Tcu1AY+UkyLAPORULYcwalOCVpIrBWQxr8eVwb/obwt79qudQGmFoEKUIIQQcw+5GxL0L3pqXvvHrYXctGt+biEiwFbnTFsCYPVTBSlLJUXYHEPbF9R45RF/H0CBKIeYe02HMLYdrzYOQkdDhnw979qChchIi/gbkTJsPY+4IhStJFUNmf1j6fB+ejU8i4mtQPYfSRFShIYSYLoTYJISoFkLcHqtRRNQ57Wc1ZiHcthXebf8BAIS9dWionIyIZy9ypr4LU94oxStJpXjdI4fo23Q6NIQQegBPAjgDwCAAlwoh+IIvkWLmnufAkD0UrlW/R9hbh8bKKYi4dyJ7ytswdRmveh4pZsweAkvvGXBveDxm98gh+i7RnNEYDaBaSrlVShkA8F8A58VmFhF1lhA6OMpmIdy6CQfmlSHUVo3sSW/A3O001dNIIxxldx28R84TqqdQGogmNHoA2HXEj2sP/hwRKWbpfSEMmQMRCTQjZ9LrMBdMVj2JNMSYexLMhWfBvf7PUd0jh6gjDFH8XnGMn5Pf+CIhZgKYCQC9evWK4nBE1FFCp0f2lHcgQ218dwkdk3PEfQi1VUPobaqnUIqLJjRqAfQ84seFAL7xkXNSymcBPAsA5eXl3wgRIooPg7NI9QTSMGPuCL77iBIimpdOlgEoFUL0EUKYAFwC4I3YzCIiIqJU0OkzGlLKkBDiWgCVAPQAnpdSrovZMiIiIkp60bx0AinlOwDeidEWIiIiSjFCysRdNiGEqAewI8qHyQNwIAZzUgWfj6Px+Tgan4+j8fk4Wiyej95SyvxYjKHUlNDQiAUhRJWUslz1Dq3g83E0Ph9H4/NxND4fR+PzQYnAe50QERFR3DA0iIiIKG6SMTSeVT1AY/h8HI3Px9H4fByNz8fR+HxQ3CXdNRpERESUPJLxjAYRERElCYYGERERxQ1Dg4iIiOKGoUFERERxw9AgIiKiuInqXicnShisUpiciTykphkcGaonAAB0Om30ptGojR1aYTLqVU84LByOqJ4AAAiGtLFDK4wG9f+b8TTsRcDVLDr69Xl5ebKoqCiOiyhRli9ffqAjHz+f2NAwOWHuf3EiD6lpeSdPVT0BAGCzm1VPAAB0784IPVLPPIfqCYc1ewKqJwBtP3FSAAAgAElEQVQA9tS5VE/QlIKu6v8d+fj3Pzyhry8qKkJVVVWc1lAiCSE6dO8y9TlMREREKYuhQURERHHD0CAiIqK4YWgQERFR3DA0iIiIKG4YGkRERBQ3DA0iIiKKG4YGERERxQ1Dg4iIiOKGoUFERERxw9AgIiKiuDluaAghnhdC7BdCrD3i53KEEAuEEFsO/md2fGcSERFRMurIGY1/AJj+tZ+7HcBCKWUpgIUHf0xERBo3ofYlFDcffVOz4uYqTKh9SdEiSnXHDQ0p5UcAGr/20+cBePHg378I4PwY7yIiojjY5RiEyzbdeTg2ipurcNmmO7HLMUjxMkpVnb1NfFcp5V4AkFLuFUJ0ieEmIiKKk61Z5Xi5/wO4fNPtqM4chb4ty/Fy/wewNatcyZ7VBwLIMuvQy9nZb0ekdXH/JyuEmAlgJgDA6Ij34YiI6DjqbUUABIY1LMSHBZfFPTKO/D7Qq1evwz8vpcRL69uwqy2MUV3NmFFqQ3GmMa5bKPE6+66TOiFEdwA4+J/7v+0LpZTPSinLpZTlwmDt5OGIiCgWnIEDuGbVVbCFWlDV5SyM3P/2N67ZiLUjvw/k5+cf/nkhBO4Zm40LS2xY1xDAHZ804Q/LmrG5KRjXPZRYnQ2NNwD86ODf/wjAvNjMISKieDkUGVmBOrzR5ya8WnoPXu7/wFHXbCR8k0mHi/s78OSkXHy/nx2bm4K4+7Mm3L+0CesbAko2UWx15O2t/wGwBEB/IUStEOInAB4CMFUIsQXA1IM/JiIijXIGDmDm2qvhDDbgjT434bOCSwB8dc1GT9d6pftsRh1mlNrx5KRcXDbAjp2tIfz282bcu6QJqw8EIKVUuo8677jXaEgpL/2WX5oc4y1ERBQHhyIj01+H54Y8he0Zw4/69a1Z5couBv06i0GHc/vaUVFkw8KdXrxR48HvlzajJMuAC0vsGNHFBCGE6pl0AniZLxFRCjsyMp4f/Ng3IkOrzHqBM/vYMLWXFR/U+jCvxo0/VLWgT4YBM0rtKO9qgo7BkRQYGkREKSpZI+NIRr3A1N5WTOxpwce7fZhb7cEjy1vQ06nHBSV2jOtuZnBoHEODiCgFpUJkHMmgE5jY04rTeljw2V4/5la78fjKVry2WY/zS2w4pcACvY7BoUUMDSKiFJNqkXEkvU7g1B4WnFxgxtJ9fszZ4sFTq9owe4sb5/W1Y0KhBQYGh6YwNIiIUkgqR8aRdEJgXHcLxnQzY8X+AGZvcePZNYeCw4aJPa0w6RkcWsDQICJKEekSGUfSCYHyrmaM7GLCqvoAZld78Pw6F+ZWe3BOsQ1TelthZnAoxdAgIkoB6RgZRxJCYHgXM4blm7CuIYg51W68tMGF12vcOLvYhmm9rbAaOvsZlRQNhgYRUZJL98g4khACQ/JMGJJnwsbGAOZUe/DvjW7Mq/HgzD42nFFkhd3I4EgkhgYRURJjZHy7ATkm3DnahOrmIOZWu/HqZjfe2urB9CIrzuxjQ4aJwZEIDA0ioiTFyOiYkiwjbi3PwvaWIOZUe/B6tQfvbPNiWm8rzu5jRZZFr3piSmNoEBElocxQA2auvZWRcQKKMo24aWQmattCmFPdfnbjve0eTOllxbl9bchhcMQFQ4OIKMlkhhpwx+6bkRmuZ2R0QqHTgOtHZOJ7/UJ4vdqD+Tu8WLDTi9MLLTivrx1dbAyOWEpsaOj0gD0roYfUsi5dM1RPAAAUdHWongAA2FPnUj0BADCoT47qCQCAXQe08XwAQFubNm7XrZV/V1Vqf7nkVmSG6/HH7g9hs7cE8Kr7dyUYiig7drS62w345bAMXFhqx7waNz6o9WHxLh9O7WHBBSU2dLPzz+KxwGeRiChJfP2ajM3eEtWTUkIXmx4/G5qBGSV2vLnVg/d3evFhrQ8nF5hxQYkdhU5+q4wGnz0ioiRwzAs/FZ7JSEW5Vj2uHOzE+X1teGubF/N3ePHpHj9GdzNjRqkNRRlG1ROTEkODiEjj+O6SxMqy6HH5QAfO7WvDO1s9eG+HF0v3+VHe1YQZJXb0zWJwnAiGBhGRhjEy1Mkw6XDJAAfO7mvDe9u9eGebB3d+2oTh+e3B0T+HwdERDA0iIo1iZGiDw6jDRaV2nFlkxfwdXry1zYPfLGnC4FwjZpTYMTjXCCF4P5Vvw9AgItIgRob22Iw6nF9ix/QiGxbu9OLNrR7ct7QZ/bONmFFqw7A8E4PjGBgaREQaw8jQNotB4KxiG6b2tmLxLi/m1Xjw4Bct6JtpwIxSO0Z2YXAciaFBRKQhjIzkYdILVBTZMLmXFR/W+jCvxo2Hq1rQO8OAGSU2jO5mho7BwdAgItIKRkZyMugEJvey4vRCCz7Z48Pcag/+vKIVPRx6XFBix/juZuh16RscDA0iIg1gZCQ/vU5gQqEVp/awYMleP+ZsceOJL1vx2mY9zi+x4dQeFhjSMDgYGkREijEyUotOCJxcYMG47mZU1QUwe4sbz6xuw+wtbpzX147TCy0w6tMnOBgaREQKMTJSl04IjO5mxqiuJqzcH8Dsajf+vrYNc6rdOLe4/doOUxoEB0ODiEgRRkZ6EELgpK5mjOhiwpqGIOZsceMf612YW+PBOX1smNrbAotBp3pm3DA0iIgSYELtS9jlGIStWeUA2iPjmlVXwRlswHNDnmJkpAEhBMryTCjLM2F9QwBzqt3410YX5tW4cVaxDRW9rbAZUy84Uu+/ERGRBu1yDMJlm+5EcXPV4cjICtTh7aLrGRlpaFCuCbPGZOO+8dkozTbiv5vcuHZRA/632QVXIKJ6XkzxjAYRUQJszSrHy/0fwOWbbgcgYAu14I0+N+GzgktUTyOF+mUbcduoLGxtCWLOFg9mb/HgnW1eTOttxfl9bSlxhoOhQUSUIPW2IgAC9lALqrqcxcigw4ozjbilPBM7W0OYW+PGvBoPdAAuGeBQPS1qyZ9KRERJ4NDLJbaDkTGg8VMUN1epnkUa0yvDgBtGZMKiFwhEpOo5McHQICKKsyOvyXijz014tfQevNz/gcPXbBClMoYGEVEcHXoLqzPYcNQ1GYeu2ejpWq94IVF88RoNIqI4OfJzMo71FtatWeWH3+5KlKp4RoOIKA74YVxE7RgaREQxxsgg+gpDg4gohhgZREdjaBARxQgjg+ibGBpERDHAyCA6NoYGEVGUGBlE346hQUQUBUYG0XdjaBARdRIjg+j4GBpERJ3AyCDqGIYGEdEJYmQQdRxDg4joBDAyiE4MQ4OIqIMYGUQnLqE3VdNZrLCXDEnkITWtoKtD9QQAQJbNpHoCACCrT47qCQCAPvlW1RO0J0/1gHYZNqOyY9t89Tjv02vhDNThzfF/QyR3JHopW6MdNQb+eZW+G/8NISI6jvbI+DGc3r14c9zfsDd3pOpJREmDoUFE9B0YGUTRYWgQEX0LRgZR9BgaRETHwMggig2GBhHR1zAyiGKHoUFEdARGBlFsMTSIiA5iZBDFHkODiAiMDKJ4YWgQUdpjZBDFD0ODiNIaI4MovhgaRJS2GBlE8cfQIKK0xMggSgyGBhGlHUYGUeJEFRpCiBuFEOuEEGuFEP8RQlhiNYyIKBZGbPk/9KhfevjHNl89LvzwUmS6dzIyiBKg06EhhOgB4HoA5VLKIQD0AC6J1TAioljYnzUEFVU3o0f90sOR4fTtxaeDb2FkECWAIQa/3yqECAKwAdgT/SQiotjZnT8GleWPYPqyGyEhYAk24+Mht2NN3ytUTyNKC50+oyGl3A3gTwB2AtgLoEVKOf/rXyeEmCmEqBJCVElfW+eXEhF1UpOzGBIC1mAzNvY8j5GRYEd+H6ivr1c9R/OW7fPDH5awGlLjMspoXjrJBnAegD4ACgDYhRCXf/3rpJTPSinLpZTlwuLs/FIiok6w+epx0YeXwBJsxoae56Go7qOjrtmg+Dvy+0B+fr7qOZq2ZK8Pf17RgpIsA87sY1U9JyaiyaUpALZJKeullEEAcwCMj80sIqLoHYoMh28fPh5yOxad9AAqyx85fM0GkZZ8stuHx1a0ojTLiDtHZ8FuTPMzGmh/yWSsEMImhBAAJgPYEJtZRETRsfnqcf6nV8HmP3DUNRmHrtno0rxW8UKir3ywy4snvmzFoFwj7hidCVuKRAYQxcWgUsqlQojXAKwAEAKwEsCzsRpGRNRZhyLD4d2HeSc//413l+zOH4Pd+WMUrSM62vs7vXhuTRvK8oy4pTwLZr1QPSmmonrXiZTyHgD3xGgLEVHUjowMfk4GaV3ldg+eX+fCiHwTbhqZCVOKRQYQ/dtbiYg0g5FByeTtrR68tMGF8q4m/GpEJowpGBkAQ4OIUgQjg5LJ69Vu/GeTG2O7m3Hd8AwYdKkZGQBDg4hSACODkslrW9x4dbMbJxeYcc2wDOhTODIAhgYRJTlGBiULKSVe2ezG3GoPJhRa8IsyJ3QitSMDYGgQURJjZFCykFLi5Y1uvLnVg8k9Lfjp0PSIDIChQURJipFByUJKiRfXu/Dudi8qeltx5WBH2kQGwNAgoiTEyKBkEZESz691YcFOL87qY8UVAx0QaRQZAEODiJIMI4OSRURKPLu6DYtrfTivrw2X9renXWQADA0iSiKMDEoW4YjE06vb8PFuHy4steF7pekZGQBDg4iShN1fj/NXzGRkkOaFIhJPrmrFZ3v8+H4/O2aU2lVPUoqhQUSaZ/fX47Kqn8Hhr2NkkKaFIhKPrWzFF/v8uHyAA+f0tamepFxCQ8NsNqJ3366JPKSm9cpL78r9un75VtUTAADZVvb31+XajMqObfLsR/nbP4c5sB8rz/gHnN1GwalsDX2dyahXPUEzgmGJP69owfL9AVw5yIEz+jAyAJ7RICINa4+MH8Ds3oeV019Ac7dRqicRHVMgLPHI8hZ8WR/AT4c4MbW3Nv7gpAUMDSLSJEYGJQtfSOLhqmasawjiF2VOTOzJyDgSQ4OINIeRQcnCG4rgD8tasLExiKuHOXFaISPj6xgaRKQpjAxKFp5gBA8ua0F1cxDXj8jA+AKL6kmaxNAgIs1gZFCycAUjePCLZmxrCeFXIzIwpjsj49swNIhIExgZlCzaAhHcv7QZta4Qbh6ZiZFdzaonaRpDg4iUY2RQsmjxR3D/0ibsdYdx68hMDO/CyDgehgYRKcXIoGTR5Avj/qXN2O8J47ZRWRiaZ1I9KSkwNIhIGUYGJYtGXxi/+7wZjb4I7hidhUG5jIyOYmgQkRKMDEoW9Z4w7lvahNaAxJ2jMzEgh5FxIhgaRJRwjAxKFnWeMO77vAmeoMSsMVkoyVL3cfzJiqFBRAnFyKBksdcdwn2fN8Mflpg1NgvFmYyMzmBoEFHCMDIoWdS2hXD/0maEpcRvxmajdwa/XXYWnzkiSghGBiWLnW0h3P95EyAE7hmbjUInv1VGQ6d6ABGlpt6r/obsPUsAfBUZlrbd2FtyPiODNGt7SxC/W9IEnRC4d2wWIyMGGBpEFBet+WUoW3Qdumx793BkRAwm1BWfpXoa0THVNAfxu6XNMBsE7h2XhQIHIyMW+CwSUVw0FYzDhpPvQ9mi6yChQ9hoxaopz6CpYJzqaUTfsLkpiAe+aIbTKHD32Gx0selVT0oZDA0iiguTZz9Kqh6BhA46GcL2wVcyMkiTNjQG8NAXLci26DBrTBbyrIyMWOJLJ0QUc0dekxE2WrF1xHUo3PDy4Ws2iLRi7YEAHvyiGTkWHX4zlpERDzyjQUQxdWRkRAymwy+XNHYfi7JF12H1pL/yzAZpwqp6Px6uakE3ux6zxmQjy8w/e8cDn1Uiipkj38K6p9+FR12T0VQwDqsn/RUZ9asVryQCVtT58ceqFvRwGPCbsYyMeOIZDSKKiY58TkZTwTiezSDllu3z488rWtA7w4C7RmfBYWJkxBNDg4iixg/jomSxZK8Pf13ZiuJMA+4cnQWbkZERbwwNIooKI4OSxce7fXjyy1b0zzbi9tGZsBoYGYnA0CCiTmNkULL4YJcXz6xuw6BcI35dngWLQaielDYYGkTUKYwMShbv7/TiuTVtKMsz4pbyLJj1jIxEYmgQ0QljZFCyqNzuwfPrXDipiwk3npQJEyMj4RgaRHRCGBmULN7a6sE/N7gwqqsJvzopEwYdI0MFhgYRdRgjg5LF69Vu/GeTG2O7m3Hd8AxGhkIMDSLqEEYGJQMpJWZXe/DqZjdOLjDjmmEZ0DMylGJoENFxMTIoGUgp8cpmN+ZWezCh0IJflDmhE4wM1RgaRPSdLN79KF/0I0YGaZqUEi9vdOPNrR5M7mnBT4cyMrSCoUFE38ri3Y9JC38Is6eOkUGaJaXEi+tdeHe7FxW9rbhysIORoSEJDQ2dTsDpNCXykNQBuXZt9Ga2VRs7qhu8qicAAEpyrUqPb3TXod87V8LkrcMKjUTG5npt/LPZecCtegIAoFeeXfUE5SJS4vm1LizY6cVZfay4YqADgpGhKdr4f3Yi0hSjuw795lwCk3svtpz7EpqzT1I9iegbIlLi2dVtWFzrw3l9bbi0v52RoUEMDSI6ytcjw1UwGvCGVM8iOko4IvH06lZ8vNuPi0ptuKiUkaFVDA0iOuyYkUGkMaGIxBNftmLJXj++38+OGaV8CUnLGBpEBICRQckhFJF4bGUrvtjnx+UDHDinr031JDoOhgYRMTIoKQTDEn9e0YLl+wO4cpADZ/RhZCQDhgZRmmNkUDIIhCX+tLwFq+oD+OkQJ6b2VvuuLOo4hgZRGmNkUDLwhSQermrGuoYgflHmxMSejIxkwtAgSlOMDEoG3lAEf1jWgo2NQVw9zInTChkZyYahQZSGGBmUDDzBCB78ohnVLSFcNyIDJxdYVE+iTmBoEKUZRgYlA1cwggeWNmN7awi/GpGBMd0ZGcmKoUGURhgZlAzaAhHcv7QZta4Qbh6ZiZFdzaonURQYGkRpgpFByaDFH8H9S5uw1x3GrSMzMbwLIyPZMTSI0gAjg5JBky+M+5Y2o94Txm2jsjA0jzfhTAW6aH6zECJLCPGaEGKjEGKDEGJcrIYRUWwwMigZNHjD+O3nzTjgjeCO0YyMVBJVaAB4DMB7UsoBAIYB2BD9JCKKRrflT8NZ+xmAIyKjrRYN/WcwMkiT6j1h/PbzJjT7I7hzdCYG5TIyUkmnQ0MIkQHgNAD/BwBSyoCUsjlWw4ioc9xdh6H43auRVf324ciQejOaSs9WPY3oG+oORkZbQGLWmCwMyGFkpJpozmgUA6gH8IIQYqUQ4u9CCN5Cj0ixtsLx2DHx9+j73jUwt+6A1JtRc9azaCscr3oa0Te8ssmFem8EQ/JMyLfqVc+hOIgmNAwATgLwtJRyBAA3gNu//kVCiJlCiCohRFXQzRMeRPFmdNehx+ePQAoddJEQ9g+7ipFBSh35faC+vv6oX7u4nx3jC8xYts+PaxcdwIvr29DoCytaSvEQTWjUAqiVUi49+OPX0B4eR5FSPiulLJdSlhvtWVEcjoiOx+iuQ7+5l8LUugsRgw17Rt2A/DX/PHzNBpEKR34fyM/PP+rXutkNuGFEJh6dkINxBRa8t92L6xY34O9r2lDvYXCkgk6/vVVKuU8IsUsI0V9KuQnAZADrYzeNiE7EkZFx5MslbYXjUPzu1dh6xlM8s0GaVeAw4OphGbiw1I551W4s2uXFol1enFZowfl9behm56cxJKto/8ldB+BlIYQJwFYAV0U/iYhO1OHIcO1Bw8Dvoan07MNR0VY4HlvPeAr2ulUMDdK8rjY9ZpZlYEapHW/UeLBolxcf1vpwcoEFF5TY0MPB4Eg2Uf0Tk1J+CaA8RluIqBOOjIxv+5yM9jMbjAxKHnlWPX48xIkLSmx4c6sH7+/04pPdPoztbsaMEjt6ZTA4kgX/SRElsY5EBlEyy7bo8cNBTpzf1463tnlQud2LJXv9GNXVhBmldhRnGlVPpONgaBAlKUYGpZMMsw4/GODAucU2vLPdg3e3ebGsrgnD8024sNSOftkMDq1iaBAlIUYGpSuHSYeL+zlwdh8bKnd48fZWD+7+rAlDco24sNTOTxXVIIYGUZJhZBABNqMOF5TYcUaRFQt2+PDmNg9++3kzBuQYMaPEhrI8E4QQqmcSGBpESYWRQXQ0i0GHc/raUFFkxcKdXryx1YMHvmhBSZYBM0rsOKkLg0M1hgZRkmBkEH07k17gjD42TOllxQe1PsyrceOPVS0oyjBgRokNo7qZoWNwKMHQIEoCjAyijjHqBab2tmJiTws+2e3D6zUePLqiFYUOPWaU2DGugMGRaAwNIo1jZBCdOINO4PSeVpzaw4Ile/2YW+3G41+24tUtelxQYsPJBRYYdAyORGBoEGkYI4MoOnqdwCk9LBhfYMYX+/yYU+3BU6va8NpmN84vsWNCIYMj3hgaRBrFyCCKHZ0QGNvdgjHdzFixP4DZW9x4dk0bZm9x49y+NkzqaYVJz+CIB4YGkQZZvPvR792rGBlEMSaEwMiuZpzUxYRVBwKYs8WDF9a5MLfag7OLbZjaywqLgcERSwwNIo2xePdj4sIfweTdx8ggihMhBIbnmzEsz4R1DUHMqXbjXxtcmFfjxll9bKjobYXNqFM9MyUkNDT0OgGnjZ/adki/fKvqCQCAokxt7Nje4lU9AQBQkqvu+TC669rPZHj3YfGE53DAPBRoUP+8NLhDqidoSq88u+oJAICdB9yqJyAQDKueEBUhBIbkmTAkz4SNjQHMrfbgv5vceHOrB2cU2XBGHyscDI6o8IwGkUZ8/ZqMA+ahqicRpZUBOSbcMdqE6uYg5la78doWN97e5kFFbyvOKrYhw8Tg6AyGBpEGHPPCTw2cySBKRyVZRtxanoXtrUHMrfZgXo0H7273YGovK84ptiHLolc9MakwNIgU47tLiLSpKMOIG0/KRG1bCHOr3Xh7mxeVO7yY3MuKc4ttyLUyODqCoUGkECODSPsKnQZcNyITF/ULYV61Bwt2eLFghxen97Tg/L52dLExOL4LQ4NIEUYGUXLpbjfgF8MyMKPUjjdqPFhc68XiXT6c2sOCC0ps6G7nt9Rj4bNCpAAjgyh5dbHp8dOhTlxQYsObWz14f6cXH9X6ML7AjAtK7Ojp5LfWI/HZIEowRgZRasi16nHlYCfO72vDW9u8mL/Di0/3+DGmmxkzSmwoyjSqnqgJDA2iBGJkEKWeLIselw904Ny+NryzzYP3tnuxdJ8fI7uYMKPUjpKs9A4OhgZRgjAyiFJbhkmHS/o7cE6xDe9u9+KdbR7c9WkTyvJMuLDUhgE56fmBlQwNogRgZBClD7tRh4tK7TizyIoFO714a6sH9yxpxqAcIy4stWNwrhFCpM/9VBgaRHHGyCBKTzajDuf1taOitw0Ld3rx5lYP7lvajH7ZRswosWF4viktgoOhQRRHjAwishgEziq2YWpvKxbv8mJejQcPLWtBcaYBM0rsKO+a2sHB0CCKE0YGER3JpBeoKLJhci8rPqz1YV6NG39a3oJeTgNmlNowppsZuhQMDt4hhihGui1/Gs7azwAcERmtu9DQfwYjg4gOM+gEJvey4s8TcnH1MCdCEYm/rGjFLR814uNaH8IRqXpiTDE0iGLE3XUYit+9GlnVbx+ODKk3o6n0bNXTiEiD9DqBCYVWPDIhBzeMyIBOAE+sasVNHzZi8S4vQikSHHzphChG2grHY8fE36Pve9dACh0iBhtqznoWbYXjVU8jIg3TCYHxBRaM7W5GVV0A/9rgwjOr2xCMSEzrbVM9L2o8o0EUI0Z3HXos+ROk0EEXCWH/sKsYGUTUYTohYDcKNPsj6G7Xo7yrWfWkmGBoEMWA0V2HfnMugamtFhGDDXtG3YD8Nf88fM0GEdHxrD4QwENfNCPfqsO9Y7OQY0mNu8LypROiKB0ZGVJvPvxySVvhOBS/ezW2nvEUz2wQ0Xdaud+PR5a3oLvdgLvHZCHDnDrnAVLnvwmRAocjw70XDQO/d9Q1GW2F47H1jKdgr1uleCURaVlVnR9/Wt6CQocBvxmbWpEB8IwGUacdGRnf9jkZ7Wc2eDaDiI7t870+PL6yFX0yDbhzdBbsxtSKDIChQdQpHYkMIqLv8uluH55Y1YqSLCPuGJUJWwpGBsDQIDphjAwiitaHtV48vaoNA3OMuG1UJiyG1IwMgKFBdEIYGUQUrYU7vXhuTRuG5Blxa3kWzPrU+9jxIzE0iDqIkUFE0arc7sHz61wYnm/CzSMzYUrxyAAYGkQdwsggomi9vc2Dl9a7MLKLCTeelAljGkQGwNAgOi5GBhFFa16NG//e6MaYbmZcPyIDBl16RAbA0CD6TowMIorW7C1u/G+zG+MLzLh2WAb0aRQZAEOD6FsxMogoGlJK/G+zG3OqPTi1hwVXD3NCJ9IrMgCGBtExMTKIKBpSSvx7oxtvbPVgYk8LZg5Nz8gAGBpE38DIIKJoSCnx0gYX3tnmxdReVvx4iCNtIwNgaBAdxeTZj37vXMbIIKJOiUiJF9a5MH+HF2cUWfGjQQ6INI4MgKFBdJjJsx/lb/8AJs8+RgYRnbCIlHhuTRsW7fLhnGIbLhtgT/vIABIcGjazHiN6ZybykJpWlGlVPQEAkGk2qp4AABjWRd0OnWsfsuZcDp2nDi2X/A85PcciR9ka7cm1hVRPII0yGfWqJ2hCREo8vaoNH+32YUaJDRf3Y2QcwjMalPZ0rn3I+tf50LXtRcslryDYc6zqSUSURMIRiSdXteLTPX5c3M+OC0vtqidpCkOD0hojg4iiEYpIPL6yFUv3+fGDAXac15eR8XUMDUpbjAwiikYwLPGXlS2oqgvghwMdOKvYpnqSJjE0KC0xMogoGoGwxCPLW/BlfQA/HuxARREj49swNCjtMDKIKBr+sMTDVc1YeyCImUOdmNxLGxf2axVDg9IKI4OIouELRfCHZS3Y0BjEL4c5MaGQkXE8DA1KG4wMIoqGJxjBQ8tasKU5iGuHZ+CUHi3RBHsAABwISURBVBbVk5ICQ4PSAiODiKLhDkbwwBfN2NYSwg0jMjC2OyOjoxgalPIYGUQUDVcggvuXNmNnWwg3npSJUd3MqiclFYYGpTRGBhFFo9UfwX1Lm7HXHcKt5ZkY0YWRcaIYGpSyGBlEFI1mXxj3LW1GnSeMX5dnoSzfpHpSUmJoUEpiZBBRNBp9Yfzu82Y0+sK4fVQWhuQxMjpLF+0DCCH0QoiVQoi3YjGIKFqMDCKKxgFvGPcuaUazP4I7RzMyohV1aAC4AcCGGDwOUafYljwO4/aPARwRGS218A39HiODiE7Ifk8Y9y5pQlugPTIG5DAyohVVaAghCgGcBfx/e3ceH1V973/89UlmycxkYxMFBALEKPpDxYC4FtS27ojyq3i1vfoo1Xr1dvlZ+9O21qW1t9f663Kt1+qPVnurV39eFPe1FmsXpSDiCpSwiqAsknUmy2S+vz8ScxMWjcnMnDkz7+fjkQfJZJjz4YSc88o5ZyYsSM84Ip9exwFHUrFoPuGVj/VEBoEQbQfP9no0EfGR91uS3PDyLuJJx3UzKjloSNDrkfLCYK/R+DnwbaAsDbOIDEjH+BNoOvUnlD8yH6wYF4zQcN5v6Rh/gtejiYiP/N83m9jZmuLz4yKMKdUljOky4CMaZnYmsM059+on3O9SM1tmZsviDbsGujiRfSpqfp/YH38EVoylOkjUfkWRIZJDeu8Htm/f7vU4+3R+TSk1Q4I8uzHB1xbv5Ml1cdo6nddj+d5gTp0cB5xtZhuAB4CTzOze3e/knLvLOVfrnKuNVgwZxOJE9tT7mgwXjNBy3FVElt/Tc82GiHiv935gxIgRXo+zTwcNCXLjMZVcd3Qlo0qL+Y+VzVz5hx08uraFRDLl9Xi+NeBjQ865a4FrAcxsJvAt59xFaZpL5BP1jgwCoZ7TJe3jjqdi0Xwa5izQkQ0R+VTMjMOGhzhseIhVH7azqC7Of65q4bG1cU6vinLq+AixYDqeR1E4dBJKfKn3U1hbp8yj7ZDZPVHRMf4EGuYsILj1NYWGiAzYwUNDXDs9RF19B4vqWnjw7y08vi7OqeMjnF4VpTyk4OiPtISGc+5F4MV0PJbIJ+nzOhnnP0DH2GP2uE/H+BMUGSKSFpMqg1xdW8mGxg4eXhPnkbo4T61P8LlxEc6cEKUyrOD4ODqiIb7Sn8gQEcmE8eVB/tdRFWxuSrKoroUn1sV5ZkOcU8ZGOHtilKElxV6PmJMUGuIbigwRyQVjygL885EVzD0oySN1cZ7bmOD5TQlmjYkwe2KUEVEFR28KDfEFRYaI5JoDYgEuP7yc86pjPLq2hcWbE/zh3QQnji7hnElR9o9pFwsKDfEBRYaI5LL9osV85X+Uc+6kGI+vi/P7TQle3NzKcaPCzJkUY0xZYe9qC/tfLzlPkSEifjEsUszFh5ZxzsQoT6xP8NzGBH/Z0sbRB4Q5d1KMceWFucstzH+1+IIiQ0T8qLKkmIsOKeXsiVGeWh/nmQ0JXtnaRu3IEOdOijGxsrB+h4pCQ3KSIkNE/K48VMS8mlLOmhDl6Q0Jnlof5zt/2cURI7qCo2ZoYQSHQkNyjiJDRPJJLFjE3OoYp4+P8PymBE+si/P9l3dx6LAg506KceiwIGbm9ZgZo9CQnKLIEJF8FQ0WMXtijM+Pi/LCpgSPr4vzgyX11AwJcm51lMOHh/IyOBQakjMUGSJSCEoCxhkTonx2XITF7yZ4dG2cf/lbAxMrApxXHWPqfvkVHAoNyQmKDBEpNKFi4/Pjo5w8NsIfN7fy6NoWblnWwLjyAOdOijJ9/zBFeRAcCg3xnCJDRApZoMg4eWyEmWNK+POWVhbVxfnZ8kbGlBYzZ1KMY0f5OziyGhqBImNYtDCusvWTWNC73rSm9ym5bw7WtJXN595L68haaOvwbB6ADQ0JT5f/kb9vz405ABrj3n5NZO/e3dHs9Qgk2pJej5A3iouMz4yJcMLoEl7e2sbDa1q4bUUjC9cUc87EKMePLiFQ5L/g0BEN8Yw1vU/Jb8/CGrfQeuF/0TrsSK9HEhHxXJEZx40q4ZgDwiz7oJ2H1rRwxxtNLFzTwuyJMWaOKSFY7J/gUGiIJ3aPjNS4Y6E5d36CFxHxWpEZ0/cPM21kiNe2tfNQXQsL3mri4boWzp7QdW1HyAfBodCQrNtrZIiIyF6ZGVNHhjlyvxBv7uzg4TUt3PNOM4vWxjmrKspnx5VQEijyesx9UmhIVikyREQGxsyYMjzElOEh3tnZzsN1Ldy7qplH17ZwxoQonx8XIRrMveBQaEjWKDJERNJj8rAQk4eF+PuuDhbVtfDA6hYeXxvn1KoIp4+PUhrKneBQaEhWKDJERNLvoCFB/ve0StY1dPDwmjgPrYnz1PoEnxsX4cyqKOVh74NDoSEZp8gQEcmsCRVBvlVbwabGJA/XtfDY2jjPbIhzytgIZ02IMqSk2LPZFBqSUYoMEZHsGVse4BtTK9jSnGRRXZynNyR4bmOCkw6McPbEKMMj2Q8OhYZkjCJDRMQbo0oDXHFEOedVR3l0bZzfb0rw+00JPjOmhHMmxRgZzV5wKDQkIxQZIiLe2z8W4LIp5ZxXHePRtXEWv5vgxc2tHD+qhDmToowqzXwGKDQk7RQZIiK5ZXikmC8fVsacSVGeWBfn+Y0J/vReK8eMCjNnUoyxZZnLAYWGpJUiQ0Qkdw0tKeZLk8uYPTHGk+vjPLshwV+3tDFtZJjzqqNUVaT/95EpNCRtFBkiIv5QES7iHw4u5ewJUZ7aEOfp9QmWftDG1P1CnDspRvWQ9AWH90+wFd8K/uUXFK1/CegVGfXvkpzyBUWGiIgPlIaK+MJBpdx+0jDOPyjGml0dfO+vu7h5yS5WftielmUoNGTAOkcdScnCSyh++5GeyCAQInnoHK9HExGRTyEaLOLc6hi/PGkYFx1cysamTm54uZ4bXt7FGzvacc4N+LEVGjJgqaoTaZ17N+FFX8F2bYBAiNbz7yVVdaLXo4mIyACUBIo4a2KUX84axsWTS3mvOcnNS+p5ddvAj24oNGRQUlUn4iJDsVQHHUdfpsgQEckDoWJjTFmA1qRjdGkxB1UO/JoNhYYMStH6l7D4DlJDJxBc9pueazZERMS/Vmxr41+X1rN/rJjrZwwZ1O9MUWjIgBWtf4mShZeQqhxPalg1rXPvpmThJYoNEREfe/WDNn7yagOjSwNcN2MIFYP8xWwKDRmw4i2v0Tr3bgiXAf99zUbxltc8nkxERAZiydZW/s+rDYwrC3DdjErK0/Dr5vU6GjJgHcd9fY/bUlUn6joNEREf+uuWVm5b0cikigDXTq8kGkzPsQiFhoiISIF7aXOCf3+9iZqhQa6ZVkEkkL4THgoNERGRArb43QR3vtHEocOCXF1bSUnA0vr4Cg0REZEC9fzGBAveauLwESG+dVQFoeL0RgYoNERERArS0+vj3PNOM1P3C/HNqZmJDFBoiIiIFJzH18a5d1Uz00aG+cbUcgJFmYkMUGiIiIgUlEV1LTywuoVjDghz5RGZjQxQaIiIiBQE5xwL17SwcE2cE0aHuXxKOcUZjgxQaIiIiOQ95xwPrG7hkbVxZo4p4bIpZRRZ5iMDFBoiIiJ5zTnH71Y28+T6BKeMLeHLh2UvMkChISIikrdSznHP2808uzHBqeMjXDy5FMtiZIBCQ0REJC+lnGPBW028sKmVM6siXHRI9iMDFBoiIiJ5J+Ucv3qjiT9ubmXOxCjn18Q8iQzIcmiUhYKcXDUim4vMaZFQsdcjAFAyyDmaA0VYsIhhQyNpmshbo0pz498xviLh9Qg9Gto7vB5B9qJup/f/V/8e0c+ruaYz5bj99Ub+sqWNudUx5lZHPYsM0BENERGRvJFMOW5b0cgrW9uYVxNjzqSY1yMpNERERPJBMuX4+fIGln7QzkWHlHLWhKjXIwEKDREREd9r73T8bHkDy7e1c/HkUk6ryo3IAIWGiIiIr7V3On6yrIE3drQz/7AyPjvO+2t3elNoiIiI+FRr0nHLsnre2dnBV6eUMevA3IoMUGiIiIj4UiKZ4sdLG1j9YQdXHF7OCWNKvB5prxQaIiIiPhPvSPEvf6unriHJ144s59hRuRkZoNAQERHxleb2FDf/rZ6NjUm+ObWC6fuHvR7pYyk0REREfKKxPcXNS+rZ3JzkqqMqOGpkbkcGKDRERER8ob4txQ+X7OL9lk6urq3giBG5Hxmg0BAREcl5u1o7uemVena2dnLNtEoOGx7yeqR+U2iIiIjksJ2JTn6wpJ5drSmumVbJ5GH+iQyAooH+RTM70MwWm9lKM3vbzL6ezsFEREQK3bZ4Jze8vIv6thTfme6/yIBBhAaQBK5yzh0CzACuMLPJ6RlL/KDthVtJrnmxz23JNS/S9sKt3gwkIpJH3m9JcuMru2jpcHzv6Epqhga9HmlABhwazrmtzrnl3e83ASuB0ekaTHJf8dha4ndfiGttAroiI373hRSPrfV4MhERf9vSnOTGV+ppTTqum1HJpEp/RgYM7ohGDzMbDxwJLEnH44k/BKpnUjLvdlI71tK5o4743RcSveQ+AtUzvR5NRMS3Njd1RUYy5bh+xhCqKvwbGZCG0DCzUuAh4BvOuca9fP5SM1tmZst2fbhjsIuTHJJq3k7bUz8AM9y2NYSOv1SRISJ76L0f2L59u9fj5LSNjV2nSwCunzGEseX+f87GoELDzIJ0RcZ9zrmH93Yf59xdzrla51ztkKHDB7M4ySGp5u20/PJUUttWQyhG+PPfof3Pd+1xzYaISO/9wIgRI7weJ2eta+jgpld2ESgyrp9RyZgy/0cGDOLprWZmwK+Blc65n6ZvJMl1fSIjGCE2/78IVM8kUP0ZnT4RERmAuvoOfrSknkjQuO7oSvaP5UdkwOCOaBwHfBE4ycxWdL+dnqa5JEf1RMaOOoLTv9gTGdB1zUb0kvvo3LTM2yFFRHxk9Ycd/HBJPbGgccOMIXkVGTCIIxrOuT8DlsZZJMf1jozYVxYRqDlpj/t0HdmYmf3hRER86J2d7fx4aQNDS4q47uhKhkWKvR4p7fIrmyRj+hMZIiLSf2/uaOeWpfUMjxTz/RmVDCnJv8gAhYb0gyJDRCS9Vmxv49ZlDewfK+Z7Rw+hMpyWV5vISQoN+ViKDBGR9Hr1gzZ+uryBMaUBvnt0JeWh/I0MUGjIx1BkiIik15KtrfzitUbGlQf47vRKSvM8MkChIfugyBARSa+/bmnlthWNTKoIcO30SqLB/I8MUGjIXigyRETS60+bW7n99UZqhga5ZloFkUBhRAYoNGQ3igwRkfRa/G6CO99o4tBhQa6uraQkUFivDKHQkB6KDBGR9Hp+Y4IFbzUxZXiIq2srCBUXVmSAQkO6KTJERNLr6fVx7nmnman7hfjm1MKMDFBoCIoMEZF0e3xtnHtXNTNtZJhvTC0nUFSYkQEKjYLnmrfTctcZigwRkTRZVNfCA6tbOOaAMFceUdiRAQqNguaat9Nx5xm4nYoMEZHBcs6xcE0LC9fEOX5UmH86vJziAo8MyHJohAJFjB4ayeYiZR86m7bz4YIzcyoyhpaGvB5B9qGiI+j1CJKjwgX0NM2P45zjgdUtPLI2zmfGlPDVKWUUmSIDdESjIHU2befDX3yO5LY1DP2nR+mcMNPrkUREfMs5x+9WNvPk+gSnjC3hy4cpMnpTihaY3SMjfPDJXo8kIuJbKee4++2uyDh1fIT5iow96IhGAVFkiIikT8o5FrzVxAubWjmzKsJFh5Riiow9KDQKhCJDRCR9Us7xqzea+OPmVs6ZGGVeTUyRsQ8KjQKgyBARSZ/OlOP21xv5y5Y25lbHmFsdVWR8DIVGnlNkiIikTzLluG1FI69sbWNeTYw5k2Jej5TzFBp5TJEhIpI+yZTj58sbWPpBOxcdUspZE6Jej+QLCo08pcgQEUmf9k7Hz5Y3sHxbOxdPLuW0KkVGfyk08pAiQ0Qkfdo6Hbcua+CNHe3MP6yMz47TC09+GgqNPKPIEBFJn9ak45Zl9byzs4OvTilj1oGKjE9LoZFHFBkiIumTSKb48dIGVn/YwRWHl3PCmBKvR/IlhUaeUGSIiKRPvCPFj/5Wz9qGJF87spxjRykyBkqhkQcUGSIi6dPcnuLmv9WzsTHJN6dWMH3/sNcj+Zp+14kPNT/3E9pWLwZ6Rcb7q4hMv1CRISIyCI3tKX6wpJ5NTUmuOkqRkQ4KDR8KjqulfsEFJF57pCcyLBQlMm2e16OJiPhWfVuKm17ZxXvNSa6ureCokYqMdNCpEx8K18yi/B9+Rf2C8wGwcClDLltIuGaWx5OJiPjXL1c0sC3eyTXTKjlseMjrcfKGjmj4UGfTdpqfvLHrg1QnsVlXKjJERAZpWzxF7ciwIiPNFBo+0/uaDAuXUnr6d4m/dGfPNRsiIjJw+tVo6afQ8JE+kRGKMuSyhZSddSOV8++nfsEFig0REck5Cg2f6P0U1siML/W5JiNcM4vK+ffTsXGZx1OKiIj0pYtBfaA/r5MRrpml6zRERCTn6IhGjtOLcYmIiJ8pNHKYIkNERPxOoZGjFBkiIpIPFBo5SJEhIiL5QqGRYxQZIiKSTxQaOUSRISIi+UahkSMUGSIiko8UGjlAkSEiIvlKoeExRYaIiOQzhYaHFBkiIpLvFBoeUWSIiEghUGh4QJEhIiKFIqu/VK3IIBIqzuYi9yrR3unZshUZIiJSSHREI4sUGSIiUmgUGlmiyBARkUKk0MgCRYaIiBQqhUaGKTJERKSQKTQySJEhIiKFTqGRIYoMERERhUZGKDJERES6KDTSTJEhIiLy3xQaaaTIEBER6UuhkSaKDBERkT0pNNJAkSEiIrJ3Co1BUmSIiIjs26BCw8xONbPVZlZnZtekayi/UGSIiIh8vAGHhpkVA7cDpwGTgQvMbHK6Buut/ulbSKxa3Oe2xKrF1D99SyYWt0/Nz/2EttVdc/RExvuriEy/UJEhIiKyF4M5ojEdqHPOrXPOtQMPALPTM1Zf4appfHDHF3piI7FqMR/c8QXCVdMysbh9Co6rpX7BBSRee6QnMiwUJTJtXlbnEBER8YvAIP7uaODdXh9vBo4e3Dh7Fzl4FiMvf5Ctt54CGLhOKAry/s9OG9DjucEMk0pRf9dcKCrGwqUMuWwh4ZpZg3lEERGRvDWY0LC93LbHPtzMLgUuBThw7NgBLyxy8CxKDjqR1tUvEp54LCUHnTDgx0p2pgb8dwHa6/5Mx7qXic26UpEhIvIJeu8Hxg5iP5Bpp4wtYVik2Osx8s5gQmMzcGCvj8cAW3a/k3PuLuAugKOOqh3wwYTEqsW0v/cWlWddR+PiOxh67g+JHDywnXyivXOgY9C2ejGJv95N6enfJf7SnYQOmqnYEBH5GL33A7W1A98PZNrZE2Nej5CXBnONxlKg2syqzCwEzAMeS89YfX10TcbIyx9k6Dk3MfLyB/tcs5EtbasXU7/gAirn30/ZWTdSOf9+6hdc0HOBqIiIiPQ14NBwziWBK4FngZXAg865t9M1WG9t65cy8vIHe45gfHTNRtv6pZlY3D51bFxG5fz7e45ghGtmUTn/fjo2LsvqHCIiIn4xmFMnOOeeAp5K0yz7VHnat/e4LXLwrAGfOhmo0s9dvcdt4ZpZOnUiIiKyD3plUBEREckYhYaIiIhkjEJDREREMkahISIiIhmj0BAREZGMUWiIiIhIxig0REREJGMUGiIiIpIxCg0RERHJGIWGiIiIZIxCQ0RERDJGoSEiIiIZY8657C3MbDuwcZAPMxzYkYZx8oXWR19aH31pffSl9dFXOtbHOOfciP7eeR/7Ab99Xfw0byZn7dfXPquhkQ5mtsw5V+v1HLlC66MvrY++tD760vroK1fWR67M0V9+mjcXZtWpExEREckYhYaIiIhkjB9D4y6vB8gxWh99aX30pfXRl9ZHX7myPnJljv7y07yez+q7azRERETEP/x4RENERER8wjehYWanmtlqM6szs2u8nsdLZnagmS02s5Vm9raZfd3rmXKBmRWb2Wtm9oTXs+QCM6s0s4Vmtqr7/8oxXs/kJTP7Zvf3y1tmdr+ZlXg9UzaZ2W/MbJuZvdXrtqFm9ryZren+c0iGZ/jY7biZhc3s/3V/fomZjc/kPPvSn22smc00swYzW9H99n0vZu01zwYze7N7lmV7+byZ2b91r9s3zGxqtmbzRWiYWTFwO3AaMBm4wMwmezuVp5LAVc65Q4AZwBUFvj4+8nVgpddD5JBfAM845w4GDqeA142ZjQa+BtQ65w4DioF53k6VdfcAp+522zXAC865auCF7o8zop/b8S8Du5xzk4CfAf+aqXk+QX+3sX9yzh3R/XZTdkfcq1nds+zt6aynAdXdb5cCd2RrKF+EBjAdqHPOrXPOtQMPALM9nskzzrmtzrnl3e830bUDGe3tVN4yszHAGcACr2fJBWZWDpwI/BrAOdfunKv3dirPBYCImQWAKLDF43myyjn3EvDhbjfPBn7b/f5vgXMyOEJ/tuO951kInGxmlsGZ9ipPt7Gzgf9wXV4BKs3sgGws2C+hMRp4t9fHm/H/Fz0tug8tHgks8XYSz/0c+DaQ8nqQHDEB2A7c3X06aYGZxbweyivOufeAW4FNwFagwTn3nLdT5YSRzrmt0LVzBfbL4LL6sx3vuY9zLgk0AMMyONMn+oRt7DFm9rqZPW1mh2Z1sD054Dkze9XMLt3L5z3bj/olNPZWtAX/dBkzKwUeAr7hnGv0eh6vmNmZwDbn3Ktez5JDAsBU4A7n3JFACxk8LJ7ruq89mA1UAaOAmJld5O1UBac/2/Gc2tZ/wjZ2OV0vwX04cBvwSLbn281xzrmpdJ0iucLMTtzt856tW7+ExmbgwF4fj6HADnvuzsyCdH0D3Oece9jreTx2HHC2mW2g63DsSWZ2r7cjeW4zsNk599FPYQvpCo9CdQqw3jm33TnXATwMHOvxTLngg48On3f/uS2Dy+rPdrznPt2nuCrY83RPVnzSNtY51+ica+5+/ykgaGbDszxm73m2dP+5DVhE16mq3jzbj/olNJYC1WZWZWYhui7ieszjmTzTfc7y18BK59xPvZ7Ha865a51zY5xz4+n6v/EH51xB/7TqnHsfeNfMarpvOhl4x8ORvLYJmGFm0e7vn5Mp4Itje3kM+Mfu9/8ReDSDy+rPdrz3PHPp+l7O+hGN/mxjzWz/j64fMbPpdO1Pd2Zvyj6zxMys7KP3gc8Bb+12t8eAL3U/+2QGXacPt2ZjvkA2FjJYzrmkmV0JPEvX1eK/cc697fFYXjoO+CLwppmt6L7tO91VLfKRfwbu696orwMu8XgezzjnlpjZQroOdyeB18iBV0zMJjO7H5gJDDezzcD1wI+BB83sy3TF2P/M1PL3tR03s5uAZc65x+jauf/OzOroOpLh1TOD9rqNBcYCOOd+RVcIXW5mSSABzPMiirqNBBZ1d08A+E/n3DNm9tVe8z4FnA7UAXGyuD3QK4OKiIhIxvjl1ImIiIj4kEJDREREMkahISIiIhmj0BAREZGMUWiIiIhIxig0REREJGMUGiIiIpIxCg0RERHJmP8PZlgMXmGGbKUAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "total_cost, pointwise_cost, accumulated_cost= dtw(x,y, 'euclidean')\n", "path = get_path(accumulated_cost)\n", "_ = plot_warped_timeseries(x, y, pointwise_cost, \n", " accumulated_cost, path)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Different length vectors\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD8CAYAAABXe05zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJztnXmQG9d9578Pc4FzYDAzHHKImQFJiRIl8RqMSVmXDymJYzuWr8SJD0W25aysLSfrlFO7cba8TrKuVLJV3iTrdSopbUTJjq1sHJ+J403sWIkt0pIG4PAQD5280HNwLs59YvD2j4cGMRgcfbzufg38PlWsOQg0HjDd337v936/749xzkEQBEH4h4DXAyAIgiDMQcJNEAThM0i4CYIgfAYJN0EQhM8g4SYIgvAZJNwEQRA+g4SbIAjCZ5BwEwRB+AwSboIgCJ9R68RBt27dynft2uXEoQmCICqSEydOTHDOO4081hHh3rVrFxKJhBOHJgiCqEgYY1eMPpZCJQRBED6DhJsgCMJnkHATBEH4DEdi3ARBEF6xtrYGTdOwvLzs9VAKEgwG0dPTg7q6OsvHIOEmCKKi0DQNLS0t2LVrFxhjXg9nA5xzTE5OQtM07N692/JxDIVKGGNhxtg3GWMvMcYuMMbutvyKBEEQDrK8vIyOjg7lRBsAGGPo6OiwvRowOuP+XwD+mXP+K4yxegCNtl6VIAjCQVQUbR0ZYys742aMhQC8GcATAMA5X+WcT9t+Zb8zdw049x2vR0FUK2e/BcyPeT0KwiOMhEpuAjAO4EnG2EnG2F8zxpryH8QYe5QxlmCMJcbHx6UPVDl+9iXg7z8GLM96PRKi2rh+BfjmI8DgV7weCeERRoS7FkA/gL/knMcALAD4bP6DOOePc84Pc84Pd3Yaqtr0N8kB8XV22NtxENWHFhdfZ4a8HQfhGUaEWwOgcc5fyPz8TQghr15SK8DIafH9LF08hMvowk2TBiWJx+M4ePAglpeXsbCwgH379uHs2bNSX6Ps5iTnfJQxlmSM7eWcvwzg5wCclzoKvzF6FlhfEd+TcBNukxVuOvfK8Yf/eA7nh+WGM++IhPD7D+4r+v9HjhzBu9/9bnzuc5/D0tISHnroIezfv1/qGIxmlfwWgK9nMkouAvi41FH4Df3CAWjWQ7jL2jIwckZ8T8KtLJ///Odx5MgRBINBfOlLX5J+fEPCzTk/BeCw9Ff3K9oAEOoG1tfo4iHcZeQ0kF4DIv3A8CCwugjUU3ZuMUrNjJ1kamoK8/PzWFtbw/LyMpqaNuVz2IK8SqygxYGeI0BrN824CXfRV3t3vEd8nRvxbixEUR599FF84QtfwEc+8hH87u/+rvTjU8m7WeauAdNXgTs/CaRTwNRFr0dEVBNaHGiNAt2Z/IDZIaDjZm/HRGzgq1/9Kmpra/HhD38Y6+vruOeee/DMM8/ggQcekPYaJNxm0Wc8PUeA6SvA5We9HQ9RXWhxoPeNIlQH0IpPQR5++GE8/PDDAICamhq88MILZZ5hHgqVmEWLA4E6YMchIBQBlmeAlXmvR0VUAzNDYobdeyfQskP8jvZYqhISbrNoCWDHQaAuSLMewl2GMu0Ae46IDcktbVSEU6WQcJthPSV28nuOiJ9DEfGVZj2EGyQHgJp6oOuA+DlEm+PVCgm3GcbOAWuLOcJNM27CRbQEsKMPqG0QP4e6adJQpZBwmyF3YxLIiTOScBMOk1oFRk7dOPcAseKjc68qIeE2QzIONG0DwlHxc10QaNxKsx7Cea6dBVLLQE9OHVyoG1icENWURFVBwm0GLS529HON0EMREm7CefTVXu+dN36n77FQEU7VQcJtlMUpYOr1jTMegDaICHfQ4iI0p++rALQ5XsWQcBslP76tQzNuwg20uJg0bFjt0ea4ivzVX/0V+vr60NfXh927d+P++++X/hpUOWkULQ6wGiAS2/j71m5g6TqZ/RDOMT8OXL8MHP7Ext/TjLs8/++zwOiLco/ZdQB4x58U/e/HHnsMjz32GNbW1vDAAw/gM5/5jNzXB824jaPFge37gPo8ly991kNxRsIpiq32GpqBYCvNuBXl05/+NB544AE8+OCD0o9NM24jpNcB7QRw8Fc3/1/urIfMfggn0OJAoFbYLORDeyylKTEzdpKnnnoKV65cwZe//GVHjk/CbYTxl4HVuc0zHoDijITzaHFg+/7CoTjaY1GOEydO4Itf/CKeffZZBALOBDUoVGKEQqlYOnoRzozm3niI6mE9BQwNFj73ACrCUZAvf/nLmJqawv3334++vj78xm/8hvTXoBm3EbQBYejTftPm/9PNfujiIZxg/AKwtlB4tQeIFd/8mKisrK13d2xEQZ588knHX4Nm3EbQEuLCyU3FyoXijIRTJAfE1/z6AZ1QBACnzfEqg4S7HEvTwPhLQE+RpSpAZj+Ec2gJYavQtrvw/2c3x2niUE2QcJdj6IT4WmzGA1CckXAOvb9p0dVej/hKE4cNcM69HkJRZIyNhLscWgIAA7rfUPwxZPZDOMHiFDD5avlJA0AThxyCwSAmJyeVFG/OOSYnJxEMBm0dhzYny6HFgW23A8FQ8cdkzX6GC29gEoQV9NVesYwSQJyX9S0k3Dn09PRA0zSMj497PZSCBINB9PT02DoGCXcpOBfCfce7Sz8ud9ZDwk3IQosDLLDZZiEfyuXeQF1dHXbvLrInUCEYEm7G2GUAcwDWAaQ45yXWbhXE5GvA8nTxVCwdKsIhnECLA9vuABpaSj+OhLvqMDPjvp9zPuHYSFQk6xFRYqkKkNkPIZ90Wtgs7H9/+ceGuoHXX3J+TIQy0OZkKbQ40BACtt5a+nFk9kPIZuIVYGWm/GoPEBOHuVFgfc35cRFKYFS4OYAfMsZOMMYeLfQAxtijjLEEYyyh6qaAaZJxkU1ixG+AinAImRRzBCxEazcADsxfc3RIhDoYFe57Oef9AN4B4FOMsTfnP4Bz/jjn/DDn/HBnZ6fUQXrCyrzo6l5qRz8XijMSMtEGgGAY6NhT/rG0x1J1GBJuzvlw5usYgO8AMKhmPmb4JMDTxmY8gBDuGRJuQhJaQuRvG1rt0R5LtVH2rGCMNTHGWvTvAbwNwFmnB+Y5WsYjolThTS6hbmAhY/ZDEHZYngXGLpibNAA0cagijGSVbAfwHSZKbmsBPM05/2dHR6UCWgLouAVobDf2+NyO2207nRsXUfkMDwLgxoU7GAbqGilUUkWUFW7O+UUABVpvVDB64c2eXzD+nNw4Iwk3YYdkZmPS6GqPMdpjqTIoHbAQ1y8DC+OlPSLyyQo3XTyETbQ4sHUvsCVs/DmU1VRVkHAXQkuIr0YzSgAy+yHkoK/2eg2GSXRIuKsKEu5CaHGgrgnovN34c7JmPzTjJmwwdRFYmjIe39YJRcT+SnrdmXERSkHCXQhtAOjuB2pMenBRnJGwi5nCm1xCEYCvizZmRMVDwp3P2hIw+qK5+LYONVQg7KLFxcqt8zZzz6M9lqqChDufkdNAOlXeWKoQrRRnJGySzKz2AjXmnkdFOFUFCXc+5ZqzliLUTWY/hHVWF4Br58yHSQAqe68ySLjz0eJAeCfQvM38c/WO22T2Q1hh+JSIU1sR7sZ2oDZIM+4qgYQ7Hy1hLg0wF33WQ6XHhBV0mwUrwp0twqEZdzVAwp3LjCb6Rlq5cACKMxL20BKi9V1Th7XnUy531UDCnUs2FctiZzYqwiGsohfeWJ00AJSOWkWQcOeiJUSccPsBa88PhkXhDgk3YZaZpNgbsS3cI6LtGVHRkHDnkhwAdvQBtfXWnk9mP4RVkjbi2zqhbiC9Jnx2iIqGhFsntSJyuK2GSXRog4iwgpYAarcA2/dZPwbtsVQNJNw6o2eB9RXrGSU6tEFEWEGLZ2wW6qwfg3K5qwYSbh07qVi5kNkPYZa1ZUmrPRLuaoGEW0eLixNfX25aJWv2Q0U4hEFGz4jYtN1JQ2MHUFNPoZIqgIRbx24qlg7NegizWHUEzCcQAFp2kHBXASTcADB3DZi+Kke4W8mljTBJcgBojQItXfaPRXssVQEJNyBvxgPQjJswj5awH9/WoXTUqoCEGxDCHagDdkjoibyljcx+COPMDgOzmv1sJh3dWphzOccjlISEGxAznh0Hgbqg/WPpRThkNEUYQeZqDxArvvVVYHFSzvEIJSHhXk8Bw4PyLhyA4oyEcbS4yATpsmizkA8V4VQFhoWbMVbDGDvJGPu+kwNynbFzwNqiZOGm6knCIFoiY7PQIOd4ZHRWFZiZcX8awAWnBuIZspeqQKYTzjCZ/RClWV8Dhk/KP/cAYVFMVCyGhJsx1gPglwD8tbPD8YBkHGjaBoSj8o4Zioi+lWT2Q5Ri9EUgtSwvowQAmjqBQC3NuCucWoOP+3MA/wVAi4Nj8Qa98IYxecfM7bjdsl3ecSdeA177V+Cux+Qd00l++kVg6pK153buBe79T3LHoxpaQnyVOeMO1GSKcCQLN+fAz/43sO+9cic5TqElgBNPAlaSa2pqgfs+A7TtlD4sWZQVbsbYuwCMcc5PMMbeWuJxjwJ4FACiUR/8YQFgcQqYeh2IPST3uLkbRN398o77/F8AiaNA34eAYKu84zrBzBDwzBeALe1AXaO5567OAae+Btz5H4C6Lc6MTwW0ONDcBbT2yD1uqFv+5uTEK8CP/ptw0XzLf5Z7bCc49mfAqz8Uq2lTcPHZte0G7vttR4YmAyMz7nsBvJsx9k4AQQAhxtjXOOcb1I5z/jiAxwHg8OHD/kgi1ePbsnJodZwqwklmxjs7or5w66ZdD30T6H6Dueeeehr47n8UZl3tN8kfmypoA0Cv5NUeICYOI6flHlP3C5/zQQhG7ya0733A+x83//w/7hXnnsKUjXFzzn+Pc97DOd8F4IMAnskXbd+ixQFWA0Rico/rhNnPyrzIgAFEwYbq2OkmpK9YKjkXfn4cuH5ZbphER89qklmEo09y/PA3sdtNKBRRfnO3uvO4tbgwrq9vknvcQEB+SuDwSYBnslT8sPGkxa13E6oG24AhPb4tebUHiM8vtQQsXZd3TD0e74e/id1MMR+k85oSbs75v3PO3+XUYFwlvQ5oJ5yZ8QDyi3D00AOg/EmF1CowfMp6tkTLDvG1kotIkgMi+0OGzUI+sotwlmeBsfMAmD/+Jsm4vW5ClSbcFcX4y2ITzDHhlmz2oyWAjj1A83b1L57RF0U3IaufbUOziOErfvHYQosD2/cD9SY3bo0ge8UyPAiAAzvvBZamgLUlOcd1Ci0uwp9WuwmFukWoZX1N7rgkUr3C7dTGpI7MOKO+2dJzxBezASmfbahH/fdplfQ6MDTo7LkHyLvB63/P2x/MHFfhv0tqRTSm6LXZdBkcmBuVNizZVLdwb2lzLmtBN/tZmLB/rOkropin57A/fFBkdBOqZHvSsQvA2oJzq72WLrHpLus8ScaBrXuBbbeLn1U+/0bOiOvOzmfrgz2W6hZu2YU3ucic9ehpgD13+kPQtAEJ/RN9sLKwSra/qcSKyVwCNUK8ZXx+G1Z76gualN6xPjDqqk7hXpoGxl9yZkdfR6bZjxYXRSzb7hDHXZ4R6YEqku0mZPOzDXUDC2Ni6VtpaAmgcaso8nAKWTf4qYsirt17BAjpm8YKp8ppcfvdhEi4FWV4UHx1asYDiBgtIOePr8WBSL8oxVV91jMkqYxbv3gUL4SwhNOrPUDeiiU3ta6+CQiG1T33ADndhIKtQF2T0u+zOoU7GQfA5Jaj5yPL7GdtSWy26Cej6rOB5ECmm9BBe8epVHvSxSlRPu7kpAEQN/iZIfub41ocqG8BOm+7cVxV/yazI6L4xu6kQW+Gouo1hmoVbi0uTkQny8YDAaBFwqxn5LRwGtQzEFSfcWsJ0RTArseI7t+h6vu0ypC+2nNoY1InFBEboMsz9o6jxcUEJ1Ajfm51wAdFFlJ7x6q9x1J9wq1vtthJFzKKjLu2fjJ2Z2Zo2eIUBU8qvZuQjDQ31VcWVtHiAAs4u9oD5NzgVxeB0bMbhVBlQdO7Cdld7QFi4qDq+0Q1Cvfka8DytPMzHkCOcCcHhI2mbg9bFxQbWyoKmsxuQg0tQENI6YvHEtqA2GRucNghWYZwD58E+HqecHeL1FQVN421uKhEldFNKBQRedzrKfvHcoDqE24nOt4UQ0YRjpbYnKGh6qwn+9lKit/6wOzHFOl0xmbB4fg2IGfFUuhaUXXvQXY3oVBE3LTmr8k5nmSqU7gbQqKgwGlae0SHE6tmPzNDwkYz/2R0wm9ZBloi001IkgG9qjcoq0y+CqzMOJuGqtPSBeEtYuPz0+KiQK2p48bvVBXua2cz3YRkCbfae0nVJ9zJuPCHDrjw1u3OeoqtDlTd8U4OyE1zqzThTkooDjFKTV3G18biiiW38CYXVQUtKXklrfgeS3UJt+5p7caFA9g/ybU4UNMgsjQ2HDciZvGri/bGJxO9m5DMMIAPzH5MocVFJlPHHndez86Nr5intaqCJrubkKo3qAzVJdy6p7VT5j75ZBsCWJz1aHEgUsDTWj85VSpO0f2aZX62PjD7MYWWEELoxmoPsCfcxVZ7DS1Ag4LOjXqmmKzV3pY2YQ2r2g0qQ3UJdza1zmQrLas0b7du9pP1tC6wOlBx1qOnucnsJqT4rMcUK3PC09qt1R5gL6WtlKe1aqG6hQng+iW5n222CEfNc6/6hLtjD9DY7s7r2TH7KeVpraKgaQPyuwllb1AVkFkydAIAdyejRCcUAVZmRSMEs5TytFZNuJ3KFCPhVoDsZotLYRIdqxkgpU5GvQhHlVS5bDch2U2XFc1gsEJ+IZUb6Dd4syG1cp7WqgmaFs90E+qTe1xVs7dQTcKd62ntJlZPci0uSuZbuzf/X32jiMGpcvFMvOJMNyEfmP0YRkuIFNQtYfde02pIbeR0aU/rUDcwPybCeSrgVDehUETc9NLrco8rgeoRbtnpQkbRTXnMFuFoA6XL8lUy+3Eqzc0HZj+GKJZa5zTZzXGTn1+50ENrZtN4XoFNY72bkBOfbSgifIIWxuUf2ybVI9xaXMzett3h7utaMfvJelqXE25FBE3vJtRxs/xjtyp0g7LK1EVgcdL91Z5VX5tyntYqhbDGLgCr8w4Jt76XpMh1lkN1CXd3xtPaTawsV414WqsUZ9TT3Jzwl1ZpZWEVJ1IljVDbIOyFzQpPMl76JqOSoGX7mzog3K0KJgFkqA7hzve0dhMrGSBaPONpfaj0cRcngLVle+Ozy/JMppuQQ2EAxc1+DKENAPXNNzyt3cTsjW92WGTxlJs06I/1Gi0ONHY4001IxeytDNUh3LqntdsxRsDajDsZL+9prUqHGKfT3BQ3+zFEvqe1m5gVbiOrg4aQuBGZjZ07gZPdhBo7hE2sCiuLPMoKN2MsyBgbYIydZoydY4z9oRsDk4qbjoD5tHSJwhSjF4/uaV1urKoU4WgJiG5CDhU1KTzrMUTW09rlMImO2c1dbUCIVb7NQi6qbBovXc90E3Loutbfpwo3qDyMzLhXADzAOT8EoA/A2xljdzk7LMkkB4RjXfM29187a/Zj8I8/dl54WpeLh6oiaMkBZ7sJqXKDskohT2s3CUWE//zqgrHHawljntYq7LFoJ8RXJz9bRfdYygo3F+gtxesy/2w2snMZffPMK8yc5JqeWlcm9KCCoGXT3JxsuqzIDcoqsj3KzWLm8zPjaR1SoEOMG92EVFhZFMBQjJsxVsMYOwVgDMCPOOcvODssieie1pJ39M8Pz+JPf/gyuJH8bDPLLS0hMgHKeVo3NItZrpcXz+TrYjbnZLaE4mY/Zcl6Wm/15vXN3OBHXzTuaR2KiDxuLzeNtbjz3YSyRThp517DAoaEm3O+zjnvA9AD4E7G2P78xzDGHmWMJRhjifFxhRLWHZrx/NEPzuNLz7yG5NRS+QebWW6Z2Wzxehnnxt6B4mY/JfGq8CYXMxkgmoE01Nzj8rR3m8bptEibdXolE+oRVaSLk86+jklMZZVwzqcB/DuAtxf4v8c554c554c7OzslDU8Cuqf19hKbLSa5MDKL46+JP+TgVQPdbUIRURJezuxncUr0xDR6oXu9jNMG3Okm5PX7tEoxT2s3MZNzbcbT2utc7slXRSqq05+tCiHJAhjJKulkjIUz328B8PMAXnJ6YNIo5mltgyePX0KwLoAtdTU4aUi4DcYZzcx4AO9notk0N4ezSr1eWVjF6/g2kGku3WFwxm3C09prQXMrU8zr91kEI1fcDgD/xhg7AyAOEeP+vrPDkkQpT2uLTMyv4LunhvErb+hBX28YJ5PT5Z9kdHZi1tPaS7Of1QXg2jl30twUNvspiZbIeFpviiy6i5Eb/Py4OU9rr4twst2EbnH2dRTdHDeSVXKGcx7jnB/knO/nnP93NwYmhVKe1hb52vNXsJpK4+P37kYsGsb54Vksr5URFKMnuRYXntYNzcYGE4pAdIjxoAhnaFDEON0IAyhs9lOS5EBxT2s3CXWX3xw3YrOQy5Y2oK7RO0FLxoVFrtOrvaZOYRnrwxm3f5G8nFpJreNrz1/B/Xs7cXNnM2LRNqTSHC8OlTGQMmL2k06LKkQzY/Vy1uNmGECPuSp28ZRE97T2MkyiY2SPwKyntZdFOG52EwoEhL2y32bcvkaLi9lGIU9rC/zDqWFMzK/ikfuEL0IsKryVy8a5a+uBpm2lO7lMvCy6lZgSbg8FTUu4103I62W5FUbOiGwEt42lChGKAEtTwrOnGMkB857WXu2xDA0C4M4YSxVCQYfKChfuAWkzHs45jh6/jL3bW3DfHpGTu7W5AdH2RgxeMRLnLnOSZ2ewJi50rwSN88xn69KFo8cZFSw9LopeSOVmx5tiZG/wRc4Tq57WRkIwTpD9bF3qHatgVlPlCrcRT2sTPHdxEhdGZvHIfbvAcnbdY9EwBq9eL1+IUy4zIjkABMPmPK2DIaC+xX3hdrubkMJmP0XR4kBrLxDa4fVIyt/gxy4Iz3jTwu3RprGWALbeKuLsbqBPusw2Q3GQyhXu7GaLnKXq0WOX0N5Uj/f0bQy7xHrDGJtbwchMGXvV1jKND6x6WnsxG9DkfrZl8WMRjtc2C7mUy4zQZ7BmQw9Z58Yx62Mzixe9Y0PdoqJ0ccq91yxD5Qp31tP6oO1DXZpYwI9fGsNDb4wiWLfRmrN/p7jrly3ECUVEwcDK/Ob/0z2trcRDPRHuuMgocLObkJ9yuWdHRPGNMsKtb44XOU+0hDVP63IhGCe4fsn9bkIK5nJXrnAb8bQ2yFPHL6E2wPDQ3Zv9Q27rCqGhNoCTV8vEuUt13Lbjae2FoCUHgIjL3YRUatVWDi9thAtR3yTCcEWF26KntReC5kXvWAVzuStTuHVPawk7+jNLa/j7ExoePBTBtpbgpv+vrw3gQHdr+cySbOPWApkldjytW7szHWLWzD/XCno3Ibd29HUUNfspiBYXMXkJqz1pFLvBL05Z97T2QtC0uGjisO12916TZtwuoXtaS7gr/138KhZX1/GJ+4ovI/t3tuHs8CxWUiU2aUptEGlx657WehGOW2Y/I2e86SYU6lbS7KcgRj2t3aRYSG1oUHy18vdsbBc+QG4KmhfdhJq3A6yGZtyOY9TTugyp9TS+8rMruOumduyLFBfVWG8Yq6k0zg+XMJFqKSLcdj2t3Z71eJXmlr3xlciFVwEzntZuUiwX2Y6ntdtFOKuLwLWz7n+2gRpRREfC7TBGPa3L8C/nrmFoegmP3Ft60yYWFRuUJePcdUGgcevmk3zyddGCyerJ6PYyTosD4SjQst2d19PxSxHOtbNAakmNislcQt0ihTO1svH32oA9T2s391hGTnnbO1ahSUOFCncmXchmA9Enjl3Ezo5G/NztpUWqqzWISGuwvOFUoZQ2fSPLajw+Gzt3S7gT3vRPVHCDqCBup0oapVBz6XRatP+yc5Nxc8bt5aavYumolSfcWU9rezOek1evY/DqND52zy7UBMrfAGLRNgxeKbdBWWB2osXteVoHw+6Z/cwMiYvUiwunqVOkdyq0QVSQ5IBxT2s3KbRimXwVWLHpad3aLdIf3dg0Tg6IlEUvugnp164iRTiVJ9xmPa2LcPT4ZbQ01OIDh3sNPT4WDWNoegljsyUKcQott7QBe57WbsYZvZzxBAIiH1mhWU9B9P0Km6s96RSyDUjqe0E2/p6hbiC9BixOWD+GEbzuJhSKiISHZQP2Fi5QgcJtv4Ho8PQSfvDiCD54Zy+aG4zlKmcNp0qFS0IREc9eXRQ/Zz2tbZ6MbsUZ9W5CXfK6CZlC9SKchQlRIKKCsVQ+hfZCZHhau7XHMqOJzCmvPttWtUJ1FSjcA8LTur7J8iG++twVcM7x8N27DD9nX6QVdTWsdAVlfhHO8MmMp7XNk9FN4ZbcTcgUoUjhPHhVUK3wJpeGFhGSyz1PtIR9T2u39lgkZYpZRrE9lsoS7vR6ZrPF+oWzuJrC3w5cxdv3d6G33bjFZbCuBndEWktnlmTv2pmTPCnpZHTD7MeBbkKmUdDsZwNmPa3dJrf6VJantVuC5nU3IcWKcCpLuCdeEU15bcxgv3VCw8zSWsmCm2LEesM4o00jtV5koyb/JNcSQPvN9j2t3TD7uSa/m5BpQt1iDAqZ/WxAi5v3tHaT3MwI3WbBbgVs41Z3No21uLfdhJq3ixCsItbClSXcNpeq6TTHk8cv41BPK/qj5i0j+3e2YXktjZdG5wo/oCXH7EffbJERs3Oj47akTV9bKDbr2YBVT2s3yRVu/Vqx62kdCDifKpdaAUZOe5sbX1MnxJtCJQ5gxdM6h39/ZQwXJxbwyH27N3huGyXWW6YjTn2j8BCeHc54Wo/JORndELTkgKj+lNRNyBJeuNEZZewCsDqvuHB3iw2+1KpcT2un91j0bkJef7YKNVSoLOG26mmd4eixy+gKBfHOA9bM73vatmBrc0PpOLfeNUTmDNaNOKMWd99YKh+VZ9xu9uC0Sm5zaZme1k4LmiqbvgoV4VSOcNvxtAbw0ugsjr02gYfv2Ym6GmsfC2MM/dFw+ZTA2aEcT+t9ll5rA43tQG3QuYtnfkysELy+cJq3Zcx+VBTujKd1+01ej6Q4+g3+ynG5ntbpIN2yAAAcWElEQVRObxqr0k0o1EPCLR07ntYQHW6CdQF8+M6orWHEom24NLGAqYXVwg/QT3ItLs/T2ukOMarMeBQ0+8mi9+BUrfAmFz3Mdf574qusv2d209gh50Y7JmwyCUVE8sNyCTM5l6gc4bbhaT0xv4LvnhrGL/f3INxoL0dZL8Q5lSwS5w71iCqzkTNyT0Yn44zZNLdDzhzfDArFGbMsXc94WisgLqXQQ02vPyPX09rJEJZK3YQUCtWVFW7GWC9j7N8YYxcYY+cYY592Y2CmseFp/fXnr2I1lcYjFlIA8znY04qaACse59b/+Ok1uVVgTgqalgC6DkrpJmQbheKMWYZOiK+qGUvl0xASgr2+KtfT2sk9Fsm9Y23hRvaWQYzMuFMAfodzfjuAuwB8ijHmYrNBA9jwtF5JreNvnr+C+/d24ubOZttDaayvxW1dLcUrKHXhBuR6Wocizpj9rKeEMKkw4wGEeZNqRTjJOMRqz7rNgivoITVA7t8zv7BMJskBdboJKWQtXFa4OecjnPPBzPdzAC4A8DAnrAA2PK3/8fQIJuZXpMy2dWLRME4nZ7CeLiAu+l1btqe1bvazMC7vmIDUbkJSUMzsB4CYNNjxtHYTJ4S7qVOE0pwQNH21p0I3oWwdhvfCbWpnjDG2C0AMwAtODAZPf1BscphlIeNMZjL0wDnH0WOXsHd7C+7bI88qsj/ahq89fxWvjc1jb1fexazvjFtY+g1evY5/OjOCz/3S7ZvzzHPjbzJvCFm/cHMXejrN8bnvnUVyatHSy97c2Yw/eHeBjJvcWY+MHGSdy8eAY38mvGPMcuU54NCvyRuLk+gTB5nC7dSmsd5N6PDH5R7XKrX1QNM2JUIlhoWbMdYM4FsAfptzvmlblTH2KIBHASAatZiZsToPpErYohajtgHY/8umPa2fvziF8yOz+JP3H7BUcFMMvSPO4NXrm4W7oQV442PAbe8yfdw/+qcLOHHlOj569y5EO/LKqnPjjDKX7FpclDWb7Cb0k1fG8fQLV3FbVwu21JuLpV5fWMWzr07gk2+5CTta8+Lqufak2yWkUuokjgJXfmbtmDsOAYc+JG8sTnLgV8QNT7antRMGYHo3IbvVnTIp1gLOZQwJN2OsDkK0v845/3ahx3DOHwfwOAAcPnzYWgDyY9+39DSrPHHsEtqb6vHemNzIz66ORoQb63Dy6nV8qFB64Tv+h+ljnkpO40SmUcPJ5PXSwi0TvSzf5I3tiWOXsD3UgH/8rftM58WfSk7jvX9xHCevTmPHgXzhdmhnPxkHbnkb8KtfkXtc1bj5AfFPNqGIyJSSiV6kppJNbqgbmLro9SgMZZUwAE8AuMA5/1Pnh+QOlycW8OOXruGhN0YRrJPbMZoxhlhvuHQFpUmeOHYJLQ21aKyvKXzcxg6xiSNT0LLdhMwtq7PFTHfvslTMdMeOEOprA4WtA5q7hNmPzBvU3Cgwc1WdOL4fcaJDjBbPdBMy1szEFUIRJYymjFxV9wL4dQAPMMZOZf690+FxOc5TP7uM2gDDQ3fZayhcjP5oG14dm8fM0prtY43MiMYOv3akF4d6woUFLRDIxBklnlTZNDdzgvbkscu2ipnqawM40F3EIremVlzMMoVbxZmd3wh1i7DGUpn2fWZQsZtQKCLava0UMZJzCSNZJcc454xzfpBz3pf59wM3BucUM0tr+EYiiQcPRbAtFHTkNfQ49+lyDYQN8JWficYOH71nF2LRMM4Nz2J5rYD3tuwinOSAmN1GYoafMjG/gu+cGsL7+3vQ1mS9mCnWG8aZoRmspgpsFsrOWdcGhDVplwIpZ35FdqrcwoQISai2CsqGJEdKP85hKqdy0gTfiCexuLqOR+6VlwKYz6HeVjAG2+ESvbHDL+4TjR1i0Tak0hxnh2Y2P1i6oMWFl0qD8fz2p1/IFDPZ/Gxj0TasptK4MFKgvFh2EY6WEHnCdc7cxKsC2Xssqq6CFKmerDrhTq2n8dTPLuONu9uxv9t8laVRWoJ1uHVbC04WK303yLcGhzY0dsj2tix0Q2iVGGdMp0WoxEQa4EpqHV997greurcTe7bZK2bq31nCIlfv5CLjfa6nMj7aigmE38gKmqTMEi0uDMVU6yakSAuzqhPufzl3DUPTS5Y63JglFhUblOlChTgGSKc5njx2CQd7WvGGnSL0srW5AdH2xsKVmaFuUc4sw+xn4mVgZdbUUjVbzCRhJbOjdQu6QkEMFrpBhSIidXRFgtmPnnKmus+I6ugdYqTNuAeALgW7CSlShFN1wn30+CVE2xvxc7dLLFIpQiwaxszSGi5NLlh6/k9eGcfFiQV8Iq+xg35D2ITMZZxJR0C9mOnW7c140y1ycoRj0XDhFYvMeKoqzod+R+amcbabkIKroLqgqGuQtbKwSFUJt54L/fF7d6Em4PxOtd7+zGqcW8+Fzm/sEOsNY3R2GSMzSxufIFvQgmGgY4+hh+vFTI/ca617UCFi0TCSU0sYn8urppVp9qMlxGwxbM/Ol0AmVCfhbzL+ktrdhBQwOqsq4dZzoT9w2J280Js7m9HSUFvccKoEL4/OFc2F7s+ETQav5N0QslWFEmYDybgpf+mjx+UXM+k3vlP5mTmtEuOMfvDR9guyBC05IL6qGr5yulWbAapGuHNzoZsbJDQvMEAgwNBXLKxRhlKNHW7rCqGhUIGKLLMfvZuQwRnP5YkF/OuFa/iI5GKm/d2tqA2wzTe+5i4AzP77XJjMpJwpKhB+Q2/LZ3fTWPVuQgp4wleNcOfmQrtJLNqGl0dnsbCSMvyccrnQ2QKV/JmoLLOfoUEA3HBGiV7M9OuSi5mCdTW4IxLafIOqrRdtzOxePEMKdK6vJEIRYG3B/qaxZm615zqhiCg0WrVmoCaDqhDu/FxoN4lFw0hz4LRmfNZtJBe6f2cbXixUoBKSEGfUdH/p8uY+2WKmg84UM/VH23BGm0FqPf99Sig91lPOTBQYESXQ91js/F2WrouMJpVXQa094uucd0U4VSHc+bnQbtLXUyLvugBGc6FjvWGsptI4n1+gIiPOqMWBzr2Guglli5kc+mxj0TAWV9fx8rW8EmMZccbkgHADrG+ydxxCICPH2aLNgqsoUIRT8cKt50IfysmFdpO2pnrctLXJsHB/32AudCybsZIXRrDbcTvbTaj8haMXM93pYDFTrLdIZo7dG1Q25UxhgfAbMgTNRu9Y18i1FvaIihduPRf6kfvkpamZJRZtw8mr18HLiCnnHE8YzIXuag0i0hosIGg2zX5MdBP64Xnni5l627dga3N94fdpx+xn/GXRsVu1kmo/07IDtjeNkwPqdxPKFuGQcDvGE8cuoSsU3JQL7SaxaBiTC6tITi2VfJzZXOhYtG1zxoXdWY+JgpQnjolipp93sJiJMYa+3rYCKwubZj+annJGM25p1NSJnHir5146LTaMVY5vA6Kac0ubpymBFS3c2Vzoe3Za8oWWRdZfpIxvidlc6Fg0DO36EsbmcroG2Y0zanGgvkXEuEugFzN97B7ni5li0TAuTizg+sLqjV/KuEFtaVc35cyv2AlhTb4mUlH9sAryOJe7ooW7VC60m+zd3oLG+hoMXiku3FZyofU496ncMIJtQRsAet4gUgtLcDRTzPSrR5wvZsoW4mgy32dC7ZQzv2Inx9lPqyAZ2Vs2qFjh1nOhf7m/B+FG677QMqitCeBgT4G86xys5ELvi4RQV8M2GjG1dIkUNyuzgdUF4Nq5sheO28VMB3taEWDAydwbnx2zn6VpUwVGhAnszES1uMhk6rhF7picwOOy94oVbj0X+uMOem6bIRZtw/kiDRBml9fw9xZyoUWBSuvG+G+gRoi3lZNq+KTocl5G0L763BWkXSxmamqoxd6u0MYbX9bsx8KsJ5typngs1Y+EIqIAZ9lCEY6WALoPi25OqhPqBhYngDULzc0l4INPyDwyfaFl0V+iAcLfDSSxYDEXOtYb3lygYnW5amBjcnE1hadfcL+YKRYN41S+Ra7Vjtt+SDnzK1aLU1bmgLHz/lkF6aG6OW9m3RUp3LovtBcFN8Xo6xUblPlZIHZzoft3tmFpbR0vjeakxVmtKkzGgfabgcb2og/5dqaYyamCm2L0R9swt5LC6+PzN35pdVmuxYFttwPBkLwBEgKrew9Dg4ZWe8ogu1WbSSpOuHN9oe/bI8cXWgadLQ3obd+yKR/Zbi50rFfPWMnduLPQCcdA4U06zXH0uGjscNjlYiY9M2fDjc/KyiKdvtGElpCPVUHLrvZ8sgryuBNOxQm3E77QsuiPtm0S7qM2c6F72rZga3PDxji3bvazXKAvZTGmrwILYyWNpX7yyjgujm9u7OAGuzua0LqlbuPnZ8XsZ+p1YHnaPzM7v6FvGptd8WkJYOutIj/aD3hc9l5xwu2EL7Qs9AYIw9OiEOd0choJm7nQjDH051vHWpkNGIhvHz0uGju8Y7/7xUyBANvc+cfW+/RBrrAfqW0Q9sJmBI3zG77ofqGhWWTA0IzbPk75QssiltcR5wlJudCxaBsu5RaoWBW02i2iq3sBXh6dw7OvisYO9bXenDax3ja8MjaH2eU18Qsrs57kANAQErM7whnMpspdvyT6pPpJuAFPi3AqSrid8oWWxe07bjRAkJkLrcd/s51irAiaFge6+0XvwAI8edz7YqZYNAzOgTPJTAjI0g0qIbJJ/JBy5ldCPeb/JoAPhTsip9uUBcqevYyxo4yxMcbYWTcGZBWnfaFloDdAGLx6XWou9MGeVtQE2I04d4vJDjFry8DImaIXzuT8Cr59snhjB7foi4bBWI4jotkb1Mo8MHbOHyXVfsbspnFyAKhvFpk+fsLDIhwj046nALzd4XHYxmlfaFnEomGcHZqVmgvdWF+L27pablRQZs1+DM4GRk4D6bWiwv31bGOHXbbHaodQsA57OptvZJbUbRF+I0YvnmGfpZz5lVBEbACvLhh7vL7aK2OzoByhbrGhn1ot/1jJlBVuzvlPAUy5MBbLuOELLYtYtA2r62npudCxaBinktNY1wtUzMwGSmxMrqTW8TfPX8Fbbu3Enm3eW23GomGcTE7fsMg1E2fU3ycV3jiLGefG1UXg2ll/3kyzRTjud8Jxp2uuQT74+HNYyW/FZYDltTSGppfw+QfvcGBUctENk2TnQvdH2/C156/i9fF53Lq9RZxUk68Ze7I2AISjQMvmlMTvnx7B+NwKPvEBNVYy/dE2fCOh4fLkInZvbcrcoAyuLLSE8MEoUWBESCAbwtKArXtKP3bkFJBO+Vu4Z4eBNnf31aQJN2PsUQCPAkA0am0Dq7mhFnU15oW7uQH44JFeR32hZdHVGsQn33ITfv727VJzofWMlcEr14Vwt/YAl35q7MlaAojetenXemOHW7aVb+zgFrmdf7LCrTf9LYVeYLTnFxweIWGqCMeE/7tyhDLl/R7kcksTbs754wAeB4DDhw9b6pv11x/14R/PAr/3DvmbMLs6GtHWKApUPnhndKPZT6nS7pkhceIVuHBeuCSKmf74/QeUKWbas60ZzQ21GLx6He/v78mY/UyKDda6EpvS1y8DC+NUMekGZjaNtTjQthtoUmNiYAoPy94pJ6pCYIyJFml6swY9zlgu/laiIOWJY6KY6X0KFTPVBBgO9bbeKMRp1d9nmYtHTzmjjBLnqdsCNHaUFzTOhT+OH2fbgJgQ1bd4MuM2kg74twCeA7CXMaYxxj7h/LAIK8R6w3h1bF4UqBid9WhxoKYB6Dqw4ddXJtUtZuqPtuGl0TksrqaMz3q0OFDXBHT6LOXMrxjZHJ/RgPlR/wo3YK9xhA3Khko45x9yYyCEfWLRNnAuSunftDUjaOU8I7QEsOMQULsxP/vJ4+oWM8WiYaynOc5oM7jLaBGONlCywIiQjJEOMfpqr4Q/jvJ4lMtNoZIK4lBva6ZAZdpYh5jUqmiekBc+sNrYwS36enOsA7KmRiUyS9aWgNEXKb7tJkashbU4UBsEtu93Z0xO4FHZOwl3BdESrMOt21pEgYoRs59rLwLrK5sE7Rtx640d3KC9qR67tzaJCkojZj8jpzMpZxTfdo1QBFiaEjfNYmhxIBITBWN+pbUbmBsF1tdcfVkS7gpDd9DjnJefDRTwiEitp/HkcfWLmWK9YQwafZ9JvQktzbhdo1wIK7Uibqh+jm8DmT0WDsxfc/VlSbgrjFg0jJmlNVyaWDAmaC07blxksN/YwS1i0TAm5legXV8qH0/V4kB4J9C8zb0BVjvlhHvkDLC+WgHCnXmfVjpO2YCEu8Loz7WOLbfjrXe8ycnRttvYwS2yhTjJ6fIbRFqC0gDdppxw+7nwJhePGiqQcFcYN3c2oyVToFLS7Gd+DJi+suHCkdHYwS1u62pBsE5Y5JY0+5kZEjnefhcIvxHSN8eLCJoWB1p7bzzOr3hUhEPCXWEEAgx9eqeYUrOeAvHto8cvobmhFh843OPCSO1RWxPAwR4R5y7ZcVuj+LYn1DcBwXAJ4U5Uxt8kGAbqGkm4CfvEom14aXQWy42ZcEehi0cbAAK1QKQPADA6s4x/OiMaO7QE/bHL3x9tw/nhGaw2dYlfFLtB1QaB7Qc2/x/hLMX2WOZGgZmrlbEKYsyTIhwS7gokFg0jzYHzCxmPkmKC1nVAlCcD+Mpzl5HmHB+T0NjBLWLRMNbWOV5dzmS/FHyfcWBH36YCI8IFiglapfX99CCXm4S7Aon1ilZm8ckG8Yv8i2c9BQwNZmc8S6vrePqFq3jbHXIaO7iF/j4TU5kiofz3mVoFhk9VxpLcj7QWETQtDtTUAzsOuj8mJyDhJmQQbqzHTZ1NSAwtFTb7GTsPrC1kZzzfGtQws7SGT7xJ7RTAfLaFgugOb8HAcEo0AM5/n6OZAiPKKPGGULdwZEytbPx9Mg50HRRFYpVAKCLM3NLrrr0kCXeFEuttw8mr18ELpcpll6qHkU5zHD1+CQe65TZ2cIv+nW2ZzJICy/JKSTnzK4U6xKyvFbRZ8DWhCMDXXS3CIeGuUESByiqWg12biwO0BNC4FWjbhZ+8Oo6L4wv4xH27lfHcNkOsN4zhmWWsNHYVuEENiFmfLiCEuxRKlbt2DkgtVVb4yqjRmURIuCuUWFTEf0fRXmAmOiBmPIzh6LFL2B5qwDsP+DOfVn+fY6yjwA0qXlkC4TcKVRVW4irIgyIcEu4KZe/2FjTW1+DiSutGs5/FKdGLsucwXh6dw7OvTuDhu3ehvtafp8K+SCvqawO4vBoWS1Xd7GfuGjB9tXIyF/xIIUHT4kBzlyi+qRRoxk3IQhSotOLFuWbxC/2kGjohvvYcwZPHLyFYF8CH77TWI1QF6msD2B8J4dx8EwAucoSBG30oK2lm5zcaWjZvGuurIB+G5YrS2C5qBWjGTcggFm1D4rrI085ePMkBgAUw2boP3z45hPf396Ctyd85zuJ9ZtIYc99noE40iSC8I9cAbGECmLpYeTdTvQjHRaMpEu4Kpj/ahqF0u/hBFzQtDmzbh6dPTmE1lcYj9+7ybHyy6I+24WoqkxGji4SWEHnCpRoIE86Tm9VUyX0/Xc7lJuGuYPp6wxjlOYKWTgNDJ7De/QZ89fkreMutndizrcXbQUogFg1jlOfcoNZTwPBg5c3s/MgG4Y4DrEZUslYaLrcwI+GuYDpbGrC1vQ0LgUwn6omXgZVZnEzvwfjcirIdbsyyozWILS1tWGGZOOPYOWBtkYRbBULdYtM4tSqymbr2A/X+qc41TCgiTM7SaVdejoS7wumPtmGEt4vZQCYV6/FLW3HLtma8+ZatHo9ODowx9O9sxyg6hHBXYsqZX9E7xMwOZWwWKjBMAogbVDolKkVdgIS7won1hpFMhbF2XQO0OFL1rfjRtWY84tOCm2LEornvMwE0bQPC/s2WqRj0VLnXnwFW5yv3ZppNCXRng5KEu8KJRdswwjuQnhkCknFcqLkV4cYGvC/WXf7JPkJ/n+vTQyKjJK+zD+ERrZnz7Px3xddKLYhyuaECCXeFc/uOEMZZBxpWJsHHX8KP5nbiI2/ciWBdjddDk8qB7lZcYx1oWLoGTL0O9FbozM5v6IJ2+ZgwPGu/ydvxOIWKM27G2NsZYy8zxl5jjH3W6UER8qivDaC2TXS0YeB4Ebfg1+/e6fGo5BOsqwELdYOBi19U6pLcbzSEgPpmgKcrexXU2CGsalURbsZYDYC/APAOAHcA+BBj7A6nB0bIo61rV/b7rjvuxfZQZeY2t24XNyTOAkAk5vFoCAA3ilOAyg2TAEAgALTsUCpUcieA1zjnFznnqwD+L4D3ODssQibdO/cAAF5Nd+PDb64Q8/oCRKI3AwCW228XPQ8JNcgKd4WvglwswjEi3N0Akjk/a5nfET7htltvAwAkm/bhQE+rx6Nxjlv37AUg3iehEKFuAAzofoPXI3EWF3tP1hp4TKGgFN/0IMYeBfAoAESjlIalEts7t+L4rt9EtP9dXg/FUbojEXyv/RF07H0fbvV6MMQNDj8i+ps2+L9KtyS73+RacRHjfJMGb3wAY3cD+APO+S9mfv49AOCc/3Gx5xw+fJgnEgmZ4yQIgqhoGGMnOOeGNgKMhEriAG5hjO1mjNUD+CCAf7AzQIIgCMI6ZUMlnPMUY+w3AfwLgBoARznn5xwfGUEQBFEQIzFucM5/AOAHDo+FIAiCMABVThIEQfgMEm6CIAifQcJNEAThM0i4CYIgfAYJN0EQhM8oW4Bj6aCMjQO4YvHpWwFMSBxOpUGfT3noMyoNfT7l8eIz2sk57zTyQEeE2w6MsYTR6qFqhD6f8tBnVBr6fMqj+mdEoRKCIAifQcJNEAThM1QU7se9HoDi0OdTHvqMSkOfT3mU/oyUi3ETBEEQpVFxxk0QBEGUQBnhpobE5WGMXWaMvcgYO8UYI8NzAIyxo4yxMcbY2ZzftTPGfsQYezXztc3LMXpJkc/nDxhjQ5nz6BRj7J1ejtFLGGO9jLF/Y4xdYIydY4x9OvN7pc8hJYSbGhKb4n7OeZ/KqUou8xSAt+f97rMAfsw5vwXAjzM/VytPYfPnAwB/ljmP+jLun9VKCsDvcM5vB3AXgE9ltEfpc0gJ4QY1JCYswjn/KYCpvF+/B8BXMt9/BcB7XR2UQhT5fIgMnPMRzvlg5vs5ABcgeuoqfQ6pItzUkNgYHMAPGWMnMj0+icJs55yPAOLCBLDN4/GoyG8yxs5kQilKhQG8gjG2C0AMwAtQ/BxSRbgNNSQmcC/nvB8ipPQpxtibvR4Q4Uv+EsDNAPoAjAD4n94Ox3sYY80AvgXgtznns16PpxyqCLcGoDfn5x4Awx6NRVk458OZr2MAvgMRYiI2c40xtgMAMl/HPB6PUnDOr3HO1znnaQD/B1V+HjHG6iBE++uc829nfq30OaSKcFND4jIwxpoYYy369wDeBuBs6WdVLf8A4KOZ7z8K4HsejkU5dEHK8D5U8XnEGGMAngBwgXP+pzn/pfQ5pEwBTiYl6c9xoyHxH3k8JKVgjN0EMcsGRK/Qp+kzAhhjfwvgrRBubtcA/D6A7wL4BoAogKsAPsA5r8oNuiKfz1shwiQcwGUAn9TjudUGY+w+AM8CeBFAOvPr/woR51b2HFJGuAmCIAhjqBIqIQiCIAxCwk0QBOEzSLgJgiB8Bgk3QRCEzyDhJgiC8Bkk3ARBED6DhJsgCMJnkHATBEH4jP8P05APXvIJB+8AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "x = np.array([1, 1, 2, 3, 1, 2, 0, 1, 2, 2, 0])\n", "y = 2*x\n", "z = np.concatenate([y,y])\n", "\n", "plt.plot(x, label='x')\n", "plt.plot(z, label='z')\n", "plt.legend()\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhQAAAIYCAYAAAA1uHWeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzs3Xl8nGW9/vHPPVuSmcme7vte6MZSWlaBUiggiygqKrihVQ6uoKICPz0CiseDK0ePVVARBY4iUNYu7CAUWqAbXaB7urdJ0yQzk9nu3x+Tpi0EmmSWZ2ZyvV+vvmgmM89cHWh78dz383yNtRYRERGRdLicDiAiIiKFT4VCRERE0qZCISIiImlToRAREZG0qVCIiIhI2lQoREREJG0qFCIiIpI2FQoRERFJmwqFiIiIpE2FQkRERNLmycZB6+rq7PDhw7NxaBERccCSJUv2WGv7OJ1D8tcRC4UxphR4Dihpf/4/rbU/eL/XDB8+nMWLF2cmoYiIOM4Ys8npDJLfunKGog2YYa1tMcZ4gReMMY9ba1/OcjYREREpEEcsFDY1jrSl/Utv+w+NKBUREZEOXdqUaYxxG2PeAHYBC6y1i7IbS/JNIrSdhicvpm37M05HERGRPNSlQmGtTVhrjwEGA9OMMRPf+RxjzGxjzGJjzOLdu3dnOqc4KBHeRcO8s2jbMpf9r3wda5NORxIRkTzTrctGrbX7gGeAczv53hxr7VRr7dQ+fbQRuFgkI3tpmD+TeMtG/OOvJt64jLYtDzsdS0RE8swRC4Uxpo8xpqr952XATGB1toOJ85JtjeydfzbxprXUnDWXimm/xF0+ipalN5HaWiMiIpLSlTMUA4CnjTHLgFdJ7aF4JLuxxGnJaBMNC2YR37eS6hkPUjJwJsblITjpe8T2LqFt6zynI4qISB45YqGw1i6z1h5rrZ1srZ1orf1RLoKJc5KxZhoWnk9s7+tUn/lPSgcfXOEqG3UF7sBQWpbpLIWIiBykW2/LYZKxVhoXXkBs9yKqTr+X0iEXHvZ94/YRmHQdsV3/JrrjaYdSiohIvlGhkA42HqbxqYuJ7nqBqtPupmz4Rzp9nn/053GVDaBl6U05TigiIvlKhUIAsPEIjU9dQnT7U1Se8ifKRl72ns81nlKCE79DdMczRHe+kMOUIiKSr1QoBJuI0vjMR2nbNo/Kk/+Af/Snj/ga/7jZuEr70rL05hwkFBGRfKdC0cvZZIzGZy+jrf4RKk78Lf6xV3bpdcbjJzDhWtq2zSO6+5UspxQRkXynQtGL2WScfc9dQdvmB6iY9ksC46/q1uv946/ClNTQskxnKUREejsVil7KJhM0vfh5Ihvvo3zqzwgc/fVuH8PlLSdw9Ddo2/Iwsb1vZCGliIgUChWKXsjaJE0vzSa87q8Ej72Z4MRv9fhYgfFfxXgrdJZCRKSXU6HoZay17H/5asJv3Ulwyv+jfMr1aR3PVVJF4KivEdl0P7HGlRlKKSIihUaFohex1rL/lW8QWvO/BCZeR/CYH2bkuIGjv4HxBGhZ/uOMHE9ERAqPCkUvYa2lefF3CK36NYGjv0n58T/BGJORY7tKa/GP/w8iG+4lvv+tjBxTREQKiwpFL2Ctpfn1G2hd+d/4x19N+Qm3ZaxMHBCYcC24SmhZprMUIiK9kQpFL9Cy9CZal/2YsrFfpGL6rzNeJgDcZf3wj01t9Iw3b8j48UVEJL+pUBS5lmW30vLGDygb/VkqT/pfjMnev/LgxG+DcdO6/KdZew8REclPKhRFrGXlz2l+7XuUjvwklSf/MatlAsAdGIR/zJWE3v4Tidb6rL6XiIjkFxWKItW66naaX72W0mGXUnXqXzAud07eNzDpOrBJWlb8V07eT0RE8oMKRREKrZnD/kVfpWTIxVSd/neMy5Oz9/YEh1E2+tOE1v6BRGhHzt5XREScpUJRZEJv/Ymml75EyeDzqT7jPozLm/MMwUnfg2SU1pX/nfP3FhERZ6hQFJHwur/R9OKV+AaeTfUZ92PcJY7k8FSMpmzEJwmt+R3JyB5HMoiISG6pUBSJ8MZ/sO+FT+PrfwY1Mx7EeEodzROY/H1sPEzLm79wNIeIiOSGCkURiGx6gH3PfgJvn5OpPuthjMfvdCS8VUdROvxSQqt+Q7Kt0ek4IiKSZSoUBS6y5REan/043roTqDn7MVzegNOROgQn34CNNdO66tdORxERkSxToShgbVvn0fj0R/BWT6Zm5uO4vOVORzqMt2YyJUMupvXNX5KM7nc6joiIZJEKRYFq2/4UDU99CE/lUdScMx9XSZXTkTpVPuVGbHQfodW/dTqKiIhkkQpFAYrufJ7GJy/EUz6a2lkLcZXUOB3pPXnrjqdk0Hm0rryNZKzV6TgiIpIlKhQFJrrrJRoWnI8rMJSaWQtxldY5HemIglNuINm2h9Da3zsdRUREsuSIhcIYM8QY87QxZpUxZqUx5uu5CCbvFt3zKg0LzsVV1p/aWU/iLuvndKQu8fU9Gd+AGbSu+Bk2HnY6joiIZEFXzlDEgWuttUcBJwJXG2OOzm4seafY3tdpmH8OrpJaas99Crd/oNORuiU4+UaS4R2E3rrD6SgiIpIFRywU1trt1trX2n/eDKwCBmU7mBwUa1zO3vlnY7zl1Mx6CndgiNORus3X/3S8fU+lZflPsYmo03FERCTDujU1yhgzHDgWWJSNMPJusX2raJg3E+MuoXbW03jKhzsdqUeMMZRPuZGGBbMIr/sL/rFfdDpSWhKtW2l89jL842bjH3WF03GkXTLaROMzHyW2d4nTUQAIjLuK8uNudjqGSE50uVAYY4LA/cA3rLXvuqmAMWY2MBtg6NChGQvYm8Wb1tIwbwZgqJ31FJ6KUU5HSotv4Nl466bRsuwnlI3+rCODyzIhEdrB3nlnkdi/hub9b1E27FKMp8zpWL1eMtZMw4Jzie1dgn/058DlczRPrOF1Wlb8DP+4L+MODHY0i0gudKlQGGO8pMrE36y1/+rsOdbaOcAcgKlTp9qMJeyl4s3r2TtvBtgENec+g6dynNOR0maMITjlBhqfvIjw+r/jH/0ZpyN1WyKym4b5Z5EMbaH8uB/T/Nr3Cb31RwJHfdXpaL1aMtZKw4Lzie15leoz/kHpsEucjkS8ZRO77x9Ny4r/onK67hYrxa8rV3kY4A5glbX259mPJPGWTTQ8MQObCFNzzkK8VcWzB7Zk8AV4ao6hZdmPscmE03G6JRnZS8O8mcSb11N91iMEJn33kH0hbU7H67VsPETjkxcS2/1vqj7w97woEwCe4DDKRn+a0No/kAjtcDqOSNZ15SqPU4ArgBnGmDfaf5yf5Vy9VqK1noZ5M0jGmqg9ZwHemslOR8ooYwzByTeQ2L+WyMZ/OB2ny5Jt+9i74BziTWuomfEQJQPO7NgXkgxtJfT2n52O2CvZeISGpy4huuMZqk69i7IRH3M60mGCk74HySitK//b6SgiWdeVqzxesNYaa+1ka+0x7T8ey0W43iYR2s7eeWeRjOym5ux5eGuPczpSVpQOuwRP1QRalt2MtUmn4xxRMrqfhgXnEm9cTvWMf1Ey6JyO7x3YF9K6/FZsMuZgyt7HJqI0PnMp0W3zqTzlDspGfcrpSO/iqRhN2YhPElrzO5KRPU7HEckq3SkzTyTCu2iYdxbJ0FZqzn4CX59pTkfKGmNcBCdfT3zfSiKbH3Q6zvtKxlpoWHg+sb1LUmvzgw8/OZfaF3IjiZaNhNf9zaGUvY9Nxmh89uO01T9K5Um/xz/mc05Hek/Byddj42Fa3vyF01FEskqFIg8kI3tpmD+TeMtGqmc+iq/vyU5HyrrS4R/DXTGGlqU3Y21+7uE9uDb/MlWn30Pp0Is7fV7J4A+m9oUsL7x9IYXIJuPse+5TtG1+kIrpv8E/brbTkd6Xp2o8pcM/SmjVb0i2NTodRyRrVCgclmxrZO/8s4k3raXmrLmU9D/d6Ug5YVxugpO+T7zhddrq828FzcYjNDx5MdEdz1J12l2UDb/0PZ97cF/IW0Q23pfDlL2PTSbY98JniWz8B+Un3EbgqK84HalLgpOvx8aaaV2lqz2keKlQOCgZbaJhwSzi+1ZSPeNBSgbOdDpSTpWN+hTu4HBalt6UV2cpbKKNxqc/THT7k1Se+ifKRn7yiK85uC/kloLYF1KIrE3S9O8vEFn/N8qP+wnBCdc4HanLvDWTKRlyMa1v/pJk9F238REpCioUDknGmtvX5l+n+sx/Ujr4XKcj5ZxxeQlO+h6xPYuIbl/odBzgwEa/j9G29XEqT/59l++VcXBfyJtENj2Q5ZS9j7WW/S9dRfjtPxM85ocEJ3/X6UjdVj7lRmx0H6HVv3U6ikhWqFA4IBlrpXHhBcR2L6LqjPsoHXKh05EcUzb6M7j8g2lZepPTUdrX5j9J25a5VEy/vdu3B+/YF7Isf/eFFCJrLfsXfY3Q2jkEJn+f4JT/53SkHvHWHU/JoPNoXXkbyVir03FEMk6FIsdsPEzjUxcR3fUCVafdTdmwDzsdyVHGXUJw0neI7nyeth3POpbDJhPse/7TRDbdT/kJPydw1NXdPsbBfSFv0Fb/SBZS9j7WWpoXf4vQ6tsJTLiW8mNvJnWvvcIUnHIjybY9hNb+3ukoIhmnQpFDNh6h8alLiG5/mspT/0zZyMucjpQX/GO+gKu0Hy1LnRmiZG2Sphc/T2TDPZQffyvBCd/s8bHydV9IIbLW0vza92ld+XP8R32V8qk/K+gyAeDrexK+ATNoXfEzbDzsdByRjFKhyJHU2vxHads2j8qT/6AJlYcwnjICE79NdPtCortezul7pzb6fYnwursIHvsjgpOuS+t4B/eFvEp02/wMpeydWt74T1qX34p/3JepmPargi8TBwQn30gyvIPQW3c4HUUko1QociB1E57LaKt/hIoTf4t/7JVOR8o7/nFfxlVSR8uy3O2lSK3Nf5XwW38kOPl6yqfcmJHjHtgX0qyzFD3WsuzHtCz9T8rGfJ6KE/+naMoEgK//6ZoBI0VJhSLLUhv9rqBt8wNUTPslgfFXOR0pL7m8AQITrqGt/jFie5Zk/f2stTS/eg2h1b8lMPHbBI/NXJFJ7Qu5jtiuF4k6uC+kULWsuI3m166nbOTlVJ40B2OK64+pgzNg6gmvu8vpOCIZU1y/U/OMTSZoeuFzRDbeR/nUnxE4+utOR8pr/vFXY3xVNC/L7l4Kay3NS75L65u/xH/U1yk//qcZ/z9g/5grU/tCcnjGpRi0rvoNzYu/Renwj1F56p8wLrfTkbIiNQPmBFqW/UQzYKRoqFBkibVJml6aTXj93QSPvZngxG85HSnvuXwVBI76Om2bHyTWuDxr79Pyxg9oXfFf+MddRcW0X2TldPrBfSFPEd3174wfvxi1rv5f9i/6GiVDL6HqA3djXB6nI2XNwRkwGwiv/7vTcUQyQoUiC6y17H/5asJv3Ulwyv+jfMr1TkcqGIGjv47xltOy9JasHL956c20LL2JsjFXUnHi7Vldm+/YF5IH99jId6G37mT/y1dRMvgCqk+/F+PyOh0p60oGX5CaAbNMM2CkOKhQZJi1lv2vfIPQmv8lMPE6gsf80OlIBcVVUo1//FeIbPw/4vtWZ/TYLcv/i5bXb6Rs1BVUnpz9tfmOfSFbnyC6Z3FW36uQhdbdTdOLX6Bk4Cyqz/wnxu1zOlJOHJwBs5bIxv9zOo5I2lQoMih1E55vE1r1awJHf5Py439SVLvTcyUw4ZsYTxkty3+csWO2rvwlzUuuo3TEZVSe8qecbfQ7sC/EqXts5LvwhvtoeuEz+AacSfWMBzDuEqcj5VRqBszRmgEjRUGFIkOstTS/fgOtK2/DP/5qyk+4TWWih9ylffCP+zLh9X8nvn9d2sdrXf1b9r/6TUqHfYSq0+7K6Ua/jn0hWx4i1rAsZ+9bCCKbHmDfc5/C1/cUqmfMxXjKnI6UcwdnwKwksvlBp+OIpEWFIkNalt5E67IfUzb2i1RM/7XKRJoCE74FxkPL8lvTOk5o7R/Z//LVlAy5kKoP/N2RtfmOfSFZvnqlkES2PELjsx/HWzeN6pmP4vIGnI7kmNLhH0/NgFmqGTBS2FQoMqBl2a20vPEDykZ/lsqT/rforpt3gts/AP/YLxJe9xcSLZt7dIzQ23+h6d+zKRl0LtVn/MOxtfmD+0L+SWzfKkcy5JO2rfNofPojeGuOoebsx3F5y52O5KiDM2Bep63+UafjiPSY/uZLU8vKn9P82vcoHflJKk/+o8pEBgUnfQeAlhU/7fZrw+vvoenFz+MbcBbVZ/7L8bX5jn0hy7Jz9UqhaNv2JA1PfQhP1QRqzp6Hy1fpdKS8oBkwUgz0t18aWlfdTvOr11I67FKqTv1L0d6ExynuwBD8oz9HaO0dJELbuvy68Mb72ff8Ffj6nkbNWQ/lxdp8al/IVUQ23EN8/9tOx3FE247naHjyQjwVY6g9Zz6ukmqnI+WNgzNgXiG6faHTcUR6RIWih0Jr5rB/0VcpGXIxVaf/vahvwuOkwKTrwMZpXfGzLj0/snku+569DG+f6VTPfATj8Wc5YdcFJlwLxkvL8p84HSXnorv+TePC8/EEh1NzzkJcpXVOR8o7B2bA6L4lUqhUKHog9NafaHrpS5QMPp/qM+7rFTfhcYqnfCRlIy+ndc3vSYR3ve9zI/WP0/jMpXhrj6Nm5mO4vMEcpeyajn0hb99FvHmj03FyJrrnVRoWnIfLP5CaWU/iLuvrdKS8lJoB8x2iO5+nTTNgpACpUHRTeN3faHrxSnwDz6b6jPsdX5vvDYKTvw/JNlpX/vw9n9O2bSGNT12Cp3oSNWc/kbdr88FJ3wFjaO3BvpBCFNv7Og3zz8FVWkftrKdw+wc4HSmv+cd8ITUDRmcppACpUHRDeOM/2PfCp/H1P4OaGQ9iPKVOR+oVPJVjKR3+cUKr/4dkZO+7vt+2/RkanrwIT+U4as/O77X5jn0hb91JonWr03GyKtawjL3zZ+LyVlAz6yncgcFOR8p7B2fAPEl010tOxxHpFhWKLopseoB9z34Cb5+TqT5rbl6tzfcGwcnXY+MttK761WGPR3e+SOOTF+ApH0HNOQtwldY6lLDrApO+CzZBSxf3hRSi2L43aZg/E+Muo+bcp/EEhzkdqWD4x30ZU1Kr+5ZIwTlioTDG3GmM2WWMWZGLQPmo4yY8tVOpmflo3q3N9wbe6gmUDvswrW/+mmS0CYDo7kU0LDwPl38QNecUztq8p3wEZaOuILT29yTCO52Ok3HxprU0zDsLjJvac5/GUz7S6UgFxeUNEJxwDW31jxHbs8TpOCJd1pUzFH8Gzs1yjrzVcROe6snta/MVTkfqtYKTb8DGmmhddTuxPUtomD8LV2nf9rX5/k7H65bgpO9BMkrrytucjpJR8f3r2DtvBtgENbOewlMxxulIBck//isYXxXNOkshBeSIhcJa+xzQkIMsB94Pm4znxY+Om/BUHkXNOfNxlVTl6mOQTnhrj6Vk8AW0rryNvQvOwVVSnSoTgUFOR+u2g/tCfksyssfpOBkRb9lEw7wZ2ESEmllP4q06yulIBatjBszmBzUDRgpG3t08IdG0ht0P5s8fRJ6qCdTOWoirpMbpKAIEp9zI3ken4/IPTm30Cw51OlKPBSdfT2TDPbS++UvKjyvs/xNNtNbTMG8Gydh+amc9hbd6ktORCl7g6K/R+ubPaVl2C9Vn3Od0HJEjylihMMbMBmYDDB3a8z/kXaV1BI/Nj0umjLuEstGf1U148oivzzSqZzyIt+bYgi4TcGBfyEdoXfUbAhO+VbBnwBKh7eydN4NkZA81sxbirT3W6UhFwVVSg3/81bQu/ynxff+Jp2q805FE3pfpyn3jjTHDgUestRO7ctCpU6faxYsXp5dMpBeI7X2DPQ8fS/DYH1E+5Uan43RbIryTvU+cQbK1nppz5uPre5LTkYpKIrKb3f8cTumwj1B12l2OZjHGLLHWTnU0hOQ1XTYq4iBv7TGpfSFv/pJkrNnpON2SjOyhYd5Mkq2bqZ75qMpEFqRmwHyZ8Pq/E9+/zuk4Iu+rK5eN3gO8BIwzxtQbY67MfiyR3iM45UZsWwOh1b91OkqXJdsa2Dv/bOLNb1N91sOU9P+A05GKVmDCt8B4euUMGCksXbnK4xPW2gHWWq+1drC19o5cBBPpLXx9puEbeA6tK28jGWt1Os4RJaNNNMyfRXzfm9TMeJCSATOcjlTUDs6A+Qvxlk1OxxF5T1ryEMkD5VNuJBnZTXjtH5yO8r6SsWYaFpxLrHEp1WfeT8mgWU5H6hU6ZsAs7x0zYKQwqVCI5AFfv1Px9T+DlhX/hY1HnI7TqWSslcaFHyS251WqT7+P0iEXOB2p13AHhlA2+rOE3rqDRGib03FEOqVCIZIngpNvIBneTujtO52O8i42HqbxqYuI7nqRqg/8ndJhlzgdqdcJts+AaS3iGTBS2FQoRPKEb8AMvH1OomX5rdhE1Ok4HWw8QsNTHyK6/WmqTr2LshEfczpSr+QpH0nZyMtpXfN7EuFdTscReRcVCpE8YYwhOOVGkq1bCK9z9p4DB9hElMZnLiW6bT6Vp9xB2ahPOR2pVwtO/j4kIkU3A0aKgwqFSB4pGXQu3tqptCz/CTYZdzSLTcZofPbjtNU/SuVJv8c/5nOO5pH2GTAjPk5o9f+QjOx1Oo7IYVQoRPJI6izFDSSa1xPecI9jOWwyzr7nLqdt84NUTP8N/nGzHcsihwtOvh4bb6V11a+cjiJyGBUKkTxTMuRCPNWTaFl6CzaZyPn722SCfS98lsjG/6P8hNsIHPWVnGeQ9+atnkjpsA/T+uavSUabnI4j0kGFQiTPGOMiOPkGEvvXENn0z5y+t7VJmv79RSLr/0b5cT8mOOGanL6/dE1w8g3YWBOtq37jdBSRDioUInmodNhHcFeOp2XpzVibzMl7WmvZ/9JVhN/+E8Fjfkhw8vdy8r7Sfd7aY9tnwPyi4GbASPFSoRDJQ8blJjj5euL7VtC2eW7W389ay/5FXyO0dg6BSd8jOOX/Zf09JT3BKTe0z4D5ndNRRAAVCpG8VTbiMtzlo2hedhPW2qy9j7WW5sXfIrT6dgITrqX8uFswxmTt/SQzfH2m4xt4Nq0rb8PGQ07HEVGhEMlXxuUhOOl7xPe+RtvWx7PyHtZaml/7Pq0rf47/qK9SPvVnKhMFJDUDZhehPJ8BI72DCoVIHisbdQXuwFBalmbnLEXLG/9J6/Jb8Y/9EhXTfqUyUWB8/U7D1+/0vJ4BI72HCoVIHjNuH4FJ3yW2+2Wi25/M6LFblv2YlqX/Sdnoz1Fx0m9VJgpUcMoNJEPbCL39J6ejSC+nQiGS5/yjP4fLP5CWpTdn7JgtK26j+bXrKRt5OZUn/wFj9EdBofINOAtvnxNpzbMZMNL76E8RkTxnPKUEJ36H6M5nie58Pu3jta76Dc2Lv0Xp8I9ReeqfMC53BlKKUw7MgEm0bia8/m6n40gvpkIhUgD8Y7+Iq7QvzUtvSus4rav/l/2LvkbJ0Euo+sDdGJcnQwnFSSWDzsNbezwty37s+AwY6b1UKEQKgPH4CUy4lui2BUR3L+rRMUJv3cn+l6+iZPAFVJ9+L8blzXBKcYoxJnV31eZ1hDfc63Qc6aVUKEQKhH/8VZiSGlp6cJYitO5uml78AiUDZ1F9xj8wbl8WEoqTSoZehKdqIi3LnJkBI6JCIVIgXN5yAkd/k7b6R4ntfa3LrwtvuI+mFz6Db8CZVM94AOMpzWJKcYoxrtSk2qbVRDbd73Qc6YVUKEQKSOCor2K8lbQsu6VLz49seoB9z30KX99TqJ4xF+Mpy3JCcVLpsEtxV4yjZVnuZsCIHKBCIVJAXL5KAkd9lcimfxFrXPG+z41seZjGZz+Ot24a1TMfxeUN5CilOMW43ASnXE+8cTltWx52Oo70MioUIgUmcPQ3MJ7g+56liNQ/QePTl+KtnkLN2Y/j8pbnMKE4qWzEJ3CXj8za3VVF3osKhUiBcZXW4h//H0Q23Ee8ac27vt+27Ukan74ET9XR1JwzH5ev0oGU4pQDM2Bie5fQtvUJp+NIL6JCIVKAAhOuBXcpLct+fNjjbTuepeHJC/FUjKH2nAW4SqodSihOKhv16azOgBHpTJcKhTHmXGPMGmPM28aY72Y7lIi8P3dZXwLjvkR4/d+IN68HILrr3zQu/CCe4HBqzlmIq7TO4ZTilNQMmOuI7X6J6PannI4jvcQRC4Uxxg38D3AecDTwCWPM0dkOJiLvLzDhW2DctCy/lejuV2hYcB4u/0BqZj2Ju6yv0/HEYf7Rn8dVNoCWZZmbASPyfrpyhmIa8La1dr21NgrcC1yc3VgiciTuwCD8Y64k/PafaVgwC1dpHbWznsLtH+B0NMkDxlNKYOK3ie54hujOF5yOI71AVwrFIGDLIV/Xtz8mIg4LTLoOAJe3gppZT+EODHY4keSTwLgv4Srtk/YMGJGu6MpkINPJY+/a5WOMmQ3MBhg6dGiasUSkKzzBYdSe/2/c/sG4/f2djiN5xnj8VJ70e1z+gU5HkV6gK4WiHhhyyNeDgW3vfJK1dg4wB2Dq1KnaViySI766qU5HkDxWOuwSpyNIL9GVJY9XgTHGmBHGGB9wGTA3u7FERESkkBzxDIW1Nm6M+QowD3ADd1prV2Y9mYiIiBSMrix5YK19DHgsy1lERESkQJls3EXNGLMb2JTmYeqAPRmIUyz0eRxOn8fh9HkcTp/H4TLxeQyz1vbJRBgpTlkpFJlgjFlsrdVus3b6PA6nz+Nw+jwOp8/jcPo8JBc0y0NERETSpkIhIiIiacvnQjHH6QB5Rp/H4fR5HE6fx+H0eRxOn4dkXd7uoRAREZHCkc9nKERERKRAqFCIiIhI2lQoREREJG0qFCIiIpI2FQoRERFJW5dmeUj+M54ya3zlTsfAV1GV9jGMST+H252Bg2SIy5V+by+mq7FcrvT+3YT3bifasi9//gXLu9TV1dnhw4c7HUMyZMmSJXu6ctt1FYoiYXzllIz7mNMxGHL2BWkfw+tN/y/g8vKStI+RKYGAL+1jxGKJDCTJD+X+9D6PF37ymQwlkWwZPnw4ixfFCJB+AAAgAElEQVQvdjqGZIgxpkuzubTkISIiImlToRAREZG0qVCIiIhI2lQoREREJG0qFCLiuDO2/pVRTYdv4hvVtJgztv7VoUQi0l0qFCLiuC3Bo7h87fWMaVwE1jKqaTGXr72eLcGjnI4mIl2ky0ZFxHHrKqdy7+gfcOWaa9gUnEjf8EbuHnsL6yqnOh1NRLpIZyhExHHeRJgztv0Vl40zsvkNXur3YZUJkQKjQiEijvImwnx+9TWM3P86bS4/CwZ9npN2/utdeypEJL+pUIiIYw4tE1FXGX8e/zPmD/0Sd4+9hcvXXq9SIVJAVChExBEHy8QbLK2dyZ/H/6xjmWNd5VTuHnsLQ1pWOZxSRLpKmzJFJOcOLRP3jP4hb/SZ9a7nrKucqn0UIgVEZyhEJKe6UiZEpPCoUIhIzvhUJkSKlgqFiOSELxHm6g3XqUz0Ui9ujbB5f5yktU5HkSzRHopiYVzgLXE6BYGAL+1jlJenf4zm5mjax6ipKUv7GACtreln8fncaR+j3J/+59pTqWWO7zKydSl/GnoDrwbOhFDPPpdkUn8hFZqWWJJfv7EfgIDHMKbay/gaL+OqvYyq8lLiNg4nlExQoRCRrHrnnolXA2c6HUlyLOAx/OqMGtY0xlI/GmLcuyZVKN0GRlR6GFftZXyNj3HVXipLdPK8EKlQiEjWdLoBs4dnJqRwGWPoH/DQP+Dh9MGpM38t0eRhBWP+pjCPbggDMCDgZmy1t71keBkYcGOMzmLkOxUKEckKXc0h7yfoc3F8vxKO75daqo0lLBv2x1ndEGVNY4zXdrbxbH0EgHKvYVyNt6NkjKr04tUySd5RoRCRjFOZkO7yug1jq1OlAcBay/bWBKsbYh1nMhbvTJ3d8rpgZOXBfRhjq72U+7RM4jQViiwxxgwB7gL6A0lgjrX2V8aYGuA+YDiwEfiYtbaxk9d/Brih/cubrbV/yUVukXSpTEgmGGMYGPQwMOhhxtDUMklTW5K1jbGOsxiPrA/xUPse3UFBd8cSybhqL/38WibJNRWK7IkD11prXzPGlANLjDELgM8CT1prbzXGfBf4LnDdoS9sLx0/AKYCtv21czsrHiJOO2PrX9kSPIp1lVMPm82xtHamyoRkVGWJixP6l3BC/9QySTRheXvfwX0Yi7a38dSWSMdzxx2yD2N4hQePSwUjm1QossRaux3Y3v7zZmPMKmAQcDFwRvvT/gI8wzsKBTALWGCtbQBoLyLnAvdkPbhIN20JHsXla6/n3tE/4Ixtf+0Y9LWo34ecjiZFzuc2HF3r4+ja1CXRSWupb0mw5sAySUOUV3a0pZ7rgtFVXsYdskwS8GqZJJNUKHLAGDMcOBZYBPRrLxtYa7cbY/p28pJBwJZDvq5vf0wk76yrnMq9o3/A59Z8G5eN0+byHzboSyRXXMYwtNzD0HIPZw9LLZM0RBIdZzDWNMZ4aF2IpAUDDCl3M67ax7gaL+OrvdSVubRMkgYViiwzxgSB+4FvWGv3d/E/1s6e9K67+RhjZgOzAfCVp5FSpOe8iTBnbPsrLhvHAM8PuExlQvJGTambkwa4OWlAKQCReJK39sU7SsYL2yIs2Bxuf+7BZZJxNV6GlXtwa5mky1QossgY4yVVJv5mrf1X+8M7jTED2s9ODAB2dfLSeg4uiwAMJrU0chhr7RxgDoAr0F+3D5ScO3TPRJvLz/MDLuOknf9iXeXxKhWSl0o9LibV+ZhUd3CZZPP+OKsPuSfGS9tTyySlbsOYak/HWYyxVR5KPVomeS8qFFliUqci7gBWWWt/fsi35gKfAW5t/+dDnbx8HvBjY0x1+9fnAN/LYlyRbju0TERdZR3LHOsqj+fytddz99hbVCok77mMYXill+GVXs4dnnpsTzi1D+NAybj/rVYs4DIwrNzTsUQyrsZLTWn6t8UvFioU2XMKcAWw3BjzRvtj3ydVJP7PGHMlsBn4KIAxZirwZWvtF6y1DcaYm4BX21/3owMbNEXywaGXhi6tncmifh/qKA/rKqdy99hbGNKySoVCClJdmZu6QW5OGZRaJgnFkqzdF2Nte8l4ekuYJzamlkn6lLk6lkjGV/sYXO7G1Uv3YahQZIm19gU63wsBcFYnz18MfOGQr+8E7sxOOpGe68p9JlJnKlQmpDj4vS6O6VPCMX1Sl6vGk5aN+w/uw1ixN8YL21LLJH6POey24b1p+JkKhYh0mW5aJQIel2F0lZfRVV4+OCJ1V89d4WTHDbfWNMR4Y/e7h5+lzmT4qCrS4WcqFCLSJSoTIp0zxtDP76afv6xLw8/6+90d98M4qsbLwGBx/FVcHL8KEckqlQmR7unO8LObTq7umGFSyFQoioXHC9UDnU5Bebkv7WN4venvmq6pKUv7GJX+9H8tmZKJLCU9vCugNxHmsje+xdD9b/DQhJtY3/88KtJO03Mu3RcgLx16X5yhQ4c6nCb/dDb8bMmuKD9b3ERzNOlwuswozoUcEcmIVJn4OkP3vcZDE25iRf/znI4kecpaO8daO9VaO7VPnz5Ox8l7xhhqSovrr+Di+tWISMaoTIhId6hQiMi7qEyISHepUIjIYVQmRKQnVChEpIPKhIj0lAqFSC920qY/M6whdYf3g2ViCW/2O1tlQkS6RYVCpBfbVj6Bj6y4jlF7XugoEzF3Ga8N/LDT0USkwOg+FCK92KaaE3jo6B/x8WXX4LJxom4/903+BZtqTnA6mogUGJ2hEOnFvIkwJ2+6C5eNY4BFQz6lMiEiPaJCIdJLHbpnIur289zwL3L81n907KkQEekOFQqRXujQMhFzl3Hf5F/w7KiruH/iT/nIiutUKkSk21QoRHqZQy8NXdlv1mF7JjbVnMD9E3/KwOaVDqcUkUKjTZkivcihZeLBCTezsv+573rOppoTtI9CRLpNZyhEeomulAkRkZ5SoRDpBVQmRPLL/rYkz9ZHnI6RUVryECly3kSYK1Z8g6FNr6tMiDhsZyjBI+tDPL0lTCwJJ/TzMa7G63SsjFChKBLG7cFb08fpGJT7fU5HAKDCn/5v0FKvOwNJoDIDWUp6mMUTD3HRK9cwsOl15h93K/WDzqcy7TTpqUzzvxGfRydWpfBsaIoxd32Il7a14TLwgUGlXDjKz6Bg8fw1XDy/EhE5TKpMXM3AvUuYf9ytrB10vtORRHoVay0r98Z4aF2IZXuilHkMF4z0c/6IMmpKM/M/LPlEhUKkCKlMiDgnaS2LdrQxd12I9U1xKktcfGJcgLOHlRHwFu8ZNhUKkSKjMiHijGjC8mx9hEfWh9gRStDf7+aLk8r5wKBSfG7jdLysU6EQKSIqEyK51xpLMn9TmMc3hmlqSzKy0sM3j6tgWv8SXKb4i8QBKhQiRUJlQiS3GiIJHl0fYuHmCJGEZUofHxeN9DOh1ovpRUXiABUKkQJ13Nt3sqtqIvV10w4pE4t5a+B5KhMiWVTfHOfh9SGe3xohaeHkgSVcNNLP8MriuPyzp1QoRArUrqqJnLfkWuYdcwtT1/2JgXsXE/f4WTnsUqejiRSlNQ0x5q5vZfHOKD4XzBxaxgUj/fT1F98VGz2hQiFSoOrrpjHvmFu48NWv4bJxYp4Aj5zwG+rrpjkdTaRoJK3l9V1R5q4LsboxRsBr+MhoP+cO91NRUrxXbPSECoVIgfLEQ0xd9ydcNo4BXh9xhcqESIbEk5YXt0WYuy5EfUuC2lIXnzk6yIwhpZTq5mqdUqEQKUCH7pmIeQK8PuIKJm+6j61101QqRNIQiSd5cnOERzeE2BtJMqTczdVTyjl5YCkeV+/baNkdKhQiBebQMhH3+DuWObbWTeO8Jdfy+PG3qVSIdFNTW5InNoaYtylMa8xyVI2XL0ws59i+vl55xUZPqFCIFJBDLw19a+B5rBx2aUd5qK+bxuPH30bffStUKES6qLNhXReNCjC2undfsdETKhRZYoy5E7gA2GWtndj+2H3AuPanVAH7rLXHdPLajUAzkADi1tqpOQktea0r95mo15KHSJdsaErN2Hh5e/EO68o1fXLZ82fgduCuAw9Yaz9+4OfGmNuApvd5/ZnW2j1ZSycFRTetEkmftZYVe2PMXdfKsj2xoh/WlWsqFFlirX3OGDO8s++Z1ILcx4AZucwkhckbD3HR619TmRDpoaS1LNrextz1vWtYV66pUDjjNGCntfat9/i+BeYbYyzwe2vtnNxFk3zijYf48GtfZWDDayoTIt10YFjXw+tD7OyFw7pyTYXCGZ8A7nmf759ird1mjOkLLDDGrLbWPvfOJxljZgOzAVyBOry+4thEVOJN/9RjaQaOkSmV/p79e/HEQ5y7+Ov0L7IyUe1P748dXbonR9ISS7JgU5jHN4RoitpeO6wr11QocswY4wE+DBz/Xs+x1m5r/+cuY8wDwDTgXYWi/czFHABP3UiblcDiCE88xLn//g/671nM0yf8lLX9z3M6kkje07AuZ6lQ5N5MYLW1tr6zbxpjAoDLWtvc/vNzgB/lMqA4651lYt2QD0Is6XQskbylYV35QYUiS4wx9wBnAHXGmHrgB9baO4DLeMdyhzFmIPBHa+35QD/ggfY27QH+bq19IpfZxTmdlgkR6ZSGdeUXFYossdZ+4j0e/2wnj20Dzm//+XpgSlbDSV5SmRA5sncO6wp6DR8Z4+fcYRrW5TQVCpE8oDIh8v7eOayrrkzDuvKNCoWIA6asvYNd1RPZ3mf6IWXiVdYNPk9lQuQQkXiSJ7dEeHT9wWFdX5lSwUkDS3TFT55RoRBxwK7qicxcdA1PT/0Jx6y9k/57XiXu8bN6xEedjiaSFzob1vXFSeUc00fDuvKVCoWIA7b3mc7TU3/CrJe+gsvGiXkCzDvpdrb3me50NBFHHTqsK56EqRrWVTBUKEQc4ImHOGbtnbhsHAMsH/1plQnp1TY0xZi7PsRL2zSsq1Dp35RIjh26ZyLmCbB89Kc5ev29bOszTaVCepWDw7pCLNsT1bCuAqdCIZJDh5aJuMffscyxrc80Zi66hoXTf65SIUUvaS2LdrQxd52GdRUTFQqRHDn00tB1g89j9YiPdpSH7X2ms3D6z+nbuEKFQorWgWFdj6wPsSOUYEDAzexJ5ZymYV1FQYVCJAcOv8/ErawbcsG7nrO9z3SVCSla8zaGuP/tEE1tSUZVerjmuApO0LCuoqJCIZJlXSkTIsXu8Y1hmtpSM2maokkW7WhjX1uScTVehpZ7VCyKgAqFSBZ5D4wgV5mQXu62D9SwqTnOmoYYaxpjvLk3xovb2gAo8xjGVnsZ1/5jdJWXUo8KRqFRoRDJEm88xIdf+yr9G15TmZBez+0yjKz0MrLSy3kjUld47A4nOwrG6sYo/1gbxQIuAyMqPIyrOVgyqnXVR95ToSgSxhi8Pudv/FLhTz9Dpd+X9jGaQtG0j9GvsrTHr/UcODPR8BqPTb6FVbWzIBTr8fEisUSPX3vA8Dp/2sfIhMZQPK3Xx5M2Q0nEScYY+vrd9PW7OW1w6vdaayzJ2sZYR8lYsCnMYxvCAPTzuxhX7esoGYOCbi2T5BkVCpEMe+eeiVW1s5yOJFIQAl4Xx/Yt4di+JUCqPG5oirO6vWS8sbuN57ZG2p9rOs5ejKvxMqrSqytFHKZCIZJBnW7ATOPMhEhv5nEZxlR7GVPt5cKRqWWSHaEEaxpiHSXjtV3R9ufCyAoP42p8HUVD48xzS4VCJEN0NYdIdhljGBDwMCDg4YwhZQDsjx5cJlndGOPxjSEeXp96/sCAu+MMxrhqLwMCbg0WyyIVCpEMUJkQcUaFz8XUfiVM7ZdaJokmLOub2jd6NsR4dWcbT9dH2p97YJkktRdjZKVHI9AzSIVCJE0qEyL5w+c2jK/xMb7Gx8WjUrf53taSYE3joSWjBQCvC0ZXHbySZGyNl6Bu/d1jKhQi3TRl7R3sqp7I9j7TD5vNsW7weSoTInnGZQyDyz0MLvdw1tDUMsm+yMGCsaYxxsPrQzzYfvHQ4KCb8e1LJONrfPQpc2mZpItUKES6aVf1RGYuuoanp/6EY9be2THoa/WIjzodTUS6oKrUzfQBbqYPSF2u2pawvL0vdfZibWPqhlsLN6eWSapLXIfdD2N4hQe3lkk6pUIh0k0HBnmd/8JsXDZGzBPomBoqIoWnxG2YUOtjQm3qHjhJa9nSnGBNY7SjZLy8va3juWOqDt50a0yVF7+WSQAVCpEe2d5nOsbGMcDy0Z9WmRApIi5jGFbhYViFh3OGpR7bG0507MFY0xjjX2+FsIABhlV4DruapK6sd97VU4VCpAcG7F6EwbK95liOXn8v2/pMU6kQKWK1ZW5OLnNz8sDUMkk4nuStxnjHWYxn6iPM25S6q2ddmeuwm271luFnKhQi3TRg9yJmLrqGJC521E1l8YSvMXPRNSyc/nOVCpFeoszjYnIfH5P7pJZJEknb64efqVCIdFPfxhXteyi+CBzcU9G3cYUKhUgv1d3hZ8MrPIyv8TKh1sfxfX1FcSWJCoVINy0de+W7HtveZ7rKhIh0eL/hZ8v2RJm3Mcz6pjiPbQjzk1OrGVnp/HDHdKlQiIiI5EDSwrqmOM9vjZCwMLbay4dG+YuiTIAKhYiISFbtDiV4dEOIp7aEaUvAcX19XDzKz/gan9PRMkqFQkREJAs2748zd30rL25rwwCnDCzlwlF+hpYX51+9xfmr6oXcbjcV1UGnY1DpT79xl2bgJjGllaVpH6PW//6nIY0xlHndR3xeJpTE0r+uvcZfHL/dNcxJ8pm1llUNMeauC/H67iglbsN5w8s4f4S/6O9PURx/woiIiDgoaS2Ld0aZu66Vt/bFqfAZPjY2wKxhZQR9veNOmioUIiIiPRRLWJ7fGuHh9SG2tSboW+bi8xOCnDmkDJ+7d51NU6EQERHpplAsycLNYR7bEKaxLcnwCg9fP7aC6f1Leu3wMBUKERGRLmqMJHh8Y5j5m8KE45aJtV7+Y0oFk+q8RXFzqnSoUGSJMeZO4AJgl7V2YvtjPwS+COxuf9r3rbWPdfLac4FfAW7gj9baW3MSWkREOrW9Nc7D60I8uzVCIgnTB5RwcRHdQyITVCiy58/A7cBd73j8F9ba/36vFxlj3MD/AGcD9cCrxpi51to3sxVUumfsyjk01k4+7LE+O16meu8y1k6Y7VAqEcmGt/elrth4ZUcbHhecObiMC0aW0T+gvz7fSZ9IllhrnzPGDO/BS6cBb1tr1wMYY+4FLgZUKPJEY+1kpj//NbBJIFUmpj//NRad9muHk4lIJlhrWbo7ykPrQrzZECPgMXxolJ9zR/ipKukdV2z0hApF7n3FGPNpYDFwrbW28R3fHwRsOeTrekBDIvLI7v4nsui0X3Pawiuo2/UqI96+j0Wn/Zrd/U90OpqIpCGRtLy0vY2560Ns2h+nptTF5UcFmTm0lDKPisSRqFDk1u+AmwDb/s/bgM+/4zmd7eqxnR3MGDMbmA3gqeibuZRyRKny4KJu9xJWTbpaZUKkgLUlLE9vCfPI+hC7w0kGBd1cNbmcUweV6kZq3aBCkUPW2p0Hfm6M+QPwSCdPqweGHPL1YGDbexxvDjAHoLT/2E5Lh2RHnx0vA0n29DmekWv/zu5+J6pUiBSY5miSJzaGmbcxRHPMMrbay2cnlHNcXx+uXn7FRk+oUOSQMWaAtXZ7+5eXACs6edqrwBhjzAhgK3AZ8MkcRZQuOLBnwho3e/qewJtTvtGxh0KlQiT/vXNY1/F9fVxUhMO6ck2LQllijLkHeAkYZ4ypN8ZcCfyXMWa5MWYZcCbwzfbnDjTGPAZgrY0DXwHmAauA/7PWrnTkFyGdqt67LLUB06R++xzYU1G9d5nDyUTkSKIJy7eea+DxjakycebgUi4ZHWB0lS7/TJfOUGSJtfYTnTx8x3s8dxtw/iFfPwa86/4Ukh86uzR0d38teUjvduierqFDhzqc5r15XTB7cjmv74qytjHK0/URnq6P4HPB6Cov42q8jKv2MrbaSyADgwp7ExUKERFJ26F7uqZOnZq3e7qMMZwysJRTBqYmEjdGEqxpjLG6IcaaxhgPrQuRtKnd8UPK3Yyr9nWUjD5lrl5/N8z3o0IhIiK9VnWpmxMHuDlxQKpgROJJ3t4X7ygZL2yLsGBzGICaUhfjqlPlYlyNl2Hlnl47t6MzKhQiIiLtSj0uJtb5mFiX2qCZtJbN+9sLRmOMNQ0xXtrelnqu2zCm2tNeMHyMqfL06vtVqFAUCWPAWyTrfWUZ+HU0huIZSJI/Sovk361IoXEZw/BKL8MrvcwannpsTzjBmvYlktWNMe5/K4QlhAGGV3g6lkjG13ipKXU7GT+nVChERES6oa7MTd0gN6cMSi2ThGJJ1u6Lsba9ZDy9JcwTG1PLJH3KXB1LJOOrfQwudxftPS5UKERERNLg97o4pk8Jx/QpASCetGzaH+9YIlmxN8YL21LLJH6PYewh+zBGV3kpcRdHwVChEBERySCPyzCqysuoKi8fHJEaNrYrnGR1Q5Q17SXjjd1RANwGRlR6Dtns6SvYAWQqFCIiIllkjKGf300/fxmnDy4DoCWaZO2BjZ6NMeZvCvPohtQySX+/u2MfxrhqL4OC7oK4XFWFQkREJMeCPhfH9SvhuH6pZZJYwrJhf5w1DTFWN0Z5bVcbz9ZHUs/1mo4lknHVXkZWevHl4TKJCoVIN41dOYfG2smHPdZnx8tU713W6V00RUSOxOtO7a0YW+3lQvxYa9neevhNt5bsSi2TeFwwqvLglSRjq72U+5xfJlGhEOmmxtrJTH/+a2CTwMFhYYtO+7XDyUSkWBhjGBj0MDDo4cwhqWWSprbUMkmqZER5dEOIuetTzx8YcDP+kMtV+/lzv0yiQiHSTQeGgZ228Arqdr3KiLfv06RREcm6yhIXJ/Qv4YT+qWWSaMLy9r5UwVjbGGPR9jae2pJaJqn0GcbV+Dr2YYyo9ODJ8l09VShEeiBVHlzU7V7CqklXq0yISM753Iaja30cXXvwrp5bWxIdSyRrGqK8siN1uWouhp+pUIj0QJ8dLwNJ9vQ5npFr/87ufpo2KiLOchnDkHIPQ8o9nD0stUzS0D78bE0Ohp+pUIh004E9E9a42dP3BN6c8o2OPRQqFSKST2pK3Zw0wM1J7xh+duCmW4cOP+vvd3PzKdU93uDp/LZQkQJTvXdZagOmSf32ObCnonrvMoeTiYi8vwPDzy4dE+DKiUFOHliCIXXGYnilJ627duoMhUg3dXZp6O7+WvIQkcKwvinG3HUhXt7ehsvAGUNKuWikn4HB9CqBCoWIiEiRs9ayfG+MuetaWb4nRpnHcMFIP+ePKMvYRFQVChERkSKVtJZF29uYuz7E+qY4VSUuPjk+wNlDy/DrKg8RERF5P9GE5Zn6MI+sD7EzlGRAwM3sSeWcNqg0a7ftVqEoEj6fm+FDqpyOwai60rSPUZaB1lzjT/8/7QHlJUd8TlWZmxE17/9rrvanfzoxE/8nUeHzpn2MTKj1x9J6fYlHe8lF3ktLLMn8jWGe2BiiKWoZVenhU8cFOaF/Ca4s3zlThUJERKTA7Q0neHRDiCc3R4gkLMf08XHhSD8Tar05uwW3CoWIiEiBqm+OM3d9iBe2RrDASQNKuGiUn+EVuT8jqUIhIiJSYFY3RJm7LsSSXVF8Lpg5tIwLRvrpm4El1p5SoRARESkASWt5fVeUh9aFWNMYI+g1fGSMn3OH+akocX5vkQqFiIhIHosnLS9uizB3XYj6lgR1ZS4+e3SQM4eUUerJ7Yjy96NCIdJN/Zb8jta+Uw57LLjl3wR2LWXn8Vc5lEpEik04nuTJzREe2xBibyTJ0HI3X5lSwUkDS7I+irwnVChEuqm17xRGPH4V2CSQKhMjHr+KDef9zuFkIlIMmtqSPL4xxPyNYVrjlqNqvHxxUjnH9PHl7IqNnlChEOmmliEns+G83zHmgY8T3LqI2hX3sOG839Ey5GSno4lIAdsZSvDwuhDP1IeJJ2Fqv9QVG2Or8+MeMkeiQiHSA6ny4CK4fTHbp31dZUJEeuydw7o+MDgzw7pyrbDSiuSJ4JZ/A0laBkylbtlfaR50skqFiHRZLoZ15ZoKhUg3HdgzYY2blkHT2Xbitzv2UKhUiMj7eeewrsoSF58YF+DsYWUEMjysK9dUKES6KbBrKRvO+x2jH7ocOLinIrBrqQqFiHQqmrA8Wx/h4fUhdoYSORnWlWsqFCLd1NmloS1DtOQhIu/W+bCuipwM68o1FYosMcbcCVwA7LLWTmx/7GfAhUAUWAd8zlq7r5PXbgSagQQQt9ZOzVVuERFJXz4M68o1FYrs+TNwO3DXIY8tAL5nrY0bY34KfA+47j1ef6a1dk92I4qISCbl07CuXFOhyBJr7XPGmOHveGz+IV++DFyay0wiIpId+TisK9dUKJzzeeC+9/ieBeYbYyzwe2vtnCMdzOt20beyNJP5eqTWn/5/UqFYMu1jDCgvSfsYlb73/z8KA5S43Ud8XiZk4j1KPYW9g/wAd5GeLpbC09mwrkvH+JmVJ8O6ck2FwgHGmOuBOPC393jKKdbabcaYvsACY8xqa+1znRxnNjAbIFA3IGt5RUTkoEIZ1pVrKhQ5Zoz5DKnNmmdZa21nz7HWbmv/5y5jzAPANOBdhaL9zMUcgLqREzo9loiIZEakfVjXowUyrCvXVChyyBhzLqlNmKdba0Pv8ZwA4LLWNrf//BzgRzmMKSIih+gY1rUpTGvMcnSBDOvKNRWKLDHG3AOcAdQZY+qBH5C6qqOE1DIGwMvW2i8bYwYCf7TWng/0Ax5o/74H+Lu19gkHfgkiIr1aZ8O6Lh7lZ0yBDOvKNRWKLLHWfqKTh+94j+duA85v//l6YEoWo4mIyPvY0BTjofZhXW4XnDaoMId15Zo+HZFuqnrldiL9jznssdLNL1C64w32TfuKQ6lEJB3WWla0D+ta1j6s68KRfs4r4OO3BLIAACAASURBVGFduaZCIdJNkf7H0P/h2WBTl7eWbn6B/g/PZseFR7y6V0TyzDuHdVWVuPjk+ABnDy3DX+DDunJNhUKkmyJDT2XHhXMY+H+XUlr/MhXL7mbHhXOIDD3V6Wgi0kW9YVhXrqlQiPRAqjy4KNv6Cg0nflNlQqRAWGt5ZH2Yh9e30hS1jK7y8KnxxTmsK9dUKER6oHTzC0CS8KBpVC79C+Ghp6hUiBSAfW1J7l7dwtgqD18/LsjRNcU7rCvXtEAk0k0H9kxg3EQGn8iOC+fQ/+HZ7SVDRPJZsv0WgGcMKWNCre4jkUkqFCLdVLrjjdQGTJP67XNgT0XpjjccTiYi4hwteYh0U2eXhkaGnqolDxHp1XSGQkRERNKmQiEiIiJpU6EQERGRtGkPRZEwBkq8xXF72Ezcna7UnT9dudKX/iChUk/+/HpERDqjP6VEREQkbSoUIiIikjYVChEREUmbCoWIiIikTYVCRERE0qZCIdJNVa/c/q65HaWbX6DqldsdSiQi4jwVCpFuivQ/JjUczCaBg8PCIv2PcTiZiIhzVChEuiky9FR2fPB3YBOU1r9M/4dns+PCOZrlISK9mm5sJdJdiRiVS/+CAcq2vkLDid9UmRDJc4mk5eUdbcxdFwKgxK2x5ZmmQiHSHYkY/R75EsG3HiPp9bPv+C9RufQvhIeeolIhkoeiCcvTW8I8uiHEzlCSgQE3X5pUzskDS5yOVnRUKES66h1lYvsldxEZeirhoado2UMkz7REk8zbFOaJjSH2Ry2jqzxcflQ5U/v5cBmdncgGFQqRrjikTLSOmsW+47/YUR4iQ09lx4VzKN3xhgqFiMP2hBM8uiHEk5sjtCUs/7+9Ow+vs67zPv7+ni05J3volq5p0o1NSi3IVmjZEWh1xBFlFBWpKLgMooPO88zMxfjM41yPl6OjMzAdQWEUwVGhYQcFamEEhSIu0JYmbeleaNL1JDnb7/njPj1NQ7rlJPdZ8nldV6+cc58753xzkt755P79fvf31NERFrbGOL4xjClIDCsFCpEj6RMm3j7/6+ya86l37NIz+RyFCZEC2rAnRVt7nOc39+CAs8dXsLClism1+jXnF73TIodzFGFCRApnZWeCpe1xVmxPUBGEi6dEuXxqjNGx8ui+XEoUKEQORWFCpChlnOPlbQnaOuKs7kpSEzY+OKOKS6ZEqYnoagiFokAhMpB0krGPfVZhQqSIpDKO5Zt6eKgjzqa9aUZHA3zixGoWTIpqGWgRUKAQ6S+dZOKjn6F6zWMKEyJFoDuV4Zdv9vDo2jidPRmm1Ib43OxazmyqIBhQkCgWChRlImBGZbg8TvVFgvmPfVYM9jn2n5kYwjBRGwnn/RyJTDrv59jVm8z7OYZCXUV+74d+gYwcO3szPL42zpPru9mXcpx4XJhPv6uGU0ZFtGKjCClQiOynORMiRWHrvhQPdcRZtrGHVAZOH1fBwtYY0+rzD+cyfBQoREBhQqQIdOxKsrQ9zotbegkG4LwJlVzREmN8tX5VlQJ9l0QUJkQKxjnHH99OsrR9H3/akSQaMha2xrisOUpDpZZ+lhIFimFiZncBVwDbnXMnZbc1AvcDzcA64C+dc10DfO61wP/K3v26c+5uP2oekRQmRAqib7OudbtTNFQEuGZWFRdOjhIrk/lgI40CxfD5IfA94J4+224FfuWc+4aZ3Zq9/zd9PykbOv4emAs44GUzaxsoeMixq//t9+gZN9u7qmW/y2krTIgMv0M165o3oZKwln6WNAWKYeKc+7WZNffbvAiYn719N/As/QIFcAnwlHOuE8DMngIuBX4yTKWOKD3jZnuNvC6/nbpX7+7TNfT6QpcmUtb6N+uaXh/io8fX8G416yobChT+Guuc2wLgnNtiZmMG2GcCsKHP/Y3ZbTIEeiafw9bLb2f8Ax/F0omDuoaKyPC59blO3urOAPCecRW8b1qMKTUhhYkyokBRfAb63+UG3NFsMbAYoGZ003DWVD7SSepevRtLJwDY+e5PK0yIDIG+x6PJkye/4/FrZlXz2229rOpM8uLWXl7c2ktl0JjeEGJmQ5iZjRGm14eIhjR/olQpUPhrm5k1Zc9ONAHbB9hnIweGRQAm4g2NvINzbgmwBGDstJMGDB3SR585E94wx6epe/VuuiefrVAhkqe+x6O5c+e+43h05vhKzhxfCXgtxld1JlnVlWRlV5KfvxHHEceA5toQMxvDzGoIM7MxTKNWepQMBQp/tQHXAt/Iflw6wD5PAP9kZg3Z+xcDX/WnvDLWL0zsH+bonny2N6fiyiUKFSI+GRUNMmpCkLMneAEjnsywemeS1Z1ewHhmQzePr+sGYHQ0kD2DEWZWQ4SJNUENkxQpBYphYmY/wTvTMMrMNuKt3PgG8FMzuw54E/hgdt+5wA3OuU855zrN7B+B32Wf6rb9EzRlkPqt5tj57utz4aFn8jlsvXIJlVt/r0AhUiCxcIDZoyuYPboC8JqArd+dYmVXklWdSf60I8lzm3u9fUPGjIZwLmRMqw+rMViRUKAYJs65Dx/ioQsG2Pcl4FN97t8F3DVMpY0sR3GdiZ7J5yhMiBSRUMBorQ/TWh/m8qnexa+2d2dY2ZlgVTZk/P4tbx5U0GBqXXYeRnYuRn2F5mEUggKFlC9dtEqkLJgZY2NBxsainDcxCnjLUFdn52Cs6kry5PpuHlnrDZOMiwWZ2RjOhYwJ1UE1E/OBAoWUJ4UJkbJWHQkwZ2wFc8Z6wyTJtGPt7hSrOpOs7EqwYnsvyzb2AFATzg6TZENGa11YF9EaBgoUUn7SScY+dpPChMgIEg56oWFGQ5grieGcY8u+tLeSJLui5OXt3jBJKACtdV64mNXofU5NRMMk+VKgkPKSTtLy+Gepbn9cYUJkBDMzxleHGF8dYsEkb5hkV683TOKFjASPrI3T1uHtP74qyKzGAyFjbEzDJMdKgaJMBAyiRdBQpzYSzvs56ioG+RzpJI1LbyLa/jhb599G5ymfgHR60HWMjlYO+nP72p1I5v0cvXl8HfsN+n0dYrt683s/0hldckUGp64iwGnjKjhtnDdMkkg71uz0AsbqriQvbunl6Q3eMEldxJjZGMkFjObaEKGAAsbhKFBIeUgnaVx6PdHVj7Lzwv9D58kfL3RFIlLkIkHjhOMinHBcBICMc2zam84NkazqTPDbrd5y1UgAptUfmIcxoyFMVRH8EVdMFCik9PULE/vmXg95/hUsIiNPwIxJNSEm1YS4aIo3TNLZk84tVV3VlWRpe5yM83okTKoJMrMhkgsZo6OBET1MokAhpW2gMCEiMkQaK4Oc2RTkzCZvCLQnlWHNzgMX3Xpucw9Pvdmd3TfQ53oYYabUhAiOoGESBQopXQoTIuKzylCAk0ZFOGnUgWGSN3encn1JVnUm+c0Wb5hkpDU/U6CQklH9wndJNJ1KYso5B4WJ7mmXKEyISEEEzGiuC9NcF+aSZm/b/uZn+y+6NVKanylQSMlINJ1K44PX07nwdqpf+SHR1V6jr72nfbrQpYmI5BxL87OxsSBfO72OcVWl/+u49L8CGTESU86h833/yaifXo2lE2TCMXZc9SPvjIWISJE6VPOzF7f2srQ9zqa96bIIFOU7mCNlKTHpTCztXe1u72k3KEyISMnZ3/zsjKaKQpcypBQopKRE1j8PQO+ks6h65W4i658rcEUiIgIKFFJCIuufo7FtMQC9zfPofN9/0vjg9QoVIiJFQIFCSkZkyyt0LlySu79/TkVkyysFrEpERECTMqWE7D3jc5A5uKdFYso5mkchIlIEdIZCRERE8qZAISIiInlToBAREZG8KVCIiIhI3jQps0yEg8aEukihy6CuIpz3c0QCR762fdACh91vdDT/6+NXRYrpGvv5v68VwfL4+2EkdW8UKSXlcYQRERGRglKgEBERkbwpUIiIiEjeFChEREQkbwoUUjKiv/lXwv36doTXLSf6m38tUEUiIoO3dV/6yDuVEK3ykJKRajqV2geuy90Pr1tO7QPXsfv9dxawKhGRo5dxjle2J2hrj7OyK0l12GiqKqYVZYOnQCElI9k8j90L76D+/g8SXvdror9bwu7330myeV6hSxMROaxUxvH85h7a2uNs3JvmuMoA155QzfmTKqkMlcdggQKFlI5MispXfwRA5M3n2Xf2LQoTIlLUelIZfrWhh0c64uzoyTCpJsiNp9Rw1vhKQmV2TRUFCikNmRQ1Dy6mcuVSXDhG/PTPEl3xA5JTzlGoEJGis7s3w2Pr4jyxvpt9ScfxjWE+dVINp46JYFZeQWI/BQopfv3CxK4P3kuyeR7JKefk5lAoVIhIMdgWT/NwR5xnNnSTzMBpYyMsbK1iRkP+V7stdgoUUtz6hIneaRfTffpnc+Eh2TyP3e+/k9CWVxQoRKSg1u5K0tYR5zebewkYnDuhkitbY0yoHjm/ZkfOVyqlp0+Y2HvBP9L9ns++Y5dk8zyFCREpCOccf96RZGl7nD+8nSAaMq5oifHeqVEaK8tj5caxUKDwmZnNBO7vs6kF+Dvn3Lf77DMfWAqszW76hXPuNt+KLAZHESZERAoh4xwvbu2lrT1Ox64UdRUBPjyzioumRKkKl8eKjcFQoPCZc24VMBvAzILAJuCBAXZd7py7ws/aiobChIgUoUTasWxjDw93xNkaTzMuFuT6k2s4d0IlkWB5TrQ8FgoUhXUB0O6cW1/oQopGJkVN22cUJkSkaOxLZnhyfTePretmV2+GlroQfz2nltPHVRAo0xUbg6FAUVhXAz85xGNnmtmrwGbgFufcn/0rq0AyKRrabqByZZvChIgUXGdPmkc64vzyzR560o5TRkdY2BLjxOPCZbv0Mx8KFAViZhFgIfDVAR5eAUxxzu01s/cCDwLTB3iOxcBigFFNExhXVTGMFR+dmvAgl0ZlUkSXfpbwyjZ6Lvo6mTNuJJ+vpjKS/4SocBGdwgwFR+64bH9BHchlmG3ck+KhjjjLN/WQcXDW+AoWtsRoriv/pZ/5UKAonMuAFc65bf0fcM7t7nP7UTP7dzMb5Zx7u99+S4AlAK0nnOKGu+Bhk0kR/cX1hF9/kJ6Lvk7ijBsLXZGIjECru5Isbd/HS9sSRAJw4eQoV7TEGBMbeSs2BkOBonA+zCGGO8xsHLDNOefM7HS8rrA7/CzONwoTIlJAuWZdHXFWdnrNuj4wPcalU2LUVujM4LFQoCgAM4sBFwGf7rPtBgDn3B3AVcBnzCwFdANXO+dK9wzEoShMiEiBjIRmXX5ToCgA51wcOK7ftjv63P4e8D2/6/KVwoSIFMBIatblNwUK8UXkf75Devwc0s3zDgoTyemXKEyIyLDb3Zvh8XVxHu/TrOv6k2uYPbp8m3X5TYFCfJEeP4fozz9B9/u/T2TF3YRffxAXjilMiMiw6tusK5WBuSOoWZffFCjEF+nmeXR/4AfE7r0KSye8FuRX3+edsRARGWJq1uU/vbPim/TE07F0AoDEGTcqTIjIkFKzrsJSoBDfBNc/B0CqeR7hl+8i1TxPoUJE8qZmXcVBgUJ8EVy3nOgDiwFItSwgPe/L3pyKD/xAoUJEBkXNuoqLAoX4Irh5Bd3vX0LVvVcBB+ZUBDevUKAQkWPSv1lXa12Im+fUcpqadRWUAoX4InHWFyDVc9C2tIY8ROQYqFlXcVOgEBGRota3WZcDzmxSs65ipEAhIiJFSc26SosChYiIFI1cs672OCu71KyrlChQiIhIwfVv1jUqur9ZV5TKkOZHlAIFijIRMKgIlsdpwNAQLPcKD8FzdCfSeT8HQDSS//clmc6/2WzPEHw9e5LJvJ+jJqxxbzlgoGZdN51Sy5njK9Ssq8QoUIiIiO9292Z4bF2cJ9Ssq2woUIiIiG/UrKt8KVCIiMiwe3N3igfa9x1o1jWxkitb1KyrnOg7KSIiw+4fX+xid8KbCzR7dIQTGsOEAoZzTkMcZUKBQnwR+Z/vkB570kHbguuWE9y8wruKpoiUtb89vZ6XtydY3ZXktc4kL2/3Og83VASY0RBmVmOYmQ1hmmtDBDUZsyQpUIgv0uPnEP3Zx3P3g+uW55qDiUj5a64L565smXGODXvSrOpKsLIzyequJC9u7QWgImhMrw8xMxswpteHialjaElQoBBfpJvn5ZqDhTqeIfLCv6nTqMgIFTBjSm2IKbUhLp7ibdvRnWZVV5JVXUlWdib5xRtxHGDAlNoQMxvCuZAxKloeS+TLjQKF+CY95WwAQuuW0zvvywoTIpJzXDTIWdEgZ42vBKA7leGNrlTuLMazG3t4Yn03AKOiAS9gZEPG5JqQuowWAQUK8U1w/fMApJrnEX75LlLqNioihxANBXjX6AjvGh0BIJ1xrN+TYlWndxbjtR1Jnt/cm93XmLE/YDSEmVYf1tU1C0CBQnwRXLec6AOLAUi1LCA978u5ORQKFSJyJMGA0VIXpqUuzGVTwTnHW92ZXMBY2ZXgv1cncHhXDm6uDeUmes5sCNNQqWGS4aZAIb4Ibl6Rm0MB2TkVH/gBwc0rFChE5JiZGWNiQcbEgsyb6A2T7EtmWN2VzIWMp9Z38+hab5hkbCzAjIZILmRMqA5qmGSIKVCILxJnfQFSPQdtS2vIQ0SGUFU4wKljKjh1TAXgNRxbuyuVm+j5h7d6Wb6pJ7uvHTQPo7UuTGQIegCNZAoUIiJSlkIBY3pDmOkNYa5o8YZJtsbTrOpMsjJ7JmNF9noYoQC01IaY2RjJBQ21Sz82ChQiIjIimBlNVSGaqkLMnxQFYHfiwDDJyq4kj62L81CHt//4quBBy1WbqoK6qudhKFCIiMiIVRsJMHdsBXPHesMkibSjY9eB62H8blsvz2zsye67f5gkwszGMC11IbVY70OBQkREJCsSNGY1RpjVGGFRq3dVz817D77o1u+27QUgHIBp9QdWksxoDFM9gq/qqUBRJsyMylDhf5ArI4dZmmXeY6FggMBh9ose7jmOUncinfdzDEUdMDS1pNJuCCrJX004/xbTe5LJvD4/7YrjvZCRIWDGxJoQE2tCXDDZGybZ2XMgYKzqSvJQR5wHsz+WE6uDB5arNkYYEw2MmGESBQoREZFjUF8Z5D1NQd7T5C1X7U071uxM5vqSPL+5l1++6Q2TjKTmZwoUIiIieagIGiceF+HE47yreh5V87PsZM9yan6mQCEiIjKEjqr52Rqv+Vk4AN84p5GJNaX/67j0v4ISZGbrgD1AGkg55+b2e9yA7wDvBeLAx51zK/yuU0TkaJnZYmAxwOTJkwtcTfHZ3/xsTCxIV0+G9btThAMwf2KUMbHyuCy4AkXhLHDOvX2Ixy4Dpmf/vQe4PfuxZAWWfxs39qSDtlnHr7FNK8jM+2KBqhKRoeKcWwIsAZg7d65mzvbhnOPVtxIsbY/zWmeSqpDxvtYYl06NUV9GF89SoChOi4B7nHMOeMHM6s2syTm3pdCFDZabMIfQ/dfm7lvHrwndfy2pD91dwKpERIZPOuP4zZZe2jrirN+dorEywF8dX82FkyuJFsGqvKGmQFEYDnjSzBzwH9lk39cEYEOf+xuz20o3ULScS+oD/0n4v/6CwJqnsee/S+pDd+Nazi10aSIiQ6o37XhmQzcPd8R5qzvD+KogN7yrhnkTKsv6QlgKFIVxtnNus5mNAZ4ys5XOuV/3eXygn7h3nELsO2Y5pmni8FQ6hNzUcwAIrP016flfUZgQkbKyJ5HhifXdPL4uzp6EY0ZDmI+fWMOcMZER0dlUgaIAnHObsx+3m9kDwOlA30CxEZjU5/5EYPMAz5Mbs5xx0uyiH7O0tc8BkJl6LoHf3klm6rkKFSJS8t6Kp3lkbZynN3TTm4Y5YyIsao0xqzFS6NJ8pUDhMzOrAgLOuT3Z2xcDt/XbrQ24yczuw5uMuauU509Ads7Ez68HIDPtfNz8r+TmUChUiEgpenN3iraOfTy/uRcDzh5fyZWtMSaXwRLQwRiZX3VhjQUeyF6KNQTc65x73MxuAHDO3QE8irdkdA3estFPFKjWIWObVuTmUEB2TsWH7sY2rVCgEJGS4Zzj9c4kbe1xXnkrQUXQuLQ5yuVTY4yKlsfyz8FSoPCZc64DOGWA7Xf0ue2AG/2sa7hl5n0Rkj0HbXMtGvIQkdKQcY6XtiVY2r6PNTtT1EaMv5xRxSVTolRHym/FxmAoUIiIiBxCMu1YvqmHhzribN6XZkw0wCdPrGbBpCiRYPlPtDwWChQiIiL9xJMZnnqzm8fWdtPVm6G5NsTnT63ljHEVZdvcK18KFCIiIlldPWkeW9fNk+u76U45TjouzGdPqeXkUeER04Z8sBQoRERkxNu8N8XDHXGWbeohnYH3NFWwsCVGa3240KWVDAWKMhEyo6Gy8Gueo5FDz3J2FiQDhING4DD7DcWV5Goqy+tHO5ku+suM+CaovxJlCK3Z6a3Y+O3WXkIBmD+xkitbYoyrKq9jiB/0jomIyIjSv1lXLGQsao1xWXOU+sqRvfQzHwoUIiIyYqQyjn/4TRdv7EwBsGBSJR87vppYWEs/86VAISIiI0bAoKkqxNZ4mj0JxzMbenhpay8zG8PMaAgzqyFCS12IsJaEHjMFChERGTECZtw4uxbnHFv2pVnVlWRVZ5KVXUle2pYA9hEOQEtdmFmNYWY2eEGjRhevOiIFCvFF5tlvQdPJB21z7ctwG14mMP/mwhQlIiOWmTG+OsT46hALJkUB2NWbYXVXMhsyEjzcEWdpdj70hOogsxrCzMyGjLGxoJaR9qNAIb6wSe8m86OP5u679mVkfvwxAtfcU8CqREQOqKsIcNq4Ck4bVwFAIu1o35VkZWeS1V1JXtjSy6829OT2ndkQzoWM5trQkKxQK2UKFOILaz0Pu/pO3F3vw73xNG75dwlccw/Wel6hSxMRGVAkaBzfGOH4bBvyjHNs2usNk3ghI8Fvt/Z6+wZgWv3BwyQjbaKnAoX4xlrm4QDal2EX3KowISIlJWDGpJoQk2pCXDjZGybp6jkQMFZ1JXmwPU7GgQGTakLeWYxsyBgVDZT1MIkChfjGdSz3brSeh3vh+7jWcxUqRKSkNVQGOaMpyBlNlQD0pDKs2ZnKhYznNvfw1JvdABxXGWBGQzgXMqbUhgiUUcBQoBBfuPZluPuuA8Cmn49d8De5ORQKFSJSLipDAU4aFeGkUQeGSd7c7QWM/SHjN1u8YZLKoDGjIcSMhgizGsNMrw9RGSrdYRIFCvGF2/Bybg4FeHMqAtfc421XoBCRMhUwo7kuTHNdmEuavW1vd6dZlR0iWdmV5Odv7MPhXSNjSk0oN0QyszFMYwlduVOBQnwRmH8zLtlD344U1nqewoSIjDijokFGTQhy9gRvmCSezLB6Z5LV2ZDx9IZuHlvnDZOMjgayASPCzIYwE2uCRTtMokAhIiJSQLFwgNmjK5g92luumso41u9OsTJ70a0/vp1k+SZvmKQqZEzvM9GztT5MRZFc1VOBQkREpIiEAkZrvRcWLp/qNTPb3p1hVWciFzLuW5UAIGgwtW7/ahLvLEZdRWHmYShQiIiIFDEzY2wsyNhYlHMnestV9ya8q3quzE72fHJ9N4+s9YZJmqqC2b4k3jyM8VX+XNVTgUJERKTEVEcCzBlbwZyx3jBJMu1YuzuV7UuSYMX2XpZt9K7qWRM2X5qfKVCUCzNCwdJdbtTXUPygdycz+T9HIp33cwD0DMHz9Kbz/3qKxe5EMq/PT2XK570QGSrhoDEje4XOK4kVpPmZAoWIiEiZOVLzs5UDND9716gIH5lVTWSQf9QpUIiIiIwA/Zuf9aYd/716Hw93xNm0N83eRA9/Mb1KgUJERESOLJVxPL+5h7b2OBv3phkVDXDF1BgLJkWpDA1+yFmBQkREZAToTmX41Zs9PLo2zo6eDJNrgtx0Si1njq8YktbrChTii8yz34Kmkw/a5tqX4Ta8TGD+zYUpSkRkBNjVm+GxdXGeXNfNvpTjhMYw159cw+zRkSFdTqpAIb6wSe8m86OP5u679mW55mAiIjL0tsXTPNQe59mN3aQyMHdsBYtaY0xvCA/L6ylQiC+s9bxcczD3xtO45d9Vp1ERkWGwdleSpe1xXtjSSzAA8yZUsrAlxvjq4f2Vr0AhvrGWeV5zsPZl2AW3KkyIiAwR5xx/3JHkofZ9/OHtJNGQcWVLjMumRn3rWKpAIb5xHcu9G63n4V74Pq71XIUKEZE8ZJzjxS29tHXE6diVor4iwEdmVXHR5CixsL8XO1SgEF+49mW4+64DwKafj13wN7k5FAoVIiLHJpF2LNvYw0MdcbbF0zRVBVl8cg3zJlQO+joS+VKgEF+4DS/n5lCAN6cicM093nYFChGRo7I3meGp9d08tjbOroSjtS7ENXNqOW1cBQEfGoAdjgKF+CIw/2ZcsgfXZ5u1nqcwISJyFHZ0p3l0bZxfvtlDT9oxe3SEK1tinHhc2JdOokdDgcJnZjYJuAcYB2SAJc657/TbZz6wFFib3fQL59xtftYpIiKFt3FPiraOOM9t8v4gO7OpgoWtMZprh2fpZz4UKPyXAr7knFthZjXAy2b2lHPutX77LXfOXVGA+kREpMBWdSZp69jHS9sSRAJw4eQoV7TEGBPzZ8XGYChQ+Mw5twXYkr29x8xeByYA/QOFiIiMIBnneGV7grb2OCu7klSHjaumx7hkSozaCn9XbAyGAkUBmVkzcCrw4gAPn2lmrwKbgVucc3/2sTQREfHJQM26Pn5Cdd7NuvymQFEgZlYN/Bz4onNud7+HVwBTnHN7zey9wIPA9AGeYzGwGGDCxMmMqa0Y5qqPLHSY5UoumWYHEI0EiVUdevwvHByKJJ7M+xmiQ7SGuzuS/ynKVDozBJUUh6pEfu9HOFD8f6mJHI2ebLOuR4apWZffFCgKwMzCeGHix865X/R/vG/AcM49amb/bmajnHNv99tvCbAE4JRT3+0QNPgXnwAAC0ZJREFUEZGit6s3w+Pr4jyxvpt9yeFr1uU3BQqfmffTcifwunPuW4fYZxywzTnnzOx0IADs8LFMEREZYtviaR7uiPPMhgPNuha2xpgxTM26/KZA4b+zgY8CfzSz32e3fQ2YDOCcuwO4CviMmaWAbuBq55zOQIiIlKBCNevyW3l9NSXAOfcccNhzWs657wHf86ciEREZas45/rQjSVsBm3X5TYFCfBF/6puEJrzroG2J1c+SWv8SsYtuKVBVIiJDq5iadflNgUJ8EZoylz13fjh3P7H6Wfbc9RFqPnlvAasSERkaxdisy28KFOKLyIz51Hzsh+y+/UoSK39J99PfpuaT9xKZMb/QpYmIDFr/Zl3T6kNcM6s4mnX5TYFCfOGcI/GnhwFIrX6W6KVfU5gQkZLV2ZPmkY6Dm3UtbI1xQmPxNOvymwKFDDvnHPt++nl6lv8HhKNEL/hrep5bQnjGfIUKESkppdSsy28KFDKs+oeJ2hseJDJzAeEZ83NzKBQqRKTY9W/WddGUKJdPLe5mXX5ToJBh0zdMhKbNI3bZ/yIycwGQnVPxyXtJrX9JgUJEitIhm3U1x6iNlPeKjcFQoJBh0TdMRC/8ErFF//SOccWIhjxEpAiVS7MuvylQyJA7mjAhIlJsBmzWNbuWM5tKs1mX3xQoZEgpTIhIqSnXZl1+U6CQIeOFiS8oTIhISdgWT/NQe5xnN3rNuk4bV8HClhjTy6RZl98UKGRIOOfo/u/Pk3h+icKEiBS1dbuSPNinWde5Eyq5sgybdflN716ZCAaM2mhhvp3OObruvYnE80uoueTL1P/F/80rTCTT+TdWrYyU11KupMZvc0zvheRhbyLDV5/voiJY/s26/KZAIXnZHyb2LrtjSMKEiMhw6k07Mg4+enw1F0yOFrqcsqKFtDJoChMiIrKfAoUMisKEiIj0pUAhx0xhQkRE+lOgkCPa/fj/o2flM8DBYSIy41yFCRERARQo5ChEmufy9pKr6V75dC5MWDhK3RX/W2FCREQArfKQo1A5awHHLf4Jb333Skj2YOEooz7XRnTW+YUuTUREioTOUMgROefofvnnkOwBoObimxUmRETkIAoUclh950xYOErt5X/L3mX/kZtTISIiAgoUchj9w8Soz7VRv+g2Ri2+j7eXXK1QISIiOQoUMqD+qzn6zpmonLWAUYvvI7HupQJXKSIixUKTMuUdjuY6E5WzFlA5a0GBKhQRkWKjMxRyEF20SkREBkOBQnIUJkREZLAUKARQmBARkfwoUIjChIiI5E2TMsuEAaHgsedD5xw7fqQwISL5MbPFwGKAyZMnF7iaQ4uGjEWtMabU6tffUNMZihHMCxM3svvZ26m79CsKEyIyaM65Jc65uc65uaNHjy50OYcUCwf4yKxqptWHC11K2VGgGKH6h4nGq76hMCEiIoOmQDECKUyIiMhQU6AoADO71MxWmdkaM7t1gMcrzOz+7OMvmlnzUL22woSIiAwHBQqfmVkQ+DfgMuAE4MNmdkK/3a4Dupxz04B/Af55sK+35dvvZeeT3wIODhOhUVMVJkREZMhomqv/TgfWOOc6AMzsPmAR8FqffRYB/5C9/TPge2Zmzjl3rC8WPeFCOu+/Becc6e3t7H72dgBqzr9RYUJERIaMAoX/JgAb+tzfCLznUPs451Jmtgs4Dnj7WF+s/uKbvetM/PSW3LaGv/wmDZd86VifSkRE5JA05OG/gU4L9D/zcDT7YGaLzewlM3vprbffOuQLVs9eBOZ9qyumn6MwISIiQ06Bwn8bgUl97k8ENh9qHzMLAXVAZ/8nOmjd96hDr/ve92obOEfF9HPofeP53JwKERGRoaJA4b/fAdPNbKqZRYCrgbZ++7QB12ZvXwU8PZj5EwA7n/wWnfffQuOHvsmEW5fT+KFv0nn/LQoVIiIypDSHwmfZORE3AU8AQeAu59yfzew24CXnXBtwJ/BfZrYG78zE1YN9ve7Xfknjh75J/cU3A+Q+dr/2y9xtERGRfClQFIBz7lHg0X7b/q7P7R7gg0PxWk1ffPQd2+ovvllhQkREhpSGPERERCRvChQiIiKSNwUKERERyZsChYiIiOTNBrkaUYqMmb0FrD/CbqMYxNU2h4HqOJjqONiR6pjinDv0hVek4A5zPCqWn7GjUUq1wvDWe1T/5xQoRhAze8k5N1d1qA7VIYVQSt/bUqoViqNeDXmIiIhI3hQoREREJG8KFCPLkkIXkKU6DqY6DlYsdcjQK6XvbSnVCkVQr+ZQiIiISN50hkJERETypkAxApjZpWa2yszWmNmtBaphkpk9Y2avm9mfzewLhaijTz1BM3vFzB4uYA31ZvYzM1uZfV/OLFAdf539nvzJzH5iZpU+vvZdZrbdzP7UZ1ujmT1lZm9kPzb4VY/k70jHGzOrMLP7s4+/aGbN/leZq+WIxyUzm29mu8zs99l/fzfQc/nFzNaZ2R+ztbw0wONmZv+afX//YGZz/KpNgaLMmVkQ+DfgMuAE4MNmdkIBSkkBX3LOHQ+cAdxYoDr2+wLwegFfH+A7wOPOuVnAKYWox8wmAJ8H5jrnTsLrgDvo7raD8EPg0n7bbgV+5ZybDvwqe19KwFEeb64Dupxz04B/Af7Z3yoPcrTHpeXOudnZf7f5W+KAFmRrGWiZ6GXA9Oy/xcDtfhWlQFH+TgfWOOc6nHMJ4D5gkd9FOOe2OOdWZG/vwfvlOcHvOgDMbCJwOfD9Qrx+toZa4Fy8VvU45xLOuZ0FKicERM0sBMSAzX69sHPu10Bnv82LgLuzt+8G3udXPZK3ozne9P3+/gy4wMzMxxpzium4NIQWAfc4zwtAvZk1+fHCChTlbwKwoc/9jRT4P0z2FOepwIsFKuHbwFeATIFeH6AFeAv4QXbo5ftmVuV3Ec65TcA3gTeBLcAu59yTftfRz1jn3BbwDvjAmALXI0fvaI43uX2ccylgF3CcL9UdxhGOS2ea2atm9piZnehrYe/kgCfN7GUzWzzA4wU75itQlL+Bkn/BlvaYWTXwc+CLzrndBXj9K4DtzrmX/X7tfkLAHOB259ypwD4KcGo/Oz9hETAVGA9Umdlf+V2HlI2jOd4U1TEJjnhcWoF36elTgO8CD/pdXz9nO+fm4A1t3Ghm5/Z7vGDvrwJF+dsITOpzfyI+ntLuy8zCeP9pf+yc+0UhagDOBhaa2Tq807Hnm9mPClDHRmCjc27/X0M/wwsYfrsQWOuce8s5lwR+AZxVgDr62rb/FG324/YC1yNH72iON7l9ssNsdbxz2Ms3RzouOed2O+f2Zm8/CoTNbJTPZfatZ3P243bgAbxhpr4KdsxXoCh/vwOmm9lUM4vgTbhr87uI7BjpncDrzrlv+f36+znnvuqcm+ica8Z7L552zvn+F7lzbiuwwcxmZjddALzmdx14Qx1nmFks+z26gMJPVm0Drs3evhZYWsBa5NgczfGm7/f3Krz/gwU5Q3E0xyUzG7d/joeZnY73e3OHf1UeVEuVmdXsvw1cDPyp325twMeyqz3OwBvG3OJHfSE/XkQKxzmXMrObgCfwZvDf5Zz7cwFKORv4KPBHM/t9dtvXsol/pPoc8OPsgbcD+ITfBTjnXjSzn+Gd1k0Br+DjFffM7CfAfGCUmW0E/h74BvBTM7sOL/B80K96JD+HOt6Y2W3AS865Nrxf4P9lZmvwzkz4uaqovwGPS8BkAOfcHXih5zNmlgK6gasLFYCAscAD2XwTAu51zj1uZjf0qfdR4L3AGiCOj8cVXSlTRERE8qYhDxEREcmbAoWIiIjkTYFCRERE8qZAISIiInlToBAREZG8KVCIiIhI3hQoREREJG8KFCIiIpK3/w+7+fp/3q6urAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\n", "total_cost, pointwise_cost, accumulated_cost= dtw(x, z, 'euclidean')\n", "path = get_path(accumulated_cost)\n", "_ = plot_warped_timeseries(x, z, pointwise_cost, \n", " accumulated_cost, path)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# TODO: fix plot dimensions for unequal length case" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python [conda env:clipseq]", "language": "python", "name": "conda-env-clipseq-py" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.2" } }, "nbformat": 4, "nbformat_minor": 2 }