atom-predict/egnn_v2/TSNE.ipynb

251 lines
54 KiB
Plaintext

{
"cells": [
{
"cell_type": "code",
"execution_count": 162,
"id": "1b9d61dc-a2c3-4400-aa3e-ead3f11a30e5",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:15.644715Z",
"start_time": "2024-06-14T15:43:15.641850Z"
}
},
"outputs": [],
"source": [
"import json\n",
"import torch\n",
"import collections\n",
"import torch.nn as nn\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"from sklearn.manifold import TSNE\n",
"from core.model import PL_EGNN\n",
"from core.data import AtomDataset\n",
"from core.data import AtomDataset_v2\n",
"from sklearn.metrics import f1_score, confusion_matrix"
]
},
{
"cell_type": "code",
"execution_count": 163,
"id": "28e77574-9a3b-4bfc-8cfa-0f18391c2153",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:15.663381Z",
"start_time": "2024-06-14T15:43:15.645894Z"
}
},
"outputs": [],
"source": [
"ds = AtomDataset_v2('/home/gao/mouclear/cc/data/end-to-end-result/tsne/')\n",
"model = PL_EGNN()"
]
},
{
"cell_type": "code",
"execution_count": 164,
"id": "61e67118-eeef-46cb-9ed5-daf70d6a526b",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:15.700330Z",
"start_time": "2024-06-14T15:43:15.664079Z"
}
},
"outputs": [],
"source": [
"ckpt_path = '/home/gao/mouclear/cc/code/egnn/logs/vor/version_1/checkpoints/last.ckpt'\n",
"ckpt = torch.load(ckpt_path)['state_dict']\n",
"ckpt = collections.OrderedDict([(k.replace('model.model', 'model'), v) for k, v in ckpt.items()])"
]
},
{
"cell_type": "code",
"execution_count": 165,
"id": "9b2b1dc0-edf8-4b23-b594-f2899f623621",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:15.703921Z",
"start_time": "2024-06-14T15:43:15.701297Z"
}
},
"outputs": [],
"source": [
"model.load_state_dict(ckpt, strict=False);\n",
"model.model.embedding_out[1] = nn.Identity()\n",
"model.model.embedding_out[2] = nn.Identity()"
]
},
{
"cell_type": "code",
"execution_count": 166,
"id": "46cfc31e-9622-4114-a613-a310193b90ed",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:17.840744Z",
"start_time": "2024-06-14T15:43:15.704511Z"
}
},
"outputs": [],
"source": [
"model.eval()\n",
"\n",
"hs = []\n",
"ys = []\n",
"names = []\n",
"\n",
"for data in ds:\n",
" with torch.no_grad():\n",
" light = data.light\n",
" pos = data.pos\n",
" edge_index = data.edge_index\n",
" name = data.name\n",
" \n",
" y = [int(data.label[0])]\n",
" h = [model(light, pos, edge_index)[0].numpy().tolist()]\n",
" \n",
" ys += y\n",
" hs += h\n",
" names += [name]"
]
},
{
"cell_type": "code",
"execution_count": 167,
"id": "46396926-df6b-4933-95f3-0adbd1245ce6",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:17.848918Z",
"start_time": "2024-06-14T15:43:17.841416Z"
}
},
"outputs": [],
"source": [
"ys = np.array(ys)\n",
"hs = np.array(hs)\n",
"names = np.array(names)"
]
},
{
"cell_type": "code",
"execution_count": 168,
"id": "1b2531f1-99fd-4933-b031-4713214be4be",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:17.858473Z",
"start_time": "2024-06-14T15:43:17.850059Z"
}
},
"outputs": [],
"source": [
"# 10\n",
"# \n",
"# idx = [name.split('_')[0] == '10' for name in names]\n",
"# hs = hs[idx]\n",
"# ys = ys[idx]"
]
},
{
"cell_type": "markdown",
"id": "3b8f5872-3e30-4f6a-b73f-accc3bbc016b",
"metadata": {
"tags": []
},
"source": [
"# TSNE"
]
},
{
"cell_type": "code",
"execution_count": 169,
"id": "83b4fc09-b7ef-4262-af53-48e743449b50",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:22.258477Z",
"start_time": "2024-06-14T15:43:17.858979Z"
}
},
"outputs": [
{
"data": {
"text/plain": "(-91.90207061767578, 101.54176483154296, -76.76430778503418, 76.35824851989746)"
},
"execution_count": 169,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"text/plain": "<Figure size 640x480 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgMAAAGFCAYAAABg2vAPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuNSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/xnp5ZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACQk0lEQVR4nO3dc5gc2RoH4F9VY2xPMpnYtm3btm17443tbJxsbNu2Oc7Yjar7x+zOTTaDRqm7z/s897mb7qpzvlH31wffoViWZUEQBEEQhMWixQ6AIAiCIAhxkWSAIAiCICwcSQYIgiAIwsKRZIAgCIIgLBxJBgiCIAjCwpFkgCAIgiAsHEkGCIIgCMLCkWSAIAiCICwcSQYIgiAIwsKRZIAgCIIgLBxJBgiCIAjCwpFkgCAIgiAsHEkGCIIgCMLCkWSAIAiCICwcSQYIgiAIwsKRZIAgCIIgLBxJBgiCIAjCwpFkgCAIgiAsHEkGCIIgCMLCkWSAIAiCICycXOwACIIgfvb5xTec3XoJPwIj4OLpjDpdqyF3iRzpXp8Yl4jHl57j/eMvcHKzR+HKBZCzqL+AEROE6aNYlmXFDoIgCIJhGKwcsgnH1pz97bnshbNi9f0FUCgVAIDEhCQcWX4K53Zew9cX3367vnCVApi0axjcs7jxHjdBmAOSDBAEIQk7Zu7H9un70n2eoij0X9IN+xcfxY/vERm2RdEUvPw9sPbRQtg52nIdKkGYHZIMEAQhuqSEZDR36QatWstpu/0Xd0OrEY05bZMgzBFZQEgQhOg2TtjJeSIAAGe2XeK8TYIwRyQZIAhCVCzL4sS6c7y0HRUaw0u7BGFuSDJAEISoJjacA41Kw0vb3tk9eWmXIMwNSQYIghDNhV3XcP/ME97ab9S3Nm9tE4Q5IXUGCIIQzfZpe3lru0CFvKjZsbLB9yfGJeLh+WdIjEuCf8GsyFMyJ4fREYS0kGSAIAjRBH8O46Xd8o1LYeLu4al1CTLCMAzun3mCG4fuICkhGf4F/BAbFYdja84iOSE59TorWyUoioJMLoNXDk/YOdpAaa2El787Ok1uDU8/d16+FoIQAkkGCIIwO9/fBoKmqUyviwiOxKjq0/D9bVCm1yYnqFL/++Pjz788d3LDBVRoWhrT/x4Dmiazr4TpIXUGCIIQTY/8Q3V6IzZU4cr50WV6Gzy68BynN19EVGg08M8rnqO7A7QaLeKjEjjrr3zjUph1dDxn7RGEUEgyQBCEaK7sv4XZ7ZaIHQantn9YCZ8cXmKHQRB6IeNZBEEIjmEY3D31CFcP3ISjh6PY4XDqyMrTYodAEHojawYIghBUUkIypjadh0cXn4sdCi9+fA8XOwSC0BtJBgiCENSkRnPx9MpLscPgjYffryclfnn1Hc+uvQJYoELT0nDzcREpMoJIH1kzQBCEYMbXn40HZ/krMiQF658uwvtHn3Fpz3U8ufQCqiT1L8/nLpED885MgpO7k0gREsTvSDJAEIQg1o3ehgNLjosdBq8c3e2hUWmREJOY4XVKGwV2fVsHJ1cHgSIjiIyRZIAgCN5FhkajrXdvscOQFLlShnlnJqNYtcJih0IQJBkgCIJfWq0WA0uPw8cnX8QORZJcvJww78xk5CyaXexQCAtGkgGCIDihVqlxdutlfH0dAJ+cXqjTtRrsHG1xYMkxrBu9XezwJK9u9xoYvWkAKCrzyokEwTWSDBAEYRSWZbF50i7sW3QUjIZJfZyiKNTvVRN3Tz5EeGCkiBGajurtKmLS7hFih0FYIJIMEARhlDUjt+LvpSfEDsNs9F/cFS2GNSJnHBCCIskAQRAGiwiORLssfVPr/RPcsHe2RbeZ7dF0YD2SFBCCIL9lBEEY7MCS4yQR4EFcVAJWDd2Mxb3XgHxeI4RARgYIgtCZVqPF4ZWncGbLJUSGRiMqJFrskMzenBMTUbZBCbHDIMwcKUdMEIROvrz6hmGVJnN65C+RMVpG4+iaMyQZIHhHkgGCIDL17U0A+hUbA61GK3YoFoXRMvj87KvYYRAWgKwZIAgiQ7GRcRhSfiJJBERi62AjdgiEBSDJAEEQGZrfbQXio8nUgFiqt68kdgiEBSDJAEEQ6Qr5EoY7xx+KHYZFe3D+CRiGyfxCgjACSQYIgkjXnA5LxQ7BKNZ2VvDI6iZ2GEZ5evklFvZYJXYYhJkjWwsJgkhXPUVbMFrTeIlw9nBE/vJ5kJyggqefO+r3rIFClfLjzsmHmNJkntjhwcbBGomxSQbfP2nPcFRvS6YMCH6QZIAgiDSxLIu6srZih5EpVx9nTDswGgUr5Ev3mu55hyDgfbCAUf3O0c0eifHJUCepDbrf3sUOe4M2QKlUcBwZQZBpAoIg0iH50/MowDuHJ9Y8WJBhIgAAy2/PRda8vpk2aeNgzVV0v4kJj4NWbfiOjLjIeLT17IVz269wGBVBpCAjAwRBpKule3fERsSLHcYv5HIZ3P3c0LB3bTQZUBf2znY63ceyLB5dfI5z2y8jPioBOYv5o2Gf2rCxt0bAuyAorZXIVjALGtt2glbD34I9iqZAURQYreF9DF7RC80G1ecwKsLSkWSAIIh0zeuyHBd2XhM7jF8cjtoGO0db3tqfUH827p99wlv7/6rVuSou7roGltH/JVhhrcD+4I28fh8Iy0KmCQiCSFefhV04aadI9QLIlj+L0e345vHm/Q1wxIb+kCtkvPYhU8jg7OEAucKwIrDqJDVOrD/PcVSEJSPJAGERIkOicOfkQ9w78xjx0dIa9pYyCgAtN/5l4v3Dz1j3ZBGGr+2LrHl9ILeSw8rOCnlL5UKb0U3glsVFp3YqtyhndCyZ8fRzx+IrM/lNOljg7f0PUCcbtpgQAHbOPkBONCQ4Q6YJCLMWExGLuR2X4eG5p6kvnHKlHA1610L/xd2gtCIrszMyr/MyXNh1nZO2Rm0cgB8BEXBwtUeVVuXg6v3/BKCxfSckJ6gybWPl3XnIVzoXJ/FkJikhGe2z9OWt+qKtow0SYhKNamPuqUkoU684NwERFo0kA4RZCf4citNbLiLofQjsXe1xbttlJMalv7dbYaVA6XrF0HZMMxSulB8Mw4CmyYAZkJJItXLvyWmbFE0BLEDJKDQfXB91ulYHy7CY1HAuIkMzPg5ZYaXA8fi/BP35nNp8EUt6r+G8XZ9cXgj6GAIY+eqbq3h2rH24kJugCItGkgHC5KlVatw6eh8rh21GZFAUJ216ZfdAn3mdUa1tRU7aM0XbZ+7Djun7BemLpikwmSykq9utOsZsGSRIPD9bMXgjjq4+w2mby27MxphaM6AysObAv2QKGY7H/WXw2gOC+BdJBgijsCyLV3fe4f7px1CrNMhXJhcqNCkNmZyfBVgsy0KVpILSWgkAOLziFLbP3Ic4nra/dZ3eFl2mtuGlbanrmnsQgj6Gih0GAMDOyRa7v62Fjb04J/g9u/EKo6pPA8tBNUZKRuGseh/G15+FB2efGt3e0uuzUahixnUWCCIzJBkgDBYRHInpLRfh1e23kMlpgKKgVWvh5uuCaQfHoEC5PEa1Hx+TgDNbLuHc9iuICIoERVOICY+FOlkDhZUc2Qtnw7sHHzn6atK3N3D9L/PblqKRXSeoEjOfx+cFlVL0iGVYeOf0wrzTk5Alt484sfwjIjgSXXINNvp7IlfIcCp5D76+DkCvgsONjmv+2SkoWbuo0e0Qlo1MjhIGSYhNwMAy4/Dq9lsAgFbDpFZXiwyJxrg6M1PmRPUQGRKF5zde4+PTLwj99gP9io/GmhFb8f7RJ0QERyE8MBLqZA0AQJ2sESQRAID1Y3YI0o/U2DuJuIedBbz9PTDzyDhsfbNM9EQAAFy9XXAsdgeKVClgVDuyf7YtZsufBXW71zCqLYqi4MfBlk2CIBNNhF5YlsWygRtwYv25dBc/MVoGqiQV/l56AoOW90y9782993h88TkYhkWhSvlQtGpBUBSFawdvY/WILfjxPSK1DVpGG1WhjUtCJR1SU7dHdez547Bo/Qd9CsWx1WdQoUlp0WL4L5qm8fnFN6PaSE5QobF9Z3j4uSHgbZBRbZVrXNLkT2UkpIEkA4TOrh64hbkdl+pUqlWrYXBh1zUMWt4TYd/DMbP1Iry++x60jAbLsmAZFg5u9nDxcMLX1wG/3S+VRABI2YpoiShK/IHDe2ce49ubAPjlk8anX4ZhEBsRZ3Q7yQnJ+P4m0Oh2Bi/ndrcHYbnE/2snJE+r1WLnnAOY1XaJXjXbE2ITcWztWQwsNRZvH3wAkPIm/2/51djwuDQTAamp1s7ydhRotVrsmX9I7DAAANun7xM7hFQURYGWSeNlU2mjhJe/p9hhEGbCMj/yEJn6ERiB/YuO4uiaM9D8M0+vL61ai+UDN3AcmbAUVnK0HdNU7DAEd2TlaU5WznPh66vvYoeQiqIoVG5ZFlf33xY7FJQiiwYJDkkjxSUk5ePTL+hVaDj+XnrC4ETAHMgVMiy5OgtyueXlzKc3XxQ7hFQfn31FXJR0Ski3HdMcmc2gjN4yEMVrFOI1jntnHiPwYzCvfRCWgyQDxC+0Wi2mNp+PhGjjyqSaMjsnW7Qe2RiHIrcif5ncYocjisTY9Ks2Co4FVg3bLHYUqfKVzoVpB8ZAbvV7kii3kmP6wdGo160GOkxoyWscGpUGA0uNM4mpNkL6LO8jj4UL/BCMUxsv4OvrANjYW6Nyy3K/FAm6f/oxQj6HiRylcKq0Koecxfxh52SHHIWzoWi1gqQcMQDvHJ4I/iyNgkMAcPXAbYzdOhgURYkdCgCgUvOy2Be4Aed3XMWz6y/BqBmUa1wKdbtXh0yW8rckRAmX+OgE9C06Eh0nt8LtYw/w8clnMFoW9q526Dy5NZoMqAuFkpy/QWSOFB2yIPsWHsGG8Sm13Rktk7p9z7+QH/44PQkeWdywedIu7P5DGgvH+ERRFFbe/QN5Swlz6I2peXbtFUZWmyp2GL84Er0dtg7iVCA0RERIJNr59BU1hqx5fbHu8cLUip0EkR7yEUjiosKi8e7hR4R+Ne7T+qW9N7Bh3F8A+/9te//+/5cX39DRrz865xqEW0fvGx2z1Pnl98XBiM0kEchAkSoFeJ/z1te3N6Y1HO7k7ih2CPj+NhB/9lsndhiECSAjAxL19fV3rBm5DfdPP059rED5POg5pyOK1yic4b2v777D4RWn8OTKC9A0DZ+cnnhy+SXPEUscBVRrXQFD1vSBk6uD2NGYBI1ag3ZZ+iHmR4zYoQAAGvSqiZEbBogdhl7qydtmegAT32gZjePxf5HpAiJDZM2AxLAsi7WjtuHvpSd+e+71nXcYW2cmZhwam25VtkPLT2L18C2QyenUmgChX3/wGrMUKa0VcPZ0gouXE+p0qYYGvWuRoVI9yRVyFCyfB3dOPIAUPjI8vfZK7BD0Zutog7ioBFFjYLQMgj6GIhspW0xkgEwTSAjLshhVfVqaiUDK8wDLsFjSZy006t+3/L28/Rarh28BAL2KAxnC2Uv8IdCMrLw7Dzs/r8HKO/PQbHADkggYqFanKpJIBAAgPCACk5v8gePrziExzjR2u1RrI42CVbRMGgsvCekiyYCEHFhyDM90+PQTFRqNu6ce/fb44RUnQdHC/NH3/qMzPLO5C9KXvmzsrZGjcDaxwzALlVqURc6i/mKHAQBIik/G3ZOPsGzgenTLMwSfnn8VO6RMNR/aUK+KhXlK5eAlDt9c3ry0S5gPkgxIyL4FR3S6jpbRCPrw+4mAD88/Sy31yyf3LK6o0b4S/Av58d6XIVx9LO+4Yb4olArMPDoe4CnHpPT8xMqyLMAC0T9iMb7uLCQnJvMTGEeyF/LD5D0jUrfuZsTRzQEdJ7XiJQ5LnCok9EOSAZFFhkRh4/i/0MSxC6LCdFuoxWgZ2Dnb/fZ4YpwwhWJ+BERgarMFUCWpBelPHxRFoUb7SmKHYVYu7bqW7gmVxjL0mGRGyyAiOAqX9tzkOCLuVWlVHjs+rkLJ2kXSv4gCpuwfCVdvfhLZkC+WUzuEMAxZQCii4M+hGF55MiKCo/T6RE/LKFRsmsYCQgFXLT8490SwvgDA098NCqUCAe/SL79K0RRs7K3RqF8dASMzfweWHOet7eR4leE3U8D9M49Qv0cN7gLiiUdWN8w/OxV3TjzAskEbEfbTJ3W/fL4YtXkgClXIB61GC3tnO87LLzu6kR00RMZIMiCixb3XICI4Eqyea/2s7axh7/LryIAqWQ11GosKzUGr4Y3Qf0l3qFVq7J1/BAeXnkBc5P+PkaVoKuVIZBd7zDkxAe6+riJGa160Wi2iedxaqEo2YnTpn+kCU1KuUSnsalQKPwIjEBUSDWcvp19+X2VyGfKVyc1psu2V3QPZJTqlR0gHSQZE8v1tIB5ffG7QvQkxibi6/xaqt/v/cPiTyy8EWS8gtAa9aqH/ku4AUuavO09pjQ4TWiDwYwjunXqEj0+/ACxQtFpBVGtbAVY2VuIGbGZomgZFUYKU1jVEvIQOMNKHu69ruklr9sJ+nCYDTQfUk0wZZ0K6SDIgkncPPxl+MwUcXnHql2QgMdY0tlrphQI6TGjx28MyuQx+eX3hl9dXhKAsC0VR8M7hgaCPup9TYGVrhVYjG2PX7IM8RpYiIiSa9z6EVq5RSRz8k7upmQ3j/sKDc08w9cBo2DkatkaDMH9kAaFI0qoToDMW+PDk8y8PZTW3N0YKqNmhMnxyeokdicXrNLm1ztfaOFhj86ulaNK/Lo8R/Z/GmGkGiSpUKT/nbT48/wxja8+Q7AgPIT6SDIjk5a23Rt2fFJ+MCQ1m4+vr7wCAnEX9kbsEP3uUxVCpWVmM3NBf7DAIAHW7VUfRagV1urZx3zrw9HOHm48LZIrMt9MZy80M14e8vfeen3bvf8TR1ad5aZswfSQZEEF4YARObbpgdDv3zzxBr0IjcHHXNQBAq5GNjW5TbP4Fs2L1/fmY/vcYMv8vERRFYcG5qajapkKG18nkMjQZUC/1HiFGB2p2rMx7H0JTq/hbCLxyyGYcXnGKt/YJ00WSAYHFRMRiSPmJ0Kq13DTIAvO6rkDgx2C84ekTBR/+Xc+UNa8Paneuiom7huGUajc2Pv8TeUrmFDc44jcyuQxT9o7E4JW9UgoQ/bQejZbRkMlpTNw17JdpnV7zOvFaBtfazgqNzXAbqVwp563IEwCsGrZZp0qnhGUhCwgFtnniLoR9D+e0TZZh0S33EE7b5IOtgw0oGQ1rWyWKVS+EZoMboGD5vGKHReih2cD6KFqlAI6sOoPHl5+DoiiUrlsMzQbV/23dyvE1Z8FouZ+jpigKVrZKzD05CXZOvxffMlUsy2Lj+J3Yt1C3SqSGkslpHPzzOIpUKcBrP4RpIUcYC0iVrEYzxy7QcDUqIBFypRyLr0xHYmwyZrRciMT4pNTtaBQoyJVyTNozHJWalRU7VEIgl/bewNwOSzltk6YpZM3niyqtyqNxvzpwz+LGaftiO7zyFFYN3SxIXzb21jgas0OQvgjTQEYGBBQVEmV2iYCTuwOW3ZyDLLl9AAB7Atbjwl9X8fDCMzBaBoUq5kP9njVJBTQL8vnFNyzps5bzdlsMa4T+i7tx3q4UaLVa7Jl3SLD+EuOToEpWQ2mlEKxPQtrIyICAYiJi0cq9p9hhcCZ3iRz44/QkOHs4iR0KIQHhQZEYWW0KAt//fogWF+xd7LA3cEOGb2CJ8UkID4yErYM1b3X++fD+8ScMKDlW0D6zF/LD+qeLSUEiAgAZGRCUo6sDZAoZd4sHRWRlZ4UlV2fCxs5a7FAICUiMS0TX3IOhSjTirIFMxEXG48yWi2jSv94vj8fHJODoqjM4vOIkIkKiUg9VKlghLzpPbQOWYXF22yWEB0bCI6sb6navgZK1i4CmpbN+Wp0sfCnxzy++4cXNNyjMQ10DwvSQkQGBTW7yB+6ceCh2GEabfWICyjUoKXYYhETM6fAnLu/l/wTBHEWyYf2TxdBqtIgMjUbotx+YUHc2EvSswFm8ZmHMOT4BSmslT5HqJy4qHm19egueFJRvXAqzjo4XtE9CmqSTGlsAlmXhm8tb7DCM1n1We5IIEL+4euC2IP1EhUZjy+TdaOPdGx2y9sOwCpP0TgQA4PHF55jbcSn3ARrI3tkOtTpVBS0T9iU5MiRK0P4I6SLTBDz6/jYQB/88jqsHbiEpPhlypRwJMaZxhsDozQNRqm4xrB62ObVaYt7SuTFgSVf45DT9hIbgDsuyYLR6Hr1poMjQaOyed4iTQ7luHL6HyLBouEhkzUufBZ3x8tYbfH8bJNj30yu7pyD9ENJHpgl48uTyC0xsOAdajRZajTB/2P8lV8oxcddQ2NjbYHrLhUhO0H0+t2rr8piybxSP0RHmgmVZ1JW1FTsMg1jbWWHawdEoVaeYJBbSxUfH4++lJ3Fs3VlEBkfx3t/aRwuRq1h23vshpI+MDOhJo9bgyKrT2DFzP+KjElIfd3RzQIlahZG7RE4Uq1EI01sthCpJ+ENUFNYKZMnjgwFLuqFkraKpj1vbWeuVDMT99LURREYoioKtg41Bw/ViS4pPxoT6c1CpRVlM3jMCcoW4L4l2TnboMq0Nukxrg+iIGLR278VbX2XqFyeJAJGKjAzoIfBDMIZVmoSo0BixQ/lFvrK50ahvbVRoUjrdbX4NrNrrVeOg1cjG6L/IPPd0E9w78OcxrBu1XewwjNJubDP0ntdZ7DB+sW7MdhxYfIzzdr2ye2DbuxWQyfg/TIowDWQBoY5USSqMrjldcokAAIzZPBANetbKcL+/g6u9Xm0KdQQtYfpUSSqc2mj8wVtiO7LqNBLjpDW60W9hV72OkNbVj4AIkggQvyDJgI6u7L+FsG/cninAlWfXX2NI+QmoK2+LOnQb1Fe2x9Tm8xEbFZd6TcM+tXVur8OEFqkVBQkiM5f23MDXVwFih2G0pPhkvLn3QewwftN9ZjtkycPt36NWrcWffdeS3QREKpIM6OjWsftih5CuZf3X4/Xd96krrLUaLW4dvY82Xr0RFpiSwDQb3ABOHo4ZtmNtZ4Whq/ug55yOvMdMmI8LO6+JHQJntBppFgTLms8384v0dHLjBXT0H4AT689Bq5Xm100IhyQDOkqKTxY7BL1p1VoMKjUOAODi6YRlN2YjZzH/367z8HPD8ltzcTRmB5keIPQWHSa9qTNDUDSF3CVyiB1Gmnyye/DSrkalwdL+6zGy2jTEx5BFw5aM7CbQkTdPf4x8iwyJxvd3QciaxwdZcvtg7cOFeHPvPV7feQ+ZnEaJ2kWRleMhSMKyZMnjg49Pv4gdhtGKVi0IJ/eMR8/E0mxwAxxeeZq39l/ffouFPVZh+sExvPVBSBsZGdBR9faVxA7BYIeXn0z9b4qikL9sHjQf0gBNBtQjiQBhtAa9a4kdAie6TG0jdgjpyprXF86e/BVHYhgWNw7fReCHYN76IKSNJAM6KlwpP+xd7MQOwyCPLz8XOwTCjJWuWwzFqhUUOwyjZM3rg6IS/xqmHhjJex8Pzj3lvQ9CmkgyoCOaptFtejuxwzBIyJcfYodAmDGKojD/3FQUr15I7FD0RlEUKJrCwGU9JVGBMCNFKheET04v3tqnKMosTlQlDEOSAT00HVQPDfvqvkVPKlQJprf4kTAtMrkMCy9OR4NepjVl4J3DE7OPTUCZesXFDkUnk3YPBy3n52WbZVjkLZOLl7YJ6SPJgB5omsbwNX0xZd9IQNofIn5Bk+IihAAenn+KU5tMo/hQvrK58Oe1Wdj2bgXKNighdjg6y1cmN+aenMR9wxSQs5g/CpTLw33bhEkgyYCeVEkqLO2/HjChIs60zIQyF8Jk7Vt4RNQkuVrbivDJpdswet6SuVC4Un7JTw2kpVTtovAvmIXTNh1d7TF5zwiT/H4Q3CDJgB60Wi1WDduM2Ii4zC+WEFWSGgwjzsmJhGVgGAYPLzwTLUnOXz4P2o1thqAPITpdf/PYPZ4j4tfYbUNA0dy8cZdrVALrniyGXz5uEwzCtFh8nYHIkCjsnH0QN4/eg0atRZ5SOVG7U1XkKZUTPjk9QdM0Xt15h0/PvmDzpN2I+RErdsh6s3e2A02TvI/gD8uwqRUwxdBmZBNEBEXqfH1cpGkX2MlbKhcWXZyOGa0XGfWa5OLljNnHJnIYGWGqLDoZuHPyIaY2mw9G+/9PzXdPPMTdEw8BAI7uDgALxISbXgKQigLqda8udhSEmZPJZchZzB8fnwhffKj1yMao0qo83j/6pPM99s6muU34Z0WrFsSBkE14du0VPj79AjsnW5RpUAKaZDUW9lyNhzpsExyxvq8AkRKmwGKPMI4MiUI7374w9y/fwdUe6x4vgkdWN7FDIczcma2XsKjnakH7VFjJcTJxNwCAZVn0LjxCp0OTus9qj06TWvEdnmgS45MwvcUCPDz/LN1r6nWvjtGbBwkYFSFlFpEMfHn5De8ffcaXl9/g7ueGnEX8sWniTjy/9lrs0Hjlk9MLc05MIHOBhCAYhkHnHAMFPd3Twc0ef4dtSf33/bNPMKH+7AzvcXRzwPb3K2DnZPqjA5l5ffcdNk/ajZc33yA5UQUgZTtl95ntUKtTVZGjI6TErJOBy/tuYFn/DYiLihc7FMHV6VoVY7cOETsMwsKEB0eio1//X6be+OSdwxM7PqxK/ffzG6+xduTWdI8i9svni3lnJsMzm2meNWIMrUYLUICMbDUm0mC2awb2zj+MjRN2ih2G8Cigx6z26DjRfIdACely83bBphdLMKTCJMRF8p+E5/9pX/yV/bcwp8Ofae5ocPdzw6Rdw1C4UgHeY5IqmZwkAUT6zHJkIPhzKLrktKy5MJmcRothDdF+XAvJnrxGWJbHl55h74KjiIuKQ+4SOXB5303ERXCbIMw/NwWu3i6ICI7CuDozM7y2SJUCmHdmMpTWSgBAfHQ83j/+jMS4JHhlc4d3Dk/Y2NtwGh9BmAqzTAaWD9qAY2vOih2GYGwcbLDi9lz4F8gqdigEka47Jx9iStN5nG1BzFnMH2HfwvWq+1GseiFM+3sMNo3/C6c2XfxlOoOiKVRvXwk9ZrWHTw7+zgAgCCkyy2SgT9GR+Pz8m9hhCKZk7SKYf3aq2GEQRKYenn+KFYM24Ps7w4/KpaiUQ3UYA5MKT393hGZyeFeeUjnRZWobVGhS2qA+CMLUmGUyMKjsOLy9/1HsMAQzdFVvNBlQT+wwCEJnh1eewqqhm8UOI1PNBtVHg961cHHXdcSEx8I7uyfqdq9OtuoSZseskgGNWoOtU/fi4J/HoVFpxA5HEPYudtj5eQ1sHchcJ2Fa9i08gg3j/hI7DJ1RFJV69kLXaW3RaXIrUsufMBsmnQyEfAnD9b/vID46AU6ejtgxYz+iw2LEDksw1rZWWHV/PrLlJ3UECNN0atMFLBuwHlqN6Z2dMXxtXzTqW0fsMAiCEyaZDKhVaizqtQYXd14TOxTeeWf3QPEahXFh1zWok1NGOxTWCtTsUBnD1/aFXGG2u0MJC8EwDO6cfITbx+4jISYBPrm8sHvuIbHDypR7Flf89Xk12bdPmAWTTAamNpuPW8fuix0GryiaQv/F3dBiaENQFAVVshqhX39Aaa2AR1Y3MjxJmK0vr76jd6ERYoehk1X35iFvqVxih0EQRjO5j5Wfnn81+0QAAFoMbYiWwxql/ltppUDWPD4iRkQQwvAvkBV5SubE+8efRD0JURfJCSqxQzA5LMuCVT0C4lYB2lcASwG0D8AkA1Q0IPMD7AeBtqoodqgWxeTOtd0+ba/YIfCOllHoMbuD2GEQhGgG/Nk95dhtCQ+A0TIaWfOSBF0fLKsCG9EHiGwPqK8BzA+ADQO0TwH2DcAEA+p7QGR3MCHlwWhN+MRYE2NyycAXHU4kM2k0MOfERFjbWokdCUGIpkiVAph/dgr88vmKHUq6KjUvCxcvZ7HDMClszFxAfVXHiyOAsFJgEk6A1bwHy2r5Dc7CmdyagYFlxuHdA/OrIUDLadRsXwntxrVA9kJ+YodDEJIQFxWPlm49JHfUuJ2jDTY8/5PUG9ADy0SADa0IwIidI3QWwHE2aOtKnMVFpDC5NQNVWpYzm2RAYaVA88H10X58Czi6OYgdDkFIzoOzTySXCNAyGhtf/An3LCQR0EvyLRiVCAAAEwBE9QADJ8DtIGhFNk5CI0xwmqBx/7qQyUwu7N84uNpjT+B69F3YlSQCBJGOxLgksUP4zZgtg0giYABWdZvD1qKB8Npg4rdx2KZlM7l3VQcXe8w4MlbsMIxiY2+NhRemwdHFXuxQCELS/CRUUIuW0Ri9eSBqd64qdiimKYmHujCxc8AkXeG+XQtkcskAAJRrWAozj4yFXMnPLIe1nRWqtq3AS9sAsPz2XOQqlp239gnCXNw+IZ1txIuvzES97jXEDsMksawWYAP5aTxqDD/tWhiTWzPwrwpNymDn59XoXXikXkeYZqZIlYKYdnAUHN0ccLl5OexfdATvHn7ipG2KplCpeVlkL2h6CwTDgyLx7XUA7F3skKtYdlAUhec3XuHMlksI+x4OdZIG0eExSIxNgpuvK+r3qIFanavAyobsiiAM8/7xJ+z547DYYQBIOfq4cMV8YodhslhtCI+tR4GJWgHaeQiPfZg/k9tN8F+7/ziELVN2c1KchJbROBa3A0or5S+Pq1VqMAyDjn4DEBNu+L7XghXzYe7JibBztDU2VMF8ef0dM1stwteftnQqrRVQWCsQH5WQ4b3Wdlbo9UcnNOhVkyQFhN7G1puFR+eeih0GAGDx5RkoWrWg2GGYLCZuCxD3B7+dOG4AbVuN3z7MmMknAzERsRhQcix+BEaA4eCwk3KNSmLmkXEpBU/+48299xhaaVK6/cit5ADLIktuH5SpXxwRQVEI+x4OV29n1OpcFWUbljCpOubf3gSgT5FR0GqM29/r5uuC+eemwr9A1t+eiwmPxalNF3Bp93VE/YiFk7sDanaojMb96sDOyc6ofgnTFRMei1YePcUOAwBg62iDSi3KomjVQqjeriKpAWIAJrxTSjEhnlEeV0HJvHnvxxyZfDIAAKHffmBux2V4ceM1J+2NWN8PDXvXTvO5728DsWzgBjy5/AIsw4KiKeQrkxvD1/ZBrmI5OOlfKvoUG4nPz75x0pabrwu2vVvxywjB67vvMKb2DCTFJf92PUVRaDa4HuQKOS7uuo6khGQ4ezqh5fBGaDqgHjmbwcx9exOAngWGix3Gbxxc7TH35ETkL5tH7FBMBssmgQ0pDUCA0s3KaqBdN/Dfjxkyi2TgXx+ffsHLW29xZd9NPL703OB2bBxscCRq2y9vOBHBkXh4/ilYFihTvzic3B3BaBnI5KbzSV8fT66+wOjq0zltc/TmgakLsOJjEtAxW38kxCTq3Y53Tk9seLoY1rbWnMZHSEdsZBxaefSU5NkECisF/vq0Cq7eLmKHYhLY5BtgI3sI1h/lcQuUjGz91JdJ7iZIT44i2SCT03h1+61R7STGJqa28f1dELrmGYx2vn0xv+tKLOi2Em28emNohUlIipfeHmguhH4Nw5ga0zlvd2m/dVjcew2+vw3E+R1XDUoEACD4Yygm1J/DcXSElDi42KNSszKSPJtAnazG8kEbxQ7DdLAZry3inPazsP2ZCbMaGVg5ZBOOrDrNSVvFaxSCf6GsOLLyTLrXOLjYY0/gul8WHMbHJODUxgs4vfkiIoKj4O7riga9a6F+r5qwsTONT7IzWi3C9UN3eO3Dyd0R0T9ijGpjT8A6uPm4chRRxtQqNd49+Ai1SoPshfzg5O4oSL+W7Mur7xhSboIkCw9RNIWzmn1ih2ESWM0nsD/qCdeh7TDQjoOE689MmE0y8OLWGwyvNFnwfovXKIxcxbPDykaJghXzYdXQzQj6FAL89F2lKArZCmbFksszJF9tUK1So5FtJ0kOz/7XoBU90XxQA4Pvf3L5BQ4tP4ln116BoimUrlcMLYY2Qr7SuRAREoVTGy/g3cOPCPwQjKAPIUiKT1nbIJPLUK1dBQz8swdJCnj27NorjKw2Veww0rQ/ZCOcPZzEDsMkMOGdAfUDAEIcNkSB8npN1hXpyWySgf4lxuDDk8+i9C2T09BqmV8SgP+iaApVWpbDlH2jhAvMAJGh0Wjr3VvsMHTSZ34XtB3T1KB798w7hE0Td6X87P7ZHUJRAMsCcqUcGpUm0zYomkKT/nXRYUILUp6WJ2Hfw9ExW3+xw0jT6gfzkadETrHDMAms5jPY8LYAGwOjzyfQhcNs0HZt+e/HjJjFmgFVkkq0RABAyptJJikVy7C4evA2fgSECxOUgSJDosQOQWe+uT0Nuu/ZtVfYNHEXAKQmAkBKIgBAp0QASPmZHl19Bn2Ljcb3d0EGxUJkzMXLCQorhdhhpMnL30PsEEwGJc8Oyv0QoKwiTIcJ24Xpx4yYRTJwcfd1sUPQDQu8vvte7CjSxLIsNk/chX7FRosdis78DazkeGj5SdAy7oYQYyPiMLDUWDCMAJ94zEDghyAs7LEK4+rOxPSWC3B5301EBEemea1cIUf9HtIrAaywUsDRVdpTflJDybIAtEDTaky0MP2YEZMtR/yvmPBYHFlxSuwwdEZL9MTFE+vOYfe8Q2KHoTOZQgafnF4G3fv06gswWm5nxxLjktAtz2Bsf7+KzFX+hypJhfePPiEqLBqrhm5B6Ncfvzx/4/A9gAIqtyyHnnM6IGse31++h73mdcKtY/fxIyBC6NDT1Wa0YdNTFk/7VZh+KMNGDS2ZySYDqiQV1o7ahpMbzv8y1Ct1hSvlFzuE3zAMg/Vjd4gdhl6qtioPucKwX9+4qHiOo0kR/CkMx9efQ5N+dXlp39RoNVrsmvM3Di49jvjoTLaXscD1g3dw/eAdUDSFCk1KY9z2wbB1sIWdoy02vVyKkdWn4cMjbs4JMYaThyO6z2wndhimiRZm9w9o03lPkAqTXEDIsiwmN52Huyceih2KXhTWCpxM2CV2GL/5+joAvQoOFzsMnTm6O2DX1zWwsta/LCzvle0oYF/QRrh4WuYqc61Gi1d33iEhJgHH153D7WP3YegrjFwpR+0uVVGqTlGUqlMMGpUGbX36cBuwniiawqYXf8Ivn3SOVjYlbOIJsNEjBOiJAuX1AhRlsp93BWeS36nHl56bXCIAANkL/V6bXwoCTGjxW5EqBTDnxASDEgEAuLCT5/UlLDC12TysuMXzoSwSwbIsXt99j9vHH+D59Vd4/+iTwcWk/kuj0uD0pos4vekiAKB6u4pw9XFBRFDa6wuEwDIswgMjSTJgKOu6QKwfwHBT5jx9LMAmAZQ9z/2YD5NMBvbMN5257Z+9e/AJwV/C4C2xVchOHoYv6tF1G56xqrQsh4HLe8Ld17hhxvQWqnHp9Z33+PDkM3IVy857X2KJDInCnvlHcGrjecGKAl3eexNypfgvWSc3nkfxGoXFDsMkUZQCcD8ANqwJwIby2JMVQJnO6bBSIM3VbBmIj47HowuGnzsgtpWDpVfG1L+gYSMWCis5KjYtzXE0v5PJaUw9MNroRAAAGK0wc4kPzz8TpB8xHFx6HG19++DvpccFrw4oROKZmcD3wWKHYNIo2gW013XAhscpH6uGoCiTe3sTlcl9t85suWwS1fHS80Ai57P/i2VZ7P7jb4PutXOyxdUDtzmO6HdcrWpJiEvExT0CbEOlYPSxz1J1cdc1rB25LdO6GubMiVQd5ATtNAZwOwGAh0WFyefAaExn+lMKxB9z09OFXdfEDsEoGpUGDMOApqWRh53efBF75x8x6N6oUOPOFtCVq7dzptewLItLu6/j9vEHoGU0ClbIh1qdKsPOyQ5AyhbUIRUmQp2o5jlaACxQsEJe/vsRiFarxeU9N3Fk1WmjDwEzB7U6CVQ4xwLQijyA920w8duBhJ2A9p/dIrQ3wBgzAhMH/KgP1vMmKNqOk1jNncklAzHhsWKHYBQbB2vJJAIsy+Kv2QfEDiNTeUrlxJv7H5CvdC4kxiXi/aPPYFkWuYpnh52jLa7su4FFvdaknh0AABd2XsOKIRvRdVpbdJrcCpMazUXQhxBB4vXLnwVFqhQQpC++aTVazGq3BDcO3ZXkCYJCy1nUH1ValRM7DLND23UF7LqCZVUAqwEoG7Axk4HE/Ua0mgg2fj0oByF2L5g+k0oGVMlqJCeqxA7DKE0HCnh6VyaCP4Ui9MuPzC8U2e1jD3Dr6H24+jgjPjoByQkpvwNKawWKVS+Ee6cfp30jC2yfvg83j9zDe4H2p9s52WLGoTFmU3jo76UncPPw3ZR/WPDUAADkK5sLs49NgEIpzfLI5oCilAD1zymwjrMAeT6wcRsB1sBRgoQ9AEkGdGJSycD60dsRFWK6ZSZdvJzQZpR0KpepJbAYSxf/lsKICIr65XFVkjr9ROAn7x99AkVR4LukRpkGxTFm8yC4eDnz2g/XAj4E48zmi3h5+y00Ki1yFsmGDhNbws3XBfsXH+VszYYp6zqjLbpMaSN2GBaFomjAritg2wls/GYgbqH+jbCm+34hNJMpOhQTHot2vn2gUZvmwqziNQtj3LbBkjrdTpWsRgvX7lCZ+GiLFFRtUwFT9o4UOwy9fHn5DYt7r8Gr2+/SfP7nEx0tDS2j4erjglK1i6D5kIbIXSKH2CFZPCZuDRD3p553yUB7v+IlHnNjMiMDjy48k34iQAFggWrtKqLHrHb4+PQbaBmF3MVzSPKEs/joBCitlSQZMIKNgzV6zu6A5kMaih2KXj6/+IYh5Sf8ss7ivywxEaBoCi5ezlh6bZbBZ18QxmGZaCDxINiEv1OOPJb5A3Z9Qdn1Bxu3GkD6v7O/oVx4i9PcmEQyEPA+CCc2nBc7jF/YOtqgxdBGsHe2xePLL6BKVCFXsexo1K8OsubxAQBkye0rcpTpYxgGczsuQ1xUnNihmLTE2CRsm7YPxWoURo7C2cQOR2erhm1GcoIeL6oWwMnDEY371kGzIQ0stpy02FjNJ7Dh7QH2p+JgTDAQdQesojxAuQFsoD4tgmVVKWsRiAxJfprg6sFbmNtxGbQSGxXIVjAL2o5uhgfnnkCrYZC/bB7U614djm7SPtZUq9Fi97zD2LvgMJIELhhjzuycbLH9/UrJ//wBIPhzKLrkHCR2GJJiZaPE8fidYodh0ViWARtWA2Ayqg/gAkC/KqKU0xJQNo2Nis0SSCYZ0Gq0uHXsPj48/gyFlQL+BbNi+8z9+Pj4s9FtN+5fB2e3XeF0OJyWUSnH4P60aFxppcDkvSNRoQn/VfkMkZyUjOGVJ+P9w89ih2KWes/rjHZjm4kdRqYeX3qOMbVmiB2GpGQv4ocNT5aIHYZFY5OvgI3kuiohDSirgnZdz3G75kcS0wSX9tzA8kEbERcZB0pGgeXwrPnuszug08SWsHW0xb4FhhXXSQvzb4w/hapKUmNaiwVY93iR5IaMI0OjMajMOIR9Cxc7FLN1ed8Nk0gG7F1IEZb/6jW3o9ghWDw26SIPrTK/TjkQ6RK1+s2Xl98wsMw4zO24FHGRKXPXXCUCFE1h9rEJ6DSxJViWxbntVzhpNzMsw2L7tL2C9KWPWW0Wk0SAZ4kcndbHN7/85MS9n/nm9kL5RtIczbMoDB8F5WSALDsP7Zof0ZKBb28CMKDkWLx78JGX9lmGxcXdKaWL75x8iMjgKF76Scv1Q3fx6o50yra+uf8Bz66R7TV8oiggR1F/scPQiSqRLBz8l72LHdY+WiR2GAQAKIrw0KgWlC2pD6EL0ZKB8fVm81705uKu63h1+y2W9hN+vmhohUnYJpERgvtnHoOWcfejlitkoGTmUWGPKywLNOlfV+wwMhUTHoNO/gONbsdczl5Y82ABbOysxQ6DAADbdty3ad0MUJThvl0zJEoy8PjyC4R+5b8Mrkwuw46Z+xEeGMF7X2n5a9YBXNgp/sFKjIbhtDyuRq1FnpI5QZOEIFWjvrVRohYfn2y4NarGdE6OHX55SzojX4Zy9nKCd3ZPscMg/kHTtgDsuW2UYWDxdbR1JEoycGTFSUH60Wq0+PLyu973yRQyTvqnKGD3vEO8l8HNTN4yuTg/UvftvQ9gtCycvZzg6C797XRcyFcuN0rUKgL5T78fHn5uGLq6D4at6Sv58wiifsTg8/NvYoehF2s7K94OSIqLiuenYcJwDqO4bU91DGziYW7bNFOi7CYI+RomTEcU8CNA/0VzXL1xsizw5cU3RIZEwdVbvEpYpesVg5e/B8K+h4PRcltVLiokWpDT7GgZBZb955wCkXKrWh2qoMXQhmBZFhHBUaAowMXLWfJJwL9uHLordgh6y6hCotHIB0bpYfQpKKSj2MWAbUvu2zUzoowMuGdxFaYj9qctgHrexyWxyyjLZDJM/3sMbOytOV07kEqAF1VGy6LX3E4YvqYvuk5vA6WN8BXFyjcuCQCgKApuPi5w9XYxmUQAAJ5cfi52CJKSNa+P2CEQP2GZOCB+Ow8Nh4FlyChQZkRJBkrUkP7cKlfsnGzh5iN+fezcJXJg/ZNFKNuwuNihGGz3Hwdh62CNnXP+FuU8hYu7bwjeJ1cOrzyFSyYcPx+6z+ogdgjEz1S3ode5A/pQP+anXTMiyjSBW1bpnNzHt8otykIm52YNgjGSE5Oxbdpe3D72UOxQDBYfnYi5nZaL1v+Lm29E69sYCbGJWDeah09cJqxEzcKo1IysMpcUls8tr2oe2zYPoiQDWi2/WwqlQq6QYdCKXmKHAZZlMbPNYtw9+UjsUEyajI8pFgHsmLEfGp638ZoS7xyemHd2ithhEP+lyM9f2/IC/LVtJgRPBjRqDVYP2yp0t6Lot6S7JPYwP7r4jCQCHChZu6jYIWRKo9bg9vEHeHP3PeJjEpG3dE4cWX1a7LBElTIyx0KrYVC4cn5M/3sMaNo0EztzRslzgVWUAdT3OG+bjRoOVu4P2A0BLScVONMieDJw7eDtlBXoEtRxUkvsmvO30e0orOToMbsjmg+qz0FUxls5dLPYIZg0igJsnWxRt1s1sUP5TURwJMIDI+Hk4YhPT79gQfeViAknx1L/LHeJ7ChZuygqNC2D/GVzm9SiT0tDOf0B9kdt7htWP0j5X+LfYJSVQblsIr8H/yF4MsDlYUFcomgK3Wa0w/F15xDzw/Aa2e5ZXFGqbjH8+B6OOyceoHT94pDJxFsz8PzGa3x7FSB4vxRFiV5fgSs2DjaYe3IS7Jykc8DPx6dfsGH8X7h/5jHZIpcJjUqDnnPIQUSmgKU9+O9EdR1seGtQ7gf578uECD5W9umFNIueuPm6gKZptB1t3KlzPwIicGbrJRxZdRqTm8xDr4IjEPA+o/O5+XV87VlRMmC5Uo6t78Rb7GcohVIO39zecPV2Rp5SOdF7Xmdsf78SBctLp/zu+0efMLTiJDw895QkAjrQaLitrUHwKH6DMP1onoGJE6gvEyFoMqBKVkOr5557uUIGVwG25nWbkVIXu/WoxihS1cjFJuz/CxcFfQzB6JozkBgnzol2r+++F+UTujpZjZ75h8HURuLUKg28s3ug9cgmmHtyItqNbQYnd0exw/rFsgHroU5Wc15AylwVLJ9H7BAIHTCJV4D4lcJ1GLcQrJb/svimQtBk4Nm1l3rfw7JARBC/51E7uTuibIMS+PLqO5LiklCzQ2XO2ma0DH4EhOP8X8KfUZCclCzqqASjZU3yg+vD88+wfuwOtPPti6NrzogdDliWxctbb7B6+BZMbT4fr+++J4mAHtyzWM5WZlPAsolg4reBiZ4OJm4jGG0imKhxQHQf4YNJJFMF/xJ0zUBEUJR+N1DclQbOSPSPGLTz7ZvSJQ2wXL/OssCVfTcFPdUu7Hs4Vg3bLP4wstj9G4HRMlgxaCOs7axRt6s4iwcT4xIxtfkCPL5Iqgca6tPzr2KHQCAlqWVjFwIJmwH89CIbt0C8mNQvhKimbhIEHRnw8td9cQhFi/Mj4jwR+EdCrDDTBBHBkZjWcgE6+Q8wyVr0UrRi0AbRFkPOaven+SQCFNB1elvU6FCZn7LY6dB3alJMapUaX15+w5dX3wX5ICQkNnYRkLARvyQCoqIASiF2EJIh6MhAoYr5IJPT0OqwoMfJ3RFRodLcgmgInxz8H5UaFxWPEVWmIPhzmNms5JeCpPhk3Dn5EOUblRK036+vA3DvlPnUh1AoFegytQ1aefQQdJojb5lcgvVlKI1ag91/HMLhFacQE56ym8nW0QbVWlfA4NW9oVSa9psWyyQACZvEDuM/WFBWNcQOQjIEHRmQyWVoPrhBptfZOdugeI1CAkQknOI1C/Pex5GVpxH8KYzMJ/Pg/pnHgve5fJB5rXa2sbfG6hFbEBeVIGi/DXrVErQ/fWm1WsxquwTbp+9LTQQAICEmEac2X0RL1+64c9J0y4gDAJt0FNIZEQAAGqB9AWvhpm6lTvCthY0H1Mv0mvioRHx4/Jn/YARUsEI+3vs4seE8GEZKf3Dm4+cXaSFotVo8vaL/glspiwmPxaFlJwVNVt2zukrioLCM3Dh0FzePpF91LzlBhSlN5uHhhWcCRsUxrcS2lNNeoFy3gKKEP/1UqgRPBoI+huh0nUwhyrEJvMmaz5f3PiKD+d11Ycku77mBJ5dfCNbfuwcfwTJkqsdYIzcMEDuETO1dmHkhNpZlsXbUVv6D4YssO4+N6zKFogAUZQDrhqCcFoLyOAtKnoPHmEyP4MmAXKlbNb7cJfx5jkRYY2vNQCTPayCkth/enLAsMKvtEqhVwpx+dn6n8FtRTQ4FFK1WAErr398MKIrC4BW9UKZeceHj0lFifBJ2zjmIt/c+6HT9p6dfceXALZ6j4gdl0wwAT5VYZfkAq+pIOylQAjbtQHmcBe22E7TzUlA2zUBRVvzEYsIoVuCVZnsXHsHGcX/pdK2Dqz3iIuPNZjEcLaNRo0Nl9JnfmZehyy2Td2PP/MNkzQCPJu8diWptKvDeT5dcAxH8KYz3fkwFLafhm8sbtTtXhbOXExQKOYpWKwjv7J7QarU4svI0Lu25Aa1ag4IV86HHnI6wc7ARO+x0JcYlYnSN6Xj78KPe229HrOuHhn14qN/PMyZ2JRDPV1VSCnCYBkpZGKw2ApB5gZL5gaLteerP/AieDGyevBu75xp/GJCu5AoZNBLbWkTRFGYeHovyjUtz2m5kaDQGlByDyNBoMKQEK+dkChnajm7Ke537yNBotPXuzWsfkkT9c6bFT9MjNvbWmLx3JMo2KCFiYNzbMO4vHFhyzKDEnaIorHuyCDkKZ+MhMn4xPzoCmvs8ta4A5XkTFO3EU/vmTfBpgiJVBD5XWoL1cFmGxbQWCxHOcWVFF08nLL0+G/nK5Oa0XSIFy7BQWvO/4OjH93De+5AcCshbKic6T26NMg1KoELT0hi6ug/2BKw3u0RAlazGifXnDB7Bo2QUjq4y0WOpXfncIaMGG7cJLENO7TSE4Kv0hJ7D06g0Ol03eGVPbJm8B/ECbXtitAz2LTyCAUu6c9quZzZ3nb9mQj+MlkH5JvzXGnByd+C9D8lhgbf3P2L+2amwd5bO6ZB8CA+IQHy04a8zjIbBkyvCLWblEk3bgbFqCCSf5KeDhLVgEzaDtW4CymEEKBn/9V3MheAjA0DKdh8pyVcmN5oNbIA/r8yEQikXrPphRtuJDHX/zBO8e/CR83YJwMbeyqgXcV15ZvNAwQrSOSVRSKokldgh8E6RxoJHfdG0KC/d3LAfzHMHKiDpMNjwVmC1uu1eI0RKBkZtHiRGt+kasb4fACBHEX/MPTUJjq4pn8xonpMCPsqNXt57A6TYNj8S45IxusZ0bJmyh/e+es7tKFpJbrE4uTvAycP8d8S4+bggZ1F/o36+peoW4zAigcVMF6ATLcD8ABu7RIC+zIMoyUDp2kVhbS+drR2Obv8fli1eozB2f1+LyXtGoO3Y5sYfZ5wBF08nHF55Ck+uvOBsx0RsZJxJHw5kCnbNOYibR7kf1flZsWqFMOvIOFjbSefvhE8UTaHpwPqQyXjafiYhFEWh0+RWRtWROLftMhITkjiMShgsEweohTozRQskHSdrCHQk2ljTljfLBT2sJCNyxa8vQAqlAtXaVkSvuR0xdf8o3uYw3z74iNXDtmB0jenoUWAY3tzXbb9xRnxzepv8J0qKotBmdFOsebgAMw+PFTucNB1Ycoz3Pso1KoXDkdvQcWJL+OT04r0/sVA0hQLl8qDduGZihyKYqq0roN+irgb/rcZGxmN83dkcR8U/VvNe4B7VgFa8Y9xNiWjvxu4+rjgYtgkVm5cR9c3L1dsZzp7pb0Vx9nDCgvNT4fzv8CXHof47IhD4Lhijqk/F19cBRrVXv1dNQSrXtRjWEJteLkXLYY2Qo0g25CiSDblLZDe6XZqmYetog7ZjmiJ38RyI42COnqIozhPP59deC1L/QiaXocfsDtj+fiXOavfhROJO5C+Xh/d++WRlo0ydgnPL4ooeszpgwfmpsLIx71GQkC9h2DRxFwaUGot+JUYj5EsY5p+dCsrAX82XN98g+kcMt0HyjRFhNIPUGtCJ4HUG0sKyLBLjk/Dmznuc++sK3j/8hE/PhDmDfNTmgajfPfOTq5ITk3Fl3y08vvQcLMuiUMX8SIhJwKEVJ/HjewSAlC0/VtZKUDSFxFjDfulpGQ1bBxvYu9jBL78vStYqirrdq6euY8hM9I8YtPbsZVDfurKyUeJw1DbI/1MyOuB9EHrkG2bQmyRFUWDBwsHFHnNPTkT+silveLeO3cfUZvM5iTtbgSwpyRZHv/GHIrbA3ln4FxqtVos98w5jqwBrF/jQf3E3NB/SAKpkNaxtrUBJcPsv1+6cfIgZrRZCq9aC+U+y7uzpiKhQw97UB6/siWYDMz/8TSoY9XsgvKFwHcrygvY4zns3rDYQbMI+QP0YgBKUdTXAulmaRY9YJhJIOg1WG5ay28G6gSRqI0giGUhLRHAkVg7ZhGsH7/DWR/vxzdFrbiej2tBqtfj+NghatRYvb73BsoEbQNM0p1UArWyVmHV0PErULPLL46okFa7su4X7Zx+D0TLIXzYP8pTOiVHVpnHWd1pW3p2HfKXTPhb23PYrWNB9pV7t5SubG55+bihRswhqda4K258qxyUnJqOtTx8kxCQaFTMfhq3pg8b9xDv1bGHPVTi79bJo/RvqcORW2DmZ9/bBn4V++4HueYdArdJwvp4nT8mcWH2fm2RZCCwTDzZUwLoRzutAW/N7TDGbeBhs9ASk/HAZ/H/42BFwGAJKngtQlgRgDTZuBRC/GiwYsCwFCiw0rAxng+vDyr4t6uQpJdrhSZJNBn726OJTbBj7F949/GR0W765vVGqTlE0HVgf2Qv5pT6uSlbj+sHbuHboDqLDYpEtvy9aj2qKrHl8Mm2TYRg8v/4aY2rP4K3yn9JagS1vlsPTzx0AcGTVKawdtf23mgJyKxk0yfxWXDwe/1eGQ7rvH3/Cpom78OTyC2jVGihtrJC/bG74FciKK3tvIOZHLCiKQqm6xdBtRtvUUYD07F98DOvHbOf6yzBalZblMfXAKNH6//TsC/qVGGNSBxrlLZ0Tq+6azpsXF/guE7720ULkKpadl7b5wITVBbSfBeiJAuX1EhSV8aJUlmVTPtFrgwDaFVCWyfSe1HtVT8BGtEWmWR5lh3g2D+zwON1Lpj2ojOshWVHNNw5jqg2DnY2w5/OYRDLwr7tnH2NS/TkG3au0UWLMlkGo3rbib8+FfAnDmFoz0jxR0TuHJ9Y8XAD7/3ySYVkWh5afxPYZ+wQpVETRFNqNbY5eczti9YgtOLSMp6IdOhi9ZRDqdatu8P2qZDXkCpnOe6VZlsWuOX9j61RpDYsXrJAXy24Y9vvIlVvH7mN6q4UmU356xPp+aNjb9OrqG2NwufF4o+NhRPqSyWVoOawh+i7sykv7fGATj4ONHilIX5TnfVD0/7ersiwDsNEAZQ2KsgGTeAGImQiwP1WDpdxAOU4CZdM40/aZyGFA8hmkjAj8KiDeHke+5MGjcE+8iHTD6fr74ahM/6CziGQrVDrWGWpGhuo+37Ci8RDY2QhXTdakzgnOUzy7ztfSchq1O1eFXCFHseoFUbV1hd/muIGUYf4J9Wene7Ry8KdQtHDpjplHx6FC49KIDI3GqU0XsHfBYSRECzd0zTIsbh29B1dvJ1ETAQDYOfuAUcmA0kq/oiv/bsU6tu4MwgOkc0yzFBZvVWhSGvuDN2LfgiO4cfguEuOS4OrjgqYD6yF/mdwYXnWKYFU1M+Po7oCaHauIHYbg+By5YVkGURL4PdSLdSNA/RZIWMtzR0qAsgUAsEwC2PgNQMKu1Dd+FvYA0th2yIb/k6xoQNk0T7d1VhuUZiLAssCiZ2Wx/nVxACxY0Kjp8znDRAAAXK2SUcY9GDdDs+JKkB82316EITX4/h79n0klAy6eznD1dkZEcFSG15VvUhJT94+GQpn5m87904/x7U1gptdNbTYfA5Z0x4Zxf4lW7jfgfTBWD98qSt8/C/4UKkq/jfrUwY4Z+yCVsayQL9I4VdDR1QG953VG73mdf3vu7x9bsGniLhxZeQrJCeJV93N0s8fK23/A2ta8dwykpWi1Qnj38CMvv7cMw6ZOHZoKiqJAOY4EoywNRI8G2Ch+OlIUA0XJUxKBiM6A5iV+fePOuP4AGz0NsG6Y5hw+o/mBB2/6411UXtjKNajq/Q3OVsmIV8sx4V41nPz+7yf6lPUDjkrd/vZ+vu74Fzf0U32GUpldp3uNZVLTBADw9v4HDCo7Pt3nfXN7Y9vbFTq3t2LwRhxdfYaL0CTHPasbL4fe0DIaZ9R7OW83MzHhsehfcgzCAyJ+W5EtluZDGqDP/M6CHGBkrMjQaLx/+BHxMQlgtCyW9FmL5IRk3vqzdbSBTw4vNOxTG4371zHtErpGCPwQjO75hvI2QrD9/UqTrUPBsgygugU2agLABnPcuhMor7tg45YB8Wtg0OpNp9Wgbf4/rcVq3uNp4GOMvvgOH2L+/zevoLXI6xiBl1FuYEHhv3vQS7gFY3+tI5l21/RsK7yMSknuKLC40SkfvN2a6B+3AUzurzNv6VxYfnM27Jxtf3uuRM3C2PJ6mV7tCfIpX6SdUw1718KOT6swYn1/9J7XCTmLcbMgJWs+X07a0ZejmwOWXpuFQpXyi9J/Wo6sPIUpTedDq5XWMdlpcfF0Qpn6JVC9bSXU7FAZ2Qpk4a8vLyfs+LAKax8tRNOB9Sw2EQAA31zeGLlxAC9tF66c32QTAQCgKBqUVSVQbn+B+4HqaLCJF4H4Hfg3EdD7o68qZTcbo/6A+697YeONcWh/9AM+xfwaq5qR4UWUR5qJAAA8CvfChxgnpLeGVMtQeBnphpdRbqmPyWkGcvr39zm+mNQ0wb8KlM+HwxHb8O1tAO6degQnd0dUbVNBp2mB/8pdMieAC9wH+TMWKb8fAn+YPb7uHBxc7VGnazXYOdqi3djmePvgA64euIX46EQkRCfg1vH7SIxNStkOyTCQW8mhSc44Qeo5p4NAX8HvPLN5YMmVmfj84hsGlhkHdVLG83B8Y1ng4fmnuH3sASo1LytqLPrqNr0tJjeZx0vbg1f2/qXMt6X79vI7L+1O2DmMl3aFRsmzAW77wYa3BcDh33TMr0mY3iUtmEhc/Xgf0y4fx5e4gjrckF4HFCber4Yd1Y4DDAsZ/f83Aw1DQcPSmPygaur9MopBNe9vcHUS7udrctMEXAv7Ho6O2fqLHQavXH1csPjyjDS3SSYlJOP2sfv4ERABV29nVGhWBjtnH8Te+YfTbKvZ4AYYvLwnzxHrZkH3lTi/46oglQAzQstolK5fHHOOTRA1Dn18exOA948+4eCfJ/DmHvclYkdtHID6PWty3q4penjhGcbVmcl5u3KlDKeSpLXDxlhM8h0gsgunbbJsyucwQwrdftb2Qd2/KbBsykJAYxVxCcXYondQwev/69SuBWfFwqflfpkeoCgW8ytp0KqkcK8pJjkywCUHV/MvVRkRFInx9WZh27sVvx0EY21rhertKv3yWO8/OqF841LYMHYHPjz5DLCAfyE/9J7X6bfCR2JqPbIJLuy8BlYrbjLAaBmEiLSoUl9Bn0KwqOdqPL3yMvUxSkZBLpdB/c+IkMJKnvrfhvr+NvNFuZbi6KrTkMlpaDne/lm7U1VO25MC2qocGEUZQM3tQWCGVrz/4y71TzLBzTTXs0hPdLnSBF428XC3SkBYkg1Ck+wBsKD+GTqW0wyGFI1CyxLCblu2+GTA2tYK2QpmxVeehvGkIuRzGO6deozyjUvpdH3hSvlF30OfmZxF/TF26yDM66L7glG+mMLhUOFBkRhacRKi/1P6ltWyUGs1qNC0NIau7gNXb2fUk7czqq8wHhaumqo3995zngjIFDK0HiXMwjLB2TTjNBkwtNp1nNYHF77x83cdkmiHkEQ7+NrEoJHfW8SprZGotYa/oxU6l2iBIr4CVmn8h8UnAwDQZ34nTGli/lXRLuy6pnMyYCpqdaqKE+vP4/l1YQ4OSo+hZ1EIacWgjYgKiU73+VtH76PjxJZw93WFe1bX1DM3DKFRS39BpRBYlkVUGLd1AKxslZh2YDT8C/plfrEpUpQROwIAwJPE8QBeZnqdIbxt7VAxWza0LFAY5bP6gZbA+RwkGQBQvlFptBrZGAeX8H+ghZhCv/4QOwRe9J7XCSOrTQPDMKKV5g35EoaosGg4e4h/4EhaXtx8gxuHMz9HflSNaVhx6w94+XsanAxQFPVLqW9Ltm3aXk52LBWtXhAeWd2Qv2we1OlS1bzPdlBdFjsCwLoj3GSlwHUy4Gtnj3Nde8JGof9id76RZACASqXGs6v8ZIBSYmrFSXRVsEI+zD05EQt6rEJ4gOGfZo3Vwa8/fHN5YdiaPihatZBocfzsR0A49s4/gsMrT+l0vSpRjX7FR8Mnl+Hb1SiaQoNeZPGgWqVOdyGuPuQKGabuGwUnd8fMLzYHiTzv7vrJv4OJsWrFPxUCHUHZ9wPseqEAKPg5OuFbTBS42h8+vnI1SSYCgAnWGeBafGwCmth1wtv7H8UOhXf1e/J7epeYClcpAEeRF4NqVBp8fRWAUdWnY3LjP8Aw4p4XEPA+CP1LjtU5EfhZ0Ie0y3PrYvCKXnDP4pb5hWbuxc03xk+XUEC9njUtJxEAAO1rwbqiKKDjpcYAHADH+aC8boKy75NS/4CikMPZGVztC7dXKNEobz6j2+GLRScDYd/D0cq9JxiRV6MLJdCIF3gpunv6EfoUGYm68rZoZNMRn559FTukVHdOPsT2GftFjWFh91WICY8VrD8rWyXmnpyIJv3FO9ZZSh5deG58IyxQoHzGp3qaE5ZlATbjMsFcYlhgcc3icPK7DNq2xW+lh0Pj/z3Xw/iRgaX1GoKSwNqA9FhsMqDVajG2zkxoLWih09ape6DVmMfXO63lAkxqOBefX3yT7BG+e+cfhlolTlGkzy++4cXNN4J+b9TJGhSpqkthFsvAVXXTRT1W4/bx+5y0JXnq+0jrBEA+sCwF2roWsnp3BUVZ//Y8w7IIijN+8ScNYG7NOqiZM5fRbfHJYpOBOyce4rsOBxSZk5gfsXh4/qnYYRht/ZjtuHmYw61HPG0L1Kg0eHOX+4I+uuCjkFBmGC2DsG/muUjVEB5+3E2VzGy9GDERwo3yiEYboN/1MsOH3SlKBsp+ULrP/3H9CqKTjTu7w8nKCjd69kX7wkWNakcIFpsM3DpyzyT2hnPN1HcUJMYlYv+SY5y2WbRqQdAyfv4UVCKVS354/pko/Vrb/f4JyxIxDIO//+Rud5JapcG5bVc4a0+yKD1349D6rqX45zWfdgflsgGUonCaVwXExGDzowd6tv0rZytr7G3dHl72plGW22J3EyQnqSD4YQES4ORh2guRzu24yvmPrcXQhnh27RW3jf7DX4QtdgzD4PEl4ZOBPCVzwCMrWTgIALvnHkIQx1Upn994jVYjGnPapuRYVQJgAyBRt+vV91KPfskQnQ+wbQUKiYA8D2BVHfFqBkeePcGZ9++QqFGjoIcnOhQuivzuHjj85hUoijK4dkk+N3fsaNEG7rbCHTRkLIsdGchZNDt0XRTSbVY7ZC+cjdd4hCCTy1C2gfCVrbj0meNFgtkKZkWl5mUxZd9IKK253fJToEJeuPm4cNqmLmIj4hARFCV4v52nthG8TymKCovG9hncH/FtCSOZFKUEbFvrd48uF1nXAW3fHZT9AFDWtfEuPBJVNq/HlEvncf3bFzwICsTOp4/RcNd2rL53B+EJ8QYXAnKzscWR9p1NKhEALDgZqNejBmQ6DA23H98cnSe1Rs0OlTK9VurylckFpbUy8wslzN6Fu2IrFEWlHrpUuUU57A3cwOmxvpN2D+esLX3IlcIM+FEyCqBS9sEPX9sXFZtKo3Kc2E5vvsTLDqUSNdIe0jY3lP1QcP3W9DbOH+sf3MOqe7dx+NVLNNq9A9GqX9cD/LtscdGt64hOTgJjwKiAnUKBkx27QvmfM2BMgcVOE7h4OmHkxgFY2H0VQOG3VddypRzzTk9Cseopf4Bizf1yqYYZJDQVm5fF7j8OcdLW9L/HpB68FPo1DKuGbsbXV3ouYMqAVzYPztrSh52jLQqUz4NXt9/x1kelFmXh4umErHl9UbtLVcvaB5+JZ9e4L2AmU8pQq7P5HUyUFop2AmvTDkjcbXRbkclWGH67Fm6EvAdNfQAFQKvDm/ybH/qvrSrtkwWL6taHh51pVoe02GQAAOp0qQYvfw/sXXAY904/Bsuw8PT3QPPBDdB8SH0olP8fNi5arSAwS8RgObBu1HaUqFUU/gWyih2KQVTJamwcu4OTtpoOrIuKzVI+yQa8D8LAUmORYALnC+iq/fgWmNZ8AeftUjQF9yyumLJv5G8nYBIp+NhLPnrTQNg62HDerlRRjuPBJh4CkPHfpEpLI0plDU+bhN+e0zAUelxthFf/HA2szyf9lz/C0LdkGax/mPmupZ7FS6FtocLI62baFV4tOhkAUlaSF61aEBq1BmqVBta2Vmn+MRevURhyhcykD2DRqLXY/cffGL99qNihGGRSo7l4csW4T10UTcHZ0wmdJqfMS0aGRGFohYmcJwIyubgzcBWblkG/xV2xbtR27hqlUkYdZh4ZRxKBDBSrXgh3TjzktM3KLcpx2p4UsdpgsAn7UmoNMEBmiQCQUkGw0rHOqO7zBZOK30J2h//XBTgfmAPPIw0fnVPSNCZWroZV927/tsVQTtNokjcfJlepARcb80jSKFbMo95MzKGVJ7F66BaxwzAKLadxMnGXyb2Yrx25FQeXnjCuEQooXbc4RqzrC89sHrhz8gEmN5nHy6aSXMWzY+3Dhdw3rKePT79gVrslnNXUkCtlsLG3QaGK+dB8SAOUrF1U0lXVxHDz6D3OR2VOJO6C0kqaNe25wCadARs1EoAW+hQdCkywQ9XjnSGjGNgrVDhY6xCyO8RAw1DoeKk5HoV7Gvzn7WxljYf9BiFZo8G9wADEq1Xwd3SGr6MD7JVWkjhpkEsWPzKgjxaDG2LH9P2IjRCuXCbXGA0DVaIKNvbSz2a1Gi22TNmN/YuOgdEaX5Ws24x26PzPiEB4UCSmNlvA2+7SOhKZ381Z1B+r7s5DM6eunLSnUWkRGxGH28cf4PbxB2g6sB4Gr+hFEoKfPLvK7TZVpY3SvBMB9TuwUcORkgTo/gepZSjs/Vgg5b9ZGnFqJRY9K4sCTj+w/k1xxGusjIorKjkJPxIS4G5ri8rZ/I1qyxRY7G4CQ1VvV1HsEIyisJKbRGEYtUqN7vmGYu/8I5wkAgCwY8Z+/AhMOdVw19yDnLWblhtHMj8uWCh8lqA+uvoMzmy5xFv7pujLi2+ctlejvWm/5mSGTfh3Kkv3REDDUAhIsMeOd//fYaFlaZz9ngMrXpYyOhH4V7xKxUk7poAkA3qq26262CEYpVi1gibxKW5Bt5UI5rhoC8MwqW9cNzgsZ5yWZ9deI1QipXntne3glZ2/nQ07Zu43uDiLObKy43D7LpVyCqRZS76ElOkB3TAscDkoG9pebI4Y9a9v+gxoaNiUAW9PW+NX9XvZm+bOAEOQZEBP+crkho2D9D9Zp6dc41Jih5AptUqNK/tvcd8wC7x/9AkAoErkP+M35hhgLlEUhZZDG/HWfujXH5jeciGOrzuHT8+lc3KkGMK+h3NzWuE/2o5uBmtb03290Y1uicCZ79kx4HodVDveCf1v1MePpIyL+jTMk9foyKzl5js9818kGdATRVGo272G2GEYzMvfU+wQMhX0MZS30/a+vwnE2wcfBFntb+sonXUZzQbXR4ma/BWtuXnkHpYNWI++RUdhVPVp+BEQzltfUrZp4k4kxnGzM8XZ0wndZ7XjpC1JUxRHZm9FLAu8jHTDucAcCEq0z7TJnsVLYnj5irAyYqG09MdPuUWSAQNUbl5W7BAMYu9sh1J1i4kdRqZYlr+5/M8vvmFQmfGICjX+aNKMuPq6IFfx7Lz2oQ+ZXIY/zkyGR1ZX3vt6fuMVRlSdivjoeN77kpL46Hhc2XuTs7UoQ1f3/qXWiTlg1a/BxMwFEzkYTPQUsKp7gE1n/LyDID6WxvcPSkSH//+N/I8n5bHqVWlk9hbtYm2NDY2bY3LVGnC0skbvkqUNjtXFWjrJvBBIMmCAIlULcHo8qVB6z+tkEquS/fJxVxJYLL3ndgRNS+vPSyaTYdGlGbBz5rdmOqNlEfwpFKc3W9bCwrDvEZzVIclWMKtZlXdmWQZM9Ayw4U2BhK1A8lkgcS/YiE5A7HxA2QDBX5XYvdwTz+/awc1bAysbBq8f2eD+Kxdsfqvbh5juxUuhVs5cqf8eUMrwD24/t2MJpPVqZSJkMpnJbafqs6AzGvWtI3YYOqFpmvc3LD6Vb1wKdbpWFzuMNPnm8sbml0vRbHADWNlxs+I6PcfXncOPgPDUxYXx0fEI/RqG5ETjzoiXqpAvYZy1NXnPCMjkplULJEPxG4DEnWk/p32NqMDzOLvPBR2GhqJUtVjY2DGwtmWRp0giSheIRP/8mRRx0qSMLCy9fQPPQ1PW6kQnJaHk+lUGhUtTFAaXKW/QvaaKFB0ywp0TD7Cw52pEh/E75Gys0nWL4Y/Tk8UOQy+3jt/H1KbzxQ5Db9XbV8L4HUNMoqjT1QO3MKvtEt77cfRwBFgWMT9iAaRUjZMpZKAoCg6uDmjcvy7aj2tm0kPiD88/xaRGczkbGTgYthmObg6ctCU2llWBDSkFIP0kkGVTfi8y0vp8czyO8Er/gmQtYCVDbhdXnOncHYVWL0OS1rCfx6I69dGyQCGD7jVVJBkwEsuymNhgDu6ffSJ2KGlS2ihxMGwzrG35/RTINZZlsbDnKpzbdkXsUDJlZaNEtbYV0aBXTRSqlN9kRozG1pmJRxeeiR0GAMDZ0xFb3iyHvZPpbeXSqDXomK0/osJiOFn4mrd0Lqy6O4+DyKSBVd1LmQ7I6JpMkgENQ+H411wYfbdWms9TKi2cT31HZLOU4kB+Dg74Fhurd6yetnZY06gpSvj46n2vqSPTBEaiKAqzj09AmQbFxQ4lTfNOTzK5RABI+b6O2TwIQ1f3gZf///fI2znZwtnLScTIfkfJaLQd2wwFK+YzmUQAAF7efCN2CKmiQmMwqto0scPQy4ubbzC7/RK0dO+ByJBoznbAtB/fgpN2JIPNfBtvZn82cppFSfd0tupqGdjf+wH7ZxGpD+mbCMgoCv1KlcHVHn0sMhEAyMgAZxiGQddcgxDyRRqFZgBg8PJeaDa4vthhGI1lWfwIiADLMHDL4orvbwIxuNwEJCUk81ZO2BCe/h7oNLElGvSuZRJJQWO7TkgWoN6CPv76tMoktr8eWXUaK4dsAi2jwGi5+yXsu6AL2oxuyll7UsBqw8CGGX98+ocYJ9Q73f7XB7UsoNLC7eQ3aNysEF1d/zfy8lmyYmOTFrBVclgsygSRkQGO0DSNjS/+lMTmVBcvJ2x+s8wsEgEgZZTAI6sbPLN5QCaTwb+gH5bfmouiVQuKHdovQr+E4c9+67BlsvHnsAuhSJUCoGXSegk4uuaM2CFk6v3jT1g5dBMAcJYItB7VBNverTC7RAAAKJkHjD0Gh9EC1x95ppQfZNmU/wcgi1GBZliEt8qB6MreqY9Dx8+4TlZW2NWqncUnAgA5qIhTNE1DoZRDnawRpf/cJbJj2t9j4G0Cn6yMlaNwNiy+NANBH0NwdPUZHFhyTOyQUu3+4xA+vfiKxn3qoHT94pJdTNhyRGPJrXX5/iZI7BAydXTVaVAUxWkJ5q7T2pjE4WEGsx0AJKww+HaGAS6OkyFb7CMkFHAGq6BBx2vwo20OsP8WEPu5kJiOI3MrG5hf8mUokgxw6NrBO6IkAv4Fs6LlsEao36um5Pa2880npxe6zWyHk5vOIyEmUTLTBrePPsDtow+Qt1ROzD01CU7ujmKH9Jsy9Yqj+6z22DplDyga4LHWk8788ku/xsSza685rZCptFbAygTX9eiDdhwCRnUB0Lz85fF/86lnt+1wepcrgr8p4eyuQa1WkShfJwayf96hVk7MiuCvVlBABaebKWeWhLbLCVZGAXTKG7/yezwcb4ZAGZQAxlqG+OJuiCvhDlaZ9mti1WzZUSlbNn6+YBNEkgEOPb3yQrC+8pbOhRmHx8DRzdEkCgnxydrWCtMPjsGkxn9AnaTmta8s+XwQoMen1/ePP2Nm68VYfHkGj1EZrtOkViheozCOrDqN28fvIzGWm1K6hqrY1PCKcUKhZdzOBTbqU9siknja/TCYH10AzZ3UxxgtsGBoNlw+7AKZjIVWS4GWsbhx0hn5S8Zj7q6PCA1Q4NTOlCJvP6dgcSXdABkNsCxcT3yDy/lAsDRAMQBLAbavouFyJgCBgwpA4/bP+Q7/bFso6e2Lrc1bCfjVS5/5/wYK6N/jcflEURQ6TmqFlXf+gLuvm8UnAv8qUbMI1j1aCL98/KwEVijlGPBnd2x5uQxype45NKNl8PTqS7y5956XuLhQqGI+TNw5DHOOTxQ7FBQob/zhMnwrVoO7Mx7cs7ii42TLeFNiWQbQPP7lsb+WeOPyEWcAgFabkmQx//z/2ye2WDTcD8e3pyQCtIxG88H10X9pN8QXdU39xG9/7wdczgcCSEkEAIBiU5ZvyaOS4bP+TWpRIjBATf8cONC2A39fqIkiyQCH/Av68dq+q48LNr1cih6z2pvEanWh+eXLgml/j+G+YQroPrs9Wg5rBIqiULhyPr1ul8lp3DmRSQU1Cbh3+pHYIaBXoeEID4oUO4wMcbUwl6IpzDw6Hs4e0toqyxc28Qh+LjwU8k2GwxvdUz7Gp4HRUrh5xgkv76dUI7VzskXX6W3RemhjFJ/TELIoFcAwcLkQkO7sIMUAypBEeG17B8/Nb9Eg0BobmprZ1k2OkGSAQ9XaVuS1/WU3ZvP2yddc+BfIiqqtuS0jSsvoX1aN1+xUVa/7GZaFOpnf6QsuXNh1VewQ8O11IEZWm4qE2ESxQ0mXf/6sKFw5v9HtsAyLnbMPcBCRiUj89Wvdv9YTCXGZLK5lgYAPKUP8sRFxaOfbBwEfgjCwTHnYfIqFLEYNZUhSppu47J9FwuFpBJoXKUI+SKWDJAMcylc6F7IX5md0oMfsDvDObv67BLgwcfdwNOhVk7NtnoyGQfZC//+5Zs3trdf9rJbFpxdfuQmGR6FfpHHscOD7YOxbeETsMDI0Yn1/KPSYLkrPjUN3kRgv7joN4fy/psW390pQOvyBUhSgtP5/Iq5RadGn8EgU9vDEiMa1YfUlTq8IyjYsodf1loQkAxybdnAMHFwzP29bVx5+bhi3fQg6TmzJWZvmTiaTYeSGATgStR2tRzaBjb21wW1RNAU3X5dfKkw6GTCse+f4Q3x48tngOPgmtaH53XMP4dubALHDSFe2/Fmw7OYcePm7G93WtQO3OIjIBCiKpP7n2b2u0GgoyBUZb2FhWQpJib++TamTNfh72QlQT8PgteM9GCtap01ElIyS5K4eqSDJAMey5vHBxudL0GZ0E9g62oCiKCitFShdrxg2v1mGSi1SjtSU/bMnVib/9Ufw70plB1d7jN48EDs/r0HtzvoNSxMpbB1s0G9RV2z/sNLgNuQKGSb8NQzPr7/GgSXHcHjlKdAyCi4GlEReN2q7wXHw7d3Dj2KH8AuGYTCj9SKxw8hQnpI5sePjar0WlKbl6sHbHEUkcfYjU//z40tr3DjhjOrNo0DL0n4rp2UsvLImQ538+9vU4eUncWrTRdAaFnSybntiyRRrxsjWQh64erug74Ku6Lug62/PTTswGi9vvcW57VcQGRIFN19X1OteHa4+Lrh55B4SYxPhlz8LyjUqCbmC/Hi44OzhhAIV8uLVrbd63Ve2QQnU61EDSwesx/c3gUbH8fHZF6Pb0EfY93Cc234FIV9C4eThhLpdqyFr3rRfEN89+MBrLAUq5MGrW+/0uufLi+84uPQ4Wg1vzFNU3NCqjastolVLoMCDAGjaHoxND2hituDdU1tER8hhZc3AP18SPr+yTvl0/89iQppm4eisQXJS2p9Xk5PUSIhJSP23LjOCNdtXNv6LMGPkbALCItw99QiTGs3V6VoPfzdM2TMKbj7O6FNsNBKiEzK/SQfOXk7YH7SRk7YywrIsts/Yh52zDv5WJU+ulMPFywnOnk4oXr0Qanepiou7ruPAkmPQavh7U/rr0yrMbLsEb+/pl3TYOtrgYNhmySbG394EoGeB4Ua10WVaG3Sd1pabgCSOZSKQ9LUCWuYrDI065Y2+WtNI+OZIxrXjzggLVMLBWYP8JRPw6r4twkPSLhNcpEp+fH8bhMiQaJ37PhS51SRPxRQKmSYgLELZBiUwatPADId03XxdMGX/SOz8uAZ5SubA6FozOEsEAKBAuTyctZWRwytP4a+ZB9Isl6tRaRD2LRzvHnzE/iXH0K/4GOxdcITXRAAA7JzsMGXvyMwv/I+EmETcO/2Y+4A4EhsZb1wDFNB8cANugjEBFO2KuLjsqYkAAFw56oLdy7wgV7DIVSgRahWF6yec0k0EAGDY2n6opcf0KUVTJBHIhDTTbYLgQf0eNVCxWWmc33EVL26+xtdXAVBaK5GzaDbU6FAFJWoWTt12NKnxHwj6kM6RqQaq0d74k9syo1FrsHXKXt0uFmhMMHfJHLB3toO9sx0qNC2DW0fv6XV/2Ddp7HJIy8/Haxti4NIecHRz4Cga0/DitjaNRyl8fv3/sxkoik33rKFyjUvBv0BWdJ/ZDifWnUNiXOa7McgugsyRZICwKI6uDmg5rBFaDmuU7jWfnn/Fw3NPOe/b1duF8zb/69Xtd7/MpUpBxwn/3wkzeEVP3DnxAIxW95GIp1deoOnAenyEZjRrOyvYOdoi3oDvedcZ7dBiSEMeopK207vkSMlE05/pZ1kKlIwC+1N9D1pOo0m/Ohi8ojcAwMrGChueLcGgMuMQ/SM2wz67z2yf4fMESQYI4jfbpun4yVoPSmsFcpfMwXm7/5UksT3rDi52eHP/A769CcDRNWcQHqD/FsYr+2+h26xAxEbEIT46Hn75skim5sa60duREKd/gSSZQoYuU1rzEJG0sSyLD8+tkdmSP4UVg82v1sDazgrvHnyEg6s98pXJ/VvBIC9/D+wP2YTtM/bhr1kHfhvtomgKE3YMRe7i/P/tmTqygJAgAESGRuPMlkv48OQz7hx/oNPQo65omkLTgfUxaHlPztpMT/DnUHTJOYj3fvRCwfgpif+0kb2wHyb8NQw5i/ob2bDhYiPj0NanDzQq/XcTZC/khw3PlvAQlbTFxSSghfO/u6zSSwhYWNsyWPt4FbLk9tG57ZiIWOxdcAT3Tz8CCwpl6hVH8yEN4JHVzei4LQFJBgiLd2ztWawcsill6JqLN67/KFa9EOacmAArG2GOqW1k2xEqnk9vlAK5Uo5Vd+eJlhA8ufwCo2tON+jeSXuGo3pb/teQSM2cDn/i8t6bmVzFAmBxKnmvZHeRmCOym4CwaFun7MHygRv+P4fNYSJg62CDCTuHYd6ZyYIlAgCQr2xuwfoSk0alwZ9914rWP0UbXu+6XKOSHEYifVqtFlOazfsnEcjsj4xClcbJJBEQGEkGCIv16OIz7JxzkLf2S9cvjpodKgv+otZ+vOWUrn599z1CvoSJ0vfnF98MvvfVbekeac2H1cO34PaxB//8K/MkqsMIaa19sQQkGSAsjlarxYNzTzCnw5+89pOruDjD12XqFYO9s+XsqRYjGWBZFgeXHDP4/vgoI+sTmJDEuEQcW3NWjztY+OcJ4i0eIm0kGSAsyt1Tj9A5xyCMrzcb0WEZb0cyVoXGpXltPz0URWHeucmi9C0GMfbph3wJQ6ARdSh8cnlxGI00BX8Oxcohm9DGuzdYRp/5NwpabSZHGxOcI5MyhFlQJasREx4LhZUCcjkNW0fb37YhPbr4DFOaztNrj7uhZAoZchQRb6X7n33Wida3kFx9XOBfMKvg/Wo1aRXO0Y1MLjP7rW5vH3zAyGpToUpUp1kJMyMyBQOltWUVYpICkgwQJi0yJAp/zTqAU5svQv3TCnoHVzu0HNYYrUc1gbVtyuK9DeP+0vuFyVAOLtwdY62vkC9h+PD4s2j9C2nQsh6/JX1C8MzmDgcXO4PKEWcv7MdDRNLx5t47DKkwSc/RgP8rUz0WFBvHcVREZkgyQJisH4ERGFxuAiICI34rXRobEY9t0/Zix8z9oOU0fHN64+ur74LF5l8wi2B9/de5HVdE65trOYpmw7fXAdCofv0kTtEU+i3qiqqtK4gSl0KpQJMB9bDrj7/13oFSrW1FfoKSgCeXn2N0zRkG3p3yjRz8x3eAsuUuKEInJBkgTIpWq8Wbex+QGJuIA0uOITwgIsPrGS0DRssImggAQPdZHQTt72fhgfpX+ZOqlsMao1qb8ji54QKu/30HGrUGRaoUQOdpbWBrb5N5AzzqNLkVHl9+jpc39Tsau3DlAjxFJK6YiFhMavSHgXenJAIT1nyBhy8LWFXjLjBCJyQZIEzGqU0XsG3aXsm/2fkXzIrClfKL1r+bD/9nIAjFM5s7bOxt0GpEY7Qa0VjscH6htFZi0cXpGFltGl7feafzfX75dK+qZ0rObr2M5ESVztfbOmiQECsHTbMoUzMG3cYFI2fBJAAUKLvuvMVJpI3sJiBMwv7Fx7Ckz1rJJwIA0GSAuIfq1OpURfA+C1cpAIW1gtM27V3sUKx6QU7b5JpCqcDckxP1uifzCnym6e6phzpfK5Mz2HbrFf5+/QzHPj3FzO2fkaNAMiiKBuU0D5SiCI+REmkhyQAhedE/YrB54k6xw9BZwQp5Re3fJ6cXWgwV9jQ8a1sl9gWuR/nGpThrc9iavpDJpL/FjGVYXeropNoxcz9/wYjov+s6MqLV0OheqQB2LfPCy/sOUGmLgnboDcr9LCibFjxGSaSHJAOE5F3cdR0aI7ZyCS1PyZxih4D+S7qh85TWnH9aT8/9M0/Qq+AIvLyj3/x5ejz93VGpeRlO2uJTbGQc+hYfrdciwpgfsTi99ZJR2xOlqEC5PKBlur+lxEfLcWCNJ8a0yoklYyqDdhgDSp6NxwiJjJCDigjJWzF4A46u1qeCmbh2fV0rmZPSEmIT8fD8UyTGJeHy3hu4e/IRb33RMprTGg5jtgxC3W7VOWuPD+vHbMf+xYZVIvQvmBXzz02V1BqPx5ee4/CKU3h65QVUyRrY2FvDN5cX6nStjlqdKsMmg0WbQR9D0D3vEDAGbCmkZRSOxe2E0kqY5JX4HRkZICTvhZ6rtcUW9u2H2CGksnWwQeUW5VCnSzU4uTvy2heXiQBFUzix/hxn7fFBq9Xi2FrDk9QvL79jTM3piIngtxKmrrZP34cxtWbgxuG7iI2MR3JCMqJCo/Hy1lssG7AeXXMPxve3gene75PTCyPW9zeob0bLIjI4ysDICS6QZICQvG+vA8QOQS9OHvy+6RrKlIalWYZF8OdQscPIUHxUApLik41q49ubQLRy74lmzl1xZNUpjiLT3/2zTzJdyxAVGoMe+YdhQsM5eHzpeZrX1OxY2eDf/4gg6S8ONmckGSAkT/VTZUFTkCW3NLeOFSgv7sJGvVDSTar+ZW3H3bHUCTGJWDlkM2a1XcxZm/rYMnm3ztfeP/0YY2rNwIE0Dmo6ufECosNiDIrhxc03Bt1HcIMkA4T0CV9t1mB8D8Ubo2hVaW/T+xkFCvW61RA7jAwprZXw8ON2bcjVA7dx/8xjTtvMTMiXMLy9/0Hv+9aN3o53Dz/+8tjRVacNjkOdrDH4XsJ4JBkgJC3gfZDe5V7FQlGU4Fv6dBUdHovZPB/Z/C+KNj5788jmhvo9pZ0MAEDJWkU5b3P5oI2ct5mRQ8tPGnSfTC7D0dVnAABqlRqX997AtwzWFGQmTynxd+FYMlKBkJC0OycepowMCJwQ2Nhbg9EysHG0gZe/B+ydbGHraIuAd0H4+PTLb9dTNAWfHJ5oNri+sIHqQJWsRo+8Qww6VMcQhh5Q87P5Z6fCzsmOg2j4Vbd7dZzZeonTNkM+hyIhNhE29taCHMJ0cdd1g+7TarR4efMNnt94jRmtFiEqNNrgGDz9PVCyNik0JCaSDBCS9vn5V1ESgb1BG2BjZ/3bcxq1Btum7sXhladSF49RNIWKzcpg2Oo+sHeW3hvY9mn7BEsEuEBRgE9OT7HD0EmeUjkhV8qhUXE3xM0wLJo5dYV3Dk+0GNIQTQfVg1zB30t1YlySwfeyLIvx9WZDnaR7GeLf0MDU/aNA02SgWkykzgAhWQ8vPMO4OjOF7ZQCWo9ogn6LumZ4WWJcIl7cfAuNSoPcJbLDPYs06gqkpY13b6M+tYlh4NLuaDG0kdhh6GTVsM04vPIUb0lrgfJ58OfVWZDJ+anG2L/kGIOOvKZlNHIVy46PTz9DqzF8W+nm10vhl1e8Uz6JFCQVIwSn1WhxYMkxdM83FG19e6N/yTG49vdt/JyXsiyL1cM2Cx5bseqF0GN2+0yvs7G3Qem6xVC+cSlJJwIAEB+dIHYIels9fCuOrTONQlM9ZneAmzd/hYNe3X6HvsVHQ5XMz66a3MVzGHSfTCFDwIcgoxKBDhNbkkRAIkgyQAgqPDgSbbx7Y93o7Qh4F4TI4Gh8ePwZM1svRo/8Q6HVpuyF//D4M768FObYYYqi4JvbG+O2DcG805OhtFYK0q9QbOx/n+4wBcsHbkDQxxCxw8iUrYMNOk9rw2sfX19+x8QGc8Aw3BV2AlIKJ907bVhVyuFr+kKVaHiCIlfK0ZXn7xuhO5IMEIIaUm4CYiPi0nwu4F0wWnn0wIVd13Blv3Anu009MArb3q5A7S5VeZ2bFUvNjpXFDsEwLDC/2wqxo9BJzQ6VYWXLXd2BtDy5/ILzctIvb75FhIGV/1YM2WjUWomR6/uZ5d+bqSLJACGYJ1dfIuxbeIbXxEclYl7n5dgz77AwQQEo16ikYH2JoecfnTgtkCOkFzfe4M6JB2KHkSlbBxv0XdCF304o4NSmC5w2mV5iroukOOOqL1ZvX8mo+wlukWSAEMzxNdKcA9aoTadMryFsbK2w8cWfcJfI4Un6Wj54k9gh6KTpwHpoNojHraUsEPyJ2xLN3jnE2bXh5OEAhZIcSiQlJBkgBKMyZvsRj2LDpXFQDJ+8snlg99e1WH5rDmq0r8R55Tw+hX4JQ9j3jEeUpKJkbe6LEKWiABcvJ06bzFnUH1nz+nLapi5GrOsneJ9ExkgyQAimWPVCYoeQpqEVJiFaIifH8e3plVe4tOdGptM1XOFqHn1x7zWctMO3ghXz8rYFECxQrEbh1EW26dGoNbh64BZWDtmEVUM34+qBW9Co05/bH7NlINeRZqjliMao1LycoH0SmSN1BgjBaLVa1Fdkvm3PKEZUK5xzcgLK1je/9QMsy+L8X1dxfN05vBT4MJjNr5fi/I5rOLnhvFG1DqzsrHA89i8OI+PP4j5rcHbLJTAcVGJMi3sWV7Qf3wJNB9b7rULhi5tvMKnRXMRHJ4CiKVCgwDAM3LO4Ys6JiXDxdsaZLZfw+flXWNkoUaFpGVjZKLFu9DZ8ePJ7ZU0uObjZY9zWwSjXqBSv/RCGIckAIZhrB29jZhv+TmVT2igAUFAlGj4dMfPIWFRoUoa7oER2ZttlLOm9BoyW2y1puijXqCRmH5sAICUhiY9OQPDnUJzfcRX3zzyGRq1B8OcwaHVcs3E4cqtJlChOSkjGpEZz8fTKS15LabcZ3QR9F3SFKkmFqwdu48SGc3h+7XWa11I0BSsbK6iT1WAYJjWJEPL3Yun12ShUMZ9g/RH6IckAIYj46Hi09urNadnWn9k62WLT8yUI/RaOyY3/MGqVdJn6xTHr6Hj+hnsFcmX/TcxuJ8zhRP+Vo3A2LLw4Ld1THEO/hmFq8wV6Vb7b+Xk1PLN5cBQhv7QaLW4du4/Tmy/i8aXnSE7gZ71M9iJ+CHwXLO1jvimgQpPSmHForCBnLRCGIckAIYjDK09h1VB+Kgo2G1wfA5f2SK1trtVqMajMeINKrP6rWtsKmLxnJEcRiqOxQ2ckxxu3/UtfShsF+szrgga9a8LKJu31AqHffqBPkZFIiEnUuV2KpnA0ZjusbU2vgFLY93B0zNZf7DBEQdMU6vWogcErepldMS9zQxYQEoJ4du0VL+3aOdn+kggAgEwmw9qHC9FtRjsorQ3bvnRl3y3ERho+uiC2r6+/CZ4IAMDmV8vQfEiDdBMBAJjddoleiQCQchLiqY0XjQ1PFB5Z3VCgfF6xwxAULQP6LOiM3d/XYeSGASQRMAEkGSCEw8MIYfMhDdI97azzlNY4FLkNfvkNq33er8ToDFdhS9mTSy8E77N8k1LwymQYPzI0Cq/uvDOo/QNLjnFejlcorUY0FjsEQc0+NhFtRzeDK49nNhDcIskAIYgiVQrwspDqyKrT2DPvEMKDItN8XmmlwNqHC1L611PY13C09uyVbttS5urrynmbtTpVSfc5F29nDFnZO9M2Lu8xvMx06Ncf+GEi9Qb+q3KLsgYnpaaEoiiM2TIIZeqXEDsUQk8kGSAEUadLVdg62nDeblxkPDZN3IUOfv2wf/GxNK9RWiux5MpMjNyo/7xtfHQCuucdYhIH5vyM6+I3rUY2xvgdQ7Hg/DTkLOYPik4Z5lHaKNBiaEOsfbgAnn7umbaTGJ9kVBymusJJJpdh+e25cM/CfZImBbSMQpn6JXAgbBPqdqsudjiEAcgCQkIwT6++xMSGc5GcmMzbdqtx24egdueq6T4fHRmHwWXGIfijfmVd3bK4YvfXtSa1Grpb3iEIfB9sVBtWtlboPqs9Wv9nmFur0SIpPgnW9taQyXTfdXHt79uY2dqw7aU29tY4HLUt3WkhU6BRazC30zJcO3Bb7FD0QtEU2P/UTbC2t0L9njVRsWkZFCifF9Y8H9RE8IskA4Sgwr6H4/jaszi95SIigqI4b9/VxwV7vq/L9E1706Td2PPH33q3r7CSo3LL8ui3qCvcfKQ9H/rm3nsMLjdBr3s6TGgOz2yeoOU0chbJhnxlcnOSALEsi8eXnmPL5D14dfutQW2Yww6Pf717+BFL+63D2wcfxQ4lQxRFoXT9YshXOjfiYxIgk8uQp2RO5Czqj+yF/MQOj+AQSQYI0STEJWLFoI04v+Mqp+1ueLYY2Qtly/Ca2Mg4tHTrYXAftIzG3FOTUIrPWvQcaGDVXueDmBr1q4Pha/pyHkPgh2BMbb4AX158M6qdvz6vglc2cQ7W4cvHp58xuuYMo+pi8ClHkWxYfHkGHFzsxQ6F4JnpjrcRJs/W3gZjtw5G7hI5OG1327R9mV5zbvsVo/pgtAzG15uFg0uOYtnADdgw9i98ev7VqDb5ULxmYZ2vPbHuHMbXn230vP7P4qLiMar6NHx7HWB0WwHvjJvykKKcRbNj3eNFYoeRplYjGmP5rbkkEbAQJBkgRPXi5ht85eCN4mfX/76Dh+efZnjN/bNPjO+IBdaO3oHja89i36Ij6Ft0FAaVGYeEWP320PPlzb33eHA24+/Dfz04+wSTG/8BrgYMz2y5hPDASE7K3r68Iey5CkLxyOqGis3KgJZJ4+XY2dMJ654sQv/F3cg6AAsijd8+wiIlJSRjarP5Rp0lkJ4lfdbizokHUCWrERkS9dsbdGRwFOd9AsDbBx8xqMw4zt5MjbFlyh6D4nh65aXBtQAA4EdAODZN2ImO/v2xfsx2zr4Xl/be4KQdKRq6ug88/Nx+SwgoCrB1tMHiqzMyXBjLFZ+cntj+fgVyFvHnvS9CWuRiB0BYrku7r/M2VxryJQyTm8z7ZRV0wQp50XRQfVRvWxFe/u54/+gTL31/fxuEm0fvo1Iz8Q48YlkWjy48M/j+s1svoaABVfM+Pv2CUTWmISEmkfNDcALeBSEhNhG2DtxvURWbm48LVt+bj0PLT+LEhvOIDI6Cg6s96veogZYjGsPd1xWFK+ZHwQp5sW/RUQR/yng3jHtWNzQdUBeqJDW8snuiauvyqd+32Mg4LOmzBtcP3f1lV0+ZBiUwZssg2Nib3/eXyBxZQEiIZmHPVTi79bLg/do52qBi87JGrxvISP6yubHi9h+8tZ+ZqwduYVbbJQbf//OJg7rSarXolmcIQr/++G0bGlfmn52MkrWL8dK2lLAsm+4ujn9PgExKSMbxtedwduslRIZGg6YpeGZzR4thjdCwdy3IFRl/1gv5EoanV16CZVkUqpQPWXL78PGlECaCjAwQoqH4qE+sg/iYRJzbfgUyBQ2tmp/ytmEiVso7v/Ma5ndZblQbWfPo/8Zw//RjhHwOM6rfzOyc/bdFJAMZbeekKAr2znawd7ZD95nt0H1mO4P68PL3QJ2u1QwNkTAzZM0AIRp9VrrzQatmeFu0ZSXSwqtXd94ZnQgAQLMhDXS6TqPW4MnlF7h+6A4u7zO81LCunl59iUATqwZJEKaAJAOEaKq2Lg+FlXiDUxRNIW+ZXCjXuCRo+v+fxIpWLYDq7SrCztnW4LbtnWwR/Fm/KodcWDF4g9Ft1O1eHT45vDK97sT6c+jg1x+ja07HjFaLOK8XkZ5Dy44L0g9BWBKyZoAQ1Zmtl7Co52rR+ndwtcffP7YgKSEZ0WExsHexg51jShIQFxWPNt69oVEZfnJhjQ6VMHHncI6izVjo1zB0yj7QqDZqd6mKsVsHpztMHR4UiTNbLuHagVt4//izUX0Zyi+fLza/WiZK3wRhrsiaAUJU9brXQPCnUPw164Ao/cdFxSMpIRnWtlaw9v/1+N2nV18alQgAwKXdN6C0VmL0JuPepHXx7W2gQfe5eDuhcKX8GLamD5zcndK97vTmi1jafx0YhuVtgaAuKNp0zocgCFNBpgkI0XWb0Q5LrsxEyTpFIVPofugNF1iGxeU9ae9fNzYR+NeZLZeQEMdvISKWZbFl0h6976vQtAz2BW7E1P2jf0kEtBotAj8EI/hzKBiGwcPzT7G4zxpoNYyoiQAooHgNcdeaEIQ5ItMEhOT0LzEaH558EaYzCqjfoyZGbRzw21OBH4LRLc8QTrrxze2FLa+X633iXlRYNJ5cegGtRou8pXMha17fNK+7tPcG5nZYqlfbShsFjsbs+OXUQY1ag51zDuLvP0+kFmpycneAjYNNpnvbhUBRFDY8Wwz/guSQHILgEpkmICRn+uGx6JJjkDCdscDnl1/BMMxvb9S+ubxRqk5RPDinX0nftAS+D8G+hUfRflzz3577ERCOs9uuIOhDMOxd7FG5VTl8fxuETeP/QmRI9C/XFqlSAOP/GgpPP3eEfQ/HvdOPERUShR0GTLNMPTD6l0RAq9FidK0ZeHH99S/XRf+IRfSPWL3b58PIDf1JIkAQPCAjA4QkNXXsgsQ47g7MyYyVjRIdJ7dC47514OjmkPr4nVMPMbkRN8WDHFztsTdwPRRKRepj+xYewaYJO1PqziJluD+zYXhnLycUq1YQV/ffNrjUb9mGJTDr6PhfEqBDK05i9bAtBrXHt3KNSqLbjHbIUzKn2KEQhFkiawYISarUoqyg/SUnqrBl0m50zzsE7x7+/4z5oA/c7WmPjYjD2/v/b/vcjivYMO4vMAwLRsuA0eo2Hx8VEo0r+28ZnAjU614D0w6O+W0kZN+CIwa1x7dV9+dj9rEJJBEgCB6RZICQpN7zOonSb1x0AibUn4OkhGQAGVeCM4Q6WQ0gZQTgr5n7YXARRgPyAIWVHAsvTMPozQOhtFL89nxEUKSBwfCDooEqrcojL0kCCIJ3JBkgJMnNxxVtRjcVvF+WYRH9IyZ1h0HxGoU4a1smp+FfKGW+++vrAAR+CDHoTd1QdXvWzHAlvpTmC2kZDaW1FbpObyt2KARhEUgyQEhW3wVd0Ht+Z8iVAq9zpYBHF1NO/PMv6IcStYsY36SMQtXWFeDimbJ9L/mfkQchXdh+BTeP3INWo03zeWfP9GsMCM2/YFYsuTID2QuRxYIEIQSygJCQPLVKjftnnmBJn7WICo3O/AYOFK1WEIsvzQAARIZGY3SNafj6KsCgtiiagk9OLyy7MRvOHilvuHFR8Wjh2p2rcPWSNa8P5p6chBc33+D0losI+/oDtIxGTHgsYsL5OVI6My2GNUSRKgWRFJcEv/y+yFcmN+dTNARBpI8kA4TJuLDzGuZxcAiPLhRWcuwL2gh7ZzsAQHJiMi7tvoGz2y4jIjgSXv4eeHnnHZJiM97xYO9ihxZDGqLFsIZwcLFPfTzoYwi65h7M69eQEYW1AuokNSiKMnghIldyFMmG9U8WixoDQVg6UmeAMBkl6xRNWXAnwHuXOlmD7TP2YeCfPQAAVjZWqN+zJur3rJl6zdUDtzCr7ZJ02yhQPg+W3ZiT5ifcgPfB3AetB3XS/xcyiq1M/RJih0AQFo+sGSAkQ5WsRmRoNFT/rLj/r3unHgm6yu30povpxgIAVVtXwISdw2Bl95/jiqmUffzzz01Nd6jbxl6cI46lqMVQ3Y5LJgiCP2RkgBDd8+uvsKjXGgS8CwKQUn+nQtMyGPBnd3hn90y9TuhFd4lxSYgIivwlhv+q2aEyKrUoi80Td+LKvluIDIkGo2Vw9+QjNHXoAgBw8nREi8EN0WRgXTi6phQ0Or7+vCBfg9R5ZfeAexY3scMgCItH1gwQovpr9gFsm7o3zeeU1gqsfbQQfvmyAACeXXuFkdWmChke9odsTF30l5bbJx/gz77rEBGY+R59N18XDFnZC7mK50DPAsNTaw5YsoUXpqJ4DeN3axAEYRySDBCiefvgAwaVGZ/hNT8vLmNZFm28eyM6LIb32GiaQr5yebD8xpw0n3986RlmtlmC2AhxVt+bg4rNSmPiruGwsiFTJgQhNrJmgBDNwu6rMr3m07Ov+PLqO4CUaoBT9o3kOywAAMOwUFopcHjlKcRHxwNIOVBox8z9GFVjGsbUmmnaiYAEdu3dPHIfrb164fLetI+QJghCOGTNACEKrUaLzy+/6XTtlxff4F8gK6LDY3Bu+xXIreTQJGt4jY+iKTy79gpPrrzApgk7kb9sbjy+9ILXPoXk6eeO0K8/xA4DSXHJmNNhKeycbMmuAoIQEZkmIEQRFRaNNl69dbqWooAsuX3w/Z8FhoRxbOytsf3jKtz4+w6OrzuHb28CkJygEjUm7+ye2PEx85EigiD4QZIBQlBarRYPzj7FpT3XcX7HVbHDMTsURaH9uObYPe9Qutf0X9wNrUY0/uWx5MRkfH0VgFd33uLvZScR8Fb4xGvbuxXwzeUteL8EQZBpAkJAD88/xcIeq/AjIELsUMxS1rw+GPBnD5RtUAI5i/ljxeBNiAmPTa0yaGNvjW4z2qHl8Ea/3WtlY4U8JXMiT8mciA6LxY6Z+3U6TplLAe+CSDJAECIhIwOEIF7efouRVadCq9VK63g8M+Hk7oD9IZt+KXKkVqlx7/RjhH79AWcPR5RrXAo2dtaZtvUjMAKdsg8Ao2H4DPk3K+7MRf4yeQTtkyCIFGRkgBDE1il7wDAMSQR40mpUk9+qHSqUClRsWibN65MSknF+xxWc3HgBMT9ikTWPD+r1rAH3LG6gKKDpwHo4vPyUEKEDSJneyFk0u2D9EQTxK5IMELy7few+Hl14JnYYZu3AkmMoUrkAClfK/8vjQR9D8PDCM7x/9BGsloWHnzvylcmJWW3/REJsYup1IV/C8OD809R/0zJhdx1Xb1cRSiuFoH0SBPF/ZJqA4I0qWYWhFSfhw6PPYodiERTWCqy4NRe5imXHzWP3sbTfOkQGR4kdVqYomsKR6O06TWEQBMEPMjJA8GZQmfH4/Fy3WgKE8dRJamyetBt1ulTFnA5LxQ5HZzU7VCaJAEGIjCQDBC9WDd9sFolA8RqF8eTKC8FX1hvq7smHeHzRtKZkNGp+C0gRBJE5Uo6Y4NytY/cFXXzGp8eXnsO/kB9qdqwMr+weYoejE1WS6RyAJJPLUk9yJAhCPGRkgODcgSXHxA6BU5+ffcXnZ1+hUMph72KHuMh4sUMyG1qNFjU7VRE7DIKweGRkgOAUy7J4fv212GHwQq3SGJUI1O1eg8NozIOHnxs0ag3IOmaCEBdJBgiT4p3DNIbqf0bRFCo1L4sR6/qidL1iYocjKWHfwjGm5gyMqj4t9XRIgiCER5IBglMURaFw5fygaH7OyA3+FMZLu3xSWivRbWY7yBVyzDg0FjU6VBJ8H7/UPbv2CtNbLRI7DIKwWOQVieBcqxGNTWb1Pd/88mfBkiszkKNwNjy69AytPXvh0u4bYLTClvo1BY8vPsfHp1/EDoMgLBJJBgjOVWxaBp2ntNb5ems7KxSukh9ypfmsZ/UrkAWzj4/HxudLkLdULjy9+gLjas9EYlyS2KFJ2unNF8UOgSAsEqlASPDm+Y3XOLLqFJ5dfY2E2ATQNI3E+KTUA3CU1go07lcHPed2hJWNFYZVmoSXt96KHDW3suTxQdGqBXFq0wWxQzEJpesWwx+nJ4sdBkFYHPP5KEZITuFK+X+rlZ+cmIyPT7+CZVnkKOwHG3ub1OcqNS9rdslAwLsgBLwLEjsMvciVcmhU4hQCcvF2FqVfgrB0ZGSAkIyY8Fi08e5N5tMFZOto808xJQrZC/mhdueqeHPvHbZP3y9KPEuuzkSRygVE6ZsgLBkZGSAkw9HNAZ2ntMb26fvEDsUi0DIas49NQJEqv775ilXO2N7VniQCBCESsoCQkJR2Y5vBO4en2GGYPb98vvjz6szfEgEAvxxlLCTPbO6i9EsQBBkZICRGaa2EVw5PBH8KFTuUVBRNmfxWSVtHW5RvXBKFKuVD6brF4ZvLO91rxfref3n+VZR+CYIgyQAhMapkNZ5efiF2GL9gGRYthzfE82uv8fbBR7HDMUiP2e3RfHADna5NTlDxHE3aGBNPuAjClJFpAkJSwgMjJPkp/O+lJzFwWQ94+JneUHbe0rnQoFdNna+3trfiMZr0ObjYidIvQRAkGSAkxtpWnDciXSwftBEbny+GjYNN5hdLROl6xbHo4jRY2ej+fa3btRqPEaWvmY4jFwRBcI9sLSQkp61PH0SGRIkdRpqcPR3hkdUd7x5Kf7pg8MpeaDawvt73BX8ORbe8Q1KLQwmBoimcSt4NmUwmWJ8EQfwfGRkgJKf7rPa8tl+jfSWc0eyFWxZXve+NCo2RfCIgV8rRZXobNOpT26D7vbN7okLjUhxHlbH+i7uSRIAgRESSAUJyGvauhSYD63HerqObA7rPbI9xO4YgMiQa4QERnPdhDK+c3BzPrFFp8NeM/ZjcZB40asMqCZaqk/lRy7SMRp8FndF3QReUqFUYBSvkRdXW5aG0UYKidD+1snT94mg2iEwREISYyDQBIVnvH3/C+jE78OnZF9A0jVL1iqHLtDb48OgzlvRZi9iIuAzvp2gKvrm9MH77UMjkMmQv7AeFUgEA+PY2ED3zDxPiy9CZs6cjVtz+A93zDYVWrTW6PYqiMODP7mgxtKHe98bHJKB9lr5ITkhGmq8QFNB9Znt0mtTqt6e+vw3EwT+P48zWy1Anq9Ptw9XHBa1HNkGLoQ0gV5CNTQQhJpIMECbr0/OvWNp/Hd7c/QCt5vc3T1cfFyy+PANZ8/j89tyxNWewfNBGIcLUS/7yefD69jvO2vPN7Y1tb1cYdO+tY/cxo9VCAID23/UDFECBQrHqhTDn5EQorRQZthEVFo0T68/hwdkniAmPg5uvC8o2KIFKLcvBK5uHXiMIBEHwhyQDhMljWRY3j9zDyQ3n8f1tEOycbVGzQ2XU61EDDi72ad4zsvpUPLv6SuBIxXEyaVfqiIi+3j74gP2Lj+L6wTvQqLXwyeWF5oMaoPGAupkmAgRBmA6SDBAWqU/Rkfj8/JvYYfCOltE4lbwbNG3c8iCWZcFoGcjkZJEfQZgjsoCQsEhZ8/qKHYIgyjcuZXQiAKSsPyCJAEGYL5IMEBapUd86YocgiLZjmokdAkEQJoAkA4RFKlWnKGp0qCx2GLwauaE/ClXMJ3YYBEGYALJmgLBYWq0WfYuOxtdX38UOhVMyOY3ZxyegdN3iYodCEISJIJt7CYslk8kgV5rXPHjhygUw88jYdHdREARBpIVMExAWzdnDSewQOGXnZEMSAYIg9EaSAcKilWtcUuwQOMVI8PhngiCkjyQDhEWr27U6YCZF8CiaQuFK+cUOgyAIE0SSAcKi2TvboX6Pmgbdm6tYdljbW3MckeEUSjka9K4ldhgEQZggkgwQFm/gsh7wy+er1wiBq48z5p2djGMxO7D17XI0H9wAMrm4f05T9o2Ci6d5rYEgCEIYZGshQQCIj47Hpgm7cGbrJaiS0j9pDwCKVi+ASbtGwNXb5ZfHQ7+FoUf+YVAlZnw/H2wcrHE0eofg/RIEYR5IMkAQP0mITcS3N4F4//AjTm++iNd33wNIqfFfpGoBDFrWAzkK+6d7f3hgBPoWG42Y8FihQgYA0DSFk8m7IZOZ11ZJgiCEQZIBgshAeFAkEmIS4J7FFTb2Njrd8+nZFwwqOx7qZA3P0f3qSPR22DroFiNBEMTPyJoBgsiAm48L/PJl0TkRAIAcRfyx6NIM+Ob2/u05ikoZZQCAis3LIG/pXJzEaetoA2s7K07aIgjC8pCRAYLgCcuyeHr1Jb6+/I6YiDiEB0YgMS4J7lncUK979dSTE4+tPYPlAzca3A8to9FyWCP0W9SVq9AJgrAwJBkgCAk4suo0Vg7ZpPd9tJyGm48LVt2bT3YSEARhMDJNQBAS0GxQfdTqWEWveyiaQoXGpbH85hySCBAEYRQyMkAQEnLt7zvYNm0PvrxIOUkxS14fZMnljQ9PPiMuKh5e/h4oXa84StYuilzFs8Pd11XkiAmCMAckGSAICVIlq8EyDKxsyKJAgiD4R5IBgiAIgrBwZM0AQRAEQVg4kgwQBEEQhIUjyQBBEARBWDiSDBAEQRCEhSPJAEEQBEFYOJIMEARBEISFI8kAQRAEQVg4kgwQBEEQhIUjyQBBEARBWDiSDBAEQRCEhSPJAEEQBEFYOJIMEARBEISFI8kAQRAEQVg4kgwQBEEQhIUjyQBBEARBWDiSDBAEQRCEhSPJAEEQBEFYOJIMEARBEISFI8kAQRAEQVg4kgwQBEEQhIX7H5SfPQnP1u72AAAAAElFTkSuQmCC"
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"tsne = TSNE(n_components=2)\n",
"X_tsne = tsne.fit_transform(hs)\n",
"\n",
"plt.scatter(X_tsne[:,0], X_tsne[:,1], c=ys)\n",
"plt.axis('off')"
]
},
{
"cell_type": "code",
"execution_count": 170,
"id": "041b1eba-cc0d-401e-94a3-cdcf73c1a952",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:22.270261Z",
"start_time": "2024-06-14T15:43:22.259355Z"
}
},
"outputs": [],
"source": [
"import pandas as pd\n",
"\n",
"df = pd.DataFrame()\n",
"df['x1'] = X_tsne[:, 0]\n",
"df['x2'] = X_tsne[:, 1]\n",
"df['label'] = ys\n",
"df.to_csv('./line_10_only_classification.txt', index=None)"
]
},
{
"cell_type": "code",
"execution_count": 170,
"id": "67cb3128-4fcf-4332-b6d7-87642e818b10",
"metadata": {
"ExecuteTime": {
"end_time": "2024-06-14T15:43:22.272408Z",
"start_time": "2024-06-14T15:43:22.271001Z"
}
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"name": "python3",
"language": "python",
"display_name": "Python 3 (ipykernel)"
},
"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.16"
}
},
"nbformat": 4,
"nbformat_minor": 5
}