{
"cells": [
{
"attachments": {},
"cell_type": "markdown",
"metadata": {},
"source": [
"# All distances between two selections\n",
"\n",
"Here we use ``distances.distance_array`` to quantify the distances between each atom of a target set to each atom in a reference set, and show how we can extend that to calculating the distances between the centers-of-mass of residues.\n",
"\n",
"**Last updated:** December 2022 with MDAnalysis 2.4.0-dev0\n",
"\n",
"**Minimum version of MDAnalysis:** 0.19.0\n",
"\n",
"**Packages required:**\n",
" \n",
"* MDAnalysis (Michaud-Agrawal *et al.*, 2011, Gowers *et al.*, 2016)\n",
"* MDAnalysisTests\n",
" \n",
"**Optional packages for visualisation:**\n",
"\n",
"* [matplotlib](https://matplotlib.org)\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"execution": {
"iopub.execute_input": "2021-05-19T05:57:18.511813Z",
"iopub.status.busy": "2021-05-19T05:57:18.511175Z",
"iopub.status.idle": "2021-05-19T05:57:19.158622Z",
"shell.execute_reply": "2021-05-19T05:57:19.159054Z"
}
},
"outputs": [],
"source": [
"import MDAnalysis as mda\n",
"from MDAnalysis.tests.datafiles import PDB_small\n",
"from MDAnalysis.analysis import distances\n",
"\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline\n",
"\n",
"import warnings\n",
"# suppress some MDAnalysis warnings when writing PDB files\n",
"warnings.filterwarnings('ignore')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Loading files"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The test files we will be working with here feature adenylate kinase (AdK), a phosophotransferase enzyme. (Beckstein *et al.*, 2009) AdK has three domains: \n",
"\n",
" * CORE\n",
" * LID: an ATP-binding domain (residues 122-159)\n",
" * NMP: an AMP-binding domain (residues 30-59)\n",
" "
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"execution": {
"iopub.execute_input": "2021-05-19T05:57:19.163708Z",
"iopub.status.busy": "2021-05-19T05:57:19.163183Z",
"iopub.status.idle": "2021-05-19T05:57:19.318656Z",
"shell.execute_reply": "2021-05-19T05:57:19.319052Z"
}
},
"outputs": [],
"source": [
"u = mda.Universe(PDB_small) # open AdK"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Calculating atom-to-atom distances between non-matching coordinate arrays"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We select the alpha-carbon atoms of each domain."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"execution": {
"iopub.execute_input": "2021-05-19T05:57:19.323326Z",
"iopub.status.busy": "2021-05-19T05:57:19.322771Z",
"iopub.status.idle": "2021-05-19T05:57:19.326254Z",
"shell.execute_reply": "2021-05-19T05:57:19.326616Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"LID has 38 residues and NMP has 30 residues\n"
]
}
],
"source": [
"LID_ca = u.select_atoms('name CA and resid 122-159')\n",
"NMP_ca = u.select_atoms('name CA and resid 30-59')\n",
"\n",
"n_LID = len(LID_ca)\n",
"n_NMP = len(NMP_ca)\n",
"print('LID has {} residues and NMP has {} residues'.format(n_LID, n_NMP))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"`distances.distance_array`([API docs](https://docs.mdanalysis.org/stable/documentation_pages/analysis/distances.html#MDAnalysis.analysis.distances.distance_array)) will produce an array with shape `(n, m)` distances if there are `n` positions in the reference array and `m` positions in the other configuration. If you want to calculate distances following the minimum image convention, you *must* pass the universe dimensions into the ``box`` keyword."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"execution": {
"iopub.execute_input": "2021-05-19T05:57:19.332303Z",
"iopub.status.busy": "2021-05-19T05:57:19.331707Z",
"iopub.status.idle": "2021-05-19T05:57:19.333831Z",
"shell.execute_reply": "2021-05-19T05:57:19.334305Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"(38, 30)"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"dist_arr = distances.distance_array(LID_ca.positions, # reference\n",
" NMP_ca.positions, # configuration\n",
" box=u.dimensions)\n",
"dist_arr.shape"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Plotting distance as a heatmap"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"execution": {
"iopub.execute_input": "2021-05-19T05:57:19.362840Z",
"iopub.status.busy": "2021-05-19T05:57:19.351895Z",
"iopub.status.idle": "2021-05-19T05:57:19.476555Z",
"shell.execute_reply": "2021-05-19T05:57:19.476918Z"
},
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, 'Distance (Angstrom)')"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcoAAAHFCAYAAAByyrkJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABpCklEQVR4nO3deXhTVfoH8G+apGm6Qlu6QWnZtwIuKIsO+1YqiCCCIBbGccVRBkZHVKQMAwUUFEVwUNYRBBHKoCBQZXUULGWKuKD0RwuFFgql+5JmOb8/mEZC29zTJtAl3w9Pnockb849d0lPzrnnvlclhBAgIiKiKrnVdQWIiIjqMzaUREREdrChJCIisoMNJRERkR1sKImIiOxgQ0lERGQHG0oiIiI72FASERHZwYaSiIjIDjaUDlq3bh1UKpX14eHhgZCQEAwYMADx8fHIzs6u9Jm4uDioVKoaLaekpARxcXE4ePCgk2pet9LT06FSqfDWW285rczMzEzExcUhJSXFaWU2VBXHZXp6eo0/eyv2TX2gUqnw/PPP13U1qAFiQ+kka9euxXfffYfExES8//77uOOOO7Bo0SJ06tQJX331lU3sn/70J3z33Xc1Kr+kpARz585tNA3lrZCZmYm5c+eyoSQip9LUdQUai6ioKPTo0cP6fOzYsfjLX/6C+++/H2PGjMGZM2cQHBwMAGjRogVatGhRV1UlcimlpaXQ6/V1XQ1qwNijvIVatmyJJUuWoLCwEP/85z+tr1c19Lp//370798fAQEB0Ov1aNmyJcaOHYuSkhKkp6ejWbNmAIC5c+dah3mnTJkCAEhNTcXUqVPRrl07eHp6onnz5hg5ciROnTpls4yDBw9CpVLhk08+wWuvvYawsDD4+vpi8ODB+PXXXyvVf8+ePRg0aBD8/Pzg6emJTp06IT4+3ibm+PHjGDVqFPz9/eHh4YE777wTn376qfQ2slgsmD9/Plq2bAkPDw/06NEDX3/9daW4M2fOYOLEiQgKCoJOp0OnTp3w/vvv26zbPffcAwCYOnWqdRvFxcVh165dUKlUSEpKssZv27YNKpUKMTExNsvp1q0bxo4da30uhMCKFStwxx13QK/Xo2nTpnj44Ydx9uzZSnX86quvMGjQIPj6+sLT0xP33XdfpXWp2Pc//fQTHn30Ufj5+SE4OBh//OMfkZ+fr7i9EhMT8eCDD6JFixbw8PBA27Zt8fTTT+Pq1auKn+3fvz+ioqJw5MgR9OrVC3q9Hs2bN8fs2bNhNpur/MzSpUvRqlUreHt7o3fv3jh69KjN+8ePH8eECRMQGRkJvV6PyMhIPProozh37pxifSrk5eVh5syZaN26NXQ6HYKCgjBixAicPn3aGjN37lz07NkT/v7+8PX1xV133YXVq1fj5ns6REZG4oEHHsD27dtx5513wsPDA3PnzrWJ+ec//4n27dtDp9Ohc+fO2Lx5c6U6/fjjj3jwwQfRtGlTeHh44I477sD69ettYmr6faIGTJBD1q5dKwCIpKSkKt8vKioSarVaDBo0yPranDlzxI2bPi0tTXh4eIghQ4aIHTt2iIMHD4qNGzeKyZMni9zcXFFWVib27NkjAIgnnnhCfPfdd+K7774TqampQgghDh06JGbOnCk+++wzcejQIZGQkCBGjx4t9Hq9OH36tHU5Bw4cEABEZGSkmDRpkti1a5f45JNPRMuWLUW7du2EyWSyxn700UdCpVKJ/v37i02bNomvvvpKrFixQjz33HPWmP379wt3d3fxhz/8QWzZskXs2bNHTJkyRQAQa9eutbvd0tLSBAARHh4u7r//frFt2zaxdetWcc899witViu+/fZba+xPP/0k/Pz8RNeuXcWGDRvEvn37xMyZM4Wbm5uIi4sTQgiRn59v3Revv/66dRtlZGSIwsJCodVqxYIFC6xlPvPMM0Kv1wsvLy9RXl4uhBDi8uXLQqVSiRUrVljjnnzySaHVasXMmTPFnj17xKZNm0THjh1FcHCwuHTpkjXuX//6l1CpVGL06NFi+/bt4vPPPxcPPPCAUKvV4quvvqq07zt06CDeeOMNkZiYKJYuXSp0Op2YOnWq3W0mhBArV64U8fHxYufOneLQoUNi/fr1onv37qJDhw7W9RDi9+MyLS3N+lq/fv1EQECACAsLE++++67Yu3eveOGFFwQAMW3atEr7JjIyUgwfPlzs2LFD7NixQ3Tt2lU0bdpU5OXlWWO3bt0q3njjDZGQkCAOHTokNm/eLPr16yeaNWsmrly5org+BQUFokuXLsLLy0v8/e9/F3v37hXbtm0TL774oti/f781bsqUKWL16tUiMTFRJCYminnz5gm9Xi/mzp1rU15ERIQIDQ0VrVu3FmvWrBEHDhwQ33//vRBCWI+3zp07i08++UTs3LlTDB8+XAAQW7dutZZx+vRp4ePjI9q0aSM2bNggdu3aJR599FEBQCxatMgaV5PvEzVsbCgdpNRQCiFEcHCw6NSpk/X5zQ3lZ599JgCIlJSUasu4cuWKACDmzJmjWCeTySTKy8tFu3btxF/+8hfr6xVf7BEjRtjEf/rppwKA+O6774QQQhQWFgpfX19x//33C4vFUu1yOnbsKO68805hNBptXn/ggQdEaGioMJvN1X624o9xWFiYKC0ttb5eUFAg/P39xeDBg62vDRs2TLRo0ULk5+fblPH8888LDw8Pce3aNSGEEElJSdU20vfff78YOHCg9Xnbtm3FSy+9JNzc3MShQ4eEEEJs3LhRABC//fabEEKI7777TgAQS5YssSkrIyND6PV68fLLLwshhCguLhb+/v5i5MiRNnFms1l0795d3HvvvdbXKvb94sWLbWKfe+454eHhYXd738xisQij0SjOnTsnAIh///vf1veqayhvjhPi+o8BNzc3ce7cOSHE7/uma9euNn/sv//+ewFAfPLJJ9XWyWQyiaKiIuHl5SWWLVumuA5///vfBQCRmJgou9rCbDYLo9Eo/v73v4uAgACbbRYRESHUarX49ddfK30OgNDr9TY/cEwmk+jYsaNo27at9bUJEyYInU4nzp8/b/P56Oho4enpaf2hIPt9ooaPQ6+3gVC45ecdd9wBd3d3PPXUU1i/fn2Vw3r2mEwmLFiwAJ07d4a7uzs0Gg3c3d1x5swZ/PLLL5XiR40aZfO8W7duAGAdLvv2229RUFCA5557rtrZuampqTh9+jQmTZpkrUPFY8SIEcjKypIafhozZgw8PDysz318fDBy5EgcPnwYZrMZZWVl+Prrr/HQQw/B09Oz0nLKysoqDQdWZdCgQfjPf/6D0tJSnDt3DqmpqZgwYQLuuOMOJCYmArg+dNqyZUu0a9cOAPDFF19ApVLhscces1luSEgIunfvbp1Y9e233+LatWuIjY21ibNYLBg+fDiSkpJQXFxsU5+q9kFZWVmVs6RvlJ2djWeeeQbh4eHQaDTQarWIiIgAgCr39c18fHwqLXvixImwWCw4fPiwzesxMTFQq9U2dQRgM6xaVFSEv/3tb2jbti00Gg00Gg28vb1RXFxsU58bt4vJZLJ+J7788ku0b98egwcPtlvv/fv3Y/DgwfDz84NarYZWq8Ubb7yBnJycStusW7duaN++fZXlDBo0yDpXAADUajXGjx+P1NRUXLhwwbqsQYMGITw83OazU6ZMQUlJSaWJeErfJ2r42FDeYsXFxcjJyUFYWFi1MW3atMFXX32FoKAgTJs2DW3atEGbNm2wbNkyqWXMmDEDs2fPxujRo/H555/j2LFjSEpKQvfu3VFaWlopPiAgwOa5TqcDAGvslStXAMDuhKPLly8DAP76179Cq9XaPJ577jkAkDpvFhISUuVr5eXlKCoqQk5ODkwmE957771KyxkxYoT0cgYPHgyDwYBvvvkGiYmJCAwMxJ133onBgwdbZyV//fXXNn+wL1++DCEEgoODKy376NGj1uVWbIuHH364UtyiRYsghMC1a9ds6qO0D6pisVgwdOhQbN++HS+//DK+/vprfP/999YfCvY+W+HGRqJCxT7IycmpcR0nTpyI5cuX409/+hP27t2L77//HklJSWjWrJlN3M3bpeJ835UrVxQntn3//fcYOnQoAODDDz/Ef/7zHyQlJeG1116rcr1DQ0OrLau64+3G9c/JyamyjIrvcG22EzVsnPV6i+3atQtmsxn9+/e3G/eHP/wBf/jDH2A2m3H8+HG89957mD59OoKDgzFhwgS7n/3444/x+OOPY8GCBTavX716FU2aNKlxnSsmDlX8wq5KYGAgAGDWrFkYM2ZMlTEdOnRQXNalS5eqfM3d3R3e3t7QarVQq9WYPHkypk2bVmUZrVq1UlxOz5494e3tja+++grp6ekYNGgQVCoVBg0ahCVLliApKQnnz5+3aSgDAwOhUqlw5MgR6x+/G1W8VrEt3nvvPfTq1avK5VfVQNXUjz/+iJMnT2LdunWIjY21vp6amipdRkWjfqOKfXDzH3wl+fn5+OKLLzBnzhy88sor1tcNBkOlHwY3TqQCft9nzZo1s3ucAcDmzZuh1WrxxRdf2Iw+7Nixo8p4e9coV3e8Ab+vf0BAALKysirFZWZmAvh9f5PrYEN5C50/fx5//etf4efnh6efflrqM2q1Gj179kTHjh2xceNGnDhxAhMmTLD7K1WlUlX6Q75r1y5cvHgRbdu2rXG9+/TpAz8/P3zwwQeYMGFClX94OnTogHbt2uHkyZOVGuia2L59O958803rH8DCwkJ8/vnn+MMf/gC1Wg1PT08MGDAA//3vf9GtWze4u7tXW5a9baTVatG3b18kJiYiIyMDCxcuBHD9B4pGo8Hrr79ubTgrPPDAA1i4cCEuXryIRx55pNrl3nfffWjSpAl+/vnnW3pBe8V+uHlf3zijWklhYSF27txpM1y4adMmuLm5oW/fvjWujxCiUn0++uijSrNob7x06kbR0dF44403sH//fgwcOLDa5Wg0Gpth4NLSUvzrX/+qUX2B66MGly9ftv5wMZvN2LJlC9q0aWPt2Q4aNAgJCQnIzMy0GQnasGEDPD09q/0xRI0XG0on+fHHH63nX7Kzs3HkyBGsXbsWarUaCQkJ1l5aVT744APs378fMTExaNmyJcrKyrBmzRoAsPZwfHx8EBERgX//+98YNGgQ/P39ERgYaJ0Ov27dOnTs2BHdunVDcnIy3nzzzVpfq+nt7Y0lS5bgT3/6EwYPHownn3wSwcHBSE1NxcmTJ7F8+XIA1/9AR0dHY9iwYZgyZQqaN2+Oa9eu4ZdffsGJEyewdetWxWWp1WoMGTIEM2bMgMViwaJFi1BQUGAzpX/ZsmW4//778Yc//AHPPvssIiMjUVhYiNTUVHz++efYv38/gOtD2Hq9Hhs3bkSnTp3g7e2NsLAw6x+7QYMGYebMmTbbVa/Xo0+fPti3bx+6deuGoKAg63Lvu+8+PPXUU5g6dSqOHz+Ovn37wsvLC1lZWfjmm2/QtWtXPPvss/D29sZ7772H2NhYXLt2DQ8//DCCgoJw5coVnDx5EleuXMHKlStrtS9u1LFjR7Rp0wavvPIKhBDw9/fH559/bj3HKiMgIADPPvsszp8/j/bt22P37t348MMP8eyzz6Jly5Y1qo+vry/69u2LN99803osHjp0CKtXr5YeyZg+fTq2bNmCBx98EK+88gruvfdelJaW4tChQ3jggQcwYMAAxMTEYOnSpZg4cSKeeuop5OTk4K233qqyl68kMDAQAwcOxOzZs+Hl5YUVK1bg9OnTNpeIzJkzB1988QUGDBiAN954A/7+/ti4cSN27dqFxYsXw8/Pr8bLpQauLmcSNQYVswsrHu7u7iIoKEj069dPLFiwQGRnZ1f6zM2zXr/77jvx0EMPiYiICKHT6URAQIDo16+f2Llzp83nvvrqK3HnnXcKnU4nAIjY2FghhBC5ubniiSeeEEFBQcLT01Pcf//94siRI6Jfv36iX79+1s9XzNK7cSq8EL/Pcrx5tuju3btFv379hJeXl/D09BSdO3e2mR4vhBAnT54UjzzyiAgKChJarVaEhISIgQMHig8++MDudqtY5qJFi8TcuXNFixYthLu7u7jzzjvF3r17q4z/4x//KJo3by60Wq1o1qyZ6NOnj/jHP/5hE/fJJ5+Ijh07Cq1WW2mW8MmTJwUA0a5dO5vPzJ8/XwAQM2bMqLKua9asET179hReXl5Cr9eLNm3aiMcff1wcP37cJu7QoUMiJiZG+Pv7C61WK5o3by5iYmJstnfFvr/50omqZqlW5eeffxZDhgwRPj4+omnTpmLcuHHi/Pnzlda1ulmvXbp0EQcPHhQ9evQQOp1OhIaGildffdVm5nLFvnnzzTcrLf/m5Vy4cEGMHTtWNG3aVPj4+Ijhw4eLH3/8UURERFiPTyW5ubnixRdfFC1bthRarVYEBQWJmJgYm0ub1qxZIzp06CB0Op1o3bq1iI+PF6tXr660jhERESImJqbK5eB/l8GsWLFCtGnTRmi1WtGxY0excePGSrGnTp0SI0eOFH5+fsLd3V1079690vejpt8narhUQihMySSiRqF///64evUqfvzxx7quClGDwlmvREREdrChJCIisoNDr0RERHawR0lERGQHG0oiIiI72FASERHZwYQDuJ5DMzMzEz4+PnbTXxER3WpCCBQWFiIsLAxubremL1NWVoby8nKnlOXu7m6TWrAxYkOJ6zkcb75TABFRXcrIyKh1di17ysrK0CrCG5eyq75Zd02FhIQgLS2tUTeWddpQHj58GG+++SaSk5ORlZWFhIQEjB49GgBgNBrx+uuvY/fu3Th79iz8/PwwePBgLFy40JqS7Nq1a5gzZw727duHjIwMBAYGYvTo0Zg3b16N0kz5+PgAALqPeR1qrf2dXdZU+ReeLk9uInG5j1zvtaS5XHkWib3pZpQqCkZ/k1Rck5BCqTiVSnkdtG4WqbLuDLwoFdfeq3Ji6yqXq1L+g9FKa//2VxWOFreTivu1WDlJejcf+8nCK5Raqs9/e6PEzI5ScW5Q3ldtm16RKis1t/rUjTfKzmoiFadvqnxHDh8Pg9wyLzSVioO73HGpzlP+Apq9lI81S1kZMv8Wb/275Gzl5eW4lG3GueRI+Po41mMtKLQg4u50lJeXs6G8VYqLi9G9e3dMnToVY8eOtXmvpKQEJ06cwOzZs9G9e3fk5uZi+vTpGDVqFI4fPw7gek8wMzMTb731Fjp37oxz587hmWeeQWZmJj777DPpelQMt6q1HlC729/Zap3ygaV2l2vY1O5yDaWbh+QVPDINpVo5BgDc9HINpdpTbvhGpqFUq+X+ILl7a6Xi9F5yh7dWYrjdy11uw+lUcnXTQrlx85BcT4tFLk7tKZcb1U3mR42XXOOsLpdcpl7uj6zaU/kYUeulipJepmxD6VamfLwJvXwv7lafBvL2UcFb8sd6dSxwjVNVddpQRkdHIzo6usr3/Pz8KiV7fu+993Dvvffi/PnzaNmyJaKiorBt2zbr+23atMH8+fOtN9rVaDiyTERUFbOwwOzgVfRmIfcjoqFrUC1Jfn4+VCqV3TsT5Ofnw9fX124jaTAYYDD8PjxTUFDgzGoSEdV7FghYJIbZlcpwBQ3m8pCysjK88sormDhxInx9fauMycnJwbx58xTv/RgfHw8/Pz/rgxN5iIioOg2ioTQajZgwYQIsFgtWrFhRZUxBQQFiYmLQuXNnzJkzx255s2bNQn5+vvWRkZFxK6pNRFRvWZz0zxXU+6FXo9GIRx55BGlpadi/f3+VvcnCwkIMHz4c3t7eSEhIgFZrf3KDTqer1U1fiYgaC7MQMDuY6tvRzzcU9bqhrGgkz5w5gwMHDiAgIKBSTEFBAYYNGwadToedO3c26inKRER0+9VpQ1lUVITU1FTr87S0NKSkpMDf3x9hYWF4+OGHceLECXzxxRcwm824dOkSAMDf3x/u7u4oLCzE0KFDUVJSgo8//hgFBQXWiTnNmjWDWi15LQQRkYvhZB55ddpQHj9+HAMGDLA+nzFjBgAgNjYWcXFx2LlzJwDgjjvusPncgQMH0L9/fyQnJ+PYsWMAgLZt29rEpKWlITIyskb10V0zQ6N1PFuF/qpcGdoSuVPEZr1cgy+TcEBTIlUUykxyh0auW9UTq2pF4vo9ADguWdzVci+puHB9rmKMRchdL3Yyv7lUXFpu5dGRm5kscseHySJ3fFzNkbyAXWI3yG6PnFxvqTjtVbnjrcSsvE/LvOWu8dRek9tuZg+5/aAtUI5TmZS3m0XiekxnsEDAzIZSSp02lP3794e922Eq3SpT6fNERESOqtfnKImI6Nbg0Ks8NpRERC6Is17lNYjrKImIiOoKe5RERC7I8r+Ho2W4AjaUREQuyOyEWa+Ofr6hYENJROSCzAJOuHuIc+pS3/EcJRERkR3sUd7A0FQNk8INesuaSty4We5exjDp5C7aNspdNw8hcf20SjKfglnumm3pGzwLs/K6Covc9nBXy62El+SOCNQWKcaEaPOlygrTy92yLbtE+eJ/P22ZVFmydHqjVJzJqHwgeWjk9ru7Ti7OKHmzc5mkFLI3AJe83zUsernyTBLJBEw+yseuxQlJT2TwHKU8NpRERC7IAhXMkPtxaq8MV8ChVyIiIjvYoyQickEWcf3haBmugA0lEZELMjth6NXRzzcUHHolIiKygw0lEZELquhROvqorfj4eKhUKkyfPt362pQpU6BSqWwevXr1csLaOoZDr0RELsgiVNL3FbVXRm0kJSVh1apV6NatW6X3hg8fjrVr11qfu7tLXqt2C7FHSUREt01RUREmTZqEDz/8EE2bNq30vk6nQ0hIiPXh7+9fB7W0xYaSiMgFOXPotaCgwOZhMBiqXe60adMQExODwYMHV/n+wYMHERQUhPbt2+PJJ59Ednb2LVn/muDQ6w00ZQIaheSF5mLl+dDuBXKZNTQauWGLsgC53SQkfvZ4XJPNgCJXt2JPuWERN6NyeSqJzCYAkOneRCrOaJH7HVhqVk7RUiaZxiW1MFAqLjvPWypOhlYyU1FZnodUnMqgvN0uq5UzCwFAeaHc8aHPk9tXJonjyCg5HKjPl4srh0TKKwDqUonsU27K66kquz39FzPcYHawr1Rx5IWHh9u8PmfOHMTFxVWK37x5M06cOIGkpKQqy4uOjsa4ceMQERGBtLQ0zJ49GwMHDkRycjJ0Op1DdXUEG0oiIhcknHCOUvzv8xkZGfD19bW+XlWjlpGRgRdffBH79u2Dh0fVP9rGjx9v/X9UVBR69OiBiIgI7Nq1C2PGjHGoro5gQ0lERA7x9fW1aSirkpycjOzsbNx9993W18xmMw4fPozly5fDYDBArbbtvYeGhiIiIgJnzpy5JfWWxYaSiMgF3e6EA4MGDcKpU6dsXps6dSo6duyIv/3tb5UaSQDIyclBRkYGQkNDHaqno9hQEhG5ILNwg1lmYoPdMuRjfXx8EBUVZfOal5cXAgICEBUVhaKiIsTFxWHs2LEIDQ1Feno6Xn31VQQGBuKhhx5yqJ6OYkNJRER1Tq1W49SpU9iwYQPy8vIQGhqKAQMGYMuWLfDxkZs8dquwoSQickEWqGBxcNarBY5lRT948KD1/3q9Hnv37nWovFuFDSURkQtiUnR5TDhARERkB3uUREQuyDmTeVzjhpRsKG9QGuAGtbv9A6e8ifJQg6ZMLpOH7LW+5fYvT/qdRHlqg9xCTXq5RQq9RSrOopbISFJ91isbWq1JKs5Ta5SKC9QVKcaEuudJldXcM18q7rJeeXJCE49SqbK8tXIb7qJ3E6k4o5tyFiJPj3KpsmT/jpq85DIfmXXKBar0cseHyVvue2r2lSvPolM+xoWXchYlS6nc8hx1/Rylg0nROfRKRERE7FESEbkgixNyvTo667WhYENJROSCeI5SHhtKIiIXZIFbnV9H2VDwHCUREZEd7FESEbkgs1DB7OBtthz9fEPBhpKIyAU558bNHHolIiJyeexREhG5IItwg8XBWa8Wznq99Q4fPow333wTycnJyMrKQkJCAkaPHg0AMBqNeP3117F7926cPXsWfn5+GDx4MBYuXIiwsDAAQHp6Olq1alVl2Z9++inGjRtXo/qoywGlXB2aEuVy3AvkstUIya2vLZY7mFXKST+gy5PMpKOQoaiCsUAuu4m6VPlchmzWoFIfD6m4LDe5L7FO7bxMKGcLA6Ti8vM8FWOMZrltq5PMQGS6JrfdtIXK+z7PInfbI1Wp3HHkc1lu31u0ynGGMp1UWZ5ZcsssK5f7oroZlcszlilvD0uZxBfZCTj0Kq9Oh16Li4vRvXt3LF++vNJ7JSUlOHHiBGbPno0TJ05g+/bt+O233zBq1ChrTHh4OLKysmwec+fOhZeXF6Kjo2/nqhARUSNVpz3K6Ojoahs0Pz8/JCYm2rz23nvv4d5778X58+fRsmVLqNVqhISE2MQkJCRg/Pjx8Pb2vmX1JiJq6CxwfNaq3PhUw9egzlHm5+dDpVKhSZMmVb6fnJyMlJQUvP/++7e3YkREDYxzEg64xnzQBtNQlpWV4ZVXXsHEiRPh61v17TRWr16NTp06oU+fPnbLMhgMMBh+v+NCQUGBU+tKRESNR4P4OWA0GjFhwgRYLBasWLGiypjS0lJs2rQJTzzxhGJ58fHx8PPzsz7Cw8OdXWUionqtIterow9XUO/X0mg04pFHHkFaWhoSExOr7U1+9tlnKCkpweOPP65Y5qxZs5Cfn299ZGRkOLvaRET1WsX9KB19uIJ6PfRa0UieOXMGBw4cQEBA9VPvV69ejVGjRqFZs2aK5ep0Ouh0clPIiYgaI+fcPaTe97Wcok4byqKiIqSmplqfp6WlISUlBf7+/ggLC8PDDz+MEydO4IsvvoDZbMalS5cAAP7+/nB3d7d+LjU1FYcPH8bu3btv+zoQEVHjVqcN5fHjxzFgwADr8xkzZgAAYmNjERcXh507dwIA7rjjDpvPHThwAP3797c+X7NmDZo3b46hQ4c6VB+TByDclWMUy9HLDUeoJK/VFXLXnTt1mTLJCwAAFrl1dTMpx7mVyy1SZZBMhiB5objBrBwnm8HETXqnKm8Pi+y2lRz9EpIJGJxJaOSWaXHiXyLZ74tZclDJLJenARZ35XU165UvqLDcposunJNwgD3KW65///4QdlIg2XvvRgsWLMCCBQucVS0iokbPIlSwOHodpYvcPcQ1fg4QERHVUr2ezENERLeGxQlDr0w4QEREjZZz7h7iGg2la6wlERFRLbFHSUTkgsxQwexgwgBHP99QsKEkInJBHHqV5xprSUREVEvsURIRuSAzHB86lc1L0tCxobyB5xULNFr7WTHKJLLC6K8apZbnVi6XgcPgq5eKk7n21+Oacw9to7fcoITaoByjLZJLMCG7zHKNQpql/7ns4aMYo9fI7dPsQskbhhcof/XKJLP8CMmLvtXFkvuqVLk8o1Eya1C5ZEYjyaxMQuYvlpMTEKkkE+WoJA4Rdany9lCV3Z6BPg69ymNDSUTkgpgUXZ5rrCUREVEtsUdJROSChBPuJyl4eQgRETVWHHqV5xprSUREVEvsURIRuSDeZkseG0oiIhfEGzfLc421JCIiqiX2KImIXBCHXuWxobyB18USaNT203BoSjwUy3G/kCe3QCGXQsTTXy7DjFArx7hfK5Mqy6KRywbkXiCZecWovK6aUqmioC6TzAojkQUFAEqLlbfvZXfl7D0AUJQnud0KletWrpH7epYZ5dZTnycXpylRjjHpJQ42ABrJfaXLk/sumJW/fjDpJZeZK5nCRyWZ+Uji+DVLHB5mw+1pfCxwc/jGy65y42bXWEsiIqJaYo+SiMgFmYUKZgeHTh39fEPBhpKIyAXxHKU8NpRERC5IOOHuIYKZeYiIiIg9SiIiF2SGygk3bubQKxERNVIW4fg5RouTb5JdX3HolYiIyA72KImIXJDFCZN5HP18Q8GG8gZGH3cIjf0sLWVNlTOSaIJ9pZanMsuNW5T6y2VBsUgl8JHLHFPWRG6ZZYFSYYBFeYhHNjOPIdB+9qQKqiCDVFxg00LFmPZNs6XK+kUqCrhq9FOM8QyQSJEDwFNnlIrLMTSVinMrlRiOC5XbtmUlcn9iVEatVJzZU/k7Ux4otz0s7nJ1K/eTO940JcrbzeKuXH9LmdzyHGVxwo2bHf18Q+EaPweIiIhqiT1KIiIXxMw88thQEhG5IJ6jlOcaa0lERPVKfHw8VCoVpk+fbn1NCIG4uDiEhYVBr9ejf//++Omnn+qukv/DhpKIyAVZoLLme631o5aTeZKSkrBq1Sp069bN5vXFixdj6dKlWL58OZKSkhASEoIhQ4agsFB5wt2txIaSiMgFif/NenXkIWrRUBYVFWHSpEn48MMP0bTp7zOxhRB455138Nprr2HMmDGIiorC+vXrUVJSgk2bNjlz1WuMDSURkQtyuDd5w91HCgoKbB4GQ/WXD02bNg0xMTEYPHiwzetpaWm4dOkShg4dan1Np9OhX79++Pbbb2/NRpDEhpKIiBwSHh4OPz8/6yM+Pr7KuM2bN+PEiRNVvn/p0iUAQHBwsM3rwcHB1vfqCme93kBdYoJaY7IbozEobzI3k9wFw0IlN2yhcmI+RbNWbpkWueu/YdE4r3ImyWEci15u++olL8TXa5XjdG5mqbJUsjtLYrvptPaPxQpqN9njTa5uFoljRKOW2x5mlVziCtmf7FLHpeQxKXtlg0ry+n83icNN5juvMt6eSy6cOes1IyMDvr6/J1rR6XSVYjMyMvDiiy9i37598PDwqLZM1U3bSAhR6bXbrU57lIcPH8bIkSMRFhYGlUqFHTt22LwfFxeHjh07wsvLC02bNsXgwYNx7Ngx6/vXrl3Dn//8Z3To0AGenp5o2bIlXnjhBeTn59/mNSEialicOfTq6+tr86iqoUxOTkZ2djbuvvtuaDQaaDQaHDp0CO+++y40Go21J3lz7zE7O7tSL/N2q9OGsri4GN27d8fy5curfL99+/ZYvnw5Tp06hW+++QaRkZEYOnQorly5AgDIzMxEZmYm3nrrLZw6dQrr1q3Dnj178MQTT9zO1SAiIgWDBg3CqVOnkJKSYn306NEDkyZNQkpKClq3bo2QkBAkJiZaP1NeXo5Dhw6hT58+dVjzOh56jY6ORnR0dLXvT5w40eb50qVLsXr1avzwww8YNGgQoqKisG3bNuv7bdq0wfz58/HYY4/BZDJBo+HIMhFRVW53rlcfHx9ERUXZvObl5YWAgADr69OnT8eCBQvQrl07tGvXDgsWLICnp2eltuB2azAtSXl5OVatWgU/Pz9079692rj8/Hz4+vrabSQNBoPNrKyCggKn1pWIqL67cejUkTKc6eWXX0ZpaSmee+455ObmomfPnti3bx98fHycupyaqvcN5RdffIEJEyagpKQEoaGhSExMRGBg1besyMnJwbx58/D000/bLTM+Ph5z5869FdUlIiJJBw8etHmuUqkQFxeHuLi4OqlPder95SEDBgxASkoKvv32WwwfPhyPPPIIsrMr3/KooKAAMTEx6Ny5M+bMmWO3zFmzZiE/P9/6yMjIuFXVJyKql5w5maexq/cNpZeXF9q2bYtevXph9erV0Gg0WL16tU1MYWEhhg8fDm9vbyQkJECrtT+HXKfTVZqlRUTkSthQyqv3DeXNhBCVzi8OHToU7u7u2Llzp93rc4iIiGqqTs9RFhUVITU11fo8LS0NKSkp8Pf3R0BAAObPn49Ro0YhNDQUOTk5WLFiBS5cuIBx48YBuN6THDp0KEpKSvDxxx9b0ycBQLNmzaBWS17sTETkYurjZJ76qk4byuPHj2PAgAHW5zNmzAAAxMbG4oMPPsDp06exfv16XL16FQEBAbjnnntw5MgRdOnSBcD1C1grEhC0bdvWpuy0tDRERkbWqD7FLfTQaO33SEuClDvhZq2n1PKEZDte1ELuYDTplTOSGJrILbS8iVx2E2N49TkdbySTscZilKtbcEieVFwXf7m0V208ryjGdPc8J1WWl0Zue5x0b64Y07VpplRZOje5DD5fm+S2r9GsHNep2WWpsi4W+UnFZVoCpOK0vsrbt3mTIqmyLqr85ZbpWS4VZ/CqfJH9zVQeyhmNLKVyx5CjBGp2eUd1ZbiCOm0o+/fvDyGq39Tbt2936PNERFQ19ijlNbhzlERERLdTvb+OkoiInI89SnlsKImIXBAbSnkceiUiIrKDPUoiIhfEHqU8NpRERC5ICBWEgw2do59vKDj0SkREjc6NGdwcxYaSiMgFVdyP0tFHfbF3715MmTIFbdq0gVarhaenJ3x8fNCvXz/Mnz8fmZlyCTyqwqHXGwjV9YfdGImfFhb7OdmtzO5yB5lZOeHH9eW6K8eY5JIGweQpl8hBpzdKxanVFsUYi0Vue4R550vFtdRfk4pr4Z6jGNNMXShVVlNNiVycTjkuxN2590n1dJfbVwLKcYG6Yqmycg1yB5xMxhoA8NIrZ8nx1ZVJlZXlLrdMjUb52AUAo8Thq3JT/l7JZLFyhsZyjnLHjh3429/+hvz8fIwYMQIvvfQSmjdvDr1ej2vXruHHH3/EV199hXnz5mHKlCmYN28emjVrVqNlsKEkIqIGa8GCBXjrrbcQExMDN7fKPZlHHnkEAHDx4kUsW7YMGzZswMyZM2u0DDaUREQuqLFM5vn++++l4po3b47FixfXahlsKImIXFBjGXq9HdhQEhG5oMbSo7yREAKfffYZDhw4gOzsbFgstueXlW60UR02lERE1Ci8+OKLWLVqFQYMGIDg4GCoVM5pyNlQEhG5IOGEodf61qP8+OOPsX37dowYMcKp5bKhJCJyQQKAo7fzrW93A/bz80Pr1q2dXi4TDhARUaMQFxeHuXPnorS01KnlskdJROSCLFBB5WBmnfqUmQcAxo0bh08++QRBQUGIjIyEVmub/eXEiRO1KpcN5Q0KI9yg1tnvZJcFKmfpKPeV66ibdXIDF8aWcplGNDqTYkxZqVzaIJ+mchlm7g65IBXn7qZcN4tM2iMA/Zucloq7x+O8VFy6sYliTJhaLm/kI37HpeIidFcVY/6gPytVllnyj1WJTOomyE35v8/nN6myTuoinLZMAGjre0UxprVeedsCgMEs9+dPJosSAJzz8FeM8dcrl2UqNkDuyHVMY5z1OmXKFCQnJ+Oxxx7jZB4iIqKb7dq1C3v37sX999/v1HLZUBIRuSCLUEHVyBIOhIeHw9fX1+nlcjIPEZELEsI5j/pkyZIlePnll5Genu7UctmjJCKiRuGxxx5DSUkJ2rRpA09Pz0qTea5dk7uj0M3YUBIRuaDGOJnnnXfeuSXlsqEkInJBjbGhjI2NvSXlsqEkInJBjXEyDwCYzWbs2LEDv/zyC1QqFTp37oxRo0ZBrVbXukw2lERE1CikpqZixIgRuHjxIjp06AAhBH777TeEh4dj165daNOmTa3KZUNZQzLXxFvc5aaCWeSu/YebRq48nUTCAVleunKpOH9tsVSct0bugn0ZIZo8qbgwyV+QZpGvGNNC4y1VVpapSCouRKO8zFC1XIIAI8xSccHaAqk4rUr5OGqtkZsUkeMut91CPeXq1tYzWzGmne6yVFmnvUOk4rzUct+FXIOnYoyPu3LyEKNRbnmOcsas1fo26/WFF15AmzZtcPToUfj7X08AkZOTg8ceewwvvPACdu3aVaty2VASEbmg6w2lo+conVQZJzl06JBNIwkAAQEBWLhwIe67775al8vrKImIqFHQ6XQoLCys9HpRURHc3eVGaKrChpKIyAVVzHp19FGfPPDAA3jqqadw7NgxCCEghMDRo0fxzDPPYNSoUbUulw0lEZELEk561Cfvvvsu2rRpg969e8PDwwMeHh6477770LZtWyxbtqzW5fIcJRERNXhCCOTn5+OTTz5BZmYmfvnlFwgh0LlzZ7Rt29ahstlQEhG5oMaWcEAIgXbt2uGnn35Cu3btHG4cb8ShVyIiV9TIxl7d3NzQrl075OTkOL9sp5dIRET1nzMm8tSjHiUALF68GC+99BJ+/PFHp5bLoVciImoUKu4e0r17d7i7u0Ov19u8z7uHOIHHVQG1QlYdi0a5Ey6ZrAYWya1frNdJxRV5S2SiMcgNImQb5bLa/OgeJhUnk5HEZJFbplYll4kGOCMVFaaufN1VbZ0oD5SKO1kSoRjjr5b7VWwRcteHnSkNkoqT2b7h7nLDWxfKA6TirpTKZfC5pPNTjPFxUz7WACC/3EMqzqiROy4LDcrfU6NZuSxTifOyWNnTGDPzvP3221CpnN/LrdOG8vDhw3jzzTeRnJyMrKwsJCQkYPTo0db34+LisHnzZmRkZMDd3R1333035s+fj549ewIA0tPT0apVqyrL/vTTTzFu3LjbsRpERA1OY5vMAwBTpky5JeXW6TnK4uJidO/eHcuXL6/y/fbt22P58uU4deoUvvnmG0RGRmLo0KG4cuUKACA8PBxZWVk2j7lz58LLywvR0dG3c1WIiKiOqdVqZGdXzgeck5PTcO8eEh0dbbdBmzhxos3zpUuXYvXq1fjhhx8waNAgqNVqhITYJjZOSEjA+PHj4e0tN5RDROSSnDEZp571KEU1Y8EGg8GhFHYN5hxleXk5Vq1aBT8/P3Tv3r3KmOTkZKSkpOD999+3W5bBYIDB8Pt5gIICuTsXEBE1Fo3pHOW7774LAFCpVPjoo49sOkpmsxmHDx9Gx44da11+vW8ov/jiC0yYMAElJSUIDQ1FYmIiAgOrnjCxevVqdOrUCX369LFbZnx8PObOnXsrqktERLfZ22+/DeB6j/KDDz6wGWZ1d3dHZGQkPvjgg1qXX+8bygEDBiAlJQVXr17Fhx9+iEceeQTHjh1DUJDtDL7S0lJs2rQJs2fPVixz1qxZmDFjhvV5QUEBwsPDnV53IqJ6yxkJA+pJjzItLQ3A9fZi+/btaNq0qVPLr/cJB7y8vNC2bVv06tULq1evhkajwerVqyvFffbZZygpKcHjjz+uWKZOp4Ovr6/Ng4jIlTTGu4ccOHDAppE0m81ISUlBbm6uQ+XW+4byZkIIm/OLFVavXo1Ro0ahWbNmdVArIiKqa9OnT7d2pMxmM/r27Yu77roL4eHhOHjwYK3LrdOh16KiIqSmplqfp6WlISUlBf7+/ggICMD8+fMxatQohIaGIicnBytWrMCFCxcqXR+ZmpqKw4cPY/fu3bd7FYiIGq56MnTqLFu3bsVjjz0GAPj888+Rnp6O06dPY8OGDXjttdfwn//8p1bl1mlDefz4cQwYMMD6vOK8YWxsLD744AOcPn0a69evx9WrVxEQEIB77rkHR44cQZcuXWzKWbNmDZo3b46hQ4c6ViHV/x5KMUoscotzM8rFySaiEWblyrlJZuaxaOUOjSKj3JRri8SGKzVqpcrK9GwiFXdFJzek7qUqV4wpsSjHAIBRyF2WVGRWzuKSZ/aUKstD8kAySKaCMkj8WSi2yGWLMgq5a9cMZrm6FZj0ijFlQu44MknWzU0llylHq1b+osrEqCRinKExJhzIycmxXjK4e/dujBs3Du3bt8cTTzxhnRlbG3XaUPbv37/a614AYPv27VLlLFiwAAsWLHBWtYiIGr9GNJmnQnBwMH7++WeEhoZiz549WLFiBQCgpKSk4SYcICIicpapU6fikUceQWhoKFQqFYYMGQIAOHbsmEPXUTa4yTxEROQMKic95K1cuRLdunWzXm3Qu3dvfPnll9b3p0yZApVKZfPo1auXdPlxcXH46KOP8NRTT+E///kPdLrrpwjUajVeeeWVGtX1RuxREhG5ojoYem3RogUWLlyItm3bAgDWr1+PBx98EP/973+tc0+GDx+OtWvXWj9T09RzDz/8cKXXYmNja1bRm7ChJCKi22LkyJE2z+fPn4+VK1fi6NGj1oZSp9NVyuEtq7oJOyqVCh4eHmjbti369u1b4/OVbCiJiFyRE3uUN+fL1ul01mHP6pjNZmzduhXFxcXo3bu39fWDBw8iKCgITZo0Qb9+/TB//vxKmdiq8/bbb+PKlSsoKSlB06ZNIYRAXl4ePD094e3tjezsbLRu3RoHDhyoUTY2nqMkInJFFXcPcfSB67c89PPzsz7i4+OrXeypU6fg7e0NnU6HZ555BgkJCejcuTOA63eU2rhxI/bv348lS5YgKSkJAwcOrDLJTFUWLFiAe+65B2fOnEFOTg6uXbuG3377DT179sSyZctw/vx5hISE4C9/+UuNNhV7lERE5JCMjAybVKD2epMdOnRASkoK8vLysG3bNsTGxuLQoUPo3Lkzxo8fb42LiopCjx49EBERgV27dmHMmDGK9Xj99dexbds2tGnTxvpa27Zt8dZbb2Hs2LE4e/YsFi9ejLFjx9Zo/dhQEhG5IGfeZqsmObPd3d2tk3l69OiBpKQkLFu2DP/85z8rxYaGhiIiIgJnzpyRKjsrKwsmk6nS6yaTCZcuXQIAhIWFobCwUKq8Cmwob1ASrIJaZ3+6c1mQctYMs7vkiLbswHeQ3LCDp145e0ypwnmDCt6+ZVJxHZtUvpt4VXTqygfvzYpMcrPb7vZOl4q7y+OCVFwbjXK2lwKLcv0BoKv7Jam4ci/lyQR36PKkysqTzATVyStLKk4mm04b7RWpssoscllyfvGSm7wRrFO+d2yIJl+qrABdsVRckE7uj2qRUfm71dS9VDGmHHJZoBxWTxIOVJe/G7ieaScjIwOhoaFSZQ0YMABPP/00PvroI9x5550AgP/+97949tlnMXDgQADXh35btWpVozryHCUREd0Wr776Ko4cOYL09HScOnUKr732Gg4ePIhJkyahqKgIf/3rX/Hdd98hPT0dBw8exMiRIxEYGIiHHnpIqvzVq1fD398fd999t3VCUY8ePeDv729Nlu7t7Y0lS5bUqN7sURIRuaIbJuM4VEYNXL58GZMnT0ZWVhb8/PzQrVs37NmzB0OGDEFpaSlOnTqFDRs2IC8vD6GhoRgwYAC2bNkCHx8fqfJDQkKQmJiI06dP47fffoMQAh07dkSHDh2sMTfmF5fFhpKIyAWpxPWHo2XURFX3Eq6g1+uxd+9exyr0Px07dnQoZd3N2FASEbmienKO0pnMZjPWrVuHr7/+GtnZ2bBYbE/g79+/v1blsqEkIqJG4cUXX8S6desQExODqKgoqFTOuQ0YG0oiIldUB+cob7XNmzfj008/xYgRI5xaLhtKIiJX1AiHXm+8RtOZeHkIERE1CjNnzsSyZcsgHM2kcBP2KG9g9hIQHvY3sPBUTjhgKXXu7w+tu9zF7jqtclyZm9xF/W6S09n0aqNkecpXxVuE3HYzO/n3nUEob7d/F0dKlXXN5C0V18EjUzEmxyw3rLWzsLtUXKBG7sL5QouHYoxWYn8CgIeb3PHh7iZ3jMvQqpS/ozVZpp9aOUkAAAToShRj/N2VkxwYjHLbzGGNsEf5zTff4MCBA/jyyy/RpUsXaLW2CS+2b99eq3LZUBIRuaJG2FA2adJEOjlBTdSqoTxz5gz+/e9/Iz09HSqVCq1atcLo0aPRunVrZ9ePiIhIyo03fHamGo9hxcfHo3Pnzvjb3/6Gbdu2YevWrXjppZfQsWNHvPXWW7eijkRE5GxOvM1WfZabm4v33nsPd9xxR63LqFFDeeDAAbz++ut47bXXcPXqVWRlZeHSpUu4cuUKXnnlFbzyyis4fPhwrStDRES3R0VmHkcf9dVXX32FRx99FGFhYVi8eDH69etX67JqNPT6wQcf4E9/+hPi4uJsXvf398ff//53XLp0CStXrkTfvn1rXSEiIqLaOH/+PNauXYu1a9eiqKgIubm5+PTTT2t8/8mb1ahH+f3332Py5MnVvj958mQcPXrUoQoREdFtIJz0qAc+/fRTDB06FJ06dcKPP/6IZcuWITMzE25ubujUqZPD5deoR3n58mVERkZW+36rVq2sN8ckIiK6HSZOnIiXX34Z27Ztk77TSE3UqEdZVlYGd/fqr8PTarUoL79NNx0lIqJaU8EJ5yjreiX+549//CNWrFiB4cOH44MPPkBubq5Ty6/x5SEfffQRvL2rvqi6sFDugmYiIiJnWbVqFZYtW4ZPP/0Ua9aswfTp0zFs2DAIISrdQaQ2atRQtmzZEh9++KFiTENlDi2D8LQfEx6Up1jOJU9fqeVZJKdWR4VmScUFeRQpxvyfV6BUWWFe+VJx0U1PSsWpJU5mXDL5SZU10POsVFygZBaiL0uUt8lX1zpLlZVXrpeK8w4uU4w5YwiRKutIjlxuy9beV6XiRjT5QTHmqyK57fFLcahU3JjAE1JxrbXK63C6PFiqrLH+x6Xitl3rIRWXK7HvUwuUjzVTsUFqeQ5rZEnR9Xo9YmNjERsbizNnzmDNmjU4fvw47rvvPsTExODhhx/GmDFjalV2jRrK9PT0Wi2EiIjqmUaYmadCu3btEB8fj/nz52PXrl1YvXo1Hn30URgMtfsRwhR2RETUKLm5uWHkyJEYOXIksrOza11OjRrKd999VyruhRdeqFVliIjoNmkkPcrvvvsOvXv3VowLCgpCcXEx0tPT0aVLlxoto0YN5dtvv60Yo1Kp2FASEdVzzsisUx8y8zz++OOIjIzEk08+iREjRlQ52fTnn3/Gxx9/jLVr12Lx4sW3tqFMS0uz+/758+crZe0hIiK6VX7++Wf885//xBtvvIFJkyahffv2CAsLg4eHB3Jzc3H69GkUFxdjzJgxSExMRFRUVI2X4dRzlLm5uVi/fj3WrFnjzGKJiMjZGsnQq1arxfPPP4/nn38eJ06cwJEjR5Ceno7S0lJ0794df/nLXzBgwAD4+/vXehmczENE5IoaSUN5o7vuugt33XWX08t17q3iiYiIGhn2KImIXFBjmcxzO9SooVTKapCXl1ejhR8+fBhvvvkmkpOTkZWVhYSEBIwePbrK2KeffhqrVq3C22+/jenTpwMArl27hjlz5mDfvn3IyMhAYGAgRo8ejXnz5sHPTy7Ly418/Uqh9rSf7qiFd55iOcXl2hov255WXjlScYFa5cw8BUYPqbIi9XLLDNfkScXJ0KpMUnFmyS/nNYtc3uFCi/I2KTHJ7dNio1w2oKtG5cTN+Wa5LD85pQrppP7HUyOXMSq9XDl7TGppkFRZ54uaSsWd9ZIrz8etVDHmhxK57GB+mhKpuLOFctmsCg06xZgSib8N5hK11PIc1sgy89xKNWoolRofPz8/PP7449LlFRcXo3v37pg6dard+4Xt2LEDx44dQ1hYmM3rmZmZyMzMxFtvvYXOnTvj3LlzeOaZZ5CZmYnPPvtMuh5ERC6nEZ6jvFVq1FCuXbvWqQuPjo5GdHS03ZiLFy/i+eefx969exETE2PzXlRUFLZt22Z93qZNG8yfPx+PPfYYTCYTNBqOLBMRuaKysjJ4eMiNoCmp15N5LBYLJk+ejJdeekn6AtH8/Hz4+vqykSQissPhW2w54Ryns1ksFsybNw/NmzeHt7c3zp69fgOF2bNnY/Xq1bUut143lIsWLYJGo5HO9JOTk4N58+bh6aefthtnMBhQUFBg8yAicinCSY965B//+AfWrVuHxYsX29w7uWvXrvjoo49qXW69bSiTk5OxbNkyrFu3DiqV8gnjgoICxMTEoHPnzpgzZ47d2Pj4ePj5+Vkf4eHhzqo2ERHVkQ0bNmDVqlWYNGkS1OrfJ0V169YNp0+frnW59bahPHLkCLKzs9GyZUtoNBpoNBqcO3cOM2fORGRkpE1sYWEhhg8fDm9vbyQkJECrtT+zbNasWcjPz7c+MjIybuGaEBHVQ84Ydq1nPcqLFy+ibdvK92e1WCwwGo21LrfensibPHkyBg8ebPPasGHDMHnyZEydOtX6WkFBAYYNGwadToedO3dKnbzV6XTQ6ZSnchMRNVqNcNZrly5dcOTIEURERNi8vnXrVtx55521LrdOG8qioiKkpqZan6elpSElJQX+/v5o2bIlAgICbOK1Wi1CQkLQoUMHANd7kkOHDkVJSQk+/vhjm/ONzZo1s+l6ExFR4zZnzhxMnjwZFy9ehMViwfbt2/Hrr79iw4YN+OKLL2pdbp02lMePH8eAAQOsz2fMmAEAiI2Nxbp16xQ/n5ycjGPHjgFApe52WlpapSFaJSaLG4TF/mh0iUn5gvJyk9xm1arNUnEWyYt6yyzOS3SgVcnVTQv7CRoqWKC8Dh4quaGRL4rkZkC31l2WipNZ1+HNfpIqK9sod1H/g74pijF5FrlRD5PCMVuhnV5uewzz+lUxppmmUKqsZu5y5/9buMsluPBxK1OMud/7N6mycsyVb8dUlW5NLkrF5RmVEz8Um5X/fhiLy1H7s2k10Ah7lCNHjsSWLVuwYMECqFQqvPHGG7jrrrvw+eefY8iQIbUut04byv79+0MI+S2dnp7u0OeJiOi6xprCbtiwYRg2bJhTy6y3k3mIiIhqIikpyTrKeKNjx47h+PHjtS6XDSURETUK06ZNq/IqhosXL2LatGm1LrfeznolIqJbqBGeo/z555+rvB/lnXfeiZ9//rnW5bJHSUTkghpjCjudTofLlytPWsvKynIorSkbSiIiahSGDBliTShTIS8vD6+++mrDnfVKRER1qJ71CB21ZMkS9O3bFxEREdYEAykpKQgODsa//vWvWpfLhpKIyBU1wnOUzZs3xw8//ICNGzfi5MmT0Ov1mDp1Kh599FHF1Kb2sKEkIqJGw8vLC0899ZRTy2RDeYPifD3cyu3nij0ncfa66KqX1PJUWrmsNr/4hEjFubspZ5i5WOgnVZZJyJ2+DnPPlYorlsgyU2JRzloCACVmuYw1Ojd/qbh8s14x5h59mlRZIZo8qTgfiX0VplbOQgMAP3lekopr4y6XmedYmXI2nc66LKmyyoTcr/jE3CipuDCPPMWY9h5y2+PLnK5Sce28sqXiiiSOyyjPTMWYMpURO6SW6JjGmnDgt99+w8GDB5GdnQ2LxfZv7BtvvFGrMtlQEhG5okY49Prhhx/i2WefRWBgIEJCQmxu0ViR0q422FASEVGj8I9//APz58/H3/72N6eWy4aSiMgFNcah19zcXIwbN87p5fI6SiIiVySc9KhHxo0bh3379jm9XPYoiYioUWjbti1mz56No0ePomvXrpUuCXnhhRdqVS4bSiIiV9QIJ/OsWrUK3t7eOHToEA4dOmTznkqlYkNJRETyGuM5yrQ0ucu4aooNJRGRK2qEPcpbhQ0lERE1GhcuXMDOnTtx/vx5lJeX27y3dOnSWpXJhvJWkEu4A2GUm3RcapLLbmJRqxRjSgxyZeWWeUrFXTN5S8UVmu1nPAKAHKNcRqMRTX6QimunzZGKKxbKX4OE/Mr3uKtKgUl5PQGgSdMSxZhMY1Opsr7JaysVl+stt3391Mp121/cUaqss6XNpOIMFrVUXJFJOftNmaX2OT2rcrncVyrOJLEOPxaGKcYYi8sVY5yiEfYov/76a4waNQqtWrXCr7/+iqioKKSnp0MIUeV9KmXx8hAiIhfUGO9HOWvWLMycORM//vgjPDw8sG3bNmRkZKBfv34OXV/JhpKIiG6LlStXolu3bvD19YWvry969+6NL7/80vq+EAJxcXEICwuDXq9H//798dNPP0mX/8svvyA2NhYAoNFoUFpaCm9vb/z973/HokWLal1vNpRERK6oDhIOtGjRAgsXLsTx48dx/PhxDBw4EA8++KC1MVy8eDGWLl2K5cuXIykpCSEhIRgyZAgKCwulyvfy8oLBYAAAhIWF4f/+7/+s7129erVmlb0Bz1ESEbmgurg8ZOTIkTbP58+fj5UrV+Lo0aPo3Lkz3nnnHbz22msYM2YMAGD9+vUIDg7Gpk2b8PTTTyuW36tXL/znP/9B586dERMTg5kzZ+LUqVPYvn07evXqVbPK3oANJREROaSgoMDmuU6ng05nf+KV2WzG1q1bUVxcjN69eyMtLQ2XLl3C0KFDbcrp168fvv32W6mGcunSpSgqKgIAxMXFoaioCFu2bEHbtm3x9ttv12LNrmNDSUTkipw46zU83PYepnPmzEFcXFyVHzl16hR69+6NsrIyeHt7IyEhAZ07d8a3334LAAgODraJDw4Oxrlz56Sq07p1a+v/PT09sWLFCskVsY8NJRGRK3JiQ5mRkQFf398vo7HXm+zQoQNSUlKQl5eHbdu2ITY21ibd3I33kASuT/C5+bXqtG7dGklJSQgICLB5PS8vD3fddRfOnj0rVc7N2FASEZFDKmaxynB3d0fbttev/e3RoweSkpKwbNky6z0kL126hNDQUGt8dnZ2pV5mddLT02E2myu9bjAYcPHiRakyqsKGkojIBan+93C0DEcJIWAwGNCqVSuEhIQgMTERd955JwCgvLwchw4dUry0Y+fOndb/7927F35+ftbnZrMZX3/9NSIjI2tdRzaUN3DL0cLNw35WjzyhnN1Ee01ys0oOe1zwksvQotFW/iV1s7IcvVRZZyUz+Hzr3lo5CECZWbm83DK5umlVyusJAMU+qVJxIZo8xZgr5T5SZRkscvv+VFm4YozsespsWwAosbhLxZ0tDVSM6eB5WaosN8lpkekFAcpBAK7plL9/Gje51FgXi/2UgwDofY1Scbnlysevt9agGKNyM0ktz2F1kJnn1VdfRXR0NMLDw1FYWIjNmzfj4MGD2LNnD1QqFaZPn44FCxagXbt2aNeuHRYsWABPT09MnDjRbrmjR48GcH3YtuI6ygparRaRkZFYsmRJzSp7AzaUREQuqC4uD7l8+TImT56MrKws+Pn5oVu3btizZw+GDBkCAHj55ZdRWlqK5557Drm5uejZsyf27dsHHx/7P1Qtlus/jlq1aoWkpCQEBir/2KsJNpRERHRbrF692u77KpUKcXFx1c6YVVLVbbby8vLQpEmTWpVXgZl5iIhcUR1k5rnVFi1ahC1btlifjxs3Dv7+/mjevDlOnjxZ63LZUBIRuapG1EgCwD//+U/rNZ2JiYn46quvsGfPHkRHR+Oll16qdbkceiUiokYhKyvL2lB+8cUXeOSRRzB06FBERkaiZ8+etS6XPUoiIhfUGG+z1bRpU2RkZAAA9uzZg8GDBwO4fglKVddXymKPkojIFTXCGzePGTMGEydORLt27ZCTk4Po6GgAQEpKijXJQW2woSQiokbh7bffRmRkJDIyMrB48WJ4e3sDuD4k+9xzz9W6XDaUNzB7myH09rvn3k1KFcspMqillqcyyeW1CPAtkYrz1ilfzHzRLLdMP+8yqbg23nL3eDNZlEf5L2qbSJXVWn9FKk4mkQAA7Cvsqhhz/IpyggAAMJrlzmb4apSPo2tG5YvrAeD01SCpuHKz3Nf9vgDlRA0phS2kysoqkbuo/67ADKm4IK3yfQllEyvc10wu7+cP+c2l4gwm5e2bUaCcPMRcovw9doa6uI7yVtNqtfjrX/9a6fXp06c7VG6dnqM8fPgwRo4cibCwMKhUKuzYsaPa2KeffhoqlQrvvPOO9bX09HSoVKoqH1u3br31K0BE1FA1kstDdu7cCaPRaP2/vUdt1WmPsri4GN27d8fUqVMxduzYauN27NiBY8eOISwszOb18PBwZGVl2by2atUqLF682Do2TUREjdfo0aNx6dIlBAUFWVPZVUWlUtV6Qk+dNpTR0dGKDdrFixfx/PPPY+/evYiJibF5T61WIyQkxOa1hIQEjB8/3jo2TURElTWWodeK9HU3/9+Z6vU5SovFgsmTJ+Oll15Cly5dFOOTk5ORkpKC999/326cwWCAwfD7eYCb785NRNToNcJZr7dKvb6OctGiRdBoNHjhhRek4levXo1OnTqhT58+duPi4+Ph5+dnfdx8d24iokavkZyjrGCxWLBmzRo88MADiIqKQteuXTFq1Chs2LABQjhW0XrbUCYnJ2PZsmVYt26d1N2tS0tLsWnTJjzxxBOKsbNmzUJ+fr71UXGBKhERNTxCCIwaNQp/+tOfcPHiRXTt2hVdunTBuXPnMGXKFDz00EMOlV9vh16PHDmC7OxstGzZ0vqa2WzGzJkz8c477yA9Pd0m/rPPPkNJSQkef/xxxbJ1Oh10Op2zq0xE1GA0lnOUALBu3TocPnwYX3/9NQYMGGDz3v79+zF69Ghs2LBBqn2oSr3tUU6ePBk//PADUlJSrI+wsDC89NJL2Lt3b6X41atXY9SoUWjWrFkd1JaIqIFpREOvn3zyCV599dVKjSQADBw4EK+88go2btxY6/LrtEdZVFSE1NTfL25OS0tDSkoK/P390bJlSwQE2N71XKvVIiQkBB06dLB5PTU1FYcPH8bu3btvS72JiKj++OGHH7B48eJq34+Ojsa7775b6/LrtKE8fvy4zS+AGTNmAABiY2Oxbt066XLWrFmD5s2bY+jQoc6uYiVCSGS2kZyhrDLKZckxmOQy/ajdtIoxpjLlGAAoUsv9VMw22L/zeAWDRFaYzCK5LC5nPeRGDbQquWum/q9YubxrBZ5SZcm6WNZEMeZKqdwlTsVFHlJx1/Ry65BR5q8YI5txJ79Mrm75Rr1UnEXi+1dqlsvM46WRy4BTapL7zhjNyt9Tme+yWfL77iiVEFA5OMnF0c87y7Vr1xAcHFzt+8HBwcjNza11+XXaUPbv379Gs5FuPi9ZYcGCBViwYIGTakVE5AIa0eUhZrMZGk31zZlarYbJZKp1+fV2Mg8REZEMIQSmTJlS7STNG6+brw02lERELqgxzXqNjY1VjKntjFeADSURkWtqREOva9euvaXl19vLQ4iIiOoD9iiJiFxQYxp6vdXYUBIRuaJGNPR6q7GhJCJyQexRyuM5SiIiIjvYo7yBm94EN0/7F6Xq3Y2K5RRr5H5mySZb8JRYJgBo3CRSArnJLVOrlbs4V6+Wq5tMZh43yZ+nssv0cpO7dkomQ4tWW7s7o1dHZh1kt4fEzXVqVJ5FKP9+Npjk/nSYzHK/xQuNcjcpMFqUs9aUSRxrAFBqlsu4U2qUi5Mhs6skd6fjOPQqjQ0lEZGLcpWhU0dx6JWIiMgO9iiJiFyRENcfjpbhAthQEhG5IM56lcehVyIiIjvYoyQickWc9SqNDSURkQtSWa4/HC3DFXDolYiIyA72KImIXBGHXqWxobyBr18p1J72xxLCfXMVyykqlcsyYjbL5eBo3zRbKk7nppw9xizklhnukycVd5/vGam4ayZvxZgLnk2lyhrs85NUXLbZRypOJkuOWi03xmSWzEQjsx/cJP8KCcnhL4NZOasNIJexxmiRW8+ycrmsNn7aMqm4UI98xZiknAipsjz0RVJxeSV6qTgZ5QblP7mW0tsznslZr/LYUBIRuSJeRymN5yiJiIjsYI+SiMgFcehVHhtKIiJXxMk80jj0SkREZAd7lERELohDr/LYUBIRuSLOepXGoVciIiI72KO8QVGRB9wsHnZjLqibKJZTlmu/jJrKKJK7EF/G1Vy5i/DdJMdUfvEOk4q7Uq6ccOBKmXIMALTxCJGKu9cjTSouoInyheeFRrl9ahRyvz1H+J9SjEn1CpYqyyS5zK5NM6XiPN3KFWN0gSapsi55yR1vF0v8pOJyDF6KMXf4X5Aq63hOS6m4js0uS8XJfGfKZJI5FJdD7sh1DIde5bGhJCJyRZz1Ko1Dr0RERHawR0lE5II49CqPDSURkSuyiOsPR8twAWwoiYhcEc9RSuM5SiIiIjvYoyQickEqOOEcpVNqUv+xoSQickXMzCONQ69ERER21GmP8vDhw3jzzTeRnJyMrKwsJCQkYPTo0db3p0yZgvXr19t8pmfPnjh69CgA4Nq1a5gzZw727duHjIwMBAYGYvTo0Zg3bx78/OQyfdzIVK6Gm1ptN6a4zF2xHJVB7veH0Mj9GiszOW83CeHcwRKjsL+9Klgkssdo3CxSZbmrJLPCmH2l4sqEcraUYRKZdACgxKKTivN0MyjGRPv8IFWW7D7wU5dKxZ03+CvG9PJJlSrrgi5AKm53ZhepOKFVPn7zjHqpsgrK5LItBXgUS8UVGZX3faT3NcWYcotyZiRn4OUh8uq0R1lcXIzu3btj+fLl1cYMHz4cWVlZ1sfu3but72VmZiIzMxNvvfUWTp06hXXr1mHPnj144oknbkf1iYgaLuGkRw3Ex8fjnnvugY+PD4KCgjB69Gj8+uuvNjFTpkyBSqWyefTq1av26+kEddqjjI6ORnR0tN0YnU6HkJCqc3tGRUVh27Zt1udt2rTB/Pnz8dhjj8FkMkGj4SlYIqL64tChQ5g2bRruuecemEwmvPbaaxg6dCh+/vlneHn9nsd3+PDhWLt2rfW5u7vySN6tVO9bkoMHDyIoKAhNmjRBv379MH/+fAQFBVUbn5+fD19fXzaSRER2qISAysHJODX9/J49e2yer127FkFBQUhOTkbfvn2tr9vrINWFet2aREdHY9y4cYiIiEBaWhpmz56NgQMHIjk5GTpd5fMBOTk5mDdvHp5++mm75RoMBhgMv58jKigocHrdiYjqNcv/Ho6Wgcp/Q3U6XZV/o2+Wn58PAPD3tz0vXtMO0q1Wr2e9jh8/HjExMYiKisLIkSPx5Zdf4rfffsOuXbsqxRYUFCAmJgadO3fGnDlz7JYbHx8PPz8/6yM8PPxWrQIRUaMXHh5u8zc1Pj5e8TNCCMyYMQP3338/oqKirK9HR0dj48aN2L9/P5YsWYKkpCQMHDjQpnNzu9XrHuXNQkNDERERgTNnzti8XlhYiOHDh8Pb2xsJCQnQau3PYpw1axZmzJhhfV5QUMDGkohcijOHXjMyMuDr+/ssc5ne5PPPP48ffvgB33zzjc3r48ePt/4/KioKPXr0QEREBHbt2oUxY8Y4VN/aalANZU5ODjIyMhAaGmp9raCgAMOGDYNOp8POnTvh4aE85Vt2WICIqNFyYq5XX19fm4ZSyZ///Gfs3LkThw8fRosWLezGVtdBup3qtKEsKipCaurv12OlpaUhJSUF/v7+8Pf3R1xcHMaOHYvQ0FCkp6fj1VdfRWBgIB566CEA13uSQ4cORUlJCT7++GMUFBRYx8qbNWsGtcI1kURELqsOMvMIIfDnP/8ZCQkJOHjwIFq1aqX4mao6SLdbnTaUx48fx4ABA6zPK4ZDY2NjsXLlSpw6dQobNmxAXl4eQkNDMWDAAGzZsgU+Pj4AgOTkZBw7dgwA0LZtW5uy09LSEBkZeXtWhIiIFE2bNg2bNm3Cv//9b/j4+ODSpUsAAD8/P+j1ehQVFSl2kOpCnTaU/fv3h7Dzi2Tv3r0Ofb6m3DN0UHvYH5It81PO4uKZJTdHyiK59S9pm8oFSlDnKtcfALKMcuvwvTpCKq64XPk6KINRboM0dW8jFdfN+4JUXDv3S4oxPXW5UmUlGeRm5u0v6KwY01GfJVVWK90Vqbgdl++UirtW5qkYU2qWO44uljSRisvMlovT6pSzMhUZ5a65u3pZbqiwtFxuXQ1lynGXfbwVY8wlt2fSSl1k5lm5ciWA63+7b7R27VpMmTIFarVasYNUFxrUOUoiInKSOhp6tUev1yt2kOpCvb48hIiIqK6xR0lE5IJUlusPR8twBWwoiYhcEe9HKY1Dr0RERHawR0lE5IqcmHCgsWNDSUTkguri7iENFYdeiYiI7GCPkojIFXEyjzQ2lDfQFgNqhcQfKqHcCdflyh08Zg+VVJyhQG43yUzVds+XW2aZXm6Z2QXKmUYAwGBQzloiLHJ1k832onNTzuICAFqVWTHGw80oVdZ/SyKl4n4rUM7gUy6ZuslN8kTR2Vx/5SAApSXKNwwwWuQGo3ILlbP8AID6ktxNCox65eMos0Quk457plxcSalczmg3g/Lxe61UeZ9aSsuklucwAcfvR+ka7SQbSiIiV8RzlPJ4jpKIiMgO9iiJiFyRgBPOUTqlJvUeG0oiIlfEyTzSOPRKRERkB3uURESuyAJAbqK5/TJcABtKIiIXxFmv8jj0SkREZAd7lDfwvCygdrf/C8kocX2912XlC9gBoNxb7ndKuY9cnFTCgUKpomDWyy2zRK+XK7BcojzJYaDsYrkkBxrJm+XJJCYwWOQuTj9VECYVl5HXRDHGaJG70N0i5DZc4VUvqTgYlffVVfhIFWXKd5eK88yVWweTQbluJqPcvvK4JrdMN6PcftBI5AmQSbxhKZP7++EwTuaRxoaSiMgVsaGUxqFXIiIiO9ijJCJyRexRSmNDSUTkinh5iDQ2lERELoiXh8jjOUoiIiI72KMkInJFPEcpjQ0lEZErsghA5WBDZ3GNhpJDr0RERHawR3kDnwtl0ChsEaOn8ibzPJcvtTxTE7msNsLNUyrOzaT8605XIJf1w80od2iUlMtlXlEblGOE5M+2K2gqFXfNRy6Dz0U/P8UYf32JVFn/dzlQKk5keyjGpPrK7XfZmYf6dMl9Va4cY2gil63GO09uWqVvutxKlHsrl2f0kjuQ/NKVMzIBgMFPrjxNqfL3r6ypcllmie3vFBx6lcaGkojIJTmhoXSROzdz6JWIiMgO9iiJiFwRh16lsaEkInJFFgGHh04565WIiIjYoyQickXCcv3haBkugA0lEZEr4jlKaWwoiYhcEc9RSuM5SiIiIjvqtEd5+PBhvPnmm0hOTkZWVhYSEhIwevRo6/tTpkzB+vXrbT7Ts2dPHD16FACQnp6OVq1aVVn2p59+inHjxtWoPsJNBeFmP/OHm1niF5RJLvuN0rJqSiYzj+wPSCFZNZXkKQqVRBIUIZc4RprFLLcS5SblLDNFRrnKmcvkvlKacuW6WQxyv2NV5ZKZY+SSC0ll5rGo5bate57cAafLk8uSozIr7yu1Qa5uHpdLpeLcynVScZoS5e+9cFM+jkzG23Tej0Ov0uq0R1lcXIzu3btj+fLl1cYMHz4cWVlZ1sfu3but74WHh9u8l5WVhblz58LLywvR0dG3YxWIiBomgd8by1o/6nolbo867VFGR0crNmg6nQ4hISFVvqdWqyu9l5CQgPHjx8PbWy7PJxERkT31/hzlwYMHERQUhPbt2+PJJ59EdnZ2tbHJyclISUnBE088YbdMg8GAgoICmwcRkUtxuDfpjFyxDUO9biijo6OxceNG7N+/H0uWLEFSUhIGDhwIg6HqW1GsXr0anTp1Qp8+feyWGx8fDz8/P+sjPDz8VlSfiKj+slic83AB9bqhHD9+PGJiYhAVFYWRI0fiyy+/xG+//YZdu3ZVii0tLcWmTZsUe5MAMGvWLOTn51sfGRkZt6L6RETUCDSo6yhDQ0MRERGBM2fOVHrvs88+Q0lJCR5//HHFcnQ6HXQ6uZlsRESNEme9SmtQDWVOTg4yMjIQGhpa6b3Vq1dj1KhRaNasWR3UjIiogWFDKa1OG8qioiKkpqZan6elpSElJQX+/v7w9/dHXFwcxo4di9DQUKSnp+PVV19FYGAgHnroIZtyUlNTcfjwYZtLR4iIiJyhThvK48ePY8CAAdbnM2bMAADExsZi5cqVOHXqFDZs2IC8vDyEhoZiwIAB2LJlC3x8fGzKWbNmDZo3b46hQ4fe1voTETVYTGEnrU4byv79+0PY6brv3btXqpwFCxZgwYIFDtdH99slaJQyZ+gkMmukn5danntxmFScrzpYKs7NoJwZxK3UKFWWWecrFaeyyM0Hk8n2YvKQKgpGb+XsLABgLpWrW4FBIjOPh1zltJe1UnG6a8rZY8ol6gUAaoksPwDgeVluhqJMtiWVSW6Z+hy5ZeqyJdMGWfSKIWrZ7ZZbLBWnk8yUo84rUoxRmZooxphMZVLLc5QQFggH7/7h6OcbigZ1jpKIiJxECMd7hC5yjrJeXx5CRERU19ijJCJyRcIJ5yhdpEfJhpKIyBVZLPK3/6mOi5yj5NArERGRHexREhG5Ig69SmNDSUTkgoTFAuHg0KurXB7CoVciIiI72KO8gfDxglDbT5Zu8VZOpq4xVH2j6UplBfpJxZUFKCc5AACVxDVRbga5ZPClAXK/oUqbyV147qacCwFmudWEMUAuaYLaRy7O17Pq27bdyMdDOQYALkpe7G7RKCcmMDeVqz+MkkkfDJJfd4nRNENTuSE3s15ue0DloxwDoNRfeV3NHnLHJBAoFWXSy5Wny1VOSmHwV97vJqMK+F5qkY6pg6HX+Ph4bN++HadPn4Zer0efPn2waNEidOjQ4YYiBebOnYtVq1YhNzcXPXv2xPvvv48uXbo4VlcHsEdJROSKLMI5jxo4dOgQpk2bhqNHjyIxMREmkwlDhw5FcfHvWZIWL16MpUuXYvny5UhKSkJISAiGDBmCwsJCZ28BaexREhHRbbFnzx6b52vXrkVQUBCSk5PRt29fCCHwzjvv4LXXXsOYMWMAAOvXr0dwcDA2bdqEp59+ui6qzR4lEZFLEuL6dZAOPa73KAsKCmweBoPcqYr8/HwAgL+/P4Drd5C6dOmSzQ0udDod+vXrh2+//dbJG0AeG0oiIhckLMIpDwAIDw+Hn5+f9REfH6+8fCEwY8YM3H///YiKigIAXLp0CQAQHGx7I4jg4GDre3WBQ69ERK5IWAA4JzNPRkYGfH1/v+OQTqc8afD555/HDz/8gG+++abSeyqV7QQqIUSl124nNpREROQQX19fm4ZSyZ///Gfs3LkThw8fRosWLayvh4Rcv2Lg0qVLCA0Ntb6enZ1dqZd5O3HolYjIBTlz6FV6mULg+eefx/bt27F//360atXK5v1WrVohJCQEiYmJ1tfKy8tx6NAh9OnTxynrXRvsURIRuSInDr3KmjZtGjZt2oR///vf8PHxsZ539PPzg16vh0qlwvTp07FgwQK0a9cO7dq1w4IFC+Dp6YmJEyc6VlcHsKHE9V85AGAyK8/UskhcOC8s5VLLtUgsDwBMRrk7nkslHDBKXiheLndomA1y5w2ETMIByR+nllK5C/FVark4M5T3l8yxAQCWUsl9Vaa8QWTXUzbhgNmJCQcsZZLHkeTxYTJKHCAAzOUSCQckz2WZjCa5OI1ceWqTcnky62n+3/dd3OI8qiYYHc43YILkMfo/K1euBAD079/f5vW1a9diypQpAICXX34ZpaWleO6556wJB/bt2wcfH7mkFLeCStzqvdEAXLhwAeHh4XVdDSIiq4yMDJvzd85SVlaGVq1aOW0WaUhICNLS0uDhoZyZqKFiQwnAYrEgMzMTPj4+1plVBQUFCA8PrzSbq6Fo6PUHuA71QUOvP9Dw1kEIgcLCQoSFhcHN7dZMIykrK0N5udzIlxJ3d/dG3UgCHHoFALi5uVX7y62ms7nqm4Zef4DrUB809PoDDWsd/Pzk8kDXloeHR6Nv3JyJs16JiIjsYENJRERkBxvKauh0OsyZM0cqw0R91NDrD3Ad6oOGXn+gcawD1S1O5iEiIrKDPUoiIiI72FASERHZwYaSiIjIDjaUREREdrh0Q7ly5Up069bNeiFy79698eWXX1rfF0IgLi4OYWFh0Ov16N+/P3766ac6rHFlSuswZcoUqFQqm0evXr3qsMb2xcfHWxMjV2gI+6FCVfWv7/sgLi6uUv0qbncENIztr7QO9X0fUP3m0g1lixYtsHDhQhw/fhzHjx/HwIED8eCDD1r/CCxevBhLly7F8uXLkZSUhJCQEAwZMgSFhYV1XPPfKa0DAAwfPhxZWVnWx+7du+uwxtVLSkrCqlWr0K1bN5vXG8J+AKqvP1D/90GXLl1s6nfq1Cnrew1l+9tbB6D+7wOqxwTZaNq0qfjoo4+ExWIRISEhYuHChdb3ysrKhJ+fn/jggw/qsIbKKtZBCCFiY2PFgw8+WLcVklBYWCjatWsnEhMTRb9+/cSLL74ohBANZj9UV38h6v8+mDNnjujevXuV7zWU7W9vHYSo//uA6jeX7lHeyGw2Y/PmzSguLkbv3r2RlpaGS5cuYejQodYYnU6Hfv364dtvv63Dmlbv5nWocPDgQQQFBaF9+/Z48sknkZ2dXYe1rNq0adMQExODwYMH27zeUPZDdfWvUN/3wZkzZxAWFoZWrVphwoQJOHv2LICGs/2B6tehQn3fB1R/uXxS9FOnTqF3794oKyuDt7c3EhIS0LlzZ+sfgeDgYJv44OBgnDt3ri6qWq3q1gEAoqOjMW7cOERERCAtLQ2zZ8/GwIEDkZycXG8ylWzevBknTpxAUlJSpfcqbgVUn/eDvfoD9X8f9OzZExs2bED79u1x+fJl/OMf/0CfPn3w008/NYjtD9hfh4CAgHq/D6h+c/mGskOHDkhJSUFeXh62bduG2NhYHDp0yPq+6qabwAohKr1W16pbh86dO2P8+PHWuKioKPTo0QMRERHYtWsXxowZU4e1vi4jIwMvvvgi9u3bZ/duBvV1P8jUv77vg+joaOv/u3btit69e6NNmzZYv369dcJLfd3+Feytw4wZM+r9PqD6zeWHXt3d3dG2bVv06NED8fHx6N69O5YtW2adMXfzzU2zs7Mr/bqua9WtQ1VCQ0MRERGBM2fO3OZaVi05ORnZ2dm4++67odFooNFocOjQIbz77rvQaDTWbV1f94NS/c3myne0r2/74GZeXl7o2rUrzpw506C+Bze6cR2qUt/3AdUvLt9Q3kwIAYPBgFatWiEkJASJiYnW98rLy3Ho0CH06dOnDmuorGIdqpKTk4OMjAyEhobe5lpVbdCgQTh16hRSUlKsjx49emDSpElISUlB69at6/V+UKq/Wq2u9Jn6tg9uZjAY8MsvvyA0NLTBfg9uXIeq1Pd9QPVMnU4lqmOzZs0Shw8fFmlpaeKHH34Qr776qnBzcxP79u0TQgixcOFC4efnJ7Zv3y5OnTolHn30UREaGioKCgrquOa/s7cOhYWFYubMmeLbb78VaWlp4sCBA6J3796iefPm9WodbnbzrNGGsB9udGP9G8I+mDlzpjh48KA4e/asOHr0qHjggQeEj4+PSE9PF0I0jO1vbx0awj6g+s2lz1FevnwZkydPRlZWFvz8/NCtWzfs2bMHQ4YMAQC8/PLLKC0txXPPPYfc3Fz07NkT+/btg4+PTx3X/Hf21qG0tBSnTp3Chg0bkJeXh9DQUAwYMABbtmypV+ugpCHsh+qo1ep6vw8uXLiARx99FFevXkWzZs3Qq1cvHD16FBEREQAaxva3tw6N5XtAdYe32SIiIrKD5yiJiIjsYENJRERkBxtKIiIiO9hQEhER2cGGkoiIyA42lERERHawoSQiIrKDDSUREZEdbCiJqjBlyhSoVCosXLjQ5vUdO3ZY75px8OBBqFQqNG3aFGVlZTZx33//PVQqlc0dNiriKx7NmjVDdHQ0Tp48eetXiIhqjQ0lUTU8PDywaNEi5Obm2o3z8fFBQkKCzWtr1qxBy5Ytq4z/9ddfkZWVhV27diE3NxfDhw9Hfn6+0+pNRM7FhpKoGoMHD0ZISAji4+PtxsXGxmLNmjXW56Wlpdi8eTNiY2OrjA8KCkJISAjuvfdeLFmyBJcuXcLRo0edWncich42lETVUKvVWLBgAd577z1cuHCh2rjJkyfjyJEjOH/+PABg27ZtiIyMxF133aW4DL1eDwAwGo3OqTQROR0bSiI7HnroIdxxxx2YM2dOtTFBQUGIjo7GunXrAFwfdv3jH/+oWHZOTg7mzp0LHx8f3Hvvvc6qMhE5GRtKIgWLFi3C+vXr8fPPP1cb88c//hHr1q3D2bNn8d1332HSpEnVxrZo0QLe3t4IDAzEL7/8gq1btyIoKOhWVJ2InIANJZGCvn37YtiwYXj11VerjRkxYgTKysrwxBNPYOTIkQgICKg29siRIzh58iTy8/Px22+/YdiwYbei2kTkJC5942YiWQsXLsQdd9yB9u3bV/m+Wq3G5MmTsXjxYnz55Zd2y2rVqhWaNGlyC2pJRLcCe5REErp27YpJkybhvffeqzZm3rx5uHLlCnuIRI0MG0oiSfPmzYMQotr33d3dERgYaJNkgIgaPpWw980nIiJycexREhER2cGGkoiIyA42lERERHawoSQiIrKDDSUREZEdbCiJiIjsYENJRERkBxtKIiIiO9hQEhER2cGGkoiIyA42lERERHawoSQiIrLj/wEVL++HZEYPCgAAAABJRU5ErkJggg==",
"text/plain": [
"