Skip to content

Instantly share code, notes, and snippets.

@lan496
Last active November 26, 2021 02:29
Show Gist options
  • Save lan496/624d5068aba04f2e34b20288003f6759 to your computer and use it in GitHub Desktop.
Save lan496/624d5068aba04f2e34b20288003f6759 to your computer and use it in GitHub Desktop.
Benchmark for creating structure objects
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "affecting-efficiency",
"metadata": {},
"outputs": [],
"source": [
"from itertools import product\n",
"from time import time\n",
"\n",
"from pymatgen.core import Structure, Lattice\n",
"from ase import Atoms\n",
"import numpy as np\n",
"from tqdm import tqdm\n",
"import seaborn\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "1db75a92",
"metadata": {},
"outputs": [],
"source": [
"seaborn.set_context('poster')\n",
"plt.rcParams['grid.linewidth'] = 0.2"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "bound-alias",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"pymatgen 2022.0.16\r\n"
]
}
],
"source": [
"!pip list | grep pymatgen"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "0c127776",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"ase 3.22.0\r\n"
]
}
],
"source": [
"!pip list | grep ase"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "a92dcefd",
"metadata": {},
"outputs": [],
"source": [
"def get_coords_and_matrix(n):\n",
" a = 2.856\n",
" matrix = a * n * np.array([\n",
" [0, 0.5, 0.5],\n",
" [0.5, 0., 0.5],\n",
" [0.5, 0.5, 0],\n",
" ])\n",
" \n",
" coords = a * np.array(list(product(range(n), repeat=3)))\n",
" return coords, matrix"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "dc4a6287",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"100%|███████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 8/8 [01:00<00:00, 7.55s/it]\n"
]
}
],
"source": [
"list_num_sites = []\n",
"times_mg = []\n",
"times_ase = []\n",
"times_ndarray = []\n",
"\n",
"for n in tqdm([2 ** i for i in range(8)]):\n",
" start = time()\n",
" coords, matrix = get_coords_and_matrix(n)\n",
" times_ndarray.append(time() - start)\n",
" \n",
" num_sites = coords.shape[0]\n",
" species = ['Al'] * num_sites\n",
" list_num_sites.append(num_sites)\n",
" \n",
" # pymatgen\n",
" start = time()\n",
" Structure(lattice=matrix, species=species, coords=coords, coords_are_cartesian=True)\n",
" times_mg.append(time() - start)\n",
" \n",
" # ASE\n",
" start = time()\n",
" Atoms(f'Al{num_sites}', positions=coords, cell=matrix)\n",
" times_ase.append(time() - start)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "d00d5d16",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x7f7f0d1badf0>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAr0AAAJ7CAYAAADulSlxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAA9hAAAPYQGoP6dpAADEPklEQVR4nOzdd5hU5fXA8e/ZxrKFDouANGkCIoooKF0RO9gQCwpiYo1BTDCGJBpLoonyM0aNWFEsCIhiFEWUqmAoCkpRQJq0pcMuy7Lt/P64M8Ps7tTdu5XzeZ77zJ25733v2ZnLcvadt4iqYowxxhhjTHUWU9EBGGOMMcYYU9Ys6TXGGGOMMdWeJb3GGGOMMabas6TXGGOMMcZUe5b0GmOMMcaYas+SXmOMMcYYU+1Z0muMMcYYY6o9S3qrGBGJF5ELReRZEfleRI6ISLaIbBCR50WkRUXHaIwxxhhT2YgtTlG1iMgFwGzP01+A5Z79s4EmwGHgIlVdXAHhGWOMMcZUStbSW/UUANOAc1W1uapeqapXAqcAE4FawLsiEl+BMRpjjDHGVCrW0luNiEhNYCdQG+inqvMrOCRjjDHGmErBWnqrEVU9CqzzPG1SkbEYY4wxxlQmlvSGICKxItJZREaIyL9FZLGIZImIeraJJaz3ChGZKiKbPYPQdovIIhH5vYjUKk28QEvP010lrccYY4wxprqJq+gAKrkpwFVuVSYiKcDbwBVFDjX0bD2B34jIUFX9pgSXGO6pZw+wqDSxGmOMMcZUJ9bSG1pskef7gfUlqcjTCjuV4wlvOvAYcANwD/C15/WTgZkicmqU9bcEnvY8Haeqx0oSpzHGGGNMdWQtvaEtAdbiTAu2XFU3icgI4PUS1HUbcJFnfw0wQFXT/Y4/LyJPAfcDdYEJQJ9IKvZ0iZgB1AOmqurLJYjPGGOMMabastkbolQk6X1DVUdEcE4szpy6J3le6qaq3wYptwzo6nlpkKp+HqbuROAzoC/wJXCptfIaY4wxxhRm3RvKRx+OJ7zzAyW8AKqaDzzr99L1oSr1zMX7Pk7C+w0w2BJeY4wxxpjiLOktHxf77c8MU/bTIOcVIiIxwFvAJcBK4BJVPVLiCI0xxhhjqjFLesvHaX77S0MVVNVdOF0hANJEpGHRMiIiwMvAUOAn4EJVPeBSrMYYY4wx1Y4NZCsf7f32N0VQfhPOLA7ec/cUOf40cKun3PmqujuaYESkOdA8mnP8NATOw5nJ4gCwEbAWZmOMMcZ4JQOt/Z5/rKq/BCtcXizpLR91/Pb3RlB+X5BzEZHBwH2ep5uBx52G32I+VNUPg9R/K/BQBHEYY4wxxrjhPxUdgCW95SPFbz87gvJH/fZTixyr67ffP0Qdm4EPI7iWMcYYY0y1Z0lvFaOqE4GJFRyGMcYYY0yVYklv+cjkeAttoud5KDX99jPKIJ7XgC9KeG53YLz3yfPPP0+XLl1cCcqYspSVlQVAUlJSBUdiTOTsvjVV0ffff8/dd9/t/9LGiorFnyW95eMgx5PeBoRPeusXOddVqroV2Fr0dREZA4wJc3qC/5MuXbrQq1cvF6MzpmxkZjr/7FJSUsKUNKbysPvWVBOVYsC7Jb3l4yeglWe/FU5/21Ba+e3/VBYBBVELaFqO1zPGGGOMKReW9JaPH4CLPPvdgbnBCopIGsenK9utqkWnKytLh4HtYcok4ExbZowxxhhTZdjiFOXjM7/9oKuseVzitx9u9TZXqep4VW0WagOuKs+YjDHGGGPcYC295WM+sAtoDPQTkTNV9duihUQkFrjX76XJ5RSf9/pR9+k1xhhjjKkKLOktB6qaLyKPAC94XnpTRAYEWEntCaCrZ/9rVZ1VXjF6WJ9eY4wxxlRLlvSGICKtgFFFXvafn+sMEXmsyPE5qjonQHUvA1cCA4FOwEoReRlYA9QDrge80yAcBG4vXfQlYn16jTHGGFMtWdIbWgtgXIjjXSicBAPkAcWSXlXNE5GrgXeAy3C6Ovw5QJ3bgOtUdXWJIi4FVR2P3xy8gYhIL2Bh+URkjDHGGOMOG8hWjlQ1Q1UvB4YA04FfgGPAXuB/wANAZ1VdVGFBGmOMMcZUQ9bSG4KqzgOkDOqdAcxwu97SsoFsxhhjjKmuLOk1/mwgmzHGGGOqJUt6jb9yHchWUFBAZmYmhw8fJicnh/z8fDeqNSYg7/0VGxtbwZEYEzm7b015iY2NJSEhgVq1apGSkkJMTPXrAWtJr/Epz4FsGRkZbN++HVUtbVXGRKSgoADA7jlTpdh9a8pLXl4ex44dIyMjAxGhadOmpKamVnRYrrKk15S7QAmviFhLhilT3vtNxPVu+saUGbtvTXnJz8/33W+qyvbt26td4mtJrylXBQUFhRLelJQU6tWrR1JSkv1SN2XKviY2VZHdt6a8qCpZWVns37+fzMxMX+Lbrl27atPVwZJe41Meszd4/yGBk/A2a9bMkl1jjDGmgokIycnJJCUlsW3bNt//15mZmdSqVauiw3OFJb3GX5nP3nD48GHffr169SzhNcYYYyoREaFevXpkZmYCzv/blvSa6qjMZ2/IyckBnH9USUlJJa3GGGOMMWXE2+VQVX3/b1cHriW9IhIDdAZOB1riLLOb7Dl8BNgFbAZWAqtUtcCtaxt3lMfsDf7906yV1xhjjKl8vIPL8/LyqtV0oqVKekWkLjAUuAzoB0TadJclIvOAj4EpqnqgNHEYY4wxxhgTSomG44lITxF5D9gJvABcgtOqKxFuyZ5zXgB2ish7InJu6X4UY4wxxhhjAouqpVdEegOPA+d5XypSZCewFqdf6F4gy1OmJk4/0KbAqThdH7znJgDXANeIyNfAOFUt9eIHJnrlMXuDMcYYY0xFiCjpFZFWwP8Bl3tf8jzuAD4E5gALVXVPhPU1BHoDA4DBHJ8x4Dxgnoj8F7hPVTdFUp9xTZnP3mCMMcYYUxEibeldg9PCJzitt1OBV1T165Jc1JMcT/ds93i6NvwKuBanX/AVwECOD4Qz5aPMZ28wxhhjjKkIkSa9NYADwL+BZ1V1v5tBqOoiYJGI3A/cC/wGqOPmNUx45TF7gzHGGGNMRYh0INuTQGtVfdjthNefqu5X1YeB1p5rGmOMMcYYU2oRtfSq6oNlHUiR6x0C/lie1zTGGLesWLGCDz/8EIAhQ4bQtWvXCo3HVH/79+/nrbfeYtasWXz//ffs37+fnJwcUlJSaNq0Ke3atePMM8+kT58+9OzZk/j4+GJ1nIj37YcffsiKFSsAGD16NHXq1KnQeEzZshXZjDHGZStWrOCvf/0rAC1btjwhkgdTcd58801++9vfcvDgwWLHDh48yMGDB1m9ejUffPABAEOHDuW9994rVvZEvG8//PBD3njjDQBGjBhhSW81Z0mvMcYYU0VNmDCBO+64w/f8tNNO44orrqBdu3YkJSWRkZHBhg0bWLJkCQsWLCAnJ6darbBlTDQs6TUmgIzsXHYdyuZITj7JCbE0rp1IamLxrwONMaai7Nq1izFjjk+t/sILL3DnnXcGLZ+RkcG0adPYvj3cJD3GVE/llvSKSA3gbpzpyNKA3cBMnNkgjpZXHCa4E31xClVl8cZ9TFq8hc/XpJNfoL5jsTHCoE5p3NSjBT1b10ek6LosxhhTvqZPn05WVhYA1113XciEFyA1NZWRI0eWR2jGVEolWobYn4i0FJGtnu2FIGUSgbnAP3EWpWgH9AL+BiwXkUaljcO4wrs4RaitWs7Ru2r7IQY9s4AbXv4fn67aVSjhBcgvUGb+sIsbXv4fg55ZwKrthyooUvfNmzcPEUFEePjhhwH46aef+M1vfkP79u1JTk6mbt269OjRg//7v//j2LFjxer4/e9/76vj7bffjui6I0eO9J3z2WefFTrWsmVLRISWLVsCkJOTw3PPPce5555Lw4YNSUlJ4fTTT+cf//gHR44cKXRueno6Dz/8MKeffjq1a9cmNTWVHj168Oqrr6Ja+HMtKjs7mxkzZnDvvff6rhUfH09qaipt27Zl+PDhzJ49O+j5EydOREQKJRb+P6f/Fkhubi7PPfcc5513HvXq1SMpKYl27dpxzz338NNPPwHw8MMP++qYN29e2J9nwoQJXHbZZZx88skkJiZSu3ZtOnfuzL333su6detCnh/oWsuXL2fkyJG0bt2axMRE6tevT//+/Zk4cSIFBQUh64vW7t27+dvf/kb//v1p0qQJNWrUIDk5mXbt2jFs2DDeeOONYp+/v/z8fN544w2uuOIKmjVrRmJiInXr1qVLly6MGTMm7M8f6N/G+vXruf/+++nUqRN16tQpdMxfad/7SPz444++/b59+5a4npLet4Huj3nz5jF8+HBat25NUlJSsfvUW75fv35h4+rXr1/Ify/+or1XRowYgYj4+vMCtGrVqtjPWzTOor+bQvFeQ0TYvHlzsePe911EmDhxIgDffvstd9xxB+3atSM1NbXQMX+HDh3i6aef5oILLvD9vPXq1aNbt248+OCD1pofjKqWagPuAwqAfOCiIGWe9CtTEGD7orRx2Fb6DaeVd1uYbTeg3m3hwoUajXXr1umaNWt03bp1UZ1Xlhas262n/vlTbfHAxxFvp/75U12wbndFh+6KuXPn+j7Phx56SCdPnqxJSUnq/zn7bx06dNDNmzcXqmP9+vUqIgponz59wl7zwIEDWrNmTQW0ZcuWmp+fX+h4ixYtFNAWLVrozp07tVu3bkHj6d69ux44cEBVVRcvXqxpaWlByw4bNkwLCgqCxtWqVaug5/pvgwcP1oyMjGLnv/766xGd7/zqLWznzp16+umnBy2flJSk7733nj700EO+1+bOnRv0Z5k3b542bdo0ZAyxsbH6t7/9LWgdRa/1xBNPaGxsbND6hgwZorm5uUHri8YzzzwT8j70bqNHjw54/oYNG7RTp04hz42Li9NHH300aAxF/21MmjTJd9/6bw899FCh89x4773y8vI0Ly8v4LE777wzaAzRKOl9W/T++M1vfhPwHP/71Pta3759w8bVt2/foP9e/JXkXrnlllsi+nmLxun/uykc/2ts2rSp2HH/9/3111/XJ598MuC/r9dff73QeVOmTNF69eqFjDsxMVEnTpwYNsZQSvP/9cKFC4vG1EsrQZ7jRveGcz2PWTjLERciIik43RrU89IrwEdAc+BhoAHQX0QGqGqx80350RNwcYpV2w9x+6TlZOVEN7AjKyef2yctZ8rtPenctHYZRVf+li9fzhNPPEFubi433HAD559/PjVr1mT16tW89tpr7Ny5kx9//JH+/fvz3XffUbu287O3adOGgQMH8vnnn7NgwQJ++ukn2rdvH/Q6kyZN4uhRp1fTr371K2JiAn/plJuby9VXX83y5csZOHAgQ4YMoUGDBmzcuJHnn3+ebdu2sXTpUkaPHs3DDz/MoEGDOHr0KCNGjKBPnz7UrFmTpUuX8p///IejR48yefJkBg4cyK233hrwellZWdSpU4cBAwZwxhln0KJFC5KSkjh8+DDff/897733Hjt37mTGjBnceuutTJkypdD5AwYM4IMPPmDOnDn8+9//BuA3v/kNAwYMCPm+Z2dnM3DgQFatWgVAgwYNGDVqFF26dCEnJ4eFCxcyadIkbr75Zi666KKQdQF8+umnDB48mNzcXGJiYrjooou44IILaNq0KdnZ2Sxbtow333yTQ4cO8cc/OrNDPvhg6JkpX375Zd555x0aNmzIiBEj6NKlCzExMSxatIhXXnmFY8eO8eGHH/KPf/zDV2dJ/e53v+Ppp5/2Pe/Tpw+XXnopzZs3Jz8/n82bN/PVV1/x5Zdfev9gL2THjh2cd955pKenA9CiRQtGjBhBhw4dyMzMZNasWbz//vvk5eXx5z//mWPHjvHoo4+GjGnRokU8/vjjiAi33HILvXv3Jjk5mQ0bNtC8eXNfubJ474Np06aNb/+NN97gt7/9LXXr1o26npLet/7+8Y9/8Omnn9KwYUNuueUWunTpAsB3331HrVq1oo4pUiW9V+69916GDBnCs88+y9y5cwFnUGCjRoW/eG7QoEGZxe5vypQpfPrpp6SkpHDzzTdz9tlnEx8fz5o1a2jcuLGv3Msvv8ztt9+OqpKQkMDgwYPp06cPaWlpZGZm8tVXX/HOO++QnZ3NiBEjSEhI4Prrry+Xn6FKKG3WDKzCacFdFOT49Rxv5X2lyLGBfscmVPRfALZF9Hn3wu+vt6rc0ltQUKADx8+LqoW36Hbh+PkhWw6rAv/WLDwtioFaEA8cOKA9e/b0lbv99tsLHf/ggw98x8aMGRPymp07d/a1tO3cubPYcW9rCqAioq+++mqxMrt27dLGjRv7Ws26du2q9evX12+//bZY2S+//NJXX6dOnYLGNXPmTM3JyQl6/MiRIzpkyJCw93/RFpxw/vKXv/jKd+7cWdPT04uVWbp0qdauXTtoC5rXjh07fK1AjRo10sWLFwe85rZt23yfQ2xsrK5du7ZYGf+WPDytXgcPHixWbt68eb4WqgYNGuixY8fC/szBvP/++77rpaSk6EcffRS07K5duwJ+BpdccomvjksuuUSPHDlSrMzMmTO1Ro0aCmhMTEzA96nov41GjRrpypUrg8bj5nvvFaqld/PmzZqQkOCLr0WLFvrUU0/pjz/+WKLfS9Het0Xvj549e+r+/ftDnuN/L4UTrqXXjXslXGtsUWXV0gtou3btdMuWLUHrW7lype/zbtu2bdD7Zs2aNdqkSRMFNDU1Vfft2xc21kCqY0tvqfv04vTzBPg5yPHz/fb/7X9AVWcDGz1Pu7sQizERW7xxH+vSM0tVx0/pGXyzscwWKawQTzzxRMD+dnXq1GHq1KmkpKQATn+0PXv2+I5ffvnlNGvWDHBanQL1/QWnxczbojl48OBCrRiB3HbbbQFbZtPS0rjnnnsAp+/mihUreP755znjjDOKlR0wYICv1Wr16tX88ssvAa918cUXB5y03yspKYk33niD5ORkgEL9AUsqJyeHF15whkPExcXx3nvvFWttAjjrrLMKtWgF889//pP9+517ctq0afTo0SNguaZNmzJ16lRiY2PJz8/nX//6V8h669Wrx/vvv+9r3ffXt29frrnmGgD27t3L0qVLw8YZiKry5z//2ff8tdde4/LLLw9aPi0tjV69ehV67YcffmDmzJkAnHTSSbz77rskJSUVO/fiiy/2zUlbUFDAk0+GXwR0woQJvhbMQMrqvQ+mRYsW/POf//Q937JlC7/73e/o0KEDdevWpX///owdO5YZM2aQkZFRomtEKjk5malTp5aopbkk3LhXKhMRYfLkyYW+NSjq4YcfJicnh8TERGbOnEmHDh0Cljv11FN9/YAzMjJ4+eWXyyLkKsmNpDfZ83g4yPHensftqroywPHvAQFauBCLMRF765stlaqeyqBOnTr86le/Cnq8adOm3HjjjQAcO3aM//73v75jsbGx/PrXvwZg3759vP/++wHreOmll3z7t99+e9iYfvOb3wQ95v+fWFpaGtdee21EZdesWRP2usHUqlWL0047DYBvvvmmxPV4ffXVV+zduxeACy+8kI4dOwYte/PNN1O/fv2gx1WVN998E4CePXvSu3fvoGUBOnTowNlnnw3ArFmzQpYNd+2BAwf69r1/1ERr+fLlvs/mzDPPDPl5BjN9+nTf/h133BHyq/V77rmH1NRUAGbOnEl2dnbQsi1atGDw4MFBj5flex/Kvffey4cffkjbtm0LvX7o0CHmzZvHP//5T4YMGUJaWhojR45ky5ay+X111VVX0bRp0/AFXeLGvVKZ9OrVK+Af7F4HDx5kxowZAFx55ZWFurYEMnDgQE466SSgdPdXdeNGn958IBaoUfSAiDQE2uI0bQfrB7rb85jqQizGRCQjO5dZq9Ndqeuz1bvIyM6tFvP49urVi8TExJBlLrjgAiZMmADAkiVLCrXC3nbbbTz66KPk5uby0ksvccMNNxQ69+DBg75+sKeccgoXXHBByGslJyfTuXPnoMf9W4m7desWtG9w0bIHDhwIWu7AgQO8/fbbfPbZZ6xatYp9+/Zx5MgRb/eeQrZt2xYy/kj4t4r2798/ZNn4+HjOO+88Pvroo4DH16xZw759+wCoW7eub0nZUGJjYwHYtGkT2dnZQT//nj17hqzH28oPod/fUBYuPP7fRKgEM5T//e9/vv0LL7wwZNnk5GR69erFp59+Sk5ODt99913Qn/O8884LOYtAWb734QwePJjLL7+cefPm8d///pdFixaxcuXKQt+2HD16lIkTJzJ9+nQmT57MxRdfXKJrBRMuyXebG/dKZRLu/fv66699s6PUqFEjovsrNTWVnTt3luqP/OrGjaR3P9AYaBngmP9v8EVBzvf+K89xIRZjIrLrUHaxaclKKr9AST+cXS2S3qKtReHK7Nixo9Cxk046iSFDhjB16lTmz5/PunXraNeune940QFs4aYiqlevXsgyNWoc/1s7VCtk0bLBWvRmzJjBqFGjfMlLOIcPB/uCK3L+7+Epp5wStnzr1q2DHvOfFmnmzJm+r/kjtX//fpo0aRLwWLgBPZG8v+H4/xERqsU7lJ07d/r2/e+9YNq1a8enn35a7Nyi/JP6QMryvY9ETExMoW48ubm5rF69mq+//prp06czZ44zTvzw4cNcc801rFy5MmxrYTTCvT9uc+NeqUyiub8mTpwYcBqzYLxdbow73RvW4HRP6CEiRTt7Xee3vyDI+d7vQ/YEOW6M645EOVtDOJnHqseynt6+qpGWCdRP0H+CfP+uDP7P4+PjI5okP1TLbWnKBrJ48WKuueYaX8Lrncv1hRde4N1332X69Ol88MEHfPDBB3Tq1AnAlXlp/ecODdT3tKhQn9HBgwdLFUtOTvC2h9K+v5Hw/yPC23c8Wv73ZCT3s/91QvV7rVmzZsh6yvK9L4n4+Hi6du3K3XffzZdffsns2bN9LclZWVkR9WGORrj3x21u3CuVSVneX7m5uSU+t7px47eYd1b5msBLIpIEICJDgSE4XRs2q+oPQc4/w1NmgwuxmFIQkTEisi3UBkwPW1EVkJwQ62p9KTXcra+ihJroP1AZb39If/379/cNsHjjjTd8/5n7D2C78sorAw7Wqkh/+ctfyMvLA+D5559n5cqVPP3009x5550MGzaMK6+8kiFDhjBkyBBX/4P3T8y8q2uFEuoz8v/Pf8yYMVGPbI5kwv2y5N//NjOzZINM/e/JSO5n/+sEup8jVdnf+wsuuIC//OUvvuehFlmpLPLzgzcmuHGvlLVQ8UfL//569tlnSzLzksGdpHci4O3AdQ2wV0T2Au/itAADPB/oRBHpBtTzPF3uQiymdE6YFdka104kNsadpYTjYoS0WiXri1fZbNgQ/m9P/zLBvo71tvbu3bvXN6DN2w8YIhvAVp5yc3N9q0Z169aNu+66K2T5QKsrlZT/e/jzz8EmwTlu48aNQY/5f0UabIaKysw//pL2Q/QO3gFn9bRw/FdGK033gqrw3vv3oQ/VlaMsJSQ4K9lH0rLtHeAZiBv3Skl4u/GUNv5oVYX7qyooddKrqvuAW3EGtAlOH916HE94vwGeDXL6ML/9arPgQRV2GNgeZqsW3VBSE+MZ1CnNlboGdWpcLfrzgjM4JNhUY15ffPGFb/+cc84JWObmm2/2fVX/0ksvcfDgQaZOnQo4fYLDDdgqb3v37vW18obr57h06dKw/5n5dwUI18rSvfvx2Rq9k+QHk5uby9dffx30eNeuXX1Tis2dOzfsZ1nZ9OnTx7fvHakeLf978vPPPw9ZNisri6+++gpwkrFQo+fDqQrvvf90fIG6BERz35aUd0qzcMvk7tu3L+RSzW7cKxD9z+yNf8+ePSET37y8vBJP3RdI7969feMbii7bbiLnSictVZ2BszLbDGAfkIsz/+4TwEBVzSt6jogkA6M8T48RYDU3U75UdbyqNgu1AVdVdJxuuamHO7PkuVVPZXDw4EFeeeWVoMd37tzJ22+/DTgtHpdddlnAcnXq1PGtAjRv3jz+/Oc/+waw/frXvw47gK28+XcxCNfa/dBDD4Wtzz+hCPcV+3nnnecbhPf555+HbLV68803Qw6yi42N9U0pt3fvXsaPD7nAYqVz5pln+vpLf/vtt74/lKJx9dVX+/b/85//hBxs+Pzzz/v68V566aWFBuNFqyLe+927d4cv5Mc/OfROu+cvmvu2pLyf79atW0O2xD/zzDMh+8y7ca9A9D+z95p5eXnMnz8/aLk33ngj4gGxkWjUqJFvxo0ffviBd99917W6TySujUxQ1WWqeqWqNlTVGqraRlX/qKrBOqnl4vTnbQW0UdWjbsViTCR6tq5Pu7TSDYBon5ZKj9b1whesQh544AEWLCg+7vTw4cMMHTrUlySMHDmShg2D93bx7yLw3HPPAU5r2ogRI9wN2AW1atXyjfRfvnw506ZNK1YmPz+f++67zzfSP5RWrVr59r/99tuQZWvUqMHdd98NOP+RXnfddQGTmWXLlnH//feHvfYf//hH6tSpA8Cf/vSnsMnDkSNHeOWVV8rlP9ERI0YgIogIDz/8cLHjIsJjjz3me37rrbcWmgu6qD179hRr+e7cuTOXXnop4PyRdsMNNwTsKz1r1ixfH9eYmBgeeOCBkvxIhbj93osIcXFxxMXFBexSM378eLp168Zbb70VNmF77733+Pvf/+57Pnz48GJlorlvS8p/qrT7778/YL/XadOm8cQTT4Ssx417BaL/mf3j/+Mf/xjw3lqwYAH33Xdf2Lqi9fjjj/u6h9x2221h/83u37+f8ePHF/p27kTnxpRlJaKqOUD1mdXfVDkiwvihXRk6YTFZJZjNISkhlqeHnl7pWi1L47LLLmP27NkMGDCAYcOGcf7551OzZk3WrFnDq6++6pteq1WrVmFHf5955pmcffbZLFmyxPfa1VdfXW5r2Udr9OjRvkR96NChXHfddfTt25e6deuyYcMG3n77bdauXUvnzp2pUaMGy5cHH4Zw2mmnkZaWRnp6Om+99RYNGjSgR48ehWZnuOiii3z7Dz74INOnT2fVqlWsWrWKTp06MWrUKE4//XRycnJYsGABkyZNIiYmhiuuuMI3T2+gGRWaNm3KlClTuPzyyzl27Bj33XcfL7zwAldeeSUdO3YkJSWFjIwMNm3axLJly5gzZw7Z2dk8+uijbr2VpTJkyBDuv/9+nn76aTIzM7niiivo27cvl156KSeffDIFBQVs3bqVr7/+mtmzZ3PHHXdw3nnnFarjpZde4swzzyQ9PZ1PPvmETp06MXLkSNq3b09mZiaff/45U6dO9X2dPW7cuKBddaJREe/9t99+y/Dhw0lKSqJv376cffbZNG/enNq1a5OVlcX69euZOXNmoft14MCBAWdPifa+LYlbb72VJ598kr179/Lf//6Xnj17cvPNN/uu+/HHHzNr1ixOPfVUEhMT+e6774LW5ca94t/PeezYsezevZv27dv7uoLUq1fPt4gIOHMCd+jQgR9//JFly5Zx5pln8qtf/YoWLVqwf/9+Zs+ezfTp02nUqBHdu3f3TRXnhq5duzJhwgRGjRpFVlYWN9xwA//4xz+4/PLLadu2LTVr1uTQoUNs2LCBJUuWsGDBAvLy8pg0aZJrMVR5Fb0Osm1VawN64beedqC1zEMpzVreZWXBut166p8/1RYPfBzxduqfP9UF63ZXdOiumDt3ru/zfOihh/S9997TpKSkouum+7b27dtHtEa9avG15efNmxfReZGub79p0yZf3bfcckvIsq+++qqv7Ouvv17seEFBgd56661Bf25ATzvtNN24caP27dvX91owr7zySsi6itqxY4d26dIlaPmkpCR97733dNy4cb7Xvv3226DXX7p0qbZv3z5kDN4tNjZWX3755WJ1PPTQQ74yc+fODfn+Fr2PArnlllvClvH6xz/+oYmJiWFjv++++wKev379eu3YsWPIc+Pi4vSRRx4p1c8UiBvvvaoWKhfo39zLL7+sKSkpEV0H0JiYGL3jjjv06NGjQWOP5r6N5v7w98UXX2hycnLQa3Tu3Dnif2eqpb9XbrrppqDn9O3bt1j5lStXaqNGjYKe07x5c12+fHmh+z3Q5+f/+zHQ76RgPvvsM23SpElEn3mNGjX0008/jbhuf6X5/3rhwoVFY+mllSCHqbCWXmMqi95tGzLl9p6MmbKCdenhp75pn5bK00NPp3PT2uUQXfkbOnQop59+Ov/+97/5/PPP2b59O/Hx8bRv357rrruOu+++O+K+j/4rYnXo0IG+ffuWVdilJiK8+uqrXHrppbz00kssW7aMw4cPU79+fdq3b8+1117LqFGjIl41a9SoUbRo0YIXX3yRpUuXsnv37pCLNpx00kksW7aMCRMm8M4777B27VqOHTtG06ZNufDCC7n33ntp3759ocFu9eoF71pz1llnsWbNGqZPn86MGTP43//+R3p6OkeOHCElJYWTTz6Z0047jX79+nHFFVcUWrGuMvj973/PjTfeyEsvvcTnn3/O+vXrOXjwIDVq1KBZs2Z069aNSy+9lKuuCjzMoE2bNqxcuZK33nqL999/n2+//Za9e/dSs2ZNTj75ZC644ALuvPPOiBawiFZ5vfe33XYbw4cPZ968eSxYsIBly5axYcMGdu/eTVZWFklJSdSrV49TTz2VXr16MWzYsLADNaO9b0vi/PPP54cffuDJJ59k9uzZbN++nZo1a9KuXTuuv/567rjjjqhWpyvtvfLGG2/Qt29f3n33XVatWsWBAwdCzm3bpUsXfvjhB/75z3/y8ccfs2XLFuLi4mjdujVXX30199xzj2/AW1kYNGgQGzdu5J133vG14u/Zs4fs7GxSU1Np2bIlp59+OgMGDODyyy/3dbkxIOq03oUuJHKmqpZNB5/g1zxDVb8rz2ua8ESkF34zbSxcuJBevXpFfP769evJy8sjLi4uotW/ypOq8s3G/Uz6ZjOzVqcXWrEtLkYY1KkxN/VoQY/WoVcJq2rmzZvnm03hoYceCtjXsqSeeeYZX9+28ePHl0k/t0h5+w56l36tqrp168a3335LnTp12L9/f7W6F01x1eW+NVVPaf6//uqrr4ourdxbVb9yNcASiLSld4mITAYeVdWfyjIgEWkPPARcC1SPeaBMlSAi9DylPj1PqU9Gdi7ph7PJPJZPSo1Y0molVptpycqLqvLiiy8CzmpDt9xySwVHVPUtXrzYN9imX79+lvAaY0wUIk16Y4DrgWEi8iHwgqq6OsWYiJwP3A1cwfE5fo2pEKmJ8ZbkltJbb73FTz85fyPffPPNIb+KN86ApLZt2wZdFWzNmjW+aeCg8HLPxhhjwos06R0DPIyzYteVwJUisgmYCnwILFXVqBahF5EYoDvOUsXX4kxdBk7CexintdcYU0Xs37+fJUuWkJOTw9KlS33zlNasWZNx48ZVcHSV32uvvcbEiRO58MILOeecc2jevDlxcXHs2rWLBQsW8OGHH/oW0Bg2bFih/tLGGGPCiyjpVdVnROQdnMR3FE63g1bAWM+WISLfAN8DPwLbcBapyMJJYmsCDYBmQHvgdOAcwL9JQ3Dm7n0F+KuqRjfrtjGmQn3//feF5rD0Gj9+PCeffHIFRFT1HDlyhA8++IAPPvggaJkbb7yRV199tRyjMsaY6iHi2Rs8SehdIvIk8DtgBOBdxqgWMNCzRcq/C8MR4HXgaVW1uXuNqeLq169Pp06deOCBB7jkkksqOpwqYezYsTRv3pz58+fz888/s2/fPg4ePEhycjJNmzblvPPOY8SIEZx77rkVHaoxxlRJUU9Z5klKfyMi44DrgBuB84Boh5bmA18DbwPvqWrwtSKNMWWmX79+RDKLS3nVc6Jq3rw5Y8eOZezYsRUdijHGVEslnqfXk6S+DLwsInWAAThdFrrgdH1I43hL8BFgF7AJpwvEEmCOqh4s6fWN+0RkDE7/7VASyiMWY4wxxhg3ubI4hSd5ne7ZTNVVC2ha0UEYY4wxxrjNVmQz/g4D28OUSQAalkMsxhhjjDGusaTX+KjqeGB8qDJFV2QzxhhjjKkKYio6AGOMMcYYY8qaJb3GGGOMMabas6TXGGOMMcZUe5b0GmOMMcaYas+SXmOMMcYYU+1Z0muMMcYYY6o9S3qNMcYYY0y1Z0mvMcYYY4yp9mxxiipIRM4ELgC6e7YWnkOtVHVzRcVljDHGGFNZWdJbNf0FGFzRQRhjjDHGVBWW9FZNi4EfgGWebTmQVqERGWOMMcZUYtantwpS1SdV9c+qOkNVt1d0PMZUJ//85z8REd82derUiM+dOHEiDz/8MA8//HDZBWiMMaZEyi3pFZEaIjJGROaJyFoRmS8iD4hIzfKKwRhjwnnttddCPg9l4sSJ/PWvf+Wvf/2r22EZY4wppVInvSLSUkS2erYXgpRJBOYC/wR6A+2AXsDfgOUi0qi0cZQFEYkVkc4iMkJE/i0ii0UkS0TUs00sYb1XiMhUEdksItkisltEFonI70Wklss/hjEmQl9//TU//vhjodc+//xztm3bVkERGWOMcYsbLb1XAs2ApsBHQcr8Fejh2ZciWwfgHRfiKAtTcPrOvg7cg/MzlLhlWkRSRGQGMAO4BmfWhRpAQ6An8A9glYj0CF6LMaasvPrqq779kSNHAlBQUMDEiRMrKCJjjDFucSPpPdfzmAXMKXpQRFKAuwH1vPQKcAVOErnX81p/ERngQixuiy3yfD+wviQViUgsMBXnZwdIBx4DbsB5L772vH4yMFNETi3JdYxLsg/D7h9h23LnMftwRUdkylhGRgZTpkwBoG3btjzzzDPUrOn8jfv666+jqqFON8YYU8m5kfSeipPQ/qCqOQGOXw4kefZfV9Vfq+rHqvoCcKNfuetciMVtS4AngGuB1qpaH6dLRkncBlzk2V8DnO4ZjPauqj6vqr2Apz3H6wITShG3KQlV2LQA3hsOT7aEF86BVwY4j0+2hCk3O8erYfKTnZ3NjBkzuPfeezn33HNp2LAh8fHxpKam0rZtW4YPH87s2bPD1nP48GGefvpp+vfvT1paGgkJCaSmptKyZUu6d+/OqFGjmDp1Kjk5gX5VHPfzzz/zhz/8ge7du9OwYUMSEhJIS0tjwIAB/Otf/yIrK8utH93nvffe48iRIwAMHz6cWrVqceWVVwKwceNG5s2bF/Tcfv36ISLMnz/f95r/YDjvNmLEiIDnHz16lOeee46BAwdy0kknkZCQQP369enevTt/+tOf2LFjR8jYJ06c6LuGt1V6xYoVjBo1itatW1OzZk2aNWvGlVdeyZIlS4qd/+mnnzJ48GBatGhBjRo1aNasGSNHjmTTpk0hrwvOZzV27Fi6d+9O3bp1iY+Pp169erRt25Y+ffowZswYFixYELYeY4wpc6paqg04AOQDk4IcfwUo8JQ5PcDxDZ5j35Y2lvLYgBE4Sb4CEyM8JxbY4XfemSHKfedX7sII69/lKd+yHH7+Xn7x6cKFCzUa69at0zVr1ui6deuiOq/Mbf9O9blzVB+qFX577hynfDXSqlUr9f9cg22DBw/WjIyMgHUsW7ZMGzduHFE9S5cuDVhHfn6+PvjggxoXFxfy/GbNmumyZcui+hnz8vI0Ly8v6PEePXoooCKimzZtUlXVWbNm+a554403Bj23b9++Ef3ct9xyS7FzlyxZoieffHLI85KSkvTVV18Nev3XX3/dV/b111/X559/XuPj4wPWFRMTo2+++aaqqubk5OioUaOCXrdWrVq6ZMmSoNd99dVXtUaNGmF/7uTk5KB1mNDC3bfGlJXS/H+9cOHCor8HemklyOHcmKc32fMY7Pvf3p7H7aq6MsDx74HWHF9VrDrqA5zk2Z+vqt8GKqSq+SLyLOAdLn498Hk5xHdi+3kOTL4Jco9EVn7PWnj9Ehj2FpxSGXvlRC8rK4s6deowYMAAzjjjDFq0aEFSUhKHDx/m+++/57333mPnzp3MmDGDW2+91dcNwP/8IUOGsGvXLgC6devGlVdeSdOmTUlOTubAgQOsXbuWuXPnsnJloF8DjltuuYW33noLgHr16nHdddfRrVs3atWqxe7du/nkk0/49NNP2bZtG/3792fZsmW0a9eu1D//mjVr+OabbwDo27cvLVu2BOCCCy6gadOmbN++nenTp3Po0CFq165d7PzHHnuMvXv38qc//YnVq1cD8MEHHxQr17x580LPv//+e/r37+9rYe7YsSPDhw+nVatW7N+/nw8//JDPP/+crKwsRo0ahaoyatSokD/LJ598wvvvv0/9+vUZNWoUXbp0ITc3l5kzZzJlyhQKCgq49dZb6dmzJ88++yyvvvoqnTp14qabbvJd94033uB///sfhw8f5oYbbmD16tUkJCQUus53333Hr3/9a/Lz84mNjWXQoEEMHDiQRo0aERMTw+7du1m5ciWzZ89m//79kX0QxhhTlkqbNQNHcVpqXwlwrCHHW3nfDnL+i54yORX9F0CEP+8Iom/p/YffOWPDlG3sV3ZXhPVbS29Jbf9O9bGTImvhLbo9dlK1afGdOXOm5uTkBD1+5MgRHTJkSNDPferUqb5jY8aMCXmt1atX6+7du4u9/uKLL/rquPzyy/XAgQMBz3///fd9LcHnnXde+B/OI1SL2ZgxYwq1lPr7wx/+4Dv2wgsvhLyGf4tvOPn5+dq5c2df+dtuu01zc3OLlXvllVdURHwtvt5WaH/+Lb2Annnmmbpnz55i5f76178WKiMi+qtf/arY+5Kbm6sDBgzwlZ0yZUqxuu6++27f8Y8++ijoz1lQUKDz588P+36YwKyl11SU6tjS60afXu+f8C0DHOvvt78oyPmJnsfQnfyqttP89peGKqiqu4BfPE/TRKRhmUV1olOFD+6IvIW3qNwj8OGd1aKP78UXX0x8fHzQ40lJSbzxxhskJztf7LzxxhuFjm/YsMG3H64lsmPHjjRsWPi2PnbsmG9u21NPPZVp06ZRp06dgOdfddVVjB07FnCmGPvf//4X8nrh5ObmMmnSJMD5Oa+55ppCx2+55RbffjRz9obzySefsGrVKgC6dOnCiy++SFxc8S/fRo0axe233w44Ler/+te/QtabkJDAtGnTaNCgQbFjY8eOJTU1FYBvv/2WTp068cILLxAbW3jMblxcXKG5hj/77LNidXk/84YNG3L55ZcHjUdE6NOnT8iYjTGmPLiR9K7BmXqsh4gU/d7Pf3BasJEMTT2Pe1yIpbJq77e/KYLy/mXaBy1VQiLSXER6lWSjcAJftW1e6HRVKI3da2DzV+7EU8nVqlWL005zPn5vVwAvbzIMsHz58qjr/vzzz9m5cycAo0ePLvZVelH+ieisWbOivp6/jz76iD17nF8/V111FSkpKYWOd+jQgbPPPhuAZcuW8f3335fqel7Tp0/37d9///3FEk9/f/jDHxCRYucFctlll9GqVauAxxITE+nWrZvv+R133BEw0Qbo0aOH7w+hNWvWFDvu/cz37dvH5s2bQ8ZkjDGVgRt9ej8DzseZv/YlERmpqlkiMhQYgtOsvVlVfwhy/hmeMhuCHK8O6vjt7w1WyM++IOcCICKXAn/2e6me5/EDETnm2f9EVR8NUv+twEMRxBFWVlYWmZmZEZfPz8+noKAAVSU/P9+NEEosZsnLiAv1FCx5GW1+bviCldyBAwd45513mDVrFqtXr2bfvn0cOXLE262lkG3bthX6/Pr374+IoKrceeedrFu3jmHDhtGhQ4eIru0/M8KhQ4fCJna5ubm+/dWrV0d0LxUUFAR83X9u3ptuuilgXcOHD/fNevDKK6/wf//3fwHr8n+vwsXk30J9/vnnhyzfrFkzOnTowNq1a9m6dSvbtm3jpJNO8h33/9nOPvvskHWlpaX59rt16xa0rIhQv359du3axYEDB4qVu+CCC5g+fToFBQX069ePBx54gMGDBxeq35ResPvWmLKmqhQUFJCfnx/V//NAmcyw4wY3kt6JwB9xkrNrgMtFJAtn2i3BSWifD3SiiHTDSdgUiL55qOrwbzrKjqD8Ub/91ADHGwLnBHi9q9/+jwGOG69jh+HHT1ypSn78GD12GGpU3cX0PvroI371q1+xb9++8IVxpibzd+qpp/LHP/6Rxx9/nCNHjvDYY4/x2GOPcdJJJ9GzZ0969erFoEGDaN8+8BcXW7Zs8e17uy5EqjSDpLZv387nnztjRZs1a8aAAYEHJl533XXcf//95OTk8M477/Dkk0+GbY0Ox9uynZqaSuPGjcOWb9u2LWvXrvWd65/0+qtfv37IemrUqBF12ezs4r+2Ro4cybRp05gzZw5btmzhrrvu4q677qJDhw707NmT3r17c8kllwTsZmGMMRWh1Emvqu4TkVtxFl6Iw+mjm+hX5Bvg2SCnD/PbX1jaWE4UqjoR54+NCpeUlFTs6+BQYmNjUVVEJOTXuWUuMx3UnZZm0Xxij+yGpLqu1FfeFi9ezNChQ8nLywOc/qUXXHABbdq0oW7dutSoUcP31bp3doKCgoJin99jjz3GOeecwxNPPMGiRU4X/p07dzJ9+nRfy+15553H008/zTnnFP6b7dChQyWOPzc3N6p7yb/spEmTfC2YN910U9B+zd5+q++//z779u3jv//9L0OHDi1Wzvs+Fb1OIBkZGYDTTSCS+L19ccFpRfE/JybmeE+1uLi4kPX5xxgfHx/xe1e0XGxsLJ999hnPP/88zz33HD///DMAP/74Iz/++COvv/46cXFxDB06lKeeeipokm4iU6G/L80JSUSIiYkhNjY2qv/nwckNKiM3WnpR1Rkici4wDmd0fy2cwVhTgMdUNa/oOSKSDHhHvBwjwGpu1UgmTss3OH8QhPuewH+p44wyiOc14IsSnnsa8IKLsVSMnBIOXgvmWHRf/VQmf/nLX3wJ7/PPP89dd90VtOzjjz8esq7LL7+cyy+/nPT0dBYuXMjixYuZP38+3377LarK119/Te/evZk5cyYXXHCB7zz/X6jff/+9r+9wWVLVQgPTnnjiCZ544omIzn3ttdcCJr3RSE1N5eDBg77pysLx/3rRPwGuSPHx8YwePZrRo0fz008/8fXXX7No0SLmzp3Lxo0bycvL45133mHevHksXbqUJk2aVHTIxpgTmCtJL4CqLgOujOKUXJz+vAC5qno0VOEq7iDHk94GhE96/b9zPOh2MKq6Fdha9HURGQOMCXN66b7TrSwSksOXiUaN6P4Krixyc3N9/Wm7desWMuEFIh6wlJaWxjXXXOObCWHr1q387ne/Y+rUqeTm5nLffffxww/Hu/k3a9bMt//LL7+US9I7b948Nm7cWKJzZ8+ezS+//MLJJ59c4uufdNJJHDx4kIyMDNLT08P2hV23bp1vvzImj+3bt6d9+/bceuutACxZsoTbbruNH374gR07dvD3v/+df//73xUcpTHmROZa0hstdZYs3hK2YPXwE+AdTt0K2BymvP/Q65/KIqAganF8No3qrVYTkFh3ujjExEFq1fzqdu/evb5W3jZt2oQsu3TpUvbujWQcZnHNmzfn7bffZt68eezZs4dVq1Zx8OBB37Rkffv25bnnngOcJXEvueSSEl0nGv4D2K6++mo6d+4c9pxFixYxe/ZsCgoKmDhxIn/+858LHffvZuDtxhPMOeec4+uj+/nnnzN8+PCgZbdu3cqPPzrd9Js3bx5RH+CKdvbZZ/Pmm29yxhlO28bChdaDzRhTsSos6T3B/ABc5NnvDswNVlBE0gBv89FuVS3PqdwOA9vDlEnAGUhXtSXWglMvgzUzSl9Xh8uc+qog/6nG/OfaDeShh0o34Ud8fDxNmzb1TQ/mTbbBmSe4YcOG7Nmzh9dee43f/va3YZPw0jh48CDvv/8+4PSB/c9//lNs7uBAvv/+e04//XQAXn/9df70pz8VSmz9u2kcOXIkZD+4q6++mokTJwLw9NNPc8MNNwTtt/nkk0/6Zoa4+uqrw8ZZWfhPneb/eRtjissvKCA3XylQJUaE+FghNsaNmWWNV5m9myLSWUQuFpEbROTmsrpOFeE/s/vFYcr6N3HNLINYglLV8araLNQGXFWeMZWp7rdVrnoqQK1atXzL+C5fvpxp06YVK5Ofn899993Hp59+GrSeZ599lqlTp5KTE3yNma+//to3x22zZs0KjepPTk7m4YcfBpxBWoMGDeK7774LGfuGDRsYM2YMu3fvLnZs4sSJiAgiQr9+/Yodf+edd3wzEngT7kh06dKFrl27ArBp0ybmzi3896t/kvfttwFXG/e55JJLfN04Vq5cyZ133hkwMZw4cSIvvvgi4AwO+e1vfxtRrGVtzJgxvgGLwbzwwvHu/973zRhznKqSmZ3Lln1HWLMjg3XpGWzYncm69AzW7Mhgy74jZGbnBpw60kTP1ZZeEWkBjAWuB4ouVPFmkbJpOLM6CLBMVf/hZiyVzHycpYIbA/1E5ExVLfY/oojEAvf6vTS5nOLzXv/E6dML0LI3NDy1dAtUNOoILXu5F1MFGD16tK8v79ChQ7nuuuvo27cvdevWZcOGDbz99tusXbuWzp07U6NGjYCLT3z77be88cYb1K5dm0GDBnHmmWfStGlTEhISSE9PZ/78+Xz00Ue+OUf/+Mc/FqvjrrvuYvny5bz22mts3LiRbt26MWjQIM4//3yaNWuGiLB//37Wrl3LwoULWbFiBeAkX9Hy79pw883R/U1+8803+6796quvFprm7IILLuDZZ53JakaNGsVvf/tbWrVq5WvBbdq0qS/RjYmJ4a233uLcc8/lyJEjvPzyyyxevJjhw4fTsmVL9u/fz4wZMwqthvbss8/SokWLqH/esjB9+nT+7//+jxYtWjBw4EC6dOlCw4YNyc/PZ/v27Xz00Ud8/fXXgNPK//vf/76CIzamcjmak8cvB46SnRu4m52iHDqay6GjuSTGx3Jy3ZrUTLAv6EvDtXdPRK4HJgDJUGzO/2J/oqhquog0AvoCl4jIC6padYfAh6Cq+SLyCMdnPXhTRAaoatEmqic4Ptfu16pauqWmonfi9OkFEIErX4TXLynZUsTxyTDkP049Vdgdd9zBsmXLeO2111BVJk+ezOTJhf/eOu2005gxYwYjR44MWIf3K/5Dhw4xZcoUpkyZErBcfHw8Dz30EHfeeWfA46+88grt27fnr3/9K1lZWXz22WcBl8D1atCgAYmJiUGPB7JixQpfK2zdunVDLqEbyI033sjYsWPJy8tj+vTphfomX3LJJfTt25f58+ezYcMGfvOb3xQ695ZbbvF1aQCn5Xju3LlcddVVbNu2jVWrVvHAAw8Uu2ZSUhLPPvts2CWey5P3M9+yZQuvvPJK0HL169dn0qRJvm4hxhjIyM5ly74sCiJswc3OzefnPUdoUT+J1MTgS8ab0FxJekXkauAt71OcGQcWA6cA7UKc+gpO0lsTp89r8e9WK5CItOL4tGpeXfz2zxCRx4ocn6OqgaZfexlndouBQCdgpYi8jLOMcz2c1nFvk+FB4PbSRV8iJ06fXq8mXWHYWzD5pugS3/hk57wmXcsqsnIjIrz66qtceumlvPTSSyxbtozDhw9Tv3592rdvz7XXXsuoUaNCJpf/+c9/GDZsGHPnzmXZsmWsW7eOPXv2kJeXR61atWjbti39+vVj1KhRtG3bNmQsY8eOZeTIkbz22mt88cUXrFmzxrdoRp06dWjTpg1nnXUWAwcO5MILLww6t24w/q281113XaHFGiLRqFEjLrroIj7++GOys7N55513fC3lsbGxzJo1i+eee44PPviAtWvXcvjw4ZD9Wbt37866det45ZVXmDFjBqtWrWL//v2kpKTQunVrBg0axN13313pZmxYvnw5s2bNYuHChXz33Xds3LiRAwcOICLUq1ePTp06cfHFF3PrrbdSt27VnMPamLJwNCcvqoTXq0CVLfuyOKVhsrX4lpCUtp+IiNQBNuKsyFYAPAo8oarHROTfwN2AqmqxERoikoqz5G4s8JKqBm7+qSAi0o8Qg86C+KuqPhykvlTgHeCyEOdvA65T1dCd5SqIiPTCbyGRhQsX0qtX5F/vr1+/nry8POLi4kImP+Vuxwr44I7Iujo06ui08FaDhPdE4l2Ewib5N1WJ3bfVi6qyfndm0C4NkUiMj6Vto5SQs8O4oTT/X3/11Vf07t3b/6XeqvqVqwGWgBt/KvwaJ+FVnITv0UhPVNUMEVmLs+BBtf/uS1UzcJZpHgzcjDOTQyOcBSh+BqYDE1S15MtTmZJp0hXuWgybv4KlL8PajwtPZxYT58zS0P02pw9vFe/SYIwxpvwdOZZXqoQXnK4OR47lk5Jorb3RcuMd8842sA94sgTn/4ST9LZ2IRZXqeo8ivdPdqPeGYALc2W564QbyFaUCLTq7WzZhyFjp7PSWo0UZx7eKjotmTHGmMph35HgM9xEV88xS3pLwI13rD1OK+9Cz4IT0TrgeSw624MpfyfWQLZQEmtZkmuMMcY1+QUFHD7qznzVh4/mkV9QYPP4RsmNpLee57H4ZJmR8XZUKnAhFlM6J95ANmOMMaYc5OYrWnwyqxJRlNx8JdZy3qi4kfQeAuoDwZceCs3bsrjPhVhMKajqeGB8qDJFB7IZY4wxJrxoZ2so7/pOBG78jbAdp99r1APRRCQe6InTPWKdC7EYY4wxxlQ6MS4PgHa7vhOBG0mvd0qvTiLSJWTJ4kbi9CMFCDS3rTHGGGNMlRcfK4hLY+MFIT7Wkt5oudG94R1gtGf/RRHpr6rHwp0kIp0B79LDeRxf3MJUkBN+9gZjjDGmjMTGxFCrZhyHjuaWuq5aNeNsEFsJlPodU9VlOPPLCnAO8KUnoQ1IRGqKyD3AVzitvAq8oqpbSxuLKTXv7A2hNhvEZowxxpRA/WR32o3qJ0e3kqRxuDXJ221AR6ADTh/dlSKyBmd5YQBE5AOgMdAVp7XQ2y7/HeFbF035sNkbjDHGmDKSXCOOxPjYUq/IllzDVugrCVeSXlU9KCL9gXeBfp6XO3oPex6v8Dz6d0KZAwyNpDuEKXs2e4MxxhhTdkSEk+vW5Oc9R0o0+0KM5/yyXoK4unKtQ4iqpgPnA7cAK3CS20AbwFpgBHChqu53KwZjjDHGmMqsZkIcLeonRT37QowILeonUTPBVmIrKVffOVVVYBIwSUQa43R1aIKz2toRIB34n6pucvO6xhhjjDFVRWpiPKc0TOaXA0cj6uqQGB/LyXVrWsJbSqV+90TkCr+nn6pqLoCq7gI+KG39pvzY7A3GGGNM+aiZEEfbRikcOZbPviPHOHw0r9CKbYJQq2Yc9ZNrkFwj1ro0uMCNPxk+xOm3+6OqfuRCfabieGdvMMYYY0wZExFSEuNISYwjv6CA3HylQJUYcebhtWnJ3OVG0psN1AC+daEuU7Fs9gZjjDGmAsTGxBBrOW6ZciPp3Qm0BGwGhirOZm8wxhhjTHXlxt8UP+DMynCKC3UZY4wxxhjjOjeS3nc9j+eKyMku1GeMMcYYY4yr3Eh6pwJfA/HARBGxtfGMMcYYY0ylUuqkV1ULgKHAcpzV2L4RkQvF5tYwxphqqWXLlogILVu2rOhQjDEmYm7M0/uaZ3c9cDrQBfgUOCAiK4A9wNEIqlJVHVXaeIwxxhhjjCnKjdkbRgBFF5AWoB7QP8q6LOmtQLY4hTHGGGOqK7fWswvWlSGaLg5FE2dT/mxxCmOMMcZUS24kvSNdqMNUDrY4hTHGGGOqpVInvar6hhuBmIpni1Mcl5mTSXpWOlm5WSTFJ5GWlEZKQkpFh2WMMcaYErIF74zxUFWW7FzCmHlj6DW5F0NmDOGGmTcwZMYQek3uxZh5Y1iycwmq1asnzrx58xARRISHH34YgB07dvDHP/6RTp06kZKSQq1atTjjjDN45JFHyMjICFqXt55+/fqFvW6/fv185QN5+OGHfcfnzZsHwPz587nuuuto0aIFNWvWpFWrVgwfPpy1a9cWOregoID33nuPgQMH0rRpUxITE2nbti2jR49m9+7dIeMqOjNBdnY2//rXv+jZsycNGzakZs2atGnThrvuuov169cHrGPVqlW+2AcOHBj2vfD+bN5zhg0bFra8m5+bv23btnHffffRvn17kpKSaNCgAeeccw5PP/00R49GMib5eD0vvPACw4YNo2PHjqSmphIfH++r78EHH+SXX34JW0/Rz+PYsWM8//zz9OvXj5NOOonY2NhCs0gUva8KCgp48803ueiii2jWrBnx8fHF7rm8vDxmz57N2LFj6du3LyeddBIJCQkkJyfTsmVLrr32WqZNm0ZBQUHAGPft20diYiIiQrt27SJ6fzZt2kRMTAwiQo8ePSI6xxhTSqpqm20Rb0AvnP7XCujChQs1GuvWrdM1a9bounXrojqvrK3eu1qHfDhEO0/sHHYb8uEQXb13dUWH7Jq5c+f6Ps+HHnpIZ82apfXq1VP/z9l/a9eunW7bti1gXd4yffv2DXvdvn37+soH8tBDD/mOz507Vx988EEVkYAx1axZU2fPnq2qqocPH9ZLL700aPxNmjTRn3/+OWhcLVq0UEBbtGihv/zyi55++ulB60pMTNSJEycGrKd3794KqIiEvJ7X9ddf76t3zpw5Ycu7+bl5ffzxx5qamhq0js6dO+vWrVsLvUfBYgv2WflvCQkJ+sorr4SMyf9amzZt0s6dOxerxz8O//tq//792qdPn4DX9te/f/+wsQJ67rnn6q5duwLGOXz48Kg+vwcffNBX/rXXXgtaLi8vT/Py8sLWZ4zbSvP/9cKFC4v+++mllSCHcWsgmzFV1qIdixg9dzRH8yJrxdpwcAMjPhvBM/2f4dwm55ZxdOVrxYoVPPXUU+Tm5jJixAh69epFamoqP/30Ey+88AK7du1i3bp1jBw5ks8//7zc4nrhhReYOnUqzZs3Z+TIkXTo0IHMzEymTZvGrFmzOHr0KNdeey2bNm3i5ptv5pNPPqFHjx4MHTqUpk2bsmPHDl566SXWrl3Ljh07GDFiBAsWLAh5zdzcXK699lpWrlxJ165dufHGG2nevDnp6elMmzaNBQsWkJ2dza233kqdOnUYPHhwofPvvPNOFi5ciKry8ssv8/e//z3otfbu3cv06dMBaNeuHf37RzfxjRuf2zfffMNVV11FTk4OAGeddRY33HADTZo0YceOHbzzzjssW7aMa6+9ltzc3JDxZGdno6q0b9+e/v3707FjRxo0aEBcXBy7du1iwYIFfPjhh+Tk5PCrX/2KtLQ0LrvsspB1Hjt2jKuuuopVq1bRo0cPrrnmGpo1a8b+/ftZvXp1wHNuvPFGFixYQKdOnbj++us55ZRTyMjIYP78+YXKZWVlkZycTL9+/ejWrRutWrUiNTWVI0eOsHbtWqZOncrPP//MokWLuPLKK1mwYAFxcYX/+7zzzjuZNGkSAC+99FLIzzA3N5fXX38dgNq1a3PdddeF/NmNMS6p6Kzbtqq1Uc1aelfvXa3d3+oeUQtv0a37W92rRYuvf4shnpbQVatWFSu3c+dObdasma/c8uXLi5XxHnO7pRfQiy66SI8cOVKs3MiRI31lunXrpoA++uijxcodPHhQO3bs6Cu7ZMmSgNf1tix6t3vuuSdgS9uTTz7pK9OoUSM9dOhQoePHjh3TRo0aKaCNGzfWnJycoO/FU0895avrqaeeClrOn5ufW15enp566qm+Mvfdd5/m5+cXKpOfn6/33Xdf0BZWf5s3b9YVK1aEjP+7777zvT9t27bVgoKCgOWKfh5PPPFEyHr97ytA77777rAtpbNnzw54b3nl5ubq3Xff7atz0qRJAct5vxWoUaOG7tmzJ2h906ZNK3R/hWItvaaiVMeW3lL36RWR11zaXi1tLMZEQ1UZ99W4iFt4izqad5RxX43z/jFQbbz11lt06tSp2OuNGzdm3LhxvueffvppucXUsGFD3n33XZKSkoode+SRR3x9NJcvX85FF13En/70p2LlUlJSeOCBB3zPP/vss7DXPeuss/jXv/5FbGxssWNjx47lyiuvBGD37t1MnDix0PGEhARGjXKmHt+1axcfffRR0Ou89NJLANSoUYNbbrklbFyBlOZz++STT3z9os8++2yefvppYmIK//cQExPD008/zdlnnx02lhYtWnD66aeHLNO1a1f+9re/AbB+/XoWLVoUtt7BgwcX+gzDOfPMM3n22WcDfn7+LrjggoD3lldcXBzPPPOMr+/wG28EHr991113AU6rdLAyABMmTPDt33777SFjM8a4x42BbCOAW1zYRrgQizERW7prKRsObihVHRsObmBZ+jKXIqp4Xbt2Dfm1rP+grFWrVpVHSAAMHz6cOnXqBDzWrFkzWrRo4Xt+zz33BK2nV69evv01a9aEve7vfve7Ysmfv7Fjx/r2p02bVuz4r3/9a9/5/omOv3nz5rFu3ToArr76aho0aBA2rqJK+7l5u1YA3HfffUEHF4oI999/f9TxBeP/eXzzzTdhy997771R1X/33XeH/PyiERcX5xtwtmRJ4AGtN954I7Vq1QKO/yFT1MaNG/niiy8AOPfcc+ncubMr8RljwnNr9gaJcgt0nqlgIjJGRLaF2oDpYSuqIib/NNmden50p57KoGfPniGPN2vWzLd/4MCBsg7HJ9zo9saNG/v2Q7VE+peLJP4LLrgg5PFzzjmH1NRUwGllLjq6v2XLllx88cUAfPHFF2zatKlYHW60+pX2c1uyZIlv//zzzw9ZV7jj/lasWMHvfvc7zjvvPBo1auSb4cC7dejQwVd227ZtIeuKjY3l3HOj60Pfu3fviMtmZWXx+uuvc9VVV9G2bVtq1arlm13Bu02e7PxbP3z4MIcPHy5WR3JyMsOHDwdg3bp1vllH/L388su+hNlaeY0pX+W5OEUMUBs4DbgcaABkA48AO12Iw5TeCbMiW2ZOJnO2znGlri+3fklmTma1mMc3XCtjjRo1fPvZ2dllHY5P/fr1Qx73jytU2Wjir1u3btjriginnHIKK1asICsri4MHD1KvXr1CZe666y4++eQTVJ0Bbd6v9MEZwPbBBx8A0KFDB/r06RPyesGU9nPbsWMHALVq1aJhw9Brz9SvX586depw8ODBoGXy8vK4++67CyV44QRKIoteNzExMaK6vPyT/VAWLVrEsGHDIppCzevw4cPUrl272Ot33nknzz//POC09vpP3+c/gK1u3boMHTo04usZY0qvQhanEJEawKPA74C7gQtVdW3os0w5OGFWZEvPSidf812pK1/z2Z21u1okvW59Fey2aOJy62dITk6OulxGRkaxpPeiiy6iZcuWbN68mddff51HHnnEN/J/4sSJHDt2DChdq19pf+bMzEyAkP1a/SUnJ4dMen/729/6vt6Pj4/noosu4uyzz6ZZs2YkJycTHx8POH2hvT93fn7of481a9aMKLZoz9m0aRODBg3yvQdt2rThoosuol27djRo0MDXOg3w7LPPMnfu3JDxdurUiT59+rBgwQKmT5/Ovn37fH88zZgxg/T0dABuvvnmqJN4Y0zpVMiUZap6DBgrzm+S+4FpItJdVbMqIh7j0BNoRbasXHdvtSO5R1yt70QQLsmpaEeORPaZ+pfzdnXwFxMTw+23386DDz7oG9B21VVXAc5X3QCJiYncfPPNLkRdMikpKRw6dIisrMj+XYR6b3755RdefPFFAJo2bcrcuXNp27ZtwLLBphorT3/72998Ce8DDzzA3//+96B9mt9+++2I6rzzzjtZsGCBb0DbmDFjgML9fH/961+XMnJjTLQqulnnL8BBoANwa8WGYk4kSfGRtWhFKjk+slbB6i4hIQHAN9drKHv37i3rcErlwIED7N+/P2QZVWXjxo2A00oabLDdqFGjfO+Ntw/v3LlzfQPYrr322mItxOWpaVOnV9Phw4fZs2dPyLL79u0L2cr7xRdf+Po2/+EPfwia8AIB+ziXN++8xY0aNeLxxx8PmvBC5PFeffXVpKWlAccTXf8BbL1796Zjx46lCdsYUwIVmvSq6lFgAc5AtusrMhZzYklLSiNWQk9jFKk4iaNRUiNX6qrq6tatC8D27aF7yezbt8+X8FVms2fPDnl8yZIlvr6oZ511VtBuBg0bNuSaa67x1bl58+ZCrX4VPaDJf/DfnDmh+7p/+eWXIY/v2rXLt9+mTZuQZctz2rtgvPG2atUq5NRmO3fuZOXKlRHVGR8f75uu7qeffmL+/Pk2gM2YSqCiW3oBvL8hQ/92NMZFKQkpDGg+wJW6BjQfUC3687rBO0/s1q1bWb9+fdByzzzzTLGZDiqj8ePHhxyI9dRTT/n2vUltMHfeeSfgtA7//e9/900T1qlTJ8477zwXoi0573zD4Hw2wX5mVeX//u//Qtbl38d5w4bgUwJu3Lgx5Fy25cUb788//xzys37kkUfIy8uLuF7/6eqef/553wC2+vXrh71XjDFlozIkvd5hx7UqNApzwhnWfpg79XRwp57qwDs9F8D9998fsN/utGnTeOKJJ8ozrBJbsmQJ9913X8AEffz48b65eRs1ahR2UYlevXpx2mmnAc5X3t4uIIFa/SZOnOibJst/9H9ZufTSSzn11FMBZ77csWPHFksAVZWxY8eGnU+3e/fuvv2nnnqKffv2FSuzdetWLr/88oj7TZclb7x79+7l6aefDljm6aef9vVTjlSLFi249NJLAZg6dapvANstt9xSaDYNY0z5qZCBbF4iUgfwTvq4uwJDMSeg7o2706ZOm1ItUNGmThvOSjvLxaiqtltvvZUnn3ySvXv38t///peePXty8803k5aWRnp6Oh9//DGzZs3i1FNPJTExke+++66iQw6qSZMmNG/enH/9618sWLCAG2+8kZNPPpndu3czbdo05s+fDzjTlr300ku+RQlCufPOO32rdoEzu4B3XteKFBsby2uvvUbfvn3JycnhqaeeYv78+dxwww00adKEHTt28M4777B06VLOOeccfvnlF980Z0X17NmTc845h//9739s2bKFDh068Otf/5pTTz2V/Px8vvnmGyZNmsSRI0cYMWJEsZXsytvo0aN9/Xp///vfM3fuXC666CLS0tLYunUrU6ZMYenSpZx00kmcdtppvrKRuPPOO/nvf/9b6DUbwGZMxamwpFdE6gHv4szdq1SDGQFM1SIiPN7rcUZ8NqJESxHXjKvJ471CD3w50dSrV4/JkyczePBgjhw5wtKlS1m6dGmhMp07d+ajjz5i5MhIp/iuGPHx8UybNo1LL72U7777LmCCXqNGDV588UUGDx4cUZ033XQTDzzwABkZGQBcd911QQe/lbcePXrw/vvvc8MNN5CRkRHwszvttNOYOnVq2EUfJk+ezIABA9i0aRN79+4tNDex129+8xvuu+++Ck96L774Yh5++GEefvhhAGbOnMnMmTMLlWnRogXTp0/n2WefjaruQYMG0bp1a99gx379+tG+fXtX4jbGRK/USa+IRDPPThxQH+gKXAb4d4T8d2ljMSZaHet35Jn+zzB67uioEt+acTV5pv8zdKxvI7CLOv/88/nhhx948sknmT17Ntu3b6dmzZq0a9eO66+/njvuuKPKzE/atGlTvvnmGyZMmMDkyZNZv349mZmZNG3alAsvvJAxY8aEnJ2gqNTUVHr06OEbIFfZBjRddtllrFmzhqeeeoqZM2fyyy+/kJSUxCmnnMJ1113HnXfeGdFcvi1btuS7777jmWeeYfr06b6+vY0bN+bcc89l1KhR9OvXj82bN5fxTxSZhx56iD59+vDss8+yePFi9u/fT506dWjdujVXXnklt99+e4n+OImJieH888/3Jb2V7fM25kQjka6WE7QCkQKcltqoT/Xbf0hVHy1VIKZcFJ2nd+HChfTq1Svi89evX09eXh5xcXFRJQtlbc2+NYz7alxEXR3a1GnD470et4S3ivH2Lw41Qh+chG3Lli20aNHC9aRs9+7dnHzyyeTk5NClS5eIZwMwVVNOTg7NmjVjz549NGzYkG3btvmmrotUpPetMW4rzf/XX331VdFvhHqr6leuBlgCbnVvKOn3uz8DD6rqNJfiMKZEOtbvyPQrprMsfRnv/vguc7bOKbRiW5zEMaD5AIZ1GMZZaWdZlwZTIq+++qpvAJt3NgdTfb3//vu+eY9vvfXWqBNeY4y73Eh6o5lzJhdnqdvNwP9UdWno4saUHxGhe+PudG/cncycTHZn7eZI7hGS45NplNTIpiUzpbJ//37fdF/16tWrFAPYTNk5duwYjz7qfIEZHx9faACjMaZilDrpVdXKPRrFmBJISUixJNeU2vz588nKymLbtm3861//8rX6Pfjgg4XmszXVw5IlS9i/fz+7d+9mwoQJrF27FnBmbGjevHkFR2eMqdApy0zlIiJjgDFhitn3c8ZE6JZbbmHLli2FXuvVqxejR4+umIBMmRo7dqxvKjuvdu3aBZy9whhT/izpNf5qAU0rOghjqpvExERatWrFsGHDuP/++4mLs1+91Vl8fDwtWrTg8ssvZ9y4cRHN4WyMKXtuTFnmXaj9A1WNetoxEbkDGAqoqp4frrwpU4eB7WHKJAANyyEWYyqEmzM2VJYpuUz5mDdvXkWHYIwJwY3mhn44U5atKuH5p/jVYSqQqo4HxocqU3TKMmOMMcaYqiCmogMwxhhjjDGmrFWGpNcbQ37IUsYYY4wxxpRQZUh6T/I8ZlZoFMYYY4wxptqq0KRXRJoBF+H0591YkbEYY4wxxpjqK6qBbCLyWojD54c57i8eZ2qsc4CaOEnv3GhiMcYYY4wxJlLRzt4wgsCzLAjQwbOVRAYQ9XRnpuqJjY0lLy+P/Px8VBURqeiQjDHGGONHVcnPd4ZaxcbGVnA07ilJ9wYpsgV7PdJtEXC+qm4t2Y9gqpKEBGdBN1UlKyurgqMxxhhjTFFZWVmoOm2c3v+3q4NoW3pHFnkuwGs4rb9fAm9HUIcC2cA+YJWqpkcZg6nCatWqRUZGBgD79+8nKSnJWnuNMcaYSkJV2b9/v+95dVpRMKqkV1XfKPqaXz/eHwMdN2VDRBKA+4GbgNY4s18sBB5T1W8rMrZQUlJSEBFUlczMTLZt20a9evUs+TXGGGMqkPcb2P3795OZ6UyoJSKkpKRUcGTucWNFtjdxWm+/caEuEwFPwjsLZyW73cB/caZ+uxK4TEQuV9VZFRdhcDExMTRt2pTt27f7Et/MzExEpFr1GzKVj/erOvvjylQldt+a8uIda+MlIjRt2pSYmMowu607Sp30quoIF+Iw0XkAJ+FdClygqocBROR64B3gLRFpraoZFRdicKmpqYUSX3B+sefl5VVwZKY6KygoAKhWv8BN9Wf3rakI3oQ3NTW1okNxlRstvaYciUgcMNrz9C5vwgugqu+KyE3AJcCtwL/KP8LIpKam0q5dOzIzMzl8+DA5OTm+kaLGlIXqOBLZVH9235ryEhsbS0JCArVq1SIlJaVa/qHlStIrIr8GEoEDqjopivNuAuoBR1T1VTdiOQGch/OebVLVZQGOv4eT9A6mEie94LRc1KpVq1p1kjeVl7ePWnXqn2aqP7tvjXFPqdN4ETkbeBH4P+CUKE9vCzwDvCQip5c2FreJSKyIdBaRESLybxFZLCJZIqKebWIJ671CRKaKyGYRyRaR3SKySER+LyLhMsCunsflQY57B7FVuvfTGGOMMaaiuNHSe6Xf/utRnvsa8GfP/jXAShficdMU4Cq3KhORFJxp3a4ocqihZ+sJ/EZEhqpqsIGBLTyP24Ic975eT0RSVDWzNDEbY4wxxlQHbnTYOM/z+JOqbonmRE/5nzxPe7kQi9uKdqLaD6wvSUUiEgtM5XjCmw48BtwA3AN87Xn9ZGCmiJwapCrvd1xHghz3T3KrVw90Y4wxxpgScqOltz3OlGU/lPD87z11tHchFrctAdbidCVYrqqbRGQE0bdoA9wGXOTZXwMMKLIwx/Mi8hTO3Lt1gQlAn5IGbowxxhhjjnMj6a3jedwfqlAI3vPqlj4Ud6nq39yox9PK+5DfS8ODrET3AHA+Tr/d3iJyoap+XqSMtyU3Ocjl/Ec7VMopy4wxxhhjypsb3RuyPY8lHVrqPa86z1fVB2fxCID5wVZMU9V84Fm/l64PUMzbhaRZkGt5X99v/XmNMcYYYxxuJL17PI+dSni+97w9IUtVbRf77c8MU/bTIOd5rfA8dgty/pmex8o2KNAYY4wxpsK4kfQuBwToIiLtojlRRNrjTK2lVO8k7TS//aWhCqrqLuAXz9M0EWlYpMjXOF1CWonIWQGquM7zOKMkgRpjjDHGVEdu9On9DLgWJ/F9TkQuUtWCcCd5+rk+7zlPKdzCWd34D9LbFEH5TTizOHjP9bWCq2qeiDwDPAK8ICJFlyG+BNiLMx1cQCLSHGgezQ/g57TwRYwxxhhjKhc3kt53cBKwJjiDsKaLyK2qGnRgm4jUw5kBYQBOwpsOvOlCLJVVHb/9vRGU3xfkXK8ncd67fsB6EZkPNAZ6A7k4A+VCDWK7lcID60osKyvLt2KQMZVZVlZWRYdgTNTsvjVVUWW9b0ud9KrqMRH5Lc4ctACXA5tF5F1gLrARZ8aBFKAVTrI2jMID336rqkdLG0sl5v+zZgctdZz/e1Fsrl1VzRGRQcDvgJtw5v7NxOnS8EiwgXLGGGNMWTpyLI+9mcfIyikgKSGGBik1SK7hRvuaMaXnyp2oqtNFZCzwD89LKTjz0t4W5BTxngo8qKpTg5QzQahqDvA3z1ZhkpKSbE14U6XY/Wqqosp836oqizfuY9LiLXy+Jp38AvUdi40RBnVK46YeLejZuj4iEqImU10kJSVVdAgBufbnl6o+LSJrgeeAlhxPbIPZAtyjqp+4FUMllsnxeYgTKbxqWiA1/fbLYq7d14AvSnjuacALLsZijDGmilq1/RBjpqxgXXrg/9byC5SZP+xi5g+7aJeWwvihXenctHY5R2mMw9XvHFR1poi0Ba7GWX2sB5CG8xV9Bk7f3W9wBq1N98xLeyI4yPGktwHhk976Rc51lapuBbYWfV1ExgBjwpye4HY8xhhjqp6F6/dw+6TlZOVE9l/5uvRMhk5YzITh3ejdtujERMaUPdc72ngS2SmezTh+wunPjOdxc5jyrfz2fyqLgIKoBTQtx+sZY4ypglZtPxRVwuuVlZPP7ZOWM+X2ntbia8qdG/P0mvB+8NvvHqqgiKRxfLqy3apanot2HAa2h9mq8yIixhhjwlBVxkxZEXXC65WVk8/9U1aiquELG+MiS3rLx2d++4FWWfN3id9+uNXbXKWq41W1WagNuKo8YzLGGFO5LN64L2gf3kj9lJ7BNxuDzmxqTJmweUTKx3xgF85cuv1E5MxA04p5Fuy41++lyeUUn/f61qfXGGNMSG99s8W1enqeUj98QWNcUuqk17O6lys8A6yqHVXNF5FHOD7rwZsiMkBVdxcp+gTQ1bP/tarOKq8YPaxPrzHGmKAysnOZtTrdlbo+W72LjOxcUhPjXanPmHDcaOndjDPfbmkplazlWURaAaOKvNzFb/8MEXmsyPE5qjonQHUvA1cCA4FOwEoReRlYA9QDrgd6ecoeBG4vXfQl4u3TG0oCYMNujTHmBLTrUHaheXhLI79AST+cbUmvKTduJpnVccbpFsC4EMe7UDgJBsgDiiW9qponIlfjLNt8GU5Xhz8HqHMbcJ2qri5RxKWgquOB8aHKiEgvYGH5RGSMMaYyOVLCwWvBZB47UWYuNZWBG0nvViJr6Y0BauN8hY7nnBycvq4nBFXNAC4XkcHAzTgzOTTCmcP4Z2A6MEFVD1VclMYYY0xgyQmxrtaXUsPd+owJpdRJr6q2jKa8iJwMXAf8AagDvKaqj5Y2jrKgqvMogxZsVZ0BzHC73tKygWzGGGNCaVw7kdgYcaWLQ1yMkFYr0YWojIlMuU9Zpqq/qOpTOAO2tgAPi8hfyzsOE5B3IFuozfrzGmPMCSo1MZ5BndJcqWtQp8bWn9eUqwqbp1dVtwE34rSkjhORkIs2mHJhi1MYY4wJ6aYeLSpVPcZEqkIXp1DVb4DvcBLfOyoyFmOLUxhjjAmvZ+v6tEtLKVUd7dNS6dG6nksRGROZyrAi2/c4SW+fig7EGGOMMaGJCOOHdiWphIPakhJieXro6YhUx0mfTGVWGZJe73wlTSo0CmOMMcZEpHPT2kwY3i3qxDcpIZYJw7vRuWntMorMmOAqQ9Lb1fN4rCKDMM7sDSKyLdSGM62aMcaYE1zvtg2ZcnvPiLs6tE9LZcrtPend1sZDm4pRoSugichVwJk4c/auq8hYDGDLEBtjjIlC56a1mTW6D99s3M+kbzYza3V6oenM4mKEQZ0ac1OPFvRoXc+6NJgKVSFJr4i0BUZQeE5Ya0GseLYMsTHGmKiICD1PqU/PU+qTkZ1L+uFsMo/lk1IjlrRaiTYtmak0Sp30isjGKK9XF0jynu553Aj8u7SxmNKxZYiNMcaURmpivCW5ptJyo6W3JZEtQ+xV9LuNlcCVqnrUhViMMcYYY4wpxq3uDdF20tkP/A94B5isqvlhyhtjjDHGGFNibiS9raIomwscVtVMF65rXCYiYyjczzqQhPKIxRhjjDHGTaVOelV1ixuBmErBZm8wxhhjTLVUoVOWmUrHZm8wxhhjTLXkxuwN+z27B4H2qppb2jpNxbDZG4wxxhhTXbmxIlstoDawwhJeY4wxxhhTGbmR9O7xPO52oS5jjDHGGGNc50bSu9Xz2MCFuowxxhhjjHGdG0nvDJx5evuISKwL9RljjDHGGOMqN5Le13AGsdUHxrpQn6kgIjJGRLaF2oDpFR2nMcYYY0y03Jind5eIjASmAI+KiAL/UNWCUkdnypvN02uMMcaYasmNKcv6AAeAh4G/Ao8Dd4vINGA5zkC3o5HUpaoLShuPKRWbp9cYY4wx1ZIbi1PMA9TvueC0Ft4bZT3qUjymhGyeXmOMMcZUV24lmRLha8YYY4wxxpQ7N5LeBRRu6TXGGGOMMaZScWMgWz8X4jDGGGOMMabMuDFlmTHGGGOMMZWaJb3GGGOMMabac2PKsr94dpeo6mclOP8C4FwAVX2ktPEYY4wxxhhTlBsD2R7GGcj2PBB10gsMAu731GFJrzHGGGOMcZ3Ni2t8RGQMMCZMsYTyiMUYY4wxxk2VIen1zudr055VPFuG2BhjjDHVUmVIeut4HrMqMggD2DLExhhjjKmmKjTpFZF4oJ/nabhky5QxW4bYGGOMMdVVVEmviNwc4nCHMMf9xeN8jT4YaI3TteF/0cRijDHGGGNMpKJt6Z1I4L63Apzv2UpCgQklPNcYY4wxpkrLzMkkPSudrNwskuKTSEtKIyUhpaLDqlZK0r1Bonw9nKPAGFVdXMLzjTHGGGOqHFVl6a6lTP5pMnO2ziFf833HYiWWAc0HMKz9MLo37o5ISdMs4xVt0vtGgNduwWmp/YnIuigokA3sA1YBn6rq4SjjMMYYY4ypstbsW8O4r8ax4eCGgMfzNZ/ZW2Yze8ts2tRpw+O9Hqdj/Y7lHGX1ElXSq6oji74mIrd4dr9Q1XtdicoYY4wxpppatGMRo+eO5mje0YjKbzi4gRGfjeCZ/s9wbpNzyzi66ivGhTq2erZ9LtRljDHGGFNtrdm3JqqE1+to3lFGzx3Nmn1ryiiy6q/USa+qtlTVVqr6VzcCMsYYY4ypjlSVcV+Nizrh9Tqad5RxX41D1dbzKgk3WnqNMcYYY0wYS3ctDdqHN1IbDm5gWfoylyI6sZTr4hQi0g24HEgDduMMYvumPGMwxhhjjKkIk3+a7E49P06me+PurtR1Iil10isidYBnPU+/V9WngpR7BvhNkZf/JCJvAbeq+s3TYYwxxhhTjWTmZDJn6xxX6vpy65dk5mTaPL5RcqOl9wrgJpypyL4MVEBERgHBZna4CTgU4rgxxhhjTJWWnpVeaB7e0sjXfHZn7bakN0puJL39PI/5wIyiB0UkBnjY81SB1cDnwMnANTiLWtwpIi+o6o8uxGNKSETGAGPCFEsoj1iMMcaY6iQrN8vV+o7kHnG1vhOBG0lvF8/jWlU9GOD4+UBTnIT3K2CgquYAiMg9OF0jYoCbgT+6EI8puVo4n5UxxhhjXJQUn+Rqfcnxya7WdyJwI+k9GSehDdZKO8hv/3FvwuvxEk4rcF2glwuxmNI5DGwPUyYBaFgOsRhjjDHVRlpSGrES60oXhziJo1FSIxeiOrG4MWVZLc/j/iDH+3geD1Okz68nAf4Gp4tDOxdiMaWgquNVtVmoDbiqouM0xhhjqpqUhBQGNB/gSl0Dmg+w/rwl4EbS661Dih4QkZpAV5yW4K+DzNCwy/NY24VYjDHGGGMqpWHtr6tU9Zxo3Eh6D3seTwpw7DyOd6FY5MK1jDHGGGOqpO5Hs2mTkxO+YAhtcnI462i2SxGdWNxIejfgtPL2EJHYIseu9NtfGOR8b6eUgy7EYowxxhhTKcmyV3l8zz5qFhSU6PyaBQU8vmcfsuxVlyM7MbiR9C7wPDbAb7orEWkLDPc8PYTTdzeQLjjdHza6EIsxxhhTKhnZuaxPz2DFLwdZn55BRnZuRYdkqoPsw7D2Yzrm5PJM+t6oE9+aBQU8k76Xjjm5sPa/Tn0mKm7M3vAaTrIbAzwhIpcDe4D+QApOQvumqhb7rSEiLYDmnjIrXYjFGGOMiZqqsnjjPiYt3sLna9LJL1DfsdgYYVCnNG7q0YKeresjUmwIizHhHd4BnqFN52ZnM3FnOuMa1mdDQvjp79vk5PD4nn1OwgtOPRk7IbFW6BNNIaVOelX1JxF5HPgLTvJ6XpEi6cDjQU6/xm//69LGYowxxkRr1fZDjJmygnXpmQGP5xcoM3/YxcwfdtEuLYXxQ7vSuamNvTZRyim8mETHnFymb9/FssQavFsrlTlJNcn3+4MqTpUBR7IYlpHJWdnHis8WcCzw/WqCc6OlF1V9WEQOAONwujl4fQOMUtU9Rc8R50/lO7xV4KzSZowxxpSbhev3cPuk5WTlRDZ36rr0TIZOWMyE4d3o3damLDdRSCi+mIQA3bOP0T37GJki7I6L5YjEkKwFNMrLJ0W1eD1eNWzKsmi5kvQCqOq/ROQ5oAPO3L2/qOq2EKfUBh7z7GcGSoyNMcaYsrJq+6GoEl6vrJx8bp+0nCm397QWXxO5Wk1AYn1dHIpKUSUlNy+yumLiIDXQpFkmFDcGsvmoar6qrlbVxWESXlT1oKq+4dnedzOO6kxEzhSRsSIyVUQ2i4h6tpYVHZsxxlQVqsqYKSuiTni9snLyuX/KSjRUS5wx/hJrwamXuVNXh8usP28JuJr0mnLxF+BJnP7QLSo4FmOMqZIWb9wXtA9vpH5Kz+CbjcEWIzUmgO63Va56TjCW9FY9i3G6hQwBmuEMFDTGGBOFt77ZUqnqMSeIlr2h4amlq6NRR2jZy514TjCu9ek15UNVn/R/blPnGGNMdDKyc5m12p32gs9W7yIjO5fUxHhX6jPVnAhc+SK8fgnkHglfvqj4ZBjyH6ceEzVr6TXGGHNC2XUou9A8vKWRX6CkH7YlYU0UmnSFYW85CWw04pOd85p0LYuoTgjVMukVkVgR6SwiI0Tk3yKyWESy/AZ9TSxhvVf4DSDLFpHdIrJIRH4vItaj3BhjqoAjJRy8FkzmMXfrMyeAUwbAyJmRd3Vo1NEpf8qAso2rmquu3RumAFe5VZmIpABvA1cUOdTQs/UEfiMiQ1U12HLLxhhjKoHkhFhX60up4W595gTRpCvctRg2fwVLX4a1HxeeziwmzpmlofttTh9e69JQatU16S36G2g/sA9oG21FIhILTAUu8ryUDrwMrAHqAdfjrEJ3MjBTRM5T1bUljNsYY0wZa1w7kdgYcaWLQ1yMkFYr0YWozAlJBFr1drbsw87SwscynYUnUk+yaclcVl2T3iXAWmA5sFxVN4nICOD1EtR1G8cT3jXAAFX1HwHxvIg8BdwP1AUmAH0CVSQibwJnR3n9D1T1wSjPMcYYE0RqYjyDOqUx84ddpa5rUKfGNojNuCOxliW5ZaxaJr2q+jc36vG08j7k99LwIgmv1wPA+UBXoLeIXKiqgZZVbg60jzIMW3LFGGNcdlOPFq4kvTf1sOnSjakqquVANhf14XjSOV9Vvw1USFXzgWf9Xro+SLl+qipRbiNc/YmMMcbQs3V92qWllKqO9mmp9Ghdz6WIjDFlzZLe0C72258ZpuynQc4zxhhTyYgI44d2JamEg9qSEmJ5eujpNle6MVWIJb2hnea3vzRUQVXdBfzieZomIg3LLCpjjDGl1rlpbSYM7xZ14puUEMuE4d3o3LR2GUVmjCkLZdKnV0TqcnxGg7pAnKo+UhbXKmP+/W83RVB+E87P7D13j+sRuUBEmuP0Ly6J08IXMcaYqqF324ZMub0nY6asYF16Ztjy7dNSeXro6ZbwGlMFuZr0isg5wF+AQUDR73weKVI2DZiB09q8WFV/62YsLqnjt783gvL7gpzrGhG5FPiz30veDmUfiMgxz/4nqvpoiGpupfAAvRLLysoiMzP8fxTGVLSsrKyKDsFUUi1rx/L+bWey4pdDzFixnYXr95Kvx6czixWhd9sGDO7alK4n10ZEyu33nt23piqqrPeta0mviPwBJ7GNpXjCW2wyRFVNF5FdOAs+dBGRR1R1X9FyFcx/lEMk60we9dtPdTkWr4bAOQFe7+q3/2MZXdsYY6olEeGM5nU4o3kdjhzLY2/mMbJyCkhKiKFBSg2Sa1TLyY6MOaG48q9YRO4G/KcJ+xGYC3QHzgpx6us4SW88cAkwyY14qjNVnQhMrOAwfJKSkkhJKd0IaGPKk92vJpyUFEirX9FRFGb3ralKkpKSKjqEgEqd9IrIScA/PE+zgTtV9Q3PsX8TOun9DDgGJAD9qXxJbyZOn2SARM/zUGr67WeUSUTueA34IsDr13m2UOIpo64bxhhjjDFlxY2W3jtwkj0F7vMmvJFQ1WMisgroRuUcIHWQ40lvA8Invf5tAwfLIB5XqOpWYGvR10XkApzuE8YYY4wx1YobSe8gz+MO4KUSnP8zTtJbGZe1+Qlo5dlvBWwOU76V3/5PZRFQGTsMbA9TJgFLjI0xxhhTxbgxT29rnFber1W12IC1CBz0PFbGBad/8NvvHqqgZzYK73Rlu1W1Uk5XFoqqjlfVZqE24KqKjtMYY4wxJlpuJL3eyQr3l/D8Gp7HXBdicdtnfvvhVlm7xG8/3OptxhhjjDGmHLnRveEAztfddcMVDMK7SEIk8+CWt/nALqAx0E9EzlTVb4sWEpFY4F6/lyaXU3yuEpExwJgwxRLKIxZjjDHGGDe5kfRuARrh9MuNiogk48w5q8AaF2Jxlarmi8gjwAuel94UkQGqurtI0Sc4Pk/u16o6q7xidFktoGlFB2GMMcYY4zY3kt4vcPq7niIivVV1YRTn/gZIwkl6A02hVSIi0goYVeTlLn77Z4jIY0WOz1HVOQGqexm4EhgIdAJWisjLOEl6PeB6oJen7EHg9tJFX6FsIJsxxhhjqiU3kt5JwFic/sETRKSXqobt3ysi5wMPe55mAW+6EItXC2BciONdKJwEA+QBxZJeVc0TkauBd4DLcLo6/LloOWAbcJ2qri5RxJWAqo4HxocqIyK9gGj+sDHGGGOMqXClHsimqj8Cr+IsPdwe+EZELhORoksRAyAirUXkKZzBXgk4rbxPV8IliH1UNUNVLweGANOBX3AW1dgL/A94AOisqosqLEhjjDHGGBOUW4uJ3wt0Bs4FTgFm4HxV7puRQUS+w2klbeR9yfP4OfBXl+IAQFXn+dXvZr0zcH62askGshljjDGmunIl6VXVHBG5EPgPMNzzcm2cVlzv3L1dKJ6IvgHcXsL5fY37bCCbMcYYY6olt1p6UdUs4BYReQWntfACIDlA0WM4fWefVNUFbl3fuMIGshljjDGmWnIt6fXyzN6w0DN3bRegCU6r7xEgHVihqtluX9eUng1kM8YYY0x15XrS66Wq+cB3ns0YY4wxxpgKU2ZJr6l6bCCbMcYYY6qrUk9ZJiIFIpIvIs+W8Px/es7PK20sptS8A9lCbdaf1xhjjDFVTmVp6XV9ejFTIjaQzRhjjDHVUmVJek0lYAPZjDHGGFNdlbp7gwviPY+5IUsZY4wxxhhTQpUh6W3jeTxYkUEYY4wxxpjqq8K6N3jm8R0CDMRZtW1tRcVijDHGGGOqt6iSXhHZGOLwzSJyWYRVxeMMhorHGcSmwCfRxGLcZ1OWGWOMMaa6iraltyVOglqUAKmeLVL+Mzb8CLwQZSzGfd4py4wxxhhjqpWSdG8INr1YtNOOFQDrgfeBf6hqVgliMe6yKcuMMcYYUy1Fm/S2KvJcgI04rb9vAg9HUIcC2cBBVc2J8vqmDNmUZcYYY4yprqJKelV1S9HXRASc5Dcj0HFjjDHGGGMqmhuzN4z0PNrsC8YYY4wxplIqddKrqm+4EYgxxhhjjDFlpTIsTmGMMcYYY0yZqrDFKUzlY/P0GmOMMaa6cj3pFZH2wAVAV6ABzty9kbQoq6qe73Y8Jio2T68xxhhjqiXXkl4RaQ28CJQkcfWuymYqls3Ta4wxxphqyZWkV0S6AvNwWnUjWaTCm+BGu6CFKUM2T68xxhhjqqtSJ70iEg9Mx/lqHOBTYCJwPTAEJ8EdgJMQtwD6AINxWgwzgfuBdaWNwxhjjDHGmGDcaOkdDrTESW4nqeoIABHp4y2gqvP9yj8vIk2Al4BLgH8CF6nqNy7EYowxxiMjO5ddh7I5kpNPckIsjWsnkpoYX9FhGWNMhXAj6b3M85gH/C6SE1R1h4hcgdNCfAUwWUROV9VDLsRjjDEnLFVl8cZ9TFq8hc/XpJNfcHy4RGyMMKhTGjf1aEHP1vW9K2oaY8wJwY2k9wycVt4lqro30pNUtUBEbgcuAk4GbgGedSEeY4w5Ia3afogxU1awLj0z4PH8AmXmD7uY+cMu2qWlMH5oVzo3rV3OURpjTMVwY3GKBp7HDUVez/fuiEjNQCeqajqwAGdA2zUuxGKMMSekhev3MHTC4qAJb1Hr0jMZOmExC9fvKePIjDGmcnAj6fUuVpBV5PUMv/1GIc7f7Hls7UIsxhhzwlm1/RC3T1pOVk5++MJ+snLyuX3SclZtt55lxh2ZOZn8fPBnftjzAz8f/JnMnMj+CDOmPLjRveEAzrytyUVe928+aAdsCXK+NyFuEOS4McaYIFSVMVNWRJ3wemXl5HP/lJV8Nrq39fE1JaKqLN21lMk/TWbO1jnk6/F7MVZiGdB8AMPaD6N74+52j5kK5UbSux4ncW1W5PUf/PYvBGYXPVFEEoHunqeHXYjFlIItQ2xM1bN4476IuzQE81N6Bt9s3E/PU+q7FJU5UazZt4ZxX41jw8GiPRwd+ZrP7C2zmb1lNm3qtOHxXo/TsX7Hco7SGIcb3Ru+9Tx2KvL6Nxzv4vArz/LERT0KnIQzEO47F2IxpeNdhjjUZquxGVOJvPVNsC/RKqYec+JYtGMRIz4bETThLWrDwQ2M+GwEi3YsKuPIjAnMjaR3juexkYj4/nxT1aPAJJxBaqnAEhF5VkR+JSK/FZEvKdyq+IYLsZjS8S5DHGqzUS/GVBIZ2bnMWp3uSl2frd5FRnauK3WZ6m/NvjWMnjuao3lHozrvaN5RRs8dzZp9a8ooMmOCcyPpnQUcwUlury9y7M/AL55jKcDdwIs4S9328yv3uaq+40IsphRUdbyqNgu1AVdVdJzGGMeuQ9mF5uEtjfwCJf1wtit1mepNVRn31bioE16vo3lHGffVOFTduXeNiVSpk15VzQbOBy4G/lvk2AGgN7AYJ/EtugG8iSVSxhgTtSMlHLwWTOYxd+sz1dPSXUsj7tIQzIaDG1iWvsyliIyJjBsD2VDVJSGObQXOE5GeOMlxE6AA2Ah8oqo/uRGDMcacaJITYl2tL6WGu/WZ6mnyT5PdqefHyXRv3D18QWNc4krSGwlVXYzT4muMMcYFjWsnEhsjrnRxiIsR0moluhCVqc4yczKZs3VO+IIR+HLrl2TmZJKSkOJKfcaE40afXmOMMRUgNTGeQZ3SXKlrUKfGpCbGu1KXqb7Ss9ILzcNbGvmaz+6s3a7UZUwkLOk1xpgq7KYeLSpVPaZ6y8otuvhq6RzJPeJqfcaEUibdG0SkK9ALOBmoC8Sq6qiyuJYxxpzIerauT7u0lFItUNE+LZUereu5GJWprpLik1ytLzm+6GKuxpQdV1t6ReQaEVkFLAf+BfwOuA0YEaBsmojsEJHdIjLdzTiMMeZEISKMH9qVpBIOaktKiOXpoafb8rAmImlJacSKOwMe4ySORkmNXKnLmEi4lvSKyIvAe8CpFJ+WrBhVTQe+BBoAV4hIU7diMcaYE0nnprWZMLxb1IlvUkIsE4Z3o3PT2mUUmaluUhJSGNB8gCt1DWg+wAaxmXLlStIrIo8Bv+Z4ojsLeACYG+bUN71VAJe4EYsxxpyIerdtyJTbe9IuLbIkon1aKlNu70nvtrayuInOsNodwxcqx3qMiVSp+/SKSFtgrOfpQeAqVZ3nOdYc6B/i9Dk4q7kl4azQ9nJp4zHGmBNV56a1mTW6D99s3M+kbzYza3V6oenM4mKEQZ0ac1OPFvRoXc+6NJgS6b5uPm1yctiQkFDiOtrk5HDWuvlwhg33MeXHjYFsv/bUo8CvvAlvJFQ1X0S+B3oCnVyIxRhjTmgiQs9T6tPzlPpkZOeSfjibzGP5pNSIJa1Wok1LZkon+zDy4yc8Hh/DiJPSOBoT/RfGNQsKeHzPPmTHx5B9GBJrlUGgxhTnRtJ7vufxZ1V9vwTnb8ZJepu5EIspBREZA4wJU6zkf9obY8pVamK8JbnGXYd3gObTMSefZ9L3MjqtQVSJb82CAp5J30vHnFznhYydlvSacuNG0tsCp5X3fyU8/7DnMdWFWEzp1AJsQKExxpjAco7Pq3tudjYTd6YzrmH9iLo6tMnJ4fE9+44nvADHSj7VnjHRciPp9U6yV9I71zvpX7YLsZjSOQxsD1MmAbCRL8YYcyJKKDyvbsecXKZv38WyxBq8WyuVOUk1yffrKx6nyoAjWQzLyOSs7GPFp3SqYbM3mPLjRtK7D2hMyROhUzyPe1yIxZSCqo4HxocqIyK9gIXlE5ExxphKpVYTkFjwW4pYgO7Zx+iefYxMEXbHxXJEYkjWAhrl5ZOiGriumDhIPal84jYGd6Ys+xnnnj872hNFpD5wFk73iJUuxGKMMcaYspJYC069LOjhFFVa5+ZxWk4OrXPzgie8AB0us/68ply5kfR+7nlsKiKDozz3DxwfGDXbhViMMcYYU5a631a56jEmQm4kvRM53h/3BRFpGclJInILzkwBChwAJrkQizHGGGPKUsve0PDU0tXRqCO07OVOPMZEqNRJr6puA57C6eLQGFgqIvd4ui4UIiKJIjJARKYBr3F8meKHVPVI0fLGGGOMqWRE4MoXIT45fNlA4pNhyH+ceowpR64sQww8BEzHSWLrAf8CdgO+7y5E5ADO6muzgSs5nvC+oarPuxSHMcYYY8pak64w7K3oE9/4ZOe8Jl3LIipjQnIl6VVVBYYCjwEFOAmt4PTX9fZir+33ugD5wF9V9VY3YjDGGGNMOTplAIycGXlXh0YdnfKnDCjbuIwJwo0pywBQ1QLgLyLyKjAauBhoF6DoL8DHwNOqutGt6xtjjDGmnDXpCncths1fwdKXYe3HhaYzIybOmaWh+21OH17r0mAqkGtJr5eqbgHuA+4TkXrASTitvEeAdFXd5fY1jTHGGFNBRKBVb2fLPuwsLXws01l4IvUkm5bMVBqlTnpFxLuYQQHwoKr61hdU1f3A/tJewzhEJB7oD1wG9MNZ2CMW2AbMAv7h+aPDGGOMKX+JtSzJNZWWGy29o3H67X7tn/CaMtEXJ7kFp5uId47ks4G7gJtE5CJVXVwRwRljjDHGVFZuDGQ77Hlc50JdJrQCYBpwrqo2V9UrVfVKnBbfiUAt4F1Pi7AxxhhjjPFwo6V3J5AKWKJVxlR1DjAnwOvZInIXzlRwLYBzgfnlHJ4xxlQLmTmZpGelk5WbRVJ8EmlJaaQkpFR0WMaYUnIj6f0KaA+c7kJdpoRU9aiIrAO6A00qOh5jjKlKVJWlu5Yy+afJzNk6h3y/GQhiJZYBzQcwrP0wujfujtgMBMZUSW4kvROBUcBpInKuqi5yoc5SEZFY4FTgLKCb5/F0oKanyBuqOqIE9V4BDMdJLBvjdO3YAHwATFDVwyFOL1Oen7ml56nNkGGMMRFas28N474ax4aDGwIez9d8Zm+Zzewts2lTpw2P93qcjvU7lnOUxpjScmMZ4q+Bl3AWnHhbRE4pdVSlNwX4AXgduAfowfGEN2oikiIiM4AZwDU4XQhqAA2BnsA/gFUi0qOUcZfGcE88e4AK/8PDGGOqgkU7FjHisxFBE96iNhzcwIjPRrBoh/2aNaaqcWsZ4t8AL+AkgytE5G8icrqIuFV/tGKLPN8PrC9JRZ4W1KnAFZ6X0nFWnrsBJ6H+2vP6ycBMEYlwaRr3iEhL4GnP03Gqeqy8YzDGmKpmzb41jJ47mqN5R6M672jeUUbPHc2afWvKKDJjTFlwY55e/1XV8oFk4AHPlisiB4BIfqOoqrrVSrwEWAssB5ar6iYRGYHT8hut24CLPPtrgAGqmu53/HkReQq4H6gLTAD6BKpIRN7EmV4sGh+o6oPBDopILZwW6HrAVFV9Ocr6jTHmhKOqjPtqXNQJr9fRvKOM+2oc06+Ybn18jaki3OjT2xJnnl4v774ACUCjCOqQInWUiqr+zY16PK28D/m9NLxIwuv1AHA+0BXoLSIXqurnAco1xxn0F42TQsSXCHwEdAG+xOniYIwxJoylu5ZG3KUhmA0HN7AsfRndG3d3KSpjTFlyq/uBBNhCHQtWtrLpw/Gkc76qfhuokKrmA8/6vXR9kHL9VFWi3EYEqsszF+/7OAtWfAMMtm4NxhgTmck/TXannh/dqccYU/bcaOlt5UIdldXFfvszw5T9NMh5rvP0lX4LuARYCVyiqkfK8prGGFNdZOZkMmdrsSnPS+TLrV+SmZNp8/gaUwWUOulV1S1uBFJJnea3vzRUQVXdJSK/4AxoSxORhqq6x+2AxOk89jIwFPgJuFBVD7h9HWOMqa7Ss9ILzcNbGvmaz+6s3Zb0GlMFuNHSW53597/dFEH5TThJr/dc15NenFkabvVc63xV3R1tBSLSHKd/cUmcFr6IMcZUXlm5Wa7WdyTXvmgzpiqwpDe0On77eyMovy/Iua4QkcHAfZ6nm4HHg4wa/lBVPwxR1a0UHqBXYllZWWRmZrpRlTFlKivL3UTHVF2xebHUi6nnWn1xeXFl9nvQ7ltTFVXW+9aS3tD8v6/KjqC8/9w3qS7HAs6UaF79Q5TbDHxYBtc3xpgqr0FiA2Il1pUuDnESR/3E+i5EZYwpa27M05sIzAESgW9V9bYozn0VOAPIwPmqPq+08VRnqjoRZ9nnSiMpKYmUFOvLZqoOu19NCil0bdaV2Vtml7quC1tcSKO6kczMWTp235qqJCkpqaJDCMiNlt4hOMv8Ks5yvNH4AhjpOfcKYLoL8bgpk+Otq4me56H4L3WcUSYRueM1nPe+qOs8WyjxlEHXDWOMKU/D2g9zJekd1mGYC9EYY8qDG0mvd7Wyozgrg0XjQ895icClVL6k9yDHk94GhE96/b/jOlgG8bhCVbcCW4u+LiIXAA3LPyJjjClf3Rt3p02dNqVaoKJNnTaclXaWi1EZY8qSG4tTnInTUvudqka1nqOn/Lc4C1R0cyEWt/3ktx/JfMT+ZX4KWqryOgxsD7OVxYwUxhhTrkSEx3s9Ts24muELB1AzriaP93rcliA2pgpxI+lt4Xn8uYTne89rEbJUxfjBbz/kOpMiksbx6cp2l8UcvWVNVcerarNQG3BVRcdpjDFu6Fi/I8/0fybqxLdmXE2e6f8MHet3LKPIjDFlwY2k1/vbIpLZDQLxnpfsQixu+8xvP9wqa5f47Ydbvc0YY0wlcG6Tc5l40UTa1GkTUfk2ddow8aKJnNvk3DKOzBjjNjf69B7A6e9a0r6g3vMOuxCL2+YDu4DGQD8ROVNVvy1aSERigXv9XqqSi7GLyBhgTJhiCeURizHGlJeO9Tsy/YrpLEtfxrs/vsucrXMKTWcWJ3EMaD6AYR2GcVbaWdalwZgqyo2kdztO4lrSP3vPxekTvNOFWFylqvki8gjwguelN0VkQIBV0J4Aunr2v1bVWeUVo8tqAU0rOghjjClvIkL3xt3p3rg7mTmZ7M7azZHcIyTHJ9MoqZEtM2xMNeBG0rsAJ+FrJCJDVXVKpCeKyHVAGk7Su9CFWLz1tgJGFXm5i9/+GSLyWJHjc1R1ToDqXgauBAYCnYCVIvIysAaoB1wP9PKUPQjcXrroK5R3IFsoCdgMD8aYaiwlIcWSXGOqITeS3qkc/2r/3yLyraqGnQNGRNoB/y5Sj1taAONCHO9C4SQYIA9nkY1CVDVPRK4G3gEuw+nq8OcAdW4DrlPV1SWKuBJQ1fHA+FBlRKQXLv6BYowxxhhTHko9kE1Vv8ZZ6EBwWgCXiMivRKRGoPIiUkNEfgV8g9MXWIEFqjq3tLGUFVXNUNXLcRbimA78AhwD9gL/Ax4AOqvqogoL0hhjjDHGBOVGSy84q6otxemqUAd4EfiniHwDbMRZ1CEFZx7bnkAqTpIMsBu42aU4AFDVeX71u1nvDKJfgKPKsIFsxhhjjKmuXEl6VXW7iAwEPgC8877UwukHW5R/MvozcJWq/uJGHKbUbCCbMcYYY6olt1p6UdXVItIN56v+u3BafIO1th7AmRHhH6qa4VYMptRsIJsxxhhjqiXXkl5w+r4Cf/JM89XDs6XhdGfIANJx+vJ+o6o5bl7blJ4NZDPGGGNMdeVq0uvlSWgXeDZjjDHGGGMqVJkkvaZqsoFsxhhjjKmuLOk1/mwgmzHGGGOqJUt6jT8byGaMMcaYasmSXuNjA9mMMcYYU12VekU2Y4wxxhhjKjtLeo0xxhhjTLVnSa8xxhhjjKn2rE+v8bEpy4wxxhhTXVnSa/zZlGXGGGOMqZYs6TX+bMoyU+lkZOey61A2R3LySU6IpXHtRFIT4ys6LGOMMVWMJb3Gx6YsM5WFqrJ44z4mLd7C52vSyS9Q37HYGGFQpzRu6tGCnq3rIyIVGKkxxpiqwpJeY0ylsmr7IcZMWcG69MyAx/MLlJk/7GLmD7tol5bC+KFd6dy0djlHaYwxpqqx2RuMMZXGwvV7GDphcdCEt6h16ZkMnbCYhev3lHFkxhhjqjrXk14RSRaRX4nIOyLyo4jsFZFjnscfPa/fJiJJbl/bGFN1rdp+iNsnLScrJz+q8/6/vTuPk6sq8z/+edLZyEYgJB2CQIAgJCBGoBmCCYaggAyyqRAcGIPgII46EXRcMrKNOA5q5KeiZEB2NaDD4oJgJMAQAhLAyJIEiBKChHRCMEt3p9OdzvP749y2b3eqqmu51bV936/XfdWte8899VT1JTx16iwtbR1cdPszvPDGpiJFJiIi1SDRpNfMZgN/Ba4HzgYOBHYHBkSPB0bH5wFvmNm/Jfn6IlKZ3J1L7lqac8LbqaWtg0vv+hPu3nthERGpSYn06TWzfsCdwJlAplEl8XO7AnOjgVFnuf5vVXKap1dK5Ym/bMi6S0M6LzVu4cm/vM2UA0YlFJWIiFSTpAayfRv4MNCZuL4F/Ax4HFgFNANDgfHAMcBMYAwhCT4zuv7ShGKR/GmeXimJO558LbF6lPSKiEgqBSe9ZjYR+Cwh4TXge8BX3b0lRfGngLvM7CvA1cDs6JrPmtkN7r6i0HikIJqnV/rcltZ2HnyxMZG6HnhxLVta2zWPr4iI7CSJPr3nA3XR/jXuPjtNwvt37r7V3S8B/js6VAd8IoFYpADuPtfd35FpI7TMiyRm7abWbvPwFqJjh9O4uTWRukREpLokkfS+P3p8G/hajtdeBmzoUY+I1JDmPAevpdO0Ldn6RESkOiSR9O5N6NrwsLu353JhVP5hQheHvROIRUQqzNCBdb0XysGwQcnWJyIi1SGJpHd49Pi3PK/vvG54xlIiUpXG7jqYun7JLCXcv59RP2JwInWJiEh1SSLpfTt63CfP6ztbeN/OWEpEqtLwwQM48ZD6ROo68ZCxGsQmhWvdDOtWwF+fCY+tm0sdkYgkIIkpy14GxgLvM7PR7p71eqBmNho4jtA94uUEYhGRCnTu0fty//NrE6lHJC/usOoxeOoGWPEb8FjfcKuDiadAw4UwfhpYMr9MiEjfSqKl94HocRBwk5lllUibWR1wY3QdwG8TiEVEKtCU/UfxzvphBdVxUP1wjt5/94QikpqyZin8cArc+iFY/svuCS+E58vuC+d/OCWUF5GKk0TSewPQuej9ycBCMzss0wVmdijwEHBKdGhzVI+I1CAzY+5ZkxmS56C2IQPr+M5Z78bUAie5+vNCuPlkWL88u/Lrl4fyf15Y3LhEJHEFd29w9w1mdimh1daB9wJ/NLOlwGLgNbpWZNuHsCLbe6LLLbrmEndXn94S0zLEUkqH7rUr8847gotuf4aWHKYxGzKwjnnnHcGhe+1axOikKq1ZCvPPhfbm3K5rbw7XnX8/jJtcjMhEpAgSWYbY3W8ys92B/6JroYrJ0ZZKZ3NMB/Bld785iTikYFqGWEpq2oGjueuiKVxy11JebmzqtfxB9cP5zlnvVsIruXOHez6Ve8Lbqb0Z7r0YLl6sPr4iFSKRpBfA3b9tZosJie80uhLbdP4P+Iq7P5FUDFIwLUMsJXfoXrvy4OxjefIvb3P7k6t48MXGbiu29e9nnHjIWM49el+O3n93dWmQ/Kx6LPsuDemsWwarFsF+05KJSUSKKrGkF8DdFxNmcTiYMCvDewgJ0jCgCVgP/BFY6O4vJfnaUjh3nwvMzVTGzKYCj3U+b9m2vdhhSQ0yM6YcMIopB4xiS2s7jZtbadrWwbBBddSPGKxpyaRwS25Mrh4lvSIVIdGkt5O7rwBWFKNuKS/n37KE09YO59yj92XK/qPU6iaJGz54gJJcSVbrZlj+62TqWv6rUN/gEcnUJyJFU5SkV2pHhzv3P7+W+59fyzvrhzH3rMnqXyki5W3zmp2nJcuXd8CWN5X0ilSAgqcsM7Obou1jeV5/VnT9jwuNRUrr5cYmzpr3BI+9kvX6JCIifa8tz8Fr6WzrfdCliJReEvP0zgI+Dhyd5/UNUR2zEohFSqylrYOLbn+GF97Y1HthEZFSGDg02foGFbawioj0jSSSXpFuWto6uPSuP+HuvRcWEelrI8aFpYWT0K8/DN8zmbpEpKjKIentHPmkDKmKvNS4hSf/ovVGRKQMDR4BE0/pvVw2Dj5F/XlFKkQ5JL17RI/qFFVl7njytVKHICKSWsOF5VWPiBRdSZNeMxsBnEBo5V1dylgkeQ+8uJYtre2lDkNEZGfjp8HoiYXVMWYSjJ+aTDwiUnQ5TVlmZpdlOH1UL+fjBhCWuz0BGEtIehflEouUv44dTuPmVs2xKiLlxwzOuB5uPjm/pYgHDIXTf6QliEUqSK7z9F5B6r63RpiFoSHPONqA6/K8VspY07aE5sKscFta21m7qZXmtg6GDqxj7K5aVUyk5MZNhpl3wPxzc0t8BwwN142bXKzIRKQI8lmcIt3X2ny/7r4OfNrdX8zzeiljwwYlNEK6Ark7T/xlA7c/8Rq/W9ZIx46u74t1/YwTD6nXSnYipXbADDj/frjnU7B+ee/lx0wKLbxKeEUqTq5J75Upjl1OaP1dAvw2izocaAU2AC8AT7nmtioLZnYJcEkvxQZmW1//fkb9iMGFBVWhXnhjE5fctZSXG1OPz+zYoZXsRMrGuMnw6Sdg1SJYckNYoji+Ylu//mGWhoYLQx9efUkVqUg5Jb3uvlPSa2aXR7tPpTovFWUEoa91Ik48ZGxN/oT/2Cvruej2Z2hpy65rR+dKdvPOO4JpB44ucnRSS5rammhsaaSlvYUhA4ZQP6SeYQO1kEJKZrDftLC1bg5LC29rCgtPDN9T05KJVIF8ujf09H+E1ttXEqhLSmsz8EYvZQYCWWVm5x69b8EBVZoX3tiUU8LbqXMlu7sumqIWXymIu7Nk7RLmvzSfhasX0hFrsayzOmbsM4OZB82kYWyDutWkM3iEklyRKlRw0uvu0xOIQ8qAu88F5mYqY2ZTgcd6q+ug+uEcvf/uSYVWEdydS+5amnPC26lzJbsHZk9TMiJ5WbZhGXMWzWHlxpUpz3d4BwteW8CC1xYwYeQErp56NZNGTerjKEVESqMcFqeQKjNkYB3fOevdNZe4PfGXDWn78GZLK9lJvhavWcysB2alTXh7WrlxJbMemMXiNYuLHJmISHlQ0iuJGjKwjnnnHVGTP9EntQKdVrKTXC3bsIzZD89m6/atOV23dftWZj88m2UblhUpMhGR8qGkVxJzUP1w7rpoSk0OxtrS2s6DLzYmUpdWspNcuDtzFs3JOeHttHX7VuYsmoMm0RGRapfEQDapYXVm/OO79uTco/fl6P13r7kuDZ3WbmrtNg9vIbSSneRiydolWXdpSGflxpU83fg0DWPzXV9IRKT8KemVgtw8q4ETjj+81GGUXHOeg9fS0Up2kq35L81Ppp4V85X0ikhVU/cGKciQQfreBDB0YLIrz9XySnaSvaa2JhauXphIXQ+tfoimtsIGYoqIlDMlvSIJGLvrYOr6JdO1o5ZXspPcNLY0dpuHtxAd3sG6lnWJ1CUiUo6U9IokYPjgAZx4SH0iddXqSnaSu5b2lkTra25vTrQ+EZFyoqS3wpjZv5jZnWa23MzeNrN2M2s0s1+b2Smljq+WJbUCXS2uZCf5GTJgSKL1DR0wNNH6RETKiZLeynMJcCawFXgcuBtYDfwj8Csz+3YJY6tpU/YfxTvrhxVURy2uZCf5q68bQl1CU431d2dM3S6J1CUiUo6U9Fae84Hd3P1wd/+Qu5/t7g3ANKAZuNTM/qG0IdYmM2PuWZMZkuegtlpdyU7yN2zrJma05Dc/b08zmlsY1ro5kbpERMqRkt4K4+5PuPtOQ6zdfRFwZ/T0+L6NSjoduteuzDvviJwT31peyU4K0NbMzM1bEqlq5pYm2KbZG0SkeiWe9JpZg5lda2ZPR31N28xse4pyI83s5Gg7LOk4alTnMl7bShpFjZt24GjuumhK1l0danklOynQwKE0tG5jQltbQdVMaGvjyNZtMKiw7jkiIuUssUlWzWw0cDPwwfjh6DFVp7Nm4AZgLPAqMCHBWOqAicCRwBHR47uBzg5rt7r7rDzqPRU4D2ggxL0ZWAncA8xz95L9Nmhmk4GzgQ7ggVLFIcGhe+3Kg7OP5cm/vM3tT67iwRcbu63Y1r+fceIhY2t+JTsp0IhxmNVx9foNzNqznq39cm/H2GXHDq5evwHr1x+G71mEIEVEykMiSa+ZjQMWA3vTlehm5O7tZnY9cCWwn5kd7e5PJhEPcBdhsFcizGwY8BPg1B6nRkfbFOCzZnZWgu+ht5jOB94HDATGA0cTWnovdvcX+yIGyczMmHLAKKYcMIotre00bm6laVsHwwbVUT9isKYlk8INHgETT2HSsvu4tvEtZtfvkVPiu8uOHVzb+BaT2tph0umhPhGRKpVU94ZfAPsQEt5lwDlAPfDDXq6Lr595YkKxAPTsUPk28Eo+FUWtxj+nK+FtBL4OfAz4DGEGBQgJ//1mNjGf18nDe4GPEz7rKYTZHD4D3NRHry85GD54ABPGDGfy3iOZMGa4El5JTsOFABzT2sotbzZm3dVhQlsbt7zZyDGtrd3qERGpVgW39JrZ6YRWRgcWAR9095boXMa5dNz9FTN7AxgHJDnjwFPAcuAZ4Bl3f9XMZhG6X+TqQuCkaH8ZMMPdG2Pnr4umCbsU2A2YBxybqiIzuw04KsfXv8fdv9LzoLtfCFxoZkOBA4HPAv8DnGlmZ7p7MkO6RaS8jZ8GoyfC+uVMamvn7jfW8vTgQfxsxHAWDtmFjljXmf7uzGhuYeaWJo5s3db1s9yYSTB+aknCFxHpK0l0bzgretwOfLwz4c3Bc8BewEEJxAKAu38jiXqiVt7LY4fO65HwdvoSYcaEycA0MzvB3X+Xotw+5P4+M3ayc/dmYClwQfQl4wLCXL5X5/g6IlKJzOCM6+Hmk6G9GQMaWrfR0LqNJjPW9a+j2fox1HcwZnsHw3rO6ztgKJz+o1CPiEgVS6J7Q2cr7xPuviqP6zsXey/HoevH0pV0Puruz6Yq5O4dwPdih85JU266u1uO26wc4r0tejwth2tEpNKNmwwz7wgJbMwwd/Zv38672trYv3176oR35h3hehGRKpdE0jsmenw5z+ujDmUMSiCWpMVnori/l7K/TXNdX1ofPZbjFwgRKaYDZsD594euDtkYMymUP2BGceMSESkTSXRv2BE95ptAd665urHwUBL3rtj+kkwF3X2tmb1OGNBWb2aj3X19pmuKYHr0uDJTITPbh9DVIh/v6r2IiJTEuMnw6Sdg1SJYcgMs/zV4R9f5fv3h4FPCoLXxU9WlQURqShJJ73pgKGHarHwcHj2uSSCWpMX7376aRflXCUlv57WJJr1mdgQhWf2Vu2/vce4Uuvrx3thLVZ+ge1/lvLVs2kBTk1ZxkvLX0pLrcIMKNvo9cPIP4fhroHk9tLXAwCEwdHTXAhTNzaWNUbJSU/etVI1yvW+TSHqfBvYDjjazEbks0GBmRwEHEPoEP95L8VIYGdt/K4vyG9Jcm5S9gbuBjWb2DGH6tJGEBPuAqMx33P3O1JcXwS8/C83/CIedBXsdqZYjkXIyaJhWWRMRiSSR9P4K+ChhtbOvAl/O5iIzGwD8v9ihexOIJWnx/1u0pi3VJT5N2PCEY4EwFduVhEUpDgKmErqXrAFuB/7H3RcV4XXTGtK+gWHLfgrLfhr6Ep5xvQbFSFkbNkxJoFQe3bdSSYYMGVLqEFJKIumdD1xBaO39gpmtcffvZbogWrL4dsLcvE6YS/f3CcRS1dx9DeGzLtRNQL6f97tIt+jI+uVh2qSZd2hwjIiIiJSVgpNed99uZhcCD0b1fdfMziMkw/t3ljOzUwnTf72XsETwLtGpFsLcsuWoibDgBMDg6Hkmu8T2txQlogS4+2pgdc/jZnYJYY7fTAZmPNveDPPPDaPC1eIrIiIiZSKJll7c/REzOxe4hZD4HU7XALXOiSHviV3S2fGzCTjH3Z9PIo4i2EhX0rsHvSe9o3pcW2lGEBYKKUx7M9x7MVy8WH18RUREpCwkMU8vAO7+c8ISuw8Tktr4RornjwBT3P03ScVQBC/F9vfLony8zEtpS5WvzcAbvWzZzUixblmYNklERESkDCTS0tvJ3V8Ejjezw4CTgSnAOGBXoJkw28AfgF+7e8Z5b8vE88BJ0X4DIaFPyczq6ZqubF0J5ugtmLvPBeZmKmNmU4HHsqpwyY2w37QEIhMREREpTKJJbyd3fw54rhh197EHgC9G+x8ErslQ9uTYfm+rt9WG5b+C1s0weESpIxEREZEaV5Skt4o8CqwFxgLTzexwd3+2ZyEzqwM+Fzs0v4/iS1QiA9nivAO2vKmkV0REREpOSW8G7t5hZlfRNUXXbWY2w93X9Sj6TWBytP+4uz/YVzEmLJmBbHHbtFqbiIiIlF5VJr1mth87T4N2WGz/PWb29R7nF7r7whTV3QCcAXwAOAT4k5ndACwDdgfOISwSAWHGhosKi76kOgeyZTIQGJ11jVoNSkRERMpA4kmvmfUDDibMZDACGJDtte5+W0Jh7AvMyXD+MLonwQDbgZ2S3mge4g8DPwVOIXR1+FqKOv8KnB0N5qtIiQ9k69cfhu+ZQGQiIiIihUks6TWzfQnJ4EfIbwleB5JKehPl7luAD5nZacA/E2ZyGENYgOLPwN3APHffVLooy9DBp6g/r4iIiJSFRJLeKBn8CWFhipKvRuDujxQjDne/D7gv6XrLReID2RouLCgeERERkaQUnPSa2cHAnXRPhl4HXgDeJnQbkMqQ3EC2MZNg/NTey4mIiIj0gSRaer9MSHgdWAF80t0XJ1Cv9L1kBrINGAqn/0hLEIuIiEjZSCLpPS56bAE+4O5rEqhTSiCRgWwDhsLMO2Dc5GSDExERESlAEknvGEIr7++V8Na4MZNCC68SXhERESkzSSS96wn9QBsTqEtKKK+BbFYHk04Pg9bGT1WXBhERESlLSSS9LxOSXk3IWvlyH8j2sZ/D8ScWJxoRERGRhCSR9N4KzACONbOh7t6cQJ1SGrkPZBs0tJjxiIiIiCSiXwJ1zCfM2jAC+O8E6pMScfe57v6OTBtwZqnjFBEREclVwUmvu7cDpxJaCC82sx+b2R4FRyYiIiIikpBEVmRz95VmdgQwDzgfOMfMHgNeBDYRZnfIpp6rkohHRKpI62bYvAbammHgUBgxTstbi4hIzhJJeiMHArtG+4OB90dbLpT0igi4w6rH4KkbYMVvwDu6zlkdTDwlmjFkmmYMERGRrCSS9JrZLOAGQncJByzacpFVa7AUT15Tlokkbc1SuOdTsH556vPeAcvuC9voiXDG9ZobWkREelVw0mtmRwM/pivJ7QAWA88DfwO2F/oa0mdyn7JMJEl/Xgjzz4X2LCeBWb8cbj45rAJ4wIzixiYiIhUtiZbeLxASXgcWAee5++oE6pW+l/uUZSJJWbM0t4S3U3tzuO78+9XiKyIiaSUxZdnU6HEjcKoS3sqlKcukZNxDl4ZcE95O7c1w78WhHhERkRSSSHpHElp5H3L3TQnUJyK1ZtVj6fvwZmvdMli1KJl4RESk6iSR9K6LHt9OoC4RqUVLbiyvekREpOokkfS+ED3uk0BdIlJrWjfD8l8nU9fyX4X6REREekgi6b2DMJDtfVqJTURytnlN93l4C+EdsOXNZOoSEZGqksTsDT8DLgCOA+aZ2UfcNZqkEmmeXimJtjwHr6WzrSnZ+kREpCoUnPS6u5vZWcDdwOnAAjP7vLs/X2jd0uc0T6/0vYFDk61v0LBk6xMRkaqQxOIUN0W7q4E2QovvUjP7M6G/7yayW23N3f2CQuORgmieXul7I8aFpYWT6OLQrz8M37PwekREpOok0b1hFjsntQYcEG25UNJbQu4+F5ibqYyZTQUe65uIpCYMHgETTwnLChfq4FNCfSIiIj0kMZANQpIb31Id620TkVrVcGF51SMiIlUniZbe8xOoQ0Rq2fhpMHpiYQtUjJkE46f2Xk5ERGpSEgPZbk0iEBGpYWZwxvVw88n5LUU8YCic/qNQj4iISApJdW8QESnMuMkw846QwOZiwNBw3bjJxYhKRESqhJJeESkfB8yA8+8PXR2yMWZSKH/AjOLGJSIiFS+JPr0i0pvWzWHlsbbmMC/tiHGaZSCdcZPh00/AqkWw5IawRHF8OrN+/cMsDQ0Xhj686tIgIiJZUNIrUizusOoxeOoGWPGb7omb1YVpuhouDIO4lLh1Zwb7TQtb6+awtPC2prDwxPA99YVBRERyllXSa2aXxZ+7+1XpzhUiXq/0PS1DnKA1S+GeT6WfjcA7wry0y+4LP+Wfcb36pKYzeISSXBERKVi2Lb1X0H0BiqsynCuEkt7S0jLESfjzQph/bvazEKxfHmYtmHmH+qam0NTWRGNLIy3tLQwZMIT6IfUMG6ilhkVEJDe5dG/o/P01VYKbxG+zSSXOkj8tQ1yoNUtzS3g7tTeH686/Xy2+gLuzZO0S5r80n4WrF9IR6xpSZ3XM2GcGMw+aScPYBkxdQ0REJAvZJr1X5nlOKoiWIS6Qe+jSkM88sxCuu/diuHhxTffxXbZhGXMWzWHlxpUpz3d4BwteW8CC1xYwYeQErp56NZNGTerjKEVEpNJklfS6+5WxvrtP9TyXeFQilWjVY4WtKAawblmYtWC/acnEVGEWr1nM7Idns3X71qzKr9y4klkPzOLa467lmHHHFDk6ERGpZLnM03sFcDlwcnFCEalwS24sr3oqzLINy3JKeDtt3b6V2Q/PZtmGZUWKTEREqoEWpxBJQuvmMJ9sEpb/KtRXQ9ydOYvm5Jzwdtq6fStzFs3BXUMDREQkNSW9IknYvKb7PLyF8I4wL20NWbJ2Sdo+vNlauXElTzc+nVBEIiJSbZT0iiShLc/Ba+lsa0q2vjI3/6X5ydSzIpl6RESk+ijpFUnCwKHJ1jeoduahbWprYuHqhYnU9dDqh2hqq60vDCIikh0lvSJJGDEuLC2chH79w1K7NaKxpbHbPLyF6PAO1rWsS6QuERGpLkp6RZIweARMPCWZug4+paaW3W1pb0m0vuZ850kWEZGqpqRXJCkNF5ZXPRViyIAhidY3dEDCXU1ERKQq5LIMcaczzOzQxCMBd/fji1CvSN8YPw1GTyxsgYoxk2D81ORiqgD1Q+qps7pEujj0t/6MGTImgahERKTa5JP0jou2JBmgCTalspnBGdfDzSfntxTxgKFw+o9qbgniYQOHMWOfGSx4bUHBdc3YZwbDBtbOIEAREclePt0brAiblAEzu8TM/pppA+4udZxlbdxkmHlHSGBzMWBouG7c5GJEVfZmHjQzmXoOTqYeERGpPvm09D4F/DbpQKQsjAD2KnUQFe+AGXD+/XDPp7Lr6jBmUmjhrdGEF6BhbAMTRk4oaIGKCSMncGT9kQlGJSIi1SSvpNfdr0w8EikHm4E3eikzEBjdB7FUtnGT4dNPwKpFsOSGsERxvM9qv/5hloaGC0Mf3hrr0tCTmXH11KuZ9cCsvJYi3qX/Llw99Wqsxj9HERFJL5+kV6qUu88F5mYqY2ZTgcf6JqIKZwb7TQtb6+awtPC2prDwxPA9a2pasmxMGjWJa4+7ltkPz84p8d2l/y5ce9y1TBo1qYjRiYhIpdOUZSJ9YfAIGH0QvOOI8KiEN6Vjxh3DLSfdwoSRE7IqP2HkBG456RaOGXdMkSMTEZFKp5ZeESkrk0ZN4u5T7+bpxqf52YqfsXD1wm7TmfW3/szYZwYzD57JkfVHqkuDiIhkRUmviJQdM6NhbAMNYxtoamtiXcs6mtubGTpgKGOGjNG0ZCIikjMlvSJS1oYNHKYkV0RECqY+vSIiIiJS9XJNetV5TkREREQqTi7dG/aLHjcXIxARERERkWLJOul199eKGYhINWtqa6KxpZGW9haGDBhC/ZB69VMVERHpQxrIVuEszNf0EHBcdGiiu68oYUgScXeWrF3C/Jfm7zTtVp3VhWm3DppJw9gGTbslIiJSZEp6K9+/EBJepwR9rlvaW/r6JSvCsg3LmLNoDis3rkx5vsM7WPDaAha8toAJIydw9dSrtaKYiIhIEWn2hgpmZu8ArgEeAFaXIoZPP/RpLnnkEp568yncvRQhlJ3FaxYz64FZaRPenlZuXMmsB2axeM3iIkcmIiJSu5T0VrZ5QB3wqVIFsMN3sOC1BVzwuws485dnsmzDslKFUhaWbVjG7Idns3X71pyu27p9K7Mfnl3zn5+IiEixKOmtUGZ2HnAy8B/lMsiw1lss3Z05i+bknPB22rp9K3MWzVGLuYiISBFUZdJrZnVmdqiZzTKz75vZE2bWYmYebbfkWe+pZvZzM1tlZq1mts7MFpvZF81sRMJvI1Mc9cC1wBLge331utmo5RbLJWuXZN2lIZ2VG1fydOPTCUUkIiIinaoy6QXuAp4HbgY+AxwN7JJvZWY2zMzuA+4DPgLsCwwCRgNTCP1qXzCzowuMO1vXASOAT7r7jj56zazVaovl/JfmJ1PPimTqERERkS7VmvTW9Xj+NvBKPhWZWR3wc+DU6FAj8HXgY4SE+vHo+N7A/WY2MZ/XySGeDwMfBr7t7n8q5msVotZaLJvamli4emEidT20+iGa2poSqUtERESCap2y7ClgOfAM8Iy7v2pmswgtv7m6EDgp2l8GzHD3xtj568zs28ClwG6EwWXHpqrIzG4Djsrx9e9x969E1+9OaOVdCVyZYz19bv6K+TSMbSh1GH2isaWx2zy8hejwDta1rNPiFSIiIgmqyqTX3b+RRD1RK+/lsUPn9Uh4O30JOB6YDEwzsxPc/Xcpyu0DHJRjGHvG9r8L1AMfc/fWHOvpc50tlrWQvCU9X3Fze3Oi9YmIiNS6qkx6E3QsXUnno+7+bKpC7t5hZt8DbooOnQPslPS6+/QC4zkNaAUuM7PLepwbGz3eZmYtwA/c/RcFvl5BaqnFcsiAIYnWN3TA0ETrExERqXVKejP7YGz//l7K/jbNdUkbDLwvw/nO/gT3FjGGrNVKi2X9kHrqrC6RLg79rT9jhoxJICoRERHppKQ3s3fF9pdkKujua83sdcKAtnozG+3u65MMxt1HpjtnZqsIs0pMdPcVmeoxs30IXS3y0a2TbuvrmXtZLB+znE3DN+X5UpXloE0HsWRtxtskK0eNPYqlTy0tPCDppqUldEEZMiTZVnmRYtJ9K5Xoueee63moLH6+VNKbWbz/7atZlH+VkPR2Xpto0pugT9C9r3Le3rz9zYznz+bsJF6mprzKq9zJnaUOQ0REJCn7lzoAqN4py5IyMrb/VhblN6S5VkRERKRW7VbqAEBJb2/iI7CymS0hvv7s8IRjEREREalEu5c6AFD3hqrh7uNzKH4T8Ps8X2oaEJ8S7lLCvMgi5exdwA9jzz9NWLVRpJzpvpVKdRTwndjzx9MV7EtKejNroqtJfnD0PJP4UsdbihJRAtx9NbA6n2vNrOehp9x9UcFBiRRRivv2ed23Uu5030qlSnHvlsUYJ3VvyGxjbH+PLMqPSnOtiIiIiJSQkt7MXort75dF+XiZl9KWEhEREZE+paQ3s3jfqYa0pQAzq6drurJ1Sc/RKyIiIiL5U9Kb2QOx/d5WWTs5tt/b6m0iIiIi0oeU9Gb2KLA22p9uZoenKmRmdcDnYofmFzswEREREcmekt4M3L0DuCp26DYzG5Oi6DeBydH+4+7+YLFjExEREZHsVeWUZWa2H3BBj8OHxfbfY2Zf73F+obsvTFHdDcAZwAeAQ4A/mdkNwDLCZMvnAFOjshuBiwqLXkRERESSVpVJL7AvMCfD+cPongQDbAd2SnrdfbuZfRj4KXAKMBb4Woo6/wqc7e4v5hWxiIiIiBSNujdkwd23uPuHgNOBu4HXgW3AW8AfgC8Bh7r74pIFKSIiIiJpVWVLr7s/Auy0HEgC9d4H3Jd0vSIiIiJSXGrpFREREZGqp6RXRERERKqekl4RERERqXpV2adXimo1cGWP5yLlTvetVCLdt1KpyvLeNXcvdQwiIiIiIkWl7g0iIiIiUvWU9IqIiIhI1VPSKyVhZgPN7Ctm9qKZbTWz9WZ2t5kdXurYRFIxswFmdoKZfc/MnjOzZjNrNbOVZnadme1b6hhFemPBQjPzaDu41DGJZGJmu5rZVdG/u01mtsXMVpjZjWa2V051qU+v9DUzGwg8CEwH1gGPAnsCU4F24EPu/mDJAhRJwczeDyyInr4OPBPtHwWMAzYDJ7n7EyUITyQrZnYRcD3ghEWcJrr7itJGJZKamU0i/Ls7DlgJLAUGAhOAScA0d1+UbX2avUFK4UuEhHcJ8H533wxgZucAPwXuMLP93X1L6UIU2ckO4BfA3Hhia2aDgR8Bs4CfmdmB7t5emhBF0jOzdwDXAA8AEwH9OiFly8xGAr8DRgH/5O4/7XH+AEJjQ/Z1qqVX+pKZ9Qcagd2BBnd/usf53wAnA7Pd/f+VIESRnJnZLsCbwK7AdHd/tMQhiewk+vf1fcAhhF/Y9kUtvVKmzOy7wGzg8+5+bRJ1qk+v9LX3EhLeV3smvJE7o8fT+i4kkcK4+1bg5ejpuFLGIpKKmZ1HaFD4D3d/rdTxiGQS/YJ2PtAC/E9S9SrprWJmVmdmh5rZLDP7vpk9YWYtsQEMt+RZ76lm9nMzWxUN5FlnZovN7ItmNqKXyydHj8+kOf9s9PjufGKT6lCm927GeIHx0dO1+dYjlauc71kzqweuJXQp+14+cUj1KtN79wjCL2fPunuLmR1vZteY2fVm9mUzOyifmNSnt7rdBZyZVGVmNgz4CXBqj1Ojo20K8FkzO8vdn0xTTWcfsr+mOd95fHczG+buTYXELBWrHO/dTM6L6lkPLC4kVqlY5XzPXgeMAD7p7juSilGqRjneu4dEj+vM7OfAR3qcv9rMrnb3y3KJTS291a2ux/O3gVfyqShqyfo5XTdxI/B14GPAZ4DHo+N7A/eb2cQ0VQ2LHpvTnI8nucPziVWqQjneu+nqHw98J3o6x9235ROnVLyyvGfN7MPAh4Fvu/uf8olHql453ru7R48fInR3/HfgHcBY4HOEmZ6+ZmYX5BKfWnqr21PAckJXgmfc/VUzmwXcnEddFwInRfvLgBnu3hg7f52ZfRu4FNgNmAccm2/gUvMq4t6NfqK7j/AP9M/d/YY84pPqUHb3rJntTmjlXQlcmUccUhvK7t6lq1F2AHCFu38rdu77ZjaA0NhwGfDjrKNzd201tBGmVfJouyXLa+qANbHrDs9Q7o+xciekKDM3OvfdNHWMjF0/rNSfl7by2Up976a4ZjDwSFT+98CgUn9G2sprK/U9C9wanZuR4tyq6NzBpf6ctJXfVgb37udi58elOD8sdn7/bN+XujdINo4lLB4B8Ki7P5uqkLt30H2QxDkpinWOGn5HmtfqPP62qz+vFC7Je/fvolaG/yVM//QkcJqrW4MkI8l79jSgFbjMzB6Jb4SfiQFui4717DMpkqsk791V0eM2d1+Too4mwhgK6LqXe6XuDZKND8b27++l7G/TXNdpafR4RJrrO5chVt8zSUKS9y4AZtYPuIMw/dOfgJPdPV0fdZFcJX3PDiZ8OUunIXq8t5fXEulNkvfuH6PHQakGtUd9h0dGT7NuIFNLr2TjXbH9JZkKuvtawhKtAPVmNrpHkccJneT3M7MjU1RxdvR4Xz6BivSQ5L2LmRlwA3AW8BLhZ7m/JRSrCCR4z7r7SHe3VBtdv7pNjI5dm9QbkJqV5L37Ol1TmB6XoopphP6+LUDWi6so6ZVsxOfDezWL8vEy3ebSc/fthPkiAX4Yn6vPwjLEJwNvATflFalId4ndu5HvAJ+Iyh3v7usKiE0klaTvWZG+kvS9+1/R47eiWXIAMLNxQOeKrTe6e1u2Aap7g2RjZGz/rSzKb0hzbaf/BmYA04FXzOxRQp+caYRpSM5z9y35BCrSw8jYfkH3rpmdBnw+erqKME9kqjrudfd7sw1QpIeRsf0k/r0V6SsjY/sF37vu/gsz+xFwMfC8mS0GOoBjCAtXPAl8JZcAlfRKNobF9luzKL81tr/TXLvu3mZmJwJfAM4lzOfXROjScFW6zu8ieUjy3t0ttp/q57ZOq1D/SMlfov/eivShxO9dd/+0mS0C/pWwqMUAwpLvPwOudfdsXufvlPRKSUQ/R3wj2kTKnrvfAtxS4jBEEuHu40sdg0g23P2nwE+TqEt9eiUb8ZGRg7Mov0tsX90UpJR070ql0T0rlars710lvZKNjbH9PbIoPyrNtSJ9bWNsX/euVIKNsX3ds1JJNsb2y/LeVdIr2Xgptr9fFuXjZV5KW0qk+HTvSqXRPSuVquzvXSW9ko3nY/sNaUsBZlYP7B09Xefu6zOVFyky3btSaXTPSqUq+3tXSa9k44HYftqVqiInx/Z7W5FFpNh070ql0T0rlars710lvZKNR4G10f50Mzs8VaFoWcDPxQ7NL3ZgIr3QvSuVRvesVKqyv3eV9Eqv3L0DuCp26DYzG5Oi6DeBydH+4+7+YLFjE8lE965UGt2zUqkq4d41d++r15I+Zmb7ARf0OHwY8KFo/zngVz3OL3T3hSnq6k/4CeID0aG1wA3AMmB34BxganRuIzDV3V8s8C1IjdK9K5VG96xUqlq6d5X0VjEzmw48nONlV7r7FWnqG06YIPqUDNf/FTjb3Rfn+Loif6d7VyqN7lmpVLV076p7g2TN3be4+4eA04G7gdeBbYQ1tv8AfAk4VP8AS7nRvSuVRvesVKpyvnfV0isiIiIiVU8tvSIiIiJS9ZT0ioiIiEjVU9IrIiIiIlVPSa+IiIiIVD0lvSIiIiJS9ZT0ioiIiEjVU9IrIiIiIlVPSa+IiIiIVD0lvSIiIiJS9ZT0ioiIiEjVU9IrIiIiIlVPSa+ISAmY2XQz82h7pNTxVCsz+6iZ/crM3jCzbeXwmZvZI7E4ppcqDpFao6RXRArS43/gbmarzWxQltdeEbtufrFjldphwU+Au4BTgHHAwNJGJSKlpKRXRJK2N3BRqYOQmvexaOv0FHALcF203VOCmLLW4wvhFaWOR6Qa9C91ACJSlb5qZje6e0upA5GadV5s/3J3v6pkkYhIWVBLr4gUQz3wuVIHITXt8Nj+j0sWRQruPt3dLdoeKXU8IrVCSa+IJOnJ2P4XzWxEySKRWrdbbP/NkkUhImVDSa+IJOkO4KVof3fg0hLGIrXt79333H1HKQMRkfKgpFdEktQBXB57/nkzG1VIhWY2PjagZ1WW16yKXTM+2zJmNsHMvmVmL5jZJjPbamZ/MrOvmtmQFHUcZGY/MLPnzWyzmW00syfN7F/NrC6P92pmdqaZ/dLMXjOzVjNba2a/M7N/NrOc/s02s1FmdqmZLTCz16P6NprZMjO7zsyOzKKOnQZUmdkuZnZBFNdqM2uLzk/O9T2neL0TzewmM3s5+ky3Rp/FPWY2y8wGZLj273/THse951ZgjA3R3/1ZM/ubmW2P4nwz+vv/yMzOMrOhaa5PO2VZ5zm6/3d0ear3YGa3ZIjRzOwMM7s1+iw3RX//183sXjP7uJllNa7HzA42s2ui9/ZW9PduNbN1ZvaMmd0c1bdb77WJlJC7a9OmTVveG/AI4NH2KcCApbFj12S49opYuflpyoyPlVmVZUyrYteMz6YMcC7QHDvWc3sW2C12/X8Qkvx05R8GhmSIcXqs7CPAcODeDPU5sBgYk+Vn8K/Axl7q20Ho7zowy7/RFcBE4IU09U0u4D4aA/y+l3gdeBk4Mou/acYtzxj7A/OyfQ3g61n8NzM9w7netlvS1H8Y8Mcsrl8BTOrlPV8BbM8ynjv68t8ebdpy3TR7g4gkyt3dzL4G/DI69Bkz+667l3O/yg8CPyD8+vUKYXqrVkLy0BCVeQ8wHzjRzL4C/Gd0/DngT4TE4CjgkOj4dGAu4YtANm4BTiMkD08By4BBwDGEpBxgCvCQmb3X3Tenq8jMrgX+LXboLeAJYC0wOHovhxK+oHwCGGdm/+i9dwMYBTwA7EP4fBYBrwHDgKOze5sp460HHgcOiB3+M/AHYBswCfiH6PiBwMNmdpK7P96jqlujGCEk/Z2uyze2Hr4F/Evs+RuEv9V6wr0zKor1oAJe4x7Cl4qj6Lr3lkSv09OTPQ+Y2bHAr4DO/vTt0fWvRPvjgamE++AgYLGZTXH35Snq+je6tzi/Fb3mm4T7dHfgYMIXoZx/2RDpc6XOurVp01bZGz1aemPHn4wd/0Gaa6+IlSllS28rsBn4SIpyZ9O9pWt29PwN4H0pyl8SK9uR4fWnx8ptix7/QopWTOBCoC1Wfl6G9/6JWLlN0bUDUpQ7DvhrrOy/Z/E3ao8efw6M7lGuX6rXyfLvdX/sNZqAmSnKHElIhDvLrQZGZqizoFbdFPWNir3/7cDHAUtTdk/gs8AFWfw3Mz2Lz/2KLGMcCzTGrrsV2DNFuXrg7li554C6HmX6E5L5zjJfTvf3JSS/56e7h7RpK5dNfXpFpFj+I7b/STPbt2SR9G4gcKa7/6LnCXe/E7g5dui7hAT0/e7+aIrycwk/00NIBM/K8vWbgQ+4+9Mp6rwRuDh26JNmdkDPcmY2HPhO9LQNOMHdb3T39hR1Pgx8gJDwA/x7qn7LPfQHfgec7e7re9S3I9Xr9MbMjiO0tHc62913Wp0v+lyOJyTyEBZB6ctp8abQNThuvrvf6u6eqqC7v+nu33f3vp4q7WpCNxGA77n7xz3FLyzu3gh8FFgYHXoX8JEexQ4G9oj2H3f3b6b7+7r72+5+s7tfU/A7ECkiJb0iUhTu/ntCixaEpO6y0kXTq19G8abzsx7P53mKn4PTlD8qyxjmuvuf052MEqhnoqdGaMHt6RPAyGj/h+7+h0wvGL2HW6Ono4CTsohztic7G0J89b5fuvtv0hV091XAN2KHPmVmlmAsmcSn31uftlSJmNloQr90CN1YvpSpvLt3AHNih/6pR5Gyfr8i+VDSKyLFFG/t/biZHViySDLbqYW3h+dzLP9CbH+/LGO4Lccyx6U4f3Js/6dZvu7C2P7UXso+10uyn4/4+7gpi/I3EwbgQehGUEj/2Vy8Hts/08zGpC1ZGu8nfLkEuNvdWzMVjvyB8AsD7Py3j7/f48zsnQXGJ1JyGsgmIkXj7o+b2W8JP1/XAVcCHyttVCm90Mv5v/V4/mIv5d+O7WezQMdb7r4yi3JPxPYnm5n1+Il9Smz/X8zs41nU+Y7Y/t69lH2ml/M5MbO96Po5HsLsFBm5+3oze5nw8zuElddWJBlXGk8SEsG9CQP5XjSzmwmDxv7g7m19EEMm8b/9YWb2gxyv383Mhrp7M4C7v25mTxIGKO4KPGNmtxMG2j3uWmJcKpCSXhEptv8g/GxuwNlm9l/u3rPltNQ2ZTrp7tt7/IqesTxhoFOntPPKxqzOokzPcoMI05xtBjCzYdHzTqm6P/Smt3lWk/6Ze3Rsf2vPfsIZrKIr6d0jQ7nEuHu7mZ0H/JowW8UewBejrdXMngb+jzAob3G6/r5FNC62P5XeW+1T2Y2ull+ACwi/BNQT3vPF0bbdzJYS3u+DwENRdwmRsqbuDSJSVO7+LKF1CMK/Of+ZoXip5JSgFCGhybbVrLnH83iSu2sCcfTWELI1gdeIGxbb7/neMomXHZ62VMKigYvvJnQziX8WgwlJ5lcJ07itMLPT+yquSOJ/f3dfRni/36f7F73+hNk0LiEkva+ZWT5fskT6lJJeEekLl9HVD/M0M2vIVDgBlfZvW2+zJnTqucLXlth+z6Rxd3e3HLfpeb+D/DTF9lOuXpZGvOyWtKWKwN3/4u4fJ7RSnwR8nbAQSTwJfidwj5ld0oehxf/+l+Txt7dooGA37t7o7p8jtPZOB74G/JboF4bIXsANZva94r09kcJV2v8YRKQCufuLdB9Y9fUcLo9Pk5Rtl6wkWr36Um99aVOV20Ys4XP3jdGxTmMLD6vo4t0ZdjGzbLsqjI/tv5VcONlz92Z3f9Ddv+buMwizX3yU7oMe/yvqt9wXGmP7if/t3X2buz/q7l9395MJ3Ts+SGjZ7vTZPvhCK5I3Jb0i0leuoKuv6wnRylHZiLco7dbbFFVmtg/ZDR4rJ6NTzbubQnyw0tIU3Sziq3a9t/Cwisvd3wDWxQ4d09s1UWIcn0ng2aTjyoe7b43meZ5OVwI6EDgx3ypzLB+fnq7of3t3b3f3BwizRsQHgn6o2K8tki8lvSLSJ6I5aOOLPGTV2uvuW+iaDWEI3ROeVLJZDKIcnZdjmYdTnP91bP/iPpzDthDx9zEri/Kz6Pp/1xrgpYTjKYi7v01YUrlTfZ5Vxaccy2Yw5IN0fak8xszenefr5sTdtxEWLOmU7/sVKTolvSLSl/6Trp/gp5F9K1i8FWtWukJm9g7gK3lFVnqXmFnaOX3NbBbQ+dOxA6lW+5oHbIz2Dwcuz/bFzWwPM6vLtnyC5sX2zzCztPdEtKpffEGFeX01S4KZjcqheLwbyrq0pTLbENvvtYtE1Gp+R/TUgNvMLKtfPMysX7S4RfzYbmaWbY6QxPsVKTolvSLSZ9z9dbonOUdneWm8P/AlZvbhngXM7GjgUcK0S6WeMzVXbYRZCBaY2eE9T5rZ+XT/3H6cal5fd98EfD526HIzuzXq8rETC95rZj8kTIe2SyFvIh/Rcsi/jR36hZl9tGc5MzuCsLzzyOjQ60BfDpz6rJktNbOLzSxln1kzG2ZmV9P15aSD7q2guYh3GTjBzLLppz4H6Fx2+DDgKTM7IV1hM3uHmX2e0Fp+do/TpwEvm9kXzGx8musHmdln6L6E8W9TlRUpB5qnV0T62jcIc8hmO2MBwHzgC4TpkwYSEqNngaWERS8OA94Tlb0COB/YN5lw+8QThC4cZwBPR4sCLCfMxTsF2D9Wdjnhs0jJ3W8xs/0Jo+wB/hn4p2he1RWEGROGERalmEx5DPo7n9Al4ABCbHeZ2SuEFv42YBLwD4QWTAgzFZwTDd7rS+8GfghcZ2Z/JiSmbxG6H+xJ6JMcn4btm9EXvXw8RddiGHsSpkH7XfR6na3bS9z9zs4L3H2NmZ1GmCt4D8JqdQ+a2RtRfeujWPcADqX31QIPAL4FfMvMVgPP0dWSO5bwpXX3WPmfuHuvC4yIlIqSXhHpU+7eGE1t9OUcrtluZmcQWvo6E8DDo+3vxQgJ9VWEJKrSzCIkJKcQEt0pKcr8ATg9atFNy90vM7MXgO8SFi2oA46ItnSeovtMGX0muifeS2jRnxEdPjDaeloJfMzdl/RVfJH41GgGTIi2VNqAq939qnxfzN13mNmngf8lfNEbS/gCE3crcGeP65aY2ZGE7i/HR4f3InyhSqcReKXHsSbCf1OdXzT2ibZUdgDXA7MzvIZIySnpFZFSuIawslPWrYzu/qqZHQZ8FjiTMKBtEGEw02PAj9z9DwCVMX6rO3ffbGanEn4q/jih9bqe0Ef3OeAnwG3uviNtJd3ru8vM7gNmEvpONxDmlh1GaCl9g9Bq/Bhwv7u/nOgbypG7NwLHm9lJhJ/apxISvQGE1sU/AvcCd7h7nyfn7v4dM/tf4AOEFt13EaZOG0FI+jYSPs+FhL/Tawm85q+jBPZfCZ/HPoS/X8YbPHrt95vZFMI0ascSWox3Iwx220BIcp8mdL94xN2396jjF2a2J3ACYTaIdxO+cI6MimwCXiZMWXZbtJCFSFmzvl8pUURERESkb2kgm4iIiIhUPSW9IiIiIlL1lPSKiIiISNVT0isiIiIiVU9Jr4iIiIhUPSW9IiIiIlL1lPSKiIiISNVT0isiIiIiVU9Jr4iIiIhUPSW9IiIiIlL1lPSKiIiISNVT0isiIiIiVU9Jr4iIiIhUPSW9IiIiIlL1lPSKiIiISNVT0isiIiIiVe//AwyvG0I5PKAgAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 675x675 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(figsize=(6.75, 6.75), dpi=100)\n",
"ax.scatter(list_num_sites, times_mg, color='C0', label='pymatgen.core.Structure')\n",
"ax.scatter(list_num_sites, times_ase, color='C1', label='ase.Atoms')\n",
"ax.scatter(list_num_sites, times_ndarray, color='C2', label='numpy.ndarray')\n",
"\n",
"ax.set_xscale('log')\n",
"ax.set_yscale('log')\n",
"ax.set_xlim(1, 1e6)\n",
"ax.set_ylim(1e-4, 1e2)\n",
"ax.set_xlabel('Number of sites')\n",
"ax.set_ylabel('Time to create structures (s)')\n",
"\n",
"ax.grid()\n",
"ax.legend()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "39997493",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"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.4"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": true
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment