openai-cookbook/examples/Clustering_for_transaction_classification.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Clustering for Transaction Classification\n",
    "\n",
    "This notebook covers use cases where your data is unlabelled but has features that can be used to cluster them into meaningful categories. The challenge with clustering is making the features that make those clusters stand out human-readable, and that is where we'll look to use GPT-3 to generate meaningful cluster descriptions for us. We can then use these to apply labels to a previously unlabelled dataset.\n",
    "\n",
    "To feed the model we use embeddings created using the approach displayed in the notebook [Multiclass classification for transactions Notebook](Multiclass_classification_for_transactions.ipynb), applied to the full 359 transactions in the dataset to give us a bigger pool for learning"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Setup"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# imports\n",
    "import openai\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "from sklearn.cluster import KMeans\n",
    "from sklearn.manifold import TSNE\n",
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "import os\n",
    "\n",
    "openai.api_key = os.getenv(\"OPENAI_API_KEY\")\n",
    "COMPLETIONS_MODEL = \"text-davinci-003\"\n",
    "\n",
    "# This path leads to a file with data and precomputed embeddings\n",
    "embedding_path = \"data/library_transactions_with_embeddings_359.csv\"\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Clustering\n",
    "\n",
    "We'll reuse the approach from the [Clustering Notebook](Clustering.ipynb), using K-Means to cluster our dataset using the feature embeddings we created previously. We'll then use the Completions endpoint to generate cluster descriptions for us and judge their effectiveness"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "df = pd.read_csv(embedding_path)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(359, 1536)"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "embedding_df = pd.read_csv(embedding_path)\n",
    "embedding_df[\"embedding\"] = embedding_df.embedding.apply(eval).apply(np.array)\n",
    "matrix = np.vstack(embedding_df.embedding.values)\n",
    "matrix.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/ted/.virtualenvs/openai/lib/python3.9/site-packages/sklearn/cluster/_kmeans.py:870: FutureWarning: The default value of `n_init` will change from 10 to 'auto' in 1.4. Set the value of `n_init` explicitly to suppress the warning\n",
      "  warnings.warn(\n"
     ]
    }
   ],
   "source": [
    "n_clusters = 5\n",
    "\n",
    "kmeans = KMeans(n_clusters=n_clusters, init=\"k-means++\", random_state=42)\n",
    "kmeans.fit(matrix)\n",
    "labels = kmeans.labels_\n",
    "embedding_df[\"Cluster\"] = labels"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 1.0, 'Clusters identified visualized in language 2d using t-SNE')"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABazElEQVR4nO29eXxcV3n//z6zaWY00kgaLZY03uIttrM4sSeOSUKAAFkwNQGXL22BJA1JaQuFtIFAaUnoThuaGtp+afozhLZQFoNr8BdnbVZInHEWx3Ecx1tsj2RrH0mj2aSZ8/vjuSONZMmSrW0knbdees3MvXfufe6dmc957nOe8xyltcZgMBgMsx/bdBtgMBgMhqnBCL7BYDDMEYzgGwwGwxzBCL7BYDDMEYzgGwwGwxzBCL7BYDDMEWaE4Cul7lNK/dd023EuKKV2KaVuGWHdIqWUVko5JunYMaXUBdZzj1LqF0qpTqXUT5RSv6OUevQ893urUuq58do0GQy9pme7/uM4xojfQ6XUNUqpg+e53/O+rnOBybg+SqkF1nfSPpH7LXQKRvCVUr+tlNpjfQinrB/s1RO4/0kV2aForW/UWn9vso+jlHpKKfWpIcf2aa2PWi83AzVAQGv9m1rr72ut3z/Zdg1liE1Tcbwpuf55x3tWa71iqo43k1FKXamUekwp1a6UarEckdqptEFrfcL6TmYmet9KqYeUUn81yjZlSqnvKKVOK6W6lVJvKaW+lLdeK6X2KaVsecv+Sin1kPU8p2exIf//52zHLQjBV0r9MfBPwN8g4rQA+Fdg0zSaNYipaigmgYXAW1rrvuk2xGCwKAceBBYh389u4LvTadA08ADgA1YCfuA3gMNDtqkDPjbKfsqshiv3/6Ozbq21ntZ/62RjwG+eZZv7gP+ynr8LiAxZ/zbwXuv5FcAeoAtoAv7RWn4C0NaxYsAGa/nvAgeADuARYGHefjXwh8Ah4BigrA+q2dr/PuCiEWx+CviU9dwO3A+0AketfWrAkXcNtgKngAbgrwC7te5W4Dnr/R2WHTda6/4ayABJ65z+Oc/upcDXgDTQa62/Pbe/PDsvBB4D2oGDwEfz1gWAn1vn+iLwl/nvHXK+u4DPDFm2F/hwvk3W85uAN5AfegNwd/65DtlH/vs+ALxi2XMSuC9vu0VDrmn+9d+b97nHrO3eZa27Evg1ELW2e1fePhcDT1t2Pgb8M9b3cJjzfxd530vkO3k38BrQCfwIcI/w3qGfyRbr/LqAl4BrhvwWfgz8h2XXfmBd3vrLrWvUDfzEOu5fjff6Wus/CRwH2oA/Z/DvzgZ8CThirf8xUDFGDbgc6D7P792g634OejDc9+UvgV9Z1+5RoHIs5z7k2Hciv7c08l37xQh2vw586CzXRAP3INqTs/GvgIeGs3+s/4Ug+DcAfWcznHMT/OeBT1jPfcCVI10g5A7iMNLKOoA/A3495KI/BlQAHuB65AdYhoj/SqB2BJufYkBwPg28Ccy39vXkkC/bduDfgGKg2vqS/17ej7QXuANpOH4faATU0OOM8CPuv3ZDf/TW8U4Ct1nnfxnSKK2y1v8Q+eEWAxch4jzSD++TwK/yXq9CRLRoGJtOYYkY4u1dPtS2Ec7lXcDFiLhcgvyAP3SWH/CnhrHzTuuzKAXqkR/wTdY+32e9rsr7Lv0jUAS8ExGCcxH8FxEvrQJxKj49wnsHnTfwcUT0HMCfAKexGgvr80xaNtuBvwVesNa5EFH6HOAEPowIz1gF/2zXdxUiYFdbx7kf+V7mfnefA14Agtb1+jfgv8eoAZ/PncN5fO8GXffz1QPr+3IEWI781p8C/m4s5z6MTQ/lrvlZzvn/Qxrr24Blw6zXwDJEb3I6Mm7BL4SQTgBo1RMXcugFliqlKrXWMa31C2fZ9tPA32qtD1jH/xtgjVJqYd42f6u1btdaJ6x9lyBesbLed2oMNn0U+Cet9UmtdTvyIwVAKVWD/Hg/r7Xu0Vo3I3cR+bdyx7XW/64l3vg9oBYJfY2XjcDbWuvvaq37tNavAD8FftPqzPoI8FXLrtetY4/EdgZfu98Bfqa1Tg2zbS+wSilVqrXu0Fq/PBZjtdZPaa33aa2zWuvXgP8Grh3bqYLVJ/RXwG9orbsQYf2l1vqX1j4fQ7zBm5RSC4AQ8Oda65TW+hngF2M9lsU3tdaN1mf+C2DNWN6ktf4vrXWb9Zl8AxHQ/P6B5yybM8B/Apday69EGolvaq17tdY/QxqdMTHK9d2MeKvPaa3TwFcRwcnxaeArWuuI9ZnfB2weLRSqlLrE2tcXrNfn+r0bjXPRg+9qrd+yfus/ZuDzGu3cz4fPAt8HPgO8oZQ6rJS6ccg2Grmb+HOllGuE/bQqpaJ5/yvPdtBCEPw2oHICY+S3I630m0qpsFJq41m2XQhsyV0sJKyhEM8vx8ncE631/yK39f8CNCulHlRKlY7Bprr8/SBeWL4NTuBUnh3/hnj6OU7n2RC3nvrGcNzRWAisz//CIEI9D6hCxGMkuwehte4G/h8DDdVvIV/o4fgI0sgdV0o9rZTaMBZjlVLrlVJPWh19nYjIVI7xvfORH/EtWuu3rMULkcYt//yvRhrUOqBDa92Tt5sRz38ETuc9jzPGz0wpdbdS6oCVWRVFQn755zl0v27r91MHNGjLBbTI//xGO+7Zru+g77D1PWzLe/tCYHvedTyAhBtHdEyUUkuRUODntNbPWovP6Xs3Bs5FD0b6vEY797NiZcblOlV3WftIaK3/Rmu9FnF6fwz8RClVkf9erfUvgQjweyPsvlJrXZb3f+BsthSC4D8PpIAPjXH7HsCbe2F5BFW511rrQ1rr30IE8+vANqVUMcO3yCeR0En+BfNorX+dt82g92mtv2l9SKuQL9IXxmDzKSSck2PBEBtSDP7gSrXWq8ew3zPsO0dOAk8POX+f1vr3gRYk1DaS3cPx38BvWQLuRkJXZxqsdVhrvQn5jP4H+bLDmZ/tvCFv/QES252vtfYD30Ya6LOilPJYx/knrfWuvFUngf8ccv7FWuu/Qz6zcuu7k2O08x83SqlrgC8id4XlWusypA9g1PNEbK5XSuVvm//5jef6nkLCNbn3ehChynES6VvKv5ZurXXDCOe5EHgc+Eut9X/mrTrX79356sG5MNq5D2WoZnxfD3SqDvXise42/wYJYS0eZn9fAf6UvPM8X6Zd8LXWncgt0r8opT6klPIqpZxKqRuVUn8/zFveQjyaDyilnEjcvSi3Uin1caVUldY6i8SQAbLIFykL5OeCfxv4slJqtfVev1LqN0eyVSkVsrwgJ/JFS1r7HI0fA3+klAoqpcqRzq3c+Z9COoi+oZQqVUrZlFJLlFJjDVU0DTmnc2EnsFwp9Qnrmjutc1xphQt+BtxnfSargFtG2d8vEU/vL4AfWZ/BIJRSLsvj8Wute5HOtNx2e4HVSqk1Sik3EhbIpwRo11onlVJXAL89xvP8DvCm1nro9+m/gA8qpa5XStmVUm6l1LuUUkGt9XEkvPM1y+argQ+O8XjjoQQRvBbAoZT6KtLfMBaeR7zqzyilHEqpTUinZY7xXN9tyLV6hxVeuI/BjdC3gb/OhfSUUlXW8c9AKVUP/C+SZPDt/HXn8b07Xz04F0Y796GM+ptUSv259VtzWZ/F5yz7zhjLobV+CunkHe33NyrTLvgAVpzyj5EPqwXxFj6DeGVDt+0E/gDp9GhAhDeSt8kNwH6lVAzJdviYdfsUR7JafmXddl6ptd6OtPo/VEp1IRf1jBY4j1Lg35FsmVyP/T+M4RT/HckA2gu8jHyh8/kk0hn0hrXvbUhYYSxsQWKlHUqpb47xPUB/GOb9SBimEbml/ToDP5jPILe1p5GOqO+Osr8Ucm7vRbzFkfgE8LZ1zT+NhJGwQi1/gXh+h5DspHz+APgLpVQ34iT8mLHxMeBmNThf+Rqt9Umk4/5PGfjefYGB38VvA+uRUN+9SGbMZPMI8DAiZMcRp2JMYRkrvvxhJIwRRfoodiJ3kOO6vlrr/Ujc+YeIxxtDstVyfTRbkLuDR633v4Bcu+H4FCKI9+V/Jnnrx/y9O189GGl/IxxjtHMfylakjyqqlPqfkXaLnFcr8tt7H/ABrXVshO3/DOn8H0p0yPf6j892LrlMD4PBMAtRSu0Gvq21PmtjfR779SGNyjKt9bGJ3HehM5PPvSA8fIPBMDEopa5VSs2zQjq3IOmVD0/Qvj9ohVmKkdTEfUgK5Kxntpy7EXyDYXaxAgkdRpEc/s16bKnDY2ETEn5oRHLEP6bnTohgVpy7CekYDAbDHMF4+AaDwTBHKKiCYJWVlXrRokXTbYbBYDDMKF566aVWrXXVaNsVlOAvWrSIPXv2TLcZBoPBMKNQSo1pNLIJ6RgMBsMcwQi+wWAwzBGM4BsMBsMcwQi+wWAwzBGM4BsMBsMcoaCydAyG2cDON3ey9ZWtNHY3UldSx+2X3c7GC89Wht1gmBqM4BsmnUgEwmFoaYGqKgiFIBgc/X0zkZ1v7uTep+8FIJaK8Xbn2/w68mvuab2HP776rIUMDYZJx4R0DJNKJAI7dkA8DjU18rhjhyyfjWx9ZSsAzT3NZHSGCncFDuXgGy98g3AkPM3WGeY6xsM3TDiRCOzaBc88A/v2QXExXHYZrF0LldaEeeHw7PTyG7sbiaViFNmLcDlkGtKSohLak+1sO7CNUDA0zRYa5jJG8A2DCIdh2zZoaID6eti8WUIw+WgNaoT5fiIR+O534bXXZB+trdDcDG1t0NgIH/wgVFRAU9Pkn8t0UFdSx9udb1PhHpirIpVJ4S/y09A17Gx/BsOUYUI6c5RwGO65Bz7+cXkMh+X//vuhuxsWLJDH+++X5Tnuuw/uuktEfzhefBF+9jN46SWIxcDrBY8H0mk4dGhgedWoVT9mJrdfdjsAXakustksid4Eib4ES8qXUF9aP83WGeY6xsOfg+SEPRAQYY9G5bXLJcsC1vTMgYCI/t/8Ddx4o4RjIhHYKmFqHnhgsKevNfzLv8Crr4LPB11dsl5rcLtl2YEDsHIlXDvWGXtnGBsv3Mg9rffwjRe+QXuyHX+RnzU1a/C7/WxeuXm6zTPMcYzgz0G2bRte2J98Ei6+WLJpaq0ZdZubxSOvqZHHNWvg9tthyxZZnxN9rcXz/9//hbIySCTAboeiIujtlfdqDX19smw288dX/zHXLLqGbQe20dDVQH1pPZtXbjbxe8O0YwR/DtLQIJ59js5OOH4cenrgxAnpZG1pkcdMBubNA5sNSktl+5tuEm89X/TvukteL18uAg/y2NcnDYJSEt5597slxLNjB2zaNDs7bgFCwZAReEPBYQR/DlJfL2GcnId/5Ih48vPmQTYLqZR46I2NUF0N73gHPP+8NAwgoZqrr4YPfEBEPif8H/iAxObfflvEvaFBBF9raSAqKmDdOmk42trgW9+CJUtmf26+wVAomE7bOcjmzSK4bW3Q0QF790J7uwhyRYWEd9rbJWf+yiulMUil5HHHDnj0Udi+/czsnQ9+UBqNkhJYulQahYsvFrGvr5ftKyslc2f/fnmcC7n5BkOhYAR/DhIKwd13i+f94osSrlm0SIQ+FoNVq0So580TL99mkwbguefkPbW1cPo0PPjg4P3+4hfyHr9f7gJyYaDKStn/2rWy3aFDEt+vqRnYpqxscDaQwWCYeExIZ44SColHffXV8MIL8NRTkqXjcIgn7/NJB20yKYL89NMSky8ulnBP7g5hyRL4whck+2bLFukH+NCHRPTfeEPe88EPyv5cLgkZnT4NTicsWzZgj89XeLn5kc4Iuw7v4uXGl9FKs652HTcsvYGg38SeDDMTI/hzmJYW8bIvv1w8/WxWYu59fbBiBWzYAHv2QF2d5NGXlUnn6+nT4sEvWiQx/upquPNO2eeWLdIofPCDcN11A7H5XD2dpiaJ2ddbKem5vgGnEy66aLquxJlEOiM89OpDHO44TLm7HIXi2RPP0hRr4pY1txjRN8xIjODPYaqqJIRTWSnCfOKExNP9fnjPe8Qjv/xyEeriYhH9SERCP5WVsqy4WN6rlGTrgIj+0qXSCOTy9IPBgU7ZSAS+9z0ZhBUIyHE6OqQhiUQKo/M23BimNdFKwBPA6/QCoJSiJdFCuDFsBN8wIzGCP4cJhaSzFKTWTSwG5eXwzneKCEejkjrZ3i6dtr/6leTQu90i+PE4fO5zAwKdL/q5u4HhCAblrqClBU6elH+QTlyXC774xck867HR0tNCui9NhWegRILb4aYj0UFLT8s0WmYwnD9G8OcwwaAIejgssferrhKRzmQkrfLaa2Wb5cvl8ZJL4Ne/FmEuKZGGYeOQMu850R9J7HNoLft86SXJDCoulgbnxz+WXP2hGUBTTVVxFS6Hi0Rfot/DT/YlcTlcVBXP0roQhlmPEfw5Tn6oZSRydwLz5sEnPynCnPP+h2M0sQcJJ/3859JwlJTIsqIiKC5N8Df/dpwb1TNUFVcRqgtNS/gkVBdif/N+DnccJquzKBRtiTaWVywnVGcGVBlmJiYt0zAquTsBr1c6Xb3e8Y+SDYUkG8jhEG8/lYLOWAIdOEjzKRc1vhrivXF2HNxBpHPqE/SD/iC3rrmVq+ZfRbI3SbwvzjULrjEdtoYZjfHwDWNiLHcC57q/tWvh2LGBEJK3toFU0oO/rg+bslFaJLUcpquTNOgPcsfaO2DtlB/aYJgUjOAbpoxIZ4RwY5iWnhaqiqv4wG9u4D//bR6BgHTyPncwSTbh58aPHut/j8/loylWYAn6BsMMxYR0DFNCpDPCjoM7iPfG+8M1x30/5RO/d5qSEkkJDZQ7eN/v7GP5JZ3974ulY6aT1GCYIIyHb5gSwo1hytxl/WGa3GNvyfN8/es3AxDp9LHj4DG6UmX4XD5i6RjRZJRrF87S4vkGwxQzbsFXSs0H/gOoATTwoNZ6i1KqAvgRsAh4G/io1rpjvMczzExaelqw2+zsb9lPZ7ITv9vPkvIl9KR7+rcJ+oNsWrGJcGOYplgTVcVVXLvwWtNJWmCEI2FT63+GMhEhnT7gT7TWq4ArgT9USq0CvgQ8obVeBjxhvTbMUWzKxjPHnyHVl6LMXUaqL8Uzx5/BpgZ/BYP+IDevvJk7193JzStvNmJfYIQjYe5//n66U90s8C+gO9XN/c/fTzhiKt/NBMYt+FrrU1rrl63n3cABoB7YBHzP2ux7wIfGeyzDzEWjUVaCvsJ6VArNCJPjGgqSbQe2EfAECHgD2G12At4AAU+AbQe2TbdphjEwoZ22SqlFwGXAbqBGa33KWnUaCfkM9547lVJ7lFJ7WlrMkPXZitaaa+ZfQ5G9iI5kB0X2Iq6Zfw16pNnQDQVJQ1cDZe6yQcvK3GU0dDVMj0GGc2LCOm2VUj7gp8DntdZdKm+4pdZaK6WG/WVrrR8EHgRYt26d+fXPUqqKq4j3xtkwf0P/sq5UV3/ZAsPZGZrSOl0jkOtL64kmowS8gf5l0WSU+tL6KbfFcO5MiOArpZyI2H9fa/0za3GTUqpWa31KKVULNE/EsQwzk1BdiB0HpVLbXMjAGa1j81wEPJfSWuYuo8ZXQywdY8fBHWxasWnKRX/zys3c//z9gHj20WSUtkQbt625bcT3FEpjNdlEwhEObDtAV0MXpfWlrNy8kmBo4Dw7I500hhvpaemhuKqYulAd/qB/Sm1U472lVuLKfw9o11p/Pm/5PwBtWuu/U0p9CajQWp+1DuK6dev0nj17xmWPoXCZKz/8XMdmwBMYJIp3b7ibUDBEpDPCN1/4Jm+0vUF3qpuSohJWBVbxR1f+0bDXY/uB7cR746QzaQ61H6Khq4G2eBvl7nI+vOrDU34dzyVLJ7+xym/op6Oxmkwi4QjP3/88noAHd5mbZDRJoi3Bhrs3EAwF6Yx0cnDHQdxlblw+F+lYmmQ0yYpNKyZE9JVSL2mt14223UR4+FcBnwD2KaVetZb9KfB3wI+VUrcDx4GPTsCxDDOYoD84q37kI5HfsQn0P247sI1QMMQP9v2Ap088TYWngmpvNbHeGE+feJrK4kq+ePWZPlEupTXcECarszTFmlAojnce52TnSRq7G6dUQEPBUL/A5xrxV/a8MqgRj3RGePjww/zkjZ/Qm+llVdUq1tatpdJbCUxfuYzJ4sC2A3gCHrwBCVHmHg9sO0AwFKQx3Ii7zE0mnaFhdwOdkU6S7UlOvnCSRdcuAi39XJPt+Y9b8LXWzwEj1Ue8brz7NxhmGg1dDSzwLxi0rMxdxonOEwA8cewJyj3lgwahaTRPHHtiWMGvKq7i6eNPU+wq5kTnCdxON2hwOpycjp1mdfXqaRHQnPeeyWY4FTvFC5EX2HV4F7+x7DfY27SXt9rfIpPN4HP5ONB6gGgiynVLrqPCUzGjy2UMF7rpaujCv2CwSLvL3HSekFHjPS092Ow2GsIN6Kymp6kHpRStB1vpS/ZRVFzE/Gvm0xvv5eCOgxPm+Q/FlFYwGCaYXMdmPvkdm4neBC6ba9B6l81Fojcx7P5CdSGaY81kdZaedA86q0llUiwqXURnqhOfyzctk7KEG8NkshneaHmDdCZNbUktDuVgy+4tHI0epdJbSbmnHJvNRomrhK7eLg61H5rR5TJyoZtUdwr/Aj+p7hTP3/88WmmS0eSgbZPRJKX10qgXVxXTvK8ZV7GLZEcSh8eBw+1AaUU2lcUT8NBxpIOi0iLcZW4aw42TYr8RfINhgtm8cjNtiTba4m1kshna4m2c6DxBwBPgwT0P4nV4OR07TSqTQmsR745kB6urVw+7v6A/yLsXvxud+1OaZRXLcDqc+Iv80yagLT0tnIqdothVjNfpxaZsVHgqiPfFaehswO1wU+erI9mbRKPpzfRyuvs00WR0xs4pkB+6sdlteANePAEPNmUj0ZYg3hYnm8nS09pDoi3Bys0rAagL1RFrjqGzmnRPGp3V9CZ7cfqcZPoyODwOUp0pAFw+Fz0tPWcz47wxgm8wTDChYIi7N9xNSVEJJzpPoNGsqlxFfWk9Nb4aNszfQKI3QUe8g+5UN4m+BP4iP7+75ncH7SfSGWH7ge08uOdBFIpKTyU3Lr2R2uJa+rJ9xNIx5vnmTYuARjojHOk4wtPHn+Z49DidSQldJPoSVBdXk9Zpkn1JSt2lrAisQGtNOpOmqrhqRnfYdjV04S5zD1rmLnOjs5oNd2+gqKSIt3a+xanwKa78kyv7s3T8QT+L370YNKRjaRp2N9Cyv4WuE110HO2g6bWm/oGJ6Via4qriSbHfFE8zGCaB/I7NXJZNLmZ/ed3lALx6+lWKXcVnZLpEOiPsOryLJ489SbWvmourL8btcKNQeJweLii/gPZkO+Xucub7509Jlk5+hpVSiuZYM/Ul9ficPqLJKLF0DLfDzdvRt+nL9JEly5v2N6krqeN0z2laE60s9i9m88rNM1bsAUrrS0lGk9iddrobu+mN96L7NGVLygiGgtSvqyeTyrB7y25++Bs/xO6246vxsfb317L0xqV0n+4mdipGvDmOu9yNcil0WhM7HcNZ4qTjWAc2u42F1y6cFPuN4BsMk0xLTws1vsEDzdfUrqG2pJY71905aHmuI/RQ+yFqfDXYlI1wQ5j1wfUsKl+E1+mVSVkmgXxRtymbhI+0xqZsnI6dZnH5Ymp8NTz59pN0JDp416J3ccPSG3jkyCNEuiM0djVS7i6npKiEecXzOB49TkN3AwFPgMvnXc6KwApebHyR2pLaGSX6g67L1Tbi343ja/bhrnCjM5pkZxJnsZPOSCf+oJ/KSytx+VzEm+IUlRXRmezk8T95nPifxjnx9Ak6Dnfgq/fhm+fD7XeDDZRD4XK76G7o5orPXlG4WToGw2xkLGMGxjquoKq4ilg61u/hg9T570h2cM9j9wzKZ490Ryhzl9Gb6aXcXd5/m3+o/RDr69dPWnZLfr683WZn16FdtCXaCJYE6Uh14HF4qC+px1ZkozfTS8AT4FD7ITYEZeT0N3d/E5ddpqYsKSpBKUXAG6CkqITfW/d7/cfpSnXNqJTMMwa9LYzxxto3WLxnMapD4Ql4uPB9F1I6r5TGcCP+oJ+Xv/0ynhoP2CAVTeEud2Nz23jua8+R7k6zbOMySueX4in3oGwKrTXJjiRL3r+EWFNsUgdjGcE3GIYwlpGt5zL6NVQX4qFXH6I10Uq6L43L4SKRTtDY3UiZu4xEb4IXIi/w3InnuLT6Uq5fdj3+Ij+JvgRepxe3w90fNpmsztn8+QqejzxPR7IDt8NNoi9Bqi9FT7qHlxpfGmRbLm7fmmjFaXNSW17LfP98AFJ9KZrSZzZOM20Gs+HmcSgrLoMPwTsq3tG/nc5qYk0xAHqae7C77BRXF2Oz20h2DGTv1K6tZfX/WU3rG5KO6fQ66Uv0UeQvmtTYfQ4j+AbDEIb7kbfF2/jWi99iSfkSqoqraI23Djuhy0jeqxoyVOVA6wEqPZW0JlopshdRU1xDe6KdX0V+RY2vhq5UF681v0ZZURk1xTUUFxVPaimK/LDTkbYj+Fw+iuxF9PT2UOYuozvdzeHoYa7nepZVLOOXh39JojfBrkO7eLPtTZx2J267dGbG03HaEm2cip1C+RSt8db+AVczLSVzuHBcRVkFzT3NUDGwLF+si6uL6TohnbvFNcWDBP/Cj1xIqjNF8bximvc305voRWc0VQurSEaTkxa7z2EE32DII9IZ4ZEjj2DDJrX6FXQluzgdO025p5x3zH8HsXSMbfu3oVA0x5uJpWP4XD4WlS1iRWAFN6+8edA+w41hFpUv4pJ5lwwsi4RpjDUyzzePIkcRAOWecpp7mnn06KNcWnMpl1RdwuGOw7xy+hU+uvqjw9495Dp4X258Ga0062rXccPSG845ZJIfdtJWncPOVCddqS4SfQmaY8343X6yOkt3uptkb5JiZzFHo0dpibWgULQn2wGIpqKkM2k8Dg+1JbU89fZTvHPhO/vvVGZS/aThwnGuRS58L/tIdaUGlUnIifXa31/L43/yOInOBJlkZtD+Ir+KcNP/vYlTe07Rm+gl2Z7EXe7GP98/JbV1jOAbDBa5MI3b4SaWinEsegytNR6Hh6zOEk1GaU+005Ho4Gj0KJlsht5sLwpFvDeOXdlp7mkmHAn318wJN4b56Rs/JegPsqJiBZXF4unW+Gp4vfn1QSNye9I9eBwegqVBytxldKY6ubz2cub55jHfP39Ysf/So1/imZPPkM6kcTvcPPf2c/zirV+wcflGblx645iFP7+43dLypbwYeZFoOspC/0K8Ti8xZwybsnGg5QAdyQ7W1a3jVOwUQX+QZeXL2Nu0V/olUh30pHooLirmIxd+hMXli9nbvJdwY5h1tetw2V388tAvZ0wtpeGK/vV6e9l400bUAUWsKUZxVTELr13YL9Zrb11LvC3eH7N3+V2s/OhKooeiHPp/h3jhGy9w/QPXk19ReKowgm8wWORCOZdWX8oPX/8hLrsLl93FsegxaotrWehfyKH2Qzx/8nkSvQma48147B7KPeXYbXZaE61cs+Aath3YRm1JbX+Mf75/Pp3JTnY37GZ9/XoqiytZM28N+1v2055op9xTTk+6h+50N9XF1SwuXzyojHSufs5Q/jX8rzxy7BE8Dg/FzmI6Uh20J9pxO9zsb9lPOpMec857/vSS80vns9e5l3JVTlu8jSbdRJm7jOsWXceKyhW09LRwLHqsf8CV1+llTc0a9rfupzfby41Lb2R5YHl/GOfdi97NgZYDpDKpQUXUpqvi57lw1mk3Lxz+PVprek72SAftB5ax7IPL8FX7qP3zWl74xgvs3rIbYFpE3wi+wWCRi9faimzM882jp7en33NfWLaQquIqDrYe5PXm1/E6vThtThQSow54ArgdboKlQU50nhjUD7A8sJzdkd3YlI09p/bgtDtpjjXzvsXvY2/zXhq7Gqnx1XBl8Eoi3RHqfHWD7Bop7v3Lw7/EZXfhc/noSnfhcXjotfVyousEvZleytxl55QRk1/cLpqM8lLjS2il8Rf5qfBU0NjdiNvhZnnlcl6IvEBtSW3/e50OJ6H6EJ2JTi6qvmjYjKT5/vln9Hk8fPhhAt5AQVdQPZeif1prHrnrEXZv2T1I7HPhmusfuB5g2kTfCL7BYJEfr60vrSeVkaHuF5RfQF+2j45EByc7T+JxesiSxamcpHWarM7SEm/hgrIL+mvm5Hf2VXorWR9cz0uNL7G7YTdXLbiK6xZfh9vpZl7JPOb55pHVWZRSuB1uno88z9HoUS6uuhi3c+S4d0+qB4/dQ0ZnyGQzOGwO7NhJZGXk7ngyYrJkWRpYOmiik7Z4G+3JdkJ1IXYd3kVHooMKTwWJvgQ96R5WVa1iSfmS/jpC+eWQK9wV+Fy+QcdI9iV54tgTbFy+cdrr/E8USimwwcJ3LWT1x1ZTVCLZN/kF0XKi7y5zT7mHb0orGAwWoboQ0WSUrlQXS8qX0JZoozXeymW1l7GqahV92T5SmRQXlF1Api9DhgyZTAaykM6kiaaivNX6FptXbu5vPHJUeitx2B39Yl/tq6a0qJTF5YsJeAPctOwm0pk09aX1vPeC9wJSVTPeGx9RAJdULEEjJQsUit5sL8m+JD6nj2UVy8aVEVPuLiejM8R742itiffGyegM5e5ygv4gd1x2B33ZPk7FTuG0O1lVtQq7zc6NS29k04pNeJ1emmJNeJ1eNq3YxPLK5YOuB8C+5n3U+GooLSrFpmyS8mjdlcxkFl6zkMtuvwy3342yKYpKi4i3xfnR5h/xj8F/5Bt136Dx5UbmrZs35bYZD99gsMiP1/ake7hq/lUoFJlshvn++Xx45Yf51ovfojvVLYXR4hkcykEqm8Jn81FZXEldaR2hYIjaztozOvuaY81ct3hwxfCcF54fAiotKuW6xdf1TwE5krf7mdBn+NITX6Iv20dfpo+edA8Om4ObV9yMy+EaV0bMisoVeJ1S5K0j2YG/yM/CqoX9efahYIjaktoRB54NZ/O5XI+ZTE9LD76agbuZtkNtvPbfr9F9shtvwIvD6aDtYBuP3f0YABduHKEzYBIwgm8w5HG2eG2kM0LAE+CJo0/Qmmil3lePsikSfQk2BDdQU1zTX/N+uM6+dy9+t9SyzyPnhQ+X7z2a+G28cCMAW1/ZSmN3I6s9q7ms9jIWly3G6/QOdC6eB6G6kOyzevWg0Ex+kbZziW2f6/WYyRRXFZOOpSkqlXTbk786SbJFSiK7iqUstrIp+uJ9vLL1FSP4BkOhkUvZrC+t52MXfYx/Df8rke4I80vnsyG4gbqSOtribYMm8x4qiLl9wJnz+oYbw8OWXxhN/DZeuLFf+CeSs2anjGOfY70eM5m6UB0HdxwEpNRx9+lu+nr78AQ8/dvYXDb60n10N3ZPqW1G8A2GMTA05PKH6/6Qn7z5E/xuPzXFNbTF20adzHs0ES008ZvsKSkno1EpBPxBPys2raAx3Ch5+tXFFPmKIDuwTTadxaZslNSVTKltRvANhjEwNOSyvGo5H+EjPHH0CU50nqC+tJ7b1tw24mTeOUYS0dkqfqNxro3KWAvWTTf+oL9/INa8y+ax665dtL/ZDlrCOb3xXtzlbi67/bIptUtpraf0gGdj3bp1es+ePdNthsFwBkNr2gP9napDSykYxkc4EmbbgW2DqojmRi7nBrPl3wXNhDTOSDjCs3/7LKdePEU2kyWwPMCGL2yYsPi9UuolrfW6Ubczgm8wjM5kic1M8VininAkzP3P30/AE6DMXUY0GaUt0cbdG+4m0h0xje4IjFXwTUjHYBgDkxFyOZcSy3OFbQe2EfAEaI+388SxJ+hOdeOyu/j6c1/n/UvfP2omk2lAz44RfINhjEx0J2a4MUwmm2F/8346U534i/zM8807p3IIWuuzjtYcbX2h0dDVQLovzTMnn8Hj8FBaVEqiN8ETbz/BisoVFLuKR8xkMg3o6JiRtgbDNHGw9SD7W/aTyqQod5eTyqTY37Kfg60Hx/T++566j7seuYuRwrJaa+565C7ue+q+CbR6cqkvrefFxhfxODx4nB5sNhtOu5PSolL2NO7pHwmd1Vm6Ul2DxgbkZ1LNppG7E4kRfINhmuhIdmBXdrxOL0opvE4vdmWnI9kx6nu11kSTUbbs3jKs6OfEfsvuLUST0REbhUJj88rNdKY7UVqhs5p0X5pUNsUF/gv6+0yGlm3Iee8tPS1n1OvxuXy09LRMx6kUJCakYzBMExXuCqKJKPHeOG6Hm2Rfkkw2Q4W7YtT3KqV44PoHANiyewsAD1z/AEqpQWL/ufWf619eSIwUaw8FQ1xceTHHu44T75PrssK/AhTUeerOGlYbae7gmT5ydyKZEA9fKfUdpVSzUur1vGUVSqnHlFKHrMfyiTiWwTBbWF65nNXVqylyFBFNRilyFLG6ejXLK5eP6f050f/c+s8N8vRngtjvOLiDeG+cGl8N8d44Ow7uINIZAeCLV32RSm8lyyqWcXH1xaDkbuj2y24/637zi98NF/IxTFBaplLqnUAM+A+t9UXWsr8H2rXWf6eU+hJQrrW+52z7MWmZhrnERKV65ot8jkIVexgY05DuS3Oo/RCdqU6cdierq1Zzx9o7ANj55s7+GkF1JXXcftntYyohMVezdKY8D18ptQjYmSf4B4F3aa1PKaVqgae01ivOtg8j+Ia5xkQJlNYa218M3LBnv5odVeynSxwf3PMgdmUn3Bim2FWMx+Ghp7eHplgTf/++v58TAj3RFEIefo3W+pT1/DRQM9xGSqk7gTsBFixYMNwmhjnGXPLSJiLVM+fh53PXI3ed1cOPdEZ46NWHaE20ku5L43K42N+8n1vX3Drp17qquIqnjz/dP0UigE3ZqPZVn1NKquHcmZIsHS23EcPeSmitH9Rar9Nar6uqMp0rc52R4rs739zJPY/dw8d/+nHueewewhGTagecEbPPfjV7Rkx/OHYd3sXhjsPYlZ0KTwV2Zedwx2F2Hd416TaH6kI0xZrQWpPVWeK9cXrSPVxcffEZGTWRzgjbD2znwT0Psv3A9v44v+H8mEwPv0kpVZsX0mmexGMZZgn5udQgc58eaT/Cd17+DpW+SmzaRkuihf3N+7n32ntHLVY2mxmpg3ak7J18Xm58mXJ3eb+H7XV6yeosLze+DGsn1+6gP8h1i6/j9ZbXiSaj+N1+Lqq+CJfd1W8PmIFUk8FkCv7PgVuAv7Med0zisQyzgEhnhEcPP4pSijJ3GcsqltGR7OA7r36H5lgzrg4XLpsLv9tPsauY77z6nTkr+GfLxhmL6GulUQxuBBQKrSY2Xz8Xnnur9S3ak+2Uu8tZUbmCS2suJZVJndFhnV8OOr/xb423cqj9EEfaj/D4kcdZH1xPwBtAocjq7KwP/U0UEyL4Sqn/Bt4FVCqlIsC9iND/WCl1O3Ac+OhEHMswO8l5c0XOImzYSGVS/OKtX3Cy6yStPa30Zfqw2+wkMgkcaQexdIwXTr4w3WZPG7lGcaRsnHzRL3OXnbF+Xe06nj3xbP/E6cm+JG2JNq5ZcM2E2Zj7TDPZDEc7jmK32Ykmo3idXhq7G7mi7goi3ZERaxPlSlK3xlvZHdlNVmeJJqIkM0n2Ne0jnU3jdXp558J39of+jPd/diZE8LXWvzXCqutGWG4wDCLnzV1acym7I7spdhVzovME6UyaTDaD0+6kyFFEX6aPZDaJz+GjIzX6iNTZzH3vuu+stXJyoj/c+huW3kBTrImWRAsdiQ5cDhfLK5Zzw9IbJsy+3Ge6v2U/viIfXqeXeG+c07HTrK5eTaQ7ctYql7mBVIfaD/V/H2w2G9Xuarp7u0HDAv8CjnQcYUNwQ/8xjeCPjBlpaygIct6cTdlYH1wv+dnJTkqKSvA6vaQyKTLZDDZspPvS9Nn7qHZVT7fZ085oqZcjrQ/6g9yy5pZJzYbKfaadyU7K3GUAeBweOpId/VUucyGfg60H6Uh2UOGuYHnlckJ1IUJ1IXYc3MHp7tPM880jmoxiUzZqfbUcbj8MgNvhJpqMArNjAvTJxgi+oSDIHxZf6a2k0lvJa6dfoy/bh1KKSGeEdDZNOpvGYXNQ5a1iw4IN0232jGaypzDMfaZ+t59kXxKv00uiL4G/yE8sHcOmbDz06kMcix7rD9G4lIuL513MG81vcMuaW9i0YhMnOk/Q1NNEaVEpAU8Av9uPw+4ADcm+JH63zCxlyiiMjhF8Q0GQ8+ZgYE7X1VWreaPlDRb5FxFNREllxcufXzqf6uJqfnfN706z1YbhyPfaj0aPEiwJ0tbTRqI3QUZnWFi1kGgySiwdY2/TXvY376ehq4EsWezKTne6m85kJzWHa/jU2k/x2Ss+298XsL95P23xNkqcJaSzaVrjrbxz4Tv7yyjM9AnQJxsj+IaCYLgJRj67/rOc6j7FtgPb6OntoSPRQZm7jDW1a/qnvZsrzJTBaPmplCurVuJ1enm9+XUCngBZspS7y5nvn0+oLsSXH/syDd0NNHQ3kNVZHDYHfbqPU7FTlHvKefrtp/nU2k8N+m4kehO0J9tZ4F/Qn6WTyWbwOr1zYg7g8WIE31AwDBdiyFVQnMtEOiNseWELB9oO0J3qpqSohF+f+DWfu/JzBSdwQ8dRLC5fTMAbwOv0EqoL9Tda4cYwJ7pOkOhNAGC32bHb7ZCBDBkSfQlaelrOaOh+55LfKbhznkmYevgGQ4Hz/X3f59mTz5LJZqj2VpPJZnj25LN8f9/3p9u0MxipJv1brW+x4+AOTnae5Gj0KDvf2klDdwNtPW3YlR2tNZlshozO4LQ5iafjeJ3es1bVNJw7RvANhgLnyaNPDszkZBuYyenJo09Ot2lnkOuozSeWjtGebCeTzfBGyxv0Znqp9dVSX1pPMpukyFEECjQau81OsbOYCk8F80rmjWkGK1N+YeyYkI7BUODEM3HKHGWDlrlsLqK90Wmx52wM1/keTUYpd5fTGGscVDDtkupLaI+3k+pL4XF6yGaz2JRNSimU1HCq+xSvN7/O8sByKr2V/fscOmn5TAl3FQJG8A2GAueiqot4rfk1vH1eYukY8b44yb4kodrz79uIdEbYdXgXLze+jFaadbXruGHpDeMWyeE6369deC3hxjA739pJra+2f1uXw9WfYaPRtPa04i3y4ra7uTJ4JY2xRjqTneyO7GZ9cD2V3sozUi9z4a5ydznV3mpivTGePfkslcWV3HP1WaffmJMYwTcYRqBQMmNuW3MbX378y5zsOonL5kLZFB6Hh7qSOiKdkXO2KVca+XDHYcrd5SgUz554lqZYE7esuWVCRH+4few6tIv2RDvlnnKSfUl60j2sqlrFfP98nMrJ1le2su/0PkrdpSwpW8KKihXsbtiNTdl4q+0tXHbXGamX+eEukGJ7Gs2TR580gj8MJoZvMAzDaNPwTSWhYIj3L30/S8qXUO4pZ0n5Ej5xySe4ZN4lZ8Szx0K4MUxropWAJ4DP5aPYVUylt5KWRMt57W8sBP1B7rj8Dkm77D6Fy+6i1lfL682v8z9v/A9/9Mgf0dzTTImrhN5sL9/f/32Oth9lWfkyTsVO8b9H/5fXml7jirorBjUm8Uwcl8016Fgum4t4Jj4p5zHTMR6+wTAMw5Vpzi2fDi+/zF3GHWvvwKbyZrXS2fMqJdDS00K6L02FZ2CydLfDTUei44x69BNB/p3SFXVXoNG0xds4Gj3KRdUX8YN9P8Dj8JDMJLHZbHjtXhSKnYd2sj64nlpfLRdWXshF1RfxYuOL1JbU9n8GuXBXpa0Sl81FOpumPdnOJdWXTPh5zAaMh28wDMNI6YWTIYhjYaTsl/MpJVBVXIXL4SLRl+hfluxL4nK4Jrw0wdA7JY/TQzqTJuANcO3Ca1lcvpj2RDv+Ij9OmxNAZuCyuTjdcxqbspHVWZYHlg+bpXPbmtsoLyon0ZegO9VNoi9BeVE5t625bULPY7ZgBN9gGAalFE++/SQPH36Y508+T2tP67TWagnVhYgmo3SlusjqbH8pgVDduXfchupCVHoqaUu0EUvH6En30BpvpcpTdV77Oxv5d0r5qZUvN77c36BWeCro6e3BYXOAoj9d02Fz4Hf7+zts4cxGNxQM8dVrv8qamjUEvAHW1Kzhq9d+dc4P1hsJE9IxGIYQ6YzQHGumI9FBwBMg0Zfg6RNPs7R8KbeuuXVabBop++VU9ym+9eK3aOhqoL60fkwlJ4L+ILeuuXVQls41C66ZkCydoeQqZubjc/nQShNLx0j3pVlQuoAnjz+J2+GmylNFojeBy+biQ8s/RKIvwZ7GPfjdfpZVLMNlP/MuJBQMGYEfI2qkOS+ng3Xr1uk9e/ZMtxmGOUykM8K3XvwWLT0tFLuK0Wi01jjtTi6quohPrf3UdJvYTzgS5v7n7yfgCVDmLiOajNKWaOPuDXcXjABuP7CdeG+8vw8EoCvVRaI3wenY6f5MoaPtR3nl9CtkdZbVVav5rdW/RSQWGZRJ1JZoY3nF8gnJJJptKKVe0lqvG2074+EbDBa5dMU9p/bQk+yhJdFCNptlZdVK3rP4PWR1drpNHMS2A9sIeAIEvAGA/sdtB7YVjOCPNBBr04pNPHz4YVoTraQzaVbXrOZDKz80aF5bp9NJfUm9zI2Q6qTcU061r9qI/Tgwgm8wWOw6vIvDHYdp7W7lcFQm2LArO2+2vUl7op0Pr/zwNFs4mIauBhb4FwxaVuYu40TniWmy6ExGCkUF/UGyOsu7F717xMyjGl8NtiIblcWVZ6wznB9G8A1znlza4H/u/U+yOsvR6FEUCqfdSSaboam7CY/dwxNHn6DMXVYw5YnrS+uJJqP9nj1ANBmlvrR+Gq06k5EGYuVPepMjv2P8bOsM54fJ0jHMafLTBj1OD41djfTpPtx2NzZlQ9kUKBHSWG9s2gdh5bN55WbaEm20xdvIZDO0xdtoS7SxeeXmabVrrJwt82gis5IMAxjBN8xp8tMGlweW94/c7NN9uB1uiuxFlLhKyOgM84rnnbVq41QTCoa4e8PdlBSVcKLzBCVFJQXVYTsauXCP1+mlKdaE1+ll04pN/XcEI60znD8mpGOY0+SnDV5eezlPHHuCtmwb3eluXFkXbpubVDaF0+bkHcF39L+vUCbMnukpiWebV3ey59ydixgP3zCnyR/BWumtZPOFm/G7/VQXV1NWVIbD7sBhc3DjkhtZXrW8/30mnmyYiRgP3zCnGZo2WFlcyarqVbTH21EoVlWt4gPLPsDxruN0pboGpRaaCbMNMw0j+IY5TX7a4Jstb3Kk4wjvWfQeFpYt7Bf2NbVrWFO7ZtjUQoNhJmEE3zDnycWKtx/YTtAfHLZC5s0rbzYCb5jxGME3GCxyHbgvnnyRx4491j9Zx9p5a7l55c3TbZ7BMG4mvdNWKXWDUuqgUuqwUupLk308g+F8qSqu4uljT/P9/d8n0Zeg0iNT6j189GF2vrlzus0zGMbNpAq+UsoO/AtwI7AK+C2l1KrJPKbBcL6E6kI8cewJPA4PJa4SsipLkb2Ii6vLeb35L4EHge3A9A64MhjOl8n28K8ADmutj2qt08APgU2TfEyD4bwI+oPYbDZKXaUkM0nsys7FNRV8YHmGeG8rUAPEgR0Y0TfMRCY7hl8PnMx7HQHW52+glLoTuBNgwYLBhaAMhqlmkX8R3enu/vo0F1VHaO3ppchRhfhHudouYcB04hpmFtM+8Epr/aDWep3Wel1VlRnIYphebr/sdjqSHf31aeyqnZZ4gvcteX/eVj5geqY6NBjGw2R7+A3A/LzXQWuZwVCQbLxwIwBbX9nKyc6TpGtL+cQl13D5oKJdMcA4J4aZx2QLfhhYppRajAj9x4DfnuRjGgzjYuOFG/uFX6KQO4AuxLOPAVHAjLI1zDwmNaSjte4DPgM8AhwAfqy13j+ZxzQYJpYgkmfgBZqsx02Y+L1hJjLpA6+01r8EfjnZxzEYJo8gRuANs4Fp77Q1GAwGw9RgBN9gMBjmCEbwDQaDYY5gBN9gMBjmCEbwDQaDYY5gBN9gMBjmCEbwDQaDYY5gBN9gMBjmCEbwDQaDYY5gBN9gMBjmCEbwDQaDYY5gJjGfQCKdEcKNYVp6WqgqriJUFyLoNzVYDAZDYWA8/Aki0hlhx8EdxHvj1PhqiPfG2XFwB5FOMxWewWAoDIzgTxDhxjBl7jJKi0qxKRulRaWUucsIN4an2zSDwWAATEhnwsi+uJt1zx7E29xBoibA6fe/g+zFS2mKNU23aQaDwQAYwZ8Ydu7knf/0P2SUIlNZhrupjQu+u53Xfud9VF1+ydTZEYlAOAwtLVBVBaEQBGdBH0I4DNu2QUMD1NfD5s1ybuPZZSTMtgPbaOhqoL60ns0rNxMKjm+fBkOhY0I64yUSgS1bKC4J0FXuJtObxtXRRUpBzSPPEaqbIhGJRGDHDojHoaZGHnfskOUzmXAY7r8furthwQJ5vP9+WX6+u4yEuf/5++lOdbPAv4DuVDf3P38/4YgJvxlmN0bwx0s4DPE43po6gv4gNrebuC1LUbqPC5MlU5elEw5DWRmUloLNJo9lZeMSxoJg2zYIBOTfbh94vm3b+e/ywDYCngABbwC7zU7AGyDgCbDtwPnv02CYCZiQznhpaRGPuqcHb0kJXqcXSjQ0NcHiZVNvRz4+n9gxk2loEM8+n7IyOHHi/HfZ1cAC/+B9lrnLONF5/vs0GGYCs0fwpyt+XVUFa9bAY4/J6+Ji6OiAvj6JNU8m+ed85Agkk7B48cD6WEzsm8nU10M0Kl59jmhUlp/vLkvriSajBLwD+4wmo9SXnv8+DYaZwOwQ/Fz8uqxMvNxYTF5v2iTrJ7MhCIWgsRHe9z549VV57vXCl7887o7FYcmJ/MGDcPQoXHSRNDKtrfDLX0J1tYhhRYWc7y23TLwNU8nmzRKzB/l8o1Foa4Pbbjv/Xa7czP3Pyz7L3GVEk1HaEm3ctub892kwzASU1nq6behn3bp1es+ePef+xu3bpZOytHRgWVeXLEunIZOBU6ckvOFwwB13TKwYT9XdRX7Dtm+fNC7Hj0MiAX4/tLfLuZaWwrJl4u3feuvMz9TJZem8+qoIvtstDdqFF8IVV5zX9TZZOobZhFLqJa31ulG3mxWC/+CD4tnbbOLpHjokwnD0KGzYAIcPQ2enhFkyGQkP/P3fDy8SkQg8/DDs2QNKweWXw403FoZobt8OJ0/C6dPw5JMi9L294uGn0yL4DoecX3GxhJpWr5YGbiaS35B2dMBzz0lj19wsd3FtbdKw1dZOfCNuMMwgxir4syOkU1Ulnu5bb8Frrw2EdtJpeOIJ6OmRDI9MRh5bWkTUP/WpwfuJROB735P9BAKgNfzqV3JnUAie8sGDcOyY2NXWJncwPT3gdEpjBtLo5Yvik08WToN1LgwN0/3859LQNTZKQxeLQTYL+/fL5/zAAyM34gaDAZgtaZnBoAjzoUPgcokX/OijInr79klj0NsLHo+Ifm8vPPXUmfvJeZOVleIh+3wi/K2thZHe2NEhYnfihHjyNpuIXVeXPO/thVRqYLu9e+WcC8H2c2Vommlrq+TgHzsmn29Tk3xWLS1yZ/Pmm7Br13RbbTAUNLPDw49E4Kqr4Kc/FWFIpUQcolHxBF0uieE7nfJfXS3bDaWlRQS0vHxgmccjgtLSMmWnMyIVFfDKKyL2uQ7MbFaEvrtbvHylxPufP1+WpVJyxzLTGJpmmmvI0mk532xWlvf2SnYSwDPPzNzw1QQQ6Yyw6/AuXm58Ga0062rXccPSG0zFVkM/4/LwlVK/qZTar5TKKqXWDVn3ZaXUYaXUQaXU9eMzcxRaWmDhQvHMS0tFGBwO6dxzOEQQ4nERieJiEXCv98z9VFVJ45ATEBBP2eUqjPTG5culMcrZHotBUZE0Srm+GKUG1tXXy/Vob58ee8dDVZWcQw63W8Te5ZJzVUoe7XYJawUCckc3R4l0Rnjo1Yf41clf4Xa68Tq8PHviWb736vdMxVZDP+P18F8HPgz8W/5CpdQq4GPAaqAOeFwptVxrnRnn8YYnJw4+n3jyIEIPIhAgwuBwiDi63TBv3pn7CYXgjTcGx/A7OmDp0unpEBxaQ+aaayQbJ3dukYg0ZCBhD5tNvHyvV7az2SSElX/HUiiMltkUDMK//7ucT65DvrRUzj2RkPMqKhpo2G02uXObo4Qbw7QmWgl4AjL4D1BK0ZJoIdwYNl6+ARin4GutD4B8sYawCfih1joFHFNKHQauAJ4fz/FGJBSSDr7ychF2p1NCAD6fCIHTKd6hzydCX1Jy5uhNEJG55ZbBWTpXXTU9nZ65GjKBgIRv9u6VDugNGyRkc+qUnIfHIx5uRYWIZy6ElStFsHq1hHcKiXBYOllPnxZb6+ulob3lFrnOkQi8+KKMMWhsFM/dZpOGobpaPudTp+Qz9njkri0QgHe+c7rPbNpo6Wkh3ZemwlPRv8ztcNOR6KClpwDCkYaCYLJi+PXAC3mvI9ayM1BK3QncCbBgOBEejZynGIuJAIB4t1oPhGNAXvt8IvTLl48sgsGgZO8MzeCZanI1ZBwOyTyKx+V8fv1reO975RwPHpT1bW0DnrDfL+e3dKmEsOz2qbk7GetYhEhExP7EiYHRs2++CQcOyOOHPyz9K7kO29zI4WPHpMHr6hKPvqRE7mQqKuDii2W7G2+c/PMsUKqKq3A5XCT6Ev0efrIvicvhoqq4AMKRhoJgVMFXSj0ODBP/4Cta6x3jNUBr/SDwIEge/jm9OT9178ILBwQmHBZhcLnEu4/H4dJL5fUbb4i3v2HDeE2fXHI1ZPbulbCSxyNi3tEhnvHq1eIZHz4s2x05IqK/eLGI4KlTcN11cMMNk3N3kh9uKikRcb7kkjNHOg89djgsGTaBgAh3PC6dy3a7ePLxuKSSXnfd4IF0CxfKZxYIwMsvD4TwFi6EFStmTyno8yRUF2J/834Odxwmq7MoFG2JNpZXLJ+6iq2GgmdUwddav/c89tsA5LvQQWvZxJKfugfyePXVInpPPCGdr5WV4s0fOCAestMpoYIXX5QBO4UqErkaMqdOiTDmGq7ychG9o0fhz/5MUhFfflnOY6oEMD/ctGCBHD8aFbEvLR34PMLhM21oaZFzydHeLo1ZMikhtNJSCdvs2yeinyMWk/O6+ebJOacZTtAf5NY1t/Zn6XSnuyktKsVhc/TPumbi+IbJCun8HPiBUuofkU7bZcCLE36UkSpElpXBpz89UG7h+edF5EE6+hYvljuA4QSpUMjVkEkmxfvNhXTWrBnIUgkGpycNMb9kMYh95eUyFmKZVSF0pEqdVVXSmL39tpxDMinnk81KGAokRPP44/IZ+Xwi9tEoXHvtVJzdjCXoD3LH2juILJX5lcvcZfhcPmLpGDsO7mDTik1G9Oc4403LvFkpFQE2AP9PKfUIgNZ6P/Bj4A3gYeAPJyVDZ2jqHgxUiAyFRCS6ugby1Xt6BgtSIeTWj0QoBHffLXcnbW3SaXnlldKYdXRIyYfpoqFB7Mjh9YrX3tY2sGykSp2hEFxwgTQQ2ayE15JJEfu1a2Ubt1u8e69XGg2vd/jwkGFYzPzKhpEYb5bOdmD7COv+Gvjr8ex/VHLZOXCmJxgMikiEwyIsWsP69RLigZlROjgUgq1b4aGHpCMzVwhu6dKxdVBOVlG3oSWLa2ulsFlpqVzrs3nkuUyomhrJhPL75ZyuvFL6HnINtBH486alp4Ua3+A7X5/LZ+ZXNsyC4mljEbX8zt38hmGmiMr5CPdknvPQlNFoVModXHONvD7XxmW2zsU7TWw/sJ14b5zSooFO765UF16nl5tXmj6Q2cjcqpY5FmaTqORi+Gdb/z//M3zJaK93Yjo+J2FiccPEEOk8M4YfTUZNDH8WM7eqZY6FYHDmCnw+990nHvUDDwwv+lrDXXdJmubQsQQTOeVhKGQEvkAJ+oNsWrGJcGOYplgTVcVVXLvwWiP2hjkk+LMBrUXst2yR10NFPyf2W7bABz4ghdb27ZPO1EBA8uQvuWRaTDdMLUF/0Ai84QyM4M8klBKRhzNFP1/sP/c5yXK57z7JhqmslIbi5z8fyIQxGAxzDiP4M42RRD9f7B94AP7gDySE094uaZzz5sno3GefhY0bp89+g8EwbRjBn4kMFf2c8OfEvqEBXnpJRD4381cqJR24DRM/4NlgMMwMZseMV3ORfNHPkQvvhMNSniA3IUpRkfwfPiwZNYY5yk7gZmC99bhzes0xTDlG8GcquZh9PnfdJctbWuBd75KiZN3dMhgqlZLBW5s3T4u5hqkkgoyHfNB6jCDifi/QjZS56rZeG9GfSxjBn4kM7aDNZuVxyxZZXlkpo19vvlkKk7W0SL2bj37UpFLOeiLADiAO1FiPO4BvAeVAALBbj+XA1ukx0zAtmBh+oTN0gNNHPgI/+MHgDtqhMf1YTIqs1dTAJz85eJStYZYTBsqA3IC7NHAI2ANUAk2AAnqsdX3We4wjMBcwgl/IDC1D3NEBH/84HDo0WOzhTNG//Xa46SYZaFVVNVBfyDDLaUE8e4BWYDeQsf6bAQ2kADdQDBQB9wN3Y0R/9mMEv5AZWoa4slKKja1dO/xI23zRLyuT2aMMc4wqIIZ4+IcQUX8N8AOdiNhrxLvPAu9HwjvbMII/+zGCX8jkZr3K533vg+PHR66lkxP9s9XaMcxiQkjMHiCKePCNwHIknv8SIvhOJIZ/IeL9n5hqQw3TgOm0LWRyZYjziUZHD80YsZ/DBIFNgBfx4DUSu/da64LI9NIrEMEHaRhMuu5cwAh+IbN5s9TBaWuTmvG55ya10nBWgkie/Z8hk80tA7qs/9wcEB1IrL/N+jffqbmAEfxCJjfrVUkJnDghj3ffPTGplTt3Strm+vXyuNPkY88+ct7+lUhM3w4sRcI+85CO2xJMh+3cwcTwC53JKEO8cyfce68UVps/X8JE994r60ydnVlGEPgUcAOSftmCePkha51hLmEEfy6ydauIfS77J/e4dasR/FlLLn5vmMuYkM5cpLFx8CTkIK8bG6fDGoPBMEUYwZ+L1NUNn/1TVzcd1hgMhinCCP5c5PbbZdRufvZPR4csNxgMsxYTw59sCnGy71ycfutWOHlSPPvPf97E7w2GWY7SWk+3Df2sW7dO79mzZ7rNmDjya+GUlUnYpK1t4lIrDQaDAVBKvaS1XjfadsbDn0xytXDa2uCJJ6Q2vcsFf/u38LOfTbd1BoNhjmFi+JNJQ4PUon/6aZlmsLRUatk/+aQZ6GQwGKaccQm+UuoflFJvKqVeU0ptV0qV5a37slLqsFLqoFLq+nFbOhOpr5ewjscj/zYbOJ0yYnarmXjCYDBMLeP18B8DLtJaXwK8BXwZQCm1CvgYsBoZ4vevSin7OI8189i8GTo7pZiZ1gOTiS9ZYnLeDQbDlDMuwddaP6q17rNevsDAUL5NwA+11imt9THgMHDFeI41IwmF4JJLZDLxeFymGVyxQh5NzrvBYJhiJjKG/7vALut5PXAyb12EEeqvKqXuVErtUUrtaWlpmUBzCoQvfEFmnFq2DC6+WJaZnHeDwTANjJqlo5R6HCmtN5SvaK13WNt8BZkc8/vnaoDW+kHgQZC0zHN9f8Fjct4NhklgJzIBeyNQB9wOmN/UaIwq+Frr955tvVLqVuRKX6cHkvobgPl5mwWtZXOTjRsnR+AjEekUbmmRu4hgUJblXodCZh5bwyxkJ3AvMoHLfGQCF6vaqxH9szKuPHyl1A3AF4FrtdbxvFU/B36glPpHpPldBrw4nmMZ8ohE4OGHJbe/pkZCRSdPwg9/CFddBQsXQiwGO3bApk1G9A3TRIQzSzIzzLJz/X5uRcTeqvLa/7gVI/hnZ7wDr/4ZmTTzMSXT6r2gtf601nq/UurHwBtIqOcPtdaZcR7LACL2O3bAoUNQWysZQOEwOBwyyOv0aUn7PHQImppk4pTPftaIvmGKiSBz65YhM2vFgIcABSzKW/Y9oBqZinGsDUAjgwMIWMc5eeamhkGY0goziUgEvvUtaG2VQV3LlknJhngcXnpJBL+hAfx+8fJdLnjjDbkLuO46uOEGI/yGSSTn0TcDR5E8jcXWulbgp0ACSdhbbi1/CvHW3400AFEkye9s39ObgW4GPHuQaRpLgO3jPouZiCmtMNvIefYtLTBvnjy+/LIIfiIBhw9L+qdSMsDr4EGphFlRIXcCr78uYwBMiMcwKeQ8+ieAXsADvALUIkIcBZKAE3gEGcJzJSLaKSRhsNTaV5izC/7tDMTsy6x9dwCfn5Azmc0YwZ8phMMi7vPmiXBXVsKbb4p373BAUZGEcRwOaG+Hri4Z1btihYzyjUbl/eGwEXzDJBAG/EgEdyewAJk/twcR4x5ru7eAY8BFwK8QwfcjdwCVgA9oGuVYuTj9ViSMU4eI/fnF73e+uZOtr2ylsbuRupI6br/sdjZeODv7AozgzxRaWmTAVnc37N0rAm63y/JUSqYsDAYlZt/TI6N6fT44ckQ6dO12qeUzdKYrg2HcRIBHkfj8auA48Jq1bimQRUT/baATCefUWu/LIqL/v4jwJxDhj3B2L38jE9FBu/PNndz79L2Uu8uZ759PNBnl3qfl7mE2ir4pnjZTsNngmWckLn/BBRLH7+iA6mrx9uPxgUYhV8ohGpUSDqdOwdKlUubhyBEJDxkME0IulFMEuJFYehqJ358ADgKnkazsnNgnkcH3biTEcwxJ4nvOWldv7XPyv6dbX9lKubucgDeA3WYn4A1Q7i5n6yuzs9aV8fALhdEmStFahLy7W+Lzdjtks5J+2dkp8fu+PmkQbDZ5jMVkG78fjh+Xx1DIhHUME0gYiaNfisTvTyHiX4R4/PmiHUCyc5oRsb8QEf6ktb1GfNASwMXosfzx09jdyHz/4IyfMncZJztnZ8aP8fALgdxEKd3dsGCBPN5/vyzPoTWsXi3C3dEh4g3i1Xs8IvA9PRLe0VrCOSUl8p+/j+JieY/BMCG8BbwO7AHaEe++F4nDLxyybQAR+DjSUXva2jaKNAAB5C7gR9Y+35p06+tK6ogmo4OWRZNR6kpmZ60rI/iFQG6ilEBAPPfc823bBrapqhLv/6KLZKCV1yvLHI6BKpxeL7jd4uE3NkqWjtbS0XvppZKeuW+fvM9gGDcR4Agi0mVIqmQbkle/Bjg0zPYupGNWA/uRGL5CYvwnkRi+zdrnESY7rHP7ZbfTkeygLd5GJpuhLd5GR7KD2y+bnbWujOAXAg0NZ3amlpXJ8hyhkAykygl4W5tk4rjdkpbpckmmTjIpr51O8fJdLrljiEblvU1NZnpFwwQRRrJtskhYpo8B8T6BhG5KkDz7CsSzb7f+KwFv3vMi63UnEmnOWvvOu8udBDZeuJGvXfs1SlwlnOw8SYmrhK9d+7VZ2WELJoZfGNTXiyAH8gaSRKOyPEcwKIOnXn9dhL6oSEI5iYSIvNstnntTk4h8bl+VlfL60CG48ELZh4nfGyaEFiRsUwK8hHj3ccTT70KE3I3ITC+SZ9+FNAwHGcjKmY+Eejqs7VYC65FGYrQUzfGz8cKNs1bgh2I8/EJg82bx2NvaJAyTe7558+DtbrhBRteWl0vNnOpqCeMsXCgNQHe3CHxVlXTgOhywfDmsWyeNx7Jlsg+DYUKoQkbHgoj4Cut5F+LVX2S9tuc9+hCR77IeqxHxzyAx/IuQQVtY+zbhx4nECH4hEArB3XdLB+uJE/J4991nhl6CQRkpm0pJyGbVKvHk582TRmDRIonv19VJY/DJT0oDcPq0PJpRtoYJJYR0uO5FRLoWCcWUI9k1bYjEuK3HLOLlVyAx/2IGausErdeLrX3ttfZ9tvBjBCml8KD1aNKNR8OEdAqFUGhssfVgEN7/fsm7Ly2V/1RKvPu2NvH429rEs1+yRDpq6+uN2BsmgSBS9+avEEEvQsolnECydToRD91nbV+FxPU9yN1AFngT8fxLkBo5bYjQZzl7TZ3hirPtGOU9BiP4M5FQSOrqgIj6s89Kh+xNN0ks/9gx8fqbmiS8c+21RuwNk0QQuB6J3e9HPPT5wD5rfTki8lGkQXAiYt6OiPUixLO3I57/ciTc4+Xswp3L/8/V3xlrHZ65jRH8mUgutBMOS+79VVfJoKxMRjz8W281Am+YQkKId32agcnx6pHUy+OIR/8OJEzTi8T709b2yxHxX4zE8l1I43DtCMfKVeT8KdKwLEc6h2FsdXjmNkbwZyrBoBF1Q4GQC+2cQAS3BngPkrmTQTzxLFLkrBVpCEoRr96BhHJakDz89YjYD/fdzg/jlCLVOH+FNBYljL0Oz9zFCL7BYJgAgsBnGRBkHxLmmY+EdV5HxPgiJAVzFeL5dyAT4tUgdfJvPssxcmGcNJL6mUb6Dl62li9goA6PieUPhxF8g8EwQeQ8/TDi6VcxMAlKKVJOAWAJUi45wUBufpSRwzg5WpCGYTeSzlmBFF3rZWBw12JrfyaWPxxG8A0GwwQSZEBocyGYLkTkn0XCOe9EBP91xDP3MnIYJ59c3n8nIvAKuZOoR7J+otZ2JpY/EkbwDQbDJJHv8fcAVyEinUFCPR/m3LzwXOew09qfzfqvQEo7WAUFzYCtETGCbzAYJpF8j38i9rUJ2AU8iYR13oukg7Yidw5jDQ/NTYzgGwyGGUQQuAO4EblzaGHwncNYw0NzEyP4BoOhwMnl3rcgoRon0h/QgMTvN3P2EgyGHEbwDQbDGAkD25hcoR0q7kHgx8ABJBWzDwnfrEfSMKPA/cDdk2DL7MMUTzMYDGMgjAhrNyK03dbriaxXn8vqiSPpl3GkTs/jSLimGmhEcvcjyMCtgPW/bZj9GYZiBN9gMIyBbQyIa05o+4DPAB8H7mH84p9fH8dmPR5jYGRuruJmMfB23vvKkLsOw2gYwTcYDGOgARHW/NdvIHVwJsrjb2GgsmaOPkTkQTz+XiQPv8t6BAnr1GMYHSP4BoNhDNQzMLAJpKyxDSmXMFGhlfwJVXLUISUUOpBGpgQR/SKk2NrbSEnlIZMFGYZlXIKvlPpLpdRrSqlXlVKPKqXqrOVKKfVNpdRha/3lE2OuwWCYHjYjwtqGxNNbEM97Zd42ZYwvtJKbUKXL2ncXsBbpuG1mIO0yVyUzZv2bDtuxMt4snX/QWv85gFLqj4CvAp9GkmSXWf/rgf9rPRoMhhlJCBHWbUhVzApEiPNDKVHGF1oZrhbPZ4FTwJ8ho2krkbz7ZUij0IQR+7EzLsHXWnflvSxGeldAPrX/0Fpr4AWlVJlSqlZrfWo8xzMYDNNJiAFxzWXttCGefdR6fts4jzHcyNwg4kfGGZjoBEwJhXNn3DF8pdRfK6VOAr+DePggzfzJvM0ijND0K6XuVErtUUrtaWlpGa85BoNhSsh5/CWIx1/C5IZWhgv3RCfxeLOTUT18pdTjDExjk89XtNY7tNZfAb6ilPoykqN177kYoLV+EJmFmHXr1ulRNjcYDAVDvsc/2QwX7jElFM6VUQVfa/3eMe7r+8AvEcFvQHpWcgQxibIGg2FcTGQhtrnJeLN0luW93ITkagH8HPikla1zJdBp4vcGg8EwvYw3S+fvlFIrkKDacaRnBcTTvwmZyyzO+HtyDAaDwTBOxpul85ERlmvgD8ezb4PBYDBMLGakrcFgMMwRjOAbDAbDHEFJ9KUwUEq1IH0B+VQiBbBnAsbWycHYOjnMJFthZtk71bYu1FqPOgqtoAR/OJRSe7TW66bbjrFgbJ0cjK2Tw0yyFWaWvYVqqwnpGAwGwxzBCL7BYDDMEWaC4D843QacA8bWycHYOjnMJFthZtlbkLYWfAzfYDAYDBPDTPDwDQaDwTABGME3GAyGOULBC75S6k+UUlopVWm9LrjpE2fSVI9KqX9QSr1p2bNdKVWWt+7Llq0HlVLXT6OZOXt+Uym1XymVVUqtG7KuoGwFUErdYNlzWCn1pem2Jx+l1HeUUs1KqdfzllUopR5TSh2yHsun08YcSqn5SqknlVJvWJ//56zlBWevUsqtlHpRKbXXsvVr1vLFSqnd1nfhR0op13TbCoDWumD/kRLLjyCDsSqtZTcBuwAFXAnsLgA7S/Oe/xHw7QK29f2Aw3r+deDr1vNVwF5kdujFwBHAPs22rgRWAE8B6/KWF6KtdsuOCwCXZd+q6f688+x7J3A58Hresr8HvmQ9/1LuuzDd/0AtcLn1vAR4y/rMC85e67fts547gd3Wb/3HwMes5d8Gfn+6bdVaF7yH/wDwRQamToS86RO11i8AZUqp2mmxzkKPYarHArL1Ua11n/XyBQYKjG8Cfqi1TmmtjyGVTq+YDhtzaK0PaK0PDrOq4Gy1jn9Ya31Ua50GfojYWRBorZ8B2ocs3gR8z3r+PeBDU2nTSGitT2mtX7aedwMHkBnzCs5e67cds146rX8NvAeZABgKxFYo4JCOUmoT0KC13jtk1ZinT5xKxjvV4zTxu8gdCBS+rfkUoq2FaNNo1OiBeSpOAzXTacxwKKUWAZchnnNB2quUsiulXgWagceQO71onmNVMN+F8dbDHxdnmz4R+FMk/FAQTPZUjxPJaLZa23wF6ENmKps2xmKrYfLRWmulVEHlaCulfMBPgc9rrbuUUv3rCslerXUGWGP1h20HLpxei0ZmWgVfjzB9olLqYiQ2u9f6kIPAy0qpK5im6RNHsnUYpn2qx9FsVUrdCmwErtNWkJECtXUECnEKzUK0aTSalFK1WutTVqixeboNyqGUciJi/32t9c+sxQVrL4DWOqqUehLYgIRvHZaXXzDfhYIM6Wit92mtq7XWi7TWi5Bbosu11qcpwOkTZ9JUj0qpG5B+kd/QWsfzVv0c+JhSqkgptRhYBrw4HTaOgUK0NQwss7IzXMDHEDsLmZ8Dt1jPbwEK4o5KiZe3FTigtf7HvFUFZ69SqiqX6aaU8gDvQ/ocngQ2W5sVhK1AYWfp5P6BtxnI0lHAvyBxsn3kZW9Mo30/BV4HXgN+AdQXsK2HkVjzq9b/t/PWfcWy9SBwYwHYejPS2KeAJuCRQrXVsukmJKPkCBKSmnab8mz7b+AU0Gtd09uBAPAEcAh4HKiYbjstW69GOj5fy/ue3lSI9gKXAK9Ytr4OfNVafgHihBwGfgIUTbetWmtTWsFgMBjmCgUZ0jEYDAbDxGME32AwGOYIRvANBoNhjmAE32AwGOYIRvANBoNhjmAE32AwGOYIRvANBoNhjvD/A/ivcKZxlxxYAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "tsne = TSNE(\n",
    "    n_components=2, perplexity=15, random_state=42, init=\"random\", learning_rate=200\n",
    ")\n",
    "vis_dims2 = tsne.fit_transform(matrix)\n",
    "\n",
    "x = [x for x, y in vis_dims2]\n",
    "y = [y for x, y in vis_dims2]\n",
    "\n",
    "for category, color in enumerate([\"purple\", \"green\", \"red\", \"blue\",\"yellow\"]):\n",
    "    xs = np.array(x)[embedding_df.Cluster == category]\n",
    "    ys = np.array(y)[embedding_df.Cluster == category]\n",
    "    plt.scatter(xs, ys, color=color, alpha=0.3)\n",
    "\n",
    "    avg_x = xs.mean()\n",
    "    avg_y = ys.mean()\n",
    "\n",
    "    plt.scatter(avg_x, avg_y, marker=\"x\", color=color, s=100)\n",
    "plt.title(\"Clusters identified visualized in language 2d using t-SNE\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Cluster 0 Theme:  All of the transactions have to do with legal deposit services from the Agency for Legal Deposit Libraries (ALDL).\n",
      "ALDL, ALDL Charges\n",
      "   Agency for Legal Deposit Libaries, Agency services\n",
      "   ALDL, Legal Deposit Services\n",
      "   ALDL, ALDL Charges\n",
      "   ALDL, Legal deposit services\n",
      "   British Library, Legal Deposit Services\n",
      "   ALDL, Legal deposit materials April to July 22\n",
      "   ALDL, Jan20-Mar 20 charge from Agency for Legal Deposit Libraries\n",
      "   ALDL, ALDL Charges\n",
      "   ALDL, Legal Deposit Services\n",
      "   ----------------------------------------------------------------------------------------------------\n",
      "Cluster 1 Theme:  These transactions all involve the purchase of goods or services from various vendors.\n",
      "GARDINER AND THEOBALD LLP, GB Feasibility Study\n",
      "   BSI, Annual subscription\n",
      "   EX LIBRIS, IT equipment\n",
      "   Private Sale, Literary & Archival Items\n",
      "   XMA SCOTLAND LTD, 2 invoices one for the replacement of obsolete hardware and the other for a new laptop\n",
      "   JISC SERVICES LTD SUBSCRIPTION ACCOUNT, Annual Subscription\n",
      "   EDF Energy, Utilities\n",
      "   JISC Services Ltd, Managed router service charge annual subscription 01/04/22 to 31/03/23\n",
      "   RICK GEKOSKI, Papers 1970's to 2019 Alisdair Gray\n",
      "   Trustmarque, Miscrosoft Software Licenses\n",
      "   ----------------------------------------------------------------------------------------------------\n",
      "Cluster 2 Theme:  All of these transactions are related to ECG Facilities Service, and involve payments for facility management fees, inspections and maintenance, maintenance contracts, electrical and mechanical works, boiler replacements, USP batteries, gutter works, and cleaning of pigeon fouling.\n",
      "ECG FACILITIES SERVICE, This payment covers multiple invoices for facility management fees\n",
      "   ECG FACILITIES SERVICE, Facilities Management Charge\n",
      "   ECG FACILITIES SERVICE, Inspection and Maintenance of all Library properties\n",
      "   ECG Facilities Service, Facilities Management Charge\n",
      "   ECG FACILITIES SERVICE, Maintenance contract - October\n",
      "   ECG FACILITIES SERVICE, Electrical and mechanical works\n",
      "   ECG FACILITIES SERVICE, This payment covers multiple invoices for facility management fees\n",
      "   ECG FACILITIES SERVICE, CB Bolier Replacement (1),USP Batteries,Gutter Works & Cleaning of pigeon fouling\n",
      "   ECG Facilities Service, Facilities Management Charge\n",
      "   ECG Facilities Service, Facilities Management Charge\n",
      "   ----------------------------------------------------------------------------------------------------\n",
      "Cluster 3 Theme:  All of these transactions are related to Kelvin Hall and involve spending money on services or goods.\n",
      "CBRE, Kelvin Hall\n",
      "   GLASGOW CITY COUNCIL, Kelvin Hall\n",
      "   University Of Glasgow, Kelvin Hall\n",
      "   GLASGOW LIFE, Oct 20 to Dec 20 service charge - Kelvin Hall\n",
      "   Computacenter Uk, Kelvin Hall\n",
      "   XMA Scotland Ltd, Kelvin Hall\n",
      "   GLASGOW LIFE, Service Charges Kelvin Hall 01/07/19-30/09/19\n",
      "   Glasgow Life, Kelvin Hall Service Charges\n",
      "   Glasgow City Council, Kelvin Hall\n",
      "   GLASGOW LIFE, Quarterly service charge KH\n",
      "   ----------------------------------------------------------------------------------------------------\n",
      "Cluster 4 Theme:  These transactions all involve construction and refurbishment work.\n",
      "M & J Ballantyne Ltd, George IV Bridge Work\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   ARTHUR MCKAY BUILDING SERVICES, Causewayside Work\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   Morris & Spottiswood Ltd, George IV Bridge Work\n",
      "   ECG FACILITIES SERVICE, Causewayside IT Work\n",
      "   John Graham Construction Ltd, Causewayside Refurbishment\n",
      "   ----------------------------------------------------------------------------------------------------\n"
     ]
    }
   ],
   "source": [
    "# We'll read 10 transactions per cluster as we're expecting some variation\n",
    "transactions_per_cluster = 10\n",
    "\n",
    "for i in range(n_clusters):\n",
    "    print(f\"Cluster {i} Theme:\", end=\" \")\n",
    "\n",
    "    transactions = \"\\n\".join(\n",
    "        embedding_df[embedding_df.Cluster == i]\n",
    "        .combined.str.replace(\"Supplier: \", \"\")\n",
    "        .str.replace(\"Description: \", \":  \")\n",
    "        .str.replace(\"Value: \", \":  \")\n",
    "        .sample(transactions_per_cluster, random_state=42)\n",
    "        .values\n",
    "    )\n",
    "    response = openai.Completion.create(\n",
    "        model=COMPLETIONS_MODEL,\n",
    "        # We'll include a prompt to instruct the model what sort of description we're looking for\n",
    "        prompt=f'''We want to group these transactions into meaningful clusters so we can target the areas we are spending the most money. \n",
    "                What do the following transactions have in common?\\n\\nTransactions:\\n\"\"\"\\n{transactions}\\n\"\"\"\\n\\nTheme:''',\n",
    "        temperature=0,\n",
    "        max_tokens=100,\n",
    "        top_p=1,\n",
    "        frequency_penalty=0,\n",
    "        presence_penalty=0,\n",
    "    )\n",
    "    print(response[\"choices\"][0][\"text\"].replace(\"\\n\", \"\"))\n",
    "\n",
    "    sample_cluster_rows = embedding_df[embedding_df.Cluster == i].sample(transactions_per_cluster, random_state=42)\n",
    "    for j in range(transactions_per_cluster):\n",
    "        print(sample_cluster_rows.Supplier.values[j], end=\", \")\n",
    "        print(sample_cluster_rows.Description.values[j], end=\"\\n   \")\n",
    "        #print(str(sample_cluster_rows.Value).values[j])\n",
    "\n",
    "    print(\"-\" * 100)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Conclusion\n",
    "\n",
    "We now have five new clusters that we can use to describe our data. Looking at the visualisation some of our clusters have some overlap and we'll need some tuning to get to the right place, but already we can see that GPT-3 has made some effective inferences. In particular, it picked up that items including legal deposits were related to literature archival, which is true but the model was given no clues on. Very cool, and with some tuning we can create a base set of clusters that we can then use with a multiclass classifier to generalise to other transactional datasets we might use."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "openai",
   "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.9.9"
  },
  "vscode": {
   "interpreter": {
    "hash": "365536dcbde60510dc9073d6b991cd35db2d9bac356a11f5b64279a5e6708b97"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}