diff --git a/notebooks/replogle_map_building.ipynb b/notebooks/replogle_map_building.ipynb index 384ba80..61be4fa 100644 --- a/notebooks/replogle_map_building.ipynb +++ b/notebooks/replogle_map_building.ipynb @@ -22,10 +22,30 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "id": "35b44634", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/Users/safiye.celik/.pyenv/versions/3.11.5/envs/eben/lib/python3.11/site-packages/scvi/_settings.py:63: UserWarning: Since v1.0.0, scvi-tools no longer uses a random seed by default. Run `scvi.settings.seed = 0` to reproduce results from previous versions.\n", + " self.seed = seed\n", + "/Users/safiye.celik/.pyenv/versions/3.11.5/envs/eben/lib/python3.11/site-packages/scvi/_settings.py:70: UserWarning: Setting `dl_pin_memory_gpu_training` is deprecated in v1.0 and will be removed in v1.1. Please pass in `pin_memory` to the data loaders instead.\n", + " self.dl_pin_memory_gpu_training = (\n", + "/Users/safiye.celik/.pyenv/versions/3.11.5/envs/eben/lib/python3.11/site-packages/tqdm/auto.py:21: TqdmWarning: IProgress not found. Please update jupyter and ipywidgets. See https://ipywidgets.readthedocs.io/en/stable/user_install.html\n", + " from .autonotebook import tqdm as notebook_tqdm\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[(0.01, 0.99), (0.02, 0.98), (0.03, 0.97), (0.04, 0.96), (0.05, 0.95), (0.06, 0.94), (0.07, 0.93), (0.08, 0.92), (0.09, 0.91), (0.1, 0.9)]\n" + ] + } + ], "source": [ "from efaar_benchmarking.data_loading import load_replogle\n", "from efaar_benchmarking.efaar import embed_by_scvi, embed_by_pca, align_by_centering, aggregate_by_mean\n", @@ -57,17 +77,35 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "id": "3e9e8847", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "2057 perturbations in the map.\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAABckAAAEhCAYAAACk8SgcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAADDi0lEQVR4nOzdeVxU1fsH8M8wrAquCCq4IiogKEKumPtGuGtuuaGWqWhlpbiimbiU+5YKhWmuueWeaaZp+dVSUXNJLcQV11R25v7+OL8ZGdlmYPb5vF8vXzB3Zs4cSJ/ufe45zyOTJEkCEREREREREREREZEVsjH2BIiIiIiIiIiIiIiIjIVJciIiIiIiIiIiIiKyWkySExEREREREREREZHVYpKciIiIiIiIiIiIiKwWk+REREREREREREREZLWYJCciIiIiIiIiIiIiq8UkORERERERERERERFZLSbJiYiIiIiIiIiIiMhqMUlOJk2SJGNPgYiIiIiIiIiITAxzRqRLTJITBgwYgFq1aqn9qV27NurXr4/u3btj586dBp9Teno6Zs2ahR9++MHgn01E5mfChAlo1apVns+3atUKEyZM0Hi8JUuWoFatWvD398eLFy9yfc2GDRtQq1atPD93wYIFqFWrFj777LN8PyP7H19fXzRs2BCjRo3CtWvXNJ4vEZm2+Ph4fPLJJ2jRogUCAgLQpk0bTJkyBbdu3crx2jNnziAiIgJNmzaFv78/WrdujcmTJ+P69es5XtuqVascccTHxwdvvPEG+vXrh19++UXt9bVq1cKSJUtynWNiYiJq1aqFbdu2qb2+Vq1amD9/fq7vUSgUaNasWY73EZHpu3r1Kj788EM0bdoUderUQUhICD744ANcvnxZ9Zpt27ahVq1aSExMzPH+/fv3491330WzZs1U7x87dizOnz+v9rrff/8dtWrVwsiRI3OdR16f8fLlSyxfvhydO3dGvXr10KBBA/Tp0webNm1CZmZmrmO8/icgIACtWrXCjBkz8jyfIyLTpEmMOnPmDN59990Cx1Jed+ma8twp+586deqgSZMmeP/99/HHH38U+Ho/Pz+8+eabmDp1Kh4/fqzzOZJ2bI09ATINvr6+mDZtmupxVlYW7t27h2+++QaffvopSpUqhebNmxtsPg8ePEBcXByio6MN9plERK/LzMzE4cOH0blz5xzP7d27N8/3KRQK7NixAzVr1sTOnTvx8ccfw8nJKdfXbtq0SfV9VlYW7ty5gwULFqB///7Ys2cPypUrV/QfhIiMZv369Zg1axYaNmyIcePGwc3NDf/++y9iYmJw8OBBxMXFoXbt2gCAVatWYf78+QgJCcHEiRNRrlw5/Pvvv9iwYQO6deuG6OhovPXWW2rjN2/eXC35lJmZiYSEBKxatQojR47E1q1bVeMXho2NDfbv34+PPvoox3P/+9//8ODBg0KPTUTGce3aNfTu3Rv16tXD5MmTUbZsWdy7dw/r1q3D22+/jbVr16JevXq5vjczMxPjxo3Djz/+iM6dO2PKlCkoXbo07ty5g82bN6NPnz744osvEBoaqva+n376Cbt27cr1nOp1d+/exZAhQ/DkyRMMGDAAQUFBSEtLw4kTJ/D5559j9+7dWL58OVxcXNTet3TpUrXzpmfPnuHYsWP49ttv8fjxYyxcuFDr3xURGZ6mMWrLli25LiJ4Xa9evdCsWTO9zff9999HixYtAABpaWm4d+8evv32W/Tv3x9LlixBmzZt8n39zZs3sWTJEvz999/47rvv9DZPKhiT5AQAcHZ2zvVE6M0330Tjxo2xbds2gybJiYhMQf369bFv374cF3T379/H6dOn4ePjg//++y/H+44fP4579+5h/vz5eOedd7B792706tUr1894PfYGBQWhQoUK6N+/P7Zv367R6ggiMk1nzpzB559/jv79+2PSpEmq4w0bNkSbNm3QtWtXTJw4Edu2bcORI0fw5ZdfIiIiAqNHj1a9tkGDBujatSvGjRuHCRMmoGbNmvD29lY9X6ZMmRxxJDg4GHXr1kVoaCh27dpVpCR5/fr1cfr0aVy6dAm+vr5qz+3Zswc+Pj7466+/Cj0+ERne119/jdKlS2P16tWwtX2VEmjTpg06dOiA5cuXY9WqVbm+d+XKldi/fz8WL16M9u3bqz3XqVMnjBo1CtOnT0erVq3g6Oioeq5EiRL4/PPP0aRJE7i6uuY5N0mSMGbMGKSmpmLHjh2oUKGC6rkWLVqgY8eOGDhwIGbMmIF58+apvdfHxweenp5qx5o3b45Hjx5h3759ePnyJYoXL17wL4iIjKooMSo35cuXR/ny5fUxVQBA5cqVc5yLdezYEe+88w4mTZqERo0awdnZOc/XN2zYEHZ2dpg4cSKuXbumdp5HhsVyK5QvBwcH2NvbQyaTARCrI1etWoW2bduiTp06aN++Pb799lu192RlZWHVqlUICwtDQEAA6tWrhz59+uC3335Te93Zs2cRHh6O+vXro1GjRvjoo49w//59JCYmonXr1gCAyMhItVIGv/76K/r164egoCDViqy7d++qnt+2bRv8/f1x+vRp9OjRA/7+/mjfvj0OHz6MGzduYNCgQahbty7atm2LPXv2qM3nzp07+Oijj9CgQQPUrVsXgwYNwqVLl3T6+yQi48ut3EBeW/BCQ0Nx/PjxHFt09+/fj2rVquWZePr+++9Rs2ZNVazKvlpcE3Xq1AEA3L59W6v3EZFpiYmJgYuLS66rsMuUKYMJEyagdevWSE5OxtKlS1G9enWMGjUqx2vt7OwwY8YMyOVyrF69WqPPLlGiBACozuEK64033oCrqyv279+vdjwzMxMHDx7MsbKdiEzfw4cPIUkSFAqF2vFixYph4sSJ6NixY67vS0lJQUxMDDp06JAjQQ6InScffPABGjZsiEePHqk99+GHHyI5ORlRUVH5zu3o0aM4f/48PvnkE7UEuVJgYCAGDRqEXbt25VqyKjcuLi6QyWRFjodEZBiaxKgJEyZg+/btuH37tqrsm7Kcyddff40OHTqgbt26+P7773Nc6w0YMACTJk3CqlWr0KJFC/j7+6NPnz45ykX9/PPP6N69OwICAtC+fXvs3r0bbdu2zbN0XXb29vaIiIjA06dPsW/fvgJfX7JkSQBFP2+jomGSnACIO/aZmZmqP2lpabhx4wYiIyPx8uVLdOnSBQAQFRWFxYsXo3Pnzli5ciU6dOiAWbNmYdmyZaqxvvjiCyxfvhy9e/fGmjVr8Nlnn+Hp06cYO3YsUlJSAACXLl3CO++8g7S0NMydOxfTp0/HhQsXMHToULi5uWHp0qUAxDYU5fc7duxAeHg4KlSogPnz5yMyMhJ//vknevfurXYSptwC2KdPH6xYsQJOTk74+OOPMWLECLRo0QIrV66Em5sbxo8fj3v37gEAHj9+jD59+uDixYuYMmUKvvzySygUCvTv31+j7TtEZBqyx7Hsfwqrffv2yMrKwuHDh9WO7927N8/E0NOnT3H48GF07doVANCtWzfEx8fj4sWLGn/uzZs3AYhVBkRkniRJwvHjx9G4ceM8yy2FhoZi1KhRSE1NxYULF9CyZcs8L45KlSqFJk2a4KeffsrxOdnjXUpKCi5fvozx48fDzs4OYWFhRfo55HI52rdvnyNJfvLkSaSlpeXbD4KITFOLFi1w584d9OnTB+vXr8f169dVze86dOiAbt265fq+EydOIDk5Od+4UqtWLSxevBgeHh5qx728vBAREYEff/wRu3fvzvP9x44dg42NTb67mJXnYK/HQ4VCoYqFGRkZePToEbZu3Yrt27ejbdu2KFasWJ5jEpHp0CRGjRw5Es2bN0e5cuWwadMmVfkSQCyAGj58OObOnYumTZvm+hkHDhzATz/9hMmTJ2P+/Pl4+PAhIiIikJWVBQD47bffMHLkSFSoUAFLlixB//79MW3aNLVFmgVp3LgxbGxsctQmzx6rUlNTcfnyZSxfvhyNGjVCjRo1tPxtkS6x3AoBEDUl/fz81I7JZDLUrFkTixYtQsuWLXHz5k1s3rwZH330kWr7f0hICGQyGb766iv069cPpUuXxoMHD/Dhhx9iwIABqrEcHBwQERGBK1euoF69eli5ciVKlSqF2NhYODg4AADc3Nwwbtw4XL9+HT4+PgBEgsjX1xcKhQJffPEFQkJC8OWXX6rGrV+/PkJDQxETE4NPP/0UgAg4I0aMUJU2+O+///Dhhx9i0KBBGDJkCACxmqBHjx64cOECypcvj7i4ODx9+hQbNmxQndC9+eabCA0NxaJFi7B48WJ9/NqJSIdu376dI44VlaurK9544w21kiu3b9/GuXPnMHfuXKxYsSLHe3744QcoFArVzcV27dphxowZ2LhxY65NPLMn8ZUnSbNmzYKLi4tGdTuJyDQ9efIEaWlpObb+50a5a+T1pNLrqlSpgp9++gnPnj1TrTjasWMHduzYofY6W1tb1KlTBzExMapzqqIIDQ3F+vXr1Uqu7N27F61bt1adxxGR+ejXrx+SkpIQExODGTNmAABKly6NkJAQDBw4EAEBAbm+T7lyu2rVqmrHFQpFjhWfNjY2sLFRX5M3dOhQ/Pjjj/jss8/QqFGjXMuuJCYmolSpUmqlCV6nXETwerPPtm3b5nitq6sr+vXrhzFjxuQ5HhGZFk1iVOXKlVGmTBnY29urSpckJycDEKVOevToke9nZGZmIiYmRhVrXr58ifHjx+Ovv/5CnTp1sGTJEnh7e2Pp0qWqBQxly5bNdXdgXmxtbVG6dGkkJSWpHZ80aZJaGT5ALIZ4vUoDGR6T5AQA8PPzw/Tp0wGIppkLFy5ERkYGFi5ciOrVqwMQd9IkSUKrVq3UkjqtWrXCihUrcObMGbRp00aVxH78+DFu3LiBf//9F0eOHAEApKenAxA1Ops3b652YRUYGKharfn6Cc/NmzeRlJSEcePGqR2vXLkyAgMDcerUKbXjgYGBqu/Lli0LAKhbt67qWKlSpQBAVUv45MmT8PHxgbu7u+pns7GxwZtvvoldu3Zp9kskIqMqV65crklrQOxKKazQ0FDMnDkTL168gLOzM/bs2QM/Pz9UqVIl19d///33aNiwIezt7VUxplWrVti9ezfGjx+f46Ivt8S+8oSMTTuJzJdcLgcA1Yqk/ChXR9nZ2Wk0pvL1ANCyZUtViZaEhATMmzcP7u7uhY4hua1kDwoKgru7O/bv3w9fX1+kp6fj0KFDOeoBE5H5GDt2LAYPHoxjx47h5MmT+P333/HDDz9g9+7dmDhxIgYOHJjjPa8nwpUWLVqElStXqh0bPXo0IiIi1I7J5XJER0ejW7dumD59eq4lCyRJUqtBnJu8nl+xYgXKlSuHjIwMbNu2DTt27MCYMWPQu3fvfMcjItNTmBilpMkCgRo1aqhdl7m7uwMQZaXS09Px559/YtSoUWrnRR06dFAtztSUJEk5zq1Gjx6tWvmemZmJu3fvYu3atejTpw++/fZbnS/8Is0xSU4AgOLFi8Pf31/1uG7duujcuTPCw8Oxbds2lClTBk+fPgWAPEsM3L9/HwAQHx+P6dOnIz4+Hk5OTqhRowYqVqwI4NVF3dOnT1XJa00oPzu31Qaurq45aofntvIgr63OyvH//fffPINRSkpKvu8nIuOzt7dXi2OvP1dYbdu2xYwZM3D48GF07twZ+/btQ6dOnXJ97aVLl1QN7N54440cz+/atQv9+vVTO7Z161bV93Z2dihXrpxW8ZGITFPJkiVRvHhx3LlzJ8/XJCcnIyMjQ7WCvKA+BLdu3ULx4sVVN/sBceNfGfv8/f1Rq1Yt9OjRA8OHD8fmzZvV4l+xYsVUCxZepzye2/mOTCZDhw4dsH//fnz00UeqcghNmzZVnf8RkfkpWbIkwsLCVOVTLl26hE8++QTz5s3L9VxHeU13+/ZttcZy/fr1Q5s2bVSPe/bsmedn1qhRA6NHj8b8+fNz9IgCxI6aX3/9Nd/rL+WKduV8lGrWrKnavVO/fn1kZmZi6tSpcHZ2Zv8EIjOkbYxS0qS00uvxRbnzRaFQ4OnTp8jKyspxTSaXy9XOwQqSkpKCZ8+e5Wga6uHhoXbdGhgYiObNm6NFixZYsmRJjpuOZDisSU65cnV1xdSpU3H37l18/vnnAF41gIqLi8PWrVtz/OnQoQNevHiBYcOGoVixYtizZw/++OMPbN26NcdWFxcXFzx+/DjH5x49ehQPHjzIcVwZiB4+fJjjuaSkJJQuXbpIP6+LiwsaNGiQ68+1devWIiXYiMj0vL6yU7k1LzdlypRBo0aNsH//fvzzzz+4fPkyQkNDc33ttm3bUKxYMXzzzTdYu3at2p+qVavm2sDT399f9ad27dpMkBNZkJCQEPz+++9IS0vL9fnNmzejUaNGuHfvHurVq4cDBw7kuVLzxYsX+PXXXwusAV6jRg2MGTMGf/31l6qvi5Krq2uu51nAq8UOuS1IAMSumn///Rd//fUX9u7di3bt2hW48p2ITM/9+/cREhKCLVu25HjO19cXH374IdLT03Ntitm0aVM4ODjk6FHg7u6udj5TkGHDhqFOnTr47LPPcjT4VO5aPnToUJ7vV35+QfFw8uTJcHNzQ1RUVK7XkURkeooSo3SlbNmysLOzyxE3lAl0TZ06dQpZWVm5Lp56XfHixVG9enX8+++/2k6XdIhJcspThw4d0KxZM+zevRunTp1CcHAwAFFjM/tJ0OPHj7Fo0SI8ffoUN27cwNOnTzFw4EDUqFFDdTful19+AfBqi15wcDB+/fVXtdVMly5dwrvvvouLFy+qthMrVatWDeXKlcvR5OXWrVs4e/Ys6tevX6SftUGDBrh58yaqVaum9rPt3LkTW7duzTEfIjJfzs7OOVY+vt5M5XWhoaE4fvw4tm7diqCgoByrAQCxCvOHH35Aq1at0LhxYzRs2FDtT9euXXH58mWcPXtWlz8OEZmw8PBwPH36FAsXLszxXFJSEmJjY1GjRg34+flh9OjRuHnzJubPn5/jtVlZWZg2bRpSU1MxbNiwAj930KBBqFmzJmJjY/HPP/+ojjdo0ADHjh3D8+fPc7xn//79OXYWZlevXj14eHhg586dOHz4MFdlEpkpV1dX2Nra4rvvvsv1Bt6NGzfg4OCQa1k5FxcXDBkyBDt27MCPP/6Y6/hXr14tcA5yuRyzZ8/Gixcv8NVXX6k917RpUwQFBWHOnDm5JsHi4+OxZs0ahIaG5qiN/jpnZ2dERkbiv//+U+trRUSmS5sY9XrfA12Ry+WoX79+jubAhw8fVis9nJ/MzEwsX74crq6uufZLeN3z589x8+bNPEt6kmGw3Arla+LEiejcuTNmzpyJ7du3o3PnzpgyZQpu376NOnXq4ObNm1iwYAE8PT1RtWpVJCcnw9nZGStXroStrS1sbW1x4MABVTmBlJQUAMDIkSPRu3dvvPfeexg4cCBSU1OxcOFCBAQEoGnTpqrk+cmTJ+Hl5YW6devio48+QmRkJMaNG4fOnTvjyZMnWLp0KUqWLKlqyFlYgwcPxs6dOzF48GCEh4ejdOnS2Lt3LzZv3ozIyMii/RKJyKS0aNECe/bsQd26dVGlShVs27atwDv2bdu2xbRp0/DNN9/kaLKidOjQITx9+lS1HfB1Xbp0waJFi7Bx40ZVcxkismz16tXD2LFjsXDhQly/fh1du3ZF6dKlce3aNcTExCAtLU2VQG/WrBkmTJiAuXPn4q+//kKPHj3g5uaGxMREbNiwAX/99Rc+//xz1K5du8DPtbW1xcSJEzF48GDMmjULq1atAgCMGDECBw8eRL9+/TBkyBBUrlwZz58/x+HDh7F161ZMmzYt30acHTp0wNq1a1GqVCk0aNBAJ78jIjIsuVyOqKgojBo1Cj169ED//v3h5eWFlJQU/Prrr1i/fj3Gjh2rag78ujFjxuDevXuIiIhAhw4d0LZtW7i5uSEpKQlHjhzBvn374O7ujsaNG+c7D29vb4waNSrHTUQbGxt8+eWXePfdd9GzZ08MHDgQ9evXh0KhwIkTJ7B+/Xr4+vqq+mkVJDQ0FN999x22b9+Ovn375tmUlIhMgzYxqkSJEnj48CGOHj2qk0bl2Y0ZMwYDBgzAmDFj0LNnT9y5cweLFi0CkLN/S0JCgmohVEZGBhITE7Fx40ZcvHgRy5Yty1HaJfvrAVExYc2aNarKDGQ8TJJTvqpXr44BAwYgNjYWGzZsQHR0NL766its3LgR9+7dQ9myZREaGooPPvgAcrkcLi4uWL58OebOnYuxY8eiePHi8PHxwbp16zB8+HCcPn0arVq1gq+vL7799lt8+eWX+OCDD+Ds7IzmzZvj448/hr29Pezt7TFkyBBs2rQJR48exa+//oru3bujePHi+OqrrzBq1Cg4OzujWbNm+Oijj4rc3M7d3R0bN27El19+iaioKKSlpaFq1ar4/PPP862pR0TmJzIyEpmZmZgzZw5sbW0RGhqKcePGYfLkyXm+p0SJEggJCcGxY8fQvn37XF+zbds2lCxZEiEhIbk+X7FiRbzxxhvYt28fb74RWZH3338fvr6+WL9+PWbNmoVnz56hQoUKaNGiBUaMGIEKFSqoXjtkyBAEBgYiLi4Oc+bMwePHj1GuXDk0bdoUn3/+OWrUqKHx5zZu3Bjt27fHgQMHcOTIEbRs2RKVKlXC999/jxUrVmDx4sV4+PAhnJ2dUbt2baxcuRLNmzfPd8zQ0FDExMSgY8eOelu9RUT616JFC2zevBkxMTFYuXIlHj9+DHt7e/j6+mLBggVo165dnu+Vy+WYM2cOwsLCsGXLFsybNw8PHz5UXfdNmjQJXbt21aif0/Dhw/Hjjz/i4sWLascrVKiATZs2YcOGDdi9ezdiYmIgl8vh5eWFCRMmoFevXlrt9J08eTK6d++OGTNmYMuWLbk2KCYi06FpjOrevTuOHj2KUaNGYcyYMXmWxCyM4OBgLFmyBIsWLcLIkSPh4eGBKVOm4MMPP0Tx4sXVXrtixQqsWLECgFioUKZMGQQHB2Pq1Km59r3L/nobGxu4uLjAz88PMTExqgoOZBwySdlJkYiIiIiIiIiIiMiK/fTTTyhfvrxakvvatWsICwvD8uXL0bp1ayPOjvSFK8mJiIiIiIiIiIiIABw/fhx79+7Fxx9/jGrVquH+/ftYsWIFqlevnufOYTJ/XElOREREREREREREBCA1NRWLFi3CgQMH8ODBA5QqVQrNmjXDuHHj4OrqauzpkZ4wSU5EREREREREREREVosdd4iIiIiIiIiIiIjIajFJTkRERERERERERERWy2oadyoUCmRmZsLGxgYymczY0yEye5IkQaFQwNbWFjY2vN9WVIxRRLrFGKVbjFFEusUYpVuMUUS6xRilW4xRRLqlrxhlNUnyzMxMxMfHG3saRBbH398f9vb2xp6G2WOMItIPxijdYIwi0g/GKN1gjCLSD8Yo3WCMItIPXccoq0mSK+8s+Pv7Qy6X5/qarKwsxMfH5/uawtDHuOY0V32Ny7nqZ1yFQoE7d+6gYsWK+d6RU342VxboBmOU5Y3LuTJGWRLGKMsbl3NljLIkjFGWNy7naty5MkbpFmOU5Y3LuVrmeZTVJMmVW1rkcnmB//E0eU1h6GNcc5qrvsblXHU7rlwuR5UqVTR+PbeL6QZjlOWOy7kyRlkCxijLHZdzZYyyBIxRljsu56rbcWUyGdzc3DQuUcAYpRuMUZY7LudqWedRVpMkJyIiIiIiIiKyVjY2NqhUqZKxp0FEZJKYJCeiVxISgNhY4OpVoGZNYOhQgCdRRGQqGKOIyJQxRhEREREVjgmcRzFJTkRCXBwwfDiQkfHq2KxZwJo1wMCBxpsXERHAGEVEpo0xiohMnQkkoIiIcmUi51HswkBE4oTp9YAEiMfDhgG3bhlnXkREAGMUEZk2xigiMnVxcUCNGsD06cCGDeKrlxewdq2xZ0ZE1s6EzqOYJCcisaLg9YCklJEhnrdCaWlpmDhxIoKDgxESEoLYfH4PV65cQd++fREQEIBOnTrht99+M+BMiSwcYxQRmTLGKCIyZSaUgCIiy6FsBFzk5pkmdB7FJDkRiS13RXneQs2dOxcXLlxAXFwcpk2bhqVLl2L//v05Xvf8+XOEh4ejRo0a+OGHH9C2bVuMHj0ajx49MsKsiUyHzk6cGKOISA8Yo4jIKphQAoqILIeyEbCNjZapZUkCkpJePTah8ygmyYlI1KQryvMWKDk5GVu2bMGkSZPg5+eHtm3bYtiwYVi/fn2O127fvh3FihVDVFQUqlSpgjFjxqBKlSq4cOGCEWZOZDoKdeL08KHYErx48atjjFFEpAeFilGZmcDRo0D28wHGKCIyZSaUgDIlmu4aHjBgAGrVqpXjT2RkpIFnTGTG7t8XpZ6GDgWqVgU++eTVcyZ0HsXGnUTWKisLsLEBZDIgPFw0RchthYGdnXjeyly+fBmZmZkIDAxUHQsKCsLKlSuhUCjULqhPnTqF1q1bQy6Xq459//33Bp0vkVm7dg3YuRPYsQM4eRJQKIAqVYDRo0WcYowiImN6+RI4eBDYtg3Yswd48kTEqL59GaOIyPSZUALKlGTfNXznzh2MHz8eFStWRIcOHdRet2TJEmRki+/nzp3DBx98gH79+hl6ykSmo6BGwMnJwM8/A4cOAT/+CLy+gPD4cXHNZ2LnUUySE1mjY8eAUaNEIAoLAypXFl2Dhw1TD0x2duK4FXY9T0pKQunSpWFvb6865urqirS0NDx9+hRlypRRHb916xYCAgIwZcoUHD58GB4eHhg/fjyCgoK0/tysrKwCn8vvNYWhj3HNaa76Gtca52pjY/OqbEE+J06SJEFx7hxk330H2a5dkF25ojaOVLUqpMaNIUtMhKxyZY1jlCRJUCgUOvlZiMjCFXRx9+gRsH27SIwfOQKkpr56ztkZqFZN3OCrVYvnUURkmiSJC6LyoNw1vHr1avj5+cHPzw/Xrl3D+vXrcyTJS5Uqpfo+KysLCxYswLBhw+Dv72/gWROZiLi4nH0OZs0S5zwDB4rHhw8DnTqpv696daB+faBRI6BVK/F+BweTOo9ikpzImiQmAuPGAZs3i8dz5gDt24vgM3Ag0LKl+gVjeLjVXtilpKSoJcgBqB6np6erHU9OTsaqVaswcOBArF69Gnv27MHQoUOxb98+VKhQQavPjY+P18lrCkMf45rTXPU1rrXM1cvL69VFRAEnTjKZDPJbt4B58wAAko0NUry88KJePTxr1AgptWsjy8UFji9eoJZy50YBMUqhUODKlStITk4u9M9ARFZCk4u7EyfEa5Tc3ICGDYHmzYEWLYAKFYBsN8x5HkVEJuX0aeDZM6B1a5NKQJkKbXYNZ7dt2zY8e/YMw7P//0ELXBBlOeOa01wlSYKbmxskSSr0uKrFUAU1Am7ZUsSU0FBIQUGQypaFpDx/qlkTKFECKF5c3MADYCNJYtwCzqOyL4YC9LcgiklyImuQmiqSUdHRQEqKCEgtW4pVU3fvisAjk4mv06YZe7YmwcHBIUcyXPnY0dFR7bhcLoePjw/GjBkDAPD19cWvv/6KnTt3YsSIEVp9rr+/v1rZluyysrIQHx+f72sKQx/jmtNc9TWutc1VdTGh4YmT9NZbkHr0AGrWhBQaCoeqVeFQqhTKOjmpvU2SJEjKk6c8YpTyNTWzbRdW/kxERGo0vbh76y0gJESUVWnVSiTIXV2BUqXEqqfc8DxKTVpaGqZPn46DBw/C0dER4eHhCC9gxWpiYiI6deqElStXomHDhgaaKZGFWbMGiIgA6tQRJQ0cHHgj7zXa7BpWkiQJa9aswcCBA1G8ePFCfS4XRFneuKY612LFiqFChQooWbIkZDIZKmVLNj979gx3797VeHFR2bJlUbVqVfFAk0bA06YBNjZIXLMGj548QVbx4oBcDjx4IP5kU6pUKVSrVk1cS+ZxHqVQKHDz5k08ffpU0x+/0JgkJzJjMpkMbm5ur8obvE6SRJ3fDz4A/v1XHKtZU1wchoQANWqICz7Kwd3dHU+ePEFmZiZsbUWoTEpKgqOjI0qUKKH22nLlyqF69epqx6pWrYq7d+9q/blyubzAJKUmrykMfYxrTnPV17hWN1cNT5xkNjaQxcYCjo7Aa7s2tCWTyfTyOyYiC6TFxR22bgVsbYGSJcVX0oqm9X6zi4qK4o4gosJKTRXJ8TVrXh373/+Axo1Fgoo38lS02TWs9Pvvv+PevXt4++23C/25XBBlOeOa+lzVymBmI5PJUKpUKZQsWVJtZXYuEwHOnoXs6FHISpYE3n1XHNeiEbBHQAAqSlKBc1UueFLOL/tx5VdVkl41Pf0siOLZHpEZs7GxUd0RzOGvv0Td8SNHxOPSpYEBA0RdqJo1gYoVc73gUygUSExMhKenZ57bzKyBj48PbG1tcfbsWQQHBwMAzpw5A39//xy/l3r16uF///uf2rEbN24gLCzMYPMlKooCb7jlR5JEt/Ly5cVjLU6c8NoNJ00wRhFZnyLFKEA0jypWTHyvTYxyd9f6oxijBG3q/Srt2rULL1++NPBMiSzEjRtAv37A77+Lx/37Az17ivj3+DFQrpxx52ditNk1rHTgwAG8+eabajXKtcUFUZY3rsnPNY8eLDkWGWVkAH/8ARw9Kv4cPw789594btCgV0lyLRoBF+Y8KLfzKEMuhmKSnMjSSJJYPTByJJCZKRLhnToBffsCPj6i0VQ+28MkSUJSUhI8PDwMOGnT4+TkhK5duyIqKgqzZs3CgwcPEBsbi+joaABiVbmLiwscHR3Rp08frFu3DkuWLEHnzp2xY8cO3Lp1C126dDHyT0GkmXxvuOVGksTKpG3bgO+/B5o2Bb75RjynxYlTYTBGEVkfrWMUANy6JXbTff+9aBQVEyOOM0YZhLb1fp88eYJ58+YhNja2SIsMWO/XcsblXDUcV6EADhyAzbBhkN2/D8nREVJkJCQ/P7ELpk4d0T/h/8eQJAl37txBxYoV873xaOkN0LXZNax07NgxjB492pDTJCoaTXqwXLsGjB4N/Por8PqN6mLFAF9fsRjKQI2AjX0exSQ5kbnK444glE0PVq0S3w8aBAQFidIqZcuqGiRQwSIjIxEVFYVBgwbB2dkZERERaNeuHQAgJCQE0dHR6N69Ozw8PLBmzRp8/vnnWLVqFby8vLBq1Sq4F2IFGpHJyswEjh0TCaft24E7d149J0niIs3GRu8nTvSKNvV+r1y5gqioKFy8eBFVqlTBpEmT0KhRIwPPmEiPJAm4cEHEp23bgHPnXj138yZjlIFpW+939uzZ6NatG7y9vYv0uaz3a3njcq55jytLTYXbd9+h4urVkGVkIN3NDbdHjkSWkxMyXr5Eqrs7FAkJ4rrxNUlJSXqZq7nQZtcwADx+/Bi3bt1CUFCQoadKVDia9mDx8gIuXxYJcmdnkRQPDBS9WJo0Eb1YXFxevd/CGwEzSU5kjgq6I+jgAHz7LXD7tkiOe3iwlmYhODk5Yc6cOZgzZ06O565cuaL2OCgoCNu2bTPU1Ih0J68bboBIOh07Bnz9NbBrl9iuq+ToCNSvD7RoAXTsCKSni2MWfuJkSjSt9/v8+XOEh4ejVatWmD17Nnbu3InRo0fjwIEDKFu2rJFmT6ShgmLU2bPivGjHjlf9VwCxKKBmTVGLt1UrcfHn4sIYZSDa1Ps9ceIEzpw5g927dxf5c1nv13LG5VwLGPfWLcjGjYPN/19/SI0bQ/7pp6iclSWu/2rWzLXni6ZztfQG6NrsGgaAa9euwcHBAZ6ensacNpHmtOnBMnOmOKdq2lSUw3RxEdd1ebHgRsDMmhGZG03vCNaqJS4ElTU4iYheV9ANN5kMePbsVSkVZ2egQQORcGrfXtyAK10650mUBZ84mQpt6v1u374dxYoVQ1RUFORyOcaMGYOjR4/iwoULaN68uZF+AiINaBKjEhOBRYvEc3Z2orRASIiIQXXqiBhVqpT6YgHGKL3TtN5vamoqpk6dimnTpuVZB1gbrPdreeNa81yV/RhsbW1frW7OyABOnQJGjBC7ZwAgPByyHj0gVyhE3KteXSS+DDhXc6TprmEAePToEUqUKFH43hhEGipyHxYlbXqw9O2r/aJKC20EzCQ5kQHoLNABmt8RlMmYICeivGl6w+2tt4DBg8VFV9u2opFd6dK5rk5SY6EnTqZCm3q/p06dQuvWrdUuhr///nuDzpdIa9rEqJ49RaK7ZUuxbbhUKbESKr8EEGOUXmla7/f8+fO4desWxowZo/b+4cOHo2vXrpgxY4ZB501kSnL0Y8jIAPbuFTtqHj0S13pRUUDt2iLBVacOULGi0eZrbrTZNRwaGorQ0FBDTY2sWKH6sNy9K8rNPX8OjB8vjmnTg4VVB1T4myAygEIFOqX0dODvv0VtKEC7O4JERHnRZgveokWAk5NYpUkmQZt6v7du3UJAQACmTJmCw4cPw8PDA+PHjy9UXU02xbOccU1xrjKZ7NUNHi1ilGLZMki2tiIxnn1BQh5zUH5G9sULkiQBABQKhdbzzk7Tn9/Sm+JpWu83ICAABw8eVHtvu3btMHPmTDRt2tSgcyYyeXZ2QIcOwJtvAn/8AcyeLXbzlSgB+PuLBp1EZB1u3xY9WDZvFk03JQmoUgX45BOD9WBRKBRITEyEp6dnrrX8zRGT5ESm6NYtUVtz717gl1+AXr1elTvQ5o4gEVmcIu1MSUkRyW5Auxtu2Vb9acMST5xMhTb1fpOTk7Fq1SoMHDgQq1evxp49ezB06FDs27cPFSpU0Opz2RTP8sY1pbm6lSuHSpUriwdaxKjE1FStm9DZ29vD2dkZcrkcWVlZePHiRY5/O0VhybV8NaFNvd8qVarkeL+7uzt7JpB1y6sfg4MDsHEjcOaMWD1asSLg5wcUL27sGRNRUeXXhwUA7t8HNmwQifGTJ9XfW7OmuIF2/z5QoYJBerBIkoSkpCR4eHgUeSxTwSQ5kb4VFOgAIDMTOH5cJMYPHBDdhbP74w9AoTDYHUEiMl1a70x5/lw03dywAShfXpwUAQa54WaJJ06mQtN6v4CoO+rj46MqZ+Dr64tff/0VO3fuxIgRI7T6XDbFs5xx9TVXSZJw9+5dVKhQQfObeVeuQLZ2LWTp6cD8+eKYFjHK09OzUHHGmL9XS2+KB2hX75eIsimoH4O9PVC1qkiM165dcAk8IjJ9Bf27B0Q/gg8/fPV8rVpAixZAaCgQHCxKYioXRAHswVIITJIT6ZOmga5NG5HIUpLJRMALCgICAsT3jx8Drq4GuSNIRGbu5Utg927gu+/Ejbe0NHG8ShXecLMQmtb7BYBy5cqhevXqaseqVq2Ku3fvav25bIpneePqY0yNbuQ9eACsXw+sXQucPSuOVakCfPGF1jGqqDtVzOX3am60qfer6XNEFk/Tfgzly4s+MdypR2T+tOnD0rq1yA+FhgKBgTkT469jDxatMElOpC+aBjrlHT9AJMUDA4G6dcUWmVKlADc3Ueoge9KDdwSJrFN+O1MkCThyBFi2TJRqSk199b7y5YFmzYD27YHkZMDZmTfczJym9X4BoF69evjf//6nduzGjRsICwsz2HyJAIj4s2OHKCF3+PCrmuE2NmJRQJs2wMOH4tyHMYqIrJGm/RhkMvUeDERkvrTpFbVliyi7VKyYYedoJZgkJ9IXbQLdsmVidWfJkmJFQOnSIinOO4JEpFTQzhSZTKwg37ZNPOfmBoSEiMR448bicZky6s03ecPNbGlT77dPnz5Yt24dlixZgs6dO2PHjh24desWunTpYuSfgixOfjfyMjJEw7ljx169vnp1sU24QwfRoNzN7dXCAYAxioisg3KXH6Bdzxgisgza/LvPfp6kIfaJ0hyT5ET6ok2ga9pUnBi5uHDLHBHlpM0WvGHDgHr1gCZNxAryMmXEaoO88Iab2dK03q+HhwfWrFmDzz//HKtWrYKXlxdWrVoFd3d3I/8EZFEKupFnZwd88onou9KsGdCxo9g9V6GCiFOv1dJXYYwiIkv25Anw9ClQrZp4bICeMURUdDKZDG5ubpr3X1FKSgK+/16U0504URzT87979onSHJPkRPqiTaArxN1AgHcEiayGNjtT5s4Vyab8dqKQRdCm3m9QUBC2KXcZEOmapjfywsJEvwQ3N6BsWbE4gIjIGqWmAv/+C/zwA7BnD/DTT+wZQ2RGbGxsNOvBAogbYdu3i35RR46IcnNVqgATJvDfvYlhVo1Il54+FdvlABHIspc1yE5HgU55R1CSpCKPRUQmJvu/a2234GmZIFcoFEhISIBCGb+IiLShyY08QJSFCg4GqlbVKkHOGEVEFiMrC0hMBLZuBQYNEjtsfv5ZNFoHXvVjeP06kv0YiMzLixfAhg1Ap06ipG54OHDokIgB1auLBpz374vX8t+9yeBKciJdyMoCVq8W22XWrhUrpdhwisiqFWoLniQB584B69eLuDJ/vjjOLXhEZGoyMwHb/7+U0OZGXiF2vjFGEZFFePwYOHMGWLwY2LdPnOvZ2ABdugClSr2qTc5+DESmLb8eLFlZwOTJwKJFQErKq/d4er7qw9KokdhRV7Lkq+f5794kMElOVFSnTgHvvQecPSsef/WVCHy2tgx0RFZMqy14V6+K7XfffQdcuyaOVakCfPEFt+ARkWlJTASWLhWNgpcsEcdYQ5eIrFiBCyNSUsS53uLFwKZNIn4Con/MmDEiefbgAfDnn6JXg40N+zEQmaqCerDI5eLf9uzZYgX5m2+KPiwhISIxXqpU3gsG+O/e6JgkJyqspCSxPS4uTjwuVgzo3Rvo2hW4dw/w8BDbihnoiCg3CQnAxo0iMX7u3KvjdnbiAqlVK7HiyNWVO1OISOe02u0iScDRo8CCBaJ2rrKW5qJFvJFHRFYvz4URWVnA7dviXO2rr0QiHAC8vYGxY4GgIHFN+fixONerWrVQu22IyEA07cHy1lvAt9+Kf+OurqIcpi3Tr+aA/5WIXlPgRWNWFrBsGTBlCvDff+JYy5bAgAGAry9Qo4a4Q0hE1qugLXjvvvuqRi8gVhzUqSMS4+3avYoj3IJHRHqi0W6Xly9FGbklS4C//np13NcXCA0FHj0CypXjjTwiotdJkmjGGRHxquSUqyswciTQvj3w5ImoR1yxokiOu7qKBVbZKBQKJCYmwtPTEzZMnhMZnyY9WKZNEze7+vQpVGKc/+6Ni0lyotfke9F4/Djw/vvAhQvicbVq4oKwYUOR1PLwyDUQMtARWRFNtuB17Qp8/TVQu7ZIfLdrB/j5AWXKiMS4XJ772NyZQkSG8PffYpV4XBzw/Lk45uAANGsGdO4smm9WqiS2DCvxRh4RWaO8FkbIZEDz5kCJEqKh+sCBQL9+ovTKgweAm5u4lnR3z3P1OPsxEJmAtDRxDgRo14OlkCvH+e/euIyaJE9LS8P06dNx8OBBODo6Ijw8HOEFbMdMTExEp06dsHLlSjRs2NBAMyWrJ0nAwoXARx+Jx8WLA++8I2qPV68OeHmJY3m+nYGOyCposwXvhx/ExVSZMiLRlFdinIhI1/Lb7ZKeDvTtC5w+LR6XLy/Od8LCRImAihVF3MotqcMbeURkTQpaGOHgACxfLnbd2NoCDx+KnYJ+fiK2svwCkUFoVWIOEPmfX38V50r29sDKleI4e7BYPKNG5blz5+LChQuIi4vDnTt3MH78eFSsWBEdOnTI8z1RUVFITk424CzJquS3EmDkSFE7uEwZsXXG21v8cXPLsTWOiKyUNlvwOnYsVN1J7kwhoiIpKKljbw9MnSpi1VtviZt7VaqI8x0XF+PNm4jIlGi6MCI4GNi/X1wv1q8vbjTa2xtnzkRWSqMScwBw65YoM/f118D16+JYlSqAQsEeLFbCaEny5ORkbNmyBatXr4afnx/8/Pxw7do1rF+/Ps8k+a5du/BS2QmaSNc0WQmwYYMIll5eov4mT3CIKDtttuAVMsHNnSlEVGiaJnXCwkQJOXd3UXNci/Md3sgjIqug6cIImUw0ZC9ZUpRdISLTkpIC7Ngh/s3+9JNYRQ4Ajo5AkyZAp06iR4uLC3uwWAGjJckvX76MzMxMBAYGqo4FBQVh5cqVUCgUOU6qnzx5gnnz5iE2NhZhYWGGni5ZOk0vGr28xEqqEiWMM08iMi03b4qteO+8Ix5zCx4RmTJtkjr16hXqZh5v5BGRVdBmYUT58vqdCxHlLb8ScxkZohH5zz+/er2yOXnXriL/U6aM+mIB9mCxaEZLkiclJaF06dKwz/aXzdXVFWlpaXj69CnKlCmj9vrZs2ejW7du8Pb2LtLnZmVlFfhcfq8pymfqclxzmqu+xi3qmDKZ7NXNGC0uGhXOzpC0/ExL/L3q+mchMiat69T98Qcwdy6wZYs4IerXj1vwiMj0GWC3CxGRVeDCCCLTV1C1ADs70Xfu4kWgTRuRGG/cWPQOKFYs73HZg8ViGS1JnpKSopYgB6B6nJ6ernb8xIkTOHPmDHbv3l3kz42Pj9fJa/T12aYwprmNW9gx3dzcXtWl0uKiMTExEUlJSYX6TGv4vRKZI43q1EkScPCgSI4fPvzqeIkSwI0bQI0a3IJHRKbn1q1XsYdJHSIi3eDCCCLTpk2Jud9/FyvGS5RgvzkrZ7QkuYODQ45kuPKxo6Oj6lhqaiqmTp2KadOmqR0vLH9/f8jl8lyfy8rKQnx8fL6vKQx9jGtOc9XXuEUdU/bixasHWlw0enp6ar2F2BJ/r8rXEVm8jAxg0yaRHFf+nbexEasM+vcHmjdX30bLLXhEZAri40UDzj//FDfyuNuFiKjoEhMBT08ujCAyddqUmKtWTevh2YPFMhktSe7u7o4nT54gMzMTtrZiGklJSXB0dESJbPWez58/j1u3bmHMmDFq7x8+fDi6du2KGTNmaPW5crm8wGSiJq8pDH2Ma05z1de4Wo+ZlCSSXbt3i201Wl40FiUAWvTvlchc5VenLjMT6NYN2LNHPHZwAFq1EjXI33gDqFABcHbOOSa34BGRDmlVEkqZHN+x49WxkyeBpk2Z1CEiKooFC4CNG4FffhHnhFwYQWS6tCkxVwjswWKZjJYk9/Hxga2tLc6ePYvg4GAAwJkzZ+Dv76+WhAwICMDBgwfV3tuuXTvMnDkTTZs2NeicycwlJQFz5gArVgDJyeLYzz+LhBcvGomsU0F16mxtgffeA44dA956C+jdG/D3F8lxJyfjzZuIrIpGJaHOnxfJ8Z07Xx1r1Ajo21dsH5YksVqKSR0iIu1IEjBpEhAdLR4vWAB88gkgl3NhBJGpyMwEfvsNCAkRj1lijgrBaElyJycndO3aFVFRUZg1axYePHiA2NhYRP///3iSkpLg4uICR0dHVKlSJcf73d3dUbZsWUNPm0yIxquqlCvHly0DUlLEsWrVRLLLwQFQKMRqcl40ElkXTevUvfWW2HlSpQrg7i7ihga4BY+IDOL15LhMBjRsKJLjDRoAVasCbm7qNTaZ1CEi0oxCAbz7LhATIx4PHgzUrg2cOwfUqQO81meNiAwsMxNYvx747DPxPUvMUREYLUkOAJGRkYiKisKgQYPg7OyMiIgItGvXDgAQEhKC6OhodO/e3ZhTJBNW4KqqBw+AefNyT463bClWj1eoIAKoEi8aiayHpnXqlPXHbbX7Xya34BGRTuRXEiojQzw+fTrv5Dhv0hERFU56urh2VJav+uAD4M03RWz19lYlyLkwgsgIsifHr18Xx0qUEGXn6tZltQAqFKMmyZ2cnDBnzhzMmTMnx3NXrlzJ8335PUeEZ89EoFy+XD053qcP0LatuGh0dweKFTPqNInIyLSpU6dlgpyISCcKKgllZydWkc+apZ4cL1dOlAHIB5M6RET5ePlS7CY8elTE08mTgYAA0bSzTh21sntcGEGkGxpVC1Amx2fMEKvGAZEc79FD7PSoWvXVa1ktgLRUqKv+58+fY9euXbh58yZGjhyJc+fOwcvLC5UrV9b1/Ihyl9+qKkdHcTKTkgJUrw706wd07ChKJZQrV+CWOF40mj/GKNII69SRkTBGkUY0LQkVFgZ4eAAVK2qUHFdiUofywhhFVu/hQ7G46uxZcW05c6ZYdFW9OuDryxIrRsYYZbnyrRaQmQl8951IjitXjru4AN26AYMGiZtYZcuql5cDWC2AtKJ1BvDq1ato164dvv/+e2zcuBEvX77EwYMH0aVLF5w6dUofcyRSFxcH1KgBTJ8ObNggvnp5AWvXiucdHESDzqgoYMsWICJCNK7y8NDohEZ50ShJkn5/DtILxijK16VLorYkIFYR2Nnl/jrWqSM9YYwijWlSEgoQF4P16gHly2ucICfKC2MUWb2EBHHtePasWJ365ZfiWrN2bdYgNwGMUVZIkoCffxb/BgcNEglyFxdgwADg++9F/7kWLQBX15wJciItaZ0knzlzJvr27Ytt27bB7v+TC9HR0ejXrx/mzp2r8wkSqSloVdWtW+Jxy5aiZlxgIOtxWhnGKMrVkyfAqFGAvz+wd684pqxT93qinHXqSI8Yo0hj2pSE4nkO6QhjFFm1ixdFb4fr18XOnAULxPlgnToiQcfye0bHGGXBEhLEQsd+/cRXZW5HJhP9oUqXFsnxd94BNm8GZs8GWrUSpXQLOA9SKBRISEiAQrlYiigPWp9Rx8fHo2vXrjmO9+nTB3///bcu5kSUN21WVZUsyTuJVogxitQoFEBMjGiutHy5ePz990BWlnh+4EBxIZT9hOz6dXGcSA8Yo0hjLAlFRsAYRVbrxAmgaVPg3j2xkOLLL8UOnbp1xUpy3ow0CYxRFkqTagFffQVs2iT6sLRuLcrMaVlijtUCqCBaR/oyZcrg5s2bOY7/8ccfKFu2rE4mRZQnbVZVkVVijCKV06fFqoNhw4BHj0TJpenTgdGjXzX1BV7VqVu/XnzlCnLSI8Yoytf16ywJRUbFGEWWLtfGgHv3Am3aAM+eAT4+QHS0WJ0aGCgS5lx4ZTJ0GaPS0tIwceJEBAcHIyQkBLHKBXe5uHLlCvr27YuAgAB06tQJv/32m9ZzpzxoWi0gMFAkxytVyvv8iKiItE6SDx8+HJMnT8b69eshSRJ+++03LF68GDNmzMCQIUP0MUeiV7iqigqgyxilzYnT+++/j1q1aqn9OXLkSFF/HHqNRh3PHz0C3nsPaNAAOHUKcHISNevi4kST38BAwNnZcJMmyobnUZSrFy+ATz8VyRmWhCIjYowiS6dsDGijXBmemAh89plYQNGgATB1qlihGhQEVKhg3MlSDrqMUXPnzsWFCxcQFxeHadOmYenSpdi/f3+O1z1//hzh4eGoUaMGfvjhB7Rt2xajR4/Go0ePdPVjWTdtqgWwJwDpmdZFtfr06QM3NzfExMTA0dERc+fORbVq1fDZZ58hNDRUH3Mka5aUJGrBzZwptriFh4vtNbkFUa6qIug2RmU/cbpz5w7Gjx+PihUrokOHDjlee/36dcybNw+NGzdWHStZsmSRfx5Sl2/H86wskTiKjBQ1yAGgSRNRizwoCKhSBXB0zPWtCoUCiYmJ8PT0fHXRRKQHPI8iNZIkmox/+CFw5444tnEj0LGj2EI8cKDosxIbK3bL1awpznWYICc9YYwiq+PpKZoCRkcDvr4iQe7vL0p3ksnRVYxKTk7Gli1bsHr1avj5+cHPzw/Xrl3D+vXrc1zrbd++HcWKFUNUVBTkcjnGjBmDo0eP4sKFC2jevLmuf0Trw2oBZEIK1XmiVatWaNWqla7nQvRKejqwdCkwY4bY9takCRAW9mpV1bBh6olyrqqibHQRo7Q5cUpPT0diYiL8/f1Rrly5In0uFdKVK6Km+B9/iMeVKonV5B06ANWri0Yv+VDWqfPw8DDAZMna8TyKAAB//SVu4il3HZUrBwweDHTqBPz3H1CqlFg1pSwJRWQguopRaWlpmD59Og4ePAhHR0eEh4cjPI8FLbt27cKyZctw9+5d+Pr6YuLEiQgICCjyHIjUJCSo33QcOlTEWAcHYMoUEZerVQOKFzf2TCkfuohRly9fRmZmJgIDA1XHgoKCsHLlSigUCrVFM6dOnULr1q0hz1b/+vvvvy/U52Yp+yLl81x+rynKZ+py3KKOKZPJXv2OtagWoFAoClVX3BR/B4Yc1xLnquufRUnrJPnSpUvzfX706NGFngwRJElsM/7wQ+DaNXGsenUgNVU8J5NxVRXlS1cxSpsTpxs3bkAmk+W9wpl0I68LGwCoWlXsNnFyAt5+W5RXqV1bbJPlynAyIbo8j9ImAfX+++/j8OHDasdWrlyJli1bavx5VDCNSkI9fy76IyxaBGRmihv9nToBgwaJciuVK4uEDZER6DJGaboj7/Tp05g0aRJmzpyJ+vXr47vvvsPw4cNx+PBhFGeyknQlLi5n3eNZs8RCq4EDxe4dHx+NGwGScegqRiUlJaF06dKwz1a+w9XVFWlpaXj69CnKlCmjOn7r1i0EBARgypQpOHz4MDw8PDB+/HgEBQVpPf/4+HidvKYw9DFuYcd0evECPs2aQaZFtQBJknDx4kWkp6cbfL6GHlNf43KuBdM6Sf7777+rPc7KykJiYiL+++8/tG/fXmcTI8ui0UXjX38BH3wAHDwoHpcsKVaCjh4tGu5lfy9XVVEedBWjtDlxunHjBpydnfHpp5/i1KlTKF++PCIiIgq1/Y6rC16xsbF5FTMKurBxcIC0cCEUN28C9eu/Kq0iSaIMiwHma4gx9TWuJc5VX6sLikqX51EsCWV68i0JJUnA5s1iIcDdu+JYYCDw7rvAG28AXl5i9XguWBKKDEVXMUqbHXlJSUkYOXIkunTpAgAYNWoUYmNjcf36da4mJ90oqDFgy5bi+pIJcpOnqxiVkpKidp0HQPX49SRscnIyVq1ahYEDB2L16tXYs2cPhg4din379qGClnXr/f391Vakv/6zxMfH5/uawtDHuEUac/9+2AweDNnXX2tdLcDX19fw8zXgmPoa1xLnqnydrmmdJP/2229zPT5r1qz8E6Bk1fK9aHz8WKyoWrZMJLNsbYEePYDx44E6ddi5mLSiqxilzYnTjRs3kJqaipCQELz77rv48ccf8f7772PTpk3w9/fXav5cXSDY29ujTp064oGmFzZNmuCyjQ1SU1KAy5cNOl9Dj6mvcTlX/dNVjGJJKDNz65a4offzz+Kxm5sorRIaCnh7A+XL57vrhSWhyFB0FaO02ZHXsWNH1fepqan45ptvULZsWXh5eRXiJyDKhSaNAbkAyyzoKkY5ODjkuKZTPnZ8rYeRXC6Hj48PxowZA0Akan/99Vfs3LkTI0aM0Gb6kMvlBSYTNXlNYehjXK3GTE8XvaPmzxePv/pKqx4sMpmsyPM3+u/AyONyrgUrVE3y3AwYMADdu3dHZGSkroYkSydJwP79wDvviEQ5ADRuLJLjbdoUWAuOq6pIG9rGKG1OnEaOHIkBAwaoVmXWrl0bFy9exObNm7VOknN1gSCTyV6d6Gp4YSOTyeDTsGGh69SZ2u/AkONa4lz1tbpAX7SNUcYqCcXdLurUdrzkUxJKcnOD7MULSHZ2kDp1gvTOO0CtWhrvejHl34EhxrXEuZrqbpe8aBujtNmRp3Ty5EmEh4dDkiR88cUXhSq1whhlOePqtOaxFo0BC1Pz2BJ/r5Yeo9zd3fHkyRNkZmbC1lakxZKSkuDo6IgSJUqovbZcuXKoXr262rGqVavirnJXGBXs2jWgd2/gzz/F47ZtRdWAx48BV1f2YCGTobMk+dGjR+HA+omUm7wuGmUyoFUrUXPcxQUYN04kzAtosKfEVVWkDW1jlDYnTjY2NjnKFlSvXh1///231vO06tUFedHiwqaoN8xM9ndgoHE5V+PRNkYZqyQUd7u84uXlhVLK8igFlISSOTggc+5cJP74I1J8fJBWqhSyCrHrxdR+B4Yel3M1Hm1jlDY78pS8vb2xbds2HDlyBBMmTICnpyfq1aun1TwZoyxv3MKO6WZri0p164oHWjQG/Pfff/FYuYBLS9bwezVV2sYoHx8f2Nra4uzZswgODgYAnDlzBv7+/jmuJ+rVq4f//e9/asdu3LiBsLCwok/cGnz7LfD++8DLl4CzMzBiBNC9u+gfpWHuh8hQtE6St2rVKsc2lpcvX+LZs2cYP368ziZGFkKDOsL45hugRAlRd5wrwqmIdBWjtDlxmjBhAmQyGaKjo1XHLl++jJoFnZCTZrS4sCEydbqKUcYqCcXdLq+o/l+gYUkoeYsWqOTtDbi7a32+Y6q/A0ONa4lzNdXdLrqKUdrsyFNydXWFq6srfHx8cO7cOWzcuFHrJDljlOWMW+gxk5MhmzsXss2bRc8rLRsDVq5cGZUrVzbMXI0wLmOU4OTkhK5duyIqKgqzZs3CgwcPEBsbq7qeS0pKgouLCxwdHdGnTx+sW7cOS5YsQefOnbFjxw7cunVL1UeB8vD8uUiIf/edeFy7NvDxx0CjRkCNGrk2KWe1ADI2rZPko0ePVgtKMpkMdnZ2qFOnDqpUqaLTyZGZ07SOsK+velNOoiLQVYzS5sSpVatW+Oijj9CwYUMEBgbihx9+wJkzZzBjxgyd/3xWScMLGyJzoKsYZaySUNztkgstSkLJPT2L9FEm+zsw0Licq/7pKkZpsyPv/PnzkMvl8PPzUx3z8vLC9evXtZ4/Y5TljavxmJIEbNwIfPIJcPu2OHbsGNC8ucaNAYta89gif68mRpf5qMjISERFRWHQoEFwdnZGREQE2rVrBwAICQlBdHQ0unfvDg8PD6xZswaff/45Vq1aBS8vL6xatQru7u46/dnMiUwmg5ubW9514E+fBt5+G7h5U+R6uncX//7q1BGLI/N4H6sFkLFpnSTv3r27PuZBlkjTBilMkJMO6TJGaXri1K5dO0ybNg0rVqzAnTt34O3tjTVr1sCziMkQq/fsGVCypFYdz4lMna5ilLFKQlEutCgJRWTqdBWjtNmRt3XrVty+fRsxMTGqYxcvXoSvr69O5kJW4MwZYPRo4LffxONy5cR5Y1YWoFCI1eQaNAYk06fLaz0nJyfMmTMHc+bMyfHclStX1B4HBQVh27ZtOvtsc2djY5N7rxuFAvjyS2DiRCAzU9Qbj4gA2rUDfHzEtR2RCdMoST5gwACNOwWvXbu2SBMiC8KLRjIQfcUobU6cevXqhV69emk8NhXgp5+AGTOAgwfFVjxe2JAZ00eMYkkoE8KSUGTm9BGjtNmR17t3b7z99tuIi4tD8+bNsWvXLpw/fx5z584t9M9E5q/AlaoAcO8eMGECsHatWEnu4AD06gWEhgJubqL3VVbWqxJXbAxolpiPMhOSBHz2GRAVJR43aACMGQMEBoryKq+VCSQyRRolyRs2bKjveZAl4kUjGQhjlIXZvFk08c3IAD79FJg/H5DLeWFDZksfMYoloUwIS0KRmdPXeZSmO/L8/PywdOlSzJ8/H19++SW8vb0RExNj1aUMKJ+VqgCQlgYsWADMnCmaAQKirEr//kDVqkC1aqKkg5OTweZL+sNrPROTkKC+cGnoUHGdJpOJm1b794uyKr17i9XjFSuyegCZDY2S5KNHj9b3PMiSSJIIgrxoJANhjLIgS5YAY8eKONKgAdCxo2j6UqqUsWdGVGj6ilEsCWVE6enAlSuAvz9LQpHZ01eM0mZHXsuWLdGyZUu9zIMsiCQBp04BffoA//wjjtWoIa4rAwNFPK5cGXB2znMINgY0P7zWMyFxcTn7zs2aJc53Bg4UuznWrweSkkSC/LUSgESmTuua5CkpKdi0aRP+/vtvZGVlqY6np6fj0qVL2Ldvn04nSGYmMRG4cQN4801eNJJRMEaZKUkCJk8WJ1kA0LatWIlQv75agpwXNmTudBmjWBLKSO7dA3r0EIny48dZEoosCs+jyCTkt1K1Xj1R5/jZM6BfP6B9e3HdWbWqRvWO2RjQvDFGGVFCQs4EOSAeDxsmzoMqVRI7OSpVEnkfIjOjdZJ88uTJOHHiBJo0aYL9+/ejY8eO+PfffxEfH887fNbuzBkgLAzw9ORFIxkNY5QZyswE3ntPxAkA6NlT1K+rVw9wcVF7KS9syNwxRpm5338HunUD7t4V2/g3bRLloWxsWBKKLAJjFBmdJitVly8XC7Nq1BDJ8TJlWM7BSjBGGVFsbO5VAgBxPDZWnAfJZEyQk9nSOkn+yy+/YNGiRWjSpAmuXbuGwYMHo06dOpg9ezauXbumjzmSOdi5U9zJT04WAfHPP4GGDUWA5EUjGRBjlJlJSQHefhvYvVskmcLDxWohf3+geHFjz45I5xijzFhMDDBypFhBXrGi2O0SECAuDB0cjD07Ip1gjCKj0nSlanCwSI6XLfuqKSdZBcYoI7p6tWjPE5kBrf+PkpaWhqpVqwIAvL29ceHCBQBA7969cfr0aZ1OjsyAJImmet26iQR53brAunWiRALv5pMRMEaZkSdPgDZtRILczg748EPg3XdFHGGCnCwUY5QZSk8HRowQCZr0dCAoSJz7dOggbuhlS5ArFAokJCRAoVAYccJEhccYRUalyUpVQFxnlivHBLkVYowyopo1i/Y8kRnQ+v8qXl5eOHHiBAARlM6cOQMAeP78OdLS0nQ7OzJtGRnionHcOJEsb99ebI9r0gSwt1d7KS8ayVAYo8zE7dtASAhw4gRQrBgwaRLQv79YlenkZOzZEekNY5SZuXsXaN4c+OorkZR5+23g88+BZs0Ab29ALld7ubIklCRJRpowUdEwRlFhyGQyuLm5QVaYRVLZ4yVXqlIBGKOMKDw87zIqdnbieSIzp3W5ldGjR2Ps2LFQKBTo0qUL3nrrLYwYMQJXrlxBs2bN9DFHMkXPnom6wYcOiYvGIUOATz8VF4y53NFnHWEyFMYoM3D5smjMmZgIlC4tGna2aSM6oLN+HVk4xigzcvIk0L27aNRZrBgwejTQqRPg5ydiF5EFYoyiwrCxsUElbXpOSZIoz7lpE5CVBXzxhTjOlapUAMYoI6pcWfQGGDZMfceHnZ04zr5zZAE0SpLPmDEDnTp1QmBgIFq3bo19+/ZBoVCgQoUK+O6777Bz507Ur18fAwYM0Pd8yRT88w/QsaNIdDk4AJ98IpLk1aqxxAoZBWOU6clzRdHvvwOhocDjx6Km75QpwJtvArVq5ViRSWQpGKPM0KpVIimekSEako8fL3bK+fhwtwtZHMYoMghJAs6fF4nxzZuB69fF8SpVgLlzX/WmmTUr95IrXKlqtRijTEBWlrhWGzhQ9AaIjRU7O2rWFP8umSAnC6FRkvzmzZt45513UL58ebz11lvo1KkTvL29AQC1a9dG7dq19TpJMiG//QZ07gwkJYlVVNOni3rknp7GnhlZMcYo05PriqITJ8SK8ZQUoEYNYOpUoEGDPHegEFkKxijTk+eNvLQ0kRxfs0Y8btBA9EuoV0/ELVutN2ESmTzGKCqShAT1hNnQoa8SZpIE3LgBfPONSIxnL5diby96PLRqJfrUlC3LlaqUK8YoIztyBFi2DFi/XiySrFQJmDbN2LMi0guNzvS//vprPHnyBAcOHMD+/fsRExMDLy8vhIWFISwsDBUrVtT3PMkUbNoEDBokLiCrVAGio8VJjbu7sWdGVo4xykwEBYlGdzdvAhMniga/1atzBwpZPMYo05PrjbykJFFO5fffX9UfHzxYlFfx9GSsIovFGEWFFhcHDB+untCeNUsktAcOFHHz0iVg5kzxnJ0dEBwsriE7dhTngSVLipJWSlypSq9hjDKi69eBHj3EjaySJYE5c8QNLZ4TkYXSeDlM6dKl0adPH/Tp0wePHz/G/v37sX//fixatAh169ZFWFgYOnbsiNKs0WjWcl1ZJUmiSdWUKeJxYKA40WnYUARIIhPAGGVC8lpR5OAAbNwo6vxWrSpWC/EEi6wEY5QZKFFClIEqXhwYM0Y0JK9Th+c6ZBUYo0hrCQk5E+SAeDxsmEh0V6oEvPWW2Hns5yduRCoT4/n1oeFKVXoNY5QR/PcfEBYmEuTVqwNduojjMhkUCgUSExPh6ekJG+4IJgtSqL/NZcqUQb9+/bB27Vr88ssv6Nq1K7Zs2YI333xT1/MjA1OurFIFOoUCiIp6lSDv0AGYPx9o2pQXjWSyGKOMKC5OlCSYPh3YsEF89fIC1q4Vz9vbi/IFVaowQU5WizHKyBISxLlNv37i661b4riDg4hbsbEioRMczHMdskqMUaSR2Njca4cD4nhsrPjexgZYt06cEzZoALi6slE7FQljlAFkZQF9+4o+dKVLA5Mni8WSrq4AAEmSkJSUBEmSjDxRIt0qUmHFixcv4uDBg/jpp59w584dtG/fXlfzIlNhYyPKIly8KO74DxsG+PsDzs7GnhlRgRijDEzTFUVsekcEgDHKKAoqDeDgIEoAODoyiUNWjzGK8pW9vnhBz2cvp6IhrlSlgjBG6VFkJLB3rzgX+vRToHFj9qEjq6B1kvyPP/7AwYMH8eOPPyIpKQkhISEYOXIkWrVqBUdHR33MkQyloBIJV68C1aoxwUUmjTHKiDRZUcSts2TlGKOMSNMbec7O3OlCVosxijRWs2bRni+AcqWqh4dHkcYhy8IYZQDffgvMmye+f+89oF07sVOY50ZkBTRKkp88eRIHDx7EoUOH8PjxYwQHB2PEiBFo3749SpQooe85kiEUtLLK1haoXVusLCcyMYxRJkKbFUVEVoQxykRoeiOPF4FkZRijSCNXrwILFwJLl4prwvBwcb2YW1y1sxPPE+kAY5QB/f67WDgAiBrkvXsDPj4iH0RkBTT6mz5kyBD4+/tj6NChCA0NhZubm77nRYak6coqJsjJRDFGmQg9rygiMleMUSaCN/KIcsUYRfm6c0fUE4+JEXWKQ0NFM7/KlcWCqmHD1K8j7ezE8UqVjDdnsiiMUQZy+7ZIjKeni74s778vGu6ykgBZEY2S5AcPHkTlypX1PRcyFpZIIDPHGGUiuKKIKFeMUSaCN/KIcsUYRbl6+hSYO1esHk9JEceaNROlOCVJ7LoZOFAsqMpesjM8nAly0inGKANISQE6dQLu3xf/fj/5RPSiK13a2DMjMiiNkuQMSBaOK6vIzDFGmQiuKCLKFWOUieCNPKJcMUZZD5lMBjc3N8jyKyuVmipKqsyaBTx5Io4FBABjxgB9+gDFi6u/vlIlLqgivWKM0jNJEje8/vwTcHEBJk4EAgOBihWNPTMig2NhIeLKKiIqOoVClGTiiiIiMlW8kUdEVs7GxgaV8op1mZnA2rUi4Z2YKI5Vrw6MGgUMGcIVpUSWauZMYOtWQC4Hxo0DGjUS//aJrBCT5MSVVURUNKmpwIQJwJw5YgsuVxQRkSnJzBQ38uzteSOPiOh1kgT873/A4MHAX3+JY+7uwIgRwHvvARUqFDiEQqFAYmIiPD09YcM+VkTmY/t2YOpU8X14ONCunWjUKZcbd15ERsIkubWSJODKFaB2ba6sIqKimTABWLQIOHcO+PZbwMND1KkkIjIFn34K/PMPsGEDb+QRkfVKSFC/QTh0qIiHMhlQt64oo1KypFg1PnIk4O2t8dCSJCEpKQkeHh56/AGISKfi44EBA8T3HToAvXuLRp0ODsadF5ERaZQkHzBgQP51y7JZu3atxh+elpaG6dOn4+DBg3B0dER4eDjC81i1vGvXLixbtgx3796Fr68vJk6ciICAAI0/i7KRJHHB+PPPwPHjIghyZRWZMX3FKNLAjz+KBDkgtua9eAHIZFxRRJQNY5QRbdwILFggvl++XCSFXFx4I48oG8YoKxAXBwwfrr4gatYssSBq4EBxPbh6NZCVBdSvzxhJJoUxSg8ePgTeegt4+VI06Hz/faBOHaBECWPPjMioNEqSN2zYUC8fPnfuXFy4cAFxcXG4c+cOxo8fj4oVK6JDhw5qrzt9+jQmTZqEmTNnon79+vjuu+8wfPhwHD58GMVfbxxC+ZMk4IMPgMWLxeOvvwbefVfUEubKKjJT+opRVIBHj8SFFQC0aQN07qyqX8cVRUSvMEYZyYULIikOAJ06iRt5Dg68kUf0GsYoC5eQkDNBDojHw4aJhVKVKonV5EyOkwlijCqaHA1709OBrl2BW7eA8uWBTz4RCXJ3d6POk8gUaJQkHz16tM4/ODk5GVu2bMHq1avh5+cHPz8/XLt2DevXr8+RJE9KSsLIkSPRpUsXAMCoUaMQGxuL69evczW5NhQK0Xhl5UrxeORIoEoVkegqV864cyMqAn3EKCqAJInk0717ovP56NFArVqi5i8RqWGMMoJnz8QFYHKyuPAbOVLU2Pz/LcS8kUf0CmOUhYuNzb33FCCOx8aKhVJMkJOJYowqGrWGvZIkeg+kpQFOTqJsZmAgULWqUedIZCo0SpJHRkZqPGB0dLRGr7t8+TIyMzMRGBioOhYUFISVK1dCoVCorerp2LGj6vvU1FR88803KFu2LLy8vDSel1JWVlaBz+X3msLQx7haj6lQQPbee7D5+mtIMhmkUaMgdeggLhzLlBFb6yAuGu/cuYOKFStqvKVJL/M10pj6GtcS56rrn6Uo9BGjqACxscDOnaKpywcfiG16ZcoYe1ZEJokxysAUCqB/f+D6dcDVVZSY8/cHSpUy9syITBJjlIW7erVozxMZGWOUDil7EBw/DmzdKsrs1qwpKgsQkfEadyYlJaF06dKwz7bq0NXVFWlpaXj69CnK5JJsOXnyJMLDwyFJEr744otClVqJj4/XyWsKQx/jajRmVhaqTp+Osnv3QpLJcL9/fzz39kaKXI6MBw+ABw9yvCUpKUnncwWM+DswkXE5V7IIf/8NjBkjvn/7baBVK7ErhYjIFMycCezZIxqQf/wx8MYbYscLEZE1qlmzaM8TkfnKq2GvgwPQty+QmsqdwETZaJQk18fduJSUFLUEOQDV4/T09Fzf4+3tjW3btuHIkSOYMGECPD09Ua9ePa0+19/fH3K5PNfnsrKyEB8fn+9rCkMf42o8ZmYmZAMGwGbvXkg2NpA++gjlWrRAuYCAXC8YLfJ3YALjWuJcla8zBVwxYEAZGeKEKjkZ8PUFBg8WJ1w6/HtNZGkYowxo3z4gKkp8P3SouIlXowbLCBDlgzHKAkmS+CqTAeHhoklnbiVX7OzE80QmjDGqkApq2GtjAxQrZrz5EZkgrVeSS5KEn376CdeuXVMrt5Ceno5Lly5hzZo1Go3j4OCQIxmufOzo6Jjre1xdXeHq6gofHx+cO3cOGzdu1DpJLpfLC0wmavKawtDHuPmOmZEB9OsHbNsGyOWQjR8PWUgIEBAAFFCD02J+ByY2Lueqf7qKUZSHGTOA06eB4sWBDz8E/PwAFxdjz4rIbOgyRqWlpWH69Ok4ePAgHB0dER4ejvACkh2JiYno1KkTVq5caXmNsG7cEDfxJEk0E+7TR9zMszXaxkkis8PzKAvw6BEwZAjw7rtAWBhQubJIig0bpp4ss7MTx5W1ionMgLHOo95//30cPnxY7djKlSvRsmXLwv0g+qZpw14iUqP1VcNnn32GrVu3wtfXF+fPn0dgYCASEhLw8OFD9O3bV+Nx3N3d8eTJE2RmZsL2/y9ekpKS4OjoiBIlSqi99vz585DL5fDz81Md8/LywvXr17WdvvVITwd69gR++EFcHE6cCDRqpFGCnMic6SpGUS5OnBCrDwBxctWwIUsYEGlJlzFq7ty5uHDhAuLi4nDnzh2MHz8eFStWzNEAPbuoqCgkJycX9ccwPcnJolHns2di5XhEhEiQF6I0H5E143mUmfvlF7FI6vZtICkJ6NBBXAsOHCiSYtnLLoSHM0lGZsdY51HXr1/HvHnz0LhxY9WxkiVLFvnn0RtNG/YSkRqtk+R79+7FF198gXbt2qFDhw6IiopCtWrVMGHCBGTk9Y8wFz4+PrC1tcXZs2cRHBwMADhz5gz8/f3VmnYCwNatW3H79m3ExMSojl28eBG+vr7aTt86pKaKC8UDB8QKgSlTgKAgJsjJKugqRtFr/vtPrNBUKIBmzYAePQBvb5YwINKSrmJUcnIytmzZgtWrV8PPzw9+fn64du0a1q9fn2eSfNeuXXj58qWufhTTIUnixl18PFCiBDB+vGhMXq6csWdGZHZ4HmWmsrKAzz8Hpk8X52qVKgEjR74quwKIY0yKkZkzxnlUeno6EhMT4e/vj3Lmcm7Bhr1EhaJ1C9sXL16gTp06AICaNWvi/PnzsLW1xXvvvYejR49qPI6TkxO6du2KqKgonD9/HocOHUJsbCwGDhwIQKwqT01NBQD07t0bv/32G+Li4vDPP/9g8eLFOH/+PAYPHqzt9C1fSorYVnfggGjAMH26SJDXrcsEOVkFXcUoes2oUWLbXrlyYoVmrVpAHqWxiChvuopRly9fRmZmJgIDA1XHgoKCcO7cOSgUihyvf/LkCebNm4cZM2YU/YcwNUuWABs2iNqaH30ENGgAVK1q7FkRmSVdnkelpaVh4sSJCA4ORkhICGJjY/N87c8//4wuXbogMDAQnTp1wk8//VSkn8Oq3L4NtG4tEuAKBdCuHfD990D//mLBFJEFMcZ51I0bNyCTyVDJnHZesGEvUaFovZK8UqVKuHTpEipWrAhvb2+cP38ePXr0gCRJeP78uVZjRUZGIioqCoMGDYKzszMiIiLQrl07AEBISAiio6PRvXt3+Pn5YenSpZg/fz6+/PJLeHt7IyYmBu7u7tpO32LIZDK4ublBln0V58uXQGio2Gbn4CBqB9epwwQ5WRVdxij6f5s3A+vWiVXjY8eKXSlubsaeFZFZ0lWMSkpKQunSpdWaoLu6uiItLQ1Pnz5FmTJl1F4/e/ZsdOvWDd7e3kWaf/b6n3k9l99rivKZuY57/Dhsxo2DDICif39IrVuLciuSJFZWmtJcTWxcztW4c9X1z6IrujyP0rSUweXLlzF69Gh8+umnaN68OY4fP46xY8di69atqF27ti5/PMuzZw8waJCoQ+7kJPrFjB2b53maQqFAYmIiPD09c+zeJjIHxjiPunHjBpydnfHpp5/i1KlTKF++PCIiItC8eXOt56/v8ygbGxuRI9KwYa8kSbkuriiIOf1/WV/jcq6WeR6ldZI8PDwcH3/8MWbNmoXQ0FB0794dtra2+PPPPxEUFKTVWE5OTpgzZw7mzJmT47krV66oPW7ZsqXpNkUwAhsbG/U7mZmZwJgxIkHu6AjMnAn4+DBBTlZHlzGKTfEAJCaKxk+AKOPUujVQvbpRp0RkznQVo1JSUtQu7ACoHr/eGP3EiRM4c+YMdu/eXeT5x8fH6+Q1uvhs24cP4duvH+SZmXhety7uNmiA5MxMZP31V6HH1BVzGpdzNa+56puuYpQ2pQx2796NRo0aqXYUV6lSBYcPH8a+ffuYJM9LWhoQGQksWCAee3mJhFhYGFCsWJ5vkyQJSUlJ8OD1IZkpY5xH3bhxA6mpqQgJCcG7776LH3/8Ee+//z42bdoEf39/reav7/MoLy8vlCpVSuOGvc+ePStSrz9z+v+yvsblXM1rrgXROkneq1cvVK1aFcWLF4eXlxeWLl2KLVu2oE6dOoiIiNDHHEkTtrbA8uWicVWjRqJWMBPkZIV0GaOsvimeQiGaPykb4Q0dCtSuza27REWgqxjl4OCQ4yJO+dgxWymk1NRUTJ06FdOmTVM7Xlj+/v6Qy+W5PpeVlYX4+Ph8X1MYuY6bkQGbFi0ge/wYUqVKKDZhAryaNNH4vMegczXRcTlX485V+TpTo6sYlVcpg5UrV0KhUKitYu7WrVuutYSteQdgrruGlf7+G+jTBzhzRjzu0gWYOlXs9LPV+vKeyKwY+jwKAEaOHIkBAwaoGnXWrl0bFy9exObNm7VOkhviPEqSJBE7CmjYK0kSXFxcUK9ePa0/w5z+v6yvcTlXyzyPKtT/RUuUKIG0tDQAQLNmzXD16lU0bdrUfJoYmLuEBPVAN3SoCHQODqIkwm+/iVqcTJCTldJFjGJTPADz5gHHjonY8uGHonxTqVLGnhWR2dNFjHJ3d8eTJ0+QmZkJ2/9PiiQlJcHR0RElSpRQve78+fO4desWxowZo/b+4cOHo2vXrlrXKJfL5QWeBGvymsJQGzciAvj9d6BYMcjGj4e8Xj1xLqRlM2GDzNXEx+VczWuuhqCLGKVNKQMvLy+19167dg0nT55Enz59tJ67yZWE0pLy5kH2XcPS/zffVCgUkG3YANn770P24gUkFxcoPv4YGDIEqFhRORmDzVXf43KullnKQBcMeR4FiH+PygS5UvXq1fH3339rPXeDnUdJkjgnyqdhr0wmK/LnmNP/l/U1LudqXnMtiNZJ8r1792LChAn46KOPEBAQAAC4cOECFi9ejC+//BJt2rTR+SQpm7g4YPhw9S0zs2aJLTMDBwJyOVC/vqhJR2SFdBWjtFkBBbxqihcbG4uwsDDd/UDG8uefwJQp4vvBg4HGjVWrDoio8HQVo3x8fGBra4uzZ88iODgYAHDmzBn4+/urxaeAgAAcPHhQ7b3t2rXDzJkz0bRpUx39VPqVY0XlN98AK1aI78eMEY06vby0TpATUU66ilHalDLI7vHjx4iIiED9+vXRunVrredvSiWhtFGqVClUq1Yt15XjymM2CgVk27YBL14gpVo1JL7/Pl40aADFgwfAgwcGm6uhx+VczWuu+mbo8ygAmDBhAmQyGaKjo1XHLl++jJqm3PzyyBHRs+6tt0RjcyLSiNZJ8sWLF2P69Ono1q2b6tiCBQuwbds2LFiwgElyfUpIyJkgB8TjYcPEVppKlZggJ6umqxhl1U3xkpNh07s3ZBkZkN54A4qePUUCqoBGeEaZqwmPy7lyBVRudBWjnJyc0LVrV0RFRWHWrFl48OABYmNjVRdwSUlJcHFxgaOjI6pUqZLj/e7u7ihbtqxufig9U+vDIkmiIV1wMFCtGtChg+jBwjJQRDqhqxilTSkDpYcPH2LIkCGQJAmLFy8uVGNJkykJpSVVsz0gz13DMjs7SBs2QJo+HfY9eqC6ry/w2o0IQ8zVUONyrpZZyqCojHEe1apVK3z00Udo2LAhAgMD8cMPP+DMmTNa78YzqKVLge3bgeHDoZgzB4nPn7NhL5EGtE6S37t3T21lpVJQUBCioqJ0MSfKS2xs7t2JAXE8NjbPrTRE1kJXMcpamuLZ29vDxcUFbm5uuHPnDp4/fw73qCi4XbuGTBcX/NupE15mZiLj2rVCjW9OK2r0NS7nal5z1TddnkdFRkYiKioKgwYNgrOzMyIiItCuXTsAQEhICKKjo9G9e3ddTNt0yGRAaCjQpg1w7hxQpQrg7GzsWRFZDF3FKG1KGQDA/fv3VY07165dm2MxgqZMpiRUYRWwa1jm4ADZlCmAo2ORds+Y9O/AAGPqa1zOVf+McR7Vrl07TJs2DStWrMCdO3fg7e2NNWvWwNPTUxc/ku799x+wd6/4vn59SI6OSLpxgw17iTSgdZLc19cX69atw+TJk9WOb968md3H9e3q1aI9T2QFdBWjrKEpnnIlQfatvZIkQTZkCJCQAFm7dqjapQvg66v1hZg5rajR17icK1dA5UaX51FOTk6YM2cO5syZk+O5K1eu5Pm+/J4zOXn1YbG3F6vJWWKFSKd0FaO0KWWQnJyMYcOGwcbGBmvXrrXePlfcNUxUIGOdR/Xq1Qu9evXSfsLGsHMnkJYGlC8PNGum9Y4TImumdZJ8woQJGDp0KI4ePQofHx8AIoA8ffoUq1at0vkEKZuCal6Zck0sIgPRVYyyxqZ4wP8nzMPCgPbtIf/7b6B6dcC2UD2eAZjXihp9jcu5mtdc9Y3nUVooqA8LE+REOqerGKVNKYOvvvoKCQkJ+Pbbb1XPAWJRgouLi45/QhPGXcNEBeJ5lAbWrRNfmzcHrPWmI1EhaZ35CAgIwIEDB7B79278888/sLW1RcOGDdG5c2frOokxhvBwcXGY28mTnZ14nsjK6SpGWU1TvLxWadrZAbVrMwlFpGM8j9KQpisqiUindBmjNC1lcODAAaSmpuZYpdmtWzfMnj1bZz+byeOuYaIC8TyqAA8fAocPi+87dgRKlzbufIjMTKGWB5YpUwbdu3dHQkICvLy8kJGRAWfWg9S/ypXF6qlhw9QvGu3sxHFeLBIB0E2MsoqmeFylSWQUPI/SAFdUEhmNrmKUpqUM9u/fX6T5WgzuGibSCM+j8rF1K5CZKXq2NGwockUm2syeyBRp3do2LS0NkyZNQoMGDdCzZ088ePBAteXl2bNn+pgj/fOPqCkFiMTV9etAVBTQr5/4ev26OE5EOo1RkZGR8PPzw6BBgzB9+vQcK6D2KhuimKOCVmneumWceRFZOJ5HaYgrKomMgjHKCCRJfA0PFwmt3HDXMBEAxqgCrV8vvjZvDpj6gi0iE6R1knzevHm4fv06tm/fDgcHBwBAREQEnjx5gpkzZ+p8glYvKQlo0QLo2xdQNhGsVEmsnlq/XnzlCnIiFV3GKOUKqD///BPHjh3D4MGDVc9duXIF3bt3z/V9V65cQcOGDQv9MxiEJqs0iUjneB6lIa6oJDIKxigDO3PmVWkE5a7h1xPl3DVMpMIYlY/bt4FffxXfs9QKUaFonSQ/ePAgJk2ahFq1aqmO1apVC5999hl++eUXnU7O6qWmAp06Af/+C5w8CRw4ADx9+mq1ARHlwBilIa7SJDIKxigNcUUlkVEwRhlQXBzQtCkwYQJ3DRNpiDEqH5s2iVxRzZpAYCBgW6jqykRWTet/NS9fvoSTk1OO4wqFAlmsdaQ7kgQMGgT8/jtQvDgwcSLg6go4OgIyGRQKBRITE+Hp6anWRJDI2jFGaYirNImMgjFKQ+zDQmQUjFEGkJEBfPghsGyZeJyZCZw+DTRpIvrBKHcNE1EOjFH5WLdOfG3RgqVWiApJ6+xqq1atsGDBArx48UJ17NatW5g5cyaaN2+u08lZtSlTgM2bAbkcGD8eCAgA/P1FkhyAJElISkqCxFXlRGoYozTEVZpERsEYpQWuqCQyOMYoPbt3TySwlAny7t2BuXOBGjXYMJ1IA4xRefj7b+DPP0Uc6dABKFXK2DMiMktaJ8mnTp0KGxsbNGjQACkpKejRowfatWuHEiVKYPLkyfqYo/WJiwM+/1x8P3KkWFVQty7Ajs1EBWKM0hDrXhIZBWOUhpSrx9mHhcigGKP06LffgPr1gRMnACcnYNw44KOPgAYNAHd3tZcqFAokJCRAoVAYabJEpokxKg8bNoivfn5AnTostUJUSFr/y3FxccGSJUtw69YtXL9+HZmZmahWrRq8vLz0MT/r88svwPDh4vsePYC33hIryHknkEgjjFFaGDgQaNlSNOm8elWUWAkPZxKKSI8YozTw6BHQubOo0xsWxtWVRAbEGKUnX30FRESIG4AeHsCnn4p65L6+ImH+GuWuYQ8PDyNMlsh0MUbl4bvvxNdWrYAyZYw7FyIzplWS/MWLF5DL5XByckKlSpVQKVsiJSkpCfPmzcPcuXN1Pkmrce0a0KWLOHlq0kQksOrUAdzcjD0zIrPAGKWFjAyxapx1L4kMhjFKQ6tXi5WWERFAtWpQ+Pgg8fZt9mEh0jPGKD1ISxM7g2NjxeOGDYEPPhClNGvW5GpPIi0wRuUhPh64fFnEk7ZtucCSqAg0utK4d+8eBg8ejDfeeAP169fHe++9h2fPngEAsrKyEBMTg/bt2+Po0aN6naxFe/xY1I56+hTw9gbGjhUJcq4eICoQY5SWHj8Wqwx++EE0CSYivWKM0kJm5qtavZ06AeXKQQLYh4VIjxijikYmk8HNzQ2y13e9JCYCzZqJBLlMBvTvL/pONWoE+PgwQU6kIcaoAnz7rfhaty5Qu7boa0dEhaJRknzGjBm4ffs25s6diwULFiApKQnR0dG4f/8+evXqhS+//BJhYWHYv3+/vudrmdLTxQryGzeAcuWAyEixuqBaNWPPjMgsMEZpaeVK4PhxsUozPh6KrCzWvSTSI8YoLfzwg0gsOTsDXbsCrq7GnhGRxWOMKhobGxtUqlRJfafLH38AgYHA//4n4tmkScCQIUBwMFC1KstIEWmBMSofkgRs2iS+b92apVaIikij29dnzpzBwoUL0bhxYwCAr68vunXrhsuXL0OSJGzatAn+/v56najFkiRRA/j4caBYMZEgDwoSq8l58kSkEcYoLaSnA4sXi++7dgXc3VWrNFn3kkg/GKO0sGCB+Nq6NVCjhlgNlZVl3DkRWTjGKD3w8xOryP/8E/jkE7HCs04doGRJY8+MyOwwRuXj1CkgIQGwtxfnTiy1QlQkGiXJ//vvP7VGCJUrV0ZGRgY8PDywcOFC2NnZ6W2CFm/mTGD9esDGRnQ3f/NNsf2OW2SINMYYpYXvvgPu3wdKlwa6deMqTSIDYIzS0IULwLFjYpHA22+zJwuRgTBGFUFCgnoD9KFDRb8XBwdg40bg55+BsmXzbNBJRAVjjMrH2rXia3Aw4OUl8kpEVGgaJcklSYL8taStXC5HRESEdQekotqwAZg6VXw/dCgQGipWHfB3SqQVxigNSRLwxRfi+9BQcSLFVZpEescYpaFFi8TX4GCgfn0mlIgMhDGqkOLigOHDRTN0pVmzgDVrgIEDxcrOxo1FLGP9caJCY4zKQ1YWsHWr+J6lVoh0oki3mYoXL66reVifkyeBwYPF92+9BbzzDuDvDzg6GnVaRJaEMeo1P/0EXLwoVjf16sVVmkRGxhiVzZMnYmcdIHa5uLsbdz5ExBiVn4SEnAlyQDweNgy4dUs8dnZmgpxIT6w+Rh09Cjx4IMr2tmzJck5EOqDx/7H37dsHZ2dn1WOFQoGDBw+ibNmyaq/r2rWrziZnKXJ0PL9xA+jUSdQGDg4GxowRjTqz/X6JSDuMURqYM0d8bdlS1MW0tzfufIisCGNUAdasAVJSRJmCtm1ZU5PIwBijtBQbmzNBrpSRIZ6fNo09poh0hDEqF8pSKw0biobALLVCVGQaJckrVqyI2NhYtWNly5bFeuWKn/8nk8msKyhpSNnxHIAod/Dvv0C1akCJEsDkyUC9erwYJCoCxigNXLwIHDokLtZ69wbKlzf2jIisBmNUAbKygKVLxfdhYYCnJxNLRAbEGFUIV68W7Xki0hhjVC7S04GdO8X3rVuLflNEVGQaJckPHz6s73lYD5lMrOI8fhw4cUI0eGHJA6IiYYzSwLx54mtwMPDGG4C1b08kMiDGqALs2SNKFxQvDnTpwobCRAbGGFUINWsW7Xki0hhjVC4OHACePhUlVt58k6VWiHSE+zH0LSEBiIoC+vUTX5X16RwcgBYtAA8PI06OiKzC/fuiUTAAvP02ULGicedDRJTdwoXia8uWQK1arN9LRKYvPBzIq2GgnZ14nohIX779Vnxt3FiUquMOPCKdYJJcn+LigBo1gOnTRYJq+nTAy+tV7SgGMiIyhCVLxJY8Ly9xc47lnYjIVPz1F3DkiDgnYkNhIjIXlSuLXgqvJ8rt7MRxZalNIiJdS04Wu/AAoE0blloh0iEu1dGXgjqet2zJkyci0r+UFGDFCvF99+5caUBEpkW5irx+fSAoCChWzKjTISIqkEIhGuQNHCiu6WJjRQ3ymjXFCnJe4xGRPu3aJRLlrq5Akyai1x0R6QRXkuuLJh3PiYj07ZtvgMePxUnUW2+x1i8RmY5nz4B168T33boBFSoYdz5ERAWRJGDyZCAtTTyuVAmYNg1Yv158ZYKciPRNWWqlaVNRRpMLoIh0hklyfWHHcyIyNoUCmD9ffB8WBlSvDsjlxp0TEZFSTIxYCeXhwe3CRGQe1q0DoqOBVq3EzmFJMvaMiMiaPHsGHDokvm/bludORDrGJLm+sOM5ERnbnj3A33+L8gU9erDWLxGZDoUCWLpUfP/WWywFRUSm79kz4OOPxfeVKwOPHgGSBIVCgYSEBCgUCuPOj4gs39atoteUh4coVefiYuwZEVkUJsn1hR3PicjY5s4VX1u3Bnx8AAcH486HiHQqLS0NEydORHBwMEJCQhCbTym3Xbt2oX379ggICECfPn1w/vx5A840F/v2ATdvipt4XbsC5coZdz5ERAWZMgV48ECUhurXT+zQs7GBJElISkqCxFXlRKRvyjJ1ISEiUc4FBkQ6xSS5vrDjOREZ059/AsePi/Iqb78NuLsbe0ZEpGNz587FhQsXEBcXh2nTpmHp0qXYv39/jtedPn0akyZNwsiRI7Fnzx4EBgZi+PDhePnypRFm/f8WLBBfW7QQu+vyWlhARGQK4uOB5cvF90OHArVrAyVLGndORGRdkpKAY8fE9yy1QqQXtsaegEXKyBAXe+x4TkTGMnu2+NqwIRAYCDg7G3c+RKRTycnJ2LJlC1avXg0/Pz/4+fnh2rVrWL9+PTp06KD22qSkJIwcORJdunQBAIwaNQqxsbG4fv06AgICDD/5q1eBn34Sq5969uRNPCIybZIEvPcekJUlzqvatRMLooiIDOm770Qcql4dCAjg9R2RHjBJrmu7d4uk+IYNorSBsuM5EZGhJCYC27aJ799+W3Q9JyKLcvnyZWRmZiIwMFB1LCgoCCtXroRCoYCNzavNgh07dlR9n5qaim+++QZly5aFl5eX1p+blZVV4HP5vQYAZPPnwwaAVLcuFPXrA05O4qKviONqQx9jmtu4nKtx56rrn4X0aN064ORJcW03bBhQowZL2BGR4X33nfjarBlQvjxLrRDpAZPkupSUJFaKJyUBkZHAzJmi1iYRkSEtWABkZoqtwCEhQKlSxp4REelYUlISSpcuDXt7e9UxV1dXpKWl4enTpyhTpkyO95w8eRLh4eGQJAlffPEFihcvrvXnxsfHF+k1Ni9fImDtWgDAnYYNkfTwIbLOntXZZ2tLH2Oa27icq3nNlQzsv/+AcePE9z16AEFB3P1CZCHS0tIwffp0HDx4EI6OjggPD0d4Ab3jEhMT0alTJ6xcuRINGzY00EwBJCQAp06J79u0YakVIj1hklxXJAkYPlwkyD09xerN/19hoFAokJiYCE9PT7WVXUREOvf8ObB6tfi+e3cRj7jKgMjipKSkqCXIAagep6en5/oeb29vbNu2DUeOHMGECRPg6emJevXqafW5/v7+kMvluT6XlZWF+Pj4fF8jW7wYNikpkCpUQPl33kH5hg2BAs6NNBlXW/oY09zG5VyNO1fl68jETZokru8qVAD69FE16yQi85e9t8udO3cwfvx4VKxYMUfZuuyioqKQnJxswFn+P2XDTh8f8YelVoj0wqhJcm3u3P38889YsGABEhIS4OnpiQ8++ACtW7c28Izz8c03wM6doknelClAnTrie0DV8dzDw8O4cyQiy7d6tUiUly8PdOgAuLoae0ZEpAcODg45kuHKx46Ojrm+x9XVFa6urvDx8cG5c+ewceNGrZPkcrm8wGRinq9RKIBlywAAstBQyKtV06phpyafrS19jGlu43Ku5jVXMqDz54EVK8T3bNZJZFG06e2itGvXLuM1Pd+4UXxt1oy7WYj0yKhJck3v3F2+fBmjR4/Gp59+iubNm+P48eMYO3Ystm7ditq1axtp9tncvAlERIjv+/cHwsJ4Z4+IDC8rC1i0SHzfqRNQtarqZh0RWRZ3d3c8efIEmZmZsLUVp3NJSUlwdHREiRIl1F57/vx5yOVy+Pn5qY55eXnh+vXrBp0zDhwArl8HHB2Brl2BcuUM+/lEZFCFKWVw+vRpjB8/Hj/99JOBZpkHSQJGjHjVrLNtWzbrJLIg2vR2AYAnT55g3rx5iI2NRVhYWKE/t1C9Xa5cgTw+HpKNDRQtWwIlSuTby0XjcYvAnHqF6GtcztUye7sYLUmuzZ273bt3o1GjRhg4cCAAoEqVKjh8+DD27dtn/CR5VhYwYADw8iVQq5aoWVehgnHnRETW6fvvRb06Z2eRgHJzM/aMiEhPfHx8YGtri7NnzyI4OBgAcObMGfj7++e4sNu6dStu376NmJgY1bGLFy/C19fXoHPGggXia4sWYkXma+ViiMiyaFvK4MqVKxg7diwcTKEp5tq1olmnvb1o1untzWadRBZE294us2fPRrdu3eDt7V2kzy1Mb5eKS5eiAoCUGjXwT2YmUq5d09tnm8KY5jYu52pecy2I0ZLk2ty569atGzIyMnKM8fz5c60/t1B37vIhmzsXNr/+CsnREYpp00S3c4WiyOMWxJzuBOlrXM7VMu/cURHMmye+tmsnbtrxYo7IYjk5OaFr166IiorCrFmz8ODBA8TGxiI6OhqAuPhzcXGBo6MjevfujbfffhtxcXFo3rw5du3ahfPnz2Pu3LmGm/D168ChQ+L7Hj24VZjIwmlbymDjxo2YM2cOKlWqhBcvXhhhxtn89x/wySfi+5492ayTyAJp09vlxIkTOHPmDHbv3l3kz9W6t4skweaXXwAAju3aodabbwJalvG19l4h+hqXc7XM3i5GS5Jrc+fOy8tL7b3Xrl3DyZMn0adPH60/tzB37vLidOUKak+dCgB40K0b7rq5Ieuvv4o8rjbM6U6QvsblXM1rrqQnJ04Ap08Dtrbigo4Xc0QWLzIyElFRURg0aBCcnZ0RERGBdu3aAQBCQkIQHR2N7t27w8/PD0uXLsX8+fPx5ZdfwtvbGzExMXA3ZJxYsECUL6hbFwgOBlxcDPfZRGRw2pYy+OWXXzBnzhy8ePECS5cuLfTn6mJBlGzCBNgkJUGqUAGKt98GqlQR8SuP91n7Aht9jcu5ckGUPmna2yU1NRVTp07FtGnT8uz5og2te7v88Qdw4wZgZwebVq1EqbpCJiPNqf+GOY3LuZrXXAtitCS5Nnfusnv8+DEiIiJQv379QjXu1PrOXV5SU2EzYABkmZmQ3ngDrlOnwjWPrTfWfidIX+NyrpZ5544KafZs8TUkBAgIYF8EIivg5OSEOXPmYM6cOTmeu3Llitrjli1bomXLloaamroXL0TpAkD0S2AjcyKLp20pg+XLlwMAtm3bVqTPLeqCKMdr1+D71VcAgLuhoXisUCDt5k2dfba2zGmBjb7G5VzNa67mQtPeLufPn8etW7cwZswYtfcPHz4cXbt2xYwZM/Q70W++EV/r1weqVQOcnPT7eURWzmhJck3v3GX38OFDDBkyBJIkYfHixTlWIGhC6zt3eZkyBfjrL6BkSchmzIC8Ro0C7+iZ090VcxqXczWvuZIe3LgB7Nkjvu/ZE6hY0bjzISLK7uuvgefPxQ6X9u2BsmWNPSMi0rPCLogqqiItiJIk2IwcCZlCAalRI7j37w/3hg0LLF9n7Qts9DUu58oFUfqkaW+XgIAAHDx4UO297dq1w8yZM9G0aVP9TlKSRM8pAGjWjP2miAzAaElyTe/cKd2/f1/VuHPt2rU5Vh8Y1OHDr5pPffwx0KSJKHFARBYlLS0N06dPx8GDB+Ho6Ijw8HCEh4fn+tpdu3Zh2bJluHv3Lnx9fTFx4kQEBAQYZqLz5oleCHXrAo0aAaVKGeZziYgKIkmAsnRCaKgoW1CIRQ5EZF4KsyBKF4q0IOqbb4BTpwB7e8iGDoW8Vi2gWDGdfra2zGmBjb7G5VzNa67mQpveLlWqVMnxfnd3d5TV903/Y8eAO3cAR0egeXPAmDkwIithtKuU7HfulHK7cweIxi/Dhg2DjY0N1q1bZ9gamq97+hQYNEh8366d+D6XpD4Rmb+5c+fiwoULiIuLw7Rp07B06VLs378/x+tOnz6NSZMmYeTIkdizZw8CAwMxfPhwvHz5Uv+TfPLkVRmDrl0BT09AJtP/5xIRaeLQIeDqVbESs0sXUUuTiCxe9gVRSvktiDK67M06e/USzTrLlzfunIhIryIjI+Hn54dBgwZh+vTpOXq77N2717gTVF7jNWwIVKokkuVEpFdGW/6szZ27r776CgkJCfj2229VzwFiFYKLoRs/RUQAiYliy/DUqayrSWShkpOTsWXLFqxevRp+fn7w8/PDtWvXsH79enTo0EHttUlJSRg5ciS6dOkCABg1ahRiY2Nx/fp1/a8mX7YMSE4GKlcGWrcGXF31+3lERNqYP198bd4cqF2bF3hEVkLTUgYmY/x44OFDUbKud2+genXueiGycNr0dtH0OZ3JygJ27BDfh4Sw1AqRgRi1RkhkZCSioqIwaNAgODs757hzFx0dje7du+PAgQNITU1Fr1691N7frVs3zFY2qzOEzZuBdevECdPUqaJ5Ak+eiCzS5cuXkZmZicDAQNWxoKAgrFy5EgqFQu0Cr2PHjqrvU1NT8c0336Bs2bLw8vLS+nPz6ySfoxt9ejpsli6FDIAiLAxSpUrKFxbqM3XZxV4fY5rbuJyrceeq65+FCuHmTeDAAfF99+5clUlkRbRZEGV0f/4JrF4tvg8PFzf0SpY07pyIyLodOAA8egS4uABNm7LUCpGBGDVJrumdu9zKGxjc7dvAiBHi+169gB492FmYyIIlJSWhdOnSak2nXF1dkZaWhqdPn+baF+HkyZMIDw+HJEn44osvULx4ca0/V5MGOcrXlNm9G9Xu30emiwtu1K6NF3fvQnr0SOvP1OazTWFMcxuXczWvuZIOLVggapL7+wMNGjDpRGRlNF0QZVSSBLz/vlhg0Lgx0KaN2J1HRGRM/19FAY0bi3KaBTQQJiLdYLdJTUgSMGSIqP1bvTowcaIot0JEFislJUUtQQ5A9fj1RlRK3t7e2LZtG44cOYIJEybA09MT9erV0+pz8+s0r9aN3sYGNv/fH8EmNBRe7doBNWpo9Vm5jqujBj76GNPcxuVcjTtX5evISJKTX9XS7NRJlDAgIqtSmFIG3bt3N2ziPDYW+P13kYAaMgTw9mYyioiMKz0d2LNHfN+sGfu5EBkQk+SaWL4c+PFHwN4emDFDbMEjIovm4OCQIxmufJzX1mBXV1e4urrCx8cH586dw8aNG7VOkmvSaV4ul0N+5Ahw8SJgbw+bHj1Ef4QiJiH10eVeH2Oa27icq3nNlXTkm2+AZ89EDc127dgvgYhMz9OnwIQJ4vuePYHgYJaFIiLj27EDeP5clFhp0AAoXdrYMyKyGiyoXZDLl191Oh86FAgNFclyIrJo7u7uePLkCTIzM1XHkpKS4OjoiBIlSqi99vz587h48aLaMS8vLzx58kSnc5LJZHBzc4NMJgPi4sTBli0BX1/A2Vmnn0VEpC1VjAKA778XBzt2BKpWLfJNPCIinWOzTiIyIarzqPv3gSpVRC1yDw/ubiEyIK4kz09GBjBgAJCSAtStC4wbx7t4RFbCx8cHtra2OHv2LIKDgwEAZ86cgb+/v1rTTgDYunUrbt++jZiYGNWxixcvwtfXV6dzsrGxQSVlc864OHFBl5EBVKig088hIioMtRj144/A/v1AsWJiNTkRkZGpLTaIjwfWrBFPsFknEZkA1XnUBx8AY8aIBuilShl7WkRWhUnybNROnABg5kzg9GmxQnPGDLESioisgpOTE7p27YqoqCjMmjULDx48QGxsLKKjowGIVeUuLi5wdHRE79698fbbbyMuLg7NmzfHrl27cP78ecydO1d/E7SxAcLCAIUCUMYsIiJTYWMjdt8pFFyZSUQmQe1Gnp8fsHMnsH07m3USkemxsQG8vER/PCIyGCbJs1E7cZIkUf8pOBho3x5o1YpbhYmsTGRkJKKiojBo0CA4OzsjIiIC7dq1AwCEhIQgOjoa3bt3h5+fH5YuXYr58+fjyy+/hLe3N2JiYuCuywa/CQmiudTVq0DNmqL8U6VKTD4RkWlgjCIic6JcbNCxo+ifwHIGRGRMeZ1HcTEUkUExSZ4XmQx46y2gbVvRNIH1fomsjpOTE+bMmYM5c+bkeO7KlStqj1u2bImWLVvqZyJxccDw4aK0itKsWWKb8MCB+vlMIiJNMUYRkanLKwEll7OcJhEZF8+jiEwGl/coJSQAUVFAv37i661b4ri9vegqTERkDAkJOU+aAPF42LBXsYqIyBgYo4jI1MXFATVqANOnAxs2iK9eXsDateJ5rtQkImPheRSRSWGSHOCJExGZrtjYnCdNShkZ4nkiImNhjCIiU8YEFBGZMp5HEZkUJsl54kREpuzq1aI9T0SkT4xRRGTKmIAiIlPG8ygik8IkOU+ciMiU1axZtOeJiPSJMYqITBkTUERkyngeRWRSmCTniRMRmbLwcMDOLvfn7OzE80RExsIYRUSmjAkoIjJlPI8iMilMkvPEiYhMWeXKorP56ydPdnbieKVKxpkXERHAGEVEpo0JKCIyZTyPIjIptsaegNGFhwOzZuVecoUnTkRkCgYOBFq2FOWfrl4VN+/Cw3nSRESmgTGKiEyVMgE1bJj69R4TUERkKngeRWQymCTniRMRmYNKlYBp04w9CyKi3DFGEZGpYgKKiEwdz6OITAKT5ABPnIiIiIiIiCwVE1BERERUACbJlXjiREQmTqFQIDExEZ6enrCxYUsJIjItjFFEREREhcPzKCLj4788IiIzIUkSkpKSIEmSsadCRJQDYxQRmTKFQoGEhAQoFApjT4WIKAeeRxEZn9WsJFcGmqysrDxfo1AocPv2bXh4eOj0zp3yM/P7bFMY09zG5VyNO1fl8/yfuG5oEqMs8e+RJY/LuTJGWRLGKMsbl3NljLIkmsaopKQklC9fHnK5XGefbe1/5/U1LufKGGVJeB5leeNyrpYZo2SSlUS99PR0xMfHG3saRBbH398f9vb2xp6G2WOMItIPxijdYIwi0g/GKN1gjCLSD8Yo3WCMItIPXccoq0mSKxQKZGZmwsbGBjKZzNjTITJ7kiRBoVDA1taWNdN0gDGKSLcYo3SLMYpItxijdIsxiki3GKN0izGKSLf0FaOsJklORERERERERERERPQ63hIkIiIiIiIiIiIiIqvFJDkRERERERERERERWS0myYmIiIiIiIiIiIjIajFJTkRERERERERERERWi0lyIiIiIiIiIiIiIrJaTJITERERERERERERkdVikpyIiIiIiIiIiIiIrBaT5ERERERERERERERktZgk1yFJknL93tTGVHr48CEkSdLpuPqcr64p56evn1+XzOn3SqaLMcq8/i0xRpG1YYwyr39LjFFkbRijzOvfEmMUWRvGKPP6t8QYZZqYJNchhUKh+l4mk+nkL096erramNk/oyi2b9+O8PBwvHz5EjKZTCdjAkBGRobqe139DrZu3Yq5c+cWeZzXPXz4EIDu5gkAWVlZao91NW5aWppqLF3N98WLF0UeI/vfTzJ9jFGMUYxRZMoYoxijGKPIlDFGMUYxRpEpY4xijGKMKjpbnY9opfbt24cjR44gPT0dJUqUwCeffAIXF5cij7lnzx7IZDJUqVIFH3/8MWxsin5fY9OmTZg2bRrs7Oxw/vx5NGnSpMhjAsDevXuxe/duODg4oHLlyvjwww+LHPA2btyI6dOnY/ny5TqZo9LFixcxceJETJgwAY0bN1b9Qy/KfA8ePIgjR45AkiTUrl0bgwcP1knA37t3L3bu3In09HSUKVMG0dHRsLe3L9J8L126hO3bt+Odd95BlSpVtH7/7t27cfz4cfz1119o27Yt6tSpgxYtWhRqLmQYjFGMUYxRLQo1FzIMxijGKMaoFoWaCxkGYxRjFGNUi0LNhQyDMYoxijGqRaHm8jquJNeBQ4cOYeLEiahRowZq1KiBy5cvIywsDEePHkVaWlqhxty/fz8mTZqE2rVro0KFCjh+/DhWrVqler6wd242bdqEGTNmYN68eQgODsalS5cKNc7rDh06hMjISPj6+sLJyQnJyclITk4u0nw3btyI6OhoLFy4EC1btszxfFHuXj19+hRXrlzB/Pnz8csvvwAo2h2xH374AePHj4erqytSUlKwf/9+HDp0qMhz3b9/P6ZOnYr69eujadOmuHnzJnr27ImsrKwiBbxbt25h69at2LhxI27cuKHVew8dOoQpU6bA19cXHTt2xNWrV7Fs2TJ8/fXXhZ4P6RdjFGMUYxRjlCljjGKMYoxijDJljFGMUYxRjFGmjDGKMYoxSocxSqJCUygUUlZWlhQZGSlNnz5d7bn33ntPevPNN6W9e/dKqampWo2Znp4uffTRR9LixYslSZJUj2NiYqSsrKxCz/frr7+W/Pz8pIMHD0qSJEnLli2TRo4cKWVkZEiZmZmFGlOhUEiSJElRUVHS559/nuO5Z8+e5XitJq5fvy7VqlVL+vrrr9WOP3r0SHrw4EGh5ppdWlqa1L17d2n48OFS3759pSNHjqie0+Z3rFAopOTkZGnQoEFSbGysJEmS9N9//0k9evSQNm/eLL148aJQ81MoFFJaWpo0atQoaenSparjiYmJUvPmzaW4uLhCjytJknT58mWpTp06UteuXaXo6Gjpxo0bGo8xd+5cacKECarH//77rxQVFSUFBQUVel6kH4xRjFGMUYxRpowxijGKMYoxypQxRjFGMUYxRpkyxijGKMYo3ccoriQvAplMBhsbG2RkZODJkycAXtWBGjx4MB4+fIjVq1cjPj5e7bmCxrSzs4OtrS2ePXuG9PR02NnZ4f79+4iNjUX37t0xZMgQ/PPPP1rNNTU1FceOHcPcuXPRtm1bAICXlxdOnz6NR48eQS6XF+rukvLuUWZmJu7cuaM6PnPmTPTp0wdDhgzBmDFjkJaWplUNK3d3d7z99ts4c+YMHj9+DACYMWMG3n33XfTt2xf9+/fHzZs3ARTurpi9vT1cXFxQsWJF+Pj4YPny5Th69CgAICUlReNxZDIZbG1tkZaWBicnJwCAi4sL0tLS8NVXX6FLly7o06cP/v77b63mqvx7kJmZqXaHrnz58ihfvrxqPG0px7px4wY8PDzQq1cvnDlzBhs3blT9PvMjSRJu3ryJBw8eqI5VrlwZ7dq1Q2pqKtavX4+4uLhCzY10jzGKMYoxijHKlDFGMUYxRjFGmTLGKMYoxijGKFPGGMUYxRil+xjFJLkOVKxYEefOncPdu3dVNZrKlSuHnj17oly5cpgwYQIkSdKqflPVqlXx22+/4dNPP0Xfvn1x7949fPrpp3jvvfeQmZmJIUOG4NGjRwA0+4vu6OiIlStXIjQ0FJmZmQCA9u3bo3bt2li7dm2Rt0t4eHjg77//hkKhwPr163H69Gn07dsXvXr1wq1bt9CzZ08oFArY2NhoNN/ixYsjNDQUMpkM27dvx+zZs3Hp0iUMGjQIkZGRsLOzQ3h4OJ4+far1thRlM4eQkBDUqlULQ4YMgbe3N9asWYOxY8di3LhxSE9P13hMOzs71KxZUzXuhg0boFAo8Mknn2DixIlwdHTEu+++i+fPn2s1V2XA++2336BQKJCZmQm5XA5fX1/Vf/vCNs6Qy+WoUKEC+vXrh969e+OPP/7QKDDJZDJ06tQJiYmJatt3qlevjqZNm6Jbt2745ZdfcOHChULNi/SDMYoxijGKMcqUMUYxRjFGMUaZMsYoxijGKMYoU8YYxRjFGKXDGFXktehW6NKlS9K5c+ekn3/+WXWsX79+UseOHaVt27ZJJ06ckFq1aiXNmjVLkiRJ6tatm7Rq1Sqtx1y6dKk0ceJEqUmTJtKuXbvUXv/2229L48aN03iuv/zyi9rxzMxMSaFQSIsXL5b69Omj2oah6daOv//+W7px44Z09uxZ1bFOnTpJ7733nrRs2TJp9+7dquMPHz6UOnXqJM2ePVvrMb/55hupZ8+e0sCBA6Vff/1V7fU9e/aUpkyZovFcz507p3Z8x44dUrdu3SRJkqTbt29LPXv2lPz8/KSFCxeqXpPXlhzlmH/88YfqmHKLUFpamvTPP/+ovb5bt27SjBkzNJ7rn3/+KUmS+O/x66+/SikpKarXREZGSv3791ebX/bnX5fb34H09HTp0qVLqsfbtm2TevbsKc2aNavArS43b96UPv30U2ngwIHSokWLpBMnTkjt27eXpkyZIt2/f1/q0aOHtGnTpgJ/VtIfxijGKMYoxihTxhjFGMUYxRhlyhijGKMYoxijTBljFGMUY5T+YpRt4dPr1mnr1q1YtGgRKlWqhGvXriEgIADvvPMOvv76a0ydOhWxsbF48eIF6tati8jISACAp6cn7t27p9WYAwYMwMiRIyGTyTBz5kzV1on09HTY29ujWbNmBd4deX3cwMBAdO3aFe3bt4dcLgcADBs2DFu3bsWSJUswYcIE1Z21/O7ibd26FStWrEDJkiVx+/ZtvPHGGxg1ahSmTJmCL774AosXL8bnn38OQNxVLFu2LOrVq6e2JaKgMRs0aICIiAgMGjQI8fHx2L17NyZPngwAyMrKglwuxxtvvIG7d+8W+DvIPm6jRo0wdOhQ+Pr6IigoCJs3bwYAHDlyBDdv3kTjxo1x+vRpHDp0CG3atMn19/D6mA0bNsSwYcPg7+8PQGydUXbozczMhK2tLYKDg1XbdLT5HYwYMQINGzaEXC5X/dzAq20qMpkMW7Zswa5du7B69Wo4OjrmGPP1vwPdunVD+/bt4ePjg4yMDNjZ2aFbt24AgO+++w42Njbo2bMnqlevnuvPX7VqVQwePBjHjh3Dt99+ix9//BE1atTAtGnTIJfLUatWLRw+fBi9evXSSTdl0g5jFGMUYxRjlCljjGKMYoxijDJljFGMUYxRjFGmjDGKMYoxSs8xqtDpdSuUkJAgtWnTRtq9e7f04sUL6fHjx9KIESOkt956S9qxY4ckSZL08uVL6cmTJ2rvmzlzpjRnzhxJknLeCcprzLCwMGndunWSQqGQZs6cKXXp0kXtfQsXLpQ++OADKS0tLde7S3mN27lzZ2ndunVqjREOHDggdenSRdq3b1+Bv4MrV65ILVq0kPbv3y/du3dPSkhIkMLCwqSePXtKe/fulfbs2SOFhoZKffr0kZKTk1XvmzNnjjRu3DgpIyMjx3zzGrNz587Sb7/9JkmSpHaXSWn27NnSxx9/LGVkZOR6xzGvcbt166a6ixURESH17t1beuONN6RTp05JV69elcaNGyf16NFDevTokcZz7datm3T48GHV6//77z+190VHR0sTJkxQ3THVZq7Zx5UkSZoyZYo0fPhwSZLEHTdfX19py5YtOcbU9O9ARkaG6j3bt2+XQkNDpZkzZ6r995MkcRfx9bn/999/Ulpamtqxzz77TIqOjs4xH9I/xijGKMYoxihTxhjFGMUYxRhlyhijGKMYoxijTBljFGMUY5T+YxRXkmvh2bNnSE5Oho+PD4oXL47ixYtj/vz56NKlCxYuXIjnz5+jT58+kMvlWLx4Mf777z+4uLhg8+bNWLFiBQDkuJOR35gxMTGwt7fHsP9r787DY7r3P4C/EyFEQoRSzYKmtaUoP3tsQUNzrbUv7S33tm61XL12UrvaKog011ZLqJBNVC2p2reH2BK1REIlpJYKIRLZZj6/P3JzmpEdwzmZ9+t5+jzM8pn3zHzPW5+ZOef88584duwYevfujQ4dOkCv12PTpk3w8vJCmTJlip117dq10Ov1GDx4MEqVKoVGjRqhadOmCA0NRZUqVdC0adN8X4P79+8DAN5//31Uq1YNAPDjjz+iY8eO2Lx5MwYOHIiJEydi7ty56NOnD5o1a4Zy5cph06ZN8PHxgYVF7iWX18zNmzejU6dOWLRoEb766it06NABDx48wNy5c9GkSRMkJiZi8+bN8Pb2znNmYXOXLVuG5ORkWFtb4+LFi1iyZAmaNWsGIOskF6VLl4adnV2xZq5cuRJ6vR7t2rXDgQMHsHnzZrRo0QJA1nGhvL29lW/enmdux44dYWZmhjJlyiAzMxPbtm3DlClTsGjRInTr1i3XN65FWQNDhgyBhYWF8q1gr169oNfrcf/+fSxZsgRWVlaoWbMmevXqpRzDLPtxREQ5KcSCBQtw79491KxZE0FBQfD29s7zeZJxsaPYUewodpSasaPYUewodpSasaPYUewodpSasaPYUewo43cUT9xZDHZ2drCxsVHODqzT6VCuXDm0a9cODg4OCA8Px7lz52BpaYnSpUsjMTER0dHRWLJkCVq3bp3nAfILmmlvb48jR44gLi4Ofn5+qFmzJs6fP48bN27Ay8sL7du3z/eg+wXNfeutt3D69GmcPn0aQNaZez/88ENlIy1I+fLlUaZMGdy8eRNA1u42FSpUQMuWLZGYmIj9+/fD2dkZP/30E1xdXZGSkoLHjx/j+++/R4cOHfLMm9dMGxsbNG/eHOnp6QgJCUF8fDzS0tIQGxuL7du34+zZs1iyZEmBr0FBc9PS0nD06FF0794dvr6+cHV1Vea89957qFOnTrFnpqSkYPv27bh//z46duyIypUrIyoqCjdu3IC3t/dzZ82ee+vWLQBA1apVcfz4cUyePBkLFy5UCulZxVkD5ubmykkXLC0tsXz5ciQlJSEqKgre3t4YOXIkEhMTARjuWgNk7cLTqlUrPHjwAPfu3cPixYvRrl275zrLM70YdhQ7ih3FjlIzdhQ7ih3FjlIzdhQ7ih3FjlIzdhQ7ih31CjrquX+DboKSk5Nl5MiRMmLECDl16pRy+cCBA2Xr1q0yfPhw+eqrrwzuk/3Tf71en+duDYXNHDZsmIwaNcrgPunp6QXOLGrW0aNHG9zn2d1y8vLw4UPp0aOHjBo1Sm7evKlc/sknn8jOnTvF3d1dPD09De6TvetJfnkLmrlr1y5xd3eXadOmKZenpaUpu14U9BoUlrVz585FOnlBcbM++/yL8n4VZW72SSHi4+OlW7dusnfv3gLnFncN6PV6SU1NlVGjRikn9khPT5eYmBjp3LmzDBgwQOLi4vJ9bXLuZlTQcyXjYUexo9hR7Cg1Y0exo9hR7Cg1Y0exo9hR7Cg1Y0exo9hRxu8oMxF+BVgU8r+f9MfHx2PSpElITU1FRkYG7t+/j1q1amHjxo34448/8NFHH2Ht2rWoX7++wf1edOYPP/wAFxeXQmcWd+769etRp04dZXeFgubq9XqYm5vj999/x/Dhw1GtWjWUL18eMTExcHR0xKZNm3D+/HlMnDgRa9asgYODA8zMzJT7vejM1atXw8nJqUivQXHmrl27Fvb29vnOetHn/zKz/vDDD6hcuTIeP36MatWqKd+OPTu/uGugbt26yn1HjRqFN998E1OnTlUuS0lJQd++fVGxYkWsWrUKNjY2yi4xV65cQaVKlZTdcuj1YEexo9hR7Cg1Y0exo9hR7Cg1Y0exo9hR7Cg1Y0exo9hRr6ajeEzyIsh+g3U6Hezt7bFs2TKcOXMG165dQ+XKldGvXz8AwNOnT5VdCrIVVkhFnVmhQoVCZz7P3PLly+faXSE/5ubm0Ol0qFWrFvz9/bFnzx7cvXsXrVu3xj/+8Q8AQEZGBkQEVlZWyrz8Cqm4M8uXL1+k16C4c589A+/Lfv4vO2u5cuWUs0vnNft51kDO+9rb2+Pq1avKmavT09NhZWWFoKAg9OjRA5MmTcL333+PUqVK4eLFi/jXv/6FNWvW8H+cXiN2VBZ2FDuKHaVO7Kgs7Ch2FDtKndhRWdhR7Ch2lDqxo7Kwo9hRr6SjCvqZuam6ffu2XL9+XVJTU5Wf62efefXZs9ampaVJeHi4JCYmys6dO+WDDz7IcxcAY8w05twHDx7I3bt3DS7Lnvfs3IyMDImJiZHU1FT5+eefpUePHnLnzp1XMpNZX/4aePjwobRu3VomT55scD8RkejoaHFzc5MdO3Yo112+fDnP50rGw47idq+lrOwo08OO4navpazsKNPDjuJ2r6Ws7CjTw47idq+lrCWto/hL8mcEBgZizZo1SEpKwrvvvou6detizJgxKFeuHDIzM5Uz1+r1epiZmSEsLAwbN27Eo0ePkJycjClTpsDR0dHoM405Nzg4GBs3bsSDBw/g4uKC1q1bo3///rC0tFR2acg5d9OmTdi9ezfS0tJw584dTJs2Lde3OMaYyawvfw3odDrY2trC29sbn332GaysrODp6amctdrJyQm1atVCXFyccp+cu8aQ8bGjuN1rKSs7yvSwo7jdaykrO8r0sKO43WspKzvK9LCjuN1rKWtJ7CgekzyHCxcuYMSIEZg6dSqcnJwQGRmJvXv3IikpCRs2bIC1tbXBGw0AqampuH37Nu7evYs33ngDzs7OBsf8McZMY849ffo0vvjiC0ycOBFvvfUWDh06hKtXr8LGxgbz58+HlZWVwQYEAI8fP8aVK1eQmJgIe3t7uLi4GMw1xkxmNd4ayLZv3z6MGzcOH374IWbMmKEU0/jx4/Hmm29i7Nixue5DxsWO4navpazsKNPDjuJ2r6Ws7CjTw47idq+lrOwo08OO4navpawltqNe6u/SNe7UqVPSq1cvefDggYhk7Rpw8eJFGTJkiHzwwQeSlJQkIn/tOvC6Zhpz7t69e2Xo0KGSkZEhIllnk/3ll19k6NChMmzYMOUsvsWZa4yZzGq8NZBTeHi4uLq6yvDhw2X69Oni5eUl7733nhw/fvy5Z9LzY0dxu9dSVnaU6WFHcbvXUlZ2lOlhR3G711JWdpTpYUdxu9dS1pLaUfkfwd4Epaam4t69e8q3J+bm5qhfvz6+/fZbVKtWDcOHD8fTp0+V67dv345Vq1ZBp9O90pnGnPvkyRPExMQo3/aULl0aHTt2xMiRI5GZmYkpU6YgLS1NmfvTTz/h559/fuUzmdV4ayCnpk2bYufOnahduzYSEhJw/fp1LF26FK1atVLOZEyvDjuK272WsrKjTA87itu9lrKyo0wPO4rbvZaysqNMDzuK272WspbUjuKH5Dm0bdsWNWvWxJgxYwxecCcnJ0yaNAnlypWDt7e3ct29e/dQvXp1g10SXsVMY8718PCAg4MDPD09odfrAQClSpVC8+bNMXToUDx8+BBBQUEAgPT0dERGRhY4z1gzmfXF1kB2hsLodDpUrFgREyZMwPLly7F48WJ06tSJ/9P0mrCjuN1rKSs7yvSwo7jdaykrO8r0sKO43WspKzvK9LCjuN1rKWuJ7aiX/dN0rdLr9SKStcvAoEGD5JtvvlEuy75+3bp18ve//105s+qz930VM405VyRrF4mQkBD55JNP5L///a/BdRkZGTJ79mz54osvlMvS09MLnWuMmaae9XnWwLFjx2TRokXK7Z4903BRFPY8yXjYUVlMebvXUlZ2lOlhR2Ux5e1eS1nZUaaHHZXFlLd7LWVlR5kedlQWU97utZS1JHcUf0n+P9kHim/UqBH69OmDa9euYfr06co3HGZmZmjUqBGuX7+Ou3fv5nnfVzHTmHOBrF0kPvjgAzRp0gQnTpzAihUrlN0hLCws0LZtW1y9ehUJCQkAsnbTKGyuMWaaetbnWQPJyclYs2YN/Pz8lExF/Qbv2celV48dlcWUt3stZWVHmR52VBZT3u61lJUdZXrYUVlMebvXUlZ2lOlhR2Ux5e1eS1lLckfxQ/IcRARlypSBh4cHevfujZiYGPzzn/9EUlISEhISEB0dDRsbG1hZWb3Wmcaea21tjWHDhuH999/HiRMnMGPGDABAWloaYmNjYWtrq5xZ9nXNZNbirwFLS0tYWFhg1apVWLZsGYDCi6k4x4si42NHcbvXWlZ2lGlhR3G711pWdpRpYUdxu9daVnaUaWFHcbvXWtYS2VFG+X26hmX/fD8tLU3Cw8Nl6NCh0rJlS/Hw8JDmzZvLzp07VTHTmHOzd3tISkqSbdu2Sf/+/aV58+bSr18/ef/992XXrl155niZM42V9XXNNNbc4qyB9evXS6dOnWTDhg3SuXNn8fb2zpUtW2RkZL7X0evFjmJHaSkrO8r0sKPYUVrKyo4yPewodpSWsrKjTA87ih2lpawlsaPMREzrrAx6vR7m5gX/gF5EDH7Gf+zYMVhaWsLa2hp169bNdX1RvMjM4lxXlLmZmZnKWW0LyysiyMzMxO7du2FtbY0qVaqgYcOGBnNv3rwJR0fHIr8GRZlZHMWZGxERAWtrazg7O6s+a3FnZstrDQDA6tWrcf/+fXz11Vf46aefsG7dOvTs2ROjR48G8Ne2ERcXh169eqFHjx7Kt4v06rCj2FHsKHaUmrGj2FHsKHaUmrGj2FHsKHaUmrGj2FHsKJV3lJiIq1evSmZmpogU/ZuIwm53+fJlOXfunJw5c6bIOYry2Ddu3JC4uDi5du2aiBTvm7GCHDx4UDlofkZGRqG3L0h2po0bN0qdOnXkzp07LzQv50wRkejoaLl69apER0e/tLkhISFSp04d2bBhw0ubKSISGxsrt27dktjY2FzXPe/c+Ph4uX//vsTFxeX7uM8qbA1cv35djh07JiIiiYmJ8uOPP0rHjh1l2bJlBrdLTk6WO3fuyDfffCPbtm17jmdBz4MdxY5iR7Gj1IwdxY5iR7Gj1IwdxY5iR7Gj1IwdxY5iR2mjo0ziQ/I///xTPv/8c1m5cqXyxr3oT/YDAwOlTZs20rt3b6lfv7588803uRbQ8wgICJDOnTtLt27dpHnz5rJ79+4XnimStfuDm5ub9OjRQzlT7YsWk7+/vzRo0OClZcwWEBAgrq6u8re//U1cXFxkz549Lzxz8+bN4uLiIgMGDJBp06aJyMvZbSMoKEjc3d3Fw8ND2rZt+9y7vjw708PDQ3r16iUtW7YULy8vuX79unL985ReXvfJWUzLly/Pddvw8PCX8o8NFY4dxY5iR7Gj1IwdxY5iR7Gj1IwdxY5iR7Gj1IwdxY5iR2mno0ziQ/InT55Ily5dpEuXLrJy5cpif4P3rIiICHF1dZWwsDB5+PChXLhwQZo1ayZLly59oZynTp2S5s2bS1hYmERERMiCBQukX79+kpqaanC751mUOp1ORo0aJXXq1JFevXq98Dd4QUFB4uLiIgcPHizwdsXNGhMTI66urrJr1y6JiYmRM2fOyOPHj58rY7aAgACpV6+eHDhwQEJDQ6VDhw7K838Rp0+flpYtW0pYWJicPn1aFi1aJB4eHpKSkvLcMyMjI6V169aya9cuiY6OlkOHDomrq6uMGDFCzp49q9yuqK9rYbdLTEwUf39/adGihSxZssTgOh6f7tVhR7Gj2FF5Y0epAzuKHcWOyhs7Sh3YUewodlTe2FHqwI5iR7Gj8qbGjir4QEAlRPny5VGuXDkAwPHjxyEi+Oyzz5QzqT57TKjsy3Q6HUqVKpVrXlxcHKpXr45OnTrB3Nwctra2GDlyJNatW4eBAweiWrVque5T2EwAuHbtGmrXrg13d3cAwK1bt3D27FkcOHAAZmZmaNKkCd544w2YmZkp84oyV0Rgbm6Ojh07wsLCAomJiejevTtCQ0OV1yWv++R8nGft3LkT9vb2aN++vXLZ3r178ejRI1SoUAFNmzaFnZ2dwYzCZgJAYmIirKys0LJlS1SqVAlA1jGNHj16BDs7O7Rs2TJXxoLmbtmyBTNmzICPjw86dOiAO3fuwNbWFhcuXMD//d//5XvMpaJkjY2NhbOzs/J+PXnyBBcvXkRUVBTMzc3RsGHDYmUFso6nZWdnBzc3N5QtWxbvvPMOVq5cib59+8LMzAz/+Mc/0LRpU4PM2bNOnjyJGzduoFKlSnByckLdunUNHjMvFStWRNeuXaHT6XIdF6uwY6XRy8OOYkexo9hRasaOYkexo9hRasaOYkexo9hRasaOYkexo7TTUSbRjLGxsUhJScGUKVPQpEkTHDt2DGvWrIFOp1M27GxHjx7FtGnT8OTJk1wbufzvgPM2NjZ4/Pgx4uLilDe7Ro0aePz4MTIzM3M9fkEzc861sLCAhYUFHjx4AADYtWsXfvvtN6xZswbjx4/H7NmzcfbsWQBZC+bYsWMFzs2WndHGxgZRUVGYMGECrK2t0adPHwDA6dOnce3aNYO8K1asQEpKSq6FmZ31u+++g06nw4IFCwAAU6ZMwcKFC7F582ZMmDABs2bNwpkzZ5SsBc3Myc7ODhkZGYiOjgYAjB8/HjNmzMB///tfDB8+HNOmTcPNmzeV51XYXEdHR/j4+KBz584AgAoVKiAjIwO7d+82eG1yKmxm9mtga2uLtLQ0HD58GAAQEBCAs2fPYuLEiRg8eHCxsmbPtLKyQmpqqnK/9PR0uLi4oEmTJrh16xZ27NiBhw8fKvfLXgP+/v7417/+hR07dsDHxwfjxo3DunXr8n2OOdna2mLAgAFo2bKlkoNeLXYUO4odlT921OvHjmJHsaPyx456/dhR7Ch2VP7YUa8fO4odxY7Kn9o6yiQ+JK9evTqcnZ1RtWpVDBs2DI0aNcLRo0cNikmn00Gv1yMyMhJly5bFqlWroNPpDOZkv7kODg6wt7fHkydPlNu88847MDc3N1gwQNZiu3jxYr4zc85t3rw53N3dYW1tjbt37+Ltt9/Gjh07EBQUhG3btuHevXsICQkBAOh0OkRERBQ491lNmzaFvb093nzzTfj4+KBixYpo0aIFRowYAWtra+h0Ouh0OuzevRsJCQkIDQ3NN6udnR2mTp2Ks2fPYty4cXjw4AFWrVqFkJAQBAcH486dOwgODgYA5ay5+c3MydbWFg4ODvjll18QGRmJe/fuYcWKFQgKCsKWLVtw8OBB+Pn5FWmuTqeDq6urUkiZmZmwsrLCqFGjEB4ejkuXLuV5n8KyZr8G7777LipWrIhvvvkGH3zwAS5cuICQkBAEBgZi/fr1OHDgADZt2lSkrNkza9WqBZ1Oh+DgYDx48ABlypQBkPXt88cff4w9e/Zg+/btStbz588DAL7//ntMnToVmzZtwsqVK/HJJ59g8eLFyj8aAPIsnJz/IObMQa8WOyoLO4od9Sx2lDqwo7Kwo9hRz2JHqQM7Kgs7ih31LHaUOrCjsrCj2FHPUmVHFXw0Fu0JCwuTefPmyahRo+TIkSPKyQtyHvMnKSlJvvvuO/n4449lxYoVyjGhsiUmJoqvr69yn5wzDx06JKmpqZKWlmZwbKaYmBhp2LChXLx4Ubns3LlzcvPmzTxn5jX34cOHBjmSkpJE5K9j8Zw7d04aN24sMTEx+WbN6zWIj49Xruvdu7eEhoaKiMj8+fOlQYMG0rZtW4PH1ev1cvnyZfH19VUeO+fMw4cPy927dyUjI0OWLl0q7dq1kylTpoher1duHxERIQ0bNpQrV67kOzOvuampqXL+/HmpX7++jBo1SsaMGWPw/hw7dkwaNGignGW4sKzPPn+RrDMV9+nTR/z9/UVEcr3/Rcl66NAhSUlJkXv37sn+/ftl7NixMm7cOIM5x48fl8aNGxucFbqgrIcOHZKMjAwJDw+XZs2ayfDhw2XChAnSvn17+fTTT0VEJDQ0VHr06CEJCQnK4/z555/Stm1bCQ4ONnj8gwcPSsOGDWXhwoXKZdmPe+fOHeWEGfRqsaPYUeyoLOwodWJHsaPYUVnYUerEjmJHsaOysKPUiR3FjmJHZdFqR5WoX5Lv2rULEyZMwNOnT5GZmYl58+bBx8cHJ06cUL4BSU9Ph7W1NUaMGIFGjRrh5MmT8PX1Nfjmq1y5cvjoo49QpkyZXDPnz5+Pb775BmfPnoWlpaVyn+TkZKSlpSnfgISGhmLgwIFISEjINTOvrAsWLMC8efNw4sQJZaa1tTWAv47FY2trC3t7e2VGUebOmzcPy5Ytw5EjRwAAbdq0gYggMDAQISEhmDFjBhwdHdGuXTukpaUByPr2platWujZsyfMzc3zfA2+++47XLp0CUOGDEHNmjXh7u4OMzMzg6xOTk4oW7ZsnjPzyzp16lSUL18evr6++OWXX3D58mXcu3dPeb7Vq1eHo6Oj8nwLy5r9/HO+ru+88w48PDywYMECxMXFoVSpUga7OBUl6/z58zFt2jTExsbCzc0Nffr0QalSpZR1JCKws7ODvb09rKysivy6Tpo0CVWqVEFgYCBcXFyQmZmJ3r17K7uqlC1bFnq9HuXLl1fyWlhYoGrVqsouQdmP3759e3h7e8PPzw/r169X1tLhw4fRv39/ZTcqenXYUewodhSUx2dHqQ87ih3FjoLy+Owo9WFHsaPYUVAenx2lPuwodhQ7Csrja7ajXsMH80YzefJk8fLyUv5+6NAhGT9+vAwZMsTgrLfZ39QkJSXJ3LlzZciQIRIbGysiuc++WtSZly5dkvfee0/i4+Pl559/FhcXF+Ubsme/GSrq3EePHkl4eLjExMRIRkaGBAcHi5ubm/JtZF5nii1obnh4uBw/flzq1KkjjRs3lmPHjomIyPXr16V3795y8uTJIr8G48aNk0GDBsn58+eVyx89eiR3794VEZGQkBCDrDm/ASto7tixY2Xw4MFy/PhxOXjwoNSpU0emTZsmkZGR8vTpUwkODpZ27dop79ezc4v6fomI/Oc//5HPP/9c7t27Z3B5UbOOGzdOhgwZIkeOHJFbt25JvXr1ZN26dZKcnCwiItu2bZMPP/xQ+eawKFnHjh0rgwYNklOnThnkuX37tvK6DhgwQB48eCB6vV55r44ePSr16tWTkJAQ5X7Z1/n7+0u3bt0kKipKue7IkSO5niMZHzuKHcWOYkepGTuKHcWOYkepGTuKHcWOYkepGTuKHcWO0n5HWbzuD+lfFp1Oh/j4eOVbEgBo164dbG1tERgYiB9++AGWlpZo2bKl8k2NtbU1Ro8ejXv37sHJyQmA4TFwijOzUqVKsLOzw9KlS7Fjxw4sXLgQ3bt3V87kW9ysVlZWqFevHubMmYPExETY2dnhzp078PT0hKOjY66shc0NCAiAt7c3unfvjqFDh8Ld3R3NmzcHkHXcoR9//BHlypXLdebZgmZu3boVixcvxogRI+Dq6orPP/8cCQkJqFq1KmJiYjB9+nQla3Feg4CAAPj6+sLT0xMbN27ErFmzcPz4cZQtWxYJCQnw9PRU3q+cc4vzfgFA37594evrCx8fH3h6eqJ06dLPlXXlypWYPn06vLy88PXXXyM0NBQ2Nja4du0apk2bhrfeeqtYWQMCArB8+XLldV29ejX27dsHS0tLXL58GXPmzFHOtJytdevWGDNmDLy8vGBpaQkPDw/lOldXV2zcuBHx8fGoXbs2gKxvcOnVYkexo9hR7Cg1Y0exo9hR7Cg1Y0exo9hR7Cg1Y0exo9hRJaSjXten88YQEBAgHh4ecuHCBYPLIyIiZNy4cTJx4kT5448/lMvz+obmeWbGx8dLQkKCNGvWTOrUqSNhYWEiIgbfrBR37oQJEyQjI0Oio6Pl8OHDsmfPHrl8+bIy93nyTpgwQTw9PZU5RVXQzPHjx8vEiRMlLS1Nzp07J/7+/uLv7y8REREvlDX7mEoZGRly+/ZtiYiIkOPHjyvHv3re1zX7/cq2ceNGg28fnzfr+PHjJTMzUyIiIsTPz0+2bNlS6GtQlKyPHz+WuLg4CQsLky1btsjZs2fznfnkyRPx8vISV1dXg2/wREQGDx4s27ZtK/R5knGxo9hR7Ch2lJqxo9hR7Ch2lJqxo9hR7Ch2lJqxo9hR7Cjtd1SJ+pA8KipKvvzyS5k5c6bcuHHD4LqDBw9Kly5d5Pjx4y995tGjR0VE5KeffpL9+/eLSMGFVJS57u7uytyXlXf//v3StWvXl/4aGCNr9tzsXXBe1sznWQNFzVrc16A4a+tZ+a2tpKQkWbFihdStW1cWLFggQUFBsmHDBnn//fclMjKyWPno5WNHsaPYUewoNWNHsaPYUewoNWNHsaPYUewoNWNHsaPYUdrvqBJ14s7atWuja9euuHLlCgICAhAbG6tc1759ezg5OWHbtm1GmSki6N69O9zc3CAiAHLvflKcuTVq1EBoaCgAKPNeNK+bmxscHR1f+mtgjKzZc431fhkja3Ffg6JkzW9mfmsr+yQga9euRUREBPz9/REaGop58+ahQYMGRX6+ZBzsKHYUO4odpWbsKHYUO4odpWbsKHYUO4odpWbsKHYUO0r7HVUijkmeswS6deuG5ORkhIaGIiUlBX379oWLiwsAoHLlyqhatapRZuZcJAWV0fNkLWieml4DZi3aXGNlzalVq1Zo2LAhLC0tkZSUhEqVKhXpH0syDjWseXYUs7KjKD9qWPPsKGZlR1F+1LDm2VHMyo6i/KhhzbOjmJUd9XKYSXG+vlCB8PBwXL9+HZUqVYKTkxPq1q0LIOuNlhwnJdixYwf27t2L33//HY0bN0bp0qURFBQEHx8ftG3b1ugztTaXWbWVtTjkmZNfkHFpaR1paS6zaitrcbCjXi0trSMtzWVWbWUtDnbUq6WldaSlucyqrazFwY56tbS0jrQ0l1m1lbU4tNZRmvqQfNu2bZg1axZcXFzw+PFj6PV69OnTB8OGDVNuo9frlTc6JiYG58+fx65du/Dmm2+ic+fO6Nixo8GbZIyZWpvLrNrKSuqlpXWkpbnMqq2spF5aWkdamsus2spK6qWldaSlucyqraykXlpaR1qay6zaylriiQbo9XpJSEiQjz76SIKDg0VE5I8//pCtW7eKi4uLzJ8/3+C2z54lOD093eD67P9e9kxjZeVrwKykblpaR1qay6zaykrqpaV1pKW5zKqtrKReWlpHWprLrNrKSuqlpXWkpbnMqq2spkITxyQ3MzODtbU1zM3NYWGRFbl69ero378/qlWrhtGjR8Pc3Bzjx4+HmZkZ5H8/jr979y4qVaqEMmXKGMzKZoyZxsrK14BZSb20tI7YUcxqrKykXlpaR+woZjVWVlIvLa0jdhSzGisrqZeW1hE7ilmNldVUmL/uAEWVkZEBCwsLXL58WblMRNC+fXt4e3vDz88P69evBwCYm5vj8OHD6N+/Px4+fPhKZ2ptLrNqKyupl5bWkZbmMqu2spJ6aWkdaWkus2orK6mXltaRluYyq7ayknppaR1paS6zaiurSSjqT87V4OjRo1KvXj0JCQlRLsv+6b+/v79069ZNoqKilOuOHDnyWmZqbS6zaisrqZeW1pGW5jKrtrKSemlpHWlpLrNqKyupl5bWkZbmMqu2spJ6aWkdaWkus2ora0mnqQ/J9Xq9rFy5Utq0aSM7d+5ULhMRiYuLEw8PD9m/f/9rn6m1ucyqraykXlpaR1qay6zaykrqpaV1pKW5zKqtrKReWlpHWprLrNrKSuqlpXWkpbnMqq2sJZ0mjkmezczMDEOGDEFycjK+/fZbpKWloXfv3gAAR0dH2Nra4tGjR699ptbmMqu2spJ6aWkdaWkus2orK6mXltaRluYyq7ayknppaR1paS6zaisrqZeW1pGW5jKrtrKWdJr6kBwAypcvj88++wxWVlaYMmUKoqOj4ezsjOTkZFy6dAnOzs6qmKm1ucyqraykXlpaR1qay6zaykrqpaV1pKW5zKqtrKReWlpHWprLrNrKSuqlpXWkpbnMqq2sJZmZyP9OZapBJ06cgI+PD9LS0qDX6/H555+ja9euqpuptbnMqq2spF5aWkdamsus2spK6qWldaSlucyqraykXlpaR1qay6zaykrqpaV1pKW5zKqtrCWNpj8kB4Dk5GRYWloiKSkJlSpVQvbTMTMzU9VMrc1lVm1lJfXS0jrS0lxm1VZWUi8trSMtzWVWbWUl9dLSOtLSXGbVVlZSLy2tIy3NZVZtZS1JNP8heTYReelvrDFmam0us2orK6mXltaRluYyq7ayknppaR1paS6zaisrqZeW1pGW5jKrtrKSemlpHWlpLrNqK2tJUGI+JCciIiIiIiIiIiIiKi7z1x2AiIiIiIiIiIiIiOh14YfkRERERERERERERGSy+CE5EREREREREREREZksfkhORERERERERERERCaLH5ITERERERERERERkcnih+REREREREREREREZLL4IbmGdezYEXXq1FH+q1u3Lpo3b44vvvgCt2/fNvrjT5o0CZMmTQIALF++HB9//HGu25w8edIg47P/TZo0SbnNq/Cij/Xxxx9j+fLl+V5fp04dnDx58rnnE5Uk7KjiY0cRvTrsqOJjRxG9Ouyo4mNHEb067KjiY0epn8XrDkAvZsqUKfDw8AAA6PV6xMTEYPr06Zg4cSL8/PxeczqgcePGOHr0qPL3Nm3aYPny5WjcuDEAoGzZsrh06dLrikdERsaOIiI1Y0cRkZqxo4hIzdhRVNLwQ3KNs7GxwRtvvKH8vVq1ahg9ejTGjx+PpKQk2NjYvMZ0QJkyZQzyAUDFihVzXUZEJRM7iojUjB1FRGrGjiIiNWNHUUnDw62UQGXKlAEAmJtnvb2PHz/G+PHj0aRJE7Rp0wazZ89GamqqcvvIyEgMGjQIjRo1QpcuXbBz507lusDAQHTt2hXvvfceWrRogZkzZ0Kn0xklt7+/P9q2bYvGjRtj8uTJSE9PB5C168zIkSMxZMgQNG/eHKdOnUJ6ejrmzJmDFi1aoEWLFhg3bhwSExOVWX5+fnBzc0ODBg3w0Ucf4fTp00V6LAA4cOAAevfujYYNG8LDwwO//PJLvpl9fHzQqlUrtGjRAoGBgQbXnThxAj179kSDBg3QqVMnbNmy5SW8SkTax45iRxGpGTuKHUWkZuwodhSRmrGj2FFaxg/JS5i4uDisWrUKbdu2Rfny5QEAU6dORVJSEvz9/eHr64sLFy5g1qxZAICEhAQMHz4c9erVw7Zt2zBixAhMnDgRV65cwalTpzBnzhz85z//wZ49ezBz5kwEBQVh3759RskeFhaGH374AT4+PtizZw+Cg4OV6/bt24du3bphw4YNaNiwIby8vPDbb79h9erV8PPzw5MnT/Dvf/8bAHDp0iUsXLgQ06dPx+7du9G0aVOMGTMGer2+0Mc6ceIERo0ahZ49e2L79u3o168fvv76a/z222+58m7duhV+fn749ttvsX79eoO8Op0OY8aMQdeuXbF79278+9//xsyZMxETE2OU145IK9hR7CgiNWNHsaOI1IwdxY4iUjN2FDtK63i4FY2bPn06Zs+eDQDIzMxE6dKl0alTJ0yZMgVAVkn9+uuvOHXqlLKry+zZs9GrVy9MnjwZO3fuRMWKFeHp6Qlzc3O8/fbbePToEVJTU2FlZYW5c+fC3d0dAODg4IB169YhOjpauexlP5datWqhdu3aaN26Na5cuaJcV6VKFQwaNAgA8PTpU2zatAnBwcHKSQ8WLlyIFi1aICoqCvHx8TAzM8Nbb70FBwcHjBkzBm5ubgallN9j/fjjj+jSpQs+/fRTAECtWrUQGRmJtWvXwsvLyyBvQEAA/v73v8PNzQ0AMGfOHPztb38DACQlJSExMRFVqlSBg4MDHBwcULVqVe7WQyaHHcWOIlIzdhQ7ikjN2FHsKCI1Y0exo0oafkiucaNHj4a7uzuSk5OxfPlyxMfHY+zYsahUqRIA4Nq1a9Dr9WjXrp3B/fR6PWJjY/H777+jfv36yq4wADBs2DDlz2XLloW3tzdiYmIQFRWF2NhYtGnTxijPxcnJSfmzjY2NwS4n9vb2yp9v3ryJjIwMDBw40OD+er0eN27cQLt27VC7dm10794d9evXR6dOndCvXz9YWPy13PN7rGvXruWa27hxY4Nv5bJdu3YNX375pfL3d955B1ZWVgAAW1tbDBo0CJ6envD19YWbmxv69OmDihUrFus1IdI6dtRf2FFE6sOO+gs7ikh92FF/YUcRqQ876i/sqJKBH5JrXOXKlVGjRg0AwLJly9C3b1+MHDkSW7duRenSpaHT6WBjY5PnRlWtWjWDDfVZR44cwZdffolevXqhbdu2+PLLLzFz5kyjPZdSpUoZ/F1ElD9bWloqf84+BtXmzZuVEshWuXJllCtXDoGBgTh16hQOHDiAkJAQ+Pv7IyQkpNDHyvk42fR6vcG3fvllBGDwes6YMQNDhgzBr7/+il9//RVbt26Fr68v2rdvn+csopKIHcWOIlIzdhQ7ikjN2FHsKCI1Y0exo0oaHpO8BClTpgzmzJmDy5cvY/369QCyds9ISkqCmZkZatSogRo1aiA1NRULFy5Eeno6atasiaioKIONa8yYMVizZg0CAwPRp08fzJo1C/369YOzszPi4uJybYivmqOjI0qVKoXExETlOVlbW2PevHlISEjAuXPnsHLlSrRs2RKTJ0/Gnj17kJaWhjNnzhQ6u1atWoiIiDC47Ny5c6hVq1au27777ru4cOGC8vdbt27h8ePHAIA///wTM2fORI0aNfDFF18gODgYLVu2xP79+1/w2RNpFzuKHUWkZuwodhSRmrGj2FFEasaOYkeVBPyQvIRp2LAh+vbtC19fX9y9exfOzs5o27Ytxo0bh8jISFy8eBGTJ09GSkoKKlSogO7duyMxMRELFy7EjRs3EBISgn379sHV1RW2trY4d+4coqKiEB0djUmTJuHPP/802O3kdbC2tka/fv0wY8YMnDx5EjExMZgwYQJiY2Ph4OCAsmXL4vvvv0dgYCBu3bqFnTt3IiUlRTleVEE+/fRThIWFYcOGDbhx4wbWr1+PvXv3Ksefymno0KHw8/NDWFgYrl69iqlTpyq7CVWsWBF79+7Ft99+i7i4OISHh+PKlSuoX7/+S389iLSEHcWOIlIzdhQ7ikjN2FHsKCI1Y0exo7SOH5KXQF9//TVKly6NRYsWAcg6iYCDgwM+/fRTDBs2DLVq1VIO+l+hQgWsXLkSp0+fRrdu3bB69WosXrwY9erVw1dffYXKlStjwIABGDZsGCwtLTFo0CBcvnz5dT49AMCkSZPQqlUrjB49Gv3794eFhQVWrVqFUqVKoV69epg7dy7WrFmDDz/8ECtWrMCiRYvg7Oxc6NxGjRph4cKF8Pf3R7du3RAcHIylS5eiVatWuW7bs2dPjB49GrNnz8bgwYPh6uqKChUqAMj6FtXX1xdXrlxBjx49MGbMGPTt2xf9+vV76a8Fkdawo9hRRGrGjmJHEakZO4odRaRm7Ch2lJaZyeveV4GIiIiIiIiIiIiI6DXhL8mJiIiIiIiIiIiIyGTxQ3IiIiIiIiIiIiIiMln8kJyIiIiIiIiIiIiITBY/JCciIiIiIiIiIiIik8UPyYmIiIiIiIiIiIjIZPFDciIiIiIiIiIiIiIyWfyQnIiIiIiIiIiIiIhMFj8kJyIiIiIiIiIiIiKTxQ/JiYiIiIiIiIiIiMhk8UNyIiIiIiIiIiIiIjJZ/JCciIiIiIiIiIiIiEzW/wMlCjg9Q+7SVgAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], "source": [ "adata = load_replogle(\"essential\", \"normalized\")\n", "metadata = adata.obs\n", "embeddings_pca = embed_by_pca(adata)\n", "embeddings_aligned = align_by_centering(embeddings_pca, metadata)\n", "map_data = aggregate_by_mean(embeddings_aligned, metadata)\n", - "metrics = benchmark(map_data, recall_threshold_pairs=recall_threshold_pairs)\n", + "metrics = benchmark(map_data, recall_thr_pairs=recall_threshold_pairs)\n", "plot_recall(metrics)" ] }, @@ -91,7 +129,7 @@ "embeddings_scvi = embed_by_scvi(adata)\n", "embeddings_aligned = align_by_centering(embeddings_scvi, metadata)\n", "map_data = aggregate_by_mean(embeddings_aligned, metadata)\n", - "metrics = benchmark(map_data)\n", + "metrics = benchmark(map_data, recall_thr_pairs=recall_threshold_pairs)\n", "print(metrics)" ] }