{ "cells": [ { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "#AUTHOR: Lisa Blum Moyse\n", "# lisa.blum-moyse@inria.fr\n", "#\n", "# REFERENCE: Blum Moyse & Berry. Modelling the modulation of cortical Up-Down state switching by astrocytes\n", "#\n", "# LICENSE: CC0 1.0 Universal" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import random as rd" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA7kAAAEnCAYAAABsVITaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAADK4klEQVR4nOydd3gUxRvHv3OX3hPSSAgQeu/SQVAUKXZQwQJWLPxUEKxYsIuKghVRFBRsCCpNUHovoQUINYSE9N7b3c3vj8nc7d3t3u3lLg3m8zz33N2W2dnd2dl533kLoZRCIBAIBAKBQCAQCASCKwFNQ1dAIBAIBAKBQCAQCAQCVyGEXIFAIBAIBAKBQCAQXDEIIVcgEAgEAoFAIBAIBFcMQsgVCAQCgUAgEAgEAsEVgxByBQKBQCAQCAQCgUBwxeDW0BUQCAQCgUAgEAgEgqZAXFxca61W+5hGoxlDKQ1u6PpcSRBC8g0Gwwa9Xv9N3759k5wqS6QQEggEAoFAIBAIBALbxMXFtXZ3d18VERERFBQUVOzh4VFNCGnoal0RUEpRVVXlXlBQ4J+ZmVlQXV19hzOCrjBXFggEAoFAIBAIBAI7aLXaxyIiIoIiIiLyPD09hYDrQggh8PT0rI6IiMiLiIgI0mq1jzlTnhByBQKBQCAQCAQCgcAOGo1mTFBQUHFD1+NKJygoqFij0Yxxpgwh5AoEAoFAIBAIBAKBHSilwR4eHtUNXY8rHQ8Pj2pn/Z1F4KkmSmhoKG3dunVDV0MgEAgEAoFAILgqmDdvHhISElo1dD2uBnJzc9GvXz+7waPi4uJyKKVhlsuFkNtEad26NQ4dOtTQ1RAIBAKBQCAQCK4KEhIS0Llz54auxlUBIUSVrEMIuSS3vEkLuYSQDgAGAYgCEAbAC0AugGwACQB2U0rLGq6GAoFAIBAIBAKBQCCoT5qckEsIGQTgMQA3AQi3s7mOEHIYwHIAP1JKC+u6fgKBQCC4sjh9GjAYgC5dGromAoFAIBAI1NBkhFxCyH0AngfQFYA0XncJ2OxtHoByACE1n1AA7gAGAOgP4H1CyM8A3qSUptRj1QUCgUDQhOGWaSKtvEAgEAgETYNGL+QSQkYA+AhAbzDhNg/AHwB2ANhPKT2vsJ8fgH5gQu4tYGbNDwO4lxCyAMC7lFIRAlwgEAgEAoFAIBAIriAavZALYEvN90YAXwNYTym1G7qbUloCYFvN5wNCSGsADwD4H9iMcBmAt1xfXYFAIBAIBAKBQCBoWlRXA5cuAbGxgFbb0LVxjqaQJ3cjgEGU0jGU0r/UCLhyUEqTKKVvAmgF4CWwGWGBQCAQCAQCgUAguOpJSwMKCoDcXPvbbt26FYQQzJ49GwcOHMCtt96KkJAQEEJw8uTJOq+rPRq9kFsj3O53YXlllNJ5lNIvXFWmQCAQCAQCgUAgEFwtHD58GABw4sQJDBs2DBqNBtOmTcPkyZPRqVOnBq5d0zBXFgiaNDodM/kgxP629cG+fcCgQcDZs0D79g1dG4FAIBAIBAJBY4CPVdUEWuRC7q5du7B9+3YMHDiwDmvmOE1SyCWEtKz5ma7GfJkQEg1ASylNrtuaCQTmlJYCfn7Aa68Bc+c2dG0Yy5ez73/+EUKuQCAQCKyhFHjvPeChh4DIyIaujUDQNHj2WeDo0YauhW169QI+/dQ1ZXEhd8GCBY1OwAWagLmyAkkALgLYTQhprmL7QwAS67RGAoEMRUXs+5tvGrYeUtzd2XdVVcPWQyAQCASNk7g44JVXgPvua+iaCASC+kSt1WFpaSnOnj2L8PBwTJ06VXab6667DpMnT3Zd5RykSc7k1kDAUgTtI4TcQik9pmJ7gaBe0dSokQyGhq2HFCVTlPx8wMMD8PWt/zoJrk7uvBO44w7g3nvNl+/bx9pinz4NUy+BQEpFBXD4MDB4cEPXpP7Q6dj35s0NW4+rgcxMICKioWshcAWumiFtDNgzVz527BgMBgPGjRsHjUZ+zvTIkSOYM2dOHdROHU11JhcAigAkAIgBsJMQcnMD10fQSNHpgDlzWLS4+sYRIffYMWbanJZWt3VSEnJDQoC2bev22AKBlFWr5GeKBg0C+vat//q4GkoBvb6hayFwliefBIYMARKboD3YiRPA9987vl9hoevrIrDm4EFmDr5sWUPXRCBwDG6qPGDAANn1Fy5cQEFBAfo24Mu8KQu5ZQAGgeXR9QOwihAyo2GrJGgs5OcDnTsDJ08Cf/4JvPMOMHOm/f3WrwdatgQqK507flwcMGuWSaBUI+QuXMh8eNevd+7Y9pCaoqxaBWRnm/5nZtbtsQWCrCygdWvg1KmGrknd8/nngJub+TMmaHocOcK+G0JR6izduzO/WkcJDXV9XQTW8Cwr9T1jfuAAGwucPq1+H0qB+fOBnJy6q5eg4VEbeIoLuf369ZNdHxcXB0IIevfu7crqOURTFnJBKS0CMBrAdwC0AD4ihHxNCGni6YsFzvLPP6zzfvtt00xKSYn9/Z55BkhJYYmwnWHoUODjj4HycvZfzWwOn/VVE9HOFRQUMHPRsWPrpny9Hnj6aeevpeDKYs0a1iY+/riha+JaysrYrJmUH35g3+IZEDQ1AgPZd2xsw9bjSqehsi7wybfp09Xvc+AA8NxzwIMPOnYsvZ5Z06nJuypoOhw+fBgeHh7o3r277Pq4uDi0a9cOgbwzaQCatJALAJRSPaX0UQAv1Sx6FMA/hJCGu6qCBkeqidLWqDzUCJoeHuzb2aBMXGDlx7QluO7axUyqeZ23bHHu2GqprolLfvFi3ZR/8CDw2WfW/pZNCUKAp55q6FrUP5s2MWuEuoA/j9zn70ph0iQ2a1ZWZlrmSCoGgcAV5OYy/2FnaSwp7xoTy5cDf/2lvP6vv4AdO+qvPs7A5Y6WLW1vJ4VbuDlq0bB+PbOm+9//HNtP0HiprKzEqVOn0L17d3jwgbMFcXFx6NPAgTWavJDLoZR+AOAuABUArgOwlxAidJBXIZ9/zgacABNWHRFy3dzUb2sLfsxqOwmuDh4Ehg1jWk4+qPjlF8eP9913wIUL6ra1HLzU1QCcH6epR3H+8suGrkH9M3o0oGCB5DSOPI9NiZ072bfU1UEIClcGTUlZ4Sp/dn6ujSloYkNz333Abbcpr7/tNuDaa2tXdm3a1vHjbAxRG+66i30ruFPKwuvoaL/G+/rSUsf2cyX79wPbtjXc8ZsSavq7EydOoLq62qa/7ZEjRxrUHxe4goRcAKCU/gFgJIAsAJ3AIi8Pa9haCeobqbZw9WrToNrWy/qFF4AOHVw3mOHH5AKeUnncB/bECfMXhyMBTigFHnkE6N/f8XrK1W3cOJOfkMAxkpOBhISGrkXtMRiA99+v22M0JSE3JYVF1lWDrSBzdSUcUQosWgTk5dVN+QJGXQu52dkmv19nOXfONeXUVqAR1I7aXOeePWv/3q8NfBbb0bo2hiwTAwcCI0fa3+74cTZuVEtVFXDmTOO1TKqqYjFqXE3fvn1BKcWiRYtk11+8eBF5eXlCyHU1lNIDAAYCOAkgDMAmACENWimBSzl9Gjh/Xv32lqbDcsybxwYHSoOZ3Fzmu2IrINWGDUxAlO6rdiBPqameAHDTTebrU1KUXxD8eGoHuvYGbOvXA088oa6s+qBjx6Yzm9qqFdClS0PXovasXQu89JL97ZSYPx+IibG9TVMRcg0GZsrHrULsIRdkrq6Fo2PHgMcfB6ZOZcFgHntMvVAuUE9dC3q9ejmXLuvgQetBtpeX+eBWrYAxdy4LglifAkllJfDaa+am/gJ5ystZe/ziC+fK4ddarelxZSXwySe1O1ZTUpT07MnS2qklPR0oLm68isazZ5mVX30rGOJq/J2EuXIdQCm9BGAwmIDrCcC9YWskcCWdOwPt28uvk9Om1WZQbTkofekl9lL5+WflfcaPZwKiXm96cfCORWmQKw02JX0RSF/258+zwfa778qXIXde+fn2TaVt1a2+tJKU2jdhOnv26vSLbQjstRl7PPcccPmyum0bu+knr9+ff6rbXi5wnCPR1WsDF2izsoAXXwQWLwZWrKibYzUUlDrfLl1FXbVZZ9LGffstm8274Qbz5ZWV5r65atvgG2+wAIz8XF1xzrfeCkRFKa9ftAh46y3ggw+cP1ZTw9Hry4Wpd95x7rhbt7LvlSvVbd/Y+2tnWbfOtim6EnKKzOpq++PNEyeYkrKu4daEhw+rv4euUM7GxcWhTZs2CAoKqn0hLqCpCrlzAdiMzUkpLQYwDsBXAC4BSK6HegkamCFDzP/36qUs5JaWMjNl6czH0aPy5fJ99XqWP9AyiioAuNeoUg4csN5PqbOQDoKlQi7fvrqazc4AwH//yZchN3gJCVEO+MSP8+uv8uul9XYGfhxbPkPz57PcwOnppmUGAzBtGrB3L5sZ57zwgroI2YLao5DP3aXccw/7buyDJkfrZ2smt6780qWDEf67sc+QS9m40b5rxrRppoCAteHaa1nf4Qz2ZqLKy5mioSFYsoR9y/kbFhWZfjvaLlz5fP79t3kfbwl/B19tM7lffMGeAUD9bKfceKY2wSP5bOX996vbXk5511gpLweef96x9jR+vO2gYkrICYTHjsmPETl6PWvzapR35887FxxU+k539JnOzKy969p7772HC2oDxdQhTVLIpZTOpZTaTUBRE3n5KUppLKVUBKG6CpAKmAATWv/+m/22FAbfe4+ZKSu4FJgh9SkZNYpFUbWEC7nSwY69gYXSy4J3RkuWmDSuSihp6H//3fYxU1OVy3TFzIl0FkHpOvz2G/uWplhJTwe++QYYPNg8Z/C8eeyeybF7Nzuv3budq/PVjrYek6/Vp5D733+Oa80drZ+c3xl/1m69FUhKYu3aFYwYAUhTD6amNg6/N0e56SagbVvb2yxerLyuqIhdY1sKux07WN/hCpTaxLXXAhERjpdXVsb8+J3BVj5bqQJ37lz2nZTErtn27bbLddVM7p13Kq/Lz2euOA0R2GvXroaP9jt9uuOBJuWuVZs2dX/tXCHk1tf9XbAA+PDD+klTp9R2bY2hkpLUl19Q4HjqpcpKlh2hrMw5hYTBYEqD2VRpkkKu4OqisJD5CNW2g+S+K5b78wGAXGck3dZgYCZhABPWDh1iv6uq5E1tX3/d9FtNCiG59fy/1GxYqQzLQa2j10lue1fMBkn9enkdCwtZnlRb2OqUlXyiN21i3//9B0yeXH+a5rQ0YMaMhg86sWKFY4FrcnOZabFl269PIbc+ueEGZtWhhqoq4M03HZ9VsuX7X1jIBNNp01wzW7V9u7nVSWoqOwbQtIRcRykuBvz9Tc87nyio62Bp9vqT2kS4LSsDxoxhfvxKZGebghMqwU0s777bep203vwaceH2u+9sl+sqgWTVKtPvzEzzPrxNG8dS2LiSYcNYJgZ+nmVldROgxxESE4FPP2X3bf9++W2Uol7X9n498wyLsMxzeishPZ6jUYp5O1y3znWB0WzB38e2Yqg0JHVtsVBQwNpDTo68heDVhBByBY2ep59mHfHmzc6V48jgT9oZSGc8pcsHDWKmtpYcP276rdZcWUnIVSOs8fNqzGkueB0nTwZuuYVp73fskBcIbJ2zGnNaW37TrmbaNDYosWybZ8+yAEy2TPRcyb33Oha45plnmKk4t3Lg8BRazqKmDTbGdgowhdbrrwOOuhLJmQtL23J2Nvuuq/MuLmbf9oSiporBAMTHM5cFrkhsrG1IDb6+9nOqhocDkZG2t+HxH+RMuqXXh/+W8x2Xw5U+uZzISOausGgR63943ev73cVNvAHTu6lrV+bm05C0bcuUpoDp2xKlMYUzyq0DB4AHH7S9jauUZx06uKYcW7jKqiU+Hli2TN22rlasU2rubuAIjXksWN+4aEjDIISEgqXwaQXAh1L6povKfcAV5VBKVTZXQWOCa1ed1X45ovlUGqhKkZrjKm1nL/AUNwWxDDzFt7enbQesU9bY69jrY5bTMgjK/v1MEOOa1S1bWERYKVVVLIryxInK5Uo7761bWUqAhtRUco2xpbD+2WcsANPKlQ1vEifH8uXse8cOc3NCV83kWrZnOdQMQAYOBN5+m7kIuIKkJCA2lpkqWvrvc2rrP6s2yFRd56WeO5cFD2rMDBjg+HMhZwbOOXqU+Y517ep01WzSGAeNzz3HvuXqJreMB7lS205TUtQ9z2r580/rYG71PSh/+WXTb72e9XvchDQnh12b8HDHyqyvuiuNKbhpqr9/3R63KeAqIbdHD/Y9ZUrd3d/8fCA42Hp5erp1QDqdjinQY2MBb2/HjyU9h7w8IDBQ/p3f2H2uHcElM7mEEDdCyMcAUgD8AuADAK9bbBNMCMkjhFQQQpo7eIgfAHzv5GeJZaGCpoG0w+KDc4DNsDoSKMDebKqU9HS23NKH7sknlcuXK8eW2W9yMnD77fLreF2lgrRc/c+eNSVzVzJjsodcuc52cpZBsl5/nQl9fDZLLkjLggVMe/3VV8rl8rawbBlw/fXAjz8qb0sIMGsW+71pk+vNpJ5+GvjnH/ZbqW1ZLi8oqLsgREq0bAk8+qj8Oktfn9oIub/9ZvKt5si1wb//tq+QmDfPfJv9+4GHHmK/ly83meTWFj7jvsTibbB8OTBzJvtdW0HfUsgtLlaOmO4oGRn2UwPVhcm8Xs+iuvNZYldx4ID6gDfSutgaaO7Z41ydbMHv48CB5stPnzb5uqqBBxhUk6/TUdQIuXo9i8KttL0U6TPsaMTuigrWHi0tRZRQI+TOn2/+/rfF0qXKM6GA+TNu+dyEhdXOv9oRIej1163HLpbvXKV3sJKQ27s3EBCgvg6O4ozAaHkuBw4whfaIESxFoBr3KJ2OmdyrSZG2di37bgqC+YUL8m1H7jyLihz345drR+XlzDRerW/w+fP1Z5Xmalxlrvw7gGcBeIDlp7V63VJK8wGsqNnm1locg7jgI2gC3HGHaWALmF5CqanmaXQGD2b+SGpfLtIOLy7O5FsrB48AO22aurKV4HW3rGNmprUvltS/V+05yaWecMVMbm21lgaD/AvLsjwepIszaBAbyAO2BRledx5t0DJ4n+W58cATo0ebm0lFRzufC/izz0y/La85F2j4OXGCg+Vnqj/6iCksXMnOneyZSUkx+ZRb4uZmHljC0hx85EimUHj/fWDCBPky7r7b2h9Qej1WrWK+h5aCpVw7lcvDSCkzG7vvPvsmdVISE9l1lcOyPd53n+nYrhByk5LYgFNO8KrNs9W8uf3cjTxKq7NQaqrj6tXAK68As2fb3ic3VzkyvauwNZPbUIwc6disObdkUfJppJQ9r7VBLtWS5TP2/ffm26emKgcglLbT++4zD3TGKSqSV9p5e7M+XilQoCVqhNznnmP1UNr/6adN/6dOZW4kSkjdMn74wbnc0tnZTBnnyHP95pv209Xs2cMUw48/bv6e5+9Xy3vLo5Tz65CfzwRpOZNXSh3Ps6vm/Lp0MVmf/fqrsuJpwACWFnD7dvbek/OdnTePuTdxlixhqRzttam4OFMASnt1TkxUFhi//NL2vhxXzPDKpR+y1cdZKh3Pn1enxLecCFHrs1xQYDtQaWPGaSGXEHI3mNCaBaAfpbQHAKW0yDze63gHDxOr8GkDIAcAtbGNdFtBE2D1avOX8YYN7Hv6dPPOgHdOatPKSF8K/frZ9oeSPvxqfYHtmSvrdCZB1jJxOKXWAxC5WUAeGXPPHla2pRAlPaYj9awNmzax2VcpXbsyXzNLLM9FzveTz5rbqj8hrLPlg0SDgb0M+fV0RDnw9dfs95o1rFw+y1wblM5Pzk/OcnajtJQJEcOH1/74cgwfzrTktli6FOjUyfTfsm1s28ZMtV56CfjjD/XHlt7DO+9kM96WMxfSbWz5/1Fqem5Wr2ZB6DiVlfIWFwAzmZ892/y+8vP7/nvml9yunXU/4Aoh11ZanNoOing/WJdQyhQd3ASW94NFRSzFhuVMJmfQIHMhKDFRedBV29kV6Uzunj1Az57m610p+G7ezMrj/atl2enp7NqotcpISVE3aNZomOWFZZYANcgJuZZtTRpYiVKgRQv2USpPipwSIzAQuO465Trt26e8Tgq3onAGqdKRc+21zPw4M5MJslu3Misj6TP+1FMms9TaMGgQey7qYtYwJob5L3Orh+Rkk5uFUj/Cr0PLlkyQlhOmt22z3Q8VFFgrkG2dX3k5U04nJACPPMKW3XOPskuIJb6+1vnVX3jBPL4GfwfYG+9Jx1ZydX7jDVPf0batcuC3p54y/1+X5ujSiQF76YfkKChw3sqJ09BBNF2NK2ZyHwQTMmdTSo/Y2fZAzbYyCViUoZReUvgkAdDb2cb4qcW5CRoB9joXtWaAvBzLnGP2yrf0HVVCbpDFtaiUsoBLfn6mqHf2sOyg4+OZVhdgL4+33wYmTTLfZtculh/OEeQ0vceP23+ZjB4NPPsse/nwMk6fltcOWgoScteKH8+WRvLPP9mgjEcI3bqVvQw//5z9j4uzXWe5srnG3zLFjK12YZmD2HJbHsBEKb+nVHvMr4Wtl9R33zkWWIojtQ4oLZW/p8nJzJ9xwgTXRaNUM+CTbvPww8rbpaaaX99nnjH95gMaaURzDj9X6QBCWs7ChWwgZzmjUdvgW1Ih15agrHYwzNPj2DLfdzX8+vBZbekM27p1yhFfLZ+rtm2VA8zYO/8DB0xmqd26mZa//bYpqjJgHuBPidoqrrjyjp+vtL8qKACiophQr1awHjOGDZpXr1a3vWWcBY7BwFxmjh9nqYOkxzcYrIVcyxl4qRuJoxH/Ofn55rkzLdO2OROZ39WCxI4dbKY2MpK9f6+7jim/LMcAzriycGGwLoUgPsb54AOT4sXW8R5/3NT/7d1rvu6ff5RTmSUksDYVHMwUgBwfH6BvX+Xj9e6tHKgvJUVdGhq12QEcuc5yfc3cuer6DjVl1RbLfuPUKfP/SumH7PU3hw4xxZuS8s3etSstZcosHhDuSsAVQi4fetnV81NKywEUAghzwXEFVwlKDzafMVQ782IwsAe4jcWcvrMvpwUL2EBcrp533WU6Bp+J+fJL+7N2lMoPFqQmiXLmicOGWQ/ap01jLymAabPfecd6PzlBnpv6ZmfbTi/QujXT6DuCrc7633+V11m+DCzNzLgvjhKbN1u/AOzlKuasWGEq39JHjVKTiTtgapNKAz5pGiUlX+rt201C+yOPOJYiSI7QUOWgJN26sZnaG2+0XcbJk9YB1zhScyZHhVypFYMc9pQXcs+wvXugtK+t/iQ+ns04EMLaUsuWpjzWUiHXlja8sJCZFPKBX0oK0KyZ+cDz11+ZiTJgOw6Aq7G8b3UREMhe2xgwwGSWKhWmPvgAeOstdcfIyWEWE2qCB1FqO4fyb7+Z5/HmJvAJCeqivQMm33dbfqJSLPuk4mKmjHrlFab06NnT2p9ebibXMvaBVElgL6e43D3X6ZibkFT5YIkzOdbVtrOnnlKOZ2EJP2db7xU1KCkeOM4KQWoUJlIFnC0l9KJFpt8VFew54u5ZY8Yo5+bdudP8v78/E3zKy63N6KX36swZ5bq0bGk7mCTnlluYKbKlFQy3CuGCvrNCbm3JymLvYXsCe23yyqp9ZqTvFbm4JgC7F1ILP8t2pXRNyspMbVwuNWZTxRVCbiCAwhoBtr6OKbiKkA46pSYt/EFU2+kdPChvouzsA/3ss8wkSK1Wf/p062WWJtGUyndGXFgF1AVaiY9ng2fe8Sp1jEuXWi/j53PPPcwX0tJ8iaNmVlqpbGdxNFXNE08Anp7qyrL8f++9wM03ywc/MRiYYMLhwricSa7lb256bSmMjRjBzOotycpi1++//5jplZyvptyssDN+Z5xu3ZQ1+lKzRzWzOY4MVpTauq12xNdJI+6eP2/7OEVFtmeVe/RgJoQAi/ickmIK5KM2unLr1sykkJuvLl/OZqSnTTMJdffcI2+hotYk7exZc8FMLZb3xNkopZaDZsA1ObjtERYG3Koy8se337IcylKlIb8OhDB/c2nQFek1smx/Sm3aVg5lW9tzAgKYWaWtmWBK7QtiUqRKqcOHWbAlaV8udy7TpzNrHVu4UsiVc9sB2LPDIzTb60ecFW45UgXLhx/K14lz5gyzqDp3jgl4zlrIcLPv2rpSdOsGXHON/e0s23NJCQv4Jcfixcrl1KbvAVjU67Ztze/psWPMKmTOHOs6HjvGTLmVrHV4oDcenFOKPcWpJS+9xCyqfHyULUR0OnPFnBwGQ+3exZSaX1duDWbZp9hqa9XV7FmXS2dnOYFwpeAKgTMfQCAhxMvehoSQFgACwPx3mxSEkGaEkAcJIT8RQk4RQkoJIZWEkMuEkD8JIYp6RULIG4QQquLTTqmMq5EhQ9jLRNqpyUX5/P571nEozTBJkZo6ctT69Nri0CHbnYujMyFKQq6jEVotfY3sBZCRwrXGXCtYGw2lEq4SctUMvu3lolQzk8ujKAPywU8s6yEdjJ4+zV6yUhMgXvalS8zXUa4MjqVfLTeh/PRTZnplOfBdscJ2jtebblJepxZ7A1k17V3ufJX2q03uYy4oSAXD99+33k4q3AwaVHuBTirkqrFs4M+TVKCR+kzK8fTT6urXsSMTplNSmJCm9tl1ZCZ3717b0c0BeYuVuop4Wts+hZsuSgO/cYsNy2BpgHn9LYVRpXPjdVPr71abc6G09n79773HlGfcKgGQDzQjnSFUwhVCrk7HlKoajf0AgZZKB3tCeG2R3ls5lyDuxw6wOAcffsiEs5Ur1QfgUoIr6pzNY962re31cu1OyfTV1njLXiwIe0jvqa3jTJgA/PQTayuffMLep9J9uYJDzsddToGstk5K4zA1fZtcHBU1KL0b7AnVUrhwza1AdDr2vNR3xof6xBVCLjf0uVbFtjxWrYJnT6MmAywN0b0AOoNdu2oA0WCBt1YRQtYTQnyUi0A1gEwbnyvM5ds59uxR51/6wgtMw2nLZ8QWrppZsDUjXJu0PnUxGHRkAMBfePzFam+AphTJ1lbZ9YE9fyv+8lqzxtzcWfpSs8zraIncTDgv4+WX2UtWLnCQVPhQEvCkg2/pzMa6dfLb24uy64oovPZmBdQIuUr+nbXBlrmyPaQ+6bXRZh84wEwnpWbnPrbeAjW8+ir7ltbT3nVbtszc1N0ezz3HzG2feII9c1yIqaw09Xt//20yEVQScleutC578GDggQfMl9VWufHLL2zA6goFjCM895wpAJ3ce0AuRZ0tIVepj6yNkCs38LRlFurM+4IHlZL2lQcP2t7HcsaM40zgGl7mrFkmn1B7grXleXfuXPvjO3IcR5g71/6zYevetm7NBDi5GWRHsBUQD5Cf9VY67+XL2ey/3P12duZaesy331bejj9XiYkseNlNN6k3o1Y6nlrUjBvlFD5qlUD2LMs4tgTU0lL5/fh1q65mkzyZmerqlJdXt77ndYErhNyfwdLzvGVLwCOE3AXgBbDAU3b0v40SN7DAWU8CaEsp9aaU+oFFbq4Jmo4xAGx1yXsopZE2Pkl1egaNFJ3OtqN7XWuZnMlfqZbazOTWhVmfWh8ywOS7K/VtnDGD+QjKCamOzBKrffnYQ03aEnsvMH5vPvuMmSPL7WdPKLflC8xnWt980/qYSqbTSuh09ttSfSgQjh1jwVuUgmGdP28/uIcrLChszTTWpyLlzz9NwrLB4NigSfpM2ru399zDfFM5lhFJleBKGD6j5OVlsiC49VZThG3p8UtK5PuLigplUz815y3Xr02axHzCpQoYexYYgPMDrvnzTYN0te1Feo6WihR7Qq5a1xhCavfOqC3cjP+VV0zL7N1L6XrpfXM28FRlpXXUflt9SX0Nui2PYxm8yh72rClsuR/p9UzQrWss853zY8tRUsJm/y3TAboCtf0nf67GjGHfp0/bt4aR46ef7G9jef/VCIaWWTQcwZl9OWVl5goHpfgfaklMlA9U2phxhZC7DEz46wtgLyHkcQDuAEAI6UcIeYgQsglMGHYDsIVSaic8TKPkOkrpAErpV5RSoz6MUppEKX0EJuH2PkJITMNUsWnywAMsml9D4UozXFfi6plcQhzLZckHZHwgl5HBTGRdEe1VLh9qXWHrOmZkKGudpf6E9mY11BxbKtiXlrJZQyU/ZyXFzv/+Z/t8cnPrR7gbOpQFcFMKhtW/v3WKl7rA1gDXEYWOK+D3Ra9X/+z++KP1IISn4JAjLMw8WqpS7mIldu0y/Zaa4PPZBWm9331Xvi09+aSyqZ8aAUfttblWhW2YI4qot9829621RE2MA8D8HJWE3IoK81lRR2c3z5xhighHUBu5WS3SVF1ySP0KpTOMzpoL84CHUmz1JXUl5BJi6rOrq62tGRYscMyqgrcBR4VjQL0yqy6QKtXqC1t9hNT9RO6Zl+YUVktt4oqowZVt01ZGCrUo1ceRaMpNLcWQ08MASqkBwC0A4sBSA30BgIss+wEsBnA92GzvfgB3O3vMhoBSutXOJt9Jfjto7X91I/W3awjfgPqYyXWUujJXdpTx401RGetzZsyVKJkSAyyCrVKU0RtuMP12NEgFR8k0KTub+VlJj8HJyzOZslqyaJF1nl0poaG2I2G7Cqmw5Coc8VWilJm48pez3Mu7voVcXochQ9QL+A88wAKacD78kAU3UcLZIHnl5baFaGmfU1Ul/8zbioStNADq2dOUv9OV/Zq0rJISFvyMB3Kz5NVXTflG5fjtN6Z4sjdwlKbhURJyp05l/pj8ftkSruVwVbAkZ7DX1mbNMv3eupXNhlVV2c6baw9KHZ8pqsuZXB6FWG7WrqqKRQRWC28bltkd1CDnH30lo6T8BUx+8vHxrrMIq+3YprLS9nPi7V27ctXiKv9zV6UObIy4ZBhAKc0CMATA/wAcBzNJJpJPAoBnAVxLKXV4Ep4Qolf6AAi3t03Np671D9J4abWMgXd1Q6lrIsA6SmM0vyguZp14QyP1+7RMIdBU4JEp1fLDD65TfCgN+pXy5wIsncy8ecrr60OIdRWOmI7xlDlq2LePmbjOnMn+5+RYPy/1LeRy8zJnzDWV/Kw5coGQHEUqRFsObqQCw7ffOn4N5c7dYGDmpk8/zYSFMBcmEJQKuc88w4SBkSOVt5cLXCila1fb+wPm7cxyMM77jS1b2PfHH5t+O4LaWeWGxNKK4/77WaYBZ1Azyyk9xjvv1F4BqQZpICxLHLVomjBBnQm+HLUxwW3K2As2BrCAenIEBDh+PDVCrlyQvfh420oxV87A1zZglZSm5k/rClw2DKCUVlNKv6CU9gYQBKADWICmEEppN0rpQkppbefpiIs+dckIyW8l8aQrIeQEIaScEFJCCDlDCFlMCOldx3VrEiilC6hrGoPWXA57g6365vXXG7oG9cODDzoeeVEJR017mupsuRIhIXVT7pQp7Fs6a8cjiT/wAAsoVN9CbmPEXn9qaRIrDbBWWOh4e5Rr71KlzLhxjpVnD0cFHIOBmfVKg7lZwi1XasO6deya8RQdr78OXH997ctzFUSjR1DriwjvYSMhsIPIPV9yKaMcYeNG+4Ku1F93zhxmOVHXqA3MY4utW9WZ4Avsp8IiRHmbxx5z/Hjr1zu+j6VZtJxAW17O+hOeh7u6Wjn9EGDeznx9zde5QvFeF+PrrVu3ghCC2bNn48CBA7j11lsREhICQghOOhL6uY5wMiC5PJTSEgB2MhI6xFwXluVyCCFBALjR2U5KqZIRRSiAEAAFYKmUOtR8HiaEvEspnVPHVW3UbN1qnmuuIfAMKMSdv07EX1OXoTQzsmEr08gYNkw+ympd0n78GlSX+SBpS/2OFB3JN2kLJVNpR3xgriaU8jhbwv0dLf3p4+JMGvcuXVxXr6uB995jkcCl2BqQyaHXW89GS/NHu3rW7X//c2z7Q4dMeaXrYsDXkLElLHH3KUV1mS+63rMCt/90r3H5gYVP47/nP4JB51zEIKlvOOfECaeKBCAf/Kih8fdv6BpcXdSVj6wSrproUPIHrq5meW3tRVaWWss1FSXt4ZocTydOnMDChQsxduxYTJs2DcnJyejEIxo2IE4LuYSQLQByKaUTVW7/M4BwSqnqUSultNEKuYQQDVi06OYAKsFMti05B+B5AH8BuEgprSaEeIDN/r4LFrTrFUJIPqVUJuyC8ViPAXgMAFq2bOnK02gUjBpVd2X3mfYVLm0fgdzTtvMLDJz5Mdrc8C8ePtAXC1tduY4wGrdqGPRagFr3pO6+JdBXecBQbW5T60ozQzVo3Ktw95/M6enEz5Pw5/3LUfcGGa5FycxMTb5JgeNIZ+CvtFnx2uBIYBxLARew7b/Lkc6K6vXWpvabNqmvQ31SFwFUaiPk+kZkwK95OqZsG4bvBhxC7hnnB4ZBsYmYfk4+KWr/pxei/9MLkXO6I77uloCm1qfWJ1wRIpdfWyCwxJZwnpVlO3+9ZayCulDC1YXLGRdyd+3ahe3bt2PgwIGuP4gTuGImdwRYDlm1DARwJUloCwDUJGPAk5RSK3sgSulymWVVADYRQnYA2AHgGgBvEEK+pZQWWm5fs883AL4BgH79+l2F1vWOMfSVtxDaKQGZx3rh+vdfAACkHeqLX8ZvQFkOk9j8o1JRkhkBqneDb3gmhs15CwAQEJ0Gk2u5OZG9D6Oq1Bd5Z53MeC6DZ2ABHtw1CHs/eh7Hlj7o8vIBJsQ+c6kFvIIKceTbR2DQa7HttbcR0fMYAlqk4MZPnkHKrmH49VbzIOi1iars7lOKFoP34OJ/oxAzZDeyTnRDZWGQ3f00btV4udyUX6fbpJ9x6MvpuLx3sOOVqMErOA8V+XVkP3uF4uZVDl2FF5riQLipaMLrkvoIKNJR0g3qdMCIEbX3PaxP6kLIddQfO7zHMTx2uJfx/xMnmRL2g4ASVJf5Kuxln8he5g6zif+Nwoqb/kWbGzZh8obRAIDQTmdw83cPIj2uHw59+RTq/BknBlmlamNn7Vp5f8wrC/mxjhzXvf88kndci/PrXex7cIUQ8/Gz8Dl7VHadmxsQodDvlAJw/YhSnrIOvZDy3KcuKYsLuQsWLGh0Ai7gQp9cB495RQhohJCPAEyv+TuDUupwWBBKaQUArkP3A4tEfVWQluZYShtH0HpWYMTc19Bt0s9GARcAovrFYUZaBDRu1bj2jdfwTHILDH9tLgCKW75njn7n1rHO+5r/LYRfJA+LSTHkxXcxR0fwyMG+ePJUJ8zRETy4tz+I1rzXanXtNszMbIYpOwfj9hV3o82NGxFhMehQYtyiRxHa+TRu/u4hKD0mRGM9kup+3zLc/N1Uq7pIcfNitp13/3UzvIKYHqX3I9+i77RFuPvv8bhv0yjcsuRBeAUWof24dbj52wcR2CpJVb2VePhAX9z7z42Yo9NgyvZhGL/4YVX73brsPuPvj0JZRJ+pO4fgnrVjalWPNqP/wazsZrh+3iz7GwvgE5aFpy9F48USHzybGgmvoKYX+eRqEnI9AwsQ0v4sRn34HCatH91g9dDrrSMON1YcE3LVDVkcmX158nQ7MwFXSs+pNkJYW6D1qMRdq29B13tWAADa3LgRE35nuaUKk1k2w80vsBw/if/eiA8Ci7GoJwsb0nPKUty08H8Y8tK76iteC0YvmI6Xyz3QeuQWaD0bILpkLfnhB5a/uTZoPSrh7lsCrWcF2oz+B3JtyN23BJ3v/B03fvK00+9aAAhokYJhc96Eb4T6eaebv30Qc3QajF/8kOzYwrL8wbM+xD1/j4dPqIO+DIIrjtLSUpw9exbh4eGYOnWq7DbXXXcdCCFWn3GuDtCgAKFOzokTQgwAMiilUSq21QIoAVBFKQ106sANDCFkHoDZNX9nU0o/cqIsX7Drorqsfv360UPORMhoBHh61i5lUGDLS5j8zw1YfuNmlGRGwMO3FBUFwZBqIzvf+Tvu/PUu5JzuiNBOZ5Cw6g4cWPgspmwbbrPs6nIvLBl4ENOOdTcuSz3QH9H9Dyjuk32yCw5+OR2Hv5kGUA1eLPWEm6f1ib3tZv2sEY0eYV1Pwi8yA8PfeA0tBuw3ritOj8SClqk12m+Ku/8ej9COpxHcNhE735mDEyvuxbhFj6LlUFM+l/MbxmDPvBdxee8g9HpwCU7+djcqCwMx8u1XMOTF94zbnfn7Fux+/2VM2T4UoARaD3lnkaLL0VjY2naIwDY3bkSn2//A+icWQaoNDu9+HI8dsc6l8u01ccg40sdqudazAjct+B+KLsfg2jdYlKt3vapg0LnjxvnPoP/TLHHjyV/vxup7f7FZJ417FR6P74KQdhewbOR2jP3qMYR2Yq7yZ/66FeufWITSrAibZdQWrWcF9JWeGDBjPtqM+hcr71oJrUdVg80ie/gXoc0Nm3B61Z1Qq63vdMdKTPjN5IGyf8Gz+Pc51yU3Du1yEmM+ewpeQfnIiu+BNY8ssfIRJBo92t70D7pM/BXuvqX44y7HRpu9eyvn8VWDu28Jbpz/LJp1PIN9859Dy+E7sOXFD2DQa3HTwunIOd0ZPs1ykXe+HU6suM9mWWxG3Btu3mXQlfvUvlIKPHKwNyJ7HzX+n988E2XZ4S4/jj0SE4Fly4A33qj3QzvMF18ATz1lf7sb5z+DsK4nsWnmpyi4GGtzhnXdOnXBtcK6xWPaURYp7eAXT2H/pzNRlBIDg84dL5V5QOtRLTub6x99GeW5zaCr8IbWswKzc4Pg5mWaro//6T50v+8nAMC219/Eng9ehLtPGSqLrIdbPad8j463r0KH8cxi59BXTyC47QUQjQG/jF9v9jyGdjmJ3g8vxn/PfwSqd4MjM3/tx68xup1w3vMth77SwUTATQitRyVeKjM/v6RtI9B6xDYsH/0vLm5mvlk3ffYk+j1hMpHi7zuOh38Rbv72IeybPwup+61nybQelRj6ytsI7ZSArXPexfhvHkHLYTuRdqgvji55BIe/mYYud/2KrPgeyEnoYrXv84V+0LqbtD2/3vYXzq013SvvZjk1765g6Cq8MSsnyKggB4ALG0ej7eiNWDH2HyRuUqdc8w7JxU2fPYXETaPRY8oPuLRtJA4sfAYVBcHQelRCX+Vpv5A6Jrz7ceSc7mTlssXZsCEBoaG2Xd/kCAqyjseh9aiERqtHdR28F0xQ+ITmoKIgSLUvPtHqQPVaWD7nsbEsAwQA7NmzB0OGDMGDDz6IJQqh/4ODgzF79mw89NBDZst9fHwQoCIUdkJCAjp3tn+tCSFxlFKrkKEOC7mEkACw6MmcJADZYLlhlXo9UrPPgwCeBnCQUjpA5fFmAfiCUlpud2OVEEKuARBKKd1Qy/0/BMCnhJ6nlH5oa3sV5V2VQm5tfeYe2DocLYeZh3FM2TsIMYP24vSft2HPBy9h6q5B0GgN+DgiG+W5ocbtPAML8PTFlvAMYLkk9n48C4OeM13un27YjKStIzFx1W3oeIt1QtJv+x9CxuG+AJgv1YxU5bwnO958vWaWmHFx8/VYMfafmkECY9L60Wh7o7nT2sq7fjcKGBv+9zn8Ikxm1K5COgj28C/CkwkdkH2yG9Y/+TVACUK7nDIOTN7xrDbVmRgAagpWzpUJAJCyewiWXssFboo5OjaV9t/zH6LlsB1I2nodbpw/g5XpoQM1mKZ8fEKzMTPDfFC+8+1Xsf2NN43/pULzlpffw555L8qem5tXOZ5M6ICAGHPh/Mh3D6P3wyx/SkFSK3zeLknt5VKFV3Ae7vx1ImKv24Lqci+4e5vPWCzqGY/sk91cekw1TDveBWFdEvD90D1I3TfI7vaegQUY/tpcDHjmUwCmwdm7XlW46bOn0Pvhb/G+f5nNQerQl9+GX2Q6sk91xdHvHzLb1jOgELPzgqz2ec+nwmyAc830hRj96TPG/5+1SUJhcisVZ+wKKKbuHmSmdAIAfbUbqF5rJlgAwL5PZ2Dvh89bBavzCsrHrBxz5cbaxxbj6BL7zq5Eq0O3e35Gmxs3YsNTX6GqRD76zbVvvGbVP/w1dSnif3pAsezoAftQkNS6iQXXY2OVFoP2orIooM6fJTevcvR66DvctNA8zMa6ad/gyHePyu7z9dfA44/bLtcnNBv3bx6BsK6n8GlMKkrSzecH+jy6CGO/ehynfp+IVZN+g3dILu78bQJ8w7IQ1vUUAGD9U19i7BdPKh7j39kfYf8nz6k4S2DQ7A9w/XvmfWlRahQWtkpFYMtLGPn2y+g2eYVx3ep7V+CGj2dg26vvwK95Oka8yZJ6b5j+BeK+NtXp3k3XI/Y6U/6kc+vHov1YFsb28v4ByD/fDuHd4rHuiW+QdkDVcNBhgtudQ3T//bi0bSSK06KdLu/25feg692/YvNL7+Pc2putBEcGxYz0cPiGKTtoSq+FlN3zXsC+j2fhrlW3wcOvGBE9jwMAMo93x+I+x8Dfu8Ftz6MsOwzXfzAbfR5drHicU79PRJeJvwMA4r5+HH0f/xp559ph3eOLcf9mU/qGtIP9EHUNG0+e/PVuXPhnDCqL/TFx5Z32LokZf01ZhvjlysmoA1slYdrxrvDwtR0q+J+nP8Phb6Zh4HMfofXILTiy+DFc3jsYlcX+qCpWlyPIPyoVpVnhRqFO61GJRw71RliXBCy/aRMu/jcK5iILRevrtmDcV48huG0idr//EhL/vRH3bx6JLa+8iz0fvGjcvrZCrhzN+xwCCFCY3FJWMUm0OngH56MstxlACfyj0lBRECSrbNO4VyG47QWUpEbDNzID5XkhKM8NZZY+7c6jPDcEBUn2EzUTjR6RvY/AUO2GzPgeZm4GUiH3888/x//+9z98/fXXmDZtmlU5Fy5cQLt27bBr1y4MqWUY9IYQcl8H8Jp0ERw3P36GUvq5yuMZAGQCmAfge0ppgYPHkpY1FMCLAMYAmEspfdPOLnJlfASAvz2cFnBryrweAE8xP4FSanfK4uoUcil6PLAUtyxR56u6uN9hZB6Vz87k3SwHVSV+xsG3f1QqitObSx5mCq1nJcZ9/RjK80KQfbIbEv6YYOVP2v2+Zbj1hylW5S9oeRnFadFoN24trn/veYR1YeF6D3/zGNY/yaIOhXU9YTZjDDAB9/SqCYqCwC+3rkFkryMYMZc9gtmnOmPf/FnIv9AOaYf64cVi5VmGb/ocRbfJy3Hx3xuNmmRbcAH22A9TsXXOu6gu88Ejh3qjuswH315zGH6RGXj6ornQwZUA138wG4Oe+wjJO4dh2UjuoEfxUpkntB7V2PHWa9gxlykAlM51XlCR1cDeKzgPs7KbGf+vvncFTv46CQDQ5a5fUF3mg55Tv0en2/7EhU03IiehMwY8swDn/7kJv9y8Dtc89TmGvToXPs3ykLJ3EJYOMyWkbN73EG5deh9+GrUVJRkm5UVo51Pofu9PqCrxw+73ZaLzgPkez8oNMtOKy/FhSIHsrEpt8A3PRN8nvkTc108oCisDZnyMGz5k+rjd77+E3R+8aBwkaD0q8fSlFkg/dA0OfP40Jq8bA12Fp1GAy7/QBksG70ebUf/i9uWTcXHLdcZB67GlU7DmkSXmPnbEgJghu1FZFGBmhnl27Xh0GL8W/z3/IfbNn4W2N23ApLVjreqaf6ENvuh4HgBBs04JeOKE+QBy3eOLcOTbWuSHqCEgJhlaz0rkn28vu77tmPW4/r3ncWDBDAx/7Q0rJYmUi5uvR+z1m62WZ8V3w7KRO1BREIzmfQ/h4f3XWG2TcaQXlo3cATevCmNsACka9yrcu/EGtBpu7tj65wM/Ws0Yaz0r8FKpNwA2E5e4aTRGzZuFgqRYrBizEdJBXJsbN6Lfk5/DM6DIWLYjMzBSPPyKYdBrazUr7eFXDI2bDp6BhSi81Bpys4LWvuAU177xupkwf/zH+7H7g5cw7JW30GXib/h39sc4+NkzcAVEo8dzWc2Ms1a5ZzqgWUdThC1La5TAlpcQ0DIZKbuG2S377r/Hof3Y9fj7wR9w/EfrdwfR6vBKJRuYLxm0Hw/ttS0ALmiVgrLsMNzw8Qz0e+IrpB/uje/6H1Z1npzeDy/GuEWPIftkF6Mg/e01cXjkYF+HyvnzgR9x4ufJ6Dvta4z53DRNvv6Jr3F48TT4hGVhZrq1BU3CH3fi6A8PIf3gNYi9/j9kHOuFh/b2x2+3rcGl7SMAsP6KaPUOtbk5OlO7etvNAFuzz/7Rl9H9vh8RGJOMvo9/jTUPL8GxZVPMxgRcccs5uuQhrH3MlHj61h/uN86kA8CKcRsw6LkPcW7deNz48Uz8NWUZbl1qrnz674V5OLzocTxfYF94WzlxJVL3D8QzyS3Mlu9+/yUMefE9FKdHYvtrb2P8YhUR42qYH5mFspwwDHzuQ4z64Hm72y/uewRZ8d1x83cPIiAmBSdW3Ivx35iUPsrXmWLSujFoO3qjcUl1mTfcfeTnr0qzQ2UVBRlHemHJ4P2Ks6yAudI87VBf5J3rgIgex4xtW8rZteNx7PuHMPGPOxTL4/D3d+2EXPN+zsO3BM06nja7VBlHasarhBonFgJaXoJvWDYqi/1hqPKAd7Nc4/aFKTEoq7FI03pUIry7dfZSfZUHtB4my8KsE92hr1SeLfePvgy/SHNzd0O1GzTuOuSc7owW4b5GIfehhx7C999/j0OHDqFvX+u+4rfffsPkyZNRVFQEH5/azVQ3lJArzZip3m4FSAXwNaX0HQeO9zaAGQC8AFQBWAfgZ7BUPTYTThBC3AH0AnALgMkAWtfU9QCAxyilx9XWo6Y8qYA7y1YkZMk+hNq4yIQQTwDbAQwA8z1voUaQv9qE3KhrDuD25fcguA1Loreo13FE9TuI48umwMO/GLNzzcNaVhQG4KNmsvG7XIrWswIDZ8zHkW8fxfDXXwfVa7HrvVeshI4+074yat4vbr6evfRqZjU/b5eIosstrMxIrp83C4NmsiZ2ef8AbHlxHpJ32ja39grKR4vBexDdfz+2z30DABDZ62hN5+mYVsHNu8ym0MxNwXe99zLivn4Cz1yKsdrG0uTOu1kOnstkA/tVk3/Gqd/uMZux+6prAkLan4VXUIHiTNSoD5/DwBnzzZadWjkBXSaY8hulx/XB0mt3QVfhbfO8frllLS5uvh6Dn/8AQ15812hmvvaxxbi041oMmvkR+jz2jXHfz9okofPE33Dk20dRWRhoNfAG2CzKDR/OQtrBfjjz923ocPPfOLduPEbMfQ2bX3ofez98AY7iF5mOJxI64N+ZnyLhjwmIGbYT9/x1s3H9kkH7kX64D6jeDTFDd2LMZ0/i6A8P4caPZwIAsk50RXg3lrPufb8yeDfLlb1fUs6tG4dfb10LotHjlSqT9YHUfH/piB3ocf9SBLdJRH5iG+NMuRLpcX3QvC8bgH8QUAJdpSdACaYd74rQTmdwbv1YuHlWGgXIzS9+gL0fzcb/ElshsGUK1j/xNY4tneqwSZu7TyleKPIDAGya+QnSD/c1E0i0HpV4vsDfymx/2cjtKE6NxoBn5yPvQjs2ULWYJe39yDcY97W5FvvvB3/A6IXT4enPDHR2vjMHw15522wbfbUb5gWWWJ3Lde+9gMGzLUIU15CyewiWjthhHHhz6wbpYPumhU+h35MsF9vFLddh+Y3/oddD35kNRKV8FJZrZUrv5lUO/6g0FKc3txAqKK7530KM/uRZxX1t4dc8Dc+mmGbUzv9zE3zDMxHe7QSqSvyw8q4/ENL2PMYtegwnfp4EojGg84TfodEabJRqYuXElTi92rHZJ47UVFKadufr7ieRk9AFGrdqdJu8HLcseRBn/roVv9/5Z82eJuHni05nFZUonBdLvKCr9LT5fgpseQnTz8eCaExDh/+e/xDxy+/DLUumou3ojUjZMxibX5yHy3tMMyTBbc+jIKm1mbWQowS2SsL/LsQa/+cnxmLn26/h+I8P4LZl98E/OhUFiW3Qc+oPAIDlo/+FZ0Ch0Q+4PC8Y3iHMh//EisnY/cFLZrPuPe5fih5TfoB/8zQzxYFalo3cDt+IDBRdjkHagf5mFkFSLAXqhD/uxOlVd+L06jsQEJOCwbM/QLd7f8LvE1Yj+0Q32f5w65x3cHzZFBSnN0er4TvMZj+l/Dx+PVL3DTSz2njbXS8baKv9+L9x95+34sTPk7Bu2mLj+7HN6H8weZ15zIm/pizD6T9vx9SdgxHRw1qAAYAvOpxHfmJbeAYWwFDtjuoyHwTFXsQTpzpC667DL7esxW3L7oVXUCH2L3gWnW5bhcBWyUg/3BtbX30XiRtvqimJotukFbj5uweN/eDKiSuhq/RCy2E7cPzHB5BzqqtsHQbN/gCDZs2DTzMWP+PM37fg9zv+Mtvmjl8mosuEldj03HwcWDDDuLz3w4tB3HQozWBjJn21B+75e7xxfdLWkWg9cqvsNb/wj/n18gnLws3fPoSWw3YYrfUskc5ay6GvdsPaR78zTl7s+2SmcbxxYOHTOL9hHL6d6w9/d/tWUYBpVhQAsk50M06shHWNNyqUC5JaI6h1EopTo+EfzbJ6VBQEwSuowG75eefbgRCK4LYXbG5HDRoQDetHM4/1tBpvegXnIbhNos0yqst84F/RxSjk9urVCwkJCSguLoaHh7XS4YUXXsC8efPga5H0d8yYMfj999/tnRqAhhFyA2EyVyYAEsHMlfvb2M0AoEgparCKY0aDpdqZDEAL08xxCoBjNcfPA0vhEwyWi7YNgJ4A+JUnAC4AeJVSatuhT74OH4ClAQKAmZRSVQ5qhJBrAbwK4AcA2yill2uWuwMYDuA9sMjKAPACpVR+dGPB1STkevgX4fl80+zXijEbkfjvjYrbN+97CKVZ4ShKaVxBvJXSOsj56jIowrqeRPbJrmio6LbSaJxyGHRafByei8qiQNz242R0m/SzcR2flbak55TvawJrmUy68863xVddzigOWCwZv/gh9HpQPjhL5rEeWDZyh80ZU61HJZ442QlBsUmqjidHxtGeiOxlHkyd38vgtudRkR+M8jw+60xx/5YRaDV8BxL/G4V9H8+uacNMRxg9YB+0npVI3nGtzJEohr8218z0XY6za8fDoHNDp9v+NFv+90PfI7RzgqLgxNn3yUwk7xiO9Lh+6PvEl9j1zhyjkmD4a29g+Gtz8VXXBJRmhZvNpsux7fU3seudVxHQIgXjvnkEPqHZaN7H3EHWvN1TPHKot+L1vPm7qeg5xZR0eNd7L2Pbq6p1pWZKJg5XOET1348xnz+B5n2O4MTPk9Bt0s/QV7njg4BSGf8lZZ1u6+s2475N5hYSJRkRWH3fz7i0jQ2OLWeo1z/1JY4sfkzS7ileKPaBu3cFVozbgMKk1sg90wljv5xmVLbE/3QfCpNboSg1GpG9jqDPo4uxqNdxZJ9gViEth2/Hnb9MhG+4dWAYaiBIj+uLpSN24pqnPseoebNx4udJyD3bEVXF/jiw8Bnc8v0UdL/XlBCAK6NADLh/80izGebMYz2wuK9VQgFFhrz0Dka+Vft08Ouf+BpJ20aix/3LMPRl0/3//c5VxpmYH4btdigKu2dgAfo98SVGvv0KANOgNuNoT3x7zWErQYUrEVaM/QfJO4bjxRKTEsCemXD0wL14cNdg/PP0Zzj05XTF7QDg1qX3ofu9y1FV6oOPmhU4ndf2ttuAP/9Ut+3EP25Dx1uZkPJxeI6kH1Omx/1Lccv3U43/1bpHtBqxFff/d526ilnALU9+ve1vFCa3RGVBEAbN/gClmZHwDCzEwBnzsWLMRgx87kO0GcWM5YpSWqA0K9yobJMj/XBvq/6K81XXBBSnRsPNu9xMiOazkpuem49DX0y3eb+8m+WYuVFJ0bhVw6Bzg2U/Y6lMe9+vrEZJKB9hL6TDGVC9FvkX2lmts+f7Gtb1BNy8y5F+yNoSRQnLvu2b3seQFd8DRKPH2K+mGRWglq5Kcox4cw7yE9vg2A/mfpx+kemYsmMIgttcRHleMBa0ugxduQ+C251D+3FrjUpdKSm7mRIoZshuozJQ41YNn7BsM4Xbmke+w7l1481Nhmsigrv7lOLhA32NcT3GntmAFu0jkHO6C/RV7gjrfArQUOSd7WC6rsQA3/AsBLQwtwgqSomBX2Q6NO46lGZGoDi9OajeDaGdEuDuWyp7PQqTW8I/KhUaNz0Kk1vCM6AIAGSF4PS4vvAOyQMo4OZdAY17FQx6NxRfboGAlsnwDctGdakPck6ze0U0engGFpoJuCUZkagq9kd1uTe8Q/JACIVPWDY0Wh2qUnqiTawWlZWV8Pf3R48ePaAki4waNQpBQUF43yIHV2BgIMJU5qSsdyFXpuBtAHIopdajWRdDCIkCyxP7EACpvYbcSfAeQgc2+7sIwEZbs6o2jtsSwKWavwYwodoWH3G/WkLICABSNVQ52IxtIADeCxoAvE8pfUVtnZq6kNu7t7rIym3HrMekNSySR+Gllvis7SU7ezRufMMzMSONaSz3fToDhz7/HwqSYu3s1bAMmv0BAqIvI/N4TxBCcfzHB4xBNT6JyjAL4EQ0etzw0Uwkbb0OZ9fcKlse0ejR++FvMfYrk/PaH/f8ioSVdzlcN5/QbDx8oC+oXotNz32Cs3/fpnpfy+BKALCwdTJCOyeYCfYrxmxEzulOGDH3VfR4YJlVOdteewv+0Zfx73OfyM4cc/o89jXGfvmE2TJqIFjU8wQej2cacubv9wgAgjt/nYDOd1p7LpRmheHy3sE48fNklGQ0txlMjc/kBcVexMP7+xlnWACgPD8IX3U5U8sARRTPpkbCL8LcmMaWSfHg59/Hde++BIBd56LL5jMnUh/33DMdsH3uXCZcgSkN+j7xJQY+a9ItlmaFYfsbb8I3IhMhbc+jWcfTSD0wABX5IRg25y18N/AA2tywyShUFaW0sDJBXtzvMB49ZDI7VZp9UYtnYAGmn4+Fd3ABAGZKWpxqblrY/d4fUZoVYWxjVaU+KMsKR+7ZjkZTvv9emId9H882K3fE3FdxzXR5Lx9Lf2YAiOwTh0cOmN757/uXWpl6Tlh5u5VSRIm985/DoJkfo6IwAJe2jzDGLLAVlMw3PBNDX34b5fkhGP4q8w7KO98Wu955FSm7h+KRuF4oSGxj9D3kHPn2EfR+5FucXXMzPAOKUHAxFvsXPous4+aB7AJapKAkIxIGnTs63rraKOgm/nsD/CLT8cPw3bh30yhEX3MQAGszX3Y6hz6PLcLRJQ8jtHMCpmyXNzH+5dY1OL9uvNXydmPXGWeafr3tL9z9J+vnKgoDZFOvcaSz2JaxIuSh6PXgEpxdc4usWbuj3HIL8Ld1mAl5iAHNOpxFVYmfVfu1DUXXu39BRVEgLmywdkmwty8IRUTPY5iyfSgW9z2K6959EZ1uW43lo/9Dyu4h+N/FllZ9jj3e9apCRM9jeGhvf7OZ8ZTdQxDc9jz8IjONy9Y8vMSYuk9OMQaYK+f8o1LNTIeZe8521JVSWuNehS4TfseplRNtmus2FF3vWQGNVm80ybZUBBsVZk7S5saNmLz+Jvx2x2oMefE9q8CgVaU++LbfEeSd66CiNHXGqFLf9bFnNqBlZ/b8Vpd7w93bZHKdfaorqIEgtONpaCTuS5QSEGIufmTF94C+it1H72Y5CIxJhkHvBkO1G9x9y0ANBBlHe5uZzEvrGhSbyARaAFSvQcZR64Celuca3PYCvIIKUJzeHBqNHr6S56mqxA8lGZGyaR49/IrRrOMZVOY0R5uwaMTFxaFfv3547LHHsGjRItmjNWvWDHPmzMGMGTNk16uhwYXchoIQ0g1sJnQAgCgAYWAmzblgQugpsPyzuyml8nYL6o/VGsBFB3aZSyl9o2bfZgCmAhgEoDuAUAABAMpqytwJ4BtKqbwtigJNXchVM4vLzXo4HwQWo7rUrw5rVT9o3KsQ1vWkor9wU8AzoBA+oTnIT7SemVZLl4m/4o6f70HO6Y74uttpF9ZOPd4hufCLSkP2iW6Qvjw8AwsQEJNinB0zh+KZ5Gj4R6WbBSKzB5/FsUdZbgh2vPkGblrwtNny1fctR8LKidYzBMSAiJ7HjMLadwMOIj3Oqq+Hxq0abt7l6DZpBZK2jUR5TqiqGRoluPnmyV/vUR0pVetRCa1HlWIQpR73L4WbVwUOL7YOYsFpMWgPpu60H8TC0qdLKjj2ffxLM5/BymI//PXAT4qKGUdw9ymFxr3abj5o7gcph1Ke1NhR/+Lef8ytWJgiw9o/mEM0esWZk9DOp/B4fFek7B6CmCG7jct/v3MVzvx1G9qO2WBUMnL4TExAixQ8ebo9qEGD74fuxR3L78EPw/bAK6gAQbEXcXnvILNZTs7hxY9i/RPcBcB80KZxqwalpNbmtnJBuNRQnh+EpdfugndIHqZsG45jP0zFmkfkLUWIVofH47ugWYdzxmULYy/h+g9mo+tdv6EkIwJfdT1tdf/5LPbWV9/G7vdU67NdxmOPAd98Y3+7xoy7bwkM1e7QelQhoGUyhr74LjKO9VL0JbWMFgwA4795GF3u+hULWqaistgfnW5fhdOr74CbV4WVEsiveRqi++9H0eUY9JzyPba8/L5s38V9ces3OF7jRBonQEreuXZYOmKnS4LdufuW4IVC8/uwfPS/SNo6EtSgZbOXLlYCEI0ePmHZCO2cgK9fiEB02yiEtDtvd7+iyy1Qlh0GjZsOPqE58GuejuK0KFQV+yu+BxlM6WNb6cq2cfcuh77aXdU5a9yrENHD2lOzOC3KKgieJV7BedDqAtG5o32Lu4sXL6JNmzbYunUrRowYYXd7Ja5aIfdq50oXcqVas3WPL0LaoWuatFAokMfDr9hORw/cfTfw66/1VKE6xM2rHJPWjcG2198y+oPe/tMkdL3nF5RmhWHn269ZRXLNONILv9yyDroKLzu+j2wG5cLGm2rSaV3BEAPmVJu/ZAsvtUTR5RgzQY2TevAaBLW+iMV9jpm9xKdsH4qYIbtx8re7sHpywzSwsK4nMPLtl9Hh5jUAmO/ugQXP2lU++IZnQl/tDjfPSpTnhbgk7UZwmwtw9yuBh1+JmZ+nb3gm7lkzFs37HrYSqKP678dDe0ypTSqL/BV94ZZdtw1eQflI/PfGOkmjxOl+749mwX3yzrXDxhkLkXe2A546a226uX/BMzUz0eylFNgqCSUZkXYVN1N2DkbMoL1I2jYCP43aipbDt+OBLSPYMc+3xZqHfwDRGHDN/xZg7SNLMGXHEBh0bvi231FXnapDlJcD3sqGJk2aNqP/gaHaHddMX4j1TyzCwOc+gq7Mxyw6v6D+4EJo0eVoHP9xCna/97LN1Fu1YexXjxkjS3/T56iVlUddwgNPuXmVI6zrSTbjeqSPcZYUAKpLfZFz2jURmF1NQEwyfMOzUFkUAINOy6wJVVow+foCKmROrFy5EhMnTkRCQgKCgoLM1oWGhsLNTZ0is1EJuYQQD7BATy0A+MKGDQCl1NruT6CapizkGgyA1oYiqHm/g3h4H3Px3v3+S9g6p3aJ6ocOBXbtUl7fowdw3IHQY8OHAzt22N9O4FqOHQN61t/7q0Fx8y7D5A2j0XLoLmz43+eI+0pFIs+rkJihO6HR6nFpx3BE9j6CjMN9IH3dtBi8G1N3DAUAfBKdLjt70LzfQUz4dQKWj9mEvLMd66vqsviEZiMgJkU2f3RjgGh1GPnOyzj582RkHutltm7oK28Zo73L8Z5vOQihNk35XY2HXzH0VR6IveFfnF8/1mwAp3GvAiG0TvJxBrRIwZTtQxHYKll2/ZpHvrPyM6wvKAWaNwcyMuxvK7hy6dIFOGUdZNjlMAW2H+oyngizUtHU6THkMIuuLE2rSAxw96qAwaCFocod1AnXl7qFQuOmq5Wfv58f0KmT/e1eeuklK19cACCEIC8vz0rwVaJRCLk1EYLfAfOXVaOuoZTS2of/EzRpIVenA9xtPFtDXnwXI99+xRjavrYUFQG2ck3v3g04krprxAhg27ZaV6dR8dBDgELu7jojNha46IjRfw0XL7J9BQL1UDwa1wtBsRfxYUhRQ1fmiscnLAsevqWYfr5NTXqwN9D/6QUovNQKZ/66vaGrV68opUQDzP2iH32UzYp8+mnd1+nmm5k/7sCBwP799rcXXLncey+w3BRXDocOAf2svVuuWlJTgWg7KZVdmSe3qaFWyHUVzgq5TguahBA3ABsBDANTp2QBCAcLppQG5oPK7X5KwHxmBVcxBoVsEB7+RRjy0rsY8vwHyDnd0elgG/7+trWWSvVQwkud2yEefRRYrJyf3Yy+fYG4OMfq4QoWLKh/IffJJ4HZs+1vZ4lKqxaBQAJxKOqvwDnKssNRlm0elOfAwmcbrkINSGVRIL69Jg6eAUUoSolBfmJb9Jy6hM1mS8y0uX/svn3sU5dcVxO8+O+/ge3bgRdfBBITge7dgXiHooEImjqUAjExQEpKQ9ekcRJl2y31qseR1J+NAVfMpT8MFgAqDUA/Sim3C8uilLYE4AdgJIA9YEL1HEqpmJe5itHr5Zff/tNkDHn+AwDAiRX3ueRYbdoor5MzYvjrL+tlnBck6U1TU83/S2nhQDDKa68F7nPNqQJggqQaCDENfOqL2naOjlxPjiuvqSsYpC6lnkAguALIONIHl7aPMAbmO/bDQzj6/cOy2z4mH3us1rz0EnBeIR5OeDgwcSLQp8YqPsZ2qmxBA/O8fDwtM4YOdaxMSgGNZOTvqLJfIGhKuELInQQWIvEVSqlV4jFKqYFSuh3AtQB2AVhCCGmcjkeCekGuUx310Uy0H7cOAHD0+wex672XXXKs++9n36tXW6+zFHK7d2dpFpSQmj5HRbFZWDnGOpA5QaMBlkm80x0RPCkFRlrkptfp5Le1hBBrofOJJ+S3tUePHuqPWV885ITb2/z5rqsHZ8UK15fZmJlrO6WvQCBoQH74gcWXqI9ZqzdF7Ce79HFiRPzxx45tT6npXXz6tLogQgJBU8UVQm63mu+VFsvNQgtRSvUAZoLlhp3lguMKmiiWgS+8gvKN+S+XDN6HtY8ukY30NmeO48e66y7WqbeVyXRjMDD/FI6PQsBPHvTIUkhTEnId8W+ZPdu8XEdd5C3rJCfkUio/ey7dd8gQ4Msvzdf/8ou6OijVOTDQ/L/GTm+jdD0d5a+/mPBfW183f9vBnmuFvXO/0rC89wKBgDF5MtCyZd0eIyYG6NXL9N+yj/b1BYYNqx/F43PPsW9X9e9XIrbe+5HOZ9yxgr+PNBrmY1mbsVVtUTrXDz+0vV8HNSlvmxj+/oCn6+Pf1SlXo7myP4AiSmmZZFklmJmyGZTSEwCKwfx3BVcp7SyyONy+giUHX9z3CNIODFDc77bb2HdttJ5+Mul1KQV++sn0P1gm80pSEjPxAqw7ozZtgFGjHK8L55NPTGXXFkvhif9v3dp8uWXHZPlfJgge7r5bXR3Uzlw/+qjyOkpZAAxncHdn5fDZeLmIzN26WS+zxEOSas5Vs7pXgpD7luPpR5ssFy82nXQr9e120JCoDMjZqFm+HLh0yfTfmdifSlkKPDyAI0dqX64cHWsRgPzkSaY8Xr8e+Ocf19bnSub669k3IfaVA3JjG1tQCtxYk2qb93G1NVmWWqE5i1yQ0IQE0281FgEbNlgv++232rXd+qB1a/tBrgTO4YqhVxYAyzmwPABehBCzITwhhADwAOBcRCHBFYPGrRptb9wEg15jlZoCMPd7dWYwEBsL/PEH8D9JGlLfmjjg3FyHm1jygEz33gu0asUGJV9/zYJYWfoweTiRb1xO6KvtTC7XxHLz5ebN2TcfFFoKtd7epmX//GPt15OXp+7477zDPpypU63rxvF1bZo8K2yZbZ09C+zdy9IR2UMa+XvGjNrX59ZbTb9tpcxqKjgyExMa6vzxZs8GSkpsR0ivK6KibJtcv+wabwpZ5s1zTCBYsEDdIO6ZZ2pfp8ZCY3yOVqywnS3AHrV9r8XEACdOMNNjJezdc3uzMjfcYPrNBS97HDrE/H4BoH179j1mjGv6hIamrvK1S9vAzTebX2t7batLFybIKWE5ZqEU+Owz4OhRU7wLpTgpHJ0OePBB6+VVVXUbaTdEkhpejSB+003WyyZMMHdD++gj5+vlKjw9zc9R4HpcIeReBuBGCJEaVZyo+bZsciMAeAIodMFxzSCEOKjPEjQG2o9fAwBYOWGV7HqDgfl8KpkSO8IddwALFwLFxWxg0p+l4jUOnLhWk8/Y8o4/LAyYNo39Tk42fyHNnWu7bpSaXiqWyA0w5AY83NR67lzgNYt0lLyMIUPYvtwMjhBgzRr2IrNk4EDzfeWOyWe15V6wUjPgqCjzgadSkCjpYKmusDVga9mSnbe9QV1AgGnW9a67nKuPtF24eiZ3wgT12770knPHGjOGfTui7Z882bljAsB33zHFSEOYR7m5mY4rZ3r99tt1d2wPD2D0aDb7lZ5uf3t7gt9jjwEFBey5lRsEyuFIajUlPvjA+TIsqax0fZn2+OIL2+s9PJyzQuHvISUuXpSfoTp4kAkYU6Yo72urj5euV+Kaa0y/1fZhffsCP/7IZqsdEf6bghAsZ+1lC6VZuvnzzfsQSs0VzdL7JqdIt3yf2nofvPii+X9KWf8mtXSyJ+RqtfIWdJWVwM6dwIEDdZMFITzcNAlhTxmkNOEgF3/EHpTKx3EB2OxyfZgY9+3Lxn7du9f9sa5kXDH02l7zLX0trgZLJ/QRIWQiIaQ9IWQCgKVgQaq2uOC4luwhhIhYgU0JYsDElXcCABL/k5eCDAY2+1Za6rrD+vkBkyYpr+eDRnsdP8D8b9XUbfp062XSjrdVK+V9+cujf38m6M6bZ/Kf5QMP/gLg3xoNMH68fLl8m9oKD126KK/jL1TL2bdNm5T3cYUCA7AehEnPj/9WOmduci2957WZYfnhB5Z/GQBur0kPOm6c+cz4LAciErz6qvxyR16yjz8uv3zlSnX+x5ZtTA2EAD//LL+uc2eWw9pSYWOJWmuCukCjMbWVhx+2DspWH4L3mDHMH8/evVYzu+moj7TUn7O2PPWU82VYUlLi+jLtIb33coEJNRp17wolundnM2KUAmvXAhcuAFlZpvWtW8u7XkREyJcn1+8podS27mSvZbM+VU0747OGnp6O+x031Si/XGlsyalT7F6WlVmv8/Mzt+Si1PRekApllMpH3x5Q49XF26Yj/ZHc+1tNwMqnnmJKc2mgyWbNmHLimmvUWWrt2qV8XKU+hwt50mdMbvbYVipGtc/ExYum39w9zhKlmWtXz2gTwpQqjdF6pSnhCiH3bzCBVhLCB9+BzeaGAvgFwGkAvwJoAaAUQF3E3mwNYD8hxI5elEEIqWPjSYE9+jzKEgUeXvyoWf5AKdIXnz3hbMSI2tWDm/pyP9aRI9lLeubM2pUHsIHx00/b3kY6gDh3jmlF5QQJy05u9mzrlxu/TnxbR0w8HRXmpOlw+L5r1rDZDF9flndRmsIiJ8f541dU2I8EainkSmejpG1m7VrrfYfVRAkwGGovwEyfzmZVBg9mQtzEiWxQ/uef5jPcaqNRA8qzJ9wc3Zky7rxTXvmitL/BYIpWroZ77pEXCghhwrVaU//6nMk9dw44c8Z8GaXmM9Pbtqkry1WBUuw9H1qtbXNl6QyLKy0KRo+2vd4Z95LGhLT9yQ3ECXFeQOMznuPGsVgPYRYOXbUVovl+Svf9zTfl33M8doF0P+lvbt1hyXvv2a7P778Dr78uv64pBLThdRwgCR2ydCl7P0kVEf7+TJnn6Wnu28/fYYSwWcrCGptGy3SEsTVJNtu0YabLlJq/0x25VtLncOdOecUpt/CQu39ceCSEnaP02NI+Tk3fwmfr77mHxfyQ7tO/P1Pu3HGH+T6bN7PvU6eYFV5ZmXV8CEJsx9tQe70sY5koIac0tPSNDgpSH73cMk6NK/Hxqb1JtJzg3tTiIrjilbcXQCwAo7cjpbQawPUAfgYLQsWb2C4AIyilp11wXEuGAKgGsJUQMlFpI0JIKCHkTQCXlLYR1A+tRmwFAPzz9OdW67gg5cjgYcKE2g2s7r6b7cdntcLD2UtLSUNri4UL2fe33zJfOVtIO3h3dzbol/oMc/hARe7ceCfKO9hBg5gZ1PffKx/XUlkgLbdDB/MXrtzLwc3N2kRu/HiTz2a3bmyQxpUOXl7KdbE8vpTVq02BwTw97Zu+WdZVOriXrrMVaVpqMufMIJ23JV9fVg9/f5NQoNasi0cGt+TQIdNgu2dP84GG3Eva1ktp1izrQYUlUkWKUqARaSorqRDvCqHKVYNfJQF98GDT73btTNdT7vkYPJjltpZjxgzg3XdN/3nkdrUDJ46l9YVcG5D2TRoNMw9du9b8+BypgM7P6c472UBaKbiWmpkdOWURJySEle3IM8RnAW0pBy2tPlxhVu0IfHZHOijVaOpeoLf3DEhjB8j1dUp9cGCgfCwDbv4qZwY7cyawapV8yjl7bgoTJii7T9SVkCvXBx08WLuyeB2l7ZAQ9n4aP960TKmP58IkLycgwKREu+kmlt/9s8+YgvS//0xuUgATiC0V3mquGT9mq1bMJFpuZnD0aPbMW5o2A9YCkrStV1Wpqws3yeezvT//zLI3WO4TFmZ9v4YPZ9/9+rFxTm0CAkrN7uXq+eOP5rFF5BgwwBQIy7IOcgFE27Vj4zM18RKCgkxKkZgY83eA5fVwNLZJly5MWeIozZrJH8tSAdfYcXoIQhmXKKWpFsuzKaX3AggAEA0ggFI6XC6XriuglMYDuAZAPICfCSFmoUEIIa0IIZ8BSAIwp2Y7QQOh9axAh/FrcHTJQzBUW0/pWJrfqqExzBzICakcacc1bJi80MYF9cuXTRpvLsjKnd+CBUyg5S8CQoBXXrEdtZn7yskN4s+cMY+0rGZmWYkVK1hgFHudslx05i1b2IBSmuLJHnImfRzpi03u2owezXxAd+wwDQxd7WPEX1Zq/dR+/FF+ed++pvvStav5zOPYseaDkqlT2cBASXAICWEB2WzhqLmyNKqr3ADTlum4nJDmKpSEfSWzaunzwQUwpRksgPnZSQP9PP88O5/Zs6233blTvoxt25TN5KRIo7oTwgbL48axY+3da+5CIScITpkCJCZa529eu5btrzRgPXHC9Fvp+Vi/HsjNddzMbuNG4PBh81QiljMmXIGydy+b0dm1Sz49nCNs2qQ+Smzv3qw99O5tWnbjjXXz7jlxAkitGVVFR7No/MnJ8tvOny8ffXbuXGYWa8tvV45u3dg5Sf0BufVI9+5MaP7ySxbQ77vvTNuoEbqU2oW9fW++WX65vVkqfjw1rhmPPgo8+6z97QDr2TfupgIoK/e4sCFnau7pyfr8li3Ztbj+eutrws+F3xc1ivi2bVlbsmXOKy1bGscjPt5aSSd9J1sqe5TYuhXYs0c5ZocUXg5/T/NZbVv3We7527HDZB59xx3Mp/j9962v6fXXs3OyF0wwKMh8ZtOWf7Z0zGNPyS/dp3t31jakgqRGY7oGABsPWqap5G4HAQGsnjExrCw12SSUiI2VfyabgsWFFKeFXELI8JpPkNx6SqmOUppOKXWhV6U8lNIsANcC+APAW4SQpYSQfoSQnwCcA/AkgE0ABlJKR9ooSlDHxAzeA3efcpz+63bZ9VzzJxfRkRD5BOYNEYXVEfgA0ceHdcC2iI4G3niDdd5cayjXkfv6MkFGbcfz1FOmGVaeeoRHwFSLVmsyNbX1YvP2ZkKYlFdeYZrFjAzTzNry5db7jpR5Oi3PcfBg07LPPpN/4fP6SfcdONAU/ZPv26sX8NBD7IXOZ7H4IP6LL9jg++OPgfx8uTNVB6+D2sG/h4dpUK90nS2vyU03mWb1MzNNM/r2tK+2okjzY3Ph35Yvbd++yiaOlqh9gUqXffKJ6beaaMFSBYrSYEPJWkQq5LZqxYJA2RsISc/Xy4vNWskpcYYOlT+u3Cyx3HMvnQ2Ttic3N9sDXyUffoAJSuPGsf2VhDbL59mSgADbigBbeHkx4VFqxs6Dv/Tpw2a71q1j5osDB5r6xfPn5WeVuVD2yCO2j3vDDawvks7OcJ9UJXjb+PlnVm+1Qi6P3K+Grl3NzR2ffdY6Uq4UOXeeZs2Ab76xP/t18qS5H7AcM2YwZcADD5iWtW/P+k1HUOoTpLOWcvz9N1Bebr3cXn8qp9SUpnCSsmiReR9jCX83uLubnl9+vaWm7Ep1euEF1lbHjbNZZUV4uddey/xH5axT5OJfdO3K2oIaevZkLlxffCEvJM2ezSzdKDUvk1+HS5es24SPj7mbky14++D34d13mSJWbkxgi2HDzBV8cXHWZuEAy5yhNH767TemrJSjbVvl51FufMqpTST24GDWtwYGyo9z3dyYUNuuHftERDDBV62ADTiXKaQx4wpz5W1ggaQahXxPKa2klN4N5it8H4D9ACYA+AlAF0rpHZTSAw1ZRwEw8h02SkveMVx2ff/+zK9R6tPHO3gfHyZ0cOHohhuAr75ipj6NGV9fZuqlNIujhL0ImWqhFPhcYhk+cyabMbblO8gHR08+aVqm0TCN6NNP2w7gJcfbb7MBVUQEGzBR6lhHDLAZKEqZ5pyb2kZGym978CATrC0HVt9/zwbNr79uHZjJUsh98kk2+J450zl/FFszmAMHymu577yT+SJZmlgrtYWbb2ZmbvPmmQ/uZs2ynWfQVj5gqU8uwK6ZXDqJS5es/VWl190yqrfa6OJSpLMsaoJvPfwwG5AlJpoGB+PHm5uQBQQwwSnewraH9ze8PURGyg/QV69mpneA/DlJ/fCkEMJmN5QGUZzPPrN+Rnx8TG1ebkZV6ZnibZ27FvBzjIoyVxpIBRlOUZFyHW+5hcUUKChQ3sae35mcWffIkaxNxMWx/j4kRD4vsJxJIL9X9oQnzqlTpnsk177lypYG+lODvXKdwZmAgl26MEVYt27Ws4x79zIFn5sbE6qcdUGQ7i+1MHn1VWWTWY6Xl/nMMWA/4vGOHey9JzUxlgrL0mwB/Nrp9fLRhPl6f3+TqS4XDh5+2LSdXDYFgF3DceNqPxvGn1dClN0gXJHJ4K23zN/5Urh5ttxygF2P774zdzVyBMv3jaenfZcaR5AGFNy40Xa/NHGi47Fe7CmUapNKUaNh47T27U31txSkPT2dezalYyhHJz4aM64QcgsBFFBKnZjjcB2EkPGEkJ0AbgULclUJluboPUrpGZs7C+oFr6B8RPc/gOxTnVFVLD/9ajBYmxf17s3Mr1asYC+7CRPYoOHbb9ngTfqAS807GguEsMGC3MtT7f6uRKOxn4icD5y4xpYHKQkJYabSrtb+vfaacj7GVauYIkMatZPPfCm9qPr0kU/34uvLBs1vvGEtJPDylcyfk5Lkl9vzH7cl3I0aBaSkyO9nSyssR/v21iayXl7KkZrtYTno0GhMA6AhQ9jx3nuPXTfL4BvSe8lNQrnJl73rxWexldq95fLFi+W38fRk/QEfIHbtago2BzDFxXXXWc9acI27Lf/U3r2ZefHdd7P/cgM63r7k1g0aZD/VzuOPmw/Kt2wxr5fcrJHSYGfsWPZM8+d+7FimADl61HyfwYNN6c84vD+Wnsfq1ey6/f476wts9VE8xyjfn6fouvdeZiHhTP8m9/xLrTjUKAhbt2bmwI88wkyQbWEZ8M+RHNJ1xfTpbBbdlsuMPeLjmZWNlIEDnQvCaAtp3ApCWFt+7z1TG5VTKlr2Gxs32s6R3qYNs2DiwQUB89gLchYhGo15HwGw9jpqFOtHv/7a1I9xoYWbeFPqfPo5JaRCrhxc4Fd6R9UlloofHuvAkSCJgOkcbb0fbPUVCQls1l8J/k564AH7z7naY3L69LF+X7u5mQu+rhqb+vrWLqKzVMju18/0LEj7SEej8TdmXCHkngfgTwiph8xRyhBC7ieExAP4C0BHAK8DaAVgFAB/sBRDw2wUIagHKAXajv4HALD2UZPt1rZt5tojuQ6OECYE8QGahwcz/7JMVXDgALB/v4sr3oA0pK8xN8HiA/66rsvcuWwmUo7evZnPknRQP306m0F0RcoTzvXXA/v2KZvDSn2rw8JMAzV710buJcnPlc+IO6rl5WWmplpHBZbj1CkWMMbSF9MWcj65ffsy38lVq5hfnjTaqJRHH2WKib172fP9998mE2q1EWPVCLllZSZBc/5806y4dLaD37cWLcwFOiUhls9oKfm3X76szipDo2GKFlf0SQEBJrO9zz5jdVSayZo/3xSZVAk3N3YfHQkmkplpiph+221MMFKj7OrTx9zP87nnmP/2kiWuidhp+R6w50su165atGDKEnd3NnOtNHvNZyD5vSDEvt/ryZO21ztLs2bMH7qpBYbhWCpmLl9myo+DB1m6M87dd5ub+7ZurU4IX7rUFFyoVy8mHPP8w4cPM1NlKTwFHGf5clbHN99kgsHvvzPTYDVmwJcuMasNZ7Flmr1li8kaxVZKwrrCUvHjqJWDUjmO0qmTsv82YFJG8Dgm9rD1XrfsQzQa+VSG3M3Dzc21KYG8vdnHlhuDJZaCcfPmbNxR2wjMjR1XCLm/AHAHUEe6K9UsBeAH4BkArSilb1FK8ymlewAMApAD4F9CiIwhlqC+2L8faDlsByqL/JF20BTy7tprmWDBcSYlwzXXNN0XvRx88NYQ2rW//wa2b3ddLltXQ4jj+RjVMGCAfc2tuzuQlmYaoKu9P9KXJk8PwXMXvvEG+5Yza3ztNWthnpcVFaUuZU3nzixgjC0z8/XrzYP/WM7kAuzazJplO8AZ3+7HH02myjffbJrtnTSJzepJzSMdGZwRwgbA//7LXvT+/ux6zJjBZjIOHzYPEDRpEkt19eST5gON6mr58m+7jZncKc2AR0fLKyU2bmS5vaW88opjqaPkuHTJPI/jPfewWTclH68ZM+RNe50lMFC9fx/AZnSk7UlKr16uswaRpi0DrIXcnTvNfWKLi9lHCX9/5WBF3G9Z6jO7eDFTAChhK7+4I8THW5vWN0X4TCmf+VLqb/v1M/eR9ve3HdlbCR8f5lNPKTPNvPFGU+Th3r3l89FKXUgshZegIPWmwS1bqvdJtYVcX8wZOdK+ZVZdYsvf3xbz5pkr42ydoxJqBdbawN9PtvpvNefapYv9mAaOotWyMi0tqGzBhW4+bvL0ZM+gLV9hW0E4t27dCkIIZs+ejQMHDuDWW29FSEgICCE4WdeaPRW4In7oAgB3APicEJJLKV3vgjJrwwMAfqaUWs0PUEoTCSEDwQJSfU8I6UApnWNVgqDOSU4GWg7fjpTdQ0EN5iqtkBAWsGDJkqabGL4umD+fvcCGDq3/YwcFsRfIoUPsf2OIYN0YOHiQaUDd3Ji5ZU6OfEoNKVyYsvXCmDWL3WdpygPO3LnsU9eMGcM+3OTZWc26Eu3aMX/Z4cOZcPDqq0wQ5dE7Ldvab7+Z/3/iCWUNtlZrLuAC7Dx4mg+pkGtrptiR6N4cR0zgHKEulDlq2LGDzZTXFmnUWcB1MQYssRykvfwysyTgFkJDh7IPD4rjiNXErbeqO354OBOWvvnGtHzoUHMFrrM4EzG1roiMtDZztoX03m/dymb0XZXDWadzXWT8t94yKRwbQ1RZNaa8nIEDzf2E65ravicsXWvUnKNU+XDsWN26p/Xty6I0K1krqaUxTRTwGWA1dOnC+jZLxS3n8GGWMOfEiRNYuHAhxo4di2nTpiE5ORmdamNP7WJc0RW8DGAHgO4A1hBCTgLYDSALgKJBGqXURhgUx6GU/mRnfSEh5CYAXwN4mRDSjlJ6jyvrILDNSsP9iL/zJ4QBOL5squw2PB2EvRmixsaffzo2u+EI3t515+Ojlsbwgm9MSEP4a7XqzOX69gX++Ye9kD/6yNxHTIqatBB1ISj89pt8wIm6Eko4c+awFE4zZrBjhYUB2dnWx7MMAKImHYUS/Jzatm1cg4/GiFI7rS1z5zKh0VFfc0eZMMF+ZGW1/Pmn+m0XLWKBEPlg3dFAg02Rs2dZ4LHaEBFhmlF1Ba40B506lfluVlU1jncg76vUCJJ799ZtXSyxvD4tW7L3pL2YA0rl2DrHAQOYu8bYsbXLAesors7HHR7OFLyNjc6dWT5mKbzNxcbKu29wIXfXrl3Yvn07BqoZwNQjrhBy3wBAYYqu3A2Amkn5Wgu5hJAXARwBsINSKhNUXh5KqQ7AI4SQ8wBkwtEIXE0xzcBSw43ItEhLHP+TTOx7sBDvPXrUPsR+Q6FG09+UmTKFmS7PmtXQNWm6vPEGU1Z07eq82RI3z3r0UaerZUSaVglgAWD27mUBuH76yTG/H0ewzDP699/mpn3vvstmx7g5uKcnG1A7M+jks0bSPKD1walTrs+/3NQYM8YUmbYuUZoZ9PdXDm5X18e+UrFl2t3U0WgczwBQV2zcyPpiR4M51QfjxrGAXLwduLsziydHUatUnT7d8bLrg/WGZ5FOj9reqMasfKfKuBSuoLwTUwAd0wPNSS+M1XxqtY2vr7KFS7Nm8hM5XMhdsGBBoxNwAdcIuTvAhNz65F2wWeJwAKqFXA6l9H1CyDmX10pgxSb6glHAbY7eaH/qe9x3H1CSId9La7Umk0JB4yEkhJmVCWqPm5vzfpmcli3r3nScp/IwGNhsnrPmWmrh/kX8/B591FyYP3SI+Q07Aw/WVN8+bEqzl+fOOWbueaXj4+P8TK/S82ErFZLgyoC7PV1ptG9fPy4rtWHhQpa/215KJ3vceSfzb3eFD7PAhFrzZEcoLS3F2bNnER4ejqlTp9rcNiEhAb1790bnzp1x5MgR11dGAaeFXErpCBfUw1Hyao5d67RFlNI/XFcdgRxZ9BSO0mXoRu7CjeQDBJPW2JYLZB1v6JoJBAK1aDT1J+ACJm2xUn7Cbt2c90u89lrg558bjwVGu3b288heTZSW1m6/N94wBW+zFUilPmjZkkUHFtQNJ06wqPJyfPcdE3IbQ2qnqwV3d9fEDRg9uunE/mjWjAWvkwagk5shbVAoZZpqV9rx13Ds2DEYDAaMGzcOGjvmK8888wxmzpyJjz76CFVVVfBwdf5JBZqqUc0xAAGEkDrQTQhcxVLDaAAwCrgAi+wKMG2fs7MxNklNZVGTjh6tw4NYcOQIs9cBWD4Xnqdl1SrzRJcNDaVAbi6L0PHii7ZDgtY1xcWNO1To2rXKCWybMseOsbDEShQX204Q6ypWrABKSswWNW/OfPw++aTuDksIi05sU7tNKZOY0tJY3pK33rJdKKXW12zRInYwpRDOtaWqytrm97//2Anl5bn2WGqorGwUEt3rr7MgWQkJEl/r7dtZPpp65uLFBrokWVnsAnD27GEhbB1tgzodi96Vm+va+rmIrl1tB3nLymJB05ocixfXT8VPnGBh3x2VKC2dNpsKBw64NgIcmNzYtm0tFWpFRVbvPlksE5Y7SlISG5s6gsHAxgB2jstNlQfY0YKvWrUKer0ezz//PKqrq3Hq1CnH6uMMlNIm9wEwGYABwEMNXZeG+vTt25c2Zi4bDtE5OtA1+ulmy3/+mT2x//3H/vfpQ2l/7KN0yRLTRtu3U3rhgul/URGlmzdTajCw/5s2UXrwoPyB9++n9MEHebdA6YQJlJ47x35v20apTkdpaSmlf/xB6T//mPZLTaX05Em2nUbDttm0ia3LzaVUq2Xrvv6alZeTwz6pqazc8nK2H0BpcjKlISHs965d7Pvxx9mxKysp1espTU+ndN8+Sisq2HEMBkrj49m2XbpQmp9vOvc336T05ZcpXbuW0oQESl99ldLqarbeYGDnkplpuj6VlaycL75g/6uqTOsopXT8eLb+++/Z9y23UHr8OPv9zDOU/v03K4Pzzz+U3n8/K0NajhwlJez8v/2W0rQ0Sk+fZsuWLaM0KYnSS5fYeVPKrsHYsey45eWUHjlC6Q03UJqdzZYtXUrp4cOUNm9O6YIFlL70EqWFhWxfvZ7VmXPgAKXTp7P6HTxI6QcfsHIs+f13SleutH0OnNRUVo+ICNOyS5co/fFHdu+/+opdl99+M63fuZPV+aefKF2zhi0zGNj9LClh5b34ovWx/vqL0oIC82X79rHrZ4nBQOmWLdb3Qqcz/79hA2sXnH79KH32WfabPx8VFex/dTWlKSmsjtOns3W3386ueV4eO8fly9k1uXzZ/DilpayN6fXseNL7QimlixZRevfd7NpRSum771Lq7c2uE0DpPfdYn+Pnn1PavTulb7xhvY4f89dfbbfHv/6i9NFH2e958yjds4fSjAz2LJaVsXN96ilKs7LYtc/LY9fwzBm2z59/svpFRpquF6WsHd9xB2uzUl54gW1jMFBaXMzOme/HOzwlzp+nNC6O/S4uZu03JYXSxERKp02zPlZAAKWxsawtAuw8eZ+zerXtY1mSlmbqTyhl/WBVlfV2773H+iC5az56tOn6UErpd9+x5zc9nT0r9jh4kPWzej2lnTtT+tFHlG7cyI717ruUeniwviI7m9JZs8z7gb/+Yt9Szp9n95G3ZT8/dn8tKS+XPx+dznz733+n9Mknrbd77jlW/pgxlP77r/L57drF2t38+awvS0pi/eTy5ZQePcp+O8LRo5SeOsV+795Nqbs7ay+UmtrcnXdSevas6f/TT5v258+99Hxzc9nvzz9nff6aNWy/u+9m11Pav9h7D6Slme4JQGmHDubvW0rZc5eWxvqTL79kz5slJ0+a3kV5eabn4ORJSo8dY3XS6Si9eJHeOqaS+qCEteVbbmHHffJJ9qxLqa5m1+K119j/nBx27r//zq7lnDmmd78lHTtSOmOG6f+pU6bzNBjY8bZvZ8fIy2Nt+L77rMuRXr+9eynt0YPVYe5c0/3auZP1mSdPsvKkzxG/vtddR+nw4axfk+P8edbPvfeeaUyxaRNrb926seNcvMj6GGn7PXqUjblOnKB04UK2rKjIVDfLPmbVKkpXrGC/d+9m/fPPP5veXydOsP1uvpm1J8t3XUEBW//LL6ZlRUWsTRcUsPIs96GUnVtiovn11OnY+9Fye2kfzjlzxtSXHDnC+lN+3U6eNN9Wr2fjL0rpKf7sUcruTVaW+firpMS8T5Xj4EH2OXTIvH/n++r1rB8+ccL0TiorY/tI++eiInZ9pOh0rL0kJ5uOIx3T6XSsb7B8r1RXs3aYns72SU1l56bAgw8+SAHQQ4cOKW5TVlZGY2Nj6eHDhymllEZFRdEl0vG+HcyutQ0AHKJy8qLcwqbwAbARQC6Arg1dFwfr7Q8WrCseQAmAQgAHATwHwENtOY1ZyK0wFNE5OtA5OtBCg/mAeNUqSvviIE2d/g6lr79Oqx+fbup87rnH9FvNZ+JE9pK57jr239/fsf35Z8MG5XWvvqq8ztPTftm+vrWr0xNPsIGII/vw69Czp2nZJ5/U7vgAE7IMBtP/2bPZ94ABbBD+5JOUhoZS44DqiSfUl71nj/K6adNs79umjfn/YcNMvwcNMl+Xl8cGCgCl33xjWr55s0nIsvy4ubEBkHTZ/v3m19Xy8+GHlI4cab38yBHT7/Xrzet1551M4OXLunVjAvj06WwwzJf/+ScbRH/2GRMO+fL77mNCR+vWTNgGKI2Olq/f55+bfltuM2OG7es9ebL1sp072Yv3jz+U93voIXVtwd2dnePp06zNp6VZb/PDD+x70iR2v0eMYP9/+onSt99mSqA+fdiA7Pnn2SBH6XwBprQaMMB6eefO6ttwy5am+2PZZrk2T/rp0IEJwj16sP9vvsnqLN2GP8NKn+XL7dfLw4N9338/u0cAU5wBTHB56y1Kp041KV2kxx440LzNZ2VZlx8Wxr5jYlibLS83rYuNpXTmTOt92rY1/R49mik1y8uZIkPt9ZZ+nniCDeykyzIz2YDw5puV9/v8c9beb7+dnZuvr0lAP36cKUhjY20f+6+/2Lvm+uut13EOHWIKpXHjWJt05NxWr2YDzEWLTG3+1VdNz760bUs/9vpNgNJWrSglhP2+6Sb2nZfHBHCAKW/4tlwBKf3s3s2+IyNZXyVt77fdxt6lW7datxXpZ8sWSvv3l6/fnj3s+l24IN+O5D41CmiDm5vyNu+8w569U6eYQM2XR0Up78ORewYA877C8p3Uu7f19oWF7Jnj1/X2283XP/CA/XOVKvCVPvPmsfZiWT7A+lpb+yq1VS4QSz9yfRVXoHfsaFo2Zw5rw5bbtmvHvnv0MB9L7dlje9xUWsoEO6nQbe8jLW/2bNb3JySYjs+FbKVP585m75FTmzczIZQr27kgefmy6bf0k57O3m16PftOTDRff/Qoa2tVVez/8eNsMkW6jVRgPXiQ1ZkLvVKBOSuLCehy9dDrmYLt2DHTsnPnmPJQbnvLc7BQGPXs2ZN6eHjQSqkAbcFrr71G77//fuP/G264gT4tVbjZwVkhl7B1TQtCyN8ALgB4GIAWwAsAvqYsenKjhRDSCsA2AK1rFpWB1d+z5v8RANdTFb7G/fr1o4d48tJGhJ7q8JWhLzJxHNHoh8e15uH1Vv6iw4RJDewsJRAIBAKBQGCLyZNNbkcCQQ0JGzagc2ioawv18XEuIXl90aMH4OGByspK+Pv7o0ePHlCSRZKSktCnTx8cO3YMMTXpGWbOnIm4uDhs375d1eESEhLQWUUUQkJIHKW0n+XypprMYDxgjOhMACwA8CYhZB2ArQAOAzjRmIReQogWwBowATcdwAOU0v8IIRoAEwEsBtAbwHIAYxuqns6yi36ITBxHBLpjkma11Xr37LQGqJVAIBAIBAKBAwgBV1BfNAUBF2DxDdq0wYkTJ1BdXY2+NqK7zZgxA/n5+YiNjTUuo5TCz8+PzbLWQ/Lpphp46hMwYZbPeBIAQWC+uosBxAEoJoQcIoR8Qwh5nBDSv0FqamIqAJ6V8U5K6X8AQCk1UEp/BTCtZt0YQkgdZ/GrGyiliKOLoYEWT2mOIZC0sNrGM0tlEB+ewPLJJ03LRoxg/+Pi1JWxYgXw6afs90MPsRwukZGm9cXF5vH4hw5lmltLEhKYkcqUKaZlSmEEx44Fli83X7ZqlSkfi5SpU1mQmNOnWQSD558Hpk2z3k6Jt94yhZnlYWLXrjUdk8fgDwhgCfZ4olEA2LDBOurOO++oO+6zz5r/X7rU9PuBB1gSyhdeMN+mbVugf3/5BMgPPsi+37ZIXX38OLvXxcXq6tW+ven37t3m6379lX23bm2/nNdfZ98eHo5nsr/rLse25wwebJ1oUvJiMGP2bOVy3n6bBS3hebgs2yJn0iTWVi3hoWmVsGwzHTpIIvxIuOkm5TL+/BN44gnbx+Fs2AD06WP6byvJbEiI8rqOHe0fKzCQlREUxJJRhoWx5aNGsfvDGTHC9Ql2+UDBxwf4/ntTMLZrr1Xehyfv3L7dFP5aKQLKbbex/B5qsXVctfSzUqqrIzzcse1btardcWwd35HcUj//LJ9AUo7ffwc++4z9vniRBTWzx+TJwPz51sulfZN0sLhhAwscM2+eaVlpKaDXA0OGsP88yejMmbaP/dtv9utnj/37nS/D1Vi+a6Q4knxcGlHWMpG4I0GOLPvVxYuB666zHx46KAiYM8d2+PsRI0xtJTLS/N2YkMDCUMu1d8s+etAg1hfffz/77+hzCrB307Fjpv+jRilv+/LLwB9/mPra4cNt5yeaM0c+KFZEhO0IZVLmzrWfe0rpvC3f3127Ai2sx8Bm9Omj7t0kxdGoxM2bm4d/toWaPEMBAQCAvn37glKKRYsWyW62adMmbN++HXFxcTh69Kjx8+uvv6KoqAgXL15UewbOIWfD3JQ+AFoCuBXMz/UvAMlgQan4R1/z0TVwPXk+4S0K6wmAxJptltorrzH65CYattI5OtAj+mWK22yZ9gulAM34YQOl69YxH47OnZmfwqFDzD/o22+tAzMkJZkXtHMnCzyzahXzJ/jhB7bN99+zdTzIghx6PfPn4JSUsGNz4uOZT4Q0+BUnPd28LufOMT+w8HDmr8GD7hQXWzv1Dx3K/DALCuSDIlHKHP9Pn2aBAfLzWcASvZ758xgMzC+jsNAUaIAHGFCDTseubUkJ+19WxvxfPvuM1csSg4H5PQMsGBXA/PQcobzcOigMJz/fOjgE95vlQZs4PPDE/v0sOMVtt7G28t9/lLZoYbqfxcWm61FWxvyStm0zL2vjRuYX9sknzC9XDt7+DAbmp/bXX+zeJCay5Tk57H98PAtSceCA+T47drA68SAkBgPz/evZk12TVq1YAJCyMlNQCUpZAJHOnZk/XnU187m9eJGtO3rU5LtTUcF8Ky9eZPeV31OOXm96BpKTmV9UdTXb/sAB8/NcvZoF8+GBZXJy2HanTjHfnbfeYvdE+sxQyvyP+DknJzNf2jVrzM9HrzcFj7EMLGMwmJ4RvZ4F5DlyhNL336dG3yxK2b3/7DNTgK+//mJ+VFu3srZQWWlqR2vXsmAln3/O/HQNBhZgg1IW2Kl3b+az9NlnJn8ry8A00rorBYwqKWF+3SUlLDgHD6qVksKCx+zebdquuJgd5/Bh5vv41ltsHQ/oYS+Iz8aNpqBqO3aw4/H9+L4VFab7cfYs67uWLGG+5dK2YTCw8zUY2H2yDFYm3e78edbeLO+7lKoqU8AiOcrL2XX/7z/W/n7/3eSHfs89rOxTp0z3iDNzJgsCYzCwPovf35wcFisgLY351EmDyJ08aYrL0KED62N5XIFDh9i7ggc8SUlRPnfOH3+w633kCPPzlga44cHkiovNj79zpyn4HS//+HG2rQ2/NVm++475YHN0OvYu3LPHVHZlpel3SQl7NqTtKSPDOiiNFIPB1A6PHGHPb3ExW84DBGZksHPT6Zjfc24uuzbS5+bgQXaP+bHLy03tVHqs3Fx2vLFj2blIOXPGFCxRyoUL7J144QIr49dfmU/wwoWsbXBWrWLP9RtvsPfE77+zumdlsWc1NZX1dTyg0FtvsetFKav7vHmmsgoL2ef775mv8zvvsOXJyaZniLfZ0lJTIC+Dgfkn8/d7aSnbdu9eVp/CQvNxgU5n8o/nx7XkzBn2LP7zD3smPv2ULT9/Xv7d+tNP7JmKjaX09dfl+xdp/TkGg3kAS0rZs8vfP5yKCtN7SMqSJSzGQGGhqc8YNIi908+ds96eU1hoasMZGWzMM2KEqQy93rwtGQxs3dixrG9MTTV/tsaPZ37XR4+yNi0lJ4eNZUpL2TjOYGDHP3jQFLhNSkYG6+Pi4lgda673qVOn2DruF6vXm64bDxglrW92NusDKipY3fPzzceovM+QC6jKx32XL7NydTpTeysuZuUZDKY68PcqD17Fyzh9mvVhRUVsPW9rVVXsGkrrwwN+paSw86yuNr3jVVBVVUU7depE3377bat1hYWFFABdqTIA6FXpk2sPQkgImOkv//QB0J5S2iDm2YQQHwDFYDPnz1NKP1TY7ksATwDIoJQ2t1VmY/TJ3Wh4Hnvpp3hZUwAPYjG7QynTkD3yCAAgPaEAzTsFMu2yXu+4dkpQP+j1QEEBm6koKjJq8QSCJk9qKut3+IytoH6glOXWsTXz7kzZgPnMpkAgELgQtX6iDnEF9V0fffQR3n//fVy8eBH+ljPcAJo3b45HHnkEb9lLy4er1yfXJpTSPACbaz4AgAbOqdsZJtPwEza24+siCSEhNefRZIhPXA1/fSsceHc5NPpqlmdPr0Ng1nl03rkIGmowbqsNqTGd1WrrJEm1wEVotSZTPCHgCq4kHDFJFbgOQupGwOVlCwQCQVPjCuq7Zs2ahVmzZimuT09Pr7e6XJFCrhyU0vIGPLzUID7VxnbSdVEAzIRcQshjAB4DgJZKPqENhI5WojTiPAYvBIYue0x2myyEYTVux0/+T2KTtXJHIBAIBAKBQCAQCJzmqhFyGxipSGcrhJp0nZUYSCn9BsA3ADNXdk3VXIMb8cT9R+Ohv16Py6NCQdzdoPFwA9zYN/H1AdEQTNAAU3xNMVMEAoFAIBAIBAKBwJUIIVfgMtoM69bQVRAIBAKBQCAQCARXOU01hVBTQ5oDRSbfhuw6lXlTBAKBQCAQCAQCgUDAETO59YM0GV40gOMK20kjodhMoBcXF5dDCLnkbMXqgFAAOQ1dCcEVi2hfgrpEtC9BXSLal6CuEW2sjvn333+76/V6XUPXoyHQ6/VuWq223s49IyPDrUuXLvEqNpVNmC6E3PohASxnrwZANwAbFLbj9r4Z9iIrU0obZd4LQsghuTDeAoErEO1LUJeI9iWoS0T7EtQ1oo3VPceOHUvq2rVrDrmCIiKr5cSJE527deuWUB/HopRCr9eHOtOehblyPUApLQOwu+bvTXLbEPa0jK75u6k+6iUQCAQCgUAgEAjUQQjJr6qqcm/oelzpVFVVuRNC8p0pQwi59cfSmu+RhJABMusnAmhT83tZ/VRJIBAIBAKBQCAQqMFgMGwoKCgQiTDrmIKCAn+DwaBk+aoKIeTWH0sBxAMgAP4ghFwPAIQQDSFkIoDFNdttoJRubqA6uoJvGroCgisa0b4EdYloX4K6RLQvQV0j2lgdo9frv8nMzCzIzMwMqaysdKe0UWX0rFNCQ0Oz67J8SikqKyvdMzMzQzIzMwv0er1T7ZlcTTenoSGEtAawFUDrmkVlYIoGnjX2CIDrKaVOTc8LBAKBQCAQCAQC1xMXF9daq9U+ptFoxlBKgxu6PlcShJB8g8GwQa/Xf9O3b98kp8oSQm79QgjxBzALwB0AYsECUp0F8DOAzyilVQ1YPYFAIBAIBAKBQCBo0gghVyAQCAQCgUAgEAgEVwzCJ1cgEAgEAoFAIBAIBFcMQsgVCAQCgUAgEAgEAsEVgxByBQKBQCAQCAQCgUBwxSCEXIFAIBAIBAKBQCAQXDEIIVcgEAgEAoFAIBAIBFcMQsgVCAQCgUAgEAgEAsEVgxByBQKBQCAQCAQCgUBwxSCEXIFAIBAIBAKBQCAQXDEIIVcgEAgEAoFAIBAIBFcMQsgVCAQCgUAgEAgEAsEVgxByBQKBQCAQCAQCgUBwxSCEXIFAIBAIBAKBQCAQXDEIIVcgEAgEAoFAIBAIBFcMQsgVCAQCgUAgEAgEAsEVgxByBQKBQCAQCAQCgUBwxSCEXIFAIBAIBAKBQCAQXDEIIVcgEAgEAoFAIBAIBFcMQsgVCAQCgUAgEAgEAsEVgxByBQKBQCAQCAQCgUBwxSCEXIFAIBAIBAKBQCAQXDE4LOQSxlRCyEZCSDohpJIQorfx0dVFxQUCgUAgEAgEAoFAILDEzZGNCSGeANYBGAmA1EmNBAKBQCAQCAQCgUAgqCUOCbkAXgBwXc3vVQD+ApAGQMzWCgQCgUAgEAgEAoGgwSGUUvUbE3IKQEcAb1JK59ZZrQQCgUAgEAgEAoFAIKgFjgq55QDcAQRRSkvqrFYCgUAgEAgEAoFAIBDUAkfNlYsBaIWAKxAIBAKBQCAQCASCxoij0ZUPAggkhITURWUEAoFAIBAIBAKBQCBwBkeF3PlgUZVn1EFdBAKBQCAQCAQCgUAgcAqHhFxK6WawCMsvEkJeJYT41E21BAKBQCAQCAQCgUAgcBxHA09tqfnZC0AggEoAJ8F8dZWglNLra1tBgUAgEAgEAoFAIBAI1OKokGuoxTEopVRbi/0EAoFAIBAIBAKBQCBwCEejK4vcuAKBQCAQCAQCgUAgaLQ4NJMrEAgEAoFAIBAIBAJBY8bR6MoCgUAgEAgEAoFAIBA0Whw1V24QCCG+AIYAGAQgCkAYAC8AuQCyASQA2E4pPdtglRQIBAKBQCAQCAQCQYNTa3NlQkgfAJMA9AMQDoCCCZwHAfxMKT3iVMUIcQcwAcBjYAKuNHgVqfm2rHwmgBUAFlNKzzhzfIFAIBAIBAKBQCAQND0cFnJrZlUXA7ibL7LYhBf4C4DHKKWlDpbvDeA5AE8DaCYpXwfgDIAcAHkAygGE1HxiwWZ3pcffAuAVSukBR47fVAgNDaWtW7du6GoIBAKBQCAQCAQCQYMQFxeXQykNs1zuaAohDYBNAEaCCZ/pYMLk5ZpNWtSsiwITNrcAuJGqPAghZCqAt2r2JwBOAPgZwA4AcZTSChv7tgYwAMCtAG4B4FNTh98APE8pTVF9ok2Afv360UOHDjV0NQQCgUAgEAgEAoGgQSCExFFK+1kud9Qn9wEA1wGoBptt/ZJSapY7t0YQfhzAJzXb3g9gmcrylwCoAvANgK8ppcfUVoxSmgQgCcCvhBAfMFPn58FmnBMAvKm2LIFAIBAIBAKBQCAQNE0cFXLvA5sdnU0p/Vxugxqh90tCiBbAAjDBWK2QuwjAu87OulJKywAsI4T8CGAiAHdnyhMIBAKBoDFSWlqK+Ph4lJaW4t9//8Xly5cxcOBABAUF4YYbbkBERERDV1HQAOh0OqSlpaG6uhqffvoptm7diuzsbDRr1gw6nQ4VFRVwc3NDixYt0Lx5c/j5+aFnz54wGAxITk5GZGQkKisrcebMGRQWFiIzMxMtWrRAu3bt4OHhgT59+iA4OBhlZWXQ6XTo2rUrWrduDTc3NxgMBlBKodVqUVhYiG3btuHs2bMoLS1FSEgIWrVqhZiYGJSWloIQgtzcXOTk5CA2Nhb9+/eHn59fQ1++KxqdTgdCCLZv347Lly+juLgYJSUl8PT0NPYXOTk50Ol08PX1hb+/PyorK5GUlIStW7ciJYUN0QkhoJSiZcuW0Gg06Ny5Mzp27IiQkBB07doVXbp0gaenp9mxi4uLkZ6eDkopDAYD9Ho9oqKiQAjBsWPHEB4ejo4dO0Kr1SI9PR0FBQWoqKhAcXExCgoK0Lp1a3Tt2hVardbqvAQCSxw1V84GEAggkFJabmdbbwCFAIoopaFO1VJghTBXFlyJVFZWws3NDVqtFsnJydixYwd27tyJoqIiUErh7++Pbt26oWPHjqCUoqSkBP7+/tBqtThy5Aj0ej0iIyOxb98+7N27F1FRURgyZAiqq6tRVlYGLy8v5OXloX379hgwYAD69+8PNzem69PpdCgrK0NFRQUIIcjJycHu3btx+PBhtGvXDoGBgfDx8UFRUREiIyPh7++P9PR048DNzc0NZWVlqK6uhoeHBzIzM9GuXTt0794d5eXlCAoKQmRkJCilIMQylEH9YzAYUFpaCn9/f7Pl2dnZyMjIwOHDh2EwGODr64vOnTsjNzcXer0ely5dQmBgIDp27Ijw8HCUlJQgPT0doaGh0Gq1KC0tRVZWFrKzs43X7OLFi9Dr9QgKCkJGRgZ27dqFzMxM4zUPDQ1FZGQk+vfvj06dOmHAgAGorq7G+vXrcfbsWfj7+6Nfv37QaDTQaDTo1q0bgoKC6uU6UkqRlpaGf//9F4WFhSguLsbx48dRVVWFsLAwaDQaxMTEwMPDA6mpqQgODkafPn3g5+eH3NxcHDhwAGlpaWjbti38/f1x8803o7KyEu7u7vD19YXBYEBZWRlyc3MBAMnJyYiPj8elS5dgMBiQnZ2N5s2bIzg4GCEhIYiJiUFRURHOnTuHiooKrFy5EmVlZQAArVaLZs2aISsrCwDg7u6ODh06oKysDAUFBfDy8kKLFi1ACEHfvn3Rr18/6HQ6uLm5obKyEn5+fggICIDBYICfnx9atGiBTp06GQezmZmZxkGvI9evMbT3xkhaWhq2bduG6upqEEKQlZUFvV4Pf39/+Pn5YdiwYcZ+JS8vDx4eHggLC0NMTAyKi4tx9OhReHl54fz58wCA6OhoHD9+HL/88gsOHTqEqqoqAKxdjBkzBlFRUcjJyYG7uzt8fHxQVVWFS5cuISsrC+np6SguLgYAeHp6orKyEgAQExNjbAvJyclITEw0CieWaDQaY/spLi5GUFAQKioqUF5uc7hoRWxsLLp37w5KKcrKyuDt7Y3Q0FB4enqioKAAHTp0MF4PnU4HjYZlw+TvhebNm6O8vNx4jgEBAYiOjkbbtm3h4+MDgPV/fD8pRUVFSE5ORtu2beHt7e1QvXU6HSorK+Hr64vS0lIUFxcjKysLSUlJaN68OXr16gV3d+s5F+m1TE9Px+bNm3H+/HlcvHgRZ86cQXh4OHJyctCpUyfExMQgICAA1dXV8PHxwbZt23Ds2DEQQuDm5gZ3d3dER0cjKioK/fv3R35+PhITE3Hx4kUUFxcjJycHly5dcui8pPTt29esT6isrERiYiK0Wi1Onz5tbEMA4Obmhs6dO6N9+/aIiorC+fPnsXXrVmPbUiIwMBDu7u7IycmRXe/n54f+/fujS5cu0Gg0xv7l4sWLMBgMCA8PByEEwcHBqK6uRnFxMXJzc3H+/HlER0cjOjoaISEhaN26NYKCgtC1a1d06tQJZWVlSEpKwvnz51FeXo7q6mr4+fnB398f/v7+cHNzQ25uLgoKCuDj44P+/fsjJiZGtn/j172srAzt27dHRESEYj9YWVlppQwQOIaSubKjQm4FgFJKaTOV2+cC8KWUeqk+iEAVQsgVNCV0Oh3OnDmDEydO4Pz58wgKCjIKLHzglJKSggMHDsDf3x9RUVE4ceIEAPbCCw0NRUlJCfR6veKLT4q/vz+GDh2KpKQkJCQkgBACb29vlJeXw9/fH0VFRcZtfX19UVVVherqatmy3N3dFdc5Srdu3XD69GkMGDAAHTt2BCEE7dq1Q35+PlJTU9GsWTNEREQgMjISVVVVSEtLQ3R0NJo3bw5KKVq0aIGCggLk5OQgNzcXAQEBCAkJAQAEBwcjPT0dWVlZyMzMRGFhIZKSklBaWorAwECEhIQgIyMDWq0W5eXl2LZtG8rLyxEZGYnOnTsDAC5cuIDk5GSXnKstmjdvjrZt2yI/Px9RUVEoKipCSkoK0tLSALBBOdf0azQaGAxmXjFwd3cHIQQdOnQwDloIIfD19UVkZCR8fX3RsWNHVFRUwNPTE0VFRUhMTESPHj0QEBCA4uJiREdHo0ePHgDYIJMP9HNzc5GXl4eEhATs2rULBw4cMBu4AUB4eDiaNWuGnJwc6PV65OXlAQA8PDxQXV0N6XuVD7b4NmrQaDSIioqCj48PfH19kZWVhcrKSuTn5xsHxIGBgfD29sawYcMwadIk+Pj4oEuXLmjRogUyMjKQmZmJn376CRcuXICfnx+CgoJQVlaGy5cvQ6fTYd++fUbh2BZeXl7w9vaGwWBAYWGh8fp7eHggKirKOHsXHx9vVDh5eXkhMzMT6enpcHd3x4gRI4wKJg8PD0RHR2PUqFFGITssLAxVVVWorKzEhQsXcPDgQcTGxgJgA9q0tDSj4F9cXAytVovY2FjExsaioqICeXl58PHxwYgRI9C7d29EREQYlVdZWVlISUmBu7s7vL290aZNG1WzQCUlJUYlhFarhcFgMA5S+TelFPv378eyZcuQlpaGa665Bh06dEC7du3Qs2dPKyEqKysLW7Zswd69e40KtNpkuPDz80Npaanivj179sQNN9yAdu3aobi4GLfffjvatm1rs0yDwYCsrCxotVqEhoaioKAAHh4e8PX1tdq2rKwMp06dQklJCby9vUEIwfHjx5GcnIzc3FwYDAaEhYUhLy8PhBBMnDgRXbt2RWBgIHJzc3Hx4kVkZGTA19cXOp0O4eHhCAwMxOnTp3Hs2DEcO3YMCQkJcHNzg4+PD8rLy5GZmYny8nJ4e3sjLS0Nfn5+0Ol08PDwAKUUer3ebnvWaDRo1aoVmjVrhqNHj6JZs2YIDg5GYWEhIiIicPnyZeMzptFo0L17d9xyyy3o2LEjMjIykJiYiNTUVKSmpsLb2xtlZWXIyMiAXq9HVVUVKioqUFZWhpiYGKSlpVkpAtzd3dGxY0cEBQVBp9PB29sbeXl5OHfunFXdufIsIiICxcXFCAgIQEpKCjIzM83KjY2NxeDBg6HVaqHT6VBeXo709HRcuHAB2dnZAICwsDB07NgRgYGBCAgIQPv27aHX69G9e3f07dvXKMQVFxcb+6ng4GB4eHigpKTEOMvL+1YlKKXIzs5Gbm4u4uPjcfToURw7dgwnTpxAVlYW2rdvjxEjRqB///7G9z8hBKmpqSgoKMDAgQORk5ODPXv2AGBKi+bNm8PLy8uogDt79iz27duH3bt3IykpCQAbX+j1erRr1w4GgwFJSUkIDAxEdnY2vLy8EBAQgMDAQERHRyMnJwfZ2dnIzs42E7a50O4oUVFRaNOmDQoKClBVVYWioiIUFRVZ3c/AwED06dMH/fr1Q8eOHeHp6Ym///4be/bsQWpqKpo3b46pU6ciLCzMqEgHgPPnzyMpKQnt27dH9+7dsW/fPlRXV8Pb2xtBQUHGsVRFRQV8fX3h5+eHDh06QKvVQqPR4MKFC/juu++QnJyMoUOHokePHvDx8UFAQABKSkoQEBCAjh07ol27dsjMzERKSgrKysqQnp6O8PBweHh4GJVFXbp0QbNmqsS/esdVQu5lAJEAoiilWXa2DQcLTJVBKY12sL4COwghV1AXUEpx+fJlFBYWwt3dHQEBAWjevDkyMjJQUlKCtm3borq62vgyTUtLQ1paGoqKilBSUmI0RUpLSzO+uIqLi3H58mXFAYhWq0Xr1q0RGRmJAQMGGAcRt912G0aNGoVu3bqZDUovXrxoNJfis2Wenp5o3749AgMDkZKSgtjYWHh4eAAAysvL4enpCY1GY5xVSktLw65du3Dy5EkUFxfD09MTPj4+8PHxgZeXFyilCAwMNGqt8/LyUFZWhtLSUnh4eCAvL884UxEQEIDS0lLo9Xr4+PjAw8MDZWVlCAoKwoULF3D69Gn4+/sjISEBGzduRK9evbBnzx5kZ2cbBwVarRbR0dEoKCgwE8CdQavVom3btggODkZBQYHxRcqvy/Dhw9GiRQucO3cOJ0+ehJubm9EUzNvbG2PGjAEhBBUVFTh+/DgiIyPh6emJFi1aoKioCKdOnUJhYSG8vLwQHh6O7OxsGAwGeHh4oGXLlggLC0NxcbFRoPT09ERpaSnCw8ONmnZLLl26hDNnzmDNmjXw9/fHnXfeie7duyMvLw/x8fHQarWoqqrC+vXrUVhYiLS0NBQWFiI5ORkajcbYBtXSqlUrlJaWmgmPHEIIWrVqZRQWRo4ciTZt2sDf3994DTmVlZWoqqqCr68vKioqcPToUVy+fBlt2rSBn58fOnbsaDQD3bRpE0JCQlBZWYmKigpotVr4+PgYlRUtW7ZE69atZTX7lZWVyMzMNA74uCBXG7jgoNFooNfrjfeHP08lJSU4d+4cTp06hcrKSlRXV6Nz586oqKhAfn4+KisrkZaWZpzZa9GiBdzd3WEwGFBRUYFmzZohOjoaubm52LdvHyilxpmV06dPGwVmtbi5uaFr165wc3MzDozT09Oh1WoRFBSE0tJSVFSw2JQ+Pj7QarWygo+7uzv8/PwQEhICf39/VFdXIzg4GOXl5SCEoKysDKmpqSgsLERQUBAKCgrMnnGuTHFzc4Onp6dxdjsmJganT582HsfDw8M409isWTPo9XpkZGQY1w0cOBCjRo3CzTffjICAAOP18fLyQnFxMTIyMrBjxw54enrC19cXISEhqKioQGZmJk6ePInQ0FD06dMHBoMBUVFR8PPzQ1JSEtq0aYOOHTvWul00BbKyshAUFGRUdgEwPl+enp5ISUkxPo+8v05LS8OpU6dw/PhxFBYWonfv3sjLy0N+fj4CAgJw+fJltGvXDmFhYejcuTPOnTuHLVu2YNeuXcbjBgUFISYmBs2bN0dlZaXx3aXVauHh4QGNRoOgoCAkJiaidevWaN68OcLCwhAVFYXU1FQcOXIEp06dQlFRkVHZGBISgsjISLRs2dIo+I4YMQL9+vWTnfWtqKgwztCXlZUpzhAaDAYkJiYiOjra4RnppoijFiMGgwEZGRkoKCjAyZMncfLkSfj7+yM2Nhbt2rUz9q/clLu4uBjV1dVGxUh+fr7RYuzy5cto1qwZPDw84O3tjYCAAERFRaF9+/bw8fEx9qOHDh3CsWPHjErz0NBQ3HjjjejYsSMOHjyItWvXWtXTzc0NISEhyMvLg06nA8AEZr1ej5KSElXnGhISgmuuuQY7d+5Updi0xbx58zB79mynyqgLXCXk/g7gDgCLKKVP2tn2KwDTAKyklN7lYH0FdhBCrsBZKKXYt28f4uPj4eXlhX///RebN29Genq62XYtWrRAWlqa0XS1tNR2VjA3Nzc0b94c0dHRCA4ORkBAACIiItC7d2/06dMH7du3R35+vnF2zt/fH4GBgXV2no0Zbo7HB8wAM5fLycmBh4cHIiIikJGRYZzhTE9PR0hICEJDQxESEoKCggIUFBSAUoqcnBxER0cjMjISYWFhRmHgaoP7bx08eBCBgYHQaDTQ6XTo0aMHDh06hKKiIoSHh2PHjh3YunUrWrdujbCwMKNvYrNmzRAUFIT27dvDy0sYIdUFOp0O8fHxiI+PN7oG+Pj4wNvbG8HBwRg8eDBSUlLg7e1tNNeOiIiw8tUsKSkxCpOVlZXYvn07EhIScPHiRQBM0cOVKtwE9+zZsygpKUFubi7y8/Ph7u6OvLw8BAUFwWAwGE26g4ODkZWVhZCQEJSWlhoVWHxwWV1djaKiIowYMQITJkxAQEAA8vPzcfnyZRw/ftw4mNVqtcjLy4Ner0fXrl0xYsQI9OnTxykFhaD+KCkpwaVLlxAVFYXg4OCGro6giVNZWYnU1FSUlJSgXbt2RvN5AEYriKCgIOOMPVfGVVRU4MyZM4iNjUVAQAAAGPugsrIy4zaFhYU4efIk9Ho9SktL0blzZ3Tp0gX+/v5GwZ4ry7jFRnx8PBITE9GiRQujBZGbmxuqq6uNArlGo8HRo0cxYsQIDBw4sEGunS1cJeSOAEsLRAGsADCXUnreYpt2AN4AMLlmu+sopdtrW3FHIIQ8wH9TStUGu2qSCCFX4Ajp6enYs2cPMjIykJWVhd9++w2JiYlGjTAABAQEYOzYsRgyZAgiIiJQXV2NjIwMbN26Fe3atUOnTp1w+PBhREVFwcvLC56enkbfn8DAQKOpTLNmzWR9nQQCgUAgEAgEAlfiEiG3pqD5AJ4FE2ABIAVAas3/GLBcuQDLczufUjqrlnV2GEKIoaYelFJ6RatJhZArUKK6uho7d+7En3/+ieTkZJw/fx4nT540rieEYPDgwRg6dChiYmIwZswYlJWVoU2bNmZaRYFAIBAIBAKBoDHjqjy5oJTOJIQkgs3WhgBoWfORkgvgDUrpF7Woq7NcffZ5gqua9PR0rFu3Dnv37sWZM2dw4MABVFdXw9fXF7GxsYiJicEDDzyAa6+9Fi1btoSXl5cwuxIIBAKBQCAQXLHUaraTUvo5IeRbADcA6AcgvGZVFoBDAP6llFa4pooOMbK2OxJCmgG4BcD1APoAaAV2fbLBzmkppXR1Lct+A8DrKjZtb2n+LRBYEh8fj7179+LgwYM4ePAg4uPjjWHz27Zti2eeeQY9evTAhAkTroqAEwKBQCAQCAQCgZRam/TWCLFraj6NAid9fzNgfj0qAFQDiK753EoI2QBgAqW0tuHJqgHYyiOhq2W5giuYjIwMHDx4EFu2bMG2bdtw9OhRACzE/zXXXIPb/9/encfZOd7/H3+9EwkhEiEaSxFSVCu0Fd8gaq+tKL4/W2whlrSV1pK0tJqkifKt2AVtUKJoqLaU2oklllpqCVpaYqdEQpImsn5+f1z3dI5j5szcZ5Yz5+T9fDzO477v677OdX8mjpn5zLXttx/77bcf/fv3XyYXGjIzMzMzK1TT81ZzWg54ArgauCsiXgOQ1Bc4HRgK7AH8Gji8zGc8GhE7tDRQq12zZs3itttu4+WXX2b27Nm888473HrrrSxatIiuXbvyzW9+k7POOosDDzyQ9ddf30mtmZmZmVkRJ7n1doqIKcWFEfE6cIykxaQtkQ6T9JOIeKu9A7TaNWPGDCZPnszYsWP58MMPkcQqq6xCz549Oe644zjkkEPYfPPNP7eFhpmZmZmZfVajSa6k32Sn70XET4vK8oiIGFpOcO2poQS3yJWkJBfSPGQnudZiL7zwAueddx7XXXcdCxcuZMCAAVx//fXstNNO3obHzMzMzKwMpXpyh5C243kZ+GlRWXPGSNbVC9JQ37JkKznnFRHRr9xnNqJwIa3Ordy2LUMWLVrEHXfcwcUXX8y9995Lt27dOOaYYxg2bBj9+/evdHhmZmZmZlWtVJJ7DSlBfa+BsvbUt5n1CpPvtohxh4LzaWW28VVJLwD9gCWk/YUfAi6NiGdaFp51ZG+99RZ//vOfmTp1KlOmTOHf//43a6+9Nr/4xS84/vjjWW211SodopmZmZlZTWg0yY2IIc0pawdHNXG/JzAQ2B+YC4wCyl39uEGSVgFOyy4fjoiXy2yqN2lv4Y+BHsBG2WuopDMj4vQWhmodzLvvvsuwYcO49da0CPkXvvAFtttuOw4//HD22GMPunTpUuEIzczMzMxqS4dfeCoiJjWnnqRNgHuAQ/lsr2uLSOoE/BZYE1gADC+jmX8CPwJuAaZHxCJJXbM4zwS2AH4qaVZEnFsiluOA4wDWXXfdMsKw9vLRRx/xve99jxtvvJEuXbowatQoBg8ezIYbbui5tmZmZmZmbUgRzR/ZK+kIYH5E/L6Z9fcHukfENWXGl4uk7wB/An4UEee0UpsXAydkl0MjopzFt0q1vwJpyPKWpJ7oL0bEJ029b8CAAfHUU0+1ZijWCubNm8fYsWO54IILWLJkCSNHjuTII49k4403rnRoZmZmZmY1RdLTETHgc+U5k9ylpNWW125m/enAOhHRLj3GkpYjJYovR8TmrdDeOcAp2eVJEXFBS9ts5Dm7kHqhAf43Iv7Y1Huc5HYs8+fP5/HHH2fo0KFMnz6dww8/nBEjRrDZZptVOjQzMzMzs5rUWJJbTvLZnJWVW1K/bBGxWNIC0sJOLSLpbOoT3JFtleBmHis436ANn2OtaMmSJdxzzz1MnDiRW265haVLl7LWWmvxwAMPsP3221c6PDMzMzOzZVJb97D2ABa28TP+S9JGwMpAk8N9m2hnPDAiu2y1oc9WGyKC22+/nbFjx/LEE0/QrVs3hg8fzsCBA9lll11YffXVKx2imZmZmdkyq82SXElbA72A6W31jKLnrQ1cTdo+qOxxvEVDlH8UEeNbHl2Ttio4b5d/LyvPI488whlnnMGdd95J3759ufLKKznkkEPo1q1bpUMzMzMzMzOaSHIlHQkcWVS8qqT7S70NWAX4KinhvLclAUpqaqGnFYB1SAs31e3HUlbPa1GCO6LUSsc52lSUmPgsaXngF9nlf4D7WvpMa33z5s3jxBNP5PLLL6d3796cddZZnHLKKd4CyMzMzMysg2mqJ7cvn9+Op2sDZY15GRiTJ6AGDCEly82Z2zsXODEi7sr7EEm/pD7BPTkizs/x3jHA6Oxy/Yh4veD2dpJ+RuplfiAi3s7e0wXYDjiLlKADjI2Ij/PGbm3rk08+Ye+992bq1KmMHDmS0aNHs9JKK1U6LDMzMzMza0BTSe4DRdejSYlkqR7OpcBs4AVSUrek7OiSa0hJbmMWA7OA54BbI2J23gdIWpe0jy2k+H8s6ccl3nJOjnm6AnbOXkiaT+qx7Ul9z/NS4P8i4uy8sVvb+uCDD9htt9148cUXmTx5MgceeGClQzIzMzMzsxJKJrkR8SDwYN21pNHA3Ij4eVsHVhDDkHZ4TKei8z5N1O+eo+1ppEWstgb6A71Jw7nnAS8BDwMTI2JajjatjS1atIgrrriCMWPGMHfuXG699VZ22223SodlZmZmZmZNyLvw1PpAS3tmO5xseHHZWx1FxBgaGZYdER9RuufbOph77rmHE088kZdeeoltt92WCRMmsPnmLd522czMzMzM2kGnpqvUi4g36uaUmtWiCRMmsOuuu7JgwQJuvvlmHnroISe4ZmZmZmZVpK33yTWrChHBhRdeyEknncS+++7L5MmTWX755SsdlpmZmZmZ5VRWkitpc+D7wLbAF4FSS81GRLRLMi1pVMFDx7bHM636zZkzh2HDhnH99dez1157OcE1MzMzM6tiuZNPSScA5wGdacE81jYyhvqVmJ3kWpOee+45DjjgAF599VXGjRvHT37yEzp1yjWK38zMzMzMOpBcSa6kgcCF2eWlwF+A24GZwIHAGsAuwGDSNkI/AN5rrWCbGyaltxwyA+Dmm29m8ODB9OrViylTprDddttVOiQzMzMzM2uhvD25PyAlkRdExMkAkgAWRsT9WZ3rJV0E3AWMA77RSrE2x/rt+CyrYjfccAODBw9myy235JZbbqFPn6Z2jTIzMzMzs2qQd1zmIFIv6YVF5Z8ZthwRzwLDgX7AyHKDyytb/fmNiHijvZ5p1eeqq67isMMOY9CgQdx3331OcM3MzMzMakjeJLcPsKAoiVwKrNBA3T8Bi4D9y4zNrFVFBKNGjeLoo49m++2357bbbmOllUqtmWZmZmZmZtUmb5I7j5S4FpoD9JD0meVoI2JRVn+98sMzax0Rwamnnsq4ceM4+uijufPOO+nRo0elwzIzMzMzs1aWN8l9B+guqTA7eDU7bllYUdJaQE863grMtoyJCE4//XTOPvtshg0bxhVXXMFyy3mLaDMzMzOzWpQ3yX0+O25cUPYAKZEdJWkFAEldgYuy+9NaEmAdSZtLmijpJUmzJS0p8VrcGs+02jB69GjOPPNMjj32WC655JK6xdLMzMzMzKwG5U1ybyMltAcVlF0CLAB2Bt6W9Aipx3c/0iJVE1oaZLY375PAUODLQPcsjlIvM8aOHcu4ceMYOnQov/rVr7wHrpmZmZlZjcv7G//twM+Bf9YVRMR00r64c4BVga2B1UgJ7tkRcV1LAizYm7czaW/ePbNbM0l78h4GXA0sBGZksezUkmdabTj//PMZPXo0Q4YMYeLEiU5wzczMzMyWAYqI1mlIWpWUgK4DfALcHRH/aoV2rwMO4bN78y4F3o+ItQrqfY20N+9s4BsRMaelz+7IBgwYEE899VSlw+iwpk2bxpZbbskee+zBTTfdROfOnSsdkpmZmZmZtSJJT0fEgOLyVlt9JyJmAte2VnsFmr03r6ThwGTS3ryj2iAWqwLvvfceu+22G6uuuiqXXnqpE1wzMzMzs2VIrvGbko6QdECO+vtLOiJ/WJ/hvXmt2WbOnMnee+/Nxx9/zF133cWaa65Z6ZDMzMzMzKwd5Z2keDVwQY765wK/yfmMYt6b15plzpw57L777kybNo0bb7yR/v37VzokMzMzMzNrZ+WsxJN35eKWrnTsvXmtSQsWLGDfffflb3/7G7///e/Za6+9Kh2SmZmZmZlVQFsvN9uDtOpxS1Rsb16rHiNHjuT+++/nqquuYp999ql0OGZmZmZmViFtluRK2hroBbzbwqYqsjevVY+JEydy8cUX88Mf/pDDDz+80uGYmZmZmVkFlVxdWdKRwJFFxatKur/U24BVgK+SEs57WxIg9Xvzvl9XEBHTJQ0GrqJ+b15IC1KNb+nevFY9br75Zr773e+y5557Mn78+EqHY2ZmZmZmFdbUFkJ9gR2Kyro2UNaYl4ExeQIqFhGzSUlucfmfJD1IG+zNa9Vh6tSpHHzwwWy55ZbceOONdOnSpdIhmZmZmZlZhTWV5D5QdD0amEtaNbkxS4HZwAvAAxGxpOzomtCGe/NaB/fCCy+w995707dvX2677TZWWmmlSodkZmZmZmYdQMkkNyIeBB6su5Y0GpgbEZ/rWTVrL2+99Ra777473bp146677qJ3796VDsnMzMzMzDqIpnpyi60PtFnPrFlTIoIhQ4Ywe/Zspk6dynrreUtkMzMzMzOrl2t15Yh4IyLebqtgJN0maYtWbK+bpBGSvttabVplXXnlldx///2cc845bLbZZpUOx8zMzMzMOpiytxCStIWk8ZKmSHoxe03JygaU2eyewBOSbpG0u6Sy4pO0nqTTgenALwGPZ60Br732GieffDI77LADxxxzTKXDMTMzMzOzDijvcGUk9QSuJO1JC2nLoDqbANsBJ0u6GTgmImblaP5bpEWt9gb2Aj6U9HvgYeCJiHi9kZhWBAYAA4F9gG2yWwuz9i7KEYN1QEuWLOHwww+nU6dOXH311XTq1GZbPJuZmZmZWRXLleRKWh64H/gaKbl9m7QC8ztZlbVJSe66wL5AX0nbRMSC5rQfEfdJ+jowBPgRsDHwveyFpDnADGAmsADoRdond3Xqe6UFfAr8DhgbEW/k+RqtYxo/fjyPPvoo1113nefhmpmZmZlZoxQRza8s/RQYR0oiTwCuigYakDQEuBRYHvhZRJxZVnDS9sBxwK7Aak1UXwI8TUpuJ0XEx+U8s1oMGDAgnnrqqUqH0S6mTZvGFltswb777ssNN9yApKbfZGZmZmZmNU3S0xHxuamyeZPcF0hDkodFxOVN1D0W+DXw94j4as54G2rvK6ThyGuRem5XAD4CPgReAh6NiLktfU61WFaS3IULFzJw4EDeffddXnzxRW8XZGZmZmZmQONJbt45uRsAi4FJzah7DTCBtO1Qi0XES6Rk1pYhv/zlL3n22We5+eabneCamZmZmVmT8ia5c4HOEbGwqYoRsUDSXLyvrpXppZde4owzzuCggw7iO9/5TqXDMTMzMzOzKpB3idqngVUkrdVURUlrkxaGerKcwGzZtnTpUo4//ni6d+/OxRdfXOlwzMzMzMysSuRNcs/Ljuc2o+45QBS8x6zZJk2axNSpUxk/fjyrr756pcMxMzMzM7MqkSvJjYh7gOHA/pLuk7SjpC519yV1ycruJe2je0JE3Ne6IVutmzFjBiNHjmTbbbdlyJAhlQ7HzMzMzMyqSN59cl/LThcBO2SvxZJmkHptVy9ocx4wUtLIBpqKiOhXTsBW+8aMGcPHH3/MZZddRqdOeQcbmJmZmZnZsizvwlN9GyjrAqzZQPlK2ashzd+3yJYpzzzzDJdddhnDhg1j0003rXQ4ZmZmZmZWZfImuUe1SRRmwJIlSxg2bBi9e/fmjDPOqHQ4ZmZmZmZWhXIluRHRnP1xzcoyceJEnnjiCa699lp69epV6XDMzMzMzKwKecKjdQhvv/02p556KjvttBODBw+udDhmZmZmZlalnORah3DKKaewaNEiJk6ciKRKh2NmZmZmZlUq75zcz1DKRnqRFphqNDOJiDdb8hyrbY899hg33ngjo0aNol8/L7ptZmZmZmblKyvJlbQX8ANga2DFJqpHuc8peubmwPeBbYEv0vjKzZC2KGrxM63tRQQjRoxgjTXWYOTIhnabMjMzMzMza77ciaCks4FTKNFzW/yWvM9o4JknAOcBnVujPes4Jk2axKOPPsoVV1xB9+7dKx2OmZmZmZlVuVxzciXtDowAFmfHr2a3PgS+ROpl/TkwE5gB7A2s35IAJQ0ELiQluJcCe2a3ZgK7AIcBVwMLs2cOBnZqwfNWljRG0jRJcyV9IulJSadI6tqCLwVJfSSdK+llSfMlzZT0sKRjtAxORP34448ZMWIEgwYN4qijvDuVmZmZmZm1XN6e3ONJw4/HRcR5QN0iQUsi4jXgNeBRSVcCU4Arga+1MMYfkHpvL4iIkwueuTAi7s/qXC/pIuAuYBzwjXIeJGk94AGgb1Y0D1geGJC9DpW0c0TMKqPtLbL4VsuK5gIrk/4wsC1wgKR9ImJBObFXozPPPJOZM2cyYcIEOnXyGmhmZmZmZtZyeTOL/8mOlxeVf6YXMiLeBk4AvgD8uLzQ/msQKbG+sIlnPgsMB/oBuSd3SuoM3EpKcN8DvhURK5HmHB8MzAG+DlxXRts9gdtICe4/gC0jYmXSvOITgEXArsD5eduuVm+88QYXXXQRRxxxBF/72tcqHY6ZmZmZmdWIvEnuasC8iPh3QdkSGl586h7gU+DbZcZWpw+wICLeKChbCqzQQN0/kRLG/ct4zhCgf3b+vxFxL0BELI2IG0i92AB7SNo5Z9sjgDWA+cCeEfFU1vbCiLgEGJ3VO07SRmXEXnVGj05f8rhx4yociZmZmZmZ1ZK8Se5sPr/w0ydAd0mfWe04IpaS5u6uXX54QBoyvKiobA7QQ9LyRc9clNVfr4znHJkdp0TEYw3cnwxMz86PyNl2Xf3JETG9gfsXk4YvdwYOzdl21Xn++ee55pprGD58OOuss06lwzEzMzMzsxqSN8l9B+gmqVdB2SvZcVBhRUkbAt1JiW5LvENKonsUlL2aHbcseuZaQE9yrsAsaUXq47+joToREcCd2eWuOdreGFi3ibbnAg/nbbtanXrqqfTs2ZPTTjut0qGYmZmZmVmNyZvkPpUd+xeU3UNKKs+UtAaApNVJ83aj4D3lej47blxQ9kD2zFGSVsie2RW4KLs/LeczNqH+3+KFEvXq7q0hadVmtr1pA+8v1fZXmtlu1YgIFixYwNSpUznppJO44447OPXUU1l11eb+E5qZmZmZmTVP3tWV/wwMBQ4BHsrKLiEt+PR14E1JH5Lm0db1po5vYYy3kbYFOgh4suCZ3wd2Bt6W9DKwEbAqKbGekPMZaxWcv1OiXuG9tUjbGLV22z0kdc96d6tKv379mD17NosXL/7cq44kDjjgAIYPH17BSM3MzMzMrFblTXLvAY4CPq4riIgPJH0b+B1pWO6a2a3/ACMi4s7iRnK6nbT37vsFz5wuaTBwFSmx3Tq7tRQYHxF5V0BeueB8Xol6hfdWbrRW67T9uSRX0nHAcQDrrrtu8e2K+/a3v83ixYtZbrnlPvfq0qULm266Kdtvv717cM3MzMzMrM3kSnIj4lNgUgPlj0nqR0o21yEtRjU1Ima3NMCsjZ83UP4nSQ8CexY88+6I+FdLn9lRRcREYCLAgAEDosLhfM5FF13UdCUzMzMzM7M2lLcnt1ERsQSY2lrtNfOZM4FrW6GpOQXnDW2H1NC9OY3WKt12Y4l/OW2bmZmZmZlZgbwLT9WqdwvOS215VHjv3UZrtazt2dU4H9fMzMzMzKwjaLWe3Cr3d9J83k6k1ZAb3OqH+pWS3896kZujcEXlTbNnlWr7peY0+vTTT8+Q9EYzY2hPvYEZlQ7CapY/X9aW/PmytuTPl7U1f8asLXXUz9d6DRU2muRKeq2VHhwR0a+V2moTETFP0iPAN4HdaWBFaEkCdssu787R9suS3iQtyrU78PsG2l4pe3az246I1ZsbQ3uS9FREDKh0HFab/PmytuTPl7Ulf76srfkzZm2p2j5fpXpy+zbx3qB+m6BS9zrcAkmNmERKNHeUNDAi/lp0/wBgg+z8mpxtXwOcDhwsaVxEvF50//tAd2AJkHdlaDMzMzMzM8uUSnKPaqS8FzAKWAV4DLgfeDu7tzawE7ANMAsYS8F2Qx3cJOCHQH/gD5KOjIj7JHUC/he4PKt3R0TcV/hGSWOA0dnl+g0ksecAxwBrAH+RdEREPC2pK2nf4XFZvYkR8Uorf11mZmZmZmbLjEaT3Ij43FZB2bDaJ0m9s7tHRENDa0dJ2gW4ATgWGNhKsbapiFgsaR9gCqkX+15J80jzdFfIqj0DHFpG259I2gu4C/gK8JSkOVm7XbJqdwMnteiL6BgmVjoAq2n+fFlb8ufL2pI/X9bW/BmztlRVny9FNH80saQzgNOAQyLixibqHghMBn4RET9rUZTtSNLKwAhgf2B90oJUrwC/Ay6OiIUNvGcMpXty6+r1AX4M7EXa2/dT0sJUk4DfRMTS1vxazMzMzMzMljV5k9x/kFaw6p7ti1uqbmdgLvB6RGzSoijNzMzMzMzMmiFvkjsP+DQiVm1m/VlA14hYqcz4zMzMzMzMzJqtU876/wF6StqwqYqSNgJ6AvPKCczMzMzMzMwsr7xJ7iOkrYEuk7R8Y5WyVYMvJS1Q9Uj54VlHJ2llSWMkTZM0V9Inkp6UdEr2OTAri6TVJB0l6VpJL0n6j6QFkt6WdLOk/Sodo9UeSadKirpXpeOx6ieph6QfS3pU0ocF38emZD8/V6l0jFadJH1L0o2S3pD0qaT5kl6TdJ2k7Ssdn3VcklaUtIek0yX9MfsM1f3sG9PMNvpIOlfSy9lnb6akhyUdI6mxbWbbTd7hylsBU0mJ7ivA2aTViN/JqqwN7EhauGkT0qJNgxrYc9ZqgKT1gAeo31N5HtAZqPsDyDPAzhExq92Ds6onaRGfXQH+U9Je0oXTH+4A/l9EeMSItZikjYFnqV9Rn4io+A9qq16SdiQtXNknK1pMWq9klYJqX4+IZ9s3MqtmWQJxGXB8QfGnpM6lbgVl50fEye0Zm1UHSTuQcriG/DwixjTx/i1Iu8aslhXNJf3srPu97W5gn4hY0NJYy5WrJzciHgeOI/2iuTFwBfAq6X+sT7PzK0gJ7hLgu05wa1O2sNitpAT3PeBb2dzrFYGDgTnA14HrKhWjVb3lgCeA7wH9IqJbRHQnrXp+ZVZnD+DXFYrPaki2J/qVpB/Sj1U4HKsBkgYBfyEluPcC2wLLR0Qv0s/KAcAvgE8qFqRVqyHUJ7g3ARtlPyNXBL4M3JLdO8mjnqyEWcB9wHjgEOD95rxJUk/gNlKC+w9gy4hYmdQJcQKwCNgVOL8NYm62XD25/32T9A3gDNIXUJwoLyVl7z+LiKdbHKF1SJKGkv6gAbBNRDxWdP8Q4PrscpeIuK8947PqJ2nHiGjsr4xI+hX1P+TXjYi32icyq0WSfghcQPrD3L/ItoVzT66VQ9KKwDRgA+APwIHeJtBai6QpwA6k71WbRMTiovtdSMnHBsDkiDik3YO0Dk1S5+KdciS9TtpFp2RPrqRxwOnAfOCrETG96P5pwJmkDs+vRMQrrRt98+SdkwtARPwtIvYkZfA7k7L/Q7Lz1SJiTye4Ne/I7DilOMHNTAbqPvRHtE9IVktKJbiZKwvOB7RlLFbbJK1P6lH7CDipwuFYbTiclGDMB4Y5wbVWtmZ2fK44wQWIiEWkqRcA3dsrKKseTW0F24S63+snFye4mYtJw5c7A4e24DktUlaSWyciPomIKRFxQ/aaEhEedlPjsr9QD8ou72ioTqQhAndml7u2R1y2zPm04LxzxaKwWnA5aZjVyRHxYaWDsZpQ90vgLRExo6KRWC16LTtuLmm54ptZT+7Xssun2isoq33Z2hXrZpeN5QBzgYezy4rlAC1Kcm2ZtQn1n50XStSru7eGpGbtrWyWww4F59MqFYRVN0nHkkYh3RsR11Q6Hqt+2e4TdaNLHpS0gaQrsxWVF0h6X9ItkvaoZJxW1S7Ljl8CfifpS3U3siTkRtJIglep8LxIqzmbFpw3Jwf4ShvGUpKTXCvHWgXn7zRa67P31mq0lllO2ZYbp2WXD0fEyxUMx6qUpLVJC27M57OrlJq1RF+gbgu9LwLPA0cDq5N2IegD7APcLumyhhowKyUibiVNrVgI/D/gn5LmSZpHmou7AykR/p+ImF2xQK0W5c0BekiqyJD5spJcSZtLmpjtXTlb0pISr8/NFbCqt3LBeamtWwrvrdxoLbMcslVwf0uak7QAGF7ZiKyK/RroCYyJiNeaqmzWTL0Kzk8jrTR6CNA9W1l5XdK6FQDDskXPzHKJiAuA/YEPsqJu1G8ftDzp966e7R+Z1biqyQFyJ7mSTgCeBIaSlinvTto3t9TLzKy1XAjslZ1/LyKeq2QwVp0kHQZ8m7Q4y3mVjcZqTKei82ERMTlbDIhsJfhDSXvJA5ze0LxKs8ZIWlHSDaRtXN4kzXvsTRotsCvwInAY8ISkzSoWqFkF5UpyJQ0k/YLZGbgU2DO7NRPYhfQ/1NWk4RMzgMHATq0Uq3UccwrOVyxRr/DenEZrmTWTpHNIe7ABnBQRv6lkPFadJH2BtF3QEuDYhlYnNWuBwp93b0XEDcUVstWWz80uewNbtEdgVjPGAwcCrwDbRcQ9EfFRRMyIiHuA7bJ7vYFLKhin1Z6qyQHy9uT+gNQze2FEDI+IutVzF0bE/RFxfUQcDWwFBDAO+FvrhWsdxLsF52uXqFd4791Ga5k1g6SzgVOyy5HZUC2zcvyStAXeROAfkroXvqifT0lBedfGGjMrUjhP7R8l6v294Hy9NorFaoyklYHjsssJETG/uE5WNiG73Db7w55Za8ibA8zOVltud3mT3EGk5PXCovLPDEmOiGdJ8+T6ASPLDc46rL8DdXv+bVqiXt299yNiZtuGZLVM0njqv5f8KCLOqWQ8VvXWz47fJf2Fufh1WkHdurKz2zNAq17Zz7u6RDdKVC383alUPbNCGwF1w9tfLVHvnwXn6zdayyyfwhWVm5MDvNSGsZSUN8ntAyyIiDcKypYCKzRQ90+kxRb2LzM266AiYh7wSHa5e0N1JAnYLbu8uz3istqUDVEekV3+KCLGVzIeM7NmqPu5t0n287AhmxScT2/jeKx2LC04LzUCoE/BuaeMWavIdrN4M7tsLAdYCfhmdlmxHCBvkjuPlLgWmkNaHnr5wsJsgYV5eAhOrZqUHXfM5moXO4C0RxuA9560smQJbt0Q5RFOcK01RMQOEaHGXsDPC+rWlZ9YuYitCl2VHdcBDiq+ma0Sf3J2+Q6e2mXN9w/StmcAxzS0aJmkztQPaZ4FeJs9a011v9cfLKlvA/e/T1qYeAlwXXsFVSxvkvsO0F1Sj4KyuqESWxZWlLQWaelyr65cmyYB00j/ff8gaWdIP7glHQBcntW7IyLuq1CMVsUk/ZL6BPfkiDi3VH0zs44iIh4GbsouL5N0kKQuAJLWIf3i9/Xs/k+zhajMmpTNt70iu/wGcKuk/tnvX52y1ZRvB7bJ6lwQEUsqEat1bJJ6Sepd96I+L1yxsLyBfW7PAd4nLS71F0lbZO11lfRd0ppMABMj4pX2+FoaoojmTwORdC1pr7etIuLJrOwc0l8j7wX2iYhPswU6ricNVf5rRGzd6pFbxWV/vZlC2vgeUs99J+qHrz8D7BwRs9o9OKtqktYF6qZFLAU+bOIt53ierrUWSWOA0ZB6cisbjVWrbMje7aSVbiHt6z2Pz+6jOzYiRrd3bFbdJHUD/shnh4suyI6FIyt/BxzuJNcaIul1mjfidlJEDCl67xbAXaRFHCGN7F0B6JJd303KCxdQIXl7cm8j9dwVDr25hPQ/1s7A25IeIfX47kdaSGFCcSNWGyLidWAzYCxpInqQhrM/TZpDuZUTXCtT8T6TfZp4Ff+V0cysoiLiP8COwLHAQ8B/SN+r3gEmA4Oc4Fo5st7cPUlTw24B3qZ+5ORbwB+AvSJisBNcawsR8TTwVeB80iJnXUjf46aSvuftUckEF/L35PYATiKtlvvrgvL9SPNPCocxLwXGR8RpmJmZmZmZmbWDXEluyYakVUl/VVoH+AS4OyL+1SqNm5mZmZmZmTVDqyW5ZmZmZmZmZpWWd06umZmZmZmZWYflJNfMzMzMzMxqhpNcMzMzMzMzqxlOcs3MzMzMzKxmOMk1MzMzMzOzmuEk18zMzMzMzGqGk1wzM7MqIOl1SSFpSKVjaSuSukp6VdICSeu0QntbZf9mD7VGfGZmVh2c5JqZmVWQpCGSxkjaodKxdADDgQ2AKyLirZY2FhGPA3cB35S0X0vbMzOz6uAk18zMrLKGAKOBHZqo9yrwMvBJG8dTEZJWBU4HFgBntWLTY7Lj/0larhXbNTOzDsrf7M3MzKpAROxc6Rja2HHAKsBNEfF2azUaEY9Leg7YHNgXuKm12jYzs47JPblmZmZWUZIEHJtdXtsGj6hr8/g2aNvMzDoYJ7lmZmYVkM3FDWD7rGh0tkhS4atvQf1GF54qqL+DpNUknZct4DRf0huSJkhavaD+epIukzRd0qeS3pR0rqSVm4i5p6SfSvqrpFnZAlFvSfqdpK1a8M+xC2ku7sfAHSWe/2VJEyW9Imle9vW9JelxSWdK+nIjb70+O+4saYMWxGlmZlXAw5XNzMwqYz7wb2BVoAvwH2BuUZ0lOdtcF/gt8MWsvU5Z2feBnSRtA2wI3A70BmYDnYF1gJOBgZK2j4jPPVfSQOAWoE9BbPOyZx0MHCTppxFRznza3bPjXyNiYUMVJH0LuBVYPitalH2NX8xeA4GF1M/B/a+IeFfSq0A/YDfgsjJiNDOzKuGeXDMzswqIiBsiYg3g0azonIhYo+iVd4XhC4EZwFYR0R3oDhxCSkY3AcYBvweeBzaNiJ7AyqRVjZcAg4CjihvNepTvJCW4NwFbACtERI+sbFz2/jMl7ZszZoDtsuMTJepcSkpw7wb6R0TXiOgFdAP6k5LbN0q8/6/ZcfsSdczMrAY4yTUzM6sdC4BdIuKvABGxKCImA+dm908g9RbvGREvZnU+jYgJ1A/pPbiBdseTFoX6bUQcEBF/i4jF2fs/iIhRwI+yumPyBCypK/D17PK5Rup8AfhSdjkkIl6ou5fF/0JE/DwirirxqGey49Z54jMzs+rjJNfMzKx2XB4RHzVQflfB+XkRsaBEnc0KC7OtffbPLv+vxLOvyY6bS+pTol6xL5CGTAN82EidOcDS7HzNHG0XmtHC95uZWZXwnFwzM7Pa0dhw338XnD/ZRJ1eReVbU/9H8fvTQshNWq/omaWsXnA+s6EKETFf0n3At4A7Jf0K+AvwTGNzeBtQ13YXSatExMfNfJ+ZmVUZ9+SamZnVjjmNlC/OUaf4D+BrFZz3aeJVZ8XmBJtZoeC8oR7mOseQhjOvDvwMeByYI2mqpJFZj3Mp8xt5ppmZ1RgnuWZmZlZK3VDi+RGhZr4eyNF+4fDq4l7k/4qIN4FvkFZivgh4mvR7zCDgbOBfknYq8ZzCJLihId1mZlYjnOSamZlZKe9nx26SvlSyZnkK5+GW7I2NiKURcVdE/DAiBmT1DwXeJCXI12cLWTWkru1PImJRS4M2M7OOy0mumZlZZdUtqNSsya4V8CgQ2XlDKy+3SETMoj6R3iDne+dExPXA0KyoD2k7oYasnx3/njtIMzOrKk5yzczMKmt2dlylkkE0JiI+AG7JLkdK2qhU/WbMjW3IQ9nxfxpps7He2TqF822XNFJnYHZ8MEdcZmZWhZzkmpmZVVbdnq97Slq7opE07hTSPNYewFRJR0vqWXdTUm9J+0v6I/C7Mtp/IDsObOT+NpKel3SSpE0kdcqeK0nbAJdl9d4GphW/WVJnYIvs0kmumVmNc5JrZmZWWZOAT4EvAW9Kel/S69nrixWODYCIeI20fc/rpNWNrwRmSZopaQ5pXu0fgP0o73eLPwALgS9L2rCROv2B84CXgE8lzcje80h2bzYwOCIa6sndGVgJ+AC4r4z4zMysijjJNTMzq6CI+CewI/BnUrK4Gmmf2fXoQPvZR8QzwFeAE4B7gRnAyqTfJf4JXE+as7t/GW1/ANycXR7aQJUngQNJPbZPZ8/uSfrjwLOk1ZU3iYiHG3lEXZtX5dhX18zMqpQioulaZmZmZm1I0nakocSvAhtGK/2CImkl4D3S3r0bZb3SZmZWw9yTa2ZmZhUXEQ8BdwP9gANasekTSD3OVzjBNTNbNrgn18zMzDoESf1Jw4//DmwWEUtLv6PJ9roD04EVgC9FxL9bHKSZmXV4HWauj5mZmS3bImKapKFAX2BN4J0WNtkXuAR4xgmumdmywz25ZmZmZmZmVjM8J9fMzMzMzMxqhpNcMzMzMzMzqxlOcs3MzMzMzKxmOMk1MzMzMzOzmuEk18zMzMzMzGqGk1wzMzMzMzOrGf8f7WV3zepTnBoAAAAASUVORK5CYII=\n", "text/plain": [ "<Figure size 1080x288 with 2 Axes>" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# With astrocyte\n", "save = 0\n", "Ttime = 10*1e3\n", "\n", "ge = 1 \n", "gi = 4\n", "ga = 1\n", "\n", "Jee = 5\n", "Jei = 1\n", "Jie = 10\n", "Jii = 0.5\n", "\n", "Jae = 0.5\n", "Jai = 0.5\n", "Jea = 1 \n", "Jia = 0.5\n", "Jaa = 0.1\n", "\n", "# without astrocytes\n", "# Jea = 0\n", "# Jia = 0\n", "\n", "thetae = 10.5\n", "thetai = 25\n", "thetaa = -3.5\n", "\n", "taue=10\n", "taui=2\n", "tauastro = 20\n", "\n", "B = 1\n", "\n", "tauAdapt=500\n", "sigmaN= 3.5*np.sqrt(2)\n", "\n", "tauN = 1\n", "nrVar = 4\n", "dt = 0.2\n", "nrSteps = round(Ttime/dt)\n", "\n", "def allderivs(y,n):\n", " dy = [0]*nrVar\n", " #RE\n", " aux=Jee*y[0]-Jei*y[1] +Jea*y[3]-y[2]+n[0]\n", " if aux<=thetae:\n", " dy[0]=-y[0]/taue\n", " else:\n", " dy[0]=(-y[0]+ge*(aux-thetae))/taue\n", " #RI\n", " aux=Jie*y[0]-Jii*y[1]+Jia*y[3]+n[1]\n", " if aux<=thetai:\n", " dy[1]=-y[1]/taui\n", " else:\n", " dy[1]=(-y[1]+gi*(aux-thetai))/taui\n", " #A\n", " dy[2]=(-y[2]+B*y[0])/tauAdapt\n", " #rA\n", " aux = Jae*y[0]+Jai*y[1]+Jaa*y[3]+n[2]\n", " if aux<=thetaa:\n", " dy[3] = -y[3]/tauastro\n", " else:\n", " dy[3] = (-y[3]+ga*(aux-thetaa))/tauastro\n", " \n", " return(dy)\n", "\n", "########################## \n", "\n", "re = []\n", "ri = []\n", "a = []\n", "nre = []\n", "nri = []\n", "rA = []\n", "nAux1=np.exp(-dt/tauN)\n", "nAux2=np.sqrt(((pow(sigmaN,2))*tauN*0.5)*(1-pow(np.exp(-dt/tauN),2)))#sqrt(((2*pow(sigmaN,2)/tauN)*tauN*0.5)*(1-pow(exp(-dt/tauN),2)))\n", "rk4Aux1=dt*0.5\n", "rk4Aux2=dt*0.166666666\n", "\n", "Tdt=int(1/dt) \n", "yAct = [0]*nrVar\n", "noiseAct= [0]*3\n", "\n", "for j in range(nrSteps-Tdt):\n", " #RK4\n", " aux2 = [0]*nrVar\n", " aux1 = allderivs(yAct,noiseAct)\n", " for i in range(nrVar): \n", " aux2[i]=yAct[i]+rk4Aux1*aux1[i]\n", "\n", " aux3 = allderivs(aux2,noiseAct)\n", " for i in range(nrVar):\n", " aux2[i]=yAct[i]+rk4Aux1*aux3[i]\n", "\n", " aux4 = allderivs(aux2,noiseAct)\n", " for i in range(nrVar):\n", " aux2[i]=yAct[i]+dt*aux4[i]\n", " aux4[i] += aux3[i]\n", "\n", " aux3 = allderivs(aux2,noiseAct)\n", " for i in range(nrVar):\n", " yAct[i] += rk4Aux2*(aux1[i]+aux3[i]+2*aux4[i])\n", "\n", " #noise update\n", " noiseAct[0]=noiseAct[0]*nAux1+nAux2*np.random.randn()\n", " noiseAct[1]=noiseAct[1]*nAux1+nAux2*np.random.randn()\n", " noiseAct[2]=noiseAct[2]*nAux1+nAux2*np.random.randn()\n", "\n", " #save result\n", " re.append(yAct[0])\n", " ri.append(yAct[1])\n", " a.append(yAct[2])\n", " rA.append(yAct[3])\n", "\n", "\n", "fig, ax = plt.subplots(2,figsize=(15,4),sharex=True, gridspec_kw={'hspace': 0.5})\n", "ax[0].plot(np.arange(int(nrSteps-Tdt))*dt/1e3,ri,color='b',label='$r_I$')\n", "ax[0].plot(np.arange(int(nrSteps-Tdt))*dt/1e3,re,color='r',label='$r_E$')\n", "ax[0].plot(np.arange(int(nrSteps-Tdt))*dt/1e3,rA,color='lawngreen',label='$r_A$')\n", "ax[0].legend(fontsize=20,loc=5)\n", "ax[0].set_ylabel('rate \\n $r_X$ (Hz)',fontsize=25)\n", "ax[0].tick_params(axis='y', labelsize=25)\n", "\n", "plt.xticks(fontsize=25)\n", "ax[1].set_ylabel('adaptation \\n a (a.u.)',fontsize=25)\n", "ax[1].tick_params(axis='y',labelsize=25)\n", "ax[1].plot(np.arange(int(nrSteps-Tdt))*dt/1e3,a,color='k')\n", "plt.xlabel('time (s)',fontsize=25)\n", "\n", "if save == 1:\n", " plt.savefig('./figures/rate.pdf',bbox_inches='tight')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.3" } }, "nbformat": 4, "nbformat_minor": 4 }