{ "cells": [ { "cell_type": "markdown", "id": "be5c940d", "metadata": {}, "source": [ "# Clustering with sample weights\n", "\n", "Like HDBSCAN, PLSCAN supports weighted samples as input." ] }, { "cell_type": "code", "execution_count": 1, "id": "0ea6577f", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from fast_plscan import PLSCAN\n", "\n", "plt.rcParams[\"figure.dpi\"] = 150\n", "plt.rcParams[\"figure.figsize\"] = (2.75, 0.618 * 2.75)\n", "\n", "data = np.load(\"data/clusterable/sources/clusterable_data.npy\")" ] }, { "cell_type": "markdown", "id": "54e2978f", "metadata": {}, "source": [ "The leaf tree without sample weights has many leaf-clusters:" ] }, { "cell_type": "code", "execution_count": 2, "id": "8dbf453b", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc8AAADuCAYAAABf/4wMAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAXEgAAFxIBZ5/SUgAAZCtJREFUeJztnQeYE9XXxk92l96lI71XEQsqihQVxAKKiIINKaJSVMCKimL5fzQBBRUFKyI2VFCUKqKioICCFOm912X7bpLvee9yw93JpM5NZiZzfzrPLMkkmcyZ3PeWU1xer9dLCoVCoVAowiYp/EMVCoVCoVAo8VQoFAqFIgrUyFOhUCgUighR4qlQKBQKRYQo8VQoFAqFIkKUeCoUCoVCESFKPBUKhUKhiBAlngqFQqFQRIgST4VCoVAoIkSJp0KhUCgUEaLEU6FQKBSKCFHiqVAoFApFhCjxVCgUCoUiQpR4WoCuXbuyTaFQKIzyxRdfUCXX+fTLL7+oixlDXKokmfk0a9aM7Tds2ECJxk0lutPWvI20Jec/s09FobAM7dq1oyVLllBycrLU983MzKRyxc+jbMqi0lSOTnmOk8vlkvoZinxSyOYcPXqUxowZQ/PmzaM9e/ZQsWLFqHbt2nTNNdfQuHHj/I7HcePHj6e1a9eyf1900UX0+OOP04033hjwMyBqL7zwAi1btozS0tKofv361K9fPxo6dCglJanBeyDcbjdlZyZTDW9zqk75HYRgtGjXhMYtGkWJjsfjoUWLFlH79u0N3T94nwULFlDHjh0jfh/cx7CPnbmj1kBKcSXT5/vfNe0ccnNzqVatWnTgwAFKSUlh/y5evDihTPLPP/9M1157LbMNL5tcqFAh9vdPP/0UE1GbOnUqE86iVJxS6SR98803dOutt0p7/9WrV7N7d9WqVWzbv38/ezxUWegPPviA3nzzTdq4cSMVLlyYLr/8cnr22WepTZs2ZFdsPfKEITt37kzHjx9no7fmzZtTamoqM9C+ffsoLy+vwPGTJk2ixx57jN3kuKmLFClCCxcuZL21N954gwYPHuz3Gb///jsTYhzTunVrJszLly+nQ4cO0e23306fffaZ4R9Boo48cf27FO4V9vHNr25CE5eNpkQHPzkIHxpVI/cOrm+HDh1YQ4x72mlcl3Q7W3halPeFqbbMzs5mbQn+Rttw9dVXM7uKNoaoQjRycnKYrWTYXwvavgplKlISJdOl1IFW0AIqRiXodN4JaSPcW265hb799lu/x4PJyKOPPkqTJ09mA5tOnTpRVlYWG3XjNV9++SV7TzuSZOcR5/XXX89EDcb8999/afbs2TR//nzatWsXrVixosDx//33H40YMYLd5LjBf/jhB9Yr+/vvv6l8+fJMVLdt21bgNbjh77rrLvYZr732Gq1cuZKJ5datW+mKK65gawsffvhhnL+5fYi4X2bfflzEoocGVNu5ixS8B9a1nCicHJcF7vE//vjD1yHCLBb2Whvjb7Qj/DEZ9tcyYcIEyqUcqktNqairONWkBpROqfTJJ59I+wy0e8899xzNnTuXDh48yNrTYCxevJgJJ9rYf/75h7W5P/74I2uDIej3338/nTp1iuyIbcVz1KhRdOzYMTY1q+dsg1GiCAyIaaoHH3yQ3QCchg0b0siRI9mNjGNEvv76a9q5cye1bNmSiSunZMmSNGXKFN8Nq9An0mlBt+TGxKqIDanR64t71+7Tr0Ywu78lCibs+eeff7K91saiyMqyv3Yw8croV6k4laKqVIs9VosaUgoVooH3PcRGvDJ48sknafTo0XTzzTdTlSpVQh6PQQfAFG2DBg18j6MNRlsM4ZwxYwbZEVuKJ268mTNnUokSJVjPJRy+//57tu/Ro4ffc/wxrIeG+xqsldatW5eNeDHSVfjjzousZcvKkvMDtzpiQ6owejHNvYKiYKIT89JLL7G91saiyMbC/v/73//ITXlUj5pRkiu/WS/kKky1qRFlUQa9++67prTTS5cujbjdtQu2FM+//vqLzpw5Q61atWLz6JiCHTZsGD388MNsXROL9yLo3cCZCOA1WmrUqEEVKlSg3bt3s3UDDqYZuFDqwR9ft26d1O+XMHg8ER2e4xDxFBtSI2Da65VXXpHusakIH1EwAZ+C1NpYFFlZ9ufs3buXJk+cTKWoLFWi8ws8V4PqU2EqQsMGD6f09PS4mva///5j68EVK1ak6tWrJ1z7aUvxhEMQqFSpEltsvuGGG2jixIn01ltvselVeMN++umnvuO5cJYrV46NVvXgxoWAal+nZ/hAr1GcwxNhz9rrltOYWB2xITUClhratm0rfe1MERlcMNGJwTIO9lobiyIry/4cTKN6yEP1qLmfA1KyK4XqUBPKoWzmFBlP9oRoP9EWly1blk6ePMkGQ3bDluKJiw2waI3FZ7hnHzlyhE2fwikI0wX33XcfcwbibvkALuSB4KIqGjHU6/ReE8qrVm/bvn07JSSRiqfHmaOVaIG3JtY8VbiUeYiCiU7Mddddx/Z6NuYiK8v+AE6OM6bPoLJUgcpTZd1jzqe6LHTl6aefpsaNG+u2QbEgLcp21y7YUjz5dAduUvS6MF2LqQHEW8GBCCEk8JTVi/NUxI+cnAgdhnIdop5CQ2oEiCY8zpV4mocomLDDoEGDfPYQbSyKrPY5IyC8zUteqka1A4a9YA20KtVkf6NdVMjBlj7u8Hbl6DkM4TGEkSBIWTw+IyMj4Hvy9YBSpUoV+ByMcgO9Tu81wQgUxxmrnp/ZeCOcTnSK1yhvSI0ixg4i+N5J5ORYQwREwcQmxiyKNmYxz126MP8MTNfK9tJ3hQjacZ0dJyFUpEWLFhQPSkbZ7toFW448McLk0wEYcWpBIgOAqVxQs2Z+rwtCGGjRHEkVxPcWX8efC+c1inN4ItRCd67bcaMVI6ARxtSYE+M8005YY5qPCyb26MwULVqU7bU2FkVWlv2tTs0Q7SfaYozA4YuixDNOcI9ZrG3Cm0vLiRMnCvR8sCjNDcnT8mm91RAzChEsXbq073HEd4I1a9bongd//IILLpDwrRROcRjSTu9FC0IdcA86MeTlxOHTZAVEwUQnBunqsNfaWBRZWfa3Oo0aNWJJFBCDytP4JVL7aUvrQQghbDx/pBb+mBiWwnPXIh2UFv4YAn9Fgr0GIrxjxw6WEpCPdBUFcUfoiu9xOUMExIbUCJjmRuozp0x3i6SdPDeDJCvkIxpEwUR7hDYBe62NRZGVZX+rU6xYMZZ3GWAZLdx21y7Y1npPPPEE28O7FmmiOPCw5esJyGDBeeSRR9ha09tvv80ClDlItYdYOdz8OEYECZXr1KnD4j0RCiNON6DnCIYPHx7Db2lzIhxJmp1qLV6IDakRsM6J/M5OW+8EmWeyfH/Lyp4TDaJgohPTp08fttfaWBRZWfa3A8OGDWP7l19+mbW1Ys7wadOmsVlBFNmwI7YVz969e7NwlPXr11PTpk3ZKBG9HGTrx7TtgAEDmNetOIUA71tM8yI2DrGh6P1hBIvE8kgjhfhQETRKyGSEHhRuArz3HXfcwdJMwfjIkIFzUOjjdkc2knTnOWMEJTakRkAjjY6fE0eeosOQx8TpflEw0V7AKRB7rY1FkZVlfzNA1jW0g5ef3XjHRXyMZ2YDKMCBQQna2AsvvNAXl4/k+Vjzff/995mA2hHbiifAhX/nnXeoXr16rFwYSuQgawXK3+BxLUiggNhQ5FVEQm1k9r/kkktYeqghQ4bofgZK5iCg+bbbbmMxVXj9eeedx8RWRkWVRMbtzo1pOj+7IjakRsF6khNJO54fQwjOnIhv5hwRUTAhoJj14un5RBuLIivT/vEG9xsKZKw8u/HvID6mvSeR9Q1tdZMmTVg5Mww8IKpIDm/Xiiq2L0mWKCRqSbI9m/dRv6bnEuqHokjJwvRdqrwKEFYF03XohWOJwYlTrjL4avJ39PZj+RWNPtz6BlWrFzpJeSyAGEIcsPYMsJyDpC1Yhw1kY5n2R0UpCFAzupSqugJ7/e/wbqQdtJGlwotXqEqiY+uRp8LaeDyR9cuSE9yBQm+0YgRMeyFJQqKHPOjhzjs3Ves22Uubr+VxnwrstTYWR6Wy7K8wl5i0VnCwQRZ/ZPrHNCcH641i4nVFYhPpGqZTJkHEhtQI8NbE2n6ie23qcerwuQw9p46YVw9SFEx0YjAK5On5tDbmIivL/gpzSZKdRR9rhFh3hKcrarghowVn1qxZLCAW+WgViU+kYuh2SnJboSE1AkQTa/FOFM8j+86tq504aF7MpyiYsAOcFLk9RBuLIqt9TmFPpP3qkGgAHlQIA0HcztixY/0az549e7J0Yl999ZWsj1VYGI8wtRYOTkmSoG1IowVrZ0is7YSQBy2pJ845DGWcMc9hSBRMbHfddRfba20siqws+ysSRDyRoB1ZeqZPn85Gm3rxj/ihY6EcHlmKxCfSWdhI0/nZFbEhNQJCHVBJyI4hD0Y5vNsaI09RMNGJKV++vC89n2hjUWRl2d8sNm3axL5z1apVWQYhJIkZPHgwa/+dhDTxxFQs0iz17ds36HG40HqpmhSJR6SZX5KSnRH2o53ek1Ez0Wl4ss/1tHKyzUuSIAomOjHwoOXp+UQbiyIr0/7xZunSpSy8b9asWSw+86abbmICCg9jZHQLlMc2EZFmPSRhRyKCUOAmC5ZlX5E45GT65x0OhlOSJIgNqREwcsH72HUEY4T09Ezf35lnzv0db0TBBNwhUmtjUWRl2T/eoN1Gchrsn3/+eTYCxRLc5s2bWaY3CKddswVFgzTr4cYIpxe8ZcsWNtxXOCucIByc4i8kNqRGQKgDGi4nhjxkCIJ56oi5SeK5YKITg0xn2GttLIqsLPvHmzlz5tDhw4fZIGnUqFG+x5Eo5tVXX2WzigsXLmTRFk5AmnheeeWVLBMPbpBAIGH7v//+S+3bt5f1sYoESgzvlFAV7WjFaIC+E0MevELh9IzUc3lu440omOjEYA2ap+fT2piLrCz7xxvkUQZwDE3SjJrxnaEBPHGDE5Amnhi2o/Hr1q0bK/iq/UFjrvyee+5hNwzPxqFIcCJMkuAkZMQ74/eGYHundDoC3VrZWZEtD8hEFEy0efCi5W2faGNRZLXP2QVeCxnhhnpgNA3UyDNCLrvsMnr99dfpwIEDbBEZi8kYzmNOHBcbxV/x3JQpU2xbv00RGREn7HaGv5BfQxot6IgiGYndRjBSEMUz07zpT1Ew0Yn57bff2F5rY1FkZdk/3lSsWJHtd+/erfv8zp07gz6faEhdsX744YdZwnXEeUI4cROdOXOGZRbq3Lkzm7Z94IEHZH6kwsLkRuoA5JABlNiQGgGNL5Ik2K0RloJwr6SfSjPvNATBRCfm448/ZnutjUWRlWX/aGB5cJs189vCAdO1AFVTjmnCUhBBgaTvAG2+E5Du7oWSNIjzPH36NFtcRq1NXMz58+f75sQVziDipAcOcRjSTu9Fi1OKKociK8O8aVtRMEXvZ62NRZGVZX+RpOpVKbl5o4Cbq1IFw5/RqVMnlj0uLS2NunTpwqpY4W9UScG/xZhWJxCzb4mRJ4b5lStXVpk0HEqk3rZOQWxIjYBGCn4ETmmsApFnosOUKJiwA5whsdfaWBRZWfaPBgxssE6u3cJt0+Fxi5HqX3/9xZbqSpUqxVKyIlTxhRdeCLommmhI+9XVrVuXGjZsGDL13uOPP87qbyoSH6fEbUaK2JAaAaEO6KDaLeRBNnkZ5n1/UTDFRP1aG4siK8v+ZlCrVi3mKfz5558zx8+BAwcyX5eNGzeygRIIdxrY7kgTT8zho1g08te++OKLAY/DXDmOVTiBCHvWDnIYkpHcAI0vwsPs2AjLJDvLvAxDomCiE1OjRg1fej7RxqLI2j25Be6322+/nSZOnMimn4cMGULnnXcerVixgj3vlFBEqfM9V111FdWpU4flucXFzcw0L/OHwny8kYaqOMRhSGxIjZKVZV6Mo1Uwc3lAFEyICpxmeHo+0caiyMq0v1U4dOgQffnllyxcpXv37uQEpFqvfv36bBEZXlmYvoWYOinXocKZGYPMqsOJkQtCwOw6gjGr6LpMRMEERYsW1bWxKLJ2rsOKJDfaDtu+fftYfD8cQ1GntFixYuQEpFsPw/fFixezufC1a9fSpZdeyryxFA7EIdOwkSI2pEZAqANKAToxPZ+heGLJcMFEJwbtHU/PJ9pYFFlZ9jeD8ePHs7XNDh06sDy31157rW/Q9Nxzz9F9991HTiEmXR/UqXvrrbdYQoTjx49Tx44d6YMPPojFRykSqKqKU9COVoxcXyRJcPx1NvE2EwUTnZijR4/60vNpbcxFVpb9zQChUVdccQXLqYxp2vXr19P1119PP/30E1uucxIxtR6SJjRp0oR69OjBsu2vW7fOkXk4nYqatg0Mb0iNANFctmwZ3X///bacAkwEuGBye3zyySfMGQhhHaKNucgiyxpEU4b9zRJPbIoYxnlyMLzHkL5x48Y0efJkmj17trruDsGJOVcjHa0YAY0wGms7jmASBQgmwk6wx4Y4Suy1NhZHpbLsrzCXuHRXEde5cuVKuuGGG9QN4yRUYnhdxIbUCGh8kSRBNcLmIQomOjFwlOTp+UQbiyIry/6KBBFP3BTvvfdewOdLlixJc+fOpeXLl7MKK4rEx6PmbfWvi9CQGgFTg0h5ib3CHETBRCcGXrTYa20siqws+yvMJa4LJfiRI3ylXbt28fxYhUl43WraVg+xITXqmPfggw86Lv2llfwmRMFE+4bsOthrbSyKrCz7K8xFeRkoYobb5BACqyI2pEZAqEPt2rVtGfJghNxc66wVioKJTgxS1mGvtbEosrLsbxbIatWzZ0+qVq0am3pG+cm2bdvS+++/7yg/hyQjuWyxlslruOHf4W4qt60ziDjDkEMQG1IjoPFFeSi7NsLRkpNlnc6CKJjoxMAxkqfnE20siqws+5sBRB+hKl988QVVrVqVZRNCpZU//viD+vbtS3fffTc5hah/dTw/Le/1qny1CitPr1kJsSE1SunSpclp5GSaV4JMC9q/Fi1asHhH0fsZU7KijbnI8t+FLPvHE4g9wg9x/p988glLksDZtGkTW5KbNWsW9e/fn0VZJDpRjzz5wjcqqYj/DndTJD652dYZIVgJcbRitDG78MILbTmCMYSFpga14UI1a9bUtbE4KpVl/3iDxAgoPdaoUaMCwgkQz89HnZjWdQJqzVMRMzyqnqcuYkNqBKw3IYOX00IerNb35oKJTgzWoHl6PtHGosjKsn+8KVKkSFjHlS9fnpyAEk9FzHCS84AZyQ14RhunzeTkZppXgkyLKJjoxKSnp/vS82ltzEXWrsktuL/Kf//9x6ZnRTBtO3PmTFYI+9ZbbyUnIE08t27dSh999JHPgYiDheTLL7+cxXk2bdqUVSJXOAPlLxQY3pAaur4eD3PccJp45mRbRzxFwYQd4FDD7SHaWBRZ7XN2AdPMH374IfOuveuuu+jiiy+mO++8k+Uuv+CCC6h69eq0ZMkSVhzECUgTT5SigbeVOIV0+PBh6ty5M0vPh9qemDO/4447aM2aNbI+VmFhvG6HrcWFibYhjRaMXBAvaLcRjFFys63jiCYKppioX2tjUWRl2T8akJcWzkraLVyQlOPnn39mo1C045999hlLCo/cyiiPh8edgjTx/PXXX5nzAnofHGQcQo23YcOGMfHEqBM31muvvSbrYxVWxqtWBfQQG1IjwOsRSRKc5tWclWUdb1tRMNGJ+fHHH33p+UQbiyIry/4iWVWKU1qDsgG3nPJyEtF/+umn1Lp1a1ZSDSlX09LSaMuWLdSnTx82gMIoNDvbOvaJJdJat4MHD1KtWrUKPIYbCYvML7zwAhUuXJj1ei677DJ20RWJj5fUmqce2uk9IzRo0ICcRl5WnmXW10XBRCcGXrTYa20siqxM+0cKZio2bNjgt4W7NId6nRUqVKDvvvuOiWiJEiXYPTht2jS66aab2Gg0WJrWREKaeKK6uOh6jd4HXJYhlljv5NSpU4eV5VEkPirDUGzrcOL3Nnz4cNuFPBgl9USa32PpZ9JNORdRMEHFihV1bSyKrF3rsKIiFjyEUb+zpNCmc5B1CCB/uROQJp6YrkW9Ts7ixYuZoGIYL4LpW/RWFImPClXRR2xIjYCGDOtVdgt5MEp2uv+0YG6OeVPXXDDRiRk5cqQvPZ9oY1FkZdk/3uzbt4/ty5Qpo/s8f/zkyZPkBKSJJ0QSw3rcIPPmzaMnn3ySZdDo1q1bgeOQiQPz5YrEx0v26lmbNVqJFjTSH3zwgeNGnu48/+vmNWndVxRMdGLggYq9no25yMqyf7ypUqUK2//111+6z/PkCHCGcgLSxPPpp59mLsxvvPEGW9vcuHEjG8a3bNnSdwzm1rdv3848thSJj1tVVQkIb0iNgM4pvBvtluYtFssBHpOKEIiCCQGdNGmSrzMj2lgUWe1zdoEPhDAt+9Zbb/mFJE6cOJH93aNHD3IC0sQTcUv//PMPPffcc8wDcPr06SwQWGTt2rXMAHxuXJHYuJS/kC7ahjRaEOpw/vnnOy49X06Wf5xnbrY510AUTHRikCSdV1gRbawVWRn2jzf4biNGjGB/P/zww9S8eXPWliOnLQZE8CB+4IEH6NprryUnIHXSHeue8KwNBHIfOinrvtNx28whIl6gAUXiEPTWjYQr4LXwK3AaGacz/R9LM+c6iIKJTgwcaXJy8sVdtLEosrLsbwbjxo2jNm3a0Ntvv02rV69m2YZKlSrFajSjUkyvXr3IKahAPEXM8NhsTSdeaKf3osWpRZWzdQoOmFWEgAsmT8+HkBmtWALtqFSG/c0C6fcWLFhAx44dYx2BEydO0NKlS00Rzn///TfsY+EfIBMlnorYoaZtdREbUiNANKdOneo48czTqedp1rStKJiwA7xosdfaWBRZWfZXEMtgF045TCwjolSaTJR4KmKGRzkM6SI2pEZAqMOiRYtsF/JglBydUWZejjniKQomtldeecWXnk+0sSiysuyvIDp69Ch16tSJpYINxDvvvEMDBw6U7qSlxFMRO7zOGhGFi9iQGgGhDkiSYLeQB8PoZRMyKcOQKJjoxPzyyy++9HyijUWRlWV/BbFiJDt27GCJG06fPu13SbA2+9BDD1GlSpVYDl6ZJIR4oqYhLg6mQerXrx9y3htppdDzQ/b/G264gVasWBH0Nb/99hs7DsfjdXg9jKYIjkeVVQlwXc41pEZBeJjT0LtubpNqx4qCiU4MvGh5ej7RxqLIyrS/07nzzjtZeCQiPW6++eYCDnRY0oBXMOJTly1bxmqoyiQhxBO9byxehwKByffffz9bZIY7NUQQ015XX301c7zQAzko4UmGmx1ld9DD4Tkeudu2Qh+Vnk8fsSE1AhxOEBpmV8eTaMnVmaLN00mcEA9EwQxmY1FkZdnfDCBCGKS4QmyjR4+O2zlhZInPQ3GS2267jdli8uTJNGTIEKpWrRqrAtOoUSPpnyttseT111+n4sWLS1+UDQXqx6HGHOKLMLcdCKQLxAVFlfPff//dl1Abf7dv356JKvZiTx5eZCizBmNARLt3784ex/w6YptQRQDJkPG6RGHN0nVUpmJpqtn4fMPvVe/CWuSJYOo2KdllKNXcgo+W0ZtD3yNy4X8Xm8rDZF73R2+k/q9YJ0SKN6RGwZrZFVdcwe5hJ6175uokhs/NMD9FIToxEEWOaGO0Ic8//zxr5GXZ3wwwisPAQQ98RxTEBm3btqV48uyzz7IZSLTxl156KRuJIgYaU7WhZiOjxiuJlJQUb9euXb3xJCMjw1uvXj1v06ZNvVu2bEE7yf6tR5cuXdjzEydO9Htu6NCh7Lnx48cXeHzMmDHs8W7duvm9Zs6cOey5m266yfD3wPljswKdUm73Ni7ekn23WG5XXXlVzD+Db3Vd1ri2nLy8PO8zzzzD9kZwu93en376ie2dxMibX/Ve6+pRYPtu+iKzT8ubm5vrveqqq9hea2Px37LsD7755ht2jze45E7vld3HBdxqNOnEjlu3bp03VsyfP599Ro0aNbwej8drBvfee6/X5XJ5a9as6d2+fXtMPytFZo+kaFE5NePC5cUXX2SLxRiWB1t8RzJ6xCEFSh2FxzByRk5eTAFzvv/++4CvufHGG9n35Qnw4/3dY+khe35GfVqXk5+nMlbAYUK2q37X8+6hnDT/EYjLm5ghAbh+CLZ3WshDelqG32NZJiVJwOi/Q4cOvoLQmMLFXlsiTRyVJqqD18yzo8677rorZvckZgKDwR2xMNp8+eWXCzyHc5oxY4a0c0mRGW/zww8/sOwaqN0Za1DBBdOmmG7FFEGwWB9kwUCJNLgqi8W6OYi54u8pgqG/+LwIviPSUyFJMorBYj3U7mSkZfpuMjt6AnrziJJc1l/G107vGWm4ixUrxn5zdrRXtHhy/ZcCvCblthUFExt8IjiijUWRxbStDPtbifT0dPr222/Z3/fcc0/MPifcRAd6nrWWFU/cDAsXLmS9DoziqlatSrECC+1YW8X65NixY0Mev2fPHrbXE06AEml4L5TSOXPmDEs3lZqa6nN9DvQ6PA7x3L17d0KI54YV//n+PnM6jUqV8a/ZZ2Xs4oChbUijBa/Fe6HRdhJ6MZ3uPHNCVUTBxHo9OtXozKChFm0siqws+1uJOXPmMAFt1aoVNW3aNGafIzvcxAgpMquqoIIKLiKmOzFaQ7J4velMoz0AuCaj/M3777/PHIBCkZaWXzwXDk2BgICeOnXKJ578NcFex+uS4jXhgNqLeqDSTL169chsdq3P72SA4/tP2k48vTZJyiA2pEbAFBWWI6655hpHTd2eOZUWcNYk3oiCCSFEu4E9bCPaWDsqlWH/aEDFK702GRWvZEzZ3hPDUSdA5EPCiac4nMYaIGInA8VPGhFPjCLhWYWL2KdPn6jPV+HP0tm/+v7+84e1VLupqrsaC7QNqZGRNgoxYBTjpNFnVrr/+mZ2pn+B7HigFcw1a9awCiNaG4siiyl2GfYXyayQRKdrBw5Zyj4Yu87VwYMHWdQDliOclBg+xW7D6UGDBrEbEJkjwgWJDUBGhr+jAQdTDgCjTvE1/HWlS5cO+ZpQBOrdBRqRxptDu86luNqzOb9qvK2wx8DTryGNFjTYSODhNLxufyHw5JqT6k4UTF7kGtVS8LdoY1FkZdk/GhDP3qJFC6nv+emnn7Lviw4BL5htBpgOR7gK/FsCgdlQy4lnvIbT3333HVufRM1QEZ5ZYv/+/b64y9mzZzNj8gu2b5++IEAEMWVbrlw5nxBCLMuUKcPWPfE6vXl8/n61atWiRCAj9dzU15G9x8l22EQ8xYbUCGiwMPLE5qRECXluf6HMMSkxvCiYEEKU6eKiKtpYFFlZ9rcKM+M0ZRsIRDzAsxY2CBYnzsvGycKW1oPQITxFD4gof44LKrJLFClShCURhrgieFYENzXQOv1gDRdV0/G8VjxhJGQqwvpBw4YNKRHwCh70WRnOqxMZL8SG1OhapZGkEnbFneMf6pGn81g8EAUTQvp///d/9NRTTzG7ijYWRRYdHVn2N5tNmzbR2rVr2Uwd1lPjDQZTKJGG64vBT506dcKeCTSK9IUSnuUBXrcIXxG9YTFtOXfu3KDTp+E0PHrbzp072fNwvOGP1a5dmz0Gd/6OHTuyv7/44gu/9/zyyy/ZHrkRtbGc4vNao0GckeYvUWI8vUIu2hMHTpp6LokMb0iNxvuhEUZj7aRRJ9AbPeSeLUAdb2BDRBpwW6KDrmdjLrLYy7K/Ffj444/ZHtnXgjlkxjLWH2v/EydOpCNHjrDoBywhBtosK54Qprp169KwYcPYPDiG05s3b/Y9j1EfegnwyI03OCeA4T1y03KQ2mzatGlsKrhfv34FXoNwGEzfIn5JPGcY6YknnmB/i0kVEolDu/IbAYV8xIbUqIhgicJppa30QlUyz5jjMCQKpljkWmtjUWRl2d9svF4vzZo1y9QpWwzIkKLykUceiXsnUpp4QoR69+7N5vGRvGDVqlV+WTbgUo91RDPEEyNEXGCMjC+88EI2xYBKKUgKj8YHYS/aChWoovLee+8xzzlkGcLo9fbbb2fTwNu2bWOCnEh5bQtgj5BJW6IdrUQL7svHH3/cUZ62+bjCShYfD0TBRDsCpxns9Wwsjkpl2N9sfvnlFxbjjmUwPrMXbzBdLNMJyJQ1z1dffZX9iFGlRC8jD8ANhuewVmgGuMkhnFOmTGHnCY83iCoqU7Rp00b3NcjSj3VPviANDzmsfw4ePDhggmQ7kpFuTpycU+ENqRHwe+NLC05CXF7gZGeZM/KEUKI4BJZxYI8BAwb40vOJNuYiCyCaMuxvFUeh3r17m9aBQ/uNqVpbiydiOjF8DiScHHi/rly5kmSD9U3tSFcPxIZGGh+KhX2kHkxkju0NXdJNIQexITXqLIT1dqy9230K0OisSHaWOWueomBiQ2ebI9pYFFnMzsmwv5lkZ2f7fEHuvtu8ikVjxoxhVVSefPJJNpqPpweztE+CExByx4YCKfAU1uPMifyYVREsxDtvSjD2aBvSaMFr4XmeKCEPRsjJME88uWCiM4OlH9gEXrSijUWRlWV/MylSpAgr2Wg2WG7r0qULjR8/npWNxDIa0qbqtVuwCWYZZSHNcpj3DpXiCSNDTNnCnVhh3RhPDhoD/EgUchEbUiPg94RiBpjxsXvIg1FyTRp5ioIJIUSBCp5tSLSxdlQqw/4KYjHOuPdxvVFhC1sgLCueWCh/6623WGKCO++8U/eY6dOn0969e9kcucJapOuUdMrOyFHiGQO0DWm0YO0MGbeQ59npDXFenjkebqJg8vShFSpU8LOxKLKYYpdhfwWxkadZSBNPBAbDbfnee+9lQbMISeHZe/Dvr7/+msV8Ymr3sccek/WxCkmcOug/nZ56LJVKl4tPwLGT0Dak0YLX/v3331LPza6YFecpCiamYxHfvn79eva4aGNRZGXZX0GmOm1K665inhnVVHATjRs3zpc9A4vKl1xyCfNWxQ2DJAmVKlWS9bEKSRze4+/9d2zfKXV94zBaiRae0cbuIQ+yiribARdM7CGEiGvXiqW2NKIs+5vJ0aNHacSIESxsD0loENYHZ1GETjkFqXM9WHtB4enXXnuNTeM2btyYpa5DDBB+5Hjusssuk/mRCkmkHvcv83T6RH49U4V8eENqFCQeUZhXjk4UTHRi4EXLOzOijUWR1T5nN5DgoUmTJiyeH9+7W7dudPnllzMHImT6cQrSuz7IK4jUU9gU9iHzjL/DUMYpFfsZC8TpPSPTdgh5QW1bBUbh5qx5QihhgyFDhjCnFThNYq+1MRdZgGlbGfY3a8SJgVFmZibLvNa1a9cCzyM5TrzB9f7kk0/Y+SB7HOor64UtYiYUtZMtJ56jR49mCQi0F1PLvHnz2Bro888/L+ujFRJIS/XPN5x+JvocxIrAiA2pEdBAo9ePRsPOU4AycOlkHYoHomDCBu+++67vOdHGosjKsr8ZjBo1io4dO0ZTp07Vbetbt24d1/NB0hokCkFR+EBx/twb17LTtnAZRq24UGDNE8l8FdbizEn/adsMHUFVGEc7vRct8OhEwg+ne9oCt0n5fblgYo/ODLxoeXo+0caiyMqyf7zBaBNZhUqUKEH3338/WQFMHaMQN+JmMepEjl2IJZI4oOILdAnni7VYxK3LJO7dVdww6sduPfSEMjPNPA/GREZsSI2A3xFyLStQytU8h6E77riDPvvsM2YP5MzGHg21aGNxVIrXyLB/vEEaPEyJXnXVVcxJCFnXkOYUGa7g29KzZ0+qVq1aXM8J1x3OSoj0gEhybcHoHs5MmOHs0KED2/Dvvn37SvvsuHdZcdOg7prCWqSf0hFPNW0b89GKEbB2Bh8DJ9b01OIxac1TFExsGPlgr7WxOCqVZf94s3HjRrZHtMQtZwtrwEEI8f0IP6xfvz6rphVPUKADU8UQTsDFUxzVt23blkV/vPnmm1I/25D1tCr+66+/BlR23DTwtkXvxYyiqYrg5Gb7jzJzMhOnUbZSqkFxtGI0PR+mquzWCMcCd645U6BcMAE6MRh5HThwgE0dijYWRVaW/aMB56BXfzhUdjgxtSqW3pKTk9m6J2Y+kJoVxTaQIg9xl/DEhf9LPMB5oFIXh4soHJuQR13MgAd/G5kYstwHH3zg+xs3C3oB2IJxwQUXsDhQhbVw6QiL26SsLYkunmJDapRDhw5R5cqVyel4dCqtxANRMCGEyPbE0/OJNhZFFsiyPye7gpfSawX+veZsMn59+JohxP+VV16hhx9+2Pcc2nSUJ0NNZ/wN79d4AFHct2+f798Y/QJUwBIHaUhjifJllhFPXpkbNwpiOeHCjOz2eqD8F26yWrVqGflIRYxw5/n33PPyEmfkaaX1JW1DGi1oxNDzx3Sa3UIeEkU8RcEEWP/Ts7EosrCVWcWj4dTZokWLqF4ris/9Og5DeAzi+fPPP1O8QHwpstfBQQh5uDGVjClkhEpihA1xfeedd5jz0M0332wd8WzXrp3vbwzXMbcsPqawD7k60145WbkJ9f2sIjDahjRa8FpM2yrMLd7OBROdmeuuu86XnFy0sSiysuwfb/jAp3jx4roVtFAWEhw5ciRu54R1ZDguLVy4kIkjRp4QTqzF8lq36DhjOhfpYWWSkggJehXG8eqNPHPMcf+PCZLd1GWOVqIFThFwgsD0GdZ+HI1JEwuiYEIIUfiCnY7X62djLrKy7B9vWrVq5QtZyT470hPhJcpkT48GAwJ58OBBv/AV1PjEKBvrtPAEHjp0KDVo0EDqZ0ubdD98+DAtX76c7UWQ0QFVVpo3b86G1JiLVlgPj46/RVaGf6UVu2KWN2YgeENqBDTQSAxvpSlpp8EFE3usCcKLlq8NijbmIsvT88mwf7ypWbMmtWzZkt1vP+tMzfLHuMiaCTQHFb4WLFjAklPIFk6p4onctYilOX36XD7U1NRUFhOEeXCsy/z44490zTXXqKkmK+INr0C2XbFSQLq2IY0WjFxmzJhhuxFMIiEKJrZly5axvdbGosjKsr8ZPPHEE2w/YsSIAiM+dOIw4gMPPvhg3M4HMZ5mLRVKE0/cNE2bNmVDZNEbFyPRXr16+RLGY8jPL7LC2mSbVGA4wWdtCzSkMhyG7NgIJwqiYKITAy9T7LU2FkVWlv3NALWY4d+yfv161t5j2hTOojwxPIp8xzNxB+59VPSytXiiukPdunULPIYSZbiRkIoKw2Ys5GLYH09vLEX0eEyKnYvXmq5ZaKf3ogUenWi8rBKC40REwURDDi9a7LU2FkVWlv3NAv4t77zzDtWrV499JySDRzkyDJbweDxp1qyZaZWFpP3qkLYJXljiNNnvv/9OF198MavxyUGZMjEuR2FdzAo8jwV5FopZFRtSI6jctuYjCiZi3XkdY62NRZGVZX+zwPcbMGAAS3iTnp5OaWlptGLFClMKUyPRPpLzYIs30hZL4HotVgrAl8FFbd++fYHjcJMh5lNhfXItNFozijvXOlObvCE1CkIeatSoYdspwERAFEx4PIvrfaKN0e4hDhKjNln2VxDzqenfvz917tyZ7RGuAscmvSxKAM9ZbuSJQtjI4oApWsyHP/vss+yG0gamIlgVgasK6+O1Uc84lMepx6RiyaFGK0ZAI/zbb78phyET4YKJPToziHXEXmtjUWRl2V9B7HpjChy+NEgRCBFFesA6der4bdplRcuMPJ9++mmaM2cODR8+3NeYwfu2TZs2vmN27drFvG779esn62MVMcRO6flCedOaVSxZD7EhNYpdp/4SBQgl/Dl4jmHu5wG7iDYWR6W4V2XZ3+lcffXVpl3HFJkLt5iqnTx5MiuWirVO1FATQcwNHIZUYnh7YKf4wVDn6vVaZwpaO70XLRi5YFkEnVI1bWsOomCC0qVL69pYFFnYKp7hHInMsmXLTPtsqW568Lj68MMP2c00evRoVi5JZODAgbR27Vpf2iSFtfHkJY54mlB9LyDi9J4R0AjD01AJp7lwwURnBtVEsNfaWBRZWfY3C3TYXC5XwA3x/E5ARVcrguBJnJGnSYnDwxmtRAumBuGAAkcUFa5iDlwwUd0GnZjjx4/77ketjbnIyrK/2dx22226qfic4tMizXp79uyJ6HiZXk+K2GAlJxvDI0+LTUHzhtSoeC5atIiFCCjxNAdRMGEPFINGUhiMwEQbiyIL0ZRhf7MZP368Lxm8E5EmnriI4S7cco8zhbXxeu0z8gyFx0LfRTtaiRY0wsjfqTAPUTDxN1KRotA1OnOijUWRxXStDPsrEkQ8A3k94YZCHBpGpvgbIS0qztMeWEhvQuL1huq4WcezUWxIjYowYtumT59u+ylAuyIKJmyASh4c0caiyMqyv8JcUuLl9bRlyxb2Q0ePDPXXFNbHY7GpTkPTthZa8xQbUiPTreisovSSk0IerOYBLgomQlAGDRpEU6dOZTYRbSyKLJBhf7OZMWMG6wTgOyCnOaIonLQcFzfL4eIiDhRxnqNGjYrXxyoMkEhtssdi4omG1GiMJsIh0Fg7qZanlarj8PNB2Ak/L176SmtjLrI8BlSG/c3m5Zdfprfeeot1Fh555BFWiPqll14ipxDXbg9y3F522WVqncYmWMlDNRSeEOn3rFTYW2xIjYC1MyTntmvIQyKMPEXBRCcGSWKw19pYFFlZ9o8GjA4Rk6/dIl2i+/jjj1mt5oyMDFYx65VXXmHf5/nnn2ex/vECNaQxqxkKxNfiWJkkmXHzawtmKxSxHpFYqRi2drQSLWY2wmZhtSlqUTDRiYEQYa9nYy6ysuwvklI+i4rUSAu4pZSRV14QMfx33303S3dXrFgxNqv4zDPP+KavX3jhBZYuL14xp2PGjAl53NixY1nGO5nE9VeHBAkoR1arVq14fqwiSizYyY/aYchrsZhVWZXtUQxYYR7ccxbFoCGgKMuFPQYJoo25yAKIpiz7RwoErkWLFjF5706dOtEll1zCqq2sXLnSryiImbMRsZixSJHZGwkEqqtgaA1HIXgIItOQwgZYS29CEEI8LfRdxIbUCPgtYaRz9OhRx4Q85FqsTJ4omBgVYzSGPZxoRBuLIgtbybC/FWnQoAETz4MHD5KVOHDggG5CB0uIJ4bquGmCKTzqfSKB/LBhw2R9rCKWJCVOyTErOZpoG9JowWtPnTpFTsJqa56iYKIzg+w6qG0MRBuLIivL/lbk5MmTbF+iRImYfcZHH31U4N/btm3ze4wDm2BNdvHixXT55ZdbUzyRJiwQiOusWrUqc6uP5UVVSMZa7VQI7JOeT2xIjQBvzc8//5x69uxpasgDzuPYgeNUrnLZ2H9WXl6BhBfes/8BOK+Ew9ql62jyw9OpVtPqNPrrJyNeM8a15rYTBRNCmJWVlX+eHk8BG4siK8v+VuPo0aP0yy+/+PKcx4o+ffr41r6xR1k+bME6XKjvCWcmS4qnGVXEFTHGOnojIT0fWQaxITUCGmgUVe7Ro4ep4rlj827q1OxG2kmb4v7ZVatUo4OHDrC/I+6Y7yUaX/zFiD8T1aLgpAL7aQVz7ty51LVrVz8biyILsZZhfzNYsWIFHTlyhNVpFsUflX3gRJSens6+f/Xq1WN2DhBBPsuJ5UKM4rt16xZw4FatWjW2HosBnEyc46anSGxCiKPbQjF12tFKtKARnjdvHplNTkYO1aWmtCFjddw9f3EtI53K7Vb2Xso7u3b6ye63qGLV8mG/FlOu5cqVY6EZsJ0omPgbMY833XQTOyfRxqLI4j1k2N8M4LuCQgRVqlRho8uyZcvS7t27afXq1ez7YQ0exaljCZYIORjBX3vttabkDlDiqUgI8kIU7rZSeTWxITUC1nGHDh1Kr7/+urlTgGdLUWH0G28xiOrz3C5KcZ2Nv8xyR/QeODY7O9v3b1Ew0XFAon6OaGNRZGXZ3wwQp//QQw8xb9o///yTrXFixI/R3+23386eQ/hKvMCI1yyiFk9MO0QLfmgIsFUoZOHJC+4Q5A7xfDwRG1Kj062WqGrhyXfWysnMoSJFipDlEfpRmWeyDNlOFEx0Zp544gkWU4g2TjxOFFkgy/7xpkmTJvTmm2+afRqWIMWOiq9QROoQ5LZQkgSxITXSeGK0ifU3q4SP5FlodB8uEHwjthMFE2AaU+84UWQx1SzD/gqijh07hn0Z0KFZsmSJ+eJp97yMisTCZaMQB+30XrRg7QzrTmvWrDF17YxfWY+FwoHCJdIQJj3bccFEZ+a5557zPS4eJ4qsLPsrKGRBEsCdi2Q7aKlujyIhyLPQyDIUaEgRJG809hSNNZwzTA95OFs0PdTUuRXJzck1ZDsumNhjRIlwPOz1bMxFVpb9FUQ7d+7U3bAsCGFF2kDkF8AMzY4dO6ReMuUwpEgIvCFmQqy05ik2pEZAT7px48amhzzwfouVKteES3ZGriHbQShRo/j3339nU7Djxo1je4x0xOPEUSlEU4b9FRQ01WudOnVYEntM7Xbu3JklSZCZGlbayHPp0qXUvXt3X5CsHshqj2OMZrdHMDRyNPbr148aNWrEAmDh8dWyZUsW94N0gIGAa3Pr1q1ZqibkBb3hhhtY7FIwEICL43A8XofXB8pooTBOTnZ+Ym1MS4a95blZ8HygDZOL4b4XPCH/XPI3bf5zC635eR2dOnZaqlnF0YoR0HBXrFiR7c2Ed0yslAIxViNPre1EwUQnBg00T4QgHieOSmXZ3wocP36cKlWqxL4zSpJZEYgncu7+3//9nzVHntOmTWPz+B9++GHAY+DOvHDhQiZ26BFEy6xZs2jAgAE+7y94raWmpjIRRLwPCs0iAT2MKvLoo4+ycjlwpUbQLBpJnDPO6csvv2TlerR89dVXrIAt1nhxziirhkVnJIVYt24djR8/nhICOEme88BnuCrH/zSyc7Lp3jceoDonqtKY/4WulhCMikWq0tHs/BybS3vMifp9Wl53Mf298C+ShThaMRIXiXVOK5QjSzovf8TpLmL+uYQFBn08q2FVjyHbiYKJ59C25OTksH+Lx4kiK8v+VmD48OF07NgxsjpI2oDc6jKRZrlVq1ZRq1atqFSpUgGPKV26NHNwQIyQEdBoPPDAA0wMIZ4cJCO+8cYbWfUWPAeR5SC3IYSzfPny7KblVQ3wN7L/I/AXe3E65cSJE9S3b182CoKIYtQMUFLtqquuogkTJjCPuXhVD4glVZeVYp4fSHXGpgHRHpowG3g65zQVvTqFHqz/EL30orHCutE6CTzx19N01H0s/1rgIkj2DBAbUiOgQ/f999+ze95Mr82kai6qsqIkeUrZY9q26g+lmHc2i02ta8x2omBCCPFvPm0rHieKrCz7m82SJUvYYAlt8TvvvENWBeXREJOKQZtMpFnv0KFDVKNGjZDH4RijGfcx6sNIVxROgPRLcAEHc+bMYTc057XXXmP7Z599tkA5IPQAUVsPCbZnzJhR4P2mT5/ORrRI/cSFE1SuXNnnmg4BTQTwo3Yl5f+w+d9mrKV5XPkjAU+Sx5eZJdoNqbmieZ032UvJhZMppXAK2yenyJ1eExtSo+I5adIk8z3fk7xn7xuyBUzAkpPYPe72uA3Zjgsm9hBMLF/xTpt4HBdZ7GXZ32xBGjhwIDVt2pRGjBhh2nns2bMn4LZx40a2vIdZxr1799L1119vzZEn1hzDKXKNvIiyewAiWPcEyAKC+XgIKgyNmxogD6gWPIYsLUh1JpYKQq8+0GvQ28f3wIgW07+x/E5Owus6OwUYYaMmkzyK7RqiOFoxmp5PZtxatOSctZXHay2nrHBwR3jOWttxwbzmmmtYJwap41B0GY+Lx4mjUln2N5MXX3yRea9ieczM74AkIaE6IbAFfGMw2rekeGLK9tdff2WKX7NmTd1j8BwcijDaixXcHRkG5YWCUZIGYgrnCr2ExbwCANYwRf75558Cz4tgVNO8eXNWuw75Hi+44IKYfB+nwRszD5nXECfFeL5abEiNkJVzmkY+PYrGjp1gqvOJ252/1pntjizhgBUQK7REYztRMPEcr+6BBls8ThRZWfY3C7STEyZMYEtdbdu2NTVhDvxQAoknr+bVrl076tWrl3WnbbE2CIHCGiAERQseQyZ+ODjg2FiBdU2AITpPFQbRBoEy/WPUjLVO5GnktfgwXXv69Omgr+OPIzGy3WlR5BiVScqiMpRFZV1Z7O9riu+L+3lUTE6lBiknqHqSeU4Ivc/bevZaZFJpVxZdXyL//pGFOL1nhOMHriTyzCazaV7kANVITqUahexRW/TyYgepNGVRxaR0al38oCHbccHEHr4R8KLFXnscF1nsZdnfDHD+/fv3Z+0lX7oyE8Ry/vTTT7rbggULWHQFRD4WM4PSRp5Q9q+//pp5rSJ5MKZP69Wrx55DwCpGcbhZbr31VrrnnnsoFsyfP5+tW2LU+dJL55xNeOgKgmUDAQHFuifEE05PYrhLoNfxEkhccEOBigN64Prwa2UWwyrHv5yUHoU8J+iZahsEd8j407rkMWpdcmvM3l8crUQ7+vDk7qeqRXNo3JOlcYcTURkyi1J0kkaf/w+RN3CImJUYWHHbuX8ktzdkOwgl/s0rfXDvZ+1x4qgUo06j9o8WRBToCcmGDfjNheaNN95gzjfvv/8+c750MlL9pGfPnk2vvvoqc85BlXRsHPRUHnvsMZbxIRZs3ryZ1ZPjXm587VNhM3LyazOS+5CJJxHbNU+xIY2aM9PP/Z02m6jMQDINz9nk6l77TdtShKM/PdtxwcTUuRhLKB4niqwU+2uoV+441T0/8AzJP6VOk9G5HMzgweES06B9+vQhs8Cofd++fSx2Ew5LwYDTEGY94aiKzoplxRO9KFzcJ598kp0wPJwATvziiy9mc9CxYP/+/WyaFtOuw4YNo0ceeaTA80hsEKrSPIq4Ah5qw1/DX4cwm1CvCUWg3l2gEakjyTs7KnD/Z+JJxNZ7VWxIo16rzPn63N+Zn5krnhg8YYn6rKe0vfAYsp0omBhRorYknAixDqe1MRdZKfaPEniftmjRIqrXDho0iDk5vf3222QW0BQ4a0JTUEM0FDgOs50Q261bt7LC2LKISYQupk3hFBRLxyAxFhOuyFh3xNy2XtIC7sCEC6gHRBBTtihyy4UQYlmmTBm27onX6fVw+PvJTPnkePL2518Cj5kjz9hPpRlPbiB2BI+TuXDnLvt520YzWhZtJwomBg/IocrjPMXjRJHl2bPsxnfffcdmEBHaJ8Jrk2IQw2PeMQuJgtmyQfggBBzrreEMWnAMZiIxXY0lPTFxv1Hs6e51FqxLdunShQ3NEYeJJNl6nldwU4bz0NGjR5mBtaAqBdB6zPKpX/68CG7+f//9l60fNGzYUOK3cjo8FZ6ZziexHUHxhlTeqMPkwspe/vn2E4RIba21nSiY2HjCCu1xEFkIC0/PJ9f+8QODDISn/CxsPOkNRJQ/Fqti38gIh6gJvWxwgUAGOsTmy84wZFvxhGcvkhcgsxGS/iIlX6CbETFVvO7bF1984fc8nJwAvIFF8EMQn9f2wnCDoNepYjxlkm2Bhji24ik2pHLwWmPd0IZxnpGub2ttJwomOtRog7DXO46LrHz7xweMpvW2nTt3sufh9Mgfi1WRdvi2IEdwpGB9FCGL5HTxxLQHvHuxcIw4I2QTCrWeirVQ8PLLL7O5bw7S8yFbEaYjkGheBC7ZmL799ttv2WeIiR5Qmw+ISRUUMuDTkZkmXs7YNmpiQ5oQeM52dOxYYivCOE+t7UTBhCMQRmY8l614nCiyCWf/OIIlNiynRQpeE6xgSDTYMivxlClTWFgMQKL2hx9+WPc4rH/ieYARIhyJEAeKBPXXXXcdmzvHNAB6SnC91pYJQpKF9957j3r27MmyDKG3CPdsrG/gRwJBToS8ttactrVH2EM08IY0YXBl5g9+PYEd8hLVdqJgoh1BAgH4emiPg7hihgqzVfAJSSj7xxH4pYSTyU4LXoPXktPFE161HC6iesCbjYsnQB5QCCfEF6KJ0SpEFYvIbdq00X2P2267jZVQw4j1jz/+YIIL56HBgwezHLsK2ZyygHjGNsOQtiG1PV4+UrfXNGQ056y1nSiYmBGDRyriIPG3eJwosgln/ziCthftMFKuYjkuHBAtgRlGlJKUiS3nDSCKgebfxU1v3h3xSQijwfAfIoxF5EDCybnyyivZcTger8OPQwlnrODej2aGPcR2DVFsSBMCHqJiYj7iqIkiPZ9oOy6Y2EMIEduuFUv2MWdFFvtEs3/t2rXZ99q2TUg+ESOQwQ5tMAYz4YJjIbZanxZHiqdCEVtiK9xiQ5oQ8K9hx+/jMmY7UTAhoPCi5en5xONEkU04+8cRhMnAcxbXGaIYrKIQnkOmORyL16AKjExi0vWBFxkqmsAjNhCBkscrFImOOL2XEE4j3rPOXV4bhqpE6CGstR3+jThCOAABHgqnPY6LLMC0bULZP44gVSpqK2O5bdSoUSw88fbbb2fFOxDCAhCSiPBCRFYgFh9T5HhNsPSsposnHGn42mCwIGBedV2hcOLIU2xIE4Ozv3VXXsLbWs92XDDhdYvcr/xv8ThRZBPP/vEFy2wrVqxgOdKRtW3ixIl+x/BRPbK3zZw5MybpWqWJJ+IekQYJNwm8murUqRN22jqFwlrwfHOxQWxI7Rgo70+ejXPbegzZThRMDAgQe47QNgwQtDYWR6WJZf/4A8fP9evX048//sjqLqMzgtlOgIgIPA+PZtkFsGMiniiOijlm9ALgiapuCoV9if1alF6mq3AJts5j8qIn2Q+PIduJgokpWDgk8vR84nGiyEI8jdjfbFD449dff2XihZh3eA0jFR8SxqNDEG3u3GiAOMZSIIMhbcIdw2e4ayOWUgmnwt7EVgR4Qxr174QnJbAM/Hxs6G2bZMx2omBiw/obT88nHgeRxUiIp+czZH+TQeWsH374gcXBo7g3vhfWFT/++GNWAASzkE5AmniiColyAlIkBrEVT7EhjQ6rTY+ezWPqNjMrVJRE6PGqtZ0omPDzwFIVT8+nPY6LrHH7mwtG2SdPnmQ5bZF5DRtS302dOpV9d2Rms+t3M0U84f2E+EmFwv7ENkmC2JBGh9Wmbc82lEl5jkjPJ9pOFEzEbSL1J0/PJx4niqxx+5sL4t6L6hTURqY35LdFNh8U60h0pFlvzJgxlJqaymp5OqHXoUhkYjvyFBvSqNAJCfGYmleWJ4a345qnMduJggkOHTqke5wosobtb2EKnc2YFKvazVZCmsMQcsOiPBjyySKmBjlfq1evrnuDwBNNZl01hcJOoAHFehFq0UaXni07QIo8s9bQ7Oww5DVsOwgmgvAxaIAo8lGXeJwossbtb00+/vhjNn3boEEDtiU6KTJT5kEU4WW2Y8cOtgVCiafCyWhHK3LwWCBUxWqOTPJn6LW2EwUTQsgrNqEd1NqYi2xs7B8eqIOpN+UKh89IGTduHHsd0uVt2rSJ/V2tWrWg5SETCakjT4U9MXfKz5nwhjQqXHqjJTNHfW5NgvjEXfPU2k4UTISgvPnmm2ztDzNu4nGiyEI0Ddlfh4tK7qJ25VMDPp9d/BD9Q3JZsGABLVmyxPfvWrVq0UcffcQ8bp2ANPFUidJtjAPWqqyEdrRif3udFSCX2xHetqLtRMHEaBPB+thrjxNFFtO2huxvgG+++UZaHObixYvZHknuEfM5evRoFuuJLHMjR46kRCfxVqwVCdAYJza8IY264fTqrXmaOWV69rO9Z0NW7ITuKD5824mCiRHljBkz2F57HEQWsZ28+ooh+1uMsmXLUtu2bWn+/Pls1Al/FuTtTXSUeCqwCK2uQhwRG9Ko0BVKM6dMcy0afyq/46i1nSiYfLSJvfY4UWQN29+iFCpUiO644w72HefNm0eJTtTTth07dmSOPx9++CHzqsW/wwWvE+fKFQonITakZo2gYoMNxSDJZch2EMpevXoxJxmscyL5AfZIoSgex0WWv0a6/S1ChQoVfJVNEp2oxXPZsmVMBFGlm/87XPA6hZWwWtB9YiM2pNHhtZhw8VGvDR2GIswTrLWdKJg8+QF/XDxOFFnj9rcuP//8M9sjWUKiE7V47ty5k+3PP//8Av9W2JHE6wFbGW1DKmXa1tSKJh4b30deQ7YTBROOQDVq1KC9e/eyAYL2OC6yhu1vIr/99hudOXOGOnXqVCCGH9/97bffZrGexYoVY9O3iU7UloNbcrB/K2xEAk4fWRmxIY0OvRGemTZ02XfkGWExbK3tRMGEEEJcsMeUrHicKLLAmP3NA45O999/P5uehXMQyn8dO3aMedsePHiQxZB+8MEH7JokOvbq9ihig0Wn0RF/mpSAwdbahlQOZtrQhmudHI/bkO1EwWRvd3YaWHucKLJwrJFv//iAUJRnnnmGTc+uW7eOCSdS8dWuXZt69OhBQ4cOpfr165MTUOKpsO50GwtgTzzx1DakUnTS1P6Py8Zr517DtuOCielYpCXdtWsX+7d4nCiyhu1vInXq1KFXXnnF7NOwt3jWrVs36g/FesD27dujfr1CNlZt9Kx6XsbQjlbkZMUx81p57GuvKByGRNuJggkh5EWuMW2rtTEXWcP2V1iCqK2Hm4Xnso0U5W0bOzyeHFq5aQS1rPs8FUopE95r3B6i3HN2zM52U2ZOMhUt5KFCpfO9qYPBi/9q/44G3E6YRWZ7k0bEmMgTv0EszoI3pFG+2nedAPs7wulHmbCU9F4ij8segeMFrl0Ugi/aThRMPI40pVgTRBsnHieKLETTmP0VVsDwvY5F44kTJ7JE8FgwDmc7cOCAnLNX+PHF+r/o7pFrqUTximwtIpytaLFSVLTmNt9W4+LDVKXZNirbcAeVKFEi5AbvOv5eSDxtpGH460hFOpxelNYeq0DZJunBjI0X0JHMInQovSgdyShKC/fUlPr+vCGNtnTfupPlaXtqaTqSUZgOZxShTSfL0b60EmQWL69uQ8eyCtPLay8jO/Dn0UrsumEbu+FiQ7bDvY6wE+yxLVq0iO21x3GRxd6o/RU2H3nOnj2bPvnkE/rxxx9p2LBh9Pzzz1P37t3p7rvv9iVQUMSfE1ku8rTvSj3vGEgze0TnLp6dnU3Npk6i2qXK0A99+kU88ozWixCNTq+fu/v+/W9Dc+6hsf9ewTaRHa3lvb84WomGrLxC1H1BrwKPLeht3trwnJ31aeb2ulS7ZCmyA72X3er7u02V8wzZjgsmcntjRIl2kaM9jo9KjdpfYfORZ8+ePenbb79lI0mkmmrWrBnLNoT4HyyGP/744yyLhiK+eD0eckG8kpN8Cakj3UqWLEnJxYqRt3AhKl68eNANrunoREE8scFNPdqRp/Z1NvbhDIo4WokGt8drqXAj99nv4bFhyJMnwkl5re24YPL0fPCixV57nDgqNWp/s0BCHCSW79evHzVq1Ij99jHz1LJlS5YUPi0tjZyE4WlbFHVFRYEVK1YwJyDU9UTF9AkTJrAp3ebNm9OYMWOYZ5ki9vgaML0GNkJywijXpG0U+N+K8K5ZNHh07IL1RrPw2rizE414irYTBROzbZdeeqlvvVM8ThRZu/5OZs2aRbfeeiu99957rKPctWtXlhAeCXJGjRrFvvuRI0fIKSTJdmNGRn0URl25ciUNGTKEjh8/zuKCcGEVscfLHSEkNKZ5cCQKgdgoiH9Hg7YZs+NIJhyMXyd/47pMjFXhn2zHXK1JriRDthMFE4IyaNAgttceJ4qsUfubBWalHnjgAVZKbePGjfT555+zZbv//vuPWrVqRZs3b6ZHH32UnELMnOOQcQjhLKgsjh+V3XpZtuWsanoljDyTw+iVi42C+Hc0uDWNL6agExGj10lPo9wRZsqRifesf4PeiNjqRNpB09pOFEzEbyKnK/ba40SRNWp/s8C67rRp06hJkyYFHq9atSpNnTqV/T1nzhzKybFhdZ0oSJE9J46LhzUwVE1ByZ0yZcrQgAED6J577pH5UY5n98njVK1kab/rUAjhQ243FRaS9p/MyqJKpc45c4QbThKO/oqNAhD/jhTt6+w3jgkP7TWL/PU6D5p4sfiI0yOhwxZ3IhRPre0glI0bN2ajLowksSbIp2bF47jIArSLRuxvRVq2bOlzNsRsIwQ10TEsnrhJFixYQDNnzqS5c+eyBhshC5gPh+ftDTfcwP6tkMf24yvogSfupaXvBU404a5Zk0r0ulv3OaxPYG06GOWKnKKLqoWO8RQbBSD+HSlJ5KHCrkzK8RamQq48KppizmiqVpmDtPt0Jd+/SxdOl/r+2msWKZVLnqAkyiUPmzhyUTLlUZlioW0VK1pU3UNrDp5PF9fcQVYHQl/UlUlZ3vw26ZLq+dmAorWdKJjcB0TvOFFkMf1pxP5WZMeOfNvju/FrkOhELZ5Y08QI87PPPmO129CLuvrqq5lgIschRpyK2JCUdIoefLIkfTflDKWkFNE9BlNCaCjWHXqLdmV+RLc1/Mv3I4YT18iRI4OmBht57XeUQmVDnovYKACxgYgUJEV4pctXviQJSUnmNDCDr8yvNSsmISB6Udr7axvSSClZOJ3G3PBZgfMrmnILmUXvi1bQnZ48KpRUnqwORn0vCfdYuaRWhm3HxQK/OUQd8FqW4nGiyBq1vxEQhw0vWS0bNmww9L6TJ09m++uvv56KFNFvkxKNqMXziiuuYILZokULGjFiBPXu3dtXnkwRW7CemZzsOvuj1P/x8R9l0aIuKuzND1sB+PEixCgcZ4Vwxn3anrf4dzRwQchv3MyZBtTOpsmeXdNes4hJ8hQ4J3atTE2N5yVER9llErLAPebyGLKdKJj4jZ06dYo9jntXa2Musobtr8PFRbLplhKBZ0jWF47dOuT8+fNZ6A2+/0svvUROwbD14HWFUQy2cIHoYm5cESWufFEJR1w8Oo4k4U+rhNewiO9nZMrG6+dwYpfmOHIMTW3pOFJ5TXQY4nby2mCVWvubiWaZVrSdKJhYwoIHKmLg0caJx4kiC9E0a2oToo0Bjyw2b97MZhtxXceNG+db+3QChrxtccG4l2UkG6YtFAau+9mgvnBGRG5XwWvNf8ThefqF/gDx/SJ7b51PO9spOPdvSkiMXic9s7hMzSprn8TwWicdV4TnrLUdBBNhJzzOGUtZPD2feBwXWZ6ez5D9LcL+/fvZNO3JkydZlrlHHnmEnETUvzh+s0S7KaLH5eKesqGvozb+T/wRh35taMT3i+S9wyMx7xOj10lvlGmFGEs7jDy1JEXYBGptJwomRpTz5s1je73juMjK/53EnxMnTrBscrt372YpB8ePH09Oww5FEBRaksJvpDSDuQI/4pCvDeP9tT3vcN/byRi9TrqvssAw3Y7i6fbmGbKdKJhwRoIXLfba40SRtfvvBGn4unTpwpbskM/83XffTaiwm3BR4mlRsnPOTXFrt7zcPMrOdFN6ejqlnknTPQY/YODWlKoSf8RB8YTnMKRtFMJ67wCw0RNvf73WGE0xJJ+G0euUlMQ8hPJPDXtWvs3c5HjMO9oizQnum89+G0Ur//vC/7eTl0N5uV625eZ4qFhyFd3fT1p6pu83FMx2omCC2rVr6x4niqxR+5sJfFW6detGq1atos6dO9Onn35qqAShnbHG3a7wo/Mz06h993t0S4g1qtyVhl2/m8qWrURlSpfSPeaaa65hP+iVfzamtSvOFS4Xf8TB+PfvGrTkx9DlmsT3C/e9A5HrTqbdOyvQqRNFaO+u8pSZaU588IrlTej0yaJ06nhRtt+8Xq4XudHrtGNHZTp6uBSdPlmEUk8WoUP7y9Khw6HDimLF0h8vobTTRWjh99ZIwXnsdBq9MHwxXd64p9/volixEnRHk3/YNnGAl9rWmaD7+ylVsjjdM+DhsGzHBRMigoIYPD2feJwoskbtbxY49169etHSpUtZTlskxHFyDL+9rOcgTnhzqGHb6ynnq490n8eo9PqX3qWTaVm0esJjfjc5j7s6ml2cvvn7KnrpnnPPDR06lF5//fWgPcYv/uhApQqHEc4ivB8I570D4nLRez90yR9VuYjukVgGLBIWrL+Ebfw8wDPnqlgZJlwbBCI9L4WmfFUwrvPyweb1/pfvaEQ/b2tEZYsVtkzKveTLbqVml91Ca8c/GnSEGmi6sdXjk+mGrleFtB0XTIAR60UXXURr1qxhWby0NuYia9T+ZjFlyhT6+uuv2d8VKlRgBUH0wPonnk90lHhaFLTbbsqP5dQDj7sKFyZXIbffMfj3k08+yf4ukuxvYv4jjnjBNADi+4X93oG8iNGYuUyeGPEFAsbuI4xcJ4Y2laGpqfFczFRmVnYpgMtFrpR8UYrWKceVnEw5Ab6QaDtRMCGEWP/DHsIsHieKLMTTsP1NAF61nK/PiqgeyF7mBPFU07YWJmTsnkffeQQ/aMRyYa99XpxaCkkYraH4fhG9t97HaWPwIgxgtwtGr5M2gT7DROGKd1WVnLw8JkB6a5V8Q35nTwCfAb5lZmXrrmty8nQS3WttJwomRrHIHMQrrIjHib9Jo/Y3C4gibBxqq23DjkE0KPG0KkW95C4Z3BPQXTKXyD/TVoEfdF7xvALZZ8QfcTA8+K9UaCcU8f1Cvffq/bvpWNrpgI1ZRnoGeXKyyZ2TTbkZ6ZSZlk5nzpxh+ZL5lpWV5ecUJZ0UVHkU/5Mr4uHaICDFCp4dO8NCJnY0iuZ72rpKxeccLnpvLF1y1+26a5XYalQqTxsmP04Zi2cGPAZb8WJF6fI+d+p+hofclFc4N6TtRMHEUkrFihV9jn7iceJv0rD9FZZATdtalNKHTlHZ1BDHbDpFnmT/H6D4gy6xPYsKl8jPgKL9EQcjpdQpKrMr9JSX9v2CvfcTA2bQicOr6O81C0K+b6HChanm288FPQZT0//73/+ku8kXyz5F6RckERyVUe6x6H65ohCuDQJRdH8uueqfpjwsMXqJCmUTpRzPI6pJplD6WCoda+2hcqviE7dYZHsmVStzEa3K+SzoccHWNEGHAS9SlZRz1YZEirtPU4l92SFtxwUTHTtMEXNBhCeteJx2VGrE/gproMTTopT600NFagQfWRXbmEvuTP+GXfxBpxzOpkpz3USjye9HHIyqX7qpePHQ03Da9wv23oX25dEFVdrRqpx5fs+tWPw3vTz8c/r2r/wE7LwREsunoTHESBSpzfAdS5cuzXJpiutaPMYuGEf2HaNylQMXLii9OJtKL85PF5iXk0HuZKJ9fff5zkNb0k38dzjP6dkgKyOL1v+ykS7sqJ86TfxeZUsVpcpTM/Ov0/795C1Tkgp/S3T69OmIz4mz+LPfaNKIT+jrHfmOZsFY9NkK+nTS9zRteX5lnhvK1qOFY1dRp35X0uHDh8O6NuGcUyAqfu+m4mX91/ojpfQaomIV9X9j5Re7KaWR/32ktZ1WML///nu68cYb/Y4r8JtMSQnrN6iwNko8wyAzM5ONcBDYvGfPHtZ4Iy0VGu5YJcN3b/iP3O5qwY/ZsZs8J/xLUYk/6KRTqeTeuNX3nLanHPC9N20j93nFQp6n+H4g2Hu7N2wlV05l3eeSczxEm3dR8eLFfY8dOXSU2tRqT/9lrPNNd9WoUYO9P47Tpjd7c9h7NGnKJNqau55kU6PGW4bfo0OHDrRo0SLWyGqv09dTvqeRTz5HO2lT5G98iqj5he+Effidd97J4gwhWJw/vvmDNu5bFlHowcwS43x/QwiWPfkGPZLvpxb1tQlHQHE/e6sYr9qUt2kLuVMrBPwMz3H/PK3a348omPh70qRJLIEAOnriceJvEvtwfoNWZPXq1cxOq1atYhtS9AHLxGTHEbXmGQKssXXs2JEJJTJrIEAYDfj7779PrVq18tWxkw4Lfg9xQ7JEAv4P40eMODI2UtE8z3/E/EerdbxIPZlKW//eTnmeXCpcMpn9wE8eP0mbV2/VTcIgvp/2vfXO1+3W/04enV748SMnyZNz7njx/Qt8x7Mc3HmIauY2ZJXsxS3tTBpLKIFtz6591I660d7de1mKNGx//byW7fkxfEtNTaV9+/bR9u3b6bLLLmMdJ2yoKHTw4EGW5Bv7K6+8ko4dO8beA68R3wMdL34eXBz0rpM7z0t1qanfuWODDYYMGeL33tiupq50d+MBvu/CN34sRqPt2rVj78HfC6ND7egcVz/Q52u3Ede/QB3o1gLngfc8cOAAHTp0iF0XvqEwMs4H1wexgTgf7XURr03YSGismYNLkOQSOv5CfrbjgsnT8y1ZssSXnk97HL9fQ/5OLAzawaeffpp523LhdCpKPEPw8ssv0x9//MEazC1btrD6pahlOmHCBNY49O3bN2bG8QQQmlDhCeIP2qMJNxF/xNhQEUF0oihzXhlq2Ko+LaNvaOGJr6hEiRJ0XoXzqMklDQscx0sPad9PK2h+5xyo0cvzb8RSkpKpnquZr1HVfhb/jueuV/4IiAs5tg+em00ty7Rm3wNbzdrVKa/NMapRqwaVLVuWbZe0a0W9u9zHRrN8Q83DZcuWUdWqValu3brsHkCnCduKFSuoSpUqzB0f+19//ZXKly9PpUqVYq/Ba/l7cFHAhtJNga4TRjTac+fbc13H0LKZvxc4vzUL19NHz31O1etUoatvuYqGXPoMlSxZktXRFc8DU9v4u1ixYuy9sEeQOxr4TSu30u7Ne1hDXq56Wb/0ehiZ6p1PsiuJkl3JBc4HI6mtv++ku2/sQ+XKlWPXBjM0v/32GzsfXJ/ly5ez8xGvC3s/4drk5uQyMQ3mJevxesib7Al6TLgOZUHv1WT/Dp3WdqJgTn/6E2rToJ1vil08rsBvMozfiVVBO/jcc8/R3LlzWcfRKbU79VDiGQT8iBEYDKZOncoaJw6qCFxwwQX0888/s6mMmKDX9SWiWf/7ivXcK9YuT0VL5fdit/69w1fmTfxBuzUCq/0RY8TARyXY7qjZ3zeqQM8S+76thrLH+DEDLx1OJ46f8DUS4vtpBc3vK3nDL53m1ghqoF6+73m3/+eePnmGanua+M4d1+3S1pcWGPl0pO50TbNrA34Wjnn0kcfY98WGMAPeMC+ZvZweGvCw37XQu9bBrlOwqJw9G/dR1bQ6BUZmM0Z+SuMmjaNZO9+hD5ZPo4+3vlmgQ8Pff/4Hi+nCKq1957vtn500fPhw9u9h7Z6lzhfezDpDg8b3pzW0vEAHqX///rq2DGTDGc/MpO2rd+tew9STZ1jd31DXqV/zx9hUfTAv2aU0h3ZX2Rj0GGxo5LWC2qv2A3T6dOq5FJS5gQU2Secp7bnf33QI3dX9bva9vpnyPR3adlT3OPF+Ded3YlXgpDd69Gi6+eabWcfRySjxDAJ6zZhmqlevHpui1dKjRw+2Ry9SFptXb6NPx3xFxcsWpTqtatPrQ95hjbfI+yNn06C+Q2j62kl0uuU+1lA0bFWP+vcayJ4XG/gkTQCg+CPGhlEIH5WwxAtnR0kYTaCzgH2RQkV8owRs2/7cRS6vS/f9tIKmBWlZ9fDk+rfI/DP0zl0rYm8Ne58aXlSPSpQtQnPf/pH+XPz32c9zFRjR3VquD53Yd64KDLYkuNS6An/W49eNpgXv/uQLFP964g++48be+wb9NX+d32uCXRvd6xQsIQU7Pc1P1eum+kn506xffvkl6wBAFLWfu/ST5ZR6OM33sgcvfJxOHcgXj7xcD9VNaszeY9Lgt+kiurrAVOr06dMLrIue+2z90/Tkuahe0rmZAvE87q0/iFbM+SvkdTp+4CTVzK1X4Dy6VbiXGhQu6Ej15+pVvr8DrdNiXVcrqLN3v0svvvjiua/ijawAg/bc9246SPu3HMn//h6Pb6ZEe9yapetoyMNDfOn5Qv1OFNZHiWcQ/vnnH7ZHBhE9+OPr1uU3njKY0PdNevaZ52nuyZm07PD39MiUgfTEE/leGAMvHE7j+k2lZm0b0jUXdaI7a/enxYsX+0ZPd97Umx13f5OhtOSjX9jfKRq1EkXnl6//oFuuue1cEnm3m0qWLUZo+8RpLy+mj4V2BGI0dMgjrJEQ32/+jIXUtMyFBabL1v+6gR64ZDgb3ba76wrq+2ov+nflRrqv8WDfMa8Pejc/pVoRF/06dyV9PnFu/udqxHNo22eoa4fu7P3H3T+F5k5Z6HtuzqT59NfhFXSyxV7q9lAX6nf7APa4RzPyzst207UX3sDO/cf3ltLXb/5ARUsXpRoXnE/LvvhV9zqdPHiKarkbsdekpWVQHU+Tcx7AHqKLq1zB/j1nyjy6qNplvu8/ru9UevSRR9m/h7R5hto26aA7ev31m5VUuW4VKlIihTb9uYX+/X0ze/z0iVSaMOBNaty6AVVtWImmPfURbVmzPd8GHhd5hU4BOgD8nMT3x/cXp75hu9633pvv9ZvkoqTk/PU5iDPew/d+SUn01FNPFbg3WGKBzExqdmVDKla6MPv3iu9W0oEdh86+d8H7RDyPjDNZ1KhQC/a5m//cQp0u7VLgvflxrA+jOY/azWtQx1bX0WsPvknb1m9n93uf5kOoUbGWbG13145ddGvn7gXWhLHW2qtVf7/1VfxOBg8aXDDJPhFtXbuDJg96p8DUcOUG+Unj//19E/23epvfufLreWWd9ux7ibea9rinOo2mlfPWsr/X/rSeurS5OXZxyoq4oLo+QYBzCKhevbru8/xx1LSTxWPTBlLbRZdTr6duZeuDU4fMoEFj+7HnOvXtSIvXfE9vL3+LLs+6girXrObrdd/22E3Uvkd+Ls57nu9JxUrkN5B3PNWd9mw8UOAzsPYHSpYtQWXKlvY9jmk/jGbbXHEle99Ro0axrCI4p1U/5P/wwc0PdaIadar7vV/LDs2pQct6BT7rq8Vf0A+HPqd3S7zG/t1ubzs21X3jFedys5avVo7m/zGXMlsfprbdLqd7brmPej7Wleo0q0HX9+/oO+6yLhcR7U5nf/cc0Y1qNKnuE4WugzpTvzF3UUpKMn3/ziI6v37+lNK9z/agshXOfcfuj91INw/sxP4uU6k0ffbx2/Td6Zl0avYuKp9SmdrffpXf9xr+7kO07e+d7O/zz69KNww4N8Xb95XevvCSC65uTrXn1PI9V65KWTrlyp/q79SnPWUtP+733mDXv3tp/em/6MxFB6hp60bUp2c/ev+K6VSoUAr9smUpfbv8K7r66qvpgzFv0GMZw+i1iybQ8PcG0aY//mOvx7pt5/s7sLVE7fsPeX0A/TQrvyMFuvTvSFd0yU/4P2DMXXRelXLs795PdafkpIIOO+I54t7AiK1WrVrsfsf58HvvscHD6LU3JtDwGYPon2Ubdd9jwP/dTU0ua8j+Ll6mOJUpVzCRve983+jH7kvxcycve4V93rS3X6ODxZ+gsRPG0EPj76cee25m66mob+vOIDZLwu+HsWPH0i+HF1GJEs/77mN2nzzUmWrVzg+I7fPyndTw4vyiCcXLFKPlm5fQo4XzZ29AzhtHqPOd+R2q/nc9QO/OnOZ3Xa7v14HueCL/Xn7szYG0a9Ne3et359O30ZW3XMrOr2S5ElS2nHFvYYW5uLxO9DEOkwceeIAFM48cOZI5DmnZtm0bNWjQgG1wJgoFqsfrAW9OTA1v2LAh5HvwNSOYjTuiyEJ8b/SoI4m9C/ae6PHjPfF+2PTem4/IZH1uJOcXq+sZ6Xnoff9Y2CTa88N5cMcmbkcQy3OK9vtHcz9pY4TFJAtmXfdQfPvtt3TLLbdQtSrJVLZ04PM7csxNx064mfOb2MnihNP26FG0aFHma+FEGVEjTwuA6SmIZzjEsoGPxXvj/bCmasZnh4OZghnOeVjt/OIdXhHt94/mdVa51pGAkXfTpk1DHlepCtGZ9O0qMYNElHgGgXvX8gQAWrCewm/gcIi2d6dQKBR6IAZdtSvmoByGglCzZv7aCALl9eCPYx1IoVAoFM5BiWcQWrbMT8+FWn168McR76lQKBQK56DEMwhIu4aMLXDo+fvv/LhBEcTXAQQMKxQKhcI5KPEMAlzxBw/OjwkbNGiQb40TvPbaayy+E3lDL7443/VfoVAoFM5AhaqEkRi+ffv2LJ8t8pwiuTXi3PBv5PNEzlO4fysUCkWigwoyPA0kQGUVhKmgcAIHaRFRZSbRUd62YcQx/fTTT6wk2axZs+ibb75hCa/79OnDbqJACRQUCoUi0UAxDAwctIiP4RgnoEaeCoVCoVBEiFrzVCgUCoUiQpR4KhQKhUIRIUo8FQqFQqGIECWeCoVCoVBEiBJPhUKhUCgiRImnQqFQKBQRosRToVAoFIoIUeKpUCgUCkWEKPFUKBQKhSJClHgqFAqFQhEhSjwVCoVCoYgQJZ4KhUKhUESIEk+FQqFQKCJEiadCoVAoFBGixFOhUCgUighR4qlQKBQKRYQo8VQoFAqFIkKUeCoUCoVCQZHx/5VU8rnDZbK3AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "c = PLSCAN().fit(data)\n", "c.leaf_tree_.plot(leaf_separation=0.15)\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "5994c87d", "metadata": {}, "source": [ "Using sample weights below 1, the leaf tree with weights has fewer leaf-clusters." ] }, { "cell_type": "code", "execution_count": 3, "id": "de849fe7", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc8AAADuCAYAAABf/4wMAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAXEgAAFxIBZ5/SUgAAXMdJREFUeJztnQecE9X2x0+yu2yjd5BeFLAgYkdEsIKKTwQF21OsT6wI+p4NRfGvAooiNkRR7AULqNgQVEBAUOkdFlgWlt1le3ry//zuMks2mZlMkum53/diliR3MzuZ3HPPuef8jiMUCoWIw+FwOByOYpzKX8rhcDgcDocbTw6Hw+FwEoB7nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3nhwOh8PhxAk3niZg6NCh7MbhcDjJMmPGDGrpOIo2bdrET6aGOHhLMuM59thj2f369evJLjxz/Uv0xpzXaCdtNPpQOBGkpaVRIBDg50VnHn30UXriiSfI4XBo9h779++no9q0oyAFaMSIEfTJJ59o9l6pDjeeJsCOxvO67ndQwbZCClHy7WIHXd+fHnhzjCrHZUd8Ph/l5uZSaWkpZWRkxHy92+2mYDBIVuOLV+fTnMfm0ns7Z1CjJg0NO9ctWrQQPdc4p0uWLKFzzjmHnE6n6Gfk9XoVfUaJcueddzLPM4tyyE3VtGrVKjrppJNU+/34fT/++COtWLGC3fLz89njsdpCz549m1555RXasGED1atXj04//XR65JFH6MwzzySrkk4W5fnnn6fff/+d1q5dS4WFhWxCaN26NQ0YMIDGjx9Pxx9/vGofIr4QkyZNoj/++INd/L169WIX6fXXX6/hX2htsupnsRU2/pcsLVs11XTCsTo4N/H0tM/JySEr0qxRc8p0ZFKDnAbUqFEjQ44B59nlclFmZmaUB+n3+9lcsmzZMkpPrzu14t8YF/m4muzcuZNemfEKNaSm1JNOouX0E13QdwgVhfar9h5PPvkkffXVV3GNuffee+nFF1+k7OxsuuCCC9hcDQP8ww8/0GeffUb/+te/yIpYds/z6aefpu+++46aNm1K5557Ll188cWUlZVFc+bMob59+9L8+fNFP8Qbb7yR1q1bR+eddx6deuqp7EM8++yz6csvvxR9n88//5wZ5AULFtAJJ5xAF110EW3dupX+/e9/07hx43T4S62Jp8Kt2u+qrvSo9rvsCDweXJ9W9CbjIXD47/P7/IYdA4wnFtFiixUYxpUrV4oaSLlxavH444+zSE83Oo4aOBpTa2pPxXSAfv31V9Xe44wzzmDh56+//poKCgrYIkKOn376iRnOZs2a0T///MPmWVyrOCZsH2A+hhdvRSxrPLH6OXToEC1fvpzmzp3Lbps3b2YhC4RIbr75ZrYSTOZDLCkpodGjR7P9IayQFi1axO6xEd+tWzeaOnUqe4wTTf3GDVQ7LQ2a5PJTLAOMJiIjdjee+dsK2P2BvIOGHQPOMSJbYuca8wQ8M7H9ZLlxaoAtn3fffZeaUktq6mjJHutCx7LIz9ABw1Qz2g8++CBNnDiRLr30UhbpUxIhBPDIu3fvXscI33777WzOnTVrFlkRyxrPfv36MU8zkjvuuIO6du1KBw4cYKHZZD7EN998k8rLy+myyy6jYcOG1T7eqlUreu6559jPMKCcaIryi1U7LYW7jZssrQA8nd9++03TkKAZWLu45vu8cdkWw45BzrsEUl5UrHHJAm8QdKXjah/LcdSnttSZyqiYvv32W9Ibl8tFCxcuZD8PHz486nnhsXnz5pEVsazxlEPYH8OeZjIf4jfffCM5RggTw6NFDJ9Tl8ryKtVOSWlRBT+9MsDTefjhh22fQRvw1Px9lRXVxh2DjHeJCBYW07iPZ1yywCh/8cUX1ILaUiNH0zrPdaae5CQnjbzkGt0jE5s3byaPx8MSrNq1axf1vJDItGbNGrIitjOe2PPEhwbvUvAwE/0QEd4Nfz4cGObjjjuOGc4tW4xbCZs5YUgtchtmq/a7ONbF7fay++oyl6HHIeVdYpvo/PPPr7NdpGRcslx06iXsvivVZO2Hk+XIpvbUjSqpTPeyld27d7N7sTkXIPu4cePGbPutosJ6C2TLx3kmT57M4v1VVVW0ceNG9nPbtm3pww8/rF0BxvshNmjQgIVry8rKZMfh8T///JPy8vJYMhHnCJUllaqdjsK96oWA7Qiuc+x52h3PYePpqjQu0iN4l2KgPGXMmDFRZSqxxiUDci5KqJBaUweq7xDPQO5Ix9Be2kFXj7qG7VeK1ZlqUSZXWVkZM7tbKLES5l0rYXnj+f3339PPP/9c+++OHTuyjXNk3CbzIQpj5MZhDFC6ahLqOSPZvn0726e1E6GAOfdP7Qg8nYEDB9Ivv/xi633P6tIaj3P3xpraQqPO9eDBg1mmf+S5htGUKruQG5cMQnSsLXWSfE09Rya1CLWl/bSbhW7FwsqcFAzbYs8RmWTwGpE5i1AtSktSYSWeKrgNDtOZHUza2PMU83jsRNAfrJN1awRy3iWy/JEHgft4xqlyXDGmcsfh57E/Ci8z8qYF9evXZ/fV1dJ71IgYAqt5ncA2y1SEXfv378+yyoRaJBTknnLKKQl9iMIYYVzDhg2T/uClLlIpj5RTg9vN6zzlwISM+uNUobzQuP0xOe8SHiUUd8Q8S7lxdqVDhw7sfu/evaLPY/5EtK9JkyaWNJ5OO2baXnXVVcwbFbJnE/kQYSwFFROpccLjCBVztMPvsXf9YrLA08E+lpjHw1H/XEt5l5hzduzYIVpTKTfOrhxzzDFMROHgwYO1Mn7hrF69mt1bNV/EdsYTNG/enN3jQ0vmQ+zdu3ed58PBlwBKRfhCHH300Zr8HRyOEuDpYI/ezvudZkHOu0QZyg033CBajiI3zq5kZ2fToEGD2M+ffvpp1PMQnAEQXLAitjSeixcvZvdCEk6iHyJqOcOfDwfyfyhTgcyfmFgDh6MX8HSwwNNS+s1ozOKxyXmXiHpha0ZMh1lunJ0ZO3Ysu3/qqaeYrKkA9H9ff/11tt120003kRWxpPGEULuYlie+YNOnT2e1njCYCN8m8yFC4g/hW0gBQv5PAEL0DzzwAPv5/vvv1+zvtCqpNkEYDTwd6DbbWiTBJJeUnHeJx1COIvWc1DgrAeEYNNM4/fANjTJA+GOCuAyAc3HPPfdQcXExnXjiiWzfd8iQIUxPHBnIb7/9Npt7rYglYwgwftCiRXgWJSnQqy0qKmIdViBWDE8Q3VPat28f9SFC3xYfIoqZ8cFDGB6TvdiHCNH5t956i6688kqmMoRWQ3gvZPhijxQGGY9x6mL1CcJqwNNBqyg7U1VunKqQmHcpRfjCPJ5xVgFbX9ATjyT8MWG7TGDatGlszn355ZfZfAuBGczHSOq0cksyS3qeKEV56KGH2F4mVIEQioU3CmN31113MSMKgxcJPkQYyZ49e7IPEV4nPkSUuEhlwl1xxRXs+QsvvJD++usvls0LUXgYZ65rK07JgRpxCY5+ixWUZtl50XKowBzXlJx3ifrJ1157TVKeT2qclYD3DGcjJHPDa8TGQVAGyZkoK0S9q5UNp2U9z86dOydcx4kPUezDjSVCjw+bowyXSbyEVCJytW83dq7Lq/Pv4v2HqFnrJoYci5R3iTAkIlTIkRBLDJIax7EmmhhPqF6gyzhCqahhHDp0KHsc+rK4idVMcuyDx1WzD8LRB3g6iKrYGZ+3rl6sp9qYa0zwLsVALeeIESMk5fmkxnGsiaphWwiwwxWHkDrafKH9V3iT6Q8++IDVUiLZh2Nf9u8qVP13Vldyb1YKeDwQSZASJLcDQX/dcKffa0z2Lc4xtnjEzjWM5jXXXCNqPOXGcVLceO7Zs4dlUKFbOko+0O8yMusS+5DYLP7888/VeluOCSk9WK7B7zTHnpcZwWR9yy232FqeL3BYmk8gGDBGOEPOu0S2PxIKpeT5pMZxrIlqnyTU+hGmRQNpeJtiJRwQUkfWlVi2Fsc+iNW5JUvQZ5JaBROCCRmJbXaemEMRtSpGlUPJeZfY5/z7778l5fmkxlkNdK+65pprqE2bNkx8plOnTnTnnXey+T+VUO2TRCgWCj2jR4+WfR1OtJjKD8c+7N6wR/XfWVxwSPXfaRfg6WBhahYhAS0oj4hmGFW6IuddArQyTGScVVi4cCGdfPLJbAsOpX2XXHIJM6AzZsygPn36SEqZ2hHVjCeEA1A6EgtcPHIC7Rzrk5ahfh6a1ScdLYGns2vXLltLv+3aUHdSzt9iTGcVOe8S+5lQJRPb15QbZxUwb1999dXs/rHHHmMeKLbgNm3aROPGjWOG06pqQYYaT6yqhKbTcmzZsoW5+xz7smeT+qvPwjze01MOJd89K1NZVnfBXW1gQ2wp7xLbFVjESG1bSI2zClBZO3DgAHOSJkyYUPs4mhI8/fTTLKr4ww8/1PYYtTuqGU/UQqJXHFZXcpqzEFPnqjz2JjNXfa1fh0P1X2kb4OlgD8rOmZzFe+suniqKjWlLJuddQgAB5ShiQghy46yCoGKFxFBnxN4tFgywAQBypqmAasYTbjs28S+77DImKBB5ASFWft1117GwBXQ4OfalvFj9FXZFcaXqv9MuYOJC6EyLRC2z4HHXret0VRjjecp5l5j/oHQmJRov55VaAaF/McoNpaKPgHuecXLaaafRSy+9RPv27WObyNhMhjuPmDhONrRk8Rz0Da3av42jjG1/71T9VO1cX1dhhnMELFQhkmB16Tc5ghH2KBDRFEIv5LxLOAZoSiHVrkxqnFVo0aIFu8/LE/8u7ty5U/Z5u6Fq3vQdd9xBv/32G6vzhOHECqyiooKpCkEbFmHbW2+9Vc235KQIAZ91Jx2twfcMouN27mZTeahumLayxJiwrZx3KRc+lxunJ6eccgpTfYu8KQHhWoCuKUURZSmooIBeOMCcnwqonvqFljSo88RFghOMtmHofiImlsyxJ2UH1P/y7N5sTHalFYCnM3PmTLIzkTbHHyGaoBeCdykG9gGRzyFVAyo1Tg2cXTpQWm576ef3biFKstrrggsuYOpx6B07ePBgVp7Sq1cv1ojjtttuq1002KGWVQma/ZXwPOHmt2rVihtOTtLs3cJrg6XApAWRBCsno8Qi4Kv7t3kj9kD1Qs67lFN6MktSF5I6EaWIvCmd05FxC0/1zz//ZFt1DRo0YJKsKFV8/PHHZfdE7YZqxrNLly509NFHx5TeGz9+PHXt2lWtt+WkCMGQMZ6GFcBkDd1UO6/4Pa66db6Fecao2ch5l6hFRg9hKXk+qXFWomPHjqyi4pNPPmGJn/A4keuyYcMG5igBpWFgq6Na2BaZZIJ+LQpow+uAwkEoV3gth6OUdKd1i8u1BhMyMtntTd24rdsgz1PwLsVAaBb7flLyfFLjrAb+vhEjRrBbOEuXLmX3qVKKqOoy6KyzzmK9NqFzixPrcrnU/PWcFMZVbVxRvNmBp4MtEjurMIUiop0eg0pV5LxLkJWVldA4q7N//37WxxTlKsOGDaNUQFXj2a1bN9bHE1lZCN/CmKaS1iFHQ3iyrawngL0sK0u/xUtlWU3Nod7IeZfYz0Q2q5Q8n9Q4KwGRG7e77sIFczzq+5FlO3XqVMrOzqZUQPUAfNOmTemnn35isfC//vqLXUzLli1T+204qQZXGJIlckKzO0Y1w5bzLiGAcPDgQUkhBKlxVmLKlClsb3PgwIFM5/a8886rdZoeffRR+ve//02pgia71yhLefXVV5kgQnFxMQ0aNIhmz56txVtxUgX7ljAmDTwdiJAYncmpJwGfMQlkct4lyvJQjoL7eMZZCSSmnXHGGUzRCmFalKmgEfsvv/zCtutSCU1jCBBN6NmzJw0fPpyp7a9Zs8bSChscA+GepyTwdNCMnqM9gncpBowmatxHjRolqv0qNc5qxhM3joZ1ngJw7+HS9+jRg1588UX66KOP+HnnxA/3PCXBpA2RBDGPh6Muct4l9jOR6yG2ryk3jmNNdCk6Ql3n8uXLaciQIZYPW3A4ZgMT8qJFi2w7MZvp7xK8S6nQLMpRpEK6UuM4KR62jXVR1K9fn77++mum78hDtxyOesDTef/99217Sn1e85R3CN6llAIPBAJwH884jjXRVe4CFxXKVwYMGKDn23J0xO328POtM/B0IJJg16iOy8DG1/F4l0iUhOqOmI633DiONbG2VhTHdLgquDCG3mBRikbEYh6PHfBUmmdBJuddQgABuR1iQghy46wGaoqvvPJKatu2LUuEQvvJ/v3709tvv2141xhLhG2hZYsLATWdUBXCv5WCcdu3b0/0rTkmprQwNdoRmQl4OrfffjvZlaryajILgncpFz4XSxiSG2clEHq+6qqr2NYbOqzAaCKLGK0of//9d2YP7LyFoIrxFPRphVUW16vlsOvBIM3RVAbfwe7du9PWrVslC/TtGM3ABK53q0Oc6+OPP57VN4qd6w4dOiQ0zgog5IzyQ5z3999/n4kkCGzcuJFtyX3wwQd08803syoLu+NMJkEIN3RSCf+30hvHnuzaqJ0cY0V5pWa/28rA00GDYqtLv0nh9YrvE/oi2pTpgZx3CePSqVMnSXk+qXFWAcIIaD12zDHH1DGcAPX81157bW1YNxXge54cVfG5tNufqjxkjJ6pFWjYsCHZFb+E8Yzs8akXUt4lPMqqqipJz1JqnFXIzMxU9LpmzZpRKsCNJ0dV/D7tFKS8Lh4SFj3nfj+deOKJts3kLNh+QPTxihL9IxFy3iUiatgTlKoBlRpnFZDXgpr9zZs3s/BsOAjbvvfee6wR9uWXX06pgGrGE/st7777Lu3cubPO43/88QedfvrprM6zV69erBM5x75Ua5jc4So3T8mCmYCnAw1pq+6lxaLikLiRLDfAeMp5l3JKT7G8UiuA/eV33nmHZddec8011LdvXxo5ciTTLj/hhBOoXbt29PPPP7PmIKmAasYTrWhGjx5d5+I4cOAAXXjhhUyeD709ETNHptbq1avVeluOyajW0MC5qrjxFAOTNfbT7JpLkF5P3OCEAvqXRch5l9jPXLBggaQ8n9Q4PYE4PUpmIm9KQUnU4sWLmReKefzjjz9movDQ8kVzgniqLqyOasYTacoIHWH1IfDWW2+xHm9jx45lxhNeJy6e559/Xq235aSQ56nl77Yy+E59+umnhk/MWlElkW3r9+v/98p5l8hCRTmKmIKa1vrDrnb1qbJ7Y8mbr6Gy/cpYfPjhh3Tqqaeyxt6QXK2srKQtW7bQDTfcwBwoeKEej3nqcrVEtdSvgoICOuecc+o8hlUYNpkff/xxqlevHlPjP+2009hJ59gTLUOrVeVcgEEMeDrQTbUr5QfLxB8vKtf9WATvUooWLVokNE4vkAnboEGDhLfm0K+zZcuWNH/+fLYVB1Am9frrr9O+ffvY43Ca/vOf/5DdcarZjDe85gqrD3xQMJbCSQYQVMBJ5tiTkFO7UFrQAE/DCsDTgUiCXTWjy4vFhTfKDupvPOW8S8x/Dz/8sGjtqdw4q4COWKhXRf/O+mFzugBUh8Cvv/5KqYBqxhPhWvTrFIDSBAwq3PhwEL7Nzc1V6205JqMwr1iz331wj3a/2+pg9W9XfB6JOk+3MYLxUt4lDAuSaMTk+eTGWYW9e2tquBs1aiT6vPD4oUOHKBVQzXjCSMKtx+pq3rx59OCDDzIZvssuu6zO66CwgXg5x574NEzFN6quz+zA07n//vt1V9vRi/15haKP79shXsKiJXLeJR6bNm2a5HNS46xC69at2f2ff/4p+rwgjoCSnFRANeP5v//9j6UwT58+ne1tbtiwgbnxvXv3rn3N+vXrmaYtMrY49qS6VDshg6oKnjAkBjwdZExKeTxWJzO7nujjDZrpH8GS8y7hLEDvVUo0Xs4rtQKCI4Sw7KuvvhpVkvjCCy+wn4cPH06pgGoJQ1DP+Oeff+jNN99kQsG4UJCBFc5ff/3FPgAhNs6xH8GgdnueoaD1O1JoAbyZ2bNnW9qrkeNA3kHRx0sKxBOJtETOu4QAAvYCvV5vVD2n3DirgIXBuHHjaMqUKUzjdsaMGax2Hzksy5YtY5nEt956K5133nmUCqgqtIh9T2TWSgHtQ0H/kGNP3NXapan7PFxhSAx4OkKXIztSUSIezdizOV/3Y5HzLmEwpVpyyY2zEpMnT6YzzzyTXnvtNVq1ahVTG0L2Lno0o1/pqFGjKFXg8nwcXSY6NSg9oL+nYQXg8Rx11FGWln6TI+AVz1A9uLdI92MRvEspeT6Uo0jJ80mNsxqQ3/v++++pqKiIhaFLSkpo4cKFhhjOdevWKX4tojNqwo0nR1XqN9ZuHypXw99tZeDxILPdytJvieD36L9/KHiXUvJ8kyZNkpTnkxrHSRwo2Clph4ntRLRKUxNuPDmq4q7WTiTBzeX5RMFkDZEEuyoMOdPEpymnAe295LxLCCGgKbSUPJ/UOE7iIL/mggsuYFKwUrzxxht02223qV4qxI0nxzJtwwp26l+aYAUwISN5w64TczAg/neVHijV/1hkvEsIIKAcRUqeT2ocJ3HQjGTHjh1MuKGsLHpbB3uzUDuCKhI0eNWEG0+OuheU06n7JJrqwNP58ccfLd1oORGCEsk5WiLnXWoxjiMPurqgPBKVHpdeeinbvhDAghJZwahPXbRoEfXo0YMo1Y1ndXU1C1PddNNNrKt5VlYWUy1CTenEiROZWLHcpjGEjbF5j9Y5Q4YMoaVLl8q+35IlS9jr8HqMw3iseDj6SuhVlOrfgsoKwNOBSIKVpd8SISDRJFvT95TxLlGGAu9SSp5PapyVgBFCxrAjxg3zsF7As8T7oTnJFVdcwc7xiy++SHfddRe1bduWdYGBnVAb1ZZBL730EuXk5Ki+KSsGGrEiLRr07NmThg4dSuXl5cwITpgwgSn/44TBVQ8H6kc4qdnZ2SxOjlUKVuw//PADffbZZ0zcIRK0EUIbNYRbzj77bGrevDnrWQeBZMgRouaJow9F+SX8VEsAgZJUwyGxF2oUyKQdOHAgCw/a1cOEF4e5TwwYLTTEBv379yc9eeSRR1hPW8zvaLsGTxQZ6PgsunXrpsl7OkJShUlxgiwyeGdfffUVaQ0assJQwhjCeIZ3drn44ouZGAPSpsO7nUNrF/3mmjVrxgp6BS1Q/IxuMDD8aOQdPgkhBRtC9jDMMKLDhg1jj2Nz+qyzzqJt27axDyeym0y8CP30oMBkdc53jtDsd+c0zqKvSuZo9vs51rqmMrLS6dvqD8ksYIGNhTgW5lpuX4QDY4F58Piz76CGzTtLvm7rqk+oMG8lm8sS7aoSi++++47ZAMiv5uXlGVLTCsM+Z84cdgyYm7XsL5qu5ooE4VO9TpDY6qdNmzYszo0iXvQOhdIHWqEBoYcoVijhItpnnHEG60gBz3nWrFks/BWe3oyLDapIguEErVq1oueee449hh52yRpPs/Ha+Hfp1It70/H9esU9NuQIUkijyG1mdqYiebNHLnma1i/dQuTA/x0Uopr1Ib7K7KcQ0bH9jqFnFzxGdgAeD65jLATt6vGI4XP7TeVdwmAicSXecXbhvcNe5zXXXKOZ4Rw9erTs80I5ELzNp556qs5zOCbM8abzPBGuxcoD3ptgsIzaDxW6tkA2CgYVnVyaNGnC2qTt2bOnTsNugI18hGShkoGYvgD+DR1HrGQilZFgmMO7CCSzcDCb54mVfn7uNtpY+bfm79W1TXfaXrCV9KYz9aQdoQ1kB+Dx4DrFNayXx2OWaMaPwU9N411iYYe5T0yeTyuv1CyeZ1VVFXMqcI95DLJ9WpDMuYPxVHPPWbUlEDbKcXFg1QEvDkbLCJC2DHDxIsEHQEIKhhN1PpGGE0A2C4S3VAOIm4c/Hw6+JMcddxzrMIBO6ieccALZgUOFNen/bSu70j/eFZq/H9Zuaq9SL6o3MuZr4JHaBZy/008/3fLSb3FjwJ8r513Co0SyophnKTfODsydO5cZzj59+mhmOIHa5SamMJ7oqoJsV5zEb775hhkciMWLeWRqu8+RKzGACzUzM5P9vHv3bnYvZjgBPFXsdcKDrKioYCszrNCEuiGpcXgcxhPxfbsYz89fmFf7GVWXu6hRs4ZkJWCMnQ5lq1N4A3bw1BASRBKcmMdja/SvVJH1LnHtrV69mnWNilzIyI3TEyTTiGUDJxv1eu9wyPa6664jLUE00HbGM1w3EFmsSOiRKgHRynh+++237Pfi4nzyySdrHxdKV5AUJAUMaGlpaa3xDC93kRonhIcxJp7wbCRo09a1a1cyA4s/X17786YV2+i0wdFet5nxepVLtvl8fsrMNG6LQS3g6cCA2mEhEInbpZ1iVSLIeZcICSKEivZckZ+F3Dg1qGzjpNBR0h1bvJu1c9MLCgpYBQKMcioJw6fbxZ3etGkT25fE6g/K/+F9RDnKKSk4Ug5yMN96HeHjEVII+gNENcEJS4NrHsLc5557ru1Ctx6XufpfynmXWLSj00i84/QEDavV3vP88MMP2cIB0T6hYbYRYAGJchVs0UmBaKjpjKeR7nR+fj774BB2HTt2LN1zzz11noewgZBMJAXi9UC4sIQxwriGDRvGHBMLqdCIlEdqBF73kckqf/NeshquCpfy11a6KTs3m6wOws9oBYhsTrt5n+WF5uqkI+dd4rlnnnmG/vvf/0aFRuXGWZ33dArZSoEyRGTW4tzKZeNj0aJmVxvL50yjFhMZbNh3vPHGG0VFC4TVxt694sYARhAhW2TkCoYQxhLZtNj3xDixTXDh93Xs2JFsQ1gyWtnBcrIargrlYb7K0mpq2qoJWR2EAqGCZUc8YYs5MyDnXQpC5YmMsyobN25kdfVwNsREZrRm/vz5rEUaFieYv1GXr1UdaySqL4EElQdk3aJdDOohwz2vr7/+WtYDjAfsIQwePJg2bNjAai5nzpwpGhKBNBOSh3Bhw0uNBOEUEJn0I4R+hefDwQoHveSQEHX00UeTHfljvvW+7MUFykPNh/brLyyuBZg4Hn30UctLv4mxY02e7PNVldo1IhAD5xiVBVLyfNOmTZOU55MaZ2XmzKkRLcH8K5dTohVPPPEEi7y88MILVFhYyBI4sYUodTOt8fz000+ZogNCp4iDw53GXqQADBdWCcjITRbEtSFesGLFCmak8X5iFy1AJuKgQYNqjzESSPMBCAuHA7Wi8OcjVzxIjDrvvPN0E4fQGl9Eso2Wja21YvPKbYpfu2XVdrILSsQjrFw6JUV5kf56x1LeJUKC2D6SCg1KjbMqoVCoVsXNqJAtHDIIhGCrTmr+1wrVjCfUTa6++moWQoLqDoxapP4CEhoQCk3WeGL1hqwuJElAQxG/L5YwAww6QGx869atdY779ddfZ6UqEJqPFH5A+BaSg+HHjBXOAw88wH4OVySyOptXWt+YuMqV73lWadg+TU8waWCvTe/JQw/SYiTX6N3jVc67xF4mNLfF9jTlxlmV3377jW2XQUNWcE70BuFiNZOADNnzfPrpp9lFA6F1MVEBgAsHzyHcmQwvv/wyffHFF+xnCLWj7YwY2P/E8wAeIlYnCCmfeOKJTOcWNVc4Xhj5t99+O0pcGyILb731Fl155ZU0fPhwJsMHbVx41NgjhUG2kzTfuiUbox7zuD2UmWWdlNTqyjgShjRs3K0n8HRwfeO6tJv0WyAoX8zp9wR0P9eXXHIJizyJyfOhq0e846yeKHT11VcblgSF6x6hWiNQ7VNETSfcZynDKYBU5uXLj9QSJgKyagUEIyoGMhAF4wmw8oPhhPGF0YS3ipOP/SLo4YqBLwOkz4RsLhhcJA/deeedkt0FrMrqn+sqLIH8bQXU5bhOZBWyc5SH0LNzrJ9pCzBxjR8/3nZZnKBkn3zYtkLn6IGcd4nQORbgWFhHCiHIjbMiHo+ndjsrUrpUT5599lkm/PDggw+yPWU9FyaqvROSgCB/F4/hSxQYRdwS4YYbbmC3eEBtFnR77c7a36K1Xpd+udJSxrP8kPI9sIpDysQtzA4mZGF/3m44nPJhW7/Oe71y3iUm7l27dknK80mNsyKZmZms0sFoEDFE0iiijOh8hUgglN/EFilIJoWjZDrjibh3LIknhEcRskU6Mcd8+N3RAgMb/9BftD0p4ilAt4kXAI8HSWtIYLOdPJ/JNB/kvEtBCjQ82qV0HCcx4ETBKMK2QNdc0DYXw7TGE1lmr776Kn300Uc0cqS4MDdafKGrCWLkHGuwb8d+shI7/9ml+LU74nitmYGng0nZLntp8djOUEhf6yrnXWJfEyV6a9eujTKQcuM4yXmeRqHaJwlVDaQtX3/99axoFiUpggAB/o29SdR8IrR73333qfW2HI0pK7KWUEJ6pvJLul49e3gAWHWjIxByDuwmz1cVI3vaU6U8QUwtpLxLGMzw0jyl4ziJY2TeiWpxK8SZ0U0FFwe0ZQUNR2wqn3zyySzhBmELiCS0bNlSrbflqIRUW9eKImuVcxTnK9+HKSooJjuA0q0xY8bYrgAf7N9VKPv8gbwi0hPBuxSr5cT5R1Ki2OcgN86KHDx4kMaNG8cEaFBHj8oEJIsicS1VUDWGgJUvemeiswmyWRGmgPoDDCtKQ2677bbaBtIca01SVqGiQrmxr4yjJtTMwOP5+2/tG5cbQWWZ/OfprvaSnsh5l1iAIu9DbCEayyu1EpAZvPDCC5maHHS5IVaDFo5QeoPSD5ynVED1ADx0BSGAjBvHOhTm2UP9xBFSHkxxSHjbVgOeDiYsrPrtVIQPvB554+jz+HQ/19OnT6e77ror6lxjPxMSofGOs5rHifwWl8vFxGOGDh1a53mI4+gNFivvv/8+Ox4I4KBFpNgCBpFQtH80nfGcOHEiq6GMPJmRzJs3j+2BPvbYY2q9NUcFdm8psMV5LD2oXK/2kAWF76UQ02y2A4EYIgi+GMZVbeS8S4Rkr7rqKvr444+jEoPkxlmJCRMmUFFREc2YMUN0rj/11FN1PR7U3aNMC2pzUudWyMY17Z4nUoa//PLLmK/DnifEfDnmorJY2pAg9G4Vgl7lX5Kg1zp/lxzwZODVWNmjkaKqQr6JhN76y4J3KVXLic4iYjWGcuOsArxNqArl5uayDlZmAFKwaMQN9SZ4ndDYhbGEiAM6vsAu4XgRlVF7HtO90A3hC7uobNgJZz3pL7XVV8t2Bx4PVt92SUYJx+eW/5t8Hn3/ZpxjiB2InWvMa5i8xeY3uXFWATJ4CIn26dOHJQlBOAYSpZBHRaLUvn37dD8mePlIVkKlR9euXWvPPfaYkcyECCcSWWFkIbWqJrpbMYQu0HeNYzJkNESt5HmmIpgwoJplx0VprD3N6kp12hsqRc67hBACSvHEOtzIjbMKSAgCqJb417/+RUOGDGEJQqjvR/lht27dWHcrPdm2bRsLFcO7BML5Dc94RvMQVH+88sorqr53UjGE0aNH1/n377//HvWYAFZcyMTF6sWIpqkceYoKpEs8/D4/V0UxMZgwRowYQXbEUyW/p1m0V99yI8G7FAMh2ZUrV0qGdI1q2xUOdGDFwvux1OHCpVWx9ZaWlsb2PXHdQZoVeuGQyEPdZc+ePVn+ix7gOMIrOAQjisQm6KiHK+Ah38Y0xnP27Nm1PyPOjFUAbnKg4XSqpDJbiQ1Ltkg+t3fzPurep6uux8NRDjwdhK6gNZpq0m+HCst0P9dt27ZlIUqxcw2JxETGJYu7VYioo3SEyF8/+a0XIQIFR2jSpEl1ullhTkd7MvRLxs/IftUDGMW9e/fW/hveL0ATj3AnDSIiaF9mGuMpdObGnhj6uSGFGer2YqCDCS6ejh07JvOWHI0oL5KehA7uKebG08TA00GyhJWTURKlskzfZthy3iWMCurZoa8qJs8nNU5PcAwoJ0yEcONzo0jCEB6D8Vy8eDHpxemnn87U65AgBLF6hJIRQkapJPSeYVzfeOMNljx06aWXqvreSX2SAwYMqP0Z7jpiy+GPcazD/t3SSi156/fSmUP1TUHnxMf+/fupVatWKXfa/G79VZWkvEsYTGh3xzvOKgiOT05OjmgHrU6darovFRbqJ7iCJCwkLv3www/MOMLzhOHEXqzQaQjOHcK5kIdVk3Q7CPRykic9I438EvJuPq81MgQTyQpGKMrKSRyCx4O9JyR0pFrYNuDTN5lNzrvEtQR1tZtuuinqmpIbZxWQZSuUrHgOe3rhCC3K1A6PygEDWVBQt0YdmbXY20XpJPZpjz76aLr77rupe/fuqr63arPGgQMHWNNo3IcDRQd0WTnuuOOYS41YNMd8yGXUelz6FqInSiJZwXYow8FkjLCtVSflZMhqUHcC1xrBuxQ717j+Fi1aJHodyo2zCh06dKDevXuz78xikdCs8JhgZI0ENgcdvr7//ntWA6224VTVeD7zzDM0cOBAKis7sncGvcOzzjqLxcGxKl6wYAGde+657IvOMRdBn7QR8eisH5ooZQk0t05kjNkQpN/sKAwfC3eFR9f3g2GE2IGYgcR+JhJlxPY15cZZiQceeIDdjxs3ro7HB21leHzg9ttv1+14kChn1FahasYTK65evXoxFzk8Gxee6KhRo1iZyvPPP89cfuEkc0yEjANWUliTom52KovjTx6pLLK+8YQngMnLDl503OgsqiTnXSI0i3IUMSEEuXFWAr2Ykd+ydu1aNt8jbIpkUSTuIGx7yy236Fo2hXONxiOWNp7Q1uzSpUudx6DsgFUY1CfgNmMjF26/ntlYnOQpP2gNA+P3xu95+XXeM9MCfMew12Z0Jqch6PzxyXmXKNcTWjHGM85qIL/ljTfeYIo+WBBADB7tyOAs4XE9QVcXo3SdVTOekG1CFpYAQkjLli2jvn371mkA26NHjzp1ORzzU7DLGqLx5cUVtm/2LZcwZGXpt4TR2dmW8y5RsI+QpZgIgdw4q4HFwS233MIEb6qqqqiyspKWLl1qSGNqdKmBOA9ueqPaMgg1nOH96vDH4KSec845dV6Hiwc1nxzrcHCP8gbTRhIIJJBtm8AYs4HMToTPrJ41HInLpe9+phLkvEsIISDCJpa8JTeOkzjIqbn55ptZf1Hco1wFiU2o8RQDz5nOeKIRNnQNEaJFUtAjjzzCLpTIwlQUq6JwlWMdQhbJQ9m4fHPcYzat2EJ9zzuB7KBtazf2btFfaDwWgncpBkKywlZVPOM4iYPaUqHlGCQCcZMCr1PT81fNeP7vf/+juXPn0v3338/+jT8G2bdnnnlm7Wt27drFsm5RB8UxDzGzNC3inIUS2L/Uu05QC+DxtG/f3vKlEJF4FWZ561mrK+ddgoYNGyY0jpMYZ599tmHefLqaG7cI1b744ousWSr2OtFDLRzU3CBhiAvDmwtvjM4VVuHA7oPxj8nTTw1FK+DpLFmyxBbJKOG4KlyKXudxeyg7J5v0QM67hFcDQXSoPYnJ80mN4yQOEpaMQtVPEhlX77zzjuTzt912G7txzEVZib76oFrhdflU7xdpFaxeAiFGSKZNXjh6/+lS3iUMZnFxcdzjrARyWOSqJb777jumcZ4K8GUQh4r2SOvaWomMrPgv57R6OhcKagA8Hkxq2BaxU0jQXa1MC9bvxaJJH89TzrvEAgZ5H6hrF5PnkxpnRa644gpRGb5UymdRzXju3r07rtermfXESY6g3x5eS0lBadxjSg/o29JKCzAZG1XrpiWZ2Qql93Tck5fzLmE8oaZ21VVXRRnPWF6p1ZgyZUqtEHyqkq521pMS1M564iSH24QlAYmQlsB+kh3KOzBpo3AdLaHs8PcIuKqUeZ4BHRd/ct4l9jMhRh7vOI41caqZ9SR2Qx0OWtkI6cSQcULrMo55cCtMzDA77ur4/w4km1gdTMw//vij7fY90UdWCVWV+l2/gncpdq6RtY5yFLHsdblxnBT3PGNlPW3ZsoUVscKAYlOZYybsUbidt166l6IUuxIYYzbg8aCDhN1wOJVdl75q/RZAct4lkOreEWuc1Zg1axYLQ8OLhp45KihSbStOt/gBTjDqQFHnOWHCBL3elqMAuwiKu73xd3/xeKzRMUYObIFAJMFuWyHbVu1Q9LriffopYMl5lxBCQJ27mDyf3Dgr8tRTT9Grr75KM2bMoHvuuYc1oX7yyScpldA1+A6N29NOO82Wq2QrU5Rvj0SGoDeoais2q4AtETT/tZv0W/72/Ypet3OdvtEDKe8SQgiod8d9POP0BNcJjjHyphRsxc2ZM4f1aa6urmbdsiZNmsQ868cee4zV+esJekgjqhkLiFPgtWriNMLLiWyYzTGWeln20Br2Vvs0U7ExM/B0xowZI+rxWJk2XVopel1u41zSCznvEo+hs4jUc1Lj1CCjVTVltq+UvKXlqiOEMnHiRLr22mtZB63s7GwWUXzooYdqQ9KPP/44azupFyjRevbZZ2O+7rnnnmOKd5at8/zrr79YgS0SiDjmwefxm0oCjRMf8HTQrQiNGcxYQ4hQZSKJMmUl5RQMxR5XXlIh6e0pAde1UqOG90G9Jvqniom/w6hIicZLjdOTlStXUoMGDVT/vRdccAGdfPLJrNPK8uXLoxqCGL3tpMXWlGrGEysSKdBdBa41EoWwL8NVhszFut82xHzNoaIyatayiS7Hw4kPIRnFjNJvE4Y9R98s+5JW7V+myu9rTM2plOqKeix8YC5d/UDiv/Phhx6mpyY9pei1ct4l5jaIBKA9Y6SBlBtnF7p3786MZ0GB+VoY7tu3T1TUIRlU+7bBXRfKUaRAv08IyI8dO1att+WoQHaj2GEvV5mLiBtP09K0aVMyI9v+2kGtKzuR1yst6RYPmF/U2NsdnDOSQocDLlddcrXicXLeJQym2+2Oe5xdOHToELvPzdU2jP7uu+/W+fe2bduiHgtf0GBf9qeffmJlkqY0nijSlgL9O9u0acM2q7U+sZz4Kd4bO2HIVamsYN3o0galeqhHxpDlwQSBpI+DBw+aLmzLPg6Hw3TH5Qg4cViM/dsL6fgzeikaJ+ddIjT99ddf09ChQ0Xl+aTG2YGDBw/Sb7/9VqtxriXILBcWIbhHUwTc5BZc6O+JhCZTGk8juohz1KFkf+xU/9KD5aY/3fEazpoxZHkwGZeWxi9NqAeOIJHD5AnNFSVVil8r513CeKJ045JLLhGV55MaZxWWLl1KhYWFrEdzWlj4GZrKSCKqqqpiC4d27dppehwwgkKUE9uF2Eu+7LLLJB23tm3bsj1ZOHBqYr5NEo7u7NsZO/u5JF+/WjpOfGDS/uSTT+jKK680XVKXx+Uhn8/c9acehQL0sbxL7DlD6SnecVYBeSuQgGzdujXzLhs3bkx5eXm0atUqtjBA9GPmzJmaHwe2CAWwj3zeeecZoh1gzU+RoyqBQGz3y+WyfkmHXcHE/P7775tS+g29YgMJ1N/qSUmB8uYAgncpJc+HchQpeT6pcVYBNfr/+c9/mCe3cuVKtmBbt24d8/ymTp3KHmvZsqWuxwSvF2UoRpCw54nN70SBy40iW445UCIU4KuyvgasXYHHM2/ePDIjIQULM6OpjkPbWc67BPDGEhlnBXr27EmvvPKK0YdhGtKTsfgcexD0xZ7gdm60vgasXYGnc/fdd9NLL71kulIIvwWMp4/1A1V+rh944AHm7USea/z70UcfjXscJ3EGDRoUl9P2888/k+HG08rhB04ESFyL4Xzm1M/ip83EmLW3YsBvfi3XQ4XxJVtJeZfIqD3jjDNo2bJlojW3UuM42jUkAUJykdplQjxhiKOomfCm5TzMblbgyYwfP57MiNPpoGACWdBaEpn1uvb3TYrHynmXSASaPHmyaEKQ3DhO4uzcuVPSuduzZw/98MMPTG/3jjvuYDc14caTowhnhrlzy9yuxMsAXNUuys7JJqsC6TdkP65evdp0NYRmM5ygMK+uQlHArdw7lvMu4dmgEF/Mw4nllXISQ07qtXPnzkzIHqHdCy+8kH02akrDqjYjLly4kIYNG1ZbKCsGVO3xGjXU7ZEe/cwzz7Dfh7oiXLBK3HKkNp966qlMqgmqLEOGDGH1S3KgABevw+sxDuOlFC3sitlFEjxJCLx7XeqIZhsFvBqUCJhyLy1o/mvFka48nCfnXcJAQixdrDWc3DgrU1xczDJsMfeiLZkZgfGE7i7shZqotgR6/fXXWTbZO++8I/kapDTDjYbaA1YEyYDecV999VVcY+69917mwuMCR9Eswjc4ZhzTZ599xhq6RvL555/TVVddxcIAOGa0VcOmM0Qh1qxZQ1OmTKFUoKpMeSG5EVSXJ97Joaq0iho1a0hWBRMXhOHtLP2mJo4IA5aVk6l8rIx3CY8ShlPMQMqNszL3338/FRXV9eTNCBwsaKuriWrLoBUrVlCfPn1kFfsbNmzIwktQ3U8WhECwh4DCYwgRZ2bKfwGgbQjD2axZM/rnn3+YkPaCBQuYF4wVO4p/I1VaSkpKaPTo0SxTDsYVm9O4R/cKrLJQ26Rkw9oOlBSYU8FGoPJQZcJjq8qqycpgwm7RooXtmmFrhau87kLQVa48qiLnXSIpBRE4MX1vuXFW5eeff2bO0i233EJmBi3SUIMKp82Unuf+/fvpzDPPjPm69u3bM+X9ZHnwwQfjev3zzz/P7h955JE6TWlhhNHhHWn+s2bNYispgTfffJPKy8uZ9BPCwwKtWrViKed4DAZUz/Y7RhFMNy5rctb7c5hSzS0yEpD1cjMVta+SGhvZ0qrIXUQPr3usJpeK6bPWPI67Z3s/TS1yWpBZwD5nMi25NAWRZJMl3O7etS/hsXLeJaJTUL9B30gx9SGpcVbE5XKx7li9evWicePG0RtvvGHIcezevTtmNy/M0UgeGjVqlDmNJwTflTS5hjai2isAJR80VoRg+PDhUc/jMRhPFJqHG89vvvlGcszFF1/M/g54tAj/6v03qQVWyc0/zqaQx1Ejko6wElbOSO+mENMmRfeJ9NbGfel/a/Mr/TV7FY1JYoV7Vr+z6Pclv4s+16nH3Lh+V+k9JTR32udkFjBp41rFNWm2ybnlFzkUMFnEv9FF9anBgxkUrAgRpTuoXl/le8WCd3nuuedGhWBhIKUEyuXGWZEnnniCduzYwfozG5mkhhKtWOcT5/6YY45he86mNJ4I2f7+++9sJdChQwfR1+A5JBTB29MTtKTxeDwstCUmWix0AcAeZjgI74Y/Hyk4fNxxxzEvGqubE044gaxIIBigeh1jX/xOA5UcQ/WDdMJ/etOfsxIP98db5/Xippfpn4q614PA6U1OJTMB4zlt2jQaPHiw6Yxnest006X0l3pLqMHliWVXy3mX2N7Bc7hFJm/JjbMamCenTp3Ktrr69+9vqGAO8lCkvtdCN68BAwYwr9O0YVvsDSIGjo4Cb731FstuCgdG5qabbmLhJbzWCNdeSu0fXjMKmNGPDi2DsG+LcG1ZWZnsODyOvwviyFY1nkGF6ZAZZGwJhDPNqesKN5QeIme6+CQXSDdX+QU8HjWVU+yOP5R4HFnOuwRS4fNY46wCFgE333wzmy+N0pQNx8icE9WMJyz7F198wRJqICDcu3dv6tq1K3sOOrbw4rD6v/zyy+m6664jPUHsW2jGLQUMKBKGBOMpjJEbJ/QmxRgloOuAGDg/wrnSn1DdTT0JMp2ppW1b6ZNOQPL4Es/s1QJ4PP/9739ZKr4py1VMRnUg8QQxOe8S/5Yqh5AbpyfoqSz2/uvXr1c0fvr06Sz55u2332bJl6mMqhGVjz76iJ5++mmWnPP333+zmwBWKvfddx899NBDar4lJ0nSgpU0uc0fFAw5Du9zYnXpJKcjdNicYueTKMdhXJbgkOzt5ArqO+Gcm72VPvZgSRFkW8AA4SGcp7MzN5OZCJa/STmhtygYnGQ643l1o81UGqxHZuL83G20iLyHl41O6pQWX69aKe8SCUFoj4U8CDEhBC2Tuno2P0CtjpJeAC/JqSLlvWOkI3hIuEQY9IYbbiCjwN7x3r17WXQTCUtybNiwgUUHkaiKkLlpjSdi+Ti5yITFASPDCeDA+/bty2LQRgBhA1BdLb3iRCNXIJTaCGOEcSiziTUmFlKrOymPVA8cDj81zzB3+vyIFolnRyZKv/obqF/9CpY/FU6NIdW37VIs0jxf04RxzYg8C4gyhpKZOL9RIZmNZrSbXuuQ2P65nHeJ+Q8yiVLyfGoX6ScCvEal81UkY8aMIa/XS6+99hoZBWwKEuNgUyCUEwu8DtFOGNutW7eydmpqocnONfamkBSE5ry44WejDCcQEphwAsWAEUTItkmTJrUXFoxlo0aNZMcJj6sp+aQ7fvNNbuGE18zp2YzA4QhGGc6ax3EzT2/ToB/X4Naaf1SOE60x5ETgPLIwjhd4lyhNk1IRksp4lhtnFebPn8+2sFDad84559TeRo4cyZ7Pz8+vfQyli1qA8kEYcOy3KlkE4DXIskXFBUoR1cRsiXCagDRliCgcPHiQfcBHHXVUneehCQoik36wbwsRBTwfGR5ACAaNYJHBdfTRR5NlcZh7sg0FDh35h7+AqF7dz0475OorEhdkUJ2KV+v8M+TfRY6MzmQ2gn4vOdPNEr5VrigUj3eJOQHzAUrXIpPb5MZZCTgZKE8RA3+38Fyk+L5aQBEOVRNianBSDB06lNXmQ2FITXF+a3+SCoGyh9D37dNPP416HklO4NJLL63zOFaR4c9HrsJwgWCPw6o1noxgfPs9+hO+UjfLqt0sx4HspV/q/tstPrEZj4lEHHyJy8nJeZfY54RxEdvvlBtnFRDVELvtPNzZBEmPwmNatciDuhuSnuIF+6MoWVQT636ScTJ27Fh2/9RTT7HYtwC6HECXFwlNKKUJBynZCN9CQ3fu3Ll1hB7Q2BaEiypYkqCyTGHDCM+MNI2hN1ORe5hnDnwmbVIfNJHxDJWIesZKgHeJ/Uux5B8YDdRAioXO5cZxlIMtNmE7LR4wJryCIqWNJxRVILQs3BAHB+GPCQpBAB7iPffcw7oAQKAebj86paDIFvsQSL2ObFaLLiqoWcVqESpD8F5HjBjBwsDbtm1jBtny0nxBk+u6hq/UHXqG/dIsstsRUbPoFVdRMhynmdqriOWdKvseyHmXKEdBUg3u4xnHUQ7yUpQo2UWCMRirJpb9JLF/KSYwH/4YXhMOVFhgOF9++WUWO0cSE4wq4uBSurxXXHEF2/eEx/rHH38wI439zzvvvJN1VrE+Jvc8g2E1lSE96yvlCunN7D3EP7HoQgh7yE3JtNd8SNkeneBdIgkyUtkG+5zh5XlKx3GUg7kX8zASgLAdpwRUSyDCiFaSamJZzxN1RlIxeOEmVouEx1BGA/cfikLYRI4laN+vXz/2Orwe45DubQ/DCW/OXHWBUbg+OPJz1Wwjj8QieMy/CDIcMeOlzHjKeZd4DOUoUs9JjbM6nTp1YvMtonFaAwU7zMFwZpSC18LYRua0pKzx5KhE8JDylxrxxXf/deRnX9jPHAoGTZS4FPN6CZp7gaFwT1bwLqWkIpHNn8g4jjJQJoPMWSxSYBTlytfwHPo+47UYgy4wpg/bYg8Re4sQY5dCSjyeozfxfJkxWevtqdYIbdRg/qa7uuLXppZOPfwmNZ4i+5shZclo8BxRN4iyEzF5PsjXxTuOoxzUmX7++edsu23ChAk0c+ZMloeC5h0oYRG261BeiMoK1OKjGgJj5ORZDTeekKUS9gblssoQ87dysbCtCMWx5xnyJlUjl3y9pZn3Gg0gZO4G5RQK+7wcZgonh5LKOpfyLjGnofcvsvPFEoOkxnHiA9tsS5cuZRrpUG174YUXol4jZDxDve29995jNftqo5rxRN0jZJCwwkJWU+fOnROWgeLoiP9AnIaMf6bmwWSNMqMIWyAHzb7Xp+z45LxLZOUjp0JKnk9qHCd+kPi5du1aWrBgAauqQEgc0U4AwXo8j7raiy66iLQiXc3mqIgxYxWATFQemrAIznhCSGafAFMMUyXhxMjLYZ3WTQzLBo6NnHcJo4kQYrzjrASafqBv89q1a1m9O4RiWrduzcTiEZI+/vjjdT0eGEctDaQcql3RcJ+Rho1aSm44LYQ/jgy5gLl1cFOOgMn3gENh+5whM4VtxaY9ZeUjct4ltqoQbRPbspIbZyXQNQuVB02bNqVzzz2XeXfYU5wzZw5r/oEIZKqg2hIIXUh4EpAFCcSh2mMahR+O/qIRCRBuME1/7SjbT5fzLuFRQr1MSp5PapyVgOcMI5kVIUn6yiuvsFIcqLIhScfK3rVSVFsGIfsJ9ZMcqxGHwpDpJ8AUw7NG9OGgV3n5kaaEZ7AGa/ajTItD2WQv510CqW4iscZZBdS8i2l533HHHUzbFko+6KGZCqhmPJ999lkqLy9nvTx5Jq1dLwGT77GpRNBvkQkuIJHsFTCJvm0oPGHITN17RMpm/Mp6xsp5l5j34F2KzX9y4+xCxuEaViPbT+qJap8ktGEHDx5MU6ZMYTU10Hxt166daIwfpSpqtobh6IXx9WnBgJ+caRpPQArKKnQ5jkT3PM3i5dXRTTZRGzfJMixlwLtE0b2Y8QhvOqF0nB2YM2cO61rSvXt3dksFVPv2P/7448woor5mx44d7CYFN55WRV9NzqBfxNMNlhGlNdP4nd0KX5N4U2V1kDDyZqmprCPqn0umJqCsBlPwLhGajFQLQpke9v4QwoxMmpQbpydo5yWW0ImEz3iYPHkyGwOpvI0bN7Kf27ZtSx9++GHKJIyq6nlyrEgcXoq/pm+ffrgkvBmNjacSqbY4PBXt2C7+sF/icb1B8/Lan6U9MlPgVCb+IeddwnFAvaFYS7JYXmmy9G+0hXo2kw49780qlrpaEuL777+nn3/+ufbfHTt2pHfffZclE6UKqhlP2wilpxxx1G4GwiZDPahaFf2YewVRZntt3zdUbZH9X1fyGdRaEghr/RWI7qFpxYQhOe8S+5mzZs2Ke5yeoKmFGuI1P/30E7tHmzXUfE6cOJHVekJh7uGHH6ZUwNpFRxx9Q7GO6Cw7TXGKhP5D2nduUFSTaKbmzpGE9pIpSAvTEnUYHeJWRydYzruUSxiSG2dlGjduTP3796dvv/2WeZ3IZYGBTgW48Ux54vgy+zfre7b8JeJ7nloTqrS2NJ4/aL5ry/Q1qco+T8G7lKrlhGiAWJKk3Dg7kJGRQVdddRVbHMybN49SgYQ/yUGDBrHEn3feeYdl1eLfSsG48Hg5xyIEdZ4AxSZchw7C9EqaFZvZgQhKJ+vpS7hIghnC3EcEw6Pwr1M0Hl7lqFGjWGKMmDyfWA/hWOPsQvPmzWu7mqQCCX+KixYtYkYQXbqFfyuFd1K3KCGd5eDEjJge8mbhe3WWFIzwmk+ezySrjZDXlZRYiJx3CQGE9u3b0549e6IyauXG2YXFixeze4glpAIJG8+dO2syL4866qg6/+bYGZ0/40Be9GP+PB3eV4FhDJi5btEkdZ51Fj/mMJ7krE4qcU7Ou4RHuWTJEsmQrtQ4q4C/raKigi644II6iwAsGl577TVW65mdnc3Ct6lAwsYTqcly/+ZwkiYgYgR8Oni/DiV9Ms2cPWoSQxUeYjfNnqdUDa8yj1DOuwToLJXIOCuAUpsbb7yRhWf79u3LWn8VFRWxbNuCggIm2zd79mz2d6YC9gy+c2whkkABgy7PoAIDbRYVH1FM0mg+WGm+MHdQyngqS7KS8y6xrwlltV27dkUZSLlxVgGlKA899BALz65Zs4YZTkjxderUiYYPH0533303devWjVIF636SnKQJ+k2yNyaJWPmAHskITnWSijQk6I9D0N8o6iTnBEy+sFA+FUp5lzCY+fn5cY+zCp07d6ZJkyYZfRjWN55dunRJ+E2RMLR9u0lUUFKEL3aspiYZ2dSv9ZGVYdBfQeQLSWYliiV2OSO6QngDfsrKqKdJ4XcV1adIUbdqaqS9KF7wEAt8SplH9lzQ6M4lbmabxGy4SYK2TBhBOBaHHiVGighJfLbKQqly3iWMI5TWENqMTAySG8dJMeOJi0DQso0Xnm2rL0sObKNb/jeWij+pyYZTwtlnZNGvy8RCXNF7V2ffOpIWv/4hqc2IVcfR0bkNqNxf856N0z20oaot/XAxacp7+9rQ0IY1Mvjs6j58iTND5SAKhIi+K2hJoxqRYeyodFGuJ5MapXvY4dUcooOcFKK/yptTv9ZkOD+VdqTT0pzkdARpbkln+rcJNNG3V2VQE28GZacFWKMXfKY4Zz8Xt6RLFZwzOe8SxvPHH39kamuRxjOWV8pJwbAtNo6vvfZauuyyy1imFcd8bCrdT81HDKDmw8+mb8+9k9o3aBJzjJTnGcnxXz5J7ZpqkyCQ53VQnrcT6c3Cogb0f5uvlvU9B7XOpVFkHKU+H13613DJ5zf2IsP5rSSD7skbyX5un9WY/n2C0UdEtLvaS0NXoSl1dCbwpQqOT867xH7mRx99FPc4jjVJ+FPERXLJJZewjeOxY8dS79696b///S+tW7eOWrZsyVrvyN04+pHpTCNHmpMc6WmUkZnOVsGxbkgEUPI6/M6QU5v9P6MUQB0B7E3J/U0OSjNYnCtQp4bSnJR7EbmocdfL3ebYo61g+/yRn+3hkIICBO9SbP8SoVmUo4jJ88mN41iThGeAK6+8kr766iuWojx9+nQ69thjmdoQaoCQqjx+/Him5cgxAWErXS2+vJVebVpgGZVisqY0tjbsP4eM1Y91i0zQpoMtQg7/aJKdWE+SmsSCdymWNYtIDVp+iUVs5MZZBQjifPnll3TTTTfRMcccw0pTcnNzmeMEYfjKSjPXPqtP0svnpk2bsk4BS5cuZUlA6OsJ1f6pU6eykO5xxx1Hzz77LKtv4hhDeZiqiiuo/qTrCvhs5Xk2y4rddaJxhrFbFAfcFbLPm8HD2Vl1pJyn0iTlM9WB5I5DzrtE0tyYMWNEk+fkxlmFDz74gC6//HJ666232N84dOhQJgoPgZwJEyawhUNhYSGlCk61U5mhqo/mqMuXL6e77rqLiouLWW0QTizHGMp9RxJ/3F71DV26RiFMo3wVtz/2OXIbXObjjuFBeZI0EmrQMutIXrRpJBICyX1uct4lhBAgTYf7eMZZBWzT3Hrrrayh94YNG+iTTz6hBQsW0ObNm6lPnz60adMmuvfeeylV0GzjBopDKGdBd3Ekn5hhJZyqVPmOTBjekPqTanqGNqEosWlGj6mnSoEnXenXJlStlHSH/FfXDEHSKt+Rc2T0HrFAVZJbDHLeJUKyCGuKhWblxlkFZBG//vrr1LNnzzqPt2nThmbMmMF+njt3Lnm9Zq8fVwen2jHx9957jwYPHsw6rdx3333Mpb/lllvYSeUYQ5vsI2HIrDT1a8yCAW12JwMGGQWnkvIrg9eC+VXydabVYdEGoyj2HGnz5Tb6hB1mT3Vy9bly3qWwjZXIOKvTu3dvdu/xeFi0MRVI2mWAR/n9998zo/n1118zA4pMTcTDUcIyZMgQ9m+Otqw7+AYd0/hG0eccAQeF/AEKobAtEBT9AiM6kO8qo86NWsT5zkFyOOy10gwqMNEhh7HGoNAlL3dX7HVR85yGZCQe5sHXyiSQGchySn1uypZlct4l9jOROImWXGLyfFLj7MCOHTVt8PB3Sy0g7EbCnyT2NN9//336+OOP2cWCWP7ZZ5/NDCZ0Dhs1MrCCPMUoql5LTz75GH0y/TbZ1x1/+snUe+RTks+ff1sX+uG1+JSfruuwlNpmYqK4n9TmzKYb6Z+yww0HmCyMg45tqH2h+dUdV9CS4pZH1HtYMX2NLk3IUdNp6+wWxq6u+zf/nvIqYeYd7DiZs8x+IMpOc1Nu2hYiMrYk7IYum2hxUc1C7egG5ihVGdhqIRW63UwYAWvJdGeQgkEntchWrr0rZRxgOEpLpZsKmMGoYN9VLHS8fv36pH7viy++yO4vuugiyszUoeeulY3nGWecwQzm8ccfT+PGjaOrr766tj0ZR1/8ARcNv7M1vfDEPGpR/8SEhA+W7XqUDoQWxP3ePRoXkZNqmuCqzcXtNrOb3nRuVEidG8lnhzspdkaulqRTMY3stEvyeV9A596rIrTJKaWh7Wra2DnJHAIqOc4MuvSoPyWjaLEEDOS8S4xHEg3K+MTk+aTGqUH/LDedlSvt3c9L1y6B7Ntvv6VZs2axv+vJJ5+kVCHpGAKyrh5++GF2UwomcMTGOSoRQkKCg9LTa4QLEiEjux6leRILrZklGUQ9Yu/hBg1uOO0n+Zq6YMj475cvTM/W6PMlUOHdLfmcP+Cmes4c2fFy3iWMJ6JxiLyJyfPJeaV6sXLlSlZKqBabNm1i0UYszCdPnly795kKJDXr4YRhRRXvza6b5kYRctQkhwST2Ht0IIaVIAGTTIzqoWSVbuzfHCD5MKjLBC3T/HXOkTm6qvhJelERkHku3EBC7ECsegD7mfPmzRPd15QbZ1Xy8/NZmPbQoUNMZe6ee+6hVCJh44mLIJkbRz2Chz/GBDT6j5CWeEIH5MhTD2OLQbDXKYsJvmKhGN6xEQRlFkahUOwFkeBdis1hgUCAlaPgPp5xVqSkpISpyeXl5TG93ilTplCqYc/UrxTDF6hpTxWQbPQbm2p3dULGl80FQW32s1Dn74zIbdBD0hVzX6xyPMOFYiCSILNmcftN0HwaNiStZlHH8q10xucL0OTXFpDH56fxt55PmZkZVM/XgqrokGgeQCiQVhsVw1i8PjK5RvAupUBjaDFijbMSkOEbPHgw27IbNmwYzZw509LiD4nCjafJwSo21mq1rKArbdnYjI49oQv5cmKHxLEfEzkpLPmmH/mbb6Er4mwE/8H0C6hBZnMaNplU56WHR5IzLURpaVBGdVAw4KRgwEHDPyNNmfXMJdSibRkFfGkUCKRRKOigYAhZrUFKzwiy4ynMb0JXvkOG8f70QdS5x17yedPJh+P0OygtI0T16vkoLT1Ia/3d6eS7yVA+mDGUWrTbi/0EKj7YhIbp7Jx8+u0qevWVF2nvph9pokQy+IABA2jxYqFVX90yrUceeSQqAQbfx7vvvpteeumlqO8Q/g1NbzHkxlkJ5Kqgg9aKFSvowgsvpA8//NDSf08ycONpYtbv2Ufnj7qJCpbEzoJt0bUHPbRdWWswTBg///xznYvek5VGP31+Nj2BTlxxsHd/C8pqrk0db1VOGmX44B2k15SNpBH5s0hzDlU1ovKN8qVWh1uMGkZeYTMqyG8m+XyHPsaXhuxzNaT85TVqNH4Dkm13+yup+UnnUYu+59Lp/brSizdfpigDvaS0ki54eCaNf+COuLxLeK0nnXQSrV69WjRxT2qcVcACYNSoUbRw4UKmaQvhm1Su4efG08Tk7S+lpqecR3m/fC37ukOHSmngpNmUP+VeRV8AoSYrHHc9oqpO8a8gXe3TKDNTmz1PX4s0MiK1rLqd+VfSnqPSZNNb0uob/9WuaOAkIyt6qn0+8reuMWI+R03GqxKys+tRsEEGOSP3DGJ4l3gOIUwxT0xunFV4+eWX6YsvvmA/N2/enDUEEQP7n3je7hj/DeNI4g55yeF0xvzS18vKJGd6TZ/OWOA1Dz74YNTjziT2pBzpTk3C06GwxAtIADL/wOEQzdYWC0XbmVCNPILk82v3HaBUwR8I0h95u+i09nUjL11aNGLXELzLTk0bRl03Ly9aSkN796KuLeuGaytdNYlD1R4f1c/NVuxdwoPt0aOHpGi8nFdqBZBVK/DFYSMqBjprcePJMZYqH3nqx06fRzJG1n75FlWxvsTZWyopc9cRLVKluHP8VC/JNdjDt88mX4PNNHnKM5KvadKyFR0qrDEI9aZHG/+BAweyZsNqGFB3jpcoHZkuh1XlhHsBB5HDa2zGkK9egNK9hxctkceIBB2N2sTFQ7PlpeRtjM/DQWlV2pWqTJz8Gc1+8xXas13Yu6xLj94n0SMvPUCPiDz3+cBh9PvCz+s8FvAEyFnpolAgFJd3iTK8Fi1aMInSSAMpN84qwCjixqmBe54mJu2Al7p8UUn0hPzrXK5qarbWK6sgFOtLnOEOUvOt8Wfrdp5XRg0bJCdCvmrpNrr7sctoknei5GswMQ09dQJdcW0/umnsEE09zy5f1SxEpM5k7ZQqrXSoOZ2+KYt9fAaTVR5gN60J7K6gntln0nbvj6LPS30vLj7pMerTLTpDLt1B1GZZNdLYo56T8y5hMKVq2OXGcaxJKhboWYaMdGXhnbS0mk378vLY3qfUl9iVYJ9P6Gw7IPiaJI7DYWepW3Z2NqU56pHfF4x6DkYTEmFq1dDhr3Ek8bzW+Hx+Ux+f3vjdMI412xtiNyS1iD3udKRRTm50wkvQV7P8CIhcToJ3KdbUGtcfylHErkO5cRxrwj1PE5NeT9naxhGsmSrT0mJ/nFKhJcFpU6LvGYlHhQbbAYU9ALG/FQmOedq0aaz2TOzYlZT71Pl9oRqPIxSC8PphPw7jI343zmE8HrBa3rHbbTdFp+RwpCXuA3g90deFTwjXilwzct6l3HUoN45jTbjxVIDL5aL/+7//Y/Jau3fvZt0RIEuFGjAtxfAd1dUUUJD4wWog3W7yK/hySn2JQwfLKFhVHZfxZGpRgQB580soGQKHysip4NgDu/ZQ4EB70QJ0lN6I8X83zqB3P5tFmyr/SuoYnY4MCobqHmNubnzFrSNHjmQqM5HnV2gWL2bg3S43pWH/NYyD+UXkC9TdC0dijCPMMDtCIU2M+671efTRs1/S/bPEMy0FSgpLyb9rB4XatmZGqF66o87xqLrIKC2hQEUCakZeLzlcrqiHQz4vBSsqKegT9y6/+eYbuvjii6M+R7nrUG6clVi1ahXLLVixYgW7QaJPuIZTDW48Y+B2u2nQoEH0xx9/sI7pKBDetWsXvf322zR//nz2eJcuXTT5cCoOVaLoLObrgn4vhbblKVLfkfoSVx8ootDOPcwzFdPmFPPe0DE+sGELBdNrEpESngjz91OwWoGQeWW1qIogju2OW8bQKzNnRL3XrjW76KjKLrTGuzxq3G197qe89C208J/vRd+uW5dutG3HNvZzuOHsmt6L1pStFN+XHTqUvvvuu6hziOdwzsUWJxdlXEW97mhP06Y/T8lIJIYiNPlyc6fG9TtuuOEG1h1DbnJ/48H36KNv3qf/zYldFsXYTKw9YVlZGeXmTktqkSFF6a79RHnxt6oLbNlJng71ox5nGbp5+RSSCM1KeZcfPP0Zzf31E/r4m2jhgFjREasAh+Grr74y+jBMgXU/RZ146qmnmIFEC7YtW7aw/qXoZTp16lTWXmj06NGavXd2lrIC5OBhUXd4gbFfW/MljjSEvsPhq6A/2gI/OXIq9es8kO0dhd/q169PC2kurcv4I+o54aa0RVGkdyWKgygzK7pX4B8/rKYf3/5VdIjX7WH7u2J7XgXbCqltdRe2CIi8wUtq16Eduw9//PWH36FO/p6Uk5MTdWvYsCEtWrSI7c+K7dmiuFxcNDzEjH/kMUy5dQYNcgyjqqqqqNuk61+gQSdcUPs7eh3dq87vPOf480XHIYoi9vfi9uabb0pO7FigYIHk8fmoE/Wo/X2bVm+JOke4zX3lWxpIl7PX7N27t84xjOxyM3V29BQ9FvyuAwcOxBVmz20k3wlFjtID0YtTv79moeTxRi/oBO9S7HN8/+nPacPvW0XfR26clcA8+Oijj9LXX39NBQUFKdO7UwxuPGXAlxmFwWDGjBnMWAigi8AJJ5zApL0QytCCgE96AoEnk79jH5vQ4B0HQ0H2s3ATE6eW+xJ7Dnt+Pn/0uL8WrqMmB9uITrjnp4+gmy+Knvhxu6v//5jHIXUs4eD4YxKCkY/2cAP+AHV1HCvq4dZDo26J7Jm0DCSM5EgmKP3yyy9RhjDNkSaZMen1eOm+e+4T/XuLCkpY31vJcyFi4KtKXeQkp6ihXjn/b6q/t2XtuV69dlXtz4NoGJ16VL+oMZjoHhj/IDOQke+Fv+leiWO/75xH6eLT/8UWQ8//8ATt77KFcnNz2a3HSUfTyAuvif59QfQaSGPvi+9N+HF4KvzUlY6lrKysuBYZCA0KRjz81rpLq6jrX7itWbKBVny/Oi7DK3zvxCKRHo+Xbhx1k+h58lb76ShPN9HrEK+HSIKS74KZQY34xIkT6dJLL6XWrVtTKsONpwxLlixhk3/Xrl2pT58+Uc+jbx/QSvBZzMsSeHzYZBrY40I2oXU5uhPtabee2rRrXevx3f5v8T0pqS9x49aN2X26I/qSyMrJpHrZmeKZjNn1qHmr5qLPbfhti+jKdN2STfTI0KdrJzhMfDkNs2MafnacLRtGP9YkV/L1OO40CRGH3AY51KRNo7jOU059aZ25/w1+kn5693fR5/7TZxytWyzd2DszOzqzusep3Sizvnj0oVGzhtSkWWPRjFJnmpPadIqe2IqLD9H3MxeJ/r5nb3hZ8rl1v26i9C0Nao3z3xtX1zHUIy++luJR3WjWuimlZUh7uFJG5pL619D4cQ9GRTcemHknraZfRSMfvc86lm4eeZvksUQmgYHMwxGfjIxoA/7ov/6PFn+0TPRXOdKgTiSuHzkkZxQVbC2UPg6O9QhxJHnhhRcwA4RGjBgh+vz8+fPZ85dffnlSZ7FXr17sJsaBPQdFH/d6vKH9uwtDXq+X3TweT+3PB/YUhjxuj+g4v98fmjBhArsPx+32hObP/FF0TEVpZahoX7Hoc8UFJaFgMCj6XKHEseO9x48dz86dcBswYEDtzzg+MXxeX0iK/G0Foo/j2KrKqyXPodSxS50nsH9XoegYj8sTKisqF30ub+OekM8nfvzVlS7Rx4XjFyMQCLCb1HtJPYdjlHqfyrIq0ecWfvhbKH/HftHn9m7dJ3mMe7fsE328vLQytPWvHXGf94P7itnjwnUefu3v2rg76nHhJnXef/9iObsGxPjtiz9EH8d5LS0qFX2uoqIq5Ha7RZ/7Yc4vIS2YNm0a+84s/rJdKFDQXfJ2w8iG7HXl5eLXZ6JkZmay35uKOPAfow24WUFo9oUXXqD77ruPnn8+Opnjn3/+oRNPPJEp9iQTuj322GPZ/fr16ylViExACi9kTzWpPQ4nUaBTfe+999LiL9vRWadJR0Vuuu8Azf6onMrLy6lBA/UEh7OyslinlVQ0I9bevdahbx3APo0Y2PMBFRUVcRnJSLZv385Cw6kEjCM3kByOOlx/537KzZHehcvfX5M5fMopp4h+71Jp4a4W3HiaAOxTpZrx5HA4yYOa81696mZai3FUOyL39u1cHlBFuPGUQciuFSs2B0i7B0rDIHx1x+Fw1OS6665jN47+8GxbGTp06MDuUacmhvB4x44dtfhsOBwOh2NSuPGUoXfv3uwe7bvEEB5HvSeHw+FwUgduPGXo168fkxZDQs/ff/8d9fxnn33G7lEwzOFwOJzUgRtPGVBkfeedd7Kfx4wZU7vHCVC6smbNGhowYAD17dtX+0+Kw+FwOKaB13nGANJ355xzDtOzhTB8//79KS8vj/0brb20FIbncDgcM4GmEuF61StWrGA1nqeddlrtY9C+RRMEu8OzbRUUAUPjFC3JPvjgA/ryyy9Zejg6UOAiateunT6fFIfD4RgMmmHAcYgk/DG8JhXgnieHw+FwOHHC9zw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw4kTbjw5HA6Hw6H4+H80nAdUQIob1QAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "c = PLSCAN().fit(data, sample_weights=np.random.rand(data.shape[0]))\n", "c.leaf_tree_.plot(leaf_separation=0.15)\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "65f1f3c3", "metadata": {}, "source": [ "Unlike HDBSCAN, PLSCAN uses the unweighted number of neighbors to compute core\n", "distances (for now)." ] } ], "metadata": { "kernelspec": { "display_name": ".venv", "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.14.3" } }, "nbformat": 4, "nbformat_minor": 5 }