Skip to content

Instantly share code, notes, and snippets.

@iwatobipen
Created March 31, 2024 12:00
Show Gist options
  • Save iwatobipen/56392e87a4ab28568bd7d4368bc746a1 to your computer and use it in GitHub Desktop.
Save iwatobipen/56392e87a4ab28568bd7d4368bc746a1 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": "88a6c977-eea2-4858-aa72-c280fa3864d5",
"metadata": {},
"outputs": [],
"source": [
"import rdkit\n",
"from rdkit import Chem\n",
"from rdkit.Chem import Draw\n",
"from rdkit.Chem.Draw import IPythonConsole\n",
"from rdkit.Chem import rdMolEnumerator\n",
"IPythonConsole.drawOptions.addAtomIndices = True\n",
"IPythonConsole.drawOptions.baseFontSize = 1.0\n",
"IPythonConsole.drawOptions.padding = 0.001"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "09d45f64-7155-408b-a461-d0e71f8f060f",
"metadata": {},
"outputs": [],
"source": [
"coresmi = 'c1ccc2[nH]ccc2c1'\n",
"core = Chem.MolFromSmiles(coresmi)\n",
"r1 = 'CO*'\n",
"r2 = 'Cl*'"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "3b587255-86fc-4a9d-94a7-e66fa9fe7fed",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVQTV/sH8Cch7IiyKohsgoKIBaVuFHCBVi1Yq4IWREWRuhF9W5XqrzZQS4WjVUTUF9emdQVR4FUq4laoilpcAFEBraIsIgICQZKQ3N8fYyMiRgiQkfB8Tk9Pcu/NzDM9nG/nzsoghABCHUcsFh88eNDCwsLZ2ZmWAs7VnltYtPAj9Y/iLeNpKQB1Nyy6C0CKhslkrl+/ftasWZ0Xo9kvsw9VHbpcd7m8sVyJoWSmYjZBe0KAXoAmUxMA6sR1BfwCHZZOJ60doWYwRlGHOX78uI2Nja2trampqaWlpVgs/uOPP+zt7U1NTTtqFY2kkf2EHfssVgxiSWPuy9yTL05ufLrxmOWxoRpDO2pdCLUSk+4CkOLYtGmTh4dHUFAQk8m8deuWq6vrggULTpw40YGrmPto7o5nOwBgkcGiG7Y3Xjq85Dnw0gekT9Ce8EjwaMnjJQTwIBWSN9wbRW1WW1t75coVExMTGxubpu2pqam3b9/OzMw8deqUUChcvXq1nZ2dmZlZR6036UXSgcoDALDDdEeQfpCk3UXLJcUqZXP55gC9AAYw3vVzHo+XlZWlqqo6bNgwFgv/8lGHwb1R1GYrV6708/OLjo5u1q6hoWFvb5+RkREcHFxeXl5bW2tubs5gvDPX2mpL+RYA8ND2aJqhFAYwvjH8RkfpncdDMzMzP/roo02bNq1atWrYsGHV1dUdVRVC+P9k1Dbp6ennz5//+uuvKyoqmnXx+fyAgICCggIul6ujo7Nq1So9PT0PD48OWW+9uD6jLgMAAvQCZPh5cXFxfHy8o6MjADg5OZ08edLPz69DCkMIYxS1AZ/PX7hw4bZt265evfp2L4vFGjVq1OTJk9XV1WfPnn358mUrKysAqKmpuXv3rpTFMnWYDUYNUgaYqJiUC8sbSSMAjNQcKUPl06ZNoz5UVVWVlpaam5vLsBCEWoQxitogLCzMycnJ3d29xRhVUlJis9nUZxUVlT179lCfs7OzfXx8pCzW0dMxZWGKlAEreq8Y12Mc9bmvcl9ZSgcQCoXjxo3Lzs5evHgxXde0IoWEMYpaKzs7e+/evdnZ2W39oba29scffyxlgImlibOWtFyzULHgi/kAwACGMkO5rQVQlJWVk5KS7ty54+fnN378eHd3d9mWg1AzGKOotXbv3i0QCKiz8w0NDWKxWFdX96effnrvD4cMGZKUlCRlQF1d3b59+3R1daUcrzxbexYACJA6UV0PpR5tLx8AQFdX19nZedKkSSkpKRijqKPgmXrUWtHR0ZX/+v777+fNm9eaDH2vY8eOOTk5/f777wkJCVKGmSibUB8K+AUyrGXdunU7d+7k8XhFRUXnzp0bPHiwLLUi1BKMUSQLIyOjjro36bPPPsvLy1u4cKH0YdZq1r2UesG/u6VtNWvWrFOnTllZWY0dO3bmzJlz586VYSEItQgn9UgWAQGyXHXUIk1NzdYMYwLTR8dnZ8XOHc92LDVYqs5Ub9NaLCwsjh07JlOBCL0H7o2iLuMbw2/Umer/CP6Z82gOn/Cb9ZYJyxrE0q6aQqiTYIyiLmOg2sBt/bYxgRlfFT8ob1B4WfjJFydP1Zza/my79wNv01zTo9VH6a4RdUc4qUddSYBegCHLcNmTZff5978v+b5pV0+lntT1+QjJGQMf24zoVVNTc/jw4YsXL2ZnZy9atOjzzz/v2/c9F9g3ksbLvMuXeZfLhGVMBtOQZThUY6iLlosqQxUAHgsen6k9Y8Ay8OzpKZctQN0dxiiiWVVV1c6dOyVfZ86c2YEPhUJIDjBGEUKoXfAUE0IItQvGKEIItQvGKEIItQte8IRoVVUFqalw7Ro8ewYqKtCvH7i6gqsrKCnRXRlCrYUximhCCGzcCD/9BDU1zbsGDYIdO8DVlY6yEGozPFOPaLJkCWzfDgwG+PjAjBlgaQkNDXDzJsTEQG4uqKrC8eMwcSLdVSL0fhijiA4JCTB9OigpweHDMH36G12NjeDnB3FxoK8PeXlgYEBTiQi1Fp5iQnT48UcAgMWLm2coALBYsHcv9OsHFRWwY4f8S0OorTBGkdzduQPUm0iWLm15gKYmBAYCAMTFya8qhGSFMYrkjnodXu/eMGDAO8e4uQEA3LkDtbVyqgohWWGMIrkrLgYAkP6KY0tLAACxGEpL5VARQu2BMYrkjscDANDQkDZG8kh83BtFHzyMUSR3VES+fCltTF3dqw89ZHwJKEJygzGK5M7ICACgqEjaGKqXwYA+feRREkLtgDGK5M7JCQCgpAQePnznmEuXAABsbEBbWz5FISQzjFEkR2fOgEgEgwe/Oke/Z0/LwxobYe9eAIBp0wAACgvh/n15lYhQm2GMIrnIzwdPT/DwgF27gMGAlSsBADZuhMzMFgZ//z3cuwfa2rBoEQDAkiVgawvLlrVw9z1CHwCMUdTJqqpg+XIYPBhOngQdHVBWBgCYPx+8vKChAdzdYd06KCgAkQhevoSLF2HaNIiMBAYDYmLA2BgaGsDYGEQiiI6GgQNh3z4Qi+neJITeRBDqJCIR4XKJoSEBIEwm8fcnT5++7uXzSWAgYTIJQPN/dHXJoUNvLCori3zyyaveoUNJerqcNwUhKfDRJKhznD8P//kP3LoFADBmDERFwUcftTAsOxvi4uDqVSgvBzU1MDaG8ePB1xd0dFoY/L//QXAwPHoEAODpCVu3vucafoTkg+4cR11JampqUFBQWVmZlDEPHz4MXrToub4+ASAWFuTo0Y6sgMcjERFES4sAEA0NEhJCams7cvkItR0eG0VtwOPxTp482atXr3f1rl+/fsyYMUcTEze7ukJ4OOTlvTrb3lE0NCAkBHJzwccH6ushMhLs7G4mJrbmpwKBIDk5ed++fYWFhR1ZEur2MEbR+5WUlPj7+9+4ccPU1NTExERJSSkxMXHFihWSAYSQ//3vf2PGjNm6datAIPDy8loYHQ1r1oCaWqcUZGYGR45AZiaMHFmvqfnxtGkjRozIbPGk/78aGxvHjh27YcOGv//+e/jw4SdPnuyUwlD3RPfuMOoC/vrrL2dn5759+06ePNnBwWHEiBFWVlZeXl719fWEkJs3b3p5eRkZGRkZGU2YMOHatWvyq0wkSuVyDQ0NAYDJZAYGBj5tehariYSEBBsbG4FAQAjZt2+fo6Oj/IpEig5jFLXWixcvfH19BwwYcODAAZFIRAipra0NDg42NjY2MjJycHCIj48Xi8XyL6y2tpbD4aiqqgKApqYmh8NpaGhoNiY4OPjbb7+lPldVVQFAdXW13CtFigkn9ai1du7c+ejRo5iYmLCwsJycHADQ0NAoLCxksVjz589PT0+fPn06g8GQf2FaWlqhoaE5OTne3t48Hi8sLMze3j4+Pr7pmLKyMmqnFQB69eqlqqpaVlYm/1KRYqI7x1HXcPjwYSMjoxMnThBCvvzySwcHh9raWkJIXl7eo0eP6K7utTNnzgwePJj62x4/fnxOTg7VPn/+/B9++IH6LBAIGAyG9OsNEGo9jFHUKgKBYO/evdSc/caNG3I9ANpGAoFg8+bNOjo6AKCsrBwREUEI2bhxo6enJzXgypUrenp61HEJhNoPL79HiqmysjIsLGz79u1HjhyZOnVqeXn5gAEDfvnlFycnp6VLl7q6uoaHh9NdI1IQGKNIkeXn5w/4941Ply5dioiIqKio8PDw+L//+z8VFRV6a0MKA2MUvXblypXQ0NCCggIbG5vIyEg7Ozu6K0KoC8Az9eiVhoaGWbNmsdnsW7dujRw50tvbm+6KOphYDEIhiETvHCAUglAox4KQosC9UfRaY2Mji8UCgGfPnhkaGtbW1mppadFdVIeJiYHgYHBxgfT0Fnr5/Fe3XN29CwMHyrk01LXh3ih6jcpQAMjPz+/Tp48iZShCnQdjFDUnFAq/++671atX010IQl0Dxih6w8uXL2fMmDFgwIDg4GC6a0Goa8AYRa89ePDA2dl54MCBu3btouW2ToS6IoxR9Mrjx4+dnJzKysru378/c+ZMHx+fe/fu0V0UQl0Ai+4C0IdCT08vLS2taUvfvn3pKqbz8Hhw82YL7XipE5IZxih6RUNDY9iwYXRX0emuXwdHR7qLQIoFYxR1L717w2eftdAuFsP+/XKvBikEjFEku8TExKbHT21tbSdPnkxjPa0xYABwuS208/kYo0hGeIoJyU5TU1PnXykpKekt3h6EkKLDvVEkOw8PD+pDfX39Dz/8sHfvXnrrQYgWuDeKOgCXy3V1de3fvz/dhSBEA9wbRe0lFoujo6P37dtHdyEI0QP3RlF7JScn6+npjRw5ku5CEKIHxihqr19++eXbb7+luwqEaIPPG0Xtcu3atRkzZhQUFCgpKdFdy3sIhfDyJbBYoKHR8oCaGgAALS1g4t4Fags8Nora5enTp9u3b//wMxQAXryA8+cBABwdwcqqhd7Tp0FJCaZOlX9pqGvDvVHULunp6WvXrgWAn3/+2dnZme5ypLl4ET75BADA3h6yskBZ+Y3e3FywtwcNDeDxaKkOdWE4e0Eyqq+vDw0NHTt2bHp6enp6uqur67Jly2pra+mu6/1yciA6+v3DDh48OHXq1AULFvz555+dXxTqwjBGUZsRQn777Tdra+uwsDBCiKmpab9+/Qgh0dHRtra2Bw4c6KgpjlgsLikp4fP5HbI0Ss+eoKICoaFQVCRt2O7du1etWjV79uzRo0d/+eWX+fn5HVgDUjQEoba4du3a6NGjqT8eJyenv/76i2q/evVq0/aLFy+2c0Xnz583Nzd3c3MzNTUNDQ1td+Hkr78IALGyIsHBBIBMmfJGb04OASAaGq++uri4xMTEUJ8XLly4cuXK9heAFBXGKGqt4uLioKAgJpMJAMbGxrGxsSKRqOkAsVgcFxdnamoKAAwGw9vb+9GjRzKvbvjw4Xv37iWEVFRUKCsrP3z4sJ31S2K0vJzo6BAAkpT0urdZjI4fP37Tpk3U5/Dw8C+//LKda0cKDGMUvR+fz4+KiurRowcAqKiosNnsFy9evGtwXV0dh8NRU1MDAE1NTQ6H8/LlSxlWOmrUqFWrVhFCnjx5oqqqWlZWJvsGEEKaxCghJCqKABBTU1Jb+6q3WYympKTo6uouWrRo2rRp5ubmkyZNaufakQLDGEXvkZycbGlpSc3WPT0979+///aYsLCwa9euNW0pKiry9/enftWvXz8ul9vW9V6/fl1dXd3Jycne3j4hIUH2DfhX0xgVCIidHQEgK1a86m0Wo4SQ+/fvHzx4MCsra/Xq1YsWLWp/AUhRYYyid8rLy/vs30cc29ranjp1qsVh1IlsJpM5b9680tLSpl3nzp0bMmQItYSxY8feunWr9Wv38fEJDAw8ffr0tGnTzMzMnjx50q6NeTNGCSEZGYTBICwWuXGDkJZilCIQCKytrQ8fPtzOtSMFhjGKWvD8+XM2m01dVK+rqxsVFdXY2PiuwdJn8SKRiMvlGhoaUlHr7+//9OnT9xbw8OFDBoNRXV1NfR03btzGjRvbuVHNYpQQ4u9PAMgnnxCxuHmM/v333ytWrPj555+dnJy8vLyaHQVGqCmMUfQGgUAQGxurr68PACwWKygo6NmzZ635YUFBgbe3N7XjaWVlFRcX17S3srIyJCRERUUFAHr16hUREcHn86Usra6uTk1N7dKlS4SQ+vp6Gxsb6jqqtnr5kiQmEn9/snVrCzFaVkZ69SIA5NCh5jFaWVm5Z8+eyMjIP/74QywWy7Bq1H1gjKLX0tLS7OzsqCh0d3fPyclp6xLOnj1rb29PLWHcuHHNZvF3796dNGkS1TtgwIATJ05IWVRcXJyZmdnIkSNNTU0XLlwoZXf4bfX1JDmZ+PsTbW0CQACIk1MLMUoIiYl5da7p6tWWJ/UIvRfGqMJKT09fuHChn5/f5s2b37s/lZ+fL9mXtLa2brYv2SZCoTA2NtbAwECyP1teXt50QFpa2qBBgyRhffv2bSlLe/78eetP9PN4r9JTS+tVegKQQYMIh0Pu3Gk5RkUiMnw4ASB+fhijSEYYo4opLi7O0NAwNjb2zJkzhw4dkjKyqqoqJCREVVUVALS0tDgcTkNDQ/sLqKysZLPZLBYLAHR0dKKiooRCoaRXIBBERUX17NkTAJSVldlstuQwqEzrIlwu8fQkqqqvopPJJMOGEQ6HFBS8HtZijBJCMjMJk0mYTIxRJCOMUcVkYWERHx8vfQx18qd3796Skz/tvzazmbt3706cOJHa8Rw4cODJkyeb9lZUVEhOZOnp6Uk/kfW20tLSHTt2TJ/OZrFepaeSEhk3jmzbRkpKWhj/rhglhAQFvVoCxiiSAcaoAqqsrGQwGFeuXFm3bt26deseP3789pgLFy44ODhQATdixIjMzMzOqyc5OVnymiZ3d/e8vLymvdevX3dxcaF6HR0d//zzT+lLKyoqio2N9fT0pHZ1AcDWtt7ZmURFkTevtmpOSow+f04MDDBGkYwwRhVQbm4uk8n09PTcuXMnm802NjZuOmV+/Pixv78/g8EAABMTEy6XK4cz0dQsXltbWzKLb3YfVHJysrm5ueQi/wcPHjRbQmFhYWRk5IgRI6jKAUBNTe2LL7747bffqqpqWlNDbi5xdyezZ7fce/gwcXcnn38u0+ah7g1jVAEVFRUpKSlJJsiDBg2i7gIqLy+fO3euuro6AGhoaHA4nPr6enkWVlJSEhQURM3i9fX1m83i6+vrIyIitLS0AEBdXT0kJKSmpub+/ftRUVHOzs6S9FRXV/f09ORyuVJuSG2P3Fyyc2dnLBgpLIxRBSQUCvX19a9evUoIEQgEvXv3Pn369MOHD6nLNhkMhp+fX4szffnIysr6hHp+MsDQoUMzMjKa9hYVFc2cOZMKTeq4LaVXr16zZs06fvx4p0b/06dES4uwWOTNohCSBmNUMW3atGnAgAEbNmzw8PAYPny4QCCYMmUKNaHev38/3dURQkhycrKZmZlkFt/sAU5XrlyxtLQ0NDTU1NT09/dPTk7ukOsHWmPNGgJATEzI8+fyWSHq8jBGFdbZs2d//PHH3bt383g8Qoi7uzsAjB8/nu66XuPxeJJZvIaGRkhISK3kgUuErF69GgB+/PFHOVclFJLRowkAmTpVzmtGXRU+/V5hjRs3bu3atfPnz9fQ0AAA6kVJktn0h4CKztzcXG9v7/r6+sjIyAULFkh6qXk9U+5v6WSx4PBh0NGBY8dg5045rxx1SRij3VRCQoKdnR2Hw6G7EDAzM4uLi8vMzBw1atSKFSvoLgcAoF+/VwG6bBlkZ9NdDfrgYYx2U1VVVXl5ecXFxXQX8sqIESMuXbo0bNgwugt5Zfp0mDcPGhrA1xdevqS7GvRhwxhFqGXR0WBjA7dvw6pVdJeCPmwYowi1TFMT4uJATQ1iYiApie5q0AcMYxShd7K3h4gIAID16288efKE7nLQBwpjFCFp2GxYvDjp5s1R/v7+IpGI7nLQhwhjFCFpGAwIDR2tp6d34cKF9evX010O+hBhjCL0HgYGBgcPHlRSUgoLC7t48SLd5aAPDsYoQu/n5ua2cuXKxsbGr776qrKyku5y0IcFYxShVlm3bt2oUaMeP34cFBREdy3ow4IxilCrsFis/fv39+zZMyEhYffu3S2OuXXr1vz587ds2SJpSUlJ8fDwGDVqVHh4OJ6hUlQsugtAqMuwtLTctWuXj4/PsmXLnJ2dbW1tm/bu2LEjMTGRxWJJ4vLevXtfffXVvn37LC0tAwICVFRUVq5cSUfhqHPh3ihCbeDt7T137tz6+nofH5+Xb94lumjRotTUVA8PD0nL7t27fXx8pk6d6uDgEB4eHhsbK/d6kTxgjCL0HiUlJUKhUPJ169atVlZWubm5vr6+0n949+5dR0dH6rODg8ODBw/4fH4nFopogjGKkDQNDQ0TJ050cXGRPMZFS0uLeh91eXm59N/yeDzqlS0AoKmpSQjh8XidWi2iBcYoQtKsWLEiOzv7+fPn1Pv4AOD8+fN5eXkAEBAQIP23ffr0kURtWVmZurq6jo5Op1aLaIExitA7paSkbN++XVVVNS4urkePHgDw7NkzPz8/QkhoaGhgYKD0n48fPz4pKUksFgNAYmLiuHHjJC/mQ4oEz9Qj1LLi4uI5c+YQQiIjI6lDnISQefPmlZaWurm5ff/9983GFxUVhYeH5+TkPH/+/Ouvv169erWvr29MTMynn35qYWGRkJCQmppKx3agTocxilALxGLx7NmzKyoqJk6cyGazqcZNmzadOHFCX1+fuje02U969Ojh7u5OvfMKALS1tdXV1TMzM1NTU2tqasLCwoyNjeW6DUheMEYRakF4ePi5c+d69+69b98+aiZ+/fr1NWvWMBiMPXv2tBiIOjo63t7ezRpVVVUnT54sj4oRffDYKELNXb16dd26dUwmc//+/b179wYAHo/n6+srEAjYbDbGImoGYxShN1RXV8+YMUMoFIaEhEhm6EuWLLl3797gwYPxWXnobRijCL1h8eLFDx8+/Pjjj0NDQ6mW+Ph4LperqakZFxcnuQ4UIQmMUYRe27Vr16FDh7S0tA4cOKCiogIADx48WLBgAQBs2bKl2U30CFEwRhF6paCg4JtvvgGA//73v9bW1gDQ2Ng4a9asFy9eTJs2bf78+XQXiD5QGKMIAQDw+XwfH5+6urqAgAA/Pz+qce3atZcvX+7Xr9/OnTvpLQ99yDBGEQIAWLly5c2bN62srCRPC/3zzz83bNjAYrEOHTqkq6tLb3noQ4Yxij4USUlJDQ0NtKz6jz/+iImJUVZWPnDggOSmT19fX5FIxOFwnJ2daakKdRUYo92Uu7v70aNHFy9eTHchAAB37tyZMGHClClTNm/eLP+1FxcXz549m7rpc/jw4QBACAkMDCwpKXF1dV29erX8S0JdC97F1F2UlpZK/g0A5ubm5ubmdBYEAADPnj1bu3bt7t27RSKRgYGBkZGRnAsQi8Vz5sypqKiYMGHC8uXLqcYtW7YkJyfr6Oj8/vvvb9/0iVBzBCmigwcP6vzLy8uLEOLm5gYA1tbWdJf2ikAgiI2N1dfXBwBlZeWgoKBnz541HTBt2jQAmDhx4sOHDzuvjJ9++gkADA0NS0tLqZbs7Gw1NTUASExM7Lz1IkWCMaqYfv7552+//bZpS1RUFPU/TkdHx+zsbLoKo6SlpdnZ2VH1uLu75+TkNO29c+fOxIkTAaBfv37UmEGDBnE4nDt37nRsGVevXlVWVmYymadPn6Za6urqbGxsAGDp0qUduy6kwDBGFVNwcPDGjRubtojFYm9vb+qSchaLFRwc/Pz5c/kXlp+f7+npSYWjtbV1XFxc096KiorFixezWCwA0NPTW7p06cyZM6lzPpQhQ4aEhYU1i13ZVFdXW1hYAEBISIikcd68eQBgZ2dXX1/f/lWgbgJjVDF5e3sPHTp0yJAhI0aMSE1NlbRXVlay2Wwqp3R0dKKiooRCoXxKqqqqCgkJoV6/oaWlxeFwGhoaJL1CoTA2NtbAwIBK+aCgoPLycqqroaEhLS2NzWYbGhpK8tTCwoLNZmdkZIjFYtnqod6k5OTkxOfzqZb4+HgAUFNTu3XrVjs3FnUrGKOKqbS0tLCwkBBy/PhxTU3NsrKypr2SWTMADBw4MCUlpVOLEYlEXC6XelQSk8n09/dvVs+ZM2fs7e2pesaNG/euYw6NjY0ZGRlsNrtPnz6SPDU1NaXyVCQStb4k6kXzWlpa9+7do1qKioqoN3zExsbKvKWoe8IYVXxmZmZNd0glkpOT+/fvT4WRp6cnFbsd7sKFCw4ODtRa3Nzcbty40bS3oKBA8oxOKyurZnP8dxGJRFSe9u3bV5KnBgYG/v7+ycnJrdm/pq6r379/P/VVKBSOHj0aAKZOnSrDNqJuDmNUMSUkJFBH9/Lz89XU1N51slsgEERFRVEva1NWVmaz2S9evOioGh4/fuzv708989jExITL5TadgNfV1XE4HGqOr6mp2WyO33q5ubkcDoe6BZ6ip6dH5alktt4iyX4oIWTNmjVUkbQcL0ZdHcaoAhIKhb6+vn369HFzczM0NNy2bZv08SUlJUFBQdQFkvr6+lFRUY2Nje0pgMfjcTgc6plyGhoaHA6n6RkbsVjM5XKpiTmDwfD395dcbNQeVJ42fQiTjo6Ov79/XFxcXV2dlB9euHBBSUmJyWSeP3++/WWgbghjVGHV1tbevHmzpqamleOzsrI++eQTKoCGDh2akZEhw0rFYnFcXJyZmRkVkd7e3s12hK9cuTJy5EhqLcOHD798+fLbC6murt6zZ8+TJ0+or42NjUlJSevXr4+Li2vNhD03NzciIqLpHZwaGhqenp5cLvft/xqVlZWmpqYAEBYWJsP2IkQwRlEzycnJVAhSB0z/+eef1v/277//loTXsGHDmgXxkydPJHP8vn37NpvjS6xZs8bc3Lxnz56Sq9+nTJni5uYWGRnp5uY2ZcqU1tfz4MGDqKgoZ2dnyWuN1dXVqTytrq4mhIhEIi8vLwBwcXFp5w446s4wRlFzPB4vIiJCS0uL2o8LCQmpra2V/hPqsACTyQQAIyOj2NjYpufNmy5QXV1d+gLz8vIaGxtdXV2pGG1sbDxy5Ah1lPPp06cA0Oxmp9b4559/Nm7cOHr0aKpCAFBVVXV0dLSysgIAbW3tTr1RCik8jFHUMukniJq6cOECFZGqqqrfffdds4lzcnKy5Ob91u/eSmK0qerqagaDUVVVJdMGEUJIeXk5l8v19PRUVlaWTPmbXn6PkAwwRpE0mZmZI0aMkH4os7a21tjY+O1LprKyslxcXCQHW9PT01u/3hZjdP369Z9++qkMW/G20tJSPz+/IUOGhIeHd8gCUXeGMYre470XzxNCJHccUSoqKthsNnXqX09PT4ZT/2/HaGJiorGxcYdc3CoWi1NSUiZOnBgdHU218A7LYsEAAAMqSURBVPn8X3/91cXFRXJzPUKth88bRe/BZDJnz55dWFjI4XCUlZV///33/v37h4aG8vl8yRjqJk4AEAqFW7Zs6d+/f3R0NJPJZLPZ9+/fX7ZsWTsfNxcbG7tixYrU1FTJ/QLtsWHDhiNHjigpKT148IBqWbJkyfXr1+vq6iQPEkSoDejOcdSV5OfnS246evvBImlpaZLLNt3d3W/fvi3DKgoLC9PS0uzt7cPCwjIyMvh8/ty5c/X19ffv35+WlpaWlibDKaYW/fDDD8uXL2/a4uXlxeVyO2ThqFvBxzajNqCi8+zZs8uXL8/NzfXx8RkzZsx//vMfJSWlrVu3pqamAsDAgQM3bdo0adIk2VaRm5ubkpIyatSo4uLiEydO2NraqqioTJ06NT09nRrQu3dv6imlCH0gMEZRm40fP/769evbt28PCwu7cOHChQsXqHZdXd3Q0NBFixZRT5CSzRdffPHFF180bYmNjW1PtQh1Njw2imShrKy8bNmywsJC6rEjDAbDzs7u7t27wcHB7clQhLoi/ItHMhKLxbq6ujdu3CgqKmIwGP369auvrxcIBNSToRHqPnBvFLVZTU3NL7/8YmFhQT0r3tTUVElJaezYsdbW1iYmJv7+/kKhkO4apcnJyfHx8Tl69OjJkyd9fHyePn2anJzs4+OTlZW1bdu2OXPm0F0g6mJwbxS1WUpKCvVKEuruTABITEwMDAz08/NraGhwdHQ8duzYjBkz6C1Sir59+wYFBUm+amtrf/TRRxoaGlQjHpRAbYV/MajNZs6cCQBbt26VxKjkffcqKipKSkofeBLp6uq6u7s3bTEzM5M8kAWhtsJJPeoYYrH4m2++cXZ2Hj169JQpU+guByH5wRhFHYPBYHz22WdfffXVqVOnbt++TXc5CMkPxijqGFSMstlsLy+vX3/9le5yEJIfjFHUAQ4dOpSZmQkAjY2NeXl5klvsEeoOGIQQumtAXcydO3c8PT1ramoaGhoMDQ0TEhKqq6sDAwPV1NRevHgxePDgI0eOUK/JQ6g7wBhFbSYSiWpqaiRfe/TowWKxCCHFxcVqamp4wzvqbv4fuj4buxKE4eYAAAD6elRYdHJka2l0UEtMIHJka2l0IDIwMjQuMDMuMQAAeJx7v2/tPQYgEABiRgYI4ANiXiBuYGRjSACJM3MwKABpJgYWDZA8ow5IpZYgWJARVZDNIQNIMTMzYmOAFDKzsDtYgGhGZtxK4ALcQPMZmYAiDKxsGUxs7AnsHBlMHJwJnFwJXNwZTNw8CTy8GUy8rAk8HAlOTEA9bKy8PBzsbKycHDzcXOKbkHzFwPfQbZk9A4OCAwMcBIDZQPH9CHEHoJoOGHs/jM3aJ3wAwoez7UFszasNQLEDYPbDY1PtYOK8PHfsz4XMBrPzi9WAZjgsQxLfj6R+P5I5YLYYAKjqNXD0Ty0vAAABSXpUWHRNT0wgcmRraXQgMjAyNC4wMy4xAAB4nI1TS26DMBDdc4q5ANaMZwx4GSCqqiogpbR3qJRl76/OkFAbtVgxLMzo8eB9XIGt6/j29Q2/y49VBYCFO8YIn4yI1QVsA/355XWCYTn122SYP6blHUiAWN/Ra489LfNlmxAMQM7HiNyBdwFtATp8bDachzlN2flwgGO4Qv0MocBw03F7J5JjxmCMzwAblVJ7xxiCcPbKX2S7QxaAnQFJx03AoD4dqokKRFXtAxcJCWFSnpYiUWfGEwcf/gNaMN41iG1sTAwefZr8GuGDsi5xcv6XJYO0OrnwuqD8PI27Qt0r1s/TmCpml09N0gfgVBjSLoRUC9Eg2xS+aFpdilg0k5hyFHWesrjEHCZKsYg5ST7zX8wx4szo9ZyQZIaKOdBkvq2QmGvOFdrzdgB1X/0A+km3DJNgNY4AAADoelRYdFNNSUxFUyByZGtpdCAyMDI0LjAzLjEAAHicXY8xTwMxDIX/CiNUieVnJ3cXd+zCgIAdMZTQrcchVAak/njuSlWnTMnnL356WW32tHraUEWtVV4+7l+Xs+LmeBuZ+hwSSQ7r5hpBUgqnILSIoKfp1TAKKeecwl+E8xlBzF1GwClh3pSsepagHgUyO4FmDWuhjnkI7Cpe3GU3/kuOS/Rd2B6m8flr+jSm9+9x/HnYvu3m/864mMfpsKNkcFCT64fiTkwdYMmBLTtk6xw66x16GxwGKw7F0DQCG5pOgKGpATG0RfT4C1WZeNzLxMIEAAAAAElFTkSuQmCC",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x73324e622a40>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"m = Chem.MolFromSmiles(f'{r1}.{r2}.{coresmi}')\n",
"m"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "a7e189c9-d262-41f2-bca4-7b6e00d748c9",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAACWCAIAAADCEh9HAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxN+f8H8Pdd2jeVSlGhlMoSEzOWDKNsoZRiRNbhm/lqIb7GDNEYyloMkmWkBRUmFTMIU0RJkkq79r1u+3bvPff3xxn317S5rbe5vZ9/fB91zufc+z7z7fFyls9C4XA4gFCfev78OYfDkZSUnDJlCl8K2Fe4L4gRtEtp17bh2/hSABpS6PwuAAma+Ph4KysrABAVFY2Pj5eWlu7zr+AA5171vXtV95KakipYFcNow3REdaxkrZbKLCUbFDOL05vTK1mVff7VCLWHMYr62L1798gfmpqa/vzzT0tLy779/NyWXIssi9iG2NYbYxtifSp9FksvDhgbIEmV7NtvRKhrVH4XgAQKh8MJDQ0FADI979y507efX8GqMEwzjG2IlafLu49yz5mQw57KLp9U7jvad7Tw6Ac1D86Wnu3bb0ToszBG0efV1dU9fvz4w4cPn20ZExOTn58vLy/v4uIiJiYWGRlZXFzch5XsLtid25IrS5N9rvXcXtFeTViNClR5ury1nHX0+GiPUR57R+zt4vCCgoLw8PC0tLQ+LAkhjFH0eXv37rW2tj59+vRnWwYHBwPA8uXLZWRkjI2NCYLg3uP3XhmrzL/SHwBcVFzGi45vs1eRrminaEcBSmeHHzt2bN68eZcvX16yZImdnV1fVYUQxij6jFevXj148MDW1vazLdlsdlhYGACYmpoCgIWFBfTpff2z2mfNnGYhitA6uXU9OJzFYr179+7GjRtPnz799ddfa2pq+qowNMThKybUlebm5s2bN585cyYpKemzjZ8/f15WVqaiomJgYAAAc+fOlZOTS0hISE1N1dbWJtuUlZV9/Pixw8MpNEqzdnNnHz5JbNLbxrcAoCuqK0OT6cG57Nu3j/whOTl5xIgRYmJiPfgQhNrDGEVdOXz4sL6+vomJCS8x+vvvvwPA8uXLqVQqAAgJCS1duvT69evBwcF79uwh2zx8+HD37t0dHj5MdtgHn04fv0ZqRVawKgBARUilBydCevbs2a5du1JTU0NCQoSEhHr8OQi1hjGKOvX+/XtPT893797x0pjJZP75558AsHz5cu5Gc3Pz69evBwUF7d69m0KhAICiouK0adM6/AQJaQk5SbnOPl+GJtPMaQYAEapIt86itZkzZ4aFhYWGhtrY2CQnJ0tJSfX4oxDiwhhFnbp69SqbzZ4wYQIANDU1EQQhJyfn6uraYeOnT59WVVVRKJQjR46QiQkA5Bi5/Pz82NhYMj2NjY2NjY3bH56VlfXbb7/9r/F/kydP7qweaao0ANSwe/5MU1hYeMSIEVu2bHF2dn779u2cOXN6/FEIcWGMok6dPn2a+3b+2LFjGRkZnWUofLqj53A4z58/b7/3zp07nV2EAsCOHTueP39eWVmpoaHRRYyOFB4JAOnN6byfAheLxbK0tNyzZ8/06dMfPnzIYDDGjRvXg89BqD2MUcSTESNGtLS0dLa3sbHx0aNHABAWFjZ27NjWu7KzsxcvXnzv3j0XF5fOHkceOXJESkrKzMys6xqmi08HgLyWvPTm9HEi3QtBOp2+ZcsWR0fHjIyMsWPHBgUFKSsrd+sTEOoMxijiiY2NTRd7Hz58WF9fr6Oj034uksmTJ0+YMCExMfHZs2cd3s4DAI/PKA0lDVWEVAqZhSdLTnqqefJYOZeJiYmJiUl3j0Los7DfKOoD3F73He41NzcHgNu3b/fyW+gUOjlI6VL5pcvll9s3yG7J7uVXINQDGKOot2pra58+fQqdx+iKFStoNNqff/7Z+x7v3yt8bzHMggDiu9zvjNKNLpZffFTz6G7VXbcSt2kp07STtMtZ5b38CoS6C2/qUW+9fftWW1tbTU1tzJgxHTZQUlKytLRMTk5++/bt119/3ZvvogL11phbh4sPHy85Hl4bHl4b3nqvjqhOKat0OH14b74Coe6i4LTNiO8SEhJevXp1/vz5KVOmzJkzZ+PGjZ89pIpdFV4bntSYVMGuEKOIjREZ86XEl/pi+uTeF3Uv0prTDMQNJopN7OfaEcIYRYNATEwM+VgAAISFhR0dHflbD0LdgjGKEEK9gq+YEEKoVzBGEUKoVzBGEUKoV7DDE+K3jAy4fx/S06GyEmRlQUMDFi+G8W0nt0do0MJXTIh/GAzYvh1u3YI2f4QUCpiZwYULoKTEp8oQ6gaMUcQnVVUwezYkJYG0NHz/PSxYALKyUFUF4eFw9ixUVYGGBrx4gUmKBj+MUcQnNjbg4wPq6vD0KbQZ/pSTA/PmwcePYGoKv//Op/oQ4hXGKOKHzEzQ0gKCgEePwMiogwaRkUDOqfzuHUyaNMDVIdQt+KYe8UNQEBAEjB/fcYYCgKEhfPEFAEBAwEDWhVAPYIwifoiOBgAwNOyqDTmJCdkSoUEMYxTxQ0EBAMDo0V21IR+Yki0RGsQwRhE/1NcDAEhIdNVGUhIAoLZ2IOpBqBcwRhE/kBHZ0NBVGzJApaUHoh6EegFjFPHDiBEAALm5XbXJywMAwIXn0KCHMYr4gVxs+eXLrtpERQEATJ8+EPUg1AsYo2gANTXBX38BAJiZAYUC797B69cdt3z//u8YNTcHAHj+/O/HqQgNPhijaKCEhICeHixYABkZoKcH5FrH330HVVVtW9bWwubNwOHA/PlgYAAlJWBiAlpa4OUFBDHwhSPUNYxR1P8SEuCbb2D5csjKAm3tv98deXrCiBHw7h1MmwZ+flBWBgBQUQG3bsH06fD6NSgowOXLAACVlaClBYWFsG0bGBp2egGLEJ/gYFDUnyor4dAhOH8eWCyQk4MDB+D774H+aXrGrCywsoI3bzo4cPJkCAgALa2/f+VwICgIdu2CvDygUGDlSjh5ElRVB+gsEOoSxijqF0wm09vbe5K///SnT0FICGxt4eBBkJVt244gIDQUQkMhJQVqakBaGrS0wMQETE2B2u5Wqb4ejh8HNzdoagIJCXBygr17QVR0YM4Ioc5gjKJuiIyMjIiIAIBly5ZN6nzGkGfPnjk7O6enp+uNHfunlBT11CnQ0emzIvLy4McfwccHAEBVFQ4fBhubPvtwhLoPYxR1g5mZWUxMDACsWLHi3Llz7Rt8/PjR1dU1JCQEAMaMGbN3795ly5b1SymPH4OjIyQmAgB76dIPR49OmDDhsweVlJQ8fPhQWFh48eLF0tixH/URfMWEeFVYWBgbGyskJESj0f7444/6f/ZAqq2tdXFxmTt3bkhIiLS09P79+589e9ZfGQoARkbw7h14e4OCwtO6usmTJ9vY2JSWlnZxxIcPH3R1dcPDw2/dujVp0qSuGyPEO4xRxKvg4GCCIObOnTtjxozGxsYHDx603tvc3Ozv789ms1euXBkZGWlrayskJNS/BVGpYGMDKSnh+vpUKtXHx0dHR+fs2bMsFqvD5q6urhs3brx27dqdO3emTJly4cKF/i0PDR0chHizcOFCZWXloKCgGzduKCsrf/vtt20aBAcHJyQk8KW2lJSUJUuWkH/S2traYWFh7duoqalFRESQP1+7dm3u3LkDWyMSWHg1iniSnZ2dkJAgIiJibGxsYmIiKioaGRnZ5r54+fLlEydO5Et5ZHQ+evRIV1c3NTXVxMTE2Ng4OTm5dZvi4mJFRUXyZ0VFxeLiYn5UigQQxijiyd27dwHAyMhIWlpaSkpq/vz5bDb73r17/K7rH4yMjOLj493d3WVkZB4/fqyvr29vb19dXU3ulZKSavg0p1RdXd2wYcP4VykSKBijiCdkYi5fvpz81dzcHADu3LnDz5o6IiQkZG9vn5KSsmnTJjabfebMGW1t7RcvXgCAjo5OXFwc2Sw+Pn78+PF8rRQJEH4/VUD/AsnJycrKyhoaGvX19eSWlpYWHR0dZWXljIwM/tbWhTdv3hgaGsrKypaVlXE4HF9f39GjR7948eLBgwfy8vIvX77kd4FIQNA/H7RoyAsODgaAhQsXiouLk1uEhISWLFni7+9/9+5dJycnvlbXqalTp/71118ZGRnDhw8HAGtr6+rq6h9++EFYWPjixYtfffUVvwtEAgK736PPmzlzZnZ2doe7yOs7CoUysBUhNIjg1Sj6jPj4+OzsbFFRUS3uRCGfZGVlZWdnx8XFfUEuhjzosdlAEECjdTBeHwA4HGCxgEL5/7lTEOIFvmJCAADNzc2d7fr9998BwNzc/I92NmzYAIPyRVNnFi4EYWFwd+94b2AgCAsDD2NKEfoHjNEhjcPhODs7q6ioaGhofPXVV1lZWW0aEARBvqM3NTVtf/jKlSsB4N69e0wmcwCqRWhwwhgd0j58+BAXF5eSkpKfnz9t2rQ9e/a0aRAdHV1cXCwvLz9jxoz2h2tpaenq6lZUVPxFLg2C0JCEMTqk6erqkjOJAICxsXFGRkabBuSl6LJly+idPC8ctB1IERow+Cwd/S0mJkZfX7/NRjU1te+//97CwqKzoywsLBgMhpiYWD9Xh9DghTGKAAASEhK8vLyiyMU4W7G1te36QCUlpR9//LHf6kLoXwBjFEF0dLSVlZWfn5+mpia/axkIhYUQH9/B9pycAS8FCQTsfj/UXb9+3cXFxdfXdyiM6jEygvDwz7TR1oaUlAGpBgkKvBod0vz9/devXz98+HBra2tyS2ZmJn9LGgAGBqCr28H2nBzAHgeoBzBGhzQLC4vFixf35hOOHDnS0tJCEASFQqHT6RYWFrodRtRg8u23sHNnB9sDAjBGUU9gjA5pIiIiIiIiPTiQw+GUlZUVFxc3NjZytxw7dmzkyJHi4uIqKirCwsJ9WilCgxfGKOq2xsbGjIwMFotFEISZmRm5MTo6WllZedKkSZWVlQwGY+TIkdyp5hESbBijqHvq6+vT0tIIgmiz3dfX18bGhkKhkDMwFhQUtLS0jBo1ii9FIjSQcBQT6gYmk5ment4+QzMzMzMzM42MjLhbCIIoKysrLy8f2AIR4gOMUdQNOTk57TMUAHx9fVevXt1mRWWCIPLy8jpb7hghgYExinjV1NRUW1vbvqNxZWVlZGTkihUr2h/C4XBKSkoGpDqE+Aa73yNeFRQUlJSUtP+DOXfuXHNz884O+xABCAkJTZo0qf+r40lDA7BYICoKHfYjYLGgoQFoNJCQGPDK0L8ZvmJCvKqpqenwH10pKaku5i5hsVhMJrPN/T6/ZGRAaioAwNKl0H42lfR0SEwENTX48suBLw39i+FNPeJVZ3Mzs1gsZ2dnOzu7mpqa9ntpNFoXU+sPMD8/sLICKys4fLiDvbdvg5UVnDkz4GWhfzmMUcSr9peiubm5jo6O58+ff/PmTVRUlKmpaVhYWJtmZP+nASyTJydOQHLyZ9rU19c7Ozubmpru2bOnsLBwQOpC/0oYo4gnzc3Nrd/R19XVeXh4rFq1KjIyUlhYWEdHR1JSsra21tnZedu2bWlpaa2P7WzKZx6/t7CwsMPuAT2moAAtLWBrC13H+6pVq16/fm1nZ9fS0mJiYjII/zFAgwTGKPoMgiAKCwuTkpLILCMIIiwszNzc3MfHh81mL1myJDQ01MfH58mTJ4cOHZKTk4uLi1u7du2BAwcqKysBgMPhiIqK9uyrDx8+PG7cuDVr1mhoaIR/dmomnn37LWhqQkQEeHt32qaioiIsLMzLy2v+/PmnTp0qKytrPxkrQiSMUdSV8vLyxMTEoqIi8losMTFx8+bNzs7OlZWVenp6V65ccXFxkZOTAwAqlWpiYnL79u3169fT6fT79++bmZl5eXmJiIj0bBX7wsLCgwcPvnnz5tmzZy4uLrt37+6rkxIWBldXAIDdu6GiouM2dDqdRqO1tLQAAIVCGTVqVHp6el8VgAQMxijqWH19fUpKSk5ODvlmqbS09MCBAxs3bnz//r2iouKhQ4euXbs2ceLENkdJSUnt2LHj5s2bRkZGDQ0NXl5ey5cvDwwM7EEB5LhSMrzq6uokJSV7f1JcFhawYAGUl0O7Rfz+JiMjY2tru2TJEjs7u3nz5mVmZpKRilB7GKOoLSaTmZ2dnZKSUl9fDwBNTU3e3t4rV668f/++iIjI+vXrg4KCTExMWl9j1tbWnjp1qqqqivxVTU3N1dX1woULWlpamZmZVlZWxsbGiYmJ3SpDWVn5yJEj8+bNMzU1PX36tL+/fx+eIwCcPAl0Ovz2G0RGdtzg7Nmzly9fNjQ09PX1VVBQUFdX79sCkMDAGEX/jyCI4uLixMTEik/3uhEREZaWlmfPnm1oaDA0NAwICNixY4e4uHibAy9duuTv729ubu7v70+O/qRQKLNnz3737t3FixeHDx/++PHjKVOmbNu2jfdR9iUlJZcuXbp8+bKZmZmIiIijo2MfnikATJgAdnbA4YCtLXTSlQtmz55taWlZVlaWm5s7ffr0vi0ACQyM0SGBw+HU19czGAwGg1FTU9PhOPfq6uqkpKSCggLyVVJKSsqWLVt27txZVFQ0fvz4S5cunT59WkVFBQCoVCqV+o+/nJUrV86ePbumpubUqVNWVlYvXrwQFRUdP368qKjo1q1bU1NT7ezsOByOl5eXtra2h4cHm83+bM0BAQHq6urr1q3buHFjRETE3bt3+/zp5MGDMHIkJCWBu3sHe728vI4ePbpv377FixefPHlSVla2b78dCQyMUQHX1NT08ePH+Pj4tLS07Ozs7OzsrKyshISE5ORkBoNBtmloaEhNTc3IyCAf/1VVVZ04ccLGxiY+Pn7YsGG7du26fv36lClTAIBCoSgqKk6aNElTU1NERIRKpZK39mpqau7u7ufPn9fQ0MjNzbW3t7e3t08lBwwByMnJeXh4vH//fuHChZWVlQ4ODhMnTvzzzz+7rlxFRSUnJ6eurg4AcnJy6HS6vLx8D/4LFBeDpycYG3ewwpKUFJw6BQDw889QVNR2r6Ghobi4uJKSUnh4+LZt23rw1WiIwDH1gqyoqKi4uLizDvBUKlVUVFRcXLyiooJswGKxAgMDPT096+vr6XT6ypUrbW1tJT6NMJeRkVFVVW09W35dXV11dXVdXR2LxaLRaKKiomJiYjdu3Dh48GB1dbWQkJCtra2Li4uMjAz3kJCQEAcHh6ysLABYunSph4fH2LFjOyyeIAhHR8fbt2+rqanl5+e7urquWbOG93PPy4MHDyAkBP74A8iL759/htpaOHYMnJzg+PH/b7lgATx6BFu2gLo67N8Pa9aAnx/v34MQxqjgSk9P9/Pze/36NZ1OX7Ro0Wcf7UVERJw6dSo/Px8Apk+fvnv37jFjxpC7REVFR40a1ToNu1ZRUeHi4nLu3Dk2my0vL79///7//ve/NBqN3NvS0nLhwoX9+/fX1tYKCwv/5z//+fnnn6WlpTv8KDabzWAwZGVluYd3LScHfv8dAgMhKurv3vWiomBkBMuWwYoVcOJEBzGakQETJwKTCStXwq1bGKOo2zBGBVNhYeHWrVvz8vLWrVsnJCQkJSU1efLkzhrn5OScPHmS7F6urq6+c+fOWbNmkbtoNJqKioqCgkIP+n6+ffvWwcEhIiICAPT19T08PObMmdO6wkOHDl2+fJkgCBUVFWdn5y1btrR55Mq7rCwICflHeoqJwfz5YGkJZmbAjej//a+DGAWA/fvh8GGgUoEgMEZRt2GMCqCmpqbg4OBNmzaFhoZ2dpVHqqmp8fLyCgwMZLPZ0tLS69evX7NmDTkbE4VCkZeXHzlyZG+GcgJASEiIvb39x48fAWDp0qVnzpzhXuQCQGxsrL29PZngBgYG7u7u3ATnxbt3727fvp2aOicg4O+J92VkYNkysLCAhQs7mMOpsxhtbIQJEyArCwAwRlG34SsmAVRYWPjhw4cJEyYkJCScP3/+zp077SdnYrPZd+7csbCwuHnzJgCYm5uTA5DIDJWSktLR0VFXV+9lhgLAsmXLkpKSXF1dpaSkQkND9fT09u7dW1tbS+41MDB4/vx5QECAmppabGysoaGhlZVVbm5u15+ZlJR08OBBHR0dfX39n3/+OSrqkJwcrFsH9+5BSQn4+ICZWQcZ2gUxMTh3rseniIY8DhIsbDY7Li5u+/bt4uLiVlZWzs7OU6dOXbhwYWwrFy5c0NTUJP8Apk2bduPGDe6u9+/fMxiM/iisoKBg3bp15MMBFRWVixcvstls7l5yOiVy9L24uLizs3NjY2Ob84qMjHRwcFBTU+P+9Y4YMcLW1vbRo0dMJsFLDV5eHCMjzsWLHe+1s+MYGXGOHOnFSaIhCWNU0NTW1r59+3bPnj1ff/01mYx37tyhUqkvX76MjY09f/78zJkzyQxSVVV1dXXlBuibN2/ITqP9Wl5MTMyMGTO4CR4VFdV6b25u7rp167jleXt7s1isyMhIOzs7sssqd9fWrVvv3bvHZDL7o0g/P05MTH98MBJMGKOCpqKiIi4u7uLFi0pKSjExMbGxsZcuXRIXF4+JidnzaQC5uLj4jh07oqKiWl+iJiUlDUyFBEF4e3uPGDECACgUiqWlZW5ubusG4eHh3HVHpKSkuOmpoaGxZ8+e6Ojofs16f38OAEdTk1NT039fggQKvmISNOXl5Xl5eWw2e9u2bZKSkvr6+kFBQUuXLt26dev8+fOrq6sVFRWvXLmirKzc5kAxMTFdXd0Bq7O+vv748eNubm5NTU0SEhJOTk579+7lTqlHEISvr++OHTs4HI6srKyZmZmlpeWsWbN6NllUtzCZMHs2xMSAlRXcutXf34YEAcaooKmurv748SObzW5paXn06FFBQcGECRPIG3lDQ8PGxsajR48aGxu3P1BKSkpLS2uAq83Nzf3pp598fHwAQFVV9fDhwzY2Nty9enp6ycnJiYmJenp6A1lVRgZMnQq1tXD9Onx6xoBQp/BNvaARExMj/2kUFhY2MTHZunVr64ehAND6FQ0XhULp25noeKSmpnb9+vVHjx7p6enl5eWtX7+e7DnAX5qa4OEBALB9O3wa0YpQpzBGBY2wsDCPy3C6ublZWVnFxcUBAIVC4ePUG0ZGRgkJCd7e3jNmzOhwvfuBt3EjWFtDXR1YWwNONIq6hjEqgJSVlXkZDlRUVJSVlUVOKkoOh+//0jpFpVJtbGyioqJaj9nnL09PGDcO3ryBn37idylocMMYFUBycnLdWrqDSqXinMTtSUqCnx8IC8OJExAWxu9q0CCGMSqAKBSKpqYmdxa7rlGp1JEjR7afiRkBwLRpcPAgcDiwaRMUF/O7GjRYYYwKJmFh4fHjx9Pp9M8mqby8vKKi4sBU9W/0v//B/PlQWgoHD4ZitxbUIYxRgSUqKqqnpycvL0+hUNo8KiW3kFPP8T793dBEpYKPDyxYYHvx4rJT5CTPCP0Txqggo9Fo6urqkyZNGjVq1LBhw8gwFRMTU1JS0tbW5ksPp38jZWVwcFhOoVB++OGH6OhofpeDBh2MUcFHp9MVFBQ0NDTIMUKjR4/Gh6HdtXjxYjs7OyaTaW1tXVNTw+9y0OCCMYoQT9zc3KZMmZKZmWlnZ8fvWtDggjGKEE9EREQCAgKkpKS8vb19fX35XQ4aRDBGEeKVpqbm6dOnAcDW1jYtLa19A4Ig7ty5Y2pqGhMTQ25hs9m//PLLjBkzjI2Nw7D3qYDCGEWoGzZv3rxmzZq6ujpra+uWdqNEV69eHRwcnJycTK4MCADu7u4BAQEXLlz4/vvv16xZk5ycPOAlo36HMYpQ95w/f37MmDGxsbEHDhxosysgIMDb25ucAoZ08eLFX375RV9f38zMbPXq1VeuXBnYYtFAwBhF6PNycnK4P8vIyPj5+dHp9OPHj9/qckZSJpOZmZk5ZcoU8ld9ff1UnDBKEGGMIvQZ8fHx48ePJ+eQJrfo6OhQKBSCIIq7HCLa0NBAEAR3zhcJCQnuWn5IkGCMItSV+vr6b7/9tqmpCQC4I2v37dvHZDKFhYVXrVrVxbEyMjLi4uIlJSXkryUlJa1XlEICA2MUoa7Y2dmlpKTo6ekdO3aM3HL//n1PT08REZEnT56QK0p14Ztvvvn9998BgMPhBAcHz58/v98rRgOut6uQIyTAgoKCrl69Kioq6u/vT96bFxQUrF+/nsPhuLq6zpo1q037e/fuhYWFpaamenp6RkdHu7m5HThwYOHChZmZmTk5OXV1dWvWrOHHeaD+hTGKUMfy8vK2bt0KAB4eHuRKpQRBrF+/vry8fNGiRfb29u0P0dDQMDIyMjIyAgAydqdNm5aUlBQeHi4tLb1gwQLumn1IkGCMItQBFou1evVqBoNhbm5OhikAHDlyJDw8XElJ6bfffutwBkI9Pb32q+8pKyuvXbu23ytG/IPPRhHqgLOzc1RU1KhRoy5dukRuiYmJcXFxoVKpPj4+n30kioYUjNGhoqmpic1mA0BjYyNBEPwuZ1D766+/3NzcyMSUk5MDgOrq6tWrVzOZzD179nS4PDUayvCmXsC1tLQUFxdXVFRQKBRy8GJubq6oqKiEhASLxeJ3dYMRg8GwsbFhs9mHDh2aO3cuudHW1vbjx48GBgaHDh3ia3VoMMIYFWQlJSWFhYUcDofsN07+L0EQHA6nrq6uoaEBAMhLVETicDgbN27Mzc01NDT88ccfyY2XL1++ceOGpKSkn5+fsLAwfytEgxDGqMDKyclhMBhd3L+TqZqXl8dkMnlc2l7gnT17Njg4WFZW1sfHh1xkJT09fefOnQDg6emppaXF7wLRYITPRgVTcXFxZWUlL1eabDY7LS0NF2sDgMTExL179wKAp6cnueJ0c3PzqlWramtrN2zYYG1tze8C0SCFMSqAmpubi4qKeHyPxOFwWlpaioqK+ruqQa6hocHKyqqxsXH79u1WVlbkxj179rx9+1ZTU/PMmTP8LQ8NZhijAoh8Hsp7e4IgSkpK+P76vqGhITQ0lF/fbm9v/+HDBz09vRMnTpBbHjx4cPbsWSEhIV9fXykpKX4VhgY/jFFBw+FwqqqqeInRTZs2HTt2TFdXFwAoFEpVVXcjuvMAAAySSURBVFX/V9epkJAQXV1dMzOzhISEgf/227dvX758ufWgz5KSko0bN5KDPr/88suBLwn9i+ArJkHT0NDQ4QAbAKirqwMA7lxt5ABHEpvNrq6uJvtIDrDY2Fh7e/uoqCgAMDAwGPiLYu6gz9OnT3MHfa5du7akpGThwoWOjo4DXA/618GrUUHT3NxM/nDgwIFvPrl9+zYAVFRUAEBcXFyHB5JzwQ2k8vJye3v7r776Kioqavjw4e7u7q9evdLX1yf3NjQ0VFdXA8Aff/zRf9N0slisb7/9trKycsWKFf/5z3/Ija6uro8fP1ZUVLx27Vpn/yYhxIUxKmi4V3N5eXknT5588uTJkydPLCws4NPlp5eXl5ubG3ll2tpAvqxnMpkeHh4aGhpnzpyhUql2dnaZmZn29vZkHyMACAkJ0dPTKygokJKScnJyUlRUXLZsmZeXV2lpad9WcujQoRcvXrQe9Pn69euDBw9SKJSrV6/ioE/EC4xRQUOn//2gpqysTEFBofWuH3/8cfTo0RQKJTAw0Nzc/O7du63voNls9sB0xQ8JCdHR0XFwcKipqTEyMnr37p2Hh4e0tDS5Ny4ubs6cOcuXL8/Ozp4yZcp33303Z84cJpMZGhq6bdu2UaNGLVq06NKlS2VlZb2vJCIi4ujRo1Qq9fr16/Ly8gBQXV29atUqJpO5e/duExOT3n8FGhI4SLA0NzfHxcW9fv1aWFh4+vTp6urq8+bNu3//fuwnfn5+U6dOJf/f19bW9vLy4u6Kj48nX9n3U20fPnxYvHgx96vDwsJa7y0vL7ezsyMvSOXl5d3d3VksFneXt7f30qVLuYOIqFTqrFmz3N3d8/Pze1ZMZWWlmpoaADg7O3M3kvOBfvHFF83NzT09SzTkYIwKoISEhNjY2KCgoCdPnkRFRa1cuXLatGmx/3Tq1CnughaGhob37t3j7kpOTq6tre3bkiorK+3s7MgrZVlZWXd3dyaTyd3b0tLi7u4uIyMDAEJCQnZ2dmRngw4/h8xTERERbp5+8cUXzs7O5CAC3pmamgLA7NmzuWFNLtspKSmZmpra4zNFQxCFg8NXBE55eXl+fj73Dj0tLW3Tpk3Pnz9v06y5ufnmzZtXrlxpaGgQERFZvXr15s2bxcXFyb0yMjKqqqrctOoxFot19erVn376qaysjE6nb9q06fDhw62fNjx+/Nje3p5cwN3IyMjDw4Psg9W1hoaG8PDwwMDAu3fvcp/z6urqWlparl69evz48Z893NLS8uXLl/Hx8eQ1aUZGxtSpU2tra318fHB6UNQtGKMCiMPhxMbGRkZGfv311wBw8+bNsLAwHx+fDhuXlpb++uuvDx484HA4CgoK//3vf5csWUK+nqZSqYqKiiNGjOC++emuJ0+eODg4vH//HgC++eYbd3f3iRMncvempqbu3Lnz/v37AKClpXXq1KkePI5sbGx8/PhxYGBgcHBwTU0NuVFXV3fZsmVLly6dPXt2ZwdyOJyMjIxx48YBQHNz84wZM96+fWtjY+Pt7d3dGtAQhzEqmBITE83NzZlMprS0dHl5+bFjx7q+QEtOTj5x4gTZ9V1XV9fJyYnbq1RISGjkyJHkGxjeZWRk7Nu3LzAwEAA0NTWPHDliaWnJ3ctgMNzc3E6fPt3S0jJs2LC9e/c6Ojr2cvKk5ubmyMjIkJCQmzdvcl/ojxkzZtmyZZaWlrNmzeqi65Kjo6O7u7uGhkZcXBz3ZRdCPMIYFVgNDQ0vX76srKwcPXo0ldpxlwwyWTifptG7f//+mTNnyMlJ58+f7+joqKSkRLaUkJBQVVWVkJBo8wkEQTQ2NrLZbCqVKiIiIiQkVF9ff/z4cVdX1+bmZgkJCScnp71793LXICIIwtfXd/fu3aWlpVQq1dra+sSJE4qKim0+NjU1NSIi4rvvviN/ZTAYQUFB1dXVs2bNmjFjRtcnzmazX758GRgYGBgYyJ0rQE1NzczMzNLScubMmW3+a/zxxx9Lliyh0+mRkZE4YAn1AMaoIGOxWLm5uWQn9tZ9mygUCoVCERISUldXFxERKSwsJHvmA0BjY6OPj8+1a9daWlrExMTWrVu3YcMG7nWivLz8yJEjyVn1qqqqiouLGxoauKlEEMSDBw/OnDlTWlpKoVDWrl177Nix1l0vnz596uDgQF7zzps3z93dvfVIKlJpaemWLVs+fPiQkZFB/nEWFhbOnDnT1NRUVVX19OnTv/zyy4YNG3g5fYIgoqKiAgMDb9++XVBQQG5UUFBYtGiRpaXl4sWL6XR6bm7ul19+WVxcfPz4cScnp+7+F0YIMEaHgpaWloqKipqampaWFg6HQ6fTJSUlZWVlW0+3UVtbm5eX19jYSP5aXFx8/vx58qmlkpLS9u3buU8tqVTq8OHD6+rqmpqaWkdzUlLSiRMnyMegenp6p0+fbr3YRl5e3o8//kg+n1VVVT18+LCNjU1n1X78+FFRUVFOTo784ywvL4+PjyeX27x8+bK3t3dkZGS3/gsQBPHq1as7d+7cvn07Ozub3KioqKiqqpqamlpXVzd//vyHDx92ds2OUNcwRtH/q6ioyM/P5y4uEhsbe/LkyfT0dAAwMDDYtWsX+UKmvQsXLly9epXD4SgpKdnZ2S1YsIBGo6moqCgpKZH3+G5ubk1NTe3v8TvDYDC4Mdra9evXr1279uTJkx6fY1JSUmBgYEBAwIcPH8gtVCo1Pj6+9bsvhLoFYxT9A5vNLi4uLikp4XxaceTBgwfu7u4MBoNKpS5atMjBwaH9DCavX792cHBo02WKQqEkJCQ4Ozvn5uZSKJSVK1eeOHGC7F30WR3GKEEQ8+bNW7Vq1fbt23t/pk+ePHFxcaHT6du3bzc3N+/9B6IhC2MUdaCpqSk/P598qAoANTU13t7e/v7+5Kv/9evXr1mzps26IwwGQ1ZWlvtr61f/BgYGHh4eM2fO5L2ADmN0586dSUlJYWFh3AGvPVZXV3ft2jV/f/8LFy5MnjwZAEpKSs6dO/fw4cPw8PD2b9IQ6gI+DEIdEBUV1dTUHDduHHn3LS0tvWPHjps3b86aNaumpubs2bOrV69+8eJF60O4GVpWVnbkyJENGzYkJCQoKCgcO3YsOjq6WxnaXn19vbW1dUpKyt27d3ufoQCwdu3a/Pz8wsJCBoMBAAwGY9WqVaKiotHR0bhgKuouvBpFXeFwOKWlpUVFRdwxUTExMSdOnMjKygKA6dOnOzk5jR07ltzFZDKDgoIuXLjQ0NAgJCRkYWGxfft2GRmZ9q/ju/bq1avi4uIVK1Y8evRIVVVVSEjIzMxMSkrqwIEDNBqN7IzVJ2enr6/v7u7OXUWZIAgajVZVVUUOS0WIRxij6PNYLFZRUVFZWRn518JisQIDAy9evFhXV0en0xcuXLh9+/bHjx/funWrsLAQAAwNDXfv3k2O2adSqRMnTuzWJeShQ4fIzwGAefPmKSoq3rp1i7uXRqOdP3++T84LYxT1CYxRxKt37961vuGtqKg4d+5caGho625PWlpaTk5O3BmkAIBGo40bN25wPm3EGEV9Ap+NIl61+RdXXl7+wIEDZ8+eJSOSQqGsWLHC19e3dYaS+L5YHkL9CtdiQryi0Wht5nUmCOLLL798+vRpbm6urKysjIxMXV2dpKRkmwP75KUQQoMWXo0iXnE7hAJAbm7u0aNH58yZExERQaVSR48enZycvHz5cmtr6wULFty9e5fbkiCI3s+21+euXr1qZWWVk5Nz6NChXbt2AcD+/ftXrVoFABs2bPj111/5XSD6N8HLBMQrWVnZ2tpa8oI0NDR08uTJaWlpTCaT3Hv//n1PT08VFZW0tLR169YtWrSIXKlYQkJiEA6ynDVrlpqaGrkgKHn5bGpq+vXXX2/btg0AlJWV+Vwf+lfBGEW8kpWVzcvLI38mxxGFhoZy9/7888/kDyIiIjQajXyQSqPRBmckaWtra2trt95iYGDAr2LQvx3GKOIVhUJRVVXNycnp7JVRZmZmQEDAq1evXFxcxMXFKRSKmJgYTt+JBN6gu9tCg5mcnNywYcM6u0mXlZWdPXv2jBkzfH192Ww2jUbj9sxHSIBhjKLuGT16dGdJKicnZ2houHfv3pKSkqSkpPHjx7cZd4+QQMIYRd1DoVDGjBmjrq7euhsTQRBXrlypqamhUChlZWV1dXUGBgaD8AU9Qv0Bn42invD19fXw8CgsLMzKyvLz87t+/TqdTre0tFRRUSkqKvrhhx9w+k40dOBgUNQTTU1N3KnyqVQqOXqyubm5sLBwxIgRZFcnhIaI/wN/VAv6AZ1HXQAAASl6VFh0cmRraXRQS0wgcmRraXQgMjAyNC4wMy4xAAB4nHu/b+09BiAQAGJGBgjgA2JeIG5gZGNIAIkzczAoAGkmBgMNkDyjJEjWsSQ/F4Qzk/1Kcx1Beh2lFaVgZsgwaAuCNTGSponNIQPIYGZmxMYAGcTMwu5gAaIZmXErgQtwA+1nZGJkB/I0jBTMFSw0mUF2+kUyszByQATNFAxNYKKsbBlMbOwJ7BwZTBycCZxcCVzcGUzcPAk8vBlMvKwJPBwJTkxAlWysvDwc7GysnBw83Fzim5CCjoHvodsyewYGBQcGOAgAs4Hi+xHiDkA1HTD2fhibtU/4AIQPZ9uD2JpXG4BiB8Dsh8em2sHEeXnu2J8LmQ1m5xerAc1wWIYkvh9J/X4kc8BsMQAuRkZbQgI5CwAAAWV6VFh0TU9MIHJka2l0IDIwMjQuMDMuMQAAeJyNU8FOwzAMvfcrfIRDK9tx0uawQ9dOgNBaxDokPgGJI/8v7G0lmVirJT0k1str/N5LATbe+9evH/gb3BcFAK58MUb4cIhY7MEWsN09vQzQTe12rnTjcZgOQALk9IzOa2w7jfu5QtABVRwjuga48mgDsMLLYsYxjKnqKvYLOActlPcQCnTfWq7PRLLM6I3xHmDQVkquHHovLjvyH1lfIVeAjQFJy8GjV50Wu4kKRO2avVslJIRBeWqKRI0JT86zvwU0Y7gKiHUM1gwu/Zr4ZOGFslzjdPkt1wTS6OSNlyud74b+KlDniG3HoU8Rs8kpSbrRnOjBt+mweWC1o3lUjqntnjft8JmiRJoSnwGDXuwWUnQGqFNixDhTLkSNjMl8Ubso81jMFqLkpZj8xJlpYjKTy9w5PS6SzAUx2UIm9gkSc6FyWWw/v1pdF79obsPEllYazgAAANZ6VFh0U01JTEVTIHJka2l0IDIwMjQuMDMuMQAAeJxdjrFuAkEMRH8lJUFey2Pv3t2akoYqpI9SRKt0SS5ClHw8HCC8SWU/P3k06+0Xr/dbbmit6dvP7n2ZDU+nVRIeC2XWQptuTWCtVTIpL4Lsev1zTMompWS6RQTfESwyFBCuCZdPLWZ3CR5RoRensGK0UR5EJpJQ6eEev+lfclqin+njOH+/HuZfF17Wl/n4ydkRYK4B6hYAzwHiJaD4EDD4GDD6FDB5DaiOrgLE0ZUAHF0NqKMvYqcziWtu4hQ2o9IAAAAASUVORK5CYII=",
"text/plain": [
"<rdkit.Chem.rdchem.Mol at 0x73324e05f6f0>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"m2 = Chem.MolFromSmiles(f'{r1}.{r2}.{coresmi} |m:4:5.13,2:6.7|')\n",
"m2"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "b8ab2c01-71af-4680-9d42-b31d578cb5d3",
"metadata": {},
"outputs": [],
"source": [
"res = rdMolEnumerator.Enumerate(m2)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "cadc6622-2deb-4df0-86b7-04901f76f726",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABUYAAAEsCAIAAACAA/4ZAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOzdeVxN+f848Ndd2tOmVSKRJVvEiCyDECpri20YWcYWhpms09gLQ8Uge3YhtCgikV1pIUqbdi23vW63u5zfH+/P3G+/snfvPcnr+ZjHPM55n+X1OvMw13mf98agKAoQQgghRIeamhp3d/fbt2+rq6uvWrXKzs6O7owQQggh9D1h050AQggh9OPavHlzcXHx1atXMzMz7ezsEhISjI2N6U4KIYQQQt8NBrbSI4QQQnQRiUQURbFYLAAYMGDAunXrJkyYQHdSCCGEEPpuYCs9QgghRBsmk0k2eDxeRkZG586d6c0HIYQQQt8XJt0JIIQQQgj+/PPPoUOHduvWje5EEEIIIfQ9wSo9QgghJDuZdZkqcSomr0zEJRRFbd68+fHjx8ePH6cxMYQQQgh9j7BKjxBCCElSZl3mqpxV5m/M1ePV5WPljV4aTc+Y/qDqATkqAlGNqKZaVE12q6qqnJycnj9/Hh4erqamRl/WCCGEWqat77fqJOh4FnjSnQiSFqzSI4QQQhJznHO8S2KXPYV74rnx6iz1zgqdK0WV50vPD307dEPehsbnT5069cqVK3l5eSNHjuzXr9/Ro0dlnzNCCKHv2u3K23My53RJ7NI6obXBSwPLZMsd73dUCivJ0SphVbGgWPwpGbU8OD0eQgghJBnXyq65ZLoAgJOm0/Y2200UTACAT/H9S/1dc1wvlV5aqbuywSWXLl0SCATiXSUlJVkmjBBC6LsQGxsbHBysoKDg5OTUvn17cXmNqGbWu1kBZQEAwGaw9dn6AhA8rX76tPrpgaIDoZ1Ceyj1oC9rJCPYSo8QQghJAJ/iL85eDADOms7nO5wn9XkAkGPIzdCaEWka+bTr09bs1g2uatWqlWY9ioqKss4bIYRQ83blypVJkyYpKyvX1NT07ds3NzdXfGjmu5kBZQFKTKU9bfcU9SrK7pmd3zM/rXuak6ZTDj/Hr8SPxrSRzGArPUIIISQBweXB+fx8OYacV1svBjAaHO2p1JOWrBBCCH3vLCwswsPDTU1NAeDJkyc3btyYP38+AASWB14tu8oAxuUOl8epjxOfb6Jgcr7D+Tmt59io2dCWNJIhrNIjhBBCEnC/6j4ADFUdqienR3cuCCGEWg5jY2Pxdl1dnXguVd9iXwCYqDGxfn2eYAAD6/M/Dux4jxBCCElAUm0SAJgrm9OdCEIIoZbp4cOHycnJY8eOBQAhJYyqigIAJ00nuvNCNMMqPUIIISQBpcJSAGjNajhaHiGEEGq65ORkJycnPz8/0kqfw88hc9qbK+Gn5B8dVukRQgghCWg8fh4hhBCSiBs3bowbN+748ePW1takhHxHBgAdtg59eaFmAav0CElYAjdBO0F7UPIguhNBCMmUBksDADhCDt2JIIQQalH27dtnZ2c3fPjw2NhYT0/PmzdvQr3vyAwGflD+0eH0eAh9QEVFxYEDB5KSkjp16uTq6iqehoR4V/fuQNGBiMqIHH5OrahWj633c6ufl+ksIyt/CigBR8ApZZfSlDtCiB7dFLuFVYTF18TTnQhCCKEWpUuXLgcPHhTvKisrA4A6S53scgQcTZYmPZmh5gGr9Ag1RFHU6NGjzc3NnZ2dT58+7eTkFBoaKj56qPjQ8uzldVQdAOjJ6WmyNN/VvTtcfPgE58SetnuW6iylL3GEEJ2Gqg7dW7j3ftX9QkGhLlv3s+c/qX7inu9+wOhAR4WOMkgPIYTQd2r06NGNC9vKtVVhqlSLqhO4CZ0UOsk+K9R8YMd7hBpiMBjnzp07dOiQjY3NP//8Ex4eXldXRw6dKzm3KGtRHVU3XWt6SveU9z3fZ/TIKOld4t3Wm81gHyo6VC2qpjd5hBBdxquPN5QzrKPqVues/pLzjxQfuVVx6+/8v6WcF0IIoRaIzWAPVBkIAAFlAXTngmiGVXqEPsDExIRsFBcXa2pqysvLA0C1qHpp9lIAmNt67lnjs+IPoipMFVdd1zumdx51eaTCVKErZ4QQveQYcvuM9gHA6ZLTv7z7JZ+fLz7EEXD2Fu5t8NblbuCuwFA4V3IugZsg61wRQgh9/+ZpzwOAC6UXnlY/pTsXRCes0iP0USKRaP369UuWLCG7l0ovlQpLVZgqe9ruaXzyQJWBaiy1xuUIoR/HJI1JB4wOyDHkTpecbvuybY/XPUakjOiU2En/pf7vOb/vLdxb/+R28u3ma88XgWhT/ia6EkYIIfT9ctB0GN5quJASjksbd7nssghE4kOva1+vzl1dI6qhMT0kMziWHqEPoyjK1dW1rq5u3bp1pORe1T0AsFGzEc9HghBCDSzSWTRSbaRXoVd4Rfhb3tvE2kQFhkIvpV6jWo1aoL0AAOQZ8p0VOmuztQFgvf7645zjV8uuPqt+9pPKT3TnjhBC6HvCBOalDpcmpE94WPXQId1Bl63bQ6kHi8F6W/s2sy4TAPoo9ZmhNYPuNJHUYZUeoQ8oLS395Zdf1NXVAwICSK97AHhT+wYA+ij3oTU1hFBz11mh8wGjA2SbR/EUGAr1jxrKGSZ3Tybb+nL6S3WW7izY6Z7vHtoptOGNEEIIoUaOcY5RFEV63bdmt75nes+vxO90yekn1U8iKiMAoBWr1VDVoVM0pkzUmAgA6ix1QzlDdSa2SLVYDIqi6M4Boealurq6d+/e5eXlkydPJiXr169v165dl8Qub3lvD7U7tFB74Scuf1HzwiLJoqti1zdmb2SSL0LoO1YqLDV5ZVImLIvsHDlMdRjd6SCEEGrWwivCx6eN51P8512f91Pu1+BohbACAHAo6I8GW+kRakhOTu7QoUP1SzQ0NMTbDGDIPCOEUIulydJcobvi7/y/N+RtiOocRXc6CCGEAABycnL2799vYGCwfPlyUhIVFXXs2LFFixYNGDCArqxe1752zHDkU/x1+usa1+cBK/M/KqzSI9SQvLy8tbV143INtgbwgCPgyD4lhFALtkpv1YGiAw+qHtyquDVa7QOLDyOEEJKlkJAQDw8PLS2tt2/fkiq9q6trbm7umzdvRo8eTVeV/j3//bjUcWXCMgdNhy1tttCSA2qecMZ7hL5UF4UuABDPjac7EYRQi6LKVF2ttxoA1uWtowBHwyGEEM3Gjx8fFRXl6OgoLvHx8bly5Urnzp3pSokr4k5Mn5hZl/mTyk8n259kYiUO1YN/GhD9RCJRUFCQp6enuKS8vNzb2zssLIzGrBobqjoUAMIqwnBFEISQZC3TWWYoZxhTE3O97DrduSCEEGpeRCCa/m760+qnHeQ7BHUMUmYq050Ral6wSo9olp2d3b9/f29vb29vb1Jy/vz5IUOGnDx5srlV6R01HVuxWpULy//O/5vuXBBCLYoiU3Gt/loA2JC3of7CwgghhNDvOb9fK7umxdYK7RSqy9alOx3U7GCVHtHM0NDw4cOH9aejmzBhQkJCQv3OTvTiU/yoqigAUGOp7TTcCQC7CnatyllF5hQlSoWl3oXedyrv0JYlQug7N197vrG8cWJt4pXMK3TnghBCqLk4XHzYu9BbniF/qcOlLopd6E4HNUdYpUc0YzKZioqK9UuUlZtRb6JaUe3k9MkjU0YGlgcCwG/av7kbuDOAsadwj26C7sDkgWNTx3Z/3V0/QX9Fzop9hfvozhch9L2SZ8hvVdzad3Xf9aPWCwQCutNBCCFEv5vRN5dlLmMA43j74yNajaA7HdRMYZUeoY+qFlXbp9sHlwerMlX12Hqk8G+Dv591fTZda7omW/NJ9ZOwirAUXkoXxS6r9VbvabsHAOQZ8kbyRm3k2tCaO0Lo++PcwZn7npuSknLy5Em6c0EIoR9XRUXFqFGjPDw8Hj58OGrUqNevXwcHB48aNerRo0c7duyYNm2abNKIjY2dOnxq73W9t2pvnaE1QzZB0fcIF7FD6MPKhGXjUsc9rn6sL6d/q9Otnko9xYf6Kfc7a3wWAHgUr0ZUo8nSrH9hD6UeWT2yZJ0uQuj7x2Kx3N3dnZ2dt2zZMmvWLAUFBbozQgihH5GysrKHh4d4t3379nw+39DQUEVFZezYsf3795dBDnl5efb29lVVVR21Oq41WiuDiOj7hVV6hD6ggF8wJnVMPDe+vXz726a3Oyl0+uBpCgwFBRa+cyOEJMbR0dHDwyMuLs7X19fV1ZXudBBC6EfEZrMtLCzEu8XFxVZWVtXV1QBw/fp1NpvduXNnCwuL7t27m5mZDRw4UFtbW7IJVFZWjhs3LicnZ+jQoSdPnmQwGJK9P2phGBSFS+Aimm3cuDE1NTUwMHDmzJkODg5dunTZunVrXFxcTU3NoEGDPD09NTQ0ZJlPVl3WqJRRb3lvuyp2DTcNbyvXVpbREUI/uKCgIHt7ex0dnfT0dFVVVbrTQQihH5pQKLSzswsNDWWxWAYGBkpKSmlpaSLR/y1NwmAwTExM+vTpY25uTv7dpk2TRl8KhcJJkyYFBQV17dr10aNHmpqan78G/diwSo/oFxoaWlVVRbZ79+6tpaV19OhRExMT8knS1tZWSUlJZslk1GVYp1in89L7KvcN6xSmw9aRWWiEECIGDRr0+PFjDw8PNzc3unNBCKEf2ooVK7y9vbW0tAICAoYNGwYAVVVVycnJiYmJMTExMTExsbGxNTU19S/R0NDo3r27xX+6du3KYrG+POKSJUsOHDigra39+PHjTp0+3FEUofqwSo+al3v37o0dO5bL5bJYrGnTptnb29vY2LRq1Uo20RNrE0enjM7j5/2k8lNox1AttpZs4iKEUH137tyxtrbW0NBIT0/H9hmEEKLL4cOHFy5cKC8vHxoaOmLEhyecFwgEycnJr1+/JpX8p0+fFhUV1T9BXl6+U6dO4hp+3759P7G6065du/78809FRcWIiIiBAwdK+HlQC4VVetSMJCYmDh48uKysrH6hkpLSyJEjbW1tbW1tDQ0NpRf9WeyzyfzJuezcMWpjAkwClJnNaC09hNCPZuTIkREREX/99demTZtkHLqurq66uho/JSCEfnChoaH29vZCodDPz2/WrFlffmFeXh5pwCf1/Ddv3tSvcJGh+GQcvoWFhaWlpY7O//qEXrlyxdHRkaKos2fPymxefdQCYJUeNRcvXrywsbEpKioyMzO7c+dOZmZmZGTk7du3IyMjxUs0m5mZ2dnZ2draDho0iMmU5BKM9+7ds7Oz69S3k+kh01NdTikwcNI7hJAUkQlE5s2bp6amBgAURd28eZPD4cyY8b9lip49e2ZpaamiopKWlqarqyuzxP7+++/Dhw9ra2vzeLxLly716tVLZqERQqj5iI2NHTp0aFVV1ebNmzdu3EgKd+7c6eDg0KFDh6+6FYfDiY2NjY2NjYuLi4uLS05OFgqF9U8wNjbu06ePrq7uqVOnuFzuzp07//jjD4k9CfoBYJUeNQtRUVF2dnbl5eXjx4+/dOlS/cHzxcXFd+/eDQoKun79ekVFBSnU0dGxsbGxs7MbM2YMeSFuipCQEAcHBy6XO336dD8/PzYbV4JACEmLSCRycnLKyMhITExMTEw0MTFJTEycPXu2SCRisVjPnz8Xn2lraxsSErJ69epdu3bJJrfi4uJNmzbt2bNHTk5u69atL168CAgIkE1ohBBqPvLy8gYMGJCTkzNt2rSzZ8+S2Z28vb1XrFjRrl275ORkRUXFr7qhQCA4cuTIkydP9PX158+fz+PxYv4TFxdH5tIHAE1NTWtra39/f8k/EmrRsEqP6BcRETFhwoSqqipnZ+dTp07Jycl98DSBQPDkyZPg4OBr164lJyeTQkVFxcGDB9va2k6aNKldu3bfEP3ChQu//PILn89ftGjR/v37Jdv4jxBCjRUVFeno6Ojr6z969MjExITL5QqFwpiYmNWrV9ev0r98+dLc3FxeXj4lJaVtW1kvvbF+/fq8vLwTJ07IOC5CCNGrsrJyyJAh8fHxQ4cOvXXrloKCAgCEhIRMmDBBJBKdOnVq5syZX3vPBQsWAMCUKVOio6MPHjyYnp4uLy9PDgmFwuTk5NjY2J07dyYkJGzYsGHLli2SfSLU8lEI0er69evkS+eCBQuEQiEpzMjIqKur+8RVaWlpXl5e1tbW9VvUzczM3NzcoqKixPf5LF9fX1KHd3Nza+qTIITQ19DT00tLSxPvRkZG9uvXr8E5Dg4OADB//nw+ny+zxHx9fYcMGWJmZlZUVCSzoAgh1BwIBAI7OzsA6Nq1a0lJCSmMiYlRUVEBgK1bt37bbd+/fy8Sicj92Wz2u3fvGp9z+/ZtAOjYsSM5E6Evh630iE5nzpz59ddfBQLB0qVLfXx8SL+m2NhYGxuboUOHXrhw4bNrfnA4nIiIiKCgoMDAwPLyclKora09duzYz3bL379/v6urK0VRf//9t7u7uwSfCyGEPkvcSk92792716CVHgCSkpK6d++uqqpKhh0pKipqamoqKSl924aamtqXLKRUXFz87t27zZs3m5qa/vPPP9J4doQQap4aLyCXm5traWmZk5Pz66+/Hj9+/JvvXFNT8+LFizNnzhQWFn5wTJNIJDIyMsrLy3v+/Hm/fv2+/RnQjwfHDCPaHDhwYNmyZSKRyM3NzcPDgxQ+evRo/PjxZWVlNTU1fD7/s2+frVu3dnBwcHBwEHfLv379elJS0unTp0+fPs1mswcMGGBnZzdp0qTOnTvXv9DT03PNmjUMBsPLy2v58uXSekiEEGqCmJgY0lzDYrGEQmFtbW1+fn5Tbqiurq6kpKSsrKyhoaGkpKSkpMRkMlNTU3/++WcfHx/SDKWtra2trb127VoHBwes0iOEfhy7du06cOCAoqJiYGAgqc9XVlaOGzcuJydn2LBhhw4dasrNi4uLT5w48eLFC9ILQKy8vLyqqsrQ0JDJZDo4OHh7e58/fx6r9OirYCs9ooe4Rr179+7ff/+dFEZGRtrb21dWVk6YMOHixYtk8NI3SE9PDwoKCg4OvnfvHp/PJ4UmJia2trZ2dnbDhg1bs2bNnj17WCzW0aNH58yZI5EnQgihr/LZVnoOh2NmZlZYWHjixAnyS8XlcktLS2tra79to7y8XCQSfSwfCwuLq1evrlmzZufOnbq6umvWrImPjycdQRFC6Hvx6tUrZWVl8U9rVVXV06dPWSzWgAED6s++3FjjBeQEAsH48eNv3brVrVu3hw8fSmRpTx6PZ2JicuXKFUtLSwC4ePHinDlznJycTp48CQBPnz61tLRs06ZNVlbWl3SqQuh/aO32j75LKSkpS5YsmTNnjrgkNDTUxsamb9++f/zxR21t7acvF4lEZGUOUqMWlwcFBZFB9TNnzpTUqNGioiI/Pz8HB4f63e9btWoFAAoKClevXpVIFIQQ+ir+/v5ubm4qKioLFy48c+YMRVEbNmyYPn26vr6+m5tbfHw8OY0sg/zzzz9LcFxlSUlJbm5uWlpadHR0VFTUrVu3Dh06RJqMWCxWZmbm7t27e/ToYWhoOGnSpKysLEnFRQghaXv06NGIESNatWq1dOlSUpKXl2dkZOTg4GBvb29qasrhcD527fPnz5WVlQFg586d4sJFixYBgLa2dkpKShNzi46O5vF4FEUJBIIOHTqEh4eT8tTUVAaDoaamVlNTQ0pMTU0B4O7du02MiH4o2EqPvk5CQsLatWs7d+58586dhIQEAMjPz+/WrdvZs2fNzc1/+eWXIUOG/P333x+7nKKoFStW+Pj4yMvLnz592tHRkZSfO3duzpw5fD5/8eLF+/btk/i080Kh8PHjx8HBwYGBgRkZGWw2+9ixY+LoCCEkS5GRkW/fviXbJiYm1tbWx48fFwgEpMTa2trExCQiIsLa2lpRUTEhIYH0/5QqsmDewoULm9izFCGE6PLmzRs5ObkrV67k5OTs27cPANauXZufn08awCdNmmRpaenm5tb4wnfv3llaWhYUFLi4uBw9epQUkv6kSkpKERERpEW9KVasWBEZGTlixIjY2FiBQBARESFe4Kl///7R0dGXL1+eMmUKAGzYsGHbtm0LFizw9fVtYlD0A6H7mwKSOpFIdPr06V9//dXd3V1S0xc/efKkZ8+eZPuff/6ZOnUq2b5//367du0+dpVAIJg9ezYAKCgoXLt2TVx+8OBBWU47P2zYMAC4fPmyDGIhhNA3qK6u7tixIwB4enqSkqKiomfPnkkv4tu3b9lstpycXGpqqvSiIISQtHl4eIhb6S0tLcVdMk+cOGFjY9P4/PLy8h49egDAmDFjxL1EL126xGQymUzmlStXJJVYfHz8+fPn79y506AvKpmyZMqUKWT39evXAKCpqUla9RH6ErgEd8u3ffv2HTt2jB07lsfjDRkyRNwQJCmpqandu3cn22ZmZtnZ2Twer/FpPB7P0dHRz89PVVWVrO1Jyj09PRcvXkxR1K5du8ST5EnViBEjAODBgwcyiIUQQt9g48aNaWlpvXr1WrlyJSlZuXLlwIEDjxw5ItlAQqGwpqYGAExNTWfNmsXn8zdt2iTZEAghRJeCggIdHR2yraOj8/79+wYn8Pn8KVOmvHr1qnv37hcvXiRLIz9//nz27NkikWjnzp2TJ0+WVDK9evVydnYeMWJE/QWYAWDatGksFiskJISs3NStW7devXqVlpbevHlTUqFRi4dV+pbv+PHju3fvdnBw2LFjh6KiosR/IOrPSy8nJ0dRVOOvBtXV1XZ2dgEBAZqamuHh4SNHjiTlpFMTAHh7e69evVqyiXG53Lq6usblQ4YMAYCoqCjJhkMIIYmIi4vz8fFhsVjHjx8nPTMjIiLOnj0rLy8/fPhwCQa6f/9+jx49xHV4d3d3BQWFc+fOkTYihBCSjZKSEok3OBEaGhpVVVVku6KiQktLq/5RiqLmz59/+/ZtAwODGzduqKurA0BGRoadnV1NTc28efNWrVoljawaMDAwGDJkSG1t7bVr10gJmZzv/PnzMoiOWgas0rd88vLy4pqtkZGReACnpBgaGmZnZ5PtrKwsLS0tsgySWFlZ2ejRo8PDw/X19SMjI8l4JIqiVq5cuWbNGhaLdeLEiWXLlkkwJQ6HY2tr26VLlw4dOkyfPr20tNTOzq5r164URQGApaWlvLx8XFxcZWWlBIMihFDTCQQCFxcXgUCwcuVKCwsLAKipqZk/fz5FUZs2bZLsoHoVFZXk5OR9+/bl5uYCQPv27efNmycUCt3d3SUYBSGEPqa4uHj48OG9evUyMjLasGGDxO/fs2fPp0+fku3nz5/37Nmz/tEtW7b4+fkpKytfu3atXbt2AFBeXm5vb19QUGBjY3Pw4EGJ5/MxDerw06dPZzAY169fF3+PQOgz6Oz1j2Ti1KlTenp68+bNGz58uKGhoXhkZlPUH0sfGxurq6tbUFBAUdTy5ctdXFzqn1lQUGBubg4A7du3f/v2LSkUCARz584FAAUFhYCAgKbn08C1a9e8vb0piqqrq+vVq5e/v3+bNm0AICkpiZwwYMAAALh165bEQyOEUFOQ8Uft27evrKwkJWSZz969e9fV1Uk8HOlTWn92aGVlZQaDIdVx+wghRCxatGjKlCkikaisrMzY2DgqKqopdysrK/P3958xY4aNjY2/v391dfXz58+1tLROnDjh6+urpaUlfg+kKOrChQsMBoPJZIpnd6qrqyPdSLt3715WVtakB/tKJSUl8vLyLBbr/fv3pGTQoEEAcO7cOVmmgb5f2Erf8s2aNevhw4eTJk06c+ZMp06dyGfIphgxYsScOXNSU1P79et348YNc3NzFxeXPn369OrV6969e/Wnu8/KyhoyZEhcXFzXrl0fPHhAluWoq6tzdnY+fvy4iopKUFDQpEmTmphPYxMmTHB1dQUABoPB5/NVVVWtrKyg3vh57HuPEGqGMjIytmzZAgCHDx9WVVUFgOjoaG9vbxaLdezYMfH0yBK0fft2Fot1+PDhjIwMADAwMFiyZAlFUZ9YuAQhhCTl1atX9vb2DAZDXV3d3t4+LCysKXfjcrkxMTFt27bt3bt3TEwMn8/v16/ftWvXnjx5Eh8ff/PmzS5dupAzHzx4MHv2bIqivLy8yOxOFEXNmzfvzp07bdq0CQ0NJZ3wZUZTU3PMmDFCofDy5cukBPveo69D7xcFJEtJSUkqKip5eXlNvE9paWnJf8SzcfJ4vMLCwvqnpaenm5iYAECfPn3Eh6qrq21sbABAQ0Pj4cOHTczk02bPnt2lS5clS5aIRCIfHx8AID/fFEWR0UrDhw+XagIIIfTlRCKRtbU1APzyyy+khM/n9+3bFwBWr14tvbgzZ84EgLlz55Ld4uJiNTU1ALh37570giKEEEVRK1assLOzKyoqSkhIsLKyEv8QSVVqaiqZM2/58uXiQvIdU1VVNTY2VgY5NHbu3DkAsLKyIruFhYVkFZLi4mJa8kHfF6zSt3wXL16cO3fuggUL9PX1Dx8+LJugiYmJpK/7Tz/9xOFwSGFZWdngwYMBQE9PTwa/mG/fvr127Vr79u0fPXr04sULAOjYsSM5xOFwmEymkpISLhCCEGomjh07BgDa2trib6A7duwAAGNjY3EnfGnIyMggHT7fvHlDSshY+iFDhkgvKEIIURRVUVExf/58MzOzqVOnTp48WaqfL4ni4uLOnTsDwPjx4wUCASk8f/48g8FgsVj1l1iWserqalVVVQaDkZGRQUpGjx4NADJ7dUffNazSt3xcLjc0NNTf3z8kJKS0tFQGEZ8/f66trU2awSsqKkghh8MhI9jbtWuXnJwsgzQIV1fXVatWCQQC0ocqNzeXlJuZmQHA48ePZZYJQgh9TFFREfnZPHv2LCl5+/atkpISANy8eVPa0RcsWAAAzs7OZLeyspI0YYWHh0s7NEIIESNGjDhy5IhUQ/B4PLJuSN++fauqqkjh/fv3FRQUAGDfvn1Sjf5ZpLP9jh07yO6JEyewSyn6Qlil/yG8e/eubdu2ZKiFoqKigYGBmZmZtbX1rFmzXJJA1HoAACAASURBVF1d3d3dvby8/P39o6KiXr16lZubKxKJvjlWZGQk6bRpZ2fH5XJJYX5+PplltEOHDmlpaRJ6rI/y8/Pz9vbmcrkcDsfc3PzQoUMURZEO/xcvXiTnLFy4EAB27twp7WQQQuiznJycAMDGxobsijvhi4cLSVVubq6SkhKDwRD3n/L09ASAfv36NeWvA4QQ+rSamprs7OyqqqqDBw8aGBhIteVJJBKRcUZkqSZSKO6Ev3LlSumF/kKBgYEA0KtXL7JbXl6upKTEZDJzcnLoTQw1fwyKoqQxRB81K5MnT7569SoAsFgsoVD42fMVFRW1tbW1tbX19PS0/6Orq6urq6tdT+MLQ0JCHBwcuFzu9OnTT548SSZzevfu3ahRo1JTU83MzMLDw0mHfKkqKCj4/fff79+/z2aznZyctm3bxmKxtm3btmHDhmXLlpFx9WfOnJk1a5a9vf3169elnQ9C6MdRVFRUW1trZGREdkUiUVJSUm1tbc+ePT82v11oaOi4ceOUlZVfvnxJpiA5fvy4i4uLtrb269evyeumtK1cudLLy0v8k8jlck1NTXNzc69fv25vby+DBL5QdXW1n5+fgYGBNKZWRQjJWHR09KJFizgcTp8+fbZv3y6evu5LCIXCiooKABCJROXl5QBAUVRZWRk5Wlpa2mDj4sWLV65cUVNTe/DgAWln4nA4AwcOTElJGT9+/PXr11kslgQf7Rvw+XwDAwMOh/Pq1avu3bsDwNSpU69cubJnz56VK1fSmxtq5rBK3/JRFNW5c+fU1NRNmzb99ddflZWVBQUFRUVFxcXFxcXFRUVFhYWFxf8hh6qrqz97WxaL1aC2z+Fwrly5wufzFy9evG/fPiaTCQBJSUmjRo3Kycnp169fWFhY69atpf/EH3b//v1hw4aZm5vHxsYCQGZmprGxsZaWVlFREUkVIYSaIj09ffv27RcuXBg1ahT5ilpbWzty5EihUKihoZGSkhIVFdX4m2Z1dXXPnj0zMjL27t27YsUKACgoKDAzMyspKTl79uz06dNlk3xxcbGJiUllZeXjx48tLS0BYP/+/cuWLevRo0d8fHwz+ZEMCgpat26doqJit27dTp06RXc6CCHJqKur+/XXX+/fv9+pUyc2m01KyLson88na7MLBILKysr6G19LR0envLw8KCiIjFEnv8+PHj2ysLC4d++eioqKJB/pWy1YsODIkSMbNmwgq59cuXJl6tSp/fv3f/bsGd2poWYNq/Qt361bt8aMGdOuXbv09PQv/ABZW1tbUlJSWlpaWlqan5+fl5dX+h/xbmFhYeMG/1atWrm4uOzdu5fsvnjxwsbGpqioaNiwYYGBgaRDPl14PJ6GhkZdXR2Hw9HQ0ACAdu3aZWdniz+FIoRQUyQlJeXl5WVnZ1+7do1U6Y8dO3bixIl79+6xWKy5c+dqaWnt3r27wVXLli3bv39///79Hz9+TH6inZyc/P39bWxsQkNDZZn/+vXrt2/fPmrUqFu3bgFAXV1d165dMzIyzp8/7+zsLMtMPobL5SoqKvr6+j569Air9Ai1GBMmTCB9zr8Kk8kk0yQxGAzyXgcAmpqaDTY0NDQYDAZFUTdv3qysrHz+/Hm/fv0AoLq6esaMGXFxcU+ePNHX15fMkzTZ3bt3R4wYYWJikpqaymAwamtr9fX1y8vLk5OTyax+CH0Ynb3+kUyQPpPbt2+X7G35fH5eXl5CQkJERMT58+f37dtHOpqGhYWRE+rq6tq3bw8Atra24kH19CKr09+4cYPskmlIyEh7hBCSiDNnzkycOJFsOzk5eXl5ke2wsDDxCEmxJ0+esFgsNpv94sULUhISEgIAysrKMph2pIHS0lItLS0AiIiIICVHjx4FAFNTUz6f/+lreTze4cOHQ0JCpJ8mdfDgwVmzZskgEEJINkj3JUtLy+PHj4eHh4eHh0dGRkZHR0dHR8fFxaWlpaWlpWVmZpLlk8vLy78tyvLlywHg999/F5cIBALxoPpmQigUGhoaAsCTJ09IyezZswFg8+bN9CaGmrlm0ZUOSU9WVlZISIi8vPzcuXNJSU5OjkTuzGazDQwMevbsOXz4cGdn56VLl5IaclRUFDlBTk7uwoULLi4uAQEBioqKEgnaRGQJvQcPHpDdIUOGQL2EEUJIst6/f6+np0e2DQwM8vPz6x8VCAQLFy4UCoV//vlnnz59AKCysvK3334DgG3btpFB9bKkoaHx+++/A8DGjRtJyZw5c7p27ZqSkvLZJvEdO3Z4enp6e3tLPUuEUMvC4XAKCgoUFBSCgoJ+/fVXa2tra2vrYcOGWVhYWFhY9O7d28TExMTEpF27dpqampqamt/c5ZO8pl64cEHcyZTFYolnj24mmEymo6MjAJw/f56UkLTJqvUIfQxW6Vu4Q4cOCYVCR0dH8lqZnp5ubGxsa2tLSWHARYMKMwBYWloePXr0YzNCyV6DOnzjhBFCSIJatWpVU1NDtquqqkgfUbEdO3bEx8ebmppu2LCBlKxbty47O7t///7Lli2Tda4AALBixQo9Pb2HDx+SPv8sFuuvv/4CgE2bNvF4vI9dlZycfPLkyfXr18suUYRQSxESEiIUCn/++ecPzrssQQMGDOjUqVNeXl4zb8shdfiLFy+STw8jR47U09NLSkqKi4ujOzXUfGGVviWrq6s7fvw4ACxZsoSUHDx4UCgU6ujoMBgMiYcbPHgwk8l8+vTpJ9786GVlZcVkMp89e1ZbWwsAPXr00NLSyszMzMrKojs1hFAL1K1btxcvXpDt2NjYbt26iQ+9fft2+/btDAbj4MGDZP35p0+fHjx4kM1m+/r60jXxsoqKipubGwBs2LCBfPl1cnLq3bt3VlbWkSNHPniJSCSaN2/ezp07xWNZEULoywUFBQGAnZ0d2ZVGm5MYWS5U3ADePPXv379z587v37+PjIwEADabPXXqVGj2aSN6YZW+JfP39y8oKDA3NyfTF3O53BMnTgDA4sWLpRFOU1PTzMystrY2JiZGGvdvOg0Nje7du/N4PJIhg8EYNGgQNGqoLygo6NixI05DghD6Knw+//bt24mJiYWFhbdv3y4vL583b965c+fOnz9/69atHTt2iL+uUhQ1b9682tpaFxeXkSNHAkBdXZ2Li0v9Tvh0WbRokZGR0YsXLwICAgCAyWRu2rQJALZt2ybucVDfv//+q6qq6uDgIIPcOBzOwoULz5079+TJk4ULF+bm5sogKEJIevh8fnh4OACMHz+elIwePXr8+PHv3r2TRrgZM2YAwKVLl+rq6qRxf0khM5I26Ht//vx5kUhEZ1qoGcMqfUv277//AoC4A+eFCxc4HM5PP/3Uv39/KUVs/qPTv6Tvvaur6+TJk8VLmyKE0JcgU8SlpqYaGhoePnz4/fv3nTt3vnr1amBg4KFDh7y8vMaMGUPO9PX1jYqK0tfX37lzJynx8PBITEys3wmfLoqKiqQL/fr16wUCAQBMmDDB0tLy/fv3+/fvb3ByVlbWxo0bXVxcYmJi0tLSKioq0tPTpZebkpKStbX1smXLtm3bZm1traqqKr1YCCEZiIyMLC8v7927t7GxMQBwOJy7d+/euXNHSmsed+vWrVevXqWlpWRdj2aL1OGvXLlC+r0OGjTI2Ng4Ozv74cOHdKeGmit6Z+dD0kNWX9fQ0KiuriYlZNEOPz8/6QU9e/YsAJCx+s0TmV9k/PjxZJf8OPbo0UN8QnBwcL9+/dLS0nR0dGjKESHUkuXl5ZE+6pcuXSIlycnJioqKDAbj9u3b9OZG1NXVdezYEQBOnTpFSsjrr6amZmlpaf0z79275/CfQYMG6enp7dmzh46UEULfJdLsREb6UBTl5+cHADY2NtKLuGPHDgCYPn269EJIhLm5OQBcu3aN7JIhUYsWLaI3K9RsYZW+xSJT3IvX6njy5AkAaGtrS3U9uezsbPLaJxQKpRelKciE/+rq6gKBgKIoHo+nrKzMYDCKi4spiiovLzc2No6JiXn37h1W6RFC0jBp0iQAGDduHNkVCoWk99C8efPoTaw+8mJtbGzM4/FIyfDhwwHA3d39Y5cEBASMHj26qYH5fCoqitq9m1qzhtq4kfr3X+r166beEyHUXHXo0AEAnj59SnbJoPF///1XehEzMzMZDIaysnJlZaX0ojSdp6cnADg5OZFdMjeetrZ2XV0dvYmh5gmr9C1TaWkpqakmJyeTklmzZgGAm5ubtEOTtegTEhKkHeibkc5d8fHxZHfYsGEAEBgYSFHUb7/9Rj6CYJUeoR/B7du39+7de+vWLZlFvHLlCgCoqallZWWRkoMHDwKAvr5+SUmJzNL4LIFAYGZmBgCHDh0iJWSAkqqqakFBwQcvSU5OvnjxYpOiBgZSHTtSAA3/GT4cK/YItTwvX74EAF1dXdIOxOPxyAJ1GRkZUo1L5lE6f/68VKM0UVZWFpPJrP/poXv37gBw48aN+qfduHGjU6dOMni9R80cjqVvmY4dO1ZTUzNmzBgyx1txcfGlS5eYTObChQulHbr5rwz3seH07969O3z4cFhYWL9+/ezs7EpLS0ljGkKoRVq7du3ixYsrKir++uuvVatWySBieXm5q6srAOzYscPIyAgA8vPz165dCwD79u3T1NSUQQ5fiMVikVnxNm/ezOVyAcDKymrs2LFVVVW7d+/+4CWdO3cmyyl/o2PHYOJESEuDESPgxAmIioKICNi9Gzp2hLt3YfBgiI399psjhJof8Vz3TCYTACIjIysqKsTj6qVHPNucVKM0kZGRkZWVlZqaWlJSEikZN24cAPzxxx/ic6qrq5cvXz5s2LCioiJ6skTNB93fFJDkiUQiUpMPCgoiJWTgkJ2dnQyiHzp0CACmTZsmg1jfxtfXFwCcnZ3JblhYGAAMHDiQx+Ol/ef+/ftaWlriZjSEUAvD5XIVFBRevXpFURSHw/lE43NT5ObmhoeHe3l5LViwwMrKSklJSUVFxdzcXDw06c8//wQAe3t7iYduOpFIROZSFQ+Pj4mJYTAYioqKOTk5Eg726hWlqEgBUDt3NjzE5VI2NhQA1aULVVsr4bgIIfoMHDgQ6g0XX7p0KdQbVy89hYWFbDZbTk6ODLpstvLy8sgoUSI4OBgAmEymuHDp0qVubm4+Pj5z586lKUfUXGCVvgW6ceMGALRv3578Py8UCslQpbCwMBlEf/XqFQAYGhrKINa3ef36df0My8vLWSyWnJyceB5BiqKysrJMTU1pShAhJHXl5eUMBqOoqIjsdu3aNTIysik3FIlE6enpwcHBnp6es2fP7tevn4qKSuPP6AwGY8CAAeKrBALB3r17s7Ozm/QwUhMSEgIA2traFRUVpGTKlCkAsHTpUglHmj2bAqD+m1+gIQ6Hat2aAqD+m64PIfS9KygoYDKZCgoK4o7lDcbVS9Xo0aMB4PDhwzKIJREikah3794AoKmpSUqePHnSoUOHqqoqrNIjCjvet0hk7brFixezWCwACAkJycjI6Nix46hRo2QQ3czMrHXr1rm5uZmZmTII9w26du2qo6OTm5tLVj1VU1Pr1asXn89//vy5+BwjI6O3b9/KPrd3795VVVXJPi5CPxo1NbUpU6Y4Ozv7+Pg4OztnZmZWVlZ+1R3y8vJu377t7e29cOHCwYMHq6mpmZiY2Nraurm5+fn5RUdHV1dXa2pqWllZLViwwMvLKzw8PDU1VVtb++nTp2QdZgBgsVgrVqxo27atFB5RAsaNGzd06NDi4mIfHx9Ssm3bNjabffjw4RcvXnxwmfpvIRRCQAAAwJIlHz5BSwtmzQIA8PeXTESEEN1CQkJEItHIkSPJapQvX77MyMjQ1dUlyzNJ23fR976+v//+Oz4+Xk5OjnSGraurc3Fx8fHx+eC3Y/QDYtOdAJKwzMzMsLAwBQWFOXPmkBJSw1+yZAkZqiRtDAbDysoqMDAwKiqKTJXX3DAYjEGDBl2/fj0qKoqM1xoyZEhsbGxUVBSZKo8WkZGRO3bsuHPnzsGDB+fPn09XGgj9OM6dO3fhwoWioqK1a9dGRUV9ul6dl5f3+vXrxMRE8u/4+PjGX980NTXNzMy6d+9O/t2zZ089Pb0G56xevdrNzW3dunXW1tYMBkPCjyQFmzZtWrRoUY8ePchuly5dRo4cGRUVZWFhQUrU1dWVlJSUlZUbb2wyNNQXCkFZGTQ0QFGx4YaSEigpgYYGvH4NlZXAYICV1UfzGDoUvLzg2TMZPDJCSAbEA+kb7MrmZXXy5MmLFy++d+9ebm6uoaGhDCI2xYULF7Zs2cJisS5fvmxvbw8Ax48fT01N9fb29vb2zs7OrqysPHTo0G+//UZ3pog2WKVvaQ4cOCAUCmfOnKmrqwsAaWlp4eHhSkpKs2fPllkOgwcPDgwMfPDgwcyZM2UW9KsMGTLk+vXrDx48IAsBDB482MfHh94p/Xg8nqen565du2jMAaEfipycHPkFiI6O5nK5Xbt2FR+SVAW+sWXLlvn4+ERHRwcFBZE3s2bu559/TkxMFL9kc7nclJQULpfbqlUrgUDA5XLLy8vLy8s/eK2bkRFkZ38mwIIFMGUKAIC6Oqirf/Q08oG4qAgEAmDjqwtC3zcej3f79m0Gg2Fra0tKGtTwpU1NTW3s2LEBAQH+/v4rV66UTdBvc//+/Tlz5lAU5e3tLf5bw9nZmYwdAIBTp069fv3awcGBvhwR/fDvxRaFx+OdPHkSAJb8133x33//FYlE06dP19LSklkaDaaUb4bILPfiDEnCDx8+FAgEbJpeFseMGUNLXIR+WP7+/mVlZWVlZT4+Pv/88w+TyRw9evTDhw8piiJzvNfXpk2b+hV4MzOzb5ugXklJac2aNcuWLVu/fr2tra1s2qOaqH6SmzdvTk9P79GjR0xMjLy8PACUlpbW1tZyudzGG7pVVVBRAbW1UFoKXO4HNrhcUFaG2loAACWlTyWhrAwAQFHA5UKrVtJ8XISQ1N29e7eysrJPnz6ke1RhYeGzZ88UFBRGjhwpsxymTZsWEBBw/vz55lylT0pKmjhxIo/HW7Vq1ZJ6Q5M0NDQ0NDTItqGhIYfDad26NU05omYBq/QtyoULFwoLC/v27UumKeZyuX5+fgCwaNEiWaZhYWGhrKz85s2boqIiHR0dWYb+Qn379lVRUeHz+TU1NcrKyvr6+p06dUpNTU1ISOjbty/d2SGEZKFjx46hoaFMJjM4ONjc3Hz69OniIe7f1gL/hRYsWLBnz55Xr175+/s7OztL6rYykJCQQL59+Pr6kvo8AEhg7b2ICACA6upPnUNmOmCxAAeOIvT9a9Am32BcvWzY2tqqq6s/f/48JSXF1NRUZnG/XHFxMVlT2dbW1tPT82On4WhNBADfQfsA+nJk2PyyZcvI7tmzZ0tKSgYOHCge9CgbcnJyZErnR48eyTLul5OTk8vOzk5JSVEmLT//tdufO3eOoihaU0MIyYiFhcWGDRvWrVtnbm4OAAkJCQBgZWWVm5tbUlLy4MEDX1/f5cuXW1tbS7A+DwDy8vIbNmwAgL/++ksgEEjwzlIlEol+++03Pp+/dOnSQYMGSfLWZBaDigr4xNLKaWn/O/N76NeAEPo0sjZT44H0ssxBUVFxwoQJAHDhwgVZxv1CtbW1EyZMSE1NtbCwuHDhApnxGqGPwb8aW44XL148f/5cU1PT0dGRlJBZMZd8bA5haSI1ZHpHp39a/Zalurq6t2/fGhoa/vPPP3p6er/88sulS5e+dvprhNB3TV1dHQA2btzYpk0baceaPXt2165dU1JSTp06Je1YkuLl5fX48eN27dpt3bpVwrc2NQXSZfTevY+ec/cuAIClpYRDI4RkLj4+/t27dwYGBqTBqfG4epkh896fO3dOxnE/i6IoFxeXR48eGRsbh4SE4LT26LOwSt9ykEWG5s2bR1qeHz16FBMTo6OjQ5YRlrHmP5xejMvlTpw48dGjR6WlpW3atCkqKjp9+rSjo6Oenp6dnZ2vr29ubq4M0uBwODExMSUlJVlZWXFxcTKIiBCqLzk5GQDIwr/SxmKx3N3dAWDTpk08Hk8GEZsoMzOTJLxv375WEh/KzmDAtGkAAD4+8MF+Uu/fA2lGmzFDwqERQjInbpMnq36QcfXm5uayX86TdMJKSkpqbu9da9asOXfunJqaWmBgoGS7iaGWCqv0LQSHw7l48SKTyRSvYEE64c+fP19RUVH2+QwcOJDNZkty4WLpqKqqsrOzCw0N1dXVffjwYW5u7qtXrzw8PKysrHg8XnBw8G+//da2bdvu3buvWbPmwYMHIpHo6wJQFISGwtKlMHo09OsHY8aAqyvcutX4tTUxMdHT07NVq1bJycne3t4Se0KE0BcoKSnhcDitWrXS19cHAC6Xe+LEiWfSXDLNycmpd+/eWVlZR44ckV4USVm6dGlVVdW0adOkNUv/qlXQqhVERcH69Q1/HisrwdkZKipgwAAYP14q0RFCMvSx5etknwmbzSbtXs1qgfpjx47t3LlTTk7u8uXLPXv2pDsd9J2gUItAps0YP3482S0sLFRQUGCxWBkZGXSlRPpTRURE0JXAZ5WUlFhaWgKAgYHBy5cvGxx9//79sWPHJk2aVH+ylmVjx1Lz51PXr1PV1Z8P8O4d1b8/BfC/fxiM/9seNIjKyZHKUyGEvt7Dhw8BoF+/fmSXtNiYmZlJNei1a9cAQF9fv/pLfk/oc/r0aQDQ0tJ6//69FMNcukTJyVEA1E8/Ufv3U2FhVFAQtXkzZWhIAVBt21IpKVKMjhCSiYKCAiaTqaSkJP7dMzY2BoDnz5/Tkg/pT2pkZCQSiWhJoIGwsDCy9NKhQ4fozgV9T7BK3xIIhcIOHToAwI0bN0gJGes4ceJEGrNasWIFAGzatInGHD4hPz+/V69eANChQ4fU1NRPnMnn86Oiotzc3Dp37hwwYMD/6uSKipS1NeXlRWVlffiyggKqbVsKgDIzoy5dosrKKIqiSkupc+eoTp0oAKp9e6q4WApPhhD6aidOnACAGTNmkN2LFy/K5ieUfFX09PSUdqBvVlxcrKurCwAnTpyQerB796hevf7v06f4Y+iUKVR2ttSjI4Sk7+jRowBgZ2dHdsn3UwMDA7pq1CKRiHxTiIqKoiWB+hITE8nSdOvXr6c7F/SdwY73LcGNGzcyMjI6duxI1jYXCoWkJyctE+OJNecZ8jIzM4cOHZqQkNCtW7eoqKiOHTt+4mQ2mz148GAPD4/k5GT7w4dh2zawtIS6Orh9G1asgHbtoE8f2LgRnj2D+t3yly2DnBywsIAnT2DqVFBXBwDQ0IBp0+DJE+jWDTIz4fffpfygCKEv8vbtWwDo3Llz/d0uXbpIO+7mzZsBwMPDo6ysTNqxvs2KFSsKCwt//vnn2bNnSz3Y0KEQFwcvXsC//4K7O2zZAn5+kJkJly+DzAfZIoSk4WO97sm4etljMBhkVmna+97n5+ePHTu2rKzM0dFxy5Yt9CaDvj90f1NAX4fH471586awsLB+4dixYwFg9+7dZDcgIAAATE1N6e1EVFhYyGAwVFVV+Xw+jWk09ubNGyMjIwCwsLAoKir6xrsUFFAnTlBTplCtWv1fa5K+PuXiQr18Sb17R7FYFAAVHf3hy+/epQAoOTkqN/ebHwQhJCmTJ08GgPPnz5PdWbNmAcCxY8dkEHr48OEA4O7uLoNYX+vOnTsMBkNJSSmF9k7vN29Srq4054AQapra2lpVVVUGg5Hz39jDn376CQCCgoJozCo2NhYAdHR06urqxIW5ubmPHz+uqqqSTQ41NTUDBgwAgMGDB3O5XNkERS0JVum/J48fP+7QocPEiRN79uz522+/kcJHjx6RUUnF//XiJjPkeXl50Zfp/5AmL7rGR33QixcvdHR0AGDo0KHl5eUSuCOfT0VFUW5uVJcu/6vYP3hA+fpSAFT37p+60MSEAqBOnpRADgihpunRowcAvHjxguySt0zZ9MMkXZlUVVULCgpkEO7LVVdXm5iYAMDOnTtpTqWyktLSogCo8HCaM0EINUFwcDDUm7Wk8bh6unTv3h0AQkNDxSV79+4lbZ8GBga2trbu7u6BgYH5+fnSiC4UCidOnAgAJiYmze0vAvS9wI7335Ndu3YtXrz46tWrz549O3PmTEpKCgC4urqKRKIePXq0Juv6Ahw8eDAqKmrOnDl05goAzW8pu+fPn1tbWxcVFY0dOzYsLExNTU0CN2WzYfBg8PCApCR4+RJ27QJLS4iNBQDo1+9TF/bvDwD/OxMhRB+RSJSamspgMDp16kRKZNbxHgCsrKzGjh1bVVW1e/duGYT7chs3bkxPT+/VqxeZGIVOqqrwxx8AAOvWfXiVO4TQ96Bxr3uRSGRtbU1WX6aRk5MT/P9971VVVXv16iUnJ5efnx8cHLxp0yZ7e3sDA4O2bdva2dlt3LjxypUr6enpEom+cuXKa9euaWlpkQWYJHJP9KPBKv33RFNTMzIysqamJj8/n8lkamlpAQCbzWYwGEwms7a2Vnzm4MGD1cngbVo1Hk7P5XKNjIwcHR1PnTol47Gjd+/eHTlyZElJyYQJE65evaqkpCT5GD16wOrVwGJBcTEAwKd/l/X1AQA4HMmngRD6GpmZmbW1tW3atCErrhcWFpaVlWloaJAePTKwdetWBoOxb9++3Nxc2UT8rLi4OB8fHxaLdfz4cTk5ObrTAVi+HNq0gefPITiY7lQQQt+CoqiQkBBoHsvXNTB9+nQGg3H16lUul0tK5s2bFx8fX11d/erVKz8/P1dXV2tray0trdzc3ODg4K1bt06dOrVjx47q6uqDBw9euHCht7f3gwcP6r+KfyFfX18fHx95efnLly+L53NB6KvR3U0AfYXCwkJtbW0tLS1DQ0N/f39SuG/fPvJ1c9iwYRwOh94MG0hNTQUAHR0d8aj+mzdviv/sycnJjRw5cu/evTIYI6mvdgAAIABJREFUpRkYGKioqAgAM2fOlMXY/okTKQBq7dpPnbNqFQVAOTlJPRmE0CeFhYUBwIgRI8ju/fv3AWDAgAGyzIGsjbxkyRJZBv0YPp/ft29fAFi9ejXdudTj7U0BUD17UkIh3akghL5aTEwMALRp04a8E/L5fG1t7frj6unVv39/ALh06dKnT8vNzQ0MDHR3d7e1tTUwMGhQq5KTkzMzM5s1a5aXl1d4eHjx5xY2CgkJYbFYDAbj1KlTknsU9CPCKv33xMHBYeXKlfn5+Tt27GjdunV8fDwpT0hIaNeuHQB07NgxKSmJ3iQbaNOmDQDUzyotLc3X19fW1lZeXl78I2hiYuLq6hoeHl5/bhJJOXv2LGllWrJkiVA274Jz5lAA1OLFnzrHxYUCoP6bEwEhRBdvb28AEE9QQtZY+uWXX2SZQ1JSEpvNlpOTS09Pl2XcD9qxYwcAGBsbV1ZW0p1LPTwe1aEDBUBduEB3KgihTykuLp4wYUKXLl26devm4eFBCskspLNmzRKfVl1dfefOHZpybMjFxYUM7LeysnJ1dfXz83v16tVnXxrz8vJu3Lixfft2BwcHU1PTBvP2MxgMU1NTBweH7du3P3nypMG1MTExqqqqALBlyxapPRb6UWCV/rshFArZbLZ4BfUpU6aIfyUpisrNzbWwsAAALS2te/fu0ZTjBzg4OADA4cOHGx+qqqoKDAxcsGCBnp6e+OdPS0vLwcHBz8+vpKREIgkcPHiQyWQCgJubm0Ru+EW2bqUAqP8a/T5s0CAKgNqzR1Y5IYQ+jKz3uee//xn//PNPANi6dauM0yDT7M+dO1fGcRtIT09XUVEBgJs3b9KbyQccOUIBUKamVDNbSAUhVN8ff/xBvoqWlpbq6OiQmUfJyKYVK1bQnV1DfD5/9erVDAaj8SAjVVVVKyurJUuWHD16NDo6ura29tO3qqioiI6OJh31rays6o/xXL58ef0zc3Jy2rZtCwBz5syR5sOhHwVW6b8nAwYM2LhxI5fLff36taGhYfj/P/dvVVWVvb09ACgoKJw+fZquJOsrLS0dNmyYsrKysrIyqaiXlpY2Pk0gEERHR7u7u5OvEgSLxbKysvLw8EhMTPzmBDw8PMhX0l27djXhOb7evXsUAKWoSH3oeSmKot6/p+TkKACqOS0HgNCPadSoUQAQEhJCdidMmAAA4sFNMpORkSEvL89isd68eSPj0GJkqioAmD17Nl05fIpA8L+1RU6coDsVhNBHLV26dOrUqSKRiM/n6+joREdHUxT1888/MxgMNTU1X19fuhP8Pzk5OWTiJzab7e7unp2dTfrVOzg4mJmZNWh1Z7PZZmZmDg4OZAL8BktKN1ZXVxcXF3fy5Mnly5fXn06/oqKid+/eADB06NDPfiZA6Etglf57kpWV5eLiYmFhMXTo0DNnzjQ+QSAQuLq6kkqsu7s7jevSi0QiPz8/0vxe/yOlePy8uLtBA+np6Z/ols/j8b48B1KfZzKZH+wjIF0iEWVqSgFQ69Z9+ARXVwqA6t1btmkhhD6ADFwST+rRrVs3ABCPbJKlhQsXAoCzs7PsQxPHjh0DAG1t7c++qtLm7FkKgGrfnsL3YISaq5ycHGNj4z59+lhbW2/bto0UpqamOjs7k/e6yZMnFxUV0ZskRVERERH6+voAYGho+ODBg8YnlJaWRkVFeXl5zZo1y8zMjMViNWjGFy9x5+/v/+rVqy958RYIBLa2tgDQrVs3SfVIRQir9C2Qr68vm80mnXm+qg4sKcnJyaSdBwCGDBny8uXLbxg/L+6WT35tCRUVFVtbW19f3/fv338iAZFIRFZdYrFYJ+la+D0oiAKgmEzK2/v/m8xJIKA8PCgGg2IyqYgIenJDCP2npqaGyWTKy8uTiTMFAoGCggKDwaBlneTc3FwlJSUGgxEbGyv76O/fvycLqZw9e1b20b+UUEj17k0BUP/+S3cqCKEP279//+DBg8PCwmbPnt2mTZv6Eyr5+/tramoCgK6u7vXr1+nKUCQSeXh4kCr68OHDv3DBeR6PJ54A38rKqvHaexoaGuKh+NHR0R98CV+8eDH5ciqDyaHRjwOr9C1TYGAgGQzp7PybLL8AVldXu7u7k3q7lpaWr69vgw+Wnx4//7Fu+VFRUX/++aeZmVn9vk/Dhg3btWtX44lSBQLB3LlzyQCEgIAAKT7tZ+3aRTEYFADVuTO1ZAn111/U4sVUx44UAMViUQcO0JkbQoiiKIqKj48nrSVkNyUlBQDatWtHVz4rV64EAHt7+/qF2dnZ//777/Hjx/39/cPCwiIjI6Ojo9++fZudnV1SUiIQCCQS2tHREQBsbGwkcjcpunpVqK19xdaWls8uCKHPMjIyiviv0WLOnDkN1s7IzMwcPnw4eZ2bNWuW7KfhLCoqGjNmDOnT6ubm9s0/oQKBIC0tTTwBfuMl5eXl5cUT4EdFRVVXV3t6epLuq48ePWpwN19fX4f/hIWFNfkp0Y+FQVEUoJYoISFh9ux5RUWXlJTah4SADJa6DAoKWrZsWWZmJpPJnDFjxp49e7S1tT92slAojIuLCwoKCg4OJuuaAACLxbK0tLSzs7Ozs6tfgRd79+7drVu3goKCSCd8AAgPDxf3CACAurq6GTNmXL58WUVF5erVq2SILJ3u3YO//4aoKBAK/1fCZsPPP8PmzTBwIK2ZIYQAAC5duuTo6DhhwoRr164BwI0bN8aPHz9q1Khbt27Rkk9xcbGJiUllZeXjx48tLS1J4e3btz/9ayYvL6+ioqKurq6oqEg2lJSUlJWVG29oaGgoKioqKysPHDiQfPklyIMrKyu/fPnSxMREug/ZRBQ1dtiwsKionTt3/vHHH3RngxBqyMrKytHRcfny5Xw+f/To0fb29uRjpRhFUUeOHFm5cmVNTU2HDh38/PyGDBkim9yePXvm6OiYmZmpra195swZUreXlMzMzLi4uLi4uNjY2Li4uMzMzPpHSRdakUjk7+9PFi6tb9KkScOGDSP/HYyNjVu3bi3BxFDLR/c3BSRFWVmiXr0oAEpHh2r0NVCS0tLSxo0bR/5E9enT5/+xd+dxMed/HMBfM92HDnKUI+SqrCL3sa5c+8s6c6xz2a2wWJLCrmLtKrmyciTWsbss61i5S7nlSORILRUpokMpnTPf3x9j2zYh1ExTr+ejP2a+n898v6/24WG953NdvHjxvT7+Aevn09PTd+/ePXny5KKtWVlZsr+aDQwMXv/6U5HS0oTQUCEwULh0SUhPV3QaIvrXkiVLAMydO1f2duXKlVD0+fALFiwAYGtrW3jlzp07Tk5OEydOtLe379u3b7du3WxsbJo0aVK3bl1DQ8NiGziVUkREROH9MzIy6tevD2DVqlWK+I3f2/HjxwHUqFEjnX+jElU8N2/e7NixY8uWLRs1ajRx4sQ37QB3+/btNm3aAFBRUXF1dS3vtaJSqbBhwxVZXd2tW7eEhIQSuxUUFNy+fTsqKurjd6R6/vy5bCm+g4NDly5d1NTUqlev3rVr1xI7t2/f/gq3TKYPxVH6Si4zE6NH49AhaGhgyxZ88UUZ3z8nB35+N9zcOmVnZxsYGPz4449OTk6yQ+M+QFZWVnBw8KFDhw4ePPjkyRPZRUNDQ1tbWzs7u4EDB8rWX5UoPT3dzs7u3LlztWvXPn78uGwr0QonMRFeXsjPx7p1io5CRAAwYcKE7du3+/v7T548GcCUKVM2bNjg4+Mj22pUIZ4/f25mZpaamhocHFw4PfXtcnNzX758+fz58+zs7Ozs7MIXaWlpOTk5shfZ2dk5OTmyF9nZ2Vu3bpUdoQRg+vTpa9eubdeu3cWLF1/f/6li6tmz56lTpxYtWrRw4UJFZyGiEhQUFMjq57fIy8tzd3f39vaWSCSjRn23cOEP5ublEiYjA5Mn46+/YG4+tUcPteXLl79+ZB2Ay5cvjxw50sTERFNT09bWdt68eWWY4erVq+3atatevfrjx4+LDmLJNGjQoHfv3klJSdbW1vPnz5cdWU9UWor+ToHKXUGB8M03AiCIRIK7e1neOThYMDcXNDSEhg372tvbv2XLOqlU+uTJk+zs7FLe+U3H2tnY2Li7u8tOQykqJSWlffv2ABo0aBAVFfXhv1J5e/ZMEIkEHR2hpB0BiUj+OnToAODMmTOyt7169QKg8HWMsrkDnTt3lsOzQkNDVVRUVFVVZWdHK4tz584B0NfXT05OLnr96tWrs2fPdnJyOnLkiKKyEdF7OXfuXMeO3Q0N0zQ1BU/P/+wpXCauXBEaNRIAwdBQOHQo/03d8vPzTU1Nt2/fXsaPL6Jly5YocmZqUSdPnrx8+XJ0dPSQIUPGjBlTfhmoUmJJX1WsXi2IxQIgfPllGZSTjx4J9vYCIACCpaVw9mzmWzofPny4UaNGFhYWZmZmgYGB7/uswmn5GhoaheV948aNHRwcDh48mJOTEx4e3qJFCwCNGjW6f//+R/xaciE7VPnyZUXnICJBEATZHu9JSUmyt3Xr1gUQExOj2FQvXryoWbMmgEGDBrm5uXl6evr4+Pj5+e3evfvQoUNBQUFXr169c+dOTExMamrqxxxrnJOTY2lpCWD+m07crMD69+8PwNXVtfBKXFycgYHBunXr/vzzT9mMLQXGI6LSy8p6dbwvIPTuLTx8WGZ33rZN0NISAKFNG+Ht/0i8dOlS3bp1JRJJbGxsOS3q+fHHHwGMHTv2LX0uX75cs2bN8ng6VWKceF+FHD2KkSPx4gV698aff8LA4ENuUlAAX18sXIiMDGhrw8UF8+fjtdlD/3r06JGlpeXhw4cL1w598LT89PR02d54R48eTU5Oll2U3U0qlZqbm4eEhBTdSL+C+vpr+PtjxQrMnq3oKERV3bNnz2rVqqWvr//8+XMAWVlZ1apVU1dXz8rKUuz88/T0dGtr65ycnMIlSG8nFov19fW1tbW1tLSKvtDR0dHU1Cz2QktLS09PT0dH5+nTp0uXLr106VLjxo1v3bqlpaVV3r9X2QoLC2vXrp2mpubff/8t+y7mu+++e/Lkib+/P4C1a9ceP348ICBA0TGJqLSOHcOkSXj8GPr6WLYMDg4fdbfMTDg4YOdOAHBwwJo1KDI2VII9e/bMnj3b2NhYU1MzMjLSx8fni7JesPrgwYNGjRppa2snJSUV3aBUIpGcOnWqd+/eADZt2rR58+bQ0NCyfTRVbu9Y4kKVyYABOHsWdnY4eRLt2+PwYZw4gdOnoaKCdevw+ip1Fxc8eICFC9Gy5asr585h6lTcvAkAdnZYuxampu946MGDBz/99NOuXbtmZmbq6up+2DZOMvr6+rKzPaRSaXh4eEBAwK5du6KiogDo6uoGBAQoQT0PoFs3+Pvj7FmW9EQKJ/sLpHnz5rK30dHRgiA0bdpUsfX88+fP+/btGxcX17Bhw0WLFqmpqaWnp2dlZeXk5MheZGdnZ2RkFH2Rl5eXlpaWlpb2YU+cPn260tXzAGxsbIYMGbJv3z4vL681a9YAuHXrVuFere3atVuxYoVCAxLR++nfHzduwMEBBw7A0RGBgdiwAR+29fvduxg+HLdvo1o1+Plh1Kh3f0QikVSrVi04OFhXVzc4OHjkyJFlXtKbmpp26NAhNDT00KFDI0eOLLz+5MmTWbNmqamp1axZMzIycu/evWX7XKr0WNJXLVZWOH8edna4eRPz5kFXF3v2AICBAdavL945MBA3bsDRES1bIjUVixZh7VpIpWjcGD//jH/+1fQOcXFxKSkpVlZW2dnZ6urq+/fvb9q06Uf+FmKx2MbGxsbGxsPDIyQk5NatWz179hSJRLJ/AX/kzctd164AcO4cBAEf8QUHEX286OhoFCnpZRV+Mzmc+flmGRkZ/fv3v3LlSpMmTU6dOiUbfH4niUSSkZGRmZmZk5NT7EV2dvaLFy9ef/HixYs7d+4IgjBhwoTPP/88JiYGgIGBgWwlgrJYvHjxX3/9tXHjxlmzZjVq1CgjI6Nw4EtXVzcjI0Ox8YjofdWsif37sWcPHBzw5584exb+/rCzA4DVq3H4MMRibN0KY+PiH/zmG0RFYckSdOiAHTswZQqysmBujj17YGlZqkebmpq+fPlS9neIqanp8+fPJRJJmX/DO3r06NDQ0J07dxYt6evWrRsREfHgwYNnz55ZWVmVuHUf0VuwpK9yGjTAuXNwccFPP8HZ+dVFPz+MH1/yQemCgO3b4eyM5GSoqcHFBR4e0NQs7eOys7Nr1qwZHBysqanp5uY2b968P//8s2x+EwBA/fr1Bw0aJNvgqmbNmr///nsZ3rxcNG6MevXw6BGiotCihaLTEFVpspK+sIYvVuHLX1ZW1sCBAy9dutSgQYPAwMBS1vMAVFRUDA0N33ImyFvMmjXLzs5O9qzJkyePKs1gVoVhaWk5evToX3/99ccff/T3969Tp87Tp09lTUlJScav/6ufiJSBvT3at8eECTh9Gp9/jq+/xsqVuHsXQUEAMGsWdu0q/pErV3D5MqZOxeTJ2LIFAL78Er6+KP0MpA4dOmhpaS1btux///vfDz/8MHDgwPKYsTVq1ChnZ+ejR4+mpqYW+wrV1NTU9J3TX4lK8oGrmkm5XL58eUQRbm5TN278dyJTr16QSuHoiPz8Ej4rlWLFCiQno2dPRETA0/M96nkA9erV09fX19TUBGBlZfXo0aOP/3WKSkhIsLGxCQwMDAwMVIJ6XqZLFwA4e1bROYiqumLD8sUqfDl7+fKlnZ3dmTNnGjRocOrUKbnNOUpMTFy0aJHsb1HlqudlPDw81NXVMzMzpVJp7969Dxw4INsk6MCBA7KFqUSkjExNcfIkPD2hpgY/Pzg5/dv0xx84frzkT6mo4O+/oamJjRuxZct71PMAxGJxYGBgVFSUk5OTsbHx1q1bPyb/m9SqVatHjx55eXn79+8vj/tT1cSSvkpo1qyZ6z9atGgRHx9ftHX4cHTogJs38fPPJXxWttJ+zx4EB3/IoLK9vX1AQMDFixcTEhJ8fX379u37ob9EyR4/fmxsbBwfHy+bNaocunUDWNITKV7FmXgvq+dPnTpVv379kJCQRo0aye3Rjx8/NjIyioqKyszMlNtDy5CZmdn9+/d37dolFou/+OKLrKys3r17jxw5cv/+/c6FU9GISAmpqMDVFZcvo1MnuLu/uihb+OnkhKyskj+ycycuX/7ArfXq1au3ZcuWc+fOrV69Wk9P70ODv8Po0aMB7JRt3EdUFjjxvkowMDAoPOB96tSpXl5exTr4+qJ9e3z/PYYOxesjQ7JB5Q/TuHHjHTt2zJ0799mzZwMHDvzuu+8+/F4lUVdXj4iImD179u3bt5s2bbp3715V1Qr/p7pwOT0RKY5EIrl//75IJGrSpInsiqenZ0REhGUpl12Wnezs7M8//zwkJKRevXohISGNGzd+vU90dHRSUpKlpWWZr3XX19d3d3c3NDS8ePHixo0bhwwZUrb3l4N69erJXmhpaV28ePHatWvZ2dlbtmwpuqE0ESmLa9eubdy4sfBtjRo1Llz4qfBtnz7Iz0dgIJYuxZIlJXy8bl2UetGSYgwfPnzatGkhISEJCQmlX2BF9BYVvvihMnXq1Knc3NwePXoUu25jg7FjsX07Zs1Cmc8DsrOzs5NtbFIOhg4dOnToUAD5+fktW7Y8dOjQ4MGDy+lZZeaTT2BggNhYxMejfn1FpyGqKnJzc+/fv1+zZk3Zke9XrlzJzc2tWbNmYeHXu3dv+U/VzsnJGTRo0MmTJ01MTIKDg83MzIp1yMjIGDJkyNOnTy0tLZ8/f37kyJEPPgq0RIXHvB0+fNjR0VEZS/qi1NTUZLurEJGSqlevnr29vex1YGDglStXinXw8YGVFby9MWYMzM3lnu+j6enpLfv662YXLtTeuxczZig6DlUGnHhftaxYscLFxaXEpmXLYGCAAwdw8KCcQ5UNNTW1hg0bFp5XX6GJxejSRVKnTtJr/5cionJy5syZFi1afP/997169ZoxYwaAwMBAAFlZWVklTt+Ui7w8jByp8fLl8Fq1agUGBpZ4IIhsCP3GjRu7du06duxY2dbzRTVv3jw5OVm2EJ2ISFFq1apl+48LFy7Mfu3QX3NzTJ2KvDxMnQol/Rtrerdu/cLCVJVlEyiq8FjSVyHR0dE3btwo/OKzmNq18cMPADBjBl6+lGuwj+Hq6url5XX79u3t27dfvXq1f//+ik5UKht79lR98mTxyZOKDkJUVSxbtszZ2Xnv3r2XLl3y8/N7+PDhwIEDxWLxy5cvGzZs6OHhIf8Dz/LyMHQoDh4UxcQ4nDp11cLCosRuAQEBM2bMePjwYWRkZJlnSE5OHjBgQGBgYGRk5Pz584cNGybi4ZpEVDFcuXIlOTm5X79+rzctXgwTE5w6he3b5Z+rLAwcCD09XLqEv/9WdBSqDFjSVyErVqyYPn26urr6mzpMmYLWrfHgAZYtk2eujzJ16tSkpKSZM2eeOHFCthJV0YlKxbJDBwDnuJyeSF4MDQ1PnjyZnZ398OFDbW1tAwMDa2vrDRs2dO7cOTk5edGiRU2aNPH29pbbiH1+PkaMwOHDqFkTx4/D3LzkNThSqfTRo0cLFixwcHCwt7cfMGBAQUFBGcYwMjL68ssv/f39HRwczM3NN23aVIY3JyL6GN7e3nPmzClxapKeHmQbQ7m4IC1N3sHKgKYmBg0CgN27FR2FKgWBqoanT5/WqlUrLS2t6MUJEwRAWLfu3yuhoYJYLGhpCQ8eCFZWAiAEBck7alWQk5OjqakpFotTUlIUnYWoSkhMTDQ0NKxRo4axsfHBgweLNp09e7Znz56y/ycaGRm5u7unp6eXa5iCAmHkSAEQDAyEq1ff1lMqlaqrq4eGhgqCkJ+fb2Zmdvjw4XLNRkRUEcTGxhobG2dnZxdecXQUAGHVqn/79OwpAMKMGYIgCO3bC4Bw6JDcg36wI0cEQDA3V3QOqgw4Sl9V6OjohIeHGxgYvL1bhw6YPBnZ2Zg/Xz65qigNDY127dpJpdKLFy8qOgtRleDg4ODk5HTr1q1p06Z9+eWXRSexd+3aNTg4WFbYy0bszczMym8qvkSCcePwxx8wMEBgIP45jaRkIpHI1NT05cuXAFRVVatXr56enl4eqYiIKpTVq1c7OTlpamq+pc/atVBTw7p1uHNHbrnKTp8+qFULkZGIiFB0FFJ6LOmrCm1tbRMTk9L0XLoURkbYuRNxceWcqWrr2rUrOPeeSC7y8vIOHz48bdq0OnXqLFiwoGPHjidOnCjWRz6FvUSCCROwcyf09XHiBNq2ffdHJk+e/N1334WGhq5cuTImJsbW1rYM8xARVUDPnz//7bffHB0d397NwgLffouCAsydK59cZUpVFcOGAQAPqKePxpKeiqtRA56ekErBoaBy1a1bNwBnz55VdBCiyk9dXd3KymrLli05OTk3b968evWqzRsGx4sV9kuWLOnZM9HDA2VS10skmDgRv/32qp5v165Un5o7d+4XX3zh6el548aNkJAQ2Ql8RESVWEFBwdGjR2vXrv3OngsXon59HD6snAP1o0cDwO+/K+vG/VRhsKSnEkyahM6dFR2isuvSpYuKisrVq1dzcnIUnYWo8jtw4EBcXFy3bt1mzpy5Zs0a2TSZN5EV9sePHx8zZvG1ay0WLUKTJli2DB+zd55UikmT8Ouv0NFBQADaty/tB0Ui0bRp0w4cOLBt27ZPPvnkwxMQESkJIyOjtqWZxQTo6mLVKgDIzCzfSOWia1c0bIiHD8FlmPRxWNJTCUQibNgAVVVF56jU9PT0WrZsmZube4Wn0xOVP1NT082bN1+5ciU4OHjEiBGl+Ujfvn23bZt/9ix69sSzZ3B1RcOGKHHEfvFitG2LTp3w5EkJ9xk7Fm3b4tAhBAVBVxfHj6Nbt4/+fYiICAAwbBg++0zRIT6MSATZ2dKce08fhyV9lbZ1KwQBU6aU0PTJJ8jPhyCgd2+5x6q88vPzIyIinvzzr37OvSeq+Lp2RXAwZIV9cjIWLYKZWfHC/sEDhIUhNBTOziXcITISYWHQ0MCpUzh6FF26yC07EVHloa0NQ0OUuF+ejw9q1YKhIdTU5B7rI8nm3v/xB8r0gFKqaljSV2nbt2P5cty/X3LrwYNYvhxXr8o3U+V18+ZNU1PTuXPndunS5auvvgJ3yCNSHqUp7AH8/jsCA994k6ZN8db5/kRE9EYrVyI1FU5OJTQ1aYKkJKSmom9fucf6SK1bY9w4LFoEiUTRUUiJiQTux1CFWVkhIgIHDmDQoBJahw3Dvn3w9sacOXJPVhkNHjy4Xbt2CxYsSE9PNzc3P3DgQP369U1MTPT09FJTU1VUVBQdkIhKJTAQixbh/HkAqFkTK1ciJARbtsDODocOoUkT3Lz5n3EkGxtcu4Zjx9Cvn6IiExEpPT8/xMdj5Ei0bFlC659/4sYN9OmDTz+Ve7KPJ5Xi+nXExCAjA0ZGaNkSjRsrOhMpE47SE8mDIAjHjh0bM2YMAH19fTs7u6NHjxobG8+bN8/X11cqlSo6IBGVVp8+OHcOhWvsq1V7db1fP9ja4t49LF2q0HxERJXR5s1YsgS3b5fcun8/liyB8q1lLCjAihUwMYGNDeztMXkyBg2CmRmsrHD4sKLDkdLgBmhE8pCampqbm2tsbCx7a2JikpiYCOCnn35SaC4ieg8SiaTwmHpLS+zdi1u3DLt2xcGDrzr4+MDaGp6eGDUK5uYKy0lEREogPx/DhiEgAGIxPv8c3bpBTw+JiQgIwLVrGDgQnp6YO1fRKUkJsKQnkgddXV2RSJSdna2hoQEgKytLT09P0aGI6P3cvXv3888/l70uKChISEjIy8sTif6d72ZhAUdHrF2Lb75BUBBEIgUFJSKiim/hQgQEwMgIAQHo2PHf6+7u8PHB7Nlwc0ObNrC1VVxEUg6ceE8kDxoaGqamptevX5dKpUhHAAAgAElEQVS9vXHjRvPmzRUbiYjel6Wl5f1/uLu7DxkyRCwu/r/RH3+EsTGCg/HbbwrJSEREyuDpU6xaBQBbtvynngcgEuHbb+HgAEHA/PkKSUfKhSU9kZxMmTLF1dU1LCxsw4YN169fHzlypKITEdGH8/HxmT179uvX9fTg6QkAzs54/lzeqYiISDns2YPcXLRogYEDS+7g6gqRCFeu4O5d+SYj5cOJ9wQXF/z4YwnX792Te5RKzdnZWSQSzZs3r3r16oGBgdUK99QiImVz9OhRXV3dTp06ldg6bhx++QWnTmHRoldjMERERP9x4QIA9O79xg6NGqFJE/z9Ny5eRIsWcstFyoglPSEmpuQFnzwgs2ypqKi4uLi4uLgoOggRfawVK1Y4Ozu/qVUkws8/o00b+PrCwUGeuYiIKrlZs7BwYQnXHz+We5SP9OgRAJiZva1P06b4+2/Ex8snESkvTrwn7N2L/PwSfoYMUXQyIqKK5+bNm7GxsYMGDXpLn5YtMWMG8vPh6iq3XEREld+LF0hOLuEnN1fRyd7XixcAoKPztj66ugDwz0krRG/Ckp6IiOg9eHt7z5o1S0VF5e3dPDxQrx4CAhAVJZ9cRESVn78/UlJK+Bk+XNHJ3pe2NgDk5Lytz8uXAKCri7AwdO+Ob7/F1q24fh35+fJISMqDJT0REVFpvXjxIj4+fuLEie/sqauLlSsBICurvEMREVHFkpubu3nzZmtrax8fn5J71K4N4B2T6mWtderg0iWcOQMfH3z5JVq3hrY2LC0xfjx8fBAUhJSUso5PSoZr6YmIiEqrWrVqISEhpexsb4/PPsORI+WaiIiIKpxt27YlJCSYm5snJyeX3KNdO+zb92qTvBKlpeHOHQBo3x6mpmjcGOHhuH4d4eG4fx937uDOHezYAQAiEczMpG3berVqZWVl1bp1a2Nj43L4najiYklPRET0UbS1YWgITc0SmlavxpUrKCiAmprcYxERkYI4ODgAWLBgwRt7DBmC+fNx8SJu3ULLliV02LwZ+fkwM4O1NcRi9O+P/v1fNWVmIioKt28jLAxhYbh2Dffu5aqozN+1S9ZuaGhoYWFh8w9zc3OxmFOzKzOW9ERERB/l55/x888lNzVtiqdP5ZuGiIgqvubNMWQI9u3DxIkICoKBwX9ar13D4sUA4OaG16txXV3Y2MDGBuPHA0B+PiIjU+/e/fbixevXr1+/fj0tLe38+fPnz5//p7uulZWVtbV169atR40apfP2PflICbGkJyIi+ihr1yIhAWPHwtKyhNadOxERgQED8Omnck9GREQV1vr1uHYNYWFo0waurujeHYaGSEjAgQNYuRJZWRg2DJMnv/s+ampo1apuq1arRoyQXUhMTAwLCwsLC7tz587t27cjIyNlFb5IJLK3ty/fX4oUgSV9lbZjB7Ky0KJFya1Ll2LOHDRsKNdIRERKZ8sWhIejXbuSS/oDB7B7NwwMWNITEVERtWrhwgV8+SWOH4eT03+a1NXh4oKffoJI9AE3NjExMTExGThwoOxtcnJyeHh4eHh4YmKinp7exwenioYlfZXWqtXbWps1k1cOIiIiIqI327ULL1+ifv2SW5ctw/z5r3aRVybGxjh2DFev4uhR3L+Ply9hYABrawwc+MZf9f0ZGRklJiaePXtWRUVl9+7dI/4ZzKdKgyU9ERERERFVaI0ava21bl3UrSuvKKVw+fLl0aNHZ2dnA/j9999v3Lihq6v7xt5t26Jt2/ILc/z48QULFuzevTs3N3fkyJH169fv1KlT+T2O5I8lPRERERERUZlp3779/fv3FZ3iFT8/P2dn586dOwNwcnLy9/dnSV/J8DwDIiIiIiKiyikqKqrVP6ttW7VqFR0drdg8VOZY0hMREREREVVOOTk56urqsteampovX75UbB4qc5x4T0REVAamToWLSwnXk5LkHoWIiOgfJiYmiYmJstePHj2qV6+eYvNQmeMoPRERURnIyUFmZgk/BQWKTkZERFWYnZ3djh07BEGQSCQ7d+4sPNyOKg2W9ERERGVgyxYkJZXwM2iQopMREVEV5uTklJKSYmNjY21tra6uPnbsWEUnojLGifdERERERESVk56e3sWLF+Pj41VUVExMTBQdh8oeS3oiIiIiIqLKrH79+oqOQOWFE++JiIiIiIiIlBJLeiIiIiIiIiKlxJKeiIiIiIiISCmxpCciIiIiIiJSSiJBEBSdgYiISIk9fIicHJiYQFe3hNbHj/HiBWrWhKGh3JMRERFRZceSnoiIiIiIiEgpceI9ERERERERkVJiSU9ERERERESklFjSExERERERESkllvRERERERERESoklPREREREREZFSYklPREREREREpJRY0hMREREREREpJZb0REREREREREqJJT0RERERERGRUlJVdACq6vbt2xcQEKCrqzt+/Ph27dopOg4REREREZHSYElf5Zw8eTIgIEBDQ2Ps2LGffPKJYsOsW7fu0KFDU6dOTU5O7tev39WrVxs3bqzYSERERERERMqCE+8rp9s5t/c/3x+ZE1ns+tatW+fMmdOlS5cGDRp069YtNjZWIfEKjR8/fs+ePXZ2dhMnTmzYsOGtW7cUm4eIiIiIiEiJiARBUHQG+kCZ0sywl2HPCp6pidRM1U2ttKxEEMma3BLcvJK8vjf+frHx4qIfSU1Nzc/Pr127NgA7O7uBAwc6OjoqIHoRgiCEhoaePn36jz/+OHPmTLVq1RSbh4iIiIiISFlw4r1Supd7b37i/APPD+QL+YUXG6o3XFp36SjDUW/5YPXq1QtfP3/+3MjIqBxTlo4gCFu3br1z5465ubmqKv9AEhERERERlRYrKOVzMevi/+79L02SVk2l2gj9Ec01mksguZB1ISgjaHTs6LjcOLc6bu+8ybFjxx49ejRgwAA5BH47sVi8ceNGALa2ttu2bXNyclJ0IiIiIiIiIuXAkl7JpBakDo8ZniZJ66zTeW/jvXXU6hQ2nXxx0uGhQ2fdzu+8yZUrVyZPnrxv3z5tbe3yDPtu0dHRNWrUqFGjBgBtbe3MzEzF5iEiIiIiIlIiLOmVzPrk9Yn5iUaqRn+Z/WWk+p9p872r9Y62iFYRqbz9Djt27Fi8ePG+ffs6dOhQnklL5cyZM56env369Xv27NmtW7f8/PwUnYiIiIiIiEhpcMd7JbMrbReAr42+LlbPy7yznl+yZMn48ePNzc23bNni6Oj422+/lUvKUvvqq69OnjzZo0eP0aNH37p1q06dOu/+DBEREREREQHgKL1yeSF5cSf7DoA+1fp82B169erVvHnzwrdmZmZlk+xDSaVSd3f33377TRCETz/99K+//uKO90RERERERKXEkl6ZPC54LIUUQFPNph92h86d373SXm5CQ0OnTp0aHh4uexsSEtKiRYvFixdPnDhRReUd0w2IiIiIiIiIE++VyQvJC9kLXbFuKT/yXPLcP9m/QCgot1AfIi0tbebMmV27dg0PDzcxMZk1a9b69es7duyYmJj41VdftWzZcs+ePYrOSEREREREVNGxpFcm2uJXG9TnCXml/MhPT376+uHXLSNb7kmrEEWyIAjbt29v3rz5mjVrxGLxjBkzoqKiVq5c6eTkdPHixYMHDzZp0uTu3bsjRozo0qXLhQsXFJ2XiIiIiIio4mJJr0wKt8SLz4sv5Uc663RuotEkKidqROyIHtE9rry8Um7p3u3GjRtdu3adMGHCs2fPevTocf36dR8fH13df2ccDBw48M6dOxs3bqxVq9aFCxe6du06YsSImJgYBWYmIiIiIiKqsFjSK5OaqjXrq9cHcD7zfCk/Mthg8B2LOxsbbKytVvt05ukOdzuMiB1xP/d+ecYsQXp6+syZM21sbC5cuGBsbLxt27bg4GALC4vCDk+ePBEEAYCampqDg8P9+/fd3d01NTX37Nljbm7u6Oj47NkzOWcmIiIiIiKq4FjSKxk7fTsAm1M2SwRJKT+iJlJzMHK4Z3nPs66njlhnT9oe8zvmjg8dnxY8Lc+k/woICGjZsuWaNWtEItGMGTPu3r07fvx4kUhU2EEqlX722Wft2rULCQmRXdHV1fXw8IiOjnZwcJBIJH5+fs2bN/fy8srJyZFPZiIiIiIiooqPJb2SmV5zurpIPSI7YkHigtdb37LGXles61rb9a7lXQcjBymkfsl+ZrfMPB57ZEuzyy9tdHR0v379Pv/880ePHnXt2jUsLMzHx0dPT69Yt9jY2GfPnoWFhfXq1WvgwIG3b9+WXa9Xr97GjRsjIiLs7OzS0tLc3NyaNWvm5+cnlUrLLzMREREREZGyYEmvZMw1zb3qegHwSvLqd6/f/uf77+Xei8mNOZ5x/NtH3xrfNH77pPq6anU3Nth40/ymvaF9pjRz0eNFzW4380v2K/2Yfym9fPnSw8OjVatWJ06cqF69+urVq0+fPt2qVasSO5uZmUVHR3t6eurr6x86dKhVq1YjRox48OCBrNXCwiIgICAwMLB169bx8fGOjo7t27cPDg4u28BERERERERKRyRbwEzK5ZeUX1wSXFIKUopd1xJr7W60207fzi3BzSvJ63vj7xcbL37TTYJeBM1NmBv+MhxAa+3WKw1X9qjdo0ziBQQETJ8+/cGDByKRaOzYsStXrjQyMirNB1NSUry9vVevXp2bm6utrT19+vR58+bp6+vLWqVS6d69e11dXWNjYwHY2touX77cysqqTDITEREREREpHZb0yipdkn404+jFrIvP8p8BMFE3aaPVZoD+AEMVQwDBL4LPZ57/tNqn3XW7v+UmAoQ/0/50TXCNy41r8XWLujXqent7W1tbf3Cq+/fvz5gx48iRIwBat269bt26jh07vu9NHjx48P333//666+CINSoUcPFxWXWrFnq6uqy1ry8vPXr13t4eDx//lwsFg8bNszb29vU1PSdt338+LG1tXV+fn5qaur7RiIiIiIiIqqAWNITcqQ5W65uWdBvgaxInjBhwuLFi+vVq/deN8nLy1u1apWHh0dOTo6BgYGHh8c333yjoqLywamuXLni4uJy+vRpAM2aNVuyZMnw4cMLN9VLTU1dtmzZm8bzSzRs2DALCwtfX1+W9EREREREVDlwLT1BU6w5tf3U+/fvu7q6qqur//LLL82aNXNzc0tPTy/lHU6ePGllZeXm5pabmztu3LioqKiZM2d+TD0PoF27dqdOnQoMDGzZsmV0dPSIESM6dep09uxZWWv16tU9PT2jo6PHjRuXnZ3t5eVlZmbm5eWVm5tb4t327t0bHx/v4ODwMZGIiIiIiIgqFI7S0388fPjwxx9/9Pf3l0qlsknv3377rYaGxpv6JyQkzJs3b8eOHQCaN2/u6+vbu3fvso1UUFCwZcsWd3f3J0+eALCzs1u5cmXTpk0LO7x9PB9Aenq6lZXV/v37a9Wq9cknn3CUnoiIiIiIKgeO0tN/NGjQYOPGjZcuXerZs2dKSoqbm1vz5s23b9/++lc/+fn5Pj4+5ubmO3bs0NHRcXd3j4iIKPN6HoCqqqqDg8O9e/c8PT2rVat26NAhS0tLR0fHpKQkWYfC8fxPPvnk9fF8AM7OziNHjmzdunWZZyMiIiIiIlIgjtLTGwUFBTk7O0dERABo166dt7d39+6vNts7ffr0tGnTZAfI29nZ+fr6NmjQQA6Rnj17tmTJknXr1hUUFOjq6k6bNu27777T1dWVtRYUFPj7+3t4eCQlJYlEonr16m3ZsqVp06aNGjUyNzfX0tLKz8+/ffv2gAEDAgIC5JCWiIiIiIioXLGkp7eRSqW//vqrm5vb48ePAdja2i5YsGDLli2y7eibNGny888/9+/fX86p7t69u3Dhwj179gCoW7fuwoULJ0+eXLh0Pysra+3atQsXLszLy9PX1793715h0+PHj7t06RIXF/f2jfSIiIiIiIiUAkt6eresrKwVK1Z4e3tnZmbq6OhkZWVpa2svWLBgzpw5hWfLyd/FixddXFzOnz8PwNzcfNGiRfb29oWt169ft7GxkUqlJ0+e7NWrl+xiUlLSgAEDrl27ppjEREREREREZYolPZXWkydPxo4de+HCBRMTk6CgoIYNGyo6EQRB2LVr14IFC2JjYwHY2touX77cyspK1mpiYvL48eMLFy506tRJoTGJiIiIiIjKBbfHo3+lpKQUfZufn1/0HLs6depMnDgxOzu7U6dOFaGeByASiUaPHh0VFbVx48ZatWoFBQXJVv7L1KxZE4COjo7iAhIREREREZUjlvQEQRD++uuv7t27GxkZxcTEAMjIyPD29jYzMyu2g72amhqA/Px82dvMzMzIyEjZCLkCqampOTg43L1719vbe8yYMYXXZYsC8vLyFBeNiIiIiIioHLGkJ4hEopiYmNWrV9epU0d2JSkpycDAYMWKFcV6FiuSjxw5YmFh4erqKs+0b2JoaDhnzhyx+N8/0rIvIFjSExERERFRZaWq6ABUIcyaNavo26ZNmzZt2vTMmTPFuhUrkmUVfuGgfUVTweMRERERERF9JI7S03soViRX8GHwCh6PiIiIiIjoI7Gkp/fAUXoiIiIiIqKKgyU9vYdiRXIF33+ugscjIiIiIiL6SCzpCQCePXsWExMjkUji4+OTkpIAxMbGJiYm5ubmxsTEZGZmyroVK5KLbYBf0VTweERERERERB+JJT0BwLZt2xwdHa2srJYsWbJx40YAM2bM2Lx5c506dRwdHa9evSrrVqxIruDD4BU8HhERERER0UfijvcEAHPmzJkzZ07RKwEBAa93K3GUvsLWzBU8HhERERER0UfiKD29B26PR0REREREVHGwpKf3wEPsiIiIiIiIKg6W9PQeOEpPRERERERUcbCkp/fAQ+yIiIiIiIgqDpb09B54iB0REREREVHFwZKe3oNyHWLHkp6IiIiIiCo3lvT0HpTrELsK/o0DERERERHRR+K59PQeVFRUxGKxRCKRSqVisbhC7T/3119/JSUlyV4PHDjQ2Ni4gn/jQERERERE9JFY0tP7UVNTy83NzcvL09TUFIvFKioqEolEIpGoqKgoNpirq2vPnj319fUBCIIA7nhPRERERESVHUt6ej/q6uq5ubn5+fmampqyt9nZ2Xl5eVpaWooNlpCQ4O3traurW3iFo/RERERERFS5cS09vZ+Kuel9ZmamRCK5cOGCv7//33//LbuoRKP0t27d+uabb+zt7X19fSUSiaLjEBERERGRcmBJT++nYm56LwjCF198cfny5fj4+M6dOwcFBVWcbO+UmZnZu3dvMzOzmTNnbtmyZcWKFYpOREREREREyoET7+n9VMxN76tVq+bv7y97raOj4+vra2trW0FmELzT7t27W7ZsOWvWLADLly+fNGnS3LlzFR2KiIiIiIiUAEfp6f0Uq+ErzuT2nJwc2QuRSCQWi6E8o/SRkZE2Njay123atImLi8vOzlZsJCIiIiIiUgos6Qn5+fk7d+7s0KHD4sWLZVekUuncuXObNGnSrFmzhQsXyjaQlylWw1eQUfqjR49aW1uvXr16+fLlK1asmDFjRsXJ9k5ZWVmyvQYBaGtry64oNBERERERESkHTrwn7N69++bNmy1btiw8133nzp1Hjx69fPlyQUFB165d27RpM3jwYFlTxRylHzBggIGBwenTpwsKCk6dOtWiRYuKk+2djI2N4+PjZa8fP36soaFRo0YNxUYiIiIiIiKlwJKeMGbMGADe3t5xcXGyK7///vv06dOrV68OwMnJ6ffffy8s6YvVyRVncnunTp06depU9IqyjNL37dt36NChL168qFat2h9//NG3b1+RSKToUEREREREpARY0lMJ4uLimjZtKnvdpEmTP/74o7CpYh5iJxMVFZWQkNCrVy/ZW2UZpe/QoUOfPn3atm1rbm5+6dKlY8eOKToREREREREpB5b0VAJBEIquny/6umIeYpedne3l5eXp6WloaHj37l19ff2Kk600tm7deuvWrZSUFBsbG11dXUXHISIiIiIi5cDt8agEpqamsbGxstexsbENGjQobCpWJy9btiwwMNDCwkL+IQv99ddfFhYWixYtys/Pt7OzK/wCokLNIHi7kJCQUaNGaWhosJ4nIiIiIqLSY0lPJRgxYsSGDRtevnz54sULf3//UaNGFTYVW6DeuXNnW1tbQ0NDheR89OiRvb394MGD4+LirK2tz549u2nTJgMDA1lrcHAwgNjY2NDQUIXEK42EhIQRI0b06tXr9u3ba9asUXQcIiIiIiJSJizpCdHR0WZmZqtXr963b5+ZmdmjR4/GjRvXunXr+vXrN2rUqFevXkOHDi3sLDtf7ebNm4rLCwD5+fk+Pj7m5uZ//vmnvr7+6tWrr1692rlzZ1lrdHR03759nZ2dVVVV09PTO3fuPGLEiHv37ik2czEFBQU+Pj4WFhZ79uzR1tZ2d3f/5ZdfFB2KiIiIiIiUiajoMmmioqRSqUgkKrr7+q+//jpx4kR9ff3U1FRbW9tVq1a1bNlS/sFCQkKmTZsWGRkJwM7Obv369fXq1ZM1FS6qz83NNTQ0nDdvXkFBwdKlS1+8eKGmpvbll18uXry4du3a8s9czNmzZ6dNmyb7ZsTOzm7t2rWmpqaKDkVEREREREqGo/T0RmKxuGg97+vrO2HCBIlEYm5urqurGxQU1KZNm2nTpj19+lRukRITE8ePH9+rV6/IyMhmzZqdOHEiICCgsJ4PCAiwtLRctGhRXl7euHHjoqKiXFxc5s2bd/fuXQcHB0EQ/Pz8mjRp4uHhkZ2dLbfMxaSmpjo6Onbv3v3mzZtmZmaHDx8OCAhgPU9ERERERB+AJT2VipeX1zfffCMIwqpVq86dOxcbGztjxgxBENatW9e4cWM3N7cXL16UawDZNHVzc/MdO3bIpqlHRET06dNH1vro0aPx48d//vnnsbGxVlZW586d2759e82aNWWtJiYmGzduvHnzpr29fWZm5qJFi5o2bern5yeRSMo1czFSqXT79u3NmjXz8/PT0NBwd3e/devWZ599Js8MRERERERUqQhEbyWVSufMmQNARUVl8+bNRZvu3r1rb28v+4MkK5vz8/PLI8OZM2c++eQT2YPs7Ozi4uIKm/Ly8lavXi3bKF5HR8fT0/PtGYKDg21sbGS3srKyCgkMLI/AJQgLk3z6aff69QH07t07MjJSTs8lIiIiIqLKi2vp6W0kEomTk5O/v7+6uvpvv/02fPjw1/uEhoa6uLicO3cOQIsWLRYvXlxY53+8lJSU+fPnb9q0SRAEMzOzNWvWFB3WPnXq1LRp0+7cuQPAzs5u3bp19evXf+c9BUH4888/582bd//+/c1du07S1MSyZWjduqwyF5eainnz4O8PqTSxf//Ljo6DBw8ur2cREREREVFVwpKe3ki2In337t3a2tr79u3r16/fWzoHBATMnj1btqt8586dvb29C/ef/zBSqfTXX3+dPXt2SkqKlpbW3Llz3dzcNDU1Za2PHz92dXXdsWMHgKZNm65du7Zv377vdf/c3Nztfn6TFi9WSU6GWIyxY7FkCUrxjcB7EATs2AEXFzx9CjU1TJmCJUtQrVpZPoKIiIiIiKowlvRUstzc3FGjRh04cMDAwODQoUNdunR550fy8/N/+eWX77///unTpyKRaPjw4Z6eno0bN/6Ap1+7dm3q1KmXLl0C0Lt377Vr17Zo0ULWVFBQ4Ovru3DhwoyMDG1tbdkGeBoaGh/wFABIS4OXF3x8kJMDdXU4OWHRIvxzsv1HuXEDU6fiwgUA6N4dvr6wtCyD2xIREREREf2DJT2VIDMzc/DgwSdPnqxVq9bx48etra3f67PLly9ftmxZdna2urr6xIkTlyxZUrhT3TulpaV5eHj4+vpKJJK6dev+9NNP48ePL2wtr7Pf4uOxZIlsbjyqV8fcufj2W7z+NUFYGLZuxfnzSEqCSARjY3TrhkmTUOwkv6ws/PADVqxAQQGMjeHpiXHjUOTsACIiIiIiojLBkp6KS0tL++yzz0JDQ42NjQMDAy0/aGz50aNHP/zww+bNmyUSiaGhoaur68yZMwunzZdIEIQdO3a4uLg8ffpUTU1typQpP/74o2zfOwCpqamLFi1au3atVCp9fVF92QgLw9y5CA4GgAYN8MMP/5biBQWYMQMbNkAQoKYGY2MIAhITIZFALIaLC5YufdUzIADTpiE+HqqqmDoVP/wAPb0yzklERERERASAJT0V8+TJk759+968ebNRo0ZBQUEfNm2+0J07d1xdXQ8dOgSgfv3633333VdffSUWl3B0YlZWVt++fS9cuACgV69evr6+hTPtZYvqnZ2dk5OT1dXVZ82a5eHh8fZvBz7KkSOYOxe3bwNAx44ICYGmJpycsHEjdHWxdCnGj39Vpaemws8P7u7Iy4ObG5YuRUoKmjZFWhq6dIGvL6ysyiskERERERERS3oqKi4urk+fPvfu3bOwsAgMDDQxMSmT2wYFBc2dOzc8PByAjY3NsmXLevXq9Xo3e3v7c+fOeXl5jRs3TvTPNPVr165NmzYtNDQUr5X65Ugqxd69cHFB16749VecOoWePaGiguBgfPpp8c5//gl7e4jFuHYNVlbYsQMFBZg4kTPtiYiIiIiovLGkp1ciIyP79u376NGjtm3bHjt2rEaNGmV4c6lUunfv3rlz58bFxQGwtbVdvny51X8HsZ89e6ahoaH3zzT158+fu7u7yxbVm5iYLF26tOiiennIykJ2NoyMMGwY9u3DhAnYurXknv/7H44cwVdfYdMmuSYkIiIiIqKqrYQp0FQFhYWFde/e/dGjRz169AgODi7beh6AWCy2t7ePiopavXq1gYFBUFBQmzZtRowY8eDBg8I+NWvWlNXzgiBs3769efPma9asEYvFM2bMuHv3rrzreQA6OjAyglSKkycBYNSoN/b84gsAOKblk88AAA5ESURBVHFCTsGIiIiIiIgAcJSeAJw9e9bOzi4jI8POzm737t1aWlrl+rjU1NRly5atXr06NzdXW1t7+vTp8+bN09fXl7VGRERMnTr1/PnzAD799NN169Z92P58ZSYuDo0aAcCTJ6hdu+Q+d+/C3BwA0tLK5gA8IiIiIiKiUuAofVV39OjRfv36ZWRkjB49et++feVdzwOoXr26p6dndHT0uHHjsrOzvby8zMzMvLy8UlNTPTw82rVrd/78+Tp16mzbtu3UqVMKrucBpKa+evGWmQtGRsU7ExERERERlT+O0ldpBw4cGDVqVG5urpOTk6+vb4l70Zer0NBQFxeXc+fOAahWrdqLFy9UVFSmTZu2ePHiwnF7BQsLQ9u2AFBQABWVkvukpLyq6u/dg5mZ/LIREREREVHVxlH6qmvHjh329va5ubmurq7r1q2Tfz0PoGPHjmfPng0MDDQyMtLT02vevPnFixd9fHwqSj0P/DuRPi3tjX1KM5JPRERERERU1ljSV1Fr166dMGFCQUGBq6urp6enSKEnrtna2tarVy8hIcHf379du3YKTFICU1Po6ADArVtv7HPzJgCYmHAhPRERERERyRNL+krowYMHAwYMqFWrVseOHS9duvR6By8vr+nTpwNYtWqVp6en3AMWl5eXFxkZKRaLix1rJ5ORkbF3797t27f//fff8s8GVVV06QIABw++sc/+/QDQo4d8EhEREREREcmwpK+EJkyYYG5unpCQ4O7uPnjw4JycnMImQRCcnZ3d3NxUVFS2bNny7bffKjBnoTt37uTm5jZp0qRatWrFmm7evNm8efP9+/dHR0dvUtSp7w4OAODnh5iYElrDw7FrFwA4Ock1FRERERERVXmqig5AZUwQhPPnz2/atElNTW3AgAFGRkbnzp2ztbUFIJFIHB0dN2/erK6uvnPnzqFDhyo67Cvh4eEArK2tX2+aNm2ai4vL7Nmz5R6qiKFD0a8fjh9Hnz7YuRPt2//bdOoURo9GQQHGjEG3boqLSEREREREVRFL+spGJBI1btz47NmzTZs2ffDgQU5OTnx8PIDMzMzWrVvfu3dPR0dn3759ffv2VXTSf12/fh1A69ati13Py8s7f/781q1bjxw5oqqq2qtXL1VVRfyJFYmwcyc+/xznzqFDB1hbw9ISgoAbN3D7NgD873/YuFEBwYiIiIiIqGpjSV8Jbd26ddKkSR4eHhYWFioqKrq6ugC+/PLLe/fuqaqqBgUFdezYUdEZ/0NW0r8+Sp+YmAjgyy+/tLKyCg8PX758+fHjxxWzk5+hIUJCsHkzNm9GWBiuXwcAFRV06YIpUzB6NBRxXgAREREREVVxPJe+MpNIJHXr1j1y5EibNm0sLS3v3LnTvn37ohvmZWVlicViLS0tBYYUBMHQ0DA9Pf3x48d16tQp2pSYmNigQYOMjAxtbe28vDxjY+OgoKDXB/Pl7eVLPHsGsRhGRlDofzoiIiIiIqriOLRYCRXuh+fp6VmnTh1ZDTxu3DgA9+7dy87OlrVu377d1NR0/fr1isopExMTk56ebmxsXKyeB1C7dm1dXd24uDgAYrFYIpEoIN/rtLVhaor69VnPExERERGRYrGkr4Q2bdrUoEGDmjVrnj59+sCBA7KZ6m5ubjY2Nqmpqdu2bZN1q169ekpKysqVK/Py8hSY9i1746moqMycOXPy5Mm//vrrqFGjWrRoUWI3IiIiIiKiqoklfSU0ffr0uLi4J0+enDhxomHDhoXXXV1dASxbtqygoADA//73v9atWyckJOzYsUNRUfHmvfFkFi1aNHPmzCtXrnTq1CkwMFAxC+mJiIiIiIgqJJb0lZNYLFZRUSl2cdiwYc2aNYuNjd29ezcAkUjk7OwMYOnSpQqc037rVhQAKyurN3UYNWqUj4+Ps7Pz66fWK0p2dvaxY8cKFzgQEREREREpBEv6KkQsFru4uADw8vKSbYs4atSoJk2a3L9/f//+/YpKFR6+p1Gj3Nat+ysqwPuSHSUwaNAg2SJ/IiIiIiIiRWFJX7VMmDChQYMGERERR44cAaCiojJ79mwAP/74o0LOPnj6FA8fIjlZ3cxMT/5P/zB9+vSJjIx8fTM/IiIiIiIiOWNJX7WoqanNnDkTwJIlS2RXJk2aZGxsfP369cDAQPnnCQ8HACsrZTrWvUuXLpqamopOQURERERExJK+6nF0dDQyMgoNDT137hwADQ2NGTNmAPD09JR/mOvXAUDhJ80TEREREREpI5b0VY6Ojs7UqVMBLF26VHZl6tSpBgYGISEhFy5ckHMYWUnPk+mIiIiIiIg+AEv6qmjmzJm6urpHjhy5du0aAD09vSlTpgBYtmyZnJPIJt5zlJ6IiIiIiOgDsKSviqpXr/7VV18BWL58uezKt99+q6WldfDgwVu3bsktRlYW/v4bamqwsJDbM8tAWFjYnj17Xr58eezYsfPnzys6DhERERERVV0s6asoZ2dndXX13bt3//333wBq1ao1adIkQRC8vb3lluHGDUilsLCAhobcnlkGHjx4EBYWNnny5CdPnty7d0/RcYiIiIiIqOoSKeToMqoIJk2a9Msvvzg6Om7YsAFAfHy8mZkZgOjo6IYNG8ohwLp1mDYNEyfil1/k8DQiIiIiIqLKhqP0VZerq6tYLN66dWtiYiKA+vXrjxw5Mj8/f8WKFfIJIFtIz73xiIiIiIiIPgxL+qqrefPmQ4YMyc3N9fHxkV2ZP3++WCz29/d/8uSJHALwBDsiIiIiIqKPwZK+SluwYIFIJFq/fn1aWhoAc3NzOzu7nJyctWvXyuHpOTkQi9GqlRweRUREREREVAmxpK/SWrdu3bt37xcvXqxbt052ZcGCBQDWrl2bnp5e3k+/eRPp6TAwKO/nEBERERERVU4s6au6efPmAVi1alVWVhaA9u3bd+/ePT09fePGjWX4lDFjUL06TEzw8OF/ruvqAoCVFapXfzUPH0D//tDRwW+/lXyradOgo4PFi8swHRERERERkVJiSV/V9erVq1OnTikpKb/8s++8rMhfuXJldnZ2WT0lKwtpaXj8GFOmlNCano60NBQUvHqbnY2XL5GfX/KtcnPx8iXy8soqGhERERERkbJiSU9wcXEBsGzZsry8PAD9+vWzsbFJSkratm1b2T5IVRVHjuCvv8r2rkRERERERFUUS3rC4MGDLS0t4+Pjd+3aJbsyd+5cAN7e3gWFQ+dl4auvAGD6dGRmluFdiYiIiIiIqiiW9ASRSCQbqF+6dKlUKgUwbNiwZs2axcTE7Nmzpwwf9OmnsLNDfDw8PMrwrkRERERERFUUS3oCgC+++MLU1PTu3bsBAQEAVFRUZEW+p6enIAhl+KBVq6ChAR+ffzfDIyIiIiIiog/Dkp4AQE1Nbfbs2QB++ukn2ZXx48c3a9asZ8+eOTk5ZfigJk0wYwYKCjB1KqTSt/WMiMDhwyX8xMeXYRwiIiIiIiIlpqroAFRRfPXVV0uWLLl8+XJISEjPnj3V1dXv3LmjoqLyMffMy8PDhzA1/c9FDw/s3o2LF7FpExwd3/jZVauwatXHPJyIiIiIiKiSY0lPr2hra0+fPn3hwoWenp49e/YE8F71fFoaYmKK/zx4AIkEERHFHoTly2Fvj3nzMHQoatYs+YbDhqFNmxKu79+Pq1dLn4uIiIiIiKjSYklP/5oxY8aKFStOnDhx9erVtm3bvt5BEITExMTY2NjY2Ni4uLiUlK8jIurExSE+HiVuja+qClNTpKcXvz58OD77DEeOYOFCrF9fchg7O0ycWML1mBiW9ERERERERABLeipKX19/0qRJq1atGjNmTGhoaMw/EhMTHz9+HBMTc/fu3aysrML+bdp0vHatjuy1oSEaNy7+06ABVFUBYPny4s9avRonT2LTJkyZIqffjoiIiIiIqJJhSU//MX78+FWrVkVHR1evXr3EDnXq1GnYsGHDhg0bNWpkYVGvVi00aoQGDaCh8X4PatoUc+fihx8wa1YZxCYiIiIiIqqCWNLTvwRBsLa2/vrrr7ds2aKnp2dsbGxiYtK4iGbNmlWrVq2sHjdvHn77DcHB+Lg9+IiIiIiIiKoolvQEAOHh4ZMnT5ZIJFlZWd9//72Pj4+WllZ5P1RLC+vWoX9/SCTl/SgiIiIiIqJKiOfSEwC4uLhMnjz5xo0bx44dc3BwEARBPs/t1w+DB8vnUURERERERJUNR+kJACQSSWJiouy1mppamd/fygpZWahTp4Sm1auRlQVBgJ7eqysdOkBLC3Xrlnyrli3Rrx+aNSvzjEREREREREpGJLfxWKrIIiIiOnXq1KpVq8zMzBUrVvTt21fRiYiIiIiIiOgdWNITAIwdO7ZatWrDhg3bsGHDrVu3Lly48KYd74mIiIiIiKiCYElPSE5OrlWrVmpqqoGBAYD27dtPnz593Lhxis5FREREREREb8Pt8QiGhoaGhoYXL14E8PTp0/j4+EaNGik6FBER/b99OzRyEIrCMPpsKkgFiBSAiEZikGkgg49ID8jgSB3MRKcFPCIdMFjMurU7WZG3d+cch3rXfvMPAAA/sNKTUkrP5/N6vS7LklJq2/ZyuXzg0a7r7vf792dd133ff+BdAACA/0HS8yecTqeqqs7nc+5DAAAAwpD05Pd6vY7H4zzPu90u9y0AAABh+Jee/G63W9u2eh4AAOAtVnoyW9e1KIppmvb7fe5bAAAAIrHSk9kwDE3T6HkAAIB3WenJadu2oigej8fhcMh9CwAAQDBWenIax7EsSz0PAADwC1Z6AAAACMlKDwAAACFJegAAAAhJ0gMAAEBIkh4AAABCkvQAAAAQkqQHAACAkCQ9AAAAhCTpAQAAICRJDwAAACFJegAAAAjpCy+ME6utMXCLAAABFnpUWHRyZGtpdFBLTCByZGtpdCAyMDI0LjAzLjEAAHice79v7T0GIBAAYkYGCOABYl4gbmBkY0gAiTNzMCgAaSZBMMXI5qABpJhZEHQGiGaGS7A7WIBoRmYkGTQlcAFuoLXMLBlMLKwJrGwZTGzsCewcCRycGUycXAlc3BlM3MwJXGwJjCwMTMwMIkxAPWzM3FxsrCys7GxcnBzis5DczcDz/1OAw6sI5gMgTvs5KwfHfeL7QewZrXwOSo9FHEDs1pSJ9jeYltmD2KVr99pfrlPbB2KrS5+wc588A6z+4ovW/YwczmD25s8iB0Jjj4LZDUdjD5xvEACzK49bHDhxPRBszgHDTftT5hqB2WE6s/bJCd0Fs8UAhsw+c/uMeMIAAAFwelRYdE1PTCByZGtpdCAyMDI0LjAzLjEAAHicfZPdTsMwDIXv+xR+gUXxT5z4cmsnhNA6CQbvgMQl7y/sjZFOikh6kbRfbJ/jdIIYr8vL5zf8DVqmCSD/85gZfHDOeTpBLOBwfHpeYb7sD/c38/l9vbwBEiD7GZ+P7P5yPt3fIMzAiVQaI+wo5UKcCXLK19GPEpyBkmFrRWGXUzGqUgcgw/zlJJFaFl8wKSsPQPHUmJBLa15HIlPTUeZy5aThb2ZUKjYKqA7mxJYbE+w8tKFZGYDVQf+eRVFrhKyKTXBANljDFabKVeJMEY86Ii1ichIt6HLCIBXfDEjMgVJqrWINoKFLtxGJt0LZCCUUq9WKMiIpSBefEdHcrkoVcZT9uC4PF+B2JQ7ndelXIib1xovP0tvrG9DeRHHfa++VuLutd0TcQuu+i/uEG3cl3EDsJkqIRtp4JSFONpZckbpR7sX2+vzKYxS51btVF/v7z+Lr6Qf7CaiqZ7Df5QAAAPt6VFh0U01JTEVTIHJka2l0IDIwMjQuMDMuMQAAeJxNjr1OBDEMhF+F8k7KRRk78U+2vIYK6BEFiuiARejKe3i8S5GtPN9oPPb1eWCcrp/nMWiM1+/Ht5h4uJ84k1TjdKFcGnFaKDvMJF1Kbk5aNS3I1fBvQag5bxa4mSVkcnHZ1ojEE2Um4eDIcgHgEVFSUFjI7ISaShZXDbHEVTOFhmWIIg+Lc5UG2x+QGnqPMSlriopW4djLSpXYjJgKrIYXJ70Y0xaLULSd0/tt/Xr5XX96yZt8Wm8fuXRMQKcJtfME7nUC9XbYQZcDla6TvNsE6z5BOw5PSMfhi3b/A3aga5m45u6KAAABGXpUWHRyZGtpdFBLTDEgcmRraXQgMjAyNC4wMy4xAAB4nHu/b+09BiAQAGJGBgjgAWJeIG5gZGNIAIkzczAoAGkmQTDFyOaQAaSYmYEMDRCDBUOA3cECRDMyI8lg6IHQ3EBrmVkymFhYE1jZMpjY2BPYORI4ODOYOLkSuLgzmLiZE7jYEhhZGJi4GUSYgFrYmLm52FhZWNnZuDg5xGchuZuB58HSKgfLkFn7QRzjdD+H540MdiB226nJ+/P+WDmA2B8aW+x7WA7bg9je/+/Yb2FwBasRi862bzPpBOtVcBbZv2ZJA5h94dqF/Ym7WA6A2HO97A7Ez9sMFo9mMjgg9mqxLYj9et6i/Wr/l4HNWW3ydl/alNNg88UAeVY/g2xWAxEAAAFselRYdE1PTDEgcmRraXQgMjAyNC4wMy4xAAB4nH1Ty07DMBC85yv2Bxrt084eaVohhJpIUPgHJI78v9hNVZxKFnYOfozHOzPOANneTq9fP/DX+DQMAPjP5+7wKYg4XCAHcDw/vywwX5+O95V5/Viu70AMJHEm+iP26bpe7isEM8joVKQqHGhktOIGOOLW2lGGNYDMhuqxTWzqpYMTmL+TiEwLTcCjMxWrHaTG1TSiGJPFwFis9hhtw1UuppQ3O8vUq7AEDkOKi3lWgJWcesAawAOOVp30hkRD7iEnWHK/CJkVOHBwytQt0pMzxFYR3IxUFK7UQRLeoNWQVBMwIUqPkzKczKSi1ixPuISXPSTfJGkRdk4zHbX0OM/L6SH/24s4rsupvYjs3HLX6NbijQmUlqGG77VFpeHu1BLRsNCb7xo+0c5dTTeImomaool3XmmK050lG6TulEexrb548RRFNooUvJeX8/vPEuPhF6vvqHwocKHWAAAA+3pUWHRTTUlMRVMxIHJka2l0IDIwMjQuMDMuMQAAeJxNj71OBDEMhF+F8k7KWv6JndhbXkMF9IgCRXTAntCV9/Akt0W2Sr4ZZzy5vDZq7XT5Pjdu7f33+aOf9HQ/CTiZlLQQMGpJqwCzoiYEYs1uaSUobJqH4ixVh4KizIlAWbSmdUHIJuwPyTH3nB5Imo0TgzOZPRQuw0MQNrIetDAURRrZWBHF5KF5ERmNMgrvUSaklrqFReq+UItT9jGGqNgfYv+Ki+5S6V3TOX3etp+3v+0aCOP6st2+AIMmUPCEHDJBIh/GKHQShx0sjDLJo06o4RNK0KGEBR1a6P0fzXVrHNOUhfkAAAEZelRYdHJka2l0UEtMMiByZGtpdCAyMDI0LjAzLjEAAHice79v7T0GIBAAYkYGCOABYl4gbmBkY0gAiTNzMCgAaSZBMMXI5qABpJhZ2BwyQDQzkgCEZnewANGMzEgqMJRCBbiB1jKzZDCxsCawsmUwsbEnsHMkcHBmMHFyJXBxZzBxMydwsSUwsjIwMTOIMAH1sDFzc7GxsrCys3FxcojPQnI3A8/WeD0HoZXKB0CcH7IL7HUPCoLZ/GJ8DrO+mTiA2Heuz7M3rbxhD9ages9+89UeOxB74kJHe5XamP0gdshC8/0T7oqA2bOXftmfvGwfmL3a2POAeY4SmC3cqXuA1dMXbM4t2dn7GVaFgtnbffL2/ff4A2aLAQB+rjwIGaQ+jwAAAXN6VFh0TU9MMiByZGtpdCAyMDI0LjAzLjEAAHicfZPBTsMwDIbvfQq/AJHt2HFy3LoJIbRWgsE7IHHk/YW9aaSDiKQHp/3i+P+dThDj5fD88QU/gw/TBID/PK01eM+IOJ0gAtgfH58WmM+7/e3NvL4t51cgBsq+x+c9uzuvp9sbghk4GYs2hgdOqkZNABNeRt/KsAIlVmykAbKZUBmAGeZPT8mM2IoHlU2VBqD40Z5RHPNKU2lEru0vpxfOsirGeZ4OSx1wxTlMpiWbS8FUC3EdSTEH/buRZdSItNI192+ywgIPlLyyjBKR1CaaB2SLnDmRZLMWOUsWKaPTCQN1z4mkhS8VqTYdkdEcP5SJqoTimqXwkOSrJM5s3h6v2LDZyKTjcrjr//VG7Nfl0G9ETO59F5/au+sLKL2H4r5bb5W4u7V3RNzC1n0X94k27kq4QdRNlBBNvPFKQpxsLLkgtlHuxfaC/MZTFLnVu1UX69u/4vH0DUo8qHj54zzBAAAA/XpUWHRTTUlMRVMyIHJka2l0IDIwMjQuMDMuMQAAeJxNjrFOQzEMRX+FsZXSyNexYzsduzABO2JAERvwKtSxH0/SJW+Kz0lyfS+vHZ3ff58/eu/c++Hyfex4uh84G4tGOnFWNUQ6I7NSYBo2E05nyqa1GKcTZa9gl3Qeo8EK6ZTqUKpDIgdQ0jjFQ3SYkiHFfL6qRaTOr2MnIJFGGsFDHx8ZcJ2qSGV9LODCBk0j1Cj80Uwo6hA1MHqNgkxInJ1Ndd5bUcUIGUzV0zF93raft7/t2ijP8WW7fWVqWIDGC7SVBbXJAmu6wFtdEM12adR8R2ixqDTsSnDDroXc/wEBD2tDMyEeSgAAARZ6VFh0cmRraXRQS0wzIHJka2l0IDIwMjQuMDMuMQAAeJx7v2/tPQYgEABiRgYI4AFiXiBuYGRjSACJM3MwKABpJkEwxcjmkAGkmJmRGRogBguMZnewANGMzISVcgOtZWbJYGJhTWBly2BiY09g50jg4Mxg4uRK4OLOYOJmTuBiS2BkZWDiZhBhAmphY+bmYmNlYWVn4+LkEJ+F5G4GHrljNQ5V1mG2IE7BZi2HyPOW+0Hsz19UDxxx5HIAsa+t6LdWesgLZt96vMK+LH2/PYi9TmitvfX9nt0gtvjVPtvzPh5gvZcVwvcdTRY4AGL7R9/c39kgA2bHLRA9sD6iG6zmQdr8/Z7FlXtB7C1TF+/frbEDbKYYAHWeP8TNqTlPAAABbnpUWHRNT0wzIHJka2l0IDIwMjQuMDMuMQAAeJx9U8tOwzAQvOcr9gdq7cuvY5NUCKEmEhT+AYkj/y92G4pTYbBzWDuz650ZewAfz/PT+yf8DJ6HAQD/+Wqt8CaIOJzBAxhPD48LTJfjeNuZ1tfl8gLEQGI5Nu+xx8t6vu0QTCChao2VAQNGFhQPrqNlMqzAIaXMyHDAkDkWxg5QYPqAA4dYBXO0HErISTtItaN9V6Wo4ZhSSamDi4ajIFQ1igXWqmjs4NKGi0wi3qJVZhPrNzBvB6eadeNSWJk7wAKL/2dCLMlJcUxYe6yrlTxQSLWYQY5U/aNLQodyEMnV7LEkLCVpTyCirSpbLc5XSn489aD8Dc0FMblKUrL2jDwt890F2K7EuC5zuxI+uTmvNmPz1xaQmolqwufmlZq8pVmiJmJtwqspRTt91fUgajKqsybeqaVOTneaXCF5x9yabQ2ZpmRNthJOeE/P17fXYvHwBSraqJ/c9JMZAAAA/HpUWHRTTUlMRVMzIHJka2l0IDIwMjQuMDMuMQAAeJxNjr1OQzEMRl+FsZVSy39xbHfswgTsiAFFbMCtUMc+PLm3Q+72nWM7+S6vnTq//z5/9N65Hy7fx97p6X4QCI0aBQErC7KXM4NZ43JCaFydy5lAKpOsBlE5Als5j2zRlHVbdFYemyMyIbqVEwNXwxiOwMJpNapV6jAMIi2GIkB3022Jx/DxxzhDerjmIxGoeNPtsIagFwYyZNta4Bgq6XBM5m2rS6F1PYsqWo7l87b8vP0t10RY48ty+wJMmkDJE2rKBEud0LJO8LQJkW33GqbviDImcdKuhCTtWuj9H+GIbFFIRvOxAAAAAElFTkSuQmCC",
"text/plain": [
"<rdkit.Chem.rdchem.MolBundle at 0x73324e05f610>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"res"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "f21373b7-05bd-46d8-8e88-3bd49b1579af",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.13"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment