Skip to content

Instantly share code, notes, and snippets.

@leelasd
Forked from greglandrum/MCS Pharmacophore.ipynb
Created October 27, 2020 22:58
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save leelasd/0a496dbee95c37057cf3e9bbfa294a45 to your computer and use it in GitHub Desktop.
Save leelasd/0a496dbee95c37057cf3e9bbfa294a45 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"from rdkit import Chem\n",
"from rdkit.Chem import rdFMCS\n",
"from rdkit.Chem.Draw import IPythonConsole\n",
"from rdkit.Chem import Draw"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAyAAAAJYCAYAAACadoJwAADB0klEQVR4nOzdeVhU1RsH8HeGHVnE\nfV/Q3HfUVMwl0Uz5ZalYLpRbqKWYplGaYqaGa6RZkmmhloVpiaUVaZp7abmbC+77hgKyM9/fH7d7\nBUWdYYYL6PfzPPOgMPfMmeXMve9Z3mMAACEiIiIiItKBMb8rQEREREREjw8GIEREREREpBsGIERE\nREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIERERERE\npBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsG\nIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIERE\nREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIERERERE\npBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsG\nIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIERE\nREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIERERERE\npBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsG\nIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIEREREREpBsGIERE\nREREpBsGIEREREREpBsGIEREREREpBsGIERERERU4GzaJJKRIXL0qEh6en7XhmzJPr8rQERERESk\nun5dZMMGkc2bRf79V8TRUaRMGREHB/3qcPv2bfnnn3/Ezc1NPDw8xNPTU9zd3cXR0VG/SjzCGIAQ\nERERUb7JyBDZu1fkt9+U28aNyu969BBp2lQkNlbf+vz++++yevVq+eijj+75m6Ojo7i5uUnRokXF\n3d1d3NzcxN3dXdzd3aVo0aLi5uYms2bNEqORk4wexAAA+V0JIiIiIno8ACL79onExCgBx+bNIklJ\nd/7u7Czi6ytSubLI55+LPPOMyIIFIjduKAFJXkpKSpIGDRpIbGysVKtWTYoUKSKJiYkSFxcniYmJ\nkv6QuWD29vbSrFkzWbZsmXh7e+dtZQsxBiBERERElKc2bVKCioMHRTp3Frl06c7fDAaRRo1E/PxE\nOnYUad1axMXlzt8PHxZp21YJTA4cEPHwyLt6vvnmmzJnzhxp0KCB7Nq1S27fvi07duzQRjscHR3F\n3t5eAMjt27clMTFREhMTJSEhQW7evClffPGF7NixQwYMGCCLFy/Ou4oWcgxAiIiIiChPjR8vUqmS\nsp5j9mxlncdTTylBR9euIuXL3/9Yk0m577ZtIkOHinz6ad7U8c8//5RWrVqJwWCQHTt2iI+Pj+zY\nsUNatmyZ4/2dnJzEzc1NPD09pVOnTvLpp5/K6dOnpUaNGpKZmSkHDx6UmjVr5k1lCzlOUCMiIiKi\nPFWypDKqceqUyE8/iVy8KBIVJRIU9ODgQ0Rk505lCpazs0hEhDJty9bS0tJk0KBBkpmZKWPGjBEf\nHx8RESlSpIh06tRJWrVqJfXr15eqVauKl5eX2NnZSWpqqly/fl1OnDghly9fFhGRypUrS//+/SUz\nM1OmTJli+4o+IjgCQkRERER5auRIEX9/kWnTRFavNn8a1dSpIhMmiEycqEzVmjRJpGpVkf37RYoU\nsV39Jk2aJO+9957UqFFD9uzZIy5Z54DdR3Jysjb9ymAwSNWqVUVE5MyZM1KjRg1JT0+XPXv2SP36\n9W1X0UcER0CIiIiIKE/t3SvSqZNIaKhlaziefloJPKZOFXn2WZH69UWuXEmUefNW26xu//77r4SF\nhYnRaJSFCxeaFXyIiLi4uEjJkiXF29tbCz5ERCpVqiSDBw8Wk8kkU6dOtVk9HyUMQIiIiIgoTyUk\nKD/d3UVWrhRZt05Z2/EwLVsq6z4yMkSGDxdZtChZihdvKBMm9JRdu3ZZXS+TySSDBg2S1NRUGTp0\nqLRp08bqMkVExo8fLy4uLhIVFSV79+61SZmPEgYgRERERJSn4uOVn66uIgEBynQsg8G8Y2fMEGnQ\n4KY4Ow+RP/74VF58sYdkZGTIwIEDJS0tzap6ffjhh7Jt2zYpV66cTUcrypYtK0OGDBEAMnnyZJuV\n+6jgGhAiIiIiylNlyohcvixy5IhIzZoinp4iN2+af/y6dT9Lly7Piqurq+zcuVO6d+8ux44dk8mT\nJ8uECRNyVadTp05J/fr1JTExUVavXi3PPfdcrsq5n0uXLkm1atUkOTlZdu7cKc2aNbNp+YUZR0CI\niIiIKE+pU7DUaVfu7pYd/+yznSUwMFCSkpJk+PDh8tlnn4nBYJApU6bI9u3bLa4PAAkKCpLExETp\n16+fzYMPEZEyZcrI66+/LgDk/ffft3n5hRlHQIiIiIgoz2RmmqRx40RJT3eQlSvtxdfXXipVEtm7\n18w5WP+5fv261K1bVwDItm3bZPz48fLdd99JZmamiIi4ublpGwZ6eHiIp6en9n83NzcpWrSoeHh4\niJubm/z9998SGRkpJUqUkIMHD0qpUqXy4qnLtWvXxNvbWxISEmTHjh3y5JNP5snjFDb2+V0BIiIi\nInp0JSTEy/79XuLp6Snx8T/LzZstpWbNJ0Vkh0XlFC9eXL7//nt54oknxNHRUf755x/JzMwUFxcX\nLSVuYmKi2eW5ublJaGioFnwkJiaKm5ubRXV6mBIlSsiIESNk2rRpMmnSJFm3bp1Nyy+s9A1A0tKU\nLTBNJhGj8c5PIiIiInokxf+3At3d3T3bv3OjZcuWYjKZpFu3bnL06FFp0KCBbNu2TYoUKSLx8fHa\nvhyJiYkSFxeX7f/x8fFy69YtSUhIkJ9//lmOHz8ue/fulfPnz8urr74q165dk507d4rB3NXxDwBA\nUlNTxdnZWd58802ZO3euxMTESJ06daR48eLaaE3RokWzjdx4enpqozTq7ytUqCClS5e2uk4FiX4B\nyKpVSsBx7pyIk5PIsGHKdpbDhulWBSJ6zF29KnLggMj58yL29kqnyMsv53etiAoGtg/KIwn/LQDx\n8PDI9u/cmjBhgvz4449SrFgxWbVqlRT5b0dCDw8Ps8s9fvy41K5dWyIjI2XUqFGyb98+OX/+vKxZ\ns8Ym60EWLVokM2fOlEWLFknNmjXFaDSKp6enHD582OKyunfvLn379pXu3btbXa+CQr8A5MIFJYFz\nZKSSAmHZMpEzZ3R7eCJ6fFy5olxHHT4s4nr+mAzYHqT84to15Q69e4t8/bXSCUL0mGH7IL2pQYct\nRkC2rF4tH3zwgdjb20tUVJRUq1YtV+VUr15dXn75ZVm8eLHMmDFD3nrrLRk5cqS8++674u/vL0Yr\nZuhcvHhR3nrrLYmLi5Pz58/L/PnzJT4+Xp566ikJDw+X27dva6MyN2/e1P6dkJAgCQkJcvPmTW06\nWVxcnKxZs0bWrVsnx48fl3LlyuW6XgWJfgFI2bIiX30lkpoqUrGiSL9+Ip9+qtvDP9YyM0UAZQTK\n0TG/a0N3W7tWSZCelqa8R/37i3z5pfKTzDJqlMiYMSKbNomsWCHyww93/lareFEZcH2j8h8PD5Ha\ntUVKlxb55hsRFxeR69dFihfPj2qTOdg+rMb2Qfkta9DRo0cPadCggRQtWtTygvbtk9YvvyxrWreW\n4z16SIcOHayq18SJE2XZsmWybNkyGT16tFSsWFH2798v33//vfTo0SPX5b7++usSFxcnXbt2FTc3\nN/nmm2/E1dVVvvzyS/H29ra4vICAAPnuu+8kLCxM5s6dm+t6FST6LcDo0UMJPiZPvtPLwulXufZf\nwgcxmZTUdgcPKufpiAiRXbM2KsPmbdqIVKminEQmTRKJjhYJD8+/SlPOTp0SeeklkeRkkdOnRZYv\nFzl7Nr9rVajUr6/M8kxLE6lUSckv37KlyODBIq+OKyn4NUYZcb11S2THDpEPPxTp2VP5XaVKIhs3\n5km90tPTZc+ePbJ06VL5888/5aYlSe9JwfZhNbYPyk9JSUnyxRdfSNGiRWX79u0ybdo08fLysnzk\n4sYNke7dReLjpWuVKjJy5Eir61a5cmUZMGCAZGZmyvTp02XcuHEiIhIaGiomc7Zpz8E333wj33//\nvXh4eMisWbNk6NChIiLywQcf5Cr4EBF57733xGg0SkREhJw6dSpXZRQ40NN77wEiwLvv6vqwj6Jh\nw4DDh4EPPlBe0qy3Bb5Lsv/CYAD69VMO/PJLIDExfytP2a1YASxfDkRGAhERyu/Un2SWxYuBI0eA\n0aOBlBQLDpwyRWkjrVrZvE5paWlo06YNWrZsCRHRbl5eXvDx8UFAQABCQ0MRFRWFXbt2ISkpyeZ1\neCSwfViN7YMe5MqVK5g/fz7mz5+PzZs327TsX3/9FVWqVIGIwMXFRXufjUYjOnbsiOXLlyM5Ofnh\nBWVkAJ07K5/HJk0AG34ezpw5AycnJxiNRvz999+oWrUqRATLly+3uKxr166hdOnSEBF89tlnCAoK\ngojgySefREZGhlX17N27N0QEQ4cOtaqcgkLfAOS115QPz7x5uj7so2jhQuDDD4FPPwXc3IDatYFn\nngEGDwZWzT6hnKx//x04fhxITQW++w5YtgwYP155H9LS8vspUFbTpwOVKwOzZ+d3TQql2bOBcuWA\njz6y8MDERKB0aeV7ae1am9ZJPfGULFkSL7zwAho3bowiRYpku9jKerO3t8cTTzyBt99+26b1eCSw\nfViF7YPuJyoqCiVLloSIwMHBASICX19fbNy40apyb968iaCgIBgMBogIGjRogJ07d+K3335Dnz59\n4OzsrL23yzp1Al5/Hdi9+/4FjhypfA7LlAHOnrWqbjkZPnw4RAQBAQH4/PPPISJ44oknkJ6eblE5\n/fr1g4igXbt22LBhAwwGA5ycnHDgwAGr6/jvv//Czs4Obm7FcfLkLavLy2/6BiDduysfoBUrdH3Y\nR9GcOcC5c8CgQRYclJamRCoiSjBCBcfbbyvvy7Rp+V2TQikoSHn5PvooF7H1nDl3etVMJpvUZ8GC\nBRARODs746+//oLJZIK/vz9ef/11hIWFYfbs2Zg2bRreeustBAQEoE6dOrCzs4OIoEaNGvDz87NJ\nPR4ZbB9WYfugu124cAEvvPCCFgT4+vpixIgR8PT01H7XtWtX7Nq1y+Ky16xZgwoVKmjvcWhoKFJT\nU7PdJy4uDvPnz0fTpk1xW70uEQHq1AHCwoArV+7cOSMDCAgAnJyArVutfeo5unDhAlxdXWEwGPD3\n33+jRo0aEBEsWbLE7DLWrl0LEYGrqyv27duH6tWrQ0Tw/vvv26yeISE/o1y5DAwcaLMi842+AUir\nVsoH7I8/dH3YR821a8rLWLNmLg7+/HPl4AYNbHYyySozMxM7d+7EoEGD8NZbb+Hbb7/F8ePHbf44\nj5zhw5X3Ze7c/K5JodS2rfLyTZ8O2NsDvXtbcHByMlCpEk76+uK377+3ui7bt2+Hk5MTRASLFy8G\nAJw7dy7HXl1nZ2fUr18fPXv2xIoVK7B7927tAsAWPWaPDLYPq7B9kMpkMiEyMhLFihWDiMDT0xMR\nEREw/Xc9cPbsWYSFhWULRPz8/MwKRC5duoTAwEDtuFatWuHw4cMPr9TevcoIR4kSdwIRJydgwwYg\nOloZGTGZgH/+se7JP8SoUaMgInj++eexZMkSiAg6dOhg1rHx8fGoVKkSRASzZ8/Gm2++qY38pNlw\nxsnJk4CjI2BnB/z7r82KzRe6BiDfvvACNrRrh1tHj+r5sI+cn39W2udTT+WiNyslRRmLF1EKsrGR\nI0fC3t4eRqMx24nEw8MDvr6+CA4ORmRkJA4cOIDMzEybP36h1b+/8p78d0Imy5Qtq7x877+v/Hz1\nVcuO/+a/Ifc6depYNU/34sWLKF++PEQEb7zxhvb727dvY+3atZgzZw6GDh2Kp59+WrufepszZw4A\n4LXXXoOIYNSoUbmuxyOH7cMqbB8EALGxsejQoUO2EY6zWaYzbdu2DUWKFEFwcDAOHTqE0NBQeHh4\nQERgMBjg7++Pv//+O8eyo6KiUKJECW0EICwszPJzfGqqEnAEBADFiinrhi9dAtasAX791ZqnbpYr\nV67Azc0NIoIdO3YgMjLynpGb+7lw4QI6deqE5s2bY8eOHbCzs4O9vT3++usvm9dz0CClHb/8ss2L\n1pWuAYg6vzM+Pl7Ph33kTJ2qfPjeeAOoVEmZVXXpkgUFqCvXzYzszTV79myICBwdHfH+++8jNDQU\n/v7+KFu2bI69W0WLFsXTTz+NsWPHYsuWLTatS6HTs6fynnz7bX7XpNCJj1fyLLi4ACNGKC/jzJmW\nlZGWloZq1apZPOR+dxlPPfUURAStW7c268SVkJCA3bt3Y/ny5Thy5AgAYO/evRARFC9eHCkWrRh+\nhLF95BrbB2VmZiIiIkK7uC5VqhQiIyPvud/kyZO1NRtubm4YN24cjh49ipCQELi6ukJdPB4QEKCN\nbJw/fx7dunXTzuudO3fGqVOnrK90fLxyrQIAcXHAZ59ZX6YZ3nrrLS04y43Lly+jcePGEBGEhITY\nuHaK06eVASI7O+DQoTx5CF3oFoDEx8dDzYJA1lGX0nz8sfLT09PC2VRxcfipTRv4NmiQq/mdOYmK\nioLRaITBYMDSpUsBALNmzUKbNm0watQozJs3D5988gkmTpyozenNGoy8+OKLiI6OtkldCqNPBg3C\nnJYtcfCXX/K7KoXOrl370bjxGLz00g8YOvRvNG16FT/9ZFmGlMzMTERGRkJEUKVKFbN7vbIaMmQI\nRAQVK1bE5cuXLT4+qyZNmkBE8M0331hVzqOC7SP32D4eb/v378eTTz6pnWsDAgJw9erV+95/3759\nCAgIyBaIhISE4PDhwxgxYoQ2fc7Ozg6tWrXSgprixYvnOji9rzVrgAULlCwKFy7Ytuz7uHr1Ktzd\n3SEi6NKlC4YMGYKQkBBMmTIF4eHh+PzzzxEVFYVff/0V27Ztw759+3DixAncuHEDGRkZmDRpEtQF\n7HmZuW3YMOX6z6LplAWMbgHIjh07ICLw9vbW6yEfWZUrKx+8Tz5Rfj79tOVljBkzBiKCl156yer6\n7NixQ+sdmZmla+255567Z9TD3d1dC0o+/vhjfPPNN+jVqxdEBC1atLC6LoWV2jO4adOm/K5KobNs\n2TLtxKqmT1R7S80RHh6OLl26IDExUVs02KZNG0ybNg3fffcd9u3b99CeVvXizNnZGX/++ae1Twmf\nfPKJNvea2D6swfbxeEpLS0NYWJgWMJQtWxbfW7CGZ8+ePQgICMh27g4JCcGhQ4cQHBwMJycnLXuW\nv78/zp07l1dPJG/KfYAFCxbgpZdeynHmxsNuRYoUgcFgsDqL2MOcPQv4+QF5/DB5Ks8DkM2bNyMg\nIEBLudYqD/KJP06uXbuGunUHws9vJebMOYEaNZLw1luWr6U4d+4cHB0dYWdnZ9Ui8WPHjmlfQkOG\nDMn2t7i4OGzevBnh4eEIDAxEnTp1tF4V9Xbs2DEkJyfDw8MDBoMBJ0+ezHVdCjO1R89WI1KPk4kT\nJ0JE8M4772jzbs1d9Dd37lwYDAYYDAYsWbIEZcuW1RZn3n0rW7Ys/Pz8EBQUhLCwMERHRyM2Nhbb\ntm3TTvKLFi2yyXO6efOmlpGFSRzYPqzB9vF42rt3L1xcXGAwGDBs2DDcupW7tK1btmzB008/rb3P\nPj4+AIDVq1dDRFC7dm1bVrvA2LdvH5YuXYpPP/0UYWFheOeddzBixAgMGDAAAQEB8PPzw5NPPom6\ndeuicuXK8PLygtFoxPPPP49vdZoqunEjkJ6u7PFTGHdWyJMAJDMzEytXrkSLFi20D62zszOqVq2K\nUqVKYd++fXnxsI+FX375RZtD+7///S9Xw9Dq8Pcrr7wCEcHw4cNzVZcrV65oPWL+/v5mLU68cuUK\n1q1bhylTpqB3795a5g11g50ZM2bkqi6F2b59+7T51f8W9rQWD5CSkoLTp0/jgo2H0l988UWICKZN\nmwYRJU2nOb744iIcHJQT9MyZM1G7dm2ICOrWrYvw8HCMGTMGzz33HGrUqKHlx8/ppo7+jRgxwqbP\nS22f48ePt2m5hQ3bh3XYPh5PatD+mY3WTqiByMKFCwEAhw4dgoigZq7ScT66LN03xBrjxikz1BYv\nVrZl6NRJmRmj02w1q9k2AElNBSIjccvHB57/ZU4oUaIExowZg/DwcBQvXlz7QlqwYIFNH7qgiY2N\nzXWPw4OoJ5GRI0eiXLly2iiCuXbu3InixYtj1qxZ2L9/PwwGA5ydnREcHIyvv/4a//77r1mZK27f\nvq0FmM2aNUOilbur//DDDxARNG3a1KpyCpOTJ08iKCgIdnZ22pzTM2fO5He1bC4tLQ2RkZHw9vZG\no0aN4OjoiKCgICQkJFhdtslkQq1atSAiiImJwbp16/DDDz889LgFC5SFuS1bnsUHH4Rpa5IaNWqE\n69ev33P/9PR0xMbGIiYmBhEREQgODoafnx+qVq0KDw8P2Nvb48KFC0hKSsL69eutfl6AMnosIihT\npoyuJ7WCgu2D7eNBHvf28TBt27aFiGDDhg24ceMGYmNjbfI6qZ2GaurkcuXKWV0m5c6HHyp7Tk+c\nqGxOLf9lMDYaAR8fIDQUsODyUHe2CUBu3FBSM6k7popgdvfumDlzJqZNm6ZtSKP23Kv/fuaZZ2ze\n25OfMjMzER0djVatWsHJyQl169ZFQEAAbty4YbPH6NGjB0QE8+bNg4iSw9tk5gr0tWvXar1RPXr0\nwIQJEyCiZMS4e53Gg1LmZmRkaFkvvL29ccmiFFw5S0lJ0fKOWxJQFUanTp3CwIEDYW9vr40O+vr6\nwtPTE76+vo9ML29mZiaWLl2qjZKJCCpUqKA976pVq+IXKxYVx8TEwMfHB05OTnB1dUXFihURHh7+\n0PnoixcrX9AGg/K19cILJ+Ho6IbGjRvnqq126dJFG70rVaqUdrFlC+qF3+OUoIHtg+3DXI9j+zCX\nOkNi9erV2qZ6tmw7t27dgoiySJ3yx/PPKyMe7dopa0KWLlV+5+JyJxipWvVfNG3aFB988AGOFrAt\nMKwPQJKTgZIl7zzbJk2ATz9F0rhx2toAEUG9evXw5ZdfIjU1Fd999502j7RUqVJYs2aNDZ5K/klO\nTkZERITWyNXAQF33UqVKFWzevNnqxzGZTFowt3XrVkybNg0TJ04069ivvwaaNx8PEcGgQYMwcuRI\niCiZLEaPHo3Jkyfj+eefR8WKFXMcRs+6ePzZZ5+FiJL1wpLFjA/z8ssva1MFHklXrgAhIYj9r2fK\nwcEBffv2xaxZs1CuXDm4uLhoI4QzZswo1L16MTExaNSokfb5qVWrFiIjI5GRkYHNmzejadOm2t8C\nAgJwJeuutw+xfv16tGzZUju+TJkyqFy5svZ/b29vLFq06L7z3JcvVzZje/ddoG5d5WsrIOAUrl27\nlqvnqo7e1apVCz179rTpZ3jWrFkQURbb5mVGlQKB7YPtw0KPVfuwUJ8+fSAiWLp0KXx8fCAiNkkA\noMrMzNQyX1qzNwzlTnr6nX0b797rMTERiIoCXnwR6NTpg2zXcg0aNMCkSZMKxEae1gcgqanA0KGA\nry8QEQEEB2vhV//WreHr64vo6Oh7eunPnDmD9u3bQ0TZ4CYoKAi3b9+2ujp6unXrFsLDw7NtmFSl\nShWEhYXh1q1bOHHiBHx9fSGi5M4ODg7OVfpCdYi+Xr168PT0hKOjI3r27IkTJ06YdfzHHys9WkWK\nmPDee19oc2ednJzw3Xff3XP/hy0eL1u2LBwdHbF9+3aLn8uD/PjjjxARNGzY0Kbl5rurV4GxYwFX\nV+Xbwt4eU157DbNmzYK3t7f2ujZt2hSvvvqq9lo3aNDApicMPcTExGS7eKpUqRIiIiK0i8Xbt2/D\n29sbPXv2xJQpU7S9gby8vBAREfHAsnfs2AF/f3+t7BIlSiAsLAxJSUnIzMxEVFSUNt3Ew8MT9evf\nQEREzovztm0D6tdX3o5atSzcR+cu6enp2nTIOXPmaBd5ud1o88iRI+jWrRuuXLmCI0eOwMnJSRsd\n9PLygq+vr7bQNyoqCgcOHCjcFwBsH2wfFrhf+zAYDKhUqRL8/Pzw2muv4Y9Fi5TNdk+cAB7DTW+H\nDh0KEcEnn3yiXWvZavqbSp0amRfTzenBYmKU9lm37oPvl5ycjOjoaAQFBWUbFGjbtq0u9XwQ6wKQ\nlSuBFSuAWbOAF15QrnLVCWg9eyL1IScHk8mE8PBwLUNG7dq177vLpiotLQ2xsbGIjo5GaGiotqfE\n559/nuPc1Lxw4cIFhIaGahcF6kXz3LlzMWHCBHh5eWkLw9PT0xEaGgo7OzvtJGruqMGtW7cwc+bM\nbFPYSpQooY2sODs74+23335g41f3CjEYgClTgJ49U1Gt2v/g4eFhUZo4dfF4SEgIDAYDXFxccOvW\nLWzcuPGhJ0ZzpaWlaSNjhwrz7jpZXbgAeHjceRN69AA+/BAnu3fX3tPatWsjKipKC9J///13PPHE\nExAR2NvbIzg4OFfzwS9duoQffvhBl3nzW7du1U5y6shmWFjYPdM91q9fr/VklyxZEnPmzEHHjh21\n45599tl7NrHav39/trz0ajrInD73GRkZ+OqrrxAQEKENynp7K9NK1A7zpCSgYcM7X95WbkkAAHjn\nnXcgIhgwYIDW25ybk318fLw2reS1117TXpuiRYvC0dExx9FJ9bugfv366NmzJ8aNG4eDBw9a/6T0\nwPaR7X5sHw9mSfv4pE2bOzMznJyUiGrsWOCnn5Rhnhw24nuUjB07FiKCsLAwLSW+JWl4zaEGlll3\nUyd9DB2qfLTNnAQDAEhNTcXatWtRrVo1GAyGHDej1JN1Aci8ecrPL78EunVTGnlg4L3jQQ+xf/9+\n1K9fXxt2Dw0NRWZmJs6fP4/o6GiEhYUhMDAQPj4+2sX33TeDwYA33njDqqfzMHv37kVgYGC2jB++\nvr6YP38+Bg0apAVS6pdjVtu3b9cyubi4uCA8PPy+j3Pp0iWEhobCy8tLK++JJ55AeHg4kpOTcebM\nGfTt21c74ZQvXxWffZaKnGYknDkDVKkChIcDrVopH9h69dKwe/eDA70HUdfxLFiwAHZ2dnB2drZZ\nD8jAgQMhIpg8ebJNyisQ/P2VhN0ffgg0bqxdbL3t56dNu7hbUlISQkJCtMD1YfPB09PTceDAAURE\nRNwzatWlSxebn3iyWrJkifY5LVmyJGbPnv3A6RCxsbH3XFTNnz8fJUqUgDrFJiwsDAcOHEBgYCCM\nRiNElPzqISEhiIuLe2idMjKAZcuU3lv1GiRrBvDp04EaNWyXLSQ2NhZGoxFFihTBuHHjICLobeEO\nUSaTScu7X6tWLW2aZOnSpbUT/Pnz5+9Z6Ovt7a29Rurt5ZdftnlvZ55h+8iG7SNnD2sfGRkZOH78\nONauXYs5c+bg8DvvKJtklS9/50m++CIwf75S4COeCGfy5MkQUbKE9evXDyK538X+ftQRtUemw7CQ\nyMwEypZVPtJ791p2bEpKCjz+SxJl7iyavGJdAPLdd8q3WEQEcPKkVV0lt2/fxrBhw7STwv3S+tnZ\n2aFGjRro0aMHJk6ciKioKHz66acQEXTs2NGqp/MgV65cQWBgIESU6VT+/v5YsmQJAgMDtZOg+vtt\n27blWMatW7cQFBSkPZfnn38+246kx44dQ3BwsNYDpgY4OU1hA4Bdu3ahbdu2aNv2B8h/Q+VRUfc+\n7okTd3q0KlcGrF2HNnfuXIgIXnjhBS3TxldffWVdof/5+eeftRNMXjh06BAmTpyIoUOH5kn52ag9\nbQsWKFMU1ZNghQrK78xI3P3PP//cMx/86tWrOHPmDKKiojB69Gi0atUqx8Dc3d1dS5354osv5slT\nPHXqFBo0aKDtlnvz5k2zj42Kisp2UTV+/PhsG1+pF1YuLi548803H7h77/1kZiptolYtZTPdN95Q\n0hV+9ZUyT9aW1Fz5H3zwAezs7ODo6GhRndULBi8vL8yfPx8GgwEODg5mjVQmJCRg9+7dWL58OTp1\n6gQRwXPPPWfN08l7bB8PxPaRnTXtAwkJwO7dwP79yqyNx2AEJDw8HCKC4OBgDBs2DCKCjz/+2KaP\n0axZM4gIduzYYdNy6cE2bVK+KqtXt/zY6OhoiBSMjKPWrwGx8ULANWvWoHTp0jAYDChdujT8/PwQ\nHByMiIgIbN68Ocd1IidPnoSIsjYhr6jpb318fHDy5MlsPVQuLi4YOnSo2dmboqKitKlGpUuX1tZa\nmBvI3G3VqjRUr37n/P3MM8r3LKAMp1erpvy+Th0lU4K1Ll68qI18zJ492yYXO+qJOT09XcvKFR0d\nnet5wlmdPHkS4eHh2noc9ZbnF2hZe9p69VJWjIWFKW+KBdLS0jBt2jTtIqpMmTL3XEwZjUbUqVMH\nAwcOxGeffYZ9+/YhIyMDBw4cgIigcuXKtn9+UHbLFVHm4+fGpUuX0KtXL+15jBw5EtHR0fDw8ICn\npyf69etnkx12MzKU5WqLFgFz5wJffGF1kff46quvIKKkpVYz/zxopDOrX375BXZ2djAajfj444+1\nuf/z1FFmC1y7dg1OTk6wt7fH+fPnLT5eN2wfD8X2obBl+wCgBCGdOwMFuX1YadGiRRAR9O/fHyEh\nIVrwZ0sdOnSAiJLimfQzevRYNG3aGWFhOy0+Vl0DbOvPQm7k+U7ouaH2aG3ZssWs+5tMJm1IKTe9\nQOZQh3tnz54NABgxYgQ8PDwQHBycq5P86dOn0aZNG4go81dFlEXhgYGBOGzhFDZAOXnMng14eSnB\nRrduyu8zM5U0bc2aKWs9baVdu3YQURa42dvbw8nJyaLevaxiY2NRpkwZbepdrVq1UKlSJYgoKf7u\nTglszoLb06dPa0FH1gX0Xl5e2v4llk4BsFjWnraLF5VeOCscPXoU9erVg7u7O5ycnODr64uQkBBE\nR0ffN0ONyWTSPl95cTG6ZcsWiAhatWqlrcfas2ePxeX89NNPqFmzphbEq3sG2bo9L16s7Bo7erRN\niwWgDG2rPdazZ8+Gvb09Ro0a9dDjjh49qr1H7777rjZVMzAwMNd1UXvKp06dmusy8hzbh9nYPmzb\nPhAQoJwoC2D72Lx5c64+I3dbsWIFRJSU+1OmTIGIYNy4cTao4R3PP/88RASrVq2yabl0fyaTSbs+\nsjQJR9Z1tgUhnXmBDED69+8PEbFos8Inn3wSIoJNmzblSZ3U3WTVqUbXr1+3evO9jIwMLT3ec889\np+1Qbo1r15REZL17K0txPv9c+b0Zsxks8sknn0BE8L///U/rBcnNgqarV69qC0q7dOmCN954Q5ty\ncL+UwG5ubmjdujVGjhyJJUuW4ODBg8jIyMCZM2dyDDqKFi2KwMBAREdHIzU1FUOGDIGIkpElz9l4\nhPD777+HiJJ20lzqaF1enCR++eUXqNMfGzZsCBHBP//8k6uyso52qb3ZhS21ptpR8dprr5l1QZuR\nkaG9bi+88IKW4rpx48ZWPXf1falatapNRhHzDNuH2dg+bNc+8Msv6iYJBSpD1sqVK2E0GtGoUSOr\n00yr05k7duyIjz76CCKC4cOHW13HrMG8mjr/yy+/tLpcMs+OHTsgIqhYsaLZe8Cp1M9E/fr186h2\nljFKAVS3bl0RETl48KDZx/j5fSrNml2Wo0eb5kmdLl++LCIiZcqUERGRYsWKSZEiRXJV1r///iuz\nZs2SDRs2SLFixURE5LXXXpNSpUpZXc/ixUU++kjk6adFfv5ZmZQlIuLgYHXR2fTs2VPs7e3ll19+\nEX9/fxER+fbbby0qIzk5Wbp16ybHjh2TJk2aSMeOHSU8PFwcHR0lOjpazpw5I1evXpWff/5Zpk2b\nJj169JAqVapIYmKibNmyRT766CN5+eWXpW7duuLh4SGVK1eWN954Q7Zu3Sru7u7Sr18/iY6OlkuX\nLsmSJUvkf//7nzg6Osru3btFRMTHx8e2L0pO7O1tWlyrVq1ERGTnzp2SmZlp1jHt2z8nvr795MiR\n0jati4hIYmKiiIi4ublp/85tuzAala+jzMxMSUlJETs7O3FxcbFNRXUyePBgERGJjIyUHTt2yPHj\nxyUjI+O+97ezs5NZs2ZJmzZtpEqVKrJu3TopXry4rFy50qrn3rFjR6lWrZqcPHlSfv/991yXk+fY\nPszG9mG79iEdO4pUqyZy8qRIAWofnTt3Fm9vb9mzZ4/MmDHDqrI8PDxERCQhIUHKlSsn9erVk+LF\ni+e6vJSUFJk0aZJUrVpVDh8+LCIi7u7u2mOQPlauXCkiyjWYwWDI1bE9evSweb1yJb8joJysW7cO\nIoL27dubfcycOUqHxl3Jp2xGzfZgi/SWalaUvn37atm/bDHkmtXChcpiwkGDbFpsNn5+fto8XAcH\nBzg4OJidCjkzMxPd/0u1WaVKFURGRsLOzg4Gg+GhmTpu3rx5zz4lZcuWhYigWLFimD9//n1Hp1JT\nU+Hk5ASj0Yj4+HiLn3NBUKVKFYgI9u3bZ9b9f/pJaRtt2ti+LpGRkdp0CHXuvbW7HN+8eRMiAg8P\nDxvVUl99+/ZF3bp1tRE4BwcHeHt7Z1vPFhMTg9jYWO2YVatWwWAwwN7eHr///rtN6qEu2n3ppZds\nUl5hwfZRsBWU9oHJkwERmPr2tU15NvL777/DYDDAycnJquuN/fv3Q0RQp04dq+u0adMmbaaCnZ0d\nIiIicOLECTzxxBOoVq0ali1bZvVjkHnUKYjmLlFQZWRkaOtr96uLhPNZgQxAzp49CxFlzwtz/fpr\n3p1EAGhzUHO7G2xWM2bMgIjgzTff1DaGuXjxog1qqa/PPvsMIkqayGeeeQZGo9HsXe3VoXhPT098\n/fXXcHV1hYhg+vTpuarL7t27tfzy1apVu+9eK7t374ZI3mXZ0kPv3r0hIvjss8/Muv/168oWCy4u\neTcVb+jQoXBzc4OIWB3YnTt3DiJ5m1Qir02bNg0dO3ZE5cqVs00HvPvm5eWFxo0ba1NqzF2Ua45z\n585p2YYs2Um7sGP7KPgKQvswnT2L5W3bolrJkjY5r9vS4MGDISJo0aJFrjcZXbt2LUSUbG+53TT4\n9u3bCAkJ0TKu1a1bF9u2bcPs2bO1c3aTJk0sngpEuXPmzBkUKVIEjo6OmDt3rkXr1n7//XeICGrU\nqJGHNbRMgQxAgDsX/Oaui7hwQQlAvLxsX5eUlBQYDAY4OjrapKG9+eabEFE2CDIajbCzsyuUOxnf\nuHEDjo6OcHBwwNatW83u2VMzZzk6OmLJkiVaVB4UFGRVfS5evKitqSlevHiOPQRq0NS3gPV6WUJN\nrzhw4ECzj6lRQ2kfu3fbti5qMD1mzBgYjUYYDAarP8tHjhyBiLL3zaMgNTVV2zw1LCwMQUFB8PPz\nQ+nSpbULrYCAAIveT3N17doVIoIPP/zQ5mUXVGwfhQvbx71u3rypbUA8d+5ci45NTk7GuHHjYG9v\nn21/ID8/P4sWLa9bt05b7Ozg4ICQkBD8/fff2npb9X15nDo3CgI1oYB6q1OnDkJDQ3H06NEHHjd8\n+HCI2D4RgTUKbACipky1ZDOtEiWUk4itk5mcOnUKIoIKFSrYpLy+fftCRLSFYaVLl7ZJufmhc+fO\nEBFMmjTJrEVzN2/eRKlSpWA0GrFw4UJtWLdr165WL7oDgMTERO2k4urqiujo6Gx/13UBeh7ZuXOn\nxUPrr7yitA0bp4FHaGio9qWmvubWUkepGjdubIMaFmyXL1/Gpk2b8mw3bnVRdu3atfOk/IKI7ePR\nkdftY9WqVQW2ffz4448QUTaYzDod7UG2bt2q7W1jMBgwYMAATJ48GZ6entkCkV27dt23jLi4uGz7\nlTVq1Ag7d+5EWFiYttly2bJl83TzTrq/27dvIzo6GoGBgXB3d88xGLl7Y0iTyaQl9dlt614WKxTY\nAES9ULQk+h82DBg8GDh1yrZ1UbMO2GrjFjVr1Pz58yEiaNiwoU3KzQ9ffPFFtl3s69Spg8DAQISH\nh2Pz5s1ISUm555ijR4/ik08+QcuWLbXX1dqMYlllZGRg6NCh2nzVrJsvqSme8ypbmh7S0tLg4uIC\ng8GAGzdumHXM558DTz4JLF1q27qoo3kTJ06EiLLTs7X++OMPiAhat25tgxo+3tLT01GuXDmIiNn7\nChV2bB9krqztI7fTlPKSmn3z6aeffuDsC3WqlLqXWPXq1bNt0Hjt2jWEhoZq2xUYDAb4+/vj77//\nvqesqVOnasHyzJkzsXv3bjRp0kQ7LigoCLdu3cqT50uWSUpK0oKRrEGmGoyEhIRg8+bN2LZtG0SU\n9bYFabpcgQ1A1N22hwwZ8tD7JicrG+8NG5Y99aytrF69GiICf39/m5SnLsBTA5BnnnnGJuXml3Xr\n1mkLo+6+ubi4oEWLFnj99dexePFi7N27FykpKVr+8KpVq+LSpUt5Uq/w8HBtfnFwcDCSk5O1BeiF\n/Qu0VatWEBH88ssvD71vZiZw+3b2XY5tRQ30wsLC4OfnhxdeeMHqMtW5y507d7ZBDWnIkCEwGAwo\nUaIE/P398eabbyIiIgK///57oVx7Zg62DzJXQW4fV69e1aYof3Gf3SE3btyI6tWrQ0Rgb2+PkJCQ\nHDv+1PJCQkK09RtqIJI1CU5qaioGDx6M/fv3ZwtqvL298dtvv+XF0yQbSElJwZo1a9C/f39tr4+s\na6lElHXHBYltcyDakJqK98CBA9rv0tPT5dixY3LsmLPs3u0tBw+K7N8vcuKECCAyb56SetbNTeTa\nNSUlrYVZynJ06dIlEbmTgjc1NVWSkpLEy8vLqvLS09OzlVtYde7cWY4fPy43b96U3bt3a7e///5b\nYmNjZceOHbJjxw7t/o6OjpKWliYlSpSQn3/+WUqXtn36SxGRkSNHipeXl7z66qsyd+5cOXTokKSm\npkqtWrW0FIWFVYsWLWTbtm2yfft26dSpU7a/xceL7NsnsnWryJYtItu2iYwYIVK/vsiqVSIuLiLH\nj4tUr259Pc6cOSMiIkWLFpWYmJhclQFAVqxYIXFxcTJkyBC5ffu2iOQ+XSndkZCQIL///rvY29vL\ntWvX5Mcff5Qff/wx2308PDykRo0aUqNGDalVq5b274YNG2qpXwsbtg8yR0FvHyVKlJDZs2dLYGCg\nvPHGG9KxY0cpX768iIjcvHlTQkJCZOHChQJAGjZsKIsWLXpgevkSJUpIWFiYvPnmmzJ79myZO3eu\n/Pjjj7J27Vrp0aOHTJkyRWrUqCGvvPKK9OzZU44cOSL29vYSHBws06ZN42euAHNychJ/f3/x9/eX\nzMxM2b59u6xYsUJWrFghFy9elAoVKoinp2d+VzO7fA6A7uvq1asQETg7O+Oll15C/fr14ejoCBFB\n27YjoYQcys3BAahTB5g2TenB6t8faNAA6NRJWZxuLbUXS921dfTo0ShTpgzWrVtncVlpaWkwGo2w\nt7fXhjrfeust6ytZQN26deuelLlGoxETJkzA1q1bdanDr7/+qg09iwj69Omjy+PmpaioKK0X9MCB\nA4iMjERQUBDq1KmDOnV2ZWsfIkC/fnd2OW7VCnB2BkJDgdTU3D3+/v37ERAQAIPBgAoVKqBEiRKI\nioqyuJyYmBhteN/NzQ2XL1/G9OnTISJ45ZVXclc5AqDM+1V3RK9ZsyZ27NiBlStX4oMPPsCAAQPQ\nsmVLbUftu2+urq4FewPDh2D7oIcpTO2jW7duEFE2/gWA6OholC9fXrtGCg0NRVouUridPXsWr732\nmnZtZW9vj/r162szBxo0aIC//vrL1k+HdJSRkYF3331Xm5pni7W2tlJgAxAA6NOnj5a6UERgNBpR\nrVo1DBkyEePGAcuXA/v23XuS2L0bKFZMObGULg38+GPuHn/z5s3w9/eHwWBApUqV4OXlhR9++AGt\nW7fW6jN27FikWnCWSktLw9dff40FCxZoqWhnz56duwoWUvkx/Wn//v0oU6YMOnfujB9++EH3x7e1\nM2fOaEPod58c/fzm4qmngLFjgVWr7g3ChwxR0o6KAA0bAhYkRsGhQ4e0Cyv1osjb21t77O7du5uV\nGnDz5s1o06aNdlyFChUwadIkvP7663B0dMSTTz6Jn3/+2bIXhbJ5//33teH3Y8eO3fd+165dw9at\nW7F48WK8/fbb6N69O7p3765jTW2P7YMepjC1j/Pnz2uZQdW1kyICX19fHD582Oryz5w5g+DgYDg5\nOaFFixZa1itLrm2o4MrIyND2slu8eHF+V0dToAMQAIiIiMAXX3yBv/76C7dv3zb7uHPngA4dlJOI\nwaBsUJiU9PDGlJ6ejq+//hqNGjXSGrm7u7u2xsFgMCA4OBhTpkyBg4MDRAT16tXD3r17LX5uHTt2\nhIjgK1tOOqb7unnzZn5XwabWr1+PypUro2zZsggICHjgwv+7xcQA1aop7cPODpg6decDEwGcOnUK\nQUFBsLe3h4iSQjkoKAiXLl2CyWRCRESENsrk6emJ8PDwHHsIt2/friVhEFH2+nn33XcxfPhwLcOK\nvb09hgwZUqAWyxU2v/zyC+zs7GA0GvFjbntgCjm2D7qfwtg+Xn31VRgMBhgMBri4uCAsLMzmozCx\nsbE4e/bsfffRosLrq6++goigcuXKBSawLPABiDVMJiA8HHByAnx9l6BWrVr3TUGWmJiIiIgILS2s\niKBUqVIIDQ3F9evXYTKZEB4erg1V1qlTB8uWLdPu7+zsjPDw8IeeFDIzMxEdHa31YvTu3RsHDhzI\ni6dPjwFLgvK7JSUp00xq106Cg4MrqlSpck+v6tmzZ7WeMREl01lQUFCOvbgXLlzQdrdXe+eypgO8\ncOGCFrQXK1YMEyZMwLhx47QLM6PRiICAAJ78rHTixAlt6sjUqVPzuzr5iu2D7lbY2sepU6fQqVOn\nbKN4PXv2zO9qUSGTmZmJBg0aQESwYMGC/K4OgEc8AFH9848JzZu31HqmZsyYofUcXLlyBaGhodnm\nelarVg3h4eFISkq6p6x9+/ZpaWednJwwefJkDBgwQDu2S5cuOU4xSk5OxoIFC7IFOCVKlMDXX3+d\n58+f6EF2796Pxo0ba5/LwMBAHDlyBCEhIXBxccl28XP8+PGHlpd1fvLdQ/ljxozB22+/jUmTJmlT\nCpRpMX45poQkyyQkJKBevXoQEXTr1o295DbA9vHoKEztQx05U/d68PLywtSpU1GkSBGISK7WoNLj\nTV0bV7FiRbNGgvPaYxGAAEq+5OHDh2tzc1u1aoXAwEDtBKL+7vvvv3/osGZSUhKCg4O1sjp06ICF\nCxeiePHiaN26dbadbm/duoXw8HAt17iIkno2PDzcqt45IltKT0/H9OnTtfagnuSMRiNeeukl/Pvv\nvxaVp25mpbaR+vXr448//kBERES2HY4t3Z2X7s9kMqFXr17aotrCnmq6IGH7KPwKU/s4fvw42rdv\nr30O/P39tVE1NQlBpUqVEB8fn881pcLEZDJpywvmzZuX39V5fAIQ1a+//opy5cppU6nUPNgxMTEW\nl/Xzzz+jbNmyEBEULVoU8+fP13ZsPXXqFEJCQrJtDtOoUSNERkYWqCwERFkdP34czzzzDN599108\n99xzuVrblNX69eu1HPXqNBURwVNPPYU//vjDRrUm4M4GYu7u7vfshEu2wfZReBWG9pGeno6wsDA4\nOztDRFCmTBl899132e6TkZGBZs2aQUQwYsSIfKopFVY//PADRJTd7PO7E/yxC0AA4PLly/jpp58w\nY8YMq7+ILl++DH9/f+3E8b///Q8vvviithjRYDDgmWeewfr1621Ue6LCJSkpCSEhIZg8eTJatWrF\n7D154Ndff9UW1a5Zsya/q0MWYPvIe4WhfezduxdNmzbVriUCAgJw/fr1HO+7b98+ODo6wmg0MlAl\nizVv3hwigjlz5uRrPQwAIGQVALJw4UIZNWqUODo6ys2bN8VoNEqXLl1k4sSJ0qxZs/yuIhE9wv75\n5x954YUXZNCgQTJhwoT8rg5RgVKQ20dKSoqEhYXJBx98IGlpaVK1alWJiIiQjh07PvC4iRMnyvvv\nvy81atSQvXv3irOzs041psJu7dq10rVrVylRooScPHlS3Nzc8qUeDEBs6ODBg3LlyhX57bffJCgo\nSCpXrpzfVSKix8T169elWLFiYjAY8rsqRAVOQWwfhw8flhdeeEGOHDkidnZ2MmLECJk6daq4uro+\n9Ni0tDRp0qSJHDx4UMaNGydTp07Vocb0qGjZsqXs2LFDQkNDpWnTpuLm5ibu7u5StGhR7d/mfA6t\nwQCEiIiISGe3bt2SunXriqurqyxcuFDatm1r0fE7d+4UX19fMRgMsmPHDvHx8cmjmtKjJiYmRjp1\n6iQeHh4SHx9/z9+reXnJ8Vu3RDw9RTw8RNzdRdzclJ/Ll4sUL251HRiAEBEREeWDgwcPSvXq1cXJ\nySlXx48aNUrCw8OlYcOG8tdff4mDg4ONa0iPqubNm8vu3bulUaNG4ubmJomJiXLz5k1JSEiQeu7u\nsuHEiZwPvHZNZNs2kfLlRZo0yfXjMwAhIiIiKoSSkpKkQYMGEhsbK9OmTZN33nknv6tEhUBmZqY0\nbNhQTp8+La+88oo0aNBAPDw8tOlXZTw8pKabmwggkpkpkpEhkpAgkpgo8tdfIoMHKz+dnEQesl7p\nfhiAEBERERVSv//+u3To0EEcHR3l77//ljp16uR3laiAmzlzprz11ltSrFgxuXHjxj1/H1K/vizY\nv//OL+zslKlYpUuLvPKKyNtvi9y8KbJihcirr+aqDva5rDsRERER5bP27dvLwIEDZcOGDTnO5yfK\n6uTJk/Lee++JiMjYsWPl+vXrEh8fLzdv3pTExERJSEiQ6iVLity+LRIfr4x6pKSIxMWJODiI1Ksn\nEhGh/L1371zXgyMgRERERIVYfHy82NnZSZEiRfK7KlSAAZBOnTrJb7/9Ji+//LJERkaad6A6BSsh\nQaRSJZH0dCUYsQIDECIiIiKiR9ylpUvl2bFj5aIoCRCK2yCbVW4xACEiIiIiepSdPy9St67AYJC9\nX34pjbp1y9fqGPP10YmIiIiIKG8NHy5y65YYWrfO9+BDhCMgRERERESPrq+/FunbV9lY8OBBZQ+P\nfMYRECIiIiKiR9H16yKjRin/njOnQAQfIgxAiIiIiIgeTcePK/t4dOggMmBAftdGw31AiIiIiIge\nBSaTyE8/KSMdTZqIODuLHDoksm+fSFqasnt5AcARECIiIiKiR8H06SLNm4tcuCASEyPyyy8iRYuK\nHDkikpSU37XTMAAhIiIiInoUGAwipUuLtG4tcuqUsoN5RITIH3/kd82y4RQsIiIiIqJHQb16SsBx\n+7ZI794iS5eKDBkiYixYYw5Mw0tERERE9KhITxdxcMjvWjwQAxAiIiIiItJNwRqPISIiIiKiRxoD\nECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIi\nIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi\n0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0D\nECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIi\nIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi\n0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0D\nECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIi\nIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi\n0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0D\nECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIi\nIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi\n0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0D\nECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIi\nIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi\n0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0D\nECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIi\nIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi\n0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0D\nECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIi\nIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi\n0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0D\nECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECIiIiIi0g0DECJ6/GRm5ncNiIiI\nHlv2+V0BIqK8kpIicuxwhtTf95XIgQMi+/aJHDwo8tJLIu3biyQkKHfq3z+/q0pERPTYYABCRI+U\nNWtEli5VYo3jx0VgspMMt2AxJMTfudORIyLe3iKvvSYSESECiBgM+Vdpur+YGJG4OBGjUSQ9XaR3\nb5Hly5WfRI+7detE4uNF0tJETCaRV14RiYxUfhIVYPpMwfrjD5HvvhP55ReRH35Qfqf+JHqc/fST\nyDffiHz5pfL/xYuz/ySzjRwpcvasyLffiqxYocQYBoNIrdoGuRE4UmTqVOV7JzZWJDpapFQpkblz\nRZYsUb6f8khKSoosX75c/vrrL4mLi8uzx3lkxcaK9OolcuOGyLFjSns5cya/a/VoAjg9sbA5dUrk\nxRdFkpKUfy9frnwREhVw+oyA/PuvSFCQyOefi1y6pPRinT+vy0M/StLTRU6fVs7HLc6vFM9D20VO\nnPjvFy1E+vYVuXJFxN1d5Jln8ru6ZI7Tp7P3wp89q5xAzp3L75oVOg0bKvFF3brKCEi9eiJ16og4\nOoqITL73gJ49RT77TGTbNpGrV0Wef17EwcGmdYqIiJDly5fLpk2btN95eXmJt7f3fW90lxIlRKKi\nlO81g0GZPrdwYX7XqtC6eVPk5Enl5nH5mPgdCL/zi1OnRD74QKRcOWXEqWfPfK4tPVTJkkpQHhcn\n4uysBCOff57ftSr01Hju1CmRFpd+kGJHtiv/OX1a+R6aOPHOyFNgYD7XtpCCHqKigE8/BWbPBubN\nU363aJEuD12YTZoExMcDX38NVK0K2NsDylUqcKZFwJ3/iAAtWwIREcqBCxfmb8XJfCtWAMuXA19+\nqfx/wYLsP8ksJ08C//uf8tF/4w0LDszIAGrXVtqQ2n5sZOfOnXBycoKIoGXLlmjcuDHc3d0hIve9\neXl5oVmzZjh//rxN61LoxcYCrVsDw4bld00KnbfeAq5cUb5qunTJftroXevv7L8QAQYNUg60cXug\nPJSeDpQurbx/Fy7kd20KpTfeAM6fB95/HyhZMnuTOOHbL/svHB2B+fOVA3muzjUDAOgS6YwZIzJ7\ntsiMGSJjx+rykIXd4sXKlE4nJ2UAKS1NpGJFkWrVRGa2XCVN3I4q89irVVN+/vabyLVrIqtXK691\n3br5/RTIHBkZIp9+qkwLGjVKpEuX/K5RobNypdJZ261bLmZ3rlihTPEpV06Z4uPqanV9rl+/Lk2b\nNpVTp07JiBEjZO7cufLJJ5/I4cOHpXz58uL632MkJibK6dOn5eDBg3LgwAG5deuWGI1GKV++vPzw\nww/SpEkTq+vySNi+XaRVK5EnnxTZsSO/a1OofPGFSHKySJEiIps2KYNJ3t4iVauKNKuVIO+WW6z8\nR/3lunXKNKzt25XhEnV6KBVc6oVCZqZIaqrNR3IfB4sXK6MeSUkiISHKYFLlyspt2pOrxcf5oEiV\nKsovqlYV2bpVuSjbulWZcdKtW34/hUJHv0XoN24oP4sX1+0hCzujUaR+fZH160UOHRIpX16dTiIi\n0v3eAwICRIYPV9balCghsmyZTetjMpnk3Llzsm/fPmnQoIFUqlTJpuU/tuztRfbsUQJILqzNlePH\nlZ/VqokMHqxcS73zjpnrynv2FGnaVOTAATkSGSk1hw2zqi4mk0n69Okjp06dkhYtWsisWbNERGT1\n6tXy66+/Zruvvb29VKpUSapXry79+/eX8ePHyxtvvCFff/21fPbZZ7JgwQKr6vLIuHZN+VmyZP7W\noxBKTVVit+++E/nkk7uXl7mLyMjsB/TsKZKYqHSGXLwo0qePSKdONq9Xenq6HDhwQOrWrSuOd05s\nlBvXrysdWcWLM/jIJUCkY0dlVu6FCyJly2b9a7f/blkEBCidV59+qpy7/f1F7Ox0rPEjQLexluee\nU4aufvhBt4cs7OrVA7y9ldkHZjtzRhketLMDjh61WV0yMzPx/PPPw8/PT5su4unpCV9fXwQFBSE8\nPBybN29GYmKizR7zsfLSS0r7+Prr/K5JoTR4sPLyTZ2q/Cxb1rLjr23YgBblyqFYsWKIi4uzqi5v\nvfUWRASlS5fGuXPntN9v2LABc+bMwWuvvYZOnTrB29sbdnZ2Wnt64oknAAD//vsvDAYD3NzccOvW\nLavq8shYtEh5YwcMyO+aFDqjRgFOTrmYmRsWprzm7drZvE6LFi3Cyy+/DBGBvb096tSpg4CAAISG\nhiI6OhqXL1+2+WM+0g4cUN6r2rXzuyaF0pUrgLMz0KGDhQdmZABPPKG89lFRNq/X4cOH8c0338Db\n2xs9e/bE1KlTsXbtWly8eNHmj5Uf9AtAfH2VN2nzZt0esjBLSVHWfNjZAadPK59zsw0YoLzWr75q\ns/q8+eabEBG4u7ujdevWKFGiRI5z2O3s7FCrVi306tULU6dOxdmzZ21Wh0dat27Ke/b99/ldk0Kp\nXTvl5Zs5U/nZpo3lZXTo0AEigvHjx+e6Hj/88AMMBgPs7e2xcePGh94/NTUVR44cwU8//YRVq1Zp\nv2/bti1EBBGch6+YPl15Y8eOze+aFDpq24iKAhISLDgwPh7w8lIO3rLFZvVZv349HBwcICKoUKEC\njEZjjueSSpUq4bnnnsOECROwcuVKnDp1ymZ1eNRkbNyY+y8+ws8/Ky9f27YWXmsByvpmEaBRI8Bk\nslmdLl++DG9vb9jb2993vaCvry+Cg4MRGRmJAwcOwGTDx9eDbgHI1a5dkVy9OtIPHtTrIQu13bvv\ndGj4+ACursA//5h58NGjgJ0dTnh745INAoBFixZBRODg4ID169cjLS0Nzz77LIKDgzFhwgSEhISg\nb9++8PHxgbOzc7ZGMnLkSBw5csTqOjzyOnVS3vCff87vmhRKzz57GjVqJGPaNBNEgIEDLS/jzz//\nhMFgQJEiRXLVw3TkyBF4enpCRPDhhx9aXoEsvv76a4gIGjZsaFU5j4pVU6diYevW2PXJJ/ldlULF\nZLoTQ7z9NmA0Au+9Z0EB774LiOC4jUaeDh8+DC8vL4gIQkJCAADz589HWFgYpk6diuHDh8PX1xdF\nihS554LLz88PO3bssEk9HjXffvstXB0cMOLll/O7KoXStGlKGxk5EihaVLnmSk018+DkZC0BQMZv\nv9mkPklJSWjZsiVEBD4+Pvjrr7/w0UcfYeTIkWjXrp12nskpKHn22WcLTSCiWwBSsmRJiAiHVs30\nxRdKg3jxRWVo0GAALJmN8eno0bC3t0dwcLBV9di4cSMcHR2z9cbu27fvng++o6MjGjVqhH79+mHU\nqFFYsGABGjRoYHWP8uMiZuBAbGvXDhe3bs3vqhQ6SUlJMBqNcHJywpgxY1CmjA/Cw1eafXxCQgLG\njx+PlJQUvPDCCxARPPPMM9iwYQPOnDlj1pd5QkIC6tatCxHBiy++aM3TAaCMjJQqVQoigj///NPq\n8gq7/v37Q0SwePHi/K5KoXLypHIeKVkS6NVL+beacM8sV6+iR+vWMBgM2Ldvn1V1uXr1KqpXrw4R\nQY8ePZCZmQkA2udc7eSqU6cO+vXrh3HjxmH8+PEYOXIk2rRpA3t7e5QrV047ju6YO3cuRASvv/56\nflelUAr4L6moGoh4e1t2/PmPPkKfVq3Q+ZlnrK5LZmYmunfvDhFBlSpVcOnSJQBAgwYNtPYRGBiI\n9957D2FhYRg3bhz8/f1RtmxZiAjKlCmDRo0aIS0tzeq65DVdAhCTyQR7e3sYDAakp6fr8ZCF3jvv\npGoReW4axMGDB2E0GuHs7IwLuUzLl7W36u2339Z+f/PmTXz//fd477330KNHD1SvXj3bMLqdnR2S\nkpKwceNGiAgqV65caCLy/NK0aVOICP7666/8rkqhowbEtWrV0gKIKDPn4yYlJaF9+/YQEQwcOBBz\n5syBi4tLtmFvR0dHeHt7w8/PT1vvFBMTg9jYWGRkZMBkMqFXr15aHWy1bmPMmDEQEQxS06I+xrp2\n7QoRQXR0dH5XpVCJjk6HCNCxI1CjhnIu2bPHsjJGjhwJEUGfPn1yXY/k5GS0atUKIoKmTZtqawXT\n09MxefJkdO/eHd7e3jAYDPd0boWGhgIAatSoARHB+vXrc12PR9WECRMgIpg0aVJ+V6VQql49ewDS\ns6dlx9+6dQtFixaFiGD79u1W1WX06NHaGtv9+/cDUIKS2rVr59g+7O3tUb9+fcyYMQPnz5/Xgvyf\nfvrJqnroQZcA5Ny5cxARFC1aVI+HeyR06NABRYtWxbvvfoG2bVdi+PB/LS6jW7du2Ya6LXHt2jU8\n8cQTEBF07979ob1OCQkJ2L59OyIiIrQvQZPJhCpVqkBEsDkP1v7cvn3b5mXmlzp16kBEcODAgfyu\nSqGzatUqiAj+97//oV69ehAR/GPGfMWUFODZZ/0hIihfvjzmzZunBR4tW7bEU089hTJlyuQ41K3e\nnJycUKFCBYgIPDw8bDrd8Pjx4zAYDHB1dbV6YXxhtnXrVjRq1MgmJ/fHzXvvvQcvr2p4993FaNfu\nB7Rpc8L8qSX/OXv2LBwdHWFnZ4ejuUhsYjKZ0KdPn3t6dHNy69Yt/PHHH5g7dy4GDRoEHx8frFmz\nBgAwceJEiAgGDx5scR0eVbdv38a8efO0NWOfcIqixW7duoVq1bqiXbsIvP/+XjRufB3Tp6dYXE5I\nSAhEBD0tjV6yiIiI0EYCf8thOld8fDw2b96M8PBwBAUFwdfXV9trSr3OmzJlCkQE/fr1y3U99JKn\nAcjly5cRGhoKDw8PiAiqV6+elw/3SFGnrL3++uvZeoHMERoailmzZmH37t0wGAxwcXHBvHnzsHv3\nbqSkPLxhpaWlab3CTZo0sSqz1dtvvw0RwdChQ3Ndxt1MJhPeffddBAYGwsPDAz4+PggMDNQycSVY\ntNIy/8XHx2uB2okTJ/K7OnnmyJEj+Oyzz7B06dIHXoRYasaMGRARjBo1Ci+99BLq1KmD+Pj4Bx6T\nlqYk5mvR4izKlauMiIgI7Yv87l7ElJQUxMbGIjo6GmFhYQgKCoKfnx+8vb21kT8vLy8MHz4cgDLV\n5GGPby51Yfx8ddOrx8jmzZvh7++vBXcigmPHjuV3tQoVdUQwNDQUIoJGjRqZfWxycjLefPNNxMXF\nYeDAgRARPP/889kyu5njnXfe0d5Da6ZxHT16VOvINOc89ii7desWwsPDtWk36nfXnDlz8rtqhY46\nU6N58+Zo164dRARr1641+/iEhARs27YNly5dgrOzM4xGIw4dOmRxPX766SdtplBkZKTZxyUnJ2Pn\nzp3ad2NsbCwMBgPc3d2RlJRkcT30lCcByKlTp/D6669nW5BcqVIleHl5PXLDpxkZGdi6dStOnjxp\ns+ll58+fh4igWLFi6NKlC0QE33333UOPM5lMCA4O1qaNrF69Gk5OTqhUqVK24Tpvb2/4+/trKQ9j\nY2O1KVImkwmBgYEQEZQrV87qLFYHDx7ULtBsddJQh5vd3Nxy7JU2Go2oXr06evbsiffeew/ff/89\nTpw4UeCmgaWmpiIiIgKlS5fWLrAelfR6WZ07dw7BwcFwdHSE0WiEnZ0dihYtivDwcKvnc+/duxc1\natSA0WhE+/btzdpBPD0d6NFDGWovUQJYuPAqPDxK3jPV0By3b9/WepyefPJJzJo1C05OTpg9e3Zu\nn1I2UVFREBE0aNDAJuUVeJmZ2B4djSZNmmRbWNmpUyd4eHigRYsWOH78eH7X0uZu3ryJgwcP2jzT\nU9WqVSEiePfddyEi6N+/v1nHJSQk4Omnn9ZGFtUA3cXFRQsCsmbg2bVrV47f74sXL9Z6dH/99Ver\nn4+Pjw9EBN/nQbbATZs2YePGjTn2PBcYZ89i29Sp2c59LVq0QEhICIoVKwYnJyeEhYU9kutk0tLS\ncOLECdy4ccOm5X744YcQEQwZMkSbcm7uefj69et48sknUaRIEWzduhWdO3eGg4MDXF1d78lQ9aD3\nZPfu3dp7OnnyZKufU/PmzSEiWLFihdVl5SXbBiCxsUBwMN7/r9fOYDCgS5cumDRpEooVK6b9Ljg4\nuMBHZg+TlpaGyMhI1KpVC3Z2dqhSpQqaN2+Ow4cPW1322rVrISJo164dOnbsCHt7+4cOfWdmAuPH\n79d6Q6ZMmaJ9oJs1a4ZevXqhdu3a2fYdyHorXrw4OnToAF9fX+3ifo+lk4XvQ50+YYuTxrfffguD\nwQA7OzusWbMGN27cuGdIUj1J3j1PMi+mgeVGRkYGFi9enC0wLF++PIoWLYqhQ4cWuhGc+7l+/TrG\njh2rvR/29vbo06eP1qsvImjdujUO5iIzXmxsLPr27auNQKifdWdnZ7z22ms4ffr0fY/dsEFJb+3l\nBUREAC4ugI/PVYwZk7tkCUlJSdr3m3oyq1mzpk0C3rS0NK2X89tvvy0UCwtzJS0NiIwEatfG3ho1\nYDAYUKpUKYwZMwYhISFwc3ODq6urNjfakh7CguzKlSsIDQ2Fl5cXnnjiCdjb2yMoKMgm00tv3boF\ng8EAZ2dnbSTdnOxsN24AnTr10TqhZs6cqaXNrV69unaRdvfN2dkZTZs2xauvvoqPP/4YH3/8sZbA\nxFZTg2bNmgURQUBAgE3KU82ZMwcGgwHFixeHiJJx6++//7bpY1jl+HEgOBhwdkaKkxPKli4NX19f\nhIWFwd/fHwaDAZ6entoagVatWiHWog3ECq6UlBRERkaievXqqFmzJry8vBASEoJUS+cS3oe6H837\n77+vfebNceHCRW3Kr7e3Nz788EM4ODhobeXum6enJ9q2bYs33ngDS5Yswf79+5Geno5z585p03j7\n9Oljk/PGnDlzIKIkeyjIbBOA7NunbKRmZweIIK54cbzcty/eeust7eQpIujQoYP2heTt7Y0//vjD\nJg+vp+TkZHz88cfZLh4rVqyoZfJwcXFBeHh4rj9E+/btQ/PmzbWFRZcuXUJKSsoDy0tPB/r0UXp0\nO3bcgalTp2oXfX369Mk2MpOWloYDBw4gMjISISEh8Pf3R+nSpbOdRIoUKYJvvvkmV/XPycyZM62e\nGwkAW7Zs0YaaP/744/veLy0tDXv37sWyZcvw1ltvoVOnTnB3d4eIeWsD8ozJBKxcid/+W7AsoqRZ\n/eSTT9CuXTvtYrpq1aoWDQEXNLdv30ZYWJi2KM9gMCAgIEBbH/Hnn38iIiJC+9J1cHBAcHCwWVP9\nrl69ipCQEO1z4OjoiKCgIKxfvx4BAQHaa+jo6Ih3312B+81o+/prZWM2Nzel3QwaZF0Kd3Wh7tCh\nQ7XvBnP2AbmfMWPG4IsvvgAAdOzYMVumoLJly2oL4sPCwhAVFYVdu3YVzjVRycnARx8BlSopb4QI\nULUqvpw+HQMHDtTeZ4PBgG7dumnZYUQEXbp0yXWCjfymzhLI2llSu3Zt7eKlTp062LVrV67LT0lJ\nwdixY7XOjVu3buH06dO4evXqA4+7fFnZzqBGjSTUr98GM2fO1NZFZV1LeP78eW1KYmBgIOrUqXPP\nfh4ODg4wGAx46623cv087nb+/HnY2dnB2dkZN2/etEmZYWFh2mfM399fS3FqNBrxyiuv4MyZMzZ5\nnFyJjweef17JnyyiXGO9+CL+/vhjtGjRQnut3dzcMGrUKPzwww8oV66cdsFbmPcRunbtGt57771s\ne46VKVNG+5w1b97c6jV3t27d0qY///rrr/jll1/w7bffPvS4kyeBhg0zULt2IGrXro1p06Zp9QoJ\nCcGpU6ewatUqTJgwAf7+/ihfvnyOQYmLi4sW9Hbo0MFmHUwXLlyweTvJKiMjA7t27UJycrJV5Vgf\ngFy6BDg4KI3D0RF45RVg1ChM+m/qkIigcePGWLFiBTIzM7Fv3z6tR9zOzg4hISGFYj5nYmIiwsPD\ns32Q6tati3feeQfHjh3DrVu3EBQUpP3N19fXovnKO3fuRLdu3bQeDPXizc3NDRMnTnzgnPIXX1Re\nfg8PYMoUoEOH9RARDBs2zOyh2DNnzuDHH3/UMih88cUXCA0NRefOnc1+DveT9aSR28W0J06c0C7C\ncpNaWO0BHJtfG5lt2QI89RQgAlORIuji44M5c+agd+/e2hdXuXLltGxYIgJ/f/9CtZFjWloaIiIi\nsnU6+Pn5ZbuQSklJQc2aNeHs7Iy3334br7/+ujYq5+3tjV9++SXHshMSEhAWFqZNVTMajQgICLin\nl+/gwYMIDAyEm5sbSpe+AKNRSbF498Dk33/f2R/hlVeUEURrHD58WJt3q8557927d67KWrp0qTaS\n+c0332gXcuXLl7/vCKb6fVqlShX4+flhxIgR1j0hvdy+DZQqpbwRdesqW9n364cd/6XwNhqN8Pf3\nz5aKOCoqSuuFL1WqFH744Yd8fAKWOX78OIKCgrRAw2AwoGvXrpg3bx4AZUqhmr7c3t4eISEhFl2U\nqOcpNbhXL26KFy+OsLCwB55rz54FatZU3opatYAPPriJIkWUqYkffPDBQx/75s2b2LRpEz766CMM\nHDhQ6/TZtWsXTpw4gS8tyv97f+r6RGvLM5lM2ga7dnZ2WLRoEQBl5DYkJESbQu7o6Ijg4OD8SwTR\ntKlybdW3r7I7fZ06uOXujqKenihRogRCQ0Nx7do17e5XrlzB888/r30vBAQEZPt7QXfx4kWEhoZm\n2+uiSZMmmDFjBjIyMrBt2zZUq1bNqg5fdeSxWLFicHBwgKOjI8qVK4eIiAhkPGQnwgMHgHLllHbi\n55eK99//RGvLM2bMuO9xcXFx2myNrEF76dKl4e7urk0ts1UQoq5nsfVo8YkTJ9C2bVs4OjqifPny\nWLnS/JT3d7PNCEhQkLLr9ogRgKcnIIIrxYujQ/v2+PHHH+/5cKSlpSE0NFQ7mdarVw+7d+8266HO\nnDmD6OhoTJ482WbDVQ8SHx+PsLAw7Ytc7bUeP3689kWYNSvH2rVrtR4IDw8PREREPLCO6kJLNfBw\ncnJCUFAQfv/9d20BoYigZMmSWLToeI4ZTL7/XsnzHhqq7J4uAkybtjVXr838+fO13kU12Nq5c6fF\n5dxNnXbz+eefW3zsrVu3tKHOzp0752qtzbZt27SL/Id9wdjUpUtKDky1d7dsWeCDD7D/3Xe1ixBn\nZ2eMHj0aV69eRXp6OsLDw7WTt6enZ67XSsTHxyMmJgb//mt5BjVLZGZmIioqSgte1d6pnNZ73bx5\nE/3799c+73Xr1sWiRYu0ud3qCVPdL0hdJ5N1BMDPz++h0wOPHbuB/v3v9I3Y2SmjhOr62dmzld/3\n6ZOLnW/v46mnnoKIYPr06bCzs4Ojo6PF+x7t3btXm2o0c+ZM7btE7X1OS0tDbGwsYmJiEBERgZCQ\nEAQEBNyzCWjp0qXxzDPPFOzOnV9/Bb79VjlvhIUB/v7Khkf/dWZNHzv2vlNPT58+ra1REBEEBgaa\ntfA/Li4OmzZtwty5czF48GC8+uqrtn5WOdqzZw8CAwO1c57RaESPHj0wbdo01K5dGyKCLf/tNp6c\nnIyQkJBsPb0Pa8PXr1/HpEmTsp2nGjVqhClTpmifSxFBtWrVsHJlTI6jffv2AcWKKSMgY8cqb0PT\nplfx8ce5S4AwYsQIiAjeeOMNuLq6wmg02mQ0YeHChRARdOrUKddlZGZm4tVXX9UCjJzmyp8+fRpB\nQUHa+1CsWDGEhYVZ3etrtrVrgW++ASZPVnaPrFDhznmkcmVsXbz4gaOekZGR2nmkTJkyZqVljYuL\nwx9//IH58+djyJAhaNWqFdavX2+z6U4Pcvz4cQQHB2f7Hrt7mpnamXV3h2+nTp3MWv8XGxuLYcOG\nZXuMZs2aaXs4qdd3D1qv1L698ha0bw+88QZQpkwGKlVqi88++8zi5xwXF6d1Fnz//ffo3LmzzTaf\nXbBgAUQEzz77rE3KM5lMmDt3rnZ+UjvJS5QokeukMtYFIJs2AStWAFFRgLv7ncbRoYMy0fohtm/f\nruX2dnBwyNbbk5GRoWWeCQ0NvWeqkHqrWrVqrlIDmmPKlClaj6uIMl99/Pjx2RZIFi9e/J7eobi4\nOPTr10+7zzPPPHNP5pDNmzdnmwvv5uaG4ODge6YU7Ny5E+3atUPFim3g6KjMVIiIuPeiae7cO6O0\nuci6q7l+/TocHR1hb2+PYcOG5XrEQaUGQepixJYtW1q0/ic9PR0dO3aEiDItwZrhRPWzpusiw4wM\nZTt7d3dg1CjlROLuDjg6on3lyggICMgx89X58+fRo0cP7fPRqlUrLSf4/ajTIkJCQuDr66tNd2zf\nvr02nScv7NmzRztJ161b16y1Pn/88Yd24WUwGNC3b1+MHz9em5JSvHhxvPLKK9lGHJ9++mmLd0I+\ndUqbOo0iRZR+kjNngKVLgdWrlemLtrJs2TLtwu9///sfROSBPWJ3i4uL03r2+vbtq104tm/f3qyg\nOy0tDceOHcO6deu0cr7++mtrnlLe+vRT5edHHynDtyLKm/TGG0p3/EOYTCZERERou2ZXrlwZv//+\nu/Z3tT2EhoYiICAAderUuSePvp2dHUqXLp1ni3b37NmDTp06aY/n7OyMoKAghIaGZpvG6+3tjXXr\n1mU7dsuWLdl6enNaXKxmmszaW+zr64vo6OhsHVAxMTHayErz5hfQrFnOp+h//lHiQTVo/29QIFey\ndvqo++TMnDkz9wX+Jy4uDk5OTrCzs8tV0o6MjAy88sorEBG4urre87rfbf/+/do+NCLKlGtzesqt\npq6b+eQT4IknlDelXj1lrZSZveQnT57UvkcMBoO2vig9PT3b9dX92of62ZsyZUqePc3du3fj+eef\n184hdnZ2eOmllzBnzhw8+eST2a6Rli5dmu3Y7777Tgu6S5YsiVWrVuX4GHv37kVgYKA2nVAdWd22\nbRsA5btk2bJlqFy5MkQEFSq0hr9/JnJK2nbxIvDaa8DAgXcm/qxYkfsATU3fO2zYMK2jzdpNP4Hs\n13LWbgB+6tSpbB0+6qha586dISLo1q1brsq1LgBR5xcuXKiklfHzAyzM056YmIjXXntN++DXrFkT\nzZo1y3EhsdoL0aFDB7z55pto3LgxRMSsOXuW+uOPP+Dn56d9oYeEhGh7NYgoQ/+hoaEPvCCOiorS\nGkfRokWxYMECREdHaxkK1FGSkJAQXL9+/YH1+fnna6hX706M16gRoM5W+eefOx2H06db/9zVi6e3\n3npLa9i5GRY8ceIEGjRogD///BPHjh2Di4sLihYtCjs7u2yZuKKionDgwIEcR2zUIKhMmTJWZ4hR\nU1EOGDDAqnLM9uOPwPLlytBUaOidOT8GA9CjB5LNGJmIjo5GxYoVISLazvaJiYlIT0/Hrl278NFH\nH+HFF1/UelGy3hwcHFCrVi0t8Msrak/I/PnzLTopp6WlISwsTJvnX7p0acyaNUv7UlM7HOrXr2/2\nxoL3c+YM8MMPygXV3LlAXmyonZKSoqXPDg8P13qczbm4NZlM2ohnw4YNMXToUO1i58qVKxbXRc0n\n37Zt21w8E52sWKGMgHz9NfDBB8CECcBD1ijk5NChQ9oImtFoRM2aNe+bIc/FxQXNmjXDq6++innz\n5mknfFvu36JKSEjAlClTYDAY4ObmhqFDh+K9997LNkWxXr16iIyMvG+AeXdPr5+fH86ePYuTJ08i\nODg423lSDTzuJz09HV9+Ga1NHxFRBp2ybj00bNidiyorZlZo1L2kpk6dChFlOnZu7dmzRxtBefbZ\nZyGiJGqJjIzEnj17zDpHpaSkaO3Mzc3Noqyc69at04I4o9H40A4hq61YoZw/liwB1qxRbrmY1ZCR\nkYH3339fG3GvVKmS1jl1983V1RVNmzbFoEGD8OGHH2rriF566aU8eILKVCj1AtzJyQl9+/bF9OnT\ns11rqdPM7neNdOnSpWwBYtbR0LtnmDg4OCAwMPC+iU+Sk5Mxffp0PP30SS0IHzQIyDq4kpoK9Oql\ntBNXV2Wgyhr//POPdp2ltvVx48ZZV+h/1HbyqdrZkwtZR9JKlSqVLcg7d+6cNh32q6++srhs6wKQ\nqCilkXz3ndkR+f3ExMSgdOnS2XqFypYt+8CLVDUd6zvvvGPVY+dEHT5+8803cfr0aa3xVq5cGR9/\n/LHZw7Dnz5/XPgRZexfKlCmDmTNnWpTxKDNTecm9vZUPf7Nmyhq1H39UdvB8wLpsi3zzzTcQUdKK\nqlOfLF0UHRcXp/Vwd+/eXevJ8PLyyrbLdNabl5cX2rZtixEjRmDhwoXajqDOzs422YDs2LFjMBgM\n8PDw0CcLm7p3w6efAs2bK2+ary9gYTaumzdv4rXXXtN6iCpUqKD1+t4dnHft2hVTp07Fpk2bcPv2\nbcTHx2vTgfLqOasBxLVr13Djxg2LpygcO3ZMC/abNWuGjIwMbdRr0KBBNp1muXgxcOyY0smeF9Q5\n5YMHD0blypXh7Oxs1uaSkyZN0t5DNZOWg4ODNi3HUgkJCdrobW6yjOnGRkNQ6enpCAsLQ7FixbRp\nTupi/aypMO8OkNV9RqwNcHOyZ88eLYiMi4vLNqLZokULrF692uzP9sqVK7Xg1tXVNdtUroCAAIsy\nNiUmAu+/f2fSgqOj0qsLAMuWKYkZHjIoYDY1/e/gwYO1zrjcbLZ65swZlCtXDuXKlcO+fftQu3Zt\n7fVQb/b29qhTpw4CAwMRFhaG6OjobMF7YmKi9r3i5eVl8WgqoEzdUjuyOnToYPHxFrPhEO3evXtR\npUoVbbRMvb4KCQm5b/vYvXs3RJTZB3lh5cqVEFESMFy6dEnbUFFE2bRy3rx5ZiXXMJlM+Oijj7SA\nvHz58qhfv75WlpubG0aPHm32msqrV5WRc0fHOwOzCxYof9uzRwk8PD0tPpXflxpwqQl7vL29c33e\nW7p0Kd5++22YTCZMmzYN9vb2cHR0vGe/tIedpy9evKh1RqujHjklsVi0aJHWyW7pmlXr14DYcBjy\nq6++goigadOmZi34+u677yBiuzluWfXu3TtbVBcaGorFixfnaiTAZDJh4sSJ2gjAxx9/bNXFYHKy\nMod9zhxl8OmuUUmrJScna19So0aNgogyJcRcWadN1a1bFwEBAdoXyqVLl7RMXFFRUdr0upx2nC5W\nrBgMBoNNLw7U0ae8GDW7h9qD9eWXwNatgJVTv3bv3o1atWppPahly5ZFQEAAwsPDsWvXrvv2tKtJ\nH/IiDXFmZiYMBgOMRqPWkxUWFmZxOSaTCYsWLdLWdvTv3x8ikqdTx/KCunu5m5sbfvvtN7MWf/72\n228wGo0wGo1YsGCBNsfW2gw2Q4YM0drw40IdCfnxxx/Nuv/48eMhIhg/PncpmB/kt99+g4gydRBQ\nplR17Ngx13thXb58GW3btkXZsmVhZ2eHwMBAq9K+X7kCDB8ODB0KvPWW8v8VK5RMWLZy+PBh7YJf\n3cxwwoQJFpURHx+Phg0bQkTQpk0bbTpX+fLlMW7cOAQEBOCJJ564JwuXeqtatSqee+45bY1auXLl\nchUEqdQ1KJacEwsKdSRK3TT1YZKTk2Fvbw97e/s8WU/26aefQkQQFBSk/V8dFczNtdahQ4fQuHFj\nFClSBEaj8aGjJw9z6hQQGKhMWnj+eWUk5KuvlOVrtszSrHZA9e/fX5t2nJsAeePGjVqH4LJly1Cp\nUqUcp9WJiBaUDB48GPPnz8e2bdu0YC8qKkpLLV+0aFEsWbLkgY+rJj3o2rWrRfXN053QLXXy5Elt\nKMocsbGxEFGmatiaOm/3559/tqqcLVu2YM2aNdqGYurJyBaWLAH+/FOZvWBrAwYMgIiygNBgMMDV\n1dXs0Zqs06bUkSR3d/eHzmu8cOECfv75Z0yfPh09e/aEnZ0dnJyccOvWLVs8JQDARx99BBHBc889\nZ7MyH8iWiwwAbNiwASLK1A1zqRnAchMYPExiYqLWK6u+1x999JHV5aoXGbZMB60XNTnFuHHjzFoY\nHR8fj4CAALz77rvaRVK/fv2srofaA1+0aNHCmaI3FwYNGgQR0bJKPYz6vezv72/zuixfvhwigl69\nelld1rVr15CWlqatM7J1fv/Fi5UBWxslqspG7QBRN+ysVq2a2b27GRkZeO655yAiqFWrlvYd4+Hh\ncc8UqISEBOzatQuRkZEIDg6Gr69vtpHiqlWrwsPDw6LslDkZPnw4RCxb31VQrFmzBiLKVD5zqdN4\n8yKF/XvvvQcRZaNMQHm/rR3x/vvvvyGiTDWz1ffe3r3K9N1584C86BP7999/tc911uQNljh8+LA2\nHSo4OFibddK8eXOcPXsWMTExmD59Ol566SXUrFkzx4Ddzs4u2xTR559/3qwF5hcuXNAClsUWzG8u\nUAGIyWTSXkBzMhqYTCY89dRlFCuWCVtvIK32pP31119WlXP3egpbb6CUV9Teu2rVqmmL2MxJtzZ7\n9myI3Fm4lnXTQEupaeRs2Qt+5coVbbOgh+XEL4iSkpLg6OgIOzs7sxfkf/vtD2je/H8YPtwGk7rv\ncvXqVYgoi8YHDx4MEclVNpC7qe1m9erVNqilvr755httB2q199fHxwcBAQEICQlBREQENm/enG2U\nNzMzM9saEFudOJs1awYRuWfx5oP89ddfuZrPWxDMnTtXm/JjjqNHr6Jp06to3Nj2GzzOmzcPIoLX\nX3/d6rLUtRTjxo3TLjBs6csvgf37gTwYCML06dMhItnWqpnbu6te7BcvXhzTpk2DiDI10dxEIhkZ\nGTh48CA++ugj2NnZwcHBweq1hK1bt4aI2GRnd72dPn3aok5eAHjjjcOoUycRy5bZvo2o7+/cuXOt\nLmvw4MHw9fXFJ58oaXHbtWtngxresXgx8O+/wJtv2rRYjbqmWZ2GVbZsWbPXVF69elXrvOratasW\ntHt7e993AXpOAbuzszNKlSoFg8GA0aNHW1R/dQaTh4fHAzcCzqpABSDAnYtOc4fQ27RR5ujZas6q\nSt2c5uTJk1aV07JlS4gIxowZAxFlo7LCIDMzUztZLFy40Kxd0deuXQs7OzsYDAZMmjTJrE0DH0Qd\n6rakt8YcXf7bo2b+/Nyll8xvas+GuSfAkyeVNuLlZd1mezk5c+YMRJR57n379tWGfq2lrgkpjCf5\nzMxM7Nq1Cw0aNNCmU93vVrJkSbRs2VLLrFesWDGrv3Oy+vzzzyGiZPB7kOTkZERGRmodL+7u7maN\n3hQ0mzZtgoiylsgcmZl3NqO09VYJEydOhIggNDTU6rLUKbHqGqPJkydbX0GdnDlzBkajEa6urhg5\nciRatmyJTZs2PfQ4dTdnZ2fnbLuqL1An41tI/X6yNqujurbK2sxC+cFkMmlJQ8xNnTp5stI+8mIL\nLVuOdKvfoWqgaouRRz2pgfpLL72kdThkzep3P0lJSdp1ZtOmTbUEJsWLF7c4uUZqaqr2HVO5cmWL\nZ5/07NkTIsr6KHNGsoxSwDRu3FhERPbs2WPm/ZWf//xj23pcu3ZNRESKFy9uVTnXr18XEZH09HSb\nlKcXo9EoL730koiI/PTTT1KyZMmHHrNlyxbJzMyUMWPGyCeffCKpqakycuRIef3113NVh169eomL\ni4ts2LBBzp8/n6sycvL888+LiMjWrVttVqaeWrVqJSLm179KFZEKFUTi4kSOHLFtXZKSkkRExMXF\nJdu/rZWcnCwiIq6urlaXpTej0Sg+Pj6yd+9euX37tty4cUN27dolUVFREhoaKgEBAeLj4yMeHh5y\n9epV2b59u/zzzz/SsmVLWbZsmVSpUsVmdendu7cULVpUtmzZIgcOHLjn7//8848MGzZMSpUqJa+8\n8ors3r1bSpQoIUOGDNG+swqTBg0aiMFgkAMHDkhmZuZD7280itStq/x7/37b1kU9h5jz3fkg6enp\nEh8fL/b29pKYmGiTMvVUsWJFad26tSQlJYmPj49s27ZN2rRp88Bj1q5dK2PHjhWDwSBTpkyRd999\nV9LS0mT8+PEyZMiQXNVj/PjxYjQa5fPPP5crV67kqozY2FiJj4+X8uXLS6lSpXJVRn4yGAxSr149\nERHZb+YHvkED5ee+fbavz+XLl0VEbPJaqu+p+r1V2N6f3r17i8FgkOjoaBk3bpx899130qJFiwce\nYzKZpF+/frJ9+3apWrWqdO3aVRYsWCDOzs6yevVqqVGjhkV1cHR0lLCwMGnevLmcPn1axowZY9Hx\nn376qZQuXVrWr18vCxYseOj9C1wA0qhRIxExPwD57+5i5t3NkpaWJomJieLg4CBubm5WlaUGICkp\nKSJSeAIQEZF+/fpJs2bNZPXq1VK+fHnx8vKS1q1by8iRI+Wzzz6TLVu2aBedIiJTp06VNWvWyLp1\n6+TKlSvSuXNnmTVrVq4f38PDQ7p06SImk0m++eYbWzwlSUtLk6VLl4rRaJQuXbrYpEy9+fr6iojI\ntm3bzD6mZUvlp61jLjVQcHFxyfZva9kymMlvXl5e4uPjIwEBATJp0iSJioqSXbt2ya1bt+TChQvy\nxx9/yFdffSVbt26VZ5991qaP7erqKn369BERkYULF4qI8l20YsUK6dixozRp0kQWLFggCQkJ4uPj\nIxEREXLmzBmZOXOmFCtWzKZ10UPRokWlUqVKkpycLMeOHTPrmIYNlZ9799q2LlevXhURkRIlSlhd\nDgApUaKEzcrUm/oZHDdunEyfPl1+/vlnuXTp0n3v36RJE2nSpImMGzdOPv30U7l586YEBATI5MmT\nc12H2rVry//+9z9JSkqSuXPn5qoM9bqkofqhKYQa/BdR7DMzosjLAEQNGqwNFgBobUM9DxW2AKRi\nxYrSqlUrSUpKEjs7O+nRo4c4Ozs/8Jhvv/1WVq1aJcWKFZOxY8fKlClTxGAwyOeff65dJ1jK3t5e\nIiMjxcXFRRYuXCjr1q0z+9jMzEypUqWKODk5SZMmTR5+gEXjKzrYu3cvRATVq1c36/7//KMMD9au\nbbs6XLhwQZuDZ43MzEzY2dnBaDRqGxNGRkbaqJb6OHHiBNq1a6etzbn7Zm9vj7p166J3796YOnWq\nlmXK2k0DVatWrYKIdfnjs1LzbJcrV87ilHEFxaVLl7RpMubOEf3wQ6WdDBxo27qom421aNFCS6Fo\nzrDxw6gLHw8dOmR9JR9z6nequ7s7Ro8erS0WFBF4enoiKCjIJhtfFRTq/Gdzp3XMn6+0DVtvD6Qm\nI8ht1iuV+v7Vr18fbdq0gYhg48aNNqqlPq5du6YtQs968/Lygq+vb46pkm/cuIEWLVpoU+pssS5q\n586d2jx1czJt3k1NK2yrfRryg5p5qn///mbd32S6s09oLrYkeiA1NXNu9jrK6saNG9m+z8SKqXr5\nad68eXB3d4fRaNRSSoeHhyMmJibHTF4mkwmTJ0/OtiGrrZLNqGtRypUrhxs3bjz0/suWLdOuE4sV\nK2bW+ukCF4Ckp6fD2dkZBoPBrAvYtDQgKEjJ22zmspGH2rdvH0QsyzSUk2vXrmlvhrru4KeffrJN\nJfPB+fPnERMTg/DwcAQGBsLHx+eeDY2qVq2KMmXKmL0I6WFSU1O1LylrN35S54a6uLjgzz//tEn9\n8ou3tzdEzM9M8uefQLVqgK23zFm/fj1ElOxuavC5c+dOq8tVd6S15XqIx1FKSgqioqLg6emp7WUk\nIvDx8UFERMQjmR1L3R/K3IvEEyeAb74BYmNtWw91DyVrgzs1IUj79u21vZWsSSObX+Li4rB48WIE\nBwejbdu22XZvz3pzdXXFk08+iZo1a2oLaa29QM1K3dH5g1ykj1T3jdEljXse2bJli/YdYK6PP1bS\nz1qwbdlDpaenw2g0ws7Ozupd5dUsUk888QS6desGEbnvrugF2e3btxEVFaXt85P1pm6y+tJLL2H6\n9OmIiYnBtWvXcOLECW3T3kGDBtmsLpmZmVqHx8svv3zf+12+fBndu3fX6tm5c2ecO3fOrMcocAEI\ncCcD1R9//HHP306eBL7/XtlYuls3oF8/JT2auh+GLRbZ/v777xC5s4vwlStXEJuLs9ORI0e0RqFe\nnOUmt3NBlpSUhL/++guff/45goODsW7duly9Vg+i7mdgzYaTK1eu1PZa+P77721XuXwSGBho0QL/\n9HRlY6UzZ2y7b4ya1rFr164IDw/HiBEjcj2ylPUkpO5QXRgXehYEBw8eRHBwcLaRS4PBgDJlyuT9\nDs75TN0f6n456TMzgSNHlE1dx49XdvzOug+GragXBRcuXPjvcXPep+dh1I1he/XqpW2+96i0i/Pn\nzyM6OhphYWEIDAxEnTp1tH0L6tSpg2rVqtn88xoTEwMRZVdnS/fjUhOzWLq4tyCJj4+HwWCAs7Pz\nfS/8TSbg+HFg1Spg1ixl41Z1DwxbUWea2GIbhT/++AMiAl9fX21Bdm43cS0IEhMTsW3bNsyfPx+D\nBw9GkyZN7rt7vZoU4dlnn0W6jdP+x8bGws3NDSKC77777p6/r1ixQvtO8vDwsHjvqgIZgKi53N9+\n+20sWbIEo0ePRvv27dGhw2CIINutWDEgMlLp4R09GqhbF8ghbrHI5MmTIaJsegQoK/tdXFwQHh5u\nUY7qf/75B+XKlUObNm1QrVo1iIjVOcgfR+qXS6VKlXJ1Et+9e7c2PDlr1qw8qKH+1GH0Pn365Pj3\n8+eVC6zgYMDHB+jSRQnU5869k07QFj744AOIPDzL0oOkpqYiIiICVapU0XpO3N3dISKFMhNTflFH\nO/z8/LJtPuXj44OPP/5Yy36TF/n8C5Jjx45BRFChQgXEx8dj69at+OSTTxAUFIRevf5EkSLZzyG9\ne9/ZByMsDJg0SRlZt4bJZIKDgwMMBgNSU1Nx7tw51KpVK1dppdV0vq+99pqWZdDWFxoFyY0bN7Bh\nwwbExMTk2fNUL1ItydCophwvUqRIroPJgkLN8nn48GFcv34dGzZswNy5czF+/BI8+eSdzHDqbe7c\nO3tgfPutbaZiqfsU1a9fH4ASkISEhORqk2a10+GFF17QrrWOHj1qfSULkPT0dBw4cEBLm+vn54di\nxYrhqaeewquvvppn50o1tXnJkiW1jo+4uDhtqpuIoGPHjjhz5ozFZRfIAGTevHk5bpJSuXJ1lCwJ\ndOyo9Fh9/TVw6NCdUY/AQKWx2NsD48YBqamWPe7mzZvRoUMHbd6beoGnphYTETzzzDNmDy9lpZ78\nzZlLR9mZTCZtXwVz0jdmdf78ea3XaqCtF0DkI3WaYOXKlZGSkoItW7ZgxowZ6NatGzp3/vqeQL1c\nOeUi69gxICBA2dk1ICD3J5LDhw8jICAABoNBG63o2rWrRTn209PTsXDhQu39ERHMnj0bGzdu1Nr/\no3yhZStHjx5FSEiI1hOlrvcICgrC31m261U3uLLFvhQFWWZmptZrd/cuwG3bzoYIULEi4O+vnCd+\n/fXOPhgVKijtpWlT5dySG5s3b0aXLl3g4eGBIkWKICwsDGPHjtXqM2TIELM3dQWUtrZ48WJttLF4\n8eK5qxhp1LWFFStWRGqWC4WkpCTExsYiJiYGERERCAkJQUBAAHx8fODq6gpPT0/4+vrmY81tQ91n\nKet6MGXaeQvtnFG+PNC5s3KtFR6udFr16wc4OAClSgHWbNF09uxZBAQEwN7eHuXLl8e1a9e0tVt1\n6tTB7t27LSrv559/RuvWrTFx4kSt88oWa1ALuszMTFyzdf7wu5hMJm1j7ueffx7r1q3Tdmt3dXVF\nWFhYrgPyAhmAJCQkYNOmTahWrRp69uyJKVOm4Mcff3zohX9KitJY7OyUBtSwIbB//4OjwszMTKxa\ntUqb9qU2ys6dO2v7WFStWhXvvfceSpQoARFloZO5Q02XL1/GxIkT4ebmhooVKxb6npP8Mm7cODg5\nOcHFxQW+vr4ICgpCeHg4Nm/ejMTExByPSUhIQMOGDbXRrFRLI9ICLDMzU5tec/fQbOvWr6BMGeUC\nKywM2LxZaRuqhQsBZ2eljZQooYwgmjuwd/r0aQwcOFCbo+rq6oouXbpo87nd3NwwZ86ch87pjYmJ\nQf369bOc+Oph9uzZCAgI0H73KEyVy0uxsbFo165dtovspk2bYuHChTm2iUOHDsFgMMDT0/O+beZR\n8euvv2pr1Jo0aYIBAwbgww8/xKZNB5DDWk7Ntm1A9epK23B2VtqPOV/ZJpMJa9asQatWrbT3wtnZ\nWft3p06dMH36dG1fmCpVqljcmfLrr79CRFCzZk2LjqN7mUwmbY1O06ZN0aRJk/smWsl68/f3t3q9\nQkGwceNG9O7dW/vObtGiBV599VV88slC/P477ttGTp8G2re/07EVGGjZupBDhw5hwIAB2jnLxcVF\n6/BdsGABGjRoABEluU1oaKjFr3VSUhJEBE5OTlbvqE53nDx5Ugvs7nTmtMWJEyesKrdABiDW2r4d\nqFEDqFjxBIoW9UJISAjS7hpTz8zMRHR0tLZ5jTrEFBoaqkXOBw8eRNOmTbWeq8DAQK3nQETQpUsX\nbX7v3U6ePInXXntNOwkZDIZCv/A5P12/fh0ff/xxjicFOzs71KpVC7169cLUqVOxZs0anDx5UluM\nVrNmzUdy5Omrr76Cr68vjEYjGjRogGHDhmHp0qU4ceLhCQCOHgU6dLhzInn55YsPHLK+du0aQkJC\ntBOGg4MDgoKCcP78eQDAxYsXtXUpIspu3jmtd4qJidHalIiyPmrevHkICgrSRj3c3NwQEhJi8SZI\nj5vbt2+jaNGicHZ2RmBgoFm9hr6+vhARfP755zrUMH9dvHjxnu99c8THA4MGKe3CzS0ePXq8fN+k\nGup5JOtnukSJEggNDcX169fx008/oWzZslrH1aRJk7RzjtFoRHBwMFKy9g7kYP/+/QgMDIS9vT3K\nlSuH5cuXW/yc6F7ffvtttg4P9cLV29sbfn5+CAoKQlhYGKKiorBr165H7vvo/PnziI2NtfhCPTMT\nmDkTcHJS2sizz+58aOKRv//+G4GBgVrHldFohL+/P1avXq1lTzQYDBg0aBBGjx6tnQtatGhh9rT1\nvXv3olevXlpmNbKthQsXolevXnBzc7Nq1COrRzIAAYDERODdd+dpvYMtWrTAv//+i7S0NERGRmoZ\nNtS1BeHh4TlmhElPT0dYWJgWsdeuXRsTJkzQFv6EhIRku39sbCyCg4O10RO1oTH4sI3z589j3bp1\nCAsLQ58+fVCvXj3Y29vnGJjY29ujePHij/S6m5MnT1p1YoyKAsqXN8HbuyscHBwQEhKS7YIoMTER\nYWFh2giHwWBAQEDAfYOVn376SZtfbDQaERQUlK1+akKB8uXL47333sPrr7+utS1HR0cEBQU9Mgts\n9bB582aL5v4uWbIEIubvFP44i44Gnn12lBY8ZE2hnpqaes95pHTp0ggNDb2nPV6+fFnrDBERdO/e\nHW+88YZ2MdawYUNcvHjxnsffsWMHunXrpp3DHBwcMGDAAK6LsqFr165h+fLl2LFjh00zbT0O9u0D\nWrRIQJUqT8De3j7Hjt7NmzdrmcPUAC8wMDDbIn6TyYTw8HDtmqlWrVqYP3++NjXX3d39gTNOfv/9\nd3Tu3DnbY6gdY2R7uUlffT+PbACi+vXXX1GxYkVtSDzrPOknnngCixYtMmtqzt69e7XpPPb29hg2\nbBheeeUVLWjZt2/fPRF+QEAADh48mNdP8bGXmpqKv//+G19++SVGjx4NPz8/lCxZEhERETlmUqPs\nrl27hQEDBmgXOrVr18b69esRERGBMmXKaO3Fz8/PrF72pKQkhIaGaoFF2bJlERUVBUAJICdNmoQ3\n33xTG01R2wpT7ua9lJQU7TvQ0nnWj6MrV67ghRde0NrAc889h6lTp2Zbt+Tt7Y3w8HAkJyc/sKzI\nyEhtGkPp0qUxa9Ys1KhRA61atco21SSni7agoCCbpTYnspW0tDSEhoZqIxbNmzfH4cOHER0djWbN\nmmmfYXd3dwQHBz8wMNi/f3+2a6y33npLW3+rJpRQqSOPajIBdeQ8ODg4V4uhKX888gEIANy6dQtB\nQUEoW7YsSpQogbp16yIyMtLiBa7JyckICQnRgoz69evj66+/1hbjqr24d0f4RIXBxo0btR7drKNK\nvr6+2Lx5s8Xl7d27V9tITERJE/j2229rCRnU0RS2FX2NGqX06g8ZMiS/q1JofPnll/D09ITBYNC+\n6xs0aGDxeeTUqVNo166d1iZ69+6No0ePwmQyITo6Olt7US/a7jfNl6ig2LBhAypVqgR1XaD6GS5X\nrhxmzJhh9qideo2lBjRPPvkkZsyYoW3mqc5gqVOnTo5THqlwMQCAPCauXbsmcXFxUr16dTEYDLku\nZ/PmzdK/f385ceKEuLu7S0JCgri6ukpQUJCMGTNGypcvb8NaE+knNTVVwsLCJCMjQzZs2CDvvPOO\n+Pv757o8ALJw4UIZO3asxMfHi4ODg6Snp4ufn5+EhYWJj4+PDWtP5jhy5IjUrl1bihQpIufPnxcP\nD4/8rlKhcPr0aVmzZo1cunRJWrVqJc8++2yuziMmk0lmz54tEyZMkNTUVKlQoYI4OzvL8ePHRUSk\ndOnSMnLkSHn99df53lChcfPmTRk+fLhUqFBBVqxYIcHBwRIUFCQuLi4Wl7V+/Xrp37+/nDt3Tlxd\nXSU0NFScnJxk9uzZcvbsWRERqVy5sowaNUpeffVVcXV1tfXTIR08VgGILSUmJsqkSZOkevXqcvbs\nWXnjjTekZMmS+V0togLp/Pnz8tVXX4mLi4vUq1dP2rdvn99Veqy1a9dONm3aJBERERIUFJTf1Xks\nHTx4UF5++WW5cuWKXLhwQSpWrMgLKnokmEwmMRqNVpURFxcnw4YNk2+//VYcHR0lLS1NRETq168v\nISEh8uKLL4q9vb0tqkv5hAEIEdFjZvny5dKnTx9p2LCh7NmzJ7+r89hKS0uTs2fPyqlTp6RNmzbi\n4OCQ31UiKlC++uoruXnzpqxZs0ZGjBghXbp0sWoGCxUcDECIiB4zqamp8s4778grr7wiDRs2zO/q\nEBHRY4YBCBERERER6ca6SXpEREREREQWYABCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERE\nRES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6\nYQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBC\nRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERE\nRES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6\nYQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBC\nRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERE\nRES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6\nYQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBC\nRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERE\nRES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6\nYQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBC\nRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERE\nRES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6\nYQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBC\nRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERE\nRES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6\nYQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBC\nRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERE\nRES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6\nYQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBC\nRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERE\nRES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6\nYQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBC\nRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6YQBCRERERES6sc/vChAREZENmUwi+/eL\nnDsnUqWKyIkTIvXrK/8mIioAGIAQEREVYunpIrt2iZz766IExASJbNkicvOmSKlSItOni/TvL7J0\nqc0DkOvXr8vIkSOlWLFiUqJECSlevLiUKFFCSpUqpf27ePHi4uTkZNPHJaLCzwAA+V0JIiIiMs/y\n5SK9eytxxsSJIjt3iiQlibg5pkmCXVGR5GSRqlVFnnpKpEMHES8vZUTEx0fkmWdsUoeFCxfK+PHj\n5erVqw+9r5ubm7Rp00Y+//xzKVu2rE0en4gKNwYgRPlt7VqRhASRlBRl6sSAASJffKH8JNu7eFG5\nXbp053V/5ZX8rhWR2ebPF6lcWYkz3nxTmWlVu7YSb0z3ixHPFrVFKlS4c8CaNSLPPSdSsaLIoUMi\nbm5WPf7JkyelYcOGkpCQIMOHDxcnJyc5e/aspKenS1JSkiQkJEhCQoJcv35drl27JmlpaVK8eHFJ\nSUmRkydPSsmSJa18Bej/7d17fI71/wfw97373snmfNqcz4dhZEI5xpLDHJLlOGcjMUkaooViElZ8\nyyJMKlbxc4qsiJyNqM0YI4UO5jQzO79+f1xdF2OH+7R7G6/n47HH1O7ruj734XNf7/fnSFTUcQhW\nYZeYKHLmjEhMjHLj6NlT5No1kXr1lDsOFX2//y4yfrxIaKgSFH/1lciffxZ0qYq0a9dEzp0TiY1V\nfk9MeE/cDnyr/EdiokiZMiJz595/3YmKkAoVRHQ65d8bNyojq8qVU//6/KMH9Ogh0rq1yOHDInPm\niLz/vtnXzszMlBEjRsidO3dkwIABsnTpUhk/fryEh4c/8thy5cpJrVq15ODBgzJgwADZtWuXfPTR\nRzJ37lyzr09kNYyvClT+JyCnTikBVr16IjduiLRpI3LggPKbRERk1SpliO6BA1nrwrlYSNw/FZW+\ndZW9vXIDCQtjBXlcVKggsn69iLOziJubMraCQbHJXntN5I03RFauFJk9O+vfBnT4W9x++UX5j3Ll\nROrWFSlV6v7rnp4uYmB7DBUNzs5KDvHSSyK+vkYcYGcnsmyZSKtWIkuWiAweLNK0qVnXXrRokezd\nu1fc3d1l2bJlIiLStGlT6dOnj8THx0t8fLzW8xEfHy8pKSlSunRpmTNnjuzatUuWLl0qU6ZMkVKl\nSpl1fXoA46s8Mb4qvPJ/CNaaNcq7HxamDHto0EAkIUFk6NB8vWxRsny5EoMmJ4uMHCmSknL/b4mN\nW4uLLkmpDA0bipQvr3StX70qUqmSSO/eBVdwsp4//xSZNk25kQQFFXRpiqRVq0Tu3lVaht9+W6RO\nHSXPqFdPpF+Ts9KoasL9xEN1/rzIm28qyUc2LbhEhdHixcrQq3HjlHtG1apK20WeAgLk38OHZV7l\nyrJk40bRqd0oRoqJiZHmzZtLcnKybN++Xbp3757jYwFIfHy8JCQkSO3atUVEpFOnTrJnzx6ZP3++\nTJs2zaRrUzYYX+WJ8VXhlf9Nfi4uIps3K2+mXi/Sp4/IunX5ftmipHRp5WVKShKZPFmkZEmlLnh4\niDjWPPzou3TsmMiLLyotuZ07WzyeV7V37145ffq0rFy5UurVqycNGzaU+vXraz/Ozs5WuQ5l459/\nRL78UqRFCyYgZkpMFOnWTeSTT5TGwKzqZ3+Qq6vIzp3K3em335SlSq0oNTVVDh06JPPnz5eKFStK\n48aNpXHjxtKoUSOpVq2aVa/1WImOvr907LVrSr04flz5TXL9uvL7xg2Rli1FBgxQRm7mJWnuXGnW\nsKH8deyYeK5eLSNHjjT6munp6TJs2DBJTk6WsWPH5pp8iIjodDopX758lvkeM2bMkD179siSJUsk\nICBAihUrZvT1KRuMr/LE+KoQgy3MnQu0agVs3WqTyxU169YBQUFAVJQJB7VrB4gA8+ZZpQynT59G\niRIlULx4cYjIIz86nQ41atTACy+8gEmTJuGTTz7Bb7/9ZpVrE4Bdu5T38/nnC7okRdLdu4CdHeDu\nDmRkmHjwhAnKaz9ggNXLNXLkSOj1+mzrVIkSJeDl5QU/Pz8EBwdjy5YtuHr1qtXLUCStWaP8XrsW\nCA4GNm4E1q8v2DIVIv7+yke2f3/l94QJxh/75ZdfQkRQpkwZ/Pvvv0YfN3PmTIgIatasiYSEBDNK\nrXjmmWcgIvjoo4/MPgc9gPFVrhhfFV62SUAGD1bezLVrbXK5osbHR3l5vv3WhIN++EE5qFQp4MYN\ni66fkJAADw8PiAh8fX3x888/Y+jQoRg4cCC6d++Ohg0bwsHB4ZFKM2zYMOzcudOia9N/NmxQ3s+X\nXy7okhRJR48qL5+nJ3D5MpCWZsLBV68CTk5KBmPFL/0PPvgAIgJnZ2d89tlnWLZsGQICAtC5c2dU\nqFAh2xuRiKB8+fKIjIxEenq61cpS5GzYAGzZAuzeDXzxBZCezgTkAaNGJaBixQy8/HImRJQAyxgb\nN27E1atX4e3tDRHBU089hXfffRehoaHYuHEj9u7di+joaPz9999ZPn+RkZGwt7eHnZ0dfvrpJ4vK\nvmnTJogIqlatipSUFIvORWB8lQfGV4WXbRKQ3r2VN3PTJptcrqipU0d5eYYNAzp0APbuNfLAjh1x\ntWJFfLdkidnXzsjIgI+PD0QEzZo1w927dxEVFfVIS22LFi3Qu3dvDBkyBBMmTNBuYJ6enmZfmx6w\nfLnyIfD3L+iSFEmrVikv36BBQPXqgKMjcOmSCSd45RVkurjg+NSpVinPzp07odfrodPp8NVXXwEA\npk6dCgcHB3h4eMDPzw9BQUEIDg7Ge++9h4CAAHh7e2uJSbly5bD1SW/RTE8HZswAfH2VrJI07du3\nh4jgxRdfRPHilbF8+Vd5HhMREQF7e3tUq1YNs2bNgpOTE3Q6XY6JsIigbNmyqF+/PsqUKQMRwVQr\n1I/MzEw0adIEIoJVq1ZZfL4nHuOrXDG+Krxsk4B06qR8An74wSaXK0pSUgCDAdDrlV5UEcDYBqZL\nhw6hWLFiKFGiBK5fv27W9WfMmKHdaC5cuAAAiI2NxahRo9C2bVuUK1fukZvS22+/jZSUFC1YOnbs\nmFnXpvuiPvkEv7Vvj9/nzy/oohRJU6YodefttwGdDnB2VuJXY9374w80dnODXq/H2bNnLSrL6dOn\nUbJkSYgI5syZo/3/wYMH5xjolS5dGh07dkR6ejqmT5+utZY98Tw9lTf21KmCLkmhoraodu3aFSKi\nJbk5OXPmOkqVKgURwYABA2Bvb6/9e8aMGRg9ejT69OmDdu3aoWHDhqhQoQLs7Oy0z2e5cuXg6OiI\nK1eu4MCBA9i8ebNF5V+3bh1EBLVq1Xqye/qsgfFVjhhfFW42SUCSunVDRtmyyDhyxBaXK1KiopRK\nUbcuULas8u+//jL++BdeeAEigunTp5t87Y0bN0Kn08FgMGD37t05Pi4+Ph779+/HypUrMXXqVOz9\nrwnh9ddfh4hg3LhxJl+bsgoMDISIYD4TELO88IJSdxYsUH57eZl+jjFjxmhd3+aKj49HnTp1ICLo\n168fMjMzs/z91q1biIyMRFhYGAIDA+Hj4wM3NzeICOrUqQMAuHLlCvR6PRwcHBAfH292WR4LLVoo\nb+jRowVdkkJFDU7atWsHEUFERESOj712TWkFbtv2HLp374PSpUsbdc/IzMzEtWvXEBMTgw4dOmh1\nQ0RQuXJlJCcnm1zuuLg4xMbGIi0tDW5ubnB1dYW9vT3c3d3RpEkTdOzYEb6+vhg/fjxmzZqFDz/8\nEOvWrcPOnTtxw8KhMI8rxlc5Y3xVuNkkAWnYsCFEBNHR0ba4XJGyceP/oWrV9vDzW4QOHTbi+edP\nG33s1atXERkZCZ1OBxcXF/z9999GH6tOirJkMmBMTAxEBCVLlsTdu3fNOsfDbt++jalTp2Ls2LEI\nCAjAX6Z8WxRhY8eOhYhg+fLlBV2UIqlhw2fh6TkeS5bsQ/v2v+H11/8x+RyXLl2Cg4MD9Ho9YmNj\nTT4+NTUVzz33HEQEzZs3N6lO/PXXXzj1QCt/ly5d+HkAgDZtlKjh558LuiSFRkZGBgwGA3Q6HQ4e\nPIhNmzbh2rVr2T723j3g2WeVl7BpU+D55/9F6dK10a9fP2SYsFrD4cOHte/6pk2bQkTw6aefmlTu\n9PR0PPPMMyhWrBhCQ0Ph4OAAZ2fnXIeAPfjz1ltv4fjx4yZd80nA+CpnjK+yKmzxlU0SkCpVqkBE\ncMmkQdlPhvfee0/rChcRtGrVyqjjVq1aBWdnZ2zduhWtW7eGvb09XF1ds6yqEx4ejqioqEe6uG/c\nuKG10vr5+VlU/tatW0NEEBYWZtF5AKVyqCukqN3/rq6umDFjxmPdEpyeno6+fftCRLBhw4aCLk6R\nc+PGDYgIXFxcEBAQABHB+++/b/TxV69excCBA3Hjxg2MHDkSIoIOHTpg9+7dJn3u1CTS3d0df/75\npzlPRaMOUXn22WctOk+R5+2NzFKlkPbjjwVdEotER0fD398fs2fPxooVK8wednTq1Ck888wzKFGi\nBPR6Pdzc3BAYGIjff/8928cnJCi9g9WqAU8/rSQiPXsmICkpyeRrq8n1wIEDteFTaSas9qDe6ypV\nqoRGjRpBRPDqq68iOTkZV65cwalTp7B7926sX78ey5Ytw+zZszFx4kQMGjQI9erVg4jAx8fH5HI/\n7hhf5Yzx1X2FMb7K9wTk33//1caesgv1UUOHDoWIaL+NGf6xevVq7QPk7+8PR0fHbFdRUH8cHR3h\n6emJfv36YcaMGXj66ae1VlpzbkQPWrFiBUQE7du3t+g8iYmJWjd/tWrVsGPHDvj6+mqTJF1dXREY\nGGj2WMzCKCMjA+Hh4WjQoAGqV68OEcGuXbsKulhFzr59+yAiaNmyJTp37gwRwXfffWfUsX/++acW\n3AwfPhyzZ8+Gk5MTDAZDlvkZbdq0gb+/P0JCQvDzzz8/sgzp4sWLIaKseHXECkMhkpKStHkkZ86c\nsfh8RVW3bt1Mej8Lm7t372LatGnanAt1SWYPDw+T5lHcuXMHr7/+uva5LF++PGrXrq19Rg0GAyZO\n/AgREcBDo/6Qmgr06aMkH9WrAyY05GYREREBEUGFChVQt25diAi++OILo449efIkHBwcoNPp4Ovr\nqyUwd+7cMer4+Ph4uLq6QkTYC/IAxle5Y3ylKKzxVb4lIAkJCQgODkaJEiW0NzM1NTW/LmcT6nMa\nPHgwevbsmeu4W2O1bNkSIoJ33nkHrVq1wocffpjr47/6CmjV6m84OpbA8OHD4eTkpFWU69ev4+ef\nf0ZoaKg2vrxWrVpZJhOKCKpUqYISJUpY3EoLKDdGdW3rmJgYs86RmJiIjh07apUjLi5O+9upU6ey\nrShF+cs2NTUVq1at0lpJ1ElqZcqUwcCBA3Hv3r2CLqLVpaSkIDQ0FC+//DImT55s8VKeD/r4448h\nIhg5ciQqVqxodGvg779noHHj1trN4p133oGIsiZ7x44dtZbm7G466rrtPj4+8PX11Va8smYP1ogR\nIyAimDVrltXOWdT06dMHIoJNRXCFn61bt6JGjRra58XPzw8rV67Mkji0bt1aG/Odky1btqBatWpa\ny6W/vz9u374NQFke19/fH05OTmjU6AhElPkewcGA2qi5d6+yMEOJEpavMq22oKq9IB4eHnkO5UpO\nToanpydERKsrdnZ22Ldvn0nXnjx5staa/aRjfGUcxleFO76yegJy7/ZtzJ49W2u9ExHUq1cPzs7O\n8PX1LZJDaRISEjBnzhytpUF9s3Q6HUaMGGHWm6VWtlKlSsHBwQFjxozJczze2rXKVgUiQEDARZQo\noaygMGHChEcmuz7ozp07iIyMxBdffIHx48drYxojIiLw0ksv4csvvzS5/A9Sh61MmzbN5GPv3r2r\nVY6qVatmqRwPOnz4sLacnYigePHiCAwMxM2bNy0qu02lpgIrVuB8r17a86hbty4WLVoEPz8/7YvM\nw8PDKq3ohUFycjL+97//aQGU+mNnZ4eJEyciMTHRovOnpaVlmSh47tw5bNq0Kdf6AAAXLwI1awKN\nG99C+/bdERQUpNXp//3vf1kee+XKFURERCAkJAT+/v5o06YNihUrluX5lCpVChNM2Q3OCHv27IGI\noHr16iaN13+cqEMT1xehPUCuXLkCPz8/7bPRrFkzHDp0CADw1VdfoVKlShgyZIi2+ICIwNvbG7/+\n+muW88TFxaF79+7aY5o3b57j98Jff/2LuXOBqlWV+4P8N/FW/WpfsACwxpYC6h4eVapU0ep0Xj05\nb775JkSUDQzV5Outt94y+dpXr16Fk5OTVVaqK6oYXxl/TsZXhT++sl4CcvcuEBKC9EqVUL9mTYgo\n4+0GDRoER0dHrQvazc0NX3/9tdUum5/u3r2LkJAQrVVVHeahBiMuLi4QEVSsWNHoMXoJCQmYO3eu\nVtkerHBq9pndhMKMDKBtW6VyjBypLDPasuXfCAiYkmew9bC2bdtq3ZEiynh3Sxw4cEB7b01phbl7\n9642rrhq1ao4f/58nsccOnQoS0UpEkOWUlKUfT6qV1feQDs7TPT2xgcffIBhw4ZpdcPV1RX169fX\nhlQEBgaatdJMYaD2eKjjk0UEjRo1wqpVqzBv3jytS7tmzZq5rhCSE3X4mvp6VahQAXq9Hn5+ftpy\nhzk5exaoUkV5K9q0AebM+RdOTqWg1+ux1sjNvNLT03H27Fl888038PLy0hKg7du3o3fv3lYZNpWZ\nmYma/32X7tmzx+LzFSXnz59H165dte/eojAZPy0tDSEhIVqLpYuLC4KDg7OMEe/Xr59WHzw8PDB8\n+HCtl81gMMDf3x+XLl1CSEiINuSoVKlSCAkJMWreSEYGEBGhbL42dy6wbJmyQXZ4uHWe44N7eAwa\nNAgdOnTAwYMHc3x8dHQ09Ho99Ho9XnzxRS0hM3cDQnWlulGjRpn7FIomxldGnZPx1X1FIb6yPAG5\ndw9YsgSoWFFretk+cCCGDBmifYHqdDoMHjwYnTp10p6Yj4+PVbqo8oMaPLm7u2vlbdWqFcaMGYOS\nJUtqleH8+fPamPO8nlN2la1Nmzb46aefEBkZmeVN79JlC956C3h4OF5CgrIvl4uL8lKPGfPoeF9j\nrF69WntOlnTv7d+/H4cPHwYAVKtWDQaDQdtk7b333sM333yD3377Ldsg+sHKUaVKFaMqx4O2bdsG\nEWWFiELt+vWszZKNGgGLFiFuyhRtPLfBYMDw4cMRGxuLe/fuITAwMMtY8aLUG5KcnIzQ0FBUrlxZ\n+zw3btwYYWFhWgC1YsUKLF++XAvcdTod/P39jRoPnpqaihUrVqBWrVra+evUqYMuXbpor6eTkxOm\nTJmCa9eybw3csQNwcFA2pXrtNeVtadXqnzz3UsiJOgelcuXKWmuVOS282Zk5cyZEBCNGjLDK+Qq7\ne/fuISgoSBv64Orqirp166JEiRIIDg4utMMTjx49isaNG2ufyf79++PKlSuPPC4zMxPh4eFZhmF5\neXmhV69e2udXXRlKHbb1zz+mr+imXEtJPLZts14CAtxfIKF27dpGJUWff/45Ro4cCZ1OB0dHxyyr\nvZkqLi4OBoMB9vb2OU68f6wwvmJ8hcc3vrI8ATlzRtnlRURZs33YMGSUKYMm/01S8/b2RmRkJADl\nyzcsLEzbVbVkyZIICQkxanhBUlIS9u/fj8WLF2P06NEWFzs7qampCAsL01od1ZvD4MGDs2TU48eP\n147J6zmpla1SpUra8c888wx+yGbToKNHj6J37/4oUyYDIsqY3bffBtResCNHgGLFlJd63DjzKgeg\njAlUW93UVilTu/f++OMPVKxYEY6Ojli1ahUcHBxyHC+v1+tRu3ZtdO/eHVOmTMGyZcu01R2qVKmC\nc+fOmfwc9u/fDxHB008/bfKxNteli5J4LFyofKv9tzOST9268PPzy3Y4weHDh7XlFU3tDUlNTcWJ\nEyewYsUKjBs3Dt9//721n9EjEhMTERISkuVz7unpifDw8CwtSOfOndN2YB4xYgRmz55tVG9ISkoK\nwsLCssybqVWrFkJDQ7WVeC5evAh/f3/Y2dnBzs4OzZqdRGAgcOvWo+f76SdgwgSlLtnbA998Y/5z\nz8zM1AJKdTJ6lSpVrLLB2rlz57Qu/cIaUFjNrl2IfOklLajq3bs3WrRokWUCaI0aNbBu3Tqzh6TF\nxMRYPYlJSUnBkiVLoNPpULt2bezYscOoY0JCQrSNyHQ6HaZPn46OHTuifPnyKFu2bKHt3U1LS0Pt\n2rXh4OAANzc3NGvWDF26dMGgQYMwceJEzJ49G8uWLcP69euxe/du7N+/X+sJXbBggcXXHzRoEEQE\nkyZNsvzJFHaMrxhfPcbxlWUJyMmTwP/9H+DvDwwffn+nFxHsHTVKG/f6sKtXr2pje0WUzZQeHrJw\n5coVhIeHIyAgAG3atIGjo2OWF71mzZpmvbC5OXDggDYO39PTE2PGjMk2o87pOb30381TRNC2bVsE\nBwdnqWxqQJaXw4eBrl3vN5p7eQGffQYkJSlDRkaNUroMLeHv7w+R+5MJTeneu3fvnrbSQ4cOHbSb\ni7+/f5YNdXr27Im6detmWVFI/SlXrhzKli1r9nu4cuVKiFi+zF2+2r4dWL8eCAkB+ve/P8DUwQHw\n90daHi14D/eGNGrUCEcf2pAtLS0NUVFRCAsL0+rKw2vrV65c2azxo8Y6dOgQypcvr13v6aefxtat\nW7Ptus7IyEBoaKg2h8Ld3R1Lly7NsTdETTwe7PHw8PDI0qPysBMnTmDcuEVa/SlXDli0SGlMVL3y\nivI3R0dgyxbLXwN1DsnQoUO1sv5o5tKxt2/fRrt27fDDDz8gIyMD7u7u2g22cuXK6NKlC6ZMmYKN\na9cqG/RZOI+mwF29qtQPEUCvx9QXXkC/fv20z32VKlXw8ccfo1mzZtpnoFGjRkZ9l6anp+Pnn39G\nYGCgltCPHz/e5AnQubl+/TrU3hpTk5s7d+4gODgYpUuXRkxMDHbu3AkRwQsvvGC18uWHY8eOYfz4\n8dkGRNn9FC9eHG3atLFKUh4dHQ07OzsUK1bM7N6hIoHxVZbnxPjq8YuvLEtAVq9Wfi9YoLTqiii7\nHhl54/3qq6+0HV2dnZ3x0ksvoVu3blq2++CPwWBA06ZNMXbsWDz11FMQEaxbt86i4j9M7fZft25d\nli67du3aGb1qz/r167VKpQaCjRs3fqQl2BiHDiljedesAT75BPj2W+Dzz83PzB905MgR7YPaoEED\niOQ9mVClDjOpXr26lmm3bds2x3G9KSkpOH36NDZu3Ij58+drXcVVq1Y1e+WOqVOnQkTw7rvvmnW8\nTaiTmZcvBzp3VhIPPz/AxO7QgwcPZpkb0r9/f/j7+8PLyyvb5QHt7OxQv359DBo0SNubok2bNvnw\nBJUhV0OGDIGzszNat26NLVu2GPU5P3/+PNq3b6+VuV+/fggKCtKeT40aNTBx4sQsQ7maNGmC8PBw\no1u/9+27P65XRJl0rq6KuHCh9ZIPQOl9UXsqpk2bBhHzdlR/cFK9l5eX9jl3dXXVhlyoPy+2aqXN\nKULt2kDv3soYgqK0Gs7u3UpTpIgy9mH4cNzs0AEOej0MBgMCAgK0JY/VeT8PBh3e3t6PLMv6zz//\nYPXq1ejXr98jrYbquHJLdrt/WGJionbuIUOGwNvbO9vhV7lRE2514YGOHTtarXz55e7du7h06RKO\nHz+OnTt3Yt26dfjwww8xa9YsjB8/Hr6+vujYsSOaNGmCrVu3WnXIVK//FvLw8fHB5s2bceDAAZw5\nc6ZITsLOEeOrRzC+elRRjq8sS0DCw5UMfdcu4P33ATOGety8eRP+/v7Q6XRaZVFbRn18fBAUFIQt\nW7ZkmY0/b948iAgmTpxoUfEflJaWBhGlOwsAdu7ciVatWmGLGRFKfHy81rW+ZMkSkytGdtavV1Yx\n+fxzi0+lUXe0VdfI7tWrV57HqENMXF1dtQmV1apVM6klKj09XWuNNHU3XZV6AyrUE+6+/lpZ1y8s\nDIiJAS5fNvtUD/aGqHuG5FRXHlzL+/r169DpdChWrJhVWh8fdvXqVYgok8BNlVNvSPPmzSEi2kov\nTZo0ybXHIy8REUor19ixwKRJwJ9/AuvWATksCGK2du3aQUSwcOFCLRkxdp8DlbqRYrly5fDBBx9A\nRGBvb48ffvgBmZmZuHDhArZs2YJ58+bh66lTle2tHRzuZ1mlS9/veVuzxrpPMD/cvg24uwPPPQe0\na6c9j53Dh+c4V+DevXt4//33Ubp0aS3h7tu3LwIDA9G6detHlsVs2rQpZsyYgYMHDyI6Ohoiyj4a\n1qoP6r3DYDBoe8qYuwjBwYMHIaIMI6GcHT58WGslfvhHr9ejQoUK8PDwQLt27fDiiy9i7ty5BV1k\n0zG+yhbjq9wVpfjK8jkgVvoSV7uqBg0alOdY5927d0NEWTHBWtTdlEuVKmX2Oa5fv46xY8di3rx5\n2pCJyxYEnfltyZIlWiti48aNsXjx4lwf/+OPP8JgMECn02HUqFEQUSb8Hjt2zORrr1+/Xqtc5qz0\npN7of7N0Yfv8ZsJOwcZQWyZatmyJffv2GRXgqnsRREdHW7UsALSArmHDhliyZAkCAgJMnnB3/vx5\nLXifPHkyUlJStCU+Fy9ebJUbTGamMlLps8+Ajz6637hoTeqmUZ07d9ZWQjF2VS1A2X1XTTiWLVum\nDYv4+OOPcz8wLQ04fVoJWFasyNrzVphFRSldUF9+CTg5KclH2bLAypVGNUPeuHEDgYGBcHZ2zjKE\nxMnJCd7e3ggJCcl2Pxh1E73cVm8ylTpcTJ2Ibu5E6+PHj0NEWXKXcnft2jUMHz4cPXv2xDPPPIN6\n9eppSenDP3Xq1EHr1q1N2rm9UGB8pWF8ZbyiEl/l+07oxvriiy+MzhITEhJQvLg7nn76DSQnW6G/\nDMClS5cgoow3NldMTAxEBPXr19e6/m9lNwO2kIiPj4ejoyPs7OyM2rhNXdt++PDh2lCZ1WZGcpmZ\nmVoLwdKlS006NjU1Ffb29rCzsyu0q+Lkl0OHDkFEGe9qLHXsrCnBsLHUJQKfffZZLeg2Z3x9eno6\nli1bhrt37wKANo8ip3XLzbVqFRAbC0yebNXTAlDmbhQrVgx2dnYIDg7Wbj7G2Ldvn1an5s+frw0z\nMKsV8uuvlSB+2DAgH5JOq1F7aNauBby9leGJ//5r8mkuXbqkte4uXrxY+wzlZNKkSdDpdHj33dw3\nJTOFOhxE7b0zJ2gAgKioKIgoc1zIPGlpafj7778RFRWFvXv34ptvvtEaYfLjO7AoYHxle4yv8lZo\nEpC4uDiThnJ4eioNZjnMwzLZgy255jp69ChEBC1atNBaxAp7i8vLL78MR0dHuLq6ok2bNvD390dw\ncDC2bNmCuLi4LOPtMzIysqx0NGXKFIuurW5q5e7unmfQ8KDTp09DRFkG8kmTlJQEg8EAg8GAJHVC\nQx6WLduJdu1WYubMP6xeHnW5vh49emitv9ZoNVFbMq8/vFZiITdgwACIKLuXe3t7G3UDuXDhgjaJ\nf9KkSdqN4/nnnzf/+yMgQPmCLMwrBW3YoPSA7N6t7JVjgVdffRUigqCgoDwfu2fPn3Bzu4ymTS26\nZBbqKj6tWrWCiGD//v1mnSc2NlZrsSfrWbNmjXZ/fxI39WR8VTAYX+Wu0CQgALRWrIsXL+b52FGj\nlAoSEmKdax8+fNjibscff/wRIqINJylWrJh1CpeP/vjjD8yaNSvbbmv1OTz11FMYMGAAZs6cqY0t\n7Ny5s1Uqf8uWLSEi+OCDD4w+Rq1Y3bt3t/j6RZG6EZix+4N8/71SV9q2tX5ZPv/8c4gIhgwZok0Y\nt7RbPDMzE3q9HjqdLl/mreSn7777DiLK7sTGDh3r2rUrRJQJtb6+vlornzk7AGsiI5U3vXz5wj0p\n3Urvr/q6t2jRIs/HpqYCJUsqL4+15kWr9y61F9DcFdDUluKqVatap2AEQGnVVXtBvv3224IuToFg\nfGV7jK9yZyeFSMuWLUVE5MiRI3k+tlUr+e+x1rl2YmKiiIi4urpafA5nZ2eLz2UrVatWlTlz5siV\nK1ckIiJCQkNDJSAgQLy9vaVWrVpy7949+eWXX2T9+vXy7rvvCgCpWbOmbNiwQQwGg8XXf+edd0RE\nJDg4WO7cuWPUMSdPnhQREXd3d7l586akpKRYXI6ipHnz5iIicvz4caMe7+Wl/D55UiQjw7pluXnz\npoiIlC5dOsu/LZGQkCAZGRlSokQJ0ev1FpfRlrp06SLu7u5y/vx56dKli8ycOVM2bNggUVFRkpaW\nlu0xa9askREjRki9evXk66+/ltKlS8uWLVssex29vESaNBG5dk1k507zz5PfrPT+durUSVxdXeX4\n8eNy5cqVXB9rby/y/PPKv7/7ziqXF0dHRxER7TvR3O8k9TxP2ndafrO3t5c33nhDRES7jz1pGF/Z\nHuOr3Fn+DK2oVatWsm3bNjly5Ij0798/18e2bq38tlYFUd+c4sWLW3wO9SZiyblsrVKlSlKpUqVH\n/v/t27fl7NmzcubMGTlz5ox06tRJPD09pWzZsla5brdu3aR9+/ayb98+Wbp0qcyYMUP7282bN+XC\nhQsSHR0tp0+f1n5fvHhRnJyc5LPPPpPPPvtMRJQb/9A9e+RayZLipNdLMTs7KabXS8vixaWak5O0\nLlFCitkVqnzbbF5eXhIWFmZ0AlK2rEj16iKXLonExoo0bGi9sqhJR/HixSUpKUkcHBykWLFiVjmn\npYlMQdDr9fLFF1/IwoULZceOHfLDDz9ofzMYDFKtWjXx8PAQLy8vadSokXh4eEiDBg2kW7du0r9/\nfzEYDPL1119LvXr1LC/MkCEigYEia9eK9Oxp+fkKMUdHR+ncubNs3rxZduzYIaNHj8718T16iHzz\njcj27SKvvGKd64uIuLm5Sd26dbUgyRh37tyR2bNny+jRo8XNzU1EmIDkh1GjRsm8efPkl19+kYiI\nCOnSpUtBF8mmGF8VHMZX2St0CYiIcRm6h4dIQIBIixYigIhOZ9m11Q+3JVm1eg57e3uLz1VYlCxZ\nUlq2bKm1nuSHoKAg6dy5s8ybN09iY2MlLi5OYmJi5Pr169k+Xq/XS3JyspQrV04yMjIkMTFR0tLS\n5K7BIFdTU7M8tqqTk9Rycnpskg8RJQERMb4HRDlGSUCOH7duAnLr1i0RuX9TKFWqlMXnLMoJiIjI\nc889J1WqVJFhw4ZJVFSUREdHy2+//SYXLlzQfrZt26Y93tnZWTIzMwWAfPjhh9K5c2frFGTIELn+\n9dey/t9/ZdDNm0X29TRWjx49ZPPmzbJ9+/Y8E5Du3UVGjhR58UXrXNvuv++XwMBAadasmdHHffvt\nt/Laa6/J5cuX5ddff5U333xTRERSH/oeI8s5OTnJpEmTJDAwUObNm/dEJiAijK8Kkyc9vipUCUjL\nli3Fzs5OTpw4IampqeLg4PDIY+LilKw8JkakaVORwYNFjh6932VoLrV7T82qFy5cKM2aNZPn1b56\nE86hVpCilKEXpE6dOkmPHj3k8uXLEhYWpv3/kiVLSp06daRWrVri4eGhtRgvX75cli1bJjNnzpRJ\nkyaJiHLDTgbkjojcy8xUfjIy5HJqqrS3QlBcmDRr1kwMBoNER0dLcnKyODk55XnMSy+J1Kyp3Fis\nSU0W1LpqjSC3qCcgIiJ169aVunXrZmlpTE1NlXPnzmVpbYqOjpaYmBgZN26clClTRsaPH2+9QlSq\nJIPKlJFdu3aJfsMGGTdunPXOXQj5+PiITqeTiIiIXOtFXJzIvn0i3t4iPj4iJ04oo9X++9o2ya1b\ntyQoKEji4uKkQoUK0qNHDwkKCpJRo0blOnzwwoULMnHiRPnuvzFgnp6e4uTkpAXF27dvN70wlKdX\nX31V3n//fdm7d6/s379f2rZtW9BFshnGV0+mQh1f5fssExOpk3COHTuGW7du4fvvv8ecOXMwZMgf\nKF/+/n5bIkBwsLJ75aefAt27K8vKm+PixYvaBlbNmzfHgQMHtEmww4cPN3olnps3byImJgZr166F\niKBbt27mFegJtWfPHixduhQRERG57iQ8c+ZMiAjmzJljw9IVLo0aNYKI4OjRo7k+7s4dZUfwMWPu\nb8BnTb1799Ymucl/qwBZ6ptvvoGIoG/fvlYoYeF369atLBuBWdO6desg8uRsbOfl5QURwY4dO7T/\nl5aWhr177+CNN4AGDe7fP955B9i6VZloW7Gicj8xdsn8zMxMhIWFaRN77e3ttUUYRARPPfUUdu3a\n9chxqampCAkJ0XZkL1WqFAYPHqytgmZvb4/AwMAnbnlxW3r77be1lfueNIyvnlyFMb4qdAnI8OHD\nISJwc3ODTqfTvtDbtPkcIoCbG9CrF/Dee0rF2LkT6NxZqTAODspNxdgVHS9fvoxXXnlFW3O5ePHi\nEBFUr14dY8eO1W4SFStWRFhYmFHn/P3339GlSxe4uLhgyJAhFrwSlJMFCxZARPDGG28UdFEKjLq7\n6ieffKL9v4SEBOzbdwCLFwODBwMNGwJ2dkrdGDr0/gZ8cXHAyZPWKYe66s9PP/0EACbv/K3KzMxE\nYmIiAGDlypUQEYwaNco6hXyCJSUlaTvKm7s7d1ESFBQEEcHo0aMRHh4OPz8/lC5dGh07ztUCqzJl\ngIEDlcRj2zagb9/7QVft2sBXX+W+D2JsbCy8vb21e1P79u0RFRWFjIwMhIeHo2bNmtrfvL29cfz4\nce1YddlqnU4HHx+fLLt5P/fccyZv4kmmi4+Ph6urK3Q6HTp27IihQ4di8uTJeO+99/Dpp59i48aN\n2Lt3L6Kjo/H3338XuZX4csP4ivJiy/iq0CUge/bswbBhwyCi7AL57LPPYvLkydi06XiOSybevg2M\nGwfodEpFadQIOHr0bI7XuHbtmraDrojAzs4Ovr6+2LZtm7Zsmfy3DJm65JvaYvLHH9nvpXD58mWM\nHz9eq2xVq1Z9rL64CpPwffvw8ooV+CCbFsYnxYcffqj1OAwaNAgNGjSAnZ0dRATly/+tBVQODoCX\nF+Dvr2zA99prQIcOgMEABAaav/3CpUuXEBAQAEdHR5QrVw7dunUza8+OzMxMbNmyBV5eXvD39wcA\nLFy48IlPMK1p5MiRWqAcHh6OmJiYQr9+vrnUvQIeDK6U1tKBePNNYO/e7Ff+3bwZ8PBQ6oyjI9C1\n65hHltJNSkpCUFCQtut6mTJlEBoa+shyyykpKQgNDdV6NXQ6HXx9fXH+/HkAwPjx4zFs2DDtXuHu\n7m50AEbW8eWXX6Jfv345Lo/68E94eDj+NWOTzMKG8RXlxZbxVaFLQADlA3zs2DGkmrh+/f79Sqvv\n009vhZ2dHfz9/ZGQkKD9PT4+HkFBQdoummor1KlTp7THZGRkIDQ0FK6urhARlC5dGoMHD9ZaEStX\nrpxla/ucKtvZszlXULLMzuvX4RUZiekXLhR0UQrMkSNH4OnpmeUm6ejoiBYtWmD69N/w6afKVhAP\nJxjJycCECfdvJk89BTzw8c9TXFwcRo0aBXt7e4gI9Hq9dlOoUqUKvvvuO6POk5GRgS+//FIbSiai\nbHx06dIl7Sb17rvvmvCKUE6OHDmiJSHqj729PTw8PODr64ugoCCEh4drrfhFWUZGBnbt2oV69eqh\nTZs2CA4ONrpXIT1d2UC+V69D2uvUtWtXnDx5Elu3btX2kdDpdPDz80N8fHyu57t+/TreeOMNODk5\nQUTg4OCAXr16aUO19Ho9Jk2ahNu3b1vjqZOJzp49i23btmHNmjVYtGgRpk+fjtGjR6NPnz5o164d\nGjZsiAoVKkCn08HBwQGTJ08u6CJbBeMryo0t46tCmYBY4t49YPbshVqAVKNGDWzatAnBwcHahzy7\nrvGHXbhwAc8//7z2+M6dO6Nbt2746KOPACg3l+wq20lrjW2hHP186xa8IiMx8dy5gi5Kgbp37x7G\njRuHFStW4MSJEybdUPbsAWrVUpIQJydg0aL1ubaKX7hwAf7+/jAYDFluBDExMYiLi0OHDh20uuLr\n65tjb4g6TKVBgwba46tVq4a5c+diypQp2o2mb9++uHXrlqkvCeUgNjYW06dPR69evVCrVq1HegjU\nH1dXV7Rs2RKjRo0yOUApTCxpHU1KSsL8+fO13c0ffK2aN29u9Aagqj///BP+/v7Q6/Vwc3Mz+zxU\nME6ePAmdTgcXF5fHohfEEoyvHn+2jK8euwRE9csvv2gTEtUuc3Xi0rFjx4w+T3h4OMqWLQsRgYuL\nC2bPnp3l5mRMZSPrOnHnDrwiIzHyCRjTnp+SkpRhWB07noCIwNPTEydOnMjymLi4uGwTj4fnE6gt\nW+q4Xnd3d2zevDnLYyIjI1G3bl2t3tSqVQsLFizAhAkTtMRDvdFERkbm+/N/kqWkpCAqKgphYWEI\nDAyEj49PlsTE3d29oItY4K5fv47AwEC0bdsWLVq0QHBwsEWJzalTp3DgwAGEh4cX+Z6mJ02PHj0g\nIpg1a1ZBF6VQYHz1+LJlfPXYJiCAsvrJ+++/j5CQEHTt2hX79+836zxXr15F3759taELD3bP57UK\nEVlfbFISvCIj0T86uqCL8ljYvXuPNnHWwcEBc+fOxa+//go/Pz8t8bC3t4efnx9iY2NzPdf58+fR\nvn37bHtD4uPjUbx4cdSsWRPz5s3DhAkTtOEpdnZ28PHx4Y2mgN24cQP79u17JHl8kmVkZDwyz4Oe\nLIcPH4aIoGTJkuyZ/Q/jq8eTLeMrHQAIGWXjxo1y5coVOXDggEyYMOGJWkO8MEnKzJTVf/0lxQ0G\n6Ve+/GO1yWBBSUxMlMDAQPnkk08EgNjZ2UlmZqY4ODjIiBEjZPr06VK9enWjzpWZmSkrV66UyZMn\nS1JSkri7u0toaKj07NlTtm3bJt9//72sWLFCUlJSxM7OTl566SWZPXu2NLTmDolERFb03HPPyU8/\n/STBwcESGBhY0MV57DC+KhxsGV8xAaEiae0//0gNR8fHbpPBgrZ7925ZsGCBuLq6iouLiwQFBUnt\n2rXNOldMTIyMGDFCjhw5IjqdTho1aiRnzpyR9PR0MRgMMmjQIHnrrbekXr16Vn4WRETWFRERIV26\ndJGKFSvKxYsXxdnZuaCLRJQvbBVfMQGhImnHjRvSrUyZgi4G5eHB3pCmTZtKZGSkDBgwQGbOnMnE\ng4iKlGeffVYOHTokS5culQkTJhR0cYjyha3iKyYgRJTv4uLixMHBQdLT06VmzZoFXRwiIpNt2rRJ\n+vbtK1WrVpXz58+Lg4NDQReJqMhiAkJERESUBwDi6ekpUVFRsnr1ahk+fHhBF4moyOLsXSIiIqI8\n6HQ6mTZtmri4uMiNGzcKujhERRp7QIiIiIiMkJ6eLrdv35ayZcsWdFGIijQmIEREREREZDMcgkVE\nRERERDbDBISIiIiIiGyGCQgREREREdkMExAiIiIiIrIZJiBERERERGQzTECIiIiIiMhmmIAQERER\nEZHNMAEhIiIiIiKbYQJCREREREQ2wwSEiIiIiIhshgkIERERERHZDBMQIiIiIiKyGSYgRERERERk\nM0xAiIiIiIjIZpiAEBERERGRzTABISIiIiIim2ECQkRERERENsMEhIiIiIiIbIYJCBERERER2QwT\nECIiIiIishkmIEREREREZDNMQIiIiIiIyGaYgBARERERkc0wASEiIiIiIpthAkJERERERDbDBISI\niIiIiGyGCQgREREREdkMExAiIiIiIrIZJiBERERERGQzTECIiIiIiMhmmIAQEREREZHNMAEhIiIi\nIiKbYQJCREREREQ2wwSEiIiIiIhshgkIERERERHZDBMQIiIiIiKyGSYgRERERERkM0xAiIiIiIjI\nZpiAEBERERGRzTABISIiIiIim2ECQkRERERENsMEhIiIiIiIbIYJCBERERER2QwTECIiIiIishkm\nIEREREREZDNMQIiIiIiIyGaYgBARERERkc0wASEiIiIiIpthAkJERERERDbDBISIiIiIiGyGCQgR\nEREREdkMExAiIiIiIrIZJiBERERERGQzTECIiIiIiMhmmIAQEREREZHNMAEhIiIiIiKbYQJCRERE\nREQ2wwSEiIiIiIhshgkIERERERHZDBMQIiIiIiKyGSYgRERERERkM0xAiIiIiIjIZpiAEBERERGR\nzTABISIiIiIim2ECQkRERERENsMEhIiIiIiIbIYJCBERERER2QwTECIiIiIishkmIEREREREZDNM\nQIiIiIiIyGaYgBARERERkc0wASEiIiIiIpthAkJERERERDbDBISIiIiIiGyGCQgREREREdkMExAi\nIiIiIrKZ/wcS/l/Gocu/agAAAABJRU5ErkJggg==\n",
"text/plain": [
"<PIL.Image.Image image mode=RGBA size=800x600 at 0x10BA95780>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import pickle\n",
"d = pickle.load(open('/Users/glandrum/rdk/benchmarking_platform/compounds/ChEMBL_II/Target_no_11631.pkl','rb'))\n",
"ms = [Chem.MolFromSmiles(y) for x,y in d[30194]]\n",
"Draw.MolsToGridImage(ms[:12],molsPerRow=4)"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"scaffold: [#6,#7]-[#6,#8]-[#6]1:[#6]:[#6]:[#6](:[#6]:[#6]:1)-[#6]1:[#7,#8]:[#7,#8]:[#6](:[#8,#7]:1)-[#6]1:[#6]:[#6]:[#6]:[#6]:[#6]:1\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElE\nQVR4nO2dd0AU1/bHv7tLLwIWFBEVRQUsGMGowRgLUXmCmih2jDGKJQnmJU83zZBn8hSj5ofPGIPG\ngiQWbAFRjCgWVCygYokVIYJGQHpny/z+uAnhscsAO7Ozxfv5C/de5p51h++eOfecc0UMw4BCoVAo\nmiLWtQEUCoVi2FAZpVAoFE5QGaVQKBROUBmlUCgUTlAZpVAoFE5QGaVQKBROUBmlUCgUTlAZpVAo\nFE5QGaVQKBROUBmlUCgUTlAZpRgtyckIDESfPhg/Hhcu6NoaivFCZZRinFy8iBEj4OmJzz6DRIIR\nI5CerhtLqJobPSLamoRilEyeDIUChw4BgFIJZ2fMnIm1a4U24+JFDB2Kjz5C//6IicHRo7h8GV5e\nQptB0SomujaAQtEKnp546aU/fxaLYWkJsRiTJkEigb09xGLY20MiQatWMDWFjQ3MzWFlBQsLWFrC\nygrm5rCxgakpWrX681dEIk3MWLsWgYFYvRoApk6FszOio6mMGhtURinGyYoVf/+8fj3y8zF5MgYN\n0vyCIpF65bW0hIUFrK1hZgZbW5iYwM4OYjFGjMDLL6tRc4mE0/ui6CFURilGzqFDWLoU8fF46SXs\n3w+lEkVFUCpRUgK5HGVlqK1FRQWqq1FVhcpK1NSgvBwyGUpLoVCguBhKJYqLwTAoKgKA58+btW5Y\nGF5+WY2aBwdr5W1SdAiNjVKMmcJCdO+OL77AP//J9VIKBUpLIZOhvBw1NaisRFUVqqsbVd7hw+Hv\n//evHzqEqVMRH4/Ro7laQtE3qIxSjJlTpzByJMrKYGOjSzN4VHOKHkITnijGjKcnkpJ0rKEA0tNR\nXIz583VsBkVLUBmlGDPZ2ThxAjU1OjZDT9ScoiWojFKMkG3bcPgwZDJcu4aVK3Uvo3qi5hQtQWOj\nFGOjthbOznj+HNev60uG5pYtCAlBSQlatdK1KRQtQBOeKMbG0aN4/hy9ezetoWvXoqoKZmZ/Z9oH\nBDR89M7NhZnZn2n5GjN/PgYMwKRJ6NoVW7Zofh2KfkJllGJsREcDwJw5Tc/csAGPH//PK1lZDWV0\n8GBkZf35M6luunYNXbr8z5z33kOHDvj8c7a1rK1x4gTatUNkJMQ0lmZcUBmlGBVFRThyBGIxpk1r\nevLSpcjN/Tvfs6QEdnYN59jZwcHh7xTRykqYmzecc/AgfHyaWMvdHd27IyMDaWkYOLDZ74diCFAZ\npRgVe/eipgajR6NTp6Ynv/de03OuX//754oK1NbC3r7hnE2b4OjY9KXGjMH33yMhgcqosUGfLihG\nBXmi11LBpbU1HBzU9CiZMAFDhjT966SoKSGBf8MouoXu1FOMh4wM9OgBKys8e6aPSZoVFWjbFrW1\n+OOPZnmvFEOBeqMU4yE6GgyDSZPUaKhM1tyWItrD2hrDhkGpxIkTOraEwi9URinGw+7dQCNP9EeP\nomNHfPSRcMaofcyjz/VGCZVRipFw/jzu30fHjhgxQs2ovz/27cNrrwlhyZkzGDQI77+v3gwAx45B\nqRTCEoowUBmlGAlkc2nWLPV9kc3MMGECxo8XwhIrK1y+jCNH1Az16oXu3fH8OVJThbCEIgxURinG\nQG0t9u8HgFmzdG0K4O2N9u2RlYW7d9WMjh0L0Od644LKKMUYOHwYBQV46SX07atrUwCx+M/ezGq1\nkoZHjQ8qoxRjQKvpohrAopUjR8LSEleuIC9PYKMo2oLKKMXgKSzEsWOQSDB9uprRq1dRXCy0SWPG\nQCJBVhYUioZDlpZ/pj0lJgptFUVLUBk1Hh4+fBgbG1tSUqJrQ4Rm9+4/C0A7dGg4xDCYOhVOTrhx\nQ1CTWrfG/fu4f1/9fhd9rjcyqIwaNtXV1SdOnPj44499fHx69Ogxbdo0Dw+Pr7/+Wtd2CQrLE31K\nCh4+ROvW6N1bYKPQrVujQ0RGf/1Vja9KMURoaxKDJCMjIyEhISEh4fTp05WVleTFNm3a9OzZ89Kl\nS8uXL2cYZvny5bo1Uhju378vFicMHTp9wgQ19ZXsWVC6omdPuLnh4UOkpmLQIF1bQ+EOQzEQqqqq\nEhMTpVKpt7d3/U/Q09NTKpUmJibW1tYyDLNr1y4TExMAy5Yt07XJQvD5558DePvtt1WHamqYNm0Y\ngLlxQ3i7muC99xiACQvTtR0UPqDeqL6TmZmZmJh44sSJY8eOlZWVkRdbt249atQoPz+/cePGOTs7\n158/ffp0ExOTmTNnfvPNN+Xl5d99951ItSWRscAwzK5duwAEq3ukj49HQQH699eLLKgGBARUX736\ny4MHd4B/69oWCmd0reMUNchksuTk5CYdTxbi4+MtLCwAhISEKBQKYcwWnjNnzgBwdnaWy+WqowsX\nbho69MH69TXCG1bH7dvMnj1qXq+srLSyshKLxbm5uYIbReEZKqN6xKNHjyIjI4OCglrVO/nMxsYm\nICAgMjIyJyenRVdLSEiwtLQEMHPmTJlMpiWbdcu8efMAfPLJJ2pHFyxY4ODg8PTpU4GtqqOwkDEx\nYRwdC8rKqlVHx44dCyA6Olp4wyj8QmVUx3B3PFk4c+aMra0tgKlTp3K5jn5SVVVlb28P4ObNm43N\n0fm7njJlkUQi+eWXX1SH1q9fD2DGjBnCW0XhFyqjuiEzM5NHx5OFc+fOkSUCAwOrq9X4RIbL3r17\nAfj4+OjaEDZI8tmCBQtUhzIyMgC0bt1abUTCKMnLy8vJyYmNjdX51xu/UBkVDq06niykpqa2adMG\ngL+/f2VlpTaW0AkBAQEA1q9fr2tD2EhLSwPg4uKidtTNzQ1ASkqKwFbphJs3b7q6uvbv3x+Ag4ND\nSEhIcnKyUqnUtV08QGVU66h1PK2trYnjmZ2dLYAN165da9euHYDhw4eXlZUJsKK2ycvLMzU1NTEx\nefbsma5tYUOpVDo5OQG4ffu26uj7778P4IsvvhDeMIE5cuQIuf/d3Nw8PT3r/hB69er11VdfZWZm\n6tpATlAZ1Qr1Hc/6+UZ1jmdNDdft42fPnu3YseP06dP1X2R5PLxz507Hjh0BvPrqq6WlpRxX1zkk\nsBgQEKBrQ5rmrbfeArB27VrVoaNHjwIYOHCg8FYJSWRkJMllnjJlCnkeunXrllQq7VCvetfb2zsi\nIuL58+e6NlYTqIzySW5u7qZNm8aPH29tbV13f9jb2wcFBW3dupX7lrFCoUhNTQ0PD/f19RWLxQCm\nTZtWN3ro0KG+ffs+efKksV+/d+9ep06dSDyxoKCAozG6xcfHB8DevXtVh86dOzdmzJjY2FjhrVLL\n3r17nZ2dv/32W9Whqqoqkvak5z61xshksnfffReASCSSSqUNHuEVCkViYmJwcLDNX4dnmZubBwQE\nxMTEGFbwlMoob+Tn569Zs6ZOPbt16xYaGsqL41lHfHx83fWtrKzGjRu3Y8cOMiSXy0nUyd3dnWWH\nKisrq3v37gAGDBiQn5/Pl2ECc+/ePQCtWrVSG+qdP38+gI8//lh4w9TCvoPk7+8PYOfOnYLZIxgF\nBQUjR44EYGFh8fPPP7PMrKysjImJCQgIkPxVtNu6dWsSPBXMWi5QGeWNqKgospmwdetWFpewmcjl\n8lu3bjV4saKiom/fvkuWLDl27FhVVVWD0cLCwoEDBwLo2rVrRkZGY1d+/Phxjx49SIRBhzmVXPjk\nk08AzJs3T3WoOVlQesV///tfANOnT9e1ITzz4MEDd3d3AE5OTpcvX27mb+Xk5ERERBCHgODu7h4W\nFsZyP+sDVEZ54+233wbwf//3f1wukp+fHxMTExIS4uTkZGpqWlxc3KJfLyoqGjx4MIDOnTs/ePCg\nsWnPnj3r06cPCfDzmFwlDEqlsmvXrgDOnDmjOhoTE0N8beEN0wyjTHtKTEx0cHAA4OXl9fvvv2tw\nhQbBU7FY7OvrGxkZWVJSwru13KEyyhuurq4Arl27pvEVxo0bRyKeBDc3t/T09JZepLi42NfX19bW\nlj2NJjc318vLi7iujx490thm4UlKSgLQpUsXtUWufn5+AFatWiW8YRpDHg4uXLiga0P4YfPmzaam\npgAmTZpUUVHB5VJyuZwET+s2GywsLPQweEpllB8eP35MfAouBeyzZs2ysLDw8/MLDw9Xmx/TTMrK\nyq5cudLktKKiokGDBjXpuuobxOv//PPPVYeqq6vNzc0BdO/enct/oMCEhoYCIO0NDRq5XC6VSone\nhYaG8tjMoaSkJCoqys/Pry7vRa+Cp1RG+WHHjh0AJk6cyD5NoVBcunQpLCzs5Zdf/uGHHxqMPn36\nVOD0+OLi4ldeeQVAhw4dVEOxekhVVZWdnR2AO3fuqJ1w584dktNuY2OzR21TEB2xadOmbt26qTUp\nISEBel+O1SSlpaWkIMLc3DwqKkpLq2RnZ0dERJAHKYKHh0dYWJhun6iojPJDk4HRtLS0mTNntm3b\ntu7jb1JzhaG8vHzUqFEA2rdvf0MPG3PW49atWzNmzADQr18/lmkVFRWz/jpnWR+arj558mTcuHHE\nnkmTJqlOMIK0p4yMDJJU37ZtW7Ux68bIycnRrG8OCZ62b9++QfBUJznRVEb5ocnA6MmTJ8nn7erq\nGhISEhMTI+TnvWvXLpabu6KiYvTo0QAcHByav6kqDBUVFYmJiaGhoV26dCH/gR07drSwsFCbMVqf\nyMhIMzMzjjt+3ImJiSGVuPb29ps2bWps2j/+8Q8A2nPitMq5c+ccHR0B9O3bl70eSbWC7tVXX+Xy\neF4XPLWysqoLngYFBcXFxQnZ1YzKKA80JzBaU1Ozfv36e/fuCWkYISUlxcTExMrKKjExsbE5NTU1\nEydOJH/t+lDiffv27bVr144aNcrMzKzOf+/QoUNwcDBJRRSJRP/617/Y/1R+++03wQxWS+Jfh38G\nBASw55Zt2LAB/1tMYShs3bqVfEZjx45lTyzZsmVL+/bt79+/X/dKRUVF/cJQT0/PVatWaRbXKi4u\nbhA8dXJyCg0NFSZ4SmWUB5oZGNUVSqWSVJKYm5ur7dhGqK2tnTx5MgBra+ukpCQhLSSoOp4AJBKJ\nt7e3VCpNTk6u+5YiniaAYcOG/fHHH8Kb2nyCg4O3b9/e5LQff/xRLBaTgjeBPSmNUSqVYWFhdRtK\nLAlbcrn8gw8+IDMjIiIajNZ/PO/YsSPHxK/Hjx+Hh4f37NmzvjqHhYVptWyfyigPsARGMzMz9SHF\nXalULlmyBICZmdn+/fsbmyaXy0kBuJWV1fHjx4WxLSMjIzIyMiAggGyyE9q1axcUFBQVFdVY0eq5\nc+dIyw9nZ+fz5883c60rV66Ul5fzZzsPFBcXk/uHuNjkB2dn56VLl+pzqLqsrGzChAkATExMNm7c\nyD4zMDCQfIuzfKPU1tYePnxYtdipqqpKM2FNSUlZvHgxiaiQ7+PRo0dHR0dr4yuKyigPsARG33nn\nHQBbtmwR3qoGKJXKf/3rX0ShWMKyCoVi7ty55KbXXlk6cTylUikpdGFxPFnIy8sjD/gmJibh4eFN\nzn/48KGDg4NXv37VWsvuysrKatH848ePu7i4kIheeHj4o0ePwsPDSRppnScVHh6uD9/E9cnJyRkw\nYAAJZLE/uMhkMlKS5OjoeO7cOQ3W+vbbb9u0acNX8LRjx44NuvnwApVRrrAHRkkBe1pamvCGqeWr\nr766ePEi+xylUkkyGc3MzA4cOMDj6modz7Zt2wYFBUVGRmrweC6TyeoSFWfOnMme7H3nzh13d/cT\nr73G2NoyjbvkmqEoLAyeNcva2vrhw4fNmV9RUREaGkrcz8GDBzfI30pNTQ0NDSW9DVFvG1ofmhym\npKSQ4qIePXrcvXu3yflbt27t06ePxglJb775Zt2t0qdPn2+++Uaz06sKCgr+85//EDdCM0tYoDLK\nFZbAKFFYOzs7g6vzUyqVH374IfEQOR4WVFlZqep4isXiFjme7OzevZtUuXh5ebGrWFlpqXLyZAZg\nRCLms88Yvj6XhASmU6ddw4ZZW1uzRJ/rSE5OJsmtxAlt7Paorq6Oi4sLCgqq22eztLTUbfB0z549\n5ICv119/vaioqJm/xfHYBRI8JckAAJr0AxpDqVSSfgu8HyNIZZQrLIFR0qxk/PjxwlvFC2QDQSKR\nNGeTpAF1jic5oJS748nOnTt3yJ5vq1atDh48yDZVqWQiIhhTUwZghg9nOKZqFhczc+YwAAPIx4xp\n0hWtrKyUSqWk5NfLy6uZpcNFRUUNtqGdnZ1DQ0OvXr3KyfiWQDaUiAEhISHC12LW1NTExsa+++67\nqkPN/yb29fUFcOLECV5NozLKGZbAKAkyqm00qT9cv36dZTQ8PJxsfWzYsKHJS9U5nh4eHlpyPFko\nLS0lmQaktWUTayUlMY6ODMBwyTH69VfGxYUBGEtLJjy8Sd/2woULZAfZ1NRUKpVqoERZWVnh4eHE\nk60fPNV2ukJVVRUpfJBIJOyntpSVlWnsLWpGZmYm+VJJTU1tcvLChQsbc3q4QGWUE+yB0RUrVnh6\negrpMrSUnTt3isViqVTKMmfNmjUikUgkEjV287E7nkJujyiVyoiICNJo3d/fv4nW1NnZzIQJjMbm\nKZXMyy8zADNkCNNUOnBVVZVUKiXNNPv27cs9Vk6Cp3VFcRKJxM/PLyoqSht5CE+ePCENGG1tbePj\n41lm5uTkeHt729raCplj8O2339bddV5eXmvXrmW55TZu3Ahg7ty5/NpAZZQTep4x2iR79uwhzXiW\nLl3KMm3Tpk3kOXTFihXkFd06nuycPn2aJCG6uLhcunRJiyvdvs18802TTujFixdJXNjExEQqlfLY\nxls1eNqqVavg4OC4uDi+wvHXr1/v3LkzgO7du7OXM1y9epWcreDm5tZYxwMt0WBH7pNPPmls5tmz\nZ6GFU1uojHKClx6juiUuLo7smy9cuJBF+7Zs2UKUdNy4cWPGjKnveLZr127mzJk///yz/rTT//33\n34kDZWlpeYy177pWqa2tDQsLI05o7969m9N2SzMKCwsjIyN9fX3rgqedOnUKDQ1lj9g0yf79+0me\n0NChQ/Py8lhmxsTEkJl+fn6FhYVcFtWYmpqagwcPTpw4kaW5V1FRkUgksrKy4vdrnsooJ7j3GNUH\n6pSUvZHH7t27xWIx2TCtczwTExP1s+Smuro6NDS0u4ODzMWFCQ5mmlNiqFAwM2cyv/6qfjQ5mWnJ\ngWvp6ekkZZI4oRx3q5vJ3bt3w8LCSJpd/eBpS5ueKJXK8PBw8sX5zjvvsHvQOTk55Gt10aJF+nkz\n1Ie4zPx2hqQyqjm89BjVE06dOuXs7MwepJfJZNbW1iKRaOPGjfrjeLLzPDqasbBgAGbQIObx4yZm\nR0UxAPPSS4zqB7ptGyMWN3c/qqZGuXx5YI8eANzd3QXeciGQ59z6NTzND55WV1fPnj2b/FZz6hoY\nhomOjm7OJqQ+QA6/aiKdo4XwI6Nnz54NCAjo3bt3YGBg8yvzDB1DD4w2QPVwpwZcuHCBeDfC2MMb\nV68yrq4MwLRty7BXuCoUzMqVjFo/5ckTxtGR+eIL5n/PtlTD9euMlxcDlPbt+9FHHwncQLYBdcFT\nEgEHYGdnFxwcnJiYqGzkjeTn57/66qsAbGxs9Od01ZYSFxcXHBysdkNs2bJlAP7973/zuBwPMpqS\nkiKRSJYtW7Zr166JEyeamZlxjMgYCiyB0RMnTnz99dcG1IC9OaxcuRLA4sWLdW1Iy3n+nBk7lgEY\niYQJC1PjbDaHJhsbymRMeDhjZsYAjKsro4WiQ40pKCggwdO6h30XFxepVNqg5diNGzfIOVedOnXS\n5wyTJiEFSx999JHqUHR0NIDJkyfzuBwPMjpp0qQ6j0yhUHTo0EGt9cZHt27dAKi920gp/bp164S3\nigsymYzlaX3MmDFo5Gh4A0CpZMLDGbGYAZjAQKbZFTjN5dYtxsfnz/qokBBGD6o21fLbb7+FhYWR\nW5fg7e0dERGRm5t79OjRVq1aARgyZAh7LPXAgQP79u0TzGYNiIuLAzB69GjVoevXrwPo1asXj8vx\nIKPLly+vH2hwdXXVh5bj2oYERh0cHAyilL6ZXLhwQSQSvfnmm6pDMpnM1tZWJBIZbod2hmGY+HjG\nwYEBmB49GL5OYJbLmfBwxtycAZiuXZmTJ/m5rDZRKBRJSUlz5syxtbUlYmpqakp2+WfPns2yG6ZU\nKv/zn/+IRCJLS0t9PgkxMzMTgJOTk+pQdXW1qampRCLhMd7C8xZTRESEjY2NoRwRzgVS6DlhwgTV\noezsbBhmKT15FFJbb2eogVFVHjxg+vVjAMbGhuF+WNPDh8yrr/7thOriBAsuVFVVxcTEBAQEmJiY\ntGnTZuHChSyTq6urSR9FkUgUFhYmlI2aoFQqiWetNk+L1A3z6OX8fZwvdw4dOrR06dIDBw6QM9CN\nmzNnzgB47bXXVIdOnTpFhkjCoAHB8qZOnz4NYPjw4cJapAXc3HDpEubORXk5pk3DggWQyTS5DsNg\n82Z4eSE5GU5OiItDZCT+cu4MBXLkxuHDh6dPn15QUNC7d+/GZhYUFIwePToqKsrGxubQoUNffvml\ngGa2GJFIRN7L7du3VUf79u0L4ObNm3wtx5uMFhYWzp07d/Xq1eRUH6OHRVZYxEifkcvlKSkpIpFo\n2LBhqqMG+qbUY2GBrVsRGQkzM2zejJEj8ccfLbtCZiZGjMCCBaioQFAQbt9GQIB2bBUIHx8fsCpL\ndnZ2Wlpaly5dzp8/Txo26zlhYWEJCQkvvfSS6hDvMsrbQ31SUhKA+v0Qhck31gkkMGpvb6/2sT09\nPf2bb77RybFLXDh//jyA3r17qw4ZSWBUlStXmC5dGIBxdGSaeW6KUslERjI2NgzAdOjAGGxKUAPI\n3++QIUNY5pw8edI4boDY2Fg0sgGlGbzJ6LNnz+r3wU5ISOjcubOxBklZAqOGy5EjR9zc3NTmMxlP\nYFSVvDzGz48BGBMTpslU88xMZuRI0haPCQpqUVGTnpOfnw/A1ta2sXxSY+LRo0doZANKM3iT0StX\nrnz66ad1HigpFXB0dExPT+drCf3BIDrgaYbaYj4DzhhtDjIZ8+GHjEjEAMyMGerLRokTamvLAEz7\n9gyvNTB6Aulpr/bot++///4f//iHvh2+rTHsG1AawFts9Nq1aytXrqypqSH/PHTo0Pjx4/Py8oYP\nH37p0iW+VtETjGe/RQXSZa4BRhUYVcXEBOvW4dAh2Nnh0SOQjcHkZAQGok8fjB+PCxewYwcWLEBZ\nGaZNw+3beOMNXRvNPywRw7S0tKNHj165ckVwo7SCSCQim/VqN6A0gDcZnT9/PsMwROMBmJub79u3\n78033ywqKhozZgx5KjQOsrOzHz16ZG9v369fP13bIgRyuZzkkxqtjBImTMCFC9i/H2ZmuHgRI0bA\n0xOffQaJBCNGoF8/vP469u3D7t34q1DdyGCRUf73ZHQNv++Iz4SnBpiZmcXExMyaNaukpGT06NEn\nT57U3lpCYrj5TJqRmppaVlbm4eFBmngaM56ecHYGgLVrERiI1asxfToOHEDr1ti9G8ePY/JkXZuo\nRYxPRvfu3evj47Nu3TrVIYORUQASiWTHjh1z5sypqKgYP3788ePHtbqcMLA84f700089e/YkHbaN\nBhLBMHJXtAGenpg9+8+fxWJYWuIF+MpkURby4HXr1i2GYYQ2iwM1NTVpaWlqYxGGJKMAJBLJtm3b\n3n333crKysDAwF9++UXbK2oblsDo6dOnHzx4UBcgNhRkMtm6detSU1PVjhp5YFQtK1b8Hf1cvx75\n+QgO1qlBQuDp6SmRSO7du6d6A7dt27Z9+/YlJSWkQs9QEO6LgZeNqiZRKpUffPABADMzs/18HxEu\nJOwZowZaSt+cjFFtH5qmpxw8yJiaNtrI2eggJ+6pza45evTo5cuXeTwBRQCqq6tNTExMTEzU5rA7\nOTkB4KUzgNa9UQI5EO3zzz+vra2dOnXqzp07hVmXd1gCozk5ORkZGXZ2dl5eXrowTXNY/GsSGHV3\ndyfZMC8WhYWYOxerV+PFKMwDq/vm7+8/cODAukOfDAJzc3M3Nze5XH737l3VUR6f6wWSUcJXX30V\nHh6uUCjmzp27bds2IZfmC1pK/wKRno7iYsyfr2s7hMNAt5JYIO/o1q1bjQ0ZnowCkEqlq1evVigU\n8+bN27Bhg8Crc4dFVu7fv2+IWUEvUCl9S/H0RFISbGx0bYdwGJ+Mrlq16vHjxzNnzlQd4vPNco8L\naMDGjRvJ0eeGVQjEHhhlGCYvL6+Js9H1D5ZCzxc9MHrlCvPpp4zxtoZQ5f79+wBcXFx0bYgQpKWl\nNXbntxShvVHC4sWLf/jhB5FI9OGHH65YsUInNmhAk4/t7dq1s7OzE9Yorjg7O3/99deLFi1SHXqh\nA6MArl3DypUwtLwLLnTv3t3a2jo7O7u4uFjXtmgdkplw//597qk1upFRACEhITt37jQxMQkLC/v4\n4491ZUaLaPIJt6Sk5JVXXtm+fbuARnGlc+fOn3322Xvvvac69EIHRgHMnw+GwV+FeS8CYrHYw8MD\njQQTV61a5erqunv3bsHt0goWFhYsG1AtQmcyCmDmzJk//fSTqanp6tWrly5dqkNLmkQul585c+bw\n4cNglZWDBw9evnx53rx5mzdvFs44rfFCB0ZfVFgihjU1NVlZWcYUOeUrPKpLGQUwderUAwcOmJub\nr127dtGiRUqlUrf2NCAvL2/fvn2zZ89u167d8OHDS0pKzMzM9uzZ09j8t99++7vvvmMYZuHChRER\nEUKayjtyufz8+fOGuGlG4YLxlYQCUCqVZWVlqq/z9o64h1e5c/ToUQsLCwDz589Xe0KckMhksjNn\nznz88ccN0j979+5NjqwBsHTpUpYr/PDDD2KxGMCKFSsEM5t3UlJSAHh4eOjaEIqgJCYmAhg6dKjq\n0L179wB06dJFcKM48fPPP1tZWak9ZurgwYMA/P39OS6hFzLKMMypU6dsbM9Ebt8AAA4sSURBVGwA\nTJ8+XW3LS22Tm5sbExMTHBxsb29fJ51WVlZ+fn4RERFZWVlk2uHDh83NzQEsXLiQRfF//PFHsVgc\nEhIilPn8s2rVKgCLFi3StSEUQXn27BkAOzs71f7NcrncyspKJBIV8X48tTYhTZF8fX1Vhx48eACg\nU6dOHJfQFxllGObs2bPkuNcpU6bU1tYKsKJcLk9NTQ0LC/P29ianyxK6desWEhISFxentobs6NGj\nlpaWAObNm8eipMnJyTr3rNmpra318fF599131eZvjR07FsAe7mdnUgwNR0dHAI8fP1Yd8vHxEYvF\nhtW/mTT2V/vFoFAorK2tAXDMU9QjGWUY5sqVK61btwYwbty4qqoqLa2Sl5dHHE8HBwcWx5OF06dP\n69Z35oUmM0YBPHnyRHjDKLpl5MiRAI4cOaI6lJmZWV5eLrxJHCE9Hn///XfVoYEDBwI4e/Ysl+vr\nl4wyDJOWlta2bVsAY8eOrVR7nINGEMczPDzc19e3+Y4nC8nJyaRHdVBQkDC+M++wHA1CA6MvMkuW\nLAEQ3uTJVIaDn58fgPj4eNUhciDQxo0buVxfzYkRumXAgAFnzpzx8/M7duyYv7//4cOHbTmc/Z2f\nn3/69OnDhw/Hx8cXFRWRF62srF555RU/P78JEya4u7trduWhQ4eePHlyzJgx+/btq6ys3L9/P9kl\nY7Fkw4YNYWFh+lNxT0vpKWox3B35xujbt+/ly5efP3+udgjc3ywXDdYed+/edXZ2BjB06NCSkpIW\n/S7vjicLV69ebabv7OvrC2DatGl6EgRgPzOZBkZfZC5fvgygX79+ujaEN1j+Nlk2oJqPnsoowzCP\nHj1ydXUdPHhwaWkpwzBnz54NCAjo3bt3YGDg+fPnVeerjXhaWlr6+fmFh4ffuXNHS3b+9ttvpHHh\nsGHDiKlquXz5Mgn7BgQEaC/s23xIk2YaGKWoUlFRIRaLzczMDKu7qGawbEA1H/2VUYZhsrKyiouL\nGYZJSUmRSCTLli3btWvXxIkTzczMrl+/zvyv40lSNRs4nsIIVp3v7Ovry+I7p6WltWnTBoC/vz+P\nYV+Nefr06ZUrV1Rfp4FRipubG/5qDq9KaWkp+cPUfzIzM5t8+mTZgGomei2jdUyaNGnixInkZ4VC\n0aFDh/nz50+dOpX4dwQrK6tx48Zt3LiRl3bWLSUzM7Nbt24AvL29nz9/3ti069evt2vXDsCmTZuE\nNK9F0IxRyhtvvAFg9+7dqkNhYWEikWj16tXCW9VSzp07165du9mzZ7NPY9mAaiY6LgZtJp6enrP/\nOmJMLBZbWlpaWlru37+/sLCwzvEsKCiIj49fvHixq6ur8BZ27dr11KlTbm5uaWlpfn5+5ElBFS8v\nr9OnT3/66acLFiwQ2EJCaWnpwYMHSce/xqCl9BSWjReSrK7/G1A//vjjyJEj8/PzS0pKmuzhJBaL\n2f8omkBjAdYVERERNjY2N2/e3Lt3b0ZGhq7N+R/++OOP3r17A3B3d9erwGJGRkZERISfnx85BGLt\n2rWNzaSBUQrDMPv27QMQGBioOnTx4kUAXl5ewlvVTORyuVQqJfoWGhrKUgUjl8vff/99AK6urlzi\nbAYmowcPHjQ1Nf1Vj48Yy83NJYcO9uzZMzs7W9fmMAzDhIeH131rmpiYDBs2jGULngZGKQzDkN5x\nXbt2VR0qLy8Xi8Xm5uZ6my5NbngzM7Pt27ezTCstLQ0ICABgbm4eFRXFZUVDktGCggJ7e3v9b5hf\nWFj48ssvA+jSpcvDhw+bnP/kyRO1WUeakZeX1+CVlJQUJyenuXPn7tu3r8mdARoYpTAMI5fLLS0t\nRSKR2i1Tsg1w+/Zt4Q1rDuXl5aNHj2YvTMrIyPD09ATQtm3bM2fOcFzRkGQ0KSkJQFlZma4NaZqi\noqIhQ4YAcHFxuX//PsvMvLw8T0/PXr16cXFdScQzJCTExcVF9ZxkpVLZ/GQOmjFKIQwYMADAhQsX\nVIcmTpzo6urKsYBSh5w/f570Dejbt29mZib3CxqSjD579iwpKUnXVjSX8vJyUpvcvn37mzdvNjat\nsLCQVPV27dpVg1BvRUXFyJEj6x976+TkpHEyCg2MUup46623AERGRqoONXYWmUGwa9cuUnA4duxY\nvtK2DGOnnpCdnX3ixAnuB6cIg7W19eHDh19//fXc3NxRo0alp6ernebg4HD8+PHBgwdnZWWNGDHi\n4cOHLVrFysrq6dOnCoXC29tbKpUmJyc/efJE4/OgyOFLHh4eHTt21OwKFKOBZbNefwqaASiVyu++\n+66ysrLJmQzDfPnllzNmzKiurg4NDY2Pj+ft5DRexFgYyMkcLa0N1S3V1dXjx48HYG9vf/Hixcam\nlZSUkGpRZ2fn3Nxc1Qm3bt1as2bNyJEj09LSGgylp6cXFhZytLOqqurXX38dOnQoaGCUwjAMw/z6\n668AXnvtNZY5KSkpcXFxOtxrKi8vf/PNNwFMnz6dfWZZWdnEiRMBmJiYcGxEooohyaiBUlNTQz5p\nOzs7tWWshIqKCj8/vyVLltR/8dy5cwsWLOjcuXPd1x6/HfUzMzMjIyODgoJa/XVw24IFC27cuMHj\nEhQD5enTpwBat27NMmfChAkAHBwcQkJCkpOTudRTasCTJ098fHyIASdOnGCZmZOT4+3tTd6ONgKD\nVEaFQC6XBwcHA7C2tmb5vKurqxvciHW5Su3atQsKCoqKiuLueFZXVx8/fvyf//xn/e5WIpFowIAB\nn332WXp6OsfrU4wGUnHHEij/9ttvSaI0wd3dXXvNKxpw7do1FxcXAG5ubr/99hvLzJSUFHJCeI8e\nPe7evasNY6iMCoRcLn/77bcBWFlZNT/v9d69e19//XVaWhr373lVxxOAjY1NQEBAZGSknqS4UvQK\n0inx2LFj7NNu3bollUo7dOhgZ2cnWLOIefPmARg1ahS7Y7F3715yVsXrr7+uvbNPqIwKh1KpJMfB\nm5mZHTp0SIAVZTJZcnKyVColTzR1eHp6SqXSxMREvc2gpugDpMJnzZo1zZksk8lUM1JKS0vj4+O1\n0Ryyurp6zZo1LDewUqkMDw8nrTJDQkK0eqtTGRUUpVL5wQcfADA1Nd23b5+WVvnjjz+ioqIaOJ7W\n1tbU8aS0CLKp22RrDxa2b99OIpIkeMqjbexUVVXNmDEDgEQiWb9+vbaXozKqA5YvX04+4B07dvB1\nzfqOZ4Nm1aGhoYmJiS9C70gKv5DKYDc3N43zK3fv3u3h4VH/MSguLo5fI1V58uQJycW2tbXl0rep\n+VAZ1Q1k70gsFm/dupXLdeocz/oZcNTxpPBCaWnpokWLTExMzM3NAwICYmJiNHs0JsFT0tYzISFB\ngyvcuHFj7NixzQluXr9+nWS2dO/enX3riUeojOqM1atXky3ylj50UMeTIhh3794dMWJEXU90R0fH\n0NBQtd2+m6S2tlZtnDQ1NZW9LOqXX34hB/F+8skn7Evs37/fysoKgK+vr2pzCe1BZVSXfP/99yKR\nSCQSrVu3rsnJz549Y3E81Z4qTqHwQk5OTkRERP/+/evnNoWFhXHvVPn8+XMzM7M2bdo0FjyNiIgg\nIj59+nT2wyzqZr7zzjsCexJURnXM5s2byWf/5Zdfqo7K5XLqeFL0h7rcJnIfisViX1/fyMhIjWsL\nL1++3KtXr7obu0+fPt9//z0Zksvlc+bMIauwn/ZcXV1N2rpLJBKdnAtNZVT3/PzzzyYmJgCkUil5\nhTqeFH1GLpcnJiYGBwdbW1uTm9PCwoJ78JR0XVqwYEHd64sXL7a2tmbPDszPzx82bBhJgo6NjdVg\nde5QGdULdu3aRZR08ODBffr0qZ/j2a9fP6lUeurUKZrjSdE3SkpKoqKi/Pz86h6VuOQ21dTUxMbG\n1j9HTyaTsZdF3bhxo2vXrgA6dep09epVDRblBRHDMKDoAfHx8dOmTXN0dMzMzLS2th4yZEhAQMAb\nb7xRv6CeQtFPcnJyDhw4sH379rpOZh4eHlOmTHnrrbe0dzbasWPHpk2bVlJSMmTIkEOHDpFMAJ1A\nZVSPyMvLu3btmrm5ua+vr6mpqa7NoVBazO3bt6Ojo3fs2JGbmwtALBYPGTJk9uzZ06dPJ61s+WL9\n+vUfffSRQqGYNm3atm3bSMWnrqAySqFQeEahUJw6dWrnzp0HDhwgnUAtLCwCAwODg4P9/f1J/Epj\n5HL5kiVLSJbLF198QQ585slwDaEySqFQtEVJSUlsbGx0dPTJkyeJ1Dg5OQUFBQUFBZHmti2lsLBw\n8uTJp06dsra2jo6OfuONN/g2WROojFIoFK2TnZ29a9eubdu23b9/n7zi6ekZFBQ0Z84cskfUHB48\neBAYGHjv3r2OHTvGxsaSZqP6AJVRCoUiHGlpaTt37tyzZ09eXh7qBU9nzJhBSpUaIzExccqUKcXF\nxf3794+LiyPNRvUEKqMUCkVoVIOnlpaWAQEBjQVPN2/e/N5778lkssmTJ0dFRZGKT/2ByiiFQtEZ\nxcXFcXFx9YOnHTt2nDx58ltvvUVOeFYoFB9++OF///tfkUi0bNmylStX1hX46w9URikUiu55+PDh\nTz/99NNPP2VkZJBX+vfvP3Xq1KSkpMTERHNz8x9//HHWrFm6NbIxqIxSKBQ9ggRPd+/enZ+fD8Dc\n3NzBweGXX34ZNGiQrk1rFCqjFApF76ipqTly5EhcXNzUqVP79OmjVxtKqlAZpVAoFE7oXbCWQqFQ\nDAsqoxQKhcIJKqMUCoXCCSqjFAqFwgkqoxQKhcIJKqMUCoXCCSqjFAqFwgkqoxQKhcIJKqMUCoXC\nCSqjFAqFwgkqoxQKhcIJKqMUCoXCCSqjFAqFwon/B2AN8s8KkDIFAAAAAElFTkSuQmCC\n",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x10c023c10>"
]
},
"execution_count": 17,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"params = rdFMCS.MCSParameters()\n",
"params.SetAtomTyper(rdFMCS.AtomCompare.CompareAny)\n",
"params.BondCompareParameters.RingMatchesRingOnly = True\n",
"params.BondCompareParameters.CompleteRingsOnly = True\n",
"mcs = rdFMCS.FindMCS(ms,params)\n",
"print('scaffold:',mcs.smartsString)\n",
"scaff= Chem.MolFromSmarts(mcs.smartsString)\n",
"scaff"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [],
"source": [
"import os.path\n",
"from rdkit import RDConfig\n",
"from rdkit.Chem import ChemicalFeatures\n",
"fdefName = os.path.join(RDConfig.RDDataDir,'BaseFeatures.fdef')\n",
"factory = ChemicalFeatures.BuildFeatureFactory(fdefName)\n",
"factory."
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"11"
]
},
"execution_count": 29,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"activeMol = ms[0]\n",
"\n",
"feats = factory.GetFeaturesForMol(activeMol)\n",
"# for this exercise strip that down to a limited set of types\n",
"includeOnly=('Donor','Acceptor','NegIonizable','PosIonizable','Aromatic')\n",
"activeFeats = [x for x in feats if x.GetFamily() in includeOnly]\n",
"len(activeFeats)"
]
},
{
"cell_type": "code",
"execution_count": 33,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0 Donor SingleAtomDonor (2,)\n",
"1 Donor SingleAtomDonor (5,)\n",
"2 Acceptor SingleAtomAcceptor (0,)\n",
"3 Acceptor SingleAtomAcceptor (2,)\n",
"4 Acceptor SingleAtomAcceptor (12,)\n",
"5 Acceptor SingleAtomAcceptor (24,)\n",
"6 NegIonizable AcidicGroup (1, 0, 2)\n",
"7 PosIonizable BasicGroup (5,)\n",
"8 Aromatic Arom5 (11, 12, 13, 14, 24)\n",
"9 Aromatic Arom6 (7, 8, 9, 10, 25, 26)\n",
"10 Aromatic Arom6 (15, 16, 17, 18, 22, 23)\n"
]
}
],
"source": [
"for i,feat in enumerate(activeFeats):\n",
" print(i,feat.GetFamily(),feat.GetType(),feat.GetAtomIds())"
]
},
{
"cell_type": "code",
"execution_count": 38,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"1 Donor SingleAtomDonor (5,)\n",
"4 Acceptor SingleAtomAcceptor (12,)\n",
"5 Acceptor SingleAtomAcceptor (24,)\n",
"7 PosIonizable BasicGroup (5,)\n",
"8 Aromatic Arom5 (11, 12, 13, 14, 24)\n",
"9 Aromatic Arom6 (7, 8, 9, 10, 25, 26)\n",
"10 Aromatic Arom6 (15, 16, 17, 18, 22, 23)\n"
]
}
],
"source": [
"# show only those features for which all atoms are contained in the MCS match:\n",
"match = set(activeMol.GetSubstructMatch(scaff))\n",
"for i,feat in enumerate(activeFeats):\n",
" aids = set(feat.GetAtomIds())\n",
" if aids.intersection(match) != aids: continue\n",
" print(i,feat.GetFamily(),feat.GetType(),feat.GetAtomIds())"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment