{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# MPSlib: Computation of entropy and self-information\n", "\n", "The self-information, and entropy (the average self-information), can be computed using MPSlib by setting\n", "\n", " do_entropy = 1\n", "\n", "This works for all algorithms except when using `mps_genesim` in direct sampling mode\n", "(i.e. when `O.par['n_max_cpdf_count'] = 1`).\n", "\n", "See details in:\n", "\n", "[Hansen, Thomas Mejer. \"Entropy and information content of geostatistical models.\" Mathematical Geosciences 53.1 (2021): 163-184](https://doi.org/10.1007/s11004-020-09876-z)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import mpslib as mps" ] }, { "cell_type": "markdown", "metadata": { "tags": [] }, "source": [ "## Setup MPSlib\n", "\n", "Setup MPSlib and select to compute entropy." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Initialize MPSlib using the mps_snesim_tree algorithm and a simulation grid of size [80, 70, 1]\n", "#O = mps.mpslib(method='mps_genesim', simulation_grid_size=[80, 70, 1], n_max_cpdf_count=30, verbose_level=-1)\n", "O = mps.mpslib(method='mps_snesim_tree', simulation_grid_size=[80, 70, 1], verbose_level=-1)\n", "O.delete_local_files()\n", "O.par['n_real'] = 1000\n", "O.par['n_cond'] = 9\n", "# Choose to compute entropy\n", "O.par['do_entropy'] = 1\n", "TI, TI_filename = mps.trainingimages.strebelle(di=4, coarse3d=1)\n", "O.ti = TI\n", "\n", "O_all = O.run_parallel()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot entropy" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'Entropy as a function of number of realizations')" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABBQAAAGDCAYAAACbREPAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAClk0lEQVR4nOzdeZxcdZ3v/9enqruzQEg6CUvIHsIiiQ6kA8Q9IKOiOCi44TI6M8ioOHOdGcdxtuhlrvPzOjrj3CtXB9FBR4kLoGCUUUHavUPoCJJAgCZJhyZhSdKBhCzdXfX9/XFOdU5X6lSdU3Wq6lT3+/l49CNd26nv+Z5zOufz+W7mnENEREREREREJI5MswsgIiIiIiIiIq1HCQURERERERERiU0JBRERERERERGJTQkFEREREREREYlNCQURERERERERiU0JBRERERERERGJTQkFEREREZGUM7MPmNlTZnbAzGY18Hv/zsxuaNT3Bb73TWb2uL+/5zb6+4vK4sxsaZO++0wzu8/M9pvZn9f5uxb5+9rmP77DzN5Th+/ZbGark96uNIcSCiI1MrPtZnbI/w+v8PP5iJ/tNrOr6l3GejKz5Wb2IzPbbWYu5D1vN7OHzOx5M3vMzF4eeO0qM+vz6+2/zezUwGufMLPhorpdEvIdF5rZA2a2z8z2mNl3zWxu8nssIiLVmOj/X9bCzNqBfwVe7Zw73jm3p07fs9rMBoLPOef+2TnXjLr/DPAhf39/24TvT4uPAnc756Y55/5PI7/YOXeJc+6rtWzDzG40s/9VtN1lzrnumgonqaGEgkgy3uD/h1f4+VASGy1kiFNuGPg28CelXjSz3wf+N/BHwDTgFcBW/7XVwD8DlwEzgW3A2qJNfKuobreGlONB4DXOuRnAqcCjwBeq3isREamHifz/ZS1OBiYDm5tdkAZayDjb3yrP08j1YB7Fd9JQOuFE6sjM3mtmvzSzz5jZoJltM7NL/Nc+Cbwc+HywlcbvanaNmT2KFxRjZu/zW/H3mtntRa34zsz+3My2+r0E/sXMMmbW4b//hYH3nmRmB83sxKT20Tn3sHPuy4T/Z/c/gWudcz3Oubxz7gnn3BP+a5cC33HObXbODQH/BLzCzE6rohxPOed2Bp7KAU3pnigiIvG00v+XZnaamf3U7w2328y+YWYzAq//jZk9YV4X9YfN7FUh+/x6M/utmT1nXtf+T4S87wzgYf/hPv+7x3RN99832oujXH36r880s/80s53+698zs+OAO4BTAz1ITjWvt+DXA5/9A/O6rO/zv/MFgde2m9lHzOx3ZvasmX3LzCaH7FfGzP7BzPrN7Gkz+5qZTTezSWZ2AMgC95vZYyGfd2b2fjN71C/LdWZm/mvFZS7uyt9tZv/LzH7t7+f3zWyWfyyfM7MNZrao6CtfV3zuBLb/x+b1xBw0r9fmwqJyjjlPS+xLyTo1s58CF3L03D+jxGe7zeyTZvYr4CCwxMzOMrOf+Of1w2b21sD7I513gW0Xzqn7bWzvImf+sAUz+46ZPekf85+b2TL/+auBdwIfLdSz//x2M7vY/32SmX3OPxd3+r9P8l9bbWYDZvZX/jmyy8z+KFC+15nZg+Zda0+Y2UfC9kXqRwkFkfq7AO9GYDbwaeDLZmbOub8HfsHR7nzBVpo3+p8728wuAv4/4K3AHKAf+GbRd7wJWAmswGvt/2M/QP8m8K7A+64E7nLOPVNcSDN7mf8fWdjPy+LuuJll/XKdaN4N3oCZfd7MpgTfVuL35YHn3uD/h7jZzD5Q4fsWmNk+4BDwEbz6FhGR1tAS/1/i/V/1/+H1hnsBMB/4BHjj3YEPAec556YBrwG2h+zv88AfAjOA1wMfMLM3Fr/JOfcIsMx/OMM5d1HI9oqVrE//tf8CpvrbPQn4N+fc88AlwM5AD5Jgor6Q3FgLfBg4Efgh8H0z6wi87a3Aa4HFwIuA94aU773+z4XAEuB44PPOuSPOueP99/yec65cI8OlwHn+97wVr76jejvwbmAucBrwG+A/8XpMPgR8vOj9x5w7AGZ2GfB3wOV4dfILju1t+Ub887S4EOXq1D/WwXP/kZB9eTdwNV5P0GeAnwA34R3btwP/z8wK3x3pvCvmnPu9wnkB/CXeubXRf/kO4HT/+zYC3/A/c73/+6f9z76hxKb/HlgFnAP8HnA+8A+B108BpuMdpz8BrjOzTv+1LwN/6l9ry4GfVtoPqQPnnH70o58afvBuFA4A+wI/7/Nfey/QF3jvVMABp/iPu4GrirbngIsCj7+M94e48Ph4vGEGiwLvf23g9Q/i3QSB95/XDsD8x/cCb61TPSz1/qSMee5Uv3z34t3czQZ+BXzSf/1iYDfejcAU4D+APHCl//rZ/jaywEuAXYXXKpRlJvA3wKpmnx/60Y9+9KMf72e8/n+JFyz+1v99KfC0//9be8z6+RxeYF/qtUV++dtKPS6uo3L16f9/nAc6S3zPamCg6LlPAF/3f/9H4NuB1zLAE8DqwDF+V+D1TwNfDNmnu4APBh6f6R+vwj46YGmZ+nLAywKPvw18rLjMIfXXDfx94PXPAncEHr8BuK/ou8LOnTuAPymqk4PAwlLnaYn9qFSno8c15PPdeD1BC4/fBvyi6D3/AXy80nkXUk/F193L8M7xM0K2N8PfxnT/8Y3A/yp6z3bgYv/3x4DXBV57DbA9cD4eYux5/jT+/R3eNfunwAlxrjX9JPujHgoiyXijc25G4OdLgdeeLPzinDvo/3o85T0e+P1UvFaWwjYOAHvwMrWl3t/vfwbn3Hq8/9RWm9lZeDc6t0fbpUQc8v/9v865Xc653XiTSr3OL9+deC0At+D957Id2A8M+K8/6Jzb6ZzLOed+Dfw78OZKX+qc2wt8FbjNxv+4WhGRVtLy/1+a2clm9k2/i/VzwNfxEuY45/rwWpo/ATztv+/UkO1cYGZ3m9kzZvYs8P7CdhISVp/zgb3OucEqtllcx3m8Og3W8ZOB3w8SfgzHbMv/vQ1vvoioon5XKU8Ffj9U4nHxtkqeO3hzHPx7oUcnsBevF0vYeVcsSp1WEtz+QuCCYC9TvGEHp0Bt552ZzcdL3LzH+b0lzCxrZp8yb9Lt5zjaIyfquVzqPAheM3uccyOBx8HjfAXePWW/mf3MzF4c8TslQUooiDRXyVURip7fifefAwDmjXGchZe9Lpgf+H2B/5mCr+J143w3cLNz7nCpLzSzlxeNjSv+eXmpz5Xj36wMFO2PK3rPdc65051zJ+MlFtqATWGbZOwQiXLa8LrenRCr0CIikkap+f8SbzJhB7zQOXeC/5nR/5ucczc5517ml8XhTUxcyk14SYv5zrnpwBeJ/n/c8/6/UwPPnRLxs48DMy0w70NAWD0XFNex4dXpE6GfiLgtvOMxwtjAvlrPU13dlBN27jyO1+0+mCib4jeEFJSr1yTqNLj9x4GfFZXneOdcYdhoVeedP1z1e8DnnHN3BF56B94QkIvxhiYsKnykRNlKKXUe7Ax57xjOuQ3Oucvw7ve+h5fskAZTQkGkuZ7CGzdYzlrgj8zsHH+Smn8G1jvntgfe89dm1ulnjv8H8K3Aa1/HG/f3LuBrYV/inPuFGzvzdvHPL0p9zjyTgQ7/8eTCZDq+/wT+zLwJrjqBvwDWBd673N/GAuB64N8LrSZmdpm/X2Zm5wN/DtwWUo7LzVurOWPeJFr/itcFdW/YPouISMtIzf+XeOPUDwDPmrc88V8XXvD/H7rI//7DeC3d+TLb2eucO+z/H/eOCvs3ynlzOzwBvMtvIf5jvHkAonx2F143/f/n10W7mb3Cf/kpYJaZTQ/5+LeB15vZq8xbyvKvgCPAr0PeX85a4C/MbLGZHY93vL5V1BpdrfvwJnle4O/L3yawzbBz54vA39rRiQinm9lbYmw3yToF7x7rDDN7t39s283sPDs6eWa1591XgC3OueL5qab55d2Dl8T556LXK127a4F/MLMTzWw2sAbvWizLvMlU32lm051zw8BzhF9rUkdKKIgk4/tFrfnfjfi5fwfebN6swCXXFvaHBfwjXuv9LrwbhrcXve02oBfvP9Af4I0jLXz+cbwJchzexD5JW4h3w1RY5eEQR2ejBm/lhg3AI3iTHP0W+KT/2mS8TPkB4B68CZH+MfDZtwN9eMMgvgb8bxdYD7mo58Rc4L/99z6A95/KmxLZQxERScp4+P/yf+JNzPesv41bA69NAj6FNz/Qk3gtp2HB7AeBa81sP14QFbd19X14yYw9eJMrxglA3403X8EWvDHpHwZwzm3BC/C2+t3lxwzXcM49jJdw+b94+/gGvKVAh2KWHbwA9b+An+MtG30Y+LMqtnMM59xP8AL+3+Ed73UJbLbkueOc+y5eL5Rv+l3+N+FNbhm1rEnWKc65/cCr8c79nXjn4f/GOzeh+vPu7cCbSvRe/RreMIUn8Jbw7in63JfxJk3dZ2bfK7Hd/4U3Z8nv8O7fNvrPRfFuYLtf7+/HG9ohDVaYeEZEWpSZOeB0f9xm2Hu+gjdr8z+EvUdERGQ80/+XIiLJ02RlIuOceesoXw6c2+SiiIiIpJb+vxQRiU9DHkTGMTP7J7yud//inNvW7PKIiIikkf6/FBGpjoY8iIiIiIiIiEhs6qEgIiIiIiIiIrEpoSAiIiIiIiIisaViUsbZs2e7RYsWNbsYIiIiqdPb27vbOXdis8sxEdTjfuT555/nuOOOS3SbE5XqMjmqy+SoLpOjukxOPeoy7H4kFQmFRYsWce+99za7GCIiIqljZv3NLsNEUY/7ke7ublavXp3oNicq1WVyVJfJUV0mR3WZnHrUZdj9iIY8iIiIiIiIiEhsSiiIiIiIiIiISGxKKIiIiIiIiIhIbBUTCmY238zuNrMHzWyzmf0P//lPmNkTZnaf//O6wGf+1sz6zOxhM3tNPXdARERERERERBovyqSMI8BfOec2mtk0oNfMfuK/9m/Ouc8E32xmZwNvB5YBpwJ3mtkZzrlckgUXERERERERkeap2EPBObfLObfR/30/8BAwt8xHLgO+6Zw74pzbBvQB5ydRWBERERERERFJh1hzKJjZIuBcYL3/1IfM7Hdm9hUz6/Sfmws8HvjYAOUTECIiIiIiIiLSYiInFMzseOAW4MPOueeALwCnAecAu4DPxvliM7vazO41s3ufeeaZOB8VERERERERkSaLlFAws3a8ZMI3nHO3AjjnnnLO5ZxzeeBLHB3W8AQwP/Dxef5zYzjnrnfOrXTOrTzxxBNr2QcRERERERERabAoqzwY8GXgIefcvwaenxN425uATf7vtwNvN7NJZrYYOB24J7kii4iIiIiIiEizRVnl4aXAu4EHzOw+/7m/A640s3MAB2wH/hTAObfZzL4NPIi3QsQ1WuFBRERERCQd+gZzbL67j1VLZtG1sLPyB0REQlRMKDjnfglYiZd+WOYznwQ+WUO5REREREQkYb39g3x6w2FG3MN0tGX4xlWrlFQQkapF6aEgIhPIoo/9oNlFiGT7p17f7CKIiIi0nJ6texjOe12Mh0fy9Gzdo4SCiFQt1rKRIiIiIiLSulYtmUV7BrIG7W0ZVi2Z1ewiiUgLU0JBREREUsPMvmJmT5vZpsBz55hZj5nd5y85fb7/vJnZ/zGzPjP7nZmtCHzmPWb2qP/znmbsi0gadS3s5KPnTeYvX32mhjuISM005EFERETS5Ebg88DXAs99Gvifzrk7zOx1/uPVwCV4q0mdDlwAfAG4wMxmAh8HVuL17O41s9udc4ON2gmRNFvamWX16qXNLoaIjAPqoSAiIiKp4Zz7ObC3+GngBP/36cBO//fLgK85Tw8ww1/W+jXAT5xze/0kwk+A19a/9CIiIhOLeiiIiIhI2n0Y+JGZfQavMeQl/vNzgccD7xvwnwt7/hhmdjVwNcDJJ59Md3d3kuXmwIEDiW9zolJdJkd1mRzVZXJUl8lpZF0qoSAiIiJp9wHgL5xzt5jZW4EvAxcnsWHn3PXA9QArV650q1evTmKzo7q7u0l6mxOV6jI5qsvkqC6To7pMTiPrUkMeREREJO3eA9zq//4d4Hz/9yeA+YH3zfOfC3tepGa9/YNcd3cfvf2akkNERAkFERERSbudwCv93y8CHvV/vx34Q3+1h1XAs865XcCPgFebWaeZdQKv9p8TqUlv/yDvvKGHz/74Yd55Q4+SCimn5I+kwXg/DzXkQURERFLDzNbireAw28wG8FZreB/w72bWBhzGn/MA+CHwOqAPOAj8EYBzbq+Z/ROwwX/ftc654okeRWLr2bqHoZE8eQfDI3l6tu7RsospVUj+DI3k6WjLaIlMaYqJcB4qoSAiIiKp4Zy7MuSlrhLvdcA1Idv5CvCVBIsmMfX2D9KzdQ+rlswaNzfQq5bMoqMtw/BInva2DKuWzGp2kSSEkj+SBhPhPFRCQUREREQSNV5b5boWdvKNq1aNu0TJeKTkj6TBRDgPlVAQERERkUSN51a5roWd42ZfxjMlfyQNJsJ5qISCiIiIiCRqIrTKSfop+SNpMN7PQyUURERERCRRE6FVTkRElFAQERERkToY761yIiICmWYXQERERERERERajxIKIiIiIiIiIhKbEgoiIiIiIiIiEpsSCiIiIiITQG//INfd3Udv/2CziyIiIuOEJmUUERERaRG9/YNVrZzQ2z/IO2/oYWgkT0dbhm9ctUoTJoqISM2UUBARERFpAbUkBXq27mFoJE/ewfBInp6te5RQEBGRmmnIg4iIiEgLKJUUiGrVkll0tGXIGrS3ZVi1ZFYdSyoiIhOFeiiIiIiItIBCUmB4JB87KdC1sJNvXLWqquESIiIiYZRQEBEREWkBtSYFuhZ2KpHQINXOdSGNp2MlUhslFERERERaRCHgKQx3UACUPpoAs3XoWInUTgkFERERkRahACj9NAFm69CxEqmdJmUUERERaRG1TMwojaEJMFuHjpVI7dRDQURERKRF1DIxozSGJsBsHTpWIrVTQkFERESkRSgAag0TeQLMVpvkcCIfK5EkKKEgIiIi0kIUAElaaY4PGW9aLUHWDEooiIiIiMi4pGCgsTTJYWvSdVKaEmTRKKEgIiIiIpG1SvAx0YKBNBwXzfHReibadRKHEmTRKKEgIiIiIpG0UvAxkYKBtBwXzfHReibSdRKXEmTRKKEgIiIiIpFauFsp+JhIwUCajksz5vhIQ++MRkl6XyfSdRKXEmTRKKEgIiIiMsFFbeFupeBjIgUDrXRckpaW3hmNUI99nUjXSTXSPAluWhJpSiiINNCij/2g2UUQERE5RtQW7lYLPtIcDCSp1Y5LktLUO6Peqt3XSoHnRLlOxpM0JdKUUBARERGZ4OK0cCv4SKeJelwmUu+MavY1TYGnJCdNiTQlFEREREQmuIncwi3pUU0X7ol07lazr2kKPCU5aUqkKaEgIiIiIhO2hVvSoZaW9PF07iY9PCFNgackJ02JNCUURERERKQlpWVSslqNl/2ohVrSx8ekizqXJx4lFERERESk5YyXseHjZT9qpZb0+iVVGtWDo5XO5VZPfKSprpVQEBEREZGWM15atMfLftQqTV24m6XVkyqtci6nKRivVprqWgkFEREREWmKWloJWz34Kmjl/Ui6lXc8zYVQjVZPqrTKuZymYLxaaaprJRREREREpOFqbSVs9eCroFX3Yzy08qZRKydVWuVcTlMwXq001bUSCiIiIiLScEm0ErZy8BXUivsxHlp5JXmtcC6HBeOtNq9CWupaCQURERERabjx0Eo4ken4SSsrDsbV46Z6SiiIiIiISMOlqcuuxKfjJ+OJetxUTwkFEREREWmKtHTZTUKrdZdOwng6fo0yEc+TVqAeN9VTQkFEREREJqSkgjt1l5YodJ6kl3rcVK9iQsHM5gNfA04GHHC9c+7fzWwm8C1gEbAdeKtzbtDMDPh34HXAQeC9zrmN9Sm+iIiIiEh8SQZ36i4tUeg8STf1uKlOJsJ7RoC/cs6dDawCrjGzs4GPAXc5504H7vIfA1wCnO7/XA18IfFSi4iIiIjUoFRwV61Cd+mskZru0r39g1x3dx+9/YPNLsq4Uku9pvE8EalVxR4KzrldwC7/9/1m9hAwF7gMWO2/7atAN/A3/vNfc845oMfMZpjZHH87IiIiIiJNl+SY6bR1l1bX+vooVa9xpO08EUlCrDkUzGwRcC6wHjg5kCR4Em9IBHjJhscDHxvwn1NCQUQSs+hjP2h2ESra/qnXN7sIIiIVpX2SuEL5Ju3LjbZkJSHp4C5N3aUnUtf6Rp6/pep1mcUrX5rOExkr7X8L0ypyQsHMjgduAT7snHvOmyrB45xzZubifLGZXY03JIIFCxbE+aiIiIiIJCDtLdnB8rUZnLtiMNHyjdfgbqLMWF/r+Rs3gCxVr/u3DdStfFEoCE5GUsdqIh6PSAkFM2vHSyZ8wzl3q//0U4WhDGY2B3jaf/4JYH7g4/P858Zwzl0PXA+wcuXKWMkIEREREald2luyg+UbcVRVvol4gz9RutbXcv5WE0CWqtfubfUpX732QUpL4lhN1OMRZZUHA74MPOSc+9fAS7cD7wE+5f97W+D5D5nZN4ELgGc1f4KIiIhI+qS9JTtYvqwRu3wT9QYf6tf7Ik0JmlrO32oDyDj1Wu/rK+0JwVaSxLGaqMcjSg+FlwLvBh4ws/v85/4OL5HwbTP7E6AfeKv/2g/xlozsw1s28o+SLLCIiIhIq+gbzLH57r6GB19Rg760t2QHyzdpX3/s8jXzBj9NgXdSoiZo6r3vwe1Xe/42IplW7+sr7QnBVpLEsZqoxyPKKg+/BMKmG3lVifc74JoayyUiIiLS0nr7B/n0hsOMuIcb2joet1U+7fMIFMrX3R0+Vj1Mo27wiwPoVu0ZUSkRECVBk+S+lypPqe1fc+HS2NtuVDKtntdXpX1IOrETZXutnEir9VilPUFbL7FWeRARERGRaHq27mE4D47Gto6nodttWCBY64123G004ga/VICbhmMQV5REQJQETVL7HlaeJOs27cm0KML2IemkVpTttWoiLUnj4ZyKSwkFERERkTpYtWQW7RnIORra/bURrfLlAvtSQQUQK9Aotf2+wRyfuSt+sFLvG/xSAW4rdn2OEqhHSdAkte9h5WnFuq1WLUm4pJNaUbbXiom0arRyL4x6UEJBREREpA66Fnby0fMmc2TGwobeeNa7Vb5SK2SpoAKIHGiEbX/L3lwqg5VSAW6zuz5XE/BEDdQrJWiS2vew8jS7bhul1tb+pBMvUbY3EZI96oVxLCUUREREROpkaWeW1avjj++uVT1b5Su1QoYFFVEDjbDtnzUzS0dbLnXBSliA26yuz9UGPEkG6tXue3EiJKw8E6Fbea2t/UknXqJsbyIkeyodl2b3Xih8/6R9OVY36DuVUBARERGRyCq1QoYFFVEDjbDtL+3MNnQCujjSFOD2bN3DkeE8DhgajheINnM/whIhaanXMPVaySWJ1v6k66+wvd7+Qa4L2edWOGa1KHdcmt17Ifj9bQbnrhhsyPcroSAiIiKpYWZfAS4FnnbOLfef+xZwpv+WGcA+59w5ZrYIeAh42H+txzn3fv8zXcCNwBS8Ja3/h78SldQoaktltYFGlJbp4oAm7o18s5IPjfjezqkdFE70vP84zQp1snPfoVQOaSmnniu5NLO1P+4cKWk/Tkkqd1yaPYdE8PtHHA37fiUUREREJE1uBD4PfK3whHPubYXfzeyzwLOB9z/mnDunxHa+ALwPWI+XUHgtcEfyxZ2Y6t0KWW77ta6q0KyAqFHfO3hwiIxB3kHGvMfVaETyY0yLasZoy2bI5ZIf0lKvfan3Si7NaO2vZo6UsFUmxuvQh7Dj0uw5JILfnzUa9v1KKIiIiEhqOOd+7vc8OIaZGfBW4KJy2zCzOcAJzrke//HXgDeihMK4UOuqCs1qRWzU9yYR1DQq+RGsk1ze8bbz5zN3xpREg9B67kuzVnKpp2rnSAkaj70YoiRIqulVkmTiJfj9k/b1N6zOlVAQERGRVvFy4Cnn3KOB5xab2W+B54B/cM79ApgLDATeM+A/JykX5ea61lUVmtWK2KjvTaKrfLOSH1esmJf499RzX7oWNmcll3qqdo6UoEadP43qBREnQRKnV0k9Ei+F7+/uHqj85oRYGoYTrly50t17773NLoZI3S362A+aXQRpoO2fen2ziyDjgJn1OudWNrscjeT3UFhXmEMh8PwXgD7n3Gf9x5OA451ze/w5E74HLAPOAD7lnLvYf9/Lgb9xzl1a4ruuBq4GOPnkk7u++c1vJrovBw4c4Pjjj090m+NV32COT284zHAe2jPw0fMms7QzO/p6sC77BnNs2ZvjrJnZMe+J8121fL5azfreYpXOy8KxGMlDW4ljkaR610m992U8XuNJXF9hdV5u23HqstLfiySte2yIWx4dxgEZ4PLT27n0tNrnJqnXdqE+5+WFF15Y8n5EPRREREQk9cysDbgc6Co855w7Ahzxf+81s8fwkglPAPMCH5/nP3cM59z1wPXgNXCsXr060XJ3d3eT9DbHq8139zHiHsbhdSE/MmPhmCU3g3W5usrvGG3RXDGLqyq0BPb2D3LLxgEMuDyhlvPVNW8hGZXOy9V4M8RXav1NooU4vBTJWE20falW2q7xuMek1PtX11iG1ZSu897+QT5zV6FFPndMi3ycuqz09yJJ0xYPsm57z2ivjSsvPi+R86he24XGnpdKKIiIiEgruBjY4pwb7cdpZicCe51zOTNbApwObHXO7TWz58xsFd6kjH8I/N+mlFoiq/eQgDjdi3v7B7ny+t8wlPN68n6nd4C172v9ceBxVOq6nbZx8qUC4+Bz11xYn2AzabUkaapZ7aRex7DU+ZPkUIhGDl0qNcwjiWRaM1fySJISCiIiIpIaZrYWr4FrtpkNAB93zn0ZeDuwtujtrwCuNbNhvBXy3u+c2+u/9kGOLht5B5qQMfXqfXMdJ5jp2bqH4dzRYcGNmryxlWbGb/YSeUGlAmMgVQmPKGoN8OMek6SOYdTzNuqEjlG21ehgPJggSTIR04yVPJKmhIKIiIikhnPuypDn31viuVuAW0Lefy+wvNRrkl6lbq4LAcakfbmaumLHadFctWQW7Vkb7aHQiMkb09biX0mzl8gLKhUYA6lJeERVa4Af95g0ekWQSkmAvsFcYEhE5Wug0cF44W/Rzn2HqjpOrZQwjEMJBREREREJ1cyb4GCw0mbeuOxGdC/uWtjJ2qtfnPgcCuWkqcU/ijR11w4LjNOS8CiodC3VGuDHPSbNWBGkXBJgy95caq+BMX+LMkZbNkMuF/04tVrCMA4lFERERESkpGbfBAeDlRFHrACjVPAWp0Wz0a2faWrxjyot3bXDAuO0JDwg2rWURIAf95jUegzLnbdxk5FnzczS0Zar6zVQbYI0+Lcol3e87fz5zJ0xJfJ2Wi1hGIcSCiIiIiLjQD16EjT7JjgYrGSNyAFGsxMh1UhTi3+axBlTXypAT0s9Rr2W0lTmKMLO22quwaWd2bpeA7X8XShOnFwRs9dSKyYMo1JCQURERKTF1SuAbvZNcDBYmbSvP/I+NTsRUq1WCyZrESVRUMt5nbbx6s2+luopyRUd6nkN1PJ3odaE33hOGCqhICIiItLi6hFAFwKyNZcuY/DgUNU3wbUGdoUAo7t7YMzz5bY7noO3WvX2D7LusSGmLa5+PookylApUdDbP8jn7nxk9LweGs7zuTsf4ZLlcyqej7Um2OqRjBjPAWUpabwGk5ijopbjNl4ThkooiIiIiLS4pG/ek+rxUK+eE5W2m5bgLW2t5IV6OzKcZ932nqYNBamUAAuW0wGGty7sLx/dzS8e3U3GKHs+1ZJgq+dwmfEaUJaSlmuwnmVK2/XdLEooiIiIiLS4pG+Uk+rxUK+hB1G22+zgLY3zOBTqzQFHhvPcsnGgKWWqlAALltOAhbOmsmPvQfLeKp4Vz6daEmytOlwmjZp9DZaSVJnSeH03ixIKIiIiIuNAkjfvSfV4qFe35zR2py6WxsB01ZJZtGWMoZzDATf3DsSeXC4JlRJgnVM7RpMHDnjtslO48TfbGRrOkwcyRtnjXkuCrRXOrfEuTS3/YWVJ4/XdLEooiIiIiMgYSfV4qEe357hzOzQrOGl2YBq2bOZbVs7nG+t3AJDLNSYQinsMBg8OYXjJhAwwbUr76HnUObUj0nGvNsGWxq76UaQpCK9Fmlr+y5UlyvU9Xo5JJUooiIiIiMgxkurxkGTPibjBRjODk2YGpsX7HUy+XL5iHt/ZsIOcK9/KX6+yfOOqVQBlj8uqJbOY1D42WEuyq3qlY5LGrvoQXvY0BeG1SlPLf7myVLq+x9MxqUQJBRERERFpCXGDjWYHJ3ED06RaNIP7PTScZ81tm8g7NxrYfPS8yRyZsbAhiY5SxwAoe1zqlYxJS5BXzXEuV/Zmn+dJalTPnijHoFJZyl3ftU4M2ko9G5RQEBEREWmSVrtxjKMe+xY32Gj2sIM4kgx2g/ttZuSdGxPYLOvMsnr10oT3oHJZshlj575DLDt1esXjUo9eAo0MvJPuTVCu7GkKwmvZ9rrHhrhyMXXv2RP1GDRjLo60JL3iUEJBREREpAla8cYxqnrtW9wb/FYaD19NsBsW4AX3u3NqB5/4/mYvoM96gc3+bQP13p1jynLLxgFu7h1g7T07jhmG0ajj0sjAO+neBOXK3ojzvJ5/r0otZ3rNhfVLeMU5BsWJrahJleJrsNAzp1KdtWJvEyUURERERJqgETeOzeoBUc99i9tyndbx8MXiBrvBAK8tY7xl5XwuD6zYUNjvm9bvIJfLex9yrt67UVLXwk56tu5hJHf0nBg8OFTXoDGsHFED71qunXr0JqhU9mDCIvg4KVGu6WrrLLhMaCOC6ErHIKneJYXX4nymuGydUzu47u6+1E4+C0ooiIiIiDRU4cavc2pHXVtL692iWO7mNWrQlNRNcCNvpuv1XXFbmcfMk5Bz3LR+B7dsHBhznHv7B1lz2yZyfh5hJO+8IQ+WWLEjS8vwkygJplqvnXr1JihX9nr3eIoShFf7/YVtDw035twodwyS7l0S9zPFPRuuXbe5Yp2GTXzaKEooiIiIiDRIudn30zZevJZWuihBU9wApG8wx+YSLXWNHDpS7++K0726EIQdGfZadku17vZs3UMuf7RXQsas4UMeCqLMip+WoSm1XjtRehOk4XqPU+eV9qmWOitse+2dG7jy4vMacvzDjkHSvUuq+UyhbH//3QdGr+9ydVqqzI1MGiqhICIiItIgxTd+9ez2XUuLcBKtdJWCpjgBSG//IJ/ecJgR93BTZ7hv9KR+5ZIXxXMT5HLHHufC8otDw3kyGePay5bTtbCT7m11KXJFYedEpeEbjZZEb4pGD7WpZchM1ORYuX2qtc66Fnay/7SOROqsUqIkSqIuqd4l1fZI6e0f5Dv3Pk4hHWgZC63TVUtm0ZYxhnOObKbxSUMlFEREREQapJHdvmvpWt2IGeXjbKdn6x6G86Vb4UutItDbP1hzYFIq6Gjk8YuSvAg+NjgmCK/UtTuNPQLChm80Ui3XTlRJ139xV/lKcykknRxr1PwUlVRKlERN1CXZu6Saz/Rs3cNw7mjvoorTn5gBzv+3sZRQEBEREWmQRgQqxd9XzXfEbaWLEyAE3xs1AFq1ZBbtGcg5QssTXEWg1mA0LOho5PGLkrwoLuflK+Yd855S50DfYI7P3JX80I1qA8UowzcaLckeBsX1ktTQmeLtFrYRZdv1SI5FqbN6Dxu6deNA2WECURN1zU6yrVoyi2zGGPGHLDnnyg55GMl5+5zLaciDiIiIyLiWhpvVSuK00sUJEEq9d9WSWZHmZPjoeZM5MmNhaHmKVxEovvmOE+yWCzoadfyiJC+qbWXesjeX+NCNWgLFKMM3WlWpeok73CfOXCZxhiQ1Y0nVYPmODOe5ZeNAYr0ZiocJFJZJDUrL5KCVdC3s5NrLlrPmtk3knaOjTFlL7ZOGPIiIiIhIIqptNY4aOMcJjkq9F4j0+aWdWVavDp9volygEDfYDdtWo4cJVDoG1QZHZ83M0tGWSzSoSmIiw66FnVyxYl5qhmIkoVS9xFkFJe5cJpW2fdP6HdyxaReXLJ/DOy5Y0PA6Loz3H8o5HHBz7wBXhMyXEfe67dm6Z7RF34A3dx273WYlUqrxjgsWcOYp0yqWtdQ+NXKeFCUURERERFpANcFsI1ZAiBPUlprvYNmp0xNpMSwXKNSydFvS3dSTVG1wtLQzm3hQlVTLb9p68NSaRCpVL1GPWzVzmZTb9k3rd/B3330AgF88uhvwgtZ6K67Dt6ycz03rd4zpop9ED5ziOrmixBAgSN85Vk7UsjZzn5RQEBEREUm5uMMKCjfvjViVIE5QW2q+gySXzwy7qa5l6baCRq7wEFQpoK02kEg6AGmllt+okkgihdVLlPqvdsWBsG3fsWnXMY/rnVAoVYeXr5jHLRsHKl6Pca/buOdgmiYmbWVKKIiIiIikXNRgtvjmfc2lyxoyXjhOcFpqvoN6Lp9Z+M5ag91mjL1OY6+IcurVStqswK+aJFKpstaS9ElyxYFLls8Z7ZlQeFxvperwmguXRroeq7luo9ZJM6+t8ZbIUEJBREREJOWiBrPFN++DB4dS2Wpcab6DepS31mC30S3wvf2DfO7OR5rSKyJNblq/Y8zEdFF75zRjGEc1QWq9eqCUUuiNcMemXSybcwKDB4cSWWK1nHJDMwpDia67uy+0buuVpGpmj6NWShJGoYSCiIiISMpFDWaLb94LyzHWKwCuZcLHUvsTvNluyxhvWTmfy0MmbGuGRo1TLtRDYfm7jB27XGa12+3ZuodJ+3KsTqSk0b+zmnOwt3+QNbdtGp1obyhG75wkArVy52nhuaC4QWozgsvCRH+N+N5CPYUNaSre/4+s6Ej0vCx37jVrtYdmJTLqSQkFERERkRYQJZgNBkCdUzu4dt3mxIOGwk16rdsvtT/Bm+2hnOOm9Tu4ZeMA37hqVc3lbiW3bhw4mkwAXrp0Nh+++Iyajt+YZI3BuSvq2zJd/J1xzpHCObZz3yFyeTf6fMbsmMAv+N4khicUKz5PywXBcYPUZgWXY66z4Tyfu/ORms+vYlGOffH+b9mbK7u9OImpSt/frDk/WmXZyjiUUBAREREZRwoB0HV39yUerARv0jNm5J1LdPuFm+1CMO04uu1lVtOmW0Zv/yDfufdxCmF0W1smkWAvGLyNOBoSvIYtE1ouiCvupdLelmFkJE8mY1x72fLRbvLFSa22jNGWzZDL1W94Qql9CgbBcYPUOMFlksM5Ct87NJwnD/yqbzcbtu9NtKdClGRJ8f6fNTNbclvVHKso3194XDgvk0y4lhvCksZhaLVQQkFERERkHKpHS1jwJh3nyGQMwyW2/cLNdmEViGBwuH/bQM3bbwU9W/eMdvE34M1dyQz5GLNkp1HV8Yob1JYaghOn1TqXd7zt/PnMnTGl5PKdwaRWLu+46AUncXg4xyXL5yQScEbZp+IgOO4EpVGCy6SHRhS+93N3PsKv+naH9lSoJYkR5e9P8f7v33Z/yW1Vc6yifH/S9Rp1e40aOtUoSiiIiIiIjEP1aAkrvkkPjo0GxkyuFhxjXs0s7VesmDfmc93bai5+XdR7IsArVsxLoJRjz4dJ+/qrms8gbvBVfA7GabUeGs5jZiw/dfqYpQ3DklrZjPGzR55hJJdnw/a9nHnKtLKBcbUJt6hBcPE8C7Ws1FCPoRFdCzv58MVnsGH73pI9FYCagu2of3+CvQTC5vYolZgq/K0pt91K3590vdZzCEuaV4ZQQkFERERknEq6JSzqZIqYMZIb+3vcoKQVWvEaORFgEgp12t0dv7dHtcFS8XEsFcQXB0trLl02urLDtes2j0kOhCW1du47xNp7dhxTvrBjVEs9B/epVKIr7Hqo9hyp17j7Qh1c+/3N3D/w7DFDU2oNjqNcw739g1z5pR6v50ym9Nwe1c4NU+n7o/ZiqMcQljjbT/vKEEooiIiIiEwAwXHnpWZcj6rSZIrDOW/2A0fR7xWCkuIb6zg38s1qvatXi2StyZR61EetQW3YjP+lgqXBg0Oh83OUS2rdsnHgmPKVO0b1SlqFXg81BOb1HHf/0K7nRn/PZo/WXSMmD7x14wBDI3kARvLe43LDBJKcG6ZSvZYL5EtdY8WJj3JzM9y0fgff2rCDB3c9Ry5ffknUnq17RueVGRoeu89p6LmghIKIiIjIOFe8DCFAW9Z4W0LLMo4Zn++3yOZyY38vFZSErRix5tJlx7RCVtq34vc24iY7jTO216s1s1wgXylRVa5MpQL+SvVaKhEQVr5mHKOw66GW769n8iNszo5GTB7oKjwulvTxDKvX3v5BPnfnIyWTF+XO58K/5a7Bm9bv4O+++8CY7yuXHOmc2jFaL3n/caGMaei5oISCiIiISBM0smWpELQFb9ZHco5vBJZljFKGsDIXB3OF7yw3frzcihF3bNp1zI182CoPxQHpLRsHRls9y91kJ1X/l6+Yh/n/pqEbcj3HcRcHX8WJqoxRss7LlalUgFhti3ycRENQ0tdiueshDedIULk5Oxox7Gj5qdPJGuQdZI2Kc4bUu7cGlD6vo/Z6ifL6HZt2HfOd5ZIjgweHyPh1lDHvcZTvaRQlFEREREQaoHiStka2LBWChsPD+WNei3ojGmVd9+LHpX4vKDW5Hs5hZiybcwIbtu8dE2SWWuWht3+QJ/YdGrNcoFF57HcSLXvF27g8ockTa5VUC25YkB18vjhRFVbn5coUFiAmGcyWSoTU+1osdz2EqTaxUUtCpBEBepje/kGuXbfZSyZkjHed1Z6KhEvwvM4AL106e8zqF6XO5+AxqHQNXrJ8Dr94dPfo41effTJ/+srTAEpONhm2vbT0kFJCQURERKTOioPPK1bMa2jLUiFouGXjAN++93FGcl4IWNzyVihrqeAi6daw4pvh9754ETf8cht557jxN9uPGWtfPPld8cR3bz9/wWhQX2osfVCpfSk8HzWoSro+kmolTyJAvGn9jtFJEYNBdvF5vObSZaMrMuTxuswHx+BHLVMjWsILmn0tRi1XpfH6UT4XVRL1X835GwzcnXMcGHYVt9WIbv7Ff5uCyQQo3QOluEzlzvfCqiV3bNrFJcvn8I4LFlQcRlH4+x3sqNXMZFCQEgoiIiIidVYcfDoqT3hWj27YXQuPLsdYasx7uZvaeoxdLl5SMDjsYfDgENdcuDT088E6zeW9Se8K26p0k11qGbq4QUqS9ZF0kFRLgNjbP8ia2zaNjqsfCgTZxefx4MGhMYmqXM6BKz0Kvh5Jg7jXSKlx8VGuxXLbW/fYENMWH7syQVxhCapK50Yaur1Xe/4WX0NnzczWZX/jnielEgbFPQeC53OpiSKvuXBp2e96xwULQpdDDduvwlCu4DC1RibjwlRMKJjZV4BLgaedc8v95z4BvA94xn/b3znnfui/9rfAnwA54M+dcz+qQ7lFREREWkbxjfPyU6cDhI69r2crXLkb0Eoz4ifdGlZcljiBXfHEdzf3DoxZmq9cMqJUMiNukJJkfaQhKAyWJZc/mhTImI126S4eXlLY756te8j7SZ1c3jWk/GG9KMKEjYu/YsW80SRbnOMY3N667T01X6NhCapK50Yaur2XS4aUq9fia2j/tvsT399q/5YW/jZF+XyUMlWqi0rbSNPfiGJReijcCHwe+FrR8//mnPtM8AkzOxt4O7AMOBW408zOcM7lEiiriIiISEvpG8yx2W/ZCltHPTj2vnDTuXPfoabcPFYzu345cVoG4wbowffv3HeItffsiBXQ1JLMCNtGteo970Hcskxq94YxZDLGtZctBxgzvORVLziZE6dNGv1M59QOMmaAS6S3RqV9KNeLIkylcfFx6yu4vSSu0bDzP8o12chu76WOT9icAlEC+eA11L0t+f2tNRCP8vlKZYpSF5W2kYbEUZiKCQXn3M/NbFHE7V0GfNM5dwTYZmZ9wPnAb6ovooiIiEjr6e0f5NMbDjPiHh7Tah62jnpv/yBXfqnHa3HP2jEtwUmXrdSNa5LBSTUtg3ED9GArYvG8CXG+v9x+12s1juLt1lrvSfVqKVWW4Dk7knP8dMvT5J3jlo0Do3Nf5PKObMZYc+mymoZbRNmHsF4U5VQaFx9XYXtDw2Ov0VonSIwbaIZ9rh7Cjk+lcyZOIJ/0/tYaiEf9fLU9v6Juo9GJozhqmUPhQ2b2h8C9wF855waBuUBP4D0D/nPHMLOrgasBFixYUOotIiIiIi2rZ+sehvMc04IZdoNaGB8LXtD26rNP4vfmz0h8xvdKXcWj3KyXWwEgOHFY0r0s4iRC4gY0pfa71CSExfNOVLsfYYFZtZKYaLIgWJbe/kHue3wfZkYGN2aJz6HhPNf/YiuF2D7v3OiSdkntQ1j38EIvCjO46KyTIu1TkgFZYXtr79zAlRefF7l7fBTF53kwcVZqFYBGqTQkKqmhGEkmSGo97kmcN0kMiSiUJU2JhIJqEwpfAP4J7//IfwI+C/xxnA04564HrgdYuXJl6dlbREREpCWZ2RnAXwMLCdxvOOcualqhGmzVklm0ZyDnxq6kEHaDWnwzNHvapLLzAJRSKaCppqt41O/o7R/kyut/w5C/gkRbjb0sim+w+wZzfOauystWFoKuzqkdNXcRDgZQQ8P5WGP2o243qWRLEhNNFis+phmDFQtmcN/As+Ryecwsdk+B4u0Hj3Gc1uDCZJA39w5w50NP8fNHn6m4j0kHZF0LO9l/WsfoNpM4ruWur0YuNVtKnCRB0r2datlOrcc9ic/XOiQizapKKDjnnir8bmZfAtb5D58A5gfeOs9/TkRERCaW7wBfBL6EN1HzhNO1sJOPnjeZrZw0ZqmvwmvFN4xXrJjHzfc+znDO0Z41rgjMrRBVpYCmmq7iUb+jZ+sehnNHt53LOd56wXzmzphSVS+L4hvsLXtzFYO1pHsUBAMoC7TM15oESGI8dKlW7FonmixWfEzzDjZsH6Q96y3TuezU6d58IIH5FgrBb2G+kLD6DwuiogahhXNuJJeeieqSOK7lrq9mT8oXN0mQRAInzcF23Dliah0SUfydhc82ewhEVQkFM5vjnNvlP3wTsMn//XbgJjP7V7xJGU8H7qm5lCIiItJqRpxzX2h2IdLg1o0DHBnO860Nj3PtZcvHLBUW1LWwk7VXv7hsK1alm8dKAU2pCffi3oiGfceqJbNoz9poa3ZhBv1qbnRL3WCfNTNLR1uubLBW/LlKS09WEgygCpNpJjEpWq2tt+WGTAS3VWtwW3xMgdGVHE6dMYV3XLCAM0+ZNmY/Sq2mEGf5v3qNj4862WOtLeH16h7fiEn50tjtPg2JlFKSTHREPbbB72zLGJiNWdmmWfUSZdnItcBqYLaZDQAfB1ab2Tl4f1O2A38K4JzbbGbfBh4ERoBrtMKDiIjIhPR9M/sg8F3gSOFJ59ze5hWp8bbszY0GViN5x5rbNnHmKdNi36zHmTG9XECTRMATto1CQqQwh0Kp5TCjKnWDvX/bQMVuw6WWNaxV8JgUB89JbLe3f5C/++4DkeosziogcY51ubkp1l79Yr74s8e466GnRudKyGbHDuEplSgopCCSXO6wUq+McomCStdPkpNa1qN7fJJDCEqJs/+NbCWPE2w3srU+yURH1GM75jtzDnDHzNFTqIdJ+3Ksrnrv4omyysOVJZ7+cpn3fxL4ZC2FEhERkZb3Hv/fvw4854AlTShL05w1M0s2MzI6Z0Heuaq7nke9ea0U0FR6vZZWyqRaL0vdYHdvi5Zwact43fFrSWiUK1cS27xp/Q7u2LSLZXNO4Cu/2jbaA+A7vQOsfV/pQK54H6MkTqKUt1Ig2bWwk3Pmz+Cuh7wRzwa8uSu8boOrH+Txeigksfxf1F4ZpUS5ftLUEl7v66uUqPvf6FbyKOdJM4ZFJN1jJMqx7ZzaMZrUc0B71sjnjy7VOubYGJy7YrAh53AtqzyIiIiIlOScW9zsMqTB0s4sV71sMdf/YivOeV2/q+16Xs/uzsEx79eu2zwaLLxl5fzIreZJtwzGCZ6CwVChO34aukWXctP6Hfzddx8A4BeP7h7zWrlArngf33Z+dfNTlNtu2PcXn3/l5vcoHiZSbg6Lao9x1IA/eF5Xun4aMaQgqmquqVqvw6j7H6WVvBal9qPSeRInGZL0Kh/BnhpxV+CIW57Bg0MYXjIhA7xl5di/AWOWd3U0LCmmhIKIiIgkzszagQ8Ar/Cf6gb+wzk33LRCNUHfYI4bN27HOchmjDWXLmva2OwwwVatMUsC5hw3rd/BLRsHQlv8oqwqUUuX+6jSFAxWcsemXWMeFwIEKN2SXxDcx2zGMP+5Ws+FKHUX9/yrR0t63GMcd4LOeg8piCrq8IxS81aUWhkizjGLsv+dUzvImHfWFnooFPeUqfS9Ya9XWsklTJRzox69GILDluJuu5rPFObACSb2woYRZY2G/R1UQkFERETq4QtAO/D//Mfv9p+7qmklaoLCygQOcM6xaeezVa8jX6/uzsHWPZwjkzGc81odK7U8lmsZjHrDHDWAWvfYENMWl+7C24hgMKlx45csnzOmZ8KfvmIJzx0ZGTOHQlgrbXCpxLX3lE/2RBW17urZ3b6gXCAa9xhXmqCzmpbwWkQN7iu1tpe6Xkp9BogdsEYZDnXtus2M5B0Zgz9+6WJ+f9kpkZIb5cpfeD3KSi5h5a50bgTraGg4z+fufIQPX3xGIsc7uO0jw3lu2TiQyDCcYpX2M/j6pH39DUuKKaEgIiIi9XCec+73Ao9/amb3N600TRJcmSCbMW7uHRj9vXjFh3IBRz0nHCtu3Vtz6TI27XyWm3sHKo7RL9cyGPWGOWoAdWQ4z7rtPWUnpKxnMFhp3HjUY1Q45nds2sUly+ccs+pHqZb1TTufHU04zJ0xJfGlEhuRLKgkSmIpTjnLnZvVtlZXOxShkAQqNddA8TYrtbaXul5KfaYec0L0bN3DkeE84E22ecMvt/H7y04Zk6ip9L3lXi/+e7lz3yF6+6PNA1Dp3Cie2+NXfbvZsH1vIj0VVi2ZRVvGWw3FATf3DlRc4abaXlWV9rPwenf3QNzdqJoSCiIiIlIPOTM7zTn3GICZLQEm3MpPSzuzoy1GO/cd4qb1O46u+PC9+1l18KcsecW76B3YHxrgxGnpr6abcVir1xUr5h3T8hhn5vmoN8yV3ldYdjPJcdpxVRo3DvFag99xwYLQ5UOLW1L/8XsPUFi18Tu9A3ziDctaZnhHHEkHwOXOzWrnY6i2W3vh/IVjZ+Qvtc1yrdClrpewzyR9nqxaMotsxspOMlvpei73euHvZdK9cODo+fC5Ox/hV327j+nNUaq+oyaQuhZ28paV80f/vudylc+proWdrLl02WhisdkJvVoooSAiIiL18NfA3Wa2FW+Y+ELgj5pbpOYIjrP95j07RoPDt2V+ypKffQWmOXr2vzw0wIkS/FQKdm5av4M1t20i71zoTP7lWoPLbT+sxSxOV/qw9/X2D/Kdex8fDcaCyxUmIWrAUDx/QfG48SSD4eB3mRm5wrTuHO2636ix/nFa5Js9D0acIQyrlsyiLesfz4jnVDXHuPCZwhE0xs6TEbbNcq3QwaEvVvR8lERhLboWdnLtZcvH/C2Ju3pHlNd7tu5JvBdOYdsfvvgMNmzfO3qedU7tOOZvG1C2V0lQ4bxbdur0MfMbhM3jEBw2VZgAd8P2vWWXE047JRREREQkcc65u8zsdOBM/6mHnXNHmlmmZuta2MmrXnAyP37QW37vvMwWAHY88yyrzg4PpqIEWpXmMlhz26bRVsWhKm7Qqw2Yo3ZRD3tfz9Y9o+U2yi9XGFecFufiIKhQtnq0Bge/q3NqB5+4fdPospLZrNcNHBjtZl6v4TBx6ieJCe/iBMBRJyUsy7mx/1ZQTcKjOBFVvGpKLUmUWzcOMDSSL9uCX8tQlrDz6h0XLODMU6aVPU5Ru+WHCauXJM714vOs+G/bLRsHxvSKguiJ3HITfxa/94oV81KzTGmtlFAQERGRxJjZRc65n5rZ5UUvLTUznHO3NqVgKfGnrzyNnz3yFNfbp3hl9ncAfPs3W7jwbEKDqSiBVrnA5NaNA6NBOUDGLFbg0ts/yH2P78PMyOAa2s2+eNbycssVxhV3KEWpFuDg70m2Bge/68xTpnHLxgF27z9C9yPPjOkGDvEn3osqThKplh4axUFiNUMP4n5/IVHl8JbgTGJCvGo+U+15U4/5EYIqBcq1JCqiKFUvSa7SUFz+4N9Og7K9SoIqTfxZ7r2O5IekwNHradK+HKsT2WJlSiiIiIhIkl4J/BR4Q4nXHFA2oWBmXwEuBZ52zi33n/sWR3s6zAD2OefO8V/7W+BP8OZn+HPn3I/8518L/DuQBW5wzn2qpr1KQOFG75Yzu1n+2O9Gn5+Rf5aerXu45sKlZVvIK908X75i3piVAgrf+Z17Hx99T2EyyDjB3pXX/+ZoC3kNS19WIxhYJDlrefFQCszonNoRexuFbueFOq9HvRS2e93dfdz50FPHjP1OIrAs1fIbp/U86nuT6FkQdVLCcuo1IV41n6lmm7UOD6mkeB6PcsOl6qW4XqpZSSHq9xT3Prpl44A/5AiWnTqdt523IHYit9J7r1gx75h5amo1ZvJYg3NXRJvQslZKKIiIiEhinHMf93+91jm3LfiamS2OsIkbgc8DXwts822BbXwWeNb//Wzg7cAy4FTgTjM7w3/rdcDvAwPABjO73Tn3YDX7lITgjd7vd/yQfdOW8uJnr+UnmT9ndmY/84tmoI9zk1kI+odzjvascXmgFb94yMDbzpsfOhlgKT1b9/iTEHryecfgwaGqy1qNQmCR5KzlwXoBr5X62nWbI49jLk60fKd3gLXvq2+gFRa81BpYhgX1cVrPo7w3iZ4FxfVQWAkAwnv4VFveNAsrf1LXY/E8HnnnEu8NEbes1aykEOf7g9sJTgz5wBPP8vBTpf82JHGNJHnuBXtdjTgaNoxCCQURERGph1uAFUXP3Qx0lfuQc+7nZrao1GtmZsBbgYv8py4DvunPzbDNzPqA8/3X+pxzW/3PfdN/b9MSCj1b98DIYf6j7f9yhg3wsynv5E2nLSD36BxeNe0I0wIBQdwW21s2DowGtkM5N6blrlSrWByrlsyiPWuj2y8ezxxnjH1hToCwMcb1FgwgCvVSaZx0WNBTnGhpxBjosICk1sC4XFAfp/W80nuT6FlQ+J5SKwGsuXRZtB1OUG//IOseG2La4sa0BAcV13fSQwKC83hcu25zor0hqilr18JOVp950ugcNFFWUqj2+7sWRp8YMslrpFC2cnO1lPtcsNdVJkPDhqYpoSAiIiKJMbOz8HoMTC+aR+EEYHKNm3858JRz7lH/8VygJ/D6gP8cwONFz18QUt6rgasBTj75ZLq7u2ss4lgHDhygu7ubSYPD/E3bt/j9bC8At+yayfcf38FZ7ady5eFufvbTu3CZLOseGxoNcoeG86y9cwP7TyvfFX/XE0eKHu+ku3vP6OOPrOhgy94cZ83Msn/b/XRvK95CeR9dOYlfPTEMGC+d2za6jahl7RvM8ekNhxnylq/3xiRn4KPnTWZpZzZyOQp1WY1CGYbzR7/7Iys6+NUTw/xiZ458HrIGk/b1j/aE6N4xzH89NETeHVveSftyZM1rBQTIZsZ+Nvi9hbovt69R31f47rV3buW3gfcuM9i/bSDysQ3W5aR9Odr8fSmug3LilBngmceHMbzjX/ie/dsGRo8DGL/duJH926KdE0N7h0bHoh/xl9csdazCyl58PoS9P2w/j27Dcftjv459Piet0vVYaj8qHcNlBhyq/W9I1LKWu8b7BnP89KHDo48txrka9fuDyl0Xcc/9qILnZda/WHIRztHCPhWSnAZccJJL5FhFoYSCiIiIJOlMvDkQZjB2HoX9wPtq3PaVwNoatzGGc+564HqAlStXutWrVye5ebq7u1l9Zier/+MV3mwOwB1zPsgPtl+AA36TfwF/mP8RrzzlAJx9GdMWD7Jue89oa+CVF59XsWVq2uJBfvWlo5/50BvOH/OZWvdoNXBVyPdGKevmu/sYcQ+PPnbASB6OzFjI6tWlJzArpbu7m2qPT6EMDsg577uvuXApV1G6F0Jv/yBf//FvRpf4HHFjy7sab3xy8RwKQb39g3zmrkIraK7ssnNR3hf3veV0d3czbfHvefu9YhZrV0RvCa2mHL39g3zzrh7yzpuH4xOXLR8dejOtf5DPbOxhaCTHb57KR96n4PlX6JZfOL5bOYkjbkro/pQ6H0qdi4X9PDKcJ5sZ4dpAuY9uw8puo1HKXY+l9uPMU6ZFPoarI3x/nCEMYWUtd41vvruPPN7fEQPeet4CrnrTCyOULPr3B63Gu8ZL/W1I4hospfi8xFHxHA3bpwsXdlT99zIuJRREREQkMc6524DbzOzFzrnfJLVdM2sDLmfskIkngPmBx/P85yjzfON9/8NHf1/8Sk5a/VHab+iBkTw/y5zPSMd02h67G86+LPbY7sJN/CfeEL5cWb1ELWup4QV5iD0JYi3Kda0v1Q25Z+sechVWxqimi3/h+WB9NWo1haC+wVwgKPK6fIfNTl9rOXr7B/ncnY+Mzpzv3Nh5OGpZkrRUt/xsxri5d4CRXHh3+uD5cHzbCC+ZV/pc7Nm65+iY9LxjzW2bRsfSF7YxNHz0nIoSVMcJvOO8t9z1WGo/3nbe/MRWiog7hKGaOSxqHb5VzfeH/W2o1wobxfODYEYuF22oSfE+7d92fyJlikIJBREREamH35rZNXjDH0aHOjjn/rjK7V0MbHHOBfu33g7cZGb/ijcp4+nAPXgNWKf7k0A+gTdx4zuq/N7atfm7/5I/g4vW0NXWMebGr+1Hp8Hg9tG3Rx2Tm+SY6XLfUct684X3fOOqVXzuzkf45aO7cUDGGBNURpbPAw4yMboZO0fX3KmRAojCyg279x+hvS3DyEieTMyVMQqKA6DOqR0lj1dSqynECT637M3VFBTFWdXhnTf0jAazGeOY4Ltzakf5beVz3r8ljnnXghl0zc5B35284oJ+9uzaRvbATn64Zw4HMpN4Eds48Xv/BjwNLgftU6F9Cl3ZSfSe9Dwjz+/lhCO7sG/kYfIMmNsF7VNgygwYep4rnx/iqbaT2Z0/noNMYgensOHRJ+ia00HXqVP4xlWrWHvnBq68+Dyg8hKeN63fEXnFhGrnGSj1nlVLZpHN2OhEpIXeHMWTW/b2HzsXRJTzqpogO87cA4X3JzGRZnB/4iTRCuq5wkbxPkK8nkPBOm3EUIcCJRRERESkHv4L2AK8BrgWeCfwUKUPmdlavN6ms81sAPi4c+7LeEmBMcMdnHObzezbeJMtjgDXOOdy/nY+BPwIb6DBV5xzmxPar+gODbK6+zLv98WvgIv/52hQNOZmeuZi2NEDznkDgyOqdRm1SoFC0pO8ffjiM9iwfW/0G/F8DvY/Cc88BJu/y2nPPAu9H4AT5sAbvwjTToEnH4Adv4GVfwyPr4ed98GkafCit8Fxs2HTLdDz/+DpLXQtexNds5bAzinwdIcXQD5+D0w6Huacw8ZJ53HlV347OgFlW9a48vz5vGX5CZxzeqDDy9Dz3vHqOB6mdHrfM3Wml+zYuxVmnQZmdC3sZM2ly7hj0y4uWT6HwYNDJYOu4pb2Qk+GsARO2Oz+cY7VWTOzdLTlqg6KogZ3hXPUARngpUtn8+GLvYVYguV974sX8fDO3bx93j66dt4EjzwNB57y6nn/kzByyEsGHDcbZiz0jsHuRyF3BHJeYmoeMK99Kocnn8hH234GwAE3hcPudPraltA5fTqzOkZg+BDkjnDczOlwyukw+wzItsNTm2FwG+w/AocGoeM4Zh45wLVtT47dqV/6P0DXzNM488gIx6+byrb8iVztTmKTLaIt79i14SnILYLjToI9j7LtyT387qcP8XYbpsNGaHM5nv9FLxxYBrlhmD4POhd5+7b/SZ6692EW5HLMYD/bRuay/rGn6DoJOPwsHNwDx53onYNDB+D53d71Mvw8PDvgJTE7jofhgzB1Jl2Tp/OvF5/A3965m5yDtrbs6LKFt/Q+zq833kffhof4zMYsH/+DF3LWovkwdIBNu3P82c2PsXt4Cv+WafeGfKw4CQ48CZh3rU05Nin2kvlTvP048BTsexxwMGWmt4/Zdq+s+Ry0TfKO66G9Xnmd3yvIORg5Au1jp92plIRoxN+0anuRVZMUKDzu7R/kurv7UrsiibnCgWuilStXunvvvbfZxRCpu0Uf+0GziyANtP1Tr292EWQcMLNe59zKZpcjLjP7rXPuXDP7nXPuRWbWDvzCObeq2WULk/j9yPZfwo3+34GXfwRe9Y+l33fPl+CHH4ErvwVnvjby5ouXLuxoy0ReujDKjfV1d/fx2R8/7I17N/jLV58Z2qIX6ab58HP0PpU79n25YS/I2P8k7N0GT22C/l97iYL9u7z3TOnEHXoWm3067H649PaLFYKt7CQvkNn7WNm356yN3+aWcH/+NJ52M1hiu3j1tG10HuqHmad5SaHndsLO38LzT4/98PGnwOF9MHIYzroUXvaXbN49zCdv+Q0duUMczB7P2179Sj75460cP/IsT7edxE3vPZcVS04ZTSL1bt/LP335WyzK7eCwTeHKCxbyyrPneoOos+1e4DV3RcmW+orHauh5r56nzADg53f9mFmzpvPbJ57nBWe/iK7Fs70yBI/j/OnePk2e4ZXRDHIjsH+nF7QOHfQC8F33ee95/hlv/w884wWIR/bz7PFL2PDEITL5HIM2nYsWT6Zz+Cn6hmby26dGmMlzLLYnOcX2MtWCk4uaFzDPXeElELLtcGS/FyDvfgQmnQAnne31JjhuNix4sZcY6DgOMlnue3gr92/dSceMOXz8B4+MBrqxl/bM52BfPw/27+Kh/ifpOu5pFk05DJbx6vSpzTyzezcnzp7NoV0PMeXZ8udYWjjLYsedCC7HyKHnaMsfqfiZYZflMB1Ms0NjX8i0QdtkcnmHjRwiN2kG7Uf2VlWuXGYS2bYOyA9759Lk6XD8yV59u/zRn3zOSzpMme5dH1Nnstcdzx1bh9mbn0omk+Wy889g3ikne0mPSdOg4zi+/btB/qPnSZ7PT2JG5hAfWHk8ly1t9/bBfw9DB/3zqA2OPAcH98K+/qOPMS+pUzgfR454ZcpkwbLe35x9/eAch3PGzue8iV2n2hFmHz+JDstDW4eX2Mq0QX7ES4p1LoJpc7ztdxznXW9TOnlq/xA/2vwkuTxkMsZrl5/CySdM8b5z+JD/czDw+/MwfIjtx53Doj/6UlXHIUzY/Yh6KIiIiEg9DPv/7jOz5cCTwElNLE/j5YaP/j779PD3/d6VXkLh6QdjJRS6FnbylpXzuWn9Dm/irhjLqEXpohy3W/sxyYmRIXjsLq+XwEPrYOQQXXN+j64FL4Yjr4Jnz4bb/9x7T+ci2LfDu0kuWPASr27yI3Dh3/HzX/ySV154MWz8Kjz0fS+YPOVFcOJZ0PMFeNFbvWD+8fWw/Rfw3BNw2kVwzru87Y4c8oKDTBb29HkJh5lLvCBg528ZXH8T7Y/8liuzP2WqHeGQ62Bo5ktg3iXQ/2tGHriV/dkZtJ2wkGkX/QMcf5IXWD+zxUuEnDAHDu2DR38MW9axDLgpy+hknNwFV/iPc21TyP7XIZh2Kpx4JmTaeNHWX/K97KGj79/o/wRNmn40iO443qu3U5Zz6XCGye2/5gS3n+02l9dkD8Ctn/aSMm1TvCE1+RFY+BKwDC/t/w3Z/BDLAO47Dk44ledtCu3PPM8V7MV1G8Nth2jPHQLMq7PJM7yWe5cbW6b2qV6Ad/zJ0DEVMu1ecDb7dKbv28FLZh4hd3CQSZkn6DhwAkyfz7z9fczM7GWnm8U25nBXbgX73HE8wYmc+7LX8Z5XX+AFWzF67ASdc+YSzjlzCX//3QcYGvHOqaGRPLdu9EZMlZtMc4xMFmYu4eyZSzg0c5AfbN3Dqvljk2a3ffcujsxYyKoXzyIz9BxbN2/g7HmzecGJHd75MHIIOhexaY/x59/ZxIHhLMOZDv7udWfzlhd1eudktt1ryd+33UuWHH8StE1h26Ob2DyY5ZzJu5h3fMar1/bJMHW2l+w5csBr5T9hLmS9wJ4T5npB7tDz3rmyfyccHPSO/3NPQH4Eyw17vQcybew+0s4Xfpdj+8iJZDPwytNn8Yr57Sw+9SR+9dDj3HHvFqbzPFPtMMdxhLOWLmHVi5Z5x+bQPq+3RG6IrHPQPpnMgadgxiKvPMedBDPmewmBg3u866WQQLSsF/yODMHkE+Dws+x87CHmz5vnn2/TvSTj8894yYNM1tuOZdhzcIRnDgwxJ/s804+fCgf3Yru3c4ntYUb2eTLmoERu+K3AW4PTZfzO/4mqkNiYPN37aZ/q1b9leP7IMIcOH2HKlKkcN3clZLLseHIfj+x7FiPPYTeZs48/gRfMnentc/sU75hk/YTGUw96CdR83ktKuDwc3sf0kTxvtNzo34VJDxu0ZQHzvrtj6uhQnudyHewdmsYJJ8zjSMeJMXasNkooiIiISD1cb2adwD/izXVwPLCmuUVqsPzI0d87F4W/b9LxMHWW16oVIqwHwOUr5nHLxoHYXdejJAvidmsfk5yYOeT1ztjjr/A5/wIvMBh6HtZ/0fsByHbAqg/CMw/Dma+HE8/wWvvndnlBcyCgdJl273HXe72foLP/4OjvL7jU+ynWPtkbogDHHo/TL6a/o4srt3iJkRMyh/mH1y/jrS89Cyga+34gwzdml2np3v0o/O5bPP30U3zuwak8MnIKc7LP8dEXH8/8qSPQPpnss0/wxPBxPP/EZuY+u5PjMiPsX3ARn+2bw4bc6ZzA85gZ7+46kT84x+8mfuBp2PYzL1g8st8Lzn73TfgdLAT+JBMow09vItd2HNkzX+O1fi59lRfAPPpjyHawb8YLmfWK93nJgV33w/4n2f3kbo64I/zSvZC8Mw7mp/CqVSuZP+WI1/J5aB+7XCdbDp3A3EVncsb8OV5L6klne98ZEvxPLfHcZGBzYA6Fz67bzLA/+dy7z17m7W8CivthP73/yJhePd/pHYjUayEsadbbP8inNxxmxD08+vwVb3xzyW0snwv/Mn3xsdfTCad61/dAllVLusaUZfGCC1hc9d77Tj677MunAH9wnjd3yM29A3RvydPR5+3LfTP28PXc0eE+bRnjWxe+GAJDbXq27mHV6cl0xX/MuplfYWWCsGOx1X9+eCTn9UZ5z4s496SMd60MHfD+9gw9z9YnnqJ/19MsPPUUliw+zUve5Ee85MyR57ykzMhhL0k26XhvqMaMBV6A3+6fzZlM+TK9wSvT/v5BPnLD0ZUXvvG6VaN1F9Xm0f3yt/GH4avFjJbhmQwfWdHBmbG+qXpKKIiIiEjinHM3+L/+DFjSzLI0TTCh0DY5/H3g3bDu21HypVI30HB0sq5qJiqrZZZzr1Bf9XoFdBzH26Yu4ZyOh5jpnqXdcszvzcDPHvfe95p/hrPfCNPnHv3sQC88eb8338Dpr/aGEqRAz9Y9jOS81uwD+cl8f8uznDZvEIA1t20andBuqNKkc7NPh4v+gZOAKwJB8+0Hh1i1eNZoIOod05VjgqJl63fwrULiIpth7sqiAGT55WO/K5/zhmH4Lbq9u4b4/274BvPzT/Dz3EquP/+iseV8lZfTe6C7m9XnrPaeO/ddAPxq/Q7+4XsPUFjcImPw/NSjQyfGnIcPZVhz6SneyiKH91UVTAbPrTNPmVbzZHulXLFiHjff+zjDOUd71jhp2iSGc0fTDFEnEAzr0dOzdQ/DeS9xEWVbpa6nRkyuWklhX0ZyY/dx1ZJZTG73VrIonpy0WeUOOxbFf9PODZSlt3+Qnu3++XVRZ6L/IQVXMalUpmDdJbF6R1BxvWzZmyv5vnpQQkFEREQSZ2YzgD8EFhG433DO/XmTitR4wSEPlVpcZyzwJoUrIbjk25A/+eKtGweqXvKvoNwEZ2E3vL39gzz3yy+xuu9TmN/1ffZxOzjhuCxtB58iN2UW7fm81+PivPfBBR84pjWPeV3eT8oElwHMA7/q282G7Xu5YsW8iktIhula2MlPNj/JZ3/8MM7BpPbMaHBQKgB5xwUL4gXXmazXpdzXs3UXG0cWc69bTDZP5CEwvf2DXLtuM4HdJGM2ZmnPYJmHhvNjVitYc2lty5bGnfG/nOJzd+3VLx4za/537n18tIdC1F49YT16Vi2ZRXsGcq70tuq1QkI9lNrHcsFsPcrdN5hjc4XJB+MuAVvPxEe5VUzCypTk6h1BxfVy1swYK+HUSAkFERERqYcfAj3AA0C+wnvHp3wwodBRPriYsQAe/m/ouxN23sd9p7yZXz0xwqols+ic2jHadTsP7N5/JNH148NWDMiOHGRb+0bmL9vPSbvu5mBmGsfv2U+X7eAX7veY9o4bOWf68/QemXfsfsVcsSINCsHT5+58hF/17R6tX4eXCCjVSlvJTet38MWfbx19fGT4aMtvnKAoqmqXtCsEh+CtuWrmLS147brNnHnKNLoWjp3J38zIO1cyuVBtwBZ3NvywbZQK1oLbW3v1i6PPoeALC6y7Fnby0fMme3MolEi+RQkca12GsFzyL+7qAmH7WI9yl9qP4uEjYS36cXpl1TNhE7aKSbnt16s8xfWyf9v9NW8zKiUUREREpB4mO+f+stmFaKqRozOn3/bA0/zVTx4ND7pmLPSWwPv6FQDsdT/ks0N/RUdbhitWzCNjkHdeC9jsaZOqu5Hf9zg8dDuccCr8+B85xCSe3TedJfks+e7nGJm0i7YXvJ7MngxXuSHOa9vCKzO/g4fgcNsJPNvWwek8zpb8fP5q+P28Z2ee3OR5pYOmFCcTygVaXQuPXd5y+anTAWIFoAV3bNo15rEZVQVFUVW73bBkQXA50uC2O6d2cO26zce8v9oAqdZW5MIx3bnvUMVgrdqETdjnlnZmWb362B5CUQPHWs6FcnM7VFOfceom6XM4bPhIuSRRcDnFwjaKf0868RFUvO1KyYRSn0myPMF6WffYENMWDzakt4sSCiIiIlIP/2Vm7wPWAaORtXOuuvXEWtHw0eXV/uUnWxnJzwTGBmmjZp8x5qMXWS8ZN8LwSBsOxtyAFtaPj3Ujv+9x+MJL4cizo09NAS4yjq4qMATcfxPnAuf6d4hfyb2O77uXcf+RhbRlMhxvh3g2N3n0Rjgt3bVLKdf7olygVSpwLrz/8hXzYpXhkuVz+MWju0cfX/3yJRVbfmtVzXaL9/kTt29iKOdwwM29A1zhJ1JKzXsQTC5UGyDVch4Fj2lbxmjLZsjlag/Wau0xESdwLGy/Z+ueMY8rKTe3QyOuyyTP4bDhI2H7UnzcMWMkN/b3wjUeN/ER9dhHTarU0sMibvmCwzDWbe9pyNwWSiiIiIhIPQwB/wL8PUcnW3dMoAkaH396L4XR7UfyR8ezFgdpgLecH0DnYnac8R4WrP8EfZP/kNflPssVK17CFSvmcd+WPq7c8iGmDn6Y3s7oy0vy1IPw9cu9ZMKLP+RN5Lf8Cn63/zj+z9rbWOie4Cf2Yr7+0j0sePEV4PI89NhWep/M8djQidx/zw7yDnJ5x2vPP4O5M6aMuZGtV2tbLcISB3FajbsWdnLd3X01BWbvuGAB4PVUuGT5nNHHaRQMDjfvfLbicqRhyYW4QTHU1mobPKa5vONt588/5hyNK6keE+998SI273qOS5bPKfv5ar+v3NwOabwuy+la2Mk7zurgseETxtRX2L6MuZZzDvASYGN+96/Zay5cWlWCKsqxqJRUiTIMJ45K5QsOw2hUklcJBREREamHvwKWOud2V3znONTbP0hu/c3M9+cjdG3t2MjRzMoxQVq2Hf7yIZh0AgssA+s/AcD3pn2GjlPfCe1T6Dr0NPxmC3zv/TzrunC5+bzrp2/l61e9ZOwN496t8F+Xwyv+Gh74Nmzt9p5f/Ep4zSdH3/Yi4ANXzaFn6x7+bcksFgS28YKu+bzA34/gspRXFHX5r1fX/VqFJQ7iBlpJBGbvuGBB4omEJOYbKCfucqSFMkQNxIrLX8t5VHyMis/ROMKGTtyycSBy2UpN1Ldh+97RuShKqbZHQbl5D9J4XZbT2z/I17cMkcvvZv22o/UVti/B4571eyXkcmN/r+aaTaJ3R/D8Trq3SKXtBSeXbVQySQkFERERqYc+4GCzC9EsP9jwEGsyW0Yfr1p6CidMn8HNvQPhN7onnHr09z/9BXznvXTsfQwevsNbLvDph0Zfvsh6uaitl77hufRsPXPsDeqv/g8MboPbPghAru04sm//r5LLM1ZqKYsSmJQaz9zsACYsERA30EpjYFbvWesL+xp3v6MGTjet31FyEscoLb2lypPUMQobOpHNGDf3DozpQl/uO4ItxECkQDIYBFrR6hqV6iGs3uo1pCauYHmB0ON068YB/HlBGRrJc2tgWFipfSk+7sDoZJvLTp1e9aojSUySGbw+11y6LNHeIpXKV6iXtXdu4MqLz9McCiIiItKyngfuM7O7GTuHwoRYNvJIdhp/cOSfuH3SPwIwc9pxXPumF4bOfXBMsDTnRfChDfAvp0HfXbDsTfDoT8Z8x4P5hXyu/ToOPPgbeOVd8PAP4fv/Aw7uZs9pl3P9o8cxkofbc6v5YvsKuiotXRlQKmip9P5mrEkfplyQGTfQalZgFhZA12t8fKljGGc50iiBWG//IGtu28SIvz7lUMTyVzq/4hyjKPWayzsuesFJHB7OMaU9y50PPRW5voPJgTyllxIs1rWwkzWXLhtNtARX14hTD2kTNs9BqbK7os8WPy4leNx7+wePWU63mrqpNUFVfH0OHhxKLOEVNdnXtbCT/ad1NOzcUEJBRERE6uF7/s+EdPmKeXyn9/TRx5d1LQJirpOeycK88+DB78GJZ8LjPfDaT8Ev/40nzng3vW0vY+mmD3D8MxvhyxfDE70wpRNe+Ba+M+Mv+NKD28k7yBpVT3IX9cY8jZMz1jMREDbhY7mb/DjDFModgzgtqHG+s9ZjGCUQ69m6h1z+aKiYMYvUYpvU+RW1XrMZ42ePPDM6yV+ciR6D9dA5tSNyS/ngwaGyq2X09g/yuTsfSd11Vk7oPAclyn7Finl8+54d3qSMWeOKoglQK53LSf4NquVvR6nrs9a/ReWSffUe/hSFEgoiIiKSKDPLAu91zl3Y7LI0S9fCTta+bxX8p/940czQ95abxfxw7gW8dOjH8JN/hCWr4fyrYdUHmAu8G+DVv/N6MTzRC+3HwVu/BotfwXn9g3R076h5kruoN+a1dBOOekPcN5hjcwqGVJS6uYfy8wcUt9S+ZeX8sktQljsGUVtQw8rZs3UPk/blWF30/iTmi6gUOK1aMotJ7V7rfSZjXHvZ8kjHMolu6KXmRiiu1zWXLuOOTbs4Mpxjw/ZBb1LKookegYpDe6oJIMvtY6l5GVphssWxQzkgkzHyeRfaVf9j50/myIyFJXtwVUpyJnWO1Pr3pR7DpKKsdNHMHitKKIiIiEiinHM5M8ub2XTn3LPNLk+zRL2x65zaQcYMOHqjXbhRHB45j5e1/QNrXjWH0172Fq/XQlBbB7z3B94SlQsu8G6K/WAnqUnuKt2YF27E11y6LPa45ag3xL39g3x6w2FG3MNN7+pd6uYeKJuECX5mKOe4af0Obtk4MGY/ggFNlHHScXsc3LJxYLRLeJvBuSvGrlHfiPkiqv2OWsoWdVnJ3v5Brl23eTRoL8hmj070GJakSUK5fQzOy5ABXrp0Nh+++IxI9dDMFuzioRwGvP38BaHJtKWdWVavPnaYTZQkZ1LnSNwJRcP2O8m6jrTSRRN7rCihICIiIvVwAHjAzH6CN58CMHHmUIiqEMTk8o5sxlhz6TK6FgaXK8zwi+Gz+cRjs/nw4oN0LZx07EbmvGh0W7WMgS+Ic2NeawtZ1Bvinq17GM4Teym0qMFUnKAr7Oa+XAKg8JlCsFq8H6XqsdbgvricxtGkx4grPQymEfNFFI97jzNBZjVli7qsZPFkigAGvLnraPBb6nxdZrGLFCpsH4uPZZxkQrNbsINDOXJ5x6kzpiTaeyMoiXOk3N+XJOqz2r9JlVa6aGaPFSUUREREpB5u9X+kjGAQ45xj8OAQcOzEbr/q282G7XvL3sA2Ywxxrd9Z6oa41A33qiWzaM/gja+OeOMcp/dD3HXnS93cl0sAFFpqv7VhBw/ueu6Ybt+l6vGaC5cm2vUaGF0KMms0vbt8vYPdwnnUObUj0rKSpSZT7PDfX/ye4Pm6f9tAYmUOU23rexpasCsN5YiyT/XuPRM1MK+1Piud88Fz9tp1m495X7leGZ1TO0Z7SzX6GCuhICIiIolzzn3VzDqAM/ynHnbODTezTGkUdiNbuFH83J2P8Ku+3bGWnitsq3NqR92XcUxi7oTgUAkoPRdB18JOPnpe6fHVYeL0fogbJITd3JfrJl0IENoydky373q1NBaXqRB8TNrX37QhIwX1DHZLLd1XaThO4Zq7ZeMAu/cf4cRpk47pml8qsO3eVl35qhn2Ua+W/aQV719hbopLls8Z0zslzvCRevaeiZqwqLU+y53zwfrImJWdoLO47A8/uf+YpVgbSQkFERERSZyZrQa+CmzH6zk838ze45z7eROLlTrlbmS7Fnby4YvPYMP2vaE3sGE37svmnFCyhauR5S+n1CSFUP6GO2x8dZioN//1CLqKj0txt/vibt+NmL+g8D1dCzvp7q5/q3ol9Qx2Sy3dF3X4T3DpwcuLVhqA2gPbRg5DaNR5FVQqmVP4W7Rh+97R5TDHzCtSZvhIPYYtlRLluNZan+XO+WB94ByZjGGUnsAyKGwp1iSH4lSihIKIiIjUw2eBVzvnHgYwszOAtUBXU0uVQuVuZMvdwJa7cf/NY3sit3DVs/xhwiYpXHPpssSCzKg3/0kHXaUCxijBc5R6TMMScUmpZ7BbbbIibq+J3v5B1j02xLTFg5HL3+hhCFGvz2B3+7iTqwYV798dm3aV3N/OqR0UVhDNO29yWg4dW6Zahy0lfc3UklAqd84Xn7NRJ7nt2Vp6KdZGDMUpUEJBRERE6qG9kEwAcM49YmbtzSxQqwq7gS134x6nhasZwiYpHDw4lGiQGfXmP8nu1GHzIZTq9h1HNS3b1QRTUT6T5BJ7aeo5EycREVzKcd32nsg9DdIykV5QqWUpq+09Ubx/lyyfU7KX1eDBIQxGV64YPDjEqUWt6rUOW0rDpJTFws754qEh77hgQdntBBNApZZirWYoTrWUUJBxYdHHftDsIoiIyFj3mtkNwNf9x+8E7m1iecadSjfu1S7j2IgW8OB49Zt7B8jl8mQzxs59XhNlNatTpEXYRJOlun3HUU3reTUJiEqfSWOQVko1yYo4iYjghKpRj0dhu40ehlBJ8QoX1fSeKLd/Z54yreREq5Pay09wWeuwpVs3DowmSYaG83zuzkfGrJBR7u9drX8Lwz5f7vmofyOqmSOknpRQEBERkXr4AHANUFgm8hfA/2tecZrjdy/8OC86c2Fdtl0q+Cl14x5V1BnIC9tOaszyFSvmjSYW1t7jDX2IE6SmbRhAqeNydBnQ6ru5x23ZrqZrfZTPpGHlgHqKmogIrgoRtTdD8ZKuvf2DdZ84tbgc5brbB1e4iNN7otKStaXqtNR1UtyqXsuwpd7+Qb5z7+NHkyTALx89uloOlJ4ANmx/4hyfsM+X226c66qWOULqQQkFERERSYyZ3eWcexVwrXPub4B/bXaZmmnvrBWwfHVN2ygXMBffqNfShTzqDOTF8zVUGtscJRjo2bqHkVz8IDWJYQD1SEgUH4ckurkHg6YoS8RF/c7g/kf5TBq77Neq2lUXvnHVKtbeuYErLz4vUm+G4PkN4QFtPZS7VorPrbit3WH7F+Xar/Qd1Q5b6tm6Z3SiwgLH0UkL8X8v9Ten1qRZ2OfLbTfOddU5tYOMeQNG0nANKqEgIiIiSZpjZi8B/sDMvom3wsMo59zG5hSrNTWye3lYV/2erXvYue9QpInWopa9OIBr5CR6URIj9ZjILYlu7oXPRTkngsNKwiZ8L3WMKpUzzr6krfdIKbV0H+9a2Mn+0zoqvq/U+Z1ET4849Vvp+yoF7qW+KziOv7B/2Yxx/+P7+Pe7HmUkV/+/W5V6XRSGPBQUJi0EQv/m1Jo0C/t8ue1Gva4KQyNyeUc2Y6y5dFnTry0lFERERCRJa4B/BOZxbO8EB1zU8BK1sEZ2Ly++oQXGLO3Yls2Qy5WfaC1K2cMSDY2YRO9zdz5SMTES3O9y8whU06KdxLGLe04UlkEsNZQkbALJJFqN650MSypZMWbFkeE8a27bRN65RMscdn7XErTGrd9aguRS3wUck4jZtPNZbu4d4CcPPjUaxJfq7VTPFVWKe10UhlONjIydtBAou2RvqSEUUcsddrwr/Z2Lcl0F57twzjF4cChutSVOCQURERFJjHPuZuBmM/tH59w/Nbs8ra44COic2lHXMdfBG9rguP9c3vG28+czd8aUyPM1hJW9uLdDIdioJuCO06pXPIt9WGIkOJFbqYC92ZMSxgkMKyUf6jl8YUygnnAyLMljEKwDM6vbcqulhieVO3crBbA9W/eMmXDwlo0DZQPgWnrJhA1pKB7HP3fGFEZyR3sEGGPnYkjiuAX3K2qviytWzAtNHEQZjlFNucO2XWtiMY1DjpRQEBGpk1ZZfWT7p17f7CLIOOSc+yczmwssJHC/4Zz7efNK1XqKxzZHmbcgKcU3rlesmBeri3RY2Yt7O9R6Qxy3VS8DvHTp7NHZ3oOJEWDMRG7Z7LHla/akhHECw7BhLOseG2La4sHEhmKU0jm1g8IQ9rzzHiclyWNQ6jxtVLAWdu5GCWA7p3aMmXDwO/c+Ti7vyg7lidqzpHguhbAgttxz2YzxlpXzuTzwd6PW41ZqeEqU4LrWID5quRsxxKee12y1lFAQERGRxJnZp4C3Aw8COf9pByihUEGplsWuhZVXCkjjuP9SZS/V2yFMcJ9qURwQBZeOK+6VUZjIzYA3d807pnxpaCGMGiCFDWM5Mpxn3faeMYFmHFHOtcGDQ3jTxnlJnCS7Zid9DIJ1UMtKKUkpN1yoULbNO58dfb8BIzk32qsm6hwnpeZEKO7JU25ujSjPBVeyqPW4lVrdoJahCWH1UCxKuRvZc6nWBEnSlFAQERGRengTcKZz7kizC9JKyt2UlruprdfNbFI3rpV6O5RSvE8fWdHB6ojfV21371LlLJbGFsJySg1jCRvOEUXUc23VkllMaq9P4qWexyANwVpYz5LgnCb5wPuzWSNjFnmOk7BjGBzuA1ScW6NUXVUaKhD3uFVahSTs+9pK9JAote04E5wGE3PFQ8+q7X0RltBohQlNC5RQEBERkXrYCrQDSijEUO6mtFwQlfTNbNKiBoDlxkdv2ZsLfW/xjXipICFKoBi1nGkIOqtRCMiGhqsP8qOea/VOvLTqMYiiVN2N6aGUO7pugQFvXTn/mHkCyvW0CJsTITjcB47ONVJtMqjaST8L4iYkxszbkXPctH5HyclIS73/iD8PRaVrPuzvS9TeF8W9rkptq9nztMSlhIKIiIjUw0HgPjO7i0BSwTn3580rUvpVuikNC6Kq6Uoc7N6c9Wc/r2d370oBYKXx0WfNzIa+N3jDXes47YkQqK69cwNXXnxeVfsZ51wbz3VZb8EA9rq7+45ZmpFAj4Tlp04vuaJAuZ44bRljOOctPVhI4AWH+1x89smcM39G1UMJCt+T1BCHoeE8n7vzET588Rlcc+FSIHz52UIvi0o9cQr1MOQPF7m5d6Bi76mwvy9REmjFf7euWDGv5LaaPU9LXEooiIiISD3c7v9IDNW26lbzueAs8SN5xz9+7wGy2UxD1o4vVmpJx+Lx0fu33T+m7GE33GmY4yAp9ehB0rWwk/2ndUQeUx936buJolA3k/blKg7FqfY4lkqyFSZKBGqbrNX8GS7MgGOvm/e/8rSyQxeifEet50qwR00e+FXfbjZs33vMspXB4Q3BpSIrTfzatbCTt6ycz03rd+CAXK5y8F7u70ulBFrx3y1H6Ykto87ZkJZrUAkFERERSZxz7qvNLkOrqrZVN+7nVi2ZRTZjo62SeQf5GsfXV6PURHDZjLFz3yGA0dbI7m1jy17upj7NAW9xIBB36EY93LR+B2tu20TeudHvAi9gOzycJ2Nw9cuX8LHXvQBoXM+DNAVNQWPG6hucu2Kw5nH6pZSahLBwPQDHDoWIeN32bN0zurxjIYi+5sKliQ+pKpSx2mNXuJY/d+cj/Kpvd+iylcXDG/75TS8MXSqy2OUr5nHLxoHICchSf1+inqel5mgpVc5Kf8PSNiSiYkLBzL4CXAo87Zxb7j83E/gWsAjYDrzVOTdoZgb8O/A6vK6O73XObaxP0UVERCRtzOwBGDMMdwzn3IsaWBwpo2thJ9detnw0kGwr6kZdS8t+nECweEnHF86dzkNP7mftPeHjnyvdcKe1q32pFuew1uVGdXvu7R9kzW2bRhNLQ4GA7fCwN/Vf3sEXf76VBbOO4x0XLEi8DGHlSlPQFBScvHDEUfbY1HIco7RUV9MjJ+wzSQ6pSkrXwk4+fPEZJSeZLDe8oXjISNjfoqgJyFITvRaej3qehn1XWLminlO3bBxoauItSg+FG4HPA18LPPcx4C7n3KfM7GP+478BLgFO938uAL7g/ysiIiITw1uAQ80uxHiX1NJo77hgwZh5E4Cab0zjBoLFwcryudN54Ilnj2mNXPfYENMWD465CW9mgFlN63lxIFBueb9GBXE9W/eQyx/NAWbMRr8rY14yoeCOTbsallCoJRCv5thE+Uxv/yC3bBzg24HJCzMZyh6bWo5jlGA36tj9alY9CX7HmkuXccemXVyyfE7Dr7uw8lYa3hBnFYdy+5TknC1RvivOajTZjHFz70BThqoVVEwoOOd+bmaLip6+DEaHC30V6MZLKFwGfM0554AeM5thZnOcc7sSK7GIiIik2U3OuRVm9l/OuXc3uzDjUdyAvdJSasU3uLXejFZzg128JFuwC3Ln1I7RIRHrtvfUdMOcVBf6cmPb4yRPyi3v14ihG739gzyx7xDtbRlGRvJk/Mk5C9919cuX8MWfbx19/yXL5yRehjDVBuLV9GyI8pnioTngTV748lOzsc7vuMcxSuKs3HtK7RvESxz29g+O9qTZsH0vZ54yrSlJhbCeSGHDG5Lq5RO2ncL105bNJNazK+4yljv3HWLtPTuaOoFjtXMonBxIEjwJnOz/Phd4PPC+Af+5YxIKZnY1cDXAggWNyXSKiIhI3XWY2TuAl5jZ5cUvOudubUKZxpW4N8lxl1KrVTWBYHGwEAzAgkMiarlhTqILfSEhsXPfoTGzzxfPPxCly3Pn1A4GDw6VTUbUsxdGMEDOmDer/58GJuID+NjrXsCCWceNtkw3qncCVB+IVxNERvlM8DwEL5kwqT3DS+e2R9qXRgV5xUmzUt3jb904EOs6SPuqA/UeqlFqO8WJ2refv2BMorYaceo5OKQjzhwQ9VDzpIzOOWdmoWMly3zueuB6gJUrV8b+vIiIiKTS+4F3AjOANxS95gAlFGoU9yY57lJqtUqiZb04QCjM9F7LDXOtQVFxAFFolTQz8s7F6pEBpdefb4Te/kHWPTbEpMGj8wDkHPx0y9P86StPO+b977hgQUMTCUHVBOJJzikQ9p5soKdPcPWRZiuVNCveN4PY10ErrZxSnFBJopdPqe0EJ8PM5R2nzphS8zVcqZ7TuupKtQmFpwpDGcxsDvC0//wTwPzA++b5z4mIiMgE4Jz7JfBLM7vXOfflZpdnPKpm/HOUscZJ3pAm1SJbKNeaS5excdMWrrz4vKq3W83NelAwIZHLO952/nzmzpjC/kPD3PDLbYCLHGzFSW5EXRUiimCvhPbs42QyNjp/Qt651LU8V+vyFfMw/98o+1PLXAXB1UeardR5Vbx6AxCrRTt4DUYZ1tMoYcuclkrURRnyUumaKt5O3CRLue8IvhZ2HpbrYdXIHjClVJtQuB14D/Ap/9/bAs9/yMy+iTcZ47OaP0FERGTiCAxzGNSQh8qqDQ6jTuxV6Fa/asms0KXUGjGbfrWT5AXL9ZEVHTWVq1zQGKUOSi35Bl5Pg1zekc0Yay5dFqmMUYOROKtCRBHstp/LOy56wcncveXp0eEajWh5rudSkMX1dbl/jKKoda6CWiVRL1FXb4iakGzWShuV6iKsXNX0Qoo6f0a5ngGdUztGJ4+NU96w14JLgxakedhJlGUj1+JNwDjbzAaAj+MlEr5tZn8C9ANv9d/+Q7wlI/vwlo38ozqUWURERNKreJhDkIY8BNR6sx520108eVzGKNtaF7xRPTLsjbFO8ka12v0svoHesjdXc1nCAsIoN+vluj07wDnH4MGhyOWIEtTFWRUiikLAWRg+8v5Xnsb7X3laXQL8OC3ISUlz0FXOTet3RJ6Ho5yo51XUVvvP3flIw+szyjkSdpyLh6Xs3HeI3v7BY/4+lptjongfK/UMgPLDlyrVY9RzNs3DTqKs8nBlyEuvKvFeB1xTa6FERESkNTnn1JgQUanJ0qIGdlGWMStMUBXlRrUtYwzlHA64uXeAK2qcXKzcfhZa8uIsjdbeluGsmdlEyhPlu8Ju1qvt9hzWwgmUbdmMsypEJeWGj9S7p0AtLchxpDnoCtPbP8ia2zYx4g89GYpQL/Xs5QFHExy5vBtNSobVZ9JliXKOlOuJERzetfaesRPQRpljongfK5Wn3OtR6jHO355mz5UQpuZJGUVERESKmdnJwD8DpzrnLjGzs4EXa16Fo2pZS7zcTWywFTpP+WAAvBvVt6ycz03rd3hd4XPJBnrFN8yFZSDjLI22asmsuk5+V+3NepTPhQXXUVpiS23/zFOmNX34SCVRWpDrEfCnOegK07N1z+g8FgAZs7L1Erf7fDXzbAQTHAa8dOlsPnzxGSWHAiTd4yTKOVLuOBcSVyO5Y8+/KHNMRE1eVHo9aj3GOWebPVdCGCUUREREpB5uBP4T+Hv/8SPAtwAlFHzBG8m4a4mXu8ktHtsbnEgtrDXx8hXzElt6rPg7im+Yq1kaDaJNfldLa2m1N+uVPhe2v1HroXj71ZSzHsNHyimVRLru7r6yk84lJa1BV5hVS2Yxqd1LAGYyxrWXLa86mZhED5DiBEc2YyWTCUl9X7GocxOUO85hfx+jzjERVp5S52zY63HqsdXO2WJKKIiIiEg9zHbOfdvM/hbAOTdiZvWNYlpQ4UYy7lriUW5y47QmVtpe1EC93Czrwc/Vo5W6WRPIVRI3uKlHd/ZGDh+BY4PC4okkS006N55EOYZRZvYvpVwysZYhOMFtRE1w1KvHSZS5CSp9vlSd1tITKW7PgbiJoqB6D2lJmhIKIiIiUg/Pm9ksvIkYMbNVwLPNLVJ6VXOjG7dVq1JrYtj2goF6W8Z4y8r5ocvxVTu5YRLSOiFfnOCm3PCIWuorbPhIPQOXwvlUmLgybcelXqKuGhBlZv9SKnX3r3YITpxtVPPeuGq9nsP+njWqN0C1dZPWxGg5SiiIiIhIPfwl3nLSp5nZr4ATgTc3t0jpVnyjW0uwV+qz1bYmBm/sh3KOm9aPnegs+J1P7DtEWzZDLhdvcsMkpHlCvqj7GzaBZdQAo9w5Uzx8pFGBS5qOSyNafqMEwvUKlsu9Vtj3nfsORUr6JdV6X1yGWzYOYBCalCxo1nmT5DmSxPCkVkjAKaEgIiIiiTGz84DHnXMbzeyVwJ8CVwA/BgaaWrgWUkuwV27YQTUtZoUb+8IylI5jb3SLezG8/fwFFQOGpNWztTSOWoaHlAqiogYYcc+ZRgUuaTouaUigxEm8Vav4HAzue8bAzMjgGh6oX3n9bxjKefMKfKd3gLXvCz8GzThv0tA7IE0JuKiUUBAREZEk/Qdwsf/7S/AmZfwz4BzgetRLIZJagr1yn62mxaxwY19Yiq1UEBT8zlzeceqMKU0JHBvVnTlMnICk0ozzhQnpOqd2RAow4p4zjQxcmn1cIB0JlEqJtzjJqHLznRSfg8F99+YJdLRljDWXLmvYcenZuofh3NFJCqMcg0afN2noHZCWBFwcSiiIiIhIkrLOub3+728DrnfO3QLcYmb3Na9YraWWYK8egWLhxv6KFfNK3ujG+c5mTzhWKRirpWxxApJyM87D2GEOay5dNma1jjjbC1OvwKVQh8UrjCS57Wq2We8ESqnVTYoFz4+RnGPH3oNjPh8lGVXpfaXOweJeRgDOOQYPDiVaB+WsWjKL9qyN9lBIY+t7WnoHpCEBF4cSCiIiIpKkrJm1OedGgFcBVwde031HRLUEe/Vs4So30VmU72x2l+Jykx4WemCM5KovW5SAJMrs/sVB4eDBoYqT9lVz3JMOXAr1WwhcM0Zix7nWc6ee10XUshXOj6HhPHngV3272bB97zG9CIZH8tyycSDSuVGctCp1DkbpZVTNPsc919Ze/eLIcyg0Qyv2DkgD/ccuIiIiSVoL/MzMdgOHgF8AmNlStMpDLLUEe81o4YrynT1b94wGm0PD1Xcprralutykh8HW22q7O1cKSKLO7l9tS2mzWzYL9VuoxyS7jidx7tSrfqL2TCmcH5+78xF+1be7ZC+C4ZE82YyFJrcqnRth52ClXkZxVJvcafb5GUUrlDFtlFAQERGRxDjnPmlmdwFzgB875wqxRQZvLgVJuXoOSeic2nE02PQfV1O+aluqy016WCiXUVt37HIBSdzAs54tpfU4zsUt8BlLrmt7EudOvQT328zKlq1rYScfvvgMNmzfW7IXQWElhrX37AidB6XSuVHuHEwiYE7DXAMTRbOHiEWhhIKIiIgkyjnXU+K5R6J81sy+AlwKPO2cWx54/s+Aa4Ac8APn3EfNbBHwEPCw/7Ye59z7/fd3ATcCU4AfAv8jkNyQEPUekjB4cIiMeS3XGaOqMdy1BDNhwViwZfgtK+fXrTt2nJ4H9Wwp7RvM8Zm7kj/OwfpNeg6FJM6deula2MmaS5ex5rZN5J3j2nWbOfOUaWWD+nK9CApDcMr1QmhmcJmWuQbGu2YPEYtKCQURERFJkxuBzwNfKzxhZhcClwG/55w7YmYnBd7/mHPunBLb+QLwPmA9XkLhtcAddSpz4hrdKhVnjfoo2wkrdxKBSK3bKA7GGjluOi1jtLfszdWthblewW69gtikrrXBg0PknYtcp5V6EaThPAlTXD6A6+7uS2VZW1mr9ARRQkFERERSwzn3c7/nQdAHgE85547473m63DbMbA5wQqGnhJl9DXgjLZJQaHSrVPFSdm3ZTFWTtkUpdxKBUj2CrUa2+Da7dRngrJlZOtpyLdXCXI/jntS11ts/yBP7DtV07YT1VkirYG+KVmhFb0Wt0hNECQURERFJuzOAl5vZJ4HDwEeccxv81xab2W+B54B/cM79ApgLDAQ+P+A/dwwzuxp/JYqTTz6Z7u7uRAt+4MCB2Ntc99jQmMnn1t65gf2n1W+8ePD7RnKOV87LMGtKO2fNzLJ/2/10b0um3H2DObbszXHWzCzLOrPs3zYQedtwbF0uM2Jvo1GC+7q0M9vs4hzjlPZDfGTFlNEyxjnOjVZcl0ke9ySutd/tPMDnf/xrhvOQNXjF3DZeOrctcp32Deb49IbDDOehPQMfPW9yKs+ZMFHrMMo1Uc3fy6Sl7dr9yIqOqq7TRtalEgoiIiKSdm3ATGAVcB7wbTNbAuwCFjjn9vhzJnzPzJbF2bBz7nrgeoCVK1e61atXJ1rw7u5u4m5z2uJB1m3vGW2VuvLi80q2+CXVVbv4+z70hguq2l65cvf2DwbG7OeqasWspi7D1HNISRL7mrTi/e3u7uaqhM/1JBXK2zm1g8/ctbludRn1Witn3Zd/zHB+GIc3v8PKZadxVYUlPoM2393HiHsYB+QcHJmxkNWro3++2aLUYdRrIslrPK5jl45Nx7W7usrPNbIulVAQERGRtBsAbvUnVbzHzPLAbOfcM0BhGESvmT2G15vhCWBe4PPz/OdaQpSu3Ul2M06qK3m57aRpLHC9u2inaV+h9P6mWbC8GbNY8xLElcS5f3y7jVl9Yv+h4Vifb5Vu7WGi1GE9rokkk4KFcy6JpWOT0AorOwQpoSAiIiJp9z3gQuBuMzsD6AB2m9mJwF7nXM7vsXA6sNU5t9fMnjOzVXiTMv4h8H+bVPaqVBo/nfQNelLjtcO2k6agqd4Bf5r2FUrv7zJrapHKCpYX58hkDMPVrS5rPfcPDDsMRgPRG365jd9fdkrNK4+0kkp1mPQ1kXRSMOmlY2vRinNSKKEgIiIiqWFma/F6ec42swHg48BXgK+Y2SZgCHiPc86Z2SuAa81sGK9x8P3Oub3+pj7I0WUj76BFJmSMKm1BayVpCprqXXdp2lcovb/7tw1U/mBAI1tMi8u75tJliS4/mbSzZmbJZkYYyXvhaN652EmqtE/AWEqccyLpayLppGDwnCu3dGwjroO09XCKQgkFERERSQ3n3JUhL72rxHtvAW4J2c69wPIEi5YqaQtao0hL0NSIuqvHvlYbzJTa3zgTGja6xbTVzu2lnVmuvWw5a27bRN45OlogwVeras6JJK+JpJOCjR5mVk6rJYtBCQURERGRlpSWAL0VtVrd1RrM1LK/zWgxTfvxCSZ3AN5xwQLOPGVayyRBatXsVvRmLB3bqH1utYQaKKEgIiIiIpJqlYKZenbFbsUW07ii1F9w5Ylr120eTe58ZEVH1TPxt6p6nBO9/YOse2yIaYsH6VrYWfGYNDrpFGefa70ei/ct7ZM0KqEgIiIiIhIiDTfz5YKZenfFbsUW0zgq1V/xcoLFK09s2ZtryYn0apH0ORFcZWHd9h7WXLpsTNImDfUZdZ+TPhda4dxSQkFEREQkpdIQzKZdLXVU6bNpuZlv9pKcaR+CUIty9VdqOcHilSfOmplt+hCAZoh7TpS71oKrLAyP5Llj065U1meUfU76XGiFc0sJBREREZEUSkswm2a11FGUz6bpZr4VluRshKSTbOXqr9Rygh3tY1ee2L/tfqYtHp/HIKm6rnStFY7B0LBXf5csn8OG7Xtbsj6LV4zYue8Qvf2DVddfK1zfSiiIiIiIpFCagtm0qqWOonw2TTfzYcHdeB+SEFSPJFu5+isODlefeRInTpvEmadMG31f97bxeQySrOtK11qh/tbeuYErLz6ProWdLTvJZWFfCsNk1t6zg1s2DlRdf61wbimhICIiIpJCaQpm0yqsjqK0rEap37g38/UaolIquAPGfFcaA42k1SvJFlZ/weMfnIyxVIA43o5BknUd9Vrbf1rH6He0cn12LeykZ+seRnLJTKSa9rpQQkFEREQkheJMApbm1qt6KlVHUVtWo9Zv1Jv5eg5RKQ7ubtk4wK0bBybccJhmJNkKx/+6u/smVI+hJOu6FVrZk9bMiVQbTQkFERERkZSqFMyOtxvTahTXUZyW1SRb/uo5RKU4ODGYUMFtQTMD03otlZjWIDvpuk5DK3sj67vZE6k2khIKIiIiIi1qvN2YJqFZQ0Xq+b3FwQnALRsHJuRwmEYGpsUBaD2WSqyUDGxm0iENSYBKCst6GnD5inkNW84xiokykaoSCiIiIiItJBhg1HpjmuYW0mo1qxW73t9bHJxMtC7kjRYWgDayR4t6IJXX2z/Ildf/hqGctw7Hd3oHWPu+0nWUpuTreBsCooSCiIiISIsoFWCsuXQZd2zaxSXL58RaMrEwC/lIbvwFK81qWW3k97ZC63ErqzYAjZqki5IMTFMQnEY9W/cwnHOjj8vVUdp6BYyn61cJBREREZEmitNLoNzkfBu27x2znF2573vnDT0cGc5TuBWvJlgZj70bRAqqCUDj9CiI0kqdtiA4bVYtmUV71kZ7KJSroyj1rb9p1VFCQURERKRJ4nZpTmJyvkJSopBMMMrfiCdRbhm/xmsQVk239Lg9Ciq1Uo+3rvFJ61rYydqrXxxpDoXC+9M0x8J4oYSCiIiISJNUE4DUOjlfMCmRzRhvWTm/4o14reWW8Wm8B2Fxu6XXo0dBLV3jx2uyJyipoQP6m1Y9JRREREREmqSaAKTWyfmSaPUcb12xJ0LgVQ8KwsZKU4+CRiR76n3dNPK6HG9/0xpJCQURERGRJkkiAKmmha7WVr00BU61Gu+t7PWkIOxYaZlsr97JnnpfN42+LsfT37RGU0JBREREpInSEoDE1arlLjZeWtmb0ctCQVh61TvZU+/rphnX5Xj5m9ZoSiiIiIiIpIi631eWZB21eit7s5cAVRCWTvVO9tT7umn163IiUUJBREREJCXU/b6ypOuolVvZk1oCNO2UZKtOPZM9tV43lY5pK1+X9Za260EJBREREZGUGC/d7+upHnVUKfBK2w18QRJLgKadkmzpVW3CIuoxVe+XY6Xxesg09dtFREREZFShm2/Wxl9gmJRG11HhBv6zP36Yd97QQ2//YF2/L45gXXRkjXdcsCAVAUaSSiWQpLXpmFYvjXWnHgoiIiIiKaFuvpU1uo7S3GtkIpwvGks//uiYVi+NdaeEgoiIiEid9A3m2Hx3X6xgT918K2tkHaXxBj5ovJ8vEyFpMtHomFYvjXWnhIKIiIhIHfT2D/LpDYcZcQ+nZqyrxJfGG/hWF3dOivGeNKm3NM4BomNavbTVnRIKIiIiInXQs3UPw3lwpK+r/ERTa0CVthv4VpbGSeXGM9W31JsmZRQRERGpg1VLZtGeQRMsNllv/yBXfqmHz/zoYa78UromVaxVb/8g193dl7p9KleuNE4qN56pvqXe1ENBREREpA66Fnby0fMmc2TGwlR1NZ5obt04wNBIHoChkTy3bhwYF8cirS3PlcqV9jkpxpvOqR1kzACn+pa6UEJBREREpE6WdmZZvXpps4sxobkKj1tVWlefqFQuzUnROL39g1y7bjO5vCObMdZcukz1LYlTQkFERESkydI4adp4ccWKedx87+MM5xztWeOKFfOaXaREpLWlP0q5NCdFYxSSOw5wzjF4cKjZRZJxqKaEgpltB/YDOWDEObfSzGYC3wIWAduBtzrn0jWwS0RERCQlynURV6Khdl0LO1l79YvHXT2mtaU/reWaiNKadJLxJYkeChc653YHHn8MuMs59ykz+5j/+G8S+B4RERGRcSesi3hax8inXakkzHhtEU/rfhWXS4mx5lByRxqhHkMeLgNW+79/FehGCQURERGRksJaEdM6Rj7NlIRJn0Yfk/GYvKhln9KadJLxo9aEggN+bGYO+A/n3PXAyc65Xf7rTwIn1/gdIiIiIuNWWCuiuivHpyRM+jTymIzHhFIr7dN4TOZIZbUmFF7mnHvCzE4CfmJmW4IvOuecn2w4hpldDVwNsGDBghqLISIiItK6SrUiqrtyfErCpE8jj8l4TCi1yj61UuJDklVTQsE594T/79Nm9l3gfOApM5vjnNtlZnOAp0M+ez1wPcDKlSvHywo+IiIiIrGFteypu3I8SsKEC55jjdTIYzIeE0qtsk+tkviQ5FWdUDCz44CMc26///urgWuB24H3AJ/y/70tiYKKiIiIjEdq2UuWkjDHKj7HPrKiY3TCs0Zo1DEZjwmlVtmnVkl8SPJq6aFwMvBdMyts5ybn3H+b2Qbg22b2J0A/8NbaiykiIiIyPqllb3xJ4zjy4nNsy95cs4tUN+MxodQK+9QqiQ9JXtUJBefcVuD3Sjy/B3hVLYUSERERmSjUsjd+pLW3SfE5dtbMbLOLJONQKyQ+JHn1WDZSRERERCJSy148aewBUJDW3ibF59j+bfc3u0hSQZrPc5EgJRREREREmiQYNFxz4dJmFyf10tgDIHgM09zbJNh63L2tyYWRstJ4nouEUUJBREREpAkUNMSXth4ApY6heptIrdJ2nouUo4SCVLToYz9odhFEpI5a4Rrf/qnXN7sIIolT0BBf2noAlDqG11y4VMdRapK281ykHCUURERERJpAQUN8aZtvQsdQ6iFt57lIOUooiIiIiDSBgobqFM8k38zJ63QMpV60YoK0CiUURERERJpEQUNt0jAPhY6hiExkmWYXQERERESkGqXmMBBphN7+Qa67u4++wVyziyLSVOqhICIiItIitDb9WJrDQOop7HoL9oxpMzh3xaCuR5mwlFAQERERaQFp6N6fNprDQOql3PUW7Bkz4tAKLTKhaciDiIiISAtQ9/7SuhZ2aqlGqagwRKG3fzDS+8tdb4WeMVmDtgzqGSMTmnooiIiIiLQAde8XqU41vXvKXW/BnjGT9vUrmSUTmhIKIiIiIi1A3ftFqlOqt0Gl66fS9VZY3aO7e6CeRRdJPSUURERERFqEligUia/a3j263kQq0xwKIiIikhpm9hUze9rMNhU9/2dmtsXMNpvZpwPP/62Z9ZnZw2b2msDzr/Wf6zOzjzVyH0Qg/ph9qZ9Cb4O/fPWZmsxUJGHqoSAiIiJpciPweeBrhSfM7ELgMuD3nHNHzOwk//mzgbcDy4BTgTvN7Az/Y9cBvw8MABvM7Hbn3IMN2wuZkArLDHZO7eDadZu1IkeKqLeBpFkrLwmshIKIiIikhnPu52a2qOjpDwCfcs4d8d/ztP/8ZcA3/ee3mVkfcL7/Wp9zbiuAmX3Tf68SClI3wYn/MmbknYs1Zl9EJqab1u9gzW2byDvXkglIDXkQERGRtDsDeLmZrTezn5nZef7zc4HHA+8b8J8Le16kboIT/+XzjowZWUMrcrQIDVFpnkLd9w3mml2UhuvtH2TNbZsYyXsJyKEWXBJYPRREREQk7dqAmcAq4Dzg22a2JIkNm9nVwNUAJ598Mt3d3UlsdtSBAwcS3+ZElfa6nLQvR5vBiIO2DLzjrHYODDvOmpll/7b76d7W7BIelfa6bLS+wRyf3nCY4Ty0Z+Cj501maWc20mdVl7UJ1n1bxgF3Ra77KNvesjfHWTOziW0zaeseG2Ik70YfGzBpX3/Nq4c08rxUQkFERETSbgC41TnngHvMLA/MBp4A5gfeN89/jjLPj+Gcux64HmDlypVu9erViRa8u7ubpLc5UaW9LlcD565ojXHQaa/LRtt8dx8j7mEckHNwZMZCVq9eGumzqsvajKn7vMWq+3J6+wf5zF09/jwmudQOI5i2eJB123sYGs6TyRjXXracd1ywoObtNvK8VEJBRERE0u57wIXA3f6kix3AbuB24CYz+1e8SRlPB+7Ba+Q53cwW4yUS3g68ownllglGE/+1pmqXlayklSfaa5Rg3WeNxOo+OAQpzfOYFFYgaeXzRAkFERERSQ0zW4vX2DvbzAaAjwNfAb7iLyU5BLzH762w2cy+jTfZ4ghwjXMu52/nQ8CPgCzwFefc5obvjIi0hHoEdcFJOltxor1GCdb9pH39idVRvZJE9dDqiUglFERERCQ1nHNXhrz0rpD3fxL4ZInnfwj8MMGiicg4lnRQ1yot5GlQqPta5w0o3mYwSQRw3d19LdsLIM2UUBARERGRlqAu5NIqWqmFfLwqJCrUW6S+lFAQERERkdTrG8wFJllTUCDpNh7Gxo8X6i1SX0ooiIiIiEjqbdmbU1AgLaXVx8aPF+otUl9KKIiIiIhI6p01M0tHW05BgYjEot4i9aWEgoiIiIik3tLOrIICEamKeovUjxIKIiIiItISFBSIiKRLptkFEBEREREREZHWo4SCiIiIiIiIiMSmhIKIiIiIiIiIxKY5FJpo0cd+0OwiiIiIiIiIiFRFPRREREREREREJDYlFEREREREREQkNiUURERERERERCQ2JRREREREREREJDZNyigiIiISYnh4mIGBAQ4fPlzV56dPn85DDz2UcKnCTZ48mXnz5tHe3t6w7xQRkYlLCQURERGREAMDA0ybNo1FixZhZrE/v3//fqZNm1aHkh3LOceePXsYGBhg8eLFDflOERGZ2DTkQURERCTE4cOHmTVrVlXJhEYzM2bNmlV1bwoREZG4lFAQERERKaMVkgkFrVRWERFpfUooiIiIiKTcJz/5SZYtW8aLXvQizjnnHNavX8/q1au59957m100ERGZwDSHgoiIpN6ij/2g2UWIZPunXt/sIsg49Jvf/IZ169axceNGJk2axO7duxkaGmp2sURERJRQEBEREUlSb/8gPVv3sGrJLM6YWfut1q5du5g9ezaTJk0CYPbs2TVvU0REJAka8iAiIiKSkN7+Qd55Qw+f/fHDvPOGHu4beK7mbb761a/m8ccf54wzzuCDH/wgP/vZzxIoqYiISO2UUBARERFJSM/WPQyN5Mk7GB7Jc2//vpq3efzxx9Pb28v111/PiSeeyNve9jZuvPHGmrcrIiJSKw15EBEREUnIqiWz6GjLMDySp70tw8qFMxLZbjabZfXq1axevZoXvvCFfPWrX01kuyIiIrVQQkFEREQkIV0LO/nGVasSnUPh4YcfJpPJcPrppwNw3333sXDhQjZt2lTztkVERGqhhIKIJC5/5HkAMpOOa3JJREQar2thJ10LOwHYv39/zds7cOAAf/Znf8a+fftoa2tj6dKlXH/99bz5zW+uedsiIiK1UEJBRBKnRIKISHK6urr49a9/fczz3d3djS+MiIhIgCZlFJHEPf/Qz3n+oZ83uxgiIiIiIlJHdUsomNlrzexhM+szs4/V63tEJH32//aH7P/tD5tdDBERERERqaO6DHkwsyxwHfD7wACwwcxud849WI/vE5F0Oektn2h2EUREREREpM7qNYfC+UCfc24rgJl9E7gMaFhCYdHHftCorxKRIpn2yc0ugoiIiIiI1Fm9EgpzgccDjweAC4JvMLOrgasBFixYkHgBtn/q9YlvU0Si+frXvw7Au971riaXRERERERE6qVpkzI65653zq10zq088cQTm1UMEamDG264gRtuuKHZxRARERERkTqqVw+FJ4D5gcfz/OdEZAL4yU9+0uwiiIiMG8cffzwHDhwYfXzjjTdy77338vnPf76JpRIREalfD4UNwOlmttjMOoC3A7fX6btEJGXa29tpb29vdjFERERERKSO6tJDwTk3YmYfAn4EZIGvOOc21+O7RCR9brzxRgDe+973NrUcIiIiE01v/yA9W/ewasksuhZ2Nrs4IjLO1WvIA865HwJaiF5kAlJCQUTGpTs+Bk8+EOsjU3IjkC1zu3XKC+GST5XdxqFDhzjnnHNGH+/du5c/+IM/iFUOmRh6+wd55w09DI3k6WjL8I2rVimpICJ1VbeEgohMXN3d3c0ugojIuDFlyhTuu+++0ceFORREivVs3cPQSJ68g+GRPD1b9yihICJ1pYSCiIiISBQVehKUcmj/fqZNm1aHwogca9WSWXS0ZRgeydPelmHVklnNLpKIjHNKKIiIiIiIjANdCzv5xlWrNIeCiDSMEgoiIiIiIuNE18JOJRJEpGGUUBARERFJsQMHDox5/N73vleT3oqISCpkml0AEREREREREWk9SiiIiIiIiIiISGxKKIiIiIiIiIhIbEooiIiIiJThnGt2ESJrpbKKiEjrU0JBREREJMTkyZPZs2dPSwTqzjn27NnD5MmTm10UERGZILTKg4iIiEiIefPmMTAwwDPPPFPV5w8fPtzQAH/y5MnMmzevYd8nIiITmxIKIiIiIiHa29tZvHhx1Z/v7u7m3HPPTbBEIiIi6aEhDyIiIiIiIiISmxIKIiIiIiIiIhKbEgoiIiIiIiIiEpulYdZiM3sG6G92OVrYbGB3swsxTqlu60v1Wz+q2/ppdN0udM6d2MDvm7DqdD+iazE5qsvkqC6To7pMjuoyOfWoy5L3I6lIKEhtzOxe59zKZpdjPFLd1pfqt35Ut/WjupU4dL4kR3WZHNVlclSXyVFdJqeRdakhDyIiIiIiIiISmxIKIiIiIiIiIhKbEgrjw/XNLsA4prqtL9Vv/ahu60d1K3HofEmO6jI5qsvkqC6To7pMTsPqUnMoiIiIiIiIiEhs6qEgIiIiIiIiIrEpoZBSZvYVM3vazDYVPf9nZrbFzDab2acDz/+tmfWZ2cNm9prA86/1n+szs481ch/SKk7dmtkiMztkZvf5P18MvL/LzB7w6/b/mJk1el/SplTdmtm3AvW33czuC7ym8zaiOHWr8za+kPo9x8x6/Dq818zO9583v+76zOx3ZrYi8Jn3mNmj/s97mrEvkg76OxaPmc03s7vN7EH//+H/4T8/08x+4l9TPzGzTv/50OtQPGaWNbPfmtk6//FiM1vv19m3zKzDf36S/7jPf31RUwueMmY2w8xu9u8RHzKzF+u8rI6Z/YV/fW8ys7VmNlnnZXQh9yqxz8XE71Wcc/pJ4Q/wCmAFsCnw3IXAncAk//FJ/r9nA/cDk4DFwGNA1v95DFgCdPjvObvZ+9bsn5h1uyj4vqLt3AOsAgy4A7ik2fvW7J9SdVv0+meBNf7vOm/rV7c6bxOoX+DHhfoBXgd0B36/w6/DVcB6//mZwFb/307/985m75t+mnI+6e9Y/DqbA6zwf58GPOL/P/Fp4GP+8x8D/rf/e8nrUD9j6vQvgZuAdf7jbwNv93//IvAB//cPAl/0f3878K1mlz1NP8BXgav83zuAGTovq6rHucA2YIr/+NvAe3VexqrDUvcqsc7FetyrqIdCSjnnfg7sLXr6A8CnnHNH/Pc87T9/GfBN59wR59w2oA843//pc85tdc4NAd/03zuhxazbksxsDnCCc67HeVfn14A31qG4LSWkbgEvUwq8FVjrP6XzNoaYdVuSzttwIfXrgBP836cDO/3fLwO+5jw9wAy/bl8D/MQ5t9c5Nwj8BHht/UsvKaS/YzE553Y55zb6v+8HHsILQC7DC+jw/32j/3vYdSiAmc0DXg/c4D824CLgZv8txXVZqOObgVep95rHzKbjBXFfBnDODTnn9qHzslptwBQzawOmArvQeRlZyL1K3HMx8XsVJRRayxnAy/1uPz8zs/P85+cCjwfeN+A/F/a8HCusbgEW+10Gf2ZmL/efm4tXnwWq28peDjzlnHvUf6zzNjnFdQs6b5PwYeBfzOxx4DPA3/rP69yVSnQu1MDv2nwusB442Tm3y3/pSeBk/3fVcXmfAz4K5P3Hs4B9zrkR/3Gwvkbr0n/9Wf/94vWgfAb4T///1BvM7Dh0XsbmnHsC7//SHXiJhGeBXnRe1iruuZj4OaqEQmtpw+uesgr4a+DbEz1Tl6Cwut0FLHDOnYvfddDMTgjfjJRxJRVa0KVqxXWr8zYZHwD+wjk3H/gL/BYqEakfMzseuAX4sHPuueBrfs8qLU9WgZldCjztnOttdlnGgTa8LuZf8P9PfR6vW/konZfR+GP7L8NL0pwKHId68SWqWeeiEgqtZQC41e+6cg9e1nk28AQwP/C+ef5zYc/LsUrWrd8dfw+A/x/zY3i9GZ7Aq88C1W0Zfte2y4FvBZ7WeZuAUnWr8zYx7wFu9X//Dl43dtC5K5XpXKiCmbXjJRO+4ZwrXHtPFbqM+/8WhiSqjsO9FPgDM9uON9zmIuDf8bo8t/nvCdbXaF36r08H9jSywCk2AAw459b7j2/GSzDovIzvYmCbc+4Z59ww3v+vL0XnZa3inouJn6NKKLSW7+FNHoiZnYE3Mcxu4Hbg7f5sqIuB0/EmXtsAnO7PntqBN6HJ7c0oeAv4HiXq1sxONLOs//wSvLrd6nctes7MVvk9Gf4QuK0pJW8NFwNbnHPB7vY6b5NxTN3qvE3MTuCV/u8XAYUhJbcDf+jPoLwKeNav2x8BrzazTr8l5tX+czLx6O9YTP7fpC8DDznn/jXw0u14yT38f28LPF/qOpzwnHN/65yb55xbhHfu/dQ5907gbuDN/tuK67JQx2/2368Wd8A59yTwuJmd6T/1KuBBdF5WYwewysym+td7oS51XtYm7rmY/L1K3Fkc9dOwWTzX4nVbHsbLjv4JXpD7dWATsBG4KPD+v8drhXyYwKzteDN8PuK/9vfN3q80/MSpW+AKYDNwn//8GwLbWem//zHg84A1e9+a/VOqbv3nbwTeX+L9Om/rULc6b5OpX+BleOM778cby93lv9eA6/w6fABYGdjOH+NNMNoH/FGz90s/TT2n9HcsXn29DK+r7u/8v133+XU4C7gLL6F3JzDTf3/odaifMfW6mqOrPCzBS9z34fW6KqxsNdl/3Oe/vqTZ5U7TD3AOcK9/bn4Pb2Z8nZfV1eX/BLb49yH/hbfSl87L6PVX6l4l9rmY9L2K+RsVEREREREREYlMQx5EREREREREJDYlFEREREREREQkNiUURERERERERCQ2JRREREREREREJDYlFEREREREREQkNiUURJrMzP7ezDab2e/M7D4zu6DMe280szf7v7/c/9x9ZvYCM9sU4bveYmYPmdndSe5Dhe98o5mdHXh8rZld3KjvFxERmajMLOffJ2wys++b2Ywqt/NeM/u8//v7zewPq9jGajN7SeBxVdtpJDPbbmazm10OkTRra3YBRCYyM3sxcCmwwjl3xP9PqyPix98J/H/Oua+b2aKIn/kT4H3OuV9GLF+bc24k4rbDvBFYBzwI4JxbU+P2REREJJpDzrlzAMzsq8A1wCdr2aBz7otVfnQ1cAD4dY3baQkJ3UOJpJ56KIg01xxgt3PuCIBzbrdzbqeZdZnZz8ys18x+ZGZzgh8ys6uAtwL/ZGbfKHrtvWZ2q5n9t5k9amaf9p9fA7wM+LKZ/YuZTTaz/zSzB8zst2Z2YeDzt5vZT4G7/MffM7Of+Jn6D5nZX/qf6TGzmf7n3mdmG8zsfjO7xcym+i0RfwD8i99CclpRL4tX+dt5wMy+YmaT/Oe3m9n/NLON/mtn1e8QiIiITAi/AeYC+P8f/7d/n/GLwv+zZvYGM1vv/998p5mdXLwRM/uEmX3EzE71/28v/OTMbGGpbfgNH+8H/sJ/78sL2/G3eY5/T/E7M/uumXX6z3eb2f82s3vM7BEze3mJ8qz233ezmW0xs2+YmfmvjfYwMLOVZtYd2Iev+vveb2aXm9mn/XuO/zaz9sBXfNR//h4zW+p//kT/XmeD//PSwHb/y8x+BfxXIkdNJOWUUBBprh8D8/3/JP+fmb3S/0/s/2/nXkKtquI4jn//amVhaYMm0UOSFMzH9RVlRS9nQS+Ni0hoSIMQrUhrVKk4aZKKYVEWGkGK5EAHPSRUtKJQUi8+IVR6OJFKNC1Mfw3W/+TudF/nXOOi/T6Tvffa6/z32hsue93/XmstA6ZIGge8R93XBEkrgPXAPEnT2onbArQCI4HWiLhR0kJgOzBN0jzKVwpJGglMBVZFRP/8/di8/r15PAJ4HJiQbTklaQylc1IbrrhO0gRJo4F9wExJX1ba2SLpu1oD81orgdZsQz/gmco9HJM0FngTmNu9x2lmZmb1IqIv8CDlnQzwNjA7+xlzgeVZvg24I9/xq4EXO4op6ad8t7cA7wAfSTrSXgxJh4G3gMX5m6114d4HXpI0CmgDXq2c6yfpduC5uvKqMXl+OHALcFfHT+NvQ4AHKB8+PgA2ZX/kNPBQpd7xLH8DWJJlS/NeJgCTgRWV+sOBSZKmdqMNZhc9T3kw60WSTkbEOOAe4H5gDbCI8g/8xkyw9wWONhj6c0nHASJiL3Az8H1dnbspiQsk7Y+II8DQPLdR0s+VupsknQBORMRxYEOWtwGjcn9ERCwCBgEDgE+7aOMw4JCkg3lcG4q5JI/X5XYHJZlhZmZmjbkyInZSRibso/QtBgATgbXZzwC4Irc3AGuijIy8HDjU1QXy6/zTlH5FwzEiYiAwSNKWLFoFrK1UqfYHBncQ5htJP2S8nVmvq+mdH0s6ExFtlL7WJ1neVnedDyvbxbk/CRheeX7X5HMFWC/pdBfXNrtkOKFg1ssknQU2A5vzpTYL2CPpzh6E/aOyf5bG/9Z/6yTeucrxuUrslcCjknZFxAzKXMmeqF2jmfabmZlZrqEQEVdREv2zKO/rX2trK9RZBrwuaX1E3AfM7yx4Jg3eBR6WdLKZGN3Qnf5AR/2ePzk/Irs//1SbbnouIs5IUpZX+zYAame/D2UUxu/VgJlgqO9DmV3SPOXBrBdFxLCIuLVS1EL5gnBdlAUbiYjLIuK2/+DyWykLOxIRQ4GbgAM9iHc1cDSnbFSnYZzIc/UOAINr8xGBJ4Et7dQzMzOzHpB0CpgDvACcAg5FxBMAUYzOqgOBH3N/emcx832/ljJV4WDlVEcx2u0P5IjKXyrrI1zI/sBhYFzuT24yRmtl+1XufwbMrlWIiJYmY5td9JxQMOtdAyhrF+yNiN2UeXevAFOA1yJiF7CTMjTxQlsO9MlREWuAGbXFIZv0MvA18AWwv1K+GpiXizMNqRVmVv8pypDLNsoXgUt6xWczM7PeIulbYDdl3aRpwMzsZ+wBHslq8ynv5R3AsS5CTgTGAwsqCzNe30mMDcBjtUUZ62JNpyzgvJvycWVhc3f5LwuApRGxnTJyoRnXZrueBZ7PsjnA+FxEci9lwUmz/6U4P7rHzMzMzMzMzKx7PELBzMzMzMzMzBrmhIKZmZmZmZmZNcwJBTMzMzMzMzNrmBMKZmZmZmZmZtYwJxTMzMzMzMzMrGFOKJiZmZmZmZlZw5xQMDMzMzMzM7OGOaFgZmZmZmZmZg37Cx152UsKKaF/AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=(18, 6))\n", "\n", "plt.subplot(1, 2, 1)\n", "plt.hist(O.SI)\n", "plt.plot(np.array([1, 1]) * O.H, [-5, 5], 'k:')\n", "plt.xlabel('Self-information')\n", "plt.title('Entropy = %3.1f' % (O.H))\n", "\n", "plt.subplot(1, 2, 2)\n", "plt.plot(O.SI, '.', label='SI')\n", "plt.plot(np.cumsum(O.SI) / (np.arange(1, 1 + len(O.SI))), '-', label='H')\n", "plt.legend()\n", "plt.grid()\n", "plt.xlabel('Realization number')\n", "plt.ylabel('Self-information')\n", "plt.title('Entropy as a function of number of realizations')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Entropy as a function of number of conditional data" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "n_cond = 1, H=2260.7\n", "n_cond = 2, H=1705.0\n", "n_cond = 4, H=1332.4\n", "n_cond = 6, H=1138.6\n", "n_cond = 8, H=1012.2\n", "n_cond = 12, H=852.4\n", "n_cond = 16, H=730.9\n", "n_cond = 24, H=615.3\n", "n_cond = 32, H=563.6\n", "n_cond = 64, H=522.3\n" ] } ], "source": [ "TI, TI_filename = mps.trainingimages.strebelle(di=4, coarse3d=1)\n", "\n", "n_cond_arr = np.array([1, 2, 4, 6, 8, 12, 16, 24, 32, 64])\n", "\n", "H = np.zeros(n_cond_arr.size) # entropy\n", "t = np.zeros(n_cond_arr.size) # simulation time\n", "i = 0\n", "SI = []\n", "for n_cond in n_cond_arr:\n", " O = mps.mpslib(method='mps_snesim_tree', simulation_grid_size=[80, 70, 1], verbose_level=-1)\n", " O.par['n_real'] = 20\n", " O.par['n_cond'] = n_cond\n", " # Choose to compute entropy\n", " O.par['do_entropy'] = 1\n", " O.ti = TI\n", "\n", " O.run_parallel()\n", " print('n_cond = %d, H=%4.1f' % (n_cond, O.H))\n", " SI.append(O.SI) # Self-information\n", " H[i] = O.H # Entropy\n", " t[i] = O.time\n", " i = i + 1" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'simulation time')" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAyoAAAFfCAYAAAChuTicAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAxOAAAMTgF/d4wjAACO7UlEQVR4nOzdd3hU1dbA4d+eNJIAoZcQSihSQkeaFwQMgg0B6YoUkWIHLFy9tmvvHQ1wQVBBujSBgHwIUhWQIhZ6CyWBAKmTNvv7YzJxiElImcmZst7nOU9yypyzdiA7s2Y3pbVGCCGEEEIIIVyJyegAhBBCCCGEECI3SVSEEEIIIYQQLkcSFSGEEEIIIYTLkURFCCGEEEII4XIkURFCCCGEEEK4HElUhBBCCCGEEC5HEhUhhBBCCCGEy/E1OoDSEhAQoKtWrZrnubS0NAICAko5ItfgzWUH7y6/s8oeExOTrrV2mx+q1A15k7JL2R1N6gb3ImX0HK5czuvVC16TqFStWpUzZ87keS46OprevXuXckSuwZvLDt5dfmeVXSkV5/CbOpHUDXmTskvZHU3qBvciZfQcrlzO69UL0vVLCCGEEEII4XIkURFCCCGEEEK4HElUhBBCCCGEEC5HEhUhhBBCCCGEy5FERQghhBBCCOFyJFERQgghhBBCuBxJVIQQQgghhBAuRxIVIYQQQgghhMuRREUIIYQQQgjhciRRsZOScoJfdg0gJeWE0aEIIYQQLu/IkSN06tSJI0eOGB2KIby9/EI4myQqAGSQmZnI0WMfkpCwl6PHPiQzM5GsrDSjAxNCCCFcjtlsJiEhgQkTJrBz505eeOEFEhISMJvNRodWKtLT0zl9+rTXll+I0uJrdACuwMd3Ips265z92NhVxMauAnyIvOWQcYEJIYQQLig4OBiLxZKzP3/+fObPn4+Pjw+ZmZkGRlY6+vXr59XlF6K0SIsKkJ76b8qX7YC2WH8c2mKifNmOtGmx1ODIhBBCCNeza9cuypcvn7Pv4+NDZGQk+/btMzCq0vP555/Tvn37nH1vK78QpUUSFeDYynoc2twAlAVLli8oC4c2N2De8+eNDk0IIYRwOTNnziQhIQGAgIAALBYL/fv3JyIiwuDISke9evXo2bMnYE1SvK38QpQWSVSAGh3iqRC+GVBcPhwJKCqEb+L28XWNDk0IIYRwKfPmzWPq1KmUK1cOk8nEE088gVKK2bNnGx1aqVq5ciUA7dq188ryC1EaJFEBgmuWx9dUjTNbHiJ2Xz/ObJ2Ar6kadZqHGx2aEF5PKXVCKfWXUmpv9jYkn+vGKKUOK6WOKqVmKKX8SjtWITzd77//ztixY6lRowbdu3dn6dKlvP322yxZsoTw8HC01te/iYMYWTdoralXrx4AHTt2NKT8QngDGUwPHFvhhyVrCMpUEUtmDEkxrTl0ui5H12/m4c9vMTo8IQQM0Vrvze+kUioceBVoC1wAlgPjgKmlEp0QXiAxMZEBAwaQlpbGwoUL6dq1a865fv360a9fPyPCMqRuUErx9ddfU7FiRVJTU40svxAeTVpUgFo9MvDxScfPfIwO++fjZz6Gj086g5670ejQhBCFMxBYobU+r60faUYBwwyOSQiPobVm3Lhx/Pnnn7z99tvXJCkuzml1Q2BgIACpqamOuJ0QIg+SqACBwZl0uftGyiT8yqMTJ1Am4Ve69r2RylUDjA5NCGH1lVLqgFJqplKqah7n6wAn7fZPZB8TQjjA1KlTmT9/Pv3792fy5MlGh2PPsLrB398fpZQkKkI4kXT9Aqo89zzJvj7Mf+BhztSsw/zIG7nhP/ey16Jpu3eP0eEJ4e1u1lqfyu5X/howB7ijuDdTSk0Gct5pBQcHEx0dnee1ZrM533OeTsouZbf5448/ePrppwkNDWX48OGsW7fOoOj+wdC6Yd26dfj7+3Pq1CmP/P/iDb8H3lBGcO9yOjVRUUqVAeYDzYBUIBZ4SGt9RCn1JfCv7ONJwESt9S/Zr5sN3ArEZd9qvdb66exzQcBMoD1gAZ7TWi8uSZx3fTgLbfq7cenXiHb0f38mymLhXEluLIQoMa31qeyvGUqpj4C8VmE9BTSw26+XfSyv+30AfGDbDwsL0717987z2dHR0eR3ztNJ2aXsABcvXmTs2LH4+vqyevVqWrVqZWB013KFuqFs2bIEBwd75P8Xb/g98IYygnuXszS6fk0HGmutW2EdxPa/7OPfAc2yj78JLMr1une11q2zt6ftjj8FpGmtGwK9gc+VUpVLEmCDI69R++xRMv88SNI3/8NyLoY6Z4/Q8OjrJbmtEKKElFLBSqkKdoeGAb/mcekS4G6lVA2llAImYP2QRAhRTFlZWdx3332cPn2aL774wqWSFFepG4KCgqTrlxBO5NRERWtt1lqv1n/P17cD66cZaK1XaK0z7Y7XUkoVpoVnCNbBcGitjwM/Av1LEuetZ1tSNzmJ9L8OkjxrKpnnzlAnOZmeMS1KclshRMlVBzYqpfYrpQ4A3YARAEqp/yml7gbQWh8DXgK2AkewtsZOMyZkITzDq6++yrp163jwwQcZNWqU0eHk5hJ1Q2BgoCQqQjhRaY9ReQJrq0pex1fbJS4ATyilHsDaRPu83fSDDh8Y91uVPRyt3gPTUettdVoaRxs2wZS2pCS3FUKUUPabjDb5nHsw1/4MYEZpxCWEp4uOjuaVV16hTZs2fPrpp0aH8w+uUjcEBgaSnJzsjFsLISjFREUp9RzQEIjMdXw4MBi42e7wf4BzWmuLUqo/sEYp1UhrnVSE5xV6UFykTy0uJSdQPvY0m4GWR37Dp144vfzruO3go8Jy5wFWjuDN5ffmsgsh8nfq1Cnuu+8+ypcvz+LFiylTpozRIbmswMBA4uLirn+hEKJYSiVRUUo9BdwD9NRap9gdH4K1STZSa33BdlxrHWP3/XdKqbeAxsBurC0sdSFnnHs94B9TkBRlUNyadU9w97ba/Gy2fipSMT6W7ocOEPav/6P3rZ8Vq8zuwp0HWDmCN5ffm8suhMhbRkYGgwcP5tKlSyxfvpz69esbHZJLk65fQjiX0wfTZ7dsDANu1VpfsTs+GOt0gj1tM3fYnQuz+74TUBlr31KwDrqfkH0uHOgOLCtJjA//37uAxlzhKgBXfS4Cmoc3vFOS2wohhBBuZcaMGezcuZMpU6Zw9913Gx2Oy5NERQjncvb0xGHA+8AxrIPewDpjV0dgLnAeWJ59HKwtK5eA2Uqp6kAW1umLB2mtr2Zf8y4wSyl1NPv8o1rriyWJM71jFXaeW0b8BWt+dCk1g53VlpFWo0dJbiuEEEK4jfnz57NixQq6d+/Oa6+9ZnQ4biEoKAiz2YzWGrv3MkIIB3FqoqK1PgPk+ZurtfYr4HU9CziXjHXmL4d5M+UknwY8w5n0WGAZpxIbERDQl7VNSjTrsRBCCOEWfv/9dx588EEqVarEt99+i6+vrAddGIGBgYB1zJ/teyGE45TGOiour3qVCljOa5RvAAA6Mw3LeU3zpjcYHJkQQgjhXElJSQwcOBCz2cyzzz5LjRo1jA7JbdiSk5SUlOtcKYQoDvnIBNAaUn2CsbSyJipZIYpUn2CDoxJCCCGcS2vNuHHj+OOPP3jnnXdo0ULWDysKW6Ii41SEcA5pUQFUFrQNy+T1LH8AOqSYaRuWiSU9y+DIhBBCCOf5/PPP+fbbb+nXrx9PPfWU0eG4HUlUhHAuaVEBmv0UxFsEk56VAUDIpVTe+q0S5w5uJ+zNrgZHJ4QQQjjezp07mTRpEg0aNODLL7+UweDFEBQUBEiiIoSzSIsKcLSDmYAGIfj5+uKjfEjNTCOgYQWqT2xrdGhCCCGEw128eJFBgwbh4+PDkiVLqFChgtEhuSVpURHCuSRRAdKCNYERVVBaEegXQGq6mcCIyvhVl3EqQgghPIvFYmH48OGcPn2azz//nFatWhkdktuSwfRCOJckKtmSd18ABUHBQaRmppG864LRIQkhhBAO99prrxEdHc2YMWMYPXq00eG4NWlREcK5JFEB0OBbqQyVhzcjuEI5Mspb97XWRkcmhBBCOMy6det4+eWXad26NZ9++qnR4bg9SVSEcC4ZTA+goPJ9TYHsVWaz0nP2hRBCCE9w+vRp7r33XsqXL8/ixYtlgUIHkMH0QjiXJCq5BAUFceGCdPsSQgjhOdLT0xk0aBCXLl1i2bJlNGjQwOiQPIK0qAjhXNL1K5egoCAZFCeEEMKjPPXUU+zcuZNnnnmGvn37Gh2Ox5BERQjnkkQlF0lUhBBCeJIFCxbw6aef0q1bN15//XWjw/EoMuuXEM4liUoutkRFBtILIYRwd3/++ScPPvggNWrUYP78+fj6So9vR5IWFSGcS2qsXGwD49LS0ihTpozB0QghhBDFk5SUxIABA0hNTWXVqlXUqFHD6JA8jgymF8K5pEUlF1ulI824Qggh3JXWmvHjx/P777/zxhtv0K1bN6ND8kjSoiKEc0mikoskKkIIIdzdF198wbx58+jbty9PP/200eF4LElUhHAuSVRykURFCCGEO/v555+ZOHEiDRo0YPbs2SiljA7JY8lgeiGcSxKVXCRREUII4a4uXbrEoEGDMJlMLF68mAoVKhgdkkfz8/PDZDJJi4oQTiKD6XORREUIIYQ7slgs3H///Zw6dYr//e9/tG7d2uiQPJ5SiqCgIElUhHASaVHJRRIVIYQQ7uj1119nzZo1jB49mjFjxhgdjtcIDAyUREUIJ5FEJRdJVIQQQrib9evX89JLL9GqVSumTp1qdDheRRIVIZxHEpVcJFERQgjhTk6fPs29995LuXLlWLx4cc4Ab1E6AgMD5T2DEE4iY1RykRk8hBBCuIv09HSGDBnCxYsX+e6772jYsKHRIXmdwMBAEhISjA5DCI8kLSq5SIuKEEIId/HMM8+wfft2nn76afr162d0OF5JBtML4TySqOQiiYoQQgh3sHDhQj7++GO6du3KG2+8YXQ4XkvGqAjhPJKo5CKJihBCCFf3559/MmbMGKpXr86CBQvw9ZWe3EaRREUI55GaLRdJVIQQQriy5ORkBg4cSEpKCitXrqRmzZpGh+TVAgMDSUtLIysrCx8fH6PDEcKjSItKLpKoCCGEcFVaa8aPH8/Bgwd544036N69u9EheT3bJDxms9ngSITwPE5NVJRSZZRSy5RSh5RS+5RS65VSDbPPVVNKrVVKHVZK/aaUutnudcU65wiSqAghhHBV06ZNY+7cudx9990888wzRocj+Pt9g3T/EsLxSqNFZTrQWGvdClgO/C/7+FvADq11I2A0ME8p5VfCcyUm0xMLIYRwRbt27eKJJ54gPDycOXPmoJQyOiTB3+8bJFERwvGcmqhorc1a69Vaa519aAdQL/v7wUBU9nW/AGeBbiU8V2K+vr74+/tLoiKEEMJlXLp0iYEDB6KUYsmSJVSoUMHokEQ2SVSEcJ7SHkz/BLBcKVUZ8NNan7c7dwKoU9xzuR+klJoMTLbtBwcHEx0dnWdQZrP5mnP+/v6cPn063+s9Se6yextvLr83l10Id2KxWLj//vs5efIkM2bMoE2bNkaHJOxITwwhnKfUEhWl1HNAQyASCHT287TWHwAf2PbDwsJ0796987w2Ojoa+3Ply5cnKCiI/K73JLnL7m28ufzeXHYh3Mkbb7zBmjVrGDVqFGPGjDE6HJGLtKgI4TylMuuXUuop4B7gdq11itb6EpCplKphd1k94FRxzzky3qCgIPlkRAghhOF++OEHXnzxRVq2bMnUqVNlXIoLkkRFCOdxeqKS3QVrGHCr1vqK3alFwITsa9oDtYBNJTznEJKoCCGEMFpMTAz33nsv5cqVY/HixTmzSwnXIrN+CeE8Tu36pZQKA94HjgEbsz8JStNadwSmAF8rpQ4D6cBwrXVG9kuLe84hgoKCiI+Pd+QthRBCiELLyMhg8ODBxMXFsWTJEho1amR0SCIf0qIihPM4NVHRWp8B8myn1lpfAHo58pyjSIuKEEIIIz3zzDNs27aNJ598knvuucfocEQBZDC9EM4jK9PnQRIVIYQQRlm0aBEfffQRXbp04c033zQ6HHEd0qIihPNIopKHoKAgzGYzFovF6FCEEEJ4kb/++osHHniA6tWrs2DBAvz8HLaesXASSVSEcB5JVPIgA+OEEEKUtuTkZAYMGEBKSgrz588nNDTU6JBEIch7BiGcRxKVPNgqHen+JYQQojRorXnooYc4ePAgr7/+Ot27dzc6JFFI0qIihPNIopIHSVSEEEKUpunTp/P1119z11138cwzzxgdjigCGUwvhPNIopIHSVSEEEKUll27dvH4448THh7OV199hckkf5rdibSoCOE8UhvmQT4dEUIIURri4+MZOHAgSikWL15MxYoVjQ5JFJEkKkI4j1PXUXFX0qIihBDC2SwWCyNGjODkyZNMnz6dtm3bGh2SKAZJVIRwHmlRyYMkKkIIIZztrbfe4vvvv2fkyJE8+OCDRocjisnPzw9fX19JVIRwAklU8iCJihBCCGfasGEDL7zwAi1atODzzz9HKWV0SKIEAgMD5T2DEE4giUoeJFERQgjhLDExMQwbNozg4GAWL16c8zdHuK/AwEBpURHCCWSMSh4kURFCCOEMGRkZDBkyhLi4OJYsWcINN9xgdEjCASRREcI5pEUlD5KoCCGEcIYpU6awdetWJk+ezD333GN0OMJBgoKCJFERwgkkUcmDLVFJirvKhal7ybwolY8QRlNKjVZKaaVUvzzO1VNKZSml9tptDQwIU4h8LVmyhA8//JAuXbrw1ltvGR2Ox3CFukFaVIRwDun6lYdAvzIAXD54jozMRK6uP0nF/g1RPiaUn+R2QpQ2pVQ9YCywo4DLErXWrUslICGK6NChQ4wePZpq1aqxYMEC/Pz8jA7JI7hK3RAYGEhMTIwzHyGEV5J33XlI+OJ3ABJjLgOQui+Osy9vJ+bFrUaGJYRXUkqZgP8BjwFpBocjRJGlpKQwcOBAkpOT+fbbbwkNDTU6JI/gSnWDtKgI4RySqOQhdEIbAFIzzQBoNH7h5aj8SAsjwxLCW00Gtmqtd1/numCl1C9KqT1KqReVUj6lEZwQBdFa89BDD3HgwAFeffVVbrnlFqND8iQuUzdIoiKEc0jXrzzMfeUxAFIz0siyZGJSPuz4ZQlHf9zH5G9XGBydEN5DKdUcGADcfJ1LzwG1tNaxSqlKwALgSeCdPO45GesbHACCg4OJjo7O86Zmsznfc55Oyu6Ysq9evZqvvvqKjh070qZNG5f/mbrLv7ur1Q0JCQlkZGSwevVqfHw84zMSd/m/UBLeUEZw83Jqrb1iq1Wrls7P2rVrr9m/cOKYBvSdjbvrdePf0aee+VH/8exKHXfqRL73cFe5y+5tvLn8zio7cEY76PcWeAjrG40T2ZsZiAUeus7rhgErC/OMotQN3kTKXnK7du3S/v7+ul69ejo+Pt4h93Q2Z/67e3LdMGrUKA3oxMTEkv2QXIg31AHeUEatXbuc16sXpOtXHqrWqUcZvwDOJp/kYOJWtsYtw7dyIJXD6hgdmhBeRWv9hda6pta6nta6HtYBs+O01l/YX6eUqqaU8sv+PgC4B/i11AMWItvly5cZOHAgAIsXL6ZixYoGR+RZXK1uCAwMBGRZAyEcTRKVPCil8PVVJKYn0PaOvpxNOcKOiytRShkdmhAim1LqFaXUhOzdLsCvSql9wB7gPPC6YcEJr2axWBgxYgQnTpzg008/pV27dkaH5FWMqBtsiYqMUxHCsWSMSh601pQJCCCgSmWiKqznsYfGcnHP72itJVkRwkBa6+52379o9/1SYKkRMQmR29tvv82qVau4//77GTt2rNHheAWj6wZJVIRwDmlRyUO6JZ2K1apxwRzPgYsHWMpP9HjkUdIt6UaHJoQQwoVt3LiR559/nubNmxMVFSUfbnkJ20LRkqgI4VjSopKHDnM7cDrtNJZUCzWowdoTa1l7Yi0+yoe9I/YaHZ4QQggXFBMTw9ChQwkODmbJkiU5b16F55MWFSGcQ1pU8rCkzxJCgkOwpFsA8FE+dKrZicV9FhscmRBCCFeUkZHB0KFDiY2NZdasWdxwww1GhyRKkQymF8I5JFHJQ8OKDalaoSqWNAt+Jj8s2sItdW6hYcWGRocmhBDCBT377LNs2bKFSZMm5cz2JbyHtKgI4RySqOQjPisena4Z3mw4CsXyI8uNDkkIIYQLWrp0Ke+//z433XQTb7/9ttHhCANIoiKEczg9UVFKfaKUOqGU0kqp1tnHKiul9tpth5RSmdmrxqKU+lEpddzu/CS7+1VTSq1VSh1WSv2mlLreqrRFprWmXHA5dJbmsZaP8UGPD6hVtpZtsSghhBACgMOHDzN69GiqVq3KwoUL8fPzMzokYQAZTC+Ec5TGYPrFwDvAFtsBrfUloLVtXyn1FNBNax1v97pJWutledzvLWCH1vo2pVR74DulVLjWOsNRASul+Fe9f/EXf5GamkpknUgi60Q66vZCCCE8QEpKCgMGDCApKYmlS5dSq1Yto0MSBpEWFSGcw+ktKlrrzVrrM9e5bAwws5C3HAxEZd/7F+As0K34EebN9ulItwWb2b3rkKNvL4QQwo1prXn44Yc5cOAAr7zyCpGR8mGWN5PB9EI4h+FjVJRSNwEVgVW5Tr2llDqglFqglKqffW1lwE9rfd7uuhNAHUfGlJCQDNr6ozkdqHnp4CFiL1y2HhdCCOH1Zs6cyZw5c7jzzjt59tlnjQ5HGExaVIRwDldYR2UM8JXWOtPu2P1a69PKulLWI1iTmGZFualSajIw2bYfHBxMdHR0nteazeZrzo3yq0ZiUvbijulp7ApvSMvfT6IsmtmZsUUJw+XlLru38ebye3PZhSiJPXv28Oijj1K3bl2++uorTCbDP/MTBpNERQjnMDRRUUqVxdqVq739ca316eyvGvhMKfWeUqqy1vpS9qD7GnatKvWAU7nvrbX+APjAth8WFqZ79+6dZxzR0dHYnxv9yEi+SbOQDGizGZPOovbZY0RuX0XvqXNKUmSXk7vs3saby+/NZReiuC5fvszAgQPRWrN48WIqVapkdEjCBUiiIoRzGP0x0BBgn9b6T9sBpZSvUqq63f4A4EL2AHyARcCE7HPtgVrAJkcGFTlgK42rWnMfkzkJjaJ76I9EDtjmyMcIIYRwIxaLhZEjR3L8+HE+/vhjbrzxRqNDEi5CZv0SwjlKY3riaUqpM0AYEK2UOmJ3Oq9B9AHA99njU/YBDwN3252fAtyklDoMzAaGO3LGL4Ab267imF8LANqmbQZgU0ZPbmy70pGPEUII4UbeeecdVq5cyfDhwxk/frzR4QgXIoPphXAOp3f90lrnW5trrW/K41gykO/HVFrrC0Avx0SXt2+ePUtQgAWAjqkb6aE1G67ezjfPneWRz5s489FCCCFc0I8//sh//vMfIiIiiIqKwjqEUggr6folhHO4wmB6lzPkuRu58NkYfgPiTzbjZn6hOX/Q6dnNRocmhBCilJ07d46hQ4cSHBzMkiVLCA4ONjok4WJ8fHzw8/OTREUIBzN6jIpLqhxWjkrVrC0nF0+0IGbrQ1So0o7KYeUMjkwIIURpysjIYMiQIVy4cIGZM2fSuHFjo0MSLiowMFASFSEcTBKVPCil0GceACCgSjrJ59pwdscEaeoXQggv89xzz/HTTz/xxBNPMGjQIKPDES4sKChIEhUhHEwSlTxorQmvb11DUp/9jtvurUr5KmWwzpYshBDCG3z33Xe89957dO7cmXfeecfocISLCwwMlMH0QjiYJCp5UJlp9BvblEplAzh1Lo76Fz7kthF1UZlpRocmhBCiFBw+fJhRo0ZRpUoVFi5ciL+/v9EhCRcnXb+EcDwZTJ+XN2qCtlCnbAanrmr4bYl1Uz7wUrzR0QnhlpRS7YBmWuuvlVIVgTJa63NGxyVEbmlpaQwcOJDExETWrVtHWFiY0SF5NE+pGwIDA7l48aLRYQjhUaRFJS8PbYPwbgRVrMCJRPgxsCWEd4OHthodmRBuSSn1MDALeDn7UCVgnmEBCVGAqVOnsn//fl555RV69uxpdDgezZPqBmlREcLxJFHJw5WydThb73aOh1ZBZ2n+HViNs+F3cKVsHaNDE8JdjQM6AQkAWuujQFVDIxIiDzNnzmTdunXcfvvtPPfcc0aH4w08pm6QwfRCOJ50/cpDl0WdUWh8K1grnKTzf9H7+CH08S/4beQ+g6MTwi2laa1Tc82cl2lUMELk5ddff+WRRx6hWrVqfP3115hM8lleKfCYukEG0wvheFIL5+HttjMoq5pS378+AFVPV6OcasbbbacbHJkQbitOKXUDoAGUUqOAU4ZGJISdK1euMHDgQLTW/Oc//6Fy5cpGh+QtPKZuCAwMJCsri4yMDKNDEcJjSItKHnYtXUtv3YxTZcqyiU34x/rR7VhTdh1fy50tOhodnhDuaCLwLdBEKXUaazePuwyNSIhsWmtGjhzJsWPH+Pzzz6lfv77RIXmTiXhI3RAYGAhAamoqfn5+BkcjhGeQFpU8jBs3htjABMqHWFeiv3r1CrGBCYwbN8bgyIRwT1rrI0BHoB3QG2iptT5ubFRCWL377rusWLGC++67jwkTJhgdjlfxpLrBPlERQjiGJCp5+G1fD64GXaZ82fIopbh6NYErQZf5bd8tRocmhDsrC/gBZYBWSqm2BscjBJs2beLZZ58lIiKCadOmkWushCgdHlE3BAUFAZKoCOFIkqjkIWyyHx0OlcNkMlEpMJCEhAQ6HipH2GRpyhWiOJRSk4AYYBmwIHubb2RMruDIkSN06tSJI0eOGB2KVzp37hxDhgwhKCiIxYsXExwcbHRIXseT6gZbi4onDaiXOkoYTRKVPIQuXkpZPwv/2rqFulqTEhdHWT8LoYuXGB2aEO7qMaCx1rq+1rpR9naD0UEZxWw2k5CQwJQpU9i5cycvvPACCQkJmM1mo0PzGpmZmQwdOpQLFy4wc+ZMmjRpYnRI3spj6gZP6vplNptJTk7miSeekDpKGEoSlTy0+foo1X1+ptaZGKr7+XElI4MKegdtvz5mdGhCuKsYrfVZo4NwFcHBwYSEhLB06VIA5s+fT0hICGXLljU4Mu/x3HPPsXnzZh5//HEGDx5sdDjezGPqBk9KVIKDgxkwYACrV68GpI4SxpFEJQ9rnuhKtT2ZaBQJVa2zv+idKax+oovBkQnhtt5QSn2ilOqllLrZthkdlFH2799PZGRkzr6Pjw+RkZHs2yfrNJWGZcuW8e6779K5c2feffddo8Pxdh5TN3hSorJ//35q1KiRsy91lDCKTE+chzqVgvilRjs+CGvIH+dOwNG9nAkM585KQUaHJoS76gyMALoAWdnHNNDBsIgMFBERQb9+/diwYQMAFouF/v37ExERYXBknu/o0aOMGjWKKlWqsHDhQvz9/Y0Oydt5TN3gSYPpIyIi0FoDYDKZpI4ShilUi4pSqoKT43ApES9FM6J2P3bUbI5PuaoAvFbmBpq/tM7gyIRwWyOAelrrtlrr9tmb270RcaQ5c+Zcsz979mxjAvEiqampDBgwgISEBObNm0dYWJjRIQkPqhs8aTB9fHw8Fy5cAKBatWoopaSOEoYobNevw0qpGUqpVk6NxkWseaIrlYP9MSnwLW9NVPxT46XrlxDFd1JrfcXoIFyF1prw8HAqVKgAWD+9DA8Pz/kEUzjHo48+yr59+3j55Ze59dZbjQ5HWHlM3eBJXb+io6NzvjeZTCxZskTqKGGIwiYqDYGDwEKl1E9KqcFKKR8nxmWoxjXK83hkI7SGwIrVAbihbBqNa5Q3ODIh3NYvSqmF2XXH3bbN6KCMopRi4cKFZGZmAnD69GkWLFgga3g40axZs5g1axa33XYbzz//vNHhiL95TN3gSYnKmjVrAGjfvj3nz5+nT58+LFy4UOooUeoKlahora9qrT/SWjcGXgfeA04ppf6jlPLIiecX7z6DUjCmZ3OUXwC/H5YZv4QogXZAVeAhYFL2NtHIgIymtSY5ORmAq1evcuyY1DHOsnfvXh555BHq1KnDN998g8kk88i4EI+pGzwlUbFYLKxdu5Z69erRoUMHLBYLsbGxRoclvFShB9MrpcoDD2CtTA4CM4AeQDTWQXAeQ2tNnUqBPHZLQ3pF1CAqrDaJCXForeXTBCGKQWvdw+gYXI3ZbEZrjb+/P+np6ezZs4cGDRoYHZbHuXLlCgMHDiQrK4tFixZRuXJlo0MSdjypbvCUwfR79uwhLi6OHj16EBoaCsDZs2epWbOmwZEJb1TYwfTTgENYu4DdrbW+XWu9VGv9GOBxtb5Siqn3taNXhHVqvqaN6pNxNc7gqIRwP0qpRtlfW+a1GR2fkZKSkgDo2LEjALt37zYyHI+ktWb06NEcPXqUjz76iA4d3HKMtkfyxLrBUwbT27p93XjjjdckKkIYobDt34ewrhz7qNb6r1znbnFwTC6nTp06JCcnc/nyZaNDEcLdfJj9dXke2zKDYiqyI0eO0KlTJ44cOeKwe9q6fbVr147AwEBJVJzgvffeY9myZQwbNoyHHnrI6HDEtTyibrDnKV2/Vq9eTfny5WnWrJkkKqLYHPV3s7BjVN4HQpRSw5VS9yml6tidO1eiCNxA7dq1ATh16pTBkQjhXrTWd2V/Dc9jq290fNeTnp5OQkICL7zwAjt37uSFF14gISEBs9lc4nvbWlTKly9P69at2bNnj8yo40CbN2/m2WefpWnTpkyfPl267boYd68b8uIJicqlS5fYuXMnPXv2xNfXl1q1agEQExNjcGTCXZjNZof+3Sxs1697gV+BAcBAYLdSamixnuiG6tSx5mWSqAhRPEqpZYU55mr69etHSEgI8+fPB2D+/PmEhIRQtmzZEt/b1qJStmxZ2rZtS3x8PCdPnizxfQWcO3eOIUOGUKZMGZYsWeKQfy/hHO5aN+TFExKVdevWobXm9ttvB5AWFVFkwcHBDv27WdiuXy8CN2qt+2ut+wPtgZeL9UQ3ZEtUTp8+bXAkQritOnkcc/mR459//jmRkZE5+z4+PkRGRrJv374S39uWqAQHB9OuXTvAOohVlExmZibDhg3j/PnzzJw5k6ZNmxodkiiYW9YNefGERMU2PsWWqFSoUIEyZcpIoiIKbf/+/YSHh+fsl/TvZmETlRSt9XHbjtb6BFCo0WJKqU+UUieUUlop1dru+Aml1F9Kqb3Z2xC7c42UUtuUUoeUUr8opSIKc85ZpOuXEMWjlBqvlPoVaKyU2mO3HQWOX+/1RqtXrx79+vXL2bdYLPTv35+IiJJXO7auX7YWFZAB9Y7w/PPPs2nTJh577DGGDBly/RcIQ7h73ZAXk8lEQECA2w6mt01L3LJly5wuX0opQkNDJVERhda4cWPi4+MB8Pf3L/HfzcImKt8rpV5WSoUppWorpV4EViqlymdPW1yQxVinL86rT8MQrXXr7G2B3fFpwHSt9Q3A28DsQp5zirCwMEASFSGKYS3WdRFO8fcaCZOwdiHtb2BchTZnzhzAmlAopZg9e7ZD7mvfotKsWTMCAgKkRaWEVqxYwdtvv03Hjh157733jA5HFMzt64a8BAYGum2Lyu7du4mLi8tpTbGRREUUxapVq7h69SpKKSZOnFjiv5uFXUflP9lfX8x1/AVAA/muUq+13gwUeiCjUqoacCPQK/vQEuAzpVRDICG/c1prx03Hk0tgYCDVqlWTrl9CFJHW+iTWDyncsv+N1prw8HCSkpK4fPkyX3/9NfPmzXPImkq2FpXg4GD8/Pxo2bIlu3fvlvWaiuno0aOMGDGCypUrs3DhQvz9/Y0OSRTA3euG/LhzopK725dNaGgoW7ZsIT09XX6vxHV98cUXAMycOZPRo0fTuXPnEv3dLOysX6YCtnyTlEL4Sil1QCk1UylVNftYbeCc1joz+9ka6ycuda5zzqlq164tLSpCeBmlFAsXLqRmzZqkp6fTr18/Fi5c6JBEwn4wPVinKY6Li5PZdYohNTWVgQMHkpCQwLx583LGFQpR2tw9USlfvjw33XTTNcdtA+rPnz9vRFjCjRw7dox169Zxzz33MHr0aIAS/90sysr0tYGu2bubtNYl/Wt6s9b6lFLKD3gNmAPcUcJ75lBKTQYm2/aDg4OJjo7O81qz2ZzvOZuAgABiYmJYvXo1Pj4lyc1cS2HK7sm8ufzeXPai8vPzIyMjw6H3tO/6BeQMqN+9e3dOd1NROI899hh79+7l5ZdfplevXtd/gRBO4q6JysWLF9m5cyf33HMPfn5+15yzn/lLPgQQBZkxYwYAEyZMcNg9C5WoKKX6AjOBLVi7en2klBqjtV5Z3AdrrU9lf81QSn2EdVFJgNNATaWUr9Y6U1lTsDpYW04SCjiX+/4fAB/Y9sPCwnTv3r3zjCU6Opr8ztmsWbOGbdu20aJFi5zB9Z6gMGX3ZN5cfm8ue1E5I1GxH0wPXDOgvm/fvg59lif78ssvmTlzJr179+aFF14wOhzh5YKCgtyy5SH3tMT2ZIpiURjp6enMnDmTBg0aXDNbZkkVtkXlJaCTbRxI9niRhUCxEhWlVDDgp7W+kn1oGNZ1WtBaxyql9gDDsQ6UHwCcsXt2vuecyX7mr5qBVbi04C8qD2mMb5VAZz9aCI+glKoJhGNX79jGsLm60mhRad68OX5+fjKgvgj27dvHww8/TO3atfnmm28wmQo7P4xwJe5cN+Tmri0qtvEpt9122z/O2WYAk0RFFOS7774jLi6Op59+2qF1cWHv5GOfDGR/X9jFIqcppc4AYUC0UuoIUB3YqJTar5Q6AHQDRti9bDwwXil1CPg3MLqQ55zG1tx54shxzq3+nfTTCZxb/TsWcyY6w1IaIQjhtpRS/wH+AD4D3s/e3GZaJj8/PywWCxaL437X7QfTg3UaxxYtWjh8iuIjR47QqVMnjhxx+uc5perq1asMHDiQrKwsFi1aRJUqVYwOSRSDu9cNuTkyUSmt313btMStWrXKSUrs2VpUZPycKMi0adPw9/dn1KhRDr1vYVtUYpVSDwKzsvdHA3GFeaHWenw+p9oU8Jq/gM5FPedMtkTl4Jdb6NapLqDg9xTOvrwdFIS92bXgGwjh3R4AGmitLxkdSHHY+mxnZGQQEBDgkHvmblEB6ziVPXv2cO7cOWrWrFmi+5vNZtLT03nhhRfYuXMnL7zwQs4fkjJlypTo3kbTWjNq1CiOHDnCp59+SseOHY0OSRSfW9cNudkSlZLM3lfav7u7du3i4sWLPPjgg3met9VF0qIi8vPnn3+yceNGhg0bRtWqVa//giIobIvKBOBBIDV7exBry4bXsHX92pf6F6mBZznZ4RVSA8/ya9AfPBz+usHRCeHyLrjzGxH7RMVRkpOTCQgIwNf378+LHLnwY3BwMCEhIcyfPx+A+fPnExISkjMmxp29//77LFu2jKFDh/LII48YHY4oGbeuG3ILDAzEYrGUqK4o7d/d/KYltilXrhzlypWTREXka/r06YBjB9HbXDdRUUr5AHdqrTsBlYHKWuvOWutjDo/GhdWoUQM/P1/OpJ0nvuFyzCHHiG+4nPO1L/L+gA+ufwMhvNt6pdRHSqkOSqmWts3ooArLtnZAenq6w+6ZlJT0jzcetpm/HDFOZf/+/f8Y0NitWzf27dtX4nsb5ciRI0RERDBlyhSaNm3KjBkzZM0Z9+fWdUNumZmZAPz222/Fvoftd9f2f1spRWRkpNN+d23TEnfunH9nFVn0UeQnNTWV2bNn07RpU7p2dXzvousmKlrrLGBk9vdJWuskh0fhBkwmE5WrQHLSMZJq7gQFSTV3Ur/uTE7t62d0eEK4uhFAX2ABsDx7W2ZkQEXhrBYV+25fAC1atMDX19chLSoRERH/mKo3IyODiIiIEt+7tJnNZhISEnjqqaf4/fffUUoxZ86ca1qjhNty67rBxvZ/dO/evQC88sorJCQkYDabi3yviIgIIiMjsS4VZ+3qePfddzvld/fixYv8/PPP9OrV6x/TEtuTREXkZ/HixVy+fJnx48c75YOjwnb9Wq+Uus/hT3czDeq3JvaCBYu17sCiISuhHh06FHuWZiG8gtY6PI+tvtFxFZYzEpWkpKR/JCplypQhIiLCYQPqo6KiAHISlm3btuV083Antq4wy5cvByArK4sOHTp4RDc2b+fudYON7f/o4cOHAVi+fHmJumvZVve29fefOnWqYwLNJTo6Ot9pie2FhoZy5coVUlJSnBKHcF9RUVGUKVOGESNGXP/iYihsojIe+FopZVZKxSulLiul4p0SkQsLD29CYlImaWYLSvljUvBnSHnKBt9gdGhCuLzsrh3PZW83Gh1PUTirRSWvNzHt2rUjJiaGCxculOj+WuucbiizZ89m+vTp+Pr68sADDxAf717V9/79+wkPD8/Z9/HxcWpXGFG63LlusLF117JNy2oymYr9f1RrnTOGbe3atYC17rG1sDhSQdMS27PN/HXu3DmHxyDc1/79+9m2bRtDhw6lYsWKTnlGYROV1ljnOG+Mdbau1hQwa5enss38FRtroU7tUYCJTmWLP7OHEN5CKTUOWAxUA6oCS7JnEizKPUYrpbRSql8+5+9SSv2plDqslFqqlCpf4sCzlVbXL/h7QH1Jx6lkZmaSkJBA+/btqVmzJmPHjmX69OmcP3+exx9/vET3Lm3nz5/n+PHjgHW8kMVioX///m7ZjU1cy93rBpuIiAj69euXk0yU5P/osWPHiI+PZ9CgQbRt25bmzZuTkJBAVlaWQ2O2WCxER0fTqlWrnEQkP7Loo8jLtGnTABg/3nnzaxU2Ublfa33SfgPud1pULiosLAyAssETaNhwCi1bTCWwTJhTPuUQwsM8CrTTWk/UWk8CbgQK/W5ZKVUPGAvsyOd8WWAm0E9r3Qg4CzhsmXJndf3Kr0UFSp6o/PTTT1y9epU+ffrkHBs1ahR33XUXc+fOZenSpSW6f2k5f/489913HyaTCZPJxMSJE1FKMXv2bKNDE47h1nWDvTlz5lzzwWVx/4/OmzcPgPvus/a4HzRoEJcuXWLjxo0ljtGebVriO+6447rXyqKPIrekpCS+/vprWrVq5dQp4gubqNxTyGMerW7dugAkHLkEr4dSNcGPFi0+kxYVIQpBax2X1/fXo5QyAf8DHgPS8rnsduBXrfWf2fufA8OKGeo/ODpRycjIID09Pc8WlZYtW2IymUo8TmXlSuvYubvvvjvnmFKK6dOnU6lSJSZMmEBsbGyJnuFsWVlZDB8+nAsXLtChQweWLl3K22+/zZIlSwgPD5cPiTyEO9cNNlprwsPDWbJkCZUqVaJx48bF+j+qtWbu3LlUq1aNnj17AtZEBWDRokUOjXn16tVA/tMS25MWFZHb/PnzSUxMZMKECU59H1zglClKqd7AbUAtpZT9HLwhTovIhdWpVgGAUz/OgZvSYcWj8MA68A+GwAqGxiaEizuslHodmJa9PxY4XMjXTga2aq13F1AZ1gFO2u2fAGoqpXy11pn2FyqlJmffE7AOgo2Ojs7zpmazmejoaE6cOAHA5s2bSzx2BP5elf7y5ct5PrtOnTps3bo137iuR2vNggULqFq1KufOneP8+fPXnB8/fjxvvvkm99xzD88//3yef2RsZTfSvHnz2LBhA3fccUdOd7Xo6GgCAwMZM2YM69atc8pzXaHsRjGg7G5dN9gbM2YMYG19OHr0KB999FGR/48ePnyYv/76i7vvvpsNGzbkHK9bty4LFy6kX79++Pj4FOme+Zk/fz7BwcEkJibmWU77MtrqkK1bt9KsWTOHPN8VeMvvujPK+c477xAYGEj16tWd+zPUWue7Ad2Al4Bz2V9t2ySgXkGvdbWtVq1aOj9r167N95y9K1PKaUAPbRWkp780Tl98qbbWL5W3bm6qsGX3VN5cfmeVHTij/1mXVAXmAxezt2+Barmvy+N1zYHtgF/2/o9Yu3Dkvu5JYJrdfhCQBfhe7xmFqRs++OADDegtW7Y45Gd05swZDejHH388z/MjRozQgL548WKx7n/w4EEN6IcffjjfawYNGqQBPXfu3DzPG/278eOPP2qTyaRbtmypU1JSSvXZRpfdSM4suyfWDXl59NFHNaBPnDhR5J/RpEmTNKB37NhxzfH//ve/GtDr168v8j3zEhsbq5VSeuDAgfleY1/G1NRUDehhw4Y55Pmuwlt+1x1dzl9++UUDety4cSW+V171gv1WYNcvrfUmrfV/gdu01v+12z7UWp8oRB7kUYLun0v5ABO/JZQlhhps5CbMAVXJGPSN0aEJ4dK01nFa66Fa6yrZ2zCtdWH6HXUF6mH91PUE0AmYrpR6KNd1p4C6dvv1gHM61yemxeXorl/JyckAeXb9gpKPU8mr21dun3/+OdWqVeORRx5xue4ccXFx3HvvvQQGBrJw4UICAwONDkk4ibvXDXlp3rw5UPRFH7Oyspg/fz4NGjSgQ4cO15xzdPevdevWFWpaYpsyZcpQqVIll6srhDFsg+idsRJ9boUdo3JQKTUke+rAF22bUyNzQW8s3kOZkCqcvZoJKH6jCW+lDeeNxSVfRVoIT6SU6pb99e68tuu9Xmv9hda6pta6nta6HtYBs+O01l/kunQt0FYp1SR7/2Gsn9I6hKMTFVvXr/zWWLDN/FXccSorVqygbNmydO/ePd9rqlSpwowZM7hy5QoPPvigy4z3sFgsjBgxgrNnzxIVFUXjxo2NDkk4gafUDXlp0aIFAAcOHCjS6zZu3Mi5c+e49957/9Eds2nTpkRERLB06dKcacdLwjY+5XrTEtuTRR8FwNWrV5k3bx4dOnSgTRvnTwBc2ERlPtbBalWAcnabVxk5fATVy/tzNSGBJa27crVMIHWJYeRw5yxyI4QHGJ79dVIe28SS3Fgp9YpSagKA1joReBBYppQ6AoQBr5bk/vZKu0WldevWKKWKlajExcWxfft2evXqRUBAQIHX3n333YwYMYI1a9Ywc+bMIj/LGd59913Wrl3LAw88wPDhw6//AuGuPKJuyIutRaWoicrcuXOBv2f7ym3QoEFcvHiRTZs2lSi+rKwsoqOjad269XWnJbZnS1Rc5UMNYYxvvvmGlJQUp05JbK/AwfR2WgBNtJf/75z/wiT8A8qRlXmGC1nwc71mhGw5xvwXJvHvecuMDk8Il6O1Hpv9tYeD7tfd7vsXc51bAaxwxHNyc1aLSn6JStmyZWncuHGxun6tXr0arXWB3b7sffzxx2zYsIFJkybRs2dP6tWrV+RnOsrWrVv5z3/+Q7Nmzfj0008Ni0M4n6fUDXkpX748devWLVLXr9TUVJYsWUK7du3ybUUcNGgQL7/8MgsXLiQyMrLY8e3atYtLly4xbty4Ir0uNDSU5ORkEhMTKV/e4UvRCDegtWbatGmEhIQwZMiQUnlmYVtUTgP+zgzEHXw85kX+6mStHLLiznO0em2iBjzEJ2O8rhecEEWilPq5MMdclb+/tfpzdItKfl2/wDpO5dixY1y+fLlI916xYgVKqUKtjQBQoUIFZs2aRVJSEg888AAWi6VIz3OUS5cuMXToUPz9/Vm4cCFBQUGGxCFKl7vXDflp3rw5f/zxR6HrjFWrVpGYmJhvawpAs2bNaNasWYm7f9lWoy/s+BQbmaJYbN++nQMHDjBixIh8P2hztMImKkeAH5VSzyilHrdtzgzMFUXXCKJq9qDOrAvnURYLtS+eZ20N+YMqxHVc03qrlPLDjbqP2lpU0tPTHXK/63X9gr8H1P/666+Fvm9aWhrr1q2jc+fOVK1atdCv69WrFxMmTGDjxo1MnTq10K9zFK01o0aN4syZM0ydOlVWnPcubl035KdFixZkZGRw6NChQl0/d+5cTCYTQ4cOLfA6R3T/Wr16NSEhIXTu3LlIr5NERURFRQHOXYk+t8ImKgHAn0BToI3d5lVaRjShvjkBAH3hLFopGh3eT8uIJtd5pRDeSSk1RSl1GWihlIq3bUACsNng8AqttAfTQ/EG1P/4448kJSUVutuXvXfffZfw8HCmTJnC4cOFXcbCMT788ENWrVrF/fffz6hRo0r12cIYnlI35Mc2oL4w3b/i4+NZvXo1t9xyCzVr1izw2pLO/hUXF8euXbvo1asXvr6F7f1vJavTe7f4+HgWLlxIly5dSvXDpEIlKlrr0Xltzg7O1WRlZRHTuhMAjU8eQWnNwUYtycrKMjgyIVxWFNYPNX7g2g85QrXWpfeRTAmV9mB6IGc2laKMU1mxwtoNv0+fPkWOqWzZsnz55ZeYzWZGjhxZavXazp07mTJlCo0bN+bzzz936grHwqV4RN2Qn6LM/LV48WIyMjIK7PZlExERUaLuX9HR0UWaltierUUlJiamyK8V7m/OnDmkpaWVypTE9gpMVJRSz9p9PyjXufedFZSrUkrROjgIpRR1Av2YqpNokZYqf1iFyIfW+qrW+oTW+nat9Um7rWgDLwxW2oPpwTogt1GjRoVuUdFas3LlSho0aEDTpk2LFVe3bt2YOHEi27dvZ+nSpcW6R1FcvnyZIUOG4Ovry8KFCwtsYRKexVPqhvw0btwYX1/fQiUqc+fOpUyZMtxzzz2FuvegQYOIi4tj8+aiNzzZxqcUZVpiG+n65b201kRFRVG5cmUGDBhQqs++XouKfXLybK5zDpmpw52YTCa+HTWY+vXrc+DAAe7peTPfDO+PyVTYHnRCeCelVDWl1BdKqW1KqT22zei4CstZLSrXe2Perl07Dh8+TEJCwnXvuX//fk6fPk2fPn1K9OHJ66+/TuPGjfnqq684ePBgse9zPVprHnjgAU6ePMnHH39My5YtnfYs4brcvW7Ij7+/P40bN75uonLq1Ck2b95Mnz59Cj2Tlq3718KFC4sUU1ZWFmvXrqVNmzbX7WKWl+rVq6OUkkTFC/34448cOnSIUaNGUaZMmVJ99vXeYat8vs9r32t0aNOeQ4cOceHQGaNDEcJdzAROYF2L6SXgLPC9kQEVhRFdv+DvcSqFGVBfkm5f9gIDA5kzZw5ZWVmMGDHCYWXO7bPPPmPZsmUMHTqUsWPHOuUZwi24dd1QkBYtWnD8+PGcFtS8zJ9vXXuyMN2+bCIiImjatGmRu3/98ssvxMfHF6vbF1jrwWrVqkmi4oVsK9EXdUprR7heoqLz+T6vfY+nMyxYzJk0D6gLwLqpS7CYM9EZxkznKYQbqa21fhtI01qvBO4BehocU6E5enriwgymh79n/ipM96+VK1cSEhJC165dSxxfx44dGTx4MHv27OGNN94o8f1y2717N0899RQNGzZk2rRp0n3Wu7l13VAQ28KPBbVMzp07l4oVKxY5eShO96/iTktsT1an9z4XLlxg6dKlREZGcsMNN5T686+XqIQppT5QSn1g/71S6kOgVinE51JiXtzK2Ze308a3IQA//7Kbsy9vJ+bFrQZHJoTLs83ra1ZKVQYysX6C6hacMT2xyWS67srxhR1Qf/bsWX755Rduv/32nFhL6t5776Vly5a89tprxVp4Mj9Xr15l8ODBgLXriiwc5/Xcum4oyPUG1P/222/s37+fQYMG5XwYUli236GizP61Zs0aKlSoQKdOnYr0LHuyOr33+fLLL8nIyCj1QfQ210tUpgJXszf7768Anzs1Mhe09sxMLqSeIKhOA0wBZdh27jfOp55gbcwso0MTwtUdyn4T8g2wE9gFFH7eXYM5YzB9cHDwdVsSKlasSP369a/bovL999aeMiXt9mXP39+fOXPmADBixAjS0tJKfE+tNWPHjuXYsWN88MEHOYmY8GpuXTcU5HpTFM+dOxcoWrcvG/vuX4WZoS82NpZffvmlWNMS2wsNDSU9PZ34+Phi30O4D4vFwvTp06levTp9+/Y1JIYCExWt9X8L2korSFdx+2svc7FSGt/UCSSsehP2n/uTCyEp3P7qS0aHJoRL01oP11pf0lp/DIwCXgSGGxtV4TljjEphZ7hq164df/31V4H93FeuXImPj0+JunTkpXXr1rz00kscPHiQl14qeT0XFRXFokWLGDhwIA8//LADIhTuzt3rhoLUrVuXsmXL5tmiYrFYmDdvHrVr16ZLly7Fuv+gQYOIjY0tVPev6OhooGTdvkBm/vI269ev5/jx44wZM8ZhrfVFJdNVFUGnw5cZ27kdvhcyaFypKZlpqYypGULnwx4xm6IQpUJrvUVrvUprXfRFAAzijETlegPpbdq2bYvWmr179+Z5PiUlhfXr19O1a1cqVqzokPjs/fvf/6Z9+/a8++67bN++vdj32bt3L5MmTSI8PJz//e9/Mi5F/IM71g0FMZlMRERE5JmobN26lVOnTjFs2LBizxxalMUfSzItsT1Z9NG7REVFoZQydMITSVSK4N9LLjNl6WWan86gfvVmAHSft4MpS62JyvGUNO7YdYjjKSXvIiGEJ1BKXbZfddpuu5y9CrVbcEbXr6K0qED+41Q2bNiA2Wx2aLcve76+vsyZMwc/Pz9GjhxJSkpKke+RmJjI4MGDsVgsLFiwgJCQECdEKtyJp9QN19OiRQvi4uKIjY295nhJun3ZRERE0KRJE5YsWVJg96+srCyio6Np06YNNWrUKPbzQFpUvElMTAwrV67k9ttvp169eobF4fRERSn1iVLqhFJKK6VaZx8ro5RappQ6pJTap5Rar5RqaPeaH5VSx5VSe7O3SXbnqiml1iqlDiulflNK3ezsMtjc8+92JJdPw4KmXnaicjj+AHc+3YbEzCzePn6OPYkpvH38HImZWZizZDYw4fVac+2q07bNdtwtOHrWr6K2qED+M3+tXLkScOz4lNyaNm3K66+/zuHDh3n22dxLahVMa82ECRM4fPgw77zzDu3bt3dSlMLNtMYD6obryWtAfXp6OosWLaJ58+YlWj9IKVWo7l8///wz8fHx3HHHHcV+lo2sTu89Zs6cSVZWFuPHjzc0jtJoUVkMdAFO5jo+HWistW4FLAf+l+v8JK116+ztQ7vjbwE7tNaNgNHAPKVUqXSc+9e0jWzPzEIBgUHBVC1fi7/O/kGn4ydo9NMBlsVeAWBZ7BUa/XSABpv3l0ZYQrisXCtO/2MzOr7CctZg+sKoXLkydevWzTNRsVgsrFy5kiZNmtCoUSOHxJafiRMn0qVLFz755BM2btxY6NfNnDmTefPm0bdvX5544gknRijciafUDddjm6LYPlFZu3Yt8fHxJWpNsSlM9y9HTEtsIy0q3iEzM5MZM2YQFhbmkAS3JIqVqCil5mS3lNS53rVa681a6zO5jpm11qv13/Pb7QDqFfLxg4Go7Pv8gnVxqG6FDr4EyjZ+iZbqIhrNgZqbqFe9KclXzhFy6Cm6VCiLT/Z1PkDXimXZ0KFxaYQlhMvLbiE9lnszOq7CcuT0xFprUlJSCt31C6ytKn/88cc/ul3t3r2b8+fPc/fdd5c4ruvx8fFh9uzZBAUFMXr0aBITE6/7mt9++43HHnuMOnXqMGvWLBmXIv7B3euG68lr5i9bt69hw4aV+P7NmzenSZMmBc7+ZZuWuGPHjiV+XtWqVfHx8ZFExcOtXr2aM2fOMHbs2BLNEucIxW1R+Q44DXzkoDiewNqqYu8tpdQBpdQCpVR9gOwpDP201uftrjsBXDdhcoSFd84jMyiD6BtmsbPuSlJbXgVgYlB/bqtQDosGvyyNRcNtFcrR2L/gNRKE8CJ3AX2yt0HACuBLQyMqAke2qKSmpqK1LnSLCljHqVgsFvbt23fNcUetRl9YDRo04N133+XkyZM8+eSTBV6bnJzM4MGDyczMZMGCBVSqVKlUYhRux63rhuupWrUq1atXz2lRSUhIYMWKFXTt2pW6deuW+P627l8XLlzgp59++sf5CxcusGvXrhJPS2xjMpmoWbOmJCoeLioqCh8fH8aMGWN0KBTrf63WepmjAlBKPQc0BCLtDt+vtT6trB+/PQKsApoV8b6Tgcm2/eDg4Jzp+XIzm835nrPn4/sIP1TpRXLIAUxakdwyHpbAwS92cLB8ZyhvYtipdL6u5883e8/QK+ovfu+RWpSwS11hy+6pvLn8pVl2rXXupZl3K6W2Aa+WSgAlZPsD74hEpbCr0tuzjVPZs2cPnTt3zjm+cuVKKleufM0xZ5swYQLfffcdM2bMoH///vl2J3nkkUf4448/eOedd0q0wJzwbO5eNxRGixYt2L59OxaLhe+++w6z2eyQbl82gwYN4tVXX2XRokV07979mnOOmpbYXmhoqIxR8WAnTpxg7dq19O3bN2eWNyMVmKgopQrsT6C1XlGShyulngLuAXpqrXP6NGitT2d/1cBnSqn3lFKVtdaXlFKZSqkadq0q9YBTecT2AfCBbT8sLEz37t07zziio6PJ75y9xMTv0ScfQinoUS6dDWFlMAX4sE8fItzkw5i9qXSPzaTlFQvrGgZSfWJbatco/JsRIxS27J7Km8tvZNmzW0dLNv1MKVJK4evr65BEJTk5GaDILSpw7YD6U6dOsW/fPkaMGIGPj09+L3U4k8nEzJkzadGiBQ8++CC//fbbP6ZFnjNnDnPmzOGOO+64bsuLEPbcrW4ojObNm/PDDz9w4sQJ5s6di5+fX87YEkfdv3HjxixZsoRPPvnkmvrAUdMS2wsNDWX37t1kZWWVat0jSseMGTNyJkFxBdfr+jXJbvsm1/7Ekjw4u8VjGHCr1vqK3XFfpVR1u/0BwAWt9aXsQ4uACdnn2gO1gE0liaWwypa9gRtrVmJMlTT6VlSMqZZOjcaV2HPgIB8El6VbbAZZlky6xWbwSUgFTJX8SyMsIVyeUupXpdSe7G0fcJR/TqDh0vz8/BzaolKURKVatWqEhYVdM0Vxacz2lZ86derw0UcfcfbsWR5//PFrzv3+++88/PDD1KpVizlz5hR7jQjhHTyhbrge2ziV9evXs2HDBm6//XaHdoXMr/tXVlYW69ato23btiWeltheaGgoWVlZxMXFOeyewjVkZGQwc+ZMwsPDufXWW40OB7j+yvQ9bBtw1H5fa31LYR6glJqmlDoDhAHRSqkjSqkw4H2gArAxewrindkvCQC+zx6fsg94GLBv2ZkC3KSUOgzMBoZrrR0zFc/1y8KoKhm0TtTU/KgKbZI0vdpWICUlhQ2zFoPWHErYBVpzatnPfDJyYGmEJYQ7mMjfH3I8DDTVWr9haERF5O/v79AWlaJ0/QJr96+DBw9iNpsBa6Li5+dHr169ShxTcYwaNYq77rqLb775hqVLlwLWxSeHDBlCWloa8+fPp0qVKobEJtzKRNy8brgeW6Ly2muvYbFYHNrtyyav2b9s0xI7stsXyKKPnmz58uVcuHCBcePGucyHTEUZo6Kvf0keL9I6vwmY85z+RWudDNxYwP0uAIb8Zc4ymwm01KTyT2VJ/2sPtX9qz40tLjObw+y9eBgdEEtM8iEupZ+jSc1O3P/2J0aEKYTL0VqXSqunM/n5+Tlk1q/idP0Ca/evFStWcODAAZo0acLGjRvp0aMH5cuXL3FMxaGUYvr06TRv3pwJEyZQq1Yt+vTpQ1xcHK+//jpdunQxJC7hXjyhbrieZs2sQ2zPnDlDcHCwU1pBW7RowQ033HBN96/Vq1cDjh2fAtdOUWwbPyc8Q1RUFH5+fowePdroUHK4RrrkJg61bYfvyE2kb7B2v0jf8CttvzgKwL60Q8QkH8JHWYhJPoTubKJqTY/qZitEsSmlmmXP4LdHKbXfthkdV1E4uutXcVpUwDpOZd26daSnpxvS7ctezZo1+eijj4iLi6NXr17ExcVRo0YNHn744ZyWHyEK4gl1Q0HMZjNZWVk5H0xUrVqVjIwMh/9+KKUYPHgwFy5cYMuWLYB1fErFihUdMi2xPVn00TMdOnSIDRs2cM8991C9evXrv6CUFJioKKVa2jYgUCnVItcxr1J/+TKCOnUC2+AxHx9Cu3ShUUUfdu/ehULTtlIMCs3BRZ/CG6HGBiyE65gP/IK1m8dDdpvbcFSiUpIWFbAmKkaOT8lt1KhRgHXaVYDz589TsWLFIidiwmu5fd1QkODgYEJCQnJ+70+cOEFISIhTfj/su39duHCB3bt3O2xaYnuy6KNnmj59OoDhK9Hndr0WleV2Wxms85vb9pc5NTIXFNCoEeV69gSLBeXnBxYL5XpGctOtdxGbZKZNyGliTRO4peZ5QiqUQ0/YYnTIQriKLK31e9kLwG61bUYHVRRGDqYHa+tFjRo1+OWXX/j+++9p2bKlQ9ZhKKn9+/fTtWvXnMUcfXx8iIyM/MeaL0Lkw+3rhoLs37+fyMjInP7+zvz9sO/+5axuXyCJiicym83Mnj2bG2644R9TXBvteoPpwwvY6pdWkK7k6nffgVJUGjkClOLqd9/RoYt1CZhN55pxIfMGYulD7363k1WunrHBCuE6NiqlbjY6iJJwdItKcT5RbdeuHfv27ePixYv861//KnEsjhAREcHgwYMBCAgIwGKx0L9/fyIiIgyOTLgJt68bChIREUG/fv3QWjv998M2+9f58+d56CFro5QjpyW2qVixIgEBAZKoeJAlS5Zw6dIlJkyYkPOhk6uQMSpFoLXGr3Zt/F9/h3GBnfF/7R38wmoTs9m6Av22MwpQHDbfzIwFTZk+2e5DoUtHYUak9asQ3mcx1tn8ziqljimljiuljhkdVFE4etavoraomM3ma97cHDlyhISEBJcYCzJnzhyUUjzxxBMopZg9e7bRIQn34fZ1w/WU1u+H2WzOaUFJS0ujYsWKBAYGOmU8TGhoqCQqHiQqKoqAgABGjhxpdCj/4NiOix4uLdNCyFvv8Nx3B9h7+hwfVarJG2+9w8O/7efDFWU4fsG6wK4ii1qBf9Hlga6QYYasdNj4OsTssn696yPw8Qe/MsYWSIjS8yXwBLALyDI4lmLx8/PLGYdREsUdTB8cHIzFYsnZX79+PSEhIfj4+JCZmVniuIpLa014eDjPP/88ffv2pXPnzsybNw+ttct9MidcktvXDQUpzd+P3HXE5cuXnVZHhIaGcvjwYYfeUxjj4MGDbNmyhfvvv9+h6/s4iiQqRdDsxbVY7CZpXrnvHCv3ncPko2jRoCG/Hv6LWt1eJXb3WGpHVCG4WQt4rTLovysOflti3ZQPvBRf+oUQwhhJWutZRgdREkZPT7x//34mTJiQM6OPj48P3bt35+OPPy5xTCWhlGLhwoU5+/369aNfv37GBSTcjdvXDQUpzd+P/fv388QTT/B///d/aK0xmUz06NHDKXVEaGgoW7duJSMjAz8/P4ffX5SeadOmAbjMSvS5SdevIlg78WZualAZH1P2oFGT4l8NK5MVWYUa5euSkZnO2SuHqdJ8Ob/9lkiTTbvgoW0Q3s2amID1a3g3eMhjxgoKURjfK6WMn6KqBIweTB8REcGQIUNQSuHv7y9jQYSncPu6wVXYxsOAtb7SWjutjrANqD9//rzD7y1KT3JyMl999RUtWrSgc+fORoeTJ0lUiuCG6uXoHVEDi0Xj72PCYtH0jqjBTH0vHfvuBuCPP8yE1NlFrTvfYoa+F6o1hSZ3WVtVfPytX5v2sR4Xwns8BixXSiUrpeKVUpeVUm7VpGj09MTwd1/3iRMnylgQ4Sncvm5wJbY6YtKkSU6tI2R1es+wYMECrl696pKD6G0kUSmixbvPoBSM7lIPpWDRrjMc+7/JVKjbGICDB81kYSImuT7HN062vmjfPFAKOj1k/bp3roElEMIQrYFwoBnQJnu/jYHxFJkjE5UyZcrgY1uPqZBsfd2XLl3K22+/zZIlSwgPD0drff0XC+G6WuPmdYOrKM06QqYo9gzTpk0jKCiI++67z+hQ8iVjVIpAa02dSoE8dktDekXUoG2diqzYG0ODpu+z0tSX2nW28fvvaSit2RPcgj7d3gf9EFSoCzc/DU3uhLAOcGARaG1NWoTwAlrrk0bHUFKO7PpVnKmJZSyI8ESeUDe4itKsI2R1eve3Z88efv75Zx588EFCQkKMDidfkqgUgVKKqfe1y9nvHVGD3hE1uLD/c85feIW/Isqwdk0iZ2PSaV9zJ82rf25NRgbP+fsmTe+ybkJ4AaXUt1rrYUqpX4F/fKyntW5rQFjF4sjpiYvT7UsIT+JJdYM3khYV9+fqg+htJFFxgKjvttO4SQB16nQC1rN4UVXuGx5A1LbtvNSyt9HhCWGk97K/TjQyCEewDU7Nysoqcrcte0lJSZKoCOFBdYM3kkTFvSUkJDB37lzatWtHu3btrv8CA8kYFQdYH/oDZ082o1rV9gQGBvLzz2c4e7IZP4RuMDo0IQyltd6d/XWTbQP2Apeyv3cbtik4SzpFcXJycrG6fgnhSTypbvBG5cqVo2zZspKouKl58+aRnJzs8q0pIImKQ0TdPYOKfqcw+fjSuEljYmNjSb10kC/unm50aEK4BKXUWqVUBaVUWWAfsEop9YrRcRWFLVEpafcv6folxN88oW7wVrI6vXvSWhMVFUX58uUZOnSo0eFclyQqDnD3R8fYk9AagAodKgKwfKeFvh8dz7km82IqF6buJfNiqhEhCmG06lrrK8AdwHKgEdDf0IiKyFGJSnEH0wvhody+bvBWkqi4p507d7Jv3z7uv/9+t/hbJImKA6x5oivpgVfYHfwHVwLvxb9sAId+38/3j/8LnWHBYs7k6vqTZJxO5Or6k1jMmegMy/VvLITnsC1dfDOwXmudAWQaGE+ROSJRSU9PJyMjQ1pUhPib29cN3io0NJTLly+TmiofwLoT2yD68ePHGxxJ4chgegeoVdaXi7WyuLirB5f9a1CxQScu7NtE/PHfifkk85r5TFL3xZG6Lw4UhL3Z1bighShdvyml1gBNgWeUUkFGB1RUjkhUSrLYoxAeyu3rBm9lW/Tx3Llz1K9f3+BoRGFcvnyZ+fPnc9NNN9GiRQujwykUaVFxgBavbmD3gTs4H1DHOh1x8zsAuGPy+1Sf2JaABiF//6RNENCwAtUnysyLwquMAqYBPbTWKUBF4FlDIyoif39/wDGJijs0twtRSkbh5nWDt5KZv9zPV199hdlsdotB9DaSqDjAt8MaEZp5EaWt3bnKhDbGv2xFgk7vwLdaEIERVcACFrLAAoERlfGrLp+oCu+htTZrrZdprY9n78dordcaHVdRSIuKEI7nCXWDt5JExb3YBtFXrFiRgQMHGh1OoUmi4gCdWjfhziaV0Ch8dBYoH9p37sa5szFs37aNpF/OodEcSUlGo0n6+RyZJZziVAhRuhwxPXFSUhIgLSpCCPcnq9O7l59++ok///yTUaNGERgYaHQ4hSaJioN8fyQRhea2qmYUmqs1rF27nn1gJMcO/8qOKyf5Pa0MO66c4tjhX/lk5BCDIxZCFIW0qAghxN+kRcW9REVFAe4ziN5GBtM7gMVioXqAZnzbEEb27cPsZZtYtiec5PBw9sdcZXdGPUzKhFIQSyixmVCmwmNGhy2EKAJHJCq2FhVJVIQQ7q5mzZqAJCru4MqVKyxevJgePXrQuHFjo8MpEmlRcQCTycR3/72PYa3rcnzwEO5tU49lrwxnyJAhXEm+SGzqPnT2+BWtLYRUzWLICx0MjloIURSObFGRrl9CCHcXGBhIxYoVJVFxA+vXrycjI8PtWlNAEhWHsKSlkZWURNzHn2Dev5+4jz8hKymJwffcA8DW/cvQWVdJS5hn/Zp0kMqh8kZFCHfiyFm/pEVFCOEJZNFH12exWFi9ejVVq1alf3/3W0tVun45wF9t2oLl7wUcE1avJmH1agJMJm644QZ2HTrAHc22QNZ5Ms3byArsTEZaGn4BAQZGLYQoCkd2/ZIWFSGEJwgNDWX79u1GhyEKsGHDBs6dO8e///3vnA/c3Im0qDhA/eXLCOrUCUzZP06TiaDOnWmwfBn1/DSJqVc5cnY7CrBk/EVi7Gw+HTXI0JiFEEUjg+mFEOJatWrVIikpicTERKNDEfmIiopCKcXYsWONDqVYJFFxhPp18O/RFW2xkOED2mLBv0cXqF+HmPYPALD//Bka9TuOX0gGZwLDmF9rsMFBCyGKwpHTE0uiIoTwBDLzl2s7e/Ysy5cvp23bttSvX9/ocIrF6YmKUuoTpdQJpZRWSrW2O95IKbVNKXVIKfWLUiqipOeM0mFuB3Z/+S5awer2Cq1g16x36Ti3Iwn396B8rVD2nz6PX8VUat4Yy7H2TTl/eyOjwxZCFIEMphdCiGtJouLaZs2aRVZWFnfeeafRoRRbabSoLAa6ACdzHZ8GTNda3wC8Dcx2wDlDLL5rMZk1KvPBAF/m9vDhw3t8yapRhUV3LeJLhjPwtlSSU7L4+ecUKjVM4PGmn/ClGm502EKIIpCuX0IIcS1Z9NF1ZWVlMX36dEJDQ+nYsaPR4RSb0xMVrfVmrfUZ+2NKqWrAjcA32YeWALWVUg2Le87Z5ShIo0qNSHnpYX5ppPAz+fHzDYrUlx+mUaVG3NjuO9q3rUtAgGL27HgyMuFibFVubPedkSELIYpIBtMLIcS1pEXFda1Zs4bTp0/z4IMP4uPjY3Q4xWbUGJXawDmtdSaA1loDp4A6JThnqOVHl6OAe8y9UMDyI8sB2LW7P3826UGve0I5cTyD7zaY2F7tJnbtdr8p4oTwZo6antjHx8ctZ14RQojcJFFxXdOmTcNkMvHggw8aHUqJeOz0xEqpycBk235wcDDR0dF5Xms2m/M9VyiZGr9EXyYljCEypg21atVkg//PRH+/lkq7/kOrzl9zeMhjqFVv8vWXl/io6x7K7f4P0RkleKaDlLjsbs6by+/NZS8OR7WoBAcHo5RyVFhCCGGYGjVqAJKouJqTJ0/y/fff06dPH2rXrs3vv/9udEjFZlSichqoqZTy1VpnKutf7TpYW0cSinnuGlrrD4APbPthYWG6d+/eeQYTHR1NfucK48yzP/GWnpCzf2tMW26NaQsKOvRsicX0PvhB8PBYkr74gCdXtaXcwBac7dmu2M90lJKW3d15c/ndpexKqXVADcACJAKPa61/zXVNd2AN8Jfd4c5a61RHxeGIWb+Sk5Ol25cQDuIqdYM38/Pzo1q1apKouJj//e9/aK2ZMGHC9S92cYZ0/dJaxwJ7ANuI8gHAGa31keKeK73o/6n6xLYENAj5+6dpgoCGFag+sS0fmRNofOwQpqwsgvoOwad6TcxfT+e1S/JLLUQhDdZat9Rat8b64cPsfK77S2vd2m5z6BsRRw2ml4H0QjiMS9QN3k5Wp3ctGRkZzJw5k7p169KrVy+jwymx0pieeJpS6gwQBkQrpWxJxXhgvFLqEPBvYLTdy4p7zhB+1YMJjKhi/UzHR4EFAiMq41c9mOZPP0zjS8fRSuGvFGVHPURmSjJbHzE8bCHcgtb6it1uCKCNiMNRXb+kRUUIx3CVusHb2RIV67BhYbSVK1dy7tw5xo0b59aD6G1KY9av8VrrMK21r9a6uta6Yfbxv7TWnbXWN2itb9RaH7B7TbHOGSl59wVQULZLLVCQvOsCAOFLv2N/zWYAdN71C2VuuY2g2vWYm5jI6dOnjQxZCLehlPpKKXUaeBW4P5/LGiil9mSvr/Swo2OQFhUhXI8r1A3erlatWqSlpXH58mWjQxFYV6L39fXlgQceMDoUh/DYwfSlSWuNb6UylL+lDoERlQmoU46UvXForfFr2IAKG3bRfv82ypRNYtD6M2zrM4g9n7/Lyy+/zMyZM40OXwiXp7UeAaCUGol1/aQ7cl2yBwjTWl9VSoUBq5VSF7XWC3Pfq7gTbdg+WPjzzz+LPQnB1atXSU1NdYtJDLx5sgUpu/uU3RXqBk9V2DKmpKQAsGjRIurVq+fkqBzL0/4dz549y/r16+nSpQv79u1j3759gHuXUxIVB1BKUfm+pjn7gRFVrF3BgLSsNN7aspRvG4RQRtfGFHCa/525wsRuNzN79mwmT55MRESEUaEL4Va01nOUUlFKqcpa60t2xxPsvj+jlPoW6Ar8481IcSfaOHbsGAB16tQp1iQEFouFtLQ06tWr5xaTGLjLZAvOIGV3v7IbWTd4qsKW8eTJk8ydO5fw8HC3GxPhaf+OU6ZMAeDll18mMjIy57g7l9OodVS8xjtvvsOyxi0I9K2DQhHoW4fljVvSqFEjLBYLzz77bJ6vy7yYyoWpe8m8KGP+hPdSSlVQSoXa7fcDLgHxua6rqZQyZX9fDrgLuGb2n5Iqadev1FTr77J0/RKi5FypbvB2sjq9a0hLS2PWrFk0atSIHj16GB2Ow0ii4mR9hvfBXN6MBQsAFiyYQ8xMmDKBwYMHs3LlSn766aec63WGBYs5k6vrT5JxOpGr609iMWeiMyxGFUEII4UAy5RSB5RS+4BHgbu01lop9T+l1N3Z1w0AbNfsANYDXzoykJJOTyyr0gvhUC5TN3g7WfTRNSxdupSLFy8yfvx4TCbPeXsvXb+c7MaGN7Kz4U5i98SSRRYmTNRpWIcbG97Ia6+9xtKlS3nmmWfYtm0bSiliXtx6zbwlqfviSN0XBwrC3uxqXEGEMIDW+iTQIZ9zD9p9/xnwmTNjKWmLSnJyMiAtKkI4givVDd5OEhXXEBUVhb+/PyNHjjQ6FIfynJTLhZ3607oeZWCDQABO/nESgEaNGjFu3Dh27NjBd999BxS8JosQwjglTVRsLSqSqAghPEnVqlXx8fGRRMVAf/zxB5s3b2bQoEFUqVLF6HAcShIVJ7NYLPgG+xIRGcFz9z9Hs8hm+Ab7YrFYu3K9+OKLBAcH8+yzz5KZmVngmixCCOM4qkVFun4JITyJj48PNWrUkETFQNOmTQPwiJXoc5NExclMJhOvP/I6g7sOBmBI1yG8/sjrOf0Hq1evzlNPPcWhQ4eYNWsWYF2TRStYH6TRdmuyCCGM4+/vD0jXLyGEyE1WpzdOamoqc+bMoVmzZvzrX/8yOhyHk0TFBTz55JNUq1aNl156iYuXr6JD/FlUK5NdaTtYXCsTHeJPanqm0WEK4dV8fHxQSpW465e0qAghPE2tWrU4d+5cTm8RUXoWLlzIlStXmDBhAkopo8NxOElUXEC5cuV48cUXOX/+PI37PUrX34/y+7kDVFHJHDx3gK6/H6XlS2uNDlMIr+fn5yctKkIIkUtoaChZWVnExcUZHYrXiYqKIjAwkPvvv9/oUJxCEhUXMXbsWBo0aEDSjoX0y/qJML94lIIwv3juK7OH4WV2Gx2iEF7Pz8+vxNMTS6IihPA0MvOXMfbu3cuOHTsYNmwYFSpUMDocp5BExUX4+/vzxhtvkJ6eztYff/57imINfmkVqHhRZv0SwmiOaFGRrl9CCE8jiz4awzaIfvz48QZH4jySqLiQgQMH0jSoHD/v2crly5fBYs1WKsZn0nHbhwZHJ4SQrl9CCPFP0qJS+hITE/nmm29o06YN7du3Nzocp5FExYWYTCYmvfUcWZYsNm1eR+uWKyhTJoErIVcp8+kzRocnhNfz9/eXwfRCCJGLJCql79tvvyUpKcljB9HbSKLiYq7s+Ikbwmuzb9+fvHRpKOVq/YXiIuumfmp0aEJ4PWlREUKIf5JEpXRprYmKiqJcuXIMGzbM6HCcShIVF9N8zDEeedYflOLP6d+xvlZvmt/2E83HHDc6NCG8XkkSFRlML4TwVJUqVcLf318SlVKya9cufv31V+677z7KlStndDhO5Wt0AOJaDzAX3dBEmV4vYo5eyf+t6cWOO7/BpLOQX38hjOWIFpWgoCBHhiSEEIZTSsmij6UoKioK8MyV6HOTFhUX8+Gs5wlPOEzIQ5MwVatB0tR3qfHHj3ww6wWjQxPC65VkeuLk5GQCAwPx8fFxcFRCCGG8WrVqSaJSCq5cucK3335Lp06daNWqldHhOJ0kKi6mfa2n6eC7E1UuhIr/fhmdlsa5D9+iVc2JRocmhNcradcvGUgvhPBUoaGhxMbGFruOFIXz9ddfk5qa6tFTEtuTRMXFzAn+P3aU6YwGyoXHEjR0FLFHLvDgjv8aHZoQXq+kXb9kfIoQwlOFhoaitebChQtGh+KxtNZMmzaNChUqMHjwYKPDKRWSqLiYZXWXcy71IuUvfkzZqwsJvSsDv4Z12bZ1Gz/++KPR4Qnh1Uo6PbEkKkIITyUzfznf1q1bOXjwICNHjvSa8Y6SqLiYJTHn6HXmXYJSdgMQnP4r/YZbCPSF+++/37oQpBDCECVtUZGuX0IITyWr0zufbRC9t3T7AklUXE7DMZtomxKEBfCzaCxA93IV+eTtVzhz5gzjx49Ha210mEJ4Jen6JYQQeZMWFee6ePEiixYt4uabb6Zp06ZGh1NqJFFxNdWasqZiMAoYnpSMAtZUDOKBJ/7DPffcw6JFi/jqq6+MjlIIrySD6YUQIm+SqDjXnDlzSE9P94opie1JouKCbtDpfBh7iclN7ueDuEvcoNNRSjF9+nRCQ0N59NFHOXr06D9el5Jygl92DSAl5UTpBy2EF7AlKkVt1UxPTyczM1NaVIQQHksSFeexWCxMmzaNKlWqcM899xgdTqmSRMXFKOD9ss255e7/wa2vEHnXDN4v2xwFVK5cma+++oqkpCTuu+++nE92s7LSyMxM5OixD0lI2MvRYx+SmZlIVlaaoWURwtP4+fkBkJmZWaTXyar0QghPV65cOYKDgyVRcYKNGzdy+PBhRo8eTUBAgNHhlCpJVFyNUjB4DjS507rf9C7rvlIAREZG8tRTT7Fz505ee+01AH7c1JxNm1sTG7sKgNjYVWza3JofN7UwpAhCeCp/f3+APLt/HTlyhE6dOnHkyJF/nLOtSi9dv4QQnkop5fRFHwuqZz3ZtGnTABg3bpzBkZQ+wxIVpVRlpdReu+2QUipTKVVJKfWjUuq43blJdq+rppRaq5Q6rJT6TSl1s1FlcKYrsSksfmsXV2JT/nHutddeo3Xr1rz22mts3bqVG9uuJKR8J87pUF7kTc7pUCqEdObGtisMiFwIz2VrUbFPVMxmMwkJCbzwwgvs3LmTF154gYSEBMxmc841tkRFWlSEEJ4sNDTUKYlKYepZT3X+/Hm+++47br31Vho2bGh0OKXOsERFa31Ja93atgHTgTVa6/jsSybZnf/Q7qVvATu01o2A0cA8pZRf6UbvPJkZWaSnZvLzimNcOJHAzyuOkZ6aSWZGVs41AQEBzJs3j4CAAO677z6mPnOQX7fUY7EazFEasVgNZvdPdZn14gnjCiKEB8orUQkODiYkJIT58+cDMH/+fEJCQq5pPbF1/ZIWFSGEJwsNDSU+Pt7hCURh6llPNWvWLDIzM71uEL2NK3X9GgPMLMR1g4EoAK31L8BZoJsT4ypV0x/fxIxJmzm8KxaAw7timTFpM9Of2HTNdU2bNuWDDz7g5MmTvHpsJi+0GsQO1RWUYofqyoutBvF23+pGFEEIj5VXorJ//34iIyNR2d0zlVJERkayb9++nGukRUUI4Q1sA+rPnTvn0Pt+//33/1jgsF27dtfUs54oKyuL6dOnU6NGDfr06WN0OIZwiURFKXUTUBFYZXf4LaXUAaXUAqVU/ezrKgN+WuvzdtedAOqUWrBONuSFDoQ1rojK/pdRJghrUpEhz3f4x7Xjx4+nT58+mDesIWTtTJS2WF+jLdyQ/heDd20ozdCF8Hi2RCU9PT3nWEREBP369cuZCUxrzV133UVERETONTKYXgjhDZwx89eWLVu4//77SUmxdoW31cMHDx4kLc2zJw1at24dJ0+e5MEHH8wpt7fxNTqAbGOAr7TWtql07tdan1bWjygfwZrANCvKDZVSk4HJtv3g4GCio6PzvNZsNud7zgjmABPaYgITaAuYAy6y60AsHPjntffddx9bN/zEsU9mUqHVbfhVr4FFmahxDIZeanLdcrla2UubN5ffm8teXHm1qIB1fnt7U6dOZeLEiTn7MpheCOENHL06/VdffcXYsWPx9/enYcOGHDt2jEmTJvHuu++SlpZG7969+emnn2jSpIlDnudqoqKiMJlMjB071uhQDGN4oqKUKou1O1d72zGt9ensrxr4TCn1nlKqstb6UvaA+xp2rSr1gFO576u1/gD4wLYfFhame/funWcM0dHR5HfOCAvf+AWlEmndsw57159CXa5A797t870+4Owp+k9+hqtvPk/3SVPYX68xv9UMpXVrTVDngucacLWylzZvLr83l7248kpUtNaEh4dz7NgxgoKCOHPmDGlpaWitc7qDSdcvIYQ3cFSLisVi4fnnn+fNN9+kTp06rFixgtdff5333nuPvn370rlzZ95//322b99Oz5492bJlC/Xq1XNACVzH6dOnWbVqFXfccQd16nhMx6Eic4WuX0OAfVrrPwGUUr5KqZzBFUqpAcAFrfWl7EOLgAnZ59oDtYBrB3C4Ma015auU4bbxLbjpnobcNr4F5auUKXCBuRumzaTmbX3J+G0vGR99wi17fqJSShInHvC+aeyEcKa8pidWSjFnzhzi4+Pp0qULnTp1+sdgUhlML4TwBo5IVJKTkxk0aBBvvvkmHTt2ZOfOnbRq1YqFCxfSt29fAPr168dPP/3E119/zdmzZ7n11ls5f/78de7sXmbOnInFYmH8+PFGh2Iow1tUsHb7mmG3HwB8r5QKACzAReBuu/NTgK+VUoeBdGC41vqfixq4KaUUt437e/2T+q2rUr911QJfE/rKFCZtPcG8A7/z44Ht1Hz1EPf07kvof6c4O1whvEp+Xb+OHTsGQMOGDenYsSM7duxg9erVDBgwAJAWFSGEd6hZsyZQ/EQlJiaGu+++mz179jB06FBmzZpFYGBgvtcPGzaMhIQEJkyYQK9evdi0aRMVK1Ys1rNdSWZmJjNmzKB27drcfvvtRodjKMNbVLTWN2mtv7TbT9Za36i1bqG1bqW1jtRa77M7f0Fr3Utr3UhrHaG13mhM5K4j6deGDAvsyfL7ZvBy5GNcTkli6qJZPPnO9zlvkIQQJZdfomJbfKxBgwYMHDgQgIULF+acl8H0QghvEBQURIUKFYqVqOzevZsOHTqwZ88eXn75ZebNm1dgkmIzfvx43n77bQ4cOMAdd9yRU9+6s1WrVnH27FnGjRuHj4+P0eEYyvBERZTckM6BZFz5A6U1Y24cRPSo/9GuegNm7V5Cq1at+Omnn4wOUQiPcL1EpWHDhoSFhdGlSxdWrVqV80GBDKYXQniL4qxOv3TpUrp27cqlS5f49ttveemll3LG+BXGM888w7PPPsuOHTvo16+f2y8EGRUVhY+PDw888IDRoRhOEhUP8E6LaqQmHwSTiUydRXjlOsyOvJ3H/vMCMTExdOvWjYkTJ+ZM7SeEKJ68picGOHr0KGBtUQEYMmQIKSkpfP/994B0/RJCeI+irE6vtebNN99kwIABlC9fnh9//JGhQ4cW67mvv/46Dz30EBs2bGDo0KFkZmZe/0Uu6NixY6xbt46+ffvmjPnxZpKoeICOtWqQ1nQgGjh89Rc0kBExmA9efJ59+/bRuXNnPv74Y1q1asWWLVuMDlcIt1VQi0pQUBA1atQAYMCAASilWLBgASCD6YUQ3iM0NJTExEQSExMLvC4tLY1Ro0bx3HPP0bJlS37++Wc6depU7Ocqpfjss8+49957Wb58OQ888AAWi6XY9zPKjBkz0Fp77Ur0uUmi4gE+GTGA1LN/sPXCUvZf3sS2C9+RevYPPhk5kBtuuIHNmzfz/vvvc+bMGW6++WYmTZokrStCFEN+icrRo0dp2LBhTleFmjVr0q1bN1avXk1iYiLJycn4+vrmzBomhBCeqjCr08fFxdGzZ0+++uor7rrrLrZs2eKQKXhNJhOzZ8+mT58+fP311zzxxBMFzprqatLT05k5cyYNGjQgMjLS6HBcgiQqHiBDa5brX4hJsfaTP5NymOV6FxnZnyT4+PgwefJk9u7dS6dOnfjoo49o3bo1W7duNTJsIdxOXtMTp6enc+LEiZxuXzaDBw/GbDazcuVKkpKSpDVFCOEVrrfo4++//07Hjh3ZsmULTz75JMuWLaNcuXIOe76fnx8LFy6kR48efPbZZ7z44osOu7ezfffdd8TFxTFu3DhMJnmLDpKoeISbp7xGcGoqoMkyaUATZE7h5imvXXNd48aN+emnn3j33Xc5deoUXbt2Zfr06aSmphb6WX+c+Ytb3/mSP8785dhCCOEG8mpROXnyJBaLhYYNG15z7YABAzCZTCxYsIDk5GQZnyKE8AoFraUSHR1N586dOX36NNOnT+e9995zyqxWZcqUYfny5XTo0IHXXnuN9957z+HPcIZp06bh5+fH6NGjjQ7FZUii4gEm/DaSVJ8LAPxeLwEAs+kCD/026h/X+vj48NRTT7F37146duzI0qVLad26Ndu2bSvwGSnmFC4nXebTtcsZ1PBrPl27nMtJl0kxSxcy4T3ySlRyD6S3qVatGj169GDt2rWcO3dOEhUhhFfIL1H57LPPuOOOO/Dx8WHdunWMHTvWqXGUK1eONWvW0Lx5c55++mlmzJhx/RcZ6M8//2Tjxo0MHDiQqlULXj/Pm0ii4gEW37WYdD/FxjZx7G5yhY1t4kj3Uyy6a1G+r2nSpAlbtmxhzJgxnDx5ki5duvD000/n27rS6pW1dH5zA+GBO6gfcpLwwB10fnMDrV9d46xiCeFy8kpU7Kcmzm3IkCGkp6dz6tQp6folhPAKuROVzMxMHn30UR577DEaNGjAjh076NGjR6nEUqlSJdatW0f9+vUZP358zgQnrmj69OkAXr8SfW6SqHiABhUaEFurIqdqpuJn8uNUzVRia1WkQYUGBb7Ox8eHQYMG8euvv9K+fXvee+892rRpw44dO/5xbVTPZ5gaOYXQmhd4Sb1JaM0LTI2cwheR/3ZWsYRwOXlNT5xfiwpA//79c7o1HD58OCepEUIIT2Wb/fCrr77i119/5c4772Tq1Kn06NGDHTt2cMMNN5RqPDVr1uSHH36gZs2aDB8+nNWrV5fq8wsjNTWVOXPm0KRJE26++Wajw3Epkqh4AJPJRHorMGFieLPhmDCR0UoVeiBW06ZN2bp1K2+99RbHjx/nX//6F88888w1Cyb9vKcvZy/XYxHDOEojFjGMs5fr8fOeu51VLCFcTn4tKn5+ftSuXfsf15ctW5bu3bsDkJiYyAsvvEBCQoLbL0YmhBB5MZvNmM1mAgICiI+Pp3v37qxbt47Ro0cTHR1NpUqVDIkrPDyc9evXExISwoABA9i8ebMhceRn8eLFxMfHM2HChCItdOkNJFHxAFprwsqG8WGPD5ncbjIf9PiAWmVrFWlKPl9fX6ZMmcKvv/5Ku3btePfdd2nTpg07d+4E4OO2w3m60vvsUF1AKXaoLjxd6X0+bjvcWcUSwuXkNevXkSNHqF+/fp4DQoODg9mwYUPO/vz58wkJCZFuYEIIjxQcHExISAhpaWkAJCRYx81+9dVXOR/0GKVZs2asXbsWPz8/7rrrLnbv3m1oPPaioqIoU6YMI0aMMDoUlyOJigdQSvF+9/e5pc4tAETWieT97u8XKytv1qwZ27Zt48033+TYsWPcdNNNTJkyhQ+qBdAo809MOgsAk87ihsw/+Kx2eYeWRQhXlrtFJSsri2PHjuXZ7Qtg//791zTj+/j4EBkZyb59+5wfrBBClLL9+/cTGRmZ8/7DZDK5VJ134403snLlSjIyMujduzcnT540OiQOHDjAtm3bGDJkCBUrVjQ6HJcjiYr4B19fX/7973+zZ88e2rZtyzvvvMNL3e8gdO8GLJjwtViwYKJNyp/cOPOU0eEKUWpyJyoxMTGkp6fnOZAeICIigkGDBqGUwt/fH4vFQv/+/YmIiCi1mIUQorRERETQr18/wNoCrbV2uTqvW7duLF68mKtXr/Lcc89x/PhxQ+OZNm0agKxEnw9JVES+IiIi2L59O6+//jonrp5h4b9nkfS/Txl8JBkFbOYWzv3rhNFhClFqcicqBQ2kt5kzZw5KKSZOnIhSitmzZzs9TiGEMIo71Hl33nkn33zzDfHx8dx6662cO3fOkDiSkpL46quvaNWqFR07djQkBlcniYookK+vL8899xzbxoRQrl59kr/9krlP9qfpf/6L3vF/3LDN8bN+xR3fyP/90Iy44xsdfm8hSiJ3olLQ1MRgHT8WHh7O0qVLefvtt1myZAnh4eFFGj8mhBDuwp3qvCFDhvD4449z9OhRevXqRXx8fKnHMH/+fBITE2UQfQF8jQ5AuIfNl29i453DWLJ/M1/vXcLGHcthx3Jq+pjo/Vsf+vTpw5133kmtWrWK/Yz01Ctkpiex6/cX8CuTxq7fX6Bbpfn4+pfFP7CC4wojRDHlnp74eomKUoqFCxfm7Pfr1y+nW4QQQngad6vzbr/9dsLCwnj66ae5/fbb+eGHHyhXrlypPT8qKorg4GDuvffeUnumu5FERRRKu1MppNXIoEV4IBOrd8SUUYXjcWn8tH8ta9asYdWqVQC0bduWPn36cNddd9G2bdtCT5EMsHlbO5SC+MAafM6bPBz4Mdt3d0Nr6Bl51FlFE6LQ8ur6ZTKZqFevnoFRCSGEKK6nnnqKK1eu8Prrr9O3b19Wr15NmTJlnP7cXbt2sXv3bsaNG0f58jIxUX6k65colPZfTuXU2f/jbMphwAeL30VuquHHlp3biYuL49tvv+Xee+/l+PHj/Pe//6V9+/aEhYUxbtw4VqxYQUpKynWfsWN3vzzXatmxu5/TyydEYeSenvjIkSPUqVMn57gQQgj38+qrr/LII4+wceNGhgwZcs0U9M4ig+gLRxIVUSjvLVhAjL8CFLpsQ0AR42/ivQULqFixIkOHDmXu3LnExsayadMmnnrqKUJCQpgxYwZ9+/alcuXK3HnnnURFRXHmzJk8nzG1/f15rtXyefv7S7WsQuTHvkVFa83Ro0cLHEgvhBDC9Sml+OSTTxg+fDgrVqxg9OjRWCwWpz3v6tWrzJs3j/bt29OmTRunPccTSKIiCuX+ESMxmQLJqNaKpNrlyKjaCpMpkPtHjLzmOl9fX26++Wbeffdd/vjjDw4fPsyHH37ITTfdxLp163jooYeoXbs2bdq04cUXX+Tnn3/OqQweWnOZJhl/XbNWS9OMP5mwpvQHuAmRF/tEJTY2lqSkpHzHpwghhHAfJpOJL7/8kr59+zJ37lwee+wxp00CMHfuXFJSUqQ1pRAkURGFsuTfj2EONuPrl0LkuvX4+qdgDjaz5N+PFfi6hg0bMnHiRDZs2EBcXBwLFixg+PDhnDp1ildffZWOHTsSGhrKmDFjqHFkK80u77tmrZaWqX8wwvTPFb9F6UqM+Ytt399GYsxfRodiKPtExTaQXlpUhBDCM/j6+jJ//nwiIyP5/PPPef755x3+DK01UVFRhISEMGTIEIff39PIYHpRKLc8+SqrF8+jza97qRwfT/MDv/FLm9b0euCJQt+jQoUKDB48mMGDB5OVlcX27dtZuXIlq1atYtasWcwC1Cp//Fpv4qbwdvwaXpWVpgB+yXiKz/fOICQkhAoVKlC+fHl8fCR5KQ2ZaalYMtLYtP1TPjjUh8mpn9Lrttcw+QXgGxBodHilzj5Rsa2hIi0qQgjhOcqUKcOyZcvo2bMnb7zxBiEhITzzzDMOu//27ds5cOAAjz32GMHBwQ67r6eSREUUStCo+xho11+z7unT1D19Gr5fDQd/K/L9fHx86NKlC126dOHtt9/m6NGjdH2pK+n7anFlz8/8+Ms2AK4Cx4E2H1zbh7Ns2bI5iUtISEiB3+d1rmzZsjJneSH8sKkNmfiw8lxvHmu8gJVnb4Rt/8JXW7it9x9Gh1fqTCYTJpOJ9PT0605NLIQQwj2VLVuW1atX0717d6ZMmUJISAjjx493yL2joqIAHHY/TyeJiiiUOkuW8suokVS9ehUTYAFiQ0LoMHuOQ+7foEEDNryxmvTPz+Pb08Ly08uprutyPu0yKe1NBPqU4+rVq1y9epUrV65c8/2xY8e4evUqqamphX6eyWSifPny+Pv7U7NmzUInOPbnAgMDPT7ZmbTpVfxUOiObr+G9io9yT8YKntv8HJnan9t6Gx2dMfz9/a9pUalfv77BEQkhhHC0SpUqsW7dOrp06cJDDz1E+fLlGTZsWInuGR8fz8KFC+nSpQsREREOitSzSaIiCmX2x2eoEVKLalevYlG+KJ1JfEgtZn9yhke+aOyQZ8yfvoSmZcM4n74Xf58YEmuZaFCpNX+mnub516/fTzQ9PT0ngckrocnr3OnTp0lKSiImJoarV68WaUpCPz+/YiU49vuuPq3t+7e8SJbyYRbjOUojtlTryOu3vImPxQIMKpUYlFLrgBpY8+NE4HGt9a95XDcG+DfWsXf/BzystXb4HJN+fn45Y1Rq1qwpTfdCGMTV6gbheWrUqMEPP/xAly5dGDFiBOXKleOuu+4q9v3mzJlDWlqaDKIvAklURKEMfbEjpwa/DyhOh/WgzukfqJ+aSvcXOjrsGeWP7CUm+2+MArJiThMTc5ryhWy18Pf3p2rVqlStWrXQz4yOjqZ3b2vTgNYas9lcqAQn9/cXL17k6NGjJCQkFGlKw8DAwEIlOPklO84erzNGzUWrv+fc2E5XttMVZbJwzmlP/YfBWusrAEqp/sBsoJX9BUqpcOBVoC1wAVgOjAOmOjoYW6Jy9OhRmjZt6ujbCyEKz6XqBuGZ6tWrx/r167n55psZNGgQa9asoXv37kW+j9aaadOmUblyZQYMGOD4QD2UJCqiUCrVDOZ8ndocqHEr8TVac7V8OM0D/qRSTcd9mnzfGx8z99Wn0SlpKEADKiiAe19412HPKIhSisDAQAIDA6lZs2ax7qG1JikpqVjJztmzZ/n9999JSkoq0jPLlStXrNacuLg4EhMTCxyv88W2S7zfOpWjQbWwKB9MOouGKTFM3hsItxTrR1Rktjci2UKw/tfIbSCwQmt9HkApFQU8h5MSlbi4OC5duiTjU4QwkKvVDcJzNW3alOjoaHr06EGfPn34v//7P9q3b1+ke2zatIm//vqLJ598kjJlyjgpUs8jiYooFKUUB5s/yKXTibS+pTZ712sO1u5ChAPHaCx55wjpQbXxSzmS89cmPag2S949wiNfuMc4AKUU5cqVo1y5ctSuXbtY98jKyiIhIaFYyc6xY8e4cuUKZrO50M+zjdfJK8Ep++dqMn5rTFKF9piCgwnocgsdA3ZwR8b3wL5ila84lFJfAT2yd+/I45I6wEm7/RPZxxzOz8+PQ4cOATKQXgijuVLdIDxb27ZtWbVqFb169eK2225j8+bNRRpnYhtEP27cOGeF6JEMT1SUUieANMA2EvpNrfUCpVQjYA5QBevkT6O01gezX5PvOeEcWmvKVynDjXfUo37rqtSoH8Khn8+jtXbYgPLUK5/i62ud/tXWouJ79RTmrM+AXg55hjvw8fGhYsWKVKxYsdj3sI3XKSjZ2bdvHxUqVPjHuZiYGK5cuUJmZqb1Zr9sA6yzsFVu0pIf63bGPPxOSnN0jdZ6BIBSaiTwNnm/ISkUpdRkYLJtPzg4mOjo6DyvNZvN/ziXmZlJeno6AElJSfm+1t3lVXZvIWV3n7K7Ut3gaaSMefvPf/7Df//7X26++Wbef//9QvXAuHz5MkuWLKF169YcP36c48ePFzfkYnHnf0vDE5VsQ7TWe3MdmwZM11rPVkoNxNr3tH0hzgknUEpx27gWOfv1W1elfuvCjwUpjPmhgxhwdiVlMzJpd+I8u+vVIEkFsiS0P0869EmerzDjdezH5+SmtSY1NZXnn3+ebWHN6HgykXrH92G+5McP5QN4f0s0/33lX84KP19a6zlKqSilVGWt9SW7U6cA+5UX62Ufy+seHwAf2PbDwsJ0fj+HvH5G5cuX5/z58wAMGDCAG2+8sRglcX0F/f/wdFJ29yu7K9QNnkbKmLfevXvTuHFjhg4dyquvvspPP/1EaGhoga956623yMzM5LnnnjPkZ+rO/5YuuTK9UqoacCPwTfahJUBtpVTDgs6VfqTCkeZOvIOAoHY0O3uRCilpNDt7kYCgdsydeLvRoXkdpRRBQUEMvLUL8+KTGZCxH0tgHP86tYoZB/6gd+fWpRVHBaVUqN1+P+ASEJ/r0iXA3UqpGsraxDcBmO+MmOxnapNV6YUwhivWDcJ7DBo0iOnTp3Ps2DFuvfVWLl26lO+1FouF6dOnU716dfr27VuKUXoGV2lR+Sq7AvkZ6xSCtYFzWutMAK21Vkqdwtqv9GoB544YEr1wCMtt3bjdbsas0CvJhG6bg+X2b4q1qKQoue2zv+Tv8amKXy7F8sul9XBYcdPt95RGCCHAIqVUINYpSOOAu7J/7/+HdZDsCq31MaXUS8DW7Nf9iLXl1eFsq9NXqlSpRN3zhBAl4nJ1g/AuY8aM4erVqzz55JPcdtttbNiwgfLly//juvXr13P8+HGee+45l1+SwBW5QqJys9b6lFLKD3gN69iTF0p6U+lrWjiuVHbfp59CT/sfNa5cxgdNFopzFSpiGv8gp5wUoyuVv7QVpux1bh3G8Z9WYTIn5Iwb0mXKU6/rXaXyc9NanwQ65HPuwVz7M4AZzo7JlqhIa4oQxnHFukF4n8mTJ3PlyhVeffVV7r77btasWUNgYOA110RFRaGUYuzYsQZF6d4MT1S01qeyv2YopT4CDgGngZpKKV+tdWZ2a0sdrP1KEwo4Z39f6WtaCK5Udq010/adJzT6a5SvH6bMTI526sP4UaOctgK8K5W/tBWm7MnmdL74ZRPKnABYJzmgTDnuGD6Q4DLe+cmQLVGRGb+EEEL897//5cqVK3z66acMHjyYpUuX5vydiImJYeXKldx2223Uq1fP2EDdlKFjVJRSwUqpCnaHhgG/aq1jgT3A8OzjA4AzWusjBZ0rpbCFkyil6H16N0opKo0cYd0/s8tpSYq4vtYvrSEr0TrDla0DWFZiOm1eWmtcUAaTREUIIYSNUoqPPvqIESNGsGrVKkaOHElWVhYAM2fOJCsrS1aiLwGjW1SqA0uUUj5YP6w9BozIPjcemK2Ueg5rK8pou9cVdE64Ka01frVrU+XhhygXGUlg69ZcXfW9Q6dAFkUz4YS1t4S/3UxsEMf4kzOAuw2NzSgZGRkAlC1b1uBIhBBCuAKTycTMmTNJTEzk22+/JSQkhCeeeILXXnuNGjVqcMcdxZ412+sZmqhorY8BbfI59xfQuajnhPtSShH20Yc5++V69qRcz54GRiRGvPEhCz/5hMY7d1AhJY2mZy/xV8fODH78MaNDK3Vms5n09HSOHTsGwLp165gwYQL+/v6yyrAQQng5X19fvv32W+644w6ioqJYunQpGRkZVKlShZSUFPlbUUxGt6gIIVzYxT79uMVuJrZaV5KoFb2eiz/8H1W9bCa24OBgLHY/iw0bNhASEoKPj8/fi2MKIYTwWgEBAfz4448AxMbGAvDbb7/J34oScMl1VIQQriF82XfEVgjB9vbcAsRWqED4d0uNDMsQ+/fvJzIyEpPJWm36+PgQGRnJvn37DI5MCCGEq9i/fz8333xzzr78rSgZSVSEEPkKaNQIFdEYk1IoPz/r14jGBDRqZHRopS4iIoJ+/fqhtSYgIACLxUL//v2JiIgwOjQhhBAuIiIigkGDBqGUkr8VDiCJihAiX0opal9JhuyZ2FCK2leTvHZygzlz5qCU4oknnkApxezZs40OSQghhIuRvxWOI4mKECJftpnYwj79hGpPPUXYJx/jF1YbrfX1X+xhtNaEh4ezdOlS3n77bZYsWUJ4eLhX/iyEEELkTf5WOJYMphdC5EtmYvubUoqFCxfm7Pfr149+/foZF5AQQgiXI38rHEtaVIQQQgghhBAuRxIVIYQQQgghhMuRREUIIYQQQgjhciRREUIIIYQQQrgcSVSEEEIIIYQQLkcSFSGEEEIIIYTLkURFCCGEEEII4XIkURFCCCGEEEK4HElUhBBCCCGEEC5Haa2NjqFUKKXSgLh8TpcFkkoxHFfizWUH7y6/s8peVWsd4IT7OoXUDfmSsnsnZ5Zd6gb3ImX0HK5czgLrBa9JVAqilDqjtQ4zOg4jeHPZwbvL781lLyxv/hlJ2aXsIn/e8HOSMnoOdy6ndP0SQgghhBBCuBxJVIQQQgghhBAuRxIVqw/+v737j72qruM4/nwlldNZVMtALb8qrrWQviK41pJB2sRcG4VGDtuglaPNWXO0Fi1Hzf5orenQaZC0MpkjXcZaRQZKmliIfPkhzYDkO7RpSoxRYxXouz/O58bh8r33e7/3+73fc+49r8d29j338znnfD/vzzm8uZ97Pud+i25AgaocO1Q7/irH3qoq95Fjr6Yqxz4SVegnx9g7ujZOP6NiZmZmZmal4zsqZmZmZmZWOh6omJmZmZlZ6VR+oCLpYkmbJe2R9IykDxbdpk6QdLqkX6Q4d0j6naQpqe5sSesl7ZX0nKRZRbe3UyQtlhSS5qXXPR+7pLdKujvFuEvSA6m8Etd+O6rUN84NGecG54YaSYOS/iJpe1oWpPKG/VL2PpO0IsUVkvpz5W3FVNZ4m8Q55DlNdV0VZ7s5u926wkVEpRfgMWBRWr8OeKboNnUoztOBT3DiuaSbgU1p/UfA8rQ+E3gJeHPRbe5AH/QBm4GngXlViR24A7grd+4npZ+VuPbb7LPK9I1zg3ODc8Mp/TII9A9R3rBfyt5nwCzgvPrY2o2prPE2iXPIc9qNcbabs9utK3opvAEFn+yzgSPAhPRawCvAlKLbNg6xzwAG0/q/av9BpddbgKuKbuMYx/smYANwGbCJE29Gejp24Mx0jb+trryy134LfVbpvnFucG6o8vWfYj7lTW2zfummPsvH1m5M3RBv/Tkc6pz2ynltNWe3W1f0UvWpX+8FXo6I4wCRnZ0DwPsKbdX4+DKwTtK7yEbNr+TqBum9PrgVeCoinq0VVCT2i4BDwDJJWyU9KelKqn3tD6fqfePcUI3YnRuauz9Nh1st6d0075du7bN2Y+rWeOvPKfRGnMPm7HbrOtrqFlV9oFJJkpaRfVrw9aLbMh4kTQXmA7cX3ZYCTADOB/4cETOAW4C1qdzsJM4NleLc0NisiJgGTAcOAj8puD02ej15TquQs6s+UHkRmCxpAoAkkY0gDxTaqg6StBT4NHBNRByNiH8AxyVNym3WR2/1wRVkMe2VNAh8GFgFfIbej/0A8AawBiAiBoD9ZG9QKnXtj0Dl8gI4Nzg3ODfURMSB9PMYcCfZddIsL3Rrzmg3pq6Lt8E5hS6OcyQ5u926Tra/VZUeqETEq8A24MZUNB94KSL2FdeqzpF0K3AD8PGIOJyreghYkraZCZwL/H7cG9ghEXFvREyOiL6I6AP+CNwUEffS+7EfBDYCVwNIugC4AHiKCl37I1G1vADODc4Nzg01ks6UNDFXdAMw0CwvdGvOaDembou30TmF9vtgXBreRJs5u926YhX9kEzRC/B+sm962QNsBS4puk0divM8IIC/AtvT8qdU9x7gUWAvsBuYU3R7O9wXmzjxwGzPxw5cCDwO7AJ2APNTeSWu/Tb7rDJ949xwUl84N0S1rv8GfTIA7Ez9sg7oG65fyt5nwEqyb3I6Dvwd2DeamMoa71BxNjun3Rhnuzm73bqil9pXm5mZmZmZmZVGpad+mZmZmZlZOXmgYmZmZmZmpeOBipmZmZmZlY4HKmZmZmZmVjoeqJiZmZmZWel4oGJmZmZmY07SoKT+urJNkua1sO9ySa9J2p5bvtHCfrMlzW2/1VYmHqhYz5M0Nf3VaTOz/3NuMCu9NRHRn1u+08I+s4GGA5XaX5q37uCBipmZmZl1DUk/lrRS0kZJeyT9XNJb0t2bJcDCdAfmNkl9kg5L+q6kbcDNkqZI2iBpZ9puXu7YIel2SQPp2AtT+VJJq3LbTZR0UNI7xzn8SvFAxcZc+ke+TNIWSfslLW5hn8UpWeyQtFVSXyr/XEokOyX9StK5qXxRSjIPStqV9rkwd7zlkvZKehb4bKdiNbPWOTeYVdLa/PQtYMYI9l2ok6d+LcjV9QOfBD5A9pfV50fEduAHnLgT8+207duB3RExPSLuBNYAD0XENOB6YLWk83PHjoi4lOzOzF0p79wHzJM0MW2zGFgXEYdGEI+NkAcq1in/iYjLgWuAFc1utUqaDdwGXBMRHwJmAa9Kmgp8L5VPAzaTJYqamcCyiLgE2AB8LR3vWrLEcxlZQuwb08jMbDScG8yqZUF++hawdQT71k/9WpureyQijkbE68AW4KImxzkGPAAg6SxgOrAaICL2An8Arshtf1+qewF4ApgVEYeBh4HPSxLwJeDuEcRibfBAxTplDUBEPA8cByY12fZa4KcR8XLa52hEHAXmAOsj4m9pu3uAj0k6Lb1+OiL219Y5kaSuBH4WEUciIoCVYxWUmY2ac4OZjYV/59ZfB5o9e3I0It5oUh/D/K5a/QqyqWVzgdciYmDYVtqoeKBinTKSBNKq+kTS6u8YLgGZ2fhxbjAzACR9StL9Y3zYI2RTvYYUEf8EtpFN3ULSFOCjZHdOamp1fWR3Wp5M+z4PvACswndTxoUHKlYGvwRulDQZQNIZks4AHgfmSjonbbcE2Jhu8zazAbhe0lnp9uxNnWq4mXWUc4NZb7uYbGDRSP0zKne0cMxHgP7aw/SNjgsskLSDbDrXFyLiQK7+NEkDwKPALRExmKv7IdmHHw+30BYbJX9FmxUuIp6Q9C3gt5IC+C9wXUQ8J+mrwPrsPQUvAl9s4Xi/lnQ52ScmR4DfdK71ZtYpzg1m3S0i+oYom517+RHgKw32XQ4sb1C3qO710tz6fuDSul0m1m2/D7hqyEZnvh8R32xQNwe4JyKONdnfxoiyabpmZmZmZtWWPhR5R3p4Pl9+DvAYcAi4Ok0hsw7zQMXMzMzMzErHU79s3EjayqnX3O6IWFhEe8ysHJwbzMxsKL6jYmZmZmZmpeNv/TIzMzMzs9LxQMXMzMzMzErHAxUzMzMzMysdD1TMzMzMzKx0PFAxMzMzM7PS8UDFzMzMzMxK53+9pHHwN0BDxAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(12, 5), dpi=80)\n", "\n", "ax1 = plt.subplot(1, 3, 1)\n", "plt.plot(n_cond_arr, SI, '*')\n", "plt.plot(n_cond_arr, H, 'k-')\n", "plt.grid()\n", "plt.xlabel('n_cond')\n", "plt.ylabel('H, Entropy')\n", "\n", "ax2 = plt.subplot(1, 3, 2)\n", "plt.plot(n_cond_arr, t, 'k-*')\n", "plt.grid()\n", "plt.xlabel('n_cond')\n", "plt.ylabel('Simulation time (s)')\n", "\n", "ax3 = plt.subplot(1, 3, 3)\n", "plt.plot(H, t, 'k-*')\n", "plt.grid()\n", "plt.xlabel('H, Entropy')\n", "plt.ylabel('Simulation time (s)')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.9.12" } }, "nbformat": 4, "nbformat_minor": 4 }