Skip to content

Instantly share code, notes, and snippets.

@iwatobipen
Created December 5, 2021 11:43
Show Gist options
  • Save iwatobipen/818fcb5059fb389a60a97c278ac829e6 to your computer and use it in GitHub Desktop.
Save iwatobipen/818fcb5059fb389a60a97c278ac829e6 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "italian-mustang",
"metadata": {},
"outputs": [],
"source": [
"import os\n",
"import sys\n",
"from rdkit import Chem\n",
"from rdkit.Chem.Draw import IPythonConsole\n",
"from rdkit.Chem import RDConfig\n",
"from rdkit.Chem import Draw\n",
"from rdkit.Chem import AllChem\n",
"from rdkit.Chem import DataStructs\n",
"aaps = os.path.join('/home/iwatobipen/miniconda3/envs/chemoinfo/share/RDKit/Contrib', 'AtomAtomSimilarity')\n",
"import numpy as np\n",
"np.seterr(all='ignore')\n",
"Draw.rdDepictor.SetPreferCoordGen(True)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "parliamentary-radius",
"metadata": {},
"outputs": [],
"source": [
"sys.path.append(aaps)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "downtown-reggae",
"metadata": {},
"outputs": [],
"source": [
"import AtomAtomPathSimilarity\n",
"def get_tc(m1, m2):\n",
" fp1 = AllChem.GetMorganFingerprintAsBitVect(m1,2)\n",
" fp2 = AllChem.GetMorganFingerprintAsBitVect(m2,2)\n",
" return DataStructs.TanimotoSimilarity(fp1,fp2)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "funny-advancement",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAyAAAADICAYAAAAQj4UaAABsDklEQVR4nO3dd1hUV/oH8C9dpIuAoKKIhSIqipjYUMS1RDQxwqaBySaBTTYB0xaSGMFNWcwvm4CpuGlgmphNNmiiCWjsBUSigICANEGRIojSZ97fH3fn6tBEHeYyw/t5Hh7lzpk779Vz5p733nPO1SEiAmOMMcYYY4ypga7UATDGGGOMMcYGD05AGGOMMcYYY2rDCQhjjDHGGGNMbTgBYYwxxhhjjKkNJyCMMcYYY4wxteEEhDHGGGOMMaY2+lIHwBhjUisvL0dSUhKGDRuGGTNmYMqUKVKHxJjKKOq3lZUVvLy8uH4zxiTHCQhjbNC6cuUKXnvtNXz88cdob2+Hvr4+ZDIZHnnkEbz99tsYMWKE1CEydtsaGhqwYcMGrt+MsQGHh2Bpo+Zm4LffgE8/BfbvB9rbr792/Dhw7pxy+epqYM8e9cbImISICN999x0mT56MzZs3o6OjA6tXr8batWuhr6+PrVu3wt3dHR988AE6OjqkDpexW0JE+Oabb+Dh4aFUv4ODg5Xq94cffsj1mzEVKG8rx6Grh6QOQ6NwAqJtDh8GnJyAt94Cjh4FXngBcHcH8vOF1//xD2DHDuX3nDwJhIaqP1bGJHDw4EF4eXnhwQcfRHl5OebOnYv09HT85z//waeffoqioiIEBQWhrq4Ozz77LFxdXbF9+3apw2asTw4cOAAvLy88/PDDKC8vx7x583DixAn85z//wWeffYbc3FwEBASgrq4OzzzzDNzc3LCj8zmBMdat4rZi/LXsr/DO88Zd+XdhU9UmtFM79l/dj40XNkodnkbhBESbtLcDDz8MvPQSsG8f8NlnQHo64OcHPPGE1NExJqmSkhL4+/tj/vz5OHnyJMaOHYvk5GQcPHgQM2bMEMuNHj0aiYmJ2Lt3L6ZMmYLCwkIEBgbCz88POTk5Eh4BYz0rLi6Gv78/fHx8lOr3gQMHMH36dLGcs7MzkpKSsGfPHnh4eKCgoAArV67E4sWLuX4z1ouq9irMyZ+Delk9No3chDcc3sCexj34qeEnqUPTSJyAaJMTJ4D6euDZZ69v09EB1q8HDh0CqqokC40xqTQ2NiIyMhJubm7YuXMnTE1NERMTg5ycHPj7+/f4voULFyIjIwPx8fGwsbHBnj174OnpifDwcDQ0NKjxCBjrmaJ+u7u7Y+fOnTAzM0NMTAzOnDnTa/329fXFyZMnxfqdmprK9ZuxXmyu3gwHAwd86/QtFpothJ+ZH34b/xvWWK6ROjSNxAmINiksFIZfGRoqb3dwAMzNhdcB4MMPAV/f6z8REeqPlbF+RkRITEyEq6srNm3ahJaWFgQFBSEvLw8REREYOnToTfehr6+PkJAQ5OfnIywsDESEzZs3w9nZGXFxcZDJZGo4Esa6ksvl2LJlC1xcXJTqd25uLiIiImBsbHzTfSjqd15eHtdvxm7ij+Y/cK/lvdCBjtShaAVOQLTJ0KFAS0vX7XI50NoKmJgIv69cCbz99vUfHp7FtMyRI0cwa9YsrF27FhUVFbj77rtx7NgxJCYmYuTIkb2+t7tOl5WVFeLi4pCVlYWlS5eitrYW69atg7e3Nw4d4omHTL0U9Ts0NBSVlZWYPXs2jh8/ftv1e9iwYVy/GbuJqvYq2OjbSB2G1uAERJuMHw+UlgKdb5/n5gIdHcC4ccLvo0cDXl7XfyZMUH+sjPWDsrIyBAYGihPLR4wYgYSEBBw+fBje3t592sc///lPeHt74+jRo11ec3Fxwa5du5CcnIxx48bh5MmTmD9/PgIDA1FaWqrqw2FMyY31+8SJE3B0dERSUhIOHTqEmTNn9mkfb7zxBmbNmoVjx451ea23+l1WVqbqw2FMozgYOOBi+0Wpw9AanIBok6lThZ8XXgDa2oRtV64Ivz/4oDAMizEt1NLSgujoaLi5uWH79u0wMjJCVFQUCgsLERwcDB2dvt0y7+jowOeff4709HTMnz8ff/vb31BbW9ulnL+/P3JzcxEbGwszMzNs374dbm5uiI6ORkt3dyEZuwOK+q1YkU1Rv8+cOYOAgIA+1+/29nZ88cUXSEtLw7x58/Dss8+irq6uS7kb67epqSm2b98OV1dXrt9sUJsxdAa+r/8eBJI6FO1ATLtUVhItXEg0bBjR1KlElpZEgYFEjY3C68uXE8XGKr9n924iZ2e1h8qYKiQkJNDo0aMJAAGgoKAgKisru+39Xb16laKiomjIkCEEgExMTCgqKoqam5u7LV9RUUEhISGkq6tLAGj06NGUkJBw25/PmIJcLqeEhAQaNWqUyuu3kZGRUv1uaWnptnzn+u3s7ExJSUm3/fmMaara9loanz2elhYspeT6ZDp89TBFV0bTqaZTtLV2K/md9ZM6RI3CCYi2ys8n2rGD6Px55e0NDUSdO1JtbUR1deqLjTEV+OOPP2jBggVix2zq1Km0d+9ele2/oKCAAgICxP2PHz++145Xeno6zZ49Wyy/cOFCOn36tMriYYNLZmYm+fj4KNXv33//XWX7v9X6nZaWRnfffbdY3tfXl+s3G3Rq22spsiKSFhcspnuL7qWNlRupoq2CDjYepLcuvCV1eBqFExBt9cwzRADR++9LHQljKnXx4kUKCgoSr8ja2tpSQkICyWSyfvm8PXv2kIeHR586Xoor1iNGjCAApK+vTyEhIVRdXd0vsTHt07l+29nZ9Wv9Tk1NpcmTJ4v1e9GiRZSVldVtWa7fjBEl1CZQSGkIHWo8JHUoGo0TEG31t78JCcgHH0gdCWMq0dLSQjExMWRpaUkAyMjIiCIiIujy5cv9/tnt7e0UHx9PNjY2fep4dR7GNWzYMIqNjaWOjo5+j5VpppaWFoqKiiJTU1PJ6vfw4cNvqX4rhnFx/WaDyV9K/kLIAH1W85nUoWg0TkC0FScgTIskJSXRhAkTxKu0K1asoPz8fLXHUVtbS2FhYaSvr6/U8Wpvb++2fOdhLq6urvTrr7+qOWo20CUlJdH48ePFehIQEEBnz55VexyK+q2np9enxKJz/Z42bRrt379fzVEzpl6PlTxGyAB9XvO51KFoNE5AtBUnIEwLZGVl0aJFi8QOzuTJkyklJUXqsCg3N5eWLl0qxuXi4kK7du3qsXznYVwrVqyg4uJi9QXMBqTTp0+Tr6+vUv1OTU2VOiw6c+YMLVmyRClx3r17d4/lOw/j4vrNtNmjJY8SMkBf1HwhdSgajRMQbfX000IC8uGHUkfC2G3ZvHmzOMTD1NSU3nrrrR5XopJKcnIyjRs3TqnjVVRU1G1ZxRAyMzMzAkBmZma0efNmNUfMBop3332XDA0NxboQExPT40pUUklOTiYnJyel+n3u3Lluy7a0tNA///lPcQiZubk5vc9zEJkW4gRENfg5IIyxASc/Px9HjhxBe3s7wsLCUFpaipdffhlDhgyROjQlNz4vwdzcHDt37oSrqyvCw8PR2NioVNbIyAgREREoKSlBWFgYmpqa8Ntvv+HUqVMSRc+kkpubi+PHj0MmkyEsLAwlJSWIiIiAkZGR1KEp8ff3R15envi8m507d8Ld3R2RkZHd1u/IyEiUlpYiLCwMV69exaFDh5CXlydR9Iz1D8VzQHTQt+fvsO5xAqKt6H8PyunjA6oYG0iKiorw3XffYd68eYiLi8OwYcOkDqlHhoaGCA8PR25uLkJCQtDR0YHNmzfDxcUFW7ZsgVwuVyo/bNgwxMXFYeXKldi5cyfS0tIkipxJpaCgANu2bYOPj4/G1O+8vDwEBQWhpaUFmzZtgouLCxITE0GKc83/KOq3v78/tm3bhsOHD0sUOWNsIOMEhDE2YJmYmEgdQp85ODggPj4ee/bswdSpU1FZWYnQ0FD83//9X7flra2tAaDPT7Fm2sfU1FTqEPrMwcEBiYmJSE1NhYeHByorK7F27Vr861//6rY812+mrRRJN9ftO8MJiLbiOyBMg2nyF/yCBQuQmZmJpKQkeHh44C9/+YvUIbEBRpPrt6+vL06dOoWEhAS4u7tz/WaM3RZOQBhjTMV0dHQQEBCAzMxM2NjYdFtGkzuhbHDT0dFBcHAwTp061ePwMa7fTFvxHBDV4ARESx2xtcUuHx9kadAQFsYUeuq87N+/H6Ghofjmm2+kCOuW6enpSR0CG4B6qt+///47QkND8d1330kR1i3rrX53nhvCmLbQadWBhY4FqJXr+J3gBERLJVZVYfn+/TjS1CR1KIypzJkzZ7BlyxYcOnRI6lDuGF8hZp1lZ2djy5YtOHLkiNShqAzXb6Zt5G/I0TC9Abp7uAt9J/hfT0tx54ZpMq6/TJsNhvo9GI6RDU5ct1WDExDGmMbQpi9+bToWphpcJxhjgwUnIFqKT2RMk3H9ZdpsMNTvwXCMjLHbxwmIluIvf6aNtKlea9OxMNXgOsHYwMftVDU4AWGMDTj8Bc+02WCo34PhGBljt48TEC3FX/5MG2lTvdamY2GqwXWCsYGP26lqcALCGBtw+AueabPBUL8HwzEyxm4fJyBair/8mSbrqf5qU73WpmNht2Yw1G/GtBW3U9XgBEQLtbS0ICsrCwCQlpbGT6RljDGmVtxJY9qqrq4OAHD58mWJI9FsnIBoEblcji1btsDZ2Rnp6emwsLBAfHw8Zs6cicOHD0sdHmN9NhiuEGvTsbBbw/WbMc2Tk5ODxYsX47fffoO9vT2ef/55REZG4urVq1KHppE4AdESf/zxB3x9fREaGorKykpMnToVTzzxBGxtbZGRkYH58+cjJCQE1dXVUofKGGOMMaYRrl69ivXr12PmzJlITU3F0KFDYWNjg5aWFmzatAnTpk3DTz/9JHWYGocTEA1XWlqKwMBATJ8+Hfv374ejoyOSkpKQmZmJd955B8XFxYiKioKhoSH+/e9/Y+zYsYiOjkZLS4vUoTPWI75CzLQZ12/GBj6ZTIa4uDiMGTMGb775Jtra2hAWFoaysjKcOnUKaWlpuOuuu1BUVIR7770Xs2bNwvHjx6UOW2NwAqKhmpubER0dDTc3N2zfvh1GRkaIiopCbm4uAgICxC/9oUOHIjo6GgUFBQgKCkJTUxM2btyIiRMnIjExUeKjYIwxxhgbWHbv3g0PDw+sW7cOdXV18PHxQXp6OuLi4mBtbQ0AmDlzJo4cOYKEhATY2dkhLS0Ns2fPRnBwMKqqqiQ+goGPExANQ0RITEzExIkTsXHjRjQ3NyMoKAgFBQWIjo7G0KFDu33fqFGjkJiYiN9//x1TpkxBeXk51q5dC19fX5w+fVrNR8FY7/gKMdNmXL8ZG5iKiorg7++PZcuWITc3F+PGjUNycjL27dsHT0/PLuV1dHQQHByMwsJCREVFwcDAAFu3boWLiws2bdqEtrY2CY5CM3ACokGOHTuGu+66C2vXrsX58+cxbdo0/P7770hMTMSoUaO6fc+ZM2ewc+dO8fcFCxYgMzMTCQkJsLW1xe+//w5PT08EBwfz/BDGGGOMDToNDQ0IDw+Hq6srdu7cCXNzc8TGxiI3Nxf+/v7dvqempga1tbUAAFNTU0RHRyMrKwsrVqxAfX09IiMj4eHhgZ9//lmdh6IxOAHRABcvXkRwcDBmz56NtLQ02NnZISEhARkZGfDx8enxfUSEkJAQ+Pv7Y8GCBTh16hQAQFdXF8HBwcjLy0NERAT09fWxdetWTJo0CXFxcejo6FDXoTHWLb56yrTZYKjfg+EYmeZTzPNwdnbG5s2bIZfLERISgry8PISHh8PQ0LDH965btw5jxoxBdHQ0WltbAQATJkzAjh07kJKSAnd3d5w9exYrVqzA4sWLcebMGXUdlkbgBGQAa2lpQXR0NMaPH4+tW7fC0NAQUVFRKCwsRHBwMHR1e//vIyIEBATA0tIS+/fvx8yZM/HCCy+goaEBAGBlZYWYmBicPn0ay5cvx+XLl7Fu3TpMnjwZu3btUschMnZLtKlTo03HwlSD6wRj6nPgwAHMnDkT69atQ21tLebNm4e0tDTEx8fD3t6+1/e2t7ejpqYG165dw8aNGzFz5kzs3btXfN3Pzw+ZmZmIjY2FhYUFUlNTMW3aNISHh+PKlSv9fWgagROQAWr79u3w8PDAxo0bce3aNQQEBCA7OxvR0dEwNTXt0z50dXURHh6O4uJiREREQEdHB++++y6cnJyU7nRMmjQJP//8M5KTk+Hs7Iz8/HwsX74c/v7+OHfuXH8eJmPd4o4Y02aDoX4PhmNkmuncuXPw9/eHj48PMjMz4eTkhOTkZBw4cADTp0/v0z4MDAywe/du7NmzB5MnT0ZWVhYWLVqExYsXIycnRywTHh6OoqIihIWFQS6XY/PmzXBxccGWLVsgl8v78zAHPE5ABpjTp0/D19cXgYGBKCwshIeHB/bs2YOkpCSMHz/+tvZpaWmJmJgYZGVl9Xqnw9/fH7m5uYiNjYW5uTl27twJV1dXztjZgKFNnRptOhamGtpUJxTHwthAceXKFaV5HmZmZoiJiUFOTk6P8zxuxtfXV7zTYWlpidTUVHh6eiI8PFwcbWJtbY24uDikpaVh7ty5uHDhAkJDQ+Ht7T2oHxLNCcgAcenSJQQHB8PT0xO///47bGxsEB8fj5MnT8LX11clnzFx4sSb3ulQZOx5eXkICQlBR0cHZ+xM7bSpI8ZYZ4Opfg+GY2QDm0wmw5YtW+Di4oLNmzejo6NDnOcREREBY2PjO9q/vr5+t3c6nJ2dERcXB5lMBgCYPn06Dhw4gKSkJIwZMwYZGRmYN28eAgMDUVZWpopD1SicgEisra0NmzZtwqRJk7B161bo6+sjIiIC+fn5CAkJgb6+vso/s7c7HY2NjQAAe3t7xMfHIy0tDXPmzBEz9lmzZuHIkSMqj4mxvtCmjps2HQtTDW2qE9p0LExzHTx4ELNmzUJoaCguXLiAOXPm4Pjx44iPj4eDg4NKP2vYsGGIi4tDeno65s2bh9raWqxbtw6zZs0S73To6OggICAAZ86cQVRUFIyMjLB9+3a4ubkNvodEE5PUypUrCQABoPnz59PJkyfV+vmVlZUUEhJCurq6BIAcHBwoISGB5HK5WKajo4M++eQTsrGxIQCkp6dHn332mVrjZIPLt99+SwDogQceUNpeUlJCKSkpdPbsWYkiU50HHniAANC3334rdShMzbZu3UoA6JFHHlHaXlxcTCkpKVRQUCBRZKqzZs0aAkDbt2+XOhQ2SK1bt07s29jZ2dGWLVtIJpOp7fOTk5NpzJgxBIB0dHQoICCASktLlcrk5ubS0qVLxX6gq6srHT9+XG0xSonvgEjo+++/R2lpKRwcHJCSkoL9+/d3+6Cb/qS403H8+HHMnj0blZWVWLt2LWbNmoVjx44BAPT09BAaGori4mJERUXB0NAQP/74IyoqKtQaK2NjxoyBn58fJkyYIHUofdLb1SziK8Ssk7Fjx8LPz++25/upG9dvNlClp6ejuLgYRkZGiImJwblz5/Dkk0/edPVQVfL397/pnQ4XFxfs2rULe/fuxZQpU1BQUIDffvtNXNZXm3ECIqHS0lKcOnUKgYGB8PPzkzQWLy8vHDhwAB999BGsra3FW4j/+Mc/xDImJiaIjo7GqFGjsHPnTnG4FmOqpumdF8VDrRwcHHDhwgWpw2EDjKbX75qaGoSGhmLkyJG4ePGi1OEw1sWhQ4fw008/4cEHH0RERASGDh0qSRxDhw5FdHQ0MjMzsXTpUnHZXi8vL6VlexcuXIj09HQYGxvjtddew7Vr1ySJV504ARkA9PT0pA4BgBDHU089hdLSUkRFRUFPTw+TJk3qUk5xBUFTT56M9RciwtatW+Hu7o7Nmzejvr4ev/zyS49lAW5HTHPI5XIkJCRgypQp2LJlCy5fvoxff/2127Jcv9lAYGFhIXUIAJTvdHh4eCAnJ6fLgwkNDQ3Feb+Dod1wAiKhgfoFrbjTcfr0aQQGBkodDhuEFG2jqKhIfF7NQLd79264u7sjODgYFRUVmD9/PjIyMvD44493KSuXy1FeXi5BlGwgUNTvgoICjanfu3btgru7Ox599FFcuHABS5YsQU5ODtauXdulrFwux/nz5yWIkjHBQO1fKe50fPTRRwgNDZU6HElxAiKh7hrIs88+Cy8vLxw9elSqsEQTJ07stvEO1IbNtMfSpUvx2GOPIS0tDa6urti+fbvUIfWoqKgI/v7+WLZsGXJzczFu3DgkJyf3OKdLkagcPXoUjz/+OHx8fCSImknpnnvuQXBwMI4dOwY3NzeNqN/Lly9HXl6eWL93794NV1fXLuV/+eUXuLu7Iy0tDSEhIZg7d64EUbPBjrp5Ds3nn3+OyMhInD17VoKIrjMyMsJTTz0FAwODLq8Npv4VJyADTGFhITIyMsQH2AxEg6mBMGkMGzYM999/P8aPH4/CwkIEBgbi3nvvRVFRkdShiTo/1Mrc3ByxsbHIzc3t9qFW5eXlePjhh8WO3KhRo7BkyRLY2dlJED2TkrW1Ne6//344OzujoKAAgYGBWL16tfg8poFAMY+pr/W7rKwMDzzwAFasWIG8vDyMHj0aS5Ysga2trQTRMya4sZ+SlJSETZs2obi4WMKImAInIBLqriPPnXvGBPfccw9yc3MRHx8PGxsb/PTTT3BxcUFoaChqamoki0smkyEuLg7Ozs7YvHkzZDKZ+FCr8PBwGBoaKpVXJCrOzs745ptvYGZmhtjYWBQVFSEgIECio2BSW7lyJfLy8hAfH4/hw4fjxx9/HJD1Wy6X91q/FYnK+PHjsW3bNqX6vXr1aomOgg12mtq/0oQYVYUTkAFGEyqfJsTItIO+vr7Y+QkLCwMAbNmyBZMmTVJ6wqy6HDhwAN7e3li3bh1qamowd+5cpKWlIT4+Hvb29kplb7UjxwYfRf3Oz88Xn6AsZf3ev38/Zs6ciXXr1qG2thbz5s27pfodFhaGoqIihIeHdzu8hDHGFDgBkRB35BnrG8UTZrOysrB06VLU1dVh3bp18PDw6HEVHlUqLi6Gv78/fHx8cPLkSYwdOxbJyck4ePAgZsyY0aX8wYMHuyQqiqfvdu7IMXZj/V6yZIlYvxXLo/e3c+fOwd/fHwsWLEBmZiacnJyQnJyMAwcOYPr06V3KHzhwoNtEJS4uDsOHD+/3eBm7Gb4DMvBxAjLAaELl04QYmXZSLGWYnJyMcePGITc3F0uXLoW/v3+/jOttbGxEZGQk3N3dsXPnTpiZmSEmJgZnzpzpdhx8SUkJ/P39MX/+fDFRSUpKwoEDB7pNVBi7kaurK3bv3o3k5GQ4OTnhjz/+gI+PD/z9/VFSUqLyz+s8j0lRv3Nycrqt3zcm4pmZmWIi3lOiwthAogl9F02IUVU4AZHQYKpojKmSv78/cnNzERsbCzMzM+zcuRPu7u6IjIxUyQMyFUNhXFxcsGnTJrS0tCAoKAi5ubmIiIiAsbGxUnlFouLm5oadO3fC1NRU7MgFBARwG2e3xN/fHzk5OYiJiRHrt5ubGyIjI3H16tU73r9MJhPr9+bNm9HR0SEOD+ytfvc1EWdMaty/Gvg4AZEQ3yJk7PYZGhoiPDwceXl5CAoKQktLCzZt2gQXFxckJiZ2uwxjXxw6dAizZs1CaGgoKisrMXv2bBw/fhyJiYkYOXKkUlkiQmJiIlxdXZUSFUVHTqqn7zLNZ2xsjIiICOTm5qq0fh88eFCs3xcuXMCcOXPE4YEODg5KZW81EWdsINOEvosmxKgqnIAwxjSag4MDEhMTcfz4cdx9992orKzE2rVrMWvWrFt+nk5DQwOWL1+OEydOwMbGBp988ok43r2zI0eOYNasWVi7di0qKipw991349ixY90mKozdrpEjRyIxMRHHjh3DXXfdhYqKCqxduxZ33XUXjh07dkv7unz5MpYtW4aMjAzY2dlhy5YtOHDgALy8vLqUPXz4cJ8TccYGmsHUkddUnIBIiO+AMKY6M2fOxOHDh5GQkIARI0YgPT0dc+bMQXBwMC5evNinfVhYWOCVV17BCy+8gIKCAoSGhkJPT0+pTFlZGQIDAzF37lykp6djxIgRSEhIwOHDh+Ht7d0fh8YYvL29ceTIESQkJMDOzg5paWli/a6qqurTPqysrBAZGYmXXnoJZ8+exZNPPgldXeVugKJ+z5s3DydOnICjoyOSkpJw6NChbhNxxgYi7l8NfJyAMMa0ho6ODoKDg1FYWIioqCgYGhpi69atGD9+PKKjo9Ha2nrTfURGRuKdd96BhYWF0vaWlhZER0eLT642MjJCVFQUCgsLERwcPChOGExaivpdVFSEqKgoGBgY3HL9Xr9+Pd5++22Ym5srbW9ubkZ0dDRcXV2V6veZM2d4HhNjTOU4AZEQZ+iM9Q8TExNER0cjOzsbAQEBuHbtGjZu3IjJkydj+/btt7y/xMRETJw4ERs3bsS1a9cQFBSEs2fPIjo6GiYmJv1wBIz1TFG/s7KyEBAQgKtXr2Ljxo3w8PDAzp07b2lfinlMivrd3NzM9ZtpPO5fDXycgDDGtNb48eORlJSE1NRUTJ48GYWFhQgMDISfnx+ys7Nv+v5Tp05h4cKFWLt2LcrLyzF16lTs3bsXiYmJGD16tBqOgLGeTZgwAUlJSUhJSYG7uzsKCgrg7++PxYsXIycn56bvV8ybWrt2Lc6fP49p06Zx/WaMqQUnIBLiDJ0x9Vi0aBEyMzMRHx+P4cOHY8+ePfD09ERoaChqamq6lK+qqkJwcDCmT5+Offv2wdbWFgkJCTh58iQWLlwowREw1jM/Pz9kZmYiNjYWlpaWSE1NhaenJ8LDw9HQ0NClvKJ+KyaW29nZISEhARkZGViwYIH6D4AxFeP+1cDHCYiENLWBMKaJ9PX1ERISgvz8fISFhYGIsGXLFkyaNAlxcXGQyWRobW0VlzrdunUrDAwMEBERgfz8fAQHB3eZsMvYQGFgYIDw8HAUFRUhLCwMcrkcmzdvhrOzs1L9jo6OhrOzs1i/b5zHxPWbMaYu/G3DbhknSUyTDRs2DHFxceIyunV1dVi3bh08PT0xZcoUREZGor6+HgsXLsTx48cRExMDS0tLqcNmrE8U9fvgwYOYOXMmamtrsW7dOkyfPh1TpkwR5zEtWrQI6enpiI6OhqmpqdRhM6ZSmnqBVxNiVBVOQCSkqQ2EMW3g7e2NY8eOITk5GWPHjkVRURHOnTuHyZMnIyUlBXv37sXUqVOlDpOx23L33XcjLS0NycnJGDNmDAoLC3Hu3Dl4eHggNTUVqamp8PDwkDpMxvoF96UGPn2pA2Cahxs20yb+/v7w9vbGiBEjYGxsjPT0dAwZMkTqsBhTCX9/f3h6emL06NEwNTVFeno6jIyMpA6LMbXThL6LJsSoKpyASKi7irZgwQLY2trCxsZGqrBuajA1EDY4KJYa1dPT4+SDaR3FECt9fX1OPtig0F0/xcHBAePGjePv+AGCE5ABJjo6WuoQRJWVlXBwcJA6DMb6HSfVTJtx/WYM+Pbbb6UOQSSTyaCnp9dl+2BqqzwHREIDtaJVV1cjNDQUY8aM6fahVgM1bsYYY4yxgdpPaWtrw6ZNmzBq1CgUFBR0eX2gxt0fOAEZANra2qQOAQAgl8vx6aefwsPDA1u2bIFcLsfJkye7lFM0EMa0xWD60meDD9dvxqR34MAB3HXXXYiMjMTFixeRmJgodUiS4gREQjNmzMBdd92F999/X1ynXSq//PIL3Nzc8OSTT6KqqgrLly9HTk4ONmzYIJZR3Bmpq6vDypUrYWxsLFm8jDHGGGPdsbe3x7hx45CYmIgjR45IGktOTg4WL14MHx8fZGZmws3NDSkpKXj99dfFMnK5HJ9//jmMjIwwbty4QXGxgBMQCfn7+8PS0hKNjY1Yt24d5s2bh/T0dLXGcObMGSxevBj33HMP8vPzxYbx888/w8XFBYDQML744gtMmTIFW7ZsQVtbGwIDAzFs2DC1xspYf+ErxEybcf1mg82TTz4JW1tblJSUYP78+XjqqadQU1Oj1hhqamoQGhoKT09PpKamwsLCArGxsfjjjz/g5+cnljt48CC8vb3x+OOP4+rVq4iNjYW+/iCYok1McsnJyeTk5EQACACtWLGCiouL+/Uz6+vrKSwsjAwMDAgAWVhYUGxsLLW1tSmV+/nnn8nFxUWMbenSpZSbm9uvsTGmbvX19WI7YEzb1NTUEACytraWOhTG1Obq1asUFRVFQ4YMIQBkYmJCUVFR1Nzc3K+f297eTrGxsTRs2DACQHp6ehQWFkY1NTVK5bKzs8nPz0/sXzk7O1NycnK/xjaQcAIyQDQ1NVFMTAyZmpoSADI2NqaIiAhqbGxU6ed0dHRQbGwsWVtb99owcnJyBnXDYIPL5cuXCQBZWlpKHQpjKlddXU0AaPjw4VKHwpjalZWVUVBQkNifGT16NCUkJPTLZ3V30fbMmTNKZS5fvtynC8DajhOQAeb8+fMUFBREOjo6BIBGjhxJCQkJJJfL73jf+/btI09PT7FhzJ8/n06ePKlUpvOdEXNz80HZMNjgwgkI02acgDBGtHfvXpoyZYrYB1q4cCGdOnVKJfvuy90MxQXgm90ZGSw4ARmgjh8/TnfddZdYmWfNmkXHjh27rX0VFRXRihUrxH05OTn12DA63xmprq5WxeEwNqDV1dURALKyspI6FMZU7tKlSwSAbGxspA6FMUnJZDJKSEggGxsbAkC6uroUFBREly5duq399fVuxq5du8jV1VXsh/n4+HS5ADzYcAIygCkaip2dnVJDuXjxYp/e39raSm+88QaZmJgQABoyZAi9/PLLXYZ17d+/X+nOyLx58ygjI6M/DomxAYkTEKbNOAFhTFldXR1FRESQoaGh+N0fExNDra2tfXp/X+d5FBYWKl0AHjduHA9n/x9OQDRAY2MjRUVFkZGRkThMJCYmhlpaWnp9n0wmI29vb9LV1aWQkBCqqKhQev3cuXNKDWPs2LHcMNigVFtbSwBo2LBhUofCmMpVVVURALK1tZU6FMYGlLy8PFq+fLnYD5o0aRL9/PPPN33f1atXaeTIkeIIlaNHjyq93tNw9r4mOIOBDhE/VU5TFBQU4PnnnxefTj5x4kS8++67uOeee3p8T1paGpqamrBgwQJxW2NjI958801s3rwZzc3NMDMzw6uvvoqwsDB+tgcblOrq6mBtbY1hw4ahtrZW6nAYU6lLly7Bzs4Otra2qKqqkjocxgac1NRUhIWFITc3FwDg5+eHzZs3w9XVtcf3/Pjjj7h27RoefvhhcYlrmUyGDz74AK+//jpqa2uhp6eHxx9/HNHR0bC3t1fLsWgKTkA0UGpqKtatW4ecnBwAQkOJi4uDm5tbr+9TPOl848aNqKyshK6uLp544glERUXBwcFBHaEzNiDV1tZi+PDhsLa2Vvta8Yz1t6qqKowYMQJ2dna4ePGi1OEwNiC1t7fjo48+woYNG3DlyhUYGBjgqaeewuuvvw5zc/Obvn///v147rnnkJmZCQCYN28eYmNjMX369P4OXSPxgwg1kJ+fHzIzMxEbGwsLCwukpqZi2rRpCA8PR0NDQ7fvOXz4MGbNmoXQ0FBUVlZi9uzZOHbsGOLj4zn5YIwxxtigZmBggPDwcBQVFSEsLAxyuRybN2+Gs7Mz4uLiIJPJun3fuXPn4O/vjwULFiAzMxNOTk5ITk7GgQMHOPnoBScgGqqvDaW0tBSBgYGYN28eTpw4AUdHRyQlJeHQoUOYOXOmxEfB2MBA/KRopsW4fjPWd8OHD0dcXBzS0tIwd+5c1NTUYN26dfD29sahQ4fEcleuXEF4eDhcXV2xc+dOmJmZISYmBjk5OfD395fwCDQDD8HSEkeOHEF4eDhOnDgBAJg9eza8vb3xxRdfoKGhAQYGBnjmmWewYcMGWFpaShssYwNMTU0NbGxsMHz4cFRXV0sdDmMqdfHiRdjb22PEiBG4cOGC1OEwplF27NiBZ599FqWlpQCAe+65B3PmzMH777+PCxcu8HD228QJiJbZsWMHwsLCUFJSAktLSzQ0NOCRRx7BW2+9hVGjRkkdHmMDEicgTJtxAsLYnbly5QreeOMNxMXFoa2tDUZGRmhtbYWnpyfee+89+Pj4SB2ixuEhWFrG398fp06dwvjx41FfX4833ngDiYmJnHww1gseosK0Gddvxu6Mubk53n77bWRnZ8POzg6tra2IiIjAiRMnOPm4TZyAaCFzc3OMHTsWAODt7S1tMIwxxhhjWmDChAlwdnYGAKxcuRK6utyNvl38L8cYG/T4CjHTZly/GWMDDScgWopPOIwxxhhjqsX9K9XgBIQxNujxCYVpM67fjLGBhhMQLcUnHMYYY4wx1eL+lWpwAsKYFnj2WWDFCiAwsOcyTz8tlHnggd731dgILF4s/Ozf332ZlBTh9e++u/2YBxI+oTBtxvWbMTbQcAKipfiEM7i4uQE//wx4evZcxt1dKDNtWu/7am8HUlOFnyeeAJqbu5aprBReLy6+o7AZY6xftchbcL79PK7Jr0kdCtMS3L9SDU5AGGPdMjQECguB//s/qSPpf3xC0UyvvQYYGwO7dnV9LSoKmDCh5/fOnw8YGAB//KG8/eBBQEcHKChQaaiSGgz1O+1aGuJr4vGvqn9hR8MOtMhbAABl7WUYnTUaZ1rOSBwhY+xGnIBoqcFwwmH9a84cYMoU4J//1K7OGNMePj5ASwvwv2X5lbS3d3/37kZyOfDXvwp/Ms3UIGvAssJluDv/brxb9S52X9mNoJIguJ1xQ15LntThMS3E/SvV4ASEMdYtPT3h7kdLizB/RJvxCWVwWrkSyM4GtmyROpL+pc31+8myJ3H42mEcnHQQ+e75SJmQgpLJJVhmsQwEkjo8xlgPOAHRUtp8wmHq86c/AffdJ8z3+PJLqaNhTLVGjADWrwf+/negokLqaNitKm0rxfeXv8czNs9gtslscbulniU+HP0hXIe4Shgd01bcv1INfakDYP2DG8jgUl0t/Ll+PfDDD92Xyc4W/rx48db2vWkT8MsvQGQksGoVYGV1+3EOVPX19QCAtrY2tLe3w8DAQNqAWJ+Ulgp/vvpq12FYBw70bR/h4cIdkIgI4KuvVBvfQHH58mUAQGtrKzo6OqCvrx2n/iPXjoBAWG25WupQGGO3iO+AMKYF/pdvQkdHGDrV3Y8iF73V8e4TJgAvvQRUVQEvv6zauKXW2tqKTZs2Yfbs2TA2NkZtbS2cnZ2RmJgodWisD4YMEf6srgbKypR/Ghr6tg9jY2DzZuDrr4G9e/svVim0trYiOjoac+bMgYmJCS5dugRXV1ds375d6tBU4lL7JQDAaIPREkfCBhO+wKsanIBoKW4gg4utrfDn668Dx451/6NYzcrB4db3v349MGkS8O9/A+npqotbSnv37sXMmTMRGRmJhoYGuLq6YuzYsSgvL8fatWuxcuVKFPDs+wHNzk74c8sW4JtvlH9Wrer7flasAAICgL/9TZi8rg1SU1Ph5eWFjRs34tq1a3BxccGYMWNQWFiIwMBA3HvvvSgqKpI6zDtioCPcqWyiJokjYYOJon/F7gwnIIyxmzIyEhIYuVwYqnIjmUyamG5XdnY2/Pz8sGjRImRlZWHy5MlISUlBRkYGCgsLkZCQABsbG+zYsQNubm4IDQ1FtWKMG9Na77wDlJcDX3xxfdu6dcChQ5KFdFuysrKwaNEiLF68GNnZ2fDw8EBqaipOnDiBgoICxMfHw8bGBj/99BNcXFwQGhqKmpoaqcO+JceuHcOBqwcw1nAsAKC4lR9IxNSPL/DeGU5AtBTfAWGq5u8P3H8/8PvvwH//e337c88Bfn7X55gMVNXV1QgNDYWnpyf27NmD4cOHIz4+HidPnoSfnx8AQFdXF8HBwcjPz0dYWBgAYMuWLZg0aRLi4uLQ0dEh5SGw29TYCOzfD+Tn91zG0VF4rsg33wi/79sHxMUJzwt54AFhWNdApqjf06dPx969e2FjYyPW70WLFgEA9PX1ERISgry8vG7rt2yAX00413oO/kX+uDv/bjxa8ihmDp0JMz0zfFWnpZN32IDE/SvV4ASEsUHsyhXhp6/eeQcYOhT46Sfh99ZWYez8nj3AzJnAK68AV6/2T6y3q6OjA3FxcXBxccGW/623GhYWhvz8fISEhHQ74dzKygpxcXHIzs7GsmXLcPnyZaxbtw4eHh7YvXu3ug+B3YHDh4Vn2vz4I/DYY8DHH/dc9oUXhGffAMDcuUBsLGBmBmzbBowfL0xYv5X2og5tbW3YtGkTJk2a1G397m7C+bBhwxAXF4esrCwsXboUdXV1Yv3+9ddf1X0IN9Uib8GbF9/ElNwp2NmwE0N0h+ChYQ/BVM8Ur454FV/WfolXKl/BpY5LkJEMeS15eK3yNTTKGqUOnTHWE2Jaac6cOQSADh06JHUoTA0++ogIIHrrrZ7LfPCBUOaf/ySqqiJaupRoyBAiY2OisWOJjhwRytXWCuX8/Lrfz6ZNwuuKz6utJQoLI9LTE7YNG0YUG0vU0aH647xVP//8M7m4uBAAAkDLli2j3NzcW95PcnIyOTs7i/tZsWIFFRUV9UPE7Fbs20cUE0NUWNj1tf37hXZRUECk+K/KzydycxP+/tNPRAcOdH1fdjZRfDxRfb3we3W1cv0ePnxg1e9JkyaJ9XL58uW3Xb/HjRunVL/PnTvXDxHfGhnJKL46nhxOOxAyQLoZuhRSGkIVbRVK5d6reo/sT9sTMiD+zM2fSwUtBZTfkk/IAKVdS5PoKJi28fT0JAB08uRJqUPRaJyAaClOQAaXW01A/vIXokceIWpsFDpSTzxBNGYMkVx+8wSkvZ3Iw6Pr52VmEs2ffz058fTsvoOnDjk5OeTn5yd2qJydnSk5OfmO9tnW1kaxsbFkbm5OAMjQ0JDCwsLoypUrKoqa9bcTJ4jmzbu992ZkEM2de71+T59OdPCgauPrq871283NjVJSUu5on62trRQbG0tmZmYDon4fvnqYvHK9xIRidt5sSr+W3mN5OcmppLWESlpLqEXeIm5vk7dRXUed0jbG7gQnIKrBCYiWmj17NgGgw4cPSx0KU4OcHKL0dKLevg+zs4UymZlC0tHaev21EyeETtWFC0KCceIE0dmzPe+rokIoc+FC19eSk4U7KoqO2ooVRCUlt31ot6S+vp7CwsLIwMCAAJC5uTnFxsZSW1tbt+XPnz9PjzzyCCUlJfX5MyorKykkJIR0dXUJADk4OFB8fDzJZDJVHQbrB21tRIsXE/3nP3e2n+RkIVkfCPXbwsKi1/pdXl5ODz30EP3www99/oyKigoKCgoiHR0dsX4nJCSQXC5X1WH0qrS1lALOBZBOhg4hA+SY5UhJdUkkJ/V8PmM3M23aNAJAmZmZUoei0QZNAtLSQlRW1nV7RQXRtWtdt9fUEJ0/r7zt6tWu2wYqTkDYrdi1SxiO1UM/5pY1NQlDY0xNhU7a0KFEUVFEzc2q2X9nHR0dFBsbS9bW1gSA9PT0KCwsjKqrq7st39zcTG+88QaZmJgQAJo0adItJxDp6eliOwNAXl5e3N4GqKYmooAAouho1eyvoYHopZeIDA2F+m1uTvT228pJvSr1VL9ramq6Ld/U1EQbN26koUOHEgByd3e/5QQiLS2N7r77brF+z5w5k44oxmn2g2ZZM0VVRpFJpgkhAzQkcwhFVUbRVdnVfvtMpjrdfX3KZMJd9c6amoQ2dKO2tq7bBipOQFRj0CQgBw8SjRjRdfu0aUTbt3fdHh5OpK8vXC1W+OEHIi+v/opQtRQnjv48YTDtIJcL80Eee0z1+z57VrhCrLhavHhxPe3atUuln7F//37xljgAmjdvHmVkZHRbViaTUXx8PDk4OBAA0tXVpZCQEKqoqOi2/M3I5XJKSkoiR0dHAkA6OjoUFBREF7q7NcQkUVQkDJfasEH4e1FR952i21FeThQURKSjQzR6dDs5O7ur/G7Bvn37lOr3/Pnzexz6oajf9vb2SvW7srLytj5bLpdTQkICjRgxol/r9y+//EL3f3i/ONxqacFSOtN8RqWfwfpPfr4wR6qzRYuIPvmk6/aoqK79q5QUIheX/opQtaZOnUoA6I8//pA6FI3GCUgvCYiTE5GPz/WTFScgTBtt3Ejk6Nj9cCpV2bdPaGuTJz9JAGjhwoV0+vTpO9rnuXPnaMWKFWLHbOzYsb3O8zh06BB5eXmJ5WfPnk3p6T2PKb8VV69epaioKBoyZAgBIBMTE4qKiqLm/rrlw/rsww+Fux83/qh6AvmvvxKtXv1vpcng+fn5d7TPoqIipfrt5OTUa/0+ePAgzZgxQyw/Z84cOnHixB3FoKCo30ZGRkr1u6XlzuZVnDlzhpYtWyYkN7o6tCx9Ge2o36GSmJn63E4CYmsrzBlU9K84ARl8OAHpJQGJiCCaMoXos8+EbZqUgNx1110EgI4ePSp1KGyAamsjevZZolGjep/voSqtrR30zjvvkIWFhTjJ9e9///stT3K9cuUKRUREkLGxMQEgMzMziomJoaampm7Ll5aWUkBAgDim3dHRkZKSkvplTHthYSEFBASIncDx48ff0vwSprlkMhklJCSQra2tePchKCiILl26dEv7aWhooLCwMDI0NOxT/S4pKVGq32PGjOm3+l1QUKCS+n3hwgUKCgoS51HZ2tpSQkICz6PSULeTgISGCouZfPmlsE2TEpApU6YQADp16pTUoWi0QZWADB0qJBY3/tjY9JyAvPIKUWoqkbU10aVLmpOA1NbW0ujRowkAffTRR1KHwwag4mJhWMqyZcKqV+pUV1dHYWFhpK+vTwBo2LBhFBsbSx03uSzdefiUYjjI+R4mZjU3N1NUVJQ4Dn7IkCEUFRVF17qb9KVie/bsIQ8PD7Gj5uvre8d3fJhmuHz5MkVERIgJhJWV1W3V75sND2xqalKq38bGxmqr36mpqTR58mSxfi9atIiysrJu+j5Fm1TMvTIyMqKIiAi6fPlyv8fM+k9+PpGuLtE77yj/TJjQcwLyt78Jd8ZtbYVzkKYkIC0tLeKQW1UPJx5sBlUCYmFB9Omnyj+Ojr0nIEREDz4oLFs60BOQzl/uiuEgfn5+lJOTI3V4bACZPFmYk2FhQWRldf1n9271xZCbm0tLly4VOzGurq7066+/dlu2u+FTaWk9r+ufkJAgJuEAKCgoiMq6W4WiH7W3t1N8fDzZ2NgQANLX16eQkJAeJ8Yz7ZKXl3dH9bun4YGKeRmjRo1SSsTLy8v783C6UNTv4cOH96l+JyQkiB03ABQQEEBn1XHrlfU7RQISHa38M25c7wkIEdHDDxOFhAz8BETR7kaOHEkAaOjQoWRgYEBhYWFUr3hoELslgyoBudUhWIoE5MIF4eFqGzYICUhpKdHvv/dntLfu119/JTc3N6Urrs8884yYjBgbG9OGDRvUcnWMDXylpdcn5N74I0X1SE5OJicnJ6WHoBUXFxNR1+ElNxs+9ccff9CCBQvEfU2dOpV+l7ix9nTHp729XdK4mHp095C/263fmZmZ5OPjI+5r2rRptG/fPjUeTVe1tbUUFhZGenp63d7RzMzMpIULF4oxT5kyhfbs2SNpzEy1bmcIliIBuXBBuPi1aZOQgDQ0EJ0ZYOsPnDp1SqkOu7q60urVq8Xv9JEjR9JXX32ltqWqtQUnIP9LQNrbic6du76M4o0JCJHwELehQ4UE5L77rj+oLTtbLeH36NSpU+Tr6ys2DA8PD0pNTRVfr6mpUTo5WFtb92k4AGPq1NjYSC+//LJ4187MzIz8/f3F+SIGBgb03HPP9ThUo6qqasCPKe+8rKm3tzcdP35c6rCYGijqt2ISd3f1+/nnn++xfl+8eFGpftvZ2Q24+n38+HGaNWuW0rLUDzzwgNIzS9555x1q7a+1iplk7iQBISKKixOSEBcXYbuurrC63MWL/RZyn3SeP2hvb6+0yl3nOVGenp60f/9+aYPWIJyATCP64gui5cuJXnhBeEpuVVXXBEQmI5o9W0hA3niDyMRESEKGDCF69VXhGSHq1LnDZWNjQ/Hx8T1eVc3IyKB58+aJDWX69Ol0UKrH+DLWg/Pnz4sPQRs1apQ4vKSn4VMtLS0UExNDlpaW4sT2gT6m/MYr4uvXr6cGTVn8nt2xnup3T8OnOtdvTZgzceMdTTMzM42Imd2Z3hKQDz4gWrNGSChWrCCqru6agMhkRLNmCQnI+vVERkZC/2r4cKKPP1b9qnU303l+lWL+4NUeOnqpqank7u6udJfz3Llz6g1aAw2aBKS+XkhCOjt+XMiyFc9zeustov/+V3hoYee5rVVV128N1tQQhYUJjQ4QJqrHxvZ/Q2ltbaWYmBiysrJS6nDV1dX16f3Jyck0ZswYcexwQEAAlajrMb6M9VFwcDABoIceeqjHMjt27KCJEycqfenf6dKn6nL16lVxPLyqlkplmuORRx4hAPTII4/0WCY5OZkmTJigVL81Zc5EY2OjOEeFV4HTfleuEH31VdftKSlCnyk3V/j9/feJEhKEkSOdH2VTUCA8EJdI6Hspnq8DEE2cSLRDDasz9zS/qi/zB9va2ig2Nla8q2lsbEwRERG3vMrjYDJoEpC+WrGC6Fb64ydPCmtZ438PWps2jai/7sB17nAtX76c8vLylMrU19fTJ93d87zBtWvXKCoqSlzGdOjQoRQVFdXjMo+MqdvGjRsJAG3YsKHb19evXy+2g7Fjx9L27iZyDXDTp08nAD0+NJFpr6ioKAJA0T08mv3VV19Vev7H999/r+YI7xwvVco6e/RR4aJvX+3bJzwKQdG/6s9h70ePHiVvb2+l+VWd5w+2tbXRb7/91ut+Og97HzlypMofTqotOAH5H7lceO7H++/f3nsTEogcHIRGoqNDFBxMVFGhmtsh2dnZ5OfnJzYMNzc3SklJUSrT+c5IXzpk5eXl4nAAADRq1ChuKGxAiI6OJgAUFRXV7esZGRniZNe2tjb1BqcinIAMXjdLQE6cOKHx9ZsTEHajrVuJXnzx1t8nkwn9KxsboX9lYCCMPlHVwlOKZ9Io+kE9za9SrKyoo6PTp+HrJ0+eVBr2PnPmTH4wdCecgJAwfyMwUJgDcuLE7U98amoiiokhMjMjMjOTk6PjTIqIiKDGxsbb2l91dTWFhIQoTeLr7oS0f/9+8vT0FCv6vHnzKDMzs8+fs2/fPvHJngBowYIF/IRPJqmbJSBEpDFjyr/88kuKiIigM52WdlG02ZOdxyIwrbdhwwYCQBs3buyxjKbU75MnT1JKSgrVKMYx/4/iOTj8/JvBraNDmCcbGkpUWEjUx9HiXVy4QLR2rXCBV0+PaP78MNq6dettXzDt7pk03c3zOHPmDC1btkzsH02cOJEOHDjQ58/pPOw9KCiILly4cFsxaxtOQEiY7xEff/3nTi/YFBURhYcfFCvs+PHjKTk5uc/v7+jooNjYWBo2bBgBID09PQoLC+vyBd+XOyN91dNTfKuqqm5rf4zdCcUV4t4SEE3h7+9PALp8B3ACMnj1JQHRFIrO2S+//KK0nRMQRiT0ryIirv/08CicPjt2jOiRR/4Q+z133313r8+E6k7n50QFBARQQUGBUpnOK8/dycqKnYe9m5iYUFRUFLW0tNzyvrQJJyD9KC0tTWlZQh8fn5veWdi1axe5uroqvafz3Yz6+noKCwu76Z2R29HdU3xjYmJ46USmVjcboqJJVqxYQQBoR6dZlNOmTSMAt3S3kmmH1157jQDQP/7xD6lDuWOKhy12fiq04knpfXlCOmO3Kjk5mcaOHau0SIPi+To9yczMVHpO1JQpU2jv3r1KZbq7M6KqVdw6D3sfP378oF6kQRes38ycORNHjhxBQkICbG1tsX//fkyfPh3BwcGorq5WKltUVAR/f38sW7YMubm5GDduHJKTk7Fv3z5MmzYNACCTyRAXFwdnZ2ds3rwZcrkcYWFhKCoqQnh4OAwMDO44ZktLS8TExCArKwvLly/H5cuXERkZiSlTpuCXX3654/0z1hdEBADQ0dGROBLGGGMDjb+/P86cOYOYmBiYmppi586dcHNzQ2RkJK5evapUtqqqCsHBwZgxYwb27dsHW1tbJCQkIDMzEwsXLhTLbd++HR4eHti4cSOuXbuGgIAAZGVlISYmBpaWlncc86hRo5CYmIi9e/diypQpKCwsRGBgIBYtWoSsrKw73r+m4QSkn+nq6iI4OBj5+fmIiIiAvr4+tm7dikmTJiEuLg51dXUIDw+Hq6srdu7cCXNzc8TGxiI3Nxf+/v7ifvbv3w8vLy+sW7cOtbW1mD9/PtLT0xEXFwdra2uVxz1x4kT8/PPPSElJgZubG/Lz83HPPffA398fRUVFKv88xrQVJ1OsM22qEz0dizYdIxuYjI2NERERgby8PAQFBaGlpQWbNm2Ci4sLEhMT0dLSgujoaIwfPx5bt26FgYEBIiIikJ+fj+DgYOjqCl3gU6dOwdfXF4GBgSgsLISHhwdSU1ORlJSECRMmqDzuBQsWIDMzEwkJCbCxscHevXsxffp0hIaGdrk4rdUkvf8yCGVkZNCcOXPEW4Dm5ubi3wMDA7s8k6OoqEgcwoH/Lcl4K/NJVKGpqYmio6PFh/IMHTqU9u3bp9YY2OCiTUNU7rnnHgJAO3fuVNquWPiBF3wYfBTLSL/++utSh3LHlixZQgBo9+7dStsVD2bL7q91Uxnr5LffflN6IODw4cPFv/v6+nZZka27eR69PdC5P1y8eJEef/xxMYYRI0YMmmGLfAdEzaZPn45Dhw4hOTkZ1tbWMDMzg7u7O06cOIFt27ZhzJgxAIArV64o3RkxMzNDTEwMcnJylO6MqIOxsTGioqJQWFiIkJAQ6Onp4YMPPlBrDIwxpi1Ii+4O9HQs2nSMTDMsXrwYp0+fRkJCAszNzWFsbAwnJyekpqZiz549mDJlCgCgtbVV6c6Ivr6+eCclJCQE+vr6aovZzs4On376Kc6cOYPFixejra0NX3zxhdo+X0qcgEjE398f/v7+qKiowLp16zBjxgzxtbKyMkydOhWbN29GW1sb7r//fpw6dQoREREwNjaWLGZ7e3t8/PHHaGxsxA8//CBZHEz7aVPnhTtojDGmHoph7/fddx/Ky8vx0ksvYdGiReLr5eXl8PLywsaNG3H16lUsXLgQaWlpiImJgZWVlWRxT5o0CQkJCairq8M333wjWRzqpL40j3WhGH+o+FPB0dERjo6OGDlyJGJjY+Hl5SVFeIwxxvqBNiWfnGCzgUixKE/nxXlGjhwJY2NjTJ48Ge+99x78/PykCK9bg62tcAIyQG3btg02NjbQ09OTNI5vvvkG5eXleOihhzB69Gg+qTC1MDR8GaNHb4CxcYfUodyx4cPj4Or6D+jrj1faPnLkt2hvb4auruonOTLGGOtKV1cX27Ztg4ODA4yMjCSNZfv27bh8+TICAgJgZWU16PpXnIBIqLfKNmLECHWH060tW7Zg//79uOuuuzB69Gipw2GDRGvrUJSXAy0td760tNQuXRqP3FxALlfeXl7uhtxcaWJi0iJ6E8Cb0IZ+hpHRD7C0bIeurqnS9uHDf8Xo0e3Q1eXzBlO/3vpXTk5O6g6nWxs3bkROTg7mzJnDCQhjnXVuEIOtgTBp/K+aaUUHjbHOtKl+Nzebor6+6/ZLl0ajvPz6sTLG2I14ErqEuDPPmPbrqbOpTZ1QxnrC9ZtJQRP6V4P9Ai8nIKxXg72BMGlw55xpM22q35xgM8ZuBycgEuLOPGPajztojDGmXprQvxrsF3g5AWG9GuwNhEmDO+dMm2lT/eYEmzF2OzgBYYyxfsQdNMYYYzcz2C7w8ipYEiL6EsCXA7oDMmTIL7C07ICurtn/thjC0rIOBpq/OiobwLSpc25qClhZgdsME2lT/eYEmw1E+vrxsLT8DPr6MqlD6VHnpaqJ7DF6dBFGjBgcXfPBcZQDlCZ8QTc1mSktsUgE1NdbQeLn9zCmMX74ofvtmtD+GWNME7W1GaC+HpDLB243t7q681LV+igvHzdozgkD93+GDUjcaWLqMBjq2WA4RtY9bfq/DwkBli4Fxo1T3q5Nx8hYfxjsbYQTEAlpQuXThBgZY4xJ48EHpY6Asa40oe/SOUZNiFmVOAFht2SwNRAmjVmzhCurU6dKHUn/uftuYPRowMRE6kiYug2G79HBcIyMsdvHCYiENOELWhNiZJpv40agshJYsgRYvfr6j0JNDfDqq4C7OxAWJl2ct6u+HoiPF/7+6KOAnR3wxRfKZU6cAPbsAf70J8DTU90RMsaY9tCEvgvfAWHsFgy2BsLU4/vvgexs4NtvhbsfI0cqv97YCGzZIow118QEpK4OiIwU/p6RASQldS1z+LBQxtycExBtFxQk1HNtvsP3978LifewYVJHwhgbiPg5IBLShM78qlXCUBg7O6kjYYNBYyMQESF1FP1r+3YgNVXqKJi6/etfgLOzkGxPnQoEBAATJwqvTZjQfVKqaZ54Ali5UmjHoaFCW7a0BI4cARYvBhoapI6QDRaa0L/685+F/pWFhfC7JsSsSpyASEgTKltkpDB0xMlJ+F0TYmaaycBA6Lx8/TWwd6/U0fSPefMAMzPg6aeBlhapo2HqdPkycO4c8Mwzwt9vdO4ccOWKNHGp0rFjwI4dwCuvKG+vrhaS7vZ2aeJig8/1pW0HrtdfF/pXtrZSRyINTkAYYwPGpk1CIvL000Brq9TRqJ6dnTA0paAA+Oc/pY6Gqdu4ccKiA9HRUkfSf2bNAj75BPjjD6kjYUyzLpYOtgu8nIBIaPJkwM8PGDFC6kj6TkcHmDEDmD5d6kiYNnJxAV54AcjPB958U+po+sff/w6MHw+89RZw6pTU0TB1GjIE+L//Az78UJgLpI0WLRKGW/31r4BcLnU0bLBychL6Kpo0B8nAQLhIMXq01JGoBycgEoiKAj77TLhNnZICLFwoDDv56SepI+tZSYkwHOvSJWG1ngMHhO2vvNJ1OAFjd+LllwF7e6GjVlh4a++VyfonJlUyNATeeAPo6ADWrdOMoQJMdVavBnx8gKee6nt9vXIFeOkloLa2f2NTlXffBU6eFBaOYEydNm4E/v1vYXjTiRPAsmXAN98M7P7V2bPC3X8zM6CoCNi3D2huFrZ1dEgdXf/hBEQCSUnCEJPs7OvbDh0C0tOli+lmysuFxvDMM8rb336bJxYy1TI3B955R5gj0bm+9ebCBWDUKGF4y0AfvvXnPwt3P/ftAxITpY6GqdsHHwhDlD75pG/l33xTaBMTJwrvHeidEhcX4PnnhYsJVVVSR8MGk++/F84bOTnXtw30/lVOjnCB98a5U9euCdsGelu/E5yASGT1as27RT1unHBF+scfpY6EabuHHgKWLwd+/RXYtq1v7/n2W+DiReEKmLe30LkfyD75BDA2Fq5s19dLHQ1TJ1dXYSjeyy8LdfZmQkKEVbPq6oBnnwXc3ITJ3gNZVBRgbS3Ub8bU6Z57hDuMmnR3edo04S7NsWNSR6I+nIBI5MknhTsHn35687KlpcCZM/0f083o6wtjl59+mjtMrP+9844wXOnll4Xb0Tfz/PPCSjvu7sDp08LQxsWLB0bb6Y6zM/Dii8IKQTExUkfD1C0yUlii9q23rm/78cfrw1tv5Ows3Dnfs0eYO1hQIKwYt3ix8pVeKbS3Cyv5NDUpbzc2FpYe/vpr7Z3vwgamv/5VSNYTEm5etrISKCvr/5huxsREGFHy6KMD/w6+qnACIhE9PSAuTjgJ3ewW9TPPAB4eQHCwMAdDSosWCaucvPGGtHEw7ae4SlxcLHRkFI4fF76kKyu7vmfRIiAzE4iNFdZWT00VriyFh0uzzOmhQ8Jn92T9emG4ynvvCWN/2eBhairMlYiPF67UXr0KPP64MD+kp8TZ1/d6/ba0FOq3p6dQx6QYCpuYKCRHf/1r94tGrFolXI2+sf0y1t/09YWhin//u5CI9CY8XFgURKo2dKOHHhLuGm7eLG0c6sIJiIR8fQF//+tPSO5OR4ewIoKODrB1q3D169//lnbo1gcfCHdueJlF1t8iIwFHR+V5Ei+/LFzZmjRJuHrc+XkaBgbCySQ7G3jwQaENbd4MzJp1Ddu3/0ctcZeVAQ88AMyfL3x2Tw8eNDQUOqFtbX27G8q0y5o1wJ/+JCQgRkbAc88JV0JTUwEvL2DDBmEs+I309a/X74cfFur3hx/K4O9/D7777juQGsadnD0rnLvWrhXmB06YAMyZ033Zjz4CdLmnwdRswQIh+f3733su09EhnC8U54jJk4U7dlIN3dLRERZueOMNoV1pO/5akFhMjDDuT3Eb/ZFHlO906OsLX+C5uUJjqq4WxgO7uQG7dvV/fNu3A0uWKD9AatQooVFr+xOrmfRMTIQrxDdOxPvPf4CwMCHxePVV4Qrsli1dk/JRo4TVT06cEDpHdnavIzBwDby8vHDo0KF+ibeh4foVtW3bhAn1sbHCybAny5YJyUpfhpkxzTVzpjCPo7N33xW+0z08gNdeE4bchoUJSenrrwNjxgh3yzuvmDVyJPDVV8Jdtgcf3IaDB3/Bgw8+CB8fH/zRT1eHLl0Szk+ursDOncID1BISgLw8Yc7WkiVCJ+5Go0YJ8QcECEkWY+ry1lvADz9cX/BH0b+qrhZ+19e/fo6YNw84f14oM2NG90MhVUkuB778Erj3XuWEx91duMO/cWP/fv6AQEztXFyI9u27/vsHHxABRM88Q6SvL/zdxobok0+IOjqU37ttG5Gjo1BGR4foxRe/pvPnz6s8xsxMIh8f4XMAoogIookTr7/e3k40ZYrwWnGxyj+eDTKTJxMZGPT8+sqVQl1buvT6ttxcoiVLrtdRLy+iw4e7f79MRvTpp1+Qra0tASBdXV16/PHH6cKFCyqJv6ODKDaWyNpaiEVPjygsjKimRni9qEjYvmZN9++/cIHI0lIo89FHKgmJabiMDKJ5867X7+nTiQ4c6Ll8cnIyjRkzhgCQjo4OBQQEUGlpqUpiaW0liom5XkcNDYVzQl2dSnbPmMpMnky0Z8/13z/6qGv/ys6O6PPPhfOCgkxG9MUXRPb21/tXkZFfUWVlpcpjPHyYaObM62375ZeJ5sy5/vqVK0QjRwqvNTer/OMHDE5AJHD//cLJRaG9nWjFCqIPPyTKzxf+rqiYEycS7dih/P62NqGz4+l5mvT09MjQ0JDCwsLoypUrdxxbSQlRQIDQ+AAh2UlKIjp+nGjRIuWyKSlE48YRlZff8ceyQe7XX4m+/77n1ysqhHrYXQcsOZnIyen6SSMggKinfte1a9coJiaGTE1NCQANHTqUIiIiqLGx8bZj372byNX1epudP1+5fRMRNTQQxccT/fZbz/vZu1cok5Nz26EwLZScTDR27PX6tWJFzxd9rl27RlFRUTRkyBACQCYmJhQVFUXNd9CL2bFDOA8pPv+ee4jy8m57d4z1q8BAorS067+3txP5+1/vXy1bdr0uu7gI554bNTUJyfaMGSdIV1eXjI2N7/gcoVBUpNy/GzdOaN+//iqct2707bfC6y0td/yxAxYnIANUSgqRu/v1iurnR5SdrVymuLiM1qxZQwAIAI0bN45+/PHH2/q8piaiqCiioUOFzxsyRPj92rU7PRLG+p/ipGFqKtTfoUOF+ttTv+vs2bO0cuVKse04OzvTTz/9dEuf2flk4uQknEwYU7XO9dvYWLgD0VOfKDc3l5YsWSLWb3d3dzp48OAtfWZ2tnDeUdRvd3fhvMSYpktOFjr3vSX1RUUldO+994ptaMKECbRz587b+ryGBuGOuKGh8Hnm5sJF5NbWOz8WTcYJyACmuNNhYSFUWgMDoRLX1yuXS09Pp7vvvltsKN7e3nTs2LE+fYZcTpSQQDRq1PUryEFBfFeDaaaCguvDtQAiN7c22rWr517T8ePHxbbz2muv9ekzmpqaKDa2mExMrifrr7zSc2eQMVUpKCBatep6/V627FFK7iXr3bt3L3l4eBAA2nPjuJReNDY20iuvvEJeXpcJIDIxIXr9dSEJYkxbtLYK/aubXbQ6fvw4zZo1S+xfLViwgE6dOtWnz5DJhLvaimFdurpEISFE/TCqSyNxAqIBzp8neuSR68OinJxa6Ouvk0gul4tl2tvb6YMPPqBhw4YRALKzs6OWm9y7y8nJoQcfbBRPZq6uwnASxjTd/v1E06YR+fj8301PGu3t7RQfH3/TIYwdHR0UHx9P9vb2ZGMzmSwt5RQSIgwPY0ydfvuNaM2aZLFTtHTpUsrNze22bEtLC3333Xc33adcLqeEhARycHAgAOTq+jAFBQnnH8a0VWGhMERL0Q/y8Gil3btTlcq0tbXRe++9R5aWlgSAHB0dqb29vdf9njx5klaubBb36+kpnJfYdZyAaJCTJ4Xx5QsW/IMAkKenZ5fb6jU1NfTUU0/Rxx9/3ON+Ll68SEFBQaSrq0szZ24kOzvhLsiNE7IY03RtbUSbN8eLJw1DQ0N68cUXqaGh4Zb3deLECZo7d67Y4fP09KTDh4tVHzRjfSSTySghIYHs7OzEhRWCgoLo0qVLt7yv3bt3k6urq1i/fXx8KDMzU/VBMzZA7dsnLKyzYMFbBIAWLlzY5aJVVVUVPfHEE/Tll1/2uJ+SkhIKCAggHR0dmjPnIxozRpi/eMP1YvY/nIBoGLmcKDHxGxoxYoR40nnyySf7dNJR3CWxtrYmAKSvr09/+9vfqKam46bvZUxTXb58mSIiIsjQ0JAAkJWVFcXGxt70ChYRUXFxMa1YsULsmI0ZM4aSkpTvPjImpTup30VFRUr128nJqdchXYxpM2FY1kdkYWGhdNGqvvO4927f20rvvPNOp/f+nerr+VzRE05ANFRTUxNFRUWRsbGxuJpPb6udJCUl0fjx48UTTUBAABUUFKg5asakk5eXR8uWLRPbgIuLC+3uYczhlStXKCIiQmxfpqamFBMTQ9d4VQY2QOXn59Py5cv7VL8bGhooLCxMTFrMzMwoNjaWWgf7rFjGiKiurq5LUh8TE9Nt+1AMXRw1apS4BHZQUBCV80Tam+IERMOVl5dTUFCQ0mo+SUlJ4uunT58mX19f8XUPD48+T0ZkTBslJyfTuHHjxDaxYsUKOnfuHBEJw1ri4+PFcfCKk0l/PGuHsf7Ql/ptb28v3kEPCQmhCp7IxFgXeXl5Skn9pEmT6OeffxZfz8zMJB8fH/H1adOm0b4bH/LGeqVDJNVD55kq7dy5E88//zwKCgoAAP7+/rC3t8eXX36JtrY2mJubY/369QgPD4ehoaHE0TImrWvXriEmJgbvvPMOWlpaYGpqisceewynT5/G/v37AQAeHh547733sGjRIomjZezWXLt2Df/85z/xr3/9Cy0tLTAzM8Njjz2GkydP4tChQwCAadOm4b333sOCBQukDZaxAYyIkJSUhJdeegnl5eXQ0dHBfffdB0tLS2zduhXt7e2wsrJCdHQ0nnrqKRgYGEgdssbgBESLdHR04PPPP8err76KmpoaGBsbQyaT4bnnnkNERASsrKykDpGxAaWyshKRkZH46quvYGVlhbq6Ojg6OuKdd97BmjVroKOjI3WIjN02rt+MqUZ7ezs++ugjbNiwAY2NjTAwMICOjg4iIyPx4osvwtTUVOoQNQ4nIFqosrISX3/9Nezs7DB58mRMnz5d6pAYG9D27t2LK1euICcnB08//TQn60yrpKam4tq1a8jOzsYzzzwDCwsLqUNiTCOVlZXh+++/h4WFBby9veHh4SF1SBqLExDGGGOMMcaY2uhKHQBjjDHGGGNs8OAEhDHGGGOMMaY2nIAwxhhjjDHG1IYTEMYYY4wxxpjacALCGGOMMcYYUxtOQBhjjDHGGGNqwwkIY4wxxhhjTG04AWGMMcYYY4ypjb7UAWiK0tJSFBQUAADGjRuHcePGdSlTXFyMoqIi6OnpwdvbGyYmJl3KZGVloaqqCnZ2dnBxcYGBgYH4WmtrK7KzszF06FC4urp2G8fJkyehq6uLadOmdXmtrKwMHR0d3cbGWH+Ry+XIyspCdXU1dHR0MH369C5PEiciZGZmoq6uDmZmZvDy8oKenp5SmdbWVmRmZqK5uRmOjo5wdnZWer26uhplZWUYMWIERo4c2SWOxsZGnD17FlZWVl3aQEFBAcrKymBlZQUXFxcMHTpURUfPWO+ICFlZWTh37hwMDAzg5eUFOzu7LmVOnTqFoqIi2NvbY9asWV3aR0tLC9LT01FbW4tJkyZ1OUdUVlbizJkzGD16NCZNmtQljvr6epw4cQLDhw/vcv5QfPbw4cPh5eXF7YOpTVNTEzIyMtDa2gpjY2PMmjUL+vr6Xcqkp6ejvb0dY8eOxfjx47vsp6qqCvn5+dDT04OLiwusra2VXi8oKMCVK1fg4uLSbd+svLwcly5dgpOTE4YNG9ZtrLm5uWhqaoK1tTXGjh17+wfNBMR6VVRURPPmzSMASj/Lly+nqqoqIiJqb2+noKAgpdfNzc1p+/bt4n4yMzNpypQpSmWmTJlC2dnZYpmjR48SABoyZAjV19d3ieX48eMEgHR0dKi5uZmIiDo6OigpKYnmzJlDAGjSpEn9/C/C2HW//PILjR07VqleGxkZ0QsvvEAdHR1ERHT+/Hny9PRUKjNx4kTKzc0V97N161aytbVVKhMQEEBXrlwRy7z++usEgBYsWNBtLG+++SYBoPnz54vbSktLaeHChUr7dXR0pN9//71//kEYu0FKSgo5Ozsr1T8dHR3661//Su3t7UREdPHiRZo5c6ZSmXHjxlFWVpa4n4SEBLKyslIq86c//YlqamrEMps2bSIANGPGjG5jiYmJIQDk5eUlbisoKCBvb2+l/To4ONDevXv76V+EMYFMJqPXXnuNjI2Nu9S/H374QSy3Y8eOLnU/MDCQWltbiYioubmZnn76aTIwMBBfNzY2ptjYWKXPmz59OgGguLi4buNRtMH33nuv29ePHz9Oenp6pKurS8HBwar5RxjkOAHpRU1NDY0aNYqsrKwoISGBWlpaSCaT0Y4dO8jJyYkSExOJiOitt94iALRhwwZqa2ujyspKWrBgARkZGVFxcTE1NjaSvb09LViwgM6dO0dEQmUePnw4ubq6ip/37bffEgDS1dWlr7/+uks8L774Iunq6hIAys/PJyKi8PBwsrKyoqCgIJo+fTonIExt9u3bR/r6+jR16lTav38/yeVyunr1KsXFxZG9vT2dOXOGiIh8fX1pyJAh9PPPPxMRUVpaGo0YMYI8PDxIJpPRsWPHSEdHh55//nmqr6+ntrY2+vTTTwkAPfvss+LnPfHEE6Srq0t6enp04cKFLvFMmzaNdHV1acyYMUREJJfLycvLi1xdXSkzM5OIiM6ePUsuLi5kaWmplNwwpmpHjhwhAwMDmjp1Kh09epTkcjk1NTXR559/TiNGjKA//viDiIiWL19OhoaG9OOPPxIR0cmTJ2nUqFE0adIkam9vpxMnTpCuri797W9/o8uXL1N7eztt3bqVdHR06MknnxQ/7+mnnyZdXV3S0dGhoqKiLvF4eXmRrq4u2djYEJFw8crFxYWWLFlCxcXFRERUVlZGHh4eZGNjI15AYKw/vPDCCwSAHn30USopKSEiouLiYnrooYdo3rx51NHRQRUVFWRiYkKTJ0+msrIyam9vpw8//JAA0Pr164mIKCIigkxNTWn79u0kk8movr6eAgMDCYBSIm1tbU26uro0d+7cLrGcO3eOdHR0SFdXl8LDw7u83tzcTJMmTaJVq1aRk5MTJyAqwglIL1599VUCQLt37+7yWktLCxEJWbyDgwPNnDlT6fXc3FwCQJGRkSSXy2n37t106dIlpTKRkZEEQLyT8tZbb5GRkRHdc889tHz5cqWyis9ZvXq1UkzV1dXU1tZGREQPPPAAJyBMbby9vcna2lqsvzdStI/Tp08TAHrllVeUXn///fcJAKWmplJTUxN9//33JJfLlcp4eXmRh4eH+Luvry/5+vqShYVFl6tU2dnZBIBWr15Nenp64tWx48ePU15enlLZzz//nADQwYMHb/vYGbuZ+fPnk6WlZZfvfaLr7SMvL48A0PPPP6/0+r///W8CQD/99BO1t7fTjz/+SDKZTKnM3LlzaeLEieLvS5cupfnz55OtrS394x//UCp75swZsX0AEJPv7OxsqqurUyqruJNSWFh4+wfPWC9KS0tJT0+PVq5c2eV7Xy6Xi9/f//jHPwgApaenK5Xx9fUla2tram9vp7Kysi53tMvLy8WLwkRE9fX1Yv3X0dERLwQrvPHGG2Rvb09z5syhlStXdol3/fr1ZGZmRmVlZZyAqBBPQu/FTz/9hEmTJmHJkiVdXjMyMgIgzLuorKzE6tWrlV53cXGBq6srjh49Ch0dHSxZsgQ2NjZKZRS/19bWAhDmkDg6OiIwMBApKSnidgA4dOgQKisr8eKLL4plAWD48OFK80gYU4cLFy4gLS0NDz74IGxtbbu8rmgfx48fB4Au7eO+++4DABw9ehTGxsa4//77oaOjo1TGxsZGqQ0UFxfDxcUFK1euxHfffadUdtu2bXBwcMDDDz8MmUyG8vJyAIC3t3eX8fCd2x1jqlZXV4eDBw9izZo1Xb73gb63j+PHj0NfXx/33nsvdHWVT9ed20dJSQnGjx+PVatW4ZtvvlEqu337dtjY2OCxxx4TywKAu7t7l/laaWlpsLKygoODw60eNmN9snPnTshkMoSFhXX53tfR0YGhoSEA4NixY3B0dISXl5dSmdWrV6O2thZnz57F6NGjsWDBAqXXFfM/FO1DUd8ff/xxWFlZYdu2bUrlv/vuO9x7771wdnYW+1YKmZmZiImJwcaNGzF69Og7Om6mjBOQXuTn52PixIm9ljl//jwAYMyYMV1eGzNmjPh6d9LS0mBmZiZOmC0pKcHYsWOxcuVK6Orq4j//+Y9Ydtu2bXB1dcXdd9+N4cOHiw2KMSnk5+cDACZMmNBruZ7ah4ODAwwMDHpsH+3t7cjMzBQny3Z0dKC8vBxjx47FmjVrcPz4cRQWForlv/nmG6xevVpsS51PIjdKS0uDjo4OpkyZ0vtBMnabzp49CyK6afuoqKgA0LV9WFtbw9TUtMf2IZPJkJGRAU9PTwDCJHbF+WPNmjXIy8tDVlaWWP7bb7/FqlWrxMm7nc8fKSkpeP/993HfffchNTUVX375JYyNjW/pmBnrK8X5oy/9q+76Vo6OjgAgXmjqLD09HQAwdepUANfPBxMnTuxyASs7OxvZ2dlYs2YNxo4dq9Q22tvb8Ze//AVTp05FWFhYH4+O9RUnID1oa2tDe3s7hg8f3mu5y5cvAwDMzc27vGZubo66urpu33fq1Cn88MMPeOGFF8SrYcXFxXBycoKlpSX8/Pzw7bffAhA6X9u3b8f9998PAHBycuq1g8VYf2tqagKAbq/u3qin9qGjowMzM7Me28fHH3+MqqoqvPzyywCEE1FHRwecnJywZMkSmJubIykpCQCQkZGBoqIirFmzBk5OTgB6TkAuXryIjz76CA899JBYljFVU7SPm50/6uvrAfR8/lC0n862bNmCsrIysX1cvHgRLS0tcHJywsKFC2FtbS1e5c3KykJeXp7YwdLR0enSPlJSUrB582bs3bsXrq6ufPeD9atr164BuHn7uHz5co9tQ/F6Z3K5HBs3bsTYsWPxyCOPABDOB3p6enB0dMSaNWtw6tQp5OTkABAu7tra2sLHxwdOTk5obGxETU0NACA2NhZZWVmIj4/vsiodu3OcgPTA0NAQRkZGuHjxYq/lFB2w7hpCXV0dRowY0WV7VVUV/vznP8Pd3R1///vfAQhXtMrKysRO0f33348DBw6goqICe/fuRXV1NdasWQOAExAmPTMzMwDCUKze9NQ+ZDIZGhoaum0fBw8eREREBJ566inMnTsXwPWEwsnJCUZGRlixYoV4FWvbtm2wt7fHvHnzYGFhASsrq27vEDY1NeGBBx7AkCFD8Pbbb9/aATN2C261fXSXiNfW1nbbPvbt24fnn38ef/nLX+Dr6wtAuX0YGBhg1apV+O6770BE2LZtG4YNGwZfX18MGTIEI0aM6NI+3n77bRQUFKCiogK2trZYuHDhTc99jN0uRQJRVVXVazlbW9tu24bifNJ5OWsAeO2113Dw4EF8/PHH4l28kpISODg4wNDQEH/6059gZWWFpKQkEBG+/fZb3HvvvdDT01O6gFVYWIioqCg888wzmDFjxh0dL+seJyC9cHNzQ15eHoioxzKKK0XddXhKS0u7XEkqKCiAj48PDAwMkJKSIjaQyspKtLW1iWtL33fffdDX18f27duxbds2TJo0Sbyd6OTkxEOwmKQUzyDIy8vrtVxP7aOiogIymaxL+9ixYweWLVuGwMBAvP/+++J2xfsV7WPNmjXIyspCVlaWOH5XMUa+uwS9uroaS5cuxdmzZ7Fnzx6+wsv6lYuLC/T09HDmzJleyynq4blz55S2V1ZWorW1tUs9/fXXX+Hv749Vq1Zhy5Yt4vbu2kdRUREyMjLw3XffYdWqVeJcwd4uYJmamuK5557D1atXkZqa2ufjZexWuLm5ARCeq9Ebe3t7lJWVdemDlZaWAoBS++jo6MBzzz2Ht99+G19//TWWLl0qvqYYXQIABgYGWLFiBb7++mulu+fA9fZTUlKCqKgoNDc349dff4WXl5f4U1lZiZ9//hne3t539o/AOAHpzZo1a1BcXIwff/yxxzKjR4+Go6Mjvv/+e6VGcurUKRQUFGDevHnith9//BEzZ87E6NGjsW/fPqXbj4oGpWgAlpaW8PHxwTfffIP//ve/WLVqlVh2zJgxqK6uFm9jMqZuw4cPx8KFC/Hdd9+J49i7M2fOHADCJNgbKX5XtA+ZTIZXX30V9957L5555hl8/vnnSpNuS0pKYGZmJk4uXLJkCUxMTPDSSy+hvLxcnLQLCO3jxoTn2LFjmD59Ourr63Ho0KFuH9LGmCqZmZlh2bJl+O9//9vr3ep58+Z1me8HQPx94cKFAK4PK1m+fDmefPJJfPXVV0pDQkpLS2FkZAR7e3sAgK+vLywsLLB+/XoUFRV1aR+K8825c+cgk8mUPru1tfUOjpyxm/P394eRkRHefffdXi/wzpkzBxUVFTh8+LDS9u3bt2PEiBHiHKuLFy/Cz88PX375JZKTkxEQEKBUvqSkRGkuyX333YeioiKsX78eVlZW4iT2kSNHQl9fHyUlJVi1ahViYmLw6KOPIiAgQPwxMTHBmDFjxKSF3QHpFuAa+K5evUoTJ04kMzMz+uSTT+j8+fNUW1tLf/zxBz333HNUUVFBREQff/wxAaCwsDCqqqqinJwcmjp1KllZWVFVVRW1t7fTSy+9RDo6OuTr60vHjh2jEydOiD/V1dWUkJCgtCQv0fWlGAFQRkaGuH337t0EgLKysqilpUXcz5/+9CcaM2aM+LviQVeM9YeTJ0/SkCFDyM3NjXbs2EHV1dV08eJF2rt3r9KyomvWrCEDAwP66quvqK6ujnbt2kUWFhbiAwUvXbpEixYtIgD04osvKrWNEydOUGtrKz3yyCM0ZcoUpc9/4IEHCADZ2toqPbPghRdeIDs7OyISlvs1NDSkyZMn0759+5T2W15eroZ/JTZYZWdn09ChQ2nMmDH0ww8/UHV1NZ0/f56Sk5PpmWeeEZcfDQ4OJj09Pfr888+prq6OUlJSaNiwYTRnzhySy+VUU1NDy5YtE88xndtHc3MzPfHEE0pL8hIRPfroowSArKysxKXaiYTl5S0sLKipqYkcHBxo2bJldObMGaqrq6PDhw/T5MmTycLCQukhh4ypmuLBsX/+85/p5MmTVFdXR+Xl5fT555/TZ599RkREtbW1ZGNjQ87OznTq1Cm6dOkSbdiwgQCIDxo8ePAg2dvbk4WFBX3zzTdKbUPxLCoTExOKiooSP7ulpYXMzc0JAD3++ONKcY0bN47++te/9hg3L8OrOpyA3ERFRQX5+/uLDwBU/EybNk0pKXj11VdJT09P6WmeR44cISKiAwcOdHmS+o0/W7ZsoejoaDIxMVH67JqaGtLX16fx48crbVesHZ+cnEyFhYU97pdPIKy/HTx4kKZMmaJU7/T09Oiee+6hxsZGIiJqaGigpUuXKpWZM2eOmGwrTig9/RQXF9PcuXNp1apVSp/9/fffEwAKCQlR2v7++++Tjo4OFRcX97rfl156SS3/RmzwOnr0KHl6eirVO11dXVq6dCldvnyZiIgaGxvpz3/+s9I5ZuHChVRZWUlE1ztqPf3k5uaSn58fLVmyROmzd+7cSQC6dJYUF7bq6uro0KFD5O7urrQ/R0dHOnDggFr+fdjg9u6779KwYcOU6p+FhQW98cYbYpm0tDRydHRUOr+8+OKLYgLv4eHRY9vw9PSkS5cuEQD64osvlD77oYceIgC0a9cupe2LFi2ipUuX9hgzJyCqo0PUy/0vJqqrqxPH8zo7O4u3um9UW1uLkpIS6Ovrw93dHfr6+gCEsYmNjY097nvo0KGQy+Voa2uDhYWF0msNDQ3Q09ODqampuE0ul6OhoQEmJibQ19dHQ0NDt/u1tLTsssY2Y/2hpKREHCY1fvz4LvUYEMbh1tXVwdzcXGl50ubmZrS0tPS4bwsLC1y9ehV6enowMTERt8tkMly5cgVDhw4VV5IDhBXsrl27BjMzs17b3ZAhQ3ipUaYWZWVlqK6uhrGxMUaNGtXtyj7l5eW4dOlSl/bR0tKC5ubmHvdtbm6Oa9eudWkfivOEsbExhgwZIm5XtA8LCwtxmGNeXh7Onz+PESNGwM3NrcszRxjrLzKZDDk5OaiursaoUaPg5OQkPgekc5n29nY4Ojoqrb545cqVLsMIFRR9p4aGBpiamio9M03Rrm5sB4CwQpdcLhcXkuisoaEB+vr6Sm2N3R5OQBhjjDHGGGNqw5c5GGOMMcYYY2rDCQhjjDHGGGNMbTgBYYwxxhhjjKkNJyCMMcYYY4wxtfl/1fCW452rXdQAAAAASUVORK5CYII=\n",
"text/plain": [
"<PIL.Image.Image image mode=RGBA size=800x200 at 0x7FECA82B0390>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"C02M1 = Chem.MolFromSmiles('c1cc(O)ccc1-c2c(N)ncnc2')\n",
"C02M2 = Chem.MolFromSmiles('c1ccccc1-c2c(N)ncnc2')\n",
"C02M3 = Chem.MolFromSmiles('c1ccccc1-c2c(N)ncnc(Cl)2')\n",
"C02M4 = Chem.MolFromSmiles('c1cc(Cl)ccc1-c2c(N)ncnc2')\n",
"Draw.MolsToImage([C02M1, C02M2, C02M3, C02M4],\n",
" molsPerRow=4, legends=['C02M1',\n",
" 'C02M2',\n",
" 'C02M3',\n",
" 'C02M4'])"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "smoking-draft",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0.7033252807348217\n",
"0.6850880906633667\n",
"0.7033252807348217\n",
"#############\n",
"0.6666666666666666\n",
"0.275\n",
"0.6774193548387096\n"
]
}
],
"source": [
"print(AtomAtomPathSimilarity.AtomAtomPathSimilarity(C02M1, C02M2))\n",
"print(AtomAtomPathSimilarity.AtomAtomPathSimilarity(C02M1, C02M3))\n",
"print(AtomAtomPathSimilarity.AtomAtomPathSimilarity(C02M1, C02M4))\n",
"print('#############')\n",
"print(get_tc(C02M1, C02M2))\n",
"print(get_tc(C02M1, C02M3))\n",
"print(get_tc(C02M1, C02M4))\n"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "difficult-central",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAyAAAADICAYAAAAQj4UaAAB4d0lEQVR4nO3dd1wT9/8H8FfCHgooDlBx4sAtbpx1K9RJ7RC1Vq21La5WtF8talvFajVaa0s3aoe4Ku7iqFsZTlw4QIYC4mKv5P374345CUNZySXwfj4eeaB3l7t3IDfenykjIgJjjDHGGGOM6YBc6gAYY4wxxhhjVQcnIIwxxhhjjDGd4QSEMcYYY4wxpjOcgDDGGGOMMcZ0hhMQxhhjjDHGmM5wAsIYY4wxxhjTGWOpA2CG6/79+wgJCUFISAhq1aqFrl27ol+/flKHxRhjjDHG9BgnIKxE7t27h1OnTiE8PBzh4eG4cuUKUlNTxfV2dnZ4+vQpvvvuO8ycOVPCSBljjDHGmD6T8USErKDMzExcuHABoaGhCA0NRUhICO7cuVNou7p166Jr167o2rUrsrKysHz5chARfv/9d0ycOFGCyBljjDFBeno6Lly4gE6dOsHKykrqcBhj+XANSBWXmZkp1mqcPn0ap06dwsOHDwttV6tWLXTr1g2urq7o1asXunTpAhsbG41tHBwc8OGHH2Lq1KmoU6cOhgwZoquPwRhjrAp72b3Mw8MDzZs3x+rVqyWOkjGmxjUgVYhKpcKNGzfEi3N4eDgiIyORl5ensZ2FhQU6deqEXr16wc3NDa6urnB0dCzRMebPn49Vq1ahWrVq+O+//9CpUydtfBTGGGNVWExMDEJCQnD+/HmEhITgwoULSEtL09hGJpOhYcOGiI2NhVKpxIYNG/Dhhx9KFDFjLD9OQCqxa9euiSVCRfXbAAC5XI6WLVtqJBstWrSAsXHZKseICFOmTMHvv/+OWrVq4fTp03B2dq6Ij8MYY6wKSkhIwMmTJ8WCs6LuZTKZDK1atYKrq6tYU9+uXTuYmJhg69atePvtt0FE+O233zBp0iSJPgljTI0TkErkzp072LNnjzgy1b179wptY2triy5duoh9N7p27Yq6deuW+ljx8fGoVasWTE1NC63Lzc3F66+/joMHD6Jp06Y4ffo06tSpU6bPxBhjrOp49uwZTp48qVF4VlSz4Lp166J3795iwVm7du1QvXr1Yve7ceNGfPjhhzAxMcGePXu4iTBjEuMEpJK4cuUKtmzZglWrVonLbG1t0atXL7FEqDRNqfJLTEzEiRMnxNKnq1evIiUlBceOHSt22N2MjAwMGDAA586dQ+fOnXHs2DFYW1uX9eMxxhirZHJychASEqKRbNy6dQtKpVJjOzs7OzHRKM+9zMfHB19//TU3EWZMD3ACUknMmDED/v7+aNu2LaZNm4auXbuiQ4cOMDMzK9V+cnNzcenSJXH0q5CQENy6dQsqlUpjO0tLS/z444945513it1XcnIy3NzcEBkZiWHDhmH37t0wMTEp0+djjDFm2Ao2C7506RLS09M1tjE3N9dINFxdXdGyZUsYGRmV+/gFmwifOnUKzZs3L/d+GWOlxwlIJZCXlwcHBwckJyfj0qVLaN++fYneV9LSp/w3BHVfkZKWPt27dw89e/ZEYmIi3nnnHWzevBkymazUn5Exxpjhun79Orp3716o70aDBg00mgW7urq+tClVUR4/fiwWmIWGhiIzMxNHjhwpctvc3FyMHDkSBw4c4CbCjEmIE5BK4NixY3jttdfQokUL3Lx5s8htiAjXr19/ZelTRXdKB4Dw8HD069cPaWlpWLhwIZYvX17mfTHGGDMsycnJcHBwgJmZGXr27Ilu3bqhS5cu6NKlCxwcHEq1r8zMTFy8eFEj4Sg4T5WRkRGePXtWbLPfjIwMDBw4EGfPnuUmwoxJhOcBqQR27NgBABg9enSR6zds2IDVq1fj/v37hdZVVKf0l3F1dcXu3bsxfPhwrFixArVr18bs2bMr9BiMMcb0086dO5GXl4fhw4dj9+7dJX5fVlYWwsLCNOb3KGpwlTp16qBPnz4aHdJfllBYWloiKCgIvXr1QlhYmFgjUtSgKowx7eAExMCpVCoxARk7dmyR28jlcty/f7/COvKVxWuvvYbffvsNEyZMwLx581CvXj14enrq5NiMMcak86p7FPBinqr8Ewnevn0bubm5GttV1OAq9vb2OHDgAHr27ImjR4/i3XffxZYtW7iJMGM6wk2wDNzp06fRq1cvNGnSBHfv3i1ym4SEBMTExKBDhw6Sl/B8/fXX8PHxgampKfbt24eBAwdKGg9jjDHtefLkCerWrQuZTIbExETY2tpqrL98+TJ8fX0RGhqKBw8eaKyTyWRwdnYWa+e7dOmCjh07lnpwlZe5evUq+vTpg2fPnsHHxwd+fn4Vtm/GWPG4BsTAvar5FSCMl17RzarKav78+UhISMDatWsxduxYHD9+HB06dJA6LMYYY1oQFBSE3NxcDB06tFDyAQAmJibYvXs3ZDIZXFxcNGo2OnbsCCsrK63G17ZtW+zcuRPDhg3DypUrUadOHcyZM0erx2SMcQ2IQSMiODk5IS4uDmfOnEGPHj2kDqlEiAiTJk3C5s2b4ejoiNOnT6NRo0ZSh8UYY6yCeXh4YO/evfjpp58wderUQuuVSiV27dqFLl26oGHDhhJEKPj777/xzjvvgIgQEBAALy8vyWJhrCrgBMSAhYaGomvXrmjQoAHu379vUG1Xc3NzMWLECAQHB6NZs2Y4ffo0ateuLXVYjDHGKsjz589Rp04dKJVKJCQkoGbNmlKH9FLffPMNPvnkE5iYmGDv3r0YPHiw1CExVmnJpQ6AlZ26+dWoUaMMKvkAhGr37du3o2PHjrhz5w48PDwKDQnMGGPMcO3btw/Z2dno06eP3icfADBv3jzMmzcPubm5GDduHC5evCh1SIxVWpyAGLDAwEAALx9ZRJ9Vr14d+/btQ6NGjRASEoLx48cjLy9P6rAYY4xVgJKMfqVvVq1ahUmTJiE1NRVDhgzB7du3pQ6JsUqJm2AZqMuXL6NDhw5wcHBAXFwc5HLDzSXv3LkDNzc3JCUlYerUqfjpp5+kDokxxlg5pKeno1atWsjJyUF8fLxBzTaem5sLDw8PHDp0iGdLZ0xLDPeptYpTlyyNHDnSoJMPAGjWrBn27t0La2tr/Pzzz1iyZInUITGmt548ARYsAE6c0Fx+9iywYQMQFwcsWwYUnHc0LExYnp2tu1hZ1bV//35kZmaiZ8+eBvfwbmJigm3btqFTp064e/cuPDw8kJaWJnVYjFUqhv3kWoVt374dgGFVbb9Mly5dsHXrVhgbG2Pp0qX49ttvpQ6JMb307BmwciXw5pvA8+cvloeHA5s2AUSAr6/wM7/iljOmDYbY/Cq/atWq4eDBg2jevDlCQ0O5iTBjFYwTEAN048YN3LhxA7Vq1UK/fv2kDqfCDB8+HL/99htkMhnmzJkj3sAYY4U5OQFffSV1FIwVlpmZiX379kEmkxlsAgIAtWrVwv79+1GnTh3s378fkydPBrdaZ6xicAJigNS1Hx4eHjA2rlxzSU6YMAFffPEFlEolJkyYgJMnT0odEmN6adUq4NtvgYgIqSNhTNO///6LtLQ0dO3aFfXr15c6nHJp2rSp2ET4jz/+wKJFi6QOibFKoXI9vVYRhl61/SqfffYZ7t27h99//x2XLl1C7969pQ6JMb3Tvj0wZQrwwQea/UHUo1l/+CHg4PBi+Y0bwk+lUncxsqqpst2jOnfujN27d2PYsGFYvnw5ateujVmzZkkdFmMGjWtADMydO3dw+fJl2NnZYdCgQVKHoxUymQxvv/02VCoVfvzxR6nDYUxvffklEBkJBAS8WJabK/xs1gzo0OHFS10QzS1ImDbl5ORgz549AIBx48ZJHE3Fee211/Ddd98BAD755BPs27dP4ogYM2xcA2Jg1CVL7u7uMDExkTga7dm5cycAYMyYMRJHwpj+srMTmmL5+ADe3sIyW1vh55w5QKNGL7bt1g0IDAQqWatNpmeOHDmCZ8+eoWPHjmjcuLHU4VSoqVOn4s6dO/j555+RmZkpdTiMGTSuATEwlalqW6lU4v3338fevXuhUqk0llemz8mYNnl5Ac2bA+vXSx0JY5XrHlUUOzs7PH78GHv37pU6FMYMGicgBuT+/fsICwtD9erVMWTIEKnDKbfjx4/jxx9/xCeffKIxl8mZM2eQmJiIFi1aoF27dhJGyJj+k8mAn38Gnj4t2fZPnwpD+L7+OjB9+osmW4yVV15eHv755x8AgKenp7TBlFNxQ+5W9gSLMV3hBMSA7Ny5E0SEYcOGwdzcXOpwyk19IS/YzEq9fPTo0TqPiTF9V726kDiYmr5Y5uwsjIg1erQw0eCECUBWlub7UlKE5XfvAh99BAQFARYWwO7duo2fVV7Hjx/H48eP0bp1azRv3lzqcMosLy8PjRs3xvjx45GRkSEuz18IOHjwYAkjZMzwcWtgA1KZSl6USqXYzyP/5yGiSvU5Gato9vZA06aAp6fQ96NXL2H5tGkvttm8ufD7BgwQXvmZmwNmZtqLlVUtleXafeTIEcTFxSEiIgKWlpbicnUh4IgRI2DGJw5j5cI1IAbiwYMHOHv2LKysrDBixAipwym3s2fPIiEhAU2aNIGrq6u4PCQkBHFxcWjSpAk6d+4sYYSM6a9Nm4C9e4XajrKKiABOnwaGDq24uFjVVVyhkiFSJ1IFR/GqLAkWY/qAExADsWvXLqhUKgwZMkSjRMZQFdfMSr181KhRug6JMYNw8yZw7RpQqxbQr1/Z9hESAsyeDfz9N1CJB9NjOlRZ+u7l5eVh165dADQTjfyFgMOGDZMqPMYqDU5ADERlKnkhInE294LNr7Zt21ZoOWPshf+/FOD11wEjo9K/f9cuoSlWixbAhg3Av/9WbHysaiquT5+hOXnyJJKTkwslUupCwKFDh1aKQkDGpMZ9QAxAUlISTpw4AXNzc3h4eEgdTrmFhoYiLi4ODRo0QPfu3cXlly5dQnR0NBo0aIAePXpIGCFj+kudgJQ1R+/bF7h8+cX/1fOGMFZWxRUqGaJXDY5i6J+PMX3BCYgB2L17N5RKJYYNG4Zq1apJHU655W9mJZPJXrmcMSa4exe4eFGYgHDgwLLto0YN4cVYRQkJCUF8fHyhPn2GRqVSFZlo5C8EdHd3lyo8xioVboJlACpbyUtxzay4+RVjL/f/TdMxYgT33WD6o7IMnX7mzJkiB0dRFwIOGjSoUhQCMqYPOAHRc7m5ubhz5w7kcnmxkw9++eWX+Oabb3QcWdlcvnwZUVFRcHBwQO/evcXlERERiIyMLLScMfZCeZtfMVbRKlPfPW5+xZjucAKi50xMTNC/f3+oVCrMnj0bKpVKY31YWBh8fX3x6aef4o8//pAoypJTX8hHjhypMfu5evnrr7+usZwxJoiJAc6fB6pV46Fzmf5Q991zcHDQ6NOnplQq8euvvxa6d+mb4vqxPHnyBEePHoWpqSlGjhwpVXiMVTr8pGcAfHx8ULt2bQQGBmLGjBka6zp37owffvgBRIRJkyaJwwfqq+JKkriEibGX++cfgAgYNkyYQJAxfdCsWTMMHDgQz58/R0RERKH106ZNw3vvvYc5c+ZIEF3JqQdHcXJyQrdu3cTle/bsQW5uLvr37w9bHrGBsQrDCYgBaNasGfbs2QMrKyv89NNPWLZsmcb6adOmwdfXF0qlEu+88w5Onz4tUaQvd+PGDVy/fh21atVCv3wTGNy8eRNXr16Fvb09+vfvL12AjOkxbn7F9JG1tTXq16+PjIwMDBgwAJGRkRrrZ8yYASsrK6xfvx5Lly6VKMpXy9+PpajBUbhwjLGKxQmIgejatSu2bt0KY2Nj+Pr64rvvvtNYv2TJEnz00UfIzMzEyJEjcePGDYkiLZ76Qu7h4QFj4xcDsKlnzy24nDEmSEpKQkJCOwwc+AV4DjSmT2QyGX788UcMHToUjx49wvDhw5GYmCiuz3/vWrJkCTZs2CBhtMUrqh9LSkoK/v33XxgbGxt8B3sGBAYCTZsWnvtowQLhdeoU0Lp14fcdOwZ07KibGKsSTkAMyIgRI/D9998DAGbNmiU+uKspFAqMGTMGjx8/xuDBgxEbGytFmMXi5leMlc3OnTsRGXkV1tYXwIPwMH1jYmKCHTt2oEePHrh79y7c3d2RlpYmrh8xYgR+++03yGQyzJ49u9C9S2qXLl0SB0dxc3MTl+/btw/Z2dno3bs37O3tJYyQVYSUFKEv3cyZQGbmi+XJycIrPR24fbvw+9LSgDt3dBdnVcEJiIGZOnUqli5dCqVSiQkTJuDUqVPiOiMjI2zZsgW9evVCXFwchg8fjqdPn0oY7Qt3797FpUuXYGdnh0GDBonL7927hwsXLsDOzg6DBw+WMELG9Bcn6UzfWVpaIigoCC1atEBYWBhGjhyJnJwccf2ECROKvXdJLX/zq6IGR+HzrvJo0waoWRP48kupI2GcgBigzz//HB9//DEyMzMxatQojeZWFhYW2Lt3L9q1a4eIiAiMHj0aWVlZEkYrUF/I3d3dYZJvAgN1SdiIESM0ljPGBI8ePcJ///0HMzMzeHh4SB0OY8Wyt7fH/v37UbduXRw9ehTvvvsuiEhcv3jxYvHepU9NhYtqfpWeno4DBw5ALpdzAlKJyGTA998Dq1cDly5JHU3VxgmIgVq7di3Gjh1bZHMrGxsb7N+/H05OTjh+/DjGjx8PpVIpYbTc/IqxstqzZw/y8vIwYMAA2NjYSB0OYy/VpEkTHDp0CDY2Nvjzzz+xcOFCjfUKhQJjx47FkydPMGjQIMTExEgUqeD69eu4desWateujb59+4rLDx48iIyMDPTo0QN169aVMEJW0Tp1AiZMALy9hZEF88vNBWrU0Hy9/bY0cVZ2nIAYKCMjI2zevBm9e/dGXFwchg0bptHcql69ejhw4ABq1KiBoKAgfPTRR5LFGhMTg9DQUFSvXl1jMsWYmBicP38e1apVw1Ce2ICVw4EDB+Di4oLZs2dj4MCBRQ4Haqg4SWeGpl27dti1axfMzMywcuVKrF27Vlwnl8uxZcsW9O7dG/Hx8ZI3Fc4/N5WRkVGh5XzeVU4rVwLXrgGbNmkuNzEBLl/WfK1fL02MlR0nIAbMwsICe/bsQbt27XDt2rVCza1cXFywb98+WFpa4ocffsDy5csliXPXrl0gIgwbNgzm+SYw+Oeff4pczlhJxcbG4q233sKIESNw48YN/PTTTzhy5Ag6d+6MhQsXIjU1VeoQy+Xp06cIDg6GsbExT4LGDEr//v3x+++/Qy6XY968edi8ebO4ztzcHHv27EH79u1x7do1jBo1SrKmwkUlGllZWdi7dy9kMhnGjRsnSVxMu+ztgRUrgIULhU7m+TVooPni8Qe0gxMQA1ewudWbb76p0dyqe/fu+Pvvv2FsbIxFixbhl19+0XmMdnZ2aNeuHUaNGqWxnEuYWFk9f/4cs2bNQtOmTfH333+jWrVqUCgUuHfvHry9vaFUKuHn54dGjRph3bp1yMvLkzrkMtm3bx9yc3PRr18/1KxZU+pwGCuVN998EytXrgQRYerUqQgODhbX2djYYN++fWjYsCFOnDghSVPhW7du4fLly6hRowZee+01cXlwcDBSU1PRuXNnNGjQQKcxMd2ZOhVwcgJKM39zejpw/z5goLcU/UKsUrh27RrZ2dkRAJoxY0ah9f7+/gSAjIyMaNeuXboPsIAHDx6QXC4nS0tLSk9PlzocZiDy8vJIoVBQzZo1xe/z9OnT6cGDBxrb3bhxg4YOHUoACAC1bNmSDhw4IFHUZTdq1CgCQN9//73UoTBWZnPnziUAVL16dbpw4YLGumvXrlGNGjWKvXdp08WLF2nQoEE0ZcoUjeWTJk0iAOTn56fTeJh2/fQTUceOmssiIohMTIjee4/o4EHh3wUFBRFZWxOdOkXUuzfRRx8RdepElJysm7grK05AKpHjx4+Tubk5AaDly5cXWv/5558TALKwsKDTp09LEOELGzduJAA0atQoSeNghuP48ePUsWNHMano3bs3hYeHv/Q9QUFB1KRJE/E97u7udPfuXR1FXD6pqalkYWFBRkZGlJCQIHU4jJWZSqWiiRMnEgBycHCgqKgojfXnzp0jS0tLAkBfffWVNEH+v+zsbLK1tSUABnOtYCVz7x7Rv/8WXr5jB1FoKFF8PFFAQOH1MTFEmzcLP9VJx/z5RFu3ajfeyo4TkEpm9+7dZGRkRDKZjH799ddC6z/88EMCQDVr1qQbN25IEKFgwIABBIA2b94sWQzMMNy7d4/c3d3FJKJRo0YUFBRU4vdnZ2eTQqGg6tWrEwAyNTUlb29vev78uRajLp/Y2Fh69913CQD16NFD6nAYK7ecnBwaMmQIAaBmzZpRYmKixvqgoCAyNjYmmUxGv/zyi0RREh08eJAAUPv27SWLgelGTg5Rr15EZmZET5+W7r3jxhGdP6+VsKoMTkAqoR9++IEAkImJCe3fv19jXV5eHo0ZM4YAUP369SkmJkbn8SUlJZGRkRGZmZnRs2fPdH58ZhhSUlLI29ubTE1NCQBVq1aN/Pz8KCMj46XvS09PpytXrhRaHh8fT9OnTye5XE4AyNHRkfz9/UmpVGrrI5RKWFgY+fr6kouLi5hs1apVi5o0aUJZWVlSh8dYuaWkpFCnTp0IAHXp0oVSU1M11utDU+Fp06YRAFq2bJkkx2e6NWAAEUD0++8lf8+mTURvvKG9mKoKTkAqqUWLFhEAsrS0pDNnzmisS01NpS5duhAAGj9+vM5j+/nnnwkAjRgxQufHZvpPqVSSv78/OTo6EgCSy+U0ffp0io+PL9H7Fy1aRDKZjDw9Pen+/fuF1oeGhlLPnj3Fh/zOnTtL1iTx2rVr9MUXX1CHDh3EeACQjY0NjR8/nho0aEAAaObMmZLEx1hFS0pKoubNmxMAGj58OOXm5mqsX7x4sWRNhXNzc8ne3p4A0LVr13R6bCaNjRuFBMTDo2TbKxREkyYRZWdrNawqgROQSkqlUtF7771HAMje3p5u3rypsT4pKYneeustSkpK0nosGRkZdPr0aVq7di299dZbZGtrSzKZjH766SetH5sZlpMnT1Lnzp3FB/GePXtSSEhIqfaxYMECMjExIQBka2tLa9eupZycHI1tVCoVBQYGkpOTEwF4acJSkXJzcyk4OJi8vb3FY6tfjo6O5O3tTSdPnqS8vDwiIrpy5YrYr+uff/7RamyM6cqdO3eoTp06BIAmTJhAKpVKXKdSqWjChAkEgIYOHarTuA4fPkwAqFWrVjo9LpNOQgKRkZHQDOtVrXKXLCEyNycaOFB4bdmimxgrK05AKrG8vDxxFJ3GjRvTw4cPtX5MpVJJERER9Ouvv9KMGTOoU6dO4sNg/temTZsoJSVF6/EwwxAdHU2enp4kk8kIADk5OVFgYKDGg0lpxMTEkJeXl/h9a9CgAQUU0bswPT2dfH19xYd8S0tL8vX1pczMzPJ+JFFOTg4FBQWRl5eXWLqqfjVr1ox8fX0pLCys2KZga9asIQBkZ2en9QSJMV0JDQ0la2trAkCfffaZxrqcnBxatGiRVvtpqe9V/v7+5OXlRS4uLmRkZER2dna0cOFCrR2X6Z8+fYRakD//lDqSqoUTkEouIyNDbG7Stm1belranlYvoVKpKCIiggICAmj69Onk4uJSZLJha2tL7u7u5OvrS0FBQSVuSsMqv9TUVPLx8RFHwLGwsCBfX98KG5r56NGj1L59e/G72K9fP7p06VKh7UqasJRUVlaWmHSohxhVv1q0aEG+vr4UERFRon2pVCry8PAgANS3b1+xdoQxQ3fkyBGxj5dCodDqsZKTk+nAgQO0dOlScnd3p9q1axe6V8lkMlIoFDw0fBWzbp2QgIwdK3UkVQsnIFVAcnIytWzZkgBQ//79y9yhNTExkYKCgsjHx4fc3NzEUYXyv0xNTcnNzY28vb0pICCAIiIi+IGJFaJSqSggIIDq168v3vi9vLwoNja2wo+lVCopICBAfOCQy+Xk5eVV5NC2R48epXbt2onf5/79+9Ply5dLdJyUlBQKCAggd3d3MaFSv1xcXEqVdISHh9OdO3fE/yclJZGDgwMBoBUrVpTsgzNmAP744w+SyWQkl8spMDCwQvb56NEjCgoKIl9fX3J3dxfPnYLJhouLC3l5eZFCoaCwsLBCTTVZ1RAXRySTEVlaEqWlSR1N1cEJSBURGxsrdmh98803XznyT3p6OgUHB7/0Aq5+sFJfwE+ePElpfPayVzhz5gx17dpV/A51796dzp07p/XjPn36lHx8fMjMzEysmfPz8yuUkBeXsBQcNpSI6NmzZ2LSoW7GpX64cXNzI4VCUeq5BPbt20dmZmbUpUsXjQeiY8eOkZGRERkbGxcaWIIxQ+bn5ycWYAUHB5fqvRkZGXTy5ElSKBTk6elZ7L3K3t5erIkPDg7mERiZhh49hFqQ7duljqTq4ASkCrl69ao4wdK8efPE5UqlksLCwkihUIhtYY2NjV96AQ8KCtJJnxJWeTx8+JC8vLzEfh5169algIAAnQ+De+vWLY15RZydnYsseX3y5An5+PiITUTs7OzIz8+PYmJiSKFQkJubGxkZGRWZdBScaK00UlNTydnZmQDQ7NmzNdb5+PgQAGrSpIlez2PCWGnNmTOHAGG29IsXLxa5TUnvVRYWFmJNfGBgIDf7Za+0erWQgLz1ltSRVB0yIiKwKuPgwYN4/fXXkZubi+HDhwMAQkJCkJycrLGdkZERWrVqha5du6Jr167o1q0b2rRpA2NjYynCZpWAu7s79u3bB2NjY8yYMQNLly5FjRo1JIvn8OHDmD17Nq5duwYAGDhwINauXYs2bdpobHfx4kXMmjULJ0+eBABYWFggMzNTXO/q6oqxY8di7NixaN68eYXEFhoaCjc3N+Tl5WH37t3w8PAAAGRnZ6Nnz564cOECZs1aAYViQYUcjzGpEREmTZqEzZs3w9HREadPn0Z6ejrCw8PF1+XLl5GWlqbxPrlcjpYtW6JXr15wc3ODq6srWrRowfcqVirR0UCTJoC1NZCUBJibSx1R5ccJSBX0119/YcqUKTA2NkZaWpp4AXd1dRUv4nwBZxXp559/xrx581CzZk3s3bsXLi4uUocEAMjNzcVvv/2G//3vf0hOToaxsTGmTJmCr776Cvb29hrbBgYGYvr06cjLy0Pnzp0xZcoUDB8+vNB2FWX58uX43//+h1q1auHSpUtwdHQEAERGRuLDDy/i9Onx+P134I03tHJ4xnQuMzMTgwcPxqlTp1C9enWkpKQU2sbW1lYsGFO/6tSpI0G0rLJxdQUuXAD271di2DAjqcOp9DgBqYKOHTuG1157DY6Ojti4cSO6dOkiPtwwpg0bNmzAxx9/jI8++gjffvut1OEU8uDBA3z22WfYtGkTiAgODg4IDg5G69atNbbr3LkzwsPDERYWBldXV63GpFKpMGTIEBw+fBj9+vXD4cOHYWQk3BQDA4Hx44XSurAwoEULrYbCmM4kJyejadOmMDc3R05ODnr16gVXV1fxxfcqpi3ffXcVGzasQM+eFvjll1+kDqfSk0sdANO9HTt2AADee+89jBw5ki/orMpzdHTE77//jmvXrmHo0KGwt7dHy5YtJY1JLpfjl19+gZ2dHf777z/89FOQuO6NNwAvLyAtDXjnHSAnR8JAGatAV65cQUpKCmxtbZGQkIA9e/ZgyZIl8PDw4HsV06ohQyxw8+Zf2LlzJ3L4oqp1nIBUMSqVSkxAxo4dK3E0jOmXVq1aYf/+/Th06JBY2yAlJycn/PTTz+jX7yjmzRuNixdfrNu4EWjeHAgPBxYtki5GxiqS+v705ptvwszMTOJoWFXSrFkztGvXDs+ePcN///0ndTiVHicgVczZs2eRkJCAZs2aoX379lKHw6qwGzdu4P3334dCoZA6FA0ymQwODg5FrlO3WJXJZDqLZ+zYMWjbtj8yMoBx4wB1s3hra+DPPwFTU2D1amDfPp2FxJhWKJVKLiBjklJ/79TfQ6Y9nIBUMXxxZ1Io6sE9Li4OP/74I/bv3y9VWAZj1SqgQwfg3j1g2rQXy11dgaVLASJgyhQgIUGyEBkrtzNnziAxMREtWrRAu3btpA6HVUHjxo0DAOzcuRN5eXkSR1PxlEollEql1GEA4ASkSiEibNu2DQAnIIwZEjMzobbD0lLogL5584t18+cDAwcKQ0dOngyoVJKFyVi5qAvIRo8eLXEkrKpycXFBq1atkJycjFOnTkkdToW5fv06Bg0aBE9PT9SvXx/r1q2TPMHiBKQKCQ0NRVxcHBo2bIjOnTtLHQ5jrBRatQLUrdVmzgRu3RL+LZcDAQGAvT1w6BCwZo1kITJWZkTENfRML4wZMwZA5WiGlZGRAV9fX3Tu3BmHDx/Gvn37kJCQgNmzZ6Nbt27i/FZS4ASkClGfTGPGjNFpG3bGWMWYOhV4/XVh9Kt33wXUBViOjsD33wv/XrYMePpUuhgZK4uQkBDExcWhSZMmXEDGJKVOgLdv3w6VgVYpK5VKrFu3Dk5OTli2bBlycnLg7e2NBw8eICgoCE2bNsWFCxfQp08fDBo0CNevX9d5jDzTXBXCza8YKx8pOqHnJ5MJtR1ubsCMGUD+uULHjQM+/xwYNkz4f16e5vqsLCA3F6hWTbcxM1YS3PyqdGJiYhAZGanTYzo5OaF58+Y6PaYUOnbsiKZNm+Lu3bs4d+4cevbsKXVIpXLixAnMnj0bF/9/2MTevXtDoVCgU6dOAAAPDw8MHToUGzduhK+vLw4fPowOHTrggw8+wLJly2BjY6ObQIlVCRcvXiQA5OjoSEqlUupwWBWzbt06AkDe3t7iskOHDhEAGjx4sISRlU7Hjh0JAF24cEHSOLKyil936RIRQPTee5rLly8nGjBAu3ExVhYqlYoaNWpEAOjMmTNSh6P3fvvtN5o6dSoB0Olr/PjxtHbtWqk/vk58+umnBIDmzp0rdSgldu/ePXJ3dxf/Xo0aNaKgoKCXvic5OZm8vb3JyMiIAFDNmjVJoVBQXl6e1uPlGpAqIn/zK7mcW94xZshKMj3CP/8A770H9Oih9XAYK5eLFy8iOjoaDRo0QPfu3aUOR68lJiZi1qxZSElJQdeuXVG9enWdHDcrKwvbtm3Djh07MGjQILRu3Vonx5XK2LFjsWrVKmzatAmDBg3CgAEDYGJiInVYRUpNTcVXX32F9evXIzMzE9bW1li0aBG8vb1hYWHx0vfWrFkT69atw6RJkzB79mycPHkSs2fPRkBAABQKBfr06aO1uDkBqSICAwMBcPMrxqoKHx+hmVZ4uGZTLMb0jbqAbNSoUdw/8RXmz5+PlJQUjBkzRuedpOfMmQOFQoEZM2bgxIkTlfpvlZaWBnNzc6SkpGDYsGGwsLDAgAED4OnpiVGjRuks8XsZIsLmzZvx2WefIT4+HjKZDF5eXlixYgXq1av30veGhobC1dVVLJDu1KkTTpw4gT179mDWrFm4ePEi+vbtC3d3d6xfvx6NGzeu8Pi5KLwKuH79OiIjI1G7dm307t1b6nAYYzrg7Q1kZgJr10odCWMvx/0TS+b8+fPYvHkzzM3NsXr1ap0f39fXF7Vq1cKpU6ewfft2nR9fVxISEvD2228jKysLvXr1grOzMzIzM7F3715MmjQJ9evXx5tvvolt27YhPT1dkhjPnj2Lbt26YdKkSYiPj0ePHj1w/vx5bNq06ZXJx61bt+Dm5gYXF5dC83B5eHjg2rVr8PPzQ7Vq1bB37160bt0aCxYsQGpqaoV+Bk5AqgD1hWLUqFEwMjKSOBrGDBdJ3Am9NExNgW+/FSYqvH9f6mgYK9rVq1dx+/ZtODg4cAHZSxARZs+eDSLCnDlztFIi/Sq2trb44osvAADz5s1DRkaGzmPQNiLC5MmTkZSUBA8PDxw+fBiRkZGIioqCQqGAm5sb0tLSsHXrVrzxxhuwsbFBr169sG7dOiToYCbYhw8fYuLEiXBzc0NoaCjq1q2LgIAAnD59Gl26dCnRPpKSkuDg4IBbt25hxIgRGD9+PO7nu0lYWFjAx8cHN2/ehJeXF7KysrBy5Uq0bNkSP/74Y4WNDMYJSBXAY6szqRX14G5ID/OGasgQYORIoTaEMX2kvj+NHDmS+ye+xJ9//olz587B0dERn332mWRxTJs2Da6uroiNjcU333wjWRzasmbNGhw6dAgODg745ZdfxPtTo0aNMGvWLJw6dQr379+HQqHAwIEDIZPJcPr0acyePRv16tUTk5H4+PgKjSsrKwtLliyBs7MzNm/eDDMzM/j6+uLOnTuYOHFiqe6jvXv3xp07d6BQKGBjY4PAwEA4Oztj1qxZeP78ubido6MjNm3ahPPnz6Nnz5548OAB3n//fXTr1g1nzpwp/4fSejd3JqmbN28SAKpRowbl5ORIHQ6rohQKBQGgWbNmicsOHjxIAGjIkCHSBVZKHTp0IAB08eJFqUMplnoULJVK+H9sLFG1akTDh/MoWEz/tGnThgBQcHCw1KHorbS0NKpXrx4BoN9++01cnp2dTffu3dP68RMTE+nJkyfi/0+dOkUymYwsLCwoOjpa68fXlXPnzpGJiQnJ5XI6duxYid7z6NEjCggIIHd3dzIxMdEYNczFxYV8fX3p1q1b5YorICCAGjRoIO7Xy8uLYmJiyrVPtZKOgqVSqSgwMJCcnJwIAMlkMvL09KT79++X+dicgFRyK1asIAA0efJkqUNhVRgnILpTMAEhIlq3TljGCQjTJzdu3CAAVKtWLZ0M+2moFi9eTACoa9euGsPor1q1iszMzOjbb7/V2rF37NhBNjY29MEHH2gsf+ONNwgAvf3221o7ti6lpKRQ06ZNCQDNnz+/0PqMjIxXPmw/fvxYTEbMzMyKTEauXbtW4pguX75M/fv3F/fRrl07Onr0aKk/W0lcuHCB+vTpIx6rY8eOdPz48ULbpaWlka+vL5mbmxMAsrS0JF9fX8rMzCz1MTkBqeQ6d+5MAGjPnj1Sh8KqME5AdKeoBCQvj6hTJ05AmH758ssvCQC9V3DSGiaKiooiCwsLkslkGnOkxMfHk7W1NQGgI0eOaO34t2/fJjMzM5LL5RQSEiIuj4mJIUtLS5LJZHTixAmtHV9XJk+eTACoe/fulJubW2j9n3/+KSYSPj4+dPLkyZfuLz09nYKCgsjLy0v8O6lfTZo0IW9vbzp58iSp8l+o/19iYiJ5eXmRXC4nAFS7dm0KCAjQyRxuQUFB1LhxYzFWd3f3ImvZYmJiyMvLS9yufv36FBAQUOTnKQ4nIJXY3bt3CQDZ2NhQ1stmLmNMyypLAtK+fXsCQJcuXZI6lCItWkSU7xlBbyUnE23dShQZWXjd0aNExc3zePkykb8/0aNHmsujo4l+/LHi42TapZ7Yc//+/VKHoreKq2lQPzCPGjVK6zHMmzePAFDPnj01HjA///xzsbTckCc43rJlCwGg6tWr0927d4vc5rvvviM7OzuNRKJ9+/a0bNmyV9ZqZGRkiMmIjY2Nxj4aNWokJiNZWVnk5+dHtra2BIBMTU3Jx8dHo/mbLmRkZJCfnx9Vq1aNAJCFhQX5+PhQSkpKoW2PHTsm3hcBUN++fUtcQMcJSCW2evXqSlVFygzX2rVrCQDNnj1bXHbgwAECQEOHDpUwstLR5wTkr7+Emo/q1YmePpU6mpeLiRFiPXu28LpevYg+/rjo961aJbyv4CXtn3+I5PKKj5Npz507dwgA2dnZcf/EYvz3338EgKysrCguLk5cfv78eZLL5WRmZlbsA3NFSklJIQcHBwJAf/31l7g8IyODGjZsSADo559/1noc2hAZGSnWUGzatOml2yqVSjp58iT5+PiIzbXULycnJzGReFkylpaWRoGBgTR+/PhCNSN169YV/92nTx+6UFxJjI7ExcWRl5cXyWQyAkCOjo7k7+9f6PPl5OTQ2rVrxQTNxMSErl69+sr9cwJSifXo0YMA0I4dO6QOhVVxnIBo182bRFZWwsO5IdQElCcBqVWLyMyMKH+rD05ADM/XX38tdqhlheXl5YlNPpcsWSIuV6lU4r3dx8dHZ/H89NNPYlObtLQ0cbm6aVLt2rXp2bNnOounImRnZ1OXLl3K/D2MiIggX19fatGihUYiUatWLfLy8qKgoKCXJtd5eXl08uRJ8vb2pho1alCNGjWodu3aejcgQ0hIiPidA0BdunSh06dPF9ru6dOn9NFHH5GtrS0NHz78lfvlBKSSiomJIZlMRlZWVpSRkSF1OKyK4wREe7KziVxdhQf6t96SOpqSKU8C0q4d0Zw5RG3aEKnv7ZyAGJ5u3boRAPrnn3+kDkUv/fzzzwSAGjZsqHEPVz/wOzg4FNkkRluUSqX4sO7r66uxTt15+ZNPPtFZPBVh4cKFBICaN29Oqamp5dqXOhlxcXHRSEZq1KghJiPZ2dnFvv/vv/8mADR27NhyxaEteXl59P3335O9vT0BICMjI/q4iAt1TEwMAaAGDRq8cp9VZtDt3Nxc5OXlITs7W+pQdOKff/4BEWH48OGwsLCQOhzGKgXSw7lL/vc/IDwcaNQI2LhR6mhKRj1f1+7dwI0bmq/MzFe/39cXSE4GFAqthqlzubm5SE9Px7Nnz6QORavu37+PkJAQVK9eHUOHDpU6HL3z/Plzca4PPz8/8R6enp6O+fPnAwC+/PJLVKtWTWcxyeVyrFu3DjKZDF9//TWio6PFdQqFAkZGRli/fj1u3bqls5jK4/Dhw1i5ciVMTU3x559/wtraulz7a926NZYsWYJr164hLCwMCxcuRIsWLfDkyRNs3rwZr7/+OpycnDBjxgwEBwcjLy9P4/3qSaL16d6Sn5GREWbMmIHo6Gj4+vrCxMQEVlZW5dpnlUhANm3ahKZNm2LSpElo2rQpNm3aJD5IVFY8+SDTJ0U9uLdu3Rr+/v7w5lnyymzfPuCbbwBjY2DLFsDWVuqISqZOHeHn2rVA166ar0uXXv1+Gxvg66+BZcuA2Fithqoz6vuUp6cn6tWrhwULFiAtLU3qsLRi165dYgGZmZmZ1OHona+++gpJSUno06cP3nzzTXH5qlWrEBcXhy5dumDy5Mk6j6tHjx548803kZmZCR8fH3F5x44d8e677yInJweffPKJzuMqrcTEREyYMAEqlQpffPEFXF1dK3T/rq6uWL58OW7evKkxg3pSUhL8/f3h7u5usOe2lZUVlixZgosXL2LBggXl21kF19LolZs3b9Lw4cPFqjD1yAIAaODAgRQRESF1iFrx4MEDksvlZGFhUe5qRcYqwpo1awgAzZkzR+pQXiknJ4fu3LlT5Lp27doRALp8+bKOoyrs4UOi2rWFpkxLl0odTemUtwmW2muvEY0cadhNsG7dukXu7u4aTTbU/3Z2dqagoCCpQ6xwvXr1IgC0bds2qUPRO7du3SJTU1OSy+UanZCjo6OLHI5X1+Li4sjKyooA0H///ScuT0xMFEd40udRzVQqFQ0bNowA0KBBg3Q6eldkZCStWLFCYzRItW3bthEAGjdunM7i0YYq3wQrKSkJEydOhIuLC/bv34/atWsjICAAycnJCAgIQN26dXH48GG0bdsWEydORGJiotQhV6igoCCoVCoMHjy43NWKjFUl6lLoQYMGISsrS+pwikUETJ4MJCUBr70GLFokdUTSUCiEWqD//hP+n5cHfPUV8PSplFGVjPo+1apVK+zduxe1atVCQEAAHj16hHPnzqFbt264ffs2Xn/9dXTv3h3nz5+XOuQK8eDBA5w5cwZWVlYYPny41OHonU8//RQ5OTmYPHkyOnbsKC5fuHAhMjMz8cYbb6BHjx6SxVevXj2xGdjs2bOhVCoBALVr18ai/78QzZ07F7m5uZLF+DLffvstDhw4gLp162LLli2Qy3X3GOzs7IwFCxZAUdnajuajbuVAJWllpP18SHeKG0P5aYExKdUzOapnqrSysiJfX99KM1fGoEGDSjSkHDMMmZmZZZplVJ988803hTqh65OCtaXNmzenyCImqWjbtq1e1ICsXi3UINSsSZRvdE6DUdIakKQkotjYF5MqFqwBIRI6pJuaCjUg33wj7NfWlsjPj0gfL+nZ2dnk5+cnDllZ3Fj/SqWSAgICqE6dOgSA5HI5eXl5UUJCgkSRV4wNGzbodWdbKannRrKxsaHExERx+YkTJ8RnldjYWAkjFGRkZFCjRo0IAP3www/i8uzsbHFEqDVr1kgYYdEuXrxIZmZmJJPJ6MCBA1KHo2H79u2V4ryIjY0VR0t7lUqTgOzZs4eaN2+uMXvjrVu3Xvqe27dvk6enZ6Wq7n706BEZGxuTqampwQ2JxzTdu3eP3N3dqWXLllS3bl3y8vLSuCkZiuPHj1Pjxo2pXr16VKtWLQoMDJQ6JFFJZ5zNysoiX19fsrS0JFtbW6pXr95LZ7LVpnPniExMiGQyor17dXroClOSBOSLL4j69BHm/Bg8WJjNvagEJD2dqGFDIQG5fJmof39h34Cw7dGjOvlIJbJ3716N+9Tw4cPp5s2bL31PamqqRoGZtbW1QReY9e/fnwDQn3/+KXUoeiUnJ0ccQcnPz09crlQqxQkbP//8cwkj1BQYGCg2GUxOThaX7927V5zU7+HDhxJGqCk9PZ1atWqlt02BOQExQBERETRw4EDxgt66detSj6F8+PBhatOmTaXoH7Jx40YCQMOGDZM6FFZGWVlZtGLFCnGSIhMTE40H5B9//NEgZp2Nj4+niRMnipMYqWdVBYTZe4vrZ6ELJa0tJSIKDg7WuD6o2z+rXy1btqRFixbpZGjelBSipk2Fh+sPPtD64bQmNlZIMM6fL7wuIIDo0CHhpf6ajxghJCvnzxP98kvh95w4QbRy5Yv/BwcTtW79IhEZOJCoBPNiaU3B+5SLi0up71ORkZGFCsz0KZkviVOnTpGRkRGZm5vrdAhZQ6CuGXJ2dtYYrvXXX38lQJjoTt+G1O/Xr1+RNdvqPhbTp0+XKLLCpk2bRoAwa7s+Ju+cgBiQR48e0fTp08nExETsYK5QKMo8o2pOTg75+/uLYxybmJiQt7e3QdQixMfHk0KhIDc3N5LL5WRnZ0dz586VOixWSiqVigICAqhevXoEgGQyGXl5eVFcXFyhJkItWrSgffv2SR1ykVJSUsjb25tMTU3FxMPPz49SUlI0zjFjY2OaPn06PXr0SKfxBQYGkrOz8ytrS69cuUIDBgwQt2vTpg0dPnxYY/Io9d9K/apfv36JZsMtq8mThQfq9u2JDLxVXqn06kVURIu4l8rJIVIohOZYgFBrNH06kS6/bs+ePSNvb2/xPmVjY1Pi+1ReXl6RywsWmA0YMKBEsw5LQaVSieeKesbsBg0aUN26dSk9PV3q8PRGcnKy2CRv165d4vLnz5+LTfC2bNkiXYDFiIiIIGNjYzI2Ntb4Dt6+fZvMzMxILpdTaGiohBEK1HOnVKtWjW7fvi11OEWqLAlIXFwcAaB69eq9cluDS0Dy8vJIoVCII4UYGRmRt7c3PX78uETvz8zMpClTphRb7f348WPy9vYmIyMjAkA1a9YkhUJR7M1AKomJieTv70+DBw8Wb27qh1YA5OjoSElJSVKHyUrozJkz4iRPAKhHjx4UEhJSaLvg4GCNiY4GDhxI165dkyDiwpRKJfn7+5Ojo6PYZn369OkUHx+vsZ36HDM2Nhar8BUKBeXm5mo1vqtXr5aotlRduKGOz97envz9/YuMLzc3l4KDg2nGjBlUt25djWSkSZMmtGDBZxQS8qIPQ3ls2yY8SJuZEV28WP79GYrffiN6442yv//xYyJvbyIjI+H3V6OGkJho85Kuvk/VrFlT4z6Vv6nKyyQnJ1P9+vXJ19e3yFLv3NxcvUjmi6JSqSgkJIR8fHyoWbNmGudE9erVxZGSPvroI6lD1Rsff/wxAaDBgwdrLF+wYAEBoF69eum8qWdJzZgxQ0yE85s7dy4BIDc3N0ljj4qKEmu6fymq+lRPcAKi5w4cOCC24QNAQ4YMoevXr5dqHytWrBBrOObOnVtsDcfFixfF2T3V1XbHjx+viI9RZjdv3ixypk1TU1Nyd3engIAAevToEb3++uvixcwQmupUZQ8fPiQvLy8xcaxbty4FBAS89IKdk5NDCoVCvJHrQ23dyZMnydXVVfxO9uzZ85UlXzdu3KChQ4dqNGXSRsfA4hKKgqXQubm5GoUbxsbGpSrcIHoxG676OuXqOoTw/53FvbyIgoKEmcvLIiaGyM2NaOPGsr3fEG3eTDRsmNDPo7zOnSPq1u1Fs6yuXYlCQiq+Q/fx48fFNvsAqE+fPhrDqZbETz/9JL6/efPmtLeYzj5SJfMF5eTkUFBQEHl5eYlJl/rl6Ogo1grm5eXRlStXyMLCggDQX3/9pdM49dGVK1fEWoT8hUmRkZHicLzh4eESRvhyjx8/Fv/m+fvQPn/+XCyU2bp1qySx5ebmkpubGwGgt99+W5IYSooTED119+5djXHSmzRpUubO4k+fPiUfHx+xeYitrS35+flptLnMLygoiBo3bqzRXCMqKqocn6Z0wsLCikw6rKysyNPTkwICAgqNnpKUlCSWQq9atUpnsbKSy8zMJF9fX7E/gbm5Ofn6+lJaWlqJ95GcnCx5bV10dDR5enqKCZSTkxMFBgaWqsQrKCiImjRponGO3b17t9yxFUwoXlZbeuDAAWrZsqUYw9ChQ+nGjRtlPra6FNjP7xw5O7946AWI6tYlmjFD6KdQ8DkxPZ3Izk5oapXfP/8QtW0rJC96WhBa4T74QOhc/t13RP7+ROX4c2gICiJq0oSoceNsMjGxJHd3d7p3716596seNEL9HWrcuHG5BjX577//qH379uL++vbtW2w/I10l8/llZWWJSYe6+ZD61bRpU/Lx8Sm2KaJCoRDvv7q8n+ojdTPPmTNnaiwfPXo0AaBJkyZJE1gpqP+eTZs21ehf4e/vLza7K829raJ8/vnn4rmo783pd+3KJlvbJ/TWW7r/PVWkSpOAPH/+XKMtefXq1UmhUBSbLJTGrVu3aMSIERqlTMW1qc/IyCA/Pz+xU7CFhQX5+PhoZZK//G1m8z+Uqdsvenl5UVBQ0Cvbz/73339kZGREJiYmdO7cuQqPk5VdYGCgRtMET0/PItul3rlzh06dOvXK/YWHh1Pv3r3F/XXq1IlOnDihjdBFqamp5OPjQ5aWluI54evrW+Z23dnZ2aRQKKh69epirZ63tzc9f/68TPvbv39/oYSiqNrSiizceJmoKKHZj5ubMHqVOhmxtCRydxc6XqekEKWlCctNTYkOH37x/sBAIienCg9Lr23fLnxu9asCclJRWhrRypX7yNzcXLy2+vn5lalzakpKCvn4+Iil+up9VUSHYfVQvLVr1xabNb5sNDxtJfNqqampFBAQQO7u7uK5r361aNGCfH19SzSAi0qlEmvqu3btWua+m4Zuz549YuFR/oIRpVJJ3t7eZGNjQzExMRJGWDLZ2dnUqlUrmj59usbgAkqlkjp37kwAaMaMGRQWFqaz1/fff09yuZxMTEwKNWdWKpV6953bsUO49o8ZI3Uk5WPwCYi6LbmDg8NL25JXhNK0qY+Li9NoLlOvXr1XNpcpiaI66qlfNjY2YtJR2rkgfHx8xFKJsj7IsYpz+fJleu2118S/bbt27ejIkSOFtsvfZKh+/folLjkKCgoSx2bXVm1dcR3lK2ps+vj4eJo+fbo46pejoyP5+/uXqinhlClTxHO0fv369McffxQ6R7VZuPEqt24RLV9O5OqqWTNSvbrQ5AgQ+iy0aPFiHouqmIDoQnx8fJmv6QX7PKnPBW3cp9Q19+qheF9Wc1/Ryfzz58/FpEOdZKlfLi4uJU46Cnr06JH4u1u0aFGZYjN069evJ2Nj42KH1y1qVD59Vdzzyf79+8nIyKjQd0cXL3Xz5PyioqKoW7du9OGHH+ri11JilSUBiY+PF+/dr6J3CUjBtuRubm5aH0WhtG3qz58/T927dxdj7Nq1a6lrGXJzcykoKIimT59O9evX1zhp6tSpo9FmtqQKtjPOzc0V4zSEatzKKikpSaMPQq1atYrs1KxSqej333/XSLwnT55cqIndy6Snp2vU1llaWlZYbd2ZM2eoa9eu4ve0e/fuWqtdCw0NpZ49e4rH6ty5M50+fbpE7/3qq6/EhKJgKVdxhRsPHjzQxsd4pUePhNoPd3dhlKaLF4Wb0PXrRC4uRMuWCdtxAqJd58+fp27dupWouRMR0enTpzUGjejZs2eRg0ZUtKJq7ovrH/LgwQONZN7BwaHEyXxSUhL5+/uTu7u7WEuU/56sUChKVbOiVCrpq6++oj/++ENjubqmXi6X0+H8VX5VxLx58wgAvfvuu1KHojVffPEFAUL/JFdXV5291MltwRq2M2fOkLGxMclkMr0aSXLnTuHaP3q01JGUj0EmIAXbkjds2LDUbcnLqzRt6ssyS212dnaJO+qV1tSpU0kul9PBgwc1lt+9e1csCdPHYfwqM3UfBHX7aHWn5qISioMHDxYaYKE8fRBiY2M1Snbr169f5tq64jrKa3uAA5VKRYGBgeTk5CSWMHt6etL9+/df+r709PQitzl16lShwo2wsDBthV9qCQkvmmBFRgoT6FlYEN2+zQmILhR3Tc8/muD9+/fL3eepIgQHB1Pr1q1LVHMfFhYmdsQFQK6urkU27Xz48KE4nLv6Hqg+79RJR1lrVP/55x+xprHg/D+fffaZWPukD6N46dLdu3fJ3NycZDJZpWwqHRUVRRYWFiSTyejMmTM6PXZaWprYDLfgPCWrVq0iAGRnZ0fR0dE6jas4nIBIICMjQ5xhGCh/W/KKUJo29QVnqVVXjavbEqelpVFgYGCRHfWaNWv20o56paEe3cve3r5QE4CtW7cSIMyg+6rZ4VnF2LdvH7Vo0UL8Wxc343HBjqsV3Qfh3LlzGiW73bp1K/GNriI6yleE9PR08vX1FUtiraysyNfXt8RNEvWhcKOk8icgREJn9IEDOQHRpYIDldjZ2dGqVato8eLFkp8L+ZWm5l6dzKub+KqT+fPnzxeZdBgbG9PAgQNJoVBUWB+Ed955hwBQ27ZtNc7d3NxcsbZz2LBhenleapM6AXN1da10o1aOHz+eANBbb70lyfGvXLkiJni7d+8Wl6tUKvG+27t3b72YZoETEB0LCwsr0YRgUilNm/rLly9T//79xW1bt25Nr7/+eqGaDvVEZcePH6/QL71SqaRBgwaJzQcK7tvLy0tsyqKLdu5V1fXr1zXmmmjatGmRCUVRHVe11QehLLV1AQEB1KBBA/FzFNdRXpdiYmLE7zEgjKwSEBBQ7Pb6WLjxKgUTkIcPiWxsiN57jxMQXbtw4YJGQZT6Vd4R0ipaaWrunz17RnPnzhXnjspfKCaTyah79+60atWqChkVrKDU1FTxfl+wRPr27dtUrVo1ksnktGnTlQo/tj5LTU0Vmwu97HpmaI4fPy4WGFVUH8GyWLNmTZGFs4mJiWJT3KVLl0oWnxonIDq2ePFisrCwoA4dOtB///0nZSjFKm2benXVuKmpqViC5urqWuaOeqWRmJgojrv91Vdfaax79uwZNWrUiLp0WUZffVWFplDWEXWnZvWNvbg+COpO3LrouFpQwdo6a2tr8vX11Rj559KlS9SvXz/xoaRdu3Z09OhRrcdWGkePHqV27dqJMfbv358uX74srlf/jtV9qyq6o7w2FUxAiIjWryeSyzkBkYJKpaJffvlFrGEorr+FPiiYMHXs2LHYmvubN2+K54eLiwutX79eJ+dHSEgImZqaFiqRJiL688+d5Or6iCwtifRkblWdCQgIIEDo/1kZBoxRKpXiXDi+vr6SxqJSqcjDw6PIwtlDhw6RXC4nuVwu+X1u1y7h2j9qlKRhlNuDBw/EPmevIlkColKpxIcDfRsOrSilaVOfmZkpPojqul3ngQMHSCaTkbGxcaFOu+fPJ5GJiTAjsMRzKlYa6hmP1TMSGxkZFdupuaSznWtbZGQkeXp6inE4OzvTpk2byMvLS+ywWrt2bZ308yir4oYnPXLkCPXt21f8bO3bt9fbwo2iFJWAKJVE3btzAiKV58+fi4UK+k6pVNKvv/4qFkTJ5XKaOnVqkfcp9TUgMDBQpzEuX76cAGEwjoIFL++9J3z/XVwqZvJJQ6FSqcQBYxYuXCh1OOWmTtobNmxYIcNRl1dSUpJY2/HFF19orJs1axbZ2dUgd/eDJOVUIZyA6FD+BMSQHDp0SGPY3h07dhS5nbqUuSzjypfXnDlzxCYqBSdd+/pr4Utep47Q6ZWVz9WrV8WHdldXVzp58mShbcoy27ku7NixQ5wzQF1bZ2RkRO+//75Gx1t9lpiYSNOmTRP/Buombba2trR27VqDKNzIT6kkCgt7MfyuWkIC0ZWq1TJFbxhSAqKWf+6qjz/+uMhtpEpA8vLyxJoaD4+xGhNrpqURtWwp3KPef1+nYUkuLCyM5HI5mZqaUmT+EggD8+zZM7Fg6M8//5Q6HNGxY8dILpeTsbGxxiAMWVlZNGhQLAFEEnVVISJOQHTKUBMQIqHT3Lp166hXr16FhlJVkzIByczMpLZt2xIAmj9fczZcpZJo0CDhiz5sWNWZVVlb1M0YVq1aVai2ICsrS6MTt5mZmeQdVwvKzMykKVOmiLUeLxt6VJ9duHCBBgwYQMuWLaMZM2YY7Gg6R44QNWumWQNCRBQSQtS8uTQxVXWGmICo3b59m5KTk4tcJ1UCQiT052rT5i1q3DiTVqzQXHflCpG5uXCP+vtvnYcmqUmTJhEAGm3AHQHmz58vdu6WupCtIHVs9evX1zgv7t4V5mECiH74QZrYOAHRIUNOQEpCygSEiOjatWvUr989AogKjr6bmEhUt67wZV+3TpLwKg11AlKw/XReXp5Gc6uBAwdqvQ9QWf37778EgAYNGiR1KFXezZvCeVlwouv//hNmUGe6V1wCEhoaSv7+/hQeHi5RZOUjZQJC9GLiNTMzooK/wjVrhHW2tkR6MkqqTiQkJIjD5hccUt8Q3Lp1i0xNTUkulxeal0wfZGVlUceOHalevYb08cea9+OtW4XvnLk5Ub4uhTrz779ETZoYfs1faRIQOfRMSEgIBg0ahP/9739Sh1IuRAQAkMlkkhzfxcUFkyY1BgC8/z5w69aLdbVrA7//DshkwPz5wMWLkoRYqRkZGcHd3R1t2rTB4cOHERwcjNatW0sdFmOsguzZswfvv/8+9u7dK3UoZSL1PWrMGMDbG8jOBjw9gefPX6ybPRtwdweePQO8vAClUpIQda5OnTpYuHAhAGDu3LnIzc2VOKLSmT9/PnJycjBp0iR07NhR6nAKMTMzw99/7wTRRXz7bWv89tuLdW+8AUycCGRlAW+/DWRm6iamnTuBPXuAQYOAu3eBH34Qvu8//ggkJ+smhoqkvp6ory8vI1kCUtzFLzk5GYcPH8alS5ckiKpymTwZeOcdID1dOKFycl6sGzJEuMhnZwPjxwOpqVJFWXnNnz8fYWFhGDBggNShMMbKSOoH9crs66+Bjh2Be/eAadNeLJfJgF9/BRwcgJMngRUrpItR1+bOnQtnZ2dcv34d/v7+UodTYocOHcLu3btRvXp1rNDjP1jz5o2wdq0dAGDmTODKlRfrNm4EWrQArl0D5s3TTTxffgmMGgWcPftiWW6uUHAcE6ObGKSidzUgrGL98APQvDlw4QKwaJHmuhUrgE6dgNu3gblzpYnP0L3s4cTc3BxmZma6DumlDhw4gDfeeAO/5S/6YXpDXWhUvz5gZvbiNWiQtHExpg1mZsDWrUC1asC2bUBAwIt1tWq9qKlfuhQ4c0ayMHXK1NQUX3/9NQDg888/R7IBFIPn5eVh7v8/RCxcuBB16tSROKKXK662w8pK+M6ZmAB//AEkJuomnubNgY8/Npyavvv372P8+PHlrv3lBKSSs7YWTiRTU+Cbb4Dg4Bfr8l/8f/4Z+Osv6eJkunH79m1s27ZNo4aRS3j1x927ws8TJ4CkpBevf/6RNCxWBD5vKoazM7BmjfDvDz/UbC48eDAwZw6Qlwe8+Sbw9Kk0MeraqFGjMGTIEDx9+hRLliyROpxX8vf3x/Xr1+Hs7CwmIvouf23Hp5++WN69O/DLL0B4OKCrPOrDD4EnT4ANG3RzvLLKyMjA559/jlatWiEwMBCLFy8uUVOr4nACoiX6dHPq3BlYtgxQqYT2tPmz+mbNgHXrhH/PnAlER0sSImMMwvkIAE2aADY2L15WVtLGxZg2TZ0qNENJTwdWr9Zc99VXQLt2gLGx7kqk9cGaNWtgYmKCH374AVfytxPSM48fP8bixYsBACtXroSpqanEEZWMldWLwtnvvtMs5PHyenEt1gVzc6GAePFiIC5Od8ctKaVSiXXr1sHJyQlffPEFcnJy4O3tjcOHDxd6xo0uxUOk3vUB0acH98rkk0+Avn2FC/jOnZrr3n1X6BDYurVQ3c0YY6xy07d77S+/CA9gBUuBzc2BXbuEZsR79wLbtxf93vxt6CsDFxcXzJgxA0qlEnPmzJE6nGItW7YMT58+xeDBgzF69GipwykVV1chwQWAKVN02+ciIQFISXnx/9GjATc3ocZPn5w4cQJdunTB7Nmz8fjxY/Tu3RshISFYt24datasKW6XnJyM999/H71794aZmVmJmnlzDUgVYWQEbNkC7NgBfPBB4fUBAcDx48JPdY1Ifrt2CaMyMMZYVaJvD+qVVY0aQk19Ud3mmjQBbG2FfoxvvikkI/kpFMCxY7qIUreWLl0Ke3t7HD16FLt27ZI6nEIiIiKwceNGGBsbY+3atVKHUybz5gkjrj19CkyYoN1+GDdvAkuWCIW9Dg5Cf5P81q8XkuxTp14sy8rSXjwvEx0dDQ8PD/Tt2xcXL15Eo0aNEBQUhBMnTqBTp07idkSETZs2oV27dvjxxx+hUqnw5ptvamxTHE5AqpD69YWajqJYWwtJyrp1wuhYBduc//sv9xEpCj+csIpUowZw6BBgaam5vF07YTljVV3DhkKHXZVK6ki0z87OTuwD8sknnyBLqqfRYsyePRt5eXmYNm0aXFxcpA6nTGQyoQZNPeLa8uUVu/+oKKFZYffugIuLMKDC9etC39uCo486OwutVdR9Up4/Bxo0EDrMJyVVbFzFSU1NxYIFC+Di4oK9e/fC2toafn5+uH79Ojw8PDS2PXToEFq3bo1Jkybh4cOHGDJkCCIiIvD777+jdu3arzwWJyCskFathGrAjAypI2EVraiEiZMo/VGrltDx1tpac7mdHY+EpW/4vJHGggVCx+FffpE6Et2YMWMG2rVrh3v37ulVLcPevXtx5MgR1KxZE19++aXU4ZRL7dqAv7+QjFy8+GI0wrIgEmowZs0CmjYVau8+/RQ4fx6oW1eY++bkSaHGpajp7hYtevHsFRoqdE7fvBlo2xb46Sft1dCoazJatWqFlStXIisrC15eXrh58yZ8fHxgYWEhbhsVFQUPDw8MHToUN27cQOPGjREUFISDBw+iVatWJT6msTY+CDPsm9OHHwodopYuBVaulDoaxhhjTFCnjtCMZcECYORI4eGxMjMyMsLatWsxYMAALF++HJMmTYKjo2Oh7ZRKJY5VcDu0zp07w9bWttDy7OxszJ49GwCwaNEi1KhRo0KPKwUPD+DIEcDRUXjQHzsWyNfFASdPCqOx9e9f+L1EwOnTwlDSQUGag/nUqyfsy9MT6NkTkBco9q9XD6he/cX/zcyAVauEQuARI4Tv+Ny5wP79wPTpQm3KN98IzcYqytmzZzFr1iyEhoYCAHr06IF169ahS5cuGtulpqbiq6++wvr165GZmYlq1arhiy++wAcffFCmwQckS0C4E7r+Uo/IMH680CaybVupI2KMMVaRDPle+9FHwkSFn36qOXdIZfXaa69hzJgx2LlzJxYuXIiAIj50RkYGBlVwNenx48fRp0+fQss3bNiAu3fvonXr1vjoo48q9JhS6t9f6Jfx/vvCUOhbtrxY9+efQs2EOgFRqYS5abZtE/rIxsa+2NbZWZhfxMNDmGizYNKR3549hZe9/rrwUtu3Dzh8WGgef+2asN+BA4W+T61bl/3zJiQkYP78+diyZQuICHXr1sXKlSvh5eVVqJXE5s2bsXDhQjx48AAymQxeXl7w8/MrMhkuKa4BYUUaPVr4gs+YIWT+jDFWFRnyg3plZWwsNMHq0UMYwrcqWL16Nfbv34/Nmzdj5syZ6Natm8Z6Y2NjDBw4sEKPaWdnV2hZQkICli5dCgBQKBQwNq58j5H16gmJwfHjwuih+Z0/Lwzms3270L9DTd3HduxYoFevlycdZTFwoNA8bONGwNdXSEg6dhRGMf3qK8DevuT7ysrKgp+fH1avXo309HSYm5vDx8cHn376KawKjPle0tqRsqh83xxWYX74Qeg09fPPUkeiv/jhhLGqic99aXXpIpQyL1hQvjb7hqJx48aYO3culi9fjg8//BAhISGQ53vKtbCwQHD+mYa1xNfXF6mpqRg+fHiFJzz6okYNoZZtxgzg8mVhrhC1t98G7t0T/u3iIjSt8vQsX01ESZmYCH1Lxo0DFi4Uamh+/FFIlvz88jBhglzjO1GcgwcPiknk4MGDsXbt2kKDCJS0dqQ8uBO6llSGm5OTkzAig69v4dEamGHiTuiMscri66+FJin5Z0+vzBYuXIh69eohPDwcmzdv1vnxw8PD8fPPP8PMzAzrihqvvxKZM0fo8F2wH+wHHwgP/2FhwndPPayuLtWrB2zaJBx/6FDg4UPg11//RevWrXGoBMMljhw5EhMnTkRQUBAOHTqkkXxkZ2djyZIlaNasGTZv3gxTU1P4+vrizp07mDhxYoU+J3AfEPZSCxYIWfb27UJ1N2OMMaYP6tYFVqwAZs6UOhLdsLa2xooVKzBx4kTMnz8f/fv3R7Vq1XRybCLCxx9/DJVKhZkzZ6KZLqcKl4CZmTAp5siRQq2H2iefSBdTQa1aCZ3T//47D4sXz8fduzcxbNgwTJgw4aX9M2QyWZH9iLZt24bPPvsMd+7cAQB4enpi+fLlWvtbcxMs9lLm5kJTLB4ClDHGKg9DKuy7e1do4/7dd0D79prrpk0TmsiUYN6zSuGdd97B+vXrcePGDTRs2FCnx7axsUGNGjWwaNEinR5XKoMHCx2+FywoXR8LXZLJgLfeMsbYsRfw/fffY/Hixdi8eTO2b98Ob29v/O9//3tlkhoREYHZs2fjyJEjAIA2bdpAoVBgwIABWo2dExAGAEhLEyY/mzVLmPQsv4EDgbVrC0+OxhhjlR3X1ksrM1Mohc7MBA4cECYizM/YGHjvPWlik4JcLscff/yBI0eO4H9FTSShRV5eXpg8eXKlGHa3pNauFWoaGjbU7yTX1NQUs2bNgqenJxYsWIAtW7Zg5cqV2LJlC5YvX15k343k5GT873//w6+//oq8vDzY29vjq6++wpQpU3QyuAAnIAxEwDvvCEPMbd4sVGsX9P9DfrMC+CGEMca058MPhbbu7doJ/RIZ0Lx5czRv3hwffPCB1KFUevXqCaNMeXvrdwKi5ujoiE2bNuHjjz/GrFmzcPbsWUyaNAnfffcd1q1bh+7du0OpVGLDhg1YtmwZnjx5AiMjI3h7e8PX11enySUnIFpiSA+ma9cKk+fUrClMtMMqL+6EzhgzFH/8Afz2G2BhIczDYG4udUSssqtZUxjdKr8PPgD+/Rdo1EiSkMqkS5cuOHHiBPz9/fH5558jJCQEvXv3xtixY3H37l2EhYUBALp3745169aha9euOo9RslGwuFpbP4SECO0bAWH2z/r1pY2HMcYYi4wUhkEFhNmfdT3SEKuaPDyAv//WXGZsLAx1+/8j1xoMY2NjfPjhh4iJiYGvry+MjIwQHByMsLAwNG7cGEFBQTh79qwkyQfANSBVWmqq0PQqN1e40I8eLXVEjDHGdEGfC/tycoR7U1qacF+qKqNcMaYNVlZWWLJkCcaOHYu4uDhER0fDy8sL1tbWksbFCUgVNmsWcOeO0LZ27VqpozFM+nwTZ4yVX3Hn+MyZMzFmzBjUq1dPirAqtc8/F+ZZaNgQ+PVXqaNhrHJo27Yt2rZtK3UYIr1LQJydneHj44NWrVpJHUq56PuDKbetrZoGDhwIf39/jYuQvn9XGdNH9erV0/vkIzw8HA0aNEDt2rWlDqXEgoOBVasAIyNhUBRbW6kjYoxpg95NRNiqVSv4+flJEVKVERUljCwCcNtabbp//z7S09M1ZhmVmr6VgDDGKl50dDQ+/vhj7N27F5MmTcLvv/8udUglkpAATJgAqFSAry/Qu7fUETHGtEWyTuiGLDMzE0uWLEGXLl2Qm5srdTilkpsLjB8PPH8OjBrFbWu16aOPPkLbtm0xceJEJCQkSB0OY6ySy8nJwapVq9C+fXvs3bsXpqamqFu3rljgp8+IgMmTgaQkoH9/oRkWY6zy4gSklPbs2YMOHTpg6dKlCAsLw6FDh6QOqVSWLAFCQ4Xx1Lltrfbk5eXByckJcrkcmzdvRuvWrfHdd98hj8c5ZsygGBsbo2PHjsjOzsaSJUuQnZ0tdUiFEBE2bdqEpk2bYv78+UhNTYWXlxfu3r0LPz+/Qi0Nzp07h5CQEFSrVg0ffPABZs2ahVOnTkmaqKxdCxw6JAyDumkTIOenE8YqN5JIVlYWASATExPKzMyUKowSO3v2LHXr1o0AEADq0KEDHTt2rMht09PTSSaTEQC6ePGiTuN8meBgIrmcyMiI6MQJqaOpHN577z2ytbWloUOH0t27dwutj4yMJE9PT/F706xZMwoMDJQg0qKdOnWKWrVqRY0bN6bly5dLHQ5jemnLli3iOdy2bVs6cuSI1CGJzp07p3Fv6tatG507d67IbWNjY2nChAni/cnW1lZ8HwBq06YN+fr60tWrV3X8GYhMTIgAop07dXpoxphEJEtAlEolKRQKMjY2JkdHR/L39yelUilVOMV6+PAheXl5kVwuJwBUp04dCggIKDLWvLw8UigUZG9vT8bGxmRkZEQAqH79+jR9+nQKCgqi3NxcCT4FUUICUZ06wgV+8WJJQqiULly4QObm5gSArKys6MsvvywyoT58+DC1bdtWvNG/9tprdOXKFQkiFiQmJtK0adPE76iDgwOlpKRIFg9j+i44OJjatGkjnsMDBw7U+YN6fgkJCSW+Nz1//py8vb3J1NSUAFC1atVIoVBQRkYGnTx5kry9vcnBwUEjGXFyciJvb286efIkqVQqrX2O58+JmjYV7k0ffKC1wzDG9IxkCQgR0cWLF6l9+/biBa9v3756U2OQmZlJvr6+ZGVlRQDIzMyMfH19KTU1tcjtw8LCyM3NTaOU7J133qH69etrXNQbNGhAs2bNouPHj1NeXp5OPotKRTRsmHCBd3MjkigHqrTi4+Np+vTp4oNAcQl1bm4u+fv7U61atQgAGRsb0/Tp0ykpKUlnsaq/15aWlgSAzM3NydfXl9LT03UWA2OGKicnh/z9/cne3l7jHH706JHOYsjKyiJfX1+ytrZ+5b0pLy+P/P39xeRCLpfT9OnTKT4+vtC22dnZdODAAZo6dap4jcpfc+vj40OhoaEV/nkmThTuTe3bExlAYwjGWAWRNAEhIlKpVBQYGEhOTk4EgGQyGXl6etL9+/cliykgIECMBwB5enrS7du3i9w2KiqK3N3dxW0bNWpEQUFBGttERESQr68vtWrVSuOibmtrS56enhQQEEBpaWla+zzr1wsX+Bo1iGJitHaYKi80NJR69uwp/n07d+5Mp0+fLrTdkydPyNvbm4yNjQkA2dnZkZ+fH2VnZ2stNpVKRQEBAWJCLJPJyMvLi2JjY7V2TMYqq8ePH5O3t7dYg1ijRg1SKBRaL1QKDAwkZ2dnjXtTZGRkkdueOHGCOnXqJG7bq1cvCgsLK/Gx1PetFi1aiPto2tSF7O2JvLyIgoKIcnLK93nS04m6diUyNye6dKl8+2KMGRbJExC19PR08vX1FZuzWFpakq+vr077h1y6dIn69+8vXmzbtWtXbFvf1NRU8vHxIQsLCwJA1tbW5Ofn98qS5JCQEJo/fz41bdpUIxmpUaMGTZ48mfbsyaGsrIr7TNeuEVlaCgnI5s0Vt19WtNIk1Ddv3qThw4eL34EWLVrQvn37KjymS5cuUb9+/cTjtG/fvtj+S4yxkrt+/ToNGTJEPLdatWpFBw8erPDjXL16lQYMGKDRV+Pw4cNFbnvv3r1ChWKBgYFlbkalUqno7Nmz9Omnn9KYMd+SMF6V8Kpfn2jWLKKTJ4mKa0H9/DmRjw/Rrl2F1y1aRBQWJryfMVa16E0CohYTE0NeXl4aTZYCAgK02gY1MTFRoy1t7dq1i21Lqy5JrlevnkZJclxcXKmPe+/ePVIoFOTm5kYymYxcXLoQIJQGubsTBQQIF+/ySEkheucdoqlTy7cfVjppaWkaCbWVlVWxCXVwcLBG7djAgQPp2rVr5Y6hNN9rxljZBQUFUZMmTcRz2N3dvchBKUrr0aNHNH36dLG21N7envz9/YvsS5iSklKmQrHSunePSKEQmvPKZC+SETu7FzUj+Stz4+KE9ZaWwnvzs7Qk+vvvCg2PMWYg9C4BUTt27JhG/5B+/frRpQquo83KyiI/Pz9xJBBTU1Py8fGhJ0+eFLn9mTNnqGvXrmJM3bt3p/Pnz1dILFeuXKG1a09QmzakUcJkbU00fjxRYKBQXV2UTz8latKEKDxcc/mxY8Ly1FSq0FoVVnLFJdQF5eTkkEKhIBsbG3F0OG9vb3r27Fmpj1nwe21mZkY+Pj709OnTCvhEjLGiZGdnk0KhoGrVqon3E29vb3pehlIk9YAmNWrUIABkZGRE3t7e9Pjx40LbKpVK8vf3J0dHx3IXipVWVFTRyYitLZGnp1CIFhkpLGvZkmj0aM33cwLCWNWltwkIkXBhDQgIoNq1a4sd6Ly8vCgxMbHc+y7Yltbd3Z1u3bpV5LbqkbDUQxfWrVtXq7UyMTHCRX3gQCJj4xcXdSMj4UKvUBDl/xVMnixc/Lt0IcrfBHnfPuF95a1FYeV39OhRateunfh969+/P12+fLnQdjExMfT222+L37X69etTcnJyiY+zZ88eat68eYm+14yxilfSQSleZs2aNSUq6Dp16hR17txZ3LZnz54UEhJSUR+lVK5fJ/riC6Ezef5CtD/+EH7u2EFkairUkKhxAsJY1aXXCYja06dPycfHh8zMzMTO235+fpRVhmL9K1euFGpLGxwcXOS2BUfCUo8YpM0O4wVFRRGtWSMkHnL5i4u6uTnRyJFE27YJCUifPkIV+Pr1L97LCYh+KU1CHR4eTn369KE333yzRPuOiIiggQMHit/r1q1bF/u9ZoxpX1GDUpw6dapE73327Bl17NiRNm/eXGRBV3R0NHl6eooFFU5OTuXq51HRbt8mWrGCaMAAouho4T4UEUE0fz5Rw4ZE6lsoJyCMVV0GkYCo3bp1S6NzXfPmzWnPnj2l2seXX35J6nHQly9fXmwn98DAQGrWrFmJRsLSlSdPhCptT08iKyvhou7lJSQgI0cSffutcEGPihK25wREPz158oR8fHzEMfmLGwVLqVS+svnGs2fPyNvbm0xMTAgA2djYkEKhoJzyDk/DGCu3ih7lMSMjw+CG0Vb3AYmIEJoDN2ggJCJEnIAwVpUZVAKiFhwcTK1bt9botBsREVGi92ZmZtK8efPowYMHRa6/fPlyiUfCklJyMtEvvxCdOfMiAcnJIXJxIRo3TtiGExD9Vp5RsNRtxGvWrKnRRrw0zbUYY7pR3lEeixtGO8YAxlXPn4AQER04IDQtjojgBISxqswgExCiF5121R1tTUxMyjUhVMERg2rVqlXsaCP6Rp2AEAkdz2Uy4SLPCYhhCA4OJhcXlxKPgnX8+HHq2LGjuH2fPn3owoULOoyYMVYWZRnlseDgJx06dKD//vtPh1GXT8EEhIho+HCisWM5AWGsKjPYBEStvBNCZWdnk5+fH9nZ2ZVoJCx9lD8BISJ6912ipk2Fcdc5ATEMxY2Clb8JVsHx/Rs3blxo0kvGmP4rySiPxQ1+YmjDaBeVgNy5Q2RhIRSWcQLCWNVk8AmIWlkmhNq7d6/GLK/Dhw+nmzdv6ijiilMwAXn4kMjGRphLhBMQw5KcnKyRUNvb29PKlStp/vz54vj+1apVIz8/P8rIyJA6XMZYGRU3KEVMTIzkg59UpKISECKi5cuF5ZyAMFY1VZoERK0kE0Jdu3ZNY8QgFxcXgx4xqGACQkT0/fcvRsziBMTwnDx5klxdXcX23urv6pgxYypkgjPGmH5ISkqi999/Xyx0sLa2Fs/3wYMHV8ikpFLKyiIKDn4x8lXB5QkJ0sTFGJOWjIgIlUxOTg6+//57LF68GKmpqTA1NcWMGTPw6aefYtWqVfj++++Rm5sLGxsbLF26FDNnzoSJiYnUYZfZu+8CT58C//zzYplKBbi5AefOAc+fA9WrSxYeKyOVSoXffvsNN2/eRGxsLKZPn47XXntN6rAYY1oQGRmJefPmwdnZGceOHcOaNWvQv39/qcMqt927gd9+A1atApydNde9/z4wdCgwerQ0sTHGpFMpExC1qKgozJ8/H9u3bwcAVK9eHSkpKZDJZHj77bfh5+eH+vXrSxxl+V27BuTlAe3bay6PiQEiI4F+/QBjY0lCY4wxVgoqlQpEBCMjI6lDqRBnzgiFYXFxQL16musaNwY++ACYP1+a2Bhj0qnUj6WNGzfGtm3bEBYWhlmzZqFHjx4IDQ2FQqFAx44dpQ6vwvzxB1CrVuEE5OZNYP16wMQE6NtXmtgYY4yVnFwulzoExhjTukqdgKh17twZJ0+ehEwmAxFVugv82bOAk1Ph5fb2wL59wIoVuo+JMcYYY4yxolSJBAR4Uaokk8kkjoQxxhirGuLihJ+ffw44Omque/ZM5+EwxvRElUlAGGOMMaZb6vFdnj4t3BdRqdR9PIwx/cAJSCVx+zbw44+ay06dEn7GxwNt2+o+JsYYY1VbnTrCz2+/LdwJ/d9/dR8PY0w/cAJSSURFAX//rbns4UPhp62tzsNhjDHGGGOsSJyAVBKDBwMBAZrLLlwAXF0BKytpYmKMMcYYY6ygyjUcFGOMMcYYY0yvcQ0IY4wxxrSidWvg7l1hWPiCTpwAqlfXfUyMMelxAsIYY4wxrbCxEV5FadBAt7EwxvQHJyCVwOzZgLV14eUWFsCGDYC5uc5DYowxxhhjrEgyIiKpg2CMMcYYY4xVDdwJnTHGGGOMMaYznIAwxhhjjDHGdIYTEMYYY4wxxpjOcALCGGOMMcYY0xlOQBhjjDHGGGM6wwkIY4wxxhhjTGc4AWGMMcYYY4zpDCcgjDHGGGOMMZ3hBIQxxhhjjDGmM5yAMMYYY4wxxnSGExDGGGOMMcaYznACwhhjjDHGGNMZTkAYY4wxxhhjOsMJCGOMMcYYY0xnjKUOwJDcvXsXUVFRAABnZ2c0bNiw0DbPnj1DfHw8WrRoAWNjzV/vw4cPERkZCRMTE7Rs2RI1atTQWB8ZGYnU1FS4uLjAwsKi0L7v37+P5ORkNGnSBHZ2dhrrMjIycPfuXTg7O8Pc3Ly8H5WxUlEqlbh8+TKePHkCuVyOzp07o3r16oW2i4mJQV5eHpo0aVJoXWRkJGJjY1GjRg20bNlS4xwgIly4cAHGxsZo3759kTFcvXoVOTk56NChA4yMjMTlT58+xbVr10BEaNasGRwcHCrgEzNWciqVCleuXEFUVBTMzc3RuXNn1KpVq9B2169fh0wmQ6tWrQqtu3HjBm7duoXatWujQ4cOsLS01Fh/+PBhGBkZoX///kXGcObMGWRkZKB3794wMzMTlycmJiI8PBxEhHbt2qFBgwbl/LSMlY4+3D8iIiKQnZ2tcf+IiorCkydPCm1b3DMaKyVir3Tjxg3q1q0bARBfMpmMRo0aRY8fPyYioitXrpCXlxeZmpoSAIqOjhbfn5GRQdOmTSNjY2Px/VZWVrRx40aN47i4uBAA8vf3LzKOtm3bEgD6/vvvxWU3b96k6dOnU7Vq1QgAnTx5Ugu/AcaKt3PnTqpXr57G+WFubk6LFi0ilUpFeXl5FBgYSG5ubgSAnJ2dNd5/79496t27t8b7GzduTKdOnRK3iY+PF9dFRkYWiiEuLo7kcjkBoNu3bxMRUV5eHn3++edkaWkpvtfY2JgWLFhAKpVKu78Uxv7fnj17yMnJSeP7LZfLac6cOaRUKkmlUlFQUJB4fjRo0EDj/VFRUdS5c2eN9zs4OND+/fvFbR49eiSuu3jxYqEY4uPjxfPj8uXLRESUm5tLc+bMEe9Z6vNjzpw5Wv19MJbfzp07qX79+lq9fzx48KDU9w8iohEjRmjsV/2Kj4/X3i+kCuEmWK/w8OFD9O3bF1FRUfjrr7+QnZ2NvLw8bN++HeHh4Th27Bj27NmDrl27IikpCa+//nqhfSxcuBCBgYHYtm0blEolnjx5gkGDBmHmzJk4e/asuN39+/chl8uxbdu2Qvu4efMmrl69CrlcjujoaABAeHg4OnTogJs3b2Ls2LFa+x0wVpwDBw5g3LhxqFevHk6fPg0iwvPnz7F06VL89ttviIqKwieffIL3338fTZo0gaurq8b7lUolRo0ahWfPnuHq1asgIty4cQNyuRyjRo1CZmYmAIjfeblcjh07dhSKY+fOneK/1dv6+/tj+fLlWLNmDbKzs5GRkYGZM2fCz88PmzZt0s4vhLF8jh49ilGjRqFOnToICwsDESE9PR0bNmzA33//jdu3b+Ozzz7DxIkT0ahRI3Tr1k3j/USEsWPHIj09HRcvXgQRITIyEtWrV8fbb7+N1NRUAJrnx/bt2wvF8c8//4j/vn//PgBg9erV2LRpE3bs2AGlUomcnBx88sknWLt2rcb2jGnLwYMHMW7cODg6OmrcP5YtW1bq+8eVK1fKfP/YtWuX+G/1toDw/Dd+/Hg8efJE48W16BVEwuTHIMyaNYtkMplGNq2WlZVFREINx7Nnz4iI6JdffilUAxIdHU0nTpzQeO/t27cJAC1fvpyIiBISEggAjRkzhuRyOcXFxWlsv3jxYmrYsCF17tyZ3njjDSISSniTkpKIiOjff//lGhCmcy4uLuTo6EhPnz4ttE59fiQnJ1NOTg4REb311luFSrBOnz5Nd+7c0Vi2ceNGAkChoaFERLRlyxbx/GjTpk2hY3Xv3p3GjBlDAOjHH38kIqFU+MCBA4ViMjc3p8mTJ5ftAzNWCp06daI6deqI94f81OfH48ePKTs7m4iI3n333UI1IJcvX6abN29qLPv11181rveBgYHi+eHk5FSohs/NzY1Gjx5NcrmcFAoFERGlpaUVqi1R16QsXLiw7B+asRJq3bo1OTo60pMnTwqtK+r+8fbbbxe6f5w5c6bQ/eP7778nABQSEkJEr75/9OjRQ7x/5G+BUrduXfLx8Snfh2TF4hqQV9i9ezdcXV3h5uZWaJ26Ha2FhQVsbGyK3UfDhg3Ru3dvjWX29vYAgMePHwN4kXV/8MEHsLa21qgFISL8+eefGDNmDJo0aSL2QzEyMiqyHTFjunD79m1cv34dEydOhK2tbaH16vOjZs2aMDExKXY/PXv2RNOmTTWWqb/X6vMjKioKFhYW+OijjxAREYGIiAhx26ioKJw/fx5TpkxBjRo1xPPD3t4eQ4cOLRRT9erVxf0ypi2xsbG4cOEC3nrrrSLvD+rzo0aNGjA1NS12P+3atUOLFi00lhV1fhgbG2POnDmIiYnB+fPnxW3v37+PM2fO4J133kHdunXFe42VlRU6dOigsV/1+9q0aVO6D8tYKd25cwfXrl2Dl5dXoT6tQMnvHz169Ch0/yjq+crc3LzY+8e5c+fE+4f6/MjLy0NSUhLXdmgRJyAvkZWVhejoaDRv3rzC9x0SEgIAYoco9UNTy5YtMXz4cPz999/ituHh4bh79y7GjRuHRo0aidsyJqVbt24BEAZkqGghISGQy+Vo27YtAOEG0rBhQ/Tt2xcODg4a58fWrVtha2uLwYMHo1GjRhpV6AVFR0cjKSmp0IMXYxVN2+eHTCYTv8fR0dGoX78+3Nzc0LBhQ2zdulXcdvv27bCyssKwYcOKvH/ExMTA398fn376KSZMmIB3330X48ePr/CYGctPfX5o6/kq//0jKiqqxPcP9fmRlJQElUqFv/76Cy1btkStWrXw2muv4d9//63weKsqTkBeIiMjA8CLbLqiKJVKLF26FM7OzuKFPioqCmZmZnB0dMS4ceNw/vx53LlzB4BwgjRo0AA9evRA48aNkZycLLb9ZUwq6enpAFDhtXBxcXH46aef8O6778LR0RGAcH40atRIbNv7119/gYgACOeHh4cHTExM0Lhx45cm6L6+vrCxscFHH31UoTEzVpD6/Kjo+0d8fDw2btyIt99+WxyJMTo6Go0aNYJMJsOYMWOwdetWKJVKAML5MXz4cFhaWqJx48aFEvSYmBh8++23+OuvvyCTydC+fXvI5fxowLRLF/ePevXqARDuH40bNy7V/aN69eqYO3cuxo0bh7Vr12Ljxo1IT0/HsGHDcPjw4QqNuariq8xLWFtbQyaTISEhoUL3++mnnyIsLAz+/v5i1Xt0dDScnJwgl8sxbNgwWFlZITAwECqVClu3bsXo0aMhk8nQuHFjcXvGpFStWjUAqNDzIy0tDW+88QZsbGywfPlycbn6BgIAY8eOxb179xAaGorr16/j0qVLGDduHAC8NAH54YcfsHnzZqxZswa1a9eusJgZK4p6GNGHDx9W2D5TU1MxZswYmJub4+uvvxaXR0dHa5wfDx8+xKlTp3Dnzh2Ehoa+9Pzo1asXIiIiEBcXh/Xr12POnDnYsGFDhcXMWFHU94+KPD+Ku38UPD8K3j/Ug/jkT9Ctra3xzTff4JNPPsGwYcPg6emJw4cPw97eHqtXr66wmKsyTkBewtTUFM2aNcONGzcqZH+5ubmYOXMmNmzYgMDAQI3x2vM/YFlaWorNsM6cOYPY2FjxBtKoUSMAnIAw6bm4uABAhZ0fCQkJGDRoEGJjY3H06FExScjLy0NsbKx4fvTr1w+1a9fG33//jcDAQNjY2GDw4MEAhPMjKSlJrL0EhD5UX3/9NWbOnIk1a9ZgypQpFRIvYy/TunVryGQyXL9+vUL29+jRIwwfPhyxsbE4cuSIWDtIRBoPWD169EC9evWwdetWbNu2DVZWVhgxYgQA4fxISUnB06dPizzGhAkT0LhxYx4Fi2mdru4fSqUSsbGx4rNTUfePIUOGACj6/pFftWrV4OLiIrZOYeXDCcgreHp64urVq+Vu9xcfH4/+/fsjMDAQBw4cwMiRIzXW379/X2Niw9GjR+Pq1av48ssvUadOHfTs2RMA4OTkBJlMxgkIk1yjRo3QpUsXbNq0CY8ePSrXvk6ePIlOnTohOzsbp0+f1pho6sGDB8jLyxNvIEZGRvDw8MD27duxbds2DBs2TOyw2KhRIxARYmJiAAApKSnw9PSEr68vfv/9d8yePbtccTJWUrVr10bfvn2xdetWPHjwoFz7OnfuHDp16oTnz5/j1KlTGp3Sk5KSkJmZKZ4f6mYmO3fuRGBgIAYNGiROWpi/ACsrKwuxsbEaxyEiZGdnlytWxkqiYcOG6Nq1KzZv3qzV+0d8fDxyc3NLdf9QD1WdkpKicZycnBzcvn27yEmoWelxAvIK8+fPR8OGDfHWW2/h119/xYMHD/D48WNcuHAB3t7e4iyZERERCA8PF7+4V69eRXh4OJ48eYJjx46hU6dOuHHjBr7//nvY2toiPDwc4eHhuHnzJpRKJaKjozVOmpEjR8LKygqHDh3CuHHjxJk5LSws4ODggHv37gEQZv8MDw/H7du3AQgdu8LDw5GYmKjLXxOrotauXYuUlBQMGjQIhw4dwuPHj5GYmIh///0XCxcuBABkZ2eL3/cnT54gKytL/H9ubi7WrFmD1157DXXq1MG6devw6NEjcf2DBw/E73r+8+ONN95AbGwsrl+/rtFhVr3NvXv3EBERga5du2Lv3r1Ys2YNWrduLe738uXLOvwtsapq1apVyMrKQt++fbF37148fvwYsbGx2LVrF+bOnQtAqBlXfy+Tk5ORk5Mj/j8nJwcbNmxA3759YWtri2+//RZPnz4V18fHx4tNqgqeH4mJibh06VKR50dUVBT+97//wdXVFf/88w+ePHmCmJgYeHt7Iz4+njuhM50o7v4RHByss/vHG2+8IS7Pf/9YvXo1OnbsiD/++APJycmIjo7G1KlT8eDBA3h7e+vwt1SJSTcCsOGIjo6mIUOGiDNlql9dunShiIgIIiJq2rRpkTNm/vXXX+Ts7FzkOgDUo0cPiomJEbfNTz0u9bFjxzSWu7m50ahRo4iIqF+/fkXu9+uvv9bJ74ax4OBgatWqlcb3z9jYmEaPHk3Z2dl07969Yr//Fy5cKHYdAPr888/FOQ+Sk5PFY+bm5lKtWrWoevXqlJmZKS7PyMggmUxGGzZsKHYWWwBUs2ZNKX5VrAr677//qE2bNhrfPyMjI/Lw8KCMjAyNWZoLvi5evPjS82P+/Pn0559/EgCNuaOUSiU5OjqShYUFpaamistzc3PJ2NiYVq9eTU+fPiUvLy8yNjYW92diYkKfffZZoXlEGNOWw4cPF3n/GDVqFGVlZZXr/rF48WLx/vHo0SPxmK+6f3z77bf04MEDGjlypMZzn729Pf36669S/JoqJRnR/w8FwF4pOTlZnGWzWbNmqFOnjrju+fPnUKlUhd5jZWWFzMzMItcBgLGxMSwtLZGSkgJra2uN8a4zMzORlZUFW1tbyGQycXlaWhqICNWqVUNqairy8vIK7dfCwgLm5ubl+biMlcq9e/dw//592NraomnTpmInXJVKhefPnxf5nurVqxeq5s7P3NwcMpkMmZmZhcaKT01NBRGJx1F7+vQpLCwskJeXh9zc3CL3K5PJipy7hDFtiYqKwpMnT2BpaYkGDRrA2toawMvPD/U1vjjm5uaQy+XIyMgodJ8o7vx49uwZzMzMYGFhAUC4r12/fh1GRkZo166d2DmYMV1S3z9sbGzQrFkzye8f6uenR48eISYmBnK5HG3atHnpnCSsdDgBYYwxxhhjjOkM9wFhjDHGGGOM6QwnIIwxxhhjjDGd4QSEMcYYY4wxpjOcgDDGGGOMMcZ05v8AZyjDXTVrTM0AAAAASUVORK5CYII=\n",
"text/plain": [
"<PIL.Image.Image image mode=RGBA size=800x200 at 0x7FEBFC9CE850>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"C12M1 = Chem.MolFromSmiles('c1ccccc1-c3n[nH]c2ccccc23')\n",
"C12M2 = Chem.MolFromSmiles('Nc3n[nH]c2cc(c1ccccc1)ccc23')\n",
"C12M3 = Chem.MolFromSmiles('c1ccccc1-c3n[nH]c2ccc(N)cc23')\n",
"C12M5 = Chem.MolFromSmiles('c1ccccc1-c2c(N)c3ccccN3n2')\n",
"Draw.MolsToImage([C12M1, C12M2, C12M3, C12M5],\n",
" molsPerRow=4, legends=['C12M1',\n",
" 'C12M2',\n",
" 'C12M3',\n",
" 'C12M5'])"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "sought-proposal",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0.38805112308175205\n",
"0.7545724621028375\n",
"0.393459707568025\n",
"#############\n",
"0.4411764705882353\n",
"0.6129032258064516\n",
"0.2894736842105263\n"
]
}
],
"source": [
"print(AtomAtomPathSimilarity.AtomAtomPathSimilarity(C12M1, C12M2))\n",
"print(AtomAtomPathSimilarity.AtomAtomPathSimilarity(C12M1, C12M3))\n",
"print(AtomAtomPathSimilarity.AtomAtomPathSimilarity(C12M1, C12M5))\n",
"print('#############')\n",
"print(get_tc(C12M1, C12M2))\n",
"print(get_tc(C12M1, C12M3))\n",
"print(get_tc(C12M1, C12M5))"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "interim-multiple",
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"id": "bulgarian-workplace",
"metadata": {},
"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.7.9"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment