{ "cells": [ { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "from pathlib import Path\n", "from brainlit.utils.Neuron_trace import NeuronTrace\n", "from brainlit.algorithms.trace_analysis.fit_spline import GeometricGraph\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from matplotlib.cm import ScalarMappable\n", "from mpl_toolkits.mplot3d import Axes3D\n", "import matplotlib\n", "from scipy.interpolate import splev" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Fitting Splines to Neuron Trace SWC Tutorial\n", "\n", "## 1) Define variables\n", "- `swc` the geometric graph\n", "- `df,_,_,_` read the x, y, and z columns in swc file\n", "- `neuron` define a new class inherited from `GeometricGraph` class\n", "- `soma` define the data on the first run as the location of soma" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "brainlit_path=Path.cwd().parent.parent.parent\n", "swc=Path.joinpath(brainlit_path,'data','data_octree','consensus-swcs','2018-08-01_G-002_consensus.swc')\n", "nt = NeuronTrace(path=str(swc))\n", "df = nt.get_df()\n", "neuron = GeometricGraph(df=df)\n", "soma = np.array([df.x[0], df.y[0], df.z[0]])\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 2) Plot the whole spline tree\n", "- `spline_tree` use the `fit_spline_tree_invariant` to locate neuron branches" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "tags": [] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQMAAADtCAYAAAC76BCNAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABw7ElEQVR4nO29eXxb9ZXG/b2SLNvyvq9JvMVx9n2HEHbKWgoFCi1LS+m0pWXetnQZ2mnfFigMU4Z5u1A6LUM7LVBaaNmhUJIQAnESx3G8JN4db7K829rX3/uHcy+yItmSLdkO6Pl8+BDLV3exdJ97fuc85zmSEIIooogiCtV8n0AUUUSxMBAlgyiiiAKIkkEUUURxGlEyiCKKKIAoGUQRRRSnESWDKKKIAgDNDN4TrUVGEUVkIc3HQaORQRRRRAFEySCKKKI4jSgZRBFFFECUDKKIIorTiJJBFFFEAUTJIIooojiNKBlEEUUUQJQMoogiitOIkkEUUUQBRMkgiiiiOI0oGUQRRRRAlAyiiCKK04iSQRRRRAFEySCKKKI4jSgZRBFFFECUDKKIIorTiJLBAoTH42F0dHS+TyOKjxmiZDDP6OnpQa/XT3rN6XTS0NAwT2cUxccVM7E9iyKMcDgcqFSTOVmS5sX1KoqPOaKRwTxDkiR8R9z5ey2KKCKNKBnMMwLd+FEyiGKuESWDeYYkSXg8njNei5JBFHONKBnMM/zlByKZM7Db7fT19UVs/1GcvYiSwTxjrnMGdrud3t7eiOw7irMbUTKIEIaHh7FYLNNuF00gRrFQECWDCKGvrw+j0TjtdtEEYhQLBVEyiBBUKlXQN7S/yCCKKOYaUTKIEPxVCfwhFNKIIopIIkoGEUKwN/lHLT/Q2tr6kbqejxOiZBAhBHuTf9TIoKurK7rMOUsRJYMIIdhlAsxtsvCjRDxRhBfRRqUIIZRlgsfjob6+nvHxcTIzM8nKyoretFHMOaKRQYiorKwMartgIwNZBKTValm7di2xsbG0tbVhMpk4dOgQbW1tGI3GKDlEEXFEI4MQYbVag9ouGDLo6+vj5MmTpKWlsXTpUgAWL17M4sWLGRsbY/ny5QwODnLixAnMZjPJyclK5KDT6WZ9LVFE4Y0oGUQIKpUKl8vl93cej4eGhgZMJhOrV6/22ysgSRJJSUkkJSVRXFyMEIKxsTEGBwepqanBbreTmppKVlYWmZmZxMbGRvqSFgQ6OjooKiqa79P4SCJKBhFCoMjAYrFQVVVFbm4uK1euZHR0NOjcQmpqKqmpqZSVleHxeBgZGWFwcJCOjg5cLhcZGRlkZmaSkZFBTExMJC5r3hElg8ghSgYzgBBi2vKZvwSiXq/n5MmTrF27lvT0dGDmpUWVSkVGRgYZGRksW7YMl8vF8PAwAwMDNDU1IUkSmZmZZGZmkp6ejlqtVt57tpb+onmTyCJKBiFCpVLh8Xgm3Vz+4H2Tey8Ldu7ciVar9bvdbKDRaMjOziY7OxuYsFMbHBxEr9dTX19PTEwMWVlZxMfHn7U3lcfjOcMiLorwIUoGISIUMvB4PMqyIC8vj5UrV57xVI6U6Eir1ZKfn09+fj4wkfgcHBykt7eXoaEhKisrlWRkUlLSWREtRMkgsoiSARNuxCqVatobHILXD6hUKoxGI5WVlZOWBb6YKwVifHw8ixYtIjk5mdbWVpYtW6YsKYxGI0lJSQo5JCQkRPx8ZgIhRJQMIogoGTChp09NTSU3N3fabeXIYCp4PB46OjoYGxvjvPPOm7Qs8MV8yJElSSIhIYGEhASKiooQQmA0GhkYGKC2thar1TqpUhEXFzen5xcIcxEZDA8Pk5SU9JFNwE6FKBkQ3A0uYzr9gLwsSE5OJicnZ0oikPc332t4SZJITk4mOTmZ0tJSZYjL4OAgVVVVOJ1O0tPTFXKYrxtlLsigtbWV5cuXR8ng4wqVSoXb7Q5620A3r3e1AKC7u3va/S0EMvCFSqUiPT2d9PR0ysvLcbvdSqWipaUFIQQZGRlkZWWRnp6ORjM3X6O5IAO32x3UcvGjiCgZEFpk4G9buVpgNpuVasHIyMiC7FqcybHUajVZWVlkZWUBEzmWoaEhDAYDJ06cQKPRKGXMSF6Lx+OJeKIzSgYfc8xmmRCoWhBsb8JCjAymQ0xMDLm5uUqOxW63Mzg4SFdXFyaTiQ8++EAhj+Tk5LDdwHMRGQRTKfqoIkoGTJCB0+kMelv5JvcnIvLeLpTIwOVyMT4+TmpqqvKFD0bctBAQGxtLQUEB+fn5jI6OsnbtWmVJMT4+TmJi4qRKxUyvaa6WCZE4hsvlwuVyLZhkrD9EyYCJMNhutwe1rZxfqK2txWKxnCEikhGKuYnT6eTAgQMkJiZiNBqJi4vDbrczPj4e1idrpCGTl06nY8mSJSxZskSpVAwODlJfX4/FYiElJUUhh/j4+JD2PxelxUj8vT/44ANeeOEF/uu//ivs+w4XPrJkID/pg8kKh7JMcLvdHD9+nMWLF7Nq1aqAX5xglwm9vb0YjUZ27dpFcnIyAGazmQMHDtDU1ITJZCI5OVkJu0O5eeYa/iIZ70pFSUkJQgilUlFdXY3D4SAtLU2pVExVfTmbRUcmk2nB6jdkfGTJoKenByCoppZgqwl6vZ6+vj7Ky8spLS2ddp9TRQZut5u6ujrsdrtys8hISEggLi6ODRs2oFKpGBsbY2BgQLl5Zlvmi1SkEcyyRpIk0tLSlLZtt9utNFy1tbXhdruVSkVGRsakSsXZTgaJiYnzfRpT4iNLBjPNA/iD7ERksVhYtGhRUAw/VWQgJx3z8/NZuXIl77//vt/3yzeX3K0o3zxyma+5uRlJkpSoIS0tbV5vlpnkONRqtVKJgIm19dDQEAMDA5w8eXLS710u11lLBmazOUoG8wW1Wo3NZgtq26nIwGw2c/ToUfLz81m1ahVNTU2zqhIYDAYaGhqUpKPH4/G7XaCbyrfMJzckdXd3U1tbS1xcnPL7ue45CEfCU6PRkJOTQ05ODvDh9fX09NDf348kSahUKjIzM0lNTT1r8inRZcI8Qq1Wz047INzUtL1GXesRLlj3WQqyS5Vtg+1N8N5OCMHJkycZHR1lx44dihlJINIINgHp25BkNpsZHBxUeg7mMt8QieqH9/XJpcvY2Fja29sZGxtDp9MpS6bZkl8kdQwmk0khuIWKjywZhKIq9CUOj8fDC+8/TF33YSxjSWzb0AYsVfYbbGQgb2e326mqqiI9PZ1t27ZN+sJN9eWbif5A7jmQM/m++YaEhAScTidOpzPskttIl0I9Hg+xsbGKNZwQArPZzMDAAI2NjUqydabWcJEUHEVzBvOImUYGZrOZg9X/pGP8AyQ1nL9rCYUpW5RtQ52UNDQ0xPHjx1m5cqXiNRAMwiFG8pdv6OzspKOjgw8++EAxQJFlxbNdj0eaDHxLi5IkkZiYSGJiomINNz4+zsDAAMePH8dqtU6qVExnDRdJwVE0ZzCPCLXfwO1209vbS1NTEyMpr2IddKJSxZCdVMaB+tfYXHEeKXGFIZUhLRYLDQ0NbN26NeSnVCRuKrVaTVpaGmNjY6xbt27Seryuri4s+YZIRwZT9UFIkkRKSgopKSkBreHkSow/a7hICY4gmjOYV6jV6qDJAGBsbAyAHTt28P/t/SGamBgSswb5+95X6W7MpPAreaTkFU5pdCrD6XRSXV2Nx+Nh586dM/qCRUqm7L1P33yDxWKZ5HEQar5hLpYJofwtg7GG8264ii4TziIcP36cwsLCgEYh3gh2mWAyGak8ug+HZpi1G7ajjdHidquxjCWQkGTEaU7kM9evpihrIzB9zmB8fJyjR4+ydOlSLBbLjJ8089Gz4KsclEPuY8eOYbfbp9U3LDQy8IU/a7ihoSHFGk7+mw8PD0+ShYcDUTIIM2bbXeiL3t5e3mj6Ac6EGqzmWE4d+g23bvormhgnktqD1ZTAV274IoWpm1BJE1/+qXIGnZ2dtLe3s3HjRpKSkmhubg7tAn0wnw1MviG3P32Db75hLsggnPvXarXk5eWRl5cHTMyxaGtr49SpU9TU1BAfH6+Q32xl4dGcQZgRSug/1bZut5uGhgYsFguxCR5M1licdi2mUS2P77mJ+CQ7SWnjuOxaTuj3sDhtu/Jef6VFuVfB7Xazc+fOsPT3L7RuxkD6BjnfEBsbS3JyMk6nM2KkEGkFolqtJiUlhZUrVwIflmmbm5vPsIbT6XQhXWM0ZxBmhEIGgRKIviKi9Z4VHO79NSe632HYqWbUkIBKAm28Dct4Im1D72F3jROrSVb26x0ZmM1mqqqqWLRoEUVFRWG7CRa6mMZfvqGrqwuj0cjevXuVEl92dnbY9A2RblTyTSD6lmlla7i6ujosFgupqakKOUzXjXg2+CScdWQwXfLOe1vfcF6uFqxdu5a0tDQAtOokdi66F2frZlbszua547fjcmoQAjQxTkxmB7/Z/xm+uvtFVJJm0jJBr9fT2Ng4aX8zgRCCtrY22tralKSd2+1eUJHBdNDpdOTm5mI2m1m/fv0Z+Ya0tDSys7NnZZsW6chgqhs2WGs4WTrt23A1Hbl//vOf55VXXiE7O5u6ujr5PenAn4EioAO4QQgxIk3s7L+BywELcLsQ4ujp99wGfP/0bu8XQvw+2Os/68ggFBMS+WZyu93U19djs9nYsWNHwM645LgChFuDwxJPQuoYcYkWnLZYHO5R3B4nKrVGiTjq6+sxGo1T7i8YuFwurFYr4+PjnH/++UpGf2xsTLEzn+1N5ItINyr5yzeMjIwoHgeA8kRNS0sL+ok5n2Tgi0DWcIODg7S2tk6yhktMTJyW2G+//Xbuvvtubr31Vu+Xvwv8UwjxkCRJ3z3983eATzChglsKbAUeB7aeJo8fApsAAVRJkvSSEGIkmGs668jA4XCE9B45jC8oKGD16tVT3ghxmhRUGie6ZA92SzxqtQuTNQ51jAuNakKw4nK50Ov1FBcXs3Xr1lknlY4cOYJGo2HdunUAynzFkZERSktLcblcZzQlZWdnhz3bHQ4EyhX4NiM5nc4z8g3BOCNFmgxmIzoKZA2n1+v57Gc/y8DAAD/+8Y+5+OKL2b59+xnv37VrFx0dHb4vXwPsPv3v3wN7mSCDa4A/iAmGOShJUqokSXmnt31LCDEMIEnSW8BlwDPBXMNZRwbBLhNg4gM5fPhw0GG8EB4cNi1CqJBUblQaFyqVGyQAicHBQerq6khOTqa8vHzmFwKKf+DatWupqak54/feDTmZmZksX74ch8PBwMAAnZ2d1NTUkJCQoJDDQkhOBZs4jImJmZTFl6Mh70SdfGN5i7Ui7YHodrvDFn15W8MdPXqUbdu2UVFRwXvvveeXDAIgRwihP/3vPkBubigAury26z79WqDXg8K8koHL5eLUqVPTegPICDaBKIfxDoeD3bt3B/0BS5KK3NQieoc68XjUqFRuNFo3EoLaxvcxDsHatWvp6uqafmcBIISgsbGR4eFhtm/fHlAi66+aoNVqKSgooKCgACEEJpNpUkJL1gFkZWXNi9X3TKsIgfQN8rRpWVIc6RZmt9sdEVsys9lMWloaN9xww4z3IYQQkiRFNIk075GBXq8PiQymyxl4Lwt0Ol3IN8WYs5kRQy5J6UY8VjVOm4aYWDvjth52bL8es9kcdN7CFw6Hg6NHj5KcnMz27dunvXGmWmd6j2wvKSlRpLf9/f20trYCzPmSIhwlxanyDUNDQxw5ckQRDoWSbwgGkcr4z6KsaJAkKU8IoT+9DOg//XoPsMhru8LTr/Xw4bJCfn1vsAebVzIINeyfLjLwrRZ0dXUF/QWVn8QakkjOMBETa8dljyUuycpQbwZZ6xJQqVQhCZ+8MTY2RnV1NeXl5Uo5LpjzCRbe0lsg4JJiumad2SAS+gLvfMPY2BgrV67EZDLR29sbUr4hGEQqJzEL9eFLwG3AQ6f//6LX63dLkvQsEwnEsdOE8SbwoCRJ8pr4EuB7wR5sXslgJq44/sjDu1qwc+dOJRoIRRUnVwlSYkuwJR/FYY0DSaCS3OBRo1ZpJu0zWAgh6O7upq2tTVEmBoPZio58lxRms1mJGoxGI0DYlxRz1cKclJSk5BusVuskJ+ZA+YZgEMnIYDoy+MxnPsPevXsZHBxEkqRuJqoCDwHPSZL0BeAUIK8zXmOirNjCRGnxDgAhxLAkST8BDp/e7sdyMjEYzPsyIRT4iwymqhbI2wfD9iqVCpvDSIehGUlKRFK5GR9OIWdJN+vXllCQvEXZLpTyZm1tLXa7fVplou+NFM6byrvVNzU1lc7OTgoLC+nv76elpSVsVYr56E2Ij4+f5G8gC4N88w3Tma1C5MjAbDZPu0x45plJCf9Cr39f6Lvt6SrCV/3tRwjxJPBk6Gd5FpKB943oT0TkDflpH6xDstlkRiV5sFviSEgxkpYzxPhwCp6k9/jtB9dwy6anSNQUBkUGNpsNk8lEXl7etCVNb79D39ciAblbLyMjQ6lSDA4OKksKnU6nrMtDWevOd6OSP2GQ3E8h1/69+yl8b/xI+RmcDU1KcBaSgdvtDrgs8Ld9sE9xh8NBfV0T+akraRrrxIMaSeUhXmcjNsGBeSyGE4PPsTX/W9PepLKhiU6no7i4eNobJNDSY64UiN7SYu8lRahVioXWqORdmoUP9Q1yl6JWq1UiouTk5HldJiwEzDsZyGF3MKGpPOzkwIEDFBYWTvvEDcbgRJ6TaLPZ2LJlC73OARrb/oRwxLMoP42y9Cup63sFdYKFDXlfnHKZIISgvb2dnp4etm3bRnV19YznLc5Xo5L3kmKqKoU/N+a5mAA1m/376ht88w12u53e3l7y8/NDzjdMhbOhYxEWABnIScFgJL16vR6z2cw555wTlIhousjAarVSVVVFbm4umZmZqNVqliZfyaK8N+gZ6GNw2MqlSy9kff6HElEhhN99ulwuampq0Gg07NixA7VaPat5iwula9FflUKeqyhHP9nZ2WRlZZ014+Bk+OYb9u3bh8fj4fjx49hstpDyDVPBZDKRkpISxjOPDOadDDQazbRPb+9lQWJiYtBNQVNFBv39/dTX17NmzRoyMjI4duwYbrcbrTqZqzfey+/fu4eBfhWHu3/FhaUPKO/z92WXZcVFRUUsWbJk0vGDvaH9kcFCRKAlRX19PWNjY8TFxZGYmDhvwqeZQu6pKCsrU6I/Wd8QTL5hKphMJgoKghYCzhvmnQym0xqYTCaOHj2qLAv27dsX0r59n8xCCJqamhgcHGT79u2K4sw7/M9KWEHFog3Udxyhvr2WHYuHiY/x767U19fHyZMnWbduHampqZN+F0pkMDY2RmNjI0lJSWRnZ58VXYu+S4quri6Gh4cZGxubdkmxEOHdIu0dEVVUVATMN2RlZZGSkjIleUdzBsGewBSRQU9PD83NzTNuEfaNDHwVgN5fTt9cQHpCPi6nBpfGjlrl3+KrsbGRkZGRgLLiYEN9u91OXV0dGzduxOFw0N/fT29vr9K9mJ2dPe0XbqEgMTFRUZROtaSYzTTm+UCgfENra6syaVpORvrmG8KVM5Ak6TImWpfVwG+FEA/NeqdemHcy8BcZeM8hnKpaEMy+ZTIYGRmhpqaGiooKcnNzz9hWJgOXx0pd/7McOrkHbawTjyMdm2sUrfpDsZAQgsrKSpKTk8+Yg+C7z6nIQE5eytcpC5IyMzNRqVSKm05raytjY2OkpKQoJb9IKglnCt8oyN+SYmBgQBlV5+1xMJs1+XwgkL5BtmiXKzBpaWlhcTmSJEkN/BK4mIkGpMOn25MbwnA5wAIgA9/IwHtZsGbNmjNuNDn0Drb64Ha7aW9vp6uri82bNwf8UNRqNXaXib8du58T7e3ExjrITC3k4tV3kRz7oQx8bGwMk8lERUWF8pQIhKmWCQ6HgyNHjpCRkUFqauoZa1BJktBoNOTn57No0SJlIEp/fz9HjhzB5XIpT6JQZx5E8okcaN++Mw78VSlk74azYUnhDX/6Bjnf8O1vf5v9+/cTExPDpz/9ac4999yZEt8WoEUI0Xb6mM8y0cr80SED76e3vCzwt/6WIZNHMF8WudSXkpLCzp07p0z6DDiraGzcR0efgdRMG9uX3sTmxbcQo/7Qsqurq4u2tjZ0Ot20RACBIwO5T0GOUg4dOjTtcsJ7IEp5efkkT4Da2toZC4XCiVCqCcFUKbzbs8+mJYX3tf3xj3/ksssuY+fOnbz00kusW7dOueYQ4a89eWtYTvg0FgQZOBwOampqcDgc0y4LglUVGo1G2tvbyczMZO3atX63cbotqCQNKklDVd/TmC0OVpWXsKHoOpZmXqxs5/F4Ji1b3nvvvaCuzV9k0NPTQ0tLy6Q+hZmUFr3XsN5Z/draWmw2G+np6UoIHg6D1mAwG4/CYJYUTqcTh8MRkSVFJJO1NpuNT37yk3zuc5+L2DHCgXknA7fbTWtrK2VlZUGp9YLxNOju7qalpYWioqKAX85R6yn+8s8/IqldFJepIWYEjTqeHUs/x6KUncp23lqE6UROvvC+oYUQNDQ0YDKZ2LFjxyQym63OwDerL1tw9ff309jYOGleQCS/9OHat78lxdDQEAaDgcrKSiD8S4pIuiiZzeZwiJgCtS2HDfNKBj09PZw6dYrc3FxKSkqCes90Fuj19fXKE3xgYEDp0PNF12ATVXVNSGoYj23FZkogI0NNQkweNuc4dvc446PjNJ/oY83qNYqkVUYwIbG8THA4HFRVVZGWlsaWLVv85kHCKTryteCSM9/Nzc2MjIygUqno7u4mOzs7rE/ZSImOVCoVqampJCYmsnPnzklj6MO1pIike3GYeh4OA0slSSpmggRuAm6e7U69MS9kIFcLHA4Hy5cvx2QyBf3eQKVIi8VCVVUV+fn5yhN8KgVibmoZhcVGknINSJKHhGQzbY25/Gf9w6QkxDM6biVGZ2TDyrVIusWM2sykxk0IioJtjZYkCZPJRF1dHcuWLQuYZwhEBuGCd+Z7aGiItrY2zGYzhw4dwuPxKDfSbJ+ykVQgevclRKJKEUkyCMffRAjhkiTpbuBNJkqLTwoh6me9Yy/MORnI1QJ5zsDg4CCjo6NBv9+fqtBgMNDQ0MDatWsnjV6bSoGYmVhEWq6FcaMObZyd9ORM1JmFeNxxOK0uUtIEJrvESy91UFn1Y7TksGplFsWFJTjValxuJ1rV1OU9s9mMXq9ny5YtJCcnB9wuUBQQjrDb4bbg9tgnFHao6Dc1QYwNkeAmPzmZ1s5TnBoZoLG3Cg1qspKKycnO9Vsvnw6RJgN/RBWoSiFrACC4JUUkJzCHa/kkhHiNCS+DiGBOycBftSAYObI3vJcJQghOnjzJ6OgoO3bsOKP2PlVkIEkSO5bdxN6GJ/G41azMuJa4uLUMjQxQXFRCYWE+A5Z6Bsf6aOtt5pS+jfcPt1Pf9S4qjRpHYz2fWPV9v/uWz2t8fJyKioopiUA+l3D3JniEi97RBl5//2V6ekdwmzMxGY20NKpJSh3EaswmK9eKx5lKYtooPafiySocZcNqI+rkQXpqUnDa3WRkZChP2elulkiSQbDJSV/lYLBLikhNYF7oKlJvzAkZeC8LfKsFM7U+s9vtVFVVkZ6eHlD4M13X4vrc2xkwNvFeVR1HtS+wyJLMzu3nKUSVl7SBvCRYVSAQwkP3+CGqOv/Myc56mgz7OK/8y+i0k8tETqeTqqoqUlJSKCwsDOoLFgkyMDsG+P8e34NT10FKYgZ5+Tq2rTkHJBduKxSXFaDTZuDBjkqKxeYc5diJ4/z1lQ+oPNLJeVtXU7Y0G5XORF+/nYYTDcRqY5VEpL+R7fMRGUyHYJcUGo0mIpGBw+FYkAIxf4g4GfguC3y/LKFGBhqNRtHxr1ixgpycnIDbBuNnsKngCzS0fx2DwUN8wQskJl92xjYTTSxqFqduJyEmm6buf2F8TMMrdT/mhg3/rWxnNBqpqqpSfA6bmppm1cI8WzilLi7Zfg471l5EgjYDjSqOoaEhenp6yExcOmlbXUw6561fxLKyxbz7QR39vcO89b9V5CxtJzVdxbqKDWwr+7pSoTCZTKSmppKTk6OszRciGXhjqiWFXq/Hbrdz4sSJsFYpzpb2ZZgDMtDr9axZsyagiCiUyEAIwdDQEOPj4+zYsWPaNe10ZUiHw0HT8QGWpV9Ci+4Zhqzj1Br+j/V5dwZ8T5quiNzEIvSeTlo6W+ku+4DC5O20dB2ltf0E2zdcpSwLgu1aDEQGM3VhnoAgL30FL7xcychQHJ+65FI0qqltwNWqGApSNnDprhSeO3Y7qxcJPELQUbcUp/MYl67UKJbm8ngxbwWhrJqMBClEovTnvaRIT0/HYDCQkpIS1irF2TBwVUbEyWDp0qVT/j7YyMDpdHLs2DFcLhdLliwJKrk1lRGJrAJctmwZazPLadj7N8aHkznU8fSUZKCS1KxK+Bd08c/R1HWc5w/+B+tyrqZ66PeoYpwkjfaxPekbSktssF2LVquV9957j5iYGLKzs3G5XLMKWxO1edz9+Wt4/3gaT//1Hc7fuZy8lNQp3yOEwOIY5P8OfRabNZ7etlwy8w1o1C6WLVk2qWHLe7yYvDavrq6mv7+fzs5OZWZkuAavzsUAFd9hst7t2fKSQi7ZBlulOFs6FmEBiI6CcSMaHx/n6NGjSq95sKXIQJGBLCvetGnT6Tl4HrRxblKyRhGeOITwIElTe+1dWPYtJOlR6lqO0657B3WMC7tNQ1XnX4lRxbO58CtBr/utVit6vZ7Nmzej0WgwGAx0dXXhdDqxWCxkZ2eTkZEREjlIkoQuJoOtq7bwl+dOsO9gHTdcsnLK9xgdvfzqz/+DMzkRbbyFzIJ+PE4NqTmj5PnoLHyh1WqVzr2srCyMRiMGg4Hq6mqcTqfiBRDqdciYizmLvvtPSEiguLhYWVLIkVBbWxtCCOVap+oNOZuWCfPeDTId23d2dlJdXc3GjRspLCyc1Vh22cWmr6+PnTt3Kh+SJKnISViDyxHD+Ci0dJ3A7OjHI/wfR6VSYbU4SLFsJiXTQe/gKDvLPk98gg23Cw52PKVsNxUZiNPTl0dGRqioqFCENaWlpSxdupTi4mKysrLo6+tj//79VFZW0t7ejtlsDur6AVSShoS0YfLyElBJ6tPX6/9vrlFpSU3WMdafxpghE7s1Fm2ck/GhFPbVvsmgpXHKY8nLA7lxZ+nSpezYsYOdO3eSmZmpXMfBgwdpbW3FZDIFnSSd76GrciRUUVHBueeey/bt20lJSaGnp4d9+/ZRWVlJW1sbRqNx0jWFskz4y1/+wsqVK5EkySNJ0ibv30mS9D1JklokSWqUJOlSr9cvO/1aizQxnFV+vViSpMrTr/9ZkqRpQ5l5jwwCwe12U1tbi9vtnmQzHkwkIcN7mTCdrPgTyx7kNcvDuJ1JPPLrX1Kxu421Sy7hwmXfOGO/NpuNEydOsHnTBVi7D9DUe5T69koklQuIw+VU8/6pRymQrgt4rjIxCSEoKCg4o9dCkiZmLebk5ChJUpPJNKn/ICMjg5ycnIBP20FTG8+98Vd6TmlJSZzeD0IXk8UXrvoqe5r/i+q2t+nrzCYtc4Tc4h7Mo0k88fZXufXc/5dFaTv9vj9QrkCj0Uy6Dt/w27uPIlDPyWz6HoJBqDqDmJgYv0uKEydOKOPUYmNj6e3tDToyWLVqFS+88AIVFRXver8uSdIKJhSHK4F84G1JkuRhn4Hamh8G/ksI8awkSb8GvsDEtOaAWJBkIM9C8FeBCCUykN83ODhIbW0tq1evPkNWLEOriWOR6kp2nL+D8tI/8c+aTvYfqGPLonES4uNRSTGKS9L4+DirV68mNTWVVeIy2vQn6Rs5RUxiDJIkEZ9o5XjfnyjI/ZTfJ5/dbufIkSOKDLuxsTGo0qJv/4Gs129oaCAuLo7s7GxycnKUJ1H38EneO9DNvd85hxWLNgb1N4tRJ3DxsvtYk38V/2z5D/pHIVZtZdydythgKm80fY8vbn3X73uDTRz6ht9yH0VzczMqlUopX3obusxFZDCbhi5/S4q3336bBx54AIvFgsvl4rOf/SwVFRUB97F8+fJAv7oGeFYIYQfaJUlqYaKlGfy0NUuSdAK4gA/lyr8HfsTZQgbyF0m2EQvkbhRKKVIIgd1u5+TJk2zbtm3KRJZ3FLFx0WcoztzEKy/X8fie68nNzOfG9Y9Tc6yWhIQE8vPzladIcepF3H7uUmyuEV6p+wFmxxBD+jQcNi3NiX+hKPbaSceRE5crVqwgOzsbmJkCUa1WKzcNcEbXYmp6IpVt1Wji3EieFLTq4DPakiSRm7yWWzb8CZvDyB+PfwabZRzzeCIxUmrA982kiuBrZ26325UKhbehS6REQTLCKUeWlxQ33HAD4+PjjI+PU1xcjMVimekuC4CDXj97T1f219acAYwKIVx+tg+IBUEGKpUKl8tFU1MTRqORHTt2BMzWBhsZuFwujh07hsfjYfv27dN+0N5kIEkqMnTLuOIqwR/3x9Nr6OEPf/49m3bks6LkKk6ePOm1rURqXBFQxM0bn+Lluu/QbmwHwOaw4ok5c+iLnLiU4e8GClV05P1kcrvdDAz20/KKlVXlyZgNbtpoU4gjFMRpk7hz0ys0Db7Oo795gU+svXLK7Web8Y+NjWXRokVnGLp0dXXhdrtxOp0zMnSZDpGcmZCVlcW11048FC666CL6+vrO2O6BBx7gmmuuCfvxQ8GCIYPKykqys7PZunXrrGchmEwmqqqqKC4uxmQyBfUh+ytDZuoquHzFd6g+0kFH2v/yeq0OtSaWWFWR3xtVF5PJjet/x4mCv/ObZ1+hPaWWwrhdCLGUxsZGRkdH/fo1+CtBzuSmcrqtqCQ1NvcoxI3S3+/mrltvoiB9KQaDgfr6eoxGI2q1GoPBQEZGRtChcWJMHtbRBBISAmsVwq0v8DZ0iY2NxeFwkJiYGBFDl0hOU/J2zH777bdnspup2pf9vT4EpEqSpDkdHQTV7jzvZDA4OMjY2Bhr1qxh0aJF024/3TLB1624ra0tqPPwJ61tbm5mYEDwycuv47UTx+kZOc5P/9+DqB3NXHRrH3mjOi5c9o0zxDwl6Rey45y36Brup3LwURyHdSQkJAQkutnKke2ucfYdf5WTjXqE20V3p4qkbAMJGU7idSplslNxcTEDAwO0t7czMDDAiRMniI39UGKcmJgY8GaOjUliSYnEnoP7KUstJDFlBcyxHNnXlNQ3oTobQ5dILUPCVFp8CXhakqRHmUggLgUOARJ+2pqFEEKSpD3A9cCzTJ7gHBDzRgZCCFpaWjAYDGRlZQVUKPoi0DLBt2nJe5kR6pfU5XJRXV1NXFyc4qJ87bqHEMLDV3dL/P3V12njb/Q2atFZN7NmRQUJ2iyFFGI1SexcegcvN9yD1lRIXlHelEQ3265Fh9vEa28cB7WTFaVLWLZcw8hYNp+8bDG5qcvOOFZ8fDyrVq0CJlq/+/v7aWhoUDL7ssTY+4bK0JWyaUs61TWVVL/+R85VLYEb9oCXHiPSjUq+T+5gDV2CGdUeyWVCsFHL3/72N772ta8BbAdelSTpmBDiUiFEvSRJzzHhd+gCvirERN17irbm7wDPSpJ0P1AN/G66488LGTidTo4ePUpCQgI7duzg+PHjQUuS/YXUsgV6SkrKGU1Lcn9CsB+0PBCluLiYxYsX+xx74ou/Zvky0mw3YrHbeP6lt6js/w8K01Zy09bHlG1jHcWUjP+I9PT0aSOe2UYGGpWOZF0qefk67rj2TjSqWJyeCUs3jWrq8rJOp6OoqIiioiLFUUi+oWJiYsjJyZlI4GlGqG2vwhMnaNFmsN5wmETLACR82BsSyQ696aoJvoYuNptNqVAYjUYlEZmVleW3cWghzFm89tpr5dzCGScohHgAeMDP637bmk9XGLb4vj4V5pwMRkdHOXbsmNLMA6FVCHwZ3ltW7M88RM4xBPNBu1wuDh06xLp166ac06BWqylKuISilUVUFL/Ny0drqa3rpvfAm3zyc8kYx4yM6+NZsWIFw8PDIV+T/FqwN5dWncDucyt47sU9jH5CT1ZSCVp16KGpSqU6wyHJYDBw4sQJ+uwHMPSkEKezMGhP40TGDjbrJickF1KjUlxc3CQrc1k92NHR4dfQ5eM+gRnmkAyEEHR2dnLq1KkzsumhtjHL6OzspL29/Yz9eSOY6oMQgtbWVmw2GxdddJEyZSkQJElS9lmSdhG371yG9px0/u/PL/PnD/4TXYLgM1t+ictmm5OuRbewk5mjYnhEzSlDI1lJwVnITYf4+HglanC4VhGX/hgDoz1ozEsxJWyhpbV1UjvzQiIDb0iSRFpaGmlpaSxbtgyn08nAwIDixJyQkIDJZMLhcISlj8IbZ5MceU7IwOVycfz4cSRJ8mtZHmobsxCCmpoanE7nJHWiP0zXxiyXILVaLQkJCdMSAUw8QZ1Op/JzStwiHA4HFWUZnHQnYRr38Pv9X0Wn07I19btT7GkCgRqagiESk93AX//5F957v4E1K5ewOKdo2vfMBFqNjktW/pvys9VqndTOnJaWhtVqDelzDAXhFB15qweFEJhMJg4ePEh9fT0Oh0MxdAml2hII0a5FLwghOHjwIIsWLZpUYvFGKJGB1WrFbDazePFiSkpKgjIlDfQFtVgsHDlyRGnL3bt3b9BGp943r9xIVVGxml25L2BxDvF/lbdjNNr49Ysvc47pj2wrv4aK7Cv87i9QaXE6MrC5RvjLO8+w70AzV1+2nou3XEmi9sxpUf6ON1vEx8dPamceGRmhurqao0ePotVqFTWkPxOUmSBSUYckSSQlJRETE8P27dsn5U1OnjypdJEGMnSZDhaLJezRRqQQcTKQJImtW7dOOecgWCGRLCuOj48Pigjkfft76g4MDFBXVzfJNzFYo1PvBiS5lOk9B0EXk8Hntj6FyWjkZfErugw99Pb9jt0l6Whz6qnImnzTSpJ0BhmOj49P+yQ02g0crurkvPPyuGL7TcRqkqbcPlKQfQESExNZu3YtkiRhMBgmRQ1y8m6mo/IiLUcGFBNdb2WnbwSUmpqqXEswbcyR7qkIJ+ZkmTDdF0Cj0WC32wP+Xu7u0+v1bNu2jaqqqqC15L6Rgfe+vKcwy9sG86WT9ylPc/anmNTFZOCO0bIy5zKWZyzGNVTO2/W/ROoYpiGjnvUVmyhI3kJaXBnp6en09vayd+9eMjIycDgcOJ1ONm6cup/A7Xbgsjsoylo1b0TgDZlI4+Lizoga+vv7aWlpmXSzBVPykzEXZOAPvhGQdxszMCkRGSgRfLZg3kVHMHXOQF7Tx8TEsGPHDlQqlRJJBEMG3pGB2+3m2LFjaDQaZV/emMoMxRsejwe9Xk92djbbtm0L+CWVJAm1J4ltRXdDERQUqdnb8BQ9+j76LL9Dl/Q4G/JvZ0fJnWzduhWHw8GhQ4dwuVwIIaiurlZKe/5CzYPVdbjtaayr8D8xaq7hL6rydhNavnz5GSU/2TptuqhhvsjAG/4MXQYGBjh16hQ1NTUkJiYqRCd/XsFWhO69915efvlltFottbW1fwPuEEKMwkT7MhNdh27g60KIN0+/7ncq82kR0rNM9ChUAZ8TQjimO4cFQQaBlgnesmLvmn8onYvytlarlSNHjiidkP4QjNTZarVy4sQJ4uLiWL169ZTb+voZlGZcSOm5F9Ixso9Xjz+IeVzDXw+/w/tSKot3HEJjT2fDopsoKipSElveBiHyOjwtLQ2Bi5HxMdJzjTjdM26ACSuCWWL5lvxGRkYwGAy0tLQoHYs5OTlnRA0LgQx8odVqKSgooKCgQJnE3N/fT3V1NXa7nWeeeQZJkrDZbNMmpi+++GJ++tOfotFokCSpCfge8J25al+GBUQGvmtmX1mx7/aheBqMjY1x4sSJM+Yq+Nt2qshgeHiYmpoaiouLA05q8kagJGBR2nn8y7nb6TMfI/m8Cv7y9m9o6qlFl2Cl5+WN5G76A1LcMFeu/hEFS9IpKS3G7fJMegolJSdRXpzD3v0t7D9UzafOXzHvIWmoST5JkpQn7fLly5WOxZaWFsbHx0lJSVGihoVIBt6QDV2Sk5MpKyvDbrfT2dnJ22+/zY4dO9i1axePPfZYwPdfcskl3j8eZEJKDHPUvgwLhAy8lwlTyYplBEsGsthkfHycnTt3TpvVnaoM2dnZSUdHhxLOj42NTXv8QCVDALVKS0HSFvr6+liSsoL8LAeZ8ctpNyagd7yNsV/HyWPfJbW0kcKkc/nk+h9PKoeNjY3RpRdkZ6k4XlfPyoJ68nIXhbQODzdmm/H37ViUcw2yg1BHRwcFBQWTfA4Wwnn7Q2xsLFdccQVPP/00b775ZlAPDy98Hvjz6X/PSfsyLCAycLlcykzC1NTUgLMQIDgycLvdHD9+HJvNRklJSVDlHX+RgcfjoaGhAavVyo4dO9BoNDidzqByC1PZnsmJzL6+Ps7dcRFa7eUAbCqBtuEfMWrtwtBnpn6gnaNtJ/ntDx7j3IsySFn/dwrTV3Px8m+SlGehvDybru5RYmKlD6W3qYnk5RaSlZV1Rl5ldo7LUyOcN5V31FBRUcH7779PfHw8ra2tStQgr89nOy9yLqTISUlJQbUvP/DAAzDRf/CnsJ/QNFgQZKBWq7HZbLz//vtUVFSQmzt1rXw6MpDzA4WFhaSnp8/IJg0meiiOHDlCenq67E2nbDdTZSFMtjyTG6G8UZJ+AQAi38MGy7mkxZVgv07NW4efpXHUhtm+lw8OniQpd5xEbR79/ano0p2kZadxqOU93mmrRtuZQHHyFjJiyxXLMXF6EvSaNWuC+nuEikgqEAEKCwspKSmZJC+WpzLLxJCamhryOUSqY9FXijxd+/JTTz3FK6+8AnCL+PCLMyfty7BAyKC3txej0cju3buDkm5ORQZDQ0McP36cNWvWkJGRQXd3Nw7HtIlUYDIZ+A5ECbTdVAjUVHXkyBGys7MpLS2d8os7YdQ6kaTUJsEnz/8sbSPZpMQu4fiiPTQO/ZWhpngaj9t56q0fo0sfYbg3E4fLBdIwmVnJXLDuOgb6h6iurmZ0dJT8/HxcLldEnoZzJUf2lRc7HA76+/tpb29nbGyM5ORkJdcQzDSjSEUGoUiR33jjDf7jP/6Dffv2kZWV5Z0RnpP2ZZhnMvB4PNTW1uJ0OklISAj6DxeIDDo6Oujq6ppkcRbMVCUZcjVBbs7ZsGGD3zmJwZKBWq0mNjaWvXv3kpWVRXJyMi0tLVRUVAScyDzl+UkaytInlhMXLl/JeZ47sW8x0t38CCmpsTg8EurYMYQjkZg4Dx0DtQiVlZiYCf/G888/H4vFovgmxsfHK1FDuFRy89GboNVqKSwspLCwUMmnGAwGOjo6FEvznJycgFHDQmhfvvvuu7Hb7Vx88cXU1NQcAw4KIf5lrtqXYR7JQHYrzsvLo6SkhH379gX9Xl8ykJ2UPR4PO3bsmPTBhmqg2tPTg8ViYfv27QGfKlMlBn2327JlCy6Xi9bWVmpra9FqtfT29uJ2u2e93tWo4nBIZtavWUvVnnEuvnA3a9YOUNP9OoPdWdisLv5Q+TkyVFtYWbYZTaybtPhEMjLT8Ag3dqtrUukyKyuL3NzcgAKa+USwQ1S83ZHkqEE2dJGjBnlJIX++C6FjsaWlxfvHdd4/zEX7MswTGQTjVjwVvEuRNpuNI0eOkJ+fT3FxsV/RSzBkIEcEsbGxftfxvvsMJRHX09NDf38/F1xwAbGxscqTq7KyEkmSlKfzTLTvupgM7vj0VRQXJ/DK3pep0B4nVuchIW0cm1mHvicOY9p+jr11nKSEp5BGd5KUKBGvyWbbpqVkZKayYfEqhEfFyNC4UrqUQ+1wJOjChZkQlK8WQPZUPHz4MB6Ph+zsbOLi4iJCfmdTxyLMMRn4yopnGpqq1WrsdrtS95+KVIJZJsiEIjsfByNHDiaBKCfsLBbLpIjF+8klK/JOnjyJ2WwmIyOD3NzckCYPxWqSuXDLpaTlWHmvrRG7NZ7YeAtqtRu7JQ6VykNimolYnQ2n7hVMHomx0TIefvQkbmEjLzeHjevy2LJxNWUr8kjUrjuDsAKJgZxuK1bXIEnaoKpX8wrvqKG8vByn06mMgxsdHcXlcilRQzDdq9PhbOpYhDkkA+9WYX9SYAg+AaVWqxkaGkKv17N169Yp5y5OFxmMjIxw7NgxVq1ahdFoDLpkON12LpeLqqoqkpOT2bRpU8Dr8lbkeTweBgcHJ63pc3NzycnJmfbLqVHFs7noRnpM79LW18D4YBoWYxyp2aN01JexeHkHVlMidlssqZmjxC+q5Vvn/ysdI/sxDFTSOw73PVDJ2o0ZXHzuJpZkr2DZsmUsW7ZMEQPJpcu0tDTSMuLoE3uobHyL4WE7KhHP6tydwO5p/34LBTExMRQUFKDRaBgaGqKgoACDwaD0vsi5hpkum0wm05RTwhca5oQMrFYrhw4d8mslJiPYfgOPx0NHRwcWi4Xzzjtv2qfnVJFBd3c3ra2tbNmyhYSEBMxm84yrBN6wWq0cPnyY4uLioExeZXgPEPGWI8tfTvnpHCgRZjKZSRu6kfPLBjjQ9Uuc2cMID+SW9OC0x6BLMqKNs2McTmR0IB23+zFG+tMAiYS0MbZcNUpn7Wqe+L9XSEn8B+ftXM/5Gz9BYmwuixYtoqAwn35TA9Vtb/DGwRqGxiyYx5NIy3YzPhzHWOHxoK91IUH+3qWkpJCSkqJEDZMUn0lJIUcNJpOJ0tLSCJ99+DAnZKDValm/fr3fzLxyIqdViFORgc1mU5628fHxQYXR/hKIQghOnDihzGiQG2R8TUsCYaqnhBxpTCd9DuYYSUlJJCUlUVZWpoS0bW1tjI+Pk5qaSm5uriIskvMwGzduJDk5mQ3Fn8ZgPM5fq75JTMYwg705OB1p6BItpGSOk5I9hFotkIQKIbmQhIQmYYTiTe/SXrOchDQXr727H5tuH7tKv4lHOPlH7S/o7BllbEiLx6MhLddCTJyDlKwxtHE2xsZSOX68htzcvBkPWJ0P+Ksm+BqgyINkZWLOzMxUooaPwtBVmCMyUKvVUxKBvM104XxNTQ0rV65Eq9UGbYHuT0h09OhRkpOT2bJly6QbO9TEoC96enpoaWmZdukyE8ghrZwIGxkZoa+vj6amJjwej9Ly7P13zklagzZWy7hZR3ZhH3ZbLB5nDKODaeiSzMQnmcnI02MeT0UlqbCY4vG44skv70GldmC259Pea6C99x5s46kMGhLJLrSQmCKwmwVtVaWULHOA24XLoWXV2pXkpuRFtHQZCUwnOvLuO1i6dOkZtmn+uhUhtJzBD37wA1588UVUKhU1NTX/AG4XQvRKE1/Q/wYuByynXz96+rxuA75/ehf3CyF+f/r1jcBTQDwT1YZ7RBBJrgUhOoKp3Y7kvoAtW7ag0+kwmUwhdy3Ch87HZWVlFBScmfAKhQyEENTV1SmluJaWFkZGRiZFGpGCLNVNS0sjJiYGg8HAokWLaGpqwm63TyoRLk49h3bexGZOI1Y3gpRgw+3W4HGrGezJJjFlHE2MG4tRR1LaOB61FeNICogEhFuNpAK3E7Q6E6nZdlTuJHITV6FNiuf8FTns3HA1T+y7BVQQr02fZAzi23Upr8HT09MXVOnS7XaH9Jn5ixr6+/s5evQoLpdLaawyGo1BRwb33nsvP/nJTwCQJOkV4N+BfwE+wYTQaCkTvQePA1slSUoHfghsAgRQdbprceT0Nl8EKpkgg8uA16c7hwVDBnJ/gjc8Hg91dXU4HA6lLwBC71p0u92Ks9H69esDzmgIRaC0a9cuBgcH6e7u5tChQ8TGxlJeXj79G8MEWbAlhFASsmVlZbhcLgYHB5WnVmbyJaxc9DkSUgUtg/vZd/J3xOpsaGKcxGgdWEw6UjKHcHvUOO1xxCdaSM8dwDyeiNMegyrGhcUYh8OmJSdpDdedfx3Dp7Q4HA7WrluD3liFUI+jS9QSq5lcgpTnGpSWluJ0OiedV1JSklK6DEYlGEnMRmfg260oX+djjz3G22+/zcDAANdddx233HLLlNGRT+ScwMQNDhNdi384/WQ/KElSqiRJeUxkat8SQgyfPo+3gMskSdoLJAshDp5+/Q/AJzmbyMD3BpcnFWdnZ58xQj1UIZHNZqOxsfEMZyNfhBIZqNVqUlNTaW5upqKigvT0dAwGAx988AFqtVqpAkRizShXKuQSmfffRqPRkJubS25uLkIIxsfH6evro729H6RCtud+nZcP/ZFFFSPokkeI1dmwmROI01lwu7W4XWoktURSqpGUjDEQEurMcTxuDWj/yYEGB/nxW4jLtvLzd76PFGPEao5HoxHkJ20KeM7e05C8z+vw4cMIIZTkqL+OxEj3PISzN0G+zocffphjx47xgx/8gMrKyqDI5r777uMPf/gDwC3A+adfLuDM7sSCaV7v9vP6tFgwZODdxizPVvCeVOyNYMlAbghyu90hD1+dDrIJqvc5ytoBed5AfX09Vqt1Utg+2y+dzWbj8OHDFBUVBTWcRc6Qy5oGvaGLNQUX8cG7x7DaNOz+5Chj6hO4HFo0WitOazxI4BGxTJjjSGjjXFiNGozDqTh1R2joaCUtbxA3KlQeQVKqGZ0mm9ykdUFdg+95ORwODAaD0pHo7Zmo0Wiora1VZjlEApHsTVi/fj07duwAph+6+sADD/DAAw8gSdKfgLuZWAbMGRYMGcg5g66uLtra2ti8eXPA5EswqkI5ssjNzQ268hAsGci9C94mqN7wnjfgG7YnJyeTm5tLdnZ2yLkFuXlq5cqVM7o54uLiKF6ylOIlS9m9s4E3332V536RwK5PG0jMHsZp1xATbyYuXmAaTSBe58Fu0eJwxxCfaEOjdeFyxZCUOYzVrCUmxo3VOBH5ZC6Kw+mxhjT6XYZWqz3Dx6Cvr4/m5masVitJSUmsWbMmYhFCpMjAarVOikRDGLr6JybW+j8kcNdiD5NFHYXA3tOvF/rZflosKDLo7OwkJiZm2lkI030h5ClL8lO7q6tryu1lBEMy3maqwax1fcP2sbEx+vr6aG1tRaPRkJOTQ25u7rRZ58HBQerq6gI2T4WK/NQV3HLFIrasq+Ht99I5dvQkcfFxLNvUBbjQaIzYrWoklcT4UDJJaRIqlQdtrAPwEKsDEKg0ToRbTVOTxFvqR7l8zff9fj5CeOg311Pb9xcGRwdYuegcilIvIFGbO2l7OTmampqK2WwmMzOThIQExVMiIyODnJycsJYuI9WbIElS0OTV3NzM0qVL5R+vAU6e/vdLwN2nnYy2AmNCCL0kSW8CD0qSJI/+ugT4nhBiWJKkcUmStjGRQLwV+Hkw57AgyMBut9PR0YFOp5tSrRcMent7aWpqmnLKUiBMFRnICTu32822bdtm9OXxlsNWVFRgtVrp6+ujtrb2jCqA99+gu7ubtrY2tm7dGtYSnVadxMrF51Dy6ZUMW5t58a39VO+JITHFTklRJglZLYw7OklKHQMknA4NMVqQVAK3Y2LmqibGjRTjIjXHxcmBN1k6tJ7yzCsnHcfpsvL7Q7cwbh3EYdMiBJxs6yIp5WnOqbierUvumLS9y+VS8kUlJRPToYqKinC73QwNDYW9dBkpP4NQZk9+97vfpbGxUT6PS5ioJMBEhHA50MJEafGO0/seliTpJ8Dh09v9WE4mAl/hw9Li6wSRPIQFQAbyUzw3N1c2g5zRfoQQNDY2Mjo6ys6dO2dU3gtUTZBNTjIzMykrKwtbqBofH6+MSne5XJMUb7L/n9FoZGRkhO3bt0esZBkfk0ZBzBa++MmVDFzYxrvV73C4th5Pr5v8Mi3DwymkZI6RlVCI0dOO26lCUnsQQoXdHEt8ghVNjMAFvN/xKyzOQVRSDIOWRlr1tVhdwzidIEmChAQV6QkFDI+ZGBy2crDpb6wvvEFZXjidTg4dOsSiRYvOUKv6zjQIVLoMNTcTqWVCKHj++ee9f7xK/sfpKsJX/b1HCPEk8KSf148Aq0I9h3klA1kOvGnTJsxmM0NDQzPaj8vlUqY6b9269YybVZYPB9OA5HA4Jn05ZG3C0qVLzzA5CSc0Gs2kbPvw8DC1tbVYrVZSUlLo6uoiNzc37GImb8SoE8hPXc2u5SkkqhLQZpjpNP4Tt9uCeSyJgVMucgoySE5W41T34nRAvM6DhAanQwWSoLfXzdDIUwiPhEBCCEBo0MY5yE4p5TPrn0SSVJgcffz23TswGqHHeIji1PNxOBxUVlZSWloa1N/au3Qpk+lMSpdnm9FLpDAvZCB385nNZkWkY7fbQ57TJ4RQLM6m6gOQn/jTkUFsbCzp6ens37+f+Ph4kpKS0Ov1bNy4MaA2IRJwu920tLSQn5/P0qVLleVETU0NDodjkmV6uL9wXV1ddHSc4uJtnzt9E92D1TlC89Cb1HdUcsrQxonaeIqW65BUdpxOQWJsCks012CJP0a/qglJ4wahQuVJIiUxlfTEfBanbaUi8xrksfaJ2lwWpa+h0VzNiLmLvDgblZWVVFRUzKi5x5dMx8fHMRgM05YuhRA4nc6wk4HVao0ocUcCc04Gsulpeno6mzdvVj6YULQD8vYDAwPU19dPO0JdTgxO1wQlSRIrVqxgxYoVNDU10dHRQVxcHHV1dUqiz1/1IJzwVzrU6XSUlJRQUlKiSGE7Ojqoqak5o0dhNmhtbaW/v5/t27dP2ld8TBprcm9idc6NdI9XcqztnzT3G9HE23C71LgT+0nIHeScjAcxGAwY+g2oJBU5eVP7NGg1sdjM8VjMLg42H2TVqlUz8rfwhXfpsry8XLFFk0uX3oNbWlpaSEpKmvXfzhdnW/syzDEZfDig9EzT01DHsjscDk6ePBmUL0IoykK5iclkMnHBBRcoo9/kpJW3biDcstpgSoe+UljvHgWtVqtULkJJqMnXbLVa2bp165QTohalbGPR+m2cGn2Xvx37HjFxdtxuqO15lfyU1SxbdtUknwZ5RqGvT4PDZaa1rwaVpGHM4GbH2rVTEvps4GuLJv/N5MngpaWliitRuD7PUFyOFgrmjAymq80HO5Zdlii7XK6gs+vBRh1y7iExMXFS1BIbG6t4DsjSZnltmpKSougGZvN08e06DAbeduIwkd/o6+ujurpaMeqYquUZJv6eNTU1aDQaNmzYEPTNsCR1FxuKrqK+7y84bPE4nTG83Xg/qXEl5CSuPMOnYWhoiL6+Purr69HpdBgT3sfQqyU1w8KWVZ+IGBH4QjZT1ev15OXlsXTpUgYGBhQTmpmYy/jD2daxCHNIBomJiVNm+f31JvhCdhbOysoKqT04GDGR7EFQVFQU0HMBUKTGsm7AWyATGxur/C4Upxy5dDgb9yeAhIQESktLlV4A75bntLQ0cnNzyczMVEhLljWnpaWxdOnSkJ+Ku4q+S4wqjiOdz6BVeZBUHl6q/zJf2LIHlfThjaRSqcjKyiIrK0vxafjd+w8Rq3OSlZTFkMGEhpEZ2ZyHCrnBTAjBmjVrkCRJGawaztJlNDKYAtOtn6Z7evsuMaqrq2fUuegPo6OjVFdXK/bqwSLQk7mqqkrx18vLywu4ZhZC0NLSwtDQUNhLh74tz8PDw8rIuri4ODIzM+nt7VVuhJli++J/pbLxLdRxw3g8GqwuJ0d7n2JTwRf8bi9JEna7HRdGYmLVLC/ZQJI6iba2NsbGxkhLSwtqEOtMIISgtrYWSZLO6HeB6UuXmZmZQcvKozmDWWCqp7f8JfZeYsykc9Efent7aW5uVtyOZgPvJ7Ost5fXzPIXKSMjQ7lWuetwy5YtEZ0jKEmSMgkZJmZLHDlyBK1Wy6lTp7DZbAGbhIJBaV4Fx9uqSUqz4LDGcqT9bwHJoK+vj5ONDUhMiJcSY7MpyDnTp6G5uZmYmBil4Wu2n40QguPHj6NWqycNxJkKsyld6vX6kKtj840FQwb+IISgubmZwcHBM+YuhpJw9JdAlJ/K8r7D/RTy1tvL4Wdvby91dXUkJiZiNpvJycmhoqJiTuvRRqOR2tpaNm3aREZGBk6nc1KTUHp6urKcCHbN3GdqQUINQoWkFthdo3SPVVKYsnXSdj09PbS1tbF0TRqHKj1IkorUuA+XZL6Rljzjoba2FpvNNskPIRTylIlAo9GwYsXMBtSGUrocHx/n0Ucf5ZZbbgn5OPOJBUsGsoFqbGws27ZtO+PDD6VCoNFo0Ov16HQ6EhMTcbvd1NTUoFarp8yehwve4afVauXgwYPodDoGBgYYGRlR8gyRrkvLlmzeScqYmBgl0+7xeJTlxIkTJ4IyZHW6rZiMFhwONUIIYuPsuF1qnq38DkmJ8Zxb9gUqsq6lq7Ob7u5utm3bxt9PfBGnXUNsvJvshMCj3nQ63SSF5uDgID09PdTW1ipP5pycnCmt3IUQ1NTUoNVqWb58eViIN1DpsqWlhS9+8YuYzWauvfZavvSlL836WHMJKRT99GmE/IZgsXfvXnbv3q0k86ZazzY3N6PT6fw6FvnC6XTS29tLX18fFosFp9NJfn5+0OFiuOCvdGixWOjr66Ovr0+pAOTm5oZ90nB/fz8NDQ2KW1Sw52swGDAYDEoOJDc3d5Jd+qitg5/+6XvEJRrJyIihOL+I1v4qXA4VkiRw2mPITM6hXPtZVq4qp3n4ZfafeBGPW6I4v5Qb1v1PyNfi7YfQ398fcPaEEIJjx44RFxc3JxGY2Wzmyiuv5KKLLsJmszEwMMCvf/3rmexqXqSLC4oM9u3bx6pVqybNSgyEtrY21Gp1SMkvo9GotDVbrVbGxsZIT08nLy+PzMzMiEYIwXQdyhWAvr4+JWTPy5u9ueiEqnDCNm6mrkJyDsRgMGA0GpUSXGpaEm+f+Dn17UcZHRWM9CVRvtJBYpYek9WKzRyLwx5DrM6BRuXB6VDj8aiIi/fwhXOfIiUuePfoQJB1IAaDQdE05OTk0N3dTUJCAsuWLYs4EVgsFq655hruuOMObr755tnuLkoGb7/9NjExMWzevHnap9epU6dwu91KV9t0kJ+M3jejd/17cHCQpKQk8vLyZuQ1MBXk0uGWLVuCLjn6nltiYqISsocy4UhWFW7evDlsKjvvcxsaGkKn05GWHYcnXk9rXw3vvt+BFN+HR0Bq1ggIFRazDo3GTU62ICu5nG0ln6MgObAz0mzOTba4czqdSknV16w0nLDZbFx77bXceOON3H777eHY5ceXDDweDw0NDXR3d3POOecEVZ/t7u7GZrNRVlY27bbt7e309PSwefPmgE9GOfTU6/UYDIYZq/l89ymXDjdt2jTjm1E23dTr9fT396NSqZRzC5Rl91YVrl+/PmJRj7eNuMFgwCPcxGcYOdm/D/1oO06nCovdSmrWGHZrPBduuIzti/+fiJwLTHyXZPfr8vJyTCYTfX19GAyGoGZPhAqHw8F1113HVVddxV133RWGKwA+rmQgtwenp6czNjbGypUrgyoj6fV6xsfHWbZsWcBtPB4P9fX1OBwO1q1bF1KobbFY0Ov19PX14fF4yMnJIS8vL2jJqnfpcM2aNWG9GWVbtb6+Pr8+CN6qwlWrVs1pXkTO9zidTiSVIDVdx4hUw5GOfxCn1XHbeQ+SGj9zXcNU8Hg8ijekl1GIAnkZZjAYGB0dnbWmwel0cuONN3LBBRdw9913h+MSZHz8yEBOqJWXl5Ofn8/Ro0cpKysLSo7b39/PwMAAK1eu9Pt7p9OpqOt8TUNDhcPhUJJ8FouFrKws8vLyAnYNyuYc6enpM1L2hQKXy6XkGcbGxkhJScFoNJKTkzMna2VvuN3uSdctj4vr6+tjYKiXhMQECnKLQl7qBAOPx6McO5ho0VvTMDAwQExMTNCuUzDxd7/lllvYunUr3/jGN8JxCd74eJGB3KvgvYavqalh0aJFQUmNh4aG6OnpYc2aM0tT081HmA3k3gS9Xs/IyIhSl8/KykKtVodkWBpu2O123n//fWJjY3E4HMTHx5OXl0dOTk7E7chdLheHDx8mNzeX4uLiM34vLyfkkF2uAOTm5s66Qcjj8XD48GEyMzNnPM5M1jQYDIZpNQ1ut5vbb7+dFStW8L3vfW/G5z0FPh5kIISgtbUVg8HApk2bJn1J6+vrycrK8uuI7IvR0VHa29tZv379pNflyczTtTWHA7LMV6/XMzg4iFarxWQysXr1avLy8iJ6bF/I8yzLy8uVY3vffEDE2rCdTieVlZUsWbIkaAK02+3Kufk2CIXqUiT3qwSbTJ4OsqbBYDAwPDysaBqysrLQarXceeedLFmyhB/+MGLmxR99MvAW+6xevfqMD72xsZHk5OSgbiSz2UxlZSUrVqxQnspdXV20t7ezefPmOR/nNTAwwPHjx8nKymJ0dFQxQs3Ly4v4ucgl06nKsfLN19fXF9Y2bLvdTmVlJWVlZTN2gnK73cpyQr755ArAVMsJt9vN4cOHycnJ8RuNhAPemoYHH3yQAwcOsGjRIn7+859HUqfy0SeDI0eOkJGREfCDa2lpIS4ujsLCQr+/94X8VO7v71ccbsNtGhoM/JUOZXcivV6Py+VSEpCBmpZmCn+qwukgL3X6+voYGRmZcRu2HI0sX748qGguGPgKilQq1aTlhPc1HDp0iLy8PIqKisJy7OnO65577sHtdrNlyxZeffVVfv3rX89quO4U+OiTwXRec+3t7UiSFNKHK3sQaDQadDod/f39ylSbUFuJQ0WwpUNZ/6/X6zGZTJOeyrOpMsxEVejvGuREWn9/P3FxcUG1YZvNZg4fPszq1atD6vQMFTabTamcWK1WMjMzyc7OpqWlhYKCgll1XAYLIQT33nsvDoeD//7v/56LpOxHnwymQ1dXFw6HI+gkkJysW7x48aQvhVwW1Ov1AEq4Hs6W0pmWDn1D4pSUFPLy8kK2Levu7qa9vX1WqkJ/kNuwvUuqcp5BvgnkZclUcysjAbfbTV9fH3V1dQBKJ2i4RWLeEELw/e9/n6GhIR5//PG5qs5EyaC3txej0TildkDG2NgYR48eZfXq1VP65slrZb1ej91uV8J1b319qAhX6dD3qRxs9j8SqkJ/kCXIfX19Sht2YmIi7e3tbNq0KSzDXEKBy+WisrJS6QaVB9L09/dPMp0JJ+n/+Mc/pqOjg9/+9rdBE35XVxe33nqrUjW56667uOeee/jRj37E//zP/yh9KQ8++CCXX345AD/96U/53e9+h1qtpqmp6TIhxJsAkiRdxsRIdjXwWyHEQ2G7OB8sKDLo7+9ncHCQFStWTLmdXq+nsbEx5EEpsuhEr9djNBoVvUAoSTQ5GikuLg46txEs5Ox/X1+fX5WhrCq0WCxs2LAh4t2W3nC73bS3t9Pc3IxWq1WMWCP5VPaGTASLFy/2W7HwFmLJpcHZJkgfeugh6urq+P3vfx+SYE2OSjds2IDRaGTjxo38/e9/57nnniMxMZFvfetbk7ZvaGjgM5/5DIcOHaK3t5eSkpI2QB7p3QRczMQA1cPAZ4QQDTO6oGmwoFqYp/MokMuS/f39Z/gbBANv9x85XO/s7OT48eOkpaUpDUuBPnhZJBUuF19fJCUlkZSUxNKlS7HZbGdMWzIajcTHx7Nx48Y59+QfGhqiu7ub3bt3ExcXpzyVW1pawiLdngpy6bKoqCggAfvOt/TnUxmK0vDRRx+lurqaP/3pTyE3icm+BzDxmS5fvpyensDjDl988UVuuukmYmNj5eR6C7Dl9K9bhBBtAKdHrF0DfPTJYCpTVFliK0mSX3+DUKFWq5W2VyEEQ0ND6PV6Ghoa/DYshXvW4XSIi4tTvtw228RMAY/Hg8Vioba2VulmnIvoQK/X09zcPGm+pPeYOLkN29uINVxt2DIRFBcXBy0g8zUiGR0dVYjLe/ZloKTrL37xC/bv389zzz0362VYR0cH1dXVbN26lQMHDvCLX/yCP/zhD2zatImf/exnpKWl0dPTw7Zt27zf5j1G3Xfs+mTHmDBiQZFBoMjA4XAo6raSkpKwPxUlSSIzM5PMzExlOKper6elpYXY2Fji4uIYHR1l27ZtEa1O+INsAltUVMSSJUsUia9er6eurk7RZczWnTkQuru76ejomNKj0Xeug2z04d3qPJMWcafTycGDB4OesOQPshtyWloay5cvP2MgjW9fx29+8xvefPNN/vrXv856+WMymbjuuut47LHHSE5O5stf/jI/+MEPkCSJH/zgB3zzm9/kySfPmI42b1hQZKDRaDCbzRiNRkUlZzKZOHLkiN9ZC5GA73DUhoYG+vr60Gq1VFVVKZWJuZiW409VqFKpFNckb+KajTtzIHR0dNDb28u2bduCJhrvpZivRXoobdjyqLWysrKwqjl951v29/dz6tQpvve979HY2MjY2BhvvfXWrCs0TqeT6667jltuuYVPfepTAJMmRX3xi1/kyisnBtQWFBT4Tgr3HqPubxx7RLCgEohCCHp6eujp6cFqtZKUlMTIyAibNm2a0xIW+C8dzpWQCD7MT4RSxzebzUqnJTCp0zJUeOsnwjF6zFtMZDAYpsz+OxwODh48SHl5+Zw8AAD++Mc/8rvf/Y7Nmzezd+9e/vSnP/ntfAwGQghuu+020tPTeeyxx5TX5VkNAP/1X/9FZWUlzz77LPX19dx8883eCcR2YCkTJcYm4EImSOAwcLMQon5WFxsAC4oMvNHa2kpHRweJiYmKfHaqTsFwIpjSoVx20+v1QXUyhoKZqAp94Ss/ltfx052fPM3abDZH1AfBtw1bPj+dTkdlZSXLli2b0czFmeAvf/kLv/rVr3jllVfCUpZ87733OPfccydJ7h988EGeeeYZjh07pgjrnnjiCYUcHnjgAZ588kk0Gg1NTU2XCyFeB5Ak6XLgMSZKi08KIR6Y9QkGwIIjAyEE9fX12Gw21q9fr1iiDwwM0Nvby9jYGBkZGeTn589awecPNpuNQ4cOUVJSEnTpUM5e6/X6WVuphUNVONPzkweMuN1u1q5dO2cVCzlc7+npwWAwKJJ1ueckknjxxRf52c9+xquvvhrxOZohICo6cjqdHD16lJSUlIC9+N4JtOHhYVJTU8nPzycrK2vWxCAPaplN6VBeJ+v1eoaGhkJK8EVKVeh7ft6dlt5NQQ0NDWg0mjk3ioWJSObgwYNUVFQobtaDg4PodDolzxDuv8lrr73GAw88wKuvvhrSMjSQqGh4eJgbb7yRjo4OioqKeO6550hLS1P6Gl577TV0Oh1PPfUUGzZsAOD3v/89999/PwDf//73ue222+DjTgYWi4XDhw9TWloa9BNZbiHu7e1lcHBwVpn1SJQO5bKW3EwVHx9Pfn6+3wTaXKkKfc9Ptnpra2tDq9VSUlIyJ52W3pBLp3IHqvf5eduWwezyIN74xz/+wb//+7/z2muvhdxsFEhU9NRTT5Gens53v/tdHnroIUZGRnj44Yd57bXX+PnPf85rr71GZWUl99xzD5WVlQwPD7Np0yaOHDmCJEls3LhRNuT5+JKB7EGwdu3aGXeB+d54Op1OufGmKxHNxLB0JpB9DPv6+pRaeG5uLu3t7fOiKoQP/QAyMzPJy8tTEqRut1tJ8EUiQSrDZrNx8ODBKSdPe2/r3bQ0U5Xhnj17+M53vsNrr70WFvHYNddcw913383dd9/N3r17ycvLQ6/Xs3v3bhobG/nSl77E7t27+cxnPgPAsmXL2Lt3r/LfE088AeC93byQwZyWFg8ePMihQ4f41Kc+pTz9T506RUdHB1u3bp3VGtm3nizfeLLzj3zjeYea8sSm4eFhduzYEfEnsqwwLC8vx2q10tvby759+wAoKSnBYrHM6bBOl8vFoUOHJnX/yXoBOUF68uRJJUEa7jH0VquVysrKoJdlcXFxyiwN2YCks7OTmpoaxQV5uoav9957j3vvvZdXX301LETgLSoyGAxKQjA3N1eJZnp6eiZJqAsLC5Wqmb/X5wtzSgalpaUcPXqU2267DZvNRmZmJhkZGTzxxBNh1bdLkkRycjLJycksW7YMk8mEXq/n0KFDqFQqhRiam5sBIj7r0B+0Wi2Dg4OUlJSwePFiDAYD9fX1WK1WJRQO9yAVbzgcDg4dOhRQ4us7Hs5b3puamqp0Ws40wScTwUxboL2VhN4NX42NjQHbsA8ePMg999zDyy+/HJZKha+oyBuSJM153mW2mFMyyMrK4itf+Qqf//znufnmmzEajdjtdnbt2sU111zDddddN+Pa7lRITExk6dKlLF26FKvVSnd3N3v37kWj0VBcXIzVap3Tibnyjbho0SLliez9xOvv71dmH8rhe0ZGRti+XHLFJNg6vu8Yeu+JzjqdTum0DLZXxGKxcOjQIdasWRMWcxDvGY0rVqxQ8gzyNOyGhgYSExN55JFHePnll2esZvRGIFGRrCXQ6/WK4YuvqKi7u1sRZu3du3fS67t37571uc0YQohQ/5s1TCaT+POf/6z8PDQ0JP73f/9XXHHFFWLz5s3ixz/+saivrw/Hoc6A1WoV+/btE11dXcJms4mOjg7x/vvvi3379omTJ0+K8fHxiBxXhsViEXv37hW9vb3Tbut2u4XBYBDHjh0T77zzjjh69KjQ6/XC5XLN6vh79uwR/f39M96HDI/HI8bGxkRjY6N49913xXvvvSdaWlqE2WwO+B6z2Sz27NkjhoaGZn38YGC328UDDzwgioqKxPLly8W//uu/CovFMqt9ejwe8bnPfU7cc889k17/1re+JX76058KIYT46U9/Ku69914hhBCvvPKKuOyyy4TH4xEffPCB2Lx5sxBi4ntfVFQkhoeHxfDwsCgqKpL/LjO5L2f937yQwVQYHR0Vf/zjH8W1114rNmzYIH7wgx+IY8eOCY/HM+t9j42NiT179oiBgYEzfudwOERnZ6c4ePCg2LNnj2hoaBAjIyNhOa6M8fFx8c4774jBwcGQ3+vxeMTQ0JCoq6sTe/bsEYcOHRJdXV3C4XAEvQ+j0RjRG9FisYi2trZJ5Do6Oqr8DU0mk9izZ48YHh6OyPH9oba2Vqxdu1Y0NzcLm80m/vGPf8z6M92/f78AxOrVq8XatWvF2rVrxauvvioGBwfFBRdcIMrKysSFF16o/J09Ho/4yle+IkpKSsSqVavE4cOHlX397ne/E6WlpaK0tFQ8+eST8svzQgYLopoQCCaTiddee43nn3+epqYmLr30Uq6//voZtfDKpcONGzdOKy6RQ3XZbCUrK4v8/PxZqQvDoSqUIcSZ05+ms3kbHx+nqqqKDRs2kJKSMqvjBwPZ6q2vrw+j0UhKSgrDw8Ns3Lgx4q7VMk6ePMlNN93Es88+S0VFxZwcM0z4+JYWg4HFYuHNN9/k+eefp7a2losuuojrr78+qJHqsykdes9JGB0dJT09nfz8/JDahyOhKvSGb0+Cr82bTESbNm2aF5Xd+Pg4lZWVpKSkYDabZ2zAGgpaWlq4/vrr+dOf/hRw0M4CxvxkHmcQTsw7bDabeOWVV8Ttt98uVq9eLb72ta+Jffv2nbGW9ng8orGxUXzwwQfC6XTO+rhut1v09/eHtIbv6uoS7777rrDZbLM+fjCwWq2ivb1dvP/++2Lv3r2iqqpKvP3228JkMs3J8X0hL41GR0eFEBOfyfDwsKivrxd79uwRBw8eFB0dHWH9+7S3t4t169aJY8eOzej9d9xxh8jKyhIrV65UXvvhD38o8vPzJy0LZDz44IOitLRUlJeXizfeeEN5/fXXXxfl5eWitLRUySUEiWjOYCZwOBzizTffFF/84hfFqlWrxL/8y7+It956S5jNZvH444+LY8eOCbfbHfbjymv42tpa8c4774jDhw+Lnp6eSaTT0tIi3n///bAQ0UzQ3d0t/vGPf4iDBw+Kd955R9TW1orBwcGw5kGmgkwEY2NjAbcxGo2iqalJ7N+/X+zfv180NzfPirg6OzvF+vXrxZEjR2a8j3379omqqqozyOCRRx45Y9v6+nqxZs0aYbPZRFtbmygpKREul0u4XC5RUlIiWltbhd1uF2vWrAklKT4vZLCg/AxmgpiYGC655BIuueQSXC4X+/fv55lnnuGWW25hw4YNLF68mOXLl4d9tp93OUuID30FmpqaiI+PRwihzHGYaw0DTJjLtrW1sWvXLmJjY5XljmzzNptmqmAg5yimy5F4l31lhaFs9ZadnR2S3qK3t5dPfvKTPP7442zcuHHG575r1y46OjqC2tbXsqysrIxDhw4BUFZWpkx5uummm3jxxRen9fecT8z9tzSC0Gg07N69m66uLh599FHuu+8+3nrrLTZt2sRtt93GSy+9hM1mC/txZUOU5cuXs2vXLiRJwuFwKMKaU6dOYbfbw37cQOjs7KS9vZ1t27YpiktZK7B+/XrOO+88CgoKMBgM7Nu3j6qqKnp7e6f0nwwFMhGE6qAsKwy3bdvGzp07SUlJobW1lb1793L8+HFlWI4/GAwGrr76av77v/+brVsj4wz2i1/8gjVr1vD5z3+ekZER4OxRFwaDsz4y8IUkSfzlL39RZL3nnHOOMpjz+eef59///d9Zvnw5119/PZ/4xCfCmtCTdf5paWmKD4Kc3PNWP0ayEaitrQ2DwcDWrVsDJudUKpVfmzc5qpEFRjPpEpQt7GebrNRoNOTn55Ofnz+pE7S+vp7k5ORJzswDAwNcddVVPPLII5xzzjkzPuZUWOiWZeHAR44MgDP0/SqViq1bt7J161Yefvhhjh07xl//+lceeOABSkpKuP7667niiitm9eX1pyoESEhIoKysjLKyMqxWK3q9nqNHj+LxeMjNzSU/Pz9s6sempiZGR0fZsmVL0DJhb5u35cuXnyHdls8xGPIaHR2lurqazZs3h7XHQqVSkZWVRVZWlkJefX19HDhwgEcffRSTycR9993H+eefH7Zj+iJYyzJZXQgEfH3BYgaJho8MPB6PqK2tFT/84Q/Fpk2bxNVXXy3+8Ic/iJGRkZD2E4qqUIasfvzggw/E3r17xcmTJ8XY2NiMknsej0fU19eLqqqqsCZLZRHRgQMHxL59+0RjY2NAhebIyIh455135rRqMTIyIs477zxx5513it27d4vvfOc7Ydt3e3v7pASi92f76KOPihtvvFEIIURdXd2kBGJxcbFwuVzC6XSK4uJi0dbWpiQQ6+rqgj18tJow3zhx4oS4//77xdatW8UnPvEJ8bvf/W5ateD4+LjYs2fPjFSFMhwOh+jq6hKVlZViz549or6+Pmj1o8fjETU1NWFTaQaC3W4XnZ2dk85xeHhYKRXu2bNnTolgbGxMnHvuueLvf/+78tpsZNreuOmmm0Rubq7QaDSioKBA/Pa3vxWf/exnxapVq8Tq1avFVVddNYkc7r//flFSUiLKy8vFa6+9prz+6quviqVLl4qSkhJx//33h3IK80IGZ43oaK7R2trK888/z4svvohOp+O6667j2muvnRQuhlNVKENWP+r1esbHx6ec+uTxeDh27BhxcXEsX758Tm3KZBu64eFhXC4Xq1atoqCgYE4qJyaTiauuuoqvfOUrfPrTn4748eYBUQXiQkVnZycvvPACL7zwAmq1mk996lPExcXR29vLvffeGzHbdHnqU29vr6J+lMuBMDHiPjU1lfLy8mn2FBnIpjRlZWUMDw/PapBssLBYLFxzzTXccccd3HzzzSG///Of/zyvvPIK2dnZygDXMNqVhQtRBeLZgJ6eHnHbbbeJ7Oxscd5554mf/exnoqOjI+LHldWPNTU14p///Kd44403RHV1ddhC41AxODgo9uzZM6kD0LeZqrKyUnR2dobUTDUVrFaruPTSS8X//u//zngf/gRF995776Ruw29/+9tCiIkw37vbcMuWLUKIiW7D4uJiMTQ0JIaHh0VxcXG4m6/OvpzBc889J1asWCEkSZrUidXe3i7i4uIU6eaXvvQl5XdHjhwRq1atEqWlpeJrX/uass4dGhoSF110kSgrKxMXXXTRnHa2hYK+vj5x7bXXCqPRKAwGg3jiiSfEJZdcIrZv3y4eeugh0dzcHNHjOxwOsX//ftHQ0CDq6uoU9WN3d/ecKR0HBgbOIAJfyO3NJ0+eFPv27RMHDhwQbW1twmq1zuiYNptNXHnlleKJJ56Y6Wkr8E0OlpeXKzmA3t5eUV5eLoQQ4q677hJPP/30Gds9/fTT4q677lJe990uDDj7yKChoUGcPHlSnHfeeWeQgfcf2xubN28WH3zwgfB4POKyyy5TEi6B2PlsweDgoHjyyScVT4af/OQnoqGhIazHsNvt4t133xU9PT3Kax6PR4yOjoqGhgZF63/q1Clht9vDemwZMhGEelObzWbR0tIi3nvvPfHuu++KpqYmYTQag3qvw+EQn/zkJ8UvfvGLmZzyGfD9fqakpCj/9ng8ys9XXHGF2L9/v/K7Cy64QBw+fFg88sgj4ic/+Yny+o9//GO/UuVZYF7IYFaLuuXLl4e0vZwUk4dM3nrrrfz973/nE5/4BC+++KLi+nLbbbexe/duHn744dmc3pwiIyODO+64gzvuuIOxsTFeeeUV/u3f/o2uri6uuOIKrrvuOlavXj3jJJ/sIFxRUTEpiSlJEikpKaSkpEzyfjx48CAxMTGKyCkcNuMDAwM0NDTMaOakTqejtLSU0tJSZcBLXV3dtLJjl8vFLbfcwrnnnstXv/rVWV/DdDgb7crChYiJjtrb21m/fj3Jycncf//9nHvuufT09Ezy2/OWaAYykzwbkZKSwi233MItt9yieDLcf//9tLS0KJ4MGzZsCPpLJ9uEBWMc6m26arFY6O3t5fDhw0iSFJKAyBf9/f2cOHFiksR5poiNjVVs3nwHtXrbvHk8Hm6//XY2bNjAN77xjVkdcyqc9XZl4cJ0oQPwNlAn/7dy5UqxcuXKSfVd32WCzWZT6u5HjhwRhYWFYmxsTBw+fFhceOGFynbvvvuuuOKKK4QQk0M1IYRITU2dWYC1gGE2m8ULL7wgbr75ZrF27VrxzW9+U3zwwQdTCoWMRqN45513Zp1DsVqtioAo1DC9r69P7N27N+Jt2C6XS/T19Ynq6mqxe/dusWHDBnHzzTfPOM8QCL7LhDDalYULZ1/OQIYvGQT6fW9vr1i2bJnyunciJlAS56MKm80mXn75ZcWT4etf/7p49913J1UHRkdHp20Bngnsdvsk9eOJEycCqh/7+vrEvn375syPQYiJdfudd94pbr/9dvGNb3xDbN68OWyE4E9QFEa7snDho0MG/f39ype6tbVV5OfnK39g3wSibBIRiJ0/DrDb7eKNN94QX/ziF8Xq1avFl7/8ZfHYY4+JW2+9Nein90whqx8PHTok9uzZI+rq6hRloV6vF/v27YtYMtIfPB6PuPvuu8U3v/nNOfNdWIA4+8jghRdeEAUFBUKr1Yrs7GxxySWXCCGE+Otf/ypWrFgh1q5dK9avXy9eeukl5T2HDx8WK1euFCUlJeKrX/2q8oEHYuePG5xOp/jP//xPkZOTIzZv3izuvPNO8frrr8/JDel0OkVvb684cuSIePPNN8Xrr78u9Hr9nN2UHo9HfOMb35hUcv6Y4uwjg0gjkI5BiNCtptra2sSWLVtEaWmpuOGGG+b0aRcK3G63uO2220RfX59wuVxi//794p577hGrV68Wt912m3jppZfCvob2RU9Pj9i3b5/o6ekR1dXV4p///Keorq4WBoMhIq5RQkwQwfe+9z3xpS99KSxEsGTJErFq1Sqxdu1asXHjRiFEYC2Lx+MRX/va10RpaalYvXq1qKqqmvXxZ4koGfgikI5hJlZTn/70p8UzzzwjhBDiS1/6kvjVr341l5cya7jdbnHw4EHxrW99S6xdu1Z85jOfEc8///yUMwpmgu7ubrF///5JqkGPxyMGBgZETU2NeOedd0RVVZXo7e0Nq/rxRz/6kbjjjjvCRjZLliw5wxI/VKXhPCJKBoHgSwYPPvigePDBB5WfL7nkEvH++++L999/X1mqeG/n8XhERkaGotDz3e5sg8fjEVVVVeJ73/ueWL9+vbj++uvFs88+O+sBMF1dXWcQgb9j+85v6O7unpXk+MEHHxSf/exnw0ou/sggVKXhPGJeyOCsNDfp6elRhEswWa/gazVVWVnJ0NAQqampSuPM2WBBNRUkSWLDhg1s2LCBBx54gLq6Op5//nkeeeQRCgsLuf7667nqqqtCmo/Q1dVFZ2cnW7dunXLupa/3ozy/oaWlhbi4OGV+Q7Cekz/72c+oqanhj3/844znNgY6z0suuQRJkvjSl77EXXfdFfJgVHnbjwvmnQwuuugixe/fGw888ADXXHPNPJzR2QVJkli9ejWrV6/mRz/6ESdPnuT555/n0ksvJSMjg+uuu45rrrlmyuGmnZ2ddHV1TWmVFujYsvqxoqICk8lEb28vBw8eVEbO5+XlBVQr/vznP+fAgQP8+c9/DnuH43vvvUdBQQH9/f1cfPHFZwxR+TgrDQNh3sng7bffDvk9oVpNZWRkMDo6isvlQqPRnB0WVDNERUUF9913H/fdd5/iyXD11VeTkJCgeDLICjuAU6dO0dPTEzIR+ENiYiLl5eWK+lGv13PkyBEAhRjkdu8nnniCf/zjHzz//PNhncAtQ/58s7Ozufbaazl06FDISsOPHWawtphz+OYMZmI1df31109KIP7yl7+cj0uZN3R0dIhHH31U7Nq1S5x//vni5z//ufjhD38ofv7zn0e829Fb/XjzzTeLG264QZx77rkRq4qYTCYlf2IymcT27dvF66+/HrLScB4RTSD6IpCOQYjQraZaW1vF5s2bRWlpqbj++uvnVFG30NDT0yNuuukmUVhYKHbt2iUeffRRcerUqTk59uOPPy62bt0qrrjiCrFp06aIJOpaW1vFmjVrxJo1a8SKFSuU78FMlIbzhCgZzBfmeHTWvOPUqVPi2muvFVardZInw44dO8TDDz8sWlpaInLcP//5z2L37t2KV+L4+HjEdAtnOaJkMF+Y49FZCwL+hD2yJ8Pll18utmzZIu6//35x4sSJsBzvb3/7mzj33HNnXf78mGBeyOAjNVEp3Ag0OuvQoUPK6CytVquMzjqb4C+TLnsyvPrqq7z55pssWbKE7373u2zatIkf/vCH1NbWTjxBQsSrr77KI488wksvvRTRKdBvvPEGy5Yto6ysjIceeihix/moIkoGp/FRH50VKlJTU/nsZz/L3//+d/bs2cOKFSv4yU9+wsaNG/m3f/s3jh49GhQx/OMf/+AnP/kJL7/8MqmpqRE7X7fbzVe/+lVef/11GhoaeOaZZ2hoaIjY8T6K+NiQwUUXXcSqVavO+O/FF1/ky1/+Mq2trRw7doy8vDy++c1vzvfpLigkJSVx44038txzz/Hee++xadMmfvazn7F+/XruvfdeKisr/RLDO++8w/e//31eeeUV0tPTI3qOH4Vobb4x7zqDuUKweoaP7OisMEGn0/GpT32KT33qU9hsNt5++20ef/xx7rrrLs4//3yuv/56duzYwYEDB/j2t7/Nq6++Oq07UzjgL1qrrKyM+HE/SvjYkMFUkIUoAH/7299YtWoVAFdffTU333wz3/jGN+jt7aW5uZktW7YghKC5uZn29nYKCgp49tlnefrpp+fzEuYFcXFxXHnllVx55ZU4HA727NnD73//e77whS9gs9n44IMPJvk1RrGwESUD4Nvf/jbHjh1DkiSKiop44oknAFi5ciU33HADK1asQKPR8Mtf/lLRz//iF7/g0ksvxe128/nPf56VK1fO5yXMO7RaLZdeeimXXnopDoeD1tZW8vPz5+z4URXh7BGdqDRPeOONN7jnnntwu93ceeedfPe7353vUzqr4XK5KC8v55///CcFBQVs3ryZp59++mwl6XlpmohGBvMAOfP91ltvUVhYyObNm7n66qtZsWLFfJ/aWQuNRhON1maJKBnMA7wz34CS+Y6Swexw+eWXc/nll8/3aZy1+NiUFhcSPg46hSjOPkTJIIooogCiZDAviGa+p8aPfvQjCgoKWLduHevWreO1115TfvfTn/6UsrIyli1bxptvvqm8HpUihwEzaGiIYpaYynchio9n45gP5qVRKZpAnAdEM98zQ6DGMSCakA0DosuEecLll19OU1MTra2t3HfffTPeT1FREatXr2bdunVs2rQJgOHhYS6++GKWLl3KxRdfrDReCSH4+te/TllZGWvWrOHo0aNhuZZIINo4NveIksFHAHv27OHYsWOK3+BDDz3EhRdeSHNzMxdeeKGyhn799ddpbm6mubmZ3/zmN3z5y1+et3OONo4tPESXCR9BvPjii8rI8Ntuu43du3fz8MMP8+KLL3LrrbciSRLbtm1jdHR0Ul/GXCLaOLbwEI0MznLI8wE2btzIb37zG4CQ5wMsNOj1euXfvo1jzz77LHa7nfb2dqVxbPPmzUrjmMPh4Nlnn+Xqq6+er9M/ezFfmcvof+H5Dyg4/f9soAbYBYz6bDNy+v+vAOd4vf5PYNN8X4Ofa/o/oBY4DrwE5Hn97j6gFWgEPuH1+uVA0+nf3Tff13A2/jeTRqUoFigkSfoRYAK+COwWQuglScoD9gohlkmS9MTpfz9zevtGebt5O+koFgyiy4SzGJIkJUiSlCT/G7gEqGPiaXrb6c1uA2TLn5eAW6UJbAPGokQQhYxoAvHsRg7wt9PmphrgaSHEG5IkHQaekyTpC8Ap4IbT27/GRDjdAliAO+b+lKNYqIguE6KIIgogukyIIoooTiNKBlFEEQUQJYMooojiNKJkEEUUUQBRMogiiihOI0oGUUQRBRAlgyiiiOI0omQQRRRRAPD/A/O2V1QfvXXBAAAAAElFTkSuQmCC\n", "text/plain": [ "