{ "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.3\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": "iVBORw0KGgoAAAANSUhEUgAAAr0AAAJHCAYAAACdLe6nAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAewgAAHsIBbtB1PgAAu79JREFUeJzs3Qd803X6B/BPRtORMsveUPYSQTYIsmTv6UIEvVPPU89151/FcXqenqfe6TlBcbH3lCUge4nsVUbZe6Yjzfi/nm9JmrRJm5a0WZ/36xXSfn+/JN82pX3yzfN9Ho3dbreDiIiIiCiMaQM9ASIiIiKiwsagl4iIiIjCHoNeIiIiIgp7DHqJiIiIKOwx6CUiIiKisMegl4iIiIjCHoNeIiIiIgp7DHqJiIiIKOwx6CUiIiKisMegl4iIiIjCHoNeIiIiIgp7DHqJiIiIKOwx6CUiIiKisMegl4iIiIjCHoNeIiIiIgp7DHqJiIiIKOwx6CUiIiKisMegl4iIiIjCHoNeIiIiIgp7DHqJiArg9OnTqFmzJu69917Y7fYC34/c9ueff0a/fv3QtWtXn25z5MgRvPDCC0hISMCxY8cK/NiUf+fOncNbb72FKlWqYNWqVYGeDhHlA4NeIgpZ06ZNQ4kSJaDRaJyXv/zlL17Pv3jxIqpXrw69Xu88Py4uDo888ki+H3vdunUq4Fy6dCkuXbpUoPlnZGTg6aefxrhx47BgwQJYrdY8b/PNN99g9OjR+Ne//oXLly8jmM2YMQOdOnVC+fLlUatWLTRp0gTPPvssVqxYgT//+c/YvXs3QsnixYvx6KOP4rXXXsOpU6cCPR0iyi87EVEIs1gs9unTp9tLlSoly63q8v333+d6m0OHDtmNRqO9e/fu9rS0tAI97tWrV+09e/a0//nPf7bfrp9++knNu1OnTj6dn5GRYY+NjVW3OXr0qD0YPfvss/bo6Gj7pEmT7FarVY3JXEePHu18nnbt2mUPRXXr1lXz/+WXXwI9FSLKB670ElFI0+l0GDp0KKZPn+4ce+yxx7B9+3avt6lduzYaNWqEAQMGIDo6ukCPKyvMsvL38ccf43aVLVs2X+fLSnXp0qURrBYtWoQPP/wQL7/8Mh566CFotZl/amrUqIFvv/0WL730EkJZfp8vIgoODHqJKCwkJiY6g+DU1FQMHDgQFy5c8Hp+bGwsjEYjgoEEsUVxm6Ly448/quvmzZt7PP73v/8dFStWRKgK5u89EXnHoJeIwsp7772nrk+cOIFhw4bBYrEEekoRx5FrPGXKFK9Bozw3RERFiUEvEYUV2cjm2Ji2evVqPPPMM/m6/YYNG1QlhQYNGqiV4NatW2P+/Pk5zluzZg1GjRqFevXqea2w8OCDD6Jq1aooV66c2sQlm9B8sX79erRv315tspPr/fv353r+77//rio/yHwldUNSLrxVlFiyZAl69uyJpk2bqrfp27Zti0mTJuU47+bNm+p+JCVBqhTs2rULd955p9qUtnLlylznc8cddzhXfMeMGaPuK7uWLVt6vG1KSgrefPNNtGjRQj2WrAg/8cQTHjft/fLLL+jevbt6vFKlSqnnYvz48WqDoCv5/skGtPr166vP//Of/6jKF3fffTdMJpPzvEOHDql0jLp166JSpUrqOfv0009z3WCYnp6u0jXkeykXqexAREEqPwnARETBSjZJOX6lpaen2zt27OjcMDVx4sQc58umsW+++cZtTD5v0KCBfc+ePerz3bt326tXr27XaDRu5z733HPOzUxyPDu5nWys69Chg/3KlStqI9crr7yizpdxua08jmMjl2yIcmxkW7p0qT0uLk7dr16vV+ONGjWy22w2t8eQ43Js6tSp9uLFi6vPZeOY42t++umnc8zrjTfesFeoUMG+efNm9fm1a9fsw4cPV+c/8MADzg1ncp916tRx3pdsFJTPHfORc3Nz6tQpe+nSpZ23r1q1qrqPvMj3qnnz5vbx48erzXpms9n+t7/9Td1H48aN7Tdv3nSe++OPP6rn5Z///Kf63GQy2fv27avOffnll53nPfroo85Nf/I9ktuVKFHCObfly5er85YtW2YvWbKk2ngn32t57AEDBqhzxowZ4zZPeZ5kXJ4r2cwoz2nZsmWd9zlt2rQ8v1YiKnoMeoko7IJecf78eXuNGjXUmASDGzduzDXo3bdvn91gMOQ4TwIYuY9ixYqpig0O69at8xr03n333erY2rVr3apM1KxZ0xlEunIEvXJ85MiR9hMnTqjx/fv3qyoTcmzbtm0eg16pQCHVKMTly5ft/fr1cwZfMkeHn3/+2eMLAHmB4AhwP/jgA+e4BJ0SrMq4vIA4cuSIffv27fZRo0bZN2zYYM/Lpk2b7OXKlXPORS7t2rWzr1+/3uttHnroIRVEupIAVIJ+uf1rr73mHG/WrJkak3k5rFq1So3JMVeOr71MmTL2xx9/XD0X//vf/+xjx461p6am2k+ePKkC3j/84Q9ut5s8ebJz7hcvXswR9LZv315VCpEXC3KRAFvG5TkgouDDoJeIwjLoFTt37rTHx8er8UqVKtnPnDnjNeh94okn7JUrV85xvxJ4elrBO3jwoMegV1ZPZQVSjsnHrv7yl7+ocQnuPAW9TZs2VQGZKwkC5djMmTM9Br2HDx92G5fVUEewKUGYQ8uWLdW8ZDU1OwkAHUGhawk3Cepk/NNPP7UXhLzwePDBB90CX7nI2KVLl9zOPX36tF2n09m/+uqrHPcjLwTkdg0bNnSODR06VK1aS9DqIN8Lx4sHV/KiQMblBYTMKTtZuZfjW7ZscRuX508e++GHH3augrsGvdlXdOUFk4w3adIkX98nIioazOklorAlOZmSVyols6SD2pAhQ2A2mz2eK3mqkjcqeZ+uF8nvlfxPuZw9e9Z5flRUlMf7kXxSb/m00hhDuN6PK8lLleoTruRxhWvuqavs50ter5RsE7/99pu6Pnz4MLZs2aLKrJUsWTLHffTp08fZvGPTpk05qhQ0bNgQBSE5rt99953Kk5bcZIfvv/9e5fSePHnSLUdacmclnzf7c/Drr7+q78P169ed50uJOtmsGBMTA5vNpsqkPffcc+qYfO7K8XWUKVPGY7kxua2oU6eO23jx4sUxefJklYvtKLuW/evLz3NFRIHFoJeIwlr//v3x9ttvOzeIPfXUUx7PS05OVkGybHpyvRw8eFAFg3LxdtvsgY8jeDp+/LjbMUcwnJ9yXY6ALXsglxvHhq2rV6+q671796pr6UDnSbVq1Zzl2+T74G9t2rTB2rVrVQc92djn2OgnnegcHI/7xRdf5HgOJDiW778EudnJ+dIK+sCBA3j11VcLND9HK2dfOuL5+7kioqLDoJeIwt5f//pXPPDAA+rjL7/8Ep9//nmOc6S0mayIelulzY9///vfamVw6tSpbuM7d+5U19JMI7/yMy+pFiEcq7pSEUFcuXJF1TD2xHGurG7eLmk1nL2CgpAyZfI9kIoR4ueff3a283WUlpMXGb44evSoCqYloJcWzvKYjpXW/HIEqUlJSfAHf/wMEZH/Megloojw1VdfqfJj4s9//nOO4EpWXyW9YdmyZR5vL8Gi61v/uenbty9++OEHVRprxYoVKgiSUmGSajF69Gh1vDBJcCuaNWvm1rhDyIqoJ46gs3Hjxrf9+PK98tYRT4Jr+d440gUcq7eO1W9ZDfbGUSotLS1Nre4KKatW0K56DlWqVMn1sWWO3n4uiCh0MOgloogguZ9z5sxRAY6sQp45c8btuNRsFZLCcOnSpRy3/+c//+nz29+yevnKK6+oFU2pG1yrVi1Vv1UCNF9r9d6Obdu2qev77rvP2RmtcuXK6mPJUc1OgkjpXnfXXXepufqzK5sn8hiSvyzpFo48544dOzpTUP773//muI3Mb+LEierj5cuXq5q6kpbhSX7TCzp37ux8YeQp31rq8ErdZiIKbQx6iSgsODYPSQDnTYUKFTBv3jzV9CE7Wf2V1UdZAW7VqhV++ukntcK3Z88eFfTISq28ne7g2BDn6W18yVWVjW5yH1u3blVpE5LT+oc//MFjXq3jPnILqj09jqe30m/cuKE2j8mmvXvuuce52c2R7yqNKK5du+Z2m6VLl6pA8Z133vEYPDrSI/JDcm0lgPVEUhzkhYVsEnSs8MpqtORfO56LP/7xj2q1WJ4DaQ4iX8vIkSPd5iXpETt27FAfy/f4hRdecOYyy/fLsXKb19chjyf5uPJ9kWYXMm+5vfwsyAsH+blwrAb78nx5e66IKMCKqEoEEVGhevPNN31uDCB1crM3nBBSpzZ7eS25SA3XAwcOuJ37xRdfqGNarVaVNXOVmJjo8X6ioqJUWbQnn3zSfuPGDef5jgYMUoLrwoULbrVypVmDHBs9erTbY3Tu3NnZ0MIxNykD1qNHD1W7V5o1ZK93K6W3HLV9z507p8Z/++03e7Vq1eyvv/662/lSl9ZR+kxKrGUvpZYbqXcrt5OmEO+9956qH+ywevVqe+3atVXt3ezlw86ePevWFMP1IiXlHOR2UjdZxqXMmdQTlnJr8+bNc56fkJBg/+mnn9T5UkvXMS5NKDz5z3/+4/FxpYmIa41eKXsnzSjk2FtvveV2H7NmzXL+TGQvJUdEgcegl4hC2ueff+7W/UsuFStWtC9YsCDX20l3suxBr1i4cKFqxiBd0aRz1+DBg3MEtRJwOmrxykWaWrz00ktuHdmkQ5cEcBIgSbCbPZgaOHCgOtfRQMNxkUYa0oThu+++U8G267Hy5curZhJCglDpHtalSxc1TwkkpQOcBOO5BajSnOKuu+5St5G6txIkL1682O2cDz/80K27m1zkfKl77AvpBieNIqSz27hx41QHulq1aqnAVDqryQsU1+5qriRwl9tXqVJFfV/r1atn/+ijj3J0pJOGExKQSmB977332pOSktT4sGHD1ON8+eWX6vNWrVrl+N7fc889Hh97/vz59jZt2thjYmLUC5CnnnrKrSGJ3KfMKftzIkG9NMSQYNf1eXz33Xd9+n4RUdHQyD+BXm0mIgonr7/+OurVq4dRo0a5jcvb7PIWutSkvf/++1X+aHx8fMDmSUQUSTKLChIRkV9II4UpU6ao+rLZSW6obOAaMGCA2hjlqJhARESFjxvZiIj8SDZTycax3N5Ek9qysbGxHrujERFR4WDQS0TkR9JpTIJaaUAhlQVcg1+p+DBz5kxVp/df//pXQOdJRBRpmNNLRORH0oXtoYcecpY0k5xd6ZAm+bynT59W5dJmzJiBrl27BnqqREQRhUEvEZGfSdezDz/8UHXxOnnyJAwGg2rIIHVpn3nmGZQpUybQUyQiijgMeomIiIgo7DGnl4iIiIjCHoNeIiIiIgp7DHqJiIiIKOwx6CUiIiKisMegl4iIiIjCHoNeIiIiIgp7DHqJiIiIKOwx6CUiIiKisKcP9AQiWYUKFWAymVCtWrVAT4WIiIgoaCQnJ8NoNOLs2bN+u0+u9AaQBLwZGRnOz5OSktSFQg+fu9DF5y508bkLXXzuQldSET13Eh9JnORPXOkNIMcK7549e9R1o0aN3D6n0MHnLnTxuQtdfO5CF5+70NWoiJ47x+P4E1d6iYiIiCjsMeglIiIiorDHoJeIiIiIwh6DXiIiIiIKewx6iYiIiCjsMeglIiIiorCnsdvt9kBPIlKxZAsRERFR0cRIXOklIiIiorDHoJeIiIiIwh6DXiIiIiIKewx6iYiIiCjsMeglIiIiorAXkUHv7t27MXLkSJQrVw7R0dGoVasWnnrqKZw5cybX2x08eBCjRo1CxYoVUbZsWYwYMQKHDx8usnkTERERUcFEXNC7ZMkStGzZElOnTsWFCxdgNptx9OhRfPLJJ7jjjjvw22+/ebzd4sWLceeddyImJgb79+/HsWPHULx4cTRv3hwbNmwo8q+DiIiIiHwXUXV6L126hNq1a6NZs2b4y1/+gvr16+P06dN4//33sXDhQnVO1apVVVAbFxfnvJ2s5spt5Lbbtm2DTqdT4xaLRdWRu3z5Mvbt24cyZcrkaz6s00tERESUE+v03qavvvoKAwcOxMqVK9GvXz/UqVMHnTp1wvz583H//ferc06cOIHJkye73e7ZZ5+FyWTCk08+6Qx4hV6vx2OPPYaLFy/ixRdfLPKvh4iIiKgoZRxIgvXUWYSiiAp6t27dik8//RQajcZtXD7/8MMPVRArtm/f7jx25MgRLFiwQH3crVu3HPfZs2dPdf3jjz+qlWQiIiKicGO3WpE652fc/PdXuPHlT7ClpSPURFTQ++6777qlLbiSjWmOpXTJ23VYtGiRujYajahZs2aO20mKhGyGk9zg2bNnF9rciYiIiALBeuESrv/9P0hb/Iv63H72PK49+wZM30yD5cRphIqICnolJzc3BoNBXTdu3Ng5tnTpUmeuryeS7lClShX18ebNm/04WyIiIqLASt/0G66/8SFsp8+5H7DZYN64HTfe+QTmzTsQCiIq6M2N7OeTVAZZtZW8XwcZE47A1pOEhAR1feDAgSKYKREREVHhsqemwfTNVKRMnApkWLyfaLOFzIpvZhIrqbJjkpP7zDPPoFSpUs5xKWsmihUr5vW2EiiLK1eu5Ptxk5KSnGkVnrCyAxERERUly9ETMH09GbaLl327gc2G9BVroX94eL4eJ7f4R+KjxMRE+BNXem/56KOPVNOJ8ePHu41LOTLhLRdY2Gw2dZ2WllbIsyQiIiIqHHabDamLf8GN9z7zPeC9xbxtl7p9MONKL4B169Zh1qxZatNayZIlc+T5Sj3e3GRkZKjr0qVL5/ux5VUMV3PDk+SDS9c+KXn32muvIVjSeFasWIEvv/wSO3bsUF0GiYiIbFeuqXQGy4HMtM58M2dkpkFEZ+6P8kVu8U9uq8AFFfErvVJ/d9y4cXjzzTfRo0cPr/m6qampXu/j2rVr6jq/zSkovM2dOxdXr17FtGnTECz+8Y9/4P/+7/8wffp0VXGkKEiXw/79+6u0odjYWNXZ8H//+58KwG+H3F5qbI8ePVq1FZda2d46KhIRkXfmHXtw/a2P3QNerRbR/XsAhij4RM6LCu611OCeXRGQgLd9+/Z4+eWXPR5v0KCBalhx7ly2XYsuHPV5PZU0C1c2mx3pFiui9Tpote51jynrZ0tqPj/++OMej69ZswZ33313kc5Jfs6l3nTr1q2L5PGkEUzv3r2Rnp6u6mDLuyaywiyNXlatWoUpU6ZAq83/a2/ZYPrggw+q+5Ta264VV4iIyDd2cwZSZyxE+uqNbuPaMqVhHDsS+lrVYD9/UVVpyIuhRRNoCvD7vCgF9+wKmbzlLHm4X3zxhddzOnTooK6PHj3qdaVYOrKJ7t27I9wlnb6K96ZtxoDxs9H/tdnqWj6XcXInK5qyQfKhhx7ymAcugV8gFCQNpyAk0H344YfxyCOPqP8/kgYk6RT33nuvOi6rzV9//XW+73fjxo1o2bIlGjZsqIJqBrxERPlnPXUW1//xSY6A19CqGYq/8mcV8Irobh3Uqm+uZFW4a2a8FMwiNuiVQHft2rWq5bBra+HshgwZoq5Pnz6Ns2dztt3btWuXs4JDly5dEM5W7kjGk58sx7Ltx5GWYVVjci2fy7gcJ9/TDHbv3h2Qx3Z0HixsS5YsUTnNkspQo0YNNSatv+fMmaOuxQ8//JCv+9y3b59aOZZ3ZyQvObf/u0RE5Dk1LO2X9bj+zifutXejDYgbM1yt8Gpis5p06atWgnHMcO+Br1arjst5wS4ig95vvvlGXSTn0rX7moO8BTtjxgz1sawm9erVS328ePHiHOcuX75cXUteYXx8PMJ9hddq85yHKeNc8fXNxIkT8eqrryLcSa67pw188n/uT3/6k1tJQF/IuzJ9+/ZV/z9lhTh7O3EiIsqd7aYJpv99h9Qp8yTYcY7ralRRq7vRbZp7vJ2s/hZ7+U8wtG2eleNriFKfq/FWzRAKIi7onTRpEt577z18//336u1XSU2Qi/zxldzdZcuWoWfPnuqYazkz2YAjK0uuZHObBDCygU02woWzmWsPeg14HeS4nBcq9u7dq16s3HHHHahQoYLK33799deRkpLiPOfnn39Wq4kSYDku2Tv7NW/e3HlMNmwJ2SQm+aqdOnXCmDFjnOfKz56s8jo2ccl9ySV7Trm8hT948GA1t3LlyqlUCfk5zKuSyO0oUaKE8+uQ1WBHak9BSVqHt/rWjpXe6tWr+3x///znP1Uu7xNPPKG+J0RE5LuMfYdx/c2PkLFzX9agRoOYnp1R7IU/Qlcu9834asX34eEo+fEbKPmfN9W1fB4KK7wRGfR+9tlnKgDZv38/6tati7Jlyzov8ke0WrVqqoLDpk2bMGjQIOft5FwJbrdu3YqXXnpJBcSS6iBv3crufFkxLl++PMJ509qvu076dK6cJ+cHO1m1b9u2rUpJkY1VJ0+exNixY/HGG2+gTZs2zs2Jkn8qPy+O4ExeEB06dMjtvmSzWp8+fVSO7rx589TPxmOPPaY+l81qrlUKpMKA6+0PHz6sLu+8845zTH7W5HEeffRR/P7770hOTlbBs5Q+k/nkVknkdsjPsQS8zz33HM6fP6/Sf6S1dv369fN16dq1a56P5ciDHzBggM+rvB988IH6uF27dup7IY8j7cElOJf8YCIiyslusSBl1mLc/HgC7NduOMc1JYoh/umxiB3UE5p8pL3JZjVNtCHoN61FdPWG2bNnqxUiX0jAm70ZhZREkj+w0ryiUqVKagVL3mqV1UJZJQxnUqXBkcObFzlPzo81BO+PlgSlDzzwgGo3LSu9QlY2n3/+eZWj/d1336kXRxLAOlYlZSVfzpVgLfvb6larVQXG8i6CkJ+Hb7/9VgWyUh4sPyTPV6o9/PWvf3Wm1Ug6gKzySnAuG7fkhdd//vMf+NPx48fV/48ff/wRo0aNco63atVKfW3+JmlB8kLT0yY/b/nBN27cUHWz5YWEBObSGlxeoA4dOhTDhw/Hu+++q743RESUyXr+IkxfT4H1uPvCVdQdDRD30FBo442IJKEXpheQBLKy4ubLRYIeT2TzjPyxllXAY8eO4ZNPPgn7gFdIWbKYKN82DMl5cn4w+/e//6067Un6QHYSbAqp/7pt2zbnuKzqyzsCstovwaerBQsWoGPHjs6azg5yfn5JcC2pEZ7m5gjoPv/8c4+bKgtKAm2pPCIrqa4Bb2GR8n9Su1ieB6PRt1+4q1evVtfyQuBf//qXCniFlF6T/HwhKSKB2hxIRBRMJJZJ37AN1//+H/eAV69H7KgBMD7+UMQFvBEV9FLBSR3ejk0yg4y8yHnBXrf3p59+UteOigKuJK/XUW954cKFznGpziHpBkJe7GSvBOKpFm9UlI8FvV3K3zlWlz3NTerryjyk9Jd0e/NXN0JJE5CVacfKcmGTFxbyIvT+++/3+TaSfiIqV67s8fsiq/FSBs6x2k5EFKnsqWkwTZiClG+nA+lZTYi0lcqj+Mt/QkznthG7EZhBL/lkSIe60OURzMpxOS+YXb9+HadOnVIfe/tPL4GvkFxaV3/4wx/UpjYJmq9cuaLGZMVfUh4kDeB2JSUlOTdQepqbBNHSttrT3ApCNsXJyq7k78oKtieSPuDYbOfrpXPnzl4fU8qVSUpQfuvzyvMmihcv7vG45FQLuW8iokhlOZKsOqtlbPndbTy6c1sU/9ufoKsc/u9O5yZ4Ey8pqCRWKokXh7fyWrZMAl45LucFM9fKDBL8Nm3aNMc50i7XU4AlGx379eunAjfZbCZ5pVLRQ4LhwpibpwDP29wKQvKYJ0yYoDbHydfSrFmzHB3iJH1ANtr5g9TYffvtt9UKevac+bw4UkUcwW92jnSH221tTEQUiuw2G9KWrELa/OWy+9w5rjHGIW70UBjuaBjQ+QULrvSSz7o0q4ZP/9QN3ZtXd+b4yrV8LuNyPNhJ8OQoo3XgwAGP5zjKgjVp0iTHMUd9WakEIquyskHSX3mwtWrVcn5ckLkVhKQ2SCAq9zts2DBnGoG/SXMXacss9a8LUm7srrvuUtd79uzxeNxRb1sqrRARRRLb5au4+eFXSJu71C3g1ddLRPFXn2bA64JBLxVoxXfuG4Mw781B6joUVngdJD1BOnoJ6cbnidRrlioBUt3BU5AoZbkkFUFWeGWVNL+rlt5IMCgrq3nNTTZPZl+RvR2yQU46D0qag2ygk/Jg/iT3Ky8MJKXBU11eby2+Xcm8JL1D0i1kE6K3+/C0AZCIKFyZt+9W6QyWgy6/R7VaVYYs/pmx0JYqEcjpBR0GvVQgsllNypIF+6Y1T1555RX11r7UoJUmEK4kP1fyW6V8mSOVIDtH6TvZNPXHP/7R6+NIFQYhG8+yk2YnInuAKXWCHTVzJV/Y1c6dO9WKqVR4cLTflc1bUl9YOgdKwxVfOObl2uhCNrJJLvOWLVvU/fkz4JUV5I8//tiZK+0gjy9tiKXphIOUJZMXG7LR7ebNm26pJbLKLt+v7OXaZMV96tSpqraxP18MEBEFK7vZDNMPs2D64gfYU7Jqt2vLJqDYi4+rhhOhWEe30NkpYBo2bKguVPQmTJhg12q19urVq9t37Nihxi5cuGDv0qWLvXfv3vbU1FSvt7127Zo9Pj5enZubUaNGSYKpvU6dOnaTyeR27K677lLHli9fro7961//ch4bP368OtasWTP7kSNH1NjRo0ftTZs2tY8dO9Zus9mc527dulWdK5eoqCj78ePH8/zaP//8c3W+Xq+3Hzt2zDk+a9Ys53099thjuX4PfCFzqVu3rj0mJsaekJDgdilVqpR6fHmsZcuWOW8zffp05xxmzJjhdn/p6enqey63k/OEfO9Gjx5tr1+/vnr+iIjCXUbyKfvV8R/YLz/2ktvl5sSpdltqmj1cNCyEGIlBbwAx6A2sNWvW2Hv16qUCMAlMW7VqZf/f//5nt1gsed728ccfdwZe2UnwVa5cOWfwJpdixYrZ58yZ4zxn06ZN6jHLli1rf/bZZ+03btxwu4+5c+fa7777bnvJkiVVQNexY0f7Tz/95DEAl+NyXxUqVMgRKGZXpUoVt3lJQProo4+qY/JYrscksJ84caK9INLS0uzVqlVzuz9Pl4oVK9qtVqvzdmfPnrXXq1dPBcvycXYSiL/88svqvitVqmRv3Lix/YUXXrBfv369QPMkIgoVsuCRumKt/fITL7sHvH9+zZ6+6Td7uGlYCDGSRv4p/PVk8qRRo0a5bs4hyg9JC5D2vNKymIiIwoft+k2YJk2HZbf7Jmddzaowjh0JXVn35kjhoFEhxEgsWUYUBq5evaoqL7Rt2zbQUyEiIj/K2HsQpm+mwX49a58DNBqVtxvTrxs0t/Z4UN4Y9BKFOOnkJuXApI2wVJ0gIqLQZ7dYkDrnZ6Qv+9VtXFOyOIyPjEBUvcxmReQ7Br1EIUwaR3z44Yd45513WKOWiChMWM9dgOnrKbAmZ3YQdYhq1hBxDw6BNt4YsLmFMga9RCFM2v5++umngZ4GERH5gWyzMq/fhpSp84D0zPKSSlQU4ob3haFjK49t6sk3DHqJiIiIAsyWkoqUH2YjY9tOt3Fd5QowjhsFXaXyAZtbuGDQS0RERBRAlsPHYJowRbUUdhXdpR1iB/eCJioqYHMLJwx6iYiIiALAbrUibdEvSFu4QnIbnOOaeCOMDw9DVJP6AZ1fuGHQS0RERFTEZFVXVndlldeVvkEdGMcMg7ZE8YDNLVwx6CUiIiIqQuZtu5Dyw0zYU9KyBnU6xA68F9HdOkCj1QZyemGLQS8RERFREbCnm1VlBvO6rW7j2nJlVGc1fY0qAZtbJGDQS0RERFTILMmnYPp6MmznLrqNG9q1QNyI/tDERAdsbpGCQS8RERFRIbHbbEhfsQ6ps5cAVmvWgZhoGB8YDEPLOwI5vYjCoJeIiIioENiu34Dp2+mw7DnoNq5LrK5aCevKlA7Y3CIRg14iIiIiP8vYfUAFvPYbN7MGNRrE9O6CmD5doNHpAjm9iMSgl4iIiMhP7BkWpM5erFIaXGlKlVCru1F1awVsbpGOQS8RERGRH1jPnFeb1awnz7iNRzVvjLgHBkNrjAvY3IhBLxEREdFtsdvtMK/dgpSp84GMjKwDhijEjegHQ/uW0Gg0gZwiMeglIiIiKjibKQUpP8xGxvZdbuO6qhVhHDcKugrlAjY3cseWHxTxfvjhB/Tp0welSpVCtWrV8Pzzz+P69esIRps2bcIDDzyAevXq3db97NmzB08//TRKly6NY8fcW2ASEZFvMg4dxfW3Ps4R8EZ3bY9iLz3JgDfIMOiliPbkk09izZo1WLBgAQ4ePIgSJUrggw8+wOjRoxFs/u///g/3338/fvzxR6Snpxf4fqZNm4ZXX30V//nPf3DlyhUEyrBhwwr0dt+pU6fw0EMPoUqVKoiNjUXr1q2xcOHCQpkjEZEndqsVqfOW4uYHX8J+5ZpzXFMsHvFPjUHc8H7QRPHN9GCjsUsiCgVEo0aNnKtuoVhsGxkWIEofsj3ClyxZgl69emHdunVo166dGjt8+DB69uyJTp06YcKECQg2a9euRceOHVG9evXbWqGV//bFihWDyWTC0aNHUaNGDRT16vqDDz7onIuvTp8+jVatWqmV+Pj4eJw5k7lZRKvVYv78+ejdu3ehzZmISFgvXoZpwhRYjyS7jesb1oVxzDBoixcL2NzCSaNCiJH4MoTyxXLiNNKXr4VZ3soxZ6gkfUPzJoju1gH6qpUQSr7//nt1nZCQ4ByrXbu2CnyDVfny5f1yP7LCWrJkSRX0FjVZqX355ZcLvDL/2muv4ZFHHoFer8euXbswePBg9Zy98sorDHqJqFCZt/wO0w+zgDSXd9t0OsQO7onoLu1DdhEoUvDZIZ+ZN+/AjXc+gXnj9syAVw1mqM/V+OYdCCX79u1T1waDAaEiKirKb/clQWNRk1VdCVglcM2v5ORkDB8+HI899phz7k2aNMEXX3yhPj5w4IDf50tEJOxp6arRhJQjcw14teXLothfn0BMt44MeEMAV3rJ5xVe0zfTAElr8MRmU8e1FcuFzIqvI5+VZWSKzmeffaZeZIwbNw6PPvpovm4rmwzlkl2bNm3UdbNmzfw2TyIiB8uxkyqdwXb+otu4of1diBvRH5ro0Fk4iXR8WUI+kZQGrwGvg82G9BVrEcxSUlJUCoNcTp48qcY6d+7sHDObzeryzTffoHHjxvj2229x4sQJleMr6QA//fST2329+eabaNGihUo7qFixIp544glcvnxZHV+/fr2qjiBBtVwefvhh523ffvttxMTEOI+tWrUqx1zlfv785z+jTp06KodXgrqlS5cW6vdnwIABzjnJRb6us2fP+uW+JQXh/fffx9dffw1/kvnJXAuyekxElNvelbSfV+PGe5+5BbyauBgYH7sPxoeGMuANMVzpJZ/+46scXh+Yt+1CnPwiCNK3eeLi4pw5u7J56/jx4yrgdGzk+uWXX/CXv/wFO3ZkpmpIlYRBgwapRPq0tDS1Unnffffh6tWr6Nq1K/r166fKiMnb9uPHj8c//vEP/Prrr9i4caPaHHfp0iW1YUsqLmSvxCABct26dXHxovvqgZBAWzasSWkyuS/JO168eDGGDh1aqN+fKVOmqMeSIP/jjz9W8xP169fP932tWLEClStXVh9brVZVceHDDz/0W16yg1SheOedd3Dvvff69X6JKHLZrl1X715a9rnv8dDXrgHj2JHQli4ZsLlRwTHopbxJlQZHDm9e5Dw5P0Rf/d5zzz347bffVMC6YcMGlS8qFwkE33jjDVVmS0iN23LlyuH11193W72dN28edu/ejffee0+dLyuQXbp0yRH0CqkLLMGkVGRwZbPZMGLECLXJbOrUqWqFWUilCakhLKvLhUGCenlcST346KOPVEUEh/3799/Wfcv3QwL4gQMHwl9kpf2tt97CxIkT8e9//1sF1jqdzm/3T0SIyMpC5p37kDJpBuw3XTb6ajSI6dsVMb3ugYa/Z0JWcC7HUXCRWoMGHzdQyXlhUJuwVq1a6rpbt25o2bKl+nzSpEno27evKpMlQeyQIUPcbiMBrmysEjNmzHCOuwaP2XkK0uS2EnCPGjXKGfC6ph8UhmvXrqmVUknVkJXT3OacXzt37lRpIrJy7M+SZ/LC5N1338WNGzdUfrAE7BaLxW+PQUSht+/k6tPjcfXPr6lrtVJ74rTP92HPyEDKlHkwfTrJLeCVVd1iz/8BsX27MeANcQx6KU/yalnKkvnC0KJJ0KY25IejOkDDhg1zHJNmFrKqKCuuslLrepHUBlkVvp2ObhJcC6lHm53kCPubBPF33323CuglRcOfJD9aGn18+eWXKF68uN/uV7rSSQqKrMpLNz0xc+ZMlYtNRJHFH5WFrKfP4cY/PkX6L+vdxqNaNEWxV59WaQ0U+kJ/SY6KhNThVb84ctvMptUiumsHhDspnSUk7UFSDvxty5Yt6rpMmTIoCs8995xajZXcZG9kk19+rV69Gv/9739V7rPkCBcG2dwn3fQk7URWyCdPnpzvqhBEFLmVhWQ/hnnNJqRMX5CZFuEQbVCVGQztWrDCTxhh0Es+kV8WxjHDvf9y0WrV8VApV3Y7HG+hS9viwgh6HaXUMjJ8zKO+TdJ2WRo8yEqpbAjz1DiioA07pFqD5CjLY3jj+IMimwilkkZByKq7BL3SsY2IIkd+KgvpHx7uPnzThJTvZyJjx163cV21yjCOGwld+bKFMWUKIAa95DNDq2bq1bL88pAqDc6ObC2aqBXeSAh4hZQmE9OmTVMb2jxZuXKl2sAm8rtKICkM58+fV5UliurrmT59uprvq6++iubNm6tWzP4g5da8cTSTkA1ujsoaBdWgQQP1fXM8N0QU/m6nslDGgSSYJk6F/ap7Klp0946IHXgvNAFo3kOFj88q5YsEtvJqWX553O4O2UCTt7UKQkqJOerwytv3Tz31lNvxCxcuqIoCjqA3NjZWXTvq97o+vmzCcuS+OrRt2xZz587F8uXL8ac//cmvc/emQ4cOajVW6gJLSbbNmzcXKKUhu9yqPjheDNxuZQjH6rtUcyisjX5EFB6Vhex6HdLmL0faklXyi9R5WFM8Xr1bGdUws0wjhafQjFYo4CTQlaLcoRrwips3b6prR+DpSt6SFxJIZZeYmIj+/furjyVI/OMf/4jt27er2rrz589XZc9GjhzpPF8aSwipB7x3b+bbaNIY4/7771cruo70AUcg6wiiJVd127ZtHueemprq9vmnn36qaupKSTPH3PPiCLQd6RryuLJBTNIr5OuTig6B8tJLL6lAXOoj+0JW3aXWsrQoJqIIkc/KQtar13Hj/c+RtvgXt4BX37geir/6DAPeCBC6EQvRbZD8UcfKq2t5MUeg62hOMWvWLFW/NjupRuB46142tEmpL2mRK8GibNqSSggOcqxmzZoquG7atKkKgiVAlQBTAmjxzDPPqOYNEojKxq8XXnhBVYiQ+5MVZSENMf7617+qj6WhhXRnkw1oQsYPHTqkVmultm9e9u3b5+y05rqBTcqKyQqsHJcUh0DkyMrXJnV9161bh88//9ytqoWUcJP8Y8dmPwnw5fn73//+h59//vm2UiSIKHwrC0me7o13/gvr0RNZg3odYkf0Q/yfHoa2eHzhTZSCBoNeijhSWkwCS9dNUBUqVFCrhBJAlS1bFrt27XJ2FZMqCtINzZV0FZNOaZLTW6VKFRgMBpWbKk0dPvnkkxzlz6RpRZs2bRAVFaVuu2TJEvTu3VvV6ZXKA5JO8P3336v7ERL0TZgwQeWpysqx5NlKYDt27FgV2ElNXWmC4ThfynbJPBs1auQMkr2R9AWpJ+xYWZZg21FfWO7fMS5fXyBWT6XkmwT78rVLgOsg3wN5sSDfO6nRK3m8w4cPVyvT8iJGXnQQUeRVFpKN1LnSANbDx4C0dOeQtkJZFPvrk4jp0p7VGSKIxu7v5EDymQQowte3cInyIkG3pCX4u94uEVGwknKauZYty8bQsRXihveF5taiAUVOjMSNbERhZNGiRX7tfEZEFJKVhaRzmtXqdp4mLhZxDw6BoXnjgM2VAotBL1EYkDdsXn/9dfTo0UO97U9EFImVhWL69YDpm6mwHjrqfrxuTRjHjFAthSlyMeglCgN/+MMfMHToUBX0EhFFIvPve5EyaQbsJpeqO1otYvp1Q0zPziFdbYj8g0EvURiQahJERJHIbs5A6oyFSF+90W1cm1AKxrEjoU/MLBtJxKCXiIiIQpL11Fnc/HoybKfPuY1HtbwDxvsHQRMbE7C5UfBh0EtEREQht48hfdVGtcKLWw12lGgD4kYNgKFNc5YioxwY9BIREVHIsN00qdzdjJ373MZ11avAOG4kdOXKBGxuFNwY9BIREVFIyNh3WFVnsF9zaR+v0SC6x92I7d8dGj3DGvKOPx1EREQU1OxWK1LnLkX60jWS2+Ac15QopkqRRTWoHdD5UWhg0EtERERBy3r+IkwTpsB67KTbeFTTBoh7aAi0xeIDNjcKLQx6iYiIKCg3q5k3bkfK5LlAujnrgF6P2KF9EN25DTerUb4w6CUiIqKgYk9NQ8pPc2DevMNtXFupPOLHjYKucoWAzY1CF4NeIiIiChqWI8kqncF28bLbeHSnNmqFV2OICtjcKLQx6CUiIqKAs9tsSFuyCmnzlwM2m3NcY4xTubuGZo0COj8KfQx6iYiIKKBsV67BNHEqLAePuI3r69VS1Rm0pUoEbG4UPhj0EhERUcCYf9uNlO9mwp6SmjWo1SJ2QA9Vf1ej1QZyehRG+JNEVEiSk5Px0ksvISEhAceOHUM47KRevnw5hg8fjrp16wZ6OkQU4uxmM0w/zobp8x/cAl5tmdIo9uLjiOnZmQEv+RV/mogKwQ8//IBRo0bhvffew+XL7psxQtU//vEP/N///R+mT58Os9mlfFAhW7lyJTp16oSaNWuifv36+Oc//wmr1Vqg+/rtt99UiaPsl9atW/t93kTkneXEaVx/5xOY12xyGze0uRPFX/kz9DWrBmxuFL4Y9BIVggceeACrVq1CTEwMwsXLL7+M//73v0X6mJ9//jl69uyJJ598EkePHsUvv/yCiRMnYtiwYQUKfCVw9+S5557zw2yJyJd3jNJWrsONd/8H25nzWQdiohH3yAiVv6uJDZ/fmxRcmNNLBWKz25CWkY6YqGhoNXzt5ElUVBRKly6N06dPI1zI11NUtm3bhj/96U947LHHVEqFqFixogq87733XnzwwQd48cUXfb6/gwcPqpXeffv25ThWp04dv86diHKy3biJlEkzkLFrv9u4rkZVGMeNhK5sQsDmRpGBQS/ly77zBzFhy49YcmAFUjPSEBsVg571umJsy/vRoBzzPD0FvuFEry+6XxmSDy2ruWPGjHEb79GjB8qXL4/XX38djz76KEqVKuXT/UlahKzoSooEERWtjL0HYfpmOuzXb2QNajQqbzemXzdodLpATo8iBJfoyGfz9i7BwEkPYvbuhSrgFXItn8u4HCfyh7Nnz6pcXkkPad68eY7jkoObmpqqcqd9cfLkSXXujRs3sH379kKYMRF5YrdYkDJjEW5+PNEt4NWULI74Z8chduC9DHipyDDoJZ9XeF9YOB4Wm+c8ShmX43JesEtPT8d//vMf9Va5VFWQ4Okvf/kLypQpgwoVKuDdd991nvv777+rt9Lj4+NRr149zJ8/3+N9yoqkvN3epEkTJCYmqk1XH374odc5yNvsAwcOxB133KEeV85/5plnVFDmSvJYn3rqKTRqlFmUXXJaW7ZsidjYWDRt2hQzZsxwnrtgwYIcG7RkPq4efPBB5zF5TH8pUaKE835lNbhDhw63dX+LFi1SuX/VqlWDzsMfRMfXJd8PX/zrX/9Sm+8kHaJFixZqtXfy5Mm3NUciyp313AXc+OdnSF+2xm08qllDFH/1aUTVSwzY3CgyMb2BfCIpDd4CXgc5PnHLT3i/z+sIVrNmzcIbb7yBnTt3OgPg3r17Y8+ePSoounTpEv72t7+hVq1aKuDq3r27evvcYrGonNChQ4di//79bgFjRkYG+vfvr3JFJfBs3LgxDh06hL59++L48eM55rBmzRp069ZN5arKfCS4e+KJJ/Dxxx/jwoUL+PHHH9V5Eoh/+eWXMJlMqF69utoYJ/cp80lLS8OuXbtUrqucL5Ui5NjatWvVfctx2ew1bdo0t8f+/vvvVZAv5dR8DRh9MXfuXHTp0kXNWTa8Se7v5s2b8dBDD+XrfipXrowVK1aoFwWiSpUqHs8rWbKkut6xY0ee9ynPjzxn8sJBnhd5ng8cOID77rtPzfvbb78Nqw2HRIEmv9PM67chZeo8IN2l0kuUHnHD+sJwd2v1ApmoqHGll3zatCY5vL5YfGC5Oj9YSYC7bt06Z67t3//+dzz//PM4f/68CngHDRrkXBmUFd+NGzeqAPHcuXNqVVYCpqlTp7rd5yuvvIIlS5aot88l4HVsjJLKA57ICrAEYkOGDIFWq1UrmRIoOlY4Hf79739j9uzZ6mMJYmV1evfu3eqtegni5DHkj4ts9kpJSVHntW/fXn09wlupNAne//rXv8JoNMIfJLCXoF2Cb/m+OTa7tWrVSs0zPxcJeIUE/67BraeVZXHx4sU85yfPtTw/8r27cuUKZs6c6dy4Js+lvPggIv+wpaTC9PVkpHw3wy3g1VWugOIvP4XoTm0Y8FLAMOilPEmVBkcOb17kPDk/WMmKnqQqSMMIIW939+nTxxkcyUqlkEBXVkkdaQUSZMnKoOOYw6lTp1QQK2+3Z39Lv3PnzihevHiOOdSuXVvNwbVigGNF89q1a27nSqqEkGD7q6++cq4wS6rFlClTnMGta9qFbNaS+5eVXAlwXckKt+TL9uvXD/4ggaSshktqh6w2+4u8ABFxcXEej8uLBceLgfyQ+xs8eLCat6zaO1a/ZYWciG6P5fAx3Pj7x8jYmvlOmkP0Pe1Q7G9PQlepfMDmRiQY9FKepCyZVGnwhZwn5we76OjMORYrVsxtvFKlSs7g2GAwuB2T3Fsh6QEOP/30k1q1lVXN7GQ1w1Nlgffff18Fqo5A99dff8W4cePUx7Jy6ym4k+DZEag7yAavu+66y1ney0FWR6Wqgc1mU8GoqwkTJqgawv6oKiEr5l27dlXpAb169YI/Ob732b8fDo7mGAUtoSb3L8+d43lzvIAgovyz22xIXbAcN/71BWyXrjrHNfFGGJ8cjbiR/aEJs0o2FJoiPuiVAEbeNq5atarPrWJl57inrk6yQSYcSR1eKUvmi171uoV03V5HkJnbMddAbMuWLW4Bsa9ks5e8tS6b5GRF9q233irQfB1B79WrWX9ohGyKk8f47rvvVGqGI1CUFISxY8fidkmOs6zsSlrI1q1bPZ6zadMmtaqdn4usjgvZUCgkn9kTx9eb3++7Kwn83377bfVxUlJSge+HKJLZLl/FzQ++RNr85fLL0Tmub1AbxV97GoamDQI6PyJXEbuRTXbJf/bZZypv0hEU+GLZsmVqg44nktsZrqQO7/y9S3LdzKbX6vBIy8wUgEghOaJCVnt9JXmokiohq8DS0tdTCoSvHCud2VesZRPeiBEjVJArL+okuJszZ45Kl3CkbNwOCahl1ViCdkmnaNasGe6+++4cLw4PHz5coPuXyhTCW2MPx/9ZybO+HdLeWL4WSQchovwxb9uFlB9mwp7ikmak1aoyZNHdO0KTyyICUSBE7E/k0qVL0bFjR5Xflx8SPMjbubJTP/tFSlCFK2k88X6fN1Rg64mMy/FIa1DhCDo9VWnwRn7mpPKC5JLeTsDruuGrQYOcqymODW3y4u7mzZsqSPXHKq+DpDbI/wdZ9ZVKEbLBzl8kmBby/8pTioMjmJaNibdDVntltbhNmza3dT9EkcSebobp+5kwffmjW8CrLZeAYi89gZh7OzHgpaAUsSu9snPekdsnQYGvOYyywUbKMEXi7tP+DXuiTplaqiyZVGlwdGSTlAZZ4Y20gFe0bdtWbXhbvXq1Wu31livrCNxkA5Xk8N555505coYdJBc3e5qFI4fVU/AnK5U9e/bMcUxWX2WTmbw78eqrr2L9+vVqZdnfXdMkxUMqIkgwL+XY/FH+S+royotS+V5J+TLXBhXy/ZF3W2R127EJ8XZeNEjZuvyWViOKVJbkUzB9PQW2c5kvuB0MbVtk5u7GBP+eDopcEf9SzFH6yBeyqiU7vr1trokEmSu+r2Pns2uw69lf1bV8HmoBryMdQVYpXUlA5Wg2kZ3jeXdNZZBmD/IzJCkLn376qdfHc2x+c9y/NL34+eef1cdnzpxRJb8c5L6yl0WT4MxR0cBBVkEl2JY2vd7q2b7wwgvq+qOPPlJBqevKssxFVn4bNmyoVp194Qi+Xb9v8s6HrDRL8OvPlWQpGSfBv+Qlu5KybvL9kHrKriXNJGVJ3m2RsnOysu0q+/fO9f/0J598grJly/pt3kThulktbfmvuPHP/7kHvDHRMI4bBePDwxjwUvCzR7ijR49KJKMu8rE327Ztc54XHx9vHzlypH316tW39dgNGzZUFypaBw8etGu1WvVc/vjjj27HJkyYoMYNBoM9KSnJOW6z2exjxoxRx+rXr283mUzOY3PmzLHr9Xp7VFSUur3VarVnZGTY3333XbtOp1O3+fvf/27fsGGDPT093V6tWjU1ptFo7FWrVrUXK1bMPnPmTHvx4sXVeOnSpe3vvPOO28+nnNu7d2/7+fPn1XhycrL9rrvust955532Gzdu5Pr1NmvWTN3HqlWr3Ma3bt3q/JmWuR8/fjzP793nn3+uzpev99ixY87xWbNmOe/rscces6emptr9Qb4P0dHR9uXLl6vPd+/eba9SpYp96NCh6vvsavr06c45zJgxwzn+1FNPqbEHHnjAfuLECTV2/fp1++uvv26fNGmSX+ZJFM6s167br388wX75sZfcLtfe/dRuuXAp0NOjMNWwEGIkBr0+Br1Dhgxxnud6kT++V69eLdBjM+gteq+99poKaB3PnwSTjRo1UscaNGjg9txKYCeB7t69e+1lypRxOyYvfJYuXeq837Vr19o7d+5sj42NtdeoUUO9KJJALTEx0d66dWv7K6+84nyRJMGmBKxybtu2bdULKvHcc8/ZS5QoYX/jjTdy/HxWr15dBZx169ZVQbN8/uKLL6rgLS9yu9q1a6vA3dW1a9dUAF+nTh17hQoV3AJFTyTYdP0exMTE2B999FF1rGTJkjm+PxMnTrT7w+TJk1XgXqtWLRXkO15YZHf27Fl7vXr11PdIPnbYuXOnvVOnTurFhdFotHfv3t0+fvx4ZwBMRN6Zd+23X3nuLfeA9w9/tafM/dlus1gCPT0KYw0LIUbSyD+IYFKmzFHw/+jRo6hRo0aOc+RbJDU9pbaq5FCuXLlS5WY6yG54eZs5ex3VvDh20UsLXKLcfj6lDbGvJfWyk5q9ch+Orm+eSFrAs88+q6oZEBHZMyxInb0E6SvcG7doSpWA8ZERiKpbK2Bzo8jQqBBipIjdyJYfsmnt/vvvdxubO3euqoUqgYg8IVLwf/Hixfm+b6kPmlsJKQbEdDskt3XGjBm5/hxJzVupvCCb8oiIrGfPq1bC1hNn3Maj7myEuAeHQGv03CmRKL9yi38kPnJ0JfWXiN/IVlADBgxQu8od9USXLFmCVatWBXpaRG4+//xzVT/X0WkuO2n+IN3gpHObt2oSRBQZ5F3N9LWbcf3t/7oHvFFRiHtgMIx/eIABL4U0rvTeBtk5LvV+5ZWK7A6fNWuWs6OUr+RVDFdzyRtHpQgpq+ULKUkmVQ/atWuHWrVqYfz48aqMmCeSqvPhhx/inXfeQd26oVV9g4j8y2ZKQcoPs5GxfZfbuK5qRVWdQVehXMDmRuFrTy7xjz8aKWXHoPc2lS9fHi+++KKqV8pWpuRvjoBV2v1KQ4smTZrkev6kSZOwfft2dRFSn9dbe2xp+5tbmTUiigwZh47CNGEK7FeuuY1Hd22P2EG9oIliqEDhgekNfiB1QYW3xgRE+SU1ZytXrqxSDxw1daWhRV5d/55++mmUK1dObVz73//+hzfffLOIZkxEocZutSJ13jLc/OBLt4BXU8yI+D89jLjh/RjwUljhT7MfyM56Ua9evUBPhcKEdBs7depUvm8nHdjOnTtXKHMiovBhvXgZpolTYU1yb6Gub1gHxoeHQ1uiWMDmRlRYGPT6wZUrV9T18OHDAz0VIiKiXJm3/I6UH2fDnpqWNajTIXZQT5XSoMnWBp0oXDDo9QNpizpkyBCvuZNERESBZk9LR8rUeTCv3+Y2ri1fFsZxI6GvVjlgcyMqChEf9KalZb3S9bZD/syZM2pDkbx1XLp0abdj8lby999/r2qhEhERBSPL8ZMwfT0FtvMX3cYN7e9C3Ij+0ESzZCGFv4h+D8NqtaomEw6zZ8+GxWLJcd5rr72GkSNHomHDhvj222+RmpqqbivNKF555RVMmTIlRzBMREQUaHabDWlL1+DGPz9zC3g1sTEwPnofjA8NZcBLESNi2xC///77qi1r9iBXr9fj448/xhNPPOEcS05OVi1aZbX3+vXrqtC/dK+SQLh///4FngPbEBMRUWGxXbsO0zfTYdl3yG1cl1gdxrEjoUsoFbC5EQUiRorYoDcYMOglIqLCkLFrP0yTpsN+w5Q1qNEgpk9XxPS+BxqdLpDTIwpIjBTxOb1EREThwp6RgdRZi5G+cr3buLZ0SbW6q69dI2BzIwo0Br1ERERhwHr6HEwTJsN68qzbeFSLpoh7YBC0cbEBmxtRMGDQS0REFMIkS9H862akTFsAZGRkHTBEIW5kfxja3QWNRhPIKRIFBQa9REREIcpmSkHKdzORscM971FXrXLmZrUKZQM2N6Jgw6CXiIgoBGUcSFKthO1Xr7uNR3fviNgB90ITxT/xRK74P4KIiCiE2K1WpM1fjrQlqyS3wTmuKR4P45jhiGpYN6DzIwpWDHqJiIhChPXCJZgmTIX1aLLbuL5xPRhHD4O2eHzA5kYU7Bj0EhERhQDz5h0w/TgbSEvPGtTrEDu4F6K7tOdmNaI8MOglIiIKYva0dKRMngvzxu1u49oKZWEcNwr6qpUCNjeiUMKgl4iIqAjYbTYgwwJE6aHRan26jeXYCZi+ngLbhUtu44aOrRA3vC80BkMhzZYo/DDoJSIiKkSWE6eRvnwtzNt3AeYMVT/X0LwJort18LpKKwFy+tI1SJ27FJBg+RZNXCziHhwCQ/PGRfgVEIUHBr1ERESFmYf7zTS3wFUCX0lVkGNSbcHQqpnbbWxXr6tSZJYDSW7j+jo1YXxkhGopTET5x6CXiIiokFZ4cwS8rmw2dVxbsZxzxdf8+16kTJoBuykl6zytFjF9uyKm1z0+p0UQUU4MeomIiAqBpDR4DXgdJI1hxVro7huE1JmLkL5qg9thbUIp1VlNn1i9cCdLFAEY9BIREfmZ5OSqHF4fmLfuRMaxk7CfOe82HtXyDhjvHwRNbEwhzZIosjDoJSIi8jep0iCb1nw81y3gjTYgbtQAGNo0Z+1dIj9i0EtERORvUXpVpcHnwPcWXfUqMI4bCV25MoU2NaJIxaCXiIjIz2TDmZQly95QIjfRPTohdkB3aPT800xUGPg/i4iIqBBIHV4pS5bnZrb4OMSPG4WoBnWKampEEYm1T4iIiAqBlCGTOrxScswbXdVKKPH6XxjwEhUBrvQSEREVEmk8Yb14GWnzl7uv+Go0mekMg+7lZjWiIsKgl4iIqBDYU9OQ8tOczBQHF5qK5VDs0fugq1whYHMjKiibzY50ixXReh202tB6wcagl4iIyM8sR5JhmjAFtouX3cajO7VB7NA+0EhlB6IQknT6KmauPYhfd51EWoYVMVE6dGxSBUM61EVipdBojc2gl4iIyI9NKdKWrEba/GVu6QwaYxziHhoCQ7NGAZ0fUUGs3JGM96ZthtVmd45J4Lts+3F17MXhrdClWTUEOwa9REREfmC7cg2miVNhOXjEbVxfrxaMY0ZAW6pEwOZGdDsrvNkDXlcyLserlyse9Cu+rN5ARER0m8w79uD6Wx+7B7xaLWIH9UT8M+MY8FLImrn2oNeA10GOy3nBjkEvERFRAdnNZph+nA3TZ9/DbkpxjmvLlEaxFx9HTM/OqlEFUahuWvt110mfzpXz5PxgxvQGIiKiArCcPAPT15NhO3PebdzQ+k7EjRoATWxMwOZG5A/pFqvK3fWFnCfnxxqCN7QM3pkREREFIbvdjvRf1iN15mLAYsk6EBONuFEDEd3mzkBOj8hvovU6VaXBl8BXzpPzgxmDXiIiIh/ZbtxEyqQZyNi1321cV6MqjONGQlc2IWBzI/I3rVajypJJlYa8yHnBXreXQS8REZEPMvYegumbabBfv5E1qNEg5t5OiOnfHRpdcK9yERWE1OGVsmS5bWbTaTXqvGDH7HoiIqJc2C0WpMxchJsfT3ALeDUliyP+mbGqQgMDXgpXiZVKqjq8Eth6IuNyPNjLlQmu9BIREXlhPXdRbVazJp9yG4+6o6FqNqGNNwZsbkRFpUuzaqoOb6h3ZNPYJSOfAqJRo8zOPHv27An0VIiIyIX8aTRv2IaUKfOAdHPWgSg94ob1heHu1tBogjt/kagw2Gx2VaVBNq0VZg5vYcRIXOklIiJyYUtJRcpPc5Cx5Xe3cV3lCjCOGwVdpfIBmxtRoGm1mqAuS5ab0Jw1ERFRIbAkHYdpwmTYLl11G4/u3BaxQ3pDY4gK2NyI6PYw6CUioohnt9mQtvgXpC1YIe/fOsc1xjjEjR4Kwx0NAzo/Irp9DHqJiCii2S5fhWnCFFgOH3Mb1zeoDePDw6EtWTxgcyMi/2HQS0REEcu8fRdSvp8Je0pa1qBWi9iB9yK6e0dotKzsSRQuGPQSEVHEsaebkTJtAcxrN7uNa8slwDh2JPQ1qgZsbkQUBkHv8ePHsWTJEmzatAmHDh3C5cuXYbVaUbJkSVSoUAHNmjVDhw4d0KVLF2j56pqIiAqB5cRpVXvXdvaC27ihbXPEjRwATUx0wOZGRCEe9C5fvhxvvfUW1q5d6zbuWiJY6h3Onz9ffZyQkIBHHnkEzz//PMqUKVMUUyQiojAnf3PSV6xD6uzFgMWadSAmGsb7B8HQqlkgp0dEodycYt++fXjsscewfv16twC3fPnyKF26NOLi4hAbG4u0tDTcuHEDFy9exKVLl5znGY1GvPnmm3j22WcRjticgoioaNiu34Bp0gxYdh9wG9fVqqbSGXRlSgdsbkShxm63F3pzlpBqTjFjxgyMGTMGMTExGDt2LDp27IjmzZujVq1aKtD15sqVKypYXrNmDRYtWoQXX3wRCxcuxKxZs1C8OHfQEhFR/mTsOQjTt9Ngv34za1CjQUyvexDTtys0Ol0gp0cUMk5fuok56w/jzKWbeOvhDgg1hbLS+9///hfjx49XKQ2y0hsVVfBi3ocPH8Ybb7yhIv0VK1agVKlSCBdc6SUiKjz2DAtS5yxB+nL31DpNqRIwPjICUXVrBWxuRKHCbrdjR9IFzF53CBv3n4YjavzsqW6oXblUZK/0yka17777Dr///juqVr393a+1a9fG999/j+nTp6s8X1nxZb9zIiLKjfXsBbVZzXritNt41J2NEPfgEGiNcQGbG1EoSM+wYuWOZBXsHj17LcfxWesO4cXhrRBK/Br0pqSkqFVeWZH1dyrCsGHDVA7w559/jscff9yv901EROGzKmVetxUpU+cB5oysA1FRiBvRD4YOLblwQpSLi9dTMX9DEhZsSsL1FHOO4zqtBp2aVsWAtrURavwa9G7cuFEFvYWVe9unTx9Mnjy5UO6biIhCm82UgpQfZyNj2y63cV2VijCOGwVdxXIBmxtRsNt/4jJmrT2INbtOwmrLmflawmhA39aJ6NsmEWWKe9+bFTFBr9TXLWyjRo0q9McgIqLQIi2EpZWwtBR2Fd2lPWIH94TmNvaWEIUri9WGtbtPYda6g9iXfNnjObUqlMCg9nXQpVk1GKJCe9NnUHRkS0pKwgcffKCaV9x555147rnnwmrDGhERFQ671Yq0RSuRtnCl5DY4xzXFjDCOHoaoJvUDOj+iYHTdlI6Fm49g3oYklc6QnWQAtW1QSQW7d9QqGzYpQUUW9A4cOBDJycnqG9eiRQt8+eWXanzHjh245557cP36dZWLtXjxYrVpbevWrShWrFhRTY+IiEKM9dIVtbprTTruNq5vWAfGh4dDW4J/Q4hcHT17DXPWH8Ly7cdhtthyHI+L1qNXy5oqX7diQjzCTZEFvW3atMG2bdvw008/qZq9wmaz4eGHH8a1a9dUHrCUJqtWrRr+/ve/491338Xbb79dVNMjIqIQYt66Eyk/zII9NS1rUKdD7KB7Ed21AzRsZU+k2Gx2bD5wRlVb+O3weY/nVEqIx6B2tdHjrhqIiw7fVKAiC3pXr16tVncdAa+QAHjnzp1q9VeaWXTr1k2NN2vWDL1792bQS0REbuxp6UiZNl9VaHClLVcGxnEjoa9eJWBzIwomKekZ+HnrMdVMQppKeNK8djmVwtCqXkVoteGRwhAUQe+WLVvQo0cPtzFZzZWAV8qROQJeUbNmTZw6daqopkZERCHAknxK1d61nbvoNm5ofxfihveDJiY6YHMjChbSLU0C3SVbjyIl3ZLjuEGvRbfm1TGwXR3UrFACkaTIgl5JZXBt/iarvHv37oVOp8Obb77pdq60ITaZTEU1NSIiCmJ2mw3pK9YidfbPgNXqHNfExiDugcEw3NU0oPMjCjSJr34/ktk1bcO+rK5prsoUj8WAdrXRu2VNFDdG5gvEIgt627Zti48++gjPP/88du/ejWeffVat8j700EOoW7eu8zyLxYKnn366qKZFRERBzHbtBkzfToNl7yG3cV1idRjHjoQugZV+KHKZXbqmHfHQNU00rJagUhg6NK4MvS6yc901dtfl10K0a9culc8bFRWFK1euqJVfaVP822+/oXTp0uqcX3/9VQW8UtFBAmKryyv6cFQYfaWJiMJFxq79ME2aDvsNl3f+NBrE9OmCmN5doNGFds1QooKSMmMLNmZ2Tbtm8t41TYLd+lUzY6xQUxgxUpGt9DZp0kTl9f773//GsWPHUL9+fbXq6wh4//jHP+LcuXOoXr26uhARUWSyZ2QgddYSpK9c5zauKVVCre5G1akZsLkRBbprmqzqrt55wmvXtD6tEtFPuqaVCM2uaWGx0ks5caWXiMid9cx5tVnNevKM23hUiyaIu38QtMa4gM2NKJBd0yTY3Zt8yeM5NV26pkWHeNe0kFnplRVcqbOrLcT6iPIYNWrUKLT7JyKioifrL+ZfNyNl2gIgIyPrgCEKcSP6qwoN4dIVisjXrmmLthzFvA2HceGa565pbepXwuAO4dU1rTD5NeiVXN1vv/0Wr7/+Ogor4J0wYQLeeuutQrl/IiIqejZTClK+n4mM39xXdHRVK8I47j7oKpQN2NyIitqxc9cwZ91hLP/tONIzrB67pvW8q6aqxCBNJShAQe+dd96Jb775Bt99952qyuDvgHf06NGYP3++X++XiIgCJ+PgEZgmToX9ivvO8+huHRA7sCc0UUW29YQooF3Tthw8q7qmbT90zuM5lRKMqrZujxY1YIwJ365phcnvv00++OADdO7cWVVlGD9+PEqWLHnb9zl9+nS8+OKLmD17tmpXTEREoc1utSJtwQqkLf5Fchuc45pi8TA+PAxRjesFdH5ERdU1bem245iz/hBOXfTcNe3O2uUwOIK6poVU0CslyZYsWYL+/fur3NvHHnsM/fr1Q/v27fOV63vmzBksWrRIpTMcOXJEBb7SnpiIiEKb9eJlmCZMgfVIstu4vnE9GEcPhbZ4sYDNjagonLlswtz1h7B4i/euaV3vrK42p0Va17TCVCjvGxUrVgxLly7Fq6++in/9619q9VdWaGUnnpQqk2BYzpGL1OtNSUlRHdhOnDihAlzpyHby5El1X3fddRe2bt2KKlXYT52IKNSZN++A6cfZQFp61qBeh9hBvRDdpR00hbgRmijQmzV3Hr2I2WsPqq5pHiqOIaF4DAa0rY3erWqhRIR2TQvpkmXbt29Xwe/ixYszHzCP3YWO6dSuXRuvvfYa7rvvvkKtBhFILFlGRJHCnpaOlClzYd6w3W1cW74sjONGQl+tcsDmRlToXdN+v9U17YznrmkNqpVWq7odG1eJ+K5pIVOyzJPmzZtj4cKFavV28uTJKviV7mxmc84OIpUrV1b5wEOHDkXfvn2hY7cdIqKQZzl2QqUz2M671xg1dGiFuOF9oYk2BGxuRIXl0vVUzN+YhIWbjuCqyeWdDZeuaXc3qaKC3QbVEgIyx0gTkOYUFosFx48fVyXO0tPTER8fr+r7lioVWT3UudJLROHMbrMhfdmvSJ3zs2xPd45r4mIQ98AQGFo0Cej8iArDgZOXMXvtIazedQIWa84Qq3icAX1a10L/NrXZNS3cVno9Pqhej8TExEA8NBEReQlQkWEBovR+yau1XbsO0zfTYNl32G1cX7uGaiWsLX37lX2IgoVVuqbtOaVKju097rlrWo3yxTG4Q92w6poWalgAkYgogllOnEb68rUwb98lyYeqA5qheRNVJ1dftVKB7tP8+16kfDcT9pumrEGtFjF9uyKm1z3crEZh43qKGYs2H8mza5qkMDRLZNe0QGPQS0QUyZUUvpnmlnogga9543Z1zDhmOAytfC8VaTdnIHXmIqSv2uA2rk0oCePYUdAnVvfn9IkC5vi565i9/hCWb/feNe3eu2piILumBRUGvUREEbrCmyPgdWWzqePaiuV8WvG1nj6Hm1/9BNtp925SUS3vQNx9A6GNY+4ihUfXNKnCsC2XrmkD2tbBvXexa1owYtBLRBSBJKXBa8DrIBvRVqyF/uHhXk+RvdDmNZuQMn1BZk6wQ7QBcSP7w9C2Bd/SpZCWmm7B0u3HMGfdIZz00jWtWeKtrmn1K6qqDBScGPQSEUXgpjWVw+sD87ZdiHtoqMc8XNtNE1K+n4mMHXvdxnXVKqvau7ryZf02Z6LAdE07jCVbj8KUlpHjOLumhZ6IDnpTU1Px1Vdf4f3338evv/6qOsXl5uDBgxg/fjxWrVqlyq516dIFb7/9tmqkQUQUMmRFVjat+ULOk/Oz1dLNOJAE08SpsF+97jYe3eNuxA7oAY0+ov+8UABTENItVkTrddAWYMXV2TVt3SFs2HvKa9e0/m1row+7poWciPytdOPGDXz22Wf497//jXPnPOflZCdNNaRpxvDhw7F//35Vdu2ZZ55RzTd+/vlntG3bttDnTUTkF1F6VaXBp8BXzpPzb7FbrUibtwxpP6+WCME5rileTG18i2pYp7BmTeRV0umrmLn2IH7ddRJpGVbEROnQsUkVDOlQF4mVSvrUNe2X30+oYDfpzFWP59Svmtk1TRpKsGtaaCryoHfOnDn4/vvvVXOKrVu3qjGbzYY//OEPqFWrFp566inVrKIwLV26FB07dsSxY8dU8JuXw4cPY9iwYahTpw6+/vprZ6c4ue2aNWvQv39/1XGuTJkyhTpvIiJ/kFQFKUsmVRryIg0kHKkN1guXYPp6CqzHTridE9WkPuJGD4W2GHepU9FbuSMZ703bDKvLsqwEvsu2H1fHXhzeStXG9eTyjTTVNW3BxiSvXdMkeJZgtyG7poW8Ig16H330UUycOFF9HB2d9ZaAVqvFl19+qVZeGzZsiOnTp6N169aFNo8hQ4aoa4PB4FPQ++yzz8JkMuHJJ590a40sq72PPfYYnn/+ebz44ovOr42IKNhJHV4pS5brZjatFtFdO6gP0zf+hpTJc4A0l8BAr0fskF6IvqcdN6tRwFZ4swe8rmRcjlcvV9xtxffgySuYte4gVu9k17RIUmTr8xMmTFAXyZfx1PlYfmE+99xz6N27N3r06KHyZwtbiRJ5J54fOXIECxYsUB9369Ytx/GePXuq6x9//BGXLnnuwkJEFGykDJmkI0hg65FWq47rypRWubsp30x1C3illFmxvz2JmC7tGfBSwEhKg7eA10GOq/OsNqzZdRLPfLYST36yHCt+S84R8ErXtGcHt8BPf+uLR+5twoA3zBTZSu/nn3+OatWqYdKkSWoVt1Ilz3UfR40apVZ9X3vtNUyZMqVQ5yQrtXlZtGiRujYajahZs2aO4/Xr11er1unp6Zg9ezbGjRtXKHMlIvI3aTwhwauUJZMqDc6ObC2aZK7wWiy4/vf/wHbxsvvt7m6NuGF9oDG4b24jKupNa5LD64tffk/GjqTzXrumta5fUaUw3JlYji/iwliRBb179+5VG746dMh8qyyv1dcVK1YgGEj+r6hatarH45LuUKVKFSQlJWHz5s0Meoko5FZ8pQ6vlCVTVRpubVqTjWqyYc01/UFjjEPcg4NhuLNxAGdMlEmqNEjuri9kRTd7wBtrkK5pNTCwXR1ULsN89EhQZEGvbE7LK+AVEjgKyaENBpLeICSw9SYhIUEFvQcOHMj3/cvtGjVq5PX4nj178n2fRET5pTarRRtgu3INpm+mwnIg83efg75uTRgfGQltKdYjpeAgZcmkSoOvga9DxdJG1R5Y2gSza1pg5Rb/SHyUmJgYmkFvgwYNcPr0aa9pDeL69et499131VsLkjYQDC5cuKCuixUr5vUcx6a8K1euFNm8iIj8zbxjD1K+mwm7KSVrUKtFTL/uiOnZyWODCqJA0d6qrCBVGnwhXdMkhUFSGdg1LTIVWdArVQ5eeOEFteHLkxMnTqgauFJGTILeRx55BMHg8uXMXLa4uDiv50jJNZGWlpbv+5dXMVzNJaJAspszkDpjIdJXb3Qb15YpDePYkdDX8lzuiSjQOjWtiuW/HXctGZ2DpOj+36g26lwKLrnFP7mtAgd90HvfffepjV6S4jB69GjV0WzZsmU4efKk6nA2Y8YMZ9DYtWtXPPHEEwgGUtZM5pqbjIzMAu+lS5cuolkREfmH9dRZ3PzqJ9jOnM+xyS3uvoHQxMYEbG5EnkgFqF23uqat33sq14BXVnSlTi8DXiryOr2TJ0/GK6+8gqeffloFuI5yX44SZrIpbOzYsfjoo49U7d5gIPm6KSkpqmWxN9euXVPXbE5BRKFCfu+mr9qA1BmLVJUGp2iDCnaj2zQP5PSIcjBbrFh1q2va4dNXva7qSkgRHaVTndN87chGkaFIg14pESY5u9LMQdr67tq1C1evXlWpA/Xq1VM1eqtXr45gIrnIknqRW7tiR31eTyXNiIiCje2mCSmTZiBj5z63cV2NKiqdQVeOL+ApeEjXNOmYNn9TEq7eTPeY23t34yoY1KEO6lcprao6yCY3GScKaBtix4rogw8+mOs569evR7t27RBoko4hZcuOHj3q8bhUmbh48aL6uHv37kU8OyKi/MnYdwimb6bBfu1G1qBGg5h7OyGmXzdofKhfTlQUpGuarOqu2pmziYQoJl3TWtVC/7aJKFsizq0UGZEnQfmTIakEnTt3htlsDvRUVMtiaZQhlSfOnj2LChUquB2X1WpHBYcuXboEaJZERLmzWyxInbsM6cvWZL7/e4umRDEYx4xAVIPaAZ0fkZCuaev3nsasdYew+1jmglJ20jVNqjB0aVYNMQxwKR+K/KdF2gvv27cPN2/ehNVq9Rjwzp071+OxQGjYsCF69eql0jHkMmbMGLfjy5cvV9eyOU9qERMRBRvruYswTZgC63H37lVRdzRQTSm08caAzY1I3EgxY/GWo5i74TDOX3UpmefaNa3era5ptdk1jYI86JXc3ZEjR6qKDb5ssCiKH2jXEmPSRtgb2VgnFSakPbJr0CsB+sSJE1W6xptvvlno8yUiyg/5XWreuB0pk+cC6S7vnOn1iB3WB9Gd2jB4oIBKPn8dc9YfxrJtxzw2mXB0TRvQrjaqlPFeL58oqILeP/7xj86WvpIKIFURoqJydkKRFd4zZ844a98WFnkcWVF2kHJqssFONttlV7duXRXcSh7ySy+9pAJcaUQhtYclmF+wYAHKly9fqPMlIsoPe2oaTD/ORsaW393GtZXKI37cKOgqu6dqERUVm82ObYfOYda6g9h60PMm8QqljRjErmnkZxq7o15YIStRooRKafjHP/6BZ5991mPA67qJrWPHjoWW4vD+++/j5ZdfzlF/VwLejz/+2GuN4HXr1mH8+PH47bffVIe2vn37qhJs2fN8feUovMzmFETkT5ak4yqdwXbJvUtkdOe2iB3SGxoDgwgqeqnpFizbfkyt7J644LKR0kWzWmUzu6Y1qMSuaRGuUSHESEUW9EpgKEGso62vL7m0e/fuRThj0EtE/mS32ZC2+BekLVghy2nOcY0xDnGjh8JwR8OAzo8i07krJpWru3jzUdxMy2zm5CpKr0XXZtUwsF0d1tSlQo2Riiy9YcCAASqFwFfhHvASEfmT7fJVmL6ZCstB9/KK+nqJMI4ZDm2pEgGbG0UeWU+T6gtScmzdnlOweVheK10sBv3bJKJP60SUjI8OxDQpwhRZ0PvWW29hyZIlaiObL/VsExMTkZSUVCRzIyIKZebtu5Hy/UzYU1w6R2q1iB3QA9E97oYmSDpcUmR0TVv9+wlVcsxb17R6VUqpFIa7m1RVq7xEYRf0litXDitXrlT5si1btkTJkt7fwti8eTOOHTtWVFMjIgpJdrMZKdMWwPzrZrdxbdkE1VlNX7NqwOZGkeWKdE3blIT5G5NwxUvXtI6NK2Nw+7poUK00q4ZQ+NfpTU5OVpvHpBqCt5q20pAit5a/REQEWE6chunrybCddd8nYWjTHHGjBkATw7eLqfAdOnWra9rvJ5BhzVl1qVhslEpf6NcmEeVKZnVNIwrroHfKlCl44IEHVJ6PXByte73hq0Aiopzk92f6ynVInbUYsLhUuImJhvH+QTC0ahbI6VEEdU2TYHeXl65p1ctldk3reie7plHwKLKfxNdff13V3pUqDtLgoUaNGjAYDDnOk5Xg6dOn+9TEgogoktiu34Rp0nRYdh9wG9fVrKrSGXRlEwI2N4qMrmlLth5VJcc8dU0Tretndk1rzq5pFMlB74kTJ6DVarF69WrUqVMn13NHjRqF4sWLF9XUiIiCXsbegzB9Mw326zezBjUaxPTsjJh+3aDR6QI5PYrgrmkxBh3ubVETA6VrWll2TaPgVWRBb6tWrXDkyJE8A14RFxeH7777rkjmRUQUzOwWC1Ln/Iz0Zb+6jWtKFofxkRGIqpcYsLlReKfRbD10DrPXHsKWg2c9nlOhVBwGtKuDnnfVQHxsznduiSI26JUuaNJl7fLlyyhdunSe58uqMBFRJLOevaA6q1mTT7mNRzVriLgHh0AbbwzY3Cg8pZotWL79uMrX9dY1rWnNshjcoQ7asGsahZgi68gmvv32W+zatQsffPBBrudJBQdp85uenrPsSThhRzYi8kR+LZvXb0XKlHmA2aWDVVQU4ob3haFjK+ZLkt+7ps3bkIRFm4947pqm0+KeZtVUvm5tdk2jIhDSHdkkXUFWbyXF4e2330bVqp7rR0qgu2DBArWhjYgo0thSUpHyw2xkbNvpNq6rUgHGsaOgq1Q+YHOj8Htxtef4JdVIYt3uk167pkm5sT6ta6FUfEwgpkkUekGvBLqHDx9WH8+bNy/P/4hcxSCiSGM5fEylM0hLYVfRXdohdnAvaKKiAjY3CrOuaTtPYva6gzh0ynPXtLq3uqZ1Ytc0CiNFFvT+5S9/weOPP66C2YSEBMTGxnrM201NTcX58+eLalpERAFnt1qRtugXpC1cIa/6neOaeCOMDw9DVJP6AZ0fRU7XtA6NKqtgt1H1BC4+UdgpspxeSVuoXr06fvnlFzRo0CDXc6Vdcffu3WG15iyNEk6Y00tE1ktXkDJxqlrldaVvUAfGMcOgLcHyjXR7Dp+6olIYcuua1rtVLfRvW5td0yhohHROb3R0NJ5//nmUKVMmz3O7dOmCHj16FMm8iIgCxbxtJ1K+nwV7alrWoE6H2IH3IrpbB2hYxYYKyGqzY/3eU5ld04567ppWrVyxW13TqiOWXdMoAhRp9Yb8iIS8Xq70EkUme7oZKVPnwbxuq9u4tlwZGMeNhL56lYDNjULbzVQzFm85irnrD+Ocl65prepVUMFuizrlw/7vLIWuRqG80ptfRqMRKSme/8MSEYUqS/IpmL6eDNs599U3Q7sWiBvRH5qY6IDNjUKX1NSVVd1l248hzey5a1qPFjUwsF0dVGXXNIpQfg16L1y4gLVr16Jnz55qo5qr5ORkn+5DNrJNnTo17Gv0ElFksdtsSF+xDqmzlwCu+xViomF8YDAMLe8I5PQolLumrTuELQc8d00rX1K6ptVGr5Y12TWNIp5fg95OnTrhwIEDGDRoEGbMmOF2rGnTprhxw3N3FyKicGa7dgOmb6fDsveg27gusbpqJawrk3eXSqLsXdPmrD+E5POe/642qVlGpTC0k65pOuaGE/k96LXZbOqVp6c04VGjRuGLL77w+b6YZ0RE4SBj9wEV8Npv3Mwa1GgQ07sLYvp0gUanC+T0KIScv5qCuRsOY/HmI7iR6rlrWuc7qmKwdE2rXCogcySKmKB39erVWLdunUpvyO4Pf/gDpkyZgsmTJ6NmzZqqmoOnwFbyeOW8v//97/6cGhFRkbJnWJA6e7FKaXClKVUCxrEjEVWnZsDmRiHaNW3PKdg8tE0rFR+Nfm1qo690TSvGrmlEQVG94f7778ePP/6Y53kypbi4OJXfG85YvYEoPFnPnFeb1awnz7iNRzVvgrgHBkFrZC1Uyl2GxYZVO0+ofN1Dp654PKdO5ZIY1L4uOjWtAoOe7xhQeGkUCtUbdu507xdfv359GAyZyfOfffaZT/chK8Amk8nfUyMiKlTygt28dgtSps4HMlzefjZEIW5EPxjat2TqFuXqys00LNh4RHVOu3zDpX7zLVoN0L5xFZXCwK5pRAEOeps1a6b+Ew4dOhTDhg1D3bp1nceKF/e9s5CnFsVERMHKZkpByg+zkLF9t9u4rmpFGMeNgq5CuYDNjYLf4dNX1aruLzuSPXZNi49xdE1LRPlSxoDMkSjUFUqd3r/97W/MySWiiJFx8AhME6fCfuWa23h01w6IHdQTmqigLYlOAe6atmHvacxad9Br1zSpqStVGLo1Z9c0ottVKP+DnnvuucK4WyKioGK3WpG2cAXSFv0iuQ3OcU2xeBgfHoaoxvUCOj8K7q5p8zYcxtkrnpswtXR0TatdHlrJaSCi4At64+PjUarU7ZdK+d///ocnnnjCL3MiIvI368XLME2YAusR98Y7+kZ1VcCrLc6uV+TupHRNW38IS7d56ZoWpUN31TWtNqqV8z0dkIgCFPTq9Xq/1PuV1WIGvUQUjMxbfofph1lAmkvnSJ0OsYN7IrpLe2i4J4FcNjduu9U1bXMeXdN63lUTxeLYNY2osARlgtD3338Ps9kc6GkQEbmxp6UjZco8mDdscxvXli8L47iR0FerHLC5UXBJk65pvx1Xwa7Xrmk1ymBQB3ZNIwrZoPfatWuoVatWgV8RX716FdevX/f3tIiIbovl2EmVzmA7777hyNChJeKG94Mmmit0lNk1TXJ1F+XRNU3ydeuwaxpRaAe9ErgeO3bstu+HtQeJKBjYbTakL/sVqXOXAtasPExNXAziHhgCQ4smAZ0fBZ783dubfAmz1h7CWi9d00rGR6N/m0T0bZ3IrmlE4RL06nQ6dOjQocC/OK5cuYK9e/eqvF4iokCyXbsO0zfTYNl32G1cX7uGaiWsLV0yYHOj4OiatnpXZte0gyc9d02rXamkaiTR6Y6q7JpGFI7VG3755Zfbuo85c+ZgyJAhfpsTEVF+mXfuQ8qkGbDfdOkOqdEgpm9XxPS6BxodA5hI7pq2cNMRzN+YS9e0RpVVCkPjGmWc71zKCnC6xYpovY5lyIjCIej1R1rCgAEDnK2LiYiKkj0jA6kzFyP9l/Vu49qEkjA+MlKt8lJkSrrVNW3l78lqlddT17RerWpiQNvabl3T5HYz1x7Er7tOIi3DqkqTdWxSBUM61EViJb5bQBSyQW96uksJn9sInL/77ju/zIeIyFfW0+dg+noyrKfcS0tF3dUUcfcPgjYuNmBzo8B1Tdu477QKdn8/ciH3rml3VkdstPuf1ZU7kvHetM3qfhwk8F22/bg69uLwVujSrFqhfx1EVAhBb2pqKi5duoSEhITbup9hw4b5bU5ERHntJzCv2YSU6QskUTPrQLQBcSP7w9C2BTfXRhhTWgaWbDmKOdI17bJLiouLu+qWx+D2ddGijueuabLCmz3gdSXjcrx6ueJc8SUK1Tq9L7zwAiZOnFgYd01E5Fe2myakfD8TGTv2uo3rqlVWtXd15csGbG5U9E5evIE56w6rrmmpZpcXQNm6pkkzCQlWcyMpDd4CXgc5LufJii8RhWDQ++2332Lz5s0YNGiQ6qxWsiRfwRJR8Mk4kATTxKmwX3WvDR7d/W7EDuwBjR86TFJorPRvP3z+Vte0M7B7iFPLSde0trXRq6VvXdNk05rk8PpCznt+aEtubiMqZH7/jc5SY0QU7OxWK9LmL0faklUS8TjHNcXjYRwzHFEN6wZ0flR0XdNW/JaMOesP4dg5z02RpPqClBxr1zB/XdOkSoPk7vo0jwyrOj/WwBdZRIWJ/8OIKKJYL1xSndWsR0+4jUc1qY+4h4ZCWzw+YHOjonHhmnRNS8JC6ZqWkrPlvV6nQeem1dTmtLpVCtY1TcqSSSqEL4GvnCfnE1HhYtBLRBEjfdNvSPlpDpDmUmVGr0fskF6IvqcdN6uFeQrDvuTLmLXuIH7d7b1rWr/WiejbJhGlb7NrmqQqSFkyqdKQFzmPqQ1EhY9BLxGFPXtqGlKmzIV5429u49qK5WAcNwr6KhUDNjcqXFJP99fdJzFr7UEcyKVrmqzqdvZz1zSpwytlyXLbzKbTatR5RFT4GPQSUVizHD2hau/aLl52Gzfc3Rpxw/pAw0Y4YenqzXQs3JSEebl0TWvXMLNrWpOaWV3T/EnKkElVBm9lyyTgleMsV0ZUNBj0ElFYsttsSF+6Bqlzl8oOW+e4Ji4WcQ8OgaF544DOjwrHkTOZXdNW7PDcNc0oXdNaZnZNq1A6q2taYZHGE1LajB3ZiAJPY5dEJwqIRo0aqes9e/YEeipEYcV25RpM30yD5UCS27i+bk3VSlhbqkTA5kb+J6uom251TdvhpWtalTLxalW3e/MaObqmFRXJI5YqDbJpjTm8REUfI3Gll4jCinnHHqR8NxN2U0rWoFaLmH7dEdOzEzRa38tOUfB3Tft561HMWX8YZ3LpmibB7l11KgQ80JTHZ1kyosDh/z4iCgt2cwZSZyxE+uqNbuPaMqVhHDsS+lrVAjY38n/XtLnrD+Pnrd67pnVrXh0D29fJs2saEUUOBr1EFPKsp87ipmxWO33ObdzQqhni7hsITeztlZ+iwJNMvN+SzmP22kPY5KVrWtkSsao9cK+WtVDch65pRBRZGPQSUUgHQumrNqoVXlhcVvyiDYgbNQCGNs1ZezfEpWdYseK34ypf11vXtEbVEzC4Q120z2fXNCKKLAx6iSgk2W6akDJpBjJ27nMb19WootIZdOXKBGxu5L+uaYs2H8H1QuqaRkSRhUEvEYWcjH2HYfpmKuzXbmQNajSIubcTYvp1g0bPX22ham/yJbWqK+W9rF66pvWVrmmtayGheGxA5khEoYl/GYgoZNgtFqTOW6bq77omdWpKFINxzAhENagd0PlRwVisNqzZdVIFu/tPuDcRcUismNk17R7pmhblv65pRBQ5GPQSUUiwnr8I09dTYD1+0m08qmkDxI0eCm184TcaIP+6ZkrHgk1HMH/jYVy67rlrWtuGlTG4ELumEVHkYNBLREG/Wc28cTtSJs8F0l1yO/V6xA7tg+jObRgMhZijZ69h1tpDWLnjOMxB0DWNiCIDg14iClr21DSk/DQH5s073Ma1lcojftwo6CpXCNjcqABd0/afyeyalnTea9c0qa3bI4Bd04gofPG3ChEFJcuRZJgmTIHtonuOZ3Tntogd0hsaQ1TA5kb57Zp2DHPWH/LaNa1FncyuaS3rBr5rGhGFLwa9RBRU7DYb0pasQtr85YAt661vjTEOcQ8NgaFZZj92Cm6nL91U7YGlTXBKes6uadG3uqYNalcH1cuzaxoRFT4GvUQUNGxXrsE0cQosB4+6jevrJcI4Zji0pUoEbG7kW/71jqQLmLXuoEpl8NY1rX/b2ujdil3TiKhoMegloqBg/m03Ur6bCXtKatagVovYAT0Q3eNuaLTstBXqXdMkhaFDo8rsmkZEAcGgl4gCym42I2X6QpjXbHIb15ZNUJ3V9DWrBmxulLuL11Ixb+NhLNzkvWtap6ZVVbBbr0rpgMyRiMiBQS8RBYzlxOnMzWpn3HfzG9o0R9yoAdDERAdsbuTdvltd09Z465pmjEaf1rXQr00iu6YRUdBg0EtEAcn9TP9lPVJnLpJ2XFkHYqJhvH8QDK2aBXJ65KVr2q+7M7um7Uv23DWtVsUSGNy+LrumEVFQYtBLREXKdv0mTJOmw7L7gNu4rmZVlc6gK5sQsLlRTtdN6Vi4+QjmbUjCxesu+da3SF+Qdg0qYVCHumjKrmlEFMQY9BJRkcnYexCmb6bDfv1G1qBGg5ienRHTrxs0Oq4OBlPXNFnVlQ1qnrqmxUXr0atlLQxoVxsV2TWNiEIAg14iKnR2iwWpc5Yifdkat3FNyeIwPjICUfUSAzY3ymJzdE1bfwi/HfbcNa2ydE1rVwc9WlRHXDQbhBBR6GDQS0SFynruAkxfT4E1+ZTbeFSzhoh7cAi08VwlDLSUdEfXtMOqqYQnzeuUx2B2TSOiEMagl4gKbbOaef02pEydB6S7lLOK0iNueD8YOrZi/meAnbnVNW0Ju6YRUQRg0EtEfmdLSUXKD7ORsW2n27iucgUYx42CrlL5gM0t0smLkd+PXMCstYewcf9pdk0joojBoJeI/Mpy+JhqJWy7dNVtPLpLO8QO7gVNFPNAA9U1beWOZMxZdwhHzl7zeE7D6gkqhaF9o8rQs2saEYUZBr1E5Bd2qxVpi39B2oIVspzoHNfEG2F8eBiimtQP6PwilZQZm78hCQs3J+GaKWfXNJ1Wg86OrmlV2TWNiMIXg14ium22y1dVZzVZ5XWlb1AbxjHDoS3BfNCitv/EZVVybPXOEx67ppUwGtC3dSL6tklEGXZNI6IIwKCXiG6LedsupPwwE/aUtKxBrRaxA+9FdPeO0Gj5NnlRdk1bu/sUZq076L1rWoUSGNShDrrcUY1d04goojDoJaICsaebVWUG87qtbuPacmVUZzV9jSoBm1skdk1btOUo5q4/nHvXtPZ10LRWWVbNIKKIxKC3gFq3bo3NmzfnGG/evDm2bdsWkDkRFRVL8ilVe9d27oLbuKFdC8SN6A9NTHTA5hZJjp1zdE1LVhvVPHdNq4kBbWujYkJ8QOZIRBQsGPQWwLJlyzwGvOKVV14p8vkQFRW7zYb0leuQOnuJvJeedSAmGsYHBsPQ8o5ATi9iuqZtPnBGBbvbvXRNq5QQj0HtaqPHXTXYNY2I6BYGvQXw9ttv49tvv1WrvdnVq1cvIHMiKmy26zdg+nY6LHsOuo3ralVT6Qy6Mtz5X9hd05ZuPYbZuXVNq11OpTC0qleRXdOIiLJh0JtP69atw6VLl/DQQw8xL44iRsbuAyrgtd9wCbY0GsT0vgcxfbpCo+OGqEB1TTPotapr2sB2dVCzQomAzJGIKBQw6C3AKu/QoUNVVyMGvRTu7BkWpM5ejPQV69zGNaVKwPjICETVrRWwuYUz+f2yU7qmrTuEDfs8d02TMmMD2tVG75Y1UdwY7Ze0iXSLFdF6HVeJiSgsaezy25V8sn37drRo0UJ9HB8fj759++Lxxx/H3XffXaD7a9Sokbres2ePX+dJ5A/Ws+dh+noyrCfOuI1HNW+MuAcGQ2uMK/J8YmRYgCh92JZBM9/qmjY7t65p1RJUCkOHxv7pmpZ0+ipmrj2IX3edRFqGFTFROnRsUgVDOtRFYqWSt33/RETBEiMx6M0HWeGdOXOmx/Gvv/4aJUrk761FBr0UjORXgnndFqRMnS9RWNYBQxTiRvSDoX3LIn2Xw3LiNNKXr4V5+67M+RiiYGjeBNHdOkBftRLCgZQZW7AxCQs2ee+a1ulW17T6fuyaJgH2e9M2e2xeIY/54vBW6NKsmt8ej4jIVwx6A0i+TT/99BMuX76Mw4cPY+XKldi9e7fbk7N69WokJCT4fJ8MeinY2EwpSPlhNjIkwHShq1oRxnGjoKtQrkjnY968A6Zvpsl77zkParWq25uhVTOEqgO3uqatyqVrWp9WiejX1v9d02SF98lPlnt8XNfA99M/deOKLxEVucKIkZjT6yNZ2br//vvdxubOnYtnnnkGx44dU0/KAw88gMWLF+frfpOSkpxPrCcMiKmoZBw6qloJ26+4v60e3bU9Ygf1giaqaH9dyAqv14BX2GzquLZiuZBa8ZWuaev2nMKstYewN/mSx3NkQ5qs6soqa3QhdU2TlIbcAl4hx+U8WfElIvK33OIfiY8SExP9+ngMem/DgAED0KlTJ3XZuXMnlixZglWrVqFz586BnhqRz+xWK9IWrkDaol/kLQ3nuKZYPIwPD0NU48CU4ZOUBq8Br4PUDV6xFvqHhyNUuqbN23AYF6557prWpn4lDO5QB3cUctc02bQmOby+kPOeH9qSm9uIKOQx6L1NJUuWxNKlS9WrFSllNmvWrHwFvfIqhqu5FCjWi5dhmjgV1qTjbuP6hnVhHDMM2uLFAjIv2bSmcnh9YN62C3EPDQ3azW3SNW3OusNY/ttxr13Tet5VU1VikKYSRUGqNMimNV/IeXJ+rIF/LojIv3KLf3JbBS4o/hbzg/Lly+PFF1/ESy+9pJbjiUKBecvvSPlxNuypaVmDOh1iB/dEdJf2gQ0ipUqD6ya63Mh5cn60AcFCVlK3HDyrSo5tP3TO4zmVEoyqtm6PFjVgjCnarmlSlkyqNPgS+Mp5cj4RUahj0OsngwYNUkFvVBRbflJws6elI2XKPJg3bHMb15YvC+O4kdBXq4yAk/xhQ5Rvga+cV8T5xrl2Tdt2HHPWH8Kpi567pt0pXdPa1UGr+hXVRrFAkFQFKUu2bLv7Cr8nch5TG4goHATHX4owUL16dXXNNsQUzCzHTqrNarbzF93GDe3vQtyI/tAEyWqprDJLWTLzxu15nmto0STgqQ1nLpswd/0hLN7ivWta1zurq81pwdI1TerwSsmyvKo3yHlEROGAQa+fXLlyRV0PHx78G2oo8kiObPqyX5E6dylgzXpLWxMXoxpNGFo0RbCROrxSsizXzWxaLaK7dkAgu6bNvtU1zVPsmFA8Bv3b1kafVrVQwg9d0/xJypBJVYa86vSyXBkRhQsGvX6yaNEiDBkyxNmxjShY2K5dh+mb6bDsO+Q2rq9dA8axI6EtHZxBjZQhkzq8edXpLepyZapr2u+3uqad8dw1TRpISBWGjo2r+KVrWmGRkmjVyxVnRzYiighsTuGjM2fOYM2aNejevTtKl3bviHTu3DmMGjUKM2bMyHEsN2xOQYUtY9d+mL6dDvtNU9agRoOYvl0R0+seaHTBv0FJdWRbsVZVaXB2ZGvRRK3wFmXAe+l6KuZvTMLCTUdw1ZTucWX07iZVVApDg2q+N6kJps13UqVBNq0xh5eIAo3NKQLotddeU62GpVLDu+++ixEjRsBgMKhyZVKmbMqUKfkKeIkKkz0jA6kzFyP9l/Vu47KqK6u7ssobKiSwlTq8UpZMVWmI0hdpDu+Bk5cxe+0hrN51AhZrzjWC4nEG9GldC/3b1EaZEv7tmlaUJNBlWTIiCmdc6fVRcnIynn32WbXae/36dVSqVAlt27bFyJEj0b9//wLdJ1d6qTBYT5+D6evJsJ466zYe1aIp4h4YBG1c6AZmRcVqtWGtdE1bdwh7j3vumlajfHG1qisb1AqraxoRUaRqxJXewKlWrRpmzpwZ6GkQeSWvX81rNiFl+oLMFVGHaIOqzGBo16JQu3yFg+spZizafCTXrmmt61fE4PZ10SyxcLumERGRfzHoJQoDtpsmpHw/Cxk73F8R66pVVrV3deXLBmxuoeD4ueuYvf4Qlm/33jXtXuma1rY2Kpcpmq5pRETkXwx6iUJcxoEk1UrYfvW623h0946IHXgvNHr+N7+drmkD2tbBvXcVfdc0IiLyL/41JApRdqsVafOXI23JKsltcI5riserUl5RDdlUwJPUdAt+3nZMNZM46aVrWrPEcipft3UAu6YREZF/MeglCkHWC5dUZzXr0RNu4/rG9WAcPQza4nwL3nPXtMNYsvUoTGk52xtHSde0ZtI1rTZqVWR9WiKicMOglyjEpG/6DSk/zQHSXGrF6nWIHdIb0fe04+aq7F3Tjl7M7Jq295TXrmn92mR2TSsZH1xd04iIyH8Y9BKFCHtaOlImz4F5429u49oKZWEcN6rIO5MFM+ma9svvJzBr3cFcu6ZJCoN0TZNVXiIiCm8MeolCgOXoCZXOYLvgXjPW0LEV4ob3hcZgCNjcgkleXdO0Ll3TGoZg1zQiIio4Br1EQcxusyF96Rqkzl0q5Qac45q4WMQ9OASG5o0DOr9gobqmrTuE1Ts9d00rJl3TWtVC/7aJKFsiLiBzJCKiwGLQSxSkbFeuwfTNNFgOJLmN6+vWhHHMCNVSOJI5uqZJsLsnj65pXZpVQwxb7BIRRTT+FSAKQubf9yJl0gzYTSlZg1otYvp1Q0zPztBotRHdNW3xliOqEoPXrmn1Kqpg987a5bixj4iIFAa9REHEbs5A6oyFSF+90W1cm1AKxrEjoU+sjkiVV9e0WIN0TauBAe1qo0qZYgGZIxERBS8GvURBwnrqLG5+PRm20+7dwaJa3gHj/YOgiY1BJHZN23roLGatPYRtXrqmVShtxKB2tVWbYHZNIyIibxj0EgVBLdn0VRvVCi8slqwD0QbEjRoAQ5vmEfcWvXRNW7r9GOasy6VrWq2ymV3TGlRi1zQiIsoTg16iALLdNKnc3Yyd+9zGddWrwDhuJHTlyiCSnJWuaRsOY/GW3LqmVcPAdnWQWCmyN/IREVH+MOglCpCMfYdh+mYq7NduZA1qNIjucTdi+3eHRq+PmJXuXbe6pq330jWtdLEY9G+TiD6tE9k1jYiICiQy/qoSBRG71arq7kr9XdizIjxNiWKqFFlUg9qIBGaLFauka9raQ0g6c9XjOfWqlMLgDnXZNY2IiG4bg16iImQ9f1F1VrMeO+k2HtW0AeIeGgJtsXiEu8s30lTXtAWbknD1pueuaR0bV8bg9nXRoFrpiMtnJiKiwsGgl6iI3sI3b9yOlMlzgXRz1gG9HrFD+yC6c5uwD+4OnryiUhhW7Uz23DUtNkqlL/Rrk4hyJdk1jYiI/ItBL1Ehs6emIeWnOTBv3uE2rq1UHvHjRkFXuQLCuWvaur2nVbC7+9hFj+dUL5fZNa3rneyaRkREhYd/YYgKkeVIskpnsF287DYe3amNWuHVGMKzruwN1TXtqKrEcP6qS1c5F63rZ3ZNa86uaUREVAQY9BIVArvNhrQlq5A2f7l0WHCOa4xxKnfX0KwRwtHx89dVbV3pmpbmoWtajEGHnnfVZNc0IiIqcgx6ifzMduUaTBOnwnLwiNu4vl4tVZ1BW6oEwrFrmqQwbD3opWtaqThVW7dnS3ZNIyKiwGDQS+RH5t92I+W7mbCnpGYNarWIHdBD1d/VaLVh1TVtmXRNW38YJy641Bp2ccetrmlt2DWNiIgCjEEvkR/YzWakTF8I85pNbuPaMqVhHDcK+ppVES7OXcnsmrZos/euaV3uqKaCXXZNIyKiYMGgl+g2WU6egenrybCdOe82bmh9J+JGDYAmNgbhUHJNqi9II4ncuqZJubE+rWuhVHzof81ERBReGPQS3UYgmP7LeqTOXAxYLFkHYqIRN2ogotvciXDomrZauqatO4TDpz13TasrXdPa18HdTaqyaxoREQUtBr1EBWC7cRMpk2YgY9d+t3FdjaowjhsJXdkEhHrXtAW3uqZdyaVrmqQwNKyWwJJjREQU9Bj0EuVTxt6DMH0zHfbrLpu3NBrE3NsJMf27Q6PTIVQdOnWra9rvJ5BhzSq15to1rXerWujftja7phERUUhh0EvkI7vFgtQ5S5G+bI3buKZkcRgfGYGoeokI1a5p6/eeVikM3rqmVStXTK3qdruzesC6pklptHSLFdF6nVppJiIiyg8GvUQ+sJ67ANPXU2BNPuU2HnVHQ9VsQhtvRCh2TVuy9agqOeata1qrehUwuEPdgHZNSzp9FTPXHsSvu06qhhcxUTp0bFIFQzrUZXUIIiLyGYNeojw2q5nXb0PK1HlAujnrQJQeccP6wnB365DLZ02WrmnrD2PZtmNeu6bd26ImBkrXtLKB7Zq2ckcy3pu2GVaXchEy52Xbj6tjLw5vhS7NqgV0jkREkcRmtyEtIx0xUdHQakJr8zKDXiIvbCmpSPlxNjK27nQb11WuoGrv6iqVd2s7jAyLCoaDsQGFpAZsO3RO5etuOXjWa9e0AdI17a4aiI81INBkhTd7wOtKxuV49XLFueJLRFSIiz/X0q5jU/I2/LRjJjaf2A6zNQMx+mj0qt8NY1vejwbl6iIUMOgl8sBy+BhME6fAdsm9TFf0Pe0QO6QXNFGZrXQtJ04jfflamLfvAswZgCEKhuZNEN2tA/RVKyHQUs0WLN9+XAW73rqmNa2Z2TWtbcPg6pomKQ3eAl4HOS7nyYovERH5HsjeNJtw0XQJF+Ry8xIupmR+fMl0OXPMdAkXTZdxyXQJGTaXspy3pFnSMXv3QszfuwTv93kD/Rv2RLBj0EvkQlZs0xatRNqCFfJbwTmuiTcibvRQGJo2cI6ZN++A6ZtpsoyadQfmDJg3blfHjGOGw9CqGQLVNW3ehiQs2nwENz11TdNpcU+zzK5ptYNwlVRWpiWH1xdy3vNDW3JzGxFFvBRz6q2A9aIzeJXAVT6/mHIZF2/eCmZTLiPdkrMcZUFYbFa8sHA86pSpFfQrvgx6iW6xXb4K04QpapXXlb5BbRXAaksUd47JCm+OgNftzmzquLZiuSJb8ZVX7nuOX8KstQexbk9od02TKg2e8o09kfPk/NgAVZUgouARyvmm3khwmhm4Zgaz8vFFl2DWEdjKWEpGakDmaLFZMXHLT3i/z+sIZvwrQSQLtNt2IeWHmbCnpGUNarWIHXgvort3zJGnKykNXgNeB5sN6SvWQv/wcBRF17TZ6w/h0CnvXdNkVbdTiHRNk7JkUqXBl8BXzpPziShy7Tt/EBO2/IglB1YgNSMNsVEx6Fmva9Dmm0pOrKzEqtVXlV6QuRLrGsA6Atwb6TeLZE6lYkugjDEBZY1lUMZYGmWNCSgdVwofrf0CGdac7xhmt/jAcvyz92tB/WKDQS9FNHu6GSnT5sO8dovbuLZcAoxjR0Ffo0rO29hsmTm8PgbTcQ8NLZTNbVeka9qmJMzf6L1rWodGmV3TGlUPra5pquNbkyqqSkNe5DymNhBFrnl7l6i312W10UEC36LON7XarLiccsUlpUCuMwNXR0qBI73gato1FIXi0cWcAWyZW5fMj0u7BbgJcaURpdN7TJd4f/UnPj2WfM9llT3OEItgxaCXIpYl+ZSqvWs7d8Ft3NC2BeJG9ocmJtrzDaVKg2xa84WcJ+dH+68awuFTV1Qjidy6pvVqVQsDQrxrmtThlbJkuW1mk413ch4RRe4Kb/aA15/5ppIucSX12q10gqxgNjOgdR+TgNeO3Dff+oPREJcZwMZlBrNl4xOQEJcVzGaOlUFCXClE6738HfORpInIqrkEtHmR8+T8YMaglyKOrNSmr1yH1NlLAIvLL8qYaBgfGAxDyztyv4MovarS4FPgK+fJ+bdJAr/1e09h9tpD2JVH17Sud1YPi/xWKUMmVRm8lS2TgFeOs1wZUeSSlAZvAa+3fFPZ/3A9/UZWxQJ17UgzyBxzBLNysdp9219wOyQ4VcGqMQEJjpVZCWTjE7ICXPk4LqFIV1K1Gq1KE5FV87z0qtctqFMbROj/ZSTKB9v1GzB9Ox2WPQfdxnW1qsE4diR0ZUrneR+SqiBlyaRKQ14MLZrcVmrDzVQzFm85irnrD+Ocl65pLaVrWvs6aFGnfEilMPhCGk9IHV52ZCMiT6uwksPrizl7FuHQxaTMElwpl1VObWGL0uqRcCuQdU0pyFyNLeM2Fm8wBu3v77Et71dpIrm9uNBrdXik5X0Idgx6KWJk7D6gAl77DZdNARoNYnrfg5g+XaHR+b4ZSurwSlmyXDezabWI7tqhQHOVmrpSW3fZ9mNIM3vumtajRQ0MbFcHVQPcNa2oVnylLJlUaZBNa8zhJYoM8ra6s5asqiGbVT/27I3zPr3t7giQd53dd9vz0Wl0SDCWUiuvZVw2fGWt0GYFsyViigdtIJsfkhYiedHe0kgk4JXjwbhhMDsGvRT27BkWlcoglRRcaUqVgPGREYiqWyvf9yllyKSMmdeyZVqtOp6fcmXylttWR9e0A567ppUvKV3TaqNXy5pB0TWtKEmgGw5pGxS5wrGcVkGkW8xqtTV7Y4TsObNS3UAaKBQ2DTQoFVfyVkpBabcc2cyUgtK3cmRLqwoHOm3kVYvp37CnyouWNBGp0uCokCEpDbLCGwoBr9DY5S8tBUSjRo3U9Z49ewI9lbBlPXsBpq8nw3ritPsBrQZRzZsgpmfn26qjqzqyrVirqjQ4O7K1aKJWeH29X0fXtDnrDyH5vOeuaU1qllH5uu0aVIJOF7l/LIlCUaiV0yqIDKvFpXLBrYuU4FLVCtwbI0g+bVGrUqISutfplKOCgazMlo4rCb2WL6iD7cVbYcRIDHoDiEFv4ZEfa/O6LUiZOj/3DWe3VmRvt3OabI5TVRqi9D7n8J6/moK5Gw5j8eYjuJGaS9e0drVRu3Kp25ofEQVPOa3sbwsHa/tWKcHlqFzgCGbd68hmBbhyXlFULpDcV9egVavRqZVHCcS8ke/znNHfh80LjEjRqBBiJL60obBjM6Ug5YfZyPCllq6fOqepQNeHsmTOrmnrDmV2TfNQlaBUfDT6tamNvtI1rVjwdk0josCW0yoI+R0kNWKz6si6B7OuebOXUq7kGkz6S4w+GuVU+oBLOoEEtbeqFbhu+JJV8uzm7b07LPJNqfAx6KWwIi2Eb349GfYr+Sj8XQSd0zIsNqzaeULl6x46dcXjOXUq3+qa1rQKDOwwRhSR5bQKGshK7quzq5ejFJcEsrfKcTmCWgloM2wWFDaDLsqlUoFcl3Epu5WZI+soxSV1Z29nw1e45JtS4WN6QwAxvcF/7FYr0hatRNrClfIXIP93YIhCyY/f8HvntCs307Bg4xHVOe3yjZy7jKUIQfvGVVTJsVDrmkZE3skKadMP7/a5qP/OZ9fkyI+UblieV2Pdc2Tl43RLzq6M/iarprKZyzWYlZxYZ11Zl+oFxaLjA/L7jJsFw0cjpjcQ5WS9dAWmCVNgTXJvWaurnwjr/qSAdE47fPqqWtX9ZUeyx65p8TFR6N2qFvq3TUT5Uka/PCYRBQ8JvHwtpyXnvbrkH7iSdtUt7SAlI7VIKheUjiuVszGCM282K5gtGVs86ANJmV8wt8GlwGLQSyHNvOV3pPw4G/ZUlz8uOh1iB/WE4Z62uPbsG0XWOU26hm3Yexqz1h3ErqOeu6ZJTV1JYejWPDy6plH44ArZ7VUukBJcrquy529eVCujeaU3OEzZOduvc5LSWpIPK3mxjlJcrjmyjjQDKdXFygUUKfiTTiHJnpaOlKnzYF6/zW1cW74sjONGQl+tsvq8KDqnSde0RZuPqEoM56+meu2aJsFui9rl2ViBgkoklNMqaOUCKcGVVXbrVi1Z1xzZW7mzUrmgKEjKQNbqq2uKgctYXGm1Wis5tUTkjkEvBaXcSoBZkk/B9NVk2M67r6Ya2t+FuBH9oXFJUSjMzmknL9zAt8t249fdnqswSLvc7qprWm1UK1c83/dPFIhyWhL4zt69ULUdDeZyWgVdzb6aej0zJ/ZW2S33bl9ZjRGupF4tksoFokKxcqhcvKIzjcBRqcA9zaA0ovXRRTIfonDFoJeCimr2sHwtzNtdmj00b6KCV13lCqrKQursnwFr1h9pTWwM4h4YDMNdTQu9c5rs+9x2q2vaZi9d09ScNMAT/ZqhV6v8d3sjitRyWgUh/ydvpN+8FaxmBrNZAexFt7qykoLga7rB7TDoDM4VWIvNgj3nDnisYRvsdXqJwg2DXgoasiKbIzg1Z6j0BDkmtXRtp9wDTV1idRjHjoQuwXvzBmk8Ibe9nc5padI17bfjmLPuMI6fv57n+VJA4uM521G3SmkkViqZ5/lE4VpOq6CBrMmcklmhwHU19lZqwUWX4FaOma0+5O3fJglQVdktx0avWzmymR+7t6uVBgqulQvkBQbLaREFHkuWBRBLlrmv8N5455Pc0xBcaTSI6dMVMb3vgUanK7TOadI1bd6Gwypn11PXtLx0b14dLw5vle/bEQV7Oa2CSMtIc2uC4MiRdZTdkqD20q36sr5WPrgd8jUl3KpckODWmjYrR9YxViLm9isXcLMgke9YsozClqQ0+BrwakuXVKu7+to1CqVzmrwO3Jt8SaUweMvX9dWvu07i+aEtuXmNQrqclpzvrQyUrLJe9NCaNnuOrIxJA4WiKMElFQmcXb2yBa+uY1LhQKctukYwLKdFFFgMeingZPVV5fD6QqtB/MtPQVesYLVtJYBNt1gRrdflCESla9qaXSdUi+CDJ711TSuJPq1q4aPZeVeEEGkZVvV4LE9GwURWGmUF15fAV/JTF+xfhsu3Ug1cc2Tl42tpeaf7+IOstLrVkXWkFNxKNchML0hQNWejdPz/RkQ58TcDBZ6kG/hSS1fY7NBKTd18Sjp9FTPXHlQrrxKISmWFjk2qYEiHukgoHouFm5Iwb2MuXdMaVcbgDnVV1zRJCPp8we/qfvIijyMBNlGgyVvrUlpLBas3L6JR+frYenJHnrczW8342+I3C2VOkvvqugrr2gjBdUy6gEXr/dM4hogiF4NeCjxpCiGBbCE1kVi5IxnvTdusmkc4SMC6bPtxLN9+XK34uh7Lq2ua7E+RgFlunxc5j6kNVFgkFUdWWl3b0zpSCrI2emW2rL1kugKrvfArF8Too2+V3XLZ3OUMYN2DWVltJiIqKgx6KeAkz1bycy17D/m9iYSs8GYPeF3JaPZjvnRNkxViCaa93a/QaTXqPKL8BrKS++qaSuBaV9Y1oL1kuoQMm6XQ5ySNDlybIagqBs6UgsyKBY40A6Mhzq1yARFRsGDQSwFlt1qRtmCFTwFvQZpISEpDboGpq5Z1b3VNq5N31zQpQyZVGbwF1BLwynGWKyOHFHOqqlagym65VCtQY87uXpnBbLolvUhKcEnagKQYZAbZl9VKcJRWjzsrN8XARr1xV5VmKsiVTmAMZIko1DHopQLLb/mv7KwXL8P09RRYjybnfXI+m0g4Nq2t3nnSp3MNeh3+/nCHfKUidGlWDdXLFfeaK8yAN/xJcOq6Auu6ycut01fKZVV3tigqF8hGL0dKgWPTV87WtWVQMta9BBfLaRFRuGPQS37tmuZrUKoaUfw4G0hzWdHS65DasgVOHD2L6mdPIho2pEOL4xWroHz/rijVvJ7Pbw//dvg8pv96AGaLbzmMcl5Bqiw4VnylLJm3qhAUWjKsFhWkOspseUopcASz0gmsKEhpLalMkNUMwXNjBCnVpdcW7Nc6y2kRUbhj0Et+7Zomq7HSAc0be1o6UibPVee70lYoi70d78EbK4/AaisLTckyMMAGM7Swp2mgm7ELL2pj1epqbl3TJM9W6useO5e/Mkq3W2VBAl2WJQteVpsVl1Ou5JpS4AhwpcJBUZCUgbzqyMrHUoJLcmqJiOj28K805WuFN0fA68pmU8el5a+nFV/LsRMwTZgC2/lLbuOGjq1wtkMHvPHlGmd+rF2jQTqyglAZl/xZSSfInjZw4Zp0TUvCQumalmIu0NfGKguhR96Ov5p63W2Tl1tKgbPr12UV8Mr5hS0uKjZbAJvZtjZrLOtYtD660OdDRERZGPSSf7um2WxIX7EW+oeHu+X+pi/7Falzfna7vSYuFnEPDlapETNzqbDgIMclf1bSCSSFYV/yZbWqu2b3SY9d00rGR6N9w8pYsvUoqyyECHleJWUgM3DNDGazAtiLbjmzkoJgsRV+CS5pzuDMi3WW4MoMZrPnykrlAiIiCk4MesnvXdPM23Yh7qGhanOb7ep1tfpr2X/Y7Rx9nZowPjJCtRSWgFU2gvliza6TqrqCBLsHvHRNq12ppKrC0PmOqmqDWtNaZVllIcCBrGzikhVXR2OErDQDGXMPZqWtbVFULlBlt5zdvTznyEopLqluwMoFREShj0Ev+b9rmpyXYYF5/2GkTJoBu8ll17pWi5i+XRHT6x5nxQfZAOZLdzN1boYV707dnGNcMhPaNZSuaXXQuEYZtyCFVRYKR1pGmnsThJTMHFlHSoEEtZdu5c760u72dslGrIS4UreqFnjPkZUxaWnLCgVERJGFQS/5JkoPe1QUNBl5B772KD1SZi6CefVGt3FtQikYx46EPrG627hsIJMg1NfA15VRdU2rif5taqNC6ayuadmxyoJvZJXVkRub1dHLUb3AfUxquwaSBK2PtX4I/RrcqwJaqXCg07LlMxERecagNwLcbj1dIbc7lFARdc/mXVPXbNdAky3gjWp5B4z3D4ImNmfbUQk+fW3r69o1bWC72ujevAZio33/MY7EKgsWmwWXU6665chmL8Xl+Fha2hYFWWl1pBQ4cmQdKQWOVdmrqdcwdsafvebtysa0rzd/j74NeqjVWyIiotxE1l//COOPeroOknc7Kb043sjjh0ayZqMtLqvB0QbEjRoAQ5vmueZF+tLWVzSsXhoPdGnkU9e0cCYBnwSyqp6sS45sZiOEW12/bpXiupJyFXb1zBQuyX11b1XrqFrgvulLuoBF6w153t/zC8fnuVFNjk/c8hPe7/O6H78SIiIKRwx6w5TUzL35zTRoPNTTTd+8A/F51NPNTlICDiEanxpr4knTUa8/OK5hqK56FZXOoCtfJtf7PnnxhqqwoNVoYPUSnEm8PK5XUwy/27cGFUXNH92sZMOXrLQ6VmNlZdZZdst5nRngXjJdUS1jC1uMPvpW0OoezGZeuwezsVE5V/Fv5/u55MAKn85dfGA5/tn7NeboEhFRrhj0hukKb46A14WMy/HiXurpeuLIu12H0jipjUGf9PNoY76CGNhUmJp9zdXQ427EDegBjV7vvWta0nnMXnsImw6cgd3LQqRep0HHxlUwolN9v2w483er1X3nD2LClh9VgCabtSTw61mvK8a2vB8NytVVX6fkvuaVI+vo9JVhs6CwSaMD1yA2QaoVGN0rFshqrIxJCa5AVC6Q58jXzW9ynpzPbmJERJQbBr1h6OL8lYjOo56uBL5yXoUnHvDpPl3zbo/r4/CFrjquaPQYkH4OrqGjjK1v1gYPDOnttfrCit+O59o1Taov9G9bC41rlkQpYzz0Ol2hB6f5lWJOxdSds/GPlR+7rbjKfc/evVBdSseWREpGKtIsLq2WC4lOIyW4HPmwnisWONIMikcXC/oSXPKiRJ4jXwJfOU/OJyIiyg2D3nw6ePAgxo8fj1WrVsFisaBLly54++23Ubt2bQTLpjXtrn0+nSvnyfm+bm5z5N2WyUjDn01HUMfqUooMwLaoEvgivgbe7d8hx23z6pomK7qdm1bDHQ21WJ48H0+teNUvwamYt3cJXsiWH+oITufvXYL3+7yB/g17It2S7rapy+Omr1u1ZqXubF4up17F7dBAo1rQOhojJLgEs1lBbWYwWzI2vEpwydciz7s8R3npVa9bWH3tRERUOBj05sPixYsxdOhQDB8+HPv374der8czzzyD5s2b4+eff0bbtm0DPUXY0jMQlW3zjyQgpGtsiLZroXVJRJDz5HxdrG+rZJJe8I9mJVBu5XbE2W3O+9XYdfgxtiqWx5bDiyNau6Uh7E2+pFIYcuua1rd1Ivq1ScTaE6vwtA/BaX7sOrNPbYiyetkQJY/17PxX8H9L3larskVBSmtJA4SsZgiOKgbuaQal4kpCr43c/6LyQkee99w2s0mTiUda3lek8yIiotCksUvSIeXp8OHDaNasmVrR3bZtG3S33nKX1d5GjRrh8uXL2LdvH8qUyX3Tliu5ndizZ4/f5pmSZsa1p19XubYHDTfwY4lkrDCeQ5rWhhibFl1N5XH/tWqoay6GNGhR4uPXEReT9056e2oaUn6aozbIZb9fnV2PxGIt8GynR9CjcQtkWGz4dfdJlcKw/8Rlj/eXWLGkaiTRuWlVGKJ0Kv1g4KQH8wxw5oz+HnXLJOJKqpTgclQpuLUK61KxwLE6eyX1GgIpSheFTwf+E+Xiy6qAVlZuJaeWCr5K7/rzUJAXQkREFPwKI0aK3GWkfHr22WdhMpnw5JNPOgNeIau9jz32GJ5//nm8+OKLmDhxYkDnGWOIwoqYUrip24Px5fbAqsl6TSMB6sJiZ7Ak/izeON8IRmtj9DXkHYBZjibD9PUU2C5exhLj2Rz3a9VYcPDmJjy1eCsG7vsDjh8qi0vX07x2TZMWwU1qundNm7D5B5/KUw35fgwyrBlqQ1phi4uK9ZoT++aKf/l0HzLXttVacpNVAUlAW6dMLVWWTKo0OFJeJKVBVngLmvJCRESRhyu9Pjhy5AgSExOdH9esWdPtuLwKady4MQwGA06fPo2EhISAvYoRb38xCZOufOIWmGans2swutSf8H9/GO31HMn3Tft5NdLmLZNCvWqF98HKm3O9X9i1qHxjNKKt5VRtWJsmHYaYdLSoXwwNasUhQ3PzVhmui1l5szcv4rzpIgLtw75/R+USFZ25slK5wBMJuJt+eLfPm6x2PruGOad+4O/KG0REFLy40hsgixYtUtdGozFHwCvq16+P6OhopKenY/bs2Rg3bhwCKbn0fliv5v5aRgLXE6UPeD1uu3INpm+mwnLgiHNMUhpyDXiFxoYz8VOgtUfBqjPBjszV20PHAMilEMjb3J5yZBNiS+O91f9RrXV9CU77NuzhUzDFTVaBId9HrpgTEVFBMej1wdKlS9V11apVPR6XdIcqVaogKSkJmzdvDmjQK6th60786tO5a0+sUedLMCGrlo4KBWd378DpNatwyXITl8qk45LOjEu6dOyO8a1FrU2bChsKZ1NYlFaPd3u/pnJkM2vKJqiWtt4Cyz3n9xdKcMpNVkRERKGFQa8PJKVBSGDrjaQ0SNB74ID31dNgLOrf9cvBuJxyRTVQcFO8cOYntWuz14+Vj5cfWo1tp37P8/b9GvbEwEaeawAXZXAquaSyiSqvTVbMOSUiIgoODHp9cOHCBXVdrFgxr+dIeoO4cuVKvu5bAmVH3oon+c1lyU9Rf5F89ST8TQK+Zzr8MUeDBKlcEKXz/CPXsWYbn6o3BFNwyk1WREREBZdb/CPxkWM/lb8w6PWBlCMTcXGeNzYJ260OaGlpvgWbhSU/+aa5ibPpEI1iSLXGQ2c3QmczIk1/Emb9+Txv279hLzzedkxEBKeZ834d/+z9GjdZERERBTEGvT6QqgxSjzc3GRmZm6VKly6dr/uWVzH+rt7gy1v60u3r3sS7USM5BSVP30AZqwGlrQYkWKNxXVsK/4upjXO6aJRyuY016jyS4ydJTYZCyWMN5eCUm6yIiIjyJ7f4J7dV4IJi0OsDyddNSUlBaqr3zVnXrmU2QchPc4rC4suq6bvNHsc9y07Bfk3q9JZQ4xLKzo2pgGkxlWB1qaFb0hiNPq1roV+bvlh3MrFQ81gZnBIREVFhYNDrgwYNGuDEiRM4d+6c13MuXbqkrj2VNAsEr6umdbvgvus1UGPGfrgWH7usicJ/jTWwJyprB1utiiUwuH1d3HNHZte0osxjZXBKRERE/sTmFD5466238Nprr6FSpUo4depUjuPSqS0+Pl59PHfuXPTv3z+gzSm8FfWPunwDN7+aApw87XZ8S1QJfBZXAze1esgCb7sGlTCoQ100zdY1zdv9Mo+ViIiI/InNKQJkyJAhKuiVbmtnz55FhQoV3I7v2rXLWcGhS5cuCDaSv3tlxVZEzV+CKGtWWoIZGkyKq4plhjKIi4nCkJa1MKBdbVQsbfTpfrkaS0RERKGCQa8PGjZsiF69emHx4sXqMmaMe2WC5cuXq+vRo0c7V3yDxZbfj8E8ZS4aXj7jNp6sjcHH8bVgK18WT7argx4tqiMuWvJ7iYiIiMIP0xt8dPDgQTRr1gx33HEHNmzY4ByXzW2yBH/jxg3s3r0b5cuX9/k+Czu9wW6z4dDz76Ksyb2T2uLostjTpDn6d6yPlnUrQKv1nsJAREREVNQKI0ZiIqaP6tati4kTJ2Lr1q146aWXkJ6erlIdRowYgatXr6pc3vwEvEVBo9VC06WD8/MbGj1W3Nke7f5vHN5+9B60rl+RAS8RERFFBKY35MPIkSNRtWpVjB8/Xm1qkw5tffv2xd69e3Pk+QaL2n06Yt3GnSivtaLcH+/D0EplAz0lIiIioiLH9IYAKqrqDXZzBqDXqZVfIiIiomDH6g1UIBoDN6gRERFRZOPSHxERERGFPQa9RERERBT2GPQSERERUdhj0EtEREREYY9BLxERERGFPQa9RERERBT2GPQSERERUdhj0BtkhZgdxZgptPC5C1187kIXn7vQxecudDUK4eeOQS8RERERhT0GvUREREQU9hj0EhEREVHYY9BLRERERGGPQS8RERERhT0GvUREREQU9jR2u90e6ElEqmLFiiEjIwOJiYnq86SkJHXt+JxCB5+70MXnLnTxuQtdfO5CV1IRPXfyOFFRUbhx44bf7pNBbwBVqFABJpMJ1apVC/RUiIiIiIJGcnIyjEYjzp4967f7ZNBLRERERGGPOb1EREREFPYY9BIRERFR2GPQS0RERERhj0EvEREREYU9Br1EREREFPYY9BIRERFR2GPQS0RERERhj0EvEREREYU9Br1EREREFPYY9AaJgwcPYtSoUahYsSLKli2LESNG4PDhw4GeFuVDamoq/vOf/6Bq1ao4duxYoKdDedi9ezdGjhyJcuXKITo6GrVq1cJTTz2FM2fOBHpqlIejR4/igQcecD53derUwfjx45GWlhboqVE+/eMf/4BGo8G3334b6KmQj1q3bq2es+yXFi1aINgx6A0Cixcvxp133omYmBjs379fBUzFixdH8+bNsWHDhkBPj/Jw48YNvPfee6hZsyaefvppnDx5MtBTojwsWbIELVu2xNSpU3HhwgWYzWYVSH3yySe444478NtvvwV6ipTLAoE8d9OmTYPNZlPPnSwQvPnmmxgwYECgp0f5sGPHDrz++uuBngblw7Jly7B582aPx1555RUEOwa9ASa/rIcNG6ZWKr7++muUKFECRqMRn332mVr17d+/Py5evBjoaVIuli5dio4dO2Lw4MGBngr54NKlS+pdlTZt2mDevHkqiFq1ahX69OmjjksQLMFTSkpKoKdK2WRkZKjfl2+88YZ6sSm/G+X5a9eunfP/4qxZswI9TfKBrMrLar08pxQ63n77bbUqv2/fvhyXgQMHItgx6A2wZ599FiaTCU8++SR0Op1zXK/X47HHHlO/1F988cWAzpFyN2TIELRt2xZjx44N9FTIB1999ZX65bxy5Ur069dPveDs1KkT5s+fj/vvv1+dc+LECUyePDnQU6VsvvnmG7zzzjvq96WkNQh5/mbPno2SJUuqz9esWRPgWZIv/va3v6nnrlq1aoGeCvlo3bp1atHgoYceQv369XNcJMUh2DHoDaAjR45gwYIF6uNu3brlON6zZ091/eOPP6ofNApuskpPwW/r1q349NNPc/yCls8//PBD9YJTbN++PUAzJG8efvhh54q8K8ntlZV74Qh+KXjJC05JLZIXoBRaq7xDhw6F3W5HqGLQG0CLFi1S15LOIPmg2ckrJ1nNkJw1Wcmg4OYIlii4vfvuu4iLi/N4TDaRNmrUSH0sOfYUXAwGg9dj8nzJu2X33Xdfkc6J8ufq1asYM2YMJkyYgDJlygR6OuQjWQSQ/UeSgy0LPJIiForvqjDoDSDJPxOy298T+QVepUoV9bG3xHEiyp/atWv7FFg1bty4iGZEt0vyQmXTr1RwqFu3bqCnQ7mQ1JS+ffuiV69egZ4K5YOkFTncvHkTU6ZMUWlhkmN/7do1hAouTQU4vUE4AltPEhISkJSUhAMHDhThzIgik7xtJ/8v5R2WUNiUQZneeustDB8+HK+++mqgp0K5kIob27ZtY+pQCP5eHDRokApyZfO9pKdIyUcxY8YMtYlt9erVKl4Jdgx6A0h2iYtixYp5PcexWePKlStFNi+iSCWrhZI//8wzz6BUqVKBng7l4ezZsyrQlco3TZs2xS+//IJ77rkn0NMiD06fPo0///nPah+Lt/QiCk4ajca5yddh7ty56veklFjds2ePqsQh6Q/BjukNAXT58mV1ndsvAKlDKVh0najwffTRR6pUoLxNTsH/XN17770q4BU7d+5Ejx49MGfOnEBPjTysFD7yyCP405/+hLvuuivQ0yE/kLKOUs9cXmw6ap9L6cdgx6A3SDdlODhqGJYuXboIZkQU2eV4pMar1KBkBYDgJ6tMv//+u0r/khJKwmKxqFKPUgaSgodUS5G6ylKmjMJHyZIl1d4kR1pDKNTIZtAbQI4fFGlf640jQZy7XIkKjwRJ48aNU129ZLWQQoe0j540aZJ67hxpY6HwNmukkP0oknP9/fffu9Wip/BQvnx5Zy8BeQEa7Bj0BlCDBg3U9blz57ye46jP66mkGRH5hwS87du3x8svvxzoqVAB/fWvf0WFChXcNglT4P3rX/9Sf8fkbfD4+Pgcl+TkZHXeH//4R/W5o2QghY5Bgwap66ioKAQ7Br0B1KFDB3V99OhRr6tPjhbE3bt3L9K5EUWK1157TeXMf/HFF4GeCt0G+YPrKIMVGxsb6OmQS4qe1WpVf888XRyNDtLT051jFFqqV6+uruvVq4dgx6A3wO1rHbtaZRdydrt27XJWcOjSpUuRz48o3Emgu3btWtVymG+9hj7HSm/Lli0DPRW6RXLkJbD1dnEETNJiWj6XagAUWq7cqi4lZQODHYPeAGrYsKFzZcJTDtry5cvV9ejRo9XbPkTkP/JHVi5SesdT9zXZFCU1KCl0HDx4EHfccYezJTERFU132SFDhqBFixYIdgx6g6DsjrwV9+WXX7qNy+a2iRMnqg1sjg0aFNxcy8rJW3UUvGTj03vvvac218hzJWlEcpFNUCdOnMCyZcvQs2dPPo9B5vr16143y0glB9lJLu1tich/zpw5g6lTpzrLrLqSPUnyezR7DBOsGPQGmLTMlOB269ateOmll9QfWUl1GDFihOpRLqtQsjuSgpvkrMlz5TB79my1UkjB57PPPsOYMWOwf/9+9f+vbNmyzku5cuVQrVo1VcFh06ZNzg0aFDz7IKSNtKzkLly4UOWLSi3z+fPnq3fE5B2zUFhtIgq1fQ8jR45U705LuoosysnfPPn/9sorr6iWxKFSVpVBbxCQHyYp6iztGStVqqR+of9/e3cCXNP1xwH8xC72fS1Bqdoj9l1RVFGUEhPLUCmmtIouxvavXQ1qp0Wrto7WUoLalVRbS1satbRKqX3flzj/+Z6Zc+e9l7fcJDfy3n3fz8yT5O733Rf53XPP+f3whzcuLk7UqVMntQ+PfJg8ebJ6PI7R4xryUaIFf/bs2al6bOQMNyP9+vUzBs94g4CXlaP8y+DBg9VgGSTFx/UJCwsTrVu3VtkakGcZGTiIyFrDhw8X7du3V4FudHS0Cn6joqLUTeeCBQtUY0GgCJFm/vcnIiIiIgpgbOklIiIiIttj0EtEREREtsegl4iIiIhsj0EvEREREdkeg14iIiIisj0GvURERERkewx6iYiIiMj2GPQSERERke0x6CUiIiIi22PQS0RERES2x6CXiIiIiGyPQS8RERER2R6DXiIiIiKyPQa9RERERGR7DHqJiIiIyPYY9BIRERGR7THoJSIiIiLbS5faB0BERETWkFKKdevWiRUrVogiRYqIihUriu7du6f2YRH5BQa9RERENnD16lURFRUl4uPjxRdffCEKFiyY2odE5FfYvYGIiCjAXb9+XTRq1Eg8ffpUbNiwgQEvkRts6SUiIgpwkZGR4vz582LHjh0iXTr+aSdyhy29REREAQx9eDdt2iQGDBgg8ubNm9qHQ+S3GPQSEREFsPnz56uvGLjWt29fUatWLdG0aVOxevXq1D40Ir/CoJeIiCiAszVs27ZN5MuXTxQvXlzMmTNHbN++XeTKlUu0b99eLFq0KLUPkchvhEj8xhAREVHAuXz5ssifP79o2bKliImJMaZfunRJPPfccyJnzpzi4sWLqXqMRP6CLb1EREQBSg9ay549u9N0BMK1a9dWwS9eRMSgl4iCBFq7Ro0apVI57dy5UwSDe/fuiUmTJomIiAiRJUsWUbp0afUeYDpZ6+bNm2LLli3PfL/oxoC+vGfOnEkwD591BMVo7YXffvtNnD59+pkfI5G/YNBLRGLZsmXqj2NISIjxyp07t/jrr7+8rvf666+LtGnTOq2XJk0a8dprrwl/8eTJE9GrVy9RokQJMXr06KB51Iu8rY0bNxYZM2YUBw4cUAEZzh3vQdu2bVP78GyXPQFFIcqUKeN2/v3798Wnn36quhv8888/prZ5/Phx0aVLF1GoUCHVX/eNN94QJ0+edLtsjx49xC+//CLOnj3rNP3vv/8WLVq0EBkyZFA/P//882LkyJHqRogoKKFPLxHRw4cP5cqVK2WhQoXQz1+9KlasKO/cueN1vXv37sk333xTLd+tWzd5+/Zt6W9u3bold+/ebZzXjh07pN316NFD5smTR8bHxxvTli5dqs4/c+bM8vHjx6l6fHYxceJE2axZM7fvJz53mF+gQAHjs3fq1Cmf24yJiZGhoaHqGt64cUP9Dvbu3Vtmy5ZNxsbGJlgev3OVKlWSbdq0kY8ePVLTli1bJnPmzCn//PNPp2WfPn0qe/bsKSMjI9XvPFEwYdBLRE5OnDihgiX9R7pTp04+1zl48KBa9tChQ9Jf3b17N2iC3kuXLsk0adLIKlWqJJi3efNmt4FTchw7dkyOHj1aBpt58+bJYsWKyevXr7udv2rVKvVe9+3b13TQi9+/LFmyyMqVK8snT54Y0xFUlylTRubNm1devnw5wXqYhiAZ17xOnTqydevW8o8//nC7D2yrRo0a6riIggnLthCREzwCrVChgti1a5f6+euvvxZVq1YV77//vsd1dJ9B/dUfZcqUSQQLpKxCOdrQ0NAE815++WXL9zdx4kT16D6Y/P777+Ltt98Ws2bN8vi579Chg/qK7gVIJWbGu+++K+7evSv69++vug5p6H7Up08fMXjwYDF06FCxcOFCp/VQlMJsejJsa8GCBSI8PFzUrFlTdO/e3dR6RIGOfXqJyK0aNWqoProwbNgw8f3333tcVi/nz9DXOFicO3dOfX0W5Wg3b94cdLlg8ZQU/WgxiMxMwJgjRw5T20Uf3PXr16vvUVzCFfrnwtKlS8XVq1dFclSqVEm0atVKBe7J3RZRoAievwJElCjI+4lBTxAfH68G1eCPMvm/27dvP5ObEbQoYzBjsKV7X7VqlTh06JD6HUmfPr3P5c3efOg8u8i0gYGXrsqWLasGJj569MiSamsYHIfPypQpU5K9LaJAwKCXiDwaPny46NSpk/r+2rVrol27dqbTXaEVzDGrw+LFi50ehzvOa9SokdO6CKIQACD7QM+ePdW08+fPi27duql8pBjR/tFHH6lgHBAEjBs3TnXNwCP9evXqiV9//dXr8eERMh4VFy5cWGTNmlW1rP34448el8f5o4tHuXLlVFCSJ08e9fgaaaBcIZCYPXu2qFixokoR9vjxYzFgwAD1GBznqoNSX/So/+rVq6u8q1i/bt26Yt68eca5O3rhhRfUMuPHj1c/79mzR/2MF6p1mfHw4UMxffp0Ua1aNZXyCu9N+fLlxYgRI9TxaFOnTlWfDT1twoQJxr5cW34Tex4pdf2RrguBHq4djqFBgwbi448/TvD582Xy5MnqK4JeK+mnKZ66iqC7Q9GiRdX3P//8c7L3V79+ffV1xowZ4saNG8neHpHfS+1OxUTkfxo2bChHjhxpDAALDw83BuJ07tw5wfIYnONukM7GjRtlSEiImrdo0SKneVeuXJHVqlVT87A/bf369bJChQrG/rp3764GyuXOnVsWLVpUZsiQwZg3YsQIlT2iSZMmavBP4cKFjXn58uVTI98d6Xlr1qxR54Rjw4AvPT19+vRy9erVCc7vyJEjsnTp0nLBggVqf3hPxo8fr9bPlCmTGhymjRs3ThYpUsTYJt7HXr16yezZsxvTMMDJl3PnzsmIiAhZq1YtefToUWOQU+PGjY33zFOmDOzT9X0148GDB2p/eI83bdpk7BODqrC9jh07JlgH+9DnacV5pNT1x+cN6w8YMEBlQ8AgMWT0KFeunMyRI4fp9+jkyZPGPjCAzwz9++FrIFv58uXVMk2bNvW4DAagYZkGDRpIK+hBqwsXLrRke0T+jEEvEXkNeuHMmTNOaZcmTZpkKugFHYi4Br3wwQcfJAjOEJxcvXpVBZmYV7duXZVe6fDhw0YaKB1oYSR7165d5cyZM1XABrt27VLBK+bPmTPHaX/6+MuWLStnzZqlgp/79+/L2bNnG8EUgtPz588b62CZUqVKqWVcIUUb1sF7g+BLw/f58+dX8xDsTZs2TQXKgwYNkq1atZLXrl3z+v4jrRTOESmqLly44DQP28HxYNvt2rWzNOjFNcJ6GP3vCDcCmI4bBMfz9BX0JuU8Uur6T5gwQU3D+q5BLAJks2bMmKG2g8+L2bRvZoNe/ZnxdF2hfv36RjpBK1SvXt10lhaiQMfuDUTkEx63fvvtt0aS+w8//FBs3brV1Lre+jyif6IrPHpGYYzKlSurnzNnziyWLFmiMkpAtmzZ1KNtuHLliujdu7ca6a63hUfWeLwNR48edbtf9FXu16+f6qaArA59+/Y1RtffunVLfP7558ayGCWPalcoPuAKZV4BRR8cB/rhmEuWLGmc48CBA9Vjd/SdxEAldP3wBgOVkD0DA40KFCjgNA/b+d///qe+R79O9Ku1yp07d9RX132iLykgIwS6eZiVlPNIqet/6tQp9VVnJdFKlSol2rRpY/qcUOgDUAXN6oGC+r11l3VDwzWABw8eWLJP/Tn96aefLNkekT9j0EtEptSpU0fMnTtXfY++lJ07dzZdXSopEOzo4MI180JYWJjb7zX0+XQM4lyhX6kr9B3Vg4ccyxQjIENVN/Sl1H1W9WvIkCEq2MLLtXqdvkFIbH9R0EF3lSpV3M5v3769kYINAaFVEEAiyEefXg3Xev/+/cbPeC+exXlYff114I7PLfqwop+19tlnn5k+pxMnThjBt9X0Z8Ybfdy4MbCC3s6///6r+nMT2RmDXiIyDYEh8ogC0hwlZmCblSnGfAUHunXZ3SApTxwH1OmUX3D48GEVJGOgj+sLA9LQ4obXoEGDEmwvKdCSt3fvXvW9p/yvCBQxSA58DdhLDGwX1xgt+xg4hoGMuNnRxwNmMzUk9zysvv4oRY3WYgxgxKBCDHqbOXNmoltMUd45pYJetHKD44BBVzdv3jTy8lpBtyrjenEwG9kdg14iSvTIdV3gAIEKWgftQrcQOgZcCAQQ4D+rVrBLly4ZrXnegh99rFYfF25iEMA3bNhQ5XKNjY31WpjEX8/DFYLUffv2iffee0+1zKPLCnLUIhvHtm3bTG9HH6fuZmClF1980egu44nOqesupVlS6Bb1xN4kEgUiBr1ElChIm7Ry5UpRpkwZ9fPy5cttk+dTB2l4pO7YEoZH+kj/5Y1VAQP6Gbv2Q/XWQmdlJTQEqkgrhv7a6OPZsWNHp6pggXIe3o7pk08+UV0UUN0M54Zje+WVV8QPP/xgahtI4QYp8YRD90X29H6hlRr9mKFZs2aW7NOxpdufKyoSWYFBLxElGv44rlu3zqg0hcfEnujBPt5axvyluAH6NUKTJk2cct+Ct8AeLYiuZWGT0yKpbyi8Bdo61+9LL70krBIdHS3i4uLUQDFfg+38+TzcQf9r/fgeATZyBCOwxw0O8vyOHTvW1HaKFStmBKBW02WL//vvP3HhwoUE89HVBtBSbdX7pd9/XC9vA+iI7IBBLxElKQhFMLhixQqfLYG676O7P+IoFgAIOlIbWmoxgA3HixKzmh7Zv3HjRqPog7vuAM2bN7fsWFCEAQ4ePCiOHDnidhlkJkDf1sjISLfnkpSbCZ2Bwlt/WtebF72su76xyT0Pq1vx165d6zQtIiLCuGHTn0Vf0OVD9/u2+mYNXS10wQt83lzpjCkofaxbnJNLtxx7GmxIZCcMeonI7WAZPWDGmxYtWohJkyZ5XUaPmv/yyy+NAWJoyUIKMN1yhcwH6CvpGFDpfqDu+no6ZhBwFzDroM9TX1J3rXSooIbAHCP7HbM7ILUZKpMBWkDbtm0rNmzYoPozo5sHUmSFh4cbLYCux5iYFF/aO++8Y6SSQgulu5ZllITG8bjrFqAHW+mAxizduouKebrvKFpDHW8C0AVi9+7dRiU6fVOD1ly87ydPnlRVzpJ7Hilx/YcNG5bgeuhKdWaDPj3YEdvGYD8zHG8IfPVdnjZtmupnO3/+fKfp2B+eJmAAm071ZoVjx445da0gsrXUThRMRP4jPj5erl27VqZNm1YVDjh+/Lip9Xr06OEx8f6WLVuMqmwoGlCsWDGZMWNGuXjxYqOIAl4lS5aUU6dOVeucPXtWFipUSE3H8o7FInCMjutNmTJFTXOsAKaLHqAy2unTp415KLqA6ajuhQpUKLSAAgNz586VoaGhbgtQwL59+2SuXLmMfTq+UDzBtWADKoihQpg+hri4OJlYqPalC3tERUUZxR327t0rw8LCVJU3d8URcM10RTgUk4iJiTG9TxTP0OeFSnO4BiVKlJA7d+40pqN6Wdu2bY33HBXo9DwUV8Axo5hJcs4jJa7/wIED1bSqVauq6wk4lubNm6siF2Y/6zhWLI9tbd261efyqPymC2PghUp+vopaLF++XKZLl04OHTpUFd3A+bdu3Vp9BvG+WQXHhuuM4zp06JBl2yXyVwx6iUj55ptvVDDqGtRlzZpVXrx40eu6+MNcu3Ztj9WmlixZogKRzJkzy3r16skdO3ao6QheUHoVlbNQ+UwHUTpI1i/8HB0drQIYXW3L8YXqWAha3nrrrQTzsO7YsWOdytx26NBBFixYUJ0vgu2ePXsaJXI9QSDXu3dvFcBhfzgfHL9rwItqWu6CY1RvSyxUbkPgg+pkCE4QgLZs2VKuW7fO7fI1a9Z0u28E9EOGDPG5P5xL//79VVCHG4M+ffrI69evq3koB4zKaniPdfUzvQ6qomEf2P+BAweSdR4pdf110Ktf+Cwi6MY6jkG1GaNGjVLbGDNmjNflULkQwavrMWEaKgJ6s2fPHlVeGdehePHi6rok9jh9iY2NVceDcuBEwSAE/6R2azMREVGgQNU+FMVAijHHHMaBZsSIEaorynfffSdeffXV1D4cohTHoJeIiCiR0Pcb5aUxQA8D0AIN+kWjvz1eKI1NFAwY9BIRESVBly5d1IDPmJgYEWhQIhoDCDEgUxcIIbI7Zm8gIiJKYuCIbBDTp08XgeT48eNizJgxYtOmTQx4Kagw6CUiIkoCFHPYvn27StdmVXGSlIYUZSi/vGbNGpVqjyiYsHsDERFRMn311VciJCREdO3aVfir/fv3q2B38ODBLDlMQYlBLxERERHZHrs3EBEREZHtMeglIiIiIttj0EtEREREtsegl4iIiIhsj0EvEREREdkeg14iIiIisj0GvURERERkewx6iYiIiMj2GPQSERERke0x6CUiIiIi22PQS0RERES2x6CXiIiIiGyPQS8RERER2R6DXiIiIiKyPQa9RERERGR7DHqJiIiIyPYY9BIRERGR7THoJSIiIiLbY9BLRERERLbHoJeIiIiIbI9BLxERERHZHoNeIiIiIhJ293+NiAY+XHPwzQAAAABJRU5ErkJggg==", "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": "ovi", "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 }