Skip to content

Instantly share code, notes, and snippets.

@ischurov
Created July 7, 2020 23:53
Show Gist options
  • Save ischurov/3bfc8e317725c8331940111a650ec12c to your computer and use it in GitHub Desktop.
Save ischurov/3bfc8e317725c8331940111a650ec12c to your computer and use it in GitHub Desktop.
peaks.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "import pandas as pd\nimport matplotlib.pyplot as plt\nimport seaborn as sns\nimport numpy as np\n%matplotlib inline",
"execution_count": 2,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "def mkpct(df, candidates, allvotes, og=False):\n if og:\n allballots = ['Число бюллетеней, выданных участникам голосования']\n allvoters = 'Число участников голосования, включенных в список участников голосования на момент окончания голосования'\n \n else:\n allballots = ['Число избирательных бюллетеней, выданных избирателям, проголосовавшим досрочно',\n 'Число избирательных бюллетеней, выданных в помещении для голосования в день голосования',\n 'Число избирательных бюллетеней, выданных вне помещения для голосования в день голосования']\n allvoters = 'Число избирателей, включенных в список избирателей'\n\n\n df = df.copy()\n allvotes_values = df[list(allvotes)].sum(axis=1)\n for candidate in candidates:\n df[f\"{candidate}_pct\"] = df[candidate] / allvotes_values\n \n allballots_values = df[allballots].sum(axis=1)\n return df.assign(turnout=lambda x: (allballots_values / \n x[allvoters]))\nresults = (pd.read_csv(\"results.txt\", sep=\"\\t\")\n .rename(columns={\"ДА\":'yes', \"НЕТ\": 'no'})\n .pipe(lambda x: mkpct(x,\n candidates=['yes'],\n allvotes=['yes', 'no'],\n og=True))\n .dropna())\n\nresults_mc = results.assign(\n **{'Число бюллетеней, выданных участникам голосования':\n np.random.binomial(\n results['Число участников голосования, включенных в список участников голосования на момент окончания голосования'],\n results['turnout'])}).pipe(\n lambda x: mkpct(x, candidates=['yes'],\n allvotes=['yes', 'no'],\n og=True))",
"execution_count": 3,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "plt.figure(figsize=(10, 10))\nsns.distplot(results['turnout'], ax=plt.gca())\nfor v in [0.7, 0.75, 0.8, 0.85, 0.9, 0.95]:\n plt.plot([v, v], [0, 3], label=str(v))\nplt.legend()",
"execution_count": 12,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 12,
"data": {
"text/plain": "<matplotlib.legend.Legend at 0x7fa46801f3c8>"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 720x720 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlMAAAJNCAYAAAAVnfADAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeXSU12H38d+dGc2MNJLQjmAkFiGQgQDCgG28Am5MTFqc1ITgtE5ct3EW0iVN2mY5h6Z+08Ztkrqp/bZ53Tp1NkySOjbUtXGIY8cbBoORbRAYAUJoAW1o3zXzvH9IyMhskmZGzyzfzzkcazSPZn46FtbP997nXmNZlgAAADAxDrsDAAAAxDLKFAAAQAgoUwAAACGgTAEAAISAMgUAABACyhQAAEAIXHa9cU5OjjVr1iy73h4AAGDM9u/f32RZVu7FnrOtTM2aNUv79u2z6+0BAADGzBhTdannmOYDAAAIAWUKAAAgBJQpAACAENi2ZgoAAMSWgYEB1dTUqLe31+4oEeP1elVQUKCkpKQxfw1lCgAAjElNTY3S0tI0a9YsGWPsjhN2lmWpublZNTU1mj179pi/jmk+AAAwJr29vcrOzo7LIiVJxhhlZ2ePe+SNMgUAAMYsXovUORP5/ihTAAAgpuzcuVMlJSUqLi7WAw88cMHzX/ziF1VaWqrS0lLNmzdPGRkZEc3DmikAABAzAoGANm/erF27dqmgoEArVqzQ+vXrtWDBgpFrHnzwwZGPH3roIR04cCCimRiZAgAAMWPv3r0qLi5WUVGR3G63Nm3apO3bt1/y+scff1x33XVXRDNRpgAAQMyora1VYWHhyOOCggLV1tZe9NqqqipVVlZqzZo1Ec3ENB8AABi/Z78inXknvK+Zv0i6/cI1UBO1bds2bdiwQU6nM2yveTGMTAEAgJjh9/tVXV098rimpkZ+v/+i127bti3iU3wSI1MAAGAiwjiCNB4rVqxQRUWFKisr5ff7tW3bNm3duvWC644cOaKWlhatXLky4pkYmQIAADHD5XLp4Ycf1tq1azV//nxt3LhRCxcu1JYtW7Rjx46R67Zt26ZNmzZNyr5YxrKsiL/JxSxfvtzat2+fLe8NAADG7/Dhw5o/f77dMSLuYt+nMWa/ZVnLL3Y9I1MAAAAhoEwBAACEgDIFAAAQAsoUAABACK5YpowxXmPMXmPMW8aYQ8aYv7vINR5jzM+MMceMMXuMMbMiERYAACDajGVkqk/SGsuylkgqlfQhY8x177vmjyW1WJZVLOlBSf8Y3pgAAADR6YplyhrSOfwwafjP+/dTuEPSD4c//m9Jt5rJ2NgBAAAknJ07d6qkpETFxcV64IELNw/94he/qNLSUpWWlmrevHnKyMgYec7pdI48t379+rDkGdMO6MYYp6T9kool/V/Lsva87xK/pGpJsixr0BjTJilbUlNYUgIA7PfsV4b+GYGdr/9x79CExt9c8zdhf+0z//APkqT8r30t7K/98s+PSpJu2jgv7K/9wmOPSJJW33Nf2F+79X+OS5Iyfm9O2F870gKBgDZv3qxdu3apoKBAK1as0Pr167VgwYKRax588MGRjx966CEdOHBg5HFycrLKysrCmmlMZcqyrICkUmNMhqQnjTEfsCzr4HjfzBhzn6T7JGnGjBnj/XIAgJ3CfajteY6cPRKx1+47HLnXbqruvPJFE9RQdSJir91f1xWx1460vXv3qri4WEVFRZKkTZs2afv27aPK1Pkef/xx/d3fXbDcO6zGdTefZVmtkl6Q9KH3PVUrqVCSjDEuSVMkNV/k6x+xLGu5ZVnLc3NzJ5YYAAAkrNraWhUWFo48LigoUG1t7UWvraqqUmVlpdasWTPyud7eXi1fvlzXXXednnrqqbBkuuLIlDEmV9KAZVmtxphkSR/UhQvMd0j6lKTdkjZI+o1l1zk1AAAg4v5x7z+GfUTxqqyrwjrVu23bNm3YsEFOp3Pkc1VVVfL7/Tpx4oTWrFmjRYsWac6c0KY7xzIyNU3SC8aYtyW9IWmXZVlPG2PuN8acW7n1qKRsY8wxSX8p6SshpQIAALgIv9+v6urqkcc1NTXy+/0XvXbbtm266667Lvh6SSoqKtKqVatGraeaqCuOTFmW9bakpRf5/JbzPu6V9LGQ0wAAgJgQiZsFxmLFihWqqKhQZWWl/H6/tm3bpq1bt15w3ZEjR9TS0qKVK1eOfK6lpUUpKSnyeDxqamrSq6++qr/+678OOdOYFqADAABEA5fLpYcfflhr165VIBDQvffeq4ULF2rLli1avnz5yHYH27Zt06ZNm3T+Tk2HDx/WZz7zGTkcDgWDQX3lK1+55ML1cWUK+RUAAAAm0bp167Ru3bpRn7v//vtHPf7GN75xwdddf/31eued8N+Vytl8AAAAIaBMAQAAhIAyBQAAEALKFAAAQAgoUwAAACGgTAEAAISAMgUAAGLKzp07VVJSouLiYj3wwAMXPH/q1CmtXr1aS5cu1eLFi/XMM89ENA9lCgAAxIxAIKDNmzfr2WefVXl5uR5//HGVl5ePuuab3/ymNm7cqAMHDmjbtm36/Oc/H9FMlCkAABAz9u7dq+LiYhUVFcntdmvTpk3avn37qGuMMWpvb5cktbW1afr06RHNxA7oAAAgZtTW1qqwsHDkcUFBgfbs2TPqmm984xu67bbb9NBDD6mrq0u//vWvI5qJMgUAAMbtzD/8g/oOHwnra3rmX6X8r30t5Nd5/PHHdc899+hLX/qSdu/erbvvvlsHDx6UwxGZCTmm+QAAQMzw+/2qrq4eeVxTUyO/3z/qmkcffVQbN26UJK1cuVK9vb1qamqKWCZGpgAAwLiFYwRpIlasWKGKigpVVlbK7/dr27Zt2rp166hrZsyYoeeff1733HOPDh8+rN7eXuXm5kYsEyNTAAAgZrhcLj388MNau3at5s+fr40bN2rhwoXasmWLduzYIUn67ne/q//4j//QkiVLdNddd+mxxx6TMSZymSL2ygAAABGwbt06rVu3btTn7r///pGPFyxYoFdffXXS8jAyBQAAEALKFAAAQAgoUwAAACGgTAEAAISAMgUAABACyhQAAEAIKFMAACCm7Ny5UyUlJSouLtYDDzxwwfOnTp3S6tWrtXTpUi1evFjPPPOMJOnkyZNKTk5WaWmpSktL9dnPfjYsedhnCgAAxIxAIKDNmzdr165dKigo0IoVK7R+/XotWLBg5JpvfvOb2rhxoz73uc+pvLxc69at08mTJyVJc+bMUVlZWVgzMTIFAABixt69e1VcXKyioiK53W5t2rRJ27dvH3WNMUbt7e2SpLa2Nk2fPj2imShTAAAgZtTW1qqwsHDkcUFBgWpra0dd841vfEM/+clPVFBQoHXr1umhhx4aea6yslJLly7VLbfcopdffjksmZjmAwAA4/byz4+qqbozrK+ZU5iqmzbOC/l1Hn/8cd1zzz360pe+pN27d+vuu+/WwYMHNW3aNJ06dUrZ2dnav3+/PvKRj+jQoUNKT08P6f0YmQIAADHD7/erurp65HFNTY38fv+oax599FFt3LhRkrRy5Ur19vaqqalJHo9H2dnZkqRly5Zpzpw5Onr0aMiZGJkCAADjFo4RpIlYsWKFKioqVFlZKb/fr23btmnr1q2jrpkxY4aef/553XPPPTp8+LB6e3uVm5urxsZGZWVlyel06sSJE6qoqFBRUVHImShTAAAgZrhcLj388MNau3atAoGA7r33Xi1cuFBbtmzR8uXLtX79en33u9/Vpz/9aT344IMyxuixxx6TMUYvvfSStmzZoqSkJDkcDn3/+99XVlZW6JnC8H0BAABMmnXr1mndunWjPnf//fePfLxgwQK9+uqrF3zdnXfeqTvvvDPseVgzBQAAEALKFAAAQAgoUwAAACGgTAEAAISAMgUAABACyhQAAEAIKFMAACCm7Ny5UyUlJSouLtYDDzxwwfNVVVW69dZbtXjxYq1atUo1NTURzUOZAgAAMSMQCGjz5s169tlnVV5erscff1zl5eWjrvnyl7+sT37yk3r77be1ZcsWffWrX41oJsoUAACIGXv37lVxcbGKiorkdru1adMmbd++fdQ15eXlWrNmjSRp9erVFzwfbpQpAAAQM2pra1VYWDjyuKCgQLW1taOuWbJkiX75y19Kkp588kl1dHSoubk5Ypk4TgYAAIzbC489ooaqE2F9zbyZRVp9z30hv853vvMdfeELX9Bjjz2mm2++WX6/X06nMwwJL44yBQAAYobf71d1dfXI45qaGvn9/lHXTJ8+fWRkqrOzU0888YQyMjIilokyBQAAxi0cI0gTsWLFClVUVKiyslJ+v1/btm3T1q1bR13T1NSkrKwsORwOfetb39K9994b0UysmQIAADHD5XLp4Ycf1tq1azV//nxt3LhRCxcu1JYtW7Rjxw5J0osvvqiSkhLNmzdP9fX1+vrXvx7ZTBF9dQAAgDBbt26d1q1bN+pz999//8jHGzZs0IYNGyYtDyNTAAAAIaBMAQAAhIAyBQAAEALKFAAAQAgoUwAAACGgTAEAAISAMgUAAGLKzp07VVJSouLiYj3wwAMXPF9VVaVbb71Vixcv1qpVq1RTUzPynNPpVGlpqUpLS7V+/fqw5GGfKQAAEDMCgYA2b96sXbt2qaCgQCtWrND69eu1YMGCkWu+/OUv65Of/KQ+9alP6Te/+Y2++tWv6sc//rEkKTk5WWVlZWHNxMgUAACIGXv37lVxcbGKiorkdru1adMmbd++fdQ15eXlWrNmjSRp9erVFzwfbpQpAAAQM2pra1VYWDjyuKCgQLW1taOuWbJkychBx08++aQ6OjrU3NwsSert7dXy5ct13XXX6amnngpLJqb5AADAuLX+z3H113WF9TXd033K+L05Ib/Od77zHX3hC1/QY489pptvvll+v19Op1PS0Hoqv9+vEydOaM2aNVq0aJHmzAntPSlTAAAgZvj9flVXV488rqmpkd/vH3XN9OnTR0amOjs79cQTTygjI2Pk6yWpqKhIq1at0oEDByhTAABg8oVjBGkiVqxYoYqKClVWVsrv92vbtm3aunXrqGuampqUlZUlh8Ohb33rW7r33nslSS0tLUpJSZHH41FTU5NeffVV/fVf/3XImVgzBQAAYobL5dLDDz+stWvXav78+dq4caMWLlyoLVu2aMeOHZKkF198USUlJZo3b57q6+v19a9/XZJ0+PBhLV++XEuWLNHq1av1la98ZdRdgBPOFPIrAAAATKJ169Zp3bp1oz53//33j3y8YcMGbdiw4YKvu/766/XOO++EPQ8jUwAAACGgTAEAAISAaT4AAMZp655TkqSOjt5Rj8/3iWtnTGom2IeRKQAAMGaWZdkdIaIm8v1RpgAAwJh4vV41NzfHbaGyLEvNzc3yer3j+jqm+QAAwJgUFBSopqZGjY2NdkeJGK/Xq4KCgnF9DWUKAACMSVJSkmbPnm13jKjDNB8AAEAIKFMAAAAhoEwBAACEgDIFAAAQAsoUAABACChTAAAAIaBMAQAAhIAyBQAAEIIrliljTKEx5gVjTLkx5pAx5s8vcs0qY0ybMaZs+M+WyMQFAACILmPZAX1Q0pcsy3rTGJMmab8xZpdlWeXvu+5ly7J+N/wRAQAAotcVR6YsyzptWdabwx93SDosyR/pYAAAALFgXGumjDGzJC2VtOciT680xrxljHnWGLMwDNkAAACi3pgPOjbGpEp6QtJfWJbV/r6n35Q007KsTmPMOklPSZp7kde4T9J9kjRjxowJhwYAAIgWYxqZMsYkaahI/dSyrF++/3nLstoty+oc/vgZSUnGmJyLXPeIZVnLLctanpubG2J0AAAA+43lbj4j6VFJhy3L+udLXJM/fJ2MMdcMv25zOIMCAABEo7FM890g6W5J7xhjyoY/9zVJMyTJsqzvS9og6XPGmEFJPZI2WZZlRSAvAABAVLlimbIs6xVJ5grXPCzp4XCFAgAAQ7buOXXFaz5xLeuQ7cQO6AAAACGgTAEAAISAMgUAABACyhQAAEAIKFMAAAAhoEwBAACEgDIFAAAQgjGfzQcAQKxjzyZEAiNTAAAAIaBMAQAAhIAyBQAAEALKFAAAQAgoUwAAACGgTAEAAISAMgUAABACyhQAAEAIKFMAAAAhoEwBAACEgDIFAAAQAsoUAABACDjoGACACWjq6FNP94C6+gJ6/oVjKslPU8nUNPkzk+Uwxu54mESUKQAAxigQtPST16v07y8e15n2Xn28yy1vklNOh9ELRxr0myMN8nlcumZWpj62vEBJTiaAEgFlCgCAMahs6tJf/eIt7atqUWFmsj68aJpmHOmRy2n02VsK1N03qKMNHTpY264X3m3Ux76/Ww/dtVSFWSl2R0eEUaYAALiMYNDSf712Ut9+7ojcToe++7El6h0IyBijjoqaketSPC6VFmaqtDBT79S26em367Tuey/rgTsX68OLp9n4HSDSGH8EAOASTjZ16eOP7Nb/ebpc18/J0a6/vEV3LiuQucKaqEX+KXrmz27SnLxUbd76pv7y52Vq6Ogd03vWtvbopaONOtvZr8aOPu0+0awzbWP7WtiDkSkAAN4nGLT02Gsn9U/PHVGS06HvfGyJ7rzaf8USdb7CrBT94rMr9S+/PqpHXjqhXx2q1xfWFOuPbpglj8s56tpA0NJvjzbox7ur9OLRRlmW9NGGTknS/7xVpySn0adWzlJRbmpYv0+EB2UKAIDzVDV36a/++23trTyr1SW5+tbvL1b+FO+EXivJ6dBfrb1KG5YV6ptPl+uBZ49o295T+vDiaRoMWhoYtNQ3GNBvjzaqpqVHeWke/emaubppbo4OPvKinA6j2avm6Sd7qvTD3Sd193WzVJxHoYo2lCkAAIYdb+zUA88eliXp2xsWa8MYpvTGYnaOT4/es0K/Pdqov//fcv3bi8fldjrkdjrkchpdlZ+ur94+X7ctnDpyB+AJz9DoVU6aR39yU5F+8EqlfrT7pP7wupmaNzUt5EwIH8oUAACS9p08q6fKajUnN1U/uGdFRO7Cu2Verm6ee/O4C1qqx6U/vnG2fvBqpX78epXuplBFFRagAwASWtCy9NyhM/rlgVoV5abqic9fH9HtDCY60uUbLlQ5qW49daBWgaAV5mSYKMoUACChPXWgVr892qgVs7L0qZWzlO5NsjvSJaW4XfrQwmlq7RlQWXWL3XEwjDIFAEhYh0+3a19Vi26am6OPlE6X0xH9x8DMm5qq6Rlevfhuo4IWo1PRgDIFAEhIPf0BPVVWq/x0rz64YGpYFppPBmOMVs3LU3NXv96pabM7DkSZAgAkqGcOnlZX36DuvLpALkds/TpcMD1deWkevfBuA6NTUSC2fnoAAAiDo/Ud2l/Vopvm5sqfmWx3nHFzGKNVJXlq6OjTkdPtdsdJeJQpAEBC6R0I6MkDtcpN9WjNVXl2x5mwRf4pyva59cK7jbIYnbIVZQoAkFB2ldervWdAd17tH9kgMxY5HUa3zMtVbWuPfnu00e44CS12f4oAABinnv6A9lWd1dUzMjUj22d3nJCVzsjQlOQk/efLlXZHSWiUKQBAwiirbtFAwNJ1Rdl2RwkLl8Oh5TMz9erxJtW19tgdJ2FRpgAACcGyLL1xskXTM7wxuej8UpbOyJRlSU8eqLU7SsLibD4AQEI4UN2qM+29+kip3+4oI7buOXXJ55rb+yRJVxpDy/K5de3sLD2xv0afXzUnZvbLiieMTAEAEsLWPafkdjm0pGCK3VHC7s5lBTrR1KU3T7XaHSUhMTIFAIh7g0FLT79dpyUFGfIkOSflPS836hRu6xZN099uP6T/3l+jZTMzJ+19MYSRKQBA3Gvq7FPvQFDXzMqyO0pEpHpcuv0D+Xr67Tr1DgTsjpNwKFMAgLhmSWpo79Mi/5S4Wnj+fhuWFaijd1C/Kq+3O0rCoUwBAOJaZ9+guvsHddc1M+yOElHXFWXLn5GsJ/bX2B0l4VCmAABxraG9V06H0frS6XZHiSiHw+ijS/16uaJR9e29dsdJKJQpAEDcGggEdbarX1k+t1I98X/P1Z3LChS0pF++yZ5Tk4kyBQCIW3tOnFUgaCnT57Y7yqSYnePT1TMytL2MMjWZ4r+mAwAS1nOHzmixwygjOcnuKJPmI0v92rL9kI6caddV+ekjnz9/q4bFHUPTgL9+3/YNn7g2vteVRQojUwCAuBQMWtpVXq+M5CQ5EmhX8HWLpsnpMNpRVmd3lIRBmQIAxKV3att0pr03Yab4zslJ9ejG4hxtL6uTZVl2x0kIlCkAQFx67tAZOR1GmSmJVaYk6Y7S6apt7dGbp1rsjpIQKFMAgLj0q/J6XVeUJZcjcab4zrltYb48Loe2M9U3KShTAIC4c7yxU8caOnXbgny7o9gi1ePS7yyYqv99+7QGAkG748Q97uYDAMSdXx0aOlLlgwumamCcXzuZBxRH0h1Lput/3z6tV481aVVJnt1x4hojUwCAuPPcoTNaXDBF0zPi9yy+K7mlJFfpXhd39U0CyhQAIK7Ut/eqrLpVty2YancUW3lcTq1bNE3PHTqjnv6A3XHiGmUKABBXdpUPTfGtXZiY66XOt750urr6A/rNkQa7o8Q11kwBAOLCubVOP3ztpLJ9bu2tPKs3TraouKNPkvRqnKyFGo9rZ2drarpH28tqWTcVQYxMAQDiRk9/QMcbO7VwerpMAu16filOh9HtH5imF482qm+Aqb5IoUwBAOLGu/UdClrSgmnpV744QXx48TT1DwZ15EyH3VHiFmUKABA3yk+3K83jUkFWit1RosayGZmamu7RO7VtdkeJW5QpAEBcGAgEdbS+Q/OnpSfUwcZX4hie6jta38FUX4RQpgAAceF4Y6f6B4NaMJ0pvvdbt2iaBoMWU30RQpkCAMSF8rp2eVwOFeX67I4SdZbPzFSa18VUX4RQpgAAMS8QtHT4dLtK8tPkcvCr7f0cDqOF06cMTfUNMtUXbvzEAQBi3punWtTVH+AuvstY5J/CVF+EUKYAADHvuYNn5HQYzZuaZneUqDUzO0VpHpcOMtUXduyADgCIaZZl6Vfl9SrOTZU3yWl3HFtsHcPu7g5jtNCfrn0nWxTMdHPHYxgxMgUAiGnv1nfo1NlupvjG4APDU30cfBxejEwBAGy1dc8pNbT3jXx8MZ+4dsYlv/65g/UyRrpqGlN8VzIr26c0j0tdfQH5PBOrAGMZBbvcv694xMgUACBmWZalp9+uG771P8nuOFHPYYbWlfUMBGRZdqeJH5QpAEDMequmTRUNnfr9qwvsjhIz5k5NVdCy1M8WCWFDmQIAxKxf7KuWN8mh3108ze4oMaM4L1WS1M26qbChTAEAYlLvQEA73qrTug9MY4pvHFLcLnlcDvVwTl/YUKYAADHpuUNn1NE7qA3LmeIbr+Qkp/oGg+ruG7Q7SlygTAEAYtLP91WrMCtZ183OtjtKzElxD93Jd6yx0+Yk8eGKZcoYU2iMecEYU26MOWSM+fOLXGOMMf9qjDlmjHnbGHN1ZOICACBVn+3Wa8ebteHqQjkcbD45Xh6XQw5jVFFPmQqHsWwyMSjpS5ZlvWmMSZO03xizy7Ks8vOuuV3S3OE/10r69+F/AgAQdk+8WSNJunOZ3+YkMcoMTfVVNHTIsiwZdkMPyRVHpizLOm1Z1pvDH3dIOizp/T+9d0j6kTXkdUkZxhhurQAAhF0waOm/99fohjk5KshMsTtOzEp2O9XeO6j64Q1TMXHjWjNljJklaamkPe97yi+p+rzHNbqwcAEAELLXK5tV09Kjj7HwPCTJw+cYHq3vsDlJ7BtzmTLGpEp6QtJfWJbVPpE3M8bcZ4zZZ4zZ19jYOJGXAAAkuJ+9Ua00r0trF+bbHSWmuZxGU9M9qmigTIVqTGXKGJOkoSL1U8uyfnmRS2olFZ73uGD4c6NYlvWIZVnLLctanpubO5G8AIAEVn22W0+/fVofW1Yo7/DICiZubl6aTjZ3q38waHeUmDaWu/mMpEclHbYs658vcdkOSZ8cvqvvOkltlmWdDmNOAAD0/146LoeRPn3zbLujxIV5U9MUCFo60cRdfaEYy918N0i6W9I7xpiy4c99TdIMSbIs6/uSnpG0TtIxSd2S/ij8UQEAiayhvVc/31ejDcsKNG1Kst1x4sLM7BQlOYe2SLgqP93uODHrimXKsqxXJF32nknLsixJm8MVCgCA9/vPVyo1GAjqMzfPsTtK3EhyOlSYlaKqs112R4lp7IAOAIh6rd39+snrVfq9JdM1K8dnd5y4UpiZojNtvRoIsG5qoihTAICo91+vnlR3f0CfX1Vsd5S4U5iZoqAl1bX22B0lZlGmAABRrW8goMdeO6kPLpiqkvw0u+PEncKsofVn1We7bU4SuyhTAICotqfyrNp6BrR5NaNSkZDmTVJGcpKqWxiZmijKFAAgag0EgnrlWJNuLM5RaWGG3XHiVmFWiqpbGJmaKMoUACBq7a9qUWffoD6/mjv4IqkwM1mt3QNq7OCcvomgTAEAolIgaOmlikbNyErRyqJsu+PEtcKsoQOjy6pbbU4SmyhTAICo9FZ1q1q7B7RqXq6GDuNApEzPSJbDSGXVLXZHiUmUKQBA1Alall482qj8dC938E2CJKdD+VO8jExNEGUKABB1DtW1q6mzT6tKGJWaLIWZKXq7uk3BoGV3lJhDmQIARBXLsvTbdxuU7XPrA/4pdsdJGIWZKeroG9TxRg49Hi/KFAAgqlQ0dKqurVe3zMuVg1GpSVMwvHnnAab6xo0yBQCIKr892qgpyUkqncG+UpMpJ9WjNK+LdVMTQJkCAESNM+29qmzq0sqibLkc/IqaTA5jVFqYobJTlKnx4icVABA1Xj/RLJfDaPnMTLujJKTSwgy9W9+hnv6A3VFiisvuAAAASEPbIZSdatXiggyleEb/etq655RNqRJLaWGGAkFL79S26ZrZWXbHiRmMTAEAokJH36D6A0F2O7fRufMP2bxzfChTAADbWbLU3jOgwsxk+TOT7Y6TsLJTPSrITNbbNW12R4kplCkAgO16BgIaCAR1HaNStrsqP03HGthrajwoUwAA27X3DMrpMFrEJszCFWIAACAASURBVJ22m5OXqhONXRoMBO2OEjMoUwAAW7V09aunf1Bp3iS5nPxastvcvDT1B4I6dbbb7igxg59aAICt9lQ2S0ZK93KDeTSYm5cqaWgneowNZQoAYJvBQFBvnmpVstvFJp1RYs5wmWLd1NjxkwsAsM1LFY3q7BtUmodRqWiR6nFp+hQvZWocKFMAANs8sb9WKW6nUtxOu6PgPMVT01TR0GF3jJhBmQIA2KKte0C7yuu1pDBDRsbuODjP3LxUHWvoVDBo2R0lJlCmAAC22PF2nfoDQV09g3P4os3cvFT1DgRV29pjd5SYwCQ1AMAWT+yvUcnUNE2f4pXYcDsqnDsDsaq5S5L0ny+fUEl+up2RYgIjUwCASXe8sVNl1a26c5lfxjDFF21y0zySpIaOPpuTxAbKFABg0j2xv0ZOh9FHSv12R8FFpLhdSvO41NBOmRoLyhQAYFIFgpaePFCrm+fmKC/da3ccXEJuukcNHb12x4gJlCkAwKR67XiTTrf16s5lBXZHwWXkpXnU0NEny+KOviuhTAEAJtUv36xVutel35k/1e4ouIy8NK/6BoNq7x20O0rUo0wBACZNd/+gnjt0Rh9ePE3eJDbqjGZ5I4vQmeq7EsoUAGDS7CqvV3d/QHew8DzqnVvPxiL0K6NMAQAmzfayOk2b4tU1s7LsjoIr8LmdSk5yqpHtEa6IMgUAmBRnu/r10tFGrV8yXQ4He0tFO2OM8rijb0woUwCASfG/75zWYNBiii+G5KV5Vd/OHX1XQpkCAEyK7QdqNW9qquZPS7M7CsYoL82jnoGAuvoDdkeJapQpAEDEVZ/t1r6qFt1RyvExsYQ7+saGMgUAiLgdb9VJktYvmW5zEowHd/SNDWUKABBRlmVpe1mtls3MVGFWit1xMA7pXpc8LgcjU1dAmQIARNSRMx06Wt+pj5QyKhVrjDHKSfXobFe/3VGiGmUKABBRT5XVyuUw+vBiylQsykhJUkvXgN0xohplCgAQMcGgpf8pq9NNc3OU5XPbHQcTkJniVkt3P9sjXIbL7gAAgNizdc+pMV03J9enurZe/c3tV0U4ESIlMyVJg0FLnX2DSvMm2R0nKjEyBQCImKfK6pSc5NTvzJ9qdxRMUGbK0IhiazdTfZdCmQIARMRgMKhn3jmt2xZOlc/DREisyhienm3pZhH6pVCmAAARUVHfqbaeAX2E42NiWmbK0NReC3f0XRJlCgAQEWXVrcryuXXj3By7oyAEHpdTKW6nWpjmuyTKFAAg7PoGAjp8ul0fXjRNSU5+1cS6c3f04eL4CQcAhN2h0+0aDFq6g40640JmShIjU5dBmQIAhN1b1a3KTEnSspmZdkdBGGSmuNXKXlOXRJkCAIRVR++AjjV0aklBhowxdsdBGGT43CN7TeFClCkAQFi9U9smS9KSwgy7oyBMRu7oY6rvoihTAICweqemTfnpXk1N99odBWFybuNOtke4OMoUACBsOnoHdOpstxZOT7c7CsIoY2RkijJ1MZQpAEDYHDnTIUvSAspUXPG4nPKx19QlUaYAAGFz+HS7MlOSlM8UX9zJ9A3d0YcLUaYAAGHRNxjQsYZOzZ+Wzl18cSiDjTsviTIFAAiLivpODQYtLZjGFF88ykxJUmv3gILsNXUByhQAICzKT7crOcmpmdk+u6MgAjJThvea6mWvqfejTAEAQhYIWnr3TIfmT0uT08EUXzzK5I6+S6JMAQBCdrK5Sz0DAc1nii9uZZzba4o7+i7gsjsAACD2lde1y+UwmpuXNurzW/ecsikRwu3cxp3c0XchRqYAACGxLEuHT7drbl6q3C5+rcQrt8shn8fFNN9F8FMPAAjJ6bZetfYMsFFnAshMSWKa7yIoUwCAkJSfbpeRVJJPmYp3mSluzue7CMoUACAkR+s7VJiVolQPy3DjXWZKklp72Gvq/ShTAIAJ6+kPqLalR8V5qXZHwSTISHErELTUwV5To1CmAAATVtnUKUvSnFzKVCLgjr6Lo0wBACbsWGOXkpxGhVnJdkfBJGDjzoujTAEAJux4Y6dmZfvkcvDrJBFk+ti482L46QcATEh7z4AaO/qY4ksgSU6HUj0u7uh7H8oUAGBCjjd2SpLmsPg8oaQnu9Tey8jU+ShTAIAJOd7YpeQkp6ZN8dodBZMo3ZvE3XzvQ5kCAIybZVk63tipolyfHMbYHQeTKM3rUjtlapQrliljzA+MMQ3GmIOXeH6VMabNGFM2/GdL+GMCAKJJc1e/2noGWC+VgNK8SerqG1QgyMad54xlZOoxSR+6wjUvW5ZVOvzn/tBjAQCi2bn1UsWUqYST5h3a6b6zj9Gpc65YpizLeknS2UnIAgCIEccbOjUlOUnZqW67o2CSpXuH9ppq72ER+jnhWjO10hjzljHmWWPMwjC9JgAgCgWDlk40dWlOrk+G9VIJ59zIFIvQ3xOOUynflDTTsqxOY8w6SU9JmnuxC40x90m6T5JmzJgRhrcGAEy28tPt6u4PsF4qQY2MTLE9woiQR6Ysy2q3LKtz+ONnJCUZY3Iuce0jlmUttyxreW5ubqhvDQCwwWvHmyRJRZSphOTzuGTEyNT5Qi5Txph8MzzOa4y5Zvg1m0N9XQBAdNp9vFk5qR5NSU6yOwps4HQYpXpc6mBkasQVp/mMMY9LWiUpxxhTI+lvJSVJkmVZ35e0QdLnjDGDknokbbIsi/slASAOWbK0v6pFJflpdkeBjdK8LkamznPFMmVZ1l1XeP5hSQ+HLREAIGr19AfU3juomVk+u6PARmneJEamzsMO6ACAMTs3GjEzO8XmJLATu6CPRpkCAIxZR++ActM8yvKxv1QiS09mF/TzUaYAAGPW3juoFbMy2V8qwaV5XbLELujnUKYAAGPSFwioPxDU8plZdkeBzc7tNcW6qSGUKQDAmJxbL7ViFmUq0bEL+miUKQDAmHT0DsphjOZPY1uERJfGLuijUKYAAGPS0TugVI9LLie/OhJdKrugj8LfCADAFXX0Dqi7PzAyvYPE5nQY+dgFfQRlCgBwRQdOtUp6b+ExkOZ1qb2HkSmJMgUAGIN9J89KklIZmcKwdG+SOvoYmZIoUwCAMXjjZIt8bqec7C+FYWlelzoYmZJEmQIAXMFAIKiy6taRO7gAaeiOvk52QZdEmQIAXEF5Xbt6Blh8jtHO7YLexS7olCkAwOW9MbxeipEpnO+9XdApU5QpAMBl7TvZosKsZLnZXwrneW8XdBahM2YLABhl655TIx9blqWXKxpVkp+m+o5eSdLz5z2PxJWefG4XdEam+N8MAMAlNXb2qas/oFnZPrujIMq8tws6I1OUKQDAJZ1s6pYkzcqhTGE0p8MoxeNiZEqUKQDAZZxs7lKax6Vsn9vuKIhC6V6OlJEoUwCAy6hs6tKsHJ8Mm3XiItK8Lu7mE2UKAHAJLd39ausZ0KzsFLujIEqleZMYmRJlCgBwCSebuiSxXgqXlj48MhW0EnsXdMoUAOCiKpu65E1yaGq61+4oiFJp3iR2QRdlCgBwCSebuzQr2ycH66VwCenDG3cm+h19lCkAwAU6egfU1NnP/lK4rLSRI2USe90UZQoAcIGTzUP7S81mvRQuY+RImR5GpgAAGOVkU5eSnEbTM5LtjoIolnpumq+PkSkAAEY52dylGVkpcjpYL4VLczkc8rmdjEzZHQAAEF16+gM609bLlggYk1SvS53czQcAwHuqznbJkjSbxecYA5/HxdYIdgcAAESXk01dchqjwix2PseV+dwudfVTpgAAGFHR0KnCrGQlOfkVgSvzeZjm428KAGBEXWuPTrf16qr8dLujIEakepzqHQhqMBi0O4ptKFMAgBHPH2mQJF2Vn2ZzEsQKn2doe4TuvoDNSexDmQIAjHj+cL2yfG7lpnnsjoIY4XMPlalEXjdFmQIASJK6+wf12vFmXZWfJsN5fBijcyNTibxuijIFAJAkvVLRpP7BIOulMC6pw2Wqi2k+AECie/5wg9I8Ls3KYUsEjJ3P45SkhN5rijIFAFAwaOn5Iw26uSRXLge/GjB23iSnHIYyBQBIcG/Xtqmps0+3XpVndxTEGIcx8rkTe68pyhQAQM8frpfDSKtLKFMYP5/Hpa5+1kwBABLY84cbtGxmpjJ9brujIAb5PE6m+QAAiauutUflp9t16/ypdkdBjEr0w44pUwCQ4M7tev4785niw8Qk+vl8lCkASHA7ymo1KztFc3JT7Y6CGJXqcalvMKjBQGKez0eZAoAEduBUi9442aK7V85i13NM2HtHyiTmInTKFAAksP98uVJpXpc+vqLQ7iiIYYm+cSdlCgAS1Knmbj178LT+4NqZI0eCABORmuDn81GmACBB/eDVSjkdRn90wyy7oyDGjUzzUaYAAImipatfP3ujWneU+jU13Wt3HMQ4n4cyBQBIMD/dU6WegYA+fVOR3VEQB7xJDjmNYQE6ACAx9A4E9NhrVbplXq5K8tPsjoM4YIyRz+NkzRQAIDE8daBWTZ19uu9mRqUQPom8CzplCgASyMmmLv3DM4e1pDBD18/JtjsO4ghlCgAQ9zr7BvXpH+2Tw2H08F1L2aQTYeVzOxN2zRQbiwBAAggGLX3xZ2U60dSlH917jQqzUuyOhDiTmsDn8zEyBQAJ4HvPV2hXeb2+vm6+bijOsTsO4pDP41L/YFADCXg+H2UKAOLcUwdq9b3nK7RhWQEbdCJiEnmvKab5ACBO9fQH9M3/LddP95zSilmZ+uZHPsA6KUTMe7ugJ966KcoUAMSArXtOXfGaT1w7Y+Tj8rp2/dm2AzrW0KnP3FykL91WIreLyQhETuq5w477GZkCAMSwQNDSD16p1Lefe1cZKUn6yR9fqxvnskYKkedL4MOOKVMAECf+ZddRPfFmjapbejR/Wrp+f6lfp852jxrVOn/0Cggn1kwBAGJWIGjplYpGPX+kQUlOhzYuL9SSgimsj8Kk8rgccjoMZQoAEFt6+gP6yZ4qVTZ1aeH0dK1fMl1p3iS7YyEBGWOGNu5kAToAIFa0dvfrsddOqrmzXxuuLtDVMzPtjoQEl6gbd1KmACAGnWnr1WOvVapvMKhPXT9LxXmpdkcChs7n424+AEC0q2zq0o92n5TH5dB9Nxdp2pRkuyMBkobKVFNnn90xJh1lCsAljXdvI0ReQ0evfvz6SaV7k/RHN8xSRop7XF8/ln+nwEQl6popdnADgBjR1TeoH+2uktPh0D0TKFJApKV6XOoPBNXTn1iFipEpACFh9GpyDAaD+umeU2rvGdCf3FSkTIoUotC5vaaau/pU4E6xOc3koUwBGJOO3gHVtPSopqVbNS096uwbVH66V9MzkjUtwyv/lGR5kpx2x4xLlmVp+4E6nWzu0sblhZqRlTi/pBBbRspUZ78KMhPn55QyBeCyOvsGtaOsVgfr2iVJDiNNTfcqzevS8cZOHahulSQlOY2un5Ojm+fmKtlNqQqnR1+p1P5TLVpdkqfSwgy74wCXdK5Mne3qtznJ5KJMAQnqStNzlmXpndo27XirTn2DQa0uydW8qWmaNiV51IG5Hb0DqmvtVVl1i1462qg9lc26eW6urp+Tw8G6YXCork3/uPOIFkxL163z8+yOA1xW6nCZSrQ7+ihTAC7Q3TeoJ8tqdaiuXQWZybrz6gJNTfde9No0b5JK8pNUkp+mm+f1aFd5vX5VXq/XTzRr44pCFeWw/9FE9Q4E9MWflSkjxa2PLvXLwfEwiHK+4VHpRBuZ4n8bAYwyEAjqh7tP6siZDq1dMFWfuXnOJYvU+02bkqxPrpylz9xcJLfLoUdfrtQL7zYoGLQiGzpO/dPOd3W0vlPf3rB4ZPoEiGZul0Muh1EzZQpAorIsS08eqFV1S48+vrxQt5TkyekY/2jIzGyfNq8q1qKCKdpVXq9P/dfehBv2D9UrFU36wauV+uTKmVpVwvQeYoMxRj6PS82dlCkACerFo40qq27VbQum6gP+KSG9lifJqY8vL9RHS/3aW3lWH/7Xl/V2TWuYksa31u5+ffkXb2lOrk9fvX2+3XGAcfF5nGrppkyNYoz5gTGmwRhz8BLPG2PMvxpjjhlj3jbGXB3+mAAi7WBtm3aV16u0MEO3zMsNy2saY7Ridpae/PwNSnI69LHv79bTb9eF5bXj2Td2HFJTZ5/+5eNLuTMSMcfndrFm6iIek/Shyzx/u6S5w3/uk/TvoccCMJnqWnv0i/3VKsxM1keX+mXCvNB5wfR0PbX5Bi3yT9EXth7Qg7uOso7qEp4/XK+nyuq0efXQNCkQa1LczoQrU1dc0WhZ1kvGmFmXueQOST+yLMuS9LoxJsMYM82yrNNhyggggizL0o636uRxOfWH181UkjMys/85qR799NPX6utPHtT3nq9QRUOHvr1hCQur9d42FT39AX3v+aPKT/cqO9XNOXqISSkel040dtkdY1KF479ifknV5z2uGf4cZQqIAQfr2nXqbLc+utSvNG9SRN7j/FKwtDBDXX2DevadM3rjZIv+4JoZykv3cuSMpGcPnlZH76D+8LqZcjlY0orY5HM71dE3qP7BYMLsNTep36Ux5j5jzD5jzL7GxsbJfGsAFzEYDOq5Q2eUl+bRspmZk/KexhjdNDdX9944W919g/q3F4+zMF3SsYZO7atq0U1zcxLqGA7En3Ojza0JtAg9HGWqVlLheY8Lhj93AcuyHrEsa7llWctzc8OzwBXAxO2tPKuzXf26/QPTJn1DyDm5qfrCmrnKn+LVtjeq9bfbD6qrb3BSM0SLvsGAnjxQo2yfW7fOn2p3HCAkKe7hI2UoU+OyQ9Inh+/qu05SG+ulgOjX0x/Q84cbNCfXp3lT7dmlfEpykv7kptm6fk62fri7Sh/8599q58EzGlqCmTh+daheLd0D+v2rCyK2Zg2YLIm4C/pYtkZ4XNJuSSXGmBpjzB8bYz5rjPns8CXPSDoh6Zik/5D0+YilBRA2vz3aoN6BgG7/wLSw3703Hi6HQ7+7eLqe+NxKpScn6bM/2a97H3tDVc2JsYD11WNN2n2iWSvnZGt2js/uOEDIUhLwsOOx3M131xWetyRtDlsiABFX09Kt1443q7QwQ9Mzku2OI0laNjNLT//pjfrh7ir986/e1Zrv/la/Mz9Pf3DtTN1YnCPHBHZij3ZtPQP68i/eUk6qR2sX5NsdBwiLcyNTLZQpAPHskZdOyJL0wQXRtT7H5XToj2+crd9dPE3/9epJ/XxftZ47VK8ZWSnauLxAPQNBTZvivez6rli6K/Dv/ueQGjr6dN9NRQlz1xPi38iaqa4Bm5NMHsoUkGC6+wf15Ju1WuSfoowUt91xLmpquldfuf0qffGDc7Xz4Bn99PVT+s6vjkoa2hBwTm6q5ualqjgvNWq/hyvZefC0fvlmrf7s1rnKH+NB0kAscDqM0r2uhDpShjIFJJin3zqtjr5BXTMry+4oV+RxOXVHqV93lPrV0N6rbz/3ro41dOpYQ6feqW2TNLQZaHFeqq7KT1Nxnj0L6cersaNPX3vyoBb5p+hP1xTrF/tq7I4EhFWWz82aKQCx7XI7Z//bi8eUl+bRzOzY2ssoL92rpTMytXRGpizLUn1Hn47Vd+hYY6f2V53V6yealZGSpJeONmrZzMzLbkBq51TgYCCoL/6sTJ19g3rw40u4ew9xKdPnZmQKQHyqa+1RTUuPfnexvXfwXcx4jk4xxig/3av8dK9unJurwUBQh890aE9ls35VXq9fH67XIv8UfXBBvrJ80TUN+H+eLtcrx5r0TxsWqzgvze44QERkpbh1uq3X7hiThjIFJJC9lWflchgtLZyc3c4ni8vp0CL/FC3yT1FjR5/2VjZr78mzOlTXrhuKc7RqXq48SU67Y+onr1fph7ur9OmbZmvj8sIrfwEQozJ9bpWfbrc7xqShTAEJom8goLKaVi0uyFCy2/5iESm5aR59ePF03Tg3V88dOqPfHm3Um1Utum1hvq6ekWHbiNxrx5r0tzsOac1VefrK7fNtyQBMluzhNVOWZUXdKHgkMFkPJIiymlb1DwZ1zezoX3geDlOSk7RxeaE+d8scZaQk6Yk3a/ToK5Vq6uyb9CyVTV363E/fVFGOT9/bVCpnHO6ZBZwv0+dW32BQPQMBu6NMCsoUkAAsy9IblWeVn+5VYWZ0bNI5WQqzUvSZW+boI6V+1bX16F+fr9C/vXhMA4HgpLz/m6da9LHv75bDSI9+asVlF8YD8SJreMuSRLmjjzIFJIDa1h7VtfXqmtlZCTHk/n4OY3TN7Cz9xa3zVJKfpn/a+a5+76FX9HJFY0Tfd3tZrTY98rpS3E794rMrNSPG7qAEJipz+MaPlgTZuJMyBSSAN0+1yuUwKi3MsDuKrdKTk/QH187U/7t7mbr6B3X3o3t1z3/t1dH6jrC+TzBo6Z9/9a7+fFuZSgsz9NTmG7hzDwklyzc0AtvcNfnT6nZgAToQ54KWpfK6Ns2bmiZvFNzRFg3WLszXqpJc/fC1k3roN8f0oX95SRuWFegT187UkoIpEx69syxLLx5t1P/9zTHtq2rRx5YVaFHBFO08eCbM3wEQ3bJ8HklKmL2mKFNAnKtt6VF776AWTk+3O0pU8bicuu/mOdqwrFD/+nyFtr1xSj/fV6O5eanauLxQv7dkuvKnjO2Yl4FAUDsPntG/v3hc5afbNX2KV9/6/UXatKJQj++tjvB3AkSf99ZMJcY0H2UKiHOH6trkMNJV+ZSpi8nyufWN9Qv1pdvm6em3T+sX+6r1988c1t8/c1ipHpemTfFq2pRkZfvccjiGNgw1knoGAkpxO1V+ul1Hz3SqPxDUnFyfvr1hse4o9XNwMRJamtclp8OoJUEWoFOmgDhmWZYO1rVrTm5qXO8tFQ5p3iTddc0M3XXNDB1r6NCDuyp0uq1Xp9t69OqxJgUs64Kvyfa5tWB6uv7oxlm6ZlaWVpfkycG2B4AcDqPMlCSdZZoPQKw7096rs139umVurt1RYkpxXppuKM4ZeTwYDKqzd1CWJVkaWofmdjmU5nGNrK+qb+/TtjeY0gPOyUxxMzIFIPYdqmuXkTSf9VIhcTkcykiJrjP+gGiX5XOrOUHKFJP6QBwrr2vXzGyfUj38fxOAyZXlY2QKQIxr6uzTmfZefXjRNLujRJ2te07ZHQGIe5k+d8JsjcDIFBCnDtUNndjOlggA7JCV4lZL94CCwQtv3og3lCkgTh2qa1NBZjJrfQDYItPnViBoqaN30O4oEUeZAuJQa3e/alp6tHAao1IA7JE9fD5fIhwpQ5kC4tDh0+em+KbYnARAoho57DgB1k1RpoA4VNHQqWyfWzlpHrujAEhQiXSkDGUKiDP9g0GdaOpScV6q3VEAJLBMX5IkJcT2CJQpIM4cONWi/sEgZQqArbKGp/kS4UgZyhQQZ1451iQjqSiHMgXAPilul7xJDkamAMSelyuaVJiVwsHGAGyXlZIYR8pQpoA40tY9oLdrWpniAxAVMhPkSBnKFBBHdp9oUtCSinMpUwDsl+Vzs2YKQGx5uaJJqR6XCrNS7I4CAMpMYWQKQIx5uaJJ1xVlyekwdkcBgKGRKcoUgFhxqrlbp85268biHLujAICkoTLV3juogUDQ7igRRZkC4sTLxxolSTfOzbU5CQAMSZQjZShTQJx4paJJ06d4NSfXZ3cUAJD03pEyLXF+pAxlCogDgaCl144368a5OTKG9VIAosO5I2Xifd0UZQqIA+/UtqmtZ4ApPgBRJYtpPgCx4pWKofVSN8zJtjkJALxn5Hw+RqYARLuXKpq0cHq6slM9dkcBgBGZI2umKFMAolhX36AOnGrRjXPZEgFAdElyOpTmdcX9+XyUKSDG7als1kDA0k3FrJcCEH2yE2DjTsoUEONermiSx+XQ8lmZdkcBgAtkp3ooUwCi2ysVTbpmdpa8SU67owDABbJ9bjV19tkdI6IoU0AMO9PWq4qGTt3EeikAUSo71R33a6ZcdgcAMD5b95wa+Xh/VYskqaN3cNTnASBaZPuGpvmCQUuOOD2EnZEpIIYda+iQz+PS1HSv3VEA4KKyU90KBC219cTvkTKUKSBGBS1Lxxq7VJzrk4MjZABEqXP73zV3xe+6KcoUEKPq23vV1Teo4rw0u6MAwCXlDO+C3tQZv+umKFNAjDrW0ClJKs5LtTkJAFzayMgUZQpAtDnW0Km8NI+mJCfZHQUALik7dWhkimk+AFFlIBBUZVMXo1IAol5milvGMM0HIMpUNXdrMGhRpgBEPafDKCvFreY43riTMgXEoGMNHXIao9k5PrujAMAVZae6WTMFILoca+hUYVaKPC6OkAEQ/bJ9HtZMAYgenX2DqmvrZYoPQMxgZApAVDneOLQlwlzKFIAYkZPqievDjilTQIw51tCp5CSn/JnJdkcBgDHJ8rnV3juo/sGg3VEigjIFxBDLsnSsoVNFHCEDIIac22vqbFd8TvVRpoAYcryxS209A6yXAhBTsn1Du6DH61QfZQqIIa9UNEqS5nIeH4AYkjOyCzojUwBs9sqxJmX53MoaPjgUAGLBe+fzMTIFwEYDgaBeP3GWKT4AMWfkfL443R6BMgXEiLLqVnX2Dao4lzIFILakeVxyOx1M8wGw18sVTXIYaQ5lCkCMMcYMb9zJNB8AG71c0ajFBRlKdnOEDIDYk53qZmQKgH3aegb0VnWrbpqbY3cUAJiQbJ+HkSkA9tl9vFlBS7qxmDIFIDZlp7rVxAJ0AHZ5qaJRPrdTS2dk2h0FACYkJ9Wj5q4+WZZld5Swo0wBUc6yLL14pEE3FOfI7eKvLIDYlO1zq3cgqO7+gN1Rwo7/MgNRrqKhU3VtvVpVkmd3FACYsPc27oy/qT7KFBDlXny3QZK0qiTX5iQAMHHZwyc3NP3/9u49Ou6yzuP4+5tMMrknTZMm6b2UpFBKqbRcpShyR2z3yNUjIiKyuIurrnqOIsdFzH8KxAAAFUxJREFU9qKuuusqihfgAGoBAWWLCwcUqIiFUHqh9GJL2rRp2jT3tM2lmVye/WMmJcSWTJPJ/GZ+83mdk5OZ+f068+nT6ZzvPM/ze54u/01CVzElkuBWbWumqiyPqUXZXkcRERkzP6+CrmJKJIF19vazZlebhvhEJOn5eX8+FVMiCWx1TQt9A44PVmmIT0SS29Awnx8X7lQxJZLAVm0PL4mwZHax11FERMYlKyOdvGCAFvVMiUi8OOf407ZmLYkgIr4R3p9PPVMiEic1TZ3s7ejRfCkR8Y3JuZm06mo+EYmXVduaAS2JICL+MTkvmLo9U2Z2mZltM7MaM/vqUY7fZGbNZrYh8nNL7KOKpJZV25u0JIKI+EqJT/fnG7WYMrN04MfA5cB84GNmNv8opz7mnFsU+bkvxjlFUkpXbz9rats1xCcivjI5N0hbVy+Dg/7any+anqkzgRrn3E7nXAh4FFg+sbFEUtvqHa2EBga1JIKI+MrkvEwGHXT09HkdJaaiKaamAXuG3a+PPDbSVWa20cyeMLMZMUknkqJe2NpIXjCgJRFExFf8unBnIEbP8zTwiHOu18z+HngI+NDIk8zsVuBWgJkzZ8bopUX8Y0V1HYPO8fTGBuaW5vLE2nqvI4mIxEzJ0P58nSEqyzwOE0PR9EztBYb3NE2PPHaEc67VOTdUZt4HLD7aEznnfu6cW+KcW1JaquELkaOpa+2mq7ef+RUFXkcREYmp4qH9+Xy2PEI0xdQaoNLM5phZJnA9sHL4CWZWMezuMmBr7CKKpJYtDQdJN6OqLN/rKCIiMTU5d2iYz19X9I06zOec6zez24HngHTgAefcZjO7G3jDObcS+CczWwb0A23ATROYWcS3nHNsaTjI3Cm5ZGWkex1HRCSmJuVkYJaic6acc88Az4x47BvDbn8N+Fpso4mknsZDvbR1hVhaWeJ1FBGRmAukpzEpJ5MWn212rBXQRRLIln0HMdB8KRHxrSn5QZoOHvY6RkypmBJJIFsbDjKjOIf8rAyvo4iITIiKwiz2daiYEpEJsLejh70dPeqVEhFfqyjKZr96pkRkIvxh835AQ3wi4m8VBVm0dYU43DfgdZSYUTElkiCe39JIaX6Qkvyg11FERCZMeWEWAI0+6p1SMSWSADq6Q1TXtqlXSkR8r6IwG4CGAyqmRCSGXtjaxMCgUzElIr431DPVcKDH4ySxo2JKJAE8vXEf04qymTYp2+soIiITquJIMaWeKRGJkbauEK+83cKVp1WQZuZ1HBGRCZUbDFCQFWC/iikRiZVnNzXQP+hYdtpUr6OIiMRFRWG2eqZEJHZWbtjH3NJczZcSkZRRXpilnikRiY39Bw7z+q42lp02DdMQn4ikiKlFWZqALiKx8fuN+3AOPnJahddRRETiprwgm5bOEL39/li4U8WUiIdWvrmPBdMKOKE0z+soIiJxM3RFX9PBXo+TxIaKKRGP1LZ0sbH+gCaei0jKKffZ8ggqpkQ88vSb+wC4cqGKKRFJLVOL/LVwp4opEQ8451j55j7OnF3M1CIt1CkiqaXcZ1vKqJgS8cDWhkPUNHXykUXqlRKR1JMXDJAf9M/CnSqmRDzwu/X1BNKMKxaUex1FRMQT5YX+WR5BxZRInIX6B3ly3V4uOrmMyXlBr+OIiHiioihbPVMiMjYvbG2krSvEdWfM8DqKiIhnKgqy2OeTYirgdQCRVLGiug6AB1fXUpAVYG9Hz5HHRERSTXlhFi2dvYT6B8kMJHffTnKnF0kyHd0h3m7sZPGsSaRp+xgRSWEVhVk4B02Hkr93SsWUSBytq2vHAYtnFXsdRUTEUxWRZWH8MG9KxZRInAw6x9rd7cwtzaU4N9PrOCIinqrw0SroKqZE4mRncxft3X0sUa+UiMiwLWWSf3kEFVMicbJmVxvZGenMn1rgdRQREc/lBwPkZqarZ0pEotPeFWJLw0EWzSgiI13/7UREzMw3a03pU10kDp5cV8/AoGPJ7EleRxERSRgVhVnqmRKR0Q0MOh5cvYtZk3OoKNSmxiIiQ8oLstQzJSKje37zfurbe3j/3BKvo4iIJJSKwiyaDh2mf2DQ6yjjomJKZILd90otM4qzNfFcRGSEiqJsBh00Her1Osq4qJgSmUAb9nSwdnc7nzp3jlY8FxEZodwna02pmBKZQPe/Ukt+MMC12tRYRORvDC3cmezzplRMiUyQfR09PPNWA9edMYO8oPYUFxEZqaIgfFFOsi/cqWJKZII89OounHN88tzZXkcREUlIBdkBcnywcKeKKZEJ0NXbz4rqOi5fUMGM4hyv44iIJCQzY8akHHa1dHkdZVxUTIlMgEfX7OHQ4X5uPm+O11FERBJaVXk+25sOeR1jXFRMicRYd6ife1fVcM4Jk1k8Syuei4i8l6opeexp66Grt9/rKGOmYkokxh5cvYuWzhBfvrTK6ygiIgmvqjwfgJqmTo+TjJ2KKZEYOni4j5/9aScXzCtl8axir+OIiCS8qrJwMbWtMXmH+lRMicTQfX+u5UBPH1+6ZJ7XUUREksLM4hyCgTTeVjElIm1dIR54pZbLF5SzYFqh13FERJJCeppx4pQ8tjUm7zCfVhIUiYEV1XU8u6mBrt5+qsryWVFd53UkEZGkMa8sn1d3tnodY8zUMyUSAwcP9/HazlZOm1FEWUGW13FERJJKZVk+DQcOc6Cnz+soY6JiSiQGnt+8n4FBx4UnTfE6iohI0qkqywNI2nlTKqZExmn1jhbW1XVwfmUpk/OCXscREUk6Q1f0bU/SeVMqpkTG4XDfAHf+bhPFuZlcoF4pEZExmVaUTU5mOtvVMyWSeu5dtYOdLV0sP20qGen67yQiMhZpaUZlWb6KKZFUU9PUyb2rdrB80VQqI13UIiIyNlVT8lRMiaQS5xxf/91bZGWkceeH53sdR0Qk6c0rz6elM0RrZ6/XUY6biimRMfh1dR3VtW3cccXJlOZr0rmIyHhVJvEkdBVTIsdpY30Hdz+9hfOrSrl2yQyv44iI+MK8SDH1dlPyDfWpmBI5Dh3dIT77q3WU5GXyg+sWkZZmXkcSEfGFsoIg+VkBtu1PvmJK28mIRGlw0PHFxzbQdOgwj992LsW5mV5HEhHxDTNjXlk+b2uYT8S/frKqhpe2NfONK+ezaEaR13FERHynsiyfbY2HcM55HeW4qJgSicJLf23i+3/YzvJFU7nh7FlexxER8aV5ZXkc6Omj+VByXdGnYkpkFK/XtvHZX69lfkUB3/roqZhpnpSIyEQY2lZmW5KtN6U5UyLvYdPeA3zi/mryszJYvmgaT63f53UkERHfqip/Z3mEpZWlHqeJnnqmRI6hpqmTGx94neyMdG5+/2zygvruISIykUrygpTkZfJWfYfXUY6LiimRo6hr7eYT91eTZsbN582hKEdX7omIxMPSylJe2tZM38Cg11GipmJKZIStDQe56qer6Q4N8PDNZ1KSpxXORUTi5dJTyjjQ08ea2javo0RNxZTIMGt2tXHtz14lzeDx285h/tQCryOJiKSU86tKCQbSeG7zfq+jRE3FlEjEC1sbueG+akrzgjz52XOPXFUiIiLxk5MZ4PyqUp7f0pg0601pRq2ktBXVdUB4+YOVb+6lojCb68+cycvbWzxOJiKSui6ZX8YftjTy1t4DLJye+Iskq2dKUlr/4CD/u2EvT23Yy4lT8rjlvDm6ak9ExGMXnVxGepolzVCfiilJWa2dvTzwyi6qa9tYWlnCjefMJpiR7nUsEZGUNyk3kzNnF/Pc5kavo0RFxZSkpI31HSy75y/Ut3dz7ZLpXL6ggjStbC4ikjAuPaWMmqZOdjQn/sbHKqYkpQwOOn7x8k6uunc1g85x6/knsGjGJK9jiYjICBefUg7A80nQO6ViSlJG86FebnpwDf/+zFY+dNIUnv38UqZPyvE6loiIHMW0omxOnVaYFPOmNNNWfM85x3ObG7nzqbc4dLiff/u7BXz8rJnasFhEJMFdekoZ33t+O40HD1NWkOV1nGNSz5T4Wn17N7c89Aa3/WotpflZrLz9PG44e5YKKRGRJHDpkaG+xO6dUs+U+FKof5DPP7qeP24Nj7VfvqCcc+eWsHZ3O2t3t3ucTkREonHilDzmVxTwoxdr+PDCqRTnJuY+qeqZEl8J9Q+yorqOC763imc37WduaR5fuKiKpZWlpKepN0pEJJmYGd+9ZiHt3SHu+O1bCbsiunqmxBd6QgM8ua6ee1ftYG9HD++bWcTF88uonJKnIT0RkSR2ytRCvnzJPL717F95/I16rj1jhteR/oaKKUlazjne2N3Ok2vr+f3GBjp7+zl9ZhHf+uipLK0s4ZHX93gdUUREYuAzS09g1bZm7np6M2fOKWZ2Sa7Xkd5FxZQkla7efl6vbeOVmhZe2NrIrtZucjLTueLUCq5ePJ2z5hSrJ0pExGfS0ozvX3sal/3gZb7w2AYev+0cMtITZ6ZSVMWUmV0G/A+QDtznnPv2iONB4GFgMdAKXOec2xXbqJJq+gYGqWnqZGvDQbbsO8ib9R2sr+ugf9ARSDNml+Ry9enTOWVaAcFAOjubu9jZ3OV1bBERmQBTi7L5j4+eyu0r1nP7inXcccXJzJqcGD1UoxZTZpYO/Bi4GKgH1pjZSufclmGnfRpod86daGbXA98BrpuIwOIPPaEBWrt6aesK0doVoq0zRHNnL3vautnT3sOetm7q27vpGwhPNgwG0jipooBblp7A0soSapo6E+pbiYiITLwrF05ld2s3P3rxbV7Y2sR1Z8zgcx+qpLzQ2zWooumZOhOocc7tBDCzR4HlwPBiajlwV+T2E8A9ZmbOw2n3h/sGqG/vJisjnZzMANkZ6WRlpE34EJBzDufADb9/5DY4wscZ9hi8+3F3lD9HlOe5d06MHD/2n3vX6492XuR2b/8goaGfgUF6+8K/Q/2D9PYP0Nk7QFdvf/gn1E9X7wCdkfudvf20doZo6wrR0zdw1PYryskgNzPApJwMzjmhhPLCLCoKsyjJCx65Gm93a7cKKRGRFPWPF5zINYunc89LNTzyeh1PrK3nixdXcdsH5nqWKZpiahowfCZvPXDWsc5xzvWb2QFgMtASi5BjsbO5iyt++Oe/eXxkLTWytBpZbI08fqwiSd4tM5BGXjBAbjCd3MwAucEAk3IySTNj9uQccoOByPEAuZnp4d/BAFkZ6V5HFxGRBDelIIu7ly/gM0tP4Ad/fJtJORme5rHROo/M7GrgMufcLZH7nwDOcs7dPuycTZFz6iP3d0TOaRnxXLcCt0buzgO2xeovEgcleFgcpiC1d3ypveNL7R1/avP48mN7z3LOlR7tQDQ9U3uB4Ys6TI88drRz6s0sABQSnoj+Ls65nwM/jyZxojGzN5xzS7zOkSrU3vGl9o4vtXf8qc3jK9XaO5qJJ2uASjObY2aZwPXAyhHnrAQ+Gbl9NfCil/OlREREROJl1J6pyByo24HnCC+N8IBzbrOZ3Q284ZxbCdwP/NLMaoA2wgWXiIiIiO9Ftc6Uc+4Z4JkRj31j2O3DwDWxjZZwknJ4MompveNL7R1fau/4U5vHV0q196gT0EVERETk2LRYj4iIiMg4qJgaxswuM7NtZlZjZl89yvGgmT0WOV5tZrPjn9JfomjzfzazLWa20cxeMLNZXuT0i9Hae9h5V5mZM7OUuRpnIkTT3mZ2beQ9vtnMVsQ7o59E8Xky08xeMrP1kc+UK7zI6Rdm9oCZNUWWRzracTOzH0b+PTaa2enxzhgvKqYihm2bczkwH/iYmc0fcdqRbXOA/ya8bY6MUZRtvh5Y4pxbSHh1/f+Mb0r/iLK9MbN84PNAdXwT+ks07W1mlcDXgPc7504BvhD3oD4R5fv7TuA3zrn3Eb5Q6ifxTek7DwKXvcfxy4HKyM+twL1xyOQJFVPvOLJtjnMuBAxtmzPccuChyO0ngAttoven8bdR29w595Jzrjty9zXC65zJ2ETzHgf4V8JfFA7HM5wPRdPenwF+7JxrB3DONcU5o59E094OKIjcLgT2xTGf7zjnXiZ8Bf+xLAcedmGvAUVmVhGfdPGlYuodR9s2Z9qxznHO9QND2+bI2ETT5sN9Gnh2QhP526jtHemGn+Gc+794BvOpaN7fVUCVmf3FzF4zs/f6li/vLZr2vgu4wczqCV+h/rn4REtZx/sZn7SiWhpBxGtmdgOwBPiA11n8yszSgP8CbvI4SioJEB4C+SDhXteXzexU51yHp6n862PAg86575vZOYTXR1zgnBv0OpgkN/VMveN4ts3hvbbNkahF0+aY2UXA14FlzrneOGXzo9HaOx9YAKwys13A2cBKTUIfs2je3/XASudcn3OuFthOuLiS4xdNe38a+A2Ac+5VIIvwHnIyMaL6jPcDFVPv0LY58Tdqm5vZ+4CfES6kNJ9kfN6zvZ1zB5xzJc652c652YTnqC1zzr3hTdykF81nylOEe6UwsxLCw3474xnSR6Jp7zrgQgAzO5lwMdUc15SpZSVwY+SqvrOBA865Bq9DTQQN80Vo25z4i7LNvwvkAY9H5vrXOeeWeRY6iUXZ3hIjUbb3c8AlZrYFGAC+4pxTb/cYRNneXwJ+YWZfJDwZ/SZ9IR47M3uE8JeBksg8tH8BMgCccz8lPC/tCqAG6AY+5U3SiacV0EVERETGQcN8IiIiIuOgYkpERERkHFRMiYiIiIyDiikRERGRcVAxJSIiIjIOKqZEJCGYWZGZ/UOqvr6IJC8VUyKSKIqA4ypmzCzdy9cXEQEVUyKSOL4NzDWzDWa2xsx+P3TAzO4xs5sit3eZ2XfMbB1wTeT+N81snZm9ZWYnRc4rNrOnzGxjZBPhhZHH7zKzLw977k1mNnvE6383bn9rEUl6KqZEJFF8FdjhnFsEfGWUc1udc6c75x6N3G9xzp0O3AsMFUrfBNY75xYCdwAPR/v6zrnRXl9E5AgVUyKSjB4bcf+3kd9rgdmR2+cBvwRwzr0ITDazgrikE5GUomJKRBJRP+/+fMoacbxrxP3eyO8BRt9zdLTnFhE5LiqmRCRRHALyI7d3A/PNLGhmRcCFY3i+PwMfBzCzDxIeCjwI7AJOjzx+OjDnKK8vIhK10b7BiYjEhXOu1cz+YmabgGeB3wCbgFpg/Rie8i7gATPbSHjH+k9GHn8SuNHMNgPVwPajvb7mTYlItMw553UGERERkaSlYT4RERGRcVAxJSIiIjIOKqZERERExkHFlIiIiMg4qJgSERERGQcVUyIiIiLjoGJKREREZBxUTImIiIiMw/8DhFK7J2RpfgUAAAAASUVORK5CYII=\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "plt.figure(figsize=(10, 10))\nsns.distplot(results_mc['turnout'], \n axlabel='turnout (Monte Carlo simulation)')\nfor v in [0.7, 0.75, 0.8, 0.85, 0.9, 0.95]:\n plt.plot([v, v], [0, 3], label=str(v))\nplt.legend()",
"execution_count": 16,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 16,
"data": {
"text/plain": "<matplotlib.legend.Legend at 0x7fa46dbc16d8>"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 720x720 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlMAAAJNCAYAAAAVnfADAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeZSUV4H//8+ttbur9wXoru4Gmg7NvgRISDIxQNQoKuqISNTETEad0fg7o6Nfx234Ovk6Gs+oGU1mdJxJJjNGQMfsZlGykIUkIISdBhroQHfT+75VL1XP7w+annSA0FBV/dTyfp3DSVfXU1Ufclg+3Hufe41lWQIAAMDlcdgdAAAAIJ5RpgAAAMJAmQIAAAgDZQoAACAMlCkAAIAwUKYAAADC4LLrg/Pz861p06bZ9fEAAADjtmvXrhbLsgrO95xtZWratGnauXOnXR8PAAAwbsaYkxd6jmk+AACAMFCmAAAAwkCZAgAACINta6YAAEB8GRoaUm1trQKBgN1RoiYlJUXFxcVyu93jfg1lCgAAjEttba0yMjI0bdo0GWPsjhNxlmWptbVVtbW1mj59+rhfxzQfAAAYl0AgoLy8vIQsUpJkjFFeXt4lj7xRpgAAwLglapE663J+fpQpAAAQV5555hlVVFSovLxcd9111znPf+UrX9GiRYu0aNEizZw5U9nZ2VHNw5opAAAQN4LBoO644w5t2bJFxcXFWrZsmdasWaM5c+aMXnP33XePfn3PPfdo9+7dUc3EyBQAAIgbO3bsUHl5ucrKyuTxeLR+/Xo99thjF7x+06ZNuvnmm6OaiTIFAADiRl1dnUpKSkYfFxcXq66u7rzXnjx5UtXV1Vq1alVUMzHNBwAALt3T35Aa9kf2PafMl95/7hqoy7V582atXbtWTqczYu95PoxMAQCAuOH3+1VTUzP6uLa2Vn6//7zXbt68OepTfBIjUwAA4HJEcATpUixbtkxVVVWqrq6W3+/X5s2btXHjxnOuO3z4sNrb23XNNddEPRMjUwAAIG64XC7de++9uummmzR79mytW7dOc+fO1YYNG/T444+PXrd582atX79+QvbFMpZlRf1Dzmfp0qXWzp07bflsAABw6SorKzV79my7Y0Td+X6exphdlmUtPd/1jEwBAACEgTIFAAAQBsoUAABAGChTAAAAYbhomTLGpBhjdhhj9hpjDhpj/uE813iNMb8xxhwzxmw3xkyLRlgAAIBYM56RqQFJqyzLWihpkaT3GWOWv+2av5TUbllWuaS7Jf0wsjEBAABi00XLlHVGz8hD98iPt++n8GFJ/zXy9e8k3WgmYmMHAACQdJ555hlVVFSovLxcd9117uahX/nKV7Ro0SItWrRIM2fOVHZ29uhzTqdz9Lk1a9ZEJM+4dkA3xjgl7ZJULulfLMva/rZL/JJqJMmyrGFjTKekPEktEUkJALDX0984898o7Xr9wx1nJjT+7qq/i/h7N3z/+5KkKd/6VsTf++XfHpUkXb9uZsTf+4UHfilJWnnb5yP+3h1PHJckZX9oRsTfO9qCwaDuuOMObdmyRcXFxVq2bJnWrFmjOXPmjF5z9913j359zz33aPfu3aOPU1NTtWfPnohmGleZsiwrKGmRMSZb0iPGmHmWZR241A8zxnxe0uclqbS09FJfDgCwS6QPtH2bw22Ho/beA5XRe++Wmp6LX3SZmk6eiNp7D57ujdp7R9uOHTtUXl6usrIySdL69ev12GOPjSlTb7Vp0yb9wz+cs9w7oi7pbj7LsjokvSDpfW97qk5SiSQZY1ySsiS1nuf1v7Qsa6llWUsLCgouLzEAAEhadXV1KikpGX1cXFysurq681578uRJVVdXa9WqVaPfCwQCWrp0qZYvX65HH300IpkuOjJljCmQNGRZVocxJlXSe3TuAvPHJX1G0muS1kp63rLrnBoAABB1P9zxw4iPKM7KnRXRqd7Nmzdr7dq1cjqdo987efKk/H6/Tpw4oVWrVmn+/PmaMSO86c7xjEwVSnrBGLNP0p8kbbEs6/fGmDuNMWdXbt0nKc8Yc0zS30r6RlipAAAAzsPv96umpmb0cW1trfx+/3mv3bx5s26++eZzXi9JZWVlWrFixZj1VJfroiNTlmXtk7T4PN/f8JavA5I+HnYaAAAQF6Jxs8B4LFu2TFVVVaqurpbf79fmzZu1cePGc647fPiw2tvbdc0114x+r729XWlpafJ6vWppadG2bdv09a9/PexM41qADgAAEAtcLpfuvfde3XTTTQoGg7r99ts1d+5cbdiwQUuXLh3d7mDz5s1av3693rpTU2Vlpf7qr/5KDodDoVBI3/jGNy64cP2SMoX9DgAAABNo9erVWr169Zjv3XnnnWMef/e73z3ndddee63274/8namczQcAABAGyhQAAEAYKFMAAABhoEwBAACEgTIFAAAQBsoUAABAGChTAAAgrjzzzDOqqKhQeXm57rrrrnOeP3XqlFauXKnFixdrwYIFeuqpp6KahzIFAADiRjAY1B133KGnn35ahw4d0qZNm3To0KEx13zve9/TunXrtHv3bm3evFlf/OIXo5qJMgUAAOLGjh07VF5errKyMnk8Hq1fv16PPfbYmGuMMerq6pIkdXZ2qqioKKqZ2AEdAADEjbq6OpWUlIw+Li4u1vbt28dc893vflfvfe97dc8996i3t1fPPvtsVDNRpgAAwCVr+P73NVB5OKLv6Z09S1O+9a2w32fTpk267bbb9NWvflWvvfaabrnlFh04cEAOR3Qm5JjmAwAAccPv96umpmb0cW1trfx+/5hr7rvvPq1bt06SdM011ygQCKilpSVqmRiZAgAAlywSI0iXY9myZaqqqlJ1dbX8fr82b96sjRs3jrmmtLRUzz33nG677TZVVlYqEAiooKAgapkYmQIAAHHD5XLp3nvv1U033aTZs2dr3bp1mjt3rjZs2KDHH39ckvTjH/9Y//7v/66FCxfq5ptv1gMPPCBjTPQyRe2dAQAAomD16tVavXr1mO/deeedo1/PmTNH27Ztm7A8jEwBAACEgTIFAAAQBsoUAABAGChTAAAAYaBMAQAAhIEyBQAAEAbKFAAAiCvPPPOMKioqVF5errvuuuuc50+dOqWVK1dq8eLFWrBggZ566ilJ0ptvvqnU1FQtWrRIixYt0l//9V9HJA/7TAEAgLgRDAZ1xx13aMuWLSouLtayZcu0Zs0azZkzZ/Sa733ve1q3bp2+8IUv6NChQ1q9erXefPNNSdKMGTO0Z8+eiGZiZAoAAMSNHTt2qLy8XGVlZfJ4PFq/fr0ee+yxMdcYY9TV1SVJ6uzsVFFRUVQzUaYAAEDcqKurU0lJyejj4uJi1dXVjbnmu9/9rh588EEVFxdr9erVuueee0afq66u1uLFi3XDDTfo5ZdfjkgmpvkAAMAle/m3R9VS0xPR98wvSdf162aG/T6bNm3Sbbfdpq9+9at67bXXdMstt+jAgQMqLCzUqVOnlJeXp127dukjH/mIDh48qMzMzLA+j5EpAAAQN/x+v2pqakYf19bWyu/3j7nmvvvu07p16yRJ11xzjQKBgFpaWuT1epWXlydJWrJkiWbMmKGjR4+GnYmRKQAAcMkiMYJ0OZYtW6aqqipVV1fL7/dr8+bN2rhx45hrSktL9dxzz+m2225TZWWlAoGACgoK1NzcrNzcXDmdTp04cUJVVVUqKysLOxNlCgAAxA2Xy6V7771XN910k4LBoG6//XbNnTtXGzZs0NKlS7VmzRr9+Mc/1uc+9zndfffdMsbogQcekDFGL730kjZs2CC32y2Hw6Ff/OIXys3NDT9TBH5eAAAAE2b16tVavXr1mO/deeedo1/PmTNH27ZtO+d1H/vYx/Sxj30s4nlYMwUAABAGyhQAAEAYKFMAAABhoEwBAACEgTIFAAAQBsoUAABAGChTAAAgrjzzzDOqqKhQeXm57rrrrnOeP3nypG688UYtWLBAK1asUG1tbVTzUKYAAEDcCAaDuuOOO/T000/r0KFD2rRpkw4dOjTmmq997Wu69dZbtW/fPm3YsEHf/OY3o5qJMgUAAOLGjh07VF5errKyMnk8Hq1fv16PPfbYmGsOHTqkVatWSZJWrlx5zvORRpkCAABxo66uTiUlJaOPi4uLVVdXN+aahQsX6uGHH5YkPfLII+ru7lZra2vUMnGcDAAAuGQvPPBLNZ08EdH3nDS1TCtv+3zY7/OjH/1IX/rSl/TAAw/oXe96l/x+v5xOZwQSnh9lCgAAxA2/36+amprRx7W1tfL7/WOuKSoqGh2Z6unp0UMPPaTs7OyoZaJMAQCASxaJEaTLsWzZMlVVVam6ulp+v1+bN2/Wxo0bx1zT0tKi3NxcORwO/eAHP9Dtt98e1UysmQIAAHHD5XLp3nvv1U033aTZs2dr3bp1mjt3rjZs2KDHH39ckrR161ZVVFRo5syZamxs1Le//e3oZorquwMAAETY6tWrtXr16jHfu/POO0e/Xrt2rdauXTtheRiZAgAACANlCgAAIAyUKQAAgDBQpgAAAMJAmQIAAAgDZQoAACAMlCkAABBXnnnmGVVUVKi8vFx33XXXOc+fPHlSN954oxYsWKAVK1aotrZ29Dmn06lFixZp0aJFWrNmTUTysM8UAACIG8FgUHfccYe2bNmi4uJiLVu2TGvWrNGcOXNGr/na176mW2+9VZ/5zGf0/PPP65vf/KZ+9atfSZJSU1O1Z8+eiGZiZAoAAMSNHTt2qLy8XGVlZfJ4PFq/fr0ee+yxMdccOnRIq1atkiStXLnynOcjjTIFAADiRl1dnUpKSkYfFxcXq66ubsw1CxcuHD3o+JFHHlF3d7daW1slSYFAQEuXLtXy5cv16KOPRiQT03wAAOCSdTxxXIOneyP6np4in7I/NCPs9/nRj36kL33pS3rggQf0rne9S36/X06nU9KZ9VR+v18nTpzQqlWrNH/+fM2YEd5nUqYAAEDc8Pv9qqmpGX1cW1srv98/5pqioqLRkamenh499NBDys7OHn29JJWVlWnFihXavXs3ZQoAAEy8SIwgXY5ly5apqqpK1dXV8vv92rx5szZu3DjmmpaWFuXm5srhcOgHP/iBbr/9dklSe3u70tLS5PV61dLSom3btunrX/962JlYMwUAAOKGy+XSvffeq5tuukmzZ8/WunXrNHfuXG3YsEGPP/64JGnr1q2qqKjQzJkz1djYqG9/+9uSpMrKSi1dulQLFy7UypUr9Y1vfGPMXYCXnSnsdwAAAJhAq1ev1urVq8d878477xz9eu3atVq7du05r7v22mu1f//+iOdhZAoAACAMlCkAAIAwUKYAAADCQJkCAADjZlmW3RGi6nJ+fpQpAAAwLikpKWptbU3YQmVZllpbW5WSknJJr+NuPgAAMC7FxcWqra1Vc3Oz3VGiJiUlRcXFxZf0GsoUAAAYF7fbrenTp9sdI+YwzQcAABAGyhQAAEAYKFMAAABhoEwBAACEgTIFAAAQBsoUAABAGChTAAAAYaBMAQAAhOGiZcoYU2KMecEYc8gYc9AY8zfnuWaFMabTGLNn5MeG6MQFAACILePZAX1Y0lcty3rDGJMhaZcxZotlWYfedt3LlmV9MPIRAQAAYtdFR6Ysy6q3LOuNka+7JVVK8kc7GAAAQDy4pDVTxphpkhZL2n6ep68xxuw1xjxtjJkbgWwAAAAxb9wHHRtj0iU9JOnLlmV1ve3pNyRNtSyrxxizWtKjkq44z3t8XtLnJam0tPSyQwMAAMSKcY1MGWPcOlOkfm1Z1sNvf96yrC7LsnpGvn5KktsYk3+e635pWdZSy7KWFhQUhBkdAADAfuO5m89Iuk9SpWVZP7nANVNGrpMx5qqR922NZFAAAIBYNJ5pvusk3SJpvzFmz8j3viWpVJIsy/qFpLWSvmCMGZbUL2m9ZVlWFPICAADElIuWKcuyXpFkLnLNvZLujVQoAACAeMEO6AAAAGGgTAEAAISBMgUAABAGyhQAAEAYKFMAAABhoEwBAACEgTIFAAAQBsoUAABAGChTAAAAYaBMAQAAhIEyBQAAEAbKFAAAQBgoUwAAAGFw2R0AAIB4s3H7KUlSd3dgzOO3+uTVpROaCfZhZAoAACAMlCkAAIAwUKYAAADCQJkCAAAIA2UKAAAgDJQpAACAMLA1AgAAMex82y68Hdsw2IsyBQBAGEIhqaVnQC6HkcvpkMth5HUx8ZNMKFMAgKQRqVGegaGgKhu65OsMqH8opM1bjo55PjvVrdbeQa1bWqIpWSmXnRfxgTIFAMA4NXUF9E9/OKJHdtdpOGTpU8EUZaa49PElBQqGLA2FLA0HQ6pq7NFPthzVT5+r0sqKSfrc9dN1dVme3fERJZQpAAAuYnA4pP/cVq2fPVeloaClJVNztKgkWznb22WMNK00Z8z1119RoNaeAe082a7XT7Tq2cpGLSjO0ur5hcpMcY9ex1qnxECZAgDgHbx6rEXfeeyATjT36sZZk/T3H5yjV4+3SpK6TfsFX5eX7tVNc6do1axJevFos1462qwjDd169+zJWl6WJ6fDTNRPAVFGmQIA4DxCIUv3PH9M//zcUU3NTdN/3rZMK2dNkqTRMjUebqdD7549WYtLsvXEvtN6cn+9dp1s14cXFUUrOiYYZQoAgLdp6x3Ul3+zRy8dbdafL/brHz86X6keZ1jvmZfu1WeumaaDp7v05P56/dtLJ9TSM6hvrp6l/HRvhJLDDpQpAADeYl9th77w4Btq7h7QP350nj55VamMicyUnDFG8/xZmjk5Q88fbtJje+q05VCD/s9NFfrk1VOZ+otTbIQBAMCI2vY+ffLft0uSfveFa/Spq6dGrEi9lcfl0PvmTdEzX75e8/xZ+vvHDuoj/7JNe2o6Iv5ZiD7KFAAAkhq7Anrg1TeVnebWQ1+4VguKs6P+meWTMvTrz16tn928WI1dAX30X7fpmw/vU3vvYNQ/G5HDNB8AIOm19Q7q/m3VcjqMfv3Zqyd0o01jjNYsLNLKigL987NVeuDVN/XS0RZ9diCoNG9467QwMRiZAgAktc7+Id33ygkNBy3dft10Tc3z2ZIjI8Wtv//gHD38hWsVDFk6eLpTbYxQxQXKFAAgaQVDlh58/aT6BoP6i+umaXKm/Ue/LCzJ1mNfuk6pHqeONvZo65EmWZZldyy8A6b5AABJ69XjLarr6NfNV5WqOCdN0vjO74u2yZkpmlOYpRPNPXrkUKO6AkNas9BvdyxcACNTAICk1NozoGcrGzW7MFPzijLtjnMOh0Mqn5Su62bk6fUTbfpTdZvdkXABjEwBAJKOZVl6dE+dHCOLv6Ox/cF4vNMoWGvXgCTp/fML1dQ9oMf3ndaUrBSV5KZNVLzzGs/IXbKdOcjIFAAg6ew+1aHjzb26ae4UZaW6L/4CGzmM0SeWligzxaVfbz+p7sCQ3ZHwNpQpAEBS6Q4M6cn99Zqam6arpufaHWdc0rwuferqqeofCmrTjhoFQyxIjyWUKQBAUnnmQIMGgyF9dLFfDpum9y5HUXaqPrrYrzdbe7XlUKPdcfAWlCkAQNJo6x3UnpoOXVOWp0kxsA3CpVpUkqOlU3P0yrFmNXUH7I6DEZQpAEDS2HasRcZI15Xn2x3lsr137hS5nQ49vb/B7igYQZkCACSFjr5B7TzZpoXF2TG/6PydpHtdWjVrko40dutIQ7fdcSDKFAAgSTz4+kkNBS1df0WB3VHCds2MPOX5PHpqfz2L0WMAZQoAkPBClqUHXj2pmZPTJ/QQ42hxORxaPb9QzT0Dev1Eq91xkh5lCgCQ8Fp6BtXSM5AQo1JnzZqSofJJ6XrucCMHItuMMgUASGiWpNOd/Zrnz1RZvs/uOBFjjNEH5hdqcDiknz571O44SY0yBQBIaO19gwoMBvW568tsOzYmWiZnpujK0hxt+lMNWyXYiDIFAEhoDZ0BeV1OfWB+od1RouJdMws0HAzpvleq7Y6StChTAICEVdPWp67+IU3K9MrlTMy/8vLTvfrAgiI9+NpJdfZxbp8dEvNXFgAAkh7dXSdJykv32pwkur64YoZ6B4P6r9fetDtKUqJMAQASkmVZenh3nTJT3UpxJfZfd7MLM3XjrEn6z23V6hsctjtO0nHZHQAAgEjYuP3UmMen2vpU3dIrr8uhxu4BbXvb84nmiyvL9bGfv6pNO2r0l3823e44SSWxqzoAIGntPtUut9MozZMc4wZLpuZoeVmu/v2lExoYDtodJ6lQpgAACWc4GNK+2k7NKcyUw5FY2yG8kztWlquhK6CH36izO0pSoUwBABLO4YZu9Q8Ftbg0x+4oE+rPyvM135+l+16plmVxZt9EoUwBABLO7poOZaS4VD4p3e4oE8oYo9uunaZjTT167Thn9k0UyhQAIKH0DgzrSEOXFpVky5FgO56PxwcWFCrX52GbhAlEmQIAJJR9tR0KWUq6Kb6zUtxOfWJZibYcalRdR7/dcZICZQoAkFB213SoMCtFUzJT7I5im09dXSpJ+vXrJ21OkhwoUwCAhNHWO6ja9n4tLM62O4qtinPS9O7Zk7X5TzUKDLFNQrRRpgAACWN/bYckaX5xls1J7PeZa6eprXdQT+2vtztKwqNMAQASxr66TpXmpiknzWN3FNtdOyNPMwp8+q/XmOqLtuTYFhYAkPCaugKq7wzogwsK7Y4y4d5+lM5Zc4qy9MTe0/rh04f1d++fNcGpkgcjUwCAhLCvrlNG0jw/U3xnXVmSLa/Loe3V7DkVTZQpAEDcsyxL+2s7NT3fp8wUt91xYobX7dSC4iwdqOtS78Cw3XESFmUKABD3Kuu71dwzoAVJfhff+SwuydFgMKQ/HGywO0rCYs0UACDuPbHvtBxGmluUaXeUmDM1L005aW796wvHFRgKSZIWdAckSc++ba3VJ0f2p8KlYWQKABDXLMvS7/edVvmkdPm8jBG8nTFGi0tzdLy5R539Q3bHSUiUKQBAXNtb26matn4t8DPFdyGLS7JlSdpT02F3lIREmQIAxLUn9p6Wx+nQHKb4Ligv3aupuWnafapdlmXZHSfhUKYAAHErFLL05L563VBRoBS30+44MW1xaY6augd0uiNgd5SEQ5kCAMStnSfb1dAV0IcWFtkdJebN92fJ5TB6o6bd7igJhzIFAIhbT+w9rVS3U++ePcnuKDEv1ePU7MJM7a3pEDN9kUWZAgDEpeFgSE/tr9eq2ZOU5uEuvvFYXJqtvsGg+geDdkdJKJQpAEBcev1Em1p7B/WhBUzxjdcVkzLk87rUw27oEUWZAgDEpSf2nla616UVFQV2R4kbTofRnMIM9Q0GmeqLIMoUACDuDA6H9PSBer13zmTu4rtEs6dkypKlwBBTfZFy0TJljCkxxrxgjDlkjDlojPmb81xjjDE/M8YcM8bsM8ZcGZ24AABIL1c1qyswzF18l2HGpHQZSX2DTPVFynhW7A1L+qplWW8YYzIk7TLGbLEs69Bbrnm/pCtGflwt6ecj/wUAIOJ+v69e2WluXVeeb3eUuON2OpTqcal3MCjLsmSMsTtS3LvoyJRlWfWWZb0x8nW3pEpJ/rdd9mFJ/22d8bqkbGNMYcTTAgCSXmAoqD8ebND7502Rx8VqlcuR5nEqGLJU38kGnpFwSb8KjTHTJC2WtP1tT/kl1bzlca3OLVwAAITthcNN6h0M6oPcxXfZ0jxn1plVNnTZnCQxjLtMGWPSJT0k6cuWZV3W/31jzOeNMTuNMTubm5sv5y0AAEnu8b2nlZ/u1fKyPLujxC2nw8jrcqqynjIVCeMqU8YYt84UqV9blvXweS6pk1TylsfFI98bw7KsX1qWtdSyrKUFBdzKCgC4NB19g3quskkfWlgop4O1PuFI8zh1uiOgzv4hu6PEvfHczWck3Sep0rKsn1zgsscl3TpyV99ySZ2WZdVHMCcAAHpiX70GgyGtXVJsd5S45zs71cfoVNjGczffdZJukbTfGLNn5HvfklQqSZZl/ULSU5JWSzomqU/SX0Q+KgAg2f1uV61mTcnQ3KIsu6PEPbfToTyfR4cbupgyDdNFy5RlWa9IesexVMuyLEl3RCoUAABvd6ypW3trOvSdD8y2O0piMNLswky9dqJVA0NBedn89LJxTykAIC78blednA6jDy/iZvFImVWYoWDIUlVTj91R4hplCgAQ84IhS4/srtXKigIVZHjtjpMwpub6lOp26jBbJISFMgUAiHmvHGtRY9eAPnYlC88jyekwKivw6URLr91R4hplCgAQ8x7aVavsNLdWzZ5kd5SEMz3fp46+IbX3DdodJW5RpgAAMa0rMKQ/HGzQmoVF8rpYJB1p0/N9kqQ3GZ26bJQpAEBMe3JfvQaG2VsqWiZnpijF7VA1ZeqyUaYAADHtf3bW6IpJ6ZrvZ2+paHAYo2l5PspUGChTAICYdbSxW2+c6tAnlpXozIEciIbp+T619g6qsStgd5S4RJkCAMSs3/ypRm6n0UcXs7dUNJ1dN7W9us3mJPFpPMfJAAAw4QaGg3r4jVq9d84U/eFgo91xElphVqq8Loe2n2jVmoVFdseJO4xMAQBi0pZDjWrvG9InlpXYHSXhOR1GU/PStIORqcvCyBQAwDYbt5+SJDV1DYx5LEn3b6tWdppbp9r65GC9VNRNz/PpD4ca1dozoLx0dpm/FIxMAQBiTnvvoI439WjJ1ByK1ASZNrJu6k9vMjp1qShTAICYs+tUuyRpSWmOzUmShz8nVSluh14/QZm6VJQpAEBMCVmWdp1s1xWT05Wd5rE7TtJwORxaMjWHO/ouA2UKABBTqhp71Nk/pKVTc+2OknSumpanww1d6uwbsjtKXKFMAQBiys6TbfJ5XZpVmGF3lKRzdVmuLIt1U5eKMgUAiBndgSFV1nfpytJsuRz8FTXRFpVky+N0aHt1q91R4gq/UgEAMWP3qQ6FLDHFZ5MUt1Pz/JnaW9tpd5S4QpkCAMQES5Z2nmzTtLw0FWSwz5Fd5hZlqfJ0l0Ihy+4ocYMyBQCICYGhkFp6BrV0GqNSdppblKnugWHVtPfZHSVuUKYAADGhe2BIXpdD84qy7I6S1OYUZUqSDp7usjlJ/KBMAQBsF7Qs9e5YEcQAACAASURBVA0EzyyAdvFXk51mTs6Q02F08DTrpsaLs/kAALbrHRhWyLKY4rPZ2bMRC9K92nKoUf7sNJsTxQfqPwDAdt2BM1N8/uxUu6NAUmFWiuo7AnbHiBuUKQCAreo6+jUwHFJ6itvuKBhRlJ2q7oFhdQfYCX08KFMAAFvtfLNNxhile1l5EisKs1MkSfWdjE6NB2UKAGCboWBIe2s75PM45TTG7jgYUZh5Zrr1dEe/zUniA2UKAGCbyvouBYZCymCKL6akepzKSXPrNCNT40KZAgDY5o1T7cpKdSvFzV9HsaYoO1X1jEyNC796AQC2aOoKqKqxR4tLs2XEFF+sKcxKVWvvoAJDQbujxDzKFADAFo/srpMl6cqSHLuj4DyKWIQ+bpQpAMCEsyxLv9tVq9LcNOVzqHFMKso6swi9vpOpvouhTAEAJty+2k5VNfVoSSmjUrEqI8Uln8fJ5p3jQJkCAEy43+2qldfl0PxiDjWOVcYYFWWn6jQjUxdFmQIATKjAUFCP7z2t982bohS30+44eAeFWalq6hrQcChkd5SYRpkCAEyo5yqb1Nk/pLVLiu2Ogosoyk5R0LLU1DVgd5SYRpkCAEyoh96o1ZTMFF07I9/uKLiIs4vQ2Qn9nVGmAAATprVnQC8ebdZHFvvldLC3VKzLTffI43SwE/pFUKYAABPm9/vqFQxZ+uhiv91RMA4OYzQp06uWbqb53gllCgAwYR7ZXafZhZmqmJJhdxSMU57Po5ZeytQ7oUwBACZEdUuv9tR06KOLi+yOgkuQl+5VZ9+QhoLc0XchlCkAwIR4dHedjJHWLGSKL57kp3tkSWrrHbQ7SsyiTAEAos6yLD26p07XzsjTlKwUu+PgEuT5zhz309pDmboQyhQAIOp213ToZGufPrKIUal4k58+UqZYN3VBlCkAQNQ9urtOXpdD75s3xe4ouESpHqfSPE61MDJ1QZQpAEBUDQVD+v2+er1nzmRlpLjtjoPLkOfzqLWHkakLcdkdAACQmDZuPyVJOtzQpbbeQeWmeUa/h/iSn+7ViZZeu2PELEamAABRtaemQ2kep66YzN5S8Sov3aPO/iENDrM9wvlQpgAAUTM4HFJlfZfm+bM4PiaO5Y0sQmd7hPOjTAEAouZwQ5eGgpYWFGfZHQVhyB/ZHqGFdVPnRZkCAETNvtpOZaa4NC3PZ3cUhCEv3SNJamVk6rwoUwCAqAgMBXW0sVvz/FlyGKb44lmK2ymf18UdfRdAmQIARMWh+i4NhywtKM62OwoiIN/nYa+pC6BMAQCiYl9th7LT3CrJSbU7CiIgL93LLugXQJkCAERce++gjjX1aIE/S4YpvoSQn+5Rd2BYA8NBu6PEHMoUACDinjnYoJAlpvgSyNntETjw+FyUKQBAxD2x97Ty0z0qzEqxOwoiJM/HHX0XQpkCAERUU1dAr51o1Xx/NlN8CWR0ewTu6DsHZQoAEFFP7a+XZYmNOhOM1+VURoqLO/rOgzIFAIioJ/fXq2JyhiZnMsWXaPJ8XkamzoMyBQCImKbugHaebNf750+xOwqiID/doxbWTJ3DZXcAAED82bj91Hm/v726VZYlhUITHAgTIi/dq96BdgWGgkpxO+2OEzMYmQIARMyh013K83k0OdNrdxREwegdfaybGoMyBQCIiP7BoI4392huERt1Jqr8s3tNsRP6GJQpAEBEVDZ0KWRJc4sy7Y6CKMkdGZnijr6xKFMAgIg4eLpLWalu+TmLL2F5XA5lpbrVxsjUGJQpAEDYBoaDqmrs1pzCTDmY4ktoWaludfQP2R0jplCmAABhO9rYo+GQpbl+pvgSXWaqW12UqTEoUwCAsB083Smfx6lpeT67oyDKslJc6uoflmVZdkeJGZQpAEBYhoIhHW7o1mym+JJCVqpbg8GQAkNsJnYWZQoAEJbjzT0aHA5pbhFn8SWDzFS3JKkzwFTfWZQpAEBYDp7uktfl0IxJTPElg6yRMsW6qf9FmQIAXLaQZelwfZcqpmTI5eCvlGSQSZk6B7/yAQCXrba9X72DQc2awl18ySIjxSUjqZMyNYoyBQC4bIcbuuQw0szJ6XZHwQRxORzyeV2UqbegTAEALtvh+m6V5vqU5nHZHQUTKCvVrS4WoI+6aJkyxtxvjGkyxhy4wPMrjDGdxpg9Iz82RD4mACDWdPQNqqEroNmFGXZHwQQ7s3HnsN0xYsZ4RqYekPS+i1zzsmVZi0Z+3Bl+LABArDvc0C1JqphCmUo2WalM873VRcuUZVkvSWqbgCwAgDhyuKFLeT6PCtK9dkfBBMtMcat/KKjBYTbulCK3ZuoaY8xeY8zTxpi5EXpPAECMGhwO6URzr2ZNyZBh1/Okw15TY0WiTL0haaplWQsl3SPp0QtdaIz5vDFmpzFmZ3NzcwQ+GgBgh2NNZw42nlXIlgjJiF3Qxwq7TFmW1WVZVs/I109Jchtj8i9w7S8ty1pqWdbSgoKCcD8aAGCTww1dSnE7ONg4STEyNVbYZcoYM8WMjPEaY64aec/WcN8XABCbQiFLhxu6dcWkDDkdTPElo8yUkZEpypQk6aIbgxhjNklaISnfGFMr6f9KckuSZVm/kLRW0heMMcOS+iWttyzLilpiAICt9td1qmdgWLO4iy9peVwOpbqdlKkRFy1TlmXdfJHn75V0b8QSAQBi2nOVjTKSKiZTppLZmY072WtKYgd0AMAlev5Ik0pz05TmZdfzZJaZ6mLN1AjKFABg3Jq6AzpQ18VGnVBmiptpvhGUKQDAuL145My2NpQpZKW61TMwrOEQG3dSpgAA47b1SLMmZ3o1JTPF7iiw2dntEbpZN0WZAgCMjyVLL1U1a8XMSex6jtGNO1k3RZkCAIxTd2BY3YFhrZzFpst4yy7olCnKFABgfDr6huRyGF1Xft5DLpBkslIYmTqLMgUAGJeOvkEtnZajjJG/RJHcUtwOeZwORqZEmQIAjMNgMKS+oaBWVEyyOwpihDFGmaludbIAnTIFALi4jr5BSdJKyhTego07z6BMAQAuqr1vSB6nQzMnp9sdBTEkK8VNmRJlCgBwEYPDIXX2Dyk7zc2WCBjjzPl8QwpZlt1RbEWZAgC8o50n2xSyLGWneeyOghiTmepWyJJ6BpJ73RRlCgDwjrYeaZbR/+54DZyVxcadkihTAICLeP5wkzJS3HIyxYe3YRf0M1x2BwAAxJaN20+Nft3aM6BjTT1y5Rg1dgf03FueA7LYBV0SI1MAgHdwuKFbkpTmcdqcBLEozeOU02HU2c+aKQAAzutIQ7cKMrxyO/nrAudyGKPMFJe6AoxMAQBwjsBQUNUtvZo1JcPuKIhh6V6XermbDwCAc1U19ShoWZo1JdPuKIhhPq+LrRHsDgAAiE1HGrqU6naqNDfN7iiIYYxMUaYAAOcRsiwdaejWzMnpcjrYEgEXdnZkykriXdApUwCAc9S29al3MMgUHy4q3etSyJICQyG7o9iGMgUAOMfhhm45jDRzMovP8c583jNbVibzuinKFADgHIcbujU1z6dU9pfCRaRTpihTAICxOvoG1dAVYEsEjIvPe6ZwJ/MidMoUAGCMs7ues14K48HIFGUKAPA2lfVdyvN5lJ/usTsK4kCahzJFmQIAjGrtGdDx5h7N82fJGLZEwMU5HUZpHifTfAAASNLTBxoUsqQFxVl2R0EcSfZd0ClTAIBRv993WgXpXk3JTLE7CuJIsu+CTpkCAEiSGrsC2l7dpvnFTPHh0pwZmQraHcM2lCkAgCTpqf31sixpgZ8pPlyadC9rpgAA0O/31WvWlAxNYooPl8jndal/KKjhUHIeKUOZAgCorqNfu06260MLi+yOgjh0dq+pviSd6qNMAQD05L7TkqQPLaBM4dIl+8adlCkAgJ7YW68FxVkqzUuzOwri0NkylazrpihTAJDk3mzp1f66TkalcNl8jEwBAJLZk/vrJUkfWFBocxLEK0amAABJKxiytPlPp3TV9FwVZafaHQdxyutyyOkwSbvXFGUKAJLYlkONqmnr1+3XTbM7CuKYMSapd0GnTAFAErt/W7WKc1L1njlT7I6COOfzOlkzBQBILgfqOrWjuk23XTtNTgfHxyA86V6XegcpUwCAJHL/tmr5PE6tW1ZidxQkAJ/HxcgUACB5NHUH9MTe0/r40hJlprjtjoMEcHbNlGVZdkeZcJQpAEhCD752UsMhS5+5dprdUZAgfF6XhoKWBoeT73w+yhQAJJnAUFAPbj+lG2dN0vR8n91xkCCS+UgZyhQAJJnH95xWW++gbr9uut1RkEB8SbxxJ2UKAJJI3+CwfvpcleYUZuqaGXl2x0EC+d+RqeTbuNNldwAAwMT52XPHVNfRr39ev0jGsB0CIsfndUpiZAoAkMCONnbrP14+oY8vKdayabl2x0GCGT3sOAn3mqJMAUASsCxL33n0gNJTXPrm6tl2x0ECcjsd8rocLEAHACSmh9+o047qNn3jfbOU6/PYHQcJKlnP56NMAUCC6+gb1PefqtSVpdlat5TdzhE9Pm9y7oLOAnQASGCWZem7jx9UR/+QfvWR+XJwBh+iKN3rUmvvgN0xJhwjUwCQwP5163E9uue0/ubGKzSnKNPuOEhwZ0am2BoBAJAgntpfr3/6wxF9eFGR/r9V5ZKkjdtP2ZwKiSzd61LfwLCCIUvOJBoFZWQKABLQnpoOfeU3e7Rkao5++LEF7CmFCZHudcrSmXV6yYQyBQAJpq6jX5/9r50qyPDq325ZohS30+5ISBJn95pq7U2uMsU0HwAkiI3bT6k7MKT/eLla3YEhferqGfrjwUa7YyGJnD1SpqVnQDMnZ9icZuIwMgUACaJnYFj3vVKtjv5B3XrNNE3OTLE7EpLM6MhUT3KNTFGmACABtPcO6v5XqtXeN6jPXDNN0/N9dkdCEkofLVPJtT0CZQoA4lxn35A+fd92tfQM6Jbl01RWkG53JCSpVI9TRlILI1MAgHjRNzis2x7YoarGHn16+VSVT6JIwT4OY+RLwo07WYAOAHFqKBjSF3/9hvbWdOjnn16SdOtUEJvSPE61cTcfACDWvH2zzZBl6aFdtdpd06GPLvJTpBAzfF6X2nuH7I4xoZjmA4A49IcDDdpd06F3z56sZdNz7Y4DjErzONXGpp0AgFi27ViLXj7WouVleVpZUWB3HGAMn8fFDugAgNh1vLlHT+2v15zCTH1wQSHHxCDmpHmcau8bUihk2R1lwlCmACBOdPYPafOOU8pP9+rjS4rloEghBqV5XQqGLHUHhu2OMmFYgA4gLG9fGH0+n7y6dAKSJLbhYEgbt5/UUMjSp5aXyst5e4hRPs+ZX5ttfYPKSnPbnGZiUKYAjJtlWapt79cbp9p1uiOgggyvjjX1KDPFpew0jzwuBruj5cn99app79cnryrVpAyOiUHsSvOcqRbtfYOaruTYiZ8yBSSp8YwonVXV1K3Xjreqpq1PvYPB817jNEblk9I1z5+p2YWZo3+gInwPv1Gr7dVtuv6KfM3zZ9kdB3hHPu+Zkan2JNprij/tAFxQ/2BQTx2o166T7cpKdatiSoZKctNUkpOmPJ9HPQPD6goMq6t/SHUd/Tp4ulMPvdEth6nTFZMydP3MfE3PS45/mUZLdUuvvvPoAU3P9+m9c6bYHQe4qLP/kEqmjTspUwDO63B9lx7dU6eegWHdMLNAq2ZNkts5dhrP63YqL90rSVpYkq33z5ui0x0B7a/r1Bun2vUfL1erNDdNhdkpWlkxiTvPLtHgcEh/s3m3PC6H1i0tkdPB/z/EvrSRNVMdfcmzcSdlCsA5th5p0h8PNWpKZoo+vXyqinPSxvU6Y4z8Oany56TqxtmTtOtku16qatbtD+xUUVaK3j+/UDMucAgvi9TP9ZMtR7WvtlO/+PSSpPpXPuKb1+WQ22mSauNOVosCGOPQ6U798VCjFhZn6YsrZ4y7SL2d2+nQ8rI8ffU9FVp7ZbH6h4K675Vq/eq1N9XSk1yHoF6Obcda9G8vHdcnry7V++YxvYf4YYxRTponqdZMXbRMGWPuN8Y0GWMOXOB5Y4z5mTHmmDFmnzHmysjHBDARGroC+u2uWhXnpOrPryyWyxH+v7ecDqMrp+boy++eqZvmTNaJll799NkqPbnvtPovsJg92bX1Duorv9mjsnyf/v4Dc+yOA1yynDRPUo2mjudPygckve8dnn+/pCtGfnxe0s/DjwVgovUNDuvB10/K63ToU1dPPWd9VLjcToduqJikv33PTF05NVuvHm/Vj7cc0WvHWxRMop2SL8ayLH39d/vU0Tekn928WKke9pNC/MnxuZNqzdRF/7S0LOslSW3vcMmHJf23dcbrkrKNMYWRCggg+oIhS5t31Kizf0ifWj5VWanR22gvI8Wtjy4u1pdWlWtKVoqe2Fevnz1XpReONEXtM+PJg9tP6dnKRv3d+2dpbhHbICA+5fo8rJm6RH5JNW95XDvyPQBx4uWqZh1r7tFHFhWpNPfy1khdqsKsVP3lddN1y/KpClmW/uI//6Rb79+ho43dE/L5sehoY7e+9/tDumFmgf7i2ml2xwEuW7KtmZrQu/mMMZ/XmalAlZZy5w4QC/oGhvXi0WbNKczUkqm5E/rZxhjNLszUFZPTNTgc0s+eq9L7/vklffLqUn3l3TNHt11IdBu3n9JQMKSfbz0ul9Oha2fkafOfai7+QiBG5aR51N43qFDIkiMJtvSIxMhUnaSStzwuHvneOSzL+qVlWUsty1paUFAQgY8GEK4XjzZrcDik98yZbFsGl8Ohz15fpq3/Z6VuWT5Vm3bU6IZ/2qqfbDmqrkByrLt45mCDGroCWnulXxkpyXGeGRJXjs+jkKWkOew4EmXqcUm3jtzVt1xSp2VZ9RF4XwBR1tk/pNdOtGpxabYmZ9p/3luuz6N/+PA8/eHL1+v6K/L1s+eqdP0PX9C/vHBMvQOJ+4fy4YYuvXa8VdfOyFPFlEy74wBhy/Wd+QdBsqybuug0nzFmk6QVkvKNMbWS/q8ktyRZlvULSU9JWi3pmKQ+SX8RrbAAIuv5w42yLOnGWfaNSp1P+aQM/fzTS3SgrlN3bzmqf/rDEd33SrWunp6rq6fnveOByvG2+eeJ5h79dmeNCrNSdNNc9pNCYshJ80g6s83H9PzEP1LqomXKsqybL/K8JemOiCUCMCFauge062S7rp6epxyfx+445zXPn6X7blumN0616+4tR/X0gQa9UtWiGyoKdNW0XLkivH3DROsKDOlz/71TDmOish0FYJezZSpZFqHzOxdIUlsqG+VyOLSiIvbXL15ZmqNf/eXV+tz1ZcpL9+r3++r14y1HtaO6LW73qAqGLH158x6dbO3TJ68qVW6MFlrgcpz99dzONB+ARHWgrlP76zq1oqIgrhY7T8/36XPXT9fx5l5tOdSgR/fU6cWjTbpx1mQtLMmOq4OAf/THI3r+cJP+34fnyhmBneaBWJJDmQKQ6O5/pVpel0PXl8fOqNTG7afGdZ0xRuWT0jWjYIaONnZrS2WjfvdGrV6sataHFhRFOWVkPLK7Vj/felw3X1WqT4/cvQgkEp/HKY/Tobbe5LgblzIFJJnO/iE9ub9eC0uy4/qoEmOMKqZkaubkDB083aVnDjbo/m3Vau4Z0Hc+MDsm7k48n8f21Omrv92r5WW5+oc1c2VM/IymAeNljFGOz500a6YoU0CSeXxPnQaGQ1o2wRt0RosxRvP8WaqYkqEXjzbrDwcb9MLhJv3te2bqM9dOi6mpv4ffqNXX/mevrpqeq/tvW/aOdyUC8e7sxp3JgN/JQJL5zc4azS7MVFF2bI7cXC6306F3z56sP375XVoyNUd3/v6QPvbzV3WkITaOp/mfnTX66v/s1fKyPP3nbVcpzcO/ZZHYkqlM8bsZSEAXWn9U19GvA3Vd+tCCwoSdXnr1eKveO2eyJmem6Pf7Tmv1T1/WDRUFWjGzYHQrhYnci+rXr5/Uq8db9dT+es2YlK73zpmiR3af95AIIKHk+jw63NBld4wJQZkCksjON9vkchgtKsmxO0pUGWO0qCRb5ZPS9eS+03r+cJMO1HVqzcIilRWkT1iOjr5BPbj9lCrruzRrSoZuvqqUvaSQNHJ8brX3JccCdH5XA0liKBjS3toOzfNnxfXC80uR7nXpE8tKdes1UzUYDOk/XqnWb3fWqKkrEPXP/tObbVr905d1tKFbq+dN0aeXsyknkktOmkcdI4cdJzpGpoAkcaCuU4GhkJZMTexRqfOZNSVTZfnpevFos16qataqH7+ov7nxCn1qeWnE1y619AzoF1uP6/5t1SrJTdNf3VCm4py0iH4GEA9y0s4cdtwVGFJ2WmJvSkuZApLEzpPtyvV5kuKcrPPxuBx6z5zJWlyarV0n2/WPT1XqX7ce023XTtet10wN+0idlp4B/fKlE/rv197U4HBIH19Sou98cLae2Mu570hOZ3dBb+sdpEwBiH8t3QOqbunVe+dMliNBF56PV366V/91+1XadbJNP996Qnc/e1T/9tJx/fmVft04a7KWl+WNexq0KzCkV6pa9Fxlk57cf1qDwyF9ZJFfX1pVPqFrs4BYlEy7oFOmgCSwu6ZdRmfOuMMZS6bm6j8+k6ujjd36xdbj+t2uWj34+il5XA5dNS1XS6flqCDDqzyfR7k+r5wOo8augBo6A2roCuhAXad2VLdpOGQpK9WtNQuL9Nc3zKBEASNy0s4cVdWeBLugU6aAJHDwdJem5fuUmRo/5/BNlJmTM/STTyzS9/98vv70ZptePNKsF48265+frbrgazwuh8ryffrs9WW6cfYkLS7JHt12AcAZOSNTe22MTAGId03dATV1D+iD0xNjx/NoSXE7df0VBbr+igJ9R9J/v/qmegeD6h0YVu/gsEIhS5mpbmWmuJXmcY7u01XV2KOqxh57wwMx6OyaqWQ4UoYyBSS4Q6fPbJo3pzDT5iTxxeV0KCvVoSxG84DLkuZxyuNyMDIFIP4dPN2l4pzUhL+b5lJcaId4AJFjjFFOmlsdSbBmikl+IIF19A2qrqNfc4uy7I4CIAnlpHmSYmSKMgUksIMjU3xzmeIDYINcnycp1kxRpoAEdqi+S5MzvcrP8NodBUASyvExMgUgjvUMDOvNll6m+ADYJifNrY4kOOyYMgUkqMr6LlmS5hYxxQfAHrkjhx0HE/ywY8oUkKAOnu5Urs+jKZkpdkcBkKRyfCOHHfcn9ugUZQpIQIGhoI439WpOYebo5pIAMNFyk+R8PsoUkIAON3QraFlM8QGw1dkjZShTAOLOkYYu+TxOleSm2R0FQBIbPZ8vwTfupEwBCSYYslTV1KOZkzPkYIoPgI1yfGeOY0r0vaYoU0CC2Vfbob7BoGZOzrA7CoAkx5opAHHpxaPNMpLKJ6XbHQVAkkt1O+VNgsOOKVNAgtl6pFnFOanyeTnHHIC9zhx27FFbD2UKQJxo6x3U3toOpvgAxIxcn4dpPgDx4+WqZlmWKFMAYkauz6NWFqADiBcvHmlWTppb/pxUu6MAgKSRkSnKFIB4EApZevFos941s4AtEQDEDEamAMSNA6c71do7qBtmFtgdBQBG5fk86g4Ma3A4ZHeUqKFMAQnixSPNkqR3UaYAxJCcJNhrijIFJIitR5u1oDhL+eleu6MAwKi8kTLVmsDbI1CmgATQ0Teo3afameIDEHOSYRd0yhSQAF451qKQJa2ooEwBiC156SMjUwm8CJ0yBSSArUealZXq1sLibLujAMAYOWlnylRbz4DNSaKHMgXEubNbIvzZFflyOfktDSC2ZKd5ZMyZExoSFX/yAnGusqFLzd0DWsF6KQAxyOkYOZ+PNVMAYtXWkS0RWHwOIFbl+jyMTAGIXS8eadacwkxNykyxOwoAnFdumoetEQDEpq7AkHadaucuPgAxjZEpADFrW1WLgiFLKyom2R0FAC4oN93DPlMAYtPWI83KSHHpylK2RAAQu3LTPGrvG1IoZNkdJSpcdgcAcGk2bj8lSbIsS08fqFdpbpp+u7PW5lQAcGG5Po+CIUud/UOjZ/UlEkamgDjV2DWgrsCwKiZn2B0FAN5Rou+CTpkC4tTRxm5J0hWUKQAxLtHP56NMAXHqSGO3pmSmKCvVbXcUAHhHZ4+USdTtEShTQBwaGArqZGuvZk5OtzsKAFzU2Wm+RN0egTIFxKHjzT0KWdJMpvgAxIGz03xtvYl52DFlCohDRxp75HE5VJqXZncUALgor8updK+LBegAYoNlWapq7FZ5QbpcDn4LA4gPOT632ilTAGJBU/eAOvqHmOIDEFdyfV5GpgDEhrNbIrD4HEA8yUvg8/koU0CcqWrs0aQMr7LTEm8XYQCJK5EPO6ZMAXGkd2BY1a29TPEBiDtny5RlJd75fJQpII68drxVwZBFmQIQd3J9Hg0Mh9Q3GLQ7SsRRpoA48uLRZnmcDk1jSwQAceZ/95pKvKk+yhQQJyzL0tajTSor8Mnl5LcugPiS50vcw475ExmIEydaelXT1s8UH4C4lHP2sGPKFAC7bD3SLIkjZADEJ0amANjuxaPNKivwja47AIB4ksjn81GmgDjQPxjU6ydatWLmJLujAMBlSfe65HYaRqYA2OP1E60aHA7phooCu6MAwGUxxijX52HNFAB7PH+4SSluh66enmt3FAC4bLk+L1sjAJh4lmXpucpG/Vl5gVLcTrvjAMBly/N5mOYDMPEq67t1ujOgd89mvRSA+JaToOfzUaaAGPdcZaMk6f9v786j46zu+4+/vxqt1uZNkm0Zr9iADcaYfXEMIbSEpDhpFqDQAoEQ8vtlIWk5DU1CgeS06S8/kpSwpIQAIQdCFigY6oSEYsxu8ApeMHjHlrHkRYu1L9/+MY/ssSpZo20ezczndY6PZ3k0z2fuyPJX997n3o8er2JKRJLbmPxs9h9UMSUiCfb8u5WcPLGY0qLcsKOInZMgXwAAGZJJREFUiAzI6Pxs6prbaGnrCDvKoFIxJTKMVdY1seaDai48oSzsKCIiA9a51tSBhtTqnVIxJTKMLXm3EoALNV9KRFJAZzG1L8WG+lRMiQxjz2+oZEJxLrPGF4UdRURkwA6vgq5iSkQSoKm1nVfe38tHTyjFzMKOIyIyYJ378+3XMJ+IJMLrm/fR2Nqu+VIikjIO9UwdTK39+VRMiQxTz2/Yw4jsCGdPGxN2FBGRQTFyRDZmaTrMZ2YXm9lGM9tkZt/q5vlrzKzKzFYHf64f/Kgi6cPdeeHdSs47dqxWPReRlBHJMEbmZaXcKui9FlNmFgHuAT4OzAKuMLNZ3Rz6G3efG/x5YJBziqSVdRW17K5p4mMa4hORFFNSmENVXWoN82XGccwZwCZ33wJgZo8DC4H1QxlMJB09tmwHEB3iM6JrsXQ+JiKSCsqKctmTYsVUPMN85cAHMfd3Bo919Rkze9vMfm9mxwxKOpE0tb6ilkljRlCYmxV2FBGRQVVamEtlbVPYMQbVYE1AfwaY4u5zgD8Dv+zuIDO7wcyWm9nyqqqqQTq1SGrZe7CZD2ubOHFCcdhRREQGXVlRDpV1zXR0eNhRBk08xdQuILanaWLw2CHuvs/dO/vsHgBO7e6F3P1+dz/N3U8rKSnpT16RlLduVw0AsydooU4RST1lRbm0d3hKTUKPp5h6C5hhZlPNLBu4HFgUe4CZjY+5eymwYfAiiqSXtRW1TByVx8gR2WFHEREZdGVFOQDsSaGhvl4noLt7m5l9BXgOiAAPuvs6M7sDWO7ui4CvmdmlQBuwH7hmCDOLpKzqhhZ2VTfyl7PHhR1FRGRIlBXlAtGN3CE1pjPEczUf7r4YWNzlsVtjbt8C3DK40UTSz7qKWkBDfCKSujqLqT21qXNFn1ZAFxlG1lbUMK4ol7EFOWFHEREZEiWFqTfMp2JKZJiorGtix74G9UqJSErLimQwtiBbPVMiMvieW7cHB2aXp8YcAhGRnqTaWlMqpkSGiefWfsjYgmzKCjXEJyKprawohz11KqZEZBAdqG/h9S37mD2hGDMLO46IyJAqK8rVMJ+IDK4/b9hDe4drvpSIpIXSolz2Hmymrb0j7CiDQsWUyDCwaHUFk0aPoHxkXthRRESGXFlRDu6w92BqrIKuYkokZJW1Tby2eS8L507QEJ+IpIWywuhaUx+myCR0FVMiIXv27d10OCycOyHsKCIiCXF44U4VUyIyCJ5evYvZE4o4trQw7CgiIgnRuT9fqiyPoGJKJERb99azZmeNeqVEJK2MKcghw1JnSxkVUyIhWrS6AjP4q5NVTIlI+ohkGCWFORrmE5GBcXeeXr2LM6eOZnyxruITkfRSVpTLnjr1TInIAKzdVcuWvfUsnFsedhQRkYQrK0qdLWVUTImE5OnVu8iKGJecOD7sKCIiCVdWpGE+ERmA9g5n0ZoKzj+ulOIRWWHHERFJuLLCXA40tNLc1h52lAFTMSUSgmVb9lFZ16yr+EQkbXWuNVWZAlf0qZgSCcHvVuykMCeTC48vCzuKiEgoSjvXmqpL/qE+FVMiCVbT0Mrid3az8JQJ5GVHwo4jIhKKw6ugJ3/PVGbYAUTSxWPLdgDw+pZ9NLd1MDIv+9BjIiLpJpW2lFHPlEgCuTvLt+1nwshcJozU2lIikr5GjcgiK2Ip0TOlYkokgXZVN7K7ponTJo8OO4qISKjMjNLC1FhrSsWUSAK9te0AWRFj7jEjw44iIhK6sqIc9mgCuojEq7mtnTU7qzmpvJjcLE08FxEpK8rVMJ+IxG/trhpa2jo4fYqG+EREoLOYUs+UiMTprW0HKCnIYdLoEWFHEREZFkqLcqhraqOhpS3sKAOiYkokAd7bU8eO/Q2cNmUUZhZ2HBGRYaGsMDVWQVcxJZIAv3p9O5EM45RJo8KOIiIybHSuNfVhkg/1qZgSGWLVDS38bsUHzD1mJAU5WidXRKTTuOLoljLJPm9KxZTIEHt02Q6aWjs499ixYUcRERlWSlNks2MVUyJDqKWtg1++to35M8YyLvihISIiUYU5meRlRTTMJyI9e/btCirrmrl+/rSwo4iIDDtmxoSRuew80BB2lAFRMSUyRNydX7yylRmlBXxkhob4RES6M62kgC1V9WHHGBAVUyJD5I0t+1lXUct1503VcggiIj2YXlLAtn31tLV3hB2l31RMiQyRX7yyhTH52XzqlPKwo4iIDFvTSvJpbXd2HmgMO0q/qZgSGQJbqg7y/IZKrjprsvbhExE5iuklBQBsrjoYcpL+UzElMgTuXrKJnMwMrjprcthRRESGtekl+QBJPW9KxZTIINtUWcdTq3Zx9TlTKCnMCTuOiMiwNnJENmPys9UzJSKH/ejP7zEiO5MbF0wPO4qISFKYnuRX9KmYEhlEa3fVsPidD/nCeVMZnZ8ddhwRkaQwrSRfPVMiEnXnnzZSnJfF9fOnhh1FRCRpTC8pYF99C9UNLWFH6RcVUyKDZMX2/SzZWMWNC6ZTlJsVdhwRkaQxLZiEvjlJh/q0hb3IIHj0je088MpWCoJ9ph5btiPsSCIiSaNzeYQtVQc5dfKokNP0nXqmRAbB+5UH2bq3nvOPKyE7U/+sRET6YuKoPLIilrQ9U/qpLzJATa3tPLOmgtH52ZwxZXTYcUREkk5mJIMpY/LZkqST0FVMiQzQfyzdwr76FhaePIHMiP5JiYj0RzJf0aef/CIDsHVvPfe8uIk5E4uZUVYYdhwRkaQ1vaSAHfsbaE3CDY9VTIn0k7vz3afWkhPJ4JKTxocdR0QkqU0rKaC13flgf0PYUfpMxZRIPy1aU8Erm/Zy88XHaSkEEZEBSuY9+lRMifRDTWMr33t2A3MmFnPlmdrMWERkoKYFyyMk47wprTMl0kfuzq1Pr2V/fTMPXXM6kQwLO5KISNIrzstibEGOeqZE0sHvlu/k6dUV3PSxmZw0sTjsOCIiKSNZr+hTMSXSB+/vqePWRWs5Z/oY/u8Fx4YdR0QkpUwvKWDLXvVMiaSsptZ2vvLYKvKzM/nJZXM1vCciMsiml+Szv76FA/XJteGxiimRON3+zHo27qnjR5fNpbQoN+w4IiIp59AefXuTa6hPxZRIHJ5cuZNfv7mDGxdMZ8HMkrDjiIikpM5ianNlcg316Wo+kV5879n1PPzqNqaNzad8ZB6PLdsRdiQRkZRUPiqP7MyMpJuErp4pkaN4f08djy7bzpiCbK48c7LmSYmIDKFIhjF7QhHLtu4PO0qfqJgS6UFlXRPXPPQWmRkZXH32FPKyI2FHEhFJeQtmlrBmZ3VSTUJXMSXSjfrmNq57eDn761v4u7MnMyo/O+xIIiJp4fzjSnGHlzftDTtK3FRMiXRR39zGtQ+/xbqKGu7+m1OYOGpE2JFERNLGSeXFjBqRxYsbK8OOEjcVUyIxDja3cc1Db7J8235+fNlcLjyhLOxIIiJpJZJhzJ9Rwkvv7aWjw8OOExcVUyKBuqZWrn7wTVbuqOauK05h4dzysCOJiKSl848rYe/BZtbvrg07Sly0NIKktc5lDppa23no1a3sqm7k8tMnUdvYpiUQRERCMn9GdD2/pe9VcWL58N8DVT1TkvYO1Lfws6Wb2VXdyBVnTEqKf7giIqmspDCHE8uLWLqxKuwocVExJWltx/4G7l26mdqmVq49dyqzJ6iQEhEZDs6fWcqKHQeoaWwNO0qvVExJ2npmTQUPvLyFnMwMvrzg2EPbGIiISPgWHFdCe4fzWhIskaBiStJOS1sHP/jDu3z116soH5XHlxdMp6QwJ+xYIiIS45RjRlKYm8mLSTDUpwnokla27a3na4+v4u2dNVxxxiROGFdIZkS/U4iIDDeZkQzmzxjL0veqcHfMhu92XvpfRNKCu/PEip184q6X2b6vgZ9dNY9//euTVEiJiAxjC2aW8GFtExv31IUd5ajUMyUpb+veeu54Zh1LNlZxxtTR/OSyuUwYmRd2LBER6cWCmaUAvPBuJcePKwo5Tc9UTEnKamhp454lm/j5S1vJzszgO584gWvPnUokY/h2FYuIyGHjinM5Y+poHnh5K1ecPmnY7pOqMQ5JOc1t7Ty2bAcX3rmUe5Zs5pNzxvPC3y/g+vnTVEiJiCSZOxbOpraxlX9ZvCHsKD1Sz5SkjPrm6KrlD7yyhT21zUwclceXPjKNyWPyeX5D8myYKSIihx0/rojr50/jZ0s389fzJnL29DFhR/pfVExJ0ltXUcMTK3bx5KqdVDe0cs70Mdz5ubls31c/rK/+EBGR+Hz9whn81zsVfPupd/jD1+eTkxkJO9IRVExJUtp5oIE/rv2QJ1buYsPuWrIjGVw0q4zr50/llEmjgOjq5iIikvzysiN8/1MncfWDb3Lvks1846KZYUc6goopSQrNbe2s2lHNko2VLHm3kvf2HARg4qg8Lj15AnMmFjMiO5MNu+vYsHt4X0IrIiJ9t2BmCZeePIH7XtzMX508nmNLC8OOdEhcxZSZXQz8OxABHnD3H3R5Pgd4BDgV2Adc5u7bBjeqpIvW9g6272vg3Q9rWbm9mpU7DrC+opaW9g4yM4wzpo7m2JICjhtXpJXLRUTSyHc/OYuX3q/i0/e+xjcvmsnfnjV5WKwX2GsxZWYR4B7gImAn8JaZLXL39TGHXQcccPdjzexy4N+Ay4YisCQ3d6euuY2ahlb217ewu6aJiupGdtc0squ6kU2VB9m6t57WdgcgJzODOROLufbcKcybPIpzpo+hMDeLx5btCPmdiIhIopUU5vDkl8/hnxet4/Zn1vObtz7gtktnc9a0cCelx9MzdQawyd23AJjZ48BCILaYWgjcFtz+PXC3mZm7+yBm7ZOm1nY2VR4kPyeT/OwII3IyycuKDNml8e5O57v1zvsQ81jM8x69f/j24eM7H8O7fE0Pr8mhr+/hNWM+AY/jNQme7+41O9xpau2gqbWdxpZ2mtqCv1vbaWxtp7Gl44jHGlraqWlspbqxldrGVqobWqhpbKWjm++KnMwMykfmMa2kgPHFeZQW5lBamEtZcQ6ZGdHfOvYdbOGZNbv7+tGIiEgKmVZSwCNfOIPn1u3he8+u5/L73+BLH5nGLZecEFqmeIqpcuCDmPs7gTN7Osbd28ysBhgDhLbV8+aqg3zyp68c9RgzMDh0xZcdeiwouAzoUgx1W/jIIZEMIytiZEcyyIpkkJcdIS8rQmFuJqWFOeRlRxiRFYn+nZ1JUV4WxXlZ5GdHdOWdiIjExcy4+MRxLJhZwn1LNzOnvDjUPAmdgG5mNwA3BHcPmtnGRJ5/gMYSYnGYhtTeiaX2Tqwkbu/JQ/rqi/n1UL30WJ7+7dC1+e1D9srwo58M4YsPjSuT+nu8Rz1+88dTTO0Cjom5PzF4rLtjdppZJlBMdCL6Edz9fuD+OM457JjZcnc/Lewc6ULtnVhq78RSeyee2jyx0q2945kC/xYww8ymmlk2cDmwqMsxi4Crg9ufBV4Ic76UiIiISKL02jMVzIH6CvAc0aURHnT3dWZ2B7Dc3RcBvwB+ZWabgP1ECy4RERGRlBfXnCl3Xwws7vLYrTG3m4DPDW60YScphyeTmNo7sdTeiaX2Tjy1eWKlVXubRuNERERE+i/8ZUNFREREkpiKqS7M7GIz22hmm8zsW908n2NmvwmeX2ZmUxKfMnXE0d7fNLP1Zva2mf23mQ3tddkprrf2jjnuM2bmZpY2V+MMhXja28w+H3yPrzOzxxKdMZXE8fNkkpktMbNVwc+US8LImSrM7EEzqzSztT08b2Z2V/B5vG1m8xKdMVFUTMWI2Trn48As4Aozm9XlsENb5wA/Jrp1jvRDnO29CjjN3ecQXV3//yU2ZeqIs70xs0Lg68CyxCZMLfG0t5nNAG4BznX32cBNCQ+aIuL8/v4O8Ft3P4XohVL3JjZlynkYuPgoz38cmBH8uQG4LwGZQqFi6kiHts5x9xagc+ucWAuBXwa3fw9caFq6u796bW93X+LuDcHdN4iucyb9E8/3N8D3iP6S0JTIcCkonvb+InCPux8AcPfKBGdMJfG0twNFwe1ioCKB+VKOu79E9Ar+niwEHvGoN4CRZjY+MekSS8XUkbrbOqe8p2PcvQ3o3DpH+i6e9o51HfCHIU2U2npt76Ab/hh3/69EBktR8Xx/zwRmmtmrZvaGmR3tt3w5unja+zbgKjPbSfQK9a8mJlra6uvP+KSV0O1kRPrLzK4CTgMWhJ0lVZlZBvAj4JqQo6STTKJDIOcT7XV9ycxOcvfqUFOlriuAh939TjM7m+j6iCe6e0fYwSS5qWfqSH3ZOoejbZ0jcYmnvTGzjwHfBi519+YEZUtFvbV3IXAi8KKZbQPOAhZpEnq/xfP9vRNY5O6t7r4VeI9ocSV9F097Xwf8FsDdXwdyie4hJ0Mjrp/xqUDF1JG0dU5i9dreZnYK8B9ECynNJxmYo7a3u9e4+1h3n+LuU4jOUbvU3ZeHEzfpxfPz5CmivVKY2Viiw35bEhkyhcTT3juACwHM7ASixVRVQlOml0XA3wVX9Z0F1Lj77rBDDQUN88XQ1jmJFWd7/xAoAH4XzPPf4e6XhhY6icXZ3jJI4mzv54C/MLP1QDtws7urp7sf4mzvvwd+bmbfIDoZ/Rr9Mtx/ZvZror8MjA3mof0zkAXg7j8jOi/tEmAT0ABcG07SoacV0EVEREQGQMN8IiIiIgOgYkpERERkAFRMiYiIiAyAiikRERGRAVAxJSIiIjIAKqZEhhEzG2lm/2c4ZzCzPDNbamYRM5tiZm5m3495fqyZtZrZ3f08/xQz+5t+fN04M3vczDab2QozW2xmM/v4GtuC9Z76zczuCBaaHTAze7G3RVPN7CYzGxFzf7GZjezn+f6/mX20P18rks5UTIkMLyOBPhdTZhZJYIYvAE+6e3twfyvwiZjnPwesG8D5pwB9KqaCzcb/E3jR3ae7+6nALUBZvF8fbKczYO5+q7s/PxivFaebgEPFlLtfMoDtaH4KfGtQUomkERVTIsPLD4DpZrbazH5oZueb2bOdT5rZ3WZ2TXB7m5n9m5mtBD4X3L/dzFaa2Ttmdnxw3Ggze8rM3g42050TPH6bmf1DzGuvNbMpXTN0k/FK4OmY+w3AhpgelMsItuwIXneKmb0QnP+/zWxS8PjDZnaXmb1mZlvM7LMxbTA/OP83gh6wH5rZW8FrfKmbTBcArcFCgQC4+xp3f9nMCoLzdrbLwphcG83sEWAtR257gZl9M2iTtWZ2U9cTBrkeDp5/J1gIsvN9fTbmM/rX4L0sN7N5ZvZc0Ht2Y3BMj59xl/PdF7zGOjO7PXjsa8AEYImZLYk559ie3kPwvjeY2c+D1/qTmeUFbbYdGGNm47ppYxHpgYopkeHlW8Bmd5/r7jfHcfw+d5/n7o8H9/e6+zzgPqCzULodWOXuc4B/Ah7pbwaLbtMxzd23dfmax4HLzewYoit5V8Q891Pgl8H5HwXuinluPHAe8EmiRVTn+V8Ozv9jovup1bj76cDpwBfNbGqX858IrOjh/TQBnw7a5QLgzqAnC6L74N3r7rODQqLzfZ5KdLXmM4nuUfhFi25tFGsuUO7uJ7r7ScBDPZx/h7vPBV4GHia6DdVZRD+Xvvi2u58GzAEWmNkcd7+LaFtf4O4XxB7cy3uYAdzj7rOBauAzMV+6Eji3j9lE0pqKKZHk9psu958M/l5BdLgMosXKrwDc/QWiPQ9F/TzfWKL/+Xb1R+Aiotsrdc10NvBYcPtXQZ5OT7l7h7uvp+chub8gur/XamAZMIa+bQZswL+Y2dvA80B5zLm2u/sb3XzNecB/unu9ux8k2q7zuxyzBZhmZj81s4uB2h7O37lNzzvAMnevc/cqoLmPc5s+H/RCrgJmA7N6Of5o72Gru68Obsd+rwBUEu3tEpE4aW8+keGtjSN/6cnt8nx9l/vNwd/t9P7vu7fX7k5jd8e5e4uZrSC699ksIN79E5tjblsPxxjwVXd/7iivs45oj093rgRKgFPdvdXMtnH4PXRtv7i5+wEzOxn4S+BG4PNE55N11fkeOzjy/XYQ/Yx6/RyCnrh/AE4Pzvtwd8f1QWyOdiCvy/kbB/DaImlHPVMiw0sdUBhzfzswy8xygl6MC/vxmi8TLSgws/OJDgXWAtuAecHj84DOobOuGQ5x9wNAxMy6+4/8TuAf3X1/l8df4/CG4FcGeY6m6/mfA75sZllB1plmlt/la14Acszshs4HzGyOmc0HioHKoJC6AJjcy/kJMn7KzEYE5/p019zBvKQMd38C+A5BW/ZDPJ9xEdHCr8bMyoCPxzzX0+fV63vowUyic8hEJE7qmRIZRtx9n5m9amZrgT+4+81m9lui/7ltJTrE01e3AQ8Gw1wNwNXB408QHT5bR3T47L2eMnR5vT8RHUI64oo1d19H91fxfRV4yMxuBqrofef4t4F2M1tDdI7RvxMdhloZzHWqAj7V5dxuZp8GfmJm/0h0ntQ2ole6PQo8Y2bvAMuBd3s5P+6+Muj9eTN46AF379r25cH76vyl9JbeXreHc33Q22fs7mvMbFWQ/QPg1Zin7wf+aGYVsfOmenoPFr3IoFtBwXos0XYSkTiZu4edQUSSSNCL9Q13/9uws8jgCgrSee7+3bCziCQTDfOJSJ+4+0qil+IP5tpWMjxkEh2uFZE+UM+UiIiIyACoZ0pERERkAFRMiYiIiAyAiikRERGRAVAxJSIiIjIAKqZEREREBkDFlIiIiMgA/A/kFz1alCYFnAAAAABJRU5ErkJggg==\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "",
"execution_count": null,
"outputs": []
}
],
"metadata": {
"kernelspec": {
"name": "py3.6",
"display_name": "Python 3.6",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.6.1",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"gist": {
"id": "",
"data": {
"description": "peaks.ipynb",
"public": true
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment