Skip to content

Instantly share code, notes, and snippets.

@bnaul
Created February 21, 2020 06:54
Show Gist options
  • Save bnaul/fffa8f526291966f4693533f90f896fa to your computer and use it in GitHub Desktop.
Save bnaul/fffa8f526291966f4693533f90f896fa to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 57,
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"%matplotlib inline\n",
"import numpy as np\n",
"import pandas as pd\n",
"from sklearn.neighbors import KDTree\n",
"\n",
"polling_centers = pd.DataFrame(10 * np.random.random((20_000, 2)), columns=['lat', 'lng'])"
]
},
{
"cell_type": "code",
"execution_count": 58,
"metadata": {},
"outputs": [],
"source": [
"people = pd.DataFrame(10 * np.random.random((3_400_000, 2)), columns=['lat', 'lng'])"
]
},
{
"cell_type": "code",
"execution_count": 60,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x12e258710>"
]
},
"execution_count": 60,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOydebxVVdnHv2sPZ7hcIUFEEQ1FRREJBCUhJZTBCVMzhfQtIqdEM30dSs0pzTLT1NRyyiwcSnECSQSnFMURQTBRBJVUQMALl3uGPaz3j3P3Yd999njOueCL9/f57A/cfdZe61nTs571rGc9j5BS0oEOdKADHfjqQNncBHSgAx3oQAc2LToYfwc60IEOfMXQwfg70IEOdOArhg7G34EOdKADXzF0MP4OdKADHfiKQdvcBMTBNttsI3v37r25yehABzrQgf9XeP311z+XUnb3vv9/wfh79+7Na6+9trnJ6EAHOtCB/1cQQnzo975D1dOBDnSgA18xdDD+DnSgAx34iqGD8XegAx3owFcMHYy/Ax3oQAe+Yuhg/B3oQAc68BVDB+PvQAc60IGvGDoYfwc60IEOfMXQboxfCHGXEGKlEOJt17uuQoinhBDvtf67dXuV3154f8V6HnztY95fsb7yx1Wr4NVXS//WEaubC7z18Resbi5Ul0E70eXF6uYCzy9eyfOLV1VPaxVl1tQ29cYmbOvQenvoCErvvH9/xfrN1o7V9mHd+34T9d2XAe15getu4I/APa53PwdmSyl/I4T4eevfF7QjDYmwurnA4reXkF6+nJ1792BrWeSDhm14I6exT9Zk2uNzuPsTwZqGLnRtaeKUnQTHffcAPtYa6TFtKtuecwakdBTDYP3Nf+aDUeMwTItlq1vYNt9E03/eJ92nD9ubzeReeIltRx1I7wOGtKFhzdLlrJ3zCts0pugybCiP/rfIBQ/NR1cUipbFGSN34/tDd6JbY7qyAqtWwbJlrBUpli77jEKvnej/zqs0Tj4NS1GgaDD3rEvYfuL3SS94iyWrNtDroGF0+XpPlq/N0WvrbCnf1nzo3Ru6d2d1c6Ht761t5X736Lz/cu4/38KwSvEduueauGJgI4cesX9FHgAvLVnNh6s38PVundi/Tze65daVy1yd7czytTk6pVQ+acoBgr16dm5T59XNBe749wfc8cJS0prCjis+4qqeG+g68gBey/Zg4I5fY9ceW7Vp188XLmabvXZHdu/ehpaFn6wj/9/P2O6Dd+itGXTevjsMGgTdSxce31+xnhfe/5xtGlPs32cbujWmy/WRK1aw7qXX6CnyyDfepPeUO7B1Hc2y+O+lv6bz/vvytX67t2mDLsveR7z6Cp0PHM7WQ75Rbu8122zP20YakPTskuWTpjz/XdtCwbT51q7blOqzahUzHn+Jy+a3sG6rLtgSfnfsAIbvuk25zYwpU9j9onOQugZFg2fPv5ozZF80RVC0JOeO3p1hXSQfv7GIyxe0sKahC0XLJq0KhCK45rsb83P6d3VzgU/f/5heTSsq6uO048pX5tHrvQVsdeBw2HNPVjcXWPhJE+tyJp2zGnv17ALQ5pspcz/i3mmv0Xv9Kj7svC2/+MEBHDlwh41jsLERmpvLY9HBlJc/5LLHF6IKsKTkvLF7MHTnbiV6c+v4YtFilnfpwfa77ug7ttzj8PPmPIPnPMmevzwbRVFRbAtuuIH3j/o+L7y/im0aMxVj1E1LNVizdDnLXl/E5922I7PDdvTskmVD0eLz9XneWt7Egbttw5Cdu9VURhhEewZiEUL0BqZJKfu3/v0u8G0p5adCiO2BZ6WUfaPyGTJkiKzXzV1nsHZKqeQ/XVEeyFOWtPDK1TfzmyduACBjFinqKWwJ9w8Yzfi3ZmKjoEqTqXsdxDHvPIclFFRpc+XoU7ho5p/JWka5nJyWZsTku1iZ6cK4Rc9yzYwbMYVK1sihuuh5+bAJ7D/9XgDevOg39P/NxWi2BYCtalw25jTe6r4Ly7v0YE1DaeLoCvz+uIEcOXCHcn0633MnO15+IVIItGKBnKqjCIFiWaSkVS5PAhYClVK/G8ADAw/hH/sfw8pUI79aOYdRD9+Bkk5DocCiSWfyv9audDZzLGvclp+OH0ZjRisvRoZt88sj+nHF44somDZAub6GopEVNg/95HIuaeiPrirkDQsJ2K3DrmtLE9+f9wQ/m/sgMpVCFoucf+hPmd5vBIa9Mc1O61Zw+LhhHDN2IK+88i633fM0TXqGxmKOCW8+wfgFs8p1fLLPfvx98BH0P/xAfv6DEbxwxY0M/tV5mKqGZhr8Zb+jeHPngXzUqRuZQgv9Pnufy5/6M7q7nVSNd6/5I3/Z6Zv847Xl5ffdWpqYsK3JPz/XGfrBG1w77fo23wlX38rWcZBSYPrPruRcdU+unH4Dx81/qpzmw8HD+frC1zFRMA2D8w77GXN6D6RX0wqaU1kaizmWd+kBwFUrX+TgqXeSFyq6bXLeoWcxrd8IFAGqAFVVaGhay4u3/oisWXSNxRTDf/KX8vhx949uGdw07HjuH3hoWaDZad0KPv3adrR06Yph2xw3uBfNd/+Nq6bfiKUoZKXF0mNP5MZcDzY0dOLt7rtw+r+n8IM3ppfLfP+I4zhswA/LgoDTNk772JTUDce/+QSXzr4NQ9HRpcXN+x/HqAP3YsB1VyClROTzkMkgpIQbboBTT+W255bw6xn/wQsVGPefZ7lmxk0UUdAtk0f2+jbH/OffiJQOhsH5h/2Up/Y+iLxplWnr2tLEyzf/kJRttum7X4yZzP2DDgXgyEXPcu2TNyFSKTTTQFx0EZx6apsFcKt1azGXLkXbeWeW653Q16ym1xcraNquV2kBahWqFl12DbvcdiNFRWvTjw66tjTRq2kF2+69O3ece3hFPZNACPG6lHJIxftNzPi/kFJ+rfX/Aljr/O3z7SnAKQA77bTT4A8/9L15HBurmwtMmfsRf3x6MbaE773eOuBUnTQ2l3z7JC59+vY2E8aBpHJCu/82Wxmp+11OS3H893/D8i49KiaiN+/ZD86m4eWX2P/ai9vk4fy+QcugSZPLDj61PBABLjx0D66btZjvvjqdK5+8ueJbP1qDaACwKE0eb10B8ooOiuDyg07m0623oyHXXPpBQEumkQXb7sKahi4MWr6I++/7BWl7IzP0Mp5dVn/MwE/epWvLOs594e+kzWKbMvNaipOOuZhFPfow9t0Xy4xBkxb//MYYvvfWTKSUZC2DvKqTsQxfmk0Upk04k0P+eWub9nePeL/vy3SrOsNPv9uXYaZsE8UySUk7onVLKAqVf/Q/mBMWzIwcS6amIykJHzlVR7NtEALdNivGmLtdAQZ8upi/338xnYst5Xfr9SwnTLiK+dvvTteWporxKIGClmLqngdyzKJny23tMKSgb8o0KyqabVXU6659xnHz8PFt6IONjG2/j97momfv8u27oLH80MkXcW7X/X1+JZDOqDb71tI3+Ns/Lqkos6ioHDLpj3yR7eybr0xn+Pfki7j7Ezhg8Vy+P28GhppCk1ZJWJz/FBYKum1w3ciJjD54EIMu/1+UQiGQpjaLsm3yzHlXc9ivz/Gtbxx86Rh/699rpZSRev5aJf5H5/2X8x+cX5ZIT5o7tWLAFRQVQ9VpNPIV30cxT7/fJXD6dy5gfboTtz/4K7K24fNlKd2lo07lwmfuJGOZvmncaX8xdjL3D9zI/Lu2NPHSLRNJW8H5RzH+JAgaLYZQmN9jVwZ/thg8Za5LNXDi+CuZv/3uXPrUrUx0SYZBE3yDliZlGejSDmWUcejNK3pg+4dhg5rml2N/wryefem5bhV3TL2STAhTiaIDn/RRQkUY8qrGSd+9hBd23qe8mC7duif3PnBxBZNyxs23lr7Bn6ZeRaNZqRf3lu0swAC3PnI1WxVzofXzq1te1TnvsJ+xqMcurYt9E//7whRModJo5BKPzaKiccikm/gi25leTSva7IQHfLqYKfdfFEqneyw6OPSdf3PLY7/1pb+gatw0bDynzp3aZjF1p3EQ1o+lnba/bj1MSMxpKXKLl9B1516BdQpDEOPf1E7aVgghtnepela2d4Grmwuc/+B8Oq1bS9+mFez30QIuevYvFZ2csi1SLim1Hrhu+vUoUpKKYDqrGr5WVn2EQQC/evIWZu4+rDzYezWtwFDUQMYP9WX+QfmkpM3gzxb7/p6yTXb4YgV9Vy5l4hvTI2kR4MuYqoEAMlUwfYAGq8BVT95CxipSECppWf34qOfi6yBtmdzx0K94ece9GLFsXvn9izv2Z/jHb5fLFMCls29HSLj06dtJB+w+vTSmzSJ/mnoVqrTRYu5svPllLYM/PH5tWbXoV04S6LbJzDtOx1ZU8lqqjapkeZcepCKEp5RtllVnDtZnOgXSn7FMzpzzAIrl3/dx6yKgjYrXjYxZpDmVbZ3LGlk29o+haCx7fWHVjD8Im5rxPwb8EPhN67+PtneBU+Z+xJj5T5d07IpKY9FfynDehW01gxCUn7NyB/F0Cdw3YDTvbtubrF05YP0Ytobk9Dn389heI2lOZRnz7hw6+exSgr5vT/hJTBYKmmlwy2O/3YSUtEVBTYG0SHvUEW4ESd5Zq9SHmRhMv5r2rqV/SozJYMSyeW3yGf7x2xVpTaFy2aw/k/boscPKdy/AMiK9CPndqwaNAwkUUEjTdscnKM0BbJNUsVSX3824gUU9dqGxmOOaA/+Hi5/xVx8BaKbBmHfnsGi7PuXdwqIefbCEEri4pS2DpMteGI/xwzELZnPXfkeje/iAbpsUeu2UsPQYkFK2ywPcB3xK6fxwOfBjoBswG3gPmAV0jZPX4MGDZTX4fH1eDv3ZvbJFS0kJdXvs1ievqNIO+D1uPr8cdao84bgrZE7VAtP4vdugaGU66lm3KHqTpm9BqXu7J613i5aSI0+6VV5zwImygFLxbZI+dMo2hahIUxSKzCmabNbSofQFlVdtf8apT9GH3lr7Osk4r6ZOD/c9QBZKx7qhT05NyZyqyaZUg2zRUnLK3qNlTtXkei0T2NbrUlnZoqXk5HHnyUFnTpEFoda9DknyKiiqHHTmFHnv3qPbjIW79jlCvvrB51XxP1lixK/58uc4jHdzP9Uy/nkfrZXfm/QH2ZRqqFsHtRlwiipzql7TAFivpmSLqgcOvE3J2OvRJk4aG+SjfYfL9Vo68Pe45VXLOGyQBUWTk8edJ79+wTT59QumyV+OOjUyvzhl+jHbGX32lTlVly1aqirG3159VFqoohfgeo+1WvOL2/feNM5Cf8/AQyK/z6u6PO3IC+S6VLaudU/aLi2qLm8aeqxvXR6d+WZV/E9KKYMY/xZ9c7fX1lmWbdW9Yvsk65R/UUtz07DjKSpaVXkKoNEqls1ADURFPptSVeOglvZxvp222/5cNuZ0UgH6ZNunrkEoKMk1ko7poCUEc3oPLL9/wfX/MBii7dSI0y9jl7xKxjLIeqyU6gXpQ4dfGm/ZJpT2ORHfVUuT92/nqfXEzG3+GZXODVOo9Fy3iu+9/XTk9ynL4A/TriVrbLqLa75qICk5+bVHKuuiqAyyvqg7DVs04+/WmOako4dy3qFnkdNSrEs1kNNSTN1zROyBHpZOt01m9B3Or0f8MBFdfhM4r6f56ZHnU90x5OaFuz7OZD3ivZe44l83o3tq6jCmJHpfpYalyFQ0vr3kVbq2NNG1pYnGYo6p/UZEMlE94WGmH5MKyj+pVZLzxKHIL++4k7yaxcr7jXA9QYeZ9UBY3zUaOb750XwMJZoCAaRtC1XagYtYe9AoPc/M3YZSVPSKdGnbZMfBe9VAhT/a1ZyzXqjVnHP8n+eweOEy+q1YAsCiHn0Ys3gOVzz1J/SQAz+olKDcrfVc70EM+2hBmx1FnMmTU3WEbbc5MHRsece8O4erZ/rb5AfRtKngDNK4jMSPzqS0S+C+/qM4duEz5YtSScwfJdCcypI2iggpMRSFrG1SEAppj6moX34SKKg6c3bqz8ilb1bV7nHqXBQKCqD6LDi1jgX3ohw3H+eboqqRssyaxltQ+e65VM3hb9g3JpV3UpJ8XyskYLde8IyDvFoSkdyXQCVw0djJ/O+Dv/e/qR8DQeacW7TEDyVzzjc/bmL4sje5/aFf8aepV/HiLRNpTjVwyKQ/YiVsArdEM2LZm6RaL9UESXx+UkPGMtCkVf6tKFTOO/QsAIqazgY9k7ie7u11rZJKEHKKhhlDiqo3jl/4NLq02rRxErXHVsUcKWmhY9PQ2l8ZH6YfBM22GLn0zTbl1qON3X2mS7vMJGwIHFNexFWFOFY3cfMply1lu+1C3XWM26YSyCtapBrJj+lLoFlNJa5PlOQeBEvK2GMkbRnM2nUoOS3Fej1LXtX4xZjJ3DvwUBZ+0pSA2nj4fxFsvRZMmfsRndat5ffTrm/juuC6adfxzTPu4R97H8SEBbPqvvpL4Oahx/JFdit+/uxf2gzEkknaRthCMOKD1/nd9OsxVb3CPDOOdOL8rsRIWy38TE7diCN9Gwh0JHkthWqXjrJsRSVt+evFBcml4CSI07ZeMz9JaRGMag8HYSaQUNln7SWNJdkpOXDutoSdcUhKErYW8HvUouNOb/u8c5f/eN/hzN++Lxc8/1dkxG7dD51azXPrNX7C+jaVQDQQwKj353L4j26k57qSk7hFPfpElFI9tmjGv7q5wB+ffo/9Vixp408FQJcWQz+cz9HvPO87wNxMOmiShHVHcyrLy18fwO0P/SqykdO2ybELWw+iWplJNYd4taSrRz4WgGt76ydFOpNBSJvjJ/yaTkaeXdb8l/Ofvbtul7aSoFq9dkNMpl8PGjaXao+Qct3j00KAEIgYB8hxBJigNBO/ewk913/Or2feUtUiVo3KbJMqwqXkR689yrELZmOoOrptcsFhZ7FXz1F1L2qLZvzL1+bQ1GD56chFz1XclPODW72QRD0ArZMiAlESod/7ehwcJkEsBpBEHywl9993IUVVJ2Pkqc4uKpwmQ6i+5wJ+6XOqjiYlim2Vb5k6jsSqKdvB5mLYSZF0cXGn1Sjf3QnMzz2HighSCQ73oXTW0GAUuHT2bb4ChbdsC4HiKqPafghaDKpZeKLyydgmJ877F4KNuv7fP3kTeu48aKzNG6gXW7SOv9fWWSxbsqhHHwyPSaAADnn/ZRo8/jfCJI64EnZRqFx+8Mks6tEHrYZr/u2JIJ1qUvbr1tOmExwCpm2LtG2ylZFDT8gE3AhbAFMyvirg7wMPbfWlIsv1qVZtVlB1iopKQQ2Wq4L02fWWMIPK2ZwL0ktfH5D4m5RlctNj1xBmi+Lerbv7sV5jy1FpWYg27Rr3cN2LwHMVF1TLLLmCrjO2aMbfrTHNd76xPWsaunDD8PG+0phKcMdUOwmLqsals29n2LJ5XDbqtE27XfSg1gneXiZttdDl0FQ/ZQuc/PpjZDzqm2rVQBnLKC1sVuX9kSQMoxY63N9uCibvpy4NSjfiw7faME9nF5ATwYYDDjP39lE9ETXOC0LBVLQ2i0p7QzFNKIZrJKrKt+45fomwurnAI/M+AWBG3+FV5RHG+Px+E5T8m2TNIr+bcQMz+w7j/r1HJ2agzsQIKnNuzz18JZIki1hcq5FaB3hR1GKJH1ynMHO9JKhVMnQQ1l/ecsLKcltmxW03ScmNcNJ2lkBOqBU74rD03r+r2SVqSAqqzm8P+AHnjz2TScddztUHnxTrWwk0a5nIdqpWdeWG4VqgNCnb+DqqJe9EbfbGG4nLjMIWzfiXr82ht+r4nZNyL6JUO87jNeOLI7WZQqXfiiUc9c5zVTAWyQYlVfHWQnD9sAl8o/VOghuGULn5m8dixGAAQYNx+u7DKKg6OS2JTUIwDKFy8jEXxzrrCIOblijmWYv+tZYdTtDWPa6K0KFVS/itg3/tOjRB6o3QpcXjfYdXVe8go4d4dAt+9uIULp19G7c/9Cu6tKyL9dUGPcMfhx1XVsVFlVfLOLaFUs7fkfTrVU7sxb1Hj6gUibFFM/5eW2cx4/g7jgH34Io7KTNGnsMXPYuwk7u0VYFOdtstngSuHvkjfvLKQ2R9gofk9TQv7zQAs3Wwur+Ly0AOXzwHW0qUGLEBwn7LKTo5RePB/iO57eGr8C5FSXplU22r61lWEkndOQWqVaUzbvEc3tt6+8RtqwHHvPOcL202lYuuXx4SaNHSicpOW0UytlXaIVsG58y5r42EHQTNtvharnKRqLXfvN9bQmljAu5FLZwlLg+RAHvvXUNJ/tiiGX+3xjSThvcGSjax9ThmjTu4JCUJYfzbT1e9PXQm1AahUVRUrvz2JF7dsT8ELCQpx+ePWnn1O0mZWdskXYX/dW9mGpIJC2b5ukMuqhom9TnM3JxnKEFIor+v190LAey29tPQNEnUIkl2HTk1xY3DjicfoqePKlNA2donqEwJKK3WL7UganGRgC3DJfw4TLt8HiUUDOGvvg3LJ6/qfLFyTURJybFFM36gHOB57Ltz2tV3iBuOhF3Ndt1vYDRIE922uOD5ezh+3pO+B1wS+OfeB7OoR58Kp3Rxy/HCEApmiH4+aKsqKJmjhbnDSFlmVTr6pFZIXom1XoiTZ9xya20Db15RzKo9dk9Zq8g5L9zL7N2Glhf0Wg0DgpCSkkYj30YFGwW/NGHqRwt8/UyFwcnNockSCi1amqKq8Y/+B2FGLIq+dRFKReCYemCLZ/zr8wZdW5q4zMf+Ny6qOTStF9yLR8o2+f78JwPTruxUimJ50/7HkVP1klM6Va8Y4JLSbeEo6NLGIloadA/4uKiGATmTyY2ij8lk2VJEKFhCbRdGFydPb5r2ZIZJEbRgV5OP823KNjl88UYBa1G3+gQQCVM1xRlH7t2z+51KsEsFP6Ek7jjPKRqGoqFJm05mgbRlMmHBLLKumMlBeTiCivPcP2A0me07GH9iFEybXk0rAhlEXtEoKhpFodLS6ijJkXKdxwhpJpP47oXjIO4g9nv/0zn/4MVbJnLa3KkIIfjz0GM4+bu/JOfx/bNBS3Pr0O9G0l2gZL7mRdB3m0IPP3WvkWV/JoZQET5nEeWJrigUterUXnERR2Xgtepxv6+Xqssvn7D8k6iiwsr1S+88/VZ/hCnCWUycNmivcVWvfL30K1JW5dOqtBi1bcPx82eS/3RFzTR6scUz/qJps7xLDzQfm+q8UEHAJaNO5Ztn3MN93xhLkZJvmDaSNsEOvep90l8LNGmRtQw6F1vImEXOeOkf7LFyGY1G2+DTncwCa7OdmTJgTOi2XEHSYFXaENd7IiZhgEcueo6H9xxB2i45ukuF0FNU9cB4AFHM2u/8wV9dEA4/dUJcNaBbP+wIGLI1z6LLWXWYMBBHR12LFVPUWJiz495tBKkgGr5sSKLbd6ucJJCWFtkI9yMmCoUYZs6GotGrqf6MvyIyy5fxqSX04u4XPSHPGHeuLLhCFXqj3ORUXc7uvc+XJtqVH43VRCJap2dlPiBCWFBbxA1vWA193vTrtbRsUXV549BjZS5BeMy4ZRkl7zFV0Ra3vd/Ydpe6h1J0nryqyV+OOlUedcI1suiJ0NZeYTej+rya6Gk2SEMoslgjXVFj1e+JqpshKkNx1rs9/Z5CzNCpLVpKrv7g46r4n5RS8lWMwLV8bY5tck1cM+PGUPfJactg5LI32kWSreabpZ170KLG87/tLsNLv24HB4kOaotqD6XjSK5e3D1kHMNPv5vff3si57YGy9mQ0CQwrDwn2EvSfo0riU/bfX/6r/ygZouYIKQsk5Rp8MB9F1a4/qjGaCBOu8Y5GE6iJnIeTdpQ5UU+CVw3fAJFj7VaQUtx/tgzy+bLXvqcnUZe1Sj6sDoJPLrnAVVQFB9eutzndX6qP++ltKmDxvKx1lh3urZoxt9r6yzbr10ReSuxXpYOUYM67sTbed0KGqx4nir9DhA36BlyWorfHfADMlb9vKm7B+RzvQdRUOO7VvNbYE565eGyPfa0fiP4/vFX8rdBh0WGWvSrs/f3OP0Zlxl6YSgq9w4Yw9j3X63ay2GcsguKxgXP/7XMJGpBe6hUkrad1711VN6S0vnald+exHN99uXSg09pE0nv3EPP4tMu3X3zFZSMEwSg2hLdo66VlC4rHrb4pS+FuskUKkVExaW0Y96cyY5mc93L26K9c3ZrTHPSiSPR/9Z+/j2ikFd1DFUjbRRih/NzSwHNehZNWjy85wiOWfQMQkpSIWaSBRTOPfQs3t22NwM/eZe8liIboOd2w5nEQflKz2/Dlr2FjGEZFMYc0rbFzDtO5+Ixp9P382VMfGN6+TcTkSg8o5c+fP72w9Iu27Fz02e+af3ylAC2zffnz0xEm5cer5WJHzSXp9B6wUYgPO3qV884cOpgoKB79pZBu6C4yCk6vzloEoaiccnTt2MoGinb5I4hR/LyTgNY1KMPaxq68K2l0e4MNB9nfQIY8/5cCmrb2/FR88BJU+/Fwo9GACWl0/XzT2HnXvUtsB46+PZ+qtXxS1nS85955HmyRUvJdalspE6tFt1p5dmBJk878gI57gfXyZEn3SoLHj1t1LNez8izD/uZ/PmYyWX686omjQjaDaHIglDlOj26vs5jIOTtg8fJnJYK1KWG1bWWNgzK//bB42ReUdutzwyEPOqEa2QB4ZvnR5261XSW4fTFFSMnSSOgjHD66q/HT6qjj5OuoKjyqBOukSccd4U87cgL5B2Dx9XcZzbIK0ZOqjj7sUG2qLqcPO48+fULpslBZ06RBZSq6mG3joGkda61T5LMJzublXLlyqr5H19FHT/Awk+aeOHrAzn5mIv5yVG/4KldhoTKUEmkOG8+lWcHJtc/cT07rf2UD7rtyDlHnENO1cm3qkjC6Ch9bzCvZ18uefp2smaRrYo50q0Xn4Jod/SpKWmxlZGLJVlCSR9+4rx/8fCeIzCFGmiGGFTXqN+qafNFPfrws3Hn+ublzq8a6UsCs3fdlwfuuxDdx55bAL02rE6UnxeCkgXOu917k0slD6dZjwtufjQl0dFL/Me6GwU1haVq7PjFZ1w/7VqOn/dk2fVCECxgfvedQ81Nf/7MX0h7dquC0uXA3824gV1Wf0y/FUv489BjKNLW26cE8q0+q1WNrv4AACAASURBVILqK6DqHZXbGs4Nv9yCxkYYJGCk0og774Tu9fXFD1u4qgeg+7SHefHWn2IoGrploCWItRqGMIbqDDYBZFq9dM7pPZBp/UawqMcuDPzkXVZnO3PH1CtD1T9SSnquWxUrWEwUrUF0Or8LSu6Exy946kuh8wSY17MvX2Q7l7f5DqxW80Z3sHqnP6xWFZHzLgwjP3ijTb5eeFVuurRQLKuN/xZJyeROsU3fyWS12nPrdqXhp58aKemBbdSCnDRPL5xviwhsVSPt4yNKswy+s/AZJr3+eOtvG8M1Bi3QKrDHqmWhtEWZSj95x2S0VhWTgWDWLkN4sfdAttnwBSe/9mgptGfM+gX97cDdjqKVfr/0BaGiSauNcGa70gflbbWmyas6ipTcvt9R/HW/o/nXuGPoFlmL5NiyGf+qVfS96GwUs1hmnHHX9zgTxs0YwtIaika/FUv4xieLOfOlB7AUFTVGvNB868WrOC4YqkE1Eabiwt02cWC5QjZCiYGf8OYT5LQ0yI0SuaFoXDL6VC6fdRtYPsxUgIhZaNyjxg16hktHn8q8nn055N0XOXPOA0hKkmdB1ZAShFBKG3YPFGmx4xefoZhm5DjxtkG1SDLGHUSNRR3JY7sOZcz7r5B13e2QlBa1jUzfP08/mqL8v4bRlDHbxmhOIRn1wWuM+uC1SEZbDwTRpkobU9XRXEYV3jnm104OI85YBhu0DD9+7TH+031n7vj3IC44dM/6EO1GUn375niq1vG/8oo0t+ocqktz29a2CK0m/Z2jM/TmURCqzIlKXXUcG95BZ06Rk8eVziiaUg0yp+ryuZ0GyLxPft76xNGz1ltX6X3yQonUVdsgCz7t5pd/i6rLQWdOkdcOnxBY/zi0W1Tqd4PKtEFOGTC6TR8URVu9csHHHtzRU7dE3FGwQeYUrSLPL9uT9KynPcdV1HiqVxlJz4+M1r6sRxu1aCm531n3ys/X56vjf1JKvpI6/t69Ucy25owGok1gEEfNUVB1/jrkiFjZyoD3OUWvsESRlPydZwIsC/zyLiDIaSnOO/Qs1jR0YU7v0hnFSzv2J20ZHPDR/FB3sW5zsCBag8p30xEXYWlT0i77lw+DXwxWP320ragM/Wg+e65aFps2P/qC9LsSwaN9v1UhDU+Y/xRZs1i6Fd2qMnRDl3aFn5X79h7Nqzv2DzUnlsAzvQchpJ3I3LHekAGPG0nPesLKCEvzQZftEuRaibg0Re25qzHzVsHXG29Uvf1gKBq9m1eyfG0uOnFCbNmqnu7dEXfeiTnxR+RtyJiFksc9WdkFRVWnMd9S8V5S2fmOPs77PkhfnGTwCEpM8PyDT2VO74GMW/Qs18y4EQuFTmY+klkbgNtArVr9bpCtvNckMQwSKvTztaLByHPLo7+NVT6UDuHCHG55IZCMfm9uYrrcel8HRy98hn/vsg8ZI9/mvdOWpqLyzNf3YfTSVzf7uUpe0VBF6dJYLfAbJ4S880uzS4B5bbUIGrthMSKSjHM3wuqdlPHrtsmyxm3ptXU24ZfR2LIl/lYoQiCkrHCA5EbWLHDc27N9O86grbTipz+UlHTG9Rqwl8/6Ey/eMpHfT/8DWbNIYwTTp7Xs9nBJFjaZwxAnFKCkZLMfF0H955aoSnmW/o2z2/Dm38muPMD0SxeFjG1y7fQ/+B7g59QUwrZrZvp+kmQcqdpQVIpCZb2eJa9q3LnvUeTVyohvSVGNlFxtORJC3YY7MAErxp0TQ1F5eud96kBd9ZCUwko6O/5V2c7tUs6WzfhXrYIf/xilkKeTVQxkFpLSZRk9QH3y+wN+wJS9xwDBjMcSCkUteuLEWfUF0GAZrT7tk0lgbhO8MBqSSB/VTGYJzNtuN1K2GVnWazvsWaVR3UYIzxO2yFebfzXfdDILvruNBqsYabUStwy/PMyQqV1QVMZO+iOXjj6NlG1iKhonv/oImRgX/eKi1v6Mk6egdJha9HgA9aYraGmKAYuamwdYispBS4Ndt1SjrgmDX34FLcUlY05j+E/+wrR+I0hrSruoejYL4xdCnC2EWCiEeFsIcZ8QIrmRcxwsW4ath8vAXobhh5+9cC8nLAi/qanKku9tP0ggr6UoqDpT9xxRti8OSlsr6sXsqqHFGcwFobLfJ/+JRcuQ1nT1RDXnF+3BrOqNOLQLQAqBQaU0LIGndtufL7KduWT2baQtg0YjT9o20W2ToqKRb423XAuji6MG9JPWJTC13wjyqjcMSrg9fpCaBkoqE8Wz63JUuO4na1YKh+70tUbwi3vGN69nX3o1raBrSxOWLbcMVY8QYgfgp8AQKWV/SsLZ+HYprHdvKARLMUEN702TkdHKiCAGV5ZATZOrR/yQq0adwrDJ93DtASdS8JmYQTTWmynF1bf6IQ7NfofZQWmT+nCpFUE7v1d67lH29R+3nLDDUG+6uPlVA2+d8nqaC8dO9k03+v25THrl4Qo/ToJSWENhW2zQMhQVNZEaLgzeehkIXtnB30xxVaeunH34ObEYraBkQBA2VhXbRrEtiorG+lS2tLgp8ZSi7n5NqjYMyi+oDAnM7bUX0+8+i7/ffzEv3vojrmh5m26N8Rw2JiPEx9SnPR9gB+BjoCultpwGjAn7phaXDc2/vDTxNfu4aZOkL7tTQMjJ486TZ4w7VxZjuCPwfl8PE832cANQz/yCyihW2T9xTfKKiiqvPvAH8tyxZyZ2r+HXDlGmtUHmkWYIvQYiliuHnKLKP3zzOFnwMRFt1tIyr2qx6uB2bV2NeaOT3opRd/d7P7PoWp8CijzhuCvkyJNujTSxDWqPWp+k89fKtI/LBtHKjDcphBBnAVcBOWCmlPIEnzSnAKcA7LTTToM//PDDqspas3Q52d12IRvhpdLdClGSrvC8SyoF5BUNKUQkTdXAQlBUNYRtkw6QuquhOQxh+dValtPmFgIpQI85XiUlXfaD/Q9mlzX/Zf/lC2OpHxwYKKQirnh569ai6qiUwlo6+vJq6u43zhxYQqBI/7MBSenmp58k70ZRKCUrtohgIU6e1dah1nlSjzy8+V1y8ClM7zeCYcvm8bsn/kDG5yZyGA3e30hIU1QdvL/bnTujzJoF++6boJSNEEK8LqUc4n2/OVQ9WwPfAXYGegKdhBAnetNJKW+TUg6RUg7pXoOvio+1Rm789g8it9CWUCj4xGeVQB6FnKrzizGTufmbx1JQS1vGvJbyjWcbyZoEJXmmDvCWZQmFw390Iyd979JI98bVluH9LezQPIhBxRU3HFWZhkRLIKSYisq1B/6Q7y58JhbTd5dVso5K3j9CCMaPvwohN0Zw80NULZxv/dpJ8WkDScky6spvT0JBtqmHW1XoPHrIeVQ94Fe/ag/H4+SdBL949m5evPVHAKWwpFq4GiWK7le33z22ui9OfhW/F42SyrrO2ByHu6OApVLKVVJKA5gKDGuvwnptnWXXFUtD00hKh7OKqHxvIrD0FEII+n6+jJNeeQRpS1JGkb8NGEvBY8kT5zBTtWzSPr5b3OWGDaSw3YkuLQ5590Ve2Hkfbt/3qLpMwjAG5g3aLiktPkHtIIFb9z0aU1QXkzQ2bJtz//03sj4SXZxJ6jDeZrXSsUDQzu+J3fdn57WfIGU4pXHr4deGQe165chJvLpTf98Yyd5v4x7+O32ZFNUYF8Rh6BaCK0dOCg0XGkVXg1Uk2+o/a4OeqdkdyjdWLGHKgDFccvApnHnE/2J4Yu1Ws1BJoEVLk9NSzLnwt+3ipK1C99PeDzAUWAg0UOqLvwJnhn1Ti45frlwpDc0//GCUfi3O39XonP1cBdSiB/Q+uVa3Biccd0W76t+DdNQFJVg/nlN1+ctRp8qcjw69mvZMQlvSp1lPywf2OijROcx6LR3sYrcd+6FF1eUFYyYHhrCspmwn36B+zrX+Vo86GAS77WhRdZlrdcV8xrhzZYuqy2Y1JQsIaYrwuRT0bFD18hlb2NlM3HYyELKgaOU22aDqskWogS7Uw8apDXLmLkPkyJNulbtdOH3LcNkgpZwLPAi8ASygtOu4rd0KXLYMOxXvYooj5eXVkullTo1nCpoUXim53pAI+q1YwqIefTBCJGtZVd4bEVQLP0+UDlRpc/GsP5P2uTNRq7mcG359k7S+KcvkO+88n0hCb/Sx23f/7oX0oSsunU46Qclh3OWzb+N3B5xIUal0+x23DtL1OPkGWatpllV2v1zNWHKXpwIIwZN99mvjYtlQNP447HiGnX43c3oP5JoZN5K1DDpZRVJIX9VXHJTuyFiJdlVBcFSRKdssnxmkbavspVN6nqiyBDDqg9d44u6zOHzRc1uOHb+U8lIp5R5Syv5Syv+RUrabwnHNNttDIX72AhDSKulqPQw66lAmbv7em5xJh65XxeBFxipy50NXMObdOdyzz2GB6iLHt0x7Q7LRZluzLVL4t2U1/uf9ykryPux3zbbq6m4iCH5H/EGMwv2bt61SlsEFz/wFIW2MAOEiikknUQWp2GUfS9X2m7u8lLQ5aMlr2Kq28Z1tcsZL/wCg34ol1OpUXVLy01+IEOpqhSJtdOK3p9+CkDGL/HraH9ol9OKWfXMX+KQpD3bbZo1iAsXWwBKXH3TyJmGMtcB9COiWANOWydUzb+bHrz8e+G01jDbuIWlFWREHnknhJyWHsehqDxdrtSKJGj8FVfd12xvGMMLep5Do0iYla78V7M07CmGLVVB6L3RsXyl+wrwnuOufl9FQxaG0BJr1DHktxbUHnMiECVcTHiamdiTpt6D00Brr4aMP60OUC1s84+/y2XLyetuT+6hBrEmL5V16sGi7PmzQ/S8Ve7dvXs+MYaiVmXixQc/w532PJqf6HzQHDcI4eYfR4Sd1+qWLYqDuBSsqjQO/7XHUQta+CrZoGv2gSJnId3yt7CrpYha240iSt6SkyoszR7wHrhmzyBkv3O97UStOXYpC5ZLRpzHsJ39hTbYz991/IcKyK+ZwNf1XD0S1xTvprnUvc4tn/E3b9Yp1cu9IjDktxeUHnUyvphU0p7KBgTG8UpnbN0xSJBnMQdL0P74xBiXEVXM1dITBJJ5OvtodQj3SJFHv1GOSB+XhZnje8h/pNyJW3kGLY1wJuxY9vBfucZ8kX4XKeRNH510QKlKpnlUp0ubZPvsy9t05/HrmLWQskxR2BS2bitE7CNulSkpuXs479CyMrvWPwbXFM/7td92RCw4/y9f3hxsO875/71Fc+vTt/P3+i5l+91ncP2AMOS1FXsSPzlnNAKpmm5xr9f9z+cEn80G3Hfn7IH99fj3gzVej/aMc1QNRfV7tt0ngZnhu2MDHnbeNLFdChSMyN+q1cPrREfVdi56uEACSqjmiIJSSd91qIRWVr+XWcdns2+rO3GsdI0G71Jyqc/ZhZzOt3wh6dqm/r54t2x9/K4RdCt1XECppaVFQU6R9vHUCTHzziZI1Q2uoxvHzZ3L4xBs4Ze5DHL9gVvvSGfKbdxKaCBTbwlB1Lp91G9/4ZDFHL3q26rL9JrkEiihIpXRbNOXa/fhJa7WU1V75FFtv4MYtL6m6y424KjUHCnDOnPsiD+sBZADjSzJmkqp4gnZQ7vcNRvvd/nW+fajfSAp6iolvTC+/97ZZ2N+monLcWzMxFJV0wtvyQSqtauFuw7B8spbB9U9cT0ZINhSH11CiP7Z4if/T9z/mN60mYI7jMCUkYIq3MwxFo7GY4z/b9K4LPRJoUcJ3H2HfOo+GJG1bJc+KrUHSo66fhyHou6n9RzJhwtXYPlvtpAd+7vOQsPRxEEeNI4B0K9NPsqMKU0UEpXPKTLJr8xtvQQwhjUy8I6zm8lUYktCbJE83vHUUwNGLnmXC/Kdiq2S8v2fNAqe++jCdPAFxwuDQ4YiH9VIHxd1lClqteqbf0GHVUw26fLa8QlryC4wRhJRtstdnS/jF83fXTIsE7tpnHDcPOx7DY2sdR3IL00fWQwrxwzGLnuXeBy5mar+RVUv4XrrtkOAZcesRNAmjDh6TMmlbxLP/iKMrrlf7xcmrWc/ywN6jayqzXvBbQMPo8vrTtYRS0yIWdyH3+y5NpXVUvc+G/PJ1oGVSdP380zqWUsIWz/if/SRf4QwtauV2MwXFNLl81p9JW7U7qLWAKfscxhkv/4OUbbb7gVKYxUIcCEoXUbJmkaPfeY4rR07CaI0lEIdxB6Vpz9iySaTCOOmCHKJVg2r17NXkpdkmM/YYjrGZpriX7oKqly97QTD9grahQ6FUl1SdgsTUe1F2vg2aE0lVRd7flFyOLcVXzybD6uYCj77wn9CgyhIwhKCoaKxLNbQ6XtvIFFLYkVZBsS0rFJW+q5aFBt+uJ+JKil5m7VcfQ9F4t3tvXz88QfXPBdwajsNw/VREcdQ7YX/HKSspc652UfX7zu/vJBppdz1UKenVtAJLix5r1bZtovMQIfjpkef7Og4MYphltaZtRQoLcZlu0G+17KDDdhRhu4U4sC2bNc31v9+6RTP+5WtzpKVVEVTZjYKi8tMjz+eQSTdx4vgrOfuwsxNbq8QdJJptcf3j18YKcVcLMwqDW98dlsZbpt66a/JTfAQxcpXSwVoc+E0Yd77V7hGS6tzj5icp7eB+d8CJ/GLMZEySuS4IU9c5+dw3YDS/HTmpTd5x+y5lm1w6+3YuP+jkckStIMS5m+yl1URw5+BxFFvVds6TV1Ryqs5zvQe2ef9Un335bKtupH2EKL92d/d/VACUeCF/ostM+m0tuwQvgvIqpNJ8vnBxHUsqYYtm/L22zrLD559UvG8jTUjJNTNuZPrdZ7HTWn9dWhzdYJC+2ZsmY1sorXFowyZzvQ6T6gEJ3Lnvd9hj5VJSVrw7ERKwFAUpZSmwd0Tg9ajtb9xbxoaisj6VrTgXqab/3O/di7DzKMC/+g7nlZ36YyrBHkmTwsnnuAWzufiZu8qCy7933Jubh34Xw8Vsw+pgKiqLtuvD2Yed7VuOpMT0rz/gxAqvklHQkEyY/xQIBYOSxVyJiFI7fGvZW23aaty7L/LAvT+v0N/7IWkb2kL4Cgb16IukO5soFXJQ/oHCk2WyzV67J6AgHrZoxt+tMc23Tjis4n0baULabFXMkTWLXDvjhrpy2iDp2hsIvD0ZfNii4vd7EB0/fuURLnjunlh0Ovk3WAYpaZcC2duVclk9JSanXM22eHqXwb6/V6svD8O1j/+eGXeeEepmu1q43VwI4ICPF3DK3KnlG6xRdKfMIp3zwRYhzoK6cNtdEh+eCqCTkSdlm6TYGGozYxlkLKNily0o+eKJ5y4xGXRp15WReeeDRfKxGjXv4uaxYOIZdN25V8LS42Tu47Lzy/bU5JZZSvnKuO/HcrVqg7xx6LHS8AlXV+1TqEN4xVpC34WF8Uv6bIowdEnaKkmb2yDzIe6iN3eb1KPe3nQGyKZUgyyI4DFoUwo5OXvnQXXt36h5Vk8X3LXQFTTf4tTFiqhLnPrlIvpG/ulPNfE+vixumTcHtr79T1w/bAJRVu4COGnu1IrrPrKGsv0k3bDyg3S/YbsDL30SyKsaeVXnorGTQ10zJ0E9JXYJkf3hTe9+wr40hYrpkWBNReWaA39ITtUxqI12N74MqjiorI+jG+9cbCEVEvheUBqjI5e+2a70ecuMoxqtBRJo1tKR6rANeoaZffZLdKbmThelTiwQbLoMJY1DkPJUAJx9NqxaFUFRcnwlGH/qrjs4e859ZReyEDzIMtgV/vLrcRDkB2cAVYug73OKxoVjz2DY6Xdz/8BDOeeIcygKtYJ5BpUdl6ZqLWAkJUbsvcYmgSf77EfOxz1Gkm1y1ippkg0hKKha60Gs4ILn/oom7cjDQoeWWhE3Dwk83XtQrPRmyFlJLTQnVUF4/3b70JfAc70HkVP1UFcnDrPcoKYpInzvTCStk1P+1D1H8HD/kW3iBfjllTYLHLjszVChKkwQi8KDe48qxYsOKF9FhrpPMG0Jy5bFKCkZtnzGv2oVO1z280ALCi8E+FoetAdqlRYDdwiKwrN99mVNQxcApvUbwSE//iOGorX5xgq4SNVeunh32c5OyD0pLQQjPpxHRlptbL798vDCPbEEJb2vJiWp1vsXGdskJf3PGrz5OHkE/R53wU7isXX4x2+H/i6BoqLxy9GnVYT7lMArPfeIQdHG9GGLRxStfm2j0zbW79DlC8uxn42Q84OCUPntyInYmo4qa78zIQBDKBy+eA4nzvtXm3gBvnkLJTJcZS049u1ZaATHLIhaQNRCnrWi/icjWzzj/2LRYloCDDTDJmXUBEgirceRsusBR8Vz3qFnlZm+g+8umFVxHyGvBQ8oLyOthaawOrc5aEeSNUuX5N27s6i83dtuv7zrjX/2P4h/9j8osi+93iiD4ByKRtF61hH/y/0DD+XcQ88ir6XYoKYoCoUpA8YwYMWSRHV9apd9A3dVQWoLN7zjw1u24+pkXaaRfIBrc4C0tFjV8LVAq6KoRd4PurRjuy+xhIJu+dsaxZmrEmhpjcvsV57bv1U1yKs6qz5bXVMeftjiGf9fPy1Zevghjk49DHHtOJy8iooaKv3UClNRmTD+13y09fZ0bWkCoGtLE+c9ezenz32ooj6abXFf67X+JAuSBPI+tyOC8mkvPXiSLXeSPKN+H/fOv/nTN78XekchSJiI8y4I6zOdWmkoHdCptoUubY5d+AyZBM7HBDDiw3lcN3wCBv6moXHoNEPYR8oo0JzKsrxLj9B7K3lVZ68V79NYTBZesF42VGnLQLGT+UFy4Kgs/YKdBQkj7t9jvReiXcw5K057v4xPtVY9n6/Py94XTJOTx50n8zEDQyexmggKGB70zXotLYsey5JqLRyCrBMMhGxRdZlXVDll71Eyp6UCLQ4MhJw87jx5y75Hxw4kb4N8bPf95fTd9q9re7bX9+3xNKUa5AnHXSGLIdZfca084tavoGhy0JlT5KAzp8iWgIDqScdLXlFlQdFip48zPty/FRVVXjFykjQIDoheUDTZouqx83Z+v9lnzFYzZqLSRllDVVuOCaE8yQky/9Qvr6+K9zngq2jVs/CTdUhgTu+B/OzwczZeMqkDBPgGDA+TGDuZhYrdR7VSa9AuRUOStQzStsWEBbPImP7up5201027jpNefSRSTSJbH0PRGL3kNUa991Ismuuhs5Ux0nnTxPmmWjQUW9hlzX+xEl56Sir9Ob8VFY0bho8HoFfTikCXH2H5+PVv2rZiBSkK+j7snMw5x3FfQPODYltYMealdD0AJ73+OFP7jYg8qzCFgkn7jAfvXI4qQ1I608irOmeNO49Df3Qjfxt4iK+VWYue4SfHXkL3kyfWj2AXtnB//JJxi57lmhk3lmJXCoEpNwYQqZUpFdQUulUs63Kj0B4qDxmSb5zy9ARRuwSUg4+3J2OtBkG64Gosj6K+UYHzn/lLoG/3OAzAQVRZJgq2EJw29yF+OucBrjnwf8ruM7yIa6lU7TiUQF5NkQmIZeHXdlELt4qkk9nWXXKQoOJGyjY54t0X2aBnaHS5W/bSYAgFRVHQAm6cx1HtVfNdUHohBOPHX0WvdSuZfvdZIP1dUmi2yaedu9Nr6/oHYQGo2AJ8GZ9qVT2rP/i4Yltsg8yrWl1UCO15CWVTXnKph8pmU7RFkm+KQpF5VavIZ4OW3mTqo1raJEgVY4P8V5+hdVOzRalU3H/nFU2ee8iZMudRzThpjRrpSNpWG7S0zKn+qqoyzSEXpKqhIaq9iojINAWhBKrYnDQ5VZfFdEbKe++tivc54Kuo6un6+aeIVGUAcsfFsmx94sJJ7/idqefBYhAdfrKdxJ/2JHWpJyQl00U/M1AjwkdPXBQ96oCoPDVpM3334W3a6r69R3PKMReRV/WKvNqj7Vr0NE/stn8bV9ZJyvnXbsMqXCkIYMQHr8X63umPautXoeqzba546k9M7fftin4VlHZCfuXE3Q0nnUspy2Bqv2+H1i/s8lo9aPB+ryErLg9606RkuMdfx8pLL+Thxz/uuMCVGL17++rhHSTt6HJ6uzoTLQmx9Y1lPWnIb0lolwQvGN7fg+gJgsB/IJlC4Y/f/B4b1Hhmo2GY03tAGxpf3LF/qNfJvKJx2OIX27TV0e88x6IefXjExSwk8MKOe1MMseUOo0+2luVXjwajwBs77MmPjr2U079zAacfeQEFn0XHDwI45L2XyPiZGvqMaQltLivFgbNgxx2PGWyylsGEBU/5Wqe1pyDkVy9V2hz79tOB5cqA97XQE0s1JEvWSi0+N4fbpAsppwxdb5cLXBVbgC/jU5OvnnvvlYVURjbH3OLH3e75WTkYQgl871g4HHXCNb7b9yg/PXHofm6nAbLYSkNcfyPexxSiLvQ4aY864ZqqrFCiymzRUnLkSbfKm4YeK4sBdLYobVUSTakGedPQY2uqk5N+nZ6VLVqqZDmlarJF0QLba72ekS1aSt61z+Gy4EpXEKqcMmC0zAf4c7JBFnwsYsIsQXIRqoag72rpn3r1aZtxSPAYjEt7HCuqJGM5bt0KKDKnqHJDTEvC0HJ0XcqVK6tmf3wVVT0Aq8cdw4Gn38Wpx1wUyzdMHKklKM0D/Uf5SreqtJF2KTDGvQ9cjLStim1yUJ5JpKh9Pl2MJuN5bgyCE3GqHvQUFA1L1bj2gBMpKhq5Vsk4DmTrE1SmoWgc+p8XOen1x3wPxwSQsdsqynTL4OTXoi2YwmjKayl+MXYyJ0y4isMn3sDR7zxPxjLJBkRUE0CjkSdrFpn4xvS2kdeE4NoRE/nxsZeS09IV5ZV2fJUtFnbgmCH57dcw6bMY4WvG3U9R76IgXU/cy29J4dc/fnTEgZteN3RsMrZFQw0xsMvwuyRQB2zxjH/52hzrGrdm69w6hOI/iP06LwhB6QTw/QUz6eTZmjsDN4XNxDemkzWLvnE8o8oIKteh3WEw9TCfjFt2FBQk4998goufuQvdNskkcIVREAotPszQneSa3QAAIABJREFUQcoy+OlL9weaqwIUVY2CUFivZ8lpKW4adnzNDuse6/stZu4+jPnb705jMVdTNDXdNum3Ygkb9AypkAtYScZnLfArQygiUs0Xxwon7mKblNlH0VbPdvPmFZdWCeQjFlAnXcUNZlXt8NVTDXptnaVz81p+P/0PgX5akur5w35rn/W5EpISY/tn/4PYEHIl3kG1HimTLlDu3x7b40AmLJjVZoJ4DxyDFrk79jsmcKpISgEqooLCpCyTlLTJGnmuPeBE1mY708msPoydAL638Gnm3PxDzn32bvquXEq6hvwA/ue1aUydcj6q9D83qseYitvvfuVEhS2tFWHjICxtMebusRojjlp3un6LXDoihpwE/jLocHTv3YBcDhobY5YcH1u4HX8pGMt1u8mKU3RHSnaQU3V0y4jlubG9Ebf8oqJx1MLnQm3xnQGvU5ImvAMrTpnVtIcFvL/NjqFpbKH4MjxD0VjYow93DjmSU159BN0n0L2fCsSB07fONxqSi5+5C0OoiQ/EoW39Syokk8lzH6worxqM+eCVdh1vtdBmINrdYWE1goWB4JF+B/K9t5/2TR+0gMVXNMajyVtOWFvH6YOm7FYUhULaNScModCycg1f2zMZnVHY4iV+gGG7blPxztsRmm1BwE3MODrLKOnXIH78WT8E6Y7DzNXcDFAAqdYdj99TT0jg3oGHhDqBE4Diw/Qlpductzz2Wya//CC6j1RfDb2CZJfVoORq2IwIGF+rJFgPhI29asqQlG6Y3jP4iNB0psvVtxuWUMpnA94nCUyhcPvg75QX7I1qU8m4Rc9XqNnC+kIqbVldtTS56fB+X2t/Lum6Q4VTN13avNdSY8Y+2CyMXwjxNSHEg0KI/wgh3hFC7N+e5c3Qtgt06lQ6wBLo0g512etlUXH0ejmhlfV2lqKV7DNC0sfJNw4kUBQKecXPGHTT4HtvzeLSp+9I9I2ktFNwu7GtZmGqh2pEAlarL38n2PnmhGMG7MewatVle781hODm/Y/jf96YHvqdJq2K/pHAP/sfjKWqFFyCTtJ+dOr4ozem+S7YlqJyw/Dx5XgLQd+XxpTwvQ8BpZ1+tQiqU5xFxe+366f/oeJdXkuhbggOn1ktNpfEfwPwLynlHsA3gHfaq6DVzQUue2xh6HGq7vOrd2IFNZTTwU7giWYtTU7VuXLkJLTWWF4p2yJjmxVqDa/lSlDe1Uzqd7b5eoVVS1zE3c34tZmzy8hKfyuXMAgIcKBdiTjtEsUQgr6xWmkpxQy2UCMO42tBXB20o4IMatOogCNRyKGUx+6De4/i7Dn3kZbeWHQlRJVx7NtPk7UMMnblwhAXgtaLTgE7Wt0y+GDrHbh56PdC8xCUBIkg9xrZBF5N4yKOwOK3+wuqa+/Be9WLtDI2uY5fCNEFOBCYCCClLALBfltrxMJPmuix9rPWqVs5XGu1HnDeD12+kO+Pv4qd137CvJ596bluFVrAgZ372yTbdK8OMWibKYABq5bGqpufXrJaJlcPyTwuzY/1/Raj3n+lwoqqmjLdMBDYQkFzSZntfeYTJ/+oPvIKL0l1+5qQ/HH/7zGj73Bm3vGTqs98BLRpu3rDGfMpaXPLY78NpCdowfLOky8rSsYbOgt/dR1D2iHY+uY43N0ZWAX8RQjxDeB14Cwp5QZ3IiHEKcApADvttFPVha3LGaiWiR5xqu5Gtdvm++6/kKKaQrdN7hx8ZKxvkkrEtXxfLaIO4NqDhqh8x773csXBo983RaFiaDopywiMy+CGjsTvZqwffV5Uu9DWki5J2WHf6lJyxpwH2P6LlbF3XUlRj7FSq1DSHuO1nmpaB0VV48hJNzH5kPBzlmqxOVQ9GrAPcKuUchCwAfi5N5GU8jYp5RAp5ZDu3btXXVjnbIqd134SO73XIiQuMmaRjGXSudhC1izy41cfqfDZUasu1k2jhAq9pbesOEhazzYmdUKl0t6mbdq4eSahS1CyrPFKuH6HfTZw0/7HYRFujx6nXG8675Y+rnmhgzCVyeZgkDYKx789q13VWl8VVKtyc769cdh43u3ai/MenM/q5tpMhv2wORj/cmC5lHJu698PUloI2gV79ezMx117RqZzDnmTQgIFVas4JMrYJrbPoWCxDk7LHGYTdmbQXmhj1aCI0BCEQQudc/i8LtVAPqalU5T+OuguQ1panDXn/kRmiYGWIRE0eO9pxumL6rw+hSNOuX5pVGlW+BJKWu7mOgRPUrYVkjbu+ZYX3jETJpTFwYy+wwGwbcnytcmik8XBJmf8UsrPgI+FEH1bXx0MLGqv8ro1pjn4G70ib2wWhcpFh5wRO19JyTnX3wcewvjxv0bzOG4TlJiOd0CkbCs0bFwSSwCv9Pjc1wdSrdV1NZM2ZZlMeHt2aJogJnrlwSdz4vgrOXvcubHLC9uJaZbBHUOO9G2jbMjtXkg2uYNoEFTeLXAf3AbtamqdgEkZlVct4d69/e7AHwbq5+PuZqo9zE0y/mpZXJw2j3NWFrdsv7M4NeBgPC6OWVCaV4Yt6ZSqv/Jtc1n1nAlMEULMBwYCv26vglY3F7jzY4mlBFdVAo/uNYIu+WCzKW+HF1QNkHxn0fPce9+Fvjbpft9GWa6ESZtR6o8RH85D86G1vSSxWib6C70HMn/73flarnZTNUmpTX/yytRyDNlCjCvy7QG/hSesjeql/vOWlwSGojJtj29x3r/vQbOsNotVrSqLJO/jImjhrVa96ce4W1SdnKLxeN/hFFS9YsGMyj9JH/iNmclzH2T8vBnoimBDsf6H5Zvl5q6Uch4wZFOUtfCTdazKduG53gMZu+TVwA5xbgEGWQOU7Ms3wvHpnym2lNMEIe4gCGPutein4xzOVrvgVANLKHyR7UzXliYun/XnuuiznRvJqVblSSqmEqUWJlQvK5FNpfsOklZTtsXR7zwXS9qtx4FzLWPcD+0xRhUpSdsmo95/BWwbE38X6UkRRKvfWLp81m08u+e32iUK1xbvsmFdzmDCmzNCmX4sKwtFpSihoKdImwaWEG1sgGsZeGGSUT3yjdopRH1fa/neclr0DL2aVtA531x1zFd3GbVYOyVldnHVApvKeiSJaqLaceCgIBQUKZFIUjG/iYsgZhj2uzeN+13c8RL0t+NQsKF1jlczFyQlk8yUZZBXdDK24TuWgmgtqhpn76rTrTHYWWG12OJdNix6632ueOrWmgapI1EKAbpRoKhoZGqIt+o30Pyk7v8PVhBRlil+9UrZJsu79GCr/AafrzYP4k5sA4FJtIVQGONJGjDFDSfoi5NHFGpRx3mRljZCiFhMvx76+ii6nTZs7wUo7H0UFMvkuuET+KuP+4uo+qWNIqMPG1plyRF0tUuuXxKsbi6waMbzkRep3AgbiJptkZE2jWY+UDeblCG0J3OPO+GrMSksHW6rTNl7TOz8nN/e6boj/VYsYcjy2s70JTBt9/3JaymKoj4hHr2whEJO0cq3Wv9wwAkUQ3wQxcGM3fenqKgUROkcwuuOwQ33ewFlv/9+jgqqlUqj+sstmeoxDy3j7iRrFciS7vbqMd+SGAPoSM558T4mvv544nIURdC1HaR92MJVPcvX5lBjBjKoRi3idyjUnoy83vlLSt7/FCkDPReGlTdhwtV82HUHjnt7NqmYN10FMHDlEv72j0si07rp9KZ1GObY915BkZbvgXlYn3pr6+e9UQJXf3siD/c/mF5NK2hOZem7chmaj3oqqm+cnG1g3OI5bdKqVNbNoTtKAq31jOLvAw9lYY9duHT27dhC0GAUKlRJm+OQPC7itHlQ/xdUjaZUJ7bNNbXrGZdbdZQEJoKnHn+JQyfFuwyaBFs04++1dZaFPfpgKBqpkIZPcmIfhqTfx9U3RuUfd3D6pbMkKEJAiAO5IOy9YglHvlNyC53U6iEufQ5s2lpDlaQpCLtlG8W02jJfiSUEiqsdnuyzHw/3P5g1DV0YvuxNfj/t+rLDMD/J3JunX3l+C5T3mw1qChXIxnBH4Sd8mEKgyWj/QgVF5foDTmRNQxeElFwx80+xyvBDVP0l0Kxn0WwTVUrf+SgpxarN1CNyVQQ9RVXjNwdO5NJn7ki8Q9hUKti0ZfCHlz9lv+MKddfzb9Gqnm6NaS784QGce8TZ5NVgD+5JdaD1koAsoZBTddbrWfKKypN99iOnaGxQ0+SFUlO0qDCVgXuSZrEDg4BEYfBHC5j4xvS66ZHDpNtqB2pcna2g5Gb4/EPO5J/9DyKvaAz7+G1evPVHjH9zBtc8cUPZBbbf47UNd8ORLuMiZZsoAVZJUWNPAJqUkfc5JDB1r4NY09CFri1NXPr0HehUb3teUHXuHDwOI+CCoongJ0f/guGn3805h5/t67YZ8D078wpDQWqxsO/c35sIhGVzwfP3ROQQjThnerVg7H9eaJcLXFu0xA8wfNdtOLvfCN7edheevGNyG5891WzjqtWj+pYjJZeMOo2sVeSF3gP5ItuZfiuWALCoRx/GvDuHq2feHLmdjStl1/sQ7PDFcwJ/b2+1V3vAUDWyZpFxi54na5tkiyX2ednsP2PXKCMpVvzFVZc2OVv479CoVA15UToDCIcAvrvoGd7suQdZq1izZVnaMvjh69MC88mnMqzLNLKmoQtbFVpqWmS89BGTRgcqrarNhA7+gpBXdWyhkDULFYtUULAhdxoI3ln89N9T+GjdFcDX6kKrgy2e8S/8ZB2WhMZijlwqg14Mj2oQNdDrycxUJNc8eRNQ0jFLRaGg6ui2xTUH/g+914T7GKrXzqPaOoUxIAvRru6Mq0HUoWInI89Fz9xZoYYwFL3mEItBErwfBNAQoMLKqynuHziWH7/+eKxxGjaeU5ZZHn9R+Th5+b13/g6LiqbZFsu79KBrSxOXzb4tclwElRXnvCYsn7jfW0BeSdHJjnbk8szOg3mr5+6c9cK9ZD1j55E9D+CYRZV3JJxF4V+7DuWw914KzFtFIp5/FvbaOYKKZNiiVT0Aiz9bB0BzKkvGM3nj6O/iMNdqtpzOe+fRkOi2RaORJ20ZXPzMXZz4f+ydeZgVxfX3P9XLvTOAjICIIiiKcUFURBQFDeICgqDGEAMuiWsSJcQlojGLicbEBI1xiZr4c43r65ooivsucUEhqLiiKIgiAg7MzF36dtf7R88d7vTtpbpvX/Sd9/d9nn5gbledOl3LqapTp8757yORnU5VlbEhUN5+FzWDm3efREGRCxWrkrL5YpjKLorOpaOm8fOJp3Ph2BPIG5mq9IL1F/MqkbXyXLfHYYkn2iDdfhI02EVEDEaiJoc4arqcMHhim90piXi3oiXw2LZ7snFuLQOaVyit9b18hal2wlRsTw8aHhn5riB0Sp4oXwiN7o7aGcvBH7zEOc/9E9MpVU00k955oVNAmsp3tpTM2X5UZBl9W9dEpokNKWXkA/xZ5bd6PbvvvrtMit8/+Jb86eSzZEHo0gHpuMeYkY8DMq/psqSYVoVmrY+3nKByw/hxQFp14jkP8tadx3Wqa5Vyovj1/u2XPqosB+SVe02RW50zW251zmx59JEXyHVGVrnec0KTRffot25tq1JfcfpwLXnSpFP5bY9sMyK03ybp00n5dceCkEWE8thSKTPON+V9yq6it2hRYvkHzPOTqaor/oN8fpuQ0txTV0zaTK86mFPFRfsdzx/3O+EbY87mt71OQkOHmg6O/SABdIOj3nis6hBUhadKOpamsy7TSNGHxzCaQWqJ8rvpL93D1Plz6N3WDFR7Nw3jr0E6GCn60pSeR7W+kp7VxDEmLPPk/U21bL+6Lz/jPpwXmC+sTybt61H5dKRvBD4vVGWA3xgNo50NKVsCbx48BXZMOdI6ETp+IcQpwKnANu0O1crYCHgxdW7qgGFyLTldhwQh1n771HWUNJ2HthvFIR7b60p4BY5qJy13pnqoZcL4EICQkoJukPFRbYTRrKThpZn1CYyeBMJxaCgWsBXvYATBWwcCuPDxa7jg8b9TNEx06Vq/+NnRh6kB47RxEOqligviLeowrzKfigo0jI7X9DaKjsS9T4IQ6hI2IcL6sPd90Lu02k6Fzl+aduUvLRvenPN2YDLwQPu/5Wd3KeUxqXJSLwwaRGOEUtTvbXmmNh2bQ96bW+VrPSn8Vnne98V2HWo5WHsSRHFb0A3uGbJ/KI9+NOt9dtBx3oFDQ0AM0jIq69KPjl9eXTpkpHuWknFKvgfUSc9VgnhMChUh5AdLaJ3qpR7tFlR+UFyEMDolAKGRcdQXIbUgqj783sfZwUYhTp8oSclby5tTKLUzQgW/lLJZSrlESjlNSvkxkMPlu4cQInk8xA2JJ57AdpyqrXWY0PAiymLBmzZpOglcP3wyU469hD/sdwKa5l40ScJvFLrbRY5647HA1d2GEhRR7/zSegVaHH7T/C6Jexej7HKhVlp+NILoloOJeN9bQuNnk2dyxV5HKgcWclfc8VR/IqD8rFWIZTHiji2UnfXViriTehobED/VlwosobOo3+AYOdSh1EZCiMlCiPeBj4BngSXAnNS5SRsrVyJPPBHdsqqERJDVQL31+WFCytIMbh8+kWVN/Tj7uX+StUsdViYlTef8sSdRVIxY5QeVVWDSiWVDnYMkGQJp8eZHR5NOoMsFCBaQlf2thODCsSdUhbEM+9aSbvD3Pb7TaQKQuHcALp/9F3720l1U2y0F47o9qoPYxO0LEne3lnSX9nVgQ435uJDAn/f7Ic3dm9ipf8+0WVKenC8E9gLek1JujRs166XUuUkbS5bgKAjKykmg7H2xHp0hatVrI3joptP448OXVa2ADMdmRc8+XDF6mjJvcXcKUbTC3kWtlCzN6HCk5uXHO/mWhNbRBmm0Q9jhmZ+1fFidBe04ospzqP5mS9M5e/wMxp10Na8OHMrNwycprfol8NyWw/jJq/djUL2QKatM4qxeT3713+RjOp+Le5Dpl6cMG0FRi1ZtRrVNUPog2AjyWvgJSNx+mJbsOPu5W7i58cO6uGWuMvPxe2g3CQL+C2jl/6vkTeNJbM75xRfSyaqZ7JWf5kw3+ZNDz5FFRCJzLkntJmlB5opFhDxn3HTZppuB+RyQrbop2zRDWjG+oUz/huGTZEEzOpnglYQWy+RSgizR2Yzvnh3HyCtHTpF5TZdtuumbJ6eb8jcH/liOPemawG/0fmvSNirTaNUzidoojac5003+beR3ZU435DqjQRY03dccMC/0VL5dxVRU9be066kUQbf8zorBf1SfLYG8YOwJ8qlBu6X6PanWTWOjlF98kUz+SSmp0ZzzKyFED+A53JCJlwPfHGfqQejbF3H55bFm4AarwLubDuKv+x5NTjNoC7kwlASy4omD8jnDb5+8lku+fWzEgXSJ8w46hbiOio88+s88sOO3od31bvmJCivpt7Is+64pP0e8/SzTX76HrGPTGOKE68PeW9B/7UqsiKDfUStLvy/3+83rCK0WbWrcNm0oFTn15XtpsEv0KOXJONX7D0szcDxhQ5OcwUj8dzeVUHFtLYE52450rXBSRGV/8UP5nXf/XjmegtRtYWX++ukb2G/J/A1yqOyHSD5NE5YsqUPBPrOB9wG649a5AfwQ+BnQRyVvGk8tF7i+XJeXdw3dX3kWtkDmjIxsznSTOSPTsUpNc3VTy+ogLzT5mwN/LNeZDaE0bxw2MdYO47ZdDpI3DD8k9iW3NC++WAjZnOkm87oZurKrB29xV7F+tL27nLC8RU2XReF/IaySxg3DJ8nbdj6w5lV20sth3r8tROzLkEnqM266OBen6vGUyysKXT6w3d6ple2YZl1W/FU/fBOfWgT/gws+lWNPuiZxx2ozMvKS0dPq1omS3MLMC00WNCM0TZvwVxv4ddZ7dhyjVEeqwiPNgVSvOlQRFA7IQgD9WgVLEaF0c7joo/7x8lDPurIQsk0zZIuelTmhyUKdbi/H4b+omH5DCv5OfUY35BubDIqtHvV7SkZ9BH/ofk0IsU4IsdbnWSeEWJv+/iN9fNmS58M+A7mp/eCs/ATBu+WzNINPm/qRj1A9bCgI3BB4plMKdIErgMZ2G3jp+d2bTgAT3/8Ph77xVGTZ4Vb1/mUkxYa2mQ4qz6TaLp6AtHFgIOmu4PjNcMJrPaebWEKvUuOofncULM3g+j0Ox8DB0c1O3m39kLTcMOsnL1Q9S1pUt1u9IXCd3+305ZIqKy1vOhW06SZfLXovFd4qEWXHv5GUsqfPs5GUMn0bozpgaP8mAM4/6Cf8ddS0jm6rMgmAqy9f0H9733cSdTphSJJf4Opsr9prSiff5n7poqA7Dqe9fHckj0SY6nU1dFjJSH+3Yqp1EdQu3ok5LgRuMPBMQASyoDxx0OBY/OjVf5G1LbqVVN3uJUOQJVJl345zb6Ps+lklbVg7JB2fVFgq+eVXoWvYJd7O9opZejS6vHdO09DRBUybP4cz5t7RYXMd1oEksM5sJGdkmDnhNACeHjyiqpEq8ycdEN67BHHgIHhsu1EcfNJV3DxsYkIO3MktaoCUhJaah8k0EVRnlYMq6YSYhJe8x2S1KLTIQ9UgWg5C+RA1yqS0lsVJnMtVcSZDFWGbdHxJqOnOSyWKQuORwSOV+lolLN3ghCnnceyRF3D7LuP80wSYkkog3y5/rN59kzEegi7vj39Ar0Y2yan5AAe3wp/ceneu2PdoljX1Y8aLt3PlgxeH5klDiCSh0eiUGPPBq/z05XsS+bWJIwhKuoFRinZTWzlYw36XqAUVCSuHkLxhNMMuoNXSlkXd4IxDft6xlF/X0B2Aa+6/iI0s/yhK4fxLXtxqGPstmR+ZPgpRl6qC0sTd1ahYW/m1u7d/xPnOoLL9rKSCEHbXQ5MwpqINVKHjsKjfYFZ3c7UORy18rKpM76RaXjzcvetB3DTiMD7eZCDnf40XuP6fRZ8eWQ6TK7FjmJ/ts2QBW675jI1zaxOFFowSqKqrL5U0Z869o30brnpBfz3ifJcRYxBVrnjzmk4eAjTD4SVHfX+r3vliS9IVrcS9uFcU1eEAvWqGMBoPbT+ayx66lFlzruDSh//Kxrl1LOo3GDMg0pOKPrtsapjWmUcQaqWtyp+fSWYtKE8YSS6T+dHyawMTh2528PgKUlH9zx6HA9C7rZlDFz3jm7eEoKCbrM10I6+bvLDlrghN4/BFz/HQTadxqfP213eB6+t+arHqWXvDzbKt/UJSnNP1nJGRvznwx6lcevH720KEWmyoPKplp2HxYUFsi44CQhaFqMuFo4JmVF2yS1pGmOloXmgyF3GZrFw/fhZhv9/veFnE/wJcmqawSfNtaLNHGVDmhv72DfGNbZopC0KTxZALkA7Iw4+eJWfte7SvnCo1fL0XuP7fxMqVdP/RiTRWeP2TPsmkz+/ZUpFvrfgoUbFR+laBa9WBlFw1cgo5havqSZHWSlEH/rv5t2LxaSIxpazLalWXDiq1psJv2Ao0Ix2yCi69/VRWuuPwq2duVI4vm+blMb8+7Vde0Go1ablJ4VX7RNEtv6+1zvzKCvo9DgTQ6FhkpIMZYBwAbkS57laen/7n7k5yqgN1usDVtQX//PmIUvXBlMR1SywJ7vwCOOaNx3hrk0GhnSBMNxgFQzps++UnjJ5+M5fse0xoWEGVsus1eZQx4rN3N6hVT1hZms9gStIWUd/jPXyPNfH5DORKGkn02H6/lwKG8Wfdetek/vL+HcZDHKSxEFA9r4tq+7TUREnRYFvs9fHCwENep1iEQYNSL7dLC/7mtmrdarlRddumVQ93SiWAnb5cgs16J1vepxYIYPziV9g4t5arRk1l1Kk3cdVeU2KvuJIKEy+d+ZtuEzrBxRlsRaFh1bF7Renby4/KilZ1wvBOAOXvzBsZX7fGfnSLQsNKYG0i8O9zAnx3FALYvG11IuueuMIwbSFZq+Atf6OD2t2TKFphZZQQ5HziN6tCACe+cj9ZzzmQpN2q5+AZrGrsQoe7QghdCDFfCDG7XmUs2mybwJlUR9I94NCtEq5ahiozUEto5EOqL87W9bh5/6Z3WzMDmldwwx7f4dxx08kHXM4KKitosMRZpQ1Z9Yly/jC0au5lNw2Z+mSpgqj6UD2mDrJy6UxfMO37F3LbbgdHfqME7h56IGaMg/Ja4Md/GiWrfGcaSNp3OhZ3ENtflZeOH0+VZZQMk/MO/LFvfr9dkx8apI2w7U7fOXv7fRh1yo3MGbofy9b4W4TVgq9zxX8a8HY9C7B69+XMQ86gEHLDNSkcoWEINdfBUaulqQseZe7Vx3H77efy4tXH0ZLtxqjp/+SSfY+h6MN7uUxL6BRCXMpKwHvlJmhnIICGdt//YTuIqBWQAHq06zbLvtkF7sro7PEzsPRqfjfUhABuXAOZooMxU9rcfds5HL3gUSXrm++99WTVb3EFWlKkqY8P+9a0dgBR5XwTUNQMbE33vdnvVMgH9/H/CgFkKsaKAA5c/AoAJUcyoFdj6nx/LYJfCDEAOAS4rp7l7NS/J48MHcPUaRelTjvrlDATBHD3g4mkwbboUSrQaFtcOvtSAK4aNZXv+/De0UGkTTbick02wXAPUo/khRa4g4qqBwPJVl8tR7dLvuoKPzWKKlTTSlyzVMPH22gtW3UD6WuP7Zc26+OCIar/hPEW90yo7L8/ChIoKtKPU/+1oLKPfFNgOiU+6tUfbwQPd2cpPJNXDM6lZEDzCs45eIe6mHN+XSv+y4CzCTLvBoQQPxJCzBNCzFu5cmWiQvr0yHLpkcOQhkHOI7BqXTnlQkLV1boNNqXNkBWLAZg/YAh37HyQL82oy09+ArWmSUrTeX2zbRMPvJNfvl/pwpaF4L4dxwT6IvJC9ZvirB5Vt+lBPNRy1hJ1BiGBNs0gp5uUIs4Lku6mLEDqRofBQdQhqQrS3AlEIZX+TuezlaIQ3LHzQeSNDK1GloKm8+xWu3L7nb9Ctt8BadMM8rrBVSOnkPfEH47DS4MwbtceAAAgAElEQVRtscfSN1m2pq3GL/DHBhf8QohJwBdSytfC0kkpr5VSjpBSjujbN/mV5dHbbsKypn5oKfS6cgcoaAaN0l9b6h74rI+J6pdfZSBuvnYlR7zxBCe9fB9T3ox2oOYtx4u8MCh4tqNxhJurCrIYufydRINJApkIB1/lcjJIJr/7ItLxT1+PFV+SyaAeUOFD4FqEHTX1D/zmwJ+ERoyzUXf7UAkTyNqljtgJttDcySA2pdqxIcsM2qkJwJSSKW88iVEq0q1UIOPYjF/8Ko22RWP7BUoDmDb1j9yw53dqiiMscCNwzX7yDT5YsS4xnUCoXKBK8wEuApbhxu39HGgDbg3LU8sFrmffXSG3Ome2PGfc9G+sz++qSxsVl57SuvhUz4s69XCRG5beprb6qSyjVfhf7lP9TbWcIFfCSS8zqV4CTFJPQTzl0KQVEEPg6+hT9cxfS7k53ZDTJ8+U0yfPlG26KVv0bCJ+1pkNcvIPLpV3v/pJYvnHN+UCl5TyXCnlACnlIGAq8JSU8pj6leiuoR7bfhQ37D5ZeZUrfd7FMXGLA68FkCalslpCVjxFTfc1LZQJeA2ql6Tp4iJchRVdalQKiet4zhT+KWtR1fiVE6RbDzKdjML6s4Vwg8Ukh6NBPDXgoEsn0TkMCdLXgrDvjZIBtfBZNpK45OG/stm6VQghQBOuUYEiP2V0t/IM+XwxwwZuXANH/ujSdvwA/ZsamLzoGeZeczxTFz6OJXQsBK16hvJGzK/i/7zvD3h4u73r2llVJpcoVA5s07G5ddfxnHfQT8jpJq1GlvgefOLxkdbkF0eQpGHt4W7dnVRNKyWumk8CxYrpKezmZhzaSYRt3DL8yqyEt+7T5KeSVr2teAoIWvRsKucAUF0PWbvEr565gYZSke5WPrSfhVn8XfjUtWyr5VPgsDO+Vu+cUspngGfqWcbKj5bxl4cuI1Ohb5OAhQw9aHxs+725a9h4Dlg8j4aKK/tJOopfHgm8tvl27P5ZdZCFMBNKiei04q1MK4ATX3uQUrt/fqG7FkNx+KoHJK4ZZVTnr0WIlPMKz79R6eIgKo/A1Y2DW++1GhB427beSLrzqERYG6b5DbW0I8BrWwzBdCxGeMafxD3P8LP8CkNah/te6CUL5s+HceNSouiiy7tl/vjJF9nLc8gicANYeFHuTHfsfBAf9hnI5EXPoDlOp0GoonpRSWOhscvni5W2pPftOIb7dz6AjQqtXDr7UhpDfMe4wqe909rhnddun0Qqt31+/KQxQQiZXKyrlF/SdCUPoioTTJzvlbjumDN2KVVBncYOYUNMFl7Us8xKYe817Yyrytr707cC85Xkev9NSdSjcfN4Vb3e/M25Ik0xaKqgS6t6VrUUeOytz2Pnm/LmU0ydP4eL51xBRsa79B2kx4XOHSODgxlgGeTNN/G9uWycWwcStBRVE2WhH7V1j1o9q8BRWEEltRaSwLu9BybI7cIGcprB2kw3CpoR+3brB722iF1mmnpl1fL8VDdJVEj1Ujmp7KhUflPhzU9NWBa6De0uMGxEx12GOGrIMCRR7S7f6H8DscTCsjU5Puj/LSyhk4kQsrC+ETLS5vePXa3k/TEOkqoXGm2Lyx68pKbr535IMutLICcMMrKk5FfdHTQiVOVUTleL+mWnL5fEXmlZ7Xw5moGB5L4d9mVB/+2Y9ciVscof8uXHic5mKnkp6GYnlWKtUFU9hE3qcfOoQpI8CE8c4V8rbT1CHRyFNHjI6yab68nNQoPQpVf8A3o18kVjT86cdGYs3zdAh7uBesK7EwiDTu3ujdOYNmzAwIk1ILSIuoziq7KegtL6rd7C0roqMZevjFPCdGyOWvgof3zkbxHcpI9WzcTv5meYFUhe0ZV30ja3qjiKj7DdhJOi2wxveUHjJMlOJarfxtkx2bguHuLA1AQbD9kuVh4VdGnB36dHlt8cMoTZQ8Yw8YS/xar0tIS+SocoonWyAqlXOX7f5FXARB3W6QQHH/fjSUUgq1qIxJ34ijGFS+VkEDdfJeIKmO6OhXBkx0Us1zoomIuipiOkGp+1WD3VAonbty4dNY1bh02gUOFyXOD2Ib88tUAlyl4a31VJq/IJK0sC/x4yBqlV8yiBFrPBVw338EnnQA0XWIPQpQU/wNAtmsjqgg/7DOS3B/24rqukpMjgRKpComC1u/uNQ0XgdoA4K+lYh54B3StKhaC6iop8H+FEr547Ou+gD9Kvl9NmpN1pV+fXHzpWtI6tpLr05vMr24/HIhpS02NNaH6TuQ6c+sq9/HXfY/j1gT9R5jeMblAaCRCyIElD7RL0myqtSe+8yO8OOJli+25N4i5Ozh03naOm/ZFzx0/v9M4SGk8uL7CqpZCw1GB0aR0/QPeMTsF2u8+ifoNpNRroUUrfLrYWqK2eO6+JvYLlke32Ztz7L8Vu0DQHihdeFw2qk1JBN9BtO3IyDONPAA3tB+EbeiL3osxnic567SSTqmsYEL98VXWapek8uu1I9v/o9U4m0OU0JaGhBwhYvz4jbIcznr+FIxc+HpNrdZTLTDOWbxKU6zmnmTQ6VvXE6Tgsa+rHCVPOY6NCa0eFrct2Z1lTP5Y19eO3T1zb0b4Z6fCnhy5n8Qcn02fYtqny2uUF//Lm9b6slzX1w3CiD9D8BkqUrs9vdVQ284uTL5gniYXgnb6D2HnlR53eFYXGQe+/3CHowpCGQA+jbQkt8MKS+pkAFHXTN0h5LSs36fk7LpJYZFQi7kFhSehIKTvMc+t95gSu99JJ780NfJ8zs7yw5S4c9MHLnQRtEG9ZWeKYBY8oLTCCrGxUkUb9eMsME+R+aDOy3LLbBH786r+q3mVxuO7e31PUTRpKRTTH7jDYsITO5ftMw9LNTubaJV1nQPMKIF3B3+VVPZXNOHrJfDQfXXrQFjZIf+c91FnoE7mqoBv8dfRRtTBexUsGyc4rP6rqmBnp0BDDIVS9VsAC95ZqLZAAtk03n4hE5TKC8qVxRhCkFqk8NIziIQ1rGAk8tfVuIIQS3358+P0exUvlN/rR6VHMcbCi0K+kF/ddFN2o9g76fr/HFqKqf1WqW3QB9w0ZozRuNOlw167jA72mNtgWPYttZJxSp3gVGWnzsxfvIOOx7MpK+38Pd5Ogf5PrGrV3WzOz5lxBxhN4XZIsKlHlANlu1dKqhtYkfLzx5oH5JW5nyhkZcpqr1FDR0fupB+IKBj9e0poMkggpL8pBKbx0w2iUEBtMZx/ET06r9oBaRtz6FcC4D+eF3vUIO7soIShqOq16ptOkVSu8B5obYhdShgTWZRrJGxmu2msK1+8+GavdHXIU/Hgu6CZFzcDxOShfL5AdsrbFIe/9h6cHDQ+dVItCY+bE0/mwz0Du2Wn/2G2ecWwM2+o0Kd25y7i6hF78/0bVM6B5BZZm0Mj6lWS5sQ3UtvESN6qO5rmF6ggN6VXpSMkWa8PjCHzYZyAX7D+U3z3xdwpC8w3SEQd+21Q8v6mqYOo9oMsT7tyBQxm19K2a7igUhI6jG5gl/0MwFRVPGsLRkBIp0qu5KN18CUFJaKBpNNgeVYSmIaVEV3CFHVZ+EtRDnSiBC8eewKsDh7LT5x9w/uN/D9xdyop/w3Yb5XsTKj0vWyoydsnr4X2o3WKnd1szh7/9bNVYjKoTQbVA/t6Cx1j8wdLUdfxdfsVfru5lTf1C/WP7+c/3osR6z5mV6FYqVB1EFowMM5+/JbTjzZj7f/j9Y1eTdWwaYjryUhWTZXvsAmorow0Fgavz3nfpm52sWZIIDE1KGgOEPrhtW66HetaBLm2EY4euCuMiSM1YQnDr8IkYAt+LX4Zjk21XAdajX4VBRUVT+YS9z2kmed3g3HHTuX7PI1jW1I/zH/t7h0mxVy2lopKLw29lmii1VNYucfGcyxn58cKqPWvSidAWol3Hny66vODfqX9PdAGruzUxc8Jp5I2Mb+eOsgiQwBPbjvR959cpTLtI0RNf1luuLZMFgw7K4Xjelc8F3H9r073H4SdoQHsRdobSYmQ72bUHlZvXTdA0351OGTqdY5r6IUxvr+rhVINOAqmS97TUcW/3HsjNwyZw4nd/w9ELHu1wyxxHNealmyRfGI0oeNvd22crHwOHicdfyZ27TQBgyIrF631RhdBMk9840B2Hv86+JHQhEgcNTul/dfxJ0KdHltMPdCtu9pAxXLf75Ko0Ktt8C42xi18NLaty1WE4Dlmr8wFl1U5BwcLID0HCUvP53S+PxL1BWBQ6azPdKAqdUruuNE1dfyVUdlRlWELj/AN+xIxDz+b8/U+qErzuStDgbyOncNnoaVgR4Qf96irq5qsECprOhWNPwDGqI5fF0SungTKtHVcv5YcL5nDTvRckjvAUteqOQ6PMm2pZ3kt1UXVkODb9K1SmG+VbE3Crzm9cWl6YTomGhCpbv8XSglPPrssFri6v4wfo3d21jO3d1syJ8x6IbBS/DmLiEBCzowOV6Q0kUtpYaBgEmzfWS9CGoaDpnDTlPBb1G8yA5hUsa+oHwImv3M+PX76nJv8kfpC4q24HlDyBmtJh1qPr/eU4QsOSkM80YDo21+1+KAA/efV+X4+cKrwv3HQbdl35EdkAL60vbDmU0w47l9XdmhjQvILjXn+og3aa5zBxEbazUaUfdeajQi9K2HthI7h89FRmzL0rosRwjPh0UaJ8afXntCYPP1QvlgTTtaE83FJIPeB6lxf8q1oK/ObfbwLtB7weO1kJWEKQqTiwLRGvYsIGh4ET6d87je12HGQdm+U9+7K6WxOru7kOX3u3NXPCvH8n6hBRaoby3yoXbPy+XZcOGpAt5rhn54M4cd6/qw8zY0AAe3xeHQehEnsse5sTXrmft/oNZup/H1MSkPftOIYjPId6fmkOe+f5wEtQtSLqcH9DoIiosswykJz6n7tBiNirnVazgV0+e4+WTCPHzH9YWZ2Txnf7qaHSoKlCxzJMBq79gmVrcv8r+OPiP4u/xGlvvWVN/TA9K7yCpnPxt3/IWS/ciqUZmLaFIR1ExGUoCVisD7wRBAEd4erKf3vfFzSDO3Ydx1ubDmbvpQv5zqJnO72vB46f9wA3jjiUHsUcy5r6MW3+nEDvkN7O/3XAPaeAaW88nhoPNsLXA6vAvfU7/eV7QvNLoFXPYmDzuwN+zJ27TWBJ7/6c+eIdvuqlu4cewNmHnMEXPXrzk1fv96X3ddVvUPlJ+PEzx3WJiSrtfNQ3WwjuvPOXFPQMDVYOs4a4DklhIZC6TjbkMmYZlerUWtuyoVTkKz3LgF6NNVKqRpcX/B+vagPcFe20+XMwKoQwgBQ6Zz9/C7O+fSyvDhxKz3wL/3PvhZg+1v3eLmeibhEQhoxT4nsLn+Qo+2Flml6+vCu9KAuEYxbM4ZgFc8gbGQSg2cF6Sb9DuFp5+CbAe7Du5TmK/4KR4bxxp/DM4D0A2OWz9/igz5aBbsDv2HU8Fz18GVPfeCKVuomj0vEiLSEfB2ELi6CyTSTCLikJXQnkdBPdcRDSwc+MI+4ipmwg4TXXDjJWKMMvOE9UuX71sMW6cJPwpOjygn+rPt2YNv9hfvv438n6bK8b22+I/vrpGzh33HRe2XIoDT63Rus5KATQ3Q62QY8awFb71rqSnkqZAI2lYgcdFYTxI3EvoTkINOn4emDcEAhTnYUNfBW9dWUazXH4qFd/Ri+Zz6w5VwDuKi0Id9x+LlmfwD5l2os22ZIhX34SyVNak2wa/dqvXsIWEEksgOJAR0Y6sAvrvyrv/ej4fVerkaWHgnWPpPr8qwxd0+qi6kFK+Y1/dt99d5kU6y67UjogpcKT03T5mwN/LNs0Uyn91/E4nueeHcfIP+97jPI3htFNI12ZrxcGDpU53ZBtmtHxW9xv9P5WQMjnthyqxMNF3/6BPGv8DNmqd27LVs2QBaHFrp8Cmrxyrykyp5uyqOmd2sCKQS+Md5XvsmLW5Ybqk6rv0uTdOxai0haELnOaHkij3LZFIWLz2erpF6WY3+qXtoAm9zztdvnlunxi+QfM85OpXducc+VKepz9c+VVQ9ax+cUzN9GQ0MwSCLU5Twui4jni7WfRY4RjDONDhUcV9YIARi19E0fo6MClo6dx1cgpSvRlBY3y3xL3dq4G7PPJmwpUoG/rGi584pqq+MSagKIR178lIASf9uyHdByMCtv58hnOhoCgdg+UUW3gfS/p7MsmLF/5CYPqwWya4wXcXbEQ7hjv9LvQsDTDPdexLUzHRgqdEpryN0G7e3Nn/V0Rr2WcxDWhLrtniYIEhBD8oeeK9Ff7QNVM8E18Eq/4X3lFysbGmlYT9ciX9oqtiJAFz0omqFw7pPw4q7M49eOAfGTwnkrp01ghBu4YhCbPGTddthmZ2DRb9KzM6/47wSS7mQ2xao9bRnlH4f09hyYLIrp/hZWp+s3eVXwaO9mg/uCAzLXvSGvpg+5uNDxdi56VZ42fEbsP2Q2NUn7xRTL5J6Xk/8sV/6BBSDt5cPIk+k8Vky/varYS3r9VYCAxFVb9guhLXqqIY9omgHGLX0nlIFwVvvp7TaO5oQfXjziUvB4v9I0uS1U3sYPKiuLLr+78+kEUf5Lwi3Gqq+u80CnoBhfudwLS5zJcFofXttje14mg9+9SAD8FzQisPy/P3sePZz/4/R5Gx7XeiufSwg9li7MwdLMLXPj41TghvpxKfg7nMiYsWVIbgz7o2oK/b1+Wn/WbRMI1iWVNXKFdRGClEHvUTzWSFtKmF4euSrlx+DPtElc/8Gemv3QPunQoCDXFiQRuGHFolSlwmkhqYRM2kUdB4qrhjjz2YvY+9WYa7ELghbi9li3yfeeFEcBP1ilx75Cx5HSTNt3fbUoYn2kjjmVN1CStWl7WsQMP/0vg6yVUsywYNChBieHo2oIf6D5urOvPJQISNbfIYYh7JSeDrPIwGIdCkNBLc9VcrxU4wCPbjqxZJ+yd8PwGbmXa8mM6NtkY4QtPeeV+dEfG0vumgbBykrZPmf/7dhzDs4P36Li5/dO5dynVnWR9jFjVCUsA33vrKdfJoYCSpgfuDvzyqpYB9WmfSnq17pgtTSfX7jOs/BhUh9uUQOtZ9XHZ0OUF/+pN+oOiq9zntty16jfVge6qUWTsTqey/fXjqYTAz7K53rbYlTzUks4GLt7vuFR3PCoqlEqU0CggKBB+eNkxWUg7sBxVxJk40px4/ege8faz3H7HL3nx6uM44/lbkI7aIXWbkeXWYQfTpsU7JDccm6xTorFUxHTs1F2DlGGhsbxbrw2y0y/TVR0PArAdh2tHHN7xd1if+mS/CQk4ikaXF/yv5wxmTjiNnJEJXKWUsb+Pv+04Az2JUDAiLqb4daqibrghDkm2yglKW4uuOAhheb9q7MmZk35OTjfJt++2VARjWqs5A/eSTzzFQ23wrpyDkJSjOPkE0MPK02hbHLPgERp9dkB+9DJ2kRNfe5BuTvCdhaDyVBBXDeRNb+LQv21NVVoHNzJeQQ930BdF34u4C44sklNfvjeybAuNzfVkjvii0OUF/7CBGzN7yBgu2P+kUB1tvVZXUYhqgLxu8u8d9qWgG6zLNJIzMjw6eCQNnkEaxb8E2vQMBd3g6UG71ch17XWl47rYnT1kDKNPvYkjj7mYI46eRUHTq4SjdxDFEZp+A1JWvFMVxGFImi/p6rIeCOo/ttAoekwbDSk73EHXixdVlBBctdeUTof13jZt0zPkdJN/Dj/ENZWhPjsvFXVQ+Z3pc5HPCxOHjTftrVByfHR5wb9tv434yZAenPfUdTVHuKoFfp1NpUM12BYT351LUTMxbYtL9jmGcR++Gno45VfufUPGMPXoPzHh+CvZe+mbdamHuEJpo3wru3zmOktbuPl22LpB0ehss5xEKPsJ+/L5TZhW3yK+UKgsI84uKopWUPsUIs6h0mxXCVy033GccOTvOHv8jE6Tsl/asL/9YCPIG5mqiSUODCT3DT2AB3bYNzDNP0YewSHHX8Gx8+fQ4JRosEsdh+KqZ0gqiFp4FYjZpzQdWlpi5FBHl3fZAPCLHbuRM00IuU5fbwhc4XP1Ht/lx/PuD/Rb4pcvI20ylhtCcubz/8TSTRpilAtwxKJneXXAUBZtNrjKQ2lasKm+uBKcVvDXh/9KUTPIOCWuG3Eob/UbHLgrq1XdVL50FdbhvX6SikKg4x91zY83v8Plyv9Ln9/iQpOS23YZx9ELoz2GViJsMvGmy6Gha4KHtt+Hmc/f6nqXdUqRbskrESUEAV4YNIwzJ5/FtPlzmDH3TqQQVeNChe8j3nySQ995PjDdc1sPZ9x7/0n9op1qnZbxyPb7cOi7LyinF47NV1+1sHFszhRQ6+WquA8wEHgaWAS8BZwWlacWlw1SSnn3w6/JnG7UdBEkjWdtplEefeQFMpfgAlH5adPN0EsgTw4aHngZJKcbcuxJ14ReYEpyyary96hLXwWhybymy2LA1XlLaLIgtNQvOAVd4Ck/BSGq6mWd2VAXFwNRdRlVh616JtHlLL/f/FwLONDhakOlDgtCk89uuUtgGWH95QffPa+q3p32ui9ourQQkd9mIeRas9G33BLIgmb4XtSKU3cFRE0XywpC7Vu85S75y1U1yT6+QRe4SsDPpZRDgL2A6UKIIfUqbFVLgV++8DkX7/uDDWqCV34qkS25q5mL9z0mMS8NtsWsbx9Lzsh0mIQV28sqCo1RS98IVGdYmkmPYq7jsDsvqtUG3hVqHD6jts6ubtNBd+wqm/ByXkM6OJrOVSOnRKo1guDHtx9fjhDcO2Q/jj3yAg4+8eqqNBkFN9VxUKm2WpdppBAQOSyqDrvZqsEgw2kKXB1+q496rdHnYpMfjUe2HYmj6YwMCJASxed19/8BVz6tR0HPYNoWRWH4elD1QkfSYHV2hlauax3X+63fRS3V/i1wTS3zQue64ZP5xx7fodXsvOcOoiNxo8X9c/jERGFW++TWxs6jgg0u+KWUn0kpX2///zrgbWCLepW3bE2OQ99+lrOe+2c1L1Q3mGrThKULsu4xnRLX3H8RM5+9mVpOGzZfu5LRp9zI94/6E89uOaxDRZGVDg22FejPxXQsWjKNzB4yhkOOu5w7dx0XWk5eM8gl8WsTgrK6JezrNdvmwz4DmDrtIkoRYRWDhHyl+iWID11KvrvoGbb/YglHvPEEesnqoFfUDGZvv49vXlX78yC0mI389sAfc/qkswLTbKhFii00utcQH/bgD16mwbbIBrgghuC2ELj3Kbxqx6xdJOvY9LALSpOP294OVszwoXEMOgTQKG1Oev1Bjlkwh+5WPjKPBO4euj+jp9/MJ736K5bUucwef7gAVtbBNbPfNmBDPcAg4BOgp8+7HwHzgHlbbrll4q3Oqg+XBqo2XF8d1SqHJNvBJNtHv99UVS0XjD1BPj9wZyVVQZlum5GRbUZG3jD8ENlmZGSLkQ0tL46PlrTrZp3Z0M7rJGm1q36CVCXFGPWsWu853ZRtASq1XLsXxqTf12Zk5G4zbpNjT7pG2W9S0jQbgkbaPNTLR5YqjSQeUC2qPXo6IK8YOSWyrcOeYo+NXJ9jCcE3SNUDgBCiB3AvcLqUsmo/I6W8Vko5Qko5om8NN9e0Tz6mFHo1X4T85eEpMRdRpcYrR+DGDxi99I3AFZDfb42lIo2lIse9/hCNpSLdS9UrqkoedCQ+3kPqjg778lKRqf99lJKPiWdlWqkZ5CKCp0eV560HW2i+uzIJXPztH1a9k0Chvbb87MRl+1MUOjMnnMbqbk30KObI+fTNoO8o08gZmcgA82GQFY9K2lYjG3uXo0ofqus+yV44bPUeh4+yeqjsA0k1r46k6FGdCuCnL9/D3GuOZ8iKD7lp+KT41mK5HKs32TxGDjV8LYJfCGHiCv3bpJT31bOst7O9MO1qax6JG3Bbc5zAgSCJ3wGCoJI/SEUUllalTJV0Xv7Kes0N5XI4CLbQsCOEnCMlrwwY6vvO++0lol0MAxiy+hyijOVNfXE8N45tBFITtOpZcBwsj0Avt5ej68wdNAxwQ4F6rYzC2qxMwywVsZ3aeqRKXyuX8PxWu1Jqj9QWBYnrQtsWWs2mpWHjMs7XJzkPKasj46iC/ILsCNzgPJfMuZwnt92T4757Xuht9arvciTLm6PVSnGxwQW/EEIA1wNvSykvrXd5Vu++/M8e3/HtKDt98SGZdg87QY1cae+bpCOXO6kdU//opZEkXZxOW86/4df34dDtEmbE7eastPn2JwuUJribh0+qEsp+6f428nv8z57+/WbsB69SNDr7f9KRNLTrpbPSwZS2r7fFomYwoHlFx9+aj6+mqO/QgYbYnqE6l6GaTgDjF7/iKh4qEL5riF4w1LIQKqJTVGjDoEkjLcTZOWVLRf5+3x/4+78u4u2+g5TzFjLZTv0lLXwdK/7RwLHA/kKIBe3PxHoVtlP/nty053eqZtnyli4MSYW9Hw1DugO1pOnK1iqqK/Zy2iiarqpBY+Gm2wQeiKax7U4LEnflpTmdBb8fj34WG354Zts9OG/cTyLr6qcv38PWq5dX/S6A7y56JvJwzz08rhZ+plPqcIo28uOFiSw9vg5kfSycgsaHX4hTv7xlxK2BDLZSeMWg5UKaC5ygvuc3vnqUCjTaFrt88aHyYrIRh42HbJcOsxX4Oqx6XpBSCinlLlLKYe3Pw/Uss7lHE9eOPCLxijhNGI7NLbtNDDTliwuJq764dNQ0Slr4fTyBa045pL3jfdPh8mvXHHWqDEszWNRvMIv6DabFbAwtN2tbHPL+XP9doEewBQkR78CWwPkHnMzqbk1MXvQMf519SSi/G3pKCPuObxJU+fEbDX4C129XYCH8/eMr8hfGo8o7iWsGmrvm2v/1zpkEy9bkOPzt5zjp1X/XRGe9ysb1od9iZBN3iqMXPMJ9O42N3O6pqC6eHjScPWfcyrPb7kHeY3oZtPJI87p2WLbKYZkAACAASURBVB0kOQxMazXmeOiVhMaZh5wBwPYrPqqy+/ZDXJWICgY0r2Cfj17n4jlX0FBnFyIqO8BvCoqaTi5mcJwoxFV1VuLfQ/bjir2/n9oZXxLoAla31snbgJ+pzzftqeXm7qoPl8qC59Zu1M3IsL/bdFMefvQs+dC39ooVSq7yt7WZxtCbxHHMvsqmgbvNuK3qRnBSU8G4pnG1mCSmkcePRtETAL1NN+Xv9zu+w+wuqflfrWaDcW7epmGemNZ3ed+pmjyqpLEQMo8mc6LaTDaorQpClzkjI3O6GYvverRzrf0h7H29Qi92eV89vd9fhIw4HKxEpI5bwp13/JJMjJBt3nSZUjHygDEOBjSvYKs1y9Fsu2NlYgsNW0r89iVx9K+qkAnz1atsR9PAXq9j1x2HXz1zo5IVSxQftXynwL15q5o2CFH1Xd3q0WUF5ZG4O6bKoEFxem8YrxL3YNxA+jJQzldEIHWdom5iODYzJ5zG3EHDOOD9l5j1yJWBZUeZeAa99/s9Kk9ShNF1DANtyZLU1T1dXvB/urqN+HfmghujwbFCO4VfHm863bHJhPiJjNPpGkpFdLvErDlXdDrw0qWTmm48CmWh4eUx6MC13jrjEoKMZ7JXcYNbEoJ7dtqfqW8+mRovUUKvoJsdB6dx6qWk6R2O59JCKC0hkDK4betSZgXymUZOPfwXrG3oQUumkR7FHBvn1tJ/rf+t1qi+1opOg5AYMcyV40wSSWl6+4uoU+jFLi/439B7Kgl+vwHq1yh5PUOjz4pNtVN4dewqK2XJes+XfjjsradxZGcqaQoEPx6rOqgCHQtRFV4uDg9B5QhPGoQGMcIqlmFIqRRXVgVlXhzCV8cP7LAvdww7mP0/eIWTX70fKQUN0n83WZ4ohJTMGvMDzn7uliprm3K7JF3p+kHgulYI2xHUa4dUhmlbLOo3mNFL5jNrzhVIKWm0LQoBYVWjyuyOTRE9tShgquNYpayC0LF0A106LPrdX9i9Doe7Vbqfb+JTi47/kzlP+3r/i9K1+f19yy7jYnvYS0M/GKXDrEUf6YC0RfV187g8qjw53ayrC4g0vkHVg6Zq3aucpxSELguaEekNtLK8gtDlnMF7Kp/XpOUWoV5tF1VuESFv2XmczHnObpLwrFJntdZNUk+eDsi80GSL2SBzRkYu+svfE8s+KaXkm+ayYUPhzYxaBBuVmfi//bdP3e7aEjr5Sk+bQj0EIURfLqukYwmNkk+TC6nmb75MLy4kbiyCsyaeznkHhdvQq3xzEKK+QbU+45QXx6IHn/IF7fEWnBI9rHzkN4iKPPsvfkWJxxKQ0/xXxnF4LSX0luql6aUR1dcF7i3yo994jOwGuEkepaqVrHfUV+ZZeN4LkvUPcO9BdLfyNJSK7PDrM+ripK1LC/5VLQXOePrTDjfEazPdKOgmuYDtYSW8jWVpRqAJoATywqgS3EFY33kEZ046k1HtnjYPOOkaphx7MQecdA1nj5/Rce09atseVk5JCK4aOYVjj7yAI4/6Ewbxb4p600eV6YfbdhnH3EHD6JP7KrKMsimmCuJMFKrfmvSAOe2JJQrewDFBMIBGp7NKKG49SOBPY34YqFqpTBdWB35tYCG4a6exSvzEnWRnf2vv1E0xzx97IrZmBPITt439JpBOtJYsiUkxGl1ax79sTQ4pJbOHjGFRv20YtvxdPurVn/9z+7nKNCRucPOfTzyDzdatCkwnpENOz2DYRWxE6M6g3Lg6kk979mVA8wqWNfVjdbemjjT9166s2U+OAEwpOfG1B7hhz+8woHkFOd2kWx2ib0XhhwvmcNSCOQjCfbgIXJ2436qw/L7ytygdehjKK7MopDkxfl1IIowq87SYjTRYxQ7fVkE0k5RjIjnyracjeVDhs/L3x7cZwUEfvkpeNzsieyVZgXux2/J3yDj+loJJaIdaVOVyiB49ElANR5cW/N0zOoWSZPKiZ5g15wpKmk7WKlad5Ed1sFt2ncCifttw8ZzLA9NlcaDj0Fd9jXHX7b8gb2QxpGuiNnvIGAD2+0BtG68Cq90/zLKmfrFXTGkIsjKNDCBRm8yiDtrL8G5Z4/CcppD2O3xXFVxBiLLu2pDfadgWP/vPnZ0sx8ImAFWkcSDqV4c2cNCH8xKrKSXBoUQnfPCy8jeXd/YG4erUIBSETv6L1Wy8Y4LMIejSqp7Wok2/fDOz5lxBY6nIRsUcGWnj9c8YtQI98fUHefiGn9aFR1M6bGTlaCwVufShS9lt2SK2WbWUY+c/7GtJk0QH3mjl6ZlvYePcWu7e+YDEW980tsxpmx/GnciS1qEfrUo9b9Jdh1cXXPkE8Rm2X0vj28p8VdIzHbvKRFZFZeqt8zj1r7ob83uigv1EwQG+d/SsqkBAJaFHqrvKkMCFY09g/ElXU4zIUy2V2mloGq84PdWYjoEuLfgH9Gpk4LqVVY0Xt0MIoMGxafAEa1dZ8QRu4XzymY7NfbedzWPXnUomQM2TZKtqSIcb7zmfJ687hSMXPqGcT3XVXUbautQ4KAsTRwRzmfSwzVuGBCw0nAo9rx9dv4k7DYSJkCSqFgf38N0vjsB6QVq9YvVGzvLS9R5ygruCfWAH/8hmUbRqrT+V/lv5bP/lx5xxyJnkNaPj+c1BP8FUNPltNRpYvlFfhi1/l1nfPpa8kaHQfgZoeYw4dE/9ls8LZ048neWZ7nE/NRJdWtXTp0eWk44eS+Zm9Zu7UB/drwrK9PxCiqh0+iC+BXR01oZ23WTQpFVZTpIJMowPVSTJX06vSYnVbn1SOUDTaqv1aisH6cQ7g0n6TV5Ymk7G6axySfp95XzCsVENsqmy4Anqhw3SZtz7L2MJrdPiJvRw05MmbfiVbQAXPXoVcwcOJVuhzz/4vRc73ZAPq4OGUp6rH/hzx98OFStt6WChkQlQfRaEzknf+y0vbD2c4/o3+aapBV16xQ9w8P67MO+M80JX3n7wizgUV62QFzo3DZ9Ezsj4RueJgzQOpVTppV1WXN1qLbQFrgdUu90Nd62WFxJqiu6livLKuyj0UAeAlZN45W9hf6vwHifwiF+6uG1cqcqyhMaFY0/gqr2m1KWeVSytgvrIqKVvdtq1jFmyoF1dvP4b2oR/xDW9Ip/w/G0iMUPOu7LSZnlP9+KWaaR/B7/LC36AtyZN49xx05XDx1U2UhlR+fze37PzATy57Z4cctzl/GLCDHJG1jdfGvpmW2iKx6bRdKLSRKFy5V/5twriTjpBA7bR9netoYpyuxTQaWj3yxRX9VWmo/p7m9mArWncuPtkCiFB7osh7oKT9KU0JnnvmUAYGm0LUzqdhOfr/XfgpS13SYGTYN4qkZb6CMDA8VUx1lKveSNDj2IOgauyThtdXvCvaikw69F3aM02Yukmxfa4qHnNIB9y4FLrtlwARy98jFvuOo9Hr59Or9zaQPPMOHbrQbClTKUTR5mQpq3e2pBQGeyV7wXQ4Bt5V72soHeO0DrCQJZXwD2sPI22xcmvPcDvDjiZnJFhXaaxik7Q+U8UP/VMX0alMPfSKwYcYJrS4b7bzmbmMzdF0laFX/172zYNda7AbQ+tffxJXFWcqgvHoH5iODbLmvqha/UZcV1e8L+1fC09W9ote2yLbHvIuqxT4tFvjQwV/lFQ3UKa0uHXT9/A44P38FVNlDdyFtCiZWIPOgFkkDU5Zas8iFNJm+aKKU5ZtZYZ9Y21HgCrlCUoO9GTvmmypSK925oZfcqNHD31D1y43wnkK27NhvFnaUaVgE3yLfUQN2bA95Z/2yVmgKCwfhBUR0WfSHyV9ML6WVjf69RvHCfwnM7bLkF8Go7NuHfnkjE0lq3JBZSaHF1e8INkQPMKSj7Brw9954WOix3hFNR+C4MAxr3/Uuh7RzeYOemMxCqbDbkaL4eydERw8HK/lV8cVE5GAnfQ+gXxrvfkE4Q01S2VEMBpL9zBAe+/xLh35zLz+VvQQoKrVwqqp7YezqXfPpacj6poQ0zUZXgFaq3nAl5I4L4hY9yb+Ib6YilstyR8Hr/3UUhqs+8t68LHrmKjtWvqourp0lY9AN1MnWVN/XzdLVTqJaNWgUEHvXEsKkqajmkHd7yiZrJJ6xrl2bhW65ko2lA9gMu/lQeQLuMNXwkUNJ2sgkth7/swE9cNjfJ5SJBb31p5yki7w898nHoav/gVNzh6DH6izmLq2c+SoiXTyM0jDuOqUVP5wbwH+eGC6OitQX05Dvwmr1rPpIKgA7/677306TEtRi41dPkV/5JVbQDIEPtuFQTljtWQFdYmfshaBb7s0Ss2T0lWTpWrREtovNNny06/vdNnIFb7IWKlHrrWgWKh8Y89vpuA42Tll3kvIbhj54OU6ypqdVxW1/jlU/0tzBo8SuUUxF9lPtUVvor6S6XsNKCqzjMcm50+X8wj1//UV+irtF+tk1lcTUDcehTAwc/9iw/f+igua5Ho8oJ/6eo2BjSvCLWS8Fok1KNjS8DRdJyQQ1gNybt9B2FFBE33QrUDB+l9Demww6pPOgmNHVYtJVNheZEGJHDe+FNYtVH15OYK5/Trv8y/ret8vtEmygfpJU2n4Ln4pzLQ49SXV/0YByrlpNl2SelK4I6dD1KyqFM5Xyn3jdc3347znry2k2ll+X0lrTj9KI3zJC//YXS8fPulKWomH817KyYX0ejSgn9VS4Frnl3MsqZ+mD7BU4IGreoBp0qa1vYA0gLXjC1MG6kjOeKNJ6HCQkC1zCh+JOBVdkXpM1URZ3D0aWvmhUHDfN/VM2JYxi5xxtw7QgNvlNVQlqZzxd5HcsXoaR2mk0EDNwwSKAZM4gLXVruro1XPMm/gTsw47BwKPuLGEiLQ2idMWI5a+mYKmnR/2qpyIM3JN6+bLO/Wq+qbdRy2HrGTAoV46NKCf9maHKausbpbE7cOmxiop/dDmHCwEVWrwfK7SlhC59qR36XVbFDm+eR5/6q6IKLCr+VTvhfVtwjUYSM6bsR6oTr8BDDjP3e59HysK1QGXdLdgOoAzDg2hmNz5ot3cPoLt3eyN/dbQVohEwPAZaO+T9i98a/rYDoKcVRiYehuF/j9o1dz+b9nuY4MKyBwo56ZPiJcpb1MHw+ZX9d5T1Q9BI2dMhpsi/5ta6p2C//acQxNWyUJHhuOLi34B/RqxG4/fHxtwJBUaArcVQo+agBvAHVb12nNNNLdyivRLqHhRDRJ0E7g1S13juS71kHxs0lncdsu42raeRQ1g32WLKjynxQFiWsHftVeUygoqkjiTBKVZxjlp6zqCoMZoQ6b+cJtgRYUURNcEqioKGqdbNafmwS/z7fflxFAd7tIJmB9Xku/DHJs5qVfC5Lm9/IV5t8n6PxMAN9Z9DTvvbk4IRchiAp7+E14agm9+O/5y+Rph82UbUJPLXTcWrNRztr3GNlmZOTaTKPM66Y8Z/x0ecPwQzqFW7tnyBilUINuKD1NXjD2hEQ8tumm/MF3z0scMk4lNJ8DskXLyDYjI+/ZcUziurRA5hTbwgHZYmZlXjfklXtNkbvNuE1udc5s+UafrUJD3VlCiwx7lzQU4dcVflC1/FpDeao8BaHLy0YeKfNadd92QwemN9bS4L/e6YNoFBBybaZRuQ/6PWvNRvna3Y8mln8EhF7s8uac5upVzHrg0k5+xMuQJJvRDWlzx7AJ3DFsQoefe4ALHv9HJ3pHLHpWiZ4AEBrv9h2kHChF4vqQEUIwc+Lp7PfRax2/l+Eg0BQ0oTMm/ZyRy97iewufwNYNsla+St8ugO5OERyY9O4LtJkNyjuZSui49RcFCRTQyJQsrhj1fa7c52h6tzUz4e3n2WnVx77fdPOwiTw9eAQ33XtB5Eo66Uru61AlAJErW9XvksDTg4YzeukbZDxuLcLGgwRazQbXNbOUnPT6A50cxZWhenYRVJb3d7/vijr8jUqTBqJkh8C9sPbapluz17JFgWmiYEibQbuno62oRJdW9axqKXD3dQ9ihnTEApAT0bFEJdCiuZdFZk44jY1za9lv8avodokBzSsY+fHCKp1jnM5nOiX2+mQhWgyzUwEIIejXsorjXn+oSlWhISkQ7NeljJHLFvHKwKGuyasjQw9AwdWDd4sQ+hJ/YRU1aIusH1QNOJjS4cwX7+DW23/Bi9ccz2UPXhyY/58jJtMntzaUL+/NTT+eo+orKF9cxMkTdfYR9M6vjNFLF/K3vaZQ0NXWfRJ4ZPBIbh02AQf3jkH3iPjAcVGuQ6/QD2svP/hZ1bzcf4dYbeS1EApLE0Vnr2WLYhtKlITGOrORnJHh3ENOR/btG4OCakE+24Bv2pNU1bPgkzXyxKMuDN1qBakD/FQBed3wVekUNV0WNCOVLXSB9WoKC6FEs6DVtrXOa7rM+WzbN/RTQJPX7PEd2WJkI9vDrx1vGD5J/n6/cHWZSlvHeRyQRTQ5Z/CeMq/pssXIdrRbUFlORVunxUcQb3MG7+mrWnNA3rbzQR39NkotpqoSVOHJ77f7d9hXrs001kzLL01eN+Qlo6dJC5E6/bT7lgPygv1OkJN/cKncbcZtcsffzJELPlmTSP5JKSX/P6p6BvRqZGHfbbA0wzdGZthM7AjhDkvWryKydonfPXEtWc8WuXxwIyP4kdARj9evbFN2vs1aEhqOwHdL3YnXqIIjEEVfBZLaDsJsIUBoHL3gEbqV/IPaB6Fc7nGvz6bkceng1jlYRqYqkE75fUHoNCQ0rXS39A7jF79Cq9mA6ZRwEJ18tVQ2z3u9tuCyfY/h5a12Ydy7c/ndk/+gJAwMx+KZrXdn/49e77RzrKVOr999Mmu69eSAgJu8R7z1FI6PC4yg71T5XeIaOXj7cvldEJ0lTZuRtarbpzKvX3kl/MMjVvKXtUv87D93hcbBrqX/JqVTyY2o+Pes5/7J6Ok3s7pbE6bt1MVlg5CyRqmxATBixAg5b968RHn/+ti7bHbWDKa+8XjsBvFLn9MzNNjF0G11lDrjoW+NYtziVzCdcB9+rkDU0COsS1Q7W9LOrfJdKmUEvXOjWWmdAl5406rw7senBM4feyK/fPoG3zsU+fZ7FkHRpJLUWVQd5HWT+4eM4YhFzyClGxwn337B8Pz9T6a5sQebtKxh4/w6znzxjkT1XkJg6zoZO7iPtegZhKZFntXEqfuCblYtjKLo+AlAFT5U84WVrZImzviy6ewHp4RACtHh2kOVzoPbj+Znh5+LAOb9+kD69EhmjC2EeE1KOaK6EJ9tQL0f4GDgXeAD4BdR6Wux6nlj/vsyr5upbfHKqpigfCrbyaittV/6Ash1RkOgCiqNrWpYnrDtftLHAZnTTTlv8+2qaBY8ljlPD9pNFtpVFnH4cNvEX0VRVg9NnzxTthkZ2eqjrkuzDlXSFDRDtummbM50k21GRr7UfwclS54nB+0mc7rRqd6ieGjTTZnzGRvl/Dkj842ydopTRtp9XzV/kDVcEpXP2JOukVudM1s+++6KxPKPb4qqRwihA1cBBwHLgFeFEA9IKf2PvmtE0+fLsDSdrM+KTpZ5ikEvLGqOJXS3ySIQd+UhAKkb/OSIXzKgeQUXPP6PTuoAbxDsOLRVeUvbSkICV+01hScH78l9t51dRd/RdI6Y+ie2XrOcBf2356vGnsy95nhEab1KprzNj+Lbr5OX009d+BhXjp7G6FNuZEDzCvZY+iZnP/dPbKGjSxsch2xFm3rrL6gP5TSDkm7Sw8oF2q/7wXRKZFi/A9lz+TuRKhNLMzhr8llstfpT37r0QuLeUJ458XS2/OozZj5/a6c8bWYDZ004jcZSgfMf/wcbWbW5BU4yzmqFSFCuJXR35+nZsVQe9sYdOx3WcAkwbPm7fNhnYESpyfB1WPXsCXwgpfxQSlkE7gQOq1dhzZsNCLw8YeO6sZW+b/3hZz0hgbymc/k+0yjEuKVbmT8KRc1kbUMPHttuFLYQVRY8ZTqVjxXTIqKSnzh1koROq9nAY9uNYus1y33z3b3zAcwfMIT7dj6QD/sMZEDziir3B22Zbpw9fgZFBZ86QbA0gwHNK1jdrYmFm2/H9Xsewd6n3szUoy5iwvFX4ngsXyw0CrrhWl3oJrcMO7hDTVNG3shw8pTzOHraHzh3/HTyMdwGq6IgdFrNBvJGhjMPOYPV3Zp869IPJU1n4gl/Y/aQMdw5bEKVHytNOry81S48M3gPZdNbr38pCRS7dSenm1w1ckrNcS+ihK4f2vQshYhyJdBiZMkZGc6cdCajTr2Ji/c9hoJPuM0kKr9a3i/ovz2GBjv17xmz5Gh8HYJ/C2Bpxd/L2n/rBCHEj4QQ84QQ81auXJm4sM23Hcjvx/3Yt5LzmUZemnxM9Y1bBHkjQ4vZoCaUdZOJJ/yNO4ZNUBooXpRwb7Suy7gmXH5eJA3pRuQZ0LwCK6AztxpZzh4/g59PPJ0jjp6FrVcPRhUHcHnd5NHBe0aaN0a9u3745MCBp0uHZU39WNB/e9/3N47ovBZY1tSvylzWdEo8ud1e/PyQM8kZGdZmupE3MrEmPNMpddzDKKM8CXzYZyAzJ5zWQTtnZDhz8s/Z+9SbOXraHxh96k1ctu+xvnQX9RvMws23485hExh1yo1cvO8x5HSTdWZ1RC0J5HSznffoW8mWZjDhxL8xbdofGXXKjcweMgYgsC69ZV02elr7StL91rM83zhzwmms7tbE6m5Nnb6/KDRKHo+tlqZz7vjpnHXIGZ1onDfhp6y4dzb7/+yfXLLfccyceDo5I9PhO79NNykKvYpeXs9gaTpFobkTm25y67AJgW5PvIudSmg44DGP9qaxhM55405hdHs9ru7WxNWjpjLhhCuVzV1V4a23kg9v5eem4ZNY2ncglx45LLF+PwzfWKseKeW1wLXgHu4mpdOnR5Y9//gL/lho45dP3dBp1m4Qkl/9+hge6i7Zf9YvOg5SZ048nbmDhnHExkW2Xf4+h918CSWhV23b3Y5qMnPi6R0DaeaE07h4zuXgSBocC0vTO+04vPktoXP2YT9n+HHfZai1hpOfW8WqdsHzuyf/gaWZGNLmV5NPZ3W3JsDfRwm4Tt6e3G6vjnRlXkqae8j31Cm/ZJ/hg9FO+RGW7VRZuZRXrjMnnMbsIWPYZtVShi1/l01yaznrxdvR7BJ6hVrpwe1HM3erYfz2qf/B0ARGsUBeN9GA3x74I+4cNoEFW+zAxXMuJ6ML9EIBslmKUjBz/IwO4XLz8En88PXZHXRvHj6J7047gK36dGeHzTbi3teXce3zH3V8j6UZmE6po04e33k/9t92ONM2tblpuWDUkgVc+tBfO+rJRuDoOgXNJFvKg9DJGxlMp9Qh5HQBpx+4Hb27Z/j1v97sEBCzh4zhP4OGsUX7Rb2ONtisL2vWFpEV9Vzm663fX8rvDh7D6XcuwMEVrn8fPZXNzpqB/GgJ7z/yHL964n862mXWt4/l1YFDWdbUj1FLFnSiddeu45j2+pyOuyiWZjBz8hks3XRLLLvzsPiwz0BuGj6J4yrq0tvn8rrJHcMmAHDYrpvz5boisxnD3EHDOi4jrunWhOZe6WD2kM7vAIasWMxG+VbWNXRnUb/BrOnWxIXfGcr+2w5nq7Vf8HHPTTn3B/sycNgWnLvZp5x970Ke2HUsowcNY8u1K1hrNtKjmOtED2B5z76dfh/QvILPem2G7UimvPlkp2+SuP6p7hx2MPfvtD9NdoGRHy/k+Ff/TbG97mZOOA2gYgxYCMfpFNfB1nWe33YPrI37YJZsbFsihFuXZ008w80bMPaLjd0RxQKabXeyGMobGUxdY83UY+l15y20oWHaJc7f/2Q+3bgfmhC8sek2VW19yT7HsKZbTxb0354dx+zBS4cPrYvQh6/BqkcIsTfwOynl+Pa/zwWQUl4UlKcWq54yVrUUWHXplQz+w6+QZgbdsRHXXw/T3CAH9zzyOnfc9TyfbdyPLxp7cvb47fnRtwcD8OFbH/HRvLcYvOw9trzw15Ax0SyL1rPO5qMjjiHf1Is3lzezSY8GdthsI/KfrWBA8wpk9x6s/HwVfTfrg/bpUvKPPcWmN1yDMAwoFmn96Wl8dOyP2HzbgR0N/MCCTznzrgWUHOjd1syWa1dwzPfH8K2h2zCgVyNrWot8ed3NjPj9TKySTYNtkTcyGLrGi7/4E59NOJwtejXSv6mRdz5fx4rFS9ku9yU7jdqV3lsPcCtj5Uq+WvQen9kGPb/4jLzt0Hv7wYjWFu7+0uCS+avI6DpF2+GnY7flqJFbuhejlizhq69aePXxV/hLcy8+3WwrLMfhsrFbcPBGBR5YvJabH32DFb03Y2VjT04YvTV7D+7DULNI7y8/gx49oKUFBg3iA6eBBUu/YtjAjVn02VquvuYhdvv8XeZvtj2nnnIIhw7boqr93lrejLl6FTsWVuNsuRVLjR50z+i0Fm0G9GqkT48st730Mb978E02aV3Ldp9/wISdN2fPww+gzbJZ+voiLnizDVMTbLrmc6ZO2ZdNBw8ABDv179nRBqtaCvxn8So+XtXCVn16sPfgPsD/be/uY6u66ziOvz+0wIQRKLRhPGnZAGdZlCcJGT4sssRtmrHE+DAfQoyZiY+bM1mmkeg/JMbghkZiXCa6ZIRhcMkWg5NlmyFLDEMYiXRlUAujpQ+UB8tgBfrw9Y972hUGFNZze+g9n9c/Pff03nO+v3PP/ZzTX889P9hW20rDsTN8tmYqS2ZPSWo6xanO84w5cYLK4y1UL67p3859v4dg/vSJFyy/pb6Ria1NdNw0kxumTaW5o5MjJzs5193DgrFdjGs9QuX8eURVFS31jVTs28vZnl4qbl9KVFXRdLKT8WPKaO44y8t1bTy94y1GjxpFTwSrbxnFfV1NdC/+OO1bX+SWNT9DY0YzqqeH1sfW8+qSO1kwaxJzpk4AoL7tbV6tP8YN5aOYXjGuv2uhtvkU//rvMZ7Y3kDfMWZ0mbh/6Syeee0whe+Fc7jAwQAABelJREFU97L2iwu4d8EMjp8+R9PJzv73YuB711fvmfM9jB9Txt/3tvK7V+rp7ull4PHrk3OnsGz25Au2e21zB2WbN7N0zSNQVk5593kOfO0BVk36BKcnVNDV28uvvvBRls+ppKW+sdB1Vzmd3Z3lVE8ZR1dbG2ObmqheXEPD5ueYv/phusvKKe/pZtfqtXzk4W/31w1cUGu0tdG4p44PHtrPbb/+BT3J60atexwWLeJE5TSaO85SsW8vx985z7lp05k7TkyqmQdVVf2ftbqxFXRNrhqwbTs41dlNQ91Btm/bSevkm2gdM4GVC6bznU/P6X9vhupyV/VkEfzlwH5gBXAE2Al8NSIue9PpNIK/X3s7HDoE1dWFN2aAy+24V/v6oa5/YB21zR1cHEgXL6cvvKeVdb+7o6XgarbDYB/yK27D97nONJaV5nquJ1ds1xD32UsdwNLYjn3L6Oru4dDxdy44GF3SRe14vzWcONjEsdr9VM6f9+7J0NUY6mf/Moq5T143wZ8Ucw+wjsJFGRsiYs2Vnp9q8JuZ5cTlgj+TPv6I2AoMPkimmZmlrqRv0mZmZu/l4DczyxkHv5lZzjj4zcxyxsFvZpYzDn4zs5wZEffjl9QOvDXExVQCx1IoZ6Rwe0ub21va0mrvhyLiPd82GxHBnwZJ/77UFxlKldtb2tze0lbs9rqrx8wsZxz8ZmY5k6fgfyLrAoaZ21va3N7SVtT25qaP38zMCvJ0xm9mZjj4zcxyJxfBL+kuSW9Kqpf0aNb1FJOkWZJekfSGpFpJD2Zd03CQVCbpdUl/G/zZI5ukSZK2SNonqS4Z1a5kSfpRsi/vlbRJ0qUH4B2hJG2QdFTS3gHzJkt6UdKB5GdFmuss+eCXVAasB+4GaoD7JdVkW1VRdQM/jogaYBnwvRJvb58HgbqsixgmvwFeiIhbgY9Rwu2WNAP4IbAkIm6jMHjTV7KtKnV/Bu66aN6jwEsRMRd4KXmcmpIPfmApUB8RDRFxHngGWJlxTUUTES0RsTuZfptCKMy48qtGNkkzgc8BT2ZdS7FJmgh8CvgjQEScj4j/ZVtV0ZUDH0iGbR0HNGdcT6oiYjtw4qLZK4GnkumngPvSXGcegn8G0DjgcRMlHoR9JFUDC4Ed2VZSdOuAR4DerAsZBrOBduBPSdfWk5LGZ11UsUTEEWAtcBhoAToiYlu2VQ2LqRHRkky3AlPTXHgegj+XJN0I/BV4KCJOZV1PsUj6PHA0InZlXcswKQcWAb+PiIXAGVLuBrieJH3bKykc8KYD4yV9PduqhlcUrrlP9br7PAT/EWDWgMczk3klS9JoCqG/MSKezbqeIlsO3CvpEIVuvM9IejrbkoqqCWiKiL6/4rZQOBCUqjuBgxHRHhFdwLPA7RnXNBzaJE0DSH4eTXPheQj+ncBcSbMljaHwj6HnM66paCSJQv9vXUQ8lnU9xRYRP4mImRFRTeG9fTkiSvaMMCJagUZJH05mrQDeyLCkYjsMLJM0Ltm3V1DC/8we4HlgVTK9CnguzYWXp7mw61FEdEv6PvAPClcEbIiI2ozLKqblwDeA/0jak8z7aURszbAmS9cPgI3JiUwD8M2M6ymaiNghaQuwm8IVa69TYrdvkLQJuAOolNQE/Bz4JfAXSd+icEv6L6W6Tt+ywcwsX/LQ1WNmZgM4+M3McsbBb2aWMw5+M7OccfCbmeWMg99sEJJOZ12DWZoc/GZmOePgN7tKku6Q9M8B98LfmHybFEn3JPN2SfptHsYFsJHLwW92bRYCD1EY2+FmYHkyMMgfgLsjYjFQlWF9ZoNy8Jtdm9cioikieoE9QDVwK9AQEQeT52zKqjizq+HgN7s25wZM95CD+11Z6XHwmw3dm8DNycA3AF/OrhSzwflsxWyIIqJT0neBFySdoXArcLPrlu/OaZYCSTdGxOnkKp/1wIGIeDzruswuxV09Zul4IBn/oBaYSOEqH7Prks/4zcxyxmf8ZmY54+A3M8sZB7+ZWc44+M3McsbBb2aWM/8Hn1fGrkMOm5wAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(1, 1)\n",
"polling_centers.plot.scatter(x='lng', y='lat', ax=ax)\n",
"people.head(10000).plot.scatter(x='lng', y='lat', color='red', ax=ax)"
]
},
{
"cell_type": "code",
"execution_count": 62,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 5.3 ms, sys: 117 µs, total: 5.42 ms\n",
"Wall time: 5.4 ms\n"
]
}
],
"source": [
"%%time\n",
"tree = KDTree(polling_centers)"
]
},
{
"cell_type": "code",
"execution_count": 63,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 6.64 s, sys: 123 ms, total: 6.76 s\n",
"Wall time: 6.8 s\n"
]
}
],
"source": [
"%%time\n",
"dists, i = tree.query(people, k=3, return_distance=True)"
]
},
{
"cell_type": "code",
"execution_count": 64,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"lat 8.960802\n",
"lng 5.894139\n",
"Name: 0, dtype: float64"
]
},
"execution_count": 64,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"people.iloc[0]"
]
},
{
"cell_type": "code",
"execution_count": 65,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>lat</th>\n",
" <th>lng</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>19201</th>\n",
" <td>8.952971</td>\n",
" <td>5.902439</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1215</th>\n",
" <td>8.933707</td>\n",
" <td>5.919021</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2580</th>\n",
" <td>9.013949</td>\n",
" <td>5.897545</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" lat lng\n",
"19201 8.952971 5.902439\n",
"1215 8.933707 5.919021\n",
"2580 9.013949 5.897545"
]
},
"execution_count": 65,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"polling_centers.iloc[i[0]]"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"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.7.6"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment