{ "cells": [ { "cell_type": "markdown", "id": "45511424-de0d-4620-88dc-1631789cfe30", "metadata": {}, "source": [ "# Benchmark" ] }, { "cell_type": "markdown", "id": "1edd243a-739b-4a67-b29d-8788eca0932d", "metadata": {}, "source": [ "### Benchmark for CutoffNeighborSearch in mdapy, freud and ovito" ] }, { "cell_type": "code", "execution_count": 1, "id": "d5b27ac4-c8fe-44ca-8e37-8b06148a17eb", "metadata": {}, "outputs": [], "source": [ "import mdapy as mp\n", "import freud\n", "import numpy as np\n", "from time import time\n", "from ovito.data import CutoffNeighborFinder\n", "import ovito" ] }, { "cell_type": "code", "execution_count": 2, "id": "f0d1f0ec-76e9-4a10-a8d1-85082031ce81", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "mdapy version is: 1.0.0\n", "freud version is: 3.4.0\n", "ovito version is: (3, 14, 0)\n" ] } ], "source": [ "print('mdapy version is:', mp.__version__)\n", "print('freud version is:', freud.__version__)\n", "print('ovito version is:', ovito.version)" ] }, { "cell_type": "code", "execution_count": 3, "id": "f6877409-471f-4b4b-88b9-8811508c3e2e", "metadata": {}, "outputs": [], "source": [ "def test_neighbor_average_time(ave_num=3):\n", " time_list = []\n", " cutoff = 5.0\n", " for num in [5, 10, 15, 20, 25, 45, 65, 85, 105, 125]:\n", " FCC = mp.build_crystal('Cu', 'fcc', 3.615, nx=num, ny=100, nz=100)\n", " ovi_data = FCC.to_ovito()\n", " print(f\"Testing system with {FCC.N} atoms...\")\n", " freud_t, mdapy_t, ovito_t = 0.0, 0.0, 0.0\n", " for turn in range(ave_num):\n", " start = time()\n", " FCC.build_neighbor(cutoff, max_neigh=50)\n", " end = time()\n", " mdapy_t += end - start\n", "\n", " start = time()\n", " finder = CutoffNeighborFinder(cutoff, ovi_data)\n", " finder.find_all(sort_by='index')\n", " end = time()\n", " ovito_t += end - start\n", "\n", " box = freud.box.Box(Lx=FCC.box.box[0, 0], Ly=FCC.box.box[1, 1], Lz=FCC.box.box[2, 2])\n", " pos = FCC.get_positions().to_numpy()\n", " shift_pos = (\n", " pos\n", " - np.min(pos, axis=0)\n", " - np.array([np.linalg.norm(FCC.box.box[i]) for i in range(3)]) / 2\n", " )\n", " start = time()\n", " aq = freud.locality.AABBQuery(box, shift_pos)\n", " nlist = aq.query(\n", " shift_pos, {\"r_max\": cutoff, \"exclude_ii\": True}\n", " ).toNeighborList()\n", " end = time()\n", " freud_t += end - start\n", "\n", " time_list.append([FCC.N, ovito_t / ave_num, freud_t / ave_num, mdapy_t / ave_num])\n", " time_list = np.array(time_list)\n", " return time_list" ] }, { "cell_type": "code", "execution_count": 4, "id": "3244fcd1-e870-4230-adc1-0bd6b0d770e6", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Testing system with 200000 atoms...\n", "Testing system with 400000 atoms...\n", "Testing system with 600000 atoms...\n", "Testing system with 800000 atoms...\n", "Testing system with 1000000 atoms...\n", "Testing system with 1800000 atoms...\n", "Testing system with 2600000 atoms...\n", "Testing system with 3400000 atoms...\n", "Testing system with 4200000 atoms...\n", "Testing system with 5000000 atoms...\n" ] } ], "source": [ "time_list = test_neighbor_average_time()" ] }, { "cell_type": "code", "execution_count": 5, "id": "1335ef0b-fc3f-47e7-afed-94fd4d8f3008", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAr0AAAJHCAYAAACdLe6nAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAewgAAHsIBbtB1PgAAwqpJREFUeJzs3QV4XGXaBuBnPJmkaeruLqnS0lKl7qXUKFRxWVhYZPdnd9EVFhZZdmGRpS0tLfVC3Y0W6u7uLmnT2Oh/vV+YNCcTmaRJxp77ukKSc05mvkxK8s43r+jcbrcbREREREQhTO/vBRARERERFTYGvUREREQU8hj0EhEREVHIY9BLRERERCGPQS8RERERhTwGvUREREQU8hj0EhEREVHIY9BLRERERCGPQS8RERERhTwGvUREREQU8hj0EhEREVHIY9BLRERERCGPQS8RERERhTwGvUREREQU8hj0EhEREVHIY9BLRERERCGPQS8RERERhTwGvUREREQU8hj0EhEREVHIY9BLRJQP58+fR40aNdCzZ0+43e5834587dKlS9G/f3907drVp685fvw4Xn31VZQqVQonT57M931T3l26dAnvvvsuKleujDVr1vh7OUSUBwx6iShozZgxA8WLF4dOp0t/+93vfpft9VevXkW1atVgNBrTr7darXj00UfzfN8bNmxQAeeyZctw7dq1fK3fbrfjt7/9LR5//HEsWLAATqcz16+ZMGECxowZg3/+85+4fv06AtmsWbPQqVMnlCtXDjVr1kRcXBxeeuklrFy5Ei+88AL27t2LYLJ48WI88cQTeOONN3Du3Dl/L4eI8spNRBTEHA6He+bMme4SJUrIdqt6mzx5co5fc+TIEXdUVJS7e/fu7pSUlHzdb3x8vLtXr17uF154wX23pk6dqtbdqVMnn6632+3uyMhI9TUnTpxwB6KXXnrJbbFY3N9++63b6XSqY7LWMWPGpP+c9uzZ4w5GdevWVetfvXq1v5dCRHnAnV4iCmoGgwFDhgzBzJkz0489+eST2L59e7ZfU7t2bTRq1AgDBw6ExWLJ1/3KDrPs/P3rX//C3SpTpkyerped6pIlSyJQLVq0CB9//DFef/11jB49Gnp92p+a6tWrY+LEifj973+PYJbXnxcRBQYGvUQUEmrVqpUeBCcnJ+OBBx7AlStXsr0+MjISUVFRCAQSxBbF1xSVKVOmqPctWrTI8vxf/vIXVKhQAcEqkB97Isoeg14iCinvv/++en/mzBkMHToUDofD30sKO55c42nTpmUbNMrPhoioKDHoJaKQIoVsnsK0tWvX4sUXX8zT1//yyy+qk0KDBg3UTvC9996L+fPne123bt06jBgxAvXq1cu2w8KoUaNQpUoVlC1bVhVxSRGaL37++We0a9dOFdnJ+4MHD+Z4/a5du1TnB1mvpG5IykV2HSWWLFmCXr16oUmTJupl+rZt2+Lbb7/1uu727dvqdiQlQboU7NmzB82bN1dFaatWrcpxPU2bNk3f8R03bpy6rcxatWqV5dcmJSXhnXfeQcuWLdV9yY7ws88+m2XR3urVq9G9e3d1fyVKlFA/izfffFMVCGYkj58UoNWvX199/umnn6rOFx07dkRiYmL6dUeOHFHpGHXr1kXFihXVz+yzzz7LscAwNTVVpWvIYylv0tmBiAJUXhKAiYgClRRJeX6lpaamujt06JBeMDV+/Hiv66VobMKECZpj8nmDBg3c+/btU5/v3bvXXa1aNbdOp9Nc+/LLL6cXM8n5zOTrpLCuffv27hs3bqhCrj/96U/qejkuXyv34ynkkoIoTyHbsmXL3FarVd2u0WhUxxs1auR2uVya+5Dzcm769OnumJgY9bkUjnm+59/+9rde63r77bfd5cuXd2/evFl9fvPmTfewYcPU9SNHjkwvOJPbrFOnTvptSaGgfO5Zj1ybk3PnzrlLliyZ/vVVqlRRt5EbeaxatGjhfvPNN1Wxns1mc//f//2fuo3GjRu7b9++nX7tlClT1M/lH//4h/o8MTHR3a9fP3Xt66+/nn7dE088kV70J4+RfF3x4sXT17ZixQp13fLly92xsbGq8E4ea7nvgQMHqmvGjRunWaf8nOS4/KykmFF+pmXKlEm/zRkzZuT6vRJR0WPQS0QhF/SKy5cvu6tXr66OSTC4cePGHIPeAwcOuM1ms9d1EsDIbRQrVkx1bPDYsGFDtkFvx44d1bn169drukzUqFEjPYjMyBP0yvmHHnrIfebMGXX84MGDqsuEnNu2bVuWQa90oJBuFOL69evu/v37pwdfskaPpUuXZvkEQJ4geALcDz/8MP24BJ0SrMpxeQJx/Phx9/bt290jRoxw//LLL+7cbNq0yV22bNn0tcjbfffd5/7555+z/ZrRo0erIDIjCUAl6Jevf+ONN9KPN2vWTB2TdXmsWbNGHZNzGXm+99KlS7ufeeYZ9bP4/PPP3Y899pg7OTnZffbsWRXwPvXUU5qv+/7779PXfvXqVa+gt127dqpTiDxZkDcJsOW4/AyIKPAw6CWikAx6xe7du93R0dHqeMWKFd0XLlzINuh99tln3ZUqVfK6XQk8s9rBO3z4cJZBr+yeyg6knJOPM/rd736njktwl1XQ26RJExWQZSRBoJybPXt2lkHv0aNHNcdlN9QTbEoQ5tGqVSu1LtlNzUwCQE9QmLGFmwR1cvyzzz5z54c88Rg1apQm8JU3OXbt2jXNtefPn3cbDAb3119/7XU78kRAvq5hw4bpx4YMGaJ2rSVo9ZDHwvPkISN5UiDH5QmErCkz2bmX81u2bNEcl5+f3PfYsWPTd8EzBr2Zd3TlCZMcj4uLy9PjRERFgzm9RBSyJCdT8kqlZZZMUBs8eDBsNluW10qequSNSt5nxjfJ75X8T3m7ePFi+vUmkynL25F80uzyaWUwhsh4OxlJXqp0n8hI7ldkzD3NKPP1ktcrLdvEjh071PujR49iy5Ytqs1abGys12307ds3fXjHpk2bvLoUNGzYEPkhOa6TJk1SedKSm+wxefJkldN79uxZTY605M5KPm/mn8FPP/2kHodbt26lXy8t6qRYMSIiAi6XS7VJe/nll9U5+Twjz/dRunTpLNuNydeKOnXqaI7HxMTg+++/V7nYnrZrmb+/vPysiMi/GPQSUUgbMGAA/vrXv6YXiD3//PNZXnf69GkVJEvRU8a3w4cPq2BQ3rL72syBjyd4OnXqlOacJxjOS7suT8CWOZDLiadgKz4+Xr3fv3+/ei8T6LJStWrV9PZt8jgUtDZt2mD9+vVqgp4U9nkK/WQSnYfnfr/88kuvn4EEx/L4S5CbmVwvo6APHTqEP//5z/lan2eUsy8T8Qr6Z0VERYdBLxGFvD/84Q8YOXKk+virr77CF1984XWNtDaTHdHsdmnz4qOPPlI7g9OnT9cc3717t3ovwzTyKi/rkm4RwrOrKx0RxI0bN1QP46x4rpXdzbslo4Yzd1AQ0qZMHgPpGCGWLl2aPs7X01pOnmT44sSJEyqYloBeRjjLfXp2WvPKE6QeO3YMBaEg/g0RUcFj0EtEYeHrr79W7cfECy+84BVcye6rpDcsX748y6+XYDHjS/856devH7777jvVGmvlypUqCJJWYZJqMWbMGHW+MElwK5o1a6YZ3CFkRzQrnqCzcePGd33/8lhlNxFPgmt5bDzpAp7dW8/ut+wGZ8fTKi0lJUXt7gppq5bfqXoelStXzvG+ZY3Z/bsgouDBoJeIwoLkfv7www8qwJFdyAsXLmjOS89WISkM165d8/r6f/zjHz6//C27l3/605/Ujqb0Da5Zs6bq3yoBmq+9eu/Gtm3b1PuHH344fTJapUqV1MeSo5qZBJEyve6ee+5Ray3IqWxZkfuQ/GVJt/DkOXfo0CE9BeXf//6319fI+saPH68+XrFiheqpK2kZWclrekHnzp3TnxhllW8tfXilbzMRBTcGvUQUEjzFQxLAZad8+fKYN2+eGvqQmez+yu6j7AC3bt0aU6dOVTt8+/btU0GP7NTKy+kenoK4rF7Gl1xVKXST29i6datKm5Cc1qeeeirLvFrPbeQUVGd1P1m9lJ6QkKCKx6Ro7/77708vdvPku8ogips3b2q+ZtmyZSpQ/Nvf/pZl8OhJj8gLybWVADYrkuIgTyykSNCzwyu70ZJ/7flZPP3002q3WH4GMhxEvpeHHnpIsy5Jj9i5c6f6WB7jV199NT2XWR4vz85tbt+H3J/k48rjIsMuZN3y9fJvQZ44yL8Lz26wLz+v7H5WRORnRdQlgoioUL3zzjs+DwaQPrmZB04I6VObub2WvEkP10OHDmmu/fLLL9U5vV6v2pplVKtWrSxvx2QyqbZozz33nDshISH9es8ABmnBdeXKFU2vXBnWIOfGjBmjuY/OnTunD7TwrE3agPXo0UP17pVhDZn73UrrLU9v30uXLqnjO3bscFetWtX91ltvaa6XvrSe1mfSYi1zK7WcSL9b+ToZCvH++++r/sEea9euddeuXVv13s3cPuzixYuaoRgZ36SlnId8nfRNluPS5kz6CUu7tXnz5qVfX6pUKffUqVPV9dJL13NchlBk5dNPP83yfmWISMYevdL2ToZRyLl3331Xcxtz5sxJ/zeRuZUcEfkfg14iCmpffPGFZvqXvFWoUMG9YMGCHL9OppNlDnrFwoUL1TAGmYomk7sefPBBr6BWAk5PL155k6EWv//97zUT2WRClwRwEiBJsJs5mHrggQfUtZ4BGp43GaQhQxgmTZqkgu2M58qVK6eGSQgJQmV6WJcuXdQ6JZCUCXASjOcUoMpwinvuuUd9jfS9lSB58eLFmms+/vhjzXQ3eZPrpe+xL2QanAyKkMlujz/+uJpAV7NmTRWYymQ1eYKScbpaRhK4y9dXrlxZPa716tVzf/LJJ14T6WTghASkElj37NnTfezYMXV86NCh6n6++uor9Xnr1q29Hvv7778/y/ueP3++u02bNu6IiAj1BOT555/XDCSR25Q1Zf6ZSFAvAzEk2M34c3zvvfd8eryIqGjo5D/+3m0mIgolb731FurVq4cRI0ZojsvL7PISuvSkfeSRR1T+aHR0tN/WSUQUTtKaChIRUYGQQQrTpk1T/WUzk9xQKeAaOHCgKozydEwgIqLCx0I2IqICJMVUUjiW04to0ls2MjIyy+loRERUOBj0EhEVIJk0JkGtDKCQzgIZg1/p+DB79mzVp/ef//ynX9dJRBRumNNLRFSAZArb6NGj01uaSc6uTEiTfN7z58+rdmmzZs1C165d/b1UIqKwwqCXiKiAydSzjz/+WE3xOnv2LMxmsxrIIH1pX3zxRZQuXdrfSyQiCjsMeomIiIgo5DGnl4iIiIhCHoNeIiIiIgp5DHqJiIiIKOQx6CUiIiKikMegl4iIiIhCHoNeIiIiIgp5DHqJiIiIKOQx6CUiIiKikGf09wLCWfny5ZGYmIiqVav6eylEREREAeP06dOIiorCxYsXC+w2udPrRxLw2u329M+PHTum3ij48GcXvPizC1782QUv/uyC17Ei+tlJfCRxUkHiTq8feXZ49+3bp943atRI8zkFD/7sghd/dsGLP7vgxZ9d8GpURD87z/0UJO70EhEREVHIY9BLRERERCGPQS8RERERhTwGvUREREQU8hj0EhEREVHIY9BLRERERCFP53a73f5eRLhiyxYiIiKioomRuNNLRERERCGPQS8RERERhbywDnqTk5Px6aefokqVKjh58mSO1+7duxcPPfQQypYtC4vFgpo1a+L555/HhQsXimy9RERERJQ/YTmGOCEhAf/973/x0Ucf4dKlS7lev2TJEgwaNAgpKSnpx06cOIH//Oc/mD59OpYuXYrmzZsX8qqJiIiIKL/Ccqd32bJl6NChAx588MFcr7127RpGjBiBNm3aYN68eTh8+DDWrFmDvn37qvNXrlzBwIEDkZSUVAQrJyIiIqL8CMugd/DgwWjbti0ee+yxXK/9+uuv8cADD2DVqlXo378/6tSpg06dOmH+/Pl45JFH1DVnzpzB999/XwQrJyIiIvIft8sFd6pNvQ82YZne4FG8ePFcr9m6dSsmTZoEnU6nOS6ff/zxxyq9weFwYPv27T4F0URERETBxnHmPFJXrIdt+x7AZgfMJphbxMHSrT2MVSoiGITlTq+H0Zh7zP/ee+/BarVmea5MmTLpfeQiIiIKfH1ERERE/mbbvBMJf/sPbBu3pwW86qBdfa6Ob96JYBDWQa8vateuneN5s9ms3jdu3LiIVkRERERUdDu8iRNmANmlM7hc6rxcF+gY9N4FGWZ3/Phx1cJM8n6JiIiIQknqivXZB7weLhdSV65HoAvrnN679csvv6juDi+++CJKlCiRr9s4duxYeopEVjiimIiIiPzB7XKl5fD6wLZtD6yjh0Cn930/Naf4R+KjWrVqoSBxp/cufPLJJ6hQoQLefPNNfy+FiIiIqGDZHXdyeHMj18n1AYw7vfm0YcMGzJkzB4sWLUJsbGy+b0eexXA3l4iIiAKOyai6NPgU+Mp1cn0e5BT/5LQLnF/c6c2HxMREPP7443jnnXfQo0cPfy+HiIiIqMDp9HrVlswX5pZxeUpt8IfAXl2AkoC3Xbt2eP311/29FCIiIqJCY+nWHsgtmNXrYenaHoGOQW8evfHGG0hJScGXX37p76VQEIy7lgJHeUUgkDqOrFixAsOGDUPdunX9vRwiIgpwhkrlYWrdNPsL9HpEjRsWFAMqGPTmgQS669evVyOHDQaDv5dDAe7HH39EfHw8ZsyYgUDx97//HX/84x8xc+ZM2Gy2Ir3vEydO4Omnn0avXr3yfRt79+5VAXvNmjVRrlw53HPPPfjuu+8KdJ1ERJTGdT0etz/5BvaNO+DFZIS5bQsUe/03MLduhmDAoNdHEyZMUG8SyGQ1fU1GEc+aNQvhwuVyI9nmUO8p+zSYNm3a4LXXXsvy/Lp164p8TZKS8+9//7tI7/PgwYPqsZCdZXniKK+U5MeSJUvQqlUr3HfffThy5AguXLiAkSNHYtSoUXjllVcKfN1EROHK7XYj9ZdtuPn2x3AcOnbnhE4Hc9d2KP7hnxH76TuIGhscO7weYd29IeMf39TU1Gyv+/bbb/H+++9j3rx56jrPtfKPQm5D/qj/4x//wLhx4xDqjp2Px+z1h/HTnrNIsTsRYTKgQ1xlDG5fF7Uq5r+LRShq3ry56uWcFZfLheeeew579vjW/7AglSxZskjvT8Z4S7BbqVKlfKd6nD9/Hg899BAGDBig+mJ7yMdS/fvhhx+ibdu2GDx4cAGunIgo/LgSbiNpylzYd2g7K+hLxsI6dihM9Qq2d25RCtug1+l0ql1bj7lz56rdIqNR+5D897//VcGJBLg55UBGR0eroDiUrdp5Gu/P2Axnht1dCXyXbz+lzr02rDW6NKvq1zUGC0kzkJfq/SHzv/HCVrVq2r8J2fXOr//85z+4efMm+vTp43VOAt///e9/ahd70KBB0Ad49TARUaCy7dqPpMlz4E64rTluvq8lrMP6Qxfp/Up3MAnLvw4ffPCBSlH4wx/+kH7s//7v/xAZGYnPP/9cEwg/++yzKuDNjfyxlR2tUN7hzRzwZiTH5bxcRzkbP348/vznPyPcZJUW5KuVK1eq9xUrVsyyl2P58uVx+PBhbN++/a7WSEQUjtzJKUicNAuJn0/SBLy6YlGIemYUosYMDfqAN2yD3ldffRV2u10Fsxnf5JgEuRkD2czXZPc2adIkhDJJacgu4PWQ83JdsNi/fz/GjBmDpk2bqqCpQYMGeOutt5CUlJR+zdKlS1XRok6nS3+rXbu25nZatGiRfk5efhdSJDZt2jR06tRJk/YiaTKyy+t5IiW3JW+Z299t3LgRDz74oFpb2bJlVaqETACU3PHCUrx48fTvQ3aD27cv2PYzcrv5df36dfVednuzUq1aNfV+y5Yt+b4PIqJwZD98HLfe/RdsG7ZqjpuaNUTMGy/B3Kzgh0T4S1gGvZQ3UqwmOby+kOuCobht8eLFKge0S5cu2LlzJ86ePYvHHnsMb7/9tnoZ/tq1a+q6nj17qpxtT1AlnQekiCoj2V3s27evSoORFJeLFy/iySefVJ9LsVrGVwqkqC3j1x89elS9/e1vf9PsBMv9PPHEE9i1axdOnz6tgueXXnpJrSc5OblQHhNJ95HA9OWXX8bly5dVp5LNmzejfv36eXrr2rVrga9N8oFzCmo9j/HVq1cL/L6JiEKR225H0qyFuP3R13Bdu3HnRIRF5e5GPT0K+phohJKwzekl36U6nCp31xdynVwfaQ7cf1oSlErV/wMPPKB2eoXsbEpOtxSWya697M56crTr1KmjCrDkWgmqMu9YSn64BMZS8Chk13jixIkqkJX2YHkheb7PPPOMSr3p3bt3elqA7PJKcL5q1Sr8/ve/x6effoqCdOrUKfUqx5QpUzBixIj0461bt1bfm7/JrvfatWvVz0Ye05iYGM35c+fOqfdms9lPKyQiCh6O0+eQOGEGXOcvaY4b69WEdcxQGEqVQCjiTi/lymI0qC4NvpDr5PpA9tFHH6mXyyWQysyT5z1//nxs27Yt/fjw4cNRpkwZbN26VQWfGS1YsAAdOnRAqVKlNMfl+ryS4FpSI7JamwS74osvvlCBe0GRQLt79+6qA0LGgDeQyBMB2ZWX73vgwIFq91vIrvlTTz2lujsIz448ERF5czudSF60Ggnvfa4NeI1GRA7th+gXHw/ZgFcw6KVc6fU61ZbMF3KdXB/Ipk6dqt5Xr17d65zk9daoUUN9vHDhwvTjFotFpRt4OglkJO24JCjLzGQy5WldiYmJ6bvLWa2tW7duah2Sey7T3grChg0bVDqC7Ex7dpYDkTyWkpIiu+fyhEXGgMt6JW+6Vq1a6ekNd9MhgogolDkvXUXCB18i5cel8hJl+nFD1UqI+ePziOjWHroQ734TuK9BU0CRPrzSliynYjaDXqeuC2S3bt1Kfyk8u8IqCXxlephnN9FDdhSlH7MEzdIBREYMnzx5UqU8SBrA3Tp27Fh6D+is1iaBnwR4UoCXeW35IUVxsrMr+buygy1DHzLbtGkTHnnkkTzdbuXKlbFmzRoUNEldkC4r8paR7PyKe++9N8snC0RECPdBE2s3Inn2IsBmv3NCr0dE7/sR0bcLdGEyZZZBL/lEBk9IH97s2pZJwCvnA31ARcbODBL8NmnSxOsaCWZF5rxR6Tfbv39//PDDD6rYTAq+vvrqKxUMF8baMt9/TmvLD8lj/uabb1RxnHwvzZo1Q8eOHTXXSCAphXaB6sCBA1i0aJH6OLvJd0RE4cp146ZqRebYry3A1pcrjahxw2GsUQXhJLT3salAyeCJz37TDd1bVEvP8ZX38rkcD4bBFJJnW6xYMfXxoUOHsrzG0xYsLi7O69xvfvOb9KElsisrvZwLKg+2Zs2a6R/nZ235IakNf/3rX9XtDh06VHWxCCbS0ULWLt0ussqDJiIKV7bNO3HrnY+9Al7L/fch5k8vhF3AKxj0Ur52fH98exDmvTNIvQ+GHV4P6bnrmer1/fffZ3nNmTNn1Evp0t0hqyBR2nJJKoLs8MouaUENJZF+vLKzmtvapDtE5h3ZuyEFcjK+V9IcJHDMOJ47kMkgGemjLDnYod4nm4jIV67bibj99VQkfjMN7qQ7v891JYoj+sXHYH1oAHRh2umGQS/lixSrSVuyQC9ay8qf/vQn9dK+9KCVIRAZSX6u5LdK+zJPKkFmngEm0qLs6aefzvZ+pAuDkMKzzGT6n8gcYEqfYE/PXMkXzmj37t2qS4F0eJDgXbhcLtVfuGHDhpg8ebJP379nXRkHXUghm+QySx9cub3C4Lk/afGWnYSEBPVkQwbD3L6tHYOZ0XfffYcXXnhBjQZfvXp1vjplEBGFGvveQ7j1ziewb92tOW6+tzli3ngRpgZ1EM4Y9FLYady4seq4oNfr8dBDD6kBEJ6AV1qTSZeEnMYES7/e6Oho1UJLdn2z4ynmkkAyY76uZ3Sup3uCnJN2YUJ2jt98802VOiGBnxTUCQmAR40apQLSxx9/PP12duzYofKLJbdVzvlS4CZBorh06ZLqzyvk+5E0ByGFerKLXdA7vp4nGNJmzDP8IzPZuZWAX/Km5ePMxRjyeMnPbPTo0SqtRB5btikjonDnTklF4pS5uP3vCXDfTEg/rouyIurJRxD16HDorWmbLWHNTX7TsGFD9Ub+sW7dOnfv3r3dJUqUcNepU8fdunVr9+eff+52OBy5fu0zzzzjnjlzZpbnrly54i5btqxU+6W/FStWzP3DDz+kX7Np0yZ1n2XKlHG/9NJL7oSEBM1t/Pjjj+6OHTu6Y2Nj3fXr13d36NDBPXXqVK/7unnzpjovt1W+fHn3rFmzclx35cqVNeuKiIhwP/HEE+qc3FfGc9HR0e7x48e775asKfPjYbFY3F26dPG69uLFi+569eq569atqz72ePfdd9V6atas6R43bpx748aNd70uIqJQYD960h3/x/fd15/8veYt4d8T3M74m+5g1bAQYiSd/MffgXe48uz27du3z99LoRAgaQFS2CUji4mIKLS57Q6kLFiBlKVrZQfzzgmLGdZh/WBu1yrb1pzhGiOxZRlRCIiPj1edF9q2bevvpRARUSFznruIxPHT4Dyrnc5prF0d1rFDYSijnRBKaRj0EgU5meQmeb6SFyxdJ4iIKDS5XS6kLv8JyfOWAY4MRcFGAyIH9ICle4eQn6p2Nxj0EgUxGRzx8ccfq/G80smAiIhCk/PKNSRNnAnHUW1nH0PlCqpQzVCpvN/WFiwY9BIFsdq1a+Ozzz7z9zKIiKiQSOmVbf0WJM1cAKSmtZxUdDpE9OyEiH7doDMxnPMFHyUiIiKiAOS6mYCkybNh33NQc1xfphSixg2DsRZbNuYFg14iIiKiAGPbtgdJU+bCnajt827ueC+sg/tAF2Hx29qCFYNeIiIiogDhSkxC8rR5sG3eqTmuK14MUaOHwNS4nt/WFuwY9BIREREFAPuBI0icOBPu+Fua46Z7msD68APQR1n9trZQwKCXiIiIyI/cNhuS5yxB6uqfNcd11ghYHx4Ec6umfltbKGHQS0REROQnjhNnkDhhOlyXrmqOGxvWUekM+hLF/ba2UMOgl4iIiKiIuZ1OpCxciZTFawCX684JswnWIX1VwVowjxEORAx6iYiIiIqQ8/wlJE6YAefpc5rjhhpVVSsyQ7nSfltbKGPQS0RERFRUY4RXbUDy3KWAw3HnhF6PiP7dEdGzI3QGgz+XGNIY9BIREREVMue1G2ljhA8f1xzXVyyXNmiiaiW/rS1cMOglIiIiKswxwr9sQ9L0+UBK6p0TOh0s3TogcmB36Ewmfy4xbDDoJSIiIioErlu3kTRlDuw792uO60uVgHXsUJjq1vTb2sIRg14iIiKiAmbbuQ9J382BOyFRc9zc7h5Yh/aDLjLCb2sLVwx6iYiIiAqIOzlFpTJISkNGumLRsI56EOamDf22tnCn9/cCiPztu+++Q9++fVGiRAlUrVoVr7zyCm7d0o6ADBSbNm3CyJEjUa/e3c1e37dvH37729+iZMmSOHnyZIGtj4gonNkPHcOtdz/xCnhNzRsh5s0XGfD6GYNeCmvPPfcc1q1bhwULFuDw4cMoXrw4PvzwQ4wZMwaB5o9//CMeeeQRTJkyBampGYoh8mjGjBn485//jE8//RQ3btxAUbh8+TKefvppVKxYEWazGTVr1sRrr72Gmzdv3vVtJyQkqNvr3LlzgayViCiv3HY7kmYswO2PvobrWvydExEWWMcNQ9RTI6EvFu3PJRLTG+hueg3C7gBMRuj0wfncacmSJfj888+xYcMGNfWmTJkymDt3Lnr16qV2QAPNX//6V/Tu3RsdOnS4q9sZNmwYhg4dimLFiiExUZtrVhguXLiAtm3b4tSpUzAYDHA6nThx4gQ++OADzJ8/Xz3pkMc+v1566SV1e7JLT0RU1BynzqpBE64LlzXHjfVqIWrsUOhLxvptbaTFoJfyxHHmPFJXrIdt+x7AZlfjEs0t4mDp1h7GKhURTCZPnqzelypVKv1Y7dq1cfToUQSqcuXKFcjtSJAfGxtbJEHviy++qHZiZ86ciRYtWqjd3b///e/45z//iYMHD6o0i6lTp+brthcuXIjp06cX+JqJiHwaI7xkDVIWrNSOETYZEflgb1g6tw3aTaFQxaCXfGbbvFM9m9X8z22zw7ZxuzonzbXNrZshWBw4cEC9l5fbg4WpAHs5Go2F/7//tWvXcOnSJSxbtiz9cZZddNnllePyxGP27NlISUlBREREnm/7N7/5jQqgn3/++UL6DoiIvDkvXkkbI3zyjOa4oVplRD06DIbyZf22Nsoen4KQzzu8XgFvRi6XOi/XBQtPPqvselLhkCK59957L8snFlIwKGw2W75ye5955hkV9DZu3LhA1kpE5EtqX8rqn3HrL59qA14ZI9yvG4r9/hkGvAGMQS/5RFIasg14PWSm+Mr1CGRJSUkqhUHezp49q45JAZTnmARg8jZhwgQVTE2cOBFnzpxBp06dVDpAxpfh5bbeeecdtGzZUqUdVKhQAc8++yyuX7+uzv/8889qV1OCankbO3asJj9XdjY959asWeO1VrmdF154AXXq1EG1atXQrFkztWNamAYOHJi+JnmT7+vixYv5vj15bNq0aZPlOfm+RGRkZJ5zeuXnIMVxks9LRFQUXDdu4van45E8bR5gt6cf15cvg2K/fxaR/btBZzD4dY2UM6Y3kE/PbFUOrw9s2/bAOnpIwOYxWa3W9Jzd6tWrq+IqCTjlY7F69Wr87ne/w86dO9Xn0iVh0KBBqsWXvAT/3//+Fw8//DDi4+PRtWtX9O/fX7URkzGTb775pnqp/aeffsLGjRtx3333qZfgR40apTouZO7EIAFy3bp1cfXqVa91SqAtBWvSmkxuS/KOFy9ejCFDhhTq4zNt2jR1XxLk/+tf/1LrE/Xr18/zba1cuRKVKmU/S97zfctjqM/Dv5fz58+rx09+bnn5OiKifI8R3rwTyd//qHrwZmTp0g6Rg3pBZ+YY4WDAoJdyJ10apGjNF3KdXG8JnjzZjO6//37s2LFDBay//PILvvzyS/UmgeDbb7+tuh4IKb4qW7Ys3nrrLc3u7bx587B37168//776nrZLe3SpYtX0CukL7AEk+vXa3fHXS4Xhg8frorMpEhLdpiFdG6QlADZXS4MEtTL/T7++OP45JNPNAGlFJwVtBUrVmjSHHz12GOPqcdddr+JiAqT63Yikqb8AHumjR9dieKqM4Opfm2/rY3yjtsklDuTUXVp8IlcJ9cHOek2ILp164ZWrVqpz7/99lv069dPteCSIHbw4MGar5EANy4uTn08a9as9OM57UZKC6/M5Gsl4B4xYkR6wJsx/aAwSE5tz549VTqC9O8tih3U//znPyrAlsfXV1988YVKhwjEPspEFFrsew7i1tufeAW85jYtUPzNlxjwBqHgj06o0EmqgrQlky4NuTG3jAvY1Ib8dDZo2NB7eo70lZVes7LjKm23Mrp9+7baFb6biW4SXIvWrVt7nSuM/sESxEswLZPeXn31VRQF+R7lsfroo498/ppjx46px1vSPYiICos7JRVJMxfCtn6z5rguOgrWRwbB3ILFs8GKQS/5RPrwSk5TjsVsej0sXdsj1J0+fVq9l7QHSTkoaFu2bFHvS5cujaLw8ssvY/fu3So3OTtS5JdXa9euzTKnVwZJ/OUvf1FT8GRAhi8k5UN2dz/++OMie1yIKPw4jp5MGzRxNa0g2cPUpAGsox6EPsa331kUmBj0kk9k8IT04c22bZler84H24CK/HA4HOq9jC0ujKDX00rNnqE6uDDJ2OUHH3xQ9cv929/+htdff93rmoIa2CF5ylLYJzu9UqTnK5maJ28DBgzIMcj2tJ+TwhMiIl+57Q4kz1uO1OXr5BfInRMWM6zD+8N83z1sbxkCGPSSz2TwhL5CWdWWTLo0pE9kaxmndnjDIeAV0ppMzJgxQxW0ZWXVqlWqgE3k9RelpDBIOy7pLFFU349MS5P1/vnPf1ZT02QUc2E8WZDOF9J5QQoF88JisWQbJEvrOOl2Ibm+HEVMRHkl/eWTZNDEOW17RmOdGrCOHQpD6cAbS0/5E/zJl3chOTlZFe1UqVJFNdHPjezsSXGRBAnSV1Qq3QN5ZG2h7fiOHYbYf72N2E/fUe/l82AMePO7GyitxDx9eP/97397nb9y5QrGjx+f/rkEY8LTvzfj/SckJKiPpTewR9u2bTXdDQpy7dlp37692vGVNAIJTAv637XkQI8ePVrddla74+fOndM8BplJfrN0kMjqbdKkSV7XEBH5NGhiyRok/P0zbcBrNCBySB9E/+4JBrwhJiyDXgk0pKVUjRo11E6dZ0hBTqRHavPmzdVAAfmjKkFyTEyM2hWTSvtwI8VqOos5qIvWpJBKeALPjCT48+wiZlarVq30l9lleMTTTz+N7du3q93G+fPnq7ZnDz30UPr1ntZa0ld2//796mP5N/fII4+oHV0hQaYnkPWM1JWc123btmX7hC2jzz77TPXUlfZfnrXnxhNketI15H6lmE3SK+T7y8+UtKzI7UvA26NHD/VEMTN58iD3m3Fq2+9//3sViEt/ZCKigua8cg0J//wSyXOXyLPy9OOGKhUQ8/rziOjeMaj/vlE23GFo1qxZ7p9//tn9zDPPSJSh3k6cOJHt9UeOHHFHRUW5mzZt6nY4HOnH7Xa7u27duu7SpUu7r1y5kud1NGzYUL1R0Vu1alX6z/6NN97QnEtMTHTHxcWpc127dnUnJyd7ff3FixfdderUSb+NjG/PPvus5lr5d1KjRg11zmAwuKtWreqOjIx0L1y40N2pUyd13GQyuUeOHOlOTU1VX/Pqq6+q4xUrVnRv2LBBHbtx44Z7+PDh6rhOp3MvXbrUvWvXLnUuOjo6/f6nTp2a6/e/f/9+dRty/XfffZd+/Nq1a+nH27Rp4z537pz7bqSkpLgHDhyovu9SpUpp3kqWLOm2WCzqvv74xz+mf438v+T5Xn7zm9/keh+rV69W18pjSUSUE5fL5U5Zu9F9/fk/u68/+fs7b0/9wZ00d4nbZbf7e4lUiDFSWAa9Hlu3bvUp6O3Xr5+65quvvvI6989//lOdGzduXJ7vn0Gvf9SrVy89sPO8lStXzv3EE0+4Z86c6bZarZpz8oRn0aJFXrcjAeJvf/tbd+XKld1ms1nd7ieffKJ+qWa2Z88eFURGRES4W7Vq5V67dq063qVLF/fQoUPdO3bs8Pqab775xt24cWN1282bN3c/9dRT7mXLlqn19erVy/3hhx+6Dxw4oK6VYFiefDVq1CjXQHHEiBEqCPV8f3q9Xt2PeOCBBzTfuwTj8rjk17Bhw7J8YpD5TYJwD3n8BgwYoIJieXKSGwa9ROQLZ/xN961Px2uD3Sd/747/0wdu+7FT/l5e0HA6Xe6kVLt6X5gKI0bSyX8QpuQl5Tp16qS3UfKMos3o+PHj6uVsz8eSEpGRvPzauHFj9dKsjEeVHq2+atSoUfptEBUEmaQmaQkyEpmIiNLYtu5G0tQf4E7UpqxZOrVB5OA+Kl2PcnbsfDxmrz+Mn/acRYrdiQiTAR3iKmNw+7qoVVE7SKkgFEaMFNYJK54BBDlZtGiReh8VFeUV8AoZIyuV5ZIfOXfu3EJZJ5Gv5N/rsGHD/L0MIqKA4EpMwu3/fY/Er6dqAl5dbAyiX3gU1ocfYMDrg1U7T+O5/6zA8u2nVMAr5L18LsflfDAI66DXF8uWLVPvpcNDdmNkK1eurD7evFk7vYWoqMgLNrK7K8ViDRo08PdyiIj8zr7/MG698wnsW3Z5td+MeeNFmBrV9dvagm2H9/0Zm+F0ZZ0YIMflvFwX6Bj05kJSGoQnsM2KJ6Xh0KFDRbYuooyeeuoptGvXTnVvICIKZ+5Um0pluP2v8XDH3xkJr7NGIuqJhxH12EPQR1n9usZgMnv94WwDXg85L9cFOg6nyIX0XBU5jUuV9IaMk7Ty4tixY+l5K1lhvi/54quvvvL3EoiI/M5x/HTaGOHLVzXHjY3rIWrUYOhjY/y2tmDkcrlVDq8v5LpXhrSCXu/7QKac4h+Jjzw1VQWFQW8uPAMFrNbsnxV6+qKmpKQU2bqIiIgojdvhQMqClWrYhGaMsNkE69B+MHdozTHC+ZDqcKbn8OZGrpPrI82BG1oG7soChHRl8DTvz47dbk8fH5tX8iyGu7lERET54zx/CYnjp8N55rzmuKFWNUTJGOGypf22tmBnMRpUlwZfAl+5Tq7Pi5zin5x2gfOLQW8uJF9XpnJlnoCVkWdyVenS/B+LiIioqMYIp65cj+QflsnoxzsnDAZEDugOSw9OVbtber1OtSWTLg25kevyktrgDwx6cyGV8DJe9tKlS9lec+3aNfU+q5ZmREREVLCcV68jaeJMOI6c0Bw3VCoP67hhMFap6Le1hZrB7euqlmQ5FbMZ9Dp1XaDjU6BctG/fPn14RVYSExNx9Wpawnz37t2LdG1ERETh1p4xdcMW1YpME/DqdLD06IRi//cbBrwFrFbFWLw2rLUKbLMix+V8YQyoKGjc6c3F4MGD8cYbb6hpaxcvXkT58uU15/fs2ZPewaFLly5+WiUREVFoc91KQNJ3c2DfdUBzXF+6JKJkd7d2dZXyIC3LYDIytaEAdWlWFdXKxhTpRLbCwKA3Fw0bNkTv3r2xePFi9TZu3DjN+RUrVqj3Y8aMQXR0tJ9WSUREFLpsO/Yi6bu5cN9O1Bw3t28N69C+cF65plqV2bbvAWx21bXB3CIOlm7tufNbwDu+rwxppbo0SNFaoOfwZqZzy2sFYergwYPp06vk43r16mV53eHDh9GsWTM0bdoUv/zyS/pxKW6T6sKEhATs3bsX5cqV8/tcaSIiolDhTk5B0rR5sG3crjmui4lG1OghMMXVh23zThXw4tf2oRp6vdoFlilsFFwKI0YK271/p9OJH3/8Mf3zuXPnZtuarG7duhg/fjy2bt2K3//+90hNTVWpDsOHD0d8fLy6nbwGvERERJQ9+8GjuPn2x14Br6lFHGLefEkFvI4z57MPeIXLpc7LdURhGfR+8MEHiIiIwB/+8If0Y//3f/+HyMhIfP7551l+zUMPPYQ1a9Zg27ZtqFixItq0aYOqVati//79uO+++4pw9URERKHLbbMjafp83P74f3DfSGsJKnSREbA+OhxRTz4MfXSUOpa6Yn32Aa/Hr63NiMI6vcHfmN5ARER0h+PkWSROmA7XxSua48YGtVU6g77knYIpKVqL/+2baTm8uTGbEPuvt1ncFuYxEgvZiIiIyK/cTidSFq1GyqJV2p1bkwmRg3vD0qmNd8Bqd/gW8Aq5Tq63mAt24RRU+JSHqJCcPn1a5YDLVL+TJ08i2MmLQtKtZNiwYSrPnYioIDgvXkbC+/9FyoIVmoDXUL0KYv70AiLuvy/rHVqTUe3g+kSuk+sprDHoJSoE3333HUaMGIH3338f169fRyj4+9//jj/+8Y+YOXMmbDZbkd3vqlWr0KlTJzXxsH79+vjHP/6hClHzatOmTWqATMmSJREbG4sBAwbgwAFtv08iKjqSnpCyagNu/eVTOE+evXNCr0fEgO4o9trTMJQvk+3XSyAsbcl8YW4Zx9QGYtBLVBhGjhypCh+lYDJUvP766/j3v/9dpPf5xRdfoFevXnjuuefUVMTVq1erTipDhw7NU+C7cuVKNV1x9+7d6vObN29i/vz56pgMniGiouW6Ho/b/xqP5Onz09IOfqWvUBbF/vAsIvt2hc5gyPV2pA+vBMk50uth6Zo2XZXCG4NeyheX24UkW7J6T1kzmUxqVzGUFOX3I51SfvOb3+Dxxx9XKRWiQoUKKvCWFoMffvihT7dz+/ZtvPbaa1i+fDkuXbqkdt7l6+UJiXz8z3/+s5C/EyLSjBHeuD1tjPDBo5pzEpjGvP48jNUq+3x7MnhC+vBmG/j+2qeXAypIMMGF8uTA5cP4ZssULDm0Esn2FESaItCrXlc81uoRNCjLPM+sAt9QYjQW3a8MyYeW3dzMUxB79Oih+mK/9dZbeOKJJ1CiRIlcd3mnTJmiUiM8HnjgAbz66qt49913cejQoUL7HojoDlfCbSRNmQv7Dm01vnRksI4dClO9Wvm6XRk8ITvE0pbMti3DRLaWcSqQZsBLHgx6yWfz9i/BqwvfhMN152VlCXzn7l2I+fuX4IO+b2NAw15+XSOFBhn+Irm8shvbokULr/P33nsv5s2bp3Knn3/++Rxva+DAgVkel17bQqYtElHhsu3aj6Tv5sB967bmuLltS1iH91c9eO+GBLbGscNgHT0kLV3CZGQOL3nhvwjyeYc3c8CbkRyX83JdoJOJep9++ql6qVy6Ksg46d/97ncoXbo0ypcvj/feey/92l27dqFnz56Ijo5WY6olDzQrsiMpL7fHxcWhVq1aqujq448/znYNO3bsULuNMtpa7leuf/HFF9VI64wkj1WCOk+/QslpbdWqlRqk0qRJE8yaNSv92gULFkCn02neZD0ZjRo1Kv2c3GdBKV68ePrtym6w5MrejUWLFqmXQWUAjCGLvD7P9yWPx90E1tJZQ1IoiKhwuFNSkThpFhI/n6QJeHXFohD1zChEjR161wFvRhLo6ixmBryUJe70kk8kpSG7gNdDzo/fMhUf9H0LgWrOnDl4++230wuaJADu06ePan4tHQmuXbumpvPVrFlTBVxS7S8vn8uI6sOHD2PIkCE4ePCgJmC02+3pnQAk8GzcuDGOHDmCfv364dSpU15rWLduHbp164Ynn3xSrUeCu2effRb/+te/cOXKFfVSvJBA/KuvvkJiYiKqVaumCuPkNmU9KSkp2LNnj8p1leulU4ScW79+vbptOS/FXjNmzNDc9+TJk1WQL+3U7iZgzExGcXfp0kWtWQreJPd38+bNGD16dJ5up1KlSiodQZ4UiMqVs87tk+4LYufOnflar/zMvv32W/X4yJMfIip49sPHkTRxJlzXbmiOm5o2hHXkg9DHRPttbRSe+FSIciXFapLD64vFh1YEdHGbBLgbNmxIz7X9y1/+gldeeQWXL19WAe+gQYPUcSlukh3fjRs3qgBRCqBkV1YC4+nTp2tu809/+hOWLFmiXmqXgFfUqVNHdR7IiuwAS9A1ePBg6PV6tZMpgaJnh9Pjo48+UgVXQoJY2Z3eu3cvzp49qwJvuQ8JmGWnMikpSV3Xrl079f2I7FqlSfAuI7ijotLGeN4tCewlaJfgWx43T7Fb69at1Trz8iYBr5DgP2Nwm9XOsrh69Wqe1yvdGmSXXXZ6iajgue12JM1ahNsffa0NeCMssI4ZonZ4GfCSPzDopVyl2FNV7q4v5Dq5PlBJjqikKsjL2kKq+vv27as+lkBYdiqFBLqyC+hJK5Ag6+GHH04/53Hu3DkVxMrL7Zlf0u/cuTNiYmK81lC7dm21BglaPTw7mtJKKyNJlRASbH/99dfpO8ySajFt2rT04DZj2sXLL7+sbl92ciXAzUh2uCXY69+/PwqCBOGyGy6pHbLbXFDkCYiwWq1ZnpcnC54nA75yuVx444030LJlS/XkQh6brl275piGQkR54zhzHrf+9h+kLl8nrRrSjxvr1kTMGy/Cct89Kg2KgpPd4cLcDUcwbc1BBCMGvZSrCJNFdWnwhVwn1wc6iyVtjcWKFdMcr1ixYnpwbDZrx1VK7q2Q9ACPqVOnql1b2dXMTH6xZ9VZ4IMPPlCBqifQ/emnn1RbLiE7t1kFdxI8ewJ1Dynwuueee9Lbe3nI7qh0NZAgL3Nbr2+++Ub1EC6IrhKyYy5B48SJE9G7d28UJM9jn/nx8PAMx8hLCzV5LN955x21Uy5pKA0aNFDHpYuD7DIT0d2NEU5etBoJf/8MrvOX7pwwGhE5tC+iX3ochlI5d1qhwOV2u7F+7zk8/vFSfD5/Jyat2IdLNxIRbBj0Uq70Or1qS+aL3vW6qeuDlSfIzOlcxkBsy5YtmoDYV1LsJWkSUiQnO7LSOis/PEFvfHy85rgUxcl9TJo0SaVmeAJFSUF47LHHcLckx1l2diUtZOvWrdlOQJNd7by8ye64kIJCIfnMWfF8v3l93IWkk8juvqSuSO62FCFmzn0mIt85L11Fwj+/RMqPS6WqN/24oWolxPzxeUR068DCsiB28Mx1vPzlGrz93c84f+12+o7v+KV7EWxYyEY+kT680pYsp2I2o96AR1ulpQCEixs30vLVZLfXV5KHKqkSsgssI32zSoHwlWenM/OOtRThDR8+XAW5kgv817/+FT/88INKl/CkbNwNCahl11iCdkmnkLZfHTt29GordvSotvm8r6QzhchuWponkJc86/ySx13ynyUfmVPZiPJONgBs6zYhadbCtN64GccI9+6MCB+nqlFgung9EeOX7sHqXWe8ztUsXxw9WlZHsOFTL/KJDJ6QPrwS2GZFjsv5cBtQ4Qk6s+rSkJ0HH3xQdV6QTgp3E/BmLPjyvFSfkaeg7b///a+aSiZBakHs8npIaoME07LrK50iJG2goEgwLaQjRlYpDp5gWgoT74YnD5sdHIjyxnXjJm5/OgFJU3/QBLz6cqVR7LWnETmgBwPeIJWYYsfXi3fj0Y+WeAW8JYtF4OXB9+DzF7qjZZ1yCDYMeslnMnjihzGT8WDjfuk5vvJePpfj4TiYom3btur92rVrc9zt9QRuUvglObwSZGXOGfaQXNzsclizCv5k17VXL+/HXnZfpchMdqP//Oc/4+eff8ZDDz2Egp6aJl0oJM1Bgvm8FJblRKandejQQfUt9rQvy/j4SDs02d32FCHml6dwUFrOEZFvbFt2pY0R3q8tlLV0bouYP70AY42qflsb5Z/D6cIPPx/FmA8WYcbaQyqFwSPCZMDobg0x8ZXe6NWqBgz64CxGZNBL+djxfQu7X1qHPS/9pN7L58G2w+sJUGWXMquAU/I8swtcMwa3MuxBOjtIysJnn32W7f15it88ty9DL5YuXao+vnDhgnqJ3UNuK3NbNNnR9XQ08JBdUAm2ZUxvdv1spUhLfPLJJyoozbizLGuRnd+GDRuqXWdfeILvjI+bFLLJTrPkNxfkTrK0jJM8aslLzkg6L8jjIf2UM7Y0kwBZWpFJ2znZ2faF3Lb0Em7evHmBrZsoVLkSk3D766lI/N/3cCfdKejVxcYg+rePwTpiIHTZPJmnwOV2u/HzvnN44uOl+GzeDtxMvLPJIo02et1TAxNe7Y1R3Roh0hLcWbEMeilfpFjNao4MyqI1GRzh6dEqO4YZyQAIITuXx48f1/xSkB1asX379vS+uNJRQYYcyG6rtD8bP368CiYlKPzHP/6R/pK/9NuVwinZwZR8W7lGOh7Ix9J+TAZKeAJSCSAz3rcnCJfgzJPOcObMGfW55LRKP9/syE6vZ8zuo48+qjknO6iyXgmeJVjN2IotO56BFpJT60npkPZokubg6Wbx1FNPFciO73333af6KEu/Y0//Xhki8swzz6ghIbLLnJE8iZAhGZK77HlCIaSTRZkyZVRxnzxuQtb31ltv4datW/jyyy/veq1Eoc6+9xBuvf0x7FvTBvt4mO9trlqRmRreacFIwePQ2et45au1eHPyzzh7VbtZ0KJOOXzxQne8POQelI6JREhwk980bNhQvVHReeONN9xms1m2bNWbTqdzN2rUSJ1r0KBB+nF5MxqN7nHjxrn379/vLl26tOZcdHS0e9myZem3u379enfnzp3dkZGR7urVq7sfeugh94oVK9y1atVy33vvve4//elP7rVr16prt27d6r7nnnvUtW3btnVv27ZNHX/55ZfdxYsXd7/99tvpt3vixAl1f9WqVXN/8cUX7rp167qrVq2qPn/ttdfct27dyvV7lq+rXbu22+VyaY7fvHnTXb9+fXedOnXc5cuXd8+aNSvH26lcubLmMYiIiHA/8cQT6lxsbKzX4zN+/Hh3Qfj+++/dzZo1c9esWdPdvHlz9zfffON2Op1e1128eNFdr1499RjJxx7Tpk1T37/JZFJrbt26tXvkyJHuJUuWFMj6iEKZKyXVffu7Oe7rT/5e83bjpbfdqVt3+3t5lE+XbiS6/z5to7vb72d4vT3+0RL35oMXvP5mhEKMpJP/+DvwDleeKnrZvSLKysmTJ9VAChlDLB/nh+x0ym14pr5lRdICXnrpJXTq1OkuVktEocRx7BQSJ8yA64o2tcrYuB6iRg+GvvjdFeKSf4rUpq05iDnrD8OWIWdXlIi2YGyPxujZsjoMBn1IxkjBnZxBRDmS3NZZs2bl+EtDet5KGoanKI+Iwpvb4UDK/BVIWbpWM1UNFjOsQ/vB3L4Vp6oFGafThUVbTmDS8n2IT9ROTbWYDBjSoS6GdaoHq+XuBxcFMga9RCFM8mGlf65n0lxmMvxBpsHJ5LbsukkQUfhwnruIxPHT4Tx7QXPcUKsaosYNg6GMdjIkBTZ5MX/TwQv4atFunLmSoDknz1u6t6iOsT0aoUzxrEe+hxoGvUQBzNMpIjVV+8w8OzLsQroeSBGYTBt78803sW7dumzbnX388cf429/+hrp1g6v7BhEVLLfLhdTlPyF53jLpXXXnhNGgeu5aunOqWrA5eu4Gvly0GzuPXfY616xWWTzVpwlqVwqv0dAMeokCmCdglW4SMtAiLi4ux+ulk4R0l5A3If15W7ZsmeW1MvY3pzZrRBQenFevI2nCDDiOausGDJXLwzpuOIyVObwlmFy5mYQJS/dixY5TmuwUUbVsMTzZpyla1ysflikqDHqJApD0nJX2Zp7xuNLiTHrJ9uvXT7Xkys5vf/tb1S83KipK9eiV9l5ERNmOEd6wBUkzFgCpGQbg6HSI6NkJEf26QWdimBAsklLtaqjErJ8OI9Wu7TUfG2XB6O6N0EcGSwRAkZq/8F8zUQCSaWPnzp3L89dJX17poUtElBPXzQQkTZ4N+56DmuP60iVV7q6xdnW/rY3yXqS2ZOtJfLt8L27c1qbCmY16DO5QF8M71UdURGgXqfmCQS8REVEYsW3fg6Tv5sKdmDZkx8Pc8V5YB/eBLsLit7VR3nbqtxy+iK8X7cbJS7e8zndrXg3jejZG2djwKFLzBYNeIiKiMOBKSkbytHmwbdqhOa4rXgxRo4fA1Lie39ZGeXPsfDy+Wrwb2494v7LXpEYZPNW3KepWDq8iNV8w6CUiIgpx9gNHkPjtLLhv3NQcN7VsAuvDA6GPjvLb2sh3V28l49tle7F020mvIrXKpaPxRJ+maNugQlgWqfmCQS8REVGIcttsSJ6zBKmrf9Yc11kjYB3xAEytmjJACgLJqQ7MXHdIvaVkKlIrHmXGqG6N0Ld1TRjDuEjNFwx6iYiIQpDj5Bkkjp8B16UrmuPGhnVUOoO+RHG/rY1843S5sWzbSUxcthfXE1I050xGPQa1q4OH72/AIjUfMeglIiIKIW6nEykLVyFl8Wrpd3jnhMmEyCF9YOnUhru7QWDr4YtqktqJi9qUFHF/0yp4tGccypdkWkpeMOglIiIKEc4Ll9PGCJ/Wtjw01KiSNka4XBm/rY18I0GudGSQzgyZNa5eWhWp1a9S0i9rC3YMeomIiEJhjPCqn5E8dwngcNw5odcjon83NWxCZzD4c4mUC0lfkF67S7acgCtTkVrFUtF4oncTtGtUkbv0d4FBLxERURBzXY9H4sQZcBw6rjmur1AWUY8Oh7FqJb+tjXKXYnOoKWrT1x5Eik1bpFbMasaorg3R795aKoeX7g6DXiIiomAdI7xxO5KmzQNSMkzi0ulg6dYekQN7QGdigVMgF6mt2H4KE5btwbVbmYrUDHo80K42RnRuoAJfKhgMeomIiIooBQF2h5TdQ6e/u107163bSJoyB/ad+zXH9aViYR07DKa6Ne9ytVSYth+9hK8W7saxC/Fe5zo3kSK1xqhQKtovawtlDHqJiIgKkePMeaSuWK/G/8JmB8wmmFvEqd1YY5WKeb492679SJo8B+6E25rj5nb3wDq0H3SREQW4eipIpy7dwteLd2PTwQte5xpWLYWn+jVV76lwMOglIiIqJLbNO5E4YYa2dZjNrtIS5Jx0VDC3bubTbbmTU5A0cwFsG7ZqjuuKRcM66kGYmzYs6OVTAbmRkIJJK/ZhkRSpZapSq1AyCo/3boIOjSuxSK2QMeglIiIqpB1er4A3I5dLnZeCs9x2fO2HjyNp4gy4rmlfDjc1awTryEHQF+NL4YEo1e7E7PWHMX3NQSSlOrRFapEmPNKlIfq3rQWzkZ01igKDXiIiokIgKQ3ZBrwe0mps5XoYxw7L8rTbbkfyD0uRunKDVK7dORFhgfWhATC3acHdwQAku7krd57GhKV7cOVmsuac0aDDgDa18UjXhohhkVqRYtBLRERUCEVrKofXB7Zte2AdPcSruM1x+pwaNOG6cFlz3FivJqLGDoO+ZGyBrpkKxq5jl/Hlol04cs67SE1SGB7r1QSVSnNn3h8Y9BIRERU06dIgRWu+kOvkeov5zhjhJWuQsmBlpjHCRkQO6gXL/ffddfcHKninL6cVqW084F2kJhPUZJKaTFQj/2HQS0REVNBMRtWlwafAV66T66V366UrKs/XeeKM5hJD1Upq0IShQtnCWjHlU/ztVExeuQ8LNh33KlIrX8KqdnY7NanMNJQAwKA3jy5duoT33nsPCxYswIULF2AwGBAXF4cnnngCo0eP5j9qIiJSO7HSlky6NOTG3DJODZRIWf0LkmcvAux27RjhPvcjok8XjhEOMDa7E3M2HMH3qw94FalFRUiRWgMMbFsbZhN/boGCQW8eHD58GJ06dUJsbCz+97//4b777kNqaiqmTp2Kp59+GkuWLFEfM/AlIiLpwyttyXIsZpPguFUz3P7XeDgOHNGeKlcGUY8Og7F6lcJfLPlMdnPX7D6Db5bsweX4JM05g16H/m1qYWTXhigeZfHbGilrDHrzQHZyr1+/jp9//hk1atRQx0wmE5588klcvHgRb775Jnr06IFx48b5e6lERORn0oZM+vBm27ZMr4el471I/N9UuJO0Y2gtXe5T+bs6M6v7A8meE1fw5cJdOHT2hte5do0q4fFecahcpphf1ka5Y9Dro2PHjmHTpk1o0aJFesCbUa9evVTQO3fuXAa9RESkyOAJ6cMrbcmkS4NnIpupSQO4E5OQuuYXzfW6EsURNWYoTA1q+23N5O3slQT8b8kebNh3zutc3col8FSfpmhSs4xf1ka+Y9Dro6tXr6r3sqPrdru9UhgiIyPV+2LF+AyPiIi0O77Sh1fakkmXBvuhY2ljhG8laK4zt2mOyOEDoLem/T0h/7uZmIrvVu7H/I3H4MxUpFY2VorU4tC5SRXo9UxrDAYMen1Uu3ZtGI1GnD9/HpMmTcKYMWM052UXWDzyyCN+WiEREQU0mx1JsxbC9tNmzWFdlFVNVZPCNwqcIrUffzmKKasOIDFF24HDajFixP0NMKhdHVhYpBZUGPT6qFSpUhg+fDimTJmCZ599FpUrV0bXrl3VuYSEBPzjH//A448/jj59+vh7qUREFGAcR0+q3F7X1eua46a4+rCOGgx9cb5KGAjkldy1u8/imyW7cfGGtkhNdnP7tq6J0d0aITaaRWrBSOeWnzD5RIrY2rVrh4MHD6pd388++wyjRo3CwIED1fE33ngjT50bGjVqpN7v27evEFdNRET+4rY7kLJgBVKWrtWOEbaYYR3WH+Z297DjT4DYe/KqKlI7eEb7xES0aVABT/RugqplY/yytnDUqBBiJAa9eSQ5vT179sTu3bvV59WqVcO7776rgt/8/EClQK5WrVrZXsOAmIgoODnOXkDShOlwnr2oOW6sXR3WsUNhKFPKb2ujO85fu43/Ld6Nn/Z6F6nVqRSLJ/s0RbNaHApSmIFtVjzxUUHGQUxvyKPy5ctj+vTp6N27N5KTk3Hq1Ck8//zzKFu2rAqGiYgovLldLqQuW4fkecsBp/POCaMBkQN7qv69HCPsf7eSbJiycj/mbTwKh1O7/1emeCTG9YxD12ZVWaQWQrjTm0dbt25VxWo//PADoqKi0K9fP+zZs0elO0yYMAEjR470+baY3kBEFFqcV66ljRE+dkpz3FClAqLGDYehUnm/rS1cBkekOpywGA3ZBqs2hxPzfjmmAt7bmYrUIs1GDO9cH4Pb10GEmfuC/lQYMRJ/onkgD3y3bt1UOkODBg3UsXXr1qFv375qYIX0561SpYqa2kZEROFD9o+kK4N0Z0Cq7c4JnQ4RvTojol9X6Iz8k1tYjp2Px+z1h/HTnrNIsTsRYTKgQ1xlDG5fF7Uqxqb/jCSFQVIZLlxP1Hy9xMe9W9fEmG6NUKJYhJ++Cyps3OnNgzZt2mDXrl24cuUKoqOj04/funVLBbo7d+5E48aN1c6vL7jTS0QU/Fw3byFx0mw49h7SHNeXLYWoscNgrFXNb2sLB6t2nsb7MzZ79dH1jAV+bVhrlC8ZpYrU9p+65nVN63rlVd5utXIsUgsk3On1o+3bt6tevE2aNNEEvCImJgazZ89G/fr1sXfvXhw5cgR16tTx21qJiKho2LbtRtKUH9R0tYwsndogcnAf6CwcI1zYO7zZBbxCjr83bROyOluzQnE1Sa1FnXKFvk4KDAx6fXT48GH1XnJ3s1KzZk107twZy5cvx7Vr1xj0EhGFMFdiEpKnzYNt807NcV3xYmljhBvV9dvawomkNGQX8HpkPlsqJgJjezRG9xbV1U4whQ8GvT6S1mRCdnGdTicMBu8pLBEREdDr9Wp6GxERhSb7/iNI/HYm3PG3NMdNrZrCOmIg9FFWv60t3IrWJIfXVxaTHsM7NcCQjnVVwRqFH/7U85DP26JFC5Xm8MUXX+C5557TnL906RLWr1+PYcOGoXTp0n5bJxERFQ63zYbk2YuRuuYXzXGdNRLWhx+AuVVTv60tHEmXBila89WXv+2BSqU5+S6cMej1kUzMkf68Xbp0wWuvvYbIyEjVukx2fKWNmQTBkuLw3//+199LJSKiAuY4cTptjPClq5rjxkZ1ETV6CPSxLIIqatKWTLo0+BL4ynUVSmrrcSj8MOjNA0lb2LFjB95//3319sILL8BkMqlgd8SIEWpIhcXCedxERKHC7XAgZeEqpCxerR0jbDbBOqQvzB3v5RhhP7lyMwmx0RZcvKEtIsyKtC/jkAliyzI/YssyIqLA5Tx/KW3QxGnteFpDzaqIGjcMhrJMZfOHxBQ7pq4+gLkbjsDucOV6vRSrffabbun9eik4sGUZERFRUYwRXrkByT8sBRyOOycMBkT074aIHh2hy6KYmQqXw+nCgk3H8d3KfbiZmGEASA48fXoZ8JJg0EtERPQr59XrSPp2JhyHT2iO6yuWQ9Sjw2GsUtFvawtX8oL0L/vP4+vFu3H26m3NOcks6dGiOu5vWhUrd57KcSIbEYNeIiIKe2qM8M/bkDRjPpCSeueETgdL9w6IHNAdOpPJn0sMS4fP3lCT1HafuOJ1rkXtsmqSmieobVm3HF4Z0kp1dZAiN+bwUmYMeomIKKy5bt1G0ndzYN+1X3NcX7okrGOHwlSnht/WFq4uxydh/NI9WLnjtNe5amVj8ESfJmp8cOYiQgl02YOXssN/GUREFLZsO/epgNedkKg5bm7fGtahfaGLYEeeoi5Sm7bmIOasPwxbpiI16dQwpnsj9L6nBgwGvd/WSMGLQS8REYUdd3IKkqbPg+2X7ZrjuphoWEcNhrlJA7+tLRw5nS4s2nICk5bvQ3xiqvYJiFGPIR3qYXjnerBamGJC+cegl4iIwor90DEkTZwJ1/V4zXFTi8awPjII+ugov60tHHOpNx+6iK8W7cLpywmac5K50K15NYzr2RhlinO0M909Br1ERBQW3Da7akOWunK95rguMgKRDw2E+d5mHDRRhI6ej1dFajuPXfY616xmGTzZtynqVCrhl7VRaCrSoPfUqVNYsmQJNm3ahCNHjuD69etwOp2IjY1F+fLl0axZM7Rv316N+tXrma9DREQFw3HqbNoY4QvaAMtYvzaixgyBviTbWhWVqzeTVZHaih2nNEPuRJUyxfBE7yZo06ACn4BQcAa9K1aswLvvvov167XPrjMOg5N/3PPnz1cflypVCo8++iheeeUVlC7NiTdERJQ/bqdTjRCWUcJwZSiMMhkR+WBvWDq3hY6bLEUiKdWOGWsPYdZPh5Fqd2rOxUZZMFqK1FrVgJFFahSMY4gPHDiAJ598Ej///LMmwC1XrhxKliwJq9WKyMhIpKSkICEhAVevXsW1a9fSr4uKisI777yDl156CaGIY4iJiAqP8+KVtDHCJ89ojhuqV04bI1y+rN/WFm5Faku3ncTEZXtx47a2SM1k1KsBEg91ro+oCBapUZCOIZ41axbGjRuHiIgIPPbYY+jQoQNatGiBmjVrqkA3Ozdu3FDB8rp167Bo0SK89tprWLhwIebMmYOYmJjCWi4REYXSGOE1vyB5zhLAbr9zQq9HRN8uiOh9P8cIFwHZ7Npy+CK+XrQbJy/d8jrftVlVPNorDmVjWaRGQbzT++9//xtvvvmmSmmQnV7TXUyxOXr0KN5++20V6a9cuRIlSoROUjt3eomICpZ0ZEicNAuOA0c1x/XlyyBq3HAYq1f229rCyfEL8fhy0W5sP3LJ61xcjdJ4qm9T1Ktc0i9ro+AQFDu9Uqg2adIk7Nq1C1WqVLnr26tduzYmT56MmTNnqjxf2fFlcjsREXmNEd60E8nTflQ9eDOydG2HyAd6QWfmy+eF7eqtZHy7bK9KZ8i8pVa5dLQqUmvbsCL/jpNfFGjQm5SUpHZ5ZUe2oFMRhg4dqnKAv/jiCzzzzDMFettERBS8XLcTkTTlB9i379Ecl44MaoxwvVp+W1u4SLY5MHPtIcxcdwgpmYrUYqxmjOraEP3a1GKRGoVO0Ltx40YV9BZW7m3fvn3x/fffF8ptExFR8LHtPoCkybPhvnVbc9zctiWsw/urHrxUeJwuN5ZvO4kJy/bieoJ2h91k0GNQuzoYcX99REea/bZGokIJeqW/bmEbMWJEod8HEREFNndKKpJmLoBt/RbNcV10FKwjB8HcvLHf1hYuth25hK8W7sLxize9znVuWgWP9YxD+ZKcbkeBIyAmsh07dgwffvihGl7RvHlzvPzyyyFVsEZERAXHfuRE2hjhq9c1x01NG8A68kHoY4r5bW3h4MTFm6ojg3RmyKxRtVKqSK1B1VJ+WRtRQAS9DzzwAE6fPq2S11u2bImvvvpKHd+5cyfuv/9+3Lp1SxUiLF68WBWtbd26FcWK8RcXERGlcdsdSJ63DKnLf5LKtTsnIiywDusP830tWSBViCR94dvle7Fkywm4MhWpVSwVhcd7N0H7RpX4M6CAVWRBb5s2bbBt2zZMnTpV9ewVLpcLY8eOxc2bN1UesLQmq1q1Kv7yl7/gvffew1//+teiWh4REQUwx5nzSBw/Ha7z2hZYxro1YB0zFIbSbH9VWFJsDjVFbfrag0ixaYvUikWaMLJrQ/RvU1sNmiAKZEUW9K5du1bt7noCXiEB8O7du9WzQhlm0a1bN3W8WbNm6NOnD4NeIqIwpwZNLFuH5HnLZbTXnRNGIyIf6KnakXGMcOFwudxYseMUJizdq1qRZWQ06DCwbR080qUBillZpEbBociC3i1btqBHjx6aY7KbKwGvtCPzBLyiRo0aOHfuXFEtjYiIApDz8lUkTpwJ57FTmuOGKhUR9ehwGCqW89vaQt2Oo5fx5cJdOHYh3utcpyaVVZFahVLRflkbUcAHvZLKkHH4m+zy7t+/HwaDAe+8847mWhlDnJiYWFRLIyKiQBs0sW4TkmYtBGyZxgj37oyIPl2gMwZEHXbIOXX5lipS23Twgte5BlVL4um+zdCwGovUKDgV2W+Ntm3b4pNPPsErr7yCvXv34qWXXlK7vKNHj0bdunXTr3M4HPjtb39bVMsiIqIA4oq/lTZGeN9hzXF92dKIenQYjDWq+m1toezG7RRMXrEfCzcfV2kNGUnbscd7xaFjXGUWqVFQ07kzbr8Woj179qh8XpPJhBs3bqidXxlTvGPHDpQsmVaA8NNPP6mAVzo6yP9Yzoz5WyGoMOZKExEFK9uWXUia+gPcSdr8UUvntoh8sDd0FuaOFrRUuxNz1h/GtDUHkZTq0JyLjjDhka4NMaBtLZiNBr+tkcJTo0KIkYpspzcuLk7l9X700Uc4efIk6tevr3Z9PQHv008/jUuXLqFatWrqjYiIwoMrMQlJ3/8I+5ZdmuO62BhEjRkCU8M7rwZSwZDd3FW7TmP8kj24clP7JMOg12FA29oY2aUBYqIsflsjUdDu9JI37vQSUbiz7zus0hnc8bc0x82tmyHyoQHQR1n9trZQtev4FVWkduTcDa9z0mf3sd5xqFyaffLJvwJ+p1d2cKXPrr4Q28fIfVSvXr3Qbp+IiAqfO9WG5NmLkLp2o+a4LsoK68MPwHxPE7+tLVSduZKgitR+OXDe61z9KiXVJLXG1Uv7ZW1ERaFAg17J1Z04cSLeeustFFbA+8033+Ddd98tlNsnIqLC5zh2CokTZ8B1+ZrmuLFxPUSNHgx98Ri/rS0Uxd9OxeSV+7Bgk3eRWrlYq9rZ7RRXBXo9i9QotBVo0Nu8eXNMmDABkyZNUl0ZCjrgHTNmDObPn1+gt0tEREXD7XAgZcFKpCxZox0jbDHDOrQfzO1bsTtAAbJJkdqGI/h+9QGvIjWrxYiHuzTAoPvqwGxikRqFhwIvZPvwww/RuXNn1ZXhzTffRGxs7F3f5syZM/Haa69h7ty5alwxEREFF+e5i0icMB3OM9r+r4Za1RA1bhgMZdj7taDIbu6a3WdUkdql+CSvIrV+99bCqG4NUZxFahRmCjzolZZkS5YswYABA1Tu7ZNPPon+/fujXbt2ecr1vXDhAhYtWqTSGY4fP64CXxlPTEREQTZGeMV6JP+4FHBkaENpMCByQHdYenTkGOECtOdEWpHaobPeRWr3NayIx3s3QZUyLFKj8FQoLcuKFSuGZcuW4c9//jP++c9/qt1f2aGVSjxpVSbBsFwjb9KvNykpSU1gO3PmjApwZSLb2bNn1W3dc8892Lp1KypXrlwYSyUiokLivHodSRNnwnHkhOa4oXJ5WMcNh7FyBb+tLdScvZqA/y3egw37znmdq1OphCpSa1qzjF/WRhQ2Lcu2b9+ugt/Fixen3WEu+Vqe5dSuXRtvvPEGHn744ULtBuFPbFlGRCE7RvjnrUiaPh9Itd05odMhomcnRPTrBp2JY4QLwq3EVHy36gDm/XIUzkxFamWKR+KxXnG4v6l0VWKuNAWXgG9ZlpUWLVpg4cKFavf2+++/V8GvTGez2TL8IvxVpUqVVD7wkCFD0K9fPxgMTK4nIgomrlsJSJo8B/bdBzTH9aVLqtxdY222nCwINocTP/58FFNXHcDtFLtXkdpDnevjwfZ1YWGRGpF/h1M4HA6cOnVKtThLTU1FdHS06u9bokQJhBPu9BJRKLFt34ukKXPhvp2oOW7ueC+sg/tAF8HCqbslf7LX7TmL/y3Zg4vXtY+z7Ob2bV1TFamViI7w2xqJwnanN8s7NRpRq1Ytf9w1EREVMHdyCpKm/Qjbxh2a47qYYqrvrimuvt/WFkr2nbqqitQOnL7uda5N/Qp4vE8TVCvLDkdE2WFSFRER5Zv9wFEkfjsT7hs3NcdNLePUZDV9dJTf1hYqzl+7jW+W7FE7vJnVrhiLJ/s0RfPaZf2yNqJgwqCXiIjyzG2zI3nuEqSu2qA5rrNGwDriAZhaNeWgibt0K8mGqav248dfjsLh1GYilo6JxKM9G6Nr82osUiPyEYNeIiLKE8fJM0icMAOui1c0x40N6iBqzBDoSxT329pCgd3hwvyNR/Hdyv1ISNYWqUWajRjeuT4Gt6+DCDP/hBPlBf+PISIin7idTqQsWo2URatk7NedEyYTIof0gaVTG+7u3mWR2vq95/C/Jbtx/lqmIjUd0LtVTYzu3ggli7FIjSg/GPQSEVGunBcvI3H8DDhPafNKDTWqpI0RLsfBB3fjwOlrqkht36lrXuda1SuPJ/s0QfVy3EEnuhsMeomIKOcxwqt/Vvm7sDvunNDrEdG/mxo2oWNP9Xy7cD0R45fswZrdZ7zO1SxfHE/2bYqWdcr5ZW1EoYZBLxERZcl1PR6JMkb40DHNcX2Fsoh6dDiMVSv5bW3B7nayDVNXH8APG47C7syQKgKo9IVxPRuje4vqMLBIjSh4g94ffvgBkydPVsMptm7dqo65XC489dRTqFmzJp5//nk1rCIYyFS5WbNmYf78+Th//jyqVauGRx55BD179vT30oiI7m6M8KYdSPr+RyAl9c4JnQ6Wru0Q+UBP6Ewmfy4xaDmcLizYeAyTV+5X3RkyijAZMKxTfQzpWFcVrBFRwSrS/6ueeOIJjB8/Xn1ssdyZzKPX6/HVV1/ho48+QsOGDTFz5kzce++9CGTLli1TgXrJkiXxz3/+E/fff7+/l0REdNdcCbfVVDX7Du0UJH2pWFjHDIWpHgcL5feJxM/7z+Prxbtx7uptzTmp/evZsgbG9miEUjGRflsjUagrsqD3m2++UW/ZkYrfl19+GUeOHEGPHj2wZcsW1K1bF4Ho73//O/74xz/isccew+effw4TdzyIKATYdu1H0uQ5cCdogzJzu3tgHdoPukh2DciPQ2evqyK1PSeuep2TfF0pUqtZIdYvayMKJ0UW9H7xxReoWrUqvv32W7WLW7FixSyvGzFihNr1feONNzBt2jQEmvfffx+vv/46xo4di6+//trfyyEiKpgxwjMXwLYhLeXMQ1csGtZRD8LctKHf1hbMLt1IxISle7Fy52mvc9XLxagitVZ1y/tlbUThqMiC3v3792Pp0qVo3759jtcVL57WkmXlypUINJJ28Yc//AGNGjVSQTwRUbCzHz6OpIkz4bp2Q3Pc1KwhrCMfhL5YcNRYBJLEFDumrTmI2esPq0ETGZWItmBsj8bo2bI6DAa939ZIFI6KLOiV4rTcAl6xefNm9T4xUduY299u3LiBZ555RuVlffLJJ5qcZCKiYOO225H84zKkrlgvCad3TljMsI4YCHObFhw0kY8itUWbj2PSin24magtUrOYDBjasR6GdayHSAuL1Ij8ocj+z2vQoIHqcJBdWoO4desW3nvvPfWLtn79+ggkb7/9Nq5du4amTZuiW7du/l4OEVG+OU6fSxsjfP6S90mXG46Dx2CoXAHGKtn/vqY7ZDNk44ELqkjtzJUEzTl53tCjRXW1u1u6OIvUiMIi6H3yySfx6quvYsqUKVmeP3PmDIYNG4aTJ0+qoPfRRx9FoEhJSUkvwuvfvz8+/fRT/Pjjj9i7dy/sdjuaNWum8nzzEwwfO3ZMpUtkZ98+bQU1EdFdjRFeug4pC1YATmfWF9ntsG3cDtvmnWrSmrl1s6JeZlA5cu4Gvlq4CzuPX/E617x2WTzZpylqV2SRGlFWcop/JD6qVatWcAa9Dz/8MObOnatSHMaMGQOHw4Hly5fj7NmzWLNmjep3K8Gl6Nq1K5599lkEUnuy27fTqpnnzZunOjfMmTMnPdXhnXfeUd+D5PlKcE9EFGicl66q3V3nCe+iqiy5XOp6GUTBHV9vl+OTVJHaih2nvM5VLVtMBbut65VnighRANG5JXIrIhLo/ulPf1I7pRLgen4ZeJZgMBhUGzAJJCMiAqc1juxQSy9eWVNCQgKMRu1zheeeey69ddnu3bt9Ts3wPMPhbi4RFRb5/Zq6diOSZy8CbPY8f725bQtEjR1WKGsLRkmpvxap/XQYtkxFarHRFozp3gi976nBIjWiu1QYMVKRBr0eV69exeLFi7Fnzx7Ex8fDarWiXr166NOnj5pqFmhkl/r7779X2+xHjx7Ncgu+Tp066o+LTJSToN4XDHqJqDC5btxE4qTZcOw/rDmuK1sa7hvxgN2R+42YTYj919vQ6cM7iHM6XVi89QS+Xb4P8bdTtQ+RUY/BHepieKf6iIpg33aiglAYMZJfSkhLly6NUaNG5XjNzz//jPvuuw+BQIJ0UaJEiSzPSzDcunVrbNq0Sb0REfmb5OQmff8D3ElpaWMeli73IaJPF9x85S8+3pA9LTi2mAtnoQFONjM2H7qIrxftxqnLt7zOd2teDeN6NkbZWKtf1kdEvgvIvinJycno3LkzbDZtyxd/iY1NK0LIaVNcpsdJwCutzYiI/MWVmISkqT/AvnW35riuRHFEjRkCU4M6cLtcagfXp3QHuc4UkH8qCt3R8/H4atEu7Dh62etc05pl8FTfpqhTKevNECIKPEX+m+zw4cM4cOCAKgxzZlE9LAGvdEbI6py/1KxZU70/d+5cttd4WrGVKlWqyNZFRJSRfe8hJE6aBfdNbdssc5vmiBw+AHprWsssSVUwt4hTXRpyY24ZF3apDVdvJmPCsr1Yvv2kpoWxqFKmGJ7o3QRtGlRgkRpRkCmyoFdydx966CHVsSE3sqMaSL9MevXqhX/84x+4ePEiTpw4gRo1anhd49nhbdOmjR9WSEThzJ2SiqTZi2Bbp02v0kVZYR05SAW4mVm6tVcpENKlIVt6PSxdcx8qFCqSUx2Yse4QZq07hBS7duOleJQZo7s1Qp/WNWFkkRpRUCqyoPfpp59Wrb+ETDOTHVHpdpCZ7PBeuHABrpx+ERexjh07Ii4uThXeTZ06VbUsy2rMcqD1Fyai0Oc4dipt0MSVa5rjprj6sI4aDH3xYll+nbQhkz688rVZBr56vTofDu3KnC43lv5apHY9QZsDbTLq8WC7Ohhxf4O7LlJzudxIdThhMRqg1wfOxg5RuCiy7g3FixdXKQ1///vf8dJLL2UZ8GYsYuvQoUNApThIH14ZPiHBurQlK1eunGa97dq1U2OKpXWZr9i9gYjyy+1wIGX+CqQsXes9RnhYf5jb3ePTK2aOM+eRunI9bNv2pOX4mk0qpUF2eMMh4N1y+KIaLnHykneRWtdmVVWRWrkSUXd1H8fOx2P2+sP4ac9ZtYMcYTKgQ1xlDG5fF7U4uIIo9FqWlS9fXgWxV654T63JSsOGDdXuaSD56quv1I71Pffcg8mTJ6s2axs2bFAtzZo0aaIGVuQUzGfGoJeI8sN57iISx0+H8+wFzXFj7eqwjh0KQ5m81xao4jbp0mAyhkUO74mLN/Hlwl3YdsR7FHNc9dKqSK1elZJ3fT+rdp7G+zM2q93kzAx6HV4b1hpdmlW96/shCjWNgrll2cCBA9VENl8FWsArZNqaFLX97W9/Q8uWLREZGanalcnAjccffzyg8pCJKPRIYJq6/Cckz1sGODK8EmY0IHJgT5Wnm9+AVX1dGLQlu3YrWaUxSDpD5ji0UuloVaR2X8OKBfL7XHZ4swt4hRyX89XKxnDHl6gIFFnQ++6772LJkiWqkK179+65Xi/BpAx9CDSS4iBvRERFyXnlGpImzoTj6EnNcUPlCoh6dDgMlcr7bW3BINnmUAVqUqiWYtOmzsVYzRjVtSH63ltL5fAWFElpyC7g9ZDzcp3s+BJRiAS9ZcuWxapVq/Dss8+iVatW6b1vs7J582acPKn9xU5EFI4kA822YQuSZiwAUjP0LtfpENGrMyL6dYUu02h00gaV0nps4rK9uHYrU5GaQY8H2tXGw/c3QHRkwe5yS9Ga5PD6Qq57ZUgrFrcRFbIi/U15+vRpOBwONcghOjo6y2tkIMWlS945VkRE4cZ1MwFJk2fDvueg5ri+bClEjR0GY63AG9seSLYfuYQvF+3C8Qs3vc51blIFj/aKQ4WSd1eklh3p0pC57Vl25Dq5PtLMJy9EhanI/g+bNm0aRo4cqXYt5M0z2jc7zI8lonAm3RSSpsyFOzFJc9zSqQ0iB/eBLgzyb/Pr5KWb+GrRbmw5dNHrXMNqpVSRWsOqhTtISNqSSZcGXwJfuU6uJ6IQCXrfeust1XtXujiMGzcO1atXh9ns/UtbdoJnzpzp0xALIqJQ40pKRvL3P6YNjshAV7wYokYPgalxPb+tLdDdSEhRRWqLtxz3KlKrWCoKj/dqgvaNKxXJpoqkKkhbsuXbT+V6rVzH1AaiEAp6z5w5A71ej7Vr16JOnTo5XjtixAjExMQU1dKIiAKC/cARJH47C+4b2pfjTa2awjpiIPRRVr+tLZCl2ByY/dNhTF97SBWsZVQs0oRHujTEgLa1C7RIzRfSh1daluVUzCZty+Q6IgqhoLd169Y4fvx4rgGvsFqtmDRpUpGsi4jI39w2G5LnLEHq6p81x3XWSFgffgDmVk39trZAJsViK3ecwoRle3HlZrLmnNGgw8C2tfFwl4aqO4M/SBsy6cqQW59etisjCrGg94MPPlBT1q5fv46SJXNv+C27wkREoc5x4gwSJ0yH65K2zsHYsC6iRg+GvkRxv60tkO08dlkNlzh6Pj7LdIHHe8WhYqmsC6aLkgyekD68nMhG5H9FNpFNTJw4EXv27MGHH36Y43XSwaFYsWJITU1FKONENqLw5XY6kbJwFVIWr5YtyzsnzCZYh/SFueO9LOjNwqnLt/C/Rbux8aB2Gp1oULWkKlJrVK00AnVnWro0SNEac3iJQngim6QryO6tpDj89a9/RZUqVbK8TgLdBQsWqII2IqJQ5Dx/CYkTZsB5+pzmuKFmVUSNGwZD2cAM2vzpxu0UTF6xHws3H1fBY0blS1jxWK8m6NSkckA/UZBAl23JiPynyP7vk0D36NGj6uN58+bleK1sPgfyLy4ionyPEV61Aclzl0qrmjsnDAZE9O+GiB4doTOwdVVGqXYn5m44gu9XH0BSqnYzJCpCitQaYOB9tWFmyy8iCpSg93e/+x2eeeYZFcyWKlUKkZGRWebtJicn4/Lly0W1LCKiIuG8diNtjPDh45rj+orl1BhhY5WKfltbIJLd3NW7TmP80r24HJ/kVQAm3RhGdmmAmCiL39ZIRMGlyHJ6JW2hWrVqWL16NRo0aJDjtTKuuHv37nA6fZtmE6yY00sUJmOEf9mGpOnzgZQMdQo6HSzdOyByQHfoTCZ/LjHg7D5+RU1SO3z2hte5do0q4fHecahcuphf1kZERSOoc3otFgteeeUVlC6de65aly5d0KNHjyJZFxFRYXHduo2kKXNg37lfc1xfqgSs44bBVKeG39YWiM5eScD/Fu/Ghv3nvc7Vq1xCFanF1Sjjl7URUfAr0u4NeREOeb3c6SUKXbad+5D03Ry4ExI1x83tW8E6tB90EXxZ3uNmYqoqUluw6ZhXP9tysVY82isOnZtUYccDojDSKJh3evMqKioKSUnaPC4iokDnTk5RqQyS0pCRLiYa1lGDYW6Sc3pXOLFJkdrPRzB1lXeRmtVixIj7G+DBdnVgNrFIjYgCLOi9cuUK1q9fj169eqlCtYxOnz7t021IIdv06dNDvkcvEYUe+6FjSPp2JlzXtAMTTC0aw/rIIOijo/y2tkB7JW/N7jP4ZvEeXMpUpCa7uf3urYlRXRshNpq74UQUoEFvp06dcOjQIQwaNAizZs3SnGvSpAkSEhIK8u6IiAKC225XbchSV67XHNdFRiDyoYEw39ss5NO1fLX35FU1Se3gmete59o2qIgn+jRBlTIsUiOiAA96XS6XegafVZrwiBEj8OWXX/p8W/wDQUTBwHH6HBLHT4frgrbVorF+bUSNGQJ9SY6ZFeeu3lZFauv3aQdyiDqVYvFUn6ZoWqusX9ZGROGhQIPetWvXYsOGDSq9IbOnnnoK06ZNw/fff48aNWqobg5ZBbaSxyvX/eUvfynIpRERFfwY4SVrkLJgpXaMsMmIyAd7w9K5LXRZ9CIPN7cSU/HdqgOYv/EoHE7thkiZ4pEY1zMOXZtVZZEaEYVW94ZHHnkEU6ZMyfU6WZLValX5vaGM3RuIgpPz4pW0McInz2iOG6pXThsjXJ47ljaHEz/+fFQVqd1OsWvOySjeh+6vj8Ht68LCIjUiCtbuDbt379Z8Xr9+fZjNZvXxf//7X59uQ3aAExO1bX6IiAJijPDajUievRiwZwjk9HpE9O2CiN73h/0YYdm0WLfnLL5ZsgcXrmt/j8tubp9WNTC6WyOUKBbhtzUSUXgq8KC3WbO0go0hQ4Zg6NChqFu3bvq5mJgYn28nqxHFRET+4rpxE4nfzoTjwFHNcX35MogaNxzG6pUR7vafuqaK1PafvuZ17t76FfBE7yaoVs73vwNERAWpUPr0/t///R9zcokodMYIb96J5O9/VD14M7J0bYfIB3pBZw7vMcIXrt3G/5bsUTu8mdWsUFxNUmtRu5xf1kZEVKhB78svv1wYN0tEVKRctxORNOUH2Lfv0RzXlSiucndN9WohnCUk2TB19QGVu2t3ZijmA1A6RorUGqNr82owsEiNiEIx6I2OjkaJEiXu+nY+//xzPPvsswWyJiKivLLvOYjESbPhvqXtL25u2wLW4QNUD95wZXe4MH/jMXy3ar8KfDOKMBswvFN9DOlQFxHmgB36SURhqMB/IxmNxgLp9yu7xQx6iaiouVNSkTRrIWw/bdYc10VHwTpyEMzNGyOcUz3W7z2H/y3ZjfPXMhWp6YBerWpgTPfGKMkiNSIKQAH5NHzy5Mmw2bS7B0REhc1x9KRqRea6qp0WZmraANaRD0IfE76TwmSCmhSpyUS1zFrVLa8mqdUoX9wvayMi8kvQe/PmTdSsWTPfuwjx8fG4detWQS+LiChbbrsDyfOWI3X5OvlFdOdEhAXWYf1hvq9l2E6JvHg9Ed8s3YM1u7Q9iYUEuU/2aYJ76pb3y9qIiPwa9ErgevLkybu+nXD9A0NERctx5jySZNDEuYua48a6NWAdMxSG0iURjm4n2/D96oOYu+GIV5GapC+M7dEYPVpWZ5EaEYVv0GswGNC+fft8B8w3btzA/v37VV4vEVGhDppYtk7t8MLpvHPCaETkAz1VO7JwHCPscLqwYOMxTF65H7cyF6mZDBjasZ56i7QEZHYcEVHRdm9YvXr1Xd3GDz/8gMGDBxfYmoiIMnJeuZY2RvjYKc1xQ5WKiHp0OAwVw6+nrGw6/Lz/PL5evBvnrt7WnJMX3nq2rI4xPRqrVmRERMGowIPegkhLGDhwYProYiKiAh008dNm1Z0BqRl2MXU6NUJYRgnrCqADTbA5dPY6vlq4G7tPXPE616JOOTzVpwlqVoj1y9qIiApKgf92T01NLZDAedKkSQWyHiIi4bp5S/Xddew9pDmuL1taDZow1qyKcHM5Pgnjl+zByp2nvc5VKxuDJ/s2UZ0ZWGNBRKGgwIPe5ORkXLt2DaVKlbqr2xk6dGiBrYmIwptt624kTf0B7sQkzXFL57aIfLA3dJbwemUpMcWOaWsOYvb6w2rQREYloi2q126ve6rDYAi/nGYiCl2F8jreq6++ivHjxxfGTRMR+cyVmISk73+EfcsuzXFdbAyixgyBqWFdhBOn04VFW05g0vJ9iE/UvipnMRnUFLVhnerBajH5bY1EREEV9E6cOBGbN2/GoEGD1GS12FjmghFR0bLvP4zEb2fBHa/t+21u3QyRDw2APsqKcMpl3nTwAr5atBtnrmjHKkvmQrfm1TCuZ2OUKR4+jwkRhZ8CD3rZaoyI/MmdakPy7EVIXbtRc1wXZYX14QdgvqcJwsnRczfUJLWdx72L1JrVKquK1GpXKuGXtRERFaXwK1MmopDlOH46bYzwZe2oXGPjeogaPRj64jEIF1duJmHC0r1YseOUZsicqFq2GJ7s0xSt67FIjYjCB4NeIgp6bocDKQtWImXJGu0YYYsZ1iF9Ye7QOmyCu6RUO6avOaSK1FLtGYZuAIiNsmB090bo06oGi9SIKOww6CWioOY8fwmJ46fDeea85rihVjXVisxQ5u46yQRTkdrirSfwrRSp3dYWqZmNejzYvi4e6lwfUREsUiOi8MSgl4iCZmww7A7AZFTjgdUY4ZXrkfzDMsDhuHOhwYDIAd1h6dExLMYIS5HalsMX1XCJU5e1RXuia/OqeLRnHMrGskiNiMIbg14iCmiOM+eRumI9bNv3ADY7YDbB1KguXNduwHk60+5upfKwPjocxsoVEA6OnY/HV4t2YfvRy17nmtQog6f6NkXdyixSIyISDHqJKGDZNu9UhWnI2BXGZod9xz7thTqd2tmN7N8dOlPo/1q7eisZE5fuxbLtJ72K1CqXjsYTvZugbcOKYZPHTETki9D/61DI/v73v+P111/HhAkTMHbsWH8vhyikdni9At4syKCJ6CcehrF2dYS65FQHZqw7hFnrDiElU5FajNWMUd0aod+9NWFkkRoRkRcGvXdh586deOutt/y9DKKQJCkNuQW8wli3ZsgHvE6XG8u2ncTEZXtxPSFFc85k1GNQuzp4+P4GLFIjIsoBg958SklJwciRI2G32/29FKKQI0VqKofXB/ad+9T1oVq0tlWK1BbtxomLN73O3d+0iipSK18y6q7vx+VyI9XhhMVogF7PtAgiCj0MevPp//7v/1CnTh3cvn0bp06d8vdyiEKLdGmQojVfyHVyvcWMUCJBrhSpbT18yetc4+qlVZFa/SolC6QYTnr6/rTnrEqZiDAZ0CGuMga3r4taFTlCnohCB4PefFi1ahWmT5+O3bt345577vH3cohCittmR9LcJb5/gdmk2piFCklfkDSGpVtPwJWpSK1iqWg83jsO7RtVKpAitVU7T+P9GZtV+oSHBL7Lt59S514b1hpdmlW96/shIgoEofOXoojEx8dj3Lhx+Oabb1C6dGl/L4copDhOnkXihOlwXbzi89eYW8aFRGpDss2B2T8dxvS1B5Fi0xapFZMita4N0e/eWiqHtyDIDm/mgDcjOS7nq5WN4Y4vEYUEBr159Nxzz6Ffv37o3bu3v5dCFDLcTidSFq1GyqJVPhWvpdPrYenaHsFMgssV209hwrI9uHYrU5GaQY+B99VWRWoS+BYkSWnILuDNuDa5TnZ8iYiCHYPePJgxYwa2bduG7du3+3spRCHDefGyak3mPHlWc9xQo4raxU2esyTrQFivV2OGjVUqIlhtP3oJXy7cheMXvIvUOjWpjMd6xqFCqegCv18pWpMcXl/Ida8MacXiNiIKegx6fXT+/Hm88MILWLBgAazWghvneezYMTRq1Cjb8/v2ZWrCTxQi1BjhNb8gec7itEI0D70eEf26IqJXZ+gMBhjr11bjhm3b7kxkk2BYdniDNeA9dekWvl68G5sOXvA617BqKVWk1rBaqUK7f+nSkLnPb3bkOrk+0sw/F0RUsHKKfyQ+qlWrVoHeH3+L+Tjb/tFHH8VvfvMbFq4RFQDX9XgkfjsLjoNHNcf1FcoiSsYIV62UfkwCW+PYYbCOHpIWHJuMQZvDeyMhBd+u2IfFm497FalVKBmFx3s3QYfGBVOklhNpSyZdGnwJfOU6uZ6IKNgx6PXBZ599hoSEBNWmrKDJsxju5lI4PYG0bdqB5Gnz4E5O0Y4R7toOkQ/0hM6U9YAFFegGaVuyFClSW38E09ccVAVrGUVHmPBI14YY0LYWzEUUXEqqgrQlky4NuZHrmNpARIUhp/gnp13g/GLQm4tDhw7h3XffxS+//AKDgbsdRPnlSriNpClzYd+h/SWnLxUL65ihMNUr2JexAoHkzq7ceRoTlu7BlZvJmnNGgw4D2tRWAa+MEC5q0odX2pLlVMxm0OvUdUREoUDnlq0XytYTTzyBCRMmICIiIsvzSUlJavfKYrHAaDSiWrVqPu/cep7FcKeXQp1t134kfTcH7lu3NcfN97WEdVh/6CKz/v8rmO06dhlfLtqFI+fivc5JCsNjvZqgUumCL1K72z69GQNe9uklIn8pjBiJO725kDHDTqcTiYmJOV6Xmpqq3nK7jiicuFNSkTRjPmwbtmqO64pFwTryQZibFfzLV/52+nJakdrGA95FajJBTYrUZKJaIJCAVvrwciIbEYUD7vTeperVq6sxxLIbPHbs2Dx9LXd6KZTZDx9H0sSZcF27oTluatYQ1kcehD7Gv7ucBS3+diomrdiHhVKklmnntHwJq9rZlTZkhV2kll+yZunSIEVrzOElIn/jTi8RBTy33Y7kH5cjdcVPUrl250SEBdaHBsDcpkXABn75YbM7MWfDEXy/+gCSUrVFalERJjVY4oH7asNsCuyaAAl02ZaMiEIZf8MRUYFxnDmPxPHT4Tp/SXPcWK+mKlYzlCqBUCE7o2t2n8E3S/bgcnySVz5s/za1MLJrQxSPsvhtjUREdAeDXiIqmDHCS9chZcEKwJmh96vRiMhBvWDpcl/Q9tbNyp4TV9QktUNntakbol3DiqrfbuUyxfyyNiIiyhqDXiK6K85LV5E4cQacx09rjhuqVlJjgg0VyyFUnL2SgP8t2YMN+855natbuQSe7NMUTWuW8cvaiIgoZwx679LJkyf9vQQi/w2aWLcJSbMWpo0HzjhGuPf9iOjbRY0RDgU3E1Px3cr9mL/xmFd7r7KxUqQWh85NqrAAjIgogDHoJaI8c924icRJs+HYf1hzXF+uNKLGDYexRhWESpHaj78cxZRVB5CYkiGwB2C1GDHi/gYY1K4OLAFepEZERAx6iSiPbFt2IWnqD3AnaSeMWe6/D5EP9oLOHJyjgjPvYq/dfRbfLNmNize0RWqym9u3dU2M7tYIsdEsUiMiChYMeonIJ67EJBXs2rfu1hzXlSiOqDFDYGpQB6Fg78mrqkjt4JnrXufaNKiAJ3o3QdWyMX5ZGxER5R+DXiLKlX3vISROmgX3zQTNcfO9zRH50ADorZEIduev3cb/Fu/GT3u9i9TqVIpVRWrNapX1y9qIiOjuMeglomy5U22qUE0K1jLSRVlhfWQQzC3jEOxuJdkwZeV+zNt4FA6ntkitTPFIjOsZh67NqrJIjYgoyDHoJaIsOY6dQuKEGXBduaY5boqrD+uoB6EvHtwv8dscTsz/5RimrNqPhGRtkZpMJhveuT4Gt6+DCE4pIyIKCfxtTkQabocDKfNXIGXpWu0YYYsZ1mH9YG7XKqjHCEuRmqQwSCrDheuJmnOymdu7dU2M6dYIJYpF+G2NRERU8Bj0ElE657mLaoyw8+wFzXFj7eqwjh0KQ5lSCGb7T19TRWr7T2l3r0XreuVV3m61csG9g01ERFlj0EtEcLtcSF3+E5LnLQMcGccIGxA5oAcs3TsE9Rhh2dGV9mPShiyzmhWK46k+TdGiTuhMjiMiIm8MeonCnPPqdSRNmAHHUe10QUPlCoh6dDgMlcojWCUk2fD9mgP4YcNR2J0uzblSMREY26MxureoDgOL1IiIQh6DXqJwHiO8YQuSZiwAUm13Tuh0iOjZCRH9ukFnCs5fEXaHCws2HcPklftV4JtRhNmAYR3rY0jHuqpgjYiIwgN/4xOFIdfNBCRNng37noOa4/oypRA1dqjK4Q3WQH7DvvP4evFu1Xc3I9nM7dWqhpqkViom+PsKExFR3jDoJQoztu17kPTdXLgTteN1zR3vhXVwH+gignO07qEz11WR2p6TV73O3VO3nCpSq1G+uF/WRkRE/seglyhMuJKSkTxtHmybdmiO64oXQ9ToITA1rodgdOlGIsYv3YtVO097nZMg94k+TdCqbvDmJRMRUcFg0EsUBuwHjiDx21lw37ipOW5q2QTWhwdCHx2FYJOYYsfU1Qcwd8MRlcObUcliaUVqPVqySI2IiNIw6CUKYW6bDclzliB19c+a4zprBKwPD4K5VVMEG4fThYWbj2Pyin24mZipSM1kwNCO9dRbpIW/3oiI6A7+VSAKUY6TZ5A4fgZcl65ojhsb1lHpDPoSxYOuSO2XAxfUJLUzVxI052RAXI8W1dXubuniLFIjIiJvDHqJQozb6UTKwlVIWbwacGV42d9kQuSQPrB0ahN0Y4QPn72Brxbtwq7j2gBetKhdVhWp1aoY65e1ERFRcGDQSxRCnBcup40RPn1Oc9xQoyqixg2DoVxpBJPL8UkYv3QPVu7wLlKrVjYGT/ZNK1ILtiCeiIiKHoNeolAZI7zqZyTPXQI4HHdO6PWI6N9NDZvQGQwIpiK16WsPYvZPh2HLVKQWG23BmO6N0PueGjAYgnc0MhERFS0GvURBznU9HokTZ8Bx6LjmuL5iObW7a6xaCcHC6XRh0ZYTmLRiH+Jvp2rOmY16DOlQD8M714PVYvLbGomIKDgx6CUK5jHCG7cjado8ICVDgKjTwdKtAyIHdofOZAqa72XzoYsqb/f0Ze8itW7Nq2Fcz8YoU9zqtzUSEVFwY9BLFIRct24jacoc2Hfu1xzXl4qFdewwmOrWRLA4ej5eTVLbeeyy17lmNcvgyb5NUadSCb+sjYiIQgeDXqIgY9u1H0mT58CdcFtz3NzuHliH9oMuMgLB4OrNZExYthfLt5+E2609V6VMMTzZpwnurV+BRWpERFQgGPQSBQl3cgqSZi6AbcNWzXFdsWhYRz0Ic9OGCAZJqXbMWHsIs346jFS7U3MuNsqC0VKk1qoGjCxSIyKiAsSglygI2A8fR9LEGXBdi9ccNzVrBOvIQdAXi0YwFKkt3XYSE5ftxY1MRWomox6D29fFQ53rIyoiOPKQiYgouDDoJQpgbrsdyT8sRerKDVLtdedEhAXWhwbA3KZFULz8v+XXIrWTl255nevavCoe7RmHsrE5F6m5XG6kOpywGA3Q6wP/eyYiosDCoJcoQDlOn1ODJlwXtAVexnq1EDV2KPQlA38C2fEL8fhq0W5sO3LJ61xcjdJ4qm9T1KtcMsfbOHY+HrPXH8ZPe84ixe5EhMmADnGV1c4wp7AREZGvGPQSBeIY4SVrkLJgZaYxwkZEPtgbls5todMHdr7r1VvJ+HbZXpXOkLlIrXLpaDzRuwnaNqyY6y71qp2n8f6MzXC67tyIBL7Lt59S514b1hpdmlUtrG+DiIhCCINeogDivHQFiRNmwHnijOa4oWolRD06HIYKZRHIkm0OzFp3SBWqSXCaUYzVjFHdGqHfvTV9KlKTHd7MAW9GclzOyzhi7vgSEVFuGPQSBchwhtQ1G5E8exFgt2vHCPe5HxF9ugT0GGEJQJdvO6lakF1PSNGcMxn0GNSuDkbcXx/RkWafb1NSGrILeDPer1wnO75EREQ5YdBL5GeuGzeR+O0sOA4c0RzXlyuDqEeHwVi9CgKZ5Ot+tXAXjl+86XWuc9MqeKxnHMqXjMrTbUrRmuTw+kKue2VIKxa3ERFRjhj0Evlxd9e+ZReSvv8B7iTt7qily32IHNQbOnPgtu86cfEmvl60G1sOX/Q616haKVWk1qBqqXzdtnRpyJwekR25Tq6PNPPXGRERZY9/JYj8wHU7EUlTf4B92x7NcV2J4ogaMxSmBrURqCR94dvle7Fkywlkzj6oWCoaj/eOQ/tGle6qlZq0JZMuDb4EvnKdXE9ERJQTBr1ERcy+5yASJ82G+1aC5rj03I0c3h96ayQCUYrNgdk/Hcb0tYdUwVpGxSJNGNm1Efq3qaUGTdwtSVWQtmTSpSE3ch1TG4iIKDcMeomKiDslFUmzFsL202bNcV2UVU1VM7eIQyCS/NoVO05hwtK9qhVZRkaDDgPb1sEjXRqgmNX3IjVfSB9eaUuWUzGbQa9T1xEREeWGQS9REXAcPalakbmuXtccN8XVh3XUYOiLF0Mg2nH0Mr5cuAvHLmjHH4tOTSqrIrUKpQpnBLK0IZOuDNm1LZOAV86zXRkREfmCQS9RAXDLEAm7Qw2QyDg4wm13IGXBCqQsXasdI2wxwzqsP8zt7gnIMcKnLt/C/xbtxsaDF7zONayaVqTWsFr+itTyQgZPSB9eTmQjIqK7pXNLCTn5RaNGjdT7ffv2+XsplE+OM+eRumI9bNv3ADY7YDapNAVLt/aAToekCdPhPKvtbmCsXR3WccNgKJ3z+F1/uHE7BZNX7MfCzcdVWkNG0nbs8V5x6BhX2S+BuqxHujRI0RpzeImIQlujQoiRuNNLlE+2zTtVyoJmVLDNDtvG7bBt2pH2ecbnlEYDIgf2VAFxoI0RTrU7MWf9YUxbcxBJqdoitegIEx7p2hAD2taC2Y9dEiTQZVsyIiLKL/4FIcrnDq9XwJtRphdQDFUqIGrccBgqlUcgkd3TVbtOY/ySPbhy07tIbUCb2qpILSbK4rc1EhERFQQGvUT5ICkN2Qa8mUT0vh8R/bpCZwys/912Hb+iitSOnLvhdU767D7euwkqlS6cIjUiIqKiFlh/hYmCpGhN5fD6wmRExIDuAZXOcOZKgpqk9suB817n6lcpqYrUGlcv7Ze1ERERFRYGvUR5JV0apGjN12vlzVKwPWzzI/52Kiav3IcFm7yL1MrFWvFY7zh0iqvCIjEiIgpJDHrzYO/evfjLX/6CVatW4ebNm6hUqRL69u2L119/HRUqVPD38qiomIyqS4NPga9cJ9f7kU2K1DYcwferD3gVqUVFmPDw/Q3wwH21YTZxlC8REYUuBr0+WrJkCQYNGoSUlJT0YydOnMB//vMfTJ8+HUuXLkXz5s39ukYqGo6DxwAfG/2ZW8b5LbVBdnPX7D6jitQuxSd5DXbod28tjOrWEMVZpEZERGGAQa8Prl27hhEjRqBNmzb43e9+h/r16+P8+fP44IMPsHDhQly5cgUDBw7EwYMHYbVa/b1cKiTuVBuS5yxG6ppffPsCvR6Wru3hD3tOpBWpHTrrXaR2X8OKqkitSpnAnAJHRERUGBj0+uDrr7/GAw88gPHjx6c35a9Tpw46duyIUaNGYcqUKThz5gy+//57PPbYY/5eLhUCx4nTSBw/A67LV337Ar0eUeOGwVilIorS2asJ+N/iPdiw75zXuTqVSqgitaY1yxTpmoiIiAIBg14fbN26FZMmTfKaQiWff/zxxyq9weFwYPv27Qx6Q4zb4UDKwlVIWbxa23vXbIJ1SF8YalZF6sr1sG3LMJGtZZza4S3KgPdWYiq+W3UA8345CmemIrUyxSPxWK843N+0KovUiIgobDHo9cF7772XbdpCmTJl1Ki8Xbt2ISIiosjXRoXHef6SGkDhPK3dNZVAV3ZxDWXT2noZxw6DdfSQtC4NJmOR5vDaHE78+PNRTF11ALdTtIV1VosRI+5vgEHt6sDCIjUiIgpzDHp9ULt27RzPm81p7agaN25cRCuiwu7Dm7pyA5J/WAo4MnQ7MBgQ0b8bInp0hM6gDSJVoFuEbcncbjfW7TmL/y3Zg4vXEzXnZDe3b+uaqkitRDSfiBEREQkGvQUQfBw/fhwWi0Xl/VJwc169jqRvZ8Jx+ITmuL5iOUQ9OrzIc3Szsu/UVVWkduD0da9zbRpUUEVq1crG+GVtREREgYpB71365ZdfVHeHF198ESVKlMjz1x87dkylR2Rn3759d7lC8vXJi+3nbUiaMR9ISb1zQqeDpXsHRMpUNZPJn0vE+Wu38c2SPWqHN7PaFWPxZJ+maF67rF/WRkRElFc5xT8SH9WqVQsFiUHvXfrkk0/UYIo333zT30uhfHLdSkDSd3Ng33VAc1xfuiSsY4fCVKcG/OlWkg1TV+3Hj78chcOpLVIrHROJR3s2Rtfm1VikRkRElAMGvXdhw4YNmDNnDhYtWoTY2Nh83YY8i+FubtHm62YsOLPt3KcCXneCNi/W3L41rEP7Qhfhv8ENdodLdWOYsmo/EpK1RWqRZiOGd66Pwe3rIMLM/42JiCj47Msh/slpFzi/+NcynxITE/H444/jnXfeQY8ePfy9HMqF48x5pK5YD9v2O63F9MWLwXVFmxeri4mGddRgmJs08Guqxfq95/C/Jbtx/lqmIjUd0Lt1TYzu1ggli7FIjYiIyFcMevNJAt527drh9ddf9/dSKBe2zTtV6zHILm/6QbtXwGtq0RjWRwZBHx0Ffzlw+poqUtt36prXuVb1yuPJPk1QvVxxv6yNiIgomDHozYc33ngDKSkp+O677/y9FPJhh9cr4M1CxMAeiOh9v9cAkqJy4Xoixi/ZgzW7z3idq1m+OJ7s2xQt65Tzy9qIiIhCAYPePPryyy+xfv16lcdryNSrlQKPpDTkFvAKGS/sj4D3drINU1cfwA8bjsLu1K5T0hfG9WyM7i2qw8AiNSIiorvCoDcPJkyYoN6WL1+e5fQ1GUX8ww8/YMiQIX5ZH3kXrakcXh/IGGGZqlZU09QcThcWbDyGySv3q+4MGUWYDRjWsT6GdKyrCtaIiIjo7vEvqo++/fZbvP/++5g3bx5SU1PVm6foSFIdDh48iH/84x8YN26cv5dKv3KevZhWtOYLuU66OhTyVDX59/Lz/vP4evFunLt6W3NONnN73lMDY7o3QqmYyEJdBxERUbhh0OuD//73v3juuedUwFK3bt1sr4uOjlZBMQXAGOE1vyB59iLfv8hsUm3MCtOhM9fx5aJd2HPiqtc5ydeVIrWaFfLX+o6IiIhyxqA3F3PnzsWzzz7r07WDBg2C1Wot9DVR9lzX45H47Sw4Dh7N09eZW8YVWmrDpRuJGL90L1btPO11rnq5GFWk1qpu+UK5byIiIkrDoNeHQFZ2eCkIxghv2onkaT/CnZyiPSkFajn9DPV6WLq2L/A1JabYMW3NQcxef1gNmshcpCZpDJLOwCI1IiKiwsegl4Ke63YikqbMhX37Xs1xfclYNUbYfTMh+7Zlej2ixg2DsUrFu1uDy41UhxMWowEutxuLNh/HpBX7cDNRW6RmMRkwtGM9DOtYD5EW/u9HRERUVPhXl4KabfcBJE2eDfctbVGYuW1LWIf3hy4yrcuGvkJZpK5cr7o0eCaySUqD7PDeTcB77Hy82sn9ac9ZpNidMBn1MBn0SEp1eG0292hRHWN7NEbp4ixSIyIiKmoMeikouVNSkTRzAWzrt2iO66KjYB05CObmjTXHJbA1jh2m2pKpLg0m413n8EqO7vszNsPpupM6IWkMmVMZmtcuiyf7NEXtiixSIyIi8hcGvRR07EdOIGniTLiuZhoj3LQBrCMfhD6mWLZfqwLdAmhLJju8mQPerDw3oBkGtq3tt0lvRERElIZBLwUNt92O5HnLkbr8J21hWoQF1mH9Yb6vZZEFl9PXHsw14BWHz95gwEtERBQAGPRSUHCcOY/E8dPhOn9Jc9xYtwasY4bCULpkkazD6XRh4ebjWL3rjE/XS67vK0NaQc8ODURERH7FoJcCdsCE5N669TqkrliPlPkrJOK8c4HRiMgHesLStV2RjA6WlmibD13E14t249TlWz5/nRS3SVcHjhMmIiLyL/4lpoDb0ZUg17b91y4LWfTYNVSpiKhHh8NQsVyRrEnyd2WS2o6jl/P8tREmg2pjRkREFApcbhdS7KmIMFmg1xX+plNBYtBLAcO2ead3P92MAa9Oh4g+9yOiTxfojIX/T/fqzWRMXLYXy7af9JptYbUYvdqSZaVDXGWmNhARUdA7cPkwvtkyBUsOrUSyPQWRpgj0qtcVj7V6BA3K1kUwYNBLgZOzm90ACQ+dDqbmje8q4M04RCK7YDQ51YEZ6w5h1rpDKj0ho+JRZozu1gj1q5TEC5+vyrGYTSatDW4fHL8IiIiIsjNv/xK8uvBNOFx3/iZK4Dt370LM378EH/R9GwMa9kKgY9BLAUFSGnIMeIXLpQZMSL/dux0iIWkHsgsrQWmtX/vnSgC7dOsJfLt8H64naEcZy9AJufahzvURFWFSx14b1jrbtmUS8Mp5z20TEREFGpvTjpvJNxGfcgvxyTdxQz5Ojtd8fubmOWw8tQXZbfFIICwBcZ3SNQN+x5dBL/mdMyFRpTb4QiaqyYCJvBSvZTVEQgLf5dtPqXMSnBazmvHVwl04ecm7SK1rs6oY17MxypWI0hzv0qwqqpWNyTWYJiKiohHM+aZ3+30npN7+NWj99S3FE8TexM0MH99Qn6cFtbdtiQVy/xL4jt8yFR/0fQuBTOeWsnTyi0aNGqn3+/btQ7iy7z2E29/OBDKNEc5J7KfvQOfjgAnZ4X3uPytyTEOQJIeszsZVL42n+jZFvSolCyRtgoiICkco5JsKCclSHKm4IbutmgA1/k5Am3FXNiUeN5NvqWMS+PpTpCkCu19aV2BPNgojRuJOL/mFO9WG5NmLkLp2Y96+0GxSI4R9JbuwuQ2RyHy2UuloPNG7Ce5rWNHnwRIS6LItGRFR0QvUfFO705G+o5q265oWyMarINU7lcATyNqcNvhThNGC2MhYlIgsjmKWaGw+s92nr5PHXHbZreZIBCr+laYi5zh2CokTZ8B1+Vqev9bcMs7n1AbZfZW0A18VizSpIrV+bWrBaAifl8WIiIJ5hzdzwFvQ+aay++pJHfBOE7jzuXYXNr7AUgfyy6AzIDayOGIjYtT7EvJxZHEUj7jzcYlMn8u1EaaI9NuQ3eMmH3dUAa0vO72SVhLIGPRSkXE7HEhZsBIpS9ZoW5FZzIjo1h4pi9fkXMym18PStb3P9yfpBpm7L+Tkqxd7oHRxq8/XExGRf0lKQ3YBb1b5pin2lF/TAn7dcZW0AfWx5/O0lIGMqQSSPuB0+/63pDDIjqsnQM0YwN4JaGMRq87d+TzaHOXzq5XZkVQFSRORXfPc9K7XLeDzqBn0UpFwnruIxAnT4TxzQXPcUKsaosYNg6FMKRjKl82+bZler64zVqno831Kfq0UlvkS+Mp1JYsF7ksyRER0h8PlUIHpooMrfLp+7t4FWHhgGVL9nDpgMVpQ4tfAtXhkTPrHKmjNEMR6jhVXn8fAqPdfuPZYq0dUmkhOTy6MegMebfUwAh2DXir0ccLSjiz5x6WAI+MYYQMiB/SApXuH9HQFc+tm0Fcoq9qSSZcGNZHNbFIpDbLDm5eAV5y+cgtRkSafgl4OkSAiKnqSOiBpAHfaZXl2YX9NIci4C5typ7hL0g3ydD/y6l8BBryyo6kNUDOlCagANsOu7K8fSwpAsGlQtq7Ki84ujUQCXjkfDAWDDHqp0DivXkfShBlwHD2pOW6oXB7WccNhrFzB62sksJU+vNKWDHaHKlrLS3sycSMhBd+u2IfFm48jlxq2tPVwiAQR0V1LVV0H0roNZO4+kF60lemY5MTmlp5Q2CQNID14/TU4VWkCEbEZAtqYXwPatAKvaEtUwL+UX5AGNOyl8qIlTWTxoRXpHTIkpUF2eIMh4BVsWeZHodqyTP5J2TZsRdKM+UCqTTtGuGcnRPTrBl0uHRjy0wIsxebA7PVHMH3NQSTbch8RnHGIhPTcJSIqTMHSQ9bpcuJmSoLXkAJP/qtqkaXZhU3bnfWl2MlfGpdrgJEthmjSBjy7syYD9/8C8d8xW5ZRwKUuZN6Ndd1MQNLk2bDvOai5Vl+6ZFpObu3qdz05LasAeeWOU5iwbC+u3EzWnDMadBjYto5qP7Zk6wkOkSCisOkhKxsQibYkzZCCjOkDWXcjuIVbKQlwZzt/q/BJIFU8otivKQKxd4qz0ndh76QLSCD+0oI/q0A9O/Ly+3t9/hw0u5GBTq/TB3Rbspxwp9ePgnWn13HmvMrTtW3PkHfbIg76yuWRungN3IlJmuvNHe+FdXAf6CIseZ6cltuO7M5jl/Hlwl04ej7e62s6xlXGY73iULFUdPoxDpEgIn/3kM2cC+lLD9lUh02TNnCnPZb30ALPLqyMl7W7fHvVq7BEma2a3Na0HdZfi7Yydx/49RrpVJCXHcSCeowp9GMkBr1+FIxBr4wLzrbDQia64sUQNXoITI3rFcjkNAl8P/tNN7Uze/ryLXy9eDc2HtB2gxANqpZUk9QaVSvtw3dERFR4O7wPfDsqx5xVg06PF9o/CaspMtPQAu0ubJJd+ypWUTPpjZquAl6dBlRg++uu7K8FXcUji8NsMBXZYx3s+aakxaA3xARb0Cs7vAl/+49PAa+pZRNYHx4IfXSUT7ctO7zLt5/K9bpOTSojxmrBQilSyxQgly8Zhcd7xakd3rvtTUhElBP50ymBaOZ0gYxpAyuOrMOZm+cQSHTQqd3WjO2ysuo+kPlzCcqD4fdqsORNU+6Y00t+JSkNvgS8xjrVEfXECJ9/QeZlctra3d7XRUeY8HCXBhh4X22YjQafboeIyMPmtKtUgByHFqQHtvHpo2Xl6/xJAtGMAaq2PVbWQwtiLMVg0Ifu78lgzjelwsegl3wuWlM5vD5wnDqXNnHNx6A3r5PTMqY7DGhbGyO7NEBMVGCPPiQKdKGwQybfg2dcbMYANbsiLs/nUuzlb6WsJVHy17QB7dAC7/6vnu4DFqPZ38smCioMesk30qVBitZ8IdfJ9RZzgU9O82jfqBIe6x2HyqWL+fw1RBQ4nQVySx2QtWQOUDWdBrIYWiBttiTw9SfZSfXku+69dNCn9chjvvE3S4L2yQZRsGDQSz5xXr2RtnPrSwq42aTamPlKuihI+zBfcnpjrGa8Neo+xNUo4/PtE5HvVe8SbM7du1CNHS2Iqne706GC1TtFWtkPLcjYE9bm53GxEoh6ug5kbJd1Zxc21mtogbTZypg68MrCN9VjmRspuGLAS1T4GPRS7mOEV/+M5LlLfAt4JeZtGZfnKWrSL1daluXUvUGv0+G9xzqiTqUSebptIsp6hze7Nk9Cjst5mcIkO76y+3ondSDz0ALv0bEywEACXBkx608GnSGLqVq/dhpIb5flPTbWYrz7lCnZLZcnDzl1b5CWWtJhgIgKH4NeypbrejwSJ86E49Ax379Ir4ela/s83Y/N7sS2I5dgNOizbTDu6dPLgJfo7qXYU/CfDd/kOv5Vzj805QnVdkpSB5xu/46Llf6tngBV02kgi6EFnnMyYtZfXQfkyYLslufWQ5YttYiKBoNeynqM8MbtSJo2D0hJvXNCp8PNKlUQdfp0lv9wHNDhbJeuaF6los/3s2b3GYxfsgcXb2RdSGIxGVQLMk5OI/LmcDl+HRd7J21AM6Agi0Iu2YVNdWT4/zoXhbFTK7uo2kBVO6Agq2OSOmDUB9+fLEkPkd1y9pAl8j/26fWjQOzT60q4jaTv5sK+U7smfalYxPfrg6cWHEZlWyL6pl5GG9sNRMCFFOix0VwCCy1lcdYclT5AIid7T15Vk9QOnrnuda5tg4qqSK1srJWT0yikugtkR34NS3CZOV1ABa2SJpBpbKznY0k38HfqgOS4atIEMoyN9ezCFs/QSkveR5giEI5C+d8wUUFjn17Kd16u6qZgMuaYa2vbtR9Jk+fAnaD9Q2pudw+sQ/vhy/m7Vc7tKaMVnxur47/WajDDBRv0cHtePnS5MXv9YZWKkJVzV2/jmyW78dNe74btdSrF4qk+TdG0Vtm7/ZYphARid4GcyC5q5jxX7ZhY70IuCWRzSzXwF9ld/UPn36KkNeMUrrQirmhLFIO3PGAPWSL/YtAbwmSCmgyUUP11pY2Y2QRzizhYurWHMUMKgjs5BUkzF8C2Yavm63XFomEd9SDMTRtmOUBCAt1UeDc5l+teGdJKs0N7K8mGKSv3Y97Go3A4tS8ulCkeiUd7xaFL06rc1aUi7y6QHckvT0sdiPfqNOAp1MrYOssTvMr6/MlsMGcozPLsuKYFqRmHFEzdMRvrTvyS6+3J4zuu1YgiWTsRUWFi0BuibJt3InHCDO0ENZtd5erKuahxw2Bu3Qz2w8eRNHEmXNduaL7e1KwRrCMHQV8sOs8DJOQ6uT7SbITN4cS8X46pgPd2irbPr9VixEOd6+PB9nVV7i7R3XQXyCl1QIYPpKULeE/Vym5X9lZqAvy9Kyh5rGnB66/tsbIYWpC5+4DshPtSuFW5eEX8fGozOwsQUdhg0BuiO7xeAW9GLpc6b9t3CPZNO7WtyCIssD40AOY2LTR/OPMyQEKuMxv0WLv7DL5ZsgcXrmsLYWQ3t2/rmhjVrSFKRIdnbh/lTlIafOku8Oel76F7nY6aoQUZd2FlvKzd5YA/SQeB4tm0x9K204pN/1w6FRRm6gA7CxBRuGHQG4IkpSHbgNfD5YJ94w7NIWO9mrCOGQpDqRJ3NUCiSc0y+N2Xa7D/9DWvc23qV8DjvZugWrkYX74VClGSOiA7qXcmamXoPvBrf9d5+xb7dFs7zu9Wb0VBWnfFZgpQM6YMeIYWqM9/LdqSIi75ukDEzgJEFE7YvSHEKhOlaC3+t2/6PjJYmIyIHNQLlvvvy7HQ7dj5eDz3nxU5DpCQveGsztaqEIsn+zZBi9rlfF8XBTz59ZFkT9bktN5JIch+aIGkFriz/JdSNHTQpe+83sl7vZMykN7/NdOurNUU6beer4WNnQWIKJCwewPlTro05CHg1VeugOjHR8BQIfeOCdKGTLoyvD9jc7aBb+ajpWMiMa5nY3RrXo1FagHO5rSrVAAVrHrGxv5atOUJWu9M4UrblZXgVb7O3yoUK6e6C2hyXNN3Wn/dlU1PJ4hBjEU7LpbYWYCIQh+D3lBjMsJtMkFnzz0Qcev1KPaHZ6E3+f7Sa5dmVVGtbIxqS7Zuz1mkZpPjG2E24KFO9TG4Q11EmEP/n1kg7ZLJWm5J14FMAWqOQwtSbqpiL3+SHFLVZeDX4PTszfO4kHAp1697sHE/fND3rSJZIxERBa/Qj0bCjKQnHClVAXUvns712iPlKuPePAS8HjUrFMd9DSupARNeRWo6oHermhjdvRFKFosI+eC0MHvISuqA3KZ3r1dtusCdvFjP7muC+j79SVIH0vJes04byFzEJZ9Hma2a1AF5bB/4dhS7CxARUYFg0BtipJ/ut6kxeDuXH67Usn+bEoNWLnee0g5kgtpXC3dhz8mrXuda1SuPJ3o3QY3yxREOwWleesjanTIuNuNAAu8BBVkNLbA5bfAneZw0AeqvH99pnfVr/9cMQwukzVZBpA6wuwARERUkFrKFWJJ2ss2BAW/MRTvbdTyXeCLLwFcC3s+iamCDuSTmvTNI9dPNzcXriRi/dA9W7zrjdU6C3Kf6NkXLOuWCJjjNHDj5OuBAgvHbqYnYenYnnp7zCpzunFtqlY0ugyRbkhox60/yfWp3WrVDC7S7srHpn1uMFvib/JtgdwEiovDSqBBiJAa9IfYDlZ3egW/OVf10qzmS0Df1MtrarsMCN1Kgx0ZzCSy0lFWjhKWf7o9vD8pxp/d2sg3frzmIuRuOwO7QvmReopgFj3Stgz6t6sJkuPudvYIMTvPyErlBZ8DbPf6gdig1qQNZpBNI6kBugW5hk/6td4LXX9tjZdiFzZg+4LlO+sQGe9eBQMqbJiKiwsWgNwAcPnwYb775JtasWQOHw4EuXbrgr3/9K2rXrh0QP1Ah3RUy9tPVud0wwwUb9Gp0sEf3FtVUN4asOJwuLNh0HJNX7FMjhDUsV1Gs8kEcvrUNyY6C2Y31NX/zhzGTve7D4XKonNaMaQOewq3Ze+bj8NVjCEQRRkuOAwqyOiaBuVHPrCQiIgptjdiyzL8WL16MIUOGYNiwYTh48CCMRiNefPFFtGjRAkuXLkXbtm0RCAa3r4tVO0+ntxVz6YBkOKGDXvXRFQa9Tl2XmTwH+mX/eXy9eDfOXr2tOSfxcs26V7Hq6kQ4r+eex1oY07eemfMKqpespvJjPbuxCanadQaaP3d9GeWiy3gNLYgwBV6hHxERUajiTq+Pjh49imbNmqkd3W3btsHw68v5stsrz0auX7+OAwcOoHTp0n7f6RUS9L47az6um7cg0XwIbp0dOrcJUbZ6KGlrhT8P6a/aj2V06KwUqe3G7hNXvG6vRZ1y6NYmCr9Z9KxPu7E1S1bzqWhLdmNvJMXj2PWTCASSBpC5aCstUL0zHjYmohiem/saUn0oMpNd8N0vrePL8URERHnAnV4/eumll5CYmIjnnnsuPeAVstv75JNP4pVXXsFrr72G8ePHIxDcNu/HuZhJmvxTCXxvW/YiOeIAbptrAUgLei/HJ6kitZU7vNucVSsbjRHdaqJGJQveXfmRT7uxAyeOhNPPLbNk7KukBVxNvO7T5C9JNdj+29WwGM0+3X6fBt3V7nZupOCKAS8REZH/cafXB8ePH0etWrXSP65Ro4bmvDwLady4McxmM86fP49SpUr5dafXt+ItPZ5r+xQOnkrA9uOnYHP/f3v3ASVFlf1x/A055yxJFCQJAiJRERQBEREQSSthYcGwYljETFKCqAcMgICCYRXwsBIkCUhaQGUJKoqAIAooOYnkUP/ze/6rz4TumR7oobprvp9zmp5O1a+7mpnbr+6796S5GHfKXMhw2p6bjKdNhsxnzanz3lYdSEwBpNInCqjCQJBOW+5lzbBq4Va/uQPDCk5T2+DgcnKQAQBA8pjp9ci8efPsec6cOZMEvFKxYkWTNWtWc+bMGTNjxgzTq1cv46Vw8mM1E/vG6nF/XQhVlUq1zSJIzQfiVxWIv0Br+c+rzca9m1Lcxj1V7kxVcKrFdco1jnSDA2rIAgAQWwh6w7Bw4UJ7XqpUqaC3K92hZMmSZvv27WbNmjWeBr0q66Qat17KnDGzGd1qqCmQ469ZVxvgZs9rUw5CaX5dkzTpvpWWwalmnMsXKkcNWQAAYgBBbxiU0iAKbENRSoOC3i1bthgvqY6pgq9wxV3ManJmym1KFShsiuYpkKC2a+Jar2NWTzKfb12S4jZbVWpmg9j0EJz+Ne5B5uU7B1BDFgCAKEbQG4YDB/6qZpA7d+6Q91F6gxw5ciRV21ag7OatBJPaXBYFXQrowgl8M5jM5oN7PjF1K5YIq3HBIw16mS+2LY/4bKwfglNtK0eW7BHbHgAAflclmfhH8ZG7nipSCHrDoHJkkiNHjpD3uXjxr2oFp0+HP8uaFhR8qVFEOIu3Wle5w9SrdFVU5bESnAIAgLRA0BsGVWVQPd7knDt3zp4XKFAgVdvWt5hIV28Id/FWz5u6RG0eK8EpAAD+9kMy8U9ys8CXiqA3DMrXPXnypDl16lTI+xw7dsyep6Y5RVpJ6xlZ8lgBAECsIegNQ6VKlcyuXbvMvn37Qt7n0KFD9jxYSTMvJJiR3bzYnDof+RlZZmMBAECsIOgNQ8OGDW3Zsh07dgS9XZ3aDh48aH9u2rSpiRbMyAIAAPyFCCgM7dq1s+fqtrZ3794kt2/cuDFQwaFJk9SV6roS3BlZAl4AAJBeEQWFoXLlyqZFixb25/nz5ye5ffHixfa8W7duJleuXFd8fAAAAEgeQW+YRo8ebbJnz24mTJiQ4Hotbps0aZJdwDZkyBDPxgcAAIDQCHrDVKFCBRvcrl271jz11FPmzJkzNtWhQ4cO5ujRo2bWrFmmaNGiXg8TAAAAQRD0pkLHjh3NsmXLzLp160yJEiVM3bp1TenSpc2mTZtM/fr1vR4eAAAAQqB6Qyo1aNAgkMMLAACA2MBMLwAAAHyPoDeKqOVeWrTdQ9pj38Uu9l3sYt/FLvZd7KoSw/uOoBcAAAC+R9ALAAAA3yPoBQAAgO8R9AIAAMD3CHoBAADgewS9AAAA8L04x3EcrweRXuXOnducO3fOXHPNNfby9u3b7bl7GbGDfRe72Hexi30Xu9h3sWv7Fdp3ep7MmTOb48ePR2ybBL0eKlasmDlx4oRtZQwAAIC/7Ny50+TMmdPs3bvXRApBLwAAAHyPnF4AAAD4HkEvAAAAfI+gFwAAAL5H0AsAAADfI+gFAACA7xH0AgAAwPcIegEAAOB7BL0AAADwPYJeAAAA+B5Bb5TYunWr6dSpkylevLgpXLiw6dChg9m2bZvXw0IqnDp1yrzxxhumVKlS5pdffvF6OEjB999/bzp27GiKFClismbNasqVK2ceeeQRs2fPHq+HhhTs2LHD/O1vfwvsu/Lly5uBAwea06dPez00pNLw4cNNXFycee+997weCsJUp04du88Sn2rVqmWiHUFvFJg/f76pUaOGyZYtm9m8ebMNmPLkyWNq1qxpvvzyS6+HhxQcP37cjBw50lx99dXm0UcfNbt37/Z6SEjBggULTO3atc20adPMgQMHzNmzZ20g9dZbb5nq1aubDRs2eD1EJDNBoH33ySefmIsXL9p9pwmCIUOGmNatW3s9PKTCN998YwYNGuT1MJAKixYtMmvWrAl62/PPP2+iHUGvx/TLun379nam4p133jF58+Y1OXPmNOPGjbOzvnfffbc5ePCg18NEMhYuXGhuvvlm07ZtW6+HgjAcOnTIHlWpW7eumT17tg2ili1bZlq2bGlvVxCs4OnkyZNeDxWJnDt3zv6+HDx4sP2yqd+N2n/169cP/F/89NNPvR4mwqBZec3Wa58idgwdOtTOyv/4449JTvfcc4+JdgS9Hnv88cfNiRMnzMMPP2wyZswYuD5Tpkymd+/e9pd6//79PR0jkteuXTtTr14907NnT6+HgjBMnDjR/nJesmSJadWqlf3C2ahRI/PZZ5+ZLl262Pvs2rXLTJkyxeuhIpHJkyebYcOG2d+XSmsQ7b8ZM2aYfPny2csrVqzweJQIxzPPPGP3XenSpb0eCsK0atUqO2nQtWtXU7FixSQnpThEO4JeD/38889mzpw59ufbb789ye3Nmze35x999JH9oCG6aZYe0W/t2rVmzJgxSX5B6/KoUaPsF05Zv369RyNEKN27dw/MyMen3F7N3Isb/CJ66QunUov0BRSxNct77733GsdxTKwi6PXQvHnz7LnSGZQPmpi+OWk2QzlrmslAdHODJUS3ESNGmBw5cgS9TYtIq1SpYn9Wjj2iS5YsWULepv2lo2WdO3e+omNC6hw9etT06NHDvPvuu6ZQoUJeDwdh0iSA1h8pB1sTPEoRi8WjKgS9HlL+mWi1fzD6BV6yZEn7c6jEcQCpc+2114YVWFWtWvUKjQiXS3mhWvSrCg4VKlTwejhIhlJT7rrrLtOiRQuvh4JUUFqR688//zRTp061aWHKsT927JiJFUxNeZzeIG5gG0zBggXN9u3bzZYtW67gyID0SYft9P9SR1hiYVEG/vLiiy+a++67z7zwwgteDwXJUMWNdevWkToUg78X27RpY4NcLb5XeopKPsr06dPtIrbly5fbeCXaEfR6SKvEJXfu3CHv4y7WOHLkyBUbF5BeabZQ+fOPPfaYyZ8/v9fDQQr27t1rA11VvqlWrZpZunSpady4sdfDQhC///676du3r13HEiq9CNEpLi4usMjXNWvWLPt7UiVWf/jhB1uJQ+kP0Y70Bg8dPnzYnif3C0B1KIWi60DaGz16tC0VqMPkiP591axZMxvwynfffWfuuOMOM3PmTK+HhiAzhX//+9/NP//5T3PjjTd6PRxEgMo6qp65vmy6tc9V+jHaEfRG6aIMl1vDsECBAldgRED6LsejGq+qQUkFgOinWaZvv/3Wpn+phJKcP3/elnpUGUhED1VLUV1llSmDf+TLl8+uTXLTGmKhRjZBr4fcD4ra14biJoizyhVIOwqSevXqZbt6abYQsUPto99//32779y0sVg4zJpeaD2Kcq4//PDDBLXo4Q9FixYN9BLQF9BoR9DroUqVKtnzffv2hbyPW583WEkzAJGhgLdBgwbm2Wef9XoouERPP/20KVasWIJFwvDeq6++av+O6TB4rly5kpx27txp7/fAAw/Yy27JQMSONm3a2PPMmTObaEfQ66GGDRva8x07doScfXJbEDdt2vSKjg1ILwYMGGBz5sePH+/1UHAZ9AfXLYOVPXt2r4eDeCl6Fy5csH/Pgp3cRgdnzpwJXIfYUqZMGXt+3XXXmWhH0Otx+1p3VatWISe2cePGQAWHJk2aXPHxAX6nQHflypW25TCHXmOfO9Nbu3Ztr4eC/6cceQW2oU5uwKQW07qsagCILUf+v7qUygZGO4JeD1WuXDkwMxEsB23x4sX2vFu3bvawD4DI0R9ZnVR6J1j3NS2KUg1KxI6tW7ea6tWrB1oSA7gy3WXbtWtnatWqZaIdQW8UlN3RobgJEyYkuF6L2yZNmmQXsLkLNBDd4peV06E6RC8tfBo5cqRdXKN9pTQinbQIateuXWbRokWmefPm7Mco88cff4RcLKNKDlpJrva2ACJnz549Ztq0aYEyq/FpTZJ+jyaOYaIVQa/H1DJTwe3atWvNU089Zf/IKtWhQ4cOtke5ZqG0OhLRTTlr2leuGTNm2JlCRJ9x48aZHj16mM2bN9v/f4ULFw6cihQpYkqXLm0rOHz99deBBRqInnUQaiOtmdy5c+fafFHVMv/ss8/sETEdMYuF2SYg1tY9dOzY0R6dVrqKJuX0N0//355//nnbkjhWyqoS9EYBfZhU1FntGUuUKGF/oesP76ZNm0z9+vW9Hh5S8Morr9jD41o97lI9Ss3gjx071tOxISF9GXnooYcCi2eSo4CXzlHRpV+/fnaxjIria/+ULVvWtGrVylZrUJ1lVeAAEFkvvPCCadu2rQ10+/TpY4Pf+++/337pnDhxop0siBVxTji//QEAAIAYxkwvAAAAfI+gFwAAAL5H0AsAAADfI+gFAACA7xH0AgAAwPcIegEAAOB7BL0AAADwPYJeAAAA+B5BLwAAAHyPoBcAAAC+R9ALAAAA3yPoBQAAgO8R9AIAAMD3CHoBAADgewS9AAAA8D2CXgAAAPgeQS8AAAB8L5PXAwAAAJHhOI6ZPXu2mTp1qrnqqqvM9ddfb7p16+b1sICoQNALAIAPHDp0yNx///3mwoUL5v333zfFihXzekhAVCG9AQCAGHfkyBFz6623mosXL5q5c+cS8AJBMNMLAECM69y5s9mzZ49ZunSpyZSJP+1AMMz0AgAQw5TDu2DBAtO3b19TqFAhr4cDRC2CXgAAYtiECRPsuRauPfjgg6Zu3brm9ttvNzNmzPB6aEBUIegFACCGqzV88cUXpnDhwqZMmTJm3LhxZsmSJSZ//vymbdu2ZvLkyV4PEYgacY7+xwAAgJhz4MABU6RIEdOiRQszb968wPX79+83pUqVMvny5TP79u3zdIxAtGCmFwCAGOUuWsuTJ0+C6xUI16tXzwa/OgEg6AWQTmi2a9CgQbaU07Jly0x6cPLkSTNy5EhTq1YtkzNnTlO+fHn7Huh6RNaxY8fMokWLrvjzKo1Bubw7d+5Mcps+6wqKNdsr3377rfn111+v+BiBaEHQC8B8/PHH9o9jXFxc4FSgQAGzffv2ZB937733mowZMyZ4XIYMGcw999xjosX58+dNz549zdVXX20GDx6cbg71qm5r48aNTdasWc26detsQKbXrvegdevWXg/Pd9UT1BSiQoUKQW8/deqUeeONN2y6wS+//BLWNrdu3Wo6depkihcvbvN1O3ToYLZt2xb0vt27dzf/+9//zO7duxNc//PPP5vmzZubLFmy2MvXXnutGThwoP0iBKRLyukFgDNnzjjTpk1zihcvrjx/e7r++uudP//8M9nHnTx50vnHP/5h79+1a1fn+PHjTrT5448/nBUrVgRe19KlSx2/6969u1OwYEHnwoULges++ugj+/qzZ8/unDt3ztPx+cXLL7/sNG3aNOj7qc+dbi9atGjgs7djx44Utzlv3jwnR44cdh8ePXrU/h/s1auXkzt3bmf16tVJ7q//c9WqVXPuvvtu5+zZs/a6jz/+2MmXL5+zefPmBPe9ePGi06NHD6dz5872/zyQnhD0Akjgp59+ssGS+0f6vvvuS/Ex69evt/fdsGGDE61OnDiRboLe/fv3OxkyZHBuuOGGJLd9/vnnQQOny7FlyxZn8ODBTnozfvx4p3Tp0s6RI0eC3j59+nT7Xj/44INhB736/5czZ06nevXqzvnz5wPXK6iuUKGCU6hQIefAgQNJHqfrFCRrn9evX99p1aqV88MPPwR9Dm3rpptusuMC0hPatgBIQIdAq1atapYvX24vf/LJJ6ZmzZrmqaeeCvkYN2fQPY9G2bJlM+mFSlapHW2OHDmS3HbHHXdE/Plefvlle+g+Pfnuu+/MI488YsaMGRPyc9+uXTt7rvQClRILx+OPP25OnDhhHn74YZs65FL6Ue/evU2/fv1M//79zaRJkxI8Tk0pwi1Ppm1NnDjR1KhRw9SpU8d069YtrMcBsY6cXgBB3XTTTTZHV5577jmzcOHCkPd17xfNlGucXvz222/2/Eq0o/3888/TXS1YHSVVHq0WkYUTMObNmzes7SoHd86cOfZnNZdITPm58tFHH5lDhw6Zy1GtWjXTsmVLG7hf7raAWJF+/goASBXV/dSiJ7lw4YJdVKM/yoh+x48fvyJfRjSjrMWM6a3c+/Tp082GDRvs/5HMmTOneP9wv3y4dXZVaUMLLxOrWLGiXZh49uzZiHRb0+I4fVZee+21y94WEAsIegGE9MILL5j77rvP/nz48GHTpk2bsMtdaRYsflWH9957L8Hh8Pi33XrrrQkeqyBKAYCqD/To0cNet2fPHtO1a1dbj1Qr2p999lkbjIuCgGHDhtnUDB3Sb9iwofnmm2+SHZ8OIetQcYkSJUyuXLnszNqXX34Z8v56/UrxqFy5sg1KChYsaA9fqwxUYgokxo4da66//npbIuzcuXOmb9++9jC4XqsblKbEXfVfu3ZtW3dVj2/QoIEZP3584LXHd91119n7DB8+3F5euXKlvayTunWF48yZM+b11183N954oy15pfemSpUqZsCAAXY8rlGjRtnPhnvdiBEjAs+VeOY3ta8jrfa/ynUp0NO+0xhuueUW8+KLLyb5/KXklVdesecKeiPJPZoSKlVE6Q4lS5a0P69Zs+ayn+/mm2+252+++aY5evToZW8PiHpeJxUDiD6NGjVyBg4cGFgAVqNGjcBCnI4dOya5vxbnBFukM3/+fCcuLs7eNnny5AS3HTx40LnxxhvtbXo+15w5c5yqVasGnq9bt252oVyBAgWckiVLOlmyZAncNmDAAFs94rbbbrOLf0qUKBG4rXDhwnble3zubTNnzrSvSWPTgi/3+syZMzszZsxI8vq+//57p3z58s7EiRPt8+k9GT58uH18tmzZ7OIw17Bhw5yrrroqsE29jz179nTy5MkTuE4LnFLy22+/ObVq1XLq1q3r/Pjjj4FFTo0bNw68Z6EqZeg5E7+v4Th9+rR9Pr3HCxYsCDynFlVpe+3bt0/yGD2H+zoj8TrSav/r86bH9+3b11ZD0CIxVfSoXLmykzdv3rDfo23btgWeQwv4wuH+/0hpIVuVKlXsfW6//faQ99ECNN3nlltucSLBXbQ6adKkiGwPiGYEvQCSDXpl586dCcoujRw5MqygV9xAJHHQK08//XSS4EzByaFDh2yQqdsaNGhgyytt3LgxUAbKDbS0kr1Lly7OW2+9ZQM2Wb58uQ1edfu4ceMSPJ87/ooVKzpjxoyxwc+pU6ecsWPHBoIpBad79uwJPEb3ueaaa+x9ElOJNj1G742CL5d+LlKkiL1Nwd7o0aNtoPzEE084LVu2dA4fPpzs+6+yUnqNKlG1d+/eBLdpOxqPtt2mTZuIBr3aR3qcVv/Hpy8Cul5fEOK/zpSC3kt5HWm1/0eMGGGv0+MTB7EKkMP15ptv2u3o8xJu2bdwg173MxNqv8rNN98cKCcYCbVr1w67SgsQ60hvAJAiHW799NNPA0Xun3nmGbN48eKwHptczqPyExPToWc1xqhevbq9nD17dvPhhx/aihKSO3due2hbDh48aHr16mVXurvb0iFrHd6WH3/8MejzKlf5oYcesmkKqurw4IMPBlbX//HHH+bdd98N3Fer5NXtSs0HElObV1HTh/gL/TTmcuXKBV7jo48+ag+7K3dSC5WU+pEcLVRS9QwtNCpatGiC27SdIUOG2J+V16m82kj5888/7Xni51QuqagihNI8wnUpryOt9v+OHTvsuVuVxHXNNdeYu+++O+zXpEYfoi5okV4o6L63wapuuLQP5PTp0xF5Tvdz+vXXX0dke0A0I+gFEJb69eubt99+2/6sXMqOHTuG3V3qUijYcYOLxJUXypYtG/Rnl3I+4wdxiSmvNDHljrqLh+K3KVZApq5uyqV0c1bd05NPPmmDLZ0Sd69zvyCkNl9U3KD7hhtuCHp727ZtAyXYFBBGigJIBfnK6XVpX69duzZwWe/FlXgdkd7/buCuz61yWJVn7XrnnXfCfk0//fRTIPiONPczkxx33PpiEAnudnbt2mXzuQE/I+gFEDYFhqojKipzlJqFbZEsMZZScODOLgdbJBVK/AV1bskv2bhxow2StdAn8UkL0jTjptMTTzyRZHuXQjN5q1atsj+Hqv+qQFGL5CSlBXupoe1qH2tmXwvHtJBRX3bc8Ui4lRou93VEev+rFbVmi7WAUYsKtejtrbfeSvWMqdo7p1XQq1luib9gMLFjx44F6vJGgjurrP3FYjb4HUEvgFSvXHcbHChQ0eygX7gzhPEDLgUCCvCv1CzY/v37A7N5yQU/7lgjPS59iVEA36hRI1vLdfXq1ck2JonW15GYgtSvvvrK/Otf/7Iz80pZUY1aVeP44osvwt6OO043zSCSKlWqFEiXCcWtqRuspNmlcGfUU/slEYhFBL0AUkVlk6ZNm2YqVKhgL0+ZMsU3dT7dIE2H1OPPhOmQvsp/JSdSAYPyjBPnoSY3QxfJTmgKVFVWTPnayvFs3759gq5gsfI6khvTq6++alMU1N1Mr01ju/POO81///vfsLahEm6SFkc43FzkUO+XZqmVxyxNmzaNyHPGn+mO5o6KQCQQ9AJINf1xnD17dqDTlA4Th+Iu9kluZixamhsor1Fuu+22BLVvJbnAXjOIidvCXs6MpPuFIrlA263126RJExMpffr0MZs2bbILxVJabBfNryMY5V+7h+8VYKtGsAJ7fcFRnd+hQ4eGtZ3SpUsHAtBIc9sW//7772bv3r1JbleqjWimOlLvl/v+a38lt4AO8AOCXgCXFIQqGJw6dWqKM4Fu7mOwP+JqFiAKOrymmVotYNN41WLW5a7snz9/fqDpQ7B0gGbNmkVsLGrCIOvXrzfff/990PuoMoFyWzt37hz0tVzKlwm3AkVy+bSJv7y49w2WG3u5ryPSs/izZs1KcF2tWrUCX9jcz2JKlPLh5n1H+suaUi3chhf6vCXmVkxR62N3xvlyuTPHoRYbAn5C0Asg6GIZd8FMcpo3b25GjhyZ7H3cVfMffPBBYIGYZrJUAsyduVLlA+VKxg+o3DzQYLme8SsIBAuY3aAvVC5psFk6dVBTYK6V/fGrO6i0mTqTiWZAW7dubebOnWvzmZXmoRJZNWrUCMwAJh5jakp8uR577LFAKSnNUAabWVZLaI0nWFqAu9jKDWjC5c7uqmOemzuq2dD4XwKUArFixYpAJzr3S41mc/W+b9u2zXY5u9zXkRb7/7nnnkuyP9xOdeEGfe5iR21bi/3CEf8LQUq5y6NHj7Z5thMmTEhwvZ5PRxO0gM0t9RYJW7ZsSZBaAfia14WCAUSPCxcuOLNmzXIyZsxoGwds3bo1rMd17949ZOH9RYsWBbqyqWlA6dKlnaxZszrvvfdeoImCTuXKlXNGjRplH7N7926nePHi9nrdP36zCI0x/uNee+01e138DmBu0wN1Rvv1118Dt6npgq5Xdy91oFKjBTUYePvtt50cOXIEbUAhX331lZM/f/7Ac8Y/qXlC4oYN6iCmDmHuGDZt2uSklrp9uY097r///kBzh1WrVjlly5a1Xd6CNUfQPnM7wqmZxLx588J+TjXPcF+XOs1pH1x99dXOsmXLAtere1nr1q0D77k60Lm3qbmCxqxmJpfzOtJi/z/66KP2upo1a9r9KRpLs2bNbJOLcD/rGqvur20tXrw4xfur85vbGEMndfJLqanFlClTnEyZMjn9+/e3TTf0+lu1amU/g3rfIkVj037WuDZs2BCx7QLRiqAXgPWf//zHBqOJg7pcuXI5+/btS/ax+sNcr169kN2mPvzwQxuIZM+e3WnYsKGzdOlSe72CF7VeVecsdT5zgyg3SHZPutynTx8bwLjdtuKf1B1LQcsDDzyQ5DY9dujQoQna3LZr184pVqyYfb0Ktnv06BFokRuKArlevXrZAE7Pp9ej8ScOeNVNK1hwrO5tqaXObQp81J1MwYkC0BYtWjizZ88Oev86deoEfW4F9E8++WSKz6fX8vDDD9ugTl8Mevfu7Rw5csTepnbA6qym99jtfuY+Rl3R9Bx6/nXr1l3W60ir/e8Gve5Jn0UF3XpM/KA6HIMGDbLbeOmll5K9nzoXKnhNPCZdp46AyVm5cqVtr6z9UKZMGbtfUjvOlKxevdqOR+3AgfQgTv94PdsMAECsUNc+NcVQibH4NYxjzYABA2wqymeffWbuuusur4cDpDmCXgAAUkm532ovrQV6WoAWa5QXrXx7ndQaG0gPCHoBALgEnTp1sgs+582bZ2KNWkRrAaEWZLoNQgC/o3oDAACXGDiqGsTrr79uYsnWrVvNSy+9ZBYsWEDAi3SFoBcAgEugZg5Lliyx5doi1ZwkralEmdovz5w505baA9IT0hsAALhM//73v01cXJzp0qWLiVZr1661wW6/fv1oOYx0iaAXAAAAvkd6AwAAAHyPoBcAAAC+R9ALAAAA3yPoBQAAgO8R9AIAAMD3CHoBAADgewS9AAAA8D2CXgAAAPgeQS8AAAB8j6AXAAAAvkfQCwAAAN8j6AUAAIDvEfQCAADA9wh6AQAA4HsEvQAAAPA9gl4AAAD4HkEvAAAAfI+gFwAAAL5H0AsAAADfI+gFAACA7xH0AgAAwPcIegEAAGD87v8A+B8Vk9ImdWAAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = mp.set_figure(figdpi=200)\n", "colorlist = color_cycler = [\"#4477AA\", \"#EE6677\", \"#228833\"]\n", "N_max = time_list[-1, 0]\n", "exp_max = int(np.log10(N_max))\n", "\n", "x, y = time_list[:, 0] / 10**exp_max, time_list[:, 1]\n", "popt = np.polyfit(x, y, 1)\n", "ax.plot(x, np.poly1d(popt)(x), c=colorlist[0])\n", "ax.plot(x, y, \"o\", label=f\"ovito, k={popt[0]:.1f}\")\n", "\n", "y1 = time_list[:, 2]\n", "popt = np.polyfit(x, y1, 1)\n", "ax.plot(x, np.poly1d(popt)(x), c=colorlist[1])\n", "ax.plot(x, y1, \"o\", label=f\"freud, k={popt[0]:.1f}\")\n", "\n", "y2 = time_list[:, 3]\n", "popt = np.polyfit(x, y2, 1)\n", "ax.plot(x, np.poly1d(popt)(x), c=colorlist[2])\n", "ax.plot(x, y2, \"o\", label=f\"mdapy, k={popt[0]:.1f}\")\n", "\n", "\n", "ax.set_title('Neighbor Search')\n", "ax.legend()\n", "ax.set_xlabel(\"Number of atoms ($\\mathregular{10^%d}$)\" % exp_max)\n", "ax.set_ylabel(\"Time (s)\");" ] } ], "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.11.13" } }, "nbformat": 4, "nbformat_minor": 5 }