examples/LMF/LMF.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Load Data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(400, 15)"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import os\n",
    "X = np.load(os.path.join(\"..\", \"..\", \"data\", \"dog.npz\"))[\"X\"]\n",
    "X.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sklearn.metrics import root_mean_squared_error\n",
    "from sklearn.metrics import f1_score\n",
    "from sklearn.metrics import precision_score\n",
    "from sklearn.metrics import recall_score\n",
    "from sklearn.metrics import auc\n",
    "from sklearn.metrics import accuracy_score\n",
    "from sklearn.metrics import matthews_corrcoef\n",
    "from sklearn.metrics import precision_recall_curve\n",
    "from sklearn.metrics import roc_curve\n",
    "\n",
    "def get_scores(y_true, y_pred, binary=True):\n",
    "\n",
    "    fpr, tpr, _ = roc_curve(y_true, y_pred, pos_label=1)\n",
    "    roc_auc= auc(fpr, tpr)\n",
    "\n",
    "    precision_line, recall_line, _ = precision_recall_curve(y_true, y_pred, pos_label=1)\n",
    "    pr_auc = auc(recall_line, precision_line)\n",
    "    rmse = root_mean_squared_error(y_true, y_pred)\n",
    "\n",
    "    if binary:\n",
    "        acc = accuracy_score(y_true, y_pred)\n",
    "        f1 = f1_score(y_true, y_pred, average=\"weighted\")\n",
    "        precision = precision_score(y_true, y_pred, average=\"weighted\")\n",
    "        recall = recall_score(y_true, y_pred, average=\"weighted\")\n",
    "        mcc = matthews_corrcoef(y_true, y_pred)\n",
    "    else:\n",
    "        acc, f1, precision, recall, mcc = None, None, None, None, None\n",
    "\n",
    "    return {\n",
    "            \"roc_auc\":roc_auc, \n",
    "            \"pr_auc\":pr_auc,\n",
    "            \"accuracy\":acc,\n",
    "            \"f1\":f1,\n",
    "            \"precision\":precision,\n",
    "            \"recall\":recall,\n",
    "            \"mcc\":mcc,\n",
    "            \"rmse\":rmse\n",
    "            }"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create Training Set and MASK"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys; sys.path.append(os.path.join(\"..\", \"..\", \"scripts\"))\n",
    "from sample_with_mask import sample_matrix\n",
    "\n",
    "X, Xtrain, MASK, removed_coords= sample_matrix(X, sample_ratio=0.2, random_state=42, stratify=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgMAAAGzCAYAAACy+RS/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABKXUlEQVR4nO3deVxUZfs/8M+wDYiICwIDCoiYuKJiEiqSSeKSgZWZmWtZmj2JPKGSKdriqJVaapLmQplL9UUjLYwMTQszQXJJUQRxYXEhAVEHZM7vj37O0yAMnOEMM+N83r7u1yvOcp1rjHGuue/73EcmCIIAIiIislhWxk6AiIiIjIvFABERkYVjMUBERGThWAwQERFZOBYDREREFo7FABERkYVjMUBERGThWAwQERFZOBYDREREFo7FAFEjW7BgAWQymV7nbtq0CTKZDOfPn5c2qX85f/48ZDIZNm3aZLBrEJFpYTFAJMLJkyfxwgsvwNPTE3K5HB4eHhg7dixOnjxp7NSIiPQm47MJiOonMTERY8aMQcuWLfHiiy+iXbt2OH/+PNavX4/r169j27ZtGDlyZJ1x7t69i7t378Le3l50DlVVVaisrIRcLte7d6Eu58+fR7t27bBx40ZMnDjRINcgItNiY+wEiMzBuXPnMG7cOPj6+uKXX35B69atNftmzJiBkJAQjBs3DseOHYOvr2+NMcrLy+Ho6AgbGxvY2Oj31rO2toa1tbVe5xIR1YbDBET18P777+PWrVtYu3atViEAAC4uLvj0009RXl6OpUuXAvjfvIC//voLzz//PFq0aIH+/ftr7fu327dv4/XXX4eLiwucnJzw5JNP4vLly5DJZFiwYIHmuJrmDPj4+OCJJ57AwYMH0adPH9jb28PX1xeff/651jWKi4vxxhtvoFu3bmjatCmaNWuGoUOH4s8//5Twb4qIzBF7Bojq4bvvvoOPjw9CQkJq3D9gwAD4+Phg9+7dWttHjRqFDh06YNGiRdA1Ijdx4kR89dVXGDduHB555BHs378fw4cPr3d+2dnZeOaZZ/Diiy9iwoQJ2LBhAyZOnIjAwEB06dIFAJCTk4OdO3di1KhRaNeuHYqKivDpp58iNDQUf/31Fzw8POp9PSJ6sLAYIKpDSUkJ8vPzERERofO47t27IykpCWVlZZptAQEB2LJli87zMjIy8NVXXyEqKgrLly8HALz66quYNGlSvb+1Z2Vl4ZdfftEUK88++yzatm2LjRs34oMPPgAAdOvWDWfOnIGV1f86BMeNGwd/f3+sX78e8+bNq9e1iOjBw2ECojrc+3B3cnLSedy9/aWlpZptU6dOrTN+cnIygH8KgH/7z3/+U+8cO3furNVr0bp1a3Ts2BE5OTmabXK5XFMIVFVV4fr162jatCk6duyIjIyMel+LiB48LAaI6nDvQ/7f3/hrUlPR0K5duzrj5+XlwcrK6r5j/fz86p2jl5fXfdtatGiBv//+W/OzWq3G8uXL0aFDB8jlcri4uKB169Y4duwYSkpK6n0tInrwsBggqoOzszMUCgWOHTum87hjx47B09MTzZo102xzcHAwdHoAUOsdBv+ep7Bo0SJER0djwIAB2Lx5M/bs2YOUlBR06dIFarW6UfIkItPEOQNE9fDEE09g3bp1OHjwoOaugH87cOAAzp8/j1deeUV0bG9vb6jVauTm5qJDhw6a7dnZ2Q3KubpvvvkGAwcOxPr167W237hxAy4uLpJei4jMC3sGiOohJiYGDg4OeOWVV3D9+nWtfcXFxZg6dSqaNGmCmJgY0bHDw8MBAJ988onW9pUrV+qfcA2sra3vu6Ph66+/xuXLlyW9DhGZH/YMENVDhw4dkJCQgLFjx6Jbt273rUB47do1bN26Fe3btxcdOzAwEE8//TRWrFiB69eva24tPHPmDABIttLgE088gbfffhuTJk1C3759cfz4cXz55Ze1LpJERJaDxQBRPY0aNQr+/v5QKpWaAqBVq1YYOHAg3nzzTXTt2lXv2J9//jnc3d2xdetW7NixA2FhYdi+fTs6duyo17LFNXnzzTdRXl6OLVu2YPv27ejVqxd2796NOXPmSBKfiMwXn01AZKIyMzPRs2dPbN68GWPHjjV2OkT0AOOcASITcPv27fu2rVixAlZWVhgwYIARMiIiS8JhAiITsHTpUqSnp2PgwIGwsbHBDz/8gB9++AEvv/wy2rZta+z0iOgBx2ECIhOQkpKChQsX4q+//sLNmzfh5eWFcePGYe7cuXo/4ZCIqL44TEBkAh5//HEcPHgQxcXFqKioQHZ2NuLi4lgIEFmYX375BSNGjICHhwdkMhl27txZ5zn79u1Dr169IJfL4efnh02bNom+LosBIiIiE1FeXo6AgACsXr26Xsfn5uZi+PDhGDhwIDIzMxEVFYWXXnoJe/bsEXVdDhMQERGZIJlMhh07diAyMrLWY2bPno3du3fjxIkTmm3PPfccbty4oXkIWn2wZ4CIiMiAVCoVSktLtZpKpZIkdlpaGsLCwrS2hYeHIy0tTVQckxmQtLHzNHYKRERkJu5WGHYZ7cprOXUfVE/KVZ9j4cKFWtvi4uKwYMGCBscuLCyEm5ub1jY3NzeUlpbi9u3b9X5YmskUA0RERCZDXSVZqNjYWERHR2ttk8vlksWXAosBIiIiA5LL5Qb78Hd3d0dRUZHWtqKiIjRr1kzUI9RZDBAREVUnqI2dQb0EBwfj+++/19qWkpKC4OBgUXE4gZCIiKg6tVq6JsLNmzeRmZmJzMxMAP/cOpiZmYkLFy4A+GfIYfz48Zrjp06dipycHMyaNQunT5/GJ598gq+++gozZ84UdV32DBAREVUjGKln4MiRIxg4cKDm53tzDSZMmIBNmzahoKBAUxgAQLt27bB7927MnDkTH330Edq0aYPPPvsM4eHhoq5rMusM8G4CIiKqL0PfTVCRf1KyWHYeXSSLZSiiewauXbuGDRs2IC0tDYWFhQD+mcDQt29fTJw4Ea1bt5Y8SSIiokYlsnvf3InqGfjjjz8QHh6OJk2aICwsTHNvY1FREfbu3Ytbt25hz5496N27t844KpXqvgUXWrTyh0wm0+MlEBGRpTF4z8DFPyWLZdc2QLJYhiKqGHjkkUcQEBCA+Pj4+z64BUHA1KlTcezYsTpXPlqwYMF9CzDIrJrCyrqZiNSJiMhSsRiQlqhiwMHBAUePHoW/v3+N+0+fPo2ePXvi9u3bOuOwZ4CIiBrC4MVAXoZksey8e0kWy1BEzRlwd3fH4cOHay0GDh8+fN+yiDWpaQEGFgJERGQyzGSdAamIKgbeeOMNvPzyy0hPT8egQYPumzOwbt06fPDBBwZJlIiIiAxDVDEwffp0uLi4YPny5fjkk09QVfXP2s3W1tYIDAzEpk2b8OyzzxokUSIiokbDuwnqp7KyEteuXQMAuLi4wNbWtkGJcJ0BIiKqL0PPGVCdOyRZLHn7RySLZSh6r0Boa2sLhUIhZS5ERERkBFyOmIiIqDoLGyZgMUBERFQd7yYgIiKycOoqY2fQqPgIYyIiIgvHngEiIqLqOExARERk4SxsAiGHCYiIiCwcewaIiIiq4zABERGRheMwAREREVkSyYuBixcvYvLkyTqPUalUKC0t1Wp6PiKBiIhIcoJQJVkzB5IXA8XFxUhISNB5jFKphLOzs1YT1GVSp0JERKQfQS1dMwOin1qYlJSkc39OTg7++9//ah5vXBOVSgWVSqW1rUUrf8hkMjGpEBGRhTL0UwvvZO6SLJZ9jycki2UooicQRkZGQiaT6ezWr+tDXS6XQy6XizqHiIio0XACoW4KhQKJiYlQq9U1toyMDEPkSURE1HgsbJhAdDEQGBiI9PT0WvfX1WtARERk8tRV0jUzIHqYICYmBuXl5bXu9/PzQ2pqaoOSIiIiosYjegKhodjYeRo7BSIiMhMGn0B4+GvJYtn3GSVZLEPhCoRERETVcQIhERERWRL2DBAREVVnJncBSIXFAN3ndv4BY6dABuLgEWLsFIjMA4cJiIiIyJKwZ4CIiKg6C+sZYDFARERUjbk8bVAqHCYgIiKycOwZICIiqo7DBERERBaOtxYSERFZOAvrGRA9Z+D27ds4ePAg/vrrr/v23blzB59//nmdMVQqFUpLS7WaiTwigYiIyOKIKgbOnDmDTp06YcCAAejWrRtCQ0NRUFCg2V9SUoJJkybVGUepVMLZ2VmrCeoy8dkTEREZgqCWrpkBUcXA7Nmz0bVrV1y5cgVZWVlwcnJCv379cOHCBVEXjY2NRUlJiVaTWTmJikFERGQwarV0zQyImjPw22+/4aeffoKLiwtcXFzw3Xff4dVXX0VISAhSU1Ph6OhYrzhyuRxyuVxrm0wmE5MKERERSURUz8Dt27dhY/O/+kEmk2HNmjUYMWIEQkNDcebMGckTJCIianQWNkwgqmfA398fR44cQadOnbS2r1q1CgDw5JNPSpcZERGRsZhJ975URPUMjBw5Elu3bq1x36pVqzBmzBjeFUBERGRmZIKJfHrb2HkaOwX6//gI4wcXH2FMD4q7FZcNGv/27hWSxXIYHiVZLEPhokNERETVmclYv1T4oCIiIiILx54BIiKi6ixsAiGLASIiouosbJiAxQAREVF1FtYzwDkDREREFo49A0RERNVxmICIiMjCcZiAiIiILAl7BoiIiKqzsJ4B0cXAqVOncOjQIQQHB8Pf3x+nT5/GRx99BJVKhRdeeAGPPfZYnTFUKhVUKpXWNkEQ+BhjIiIyDaaxUn+jETVMkJycjB49euCNN95Az549kZycjAEDBiA7Oxt5eXkYPHgwfv755zrjKJVKODs7azVBXab3iyAiIiL9iSoG3n77bcTExOD69evYuHEjnn/+eUyZMgUpKSnYu3cvYmJisHjx4jrjxMbGoqSkRKvJrJz0fhFERESSUqula2ZAVDFw8uRJTJw4EQDw7LPPoqysDM8884xm/9ixY3Hs2LE648jlcjRr1kyrcYiAiIhMBosB3e59aFtZWcHe3h7Ozs6afU5OTigpKZEuOyIiIjI4UcWAj48Pzp49q/k5LS0NXl5emp8vXLgAhUIhXXZERETGIKila2ZA1N0E06ZNQ1VVlebnrl27au3/4Ycf6nU3ARERkUkzk+59qYjqGZg6dSqGDx9e6/5Fixbhs88+a3BSRERERiUI0jWRVq9eDR8fH9jb2yMoKAiHDx/WefyKFSvQsWNHODg4oG3btpg5cybu3Lkj6ppcgZCIiMhEbN++HdHR0YiLi0NGRgYCAgIQHh6OK1eu1Hj8li1bMGfOHMTFxeHUqVNYv349tm/fjjfffFPUdVkMEBERVWekuwmWLVuGKVOmYNKkSejcuTPi4+PRpEkTbNiwocbjf/vtN/Tr1w/PP/88fHx8MHjwYIwZM6bO3oTqWAwQERFVJ2ExoFKpUFpaqtWqr8ILABUVFUhPT0dYWJhmm5WVFcLCwpCWllZjmn379kV6errmwz8nJwfff/89hg0bJurlshggIiIyoJpW3VUqlfcdd+3aNVRVVcHNzU1ru5ubGwoLC2uM/fzzz+Ptt99G//79YWtri/bt2+PRRx/lMAEREVGDSXhrYU2r7sbGxkqS5r59+7Bo0SJ88sknyMjIQGJiInbv3o133nlHVBw+tZCIiKgaQS3dg4rkcjnkcnmdx7m4uMDa2hpFRUVa24uKiuDu7l7jOfPmzcO4cePw0ksvAQC6deuG8vJyvPzyy5g7dy6srOr3nZ89A0RERCbAzs4OgYGB2Lt3r2abWq3G3r17ERwcXOM5t27duu8D39raGsA/TwOuL/YMEBERVWekRYeio6MxYcIE9O7dG3369MGKFStQXl6OSZMmAQDGjx8PT09PzZyDESNGYNmyZejZsyeCgoKQnZ2NefPmYcSIEZqioD5YDBAREVVnpGWER48ejatXr2L+/PkoLCxEjx49kJycrJlUeOHCBa2egLfeegsymQxvvfUWLl++jNatW2PEiBF47733RF1XJojpR6iFIAgNfuqgjZ1nQ9MgidzOP2DsFMhAHDxCjJ0CkSTuVlw2aPxba/4jWawm01ZKFstQJJkzIJfLcerUKSlCERERGZ9akK6ZAVHDBNHR0TVur6qqwuLFi9GqVSsA/6ygpItKpbpvwQUpeheIiIgkYWEPKhJVDKxYsQIBAQFo3ry51nZBEHDq1Ck4OjrW6wNdqVRi4cKFWttkVk0hs24mJh0iIiLDsLBiQNScgcWLF2Pt2rX47LPPtB5VbGtriz///BOdO3euV5yaegZatPJnz4CJ4JyBBxfnDNCDwuBzBj6aKlmsJjPiJYtlKKJ6BubMmYNBgwbhhRdewIgRI6BUKmFrayv6ojUtwMBCgIiITEbD59abFdETCB9++GGkp6fj6tWr6N27N06cOMEPciIierAY6amFxqLXOgNNmzZFQkICtm3bhrCwMFRVVUmdFxERETWSBi069Nxzz6F///5IT0+Ht7e3VDkREREZl5ncEiiVBq9A2KZNG7Rp00aKXIiIiEyDkVYgNBY+qIiIiMjC8dkERERE1XGYgIiIyLIJZnIXgFQ4TEBERGTh2DNARERUHYcJiIiILJyF3U3AYoCIiKg6C+sZ4JwBIiIiC8eeASIiouos7G4CFgNERETVWdgwQYOKgfLycnz11VfIzs6GQqHAmDFj0KpVqzrPU6lUUKlUWtsEQeDTD4mIiIxA1JyBzp07o7i4GABw8eJFdO3aFTNnzkRKSgri4uLQuXNn5Obm1hlHqVTC2dlZqwnqMv1eARERkdQEtXTNDIgqBk6fPo27d+8CAGJjY+Hh4YG8vDwcPnwYeXl56N69O+bOnVtnnNjYWJSUlGg1mZWTfq+AiIhIampBumYG9B4mSEtLQ3x8PJydnQEATZs2xcKFC/Hcc8/Vea5cLodcLtfaxiECIiIi4xBdDNz70L5z5w4UCoXWPk9PT1y9elWazIiIiIzE0p5NILoYGDRoEGxsbFBaWoqsrCx07dpVsy8vL69eEwiJiIhMmpl070tFVDEQFxen9XPTpk21fv7uu+8QEhLS8KyIiIio0cgEQTCJ8sfGztPYKdD/dzv/gLFTIANx8GCxTg+GuxWXDRr/ZsxIyWI1fX+HZLEMhYsOERERVWcmtwRKhcUAERFRdRY2Z4APKiIiIrJw7Bmg+5j6uLIpz2kw9b87IqofwcJ6BlgMEBERVWdhxQCHCYiIiCwcewaIiIiq4wqEREREFo7DBERERGRJ2DNARERUnYX1DLAYICIiqsZEVupvNKKGCTIyMpCbm6v5+YsvvkC/fv3Qtm1b9O/fH9u2batXHJVKhdLSUq1maX/xREREpkJUMTBp0iScO3cOAPDZZ5/hlVdeQe/evTF37lw8/PDDmDJlCjZs2FBnHKVSCWdnZ60mqMv0ewVERERSUwvSNTMg6qmFTZo0walTp+Dt7Y1evXph2rRpmDJlimb/li1b8N577+HkyZM646hUKqhUKq1tLVr5QyaTiUyfLBFXICQiQz+1sPTFxyWL1Wx9imSxDEXUnIEmTZrg2rVr8Pb2xuXLl9GnTx+t/UFBQVrDCLWRy+WQy+Va21gIEBGRqbC05YhFDRMMHToUa9asAQCEhobim2++0dr/1Vdfwc/PT7rsiIiIyOBE9QwsWbIE/fr1Q2hoKHr37o0PP/wQ+/btQ6dOnZCVlYVDhw5hx44dhsqViIiocbBnoHYeHh44evQogoODkZycDEEQcPjwYfz4449o06YNfv31VwwbNsxQuRIRETUOtYTNDIiaQGhINnaexk6BzAQnEBKRoScQlowbJFks5y/2ShbLULjoEBERUTWWNoGQxQAREVF1FlYM8EFFREREFo49A0RERNWZycQ/qbAYICIiqsbS5gxwmICIiMjCsWeAiIioOg4TEBERWTZLGyZgMUBERFSdhfUMcM4AERGRhRNVDPznP//BgQMNXwpWpVKhtLRUq5nIqshEREQQ1NI1cyCqGFi9ejUeffRRPPTQQ1iyZAkKCwv1uqhSqYSzs7NWE9RlesUiIiKSnIU9qEj0MMGPP/6IYcOG4YMPPoCXlxciIiKwa9cuqNX1f8WxsbEoKSnRajIrJ7GpEBERkQREFwPdunXDihUrkJ+fj82bN0OlUiEyMhJt27bF3LlzkZ2dXWcMuVyOZs2aaTWZTKbXCyAiIpIahwnqydbWFs8++yySk5ORk5ODKVOm4Msvv0THjh2lzI+IiKjxcZhAPC8vLyxYsAC5ublITk6WIiQRERE1ElHrDHh7e8Pa2rrW/TKZDI8//niDkyIiIjImc+nel4qonoHc3Fy0atXKULkQERGZBGPOGVi9ejV8fHxgb2+PoKAgHD58WOfxN27cwPTp06FQKCCXy/HQQw/h+++/F3VNrkBIRERUjbF6BrZv347o6GjEx8cjKCgIK1asQHh4OLKysuDq6nrf8RUVFXj88cfh6uqKb775Bp6ensjLy0Pz5s1FXZfFABERkYlYtmwZpkyZgkmTJgEA4uPjsXv3bmzYsAFz5sy57/gNGzaguLgYv/32G2xtbQEAPj4+oq/L5YiJiIiqE2SStZpW3VWpVPddsqKiAunp6QgLC9Nss7KyQlhYGNLS0mpMMykpCcHBwZg+fTrc3NzQtWtXLFq0CFVVVaJe7gPbM3A7v+HLJhOJxd87qi8HjxBjp0A6SDlMoFQqsXDhQq1tcXFxWLBggda2a9euoaqqCm5ublrb3dzccPr06Rpj5+Tk4Oeff8bYsWPx/fffIzs7G6+++ioqKysRFxdX7xwf2GKAiIjIFMTGxiI6Olprm1wulyS2Wq2Gq6sr1q5dC2trawQGBuLy5ct4//33WQwQERE1hKCWblVcuVxerw9/FxcXWFtbo6ioSGt7UVER3N3dazxHoVDA1tZW67b/Tp06obCwEBUVFbCzs6tXjpwzQEREVI0xbi20s7NDYGAg9u7dq9mmVquxd+9eBAcH13hOv379kJ2drfV8oDNnzkChUNS7EABYDBAREZmM6OhorFu3DgkJCTh16hSmTZuG8vJyzd0F48ePR2xsrOb4adOmobi4GDNmzMCZM2ewe/duLFq0CNOnTxd1XQ4TEBERVSMIxnl43ujRo3H16lXMnz8fhYWF6NGjB5KTkzWTCi9cuAArq/99j2/bti327NmDmTNnonv37vD09MSMGTMwe/ZsUdeVCYIgSPpK9GRj5ylpPM7qJiJTxrsJGuZuxWWDxr8U9Jhksdr8/rNksQyFwwREREQWTnQxsGrVKowfPx7btm0DAHzxxRfo3Lkz/P398eabb+Lu3bt1xqhpAQYT6aAgIiKCoJZJ1syBqDkD7777LpYuXYrBgwdj5syZyMvLw/vvv4+ZM2fCysoKy5cvh62t7X2LK1RX0wIMMqumkFk3E/8KiIiIJGZp309FzRnw8/PD0qVL8dRTT+HPP/9EYGAgEhISMHbsWADAjh07MGvWLJw9e1ZnHJVKdd9SjC1a+UMmk66C4pwBIjJlnDPQMIaeM5DXK6zug+rJO+MnyWIZiqiegfz8fPTu3RsAEBAQACsrK/To0UOzv1evXsjPz68zTk0LMEhZCBAREVH9iZoz4O7ujr/++gsAcPbsWVRVVWl+BoCTJ0/W+IhFIiIic8I5AzqMHTsW48ePR0REBPbu3YtZs2bhjTfewPXr1yGTyfDee+/hmWeeMVSuREREjcLS5gyIKgYWLlwIBwcHpKWlYcqUKZgzZw4CAgIwa9Ys3Lp1CyNGjMA777xjqFyJiIjIALjoEBGREXACYcMYegJhTrfBksXyPf6jZLEMhcsRExERVWOs5YiNhSsQEhERWbgHtmeAXXAPLlMeAuLvHdGDQcyjhx8ED2wxQEREpC81hwmIiIjIkrBngIiIqBpLm0DIYoCIiKgac1k5UCosBoiIiKoxjRV4Gg/nDBAREVk40T0DBQUFWLNmDQ4ePIiCggJYWVnB19cXkZGRmDhxIqytrQ2RJxERUaOxtGECUT0DR44cQadOnfD999+jsrISZ8+eRWBgIBwdHfHGG29gwIABKCsrqzOOSqVCaWmpVjORVZGJiIigFmSSNXMgqhiIiorCzJkzceTIERw4cACbNm3CmTNnsG3bNuTk5ODWrVt466236oyjVCrh7Oys1QR13UUEERERSU/Ug4qaNGmCEydOwNfXFwCgVqthb2+Pixcvws3NDSkpKZg4cSIuX9b9AAmVSgWVSqW1rUUrf8hk5lFBkXFxBUIiMvSDio63GyFZrG6530kWy1BEzRlwdXVFQUGBphgoKirC3bt30axZMwBAhw4dUFxcXGccuVwOuVyutY2FABERmQpLG7kWNUwQGRmJqVOnIjk5GampqRg7dixCQ0Ph4OAAAMjKyoKnp7SPIiYiIiLDEtUz8O6776KgoAAjRoxAVVUVgoODsXnzZs1+mUwGpVIpeZJERESNyVwm/klFVDHQtGlTbN++HXfu3MHdu3fRtGlTrf2DBw+WNDkiIiJj4HLE9WBvby91HkRERGQkXI6YiIioGkubQMhigIiIqBrOGSAiIrJwljZngA8qIiIisnDsGSAiIqqGwwREREQWzsLmD3KYgIiIyNLp1TNQUVGBnTt3Ii0tDYWFhQAAd3d39O3bFxEREbCzs5M0SSIiosZkacMEonsGsrOz0alTJ0yYMAFHjx6FWq2GWq3G0aNHMX78eHTp0gXZ2dmGyJWIiKhRCIJMsmYORPcMTJs2Dd26dcPRo0c1Tyu8p7S0FOPHj8f06dOxZ88eyZIkIiIiwxFdDPz66684fPjwfYUAADRr1gzvvPMOgoKCJEmOiIjIGNTGTqCRiS4GmjdvjvPnz6Nr16417j9//jyaN2+uM4ZKpYJKpdLaJggCZDLz6E4hIqIHmwDL+jwSPWfgpZdewvjx47F8+XIcO3YMRUVFKCoqwrFjx7B8+XJMnDgRL7/8ss4YSqUSzs7OWk1Ql+n9IoiIiEh/MkEQ/ziGJUuW4KOPPkJhYaHm27wgCHB3d0dUVBRmzZql8/yaegZatPJnzwDVy+38A8ZOoVYOHiHGToHIItytuGzQ+PvcRkkW69GiryWLZSh6FQP35Obmat1a2K5dO70TsbHz1PtcsiwsBojI0MXAz27PShbrsaKvJItlKA1adKhdu3YIDg5GcHCwphC4ePEiJk+eLElyRERExiBAJlkzB5KvQFhcXIyEhASpwxIREZGBiL6bICkpSef+nJwcvZMhIiIyBby1sA6RkZGQyWTQNdWAEwGJiMicmUv3vlREDxMoFAokJiZqliGu3jIyMgyRJxERERmI6GIgMDAQ6enpte6vq9eAiIjI1KklbOZA9DBBTEwMysvLa93v5+eH1NTUBiVFRERkTObyIS4V0cVASIju+6gdHR0RGhqqd0JERETUuEQXA0RERA86S5tAyGKAiIioGrVl1QLSLzpERERE5kXyYqCoqAhvv/221GGJiIgajRoyyZo5kLwYKCwsxMKFC6UOS0RE1GgECZs5ED1n4NixYzr3Z2Vl6Z0MERGRKeCthXXo0aNHrQsL3dte13LEKpUKKpVKa1t9ziMiIiLpiS4GWrZsiaVLl2LQoEE17j958iRGjBihM4ZSqbxvKEFm1RQy62Zi0yEiIpKc2sK+nIouBgIDA5Gfnw9vb+8a99+4caPO5YhjY2MRHR2tta1FK3+xqRARERmEuYz1S0V0MTB16lSdyxF7eXlh48aNOmPI5XLI5XKtbRwiICIiMg7RxcDIkSN17m/RogUmTJigd0JERETGZmkTCCW/tfDixYuYPHmy1GGJiIgajVomXTMHkhcDxcXFSEhIkDosERERGYjoYYKkpCSd+3NycvROhoiIyBQYc+XA1atX4/3330dhYSECAgKwcuVK9OnTp87ztm3bhjFjxiAiIgI7d+4UdU3RxUBkZGSt6wzcw8mARERkzox1N8H27dsRHR2N+Ph4BAUFYcWKFQgPD0dWVhZcXV1rPe/8+fN44403EBISotd1RQ8TKBQKJCYmQq1W19gyMjL0SoSIiOhBpFKpUFpaqtWqL7x3z7JlyzBlyhRMmjQJnTt3Rnx8PJo0aYINGzbUGr+qqgpjx47FwoUL4evrq1eOoouBwMBApKen17q/rl4DIiIiUyflBEKlUglnZ2etplQq77tmRUUF0tPTERYWptlmZWWFsLAwpKWl1Zrr22+/DVdXV7z44ot6v17RwwQxMTE61xnw8/NDamqq3gkREREZm5S3Fta00F71tXYA4Nq1a6iqqoKbm5vWdjc3N5w+fbrG2AcPHsT69euRmZnZoBxFFwN1jUc4OjoiNDRU74SIiIiMTcr+7ZoW2pNCWVkZxo0bh3Xr1sHFxaVBsUQXA0RERCQ9FxcXWFtbo6ioSGt7UVER3N3d7zv+3LlzOH/+vNbzgNTqf/o0bGxskJWVhfbt29fr2pKvM0BERGTujLHokJ2dHQIDA7F3797/5aFWY+/evQgODr7veH9/fxw/fhyZmZma9uSTT2LgwIHIzMxE27Zt631tvXsGLl26hObNm6Np06Za2ysrK5GWloYBAwboG5qIiMiojLUccXR0NCZMmIDevXujT58+WLFiBcrLyzFp0iQAwPjx4+Hp6QmlUgl7e3t07dpV6/zmzZsDwH3b6yK6GCgoKEBERATS09Mhk8nw/PPP45NPPtEUBcXFxRg4cCCqqqrEhiYiIrJoo0ePxtWrVzF//nwUFhaiR48eSE5O1kwqvHDhAqyspO/Ulwki7wOcMGECsrKysGrVKty4cQNz5syBTCbDjz/+iBYtWqCoqAgKhUIzblFfNnaeoo4ny3U7/4CxU6iVg4d+C34QkTh3Ky4bNP6nbV6QLNYrlzZLFstQRJcXP/30Ez7++GP07t0bYWFh+PXXX6FQKPDYY4+huLgYAFcgJCIi8ybIpGvmQHQxUFJSghYtWmh+lsvlSExMhI+PDwYOHIgrV67UGaOm1Zi4UBEREZFxiC4GfH19cezYMa1tNjY2+Prrr+Hr64snnniizhg1rcYkqMvEpkJERGQQagmbORBdDAwdOhRr1669b/u9gqBHjx51fsuPjY1FSUmJVpNZOYlNhYiIyCAsrRgQfTfBe++9h1u3btUczMYG//d//4fLl3VP7KhpNSbOMyAiIjIO0T0DNjY2aNasWa37CwoKsHDhwgYlRUREZEyChM0cSH6zYnFxMRISEqQOS0RE1GiMsQKhMYkeJkhKStK5PycnR+9kiIiITIG5jPVLRXQxEBkZCZlMpnOSIMf/iYiIzIfoYQKFQoHExESo1eoaW0ZGhiHyJCIiajSWdjeB6GIgMDAQ6enpte6vq9eAiIjI1FnaBELRwwQxMTEoLy+vdb+fnx9SU1MblBSRuZL6uQl81gERNQbRxUBIiO5/nBwdHREaGqp3QkRERMZmLncBSEV0MUBERPSgM5exfqlI/1BkIiIiMivsGSAiIqrGXCb+SUWvYuD69es4duwYAgIC0LJlS1y7dg3r16+HSqXCqFGj0KlTJ6nzJCIiajRqCysHRBcDhw8fxuDBg1FaWormzZsjJSUFo0aNgo2NDdRqNRYvXoyDBw+iV69ehsiXiIiIJCZ6zsDcuXMxatQolJSU4M0330RkZCQGDRqEM2fOIDs7G8899xzeeecdQ+RKRETUKLjoUB3S09MRHR0NJycnzJgxA/n5+ZgyZYpm/2uvvYY//vhDZwyVSoXS0lKtxoWKiIjIVFjaokOii4GKigo4ODgAAGxtbdGkSRO4uLho9ru4uOD69es6YyiVSjg7O2s1QV0mNhUiIiKDYM9AHdq2bav1ZMJt27ZBoVBofi4oKNAqDmoSGxuLkpISrSazchKbChEREUlA9ATC5557DleuXNH8PHz4cK39SUlJ6NOnj84YcrkccrlcaxufdEhERKaCKxDWIS4uTuf+uXPnwtraWu+EiIiIjM3Sbi2UfAXC69evY9q0aVKHJSIiIgORvBgoLi5GQkKC1GGJiIgajaXdTSB6mCApKUnn/n9PLiQiIjJH5nIXgFREFwORkZGQyWQ61wXgZEAiIiLzIXqYQKFQIDExEWq1usaWkZFhiDyJiIgajRqCZM0ciO4ZCAwMRHp6OiIiImrcX1evAVme2/kHjJ1CrRw8QiSNJ/VrlTqe1K+X6EFlaZ9ioouBmJgYlJeX17rfz88PqampDUqKiIiIGo/oYiAkRPc3C0dHR4SGhuqdEBERkbFxAiEREZGFM5exfqmwGCAiIqrGskoBCRcd8vX1xdmzZ6UKR0RERI1EdM/Axx9/XOP2CxcuYOPGjXB3dwcAvP766w3LjIiIyEg4Z6AOUVFR8PT0hI2N9qlqtRqff/45bG1tIZPJWAwQEZHZEixsoEB0MfDyyy/j999/x5YtW9CpUyfNdltbW/z444/o3LmzpAkSERGRYYmeMxAfH4/58+cjPDwcq1at0uuiKpUKpaWlWo0LFRERkalQS9jMgV4TCEeOHIm0tDTs2LEDQ4cORWFhoajzlUolnJ2dtZqgLtMnFSIiIslZ2nLEet9N4OnpiZ9++gkDBgxAz549RX2zj42NRUlJiVaTWTnpmwoRERE1QIPWGZDJZIiNjcXgwYNx8OBBKBSKep0nl8shl8vvi0VERGQKzOP7vHQkWWcgMDAQM2bMQIsWLXDx4kVMnjxZirBERERGwWGCBiouLkZCQoLUYYmIiMhARA8TJCUl6dyfk5OjdzJERESmwFzuApCK6GIgMjISMplM54RBjv8TEZE5s7RFh0QPEygUCiQmJkKtVtfYMjIyDJEnERFRo+E6A3UIDAxEenp6rfvr6jUgIiIi0yJ6mCAmJgbl5eW17vfz80NqamqDkiLjup1/wNgp6OTgEWLsFGoldW5S/7+QOp4p/78gaghLGyYQXQyEhOh+8zs6OiI0NFTvhIiIiIzNXLr3pSL5rYVERERkXhq0AiEREdGDSG1hc99YDBAREVVjWaWABMWAIAjYt28fsrOzoVAoEB4eDltbWylyIyIiokYguhgYNmwYtm7dCmdnZxQXF2PYsGE4fPgwXFxccP36dTz00EP45Zdf0Lp1a0PkS0REZHDm8kwBqYieQJicnAyVSgUAeOutt1BWVoZz587hypUryMvLg6OjI+bPny95okRERI1FkPCPOWjQMMHPP/+MpUuXol27dgCANm3aYMmSJZgyZYrO81QqlaaguEcQBC5jTEREZAR63Vp470P777//Rvv27bX2+fn5IT8/X+f5SqUSzs7OWk1Ql+mTChERkeS4HHE9TJw4EU899RQqKyuRm5urta+wsBDNmzfXeX5sbCxKSkq0mszKSZ9UiIiIJKeGIFkzB6KHCSZMmKD574iICNy6dUtr///93/+hR48eOmPI5XLI5XKtbRwiICIiU2EuY/1SEV0MbNy4Uef+uLg4WFtb650QERERNS7JlyMuLi7Gq6++KnVYIiKiRsM5Aw1UXFyMhIQEqcMSERE1GkEQJGvmQPQwQVJSks79OTk5eidDRERk6VavXo33338fhYWFCAgIwMqVK9GnT58aj123bh0+//xznDhxAgAQGBiIRYsW1Xp8bUQXA5GRkZDJZDqrHU4GJCIic2asuwC2b9+O6OhoxMfHIygoCCtWrEB4eDiysrLg6up63/H79u3DmDFj0LdvX9jb22PJkiUYPHgwTp48CU9Pz3pfV/QwgUKhQGJiItRqdY0tIyNDbEgiIiKTIuWcAZVKhdLSUq1WfeG9e5YtW4YpU6Zg0qRJ6Ny5M+Lj49GkSRNs2LChxuO//PJLvPrqq+jRowf8/f3x2WefQa1WY+/evaJer+iegcDAQKSnpyMiIqLG/XX1GpDlcfAIMXYKZkvqv7vb+QckjUdEdVMqlVi4cKHWtri4OCxYsEBrW0VFBdLT0xEbG6vZZmVlhbCwMKSlpdXrWrdu3UJlZSVatmwpKkfRxUBMTAzKy8tr3e/n54fU1FSxYYmIiEyGlOsMxMbGIjo6Wmtb9bV2AODatWuoqqqCm5ub1nY3NzecPn26XteaPXs2PDw8EBYWJipH0cVASIjubyqOjo4IDQ0VG5aIiMhkSDlnoKaF9gxh8eLF2LZtG/bt2wd7e3tR5zboQUVEREQkDRcXF1hbW6OoqEhre1FREdzd3XWe+8EHH2Dx4sX46aef0L17d9HXlnydASIiInNnjHUG7OzsEBgYqDX5795kwODg4FrPW7p0Kd555x0kJyejd+/eer1e0cXApUuXcO3aNc3PBw4cwNixYxESEoIXXnih3pMciIiITJWxViCMjo7GunXrkJCQgFOnTmHatGkoLy/HpEmTAADjx4/XmmC4ZMkSzJs3Dxs2bICPjw8KCwtRWFiImzdvirqu6GLg6aefxqFDhwAA3377LR599FHcvHkT/fr1w61btxAaGopdu3aJDUtERGQyBAn/iDF69Gh88MEHmD9/Pnr06IHMzEwkJydrJhVeuHABBQUFmuPXrFmDiooKPPPMM1AoFJr2wQcfiLquTBB5H2DTpk1x/PhxtGvXDo888ghGjhyJ2bNna/avWrUKGzZsEL3egI1d/RdHIMOS+vYz3lpoOvj/lh4UdysuGzT+4LZDJIv148VkyWIZiuieARsbG5SVlQEAcnNzMXToUK39Q4cORVZWljTZERERGYEagmTNHIguBkJDQ7F161YAQM+ePbFv3z6t/ampqXUugVjTakxcqIiIiEwFH1RUh8WLFyMkJAT5+fno378/5s6diz/++AOdOnVCVlYWtm/fjvj4eJ0xalqNSWbVFDLrZmLTISIiogYSPWcAAM6dO4e33noLu3fv1sxYtLGxwcMPP4yYmBhERkbqPF+lUt23LnOLVv58wJGJ4Ljyg4v/b+lBYeg5AwPbPC5ZrNRLKZLFMhS9Fh1q3749tm7dCkEQcOXKFajVari4uMDW1rZe59e0GhMLASIiMhVSLkdsDhq06JBMJoObmxsUCoWmELh48SImT54sSXJERERkeJKvQFhcXIyEhASpwxIRETUatSBI1syB6GGCpKQknftzcnL0ToaIiMgUmMdHuHREFwORkZGQyWQ6b5fg+D8REZH5ED1MoFAokJiYCLVaXWMTu/IgERGRqeGiQ3UIDAxEenp6rfvr6jUgIiIydZZWDIgeJoiJiUF5eXmt+/38/JCamtqgpEgcqe8dJyKydJb2pVZ0MRASonuREUdHR4SGhuqdEBERETUuvRYdIiIiepCZS/e+VFgMEBERVcMVCImIiMiiiC4GPvzwQ+Tl5RkiFyIiIpNgaY8wFl0MxMTEoH379nj88cexfft2VFRUGCIvIiIio7G0Wwv1Gib47LPP4OjoiHHjxsHDwwNRUVE4ceKE1LkRERFRI9CrGBg2bBh27tyJS5cuYdasWdizZw8CAgLQp08frFu3DmVlZTrPV6lUKC0t1Wrm0pVCREQPPg4TiODq6opZs2bh1KlT2LdvHzp37oyZM2dCoVDoPE+pVMLZ2VmrCWrdBQQREVFj4TBBHWp7CFFISAg2bdqE/Px8LF++XGeM2NhYlJSUaDWZlZPYVIiIiEgCotcZqKvLo1mzZpgyZYrOY+RyOeRyudY2PumQiIhMhaWtMyC6GFCr1YbIg4iIyGSozWSsXyqSLzp08eJFTJ48WeqwREREjUaQ8I85kLwYKC4uRkJCgtRhiYiIyEBEDxMkJSXp3J+Tk6N3MkRERKbA0oYJRBcDkZGRkMlkOicScjIgERGZM3Pp3peK6GEChUKBxMREqNXqGltGRoYh8iQiIiIDEd0zEBgYiPT0dERERNS4v65eAwJu5x8wdgo6OXiEGDsFMhNS/y7zd49MBYcJ6hATE4Py8vJa9/v5+SE1NbVBSRERERmTpQ0TiC4GQkJ0V+6Ojo4IDQ3VOyEiIiJqXKKLASIiogcdhwmIiIgsnKUNE+i16NCuXbswf/58/PrrrwCAn3/+GcOGDcOQIUOwdu1aSRMkIiIiwxJdDHz66acYOXIkvv/+ewwbNgybN29GZGQkPD094ePjg6ioKHz00UeGyJWIiKhRCIJasmYORA8TfPzxx/jkk08wZcoUpKamYtiwYfjwww/x6quvAgAeeeQRLF26FDNmzJA8WSIiosag5jCBbrm5uQgPDwcADBw4EFVVVRgwYIBm/6OPPoq8vDzpMiQiImpkgiBI1syB6GKgVatWmg/7/Px83L17FxcuXNDsz8vLQ8uWLXXGUKlUKC0t1Wrm8hdGRET0oBE9TBAREYEXX3wREyZMQFJSEsaPH4///ve/sLKygkwmQ0xMDAYPHqwzhlKpxMKFC7W2yayaQmbdTGw6REREkrO0YQKZIPIreXl5OWbOnIm0tDT07dsXK1euxMcff4y5c+eisrISoaGh2L59O1xdXWuNoVKpoFKptLa1aOVvMQ844nLEZCz83aMHxd2KywaN79mii2SxLv99UrJYhiK6GKjNnTt3UFlZCScnJ73Ot7HzlCINs8B/kMlY+LtHDwoWA9LSa52Bmtjb28PJyQkXL17E5MmTpQpLRETU6NSCIFkzB5IVA/cUFxcjISFB6rBERESNRpDwjzkQPYEwKSlJ5/6cnBy9kyEiIqLGJ7oYiIyMhEwm03kroKVMBCQiogeTpd3uLnqYQKFQIDExEWq1usaWkZFhiDyJiIgajRqCZM0ciC4GAgMDkZ6eXuv+unoNiIiIyLSIHiaIiYlBeXl5rfv9/PyQmpraoKQedLx9ioyFv3tE9WNpX2olW2egoSxpnQEiImoYQ68z0NKpg2SxisvOShbLUET3DBARET3oTOR7cqORfJ0BIiIiMi/sGSAiIqrGXO4CkIpexcDt27exdetWHDx4EAUFBbCysoKvry8iIyMxaNAgqXMkIiJqVJY2TCB6AmF2djbCwsJw+/ZtyOVyXLp0CcOGDcO1a9dw5MgRPPXUU9iyZQtsbMTVGZxASERE9WXoCYTNHH0li1Vabvor84qeM/D6669jyJAhKCwsxIULF6BUKqFWq3Ho0CGcOnUKf/zxB959911D5EpERNQoLO1BRaJ7BhwdHZGZmYkOHf657aKiogJNmzZFQUEBWrVqhW+//RZRUVHIzc2tNYZKpYJKpdLa1qKVP5cxJiKiejF0z4BjEx/JYpXfOi9ZLEMR3TPQvHlzlJWVaX6+desW7t69Czs7OwBA9+7dUVBQoDOGUqmEs7OzVhPUZTrPISIiIsMQXQw8/vjjiI6OxunTp5Gbm4upU6eiR48ecHJyAgBcuHABrq6uOmPExsaipKREq8msnPR7BURERBKztGEC0XcTLF26FBEREejcuTNkMhnatm2LHTt2aPZfvXoVMTExOmPI5XLI5XKtbRwiICIiU8G7Cerp7NmzUKlU8Pf3F33nQE14NwEREdWXoecM2Nt7SRbrzp0LksUyFL1XIOzQoQO6du16XyFw8eJFTJ48ucGJERERGYsg4R9zIPlyxMXFxUhISJA6LBERUaMRBEGyJtbq1avh4+MDe3t7BAUF4fDhwzqP//rrr+Hv7w97e3t069YN33//vehriu7fT0pK0rk/J8f0F1cgIiLSxVhzBrZv347o6GjEx8cjKCgIK1asQHh4OLKysmqcnP/bb79hzJgxUCqVeOKJJ7BlyxZERkYiIyMDXbt2rfd1Rc8ZsLKygkwm0/kXJZPJUFVVJSYs5wwQEVG9GXrOgK2En0mVInINCgrCww8/jFWrVgEA1Go12rZti//85z+YM2fOfcePHj0a5eXl2LVrl2bbI488gh49eiA+Pr7e1xU9TKBQKJCYmAi1Wl1jy8jIEBuSiIjIpAgSNpVKhdLSUq1WfeE94J9F/NLT0xEWFqbZZmVlhbCwMKSlpdWYZ1pamtbxABAeHl7r8bW/YJFGjBghzJs3r9b9mZmZgkwmExu2Xu7cuSPExcUJd+7cMalYjPdgxzPl3BjPdGIxnunFMxVxcXH31QhxcXH3HXf58mUBgPDbb79pbY+JiRH69OlTY2xbW1thy5YtWttWr14tuLq6ispRdDHwyy+/CD/88EOt+2/evCns27dPbNh6KSkpEQAIJSUlJhWL8R7seKacG+OZTizGM714puLOnTtCSUmJVqup4DFmMSB6AmFISIjO/Y6OjggNDRUbloiI6IFU00J7NXFxcYG1tTWKioq0thcVFcHd3b3Gc9zd3UUdXxvJby0kIiIi8ezs7BAYGIi9e/dqtqnVauzduxfBwcE1nhMcHKx1PACkpKTUenxtGr50IBEREUkiOjoaEyZMQO/evdGnTx+sWLEC5eXlmDRpEgBg/Pjx8PT0hFKpBADMmDEDoaGh+PDDDzF8+HBs27YNR44cwdq1a0Vd16yKAblcjri4uHp1tzRmLMZ7sOOZcm6MZzqxGM/04pmj0aNH4+rVq5g/fz4KCwvRo0cPJCcnw83NDcA/DwO0svpfp37fvn2xZcsWvPXWW3jzzTfRoUMH7Ny5U9QaA0ADnk1AREREDwbOGSAiIrJwLAaIiIgsHIsBIiIiC8digIiIyMKxGCAiIrJwZlMMiH2+sy6//PILRowYAQ8PD8hkMuzcuVPvWEqlEg8//DCcnJzg6uqKyMhIZGVl6R1vzZo16N69O5o1a4ZmzZohODgYP/zwg97x/m3x4sWQyWSIiorSO8aCBQsgk8m0mr+/v97xLl++jBdeeAGtWrWCg4MDunXrhiNHjugVy8fH577cZDIZpk+frle8qqoqzJs3D+3atYODgwPat2+Pd955p0GPNi0rK0NUVBS8vb3h4OCAvn374o8//qjXuXX93gqCgPnz50OhUMDBwQFhYWE4e/as3vESExMxePBgtGrVCjKZDJmZmXrFqqysxOzZs9GtWzc4OjrCw8MD48ePR35+vt65LViwAP7+/nB0dESLFi0QFhaG33//Xe94/zZ16lTIZDKsWLFC73gTJ0687/dwyJAhDcrv1KlTePLJJ+Hs7AxHR0c8/PDDuHDhgl7xanqfyGQyvP/++6Jj3bx5E6+99hratGkDBwcHdO7cWefT8uqKV1RUhIkTJ8LDwwNNmjTBkCFDdP4ekzTMohi493znuLg4ZGRkICAgAOHh4bhy5Ype8crLyxEQEIDVq1c3OLf9+/dj+vTpOHToEFJSUlBZWYnBgwejvLxcr3ht2rTB4sWLkZ6ejiNHjuCxxx5DREQETp482aA8//jjD3z66afo3r17g+IAQJcuXVBQUKBpBw8e1CvO33//jX79+sHW1hY//PAD/vrrL3z44Ydo0aKFXvH++OMPrbxSUlIAAKNGjdIr3pIlS7BmzRqsWrUKp06dwpIlS7B06VKsXLlSr3gA8NJLLyElJQVffPEFjh8/jsGDByMsLAyXL9f9iNO6fm+XLl2Kjz/+GPHx8fj999/h6OiI8PBw3LlzR6945eXl6N+/P5YsWdKg3G7duoWMjAzMmzcPGRkZSExMRFZWFp588km9X+tDDz2EVatW4fjx4zh48CB8fHwwePBgXL16Va949+zYsQOHDh2Ch4eHzuPqE2/IkCFav49bt27VO965c+fQv39/+Pv7Y9++fTh27BjmzZsHe3t7veL9O6+CggJs2LABMpkMTz/9tOhY0dHRSE5OxubNm3Hq1ClERUXhtddeQ1JSkujcBEFAZGQkcnJy8O233+Lo0aPw9vZGWFiY3v+mUj2JepKBkfTp00eYPn265ueqqirBw8NDUCqVDY4NQNixY0eD49xz5coVAYCwf/9+yWK2aNFC+Oyzz/Q+v6ysTOjQoYOQkpIihIaGCjNmzNA7VlxcnBAQEKD3+f82e/ZsoX///pLEqsmMGTOE9u3bC2q1Wq/zhw8fLkyePFlr21NPPSWMHTtWr3i3bt0SrK2thV27dmlt79WrlzB37lxRsar/3qrVasHd3V14//33Ndtu3LghyOVyYevWraLj/Vtubq4AQDh69KheudXk8OHDAgAhLy9Pknj3HnDz008/6R3v0qVLgqenp3DixAnB29tbWL58eZ2xaos3YcIEISIiol7n1yfe6NGjhRdeeEGyeNVFREQIjz32mF6xunTpIrz99tta2+r7O109XlZWlgBAOHHihGZbVVWV0Lp1a2HdunV1xiP9mXzPgD7PdzamkpISAEDLli0bHKuqqgrbtm1DeXm56HWm/2369OkYPnz4fc+81tfZs2fh4eEBX19fjB07ttauyrokJSWhd+/eGDVqFFxdXdGzZ0+sW7dOkhwrKiqwefNmTJ48GTKZTK8Yffv2xd69e3HmzBkAwJ9//omDBw9i6NChesW7e/cuqqqq7vs25+DgoHfvyj25ubkoLCzU+n/s7OyMoKAgk32fyGQyNG/evMGxKioqsHbtWjg7OyMgIECvGGq1GuPGjUNMTAy6dOnS4JwAYN++fXB1dUXHjh0xbdo0XL9+Xe/cdu/ejYceegjh4eFwdXVFUFBQg4Y3/62oqAi7d+/Giy++qNf5ffv2RVJSEi5fvgxBEJCamoozZ85g8ODBomOpVCoA0HqPWFlZQS6XN/g9QrqZfDFw7do1VFVVaZZivMfNzQ2FhYVGyqpmarUaUVFR6Nevn+ilIP/t+PHjaNq0KeRyOaZOnYodO3agc+fOesXatm0bMjIyNOtYN1RQUBA2bdqE5ORkrFmzBrm5uQgJCUFZWZnoWDk5OVizZg06dOiAPXv2YNq0aXj99deRkJDQ4Dx37tyJGzduYOLEiXrHmDNnDp577jn4+/vD1tYWPXv2RFRUFMaOHatXPCcnJwQHB+Odd95Bfn4+qqqqsHnzZqSlpaGgoEDvPAFo3gvm8D65c+cOZs+ejTFjxqBZs2Z6x9m1axeaNm0Ke3t7LF++HCkpKXBxcdEr1pIlS2BjY4PXX39d73z+bciQIfj888+xd+9eLFmyBPv378fQoUNRVVUlOtaVK1dw8+ZNLF68GEOGDMGPP/6IkSNH4qmnnsL+/fsbnGtCQgKcnJzw1FNP6XX+ypUr0blzZ7Rp0wZ2dnYYMmQIVq9ejQEDBoiO5e/vDy8vL8TGxuLvv/9GRUUFlixZgkuXLjX4PUK6mdWzCUzd9OnTceLEiQZXsB07dkRmZiZKSkrwzTffYMKECdi/f7/oguDixYuYMWMGUlJSah1bFOvf34q7d++OoKAgeHt746uvvhL9zUKtVqN3795YtGgRAKBnz544ceIE4uPjMWHChAbluX79egwdOrTOsV9dvvrqK3z55ZfYsmULunTpgszMTERFRcHDw0Pv/L744gtMnjwZnp6esLa2Rq9evTBmzBikp6frnac5qaysxLPPPgtBELBmzZoGxRo4cCAyMzNx7do1rFu3Ds8++yx+//13uLq6ioqTnp6Ojz76CBkZGXr3IlX33HPPaf67W7du6N69O9q3b499+/Zh0KBBomKp1WoAQEREBGbOnAkA6NGjB3777TfEx8c3+JHxGzZswNixY/X+N2LlypU4dOgQkpKS4O3tjV9++QXTp0+Hh4eH6N5IW1tbJCYm4sUXX0TLli1hbW2NsLAwDB06tEETd6luJt8zoM/znY3htddew65du5Camoo2bdo0KJadnR38/PwQGBgIpVKJgIAAfPTRR6LjpKen48qVK+jVqxdsbGxgY2OD/fv34+OPP4aNjY1e31Kqa968OR566CFkZ2eLPlehUNxX4HTq1EnvYYd78vLy8NNPP+Gll15qUJyYmBhN70C3bt0wbtw4zJw5s0G9LO3bt8f+/ftx8+ZNXLx4EYcPH0ZlZSV8fX0blOu994Ipv0/uFQJ5eXlISUlpUK8AADg6OsLPzw+PPPII1q9fDxsbG6xfv150nAMHDuDKlSvw8vLSvE/y8vLw3//+Fz4+Pg3K8R5fX1+4uLjo9T5xcXGBjY2NQd4rBw4cQFZWlt7vldu3b+PNN9/EsmXLMGLECHTv3h2vvfYaRo8ejQ8++ECvmIGBgcjMzMSNGzdQUFCA5ORkXL9+vcHvEdLN5IsBfZ7v3JgEQcBrr72GHTt24Oeff0a7du0kv4ZardaMpYkxaNAgHD9+HJmZmZrWu3dvjB07FpmZmbC2tm5wbjdv3sS5c+egUChEn9uvX7/7bsM8c+YMvL29G5TTxo0b4erqiuHDhzcozq1bt7SeDgYA1tbWmm9qDeHo6AiFQoG///4be/bsQURERIPitWvXDu7u7lrvk9LSUvz+++8m8T65VwicPXsWP/30E1q1aiX5NfR9n4wbNw7Hjh3Tep94eHggJiYGe/bskSS3S5cu4fr163q9T+zs7PDwww8b5L2yfv16BAYG6j3XorKyEpWVlQZ5nzg7O6N169Y4e/Ysjhw50uD3COlmFsMEdT3fWaybN29qVei5ubnIzMxEy5Yt4eXlJSrW9OnTsWXLFnz77bdwcnLSjM86OzvDwcFBdG6xsbEYOnQovLy8UFZWhi1btmDfvn16/aPk5OR039wFR0dHtGrVSu85DW+88QZGjBgBb29v5OfnIy4uDtbW1hgzZozoWDNnzkTfvn2xaNEiPPvsszh8+DDWrl0r+jnc/6ZWq7Fx40ZMmDABNjYN+/UeMWIE3nvvPXh5eaFLly44evQoli1bhsmTJ+sdc8+ePRAEAR07dkR2djZiYmLg7+9fr9/lun5vo6Ki8O6776JDhw5o164d5s2bBw8PD0RGRuoVr7i4GBcuXNCsB3Dvw8jd3f2+3gZdsRQKBZ555hlkZGRg165dqKqq0rxPWrZsCTs7O1G5tWrVCu+99x6efPJJKBQKXLt2DatXr8bly5drvY20rtdavTixtbWFu7s7OnbsKDpey5YtsXDhQjz99NNwd3fHuXPnMGvWLPj5+SE8PFyv/GJiYjB69GgMGDAAAwcORHJyMr777jvs27dPr3jAP8Xi119/jQ8//LDGGPWNFRoaipiYGDg4OMDb2xv79+/H559/jmXLlukV7+uvv0br1q3h5eWF48ePY8aMGYiMjNRrQiKJYNR7GURYuXKl4OXlJdjZ2Ql9+vQRDh06pHes1NRUAcB9bcKECaJj1RQHgLBx40a9cps8ebLg7e0t2NnZCa1btxYGDRok/Pjjj3rFqklDby0cPXq0oFAoBDs7O8HT01MYPXq0kJ2drXe87777Tujatasgl8sFf39/Ye3atXrHEgRB2LNnjwBAyMrKalAcQRCE0tJSYcaMGYKXl5dgb28v+Pr6CnPnzhVUKpXeMbdv3y74+voKdnZ2gru7uzB9+nThxo0b9Tq3rt9btVotzJs3T3BzcxPkcrkwaNAgnX8PdcXbuHFjjfvj4uJExbp3a2JNLTU1VXRut2/fFkaOHCl4eHgIdnZ2gkKhEJ588knh8OHDer/W6uq6tVBXvFu3bgmDBw8WWrduLdja2gre3t7ClClThMLCwgblt379esHPz0+wt7cXAgIChJ07dzYo3qeffio4ODjU+ftXV6yCggJh4sSJgoeHh2Bvby907NhR+PDDD2u9pbeueB999JHQpk0bwdbWVvDy8hLeeuutBr3nqH5kgsBZGURERJbM5OcMEBERkWGxGCAiIrJwLAaIiIgsHIsBIiIiC8digIiIyMKxGCAiIrJwLAaIiIgsHIsBIiIiC8digIiIyMKxGCAiIrJwLAaIiIgs3P8DYFRNXIjuUOMAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgMAAAGzCAYAAACy+RS/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABJb0lEQVR4nO3deVxU9f4/8NdhGxARFQQGVFyw0FRUXEJDMklcrkGrmiVq2dX0pnJDIxdcStwyS03S1KiuS92LRmqYGlgWpoJkehUlCBcWFxIQdEDm/P64P+fbsAycmTPMjPN6+vg8HvE5c97zHp2J93w+n/M5giiKIoiIiMhq2Zg6ASIiIjItFgNERERWjsUAERGRlWMxQEREZOVYDBAREVk5FgNERERWjsUAERGRlWMxQEREZOVYDBAREVk5FgNEJjZx4kR06NBBr3MXLVoEQRDkTYiIrA6LAaJ6CILQqJaammrqVImIDCLw3gREdfviiy+0fv7ss89w8OBBfP7551r9Tz75JDw9PfV+nqqqKqjVaigUCsnn3rt3D/fu3YOjo6Pez09ExGKAqJFmzJiBDRs2oKGPTEVFBZo1a9ZEWRERGY7TBEQGePzxx9G9e3ekp6dj8ODBaNasGd5++20AwNdff41Ro0bB29sbCoUCnTt3xtKlS1FdXa0Vo+aagT/++AOCIGD16tXYtGkTOnfuDIVCgX79+uHEiRNa59a1ZkAQBMyYMQN79uxB9+7doVAo8MgjjyA5OblW/qmpqejbty8cHR3RuXNnfPzxx1yHQGSF7EydAJGlu3nzJkaMGIGxY8fipZde0kwZfPrpp2jevDmioqLQvHlzfP/991i4cCFKS0uxatWqBuNu374dZWVl+Pvf/w5BELBy5Uo888wzyMnJgb29vc5zjx49isTERLz++utwcXHBhx9+iGeffRaXLl2Cm5sbAODUqVMYPnw4lEolFi9ejOrqaixZsgRt2rQx/C+FiCyLSESNMn36dLHmRyYkJEQEIMbHx9d6fEVFRa2+v//972KzZs3Eu3fvavoiIyNFX19fzc+5ubkiANHNzU0sLi7W9H/99dciAPGbb77R9MXGxtbKCYDo4OAgZmdna/p+/fVXEYC4bt06Td/o0aPFZs2aiVevXtX0Xbx4UbSzs6sVk4gebJwmIDKQQqHApEmTavU7OTlp/rusrAw3btxAcHAwKioqcP78+QbjjhkzBq1atdL8HBwcDADIyclp8NzQ0FB07txZ83PPnj3RokULzbnV1dU4dOgQIiIi4O3trXmcn58fRowY0WB8InqwcJqAyEA+Pj5wcHCo1X/27FnMnz8f33//PUpLS7WOlZSUNBi3ffv2Wj/fLwz+/PNPyefeP//+udeuXcOdO3fg5+dX63F19RHRg43FAJGB/joCcN+tW7cQEhKCFi1aYMmSJejcuTMcHR2RkZGBuXPnQq1WNxjX1ta2zn6xERcAGXIuEVkfFgNERpCamoqbN28iMTERgwcP1vTn5uaaMKv/4+HhAUdHR2RnZ9c6VlcfET3YuGaAyAjufzP/6zfxyspKfPTRR6ZKSYutrS1CQ0OxZ88e5Ofna/qzs7Px7bffmjAzIjIFjgwQGcHAgQPRqlUrREZG4o033oAgCPj888/Naph+0aJF+O677zBo0CBMmzYN1dXVWL9+Pbp3747MzExTp0dETYgjA0RG4Obmhr1790KpVGL+/PlYvXo1nnzySaxcudLUqWkEBgbi22+/RatWrbBgwQJs2bIFS5YswdChQ7m9MZGV4XbERKQlIiICZ8+excWLF02dChE1EY4MEFmxO3fuaP188eJF7N+/H48//rhpEiIik+DIAJEVUyqVmDhxIjp16oS8vDxs3LgRKpUKp06dQpcuXUydHhE1ES4gJLJiw4cPx44dO1BYWAiFQoGgoCAsW7aMhQCRleE0AZEV27ZtG/744w/cvXsXJSUlSE5ORp8+fUydFpHV+uGHHzB69Gh4e3tDEATs2bOnwXNSU1PRp08fKBQK+Pn54dNPP5X8vCwGiIiIzER5eTkCAgKwYcOGRj0+NzcXo0aNwpAhQ5CZmYlZs2bh1VdfxYEDByQ9L9cMEBERmSFBELB7925ERETU+5i5c+di3759OHPmjKZv7NixuHXrFpKTkxv9XBwZICIiMiKVSoXS0lKtplKpZImdlpaG0NBQrb6wsDCkpaVJimM2CwjtHHxMnQIREVmIe5VXjRq/6kbDtwpvrLj1n2Hx4sVafbGxsVi0aJHBsQsLC+Hp6anV5+npidLSUty5c6fOG6nVxWyKASIiIrOhrpYtVExMDKKiorT6FAqFbPHlwGKAiIjIiBQKhdF++Xt5eaGoqEirr6ioCC1atGj0qADAYoCIiKg2UW3qDBolKCgI+/fv1+o7ePAggoKCJMXhAkIiIqKa1Gr5mgS3b99GZmam5s6hubm5yMzMxKVLlwD8b8phwoQJmsdPnToVOTk5mDNnDs6fP4+PPvoIX375JWbPni3peTkyQEREVINoopGBkydPYsiQIZqf7681iIyMxKeffoqCggJNYQAAHTt2xL59+zB79mx88MEHaNu2LT755BOEhYVJel6z2WeAVxMQEVFjGftqgsr8s7LFcvB+RLZYxiJ5ZODGjRvYunUr0tLSUFhYCOB/CxgGDhyIiRMnok2bNrInSURE1KQkDu9bOkkjAydOnEBYWBiaNWuG0NBQzbWNRUVFOHz4MCoqKnDgwAH07dtXZxyVSlVrw4VWbv4QBEGPl0BERNbG6CMDl3+VLZZDuwDZYhmLpGLg0UcfRUBAAOLj42v94hZFEVOnTsXp06cb3Plo0aJFtTZgEGyaw8a2hYTUiYjIWrEYkJekYsDJyQmnTp2Cv79/ncfPnz+P3r17486dOzrjcGSAiIgMYfRiIC9DtlgOvuZ/J1BJawa8vLxw/PjxeouB48eP19oWsS51bcDAQoCIiMyGhewzIBdJxcCbb76J1157Denp6Rg6dGitNQObN2/G6tWrjZIoERERGYekYmD69Olwd3fH+++/j48++gjV1f/bu9nW1haBgYH49NNP8cILLxglUSIioibDqwkap6qqCjdu3AAAuLu7w97e3qBEuM8AERE1lrHXDKh+PyZbLEXnR2WLZSx670Bob28PpVIpZy5ERERkAtyOmIiIqCYrmyZgMUBERFQTryYgIiKycupqU2fQpHgLYyIiIivHkQEiIqKaOE1ARERk5axsASGnCYiIiKwcRwaIiIhq4jQBERGRleM0AREREVkT2YuBy5cvY/LkyTofo1KpUFpaqtX0vEUCERGR7ESxWrZmCWQvBoqLi5GQkKDzMXFxcXB1ddVqorpM7lSIiIj0I6rlaxZA8l0Lk5KSdB7PycnBP//5T83tjeuiUqmgUqm0+lq5+UMQBCmpEBGRlTL2XQvvZu6VLZZjr7/JFstYJC8gjIiIgCAIOof1G/qlrlAooFAoJJ1DRETUZLiAUDelUonExESo1eo6W0ZGhjHyJCIiajpWNk0guRgIDAxEenp6vccbGjUgIiIye+pq+ZoFkDxNEB0djfLy8nqP+/n5ISUlxaCkiIiIqOlIXkBoLHYOPqZOgYiILITRFxAe/0q2WI79n5ctlrFwB0IiIqKauICQiIiIrAlHBoiIiGqykKsA5MJigGq5k/+jqVOg/8/JO9jUKRBZJ04TEBERkTXhyAAREVFNVjYywGKAiIioBku526BcOE1ARERk5TgyQEREVBOnCYiIiKwcLy0kIiKyclY2MiB5zcCdO3dw9OhR/Pe//6117O7du/jss88ajKFSqVBaWqrVzOQWCURERFZHUjFw4cIFdO3aFYMHD0aPHj0QEhKCgoICzfGSkhJMmjSpwThxcXFwdXXVaqK6THr2RERExiCq5WsWQFIxMHfuXHTv3h3Xrl1DVlYWXFxcMGjQIFy6dEnSk8bExKCkpESrCTYukmIQEREZjVotX7MAktYM/Pzzzzh06BDc3d3h7u6Ob775Bq+//jqCg4ORkpICZ2fnRsVRKBRQKBRafYIgSEmFiIiIZCJpZODOnTuws/u/+kEQBGzcuBGjR49GSEgILly4IHuCRERETc7KpgkkjQz4+/vj5MmT6Nq1q1b/+vXrAQBPPfWUfJkRERGZioUM78tF0sjA008/jR07dtR5bP369Rg3bhyvCiAiIrIwgmgmv73tHHxMnQL9f7yFsfngLYyJ6nav8qpR49/Zt1a2WE6jZskWy1i46RAREVFNFjLXLxfeqIiIiMjKcWSAyIzJPWXDaQeiRrKyBYQsBoiIiGqysmkCFgNEREQ1WdnIANcMEBERWTmODBAREdXEaQIiIiIrx2kCIiIisiYcGSAiIqrJykYGJBcD586dw7FjxxAUFAR/f3+cP38eH3zwAVQqFV566SU88cQTDcZQqVRQqVRafaIo8jbGRERkHsxjp/4mI2maIDk5Gb169cKbb76J3r17Izk5GYMHD0Z2djby8vIwbNgwfP/99w3GiYuLg6urq1YT1WV6vwgiIiLSn6RiYMmSJYiOjsbNmzexbds2vPjii5gyZQoOHjyIw4cPIzo6GsuXL28wTkxMDEpKSrSaYOOi94sgIiKSlVotX7MAkoqBs2fPYuLEiQCAF154AWVlZXjuuec0x8ePH4/Tp083GEehUKBFixZajVMERERkNlgM6Hb/l7aNjQ0cHR3h6uqqOebi4oKSkhL5siMiIiKjk1QMdOjQARcvXtT8nJaWhvbt22t+vnTpEpRKpXzZERERmYKolq9ZAElXE0ybNg3V1dWan7t37651/Ntvv23U1QRERERmzUKG9+UiaWRg6tSpGDVqVL3Hly1bhk8++cTgpIiIiExKFOVrEm3YsAEdOnSAo6MjBgwYgOPHj+t8/Nq1a/Hwww/DyckJ7dq1w+zZs3H37l1Jz8kdCImIiMzErl27EBUVhdjYWGRkZCAgIABhYWG4du1anY/fvn073nrrLcTGxuLcuXPYsmULdu3ahbffflvS87IYICIiqslEVxOsWbMGU6ZMwaRJk9CtWzfEx8ejWbNm2Lp1a52P//nnnzFo0CC8+OKL6NChA4YNG4Zx48Y1OJpQE4sBIiKimmQsBlQqFUpLS7VazV14AaCyshLp6ekIDQ3V9NnY2CA0NBRpaWl1pjlw4ECkp6drfvnn5ORg//79GDlypKSXy2KAiIjIiOradTcuLq7W427cuIHq6mp4enpq9Xt6eqKwsLDO2C+++CKWLFmCxx57DPb29ujcuTMef/xxThMQEREZTMZLC+vadTcmJkaWNFNTU7Fs2TJ89NFHyMjIQGJiIvbt24elS5dKisO7FhIREdUgquW7UZFCoYBCoWjwce7u7rC1tUVRUZFWf1FREby8vOo8Z8GCBXj55Zfx6quvAgB69OiB8vJyvPbaa5g3bx5sbBr3nZ8jA0RERGbAwcEBgYGBOHz4sKZPrVbj8OHDCAoKqvOcioqKWr/wbW1tAfzvbsCNxZEBIiKimky06VBUVBQiIyPRt29f9O/fH2vXrkV5eTkmTZoEAJgwYQJ8fHw0aw5Gjx6NNWvWoHfv3hgwYACys7OxYMECjB49WlMUNAaLASIioppMtI3wmDFjcP36dSxcuBCFhYXo1asXkpOTNYsKL126pDUSMH/+fAiCgPnz5+Pq1ato06YNRo8ejXfffVfS8wqilHGEeoiiaPBdB+0cfAxNg2RyJ/9HU6dARuLkHWzqFIhkca/yqlHjV2z8h2yxmk1bJ1ssY5FlzYBCocC5c+fkCEVERGR6alG+ZgEkTRNERUXV2V9dXY3ly5fDzc0NwP92UNJFpVLV2nBBjtEFIiIiWVjZjYokFQNr165FQEAAWrZsqdUviiLOnTsHZ2fnRv1Cj4uLw+LFi7X6BJvmEGxbSEmHiIjIOKysGJC0ZmD58uXYtGkTPvnkE61bFdvb2+PXX39Ft27dGhWnrpGBVm7+HBkwE1wz8ODimgF6UBh9zcAHU2WL1WxmvGyxjEXSyMBbb72FoUOH4qWXXsLo0aMRFxcHe3t7yU9a1wYMLASIiMhsGL623qJIXkDYr18/pKen4/r16+jbty/OnDnDX+RERPRgMdFdC01Fr30GmjdvjoSEBOzcuROhoaGorq6WOy8iIiJqIgZtOjR27Fg89thjSE9Ph6+vr1w5ERERmZaFXBIoF4N3IGzbti3atm0rRy5ERETmwUQ7EJoKb1RERERk5XhvAiIiopo4TUBk3uS8Vp57KhBRXUQLuQpALpwmICIisnIcGSAiIqqJ0wRERERWzsquJmAxQEREVJOVjQxwzQAREZGV48gAERFRTVZ2NQGLASIiopqsbJrAoGKgvLwcX375JbKzs6FUKjFu3Di4ubk1eJ5KpYJKpdLqE0WRdz8kIiIyAUlrBrp164bi4mIAwOXLl9G9e3fMnj0bBw8eRGxsLLp164bc3NwG48TFxcHV1VWrieoy/V4BERGR3ES1fM0CSCoGzp8/j3v37gEAYmJi4O3tjby8PBw/fhx5eXno2bMn5s2b12CcmJgYlJSUaDXBxkW/V0BERCQ3tShfswB6TxOkpaUhPj4erq6uAIDmzZtj8eLFGDt2bIPnKhQKKBQKrT5OERAREZmG5GLg/i/tu3fvQqlUah3z8fHB9evX5cmMiIjIRKzt3gSSi4GhQ4fCzs4OpaWlyMrKQvfu3TXH8vLyGrWAkIiIyKxZyPC+XCQVA7GxsVo/N2/eXOvnb775BsHB8t1RjoiIiIzPoGKgplWrVhmUDBERkVngyAAREZGVs5BLAuXCYoCIiKgmKxsZ4I2KiIiIrBxHBqgWJ2/zXgR6J/9HU6dQL3P/uyOixhGtbGSAxQAREVFNVlYMcJqAiIjIynFkgIiIqCbuQEhERGTlOE1ARERE1oQjA0RERDVZ2cgAiwEiIqIaRNG6igFJ0wQZGRnIzc3V/Pz5559j0KBBaNeuHR577DHs3LmzUXFUKhVKS0u1mrX9xRMREZkLScXApEmT8PvvvwMAPvnkE/z9739H3759MW/ePPTr1w9TpkzB1q1bG4wTFxcHV1dXrSaqy/R7BURERHJTi/I1CyCIEr6SN2vWDOfOnYOvry/69OmDadOmYcqUKZrj27dvx7vvvouzZ8/qjKNSqaBSqbT6Wrn5QxAEiemTNeIOhER0r/KqUeOXvvKkbLFabDkoWyxjkbRmoFmzZrhx4wZ8fX1x9epV9O/fX+v4gAEDtKYR6qNQKKBQKLT6WAgQEZG5sLbtiCVNE4wYMQIbN24EAISEhODf//631vEvv/wSfn5+8mVHRERERidpZGDFihUYNGgQQkJC0LdvX7z33ntITU1F165dkZWVhWPHjmH37t3GypWIiKhpcGSgft7e3jh16hSCgoKQnJwMURRx/PhxfPfdd2jbti1++uknjBw50li5EhERNQ21jM0CSFpAaEx2Dj6mToEsBBcQEpGxFxCWvDxUtliunx+WLZaxcNMhIiKiGqxtASGLASIiopqsrBjgjYqIiIisHEcGiIiIarKQhX9yYTFARERUg7WtGeA0ARERkZXjyAAREVFNnCYgIiKybtY2TcBigIiIqCYrGxngmgEiIiIrJ6kY+Mc//oEffzR8K1iVSoXS0lKtZia7IhMREUFUy9csgaRiYMOGDXj88cfx0EMPYcWKFSgsLNTrSePi4uDq6qrVRHWZXrGIiIhkZ2U3KpI8TfDdd99h5MiRWL16Ndq3b4/w8HDs3bsXanXjX3FMTAxKSkq0mmDjIjUVIiIikoHkYqBHjx5Yu3Yt8vPz8cUXX0ClUiEiIgLt2rXDvHnzkJ2d3WAMhUKBFi1aaDVBEPR6AURERHLjNEEj2dvb44UXXkBycjJycnIwZcoU/Otf/8LDDz8sZ35ERERNj9ME0rVv3x6LFi1Cbm4ukpOT5QhJRERETUTSPgO+vr6wtbWt97ggCHjyyScNToqIiMiULGV4Xy6SRgZyc3Ph5uZmrFyIiIjMginXDGzYsAEdOnSAo6MjBgwYgOPHj+t8/K1btzB9+nQolUooFAo89NBD2L9/v6Tn5A6ERERENZhqZGDXrl2IiopCfHw8BgwYgLVr1yIsLAxZWVnw8PCo9fjKyko8+eST8PDwwL///W/4+PggLy8PLVu2lPS8LAaIiIjMxJo1azBlyhRMmjQJABAfH499+/Zh69ateOutt2o9fuvWrSguLsbPP/8Me3t7AECHDh0kPy+3IyYiIqpJFGRrde26q1Kpaj1lZWUl0tPTERoaqumzsbFBaGgo0tLS6kwzKSkJQUFBmD59Ojw9PdG9e3csW7YM1dXVkl7uAzsycCff8G2TLYWTd7CpU9DJmv4tzP218r2iP3P/uyN5yTlNEBcXh8WLF2v1xcbGYtGiRVp9N27cQHV1NTw9PbX6PT09cf78+Tpj5+Tk4Pvvv8f48eOxf/9+ZGdn4/XXX0dVVRViY2MbneMDWwwQERGZg5iYGERFRWn1KRQKWWKr1Wp4eHhg06ZNsLW1RWBgIK5evYpVq1axGCAiIjKEqJZvV1yFQtGoX/7u7u6wtbVFUVGRVn9RURG8vLzqPEepVMLe3l7rsv+uXbuisLAQlZWVcHBwaFSOXDNARERUgykuLXRwcEBgYCAOHz6s6VOr1Th8+DCCgoLqPGfQoEHIzs7Wuj/QhQsXoFQqG10IACwGiIiIzEZUVBQ2b96MhIQEnDt3DtOmTUN5ebnm6oIJEyYgJiZG8/hp06ahuLgYM2fOxIULF7Bv3z4sW7YM06dPl/S8nCYgIiKqQRRNc/O8MWPG4Pr161i4cCEKCwvRq1cvJCcnaxYVXrp0CTY2//c9vl27djhw4ABmz56Nnj17wsfHBzNnzsTcuXMlPa8giqIo6yvRk52Dj6zxzHlVstzMfZWzNf1bmDu+V/Rn7n931uZe5VWjxr8y4AnZYrX95XvZYhkLpwmIiIisnORiYP369ZgwYQJ27twJAPj888/RrVs3+Pv74+2338a9e/cajFHXBgxmMkBBREQEUS3I1iyBpDUD77zzDlauXIlhw4Zh9uzZyMvLw6pVqzB79mzY2Njg/fffh729fa3NFWqqawMGwaY5BNsW0l8BERGRzKzt+6mkNQN+fn5YuXIlnnnmGfz6668IDAxEQkICxo8fDwDYvXs35syZg4sXL+qMo1Kpam3F2MrNH4IgXwVlznOPcjP3uUxr+rcwd3yv6M/c/+6sjbHXDOT1CW34QY3km3FItljGImlkID8/H3379gUABAQEwMbGBr169dIc79OnD/Lz8xuMU9cGDHIWAkRERNR4ktYMeHl54b///S8A4OLFi6iurtb8DABnz56t8xaLREREloRrBnQYP348JkyYgPDwcBw+fBhz5szBm2++iZs3b0IQBLz77rt47rnnjJUrERFRk7C2NQOSioHFixfDyckJaWlpmDJlCt566y0EBARgzpw5qKiowOjRo7F06VJj5UpERERGIKkYsLGxwdtvv63VN3bsWIwdO1bWpIiIiEzJUob35cLtiImIiGow1XbEpsIdCImIiKzcAzsywGuCzYfc/xa8Fv3Bxb8/MhdSbj38IHhgiwEiIiJ9qTlNQERERNaEIwNEREQ1WNsCQhYDRERENfDSQiIiIitnbTsQcs0AERGRlZM8MlBQUICNGzfi6NGjKCgogI2NDTp16oSIiAhMnDgRtra2xsiTiIioyVjbNIGkkYGTJ0+ia9eu2L9/P6qqqnDx4kUEBgbC2dkZb775JgYPHoyysrIG46hUKpSWlmo10drGZIiIyGypRUG2ZgkkFQOzZs3C7NmzcfLkSfz444/49NNPceHCBezcuRM5OTmoqKjA/PnzG4wTFxcHV1dXrSaqGy4iiIiISH6CKOErebNmzXDmzBl06tQJAKBWq+Ho6IjLly/D09MTBw8exMSJE3H16lWdcVQqFVQqlVZfKzd/CIJlVFBkWtyBkIjuVer+PWOo3zqOli1Wj9xvZItlLJLWDHh4eKCgoEBTDBQVFeHevXto0aIFAKBLly4oLi5uMI5CoYBCodDqYyFARETmwtpmriVNE0RERGDq1KlITk5GSkoKxo8fj5CQEDg5OQEAsrKy4OPjY5REiYiIyDgkjQy88847KCgowOjRo1FdXY2goCB88cUXmuOCICAuLk72JImIiJqSpSz8k4ukYqB58+bYtWsX7t69i3v37qF58+Zax4cNGyZrckRERKbA7YgbwdHRUe48iIiIyES4HTEREVEN1raAkMUAERFRDVwzQEREZOWsbc0Ab1RERERk5TgyQEREVAOnCYiIiKycla0f5DQBERGRtdNrZKCyshJ79uxBWloaCgsLAQBeXl4YOHAgwsPD4eDgIGuSRERETcnapgkkjwxkZ2eja9euiIyMxKlTp6BWq6FWq3Hq1ClMmDABjzzyCLKzs42RKxERUZMQRUG2ZgkkjwxMmzYNPXr0wKlTpzR3K7yvtLQUEyZMwPTp03HgwAHZkiQiIiLjkVwM/PTTTzh+/HitQgAAWrRogaVLl2LAgAGyJEdERGQKalMn0MQkFwMtW7bEH3/8ge7du9d5/I8//kDLli11xlCpVFCpVFp9oihCECxjOIWIiB5sIqzr95HkNQOvvvoqJkyYgPfffx+nT59GUVERioqKcPr0abz//vuYOHEiXnvtNZ0x4uLi4OrqqtVEdZneL4KIiIj0J4ii9NsxrFixAh988AEKCws13+ZFUYSXlxdmzZqFOXPm6Dy/rpGBVm7+HBmgRrmT/6OpU6iXk3ewqVMgsgr3Kq8aNX6q5/OyxXq86CvZYhmLXsXAfbm5uVqXFnbs2FHvROwcfPQ+l6wLiwEiMnYx8L3nC7LFeqLoS9liGYtBmw517NgRQUFBCAoK0hQCly9fxuTJk2VJjoiIyBRECLI1SyD7DoTFxcVISEiQOywREREZieSrCZKSknQez8nJ0TsZIiIic8BLCxsQEREBQRCga6kBFwISEZEls5ThfblIniZQKpVITEzUbENcs2VkZBgjTyIiIjISycVAYGAg0tPT6z3e0KgBERGRuVPL2CyB5GmC6OholJeX13vcz88PKSkpBiVFRERkSpbyS1wukouB4GDd11E7OzsjJCRE74SIiIioaUkuBoiIiB501raAkMUAERFRDWrrqgXk33SIiIiILIvsxUBRURGWLFkid1giIqImo4YgW7MEshcDhYWFWLx4sdxhiYiImowoY7MEktcMnD59WufxrKwsvZMhIiIyB7y0sAG9evWqd2Oh+/0NbUesUqmgUqm0+hpzHhEREclPcjHQunVrrFy5EkOHDq3z+NmzZzF69GidMeLi4mpNJQg2zSHYtpCaDhERkezUVvblVHIxEBgYiPz8fPj6+tZ5/NatWw1uRxwTE4OoqCitvlZu/lJTISIiMgpLmeuXi+RiYOrUqTq3I27fvj22bdumM4ZCoYBCodDq4xQBERGRaUguBp5++mmdx1u1aoXIyEi9EyIiIjI1a1tAKPulhZcvX8bkyZPlDktERNRk1IJ8zRLIXgwUFxcjISFB7rBERERkJJKnCZKSknQez8nJ0TsZIiIic2DKnQM3bNiAVatWobCwEAEBAVi3bh369+/f4Hk7d+7EuHHjEB4ejj179kh6TsnFQERERL37DNzHxYBERGTJTHU1wa5duxAVFYX4+HgMGDAAa9euRVhYGLKysuDh4VHveX/88QfefPNNBAcH6/W8kqcJlEolEhMToVar62wZGRl6JUJERPQgUqlUKC0t1Wo1N967b82aNZgyZQomTZqEbt26IT4+Hs2aNcPWrVvrjV9dXY3x48dj8eLF6NSpk145Si4GAgMDkZ6eXu/xhkYNiIiIzJ2cCwjj4uLg6uqq1eLi4mo9Z2VlJdLT0xEaGqrps7GxQWhoKNLS0urNdcmSJfDw8MArr7yi9+uVPE0QHR2tc58BPz8/pKSk6J0QERGRqcl5aWFdG+3V3GsHAG7cuIHq6mp4enpq9Xt6euL8+fN1xj569Ci2bNmCzMxMg3KUXAw0NB/h7OyMkJAQvRMiIiIyNTnHt+vaaE8OZWVlePnll7F582a4u7sbFEtyMUBERETyc3d3h62tLYqKirT6i4qK4OXlVevxv//+O/744w+t+wGp1f8b07Czs0NWVhY6d+7cqOeWfZ8BIiIiS2eKTYccHBwQGBiIw4cP/18eajUOHz6MoKCgWo/39/fHb7/9hszMTE176qmnMGTIEGRmZqJdu3aNfm69RwauXLmCli1bonnz5lr9VVVVSEtLw+DBg/UNTUREZFKm2o44KioKkZGR6Nu3L/r374+1a9eivLwckyZNAgBMmDABPj4+iIuLg6OjI7p37651fsuWLQGgVn9DJBcDBQUFCA8PR3p6OgRBwIsvvoiPPvpIUxQUFxdjyJAhqK6ulhqaiIjIqo0ZMwbXr1/HwoULUVhYiF69eiE5OVmzqPDSpUuwsZF/UF8QJV4HGBkZiaysLKxfvx63bt3CW2+9BUEQ8N1336FVq1YoKiqCUqnUzFs0lp2Dj6THk/W6k/+jqVOol5O3fht+EJE09yqvGjX+x21fki3W3698IVssY5FcXhw6dAgffvgh+vbti9DQUPz0009QKpV44oknUFxcDIA7EBIRkWUTBfmaJZBcDJSUlKBVq1aanxUKBRITE9GhQwcMGTIE165dazBGXbsxcaMiIiIi05BcDHTq1AmnT5/W6rOzs8NXX32FTp064W9/+1uDMerajUlUl0lNhYiIyCjUMjZLILkYGDFiBDZt2lSr/35B0KtXrwa/5cfExKCkpESrCTYuUlMhIiIyCmsrBiRfTfDuu++ioqKi7mB2dvjPf/6Dq1d1L+yoazcmrjMgIiIyDckjA3Z2dmjRokW9xwsKCrB48WKDkiIiIjIlUcZmCWS/WLG4uBgJCQlyhyUiImoyptiB0JQkTxMkJSXpPJ6Tk6N3MkRERObAUub65SK5GIiIiIAgCDoXCXL+n4iIyHJIniZQKpVITEyEWq2us2VkZBgjTyIioiZjbVcTSC4GAgMDkZ6eXu/xhkYNiIiIzJ21LSCUPE0QHR2N8vLyeo/7+fkhJSXFoKSIiIio6UguBoKDdd+IxdnZGSEhIXonREREZGqWchWAXCQXA0RERA86S5nrl4v8N0UmIiIii8KRASIiohosZeGfXPQqBm7evInTp08jICAArVu3xo0bN7BlyxaoVCo8//zz6Nq1q9x5EhERNRm1lZUDkouB48ePY9iwYSgtLUXLli1x8OBBPP/887Czs4Narcby5ctx9OhR9OnTxxj5EhERkcwkrxmYN28enn/+eZSUlODtt99GREQEhg4digsXLiA7Oxtjx47F0qVLjZErERFRk+CmQw1IT09HVFQUXFxcMHPmTOTn52PKlCma4zNmzMCJEyd0xlCpVCgtLdVq3KiIiIjMhbVtOiS5GKisrISTkxMAwN7eHs2aNYO7u7vmuLu7O27evKkzRlxcHFxdXbWaqC6TmgoREZFRcGSgAe3atdO6M+HOnTuhVCo1PxcUFGgVB3WJiYlBSUmJVhNsXKSmQkRERDKQvIBw7NixuHbtmubnUaNGaR1PSkpC//79dcZQKBRQKBRafbzTIRERmQvuQNiA2NhYncfnzZsHW1tbvRMiIiIyNWu7tFD2HQhv3ryJadOmyR2WiIiIjET2YqC4uBgJCQlyhyUiImoy1nY1geRpgqSkJJ3H/7q4kIiIyBJZylUAcpFcDEREREAQBJ37AnAxIBERkeWQPE2gVCqRmJgItVpdZ8vIyDBGnkRERE1GDVG2ZgkkjwwEBgYiPT0d4eHhdR5vaNSArM+d/B9NnUK9nLyDZY1nzq8VkP/1Ej2orO23mORiIDo6GuXl5fUe9/PzQ0pKikFJERERUdORXAwEB+v+ZuHs7IyQkBC9EyIiIjI1LiAkIiKycpYy1y8XFgNEREQ1WFcpIOOmQ506dcLFixflCkdERERNRPLIwIcfflhn/6VLl7Bt2zZ4eXkBAN544w3DMiMiIjIRrhlowKxZs+Dj4wM7O+1T1Wo1PvvsM9jb20MQBBYDRERksUQrmyiQXAy89tpr+OWXX7B9+3Z07dpV029vb4/vvvsO3bp1kzVBIiIiMi7Jawbi4+OxcOFChIWFYf369Xo9qUqlQmlpqVbjRkVERGQu1DI2S6DXAsKnn34aaWlp2L17N0aMGIHCwkJJ58fFxcHV1VWrieoyfVIhIiKSnbVtR6z31QQ+Pj44dOgQBg8ejN69e0v6Zh8TE4OSkhKtJti46JsKERERGcCgfQYEQUBMTAyGDRuGo0ePQqlUNuo8hUIBhUJRKxYREZE5sIzv8/KRZZ+BwMBAzJw5E61atcLly5cxefJkOcISERGZBKcJDFRcXIyEhAS5wxIREZGRSJ4mSEpK0nk8JydH72SIiIjMgaVcBSAXycVAREQEBEHQuWCQ8/9ERGTJrG3TIcnTBEqlEomJiVCr1XW2jIwMY+RJRETUZLjPQAMCAwORnp5e7/GGRg2IiIjIvEieJoiOjkZ5eXm9x/38/JCSkmJQUmRad/J/NHUKOjl5B5s6hXqZc26A/P+25v56ifRlbdMEkouB4GDdH35nZ2eEhITonRAREZGpWcrwvlxkv7SQiIiILItBOxASERE9iNRWtvaNxQAREVEN1lUKyFAMiKKI1NRUZGdnQ6lUIiwsDPb29nLkRkRERE1AcjEwcuRI7NixA66uriguLsbIkSNx/PhxuLu74+bNm3jooYfwww8/oE2bNsbIl4iIyOgs5Z4CcpG8gDA5ORkqlQoAMH/+fJSVleH333/HtWvXkJeXB2dnZyxcuFD2RImIiJqKKOMfS2DQNMH333+PlStXomPHjgCAtm3bYsWKFZgyZYrO81QqlaaguE8URW5jTEREZAJ6XVp4/5f2n3/+ic6dO2sd8/PzQ35+vs7z4+Li4OrqqtVEdZk+qRAREcmO2xE3wsSJE/HMM8+gqqoKubm5WscKCwvRsmVLnefHxMSgpKREqwk2LvqkQkREJDs1RNmaJZA8TRAZGan57/DwcFRUVGgd/89//oNevXrpjKFQKKBQKLT6OEVARETmwlLm+uUiuRjYtm2bzuOxsbGwtbXVOyEiIiJqWrJvR1xcXIzXX39d7rBERERNhmsGDFRcXIyEhAS5wxIRETUZURRla5ZA8jRBUlKSzuM5OTl6J0NERGTtNmzYgFWrVqGwsBABAQFYt24d+vfvX+djN2/ejM8++wxnzpwBAAQGBmLZsmX1Pr4+kouBiIgICIKgs9rhYkAiIrJkproKYNeuXYiKikJ8fDwGDBiAtWvXIiwsDFlZWfDw8Kj1+NTUVIwbNw4DBw6Eo6MjVqxYgWHDhuHs2bPw8fFp9PNKniZQKpVITEyEWq2us2VkZEgNSUREZFbkXDOgUqlQWlqq1WpuvHffmjVrMGXKFEyaNAndunVDfHw8mjVrhq1bt9b5+H/96194/fXX0atXL/j7++OTTz6BWq3G4cOHJb1eySMDgYGBSE9PR3h4eJ3HGxo1IOvj5B1s6hSIiEwmLi4Oixcv1uqLjY3FokWLtPoqKyuRnp6OmJgYTZ+NjQ1CQ0ORlpbWqOeqqKhAVVUVWrduLSlHycVAdHQ0ysvL6z3u5+eHlJQUqWGJiIjMhpz7DMTExCAqKkqrr+ZeOwBw48YNVFdXw9PTU6vf09MT58+fb9RzzZ07F97e3ggNDZWUo+RiIDhY97c8Z2dnhISESA1LRERkNuRcM1DXRnvGsHz5cuzcuROpqalwdHSUdK5BNyoiIiIiebi7u8PW1hZFRUVa/UVFRfDy8tJ57urVq7F8+XIcOnQIPXv2lPzcsu8zQEREZOlMsc+Ag4MDAgMDtRb/3V8MGBQUVO95K1euxNKlS5GcnIy+ffvq9XolFwNXrlzBjRs3ND//+OOPGD9+PIKDg/HSSy81epEDERGRuTLVDoRRUVHYvHkzEhIScO7cOUybNg3l5eWYNGkSAGDChAlaCwxXrFiBBQsWYOvWrejQoQMKCwtRWFiI27dvS3peycXAs88+i2PHjgEAvv76azz++OO4ffs2Bg0ahIqKCoSEhGDv3r1SwxIREZkNUcY/UowZMwarV6/GwoUL0atXL2RmZiI5OVmzqPDSpUsoKCjQPH7jxo2orKzEc889B6VSqWmrV6+W9LyCKPE6wObNm+O3335Dx44d8eijj+Lpp5/G3LlzNcfXr1+PrVu3St5vwM6h8ZsjkHHdyf9R1ni8tNB88N+WHhT3Kq8aNf6wdsNli/Xd5WTZYhmL5JEBOzs7lJWVAQByc3MxYsQIreMjRoxAVlaWPNkRERGZgBqibM0SSC4GQkJCsGPHDgBA7969kZqaqnU8JSWlwS0Q69qNiRsVERGRueCNihqwfPlyBAcHIz8/H4899hjmzZuHEydOoGvXrsjKysKuXbsQHx+vM0ZduzEJNs0h2LaQmg4REREZSPKaAQD4/fffMX/+fOzbt0+zYtHOzg79+vVDdHQ0IiIidJ6vUqlq7cvcys2fNzgyE5xXfnDx35YeFMZeMzCk7ZOyxUq5clC2WMai16ZDnTt3xo4dOyCKIq5duwa1Wg13d3fY29s36vy6dmNiIUBEROZCzu2ILYFBmw4JggBPT08olUpNIXD58mVMnjxZluSIiIjI+GTfgbC4uBgJCQlyhyUiImoyalGUrVkCydMESUlJOo/n5OTonQwREZE5sIxf4fKRXAxERERAEASdl0tw/p+IiMhySJ4mUCqVSExMhFqtrrNJ3XmQiIjI3HDToQYEBgYiPT293uMNjRoQERGZO2srBiRPE0RHR6O8vLze435+fkhJSTEoKZJG7mvHiYisnbV9qZVcDAQH695kxNnZGSEhIXonRERERE1Lr02HiIiIHmSWMrwvFxYDRERENXAHQiIiIrIqkouB9957D3l5ecbIhYiIyCxY2y2MJRcD0dHR6Ny5M5588kns2rULlZWVxsiLiIjIZKzt0kK9pgk++eQTODs74+WXX4a3tzdmzZqFM2fOyJ0bERERNQG9ioGRI0diz549uHLlCubMmYMDBw4gICAA/fv3x+bNm1FWVqbzfJVKhdLSUq1mKUMpRET04OM0gQQeHh6YM2cOzp07h9TUVHTr1g2zZ8+GUqnUeV5cXBxcXV21mqjWXUAQERE1FU4TNKC+mxAFBwfj008/RX5+Pt5//32dMWJiYlBSUqLVBBsXqakQERGRDCTvM9DQkEeLFi0wZcoUnY9RKBRQKBRafbzTIRERmQtr22dAcjGgVquNkQcREZHZUFvIXL9cZN906PLly5g8ebLcYYmIiJqMKOMfSyB7MVBcXIyEhAS5wxIREZGRSJ4mSEpK0nk8JydH72SIiIjMgbVNE0guBiIiIiAIgs6FhFwMSERElsxShvflInmaQKlUIjExEWq1us6WkZFhjDyJiIjISCSPDAQGBiI9PR3h4eF1Hm9o1ICAO/k/mjoFnZy8g02dAlkIud/LfO+RueA0QQOio6NRXl5e73E/Pz+kpKQYlBQREZEpWds0geRiIDhYd+Xu7OyMkJAQvRMiIiKipiW5GCAiInrQcZqAiIjIylnbNIFemw7t3bsXCxcuxE8//QQA+P777zFy5EgMHz4cmzZtkjVBIiIiMi7JxcDHH3+Mp59+Gvv378fIkSPxxRdfICIiAj4+PujQoQNmzZqFDz74wBi5EhERNQlRVMvWLIHkaYIPP/wQH330EaZMmYKUlBSMHDkS7733Hl5//XUAwKOPPoqVK1di5syZsidLRETUFNScJtAtNzcXYWFhAIAhQ4aguroagwcP1hx//PHHkZeXJ1+GRERETUwURdmaJZBcDLi5uWl+2efn5+PevXu4dOmS5nheXh5at26tM4ZKpUJpaalWs5S/MCIiogeN5GmC8PBwvPLKK4iMjERSUhImTJiAf/7zn7CxsYEgCIiOjsawYcN0xoiLi8PixYu1+gSb5hBsW0hNh4iISHbWNk0giBK/kpeXl2P27NlIS0vDwIEDsW7dOnz44YeYN28eqqqqEBISgl27dsHDw6PeGCqVCiqVSquvlZu/1dzgiNsRk6nwvUcPinuVV40a36fVI7LFuvrnWdliGYvkYqA+d+/eRVVVFVxcXPQ6387BR440LAL/h0ymwvcePShYDMhLr30G6uLo6AgXFxdcvnwZkydPlissERFRk1OLomzNEshWDNxXXFyMhIQEucMSERE1GVHGP5ZA8gLCpKQkncdzcnL0ToaIiIianuRiICIiAoIg6LwU0FoWAhIR0YPJ2i53lzxNoFQqkZiYCLVaXWfLyMgwRp5ERERNRg1RtmYJJBcDgYGBSE9Pr/d4Q6MGREREZF4kTxNER0ejvLy83uN+fn5ISUkxKKkHHS+fIlPhe4+ocaztS61s+wwYypr2GSAiIsMYe5+B1i5dZItVXHZRtljGInlkgIiI6EFnJt+Tm4zs+wwQERGRZeHIABERUQ2WchWAXPQqBu7cuYMdO3bg6NGjKCgogI2NDTp16oSIiAgMHTpU7hyJiIialLVNE0heQJidnY3Q0FDcuXMHCoUCV65cwciRI3Hjxg2cPHkSzzzzDLZv3w47O2l1BhcQEhFRYxl7AWEL506yxSotN/+deSWvGXjjjTcwfPhwFBYW4tKlS4iLi4NarcaxY8dw7tw5nDhxAu+8844xciUiImoS1najIskjA87OzsjMzESXLv+77KKyshLNmzdHQUEB3Nzc8PXXX2PWrFnIzc2tN4ZKpYJKpdLqa+Xmz22MiYioUYw9MuDcrINsscor/pAtlrFIHhlo2bIlysrKND9XVFTg3r17cHBwAAD07NkTBQUFOmPExcXB1dVVq4nqMp3nEBERkXFILgaefPJJREVF4fz588jNzcXUqVPRq1cvuLi4AAAuXboEDw8PnTFiYmJQUlKi1QQbF/1eARERkcysbZpA8tUEK1euRHh4OLp16wZBENCuXTvs3r1bc/z69euIjo7WGUOhUEChUGj1cYqAiIjMBa8maKSLFy9CpVLB399f8pUDdeHVBERE1FjGXjPg6Nhetlh3716SLZax6L0DYZcuXdC9e/dahcDly5cxefJkgxMjIiIyFVHGP5ZA9u2Ii4uLkZCQIHdYIiKiJiOKomxNqg0bNqBDhw5wdHTEgAEDcPz4cZ2P/+qrr+Dv7w9HR0f06NED+/fvl/ycksf3k5KSdB7PyTH/zRWIiIh0MdWagV27diEqKgrx8fEYMGAA1q5di7CwMGRlZdW5OP/nn3/GuHHjEBcXh7/97W/Yvn07IiIikJGRge7duzf6eSWvGbCxsYEgCDr/ogRBQHV1tZSwXDNARESNZuw1A/Yy/k6qkpDrgAED0K9fP6xfvx4AoFar0a5dO/zjH//AW2+9VevxY8aMQXl5Ofbu3avpe/TRR9GrVy/Ex8c3+nklTxMolUokJiZCrVbX2TIyMqSGJCIiMiuijE2lUqG0tFSr1dx4D/jfJn7p6ekIDQ3V9NnY2CA0NBRpaWl15pmWlqb1eAAICwur9/H1v2CJRo8eLS5YsKDe45mZmaIgCFLDNsrdu3fF2NhY8e7du2YVi/Ee7HjmnBvjmU8sxjO/eOYiNja2Vo0QGxtb63FXr14VAYg///yzVn90dLTYv3//OmPb29uL27dv1+rbsGGD6OHhISlHycXADz/8IH777bf1Hr99+7aYmpoqNWyjlJSUiADEkpISs4rFeA92PHPOjfHMJxbjmV88c3H37l2xpKREq9VV8JiyGJC8gDA4OFjncWdnZ4SEhEgNS0RE9ECqa6O9uri7u8PW1hZFRUVa/UVFRfDy8qrzHC8vL0mPr4/slxYSERGRdA4ODggMDMThw4c1fWq1GocPH0ZQUFCd5wQFBWk9HgAOHjxY7+PrY/jWgURERCSLqKgoREZGom/fvujfvz/Wrl2L8vJyTJo0CQAwYcIE+Pj4IC4uDgAwc+ZMhISE4L333sOoUaOwc+dOnDx5Eps2bZL0vBZVDCgUCsTGxjZquKUpYzHegx3PnHNjPPOJxXjmF88SjRkzBtevX8fChQtRWFiIXr16ITk5GZ6engD+dzNAG5v/G9QfOHAgtm/fjvnz5+Ptt99Gly5dsGfPHkl7DAAG3JuAiIiIHgxcM0BERGTlWAwQERFZORYDREREVo7FABERkZVjMUBERGTlLKYYkHp/Z11++OEHjB49Gt7e3hAEAXv27NE7VlxcHPr16wcXFxd4eHggIiICWVlZesfbuHEjevbsiRYtWqBFixYICgrCt99+q3e8v1q+fDkEQcCsWbP0jrFo0SIIgqDV/P399Y539epVvPTSS3Bzc4OTkxN69OiBkydP6hWrQ4cOtXITBAHTp0/XK151dTUWLFiAjh07wsnJCZ07d8bSpUsNurVpWVkZZs2aBV9fXzg5OWHgwIE4ceJEo85t6H0riiIWLlwIpVIJJycnhIaG4uLFi3rHS0xMxLBhw+Dm5gZBEJCZmalXrKqqKsydOxc9evSAs7MzvL29MWHCBOTn5+ud26JFi+Dv7w9nZ2e0atUKoaGh+OWXX/SO91dTp06FIAhYu3at3vEmTpxY6304fPhwg/I7d+4cnnrqKbi6usLZ2Rn9+vXDpUuX9IpX1+dEEASsWrVKcqzbt29jxowZaNu2LZycnNCtWzedd8trKF5RUREmTpwIb29vNGvWDMOHD9f5PiZ5WEQxcP/+zrGxscjIyEBAQADCwsJw7do1veKVl5cjICAAGzZsMDi3I0eOYPr06Th27BgOHjyIqqoqDBs2DOXl5XrFa9u2LZYvX4709HScPHkSTzzxBMLDw3H27FmD8jxx4gQ+/vhj9OzZ06A4APDII4+goKBA044ePapXnD///BODBg2Cvb09vv32W/z3v//Fe++9h1atWukV78SJE1p5HTx4EADw/PPP6xVvxYoV2LhxI9avX49z585hxYoVWLlyJdatW6dXPAB49dVXcfDgQXz++ef47bffMGzYMISGhuLq1YZvcdrQ+3blypX48MMPER8fj19++QXOzs4ICwvD3bt39YpXXl6Oxx57DCtWrDAot4qKCmRkZGDBggXIyMhAYmIisrKy8NRTT+n9Wh966CGsX78ev/32G44ePYoOHTpg2LBhuH79ul7x7tu9ezeOHTsGb29vnY9rTLzhw4drvR937Nihd7zff/8djz32GPz9/ZGamorTp09jwYIFcHR01CveX/MqKCjA1q1bIQgCnn32WcmxoqKikJycjC+++ALnzp3DrFmzMGPGDCQlJUnOTRRFREREICcnB19//TVOnToFX19fhIaG6v3/VGokSXcyMJH+/fuL06dP1/xcXV0tent7i3FxcQbHBiDu3r3b4Dj3Xbt2TQQgHjlyRLaYrVq1Ej/55BO9zy8rKxO7dOkiHjx4UAwJCRFnzpypd6zY2FgxICBA7/P/au7cueJjjz0mS6y6zJw5U+zcubOoVqv1On/UqFHi5MmTtfqeeeYZcfz48XrFq6ioEG1tbcW9e/dq9ffp00ecN2+epFg137dqtVr08vISV61apem7deuWqFAoxB07dkiO91e5ubkiAPHUqVN65VaX48ePiwDEvLw8WeLdv8HNoUOH9I535coV0cfHRzxz5ozo6+srvv/++w3Gqi9eZGSkGB4e3qjzGxNvzJgx4ksvvSRbvJrCw8PFJ554Qq9YjzzyiLhkyRKtvsa+p2vGy8rKEgGIZ86c0fRVV1eLbdq0ETdv3txgPNKf2Y8M6HN/Z1MqKSkBALRu3drgWNXV1di5cyfKy8sl7zP9V9OnT8eoUaNq3fNaXxcvXoS3tzc6deqE8ePH1ztU2ZCkpCT07dsXzz//PDw8PNC7d29s3rxZlhwrKyvxxRdfYPLkyRAEQa8YAwcOxOHDh3HhwgUAwK+//oqjR49ixIgResW7d+8eqqura32bc3Jy0nt05b7c3FwUFhZq/Ru7urpiwIABZvs5EQQBLVu2NDhWZWUlNm3aBFdXVwQEBOgVQ61W4+WXX0Z0dDQeeeQRg3MCgNTUVHh4eODhhx/GtGnTcPPmTb1z27dvHx566CGEhYXBw8MDAwYMMGh686+Kioqwb98+vPLKK3qdP3DgQCQlJeHq1asQRREpKSm4cOEChg0bJjmWSqUCAK3PiI2NDRQKhcGfEdLN7IuBGzduoLq6WrMV432enp4oLCw0UVZ1U6vVmDVrFgYNGiR5K8i/+u2339C8eXMoFApMnToVu3fvRrdu3fSKtXPnTmRkZGj2sTbUgAED8OmnnyI5ORkbN25Ebm4ugoODUVZWJjlWTk4ONm7ciC5duuDAgQOYNm0a3njjDSQkJBic5549e3Dr1i1MnDhR7xhvvfUWxo4dC39/f9jb26N3796YNWsWxo8fr1c8FxcXBAUFYenSpcjPz0d1dTW++OILpKWloaCgQO88AWg+C5bwObl79y7mzp2LcePGoUWLFnrH2bt3L5o3bw5HR0e8//77OHjwINzd3fWKtWLFCtjZ2eGNN97QO5+/Gj58OD777DMcPnwYK1aswJEjRzBixAhUV1dLjnXt2jXcvn0by5cvx/Dhw/Hdd9/h6aefxjPPPIMjR44YnGtCQgJcXFzwzDPP6HX+unXr0K1bN7Rt2xYODg4YPnw4NmzYgMGDB0uO5e/vj/bt2yMmJgZ//vknKisrsWLFCly5csXgzwjpZlH3JjB306dPx5kzZwyuYB9++GFkZmaipKQE//73vxEZGYkjR45ILgguX76MmTNn4uDBg/XOLUr112/FPXv2xIABA+Dr64svv/xS8jcLtVqNvn37YtmyZQCA3r1748yZM4iPj0dkZKRBeW7ZsgUjRoxocO5Xly+//BL/+te/sH37djzyyCPIzMzErFmz4O3trXd+n3/+OSZPngwfHx/Y2tqiT58+GDduHNLT0/XO05JUVVXhhRdegCiK2Lhxo0GxhgwZgszMTNy4cQObN2/GCy+8gF9++QUeHh6S4qSnp+ODDz5ARkaG3qNINY0dO1bz3z169EDPnj3RuXNnpKamYujQoZJiqdVqAEB4eDhmz54NAOjVqxd+/vlnxMfHG3zL+K1bt2L8+PF6/z9i3bp1OHbsGJKSkuDr64sffvgB06dPh7e3t+TRSHt7eyQmJuKVV15B69atYWtri9DQUIwYMcKghbvUMLMfGdDn/s6mMGPGDOzduxcpKSlo27atQbEcHBzg5+eHwMBAxMXFISAgAB988IHkOOnp6bh27Rr69OkDOzs72NnZ4ciRI/jwww9hZ2en17eUmlq2bImHHnoI2dnZks9VKpW1CpyuXbvqPe1wX15eHg4dOoRXX33VoDjR0dGa0YEePXrg5ZdfxuzZsw0aZencuTOOHDmC27dv4/Llyzh+/DiqqqrQqVMng3K9/1kw58/J/UIgLy8PBw8eNGhUAACcnZ3h5+eHRx99FFu2bIGdnR22bNkiOc6PP/6Ia9euoX379prPSV5eHv75z3+iQ4cOBuV4X6dOneDu7q7X58Td3R12dnZG+az8+OOPyMrK0vuzcufOHbz99ttYs2YNRo8ejZ49e2LGjBkYM2YMVq9erVfMwMBAZGZm4tatWygoKEBycjJu3rxp8GeEdDP7YkCf+zs3JVEUMWPGDOzevRvff/89OnbsKPtzqNVqzVyaFEOHDsVvv/2GzMxMTevbty/Gjx+PzMxM2NraGpzb7du38fvvv0OpVEo+d9CgQbUuw7xw4QJ8fX0Nymnbtm3w8PDAqFGjDIpTUVGhdXcwALC1tdV8UzOEs7MzlEol/vzzTxw4cADh4eEGxevYsSO8vLy0PielpaX45ZdfzOJzcr8QuHjxIg4dOgQ3NzfZn0Pfz8nLL7+M06dPa31OvL29ER0djQMHDsiS25UrV3Dz5k29PicODg7o16+fUT4rW7ZsQWBgoN5rLaqqqlBVVWWUz4mrqyvatGmDixcv4uTJkwZ/Rkg3i5gmaOj+zlLdvn1bq0LPzc1FZmYmWrdujfbt20uKNX36dGzfvh1ff/01XFxcNPOzrq6ucHJykpxbTEwMRowYgfbt26OsrAzbt29HamqqXv9TcnFxqbV2wdnZGW5ubnqvaXjzzTcxevRo+Pr6Ij8/H7GxsbC1tcW4ceMkx5o9ezYGDhyIZcuW4YUXXsDx48exadMmyffh/iu1Wo1t27YhMjISdnaGvb1Hjx6Nd999F+3bt8cjjzyCU6dOYc2aNZg8ebLeMQ8cOABRFPHwww8jOzsb0dHR8Pf3b9R7uaH37axZs/DOO++gS5cu6NixIxYsWABvb29EREToFa+4uBiXLl3S7Adw/5eRl5dXrdEGXbGUSiWee+45ZGRkYO/evaiurtZ8Tlq3bg0HBwdJubm5ueHdd9/FU089BaVSiRs3bmDDhg24evVqvZeRNvRaaxYn9vb28PLywsMPPyw5XuvWrbF48WI8++yz8PLywu+//445c+bAz88PYWFheuUXHR2NMWPGYPDgwRgyZAiSk5PxzTffIDU1Va94wP+Kxa+++grvvfdenTEaGyskJATR0dFwcnKCr68vjhw5gs8++wxr1qzRK95XX32FNm3aoH379vjtt98wc+ZMRERE6LUgkSQw6bUMEqxbt05s37696ODgIPbv3188duyY3rFSUlJEALVaZGSk5Fh1xQEgbtu2Ta/cJk+eLPr6+ooODg5imzZtxKFDh4rfffedXrHqYuilhWPGjBGVSqXo4OAg+vj4iGPGjBGzs7P1jvfNN9+I3bt3FxUKhejv7y9u2rRJ71iiKIoHDhwQAYhZWVkGxRFFUSwtLRVnzpwptm/fXnR0dBQ7deokzps3T1SpVHrH3LVrl9ipUyfRwcFB9PLyEqdPny7eunWrUec29L5Vq9XiggULRE9PT1GhUIhDhw7V+ffQULxt27bVeTw2NlZSrPuXJtbVUlJSJOd2584d8emnnxa9vb1FBwcHUalUik899ZR4/PhxvV9rTQ1dWqgrXkVFhThs2DCxTZs2or29vejr6ytOmTJFLCwsNCi/LVu2iH5+fqKjo6MYEBAg7tmzx6B4H3/8sejk5NTg+6+hWAUFBeLEiRNFb29v0dHRUXz44YfF9957r95LehuK98EHH4ht27YV7e3txfbt24vz58836DNHjSOIIldlEBERWTOzXzNARERExsVigIiIyMqxGCAiIrJyLAaIiIisHIsBIiIiK8digIiIyMqxGCAiIrJyLAaIiIisHIsBIiIiK8digIiIyMqxGCAiIrJy/w9kSZu+rR9rfgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgMAAAGzCAYAAACy+RS/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABKPUlEQVR4nO3deVxUZf8//tewDYiKKQKDCy6YaCoYKqEimSQuIdhdruVC2VfTErlDxQ1tcaTStFy4NU1aXKobjdQwNDBNTAW5zVLUwJVFkQREHZA5vz/6OZ+GZeDMnIEZ5/XscT0ecZb3vA/MyJvrus51ZIIgCCAiIiKLZdXYCRAREVHjYjFARERk4VgMEBERWTgWA0RERBaOxQAREZGFYzFARERk4VgMEBERWTgWA0RERBaOxQAREZGFYzFA9AhJTU2FTCbDt99+29ipEJEZYTFAJKGtW7dCJpNBJpPhyJEj1fYLgoB27dpBJpPhueeea4QMiYiqYzFAZAT29vbYtm1bte2HDh3CtWvXIJfLGyErIqKasRggMoIRI0bgm2++wYMHD7S2b9u2Db6+vnBzc2ukzIiIqmMxQGQE48ePx61bt5CcnKzZVl5ejm+//RYTJkyodvyHH36I/v37o1WrVnBwcICvr2+N4/7JyckYOHAgWrRogaZNm6Jr165YsGCBzlxUKhWee+45ODk54ejRo4ZfHBE9clgMEBlBhw4d4O/vj+3bt2u2/fDDDyguLsa4ceOqHb9mzRr07t0bb7/9NpYvXw4bGxu8+OKL2Lt3r+aY33//Hc899xxUKhXefvttrFy5EqNGjcIvv/xSax737t1DSEgIjh49igMHDqB///7SXigRPRJsGjsBokfVhAkTEB0djXv37sHBwQFfffUVAgMD4e7uXu3Y8+fPw8HBQfP1rFmz8OSTT2LVqlUYOXIkgL97BcrLy/HDDz/A2dm5zte/c+cOnnvuOfz+++/46aef4OPjI9m1EdGjhT0DREYyZswY3Lt3D3v27EFpaSn27NlT4xABAK1C4K+//kJxcTECAgKQkZGh2d6iRQsAwHfffQe1Wq3ztYuLizF06FCcO3cOqampLASISCcWA0RG0rp1awQFBWHbtm1ISEhAZWUlXnjhhRqP3bNnD5566inY29ujZcuWaN26NTZs2IDi4mLNMWPHjsWAAQPw6quvwtXVFePGjcPXX39dY2EQERGBEydO4MCBA3jiiSeMdo1E9GhgMUBkRBMmTMAPP/yAuLg4DB8+XPPX/T8dPnwYo0aNgr29PdavX499+/YhOTkZEyZMgCAImuMcHBzw888/48CBA3j55Zdx+vRpjB07Fs8++ywqKyu1YoaGhkIQBKxYsaLOXgQiIhYDREY0evRoWFlZ4dixY7UOEfz3v/+Fvb099u/fj/DwcAwfPhxBQUE1HmtlZYUhQ4Zg1apV+OOPP/Dee+/hp59+QkpKitZxYWFh2LJlC7Zt24aZM2dKfl1E9GjhBEIiI2ratCk2bNiAS5cuISQkpMZjrK2tIZPJtP66v3TpEnbv3q11XFFREVq2bKm17eFcAJVKVS3upEmTUFJSgjfeeAPNmzdHbGysYRdDRI8sFgNERjZ58mSd+0eOHIlVq1Zh2LBhmDBhAm7cuIF169bB09MTp0+f1hz39ttv4+eff8bIkSPh4eGBGzduYP369Wjbti0GDhxYY+xZs2ahpKQECxcuhJOTU51rEhCRZWIxQNTInnnmGWzevBkrVqxAREQEOnbsiNjYWFy6dEmrGBg1ahQuXbqELVu2oLCwEM7OzggMDMSyZcvg5ORUa/wFCxaguLhYUxBw2ICIqpIJ/5yhRERERBaHEwiJiIgsHIsBIiIiC8digIiIyMKxGCAiIjIRP//8M0JCQuDu7g6ZTFbtFuOapKam4sknn4RcLoenpye2bt0q+nVZDBAREZmIsrIyeHt7Y926dfU6PicnByNHjsTgwYORmZmJiIgIvPrqq9i/f7+o1+XdBERERCZIJpNh165dCAsLq/WYefPmYe/evThz5oxm27hx43D79m0kJSXV+7XYM0BERGREKpUKJSUlWq2mVUP1kZaWVm358uDgYKSlpYmKYzKLDlUUZksaz8E9QLJY93IPSxaLDMefLTUGKd93xiD1e1nq65U6P1vnTpLGq0rK30nKtZ9j2bJlWttiYmKwdOlSg2Pn5+fD1dVVa5urqytKSkpw7949rcej62IyxQAREZHJUFfWfUw9RUdHIzIyUmubXC6XLL4UWAwQEREZkVwuN9ovfzc3NxQUFGhtKygoQPPmzevdKwCwGCAiIqpOUDd2BvXi7++Pffv2aW1LTk6Gv7+/qDicQEhERFSVWi1dE+HOnTvIzMxEZmYmgL9vHczMzMSVK1cA/D3kMGnSJM3x06dPR3Z2NubOnYtz585h/fr1+PrrrzFnzhxRr8ueASIioiqERuoZOHnyJAYPHqz5+uFcg8mTJ2Pr1q3Iy8vTFAYA0LFjR+zduxdz5szBmjVr0LZtW3z66acIDg4W9bosBoiIiEzE008/DV3L/9S0uuDTTz+NU6dOGfS6oouBwsJCbNmyBWlpacjPzwfw9wSG/v37Y8qUKWjdurVBCRERETU6kd375k5UMXDixAkEBwejSZMmCAoKwuOPPw7g75mLH3/8MVasWIH9+/ejT58+OuOoVKpqCy5YqVQmd6sFERFZKDOZQCgVUcXAG2+8gRdffBFxcXGQyWRa+wRBwPTp0/HGG2/UufKRUqmstgDDoqg3sWTubDHpEBERkQREFQP/+9//sHXr1mqFAPD3Gspz5sxB796964xT0wIMVqXXxaRCRERkPBIuOmQORBUDbm5uOH78OLy8vGrcf/z48WrLItakpgUYKsoLxaRCRERkPBwmqN1bb72F1157Denp6RgyZIjmF39BQQEOHjyITZs24cMPPzRKokRERGQcooqBmTNnwtnZGR999BHWr1+Pysq/u1Gsra3h6+uLrVu3YsyYMUZJlIiIqMHwbgLdxo4di7Fjx6KiogKFhX937Ts7O8PW1lby5IiIiBpDYy061Fj0XnTI1tYWCoVCylyIiIioEXAFQiIioqo4TEBERGThOEzwaLiXe7ixU6D/n4N7gKTxpPzZmnJuZFr4s7UwFrbOAB9hTEREZOEe2Z4BIiIivXGYgIiIyMJZ2ARCDhMQERFZOPYMEBERVcVhAiIiIgvHYQIiIiKyJJIXA1evXkV4eLjOY1QqFUpKSrSaSqWSOhUiIiK9CEKlZM0cSF4MFBUVIT4+XucxSqUSTk5OWi12TZzUqRAREelHUEvXzIDoOQOJiYk692dnZ9cZIzo6GpGRkVrbrEqvi02FiIiIJCC6GAgLC4NMJoMgCLUeI5PJdMaQy+WQy+Va2yrKC8WmQkREZBycQKibQqFAQkIC1Gp1jS0jI8MYeRIRETUcCxsmEF0M+Pr6Ij09vdb9dfUaEBERmTx1pXTNDIgeJoiKikJZWVmt+z09PZGSkmJQUkRERNRwRBcDAQG6H/nq6OiIwMBAvRMiIiJqdGbSvS8VrkBIRERUFScQEhERkSVhzwAREVFVHCZ4NDi4657bIMa93MOSxbJEpvz9M+XcAGnfx8Yg9ffPkj63/Nka5kG5kReq4zABERERWZJHtmeAiIhIbxbWM8BigIiIqApzedqgVDhMQEREZOHYM0BERFQVhwmIiIgsHG8tJCIisnAW1jMges7AvXv3cOTIEfzxxx/V9t2/fx+ff/55nTFUKhVKSkq0mkqlEpsKERERSUBUMXD+/Hl069YNgwYNQs+ePREYGIi8vDzN/uLiYkydOrXOOEqlEk5OTlotdk2c+OyJiIiMQVBL18yAqGJg3rx56NGjB27cuIGsrCw0a9YMAwYMwJUrV0S9aHR0NIqLi7XavNnTRcUgIiIyGrVaumYGRM0ZOHr0KA4cOABnZ2c4Ozvj+++/x+uvv46AgACkpKTA0dGxXnHkcjnkcrnWtoryQjGpEBERkURE9Qzcu3cPNjb/Vz/IZDJs2LABISEhCAwMxPnz5yVPkIiIqMFZ2DCBqJ4BLy8vnDx5Et26ddPavnbtWgDAqFGjpMuMiIiosZhJ975URPUMjB49Gtu3b69x39q1azF+/HgIgiBJYkRERNQwRBUD0dHR2LdvX637169fD7WFVVNERPQI4gRCIiIiC2cmY/1S4YOKiIiILNwj2zNwL/dwY6dAZDBLex9b0vVKfa0O7gEmHc/sfrZm0r0vlUe2GCAiItKbhQ0TsBggIiKqysJ6BjhngIiIyMKxZ4CIiKgqDhMQERFZOA4TEBERkSVhzwAREVFVFtYzILoYOHv2LI4dOwZ/f394eXnh3LlzWLNmDVQqFV566SU888wzdcZQqVRQqVRa26xUqmqPNSYiImoUFvacHVHDBElJSfDx8cFbb72F3r17IykpCYMGDcLFixdx+fJlDB06FD/99FOdcZRKJZycnLRa7Jo4vS+CiIiI9CcTRDxmsH///njmmWfw7rvvYseOHXj99dcxY8YMvPfeewD+fpBReno6fvzxR51xauwZKL3OngEiIj1JvWKg1KRegdDWuZOk8aq6tz1GslgO45dJFstYRPUM/P7775gyZQoAYMyYMSgtLcULL7yg2T9x4kScPn26zjhyuRzNmzfXaiwEiIjIZFjYUwtF300gk8n+PtHKCvb29nByctLsa9asGYqLi6XLjoiIiIxOVDHQoUMHXLhwQfN1Wloa2rdvr/n6ypUrUCgU0mVHRETUGAS1dM0MiLqbYMaMGaisrNR83aNHD639P/zwQ73uJiAiIjJpZtK9LxVRPQPTp0/HyJEja92/fPlyfPrppwYnRURE1KgEQbom0rp169ChQwfY29vDz88Px48f13n86tWr0bVrVzg4OKBdu3aYM2cO7t+/L+o1uQIhERGRidi5cyciIyMRExODjIwMeHt7Izg4GDdu3Kjx+G3btmH+/PmIiYnB2bNnsXnzZuzcuRMLFiwQ9bosBoiIiKpqpLsJVq1ahWnTpmHq1Kno3r074uLi0KRJE2zZsqXG448ePYoBAwZgwoQJ6NChA4YOHYrx48fX2ZtQFYsBIiKiqiQsBlQqFUpKSrRa1bV2AKC8vBzp6ekICgrSbLOyskJQUBDS0tJqTLN///5IT0/X/PLPzs7Gvn37MGLECFGXy2cT1IOlLeZB1Fik/KxZ2udC6us19X/3zIlSqcSyZdoLD8XExGDp0qVa2woLC1FZWQlXV1et7a6urjh37lyNsSdMmIDCwkIMHDgQgiDgwYMHmD59OocJiIiIDCbhrYXR0dEoLi7WatHR0ZKkmZqaiuXLl2P9+vXIyMhAQkIC9u7di3feeUdUHPYMEBERVSGopXtQkVwur9cqu87OzrC2tkZBQYHW9oKCAri5udV4zuLFi/Hyyy/j1VdfBQD07NkTZWVleO2117Bw4UJYWdXvb372DBAREZkAOzs7+Pr64uDBg5ptarUaBw8ehL+/f43n3L17t9ovfGtrawCAiEcPsWeAiIiomkZadCgyMhKTJ09Gnz590K9fP6xevRplZWWYOnUqAGDSpElo06YNlEolACAkJASrVq1C79694efnh4sXL2Lx4sUICQnRFAX1wWKAiIioqkZaRnjs2LG4efMmlixZgvz8fPj4+CApKUkzqfDKlStaPQGLFi2CTCbDokWLcP36dbRu3RohISGapwnXl6hHGNdGEATNA4z0VVGYbWgaRmPqs2otbdY0Pbp4N4HpkPrfPXN7hPHdDW9IFqvJjE8ki2UskswZkMvlOHv2rBShiIiIGp9akK6ZAVHDBJGRkTVur6ysxIoVK9CqVSsAf6+gpItKpaq24IKVSlWv2ZZERERGZ2EPKhJVDKxevRre3t5o0aKF1nZBEHD27Fk4OjrWa7igpgUYFkW9iSVzZ4tJh4iIyDhYDNRu+fLl2LhxI1auXKn1qGJbW1ts3boV3bt3r1ec6Ojoar0MVqXXxaRCREREEhFVDMyfPx9DhgzBSy+9hJCQECiVStja2op+0ZoWYKgoLxQdh4iIyCgMn1tvVkRPIOzbty/S09Nx8+ZN9OnTB2fOnDH4TgIiIiKT0khPLWwseq0z0LRpU8THx2PHjh0ICgpCZWWl1HkRERFRAzFo0aFx48Zh4MCBSE9Ph4eHh1Q5ERERNS4zuSVQKgavQNi2bVu0bdtWilyIiIhMQyOtQNhY+KAiIiIiC8dnExAREVXFYQKqimucmxauX68/U19vXsp4fKaIYaTOT+qfx4Ny465NI5jJXQBS4TABERGRhWPPABERUVUcJiAiIrJwFnY3AYsBIiKiqiysZ4BzBoiIiCwcewaIiIiqsrC7CVgMEBERVWVhwwQGFQNlZWX4+uuvcfHiRSgUCowfPx6tWrWq8zyVSgWVSqW1zUqlqvZYYyIiIjI+UXMGunfvjqKiIgDA1atX0aNHD8yZMwfJycmIiYlB9+7dkZOTU2ccpVIJJycnrRa7Jk6/KyAiIpKaoJaumQFRxcC5c+fw4MEDAEB0dDTc3d1x+fJlHD9+HJcvX0avXr2wcOHCOuNER0ejuLhYq82bPV2/KyAiIpKaWpCumQG9hwnS0tIQFxcHJycnAEDTpk2xbNkyjBs3rs5z5XJ5tSGBivJCfVMhIiIiA4guBmQyGQDg/v37UCgUWvvatGmDmzdvSpMZERFRI7G0ZxOILgaGDBkCGxsblJSUICsrCz169NDsu3z5cr0mEBIREZk0M+nel4qoYiAmJkbr66ZNm2p9/f333yMgwLSfFEZERETaDCoGqvrggw8MSoaIiMgksGeAiIjIwpnJLYFSYTFARERUlYX1DPBBRURERBaOPQP14OAu7aTIe7mHJY1naUz5+2fq7xVT/t5JzZKulaQnWFjPAIsBIiKiqiysGOAwARERkYVjzwAREVFVXIGQiIjIwnGYgIiIiCwJewaIiIiqsrCeARYDREREVQiCZRUDooYJMjIykJOTo/n6iy++wIABA9CuXTsMHDgQO3bsqFcclUqFkpISraZSqcRlTkRERJIQVQxMnToVf/75JwDg008/xf/7f/8Pffr0wcKFC9G3b19MmzYNW7ZsqTOOUqmEk5OTVotdE6ffFRAREUlNLUjXzICoYYILFy6gS5cuAID169djzZo1mDZtmmZ/37598d577yE8PFxnnOjoaERGRmptsyq9LiYVIiIi4zGTX+JSEVUMNGnSBIWFhfDw8MD169fRr18/rf1+fn5awwi1kcvlkMvlWtsqygvFpEJERGQ0lrYcsahhguHDh2PDhg0AgMDAQHz77bda+7/++mt4enpKlx0REREZnaiegdjYWAwYMACBgYHo06cPVq5cidTUVHTr1g1ZWVk4duwYdu3aZaxciYiIGgZ7Bmrn7u6OU6dOwd/fH0lJSRAEAcePH8ePP/6Itm3b4pdffsGIESOMlSsREVHDUEvYzIDodQZatGiBFStWYMWKFcbIh4iIiBoYFx0iIiKqwtImELIYICIiqorFwKPBwT1Aslj3cg9LFosebXyvUH1J+W8UYPrvPVPPz9I9ssUAERGR3sxk4p9UWAwQERFVYWlzBkTdWkhERESPHvYMEBERVcVhAiIiIstmacMELAaIiIiqsrCeAc4ZICIisnCiioE33ngDhw8bfq+oSqVCSUmJVlOpVAbHJSIikoKglq6ZA1HFwLp16/D000/j8ccfR2xsLPLz8/V6UaVSCScnJ60WuyZOr1hERESSs7AHFYkeJvjxxx8xYsQIfPjhh2jfvj1CQ0OxZ88eqNX1v+Lo6GgUFxdrtXmzp4tNhYiIiCQguhjo2bMnVq9ejdzcXHz55ZdQqVQICwtDu3btsHDhQly8eLHOGHK5HM2bN9dqcrlcrwsgIiKSGocJ6snW1hZjxoxBUlISsrOzMW3aNHz11Vfo2rWrlPkRERE1PA4TiNe+fXssXboUOTk5SEpKkiIkERERNRBR6wx4eHjA2tq61v0ymQzPPvuswUkRERE1JnPp3peKqJ6BnJwctGrVyli5EBERmYTGnDOwbt06dOjQAfb29vDz88Px48d1Hn/79m3MnDkTCoUCcrkcjz/+OPbt2yfqNbkCIRERURWN1TOwc+dOREZGIi4uDn5+fli9ejWCg4ORlZUFFxeXaseXl5fj2WefhYuLC7799lu0adMGly9fRosWLUS9LosBIiIiE7Fq1SpMmzYNU6dOBQDExcVh79692LJlC+bPn1/t+C1btqCoqAhHjx6Fra0tAKBDhw6iX5fLERMREVUlyCRr9V11t7y8HOnp6QgKCtJss7KyQlBQENLS0mpMMzExEf7+/pg5cyZcXV3Ro0cPLF++HJWVlaIu95HtGbiXa/iyyebCwT2gsVPQyZJ+FqaO7xX9Sf29M+VrJWmHCZRKJZYtW6a1LSYmBkuXLtXaVlhYiMrKSri6umptd3V1xblz52qMnZ2djZ9++gkTJ07Evn37cPHiRbz++uuoqKhATExMvXN8ZIsBIiIiUxAdHY3IyEitbVIttKdWq+Hi4oKNGzfC2toavr6+uH79Oj744AMWA0RERIYQ1DLJYsnl8nr98nd2doa1tTUKCgq0thcUFMDNza3GcxQKBWxtbbVu++/WrRvy8/NRXl4OOzu7euXIOQNERERVNMathXZ2dvD19cXBgwc129RqNQ4ePAh/f/8azxkwYAAuXryo9Xyg8+fPQ6FQ1LsQAFgMEBERmYzIyEhs2rQJ8fHxOHv2LGbMmIGysjLN3QWTJk1CdHS05vgZM2agqKgIs2fPxvnz57F3714sX74cM2fOFPW6HCYgIiKqQhCkGyYQY+zYsbh58yaWLFmC/Px8+Pj4ICkpSTOp8MqVK7Cy+r+/49u1a4f9+/djzpw56NWrF9q0aYPZs2dj3rx5ol5XJgiCIOmV6KmiMLuxUzBbnCFO9cX3iv54N4FpsXXuZNT41/yekSxW219/kiyWsXCYgIiIyMKJLgbWrl2LSZMmYceOHQCAL774At27d4eXlxcWLFiABw8e1BmjvgswEBERNQZBLZOsmQNRxcC7776LBQsW4O7du5gzZw5iY2MxZ84cTJw4EZMnT8ann36Kd955p844SqUSTk5OWi12TZzeF0FERCQlQZCumQNREwi3bt2KrVu34vnnn8f//vc/+Pr6Ij4+HhMnTgQAeHl5Ye7cudVWWqqqpgUYrEqvi0ydiIjIOMzlL3qpiCoGcnNz0adPHwCAt7c3rKys4OPjo9n/5JNPIjc3t844NS3AUFFeKCYVIiIikoioYQI3Nzf88ccfAIALFy6gsrJS8zUA/P777zU+YpGIiMicWNqcAVE9AxMnTsSkSZMQGhqKgwcPYu7cuXjrrbdw69YtyGQyvPfee3jhhReMlSsREVGDMJexfqmIKgaWLVsGBwcHpKWlYdq0aZg/fz68vb0xd+5c3L17FyEhIfWaQEhERESmQ1QxYGVlhQULFmhtGzduHMaNGydpUkRERI3JXLr3pcLliImIiKporOWIGwtXICQiIrJwj2zPgJTriJv6GuKmnh89uizpvSf1tfJZB6ZNzKOHHwWPbDFARESkLzWHCYiIiMiSsGeAiIioCkubQMhigIiIqAreWkhERGThLG0FQs4ZICIisnCiewby8vKwYcMGHDlyBHl5ebCyskKnTp0QFhaGKVOmwNra2hh5EhERNRhLGyYQ1TNw8uRJdOvWDfv27UNFRQUuXLgAX19fODo64q233sKgQYNQWlpaZxyVSoWSkhKtplKp9L4IIiIiKakFmWTNHIgqBiIiIjBnzhycPHkShw8fxtatW3H+/Hns2LED2dnZuHv3LhYtWlRnHKVSCScnJ60WuyZO74sgIiIi/ckEof7TJJo0aYIzZ86gU6dOAAC1Wg17e3tcvXoVrq6uSE5OxpQpU3D9+nWdcVQqVbWeAKvS65DL5XpcQs0saQVCovriqnemgz8Lw9g6dzJq/N86hkgWq2fO95LFMhZRcwZcXFyQl5enKQYKCgrw4MEDNG/eHADQpUsXFBUV1RlHLpdX+8VfUV4oJhUiIiKj4d0EOoSFhWH69OlISkpCSkoKJk6ciMDAQDg4OAAAsrKy0KZNG6MkSkRERMYhqmfg3XffRV5eHkJCQlBZWQl/f398+eWXmv0ymQxKpVLyJImIiBqSuUz8k4qoYqBp06bYuXMn7t+/jwcPHqBp06Za+4cOHSppckRERI2ByxHXg729vdR5EBERUSPhcsRERERVWNoEQhYDREREVXDOwCPC0u65pUeTpd2Lbknrg5h6fpbO0uYM8EFFREREFu6R7RkgIiLSF4cJiIiILJyFzR/kMAEREZGl06tnoLy8HLt370ZaWhry8/MBAG5ubujfvz9CQ0NhZ2cnaZJEREQNydKGCUT3DFy8eBHdunXD5MmTcerUKajVaqjVapw6dQqTJk3CE088gYsXLxojVyIiogYhCDLJmjkQ3TMwY8YM9OzZE6dOndI8rfChkpISTJo0CTNnzsT+/fslS5KIiIiMR3Qx8Msvv+D48ePVCgEAaN68Od555x34+flJkhwREVFjUDd2Ag1MdDHQokULXLp0CT169Khx/6VLl9CiRQudMVQqFVQqldY2K5UKcrlcbDpERESSE2Ae3ftSET1n4NVXX8WkSZPw0Ucf4fTp0ygoKEBBQQFOnz6Njz76CFOmTMFrr72mM4ZSqYSTk5NWi10Tp/dFEBERkf5kgiD+cQyxsbFYs2YN8vPzIZP9XT0JggA3NzdERERg7ty5Os+vsWeg9Dp7Boiq4HLE+jP1ayXD2Dp3Mmr8VNcXJYv1dME3ksUyFr1uLZw3bx7mzZuHnJwcrVsLO3bsWK/z5XJ5tV/8FeWF+qRCREQkOTWHCeqvY8eO8Pf3h7+/v6YQuHr1KsLDwyVJjoiIqDEIkEnWzIHkKxAWFRUhPj5e6rBERERkJKKHCRITE3Xuz87O1jsZIiIiU8BbC+sQFhYGmUwGXfMOH04qJCIiMkfm0r0vFdHDBAqFAgkJCZpliKu2jIwMY+RJRERERiK6GPD19UV6enqt++vqNSAiIjJ1agmbORA9TBAVFYWysrJa93t6eiIlJcWgpEyNpd3rTfqT+r1iafjZIFNhLr/EpSK6GAgI0P2PnaOjIwIDA/VOiIiIiBqWXosOERERPcosbQIhiwEiIqIq1JZVC0i/6BARERGZF8mLgYKCArz99ttShyUiImowasgka+ZA8mIgPz8fy5YtkzosERFRgxEkbOZA9JyB06dP69yflZWldzJERESmgLcW1sHHx6fWhYUebq9rOWKVSgWVSqW1zUqlqvZYYyIiIjI+0cMELVu2xKZNm5CTk1OtZWdnY8+ePXXGUCqVcHJy0mqxa+L0ugAiIiKpqWUyyZo5EN0z4Ovri9zcXHh4eNS4//bt23UuRxwdHY3IyEitbVal18WmQkREZBTmMtYvFdHFwPTp03UuR9y+fXt89tlnOmPI5fJqQwIV5YViUyEiIiIJiC4GRo8erXP/Y489hsmTJ+udEBERUWOztAmEkt9aePXqVYSHh0sdloiIqMGoZdI1cyB5MVBUVIT4+HipwxIREZGRiB4mSExM1Lk/Oztb72SIiIhMQWOuHLhu3Tp88MEHyM/Ph7e3Nz755BP069evzvN27NiB8ePHIzQ0FLt37xb1mqKLgbCwsFrXGXiornUGiIiITFlj3U2wc+dOREZGIi4uDn5+fli9ejWCg4ORlZUFFxeXWs+7dOkS3nrrLQQEBOj1uqKHCRQKBRISEqBWq2tsGRkZeiVCRET0KFKpVCgpKdFqVRfee2jVqlWYNm0apk6diu7duyMuLg5NmjTBli1bao1fWVmJiRMnYtmyZejUqZNeOeq1zkB6ejpCQ0Nr3F9XrwEBDu76VW61uZd7WNJ4UrO065WS1Ncq9c9Capb0syXTJuXEP6VSWe2ZPTExMVi6dKnWtvLycqSnpyM6OlqzzcrKCkFBQUhLS6s1/ttvvw0XFxe88sorOHxYv8+Q6GIgKipK5zoDnp6eSElJ0SsZIiIiUyDlrYU1LbRX0/L7hYWFqKyshKurq9Z2V1dXnDt3rsbYR44cwebNm5GZmWlQjqKLgbrGIxwdHREYGKh3QkRERI1Nyv7tmhbak0JpaSlefvllbNq0Cc7OzgbFEl0MEBERkfScnZ1hbW2NgoICre0FBQVwc3Ordvyff/6JS5cuISQkRLNNrf67T8PGxgZZWVno3LlzvV5b8nUGiIiIzF1jLDpkZ2cHX19fHDx48P/yUKtx8OBB+Pv7Vzvey8sLv/32GzIzMzVt1KhRGDx4MDIzM9GuXbt6v7bePQPXrl1DixYt0LRpU63tFRUVSEtLw6BBg/QNTURE1KgaazniyMhITJ48GX369EG/fv2wevVqlJWVYerUqQCASZMmoU2bNlAqlbC3t0ePHj20zm/RogUAVNteF9HFQF5eHkJDQ5Geng6ZTIYJEyZg/fr1mqKgqKgIgwcPRmVlpdjQREREFm3s2LG4efMmlixZgvz8fPj4+CApKUkzqfDKlSuwspK+U190MTB//nxYWVnh119/xe3btzF//nwMHjwYP/74Ix577DEA4K2FRERk1hrzQUWzZs3CrFmzatyXmpqq89ytW7fq9Zqiy4sDBw7g448/Rp8+fRAUFIRffvkFCoUCzzzzDIqKigBwBUIiIjJvgky6Zg5EFwPFxcWaHgDg71smEhIS0KFDBwwePBg3btyoM4aY1ZiIiIjIuEQXA506dcLp06e1ttnY2OCbb75Bp06d8Nxzz9UZQ6lUwsnJSavFrokTmwoREZFRqCVs5kB0MTB8+HBs3Lix2vaHBYGPj0+dcwaio6NRXFys1ebNni42FSIiIqOwtGJA9ATC9957D3fv3q05mI0N/vvf/+L69es6Y9S0GlNFeaHYVIiIiEgConsGbGxs0Lx581r35+XlVXsgAxERkTkRJGzmQPKbFYuKihAfHy91WCIiogbTGCsQNibRwwSJiYk692dnZ+udDBERkSkwl7F+qYguBsLCwiCTyXROEuQ6A0REROZD9DCBQqFAQkIC1Gp1jS0jI8MYeRIRETUYS7ubQHQx4Ovri/T09Fr319VrQEREZOosbQKh6GGCqKgolJWV1brf09MTKSkpBiVlau7lHpY0noN7gKTxTJ3U3z9LIvV7xdR/FlJer6lfK5EpEV0MBATo/rA6OjoiMDBQ74SIiIgam7ncBSAV0cUAERHRo85cxvqlIv1DkYmIiMissGeAiIioCnOZ+CcVvYqBW7du4fTp0/D29kbLli1RWFiIzZs3Q6VS4cUXX0S3bt2kzpOIiKjBqC2sHBBdDBw/fhxDhw5FSUkJWrRogeTkZLz44ouwsbGBWq3GihUrcOTIETz55JPGyJeIiIgkJnrOwMKFC/Hiiy+iuLgYCxYsQFhYGIYMGYLz58/j4sWLGDduHN555x1j5EpERNQguOhQHdLT0xEZGYlmzZph9uzZyM3NxbRp0zT7Z82ahRMnTuiMoVKpUFJSotVUKpX47ImIiIzA0hYdEl0MlJeXw8HBAQBga2uLJk2awNnZWbPf2dkZt27d0hlDqVTCyclJq8WuiRObChERkVFYWs+A6DkD7dq1Q3Z2Njp06AAA2LFjBxQKhWZ/Xl6eVnFQk+joaERGRmptsyq9LjYVIiIikoDoYmDcuHG4ceOG5uuRI0dq7U9MTES/fv10xpDL5ZDL5VrbKsoLxaZCRERkFFyBsA4xMTE69y9cuBDW1tZ6J0RERNTYLO3WQslXILx16xZmzJghdVgiIiIyEsmLgaKiIsTHx0sdloiIqMFY2t0EoocJEhMTde7Pzs7WOxkiIiJTYC53AUhFdDEQFhYGmUwGQai93pHJLGzmBRERkRkTPUygUCiQkJAAtVpdY8vIyDBGnkRERA1GDUGyZg5E9wz4+voiPT0doaGhNe6vq9eAgHu5hyWN5+AeIGk8qfOzJKb+vTP194qU8aS+VqmZ+ntFalL/PB6UG3dtGkv7LSa6GIiKikJZWVmt+z09PZGSkmJQUkRERNRwRBcDAQG6qztHR0cEBgbqnRAREVFj4wRCIiIiC2cuY/1SYTFARERUhWWVAhIuOtSpUydcuHBBqnBERETUQET3DHz88cc1br9y5Qo+++wzuLm5AQDefPNNwzIjIiJqJJwzUIeIiAi0adMGNjbap6rVanz++eewtbWFTCZjMUBERGZLsLCBAtHFwGuvvYZff/0V27ZtQ7du3TTbbW1t8eOPP6J79+6SJkhERETGJXrOQFxcHJYsWYLg4GCsXbtWrxdVqVQoKSnRaiqVSq9YREREUlNL2MyBXhMIR48ejbS0NOzatQvDhw9Hfn6+qPOVSiWcnJy0WuyaOH1SISIikpylLUes990Ebdq0wYEDBzBo0CD07t1b1BLE0dHRKC4u1mrzZk/XNxUiIiIygEHrDMhkMkRHR2Po0KE4cuQIFApFvc6Ty+WQy+Va2yrKCw1JhYiISDLm8fe8dCRZZ8DX1xezZ8/GY489hqtXryI8PFyKsERERI2CwwQGKioqQnx8vNRhiYiIyEhEDxMkJibq3J+dna13MkRERKbAXO4CkIroYiAsLAwymUznhEGZTGZQUkRERI3J0hYdEj1MoFAokJCQALVaXWPLyMgwRp5EREQNhusM1MHX1xfp6em17q+r14CIiIhMi+hhgqioKJSVldW639PTEykpKaITcXAPEH2OLvdyD0saT0qWdK1kGFN/r0idn5T4uSBDWNowgehiICBA94ff0dERgYGBeidERETU2Myle18qkt9aSERERObFoBUIiYiIHkVqC5v7xmKAiIioCssqBSQoBgRBQGpqKi5evAiFQoHg4GDY2tpKkRsRERE1ANHFwIgRI7B9+3Y4OTmhqKgII0aMwPHjx+Hs7Ixbt27h8ccfx88//4zWrVsbI18iIiKjM5dnCkhF9ATCpKQkqFQqAMCiRYtQWlqKP//8Ezdu3MDly5fh6OiIJUuWSJ4oERFRQxEk/M8cGDRM8NNPP+H9999Hx44dAQBt27ZFbGwspk2bpvM8lUqlKSgeEgSByxgTERE1Ar1uLXz4S/uvv/5C586dtfZ5enoiNzdX5/lKpRJOTk5aTVCX6pMKERGR5LgccT1MmTIFzz//PCoqKpCTk6O1Lz8/Hy1atNB5fnR0NIqLi7WazKqZPqkQERFJTg1BsmYORA8TTJ48WfP/oaGhuHv3rtb+//73v/Dx8dEZQy6XQy6Xa23jEAEREZkKcxnrl4roYuCzzz7TuT8mJgbW1tZ6J0REREQNS/LliIuKivD6669LHZaIiKjBcM6AgYqKihAfHy91WCIiogYjCIJkzRyIHiZITEzUuT87O1vvZIiIiCzdunXr8MEHHyA/Px/e3t745JNP0K9fvxqP3bRpEz7//HOcOXMGAODr64vly5fXenxtRBcDYWFhkMlkOqsdTgYkIiJz1lh3AezcuRORkZGIi4uDn58fVq9ejeDgYGRlZcHFxaXa8ampqRg/fjz69+8Pe3t7xMbGYujQofj999/Rpk2ber+u6GEChUKBhIQEqNXqGltGRobYkERERCZFyjkDKpUKJSUlWq3qwnsPrVq1CtOmTcPUqVPRvXt3xMXFoUmTJtiyZUuNx3/11Vd4/fXX4ePjAy8vL3z66adQq9U4ePCgqOsV3TPg6+uL9PR0hIaG1ri/rl6D2tzLPSz6HHNlSddKhjH194qp5yclB/cASeOZ+vfO0q7XmJRKJZYtW6a1LSYmBkuXLtXaVl5ejvT0dERHR2u2WVlZISgoCGlpafV6rbt376KiogItW7YUlaPoYiAqKgplZWW17vf09ERKSorYsERERCZDynUGoqOjERkZqbWt6lo7AFBYWIjKykq4urpqbXd1dcW5c+fq9Vrz5s2Du7s7goKCROUouhgICNBdLTo6OiIwMFBsWCIiIpMh5ZyBmhbaM4YVK1Zgx44dSE1Nhb29vahzDXpQEREREUnD2dkZ1tbWKCgo0NpeUFAANzc3ned++OGHWLFiBQ4cOIBevXqJfm3J1xkgIiIyd42xzoCdnR18fX21Jv89nAzo7+9f63nvv/8+3nnnHSQlJaFPnz56Xa/oYuDatWsoLCzUfH348GFMnDgRAQEBeOmll+o9yYGIiMhUNdYKhJGRkdi0aRPi4+Nx9uxZzJgxA2VlZZg6dSoAYNKkSVoTDGNjY7F48WJs2bIFHTp0QH5+PvLz83Hnzh1Rryu6GPjXv/6FY8eOAQC+++47PP3007hz5w4GDBiAu3fvIjAwEHv27BEbloiIyGQIEv4nxtixY/Hhhx9iyZIl8PHxQWZmJpKSkjSTCq9cuYK8vDzN8Rs2bEB5eTleeOEFKBQKTfvwww9Fva5MEHkfYNOmTfHbb7+hY8eOeOqppzB69GjMmzdPs3/t2rXYsmWL6PUGKgq5ciERmS5Lu9XO1K/X1rmTpPGqGtpumGSxfryaJFksYxHdM2BjY4PS0lIAQE5ODoYPH661f/jw4cjKypImOyIiokaghiBZMweii4HAwEBs374dANC7d2+kpqZq7U9JSalzCUQxqzERERE1ND6oqA4rVqxAQEAAcnNzMXDgQCxcuBAnTpxAt27dkJWVhZ07dyIuLk5njJpWY1oU9SaWzJ0tNh0iIiIykOg5AwDw559/YtGiRdi7d69mxqKNjQ369u2LqKgohIWF6TxfpVJV6wmwKr3eIIsyEBHpw9TH0KVm6tdr7DkDg9s+K1mslGvJksUyFr0WHercuTO2b98OQRBw48YNqNVqODs7w9bWtl7n17QaU0V5YS1HExERNSwplyM2BwYtOiSTyeDq6gqFQqEpBK5evYrw8HBJkiMiIiLjk3wFwqKiIsTHx0sdloiIqMGoBUGyZg5EDxMkJibq3J+dzfUCiIjIvJnHr3DpiC4GwsLCIJPJdN4uIZPJDEqKiIiIGo7oYQKFQoGEhASo1eoam9iVB4mIiEwNFx2qg6+vL9LT02vdX1evARERkamztGJA9DBBVFQUysrKat3v6emJlJQUg5IyNaZ+vy09uvjeMx1Sf+9M/Wdr6tf7oPy6pPGqsrQ/akUXAwEBun+gjo6OCAwM1DshIiIialh6LTpERET0KDOX7n2psBggIiKqgisQEhERkUURXQysXLkSly9fNkYuREREJsHSHmEsuhiIiopC586d8eyzz2Lnzp0oLy83Rl5ERESNxtJuLdRrmODTTz+Fo6MjXn75Zbi7uyMiIgJnzpyROjciIiJqAHoVAyNGjMDu3btx7do1zJ07F/v374e3tzf69euHTZs2obS0VOf5KpUKJSUlWk2lUul1AURERFLjMIEILi4umDt3Ls6ePYvU1FR0794dc+bMgUKh0HmeUqmEk5OTVotdE2dIKkRERJKxtGEC0bcW1vYQooCAAAQEBODjjz/Gzp07dcaIjo5GZGSk1jarUuOuJkVEREQ1E10M1NXl0bx5c0ybNk3nMXK5HHK5XGtbRXmh2FSIiIiMwtLWGRBdDKjVamPkQUREZDLUZjLWLxXJFx26evUqwsPDpQ5LRETUYAQJ/zMHkhcDRUVFiI+PlzosERERGYnoYYLExESd+7Ozs/VOhoiIyBRY2jCB6GIgLCwMMplM50TC2u44ICIiMgfm0r0vFdHDBAqFAgkJCVCr1TW2jIwMY+RJRERERiK6Z8DX1xfp6ekIDQ2tcX9dvQa1cXAPEH2OLvdyD5tkLHNgyj8LqUl9rVIz5e8dGUbqn62pf27N7b3MYYI6REVFoaysrNb9np6eSElJMSgpIiKixmRpwwSii4GAAN3Vp6OjIwIDA/VOiIiIiBqW6GKAiIjoUcdhAiIiIgtnacMEei06tGfPHixZsgS//PILAOCnn37CiBEjMGzYMGzcuFHSBImIiMi4RBcD//nPfzB69Gjs27cPI0aMwJdffomwsDC0adMGHTp0QEREBNasWWOMXImIiBqEIKgla+ZA9DDBxx9/jPXr12PatGlISUnBiBEjsHLlSrz++usAgKeeegrvv/8+Zs+eLXmyREREDUHNYQLdcnJyEBwcDAAYPHgwKisrMWjQIM3+p59+GpcvX5YuQyIiogYmCIJkzRyILgZatWql+WWfm5uLBw8e4MqVK5r9ly9fRsuWLXXGUKlUKCkp0Wrm8g0jIiJ61IgeJggNDcUrr7yCyZMnIzExEZMmTcK///1vWFlZQSaTISoqCkOHDtUZQ6lUYtmyZVrbZFZNIbNuLjYdIiIiyXGYoA6xsbF4+umnsWPHDvj4+GDjxo145ZVXEBoaiuHDh6NVq1ZQKpU6Y0RHR6O4uFiryaya6X0RREREUrK0YQKZIFGm9+/fR0VFBZo10++Xuo1dGynS0DC3dbBNiamvcS4lPpuAHhWW9LkFAFvnTkaN3+axJySLdf2v3yWLZSx6rTNQE3t7ezRr1gxXr15FeHi4VGGJiIganFoQJGvmQLJi4KGioiLEx8dLHZaIiKjBCBL+Zw5ETyBMTEzUuT87O1vvZIiIiKjhiS4GwsLCIJPJdE6KkMlkBiVFRETUmMxl4p9URA8TKBQKJCQkQK1W19gyMjKMkScREVGDUUOQrJkD0cWAr68v0tPTa91fV68BERERmRbRwwRRUVEoKyurdb+npydSUlJEJ2Lqt7FYEkv6WVjStQKWd/sZmQ6p33sPyq9LGq8qS/ujVnQxEBCg+wfq6OiIwMBAvRMiIiJqbOZyS6BURBcDREREjzpL6xmQfJ0BIiIiMi/sGSAiIqrCXO4CkIpexcC9e/ewfft2HDlyBHl5ebCyskKnTp0QFhaGIUOGSJ0jERFRg+IwQR0uXryIbt26ITo6GgcOHMD+/fshk8lw4sQJBAcHY8yYMXjw4IExciUiIiIjEF0MvPnmmxg2bBjy8/Nx5coVKJVKqNVqHDt2DGfPnsWJEyfw7rvvGiNXIiKiBmFpDyoS/QhjR0dHZGZmokuXLgCA8vJyNG3aFHl5eWjVqhW+++47REREICcnp9YYKpUKKpVKa5tV6XXI5XI9LoGI6ovrDDy6TP1na27rDDg26SBZrLK7lySLZSyiewZatGiB0tJSzdd3797FgwcPYGdnBwDo1asX8vLydMZQKpVwcnLSarFr4sSmQkRERBIQPYHw2WefRWRkJOLi4iCXyxEdHQ0fHx80a9YMAHDlyhW4uLjojBEdHY3IyEitbValxq3yiIiI6stcuvelIroYeP/99xEaGoru3btDJpOhXbt22LVrl2b/zZs3ERUVpTOGXC6vNiRQUV4oNhUiIiKjsLS7CUQXAy4uLkhLS8OFCxegUqng5eUFG5v/C/PCCy9ImiAREREZl94rEHbp0gU9evTQKgQA4OrVqwgPDzc4MSIiosYiSPifOZB8OeKioiLEx8dLHZaIiKjBCIIgWRNr3bp16NChA+zt7eHn54fjx4/rPP6bb76Bl5cX7O3t0bNnT+zbt0/0a4oeJkhMTNS5Pzs7W3QSREREpqSx5gzs3LlTM0nfz88Pq1evRnBwMLKysmqcnH/06FGMHz8eSqUSzz33HLZt24awsDBkZGSgR48e9X5d0esMWFlZQSaT6fxGyWQyVFZWigmLikIWEUTGZur3opP+TP1na27rDNjatZEsVoWIXP38/NC3b1+sXbsWAKBWq9GuXTu88cYbmD9/frXjx44di7KyMuzZs0ez7amnnoKPjw/i4up/y77oYQKFQoGEhASo1eoaW0ZGhtiQREREJkWQsKlUKpSUlGi1qgvvAX8v4peeno6goCDNNisrKwQFBSEtLa3GPNPS0rSOB4Dg4OBaj6/9gkUKCQkRFi9eXOv+zMxMQSaTiQ1bL/fv3xdiYmKE+/fvm1Qsxnu045lyboxnOrEYz/TimYqYmJhqNUJMTEy1465fvy4AEI4ePaq1PSoqSujXr1+NsW1tbYVt27ZpbVu3bp3g4uIiKkfRxcDPP/8s/PDDD7Xuv3PnjpCamio2bL0UFxcLAITi4mKTisV4j3Y8U86N8UwnFuOZXjxTcf/+faG4uFir1VTwNGYxIHoCYUCA7nEfR0dHBAYGig1LRET0SKppob2aODs7w9raGgUFBVrbCwoK4ObmVuM5bm5uoo6vjeS3FhIREZF4dnZ28PX1xcGDBzXb1Go1Dh48CH9//xrP8ff31zoeAJKTk2s9vjaiewaIiIjIOCIjIzF58mT06dMH/fr1w+rVq1FWVoapU6cCACZNmoQ2bdpAqVQCAGbPno3AwECsXLkSI0eOxI4dO3Dy5Els3LhR1OuaVTEgl8sRExMjyaOOpYzFeI92PFPOjfFMJxbjmV48czR27FjcvHkTS5YsQX5+Pnx8fJCUlARXV1cAfz8M0Mrq/zr1+/fvj23btmHRokVYsGABunTpgt27d4taYwDQY50BIiIierRwzgAREZGFYzFARERk4VgMEBERWTgWA0RERBaOxQAREZGFM5tiQOzznXX5+eefERISAnd3d8hkMuzevVvvWEqlEn379kWzZs3g4uKCsLAwZGVl6R1vw4YN6NWrF5o3b47mzZvD398fP/zwg97x/mnFihWQyWSIiIjQO8bSpUshk8m0mpeXl97xrl+/jpdeegmtWrWCg4MDevbsiZMnT+oVq0OHDtVyk8lkmDlzpl7xKisrsXjxYnTs2BEODg7o3Lkz3nnnHYMebVpaWoqIiAh4eHjAwcEB/fv3x4kTJ+p1bl3vW0EQsGTJEigUCjg4OCAoKAgXLlzQO15CQgKGDh2KVq1aQSaTITMzU69YFRUVmDdvHnr27AlHR0e4u7tj0qRJyM3N1Tu3pUuXwsvLC46OjnjssccQFBSEX3/9Ve94/zR9+nTIZDKsXr1a73hTpkyp9j4cNmyYQfmdPXsWo0aNgpOTExwdHdG3b19cuXJFr3g1fU5kMhk++OAD0bHu3LmDWbNmoW3btnBwcED37t11Pi2vrngFBQWYMmUK3N3d0aRJEwwbNkzn+5ikYRbFwMPnO8fExCAjIwPe3t4IDg7GjRs39IpXVlYGb29vrFu3zuDcDh06hJkzZ+LYsWNITk5GRUUFhg4dirKyMr3itW3bFitWrEB6ejpOnjyJZ555BqGhofj9998NyvPEiRP4z3/+g169ehkUBwCeeOIJ5OXladqRI0f0ivPXX39hwIABsLW1xQ8//IA//vgDK1euxGOPPaZXvBMnTmjllZycDAB48cUX9YoXGxuLDRs2YO3atTh79ixiY2Px/vvv45NPPtErHgC8+uqrSE5OxhdffIHffvsNQ4cORVBQEK5fr/sRp3W9b99//318/PHHiIuLw6+//gpHR0cEBwfj/v37esUrKyvDwIEDERsba1Bud+/eRUZGBhYvXoyMjAwkJCQgKysLo0aN0vtaH3/8caxduxa//fYbjhw5gg4dOmDo0KG4efOmXvEe2rVrF44dOwZ3d3edx9Un3rBhw7Tej9u3b9c73p9//omBAwfCy8sLqampOH36NBYvXgx7e3u94v0zr7y8PGzZsgUymQz/+te/RMeKjIxEUlISvvzyS5w9exYRERGYNWsWEhMTRecmCALCwsKQnZ2N7777DqdOnYKHhweCgoL0/jeV6knUkwwaSb9+/YSZM2dqvq6srBTc3d0FpVJpcGwAwq5duwyO89CNGzcEAMKhQ4cki/nYY48Jn376qd7nl5aWCl26dBGSk5OFwMBAYfbs2XrHiomJEby9vfU+/5/mzZsnDBw4UJJYNZk9e7bQuXNnQa1W63X+yJEjhfDwcK1tzz//vDBx4kS94t29e1ewtrYW9uzZo7X9ySefFBYuXCgqVtX3rVqtFtzc3IQPPvhAs+327duCXC4Xtm/fLjreP+Xk5AgAhFOnTumVW02OHz8uABAuX74sSbyHD7g5cOCA3vGuXbsmtGnTRjhz5ozg4eEhfPTRR3XGqi3e5MmThdDQ0HqdX594Y8eOFV566SXJ4lUVGhoqPPPMM3rFeuKJJ4S3335ba1t939NV42VlZQkAhDNnzmi2VVZWCq1btxY2bdpUZzzSn8n3DOjzfOfGVFxcDABo2bKlwbEqKyuxY8cOlJWViV5n+p9mzpyJkSNHVnvmtb4uXLgAd3d3dOrUCRMnTqy1q7IuiYmJ6NOnD1588UW4uLigd+/e2LRpkyQ5lpeX48svv0R4eDhkMpleMfr374+DBw/i/PnzAID//e9/OHLkCIYPH65XvAcPHqCysrLaX3MODg569648lJOTg/z8fK2fsZOTE/z8/Ez2cyKTydCiRQuDY5WXl2Pjxo1wcnKCt7e3XjHUajVefvllREVF4YknnjA4JwBITU2Fi4sLunbtihkzZuDWrVt657Z37148/vjjCA4OhouLC/z8/Awa3vyngoIC7N27F6+88ope5/fv3x+JiYm4fv06BEFASkoKzp8/j6FDh4qOpVKpAEDrM2JlZQW5XG7wZ4R0M/lioLCwEJWVlZqlGB9ydXVFfn5+I2VVM7VajYiICAwYMED0UpD/9Ntvv6Fp06aQy+WYPn06du3ahe7du+sVa8eOHcjIyNCsY20oPz8/bN26FUlJSdiwYQNycnIQEBCA0tJS0bGys7OxYcMGdOnSBfv378eMGTPw5ptvIj4+3uA8d+/ejdu3b2PKlCl6x5g/fz7GjRsHLy8v2Nraonfv3oiIiMDEiRP1itesWTP4+/vjnXfeQW5uLiorK/Hll18iLS0NeXl5eucJQPNZMIfPyf379zFv3jyMHz8ezZs31zvOnj170LRpU9jb2+Ojjz5CcnIynJ2d9YoVGxsLGxsbvPnmm3rn80/Dhg3D559/joMHDyI2NhaHDh3C8OHDUVlZKTrWjRs3cOfOHaxYsQLDhg3Djz/+iNGjR+P555/HoUOHDM41Pj4ezZo1w/PPP6/X+Z988gm6d++Otm3bws7ODsOGDcO6deswaNAg0bG8vLzQvn17REdH46+//kJ5eTliY2Nx7do1gz8jpJtZPZvA1M2cORNnzpwxuILt2rUrMjMzUVxcjG+//RaTJ0/GoUOHRBcEV69exezZs5GcnFzr2KJY//yruFevXvDz84OHhwe+/vpr0X9ZqNVq9OnTB8uXLwcA9O7dG2fOnEFcXBwmT55sUJ6bN2/G8OHD6xz71eXrr7/GV199hW3btuGJJ55AZmYmIiIi4O7urnd+X3zxBcLDw9GmTRtYW1vjySefxPjx45Genq53nuakoqICY8aMgSAI2LBhg0GxBg8ejMzMTBQWFmLTpk0YM2YMfv31V7i4uIiKk56ejjVr1iAjI0PvXqSqxo0bp/n/nj17olevXujcuTNSU1MxZMgQUbHUajUAIDQ0FHPmzAEA+Pj44OjRo4iLizP4kfFbtmzBxIkT9f434pNPPsGxY8eQmJgIDw8P/Pzzz5g5cybc3d1F90ba2toiISEBr7zyClq2bAlra2sEBQVh+PDhBk3cpbqZfM+APs93bgyzZs3Cnj17kJKSgrZt2xoUy87ODp6envD19YVSqYS3tzfWrFkjOk56ejpu3LiBJ598EjY2NrCxscGhQ4fw8ccfw8bGRq+/Uqpq0aIFHn/8cVy8eFH0uQqFolqB061bN72HHR66fPkyDhw4gFdffdWgOFFRUZregZ49e+Lll1/GnDlzDOpl6dy5Mw4dOoQ7d+7g6tWrOH78OCoqKtCpUyeDcn34WTDlz8nDQuDy5ctITk42qFcAABwdHeHp6YmnnnoKmzdvho2NDTZv3iw6zuHDh3Hjxg20b99e8zm5fPky/v3vf6NDhw4G5fhQp06d4OzsrNfnxNnZGTY2Nkb5rBw+fBhZWVl6f1bu3buHBQsWYNWqVQgJCUGvXr0wa9YsjB07Fh9++KFeMX19fZGZmYnbt28jLy8PSUlJuHXrlsGfEdLN5IsBfZ7v3JAEQcCsWbOwa9cu/PTTT+jYsaPkr6FWqzVjaWIMGTIEv/32GzIzMzWtT58+mDhxIjIzM2FtbW1wbnfu3MGff/4JhUIh+twBAwZUuw3z/Pnz8PDwMCinzz77DC4uLhg5cqRBce7evav1dDAAsLa21vylZghHR0coFAr89ddf2L9/P0JDQw2K17FjR7i5uWl9TkpKSvDrr7+axOfkYSFw4cIFHDhwAK1atZL8NfT9nLz88ss4ffq01ufE3d0dUVFR2L9/vyS5Xbt2Dbdu3dLrc2JnZ4e+ffsa5bOyefNm+Pr66j3XoqKiAhUVFUb5nDg5OaF169a4cOECTp48afBnhHQzi2GCup7vLNadO3e0KvScnBxkZmaiZcuWaN++vahYM2fOxLZt2/Ddd9+hWbNmmvFZJycnODg4iM4tOjoaw4cPR/v27VFaWopt27YhNTVVr3+UmjVrVm3ugqOjI1q1aqX3nIa33noLISEh8PDwQG5uLmJiYmBtbY3x48eLjjVnzhz0798fy5cvx5gxY3D8+HFs3LhR9HO4/0mtVuOzzz7D5MmTYWNj2Ns7JCQE7733Htq3b48nnngCp06dwqpVqxAeHq53zP3790MQBHTt2hUXL15EVFQUvLy86vVerut9GxERgXfffRddunRBx44dsXjxYri7uyMsLEyveEVFRbhy5YpmPYCHv4zc3Nyq9TboiqVQKPDCCy8gIyMDe/bsQWVlpeZz0rJlS9jZ2YnKrVWrVnjvvfcwatQoKBQKFBYWYt26dbh+/Xqtt5HWda1VixNbW1u4ubmha9euouO1bNkSy5Ytw7/+9S+4ubnhzz//xNy5c+Hp6Yng4GC98ouKisLYsWMxaNAgDB48GElJSfj++++RmpqqVzzg72Lxm2++wcqVK2uMUd9YgYGBiIqKgoODAzw8PHDo0CF8/vnnWLVqlV7xvvnmG7Ru3Rrt27fHb7/9htmzZyMsLEyvCYkkQqPeyyDCJ598IrRv316ws7MT+vXrJxw7dkzvWCkpKQKAam3y5MmiY9UUB4Dw2Wef6ZVbeHi44OHhIdjZ2QmtW7cWhgwZIvz44496xaqJobcWjh07VlAoFIKdnZ3Qpk0bYezYscLFixf1jvf9998LPXr0EORyueDl5SVs3LhR71iCIAj79+8XAAhZWVkGxREEQSgpKRFmz54ttG/fXrC3txc6deokLFy4UFCpVHrH3Llzp9CpUyfBzs5OcHNzE2bOnCncvn27XufW9b5Vq9XC4sWLBVdXV0EulwtDhgzR+X2oK95nn31W4/6YmBhRsR7emlhTS0lJEZ3bvXv3hNGjRwvu7u6CnZ2doFAohFGjRgnHjx/X+1qrquvWQl3x7t69KwwdOlRo3bq1YGtrK3h4eAjTpk0T8vPzDcpv8+bNgqenp2Bvby94e3sLu3fvNijef/7zH8HBwaHO919dsfLy8oQpU6YI7u7ugr29vdC1a1dh5cqVtd7SW1e8NWvWCG3bthVsbW2F9u3bC4sWLTLoM0f1IxMEzsogIiKyZCY/Z4CIiIiMi8UAERGRhWMxQEREZOFYDBAREVk4FgNEREQWjsUAERGRhWMxQEREZOFYDBAREVk4FgNEREQWjsUAERGRhWMxQEREZOH+P4WljwVR+uyVAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import seaborn as sns\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "target_figure = 1\n",
    "sns.heatmap(np.reshape(X[:,target_figure], (20,20)))\n",
    "plt.title(\"Original\")\n",
    "plt.show()\n",
    "\n",
    "sns.heatmap(np.reshape(Xtrain[:,target_figure], (20,20)))\n",
    "plt.title(\"Training\")\n",
    "plt.show()\n",
    "\n",
    "sns.heatmap(np.reshape(MASK[:,target_figure], (20,20)))\n",
    "plt.title(\"Mask\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Factorize"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "from TELF.factorization.decompositions.lmf import LogisticMatrixFactorization\n",
    "\n",
    "model = LogisticMatrixFactorization(\n",
    "    k=5,\n",
    "    l2_p=1e-6, \n",
    "    epochs=3000, \n",
    "    learning_rate=0.001, \n",
    "    tolerance=1e-3, \n",
    "    device=2,\n",
    "    random_state=1\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "100%|██████████| 3000/3000 [00:01<00:00, 1522.40it/s]\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAHHCAYAAABeLEexAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABSC0lEQVR4nO3deVxU9f4/8NfMwAzrsMOAoiIYiIrmhnPNpSTRqKtmtzS3zPJrYb/UMvNWpnbvtfS2eFv09u0WVpotX63UXBBFS3HjigsquaCYOuDGDCDrzOf3B8zREVTUgTMwr+fjcR7OOeczZ97nCM7L8/mccxRCCAEiIiIiJ6aUuwAiIiIiuTEQERERkdNjICIiIiKnx0BERERETo+BiIiIiJweAxERERE5PQYiIiIicnoMREREROT0GIiIiIjI6TEQEZFDeuqpp9CmTZs7eu/s2bOhUCjsWxARNWsMRER0WxQKRb2m9PR0uUuVxVNPPQUvLy+5yyCi26Tgs8yI6HZ8/fXXNvNffvklUlNT8dVXX9ksf/DBBxESEnLHn1NZWQmLxQKNRnPb762qqkJVVRXc3Nzu+PPv1FNPPYUffvgBxcXFjf7ZRHTnXOQugIialtGjR9vM79ixA6mpqbWWX+/KlSvw8PCo9+e4urreUX0A4OLiAhcX/vNGRPXHLjMisrv+/fujY8eOyMzMRN++feHh4YG//vWvAICffvoJSUlJCAsLg0ajQWRkJN566y2YzWabbVw/hujkyZNQKBT45z//iU8//RSRkZHQaDTo0aMHdu/ebfPeusYQKRQKTJ48GT/++CM6duwIjUaDDh06YN26dbXqT09PR/fu3eHm5obIyEj8+9//tvu4pO+//x7dunWDu7s7AgMDMXr0aJw5c8amjcFgwPjx49GyZUtoNBqEhoZiyJAhOHnypNRmz549SExMRGBgINzd3REREYGnn37abnUSOQv+F4qIGsTFixcxePBgjBgxAqNHj5a6z1JSUuDl5YVp06bBy8sLmzZtwqxZs2AymbBgwYJbbnfZsmUoKirC//zP/0ChUGD+/Pl49NFHceLEiVueVfrtt9+wYsUKPP/88/D29sa//vUvDB8+HHl5eQgICAAA7N27F4MGDUJoaCjmzJkDs9mMuXPnIigo6O4PSo2UlBSMHz8ePXr0wLx585Cfn4+FCxdi27Zt2Lt3L3x9fQEAw4cPR3Z2Nl544QW0adMGBQUFSE1NRV5enjQ/cOBABAUF4dVXX4Wvry9OnjyJFStW2K1WIqchiIjuQnJysrj+n5J+/foJAGLx4sW12l+5cqXWsv/5n/8RHh4eoqysTFo2btw40bp1a2k+NzdXABABAQHi0qVL0vKffvpJABCrVq2Slr355pu1agIg1Gq1OHbsmLRs3759AoD48MMPpWWPPPKI8PDwEGfOnJGWHT16VLi4uNTaZl3GjRsnPD09b7i+oqJCBAcHi44dO4rS0lJp+erVqwUAMWvWLCGEEJcvXxYAxIIFC264rZUrVwoAYvfu3besi4hujl1mRNQgNBoNxo8fX2u5u7u79LqoqAgXLlxAnz59cOXKFRw5cuSW233iiSfg5+cnzffp0wcAcOLEiVu+NyEhAZGRkdJ8XFwctFqt9F6z2YyNGzdi6NChCAsLk9pFRUVh8ODBt9x+fezZswcFBQV4/vnnbQZ9JyUlISYmBmvWrAFQfZzUajXS09Nx+fLlOrdlPZO0evVqVFZW2qU+ImfFQEREDaJFixZQq9W1lmdnZ2PYsGHw8fGBVqtFUFCQNCDbaDTecrutWrWymbeGoxuFhpu91/p+63sLCgpQWlqKqKioWu3qWnYnTp06BQCIjo6utS4mJkZar9Fo8M4772Dt2rUICQlB3759MX/+fBgMBql9v379MHz4cMyZMweBgYEYMmQIvvjiC5SXl9ulViJnwkBERA3i2jNBVoWFhejXrx/27duHuXPnYtWqVUhNTcU777wDALBYLLfcrkqlqnO5qMcdRO7mvXKYMmUKfv/9d8ybNw9ubm5444030L59e+zduxdA9UDxH374ARkZGZg8eTLOnDmDp59+Gt26deNl/0S3iYGIiBpNeno6Ll68iJSUFLz44ot4+OGHkZCQYNMFJqfg4GC4ubnh2LFjtdbVtexOtG7dGgCQk5NTa11OTo603ioyMhIvvfQSNmzYgIMHD6KiogLvvvuuTZtevXrh73//O/bs2YOlS5ciOzsby5cvt0u9RM6CgYiIGo31DM21Z2QqKirwySefyFWSDZVKhYSEBPz44484e/astPzYsWNYu3atXT6je/fuCA4OxuLFi226ttauXYvDhw8jKSkJQPV9m8rKymzeGxkZCW9vb+l9ly9frnV2q0uXLgDAbjOi28TL7omo0fzpT3+Cn58fxo0bh//3//4fFAoFvvrqK4fqspo9ezY2bNiA3r1747nnnoPZbMZHH32Ejh07Iisrq17bqKysxN/+9rday/39/fH888/jnXfewfjx49GvXz+MHDlSuuy+TZs2mDp1KgDg999/x4ABA/D4448jNjYWLi4uWLlyJfLz8zFixAgAwJIlS/DJJ59g2LBhiIyMRFFREf73f/8XWq0WDz30kN2OCZEzYCAiokYTEBCA1atX46WXXsLrr78OPz8/jB49GgMGDEBiYqLc5QEAunXrhrVr1+Lll1/GG2+8gfDwcMydOxeHDx+u11VwQPVZrzfeeKPW8sjISDz//PN46qmn4OHhgbfffhszZsyAp6cnhg0bhnfeeUe6ciw8PBwjR45EWloavvrqK7i4uCAmJgbfffcdhg8fDqB6UPWuXbuwfPly5Ofnw8fHBz179sTSpUsRERFht2NC5Az4LDMionoYOnQosrOzcfToUblLIaIGwDFERETXKS0ttZk/evQofvnlF/Tv31+egoiowfEMERHRdUJDQ/HUU0+hbdu2OHXqFBYtWoTy8nLs3bsX7dq1k7s8ImoAHENERHSdQYMG4ZtvvoHBYIBGo4Fer8c//vEPhiGiZoxniIiIiMjpyTqGaNGiRdKzhLRaLfR6vc29Pvr37w+FQmEzTZo0yWYbeXl5SEpKgoeHB4KDgzF9+nRUVVXZtElPT0fXrl2h0WgQFRWFlJSUxtg9IiIiaiJk7TJr2bIl3n77bbRr1w5CCCxZsgRDhgzB3r170aFDBwDAs88+i7lz50rv8fDwkF6bzWYkJSVBp9Nh+/btOHfuHMaOHQtXV1f84x//AADk5uYiKSkJkyZNwtKlS5GWloZnnnkGoaGhDnOZLxEREcnL4brM/P39sWDBAkyYMAH9+/dHly5d8MEHH9TZdu3atXj44Ydx9uxZhISEAAAWL16MGTNm4Pz581Cr1ZgxYwbWrFmDgwcPSu8bMWIECgsLsW7dunrVZLFYcPbsWXh7e0OhUNz1PhIREVHDE0KgqKgIYWFhUCpv0SkmHERVVZX45ptvhFqtFtnZ2UIIIfr16ycCAwNFQECA6NChg3j11VdFSUmJ9J433nhDdO7c2WY7J06cEADEf//7XyGEEH369BEvvviiTZvPP/9caLXaG9ZSVlYmjEajNB06dEgA4MSJEydOnDg1wen06dO3zCGyX2V24MAB6PV6lJWVwcvLCytXrkRsbCwA4Mknn0Tr1q0RFhaG/fv3Y8aMGcjJycGKFSsAAAaDQTozZGWdNxgMN21jMplQWlpa5xO5582bhzlz5tRafvr0aWi12rvfaSIiImpwJpMJ4eHh8Pb2vmVb2QNRdHQ0srKyYDQa8cMPP2DcuHHYsmULYmNjMXHiRKldp06dEBoaigEDBuD48eOIjIxssJpmzpyJadOmSfPWA2od/E1ERERNR32Gu8h+p2q1Wo2oqCh069YN8+bNQ+fOnbFw4cI628bHxwOofvI0AOh0OuTn59u0sc7rdLqbttFqtXWeHQIAjUYjhR+GICIiouZP9kB0PYvFgvLy8jrXWZ80HRoaCgDQ6/U4cOAACgoKpDapqanQarVSt5ter0daWprNdlJTU6HX6xugeiIiImqKZO0ymzlzJgYPHoxWrVqhqKgIy5YtQ3p6OtavX4/jx49j2bJleOihhxAQEID9+/dj6tSp6Nu3L+Li4gAAAwcORGxsLMaMGYP58+fDYDDg9ddfR3JyMjQaDQBg0qRJ+Oijj/DKK6/g6aefxqZNm/Ddd99hzZo1cu46ERERORBZA1FBQQHGjh2Lc+fOwcfHB3FxcVi/fj0efPBBnD59Ghs3bsQHH3yAkpIShIeHY/jw4Xj99del96tUKqxevRrPPfcc9Ho9PD09MW7cOJv7FkVERGDNmjWYOnUqFi5ciJYtW+Kzzz7jPYiIiIhI4nD3IXJEJpMJPj4+MBqNHE9ERETURNzO97fDjSEiIiIiamwMREREROT0GIiIiIjI6TEQERERkdNjICIiIiKnx0BERERETo+BiIiIiJye7A93dWZVZgvyi8phsQiE+3vIXQ4REZHTYiCS0cWSCvR+exNUSgWO/+MhucshIiJyWuwyk5FSoQAAmC0CvGE4ERGRfBiIZOSiVEivLcxDREREsmEgkpFKdTUQVVksMlZCRETk3BiIZKRSXA1EZp4iIiIikg0DkYxUSgYiIiIiR8BAJCMXBiIiIiKHwEAko2vPEFUxEBEREcmGgUhGCoUC1kxkYSAiIiKSDQORzFyU1X8FPENEREQkHwYimVm7zTiGiIiISD4MRDJjICIiIpIfA5HMrIGIXWZERETyYSCSmQvPEBEREcmOgUhmSgYiIiIi2TEQyYxniIiIiOTHQCSzq2OI+HBXIiIiuTAQycwaiCyCZ4iIiIjkwkAkM+kMkZmBiIiISC4MRDLjGCIiIiL5MRDJTFXz6A4zu8yIiIhkw0AkM1XN3wBvzEhERCQfBiKZSWeIOIaIiIhINgxEMpPGELHLjIiISDYMRDJTKTiomoiISG4MRDLjw12JiIjkx0AkMxdVzY0ZGYiIiIhkw0AkM6WCZ4iIiIjkxkAks6s3ZuSzzIiIiOTCQCQzlRSIZC6EiIjIiTEQycw6hohniIiIiOQjayBatGgR4uLioNVqodVqodfrsXbtWml9WVkZkpOTERAQAC8vLwwfPhz5+fk228jLy0NSUhI8PDwQHByM6dOno6qqyqZNeno6unbtCo1Gg6ioKKSkpDTG7tULxxARERHJT9ZA1LJlS7z99tvIzMzEnj178MADD2DIkCHIzs4GAEydOhWrVq3C999/jy1btuDs2bN49NFHpfebzWYkJSWhoqIC27dvx5IlS5CSkoJZs2ZJbXJzc5GUlIT7778fWVlZmDJlCp555hmsX7++0fe3Li582j0REZHsFEI41i2S/f39sWDBAjz22GMICgrCsmXL8NhjjwEAjhw5gvbt2yMjIwO9evXC2rVr8fDDD+Ps2bMICQkBACxevBgzZszA+fPnoVarMWPGDKxZswYHDx6UPmPEiBEoLCzEunXr6lWTyWSCj48PjEYjtFqtXfd3+vf78H3mH3hlUDSe7x9l120TERE5s9v5/naYMURmsxnLly9HSUkJ9Ho9MjMzUVlZiYSEBKlNTEwMWrVqhYyMDABARkYGOnXqJIUhAEhMTITJZJLOMmVkZNhsw9rGuo26lJeXw2Qy2UwNxdWl+q+gssqhcikREZFTkT0QHThwAF5eXtBoNJg0aRJWrlyJ2NhYGAwGqNVq+Pr62rQPCQmBwWAAABgMBpswZF1vXXezNiaTCaWlpXXWNG/ePPj4+EhTeHi4PXa1Tuqax91XcVA1ERGRbGQPRNHR0cjKysLOnTvx3HPPYdy4cTh06JCsNc2cORNGo1GaTp8+3WCf5VpzlVkFr7snIiKSjYvcBajVakRFVY+d6datG3bv3o2FCxfiiSeeQEVFBQoLC23OEuXn50On0wEAdDoddu3aZbM961Vo17a5/sq0/Px8aLVauLu711mTRqOBRqOxy/7diquKXWZERERyk/0M0fUsFgvKy8vRrVs3uLq6Ii0tTVqXk5ODvLw86PV6AIBer8eBAwdQUFAgtUlNTYVWq0VsbKzU5tptWNtYtyE3KRDxDBEREZFsZD1DNHPmTAwePBitWrVCUVERli1bhvT0dKxfvx4+Pj6YMGECpk2bBn9/f2i1WrzwwgvQ6/Xo1asXAGDgwIGIjY3FmDFjMH/+fBgMBrz++utITk6WzvBMmjQJH330EV555RU8/fTT2LRpE7777jusWbNGzl2XqF0YiIiIiOQmayAqKCjA2LFjce7cOfj4+CAuLg7r16/Hgw8+CAB4//33oVQqMXz4cJSXlyMxMRGffPKJ9H6VSoXVq1fjueeeg16vh6enJ8aNG4e5c+dKbSIiIrBmzRpMnToVCxcuRMuWLfHZZ58hMTGx0fe3LhxDREREJD+Huw+RI2rI+xB9/lsu5q4+hEc6h+HDkffaddtERETOrEneh8hZXb0PEc8QERERyYWBSGbqmi4zjiEiIiKSDwORzKxXmXEMERERkXwYiGTGy+6JiIjkx0AkM2sg4tPuiYiI5MNAJDO1C8cQERERyY2BSGZXxxDxDBEREZFcGIhkxjFERERE8mMgkhkDERERkfwYiGSmtnaZ8caMREREsmEgkpnGtfqvoJyBiIiISDYMRDJzd1UBAEorzDJXQkRE5LwYiGRmPUNUVmUGn7NLREQkDwYimVnPEAnBbjMiIiK5MBDJzK0mEAFAWSW7zYiIiOTAQCQzV5USrjVPvC9lICIiIpIFA5EDcOPAaiIiIlkxEDkAayAqq+QYIiIiIjkwEDkA6dJ7dpkRERHJgoHIAbhLZ4gYiIiIiOTAQOQA3NQcQ0RERCQnBiIH4OZS/dfALjMiIiJ5MBA5AHc1u8yIiIjkxEDkADiGiIiISF4MRA6AV5kRERHJi4HIAVgHVV/hoGoiIiJZMBA5AG83FwBAUVmVzJUQERE5JwYiB6B1cwUAmEorZa6EiIjIOTEQOQAf95pAVMZAREREJAcGIgegtQaiUnaZERERyYGByAFoa8YQ8QwRERGRPBiIHID1DJGRY4iIiIhkwUDkAKQxRAxEREREsmAgcgDWq8yKyqtgsQiZqyEiInI+DEQOwHofIiGqQxERERE1LgYiB+DmqoJHzd2qL5dUyFwNERGR82EgchDB3hoAwPnicpkrISIicj4MRA4iqCYQFZgYiIiIiBobA5GDCPZ2AwAUFJXJXAkREZHzYSByENYzROeLeIaIiIiosckaiObNm4cePXrA29sbwcHBGDp0KHJycmza9O/fHwqFwmaaNGmSTZu8vDwkJSXBw8MDwcHBmD59OqqqbK/WSk9PR9euXaHRaBAVFYWUlJSG3r3bInWZMRARERE1OlkD0ZYtW5CcnIwdO3YgNTUVlZWVGDhwIEpKSmzaPfvsszh37pw0zZ8/X1pnNpuRlJSEiooKbN++HUuWLEFKSgpmzZoltcnNzUVSUhLuv/9+ZGVlYcqUKXjmmWewfv36RtvXWwnRVneZGYzsMiMiImpsLnJ++Lp162zmU1JSEBwcjMzMTPTt21da7uHhAZ1OV+c2NmzYgEOHDmHjxo0ICQlBly5d8NZbb2HGjBmYPXs21Go1Fi9ejIiICLz77rsAgPbt2+O3337D+++/j8TExIbbwdsQ7ucOAMi7dEXmSoiIiJyPQ40hMhqNAAB/f3+b5UuXLkVgYCA6duyImTNn4sqVq6EhIyMDnTp1QkhIiLQsMTERJpMJ2dnZUpuEhASbbSYmJiIjI6POOsrLy2EymWymhtY6wBMAcKawFJVmS4N/HhEREV0l6xmia1ksFkyZMgW9e/dGx44dpeVPPvkkWrdujbCwMOzfvx8zZsxATk4OVqxYAQAwGAw2YQiANG8wGG7axmQyobS0FO7u7jbr5s2bhzlz5th9H28m2FsDjYsS5VUWnC0slQISERERNTyHCUTJyck4ePAgfvvtN5vlEydOlF536tQJoaGhGDBgAI4fP47IyMgGqWXmzJmYNm2aNG8ymRAeHt4gn2WlVCrQyt8DRwuKceriFQYiIiKiRuQQXWaTJ0/G6tWrsXnzZrRs2fKmbePj4wEAx44dAwDodDrk5+fbtLHOW8cd3aiNVqutdXYIADQaDbRarc3UGFoHeAAATl0suUVLIiIisidZA5EQApMnT8bKlSuxadMmRERE3PI9WVlZAIDQ0FAAgF6vx4EDB1BQUCC1SU1NhVarRWxsrNQmLS3NZjupqanQ6/V22hP7iAr2BgAcMRTJXAkREZFzkTUQJScn4+uvv8ayZcvg7e0Ng8EAg8GA0tJSAMDx48fx1ltvITMzEydPnsTPP/+MsWPHom/fvoiLiwMADBw4ELGxsRgzZgz27duH9evX4/XXX0dycjI0mup7+0yaNAknTpzAK6+8giNHjuCTTz7Bd999h6lTp8q273WJDas+E3XoXMMP4iYiIqKrZA1EixYtgtFoRP/+/REaGipN3377LQBArVZj48aNGDhwIGJiYvDSSy9h+PDhWLVqlbQNlUqF1atXQ6VSQa/XY/To0Rg7dizmzp0rtYmIiMCaNWuQmpqKzp07491338Vnn33mMJfcW8WGVgeiI+eKYLYImashIiJyHgohBL95b8FkMsHHxwdGo7FBxxOZLQId3lyHskoL0l7qh8ggrwb7LCIioubudr6/HWJQNVVTKRWI0VX/hR08Y5S5GiIiIufBQORg7m3lCwDYc/KyvIUQERE5EQYiB9OzTfVdunflXpK5EiIiIufBQORgekRUB6Kc/CIUXqmQuRoiIiLnwEDkYAK9NGgbVH2XanabERERNQ4GIgcUX3OWaMeJizJXQkRE5BwYiBxQ76hAAMCW38/LXAkREZFzYCByQH2igqBUAEcLivHH5Styl0NERNTsMRA5IB8PV3Rt5QcASM/hWSIiIqKGxkDkoPpHBwFgICIiImoMDEQOqn90MABg+/ELKK8yy1wNERFR88ZA5KA6hGkR5K3BlQozdpzgTRqJiIgaEgORg1IoFHgwNgQAsD7bIHM1REREzRsDkQMb1EEHANiQbYDZImSuhoiIqPliIHJg+sgAaN1ccKG4ApmneNdqIiKihsJA5MBcVUok1HSbrTvIbjMiIqKGwkDk4KzdZuuzDRCC3WZEREQNgYHIwfW9JwgeahXOFJbiwBmj3OUQERE1SwxEDs7NVYX7a+5JtJbdZkRERA2CgagJGNyputts9f6z7DYjIiJqAAxETcADMcHwUKtw+lIpsk4Xyl0OERFRs8NA1AR4qF2Q0L76arNV+87JXA0REVHzw0DURPy5cxiA6m4z3qSRiIjIvhiImog+9wRC6+aCgqJy7Mrls82IiIjsiYGoidC4qDC4YygAYNX+szJXQ0RE1LwwEDUhj9R0m609cA6VZovM1RARETUfDERNSK+2/gj0UuPylUr8duyC3OUQERE1GwxETYiLSomkTjXdZvvYbUZERGQvDERNjLXbbEN2PsoqzTJXQ0RE1DwwEDUxXVv5oYWvO4rLq7D5SIHc5RARETULDERNjFKpwMNx1d1mq/fzJo1ERET2wEDUBFm7zdKO5KO4vErmaoiIiJo+BqImqEOYFhGBniirtCDtcL7c5RARETV5DERNkEKhwCNxvNqMiIjIXhiImihrt9mW38/DeKVS5mqIiIiaNgaiJqpdiDdidN6oNAuszzbIXQ4REVGTxkDUhFnPEvHZZkRERHeHgagJs15+v+3YBVwoLpe5GiIioqaLgagJax3gic4tfWAR1Q98JSIiojsjayCaN28eevToAW9vbwQHB2Po0KHIycmxaVNWVobk5GQEBATAy8sLw4cPR36+7aXmeXl5SEpKgoeHB4KDgzF9+nRUVdnenyc9PR1du3aFRqNBVFQUUlJSGnr3GoXUbbaPgYiIiOhOyRqItmzZguTkZOzYsQOpqamorKzEwIEDUVJSIrWZOnUqVq1ahe+//x5btmzB2bNn8eijj0rrzWYzkpKSUFFRge3bt2PJkiVISUnBrFmzpDa5ublISkrC/fffj6ysLEyZMgXPPPMM1q9f36j72xCSarrNdp28hHPGUpmrISIiapoUQgghdxFW58+fR3BwMLZs2YK+ffvCaDQiKCgIy5Ytw2OPPQYAOHLkCNq3b4+MjAz06tULa9euxcMPP4yzZ88iJCQEALB48WLMmDED58+fh1qtxowZM7BmzRocPHhQ+qwRI0agsLAQ69atu2VdJpMJPj4+MBqN0Gq1DbPzd+HxxRnYdfISXk9qj2f6tJW7HCIiIodwO9/fDjWGyGg0AgD8/f0BAJmZmaisrERCQoLUJiYmBq1atUJGRgYAICMjA506dZLCEAAkJibCZDIhOztbanPtNqxtrNu4Xnl5OUwmk83kyB7pzJs0EhER3Q2HCUQWiwVTpkxB79690bFjRwCAwWCAWq2Gr6+vTduQkBAYDAapzbVhyLreuu5mbUwmE0pLa3czzZs3Dz4+PtIUHh5ul31sKIM7hUKpAPb9YcSpiyW3fgMRERHZcJhAlJycjIMHD2L58uVyl4KZM2fCaDRK0+nTp+Uu6aYCvTToHRUIAFi9n4OriYiIbpdDBKLJkydj9erV2Lx5M1q2bCkt1+l0qKioQGFhoU37/Px86HQ6qc31V51Z52/VRqvVwt3dvVY9Go0GWq3WZnJ0j8RVX232cxa7zYiIiG6XrIFICIHJkydj5cqV2LRpEyIiImzWd+vWDa6urkhLS5OW5eTkIC8vD3q9HgCg1+tx4MABFBQUSG1SU1Oh1WoRGxsrtbl2G9Y21m00B4kddHBRKpCTX4TcC+w2IyIiuh2yBqLk5GR8/fXXWLZsGby9vWEwGGAwGKRxPT4+PpgwYQKmTZuGzZs3IzMzE+PHj4der0evXr0AAAMHDkRsbCzGjBmDffv2Yf369Xj99deRnJwMjUYDAJg0aRJOnDiBV155BUeOHMEnn3yC7777DlOnTpVt3+3Nx8MVvdoGAAA28NlmREREt0XWQLRo0SIYjUb0798foaGh0vTtt99Kbd5//308/PDDGD58OPr27QudTocVK1ZI61UqFVavXg2VSgW9Xo/Ro0dj7NixmDt3rtQmIiICa9asQWpqKjp37ox3330Xn332GRITExt1fxtaYofqgeN82CsREdHtcaj7EDkqR78PkZXBWIZe86q7Bnf9dQCCtW4yV0RERCSfJnsfIro7Oh83dA73BQCkHs6/eWMiIiKSMBA1M1e7zRiIiIiI6ouBqJkZGFt9q4GM4xdgKquUuRoiIqKmgYGomYkK9kJkkCcqzQKbjxTc+g1ERETEQNQcJXaoPku04RC7zYiIiOqDgagZGtA+GADw6+/nUWW2yFwNERGR42Mgaoa6hPvB18MVprIqZJ0ulLscIiIih8dA1AyplAr0bRcEANicw3FEREREt8JA1Ez1j64JREfOy1wJERGR42Mgaqb63hMEhQI4dM6EfFOZ3OUQERE5NAaiZirQS4O4lr4AgC05PEtERER0MwxEzVj/eziOiIiIqD4YiJqx+2OqL7//7egFVPLyeyIiohtiIGrG4lr4wN9TjaLyKuzNK5S7HCIiIod1R4Ho9OnT+OOPP6T5Xbt2YcqUKfj000/tVhjdPaVSgT9FBgAAth27IHM1REREjuuOAtGTTz6JzZs3AwAMBgMefPBB7Nq1C6+99hrmzp1r1wLp7twXFQiAgYiIiOhm7igQHTx4ED179gQAfPfdd+jYsSO2b9+OpUuXIiUlxZ710V3qXROIsk4Xori8SuZqiIiIHNMdBaLKykpoNBoAwMaNG/HnP/8ZABATE4Nz587Zrzq6a+H+Hgj3d0eVRWBX7kW5yyEiInJIdxSIOnTogMWLF+PXX39FamoqBg0aBAA4e/YsAgIC7Fog3b2r3WYMRERERHW5o0D0zjvv4N///jf69++PkSNHonPnzgCAn3/+WepKI8fxp0iOIyIiIroZlzt5U//+/XHhwgWYTCb4+flJyydOnAgPDw+7FUf2Yb3S7IihCBeKyxHopZG5IiIiIsdyR2eISktLUV5eLoWhU6dO4YMPPkBOTg6Cg4PtWiDdvQAvDdqHagEA24+z24yIiOh6dxSIhgwZgi+//BIAUFhYiPj4eLz77rsYOnQoFi1aZNcCyT56W+9HdJTdZkRERNe7o0D03//+F3369AEA/PDDDwgJCcGpU6fw5Zdf4l//+pddCyT7sF5+v4NXmhEREdVyR4HoypUr8Pb2BgBs2LABjz76KJRKJXr16oVTp07ZtUCyj+5t/KBUAKcuXoHBWCZ3OURERA7ljgJRVFQUfvzxR5w+fRrr16/HwIEDAQAFBQXQarV2LZDsw9vNFR3CfAAAO3mWiIiIyMYdBaJZs2bh5ZdfRps2bdCzZ0/o9XoA1WeL7r33XrsWSPbTM8IfALAz95LMlRARETmWOwpEjz32GPLy8rBnzx6sX79eWj5gwAC8//77diuO7CveGohO8AwRERHRte7oPkQAoNPpoNPppKfet2zZkjdldHDWM0THz5fwfkRERETXuKMzRBaLBXPnzoWPjw9at26N1q1bw9fXF2+99RYsFou9ayQ78fVQI0ZXPRh+F7vNiIiIJHd0hui1117Df/7zH7z99tvo3bs3AOC3337D7NmzUVZWhr///e92LZLsJz7CH0cMRdh54iIe6hQqdzlEREQO4Y4C0ZIlS/DZZ59JT7kHgLi4OLRo0QLPP/88A5EDi28bgCUZpziwmoiI6Bp31GV26dIlxMTE1FoeExODS5f4RevIerSpHkeUk1+EwisVMldDRETkGO4oEHXu3BkfffRRreUfffQR4uLi7rooajhB3hpEBnlCCI4jIiIisrqjLrP58+cjKSkJGzdulO5BlJGRgdOnT+OXX36xa4Fkf/FtA3D8fAl25V7CwA46ucshIiKS3R2dIerXrx9+//13DBs2DIWFhSgsLMSjjz6K7OxsfPXVV/aukewsnjdoJCIisqEQQgh7bWzfvn3o2rUrzGazvTbpEEwmE3x8fGA0GpvFo0kMxjL0mpcGpQLIenMgtG6ucpdERERkd7fz/X1HZ4ioadP5uKF1gAcsAsg8eVnucoiIiGTHQOSketZcbbaDD3olIiKSNxBt3boVjzzyCMLCwqBQKPDjjz/arH/qqaegUChspkGDBtm0uXTpEkaNGgWtVgtfX19MmDABxcXFNm3279+PPn36wM3NDeHh4Zg/f35D75rDi28bAADYeYLjiIiIiG7rKrNHH330pusLCwtv68NLSkrQuXNnPP300zfc9qBBg/DFF19I8xqN7fO3Ro0ahXPnziE1NRWVlZUYP348Jk6ciGXLlgGo7j8cOHAgEhISsHjxYhw4cABPP/00fH19MXHixNuqtzmxDqw+cMaIkvIqeGru+LF2RERETd5tfQv6+Pjccv3YsWPrvb3Bgwdj8ODBN22j0Wig09V9afjhw4exbt067N69G927dwcAfPjhh3jooYfwz3/+E2FhYVi6dCkqKirw+eefQ61Wo0OHDsjKysJ7773n1IEo3N8DLXzdcaawFJmnLqPvPUFyl0RERCSb2wpE156paSzp6ekIDg6Gn58fHnjgAfztb39DQEB1d09GRgZ8fX2lMAQACQkJUCqV2LlzJ4YNG4aMjAz07dsXarVaapOYmIh33nkHly9fhp+fX6Pvk6OIb+uPFf89g525FxmIiIjIqTn0oOpBgwbhyy+/RFpaGt555x1s2bIFgwcPli7rNxgMCA4OtnmPi4sL/P39YTAYpDYhISE2bazz1jbXKy8vh8lkspmao14RHEdEREQE3OGdqhvLiBEjpNedOnVCXFwcIiMjkZ6ejgEDBjTY586bNw9z5sxpsO07ivi21eOI9v1RiNIKM9zVKpkrIiIikodDnyG6Xtu2bREYGIhjx44BAHQ6HQoKCmzaVFVV4dKlS9K4I51Oh/z8fJs21vkbjU2aOXMmjEajNJ0+fdreu+IQWvl7QKd1Q6VZYG8e70dERETOq0kFoj/++AMXL15EaGgoAECv16OwsBCZmZlSm02bNsFisSA+Pl5qs3XrVlRWVkptUlNTER0dfcPxQxqNBlqt1mZqjhQKhXSWaAcf40FERE5M1kBUXFyMrKwsZGVlAQByc3ORlZWFvLw8FBcXY/r06dixYwdOnjyJtLQ0DBkyBFFRUUhMTAQAtG/fHoMGDcKzzz6LXbt2Ydu2bZg8eTJGjBiBsLAwAMCTTz4JtVqNCRMmIDs7G99++y0WLlyIadOmybXbDiVeGkfEGzQSEZHzkjUQ7dmzB/feey/uvfdeAMC0adNw7733YtasWVCpVNi/fz/+/Oc/45577sGECRPQrVs3/Prrrzb3Ilq6dCliYmIwYMAAPPTQQ7jvvvvw6aefSut9fHywYcMG5Obmolu3bnjppZcwa9Ysp77k/lrWM0R7TxeirLJ5PYOOiIiovuz6cNfmqrk93PVaQgj0+HsaLhSX49uJvaQ7WBMRETV1fLgr1du144h2chwRERE5KQYiQq8IayDiOCIiInJODEQkdZNlnrqMiiqLzNUQERE1PgYiQrtgL/h7qlFWacGBM4Vyl0NERNToGIioehxRTbfZDj7Gg4iInBADEQHANYGI44iIiMj5MBARAKB3VCAAYFfuJd6PiIiInA4DEQEAooK9EKLVoLzKgj0n+VwzIiJyLgxEBKB6HNF9UUEAgF+Pnpe5GiIiosbFQESSvvdUd5v9evSCzJUQERE1LgYikljHER06Z8KF4nKZqyEiImo8DEQkCfTSIDa0+lkv247xLBERETkPBiKy0acdu82IiMj5MBCRjfukQHQeQgiZqyEiImocDERko0cbf2hclMg3leNoQbHc5RARETUKBiKy4eaqkh72uvlIgczVEBERNQ4GIqplQEwwACDtMAMRERE5BwYiquWBmkC059QlFF6pkLkaIiKihsdARLWE+3sgOsQbFgGk5/Cu1URE1PwxEFGdBrSvPku08XC+zJUQERE1PAYiqtOA9iEAgC2/n0el2SJzNURERA2LgYjq1CXcFwGeahSVVWH3yUtyl0NERNSgGIioTiqlAv2ja7rNDvFqMyIiat4YiOiGHoyt7jZbn23gXauJiKhZYyCiG+ofHQQPtQpnCkux/w+j3OUQERE1GAYiuiE3V5V0T6JfDpyTuRoiIqKGw0BEN/VQp1AAwJoD59htRkREzRYDEd3U/dHBcHdV4Y/LpTh4xiR3OURERA2CgYhuyl19tdtsDbvNiIiomWIgoluydpv9wm4zIiJqphiI6JbujwmCm6sSeZeuYB+vNiMiomaIgYhuyUPtgoGxOgDAj3vPyFwNERGR/TEQUb0M69oCALBq31k+24yIiJodBiKqlz5RgQj0UuNiSQV+PXpe7nKIiIjsioGI6sVFpcQjncMAACv+y24zIiJqXhiIqN4evbclACD1UD5MZZUyV0NERGQ/DERUbx1baBEV7IXyKgvWHTDIXQ4REZHdMBBRvSkUCgy7t3pw9f/99w+ZqyEiIrIfBiK6LcPubQGFAtiZewknL5TIXQ4REZFdMBDRbQnzdUe/e4IAAN/uOS1zNURERPYhayDaunUrHnnkEYSFhUGhUODHH3+0WS+EwKxZsxAaGgp3d3ckJCTg6NGjNm0uXbqEUaNGQavVwtfXFxMmTEBxcbFNm/3796NPnz5wc3NDeHg45s+f39C71qyN6BEOAPgh8w/ek4iIiJoFWQNRSUkJOnfujI8//rjO9fPnz8e//vUvLF68GDt37oSnpycSExNRVlYmtRk1ahSys7ORmpqK1atXY+vWrZg4caK03mQyYeDAgWjdujUyMzOxYMECzJ49G59++mmD719zNaB9CAK91DhfVI5NRwrkLoeIiOjuCQcBQKxcuVKat1gsQqfTiQULFkjLCgsLhUajEd98840QQohDhw4JAGL37t1Sm7Vr1wqFQiHOnDkjhBDik08+EX5+fqK8vFxqM2PGDBEdHV3v2oxGowAgjEbjne5es/OPXw6J1jNWi/Ff7JK7FCIiojrdzve3w44hys3NhcFgQEJCgrTMx8cH8fHxyMjIAABkZGTA19cX3bt3l9okJCRAqVRi586dUpu+fftCrVZLbRITE5GTk4PLly/X+dnl5eUwmUw2E9l6ont1t1l6TgHOGUtlroaIiOjuOGwgMhiq73MTEhJiszwkJERaZzAYEBwcbLPexcUF/v7+Nm3q2sa1n3G9efPmwcfHR5rCw8PvfoeambZBXugZ4Q+LAL7fw0vwiYioaXPYQCSnmTNnwmg0StPp07yaqi4je1YHxW93n4bZImSuhoiI6M45bCDS6XQAgPz8fJvl+fn50jqdToeCAttBvVVVVbh06ZJNm7q2ce1nXE+j0UCr1dpMVNvgjqHw9XDFmcJSbObgaiIiasIcNhBFRERAp9MhLS1NWmYymbBz507o9XoAgF6vR2FhITIzM6U2mzZtgsViQXx8vNRm69atqKy8+uyt1NRUREdHw8/Pr5H2pnlyc1Xh8ZqxRF/uOCVzNURERHdO1kBUXFyMrKwsZGVlAageSJ2VlYW8vDwoFApMmTIFf/vb3/Dzzz/jwIEDGDt2LMLCwjB06FAAQPv27TFo0CA8++yz2LVrF7Zt24bJkydjxIgRCAurfjL7k08+CbVajQkTJiA7OxvffvstFi5ciGnTpsm0183L6PjWUCiArb+fRy7vXE1ERE1VI1z1dkObN28WAGpN48aNE0JUX3r/xhtviJCQEKHRaMSAAQNETk6OzTYuXrwoRo4cKby8vIRWqxXjx48XRUVFNm327dsn7rvvPqHRaESLFi3E22+/fVt18rL7mxv/xS7ResZqMefnbLlLISIiktzO97dCCMHRsLdgMpng4+MDo9HI8UR12JxTgPFf7Ia3mwt2/nUAPNQucpdERER0W9/fDjuGiJqOfu2C0DrAA0VlVfhx71m5yyEiIrptDER015RKBcb0ag0A+DLjJHjSkYiImhoGIrKLv3QLh5urEkcMRdh9su47gBMRETkqBiKyCx8PVwzt0gIA8MW2XJmrISIiuj0MRGQ343tHAADWZxuQd/GKzNUQERHVHwMR2U20zht97wmCRQCf8ywRERE1IQxEZFcT+7QFUP18s8IrFTJXQ0REVD8MRGRXvaMCEKPzRmmlGUt35sldDhERUb0wEJFdKRQKPFtzlmjJ9pOoqLLIXBEREdGtMRCR3T3SOQwhWg0Kisrx8z7eqJGIiBwfAxHZndpFiaf+VH3F2f9uPcEbNRIRkcNjIKIG8WR8K3iqVcjJL8KmIwVyl0NERHRTDETUIHzcXTFG3wYA8OGmYzxLREREDo2BiBrMM30i4OaqRNbpQmw7dlHucoiIiG6IgYgaTKCXBiN7tgIAfLjpqMzVEBER3RgDETWoiX3bQq1SYmfuJezKvSR3OURERHViIKIGFerjjse6twTAs0REROS4GIiowT3XLxIqpQK/Hr2AvXmX5S6HiIioFgYianDh/h4Ydm8LAMB7qb/LXA0REVFtDETUKF4c0A6uquqzRBnHecUZERE5FgYiahTh/h4Y0aP6irN/bsjhfYmIiMihMBBRo3nhgSi4uSqReeoy0nPOy10OERGRhIGIGk2w1g3jau5evWB9DiwWniUiIiLHwEBEjWpSv0h4aVxw6JwJaw8a5C6HiIgIAAMRNTI/TzWe6RMBoHosUaXZInNFREREDEQkg2f6tEWglxq5F0qwdMcpucshIiJiIKLG56VxwdQH7wEALEw7CmNppcwVERGRs2MgIlk80T0c7YK9cPlKJT7ZfEzucoiIyMkxEJEsXFRK/PWh9gCAL7adxOlLV2SuiIiInBkDEcmmf3QQekcFoMJswfz1OXKXQ0REToyBiGSjUCjw14faQ6EAVu07ywe/EhGRbBiISFYdwnwwvGtLAMDc1Yd4s0YiIpIFAxHJbnpiNDzVKuzNK8SKvWfkLoeIiJwQAxHJLkTrhhcGtAMAvL32CExlvAyfiIgaFwMROYSne0egbaAnLhSXY+HGo3KXQ0REToaBiByC2kWJN//cAQCQsv0kfs8vkrkiIiJyJgxE5DD63ROEgbEhMFsEZv+cDSE4wJqIiBoHAxE5lDcejoXGRYntxy9izYFzcpdDREROgoGIHEq4vwcm9YsEALy1+hAHWBMRUaNw6EA0e/ZsKBQKmykmJkZaX1ZWhuTkZAQEBMDLywvDhw9Hfn6+zTby8vKQlJQEDw8PBAcHY/r06aiqqmrsXaHb8Fz/SLQJ8EC+qRwL1vEO1kRE1PAcOhABQIcOHXDu3Dlp+u2336R1U6dOxapVq/D9999jy5YtOHv2LB599FFpvdlsRlJSEioqKrB9+3YsWbIEKSkpmDVrlhy7QvXk5qrCP4Z1AgB8vfMUMk9dkrkiIiJq7hw+ELm4uECn00lTYGAgAMBoNOI///kP3nvvPTzwwAPo1q0bvvjiC2zfvh07duwAAGzYsAGHDh3C119/jS5dumDw4MF466238PHHH6OiokLO3aJb+FNUIP7SrSWEAF79vwMorzLLXRIRETVjDh+Ijh49irCwMLRt2xajRo1CXl4eACAzMxOVlZVISEiQ2sbExKBVq1bIyMgAAGRkZKBTp04ICQmR2iQmJsJkMiE7O/uGn1leXg6TyWQzUeN7Lak9Ar3UOFpQjH9vOSF3OURE1Iw5dCCKj49HSkoK1q1bh0WLFiE3Nxd9+vRBUVERDAYD1Go1fH19bd4TEhICg8EAADAYDDZhyLreuu5G5s2bBx8fH2kKDw+3745Rvfh6qDHrkep7E3206RiOFfDeRERE1DAcOhANHjwYf/nLXxAXF4fExET88ssvKCwsxHfffdegnztz5kwYjUZpOn36dIN+Ht3YI3GhuD86CBVmC6Z+uw+VZovcJRERUTPk0IHoer6+vrjnnntw7Ngx6HQ6VFRUoLCw0KZNfn4+dDodAECn09W66sw6b21TF41GA61WazORPBQKBd4eHgcfd1ccOGPEh5uOyV0SERE1Q00qEBUXF+P48eMIDQ1Ft27d4OrqirS0NGl9Tk4O8vLyoNfrAQB6vR4HDhxAQUGB1CY1NRVarRaxsbGNXj/dmRCtG/4+rCMA4OPNx7A377LMFRERUXPj0IHo5ZdfxpYtW3Dy5Els374dw4YNg0qlwsiRI+Hj44MJEyZg2rRp2Lx5MzIzMzF+/Hjo9Xr06tULADBw4EDExsZizJgx2LdvH9avX4/XX38dycnJ0Gg0Mu8d3Y6H48IwpEsYzBaBad/tw5UK3kuKiIjsx6ED0R9//IGRI0ciOjoajz/+OAICArBjxw4EBQUBAN5//308/PDDGD58OPr27QudTocVK1ZI71epVFi9ejVUKhX0ej1Gjx6NsWPHYu7cuXLtEt2FuX/uiFAfN+ReKMHsn298lSAREdHtUgg+QfOWTCYTfHx8YDQaOZ5IZtuPX8Doz3bCIoAFj8XhL915BSAREdXtdr6/HfoMEdH1/hQZiGkP3gMAeOOngzh8jveIIiKiu8dARE3O8/2j0D86CGWVFjy/9L8o4gNgiYjoLjEQUZOjVCrw/uNdEFYznuiFb/aiivcnIiKiu8BARE2Sn6ca/x7THW6uSqTnnMdbqw/JXRIRETVhDETUZHVq6YMPnugCAFiScQpLtp+UtR4iImq6GIioSRvUMRQzBsUAAOasysa6g+dkroiIiJoiBiJq8ib1a4sRPcJhEcAL3+zF1t/Py10SERE1MQxE1OQpFAr8fVgnPNRJh0qzwMSv9mD3yUtyl0VERE0IAxE1CyqlAh88cS/63VN9Of7TX+zGHoYiIiKqJwYiajbULkosHt0N8RH+KCqvwpj/7MK2YxfkLouIiJoABiJqVtzVKqSM74k+7QJRWmnG+JTd2HQkX+6yiIjIwTEQUbPjrlbhs3Hd8WBsCCqqLJj4ZSa+33Na7rKIiMiBMRBRs6RxUeGTUV0xtEsYqiwC03/Yj3c35IDPMiYiorowEFGz5apS4r3HuyD5/kgAwIebjmHKt1korzLLXBkRETkaBiJq1pRKBaYnxmD+8Di4KBX4KessRny6A+eMpXKXRkREDoSBiJzC4z3CseTpntC6uWBvXiEe+fA3ZBy/KHdZRETkIBiIyGn0jgrE6hf6oH2oFheKKzD6Pzvx6dbjHFdEREQMRORcWgV4YMVzf8Kj97aA2SLwj1+O4Jkle3ChuFzu0oiISEYMROR03NUqvPt4Z7w1pAPUKiXSjhRg0Ae/YgufgUZE5LQYiMgpKRQKjNG3wU+Te6NdsBcuFJdj3Oe7MHfVIZRV8io0IiJnw0BETq19qBarXrgP4/StAQCfb8vFQ//6lc9BIyJyMgxE5PTcXFWYM6QjPn+qO4K8NThxvgR/+XcGZv+cjZLyKrnLIyKiRsBARFTjgZgQbJzaD493bwkhgJTtJzHw/a3YfKRA7tKIiKiBMRARXcPHwxXzH+uMryfEo6WfO84UlmJ8ym48nbIbuRdK5C6PiIgaCAMRUR3uaxeI9VP64n/6toWrSoFNRwow8P0teHvtERSzG42IqNlRCN6V7pZMJhN8fHxgNBqh1WrlLoca2fHzxZiz6hC21lyWH+CpxvP3R2FUfCu4uapkro6IiG7kdr6/GYjqgYGIhBBIO1yAv605hJMXrwAAQn3c8OKAdhjerSVcVTzZSkTkaBiI7IyBiKwqzRb8X+YfWJh2FOeMZQCAFr7uePq+CIzoEQ5PjYvMFRIRkRUDkZ0xENH1yirNWLozD4vSj+FCcQUAwMfdFaN7tcLoXq0R6uMuc4VERMRAZGcMRHQjZZVmrPjvGXz26wmcqLkKTakAHogJxogerdA/Oggu7E4jIpIFA5GdMRDRrVgsAqmH8/Gf33KxK/fqXa51Wjc80jkUD8eFIa6lDxQKhYxVEhE5FwYiO2MgottxrKAY3+7Oww+Zf+DylUppeSt/DyTFhWJgbAjiWvpCpWQ4IiJqSAxEdsZARHeivMqMzUcKsHr/OaQdLkDpNQ+N9fNwRZ92QegfHYS+9wQh0EsjY6VERM0TA5GdMRDR3bpSUYVNRwqw9oABW4+eR1GZ7c0d2wZ5omcbf3Rv448ebfzQyt+D3WtERHeJgcjOGIjInqrMFuw9XYj0nAKk55xH9llTrTYBnmp0aOGDDmFadAyr/rOVvweU7GYjIqo3BiI7YyCihlR4pQKZpy5j98nL2H3yEvb/UYhKc+1fSy+NCyKDPBEZ5IW2NX9GBnuhdYAHNC68YzYR0fUYiOyMgYgaU1mlGUcMRcg+a8TBMyYcOmvEYUMRKqosdbZXKoAQrRta+rmjha87Wvp5oIX02h0hWjfeMJKInBIDkZ0xEJHcKs0W5F4owYnzxTh+vgTHC4px/HwxTpwvQVE9HjbrqVYhyFsjTcHebtWvvTQI8FLD10MNPw9X+Hqo4ePuyivgiKhZuJ3vb/63kagJcFUpcU+IN+4J8bZZLoTA+eJy/HG5FGcul+JMYSn+uHxFen3mcilKKszV08Ur0nPYbkXr5gI/TzV83atDkq+HK3zcXeGlcYGXm0v1n9bJzfa1t8YVbq5KDgonoibFqQLRxx9/jAULFsBgMKBz58748MMP0bNnT7nLIrpjCoUCwd5uCPZ2Q9dWfnW2KSmvwvmichQUldf8WYbz0utyXCqpwOUrFTBeqZTONpnKqmAqq8KpO6xLqQA8NS7wUKvg7qqCW83k7qqC+zXL3NVKuLlUL7t+vcZFCVeVEmqXayaVEprr5tU17TQuDGFEdOecJhB9++23mDZtGhYvXoz4+Hh88MEHSExMRE5ODoKDg+Uuj6jBeGpc4KlxQZtAz1u2rTRbYCytROGVShReqUDhlcrqsFRaCWNpJYrLq1BcVoWSiioUlVVdnS+vQlF59bwQgEUARWVVtW4v0NBcVQopJNmGJhXUKgVUSgVcVEq4qhRQKZVwVVYvc1Up4VKz3lVZ/dqlpq30WnnN+2qWqVTV23BRKWvaVy9XKqq3pVQqoLK+rvlTpcTV9dKya14rFFAqcc3rq38qFaijLUMgkT04zRii+Ph49OjRAx999BEAwGKxIDw8HC+88AJeffXVm76XY4iI6kcIgdJKM4rLqgNSWaUZZZVmlFZYUFppRmmlGWUVZul1aUXNemu7SgtKK8woraxCeaUFFWYLKqpqJvN1f1ZZUGVxin++bskmSFlD1HWBSqGoDmLWP5WK6jOM184ra86wKWu2Vd1eAQUgrVcqFIDCdl5Rsy3l9du+ZlsKa9trtqWQ2trOK+uoVwFAqaz5LFj/RK151LSts13N/t1sG9Xrr93G1eN0y+2jekZx/TZuun3bbcBmvvY2rNuXjrvNvl+/vdrbxjX7dW0t1x4XSO1h89p2H2zbKmza3tnnuKgUdn8wNscQXaeiogKZmZmYOXOmtEypVCIhIQEZGRkyVkbUvCgUCnioXeChdkFjnHe1WAQqzBaUXxOaKq8JTdcurzJbUGkWqLJYYLYIVJoFzJaaZebqcFVlETXrLKgyV89fXWe7rNIiYK7ZXvWy6veZLQJmIWCp+dNsgfT66rLr1tdaJq4uq0fmM1sEzBCA+dZtiRxVsLcGu15LkO3znSIQXbhwAWazGSEhITbLQ0JCcOTIkVrty8vLUV5eLs2bTLVvnEdE8lMqFXBTVo8/aq6EqA5F1pB0NTyhVniyDVLVYcv6WghAoHpb1fM1ry0CApDaWGqWi+vmpW1cM2+zXYvtNqzbtG7r6udYt3Gz7QK4bl7UbNNsufqZ1v4N67prl1vnYZ2vY51A9YyotY2r87j2fTfb/g22gevrum4bqLVN222grn2rOS61tl2rbtvt1VR79X3XhG2p1hu1ldpZW12/vI73W2u5+hbb43LdZ8n9e+wUgeh2zZs3D3PmzJG7DCIiKBQKqGrGDhFRw1HKXUBjCAwMhEqlQn5+vs3y/Px86HS6Wu1nzpwJo9EoTadPn26sUomIiEgGThGI1Go1unXrhrS0NGmZxWJBWloa9Hp9rfYajQZardZmIiIioubLabrMpk2bhnHjxqF79+7o2bMnPvjgA5SUlGD8+PFyl0ZEREQyc5pA9MQTT+D8+fOYNWsWDAYDunTpgnXr1tUaaE1ERETOx2nuQ3Q3eB8iIiKipud2vr+dYgwRERER0c0wEBEREZHTYyAiIiIip8dARERERE6PgYiIiIicHgMREREROT0GIiIiInJ6DERERETk9BiIiIiIyOk5zaM77ob1Zt4mk0nmSoiIiKi+rN/b9XkoBwNRPRQVFQEAwsPDZa6EiIiIbldRURF8fHxu2obPMqsHi8WCs2fPwtvbGwqFwq7bNplMCA8Px+nTp/mctFvgsao/Hqv647G6PTxe9cdjVX8NdayEECgqKkJYWBiUypuPEuIZonpQKpVo2bJlg36GVqvlL0w98VjVH49V/fFY3R4er/rjsaq/hjhWtzozZMVB1UREROT0GIiIiIjI6TEQyUyj0eDNN9+ERqORuxSHx2NVfzxW9cdjdXt4vOqPx6r+HOFYcVA1EREROT2eISIiIiKnx0BERERETo+BiIiIiJweAxERERE5PQYiGX388cdo06YN3NzcEB8fj127dsldUqObPXs2FAqFzRQTEyOtLysrQ3JyMgICAuDl5YXhw4cjPz/fZht5eXlISkqCh4cHgoODMX36dFRVVTX2rtjd1q1b8cgjjyAsLAwKhQI//vijzXohBGbNmoXQ0FC4u7sjISEBR48etWlz6dIljBo1ClqtFr6+vpgwYQKKi4tt2uzfvx99+vSBm5sbwsPDMX/+/IbeNbu71bF66qmnav2cDRo0yKaNsxyrefPmoUePHvD29kZwcDCGDh2KnJwcmzb2+r1LT09H165dodFoEBUVhZSUlIbePbuqz7Hq379/rZ+tSZMm2bRxhmMFAIsWLUJcXJx0c0W9Xo+1a9dK6x3+50qQLJYvXy7UarX4/PPPRXZ2tnj22WeFr6+vyM/Pl7u0RvXmm2+KDh06iHPnzknT+fPnpfWTJk0S4eHhIi0tTezZs0f06tVL/OlPf5LWV1VViY4dO4qEhASxd+9e8csvv4jAwEAxc+ZMOXbHrn755Rfx2muviRUrVggAYuXKlTbr3377beHj4yN+/PFHsW/fPvHnP/9ZREREiNLSUqnNoEGDROfOncWOHTvEr7/+KqKiosTIkSOl9UajUYSEhIhRo0aJgwcPim+++Ua4u7uLf//73421m3Zxq2M1btw4MWjQIJufs0uXLtm0cZZjlZiYKL744gtx8OBBkZWVJR566CHRqlUrUVxcLLWxx+/diRMnhIeHh5g2bZo4dOiQ+PDDD4VKpRLr1q1r1P29G/U5Vv369RPPPvuszc+W0WiU1jvLsRJCiJ9//lmsWbNG/P777yInJ0f89a9/Fa6uruLgwYNCCMf/uWIgkknPnj1FcnKyNG82m0VYWJiYN2+ejFU1vjfffFN07ty5znWFhYXC1dVVfP/999Kyw4cPCwAiIyNDCFH9RahUKoXBYJDaLFq0SGi1WlFeXt6gtTem67/kLRaL0Ol0YsGCBdKywsJCodFoxDfffCOEEOLQoUMCgNi9e7fUZu3atUKhUIgzZ84IIYT45JNPhJ+fn82xmjFjhoiOjm7gPWo4NwpEQ4YMueF7nPVYCSFEQUGBACC2bNkihLDf790rr7wiOnToYPNZTzzxhEhMTGzoXWow1x8rIaoD0YsvvnjD9zjrsbLy8/MTn332WZP4uWKXmQwqKiqQmZmJhIQEaZlSqURCQgIyMjJkrEweR48eRVhYGNq2bYtRo0YhLy8PAJCZmYnKykqb4xQTE4NWrVpJxykjIwOdOnVCSEiI1CYxMREmkwnZ2dmNuyONKDc3FwaDwebY+Pj4ID4+3ubY+Pr6onv37lKbhIQEKJVK7Ny5U2rTt29fqNVqqU1iYiJycnJw+fLlRtqbxpGeno7g4GBER0fjueeew8WLF6V1znysjEYjAMDf3x+A/X7vMjIybLZhbdOU/427/lhZLV26FIGBgejYsSNmzpyJK1euSOuc9ViZzWYsX74cJSUl0Ov1TeLnig93lcGFCxdgNptt/tIBICQkBEeOHJGpKnnEx8cjJSUF0dHROHfuHObMmYM+ffrg4MGDMBgMUKvV8PX1tXlPSEgIDAYDAMBgMNR5HK3rmivrvtW179cem+DgYJv1Li4u8Pf3t2kTERFRaxvWdX5+fg1Sf2MbNGgQHn30UUREROD48eP461//isGDByMjIwMqlcppj5XFYsGUKVPQu3dvdOzYEQDs9nt3ozYmkwmlpaVwd3dviF1qMHUdKwB48skn0bp1a4SFhWH//v2YMWMGcnJysGLFCgDOd6wOHDgAvV6PsrIyeHl5YeXKlYiNjUVWVpbD/1wxEJGsBg8eLL2Oi4tDfHw8Wrduje+++65J/SNAjm3EiBHS606dOiEuLg6RkZFIT0/HgAEDZKxMXsnJyTh48CB+++03uUtxeDc6VhMnTpRed+rUCaGhoRgwYACOHz+OyMjIxi5TdtHR0cjKyoLRaMQPP/yAcePGYcuWLXKXVS/sMpNBYGAgVCpVrdH1+fn50Ol0MlXlGHx9fXHPPffg2LFj0Ol0qKioQGFhoU2ba4+TTqer8zha1zVX1n272c+QTqdDQUGBzfqqqipcunTJ6Y9f27ZtERgYiGPHjgFwzmM1efJkrF69Gps3b0bLli2l5fb6vbtRG61W2+T+s3OjY1WX+Ph4ALD52XKmY6VWqxEVFYVu3bph3rx56Ny5MxYuXNgkfq4YiGSgVqvRrVs3pKWlScssFgvS0tKg1+tlrEx+xcXFOH78OEJDQ9GtWze4urraHKecnBzk5eVJx0mv1+PAgQM2X2apqanQarWIjY1t9PobS0REBHQ6nc2xMZlM2Llzp82xKSwsRGZmptRm06ZNsFgs0j/aer0eW7duRWVlpdQmNTUV0dHRTbILqL7++OMPXLx4EaGhoQCc61gJITB58mSsXLkSmzZtqtUNaK/fO71eb7MNa5um9G/crY5VXbKysgDA5mfLGY7VjVgsFpSXlzeNn6u7HpZNd2T58uVCo9GIlJQUcejQITFx4kTh6+trM7reGbz00ksiPT1d5Obmim3btomEhAQRGBgoCgoKhBDVl2m2atVKbNq0SezZs0fo9Xqh1+ul91sv0xw4cKDIysoS69atE0FBQc3isvuioiKxd+9esXfvXgFAvPfee2Lv3r3i1KlTQojqy+59fX3FTz/9JPbv3y+GDBlS52X39957r9i5c6f47bffRLt27WwuJS8sLBQhISFizJgx4uDBg2L58uXCw8OjyV1KfrNjVVRUJF5++WWRkZEhcnNzxcaNG0XXrl1Fu3btRFlZmbQNZzlWzz33nPDx8RHp6ek2l4pfuXJFamOP3zvr5dHTp08Xhw8fFh9//HGTu5T8Vsfq2LFjYu7cuWLPnj0iNzdX/PTTT6Jt27aib9++0jac5VgJIcSrr74qtmzZInJzc8X+/fvFq6++KhQKhdiwYYMQwvF/rhiIZPThhx+KVq1aCbVaLXr27Cl27Nghd0mN7oknnhChoaFCrVaLFi1aiCeeeEIcO3ZMWl9aWiqef/554efnJzw8PMSwYcPEuXPnbLZx8uRJMXjwYOHu7i4CAwPFSy+9JCorKxt7V+xu8+bNAkCtady4cUKI6kvv33jjDRESEiI0Go0YMGCAyMnJsdnGxYsXxciRI4WXl5fQarVi/PjxoqioyKbNvn37xH333Sc0Go1o0aKFePvttxtrF+3mZsfqypUrYuDAgSIoKEi4urqK1q1bi2effbbWfz6c5VjVdZwAiC+++EJqY6/fu82bN4suXboItVot2rZta/MZTcGtjlVeXp7o27ev8Pf3FxqNRkRFRYnp06fb3IdICOc4VkII8fTTT4vWrVsLtVotgoKCxIABA6QwJITj/1wphBDi7s8zERERETVdHENERERETo+BiIiIiJweAxERERE5PQYiIiIicnoMREREROT0GIiIiIjI6TEQERERkdNjICIiukMKhQI//vij3GUQkR0wEBFRk/TUU09BoVDUmgYNGiR3aUTUBLnIXQAR0Z0aNGgQvvjiC5tlGo1GpmqIqCnjGSIiarI0Gg10Op3NZH3yvEKhwKJFizB48GC4u7ujbdu2+OGHH2zef+DAATzwwANwd3dHQEAAJk6ciOLiYps2n3/+OTp06ACNRoPQ0FBMnjzZZv2FCxcwbNgweHh4oF27dvj5558bdqeJqEEwEBFRs/XGG29g+PDh2LdvH0aNGoURI0bg8OHDAICSkhIkJibCz88Pu3fvxvfff4+NGzfaBJ5FixYhOTkZEydOxIEDB/Dzzz8jKirK5jPmzJmDxx9/HPv378dDDz2EUaNG4dKlS426n0RkB3Z5RCwRUSMbN26cUKlUwtPT02b6+9//LoSoflL5pEmTbN4THx8vnnvuOSGEEJ9++qnw8/MTxcXF0vo1a9YIpVIpDAaDEEKIsLAw8dprr92wBgDi9ddfl+aLi4sFALF27Vq77ScRNQ6OISKiJuv+++/HokWLbJb5+/tLr/V6vc06vV6PrKwsAMDhw4fRuXNneHp6Sut79+4Ni8WCnJwcKBQKnD17FgMGDLhpDXFxcdJrT09PaLVaFBQU3OkuEZFMGIiIqMny9PSs1YVlL+7u7vVq5+rqajOvUChgsVgaoiQiakAcQ0REzdaOHTtqzbdv3x4A0L59e+zbtw8lJSXS+m3btkGpVCI6Ohre3t5o06YN0tLSGrVmIpIHzxARUZNVXl4Og8Fgs8zFxQWBgYEAgO+//x7du3fHfffdh6VLl2LXrl34z3/+AwAYNWoU3nzzTYwbNw6zZ8/G+fPn8cILL2DMmDEICQkBAMyePRuTJk1CcHAwBg8ejKKiImzbtg0vvPBC4+4oETU4BiIiarLWrVuH0NBQm2XR0dE4cuQIgOorwJYvX47nn38eoaGh+OabbxAbGwsA8PDwwPr16/Hiiy+iR48e8PDwwPDhw/Hee+9J2xo3bhzKysrw/vvv4+WXX0ZgYCAee+yxxttBImo0CiGEkLsIIiJ7UygUWLlyJYYOHSp3KUTUBHAMERERETk9BiIiIiJyehxDRETNEkcDENHt4BkiIiIicnoMREREROT0GIiIiIjI6TEQERERkdNjICIiIiKnx0BERERETo+BiIiIiJweAxERERE5PQYiIiIicnr/H0vLd8Ye+6tHAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "W, H, row_bias, col_bias, losses = model.fit(Xtrain, MASK, plot_loss=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "Xtilda = model.predict(W, H, row_bias, col_bias)\n",
    "Xtilda_bool = model.map_probabilities_to_binary(Xtilda, threshold=0.5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "test_mask = np.zeros(Xtrain.shape).astype(\"bool\")\n",
    "for x, y in removed_coords:\n",
    "    test_mask[x,y] = True"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "roc_auc = 1.0\n",
      "pr_auc = 1.0\n",
      "accuracy = 1.0\n",
      "f1 = 1.0\n",
      "precision = 1.0\n",
      "recall = 1.0\n",
      "mcc = 1.0\n",
      "rmse = 0.0\n"
     ]
    }
   ],
   "source": [
    "for key, value in get_scores(X[test_mask], Xtilda_bool[test_mask]).items():\n",
    "    print(key, \"=\", value)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "roc_auc = 1.0\n",
      "pr_auc = 1.0\n",
      "accuracy = None\n",
      "f1 = None\n",
      "precision = None\n",
      "recall = None\n",
      "mcc = None\n",
      "rmse = 0.009379555872214364\n"
     ]
    }
   ],
   "source": [
    "for key, value in get_scores(X[test_mask], Xtilda[test_mask], binary=False).items():\n",
    "    print(key, \"=\", value)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "TELF",
   "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.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}