{ "cells": [ { "cell_type": "markdown", "id": "sitting-budapest", "metadata": {}, "source": [ "# USP test of on continuous data\n", "\n", "We demonstrate how the [USP test of independence](https://doi.org/10.1098/rspa.2021.0549) can be applied to continuous data.\n", "\n", "A test of independence is stronger than a test for zero correlation. A test of independence can also detect dependencies which give zero correlation. " ] }, { "cell_type": "code", "execution_count": 1, "id": "embedded-warner", "metadata": {}, "outputs": [], "source": [ "from resample import permutation as perm\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "rng = np.random.default_rng(1)\n", "\n", "x1 = rng.normal(0, 2, size=100)\n", "y1 = rng.normal(0, 3, size=100)\n", "\n", "cov = np.empty((2, 2))\n", "cov[0, 0] = 2 ** 2\n", "cov[1, 1] = 3 ** 2\n", "rho = 0.5\n", "cov[0, 1] = rho * np.sqrt(cov[0, 0] * cov[1, 1])\n", "cov[1, 0] = cov[0, 1]\n", "\n", "xy2 = rng.multivariate_normal([0, 0], cov, size=500)\n", "\n", "d = {\"x,y are independent\": (x1, y1), \"x,y are correlated\": xy2.T}" ] }, { "cell_type": "code", "execution_count": 2, "id": "prompt-military", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1AAAAGHCAYAAAC+v7OPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAABqz0lEQVR4nO3dd1gU1/s28Hspu4D0uhARbMEGNtSgiWIkYI+xJcYoqMESsFcSo6BGjL0bNQrmG40lieVn7yWKXbATNBA0AmpUEAsInPcPXyaOgCzqsoven+uaS+fMmXOeMwMsDzNzRiGEECAiIiIiIqJiGeg6ACIiIiIiorKCCRQREREREZGGmEARERERERFpiAkUERERERGRhphAERERERERaYgJFBERERERkYaYQBEREREREWmICRQREREREZGGmEARERERERFpiAkUEdELuLu7IygoSO/b1KU3bTylJSkpCQqFAtHR0boOhYiISoAJFBER6b2tW7ciPDxc12G8dY4cOYL3338fZmZmUKvVGDRoEDIzMzXaNz09HaNGjULVqlVhamoKNzc39OnTB8nJyQXq/vPPP+jatSusra1haWmJjz/+GH/99VeBemlpaejVqxccHR1hamqKevXqYd26da88TiKikjDSdQBERPosPj4eBgb8W5Oubd26FQsWLGASVYpiY2PRokULVK9eHTNnzsT169cxffp0JCQkYNu2bS/cNy8vDx999BEuXryIr776Cu+++y6uXLmChQsXYseOHbh06RIsLCwAAJmZmWjevDnS09Px9ddfw9jYGLNmzUKzZs0QGxsLOzs7AEBGRgbef/99pKWlYfDgwVCr1Vi7di26du2KlStX4vPPP9f6MSEiAphAEdEbKi8vD9nZ2TAxMXmldlQq1WuKiKhs+frrr2FjY4P9+/fD0tISwNPbNYODg7Fz5074+/sXue/Ro0dx4sQJzJ8/HyEhIVK5h4cHevfujd27d+OTTz4BACxcuBAJCQk4fvw4GjRoAABo1aoVatWqhRkzZmDy5MkAgMWLF+PKlSvYs2cPPvzwQwDAgAED8N5772H48OHo3LkzlEqlVo4FEdGz+GdVInqhR48eoVq1aqhWrRoePXokld+5cwfOzs5o3LgxcnNzC903MzMT5cqVw+DBgwtsu379OgwNDREZGfnC/qdPn47GjRvDzs4OpqamqF+/Pn799dcC9RQKBUJDQ7Fy5UrUrFkTKpUK27dvB/D09qDevXvDyckJKpUKNWvWxPLlyzUa//PP90RHR0OhUODw4cMYNmwYHBwcUK5cOXzyySe4deuWbF8hBCZNmoTy5cvDzMwMzZs3x4ULFwrt5969exgyZAhcXV2hUqlQpUoVfP/998jLy5Pq5D8zM336dMyaNQtubm4wNTVFs2bNcP78+QJtXr58GZ07d4atrS1MTEzg7e2NTZs2yeroy3iWLFmCypUrQ6VSoUGDBjhx4oRULygoCAsWLADw9DznL/lSUlJw+fJlPHnypNBYXvb4PSstLQ1GRkaIiIgosC0+Ph4KhQLz588H8PR7Y8SIEfD09IS5uTksLS3RqlUrxMXFvbAPAPD19YWvr2+B8qCgILi7u8vK8vLyMHv2bNSsWRMmJiZwcnJCv379cPfu3WL7KU5GRgZ27dqFL774QkqeAKBnz54wNzfH2rVri90fAJycnGTlzs7OAABTU1Op7Ndff0WDBg2k5AkAqlWrhhYtWsj6OXToEBwcHKTkCQAMDAzQtWtXpKam4sCBAy8xUiKikuMVKCJ6IVNTU6xYsQJNmjTBN998g5kzZwIAQkJCkJ6ejujoaBgaGha6r7m5OT755BOsWbMGM2fOlNX75ZdfIIRA9+7dX9j/nDlz0L59e3Tv3h3Z2dlYvXo1unTpgs2bN6NNmzayunv37sXatWsRGhoKe3t7uLu7Iy0tDe+9956UYDk4OGDbtm3o06cPMjIyMGTIkJc6LgMHDoSNjQ3Gjx+PpKQkzJ49G6GhoVizZo1UZ9y4cZg0aRJat26N1q1b4/Tp0/D390d2drasrYcPH6JZs2b4559/0K9fP1SoUAFHjhxBWFgYUlJSMHv2bFn9n376Cffv30dISAgeP36MOXPm4MMPP8S5c+ekX1gvXLiAJk2a4J133sGYMWNQrlw5rF27Fh06dMBvv/0m/fVfH8azatUq3L9/H/369YNCocDUqVPRsWNH/PXXXzA2Nka/fv1w48YN7Nq1C//73/8KnIuwsDCsWLECiYmJBZKMwmhy/J7n5OSEZs2aYe3atRg/frxs25o1a2BoaIguXboAAP766y9s2LABXbp0QcWKFZGWlobFixejWbNmuHjxIlxcXIqNURP9+vVDdHQ0evXqhUGDBiExMRHz58/HmTNncPjwYRgbGwN4+oeMx48fF9uesbExrKysAADnzp1DTk4OvL29ZXWUSiXq1KmDM2fOvLAtb29vlCtXDt9++y1sbW3h4eGBK1euYNSoUWjQoAH8/PwAPE0Cz549i969exdoo2HDhti5cyfu378PCwsLZGVlyRKvfGZmZgCAU6dO4aOPPip2nEREr0wQEWkgLCxMGBgYiIMHD4p169YJAGL27NnF7rdjxw4BQGzbtk1W7uXlJZo1a1bs/g8fPpStZ2dni1q1aokPP/xQVg5AGBgYiAsXLsjK+/TpI5ydncXt27dl5Z999pmwsrIq0P7z3NzcRGBgoLQeFRUlAAg/Pz+Rl5cnlQ8dOlQYGhqKe/fuCSGEuHnzplAqlaJNmzayel9//bUAIGtz4sSJoly5cuLPP/+U9T1mzBhhaGgokpOThRBCJCYmCgDC1NRUXL9+Xap37NgxAUAMHTpUKmvRooXw9PQUjx8/lsry8vJE48aNRdWqVfVqPHZ2duLOnTtSvY0bNwoA4v/+7/+kspCQEFHUR1ZgYKAAIBITEwvdnq8kx68wixcvFgDEuXPnZOU1atSQfT0+fvxY5ObmFuhbpVKJCRMmFIgnKipKKmvWrFmh3xeBgYHCzc1NWj906JAAIFauXCmrt3379gLl+cenuOXZfvO/xw8ePFggli5dugi1Wl3oMXrW5s2bhbOzs6yPgIAAcf/+fanOrVu3BADZccm3YMECAUBcvnxZCCHEwIEDhYGBgUhKSpLV++yzzwQAERoaWmxMRESvA2/hIyKNhIeHo2bNmggMDMRXX32FZs2aYdCgQcXu5+fnBxcXF6xcuVIqO3/+PM6ePYsvvvii2P2f/Yvz3bt3kZ6ejg8++ACnT58uULdZs2aoUaOGtC6EwG+//YZ27dpBCIHbt29LS0BAANLT0wttRxN9+/aV3Ub2wQcfIDc3F3///TcAYPfu3cjOzsbAgQNl9Qq74rVu3Tp88MEHsLGxkcXo5+eH3NxcHDx4UFa/Q4cOeOedd6T1hg0bolGjRti6dSuAp7eQ7d27F127dsX9+/el9v79918EBAQgISEB//zzj96M59NPP4WNjY2sbwCFzsJWmOjoaAghNLr6BBR//IrSsWNHGBkZya7KnT9/HhcvXsSnn34qlalUKmnikdzcXPz7778wNzeHh4fHS3+9PW/dunWwsrLCRx99JDvG9evXh7m5Ofbt2yfVHTVqFHbt2lXsMmPGDGmf/Nt1C3sG0MTERHY7b1EcHBxQt25dfPfdd9iwYQPCw8Nx6NAh9OrVS+N+nq3z5ZdfwtDQEF27dsWRI0dw9epVREZGYv369bJ6RETaxlv4iEgjSqUSy5cvR4MGDWBiYoKoqCjZL9JFMTAwQPfu3bFo0SI8fPgQZmZmWLlyJUxMTKRbnl5k8+bNmDRpEmJjY5GVlSWVF9Z3xYoVZeu3bt3CvXv3sGTJEixZsqTQ9m/evFlsDIWpUKGCbD0/Ach//iQ/8ahataqsnoODgyxZAICEhAScPXsWDg4OGsX4fJsA8O6770rPi1y5cgVCCHz77bf49ttvi2zz2SRCl+Mpru/Xrbjj9+jRI6Snp8u2q9Vq2NvbS8/lTJw4EcDT2/eMjIzQsWNHqW5eXh7mzJmDhQsXIjExUfaMYP6Mcq8qISEB6enpcHR0LHT7s8e4Ro0asj8saCL/DxfPfs/le/z4caG30j3rr7/+QvPmzfHTTz+hU6dOAICPP/5YeqZw27ZtaNWqVbH9PBuLl5cXVq1ahf79+6NJkyYAnp6X2bNnY8CAATA3Ny/RGImIXhYTKCLS2I4dOwA8/cUmISGhQMJSlJ49e2LatGnYsGEDunXrhlWrVqFt27bS8xZFOXToENq3b4+mTZti4cKFcHZ2hrGxMaKiorBq1aoC9Z//pS5/woIvvvgCgYGBhfbh5eWl0RieV9RzX0KIEreVP+XzqFGjCt3+7rvvlrg9ABgxYgQCAgIKrVOlShXZui7H8zr7fh3WrFkju0rybCyfffYZevXqhdjYWNSpUwdr165FixYtYG9vL9WdPHkyvv32W/Tu3RsTJ06Era0tDAwMMGTIENkkGoVRKBSFjvv5iVry8vLg6Ogou7L7rGeT1/T0dI2uziiVStja2gL4b7KHlJSUAvVSUlKKfY4rOjoajx8/Rtu2bWXl7du3BwAcPnwYrVq1gq2tLVQqVZH9AJD11blzZ7Rv3x5xcXHIzc1FvXr1sH//fgAl/z4hInpZTKCISCNnz57FhAkTpF8ev/zyS5w7d67YJAgAatWqhbp162LlypUoX748kpOTMW/evGL3++2332BiYoIdO3bIbvGJiorSKGYHBwdYWFggNzdXemi9tLi5uQF4eqWgUqVKUvmtW7cKXFmpXLkyMjMzNY4xISGhQNmff/4p3cKW35+xsfFrG7c2x6MJTa52aqq44xcQEIBdu3YVum+HDh3Qr18/6Ta+P//8E2FhYbI6v/76K5o3b45ly5bJyu/duydLtApjY2NT6K2L+VcA81WuXBm7d+9GkyZNir0aNHjwYKxYseKFdYCnt8DmJyO1atWCkZERTp48ia5du0p1srOzERsbKysrTFpaGoQQBRK//JkSc3JyADy9Qu3p6YmTJ08WaOPYsWOoVKmS9L6ofEqlUjZj3+7duwGg1L/HiejtxWegiKhYT548QVBQEFxcXDBnzhxER0cjLS0NQ4cO1biNHj16YOfOnZg9ezbs7OzQqlWrYvcxNDSEQqGQ/RKWlJSEDRs2aNSnoaEhOnXqhN9++63Qaaqfn6b7dfLz84OxsTHmzZsnu6Lw/Ax0ANC1a1fExMRIV/iede/ePemXzXwbNmyQPcN0/PhxHDt2TDqmjo6O8PX1xeLFiwv9y/7LjFub49FEuXLlpP2fp+k05vmKO37Ozs7w8/OTLfmsra0REBCAtWvXYvXq1VAqlejQoYOsfUNDwwJXkdatW1fgubPCVK5cGZcvX5ado7i4OBw+fFhWr2vXrsjNzZVuJXxWTk6O7Di9zDNQVlZW8PPzw88//4z79+9L5f/73/+QmZkpu/324cOHuHz5Mm7fvi2VvfvuuxBCFJju/JdffgEA1K1bVyrr3LkzTpw4IUui4uPjsXfv3mJv801ISMAPP/yAtm3b8goUEZUeXcxcQURly7hx44RCoRB79+6VyiZNmiQAiC1btkhl+bOKPTsjW77U1FRhZGQkAIgBAwZo1O+ePXsEAPHBBx+IRYsWiYiICOHo6Ci8vLwKzMgGQISEhBTar5ubmzAzMxODBw8WixcvFpGRkaJLly7Cxsam2BiKmoXvxIkTsnr79u0TAMS+ffuksrCwMAFAtG7dWsyfP1/06dNHuLi4CHt7e1mbDx48EPXq1RNGRkbiyy+/FIsWLRLTp08XgYGBoly5cuLWrVtCiP+Or6enp3B3dxfff/+9mDBhgrC1tRV2dnbixo0bUpsXLlwQNjY2ws7OTowZM0YsWbJETJw4UbRu3Vp4eXnp1XimTZtW4LgDEOPHj5fW165dKwCIHj16iJ9//ln88ssv0raSzsKnyfF7kZ9//lkAEBYWFqJdu3YFto8bN04AEEFBQWLJkiVi4MCBwtbWVlSqVEk2011hs/BdvHhRGBgYiLp164r58+eLcePGCUdHR+Hp6SmbhU8IIfr16ycAiFatWolZs2aJ+fPni8GDBwsXFxexbt06jcbyIqdOnRIqlUrUrVtXLFq0SHzzzTfCxMRE+Pv7y+rlf608e75u374t1Gq1UCqVYtCgQWLx4sWiX79+wtDQUNSsWVNkZWVJdTMyMkTlypWFo6OjmDp1qpg1a5ZwdXUVLi4u4ubNm7K+qlevLsaNGyd+/PFH8c033whbW1vh5uYmm1WRiEjbmEAR0QudOnVKGBkZiYEDB8rKc3JyRIMGDYSLi4u4e/euEEKIc+fOCQBizJgxhbbVunVrAUAcOXJE4/6XLVsmqlatKlQqlahWrZqIiooS48eP1ziBEkKItLQ0ERISIlxdXYWxsbFQq9WiRYsWYsmSJcX2/yoJVG5uroiIiBDOzs7C1NRU+Pr6ivPnzxdoUwgh7t+/L8LCwkSVKlWEUqkU9vb2onHjxmL69OkiOztbCCFPOGbMmCFcXV2FSqUSH3zwgYiLiysQ+9WrV0XPnj2FWq0WxsbG4p133hFt27YVv/76q96N53nP/0Kek5MjBg4cKBwcHIRCoZCd/5ImUJoev6JkZGQIU1NTAUD8/PPPBbY/fvxYDB8+XDpOTZo0ETExMQWmKC8sgRLiaYJWqVIloVQqRZ06dcSOHTsKTGOeb8mSJaJ+/frC1NRUWFhYCE9PTzFq1CiNk8HiHDp0SDRu3FiYmJgIBwcHERISIjIyMmR1CkughBDi+vXronfv3qJixYpCqVQKZ2dnERwcLCXQz7p27Zro3LmzsLS0FObm5qJt27YiISGhQL3PPvtMuLq6CqVSKVxcXET//v1FWlraaxkrEZGmFELo6CldInrjLFy4EKNGjcLVq1cLfSHpJ598gnPnzuHKlSs6iK7sS0pKQsWKFTFt2jSMGDFC1+GUOTx+RET0OvAZKCJ6bfbt24dBgwYVmjylpKRgy5Yt6NGjhw4iIyIiIno9OAsfEb0269atK1CWmJiIw4cP48cff4SxsTH69eung8iIiIiIXg9egSIirTpw4AB69OiBxMRErFixAmq1WtchEREREb00PgNFRERERESkIV6BIiIiIiIi0hATKCIiIiIiIg0xgSIiIiIiItIQEygiIiIiIiINMYEiIiIiIiLSEBMoIiIiIiIiDTGBIiIiIiIi0hATKCIiIiIiIg0xgSIiIiIiItIQEygiIiIiIiINMYEiIiIiIiLSEBMoIiIiIiIiDTGBIiIiIiIi0hATKCIiIiIiIg0xgSIiIiIiItIQEygiIiIiIiINMYEiIiIiIiLSEBMoIiIiIiIiDTGBIiIiIiIi0pCRrgN4kby8PNy4cQMWFhZQKBS6DoeI6K0hhMD9+/fh4uICAwP+re1Z/GwiItINffls0usE6saNG3B1ddV1GEREb61r166hfPnyug5Dr/CziYhIt3T92aTXCZSFhQUA4P16I2BkqNJxNABOnNd1BPqrQS1dR/AfnieiV5aDJ/gDW6Wfw/Sf/GNy7do1WFpa6jiasu3BgwdwcXEB8DQxLVeunI4jIiJ9lpGRAVdXV51/Nul1ApV/a4SRoQpGRiY6jgaAwljXEegvfTg/+XieiF6dePoPb1ErKP+YWFpaMoF6RYaGhtL/LS0tmUARkUZ0/dnEG9uJiIiIiIg0xASKiIiIiIhIQ0ygiIiIiIiINKTXz0ARERERkWby8vKQnZ2t6zCIXpqxsbHs2Uh9xQSKiIiIqIzLzs5GYmIi8vLydB0K0SuxtraGWq3W+UQRL6LVBCo3Nxfh4eH4+eefkZqaChcXFwQFBWHs2LF6fVCIiIiIygohBFJSUmBoaAhXV1e+/JrKJCEEHj58iJs3bwIAnJ2ddRxR0bSaQH3//fdYtGgRVqxYgZo1a+LkyZPo1asXrKysMGjQIG12TURERPRWyMnJwcOHD+Hi4gIzMzNdh0P00kxNTQEAN2/ehKOjo97ezqfVBOrIkSP4+OOP0aZNGwCAu7s7fvnlFxw/flyb3RIRERG9NXJzcwEASqVSx5EQvbr8PwI8efJEbxMorV7jbdy4Mfbs2YM///wTABAXF4c//vgDrVq1KrR+VlYWMjIyZAsRERERFY+PR9CboCx8HWv1CtSYMWOQkZGBatWqwdDQELm5ufjuu+/QvXv3QutHRkYiIiJCmyERERERERG9NK1egVq7di1WrlyJVatW4fTp01ixYgWmT5+OFStWFFo/LCwM6enp0nLt2jVthkdERERERFQiWk2gRo4ciTFjxuCzzz6Dp6cnevTogaFDhyIyMrLQ+iqVCpaWlrKFiIjeDgcPHkS7du3g4uIChUKBDRs2yLYrFIpCl2nTpkl13N3dC2yfMmWKrJ2zZ8/igw8+gImJCVxdXTF16tTSGB4RvaUK+3lWEuHh4ahTp85ri+d1SUpKgkKhQGxsrK5DKXVavYXv4cOHBabSNDQ05DsKiIiogAcPHqB27dro3bs3OnbsWGB7SkqKbH3btm3o06cPOnXqJCufMGECgoODpXULCwvp/xkZGfD394efnx9++OEHnDt3Dr1794a1tTX69u37mkdE2uI+Zssr7Z80pc1rioReha+vL+rUqYPZs2e/tjaDgoJw7969V0pYXraN8PBwbNiwoUBCkZKSAhsbG43aUCgUWL9+PTp06CCVjRgxAgMHDixRLKRdWk2g2rVrh++++w4VKlRAzZo1cebMGcycORO9e/fWZrdERFQGtWrVqshJhgBArVbL1jdu3IjmzZujUqVKsnILC4sCdfOtXLkS2dnZWL58OZRKJWrWrInY2FjMnDmzyAQqKysLWVlZ0jonOCKikijq55GmzM3NYW5u/pqi0S/Z2dllcvZIrd7CN2/ePHTu3BlfffUVqlevjhEjRqBfv36YOHGiNrslIqI3XFpaGrZs2YI+ffoU2DZlyhTY2dmhbt26mDZtGnJycqRtMTExaNq0qewDOyAgAPHx8bh7926hfUVGRsLKykpaXF1dX/+AiN4yQUFBOHDgAObMmSPdbpuUlAQAOH/+PFq1agVzc3M4OTmhR48euH37trTvr7/+Ck9PT5iamsLOzg5+fn548OABwsPDsWLFCmzcuFFqc//+/YX2/zJtjB49Gu+++y7MzMxQqVIlfPvtt3jy5AkAIDo6GhEREYiLi5P2i46OBiC/hS87OxuhoaFwdnaGiYkJ3NzcpEdb3N3dAQCffPIJFAqFtF7YLXzLly9HzZo1oVKp4OzsjNDQ0CKPta+vL4YMGSIr69ChA4KCgqR1d3d3TJ48Gb1794aFhQUqVKiAJUuWyPY5fvw46tatCxMTE3h7e+PMmTMF+iru3Pn6+iI0NBRDhgyBvb09AgICioxbn2k1gbKwsMDs2bPx999/49GjR7h69SomTZpUJjNNIiLSHytWrICFhUWBW/0GDRqE1atXY9++fejXrx8mT56MUaNGSdtTU1Ph5OQk2yd/PTU1tdC+OMER0es3Z84c+Pj4IDg4GCkpKUhJSYGrqyvu3buHDz/8EHXr1sXJkyexfft2pKWloWvXrgCe3g7XrVs39O7dG5cuXcL+/fvRsWNHCCEwYsQIdO3aFS1btpTabNy4cYG+X7YNCwsLREdH4+LFi5gzZw6WLl2KWbNmAQA+/fRTDB8+HDVr1pT2+/TTTwv0PXfuXGzatAlr165FfHw8Vq5cKSVKJ06cAABERUUhJSVFWn/eokWLEBISgr59++LcuXPYtGkTqlSp8srnZMaMGVJi9NVXX2HAgAGIj48HAGRmZqJt27aoUaMGTp06hfDwcIwYMUK2f3HnLt+KFSugVCpx+PBh/PDDD68cty5o9RY+IiIibVi+fDm6d+8OExMTWfmwYcOk/3t5eUGpVKJfv36IjIyESqV6qb5UKtVL70tEhbOysoJSqYSZmZnsFrf58+ejbt26mDx5slS2fPlyuLq64s8//0RmZiZycnLQsWNHuLm5AQA8PT2luqampsjKynrhbXMpKSkv1cbYsWOl/7u7u2PEiBFYvXo1Ro0aBVNTU5ibm8PIyOiFfScnJ6Nq1ap4//33oVAopP4BwMHBAQBgbW39wjYmTZqE4cOHY/DgwVJZgwYNiqyvqdatW+Orr74C8PRq26xZs7Bv3z54eHhg1apVyMvLw7Jly2BiYoKaNWvi+vXrGDBggLR/cefu3XffBQBUrVq1zE/eo9UrUERERK/boUOHEB8fjy+//LLYuo0aNUJOTo50a5BarUZaWpqsTv76qz6nQESvLi4uDvv27ZOe+zE3N0e1atUAAFevXkXt2rXRokULeHp6okuXLli6dGmRt98W5WXbWLNmDZo0aQK1Wg1zc3OMHTsWycnJJeo7KCgIsbGx8PDwwKBBg7Bz584S7X/z5k3cuHEDLVq0KNF+mvDy8pL+r1AooFarcfPmTQDApUuX4OXlJfujlY+Pj2z/4s5dvvr167/22EsbEygiIipTli1bhvr166N27drF1o2NjYWBgQEcHR0BPP3AP3jwoPTcAgDs2rULHh4eGs+SRUTak5mZiXbt2iE2Nla2JCQkoGnTpjA0NMSuXbuwbds21KhRA/PmzYOHhwcSExM17uNl2oiJiUH37t3RunVrbN68GWfOnME333yD7OzsEo2vXr16SExMxMSJE/Ho0SN07doVnTt31nh/U1PTEvUHAAYGBhBCyMqe/RmYz9jYWLauUChKNHN2cecuX7ly5Uo4Av3DBIqIiPRCZmam9IELAImJiYiNjZX9hTcjIwPr1q0r9OpTTEwMZs+ejbi4OPz1119YuXIlhg4dii+++EJKjj7//HMolUr06dMHFy5cwJo1azBnzhzZrX9EVDqUSiVyc3NlZfXq1cOFCxfg7u6OKlWqyJb8X7wVCgWaNGmCiIgInDlzBkqlEuvXry+yzcKUtI0jR47Azc0N33zzDby9vVG1alX8/fffxY6nMJaWlvj000+xdOlSrFmzBr/99hvu3LkD4GkS86I2LCws4O7ujj179hTbTz4HBwfZayByc3Nx/vx5jfcHgOrVq+Ps2bN4/PixVHb06FFZHU3O3ZuCCRQREemFkydPom7duqhbty6Ap88z1a1bF+PGjZPqrF69GkIIdOvWrcD+KpUKq1evRrNmzVCzZk189913GDp0qGwmKSsrK+zcuROJiYmoX78+hg8fjnHjxvEdUEQ64O7ujmPHjiEpKQm3b99GXl4eQkJCcOfOHXTr1g0nTpzA1atXsWPHDvTq1Qu5ubk4duwYJk+ejJMnTyI5ORm///47bt26herVq0ttnj17FvHx8bh9+3ahV1pepo2qVasiOTkZq1evxtWrVzF37lwp4Xp2PPl/+Ll9+7bs9Qf5Zs6ciV9++QWXL1/Gn3/+iXXr1kGtVsPa2lpqY8+ePUhNTS3ytsLw8HDMmDEDc+fORUJCAk6fPo158+YVeZw//PBDbNmyBVu2bMHly5cxYMAA3Lt3T5NTJPn888+hUCgQHByMixcvYuvWrZg+fbqsTnHn7k3CSSSIiEgv+Pr6FrjN5Hl9+/YtMtmpV69egb+IFsbLywuHDh16qRiJypJXfeFwSZX0BcUjRoxAYGAgatSogUePHiExMRHu7u44fPgwRo8eDX9/f2RlZcHNzQ0tW7aEgYEBLC0tcfDgQcyePRsZGRlwc3PDjBkzpHfIBQcHY//+/fD29kZmZib27dsHX19fWb8v00b79u0xdOhQhIaGIisrC23atMG3336L8PBwqd1OnTrh999/R/PmzXHv3j1ERUXJpgoHnl5Bmjp1KhISEmBoaIgGDRpg69atMDB4ek1jxowZGDZsGJYuXYp33nlHen7zWYGBgXj8+DFmzZqFESNGwN7e/oW3Afbu3RtxcXHo2bMnjIyMMHToUDRv3rxE58rc3Bz/93//h/79+6Nu3bqoUaMGvv/+e9mLzF1cXF547t4kClHcp5UOZWRkwMrKCr4NvoGRkUnxO2jb0ThdR6C/3iv+WYRSw/NE9MpyxBPsx0akp6fD0tJS1+HolfzPJh6bV/fgwQPpBaGZmZka3+bzqolBSX/R13ePHz9GYmIiKlasKHvIX98TKKLCFPX1DOjPz983Kx0kIiIiIiLSIiZQREREREREGmICRUREREREpCEmUERERERERBpiAkVEREREZdL+/fuhUCikabmjo6OlKcELk5SUBIVCIb1vTt88Px5tcHd3x+zZs7XW/tugbExjfuI8oDAuvh7pDme+Kxs4WyIREZFe8PX1RZ06dUo9mTlx4sQb92Lb0lY2EigiIiIiInplDg4Oug6hzGMCRURERPQG0vf3Mvn6+sLLywsmJib48ccfoVQq0b9/f+nltElJSahYsSLOnDmDOnXqAADu3bsHGxubQl+QWxJ//fUXhg4dimPHjqFq1ar44Ycf4OPjgwcPHsDZ2RnLly+XvZx2w4YN6N69O1JTU/Hvv/+iYsWK+OWXXzB37lycPn0aVapUwYIFC9CsWTNpnwMHDmDkyJGIi4uDra0tAgMDMWnSJBgZGSEoKAgHDhzAgQMHMGfOHABAYmKitO+pU6cwevRoXLx4EXXq1EFUVBQ8PDyk7Rs3bkRERAQuXrwIFxcXBAYG4ptvvoGRkRGEEIiIiMDy5cuRlpYGOzs7dO7cGXPnzgXw9Ba+IUOGYMiQIcXWpcLxGSgiIiIi0okVK1agXLlyOHbsGKZOnYoJEyZg165dWu/3m2++wYgRIxAbG4t3330X3bp1Q05ODsqVK4fPPvsMUVFRsvpRUVHo3LkzLCwspLKRI0di+PDhOHPmDHx8fNCuXTv8+++/AIB//vkHrVu3RoMGDRAXF4dFixZh2bJlmDRpEgBgzpw58PHxQXBwMFJSUpCSkgJXV1dZfDNmzMDJkydhZGSE3r17S9sOHTqEnj17YvDgwbh48SIWL16M6OhofPfddwCA3377DbNmzcLixYuRkJCADRs2wNPTs9DjUJK69B9egSIiIqK3ivuYLa+0v75f2SlLvLy8MH78eABA1apVMX/+fOzZswcfffSRVvsdMWIE2rR5eh4jIiJQs2ZNXLlyBdWqVcOXX36Jxo0bIyUlBc7Ozrh58ya2bt2K3bt3y9oIDQ1Fp06dAACLFi3C9u3bsWzZMowaNQoLFy6Eq6sr5s+fD4VCgWrVquHGjRsYPXo0xo0bBysrKyiVSpiZmUGtVheI77vvvpOuZo0ZMwZt2rTB48ePYWJigoiICIwZMwaBgYEAgEqVKmHixIkYNWoUxo8fj+TkZKjVavj5+cHY2BgVKlRAw4YNCz0OJalL/+EVKCIiIiLSCS8vL9l6fsJSmv06OzsDgNRvw4YNUbNmTaxYsQIA8PPPP8PNzQ1NmzaVteHj4yP938jICN7e3rh06RIA4NKlS/Dx8YFCoZDqNGnSBJmZmbh+/forxRcXF4cJEybA3NxcWvKvZD18+BBdunTBo0ePUKlSJQQHB2P9+vXIyckptJ+S1KX/MIEiIiIiIp0wNpbPsqxQKJCXlwcAMDB4+muqEELa/uTJk9feb36Sk98vAHz55ZeIjo4G8PT2vV69esmSIW17UXyZmZmIiIhAbGystJw7dw4JCQkwMTGBq6sr4uPjsXDhQpiamuKrr75C06ZNCz12JalL/2ECRURERER6J3+2uJSUFKmstN7f9MUXX+Dvv//G3LlzcfHiRel2uWcdPXpU+n9OTg5OnTqF6tWrAwCqV6+OmJgYWfJ3+PBhWFhYoHz58gAApVKJ3NzcEsdWr149xMfHo0qVKgWW/KTT1NQU7dq1w9y5c7F//37ExMTg3LlzhbZXkrr0FJ+BIiIiIiK9Y2pqivfeew9TpkxBxYoVcfPmTYwdO7ZU+raxsUHHjh0xcuRI+Pv7S0nPsxYsWICqVauievXqmDVrFu7evStN9vDVV19h9uzZGDhwIEJDQxEfH4/x48dj2LBhUpLj7u6OY8eOISkpCebm5rC1tdUotnHjxqFt27aoUKECOnfuDAMDA8TFxeH8+fOYNGkSoqOjkZubi0aNGsHMzAw///wzTE1N4ebmVqCtktSl/zCBIiIiolKVP4lDXvZjqaz6t9thoDTRVUikp5YvX44+ffqgfv368PDwwNSpU+Hv718qfffp0werVq2SzYD3rClTpmDKlCmIjY1FlSpVsGnTJtjb2wMA3nnnHWzduhUjR45E7dq1YWtriz59+sgSwBEjRiAwMBA1atTAo0ePZNOYv0hAQAA2b96MCRMm4Pvvv4exsbE0+QUAWFtbY8qUKRg2bBhyc3Ph6emJ//u//4OdnV2BtkpSl/6jEM9eW9QzGRkZsLKygi8+hpHCuPgdiOjF3qut6wj+czRO1xHQC+SIJ9iPjUhPT4elpaWuw9Er+Z9NPDYv79kE6tqsp+/acR36a5lJoPRtFr7Hjx8jMTERFStWhIlJ2TiGZcH//vc/DB06FDdu3IBSqZTKC3s/Fb0+L/p61pefv7wCRURERET0/z18+BApKSmYMmUK+vXrJ0ueiABOIkFEREREJJk6dSqqVasGtVqNsLAwXYdDeohXoIiIiIiI/r/w8HCEh4cXud3d3R16/AQMlQKtX4H6559/8MUXX8DOzg6mpqbw9PTEyZMntd0tERERERHRa6fVK1B3795FkyZN0Lx5c2zbtg0ODg5ISEiAjY2NNrslIiIiIiLSCq0mUN9//z1cXV0RFRUllVWsWLHI+llZWcjKypLWMzIytBkeERERERFRiWj1Fr5NmzbB29sbXbp0gaOjI+rWrYulS5cWWT8yMhJWVlbS4urqqs3wiIiIiIiISkSrCdRff/2FRYsWoWrVqtixYwcGDBiAQYMGYcWKFYXWDwsLQ3p6urRcu3ZNm+ERERERERGViFZv4cvLy4O3tzcmT54MAKhbty7Onz+PH374AYGBgQXqq1QqqFQqbYZERERERET00rR6BcrZ2Rk1atSQlVWvXh3Jycna7JaIiIiISC8oFAps2LDhpfcPDw9HnTp1Xls89Oq0mkA1adIE8fHxsrI///wTbm5u2uyWiIiIiPScr68vhgwZ8lrbDAoKQocOHXTSRlGJTkpKClq1aqVRG4UlWyNGjMCePXtKHA9pj1Zv4Rs6dCgaN26MyZMno2vXrjh+/DiWLFmCJUuWaLNbIiIiIiK9oFarX2l/c3NzmJubv6Zo6HXQ6hWoBg0aYP369fjll19Qq1YtTJw4EbNnz0b37t212S0RERER6bGgoCAcOHAAc+bMgUKhgEKhQFJSEgDg/PnzaNWqFczNzeHk5IQePXrg9u3b0r6//vorPD09YWpqCjs7O/j5+eHBgwcIDw/HihUrsHHjRqnN/fv3F9r/y7QxevRovPvuuzAzM0OlSpXw7bff4smTJwCA6OhoREREIC4uTtovOjoagPyqUnZ2NkJDQ+Hs7AwTExO4ubkhMjISAODu7g4A+OSTT6BQKKT1wq5sLV++HDVr1oRKpYKzszNCQ0Nf6XxQyWg1gQKAtm3b4ty5c3j8+DEuXbqE4OBgbXdJRERl0MGDB9GuXTu4uLgUehtLUFCQ9ItJ/tKyZUtZnTt37qB79+6wtLSEtbU1+vTpg8zMTFmds2fP4oMPPoCJiQlcXV0xdepUbQ+NiJ4zZ84c+Pj4IDg4GCkpKUhJSYGrqyvu3buHDz/8EHXr1sXJkyexfft2pKWloWvXrgCe3g7XrVs39O7dG5cuXcL+/fvRsWNHCCEwYsQIdO3aFS1btpTabNy4cYG+X7YNCwsLREdH4+LFi5gzZw6WLl2KWbNmAQA+/fRTDB8+HDVr1pT2+/TTTwv0PXfuXGzatAlr165FfHw8Vq5cKSVKJ06cAABERUUhJSVFWn/eokWLEBISgr59++LcuXPYtGkTqlSp8srnhDSn1Vv4iIiINPXgwQPUrl0bvXv3RseOHQut07JlS9nL2Z+fubV79+5ISUnBrl278OTJE/Tq1Qt9+/bFqlWrADx9Qbu/vz/8/Pzwww8/4Ny5c+jduzesra3Rt29f7Q2OiGSsrKygVCphZmYmu8Vt/vz5qFu3rjSDM/D0aourqyv+/PNPZGZmIicnBx07dpSeqff09JTqmpqaIisr64W3zaWkpLxUG2PHjpX+7+7ujhEjRmD16tUYNWoUTE1NYW5uDiMjoxf2nZycjKpVq+L999+HQqGQzQvg4OAAALC2tn5hG5MmTcLw4cMxePBgqaxBgwZF1qfXjwkUERHphVatWhX7oLVKpSryF4tLly5h+/btOHHiBLy9vQEA8+bNQ+vWrTF9+nS4uLhg5cqVyM7OxvLly6FUKlGzZk3ExsZi5syZTKCI9EBcXBz27dtX6DM/V69ehb+/P1q0aAFPT08EBATA398fnTt3ho2NjcZ91K5d+6XaWLNmDebOnYurV69KiZylpWWJxhcUFISPPvoIHh4eaNmyJdq2bQt/f3+N97958yZu3LiBFi1alKhfer20fgsfERHR67J//344OjrCw8MDAwYMwL///itti4mJgbW1tZQ8AYCfnx8MDAxw7NgxqU7Tpk2hVCqlOgEBAYiPj8fdu3cL7TMrKwsZGRmyhYi0IzMzE+3atUNsbKxsSUhIQNOmTWFoaIhdu3Zh27ZtqFGjBubNmwcPDw8kJiZq3MfLtBETE4Pu3bujdevW2Lx5M86cOYNvvvkG2dnZJRpfvXr1kJiYiIkTJ+LRo0fo2rUrOnfurPH+pqamJeqPtIMJFBERlQktW7bETz/9hD179uD777/HgQMH0KpVK+Tm5gIAUlNT4ejoKNvHyMgItra2SE1Nleo4OTnJ6uSv59d5XmRkJKysrKTF1dX1dQ+N6K2kVCql79989erVw4ULF+Du7o4qVarIlnLlygF4OilDkyZNEBERgTNnzkCpVGL9+vVFtlmYkrZx5MgRuLm54ZtvvoG3tzeqVq2Kv//+u9jxFMbS0hKffvopli5dijVr1uC3337DnTt3AADGxsYvbMPCwgLu7u6c1lzHeAsfERGVCZ999pn0f09PT3h5eaFy5crYv3+/Vm9nCQsLw7Bhw6T1jIwMJlFEr4G7uzuOHTuGpKQkmJubw9bWFiEhIVi6dCm6deuGUaNGwdbWFleuXMHq1avx448/4uTJk9izZw/8/f3h6OiIY8eO4datW6hevbrU5o4dOxAfHw87OztYWVnB2NhY1u+xY8dK3EbVqlWRnJyM1atXo0GDBtiyZYuUcD07nsTERMTGxqJ8+fKwsLAo8JzmzJkz4ezsjLp168LAwADr1q2DWq2GtbW11MaePXvQpEkTqFSqQm8rDA8PR//+/eHo6IhWrVrh/v37OHz4MAYOHPi6Tg0Vg1egiIioTKpUqRLs7e1x5coVAE/ftXLz5k1ZnZycHNy5c0d6bkqtViMtLU1WJ3+9qGerVCoVLC0tZQsRvboRI0bA0NAQNWrUgIODA5KTk+Hi4oLDhw8jNzcX/v7+8PT0xJAhQ2BtbQ0DAwNYWlri4MGDaN26Nd59912MHTsWM2bMkJ6fDA4OhoeHB7y9veHg4IDDhw8X6Pdl2mjfvj2GDh2K0NBQ1KlTB0eOHMG3334ra7dTp05o2bIlmjdvDgcHB/zyyy8F+rawsMDUqVPh7e2NBg0aICkpCVu3boWBwdNfyWfMmIFdu3bB1dUVdevWLfS4BQYGYvbs2Vi4cCFq1qyJtm3bIiEh4ZXOBZWMQgghdB1EUTIyMmBlZQVffAwjhXHxO5DuvFdb1xFI/q1dTtchSOziHug6BBkemyIcjdN1BHonRzzBfmxEenq6ThIGhUKB9evXo0OHDkXWuX79OipUqIANGzagffv2uHTpEmrUqIGTJ0+ifv36AICdO3eiZcuWuH79OlxcXLBo0SJ88803SEtLk/4q/fXXX+P333/H5cuXNYot/7NJV8fmTeA+ZgsAIC/7Ma7Nevr8h+vQX2GgNNFlWBpLmtJG1yHIPH78GImJiahYsSJMTMrGMSQqyou+nvXl5y+vQBERkV7IzMyUHhgHIN0Kk5ycjMzMTIwcORJHjx5FUlIS9uzZg48//hhVqlRBQEAAAKB69epo2bIlgoODcfz4cRw+fBihoaH47LPP4OLiAgD4/PPPoVQq0adPH1y4cAFr1qzBnDlzZLfoERERvQifgSIiIr1w8uRJNG/eXFrPT2oCAwOxaNEinD17FitWrMC9e/fg4uICf39/TJw4UfaMwcqVKxEaGooWLVrAwMAAnTp1wty5c6XtVlZW2LlzJ0JCQlC/fn3Y29tj3LhxnMKc3ihCCDx8+FAnfZuZmUGhUOikb6LSwgSKiIj0gq+vL150V/mOHTuKbcPW1lZ6aW5RvLy8cOjQoRLHR1RWPHz4sND3KJWGzMxMabY8ojcVb+EjIiIiIiLSEK9AEREREb1BzMzMkJmZqbO+NSWEQL9+/fDrr7/i7t27OHPmDOrUqaO94DTg7u6OIUOGYMiQITqNg/QbEygiIiKiN4hCoSgTt9Ft374d0dHR2L9/v/RaAl07ceKETo5dUFAQ7t27hw0bNpR631RyTKCIiIiIqNRdvXoVzs7OaNy4cZF1srOzoVQqSy0mBweHUuuLyi4+A0VEREREpSooKAgDBw5EcnIyFAoF3N3dATydTCY0NBRDhgyBvb299JqCAwcOoGHDhlCpVHB2dsaYMWOQk5Mjtefr64uBAwdiyJAhsLGxgZOTE5YuXYoHDx6gV69esLCwQJUqVbBt27YXxuXu7o7Zs2dL6wqFAj/++CM++eQTmJmZoWrVqti0aZO0ff/+/VAoFNiyZQu8vLxgYmKC9957D+fPn5fqhIeHF7g1cfbs2dKYw8PDsWLFCmzcuBEKhQIKhQL79+8v+UGlUsMEioiIiIhK1Zw5czBhwgSUL18eKSkpOHHihLRtxYoVUCqVOHz4MH744Qf8888/aN26NRo0aIC4uDgsWrQIy5Ytw6RJk2RtrlixAvb29jh+/DgGDhyIAQMGoEuXLmjcuDFOnz4Nf39/9OjRo8RTvEdERKBr1644e/YsWrduje7du+POnTuyOiNHjsSMGTNw4sQJODg4oF27dnjy5IlG7Y8YMQJdu3ZFy5YtkZKSgpSUlBdelSPdYwJFRERERKXKysoKFhYWMDQ0hFqtlt06V7VqVUydOhUeHh7w8PDAwoUL4erqivnz56NatWro0KEDIiIiMGPGDOTl5Un71a5dG2PHjkXVqlURFhYGExMT2NvbIzg4GFWrVsW4cePw77//4uzZsyWKNSgoCN26dUOVKlUwefJkZGZm4vjx47I648ePx0cffQRPT0+sWLECaWlpWL9+vUbtm5ubw9TUFCqVCmq1Gmq1ulRvW6SSYwJFRERERHqjfv36svVLly7Bx8dH9oLeJk2aIDMzE9evX5fKvLy8pP8bGhrCzs4Onp6eUpmTkxMA4ObNmyWK59l2y5UrB0tLywJt+Pj4SP+3tbWFh4cHLl26VKJ+qOxgAkVEREREeuNlZ8EzNjaWrSsUCllZfgL27FWrl223JG0YGBgUeEm4prf3kX5iAkVEREREeqt69eqIiYmRJSGHDx+GhYUFypcvr8PI/nP06FHp/3fv3sWff/6J6tWrA3g6s19qaqos/tjYWNn+SqUSubm5pRIrvTomUERERESkt7766itcu3YNAwcOxOXLl7Fx40aMHz8ew4YNg4GBfvwqO2HCBOzZswfnz59HUFAQ7O3t0aFDBwBPZwi8desWpk6diqtXr2LBggUFZgN0d3fH2bNnER8fj9u3b/MKlZ7Tj686IiIiIqJCvPPOO9i6dSuOHz+O2rVro3///ujTpw/Gjh2r69AkU6ZMweDBg1G/fn2kpqbi//7v/6SJIKpXr46FCxdiwYIFqF27No4fP44RI0bI9g8ODoaHhwe8vb3h4OCAw4cP62IYpCG+SJeIiIiISt2QIUMwZMgQWVlR7z9q1qxZgZnvitsvKSmpQNnzzyIVt09h9e/du1eg7P3335e9++l5/fv3R//+/WVlX3/9tfR/BwcH7Ny584Wxkf7gFSgiIiIiIiINMYEiIiIiIiLSUKklUFOmTIFCoShwqZaIiIiIqCzy9fWFEALW1ta6DoVKUakkUCdOnMDixYtlLyIjIiIiIiIqa7SeQGVmZqJ79+5YunQpbGxstN0dERER0VupuAkSiMqCsvB1rPUEKiQkBG3atIGfn1+xdbOyspCRkSFbiIiIiKhohoaGAIDs7GwdR0L06h4+fAgAMDY21nEkRdPqNOarV6/G6dOnceLECY3qR0ZGIiIiQpshEREREb1RjIyMYGZmhlu3bsHY2FhvXi5LVBJCCDx8+BA3b96EtbW19IcBfaS1BOratWsYPHgwdu3aBRMTE432CQsLw7Bhw6T1jIwMuLq6aitEIiIiojJPoVDA2dkZiYmJ+Pvvv3UdDtErsba2hlqt1nUYL6S1BOrUqVO4efMm6tWrJ5Xl5ubi4MGDmD9/PrKysgpkliqVCiqVSlshEREREb2RlEolqlatytv4qEwzNjbW6ytP+bSWQLVo0QLnzp2TlfXq1QvVqlXD6NGjy8TBISIiIiorDAwMNL7rh4hentYSKAsLC9SqVUtWVq5cOdjZ2RUoJyIiIiIiKgv4lCEREREREZGGtDoL3/P2799fmt0RERERERG9VrwCRUREREREpKFSvQJFr9l7tXUdwX+Oxuk6Aokd9Oe4tF2+X9chyKyY00bXIUiufmaq6xAklY/qOgIiIiIqK3gFioiIiIiISEO8AkVERERUAu5jtrzS/klT9OduACIqOV6BIiIiIiIi0hATKCIi0gsHDx5Eu3bt4OLiAoVCgQ0bNkjbnjx5gtGjR8PT0xPlypWDi4sLevbsiRs3bsjacHd3h0KhkC1TpkyR1Tl79iw++OADmJiYwNXVFVOnTi2N4RER0RuCCRQREemFBw8eoHbt2liwYEGBbQ8fPsTp06fx7bff4vTp0/j9998RHx+P9u3bF6g7YcIEpKSkSMvAgQOlbRkZGfD394ebmxtOnTqFadOmITw8HEuWLNHq2IiI6M3BZ6CIiEgvtGrVCq1atSp0m5WVFXbt2iUrmz9/Pho2bIjk5GRUqFBBKrewsIBarS60nZUrVyI7OxvLly+HUqlEzZo1ERsbi5kzZ6Jv376vbzBERPTG4hUoIiIqk9LT06FQKGBtbS0rnzJlCuzs7FC3bl1MmzYNOTk50raYmBg0bdoUSqVSKgsICEB8fDzu3r1baD9ZWVnIyMiQLURE9PbiFSgiIipzHj9+jNGjR6Nbt26wtLSUygcNGoR69erB1tYWR44cQVhYGFJSUjBz5kwAQGpqKipWrChry8nJSdpmY2NToK/IyEhERERocTRERFSWMIEiIqIy5cmTJ+jatSuEEFi0aJFs27Bhw6T/e3l5QalUol+/foiMjIRKpXqp/sLCwmTtZmRkwNXV9eWCJyKiMo8JFBERlRn5ydPff/+NvXv3yq4+FaZRo0bIyclBUlISPDw8oFarkZaWJquTv17Uc1Mqleqlky8iInrz8BkoIiIqE/KTp4SEBOzevRt2dnbF7hMbGwsDAwM4OjoCAHx8fHDw4EE8efJEqrNr1y54eHgUevseERHR83gFioiI9EJmZiauXLkirScmJiI2Nha2trZwdnZG586dcfr0aWzevBm5ublITU0FANja2kKpVCImJgbHjh1D8+bNYWFhgZiYGAwdOhRffPGFlBx9/vnniIiIQJ8+fTB69GicP38ec+bMwaxZs3QyZiIiKnuYQBERkV44efIkmjdvLq3nP3cUGBiI8PBwbNq0CQBQp04d2X779u2Dr68vVCoVVq9ejfDwcGRlZaFixYoYOnSo7PklKysr7Ny5EyEhIahfvz7s7e0xbtw4TmFOREQaYwJFRER6wdfXF0KIIre/aBsA1KtXD0ePHi22Hy8vLxw6dKjE8REREQF8BoqIiIiIiEhjTKCIiIiIiIg0xASKiIiIiIhIQ0ygiIiIiIiINMQEioiIiIiISENMoIiIiIiIiDTEBIqIiIiIiEhDTKCIiIiIiIg0xASKiIiIiIhIQ0ygiIiIiIiINKTVBCoyMhINGjSAhYUFHB0d0aFDB8THx2uzSyIiIiIiIq3RagJ14MABhISE4OjRo9i1axeePHkCf39/PHjwQJvdEhERERERaYWRNhvfvn27bD06OhqOjo44deoUmjZtqs2uiYiIiIiIXjutJlDPS09PBwDY2toWuj0rKwtZWVnSekZGRqnERUREREREpIlSm0QiLy8PQ4YMQZMmTVCrVq1C60RGRsLKykpaXF1dSys8IiIiIiKiYpXaFaiQkBCcP38ef/zxR5F1wsLCMGzYMGk9IyODSVQZ8W+/xroOQXKvep6uQ5AMtE7WdQhyg7foOgLJyqSGug5BP71XW9cRPJXzGDixUddREBER6Z1SSaBCQ0OxefNmHDx4EOXLly+ynkqlgkqlKo2QiIiIiIiISkyrCZQQAgMHDsT69euxf/9+VKxYUZvdERERERERaZVWE6iQkBCsWrUKGzduhIWFBVJTUwEAVlZWMDU11WbXREREREREr51WJ5FYtGgR0tPT4evrC2dnZ2lZs2aNNrslIiIiIiLSCq3fwkdERERERPSmKLVpzImIiIiIiMo6JlBEREREREQaKrX3QBEREdGbwX2M/rxTjoiotPEKFBERERERkYaYQBEREREREWmICRQREREREZGGmEARERERERFpiAkUERERERGRhphAERERERERaYgJFBER6YWDBw+iXbt2cHFxgUKhwIYNG2TbhRAYN24cnJ2dYWpqCj8/PyQkJMjq3LlzB927d4elpSWsra3Rp08fZGZmyuqcPXsWH3zwAUxMTODq6oqpU6dqe2hERPQGYQJFRER64cGDB6hduzYWLFhQ6PapU6di7ty5+OGHH3Ds2DGUK1cOAQEBePz4sVSne/fuuHDhAnbt2oXNmzfj4MGD6Nu3r7Q9IyMD/v7+cHNzw6lTpzBt2jSEh4djyZIlWh8fERG9GfgiXSIi0gutWrVCq1atCt0mhMDs2bMxduxYfPzxxwCAn376CU5OTtiwYQM+++wzXLp0Cdu3b8eJEyfg7e0NAJg3bx5at26N6dOnw8XFBStXrkR2djaWL18OpVKJmjVrIjY2FjNnzpQlWkREREXhFSgiItJ7iYmJSE1NhZ+fn1RmZWWFRo0aISYmBgAQExMDa2trKXkCAD8/PxgYGODYsWNSnaZNm0KpVEp1AgICEB8fj7t37xbad1ZWFjIyMmQLERG9vZhAERGR3ktNTQUAODk5ycqdnJykbampqXB0dJRtNzIygq2traxOYW0828fzIiMjYWVlJS2urq6vPiAiIiqzmEARERG9QFhYGNLT06Xl2rVrug6JiIh0iAkUERHpPbVaDQBIS0uTlaelpUnb1Go1bt68Kduek5ODO3fuyOoU1sazfTxPpVLB0tJSthAR0duLCRQREem9ihUrQq1WY8+ePVJZRkYGjh07Bh8fHwCAj48P7t27h1OnTkl19u7di7y8PDRq1Eiqc/DgQTx58kSqs2vXLnh4eMDGxqaURkNERGUZEygiItILmZmZiI2NRWxsLICnE0fExsYiOTkZCoUCQ4YMwaRJk7Bp0yacO3cOPXv2hIuLCzp06AAAqF69Olq2bIng4GAcP34chw8fRmhoKD777DO4uLgAAD7//HMolUr06dMHFy5cwJo1azBnzhwMGzZMR6MmIqKyhtOYExGRXjh58iSaN28urecnNYGBgYiOjsaoUaPw4MED9O3bF/fu3cP777+P7du3w8TERNpn5cqVCA0NRYsWLWBgYIBOnTph7ty50nYrKyvs3LkTISEhqF+/Puzt7TFu3DhOYU5ERBpjAkVERHrB19cXQogitysUCkyYMAETJkwoso6trS1WrVr1wn68vLxw6NChl46TiIjebryFj4iIiIiISENMoIiIiIiIiDTEBIqIiIiIiEhDTKCIiIiIiIg0xASKiIiIiIhIQ5yFryw7GqfrCP5Tu7GuI5AM9t+m6xAk78V21nUIMkfr/KrrECRz/myl6xAkVu+Z6ToEIiIiKiNK5QrUggUL4O7uDhMTEzRq1AjHjx8vjW6JiIiIiIheK60nUGvWrMGwYcMwfvx4nD59GrVr10ZAQABu3ryp7a6JiIiIiIheK60nUDNnzkRwcDB69eqFGjVq4IcffoCZmRmWL1+u7a6JiIiIiIheK60mUNnZ2Th16hT8/Pz+69DAAH5+foiJiSlQPysrCxkZGbKFiIiIiIhIX2g1gbp9+zZyc3Ph5OQkK3dyckJqamqB+pGRkbCyspIWV1dXbYZHRERERERUIno1jXlYWBjS09Ol5dq1a7oOiYiIiIiISKLVaczt7e1haGiItLQ0WXlaWhrUanWB+iqVCiqVSpshERERERERvTStXoFSKpWoX78+9uzZI5Xl5eVhz5498PHx0WbXREREREREr53WX6Q7bNgwBAYGwtvbGw0bNsTs2bPx4MED9OrVS9tdExEREekd9zFbXmn/pCltXlMkRPQytJ5Affrpp7h16xbGjRuH1NRU1KlTB9u3by8wsQQREREREZG+03oCBQChoaEIDQ0tja6IiIiIiIi0Rq9m4SMiIiIiItJnTKCIiIiIiIg0xASKiIiIiIhIQ0ygiIiIiIiINMQEioiIiIiISENMoIiIiIiIiDTEBIqIiIiIiEhDTKCIiIiIiIg0xASKiIiIiIhIQ0ygiIiIiIiINMQEioiIiIiISENMoIiIiIiIiDTEBIqIiIiIiEhDTKCIiIiIiIg0xASKiIjKBHd3dygUigJLSEgIAMDX17fAtv79+8vaSE5ORps2bWBmZgZHR0eMHDkSOTk5uhgOERGVUUa6DoCIiEgTJ06cQG5urrR+/vx5fPTRR+jSpYtUFhwcjAkTJkjrZmZm0v9zc3PRpk0bqNVqHDlyBCkpKejZsyeMjY0xefLk0hkEERGVeUygiIioTHBwcJCtT5kyBZUrV0azZs2kMjMzM6jV6kL337lzJy5evIjdu3fDyckJderUwcSJEzF69GiEh4dDqVRqNX4iInoz8BY+IiIqc7Kzs/Hzzz+jd+/eUCgUUvnKlSthb2+PWrVqISwsDA8fPpS2xcTEwNPTE05OTlJZQEAAMjIycOHChSL7ysrKQkZGhmwhIqK3F69AERFRmbNhwwbcu3cPQUFBUtnnn38ONzc3uLi44OzZsxg9ejTi4+Px+++/AwBSU1NlyRMAaT01NbXIviIjIxEREfH6B0FERGUSE6iSeK+2riOQufqZqa5DeEaergOQrExqqOsQJDkbHYqvVIreQ2ddhyCpPOSorkOgFxFPdB3BCy1btgytWrWCi4uLVNa3b1/p/56ennB2dkaLFi1w9epVVK5c+aX7CgsLw7Bhw6T1jIwMuLq6vnR7RERUtjGBIiKiMuXvv//G7t27pStLRWnUqBEA4MqVK6hcuTLUajWOHz8uq5OWlgYART43BQAqlQoqleoVoyYiojcFn4EiIqIyJSoqCo6OjmjTps0L68XGxgIAnJ2dAQA+Pj44d+4cbt68KdXZtWsXLC0tUaNGDa3FS0REbxZegSIiojIjLy8PUVFRCAwMhJHRfx9hV69exapVq9C6dWvY2dnh7NmzGDp0KJo2bQovLy8AgL+/P2rUqIEePXpg6tSpSE1NxdixYxESEsIrTEREpDEmUEREVGbs3r0bycnJ6N27t6xcqVRi9+7dmD17Nh48eABXV1d06tQJY8eOleoYGhpi8+bNGDBgAHx8fFCuXDkEBgbK3htFRERUHCZQRERUZvj7+0MIUaDc1dUVBw4cKHZ/Nzc3bN26VRuhERHRW4LPQBEREREREWlIawlUUlIS+vTpg4oVK8LU1BSVK1fG+PHjkZ2dra0uiYiIiIiItEprt/BdvnwZeXl5WLx4MapUqYLz588jODgYDx48wPTp07XVLRERERERkdZoLYFq2bIlWrZsKa1XqlQJ8fHxWLRoUZEJVFZWFrKysqT1jIwMbYVHRERERERUYqX6DFR6ejpsbW2L3B4ZGQkrKytp4ZveiYiIiIhIn5RaAnXlyhXMmzcP/fr1K7JOWFgY0tPTpeXatWulFR4REREREVGxSpxAjRkzBgqF4oXL5cuXZfv8888/aNmyJbp06YLg4OAi21apVLC0tJQtRERERERE+qLEz0ANHz4cQUFBL6xTqVIl6f83btxA8+bN0bhxYyxZsqTEARIREREREemLEidQDg4OcHBw0KjuP//8g+bNm6N+/fqIioqCgQFfO0VERERERGWX1mbh++eff+Dr6ws3NzdMnz4dt27dkrap1WptdUtERERERKQ1Wkugdu3ahStXruDKlSsoX768bJsQQlvdEhERERERaY3W7qkLCgqCEKLQhYiIiIiIqCziQ0lEREREREQaYgJFRERERESkISZQREREREREGmICRUREREREpCEmUERERERERBpiAkVERERERKQhJlBEREREREQaYgJFRERERESkISZQREREREREGmICRUREREREpCEmUERERERERBpiAkVERERERKQhJlBEREREREQaMtJ1AGXK0ThdRyBTGbV1HYJ+Opqg6wgkV2fb6ToEuT/1KJ6tug7gP1at9edrhoiIiPQbr0ARERERERFpiFegiIiIiMoQ9zFbXmn/pCltXlMkRG8nXoEiIiIiIiLSEBMoIiIiIiIiDTGBIiKiMiE8PBwKhUK2VKtWTdr++PFjhISEwM7ODubm5ujUqRPS0tJkbSQnJ6NNmzYwMzODo6MjRo4ciZycnNIeChERlWF8BoqIiMqMmjVrYvfu3dK6kdF/H2NDhw7Fli1bsG7dOlhZWSE0NBQdO3bE4cOHAQC5ublo06YN1Go1jhw5gpSUFPTs2RPGxsaYPHlyqY+FiIjKJiZQRERUZhgZGUGtVhcoT09Px7Jly7Bq1Sp8+OGHAICoqChUr14dR48exXvvvYedO3fi4sWL2L17N5ycnFCnTh1MnDgRo0ePRnh4OJRKZWkPh4iIyiDewkdERGVGQkICXFxcUKlSJXTv3h3JyckAgFOnTuHJkyfw8/OT6larVg0VKlRATEwMACAmJgaenp5wcnKS6gQEBCAjIwMXLlwoss+srCxkZGTIFiIiensxgSIiojKhUaNGiI6Oxvbt27Fo0SIkJibigw8+wP3795GamgqlUglra2vZPk5OTkhNTQUApKamypKn/O3524oSGRkJKysraXF1dX29AyMiojKFt/AREVGZ0KpVK+n/Xl5eaNSoEdzc3LB27VqYmppqrd+wsDAMGzZMWs/IyGASRUT0FuMVKCIiKpOsra3x7rvv4sqVK1Cr1cjOzsa9e/dkddLS0qRnptRqdYFZ+fLXC3uuKp9KpYKlpaVsISKit1epJFBZWVmoU6cOFAoFYmNjS6NLIiJ6w2VmZuLq1atwdnZG/fr1YWxsjD179kjb4+PjkZycDB8fHwCAj48Pzp07h5s3b0p1du3aBUtLS9SoUaPU4yciorKpVBKoUaNGwcXFpTS6IiKiN9SIESNw4MABJCUl4ciRI/jkk09gaGiIbt26wcrKCn369MGwYcOwb98+nDp1Cr169YKPjw/ee+89AIC/vz9q1KiBHj16IC4uDjt27MDYsWMREhIClUql49EREVFZofVnoLZt24adO3fit99+w7Zt27TdHRERvaGuX7+Obt264d9//4WDgwPef/99HD16FA4ODgCAWbNmwcDAAJ06dUJWVhYCAgKwcOFCaX9DQ0Ns3rwZAwYMgI+PD8qVK4fAwEBMmDBBV0MiIqIySKsJVFpaGoKDg7FhwwaYmZkVWz8rKwtZWVnSOqeKJSKifKtXr37hdhMTEyxYsAALFiwoso6bmxu2bt36ukMjIqK3iNZu4RNCICgoCP3794e3t7dG+3CqWCIiIiIi0mclTqDGjBkDhULxwuXy5cuYN28e7t+/j7CwMI3bDgsLQ3p6urRcu3atpOERERERERFpTYlv4Rs+fDiCgoJeWKdSpUrYu3cvYmJiCjyY6+3tje7du2PFihUF9lOpVHyQl4iIiIiI9FaJEygHBwfpgd0XmTt3LiZNmiSt37hxAwEBAVizZg0aNWpU0m6JiIiIiIh0TmuTSFSoUEG2bm5uDgCoXLkyypcvr61uiYiIiIiItKZU3gNFRERERET0JtD6e6Dyubu7QwhRWt0RERERERG9drwCRUREREREpCEmUERERERERBpiAkVERERERKShUnsGioiIiPSD+5gtug6BiKjM4hUoIiIiIiIiDTGBIiIiIiIi0hATKCIiIiIiIg0xgSIiIiIiItJQ2ZhEokEtwMhE11EAR+N0HYGcPsXzXm1dR6CXKg85qusQ5PTpPB1N0HUERERERCXGK1BEREREREQaYgJFRERERESkISZQREREREREGmICRUREREREpCEmUERERERERBpiAkVERERERKQhJlBEREREREQaKhvvgSIiIiKi18J9zJZX2j9pSpvXFAlR2cQrUERERERERBpiAkVERERERKQhJlBEREREREQaYgJFRERERESkISZQREREREREGmICRUREREREpCEmUERERERERBpiAkVERGVCZGQkGjRoAAsLCzg6OqJDhw6Ij4+X1fH19YVCoZAt/fv3l9VJTk5GmzZtYGZmBkdHR4wcORI5OTmlORQiIirD+CJdIiIqEw4cOICQkBA0aNAAOTk5+Prrr+Hv74+LFy+iXLlyUr3g4GBMmDBBWjczM5P+n5ubizZt2kCtVuPIkSNISUlBz549YWxsjMmTJ5fqeIiIqGzS6hWoLVu2oFGjRjA1NYWNjQ06dOigze6IiOgNtn37dgQFBaFmzZqoXbs2oqOjkZycjFOnTsnqmZmZQa1WS4ulpaW0befOnbh48SJ+/vln1KlTB61atcLEiROxYMECZGdnl/aQiIioDNJaAvXbb7+hR48e6NWrF+Li4nD48GF8/vnn2uqOiIjeMunp6QAAW1tbWfnKlSthb2+PWrVqISwsDA8fPpS2xcTEwNPTE05OTlJZQEAAMjIycOHChUL7ycrKQkZGhmwhIqK3l1Zu4cvJycHgwYMxbdo09OnTRyqvUaOGNrojIqK3TF5eHoYMGYImTZqgVq1aUvnnn38ONzc3uLi44OzZsxg9ejTi4+Px+++/AwBSU1NlyRMAaT01NbXQviIjIxEREaGlkRARUVmjlQTq9OnT+Oeff2BgYIC6desiNTUVderUwbRp02QfdM/LyspCVlaWtM6/8hERUWFCQkJw/vx5/PHHH7Lyvn37Sv/39PSEs7MzWrRogatXr6Jy5cov1VdYWBiGDRsmrWdkZMDV1fXlAiciojJPK7fw/fXXXwCA8PBwjB07Fps3b4aNjQ18fX1x586dIveLjIyElZWVtPADioiInhcaGorNmzdj3759KF++/AvrNmrUCABw5coVAIBarUZaWpqsTv66Wq0utA2VSgVLS0vZQkREb68SJVBjxowpMD3s88vly5eRl5cHAPjmm2/QqVMn1K9fH1FRUVAoFFi3bl2R7YeFhSE9PV1arl279mqjIyKiN4YQAqGhoVi/fj327t2LihUrFrtPbGwsAMDZ2RkA4OPjg3PnzuHmzZtSnV27dsHS0pK3mRMRkUZKdAvf8OHDERQU9MI6lSpVQkpKCgD5M08qlQqVKlVCcnJykfuqVCqoVKqShERERG+JkJAQrFq1Chs3boSFhYX0zJKVlRVMTU1x9epVrFq1Cq1bt4adnR3Onj2LoUOHomnTpvDy8gIA+Pv7o0aNGujRowemTp2K1NRUjB07FiEhIfz8ISIijZQogXJwcICDg0Ox9erXrw+VSoX4+Hi8//77AIAnT54gKSkJbm5uLxcpERG91RYtWgTg6ctynxUVFYWgoCAolUrs3r0bs2fPxoMHD+Dq6opOnTph7NixUl1DQ0Ns3rwZAwYMgI+PD8qVK4fAwEDZe6OIiIheRCuTSFhaWqJ///4YP348XF1d4ebmhmnTpgEAunTpoo0uiYjoDSeEeOF2V1dXHDhwoNh23NzcsHXr1tcVFhERvWW0kkABwLRp02BkZIQePXrg0aNHaNSoEfbu3QsbGxttdUlERERERKRVWkugjI2NMX36dEyfPl1bXRAREREREZUqrUxjTkRERERE9CbS2hWo1+rEeUBhrOso6EWOxuk6Av30Xm1dRyDH81Q4fTpPPEdERER6jVegiIiIiIiINMQEioiIiIiISENl4xY+IiIiItIL7mO2vNL+SVPavKZIiHSDV6CIiIiIiIg0xASKiIiIiIhIQ0ygiIiIiIiINMQEioiIiIiISENMoIiIiIiIiDTEBIqIiIiIiEhDTKCIiIiIiIg0xASKiIiIiIhIQ0ygiIiIiIiINMQEioiIiIiISENGug6AiIiISsZ9zBZdh0BE9NbiFSgiIiIiIiIN8QoUEREREZWaV72CmjSlzWuKhOjl8AoUERERERGRhphAERERERERaYgJFBERERERkYaYQBEREREREWmICRQREREREZGGOAsfEREREZUZnMWPdI1XoIiIiIiIiDTEBIqIiIiIiEhDWkug/vzzT3z88cewt7eHpaUl3n//fezbt09b3REREREREWmd1hKotm3bIicnB3v37sWpU6dQu3ZttG3bFqmpqdrqkoiIiIiISKu0MonE7du3kZCQgGXLlsHLywsAMGXKFCxcuBDnz5+HWq3WRrdEREQaWbBgAaZNm4bU1FTUrl0b8+bNQ8OGDXUdFhGVAk5CQa9KK1eg7Ozs4OHhgZ9++gkPHjxATk4OFi9eDEdHR9SvX7/I/bKyspCRkSFbiIiIXqc1a9Zg2LBhGD9+PE6fPo3atWsjICAAN2/e1HVoRERUBmjlCpRCocDu3bvRoUMHWFhYwMDAAI6Ojti+fTtsbGyK3C8yMhIRERHaCImIiAgAMHPmTAQHB6NXr14AgB9++AFbtmzB8uXLMWbMGB1HR0T6jlewqEQJ1JgxY/D999+/sM6lS5fg4eGBkJAQODo64tChQzA1NcWPP/6Idu3a4cSJE3B2di5037CwMAwbNkxaT09PR4UKFZCDJ4AoSaREeiLnsa4jkBNPdB2BftKn86Qn5ygHT+MQ4s364ZudnY1Tp04hLCxMKjMwMICfnx9iYmIK3ScrKwtZWVnSenp6OgC80l0StcbveOl93yR52f997+VlPQREng6jISodFYaue6X9z0cEvKZIyp78n7u6/mxSiBJEcOvWLfz7778vrFOpUiUcOnQI/v7+uHv3LiwtLaVtVatWRZ8+fTT+C9/169fh6uqqaXhERPSaXbt2DeXLl9d1GK/NjRs38M477+DIkSPw8fGRykeNGoUDBw7g2LFjBfYJDw/n3RFERHpE159NJboC5eDgAAcHh2LrPXz4EMDTv+o9y8DAAHl5mv91ycXFBdeuXYOFhQUUCkVJQtVLGRkZcHV1xbVr12SJJf2Hx6h4PEbF4zEqXnHHSAiB+/fvw8XFRQfR6Zfn747Iy8vDnTt3YGdnp7efTW/a9wDHo7/epLEAb9Z43qSxAP+N5+LFizr/bNLKM1A+Pj6wsbFBYGAgxo0bB1NTUyxduhSJiYlo00bz+z4NDAzeqL985rO0tHwjvpC1iceoeDxGxeMxKt6LjpGVlVUpR6N99vb2MDQ0RFpamqw8LS2tyBliVSoVVCqVrMza2lpbIb5Wb9r3AMejv96ksQBv1njepLEAwDvvvFPgIk1p00rv9vb22L59OzIzM/Hhhx/C29sbf/zxBzZu3IjatWtro0siIqJiKZVK1K9fH3v27JHK8vLysGfPHtktfUREREXRyhUoAPD29saOHXxIloiI9MuwYcMQGBgIb29vNGzYELNnz8aDBw+kWfmIiIheRGsJFBWkUqkwfvz4AreC0H94jIrHY1Q8HqPivc3H6NNPP8WtW7cwbtw4pKamok6dOti+fTucnJx0Hdpr86adX45Hf71JYwHerPG8SWMB9Gs8JZqFj4iIiIiI6G2m2yewiIiIiIiIyhAmUERERERERBpiAkVERERERKQhJlBEREREREQaYgJFRERERESkISZQOpaVlYU6depAoVAgNjZW1+HojaSkJPTp0wcVK1aEqakpKleujPHjxyM7O1vXoenUggUL4O7uDhMTEzRq1AjHjx/XdUh6JTIyEg0aNICFhQUcHR3RoUMHxMfH6zosvTVlyhQoFAoMGTJE16G89Ur6vb1u3TpUq1YNJiYm8PT0xNatW2XbhRAYN24cnJ2dYWpqCj8/PyQkJEjbNfkZm5SUBIVCUWA5evSo3o0HANzd3QvEOmXKFFmds2fP4oMPPoCJiQlcXV0xdepUvRvL/v37Cz3uCoUCJ06cAKBf5+b333+Hv78/7Ozsivxd5vHjxwgJCYGdnR3Mzc3RqVMnpKWlyeokJyejTZs2MDMzg6OjI0aOHImcnBy9G8+dO3cwcOBAeHh4wNTUFBUqVMCgQYOQnp4uq1fY+Vm9erVejQUAfH19C8TZv39/WZ2ycm6K+r5QKBRYt26dVO9lzk0BgnRq0KBBolWrVgKAOHPmjK7D0Rvbtm0TQUFBYseOHeLq1ati48aNwtHRUQwfPlzXoenM6tWrhVKpFMuXLxcXLlwQwcHBwtraWqSlpek6NL0REBAgoqKixPnz50VsbKxo3bq1qFChgsjMzNR1aHrn+PHjwt3dXXh5eYnBgwfrOpy3Wkm/tw8fPiwMDQ3F1KlTxcWLF8XYsWOFsbGxOHfunFRnypQpwsrKSmzYsEHExcWJ9u3bi4oVK4pHjx4JITT7GZuYmCgAiN27d4uUlBRpyc7O1rvxCCGEm5ubmDBhgizWZ7/309PThZOTk+jevbs4f/68+OWXX4SpqalYvHixXo0lKytLNoaUlBTx5ZdfiooVK4q8vDy9Ozc//fSTiIiIEEuXLi3yd5n+/fsLV1dXsWfPHnHy5Enx3nvvicaNG0vbc3JyRK1atYSfn584c+aM2Lp1q7C3txdhYWF6N55z586Jjh07ik2bNokrV66IPXv2iKpVq4pOnTrJ6gEQUVFRsvPz7NerPoxFCCGaNWsmgoODZXGmp6dL28vSucnJySnwvRMRESHMzc3F/fv3pXolPTeFYQKlQ1u3bhXVqlUTFy5cYAKlgalTp4qKFSvqOgydadiwoQgJCZHWc3NzhYuLi4iMjNRhVPrt5s2bAoA4cOCArkPRK/fv3xdVq1YVu3btEs2aNWMCpWMl/d7u2rWraNOmjaysUaNGol+/fkIIIfLy8oRarRbTpk2Ttt+7d0+oVCrxyy+/FBnH8z9j839JL+lnk67G4+bmJmbNmlVkXAsXLhQ2NjYiKytLKhs9erTw8PDQu7E8Kzs7Wzg4OIgJEyZIZfpybp5VVEz37t0TxsbGYt26dVLZpUuXBAARExMjhHj6+5CBgYFITU2V6ixatEhYWlrKzpc+jKcwa9euFUqlUjx58kQqAyDWr19f7L75dDWW4j4Dyvq5qVOnjujdu7esrKTnpjC8hU9H0tLSEBwcjP/9738wMzPTdThlQnp6OmxtbXUdhk5kZ2fj1KlT8PPzk8oMDAzg5+eHmJgYHUam3/JvqXhbv26KEhISgjZt2si+nkg3XuZ7OyYmpsC5CwgIkOonJiYiNTVVVsfKygqNGjV64c+Lon7Gtm/fHo6Ojnj//fexadMmvR7PlClTYGdnh7p162LatGmy24xiYmLQtGlTKJVKWT/x8fG4e/eu3o0l36ZNm/Dvv/+iV69eBbbp+txo4tSpU3jy5ImsnWrVqqFChQpSOzExMfD09ISTk5Osn4yMDFy4cEGvxlOY9PR0WFpawsjISFYeEhICe3t7NGzYEMuXL8fT3931bywrV66Evb09atWqhbCwMDx8+FDWT1k9N6dOnUJsbCz69OlTYJum56YoRsVXoddNCIGgoCD0798f3t7eSEpK0nVIeu/KlSuYN28epk+frutQdOL27dvIzc2V/QADACcnJ1y+fFlHUem3vLw8DBkyBE2aNEGtWrV0HY7eWL16NU6fPi09S0G69TLf26mpqYXWT01NlbbnlxVV53mF/Yw1NzfHjBkz0KRJExgYGOC3335Dhw4dsGHDBrRv317vxjNo0CDUq1cPtra2OHLkCMLCwpCSkoKZM2dK7VSsWLFAG/nbbGxs9GYsz1q2bBkCAgJQvnx5qUxfzo0mUlNToVQqYW1tXWQ7RfWTv02fxlNYHBMnTkTfvn1l5RMmTMCHH34IMzMz7Ny5E1999RUyMzMxaNAgvRrL559/Djc3N7i4uODs2bMYPXo04uPj8fvvv7+wn/xthdGXc7Ns2TJUr14djRs3lpWX5NwUhQnUazRmzBh8//33L6xz6dIl7Ny5E/fv30dYWFgpRaY/ND1G1apVk9b/+ecftGzZEl26dEFwcLC2Q6Q3REhICM6fP48//vhD16HojWvXrmHw4MHYtWsXTExMdB0O6Ymifsba29tj2LBh0nqDBg1w48YNTJs2rchf0nXp2Vi9vLygVCrRr18/REZGQqVS6TCyl3f9+nXs2LEDa9eulZWXtXPzpsrIyECbNm1Qo0YNhIeHy7Z9++230v/r1q2LBw8eYNq0aSX6Jb00PJv4eXp6wtnZGS1atMDVq1dRuXJlHUb2ah49eoRVq1bJzkO+13FueAvfazR8+HBcunTphUulSpWwd+9exMTEQKVSwcjICFWqVAEAeHt7IzAwUMej0C5Nj1G+GzduoHnz5mjcuDGWLFmiw8h1y97eHoaGhgVmLUpLS4NardZRVPorNDQUmzdvxr59+2R/tX3bnTp1Cjdv3kS9evVgZGQEIyMjHDhwAHPnzoWRkRFyc3N1HeJb52W+t9Vq9Qvr5/+rSZsl/RnbqFEjXLlyRW/H83ysOTk50l0eRfXzbB/6NpaoqCjY2dlplBTp4txoQq1WIzs7G/fu3SuynZKeG12OJ9/9+/fRsmVLWFhYYP369TA2Nn5h/UaNGuH69evIysoqsE3XY3k+TgDS11JZPDcA8Ouvv+Lhw4fo2bNnsXVfdG6KwgTqNXJwcEC1atVeuCiVSsydOxdxcXGIjY1FbGysNE3jmjVr8N133+l4FNql6TECnv5V1NfXF/Xr10dUVBQMDN7eL1elUon69etjz549UlleXh727NkDHx8fHUamX4QQCA0Nxfr167F3794Ct+u87Vq0aIFz585JP3tiY2Ph7e2N7t27IzY2FoaGhroO8a3zMt/bPj4+svoAsGvXLql+xYoVoVarZXUyMjJw7NgxWZsv8zM2NjYWzs7OejmewmI1MDCAo6Oj1M/Bgwfx5MkTWT8eHh4Fbt/Th7EIIRAVFYWePXsW+8t5/nhL+9xoon79+jA2Npa1Ex8fj+TkZKkdHx8fnDt3Djdv3pT1Y2lpiRo1aujVeICn58zf3x9KpRKbNm3S6Ip+bGwsbGxsCr0aqsuxFBYnAOlrqaydm3zLli1D+/bt4eDgUGzdF52bIr3SFBT0WrzsbDpvsuvXr4sqVaqIFi1aiOvXr8ummnxbrV69WqhUKhEdHS0uXrwo+vbtK6ytrWUz47ztBgwYIKysrMT+/ftlXzMPHz7UdWh6i7Pw6V5x39s9evQQY8aMkeofPnxYGBkZienTp4tLly6J8ePHFzpVtrW1tdi4caM4e/as+Pjjj2VTZWvyMzY6OlqsWrVKXLp0SVy6dEl89913wsDAQCxfvlzvxnPkyBExa9YsERsbK65evSp+/vln4eDgIHr27Cm1ce/ePeHk5CR69Oghzp8/L1avXi3MzMyKnca8tMeSb/fu3QKAuHTpUoG49Onc/Pvvv+LMmTNiy5YtAoBYvXq1OHPmjOxrqX///qJChQpi79694uTJk8LHx0f4+PhI2/Onyvb39xexsbFi+/btwsHBQaOpskt7POnp6aJRo0bC09NTXLlyRfa9k5OTI4QQYtOmTWLp0qXi3LlzIiEhQSxcuFCYmZmJcePG6dVYrly5IiZMmCBOnjwpEhMTxcaNG0WlSpVE06ZNy+S5yZeQkCAUCoXYtm1bgbhe5twUhgmUHmACVVBUVJQAUOjyNps3b56oUKGCUCqVomHDhuLo0aO6DkmvFPU1ExUVpevQ9BYTKP3wou/tZs2aicDAQFn9tWvXinfffVcolUpRs2ZNsWXLFtn2vLw88e233wonJyehUqlEixYtRHx8vLRdk5+x0dHRonr16sLMzExYWlqKhg0byqai1qfxnDp1SjRq1EhYWVkJExMTUb16dTF58mTx+PFjWTtxcXHi/fffFyqVSrzzzjtiypQpejeWfN26dZO9K+lZ+nRuivpaGj9+vFTn0aNH4quvvhI2NjbCzMxMfPLJJwV+6U1KShKtWrUSpqamwt7eXgwfPlw2Lbi+jGffvn1Ffu8kJiYKIZ6+Z61OnTrC3NxclCtXTtSuXVv88MMPIjc3V6/GkpycLJo2bSpsbW2FSqUSVapUESNHjpS9B0qIsnNu8oWFhQlXV9dCj/fLnpvnKYQo4bx9REREREREb6m396ESIiIiIiKiEmICRUREREREpCEmUERERERERBpiAkVERERERKQhJlBEREREREQaYgJFRERERESkISZQREREREREGmICRUREREREpCEmUERERERERBpiAkVERERERKQhJlBEREREREQa+n95DTxdKCaOigAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0cAAAGHCAYAAAB7zSkeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAABy90lEQVR4nO3de1zO9/8/8MfVuXTW4aotFZpj1EIrRqbJYQdmjE9D1mJWhhz7IIUtsyFymn0o+4wx2xyGITltJEQOoeFTMjqwVCt0ul6/P/x6f106KLq6Ko/77fa+6Xq/X4fn67q6XD2v1/v9esuEEAJEREREREQvOA11B0BERERERNQQMDkiIiIiIiICkyMiIiIiIiIATI6IiIiIiIgAMDkiIiIiIiICwOSIiIiIiIgIAJMjIiIiIiIiAEyOiIiIiIiIADA5IiIiIiIiAsDkiIiInuLw4cOQyWQ4fPiwukNR4uDgAD8/P3WHUSdiYmIgk8mQlpam7lCIiF5oTI6IiEhl9uzZg7CwMHWHQXVg3bp1aNeuHfT09ODk5ISoqKga1y0qKsKMGTNga2sLfX19uLu7IzY2ttKyx48fR48ePWBgYAC5XI7PPvsMBQUFz9UmEVFNMTkiIiKV2bNnD8LDw9UdBj2nb775Bh9//DE6dOiAqKgoeHh44LPPPsOXX35Zo/p+fn5YsmQJfH19sWzZMmhqamLAgAH4448/lMolJSWhT58+uH//PpYsWYKPP/4Ya9euxdChQ5+5TSKi2tBSdwBERI2NQqFAcXEx9PT01B1Kte7fvw8DA4MK+0tLS6FQKKCjo6OGqKixefDgAWbNmoWBAwfip59+AgAEBARAoVBg/vz5GDt2LMzMzKqsf/LkSWzevBlfffUVpk6dCgAYNWoUOnbsiOnTp+P48eNS2X//+98wMzPD4cOHYWxsDODR6ZMBAQHYv38/+vbtW+s2iYhqgzNHRKQ2Dx48QNu2bdG2bVs8ePBA2p+TkwMbGxt4enqirKys0roFBQVo1qwZJk6cWOHYX3/9BU1NTURERFTb/9dffw1PT080b94c+vr6cHNzk/74e5xMJkNQUBA2btyIDh06QFdXF3v37gUA3Lp1Cx999BGsra2hq6uLDh06YP369TV+Dr7//nt069YNBgYGMDMzQ8+ePbF//36lMqtWrZL6tbW1RWBgIHJzc5XKeHl5oWPHjkhMTETPnj1hYGCAf//730hLS4NMJsPXX3+NyMhItGrVCrq6urh06RIA4MqVK3j//fdhbm4OPT09dOnSBTt37nxq3L///juGDh2KFi1aQFdXF3Z2dpg8ebLS6+jn54eVK1dKz2H5Vk6hUCAyMhIdOnSAnp4erK2tMW7cONy7d0+pLyEEFixYgJdffhkGBgbo3bs3kpOTK43r+vXruH79+lPjL7/G5+jRoxg3bhyaN28OY2NjjBo1qkL/T/rpp58gk8lw5MiRCse++eYbyGQyXLx4EQBw/vx5+Pn5oWXLltDT04NcLsdHH32Ev//++6kxymSySk9JrOxaq9zcXEyaNAl2dnbQ1dVF69at8eWXX0KhUDy1n6c5dOgQ/v77b3z66adK+wMDA1FYWIjdu3dXW/+nn36CpqYmxo4dK+3T09ODv78/4uPjcfPmTQBAfn4+YmNj8eGHH0qJEfAo6TE0NMSPP/5Y6zaJiGqLM0dEpDb6+vrYsGEDunfvjlmzZmHJkiUAHv3RlZeXh5iYGGhqalZa19DQEIMHD8aWLVuwZMkSpXI//PADhBDw9fWttv9ly5bhnXfega+vL4qLi7F582YMHToUu3btwsCBA5XKHjx4ED/++COCgoJgYWEBBwcHZGVl4bXXXpOSJ0tLS/z222/w9/dHfn4+Jk2aVG3/4eHhCAsLg6enJ+bNmwcdHR0kJCTg4MGD0jfkYWFhCA8Ph7e3N8aPH4+UlBSsXr0ap06dwrFjx6CtrS219/fff6N///4YPnw4PvzwQ1hbW0vHoqOj8fDhQ4wdOxa6urowNzdHcnIyunfvjpdeegkzZ85Es2bN8OOPP2LQoEH4+eefMXjw4Cpj37p1K+7fv4/x48ejefPmOHnyJKKiovDXX39h69atAIBx48bh9u3biI2NxX//+98KbYwbNw4xMTEYM2YMPvvsM6SmpmLFihU4e/as0thCQ0OxYMECDBgwAAMGDMCZM2fQt29fFBcXV2izT58+AFDjhQ2CgoJgamqKsLAw6bm9ceOGtAhFZQYOHCj9sd6rVy+lY1u2bEGHDh3QsWNHAEBsbCz+97//YcyYMZDL5UhOTsbatWuRnJyMEydOVNlHbdy/fx+9evXCrVu3MG7cOLRo0QLHjx9HSEgIMjIyEBkZKZW9d+9elV84PM7AwECadTx79iwAoEuXLkpl3NzcoKGhgbNnz+LDDz+ssq2zZ8/ilVdeUUp4AKBbt24AHp1KZ2dnhwsXLqC0tLRCPzo6OnBxcZHiqE2bRES1JoiI1CwkJERoaGiIo0ePiq1btwoAIjIy8qn19u3bJwCI3377TWl/p06dRK9evZ5a//79+0qPi4uLRceOHcUbb7yhtB+A0NDQEMnJyUr7/f39hY2Njbh7967S/uHDhwsTE5MK7T/u6tWrQkNDQwwePFiUlZUpHVMoFEIIIbKzs4WOjo7o27evUpkVK1YIAGL9+vXSvl69egkAYs2aNUptpaamCgDC2NhYZGdnKx3r06ePcHZ2Fg8fPlTq29PTUzg5OUn7Dh06JACIQ4cOSfsqG1tERISQyWTixo0b0r7AwEBR2UfN77//LgCIjRs3Ku3fu3ev0v7y52DgwIHS8yKEEP/+978FADF69Gil+vb29sLe3r5Cf0+Kjo4WAISbm5soLi6W9i9atEgAEDt27Ki2/ogRI4SVlZUoLS2V9mVkZAgNDQ0xb948aV9lz9MPP/wgAIijR49WiCc1NVXaB0DMnTu3Qn17e3ulcc+fP180a9ZM/Pnnn0rlZs6cKTQ1NUV6erpSXQBP3R7vNzAwUGhqalb6PFhaWorhw4dXeqxchw4dKrynhBAiOTlZ6Xe2/L3/+PNSbujQoUIul9e6TSKi2uJpdUSkdmFhYejQoQNGjx6NTz/9FL169cJnn3321Hre3t6wtbXFxo0bpX0XL17E+fPnq/0mu5y+vr70871795CXl4fXX38dZ86cqVC2V69eaN++vfRYCIGff/4Zb7/9NoQQuHv3rrT5+PggLy+v0nbKbd++HQqFAqGhodDQUP6vuHw24cCBAyguLsakSZOUygQEBMDY2LjC6Uy6uroYM2ZMpf0NGTIElpaW0uOcnBwcPHgQw4YNwz///CPF/vfff8PHxwdXr17FrVu3qoz/8eeusLAQd+/ehaenJ4QQSt/wV2Xr1q0wMTHBm2++qfTcubm5wdDQEIcOHVJ6DiZMmKA0y1LVrFxaWlqtlsMeO3as0uzb+PHjoaWlhT179lRb74MPPkB2drbS8uY//fQTFAoFPvjgA2nf48/Tw4cPcffuXbz22msAUO3vR21s3boVr7/+OszMzJSeS29vb5SVleHo0aNS2Y0bNyI2Nvap26hRo6Q6Dx48qPL6ND09PaVTKSvz4MED6OrqVlq3/Pjj/1ZV9vF+atomEVFt8bQ6IlI7HR0drF+/Hl27doWenh6io6NrdLqRhoYGfH19sXr1amnxgY0bN0JPT6/S1a2etGvXLixYsABJSUkoKiqS9lfWt6Ojo9LjO3fuIDc3F2vXrsXatWsrbT87O7vKvq9fvw4NDQ2lhOtJN27cAAC0adNGab+Ojg5atmwpHS/30ksvVflH7JPxX7t2DUIIzJkzB3PmzKky/pdeeqnSY+np6QgNDcXOnTsrXKOTl5dX5ZjKXb16FXl5ebCysqqyb+D/ngMnJyel45aWltUuAlBTT7ZraGgIGxsbKcHKy8tT+kNbR0cH5ubm6NevH0xMTLBlyxbpVL4tW7bAxcUFr7zyilQ+JycH4eHh2Lx5c4Xfh5o8TzVx9epVnD9/Xin5fdzj/Xbv3r3W7evr61d6CiPwKOF7PAGsqv7j76/H65Yff/zfqso+3k9N2yQiqi0mR0TUIOzbtw/Aoz9url69WuGP+aqMGjUKX331FbZv344RI0Zg06ZNeOutt2BiYlJtvd9//x3vvPMOevbsiVWrVsHGxgba2tqIjo7Gpk2bKpR/8o+t8gvdP/zwQ4wePbrSPjp16lSjMdSV6v4grCr+qVOnwsfHp9I6rVu3rnR/WVkZ3nzzTeTk5GDGjBlo27YtmjVrhlu3bsHPz69GiwAoFApYWVkpzfo9rqo/9OvbxIkTsWHDBulxr169cPjwYejq6mLQoEHYtm0bVq1ahaysLBw7dgxffPGFUv1hw4bh+PHjmDZtGlxcXGBoaAiFQoF+/fo982IJT14zpFAo8Oabb2L69OmVln88Wbtz506NrjkyNDSEoaEhAMDGxgZlZWXIzs5WSmaLi4vx999/w9bWttq2bGxsKp2FzMjIAACpvo2NjdL+J8s+3k9N2yQiqi0mR0SkdufPn8e8efMwZswYJCUl4eOPP8aFCxeemuAAQMeOHeHq6oqNGzfi5ZdfRnp6eo1uTvnzzz9DT08P+/btUzo9Jzo6ukYxW1pawsjICGVlZfD29q5Rnce1atUKCoUCly5dgouLS6Vl7O3tAQApKSlo2bKltL+4uBipqanP1G+58va0tbVr3c6FCxfw559/YsOGDUqnX1V2A86qZgBbtWqFAwcOoHv37tUmdeXPwdWrV5Wegzt37jx1VbmauHr1Knr37i09LigoQEZGBgYMGAAAmD59utIpmo/PVn3wwQfYsGED4uLicPnyZQghlE6pu3fvHuLi4hAeHo7Q0FClPmvCzMyswqqExcXFFZKHVq1aoaCgoEavY9euXSvMOFZm7ty50kp55b+fp0+flp6X8scKhaLK399yLi4uOHToEPLz85UWUEhISFBqv2PHjtDS0sLp06cxbNgwqVxxcTGSkpKU9tW0TSKi2uI1R0SkViUlJfDz84OtrS2WLVuGmJgYZGVlYfLkyTVuY+TIkdi/fz8iIyPRvHlz9O/f/6l1NDU1IZPJlL5FT0tLw/bt22vUp6amJoYMGYKff/5ZWrb5cXfu3Km2/qBBg6ChoYF58+ZVmEEQQgB4dE2Vjo4Oli9fLu0DgHXr1iEvL6/Cinq1YWVlBS8vL3zzzTeVflNfXfzlKwM+HpMQAsuWLatQtlmzZgBQ4Y/8YcOGoaysDPPnz69Qp7S0VCrv7e0NbW1tREVFKfX3+Apsj6vpUt7l1q5di5KSEunx6tWrUVpaKv0OtW/fHt7e3tLm5uYmlfX29oa5uTm2bNmCLVu2oFu3bkoznpU9T9XF/qRWrVopXS9UHu+TMz/Dhg1DfHy8NPv6uNzcXJSWlkqPn+WaozfeeAPm5uZYvXq1UturV6+GgYGB0u/h3bt3ceXKFdy/f1/a9/7776OsrEzp9NOioiJER0fD3d1dWlXOxMQE3t7e+P777/HPP/9IZf/73/+ioKBA6VTZmrZJRFRraloIgohICCFEaGiokMlk4uDBg9K+BQsWCABi9+7d0r7yVdeeXJ1MCCEyMzOFlpaWACDGjx9fo37j4uIEAPH666+L1atXi/DwcGFlZSU6depUYXU1ACIwMLDSfu3t7YWBgYGYOHGi+Oabb0RERIQYOnSoMDMze2oMc+bMEQCEp6en+Prrr0VUVJQYNWqUmDlzplRm7ty5AoDo27evWLFihZgwYYLQ1NQUXbt2VVplrVevXqJDhw4V+ih/3r766qsKx5KTk4WZmZlo3ry5mDlzpli7dq2YP3++GDBggOjUqZNU7snV6oqLi0WrVq2EhYWF+Pzzz0VUVJTw8vISnTt3FgBEdHS0VPfHH38UAMTIkSPF999/L3744Qfp2Lhx4wQA0b9/f7F06VKxYsUKMXHiRGFrayu2bt0qlQsJCREAxIABA8SKFSuEv7+/sLW1FRYWFs+9Wp2zs7N4/fXXRVRUlAgKChIaGhqiR48eSivjVefjjz8WhoaGQiaTicWLF1c43rNnT2FgYCBmzZolVq1aJQYNGiQ9T4+vCFfZanVr1qwRAMR7770nVq9eLT755BPh6OhYYdyFhYXi1VdfFVpaWuLjjz8Wq1evFl9//bUYPXq0aNasmbhz506NxlKdlStXCgDi/fffF99++60YNWqUACA+//xzpXLlv6+Pr2woxKPV5rS0tMS0adPEN998Izw9PYWWlpY4cuSIUrnExEShq6srXF1dxerVq8WsWbOEnp6e6Nu3b4WYatomEVFtMDkiIrVJTEwUWlpaYsKECUr7S0tLRdeuXYWtra24d++eEEKICxcuCABKicPjBgwYIACI48eP17j/devWCScnJ6Grqyvatm0roqOjpT/uHldVciSEEFlZWSIwMFDY2dkJbW1tIZfLRZ8+fcTatWtrFMP69euFq6ur0NXVFWZmZqJXr14iNjZWqcyKFStE27Zthba2trC2thbjx4+Xnpdyz5IcCSHE9evXxahRo4RcLhfa2tripZdeEm+99Zb46aefpDKVLeV96dIl4e3tLQwNDYWFhYUICAgQ586dq5AclZaWigkTJghLS0shk8kqPLdr164Vbm5uQl9fXxgZGQlnZ2cxffp0cfv2balMWVmZCA8PFzY2NkJfX194eXmJixcvVljSWojaJ0dHjhwRY8eOFWZmZsLQ0FD4+vqKv//++6n1y8XGxgoAQiaTiZs3b1Y4/tdff4nBgwcLU1NTYWJiIoYOHSpu375do+SorKxMzJgxQ1hYWAgDAwPh4+Mjrl27Vum4//nnHxESEiJat24tdHR0hIWFhZR0P55EP4+1a9eKNm3aCB0dHdGqVSuxdOnSCklkVcnRgwcPxNSpU4VcLhe6urqia9euYu/evZX28/vvvwtPT0+hp6cnLC0tRWBgoMjPz69QrjZtEhHVlEyIJ+b7iYgaoFWrVmH69Om4fv260s1Nyw0ePBgXLlzAtWvX1BAdNTblN589depUhZuOEhHRi4vXHBFRo3Do0CF89tlnlSZGGRkZ2L17N0aOHKmGyIiIiKip4Gp1RNQobN26tcK+1NRUHDt2DP/5z3+gra2NcePGqSEyIiIiaio4c0REjdaRI0cwcuRIpKamYsOGDZDL5eoOiYiIiBoxXnNEREREREQEzhwREREREREBYHJEREREREQEgMkRERERERERACZHREREREREAJgcERERERERAWByREREREREBIDJEREREREREQAmR0RERERERACYHBEREREREQFgckRERERERASAyREREREREREAJkdEREREREQAmBwREREREREBYHJEREREREQEgMkRERERERERACZHREREREREAJgcERERERERAWByREREREREBIDJEREREREREQBAS90B1DWFQoHbt2/DyMgIMplM3eEQEb0whBD4559/YGtrCw0Nfvf2OH42ERGpR20/m5pccnT79m3Y2dmpOwwiohfWzZs38fLLL6s7jAaFn01EROpV08+mJpccGRkZAQB6YAC0oK3maIiIVE9DX0/dIQAASkUJjj78Rfp/mP5P+XNy8+ZNGBsbqzkaIqL6VVhYCFtbWwCPvixq1qxZvfWdn58POzu7Gn82NbnkqPx0BS1oQ0vG5IiImj4NmY66Q1DC08YqKn9OjI2NmRwR0QtHU1NT+tnY2Lhek6NyNf1s4knhREREREREYHJEREREREQEgMkRERERERERgCZ4zRERERER1YxCoUBxcbG6w6AmrqioCPb29tLPj1+D9Ly0tbXrtD0mR0REREQvoOLiYqSmpkKhUKg7FGriFAoF1qxZAwDIyMio83vhmZqaQi6X18mCQEyOiIiIiF4wQghkZGRAU1MTdnZ2vHEzqVRZWRkePHgAAHBwcKizmR4hBO7fv4/s7GwAgI2NzXO3yeSIiIgaraNHj+Krr75CYmIiMjIysG3bNgwaNEipzOXLlzFjxgwcOXIEpaWlaN++PX7++We0aNECAPDw4UNMmTIFmzdvRlFREXx8fLBq1SpYW1tLbaSnp2P8+PE4dOgQDA0NMXr0aEREREBLix+j1DiVlpbi/v37sLW1hYGBgbrDoSaurKxM+llPT69OT4PT19cHAGRnZ8PKyuq52+bXBERE1GgVFhaic+fOWLlyZaXHr1+/jh49eqBt27Y4fPgwzp8/jzlz5kBP7/9unDt58mT8+uuv2Lp1K44cOYLbt2/jvffek46XlZVh4MCBKC4uxvHjx7FhwwbExMQgNDRU5eMjUpXyP1Z1dBrWfdKInkV5gl9SUvLcbT1zcnT06FG8/fbbsLW1hUwmw/bt25WOCyEQGhoKGxsb6Ovrw9vbG1evXn1quytXroSDgwP09PTg7u6OkydPPmuIRETUxPXv3x8LFizA4MGDKz0+a9YsDBgwAIsWLYKrqytatWqFd955B1ZWVgCAvLw8rFu3DkuWLMEbb7wBNzc3REdH4/jx4zhx4gQAYP/+/bh06RK+//57uLi4oH///pg/fz5WrlzJC9mp0eNNm6kpqMvf42dOjp72bd2iRYuwfPlyrFmzBgkJCWjWrBl8fHzw8OHDKtvcsmULgoODMXfuXJw5cwadO3eGj4+PdB4hERFRTSkUCuzevRuvvPIKfHx8YGVlBXd3d6Uv8xITE1FSUgJvb29pX9u2bdGiRQvEx8cDAOLj4+Hs7Kx0mp2Pjw/y8/ORnJxcad9FRUXIz89X2oiIqOF75uSoum/rhBCIjIzE7Nmz8e6776JTp0747rvvcPv27QozTI9bsmQJAgICMGbMGLRv3x5r1qyBgYEB1q9f/6xhEhHRCyo7OxsFBQVYuHAh+vXrh/3792Pw4MF47733cOTIEQBAZmYmdHR0YGpqqlTX2toamZmZUpnHE6Py4+XHKhMREQETExNps7Ozq+PRERGRKqjkmqPU1FRkZmYqfRNnYmICd3d36Zu4JxUXFyMxMVGpjoaGBry9vausA/DbOSIiqlz58sTvvvsuJk+eDBcXF8ycORNvvfWWtKSsqoSEhCAvL0/abt68qdL+iKhpqeySldoICwuDi4tLncVTV27fvg0tLS0kJSWpO5QqqWSZnfJv0ir7pq2qb9nu3r2LsrKySutcuXKlyr4iIiIQHh7+nBETEVFTY2FhAS0tLbRv315pf7t27fDHH38AAORyOYqLi5Gbm6s0e5SVlQW5XC6VefL616ysLOlYZXR1daGrq1tXQ6E64DBz93PVT1s4sI4ioefh5eUFFxcXREZG1lmbfn5+yM3Nfa5k5FnbCAsLw/bt2yskCxkZGTAzM6tRGzKZrMJKnVOnTsWECRNqFQs90uhXq+O3c0REVBkdHR107doVKSkpSvv//PNP6U7tbm5u0NbWRlxcnHQ8JSUF6enp8PDwAAB4eHjgwoULSte/xsbGwtjYuELiRURUF+Ry+XN9wWJoaIjmzZvXYUQNh6oXwlFJclT+TVr5N2vlHv8m7kkWFhbQ1NSsVR3g0bdzxsbGShsREb0YCgoKkJSUJH3rmpqaiqSkJKSnpwMApk2bhi1btuDbb7/FtWvXsGLFCvz666/49NNPATw65dvf3x/BwcE4dOgQEhMTMWbMGHh4eOC1114DAPTt2xft27fHyJEjce7cOezbtw+zZ89GYGAgZ4eI6pGfnx+OHDmCZcuWQSaTQSaTIS0tDQBw8eJF9O/fH4aGhrC2tsbIkSNx9+5dqe5PP/0EZ2dn6Ovro3nz5vD29kZhYSHCwsKwYcMG7NixQ2rz8OHDlfb/LG3MmDEDr7zyCgwMDNCyZUvMmTNHWm46JiYG4eHhOHfunFQvJiYGgPJpdcXFxQgKCoKNjQ309PRgb2+PiIgIAI9uqAoAgwcPhkwmkx5Xdlrd+vXr0aFDB+jq6sLGxgZBQUFVPtdeXl6YNGmS0r5BgwbBz89Peuzg4IAvvvgCH330EYyMjNCiRQusXbtWqc7Jkyfh6uqKZs2aYdSoURW+rAKe/tp5eXkhKCgIkyZNgoWFBXx8fKqMuy6oJDlydHSEXC5X+iYuPz8fCQkJ0jdxT9LR0YGbm5tSHYVCgbi4uCrrEBHRi+306dNwdXWFq6srACA4OBiurq7SPYgGDx6MNWvWYNGiRXB2dsZ//vMf/Pzzz+jRo4fUxtKlS/HWW29hyJAh6NmzJ+RyOX755RfpuKamJnbt2gVNTU14eHjgww8/xKhRozBv3rz6HSzRC27ZsmXw8PBAQEAAMjIykJGRATs7O+Tm5uKNN96Aq6srTp8+jb179yIrKwvDhg0D8OgUtREjRuCjjz7C5cuXcfjwYbz33nsQQmDq1KkYNmwY+vXrJ7Xp6elZoe9nbcPIyAgxMTG4dOkSli1bhm+//RZLly4FAHzwwQeYMmUKOnToINX74IMPKvS9fPly7Ny5Ez/++CNSUlKwceNGKQk6deoUACA6OhoZGRnS4yetXr0agYGBGDt2LC5cuICdO3eidevWz/2aLF68GF26dMHZs2fx6aefYvz48VICVFBQgLfeegvt27fHyZMnERAQgGXLlinVf9prV27Dhg3Q0dHBsWPHVH7N6DNfc1RQUIBr165Jj8u/rTM3N0eLFi0wadIkLFiwAE5OTnB0dMScOXNga2urdD5knz59MHjwYClzDQ4OxujRo9GlSxd069YNkZGRKCwsxJgxY559hERE1GR5eXlBCFFtmY8++ggfffRRlcf19PSwcuXKKm9NAQD29vbYs2fPM8dJRM/PxMQEOjo6MDAwUDqraMWKFXB1dcUXX3wh7Vu/fj3s7Ozw559/oqCgAKWlpXjvvfekU2qdnZ2lsvr6+igqKqr2TKWMjIxnamP27NnSzw4ODpg6dSo2b96M6dOnQ19fH4aGhtDS0qq27/T0dDg5OaFHjx6QyWRS/wBgaWkJADA1Na22jQULFmDKlCmYOHGitK9r165Vlq+pAQMGSDPxM2bMwNKlS3Ho0CG0adMGmzZtgkKhwLp166CtrY3i4mJkZ2dj4cKFUv2nvXavvPIKAMDJyQmLFi167nhr4pmTo9OnT6N3797S4+DgYADA6NGjERMTg+nTp6OwsBBjx45Fbm4uevTogb179yrdlfz69etK02YffPAB7ty5g9DQUGRmZsLFxQV79+6tsEgDEREREREAnDt3DocOHYKhoWGFY9evX0ffvn3Rp08fODs7w8fHB3379sX7779f4wUPAKBz587P1MaWLVuwfPlyXL9+XUrSansJiJ+fH9588020adMG/fr1w1tvvYW+ffvWuH52djZu376NPn361KrfmujUqZP0s0wmg1wul67PvHz5Mjp16gQ9PT2UlZUBUE4ogae/duXJkZubW53HXpVnTo6e9m2dTCbDvHnzqj3toPw80ccFBQVVew4kEREREVG5goICvP322/jyyy8rHLOxsYGmpiZiY2Nx/Phx7N+/H1FRUZg1axYSEhLg6OhYoz6epY34+Hj4+voiPDwcPj4+MDExwebNm7F48eJaje/VV19FamoqfvvtNxw4cADDhg2Dt7c3fvrppxrV19fXr1V/wKPb6Tz5d375tVKP09bWVnosk8mk2yjUxNNeu3LNmjWrcZvPq9GvVkdERERELwYdHR1pFqLcq6++iuTkZDg4OKB169ZKW/kf1TKZDN27d0d4eDjOnj0LHR0dbNu2rco2K1PbNo4fPw57e3vMmjULXbp0gZOTE27cuPHU8VTG2NgYH3zwAb799lts2bIFP//8M3JycgA8SlCqa8PIyAgODg5K1/U/jaWlJTIyMqTHZWVluHjxYo3rA49um3D+/Hk8fPhQ2vdkGzV57eobkyMiIiIiahQcHByQkJCAtLQ03L17FwqFAoGBgcjJycGIESNw6tQpXL9+Hfv27cOYMWNQVlaGhIQEfPHFFzh9+jTS09Pxyy+/4M6dO2jXrp3U5vnz55GSkoK7d+9WOkPyLG04OTkhPT0dmzdvxvXr17F8+XIpmXp8POXX7d+9exdFRUUV+l6yZAl++OEHXLlyBX/++Se2bt0KuVwu3ZutPPHJzMzEvXv3Kn3ewsLCsHjxYixfvhxXr17FmTNnEBUVVeXz/MYbb2D37t3YvXs3rly5gvHjxyM3N7cmL5HkX//6F2QyGQICAnDp0iUcO3YM33//vVKZp7126qCSm8ASERERUePzvDfLra3a3lx36tSpGD16NNq3b48HDx4gNTUVDg4OOHbsGGbMmIG+ffuiqKgI9vb26NevHzQ0NGBsbIyjR48iMjIS+fn5sLe3x+LFi9G/f38AQEBAAA4fPowuXbqgoKAAhw4dgpeXl1K/z9LGO++8g8mTJyMoKAhFRUUYOHAg5syZg7CwMKndIUOG4JdffkHv3r2Rm5uL6OhopeWygUczP4sWLcLVq1ehqamJrl27Ys+ePdDQeDTHsXjxYgQHB+Pbb7/FSy+9VOllK6NHj8bDhw+xdOlSTJ06FRYWFnj//ferfJ4/+ugjnDt3DqNGjYKWlhYmT56stNZATRgaGuLXX3/FJ598gi5dusDBwQFBQUGYMWOGVMbW1rba104dZOJpy/w0Mvn5+TAxMYEX3oWWTPvpFYiIGjmNZzifXBVKRTEOPtiCvLw83nPuCeWfTXxu1Od5/+iv7R/xDd3Dhw+RmpoKR0dHpcWyGnpyRI1TWVkZzp49CwBwdXWFpqZmnbZf1e8zUPv/f3laHREREREREZgcERERERERAWByREREREREBIDJEREREREREQAmR0RERET0Ajl8+DBkMpm0NHVMTIy0LHZl0tLSIJPJkJSUVC/x1daT41EFBwcHREZGqqz9hoRLeRNRtTQ7vKLuECTifzfVHYJEpq/39EL1pCyn8vta1DeFqHhvECIiqjteXl5wcXGp90Tl1KlTarspa31jckRERERERFWytLRUdwj1hskREREREQFo+Pcd8vLyQqdOnaCnp4f//Oc/0NHRwSeffCLdWDUtLQ2Ojo44e/YsXFxcAAC5ubkwMzOr9OautfG///0PkydPRkJCApycnLBmzRp4eHigsLAQNjY2WL9+vdKNVbdv3w5fX19kZmbi77//hqOjI3744QcsX74cZ86cQevWrbFy5Ur06tVLqnPkyBFMmzYN586dg7m5OUaPHo0FCxZAS0sLfn5+OHLkCI4cOYJly5YBAFJTU6W6iYmJmDFjBi5dugQXFxdER0ejTZs20vEdO3YgPDwcly5dgq2tLUaPHo1Zs2ZBS0sLQgiEh4dj/fr1yMrKQvPmzfH+++9j+fLlAB6dVjdp0iRMmjTpqWUbO15zRERERESNxoYNG9CsWTMkJCRg0aJFmDdvHmJjY1Xe76xZszB16lQkJSXhlVdewYgRI1BaWopmzZph+PDhiI6OViofHR2N999/H0ZGRtK+adOmYcqUKTh79iw8PDzw9ttv4++//wYA3Lp1CwMGDEDXrl1x7tw5rF69GuvWrcOCBQsAAMuWLYOHhwcCAgKQkZGBjIwM2NnZKcW3ePFinD59GlpaWvjoo4+kY7///jtGjRqFiRMn4tKlS/jmm28QExODzz//HADw888/Y+nSpfjmm29w9epVbN++Hc7OzpU+D7Up2xhx5oiIiIiIGo1OnTph7ty5AAAnJyesWLECcXFxePPNN1Xa79SpUzFw4KOZtfDwcHTo0AHXrl1D27Zt8fHHH8PT0xMZGRmwsbFBdnY29uzZgwMHDii1ERQUhCFDhgAAVq9ejb1792LdunWYPn06Vq1aBTs7O6xYsQIymQxt27bF7du3MWPGDISGhsLExAQ6OjowMDCAXC6vEN/nn38uzULNnDkTAwcOxMOHD6Gnp4fw8HDMnDkTo0ePBgC0bNkS8+fPx/Tp0zF37lykp6dDLpfD29sb2traaNGiBbp161bp81Cbso0RZ46IiIiIqNHo1KmT0uPyZKQ++7WxsQEAqd9u3bqhQ4cO2LBhAwDg+++/h729PXr27KnUhoeHh/SzlpYWunTpgsuXLwMALl++DA8PD8hkMqlM9+7dUVBQgL/++uu54jt37hzmzZsHQ0NDaSufgbp//z6GDh2KBw8eoGXLlggICMC2bdtQWlpaaT+1KdsYMTkiIiIiokZDW1tb6bFMJoNCoQAAaGg8+tNWCCEdLympm5U0H++3PIEp7xcAPv74Y8TExAB4dErdmDFjlBIdVasuvoKCAoSHhyMpKUnaLly4gKtXr0JPTw92dnZISUnBqlWroK+vj08//RQ9e/as9LmrTdnGiMkRERERETUJ5auqZWRkSPvq6/5EH374IW7cuIHly5fj0qVL0ilsjztx4oT0c2lpKRITE9GuXTsAQLt27RAfH6+U2B07dgxGRkZ4+eWXAQA6OjooKyurdWyvvvoqUlJS0Lp16wpbeUKpr6+Pt99+G8uXL8fhw4cRHx+PCxcuVNpebco2NrzmiIiIiBo8h5m71R0CNQL6+vp47bXXsHDhQjg6OiI7OxuzZ8+ul77NzMzw3nvvYdq0aejbt6+U0Dxu5cqVcHJyQrt27bB06VLcu3dPWjjh008/RWRkJCZMmICgoCCkpKRg7ty5CA4OlhIYBwcHJCQkIC0tDYaGhjA3N69RbKGhoXjrrbfQokULvP/++9DQ0MC5c+dw8eJFLFiwADExMSgrK4O7uzsMDAzw/fffQ19fH/b29hXaqk3ZxogzR0RERETUZKxfvx6lpaVwc3PDpEmTpNXe6oO/vz+Ki4uVVop73MKFC7Fw4UJ07twZf/zxB3bu3AkLCwsAwEsvvYQ9e/bg5MmT6Ny5Mz755BP4+/srJXdTp06FpqYm2rdvD0tLS6Snp9coLh8fH+zatQv79+9H165d8dprr2Hp0qVSQmNqaopvv/0W3bt3R6dOnXDgwAH8+uuvaN68eYW2alO2MZKJx+fumoD8/HyYmJjAC+9CS6b99ApEVC3NDq+oOwSJ+N9NdYcgkenrqTsESVnOPXWHAAAoFSU4jB3Iy8uDsbGxusNpUMo/m/jcPLvGPnPU0O4f9PDhQ6SmpsLR0RF6eg3n/7PG7r///S8mT56M27dvQ0dHR9pf2f2XXiRlZWU4e/YsAMDV1RWampp12n51v8+1/f+Xp9URERERET2H+/fvIyMjAwsXLsS4ceOUEiNqXHhaHRERERHRc1i0aBHatm0LuVyOkJAQdYdDz4EzR0REREREzyEsLAxhYWFVHndwcEATu5KlyeLMEREREREREZgcERFRI3b06FG8/fbbsLW1hUwmw/bt26ss+8knn0AmkyEyMlJpf05ODnx9fWFsbAxTU1P4+/ujoKBAqcz58+fx+uuvSzdLXLRokQpGQ0RE6sbkiIiIGq3CwkJ07twZK1eurLbctm3bcOLECdja2lY45uvri+TkZMTGxmLXrl04evQoxo4dKx3Pz89H3759YW9vj8TERHz11VcICwvD2rVr63w8RESkXrzmiIiIGq3+/fujf//+1Za5desWJkyYgH379mHgQOXllC9fvoy9e/fi1KlT6NKlCwAgKioKAwYMwNdffw1bW1ts3LgRxcXFWL9+PXR0dNChQwckJSVhyZIlSkkUERE1fiqdOXJwcIBMJquwBQYGVlo+JiamQlmuvU9ERM9KoVBg5MiRmDZtGjp06FDheHx8PExNTaXECAC8vb2hoaGBhIQEqUzPnj2Vlub18fFBSkoK7t2r/B5TRUVFyM/PV9qIiKjhU+nM0alTp1BWViY9vnjxIt58800MHTq0yjrGxsZISUmRHstkMlWGSERETdiXX34JLS0tfPbZZ5Uez8zMhJWVldI+LS0tmJubIzMzUyrj6OioVMba2lo6ZmZmVqHdiIgIhIeH18UQiIioHql05sjS0hJyuVzadu3ahVatWqFXr15V1pHJZEp1yj+AiIiIaiMxMRHLli2TzkqoTyEhIcjLy5O2mzdv1mv/RNS4PW2BmacJCwuDi4tLncXzIqm3BRmKi4vx/fff46OPPqr2Q6qgoAD29vaws7PDu+++i+Tk5Grb5akLRERUmd9//x3Z2dlo0aIFtLS0oKWlhRs3bmDKlClwcHAAAMjlcmRnZyvVKy0tRU5ODuRyuVQmKytLqUz54/IyT9LV1YWxsbHSRkTPz8vLC5MmTarTNv38/DBo0CC1tFFVEpORkfHU6ynLVZZITZ06FXFxcbWOh+oxOdq+fTtyc3Ph5+dXZZk2bdpg/fr12LFjB77//nsoFAp4enrir7/+qrJOREQETExMpM3Ozk4F0RMRUWMzcuRInD9/HklJSdJma2uLadOmYd++fQAADw8P5ObmIjExUap38OBBKBQKuLu7S2WOHj2KkpISqUxsbCzatGlT6Sl1RETPSy6XQ1dX95nrGxoaonnz5nUY0Yuj3pKjdevWoX///pUuo1rOw8MDo0aNgouLC3r16oVffvkFlpaW+Oabb6qsw1MXiIheXAUFBVLiAwCpqalISkpCeno6mjdvjo4dOypt2trakMvlaNOmDQCgXbt26NevHwICAnDy5EkcO3YMQUFBGD58uPR59a9//Qs6Ojrw9/dHcnIytmzZgmXLliE4OFhdwyZ6Ifn5+eHIkSNYtmyZtHBXWloagEfXtffv3x+GhoawtrbGyJEjcffuXanuTz/9BGdnZ+jr66N58+bw9vZGYWEhwsLCsGHDBuzYsUNq8/Dhw5X2/yxtzJgxA6+88goMDAzQsmVLzJkzR/qiJSYmBuHh4Th37pxULyYmBoDybFBxcTGCgoJgY2MDPT092NvbIyIiAgCkWfDBgwdDJpNJjyubkVq/fj06dOgAXV1d2NjYICgo6Llej6aqXpbyvnHjBg4cOIBffvmlVvW0tbXh6uqKa9euVVlGV1f3uTJrIiJqvE6fPo3evXtLj8sTltGjR0t/ZDzNxo0bERQUhD59+kBDQwNDhgzB8uXLpeMmJibYv38/AgMD4ebmBgsLC4SGhnIZb6J6tmzZMvz555/o2LEj5s2bB+DR9e25ubl444038PHHH2Pp0qV48OABZsyYgWHDhuHgwYPIyMjAiBEjsGjRIgwePBj//PMPfv/9dwghMHXqVFy+fBn5+fmIjo4GAJibm1fo+1nbMDIyQkxMDGxtbXHhwgUEBATAyMgI06dPxwcffICLFy9i7969OHDgAIBH/988afny5di5cyd+/PFHtGjRAjdv3pQmA06dOgUrKytER0ejX79+0NTUrPS5W716NYKDg7Fw4UL0798feXl5OHbs2HO+Ik1TvSRH0dHRsLKyqnB/iacpKyvDhQsXMGDAABVFRkREjZmXlxeEEDUuX/4t8+PMzc2xadOmaut16tQJv//+e23DI6I6ZGJiAh0dHRgYGChd77dixQq4urriiy++kPatX78ednZ2+PPPP1FQUIDS0lK89957sLe3BwA4OztLZfX19VFUVFTlNYTAo+ToWdqYPXu29LODgwOmTp2KzZs3Y/r06dDX14ehoSG0tLSq7Ts9PR1OTk7o0aMHZDKZ1D/wKDkEAFNT02rbWLBgAaZMmYKJEydK+7p27Vpl+ReZyk+rUygUiI6OxujRo6GlpZyLjRo1CiEhIdLjefPmYf/+/fjf//6HM2fO4MMPP8SNGzfw8ccfqzpMIiIiImqEzp07h0OHDsHQ0FDa2rZtCwC4fv06OnfujD59+sDZ2RlDhw7Ft99+W+U9yqryrG1s2bIF3bt3h1wuh6GhIWbPno309PRa9e3n54ekpCS0adMGn332Gfbv31+r+tnZ2bh9+zb69OlTq3ovKpUnRwcOHEB6ejo++uijCsfS09ORkZEhPb537x4CAgLQrl07DBgwAPn5+Th+/Djat2+v6jCJiIiIqBEqKCjA22+/rbT4SlJSEq5evYqePXtCU1MTsbGx+O2339C+fXtERUWhTZs2SE1NrXEfz9JGfHw8fH19MWDAAOzatQtnz57FrFmzUFxcXKvxvfrqq0hNTcX8+fPx4MEDDBs2DO+//36N6+vr69eqvxedyk+r69u3b5WnPDx5wdvSpUuxdOlSVYdERERERI2Qjo4OysrKlPa9+uqr+Pnnn+Hg4FDhLKVyMpkM3bt3R/fu3REaGgp7e3ts27YNwcHBlbZZF20cP34c9vb2mDVrlrTvxo0bTx1PZYyNjfHBBx/ggw8+wPvvv49+/fohJycH5ubm0NbWrrYNIyMjODg4IC4uTukaTapcva1WR0RERET0PBwcHJCQkIC0tDTcvXsXCoUCgYGByMnJwYgRI3Dq1Clcv34d+/btw5gxY1BWVoaEhAR88cUXOH36NNLT0/HLL7/gzp07aNeundTm+fPnkZKSgrt37yot21/uWdpwcnJCeno6Nm/ejOvXr2P58uXYtm1bhfGUr7J59+5dFBUVVeh7yZIl+OGHH3DlyhX8+eef2Lp1K+RyOUxNTaU24uLikJmZWeWpfmFhYVi8eDGWL1+Oq1ev4syZM4iKinqel6LJYnJERERERI3C1KlToampifbt28PS0hLp6emwtbXFsWPHUFZWhr59+8LZ2RmTJk2CqakpNDQ0YGxsjKNHj2LAgAF45ZVXMHv2bCxevFi6yWpAQADatGmDLl26wNLSstJV3J6ljXfeeQeTJ09GUFAQXFxccPz4ccyZM0ep3SFDhqBfv37o3bs3LC0t8cMPP1To28jICIsWLUKXLl3QtWtXpKWlYc+ePdDQePRn/OLFixEbGws7Ozu4urpW+ryNHj0akZGRWLVqFTp06IC33noLV69efa7XoqmSidos89MI5Ofnw8TEBF54F1oybXWHQ1RrGg3t3OB2LdUdgUShWy8LbNaIVlaeukOQlP4vTd0hAABKRQkOYwfy8vJgbGys7nAalPLPJj43z85h5m51h/Bc0hbWbsVeVXv48CFSU1Ph6OgIPT09dYdDTVxZWRnOnj0LAHB1da1yyfFnVd3vc23//+XMEREREREREerpPkdERERE1HAJIXD//n219G1gYACZTKaWvomexOSIiIiI6AV3//59GBoaqqXvgoICNGvWTC19Ez2Jp9URERERERGBM0dERERELzwDAwMUFBSore+aEkJg3Lhx+Omnn3Dv3j2cPXsWLi4uqguuBhwcHDBp0iRMmjRJrXFQ3WByRERERPSCk8lkjeLUtr179yImJgaHDx9Gy5YtYWFhoe6QcOrUKbU8d35+fsjNzcX27dvrve+mjMkRERERETUK169fh42NDTw9PassU1xcDB0dnXqLydLSst76ItXjNUdERERE1OD5+flhwoQJSE9Ph0wmg4ODAwDAy8sLQUFBmDRpEiwsLODj4wMAOHLkCLp16wZdXV3Y2Nhg5syZKC0tldrz8vLChAkTMGnSJJiZmcHa2hrffvstCgsLMWbMGBgZGaF169b47bffqo3LwcEBkZGR0mOZTIb//Oc/GDx4MAwMDODk5ISdO3dKxw8fPgyZTIbdu3ejU6dO0NPTw2uvvYaLFy9KZcLCwiqcLhgZGSmNOSwsDBs2bMCOHTsgk8kgk8lw+PDh2j+pVAGTIyIiIiJq8JYtW4Z58+bh5ZdfRkZGBk6dOiUd27BhA3R0dHDs2DGsWbMGt27dwoABA9C1a1ecO3cOq1evxrp167BgwQKlNjds2AALCwucPHkSEyZMwPjx4zF06FB4enrizJkz6Nu3L0aOHFnrZc7Dw8MxbNgwnD9/HgMGDICvry9ycnKUykybNg2LFy/GqVOnYGlpibfffhslJSU1an/q1KkYNmwY+vXrh4yMDGRkZFQ7m0Y1x+SIiIiIiBo8ExMTGBkZQVNTE3K5XOl0NicnJyxatAht2rRBmzZtsGrVKtjZ2WHFihVo27YtBg0ahPDwcCxevBgKhUKq17lzZ8yePRtOTk4ICQmBnp4eLCwsEBAQACcnJ4SGhuLvv//G+fPnaxWrn58fRowYgdatW+OLL75AQUEBTp48qVRm7ty5ePPNN+Hs7IwNGzYgKysL27Ztq1H7hoaG0NfXh66uLuRyOeRyeb2eStiUMTkiIiIiokbNzc1N6fHly5fh4eGhdHPZ7t27o6CgAH/99Ze0r1OnTtLPmpqaaN68OZydnaV91tbWAIDs7OxaxfN4u82aNYOxsXGFNjw8PKSfzc3N0aZNG1y+fLlW/VDdY3JERERERI3as64Wp62trfRYJpMp7StPrh6fbXrWdmvThoaGBoQQSvtqesodPR8mR0RERETUpLRr1w7x8fFKCcaxY8dgZGSEl19+WY2R/Z8TJ05IP9+7dw9//vkn2rVrB+DRCniZmZlK8SclJSnV19HRQVlZWb3E+iJhckRERERETcqnn36KmzdvYsKECbhy5Qp27NiBuXPnIjg4GBoaDePP33nz5iEuLg4XL16En58fLCwsMGjQIACPVtK7c+cOFi1ahOvXr2PlypUVVs1zcHDA+fPnkZKSgrt373JmqY40jN8OIiIiIqI68tJLL2HPnj04efIkOnfujE8++QT+/v6YPXu2ukOTLFy4EBMnToSbmxsyMzPx66+/SosqtGvXDqtWrcLKlSvRuXNnnDx5ElOnTlWqHxAQgDZt2qBLly6wtLTEsWPH1DGMJkcmnjyhsZHLz8+HiYkJvPAutGTaT69A1MBo6OurOwRl7VqqOwKJQrfh3LdaKytP3SFISv+Xpu4QAAClogSHsQN5eXkwNjZWdzgNSvlnE5+bZ+cwc7e6Q3guaQsHqjsEJQ8fPkRqaiocHR2hp6en7nBeKIcPH0bv3r1x7949mJqaqjucelFWVoazZ88CAFxdXaGpqVmn7Vf3+1zb/385c0RERERERAQmR0RERERERACAhnOOChERERFRE+fl5VVhmW5qODhzREREREREBCZHRETUiB09ehRvv/02bG1tIZPJsH37dulYSUkJZsyYAWdnZzRr1gy2trYYNWoUbt++rdRGTk4OfH19YWxsDFNTU/j7+6OgoECpzPnz5/H6669DT08PdnZ2WLRoUX0Mj0jlOINBTUFd/h4zOSIiokarsLAQnTt3xsqVKyscu3//Ps6cOYM5c+bgzJkz+OWXX5CSkoJ33nlHqZyvry+Sk5MRGxuLXbt24ejRoxg7dqx0PD8/H3379oW9vT0SExPx1VdfISwsDGvXrlX5+IhUpXy1sOLiYjVHQvT87t+/DwDQ1n7+lap5zRERETVa/fv3R//+/Ss9ZmJigtjYWKV9K1asQLdu3ZCeno4WLVrg8uXL2Lt3L06dOoUuXboAAKKiojBgwAB8/fXXsLW1xcaNG1FcXIz169dDR0cHHTp0QFJSEpYsWaKURD2uqKgIRUVF0uP8/Pw6GjFR3dDS0oKBgQHu3LkDbW3tBnNjVGqaysrKpJ8fPnxYZ0t5CyFw//59ZGdnw9TUtE7aZXJEREQvjLy8PMhkMuneIvHx8TA1NZUSIwDw9vaGhoYGEhISMHjwYMTHx6Nnz57SzRkBwMfHB19++SXu3bsHMzOzCv1EREQgPDxc5eMhelYymQw2NjZITU3FjRs31B0ONXEKhQJ3794FAKSlpdV5Mm5qagq5XF4nbak0OQoLC6vw4dCmTRtcuXKlyjpbt27FnDlzkJaWBicnJ3z55ZcYMGCAKsMkIqIXwMOHDzFjxgyMGDFCuhFgZmYmrKyslMppaWnB3NwcmZmZUhlHR0elMtbW1tKxypKjkJAQBAcHS4/z8/NhZ2dXp+Mhel46OjpwcnLiqXWkcvfv38fAgY9uhHzmzBkYGBjUWdva2tp1elNZlc8cdejQAQcOHPi/DrWq7vL48eMYMWIEIiIi8NZbb2HTpk0YNGgQzpw5g44dO6o6VCIiaqJKSkowbNgwCCGwevVqlfenq6sLXV1dlfdD9Lw0NDSgp6en7jCoiSsrK5NmKHV1dRv075zKTzDV0tKCXC6XNgsLiyrLLlu2DP369cO0adPQrl07zJ8/H6+++ipWrFih6jCJiKiJKk+Mbty4gdjYWGnWCADkcjmys7OVypeWliInJ0c6RUMulyMrK0upTPnjujqNg4iIGgaVzxxdvXoVtra20NPTg4eHByIiItCiRYtKy8bHxyudhgA8Oq/78aVZn8SLXomIqCrlidHVq1dx6NAhNG/eXOm4h4cHcnNzkZiYCDc3NwDAwYMHoVAo4O7uLpWZNWsWSkpKpJWQYmNj0aZNm0pPqSOqjMPM3c9VP23hwDqKhIiqo9KZI3d3d8TExGDv3r1YvXo1UlNT8frrr+Off/6ptHxmZqZ0Hnc5a2tr6bzvykRERMDExETaeE43EdGLo6CgAElJSUhKSgIApKamIikpCenp6SgpKcH777+P06dPY+PGjSgrK0NmZiYyMzOlayzatWuHfv36ISAgACdPnsSxY8cQFBSE4cOHw9bWFgDwr3/9Czo6OvD390dycjK2bNmCZcuWVfgyj4iIGj+Vzhw9vrxqp06d4O7uDnt7e/z444/w9/evkz540SsR0Yvr9OnT6N27t/S4/PNg9OjRCAsLw86dOwEALi4uSvUOHToELy8vAMDGjRsRFBSEPn36QENDA0OGDMHy5culsiYmJti/fz8CAwPh5uYGCwsLhIaGVrmMNxERNV71upS3qakpXnnlFVy7dq3S41Wd113dOd286LXx0jRvQKej2FiqOwLJndeqvi5PHcwvFao7BElO+2bqDkFilZWn7hAkGvr66g4BAKAhNIEH9dunl5dXtXdGr8ld083NzbFp06Zqy3Tq1Am///57reMjIqLGpV7v+FVQUIDr16/Dxsam0uMeHh6Ii4tT2hcbGwsPD4/6CI+IiIiIiF5gKk2Opk6diiNHjiAtLQ3Hjx/H4MGDoampiREjRgAARo0ahZCQEKn8xIkTsXfvXixevBhXrlxBWFgYTp8+jaCgIFWGSUREREREpNrT6v766y+MGDECf//9NywtLdGjRw+cOHEClpaPTmFKT09XukOup6cnNm3ahNmzZ+Pf//43nJycsH37dt7jiIiIiIiIVE6lydHmzZurPX748OEK+4YOHYqhQ4eqKCIiIiJSh+ddypqIqD7U6zVHREREREREDRWTIyIiIiIiIjA5IiIiIiIiAsDkiIiIiIiICACTIyIiIiIiIgBMjoiIiIiIiAAwOSIiIiIiIgLA5IiIiIiIiAgAkyMiIiIiIiIATI6IiIiIiIgAMDkiIiIiIiICwOSIiIiIiIgIAJMjIiIiIiIiAEyOiIiIiIiIADA5IiIiIiIiAsDkiIiIiIiICACTIyIiIiIiIgBMjoiIiIiIiAAwOSIiIiIiIgLA5IiIiIiIiAgAkyMiIiIiIiIATI6IiIiIiIgAMDkiIqJG7OjRo3j77bdha2sLmUyG7du3Kx0XQiA0NBQ2NjbQ19eHt7c3rl69qlQmJycHvr6+MDY2hqmpKfz9/VFQUKBU5vz583j99dehp6cHOzs7LFq0SNVDIyIiNWByREREjVZhYSE6d+6MlStXVnp80aJFWL58OdasWYOEhAQ0a9YMPj4+ePjwoVTG19cXycnJiI2Nxa5du3D06FGMHTtWOp6fn4++ffvC3t4eiYmJ+OqrrxAWFoa1a9eqfHxERFS/tNQdABER0bPq378/+vfvX+kxIQQiIyMxe/ZsvPvuuwCA7777DtbW1ti+fTuGDx+Oy5cvY+/evTh16hS6dOkCAIiKisKAAQPw9ddfw9bWFhs3bkRxcTHWr18PHR0ddOjQAUlJSViyZIlSEvW4oqIiFBUVSY/z8/PreORERKQKTI5eMBr6+uoO4f/YWKo7Asmd1yzUHYIkx+vh0wvVoxwvTXWHIDEwzFN3CJLi6+bqDkGindswnheZQgN4oO4o/k9qaioyMzPh7e0t7TMxMYG7uzvi4+MxfPhwxMfHw9TUVEqMAMDb2xsaGhpISEjA4MGDER8fj549e0JHR0cq4+Pjgy+//BL37t2DmZlZhb4jIiIQHh6u2gESEVGd42l1RETUJGVmZgIArK2tlfZbW1tLxzIzM2FlZaV0XEtLC+bm5kplKmvj8T6eFBISgry8PGm7efPm8w+IiIhUTqXJUUREBLp27QojIyNYWVlh0KBBSElJqbZOTEwMZDKZ0qanp6fKMImIiOqUrq4ujI2NlTYiImr4VJocHTlyBIGBgThx4gRiY2NRUlKCvn37orCwsNp6xsbGyMjIkLYbN26oMkwiImqC5HI5ACArK0tpf1ZWlnRMLpcjOztb6XhpaSlycnKUylTWxuN9EBFR06DS5Gjv3r3w8/NDhw4d0LlzZ8TExCA9PR2JiYnV1pPJZJDL5dL25OkMRERET+Po6Ai5XI64uDhpX35+PhISEuDh4QEA8PDwQG5urtLn0sGDB6FQKODu7i6VOXr0KEpKSqQysbGxaNOmTaXXGxERUeNVr9cc5eU9umjY3Lz6C5kLCgpgb28POzs7vPvuu0hOTq6ybFFREfLz85U2IiJ6MRQUFCApKQlJSUkAHi3CkJSUhPT0dMhkMkyaNAkLFizAzp07ceHCBYwaNQq2trYYNGgQAKBdu3bo168fAgICcPLkSRw7dgxBQUEYPnw4bG1tAQD/+te/oKOjA39/fyQnJ2PLli1YtmwZgoOD1TRqIiJSlXpLjhQKBSZNmoTu3bujY8eOVZZr06YN1q9fjx07duD777+HQqGAp6cn/vrrr0rLR0REwMTERNrs7OxUNQQiImpgTp8+DVdXV7i6ugIAgoOD4erqitDQUADA9OnTMWHCBIwdOxZdu3ZFQUEB9u7dq3Qt68aNG9G2bVv06dMHAwYMQI8ePZTuYWRiYoL9+/cjNTUVbm5umDJlCkJDQ6tcxpuIiBovmRBC1EdH48ePx2+//YY//vgDL7/8co3rlZSUoF27dhgxYgTmz59f4Xhl95Kws7ODF96Flky7TmJvShrSUt6ylg0nkeVS3o2DgWHR0wvVE/kyXXWHINE+n6ruEAAApYpixN3bgLy8PC5A8IT8/HyYmJi80M+Nw8zd6g6hUUtbOFDdIRA9s8LCQhgaGgJ4NOPfrFmzeuu7tv//1st9joKCgqS7jtcmMQIAbW1tuLq64tq1a5Ue19XVha5uw/kjhYiIiIiIGieVnlYnhEBQUBC2bduGgwcPwtHRsdZtlJWV4cKFC7CxsVFBhERERERERI+odOYoMDAQmzZtwo4dO2BkZCTdLM/ExAT6///0rlGjRuGll15CREQEAGDevHl47bXX0Lp1a+Tm5uKrr77CjRs38PHHH6syVCIiIiIiesGpNDlavXo1AMDLy0tpf3R0NPz8/AAA6enp0ND4vwmse/fuISAgAJmZmTAzM4ObmxuOHz+O9u3bqzJUIiIiIiJ6wak0OarJWg+HDx9Werx06VIsXbpURRERERERERFVrl7vc0RERERERNRQMTkiIiIiIiICkyMiIiIiIiIATI6IiIiIiIgAMDkiIiIiIiICwOSIiIiIiIgIAJMjIiIiIiIiAEyOiIiIiIiIADA5IiIiIiIiAsDkiIiIiIiICACTIyIiIiIiIgBMjoiIiIiIiAAwOSIiIiIiIgLA5IiIiIiIiAgAkyMiIiIiIiIATI6IiIiIiIgAMDkiIiIiIiICAGipOwAiIiIiqp7DzN3PVT9t4cA6ioSoaWNypGIa+vrqDkGJTF9P3SFIHtiZqDsEidXINHWHIGmpVazuEJRsaRmn7hAkrX4ap+4QJBpFD9QdgqQs5566QwAAlIkSdYdARET0XHhaHREREREREZgcERERERERAWByRERETVhZWRnmzJkDR0dH6Ovro1WrVpg/fz6EEFIZIQRCQ0NhY2MDfX19eHt74+rVq0rt5OTkwNfXF8bGxjA1NYW/vz8KCgrqezhERKRiTI6IiKjJ+vLLL7F69WqsWLECly9fxpdffolFixYhKipKKrNo0SIsX74ca9asQUJCApo1awYfHx88fPhQKuPr64vk5GTExsZi165dOHr0KMaOHauOIRERkQpxQQYiImqyjh8/jnfffRcDBz5aqcvBwQE//PADTp48CeDRrFFkZCRmz56Nd999FwDw3XffwdraGtu3b8fw4cNx+fJl7N27F6dOnUKXLl0AAFFRURgwYAC+/vpr2NraVui3qKgIRUVF0uP8/HxVD5WIiOoAZ46IiKjJ8vT0RFxcHP78808AwLlz5/DHH3+gf//+AIDU1FRkZmbC29tbqmNiYgJ3d3fEx8cDAOLj42FqaiolRgDg7e0NDQ0NJCQkVNpvREQETExMpM3Ozk5VQyQiojrEmSMiImqyZs6cifz8fLRt2xaampooKyvD559/Dl9fXwBAZmYmAMDa2lqpnrW1tXQsMzMTVlZWSse1tLRgbm4ulXlSSEgIgoODpcf5+flMkIiIGgEmR0RE1GT9+OOP2LhxIzZt2oQOHTogKSkJkyZNgq2tLUaPHq2yfnV1daGrq6uy9omISDXq5bS6lStXwsHBAXp6enB3d5fO9a7K1q1b0bZtW+jp6cHZ2Rl79uypjzCJiKiJmTZtGmbOnInhw4fD2dkZI0eOxOTJkxEREQEAkMvlAICsrCylellZWdIxuVyO7OxspeOlpaXIycmRyhARUdOg8uRoy5YtCA4Oxty5c3HmzBl07twZPj4+FT5oyh0/fhwjRoyAv78/zp49i0GDBmHQoEG4ePGiqkMlIqIm5v79+9DQUP6o09TUhEKhAAA4OjpCLpcjLi5OOp6fn4+EhAR4eHgAADw8PJCbm4vExESpzMGDB6FQKODu7l4PoyAiovqi8uRoyZIlCAgIwJgxY9C+fXusWbMGBgYGWL9+faXlly1bhn79+mHatGlo164d5s+fj1dffRUrVqxQdahERNTEvP322/j888+xe/dupKWlYdu2bViyZAkGDx4MAJDJZJg0aRIWLFiAnTt34sKFCxg1ahRsbW0xaNAgAEC7du3Qr18/BAQE4OTJkzh27BiCgoIwfPjwSleqIyKixkul1xwVFxcjMTERISEh0j4NDQ14e3tLqwA9KT4+XukiVgDw8fHB9u3bKy3P5VKJiKgqUVFRmDNnDj799FNkZ2fD1tYW48aNQ2hoqFRm+vTpKCwsxNixY5Gbm4sePXpg79690NPTk8ps3LgRQUFB6NOnDzQ0NDBkyBAsX75cHUMiIiIVUmlydPfuXZSVlVW6CtCVK1cqrZOZmVntqkFPioiIQHh4eN0ETERETYqRkREiIyMRGRlZZRmZTIZ58+Zh3rx5VZYxNzfHpk2bVBAhERE1JI3+PkchISHIy8uTtps3b6o7JCIiIiIiaoRUOnNkYWEBTU3NalcBepJcLq9VeS6XSkREREREdUGlM0c6Ojpwc3NTWgVIoVAgLi5OWgXoSR4eHkrlASA2NrbK8kRERERERHVB5TeBDQ4OxujRo9GlSxd069YNkZGRKCwsxJgxYwAAo0aNwksvvSTdc2LixIno1asXFi9ejIEDB2Lz5s04ffo01q5dq+pQiYiIiIjoBaby5OiDDz7AnTt3EBoaiszMTLi4uGDv3r3Sogvp6elK96Dw9PTEpk2bMHv2bPz73/+Gk5MTtm/fjo4dO6o6VCIiIiIieoGpPDkCgKCgIAQFBVV67PDhwxX2DR06FEOHDlVxVERERERERP+n0a9WR0REREREVBeYHBEREREREYHJEREREREREQAmR0RERERERACYHBEREREREQFgckRERERERASAyREREREREREAJkdEREREREQAmBwREREREREBYHJEREREREQEgMkRERERERERAEBL3QE0dYoHD9QdghKNdi3VHYLkoUXD+fW7e89M3SFIvFukqDsEJW3/GKnuECQvxwl1hyDRSPpT3SFIFOoOgIiIqIngzBERERERERGYHBEREREREQFgckRERERERASAyREREREREREAJkdEREREREQAmBwREREREREB4FLeREREVAMOM3erOwQiIpXjzBERERERERGYHBERURN369YtfPjhh2jevDn09fXh7OyM06dPS8eFEAgNDYWNjQ309fXh7e2Nq1evKrWRk5MDX19fGBsbw9TUFP7+/igoKKjvoRARkYoxOSIioibr3r176N69O7S1tfHbb7/h0qVLWLx4MczMzKQyixYtwvLly7FmzRokJCSgWbNm8PHxwcOHD6Uyvr6+SE5ORmxsLHbt2oWjR49i7Nix6hgSERGpEK85IiKiJuvLL7+EnZ0doqOjpX2Ojo7Sz0IIREZGYvbs2Xj33XcBAN999x2sra2xfft2DB8+HJcvX8bevXtx6tQpdOnSBQAQFRWFAQMG4Ouvv4atrW2FfouKilBUVCQ9zs/PV9UQiYioDnHmiIiImqydO3eiS5cuGDp0KKysrODq6opvv/1WOp6amorMzEx4e3tL+0xMTODu7o74+HgAQHx8PExNTaXECAC8vb2hoaGBhISESvuNiIiAiYmJtNnZ2alohEREVJeYHBERUZP1v//9D6tXr4aTkxP27duH8ePH47PPPsOGDRsAAJmZmQAAa2trpXrW1tbSsczMTFhZWSkd19LSgrm5uVTmSSEhIcjLy5O2mzdv1vXQiIhIBXhaHRERNVkKhQJdunTBF198AQBwdXXFxYsXsWbNGowePVpl/erq6kJXV1dl7RMRkWpw5oiIiJosGxsbtG/fXmlfu3btkJ6eDgCQy+UAgKysLKUyWVlZ0jG5XI7s7Gyl46WlpcjJyZHKEBFR06CS5CgtLQ3+/v5wdHSEvr4+WrVqhblz56K4uLjael5eXpDJZErbJ598oooQiYjoBdC9e3ekpKQo7fvzzz9hb28P4NHiDHK5HHFxcdLx/Px8JCQkwMPDAwDg4eGB3NxcJCYmSmUOHjwIhUIBd3f3ehgFERHVF5WcVnflyhUoFAp88803aN26NS5evIiAgAAUFhbi66+/rrZuQEAA5s2bJz02MDBQRYhERPQCmDx5Mjw9PfHFF19g2LBhOHnyJNauXYu1a9cCAGQyGSZNmoQFCxbAyckJjo6OmDNnDmxtbTFo0CAAj2aa+vXrh4CAAKxZswYlJSUICgrC8OHDK12pjoiIGi+VJEf9+vVDv379pMctW7ZESkoKVq9e/dTkyMDAgKcpEBFRnejatSu2bduGkJAQzJs3D46OjoiMjISvr69UZvr06SgsLMTYsWORm5uLHj16YO/evdDT05PKbNy4EUFBQejTpw80NDQwZMgQLF++XB1DIiIiFaq3BRny8vJgbm7+1HIbN27E999/D7lcjrfffhtz5sypdvaI95IgIqLqvPXWW3jrrbeqPC6TyTBv3jylsxaeZG5ujk2bNqkiPKJ64TBz93PVT1s4sI4iIWrY6iU5unbtGqKiop46a/Svf/0L9vb2sLW1xfnz5zFjxgykpKTgl19+qbJOREQEwsPD6zpkIiIiIiJ6wdRqQYaZM2dWWDDhye3KlStKdW7duoV+/fph6NChCAgIqLb9sWPHwsfHB87OzvD19cV3332Hbdu24fr161XW4b0kiIiIiIioLtRq5mjKlCnw8/OrtkzLli2ln2/fvo3evXvD09NTuvi1NspXAbp27RpatWpVaRneS4KIiIiIiOpCrZIjS0tLWFpa1qjsrVu30Lt3b7i5uSE6OhoaGrVfNTwpKQnAo/tUEBERERERqZJK7nN069YteHl5oUWLFvj6669x584dZGZmIjMzU6lM27ZtcfLkSQDA9evXMX/+fCQmJiItLQ07d+7EqFGj0LNnT3Tq1EkVYRIREREREUlUsiBDbGwsrl27hmvXruHll19WOiaEAACUlJQgJSUF9+/fBwDo6OjgwIEDiIyMRGFhIezs7DBkyBDMnj1bFSESEREREREpUUly5Ofn99RrkxwcHKRECQDs7Oxw5MgRVYRDRERERET0VCo5rY6IiIiIiKixYXJEREREREQEJkdEREREREQAmBwREREREREBUNGCDNRwaeQWqjuExxirOwCJ5mETdYcg2W/YTd0hKGm5/a66Q/g/GTfVHYGk7MEDdYdAREREdYwzR0RERERERGByREREREREBIDJEREREREREQAmR0RERERERACYHBEREREREQFgckRERERERASAyREREREREREAJkdEREREREQAmBwREREREREBYHJEREREREQEgMkRERERERERACZHREREREREAJgcERERERERAWByREREREREBIDJERERvUAWLlwImUyGSZMmSfsePnyIwMBANG/eHIaGhhgyZAiysrKU6qWnp2PgwIEwMDCAlZUVpk2bhtLS0nqOnoiIVI3JERERvRBOnTqFb775Bp06dVLaP3nyZPz666/YunUrjhw5gtu3b+O9996TjpeVlWHgwIEoLi7G8ePHsWHDBsTExCA0NLS+h0BERCrG5IiIiJq8goIC+Pr64ttvv4WZmZm0Py8vD+vWrcOSJUvwxhtvwM3NDdHR0Th+/DhOnDgBANi/fz8uXbqE77//Hi4uLujfvz/mz5+PlStXori4uNL+ioqKkJ+fr7QREVHDx+SIiIiavMDAQAwcOBDe3t5K+xMTE1FSUqK0v23btmjRogXi4+MBAPHx8XB2doa1tbVUxsfHB/n5+UhOTq60v4iICJiYmEibnZ2dCkZFRER1jckRERE1aZs3b8aZM2cQERFR4VhmZiZ0dHRgamqqtN/a2hqZmZlSmccTo/Lj5ccqExISgry8PGm7efNmHYyEiIhUTUvdARAREanKzZs3MXHiRMTGxkJPT6/e+tXV1YWurm699UdERHWDM0dERNRkJSYmIjs7G6+++iq0tLSgpaWFI0eOYPny5dDS0oK1tTWKi4uRm5urVC8rKwtyuRwAIJfLK6xeV/64vAwRETUNTI6IiKjJ6tOnDy5cuICkpCRp69KlC3x9faWftbW1ERcXJ9VJSUlBeno6PDw8AAAeHh64cOECsrOzpTKxsbEwNjZG+/bt631MRESkOipLjhwcHCCTyZS2hQsXVlunJveaICIiqikjIyN07NhRaWvWrBmaN2+Ojh07wsTEBP7+/ggODsahQ4eQmJiIMWPGwMPDA6+99hoAoG/fvmjfvj1GjhyJc+fOYd++fZg9ezYCAwN56hwRUROj0muO5s2bh4CAAOmxkZFRteUnT56M3bt3Y+vWrTAxMUFQUBDee+89HDt2TJVhEhHRC2zp0qXQ0NDAkCFDUFRUBB8fH6xatUo6rqmpiV27dmH8+PHw8PBAs2bNMHr0aMybN0+NURMRkSqoNDkyMjKq8fnY5fea2LRpE9544w0AQHR0NNq1a4cTJ05I3+ARERE9j8OHDys91tPTw8qVK7Fy5coq69jb22PPnj0qjoyIiNRNpdccLVy4EM2bN4erqyu++uorlJaWVlm2JveaqAxvtEdERERERHVBZTNHn332GV599VWYm5vj+PHjCAkJQUZGBpYsWVJp+Zrca6IyERERCA8Pr8vQiYiIiIjoBVSrmaOZM2dWWGThye3KlSsAgODgYHh5eaFTp0745JNPsHjxYkRFRaGoqKhOB8Ab7RERERERUV2o1czRlClT4OfnV22Zli1bVrrf3d0dpaWlSEtLQ5s2bSocl8vl0r0mHp89evxeE5XhjfaIiIiIVMth5u7nqp+2cGAdRUKkWrVKjiwtLWFpaflMHSUlJUFDQwNWVlaVHndzc5PuNTFkyBAAFe81QUREREREpCoqueYoPj4eCQkJ6N27N4yMjBAfH4/Jkyfjww8/hJmZGQDg1q1b6NOnD7777jt069ZN6V4T5ubmMDY2xoQJE5TuNdEYaZqbqTsEJYqMhnPfKPODdXuKJamGePhQ3SFIynLuqTsEIiIiasJUkhzp6upi8+bNCAsLQ1FRERwdHTF58mQEBwdLZUpKSpCSkoL79+9L+552rwkiIiIiIiJVUUly9Oqrr+LEiRPVlnFwcIAQQmlfTe41QUREREREpAoqvc8RERERERFRY8HkiIiIiIiICEyOiIiIiIiIADA5IiIiIiIiAsDkiIiIiIiICACTIyIiIiIiIgBMjoiIiIiIiAAwOSIiIiIiIgLA5IiIiIiIiAgAkyMiIiIiIiIATI6IiIiIiIgAAFrqDoCIiIiImjaHmbufq37awoF1FAlR9ThzREREREREBCZHREREREREAJgcERERERERAWByREREREREBIDJEREREREREQAmR0RERERERACYHBERURMWERGBrl27wsjICFZWVhg0aBBSUlKUyjx8+BCBgYFo3rw5DA0NMWTIEGRlZSmVSU9Px8CBA2FgYAArKytMmzYNpaWl9TkUIiKqB0yOiIioyTpy5AgCAwNx4sQJxMbGoqSkBH379kVhYaFUZvLkyfj111+xdetWHDlyBLdv38Z7770nHS8rK8PAgQNRXFyM48ePY8OGDYiJiUFoaKg6hkRERCrEm8ASEVGTtXfvXqXHMTExsLKyQmJiInr27Im8vDysW7cOmzZtwhtvvAEAiI6ORrt27XDixAm89tpr2L9/Py5duoQDBw7A2toaLi4umD9/PmbMmIGwsDDo6OhU6LeoqAhFRUXS4/z8fNUOlIiI6gRnjoiI6IWRl5cHADA3NwcAJCYmoqSkBN7e3lKZtm3bokWLFoiPjwcAxMfHw9nZGdbW1lIZHx8f5OfnIzk5udJ+IiIiYGJiIm12dnaqGhIREdUhJkdERPRCUCgUmDRpErp3746OHTsCADIzM6GjowNTU1OlstbW1sjMzJTKPJ4YlR8vP1aZkJAQ5OXlSdvNmzfreDRERKQKPK2OiIheCIGBgbh48SL++OMPlfelq6sLXV1dlfdDRER1i8kRERE1eUFBQdi1axeOHj2Kl19+Wdovl8tRXFyM3NxcpdmjrKwsyOVyqczJkyeV2itfza68TGPgMHO3ukMgImrweFodERE1WUIIBAUFYdu2bTh48CAcHR2Vjru5uUFbWxtxcXHSvpSUFKSnp8PDwwMA4OHhgQsXLiA7O1sqExsbC2NjY7Rv375+BkJERPWCM0dERNRkBQYGYtOmTdixYweMjIyka4RMTEygr68PExMT+Pv7Izg4GObm5jA2NsaECRPg4eGB1157DQDQt29ftG/fHiNHjsSiRYuQmZmJ2bNnIzAwkKfOERE1MSpJjg4fPozevXtXeuzkyZPo2rVrpce8vLxw5MgRpX3jxo3DmjVrah2Dhr4eNGQVl1etb2U599QdghINfX11hyApvZ2h7hAapIb0GhE1dqtXrwbw6PPlcdHR0fDz8wMALF26FBoaGhgyZAiKiorg4+ODVatWSWU1NTWxa9cujB8/Hh4eHmjWrBlGjx6NefPm1dcwiIionqgkOfL09ERGhvIfvnPmzEFcXBy6dOlSbd2AgAClDxwDAwNVhEhERC8AIcRTy+jp6WHlypVYuXJllWXs7e2xZ8+eugyNiIgaIJUkRzo6OkoXqZaUlGDHjh2YMGECZDJZtXUNDAwa1QWuRERERETUNNTLggw7d+7E33//jTFjxjy17MaNG2FhYYGOHTsiJCQE9+/fr7Z8UVER8vPzlTYiIiIiIqLaqpcFGdatWwcfHx+l5VMr869//Qv29vawtbXF+fPnMWPGDKSkpOCXX36psk5ERATCw8PrOmQiIiIiInrB1GrmaObMmZDJZNVuV65cUarz119/Yd++ffD3939q+2PHjoWPjw+cnZ3h6+uL7777Dtu2bcP169errMO7kBMRERERUV2o1czRlClTpNV9qtKyZUulx9HR0WjevDneeeedWgfn7u4OALh27RpatWpVaRnehZyIiIioaXvemxinLRxYR5FQU1er5MjS0hKWlpY1Li+EQHR0NEaNGgVtbe1aB5eUlAQAsLGxqXVdIiIiIiKi2lDpggwHDx5EamoqPv744wrHbt26hbZt2+LkyZMAgOvXr2P+/PlITExEWloadu7ciVGjRqFnz57o1KmTKsMkIiIiIiJS7YIM69atg6enJ9q2bVvhWElJCVJSUqTV6HR0dHDgwAFERkaisLAQdnZ2GDJkCGbPnq3KEImIiIiIiACoODnatGlTlcccHByUbs5nZ2eHI0eOqDIcIiIiIiKiKtXLfY6IiIiIiIgaOiZHREREREREYHJEREREREQEgMkRERERERERABUvyEBEREREpG68iSzVFGeOiIiIiIiIwOSIiIiIiIgIAJMjIiIiIiIiAEyOiIiIiIiIADA5IiIiIiIiAsDkiIiIiIiICACX8iYiIiIiqhaXAn9xcOaIiIiIiIgITI6IiIiIiIgA8LQ6IiKiRuF5T+shIqKn48wREREREREROHNERERUIytXrsRXX32FzMxMdO7cGVFRUejWrZu6wyKiRoALOjQenDkiIiJ6ii1btiA4OBhz587FmTNn0LlzZ/j4+CA7O1vdoRERUR3izBEREdFTLFmyBAEBARgzZgwAYM2aNdi9ezfWr1+PmTNnqjk6Imrq1H3N4Ys0c9XkkiMhBACgVJSoOZJHFA0kjnIaQlPdIUga2nPTUDSk16ih4e9Mw1aKR69P+f/DTUVxcTESExMREhIi7dPQ0IC3tzfi4+MrrVNUVISioiLpcV5eHgAgPz//meNQFN1/5rpERM/jef7vAoDCwkKltsrKyp43pBorj72mn01NLjn6559/AABHH/6i5kgaqAfqDoCeiq8RNXL//PMPTExM1B1Gnbl79y7KyspgbW2ttN/a2hpXrlyptE5ERATCw8Mr7Lezs1NJjEREqmQSWXdt2dra1l1jtVDTz6YmlxzZ2tri5s2bMDIygkwmU3c4tZKfnw87OzvcvHkTxsbG6g6nzjTFcXFMjUdTHFdDHZMQAv/884/aPvgakpCQEAQHB0uPFQoFcnJy0Lx5c8hksgb7Gj4vjqvxaapj47gaH1WNrbafTU0uOdLQ0MDLL7+s7jCei7GxcZP7hQea5rg4psajKY6rIY6pKc0YlbOwsICmpiaysrKU9mdlZUEul1daR1dXF7q6ukr7TE1NK5RriK9hXeC4Gp+mOjaOq/FRxdhq89nE1eqIiIiqoaOjAzc3N8TFxUn7FAoF4uLi4OHhocbIiIiorjW5mSMiIqK6FhwcjNGjR6NLly7o1q0bIiMjUVhYKK1eR0RETQOTowZEV1cXc+fOrXAqRmPXFMfFMTUeTXFcTXFMDd0HH3yAO3fuIDQ0FJmZmXBxccHevXsrLNJQU031NeS4Gp+mOjaOq/FpKGOTiaa25ioREREREdEz4DVHREREREREYHJEREREREQEgMkRERERERERACZHREREREREAJgcERERERERAWBy1ODt3r0b7u7u0NfXh5mZGQYNGqTukOpEUVERXFxcIJPJkJSUpO5wnllaWhr8/f3h6OgIfX19tGrVCnPnzkVxcbG6Q6u1lStXwsHBAXp6enB3d8fJkyfVHdIzi4iIQNeuXWFkZAQrKysMGjQIKSkp6g6rTi1cuBAymQyTJk1SdygvhNq+P7Zu3Yq2bdtCT08Pzs7O2LNnj9JxIQRCQ0NhY2MDfX19eHt74+rVq0plcnJy4OvrC2NjY5iamsLf3x8FBQVKZc6fP4/XX38denp6sLOzw6JFixr9uNLS0iCTySpsJ06caNDj+vzzz+Hp6QkDAwOYmppW2k96ejoGDhwIAwMDWFlZYdq0aSgtLa3xuBry2Cp7zTZv3txgx1XTz+/G9h6rybjq4j2mjrEBwDvvvIMWLVpAT08PNjY2GDlyJG7fvq1U5nlfMwhqsH766SdhZmYmVq9eLVJSUkRycrLYsmWLusOqE5999pno37+/ACDOnj2r7nCe2W+//Sb8/PzEvn37xPXr18WOHTuElZWVmDJlirpDq5XNmzcLHR0dsX79epGcnCwCAgKEqampyMrKUndoz8THx0dER0eLixcviqSkJDFgwADRokULUVBQoO7Q6sTJkyeFg4OD6NSpk5g4caK6w2nyavv+OHbsmNDU1BSLFi0Sly5dErNnzxba2triwoULUpmFCxcKExMTsX37dnHu3DnxzjvvCEdHR/HgwQOpTL9+/UTnzp3FiRMnxO+//y5at24tRowYIR3Py8sT1tbWwtfXV1y8eFH88MMPQl9fX3zzzTeNelypqakCgDhw4IDIyMiQtuLi4gY9rtDQULFkyRIRHBwsTExMKvRTWloqOnbsKLy9vcXZs2fFnj17hIWFhQgJCanRuBry2IQQAoCIjo5Wes0eb6Ohjasmn9+N8T1Wk3E973tMXWMTQoglS5aI+Ph4kZaWJo4dOyY8PDyEh4eHdPx5XzMhhGBy1ECVlJSIl156SfznP/9Rdyh1bs+ePaJt27YiOTm50SdHlVm0aJFwdHRUdxi10q1bNxEYGCg9LisrE7a2tiIiIkKNUdWd7OxsAUAcOXJE3aE8t3/++Uc4OTmJ2NhY0atXLyZH9aC2749hw4aJgQMHKu1zd3cX48aNE0IIoVAohFwuF1999ZV0PDc3V+jq6ooffvhBCCHEpUuXBABx6tQpqcxvv/0mZDKZuHXrlhBCiFWrVgkzMzNRVFQklZkxY4Zo06ZNox5X+R9uz/rZoI5xPS46OrrSBGLPnj1CQ0NDZGZmSvtWr14tjI2NlV7Dxjg2IR4lR9u2bavROJ6k7nGVe/LzuzG+x2oyrud9jwnRcMa2Y8cOIZPJpMTueV8zIYTgaXUN1JkzZ3Dr1i1oaGjA1dUVNjY26N+/Py5evKju0J5LVlYWAgIC8N///hcGBgbqDkcl8vLyYG5uru4waqy4uBiJiYnw9vaW9mloaMDb2xvx8fFqjKzu5OXlAUCjel2qEhgYiIEDByq9XqQ6z/L+iI+Pr/D6+Pj4SOVTU1ORmZmpVMbExATu7u5Smfj4eJiamqJLly5SGW9vb2hoaCAhIUEq07NnT+jo6Cj1k5KSgnv37jXacZV75513YGVlhR49emDnzp3Vjkfd46qJ+Ph4ODs7w9raWqmf/Px8JCcnN+qxlQsMDISFhQW6deuG9evXQwjRqMb15Od3Y3yP1WRc5Z7lPdaQxpaTk4ONGzfC09MT2traUj/P+ppJY6lRKap3//vf/wAAYWFhmD17Nnbt2gUzMzN4eXkhJydHzdE9GyEE/Pz88Mknnyh9MDYl165dQ1RUFMaNG6fuUGrs7t27KCsrU/rABgBra2tkZmaqKaq6o1AoMGnSJHTv3h0dO3ZUdzjPZfPmzThz5gwiIiLUHcoL41neH5mZmdWWL//3aWWsrKyUjmtpacHc3FypTGVtPN5HYxyXoaEhFi9ejK1bt2L37t3o0aMHBg0aVKM/3tQ1rpp4ntcLaNhjA4B58+bhxx9/RGxsLIYMGYJPP/0UUVFRT63XUMZV2ed3Y3yP1WRcz/MeawhjmzFjBpo1a4bmzZsjPT0dO3bseGo/j/fxNEyO6tnMmTMrvQju8e3KlStQKBQAgFmzZmHIkCFwc3NDdHQ0ZDIZtm7dquZRKKvpmKKiovDPP/8gJCRE3SE/VU3H9Lhbt26hX79+GDp0KAICAtQUOT0pMDAQFy9erNWFwQ3RzZs3MXHiRGzcuBF6enrqDodIZSwsLBAcHAx3d3d07doVCxcuxIcffoivvvpK3aFRNebMmYPu3bvD1dUVM2bMwPTp0xvNa9ZUP7+rGldjf49NmzYNZ8+exf79+6GpqYlRo0bVaJayprTqrCWqkSlTpsDPz6/aMi1btkRGRgYAoH379tJ+XV1dtGzZEunp6aoMsdZqOqaDBw8iPj4eurq6Sse6dOkCX19fbNiwQYVR1k5Nx1Tu9u3b6N27Nzw9PbF27VoVR1e3LCwsoKmpiaysLKX9WVlZkMvlaoqqbgQFBWHXrl04evQoXn75ZXWH81wSExORnZ2NV199VdpXVlaGo0ePYsWKFSgqKoKmpqYaI2yanuX9IZfLqy1f/m9WVhZsbGyUyri4uEhlsrOzldooLS1FTk6OUjuV9fN4H41xXJVxd3dHbGxstWNS57hqQi6XV1jNq6avF9Cwx1YZd3d3zJ8/H0VFRRU+9x+n7nFV9/ndGN9jNRlXZWr6HgPUPzYLCwtYWFjglVdeQbt27WBnZ4cTJ07Aw8PjuV6zcpw5qmeWlpZo27ZttZuOjg7c3Nygq6urtPxwSUkJ0tLSYG9vr8YRVFTTMS1fvhznzp1DUlISkpKSpCUct2zZgs8//1zNo1BW0zEBj76Z8fLykmb3NDQa19uq/PctLi5O2qdQKBAXFwcPDw81RvbshBAICgrCtm3bcPDgQTg6Oqo7pOfWp08fXLhwQXr/JCUlSV8sJCUlMTFSkWd5f3h4eCiVB4DY2FipvKOjI+RyuVKZ/Px8JCQkSGU8PDyQm5uLxMREqczBgwehUCjg7u4ulTl69ChKSkqU+mnTpg3MzMwa7bgqk5SUpPQHU0MbV014eHjgwoULSslhbGwsjI2Nlb4IbYxjq0xSUhLMzMyqTYwA9Y7raZ/fjfE9VpNxVaam7zF1j+1J5WdaFRUVSf0862smqfHSDVTvJk6cKF566SWxb98+ceXKFeHv7y+srKxETk6OukOrE3WxWoq6/fXXX6J169aiT58+4q+//lJaErMx2bx5s9DV1RUxMTHi0qVLYuzYscLU1FRpVaXGZPz48cLExEQcPnxY6TW5f/++ukOrU1ytrn487f0xcuRIMXPmTKn8sWPHhJaWlvj666/F5cuXxdy5cytdstbU1FTs2LFDnD9/Xrz77ruVLnnt6uoqEhISxB9//CGcnJyUlrzOzc0V1tbWYuTIkeLixYti8+bNwsDAoFbLDDfEccXExIhNmzaJy5cvi8uXL4vPP/9caGhoiPXr1zfocd24cUOcPXtWhIeHC0NDQ3H27Flx9uxZ8c8//wgh/m8p7759+4qkpCSxd+9eYWlpWeulvBvi2Hbu3Cm+/fZbceHCBXH16lWxatUqYWBgIEJDQxvsuGry+d0Y32M1GdfzvsfUNbYTJ06IqKgocfbsWZGWlibi4uKEp6enaNWqlXj48GGdvGZCcCnvBq24uFhMmTJFWFlZCSMjI+Ht7S0uXryo7rDqTFNIjqKjowWASrfGJioqSrRo0ULo6OiIbt26iRMnTqg7pGdW1WsSHR2t7tDqFJOj+lPd+6NXr15i9OjRSuV//PFH8corrwgdHR3RoUMHsXv3bqXjCoVCzJkzR1hbWwtdXV3Rp08fkZKSolTm77//FiNGjBCGhobC2NhYjBkzRvpjtNy5c+dEjx49hK6urnjppZfEwoULG/24YmJiRLt27YSBgYEwNjYW3bp1E1u3bm3w4xo9enSl/+8cOnRIKpOWlib69+8v9PX1hYWFhZgyZYooKSlp9GP77bffhIuLizA0NBTNmjUTnTt3FmvWrBFlZWUNdlw1/fxubO+xmoyrLt5j6hjb+fPnRe/evYW5ubnQ1dUVDg4O4pNPPhF//fWXUjvP+5rJhKjDK5iIiIiIiIgaqcZ1cQQREREREZGKMDkiIiIiIiICkyMiIiIiIiIATI6IiIiIiIgAMDkiIiIiIiICwOSIiIiIiIgIAJMjIiIiIiIiAEyOiIiIiIiIADA5IiIiIiIiAsDkiIiIiIiICACTIyIiIiIiIgDA/wOpmeDUauApXgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "for label, (x, y) in d.items():\n", " # input is a histogram\n", " w, xe, ye = np.histogram2d(x, y)\n", "\n", " # apply USP test\n", " r = perm.usp(w, random_state=1)\n", "\n", " fig, ax = plt.subplots(1, 2, figsize=(10, 4))\n", " plt.sca(ax[0])\n", " plt.pcolormesh(xe, ye, w.T)\n", " plt.sca(ax[1])\n", " plt.hist(r.samples, bins=20, label=\"test statistic under\\nnull hypothesis\")\n", " plt.axvline(r.statistic, color=\"k\", label=\"test statistic\\nfrom input\")\n", " plt.suptitle(f\"{label}: p-value={r.pvalue:.3f}\")\n", " plt.legend()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.9" } }, "nbformat": 4, "nbformat_minor": 5 }