{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from matplotlib import pyplot as plt\n", "from scipy.cluster.vq import vq, kmeans, kmeans2, whiten\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ ":36: RuntimeWarning: covariance is not positive-semidefinite.\n", " datasetA = np.random.multivariate_normal(muA, covarianceA, size= sizeA)\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAqh0lEQVR4nO3dcZAc1Z0f8O9vRi2YlRPNqrw+ozFr6SgiBR3WrtnC5JSkLM6HsDF4DQaZ2ImrLlWcq85VEXFt3SqmLOnOd2xlC0PqcpdEzlG+FDJeG8FaWFyEbZGQkk/YK+/KQgeKkYGFkcrIJy12tGM0O/vyx0yPenpe93TPdPd093w/VSppZ2ZnulfSr1//3u/9niilQERE6ZTp9gEQEVF4GOSJiFKMQZ6IKMUY5ImIUoxBnogoxVZ0+wCs3v3ud6t169Z1+zCIiBLl2LFjv1RKDeiei1WQX7duHWZmZrp9GEREiSIirzs9x3QNEVGKMcgTEaUYgzwRUYoxyBMRpRiDPBFRisWquobIi+nZIiYPncKZhRLW5nMY27YBo8OFbh8WUSwxyFOiTM8WsfPJEyiVKwCA4kIJO588AQAM9EQaTNdQokweOlUP8KZSuYLJQ6e6dERE8cYgT4lyZqHk63GiXscgT4myNp/z9ThRr2OQp0QZ27YBOSPb8FjOyGJs24aW3zs9W8SWicNYP34QWyYOY3q2GNZhEsUGJ14pUczJVb/VNZywpV7FIE+JMzpc8B2Y3SZsGeQpzZiuoZ7ACVvqVYEEeRF5VETeEpEXLY/tFpGiiMzVfn0siM8iaofbhC1z9ZRmQY3kvw7gVs3jDyulhmq/ngnos4h8c5qw3bpxADufPIHiQgkKl3P1DPSUFoEEeaXU8wDOB/FeREGxjtAnD53CXTcUUMjnIAAK+RwevPN6PPfyOS6uolQLe+L1CyLybwDMAPiiUuqC/QUich+A+wBgcHAw5MOhXqGrptl/rIgH77y+YaJ1x9Sc9vuZq6e0CHPi9b8AuAbAEICzAB7SvUgptVcpNaKUGhkY0G5RSOSbl/YHD0yfcPx+Lq6itAhtJK+U+oX5ZxH5GoDvhvVZRHZOI/HiQqmeb993dF77GgGaFlex8yUlVWhBXkSuUkqdrX35SQAvur2eepfXANrqddbnMyKoKKX9vJ1PnsBSpQL9s4DC5QVS07NF7D5wEgulcv15LqSiJAkkyIvI4wA+DODdIvImgF0APiwiQ6j+n3kNwB8G8VmULl5XorZ6nf15pwAPoCmNY9ffZ2g/0/4eXEhFSRBIkFdK3at5+K+DeG9KN68rUVu9Tvd8u8zrQ6v35OQsJQHbGlBXeVmJOj1bRLHF65yeb8dCqYx14wdbvo6Ts5QEbGtAXdWqdbCZMnGiAAzteTaMQ3PltfMlUbcxyFNXtWod7CUNY50UjUJG0FRvTxRXolwmqKI2MjKiZmZmun0YFDGzKqa4UEK2VhUjcjk3HjcC4NWJ27p9GER1InJMKTWie44jeeq60eFCfURvVsXENcADzMVTsjDIUywEWR0TJubiKWkY5CkWklCOmJHLZZvsUklJwSBPsZCEFMhyLYVUXCjh/qk51943RHHBIE+xoKuyiTOFau8bjugp7hjkqSvsuzEB1bLEJFFAy77z3HWKuo0llBQ5XU8YQTVoZl0ai8VVf5+BXbdvaqqb151nzsjGosaeXTXTxa2Ekm0NKHK6ShozrCctwAPAhcUyxp44DgC++u0ExW/A9toUjtKB6RqKXBIqafwqV1RT6sZLX55OmQHbzx61XjZUofRgkKfIJaGSph324N2qL08Q2gnYUVx8KD4Y5ClyY9s2wMhKtw8jcPbg3aovTxDaCdhRXHwoPhjkKXKjwwUYmXQFeSMrTcF7dLiAB++8HoV8DgKgkM8FPunaTsCO4uJD8cGJV4rc9GwRi+Xlbh9GoG5c14/JQ6dw/9Rcw+Sn+atTTpOrY9s2aCt43AK2eTysrukNDPIUuTRO8P3w9Pl6hVDQ1SpeqmH8BuygLj4Uf6yTp8itHz/ouIl2mvT3GehbuaLj0fKWicPana8K+RyOjN8cxKFSwrHVMMVKr0zwXVgs+yptdMJqGOoEgzxFxlziX1woIV3Trt60W4vezWoYtmVIPubkKVTWXZ/M1gWo/W5+nZHLHR7Trp3R99aNA3js6Lz28TBxZWw6cCRPobGuxgTQlIdXqOaVf/5g72yl187o+7mXz/l6PChcGZsOgQR5EXlURN4SkRctj60Rke+JyM9qv/cH8VmUHF52ezJHtoUU5emdUlECtFWL3q2cvNP7FxdKTOEkSFDpmq8D+M8A/oflsXEAP1BKTYjIeO3rPw7o8yhGdDXcALQVIXbmyHZs2waMffs4yinI2/zuNWvwk/m3my5wv3vNmnqaw09TsbX5nPZnab0r6KSrpNP3On2u4PLfLVM48RdYCaWIrAPwXaXU79S+PgXgw0qpsyJyFYD/pZRyHcawhDJ5dO10jYwAUm3a5cbadnd6toixJ463/J4kKORz2LpxAPuOzjekqMzzBeCrBXGrlsVeWhpPzxax5+mTuLBYBgDkcwY+vvkqfPf4WSyUyh2fcz5nYNUVnZeLUnvcSijDDPILSqm85fkLSqmmlI2I3AfgPgAYHBy84fXXXw/keCgaTjXcTszJ1v4+A0oBb5fKWJvPYfHSUj0AJZ3AefRdcDlXt7p3XZDefUe1h32rOvpuXEDj0je/V8S6n7xSai+AvUB1JN/lwyGf/OaFH94+BAANQcfPRSIJ1uZzrvlsJ7rvsQd300KpXE+TtMrZTx46FfkdUhh986k9YVbX/KKWpkHt97dC/CzqEj/VIoV8DqPDBex5+mQq0jJGRpq6aZp9Y5x+LllxXiFg/x5zBO50h2MG0lZ19N1aNMXFWvEQZpA/AOBztT9/DsB3Qvws6hJdR0Nd8DMygsVLS1g3fjA1aZnJuzdj8lObtV0mnTo9uu18Za+88TICP7NQctwEvbhQwrqQW0gU8jn09xna53plZXPcBVVC+TiAvwOwQUTeFJF/C2ACwO+LyM8A/H7ta0oZXTtde/DL5wxA4Cm4J62U0gzoZopm8tApTM8WHdsMO51fPmc0pTa8jITX1u6O7roh2rRIzsjike1DODJ+M3bdvomti2MskJy8Uupeh6d+L4j3p3hz6mhoPrZl4rDnCo7iQikxm3mbi4LcVoXqfi66Spjdd2xqep3T5K31+8xAevCnZ9s/EY/MSfOCrXqGrYvjresTr5R+fnOzSQjwAOojdz+bdbsFRPska87IIJsRVDRrB/r7DOy6fVP9/aJIgZkB3loBZK+xf3j7EIN7zDDIU+hajUiTyq2Kxv6422Kl6dkihvY823S3U6ptrGLt+WMtnTQ9MH0imBPywHpe7G2TDAzyFLqxbRuwY2qu24cRKDNVoitvBJpXo+qC4czr53Hwp2dbjsIVnOvOp2eL2uZlYbGel9+7GOoONiij0KXhP3xGqqNo6yQqALytCdD2/V6dguG+o/Oe0yylcgV7nj7Z9LjusSCszErLyVT2uU8GjuQpEkmZTHWyrAAR4NWJyx0zh/Y8C91OtSsy0nBhcwp6fn8aFxbLeGD6BJ57+Vw97RNWLn5FRvDnd15fTzHlayuU75+aw+ShU/WKolY9daj7OJKnSCQ5wJt0q051SrZNyoMMevuOzjfsNhWWxfJyvTx0dc7AhcUyFkrl+ueOPXEcWzcOsHQyATiSp0A5TTAWUjr56sXYtg1NZZPWyVQ/orxUfuZrf6ftpglUm8899ZMirjQy9ed1k8LUfRzJU2Csm4TY9zV1WgGaJAI09FB3Wulpf1y3MOozNw1Wu3XG2JHT5133A7h4qdJwd/POki55Rd3GkTwFxq3awqytto/yza0Bk8AcRZsXr7tuKGDqx280tB4wsoJdt1cXNrXq8e6lsiZJWFkTTwzyFJhW1RZeV4AmQalcweMvvIF7P3R1w0SodWFTqxryhRQFeBMra+KHQZ4C0061hX0F6OqcgUtLFSyW43/rX1EK+48VtfXrXmrI07hIbG0+19EuVRQ8BnkKjG6C0Uu1hW6EPz1bxBe/dTz2VTlm/brX5mLWx9O4SGxh8VLDNo5cBdt9nHilwDh1XvSz1+iWicNYN34Q90/NeQrw+Zx+8jNKFxbLTZtZt+rxDlR/XnE4fq+8TBNfvFRp2qfXaSEXRYMjeQqUU969FXsO28v4vb/PiDSv7Vb2aHaktC4eMjLSEPB0dzW779iUmNF8J/dUFxbLGNrzLEssu4AjeYoFXQ67lds+cFVkqyvzOQMrss5jWTMtYZaPXlgsA5pWCPYANzpc8DRCTgNzy0L7XQ+FiyN5ioV2qjKee/mcdh4gaALn1a2mrEjTMZQrCquuWIG5Xbc0PG6dmMz3GZEucOq2UrmCL37rOADm6KPCkTzFQjsj8jMLpYZ5gLC0CsJu2/rpWg43jfgTrJ31XBWlOKKPEIM8xYJuRawZP5w2vzYvDKPDBRwZv7krWwdmRVwvMuYxmpPKO6bmErcmwEkhn8NX7xnCZ28a9J1yMstJKXwM8hQLusqch7cP4bWJ2/DQPZs9NcLaunEg0vx2zsjioXs2u27cPbZtQ8PoPS0E1Z/36HABXxm9Hg9vH2r4u/vsTYOObR9MXDgVDVExqkMeGRlRMzMz3T4MiqFWC2zs1Tlhs2+/Zz1G6z61hXwOi5eWEp+W0RGg5XZ/busd7FsJUvtE5JhSakT3HCdeKRHsK2Ott/rd6H/Tt3KFtlIGaN7YO60U0LJXjfmcdYEUABgZYUviiDDIUyLoesGMffs4IGhoEBYV+16n5sUn43NzlIxUNyRJKs8pF3serVfqRmOAOXlKBF0dfXlZBR7gvcYi64SqtVrGbxuGJAd4wFtV1OShU01/T+WK4sRrREIfyYvIawB+DaACYMkpb0TkptNJOnO1qluOvGBpf2zm/rduHMD+Y0XHfjztLOJKE13KxT5/4pSy4sRrNKJK12xVSv0yos+iFOq0Y6MZ4I+M36ydpDUDt64tw8j71zRMqFrL/3o5UBmZ6kXu/qm5+mQ40Dwn4dQOgnvBRoPpGkoEXYmikREYtlYDbrtNWfva+2mkZi2RNNMxZhuD1Q4NxmK+6VMgysto2gVs94GTTXc2Cs1pL+4FG50oRvIKwLMiogD8N6XUXuuTInIfgPsAYHBwMILDoSSyV9dYR45ed5uyd4D0s6zeqT/8lUYGOSPbdFdwxYpMy1YIaVMqVxxTV+adFHvMRy+KIL9FKXVGRN4D4Hsi8rJS6nnzyVrQ3wtU6+QjOB5KKKfA7GW3qU5Hjk5pmYXFMh7ePtR0obnfZ2fJK1ZkUr1HKmviuyf0IK+UOlP7/S0ReQrAjQCed/8uSrOwdw5yGvV38hluu17pLj5+ave3XLMGPzx9vu1ji5P+PgO/KS8HeoGlzoQa5EVkFYCMUurXtT/fAuBPwvxMijcve58God2+9k787nrltTtmf5+BH54+n5pOlNdd9Y9w98ggt/+LkbBH8r8F4CmpNphaAeAbSqn/GfJnUox52fs0jvzeHVhf7zSiFwBKdbYZR9wcOX0e6wfexdRMjLB3DUVq/fhBbVATAK9O3Bb14UTC6ZwB992mkiorgtMPfqzbh9FT3HrXsISSIuVl79O0cTq3Qj6XyvOO++brvYZBniLl1pI3Lsze7+vHD2LLxOGON7dwO2fdc0nn1P+fuoMNyihSYVS+BCmMiWEv57z7wMl6XX3SUzj3fujqbh8CWTDIU+SCrnwJUlgTw63O2Vojn+QADwBfGb2+4euwS2bJHYM8kYXToie/PWqaNutWwNulsjbIpanJWT5nYMvEYccGb2GVzJIz5uSJLIKYGNZt1r1QKjf0eLHm+dPS5CwD4OKlpYZ+NvuOzjveGVE0GOSJLIKYGG41MrcHudRU2Gg2cHFKPaV5x6y4YZCnwAVdnRIlvx0qdbyMzK2vGdu2IRUbJfnZAEWARP27SDLm5ClQbtUpQHyraqw6nRj20vt+bT7XkLe/0sigVI5vg7KgK3687A9LwWCQp0A5VafsefpkQ+OqoCbg4li50apvTc7IYuvGgYbXlMrLMDLSsNl1nFhbBbsdob3tspu0zEXEHYM8eeI1mDqNYHXb7VlLE9sJ1mE3O2v3AmKvi9dV1zjtWZv1uRF4VKytgrdMHNb+PZvbJ37xW8c9nUNq5iJijkGeWvITTDPiLzd7ZqHUdrAOs9lZpxcQe8rHesFwa1pWUcpxRB/FIind35914nl6tojFS0tN32fdPtFrL/04rXJOM0689jCvE6RuwdT+fn6zDWvzOc/vbxdUTbtOu8ekYy+pNPc91REA22+8Gn1G839NhWow3XLNmkBbB5jvVcjn8NV7hvDI9qGmiWcAGP6TZ7Fjaq7priyfMxomp72M0Pv7jK6n1XoFR/I9ys9I1WswdQuAq1ZmsXip0jASNUd/TiO/VsHaaYIzI4Lp2WJHQSTIC4juguF0LVQA9h97E79xmIQtlSt47R9KOP3gx1y7W3qRM7KOlUP2uxC3OYZVV6xoeH2rOQkBsOv2TR0cOfnBkXyP8jNS9bpAyC0AXlpabgpIGQHun5pDxmFU2mpE6NTcq6JU04Ijv4Lslun3wlAqN/+sdO/n91gEQM7I1H9/Z6mCHVNzuGbnM3hg+oTj97Wq+7efn1mG2t/XvMm5APjMTYMcxUeIQb5HOQUe3cjY6wIhp6AjgDbHfLE2stdN0nlZgGQGE13qotNVlUF2ywx6gtF8P78dLNfmc3jpTz+Kz9w0iFJ5uZ5aqyiFx47OOwZ6L3dUdqPDBcx++Zam1M/D24eaettQuBjke5RbQLaPgN0WCFnz+ouXlmBkGgNuzsh6TilkRXwvQBodLmDZoZKjk9x8EIuiTEFOMJoXGnMit1SueM7Pmxfwx194Q/v8vhfmtY+7XaRaXfhGhws4Mn4zXp24DUfGb+YIvgu4M1SPmp4t4v6pOW0AtpbLtXoPe+7VyApWrVzRVC7oZRl7u7tDuZX0xWUbuk7z50A1xfLgnR8AAIw9cbyphUAQHtk+1FTSujpn4OKlpabPy+cM7L5jk2PgjuMahrRy2xmKE689anS4gB0+Jjx1/2G1td4VhVVXrMDcrlsaHveyqbV1xOgnQGzdOIB9R+e1k7px4TUc53MG3lladvhZVUfse54+GUqAB4AdU3PYfeBkQ1BfKJVhZAT9fQYWFvWdNO2i2rCdWmOQ72EFh+oU++25039Yp6Ctm4gD0DAy/NVvyg3llkZGGmqxvQaI6dki9h8rNgRRAfDBwdWYPHQK90/NxWIU6fSztsoZWey+o1p1oltQZM4z6BaWBcncvMSqvKzQt3IFZr98i+Y7miV1w/Y0YpDvYbpSN90I2Ok/rNPqTPMi4TQan54tYuyJ41i2jEaXUd0dyay2cQpw9lSC7rUKwA9Pn68H/jiMInU/a11qyzy+dstKw+Tns8Ncw0D+MMj3MK9b8Tn9x6wo1dSrxDox6DQanzx0qindUFlW9RGk05L4omZ1rNNr7Y92exTpd9tDpzUACoAI0I2pND9VQk7Hz1YG0Qs9yIvIrQD+E4AsgP+ulJoI+zPJOy8dF/N9hjZFULDk5u2Ba8vEYcfb9XZHc1mRjnZR6vYo0k93S908gynIAJ/PGdr0jJ3fOQ6vd4kUvlBLKEUkC+AvAXwUwHUA7hWR68L8TArWA9MnHHPAxVoflq0bB7C21qFw8tApTM8WHfPPxdrFoB0VpToK1EkZRermGcKQzxmY21WtZbfX25sTre2WjwZZgkqdCXskfyOAV5RSPwcAEfkmgE8A+PuQP5cCMD1bxL6j+tppU3GhhMcsrzHTMm7NtC6+swQjK74rRLIieO/qK9vaVShJo8go9ny1TvL6TSV5FecN23tJ2IuhCgCsKy/erD1WJyL3iciMiMycO3cu5MMhPyYPnWprNFkqV1y/b6FUBlS1n40fFaV8r/IEqtU2SRpFer1babdJmW5UPTpcwNi2DU13ZJR8YY/kdf8KG/7/K6X2AtgLVBdDhXw85EOYOezyssLSpcbRqjn6d6raKeRz9cC0+8BJT7lk1N7z2zPziVmY42VnqZyRxV03FLD/WNHzqN+tIRnr2tMr7JH8mwCutnz9PgBnQv5MCkgnOWwvY0x7GDd3H3rons3avjFbNw5gy8Rh3D81h1/9xl+t+JHT5xta/XbawCxMursVXY78K6PXN+S9V63M1n/uWRFsuWaN55x4kK2VKV7CHsn/GMC1IrIeQBHApwH8q5A/kwLSqmWsE3OUOfWjN3xvZ3dmoaQdrQsUpn78Rj2P32mFSbdLKt34yZGbeW9zJG7+WCpK4Sfzb7dMU5lrDpzuHLpdkUSdCzXIK6WWROQLAA6hWkL5qFLqZJifScGxBptW6YOsCJaVqgckAPju8bOeUyom693DO0uXe6ovhrDJtRnA3BZtdSvF43fSsp0Vpq36xAPJqUgiZ6HXySulngHwTNifQ+GwjxSdAkJFqYbmVu3cARjZy60N9jx9sq0KEyMDeL0erM3nHHPRM6+fb7hzKC6UMPbEcQDxzFG3s8K0VRVPkiqSyBlXvKZYkCNR8/vcNmm2rmhtJ0CXKwqTh05h5vXzbfdn2X7jIIBqO92Kqm6MfdNv9+Mn829rF+Y4jYC/8cJ801aG5YrCnqdPxjLIt7PC1O0CUIj55DR5xyCfUmFUS5jf5zRK73RFq3mcj7WozXfz3MvncGT85qaNKZwueE49YpymEsJuDtaudlaYOl0Y4tSimTrHTUNSKqhqCftm3wDqGzvrFGtNw7rFeoGxHvvuAyexeGmp6fWrc81b1CVROytMg9z9iuKLI/mUCqILoNPdgLnlnlPaxunxKCgA68YPNj1unQC25t0vagK/m3yMLwp+J2vDWulK8cIgn1JBdAF0uxvoZiAPQqlcqeftvcoA9VYAacHWA+nHdE1KBXEr7nY3UEhBaZ3fC9XqPoMBkRKHQT6lgugC6DTqN2/rjWz3cu9B8Nv7Ja6TrkRuGORTbHS4gCPjN+PVidtwZPxm36NQt7uB0eECVq0MN9t37XtWIaw53JyRxb0fulp7fk4f2W5DMKJuYpAnR63uBt72uZrVr5+9dTG0HZA+OLgaI+9fgyuNy/8F8jkDD955vWMHzaTPQ1Bv4sQruXKbmPPSLTGujpw+j6OvXkDFUhBvtlFw2nS7UFshy2oUShKO5Klt7fR2j5OKbcVTqVzBjqm5+qYmVmYXzJ1PnkhMN0sigCN58sBp9GqOYHc4rBpNqoVSud7ad2GxXD9nP03AOOKnuGCQJ0fTs8WmzTl07RHcFkYlVXlZoW/lCsx++Zb6Y04tEOylptyAg+KE6RrSMgOVrlWwOXo1X5O2AG+yB2+3klIrbsBBccKRPGm16iRZrO0DGvaG092UEcH0bLHePvniO80tEHQLzFq1lGAqh6LEIE9arXrcCJDYyhqvKkrVe9zo9lLt7zOw6/ZNTQHaraUEUzkUNaZrSKtVjxtzw+1WrK9YtTILI5OsBUVmjxvdHcuvSku4f2oOWyYON1TYuC0iYyqHosYgT1peyiMrSrV8jTVbv6yA7TdeXV9clc8ZyCYg6Lt129SVUrotIguiOyiRH0zXkJaX/V0LltLCMwslrM4ZEAEWFsvIaCpuSuUKHjs6j3zOwMO1rQKH9jzrex/YOLKXUjotIguiOyiRH6JiVBkxMjKiZmZmun0YZKPbszVnZF0bnq0fP+jYHsD02ZsGse/ofMvXJYl54QP0fdrb+VkStSIix5RSI7rnOJKnltrZXMJLy4N9R+eR7zNS1d2xvuG3qtbam4/ZJ1dZXUNR4UieQjE9W/S0EjafM1KRrvGCe6dSWNxG8px4pVCMDhfQ39d6q7xeCfBA55Or9v122TOHvAgtyIvIbhEpishc7dfHwvosiqfbPnBVtw8hVjqZXDVz+WyORn6FPZJ/WCk1VPv1TMifRTEyPVvE/mMMQCa/Wy/asb6e2sWJVwqUuWQ/7ath/XBaGesH6+upXWGP5L8gIj8VkUdFpF/3AhG5T0RmRGTm3LlzIR8OhcmaUqDLlKqOxDvJpXttjkZk11F1jYh8H8B7NU99CcBRAL9EddHjnwK4Sin1B27vx+qa5LE229ItgCK9z940iK+MXu/59ayvJzeh1ckrpT7i8QC+BuC7nXwWxY898PgJ8Dkjiw8OrsaR0+fDOrxYe+zoPEbev8ZzgGZ9PbUrtJy8iFyllDpb+/KTAF4M67OoO/y2Gs5INXVhXwFqDVxbNw5g6sdvoFxJ/x2BbkcpN2777RI5CXPi9T+KyBCq6ZrXAPxhiJ9FXeB30m9ZAa9N3NbwmBm4zGC/7+g8VidsgVQGwHIb38dJU4pCaEFeKfWvw3pvigcvrQusBKhvwgE0VuIILnesTFKANytn9jx90nd7Bk6aUhS44pXa5rfuWwH1um57JU5SkzNKVe9Gdt2+qWXbZSsjKx3VzTvhqliyY5Cnto0OF5DPtW5dYGWmKNKydeBCqYwHpk/Uz8fLRir9fQYmP7U58Pw6V8WSDhdDUUd237GpqbTPjZmiSFM++rGj8/U/V5RCBoBkBJXly/cnRlZCCexWbqtiOWHbuziSp45Yd0Fqxbq0P8356GWgIcADiCQfxVWxpMORPHXMHCW2GtGbC3emZ4tYvLQU1eHFQnlZYfeBk45lo8+9fK7j+nfuOkU6HMlTIFrl2Av5XMPOSGnaKMQrM39vz5s/dnQ+kDy62wbi1LsY5CkQbikBa6BJy4Rrux5/4Y2W599ud0m3DcSpdzFdQ4FwShVkRRoCTa/nh722fmj358RVsWTHkTwFwilV8NA9jRUl7eaH/dSgx5mHCksAzKNTcBjkKRBeUwW6i4FOzsg0vZeXCp64y63IeDr/4kIJ1+x8Bg9Mn4jgqCjNmK6hwHhJFdi7KYpUe9rYvbO0rK028VOTH0eL5WX09xn1c8jnDHx881V47uVzTemuilL1Gnx7W2J7hQ47UpKTjvrJB4395NPDaxBaP36wZQl5zsjirhsK9UCYrfWtz+cMXLy0lOiOldae8NfsfEabs8+K4PSDl7dIZm95sgutnzyRjj0ImWWBAJqCkJcmZ6VyBfuOztcvBhWlkDOy2H3HJgCt7wrizLoi1WlStqJUy81ZuLKVnDAnT4Hzs+m01xy9PfxZg9rYtg1Ym88lLsCbzEoap743ImiorXe6GHDbRdJhkKfA+Vleb5+w9dLgy/p+adhXdm0+h+nZIq5YoT/33IqMp3kIs5UzkRWDPAXO76bTo8MFHBm/Ga9O3IaH7tkMr2FeAdgxNRfJRKzfbpte5Ywstm4cwM4nT2Cx3Lj1iEh1L9hS2duWJNZWzkQmBnkKXCfL60eHC6H08jIy3u8QdHzcYPiSkWoXS92Fau3qHL4yer2vmvleX2xGzRjkKXCdLq93qof3k8qxyucMvOvKzmoMFhbLoYzmL15yvgsxA7buoun0k+AiKrJjkKdQWFMwR8Zv9lX14RTU7v3Q1Z5TOSYjI7h4aamhIZqR9X+xWJvPYfcd/nZ/6pQZsHUXzc/cNMhmZOQJSygpdkaHC5h5/XxD2aQCsP9YEfk+w7WDZX+fAaWAt0tlrM3nsGgL8ABQrij09xn4f79ZQtlDSY4ZPM0L1Y6puTbPzDt7wNYtNBt5/xouiKKWuBiKYmnLxGFtxYzTAigjI5i8u3nnJbfFVo9sH6pvJO7E3Kjb+r5Ox9ZnZFBeVh0vzioE2GOeeoPbYiimayiWnCYQ3y6VMfmpzejvu5wfz+cMbYAHgHyfex59bNsG11x/38oVnvrv5Iws/vzOD2DyU5tdP89Nzsjike1DGNu2AfuPFblXKwWC6RqKJaeVsKtzhq92um43ql7SLk61/QAcUyWt7g50Cpb32DJxmHu1UmA6GsmLyN0iclJElkVkxPbcThF5RUROici2zg6Tes3Ytg3asseLl5Z8jWjfLnW2A9Vqh4oa68Ty2LYNmDx0CuvHD2LLxGFs3Tjga4K2kM81TE5zr1YKUqfpmhcB3AngeeuDInIdgE8D2ATgVgB/JSLpaAhOkRgdLmjLHssV5WvBT6clhfZMzvRsEVsmDtcDum47v/3HirjrBu8jbnvw9ruYjMhNR0FeKfWSUkr3P+4TAL6plHpHKfUqgFcA3NjJZ1HvWXCoovEzonW6I2jnGKwtFMyAvk+zkKlUruCxo/Oe6/rtwZt7tVKQwpp4LQB4w/L1m7XHmojIfSIyIyIz586dC+lwKImCGNGODhcweffmthcyWT9L13jNrY7Gy1Z/uuDNvVopSC0nXkXk+wDeq3nqS0qp7zh9m+Yx7b94pdReAHuBagllq+Oh3jG2bYO2b7rfEW27G45YP2t6tth2E7SsCJaVwlofpZHcq5WC0jLIK6U+0sb7vgngasvX7wNwpo33oR7WqorFD90o3MqtNt1M0zgRuI/ol5XCqxO3+T5moiCEVUJ5AMA3ROSrANYCuBbAj0L6LEqxoEa0Tnl8ARoCsHVzDnOC1+0CYd+1SocTptRNHQV5EfkkgL8AMADgoIjMKaW2KaVOisi3APw9gCUAf6SUSu7GnJR4TnX31gDstKOV2x2ANVf+wPSJhlYMACdMqfs6CvJKqacAPOXw3J8B+LNO3p+oHbr9Zb3k9512tMpqttsDGrtiTs8Wsf9YsSHAC4C7bmBunbqLbQ0oVXRljmY+vVXFilNKx9xTVve42W7AqfLmuZdZMUbdxbYGlCpu+8u2annslNLJijimbMz35ipViiuO5ClVOgm2TpuKt6p3N9NCOpx0pW5jkKdU6STYtrupuJn3D3OVqr2dAjtSklcM8pQqnQZba+OxZR8rVsNcpeo0z8BAT14wJ0+pEuQCKrccvbmC1freYa1SdZtnYOUOtcIgT6kTVLB1KruMuo8MJ3WpE0zXEDmIS6MwTupSJziSJ3IRh0ZhQTVqo97EIE8Uc0HOM1DvYZAnSoA43FFQMjHIEwVM1zvHHqC9vIYoCAzyRAFy6mQJXE67eHkNUVBYXUMUILeadj+vIQoKgzxRgLzUtLPunaLEIE8UIC817ax7pygxyBMFyEvvnLCbmRFZceKVKEBeatpZ905REuWh015URkZG1MzMTLcPg4goUUTkmFJqRPcc0zVERCnGIE9ElGLMyRN1CVe9UhQ6GsmLyN0iclJElkVkxPL4OhEpichc7dd/7fxQidKDuz1RVDpN17wI4E4Az2ueO62UGqr9+nyHn0OUKlz1SlHpKF2jlHoJAMTjhsdEVMVVrxSVMCde14vIrIj8bxH5F04vEpH7RGRGRGbOnTsX4uEQxQdXvVJUWgZ5Efm+iLyo+fUJl287C2BQKTUM4N8D+IaI/GPdC5VSe5VSI0qpkYGBgfbOgihhuOqVotIyXaOU+ojfN1VKvQPgndqfj4nIaQD/BABXOhGBq14pOqGUUIrIAIDzSqmKiPw2gGsB/DyMzyJKKu72RFHotITykyLyJoB/BuCgiByqPfUvAfxURI4DeALA55VS5zs7VCIi8qvT6pqnADyleXw/gP2dvDcREXWObQ2IiFKMQZ6IKMUY5ImIUixW/eRF5ByA17t9HCF4N4BfdvsguqiXz5/n3ruiPP/3K6W0C41iFeTTSkRmnBr694JePn+ee2+eOxCf82e6hogoxRjkiYhSjEE+Gnu7fQBd1svnz3PvXbE4f+bkiYhSjCN5IqIUY5AnIkoxBvkQicikiLwsIj8VkadEJG95bqeIvCIip0RkWxcPMxRO+//Wnkv1uZtE5NbaOb4iIuPdPp4wicijIvKWiLxoeWyNiHxPRH5W+72/m8cYFhG5WkSeE5GXav/m/13t8VicP4N8uL4H4HeUUh8A8H8B7AQAEbkOwKcBbAJwK4C/EpGs47skk3b/3x45d9TO6S8BfBTAdQDurZ17Wn0d1b9Pq3EAP1BKXQvgB7Wv02gJwBeVUv8UwE0A/qj2dx2L82eQD5FS6lml1FLty6MA3lf78ycAfFMp9Y5S6lUArwC4sRvHGBal1EtKKd2u1Kk/95obAbyilPq5UuoSgG+ieu6ppJR6HoC9nfgnAPxN7c9/A2A0ymOKilLqrFLqJ7U//xrASwAKiMn5M8hH5w8A/G3tzwUAb1iee7P2WC/olXPvlfN081tKqbNANRACeE+Xjyd0IrIOwDCAFxCT8w9lZ6heIiLfB/BezVNfUkp9p/aaL6F6S7fP/DbN6xNXy+rl3HXfpnkscefuQa+cJ9WIyLtQ3Udjh1LqVyK6fwLRY5DvUKs9cEXkcwA+DuD31OVFCW8CuNrysvcBOBPOEYannf1/kZJz96BXztPNL0TkKqXUWRG5CsBb3T6gsIiIgWqA36eUerL2cCzOn+maEInIrQD+GMAdSqlFy1MHAHxaRK4QkfWo7oH7o24cYxf0yrn/GMC1IrJeRFaiOtl8oMvHFLUDAD5X+/PnADjd3SWaVIfsfw3gJaXUVy1PxeL8ueI1RCLyCoArAPxD7aGjSqnP1577Eqp5+iVUb+/+Vv8uySQinwTwFwAGACwAmFNKbas9l+pzN4nIxwA8AiAL4FGl1J9194jCIyKPA/gwqu11fwFgF4BpAN8CMAhgHsDdadzrWUT+OYD/A+AEgOXaw/8B1bx818+fQZ6IKMWYriEiSjEGeSKiFGOQJyJKMQZ5IqIUY5AnIkoxBnkiohRjkCciSrH/DzVYhbfMriwxAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Data generation\n", "# if we want to repeat our results\n", "np.random.seed(1234) \n", "\n", "# First we set mean for datasetA and datasetB\n", "muA = [-5,-5]\n", "muB = [5,5]\n", "muC = [-5,10]\n", "\n", "# Then we set standard deviation for datasetA and datasetB\n", "sigmaAx = 10.0\n", "sigmaAy = 3.0\n", "# covariance\n", "sigmaAxy = -10.0\n", "# covariance matrix\n", "covarianceA = [[sigmaAx, sigmaAxy], [sigmaAxy, sigmaAy]]\n", "\n", "sigmaBx = 5.0\n", "sigmaBy = 5.0\n", "# covariance \n", "sigmaBxy = 0.0\n", "# covariance matrix\n", "covarianceB = [[sigmaBx, sigmaBxy], [sigmaBxy, sigmaBy]]\n", "\n", "sigmaCx = 4.0\n", "sigmaCy = 4.0\n", "# covariance \n", "sigmaCxy = 0.0\n", "# covariance matrix\n", "covarianceC = [[sigmaCx, sigmaCxy], [sigmaCxy, sigmaCy]]\n", "\n", "sizeA = 500\n", "sizeB = 500\n", "sizeC = 500\n", "\n", "datasetA = np.random.multivariate_normal(muA, covarianceA, size= sizeA)\n", "datasetB = np.random.multivariate_normal(muB, covarianceB, size= sizeB)\n", "datasetC = np.random.multivariate_normal(muC, covarianceC, size= sizeC)\n", "\n", "labelsA = [\"A\" for i in range(sizeA)]\n", "labelsB = [\"B\" for i in range(sizeB)]\n", "labelsC = [\"C\" for i in range(sizeC)]\n", "labelsABC = labelsA + labelsB + labelsC\n", "\n", "dataT = np.concatenate((datasetA, datasetB, datasetC),)\n", "\n", "plt.scatter(dataT[:,0], dataT[:,1])\n", "plt.axis('equal')\n", "plt.show()\n", "\n" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[ 0.00155107 1.15002859]\n", " [-0.96014479 -0.72213698]]\n" ] } ], "source": [ "# Compute k-means\n", "# to learn other options: \n", "# https://docs.scipy.org/doc/scipy/reference/cluster.vq.html#\n", "\n", "dataTW = whiten(dataT)\n", "\n", "codebook, distortion = kmeans(dataTW, 2)\n", "\n", "print(codebook)\n" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEICAYAAABCnX+uAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAA/k0lEQVR4nO2df5QbZ3nvv4+0Y1triGWTpWAljkPg2sUYe5stuA2X4kBxaEjYxiQuTS79cXpTei7nNG7u0g3JiW0IjYubJj0XenpzC4fSuLBJDFsHh9pw4hYaamDdXWPc2D1AsB05hQVbDnhlW6t97h/SyKPR+868MxppRtLzOccnWWn0zjsazfM+7/OTmBmCIAhC55KKewKCIAhCc4ggFwRB6HBEkAuCIHQ4IsgFQRA6HBHkgiAIHY4IckEQhA5HBLkgeEBEbyOiF+KeBwAQ0ZeJ6HfinoeQPESQC4Egoh8S0Tscf/8WEZ0hol+Lc17NQERvIqKniahARKeJ6FtE9HsRn6PuewsDM7+Lmf8uqjkJ3YMIciE0Ve3wkwBuZOZ/iXs+YSCiXwHwDIB/AfBaAK8A8EcA3hXnvJxQBXlWBS3y4xBCQUR3AngIwAZm/kb1teVExET0e0R0sqqpf4CIfpmIvlPVeD/hGuf3iei56rF7iegqx3t/VR3nJSI6SET/3fHeViJ6nIg+S0Q/I6IjRDTkeP9PiShffe8YEb1dcyk7APwdM/85M/+EKxxk5ts0181E9FrH358hogeq/385EX3Jodl/nYhSRPT3AJYBeIqIfk5EH6oev46IvlE9/hARvc0x7j8T0ceI6FkAMwBeU33tD6rv/y4R/SsR/UX1u3ueiN7l+PzVRPS16vV/lYg+SUSPed1ToYNhZvkn/4z/AfghgF0AfgRgjeu95QAYwN8AWADgnQDOAxgH8EoAOQA/BvBr1eOHAXwPwC8C6ANwH4BvOMa7AxUNuQ/A3QD+C8CC6ntbq2P/BoA0gAcBHKi+twLASQBLHfO6RnEt/QDKANZ7XO/bALzg+JsBvNbx92cAPFD9/wer125V//13AOT43t7h+FwOwE+r808B+PXq3wPV9/8ZwAkAq6rXb1Vf+4Pq+78LoATgf1av/48AnHKc798A/AWAeQDeAuAlAI/F/fuRf635Jxq5EIZfB3AAwGHN+x9l5vPMvA/AOQCfY+YfM3MewNcBDFaP+0MADzLzc8w8C+DPAKy1tXJmfoyZf8rMs8z8EID5qAhpm39l5qeZuQzg7wGsqb5erh77eiKymPmHzPx9xTwXoyJEXwz3NTRQAvBqAFcxc4mZv87MumJGdwB4ujr/OWb+CoAJVAS7zWeY+Uj1+kuKMY4z8/+rXv/fVc/9C0S0DMAvA7ifmS8y878C2B3RNQoJRAS5EIYPAPhvAP6WiEjx/o8c/19U/P2y6v9fBeCvqqaFAoDTAAgVbRVEdHfV7HK2+v4iAJc7xvovx//PAFhARH3M/D0Ad6Gitf+YiD5PREsV8zwDYA4VARgFO1DZYewjoh8Q0ajHsVcBuNW+9ur1vcU1l5M+56tdPzPPVP/3ZQCWAjjteM1kLKGDEUEuhOHHAN6Oiungr5sY5ySAP2TmrONfhpm/UbWH/ymA2wAsZuYsgLOoCHpfmPkfmPktqAhMBvDnimNmUDFBbAww5xlUTDI2r3KM9zNmvpuZXwPgJgB/4rDNuzXzkwD+3nXtC5l5u3OKAebl5EUAS4jIOc8rQ44ldAAiyIVQMPMpANcDuIGIHg45zN8AuIeIVgEAES0iolur770cwCyAaQB9RHQ/gMtMBiWiFUR0PRHNR8WOXkTF3KLiQwB+l4hGiOgV1c+vIaLPa46fAvDbRJQmohsA1MIuiejdRPTa6i7lpeo57fP+CMBrHOM8BuAmItpQHWsBVWLWrzC5Ri+Y+TgqZpqtRDSvGplzU7PjCslFBLkQGmY+iYowfy8RPRji819ERVP+PBG9BOC7uBT2txfAlwH8J4DjqAhkU/PAfADbAfwEFfPDKwF8WDOHb1Sv4XoAPyCi0wAeBfC0Zuw/RkUoFgDcjooj1+Z1AL4K4OeoaPp/zcz/XH3vQQD3Vc0o/7v63b2nOq/p6rWNILpn8nYAv4KKA/UBAGMALkQ0tpAwbA+3IAhdDBGNATjKzFvinosQPaKRC0IXUo3dv6Yax34DKtr/eMzTElpEX9wTEAShJbwKwBdQicN/AcAfMfNkvFMSWoWYVgRBEDocMa0IgiB0OLGYVi6//HJevnx5HKcWBEHoWA4ePPgTZh5wvx6LIF++fDkmJibiOLUgCELHQkTHVa+LaUUQBKHDEUEuCILQ4YggFwRB6HBEkAuCIHQ4IsgFQRA6HMnsFEIzPpnHjr3HcKpQxNJsBiMbVmB4MBf3tASh5xBBLoRifDKPe75wGMVSpUprvlDEPV+oNAwSYS4I7UVMK0Ioduw9VhPiNsVSGTv2HotpRoLQu4ggF0JxqlAM9LogCK1DBLkQiqXZTKDXBUFoHSLIhVCMbFiBjJWuey1jpTGyYYXmE+EYn8zjuu3P4OrRPbhu+zMYn8xHOr4gdAPi7BRCYTs0Wxm1Ig5VQTBDBLkQmuHBXEsFqpdDVQS5IFxCTCtCYhGHqiCYIRq5kFiWZjPIK4S2l0NVkpSEXkQ0ciGxBHWo2jb1fKEIxiWbujhIhW5HNHIhEXhp0qYattjUhV5FBLkQO37RKSZCeHwyrzTDAGJTF7ofEeRC7DSjSY9P5rHtqSM4M1PSHiNJSkK3IzZyIXZ0GnO+UPRMArI1eS8hHiZJSZKQhE5DNPIuJ8oojrBj+X1OF50CeCcBffgL30GxNOd57gdvWW2k1dvzW5SxcO7iLEpl9j2/ICQF0ci7mCijOMKOZfI5VXSKE2dVRVtbXj66BzM+Qnxxv2UkxJ3zKxRLNSGuOr8gJBER5F1MlKVmw45l8rnhwRwevGU1ch627FOFYp3QNYHZ/xjV/HTnF4SkIqaVLqbZzEinyUEnE/3G8jKZDH5kH5iBs8USFmUsEOnHYQB3P34IZRPpXOVssWI7t68jXygiTYQyc+2/pojDVEgyIsi7mDCZkTbukECvc+hs4H5mF6eTslDUOyxtgghe59xGnjxUM5fYYwQZqxVVHQUhSpo2rRDRlUS0n4ieI6IjRPTHUUxMaJ5mSs2amBwyVhrrVw402MA3j03hvvHDsdqV7evc9tSRBpu3CSkCCEAumzFymApCnEShkc8CuJuZ/52IXg7gIBF9hZn/I4KxhSbQZUYCwHXbn8GpQhHZfgvnS+Va9AcBWjOKk8X9FrbctEop8BnAzgMnjMZpBdmMha03r8LwYA53jU2FGoMZeH77jdFOTBBaRNOCnJlfBPBi9f9/RkTPAcgBEEGeANyZkW6TiTsG21T4nq8Kfp2NnIHAdmhT/BabQrGEieOnm9KixSYudBKRRq0Q0XIAgwC+qXjvTiKaIKKJ6enpKE8rBMA0SsOPYqmMu8amkPLwUJaZPcMKw5Ltt/DIprWexzx24ATuGz+MbMYKPL7YxIVOIzJBTkQvA7ALwF3M/JL7fWZ+lJmHmHloYGAgqtMKAYk6jM5L47bty2GEqReFmRKGB3O+4z524ISREzWFiqnItolvvDaHHXuPSWan0DFEIsiJyEJFiO9k5i9EMabQGtplMnBqtQvnRxscZV/Du9e82vgz5PqvkzlUTExLsxmsXzmAXQfzSuetICSVKKJWCMCnADzHzH/Z/JSEVuKXRdkM6aqZJU2EjddW7NNBEnhMcC4Q+4+am+gYQL+V8rSt5wtF7DxwQuu8Fc1cSCpRqErXAfgfAA4T0VT1tQ8z89MRjC0EwKQWijOSJUoBS6iP0d51MI8vHXoxEnu8E2dWaND5+6X0A3onKgMtqWsuHY2EKCBuQVSBH0NDQzwxMdH283YzqgSejJXGxmtz2H90Wikortv+TCTC3DRkMUriOCdQCW207e52CGYzRchU90zi1ivIItcIER1k5iH365LZ2SXoapo447mdlfwAYObibNPnjUugxhWj7nSenpkpYeTJQwDCVUbs1I5G7RCwfs1GhHqkaFaX4BXP7aRYKmPr7iMYefKQZx1vUxhqByIALJyX1r7XTryKcTVLqcyhM1ibrYUTB+3qixplwbdeQAR5lxAkGkVVqrUZdCMVS+XYNGebNBFOFYpIGa4oKQBWOtjyE1bw6u5ZkpOR2iVgO3GRixMR5F3CyIYVgQVQq5ljGAvQVlFmBlfnYsIcgIXz+gJp8WEFbzO1cOKiXQK2Exe5OBFB3kXMRqhlR4WpAI2KjNX8T/pssYRnR683OtZKU2jB66zD3ikFutolYDtxkYsTcXZ2AbbdMnlivP1cnGXkPFrHmbAoYxnbfDf98pU1wRvGCeiuhdNugs55ZMMKZaRN1AJWV/AtyYtcnEj4YRcQVRihUGHhvDTmGEYx8LlsBs+OXt+RoYRh5yxhgfEh4YddjDiAouXcRfMkJnsB3br7SKBQwiQIw7Dhj3HvIoRGRJB3AV5d6IXWs3x0j/Y95yLrbDnnjL+PK0ZaIkO6BxHkHYhbm7MLPTm1K6va4ibKMEMhOLYT0G3GUMX3tzsRqJlWgJ1EEnY/rUYEecLxE9r5QhG7DuaVqfgAlBqg0D6c98HP5t5uTXj9ygE8duCE8vVuoVcyREWQJxjVj1DVQq1YKmP/0WllyJwzoiJs27NeQbfYEYAFVqrWDs+Uxf1W7fs3EdLt1oR11SODVJVMOp1aBiEoEkeeYHT9MFX4CYrhwVxLU9W7Aa/KhwusdKCHJWOlseWmVbW//YR0HDHSvWAjN73G8ck8rtv+TMc2ExGNvI342erc7wdxYJpocyMbVmDkiUMotTtLpws4M1OClSZcNq8PZ4slLMpY+NmFWZQV36Wz+bONKv7axq7frtMQW2XjbdZGHoftOeg5Ta6xG8wvIsjbhN+PRfW+11bf+XogbS5ZWfwdRanMWDi/D1Nb3gkAWLttn7KV3ML5fb514J330K7fPnTVkobPtVLINJPc0+y8xifz2PbUkVrhNtXiZx9nC+5FGQvnLs7WHPjuap5bdx9pKDGsukaqfva67c9gZMMKI/NL0h2mYlppE37FhnRmFLfczVhp3L5uWUNaNwDfreGOvcckiqVJnFvys5p+oLrt/PBgDs+OXo9cNqONWnEyPpnH3Y8falmRKmeJAKCyM7DH9jMtNFM8a3wy31B9s1AsYeSJQ3XndVdaVBV7K5bK2PbUEYw8cUhbYth5jaqwT93O176P7ar42AyikbcJP1udVxnaXDbjqQmotKPNY1OYOH4aDwyv9p2DYI5zSx7WNGFit7Xvqa65dbP30q0RA5c6PJlo183Y13UKRWmOa1qwvYh5Nfe20ZVjtksMPzt6PYYHc8oM6GKpjDSR8jz2fewEh6lo5G3Cr9iQ7n07Bfz57TfWfpBudNq8u8+k7hxibTHDbXYIW9jJpPCUX7hiMxEuKo3YjZ923UzxLC9hf6pQ9F3EguA8l+68ZWbP+9gJTmER5G3C76E3FQpO7/rabfuw6v5/0m4N7T6TfnP41WuW1Bon9xqve+VCo4VMVZkwbPVCk3vtJSSajXAxNbF5zUF1DbbtefAj+7B22z4sH92Da+55Gstd5r5FGUs77tJsxijm3hT3DkqFfd9097ETSupK0aw2EjRqRfW+LvLBi0c2rdU6bVRZob2ElSJsetOV+Nw3T2o1wIXz0jjykRt8xwriEPM7VlcILU2Eh25b09SW/urRPUbJYfZuUIez5IAJdg/ZsW+f1C4ki/utSDpXOck5EuTCFglLSkE0XdEsEeQdRNgqh15NmKVyYiVi4myxpBVu2YxVi1TREfXD3krhYXLPg5xr8CP7IhW+rchCtq8HCFcaNylRK1L9sAsIa5PzasKcJDtfXBSKJc8a5rroFCdRO8TC1uM2Cesb2bACI08e0mrFuQDn2rH3WOQatB2tFUaYZ6w0FliphjnZ90LnZ9LhFuAPO3a3SUIEeQexKGMp45ZN0IW7SeXECiMbVmDz2JRSeDTjwDNdKHUaX5C64OtXDjSYLQrFEv5kbArbnjqCwkypNvaO964xiuPWndcd8RI17mitF88WfbtN2QvQZk0piqBKSyclCokg7yCi9keeKhTx8Ka1oezu3US/lcLwYA4Tx0831LLxcyzawlQnY0wWgSACw6sUrqoOD1DpQ2oLXXvsB29Zjcn7vc1Furl6afNR4bbPe5UKdpuBdHb7oM7JTgg7tJGolQ6iELEGtDSbqUVe9GrUCgDMlOZw3/hhPDC8Gg9vWmschTI+mcfIE4e0OxrT6BLT5BpnYgrQuMsyFa3FUhl3jU2Fqimy7akjoYR4kN+X6nvLaiJdCGi4R1H1++yEsEMb0cg7CC8ziJ25ZmomsVKEmYuzWD66p5YQ0S2lbsNcx84DJ2op8qba1tbdR7R1a0ztzIC5wIgyLA/Qa/5ejr0w5hQrRZjXlzLqvLS43wIzsHlsCjv2HqsJX51JMVPdTTnx8i8EcVp2Ur120cg7CC+N4lShqNREdDi323bYXTcI8WzGwvPbbwxc6dEdc2+Cl78iiFPNNE65FZqgW/NvRTp6aY6N2+edL82hUI0gyheKGHniUC3VXsWMprSwXQ7h+e031uqpLB/dg7vGpoyvLSrNvh2III+RoKUzhwdz2i2m00ximwa8trOqqn02nWxkKRRLuHp0TygHbsu3zDt3AsuXA6lU5b87dwIwFxhhNEErBVhp7zvqvG6dmedPHp/C+GS+5b8N97lLc+xrylk+ugfLR/dg8CP7Gp4htzlKdT6VCeu67c9g89gUFlgpZDNWoISvOBDTSkyE9YhvvXmVZ8U6p2kgbDOJTtfMw87f2ZbNZPutS15Z3K9YbHfuBO68E5iZqfx9/HjlbwDDt98OwD/U0KsUro5XXpapaaQ6YZYiwtWjezxNd3MMjDx5KNG/DWehLKd5xbQzkyoa58xMCRkrndiwQxtJCIoJXVKGXzYdECw5QVdq1Qu/BJlu5Y51yzB01RLjRBxVBIeVJux4ryLzcvnyivB28cJlA9h0z+cDJaaYFpMCKrur57ff2LZokySRNQzXtf0ZXoukyXPZDiQhKGE04xEP5JBTaPBejZkzVhpE8Wnl89KE8lzFbp8mwmsG+vG9H59ry3z2H53G/qPTStPC3Y8fwuaxKWWMt9GieqKxNyYALH3pJ4Hik4cHc9o4aeX4jgp+SRTiVoo8G50000TcRIgToK1J7iSJkSpORJDHRLs84jphY7+WLxRrUSt+CRXt4GKZa4kgr1q0ANM/u9i2RcXrYdWVeDVeVJctU2rkpy67HECw+GTTJC5bSAHJFUReQjxNhB23rsHE8dOetXD88GrQ8qvXLDGqF+P1XCYhfV8EeUw0050lKDphowo52zw2hZSmPrMXVorwsgV9ODNTQorgm4Wnw66gB5iHUkaF/bD6ndfW0IEAGX4f+1i9jRzATN98fPyt76/9bSpsTW3lt69bVlfBr5MyeJ21UXYdzDdV0tbOEnUrLctfkcE3vn/aV1HIWGmsXzmA67Y/owxnTEL2ZyQ2ciL6NIB3A/gxM7/B73ixkVeIayWPsgIiAQ1z19n/CUDWo7pdq+PYdQ0EgHrBYepQDFzEaudO4N57MXf8BE5ddjk+/tb3Y/eq9bW3g9hh/WzlGSuF5z76rrrjOy2DN5fNYObibNOlAFTf6/hkXluSwUk2Y+Hda16tfD7smHeVCSeKKpUqWlr9kIjeCuDnAD4rgjzZqB5onQD1EnyAXvB4lUnVOaBaUb7UScZKewoxValfEw1WJySCliMOU9nQLyrpEVekxX3jh/HYAbWtvpuZ35fCxdm5QFU/cwGO1WE/V0GSw3zH1AjySOLImflrAE5HMZbQWnTdhFTMMeORTWuRUgQPW2nSmoG8OhG5hfjifguPbFqLyfvfGTiJxx7T5JiN1+a04+eqMfg2w4M54+QqtznEJKEmbEMKN155BUBjgtP+o9OBxu8WLszONdwLLzOWvTjb9yOsf8FdbbSVPT7blhBERHcS0QQRTUxP9+YPKgkE+VHa3VpU9u6F8xo7xdvouseoFox+xzhBMlOBS42o/WAAjx04gZmLs5UoCAd2qQJ3UpZpOrx70TKtm+LMPPTKAtUljdmve0VmuO91Uh2e7cRZ9VOF00FsE0UAQjP1bUxomyBn5keZeYiZhwYGBtp1WsGFad9O2/Gqe/i9anSrNE6d1p8vFGuCasfeY9h4rblWOr8vhaGrlhgff2amBBBqmXrZjAVQ5fUgGpuNyjkdZaElnXZ/3/hhz2xFGzvRxxYeSawREge6chaEegexjZ+CsbjfMi4K1irtPLKEICJaDuBLYiNPNjr7bNAOQkETJLwcoO6ysSmCtjaH6vj5falASU+2I0pnB/cqQJYmwhyz1jkdZaLXuQuzWkdamKiiMnPoaKJuIpuxsHB+nzL01mtntHX3kYb7EcZJDoRPMJKEIAFA8M4zUYVJjmxYgZEnDjXEDasaXqiw0oSF8/oaHqRiqYwFVipQxEuZ2fOhyxeKuGPdsoZIBROHZNjvSxXG5jX/oHjFa/cSVopw7uKlBdL+Lmcuznp+zg7h9XNkq2rFq4jazBVV1MrnALwNwOUAfgRgCzN/Sne8aOSdRRRhks2miLeibIBfOKJul+KHKrzTb5wgkRG6eXdLGeJW4Rf+GmVDZb/Ip6g1cqm1IhhjKtBVxwXptq7CjlcPOoafcPMKS4yivoZpqKFpZ3t7gQkb9y94E3VNFV247+3rluGB4dWBx2tp+KHQ/ZjWqdYdZyqAdS4je0EIE9Wic0TZYX86wvR4dEeYmEax6ByRi/utOqexvUsolsq16+rl7k5RE7XJY3gwh43X5up+14xKtmqUDk8R5IIRpgJJd5xO2PRbqTqB9KvXLNHW5nZGw3jhjM1+YHg1HrptjeeYuvGCRHkEXcBOFYp1gv/chVll3fALpTLOXajYb89dmMXYt0/WxiwzI2Olm0pfF+ppRWTP/qPT2ubnUSHOTsEI07A63XG20KmrwpgmlMpcE0RlZvz7ibOetmnb6RQkOsTPwRuFQ9drAVMJ2my/Vef8LRRLSAFYOC9dF7EzU5qrdcFRRbB4nUMIRhgnvom5sR29P0WQC0aYVmvUHZdz2Mq9wuuKpTL2H532tVMGFb5eVQpVgn79yoFaETETR2eQBYyg7n05B2AmhN27m/qtxsXifgtbbloVqK+nacGsdlQ6FdOKYIRpOzKv49zZjLqkItOa7FGkuTvHc/Z43HUwH6hvpe6htOdlm2/8BG4YxTqXzRhluApq7li3DJP3v7OumqHJvTc1N7aj96do5IIRpvHnQeLUm9VUgjTYCILXAxom3t7PHNQMznPs+c6LLS081q2Mfesk9nznRRRmSsoSzrp7b2oyCZq7EQYR5IIxpoLT9Lh21mQPQhibpsnDGoVN1K77XpgpNdTFFjN5OEpzXFsAdb4G1b0Looi0SumwEUHeBSShQ0kY2qGphCHsTsHvYW2muYOq7rt93+8amxIbeYtR3fskKSKSENThRFXbWrhEq75Tr1rwXo2CdTXPO61RRNSYNleOCqdD1KbdSpTUWulSwthzBW9atVPwGldXklZVVtUeoxeFuHtRW7ttX9uE+ZmZEkaerG/x12qTiSkiyDucdsSo9iKtekB14+ruF0Pd+7EX76/dzCRIB6eoKZU5kUqShB92ODq7rdSe7iy8wheDHN+tLO63sOO9azBx/DQ2j02FFuKPbFqLxf36rkomJHERFUHe4bQjRrVb0HXbSQJB76OuMQJQEf53rFsWqnVe0shmLrUCBICdB06EdurOSxOGB3NNR/ckcREV00qHk9TIj6RhmoUXF0Hvo/N4VYME4FKPzk6OaCkUL9mld+w91tR1MCq/A6/uVn549aqNE4laEXqCqDodJQ2JXglGs3VpHtm0NtDCH3VUi0StCD1NUpzCqgcbCL+j6tXolbDohLiVIszrS2lbDAKVRT9I6GE7d4EiyIWeoB2Fi/xQPdgjTx4C+FIrtqAPexIdb0kgiDnJaY7SdbFy+ytMhHQ7Q4PF2Sn0BElwCqse7FKZG/ppBqlVnUTHWxLoS0FZ390NAXh29HoMD+awY+8xpRBPEzUkg5kUzNItsvlCMXJnuwjyhJDkiIpuIOpqiWEIoj2bHhu0a1KvUJoDwJWwRYK+i5JzIdSFNJaZQxXM8lpkTSpqBkFMKwkgjC2tU+urxEncWXhBaq3ohIDqvj94y2pse+pIz1U+9HNcluYY/fP6MHn/O7VlF5w7Mt14qkVAdy9TRLh6dE+tpr1Xb9UozSyikScA07rGNkFqJrcb2VnoUWnPVppgpeoFhc7ko7vvADB5/zvxyKa1tR2HCblspunkmDiZq4ZbemFryCY7Mt2ioHpdtxMqM9fuzdi3T2K27O2IjsrHIRp5SKLUiHVamu71ME6UdmjwSY/VdtPuXY0uVlz1mmoefvfduePwq31uh112cvii/V15zd+5s/HbkeU8ulu5cd9LVR1zlb3da37NIBp5CKLSiG3tVYfOrhc0lK5dGnzQnUWcxLWrcXdJsoXLs6PX4+FNawEAm8emlLsZL+eZ+1gv2znhksMNQJ2mms1YSKcaf3cpQsPOoZX4ncvdkFu1szBxZjt3kDMXZ413R0D9vZwLEZsepbO96zXyVmhdUYQVmWhCuq1e0FC6doVBJSVW24SkVY002c142djdx7ozP1ME2MEx9q/KPseDt6yuS4oan8zX2dyzmYqQbFeVwWzGwtabVzX0UPVryG3yrDuPyfZb+Pn52VrU0JmZEqw0IZuxcLZYCiQvgtaaV0XCNENXC/JWbfWjEFgmiRzuLZ2z6ps7TtZrdW+XgDVxACXFKZu0RUe3sNz7xcM1wbMoY8FKk3LL7rcI6RRG1edUJoirR/cEvKJ6fuHl8/Cjn130PS5jpbH15lWhHNN+n3HLA5VzuFRmLJzfh6kt7wx0bj8Tj5NW9AvoatNKq7b6UVQc9BMYugQEW1Ay6oskPXjLagBQOhrbVSHRxAGUFKds0qpG6n4P5y6Wa+afQrHkmeXiHEP1ewl6bifNfi8vnS/jumuW1MyFaSLcsW5ZnYPW6YC8b/wwrrnnaSwf3YNr7nka940fbur8gHkWbJjF3O1MXdxvKc1Di/utloS9drVG7heQH1Y7jKLFk9dWLKfQXFU/QobeaeXcfbSrJZWJAygpTS+S1KYLMN+au5OH3GPYBEndNxHSqu8rSPZksVTGgR+cwUO3rdE27La5b/wwHjtwovZ3mbn29wPDqw3P2IipgA67aLl3BO10pne1IPd6OJoxswSpVKe7mTpBolut/UwBXrsP2/7Zjh+V88es244nwWaetKqRIxtW4K6xqdCfd3cSMv2O/RYv5+93UcbCAitV1/h589iUsTAvMxs9c5/75knl648dOIGhq5aEvkcmi2WUi3k78xa6WpD72a2a0Q5NbpJpPQYTQeLn4PQT9HEkwyShvokXcScIORkezAUSim7cnYTMhFbKc5s/Ppmvqz1SKFacgQ87KgAG7dTjNG06Fwgi1BYIrySfZhLlVPLAShMWzusL7Nz0Io5kva4W5CY/tjDaoemNChL364efKaAZodmKH974ZB7nLsw2vC5NL/SYCHGdOcPtGDdzvnmH+G176kiDY7VUZmx76oincAQqzrc5zbi2QmN/xhkN47coFEtl3P34IWwem2r4rfopTu3YhcWVS9HVzk7gUqxnVC2zgsQfRxkZ4ZeZtn7lQMNjaRpHG3U8tT2mO1xN5eiRTNBL6H6jaaLaPb993TKj4l/u34sqJ8HP8a9L+Xe+PjyYw8Zrcw2/vXSatMtEmqipBCSd89wkuEEVxx8lceVSdLVG7iQq51aQ+OOwWrJOQ9Zp8OOTeew6mK/T1AjAxmv942t117N19xHfH7luXJ2jrX9eX4MQ76RM0FZj6jcZumqJkVbZLn/F/qPTDbuEUpmxuN/C+dJcw/VEmUXqfPaSEFIa1xx6RpBHta0KcqPCLB5hhJsuomX/0Wnf8XTXUyiWMD6Z97RF6sY1/Y6SlpQTN6a/0TC2fS+lQrcgZzOWMgmICHV5ATpzSGGmhNvXLcPnvnkSZWakibDx2lztby/syC2T2Gz7d5UEn0xcc4hEkBPRDQD+CkAawN8y8/Yoxo2aZp1b45N5ZUgdoL5RYRaPMMLNS3D6jeflFPM6p9e4pj/mJGhQSaNVDtj1KwcaGhdnrDTWrxzQLshbb16FkScONYQ82j//WmMMDdl+C7sO5mvPS5m57m8dzvR7wDucFbj0u0pCSGlcc2hakBNRGsAnAfw6gBcAfJuIdjPzfzQ7dpK4b/ywtoO3faO8TCLAJVOEylFjE7SAFuCtBfgJy/UrB+pidlXHBHnvVKGI29ctUwoN93cUZFEUwuNlett/dNo4bDVIYShCReCrxvYqP+vOoXA/P15CMgkhpXHNIQqN/E0AvsfMPwAAIvo8gPcA6BpBPj6Z1wpxu2YCAE8ThqnJxFkTwwmhkrVpGlZl/8B1ETu2sLQ7ravwEqi6xcPWwlRCA6j/jlQPs0S1RI+X6S1I2GqQNH0GtN3qy8wNtnKTtHUTIZmEkNI45hBF1EoOgDOC/4Xqa3UQ0Z1ENEFEE9PTeuGRRHbsPaYNDSszY8feY7hrbMrTW23qzdYl7tleepW3XpUePL8vhc1jUzh3Ybah5ZVTWHpp3V4CVZWOT6hENOiEhl+2YRxde3oBL2EdpFRBkJ1SLpvRHm/f5zDdmuyoE79Kkb1GFBq5KsqoQRwx86MAHgWAoaGhsHkPseAl7OySoH6fjdoe7Lab21qAW/MvFEuwUoTF/VZdRp79Oa1mnbGMtSNVEa8w17l+5UBsKc7djJfpLYhNV5dQ42we7f68buxmtFaJdmokCo38BQBXOv6+AsCpCMZNDF6aiN+KZH/WVPOxS4aaoBKOyga/1ZZXqthZXVPirTevahjbHfMNoBajb/I9+Gl0ztTsJHdB6jS8Gk/75Sc4UR27471rsOPWNXW/2wVWSnt8FDuuTqp73y6i0Mi/DeB1RHQ1gDyA3wLw2xGMmxjCFgyy0lTTTEw1H12kgAqVYAyq+bvtjna69OaxKezYe6z2sIcJN7Tx0tCcOG3mEpoYHX625SDaserY8ck8LsxeyuM8M1NSZlSGxb0zizJLu1toWpAz8ywRfRDAXlTCDz/NzEeanlmC0D0IfnUmFjoSYEwcNfYPtjTHNc9+Lqtu4kqomCLcn9WJfy9tWGeWcQrrMOGGNm4tTFccypl9GNYUFcQc00umm1Y64Fq56Kp+kzolqpejnSKJI2fmpwE8HcVYSUX3IHhpmIViqaGhgmmCje3ZdwoXZ+QMo1INbs93XsSNb3y1Z7du0wp3KmFsP5BeJYEf2bRWK5xz2UxDVMHE8dPKkMf3vfmShS5MYkUQ26m7IJQzJrpbhXmraGU+gC7iJkhjlV6gY2qtJLEmh20D1PXWBGBs3/Wz+6nSoIHKNvaxAye0QjyXzWDjtZW0edV3525AoMIruoEATBw/DVWLRStFyofrgeHVuGPdsoYmA85a0152XR1BbKdeBaGEYLSySYduMbBr8Udpe+9kOiJFP8leavv8fqnEdtU252ec+DXBCFIq1MauUe3+7jaPTeGusSnkshmcuzDrmwJt7yhUWjej4qRUmfRftqBPe38eGF7t2SQgTGKFiWZo7z5MCkIJZrQym1G3M7MbqggVOkKQx+348rOluoWOV8y5ewHys20D/qU9dSzNZrRbU9NxdVq1E12WXiGAUFR9x0EwKZ9g0vBaCE4rsxmTkHbfCRD71D1oBUNDQzwxMWF8/NWje5SCjgA8v/3GyOalQpcW7LWV89Ogde3ZosSeYzPNCoBKgSQwtELS77PZTCV+3Y6GOTNTqjly7f9mMxbOXZytM3VYKQKoPgVc9717fY/Oz5jsbLIZK3DjXaG19JJT2g8iOsjMQ+7XO0Ijj7OqWZjdgF/VNq/2bM1AdKmgkR3La9oLUoc9XlAhbn/WNlU4q+g5iyi537NRhV/qvnfd92iXT7CP93O+WSlSxs8L8ZKEtPuk0xHOzjCOr6gI45H3c4L6tWcLg5Um9DnOZ8fyOkMUdSzutzwdtklC9Z3pvsc55joB4LXwp4mw49bGxsCC0Al0hCBvVYaYCWE98sODubpwOie2cG1mR9Fvpeq+j4Xz+hq02GKpjC8detFznIyVxpabVuGh29Y0LJZJJEgNEPfrXouaW+gLQifREaYVIL7tVTPOFl1lQfv1oF3InTBQ57XXVaZTmS1sFs5L42O/Wb8g1mLGHWaapBCkBojq2LCVHgUh6XSMII+LZjzyzYYU3rFumbZWeLFU39o2jC3cSlc2ZG5n0u3rlmHs2ye1tabj4peWLVI6Om0buTMbVnWPwlZ6FISk0xFRK52KTlib1GkxOc5O399/dFqZupyx0lhgpQLHRpvOL+iYGSuFmZKut7oZzsShoBFFuvuxuN/C5P3qSBWJmBCShC5qpSNs5J2Krma3qZD0Oy5fKOKxAydqwslOXQYu+RG23BQ8CqMVSzsDWLxwftPjOCsk6iKKNj8+hcGP7GvIZNU5zXXfkVRgFDoFMa0EIKh2pjLLNBMKaIKduuy0n2976kjLMxZ1nY1sCNFE6ZSZa52SdKdzhj2qsoBN72HUiWii3QutQgS5D86CUk5tOmyZgMX9VsuFqi0wnenorTCX2FhpwqyPPZ0R3bUHXQydwjeI0zzKYlBJLjMhdD5iWvHAXVDKLar8itmrtuZn21DLY2k24zv3qFg4Lw2w2fg6IZ6xWv8zDCN8oywGJc0QhFYiGrkHJpmXdgSKapus+nxzrj4zG7tdK10198X9Fvrn9UVm4jlfmguV9elkfl+6IQonalJEDSWFVTjNH4syFqw0NZQJCBPh0o766kLvIhq5B6ZanM4JFlQLdDZPVmmpGSuN29ctQ85DI1zcX+m1qTt3YaaEZ0evxyOb1jY0ZQ5Ds0Ic8I51t7ETn8JSZvZ1WLp3UIViCeDKd9psIloY7V6crYIpIsg9CLKFVm2Tg3Ydf3b0ejy//UZsuWkVllQjPOzUeVuIPDC8uiaIvSIwdOfO9ld6K+7YeyxxceJeMIBzF2YjWXyKpTK27m6sO27S7xRAqLr4ra6vLvQ2Isg90IUP6nBrwarPWylqEEbOB9pt21Z1CgL8yxaMbFihFHo/Pz+L8cl8R/Y3dGrIUYw1Ppmva1ji1wuyGQ05TJmJVnbeEboLSQjyQWWj1LVFUxW719XZ1tk9/TI+7TA/Xfaik7Xb9inNFqlq+n3n6OONUAQlBBb3WzhfmvP1g9j3VXdvWtXkoN3nE5JPR5exjRPTXp26bbLu8zoB7Kdt2bHa+UIRd41NYevuI9h68yrleGc1tmeveO9OIQr9wyQU0r6v45N57QKbLxSNnZJBnJfSVEEwRQR5CFrZESVo0lChWNLGI7cjAambSRNh47WXWvl5MfLEoVr1SV2MeNBY8lb+zoTuQkwrCSNs1yCdWcd0rFYmDHUyGSuN+X0po8gaN+57IqYSoVnEtNIhOLWwINp03uGQc8ZBz5bNFoQwrdya4eYj+/Ghr30WS1/6CU5ddjk+/tb3Y/eq9W07vynFUjl0Fye3mSypzkuJVe98RJAnEPshClIjJU3UoIEH0SLbLcS3/9Mn0D97AQBwxUvT2P5PnwCARArzsLhDQONsWahDSgd0BxJ+mEDshytIXZIyc+Q9QFvFh7722ZoQt+mfvYAPfe2zMc3Im4Xz0p7dk9IpqjSLdpCx0li/cqAu5nz9yoHYWhbqkFj17kAEecw445jtBJMwAjmXzQTeoruFT7tY+tJPAr0eNxdn5/BLyxZp33/5/D7suHVNXXw7gTH2rZN1Mee7Duax8dpcLC0LdSTV3CMEQ0wrMaLb1gYV4rZWF8SunstmMHNxtuWVGFWcuuxyXPFSY9u1U5dd3va5mFCaY/zbD05r37fDPM876sWoGmgUS2XsPzptlGvQLuGeRHOPEBzRyFuMSuO20W1rg5Amqml1qkxSFY9sWotnR69HIQYhDgAff+v7MdNX32Ripm8+Pv7W94Pi2ST44hV7vzSbMd5FuTXd+8YPY/PYVJ3mftfYFNZu29eWmiphSgcIyUMEeQvxS+kOsn3VpfY/dNuamvbmTAPXcce6ZbXj49K6dq9aj9EbPogXLhvAHAgvXDaA0Rs+iC+/8frENXw2YWTDCuN76fzOxyfz2HnghDLs084PaLUwD1M6QEgeEkfeQvx6RJo0YAYupeMDwZJD/Lbs45N5jDx5KBHFsxbOS+Pi7FwtqSZJEIAFVkpZarffSuE/Pvou43u5uN9CYaaEpdkMzl2Y9Y0skhhzwYnEkceATks7M1Mp2KRKwfYiSHcb9/G2UN88NtUg1NvRCs5NmghzzL71a5LAAiuFjddegbFvnaxbaKwU4c9ueSMAdTq9CmcLOhPE6SiYIIK8hXilyO/Ye6ymafkJsSCxvboiXV6xwsODOWONMiqcJiEA2Dw21bZz60ihkt3a2AlqDrsO5rHpTVdi/9Fp5Q4nbCKXH+J0FEwQ00oLGZ/M4y6NgCIAz2+/se5YP43Ob5utGiNjpbHASik1bud4XnM1wU7xTxtkiNqmJee87378UFuTksKQJmpYgFSELbPgJmOlO8peLRmiraclphUiuhXAVgC/COBNzNz90jkAw4M5bN19RGkHVWla8/tSng+/3zZbFwWjG/OUo2pfs1rkw5vWYngwZyTEXv/ql+O67c/gVKGIbL+Fn5+fTbwQBypJV7piWG4B9uAtq+te09nD7dZ7dkkFItRs6J0kCCVDNF6aNa18F8AtAP5vBHPpSrbevMq3FKmpBue3zQ5qT12UsSLRHG1s4bwoY4HAylhqAHj2+5disuOIY28GO+vR6XtwC7DNY1O4fd2yut2TyrFspQlbblKXIO40vDJEu+H6kk5TgpyZnwMASmrwbwLQlSIFLgk+k4JVJrG9QcrWEirNGaIS4u4aLybx7J2K3XBbd+8YwM4DJzB01RIAHnbz5G9CjJEM0XgRZ2cbcEebuLU4LyFOgHHn93yhaFyOlgGjhCCT8QiNC0KxVE5sadxsxsLF2bJ2x+AH4VLUie7eMYCtu4/gwqy+A1FpjrF19xHf++oVPpoUm7RkiMaLb0IQEX2ViL6r+PeeICciojuJaIKIJqanG9OzewnTLMBcNlNr+uv1sDt7fDIu9RXNZTNYOE+tGeeyGd+HLJfN4OFNa2vJIjp0wjqJQtxKARdm50ILccD8ugrFku99tnuHuvFLJmumf2grkAzRePEV5Mz8DmZ+g+LfPwY5ETM/ysxDzDw0MDAQfsZdgMl20/QhUC0KjEtJRBdnGwWWlSaMbFjhO76t4T07ej2e336jtulxTLW3QlHm6MxJUaGqNOhXlTBpVQslQzRexLQSA7ptqDtJxuQh8LJN7th7TJkpuXBeX21sXVTN4n6r4fw6C9D8vhRmyxw6K7PfSjWlIdvYjam9aFfiqFfYpxvVPfSzOSfRJh00YU2IjqZqrRDRbxLRCwB+BcAeItobzbS6G9029KHb1viaUtzozCNLPcraOpsyb715lXIuW25a5fk5J+dLc3jZgnA6weJ+C392yxuVc+i3zH+eC+el8dtvXhZqDs1ipQl3rFtWp43avT5NyCp2Ol731eR9obdoSpAz8xeZ+Qpmns/Mv8DMG6KaWDcT5TbUyza5KKM2hThfDzIXL+ERtpJiYaaE4cEcNl6bQ7oa/WQ3PVbVNtFx7mIZuw7GYx+enWMMXbWkZoIa2bACuw7mjUMrVTsdP5uz1/v3jR/GNfc8jeWje3DNPU/jvnHvxtFC5yOmlZiIahvq1Wl921NHlJ9xRosGiXxQ1RMJUwvdyaKMhfHJPHYdzNciQMrM2HUwj0UZy7hdXZShlEFhRl3yS9DGILqdTqVQV2WcbMbC1psvxZzr7vvE8dN47MCJ2hhl5trfDwyvDnxtSYqMEfSIIE8QYR8a3aKg05Lt102z8ZzzyvZbmN+XwtliY/ahSshvvDaHPd95Uaudnrs4i627jygddwusFKwUGdne404MdSa/BLVT2zsdrzDSCwqnteq+3/34IeU5PvfNk4EFuWRrdg4iyBNCKx4av9hek2w897zOzFSSfeyUfBuvnYEtQAY/sq9BoJfKrNW6z8yUEhMRYxITbwvwbL8VKGN1ZMOKhu+5sXCXWZakLq7dK1dBp0BItmbnII0lEkIrwsn87Kw6zTFfrcESdF52qOLDm9YCqFQ0dHZFCmpHr0TxBPpIA3esM3eApgBt3L3JNOwFMsjuwF6n7n78kK85xkTTT2uyrHWve8Wjm/w+hGQgGnlCaEU4mZeWDHin9Nu7gaDz8tpZeJ3PSlNdHZKMlW7a5r2438LQVUvqbMZuMlYK50tzdd/N8tE9oc63/BUVQa6zeatgVL5rk6JhOmezU6POWCnMlBrHet+br1R+1muhNvl9iGaeDESQJ4RWpTh7OVW9miH4PcxLs5k6m65dvlZVxtYea2TDCmwem1Jqtwvn9WHh/L66BafZqoxbblplsKOhBjNRLkDNGiff+P5pjE/mA9W8SRMZLVi6BDH3wjlTmkOKqnXVuTL++958pdY+7rVQP7xpre/vQwR5MhDTSkKII8XZDj3UcapQ1M5r/cqButIAzogT3VjDgzmtieJssVQL37Pj6Ec2rIAV0kjeb6WMHI8qM5FpE2s3DNQWLBOslH+xNJv5fSlMHD/d0MhbpVHPMbB0UQY/3H4jvv/gb3g6Ob1CSk1+H0IyEEGeEOJKcR4ezGmbNTsfZve89h+dDmT6sAWG17mAioZpC6sde49h05uuRNYR975wXtpIuNsx6CY7GrdAsq9ZZ1f2G8v4nhHqrs2LQrGExw6caLBl6zR/UyHrp0D4/T6EZCCmlQQRV4qzzsRy7sIsxifzynkFac3mTmTRxaKr7Ou7DuYbFjSTZhi2kDHppakSSLqQSr/olSDCrVRmXJwth/YHFEtlbUcm03n4+VEA73smJAMR5IK2CXOhWNI6tfzswLq6MV6C47rtzxiFuzkXlvvGD2PngRN1wtVKE85dmMXy0T0Ntnu3IPYSSKq5el2zc6zFhiGIM6U53LFuGfYfnQ5lly8zNywEQYWsnwJhIuyFeJGenUINXQNmVa9Qr65GYXtNXj26R6vtetVldycs/fz8rDKJyE5Q0jVQNkH3Hbn7eao6AumwK1WG6VuacziGO6FuudAcLenZKXQXQUIN3V3jbY0314Sg8EqkcdqFnedXzUVHsVTG/qPTng2s/dCZGdwLl1uL9So3YLeHC6pS2Zq3l0Yt2Zm9gQhyoUbQEMiobfomymixVMbW3UfqBGSpPIdzF81szM1GWgQxM7i/H1Vmq42JEL/umiX44U+LgTRryc7sDUSQCzXWrxxQJs+sX9meRiCmiTSFYqmm3ZoW1bLRVYQMgmoBMzFfbLmpsRF3EH7402Lg3UQS65YL0SOCXKix/6i6BZ/u9agJkkgTlqARhSYC2tR8YWoC0hFG+Eovzd5A4siFGnFrb2ETcYIQpN6LaV9ML/OFMy7+uu3PAACeHb1eG5vtRRjhK700ewMR5EKNuLvOqJKP7li3rK5XqKlGrTssyLWYFgzzKi6lWwhUAtaeczZjwUrXX0FY4Su9NHsDMa0INZKQ+OG2P9tNJ2xMHKLZjIV3r3k1dh3Mh7oWv4Qjt+D26sGqWwhsW7fObBNlyKD00ux+RJALNZKY+KHrtmMnHC3KWCCqmEzc8x26aknga/GKj7dxa/W6BVA3hr0QeAlYEb5CEESQC3UkTYDozBZzzHh++42en3Vei63hbh6b8hTqfm3aVFq9bgHUafXiaBSiRgS5kGiiiLoIkhTj5dj1SnZSLRqqlm3iaBRagTg7hbbjjuTw6jYTRdRFkC5HugUiTYRThWItEkWHM9IFqAhx220pjkahVYggF9qKaUifTRRRF0HCKnUhkGVmo/mqFg3GpXo1IsSFViCmFaGthEkZb9Zu72eecUeIOAtrpTw6HqnmFHcsvtCbiEYutJU4BJ2XeUa1Q9h1sBLn/fz2GzHn0fFIRdyx+EJvIoJcaCtxCDov84yf/TzofDsxkzKIz0JIJmJaEdpKXElHOvOM3w4h6HyTGIvvhZS57Q5EkAttJWmCzs9+Hma+SYvF90LK3HYHIsiFtpMkQWeicSdpvlEjztnuQGzkQk/T60WlxDnbHYhGLvQ83axx+5GEQmlC84ggF4QeJmk+CyEcIsgFocfp5R1JtyCCXBBaTBS1xaOsTy50H00JciLaAeAmABcBfB/A7zFzIYJ5CUJXEEWctsR6C340G7XyFQBvYOY3AvhPAPc0PyVB6B6CVF5s5RhCd9OUIGfmfcw8W/3zAIArmp+SIHQPUcRpS6y34EeUceS/D+DLujeJ6E4imiCiienp6QhPKwjJJYo4bYn1FvzwFeRE9FUi+q7i33scx9wLYBbATt04zPwoMw8x89DAwEA0sxeEhBNFEa1OLMQltBdfZyczv8PrfSL6HQDvBvB2ZpMe54LQO0QRpy2x3oIf1IzsJaIbAPwlgF9jZmN7ydDQEE9MTIQ+ryAIQi9CRAeZecj9erM28k8AeDmArxDRFBH9TZPjCYIgCAFpKo6cmV8b1UQEQRCEcEhmpyB0CJLdKegQQS4IHYBkdwpeSD1yQegAJLtT8EIEuSB0AJLdKXghglwQOgDJ7hS8EEEuCB2AZHcKXoizUxA6AMnuFLwQQS4IHYJ08hF0iGlFEAShwxFBLgiC0OGIIBcEQehwRJALgiB0OCLIBUEQOpym6pGHPinRNIBzAH7S9pMnh8sh19+r19/L1w7I9Tdz/Vcxc0OLtVgEOQAQ0YSqQHqvINffu9ffy9cOyPW34vrFtCIIgtDhiCAXBEHocOIU5I/GeO4kINffu/TytQNy/ZFff2w2ckEQBCEaxLQiCILQ4YggFwRB6HBiFeRE9FEi+g4RTRHRPiJaGud82g0R7SCio9Xv4ItElI17Tu2CiG4loiNENEdEPROKRkQ3ENExIvoeEY3GPZ92QkSfJqIfE9F3455LuyGiK4loPxE9V/3d/3GU48etke9g5jcy81oAXwJwf8zzaTdfAfAGZn4jgP8EcE/M82kn3wVwC4CvxT2RdkFEaQCfBPAuAK8H8D4ien28s2ornwFwQ9yTiIlZAHcz8y8CWAfgf0V572MV5Mz8kuPPhQB6yvPKzPuYebb65wEAV8Q5n3bCzM8xc691Dn4TgO8x8w+Y+SKAzwN4T8xzahvM/DUAp+OeRxww84vM/O/V//8ZgOcARFZcPvbGEkT0MQDvB3AWwPqYpxMnvw9gLO5JCC0lB+Ck4+8XALw5prkIMUFEywEMAvhmVGO2XJAT0VcBvErx1r3M/I/MfC+Ae4noHgAfBLCl1XNqJ37XXz3mXlS2XjvbObdWY3LtPQYpXuupXWivQ0QvA7ALwF0ui0RTtFyQM/M7DA/9BwB70GWC3O/6ieh3ALwbwNu5y4L6A9z7XuEFAFc6/r4CwKmY5iK0GSKyUBHiO5n5C1GOHXfUyuscf94M4Ghcc4kDIroBwJ8CuJmZZ+Kej9Byvg3gdUR0NRHNA/BbAHbHPCehDRARAfgUgOeY+S8jHz9OJZCIdgFYAWAOwHEAH2DmfGwTajNE9D0A8wH8tPrSAWb+QIxTahtE9JsA/g+AAQAFAFPMvCHWSbUBIvoNAI8ASAP4NDN/LN4ZtQ8i+hyAt6FSxvVHALYw86dinVSbIKK3APg6gMOoyDsA+DAzPx3J+F22mxcEQeg54o4jFwRBEJpEBLkgCEKHI4JcEAShwxFBLgiC0OGIIBcEQehwRJALgiB0OCLIBUEQOpz/D2iGxovoWvUmAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot dendrogram\n", "# Read https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.cluster.hierarchy.dendrogram.html\n", "# to learn other options\n", "\n", "plt.title('Kmeans Clustering')\n", "\n", "plt.scatter(dataTW[:,0],dataTW[:,1])\n", "plt.scatter(codebook[:,0],codebook[:,1],c='r')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[[ 0.92921799 0.7727284 ]\n", " [-0.96452347 -0.72228906]\n", " [-0.92008042 1.51778476]]\n", "[1 1 1 ... 2 2 2]\n" ] } ], "source": [ "# Compute k-means\n", "# to learn other options: \n", "# https://docs.scipy.org/doc/scipy/reference/cluster.vq.html#\n", "\n", "centroid, labels = kmeans2(dataTW, 3)\n", "\n", "print(centroid)\n", "\n", "print(labels)\n" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEICAYAAABCnX+uAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAACWQUlEQVR4nOyddXgU1/eH3zvr8eDuUqS4FqelUKQUWqrUqX3r/qu7u0FL3aFQakhxd3f3YAmB+OrM/f0xm5BlZ2MEAnTe58kDGbn3zuzmzJ1zz/kcIaXExMTExOTcRSnrAZiYmJiYnBqmITcxMTE5xzENuYmJick5jmnITUxMTM5xTENuYmJico5jGnITExOTcxzTkJuYFIAQoqcQIqmsxwEghJgihLi5rMdhcvZhGnKTYiGE2COEuCTf79cKIY4LIXqU5bhOBSFEByHEZCFEmhDimBBimRDi1lLuI+S+lQQp5WVSyu9La0wm5w+mITcpMcHZ4WfAACnl3LIeT0kQQnQGZgFzgQZAeeAe4LKyHFd+hI75t2oSEfPLYVIihBB3Au8BfaWUi4Lb6gghpBDiViHE/uBM/W4hRHshxLrgjPfTk9q5TQixOXjsVCFE7Xz7Pgq2kyGEWCmE6JZv34tCiN+EED8IITKFEBuFEO3y7X9SCHEguG+rEOLiCJfyDvC9lPItKeVRqbNSSnl1hOuWQogG+X7/TgjxavD/FYQQE/PN7OcLIRQhxI9ALeAfIUSWEOKJ4PGdhBCLgsevFUL0zNfuHCHEa0KIhUAOUC+4bURw/y1CiAVCiHeD9263EOKyfOfXFULMC17/DCHEZ0KInwr6TE3OYaSU5o/5U+QfYA/wO3AEaHnSvjqABD4HnMClgAf4E6gEVAeSgR7B468AdgBNACvwLLAoX3vD0WfIVuBR4DDgDO57Mdh2f8ACvAEsCe5rDOwHquUbV32Da4kCVKBXAdfbE0jK97sEGuT7/Tvg1eD/3wheuy340w0Q+e7bJfnOqw6kBsevAH2Cv1cM7p8D7AOaBa/fFtw2Irj/FsAP3BG8/nuAg/n6Wwy8C9iBrkAG8FNZf3/Mn9PzY87ITUpCH2AJsD7C/leklB4p5TQgG/hVSpkspTwAzAdaB4+7C3hDSrlZShkAXgda5c7KpZQ/SSlTpZQBKeV7gAPdSOeyQEo5WUqpAj8CLYPb1eCxTYUQNinlHinlToNxJqIb0UMluw1h+IGqQG0ppV9KOV9KGUnMaDgwOTh+TUo5HViBbthz+U5KuTF4/X6DNvZKKb8MXv/3wb4rCyFqAe2B56WUPinlAuDvUrpGk7MQ05CblIS7gUbAV0IIYbD/SL7/uw1+jwn+vzbwUdC1kAYcAwT6bBUhxKNBt0t6cH88UCFfW4fz/T8HcAohrFLKHcBD6LP2ZCHEGCFENYNxHgc0dANYGryD/oYxTQixSwjxfwUcWxsYlnvtwevretJY9hfSX971Sylzgv+NAaoBx/JtK0pbJucwpiE3KQnJwMXoroORp9DOfuAuKWVCvh+XlHJR0B/+JHA1kCilTADS0Q19oUgpf5FSdkU3mBJ4y+CYHHQXxJXFGHMOuksmlyr52suUUj4qpawHDAIeyeebP3lmvh/48aRrj5ZSvpl/iMUYV34OAeWEEPnHWbOEbZmcA5iG3KRESCkPAr2BfkKID0rYzOfAU0KIZgBCiHghxLDgvlggAKQAViHE80BcURoVQjQWQvQWQjjQ/ehudHeLEU8AtwghHhdClA+e31IIMSbC8WuA64UQFiFEPyAv7FIIMVAI0SD4lpIR7DO33yNAvXzt/AQMEkL0DbblFHrMeo2iXGNBSCn3ortpXhRC2IOROYNOtV2TsxfTkJuUGCnlfnRjfpUQ4o0SnP8H+kx5jBAiA9jAibC/qcAUYBuwF90gF9U94ADeBI6iux8qAU9HGMOi4DX0BnYJIY4Bo4HJEdp+EN0opgE3oC/k5tIQmAFkoc/0R0op5wT3vQE8G3SjPBa8d4OD40oJXtvjlN7f5A1AZ/QF1FeBsYC3lNo2OcvIXeE2MTE5jxFCjAW2SClfKOuxmJQ+5ozcxOQ8JBi7Xz8Yx94Pffb/ZxkPy+Q0YS3rAZiYmJwWqgAT0OPwk4B7pJSry3ZIJqcL07ViYmJico5julZMTExMznHKxLVSoUIFWadOnbLo2sTExOScZeXKlUellBVP3l4mhrxOnTqsWLGiLLo2MTExOWcRQuw12m66VkxMTEzOcUxDbmJiYnKOYxpyExMTk3Mc05CbmJiYnOOYhtzklJDqQaRvOVI7VtZDMTH5z2JmdpqUCKnlINMeAt9iEHaQXqRrGCLuOczykiYmZxbzL86kRMjMl3QjjhdkJuAD9wRkjlkW0sTkTGMacpNiI6UP3JMIV0V1Q/a3ZTEkE5P/NKYhNyk+0oNeIc1oX8YZHYqJiYlpyE1KgogFi1GZSwH2jqXalZQa0rsU6Z6EVEurRrKJyfmFudhpUmyEEBD3MvL4PYAPfXZuBeFCxD5eav3IwH7ksZtApuVuQEZdjYh9FuOazyYm/03MGblJiRCOLojy48A5CKwtIOo6RIWJCGvdUutDpt0D2iGQ2foPXsgZD55IVdhMTP6bmDNykxIjbI0RCe+clrZlYA8E9hHui3cjc35CuAacln5NTM5FTENucnYic0BYwKjuicwyPkXLQma+B56JgAbOvojYxxFK4mkdqolJWWMacpOzE2sjjL+eDnBeFrZVSg157HoI7AT8+kb3n0jfcqgwCSHsp3O0JiZliukjNzkrkNKL9M5HehcgpRchrIj4twAnJwy6Cyw1EFE3hTfgWwjqPvKMOAABUFPAO/O0j9/EpCwxZ+QmZY70ztPT/fOT8BHC2Rsq/InMGQPqQYSjB7gGIYQzvJHAVpA+g9ZzkP7NCINZvInJ+YJpyE3KFKmmIo/fB3hCtx+/FyrNQVjrIeKeNj5XepBZo8H9O2ipQMDgqCiEtU5pD9vE5KzCdK2YlC0RQwk1ZMa7yMB+w71SSuSxmyH7Sz1EEaPZuADhNPSpF4ZUDyN9K5FaWrHPNTE505gz8vMc6d8CvhVgKQ+O3gjhOIW2toF/nZ7Vae+EEJainadlg3eaPmu2tQVbqxMJPTKLUL92Lj7w/I30/IN09kPEvxnSn/RMBv96jGfhuVgQ5X9DCFfB4wvsA98iEDFIe2fIeAa8C4Oqjj5k1HWI2KfMJCSTsxbTkJ+nSKkh0x8DzwxAgrACNij3I8LWuJhtBZBpD4J3PvosVwElEcr9jDBM1c93rn9TMDszgD5rtoGjIySMRAgrOLpA1ihANTg7OMv2TEPaLgRrHT28UN2rhycWAWGtVcB1SWTm25DzU/C6LCC9+v/xB/8P5IxFWuogoq8vUp8mJmca07VyvuL+Ezwz0X3P3mB2ZBoy7X9IaRScHRmZ82PQiHsAt96WegiZ9nDB50mp+79lBpCDPnt2g3cJMuc3AISthe76EFEFXQxkj9bbCmwushFHxBe837cI3L+gqzh6gtmjAcLfENyQ83XR+jQxKQPMGfl5inSPAdzhO9SjoO4Ea4OitaNlBmesnpP2qOBfj9SOIZRyxierO4P+65PxQNZINO+/4FsN2HSXi3YA1N0YZgFpKcbbCxz8Cb+5FkiC7K/AtxSUyuC6HDzTQBrcIyO09OL1bWJyBjEN+fmKNPI7o7tFDMP0TjpdPYJMfxx8KzH2YQPobUn/dr3IhBIPjksQSrTu2sl4G2OXCSCTwZcc/MUL/qDbJqKxLqYRB1Ci0AL74NiNoQ8UdSf4FwX7Kwqlr+poYlKamIb8fMU1GDJ3EjaTFk6wXlDgqXlZkupBIhpiAEslZOYn4Pkb3SWh66JIew/d8PkWFXPQJTDWEXGA40o42h/jiJbC+lPQr8cGwoGIfawUx2ZiUrqcsiEXQtQEfgCqoH/zR0spPzrVdk1ODRF1HdIzJZgokwM4QCiI+A8gsB3pWwQiDunopWdFuifr7gvpAfUAkF1A61Y9osM1DLI/J6xSkG8O+OYRsfjEaccJ9rZgcRDZiBeEFezdQTsCttaImBEIS7XSHqSJSalRGjPyAPColHKVECIWWCmEmC6l3FQKbZuUECEcUO4X8M5B+hbrfmHn5cjsT8D9N7qRFcBTJevA1lJ3u0RceDydRtxCgW8KSnlEzEPI7C9L1rxwIWJuQ9g7lOx8E5MzzClHrUgpD0kpVwX/nwlsBqqfarsmp44QFoTzYpS4Z1Fi7kAENoDnH3R3i4/wmptFJaAvUqpJpTfYIhMNFReDSIh8iHYAeWw4RFqELQzpK/JisInJ2UCphh8KIeoArYGlBvvuFEKsEEKsSElJKc1uTYqIdP9e9CiNQvHoi4YREUAUUPIEJGN0QS1iCg59BG+wQHRxcenrC/7tyEBB12dicvZQaoZcCBED/A48JGV4BV4p5WgpZTspZbuKFSuWVrcmxUEW4I4odRyQOBKirgVRoRTb1ZN1RNQ1QGE645nFbDcabC10+du0e5BHr0A7egVSTS3xaE1MzgSlYsiFEDZ0I/6zlHJCabRpUvoI1+BCEm9OBRe679oJIhHiXoDM18A9/kTNzdLAUhuhJCKEAkpscU4ECjpeAjngXwr4gtIBXghsQh69FKkdP5VRm5icVk7ZkAtdgOJrYLOU8v1TH5JJSZFaBtL9FzLnd6Rq4L5y9gN7FyhEe6RkaMEfP9g6gG8JBPbky5YsDezg7I9UU3TdF604PnqVwmfoEcIRZSYy9eZi9FV0pNSQWhpSltY9MvkvUhpRK12AG4H1Qog1wW1PSynNCrlnEM09HdIfBRQQEjI0ZOyTCEtFZPbXoB3TQ+pin0do+5DZP4B3NiULzzMid+FUBd8s/d9Sj1zxQfYoZPZnYOt5GtovAHUnmnc5eKZCYC1Y6iCibkDYW5W4SS1nLGS+r8/+hR0ZfRsi+l79bcMEKSWeHC8Olx1FMe9JQYji6m6UBu3atZMrVqw44/2er0gtDZncnfA0eiu6S8F74nclHlF+ou4qSB1I2cV6n2sYhTzaIOYRlJjbi92adE9Gpj9FqIyCC2LuRIm59xTGeXo5tPsIP708nnVzN1G+WiLXPTWEjgPalno/C/5YyqiHviP10DHsTjtXPHAZN790DRZL0RQ3z1eEECullO3CtpuG/NxH5vyOzHyliGJSdnD2Be8MPfnnlLMpczMgTyY3/f3Mf79CiTS+0sKOqDgXYSlfrLO0lH6g7grfIaIRlVYUWSL4THJo9xHuafME7iwPmqrfU2e0gxFv3sDge0uvAtOa2Rt4dtAbeHNOvC06ouwMvLMPd79/S6n1cy4SyZCb7yvnBT6QRTVWPvBMCoYhloaR1dBnqycjOTNfr4L6EKCUL8Y4BCXyNvqWFP8c7bDxduktxRDR0uWX134PMeIAnmwvXz/9Cz5vJD2e4vPDS+NCjDiAN8fHxC+m48kpae7D+Y1pyM8HHD0onlEu7RlqpLDGMzEbL+haZFA1sajXK8BSA2xtijcEJbp4xwNYG0ZoKwFECdo7A6ybuynEiOfn0M4ID6YScHCHkWImCEVw/EhaqfVzPmEa8vMAYakGrqvKehgGaJzRr5hIAGyn0IAG6mFE3Au6pEGR+nSA/aJi9yRinwBOLiLthJjHz9pKRBVqGLuPAj6VhEqFaL8Xg/ot6xhuF0JQvloJs3XPc0xDfh4g1cPB1PuzkTO4mCrTwNGXUzPmXqR6BERMEY61Icp9ixD20GEUwc0l7O0R5b7WZ/8iFqyNEQnvoUQNKeG4S4amaSRtO0hKUuFJT9f+3xAcUaGZunanjY792xBfIa7UxnTzy9fgiAq9p44oB9c/MxS741Q+2/MXc7HzPEBLHQ7+5ZT9wuLZQBR6NaJTaeJmyPmBwu9nNKLSPIQSi5QqMnsUZH8LMhMs9RFxzyEcxZ+tnylWz1rPWzd9QnZ6DpqqUad5TZ4f9xiVa0fOvJ705XS+fOInNFUj4FfpNKgtj397L67ok98uTo1NS7Yx+vEf2bF6N+WqxHPdU0Ppd1vvs/Zt5UxhRq2cp0gtA5ncmcjFH0yKj42i3U8rIu55RNS1aBmvQc5YQkNAnYhyPyLsLcPOlP514F2s+8SdlyGU0pvRFoXDe5K5o/kjIYuHikWhYs3y/LDj0wLjtv0+P4d3JxNfMY64csXJrjU5VSIZcrOwxDmPmRFY+hT1oRhAZn2J9C4G77+Ez+A9yKxPEOW+AkD6NyP963UZYf86vR9hh8w3IPFLhL19KV5DwUz+cgaBQOh3R1M1MlIzWTtnI617XxjxXJvdRs3GpsDp2YRpyM9BpG8VMvsrXUbW3hEsNY1jkk1OP9p+8O6PvD+wEyl9yOP36vVC0QjJpg2GGsrj90Glhbqy4xng8J4UAr7waCOpSY4eOHZGxnCmkFIS8Aew2qznrWvGXOw8B5DSR64LTHNPQh67BbwzIbAFcn4GLRndN5zrp3ShC0TZgWiKXpvyv4yV03KfbI11iQTfUk7owBvhC87SzwxtLr4Quyt84dCT46VWkxpnbBynmxk/z+O6mncxwHUDV1W6nQkfTaQs3MmnG9OQn8Vo7n/Qkrsjj1yITO6IlvUVZLyEbhByv4wBfVbn7K1rdLuuR8S/hKi8CFFxJiL+VYh9ugyv4hzBUp3wcMBTRSBi7oec3wiXTzDizBmYXtd1McwhE0Iw44e5Z2wcp5P5E5by4V1fkHrwOFJKMlIz+faZMUz4qCQ69Wc3piE/S5Ge6ZD+TDADUOqhdVkfBdUET0YF30qUmFtR4l9EuK5ACAfCUhnhGoASfbMuLXu2MiED0X43otp2RPvdMCFMzv70o+4lVPckPyVMl3f0RdiaUTRhMqtePu8M4fcGkFq4JZeaZM5vxS2afXby7bO/hmWIenK8/Pzq72Gz8h1rdjPr1wXsXLvnDI6w9DB95GcQqWWDugeUSghLeIiX9M5DZn2i+741N+GzuALSkwsra+bsD+7fOOuiWyZkIB5LRriDf1hJAXgsWZ+bDj2zkRyRKW5BDgVEDCIu+CbkuBTcYzFemA4WxU74xNA/LqVk/u9L+POTKWSn59Dtyo4MeXAA0XGnpisvlMhuJMVS9Pnd3s1JpOxPpX7L2iRWTjilMRUVKWWRfN3Je40rkWWn5+B1+3BGOXBne3h24BtsXb4TxSLQVI3G7Rvw6sSnSj2k8nRiGvIzgJQSmT0Ssr4AYQXpQzp6IBLeRQS1wTX3P/oMvNBX8Fw9kPwG2YWIvqPgs2IfQHpnnZjhnyWIN1JPGPHcbW4Jb6QizxpDnosFlFpgawSB/aAa1Re3grM/IuZBhKUKACL2fqRvLqip6LN+G7pRF/q/tvYRU/a/fPJH/hk1DU+2/hBP2naQmT/PZ+TKt0/J0ETHRdGsywWsn785JO3e7rTR95aehZ6feTyL5wa9yY7Vu7Harfg8fgbceQn/+/DWIhnZjYu2MvXb2ezbnERUnIuOA9vS58YeRMWe0MqXUrJ56XY2LtxKuSoJ1LygOl889j3r52/G7rBxyU09uOvdm3BFO9m7aT/bVuwivlIcDdvUI6FiHDUaVYs4w3552Htc939DmP3rAjYv2Y4/n1bM5iXb+fKJn3jgsxEh53jdXtKPZlKuSgJW29llOs048jOAdE9EZjxzkhiSA5x9UBLe1w19ShfQjhahtQSwNwffchA2kH6Ivg2URPAuAEs1XSfb1jh0DFJDpvQCzVjHoqwQ1bYjDL6CUoA8GEGPpKwQ0YhyYxG2RmhH2uqJP2FYEJVXI0SokZXSDe5JSN888Mwk9EFsBWtdRPmJIUYweV8KNzW8H9Uf+kbgiHJw1zs3Muievqd0Ocn7Unio23NkpWXj9/iRQExCNPd+fBvdr+pUoGTsc5e/yYppawn4TrxlOKMc3PPhrfQfcXGB/X724DdM/nIGPs+Je2CxWoivGMfIFW9RvmoiakDlhSFvs3bORj3ixGrVY94FefMQm8NG086NiEmMZvmU1QT8at5DqVKtCgz6X19+emkcXrexa8sRZUcNqIbRO85oB/9k/gSAGlAZ/fgPTBo9A4TAarNw88vXMOT+/gVe5+nAVD8sQ2T2FwaKdl7wTENqWSAzQEsvQktOiHsUpdw3iIrTEInfQoV/9fT8zPfBNw/cvyFTh6Clv4GU+b7A/nV6P2cb1SPMbCJtL0ukHyy59UcjjU9g9GclhAsRdRVY6hEeHRMA9QD4V+dt2b/1AHe2fCzMiAN4c7wsm7I6bHtR2bVuLy8MeZsHuz5LYqV4mnRsiESPI09PyeD9Oz7nmQFvoAaMXUpZadmsnB5qxEH3P0/4cGKBfW9buZMpX88MMeKgG8u05DS++j/deE4cPZ3VszbgyfYS8KknEpfyPfT9Xj8bFm5h6aRV+Dz+kDeL5H1H+fnl8dz9/i3UvbCWoSvJm+MzNOJAyPi+eeYXJn05A6/bhzfHS3Z6Dt889Quzxyws8FrPJKYhPxNokXQsFH1WJ6KIbBjs+o+lJsS9ghJ1DQDCUlXPGHSPAzWZEy4ZDQiA+1tkyqVINTjLl5lEDq8ru/BE+VR5pCu0f+kSyKeKp+99+nGAoycidy3CNUTfFoIVHN3DtFdCCOzCePFT6MY8yItD3yE73VhqQLEoVKpVsoLW21bu5IGLnmHR38s5mnSMbSt3sWrGelS/itR0K+nJ8rBx4RYWTFhq2EZOphsRIfMzK81oMf4Ei/9eEWbEc9FUyZJ/VqIGVH588Td8EWbS+VH9atgDJRef18/2VbsYvfY9ylUxXuw3MvBCCFr1ag5AwB/gr8+mGi6a/vTKuELHd6YwDfmZwN4Rw1utRINSGSFsEDWc8PA3FyL+TZQqG1AqzkRYqqClXoN2pD1a6jC0nD8h+xsiRkVoR5CZr+r/t7XSZ5RhOMFSv4QXVgoMjUO+WwlZw6q7U2pYke9WOosWOkE30L0R8W/nbRGxD4KtafAh7NSlZy01EHGvFtyUvR16nP9JyIDeHnBgxyGORFioA7A5rCV2q4x+4ke8Od5Cl0k82V7mjjOOXqlYozyxieGiYkIRJFSK54+PJ/FQt2cZWuFWrql2B28M/4h9W/SHlM1pK3Ch1ea08cVjP5CRmlX0i4qApmrs36r3m1jZ+PtksVqIjo/C7tRj6u1OG9HxUdz3yW2A/tCK9GZy9MDZU5D7LHx/Pf8QMQ8hvfOCFXxyvxROiH0urz6jiH0ESQByfg2eZIOYBxGugQBI7wLk8f+RN/P2r9V/CkQFz3S9OSUWGfskZL6NHv0iAZeuvy2iOP2VdApgaFwZLWxaIWqEnjIvD0Y+rMJ0FGtoSroQLig3RneHBLYisYN3ATJ1CNJSCRF9D8J5SVhTwjUUmf0laH5ORLE49Zm8VX+g+jz+iDNeRRE8+f391GlWsyQXzPaVRcsAFkKELDyevO+RL+/m5WHv4vcG8lwaUpPsXLOHkQ99F3L8rF8WsGDCMl6b9BSZqZloAePvmVAEDVrVYeLo6aWStGO1WSlfNZGjB49x7ZNDeOe2kfpDLIjNYaPLFe259+PbmPzlDLav2k3DNnXpf8clJFTUZXljEqKJTYzm+JFw12ejtnVPeYylhbnYeYaQgSTdV+5bAZaaiJg7EfawNQuk9ICWBkp5faYeRDs6AALbS9S3SPhCj7TwLkCqKaDuBC1DD4vL+U4PiSx2iN35ggMSv4Xjw4n0IBOVNxdYek0G9iJThwQf1LltuCD2UZTom8KPV1OQWR/p2bk4IepaRPRteZ+3pmlcW/2usCIKVoeVG58bxvVPDy3+ZQa5pdH9HNhReBEIR5SDt6Y9R7OLGkc8ZveGffz2zl/M+GlekQKhEirFkZ3hxh/BtaL3a9cXJ0vJLEXFufB7A/S5qQeVa1Xglzf+wGJV8HsDtO/Xiv/76YFCo39m/jKfD+78IuQh4Ihy8N7sF2ncvkHpDLSImKJZZYyw1kDEv1L4ccIJwbC1EAIl11KRaf9Dn3FbQSiA1I27Eo3MPMR/14gDSAhsQ0/6iVB71LcYHF0jt5D16UlGHMANWR8go64N85kLS0U94xZjN4yiKDz9y4M8d/mbqAEVvzeAM8ZJjYZVGfrQgAKvJvN4FssmryYjNYMmnRrRuH2DkEiY654eyod3fUHAYBEV9GgNVdW44dkrCzTiADUaVcXvDRTZ6KYlF77YfrIvuiRYrAqaJpGaJCdDDzKY+fN8bn/jesYd+YoD2w5RrmpCRL95frIzcnBGObjmicEsnbSSI3uP0rBNXW599Toatql3ymMtLUxDfq6gVADtSAlPVoM//rw/Opl2D8S/qRv2syesvAzw6XHcEW+CBdRCQjZ9KzF+CEhQ94O1aGsQyfuP4nP7qNagCq16NeebzR8x9bvZHE1KpVWv5nQd2jEkflnTNLav3IXP46dxhwYs/nsFb938CQFfIG/hMiYhij439WTBhKVkpGbSoE1d+tzSkylfzQy5ZLvLRod+rek6tBOtL24e0chpmsa07+cy9u0/Obj9EJp29nx5HFEO+tzYnanfzkY9SdnRm+Plj48mM+T+/jRoXTSXyKxf5/P+iM+x2PS3MSnhxd8fo80lLUp97KeKacjPFaLuhqyXSrFBqWePylOfAZ3rCFsNpOsacP9ssNcCNmNJV6nlIHO+iRz/LwP6A7gQDu9J5uWr3mXvpiSEohCTGM1TPz5Ay57NGP7siRJ+7mwPP708jhk/zSfgD5CT6UZKiaIoSE3i9wXCIjiy0nL44+PJeb9vXLiVHav38Na055g3fglrZm8gvnwsQx8cQPdhnQtN5nlxyDssnbwSTT29BjyxcjwVa5Qnef9RstJyIkamCKH7uivXqciVDw2k8+D2TP1utuGx2RlFLzhyaPcR3hvxuR45ky9y+IUhbzMm6Qui48+uuqqmIT9XcA2CrFcozQVJoTiRUTdCzo8UTQ/k/ESmP4eIfw3pnR3MfM29xw5wdEXYLgg/R/qRx64tIJzQAc5LEErBtSxVVeXRni9wNCk1b3brzfHy7KA3+HrTh1SqWSHkuD0b9xfoYy4KPreP3z+YyGsTiyemtnHRFhb/c/rXthwuO0MfGsi1T14BwF2tH2PX2r2Gxzbv1oTXJj2d5+eWUlK+WjkO704OOU6xKLTv17rIY5j1y3w0Ndz9JIRg4Z/LufTmnkVu60xghh+eIwglGoRxFAEAjmEU67ksfUj/FtA8YKl2yuM7p/H8jcwajajwN0TdAkpVsNSBmIcRCZ9EOGc6BHZibMRt4LwUEf9GoV2vmbWBzONZYS4K1a/q7o8gK6auJWnrwVM24qAbu2WTV/P0gNdxZ4UmqqmqypKJK/n51d+Z9esCfJ4T1zfhw+KpBlodVpzRDpp3vaBI+i1C0bMm4yvF4Yx2kJOpj61Nn8iujPs/HRGyWCmE4NGv7sER5cjr0+awEZMQza2vXMvauRv58eVx/DNqKhnHjDJzdbLT3YbJWKqq5o3rbMI05OcIQlgg6vbIB/img70TxVLqyx4Jnp+DUSvnCfbLQJSg0nr2F4BAifs/lEpzUSpOQ4m5LWKhB+meSEQBsphHURLeC0vTNyL14PE8f3Z+/L4Ah/PFkm9dtgN3VlGkcIvOmlkbeP+Oz/N+z07P5p42T/D69R/y/Ytj+fCuLxhe914O7dbXZtJSipcZHPAG8GR72bR4a0jWZSSkJlGsCsl7j/L1Uz9zc4P7OLjzMLN+nh/xnKWTVoVta9WrOaNWvsVlIy6maedGVG9QBa/by80N7+OJPi/zw0u/8fmj33NDnf+xYcFmw3Y7DWwbVmhaH6Sgfb9WhV7LmcY05GWElF60zA/RkruiHemAlv4MUi2kknn03QU0mIWIf1OfXYto9IzQgvBz3q1yWhpBzL16FmyxsYLfSAQrAqrxqz6gC6MVkSadGhoaOWe0g9bB7ELQtUOc0QaGpQhEcnv7vX4W/rk8b1b+3fNjSdp2EHeWB6lJ3Fkejh9JY0Szh5n4xbQCZ8YFURx/us+tPxw92V4yUjN56rLXOHYoLeLxP74yjkutVzMwZjgf3/sV7mz9YVezcXXueudGjuxNYd+WA3hzfKgBTY9hl3qsvifLw0tXvYuWT85XVVWm/zCX754bg91lx2rXP0sh9M/kivv7Ub1B1RLchdOL6SMvI+TxO8G3ijxpWvcfSO8CqDAFoRhLlCqKFc3aFAIGBsfWEmGpBBWmgW8BqEnIrM9PIdLlHETdBscGlvBkDSzlkeoRZOZb4J0N2MA1BBH7UJ5KZR4RZ9sCrEUPS6vZuDpdh3Zi4Z/L8uKU7U4bFWtWoOc1F+Ud131YZ0Y//kOeCmJRqVqvEp0HtWPCR5MN90tN47VrPiChSgILJizVwwlPwufxM/Kh7+h+Vcdi9X2qaJrkYCEx775guKI3x8uUr2eyb9N+3p2tBwXM+Gk+2ek5Bb4NZB3PZsfq3SgWhW+f+ZXVs9YTyCdXYHPYSKySQNs+Leh7S6+81P2zDdOQlwHSvwH8awjVFw/oiUCeSRA1LOK5Iu4FvdQbPvSQQgsIByLuWX2/sICjh96PlgNZ73LezbxPByIOlKrIo32D2jhB/2jOL/rnVe6n0IgO56WQtYVw94oNYWtVaHeZx7PYtmIn8RXjePy7/9Himyb8M2oaXreXHtd0Ydijg7A7T7xVRcW6eH/eKzxxyUuGWYaRyE53c88Ht7J/20GWT1kTtj8QUFk6ZTVCiAKzKf1eP7PHLsLusuXNms82Ar4AW5btYOp3s6lQozwb5m8q9MEX8KskbT3I+3d+gdcdLl3g9/pxZ3m49OaeZ60RB9OQlw3+TXpQahhupH81ggIMub01VJigp3n7t4CtGSL6DoS1TvixriF6BmGxIlJyfexllCRkuQDULfk25NMtPZ3INL22ppZJ6LV7wb8amdIDKVzgugoRfTMi6nqkeyyoKZx4ILsg9gGEEq5Dkp9f3/yDn14eh9VuRVM1qtSpxOtTnmHAnX0KPK9Os5okVkkoliGvWLM8qqqyY/Ue4wNy8wqKkOGtBTR8EdLrS4JiUQr2nQtwRjvxFGNtwOv28e5tIxFCIEQR+gCmfj8nJGvzZDxZHtbM3kjr3sZhqGcDpo+8LLDUAMOUb2dQ5rRghLU+SvybKBX+RIl/zdCIAwhLBYh7Dl2lz4ZupB26scSJ/vHnfgUc+VQYyyrTMwqir4foh8B5AyT+gKHA1GlBC7pTjCISAnpYorobsj7R3WIiBlH+L4i5B6wXgr0nInEkSnQBC9LAsimr+eW13/F5/ORkuPFke9m35QDPXf5mkUbpiil6MQm7y87w565i/bzNeLJLd6G0NCjMwEbHR/PIl3fnCVoVByklmiYL7ePC7k3YsqRg6Qu7y05ipchhpPu3HuDXN/5gzJt/kLStAM2e04g5Iy8L7J1AqQiql5DyX8KKiCq5joYRStQ1SPtF4PkXKX0I58UI2wVI/ybwzgMRhbR31w2Vd3Ywpry45D6UJKcW565CRm7aegA84zmjDxWlAvoDriCj5wHfImT6E4j411Bi/gcx/ytyF398NCnsdV9TNZK2HSRp20FqNCo4FPTy//Vj55o9hboMLFaFu9+7ia5DOjL3t0WIMpQqLgk2h43357xE6sFjQNFm1pFwRNnxefy63zv4gme1WbA6rGxdvjOkOpARiqJQrmoCL131Lu4sDz2v6cLFN3TFZrcx9u0/+eGlcXkKiT++Mp6bXryaax4fXKKxlhTTkJcBQihQ7hdk+pO6jgeAtREi/s0Tetel2Z+1JsTckfenLKXUq8ZH36brgPjWItPuDsrcGmfQGRMFqLo2d/QdCN9CpGemvthqaIALc5OcbJyKt7B3algg5n44vhpkrjpkAXgmImUGIvGLYvWSftQ4dtlitZBxrHDp1l7XdmHDgs1M/DyyQqBQBF9teJ8ajXTFxmZdGuP3F+dzLVtcsS6GPNAfZ7SDF698N6J+eVGJSYzhhx2fkpaczrTvZrNr3T6qNazCuHf+LvDhYLFZiEmIpl3fVrxz62d5D8+V09Yy8qFvuOudm/jhxd9Cxqf6Vb57bgytejWjcbszJ6hVKuqHQohvgIFAspSy0BWB/6L6YSSkdIMMIJTYM9Kf5p4OmS+Ddgyw6AurnsWg7ShmS1YoPw5hqRUydu3owKAIlRHR6K4Lgz8eERd8kJRFsoUDYh9Bib4VGdiBTH8mn0SwJLJRdyAq/IU4KUpFyqDGCuhKl/kWSX96ZRy/vvFHmHFyxjj5PfnrkAXOgvjyyZ8Y//7fEUP7Hhx1BwPvujTf8T/y2zt/F6ntssbhsqNYFbxuHzIofnUqtOvXiodG3Unl2icKnj/R5yVWz9xgeLzVZqHr0I4MeaA/5auX47YLHjR8mCgWRS/TaDA+xaLQ+/qu3PrqdXnZuaXB6S719h3Qr5Ta+k8hhOuMGXHpWw7pjwZDEv2AB3LGgbYzwhl2EBEU4pyDUWzNDMYeyZ8p0GfpJ3/lHGBpCLEvRVg3OFWE3kfEl08LVJiEEn2rfrS1AUr5sZD4vb6WUdDMXNjCpIWlfxPyaF/k0YHBn75I/4mkkyEP9Kdc1UQcLt1gCyFwRNm596Nbi2zEAa5/ZmhEvXCAr576OSQ+evjzkRfQzza8bh/uTA9aQDtlIw564tNtTR7kpavexef1o6oq6+dviXj8iDeH88yvD9O0c2M2zN+CxWr8vdTUyOPTVI0ZP87jpgb38fIwvd/TSakYcinlPOBYabRlcvqQWSMJ9/96iGisRBRII3+xAtF3GJ/juhrjBUol2FcgdJuzP6LCHwhnT4rvD7dS+FdY6jHfzisJr8DkAOdA8M5GS38aLfsnpJaFVA9C2l2g7iukaRUstU/8qmUhj90YzJT16D/qHuSxG/XarOgLeJ+vfoebXrqGFj2a0vPaLrwz80X63drbsIu9m5P48skf+eCuz1k6aWWecVYDKh37t4k4NL83wIYFJ4zV3o37C6zMcz4T8AXwefwsm7Kab5/9FU+2t0CXSof+JzRZYhKjT+m+qX6VhX8uY3D8TVzmvI5nBr7OgR2lXwC91ApLCCHqABMjuVaEEHcCdwLUqlWr7d69BWTGmZwWtOTeoCUZ7LGgz1zzG1knWBtBYJ3B8S5E3NOIYP3Q/EipItMeAG/QTy5sIBXy3gBORsSBcOg+e3s3yB5N0XzjNoi6GRzd4PjNRTjerme8ymwQ9mBptSa66JX0obt0XKBEgaOvXgs1Ugp+bnu2Fijlfzlx7TnjkBmvASer7EUh4p5DRF1ZhHGe4N9vZ/HpfV8T8AdQAxrOGCctezTlqkcG8dzgt5CqFrFCvNVmoeYF1WnTpwVDH+iPqmrc3PD+Upnhnsu4Yp38lfYDV1cdYaiPnlA5nnGHvsr7PeAPcE31O8mItLZhsyAEEYs4n4xQBDEJ0Xy75SPiKxS/Ktbpdq0UipRytJSynZSyXcWKFQs/waT0sV2I8UduA2sTwAUihtxCw3qtUQOXhFBAGLuDhLCgJH6GKP8LIvYJRPzrUO7bYEELA2QGaCl6glT2lxD9KIUXg7YAdsj5CXxriOzOyY9P78tSHaIfRFT4Q5+py0xO+OXdoB0H73SMjbgSHJsdnAMQiaNDd2spGPv4PaAlG2w3GKXHh6ZpZKdn88l9X+N166nloMczr52zkReueAtPlieiEQc90WX3+n389cm/jGj+CIv/WoGZGKan/m9bsZMmHRvl6YznYrNbeeaXh0K2WW1W3pr2HDEJ4dnWDped/iMu5qsNH2B3Fc0tJjWJN8fLpC9nlPgajDCjVv5DiJj7kL65wWo2ubgg5i6UmP/pIYnqfrA2QVhrIQN7kDk/Eh7JIsDZq+C+bM3A1kz/RWpIXITPVE/GA9nvEzm6JXe7CgSrtWePAqV8UH62MNRgLPi7SPUw+JYSvvCqBQ2ynfBEKhuUn4Cw1jMu/WZrpStUypOuUzjBVrCE6upZ6/n43q84uP0QNqed1r2bY7GGP/w82V4Ug+2RCPgDBPwBRj3yXZHPOZ9JqBjHo71exOfx5WnQWOwWGrSsw4g3hxtmbzZoVZfxKd/wxWM/MPlLXZFSahqX3Nidez64BZvdxhPf3ss7t31WpApHPo+fbcsjrUuVDNOQ/4cQtoZ62GPmO3pUhlIBou9CuIYG9zfNq+QOIKx1kHGvQsazQSEoCVgRiV+Ea48UhJpE0aNRIsRwK1X0OqNhDwNP5MIOEfGC+ycKLjjtQ5/55yuW7eiJYmsYuVl7Z7A2A/96cq9D4iDlUBXevnoSFWosZcgD/WnSMbSNHWt289zlb+YZAW+OlxXT1yKNIlJE5PcVi1XJm72bGJOd7s6T5s29u6pPpc6FtWhaQGk7i8XC/z64lTveGs7RA8dIqBiHK+bE30CPqy+iSt1KfP/CWFZOW6vHOkVwY9mdNuq1qlNKV6RTWuGHvwI9gQrAEeAFKeXXkY43ww/PLaSWjcz6DDx/6GGLlmoQ8xiKy7h+pNSO6zIESgWErTFa2pPg+YtTShayNATtEMjCY62LhqJn0ap7iegLVyrqC5rCBVHXI6JvjShrm4uUPmT2j+D+HVXTGPuRhXGj4shJDyCEwO6ycf9nI+h784k3mleufo/5vy8tUpq8I8qO1WYlO73o1W5MdApKKrLaLfS6titPfHffKfdzaPcRxr/3D1uX7+DIvqNkHc8OqXAUHR/FN5s/LFLN0JM5rcWXpZTXlUY7Jmcn0jtL90fnzpbVA5D+FBoSxXVCbVBKicz6GLK/Ci5yqkhrHdDSKZoRt3Oivmh+FLA2AG8hUSSGRJp1C3B0B+98UCOkaEs3SuVwvWvDQ/3rkJkf6+GI1oaI+Df45fUtjP3gr7zMQSkl3hwfnz3wDb2u6ZIXbrh3U5KhEc+VrVUUgaZJVFVjwJ19kFIycdQ0FEXB7w+UOOPxv0ZB9yngU5kzdhH3fHALsYkFa+UURtW6lbn/0xEAuLPcvHL1+yz/dw2gh5vWa1knTx63tDC1VkwKJ/N9DMMWM98L3eSdDtnfAN7gzNmtJwfJSCJPSlDfxaH/OIfoBjusOIYdEXM7xD5E0bVXrLo7xnU9xs4IK8J1JZT/g4ja7UqVIvUkfcuQqcPBN09/a/DNQx67kdS9UwzTv4UQ7F6vP5TSUtKpe2Etwwo6Po+fO98eTsdB7bigY0O6XNGBiZ9PY9IX0xGKIBBQadyhAY6oosefm0TGardw7NDxUm0z9eBx1s07kUcgpWTzkq08fdlrpdqP6SM3KRwtghDQSdtl9neE+8LVYHifg9CwwmDIoXoA3XD79f0Jn0PagxDYGvTLC4h9EWFrgbC1QNouRGaNBt98jBdELeil1vogYp9CWCqg2VpBxjOELKLGPISwNUQAWtSNkPMzoQ8rFyLm3sLujH7dGa9h9KC79n8bmfJ9nbDjVb9K5vFs7mn7BHs27EPTpOGMXFM1Pr7367wZuRHbV+4yXBQ9G+gl93E7G6hIDilE8TXNmS1qlfWwIqKpksp1KhX7vJxMNzaHFZs9PHrqj0+mEPCFPswDPpW9m5LYtW4v9VrUDjunJJiG3KRwlKrGxlw5qVKKFmE2I+x6LcycH9ANe0CfeQd2EmLcPf+CiEGpMB4Z2K+HC1ob6noweU11QJTrgHbsJvCtJNS/7dTDHm2hkQdK1OVIx0X6G4MMgLM3wlL9RJuxjyJRIWdMcIsVoq5CBjYhM3chXIMiKkwC+kPHgMo10nFG2fHki2RQLArVGlThlavfJScj1PgLYaxuHMmI5+KPUGG+LOkl9/EIK3EG3WSVyeERVoLkrDTmzigH1z41BGeUg63LdzDnt0UIAb2u7UrDNsaKpOvmbeKDu77g0M7DKBYLPa+5iPs/C60hemDbQcMFaMWicGRvSqkZ8rPzUW5ydhHzMOFZkc7g9vybLsbYTaEgYv6HqLQEUf5PRKWFupE2cte4xyOlH2GtibA1CzHi+REJnwULaNj1sSgVEAnvhxnxvOMtFRBR1yGibwwx4gBCWFHinkZUXgYVpoKrH7jH6r7+7FHIo4PQcsYathscjPFmJYHrnh6K3WkjOj4KZ7SD6DgXezcnhRlx0I14sbMIpeSiy9sX75wzwO1syDPiuThRuR1jfZOypHLtitz/2Qiuf3ooX/7fTzza60V+f38i49+fyMPdn+P7F8I/+31bDvB0/9dJ2qobar/Xz5yxi3j5qndDjmvRs5mhDK/fG6BB67qldg2mITcpFCVqMMS9cmIGrlSFuFf07fkQ0beDUg7djQL618sJcS8jhE3/sdZBKPGRZ+8EIsgChCKUWJTEkYhKixAV/0VUXIBwXlLSS9TbFE6Euk+v0iTd6G6YAOCFjFeQWgQViugRhPvuXRB9G9c/fSVjDozmud8eoec1F+Hz+PW6kREobualYrVw++vXl7ie5+miYoScgUjby5I2l1zIvPGLeeqyV/njo0l4c7x5YljeHB+/vfMX+7ceCDnn9w8mhq1/+L1+1s3bzKFdJ8orDryrD1FxUSHuL0eUg0uGd6dijfKldg2ma8WkSChRgyFqMFLK0JJn+RBKIlSYiMz5Vdc6t1RDRN+EsBlUVrG1CsrdntxRpWB2adEQShxQ/FTnSEjP5AgPEqt+Ta4rwscQfZsecpnzo57BKrVguOKdAMQmxtDmkha8dOW7BWZjFhWhCBQhsNit3PDMUKo3rErLHs1Y9u/qsyYFP4UoKhsY7RSM69GWJVO+nlXgfqlJlvyzkpqNT7zJ7d203zAKxma3cmh3MlXrVQYgrlwso1a9zfcvjGXZpFVExbm44oH+DLr70rBzTwXTkJ8n6ItlXsAR0dCWBoW1LZQ4RMxdEHNXwcfFPoE8tipoNDVyVQpF3AundfyFE1xgPXkhVQgi/bkIoSDiHkfG3gfqYVAqhxXQ1jStxFV6rHYrikUghKDZRRfQvNsFICVx5eNI2n6Ih7o+x441u0KMeCR/+5nia5qH+MgBPFj4mrO37mUkhEXB5gh1jzTt3Jity3eGxIeDPiuv3bRGyLYK1crx6Jf3nNYxmob8HEdKDZn9he7PldmgVEbGPoXiOrtVhYXtAig/AZk1CvzrwFIHEXMPwt6qbMcVNRjpHkeY/16quv5MQecKF1iN/Z4Wi4XazWqyZ8P+Yo3HGe1gxJvDEULQ9KJGNGhVF5/Xz+O9X2T3+n24I9SzLMyIxyRGk3U8u1hjKQ6zRS2QnFNRK5Hwe/ysmL6WOs1r5qXwD32wP5O/moHqV/MijhxRDnpf14XyVYuf6HOqlJr6YXEwMztLDy3zo2Dsdv6wPyci8VOEo3tZDeucRsv6FLK+QM+H16e2IuFDhNNYaraorJ27kWcGvIHP7Qu6qEAoCrHlY0hPyTCMpoyKczE++euQ0LY/Pp7E10//UiRdj/MJxSLoNKAt6+ZvxuawkZ2eg68UXFVGCEWEualsDivXPDmEm1+8GoCk7Yf48okfWTN7A9HxUVxxf3+ufHgAFsvp0NUPjitCZqdpyM9hpPQhkzuEizQBWC9EqfD7mR/UeYIMJIF3TlDL/BKEklAq7e5Ys5tfX5/A7g37adC6Ljc8M5T37/icTYuNqypd+fBA7n4vVKb3f+2fZPvKXaUynnOJ6g2r8vy4R/j22TFsXLyVrGPZRZI1KE1sDis/7PiUCtVLb6GyOJzWFH2TMkLL0F/5jSisKIJJgQhrDbAOL/V2G7Sqy3O/PRqyLSrOeAHQ5rDS56Yehtv/azijHdz04jDu6/h0ocWSI2GxWWjRvQnr5m4qsbiYlLB61gb63Bj+uZQlZvjhuYySoCfbGGEtQKXvP4iUEi3nd7SUi9EOt0BLvQbpW1nWwwLg8v/1DQsfFAIq1ChvmDAy4I4+huGGQggc0XaEEBHLk51LOKMcNO3ciJtfvoYfdnzKDy+OK7ERB73M3j0f3HpqCpFSEh3hwVuWmIb8HEYIq175PSyG2YmIfaQshnTWInO+0YtOq/sBD/hXI4/divStLfTc002ngW0ZfN9l2J02omJdRMW6KFetHK9NfMowgueSG7tz0eD2OFx27E4brlgnMQlRXPt/V/DwF3fz0cJXqNG46jlf2s2T42X7qt3kZLjRNMmB7SUvkSYENGpbj2WTV2GxldzsWR022vVrVeLzTxemj/w8QMv5A7I/BTVZT2mPfRLh6FjWwzpr0NcSOupRPSdj74JS7tszPygDjh5IZePCrcSWj6Vlz6aFLprtWreXNXM2snzKKtbN3ZSXyq+qKgIREucsLMHFu7MjzLxYKBaFctUSOLr/1MoC2xxWouOjDEu8FQWLReGDBa/QpGOjIh3v9/mZ/v1cZv26ALvTxoA7+3DR4PanFF5r+sjPY5SoIRA1pKyHcfaipUZeSwhErqZ+OpBSkpWWHZxN2/F5fIx9+y+mfjcbTZX0vr4L1z99ZZEiH+q1qM0/o6ayft5mfJ5Ql4M8yWJLVRZeQe8sRVO1UzbioKfFl9iI2yw88NmIPCPu9/n59c0/mfTFdLw5Xjr0b82IN4dTqWYFQH+YPtnnFbat3IU3R9cTWj9/M31v7cV9H99+ytdyMqYhNzn/UcpF3mc5c3HNa+du5IM7PufI3hSEIug6tBNH9iSzY82evDC6CR9OZvmUNYxc8Vahfm53todp38/Nq3hTKOfgbPxMYHfacEY7yUg1LrAMejHrCzqcWHd69ZoPWDFtbd7nNve3xayesZ5vtnxEbGIMS/5ZyfbVu/OMOOhl+qZ8NZOhDw6gWv2iSSQXFdNHbnLeI4QDooZjuJYQ88AZGcO+LQd4duAbHNhxmIBfxe8NMH/8YjYv2x4SC+33+jm06wiL/ync9ZiekoFyjvvBTxdG+u5GOGOcPDPmYX7eOwpHhALKVruVlj2a5S087996IMSIg/7W4M7y5KX7L/93NR6DZC2/L8B7t48k6RT8/UaYhvws4bjbzZgN6/hm9Up2Hkst6+Gcd4jYRyH6NhDRgEUX/op/B+Hockb6//2Df8LcHwG/aliX053lYeuyHYW2WaF6ubBK8CY6Nqc1omHOT7POjbjo8vasn7eJQMDY/Va1XiVe/OPxvN93rd2L1eBtyev2sWmRLmkcXzEOq8FnIzXJhgVbuKfN42xeGqEyVQkwXStnAXP27ObeyX8DAlVqvLNoATe2aMXT3YxjVb2BAJO2b2XZgSRqxScwrGlzKkZHn9lBn2MIYUHEPoiMuV/XdxGuM6rpsm/zgSKXZHNGOahUu2LE/Yd2H+HQrmTqNKvBra9ey+jHfwx7SPzX8Wb7qFpPLxKReSwbq92q1870n9BGcUQ5uObJKwDYsXo3qt/YkNdsXD0ks7ZKvcpoWvhnabFZOLjzMK9e9wEtezTDYrUQMGhT0ySebC8f3/slo1a8fSqXmYdpyMuYHL+f+6b8gzsQKr7z8/q19K5bj041aoZsz/B6GDL2F45kZ5Hj9+OwWBi5Yik/XnEVratWO5NDN2RD8hG+WrWCvelpdKpRk9tatT2rHjJC5JaXO7M069KYrct34PcaFIE4SaPLarfS+7rwNwV3todXhr3H2jkbsTls+Dx++tzYnSe+v4+RD37LscNpBRYYPt+w2a2oqhbxeo8eOM7Pe0aSWDmB7Iwc3rzxE1ZOW4vVZkGxKNz13k207q0rc+7ZGFkDp3az0L/BRm3rUatJDXat2xsimqX6VXav38fu9ftYOmkVVetXZu/mJMO3LoCdq/egBtRSifk3Dfkp4FNVLEJgUUruoVq4b6/hdnfAz/hNG8IM+cjlSzmQkY4vOCPwqiqoKg9Pm8zsm26POMucv28P369ZTZrHTd/6Dbn+wpZE20u31uPMXTu5/9+J+FQVTUo2pSTz28b1/HPdjVSLLT2p2dJAk5KJ27YwduN6VCm5skkzhlzQFOspfJYFMeSBAUwarYss5YYJOqLsdBrYjsO7j7Br7V4QguoNq/L0zw8QHR/+8Pvkvq9YM2cjfo8/bwY+85f51GhcjbEHv+TowWOsmr6OkQ9/S3ZawbrfQhH0vq4rCyYsLRVp3bKgSedG7F63l8wI4l+KReDO8pBYGaLjonjlrydJS0kn/Wgm1epXDpllW23GptBis9C4Xf2QbUII3pz6LB/c+TmL/1mJpmr6DD2fvfZkeQoVSLM7bUX25ReG6SMvAdtTUxk27leajvyIpiM/4oEpE0n3lEyiNDUnG7ff+LV4h4GvfNL2bXlGPD9HsrI5lGW86v75imXcPfEvZu3ZxarDh/hg6SKuGPszORH6LQmalDw9axqeQAAtmJvg1zQyvF4+XLKo1PopLR6eOpmnZ05ncdJ+lh1I4sU5M7njnz9Om3ZH+aqJfLbsTboM7UhMQjSValfk1leu5elfHuTTpW/y456RDH/uSjRV5bnBbzH6iR/IPJ6Vd77f52fOmEX4T3Kh5BY+kFJSoVo5Lr25J/1HXGJYlSY/dZrV5M53buThL+8mJjEaZ4wzKJdb9iahKB4vXWmwK+/NeYnoBOM3rNjEmDxd8FwSKsZTu0kNbHYbWWnZjHrkO66reRfL/11tuN5gs1tpe2lLw7afH/cYf6X/wD3v34zTVbzCHla7lX639y41995/YkauSYnb7yfKZjvlG3fMncOw8b+S6fUig21P3bmdXceP8c91Nxa5/QOZGdw3+R82Jh+JGBV23ODhYI8QXyyR2C3hH2eax81HSxfpM/cgnkCAg5kZjNu0gZtbti7SeAvjUFYmmb7wmZ0qJfP27SmVPkqLDclHmL5rB5587ix3IMDygwdYnLSfi2qenpDEavWr8PxJOiu5jHroO5ZMXJkXrvbHx1NY+McyRq97D4fLgd8bQFONfbhpyRnc1/H/eHvGC0THRXHjC8NYN28T+zYl4cnxGhabOLDjEDfWu5fb37yB8Ue+5uDOw8SWi2HXur08eekroaGKBvLspxNXrAtPlifvzcVis4AMyiyoGs5oB7Wb1qDPzT2xO2x8v+0T7mr1GBnHsvB7/FhtFiw2K49/d1/I3+PW5Tv44aVx7Nmwj1pNqrNv8wGOHU7Lc48oVgWhCCwWBYvNCkieHfsIzqjIRtrusBFfMR5hKZ5dadyuPne+fWPxb04EzmtDLqVk1IplfLFyGTn+AOVcTp7s0p2hTZqVuM1xGzfgDagh32u/prEnPY1Vhw/Stmr1iOfmoknJtePHcjgrE7WAGaDjJN+Z2++nfbXqIa4VAIsQNK1YiQpR4TOTNYcPY7dYQgw56IZr5q6dpWbIY+32vJn4ycTY7WxKSaZ+Yjkc1rL/yi1J2o9q8FaT4/ezaP++02bII7F3cxKL/1kREs4W8AVIPXicT+/7GsVqoXrDKlSqXZFDO48YtrF7/T6+evInHhx1J65oJx8tfJVvnx3DnLELOXrgWNhCns+tz+y/eeoXWve+kLrN9Wtuc3ELvlr/AW/e+DG71u1FsQiqN6hK0vaDqP7T73svXz2R9+e8zDfP/MqaWeuJLR/LsEcvp3H7+kz+cgbHk9PpMrgD3Yd1ynONxFeI45vNHzHt+zmsmbWBKvUqcfk9fUNitdfO2cgzA1/Pk/5N3nc0rG8toGF32ek/4mLqNKtJtys7EVc+ttAxdxrYplgPOrvTxmuTn8buLD3XZtn/VZ1GRq5YysjlS/MWElNycnh29gyi7Xb61i+ZqNTW1KN4VeOq5buPHy+SIV+ctI80j7tAI+60Whl6wYkHztQd23l02hQUIfLOsysKNouVeKeDT/oN5EBGBisOHaBCVBSdqtfEoigkOJ2G/Qgo1UXIOIeT7rXrMHfPHvzaCaOhINiXlsY148cigOd79OKqpmVbJSbR5cJmseA/yZg7LBbKGzwMTzfblu80jAf3un38++1sAOwuO0Lor/p+X/j3z+8NMPOX+Tw4Si8v9/bNn7Lor+V4sr1hx4ac5wsw48d53PHWCaXH2k1rMGrliWiK12/4iL2bkkp0baD744UiCqxVCvqawe2v30DFmuWp1qAyq2as5fDuZOaOXUjzrhdw/6cjIp4bFeviivsu44r7LjPc/9lD3xRNv11KqjesyoA7+xR+bBBXjItX/vk/Xrji7bxan5EKfjhcdm566ZpSF946bw25qml8sXJ5WDSIJxDg+dkzqRkXT9OKlYrdbovKVZi6c3tYu1JKLqgQOWQsP0eysgqs4BJls9GgXHluaaXPlg9lZvLwtMkhrgDQJwFvXXwpl9ZvwKvz5zB24/rgYp0gzuHgl6FX07JyFSpERZGUkREyY3ZardzUolWRxltU3u1zGXdP/IvVhw9htyhk+XxIJAEJAb/+R/TCnJnUSUikXbXCH3ini771G/LS3PA6jYoQDGzU+IyPp3z1coUKXBWlgEKuqt/OtXtY+OeyIhkuqWmFLna27t2cxX8X/lCI3Idu3CrVrsDxw2koikKXoR3pOKANP7zwG4d3H6FSrYrc8sq19L6uKy9d9Q6L/l6RZ/hXzVzP3a0f4/sdn1KxhDrgRa3MZLFZqFxA6GckWvZoxm+Hv2L1zPX4vX6OHT7OF4/+gFAEAZ+KpqrUblqTO9+9ifZ9WxW7/cI4bw15tt+HN2A8c07Jyeaqcb/St34D3ru0P0ox/OZDmzTjs+VL8QYjM0CfybWoXIXmlSqHHS+lZPPRFPyaRrOKlbAqCq2qVEUzeBezKQqdatTkxhat6FWnXl40zN/bNhu6LWwWC5l+H1N2bmfcpo14VTXPheIO+Llz4p9MHX4LP14xjFv/nsDhrEx9Rq9Jnu/ei5ZVqhb5uotCnMPBL1dezd403c309MxpYS4dTyDA16tXlKkhj7Hb+X7wldw16S/cfj8CgdWi8MllA6kYdeZDJVv1akZ8hTjcmSVbMAc9k7HjgDYArJu7qchFmB1RDrpdaSyw5vf5WfjHMvZvPYgrxokaUPPCJ53RDjr0b82iP5cbxkobcezgca75vyu45aVr87b1vrZryDGHdh9h0V8rwkIK/d4Ab9/8Ce/MeLFIfZ1MbLkYMo5GTsEH/R7GJETT4bKSuRvtDhsd+7fJ+73b0E4s/mclSEmnQW0pV+X0lYA7bw15jN1BrN3BMY/bcL8nEGD6zp1M3LaFyxs3KXK7cQ4Hf117A6/Nn8ucPbuwW6xc1bQZj3YOj/vdmHyEOyf+RbrXg0BfqPyo30C61qrNpfUaMH3XjryZvV2xUC0uli8GDsZpDY04yPT68BssdAVUjSyfj983b8QdCI1m0KQkKSOdXcePUS+xHNOH38KW1KNker1cWKkyLlvBUQ2nQu2EBI5kZxn65iX6G0ZZ07pqNZbcfjfrjxxGlZIWlaucttDDwlAUhXdmvcCNde8tcRvxFWL53we3AJBQKV5fICwkScgZ7aDblZ1o0b1p2L70oxnc3+lp0pLTcWd5sDn1iJbKtSsSHR/FoHv60v6yViz5ZyUU0ZAH/Crj35vI4P/1Y85vi9i2Yif1WtSmQeu67Nt8gArVdU0co2QbgDWzNpK07SA1GkXOl9A0jS1Lt+PJ8dG0c6O8hcqrHhnEz6/+HqJ9YnfaiIqLIjs9Gymhcfv6PP3zg6cc152SlMryKauxu+z0GNbJMJS0tDlvDbkiBI9f1I2X5s0Kc0nkkhPw89vGDcUy5HP27ObNBXPZnXacqrGxPNKpi+H5noCf4X+MI9174ouT7fdz18Q/mX3z7bx36WWM2bCOn9avxe33M6BRY+5q2z7MiAP0qFOH79asIuckY60ogu616jBh88aI9yA3tFEIQRMD14+UkgX79vLn1s0oAoZc0IzONWqecnRPkwoVw3zQoD/Mutc2LlB8plGEKPW3kpJSuVbFEkeHWKwWvtr4AXHl9IW5iwa34+N7CzZGiiK4+IZuPDjqTsPP+ovHfuDI3pS8mbHfo/8NxZaLyfOfSylJqJxA8t6Uog9WSkY0fxiv2483x4uiCDRNYnNYsdpt2OzWAu/Bw92fo13fVuxev4+GbetxzeOD8wz7rnV7ebr/6+Rk5iCELuP74Od3cskN3bnmicFkHc/ir0//zUuaGnxfP25/4wbSj2ZisSp59+9UGPPWH/z40jiEIlAUhQ/vHs0L4x+lfb/SCSqIxHmvRz5p21beWDiXgxFmgR2r1+DXK68pUltz9uzmf5P/DnkwuKxWXuzRm2HNLgw5duK2LTw1cxrZJ8VqOywWHup0EXe17VDka5BS8tDUyczctTPPmLusNoY1bcaLPS/m02WL+Wjp4rBFzXiHk+V33FPgTPOpGdP4e9uWvBm9y2rj6mbNeaFHyQsNewJ+3lgwl1/WrwsZk01RSHS6mHzDTZRznVjsCWgaAU01fIj9V9A0jb62ayIaMVesE4fTTk6mJ0Tt0OGy0+/23mHSqDtW7+b5K94i83g2AgzDEB1RDl7643Ha9gmPkx4Ud6Oh6BPAmANfUL6qPnse89affP3Uz0W+TkVRQBAxG1MoIi/UMBK5hZEVi4Ldaee9OS9Sr0Vtrq1xl17EOh8Ol53Plr9J7aZ6Yp0728PRpFQq1CiPK9pZ5HEXhW0rd/JIj+fD1iac0Q7GHvySqNiTRduKTyQ98rKP/j/NDGjUmDk3jyDeEf6huazWYoUivrVwXtjs3h0I8M7iBWFfvFR3DgGDGalXVUnJNs5ECzkuEMjz8Qsh+KBvfz7o259+9RsyoGFjRg24PM/Ybj92zNCH/vhFXbEqCmuPHOaOf/6gx3dfcffEv9iYrIewrT1ymL+3bQ5xy7gDfn5Zv5atqeHhWSejScnaw4dYtH9fSFLTfZMn8tvGDSFGXKCvL0y6/oQRz/R6eXjqZJqP+pjmoz5h0K8/sj7ZOLzufEdRFFr1NP4u1mtZm+fHPcaYg6O5+eWrccY4ccY4sTltXDy8O3e9e1PYOQ1a1+XnPaP4cP4rPPLVPdgd4Q9Jb46XPz6ebNinFkFACmDe+KV5/5/+wxzDY8pVSQgTjXJEOUCJbMRBXxhVrEqB2um5DyRN1fBkexj54LesmrHesAyc3xdg0pcz8n53RTup2bh6qRtxgBk/zjPUvBGKwvIpq0u9v/yct66V/FgVfSHrrol/okmJV1WJstloX606Qy4I9w9GYk9amuH2Y243XjUQMqPsWN3YPRFls9GlVngdxlz2p6fz5Ix/WX7wAACdatTkrUv6Ui02jj71G9CnfoOQ41Oys5m2c3vYRM4iBOuTj7Bo/z5G/PMH3kAACSRlpOvp+ldcxeKkfYZuJ7+mMWr5Ej7sNzDiOLelHuXWvyaQ4fWgCEFA03i51yW0q1qdhfv3hvnGrYqCy2oLCe+7/e8JrEs+gi947MaUZK77fSzTht9y1qX0nwke/PwuHuj8ND6PD2+OD0eUA1eMk1f+epJKtXS32NWPDeaK+/uTsv8oCZXiCwxjE0JQv2Ud3JluPcHFIDrl2KE0vn9xLCunraVCjfJc9fBAmnZuTJ0La7Ft+U7Ddqd9NxupaVx6c0+SthnLsR47nMZrk5/m/RGjSD14HIvVQscBrVn+75pCF3WtVgsD7ryEf0ZNK9Ki7ZZlO8g6nmU4i9dULWyWfrrwe/3G45US3ynUGi0K/wlDDtC1Vm3m3DKCv7du4Zg7h4tq1uKiGrWK5QuuHhvLrrTjYdvjHA4cJ2VVXlChIv3qN2Tqzh353BZWmleqTI8IPmK338/lY34i3Xvii744aT9X/vYrc28ZYZjVuTvtuOGiohrUOll16GCIsZbobxGvzJvNoEYXRLzWWXt2I6U0vD+qpnHjH+M5mpMd8gB5bvYMnu7Sw3A8fk1j89HkvN83H01hY0pynhHPO07V+GndGp7o0j3i2M5XajSsyvfbP2Ha93PYvW4vDdrUo89NPcKMtd1ho3qDovv2pZR4ssIX/W1OGwe2H+K3t/7C5/UjxHaWTV7Fg6PuZMQbN/Bkn1cMjeOO1bvZv/UAP70yPuLMOaFSHBM/n0Z2uq75ogZUlk1ajTPGWaghL1c1kTvfuYkNC7aQtO1QoaGX0fEuWvRoaqhe6Ix20nlQmCfitNDtqs7M+GleWJimGlBpdxpCDvNTKoZcCNEP+AiwAF9JKd8sjXZLm4pR0dzeum2Jzt2WmkJyTrhLxGW1cn+HzmT6vHy3ZjUzdu2gXFQUt7Vqy7uXXkb3rVsYs2Ed3kAAu9XCxiNHaDbyI7rXrsPzPXpTPd/M86vVK0KMOOjui+MeNzN27aR/w/BagXUSEsKMIegz8gsqVGTcpg2G17P5aApfDBzM6wvmGu73BAJkeL3EO8NfQZcc2E+O3x/2FuBTVVYfPmQ4y7cpFi6sVAVPwM/sPbtZvH+/4UPCr6lFcuucr8QmxnDlQ5HfhIrLjtW7eeqy1/PS3XOx2CzYnXa8Od68FHUpde2Wzx74ht8Of0X3qzuz6M/lhi4Lb44Pn9sfVlIulyYdG7F65voQo+bJ8aJqGkIIwweEUATOaCdP//IgDqedDxe8yj+jpjJnzEJcMU6i4qJYNWNdiPvCEWVn8H2XUaF6ea56dBATPpyU16cz2kG9FrXoOvTM1K9t3bs53YZ2Yv6EJfpCbjDV/653byKxUvxp7fuUDbkQwgJ8BvQBkoDlQoi/pZSbTrXts4WxG9bz9Kxphl/ZrjVrM/SCpgz45UdScrLzjOryA0k80LEzd7XtwODGF3DVuF9Ze+Rw3v6Zu3ax6tAhZt10G7EOPUTqj83Gt8ynqiw9sJ/Ze3ay6/hx2lerzq2t2lI5JoZK0TH0a9CQqTt24MmXceqwWrmrbXum7thOhi88kSPB4STR6YoYKGERIqI6YrrHYzgR06RkzZFDhv56h9VCh+o16PjVF2hSQ9VkyHjzjrNYaHWWRJKcD3z3wlh87vDPXwhIrBRn6BrRNMm+zUk8/fODLJiwlKnfzWH5lNVhxregBUl3tsfQXxzwBnBGO8IyHy1WhUtu6smIN64noaJu9FzRTq5+bDBXPzZYP9cf4L3bRzF33GLsTl3Gt/d1XbnhmSsBuPWV67iwW1MmfjGNnAw3va7twsXDu0dUNixthBA8/t299L2tFwv/XIYzysHFw7tTu0mN0953aVxhB2CHlHIXgBBiDDAYOC8M+fbUVF6cOzNiRNS65CNcM34MBzMzQo5xBwJ8uGQx1zVvyc5jqWxNPRoyc9aQ5Ph9/LFlEzcF9U6OuY1j3gF+Xb8Wie4y2ZB8hLEb1/P3tTdSMz6ety/pR9WYRfy0fg3ZPh/NKlWib72GLNi3l0GNL2DC5o0hmaguq5URbdqxJGm/oRsE9Jl+pGiX9tVq4NPCz7FbLOxPTwuLnlGE4I3efXhyxlQyDR4q+Y9z2WzccGF4FIVJydixapdhFrHVZiUmwTi+WfUHiCsfi6IodL+qM92v6syg2OFFzux0xTipXLsidqct7By7y8aQB/ozcfR0At4AUmoEAhqD7+3HnW8XLDpntVl58of7uePt4ezfepDDu5PJPJbFhgVbaNGjKUII2l3aknYGaoVnCiEELXs0o2WPkus5lYTSMOTVgfz5r0lA2LuMEOJO4E6AWrXOrCjRqTB+8wbDZJxckrOzOJKdZbjPblHYkHyE/RnphrMXdyDAhnxRGolOp+HsGSCQ73y/ppHp8/HOovl8fNlAbBYLT3TpxhNdurHhyGGG/zmez1cuI6BJQFI9Lp4DGelYFYWAJhneohUj2rRjxq4d2BRjQ14vMXLB4orR0dzTtgOjV60IiXjxq6rhA08A8/btKfA+ArSvVp23L+lHOVcU65OP8Pr8Oaw7cphEl4u72rZn+IWtzmhVn/OBag2qkHowfF1H0yRXPTqId28bGWJsLTYLjdrVz6sGn0vv67sy/Ye5IYUxbA4rCIHqC+S5boQicLjs3PrqtcwbvzisX0VRuOrRQQx//ipWTltH5rEsWvRoWqy0+IAvwJs3fkJ2eg5+rx+r3UqdZjV5e8bzpyUa5VygNMIPjf6ywv6epZSjpZTtpJTtKlYsvpZBWZHl81KQ1E9Ba+oBTaN8VBR1ExINDZDTag3RZ+lVt37YMZHQpGTh/tCiFFJK7pr0NxleL9l+P141gFdVOZSZyYs9ezPh6htYccc9PNW1B4oQdKpR0zBpJ8pmY6DBQujhrEyWJO3ncFYmD3a6iM8HXk6DxHIowa9ApHthURQUUXBImd5+FjXj49maepRrx49h6YEk3IEABzMzeXPBPN5fsrBI98bkBDc+PwxHVKiLzBGlK/z1GHYR1z89FLvTRnR8FI4oOw1b1+X58Y+FtXPXuzdTr2UdnNGOvJ86zWvxyeLXadS2ft5nqyiC9v1bExXr4s1/n6VclQRcMU5cMU4SK8fz+pRniE2MwWa30WlgW/rc1KPY2iZv3/IZxw4dx53pJuAL4MnysGvtHn5+ZXyJ79O5TmnMyJOA/GVsagAHS6Hds4I+9Rry19YthkUYCkrEE0Dt+AQal6+AlJI68QlsP5aaZzgF4LRYuTJfHPsdbdrx07rVIbPvgog7KTZ+89GUsMVS0GPD/9yymWuatQg7/6WevXlx7iz8qooqJXaLBafVysxdO6kWG0fLylXwqyqPTf+XaTu3Y7dY8aoB+tZvwDt9LuOY222oGxN6LwR3tG7H31s3F3jcvvQ0AD5ZuthQdvfr1Su5p11Hok6jvMD5RuveF/LEd/cx6uHvOJ6cjs1u5fL/9eW2164H4LqnhnL5//qyY80eylVJoGZjYw2cqFgXnyx+nc1LtrF3UxI1L6hOs4sa4/P4SD18HItFQQ1oqAGNOWMWsX/LQT5e9Bq/Jn3BzjV7kFLSoHVdPSGoGKiqyvTv5zJp9HQCfpXuwzqxfsHmsFh0n8fP9B/nMuLN4RFaOr8pDUO+HGgohKgLHACuBa4vhXbPCrrXrkPnGjVZnLQ/z5hbhODCypXZnpoalrmZn4c6XQTofrNfrryaF+bMZPL2bahS0rF6DV7t3SckKiTGbmdw46b8tXUTmtTjAZxWK1ViYjmQmRHiY7dbLNza6oRAj6pprDl0KCw6IRejyBaAq5tdSOsq1RizcR1/bNlEjt/PMbebP7duZvKObTzbrSdJGRlM37UjRJRr+q6dfLhkkeGDIz92i4U3Lu5DvXLleOPiS3lyxlRDVw6Q55PfkHzEcMHUIhSSMtJpVL5C2D7Q30iWHkjin21bsAjBFRc0pU2EOqZSSiZs3sSXq5dzzO2mc41aPNq5C7XiEwq8nnOR7ld1ptuVncjJdOOMcoRpiUTHRxfJpyuEoGnnxjTtfEIhcs7YRWQey8pTXgQ9nnrPhn1sXLiF5l2b0LBNvRKP/Y0bPmLppFV57p99WyIXsS5MJvd85pQNuZQyIIS4D5iKHn74jZTSWPzjHEQRgi8GXsGMXTuYuG0rLpuewt68YmU6fv15geeO37SB5OxsLqlXnyoxsXzQdwDvX9ofGWw3P1O2b+XR6f9iEbrGuE8N0Kdufe5p35F6ieUY8fcfLD+YlDf3DagqX65aTo86dfEFVG78cxxZXp9hJIjLamVohMQnTUo2phxh1q5dpHs8ee1rUuIJBHh1/hysihIWUugJBPh+7SouqFCRjSnJYe2CHoHyzeVD6Rws1HB54yZcWLkKA37+Ho+BMR8QlJCtl1iOfRnpYfsDmkqVmBjDvgCenzOTPzZvyvPb/755I7e2asNjF3ULO/aDJYv4evWKvEXgydu3MnfPbibfcNN5mYwkhCh1DWzQ09KNFkE1TbJz7V6ady26jtHJ7Fy7J1g16UQcuc/t01P0T3rQW+1Wul99UYn7OtcplRR9KeVkKWUjKWV9KeVrpdHmyexJO84Hixfy8txZzN+757TVVjRCEYJL6zfk48sG8tYlfWlbtToOq5UfrriK6Aiv+RKYu3cPry+YS6/vv+andWsA/Q/qZCN+JCuLR6b9iycQINvvwx3wo0rJnL17qBoTS4zdTqLLGXKeBhzIzOTi77/mqt9+Jjk7O0xUC3QjHmN38PP6tTwwZWJeen4uD/47kWdmzWBvRpqhg8SqKGQblHAD3d2R6fViibAAaVMU2p4kV1s3IZGpw28l9qTQxqYVK/LWxX0BuK9DJ5wnVRJyWq0MuaBpmDspl3VHDjNh80ZyAnp8e27i09erV7Hr+LGQYzO8Xj5fuSwkkkeVkpyAn9GrzowG0PlCrSY19NT7k7BYFao3PLUw0vXzNxtmSkpNYnPYcEbr/bpinFSpU4lbXi6aZtL5yDmR2fnP1i08OWMqAakR0DR+27SBi2rUYtSAy0+pgv2p0qJyFebePIKu3442nGH6NS3PJ/7y3Fm0r1adxgYKhJN3bMPI2+5VA3T/7iuGNW3GjF07DSv9aEBmBPdO9dg4Ut05HMvJJiUnm23HUpmxawdv9+lHvwaN2JSSzKzdu8KKZORHSknD8uXZlhpeCBr0Wp35qxblZ1jTCw2zUWvGx7P27vuZv3cPW46m0LNOXRrmc5e0qVqNkf0v54U5MzmYmYHdYuGGC1vx+EVdw9rKZebunYb68xLJnD27qZdYjkOZmXy/djXTdm431MEJaBrLDpS8Es5/kUtu6Mb3L4zF5/blTa4sNgvlqiTS5pILCzm7YBIr63otJ8ej2502rnt6KImV4jmw4zCN29Xnoiva55V++y9y1hvybJ+P/5s5NcRlkOP3syhpH9N27eCyBuHZjqeDlOxssvw+ascnhMyMy0VF8f2Qq/jfpH/wBPz4VdWwyn1ASq4c9yuTr78pxA/rU1XWHj4UMTTPqwYYu3F9gWXhIuEO+PM0ViDoLlFVHvh3EtG26bSvVr3QkMA4h5N3LunH1ePHGPq2/ZoWMRhlV9qxCHt0utWuQ7fadTjmzuGJ6f/y787tCASDGjVmeItWXFSzFkuS9lMtNpbutetgi1B4GvTUfiEEJwdNW4TAabWyNfUow8b9ijeghpSiO5kaceefW+V0Eh0fzccLX+O9EaPYtHgbQhF0uKw1D39xV7EXNk+m86C2WO02IHQdRrEo9B9x8Wkt1HCucdbL2M7cvZOHp04my+D1vm/9BowaMLi0hxdCSnY290+ZyJojh7AIhWi7jbcv6UfPOqF6KaqmsSElmb+2bOKX9WsNjblAF9P65cqrAb3AwpXjfuG42x1xAbCkuKw2JDKiFjvorg/AMAQxlwENG/Fa70sZ8MsPHMg0Fh+KFL1TMSqKnGAFnsYVKpCSnU2O30/HGjV4pFMXasTF49c0BvzyAwczM/LGYRUKGhIBeQ8wl9XKCz16c/VJcsFSSp6fM5NxmzYYLug6rVbm33IH907+h2UHC59tfzFgcJgwmUnR8Hn9KIoo1UzK3Rv28cIVb3P8SBoIgSvayTNjHjrjCTdnC5FkbM/6GbldsUSM8TtZqKq0kVJy05/j2XEsNWhQVNwBP/dO/pu/rh1Og3In6gdaFIWWlauQ4HDy64b1YBB9LoHlB5PwBgI4rFaenDmV5KysAuPUS4ICVI+LxRMIkJQRWfmtoNl0LpO2b2Pazh2GrohcIk0FUnJy8v6/8tCJiNRJ27cxafs2QH8InOyaCcjwvtyBAK/Pn8uQC5qGzMwnbNnEhM2bwoy4IgR2i4UPLu1P+agoVhw6UNBlArrR71W35BEW/3WMpHJPlbrNa/H99k/Yt+UAql+lTvOapzzTPx856+9Ixxo1DSuM6wUQTs0HVxjrg1mZJ7s1fKrKD2uN9YVrJyTwQIdOhbbtDQRYvH9fqRtx0B8hSRkZXFK3fpE+4PKuggXv/VphkeIlJ1d2oCgEpMb+k6JZvl+zKqzMHeiG/O9rh9O3QUNAn9EXhACe7dazzMq9mURGCEHtJjWo16K2acQjcNbfFbvFwuiBVxBtsxFts+G0WnFYLNzUUvehnk6Ss7MMCzOrUoYZlPzc074jAxsaV2MXkBePXlpuLaNZtScQYOH+fYV+8eskJLJ0xD1cUKFisYpQlwUBTaPcSQ+dLL9xRI3dYgm5L8OaXoi1gOsTQnC9qfFico5y1htygA7Va7Dk9rt5rXcfnu3Wk6nDb+HJM6BX3bxS5Yh+1y41IxeHABjRuq2hgVWEwi8b1uKwWqkZH1nasjCTmuBw0q9+Q1pUqhxxtrz9WGqeH9wIm6LwTLeeKELw9aAh1IiNI9pmOysNut1ioXedeiQ4Qw15v/oNDSNjomw26ubTi3n8oq4FPtSi/sNl5kzOfc56H3ku0XZ7sYoklwZVYmK5ptmFjN90okp9bt3Jwtw6SZkZRNlsYZmfPk1l9aFDbD2aQqCAkldxdgdxTgf7I/i4M7weRg64HIC2oz/juCc8w9JusUQMLbQpCl9fPpQLKlTk/cULWXHoAN1q1aZF5Sq8uXCeYXtlictqy8uUzeWYOwdFEQj06BRVSqxCYLNYePuSfiEPpOMed4EPqGHNmp+uoZuYnHbOGUNeVrzYozcXVqrMd2tXk+n1cmn9BtzdtgNxjvAkiPzUSUg0TDNXEMzft4fZe3YV6Hce0qQpDRLL8eycmYb7FUUhJSebNxfMM0zYcQXfGmbv2WXog/ZrGg/9O5F0rzdv/7IDSfy2acNpURiMsdnIKkDOoDDSvR4uH/MTy0bcTZzDyYGMDC4f8xM5fh9eVcUiBBYh6NegEQ91uihMvTHGbkeNsGAbZbNFjFFPyclm0ratpHs8dKlVm7ZVq5kKjCZnHWd9+OG5zJW//WJYzqwoOBQLCAoMS7QpCqqmhSyY5kZr3NyyNfe060ib0Z8ZPlDONJfUrc+M3cY1IIvDba3a8Gz3Xtw3+R/+3bE9TLArxm7nu8FXUj+xHMc8bmrExuVFudzxz5/M37snREvdYbHycs/eDDN4w5q/dw93T/oLTUp8qorLaqN77Tp82n/QWel+Mjn/OWfDD88WjubkMHrlcubs3UV5VxR3tGlP70JC1b4bfCXPz5nJ5O1bUaXMK1JcFLwFJK3kYhT/bbdY+GrQkLyF4DrxCYZ1Rs80644cLpV2Fu3fR0pONnP37jZUXczy+bh63K8oioLdYsEiBE917cG1zVvw3qX9uHPiX6w7chibouBTVW5s0ZKrmoa7VbyBAPdNmRjimsoJ+Jm3dw+Tt281lPktjBm7dvD+kkUcyEinYbnyPH5RNzrWqFn4iSYmhWAa8kLwBgJM2bGNF+bM0jM3NY0dHGPdkcN5pdwiEdA0bIqi/yAMw+ROB7vTjnNRzVr4VJUhFzTls+VLDCUEziQpBvVOS0JSZgbdvv2ywLccDdA0Le+h+cq82VSNiaVHnbqMufIadqcd51BmJhdUqEA5l7GQ1KpDBw3rUeYE/EzYvKnYhvzvrZt5aua0vAfDqsOHuPXvCXw9aEieqJiJSUkxDXkBrD18iJv/+h233x82+3UHAny0dDE3XNiKGIPaln5V5apxv7A/I6PIs/DSQAFqxcXzxcplfLpsCVLqM3cFoxSlU8OqKEgJmiw8zlxyYkHyVDDK8C0MdyDAqBXL6BHMxq2bkEjdhILTuxUhImY6FderIqXkjQVzwxaePYEAby6cx1/X/jc1tE1Kj3Mi/LAs8Ksqt/49gQyvN2IKu0CXSjUy1LP27CI5O/uMGnGAeKeLI9lZfLx0Mdl+PzlBJUWbxUqN2DjDUL0SI3WdwaKaZiMjHm214SzCmATBLN8Ssj75MNf9Ppbv1qzCXciiqyYlDcqVx2KQiBZltRm6YgrCEwhwNF+Wa362HzMWI8vFr6okZ2cVqolj8t/GnJFHYMmB/QWGB4I+03tn4Xw+WrqIrwcNoXW+IgZbjx41rCpUUhQhsAoFUcACqEUoPNmlGx8sWRQ2+/OqAY653fxfl+58t3YVydnZBeqwFIVA0O9/KmQXwd0k0Ot5LjtYeJp9JNyBAEsPJLH2yGF+3bCOP6+5AddJEsRSSr5ctYKRy5eSE/DjsFiCrjELPk3Fpij0bdCQfvUbFqtvh9VKlM1uWHi6SrSxvrqUks9XLGPkimUENA2rIrirbXvubd/JjJoxCcOckUfA4w8UnpWD7jNN83i4JeiCyaVuYmKxSpLViIunW63aXFy3Hj1r10HhxIfjstpIcDoZf/W13NOuI7Xi4g01wG0Whd5163M0gj/aqwa4skkzJlx9falFspyJiBiLorA3Pa3A5KaiouvPpPP75vDaJ9+vXc1HSxeR4fMS0DSy/X4sQnBZw0Y82rkLXwy4gjiHg94/fMOVv/3ClB3bitSnInQjfLJMgMtq5cGOnQ3P+XHdGj5dvoRsvw+vGiDb72fUimV8s3pl8S/a5LzHNOQR6FijBv5CZuT50aRk1u5deb9fWq8BMXZHyIzVIgQuqxWnJfwP+oXuvRg98ArSPB6WHTyAhj7DtikKt7Zszeybbqd5pSo80LEzU4ffwoWVq+RlIypBqdbnuvcixm6nWcXKhmNUpWTC5o0sSUoqXRfLaSagaRzJzi5QpbE4uIML2KDr5qw8dIC1Rw7zydLF4X5sVWVJ0n6GNW3OI9Mn88v6texNT2P14UM8Nu1fPl66qEh93tOuA/e060C0zY5NsZDgdPJ01x4MjlC56bPlS8PG4g4EGLVyWQmu2OR8x3StRCDO4eT5Hr14Zd5sfKqKJiUuq5Voq42jHnfY8aqUZOR7dXZYrUy4+nqenjWNBfv0avfdatXhxZ69Gbl8KX9t3YwSzEJ84qJuXFyvPiOXL2VjcnKe9rpfaiDh+3Wr+XnDWnyqSrXYOFpVqcKI1m1xBwLM2LWDRKeL6y9sSfNKugF/qmt3ho0fY+iff2PhPG6+sFWhfuKzmdJYND2QkcHsPbt4+N/JaEgCmhbR1XQkO4vv1qwKWy9xB/RZ8i2t2kSsXJSLEIL7OnTm7nYdyfJ5iXM4C3RLpbqNferH3G6klKZ7xSQEMyGoELYcTeG3jetJ93rpW78BFqHw4NRJYf5vh8XCvzfcQu2EhLA2cg1qfmW9bJ+PNI+HyjExedsv+fEbdh0vWsy3RQgqR8fw+9XXU9mgjuXFP3zN7rS0Il7lfw+BHnNfFB34huXKE223s+bwobB9sXYHowcOpnxUFL9v3ki2z8+l9RvQpWatMGN7JCuL95csZPbuXUTZbdzUojU3t2xtWOXq0p++Zcex8MIc9RISmXHTbUW/UJPzCjMhqIRcUKEiz/fonfe7JqW+8HbgQF5cuMtq4/rmLQyNOGAojRpttxN9UtiiKIpTPogqJQezMun5/Vf8fe2NNCxfPmR/rfhE05AXgEQvYF0YTquVp7p0570lCw33Z/m8LNi3l6/XrMSvqrr7astGetauy8eXDcybdad7PAwa8yNpbjcBKcEN7y1ewMaUZN679LKwdp/r1ou7Jv0V8pbgtFp5pnvPEl2vyfmN6SMvJooQfDloCK9f3IdederSr35DRg24nKe79Tjltq9q0iys6HBheFWVx6dPCds+ok3bQjW4c7ErFpoa1BI93ynMjAt0SYCkzAx2FBAmOGrFUjyBQJ67J8fvZ87e3czduzvvmDEb1pHl9elGPIg7EGDy9q0kGUgid6tdh28vH0q7qtVJdDppU6UaXw8aQq86ZuELk3DMGXkJsCoKgxs3YXApqzHe0qoNs/fuZkPyEdx+f5HjszemJJPh9RJlszF243rGbVxPTsBPjM1eYGHlXGwWhSYVK7Hj+LES6cKcq9gUS4H1OyV6BSJfQI3ogpEY5w3l+P1M2rY1z/AuPZAUUnc2bwwWCxuSk6kRFy5p3LFGTX4bdm1RLuWU2LRkGz++NI59m5Oo07wWN70wjMbtzXJ35xKmIT+LcFit/Dr0ar5ds4rPVy6LmERyMqqUZHo9PDx1EkuS9hfJeOfHp6rM3L3ztBnx/DU9c/9/+YqVPDZxCtWOp3EwMYF3B17G3+3anpb+I2GQ7xNGSnZ2iRZW9UiiE+Gn9RITWbh/b1jkjarJMi34vGrmep4f/CbeHD1jNnnfUdbO2cDrk5+hRXfjiBqTsw/TtXKWsfzgAd5dvKDIRhx0o/H16pUsTUoqthEHPYU/7TTqj8uT/n/5ipW8PmY8NY6noQA1jqfx+pjxXL7izMZIF2WhUwKJBUSk5NYcPRmHxcIVjS9ge2oqydlZ3NSyddhaiU1RqJ+YSLOKlYo79FLjswe/yTPiuXhzfIx6+NsyGpFJSTBn5GcB3kAgLxTx7UXzi51xqUnJmsOH8Bm8uhdE/pnymeSxiVOIOinqJ8rv57GJU874rLwwNCnJ9Hkj3iu7xcKLPS/mlXmzEehvR6qm0bVWbUZM/BNV0/BrGm2rVuPjfgN5Zf5sDmdlIaWka63avNvnsjILJZRSsm9TkuG+Xev2neHRmJwKpiEvQzYfTeGpGdPYkHIEixD0rd+QrUdTit1OlM1Go/IV2HYslUAR4sMtQpDodJHh85aJT7za8bRibS9rAgW4VqJsNq5pdiEX163H7X9PYHNKClaLhem7QrXXVxw8gFdVmX3T7Rx15+C0WIkNFidRNY1F+/eRnJNN6ypVw4pinC6EEMQkRpN1PDwTOK5C7BkZg0npYBry08zG5COsPnyIKjEx9KhdN6/IQUp2NteMH5On5qdJyb87thVbWs8iBAlOJ49d1I0pO7YXerzLamXsVdcya/cuPipiVmJpczAxgRoGRvtgYgL2oK7JucJxjwcpJW8tmMe21FQCUhIweKPyaxqbUpLZk54Wory46/gxrvt9LNk+Pwg956B3nXpnrHjFVY8M4tc3/sCbky+ZLcrB1Y9dftr7Nik9TB/5acKvqtzxz58MGz+G1+bP4eGpk+n+3ZfsDcZ2/7phXZiiXUBKw2xMe1DT3IjBjZvw5zXDqRAVxU9Dh1EtNjZi+r3TauXHIcNoXqkylaKjQxbjziTvDryMnJN0aHJsNj66fBCBUzTiZ9pJkeh04lNV/tm+tVCfu01RSMk+Mfv9e+tmLv3xW1JycsgJ+Mnx+/GpKv/u3E7P777iQIR6raXJdU8NYdA9l2J32XHFOHG47Ax54DKuemTQae/bpPQwMztPE9+uWck7ixaE+LsF0Lh8BSbfcDP3Tf6HyUUQXYq22hjRpi1709OZunM7fk3DplgQAr4aNIROJ1WYkVKy/Vgq4zdt4Kd1a1AUBX8gQJuq1fl84GDinfrCXYbXS9dvR5dI37s0ODlq5YNBA1jRqyf7DWKqz2Ze6dGbyxo1pvPXowsMZQR91tSichViHQ661KzN+4sX4CtAP6ZaTCxzbxlhmPlZ2rizPaQeOEb56uVwRRcsN2BSdkTK7DQN+Wni0h+/Zcfx8BRrgCEXNKVx+fJ8uHRxgQubLquVF3r05upgPckNyUdYuH8v8Q4nlzVolGeUI5Hj97Pz+DEqRkXhDgSYtG0rfk2lb/2GNK1YiXVHDnPPpL85lJVZ8gstIQ6LBQEENEm9xERqJyQwc9fOUi9+URpUj43lQGb4PYq12Vl7z/36wuW3XxZ6Hy0I1OCSaVH0YqJtNj69bFBeQQwTEzNF/wxTkFLf5O1bqRXXnmibDV9ANaw9CXrm36vz5yCAYc0upHmlynnCWEZIKZm4fSufr1hGak4OHarX4JHOXZi+cwevL5iHKjU0TdfcvqlFa/6va3cW3HoHF30zmiPZWad4xUXHoSg81bUHLpuNVpWrUjU2lvZfjjwrjbiCICUnJ6zQtd1iYeRA3Y8shODV3pdw7+R/Cnwwq/k+56LEpmtScrgMHrIm5x6mj/w0MahR44i+aq+q8uuG9fx57XD6NWiIowBJ2SyfjxfnzjLUzz6ZkSuW8n8zprL5aArJOdlM3rGNQb/+yGsL5uJVAwQ0DQ2JJxDgx3WrWXfkMEII/te+Q6lofReVCypW5KaWrRnWtDkVoqJ4b/GCs7YCjobEp6pYhKB/w8bUTyxHv/oNGT/sOrrUrJ13XK869Rh71bW0qFy51Pz0EmhRpWoptXb6UVWVDQu3sGrmerzu8CIaJqcPc0Z+mrizbQem7twRsZRXutdD9dg4Pu0/CCkl14wfy8pDBwzn5u5AgPcXL+TKJs0i9pfj9zPyJA1rTUpy/H5Dw5JbVFpKybrDh0l0ukguYYFkp8VC/XLlaVqhIhO3by00KWlE63bc/Od4UrKz2ZeeTkBTz8rZeH48QSnjKTfcHJLYo2oaI5cv5bu1usxts4qVaFyhIvvS08gJFqYAYVjXVKAXzTBa4HZZrXStVZsm54gGzraVO3l24Bt4cnwIAZqq8djX/6PH1ReV9dD+E5ySIRdCDANeBJoAHaSU57fjuxjE2O38c92N9PzuKw4buC3a5CsL9/e2LWxIPlxgck5hr9i7jh8zXBSLpAUihGBzSgrfr12dp7deUoa3aMWl9Ruy9MB+nFYrvqAKYCQemjr5lPXEy4IZu3fy8txZvNzrkrxtL8yZyYTNG/EE3yjWJR/Brig8elFX1h4+TILTRffadXho6qQwt0ui08Wwps2YtWc3sQ4HFVwuNqak4LBYuO7CFtzcss0Zvb6S4vP4+L9LXyHzpHj0d279jPqt61Kj4bnzVnGucqoz8g3AUOCLUhjLeYfdYuHDfgO45a/f84ylRQgcViuPdOrChuQjVImJ5e2F8/MMQSSqxxasx1E5OqZY7gmrorDkwP5SSQhae+QwP69fRyBY1zK3DF0kY306jPiZyFL1qSrjNm3g8Yu64bLZOJCRzriN6/GfdD0+TeOvLVuYeP2N7EtPY9K2rZRzOjmYFfpAv6ddB25v044nu566cmZZsmzKatRA+FuFGlCZ+u1sbn/9+jIY1X+LUzLkUsrNgFmtpAA6VK/BhKuv54uVy9iWmkqLSpWxW60M/2M8NouCNxAotISZ02rliS7dIu7fm5bGrxvWkuB0kZqTXWAmYi596zVk1t5dp2zIBXo0Ta47JfdaEp1O0jyeMyYBUJx+TqXCkFVReGvhfP7cugmfqka815uPJvP5imV8tHQxAc34DeWtRfO54oKmlI+KMmwj2+djyo5tHMnOomXlqlxUs1ZeklBA0/h2zSp+WrcGt9/PxfXq80inLlSMji7RdZ0KWcez0Qy+wwG/SnrKuRVOeq5i+sjPABdUqMgHfQcAMG7jel6cOwuvGsBbBBtaOz6Bxy/qSv+GjQ33z92zm/9N/ptAUNMjdzbssFjwa5qhyyTKZqN7nTrM3LMzbF9+qsTEcNztRpMy4sNGgqFPPLsYMrxnEkFw4lFCQ+4JBJiweUOhb1ACeH/JQkP/dy4BTWPgrz8ybfgteen6uWxNPcq148fgVzXcAT8um40mFSry45CrcFptPDx1EjN378pz1/y+eSNz9uxi6vBbiTuprdNNq97N0Qzq2zpjnHQaGBYpZ3IaKDRUQQgxQwixweBncHE6EkLcKYRYIYRYkZJSfD2R84XPVy4vkkKhy2rlg0svY/bNt0c04qqm8ej0KbjzzepVKbFbLNxwYUsal69geJ6UkkvrNyTaZjfcn8tVTZox+YabeaRzF2LtxTMO6lkahaIIUaBxLQxVykKNOOiLmEWJA0p15/DtmlVh2++f8g/pXi85Af2BmOP3szE5ma9Xr2R32nFm7NoV4nMPaBrpXi+/b95QjKspHarUqcTg+y7DGX3iO+KMctC4XX06Djw3/PznOoXOyKWUlxR2TFGQUo4GRoOeEFQabZ6LHHOHF24G/XU/zuEgzeOhQlQ0j3XuErHCei47jh8zjFv2qSoTt20l3RsuTWtVFP7XviMxdjsj+w/iqvFjDNtWhGB4i1ZUio7hrrYdSHN7+GLV8ojHnjzzL44ZtwpRJHdQaXAmF1ktigKFPDQCmsbE7Vt4oGPnvG0HMzPYnx7ukvCoASZs3kTt+ARsihL2RucJBFialMStrc68guQdbw2nVa/mTBo9Ha/bR69ru3DxDd2wFBBaa1J6mK6VM0z76tWZuWtnmNuhYlQ0C269g0BwRl0Uoqy2iNEmxz1uQ3eIIgR3t+0AQJtq1elbrwEz9+wKmaUqQvBij95Uij5R1PnKps0YvWq5cQQMJV9s/LTfQHYcT+WLk95Uykpit6QIwGWzIRBIJFc3acZ369YU6dzC3oxOpkZcvGESmU2xUC8x0eCM048Qgg6XtabDZa3LpP//OqeUBSKEGCKESAI6A5OEEFNLZ1jnL09c1I0omz3Ply3Q3Sgv97oYRVGKbMQBasbHUzchMUwlz1WAGJZfVUOM9nt9+3Nx3XrYLRaibTaibTae796L4S1ahZxXPTYOm2I8tvqJ5fJUHYuDw2KhT/0G3N+hMyNatyUmWIy6XmIizSoVvdiCIgSX1qtfZtltDouFLwcN4eWeF/PmxZfyf126M2ZT0V0cJ8vWVouNo2Z8eOk3p8XK0CZNaVm5CrXi4sMLVVgUbriwFeuTj/Dw1Ek8+O9EVh86WLKLMjmnONWolT+AP0ppLP8JGpQrz8TrbmTUiqWsOnyQOvGJ3NOuA63zxZUXh88HDOa6CWNJ93gBXT3xigsuYGlSErvSjocd77RaceQryhxls/FKrz4czclCIKiTkBiyPxeXzcbVzZozfvPGEHeOy2pleItWvLlwXrEjYAKaRo7fx3uLFzJu0wbsFgtOixW7YqFBufJsTkkpkitESsms3bvKLKlIlZKdx1K5o217pJR0/PrzYhUHWbw/vIjD25f0Y/iE3whIPbM0KrjYOaJ1O4QQ/DhkGI9Nn8Ki/fsQQlAtNo63L+nLO4vm8/e2LXnt/LNtKxfXrcfogVeUKLrsjy2beH/xQg5lZVIzLp4nu3SjX4NGxW7H5PRiulbKgNoJCbx5Sd+w7Sk52SxN2k+M3UGXmrWKNMutGR/PvFvuYGnSfpJzsmlTpRo14+Pp9/N3hsfrZd3cJDhdHMjI4P5/J7IpJRmBPhP8oG9/WlSuEnbeMXcO93XohEVRGLthPRJJtM3OU127M7RJM0auWErOSUUtrEF9kkimWJOSL1eu4PfNG/GqJwoc7ziWit1iwWaxoBbBIFrOoI/diICm8eHSRbSqWpXmFStzPMI6SCSOB8vspebksCvtGNtTU3lt/hwUoSCkHpt/T7sO3NOuY97bV/moKL4dfCWZXi9eVaW8y8WSpP0hRjyXmbt3MX3XTi6tX7yCyuM2rueFubPyHkp709N4ZNoURLAIisnZg2nIzxJGLV/Kx8sWY1UUhBDYFAvfX3FlgSJZuShC0LlmrZBt2T7jSkF2i4VjbjcxdgdXjx/DkeysPD/77rTj3DBhHHNuvj0vtnnHsVQemjqZ7alHAWhasRJ/XnsDiS4X5ZyuvGzSTy8bxM1/jtejOgIBomw2LqxUmVH9L+fPrZt4dd6csBmzBL5cvSLMlx+Qki2pR7m0bgMm7tha6PWfCSNemM/eEwjw07q1fNi3f5EfQLk0q1iJZ2ZNZ8LmjVgVhWyDKk8jly/l+uYtSXS5QrbHOhzk1vL5clXkxOrRq5ZFNOTrk4+wYN8e4h1OBjRsnKeq+d6ShWFvFp5AgHcWzjcN+VmGKZp1FrDsQBKfLl+CV1XJ9vvJ8vk47nFz61+/lzhUrmP1GoYVZixCoWZcPPP27iHD6w1bLM32+3jg34l4An6yfT6uHjeGzSnJ+INx6uuTj3DDhHHE2R0hkgBtqlZj/q138H9dujP0gibUTyjH4awsXp43m+6161I1QmZqQdc3ZWfheu0FYRHijImBSSDd4+GY210sF5PTaqVJxYpMCL6VGBnxXP7dWXAFqIAsIGbdIM5bk5KHp07m2vFj+GDxQl6bP4eLvhnN4v37CGgaydnG2jvnmmb8fwHTkJ8F/LphraFP1RtQWXHwQInafLDjRcTkW1QF3Z/9dNce2CwWDmZmoEYwokuT9jP8j/FM3LYFn6aGzEQ1KfEE/EzbtSPsvASni0rR0UzZsZ31KUfYm57GP9u2cPmYn0LGkR8JEfedaqjgO336MbQAobGT6VW7Dh2r1zDcV5SR1I5PYOH+vTiLsfD7cs/ejNmwrtDqQgFNI9ugCMjetDQemDKRDl+NYlcE/XuAmwx0W/7dsZ3pO3fgDgQISIk7EMAd8PO/yX8jpaRChIzTKJuNsqhjYBIZ07VyFpDp9RkbCqHPkEtCzfh4Jl5/I58tX8qSpP1Ui43lnnYd6VpLl15tWaVqxPKgGrDlaArLDiSF+b1Bf71OSk9n/t49/Ll1M1JKetSuQ6869Xhx7qyQMEJVStx+P45YK9YISn92iwUZbNciBHaLhYrR0ewziKUuKpWioulXvyHPz54Z8ZhG5crTv2FjEpxO+tZvSOWYGHyqSpPPPoxouAtK75+8YytdatUqctik02Lls+XLivTAsigKPWqHFpg4kJnB5WN+Itvvy3uzUiDMhXVhxcpccUGTsDZ/37yBnED456tqklWHDvJwpy68MHtGmOvKHQgwetVy7gqGsZqUPaYhPwsY2Kgxi5P2hWV8BlSNDtWMZ4hFoUZcPG9cfKnhvgsrVaZTjZrM27vH0JD4VQ2rxUKUzRZmzB1WKysOHeDT5UtxBw3Bn1s3n0h/PwkJJGdnR9TpVjWNRzp3ZXHSPmrGxXNzy9ZM2LKJL1YuL7Eq4/0dOrMoaX+BRTz3pqex5WgKIwecKDRst1iol1iOnRFmtwUZ3VS3m/bVqhc5YiXabiO5CAU9XFYrVze7kIbly6MF67raLRY+X7EMd8Afco809EXmVpWrIBDc2KIVAxo1Nv5cIl2K0D+z65q34O2F80j3hmqL+1SVL1aahvxswjTkZwEDGjZmzMb1bEg+kqdhbbNYeK57rzANjtLk8wGDuX/KRKbv2hE2i7RbLPSsU5cNyUfYefxYnt/XYbFQMy6exfv3haWqS4j4yl0xKor9GQHC54v6YuUtrdpwZ9v2pORk8+v6dWxKScahWHCrRV80zM9xj5uK0caugVy8qsqsPbtIzckJEa56qefF3P73hELdHScjgC1HjxY5menCSlVYeSiy6yx39u8JBNibdpzHp//LxGC5vkblypPj9xu+4TitVp7p1pOWhRSlGNqkKUsP7A+bQAhEnsxyJJ99mseDqmlnpJ6oSeGYn8JZgM1i4achw3j7kr4MbNiY65q3YPyw67iueYvT3u87ffoR53CETVxz/H4mbN7I670v5dZWbagSE0PVmFjuaNOOS+s1KLBo8Mm4rFbuatuBjjVqGE6Qm1ashN1iYXtqKpf88C2jVixl7t49aEhsikLVmFjqJiRyS8vWxNsdhou4J6NKjc41ahW6WGxXLKScVFDjopq1GDfsOmKKmXFps1iIstuL7FpZnLQPb0Gl4YIPRQnM2bsnGKYZQAtG9SRlZhjeT7+qUiUm1mBPKP0bNqZX3Xq4rDYU9AeAy2rl0/4D8xLTGpyUrJRLrbh404ifRZjFl03YnprKo9Mms/loaAJObtr5X9fcQP1y5fO2f716Je8umh9xxqoIgVVRsCkWNKlxd7sO3Ne+E7vTjjNk7C94An78moZNUbBZLPw89GpaVq7Cdb+PZdmBpBBDKIButerw3RVXAvri3qvz57Bg314UIfJcO/mxCkGdxETSPR4sQiE5Jzuii8ZltbHyzntwGmTD7ktP4/rff8vTrHEHFR0j/cVUj41j3i0jaDbyYzxFeJOwCMHARhcwbef2IgmpnYyC7srK/5nZLRa61arDl4OuKFIbUkrWHD7E/H17SXDq4Yf5304W7t/LHf/8GeIuclqtfNR3AH2KGZducupEKr5sGnITQC/91nb0yLDFL0UILmvQkE8uG5S37UhWFr1++DqiLzjB6WTaDbeQ6nFTMy6eKJst5Nzv165mzZFDXFC+Ire2akPN+HiklDT89ANDgyuA13r3YWCjC/LS+HN5b9F8Ps/nS1eEQBAeW567CJjf7eGy2nioY2fuaNs+4n1RNY0lB/aTkp3DqBVLI5buS3A6GXvltTQsX56nZkxj3KYNEYtq56d9terc16ETL8+dHdEvXxB14hPI8vvI8ukLnn3rN+CNi/uG3PPtqamMXrmMLalHubBSZe5s2546CUXXZFl2IIl3Fy9gx7FU6iYk8nCnLnmL5iZnFtOQmxTI7rTjDPr1R8MolaoxsSy87c6Qbf9u38bD0yaHzcpdVisv9ujNsGYXFqt/KSUXfPZhRN3zKKsNu8XCuGHXhrwdgG6oPlm2mK2pKew8dixiqr4AKkVHA4JK0dHc064j/RoUPbHlmVnT+W3j+rAFT5uiMP+WO6gUo4uMpXncDP3tF5KzsskJ+HFYLAW+vTStUJGNKcnFFgmzKbq2yrPde3IkK4tYhyPsQbfy0AFu+mN8Xvm93ApVv111LU0rFl3PxuTsIJIhNxc7TQAo74qK6E+uEReezNOvYSO61q7Db5vWM3HrFg5kZFA9Lp77O3aiV516JRqDVVEiGvKcgB93wM9j0//lj2tuYG9aGu6An1px8TwydTLbj6cWmogj0dPht9z7UInGd2eb9vy1dXPIw85ltTK4cZM8Iw56PP3UG25h+q6dbDmaTJ2ERF6eO5sMX3hleU1KNqQkF6n//G8TArBbrNzeui2KEFSNNfaJPz97Zlg4aI7fz8tzZzPmqmuK1K/J2Y9pyE0AiHM46FarDjN3h1YNEsAtLY2lSWPsdm5r1ZbbSkH/+mBWZqGhhhK9rNylP35LUmYGoLuEijOTdZyCPnbt/2/vzqOjqrMEjn9vLUllw0ASggQBgYRFREAURHFDaYILSg/g0qKDTh/tcY72OM40bdv26dGZP3CwbdvW4xy3dgNbXNpGEbBBRRZBgQASIB0ICUgIhCyEylJVv/mjQg2hqqCy1aui7uccziFJ1av7Qrh59Xv3d29mJn+eeTtPfrmS7344QEZyMvdcNJb7xlzMK5u+9Y9d83i4ftBgHho/kWn5BUzL9zeYWr+/gvd3bO/wJiebCHePGsPHu4qpb27mkrw8fjXpavJC/JI9wePzUXw49BCXTQe1K+LZRBO5CqgIsfXaJsLainIKw0wp6ippTmdE09e8xlBac7TD9eW3nGFYx6lqGxspqjxIVmoqw7NzGJ6dw1szZrV5zENLl7CitCRw5bto+1Y+31PKZz+5J7DU8fD4iawoLaGhuYVmX/unJ7kcDi7MzeXxq66J+Dl2EVwOZ8gbwt1Z1qqiT+uHFOBf1w21xdtrzBl7fHSFTFcK4/v1O21vFBuhpxFFyibCBe1YF35+wzomvPwiD376MTP/vJCpb73OD/X1bR6zp+Yoy/5e0mb54kSHycU7trP7yBFe3PgNS3bv5NXpP2bumLHt6jl/giDtfp6IcMeFo3Cd0pbY5XCEfZel4pMmcgX4m2mFS4+nDjDoLs/8aBrDsnNIcThIczqR1td22Gz+PuV2R6eaYPmMibjh08q9pfxxw3qavF7qm5txe1ooPVrNT//6YZvHba08GPL74/Z4eKNoM9MXvcmCdV/z9NrVzH5vEYN69uJfLplAsj30m+GwPWmMCdqiH4lHJ07i+kFDSLbbyUhKIslu5+aCYTwwbny7j6Vily6tKMD/VvuSvnl8s7+izTpust3OrBHtq0DpqF4pqXx020/YUXWIA8fqGZGdw4YD+3l85QqavV58Ph8tEVyNO8WG3W4LKo9MdTopCDOQ+mTulhaeW782qLbbawx/P1pN6dHqwFSfvhk9MCF+BdrFxr7amqAbyI+vXMHnd83lq/Iyth2qpNnrJcnm/yU6sncuI3vnYhfhjaLN2EQC70D+eMPNpCW1b4MS+OvKn516A4cajrGvtpaBmT3DNsNS8UsTuQr4nymFzH5vEdVuN16fDxFhTJ9zeWBcdHtqDM/pzfDWJZBXN39MQ0tLRMspSXY7uWnpfDj7Tma9t5CykxKpQ4SerhSmDDr9JpbXNn/H/DWraQqzocdhs1F3Uu+Ri8/tS9/0DPbWHG1Tuy4YvL7QCX51eRnvzJjFN/sr2FJ5kHMzMpgyaEibyUz3jB7Ll2V7cTkcTD5/cKfXtHunpbeZwarOLlpHrtrw+nx8Xb6P8rpaRvbO5aIQ04Kipep4A5Ne/d+QZYVJdjsDM3syIS+PFq/haKObyYMGc2P+UJIdDmoa3Tz55So+LdmFz/gTbr8ePejpSmH6sOEMy84JOuYXe/fws0/+ctpdlmlOJxv/6Wdtkm5VQwM//+wTNhyowCZC77R0Rvfpw5JdO4Nq2lOdTp646lpmjhjZ4e+LSlxaR64iYrfZuHLAQKvDAPzdH8N1VclJTWPpnXeHfW6mK4WnpxQy//qpzPt8GR/v2snaCv98y9eLNvHIhMu5d2zb/w8vfbchbBIX/DcJf3P15KCZpjlpabw5Yya1jY0cb2lh15HDLC8twW6z4TtlacXr83FtB+vslQpHE7myjDHmtAOBTzTq2ltb0+bzSXY7NxZEVg654cB+Pt61M1CCZ1q7CT69djU3Fgwj96SNPOEm4gj+JlFXDxzExH79Qz7mRFz3fLSY3dVHcLe0BJpKOWy2wE3MJ6+5rk0vE6W6glatqKjyGcMLG9Yz9qXnGfzcAqa8+Sqr95WFfKyI8MzUG0hzJgWqPFKdTvr3OCfiqovPSnbTGKKO2ibCqrI9bT43qf/AkFUxBthXV8s724uY/MYrfLo79BzRFzZ+Q/HhKo63Ntc6sT7fy5XCoxMn8bc59/JjXVJR3UCvyFVUzV/zFX/asimwhFFS7S/pe/PWmYEe2Ce7KLcPq+6+l/d3bKe8rpZxffOYOqQg4prqJLsdW4ipPjYRku12fqiv579Xf8HKvaWBK2cjEtR0y0Bgrf6R5Uu5csD5QVUkHxR/H7KnSk1TIzcNHUZOalpEMSvVXprIVdS4W1p4fcumkJPZf7d+DX+65R9CPi8rNfW0HQpP55bhI3i9aFPQVHuvMVyadx7TF71JtdsdqIrxj5lLJyMpiaONbqqOHw86pkNsrCnfp21cVczQpRUVNYcaGsIOhSg5Ero9bGcNzcrm3y67gmS7nVSHkzSnE5fDwXOFN7K8tISG5uY2pY3NXi81jW6enlLIFf0Hhj1uqPO4ddiIoF4uAgzu2Svmr8ar3cc5eKxehyrHKb0iV1GTm54WNlFEslGno+aOuZgb8oeyqmwPyXY7154/mB7JyfzrZ5+ErFIREYoPV3HrsBEsLdkd1KvEh2HiecE3PR8Ydylf7dsbuNmZ4nQGNuTEqgP1dTy0dAlbD1UiQJ/0DBZMKWRMiGUuFbs0kauocTmc/OPoi3l187dtEqjL4eDhCRO79bVz09OZfUqP9IKsbFx2R8hpPgMyMxl3bh6zL7iQhduL8Pp8ga34zxfeRIozeKJQitPJ4ll38PW+MrZUHqRvRgaFQwpCPjYWeH0+bl/8Lgfq6wL3EMpqa7jrw/f425y5uoEojmgiV1H1yGWXk+ly8dK3Gzja6KYgK5vHr7yG0WcYFNwdZo4YyQsbv6HJ+/+tcJ02GwPO8SdxEeHXV13D7SNH8UXZHtKSkpg6OJ+eKSlhj2kTYdKAgUyKkVr801lTsY9q9/GgG8Fen493t2/lwUsvsygy1V6ayFVUiQj3jR3HfWODNqdFXVZqKu/OvI1ffr6MLZUHsYtw3aAhPHXtdW3q2/OzssjPyjrNkeLTD/Whe8A3eb2U1UbWXEzFBk3kKqENzcpm8aw7aPJ4sIng7MTgiXgzKrdPyI6XqU4n4/P6RT0e1XFataIUkOxwJFQSBxiWncOk/gPb9Ct32uzkpKZFvHNWxQa9IlcqgT0/7SZe2/wdb2/dQpPXy7T8Ah68ZAIuR2zeoFWhaSJXqhvVNjbyzrYi1u8vZ1BmL+66aDQDM3t2+Hhn6k/TXg6bLWbuWaiO00SuVDepPHaMmxa+wbGmZhq9Hr4u38fC7UW8cvMMxvc7L+LjGGN4e1sRv1+/lqrjDfTr0YN5l1/Z7XNUVfzo1Bq5iMwXkWIRKRKRD0Qks4viUiruLVj3NTVud6BO3ePz4fZ4+PcVn7VrB+UbRZv5r69WUXXc352xoq6OR5YvZUVpSXeEreJQZ292LgdGGmNGAbuAeZ0PSamzw8q9pUHNtwAqG45xOEQPl1CMMTwbYuxco8fD/DWruyROFf86lciNMcuMMSd+wtYBWrOkVKs0Z+gZm8YQNNk+nEaPh7qmxpBfi3SQtDr7dWX54Vzg03BfFJGfishGEdlYVVXVhS+rVGyaM2o0KackbKfNxhX9+0c8g9PlcHBOsivk1wack9nZENVZ4oyJXERWiMi2EH+mn/SYxwAP8Fa44xhjXjLGjDPGjMvJCZ6XqNTZZs5FYygcUkCy3U56UhIpDgfDsnN4+vrCiI8hIvz8ssuDfiG4HA4enTipq0NWcarTw5dF5G7gfmCyMSaihT8dvqwSyf66OrZXVZKX0YMLeud26Bjvfb+NZ9atobLhGAPOyeQXl1+p/dATULjhy51K5CIyFVgAXGWMiXi9RBO5Ukq1X7hE3tk18j8AGcByEdksIi928nhKKaXaqVMbgowx+t5OqSipa2riw+LvKT5cxQW9c5k+dDjpSaErY1Ri0Z2dSsWBspoaZrz7No2eFtweDyk7i3l2/Ro+nH0nfTN6WB2esph2P1QqDvxq5XJqmxoDG4Pcnhaq3W5++8VKiyNTsUATuVIxzmcMayvKg4ZA+IxhVdkei6JSsUQTuVIxTgB7mI6HJ+aIqsSmPwVKxTgRoXBIAc5TknaSzcb0ocMtikrFEk3kSsWB31x9LYN79iLV6cTlcJDqdDI0O4d5V1xldWgqBmjVilJxINOVwpI75rB+fwWlR6vJz8pi3Ll5XTpkQsUvTeRKxQkRYUK/85jQjqEUKjHo0opSSsU5TeRKKRXnNJErpVSc00SulFJxThO5UkrFuU4PlujQi4pUAQ3A4ai/eOzIRs8/Uc8/kc8d9Pw7c/4DjDFBI9YsSeQAIrIxVIP0RKHnn7jnn8jnDnr+3XH+urSilFJxThO5UkrFOSsT+UsWvnYs0PNPXIl87qDn3+Xnb9kauVJKqa6hSytKKRXnNJErpVScszSRi8h/ikiRiGwWkWUi0tfKeKJNROaLSHHr9+ADEcm0OqZoEZGZIrJdRHwikjClaCIyVUR2ikiJiPzC6niiSUReEZFDIrLN6liiTUTOE5GVIrKj9ef+oa48vtVX5PONMaOMMaOBvwK/tjieaFsOjDTGjAJ2AfMsjieatgEzgC+tDiRaRMQOPA8UAiOA20VkhLVRRdVrwFSrg7CIB3jEGDMcmAD8c1f+21uayI0xdSd9mAYk1J1XY8wyY4yn9cN1QD8r44kmY8wOY8xOq+OIskuBEmNMqTGmGVgITLc4pqgxxnwJVFsdhxWMMT8YY75r/Xs9sAPI66rjWz5YQkSeAuYAtcA1FodjpbnAIquDUN0qDyg/6eMKYLxFsSiLiMhAYAywvquO2e2JXERWAH1CfOkxY8xHxpjHgMdEZB7wIPBEd8cUTWc6/9bHPIb/rddb0Yytu0Vy7gkm1Fy2hHoXmuhEJB1YDDx8yopEp3R7IjfGXBfhQ98GlnCWJfIznb+I3A3cCEw2Z1lRfzv+7RNFBXDynLZ+wAGLYlFRJiJO/En8LWPM+115bKurVvJP+vBmoNiqWKwgIlOB/wBuNsYctzoe1e02APkicr6IJAG3AX+xOCYVBeKfkv0ysMMYs6DLj2/lRaCILAaGAj6gDLjfGLPfsoCiTERKgGTgSOun1hlj7rcwpKgRkVuB54AcoAbYbIz5kaVBRYGITAN+B9iBV4wxT1kbUfSIyDvA1fjbuFYCTxhjXrY0qCgRkSuAr4Ct+PMdwC+NMZ90yfHPsnfzSimVcKyuI1dKKdVJmsiVUirOaSJXSqk4p4lcKaXinCZypZSKc5rIlVIqzmkiV0qpOPd/Gr9sBkRsfcoAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot dendrogram\n", "# Read https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.cluster.hierarchy.dendrogram.html\n", "# to learn other options\n", "\n", "plt.title('Kmeans Clustering')\n", "\n", "plt.scatter(dataTW[:,0],dataTW[:,1],c=labels)\n", "plt.scatter(centroid[:,0],centroid[:,1],c='r')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "celltoolbar": "Slideshow", "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.8.5" }, "latex_envs": { "LaTeX_envs_menu_present": true, "bibliofile": "biblio.bib", "cite_by": "apalike", "current_citInitial": 1, "eqLabelWithNumbers": true, "eqNumInitial": 1, "hotkeys": { "equation": "Ctrl-E", "itemize": "Ctrl-I" }, "labels_anchors": false, "latex_user_defs": false, "report_style_numbering": false, "user_envs_cfg": false }, "toc": { "colors": { "hover_highlight": "#DAA520", "running_highlight": "#FF0000", "selected_highlight": "#FFD700" }, "moveMenuLeft": true, "nav_menu": { "height": "12px", "width": "252px" }, "navigate_menu": true, "number_sections": true, "sideBar": true, "threshold": 4, "toc_cell": false, "toc_section_display": "block", "toc_window_display": false } }, "nbformat": 4, "nbformat_minor": 2 }