Skip to content

Instantly share code, notes, and snippets.

@iwatobipen
Created June 20, 2019 23:14
Show Gist options
  • Save iwatobipen/f0845d1b7b4eca9deab0dc833d8cfd75 to your computer and use it in GitHub Desktop.
Save iwatobipen/f0845d1b7b4eca9deab0dc833d8cfd75 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,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import os\n",
"from rdkit import Chem\n",
"from rdkit.Chem import AllChem\n",
"from rdkit.Chem import DataStructs\n",
"import numpy as np\n",
"from rdkit.Chem.Draw import IPythonConsole\n",
"from rdkit.Chem import Draw\n",
"import seaborn as sns"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"from rdkit.DataManip.Metric import GetTanimotoDistMat\n",
"from rdkit.DataManip.Metric import GetTanimotoSimMat"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"from rdkit import rdBase\n",
"from rdkit.Chem import RDConfig"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"2019.03.2\n"
]
}
],
"source": [
"print(rdBase.rdkitVersion)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"sdfdir = os.path.join(RDConfig.RDDocsDir, 'Book/data/cdk2.sdf')"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [],
"source": [
"mols = [m for m in Chem.SDMolSupplier(sdfdir)]"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA+gAAADICAIAAAD0hVwYAABfLElEQVR4nO3deVzM+R8H8NdM06VSITkq973oIFeOyK5dsdjNHjaWJSybPWXXkauf2EOuZbBL2LWb3WVjc4Rl3ZSQIuQootCdrpn5/P74jDEq6ZjvHPV+PvaxDzNNn++nmu933t/P5/15f0SMMRBCCCGEEEL0m1jXHSCEEEIIIYS8HAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAAwscF+zBg4OyMkBgKAgbNmC9esxf77yq/PmYeNG3XWOEEIIIYQQwRhY4A7AyAhz5ij/zdiz/6v/gxBCCCGEkBpGousOVNr06di+HdHRz57Ztw+ZmQBw6hT8/HTULUIIqd3WrMGSJbhyBVZWCAqCoyMKCpCSopwUnTcPTk6YOFHHnSSkFqJzsyYxvBF3iQTr1mHqVMjlymc6dcLo0Rg9Gh076rRnhBBSu9GMKCH6ic7NGsPwRtwBuLujZ09s3IhFiwDA0RF9+gDAvn267RchhNRqNCNKiH6ic7PGMLwRd5EIAIKCnntY4t8HDiAhQbvdIoSQWo9mRAnRT3Ru1hgGFrhPm4ZPPwUAKyskJWHcOEyahNxcDB6MmzexcCEmTsTMmcjIwMyZSEvTcW8JIaS2Uc2IcnxGtE8fODrqtFuE1Hp0btYMBha4l+n4cRw8iIcPlQ+nTMGrr0ImQ506Ou0WIYTUMjQjSoh+KvPcLCh49m/QuWkgRMzwVyV4euLIERw+DE9PADh1Cj/+iAUL0LKlrntGCCG1VW4uvvwS3bo9V61i5ky4ueHXX7FhAxo21F3nCKn1+vXDqVNITISTE0DnpuGoCSPufGT9yRPlw3nzYG6O779HcrIOO0UIIbXapUuQSvHjj889STOihOgJExPIZLh0SfmQzk1DYZBVZUooEbhHRuqwL4QQQgAgNhYAunR57snUVAQGYtUqWFrqpFOEEKUuXXDoEGJj4e0N0LlpOGrgiDshhBCd44F7587PPUkzooToCX5u8vMUdG4ajho44k4IIUTn+BR8iRF3mhElRE/wc1OVKkPnpqGgEXdCCCEaxljZqTKEED3RqRMkEiQkKGvLEENBI+6EkDLcv38/NTXV2dlZ1x0hBikpCZmZsLODvb2uu0JqAZlMlpCQEBMTExMTc+HCha5du3p6eg4bNkzX/dJrZmZo3RpXr+LqVdCV3oDUhMDdzu6Ru7vMxKQIcNJ1XwgxPEVFRXfv3o2Li4uPj79582ZcXNzly5ezsrLs7e1fffXVZcuWNWrUSNd9JAaGz79TNEAEkp+Py5cvnz9/ggfrsbGx+fn5qq+ePXs2JCRk165dw4cP12En9V+XLrh6FZcu0alqSGpC4K5QbDt79rNevWYAIbruCyF6jTEkJeHaNSQk4OpVXLsGC4tDf/89uPR+DvXr15fL5Vu3bk1PT9+9e7dIfSsdQl6mzAR3QqosKwuxsYiORnQ04uMRG4vu3U+cODFF9YLGjRu7PXXy5Mng4OAPP/wwOjq6RYsWOuy2nuvcGWFhz9anEoNQEwL3OnXqAHhCuTKEAADWrMGSJbhyBVZWCAqCoyNSUxEVhYQEXLsGtWEpAOjdu6VEInF0dGzZsmXHjh07derU8qmUlJQuXbr8888/69evnzx5so5+GmKQMjKO9unj5ObWALDSdV+IwShx7XJwQFISzp/HhQu4ffu5Vxobw8mpR+vW41xcXJydnZ2dna2trVVf9fb2vnr16q5du955553jx4+bmJho+QcxFN263ezefWtSUgaNexoQCtwJqYGMjDBnDlasAADG8PffOHFC+aXGjdGuHdq1Q9u2aN8e7ds7Ojk9kUjKuBQ0adJkw4YNo0aN+uKLLzw9Pdu2bavFn4AYtoiIqVeuXFm16jzgouu+EEOifu0C8PPPypDdxAStW8PN7dl/5ubOwOYyGxGJRD///PPFixfPnTs3a9asH374QSt9Nzzt24vPnZufnNyIAncDIio9RW5w/vrrr7feemvkyJF//fWXrvtCiO6tWYMnT7B9OzZswL59cHBA3bp48kQZrNetW7nWPvjgg19++aVXr17Hjh0zMjISpsukRiksLLS0tASQk5NjZmam6+4Qg1H62iUSwcQEzs5o2xbicsvg3bp1KyYmprCw8L333uPPnD17tm/fvsXFxX/++efIkSO18QMYGsaYjY1NdnZ2ampqw4YNdd0dUiE04k5IDSSRYN06TJ2q3BKv/M+sEotTb968aWZmtnv3bv7VNWvWHDt27NSpU0uWLJkzZ47wfScG7/LlyzKZrFOnThS1k8oqce0aO7bsl8lkuHoVsbH/nDt36MKFCzExMZmZmQBatWqlCtzd3d2XLFnyxRdfTJgwwdnZmZLdSxOJRK+88srJkydjY2MHDRqk6+6QCqHAnZCayd0dPXti40YsWvTc8xkZiItDfDxu3sTNm4iLg6npVzExK9VfY2VlxRjjC1Ktra23bds2YMCABQsWvPbaa927d9fmT0EM0aVLlwB0oaWppErKvHYVFuLCBcTEKP+LjUVBAXr2jD99ejl/QePGjZ2dnV1dXVXXLgCfffbZf//99/fff8+Zc3TTphaU615aly5dKHA3LBS4E1IDiURYsgTNmpV8/p9/lONY6jw9+7Rrt799+/Zt27Zt95R6GZm+ffv6+/uHhIR8++3d0NDu5uYC954YuNjYWACd+Y7qhFSGSIR9++DlhZ07n3s+Jga9ej33stat0atXv+HD/+fs7Ozi4lJm1VqRSLRp06YJE65v3+5uZ4eQEGE7b4j4eRpLlWUMBwXuhNQ006ahqAimpjA2RlHRc19q3x4NG6J9e2W+O1+l2qLFaGPj0eW3uWTJkqSkT3fsaNakCX34kZegEXdSNdOmAYCtLTIzkZEBG5tnX+rSBS4u6NoVzs5wcYGzM1+u0wPoUX6btra233zjHhGBlSvRrx9GjRKu+waJn6f8nCUGweAXpzLGfvrpp0mTJllaWm7ZsmX48OG0fo6Q1FQ0agQ7O6SlaazNmBj07IniYkREYMgQjTVLah57e/u0tLQ7d+44OdGmeKTSLCzw5Any8pTbomvE8uX4/HPY2CA6Gi1baqzZGiArK8vW1tbMzCwnJ4fCJ4NQ7iJtvXfw4MEePXpMmjRJIpHk5uaOGjXK0dFx1qxZ9+7d03XXCNGlrCwAz41XVZ+LCwIDwRgmTkR6uiZbJjXJrVu30tLSbGxsHB0ddd0XYpD4PKFm89E//RSjRiEzE++8g8JCTbZs6KytrZs0aZKfn79t27aCggJdd4e8nKEG7idOnBg4cODgwYPPnTtnZ2c3d+7c1atXd+rU6f79+0uXLm3ZsuXo0aMPHjyo624SohuZmQCgtiGJZsyahX79cO8ePvlEwy0Tg5aXl3f8+PGlS5cOHjy4Q4cO9vb2CoXi9OnTuu4XMTwKBWQyiMUoa2+JqhOJ8NNPaNECUVGYOVOTLRs0hUKxYsWKBw8eODg4fPjhh1ZWVt26dZs1a9bBgweLSuRZEv3BDM3p06e9ny6vq1+/fmBgYHZ2Nv+SQqE4cODAm2++qZrucXV13bTp1ydPdNtlQrRt/34GMC8vzbd88yazsmIA275d840TA5KenhUeHv7ll1+6u7urb+AlFoubNm0KwMTE5Pvvv1coFLruKTEk+fkMYGZmgjR+9iwzNWUA++MPQdo3LLdv3+7bty8AkUg0YMAAZ2dnsVqp/Lp16w4bNiwkJCQ29gadxHrFkAL36Gj20Uez+FvKxsZm4cKFqpC9hHv37gUGBtrZ2QHo23ertTXz82Px8VruLyE68/vvDGBvvy1I4+vXM4DZ2LCkJEHaJ3orLY2Fh7OAAObmxnr0CFd9xhsZGbm5ufn7+4eFhT1+/Li4uDgwMJAHAYMHD37w4IGuO04MRlYWA1jdukK1v2IFA5irK5PLhTqEQQgLC7O1tQVgb2//zz//8CcfPnwYHh4eEBDg5uamKizWv3+wnR3z8WFSKbt8Wbe9JowZSuAeF8d8fJhIxFxcDlpYWAQEBKSnp6u/oLi4uPR3FRQUbNmyZeTITIABTCxmQ4eyiIjafrqS2kAqZQCbOFGo9r29GcDGjxeqfaIrq1ezpk0ZHxJZvJiFhrJ799ivv7IpU1jHjkwkYvxyCrCGDdP79u07e/bs/fv35+TklG4qPDy8fv36ABwcHI4ePartn4QYpocPGcAaNBCqfYWCff89y8wUqn39l5WV5evry4PykSNHPnr0qMyXJScnb968eezYsX37XlCd9QBr3pxNmMC2bWMpKVruOFHS98A9IYG9/z4TixnA6tRhX37JHj7MUH/BqVOnvL29v/jii3IaiY5mfn6sTh3l265VKxYczF7wXiWkJli2jAGs3NOiWh48YP7+7AUzXsSArV7NnJyYvz9jjC1ezObMYeqf2XXqsEGD2IIF7MgRVpEUxOTkZD4Xb2RkFBgYKJPJhO4/MXT37jGANWkiSOOl70ulUhYYqPzq3LlswwZBjqs//vvvWLNmzQBYWlpu3Lixgt919Sr78Uf29tusfv3nLgheXiV/n0QL9HdxalISJk9Gp0749VdIJPDzw/Xr+PZbNGhgw18QHR09dOjQXr167dmzZ/v27YUvXiju6gqpFElJWLoUzZsjMRGzZsHJCdOmwcEBOTkAEBSELVu08oMRIjxeVUbji1O5NWvg5obFi2FlhaAgTJiA+fOVX5o3Dxs3CnJQojXTp+PYMURHA0Dz5mjQAF5eCAxEZCTS03HwIObNQ//+KGcfrpMnT8rlcgAODg6HDx/+5ptvGGMLFiyYOjXm/n1t/Rg68vDhw1WrVvn6+q5atSovL0/X3TE8QpSUUWdkhDlzlP/m1bBVNbENvDj2S8hkmD8fAQEt79y507179+jo6I8++oh/6cqVK+vXry/ne9u1w9Sp2LEDjx4hMRFSKXx8YG2NRo3K+H0Soelj4M4Ypk9HmzZYvx4iESZNwvXrkErRpInyBZcvXxs1alT37t0jIiKsrKzmzp0bFxdnampafrP162PmTCQmIjISPj4A4OhI7zlSMxUVLW3bdnSjRv8I1H6JE6eWfPLVEhIJ1q3D1KmQyyGRIC0NkZGYPx9eXnjZVRYAoqKiPD09Bw4cyMvySiSSoKCg/fv3Dx68YuPGbs7OOHBA8B9BV9LS0gYOHOjv779t2zZ/f38nJ6fZs2c/ePBA1/0yJMXF93v1Wu7qukOg9tXvS7l9+zBjBmbMwP79Ah1T965dQ69eWLAAZ882+e67f0+ePNm2bVsAjLG1a9d269bt448/PnnyZEWaatkSfn4IC8PDh3BzK+P3SYSmF4H7mjXPjXxv3YqsLMhk8PFBXBzWr4dqE49btzB5Mt59N2fXrl3m5ub+/v7Xr19fuHChTYULVovF8PJCWBhu34aVFb3nSM2UkhJ77doOMzOhyq2XOHFqwydfreLujp49lZMnT5eoVVRBQUGDBg3+++8/Z2fnf/5R3jp6eXn98ov/a68hLQ1DhmDGjJJ7+tYAaWlpgwYNunz5cvv27UNDQ728vNLT0//3v/81a9Zs7Nixly9f1nUHDUNh4eNTpz6/dm2hQO2r35dynTph9GiMHo2OHQU6po5t2QI3N0RFoVkzHD6ML74YwMtAPXz4cMSIER9//PGTJ09GjRrVsZI/v7ExjI3L+H0SoelF4I5SA3hBQYiLQ1gY2rRRPnnrFj78UDkMf+OG27x50sTExBUrVtjb21ftiHZ2QFnnMCE1QGZmJoCK39BWVokTp8Z/8tUqPFIPCqrit3t4eFy8eHHo0KGPHj0aNmzYjBkzeEFoOztERCAkBBIJVq6Ehwdu3tRcp3UtNTV14MCBly9f7tChw7///jt27NjIyMjo6GhfX1+FQrF169bOnTt7eHjs3r1b1z3Vd/zdYiJcrszz96UAHB3Rpw/69EHN2zHs0SOMHIlx45CbCx8fxMSgXz/llyIiIjp37hweHl6vXr0dO3aEhYVV+fOixO+TCE1fAvcSA3hOTmjfXvnvx48xZQratUNoKMRiTJ6M69cxf/6kRo0aaeTQ9J4jNU9WVhYAa4GS3AE8f+LUsE++8+fPe3p6hoSELF26VNd90baxY7FjB+bNg5UVkpIwblxVGmnQoMHu3buXLVsmkUhWrlw5adLZW7cAQCTCjBk4fhwtWuDcObi6IixMs93XDR61x8XFde3a9b///lN9Nrm6um7ZsuXatWv+/v4WFhYnTpwYPny4i4vLli1biouLddtnvSV04C4SIS7u2X2p+oRSZSeX9NzBg3B2xq5dqFsXW7ciLAy2tgBQUIAZM/Dtt3X5+/bixYtvv/12lY/Cf2lvv42sLMTFaajrhqOgoGDOnDmzZs367LPPkpKStHRUXa+OZYyx1avZDz+wM2dY9+5swQK2efNzX83NZQ0bMrGY+fiw69c1fNzlyxljLDubOTqWPC4hhuuVV14BcPHiRSEaX72aBQUxmUx54owfz+bNU37J0MsypKWlTZw4kRcg57PJH330UVFRka77pT0REQxgvXtrprVTp069/vq3YjGztmY7djx7/tEjNnw4A5hI9Kymh4G6f/8+zzFwdnZ++PDhi16WmZkZEhLCt6YC0Lhx48DAwIyMDC321DAcOXIEQL9+/QRq/7//mETCJk5kNXtTofx85ujIANa/P7tz59nz58+zjh0ZwMzN2erVf8g1VB6b1zHz89NIYwYjNzd30KBBAIyNjfn/P/jgA4E+dtXpUeDOGPvkk7ID6H/+YVevCnX0TZuYgwObM0eo9gnRPgcHBwB31C7Y+fn5mmo8P5+5u7PXXmM1KeooLi6WSqUNGjTg11++vrBOnToABg4cWHsCrIAABrBvvtFYg1lZ7J13lMXjfH1ZXt6zL0mlzNycTZ9uwBXlkpOTW7duzaP2F9XDVldYWBgaGqpKJq5bt66/v38S7WSm5sCBAwC8nt/2OTc3VyONp6Wxpk0ZwGbP1kh7eqR0pctDh9iSJc82rlEoWEiIcuPYDh3Y+fOaPPrZswxgbdtqsk09l5OT4+npCaBRo0Z//vmnr68vD98B9OnTJzw8XLhD60vgvnw5KyxkO3aw+vW1PfK9aRMD2AcfaPWgeuvSpUu0RXkNYGlpCSArK0v1jJubm5eXl0auJuPHK7fhePHwooE5cuRIly5d+DWXry/kz1+4cIHfArVu3TohIUG3ndSOHj0YwPbv13CzoaHKnTQ6dmSxsc+ev3+/ZOV4A5r5TEpK4lG7i4tLRaJ2FYVCERkZ6e3tzd9yJiYmvr6+ly5dEq6rBmTPnj0A3njjDdUzV65cMTc39/HxOXHiRHValsvZa68pB6HL2rPRsJV/Hql+dpGITZtWoU0YKkUmYzY2DGB372q4Zf2Unp7eo0cPAE5OTteuXeNP3r9/PzAwUJWh6uLiEhoaWub2oNWkF4E7l52t3OBDyyIjGcA8PbV9XD2UmJhoZmbWr1+/goICXfUhPT399OnT7du3Hzdu3IYNG+Lj4+lGorKKi4tFIpFYLFb96hISEsyflt3u3r3777//XuWtcJYuZQCztGQ1I8y4e/eur68v39zb0dExtNR47927d11dXQHUr1+/xm//mZPDjI2ZRKIct1MoWLt2bOTI54bJqywmhrVrxwAmlT73/OrVbNky5uLCoqIMKXC/c+dOq1atALi6uj5+/LhqjZw/f97X15cnZakG6mr5Fe+vv/4CMGLECNUzq1atEj1NP3d3d//555+fVCnwnD+f7/jL7t3TXHf1xkvPo0WLWMOGTIMDwX//zT74gJ08qXw4dCgD2C+/aKx9vZWamtq1a1cALVq0SExMLPHVrKyskJCQJk/rl7do0SIkJCRPI9fQp/QocGeM2doyQPBhvOxsFhjIZs5UPoyPZwBr00bYgxqE4cOHAxg3bpwO+zBz5kzR80uE6tat6+XlFRgYGBkZqcF8jxrs0aNHAGxtbdWffPjwYXBwsOpq0rx58+Dg4MxKbvy9fz8zMmIiEQsL02iPdaGoqCgkJMTKygpAnTp1AgMDX/TuysnJGTZsGABTU9OtW7dquZ/atHcvA1ivXsqHFy8ygDk6aqz9nBy2Zk3JJ8tf46Sf7ty507JlSwBubm5VjtpVbt26FRAQoBqo69q1a2hoaK1aWaHut99+AzB69Gj1J69duxYQEMAz2QBYW1v7+fnFqs/dvMy//zIjIyYWswMHNN1j/fDS80gm03Bw9dlnDHi2RuXbbxnAJk3S5CH0UEpKSqdOnQC0b9/+7ovnFwoKCkJDQ9s/rbJiZ2cXGBhYqXm5cuhX4O7szAAWFSXsUfLzmUjEjI2VuV+6GunXN/v37wdgZWWVkpKiqz7cvXu3Tp06IpHo119/Xbly5TvvvKNay8WZm5vP9vFhs2ezf/5hlQw6a4+1a9c2bNjQ1NS0dPosv5p06NBBPcW2nKuPuqtXlZOhixYJ028tioyMVF1Svb29b926Vf7rZTKZv78/AJFIFBgYKNyYaE5OTkhIiJeX17fffnv79m2BjvIis2YxgH39tfLhypUMYGPHCnvQl65x0je3b99u0aIFgG7duqWnp2uq2cePHwcFBakq0jg6Op47d05TjRuQGTNmAGjSpEl4eHiJicGCgoKwsDAvLy/VJ4Kbm1tFbnIePGCNGzOAzZ8vZNd1Svvn0d9/K/OOuHPnav4Y6O3bt/k8m7Ozc1pa2ktfL5fLw8PDe/bsyd+ulpaW/v7+6mvPqka/Avc332QA+/NPwQ9Uvz4DWGqq8qGVFQOY5q7AhqewsJDvo/YDP/V1ZMKECQDeffdd9Sfv3bsXFhbm7+/v5uYmFovD+vVTrnQDWMuWzNeXSaXsaVKyUnY2q5XrvRQKxcKFC9WnLExMTD788MMSQ1P8aqL6/OMptpdL/A6fl5WlrEUwapRhV2O4fv360KFD+Q/erl27ffv2Vfx7pVIpz2p45513ND75k5qa+s0336iXUm7duvX9+/c1e5Ty9ezJAKb6lbz1FgPYTz8Je1DDqu6litp79+6tvoZEU/jq1VdeecXc3NzT07O2rVvduHGjkZGRKq/Pyclp8eLFpc+C+Ph4f39/vpIHgL29fUBAwItuv2Uy+YABDGBeXkxDNVT0keo8unmT1a/P5s4V/IiZmczIiJmaKjPmdZjmfvPmzblz5/br18/f3/+eYIlQV69e5UueunXrVtmx82PHjnl7e/OPZrFY7O3tXZ3bcv0K3D/5hAFMC6Fjly4MYNHRyocdOjCAVWbaraZZvHgxgI4dO+pwfvbq1asSicTY2Pj6i6t+pqenp+3bx2bNYh4ezMzsWQTPZ/Tff5/t3csmTGB2dqx7d+bgwE6f1uaPoFuFhYVjx44FYGRkFBwcXJH02aioKNVrRCLRi1avymSyESNyAebszDRU3UEH8vLyAgMDTU1NAdjY2AQHBxcWFla2kX379tWtW5fHbRUZcamIW7du+fv78wo2/C/1+++/d+/eHcArr7yiqdnVlyqd4N6wIQNYqRxOzRs4kHl7G8C99q1bt5o3b87/Rtn81yQMhULx+uuvA1i2bJlwR9E3ISEhPLL59NNPpVIpT0jgFzRvb+/IyMgSl6/s7GypVMqzjXk85OXlFRYWVmKcfvbs2T17jmrTJlN3c8la9d13DGATJmjjWK6uDGCHDysfenszgG3bpo1DqyQlJfF7adW0/NSpU2/cuKHZo8THx/NE0759+1b5jv3SpUsaKT6jX4E7f8N9+qngB3rjDQY8W6Xh5cUAtnev4MfVT8nJyRYWFgAOHTqkw268+eabAKZNm1bRbyguZlFRLCSE+fiwBg2U4fvUqczTk/HR0D//ZAMGCNdhvfL48eMBAwbwybi///5b9XxF0mevX78+depUVeDYrVu3mzdvqr/gyy+/bNOmc+fON7UQwwkkPDzcycmJ35/4+vo+ePCgyk1dunSJN9WqVasrV65Up1f8Oq66cfL29lbVzcjIyHB2dgbg7u4uaIyosn//oS5dPh479r+nfdNwgvuLZGQwIyNmYqLv94TXrl3jmXseHh5a+Ivs3LkTgIuLi9AH0gcKhSIgIICfBcv5uDFjjLFjx475+PioAp02bdoEBweXvpWNiory8/Pj4/SWlpbqcdXevXvFYrFEIjl69Lh2fhadu3CBAczJSRvH+vzz59LcV6/O7ts37uuvd2rj2IwxxlJTU3nmZ48ePY4dO+br62tkZKQa1T516pRGjhIdHc3XV3h6eubk5FSzteoXn9GvwH3HDgawkSMFP9CkSQxgP/6ofDhuHAPYxo2CH1c/+fj48Nl/Hfbh9OnTIpHIwsKiirkBCgW7fJlJpeyDD9j69conZTLWpAnTStCjWzdu3GjXrh1PDI0qa40IX+de/uYvfIOYJk2a2Nvbq+eBbN26FYCxsbGB1lSJiYnp27evKiP2pKoIQjWkpKS4ubkBqFev3r///luFFtRnTo2NjctMVUpNTeV/Vk9PTy0sy/76668BBAQE8IcbN0YMGBD5+eeC5+r8+ScD9P0WOyEhgZ8+ffv2rf4nd0UUFhbWq1cPQFxcnBYOp0MymWzSpEk8Z++3334r/YL79+8HBwc3a9aMn8WmpqY+Pj6RkZElXpaWlhYcHLxw4ULVM8nJyXZ2dgD+97//Cfsz6BOFgtnbM4A9rVIooPBw1qJF4fjxyhJjUVFRPMdP8AMzxhh7+PAh32qwa9euqjXiN27c8Pf3V2VbVb9S09mzZ/mZOHToUA1ehx8/frx48eKGDRuqEiPL/Owuk34F7mfOMIC5ugp+IF4Wat485S3Ot99eHzDg35AQg4xLqunQoUMA6tSpo/2VcOr4aPE81Q6cVTZ1Ktu+XflvmYw5ObGaXpzhxIkT/MOpS5cu5WfE8vRZ1QS0lZVV6dWrT548iYmJUT2Mjo7mI/Hr1q0TqP9Cc3Fx4VmwP//8swYXlebm5vI5IhMTk9JFJF9EoVCEh4f37t2b/wksLCzKX6t0584dPrr/6quvViGxp1J4ryIiIvhDvgv6RuHHMz7+WN9XPF+9epXPkvfr1087UTvH1/wEGvrWsuUqKCjg7zQLC4u95c56y+XyyMhIHx8fPqQKwNXVVSqVvugvUlxc7OHhAeD111/X1P6ghuLdd58bmhROZmamkZGRiYkJL3cok8n4Kp3qr7+syKH56Ennzp1L71icmpoaGBhoa2vL3ypdunSpWqWmI0eO8Mpjo0ePFiKRWFV8xsTEZMyYMRVc7K5fgfuDB3xyVvC1b5s2/WFuXu/DDz/kD3/88UcAk2p8HaNSiouLO3fuDGDJkiU67EZ4eDiABg0aaGCx14EDbOhQZd3prVvZ8zXFdEDgZbK//fabmZkZgCFDhlRw+r7E5i/GxsY+Pj5lLpS5f/8+X4tTifwlPRMfH29lZWVpaSnEiiWFQhEYGMh/jf7+/uUHB0VFRep3TdbW1gEBARXJX09ISLC3twfw1ltvVbn6/kvl5eWZmJhIJBJ+DioUCj4UpPFU0dJ4cXdNTIQIQhW19+/fX5tRO2MsMjKSZ2Rp86DalJOTM3jwYAC2trYV31zpxo0bAQEBfLQCQN26df38/ErvM//VV18BcHBwKB3V1XgbNjCAvfWWNo7Fo2dVni2vnCt02dzc3Fx+V9amTZty6uDxIl2Ojo78rdKsWbNKFUGOiIjgI/djxowRYh8lFblczncAPHjwYEVer1+Bu0LBmjTpDEDo6+PevXuhtqny33//jee3aqsNHj9+PHLkSCsrq3r16glRHqGC5HI5X10UEhKigeaKiti77zIbG9a1K2vfnl25wq5f181mbsXFgi6T5VEjz7Xw9/evQkhX/urV/Px8vjOch4eH0GO9wpHL5fynmDVrlkCH2LBhA8/B9fHxKXNfmIKCAqlUygfOATRq1CgwMLBSFfSjo6Otra2HdOuWP326QDV9+D7z7u7u/OHly5d50CPEsdTdvcsAZmWlpxNjV65cady4MZ/xqNqmP9Uhl8v5PUPF59A1TMhxh8ePH/MyeY0bNy4ddr+UqjSkqoiWm5ubVCrlf6Y9e/aIRCKJRHL8eG1JbVd3+zYDmI0NE+xO/5kvvvhCfbb8+++/B/DRRx8Jd8S8vLz+/fsDcHJyqkimQImpZl4E+aVDOeHh4bySweTJk7UwY/PJJ58ACA4OrsiL9StwZ4zxooRCZ/VdunQJQIcOHfjD6OhoniYl6EFLk0qlERERK1asEPRmrkx//PEHLxjMx2s7d+6sniChTT///DOA5s2ba3K71sxMZbC+bx+ztGSDBumgDNju3cItky0oKBgzZgwAIyOjVatWVaepF61e/eijj/gQhaZqp+hKVFSUWCw2MTFJSEgQ6BAHDhzgv8CePXumqqrMPl1aoKrM3bp165CQkKplSV46caKQ11r76ivNdfyZ2bNnA5j5dF+61atXA/jggw+EOJa6TZsYwIYNE/o4VXT48GELCwtPT0/tR+0c3z3giy++0PaBBR53uH37Nl+/0bJly2rO6ly8eHHq1Kk8n4EnxU2bNo3nSHz//fea6rDBad2aAezsWcEPxCfM+/Xrxx/yaErAaaLCwsJhw0b06OHg4FCiiEL5eJqiqgiyqampr6/v1atXy3zx9u3b+ZDW9OnTtbOT8ebNm/noT0VerHeBO584Kz/XrfrS09MBWFlZ8YcPHjwAUL9+fUEPqq6goGDcuHG8DwA6dux4QFv7uT148ICvRsXTwnP8AiqRSAICArRcDjI/P5+PRAo1s5aWxho10lKR0RKmTi25TPbmTebuzoKD2YvrXVbEo0eP+GpLS0vLPXv2aKCrpTZ/4WtxLC0tqzASpofGjx/PB02FO0RsbCxfPNeiRYv4+PgSdQNcXV1DQ0Orm+gSGclMTRnAgoI01Otn+vTpA+Cff/7hD/klYsOGDRo/UAkffMAAppHJNiHwkv+rV6/WVQdOnToFoEmTJsJlSZVNyHGH+Ph4nr3g6uqqfqNbHbw0JC/ExAekhgwZop2QSz9NnswApoUc2BJp7nK5nN81CZLmXljIhg5lQKGjY1JVP0Z5EeRyis+sX79eLBZDbaW+FsTGxvL72Iq8WO8C94kTJ0IrK+F4AUSeFiyXy01MTABoZ1jl4cOH/fr1A2BhYTFv3ryOHTvya423t7fQGaVhYWG8qlGdOnWCg4P5BNCTJ08CAgL4O9Xd3f1F96BC+PbbbwF06dJFwKmo8HAGpPXoEXvpklCHKFPpZbJr1z6rOt+lC1uwoOS+URVw7dq1Nm3aAGjatKnGJ0lUU4omJib16tWTSqWabV9XUlNT+ZIpTd3nlOn+/fu8+HqdOnV48oxIJHr99dePHDmisWPs3MkkEgYwjQ4l5uXlmZqaGhkZqQoNzZgxo2nTpteEL0vh5KTXe2hs2rQJwKBBg3TYh9atWwPQ5LuoIkqPO2Rns4iI6u9Xfe7cOf4Z1L9//0oljFXQiRMn+G9MM7mXBuuPP2Rdu2b6+WnjbdOtWzcAqlUKw4cPB7BlyxYNH0YmY++8wwDWoEEVPjpLeFHxmR9//JHHQurlibRAJpNZWFiIRKKKLHzSu8B94cKFAL5WbbotGB79xMfH84d8tEwLK7EuXLjAj+Xg4MAzF4uKikJCQvjgnLGxsb+/vxAZ5/fu3eNFMAAMGTKk9N3wsWPH+F6+5ubmqpheUJmZmfXr19fCBMuRuXPFIlHXrl01mY1Tvp9+Yrt3l1wmm5/PwsOZn59yYxv+X4sWzN+fHTtWkdzlgwcP8gDU3d1duD01FQpFr169APz1118CHUL7li9fzidwBa2rmJ+f361bN3t7ez6Qc1aIierQUCYWM5FII/Vr5XL5uXPnePWSTp06Vb/BSomLi7Oyavraayv1dmA0IyOD39JUp/B/NfEynZMnT9bqUUuPO9y6xQBmZMTc3FhgIIuKqsJyi8OHD/NJ5mHDhgk3TLZmzRqU2oG7tnn06JFYLDYzM9PCcOTFixfVZ05++OEHAMOHD9dkDrBczsaMYQCztmaaW/KRkpKininKS76KRKKVK1dq6hAVx+t6la5zWpreBe6hoaEA3n//faEPFBgY6O/vr1rZwGeKha5UHRYWxovr9e7du8QnwaNHj/z9/fn0TePGjaVSqaZCZ4VCIZVK+eXSxsZGKpW+aAIxKyvLz8+Pv4MHDRokdEUn/oGkyo0TTm5uLl87oY1UUYWCffklA9iQISWXyaoUFbF9+9ikSczOThm+i8XDunf/6quvTp8+/aK/zk8//cTHcUeNGsUnJYUzb948AF8Jk06tE6oCSkJXdOaT9cLOGa5apQyhwsKq1kBiYqJUKvXx8eF3zvwTy9jYeNeuXZrtaflWrlwJYMyYMdo8aGXx+ks/aqG03gvExcXx7DWtrhEvXZ4rJob178+MjZ8NOjRuzCZMYGFh7PkdIV5k586dfEmVr6+voMu6rl69CqBhw4a1OVWGPa33UpFAULNmzpxpa2vLN2bx9vaWSqXl1yl+OYWC+fkxgNWty86c0VA3n8nOzubFZywtLc3NzXV1svMFLRUp8ad3gfuRI0cAeHh4qD8ZFxfXs2fPipeLqoLRo0cD+OWXXwRqX6FQBAcH8ymYMWPGvGjYLyoqit9CAOjWrVv1f+QbN254enqqUnHuVqC+SkREBC9lULduXeGSJVJSUvg9jEY2xHmp33//XSQStWrVysvLS8ASDQoFmzGDAczEhP3xB2Nqy2TLJJOxw4fZtGmPnv6NADg6Ovr7+x89elR151apsoMaERERUfo0NHRa2LLgzJkzABo0aCD43M6CBQxgpqasYuXDGGPs0aODO3f6+fm1bNkSalq1ajV58mRVTfo/+PtWK/hBf/75Z60dsQr4WJKnp6cO+8B3mVGtQBDc99+zhISyxx2ystiff7JJk5iDw7MIXiL5bdy4oKCg6OjoF8XKmzZt4qv9tHMF43VsL2k5PVLP8M1ohSuoVZpCoeBVOMVisaoII9elS5eZM2ceOnSo0vefCoVyr4c6dZiQCWOFhYV8cEfL4xcq/FLz9ttvv/SVehe437p1i8cu6k++//77/K0wfvx4IaYsMzIyOnXqZG1t3aJFi8ruPVshOTl/fvIJAIlE8tLEO4VCERYWptqe3cfHp2p3q8XFxSEhITyV397efseOHRX/3rS0tFGjRvHz7fXXXxeiBjYf2q/Ie7T6wsPD+YSDKvPYx8dH85WL5HL20UfKcKqSZ75cLo+KigoMDOQZXFz9+vV9fX3/+OMPflcpkUi0NhKQkZHBp1kNtxBkmfhWL8LNofMNILU0U/HVV6xZs5fsjlhczKKiWHAw8/JiEsmWpzvI1q9f38fHRyqVqsoyKBQKXtbNyMho8+bNWui+TCbji9hu3bqlhcNVWVZWlpmZmVgsLqdctNCCgoKglSI/jDG2eDEDWPPmLD//JeMOiYksJIR5eTEzs44NGvC3lp2dnY+PT2hoqGonS8ZYcHAwv/BqbTOpsWPHAli+fLl2DqefeI3Xbt26aedwpXfAvXXr1tq1a0eMGKGq+QPA0tIy7pNP2Nq1rIInPh+kMDOrxCBFVf3vf/8DMH78eKEPVCZehLdFixYvfaXeBe5FRUVGRkZisVi9CF1eXl5gYCCfZbOwsAgMDNTggNbBgwf5raGlpSV/YzVv3vy7777T2LqZ5GTm5sbMzGb07Fnx0jGlf+RK5eZeunSJr5MD4OPjU5HlDqWFhYXx6iJ2dnZ//vlnFVp4kYSEBIlEYmRkpFpjIJyQkBA+0TF69OikpKSAgAC+HkUsFvv4+GhsJa5MxsaNU44KVK9A0OXLlwMDA9u3b6+60onFYltbW9UOF9rBOyBIlrbuJCUl8VvZf//9V+ON5+Tk1K1bF2orZ4SlUDDVzjLq9bblchYdzZYuZYMHM3PzZ8OiZma3JkwIDg6Oiop60ZAnD7DEYrEWNkw9ffo0gDZt2gh9oOrji+3WrFmjqw4kJibyxIPc3Fxhj7R8uTILqzKTz7Ls7J07d06ePFm1WQG/A+zdu/fChQs//PBD/nC9arWr8Pjgpbe3t9aOqIfy8/PNzMyMjIzUb6IEUlhYyEtR1alTR7X7sopMJuMjU25ubiKRqLhJE+VFqWVL5ufHwsJYOVv3JCSwli3Z7t3C/gCMMcbi4+P5uIb2K3QzxuRyOY9CX7plmN4F7suXL5dIJJaWlvXq1QsJCVGvgXXjxg1VHcM2bdpUv0BEfn6+qppKjx49Ll26FBoaqqrxYmlp6efnV93P4OPHmb09A1ibNqzyTSUlJfn6+vL+tG7dOqwCKa1FRUXBwcG8SE7z5s33799fpX4rpaSk8IJo/AZAU+c/r0Gu2rlWIAUFBXzcRSQSBQQEqOZwU1NTAwIC+E0RD9+vV68+IysqYm+/zQBmYcE0F15fvHhx3rx5fDWqcElcL8JLKOpkjY6gFixYAKBTp04avzRv2LABWlmz8ZzS9bZ79XoWrIvFzNWVzZzJDhxgFVujtnTpUn7KCD1ayUeRp0yZIuhRNGLr1q0ABmi0JGKlyGSyxo0b29nZubq6qu+SpmGrVjGRiIlErBplQPnyCW9vb755DQ+DTExMNDv081L37t0DYGVlpeUCx/pm4MCBAITOf8vNzX311VcB2NjYvHTHqwf37rGffmI+PoxvTMH/MzdnQ4awFSsYLxBSYucvLU788urY2i7i9BTfDvalYZt+Be688oNYLOa/OwBubm4l8rwPHjyo2gHLy8uryoH1pUuX+B6zJeqXy+XyPXv2qPZjE4vFw4cPP1W18bn165mJCQPYa69VcPlOmQ4fPsxTrwAMGjQo9sW1006ePMlvPEQikZ+fHy92WU18bSu/EWzcuLFGCuqdP3/e2Ni4YcOGoaGhAn0I3bt3z93dnd+AlVkdJSkpyd/fn3+6GBsb+/r6JiYmVuFAivx85u2t3KdO09uUsKcFUlVRVEZGRkREhNCrqBlj69evB/Dee+8JfSAtKygo4PlI1dy4qrQxY6517jxF80XQyle63va0aaxxY+bjw6RSVqUkt7Vr1/LhjAULFmi4t2p4SFGpFD5dyc7ONjc3F4vFQiQNvlRxcfE777wDQBUKu7m5/f777xqu7P7zz8pqRWvXaqS9nJycXbt28R0udbJCgEcR2llDpbf47fHUqVOFO0R6ejqvQtaoUaMLFy5U4jtlMmUiX58+TCxWRvAffSTozl8vxRcGfPbZZ1o+LjdjxgxUoIKCHgXuPGoXiUR8RjI8PLx58+Z4mpGsXuGkmvUTZTKZaky6Q4cOL1qqmJCQ4O/vzxdQbujbl3XtyqTSCo5aMZmMBQQo34h+fqzaY3vFxcVSqZSXv5VIJH5+fiUmU3gtdl6UplWrVocPH67mEUu4efMm3/SH3xLklDOxVQEJCQmqNbienp4aTy04f/48n7R1dHSMjo4u55W3b9/28/Pjq6ZMTEz8/Pwqsn5XJS8vb8hrr23t04fZ2gq0SZ1UKoVanaXffvtNO1PAfD+I5s2bC30g7du5cycAW1tbDW4Ke+mSslKZwPV+Sildb/tl06wVsXXrVn4xEWgLEj6JLxaLXzoprCdGjBghxM3eSxUWFvLlRtbW1gcPHgwJCeEV6wC0aNEiJCREMwWmeI1RgH37rQZaU1O6wMvXX3/t4eGhhRqF06ZNa9588MqVLxkArtn4cvkKbutTBSkpKXxUsXnz5tXa9iE1lW3dysaMYfPnC7fzV0XwLc909cG3ZcsWAG+99Vb5L9OXwD0kJEQkEqmidu7JkyfBwcF8rLdOnTol8rxL1E+s4NhtiQD0pVe9tLS0JUFBT1q0UEbh9vYsMJCVv0D28WPm5aVcpKjRggmPHz9W/cjqqURHjx7l5Q4lEom/v79AtQL5DQ8f9WnZsuV///1XndYUCkVoaGjDhg1Vd18amR9gjP3222/8dsvDw6OC2/LdvHnTz8+P/2J5+F6RhWjZ2dl8Iy0HB4e8au8H8SIxMTFQSwVOTEwEYGdnJ9DhVORyed++iwcMOJGSIngJCO17/fXXAfj5+WmqwU8+YQCbPl1T7VVY6XrbGsoN2L59O1/M/eWXX2p2WiwuLu69994Ti8Vt27bVYLOC+uWXX7SfB1VYWMhvGGxsbE4/HXrku6Sp1sDY2dkFBgZWbRWT0h9/KHf1EmabTV6jTLWayMXFBcIsMilh504ZoP3AT7+kpqba2to2aNCgbdu2wcHBmk12T0xM5Hu/dOzYsVIDXuUpc+cvLVIoFLwkUeVmDzSEJ9k3a9as/JfpReC+bt260lG7SnJysirPu1WrViXyvKOjo3nVegDdu3cvsXVtCaGhoaqUj9LrJ8pTWMjCwljPnsrw3cSE+fiUPYlTVMQ6dmQAa9SICTNJd+HCBT7/CKBLly7e3t48q8fZ2fn8+fNCHLHE0XmKkZGR0bFjx6rZWkZGhupWpGnTpqGhodVpjZdN5L8NPz+/ylZEiYuL8/X15UkCFhYW/v7+5ZQwysjI4PODjRs31nyBGjXFxcXm5uaqDdUUCgWfeBG6yj5jytvPGrQL0zPXrl0zNTUVi8VnNFEVOD+f1avHAKbprWwroHS9bc0JDw/nN+qTJ0+ufgm//Pz8sLAwVRZinTp1TExMNL77r0BycnJ4tozGApSXKSgoGDZsGJ8aKr1GXC6Xh4eH80uQ6npVhWvCzp073/fwKDI2ZosWaajjJfGiWKqdDT777DMAWigvk5HBjIyYiQkTekGv3lIoFG+88QYAvqCLv08mTpyokYLIsbGx/JbM3d29WveNJQg2ElFxU6ZMATB//nwtH5cxJpfLeQWe8meDdR+4r1y5kkft5W9ZcvjwYR4v8jzvy2oDnAqFYvPmzY0aNeLDpWXuKJmamsorA6Caiyz/+4+NGsWMjJQRvKcn47dl6mspNm1iXbsywWpFM8by8/PXr1/P7wvr1atnbGw8YsQIjafHvEhRUVFgYGD//v01lWQZHR3do0cP1R+3apVecnJy+OiURCIJDg6ucmdiY2N9fHx4bGFpaRkQEJCenl7iNenp6TyBvlmzZlrYbZffmqpKEg0ZMkQ7ycFz5zKA1aBdmJ7z5ZdfAujVq1f1h5O3bmUAc3fXSL8qqajohft8aUJERASvwjRmzJgqL+eNiory9/fnJaoA1K1b18/Pj1caaNOmjaam2oTGU1ZWrFihhWPl5eV5eXnxJJOLFy+W88pjx46pxm74fr3nzp2r4FH++ecfnjK6bdkyTfS6bHwfU9VqmV27dkFbK33d3BhQzSpfBmzJkiX8xu/GjRuRkZGq9wkANzc3qVRa5SJFR48e5bnKAwcO1PD5K+RIRAXt3bsXgIuLi/YPzRjjKSH79u0r5zU6DtylUtarV5KxsXlFKm2p53nz5Ar1io1ZWVlfffXV559/Xvob//zzT/5dNjY2W7du1UC/b95kAQHM1paJxSwhoYy1FALfIx4/fhxAz549165de+TIkY0bNwJ49dVXBT1oCZpdGiWXy0NDQ/k+jubm5pUtf5mYmMiXLNerV++gJqq9Xrx4UVXCyMrKKiAgQPVme/DgAb+HbNu2bXU3hKsYvmBl8eLF/OHcuXMBzJw5U+jjRkQwgPXtK/RxdCM7O5uPGFVznocx1r8/A5gWi92VUn697eo5cuQIHwR65513KlWjIzMzUyqV8tSI0uFCfn5+165dAUycOFGgnmvW9u3boZVdyXJzc/nKXXt7+3JKEai7dOmSr68vT20C0KdPn/Dw8PK/JTIykg/Efvrpp5ro9QvxnV8bN27MHz5+/JjvEVGpK3zVzJzJACbMMg19d+rUKWNjY5FIpL6d0PXr1wMCAuzs7NRvocu/Myxtz549/GZ+xIgRmt9pTuCRiIp1oYgXc1Ntc6FNn376KYCgoKByXqPLwF1Ve2rz5kpsbKGe512/fv0SJSNLyMrK4hv9ABg8eLCGZzmzs1l4eBlVHYQXHR0NwNXVlT88evQogL6GH2E9evTIz8+Pjwq0atWqgulM//33H8+Vb9u2rcbqsjPGGDt58iTf8Jy/2QIDAxMSEvgdQocOHbRWYmLbtm0A3nzzTf4wPDwcgJfXEKGPm5HBxGJmbq7NYlxatW7dOh4eVWfThtxc1rkzs7TUciqmVp09e5aPl3t7e1ck3oqKivLz8+NLTfiYn5+fX+ltLOPi4vhrfv31V2E6rkk5OTl16tSRSCTff/+9cGsrMzMz+QxbFXLwbt++HRAQwIdC+ZBhaGhomZ+Px48f51mjEydOFKq45FMKhYJfn1VVd/kNmxZKY+3fzwCmrQ2I9Mjjx4+bNWuGFywuLygoUE9aU91RV+Rd/csvv/D7wylTpgi4A66QIxEV8e677wJ46XaZmpKZmam6v+LFZ0eNGlXO63UWuEulyqh99eqqfHtMTEzfp7sAli4ZyZ04cYKvnDA3Nw8ODhbqTaaLtRR8DKNDhw784dmzZwF0795d6ONqx9GjR/kW3zxQKD9xUyqV8uvI66+/rrE9s553+PBhXl0VTyuyubi4aLMaRkJCAgBHR+U69wcPstu3z6tfnwm/cThr354BrMJz74ZEJpMNHDiQb5nUoUOHKVOmLFq0KDQ09ODBg1euXKnsCm9dDM1oVXR0NJ+3HDJkyIs+4O/fvx8SEqIq1ysWi728vMLCwspZbbJ27VoA1tbWOhncqqx9+/bxoi4NGjQICAjQ+K17RkZGz549ATg5OVU5By8rKyskJKRx48b8r9CyZcuQkBD1P9mpU6f4FMr48eMFjL3U8B2LNzwtDz9zZpC7++KQEMH3KcvLY6amzMiICb8BkR5RKBQ8Mbhnz56qKbKsrKzSWcRXrlwJCAhQ5bDZ2Nj4+fmVc7u4evVqvgZMoGJT+oNPr2mnjGl4eLijo6NEIuELfng45+TkVM636CZwX7++WlG7Snh4OL+t5CUjVXkLBQUFqtqI3bt31+wobEm6WEtx8+ZNqG2Ne+nSJQCvvPKK0MfVmuLi4pCQEP7pwgsKlf7sLy4u/uSTT/gVx9/fX+hPoGPHjvXr12/ChAlvvPFG6ax3QSkUij59HllYKFRjEI6ODKjCjl6VNn48A1iN24WJMcbmzJnDJ1KmTp2Ksri45L3yCnvjDTZpEps/n/38M9u3jwUFsaZNlffmixez0FAmlTLVQru5c6uzd42+i4uL48lF/fr1U09slcvlkZGRPj4+qlSNpk2bBgQElI7F09PTpVKpm5uberknXqTc3d1d//fKkclkv//+O1/fwoeEJk+enJCQoJHG09PT+XbXzZs3r9qeEuoKCgpCQ0N5tTEADRs2DAwMTE9Pj4mJsbW1BeDr66udqJ0xtmrVqnr12nz+uXIXuT//ZAAbOFAbh+ZpbDVyhf2L8Nrt9erVu6220O777783MjLy8vIqvYEXXzWuKtAsEon4/XaJ85FvqywSiZYJuSJCT2RnZ5uamhoZGQk6Qnf79m3VlH6fPn2uXr0aHh7u4OBgZGRU/p5ZOgjcVVG7Rkri5uXlBQYG8lw9CwuLwMDA6OhoZ2dnlNpZSSi6WEuRkpICoFGjRvzh9evXAbRu3VoLh9ame/fuqQoKtWvXTj15/dGjRzwN1NTUdPPmzVrrUmUr1WjK4MEMYDt3Kh++9RYDmBZ+bqmUAazG7cLEDh48aGRkJBaLIyMjP/nkE76B7tdff+3r6ztgwIC2bduam5tbWSlU+/qp/nv7bebkxPz9GWNs8WK2eTOTStm8ecpm58ypyYE7Y+zq1at8TXy3bt0eP36cnJwcHBzMN9wAwCODsLCwEstY5XL53r1733rrLb4UEs/vyJuens53XZgzZ47Wf6AqOnbsmI+PDx8b4ktCIyMjq9NgamqqauVMcnKypvopk8m2b9+uWmZgaWlpYWEBYPTo0Rrev6lcsbGFYjFzdFQ+fPRImYOn8QTp0hYuZACbNk3wA+mJ48fzmzVrLRaLSySafvnll6r76rZt23733XelS8HwDDf+DuEBRkBAwK1btxQKBa8FZGRktHHjRi3+NLrEi0AIFF3I5XKpVMqHJq2trUNCQm7evKnapd7Dw6P84WZtB+6ajdpVEhMT33zzTf4z83dnu3btSpfQEoQu1lJkZGTwvzd/mJyczEe5tHBo7Tt06BAvWiwSiXx9fVNTUy9dutSiRQsATZo00Ug5P/33zTcMYN98o3wYHKylTyO+r1AN24UpPznZqXFjAIsXL969e7dIJJJIJKVrmz56xC5eZLt3s3Xr2Ny5bNw45uXFpk9ny5YxFxcWFfUscHd3Z/7+zN+fde9ewwN3xlhiYiKP1OvXr8/nzfn1dtmyZaW3Tbh7925wcDA/W8uJ7I8dO9asWesuXU5qYmG59ty4cUO1SR8AV1fX0NDQKtTeefDgAU8ObN++vUArZ1TFZwYOHDhixAgtT24oFKxBAwYw1UTCK68wgFW7nvDLnTjBrK3ZjBmCH0gfpKWxpk2ZjU360qW/lP7q/fv31W+zTU1NfXx8IiMjSwzAZ2RkrFixgm/Bzs9Z/i1mZmY7VUNHtQDP4mvTps3evXs1u5A6JiaGz60B8Pb2vnXrliq/wMbGJiQk5KVTYVoN3AsLWZcuTCRiFSghUxWHDh2ytrZu3Ljxe++9J9AmRC+k3bUUBQUF/MTjDx89esSnxrTWAS3Lz89XzatYWVnxJe3u7u462YFcJ3buZADz8lI+PHRIS/UH5XLWuzebMqX6m//qjaIi1rt3VufOn7/77q1bt3gho++++66C3716NfvhB3bmDOvenS1YoAzcx49nx4+z48fZuHE1P3BnjKWkpDg5OTk6OpqYmJT52S+TyXjyDN+TmM8HBgYGlrNeZcECGcCaNtXIxq9alZaWFhwcrMop5xuaVrzQ3v3793mQ1KFDh4rs+1YdMTExxcXFOpk2HDmSAc82JJw2jQHsaaEsAcnlTItTC7okl7PXXmMA69mzvHRdVWKb6txs3759mXsz8QF4MzMze3t7S0vLA7WsrOaDBw94ATeeFOfl5RUSElLN2bC8vDxVFnfTpk3/+uuvEkF8BQuoaHvEPS2N/VLGraDG8GruWtibRuf4WBe/M8vLy6tnbt7haeZMTXXjxo033njDzMysefPm77zzjrbvzXQqJYUBzNqa8QApO5uJxczERBtzzTXN558zgDk4sHv3Vr/zDoA333yz4oU1eODOGPvkE+boWOtSZVTGjx8PYHWpVUoJCQkBAQH29vblj+qVJpczT08GsDfeYAKXOREEzynv0KED/8Hr1q3r7+//0o/5O3futG7dGoCzs7M217trX0gIA9i4ccqHYWEMYIMHC37c1atLLkqpqRYvZgCzs6voECKfDXN0dOTvWDMzM19f3+jo6BIvCwsLAzBo0CDN91jvJScnBwYGurm5qcrviEQiV1fXefPmnT17trJLRHbv3s3XZPId7lNTU1VBfIsWLSq1JagggXvpU0U767d4AkmdOnW0tuZGhwrbtJHZ2ytz6+VyBjCx2CA/8SopISGh/E3FaqomTRjAVAvhpk1j//sfy8nRWPsKBSsdZtS0j73du5lIxIyN2fHjbPp0Jhb/O3JkpZYar17Nli9njLHsbGXgvn79s8C9Zi9OVcfL7EqlUtUzW7ZsUVX6AtClS5eVK1dW6nebnMzq12eAhhMptYlvaMr3TgJgYmLi6+v7olrst2/fbtmyJc+x0eTek3opJua5pLvUVCYSsWbNNPmR9fffrHRZ0dWrSy5KqZGOHmUSCROL2d69lfvGoqKisLCwgQMHqmJTd3f3TZs2qV6we/duAEOHDtVshw1LWlpaaGioj48Pr0LGNWjQwMfHJzQ0tCLl7HgtBADdu3ePiYnZvXs3X9jDg/icSn6QCxW462T91smTJ6FW3byG43usq6a3TEwYQAOwNdiIEQxg27YJ0nhiIhs4kDVvXvJOoEZ97OXlsYYNGcB++IH9/jsDmKlpzSx1KbyPP/4YgPrGeXyxuJWVla+vb5VXav71l/LPEhOjmX7qSlRUlK+vryobge+IpD7tcPPmTZ46zJf56rCr2iGXK+/KVJVOrl7VWNSem8v8/BjArKxYiT3xVq8uuSil5klNZU2bMoDNnl31Rq5duxYQEMBLvg5Uq/jz119/ARg5cqQGOmr4iouLjx07FhAQwNfdcRKJpE+fPsHBwfEvLvR28eJFa2vr4ODg5ORkVckNFxeXiu9zrE4MYUyfjmPHEB397Jl9+zBjBmbMwP79Ah0TV65cAaCarKzhzM0BID//uYcFBTrrDxEYT4Q7d07Dzcpk+PZbvPIKDh9Gfj6uXy/5gtLnsqGqUwfh4Zg+HUOHYtIkAPjhB3TrputuGSQekhYXF6ue+frrr3/66af79+9v2bJFNeRcWSNHYsoUFBZizBg8eaKZruqEm5vbli1brl275u/vb2FhceLEieHDh7u6um7ZsqW4uDghIaFv3763b9/28PA4fPiwqpB2DSYWg9cbPHpU+Uy7dng6yFst587B1RXr18PMDLNno2lT5fNyOR48AACJBOvWYepUyOUaOJy+USjg64t799CvH+bPr3o7bdq04WFlaGjo7NmzVc/zc1xVkaaWk0gkHh4ewcHBV65cSUxMDAkJ4ftYnThxYtasWR07dmzVqtWMGTMOHjxYVFSk/o1dunRJSkqytbXt1KnT1q1bLSwsgoODz507161KH0BCBe6lT5VOnTB6NEaPxtPFyppXuwJ3MzNALXAv8ZDUODxwP3tWk21evow+fTBzJvLz4eODy5eh2p/+xg38+itQwz72evTAqlVYtw7Z2Xj3XXz8sa47ZKj4B7lMJlM94+XlNWHCBFUtuSpbvhxduiA+Hp99Vs2WdK9FixYrVqy4c+fOokWL7O3tL1y4MG7cuObNm/fs2fPevXsDBw7ct28fryZRG/TvD6gF7tUnk8kXL0afPrh2DV27IioKAQHghY6SkjBoELy8wG8t3d3Rsyc2btTYofXHqlU4cAANG2L7djyd4Kk6MzOzsWPH8tkzjgL3F2nZsuWMGTMiIyMfPHgQFhbm5+dnb29/8+bNlStXDh48uF69esOGDVu/fv2DBw8AXLp06bXXXps8eXJ2dra3t3d8fLwqwb0KhArcUepUcXREnz7o0wdP10JoXu0K3MsccafAvebq0QM//YQNGzTTWlFRUWBg4Pvv5549Cycn7N2LsDA0aAAACgXWr4ezMyZMwP37QM342MvJQXKy8t/ff48ff8T69TrtkGHjI+7qgbummJnhl19gbo6DB5GZqfHmdaB+/fpz5sxJSkoKDQ3t1KlTSkpKs2bNXnvttT179lT/PseA9O8PIyNkZGimtTt37gwaNHD//usyGfz9ceYMnu7Yi61b0aULjh5FZiYeP1aO6wcFaea4+mb8eLz7LrZuRZMmgrRPgftL1atXz8fHRyqV3r1798iRI1999VWHDh3y8vL27NkzefJkBweH9u3bu7m5nT592sHBYdeuXaoE9yoTKnAvcarwh+fOQaHQzOxYma5evYraFrircmMoVaam27oV8+aBn+9BQdiypepNnTlzxtXVdeHChUVFE6ZPx+XLGDJE+aXLl9GrFyZPRl4ehg+HlZXhf+zJZPjoI7RqhbfegqMjzpyBSISpU1FrRjqFIFzgDuCVVxAejpgY2NgI0bxumJiYjB07NjY29sCBA8ePH//77795Wdvaw9kZ6en46y8NNLV58+bOnTv/999/jx6NOXRIsWIFTE0BICsLvr4YOxZZWXjrLVy8iEWL8OmnAGBlhaQkjBungaPrxJo1cHBATg7w9Pq/fj3mz0fduti+HcePCzWwwlM+KHCvCIlE0r9//2XLlsXHx9+6dUsqlXp7e0skkoKCgrp16/r5+cXFxal2HKrWgarfRGnTpin/wU8Vzt8f7u7YvBkLFwpxTBQUFNy+fdvY2LhVq1aCHEDfUKpM7WNkhDlzsGIFADBWlRby8/MXLFjw3XffyeXyVq1arV071dNT+aXiYvzwAwIDUViIJk2wZg1GjHj2jernsoHZtw+3biEpCWZm+OsvzJqFf//VdZ8MnnCB+5o1WLIEV67AygpBQXB0REEBUlKU+bv83nXiRI0fVktEItHgwYN13QvdWLeu5F927NhKN5KVlfXxxx//+uuvAN5++22pVFqvnnL88fBhjBuHu3dhZYXvvoOfn2a7r3ulr/+qT4GqfRxUBI24V03z5s39/Pz8/Pxyc3Nv3brl6Ohoo7lxCAFTZUpwdweA2bOFWnJ040ZO9+7BffqMU+2qXcPR4tTaR32dqFyO7OzKffu1a9deeeWVpUuXikSiWbNmXb582fNp2B4Tgx49MGsWiorg54crV56L2g1bRATee095Z/vmm7h2TTlsRaqBZ2cKNOLOAxROawEK0YLSf9lKefToUadOnX799Vdra+utW7fu2LGDr+stLGRffonBg3H3Lvr0wcWLNTBqh45qfvDAvbaEVQKwtLTs3LmzBqN2aDNwHzMG3brh3j0sXy5I+3FxdqdPf2ljo6EUYP33+efYsAGdOysfrlmDqKhnD0lNpL5ONC4O9erBwwNLlyIhoULf3qxZMzMzsy5dupw8eXLJkiV8J9onT54EBm7o0QMxMWjbFkePQiqFWrHaGkE9K0YiUQbxpBokkrkAk8kWCdG4TgIUogXqf9niYvj746+/KpH13qBBAy8vr549e0ZHR3/wwQf8yfj4eA+P3vv354rFCAzE0aNo0UKY3uuaTmp+iMWTHR2LrKwMN1eyBtJe4C4S4bvvACA4WFmkSbOuXAGAWpLfDgCvvYaJE9GkiXLVnbMz3NwobbfGU60TTUkBgBMnMGsW2rdH166YMwdnz0KheOH3mpqaRkREREVFqfZYPnbsmIuLy8KFfh4eUTNn4sIFqG2hU1OMHIlt25Qzfdu3o2dP0LRvtfFfoTAD7roJUIgWqP9lr1/HqlV46y00aIBu3TBjBnbsePks4tq1a48fP84TYhljK1ascHNzi4o63aDBZydPYv58VLVQh2HQfs2P3Fzz5GRjubyOUAcglae9wB1A//4YPhy5uYKkude6wB1lrbojNZr6OtEhQ5CWhm3bMHo06tbFpUsICkKPHmjaFJ9/fj88PPxJWUlpzZo149mKWVlZM2bMGDBgwLVr1zp37rx0KVu6FDVzsdyAAbCyQtOmcHZGUBAWLNB1h2oCXnhOoMAdughQiHao/rJ2dli4EAMGwNgY0dFYuRKjR6NhQwwciMWLcfp0dJmJWObm5jxNKzU1ddiwYZ9++mlBQYGvr+/u3cufDkfUWGXW/AAgkwlY84PX06SxDr0iyOLUcixbhr17sX49pk17VrxJI2pj4E6r7mqTMtd8jxmDMWMgl+PUKezZg127kJCA8+fDli//1MzMzMPDw9vb++23326q2pUEABARETFlypTk5GRjY+Ovvvpq4cKFNTmF0dgY27cjKwu5uXj+90CqTNDAXSTC/v2YPx+7dikfqn+JGC5V6LlrF+zsMG4c5s5Ffj5OnMDBgzh+HGfO4N9/ER+fMnduNwsLi169enl5eXl5ebm6uorU/vY7d+708/N79OiRnZ3dhg0bNFKpQ8+Vef0HkJyMd97B+PFCrdimwF0fVX3v16r6+GMGMG9vTbYpkzEzMyYSsexsTTar76ZOZevXK/8tk7EmTWrZz0/KcPEiW7Fie48ePcRi5XyaSCTq3r37okWLYmJi0tPT/Z6u2+rVq1dcXJyu+0sM0tq1DGBTpgjS+P79TCxmffsymUyQ9oneevyY/fEHmzfvdLt27dQDFXt7+/fff/+nn36Kj4/39/fnT7766qv37t3TdZd1bPt2BjAzM3bhgiDtf/UVA9jSpYI0TqpGB4F7aiqrW5c5ONw6ckRjb7Tr1xnAHB011Z6BmDqVbd+u/LdMxpycWFGRTjtE9Mj9+/f5WFSdOs/SE+vXrw/AwsIiJCRELpfruo/EUG3YwAA2caLmW05OZnZ2DGCLFmm+cWJAkpOTN2/e7Ovr20RtbyE+HmFubr5mzRqFQqHrPuqFyZMZwFq3ZpmZmm/8008ZwJYv13zLpMq0muPONWyI//3v8MOH7T//fIKinJV0lVEb82RAq+5IeRo1ajRx4sRdu3Y9fvw4MjLS39+/QYMGAJo1axYbGztjxgzVkDwhlSVQqoxMhvfew8OHGDQIX3+t4caJYXFwcBg3btyWLVvu3bsXHx+/atWqN99809jYWCwWnzp16uOPPxZR4hQAYMUKuLrixg1MmqT5xilVRg9pO8ed++ij3t9+2+j8+fPbtm0bq7YHw82bN2/fvp2Xl5eXl5ednZ2dnZ2fX5yb+3VmJnJzkZeHvDxkZGDo0JILzGpp4D5gAH7+GU2bolkzFBZi505dd4joIzMzM6+nhg8f3q1btxY1tV4a0RaBAvdZs3D8OBwc8NtvNbw8CKmUDh06dOjQYfr06R06dOBbpBMVU1OEhcHNDTt2YM2aZ9nwGkGBux7STeBuZma2aNGisWPHfvPNN2+//bZqKj8wMHDbtm3qrxSJxGLx16qiYFybNiUbrKWBO626I5XRunVrAImJibruCDF4Lwrcb99GSgpcXatSK3/PHvzwAyQS/PYbGjTQQCdJzdO5c+erV6/GxsZ27dpV133RI61aYcsWjBiBzz6Dmxt69qx0C0VFiI1F3bol4ysK3PWQbgJ3AB988MHKlSujoqKWL18+e/Zs/mSHDh08PT3r1KljYWFhY2NjaWlpYWFhZSW3tjaqUwcWFrCygrU1GjYs2dqCBRg1CrX0RLa2hrW1rjtBDECLFi3EYvHt27eLi4tpC2tSHc2bY+JEdOtW8vktWxAYCIkEbdvCwwN9+sDNDR07vrwazJ07GDcOjGHZMvTpI1CvicHr3Lnzjh07YmNjdd0RvTN8OKZPx6pVeP99nD8vt7F5yYxVbm7+hQvmMTHg/8XFobgYn35acotMBwe4ucHeXsCek8oSMd1tIX306NEBAwZYWlpev369UaNGVWtkzRosWYIrV2BlhaAgODpCLfWGEPIcR0fHu3fvJiYmtmzZUtd9ITXQ6tXYsAFxcc/2TrKyYhKJnZubS58+fdzc3Hr37s1XSKsrKCiYOPHI778PGTYMf/5JNR/JC+3cuXPUqFGvv/56RESErvuid4qLMXAgMzI6WrduyN9/7yyxBiArKys2Njb6qcTEJJEop6BA+VUjI7RrhzFj8M03Oug5qRzdro0dNmwYgKlTp1bt2xUKtno1c3Ji/v6MMbZ4Mdu8WZPdI6SG6d+/P4ADBw7ouiPEgK1ezZo2VdaeXbyYhYaWfEFODvv3X7ZkCXvzTTZwYKr6J45YLO7YseO///6r/vopU6YAGDJkoRBlMUhNcv36dQAODg667oieunXrjq2tLYAffvhB/flFixaViP1MTU2HDcuYNIn9+CM7dYrl5ZVs6qWnOdEVnaXKcN9+++2+ffvWr1/v7e1tb2+fnZ395MmTvLy8zMzM3NxcvkpVoQi4f982Lw85OcjOVi5R5Xndr72GIUMwfTq2b0d0tG5/FEIMQKtWrY4ePZqYmDh48GBd94UYMCMjzJmDFSsAQC5HXh4sLJ591dISAwZgwAD+qGFKyr3o6OgTJ04cP378/Pnz8fHx9erVU734999/X7dunampaVDQUEr6I+Vr2bKlhYXF3bt309PT1d9FhGve3GnLli3Dhw8PCAhwd3fv8zTtrE2bNpaWlu3atevYsaObm5ubm1u3bt3Myl2Jwthzp7nukjNISToO3Nu1azdp0qTNmzcPHTr0Ra9xdp514ULZX8rLAwCJBOvWYepUeHsL0klCaoxWrVqB1qeSalMfLrl5E9bWaNcObm5wc4OHB1xcoF5rtEmTJk2aNOHzq8XFxTExMZ2e7pt97do1viPYypUrXV1dtf+DEMPCZ2zOnTsXFxfXt29fXXdHH3l7e3/++efff//9e++9d/78eV4F+K233vLx8Sm/BHBKCqKjlf+dPQtfXxoV1VM6DtwBLFy4cPTo0R999JFqNaqFhYWtrS1fomplZdWoUZFIBAsLWFvDygp8laqtLSwsYGKCNWsAwN0dPXti40aUmg4ihDxDgTvRCPXhkocPIRYjPh7x8di6FQBsbJTX5B490KPHY/WkdmNjY3d3d/7v/Pz80aNHZ2dnv/vuu6oNfQkpX+fOnc+dOxcbG0uB+4sEBwefPn36xIkT48aN27Nnj0gkkkhKBntFRUXx8bejo9vyxakXLypHQlVSUtC0KY2K6iPdB+7169fv37//jRs3qtwCX4ARFIRduzTVKUJqprZtew4YcFwiafLylxJSLvXhklWrcOkSjh9XDtfFx+PAARw4AEfH4uTkBo0bN/bw8OCLU7t3725qaspbmDp16sWLF9u2bbt+/Xrd/izEgPTuPffs2RXXrpnquiP6SyKR/Pbbby4uLhEREcuWLQsICABQVFR0/fr1aDUiUZ3i4seqoq6NGysnzdzc0KkT9u5FURGNiuojXVaVIYRoWWYmbG1haYnsbKrdQapozRpl5bicHHTqhEWLMG7ccy9IScGZMzh9GmlpiTt2dM1TG8ozNzd3dXXt0aNHXl6eVCq1sLA4c+aMKnOGkJc6dAheXvDwwLFjuu6Kftu7d6+3t7dIJPL09Lx79+61a9fU96o3MjJq27btgAFnW7a0dHGBqytsbZ/79pee5kRXKHAnpHapXx/p6XjwgErzEm2Qy+VXr17lI3wnTpyIiYnh0YOZmZlMJtu0adMHH3yg6z4SQ/LwIRo2hLU1MjJo9OElvvnmm3///ff06dMAJBJJ27Zt3Z5ycXGxUF9RTgwHBe6E1C7u7jh3DidOoHdvXXeF1D6ZmZlnzpw5c+ZMYWHhuHHj2rZtq+seEcNjb4+0NCQlwdFR113Rb3K5PDc3948//nB1de3UqZOJiYmue0Q0gAJ3QmqX997Db78hNJS2KiOEGCQvLxw6hH/+wRtv6LorhGhdebWBCCE1T6tWAEB1ZQghBqpzZwCIjdV1PwjRBQrcCaldKHAnhBi0V14BKHAntRUF7oTULhS4E0IMGo24k9qMctwJqV34thp2dkhL03VXCCGk8vLyULcuJBLk5IDWW5LaRvcbMBFCtKlxY3h4wNERxcUwNtZ1bwghpJIsLLBkCRwdQQOPpBaiEXdCCCGEEEIMAOW4E1LDrVkDBwfk5ABAUBC2bMH69Zg/X/nVefOwcaPuOkcIIS9DFzFCVChwJ6TmMzLCnDnKf/M5NtVMG025EUL0H13ECOEox52Qmm/6dGzfjujoZ8/s24fMTAA4dQp+fjrqFiGEVAxdxAjhKHAnpOaTSLBuHaZOhbe38plOnTB6NABkZemwX4QQUiF0ESOEo1QZQmoFd3f07PksE9TREX36oE8fODoqnzlwAAkJuuodIYS8BF3ECAEF7oTUBiIRAAQFPfdQ/d8zZyIjAzNnUnF3Qog+oosYIRyVgySE4OZN2Nrigw/w+++wtNR1bwghpJLoIkZqCRpxJ4QgNRX+/li1ij7wCCEGiS5ipJagEXdCCAYPRosWMDbGrFnPEkYJIcRQ0EWM1BIUuBNCCCGEEGIAKFWGEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogBoMCdEEIIIYQQA0CBOyGEEEIIIQaAAndCCCGEEEIMAAXuhBBCCCGEGAAK3AkhhBBCCDEAFLgTQgghhBBiAChwJ4QQQgghxABQ4E4IIYQQQogB+D9Fs5XvxARY0AAAAABJRU5ErkJggg==\n",
"text/plain": [
"<PIL.PngImagePlugin.PngImageFile image mode=RGB size=1000x200 at 0x11D1A5128>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"Draw.MolsToGridImage(mols[:5], molsPerRow=5)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"morganfps = [AllChem.GetMorganFingerprintAsBitVect(m,2) for m in mols]"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
"distmat = GetTanimotoDistMat(morganfps)\n",
"simmat = GetTanimotoSimMat(morganfps)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[0.53703704 0.51851852 0.43396226 ... 0.88659794 0.92079208 0.86538462] 1081\n",
"[0.46296296 0.48148148 0.56603774 ... 0.11340206 0.07920792 0.13461538] 1081\n",
"47\n"
]
}
],
"source": [
"print(distmat, len(distmat))\n",
"print(simmat, len(simmat))\n",
"print(len(mols))"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [],
"source": [
"def tri2mat(tri_arr):\n",
" n = len(tri_arr)\n",
" m = int((np.sqrt(1 + 4 * 2 * n) + 1) / 2)\n",
" arr = np.ones([m, m])\n",
" for i in range(m):\n",
" for j in range(i):\n",
" arr[i][j] = tri_arr[i + j - 1]\n",
" arr[j][i] = tri_arr[i + j - 1]\n",
" return arr"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [],
"source": [
"distarr = tri2mat(distmat)\n",
"simmat = tri2mat(simmat)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x11e676278>"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAWQAAAD8CAYAAABAWd66AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFvBJREFUeJzt3XuUXWV9//H3x0mCwaRcBBUIctEgIrq4xGi1RSyFhtaF1V4Ea0WXmnYVtLbLWqyt/IRfW/WnWNcqbU01rros0Iq2ppSKgEXbKphwU8KlhGhljII2CFKsZGY+vz/ODh6GmXPOzJy95zmbz8u1V/bZt+9zTPjOM89+9v7KNhERsfiesNgNiIiIjiTkiIhCJCFHRBQiCTkiohBJyBERhUhCjogoRBJyREQhkpAjIgqRhBwRUYgldQfY9b3tjTwK+L7j/6iJMAA8qOaebrxh8r7GYh05tldjsS6678ZG4hywfN9G4gC8ao/DG4t11eS9jcVaptrTxCOuuPtftNBrzCXnLN3v8AXHG6b0kCMiCtHcj76IiCZMTS52C+YtCTki2mVyYrFbMG9JyBHRKvbUYjdh3pKQI6JdppKQIyLKkB5yREQhclMvIqIQbe4hSzoSeDlwEGBgB7DJ9m01ty0iYs48wrMsej4YIun3gUsAAV8BNlfrF0s6p/7mRUTM0dTU4Eth+vWQ3wA8x/au7o2SLgC2Au+Z6SRJ64H1AH/xgf/LG197xhCaGhExgBYPWUwBBwL/NW37AdW+GdneAGyA5t5lEREBtPqm3luBqyXdCdxdbXs68Ezg7DobFhExL23tIdv+rKQjgLV0buoJGAc22x7dH0MR0V4jfFOv7ywLd55DvLaBtkRELFyBN+sGlXnIEdEqo/zLexJyRLRLW8eQIyJGToYsIiIKkR5yREQhJnf1P6ZQScgR0S4ZsphdU9Wg3379+Y3EAXjZsWc1FquNlaChuWrQqQS9cFff89XGYg1FhiwiIgqRHnJERCGSkCMiyuDc1IuIKETGkCMiCjHCQxY9K4ZERIwcTw2+9CFpnaQ7JG2bqUqSpEMkXS3pq5KukbSqa9+Zku6sljMHaXoSckS0y5BKOEkaAy4ETgWOAs6QdNS0w94PfNz284DzgD+tzt0XOBd4AZ3XF58raZ9+TU9Cjoh2GV4PeS2wzfZ22w/TqS/68mnHHAVcXa3/a9f+nwOutL3T9n3AlcC6fgGTkCOiXSYmBl96O4gfV0qCTnGOg6YdczPwS9X6K4CVkp484LmPMe+ELOn18z03IqI2c+ghS1ovaUvXsr7rSprp6tM+vw14iaQbgZcA3wImBjz3MRbSQ373bDu6v+TmB7ctIERExBzNYQzZ9gbba7qWDV1XGgcO7vq8CtjRHcr2DtuvtH0s8M5q2/2DnDuTntPeJM32ELuAp852XnfV6T8+5NdSdToimjO8ecibgdWSDqPT8z0deHX3AZL2A3ZWpe7eAWysdl0B/EnXjbxTqv099ZuH/FQ6g9P3Tdsu4Ev9Lh4R0bghzUO2PSHpbDrJdQzYaHurpPOALbY3AScCfyrJwBeBs6pzd0o6n05SBzjP9s5+Mfsl5MuAFbZvmr5D0jWDfa2IiAYN8Uk925cDl0/b9q6u9UuBS2c5dyM/7jEPpGdCtv2GHvtePdu+iIhF03/2RLHy6HREtItH97ZVEnJEtMsIv8siCTki2iUJOSKiEHn9ZkREISYnF7sF81Z7Qn5QzQywv2vNH3LD5PTp0vW47MYLG4kDcOAzTm0sVlOFR6G54qMpPLpwD3xg+vt0Cpchi8XXVDKOiMIlIUdEFCJjyBERZfBU5iFHRJQhQxYREYXILIuIiEKkhxwRUYgk5IiIQozwy4X6lnCSdKSkkyStmLa9bwXViIjGzaGEU2l6JmRJbwE+A7wZuEVS9yM7f1JnwyIi5mXKgy+F6Tdk8SbgeNsPSjoUuFTSobY/xMxVVSMiFtcIz7LoN2QxZvtBANvfoFM/6lRJF9AjIXdXnb7pB6k6HRHN8dTUwEtp+iXk70g6ZveHKjm/DNgPeO5sJ3WX1j5m5TOH09KIiEG0eMjitcCjClTZngBeK+nDtbUqImK+2vouC9vjPfb9x/CbExGxQAX2fAeVecgR0S4To3tTLwk5ItqlrUMWEREjJ0MWERFlKHE626CSkCOiXdJDjogoRBLy7JoqPnrk2F6NxIFmK0HvuOtfGov1vuP/qLFYTVWDTiXohdMhzVQIH5oRfnQ6PeSIaJXU1IuIKEUSckREITLLIiKiEOkhR0QUIgk5IqIMnsyQRUREGdJDjogoQ6unvUlaC9j2ZklHAeuA221fXnvrIiLmqq0JWdK5wKnAEklXAi8ArgHOkXSs7T+uv4kREXMwxCFkSeuADwFjwEdsv2fa/g8CL60+7gk8xfbe1b5J4GvVvm/aPq1fvH495F8GjgH2AL4DrLL9gKT/B1wHzJiQJa0H1gMctfdzWLXi4H7tiIgYCk8MJyNLGgMuBE4GxoHNkjbZvvWRWPbvdB3/ZuDYrkv80PYxzEG/IqcTtidtPwTcZfuBqhE/pMfPoe4ip0nGEdGoqTksva0Fttnebvth4BKg10tEzgAuXkjT+yXkhyXtWa0fv3ujpL0Y6i8GERHD4SkPvPRxEHB31+fxattjSDoEOAz4fNfmJ0raIulaSb84SNv7DVmcYPtHAPaj6qIsBc4cJEBERKPm0FXsHl6tbLC9YffuGU6ZLYufDlxqu/tVc0+3vUPS4cDnJX3N9l292tOv6vSPZtn+PeB7vc6NiFgMc5n2ViXfDbPsHge6x1xXATtmOfZ04Kxp195R/bld0jV0xpd7JuR+QxYREaNleGPIm4HVkg6TtIxO0t00/SBJzwL2Ab7ctW0fSXtU6/sBLwZunX7udHkwJCJaxRNDuo49Iels4Ao609422t4q6Txgi+3dyfkM4BLb3V3zZwMfljRFp+P7nu7ZGbNJQo6IVvEQpxtUD8BdPm3bu6Z9/j8znPcl4LlzjZeEHBHtMsLzv5KQI6JVhtlDbloSckS0ShJyD01Vg77ovhsbiQNwwPJ9G4vVZCXot19/fmOxvnjsWf0PGoJUgn788eRM04dHQ3rIEdEq6SFHRBTCU+khR0QUIT3kiIhC2OkhR0QUIT3kiIhCTGWWRUREGXJTLyKiEKOckOf8+k1JH6+jIRERw2APvpSmX9Xp6e/+FPBSSXsDDFJFNSKiSaPcQ+43ZLGKzkuVP0KndImANcAHep3UXRblZ/Zdw9Ern7HwlkZEDGCUp731G7JYA1wPvBO43/Y1dEpbf8H2F2Y7qbvqdJJxRDRpclIDL6XpV1NvCvigpE9Wf97T75yIiMU0yj3kgZKr7XHgVyT9AvBAvU2KiJi/No8hP4rtfwb+uaa2REQsWImzJwaV4YeIaJXHTQ85IqJ0k1NzfryiGEnIEdEqGbKIiCjEVNtnWUREjIrWT3uLiBgVGbLooalq0E1Wgn7VHs1V/L1q8t7GYjVVCRrgshsvbCTOrovf30gcaHEl6MmJxW7BnGTIIiKiEJllERFRiBEesUhCjoh2yZBFREQhMssiIqIQI1x0Ogk5ItrFpIccEVGEiQxZRESU4XHTQ5b0U8Ba4Bbbn6unSRER8zfKY8g9Z1BL+krX+puAPwdWAudKOqfmtkVEzJnRwEtp+j3SsrRrfT1wsu13A6cAvzbbSZLWS9oiacv/Pvz9ITQzImIwU3NYStNvyOIJkvahk7hl+7sAtv9H0qwPuNveAGwA2H+vZ43ygzMRMWImC+z5DqpfD3kv4HpgC7CvpKcBSFoBI/ytI6K1pjT40o+kdZLukLRttmFaSb8q6VZJWyVd1LX9TEl3VsuZg7S9Zw/Z9qGz7JoCXjFIgIiIJk0Nqa8oaQy4EDgZGAc2S9pk+9auY1YD7wBebPs+SU+ptu8LnAusofN6jeurc+/rFXNer0Wy/ZDtr8/n3IiIOnkOSx9rgW22t9t+GLgEePm0Y94EXLg70dre/b7cnwOutL2z2nclsK5fwNF9T11ExAzmclOvewJCtazvutRBwN1dn8erbd2OAI6Q9B+SrpW0bg7nPkYeDImIVpnS4EMW3RMQZjDThaZ3rJcAq4ETgVXAv0k6esBzHyM95Iholck5LH2MAwd3fV4F7JjhmM/Y3lUN495BJ0EPcu5jJCFHRKsMcZbFZmC1pMMkLQNOBzZNO+YfgZcCSNqPzhDGduAK4BRJ+1RTh0+ptvWUIYuIaJVhzbKwPSHpbDqJdAzYaHurpPOALbY38ePEeyudTvfv2f5vAEnn00nqAOfZ3tkvplxzidbnPe0nG3swpKnio00WHl2m5n5mXn3PVxuL9cAHpt+srsfSM97WSByAiS99urFYjRYeHWvu3+Dy09624Gz6iQNfM3DOec2OTxT1PEVreshNVoKOiHIN8sBHqVqTkCMioMx3VAwqCTkiWmUyPeSIiDKkhxwRUYgk5IiIQoxwSb0k5Ihol/SQIyIKMcAj0cXqV1PvBZJ+olpfLundkv5J0nsl7dVMEyMiBjfMF9Q3rd+7LDYCD1XrH6JTQeS91baP1diuiIh5aXVNPdu7n9FcY/u4av3fJd1UY7siIualxEQ7qH495Fskvb5av1nSGgBJRwC7Zjup+6XPOx+6Z0hNjYjob4gVQxrXLyG/EXiJpLuAo4AvS9oO/HW1b0a2N9heY3vNvns+dXitjYjoY5THkPsVOb0feJ2klcDh1fHjttPtjYgijfIsi4Gmvdn+AXBzzW2JiFiwqSIHIwaTecgR0SqjfFMvCTkiWmV0+8dJyBHRMukhR0QUYkKj20dOQo6IVhnddJyEHBEtkyGLHlIJemHaWAkaQIc08++iyUrQS170ysZiNfm9fNvWxmJx2sIvkWlvERGFGN10nIQcES2TIYuIiEJMjnAfOQk5IlolPeSIiEI4PeSIiDKkhxwRUYhMe4uIKMTopuP+VaffIungphoTEbFQE3jgpTT9SjidD1wn6d8k/Zak/ZtoVETEfHkO/ytNv4S8HVhFJzEfD9wq6bOSzqzKOs2ou8jp5ge3DbG5ERG9Tc1hKU2/hGzbU7Y/Z/sNwIHAXwDr6CTr2U56pMjp81c8c4jNjYjobZR7yP1u6j2qLqvtXcAmYJOk5bW1KiJinkrs+Q6qX0J+1Ww7bP9wyG2JiFiwSZfX8x1Uz4Rs+z+bakhExDBkHnJERCFKHBseVL+behERI2WYsywkrZN0h6Rtks7pcdwvS7KkNdXnQyX9UNJN1fJXg7Q9PeSIaJVhDVlIGgMuBE4GxoHNkjbZvnXacSuBtwDXTbvEXbaPmUvM9JAjolWGOO1tLbDN9nbbDwOXADPVOTsfeB/wvwttexJyRLTKpD3w0v0QW7Ws77rUQcDdXZ/Hq22PkHQscLDty2ZoymGSbpT0BUk/PUjbM2QREa0ylyEL2xuADbPs1gzbHrm4pCcAHwReN8Nx3waebvu/JR0P/KOk59h+oFd7ak/ITVWDTiXohWuqEnSjJicaC9XaCtdLlzUWaxiG+GDIOND9crVVwI6uzyuBo4FrJAE8jc5Dc6fZ3gL8CMD29ZLuAo4AtvQKmCGLiGiVIY4hbwZWSzpM0jLgdDpPKnfi2Pfb3s/2obYPBa4FTrO9RdL+1U1BJB0OrKbH6yZ2y5BFRLTKsGZZ2J6QdDZwBTAGbLS9VdJ5wBbbm3qcfgJwnqQJYBL4Tds7+8VMQo6IVvEQH522fTlw+bRt75rl2BO71j8FfGqu8ZKQI6JVJkf4Sb0k5IholbzLIiKiEMMcsmhaEnJEtEp6yBERhRjlt731TMhdc+922L5K0quBFwG3ARuqCiIREcVo7QvqgY9Vx+wp6UxgBfBp4CQ6L944s97mRUTMTZuHLJ5r+3mSlgDfAg60PSnpE8DNs51UvaBjPcARex/JgU9aNbQGR0T0MsoJud+j00+ohi1WAnsCe1Xb9wCWznZSd9XpJOOIaJLtgZfS9OshfxS4nc5jg+8EPilpO/BCOu8GjYgoyij3kPsVOf2gpL+r1ndI+jjws8Bf2/5KEw2MiJiL1s6ygE4i7lr/PnBprS2KiFiASQ/xBZwNyzzkiGiVEseGB5WEHBGt0tox5IiIUdPqMeSIiFEylSGLiIgypIccEVGIzLLooalq0KkEPWKaqgY91lyfw7dtbSxWk5Wglzz/ZY3FGoYMWUREFCJDFhERhUgPOSKiEOkhR0QUYtKTi92EeUtCjohWyaPTERGFyKPTERGFSA85IqIQrZ5lIekZwCuAg4EJ4E7gYtv319y2iIg5G+VZFj1r6kl6C/BXwBOB5wPL6STmL0s6sfbWRUTM0aSnBl5K06+H/CbgmKrS9AXA5bZPlPRh4DPAsTOd1F11+qi9n8OqFQcPs80REbMa5THkflWn4cdJew861aex/U0GrDqdZBwRTZqyB15K06+H/BFgs6RrgROA9wJI2h/YWXPbIiLmbJR7yP2qTn9I0lXAs4ELbN9ebf8unQQdEVGUVs9Dtr0VaO69ghERC9DaHnJExKgpcfbEoJKQI6JVSrxZN6gk5IholQxZREQUYpSf1EtCjohWSQ85IqIQozyGjO0iF2B9m+Ik1mjFauN3anOstiyDPDq9WNa3LE5ijVasNn6nNsdqhZITckTE40oSckREIUpOyBtaFiexRitWG79Tm2O1gqrB94iIWGQl95AjIh5XikvIktZJukPSNknn1Bhno6R7Jd1SV4yuWAdL+ldJt0naKum3a4z1RElfkXRzFevddcWq4o1JulHSZTXH+Yakr0m6SdKWmmPtLelSSbdXf2c/WVOcZ1XfZ/fygKS31hTrd6p/D7dIuljSE+uIU8X67SrO1rq+T2st9ry7afMWx4C7gMOBZcDNwFE1xToBOA64pYHvdQBwXLW+EvjPGr+XgBXV+lLgOuCFNX633wUuAi6r+f/DbwD71f13VcX6G+CN1foyYO8GYo4B3wEOqeHaBwFfB5ZXn/8eeF1N3+No4BZgTzoPnl0FrG7i760NS2k95LXANtvbbT8MXAK8vI5Atr9IQ1VPbH/b9g3V+g+A2+j8R1JHLNt+sPq4tFpquVEgaRXwC3Qqy7SCpJ+g88P6owC2H7b9/QZCnwTcZfu/arr+EmC5pCV0kuWOmuI8G7jW9kO2J4Av0KlaHwMoLSEfBNzd9XmcmhLXYpF0KJ3isNfVGGNM0k3AvcCVtuuK9WfA24EmXkBr4HOSrq+K6NblcOC7wMeqoZiPSHpSjfF2Ox24uI4L2/4W8H7gm8C3gfttf66OWHR6xydIerKkPYGfp1OpPgZQWkLWDNtaMw1E0grgU8BbbT9QVxzbk7aPAVYBayUdPewYkl4G3Gv7+mFfexYvtn0ccCpwlqS6SogtoTOU9Ze2jwX+B6jtXgaApGXAacAna7r+PnR+0zwMOBB4kqTX1BHL9m10am9eCXyWzrDjRB2x2qi0hDzOo3+arqK+X60aJWkpnWT8t7Y/3UTM6lfta4B1NVz+xcBpkr5BZ2jpZyR9ooY4ANjeUf15L/APdIa36jAOjHf9VnEpnQRdp1OBG2zfU9P1fxb4uu3v2t4FfBp4UU2xsP1R28fZPoHOsOCddcVqm9IS8mZgtaTDql7D6cCmRW7TgkkSnTHJ22xfUHOs/SXtXa0vp/Mf4+3DjmP7HbZX2T6Uzt/T523X0uuS9CRJK3evA6fQ+dV46Gx/B7hb0rOqTScBt9YRq8sZ1DRcUfkm8EJJe1b/Fk+icx+jFpKeUv35dOCV1PvdWqWo12/anpB0NnAFnbvOG90psjp0ki4GTgT2kzQOnGv7o3XEotOb/HXga9XYLsAf2L68hlgHAH8jaYzOD9y/t13rlLQGPBX4h04uYQlwke3P1hjvzcDfVp2C7cDr6wpUjbOeDPxGXTFsXyfpUuAGOsMHN1LvU3SfkvRkYBdwlu37aozVKnlSLyKiEKUNWUREPG4lIUdEFCIJOSKiEEnIERGFSEKOiChEEnJERCGSkCMiCpGEHBFRiP8P8QD5aTAew90AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sns.heatmap(distarr[:10,:10])"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x11e805400>"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAV0AAAD8CAYAAADUv3dIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFkxJREFUeJzt3X+UHWV9x/H3JwlIEgJ4ABWSIMEGMQJHIAaUNqCgTfwBtccqUMuPg6znVAS0to3VgmDbU7VK6Sn+iBLEKqAgtgEj4C/EWoFEEUgCSAhIlgCJioCESnb32z9mYi/r7p17s3eenZn9vDhzuHdm7vN9LgnfffaZZ+ariMDMzNKYNN4dMDObSJx0zcwSctI1M0vISdfMLCEnXTOzhJx0zcwSctI1MxuFpGWSNklaPcpxSfo3Sesk3Snp0KI2nXTNzEb3BWBRm+OLgbn51gd8uqhBJ10zs1FExM3Ar9qccjzwxcjcAuwmaa92bU7pZQdHsvUX65Pc8vaVg89NEQaAaUNDyWLNnfpUslgbtkxPFuvw+Y8kifP0Y6X/Ff+dmze2/X+tpw6b3i4P9NbAwORksQ564FqNtY1ucs6Oe77kXWQj1G2WRsTSLsLNBDa0vO/P9436Fzzd30gzs4rJE2w3SXa4kX5ItE36Trpm1ixDgymj9QOzW97PAja2+4DndM2sWQYHOt/Gbjlwcr6K4QjgiYhoO3fmka6ZNUpE7665SLoCOBrYQ1I/cB6wQxYnPgOsAN4ArAO2AKcVtemka2bN0sML3RFxYsHxAN7dTZtOumbWLD0c6ZbBSdfMmiXthbSuOemaWbPUfaQr6QCyuy5mkq0/2wgsj4i7S+6bmVnXojerEkrTdsmYpL8FriRbAHwbsDJ/fYWkJeV3z8ysS0NDnW/joGikezrw8ojY2rpT0ieBNcA/j/QhSX3kt9Z96hP/wDtPbnsB0Mysd2o+vTAE7A38fNj+vfJjI2q9tS7VsxfMzIDaX0g7B/iOpPv4/4c67AP8AXBmmR0zM9sudR7pRsT1kvYHFpBdSBPZvcYrI6LaP07MbGKq+IW0wtULkd1Td0uCvpiZjd04XSDrlNfpmlmjVP2XcCddM2uWOs/pmpnVjqcXzMwS8kjXzCyhwa3F54wjJ10za5aJPr2Qqkrv2++8IEkcgHteeXayWE2s0AvpqvS6Qu/Y7bNwS7JYPeHpBTOzhCb6SNfMLCknXTOzdMIX0szMEvKcrplZQp5eMDNLyCNdM7OEPNI1M0vII10zs4QGqv0Q87bVgNuRdFovO2Jm1hMx1Pk2DrY76QLnj3ZAUp+kVZJWfXfLfWMIYWbWpTqXYJd052iHgBeO9rnWasBf2vsdrgZsZunUfE73hcAfA48P2y/gf0rpkZnZWNR89cJ1wM4R8dPhByTdVEqPzMzGos4j3Yg4vc2xk3rfHTOzMar46gUvGTOzZolqX0Zy0jWzZqn5nK6ZWb046ZqZJVTnC2lmZrUzODjePWir9KQ7LdFQ/9oDP8TcqU8liXXAyouSxAHY+5R0d1unKhYJ6QpGuljk2M1b/miyWA/3ohFPL6SRKuGaWcVVPOmO5dkLZmbV08MH3khaJOleSeskLRnh+D6Svifpdkl3SnpDUZtOumbWKDEUHW/tSJoMXAwsBuYBJ0qaN+y0DwFfjYhDgBOATxX1rzHTC2ZmQC+nFxYA6yJiPYCkK4HjgbUt5wSwS/56V2BjUaNOumbWLF2sXpDUB/S17FqaPyURYCawoeVYP3D4sCY+DNwo6T3AdODYophOumbWLF2MdFsfQzsCjfSRYe9PBL4QEZ+Q9CrgPyQdGDH6hLGTrpk1S++mF/qB2S3vZ/H70wenA4sAIuJHknYC9gA2jdaoL6SZWbNEdL61txKYK2mOpB3JLpQtH3bOQ8AxAJJeBuwEbG7XqEe6ZtYsPRrpRsSApDOBG4DJwLKIWCPpAmBVRCwH/gr4nKT3kk09nBrRPpsXJl1JB5BNKN8aEb9p2b8oIq7f/q9kZlaCgqVg3YiIFcCKYfvObXm9FjiymzbbTi9IOgv4L+A9wGpJx7cc/qduApmZJTE42Pk2DormdM8ADouIPwGOBv5e0tn5sZGu7GUHWqoB37hlXW96ambWgRga6ngbD0XTC5O3TSlExIOSjgaulvRi2iTd1mUY17zopGo/xt3MmqWH0wtlKBrpPirpFdve5An4TWRLIg4qs2NmZtulh89eKEPRSPdk4DlV3iJiADhZ0mdL65WZ2faq+Ei3qBpwf5tjP+x9d8zMxmhggj/E3MwsKZfrMTNLqM7TC2ZmdTNeS8E65aRrZs3ika6ZWUITPemmKhi5Ycv0JHEgbYXeXS67NFmsbx58bvFJPZKqSq8r9I7d5i1PJIvVExO9BLuZWUpFtc/Gm5OumTWLk66ZWUJevWBmlpBHumZmCTnpmpmlE4OeXjAzS8cjXTOzdLxkzMwspbonXUkLgIiIlZLmAYuAe/IqmWZm1VLtKd32SVfSecBiYIqkbwGHAzcBSyQdEhH/OMrn+oA+gHN3P4i3znhxTzttZjaaGKh21i0a6b4VeAXwPOBRYFZEPCnp48CtwIhJt7Uw5V1z3lztsb6ZNUu1c25h0h2IiEFgi6T7I+JJgIh4RlLFv5qZTUR1v5D2rKRpEbEFOGzbTkm7UvmfJ2Y2IVU8MxUl3YUR8VuAiOcUHtoBOKW0XpmZbadaj3S3JdwR9v8C+EUpPTIzG4uaj3TNzGolBsa7B+056ZpZo1S8AruTrpk1jJOumVk6HumamSU04ZNuqiq9h89/JEkcgKcfS/ezKmWF3rffeUGyWPe88uwkcVyhd+KJQY13F9rySNfMGmXCj3TNzFKKIY90zcyS8UjXzCyhCI90zcySqfpId9J4d8DMrJeGBtXxVkTSIkn3Slonacko57xN0lpJayRdXtSmR7pm1ii9upAmaTJwMfA6oB9YKWl5RKxtOWcu8AHgyIh4XNILitp10jWzRunh6oUFwLqIWA8g6UrgeGBtyzlnABdHxOMAEbGpqNGupxckfbHbz5iZpRLR+VZgJrCh5X1/vq/V/sD+kn4o6RZJi4oaLSpMuXz4LuA1knbLvlwcV9htM7OEuhnpthbRzS3NazxClu9+r/lh76cAc4GjgVnADyQdGBG/Hi1m0fTCLLKh9OfzYALmA59o96HWL3LmjPksnvqSgjBmZr3RzZKx1iK6I+gHZre8nwVsHOGcWyJiK/CApHvJkvDK0WIWTS/MB34MfBB4IiJuAp6JiO9HxPfbfZGImB8R851wzSylwUF1vBVYCcyVNEfSjsAJwPDf/v8TeA2ApD3IphvWt2u0qFzPEHChpKvyfz9W9Bkzs/HUq5sjImJA0pnADcBkYFlErJF0AbAqIpbnx14vaS0wCPx1RPyyXbsdJdCI6Af+TNIbgSfH8kXMzMrUy2cvRMQKYMWwfee2vA7gffnWka5GrRHxDeAb3XzGzCylDlYljCtPFZhZo/gpY2ZmCQ0OVfvpBk66ZtYonl4wM0toyI92NDNLx8/TNTNLaMJPL6Sq0puyQu/NG/dKFuuw6b9KFitVhV6AA1ZelCTOvnPfnCQONLdC7yRV+8LUcJ5eMDNLyKsXzMwSqvjsgpOumTWLpxfMzBLy6gUzs4QqXgzYSdfMmiVGLPhQHU66ZtYoA55eMDNLp1EjXUl/SFaWeHVE3FhOl8zMtl/V53TbriKWdFvL6zOAfwdmAOdJWlJy38zMuhao4208FN26sUPL6z7gdRFxPvB64M9H+5CkPkmrJK267KE0twGbmUE20u10Gw9F0wuTJD2fLDkrIjYDRMTTkgZG+1BrWeNfvvGoqt8gYmYNMljzOd1dyUqwCwhJL4qIRyXtnO8zM6uUilfrKSzBvu8oh4aAt/S8N2ZmYzRU8fHgdi0Zi4gtwAM97ouZ2ZhVfT7T63TNrFGqvmTMSdfMGmVIDZxeMDOrqsHx7kABJ10za5Rar14wM6ubRq5e6EYTC0amLBY5MDA5Wax9Fm5JFitVwcgH77s2SRyA6TMXJouVsljkUFT90tRzefVCIikr9JpZdXl6wcwsoaqPy510zaxRBj3SNTNLxyNdM7OEnHTNzBKqeIk0J10zaxaPdM3MEvJtwGZmCVV9nW5RYcrDJe2Sv54q6XxJ10r6qKRd03TRzKxzVa+RVnQv4TJg272hF5GV7/lovu/SEvtlZrZd6p50J0XEtgKU8yPinIj477wi8H6jfai1GvDlmx/uWWfNzIpEF1sRSYsk3StpnaQlbc57q6SQNL+ozaKku1rSafnrO7Y1KGl/YOtoH4qIpRExPyLmn7TnzKI+mJn1zJA639qRNBm4GFgMzANOlDRvhPNmAGcBt3bSv6Kk+07gKEn350F/JGk98Ln8mJlZpQx2sRVYAKyLiPUR8SxwJXD8COd9BPgY8L+d9K+oGvATwKl5Jt8vP78/Ih7rpHEzs9SGuni4o6Q+oK9l19KIWJq/nglsaDnWDxw+7POHALMj4jpJ7+8kZkdLxiLiKeCOTs41MxtP3VwgyxPs0lEOjzQB8buMLmkScCFwahchC6cXzMxqpYcX0vqB2S3vZwEbW97PAA4EbpL0IHAEsLzoYppvjjCzRunhUrCVwFxJc4CHgROAk7YdzKdf99j2XtJNwPsjYlW7Rp10zaxRBtSbgj0RMSDpTOAGYDKwLCLWSLoAWBURy7enXSddM2uUXtZIi4gVwIph+84d5dyjO2nTSdfMGmXCP2XMFXrHJmWF3nnLH00Wa/OWJ5LESVmh9+mHb04WK+X3mrnz7sli9UI3S8bGg0e6ZtYo1U65Trpm1jATfnrBzCylwYqPdZ10zaxRPNI1M0soPNI1M0vHI10zs4S8ZMzMLKFqp1wnXTNrmIGKp92iasBnSZrd7hwzsyqJLv4ZD0XP0/0IcKukH0j6S0l7dtJoa2HK7265b+y9NDPrUN2rAa8ne3DvR4DDgLWSrpd0Sl7CZ0SthSlfO21uD7trZtZe3Ue6ERFDEXFjRJwO7A18ClhElpDNzCql6iPdogtpz6kRFBFbgeVkJSmmltYrM7PtNBjVvpBWlHTfPtqBiHimx30xMxuzWq/TjYifpeqImVkv+DZgM7OEfBuwmVlCtZ5eMDOrG08vmJklVPfVC2ZmtTLhpxdSVel1hd6xS1WhN6VJKrr/p3eaWnl46t5/lCxWL/hCmplZQp7TNTNLaMJPL5iZpRS+kGZmlo5LsJuZJeTpBTOzhDy9YGaWkEe6ZmYJ1XrJmKQdgROAjRHxbUknAa8G7gaW5g81NzOrjLrfBnxpfs40SacAOwPXAMcAC4BTyu2emVl36j69cFBEHCxpCvAwsHdEDEr6EnDHaB+S1Af0AXx4z5fztl336VmHzczaqXrSLboxfVI+xTADmAbsmu9/HrDDaB9qrQbshGtmKUVEx9t4KBrpXgLcA0wGPghcJWk9cARwZcl9MzPrWtVHukU10i6U9JX89UZJXwSOBT4XEbel6KCZWTdqvXoBsmTb8vrXwNWl9sjMbAwGo9oPd0z3sFEzswR6OacraZGkeyWtk7RkhOPvk7RW0p2SviPpxUVtOumaWaMMER1v7UiaDFwMLAbmASdKmjfstNuB+RFxMNkswMeK+ueka2aNEl38U2ABsC4i1kfEs2SLB45/TqyI70XEtlIytwCzihp10jWzRhmK6HiT1CdpVcvW19LUTGBDy/v+fN9oTge+WdQ/P3vBzBqlm9ULEbEUWDrKYY3Y/EgnSu8A5gNHFcV00jWzRunh6oV+YHbL+1nAxuEnSTqW7D6GoyLit0WNlp50U1XpdYXeeklVpXco4fKhmTvvnixWygq9z2z8QbJYvTDUuzvNVgJzJc0hewzCCcBJrSdIOgT4LLAoIjZ10qhHumbWKL26OSIiBiSdCdxAdlfusohYI+kCYFVELAc+TvYgsKskATwUEce1a9dJ18wapYcjXSJiBbBi2L5zW14f222bTrpm1ii1vw3YzKxOBmNwvLvQlpOumTWKC1OamSVU60c7mpnVjUe6ZmYJ9XL1QhkKk66klwBvIbszYwC4D7giIpq5at/Maq3qqxfa3hYk6SzgM8BOwCuBqWTJ90eSji69d2ZmXRqMoY638VB0L+YZZLe3/QNZmZ55EfFBYBFw4Wgfan1yz9VP/bx3vTUzK1D3wpTbzhkkqwA8AyAiHpLUthow+ZN77prz5mqP9c2sUeo+p/t5YKWkW4CFwEcBJO0J/KrkvpmZda3Wqxci4iJJ3wZeBnwyIu7J928mS8JmZpVS+3W6EbEGWJOgL2ZmY1brka6ZWd1UvQS7k66ZNUrdL6SZmdWKpxfMzBKq+h1pTrpm1ige6ZqZJVT1Od2ubplLuQF9TYrjWPWK1cTv1ORYddrS1MHePn0Ni+NY9YrVxO/U5Fi1UeWka2bWOE66ZmYJVTnpLm1YHMeqV6wmfqcmx6oN5RPeZmaWQJVHumZmjeOka2aWUOWSrqRFku6VtE7SkhLjLJO0SdLqsmK0xJot6XuS7pa0RtLZJcbaSdJtku7IY51fVqw83mRJt0u6ruQ4D0q6S9JPJa0qOdZukq6WdE/+Z/aqkuK8NP8+27YnJZ1TUqz35n8fVku6QtJOZcTJY52dx1lT1veptfFeKDxsMfVk4H5gP2BH4A6yumxlxFoIHAqsTvC99gIOzV/PAH5W4vcSsHP+egfgVuCIEr/b+4DLgetK/m/4ILBH2X9WeazLgHfmr3cEdksQczLwKPDiEtqeCTwATM3ffxU4taTvcSCwGphGdsfrt4G5Kf7c6rJVbaS7AFgXEesj4lngSuD4MgJFxM0kKjkUEY9ExE/y108Bd5P9j1BGrIiI3+Rvd8i3Uq6WSpoFvJGsrFMjSNqF7AfyJQAR8WxE/DpB6GOA+yOirEquU4CpkqaQJcSNJcV5GXBLRGyJiAHg+8BbSopVS1VLujOBDS3v+ykpOY0XSfsCh5CNQMuKMVnST4FNwLcioqxY/wr8DZDiqdEB3Cjpx5LKvNNpP2AzcGk+bfJ5SdNLjLfNCcAVZTQcEQ8D/wI8BDwCPBERN5YRi2yUu1DS7pKmAW8AZpcUq5aqlnQ1wr7GrGmTtDPwNeCciHiyrDgRMRgRrwBmAQskHdjrGJLeBGyKiB/3uu1RHBkRhwKLgXdLKqtG3xSyaadPR8QhwNNAadcWACTtCBwHXFVS+88n+41xDrA3MF3SO8qIFRF3kxWw/RZwPdkU4UAZseqqakm3n+f+VJxFeb8GJZWXrP8a8OWIuCZFzPzX4puARSU0fyRwnKQHyaaBXivpSyXEASAiNub/3gR8nWwqqgz9QH/LbwdXkyXhMi0GfhIRj5XU/rHAAxGxOSK2AtcAry4pFhFxSUQcGhELyabw7isrVh1VLemuBOZKmpP/9D8BWD7OfRozSSKbI7w7Ij5Zcqw9Je2Wv55K9j/cPb2OExEfiIhZEbEv2Z/TdyOilNGTpOmSZmx7Dbye7NfYnouIR4ENkl6a7zoGWFtGrBYnUtLUQu4h4AhJ0/K/i8eQXVcohaQX5P/eB/hTyv1utVOp5+lGxICkM4EbyK7mLousGnHPSboCOBrYQ1I/cF5EXFJGLLJR4V8Ad+VzrQB/FxErSoi1F3CZpMlkP1S/GhGlLudK4IXA17N8wRTg8oi4vsR47wG+nP/gXw+cVlagfN7zdcC7yooREbdKuhr4Cdmv+rdT7i26X5O0O7AVeHdEPF5irNrxbcBmZglVbXrBzKzRnHTNzBJy0jUzS8hJ18wsISddM7OEnHTNzBJy0jUzS+j/AIgaNXMFJHdAAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"sns.heatmap(simmat[:10,:10])"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA+gAAAGQCAIAAACyL902AAChY0lEQVR4nOzdd1iT1xcH8G8S9hIcCFr33gM3jqo42oKjFkcrriquFrS2Uls1WmuL/dWWLhWtVlq1FjdareLGDYrWreCqIqDsPZLz++NiRFyQvDEEzufp00djcjnR5H3Pe99zz5URERhjjDHGGGOlm9zQATDGGGOMMcZejhN3xhhjjDHGjAAn7owxxhhjjBkBTtwZY4wxxhgzApy4M8YYY4wxZgQ4cWeMMcYYY8wIcOLOGGOMMcaYEeDEnTHGGGOMMSPAiTtjjDHGGGNGgBN3xhhjjDHGjAAn7owxxhhjjBkBTtwZY4wxxhgzApy4M8YYY4wxZgQ4cWeMMcYYY8wIcOLOGGOMMcaYEeDEnTHGGGOMMSPAiTtjjDHGGGNGgBN3xhhjjDHGjAAn7owxxhhjjBkBTtwZY4wxxhgzApy4M8YYY4wxZgQ4cWeMMcYYY8wIcOLOGGOMMcaYEeDEnTHGGGOMMSPAiTtjjDHGGGNGgBN3xhhjjDHGjAAn7owxxhhjjBkBTtwZY4wxxhgzApy4M8YYY4wxZgQ4cWeMMcYYY8wIcOLOGGOMMcaYEeDEnTHGGGOMMSPAiTtjjDHGGGNGgBN3xhhjjDHGjAAn7owxxhhjjBkBTtwZY4wxxhgzApy4M8YYY4wxZgQ4cWeMMcYYY8wIcOLOGGOMMcaYEeDEnTHGGGOMMSPAiTtjjDHGGGNGgBN3xhhjjDHGjAAn7owxxhhjjBkBTtwZY4wxxhgzApy4M8YYY4wxZgQ4cWeMMcYYY8wIcOLOGGOMMcaYEeDEnTHGGGOMMSPAiTtjjDHGGGNGgBN3xhhjjDHGjAAn7owxxhhjjBkBTtwZY4wxxhgzApy4M8YYY4wxZgQ4cWeMMcYYY8wIcOLOGGOMMcaYEeDEnTHGGGOMMSPAiTtjjDHGGGNGgBN3xhhjjDHGjAAn7owxxhhjjBkBTtwZY4wxxhgzAkaQuO/ahbNnAWDbNgDYvRvJyQCweTPy8gwXFmOMMcYYY6+QESTuu3fD1xf5+QgJAYDQUCQlAcDWrcjNNWxojDFm9ApPjvDMCGNlDH/ByxgjSNwBjBmDxYsf//bzzzFxIo4fN1xAjDFWVhSeHOGZEcbKGP6ClzEmhg6gWDw88PHHuHOn4LcLF6JOHYwahdxcXLqEdu0gkxk0PsYYM2aFJ0c+/xy2tjwzwljZwV/wssQ4ZtwB+Pvj9GkAuH4df/6JlBQA+PlnXLtWUPvOGGNMOx4euHwZN28CwMKFCAxE586GjokxJhH+gpclRpC4BwSgcmU4ORUk69HR+Pxz/Pcffv8dY8bg/Hn+/DHGmK40kyMaMhlycxEeDiIDxcQYk4jmC65WIyoKRPwFN1YyMrZ/sQ4dEB6OkyfRoQOGDEGXLujXD82bGzosxhgzfm5usLHB8uVwdAQAf3/UqAFrawwaZODAGGOSaNwYV6/i2jU0aMBfcKNkBDPuRVhZAUBmJgBs2oQZMzhrZ4wxCWRk4MAB7NoFe/uCR957j+9qMlam1K8PAP/+C/AX3DgZa+KelWXoOBhjrGy5cAFqNZo0gZlZwSPTpqFKFTx4YNCwGGPSadkSAM6fB/gLbpyMo6tMYYVn3BljjElFTMKJ87qwaZOhYmGM6UWLFsCjLzt/wY0RJ+6MGb2UlJRr167VqFHDycnJ0LEwIyYm4cR5nTH9uXnz5tmzZyMjIyMjI5VKZY0aNapWrWrooMoL8QUXX3ZmjIwvca9R407nznJABdQydCyMvWr5+bh5E9HRpy5eDLt69eq1a9euXLkSFxcHwMvLq2HDhrNnzzZ0jMxYPT3jzpju8vNx5QoiI3HlynfHj++IjIxMFrt3AgDCwsJcXFz27NmjUCgMF2M50rgxzM1x4wbS02FjY+hoWMkZX+Kemfn98eMBnp7fAdMNHQtjerFrF5yd0bo1tm1DnToID8e1a7h2DVeu4MYN5Oaia9cLR458rHm+lZVVzZo1165dK5PJevbs6erqasDgmfG6cgXgGXemrcIHrho1cOoUIiMRGYnz55GdDQDt26eHhx8A4OTk1KZNm9atW9etW3f27Nn79+9fsGDBvHnzDBt/OWFigtat1ampWVeuZLVrV9nQ4bASM77E3crKCkAm18qwsmv3bkRGYt8+hIQgLw9//PH4j2Qy1K6NZs1atmnzYaNGjRo1atSwYcMaNWrIZLLZs2cvXLjQy8vr3Llztra2hgufGaW7d+8mJTXp2PHdatUCDR0LM0qFD1z37mH37oLHZTLUr482bdCjx4D589u3adOmcFFf/fr13dzcFixY4Orq2qdPH8OEXs40bDjmjz/+iIxc3q7dBEPHwkrM+BJ3S0tLAFncVoaVaZodqnv3BhEaN0bDhgX/WVoCaAe0A5CXl3fz5s2TJ0926tRp3rx5e/bsCQ8P/+ijj1asWGHY+JnROXfuXG5uurV1lKEDYUZMc+Byd4ezM9q0QevWaN0adnbiz1urVC2uXbt24MCByMhIhULx9ddfv/7667Nnz54/f/7IkSMjIyOrVatmyDdQPrRo0QLAea5zN07Gl7jzjDsrDzw88PHHuHMHK1di9OiCB+PjC8pm/vsv4+zZd69cuXLz5s28vLwGDRpcu3bNxMQkKCjIxcXl119/7d+//5AhQwz6DpiR+ffffwG05Ap3poPCBy4hOxsXLhTUzNy5k33wYOWMjAzxR1WqVPn6668BzJ0798iRI/v27Rs1atzu3f9wrbu+ia+5+Mozo8OJO2Ol0YkTGD368W52ublwdkZiYsFvzcysVardKlWOXC6vU6dOo0aNiEgmkzVp0uTrr7+eNm3awoXru3UbIja/ZKw4xPRbC65wZzpITMSsWejQoeC3KhWqVEF6esFvZTILO7vKdeo4tm7dWtS4q9VquVwul8vXrVvXt+/YBw/WffEF5s83VPjlhfiac+JupIwvcRelMoXXpDNWxnz1FaytYW5esKILgJkZ7OxAhIYNC8pmmjTZVr9+tYYNG5qbmxd+rY+Pz8WLtdatG+jtja1bX33szFjxjDvTUUAALCyQk/P4wKVQoHlzpKcX1My0aYPWrS86OFg//VpHR8fvv/+7Tx9cuICuXcG17npVrVq1KlWqPHjw4O7du6+99pqhw2ElIyMiQ8dQAnv37v3ggw+Sk5MTEhIGDx48ffr0zrxXLytzYmJQvTqcnHD//uMHU1M1daIvce8eWrZEYiJ+/RXvv6+nGFmZkpOTY2trq1arU1NTxV1NxrQgl4MIKhXkj7ZlV6sf//qlvvgCSiUcHREZCa5116uePXsePHhw06ZNb7/9tqFjYSVT7O+Toe3du7dz5859+vS5evVqXl6eWq3esGFDly5dunTpsnbt2pycHEMHyJhkxP0ke/snHixm1g6genX89BMA+PoiipcasudLS0vbu3fvp59+6urqamZmZm1tzTczmdby8kAEU9MnMvXiZ+0AZs9G376Ij8e77yI/X/IAWYHbt2/fuHHD2dnZ09OzXbt2n3766d69e7M1N0pYKUel3rFjx3r37i2irVy5sr+/f2Zm5vXr1z/66CMHBwfxuKOj4//+d+rOHUPHypgUjh0jgDp21GmQ4cMJIFdXys+XKCxWJsTG0ubNaR9++GHLli3lhbIq0UK0SpUqf//9t6FjZEYpPZ0AsrbWaZC4OKpWjQCaPVuisNiTVq9erfmym5g8rpe2tLTs3bv3Tz+FnjhBeXmGjpI9X6lO3M+dI09P6tDBA0DFihWVSmVKSkrhJ2RnZwcFBbVu3VqhMKtePU8uJzc3CgkhtdpQITMmgZ07CaB+/XQaJCmJatQggL7+WqKwmPHYuZMiI4mItm4lIoqJoeBg8vEhFxeSyUguV1WsWAmAiYmJi4uLn59fSEhIdHT0m2++CUAmk/n4+OTk5Bj0HTDjk5hIADk46DrOgQNkYUH/+58UMbFCkpOT33vvPZGmDxkyJCEhISMjIzQ01M/Pz8XFRVzGt2+/XVx9ubmRvz9FRHBCVeqU0sT97FkaMIBkMgLIxeXsvHnzRMp+/fp1Ly+vnTt3Fnn+iRPXRowgMzMCCKAmTejnn+nJJJ8xo7FuHQE0bJiu44SGUt26dOiQFDExo+LrS927U14ejRtHnp4FB0bxn7U19elDAQHrDx8+nJ2dXfhVarU6ICDAzMwMQLt27aKiogwVPzNG9+8TQFWr6jrOzp0UGkpEtHUr/fMPJSUREW3aRLm5uo5cnh04ECbWodrZ2QUFBT39hAcPHmzYsMHXN7VhwyeOGI6O9OOPT0wEMMMqdYn7pUvk5UVyecEJxseH4uKIiG7fvj1+/HhTU1MAnTt3fuZrY2NpwQKqXr3g0+biUnTaiTGjEBioAsjbW6dBxIc/N7fgw8/nv3LF15dWrSJ/fxo3jj77jGxsyM2NlEoKDaUnc/XHQkJCzp07R0THjh2rVasWgH79tmzc+ErDLmN27tz57bff1q9fX6lUJomvX5l2+zYBVKOGruMUvuycMYNu3CAi8vKi9HTdYyyPcnNJqaTWrVMUCrOOHTtev36diA4cOPDxxx8/7yX371NwMHl7U61aZGJCkyc//hdhBlcqEndNer1xI1WtSgBZWtL06QUp+717NHNmvOgCaWJiMm7cuJs3b75gNJWKQkLIzY1+/PGJ7z9jxsLf31+hMPvssy91GaTIh5/Pf+WKry89eECjR1PPnpSe/vJ1DlFRUba2tpaWlkuXLiWixMREb+/F4h7mlCmUlfUqYi5jduzYYW5uLpPJRGWCvb29n5/fvXv3DB2XHkVF3apRo0e3bmN1HKfwZeeMGTRiBHl7U/36fODSxqVL1KYNAWRiQj/8cCY/Pz8nJ8fPz08UxuzYseOlI9y+/cS/CDO4UtFVZvdu+PoiPx87d2LWLEyZguvX8d13kMvx6aeoXx/ffFPFxeUdT0/Pixcvrly5snbt2i8YTS6HhwdCQ/Hhh0ChHZgZMxZJSUkqVa6Nja5fzyIf/s8/x8SJOH5cx1HLEbVabegQdOLvj9OnYW2Nl25FWb169bFjx2ZlZU2ePFm0hwsM/GjFClhZYckSuLjgwoVXEXCZsWvXrnfeeScnJ2fixIl79uxxd3dPTk5etGhRnTp1hg4dGh4ebugA9SInJ+O//w49fHhC96E8PHD5Mm7eBICFCxEYCO78rIXff0f79oiMRO3aOHAAPj5trl+P7ty586JFi+RyuVKp7Nev30sHqVkTePJfxEiRUXU/f4FSkbijUIbh64tffkH16ti2DXXqYNEiZGdjyBAsW7Y6ODi4YcOGJR25DHzaWHmTkpICoEKFCjqOU+TDz+e/4ktISPD19R09enS/fv2uXLli6HBKrEYNzJqF2FikpBTr+RYWFj/88MP69evt7Oy2bNkyZkzCqVMYNQqnTqFZM1y6hA4dsHy5noMuK3bu3Dl48ODs7Ozp06cvXbq0T58+27dvP3PmjJeXl2hk3KFDh65du27fvr3MZBJCbm4uALFGQkchIfj6a5w+/cSDj+5esJdLSYG7O0aPRkYGxo3D+fPo2hW//46336597dr9WrVq7d+/f968eYW7yrxUx444cgQZGfqLWl8OHz48c+bMxo0b+/r6xsTEGDocnRl6yp/oybu6GvfukaUlublRRIROI+/bR6tWka2t7mEy9ooMHz4cwNq1a3UZxNeX1qyhDRvIzo7oyVKZjAwpoiyj8vLyvv/+e3HVpFAoAFSqVOnw4cOGjqtk2rUjgPbtK/ELo6KiRo7cDpCZGS1eTGo1pabSyJEEUKVK9NdfvGroJUSFDICPPvro6T+9efOmn5+f5pq8VatWQUFBuWVl0cnJkycBtG/fXsdxFi8mgN57T5KgyhdN4fGWLdSjB1WoQOvWERHFxZGHR8HyPz+/y6mpqVoM3qEDAbRnj5QBvwK7d+8uvKmcpaXlpEmTRKG/kSotifu9e3T8eNH0WpK+7K1bE0Dh4RIMZVgxMTGGDsFg8spZU9n+/fsD0LTTfuedd7744ouHDx+WaJDz58nWlkxM6MwZPYRYRh04cKBFixbi+O7m5hYRETFw4EAAZmZmz+zDUDolJ5NCQWZmWl6h5eWRUlnQIcDNje7fJyJatYq2b+dVQy+xceNG0UHhBcv+iCglJSUgIECz1byTk5NSqUxMTHxlcerJkSNHALi6uorf7ty508XFpaRXJidOkJkZyWS0ZYtegizbCn9Db9+mu3eJiM6fJyengk6dwcHaDz5zJgH02WdSBfsq7Nixw8LCAsCECRPCw8O9vLzEdIxcLnd3dz9+/LihA9RGqUjciWjWLAJowQLpR3Z3J8DoDwFZWVl169bt2bNnfHy8oWN51YKCgnr06FG3bl0vL6/AwMALFy4YOiK969SpE4CjR48S0dmzZ8XZ3crKaurUqcXs0JeQQPXqEUAjR+o51rLi7t27Xl5e4q+6Xr16wY/Ob/n5+X5+fuJxHx8flUpl2DiLY8cOAqhrVyKiH36gfv3on39KPEhICFWqRBUr0u3bjx/kNWovsGHDBpG1f/LJJ8V5fk5OTlBQUPPmzcWny9bW1sfH53bhv25js3//fgA9H906f+edd8Rbc3Z2nj179p1iTMUlJlLt2gRQ8f4KWVHP/IZmZ1OrVtSzJ/33X4kHVKnou+/I05Py8ws2GOnSRcJ49euvv/4SX8kpU6aoH7Wjj4qK8vHxEf1OxHVmSEiI2qia1ZeWxH3pUgJo/HhpRjt5krp2pcmTiYgmTiSAfv5ZmpENZd68eQCaNWtW3uae//77bwBF6vBee+214cOH//PrrxQZWfb2BVWpVDVr1gTw+++/i0fCwsLc3d1FewoxTyBy+ufJy6NevQigtm25KublcnJyAgICbGxsxNWRUqnMeqqLyvLly8UJwNPTMzMzU/cfev/+/XXr1t26dUsfJ4xPPnm88WT//gTQ+vXajHPnDoWFPfHIM8saGREFBweLw9TMmTNL9EK1Wv3333/37NlTHNxMTU1Hjhx5Q5S1GZvNmzcDaNGihfgGpaSkBAYGtmzZUrw1uVzu5uYWHByc/5yDtlpNAwcSQJ06EW//pZ3nfUNjY0nrOYf69QmgiAhKSyNTUzIxIa0KbV61NWvWiK+kn5/f038aFxenVCorVqwoPpwtWrQwoqK10pK4iyu5vn2lGU3sGC8K7b74ggCaNUuakQ3i9u3b1tbWMpnswIEDho7llVKpVK1btwbw3XffRUREBAQEeHp6VqlSRXzT9r3+OgFkY0OuruTnRyEhlJxs6JB1lZ6ePmjQIJlMJspk27ZtGxQUJK7Wrl+/XnieQNyDfuYpcOpUAsjJSZv5lfImJCSkbt264q/U3d391q1bz3vmnj17RGlyp06d4kSrWq3cuHFD/DvK5XJ7e/upU6dqPdTztG9PAIWGUl4e2doSQFI1IRRpwf37BQsnmPDXX3+9IEUopsjISC8vL1NTU4VC4efnl2JsOwjGxcW1aNHC0dERQIUKFby9vTV3RyMiIry9vTXHrurVqyuVyqe/RP7+KoAqV+YDl/bEN/T6dbK2lmzM8eMJoG+/JSLq2JEA2r1bssE1EhIS1qxZM378+LVr1+o+Qbls2TLR7/LFX8m0tLSAgIAaNWpobg0plcrkUp9IlJbE/eJFAqhRI2lGu3OHAHJ2JiJatYoAGjVKmpENYtCgQQBGlr+ih6CgIAC1a9cuvL+jWq2+dOnS8uXLYz7+uKAcRPOfqSl1704qFS1dSh9/TBcvGjB4LcTExLRr1w6Avb29l5dX5cqVxdGkfv36v/zyS0ZGBj2aJ6hUqZKmqCMgICCj0Lz6mjWn5HKVhQUZZ/Heq3Pt2rU333xT/DU2btx4dzHORefPnxc7E9WtW/fSpUsl/YlnzpwZNmyYpsKya9eu4vJszpw5Wr2DZ0tJybWwIDMzSk+nEycIoMaNJRu8e3eaMIE3pX7C+vXrRdY+b9483Ue7detW48aNAaxevVr30V6ZmzdvNmjQQNwObdWqlfhayWSyPn36bNq0SeRhcXFx/v7+derUEX9qbm7+7rvvavalOnToUIMGzRs2PL9rlyHfSBmQnExmZmRl9dyt1kpqzRoCyN2diMjPTy9l7ikpKe3bt9dselC7du0ffvghXdu+/f/73/9kMplMJlu8eHFxnp+bmxsUFNSsWTPx0+3s7Hx8fErzlgulJXFPTy/Yd0mS+8Z5eaRQkFxOubm0ezcB1KuXBMMaxJ49e0T5Y2n+GOlDTk6OmAr9448/XvS82FgKCSE/P3J1JTMzcnWln3+m+fPp3j0aO1abthoG8u+//4oKmbp1616+fJmIsrOzg4KCxFlcTGJpjiZinkCzocGKFSvEIGFhYWZmZh07DgwKKl8lVSWSnp6uVCpF0mxvbx8QEFD8CZ779++3b98egIODw75if7oKFzuZmpp6eXldvHiRiEJCQkQFzqJFi7R8M0/ZuXOntbXjiBE/ENFPPyWbmtLEidKMfO1aQW8ZY6jzf0X+/PNPkbXPnz9fqjEDAwMB9OvXT6oB9e3SpUtioa2Li4uYRxdT7NbW1uIA5eTk5OfnJ25nqVSq0NBQT09PExOTunXrikUjcXFx1apVAzB79lwDv5kyoVUrAkiqO/R37xJAFSpQfj7t3Rvn6rrkjTd03WOrsIyMjO7du4tz3w8//KBJoMUpr6RtOfz9/cVF448//liiF6rV6pCQkG7duomfLprPlM68q7Qk7kRUqRIBpMMt6Cc4OxNAd+7QxYv5zZqlDRwoRYeaVy4nJ6dRo0YA/ve//xk6llftu+++E5VnJVgOmJZGUVE0aFBBZffZs/Tpp/qLUEK7d+8WZRidO3cucgdZpVKFhIR06dJFM02lSfvy8vLWr1/fv39/UXV9+/ZtcZ96xowZhnkbRkJ025TL5d7e3g8ePCjpy7OyssQIJiYmYp/R5xFngs6POufb2Nj4+PgUWaK3du1auVwuk8mWLVtW4nfyLDNnzgTw+eefE9Gbb75pbe24cWPYS19VHEuWEEDDhkkyWFmwbt06cf/kiy++kHDYxMREMzMzExOT2NhYCYfVk1OnTol7gz169ChS3iNq3DVtmorUuN++fVt0WVWpVH379hUjlLdFXHoyY8bjVS6SGDz4Dycnl/Dw8LS0NFNTUxMTE+0aSj4tOztbdFGrUaOGWNohDptubm6FT3lXrlwpzmhz5swBoFAofvvtN61DOn36tGg+Y2Ji8v3335fChgSlKHHv2FFVu3ZuRIQ0mxoPGjTHycnl+PFTycnJAKwlLPh6hebOnQugSZMmxrJmQiqpqakiB925c2eJX/zhhwWNMDZtojVrJI+tgFotVUFOYGCgmLQbOnToCxY+RkREFGllFRoaqvnTzMxMUWbTp08fPvm9QHJycosWLUxMTPbo0I5YrVYrlUpxXnlmqxlx77Vp06biOZUrV1YqlQkJCc8cbenSpTKZTC6X/71hg9YhaXTs2BHAnj178vLy7OzsAEg1aTRkCAG0fLkkgxm9lStXiiLaBXrohubu7g7g51LfVGHfvn22trYABgwY8OJjV+Ea93r16vn7+2s6pImvkqOjY+mc3TRGYtFgp06SDThhwgTNBKJoevaPFp2qnpKTkyM+6lWrVhX3mQt7+pR34sSJ5w2lVqunTZsmbmlukOJAeuXKFVGSqkVVpL6VosRdVHJL8jdeZDRxZJF8wUFOTk5cXFyELhtEPV9KSoqPj4+9vb2ZmVn37t2fd8ovqz7//HMA3bt31+bFFy7QsGE0cSJNmUIffUR6KpmUoiBH02pQJpP5+fkVp8GIaGWl2U5CrF7Nzc0Vc8ANGzYsA92g9U2Uto/XuYnVypUrRaHLkCFDNMsM0tPTC692qlWrVpFFCM8UEBDQvVWrlIoVaft2XUIqPB8mdsNpJNHKIZWKKlcmgKKjJRnPuN26dUvsD6qne6Fr165FoYboEpNo0mHz5s2iQ/aoUaOKM1kQGxtbpMbd09Pzf//7n0KhkMvlulxIsyIyMsjcnBQKkupssGbNGgDu7u5E9OmnnwKYpXvHj7y8lT4+4prtBcmxOOWJTxqe070xPz///fffFx+qLdI1/x4yZAgKtXcrPUpR4u7r6wvgW7F0WWdTp04FEBAQQESi2kTa/t8PHjx4/fXXGzZsKJPJevfuff78eQkH3759uygZNDU1FZW4NWvW3Lt3r4Q/ojSLi4sT11rHjh3TfpS0NPr114JFyiUvh3i5IgU5H31EixfTzZslCTDNw8NDHGtKemiIi4ubO3euZvWqmBioUKFCKZwbKIWuX79ubm4ul8tPnjyp41B79+61t7cH0KFDh0uXLhVeN9y8efMS9RfLnT27YKGPDqWpu3btAtClSxciWrRoEYCJElW4R0QQQLVqSTKY0du0aZOjo2OfPn30NH56erroJHazJIeU4pJi0mHJkiXihoOPj0+JWpqqVKodO3a89dZb4uXi/xKuEGBCjx5S7mATExMjTjH5+fniINNJx/n8/HwaPpyAeW+8cfbs2Zc+PTY2VqlUOjg4iKNry5YtNf3W8vPzR40aBcDKykray7+vvvoKgK+vr4RjSqIUJe6iptnHx0eS0cTfuNgIo1evXgCK0zWimCIjI0VziUqVKomb0aamptOnT9cskNdaYmKit7e3+Gi2adPmzJkz0dHRYrWETCbz9vbWep21EZk0aZKYxdR1IJWKunfPr1ZtuWjpL63CBTmBgWRqWtDZpmlTUirpqbt+Rdy7d69t27YAKlasePDgQe1C0KxeNTc3r1ix4nrtmnWXS6IQvF27drrXL168eFFMIoprbDEntH37dm0atE+fTgBZW9ML+/S/wIwZMzSTYX/88Yerq2uwLjslFvLjjwnm5pJttWHs9u3bJ25w6e9HiHtoEq5afqzIpIOfH23bRmlpxR9As/7P399f6yhu3brl5+dnY2OjUCiMrvdl6bdoUWbnznfmzQt9+VOLp2HDhk5OTtHR0RKUuavVNGECAWRnR6dOFf91qamphbccrl279uLFiwcMGADAxsZm//79WsbzHLt37wbQVWxlV5qUosR906ZNAAYOHCjJaL///juAd999l4jEhoirVq2SZOTt27eLZL1Tp073799PSEjw8fERZVgVK1YMCAh43u4SLxUcHCyalFtaWvr7+2vGycvL8/f3F3dm69atKxb0lFXXrl0TnYwlmTzOu327RY0aEv7rExHdukUdO9KWLQUFOR98QPHxtGEDjRhR0DRb/Ney5ZpvvnnmrZhz586JUor69esXc83NC2i63UtSdFhOpKWlVa9eHYAua5g07t27Z2Fh4ezs3LFjx5CQEO0HUqvp/fcJIHt7iows5ovy8/MjIiL8/f3d3NxMTEysra0l6UtYRN++fa2sKm/adETykY1Rfn5+1apVAUQW+5+ppLZt2wagdevW0g9deNLhm28KjlcmJuTqSv7+FBHxguZuarVaXBwqFIrlUix3EEsytFnLxF7o6NGjAJo0aSLVgIW7JjRo0EAmk/Xt23ft2rUlXt+vVtOUKQSQlRUdOqRFJNnZ2b/++qum35pCoXBwcDhVkguAYkpISJDJZFZWVlondXpSihL38PDwIsepZcuWBQYGajcllp6efv/+ffHrWbNmAVAqlTpGqFar/f39xa29ESNGFF6Lc+nSpX79+mlmykuaW9+/f//tt98WL+/WrdvVq1effs758+fbtGkjPqZ+fn7ZUvVoLWU8PT0BeHt7SzXgzz//LObG5s+fL8G8zpUr9NprBNCgQURUdJoqO5tCQ8nHhxwd852c5DIZgDp16vj4+ISFhYkp2J07d4pCIFdXVy1amjyTOJXqI10rw0TVpqOjo+43yrZs2QKgcePGEmyDmpdHgwYVlHi9uEzi8uXfliwZOHCg6EekOYeJG4BSTbQLOTk5onLDKPqcvBqTJ08G8JnkHa0fycnJEds6XpRwP4rsbPrjj8ergD74gK5coQULyNWVFIrHkw61as384IPNmzcXmVLNy8sbO3asuLkk1Wq0zz77DNwISw/y8vLEkeE/qbezEp1bNFUrcrncxcXFz88vNDS0WJWBM2cW1ATq1qxZpVJt3bq1bt26MplMf+c+0XNZyu+gFEpR4h4fHy8+DeK3Dx48EBPbbdq0CQvTvp1ZcnJyly5dHB0d7ezslEql9qlSenrMlCkWJiYKheJ5hfghISGalTfu7u7FLE8MDg4WB+gKFSoEBAS84EIlNzdXqVSKc3Pz5s3PnDmj3VsptcLDw2UymaWlpVTHmrNnz4qiJvGXVqlSpa+//lr7cqNLl6haNQKoW7eXbPqcm3tv/35vb2/RG0dzX69Xr17iwm/MmDE50m3qvXHjRhhV4+fSQK1W9+jRA8D06dN1HEqsdv3++++liIsoJ4f696devZ5RvRAfT8HB5O1NtWoRoGzfXny06tat6+3tHRwcnJiYKFpRKRQKCe8y7d+/H0CrVq2kGrAMEH8ndevW1d+PEOvt5s6VqLV5bi4NGEAALVxI9NSkQ1oahYSQtze99lp6kyaa60AXFxelUhkREZGVlTV48GBRkyBhJbEoOtLLjYVyT9SQSLiTl1qtnj59uvhgfPnll999912/fv00y0YBbHnjDRo6lFatouc1X9+xgwAyN5eqacSKFSvwaNWsPrzzzjsAgoKC9DS+dkpR4h4bG2tpaWlnZ/f555+LnKZIHvyC3cifZ/fu3eKGuKYLh6Wl5YQJE0q8UPX2bWrdmoBzffu++KZeZmamv7+/jY2N+KF+fn5pz68dvHHjRp8+fURgb775ZpEGz89z/Pjxhg0bink1pVJZ2m7i6KJ3794STmIVLmoKCQkRSx0AVK5c2d/f/6WNPoo6c6agrcbrrxe/HlSlUoWFhfn4+IjPYcWKFRUKhe43f4q4d+8eAHt7+1LYcbY0i4yMFM16//33X60H+e+//xQKhZmZmaa9nQTS0ykn54nuH7/8Qi1bkkz2eFq0atVzn3zy22+/PX2Vq/UWJM8jujx99NFHkoxWNqhUKmdnZwD6m0AJDQ0F4OTkJEFnhUdrAcnenl4csFodc/bsV1991a1bN9GmVhD3CR0dHaVtpJaVlWVpaSmTyaT8+jAiIvrhhx8AeHl5STJaXl7emDFjxP2WjRs3ah7PzMzcuXOnj49Pw4YNE9u0eXyAqluXfHwoNJSKTFF99plka2aJ4uLiFAqFubm5VH3li/j6668h3dpLqZSWxD02NlbT8xhAkyZNxDV9kTxYqVRmZWUVZ8DMzEw/Pz8xu9mpU6fr16/v27fPw8NDPCK2Yj6xZ0+xdmo9epScnAig+vWL2T/r7t27Xl5eYqPE1157LSgoqMg9dLVaHRgYKN6Xg4NDYGBgcYZ93rt7ZmmN0VGpVLNmzZLL5RMnTtTxS/i8oqawsDAxyQqgSpUq/v7+xfw4qU6cIHt7AsjDQ7uNpFUq1YEDB0xMTBQKhZjyz8zMLM5q+mISRfOl7Y5e6ScKHnrpsLXyDz+crVat04gRIySMqkDh7h9eXgVVoW5uLy1EJqKff/65RJt+v5jo3Pz333/rPlRZInqXSdAX71nUavXUqVNFUiKTydzc3LRfPqFS0ciRBbtfhocX/3Xp6ekhISHe3t6vvfZaxYoVra2tdV+T8zQxXyNtcRcjoosXLwJwdnbWvYQvIyND3Fd8yf2W6Gj65Rfy8CBr68cZ/Jw5Uu158kxdu3aFdJ3EixBb1+urMau2SkXiHhcX17x5c5Gv//XXX5oNb93d3aOjo+nJPLhGjRovvW1x6tQpsXDh6TnpqKgoPz8/0cFtV7t2VL8+BQTQC2on1qwhCwsCqE+fkvZEPXnypFh5A6BHjx6aLO369eua9NHT01PrmYY9e/aI5dViMWsZmG396aefRLZds2bNzZs3azdIVlbWe++9J27nPbPpQWhoaPtHNQY1atQICAh48YKBsLCwpg0aRNeoQQMGaJe1a7Rq1QrAkSNHsrOzLSwsTExMSjzx/xxDhw4F8Ouvv0oyWvmRmJgoumoWnkMqPpWKatUiuZwOHpR4mwiiJ7t/TJpEBw8Wnbt6oWXLlolvk5+fny5RJCQkmJiYmJiYcOuPIg4ePCiqZSRY2/AklUo1btw4ABYWFgMGDNDcMe7cufOWLVtKdqhXq8nbu6BbkbZFp2q1WkwNnD9/Xq1Wz58/v3fv3sWc9Xgpf///NWs2du5cHZr/sueoVq0agLVr1+oySFJSkkiOK1asePz48WK9Jjub9u6ljz+m5s1p+nTd24++wLfffgtg5MiRko9MRImJiWJ9aqna1tDwibsma2/cuLFYTpqbmxsQECCKHMzMzHx8fMT868GDB1u2bCmOX7169Xpmv47CDViaNm16+vTpZ/7QpKSkpd9/r65du+CKsGJF+vRTerquWnRnA8jXl7T6Z1OpVEFBQaL/gFwuHzly5OzZs0XbOCcnp02bNmkxZmGJiYkiSQXQr1+/MtAs8vTp05qrnd69ez+9m9qL3b17V2wgamtru23bthc8MzQ0VDRkBFCrVq3AwMBnfjNDQ0PFWfOTKVO0+wwUJvafE8XQIok/qm3jvyJWrDjbocN9Hx9pTqXlytKlS8UlnBZfn127Cu4JS525EZEEewCvWbNGVDtol7tfu3bNz8+vcuXKNWrUqFmzZqk6dZUGKpVKJEbSVo/k5+ePHj1a3GQWuyM/ePDA399fVOYAqFevXkBAwAt2KtVQq9Xb587VMWsXRHO2n376iYjEiVjrPrZFnDxZcD+bSUvMSoiPTePGjf39/bXYni82NlZ0LatWrZqW+9UUaT8qtaioKFEpKuGyscJEwba0e/XoyMCJe1zcA5G1N2/evMjE8/379729vcWMUbVq1US1SV5e3s8//yyWcpqaml6/fr3wS27cuCGuC0XL85fPZeblUXAwdelSkJ2bmtLw4ZSS8vi2zty5ZGJCOu87nZSUNG3aNLHJoomJiUwmGz58uIQ7uW7cuLFy5coDBgyQakDDElc7YiMbS0vL4tdHHTt2zMnJCUD9+vWLUzSiVqtDQkJEAg2gTp06gYGBhe/P7Ny5U6y8GT9+vCQ3NJYvX45HXUrHjx8P6VY0njhBADVrJslg5YtKpRJ3YLRYBfj22wTQ11/rIy56ovuHtrPdf/31lzjyTJkypZgTw5mZmb///nv37t3FTU4x7wtg9uzZ2sVQhn344Ye639MoLD8/f+TIkQCsra2L9KUW+zaIBU4AqlatqlQqX5yKiU3g57i56T7T+euvvwJ45513iMjHxwfS7ZqUn08ODgSUaP869hJqtXrgwIEAqlevrumRYGNj4+3tXfxVGTdv3mzQoIHI+2+LSQQt6DwB8VItWrQAIK5yJSc63Um4xld3hkzc4+OpTRt1x45faebanxYREdG5c2fxmevQocOJEyeIKDEx0cfHx9PTU/O0wiXjNWvWLHEf/tOnycuLTE2pWTP66afHt3X27iWJLrMGDRoE4K233vrll19EhxNpp6/u379fuM1qGfDw4UNvb2+ROtSrV++ljX7Xrl0r0os+ffqUaF5BpVIFBwdrTodNmjQJCgpSqVQhISHi3sikSZOkKkOKjIwU1xVEFBgYqEnidZeTQxYWJJNJtsd1uXL06FGZTGZhYXHjxo3ivyo2lkxNycSE7t3TX2hPdf8ouR07doivxoQJE178Sb548aKfn59m81dLS0tPT8/Q0NDDhw+Lfen36eFOt1E7fPiwuGUnSbVMbm6u6AtcoUKF5+0bLQ5NmmI/W1tbHx+fZ7bhEp0WzczMduzYoXts169fB1ClShW1Wi02XdFlZUgRAwcSQBJutsHECnUHB4ebN2/m5OQEBwe7ublpLsVdXFwCAwNffI/x0qVLohbXxcVFp6XDUkxAvJjoUPnBBx/oY3DxN/nhhx/qY3DtGCxxj4+n5s0JoDZt1HFxL2rRqFarg4KCxEyqXC738vISGarmDBQbGyu2jgfg6empxc2gAv/9RydP6um2jihB/uuvv4hIXGDoaRF0GXP48GFxTwaAu7v7My/68/Pz/fz8xHO8vb2Lv8l8YXl5eStXrhRNWwHUrl1bXF+JzXelkp+fL/phP3z48MyZs82ajR0y5EX1PCXSuTMBJN0GweWISqUSd/87dOiwZcuWU6dOxcTEvDQVi48nPz+SbssBPdq/f7847IwYMeLpKYOUlJTAwEBXV1c8Is7rhY9RSqVSzN5JtflA2aBSqZo1a1azZk2xB7t2Bx8hJydHzO/Y29uLKaoXCwsLc3d3F/9eZmZmXl5ehW8zin8vU1PTF1cMlohmBfyDBw9kMpmNTYXsbO3fb2E//kgAvfeeJIMxOn78uKmpqUwm27p1KxFpvvKa4jfxsalQoYK3t/cze2qdOnVKPK1Hjx7SLG7ReQLiBSIiIsTRSarVJiqVKjAwsE+fPiqVSjR36tKliyQjS8IwiXt8PLVoQQA1bvzcdp9FpKenK5VKMQNqb2/v7+8v6plElYh4UMcVGAX0c1tn1KhRmrstYntU7n5VTHl5eQEBAaIZmegsVLiULSUlRVy2mZiY/KxzUVNubm5QUFCdOnWsrKyqVq0q4R1wDS+vdd27nw8NTcvLI0tLksno4UNpRv7oIwKId2HSwhdffCEmmDWTzSLpefPNn1xdadgw+ugjCgigDRto8+aCLU23biUi+ucfEts3bdpEOuRsr8Lhw4fFwqEBAwZoVmNHRER4e3uLnF4cRb29vTW7gV66dOmjjz4SRy2VStWzZ08Ab775puRrMY1aWlqaqHQX93sXL16sxaRMdna2OI6VdA/IyMhILy8vsZJBLpe7u7sfPXr0u+++A6BQKP7888+SRvIC7733nomJxcqV24ho8OC7trZqqXbxvnhRNDjVz1qRciYhIUHsXjJz5kzxyNdff92kSRNNjXt2draYgC9yoa5ZNbF//35xwvXw8CjOUorSQEy6SbJ/6pkzZ8RKOQBr164Vq8P/97//6T6yVAyQuGuRtWtcvnxZs0Fp06ZN+/btK3795ptvPq/YpsT0c1tn4sSJAJYtW0aP5i20rxgrl+7duydWRwFo1KiRqGa7fv16kyZNAFSqVOnAgQNS/aysrKxjx47paWZx2jQC6MsviR7NkUu1mcmGDQRQ//7SjFZ+7N+/X9SBjB8/XqFQtGrVqm3btuIWX/fuWzQ9zcR/EyZQ9+6Ul0fjxhERzZhBorjGy+tFvalKifDwcHFl0qdPn59//lmz1l8mk7m6ugYGBop1QVlZWYVvrLds2VK8/Pbt22K7xF9/3WrQ91Hq5OTkBAUFafqhifKV4h/hMzIyxIYejo6O586d0yKA69evT5o0ydLSUpOKyeXyP/74Q4uhXiAo6KG1tXroUCKiqVMJoAULJBu8enUC9NQwsBxRqWjixAO2tnaurq6a+z+inav4ZE6aNEnT4O7SpUt+fn6aPVAdHBx8fHwCAwNFZd3IkSONaD26JKtNMjIylEqlWBRUrVq1sWPHiguYChUqHDlyRKpQdWeAxL17dwKoaVPSevPs0NBQka5VrFjR0tIyICBA+hkgqW/r+Pr64tFKRFFOrY+GuGXerl276tevj0ed+MWWzi1btizmJrWlwdq1BNDAgUREvr6Pk3jd/fcfAVSlCs9alURMzMm33rIwMXn//ffNzMxkMpmmD2l2dvbNm2lhYbR2LX3zDfn40JAh5ONDq1aRv//jxH3ECPL2pvr1jSBxJ6LTp09XrlxZs9mhs7Ozn59fVFSU+NOIiAgfHx+x+h+AnZ2dt7d34Y2rN23a3KNHgKWlWqv0soxTq9WhoaGa8hVTU1NPT8+TJ0+++FXp6elib7iqVavq2LkiPj5eqVRWqVJl5syZ+ugMe+UKAeToSGo1BQcTQG5ukg3u708LF+p5uUg58PXXogL52p07j5c9vLjGPTU1ddmyZaJ1DABx/+2jjz4yrhtr+/btc3BwMDU1dXV19ff312JLk7///ltM28vl8iFDhmgm3d3d3Yu5OeYrY4DE/fRp6ty5xHPtRWRnZ4vLxL1790oUl359+umnAL766it61AdQws13ypXc3Fx/f38LCwtxfHF3dzeu9tLXrhFATk5ERGvWPE7idbdzJ/32G6WnF1RxsJfLy6MePQiIHzbszdatUYwlDb6+9OABjR5NPXsSGduMu/D7778DqFy5ckhIiJhRS05ODgwM1HRH1ZzXn7nrs+gJ3rQpSbQDQRlUuHwFgKura0hIyDPToLS0tNdffx2Ak5OTVLun6akpnlCtGgF05QrFx5NMRlZWJdpa4EV27nyiCI1p4fhxMjUlmey5f4dXr14tvPq8SI378ePHxaJnCZcdvzL5+fn79+8vvNdvw4YNP/roo3379r105cn9+/c19/Nbtmw5fvx4zaS7dvt76JveE/fC30YJ60HT0tJEF4jCzftKs8Rvv81o0iTZ35+Icvr2Vdvb5xdzIwP2LFFRUbt37969e7dxzQoQkVpNlSoRQHfv0o0bNHEi/fWXlkNFRDzR9MjX94kqDvZys2YVlNa+/rqqSpUf33vvpbeGReJ+/z7Z2RE9mbgbSyJ7/PhxAJ06dRK/nTVrlqbEokqVKjNmzHjx5glZWdSyJQHGsSrXgGJiYpRKpdjsD0CLFi0CAwMLd7ZNSkoSNQw1a9Ys0tq41Bo6lABatoyIqGlTatKEStKH6bF9+2jECCrc4ogPXzpKSKBatQh4eU+NzMzM3377TbNfikwm6927t8huRQcVfSzuejU0e/2KWkfB2tra3d09MDAw5jkTxuJyxcbGZsKECZpJd29v71LbQUTviXvhb6OEs1OnTp0qXHxpBBYtIoDEYpHXXyeAStqzkpUVffoQQNruDEtElJVFSiWZmlKrVo+vgX19n6jiYC+hVtPIkWRiQuPGFdQY3b1r6JhehfDwcADt2rUTv124cKFcLndzcwsKCirmQrQLF8jSkgBat06fgZYJqampAQEBYl0THjVfT0hISExM7NChA4BatWqJDcKNwi+/EEAjRhARabdqMTubPv6Y5HIC6LffCh5UqfjwpRO1mjw8CKBOnUowK3rp0iUfHx8bG5u33npLPLJgwQKUiR0bVCpVRESEUql0cXHRFAjJ5XIXFxelUhkREVF4ym/fvn19+/YVnVgBtG7d+qUVbob1KhJ3zbdRwnrQoKAgAMOGDZMoTP374QcCSLQCfeMNAuhljclZWfXZZwTQrFlavvzQIWrQgABSKOijj0hM4f35J02e/EQVByuWvXvJ2prkcsnWCJd6YjOB1q1bi98mJCTcunWrpIP8/DMBZG/Pm+YUS05OzurVq8U2MWJuT8wINmzY8Jkt2Est0f7F2Vnrl+e0alWw2+GCBZSfT3l5pFTS228XLUJjJSLyi8qVn7H/+0ulpKRoFonNnTsX0m2tVUrcunUrMDDQ09NT0z5LXEJ7eXkFBwenpqZqNny0srLy9/cv/XUcryJx13wbJZxxFyXj84yo9d3y5QU9KYho8GACaNMmQ8fEDOPsWVq9WptbzKmpqZMnT+7WLQqg5s1J9HqOiaFBgwigli2fqOJgz6VWP94dmYjOnNF9d2Qjcv78eQDNmzfXcZxBg8jBgYxkkVFpIZqvy2Sy7t27N2rU6J6xLcZUq6lWLerVi0q68bfYJPG111pXraqqU4dEi47Ll6ltWwLIxIRGjODDl/YSE2nwYAkmA2fNmqVZjFf2pKenb926dcKECZr+rQBEk3EAHh4extLr73Ehv175+6NRI7Rtix07kJUFAI/uXWjp8uXLAERvGeMg2jiINy8qSrOzDRkPM5yYGLRqhTp1sG0bBg4s7qt27do1adKkO3fuVKmyf8GCCzNnmpiaYulS+PkhLQ0ODvjoI4iNNVJS9Bd7mbBkCRISMH06Zs/GyJHo1Qtt2hg6pldHrN/Kz8/XZZBduzBtGgICcPYsAOzejY4dYW+PzZvh4QFTUykCLYu6du3atWvXc+fO1ahRw8TERHTWNyIyGZYsQbVqqFChBIev+/fvjx07dvfu3QA++2zjZ58NtbbG779jyhRkZKB2bfz+O7p1K3gyH76KY9cuODujdWts2wYLC3TsiM2bsXkz8vJ0+vbl5eUBMC2jX2Bra+uBAwcOHDgQwMWLF3fs2LF9+3aVSmVqaurj4/POO+8YOsDi0nviHhBQ8IuUFFy5gubNoVDg4kVYWek0bEZGLUvLio0bN9Y5wFelcLJeOIln5c/u3YiMxL59CAkpOF096r/3XFOnTl2yZAmA9u3br1y5skULkxs3MGEC9u8HAHd3LFuG6tX1HXhZsXcv1q6FlRV8fbF+PXr1MnRAr5QkibvmM/zFFxg4EKGhaNgQ9vbYuhX9+nHi/hKisZiR2rPn8eGrOIl7VFRU586dHz58WKVKlRUrVgwcODAuLm/oUNOdOwFg3DgEBMDWVt9RlzWFTyIODpJ9+8p24l5Ys2bNmjVr5ufnl5uba2ZmZuhwSkb+Kn9Y48YYPRq5uZg1S6dxcnNx8OAPeXkJDRs2lyg0/XNxQWAgxo8HAF9fhIbi0U5SrBwaMwaLFwPAkiVwdES7dpg3D5cuPff5zZo1s7S09Pf3P378eJMmTb755psxY+7u3w8nJ2zYgO3bOWsviRo18PAhAERHo7nxHEMkYmpau0aN3IoVL+g4juYzLHz+OSZOxPHjOo7KjIDmn/699zBxIv76Cw8ePPfJ9erVa9++fZ8+fc6ePTtw4MDdu3d37Nji+vVse3usW4eVKzlr11LhL6BU3z5Ly/lNmqTb2b2n60DGw+iydrzixB3AggWwtsbGjThyRPtBrl9Hfj5q14aFxauOX3t16mDCBNy8iU8+gYkJ3Nzw2muGjokZjIcHLl/GzZt4+BAKBU6fxvz5aNYMTZpg5kwcPx6jUqkKP3/SpEmXL1/28/O7cOFCp06d/Pz8rl59a8KE/IsXYTz390qNiRMxcyYmTcKBA/DwMHQ0r5pCYfLff6YxMeY6jqP5DAsLFyIwEJ076xoeK/3EP31UFDZswPLlGD4cVauidWvMmIGdO5Ge/kQVqEwm++uvv3bv3l2hQoVJkyb179//9u2rHTt+d/EiRoww1DsoCwp/AaX69j18WOHyZev8/Mq6h8f051UnvtWq4eOPAWDGDBBpOciVKwBgRPXtBZYsQXw8pk/Ht98WlDiwcszfH6dPIyAASUkICYG3N6pWxZUr+N//MHHiUEdHx6FDh/7+++8pKSkA5HJ5tWrVFi1a1KFDh9OnT1erVm358i+WLzd5aY0Ne4ZmzbB+Pb79Fj/9BGMrMtad2KIkL0/XcWJiMGMGTp8u+riO65eYUfD3x9mzOH4cAQHw9ISNDc6dw3ffYdiwtCpVHNq1a/fpp5/u3bs3JycHgK2t7enTp8WuXhYWFv7+/kFBnxZaH8i0JE4iAHJz4euLBw90/faJw4IRzkGXLzLSOn3WVno6GjbE/ftYvx7DhmkzwpdfYs4cfPIJvvlG6uD0avDggsrac+ewfj2+/trQAbHSRaXCkSPYtSt327ZWV8TlKWBubt6jR49WrVrt2LHj8uXLcrl8ypQpX3/9deHOVowVX0ICKldGpUoF5ULayc1Fjx64fRt//12uVvayZ8vOxvHj2LcPN26EbdjQS7OCwtraunv37paWltu3b8/Ly2vZsuXatWubl7/6NH0LCMD06ahZE2fO4NG+qNp47z2sW4e1a/Huu9IFx6RmgFITGxt88QVatTq8YsU74nK8pC5fBoxxxr18V9ayl1Io0KMH/P3NLl++HB0dHRAQ4Obmplar9+zZ88svv1y5cqVevXp79+796aefOGtnWhMz7rqtTcUnn+DECZiY4NHOQqxcs7BAz5748kusW9ctOTk5NDTUz8/PxcUlMzNz165d4eHh+fn5Pj4+p06d4qxdHz74AF274s4deHlBrdZ+HDHjXg7Wpho3A8y4A1CpqHXrlhcuXPj2229nzJghHtyyZUtGRkZGRkZycnJmZnZ6ujI9HcnJyMiAWg2xAl1o2xaRkTh+HJ06vfrYdXDxIhYsgL09TE2xcGE5vEfPtJCQkCBaz8rl8pSUFGtra0NHxIxbRgZsbGBtjfR0LUcICcGgQTAxwcGD6NJF0uBY2RIbG/vrr7/OmTOnXr16UVFRhg6nLLt7F23a4OFDfP01Pv1Uy0EGDcK2bdiyBYMGSRkbk5ZhEncA//zzzxtvvGFvbx8VFaXZsyqrUIdEU1PSVGHKZMjPh1wOAESoVAlJSUhKgr39qw9cZ+np4OlSVkLVq1ePiYm5efNm7dq1DR0LM245ObCwgLn5EztJ+PvD2hqdOqF165fMt0VFoV07pKTgxx/x4Yf6DpYZvaysLFtbW7lcnp6ebowdPIzIvn3o2xdyOfbvf9wX/5ny8nD5MiIjERmJL798nJK89RZ27sTff+PNN19BvExLr2gDpqf179+/b9++e/bs+fLLL7///nsAnp6e+fn5NjY29vb21tbWtrb51tYmdnawtoa1dcGSC7HpwIMH+PVX48zaAc7amRbq168fExMTHR3NiTvTkZkZgoOfWMSmVsPfv2DjGwsLtGuHjh0/bd++TefOnWvWrFn4tVlZWWPGmKakmHh6ctbOisXS0rJevXrXrl27evVqixYtDB1OWda7Nz79FF99ha+/vtGokY2jo6PmjzIzM8+fP3/mzJnIyMjY2MV79thq6pQ9PeHqWvDr0aPRvbsR1iGXNwbctfXcuXMKhcLU1PTatWsvfmZGBsXHU3Q0TZ1K3btTXh6NG/dqYmSsVBg7diyAZcuWGToQZvR27qTISCKirVsLHsnJoWXLaOxYatqU5HKqX//xVLyzs/P06dM1rx0/fnyNGj369k1KSTFA5MxIDRkyBMDatWsNHUjZl59P77+/G0CfPn1UKpV48MMPP1QoFJovddeuUXI5NWxIw4aRvz/dvl3w2qePDKx0MtiMO4CWLVuOHj161apVb7/99htvvJGSkpKWlpaRkZGZmZmUlGRrezg62kqUuWvKeQYNKrrrB2PlQb169QBER0cbOhBm9ArvuWhnh9xcdOyIiRMxcSIApKTgzJnkI0e+OHny5MmTJ+/fv5+WliZe+Oeff/76668WFhbbtt20s+NWMqy4mjdvvmnTpvPnzxs6kLJPocDCha127nQODQ398ssv586dC8DR0VEulzdv3rxNmzZt2rRp107WosXjfa/y83H+PO7cQWhoyfbEZYZisBp3ISYmZuDAgdHR0UlJSUX+qFmztIsXC6pKLCxgbY0KFdCyJVaswMcf484d7oTOypH169ePGDHi7bff3rRpk6FjYcZt2jS0aoX4eFy7hps3ceAAADg7o2tXuLrCxQUdOjxu5Hz9+nWZTFa/fv2LFy927NgxIyNj5cqV48aNM2D8zOhs2rTpnXfeeeutt3bs2GHoWMqFAwcO9OnTh4j++eefPn36pKammpubm5sX7LmWnZ1//rzJmTMFNe7//ovsbFhaYvx4tGlTcGRYudKw74C9iCFn3AFUq1YtPDz8hx9+yM7OdnBwsLa2tra2flTmDisr2NnBzg6amzzTpgGAvz8aNTJc0Iy9cvXrN2/bdqapaUdDB8LKAg+PgumP7t2Rm4szZ3D/PjZswIYNANCixXt2drddXFxcXFy6detWp06d9PT0oUOHZmRkvPvuu5y1s5Jq2bJD9+5/yuUtDR1IedGzZ8/Zs2fPnz9/5MiRkZGRtra24eHhp0+fvnTp0sWLF69evZ+QEK2Zs5XJUL8+2rRBbu7jIwMrzQw8484YK47kZDg4wMYGj8oWGNPStGmYPRv5+WjUqGBBal4ezp3DyZM4eRLnzqkvXrRUqXI1z3/ttddMTExu3brVokWLEydOWFlZGSx0ZpxUKtjZISsLiYlG21XC2KjV6n79+u3du9fW1jY9Pb1wpmdiYtKjR2r16pZt2kD8J3pTP31kYKUTJ+6MGYeKFZGUhNhYVK1q6FBYmZaamvrvv/8ePXr0yJEjJ06cePjwYbVq1WrUqPHbb7814X4TTCvt2uH0aRw58riBCdO3+Pj4uXPn7ty5MzY2tkGDBi6PtG3bli+/jRon7owZh/btERGBo0d5yxv26hDR1atXY2NjX3/9dUPHwozY2LFYvRrLlhWsgWavzH///efs7GxiYuC6aCYhuaEDYIwVS716AMB9ZdirJJPJGjduzFk705Fo4M59ZV69GjVqcNZexnDizphx4MSdMWakOHFnTCqcuDNmHDhxZ4wZqebNAU7cGZMC30BhzDi0aYN33uECd8aY8XF2xoIFaNgQKhUKbeLJGCsxnnFnrPTatQtnzwLAtm1o0wbjx2PECADYvBl5eQaNjDHGXqbwEax9e/TtC4WCD1+M6YQTd8ZKr9274euL/HyEhABAaCjEFsNbtyI398UvZYwxAyt8BOPDF2OS4FIZxkq1MWOwePHj337+OWxtcfy44QJijLFiK3wE48MXY7rjGXfGSjUPD1y+jJs3C367cCECA9G5s0FjYoyx4il8BCt8+MrNRXg4eCMZxkqKE3fGSjt/f5w+XfTBvDw+7THGjMDTRzCZDN99h2vXsG2bgWJizGjxzqmMGSV/f9SoAWtrDBpk6FAYY6yE/vsPv/yC6dNRtaqhQ2HMqPCMO2NG6b33cP4818wwxozStGmoUgUPHhg6DsaMDc+4M2aUhgxBly7o169gZxPGGGOMlXmcuDPGGGOMMWYEuFSGMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s4YY4wxxpgR4MSdMcYYY4wxI8CJO2OMMcYYY0aAE3fGGGOMMcaMACfujDHGGGOMGQFO3BljjDHGGDMCnLgzxhhjjDFmBDhxZ4wxxhhjzAhw4s5YUfcPHtzUoMEqR0f89puhY2GMMcYYKyAjIkPHwFhpolLlOzubPHhQ8Nv16zFsmEEDYowxxhgDeMadsaJUKpMHD1RAuPjt3buGDYcxxhhjTODEnbFnUADtAQCxcXEGDoUxxhhjDAAn7owVZWYGf384OGSam58Cvrp3z9ABMcYYY4wBXOPO2PPcvHmzYcOGMpksOjq6Ro0ahg6HMcYYY+Udz7gz9mx16tTx9PTMy8v78ccfDR0LY4wxxhgn7qwc2r0byckAsHkz8vJAhGXL8MknuHSpyBP9/PxkMllgYGCyeD5jjDHGmOFw4s7Kn9BQJCUBwNatyM3FkiWIj8f06fj2W+zfX/iJrVq1cnNzS0tLW7ZsmWFCZYwxxhh7hGvcWfnz8ceIiYGtLfbvx9mzGDkSa9fCygrnzuH77+HkhNat0aoVGjaEQrF3794+ffpUrVr11q1bFhYWhg6d6Ud8PNLTIZejdm1Dh8IYY4w9F8+4s3Jp4UIEBqJzZwCoUQMPHwJAdDSysrBoEUaMQNOmsLNDx46qQ4dq1qwZFxe3YsUKw4bM9OXgQdSpg3r1UKcO5s0zdDSMMcbYc3Hizso3mQwTJ2LmTEyahAMHMHEi5s7FwIGoXRuZmTh16s/Nm+/cuQPAx8enfv3677zzTm5urqGDZpL66y9kZiaLX//2m0FDYYwxxl6ES2UYAwCkp8PG5olHkpNx9uyKCxeOnT79559/5uTkADA3N9+4caO7u7thgmT6MHkyli27DDQBcp2czO7fN3RAjDHG2LPxjLs24uLirl27durUKUMHwqRTJGsHYG+P11+f8MEHv/322+LFiwHY29vn5OQcPXrUAOEx/WnfHkATAMAJtdqwsTDGGGMvwDPuJaNWq9esWTNjxgw7O7sbN254enouWbKkcuXKho6L6Vd2dnatWrXi4+PlcnlUVFSdOnUMHRGT1P79m1es+HP9+r/l8vuJiRUqVDB0QIwxxtgz8Ix7CURERHTs2HH06NEPHz60tra2tLTcsGFD8+bNN27caOjQmH5ZWFjUqlXL2dn53Xff5ay9DOrV663Vq0PMzLLUaqVSqb+fExcXd+nSpfXr1+fl5envpzDGGCurOHEvluTkZF9f306dOkVERFSrVi0oKOjff/+9ePFi79694+LiPD09PTw8YmJiDB0m05ddu3ZFREQkJSXNmjXL0LEwvTA3Nx88eDCAVatW6ek+5N69e9u0adO7d+8RI0Y0b958w4YNfMOTMcZYiXDi/hJE+P13DBuW+uOPP8rl8k8++eTatWujRo1KSUmpU6dOaGjokiVLbG1td+zY0apVq40b43T5WTk5Odu3bx81atSbb74pVfxMd8nJyd7e3kT05ZdfNm3a1NDhMH358ccfZTJZWlraH3/8Ie3Iubm5M2bM6Nu37/3796tWrVq/fv1r164NHTq0S5cuYWFh0v4sxhhjZRmx57twgXr0IIAAGjVq5YULF4jo3r17Xl5eDRs2zM7OFk+7d+/egAEDXF1HAfTmm/TffyX7KXl5ef/888+YMWM0lbUKhSI2Nlbyt8O0M3z4cACurq75+fmGjoXpV7t27QDUq1dPwjGvXLnStm1bACYmJkqlMj8/Py8vb+nSpU5OTuL7PmDAgMuXUyX8iYwxxsoqTtyfLSODlEoyMyOAKlakwEBSqyk3N/ebb76xtrYGYGtre+rUKc3z1Wr16tVZ9vYEkL09rVr18h+hUtHhwzRr1qoqVaporqPatm37zTff3L59W4/vjZXE5s2bAVhbW1+/ft3QsTC9O3LkiPgmRkZGSjLgmjU7xBGjbt26J06cEA+KL3hGRoa/v3+FChWaNRurUJCnJ926JcnPZKwU2bmTxJdp61b65x9KSiIi2rSJcnMNGRVjxosT98c0x5fAQKpRgwCSy2nyZEpMJCI6dIiGDTskTuqenp537959eoT792nw4IIZ+n796Hnp94ULpFRS3boEUPfuZwE0adJEqVRevnxZb2+OaeP+/VjRMmjZsmWGjoW9Is7OzgBq1ar1559/Joovv1ZSUujdd8nZ+b+KFSt5enomiYSF6OzZsyYmJp6enlFRUUQUFxf3+ef3TU0JICsr+vxzSkmR5H0wVir4+lL37pSXR+PG0YwZdOMGEZGXF6WnGzoyxowTt4N8bNo0REZi3z54e+PSJeTmYskSdOqE+Hh88gn++AOmpujSZdKnnw7u16/fC8ZZswbTpiEhARUqYOlSNGmC1q2xbRvq18f69Vi/HlFRBc+sXRtjxmQOGnS9VatWr+IdshIaOJBSUk5aWfn//fcWmUxm6HDYq7Bx48YNGzYEBwcDUCgUnTp18vDw8PDwKNHyhqNH8d57uH0btrb49df/hg6tofmjVatWTZ06NTs729zcfOrUqZ999lmlSpWuX8dnn2HTJhChShXMmYNJk2BqKv27Y+wVmzYNrVohPh7XrsHBATExsLXF/v04exbW1oYOjjEjxIn7Y4WPLwsXokoVKBTYvBnjxyMpCRYWmDULM2fCwuLlQ8XHY+pUVKoEC4uCi4Hx47FtG5KTAcDZGZ6eGD4cnTqBs8FSa+VKjB8PBwecP0/Vq/O/Uzly586dzZs379ix4/Dhw5q+jQ0bNhwzZkPnzi27doWJyXNfq1Lh228xZw7y8tCuHdatQ4MGRZ9z9+7dBQsWrFy5ElA0b548aJClnx8sLREejpkzcfAgatXCjz+iZs2Ca/6BA/X2VhnTs2nTMHs2Pv4Yd+6gbVtMnYo6dTBqFH74AVFRaNeOT4KMlQwn7o8VPr7s31/wYEQEOnbE66/jl1/QuHHJBszLwyefPL4YsLdHWhpGjED37lAoJA+fSenWLbRqhdRUrFuHESMMHQ0zkIyMjP379+/YsSMkJCQ2NrZu3ZwbN8wcHODmBnd3eHjAweGJ59+7hxEjEBYGuRyffIIFC140a37u3LmlS28HBg4A8NprmD8fo0dDocD27ZDJsHdvwTX/xIlYuVLP75MxvREn1vx8NGiALl0wcyZ698aoUWjQAHXrwtoagwYZOkTGjAon7o9pji+NGiEl5fHjkZFo00anMYtcDLBSTq1G7944eBCDB2PzZkNHw0oBlUp1/Hj49u2dduzApUsFD5qYoGtXDB6MBg3g7IzWrfHHH/Dzg1qN1avRv3+xRt63D35+OH0aAH78ER9+WPB44RuAnLizMiAgANOnY8IELF8OAP/9h19+wfTpqFrV0JExZlS4j/tjAQGwtsb69UVvTGudtWv4+xecmJlR+OEHHDwIR0csW2boUFjpoFAounbttGgRLl7EjRsIDIS7O+RyHDyIbduwezd8fZGfj4MHERKCCxeKm7UD6N0b4eFYtw5dumDs2Cf+yMMDly/j5k3J3w1jBhARAQDt2hX8dto0VKmCBw8MGBFjRoln3J+gUsHJCQ8f4tIlNGmi62g5OfDwQP/+mD6dy/iMxtWraNMGWVnYvh3u7oaOhpViKSnYswe2tvjnH+lnx/v1w7RpaNkSTZs+cQOQMSPVqBGuXcOZMxLMhTFWnvGM+xMUCnh4AMCWLRKMFhqK0FCsWcNZuzGpVQs+PvD25qydvUSFCvD0LJhcl3Z2PCsLhw/D3R3m5py1s7IgJQVRUbCwQPPmhg6FMSPHiXtRgwcDEiXumzYBwJAhEgzF9GfXLpw9CwDbtgHAoUP49FMsW4bNm/GooQhjL+fvj1OnpBnq8GFkZ8PFBZUrSzMgY4Z1+jTUarRqxU1Oy6DC59Dduwu65/EJVH84cS+qb1/Y2iIiQtfJs/x87NgBAG+/LUlcTF80BcohIQAQGoqkJMhk2LoVubmGDo4Zg4AAVK6M4cORkSFN7h4aCgB9+kgwFGOlQXg4ALRvb+g4mB4UPoeKEyjAJ1A94sS9KHNzvPEG8Gj+VWsHDuDhQzRtKkGtPNO3MWOwePHj337+OSZOxPHjhguIGSGx6m71agmG2rMH4MSdlSExMX/16BHepUuioQNhelH4HCrNCVSlAgAiLFuGTz553MyLceL+TJJUy3CdjBEpUqC8cCECA9G5s0FjYsZG9IRZtw6ZmTqNExeHCxdgbc2fQFZ2hIR8euhQh1atYg0dCNOLwudQXU+gFy9i1CiMGgUAS5YgPh7Tp+Pbb7mjtgYn7s/w5pvo3Xt+TEyz+Ph47UZQq2n7djUeXQOw0k+07FSrcf8+Ll4seJBXFbPia9YMHTogJQVbt+o0TmgoiNCjB8zNpQmMMcN6+PDhrVu3bGxsGjVqZOhYmMQyMpCTAzyr7XV6Oo4dK8lY58/jnXfQogX++APbtuHhQ+zdi48/RrVq8PUtqCBknLg/k50dLCzCo6IubdO2XObYsaMqVY233/6K+16VfqJA2ckJKSkIDMS6dfjjDwD4/XdYWRk6OGZUxKT7b7/pNMiFC3906BDyxhu6zdszVmqEh4cDaNu2rYL3DC9boqLQuTOio+HgUHAO/fZb1KkDAL//jsWL4eoKDw/cuvWSce5fuoShQ9G6NTZtgoUFfH1x/ToqV0aNGnj4EACio7khkQYn7s82ePBgAFsKlcusWrVqw4YNu3fvDgsLO3PmzI0bN+Lj4zMyMp758o0bN8bFxTRokPqKwmUSefNNyGTYubNgCoGxEhkxAlZW2LcPN25oP8jatZ+dOjWwZ0/eeImVESJxb89LU8uWHTvQvj3On8fduwXZdRF2drC0xI4daNYM8+Y9u4bw0qVLo0aNqtO6dXZkJExM4O2N69cREABTU8yahd69MXMmJk3CgQMFvboZAGLP8vDhQxMTE1NT08TERCLKzs5+wd+hvb199erVGzZsePz4cSJSq9U1a9YEcOrUKUO/D1ZirVoRQP/8Y+g4mHEaMyazS5e13323VruXnz9/HkC1atXUarW0gTFmKB4eHgD+/PNPQwfCpKFWk78/yeUE0IABlJz83Gfevk3DhpFMRgC1b5+1bt06zZHt4sWLw4cPl8vlACwsLHbMm0f37hERPXxISiVVqEAADR1KRJSW9grelBExeSVXB8anUqVK3bp1O3DgwM6dO9977z2VSjV27Nj09PSUlJSMjIyMjIzU1NS0tLSMjIzMzMzk5OTk5GQAKpUKQHh4+J07d1577bV2ms2dmfEYMADnzmH7dvTrZ+hQmBEaOfKYm9t7MTG1fX0LzkmFZWRk3Lt3Ly4u7t69e7GxsTExMeL/Y8eOfe+99wCEhoYC6NOnj4wXWDDjl5SUFBYWdvjwYQAuLi6GDocBQFpaWlZWlqOjo3YvT0pKGjfObutWhUKBRYvwyScvWgxWsybWr8dHH2HaNBD9/O67n3z//fcBAQF///23v7+/Wq02NzefMGHCp59+Wr16dSQmYt48BAQUbDvn5oYZMwDAxka7UMssQ185lF4//vgjgCFDhrz4aWq1OjEx8c6dO1euXMnIyCCimTNnApg2bdorCZNJLDycAKpenXjGk2lBrVbXrVsXwP79+8Ujn3/+effu3Rs3bmzz/NOPn5+fePIbb7wBYM2aNYZ7B4zpJDU1NTQ01M/Pz8XFRVy7mpmZmZqauru7JyUlGTq68m7u3LmTJ08G4Ozs7O7u7u/vHxYWlp2dXcyXnzt3rl69eq6uSypVoj17SvBzVSpauXKlk5MTALlc3r17d1NTU29v7//++4+IUlNT/f39V7i5EUAAubnRiRNavLtyQkZEr+QCwfjcvXu3Zs2aJiYmHTp0sLGxsbOzs7Ozs7Kysra2tre3t7Gxsba2Fr8Wv7C1ta1QoYKVlVXLli2joqLCwsK6du1q6DfBSowINWvi7l2cPo22bQ0dDTNC8+bNmz9/vpeX1++//w7A3d3977//Fn9kaWnp/Ei1atWcnJzE/5s0aVKzZs3c3NxKlSqJWXlnZ2eDvgnGSuDhQxw+jEOHcOgQ7typm5RUsELD3Ny8Y8eONWvW3LFjR3JycsOGDTdt2tScVxkayOrVq8eOHatQKCwsLAqv0LO0tGzXrl3Hjh27dOnSsWPHatWqPfPla9eu9fb2zszMbNeu3ebNB2rUKPFEeFpa2sKFCwMCAnJycnr06HHw4MGkpKTvv//+hx9+SE1NNZHLE0eMsP3gA3TqpP2bLAc4cX+uS5cutWzZ0srKKi0trUQvNDExsbW1ffjw4dM3yplRmDkz5tSpPz08qsyYMcrQsTDjc+vWrXr16pmbm8fExNjb2588eTIjI0Pk6Pb29k8/Pz8/Pz4+PiYmJjQ09LPPPmvZsuW5c+deedSMFcuuXXB2RuvW2LYNFhbYsQMHD+LiRWhSCVfXRaam/7z++us9evTo1KmThYUFgDt37rzzzjvh4eEWFhZLliwZKxowsVfoyJEjbm5uOTk5S5cunTBhwpUrV06fPn306NEjR45cuXJFrVZrnuns7Ozi4tK1a1dXV9d27dpZWFjk5+fPnj170aJFALy8vAIDAy0tLbWOJDo6unv37llZWe7u7iEhISkpKQD69u07b968zrx7RTFw4v5cgwYN2rZt2/DhwydPnpyZmZmWlvZ0gXtGRkZycnJ6errm1xkZGZaWlqmpqWFhYV26dDH0m2Da2LVr15tvvtmmTZszZ84YOhZmlHr37r1///7ly5dPmDABQG5u7sOHD+/fvx8TE5OUlCR+ofn/nTt38vPzxQtr1apVr169ffv2GTR8xp5r2jRERmLfPkycCLW6YKtgKyu0aYOuXeHmhq5dYWEBALm5ueHh4ZmZmX369AGQlZU1ZcqU1atXA/jkk9kLFy4wNTXc2yhnbt261bFjx/j4+OnTp3/33XdF/jQpKenEIydPnhSZtGBubt62bduEhIRr166Zm5v/+OOP3t7eusfj5OQUFxfXqFGjq1evurq6LliwoGfPnroPW05w4v5sp06d6tSpk5WVVVRUlKjKKr7PP//8q6++qlu37rlz515Q1cpKrZycnCpVqqSlpd28ebN27dqGDocZnzVr1nh5ednb27/22muxsbEPn9ks7RG5XO7o6Ojk5FShQoUjR44AOHXqVFuu02Kl0rRpaNUK8fG4dg3DhyM8HD16oEMHiCw8Jwfh4VkHDnx76NCh48ePZ2ZmtmrV6uzZs5qX//7771OmTGnaNJjozY0bUauWod5HOZKWlubq6nr+/Pl+/fr9/fff//33X35+fv369Z/3/Bs3bhw5ckTMx0dGRqrV6jp16jx8+HDPnj2dpChiuX79esOGDatWrbp+/XpTU1NXV1fdxyxfDFtiX2qJi785c+Zo8dqcnJwWLVqA16cas3feeQfATz/9ZOhAmFHavXu3lZWVmZmZ5kjr4ODQtGlTNzc3Ly8vPz+/gICA4ODgsLCw6Ojo3NxczQs//vhjAO3bt1epVAaMn7Hn8fWlBw9o9Gjq2fPxg0eO0Lx51LMnWVqSiQnZ2lYHIJPJmjdv7uPjU+TDfPr0/Zo1CaCqVenAgVccfrmTn5/v7u4OoEmTJsnJySkpKc2aNatUqdLp06eL8/KUlJQVK1YAsLe3L/4a1hf77bffALz99tuSjFYOceL+DNu3bwdQuXLllJQU7UY4c+aMqampXC4/ePCgtLGxV0MsK+zTp4+hA2HGJy8vr2XLlgDGjRt37ty52NjY4jdlT01NrV69OoCVK1fqNUjGtDNhAl24QPfvk53d4wc17UDkcmrZkpTK9Zs2bXrw4MHzBnn4kPr1I4AUClIquYWXHn344YcAKlWqFBUVpVKpRE/9Jk2alKjDT5s2bQBs2rRJkpBEAeG3334ryWjlECfuRalUqtatWwP4/vvvdRlnzpw5AOrUqZPGewcYoYSEBLEDF/cvYyX1/fffA6hXr15WVpYWL1+3bp040T58+FDy2JjB5ebmrl69ulu3bkuXLj169KihwymxuXMJoPnzn3jwt99o2jTaupUSEoo7Tn4+ffZZwdY8Q4ZQZqbkkTL69ddfAZiZmYk5RF9fX3FsuX79ejFHSE5OJqLFixcDGDx4sCRRNW3aFIDYsJJpgRP3ooKCggDUqlVLx7tCOTk5LVu27NTpfzNn8gHJKIkNuj08PNasWXPu3LnC9QyMPU9sbKxoHbNjxw6tB+nVqxeAqVOnShgYM7j09PSAgACxr7aoJLG2tg4NDTV0XCUjZso3bJBmtB07yMGB+vShHTsoMpKIaOtWaUZmhw4dEtV6K1asIKKVK1cCMDU1PVC8+iS1Wv3222+bm5v/999/sbGxJiYmZmZmus8mJCYmyuVyc3NzqQpvyiFO3J+Qk5MjNk8JCgrSfbSzZ1NMTUkuJ66XMTp79+61trYu3EvbzMysVatWXl5ev/yyfc8eio01dIisVPLy8gIwYMAAXQa5ePGiqLU7efKkVIExA0pJSQkICND0Oahfv/533303ZswYcWCRqgLhFVCrqVIlAuj2bcnGvH6dEhLI15e6d6e8PBo3TrKRy7MrV644ODgA+Oyzz6hQEr98+fLiDyLWen3zzTdE1K9fPwBLly7VMbAdO3YA6Nq1q47jlGecuD9B3ONu3rx5fn6+JAMqlQRQ7dqUmirJeOxVWLVqlampKQA3N7c5c+YMGTKkfv36mq783bsHi2pOR0fq04dmzKCgILp8mXbu5Bmj8u748RsKhYmFhUV0dLSOQ/Eq1bLh/v37SqWyQoUK4ujRtm3boKAgcX5Rq9XTp08HoFAoVq9ebehIi+X6dQLIyUn6kX19adUq8vfnxF0CCQkJDRs2FMUtKpUqOjq6cuXKAD799NMSjbNt2zYAzZo1I6I//vgDgKurq46xffbZZyi0VzTTAifuj6WlpTk6OgLYvn27VGPm5VG7dgTQlClSDcn0SK1WK5VKcYot0gwhPT39xIkTgYGBSuXVbt2oQoWCxVjiv/fe4xmj8i4/n1q3pubN0775Zr3uo/EqVWMXFRXl4+Mjth8SGU9ISMjTT/P39wcgl8tLNBX6YsVfDF1Sa9cSQLrdT3q2ZzarYVrIzc3t3bs3gDZt2qSnpxdJ4ks6lMj4z507l56ebmNjI5PJoqKidAmvR48eALZt26bLIOUcJ+6PieWkkt/BuXiRzM1JJqPdu6UdmEksJydn5MiRAExMTAIDA4no3r17LzgF3rxJISH05Zfk6UmrVvGMUXn3008EUM2alJEhzYB//vknr1I1RmfOnPHy8lIoFCIjd3d3P3Xq1AueL3J3mUyme5+NuLg4pVLZtm3bDz/88OzZszqO9rRp0wigBQskH5imTKEHD4o2q2FamDRpEgBnZ+c7d+5okvjWrVunp6drMdqUKVMAfPzxx/SoDnB+kYXJJZGbm2tlZSWTyeLj47UehHHiXiA+Pt7W1hbAkSNHJB98wQICqF49ysuTfGwmjcRE6tVL3bHjVzY2NmJZ4fHjx6tUqTJ79uxijsAzRuXZw4cFtb/SFkqJVapT+IadkQgLCxM9s0XxupeX15UrV4rzwqVLl4piPK1LCG7dujV16lTNBL/oaHzu3DntRnueMWPONm2avnu39OVbbdtS9ep06ZLkA5cv58+fNzU1tbCwENeKH3zwgUji//vvP+0GPHbsGIBq1arl5+fv3r1brNDQ+pbOqVOnADRq1Ei7lzOBE/cCU6dOBTBw4EB9DJ6XR+++S7zMrNS6cYOaNCGAGjRQnT17gYj++usvcQr08PAo5oIHkbhHR5O1tZ7DZaXPuHEEkOR9/3mVqlFQqVQhISEdOnQQGbOtra2Pj8/du3dLNMiaNWtMTEwAfPjhhyVKjKKjo318fMzNzcW0vbu7+5EjRwYOHCi2/YqIiCjhu3muvLw8MV0q1S2g8ePpgw8oNpZycsjMjORy4ubJOlq/fr2dnZ0mk4mMjGzQoMGJEyd0GVNU2oSGhqpUKlG/p/WAYhnhOL4rrRtO3ImIbty4YWZmplAoLly4oL+fwosXS6eTJ6lqVQKoefOCVgkBAatkMpnox1eiZcqpqWRnRxYWxM3fy5XwcJLLycyMije7WjK8SvXl1GpaupQ+/pguXjTIz9+wYYM4Yjg5Ofn7+2u9c9+2bdtE/j1x4sTi/HP/+++/hWtyPD09Lz76G8jJyRG5e5UqVaSadz9z5gyABg0aSDJaXh5ZWpJMRsnJFB5OADVtKsnA5drRo0cB1KhRQ/P50b2R8bx58wCMGjWKiGbMmAHggw8+0G4oT09PAL/++quOIZVzZSdxL5wW//NPQea0aRM9/aHNyMi7f//+1atXw8PDQ0NDN27c2LVrVwBjx44VT3jw4EGeHopaePFiKbRlC1lZEUBubpScTPn5NHUq1auXbW9fW6lUajFg794EEK8nLFfGjiWAZs3Sy+AlXaVa/CNhWXD/Pl2+TD//TPPn0717NHYsbdlCV6++yhAuXLgAwNHRcdmyZdptuVXYzp07LS0tAYwYMeIFKdfRo0fd3d3F1YKoybn61LvOyckZMGCAyN3//fdfHQMjosDAQADvvvuu7kMR0dmz4ianGJkA8vKSZOByTa1W165dG0BYWJhUY0ZFRYk9B9LS0iIjI8XCm5ycHC2Geu211wBcvnxZqtjKJxOUFbt3IzIS+/YhJAQODmjYEPb22LoVVavi88+RmoqUFCQnIzUV9erlXr3qXPi1dnZ2VlZWYhHG9u3b33///enTp8+aNUvyIMeMweLFko/KtJSQgNGjkZmJ8eOxdCmysuDhgV27YGlpvn79uQED7Eo02sOHqFwZ772Hffuwdi3GjdNT1ABARJmZmdbW1nr8Gexldu2CszNat4aHB7p2xbBhEo9PRDKZzNbW9ttvvx0xYsRHH32UkpJS+B/dyqqjWt3K0vLxS2rVeu6RsF8/mJpKHKEWiOjKlStHjx51dHSUy+Wvv/66jY1NyYZIT8eJE9i7F3v34swZvPkmTE2xdi2srODri5kzsWcPatZE796qPn0e9Oyp6Z6uJ5s2bQIwcODAiRMn6j7aG2+88c8//3h4ePz555/p6enBwcGFy9YBHDlyZNGiRaIZtrW19fvvv//JJ5+IfKgIMzOz4ODgIUOG7Nq5U+3nh0WL0KKFLrGFh4cDEDvT6e7MGQBwcQGA06cBoG1bSQYu12Qy2bBhwxYtWrR+/XoxI6m7evXqderU6fjx4yEhIe+++26LFi3Onz+/e/duDw+PEo1z69atu3fvOjg4NGrUSJLAyi8DXzhIp3BPjxkzaMQI8vam+vVp374n2vYB1KBBrqOjY/369V1cXHr37j148OA6dergUYH7vn37ZDKZubm5jmUzajUVuT/JixdLicJTkjt20FdfERHFx1Pr1gRQlSqkxU7Mw4aRQkGXLlFKCllaklxO2q4FegnR+/mDDz6ws7Pz8fG5deuWXn4MKwa93kNLSUlp2bLl8uXLxS3vmjVriom0wnr0CClycBs9+rlHQq1aSkgjLy8vIiIiICDA09OzSpUqInhREa5QKFxcXPz8/EJDQ18wwZyenr5r164z/v7Urh0pFI/fsI0NDR1KH35YUOW2aRONHEmOjuJPTzZrBqB58+a+vr4hISF62vy4devWAHbu3CnhmOHh4ZUqVQLQs2fPtLQ0IlKr1SEhIZ06dRJ/e6KM/v79+y8dKjs7+9/JkwsObefP6xJVy5YtARw9elSXQTR8fVUALVpERAVNkw8dkmTg8k5UNFWpUkXCwoElS5YA6N+/Pz1qgjR06NDivDA7OzsxMfHWrVtXr1798ssvAbi7u0sVVbmlTeJ+8+bNiIiIvXv3JiYmSh6Q1gqnxTNm0I0bREReXhQbS/v3U0QERUfTgwf0zNs79+7dE3uMrVu3jogmTJgAoGPHjrpswxQQQAoFie5e6enk6UleXtzuqlR4Zr6Vm0v9+lH9+lreZh8/ngASHWjeeYcA0rmx2zPk5+e///77AMQiNvGLYcOG6bj2iGlHrw1ARS1ply5d1Gr14cOHxWzCu+++613InDlho0aRp+fj/3755blHwlecuCclJf3999+fffZZt27dLAvfFACqV68+fPjwyZMnd+nSRfNJBlCxYkVPT8/ly5ffvHmTHqX7/v7+bm5uovJ7mri2NjEhFxfy86PQUBK7pl+4QMOG0cSJ9MEHlJJSMGuyePGfkydrblBYWFjoXsfytJs3b4q7IpIPfvHixWrVqgHo1q3b+vXrXcTUNFC5cmWlUplUomU0OTk0YEBB7q5tzcyFCxcUCoWJiUlmZiYRqdXqyZMn67KNjqtr11q1eh84cCk3N7dTpwcVK6q1XRrAimratCmAPXv2SDXgw4cPP/roo8jISCK6e/euQqEwNTWdNm2ar6+vt7f3sGHDBg0a5Obm1rFjx9atW9etW9fJyUnkVEVUrVpVx12li6lsVwyWIHG/e/duQECAq6urTCZzc3MT/wzOzs7u7u5KpTI4OPjChQv62/fhpXx96dIlun2b7OyeOF0Vs6fyr7/+CqBSpUqxsbEpKSk1atQAsGTJBu2CuXiRLC0JILHJgLc3AdS5s3aDMYk9L99KSaGEBC3HPHCAAKpVi9Rq2rKFAGrTRvdIn5Cfnz9q1CgAVlZWu3fvPnPmjLe3t+Y2uouLS1BQkD7WZrDn0d89tGvXrpmbm8vl8lOnTuXn57dq1QrAF198UfyoxASBOBKmp5fgSKiLe/fuBQcH+/j4uLi4aHYaFurWrevl5RUYGFjkTmZycvLmzZsnT55cv379ws+vU6eOlZWV5rcKhaJDhw4LZs+mPXue+06e1ZEkJyfn0KFDc+fO9fX11cM7pu+++w7A8OHD9TH4lStXxJlIeO2113744YcM7f4hc3LIw6Ok8+75+flhYWF+fn7issHJycna2nrfvn1EdOrUKbEoVrvcPT8/39raWiaTJSQkiBlibhEoIbGNoJ6atyQkJFSsWLHwJ/MFzMzMHBwcatas2aBBg2bNmsnlcoVC8c8//+gjsMIKT88ZcP5CT16euMfFxf3yyy/du3fXHIhtbGzeeuutLl26PF2bWKVKlT59+vj5+a1fv/7KlSuvuA2Cmxs1aqT9zcD+/fvjUcHMP//scXX9w9JSmzYR2dnUqhUBNHEiEdHWrQSQubnWMx1MYvrIt1QqqlGDADp2jLKzqXfvre3bu1+SrilxTk7O22+/Lb59+/fv1zwuNlQXm9uJC2mlUpmg9fUHK4nCKTIRaaZcIyJ0XagqjkWTJ08mop9++kkkvtrN6W7bRo6O+t39LSYmZu7cuUXKrC0sLLp16zZr1qwdO3YUvjd7586d69evPz3Iv//+W7FiRZGyKxQKkbXXrVvX29s7ODi4VN3dLaxbt24A/vrrLz2NLzbhql69+ooVK7RbDviYJnd3dHzJafL+/ch164YMGWJn93idj729fYMGDQBYW1uHhoYSUUhIiLgTMmnSpJKe68+fPy/+fenRrJlUa14ZEV25ckX8k2WLW1LSycvLE5tL1KtX76uvvvr+++8DAwP//PPPzZs3h4aGHjt27PTp09HR0ffu3UtMTHx6Jnf+/PkA7OzsJDw5PlMprBiU0HMT94yMjF9//dXNzU1cVQOwtLQcMmRIcHCwuFNGRCqV6vLly3/++efMmTPd3Nw02YOGjY2Nq6vrBx98sHLlyvO6Vde91O+/E0CVKpHWG3LdunVL7MG0bVsYEY0ZQwB16kQlrZdZuHBNxYoPGzSg9HSKiytoNfjDD1pGxSQn8q2LF8nCgu7dk2xYP7/8du12KZVBRDR+/HgAxd+86cUyM/NFJlexYsVnbsGYnZ0dFBTUrFkzzffO29tb30dGppGVRd7e1KsX5edTaipVrEgAad3xbNMmddeuqxs2bPrw4cOEhARR66z1DuFffUUAVatGetp99fLly1WrVhUzsnZ2dm5ubkqlMjQ0tPBlRnR0dFBQkLe3t7iD7/VU95BLly6JJaS9evVq0qQJgKZNm5b+1hNxcXEKhcLc3Dw1NVVPP0K0TPj888+lGS4ri/r3J4B++eUZf3rhAvn7k6sryeWnGzcWB5O6dev6+PiEhobm5OSo1Wqx4Ym5ufmWLVuI6O+//xY3/SZMmFCi3H316tUAPD09iWjy5MkAdN81lhUmll5slbr5tPgAODk5abejk1qtHjp0KIBGjRolJydLG1thz6sYTE2l1q3J25v02fpb756buKekpIjqQ9Fwql27drGxsUWeExMTExMTU/iRe/fuhYSE+Pv7e3l5NW3aVLxW6Nq167vvvqunf6qEhIL1SEFBOo2zYsXvnTv/5ehIDx5QcjK99hoB9P33JRhh//79crm8evWap05lqNXk7l6wLYvhaojYsw0dSgB9841kA4o+WVWqVMnNzT1w4ACAWrVq6V48lp5OvXvT668vcHR0fOke5mLjRvG9k8vlbm5uISEhBixgKyfu3qUqVQggUcyybh0BZGFBZ86UeKjMTKpdmwBaulRNRN7e3ni0Jkw7KhX16EEADR6s9RjPdf78ebHYtFevXv/++6/mk5aZmXno0KEvv/zyzTffrFChQuHZnIoVK04U9yIfOXv2rJj00WTtrVu3fvDggfThSm358uXQ52I7lUrl7OwMQNQWSyMrizZseKLz/eXLNHq0ZkUvAWRlpXZ3X7Fs2dNr39Vq9bRp00QJxIYNG4jowIED4t77iBEjil+q99tvv9WrV+/rr78moo4dOwIofCOR6U4sIR0xYoSEY65atUrcTDuuRQMHot27dycmJqalpbVo0UIc1nRZRvhiH3xAY8ZQRETR2ulduwo+4zIZ9e5Nf/317HWPpdxzE/ctW7aIrH306NFikUHr1q3v3LlT+DmffPIJAAcHB1dXVx8fn6CgoAsXLhS+7E5ISAgNDf3mm2/effddsflW06ZNb4i/QkmJPso9euiaH6vV1KcPASQ+7SEhBJCTEz26x/ASSUlJNWvWBDB//nwiWr36sIVFVqVKUk7rMqmIEqZmzaQcs3nz5gB27typUqlECeCxY8d0GTApiTp1IoCqV6erV4t7L+natWs+Pj6aEuFWrVoFBgZmFvNDzLSybx8pFCSXk1gPJhYrN2xIqaklOy3MnVuwQCI/n8LDqUGDix06vPF0i+4SuXGD7OwkmNco4sqVKyKt7Nevn2Z+/dixY507dzYzMyucrNepU8fLy2vZsmVPr4O6dOmSyNr79Omjydql2phT39544w0Uu7++Fg4ePKipJ5FY4c73q1YV5DK1a5O3NwUHv3T/0tmzZwNQKBSrV68mosOHD4uiGk9Pz5K27snPzxe7sZZsuS17mdu3b8tkMisrqzSJdqM9cuSI+F5rt33Sb7/9JpfL+/Xrl5+ff/PmTfGtnzNnjiSxPW3ChIID6dM54fHjNHZswf4tAPXvf3XmzJnPrN8rtV5U4x4QECDq2t955x2RdlerVu306dOaJ3zyySdPLxy2s7Pr3r27r6/v6tWrz507p/ka3717t02bNgAqVap0+PBhCd/DoUMkk5G5OUlyZ/XmTbKxIYA2bSIiWry4oM9YcYwYMULcncjNzb18+bKVlVX9+s22buWl8qVRbi5VrkxA0a6duli4cCGA9957jx5d1k6dOlXr0RISqH37gjWvWhxV4uPjFyxYIFIrcXOzc+fO/fv393xk7NixojmJ6A6xb98+rXd8ZEQ0Zw4BVLUq3btHWVnk6prm6vpuMZumCbdvk5UVyWR06BCp1QXXbDr07Xjst98IIDs7unlTgtGoUNbet2/fwlUxonZZoVA0bdrU29s7KCjo5gt/ZEZGRu/evXv16tW4cWPjytpTU1MtLCwUCkW81tWZL/PBBx8AmKWPnb0GDSpY43v2LH36Kf38c0lLB8SErlwuF2lceHi4WKLw1ltvlWgxhh4vTsq9Ll26APjzzz91H+rWrVuOjo4APvnkE+1GuH37thjh448/JqK9e/eamJjIZDJ9rA/54YeCe54nTz73OcnJ9PPP1LIlNWniJU6RLi4uxjLD9ZLFqZs3bxbzdh07dhT3s2xtbYusCBblMUql0tPTU5QwFla4E0JaWpq7u7uokFuzZo0kbyA7W92oEQG0YIEk4xER/fhjwUR7idb45efnf/jhh7a2ttevX8/LyxNdiseMGSNZWExqoruxtgeiZxC94aytrTMzM8+ePaupnNFiqPv3qUULAqhRI51awufk5AQHB3fq1OnpLuBFVK1aVTSH5t7w2lGpCu7X9ehB+fl07dp1MQ35yzPriZ8lOpr69KFRo4gepdpOTiTVxdTbb4sdgrN07xlw9epVkbX37t27SJamUqn27NlToprvGzduiKy9TZs2xpK1E9G6desAvP7663oaX6VSiXaQhSfLJFO487225+JFixaJ2/IBAQFEdObMGTGNquk9/wI5OTkhISGenp6mpqa1atWysLAQvZiZhH744QcAgwYN0jzy77//njx5sqSNidLS0kQLfzFfrnU8R48eFXP2q1atokcdmWxsbM5JOHlGFBpKJiYkk9HatcV6/vHjx8eOHau5QV2tWrXSn7u/vKvM2bNnxU3/evXqDRo0yNLS8uQLrmKI4uLidu/e7e/vP3To0IYNG27fvr3wn+bn54tZBJlMplQqda++nTt3bps27l263JWwUEmlou7dCSg4g5aI2BFD7Lpap04dnsIszY4eLVi3J2Gh3bJlyzSrQqtVq2Zpafnaa6+JlqkhISFPLxR5puxsatiQAGrRgor3ipcTuxNYW1sPHjz4vffe8/b2Hj169MiRI8eOHTtnzpxZs2b17NlTHLlMTExGjBgREREhzQ8uT2JjydmZunZdpVQuIKLg4GAxT1Giv8ysLEpJIWdnAkjCZObBA3J1jaxfv5m/v78u41y7dk0klAAsLS3d3Nz8/f0jIiJKejDX9I4U0z0uLi6ltnXMM/3zzz9169bt2rWrnpqnHT58WJxE9LJGpUjne20tWbJErKhZsGABEf37779Vq1YFsHjx4ue95MSJE1OmTBHT8wBMTU3r1q0rfv3+++9LVdfBiOj+/fsKhcLMzEzTZ2z06NHir1rTyDskJCQ6OvoFgxReTqp7OdPSpUsBWFhYiC4LYluSWrVqSXXb6saNghvpJV3OnZKSEhgY2KRJkypVqkyZMkWSYPSnWH3c7927J5oGVKxYUdS0aS0pKSkoKEhcvQEYNmyYLvtWXL161cLCQiaTSb6u5do1srKitm21aR505MgRhUIhl8sPHjwobVRMWmo11atHAIWGSn9q/Pnnn8UdqiIT2wMHHvfwoDlzaNMmesEBc/VqateOJFyhl5eX9/rrrz9zrn3evHniOWfPni3cG97V1TU4OFh/64fKpLCwS+Lrv2vXLnpU7eDh4VGiQaZNI4BcXSVe1L5nzx6ZTGZqavrM3kQvlpSUtHPnzjlz5oj2Lw0aNGjVqlXh9gO1atXy9vZ+8Vx74fYymtfWqlVr3LhxxpW1E1F0dLT4powfP14fubuPjw+07ZJeXFJkycuXLxcltSLUK1euTJ8+/em/ELEPjNiRQGjatKm/v39sbKxarQ4MDBRbZdWuXTssLEz3qJjQu3dvzQw3Ec2ZM6dFixampqZFTgGOjo59+vT55JNP1qxZU+SO69y5cwE4ODhcu3ZNkpAmTpwoJrbv3buXlZUlSjl69+6t+yYkycnJgwdHAzRoEGn3jRRbCrSRfBMWqRV3AyZNlYuZmdmoUaOK1K+/VHJy8u+//+7u7q5ZtCS6aAFo167dsmXLjh07VtKOWmq1WvQT9fb2LtELi+nECdq+/fHmW8+Umppx+/bts2fPHjhwYPPmzStXrly4cKE4sX322Wf6iIpJ67vvbrVvP2/cOIk/Qlu2bFEoFDKZbPny5RcuXAgKCvLx8XF1dbWxsWnbNqHwHvUVKtDs2U98zDTbvAUHS7zN28OHD0VvwbZt22p23/Ty8tok1nM8InrDa6bE6tevHxAQoOW2L+WS2Nm7YsWKt27dys7O/uKLL0r0t3f7NpmakkJBL+shpA3Rza1JkybFuR0s8uwiGyppdkqytrbu0aPHmDFjRowYUb16dQBOTk5Fpoezs7OPHDmyaNEiDw8P8dnTsLOze+ONNxYsWGC8ExwHDx4UDVWGDx8u7d5narVa3OgODw+XcFg9WbdunehBN3PmzCJ/lJWVFRwc7O7urtki19nZ2cfH5+n6n4sXL4oOho0avT1/vpq3kpPEihUrAHTv3r1wtpabmyvOSn5+fu7u7qIxlMaiRYs0z9y0aZNMJlMoFH///bdUIeXm5vbo0QNAly5dsrOzY2JixNFDxy3S8vPz33jjDZlMPmLE0ZJekM6aNcvb2/vhw4diiU7z5s11ieQVKMHOqfn5+X5+fv369dP8A5uamjZt2tTLyysgICAsLCz9+bPTkyZN0uTrPXr0mDJlyoABA0QSX3jDPHH7xs/P7+kGNU8TuzZUrVpVf1M1hTff+vlnGj6c+vWjDh2oQQOqXJlMTKhHj+14Sr169ezt7fWxvTaT3PXr18WO5RImpgcOHBAXpU/XJKhUqqtX8//6iz79lPr3JycnAmjo0McfMyL9bvN25swZsf98kVYYs2bN8vPzu3v3ruaRtLS0gICAWrVqiU91lSpV/Pz8irR/Zc+kVqsHDBgAoGPHjtptmrNrF339teRxERFlZWWJ3kfPPE2mpaXt37//yy+/fOuttzRXboKFhYWrq+vHH3/8008/TZs2TXSA0WjUqNHw4cPnzp2bkZGRmpoaGhqqVCrd3NzEh63IEd7f3z8sLEy7tR+lTVhYmFjJ4OHhIeFmN0eOHAFQu3ZtY+nlGhwcLOZxJ0+erFarVSpVWFiYt7e35pajubm5u7t7cHDwC/7dc3JyZs78tHHjTLHR+AsrOFixREVFWVhY2NnZmZiYiGzN398/JCQkLi6u8NNu3ry5devW+fPnv/3225o7HpGRkeI2iFjDIKHY2FhxXSoWAR47dkycMVesWKH1mDNmzABQqVKlqKiokr5WFGtFR0dfvnwZQOPGjbUO49UoQeIuXL9+XdSvN2jQoMim1gqFokmTJu++++7//ve/vXv3Ft6+8cCBA926dRs7dmz//v01N+IVCsXrr7/+0UcfjR49uk2bNkU6iAGwsbHp3LnzxIkTlyxZEhYWVrgH/IMHD8Q6GElWTD9P4c233nuPCk+Uiv/69TtWs2bNVq1avf7664MHDx43bpyPj4844f3222/6C4xJSNyqk+qD9O+//9rb2wMoZp3c/fv0wQePP2ZEet/mLSgoSORhmsm85ORksTTHzMzMy8urcLd4lUoVEhLSuXNnzdnXy8uLt3Z6qYSEBLEa2MfHJ/EZMhMTqch/aWm0c+dLbvFJ4vTp02ZmZjKZbOfOnVSo1tzV1bXIQbhwnv30TMTt27eXL18+ZMgQ8YHXHLQLnxfkcnmLFi2mTJmyZs2a28Xvz2VUIiIixM2EN998U6plbaLCSusOHgaxbds2kX65uLiItsgAZDJZ165dV6xYUfwtXPbvL9iF2taWAgP1GnIZp1arRaGEra1tkWwNQM2aNT08PObMmbNx48aoqKgil4gPHjyoU6cOgNGjR+sjtsjISHHSWbp0KT06K5mbm//xxx/37t0r6VW9eLmpqemBAwe0CEZ0Tbxy5UpUVBSAevXqaTHIq1TixL2w1NTUsLCwH3/8cezYsa1bt366cKpWrVqDBg2aMmXKoEGDNPm6XC7v0aPHL7/8UmShXl5e3vnz59etW+fn5/fGG2+IuydF1KlTZ+DAgbNnz+7evTuAfv366fb2X6Lw5luHD9Off9I//9DJk3T1KsXH0/Pu5a1duxaAo6MjN6Y1CmJL+bfeekvzyPr164cNG7ZmzZqS9riIjo4WhVLDhw8vftlr4Y8Z6XnGXRDb+tSuXXuP6DpOdObMGS8vL81X2NXVtcjOTQcPHvTw8BBHf7lcPnDgQG4+82LHjh0zMTF5umEugB49tj49C9C37xO3+PRKFPPY2NiIpYQaJiYm7dq18/HxWbduXfHz7Pz8/IiICH9/fzc3tz59+piZmbm4uPj4+AQHBxtRlxhdREZGinqD119/XcfllVFRUT4+PnZ2djVr1jxx4oRUEb4a+/fvNzc3Fx2HatSo4efnp11hdGIiDRtW8L2YMEHyMMuLb775BoCDg8PNmzdzcnIKF22KqfTCbG1tXVxcRAHFgQMHunbtikfVLHoKb82aNQAqVaokenh88MEHhZeEWVhY1K1b19XV1d3d3cvLy8/PLyAgICgoKDQ09MKFC/fu3dOcnjQT9oHaXueJTccvXLhw69YtcUkj2ZvUD50S9yLy8vI0nww3NzfNnVbRgkAul7u6ugYEBNwr9nZESUlJYWFhgYGB4qOm6dcj0mJbW1t97OVUmMio7t8nO7uSvVCsApw2bZp+4mJSevjwoZmZmYmJieZKcvDgweJjplAoXFxclEplcbqCxMfHiwv3Xr16lehgV+RjVjhx11NheXZ29rvvviu2r2vZsmVgYKCYT42JiVEqlZoJ1AYNGgQEBBSeRBRZhZWVlY2NTdeuXfXSqK6sSEpKsra2Njc3r1ChgsOT+vff5uBARf57++0nbvHpVX5+/pAhQ/r27QvAzs7Ozc1NqVSGhobqXjCWn59fNmpgSury5cviTNe1a1ftmomFh4cPGTJEXBsrFIrt27cbS52MRnZ2trW1tUwm27Rpk+7BBweTgwOtXy9JaOXOiRMnxI21rc+5f6dp5O3u7q5p7FM4b65Zs2Yx26Bp7dtvv7148aL4dW5u7rFjx1q2bOns7KxZEfECFhYWNWrUaNeundibWVP7FxsbW9IFJ2LN9NmzZ+/duydSVmnfpuSkTNyLUKlUV65c+fPPP4OCgn777bfi5+vPk5eXd+nSpfXr18+aNSsgIKA074l94cIFExMTExMTaRuUMj0R9xM1lXy3bt0KDAx0d3cX1/FCnTp1vL29Q0JCnpmUp6amtm3bFkC7du1KuszaIPLy8r7++mtNX7/q1av7+/uLe0SpqakBAQGam92Ojo5KpbLw1Gl8fLy4tqlbt25CiTY7KE+USiWAPn36FP8lRe696Nu1a9d03JOVFXb16lVRuevi4lKiWw1hYWHu7u6izY4oV7ty5Yr+4tSfXbt2QbqmHDt3kih8KLJqf9MmiVftlz0JCQlieVLxuxLFx8eHhob+73//GzlyZPPmzYODgyNF3Z6BZGZmRkdHh4WFhYSEBAYGKpVKHx8fLy8vNze3pk2bOjs7F25LZWVlJar2d+7cWblyZaVSWaKfNWDApSZNIk6fzoiPz69RI7ply9J+M1mPiXs59+GHHwLo2bO3oQNhL/fzzz+bmZlZWFi4u7sHBgZqLjJTU1M3btw4ZswYseWbYGdn5+npWbihXk5Ojpi8rFevnr6nKKQltmdq3769eGs2Njbe3t4imcvNzV2zZo3Y7RjA7NmzC79w06ZN4kLFzc2N+0U+LTk5WRTJlGiXaK1v8bFS4tatW/Xq1RPJ60tbU4sFJJ06ddJ8+3x8fAovEDc6ojS/yLFCa0Uqx15BDWEZoVLNGT9e3PyRttlRqZKZmXnr1q1jx46JVWpvv/02EYWFhSkUChMTk2PHjhV/KLFH9fHjlJhIADk46C1oiXDiri+JiYmurr7NmqXrYUNfJrH+/fsD0FzBy+Xyjh07fvnll5plmiqVKiIiQqlUuri4iKdptg5QqVRifwpHR0epOt2+eoXn/ORyubu7+9GjR8Uf7du3b8CAAWJnMeHHH3+UyWRvvPGGuJ7hzqdP02K6nZUNMTExomS2SZMmz7vPnJubGxQUpOnMU6VKFaVSWQZuXomLluPHj0syWpHKMX2v2i87vvoqX6H4+q23/tNlz23jcefOHTFLInYZ8vPzE3eDi3/ru1s3AujwYUpNJYBsbPQZrhQ4cdejFSsIoNde46NMqSa2t6xQocK///4bFBTk6elZeIlM1apVvby8goODNWvObt26tWTJEk0h77Rp08Q0/JkzZwz3JqRx9epVHx8fzTpyFxeXoKCgpyfUb968KdpojB8/3sTERCaTBQcHGyTg0km76XZWZsTGxopNBhs2bHjnzp3Cf5Senh4QECAqasRd/jKzScKlS5cAVK5cWapbcC9YtZ+YSKdOSbw9WRlx6BCZmJBcTrt2GTqUV0d0lbGxsbl+/Xpubq64h/z+++8X8+U9exJA+/dTVhYBZGGh12AlwIm7HqlU1LGjNrvvslcmJSVF9C9avny55sHMzMzQ0FAfH5/XXntNk8GL3d0DAgIK38tesGCBKEvVtGcpA2JjY5VKpWa7nLp16z6dW+zevVvsDzphwgRxxLxw4YKhAi5tFixYJIqIDB0IM5jExMQOHTqI1Fw0lk5OTvb399d8rZo3bx4UFFSWKhn+97//AfDy8pJqwBes2v/iC1qzhrZskepHlRXx8VS9evnMOYYPHw7A1dU1Pz//0qVLopdJMWeU+vYlgHbvprw8Akih0HewuuLEXb8iIkguJzMz4jVgpZMoyuzYseMzuzeq1eqIiIh58+a1a9dOU0gjk8nat2//xRdfLFiwQCaTyeXyMjnfnJaW9uOPP2q6DVSuXHn8+PF+j8ybN8/NzU1cz4juNDVr1ix+q+YyLCmJnJ3VPXpEHDliZL38mLSSkpLEBgiOjo6jRo3SNOB7utdq2fD++5uqVGmu151VNO7cIT8/Mqr1RPqnUlG/fgRQjx7PbVZddiUlJYmGCl9++SU96vJcuXLlkm4aKJMRUNpv5nDirndjxxJA7u6GjoM9JTw8XCxkOVuMzeXj4+OLFNKIHROXLFnyCkI1FLF+ztXVVS6XP90yrLB27dppt0toGaNUEkA8286IKDU1tXv37ubm5uKyX6Tshg5KL1JSyNSUzMwoMfFVpIxvv03ffkvnz7+CH/WqxcTEXLt2TZtbMV9+SQA5OpLOHfyM1KFDh+RyuYmJyYkTJzSbT/Xp06dEF8mmpgSU9rZFMiJ6wcmY6S4+Ho0aITkZe/agTx9DR8Meyc+Hl9f10NAB48cP9Pf3L/4Ls7Ky9u/fv2zZsh07drRs2fLcuXP6C7L0uHLlSlhYWGJiovhtZmZmTk5OTk5OXFycs7NzUlKSp6fnG2+8YdggDS4lBXXqICkJhw6he3dDR8NKgczMzDNnzhw5cuSNN94Q7aLLpI0b4emJ7t1x6JChQzFy/v7+s2bNMjU1rVu3buPGjRs1atSwYcNGjRp1atTIpEqVF73yzBkMH46ff0bfvq8q2FLn448/Xrx4cb169c6ePZuZmdmyZcu4uLgffvjBx8enOC/ftQvOzmjdGtu2YeBAfQerPU7cX4Xly5GSAl9fPLmhODOk77/HRx+hdev8o0dzC+/tVUxnzpxxcXFp06bNmTNn9BEeM0bz5mH+fLi5ITTU0KEw9gqNHYvVq7FoEWbONHQoRu7bb7/95Zdf7ty5o1arNQ+ayOV55uYwN0e9eqhbF02bolkz1K2LRo2wZg2iozF2LJo2RV4entrAvlzJycnp2LHjuXPnJk2atHTp0m3btg0aNMjc3Dw8PFyUdL7YtGmIjMS+fZg4EStXvoJ4tcSJu94VvoazsEDHjrC3x+bN8PAo518xQ/rvPzRrhrQ0bN8Od3dtRrh161adOnVq1aolNklmjKfbWfmkVqN6dcTG4sIFNGtm6GjKhOzs7GuPXLlyRREf/9vJk0hOLvq8vn3h6orx4zF7NkaORK9eBoi1lLl06VK7du2ysrK2bds2YMCAiRMnLl++XKlUzps373kvuXkTe/YgLw9RUWjVCvHxuHaNE/fyrfA1nIMDpk5FnToYNQpLl+LRaiX2qg0ahG3bMHQo/vpLyxHS0tLs7OxsbW1TU1MlDY0ZK55uZ+VTeDg6dEDNmrh929ChlG3x8bhyBVev4to1XL2Kq1dha4vDh2FlhXPnsH49vv7a0CGWCgEBAdOnT69Spcq///5ra2u7Z88esdV3YZmZOHYMe/di716cPg0ATk4YNgyzZ+Pjj3HnDvbvN0DkxWRi6ADKhTFjsHhxwa8//xy2tjh+3KABlW87d2LbNtjZ4bvvtB/E1tbWzMwsLS0tNzfXjEugGNCgAWrXxty5ho6DsVfr778BaHnrkpWAoyMcHZ+4nefjg4cPUbMmoqPRvLnhIitdfH199+zZs2vXrjFjxuzatUuTtavV6jNnzuzZs+fGjQG//948L6/g+Q4O6NUL/frh/HkA8PdHo0YGCr14OHF/FTw8Cq7hHBywcGHBjDsziMxMTJ0KAF99herVdRrKwcEhLi4uMTHRyclJktjYq0RE7u7uAwYMGDVqlKWlpdbjaGrhbGxw/TpM+JjKyhmxPr/cL003hIkTMXMm7O1haoqFCw0dTWkhk8l+/fXXli1b7t69e9myZUOGDDl06NDevXt37NgRExMDoEMHUqubu7jAzQ1ubujRo2jdckqKYSIvJj7JvCI1amDjRrRt+/iRR23B2SuVmoomTeDoiMmTdR2qYsWKnLgbKSKqUqVKQkLCzp0758yZM2XKlClTpjg6OhZ/hPR03L6NW7ewYQNu3sS+fQgJKdWNCBiTnLhq3bIFv/yC3r0NHU051KwZ1q9HejpsbAwdSulSrVq1wMDAd955Z9q0aVOmTNE8Xrt27b59+/bv36ZnT9jbGy4+3XCN+6tw6RKaNYO9Pe7fx6Pt5NmrVniVcI8eEnxpBw3KPHrUYutWuLrKJYiPvVpt27aNjIxs0aLF+fPnAZibmw8dOtTPz6/Zk8vrHjx4cPv27Tt37iQkvHbxYofbt3H7Nu7cQUJCwRO6dsW4cUawnokxyRlLFw5WPs2fP//u3bvr1q3r0qWLm5ubm5ubi4uLoYOSAM+4vwqrVgHAiBGctb9S+fmIi0NMDGJjUa0adu8uOMdINTNKZPXwIR4+lGAo9uo1a9YsMjLS19e3UaNGixcvDgkJ+eOPP9asWdOtW7cGDRrcu3fvzp07N2/ezMrKEs/v1s07LKyD5uWWlqhVC7VqQSZ7XAvHWHlTeAUXgPh4lOTGFWN6pFQqU1NTlyxZYlq2Wvhx4q53+flYuxYAxo41dChlTk5OTmxs7L179+Li4sT/c3OdLlyYevcu4uMRFwfN/aRJk2BuXvQco6OKFQHg0ZZEzMiImfWLFy++//77Xbt2vXHjxg8//LBy5cr09PSVhSYPHRwcatWqVbNmzTZt2g4ahJo1C/J1TXYybRrwaD3TxYuwsEC9eq/+3TBmGIWvWvfuxcCB+OknjBtn6LAYA/Bog/MyhhN3vfv7b8TGolkztG9v6FDKipSUlA0bNnz66acJmnqFR5r/v737D6q6zvc4/jr8UCRAXESBFFsqU4e9upMIBboqhlhq9sNubKAe4EqOuyLJr91AG1dNxl+0joF5GfxZVoubQWyK4DqRmoatBu7soFk6izfFNBH1Jj/uH3AT7222Vb6Hr1/2+fiDOTDnvL9v/jjDi+/5fD7v4Eeqq+e2PXZykp+f/PwUEKDgYNXWGnxn1MdHIrhbVnBwsKTq6uq2b4OCgl577bWFCxeeOHGiqqoqMDDwvvvuGzRokKen5z+uk5vb/uCNNzRihKKi2g/ZAP5FfH8Kx8GDunpViYlqbNSvf212W0A3RXB3uMJCSdyBMFJkZGRVVZUkV1fX/v3733vvvW1f/fz8AgMf6ttXAQHy81P//nJ2vvmqjndGExI0a5ZGj+5UG336SAR3ywr52c/Kx4590OmW/Qk+Pj4+Pj6hoaF3UPCxx+TlpdJSlZbq8ccN6hK4i33/X+tHH+n4cWVl6Z57tGCB5s1T//8+/VxqoJnNAd0Um1Md69w5DRig1ladOSOOHjHE6tWrFyxYIGnFihWpqal3UOH11zV3rnr3Vnm5OrNTJT9fc+boxReVl3fnRWCa1lb17q2GBl240L7sqdNyc5WSogceUHW1evY0pCRwtzt1So8+qqYmlZYqJERbt6pk2bG3zo61vfBL/f73cmLvPmAk3lGO9Yc/vB0W9p/Tp18mtRuiuro6PT1dUkRExJ2ldkkvvqiYGH37raKjdfz4bbzwxg11HJPaFvb+32odWITNpqFDJammxqiSv/qVgoN14oTWrjWqJHC3CwjQo4+qvl7jxqmsTLGxenPZl7ZrV7VuneLj1dRkdoNAt0Jwd6z8/KUfffQfMTH7zG6kO7hy5crEiRObm5t79uy5Y8eOO67j5KRNmzR5surrFRWlU6d+5PnXrqm4WElJGjBAK1bc/Hl4uIqK9Nvf3nEjMFvbrEHjgruLS/vigcWLdfasUVWBu1rPntq+XS+8oMZGLUk60/R+qdO0qSopkYeHNm1SerrZDQLdCsHdgT755JPPP/+8f//+0dHRZvfSHSQkJLSNPVu3bp2vr29nSrm66p13NHas/v53PfbYD2es8+fPFxQUJCau/slPNHWq3nhD587dcof+2DEFBbUfDA9Lajuy3bjgLikyUlOnqqFBWVkGVgXuaq6u2rxZr7x0ucwW5fLMk9qyRRMmqKJCI0YoOdns7oBuheDuQIWFhZJmzJjRzc4Q7WqHDikq6syQIU+88467FBISEm/EVt9evbRzp0aOVEPDf8XHP/9Nh02meXl5Y8aM8ff3T0xMLCzM8PZuefhhLVqkv/5VRUWSdPWqzp7Vrl1KTlZTk95/v/PtwAzBwfLyUkuLsVVXrdLDDx/+9NPwQ4cOGVsZuGs5OWnRKq8es2epqUkzZ2rNGoWE6MgRDRpkdmtAt8LmVEe5fv26v7//pUuXampqhg0bZnY7VhYYqDNn2h4uc3J6uqZmyJAhRtWur9fzzz9ZXv7+qFGj9uzZ03b237PPPltUVOTm5jZ+/Php06Y9+WRMv34ekr75RiUlKinRn/6kZ56Rt7eGD2dkppU1N99y8JBxMjMzc3JyHnnkkY8//thmszniEsBdavVqpaaqd2/97W/y9dX69Tp5Una7+DsIGIHg7ijbtm2LjY0NCws7cOCA2b1YnK+v6us/dXEZ2dR0MDw8rLLS2PJ1dXWjR4/+4osvxo0bV1pa6ubmVllZefbs2ejo6LYc/9VX9X/8Y9/33lNlpZqbJcnJSY8/rvvvV1ZW+8HwFRXGNgUHa229mSfOnFFoqLy9tWOHpkyREZ+PXbly5aGHHqqrq9uyZUtsbGznCwJWUliowYMVHq5163ThghITlZWl2FiNH292Z4DlsVTGUdrWydgZl9p5bm6SRjY1SQrp5NHrPyQgIGDv3r2BgYF79+597rnnbty4ERERMX369NOnT+fk5ERERAwZ8kB2duu+fZIUHq7cXJ0+reLi9pcvX66qKsObgoO9/rrOnVNKilauVEGBLl6UpPfe03ffGVLew8Pj1VdflZSWlna541FEwL8Cu13h4ZK0Z49SUxUQoORklZWZ3RbQHTCAyVEmTJiwf/9+Dw8PsxuxvEsbNlQ98URrS8u/TZzYzzE7/gIDA8vKysaMGVNcXDx58uRhw4bt3Lnz1P8eN+Pl5RUf/1VY2H2TJqnj+OTvh498+60jmoIj7dmjbdvk7q7kZM2erZdflqenDP1wLC4ubsOGDZWVlcuXL1+2bJmBlQHLGDhQ9fUKDNTJk+2HOAHoHIK7o7S2tl67ds1ut/fq1eupp54yux0L846O9ti/v6KiYsJvfuO4qwwePLikpCQyMvLYsWO7d++W1Ldv30mTJk2fPj0qKqon03S6mY55om9fLV2qn/5UM2YYeAWbzZabmztq1KhVq1bZ7fYHH3zQwOKANSQlKT1d3t5yddXSpWZ3A3QHrHF3oOzs7CVLljg7O2/cuJF1rpZw5MiRixcv7t69e9q0aaGhoU7M/Ouuamr0u9+15wmbTSkp7cE9P1/u7gZex263v/nmmxEREW+99Va/fv0MrAxYxpUr4sNnwCAEd8datGjR4sWLnZ2dCwsL4+LizG4HQAeOzxNff/11enr65s2bBwwY8O6774aFhTn0cgCA7o3g7nA5OTmZmZnOzs4FBQUzZ840ux0AXer8+fMxMTHl5eUuLi5LlizJyMgwuyMAgFWxEsDhMjIycnJympub7XZ7Xl6e2e0A6FK+vr4ffvhhRkZGc3NzZmZmXFzc1atXzW4KAGBJ3HHvIitXrkxLS7PZbGvXrp07d67Z7QDoatu3b09MTGxsbBwxYkRRUVFQUJDZHQEALIbg3nXy8vLaIntubu68efPMbgdAV6uurn766adra2v79OlTW1vr4+NjdkcAACthqUzXmTNnTttSmfnz5+d+fwb4j2ptVX6+0tJ0/Lgk7dqlS5ckaccO3bjhkEYBOIa/v39kZKTNZvPy8iK1AwBuF8G9SyUlJa1fv95ms6WkpCz9Jw+17TjisaJCZWWGT3kE4HDXr+evWhUUFJSfn+/i4jJlypTveP8CAG4TS2VMUFBQMHv27JaWlsWLF2dnZ//AM+rrdfiwDh/WyZO6fLl9xOPRo9q+XTduqK5Onp6qqNBf/qJ77uny9gHcpuJizZ//ycCBYfv2TZgwYc2aNcFMkQQA3D4mp5ogISHB3d19xowZCxcubGxsXL58eUNDQ1VV1eHDh6/X1maXlenLL28+Oz7+lpHRn33miCmPABziwAG99JIOHpQU4uPz5/LyX4wfb3ZPAACrIribIyYmprm5edasWTk5Odu2baurq2tpaZHk5uKS3aOHPDz0858rJEQhIbr//ltGRn/22c0qNptpvwCAf8aaNTp4UP7+euUVp4SEXzg7m90QAMDCWCpjprfffnvr1q0ffPCBq6vr8OHDQ0JCRo4c+e+jRrkPGaL/8weekdHAXWvXLoWGyttbO3ZoyhS5uGj9ep08Kbtdbm7atElpabx/AQCdR3A339GjR4cOHdqjRw+zGwFwR1JTNXdu+wK2vDxt3KgLF5SYqKwsxcaKtTEAAIOwVMZ8w4cPN7sFAJ3z8svy9NSBA5K0Z0/7hvLkZG3fTnAHABiF4A4AndZxy/jAgbdsKAcAwCAEdwAwjs2mpKRbNpQDAGAQ1rgDgAOwoRwAYDSCOwAAAGABTmY3AAAAAODHEdwBAAAACyC4AwAAABZAcAcAAAAsgOAOAAAAWADBHQAAALAAgjsAAABgAQR3AAAAwAII7gAAAIAFENwBAAAACyC4AwAAABZAcAcAAAAsgOAOAAAAWADBHQAAALAAgjsAAABgAQR3AAAAwAII7gAAAIAFENwBAAAACyC4AwAAABZAcAcAAAAsgOAOAAAAWADBHQAAALAAgjsAAABgAQR3AAAAwAII7gAAAIAFENwBAAAACyC4AwAAABZAcAcAAAAsgOAOAAAAWADBHQAAALAAgjsAAABgAQR3AAAAwAII7gAAAIAFENwBAAAACyC4AwAAABZAcAcAAAAsgOAOAAAAWADBHQAAALAAgjsAAABgAf8DCKUfk1mxnRIAAAAASUVORK5CYII=\n",
"text/plain": [
"<PIL.PngImagePlugin.PngImageFile image mode=RGB size=1000x400 at 0x11E8C4F98>"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"Draw.MolsToGridImage(mols[:10], molsPerRow=5)"
]
},
{
"cell_type": "code",
"execution_count": null,
"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.6.8"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment