{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# MVRegressor.py\n", "# G. Cowan / RHUL Physics / November 2021\n", "# Simple program to illustrate regression with scikit-learn\n", "\n", "import scipy as sp\n", "import numpy as np\n", "import matplotlib\n", "import matplotlib.pyplot as plt\n", "import matplotlib.ticker as ticker\n", "\n", "from sklearn import linear_model\n", "from sklearn.neural_network import MLPRegressor\n", "from sklearn.model_selection import train_test_split\n", "from sklearn import metrics" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# read the data in from file\n", "events = np.loadtxt('trainingData.txt')\n", "nEvt = events.shape[0]\n", "X = events[:,1:] # eta, s1, s2, s3\n", "y = events[:,0] # E\n", "X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.5, random_state=1)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Test R2 score: 0.959\n" ] } ], "source": [ "# create regressor object, train and test\n", "model = \"MLP\" # choose one of the options\n", "if model == \"linear\":\n", " regr = linear_model.LinearRegression()\n", " model_label = \"Linear regressor\"\n", "elif model == \"MLP\":\n", " arch = (10,20,20,10)\n", " regr = MLPRegressor(hidden_layer_sizes=arch, activation='relu', max_iter=2000, random_state=1)\n", " model_label = f'MLP regressor ' + str(arch)\n", "elif model == \"SignalSum\":\n", " regr = linear_model.LinearRegression()\n", " model_label = \"Sum of raw signals\"\n", " X_train = X_train[:,1:2] + X_train[:,2:3] + X_train[:,3:]\n", " X_test = X_test[:,1:2] + X_test[:,2:3] + X_test[:,3:]\n", "\n", "regr.fit(X_train, y_train)\n", "y_pred = regr.predict(X_test)\n", "R2 = regr.score(X_test, y_test)\n", "print(f\"Test R2 score: {R2:.3f}\")" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc8AAAGkCAYAAABXS66yAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB4eElEQVR4nO3dd1xT1/sH8M+FQMKeogzBLaKirXVXRRFxVkGtdW9btV9n3RPrrrb1Zx0dWtxad6vWutC2KopFW8W6RdwDNSAQZJzfH7cJuRmQGxLCeN6vFy/k3pN7Ty5tHs56DscYYyCEEEKIwawsXQFCCCGkpKHgSQghhIhEwZMQQggRiYInIYQQIhIFT0IIIUQkCp6EEEKISBQ8CSGEEJEoeBJCCCEiUfAkhBBCRLJo8ExNTcXkyZPRrl07lCtXDhzHYe7cuTrLxsfHo23btnB0dISrqysiIyNx584dnWVXrlyJwMBASKVSVK5cGVFRUcjKytIq9+zZMwwaNAienp6wt7dH06ZNcfz4cVO+RUIIIaWQRYNncnIyvvvuO2RmZqJbt256y127dg0hISF4+/YtfvrpJ6xfvx43btxAixYt8Pz5c0HZBQsWYOzYsYiMjMRvv/2GUaNGYeHChRg9erSgXGZmJkJDQ3H8+HGsWLEC+/fvR/ny5dG+fXucOnXKHG+XEEJIacEsKDc3l+Xm5jLGGHv+/DkDwObMmaNVrmfPnszT05PJ5XLVscTERGZjY8MmT56sOvbixQsmk8nYiBEjBK9fsGAB4ziOJSQkqI6tWrWKAWBnzpxRHcvKymJBQUGsUaNGpnqLhBBCSiGLtjw5jgPHcfmWyc7OxoEDB9C9e3c4OzurjgcEBKB169bYu3ev6tjhw4ehUCgwePBgwTUGDx4Mxhj27dunOrZ3717UrFkTTZs2VR2TSCTo168fzp8/j4cPHxby3RFCCCmtiv2Eodu3byMjIwPBwcFa54KDg3Hr1i0oFAoAwJUrVwAAdevWFZTz9vaGp6en6ryyrL5rAkBCQoLJ3gMhhJDSRWLpChQkOTkZAODu7q51zt3dHYwxvHr1Ct7e3khOToZUKoWDg4POssprKa+r75rq99WUmZmJzMxM1c+5ubl4+fIlPDw8CmxFE0IIMQ3GGFJTU+Hj4wMrq6JvBxb74KmUX2BSP2doObFllRYtWoSoqCi9ryOEEFJ07t+/Dz8/vyK/b7EPnh4eHgB0twRfvnwJjuPg6uqqKqtQKJCeng57e3utsg0aNBBcV981Ad0tXQCYNm0aJkyYoPpZLpfD398f9+/fF4zJEkIIMZ+UlBRUrFgRTk5OFrl/sQ+eVatWhZ2dHS5fvqx17vLly6hWrRpkMhmAvLHOy5cvo3HjxqpyT548wYsXL1CnTh3Vsbp16+q9JgBBWXVSqRRSqVTruLOzMwVPQggpYpYaLiv2E4YkEgm6dOmCPXv2IDU1VXU8KSkJMTExiIyMVB1r3749ZDIZoqOjBdeIjo4Gx3GCtaQRERG4du0azp07pzqWnZ2NzZs3o3HjxvDx8THbeyKEEFKyWbzl+euvvyItLU0VGK9evYpdu3YBADp27Ah7e3tERUWhYcOG6Ny5M6ZOnQqFQoHZs2fD09MTEydOVF3L3d0dM2fOxKxZs+Du7o527dohLi4Oc+fOxbBhwxAUFKQqO2TIEKxatQo9e/bE4sWL4eXlhdWrV+P69es4duxY0T4EQgghJYuF15mygIAABkDn1927d1XlLly4wEJDQ5m9vT1zdnZm3bp1Y7du3dJ5zRUrVrAaNWowW1tb5u/vz+bMmcPevn2rVe7JkydswIABzN3dnclkMtakSRN29OhRUfWXy+UMgCCBAyGEEPOy9GcvxxhjlgvdJV9KSgpcXFwgl8tpzJMQQoqIpT97i/2YJyGEEFLcUPAkhBBCRKLgSQghhIhEwZMQQggRiYInIYQQIhIFT0IIIUQkCp6EEEKISBQ8CSGEEJEoeBJCCCEiUfAkhBBCRKLgSQghhIhEwZMQQggRiYInIYQQIhIFT0IIIUQkCp6EEEKISBQ8S4no6GhwHAeO43Dy5Emt84wxVKtWDRzHISQkRHCO4zh8+umn+V4/JCREdX2O42BnZ4d69erh66+/Rm5urgnfCTFUaGgoPvnkE8GxmTNnonPnzvD19QXHcRg0aJDe19+5cweRkZFwdXWFo6MjwsLCEB8fb3R9Tpw4gSFDhiAwMBAODg7w9fVF165d8ddff+ksHx8fj7Zt28LR0RGurq6IjIzEnTt3jL7/nj170Lt3b1SrVg12dnaoVKkS+vbti5s3b+osf+zYMTRt2hT29vbw9PTEoEGD8OzZM6Pvn5CQgFGjRqFp06ZwcHDQ+/+i0vbt21G/fn3IZDL4+Phg3LhxePPmjaDMunXr4Ovri7S0NKPrRcyDgmcp4+TkhHXr1mkdP3XqFG7fvg0nJyejr12lShWcPXsWZ8+exY4dO+Dr64vx48dj2rRphakyMcL+/ftx+vRpzJo1S3D8q6++QnJyMj744APY2trqff3z58/RokUL3LhxA+vXr8dPP/0EhUKBkJAQXL9+3ag6rVmzBomJiRg7diwOHTqEFStW4NmzZ2jSpAlOnDghKHvt2jWEhITg7du3+Omnn7B+/XrcuHEDLVq0wPPnz426/5IlS5Ceno4ZM2bg8OHDmD9/Pi5evIh3330XCQkJgrKnTp1Chw4dUL58eezfvx8rVqzAsWPHEBoaiszMTKPuf+HCBezbtw/u7u4IDQ3Nt+yWLVvQu3dvNGzYEL/++ivmzJmD6OhoREZGCsoNHDgQDg4OWLp0qVF1ImbESKHI5XIGgMnlcovW48cff2QA2LBhw5idnZ1Wffr168eaNm3KateuzVq1aiU4B4CNHj063+u3atWK1a5dW3Ds7du3rEqVKsze3p69fftWVH3T0tJElS+Mt2/fsqysrCK7nymkp6fne75Ro0bso48+0jqek5Oj+reDgwMbOHCgztdPmjSJ2djYsMTERNUxuVzOPD092YcffmhUnZ8+fap1LDU1lZUvX56FhoYKjvfs2ZN5enoK/jtNTExkNjY2bPLkySa7/8OHD5mNjQ0bOnSo4HjDhg1ZUFCQ4L+L06dPMwBs9erVRt1f/dnv3LmTAWAxMTFa5bKzs5m3tzdr166d4PiWLVsYAHbo0CHB8WXLljEXF5ci/X+mJLD0Zy+1PEuZ3r17AwC2bdumOiaXy7F7924MGTLEpPeysbFBgwYNkJ6enm9rYe7cueA4DvHx8ejRowfc3NxQtWpVAHx38urVq1G/fn3Y2dnBzc0NPXr00Oq+Y4xh4cKFCAgIgEwmw3vvvYejR48iJCRE0A198uRJcByHTZs2YeLEifD19YVUKsWtW7cAQNW6cHZ2hr29PZo3b47jx48L7vX8+XOMGDECFStWhFQqRbly5dC8eXMcO3ZMVebixYvo3LkzvLy8IJVK4ePjg06dOuHBgweqMgqFAtOmTUPlypVha2sLX19fjB49Gq9fvxbcr1KlSujcuTP27NmDd955BzKZDFFRUXqf58WLF3H+/Hn0799f65yVlWH/S+/duxdt2rRBQECA6pizszMiIyPxyy+/IDs726DrqPPy8tI65ujoiKCgINy/f191LDs7GwcOHED37t3h7OysOh4QEIDWrVtj7969ou+t7/4+Pj7w8/MT3P/hw4eIi4tD//79IZFIVMebNWuGGjVqGH1/Q599bGwsHj9+jMGDBwuO9+zZE46Ojlr379u3L1JSUrB9+3aj6kXMg4JnKePs7IwePXpg/fr1qmPbtm2DlZUVevXqZfL73b59GxKJBG5ubgWWjYyMRLVq1bBz506sXbsWAPDxxx9j3LhxaNu2Lfbt24fVq1cjISEBzZo1w9OnT1WvnTFjBmbMmIH27dtj//79+OSTTzBs2DDcuHFD572mTZuGpKQkrF27Fr/88gu8vLywefNmtGvXDs7OztiwYQN++uknuLu7Izw8XBBA+/fvj3379mH27Nk4cuQIfvjhB7Rt2xbJyckAgLS0NISFheHp06dYtWoVjh49iq+//hr+/v5ITU0FwAf7bt26YdmyZejfvz8OHjyICRMmYMOGDWjTpo1W12B8fDwmTZqEMWPG4PDhw+jevbve53jgwAFYW1ujZcuWBT5zXTIyMnD79m0EBwdrnQsODkZGRkahxh7VyeVyxMfHo3bt2qpjt2/fRkZGht7737p1CwqFwiT3v3PnDu7duye4/5UrV1T30nV/5Xlz0Xd/GxsbBAYGat2/QoUKCAwMxMGDB81aLyKSRdq7pYiluw6UlN22cXFxLCYmhgFgV65cYYzxXVSDBg1ijLFCd9tmZWWxrKws9ujRIzZ16lQGgPXs2TPf186ZM4cBYLNnzxYcP3v2LAPAli9fLjh+//59Zmdnp+q+e/nyJZNKpaxXr146X6/+fpTvvWXLloKyaWlpzN3dnXXp0kVwPCcnh9WrV481atRIdczR0ZGNGzdO7/u5cOECA8D27dunt8zhw4cZALZ06VLB8R07djAA7LvvvlMdCwgIYNbW1uz69et6r6euQ4cOLDAwsMBy+rptHz58yACwRYsWaZ3bunUrA8DOnDljUF0K0rdvXyaRSNiFCxdUx5Tdo9u2bdMqv3DhQgaAPXr0qND3zsrKYiEhIczZ2ZklJSWpjiu7R8+ePav1mhEjRjBbW9tC3zu/btsFCxYwAOzx48da59q1a8dq1Kihdbxv376sfPnyha5XaWLpz15qeZZCrVq1QtWqVbF+/XpcvnwZcXFxJumyTUhIgI2NDWxsbODj44Ply5ejb9+++P777w16vWZr6sCBA+A4Dv369UN2drbqq0KFCqhXr55qpmJsbCwyMzPx4YcfCl7fpEkTVKpUyaB7nTlzBi9fvsTAgQMF98rNzUX79u0RFxenmtHYqFEjREdHY/78+YiNjUVWVpbgWtWqVYObmxumTJmCtWvX4urVq1r3V06Q0Zzt2rNnTzg4OGh1FQcHB6NGjRo634umR48e6eyiFIvjOKPOGWrWrFnYsmULvvrqKzRo0KBI788Yw9ChQ/HHH39g48aNqFixosH3MMV7N4SY+3t5eeHZs2dGdacT85AUXISUNBzHYfDgwfi///s/KBQK1KhRAy1atCj0datWrYrt27eD4zjIZDJUrlwZ9vb2Br/e29tb8PPTp0/BGEP58uV1lq9SpQoAqLpLdZXT91pd9wKAHj166K3fy5cv4eDggB07dmD+/Pn44YcfMGvWLDg6OiIiIgJLly5FhQoV4OLiglOnTmHBggWYPn06Xr16BW9vbwwfPhwzZ86EjY0NkpOTIZFIUK5cOcE9OI5DhQoVVO9JX33zk5GRofd9G8LNzQ0cx2nVAeCfAQC4u7sbfX0AiIqKwvz587FgwQKtZVAeHh4AoPf+HMfB1dXV6HszxjBs2DBs3rwZGzZsQNeuXUXdv7DvvSDq99f8Peq7v0wmA2MMCoUCjo6OZq0fMQwFz1Jq0KBBmD17NtauXYsFCxaY5JrKiTrG0vyL2tPTExzH4Y8//oBUKtUqrzym/LBRHwNVevLkic7Wp657AcDKlSvRpEkTnfVTfpB5enri66+/xtdff42kpCT8/PPPmDp1Kp49e4bDhw8DAOrWrYvt27eDMYZ//vkH0dHRmDdvHuzs7DB16lR4eHggOzsbz58/FwRQxhiePHmChg0b5lvf/Hh6eqqCnDHs7OxQrVo1XL58Wevc5cuXYWdnp/rDxRhRUVGYO3cu5s6di+nTp2udr1q1Kuzs7PTev1q1apDJZEbdWxk4f/zxR6xbtw79+vXTKlOnTh3VvTp27Kh1f+V5c6lbt67qXkFBQarj2dnZuHbtmmrSn7qXL19CKpVS4CxGqNu2lPL19cWkSZPQpUsXDBw40NLV0alz585gjOHhw4d47733tL6UHzKNGzeGVCrFjh07BK+PjY3FvXv3DLpX8+bN4erqiqtXr+q813vvvadzXaS/vz8+/fRTvQkEOI5DvXr18NVXX8HV1VVVRrnOb/PmzYLyu3fvRlpaWoHrAPMTGBhY6Ak9EREROHHihGAWampqKvbs2YMPPvhAMAtVjM8//xxz587FzJkzMWfOHJ1lJBIJunTpgj179qgmWAFAUlISYmJitNY6GooxhuHDh+PHH3/Et99+qzWbVcnX1xeNGjXC5s2bkZOTozoeGxuL69evG31/QzVu3Bje3t6Ijo4WHN+1axfevHmj8/537twRBFpiedTyLMUWL15scNnbt29j165dWseDgoLM9j9t8+bNMWLECAwePBgXLlxAy5Yt4eDggMePH+PPP/9E3bp1MXLkSLi7u2PChAlYtGgR3NzcEBERgQcPHiAqKgre3t4GLRFwdHTEypUrMXDgQLx8+RI9evSAl5cXnj9/jr///hvPnz/HmjVrIJfL0bp1a/Tp0weBgYFwcnJCXFwcDh8+rPpQO3DgAFavXo1u3bqhSpUqYIxhz549eP36NcLCwgAAYWFhCA8Px5QpU5CSkoLmzZvjn3/+wZw5c/DOO+/oXGZiqJCQEFVSAc1x0lOnTqmWDeXk5ODevXuq32urVq1UreDPPvsMmzZtQqdOnTBv3jxIpVIsXrwYCoUCc+fOFVxz0KBB2LBhA+7evat3jBkAli9fjtmzZ6N9+/bo1KkTYmNjBefVW/xRUVFo2LAhOnfujKlTp0KhUGD27Nnw9PTExIkTtd7vqVOnwBjL97mMGTMG69atw5AhQ1C3bl3B/aVSKd555x3Vz0uWLEFYWBh69uyJUaNG4dmzZ5g6dSrq1KmjFXSV7zkxMTHf+6enp+PQoUMAoLr3qVOn8OLFCzg4OKBDhw4AAGtrayxduhT9+/fHxx9/jN69e+PmzZuYPHkywsLC0L59e8F1c3Nzcf78eQwdOjTf+5MiZpFpSqWIpWd8KanPts2Pvtm2+r7mzJnDGNOdJMFQytm2z58/13l+/fr1rHHjxszBwYHZ2dmxqlWrsgEDBghmaObm5rL58+czPz8/Zmtry4KDg9mBAwdYvXr1WEREhKqccrbtzp07dd7r1KlTrFOnTszd3Z3Z2NgwX19f1qlTJ1V5hULBPvnkExYcHMycnZ2ZnZ0dq1mzJpszZ45qkfq1a9dY7969WdWqVZmdnR1zcXFhjRo1YtHR0YJ7ZWRksClTprCAgABmY2PDvL292ciRI9mrV68E5QICAlinTp0Mfp5yuZw5OjpqzeRljP896ftdas78vHXrFuvWrRtzdnZm9vb2LDQ0lP31119a1+zevTuzs7PTqreYe+v6qLlw4QILDQ1l9vb2zNnZmXXr1o3dunVLq1yDBg1YhQoV8n8ojH+O+u4dEBCgVf7IkSOsSZMmTCaTMXd3dzZgwACdiRY8PT1ZkyZNCrz/3bt3Rd1/69atLDg4mNna2rIKFSqwMWPGsNTUVK1yx48fZwB0/m7KMkt/9lLwLCRL/wLLsjt37jBbW1u2YMECS1elyH366aesVq1aLDc31+z3Kl++PPvss8/Mfh9dUlJSmEQiYd98841F7p+QkMAAsAMHDljk/ozx2cGaNWtmsfsXV5b+7KVuW1Ii/P3339i2bRuaNWsGZ2dnXL9+HUuXLoWzs3OZ7M6aOXMmNm7ciN27d+c7g7iwEhISkJ6ejilTppjtHvn5/fff4evri+HDh1vk/jExMWjatCk6depkkfvfvn0bO3bs0MoNTCyPY6yAgQSSr5SUFLi4uEAulwtSjRHTunXrFj755BP8/fffeP36NVxcXBASEoIFCxagZs2alq6eRRw4cACvXr0q1PgpKd5iYmJw8+ZNjBgxwtJVKXYs/dlLwbOQLP0LJISQssjSn720VIUQQggRiYInIYQQIhIFT0IIIUQkCp6EEEKISBQ8CSGEEJEoeBJCCCEiUfAkhBBCRKLgSQghhIhEwZMQQggRiYInIYQQIhIFT0IIIUQkCp6EEEKISBQ8CSGEEJEoeBJCCCEiUfAkhBBCRKLgSQghhIhEwZMQQggRiYInIYQQIpJE7AvS0tJw8uRJnD59Gg8fPkRGRgY8PT0RFBSE1q1bo3bt2uaoJyGEEFJsGBw8b968ieXLl2Pr1q148+YNOI6Dq6srZDIZXr16BYVCAY7jULduXYwZMwaDBg2ClRU1bIurzMxMfPLJJzh27BjkcjmCgoLw5ZdfolmzZpauGiGEFHsGRbfx48ejTp06OHPmDGbNmoVz584hMzMTycnJePjwIdLT0/Hw4UPs3LkTwcHBGDt2LIKDg3HhwgVz158YKTs7G5UrV8bp06fx+vVrjBw5Eh988AHS09MtXTVCCCn2OMYYK6hQq1at8Pnnn6Nly5YGXTQlJQVff/01XFxcMHbs2EJXsjhLSUmBi4sL5HI5nJ2dLV2dQnF3d0dMTAzq1atn6aoQQki+LP3Za1DL89SpUwYHTgBwdnbG7NmzS33gLI4OHToEjuNUXxKJBJUqVcKECRPw5s0bva+7du0aMjIyULVq1SKsLe/NmzcYN24cfHx8IJPJUL9+fWzfvt3g158/fx7h4eFwcnKCo6MjWrdujdOnTwvKnDx5UvBc1L9iY2NFX48QUrYZFDz37duHnJwcc9eFmEB8fDwAYPfu3Th79ix+++03tGnTBl999RXGjBmj8zXp6eno378/Zs6cCUdHx6KsLgAgMjISGzZswJw5c/Drr7+iYcOG6N27N7Zu3Vrga+Pi4tCyZUtkZGRg06ZN2LRpExQKBUJDQ3H27Fmt8gsXLsTZs2cFX3Xq1DH6eqT4SpIDu67y3wkxOWYAjuNY+fLl2eTJk9m///5ryEvKDLlczgAwuVxu6aowxhiLiIhgMpmMZWdnq47l5OSwSpUqMU9PT63yb9++ZZ06dWIDBgxgubm5RVlVxhhjBw8eZADY1q1bBcfDwsKYj4+P4H3oEh4ezsqXL8/S0tJUx1JSUpinpydr1qyZ6lhMTAwDwHbu3GmS65Hi7d5rxmp+w5j/1/z3e68tXSNiapb+7DWo5fntt9+iSpUq+OKLL1C7dm28//77+PHHH5GWlmbOuE6M8Ndff6FWrVqwtrZWHbOyskK5cuUgkQgnV+fm5mLAgAGwtrbGunXrwHFcUVcXe/fuhaOjI3r27Ck4PnjwYDx69Ajnzp3L9/WnT59GSEgI7O3tVcecnJzQsmVLnDlzBo8fPxZVH1Nfj1jG+YdARjb/74xs/mdCTMmg4Dl8+HCcOXMG//77Lz777DPcvXsXQ4cOhbe3N4YNG1Yk3VkXL15Et27d4OPjA3t7ewQGBmLevHlas0Pj4+PRtm1bODo6wtXVFZGRkbhz547Oa65cuRKBgYGQSqWoXLkyoqKikJWVZfb3Yi7JyclISkpC3bp1BcefPn2KhIQErQD18ccf4/Hjx9ixY4dWYNWHMYbs7GyDvgxx5coV1KpVS+v+wcHBqvP5efv2LaRSqdZx5bHLly8Ljo8ePRoSiQTOzs4IDw/Hn3/+WajrkeKpkS9g999/UnYS/mdCTMqY5mpOTg775ZdfWGRkJJNKpczKyooFBgayL774gj158sTEjWPGEhISmEwmY/Xq1WM7duxgx48fZ3PmzGHW1tbsgw8+UJX7999/mZOTE2vRogU7ePAg2717N6tduzbz8fFhz549E1xz/vz5jOM4Nm3aNBYTE8OWLl3KbG1t2fDhw0XVzdJdB+qOHDnCALDFixezrKwslp6ezmJjY1njxo1ZREQES01NVZVNTExkAJhMJmMODg6qr99//z3feyi7Pw35unv3boF1rl69OgsPD9c6/ujRIwaALVy4MN/X169fn9WoUYPl5OSojmVlZbEqVaoIuoPj4+PZ2LFj2d69e9nvv//O1q9fz2rVqsWsra3Z4cOHRV+PFH/3XjO2M4G6bEsrS3/2GhU81b148YJ99dVXrF69eszKyorZ2tqaol4CM2bMYADYrVu3BMdHjBjBALCXL18yxhjr2bMn8/T0FDzMxMREZmNjwyZPniyos0wmYyNGjBBcb8GCBYzjOJaQkGBw3Sz9C1S3aNEinUEsLCyMvX371iT3SElJYXFxcQZ9ZWZmFni96tWrs/bt22sdVwbPRYsW5fv6devWMQBs5MiR7MGDBywpKYkNHTqUWVtbMwBs+/btel/76tUr5ufnx4KDg01yPUJI0bH0Z2+hUwA5OzvDz88P3t7eAGBwd50YNjY2AAAXFxfBcVdXV1hZWcHW1hbZ2dk4cOAAunfvLljzExAQgNatW2Pv3r2qY4cPH4ZCocDgwYMF1xs8eDAYY9i3b5/J30NRiI+Ph7W1Nc6cOYO4uDgcPnwYbdq0wdGjR/H999+b5B6Ojo6oX7++QV+2trYFXs/DwwPJyclax1++fAmAX3uanyFDhmDx4sXYtGkT/Pz84O/vj6tXr+Kzzz4DAPj66u+vc3V1RefOnfHPP/8gIyOj0NcjhJQdRgfPv//+G2PHjoWPjw969eqFCxcu4H//+x8uXbpkwurxBg4cCFdXV4wcORJ37txBamoqDhw4gG+//RajR4+Gg4MDbt++jYyMDNVYmbrg4GDcunULCoUCQN44mubYoLe3Nzw9PQscZyuu4uPjERQUhKZNm+K9995DeHg4du/eDQcHB/zwww8mucepU6dgY2Nj0FdiYmKB16tbty7+/fdfrT+6lGOL6stI9JkyZQpevHiBy5cvIzExEWfOnMGrV6/g4OCABg0a5Pta9l+OEPXJUoW5HiG0RKZsEJUY/tWrV9i6dSvWr1+PS5cugeM4hIaGYujQoejWrZtBLQ1jVKpUCWfPnkVERIRgEf+YMWPw9ddfA4Cq9aKrpeLu7g7GGF69egVvb28kJydDKpXCwcFBZ1ldLSGlzMxMZGZmqn5OSUkx9m2ZlFwux507d7Ra08pJU5s2bcKdO3dQpUqVQt2nQYMGiIuLM6isj49PgWUiIiLw/fffY/fu3ejVq5fq+IYNG+Dj44PGjRsbdC+pVKoKtElJSdixYweGDx8OOzs7va959eoVDhw4gPr160MmkxX6eoQkyYF2m/kZvnYS4Eg/wN+l4NeRkseg4HnkyBGsX78e+/fvR2ZmJgICAjB79mwMHjwY/v7+5q4jEhMT0aVLF5QvXx67du1CuXLlcO7cOcyfPx9v3rzBunXrVGXzW26hfs7QcpoWLVqEqKgoke/A/OLj48EYQ6NGjbTO9ejRA5s2bcLevXsxceLEQt3HyckJ7733XqGuoa5Dhw4ICwvDyJEjkZKSgmrVqmHbtm04fPgwNm/erFpyc+rUKYSGhmL27NmYPXu26vVXrlzB7t278d5770EqleLvv//G4sWLUb16dXz++eeqcn369IG/vz/ee+89eHp6qjY6ePr0KaKjo0VfjxBddC2RoeBZShkyMMpxHJPJZKxXr17syJEjRb6YvlevXszLy4u9efNGcHz9+vUMADt58iS7du0aA8BWrVql9frPPvuMcRzHMjIyGGOMTZ06lQEQLIRX8vT0ZL1799ZbF4VCweRyuerr/v37xWLC0LJlyxgA9tdff2mdUygUzMnJib3//vsWqFnBUlNT2ZgxY1iFChWYra0tCw4OZtu2bROUUc7ynTNnjuD49evXWcuWLZm7uzuztbVl1apVYzNnztT6b2XRokWsfv36zMXFhVlbW7Ny5cqxiIgIdv78eaOuR4gulJyh6Fh6wpBBwXPFihUsOTnZ3HXRq2bNmiwkJETr+OXLlxkA9s0337CsrCxmZ2fHPvnkE61y4eHhrHr16qqft2zZwgCw2NhYQbnHjx8zAGzBggUG183Sv0BCSPFCS2SKhqU/ew2aMDRmzBitsUS5XI7ffvsNW7ZswatXr0zWEtbFx8cHCQkJWonNlckZ/Pz8IJFI0KVLF+zZswepqamqMklJSYiJiUFkZKTqWPv27SGTyQTddQAQHR0NjuPQrVs3s70XQkjp5u8C9Aii7tpSz5iIO2/ePGZvb884jmNWVlaqrsI2bdoUuC7PGPv372ccx7EmTZqokiQsWLCAOTo6sqCgINV6wn///Zc5Ojqyli1bskOHDrE9e/awOnXq5JskYfr06ezkyZPsiy++YFKptEQnSSCEkLLC0p+9ooPnqlWrmLW1Nfvf//7HDh06xDiOUwXPL7/80mzjaidOnGDt2rVjFSpUYHZ2dqxGjRps4sSJ7MWLF4JyFy5cYKGhocze3p45Ozuzbt26aSVXUFqxYgWrUaMGs7W1Zf7+/mzOnDmikwlY+hdICCFlkaU/ew3aDFtdUFAQOnfujKVLlyInJwc2Nja4cOEC3n33XRw4cADDhg3DkydPzNFILpYsvSErIYSURZb+7BWdJOHOnTsIDw/Xec7JyQmvX78ubJ0IIYSQYk108HRxccHTp091nktMTISXl1ehK0UIIYQUZ6KDZ2hoKJYuXSrYy5PjOGRnZ2PNmjV6W6WEEEJIaSEqPR8AzJs3Dw0bNkRQUBAiIiLAcRy++eYbXLx4EUlJSfjpp5/MUU9CCCGk2BDd8qxWrRpOnz6NWrVqYfXq1WCMYePGjfD09MQff/xRJOn6CCGEEEsS3fIE+Bm3hw8fRmZmJpKTk+Hm5kYJswkhxEKS5Hwe3Ua+lJyhqBgVPJWkUqlBO2cQQggxD9rJxTIMDp4xMTHw9PRU7YHJGMPYsWMFZVxcXGjnCUIIKUK0k4tlGBQ8T548ibZt2yImJkZ1LDc3F998842gHMdxaNSoEbp06WLaWhJCCNGpkS/f4lS2PBv5WrpGZYNBE4bWrVuH9u3bo2XLllrnLly4gNzcXOTm5mLYsGHYsmWLyStJCCFEN38Xvqt2eRh12RYlg4Lnn3/+iR49ehRYrl27djh//nyhK0UIIcRwtJNL0TMoeD558gTVqlUTvtDKCl27doWbm5vqmJeXV5nKa0sIIaRsMmjM09raGpmZmYJjHMdh7969gmOZmZmwshK9dJQQUsLRUglS1hgU6fz9/REfH19gufj4eEqSQEgZo1wqMfEo/z1JbukalQxJcmDXVXpeJZVBwbNt27ZYu3atIJ+tpjdv3mDNmjVo166dySpHCCn+dC2VIPmjPzhKPoOC54QJE/D8+XO0bdsWFy9e1DofHx+PsLAwvHjxAuPGjTN1HQkhxZhyqQRASyUMRX9wlHwGjXlWqlQJW7ZsQZ8+ffDee+/B29sbAQEBAIB79+7h8ePHsLOzw7Zt21CpUiVz1pcQUswol0qUljHPwo7fGvJ6WptZ8nGMMWZo4Tt37mD58uU4duwY7t+/D4AfDw0LC8P48eNRpUoVs1W0uLL0buaEENMpbKo7Ma+nSVaFY+nPXlG5batUqYJVq1aZqy6EEGJRhU11J+b1/i4UNEsyWldCCCH/Kez4LY3/lh0GBc9Ro0aJTn6wZ88eStVHCClRCpvqjlLllR0GBc/r16+jSpUq6NevH3777Tekp6frLHfr1i0sX74cderUwdChQwXZhwghpCQobKo7SpVXNhg05nn8+HHs378fixYtQocOHSCRSFC9enV4eXlBJpPh5cuXuHPnDl6+fAkHBwcMGjQIM2fOhJeXl7nrTwghhBQ5UbNtAeDixYs4cOAAYmNj8ejRI2RkZMDT0xOBgYEICQlB165d4eTkZK76FjuWnvFFCCFlkaU/e0XNtgWAd955B++884456kIIIYSUCDTblhBCCBGJgichhBAiEgVPQgghRCQKnoQQQohIFDwJIYQQkSh4EkIIISJR8CSEEEJEouBJCCGEiETBkxBCCBGJgichhBAiEgVPQgghRCQKnoQQQohIFDwJIcRMkuTArqv8d1K6iN5VhRBCSMGS5EC7zUBGNmAnAY70ow2ySxNqeRJCiq2S3HI7/5APnAD//fxDy9aHmBa1PAkhxVJJb7k18uXrrax/I19L14iYkuiW5+7du5Gbm2uOuhBCilBRtuqMuVdJb7n5u/ABf3lYyQv8pGCiW549e/aEr68vPvnkEwwfPhxeXl7mqBchxIyKslVn7L1KQ8vN34WCZmkluuV58uRJNG3aFFFRUfD390f//v0RGxtrjroRQsykKFt1xt7L3C23kjyeSixPdPBs2bIlfvrpJ9y7dw+TJ0/G8ePH0bx5czRo0ADR0dHIzMw0Rz0JISakbNUB5m/VFeZe/i5AjyDzBM52m4GJR/nvFECJWBxjjBXmAtnZ2di5cye+/PJLxMfHw93dHcOGDcOYMWPg7e1tqnoWWykpKXBxcYFcLoezs7Olq0OIwZLkfCuwka95u2zPPwT8nIEHKea9lxi7rvKBU2l5GB+kSclh6c/eQi9VuXv3Ls6dO4ebN2/C2toadevWxYoVK1CjRg388ssvpqgjIcQMzNWqU1Jv3Q3aX3wCJ1C0LW9SOhkVPBlj+PnnnxEeHo5atWph69at+PTTT5GYmIgTJ04gMTERISEhGD9+vKnrSwgpIYrzbFmaCUsKS/Rs2yVLlmDt2rW4d+8e6tWrh++//x59+vSBVCpVlfHy8sKkSZPQunVrk1aWEFJyFPfZsjQTlhSG6OA5c+ZMfPDBB4iOjkarVq30lqtatSpmz55dqMoRQkouZevO3OOqhFiC6AlD9+7dQ0BAgLnqU+JYetCaEHMpiglFhBjL0p+9olueFDgJKf2KU2o8CuKkOBIdPIcMGaL3nJWVFVxdXdGwYUNERETA1ta2UJUjhFiGrsk+5ghcBQXG4hTECVEnOnjGxMRALpfj9evXkEgk8PDwQHJyMrKzs+Hq6grGGL788kvUrFkTJ0+eRPny5c1Rb0KIGRXFZB9DAmNRBXFCxDIqMbyTkxO2bduGjIwMPH78GBkZGdi6dSucnJzw22+/4c8//8SrV68wffp0c9SZEGJmxi7lEJPyzpClLLQekxRXolueEyZMwGeffYZevXqpjllbW+Ojjz7C06dPMWHCBPz555+YMmUKli1bZtLKEkKKjtilHGK7WA1p3eY3Y5fGQokliW55xsXFIShIdx6rOnXq4OLFiwCA+vXr48WLF4WrHSGkxBCbFMHQ1q2uTEiUm5ZYmujg6ezsjJiYGJ3nTpw4oZoynJGRAScnp8LVTsOff/6Jjh07ws3NDXZ2dqhevTo+//xzQZn4+Hi0bdsWjo6OcHV1RWRkJO7cuaPzeitXrkRgYCCkUikqV66MqKgoZGVlmbTOhJQVxnSxGpsisDhnLyJlg+hu2z59+mDJkiVgjKFnz54oX748nj59ih07dmD58uUYO3YsAOCvv/5CrVq1TFbRrVu3on///vjwww+xceNGODo64vbt23j06JGqzLVr1xASEoL69evjp59+gkKhwOzZs9GiRQtcunQJ5cqVU5VdsGABZs2ahalTp6Jdu3aIi4vDzJkz8fDhQ3z33XcmqzchZYWYpAiF7XIt7tmLSBnARMrMzGS9e/dmHMcxKysr1RfHcaxPnz7s7du3jDHGjh49ys6dOyf28jo9ePCAOTg4sJEjR+ZbrmfPnszT05PJ5XLVscTERGZjY8MmT56sOvbixQsmk8nYiBEjBK9fsGAB4ziOJSQkGFw3uVzOAAjuSQjR795rxmp+w5j/1/z3e6+Nv87OBONfT0o2S3/2iu62tbW1xdatW5GQkIBVq1Zh3rx5WLVqFa5cuYItW7bAxsYGANC2bVs0atTIJAH+hx9+QFpaGqZMmaK3THZ2Ng4cOIDu3bsLsk0EBASgdevW2Lt3r+rY4cOHoVAoMHjwYME1Bg8eDMYY9u3bZ5J6E0K0marL1dy7whCSH1HdthkZGahWrRrWrl2LLl26mLRbNj+///473N3dce3aNXTt2hVXrlyBu7s7IiMjsXTpUjg7O+P27dvIyMhAcHCw1uuDg4Nx9OhRKBQKyGQyXLlyBQBQt25dQTlvb294enqqzhNCTI+6XElpIKrlaWdnh4yMDDg4OJirPjo9fPgQ6enp6NmzJ3r16oVjx45h0qRJ2LhxIzp27AjGGJKTkwEA7u7uWq93d3cHYwyvXr0CACQnJ0Mqlep8H+7u7qpr6ZKZmYmUlBTBFyHEcLQdGCkNRE8YCg0NxbFjx9CmTRtz1Een3NxcKBQKzJkzB1OnTgUAhISEwNbWFuPGjcPx48dhb28PAOA4Tu911M8ZWk7TokWLEBUVJfYtEELU5LeGlNZvkpJA9Jjn9OnTsX37dsybNw9XrlxBcnIyXr58KfgyNQ8PDwBAeHi44HiHDh0A8MtTlGV0tRpfvnwJjuPg6uqqup5CoUB6errOsrpar0rTpk2DXC5Xfd2/f9+o90QI0UbrN0lJITp4NmjQAImJiZg7dy7q1asHLy8vlCtXTvBlarrGMQGA/bebmpWVFapWrQo7OztcvnxZq9zly5dRrVo1yGQyAHljnZplnzx5ghcvXqBOnTp66yKVSuHs7Cz4IoSYBq3fJCWF6G7b2bNn59utaQ7du3fHd999h19//RXvvPOO6vihQ4cAAE2aNIFEIkGXLl2wZ88eLF26VJWgISkpCTExMRg/frzqde3bt4dMJkN0dDQaN26sOh4dHQ2O49CtW7eieWOEEAGaTERKDIsskDFCly5dmFQqZZ9//jk7evQoW7RoEZPJZKxz586qMv/++y9zdHRkLVu2ZIcOHWJ79uxhderUYT4+PuzZs2eC682fP59xHMemT5/OTp48yb744gsmlUrZ8OHDRdXL0muNCCltaP0mMYSlP3sLFTzT09PZgwcPWFZWlqnqk++9pkyZwipWrMgkEgnz9/dn06ZNYwqFQlDuwoULLDQ0lNnb2zNnZ2fWrVs3duvWLZ3XXLFiBatRowaztbVl/v7+bM6cOaokD4ay9C+QkJLi7H3GJh/lvxNSWJb+7OUY+2/gUISYmBhMnz4dcXFxAIDz58/j3XffxejRoxEaGorIyEiTto6Ls5SUFLi4uEAul9P4JylTDNnIWnn+USrQa3feuR3dgSZ+RVdXUvpY+rNX9IShEydOoF27dlAoFPjss8+Qm5urOufp6Yno6GhT1o8QYiZi9t7U9VrlrNg2G4HYB/rPt9sMbPpHeH7vNePrTUhxIDp4zp49Gx07dsTFixcxf/58wbl69erh0qVLpqobIcRMCrskRH1WbFYu0G+v8Bqas2YrajQMIgKNr7sYhfkDgZD8iJ5te/HiRezcuROAdjKBcuXK4dmzZ6apGSHEbHQtCRGTkKCRL2BjxQdOgP+ufg3NWbN96gIhlfgWZ0Rg0XTZit2cmxAxRAdPiUSid8/LZ8+emXwPT0KI6RV2SYi/C7A5gm9xZuVqX0PX9mT+LtpB05zZhAr7BwIh+REdPBs2bIhNmzaha9euWud27dqFpk2bmqRihBDzEbP3pj5N/IATA/RfI78UfID5W4a0ZpSYk+jgOXXqVISHhyMiIgIDBgwAx3E4d+4c1q9fj127diEmJsYc9SSEmFhBwc3c19CXTchULVFT/IFAiD6ig2fbtm2xYcMGjBs3Dvv37wcAjB49Gq6uroiOjsb7779v8koSQkofzZahn7PpW6LK1ysDMwVQYiqigycA9OvXD927d8eZM2fw9OlTeHp6onnz5kW+VRkhpOTSbBkaM0ZpyFpTmjREzMGo4Anwe3uGhoaasi6EkGJMPVABpukO1ez2NWSMUlkPP2dg0P78AyNNGiLmYlTwZIwhLi4O9+7dQ0ZGhtb5AQMGFLpihBDzEjPTVb0FJ5MAjAGZOaZtzRkyRqleD/WlMvoCI00aIuYiOnjeuHEDH3zwAW7evAldmf04jqPgSUgxZ0h3pnpwVW/BKbLzymRkA7/eBDzsCw7ChgTrgiYgaSZnUAZQfYGRJg0RcxEdPEePHg2FQoEdO3YgODgYUqnUHPUihJhRQd2ZmsE1umteC0695Sm1Br48xwfU/Fqhphp71GxJRncFHqQULiATYgzRwfP8+fP4/vvv0aNHD3PUhxBSBHR1Z+praWZk8wFKvQUH8P9OzgAW/plXTt+YoqnGHqklSYoL0cHT0dGRdg8hpITTDEIAELaZb0HKJMCGrtrBVbMF5+/CB9yvYgseUzTl2CO1JElxIDp4Dh48GFu3bkX79u3NUR9CSBFRD0Lf/pU3lqnIBv5+algLz9CWILUYSWkjOnjWqVMH27ZtwwcffIAuXbrAw8NDq0xZ2s+TkFJBc+4fM7yFZ+pyhJQEojfDtrLKfxczjuOQk5NTqEqVJJbekJUQU0iSA2035U0COtafAh0p3iz92Su65Um5awkpPky1K4m/Cx8wqVuVEMOIDp6tWrUyRz0IISKZOvUcdasSYrj8+2DzIZfL8dtvv2HLli149eqVKetECDGAvl1JkuTArqv8d0KIeRiVnu/zzz/H4sWLkZGRAY7jEBcXBzc3N4SGhiIsLAxTp041dT0JIRr0rdWkROiEmJ/olufq1asRFRWFoUOH4uDBg4IUfZ07d8bBgwdNWkFCiG7K5R/Lw/KCpL7WqBgFtVxN1bKlFjIpyUS3PL/55htMmDABS5cu1ZpVW716ddy8edNklSOE5E9znLKwyQgKarkact6QSUfUQiYlnejgeefOHYSHh+s85+TkhNevXxe2ToQQI4lJRqAr0BWURk/z/K83gY/fy7ueoQGRtgojJZ3o4Oni4oKnT5/qPJeYmAgvL69CV4oQYjxDZs3qC3QFtVwb+fLrQDP/63T68hzQobr+LmNlCj/NIE1bhZGSTnTwDA0NxdKlS9G1a1fIZDIAfGKE7OxsrFmzRm+rlBBSfOgLdAW1XP1dgIlN85LBK9ReK2YCE6XrIyWd6OA5b948NGzYEEFBQYiIiADHcfjmm29w8eJFJCUl4aeffjJHPQkhJpRfy6+glmuHarqTwesKiLuu6u+epXWlpCQTnZ4PAK5evYoJEybgxIkTyM7OhrW1NVq3bo0VK1agVq1a5qhnsWXpFFGEGEJX16nymJ9zwXtiGnI9feVoYhAxB0t/9hoVPJUyMzORnJwMNzc32NnZmbJeJYalf4GEFCS/AFYUwc1UKQQJUWfpz16jMwwBgFQqhY+PT5kNnKTsKI5rEg2tU35rP02xLrQg/i5AjyAKnKR0MSrDECFlSXHsetRXJ7EzW2nWKyHGoeBJSAGK45pEfS3G/Ga2/noLWvt20qxXQoxDwZOQAhSH1plmi1JXnQoK8soZsl+dE7aeadYrIeJR8CSkAOZqnRV2xqquOikDqkwCJKfzr80vgQEhxDgUPAkxgKlbZ6ZIZaeeUk9ZR2X37LKzwMLTfAago/2KR+uZkNLEoOD5+++/i7poy5YtjaoMIWWFmJagvsAX+wDotxfIyuVT5k1syicwAIC3/6XPU2TzwfTjBjS2SYgpGRQ8Q0JCwHEcAIAxpvq3Ppq7rRBS1hTUJSumJairizZJnhc4AT7X7MI/+XHNQfU0LmD0Sm7ToHWepDQyKHjGxMSo/p2amopPP/0UNWvWRJ8+fVChQgU8efIEW7ZswfXr17Fq1SqzVZaQksCQLlmx46ia3cbnH+YFTnUZ2YCbXV7ydqk1n7jdUsttiuMyH0JMwaDg2apVK9W/R40ahZYtW2Ljxo2CMgMHDkT//v3xyy+/oHPnzqatJSEliKFdsoUZR1VvuUo4wNqKD5Z2Er7rtkM1w3PMFqQwLUeaqERKK9EThnbu3InNmzfrPNe3b1/0798f3377baErRkhJVZjJOYYGKs2WK6D9OlNskl3YliNNVCKllejgmZ6ejmfPnuk89/TpU6Snpxe6UoSUZMYubREbqDRbroYGWz9n4Qzd/BjScswv4FMSBlJaiQ6eLVq0wIwZM/DOO++gTp06quOXL1/GzJkz0aJFC5NWkJCSyJguWbFdnGK6U9V3UBm03/AAXVDL0dDxXQqapLQRHTxXrFiBli1bon79+qhdu7ZqwlBCQgI8PDywYsUKc9STkFIvv0ClGSjFtFLVy1pzQM5/s28NCdAFtRxpTJOUVaKDZ82aNXH58mV8+eWXOHnyJO7cuQMPDw9MmjQJ48aNQ/ny5c1RT0JKLfXAqCtQ6QqUYoKWelll4AQAGyu+JVqQ/FqONKZJyiqjMgx5eXlh8eLFpq4LIWWOrsDYI0hYRlegFBO0GvnygVJ9aYs1x/88aH/hlo/QmCYpq4xOzyeXyxEbG4sXL16gY8eOcHNzM2W9CCkT9LUg1VujugKlmKDl7wJsjshLqqAeSE3R1UpjmqQsMip4fv7551i8eDEyMjLAcRzi4uLg5uaG0NBQhIWFYerUqaauJyGlkq7AqKs1qitQiglaTfyAEwN0TxqirlZCxLMS+4LVq1cjKioKQ4cOxcGDB8FY3iBK586dcfDgQZNWkJDSTNmCXB6W132qrzXaI6jwLcQeQXwg1bwnIUQc0S3Pb775BhMmTMDSpUu1cthWr14dN2/eNFnlCCkLNFuQRZVkgYImIcYTHTzv3LmD8PBwneecnJzw+vXrwtaJkDLNVEkWorsCD1IM2y9UcxkMTQAiJH+ig6eLiwuePn2q81xiYiK8vLwKXSlCyjpTJFlQThDKbx2oroArJokCIWWV6DHP0NBQLF26FGlpaapjHMchOzsba9as0dsqJYTkL0nOJ3BPkhv3emV3L6B7Rq0umgF37zXt8VZCiDbRLc958+ahYcOGCAoKQkREBDiOwzfffIOLFy8iKSkJP/30kznqSUiplSQHfr0JfHmO37xaX4uvoO5Uzfy1hsyo1RxfjQgE9l+nmbiEFIRj6tNlDXT16lVMmDABJ06cQHZ2NqytrdG6dWusWLECtWrVMkc9i62UlBS4uLhALpfD2dmAdC2kxDDH2J8yUILjtw0D8rpN1S0PEyZLMGZ3E0PrT2OepCSy9GevUes8g4KCcPjwYWRmZiI5ORlubm6ws7Mzdd0IsRhzbOKcJAfabuL33QSAL2OBCY21A6euFp8xOWQNHTfVtTsLBU1C8id6zHPIkCG4e/cuAEAqlcLHx0cVOO/du4chQ4aYtoaEWICuYGWKa2aqre5SZAPg8sYppdbA9Pd1B2r18UzqTiXE8kR321pZWSE2NhaNGjXSOvfXX3+hUaNGWus/SzNLdx0Q8yio5WlI16au7lD1lqdMAhztx//bmO5VQsoyS3/2im555ufly5eQSqWmvKReP/zwAziOg6Ojo9a5+Ph4tG3bFo6OjnB1dUVkZCTu3Lmj8zorV65EYGAgpFIpKleujKioKGRlZZm7+qSY05X5R0kZWCce5b/rmh2rq4y/C3CsPzC9Od/CPNovr4u0oOxBxgTOws7eJYToZ9CY5++//46TJ0+qfv7hhx9w+PBhQZmMjAzs378fQUEaW0KYwcOHD/HZZ5/Bx8cHcrnwk+HatWsICQlB/fr18dNPP0GhUGD27Nlo0aIFLl26hHLlyqnKLliwALNmzcLUqVPRrl07xMXFYebMmXj48CG+++47s78PUrzpG/szZPxRXxl/F+Dj9/TfU1eQNHaykKnHbAkheQwKnjExMYiKigLAr+n84YcfdJYLCAjAqlWrTFc7PT755BO0bNkS7u7u2LVrl+Dc7NmzIZVKceDAAVVTvkGDBqhevTqWLVuGJUuWAACSk5Mxf/58DB8+HAsXLgQAhISEICsrCzNnzsS4ceOK5A8BUrzpCmaGpM8zJsWevoBnzGQh2qSaEPMyqNt28uTJeP78OZ49ewbGGA4fPoznz58LvlJSUnD37l20bt3arBXevHkzTp06hdWrV2udy87OxoEDB9C9e3dBH3hAQABat26NvXv3qo4dPnwYCoUCgwcPFlxj8ODBYIxh3759ZnsPpOgUputSX/dsfl26SoaU0aQZ8H69xf+7kS8/Pgrw35WBOPYBMOUY/12TrglG1I1LiOkY1PK0s7NTzai9e/cufHx8YGNjY9aK6fLs2TOMGzcOixcvhp+fn9b527dvIyMjA8HBwVrngoODcfToUSgUCshkMly5cgUAULduXUE5b29veHp6qs6TkquwXZe6gpmHXV4r1FRLRZQa+QK21sDb/yYULT/LrwV9lApk/3dMOb1v/zVgzG/8v7cnADu687ulqN9bPT8uQN24hJiS6AlDmZmZOHPmjM5zp06dMuuuKqNGjULNmjUxcuRIneeTk5MBAO7u7lrn3N3dwRjDq1evVGWlUikcHBx0llVeS1NmZiZSUlIEX6R4KuxyE/XWm0zCB7OJR4E2G3W39kxBfe57Zg6fUKHfXiCbCY+NPyJ83d5r2i1L9YlI5lh6Q0hZJjpJwoQJE1CjRg20atVK69wvv/yCGzdu4OeffzZJ5dTt3r0bv/zyCy5evAiO4/Itm9959XOGllO3aNEi1fgvKd4Ks7UXIGy9JacDC0/zx7Nygb57gJiBpkmcoEynt/daXj5apVeZwmM2VvyxHI0FZs388m9ZFvZZEEKERLc84+Li0LJlS53nWrVqhbi4uEJXStObN28wevRo/O9//4OPjw9ev36N169f4+3btwCA169fIy0tDR4eHgCgs9X48uVLcBwHV1dXAICHhwcUCgXS09N1ltXVegWAadOmQS6Xq77u379vondJTM2YcUdd1+gRBHSoDkjU/p7KZv+l2dPDkPFF9THVXrv57ldd1JO9Lw8DfryUd84KwP+F8wE2v5alKZ4FISSP6JanXC7XubYS4MdGld2ipvTixQs8ffoUy5cvx/Lly7XOu7m5oWvXrti1axfs7Oxw+fJlrTKXL19GtWrVIJPJAOSNdV6+fBmNGzdWlXvy5AlevHiBOnXq6KyLVCotsrWspPBMlWrO3wUY0QBYfUHtoJ6OC11jrYD2rF31rlR93GTCscvzD//LTPSfqe8DXQP5exbUsqS0e4SYjujg6evri/Pnz6Nt27Za586fPw9vb2+TVExdhQoVEBMTo3V88eLFOHXqFH799Vd4enpCIpGgS5cu2LNnD5YuXQonJycAQFJSEmJiYjB+/HjVa9u3bw+ZTIbo6GhB8IyOjgbHcejWrZvJ3wcp2XrXAdZf4oOXTMJP5lFfygL818WboT3R6KtY7S5V9a5UXaTW/D00g556kFQmlzd2A21CiHFEB89u3bph8eLFaNq0qWBZysmTJ7FkyRIMHTrUpBUEAJlMhpCQEK3j0dHRsLa2FpyLiopCw4YN0blzZ0ydOlWVJMHT0xMTJ05UlXN3d8fMmTMxa9YsuLu7q5IkzJ07F8OGDaM1nkSLvwufFUjXDFaZhJ/sk5nDBz2ZJG97MTD9CROUAS/xFbDygvB+Q+rndb8qg2F+QZJaloQUHdG5beVyOZo3b45///0XNWrUgJ+fHx48eIAbN24gKCgIp0+fLrI8g4MGDcKuXbvw5s0bwfG//voLU6ZMwdmzZyGRSNCmTRssW7YMVatW1brG//3f/2HVqlVITExEhQoVMHjwYMyYMcPgpTiWzq9IjGOKbbh2XeXHK3WZ3hzwsDd8mYjmtSQcILHOf39PQsoyS3/2GrWf55s3b/DVV1+pkiWUK1cOHTp0wLhx4/SOh5ZWlv4FEvE0xySjuwo3jjYkUMU+ADb9Axy5DbzNFbY8C0okD+Sfgs/GChj+rnB8VXN/T0LKOkt/9hq1n6ejoyNmzZqFWbNmmbo+hJid5prHvdcKTmWnHvwepfKzY5VaVwIiagJP0oBXGYCbjq1tlV2q+hI36Epq8OMlWlpCSHFlVPAkpCTTXPMYEQjsv647UCXJ+SUpX57L60INrSy8Xkwi/6Xuq1jdLdj8cs5qjlnSBCBCii/RwbNNmzb5nuc4DsePHze6QoSYm65JN7oCleb+mwAf8Coa0EOkrwUrJlmBrglAxuwjSggxPdHBMzc3Vyv7zosXL3D9+nV4eXmhRo0aJqscIaamHljUxxB1BapfbwoDJ8AHvD51gZBKwjFPTfmttTS2RWlIrl7aioyQoiE6eKrv66nuxo0b6Nq1K+bMmVPYOhFiFqIDi0YShHrl+U2slYHWxwk4orbHujWAr/7L9pNfYDR2Scmvt4zfR5QQYlqi0/PpU6NGDUyaNAmTJ0821SUJKTT1NHlik6PXKy/8+e+n/KxcZcq98w/zdkABgBwAT97kJWM3pSQ5n5heSaanZatrKzJCiOmZdMJQpUqVaCsvUmzoWpIiJjn6Ax0b5qi35vyc+b8+1XttX2ea8h3kOf9Q2IU8oXH++4jSmCch5mXS4Ll79274+PiY8pKEGE2zpfkgRVxgaeTLZwvSHPf0c+YD86D9wsAJAK5S0yRg0FUXQVq+6vrLUqYhQsxPdPAcMmSI1rHMzEz8888/uHr1KpYuXWqSihFSWLpmtooNLNk52scepPBduJo5aW2tgHoVCp+AQRdqURJSvIgOnidOnNCabSuTyVCpUiVMmzYNffr0MVnlCCkMsQFHs4W49TI/jqnpVBJw8Ib2cSsr4O8n4hMwiHk/FDQJKR5EB8/ExEQzVIMQ89AMOPq6UHWNj/5wUfc1f76u+7giGwBneAIGQkjJJSp4ZmRkoG3btoiKitK5JRkhxVl+S1U0x0c3/cMvOTGExArIzs3bpqxDNd0JGPyctXdJIYSUTKKCp3KjaYmEsvqRkie/NZDq46MyCZ/8QMmaAz6qA4AB2xOAHLWtFGw4QDmKkZMDLDkN9A/WTsAAUPICQkoT0es8mzZtivPnz5ujLoSYVX5rIJUtxOVh/DIQ9axBOQzYchnYcoX/t40V0LcuMLIBEF4tb61nFgMO3OSTxsc+EN5b7BpTQkjxJroJuXz5cnTt2hUVKlRAZGRkmduCjJRcBU0gUv786y1+P81sPZv1ZeUCAc7AV+e0Z9wq7b0GNPHL+1lMTltCSPEnej9PJycnvH37FtnZ/KeGvb29YPYtx3GQy+WmrWUxZuk95UjhKScR+TnnLSuxsdI/5mknAcY3ARb+mXesdSXhzio7uguDp/p9aKkJIYVn6c9e0S3P7t27ay1VIaSk0tyEWhkwNQNnY1+gb528vLUAv+2YsiU5L4Tf53PvNX6GrWbgBGipCSGlieiWJxGy9F8/xHhJcmBVHD8JSBfN7EKaE32oJUmI5Vj6s1f0hKF58+bh0aNHOs89fvwY8+bNK3SlCDEXZaL42Ad8i1M9cNpo/N8wsQnwUe28nzUn+vi7mCcJPCGk+BMdPKOiovDgwQOd5x49eoSoqKhCV4oQU1DfUUX5c7vNwMSjQL+9wsk+H9UGNkcIZ+N2qA6Mbki7lBBCtIke88yvl/fNmzewsbEpVIUIMQXNsczlYfzem8qAqTmmqRyn1DUbV/MYddcSQgwKnv/88w8uXbqk+vnQoUO4du2aoExGRga2bNmCqlWrmrSChBhDfV1lVi4w5jfhec2lKA9S9AdF9Yk+sQ/4VmtWLiU7IKQsMyh47t27V9Udy3Gc3nFNOzs7/Pjjj6arHSFGauSb/3KT4Q2A6Et5s2X9nIUt1c0RupeaKAMnwJf99SbgYU+tUELKGoNm2z5+/BiPHj0CYwyNGjXCjz/+iDp16gjKSKVSVK1aFXZ2dmarbHFk6RlfpZUpukbVW4ma/i8ceJIGgPFjm+cf8mOhSjZWwIkBwnvvuiosI+EAiTWfEJ5aoYQULUt/9hrU8vT29oa3tzcAICYmBu+++y6cnJzMWjFSduWXwN3Q1ysD74kB/L9vvwJWX8grM/EoH1RlEn4PzuQMYVduVq729mHqWYJsrIDh7+Zds7DbjRFCShbRE4ZatWqldezcuXO4ePEiWrVqhVq1apmkYqTsyi+Be0F0Bd4eQXwr9Pt4Piiqd+cqsoG+e/igaWsN2LC88UzNmbWa6f0A4MdLlHKPkLJIdPAcNmwYsrOzER0dDQDYvn07+vbtC8YYbG1tERMTg6ZNm5q6nqQMKUweWH0J2Afs44OihAM+rM0neldStjbf5gDTm+c/hqmZJUjMZtuEkNJD9DrPmJgYtGnTRvXzggULEB4ejkuXLqFZs2ZYuHChSStIyh71HU7Edtnq2jnl15t5mYKyGZD6ls8eBPCtTeW/7f7rwhVbV0qUQEjZI7rl+eTJEwQEBADgkyIkJCRg9erVCA4OxtixY/HJJ5+YvJKk7FFv4emaPKSezP3vp1BN/FEG3l9vAlCmYNZIxfzzdX7CkHqeWs3E8KaeAERrQwkpXUQHTxsbGygUCgDA6dOnIZPJ0KRJEwCAm5sbXr9+bdIKkrJNcwwzuivw9xPgy3P8eKW65bHAsf78v5XbhX1xhm/BWnPCTazPPACWtM372d+Fn01r7FirmPdAs3IJKflEd9sGBgZi06ZNSE1Nxbp169C8eXNVVqEHDx6gXLlyJq8kKbs0xzD77QUWntYOnADfNXv+oXaChIlHgRnvC8tGBGq/Pr/Nsk35HmgjbEJKPtEtz4kTJ+Kjjz7Ctm3bAAD79u1TnTt+/DiCg4NNVjlSNiXJ+Q2pX2UA4PjlJIoC9thU8vtvuZd6SzMrlx/rHPkecO0FEOgJ+OhYaVXQZtnGoo2wCSl9RAfPnj17omLFijhz5gwaNmyIFi1aqM75+fmhe/fuJq0gKVuS5EDYZmHL0tYKmP4+UK983pik1BqY2BR4pQDWqK3fPJkIRP8t7KKVSYDlZ/MmDcUk8ktMJjTmJwg9SMkLlubYc9NcQZkQYjm0n2chWTrLRWmjmcVHafr7wMcNtCfeJMmBtpvyAqNm6/Sj2kAVV76rNz80FklIyWLpz17RLU91z58/R0ZGhtZxf3//wlyWlGGNfPO6aQWY7sB5/iEw9J28TD/qSRBsrPixTR8nfjKR+sbWmihDECFEDNHBMzU1FePHj8e2bdtUs2415eTk8ylFSD78XYANXYFv/+K7YHPBd9HWq6A963bAPj4g2loJx0WnNQcWneYD6KD9fFnuv+UqttbAsHeAdReFwZTGIgkhYogOnuPGjcPWrVsxdOhQBAcHQyqVmqNepIxKkgvHNcOqAP2D+XFJ9Rmrm/7JC35vc4G+QcBPV/mAqQycyrJ7r+W1ZN/mAFXd+CUtyrWd6mOehBBiCNHB8+DBg1i8eDHGjh1rjvqQMk59WUdmDnDgJnD0DrCxm3DGqmage/QmL2Bqzsht5gfsv657tquPk/bWY4QQUhDRwVOhUKBu3brmqAshgmUdSpk5fEszumteKxEA1l/iW5QyCb/8JCZR9zWzcrUTumsmLQB0j6dSi5QQoovo4NmxY0f88ccfgvy2hJiKclnHktN8q1NJ2QKd2DSv3FGNgLj+v3FMzRm3fs7CJSiamYR+vQV8FSscTzVXmj5CSOkgOnjOnDkTPXr0gJOTE7p06QIPDw+tMu7u7iapHCl7lC2+dlWA324Lg2BmDrDwTz7QKVuLSv4ueeOYyenCpSkPUoRlNZMWgAmD6d5r5knTRwgpPUSv87Sy4jP6cRynt0xZmm1r6bVGJZmupSfK7lQlKwASK35SkLrp7wtbi+qtQ0NyyarfG9CeyUstT0KKN0t/9opuec6ePTvfwEmIIXQFOPXJQkq5AHJy+dR6P17ixzjtJHzqPvXW4ao4YHTDvO7ZgjL6FLQvJ2UEIoTkhzIMFZKl//opqTQzCU1vDlRwBMYfEabWU1oexgcy5fKSgfu1EynYWPHllFuNUdAjpPSy9GdvoTIMKRQKvHr1Cm5ubpDJZKaqEyllYh/w44gRgXnLQhr58us4lWs1l50Vds1KAChjo9RamHt211Xdu6pk5QJjfuP/Td2thBBzEr0lGQCcOXMGLVq0gJOTE/z8/ODk5IRWrVrh7Nmzpq4fKeFiHwC9dgPbE/jvsQ/44/4ueTNnAe0xzdrl+a7a6c35iUDqQVB96zCZBJDoGEWgrb8IIeYkuuUZGxuLNm3awNXVFSNGjICPjw8ePnyIPXv2oE2bNjh58iQaN25sjrqSEmjvNe2fla3PDtXyJv1YcUCuWnft30+BG8m6W4/+Lvyknr3X+AQI/74AfrgonJlL6fYIIeYkesyzXbt2eP36NWJiYuDg4KA6npaWhtatW8PNzQ2//fabyStaXFm63704S5IDW68Itwzb0V2Y0WddPDDvD/3XUG5FtukfoKIz0Oe//Byas3JtrfnlLe2q0JgnIWWBpT97jWp5rl+/XhA4AcDBwQGTJk3C0KFDTVY5UnLpWnai2b2aJAcW/Jn/db44DWSp/Xm37iLwWTPtWblv/0vld/wujXUSQsxP9JhnTk6O3mTwMpmsTK3xJPrpWnaSzYC+e/igqSyja2atuiyN829zAbC8MU9NNNZJCCkKooNnvXr1sGbNGp3nvv32W9SrV6/QlSIln/qkHnXZDPj1v7R7fgb0tFhrtFZtrYAO1fnW5fIwvht4+vv8jFyAxjoJIUVDdLft1KlT0a1bN7zzzjvo168fvL298fjxY2zduhWXLl3Cvn37zFBNUpzpS6I+vgkABtyTA1uu5B3/5xn/Gs20eQD/H2Tuf5OHlOs2Jx3jl7RYccCysLx7KL838eMnH1FSA0JIUTEqScLWrVsxefJkPHr0SHXM19cXS5cuRe/evU1aweLO0oPWlqYrUxAgPLYklE+IoD4bVsLxLcYFf+ruurWxAjZH8AFWPZmCjRVwYgAFSELKOkt/9hqVJKFPnz7o3bs3rl+/juTkZHh4eKBmzZqUtq8MUh/bVB9vVD+mDJwcAGWczGb5z7LNygX+fsL/W32XlKxcStROCLE8ozMMcRyHwMBAU9aFlECaO5Q08gUepQrLKAOfmC4OCQd8eY7PJGRrlRdAaUyTEFIciA6eS5YswYMHD7By5Uqtc//73/8QEBCAzz77zCSVI8WfriTqhZ3tagWgVx1gy2X+57e5fKYhD3sa0ySEFA+iZ9tu2LABderU0XmuXr162LBhQ6ErRUoWfxegR1BeUFOfaWtjRALIXPD5a9XVqyC8ByGEWJLoj7Z79+6hRo0aOs9Vq1YNiYmJha2TlhMnTmDIkCEIDAyEg4MDfH190bVrV/z1119aZePj49G2bVs4OjrC1dUVkZGRuHPnjs7rrly5EoGBgZBKpahcuTKioqKQlZVl8vqXBUlyPuAlyfNao8vD+Ek/+tZk5idTY7mwrpm5hBBiKaKDp42NDZ49e6bz3NOnT80yaWjNmjVITEzE2LFjcejQIaxYsQLPnj1DkyZNcOLECVW5a9euISQkBG/fvsVPP/2E9evX48aNG2jRogWeP38uuOaCBQswduxYREZG4rfffsOoUaOwcOFCjB492uT1L+2UM24nHuW/KwNojyB+GcmRfkBf3Z0VAk19+UTvAP+d1m4SQoor0UtV2rZtCwA4duyYznO5ubmCgGYKz549g5eXl+DYmzdvUK1aNdSpU0dVlw8//BAxMTG4ffu2auryvXv3UL16dYwfPx5LliwBACQnJ8PPzw8DBgzAt99+q7rmwoULMXPmTFy5cgVBQUEG1c3S06UtRX1t5/mHwuUkyr03f70JvMrkZ9muu6jdmtQktQY2duNbmcpgSWs3CSG6WPqzV3SH2meffYZOnTohJCQEo0aNgq+vLx48eIC1a9fi999/x6FDh0xeSc3ACQCOjo4ICgrC/fv3AQDZ2dk4cOAABgwYIHiQAQEBaN26Nfbu3asKnocPH4ZCocDgwYMF1xw8eDBmzJiBffv2GRw8yyLNtZ3RXYUzbm2sgNYb+OUoYmTm8IGzh9qjp6BJCCmORAfP9u3b47vvvsPEiRPx0UcfgeM4MMbg4uKC77//HuHh4eaopxa5XI74+Hi0adMGAHD79m1kZGQgODhYq2xwcDCOHj0KhUIBmUyGK1f4dDd169YVlPP29oanp6fqPNFNc23n30+A8Y2B15mAPBMYd0S4vZihqHuWEFJSGLXOc+jQofjoo49w+vRpvHjxAuXKlUOzZs20dloxp9GjRyMtLQ0zZswAwHfFAoC7u7tWWXd3dzDG8OrVK3h7eyM5ORlSqVRnfd3d3VXX0iUzMxOZmZmqn1NSyt5MFvW1nVLrvPWYhrC1Boa9k7fm003Kz6RVdtVSS5MQUhIYnSTBwcEB7dq1M2VdDDZr1ixs2bIFK1euRIMGDQTn8puwpH7O0HKaFi1ahKioKBG1LV2UY53RXfmAl5wBLMxnW7EmvoCnPb9dGMBvHeYqy9sEW5nST32PT0IIKe6MWIUHPH/+HNOmTUPTpk1Ro0YNJCQkAOB3Vbl48aJJK6gpKioK8+fPx4IFC/Dpp5+qjnt4eACAzlbjy5cvwXEcXF1dVWUVCgXS09N1ltXVelWaNm0a5HK56ks55lpaqC850XVOOat20H6+pVivfP5rOfvUAaY0z1uuYicBwHSn9COEkJJCdMvz7t27aN68OeRyOerVq4fbt2+rujH/+ecfxMbG4scffzR5RQE+cM6dOxdz587F9OnTBeeqVq0KOzs7XL58Wet1ly9fRrVq1SCTyQDkjXVevnwZjRs3VpV78uQJXrx4oTcJBABIpVK9+5mWdLqSvKt3o2qOdf56E/jqHJ82z8YKYLmAZu/ttitA10C+pbr3GhARCPg48a9TT+lHCCElieiW5+TJk+Hq6oqbN2/i999/h/pKl/fffx+nT582aQWVPv/8c8ydOxczZ87EnDlztM5LJBJ06dIFe/bsQWpqXnLVpKQkxMTEIDIyUnWsffv2kMlkiI6OFlwjOjoaHMehW7duZnkPxZ2+JO9KmpmDXmfmlc/SETgBIPYhEPuAb6luT+C/A3lJFDQDNCGElASiW57Hjx/HmjVr4OPjg5wc4cI9b29vwTZlprJ8+XLMnj0b7du3R6dOnRAbGys436RJEwB8y7Rhw4bo3Lkzpk6dCoVCgdmzZ8PT0xMTJ05UlXd3d8fMmTMxa9YsuLu7o127doiLi8PcuXMxbNiwMrtMxc9ZdwJ29TWdS0KB8Uf4Mt/9xU8Yyszhkxrk5vJ5aNUx8C1OzaBMqfYIISWZ6OCpUCj0jgmmpaXBysqoYdR8/fLLLwD49ZmHDx/WOq9s/QYGBuLkyZOYMmUKevToAYlEgjZt2mDZsmUoV66c4DUzZsyAk5MTVq1ahWXLlqFChQqYOnWqavZuWZMk51uFyi7YJaHAr7eAVxnAj3/zs2ml1vwSFOX+m9kMQC6/L2eHavyxrZeB7+Lzykit+a7a/depm5YQUnqIDp41a9bEsWPHEBYWpnXu999/z3e80FgnT540uGyDBg10Zj/SZcyYMRgzZoyRtSpd1Ltss3K1N68GdGcIUiZCULYiq7kLN7ee2DQvRR9lCyKElBaig+fw4cMxYcIE+Pj4oG/fvgCAt2/fYteuXVi9ejW++eYbk1eSmJ/62k31zafV2Vrx3bCa5+6pzcz108iSVa88/93fhYImIaT0EB08R40ahUuXLmH8+PGqccT3338fjDEMHz4cAwcONHklSdEY3wQA45MWDNyfl/jAGkCnmvzsWl1B9acE4IMaees+1Sl3Q1EfN6UgSggp6UQnhleKjY3FwYMH8fTpU3h6eqJz585o1qyZqetX7Fk6ObEpqC9RkUmACY35AHrqHvDtBaCAfO4A8lqrUmuA4/jAq1zuAuS/BIYQQsSy9GevqJZnRkYGqlWrhrVr16JLly6qWa6kZFMf71RkAwtP80G0fnnDAieQ1yLNzAGmNwc87PNambuuas+2peBJCCnJRAVPOzs7ZGRkFGkOW2J+6uOdSopsfo2modSXuHSozh9TrhNVvz7NtiWElAaixzxDQ0Nx7Ngx1W4mpOTzd+EzAG36BzhyW3utZkFsrIDNEcJ9ODW7aWm2LSGkNBEdPKdPn47u3btDJpMhMjIS3t7eWonU88sNS4oH9Qk8AL/GM0NXiiADcOBT7imTu+vqpqWkCISQ0kR08FTuYjJ37ly9u4toZh4ixUuSHGi7iR+flFoDE5voDpxWAHrXBXYmaLdGrZE3Hvo2VziOSd20hJDSTnTwnD17dr5bdpHi79ebeQkPMnOAK8/51qPmtOtcADm5wLKwvKQJEiughT8QUROYclx3gPR3oW5aQkjpZvRSFcKz9HRpY3z7V/57cGrSlTTBTpK3pycFSEJIUbP0Z6/pE9GSYq9eecAqn84DKwCdq+f9rMx3qy4jmw+cNJZJCCmLKHiWMbEPgH57+QTv+uSCD4jqG1hvjuATwEut847RWCYhpKwSPeZJSq7YB0DvPfkHTiXGtMctm/jxu6fQWCYhpKyj4FlGJMmBvgYGTgBYdwlws+ODpXqQpATvhBBC3bZlxvmHeduHGeJtDj+pqN1mPvASQgjJQ8GzlEmS80kKNAOenzO/HEUsZZIDQggheajbthRR3x1FffeSJDkwYJ/2Ok5NrSsBnzQA/n4CfHkub2cUmhhECCFCFDxLEfXdUdR3L1FPiqDJCvzsWjsJMC9EbWJQdZoYRAgh+lDwLEXU0+JJrYHbr4DFfwLXk/W/xsoKmNqUD5Y0MYgQQgxDwbOUUCZ6XxIKHLnD746y+kLBr8vO5ffepEBJCCGGo+BZCqiPdYolsaIxTUIIEYtm25YC6mOdYlhzwJYIanUSQohY1PIsBfyMyIlcrzzwTQcKnIQQYgxqeZZAmms5H6SIv8b09ylwEkKIsajlWcLoWsvZyFf3tmGA7n06ba0BH6eiqC0hhJRO1PIsYfSt5Vweprt8TQ++lam+xdjbHMoaRAghhUHBs4Rp5AvI/usvkFoDyen8bilnHuguf+cVn9x9SnPhFmM0w5YQQoxH3bYlzKNUIPu/bEGZOcDC09pllFmDAOBtLt/K7BGkvcUYIYQQ41DwLEGS5PxG1gXtjiKxBqw47dy0lDWIEEJMg4JnCaDMHnTrle5JQZre5gDTm/OZg6iVSQghpkfBs5hTn11rY+AItZ1EO1ctIYQQ06EJQ8VYkhxYFZc3uzYrF7Au4DUf1c7biowQQoh5UMuzmNKVr1ZixSdy18dOAoxuSIGTEELMjVqexZRmvlprLv/AaWMFRHelwEkIIUWBgmcxpdybE+ADZ04BM2yzco1L00cIIUQ8Cp7FlL8L35K0sSo4cAKU+IAQQooSjXkWU0lyYO+1gpemyCTAhMY0u5YQQooSBc9iqKDNrTkAK8L5wErrOAkhpOhR8CyGNCcLae6MwgA8eQN8/F4RV4wQQggAGvMslhr5AhIu72cGoHUljUIcCCGEWAgFz2LKSiM4KrIBm/+O2VrzO6UQQgixDAqexUCSHNh1lf8O8N22bzUmCp19AGT913erGVgJIYQULRrztLDYB/xOKVm5/HKT6K5A/BM+QObqWaKiUNsEmxBCSNGj4GlByi3GlMtRMrKBPnuE6zpDKwPH7wpfR2s6CSHEsih4WtD5h8J1nLoyCZWzB3Z059d8NvOj5SmEEFIcUPC0ID9nflatcnNra47vrlUPqBGBQBM//osQQkjxQMHTQpRjndlqLc23ucD094FXCuC+HOgfTEGTEEKKIwqeRSRJznfTKscq1cc6lewk/BIU6pIlhJDijYJnEVBPt2cnAcY3EQZOCQdMbk6BkxBCSgoKnmakbG0mZ+Sl28vIBsD4IJqRze+asjmCumcJIaQkoeBpJuqtTak1v/uJ4r+WZ4fqQL0K/Axa5YQgQgghJQcFTzNRT+6emQNMbw542PNjno9S88Y8d/9LLU9CCClpKD2fmTTy5VuZQF5rs0cQ/7P6ZKGsXP5nZWo+QgghxR+1PM3E3wU40i9vhq1yItC2K9qzbLNyKd0eIYSUJBQ8zUgZDM8/zDv27QXtcpRujxBCShYKnmakPmnIxgoY9g6Qo1Gmc3VgSnNqdRJCSElCY55m9OutvElDWbnA9xcBW7UnLrWmwEkIISVRmQ6eb968wbhx4+Dj4wOZTIb69etj+/btJrl2khxYdlZ4LDsXGPYuP/N2+vvAsf4UOAkhpCQq0922kZGRiIuLw+LFi1GjRg1s3boVvXv3Rm5uLvr06VOoa59/CLzV7KMF4CoFPn6vUJcmhBBiYWU2eB46dAhHjx5VBUwAaN26Ne7du4dJkyahV69esLa2NuraSXI+q5CttTCA2lrzS1YIIYSUbGU2eO7duxeOjo7o2bOn4PjgwYPRp08fnDt3Ds2aNRN9XfVJQjIJMPI9AAxws6PctYQQUlqU2eB55coV1KpVCxKJ8BEEBwerzhsTPNUzCymygWpueckRCCGElA5lNngmJyejSpUqWsfd3d1V53XJzMxEZmam6me5nE8NlJKSAgCo5QTYZgGKHEBmzf/83ylCCCEmovzMZYwVUNI8ymzwBACO40SfW7RoEaKiorSOV6xYUWf5OpOMqxshhJCCJScnw8Wl6MfDymzw9PDw0Nm6fPnyJYC8FqimadOmYcKECaqfX79+jYCAACQlJVnkF1hcpaSkoGLFirh//z6cnZ0tXZ1ihZ6NbvRcdKPnoptcLoe/v7/ez2pzK7PBs27duti2bRuys7MF456XL18GANSpU0fn66RSKaRSqdZxFxcX+g9bB2dnZ3ouetCz0Y2ei270XHSzsrJMuoIymyQhIiICb968we7duwXHN2zYAB8fHzRu3NhCNSOEEFLcldmWZ4cOHRAWFoaRI0ciJSUF1apVw7Zt23D48GFs3rzZ6DWehBBCSr8yGzwBYM+ePZgxYwZmz56Nly9fIjAwENu2bcNHH31k8DWkUinmzJmjsyu3LKPnoh89G93ouehGz0U3Sz8Xjllqni8hhBBSQpXZMU9CCCHEWBQ8CSGEEJEoeBJCCCEiUfA0kjn3Ai1uTpw4gSFDhiAwMBAODg7w9fVF165d8ddff2mVjY+PR9u2beHo6AhXV1dERkbizp07Oq+7cuVKBAYGQiqVonLlyoiKikJWVpa5347Z/PDDD+A4Do6OjlrnyuJz+fPPP9GxY0e4ubnBzs4O1atXx+effy4oU9aey8WLF9GtWzf4+PjA3t4egYGBmDdvHtLT0wXlSutzSU1NxeTJk9GuXTuUK1cOHMdh7ty5Osua4xk8e/YMgwYNgqenJ+zt7dG0aVMcP37cuDfDiFHCwsKYq6srW7t2LTtx4gQbNmwYA8C2bNli6aqZXI8ePVjr1q3Z6tWr2cmTJ9nOnTtZkyZNmEQiYcePH1eV+/fff5mTkxNr0aIFO3jwINu9ezerXbs28/HxYc+ePRNcc/78+YzjODZt2jQWExPDli5dymxtbdnw4cOL+u2ZxIMHD5iLiwvz8fFhDg4OgnNl8bls2bKFWVlZsY8++oj9/PPP7MSJE+z7779nUVFRqjJl7bkkJCQwmUzG6tWrx3bs2MGOHz/O5syZw6ytrdkHH3ygKlean8vdu3eZi4sLa9mypeozc86cOVrlzPEMFAoFq1OnDvPz82ObN29mR44cYV27dmUSiYSdPHlS9Huh4GmEgwcPMgBs69atguNhYWHMx8eHZWdnW6hm5vH06VOtY6mpqax8+fIsNDRUdaxnz57M09OTyeVy1bHExERmY2PDJk+erDr24sULJpPJ2IgRIwTXXLBgAeM4jiUkJJjhXZhX586dWZcuXdjAgQO1gmdZey4PHjxgDg4ObOTIkfmWK2vPZcaMGQwAu3XrluD4iBEjGAD28uVLxljpfi65ubksNzeXMcbY8+fP9QZPczyDVatWMQDszJkzqmNZWVksKCiINWrUSPR7oeBphGHDhjFHR0eWlZUlOL5161YGgJ0+fdpCNStarVu3ZjVq1GCM8f8R2tnZsY8//lirXLt27Vj16tVVP2/evJkBYGfPnhWUe/ToEQPAFixYYN6Km9imTZuYk5MTu3//vlbwLIvPZe7cuQwAS0xM1FumLD+X58+fC45PnjyZWVlZsTdv3pSp56IveJrrGbRt25bVrFlT65oLFy5kANiDBw9E1Z/GPI1gyF6gpZ1cLkd8fDxq164NALh9+zYyMjJUz0BdcHAwbt26BYVCASDv+dStW1dQztvbG56eniXq+T179gzjxo3D4sWL4efnp3W+LD6X33//He7u7rh27Rrq168PiUQCLy8vfPLJJ6ptpMricxk4cCBcXV0xcuRI3LlzB6mpqThw4AC+/fZbjB49Gg4ODmXyuWgy1zO4cuWK3msCQEJCgqh6UvA0QnJyss5M/gXtBVqajB49GmlpaZgxYwaAvPes77kwxvDq1StVWalUCgcHB51lS9LzGzVqFGrWrImRI0fqPF8Wn8vDhw+Rnp6Onj17olevXjh27BgmTZqEjRs3omPHjmCMlcnnUqlSJZw9exZXrlxB1apV4ezsjC5dumDgwIFYsWIFgLL534smcz0DU39ul+n0fIVhzF6gpcWsWbOwZcsWrFy5Eg0aNBCcM/S5lIbnt3v3bvzyyy+4ePFigXUuS88lNzcXCoUCc+bMwdSpUwEAISEhsLW1xbhx43D8+HHY29sDKFvPJTExEV26dEH58uWxa9culCtXDufOncP8+fPx5s0brFu3TlW2LD0XfczxDEz5vKjlaQRj9wItDaKiojB//nwsWLAAn376qeq4h4cHAN1/vb18+RIcx8HV1VVVVqFQaE3PV5YtCc/vzZs3GD16NP73v//Bx8cHr1+/xuvXr/H27VsA/D6vaWlpZe65AHn/LYSHhwuOd+jQAQC/BKEsPpepU6ciJSUFv/32G7p3746WLVti0qRJ+Prrr7F+/XqcOnWqTD4XTeZ6Bqb+3KbgaYS6devi33//RXZ2tuB4QXuBlnRRUVGYO3cu5s6di+nTpwvOVa1aFXZ2dqpnoO7y5cuoVq0aZDIZgLzxCc2yT548wYsXL0rE83vx4gWePn2K5cuXw83NTfW1bds2pKWlwc3NDX379i1zzwWAznElAGD/pdG2srIqk8/l0qVLCAoK0upibNiwIQCounPL2nPRZK5nULduXb3XBIz43BY1vYgwxhg7dOgQA8C2b98uON6+fftSuVSFMcbmzZvHALCZM2fqLfPhhx8yLy8vlpKSojp27949Zmtry6ZMmaI6lpyczGQyGfvkk08Er1+0aFGxnWKvKSMjg8XExGh9hYeHM5lMxmJiYtjly5cZY2XruTDG2G+//aZztueXX37JALA//viDMVb2nkvr1q1ZuXLlWGpqquD4d999xwCwffv2McbKznPJb6mKOZ7B6tWrGQAWGxurOpaVlcVq167NGjduLLr+FDyNFBYWxtzc3Nh3333HTpw4wYYPH84AsM2bN1u6aia3bNkyBoC1b9+enT17VutL6d9//2WOjo6sZcuW7NChQ2zPnj2sTp06+S5snj59Ojt58iT74osvmFQqLbaLuw2la51nWXwuXbp0YVKplH3++efs6NGjbNGiRUwmk7HOnTurypS157J//37GcRxr0qSJKknCggULmKOjIwsKCmKZmZmMsdL/XA4dOsR27tzJ1q9fzwCwnj17sp07d7KdO3eytLQ0xph5noFCoWC1a9dmFStWZFu2bGFHjx5lERERlCShqKWmprIxY8awChUqMFtbWxYcHMy2bdtm6WqZRatWrRgAvV/qLly4wEJDQ5m9vT1zdnZm3bp101oUrrRixQpWo0YNZmtry/z9/dmcOXPY27dvi+ItmY2u4MlY2Xsu6enpbMqUKaxixYpMIpEwf39/Nm3aNKZQKATlytpzOXHiBGvXrh2rUKECs7OzYzVq1GATJ05kL168EJQrzc8lICBA72fJ3bt3VeXM8QyePHnCBgwYwNzd3ZlMJmNNmjRhR48eNep90H6ehBBCiEg0YYgQQggRiYInIYQQIhIFT0IIIUQkCp6EEEKISBQ8CSGEEJEoeBJCCCEiUfAkJd6ZM2cwd+5cvH792tJVIYX0xx9/QCqV4t69e4Ljubm52Lx5M8LDw+Hl5QUbGxu4urqiSZMmWLZsGV68eCHqPllZWShfvjyaNGmit0xubi78/f1V6QaPHz8OR0dHPHz4UPwbI6UOBU9S4p05cwZRUVEUPEs4xhjGjRuH4cOHIyAgQHU8IyMD7du3x4ABA+Du7o7/+7//w/Hjx7F582a0adMGX3zxBSIiIkTdy8bGBv3798e5c+dw9epVnWWOHTuG+/fvY+jQoQCA0NBQNGrUSCuvMymjjEqtQEgx8sUXX2hlJ8lPenq6eStUjOXm5hbb96/MGX3t2jXB8REjRjAAbOvWrTpfl5aWxr777jvR97t69SoDwCZOnKjzfK9evZitra0g+8+uXbuYtbU1S0pKEn0/UrpQ8CQl2pw5c3Sm+YqJiWGM8anAOnXqxHbv3s3q16/PpFIpmzJlCrt79y4DwH788Ueta0JHsuobN26w3r17s3LlyjFbW1sWGBjIvvnmG4PqmJuby1atWsXq1avHZDIZc3V1Zd27d2e3b98WlGvVqhWrXbs2O3/+PHv//feZnZ0dq1y5Mlu0aBHLyckRlJXL5WzixImsUqVKzMbGhvn4+LCxY8eyN2/eaL2X0aNHszVr1rDAwEBmY2PD1qxZwxhj7I8//mBNmjRhUqmU+fj4sJkzZ7Lvv/9e8IfIkCFDmJubmyrnqLrWrVuzoKAgg56BIbp06cIaNmwoOPbo0SMmkUhYp06dRF3L0GfetGlT5uXlxbKysgTHX716xWQyGevZs6fgeGZmJnNxcWGzZs0SVR9S+lDwJCXa/fv32f/+9z8GgO3Zs0eVrF4ulzPG+ODp7e3NqlSpwtavX89iYmLY+fPnRQXPhIQE5uLiwurWrcs2btzIjhw5wiZOnMisrKzY3LlzC6zj8OHDmY2NDZs4cSI7fPgw27p1KwsMDGTly5dnT548UZVr1aoV8/DwYNWrV2dr165lR48eZaNGjWIA2IYNG1Tl0tLSWP369Zmnpyf78ssv2bFjx9iKFSuYi4sLa9OmDcvNzRW8F19fXxYcHMy2bt3KTpw4wa5cucL+/vtvJpPJWHBwMNu+fTv7+eefWceOHVmlSpUEwfPvv/9mANj3338veE8JCQkMAFu1apUhv6YCZWZmMjs7OzZ58mTB8S1btjAA7NtvvxV1PUOf+Q8//CDY0URp1apVDAA7fPiw1rU7dOjA3n33XVH1IaUPBU9S4uXXbRsQEMCsra3Z9evXBcfFBM/w8HDm5+enCshKn376KZPJZOzly5d663b27FkGgC1fvlxw/P79+1rBQpmA/9y5c4KyQUFBLDw8XPXzokWLmJWVFYuLixOU27VrFwPADh06JHgvLi4uWnXs2bMnc3BwYM+fP1cdy8nJYUFBQVrPslWrVqx+/fqC148cOZI5Oztrba9lrHPnzunc5m/x4sV6g1hWVpbgS0nMM09NTWWOjo7sgw8+EJRt0KABq1ixolaLnzHGZsyYwaysrLRa+aRsoQlDpNQLDg5GjRo1jHqtQqHA8ePHERERAXt7e2RnZ6u+OnbsCIVCgdjYWL2vP3DgADiOQ79+/QSvrVChAurVq4eTJ08KyleoUAGNGjXSqr/67NMDBw6gTp06qF+/vuCa4eHh4DhO65pt2rSBm5ub4NipU6fQpk0beHp6qo5ZWVnhww8/1HoPY8eOxaVLl3D69GkAQEpKCjZt2oSBAwfC0dEx3+dnqEePHgEAvLy8DCp/6dIl2NjYCL6UM27FPHNHR0d8+OGHOHToEJ4+fQqA35T6r7/+wqBBg2Blpf0R6eXlhdzcXDx58qSQ75qUZBQ8Sann7e1t9GuTk5ORnZ2NlStXan1Yd+zYEQDyXSbx9OlTMMZQvnx5rdfHxsZqvdbDw0PrGlKpFBkZGYJr/vPPP1rXc3JyAmNM65q63n9ycjLKly+vdVzXsa5du6JSpUpYtWoVACA6OhppaWkYPXq03vctlvL9yWQywXF/f38A0Fq6UrNmTcTFxSEuLg7Dhw8XnBP7zIcOHYrs7Gxs2rQJALB+/XpwHIfBgwfrrKuyjuq/E1L2SCxdAULMjeM4rWPKD8DMzEzB8eTkZMHPbm5usLa2Rv/+/fUGi8qVK+u9t6enJziOU61f1KTrWEE8PT1hZ2eH9evX6z2vTtf79/DwULW01OlqTVlZWWH06NGYPn06li9fjtWrVyM0NBQ1a9YUXXd9lHV++fKl4HhISAgkEgl+/vlnjBgxQnXczs4O7733HgC+pal5LTHPvFmzZqhVqxZ+/PFHjB07VrUERt/vVVlHzedMyhYKnqTEU34YimkJlC9fHjKZDP/884/g+P79+wU/29vbo3Xr1rh48SKCg4Nha2srqm6dO3fG4sWL8fDhQ51dosbo3LkzFi5cCA8Pj3wDd35atWqFQ4cO4cWLF6ogkJubi507d+osP2zYMMydOxd9+/bF9evXsWTJEqPrr0utWrUAALdv3xYc9/b2xpAhQ/Ddd99h+/bt+Oijjwq8ljHPfMiQIZg0aRJmzpyJ58+fY8iQIXrL3rlzBx4eHjpb6aTsoOBJSry6desCAFasWIGBAwfCxsYGNWvWhJOTk97XKMfE1q9fj6pVq6JevXo4f/48tm7dqlV2xYoVeP/999GiRQuMHDkSlSpVQmpqKm7duoVffvkFJ06c0Huf5s2bY8SIERg8eDAuXLiAli1bwsHBAY8fP8aff/6JunXrYuTIkaLe77hx47B79260bNkS48ePR3BwMHJzc5GUlIQjR45g4sSJaNy4cb7XmDFjBn755ReEhoZixowZsLOzw9q1a5GWlgYAWmN9rq6uGDBgANasWYOAgAB06dJFVJ0L4ufnhypVqiA2NhZjxowRnPv6669x9+5d9O3bFz///DO6du0KHx8fpKen49q1a9i+fTtkMhlsbGwAGPfMBwwYgOnTp+OLL76Aq6srIiMj9dY1NjYWrVq10tmiJ2WIhScsEWIS06ZNYz4+PszKykrnOk9d5HI5GzZsGCtfvjxzcHBgXbp0YYmJiTrXed69e5cNGTKE+fr6MhsbG1auXDnWrFkzNn/+fIPqt379eta4cWPm4ODA7OzsWNWqVdmAAQPYhQsXVGWU6zw1DRw4kAUEBAiOvXnzhs2cOZPVrFmT2draqpbSjB8/XrAUA/+t89Tljz/+YI0bN2ZSqZRVqFCBTZo0iS1ZsoQBYK9fv9Yqf/LkSQaALV682KD3LNasWbOYm5sbUygUWudycnLYxo0bWVhYGPP09GQSiYS5uLiwRo0asVmzZrEHDx5ovcaQZ64uIiKCAWCjRo3SW8dbt24xAGz37t3Gv1FSKnCMMWbB2E0IKUbatWuHxMRE3LhxQ+vcxIkTsWbNGty/f1/nxKbCevToESpXroyNGzeiV69eJr++KcyaNQsbN27E7du3IZFQx11ZRr99QsqoCRMm4J133kHFihXx8uVLbNmyBUePHsW6desE5WJjY3Hjxg2sXr0aH3/8sVkCJwD4+Phg3LhxWLBgAXr27KlzmYglvX79GqtWrcLKlSspcBIKnoSUVTk5OZg9ezaePHkCjuMQFBSETZs2oV+/foJyTZs2hb29PTp37oz58+ebtU4zZ86Evb09Hj58iIoVK5r1XmLdvXsX06ZNQ58+fSxdFVIMULctIYQQIlLx6hchhBBCSgAKnoQQQohIFDwJIYQQkSh4EkIIISJR8CSEEEJEouBJCCGEiETBkxBCCBGJgichhBAiEgVPQgghRKT/Bx5ba4+IfrC8AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAekAAAGxCAYAAABGE04NAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB+q0lEQVR4nO3deVwU5R8H8M8Cyy0g4LUg3orKZV5oeaDifaGReeRtZZaapnkfpWmWlfWzNFOxFDTv20RFK/NAPBKPUjwQNQ9QDgXkeH5/jLO7szu7Owu77ALf9+vFC3f2mZlnH9b5znOOjDHGQAghhBCrY2PpDBBCCCFEHAVpQgghxEpRkCaEEEKsFAVpQgghxEpRkCaEEEKsFAVpQgghxEpRkCaEEEKsFAVpQgghxErZWToDZVlhYSHu3buHChUqQCaTWTo7hBBCjMQYQ2ZmJhQKBWxsSr5eS0HajO7du4fq1atbOhuEEEKK6c6dO/D19S3x81KQNqMKFSoAAG7evAlPT08L58Z65eXl4eDBg+jcuTPkcrmls2O1qJykoXKShspJmrS0NNSqVUt5PS9pFKTNiG/irlChAtzc3CycG+uVl5cHZ2dnuLm50cVCDyonaaicpKFykiYvLw8ALNZlSQPHCCGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQbqUiIqKgkwmg0wmw9GjR7XeZ4yhbt26kMlkaN++veA9mUyG999/X+/x27dvrzy+TCaDk5MTgoOD8c0336CwsNCEn4RI1bFjR7z77ruCbRs2bEDfvn3h4+MDmUyG4cOH69z/xo0b6NevHzw8PODq6orw8HCcPXu2yPk5cuQIRo4cCX9/f7i4uMDHxwd9+vRBQkKCaPqzZ8+iU6dOcHV1hYeHB/r164cbN24U+fzbtm3DwIEDUbduXTg5OaFmzZoYPHgwrl27Jpr+8OHDaNWqFZydneHt7Y3hw4fj4cOHRT7/pUuX8N5776FVq1ZwcXHR+X+Rt3HjRoSEhMDR0REKhQITJ05EVlaWIM3q1avh4+ODZ8+eFTlfpGyjIF3KVKhQAatXr9bafuzYMSQlJaFChQpFPnbt2rVx4sQJnDhxAps2bYKPjw8+/PBDTJ8+vThZJkWwc+dOHD9+HLNnzxZs37VrF1JTU9G7d2/Y29vr3P/Ro0do06YN/v33X6xZswa//vorcnJy0L59e/zzzz9FytMPP/yAW7duYcKECdi3bx+WLVuGhw8fIjQ0FEeOHBGkvXr1Ktq3b48XL17g119/xZo1a/Dvv/+iTZs2ePToUZHO//nnn+P58+eYOXMmDhw4gAULFuDcuXN45ZVXcOnSJUHaxMRE9OrVC1WqVMHOnTuxbNkyHDp0CB07dkRubm6Rzn/mzBns2LEDnp6e6Nixo960GzZswMCBA9G8eXPs378fc+fORVRUFPr16ydIN2zYMLi4uGDJkiVFyhMpBxgxm/T0dAaAPX78uNjHWrt2LQPARo8ezZycnFh6errg/SFDhrBWrVqxxo0bs3bt2gneA8DGjRun9/jt2rVjjRs3Fmx78eIFq127NnN2dmYvXrwwKr/Pnj2TnPbFixdsx44dRp9Dff+8vLwi7Wspz58/1/t+ixYt2JtvvinY9uLFC7Zt2zZlObm4uLBhw4aJ7j9lyhQml8vZrVu3lNvS09OZt7c3e+ONN4qU5wcPHmhty8zMZFWqVGEdO3YUbI+MjGTe3t6C7+mtW7eYXC5nU6dONdn57969y+RyORs1apRy24sXL1i9evVYw4YNBd+L48ePMwDs+++/L9L5CwoKlP/evHkzA8Di4uK00uXn57Nq1aqxzp07C7Zv2LCBAWD79u0TbP/yyy+Zu7u7Uf9nTKG4/+/Ki8ePHzMAWtfckkI16VJm4MCBAICYmBjltvT0dGzduhUjR4406bnkcjmaNm2K58+f6639zJs3DzKZDGfPnsXrr7+OihUrok6dOgC4Zvjvv/8eISEhcHJyQsWKFfH6669rNXsyxrB48WLUqFEDjo6OaNasGWJjY9G+fXtB8/3Ro0chk8nwyy+/YPLkyfDx8YGDgwOuX78OAMrakpubG5ydnfHqq6/i8OHDgnM9evQIb7/9NqpXrw4HBwdUqlQJr776Kg4dOqRMc+7cOfTs2ROVK1eGg4MDFAoFevTogZSUFGWanJwcTJ8+HbVq1YK9vT18fHwwbtw4PH36VHC+mjVromfPnti2bRuaNGkCR0dHzJ8/X2d5njt3DqdPn8Zbb72l9Z6NjbT/stu3b0eHDh1Qo0YN5TY3Nzf069cPu3fvRn5+vqTjqKtcubLWNldXVzRq1Ah37txRbsvPz8eePXvQv39/uLm5KbfXqFEDYWFh2L59u9Hn1nV+hUIBX19fwfnv3r2La9euYfDgwbCzs1Nub926NerXr1/k80st+5MnT+L+/fsYMWKEYHtkZCRcXV21zj948GBkZGRg48aNRcoXKdsoSJcybm5ueP3117FmzRrltpiYGNjY2GDAgAEmP19SUhLs7OxQsWJFg2n79euHunXrYvPmzVixYgUA4J133sHEiRPRqVMn7NixA99//z0uXbqE1q1b48GDB8p9169fjzlz5qBr167YuXMn3n33XYwePRr//vuv6LmmT5+O5ORkrFixArt370blypWxfv16dO7cGW5ubli3bh1+/fVXeHp6okuXLoJA/dZbb2HHjh2YM2cODh48iJ9++gmdOnVCamoqAODZs2cIDw/HgwcPsHz5csTGxuKbb76Bn58fMjMzAXA3FX379sWXX36Jt956C3v37sWkSZOwbt06dOjQQatJ9ezZs5gyZQrGjx+PAwcOoH///jrLcc+ePbC1tUXbtm0NlrmY7OxsJCUlISgoSOu9oKAgZGdnF6tvWF16ejrOnj2Lxo0bK7clJSUhOztb5/mvX7+OnJwck5z/xo0buH37tuD8fNN3YGCg6PkTExNNcm5d+ONrfn65XA5/f3+t81etWhX+/v7Yu3evWfNFSic7w0mItRk5ciTCwsJw6dIlNG7cGGvWrEFkZGSx+qN5fA3r0aNH+Pbbb3H27FlERkbCycnJ4L7Dhg0T1BBPnjyJVatWYenSpZg0aZJye5s2bVC/fn189dVX+Pzzz/HkyRPs2rULkZGRWLlypTJdQEAAWrVqhfr162udq06dOti8ebPy9fPnzzFhwgT07NlTUFPp3r07XnnlFcyYMQOnTp0CABw/fhyjR4/GmDFjlOn69Omj/PfVq1eRmpqK1atXC7a/8cYbyn8fPHgQv/32G5YsWYIpU6YAAMLDw1G9enUMGDAAP//8s+D4Dx8+xOXLl0U/i6YTJ06gXr16cHV1NZhWzJMnT8AYg6enp9Z7/Db+hqS4xo0bh2fPnmHmzJnKbfyxdZ2fMYYnT56gWrVqxTp3fn4+Ro0aBVdXV3z44YeSz2+qz66LofPfunVLa/srr7wiaMkhhEc16VKoXbt2qFOnDtasWYOLFy8iPj7eJE3dly5dglwuh1wuh0KhwNKlSzF48GCsWrVK0v6atcM9e/ZAJpNhyJAhyM/PV/5UrVoVwcHBypGxp06dQl5eHl5//XXB/qGhoahZs6akc/31119IS0vDsGHDBOcqLCxE165dER8frxxB26JFC0RFRWHBggU4efIk8vLyBMeqW7cuKlasiI8//hgrVqzA5cuXtc7PD5TSHF0dGRkJFxcXrSb2oKAgSQEaAO7duyfatGssmUxWpPekmj17NjZs2ICvv/4aTZs2LdHzM8YwatQo/PHHH/j5559RvXp1yecwxWeXwpjzV65cGQ8fPixSNwQp26gmXQrJZDKMGDEC3377LXJyclC/fn20adOm2MetU6cONm7cCJlMBkdHR9SqVQvOzs6S99esGT148ACMMVSpUkU0fe3atQGoah5i6XTtK3YuAFqBXl1aWhpcXFywadMmLFiwAD/99BNmz54NV1dXREREYMmSJahatSrc3d1x7NgxLFy4EDNmzFDW+saMGYNZs2ZBLpcjNTUVdnZ2qFSpkuAcMpkMVatW1aqtGVNrzM7O1vm5pahYsSJkMplojTEtLQ2AeC3PGPPnz8eCBQuwcOFCrel9Xl5eAMRr62lpaZDJZPDw8CjyuRljGD16NNavX49169YJWjuknL+4n90Q9fNr/h11nd/R0RGMMeTk5BS5BYWUTRSkS6nhw4djzpw5WLFiBRYuXGiSY/IDtopKs4bg7e0NmUyGP/74Aw4ODlrp+W38RU29j5r333//idamxc4FAN999x1CQ0NF88dfML29vfHNN9/gm2++QXJyMnbt2oVp06bh4cOHOHDgAACuP3Pjxo1gjOHvv/9GVFQUPvnkEzg5OWHatGnw8vJCfn4+Hj16JAjUjDH8999/aN68ud786uPt7a0MpkXh5OSEunXr4uLFi1rvXbx4EU5OTsobpKKYP38+5s2bh3nz5mHGjBla79epUwdOTk46z1+3bl04OjoW6dx8gF67di1Wr16NIUOGaKXh+6cTExPRu3dvrfMHBAQU6dxS8X3hFy9eRKNGjZTb8/PzcfXqVeXgT3VpaWlwcHCgAE20lNnm7qysLEycOBEKhQKOjo4ICQmRPHoyLi4O4eHhqFy5MlxdXREUFIRvv/0WBQUFZs61dD4+PpgyZQp69eqFYcOGWTo7onr27AnGGO7evYtmzZpp/fAXsxYtWkAulwv6mAGuT/v27duSzvXqq6/Cw8MDly9fFj1Xs2bNROcV+/n54f3339e50IdMJkNwcDC+/vpreHh4KNPw82TXr18vSL9161Y8e/bM4Dxaffz9/Ys9sCsiIgJHjhwRjHrOzMzEtm3b0Lt3b8GoZ2N8+umnmDdvHmbNmoW5c+eKprGzs0OvXr2wbds25UA7AEhOTkZcXJzWXGGpGGMYM2YM1q5di5UrV2qNnub5+PigXr16iI6OFvyfPXnyJP75558in1+qli1bolq1aoiKihJs37JlC7KyskTPf+PGDUFAJ4RXZmvS/fr1Q3x8PBYvXoz69esjOjoaAwcORGFhIQYNGqRzv0OHDqFLly5o27YtVq1aBRcXF+zatQsTJkxAUlISli1bVoKfQr/FixdLTpuUlIQtW7ZobW/UqJHZLg6vvvoq3n77bYwYMQJnzpxB27Zt4eLigvv37+PPP/9EYGAgxo4dC09PT/Tu3RtbtmzBu+++i4iICKSkpGD+/PmoVq2apKkvrq6u+O677zBs2DCkpaXh9ddfR+XKlfHo0SNcuHABjx49wg8//ID09HSEhYVh0KBB8Pf3R4UKFRAfH48DBw4oL5579uzB999/j759+6J27dpgjGHbtm14+vQpwsPDAXCDxLp06YKPP/4YGRkZePXVV/H3339j7ty5aNKkiej0Kanat2+vXPxDsx87MTER2dnZsLOzQ0FBAW7fvq38u7Zr105Zq//oo4/wyy+/oEePHvjkk0/g4OCAxYsXIycnB/PmzRMcc/jw4Vi3bh1u3rypcwwAACxdulQ5Ar9Hjx44efKk4H31Foz58+ejefPm6NmzJ6ZNm4acnBzMmTMH3t7emDx5stbnPXbsGBhjestl/PjxWL16NUaOHInAwEDB+R0cHNCkSRPl66FDh2L+/PmIjIzEe++9h4cPH2LatGkICAjQCu78ZxYb0KXu+fPn2LdvHwAoz33s2DE8fvwYLi4u6NatGwDA1tYWS5YswVtvvYV33nkHAwcOxLVr1zB16lSEh4eja9euguMWFhbi9OnTGDVqlN7zk3Kq5Kdmm9/evXsZABYdHS3YHh4ezhQKBcvPz9e57+DBg5mDgwPLysoSbO/cuTNzc3MzKh/mWMwkPj5ebzpdi5no+pk7dy5jTHwxE6nmzp3LALBHjx6Jvr9mzRrWsmVL5uLiwpycnFidOnXY0KFD2ZkzZxhj3KIK27dvZ/Pnz2e+vr7M3t6eBQUFsT179rDg4GAWERGhPFZcXBwDwDZv3ix6rmPHjrEePXowT09PJpfLmY+PD+vRo4cyfU5ODnv33XdZUFAQc3NzY05OTqxBgwZs7ty5ysUkrl69ygYOHMjq1KnDnJycmLu7O2vRogWLiooSnCs7O5t9/PHHrEaNGkwul7Nq1aqxsWPHsidPngjS1ahRg/Xo0UNyeaanpzNXV1e2ZMkSwfYXL16wxo0b6/xbai6scf36dda3b1/m5ubGnJ2dWceOHVlCQoLW+fr378+cnJy08q2pXbt2er9Lms6cOcM6duzInJ2dmZubG+vbty+7fv26VrqmTZuyqlWrGiyXGjVq6Dx3jRo1BOW0Y8cOtm/fPhYaGsocHR2Zp6cnGzp0qOiCKN7e3iw0NNTg+W/evCnp/Lzo6GgWFBTE7O3tWdWqVdn48eNZZmamVrrDhw8zAKJ/G3OixUyksfRiJmUySI8ePZq5urpqrUIVHR3NALDjx4/r3Hf48OGsQoUKgtWFGGPsjTfeYJUrVzYqH6YM0mWZrovFjRs3mL29PVu4cKGFcmY577//PmvYsCErLCxUbjPXRbVKlSrso48+MukxpcrIyGB2dnbsf//7n8mOaUw5Xbp0iQFge/bsMdn5jTVkyBDWunXrEj8vBWlpLB2ky2Rzd2JiIho2bKjV78YvLpCYmIjWrVuL7vvuu+8iJiYG48ePx4wZM+Ds7Izdu3dj+/btWLRokd7z5ubmChaxyMjIAADk5eVpTfMhKnl5ebh58yamTZuGV199FW5ubvj333+xdOlSuLm5YejQoeWu/D7++GP8/PPP2LRpk3K6GV8GpiyLS5cu4fnz55g0aZJFyvjIkSPw8fHB8OHDTXZ+Y8rp0KFDCA0NRefOnS3y+ZOSkrBp0yYcPHiwxM9vju9TWWTp8pExZqAjqBSqX78+ateurRypy7t//z4UCgU+++wzvQ+N+OuvvxAZGYl79+4B4PqYFi1apFy0Qpd58+aJLvcYHR1t1FSm8uj+/fvKBzg8e/YMzs7OCAgIwJAhQ+Dj42Pp7FlEfHw8srKyEBYWZumsEDO5ePEi7t27hy5dulg6K0SH58+fY9CgQUhPTxcsc1tSymRNGij6QgoJCQmIiIhAy5YtsXLlSri4uODIkSOYNWsWcnJytJ5KpG769OmClbUyMjJQvXp1hIWFKacZEW15eXmIjY3FiRMnIJfLLZ0dq9G9e3fBa76cwsPDqZz0KE3lpPk3LkmlqZwsydwr1BlSJoO0l5dXkRdyGDduHKpUqYLt27fD1tYWABAWFgYbGxvMmzcPgwcP1jnH1MHBQXQ+ML+KF9GPykkaKidpqJykoXLSz9JlUybnSQcGBuLKlStaS+zxiyvoW8zg/PnzaNq0qTJA85o3b47CwkJcuXLF9BkmhBBCRJTJIB0REYGsrCxs3bpVsH3dunVQKBRo2bKlzn0VCgXOnDmjtXDJiRMnAAC+vr6mzzAhhBAiokw2d3fr1g3h4eEYO3YsMjIyULduXcTExODAgQNYv369spY8atQorFu3DklJScrn7n744YcYP348evXqhXfeeQfOzs44fPgwli5dik6dOiE4ONiSH40QQkg5UiaDNABs27YNM2fOxJw5c5CWlgZ/f3/ExMTgzTffVKYpKChAQUGBYKWjDz74AD4+Pvj6668xevRoZGdno2bNmpg7d67gcXiEEEKIuZXZIO3q6oply5bpXcYzKipKa31dgFtS1Nzr+xJCCCGGlMk+aUIIIaQsoCBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBNCCGEWCkK0oQQQoiVoiBtBZLTgS2Xud+EEEIIz87SGSjvktOBzuuB7HzAyQ44OATwc7d0rgghhFgDqklb2Om7XIAGuN+n71o2P4QQQqwHBWkL83UD5C//Ck52QAsfy+aHEEKI9aAgbUHJ6cDwnUBeIReoo/pQUzchhBAVCtIWpN7UnVcIpGRYNj+EEEKsCwVpC2rhwzVxA9TUTQghRBuN7rYgP3duNPfpu1yApqZuQggh6ihIW5ifOwVnQggh4qi5mxBCCLFSVJO2MsnpquZvgJrCCSGkPKMgbUXUVx9ztAMYA3ILaCUyQggpr6i524qoT8nKyecCNEArkRFCSHlFQdqKqE/JcrQDHGy5f9P0LEIIKZ+ouduKaE7JAqhPmhBCyjMK0lZGc0oWBWdCCCm/qLmbEEIIsVIUpAkhhBArRUGaEEIIsVIUpAkhhBArRUGaEEIIsVIUpAkhhBArRUGaEEIIsVIUpAkhhBArRUGaEEIIsVIUpAkhhBArVWaDdFZWFiZOnAiFQgFHR0eEhIRg48aNRh9n1qxZkMlkCAgIMEMuCSGEEN3K7Nrd/fr1Q3x8PBYvXoz69esjOjoaAwcORGFhIQYNGiTpGOfPn8eXX36JKlWqmDm3hBBCiLYyGaT37duH2NhYZWAGgLCwMNy+fRtTpkzBgAEDYGtrq/cY+fn5GDFiBN555x1cuHABjx8/LomsE0IIIUplsrl7+/btcHV1RWRkpGD7iBEjcO/ePZw6dcrgMRYvXoy0tDQsXLjQXNkkhBBC9CqTNenExEQ0bNgQdnbCjxcUFKR8v3Xr1jr3v3z5MhYsWIBt27bB1dVV8nlzc3ORm5urfJ2RkQEAyMvLQ15enlb6+HvAnn+BnvWB5gpuW0oGkHAPaKoAfN0kn7pU48tGrIyICpWTNFRO0lA5SWPp8imTQTo1NRW1a9fW2u7p6al8X5fCwkKMHDkS/fr1Q/fu3Y0676JFizB//nyt7XFxcXB2dhbdpyWAR+eBfedV2+QA/r4B/G3U2Uu/2NhYS2ehVKBykobKSRoqJ/2eP39u0fOXySANADKZrEjvffXVV7h27Rp27dpl9DmnT5+OSZMmKV9nZGSgevXqCAsLg5eXlyDt3KPA1iuq1/0bAq9UBWbGqbYtDAP6+BudjVInLy8PsbGxCA8Ph1wut3R2rBaVkzRUTtJQOUmjr1JXEspkkPby8hIt2LS0NACqGrWm5ORkzJkzB4sXL4a9vT2ePn0KgBtEVlhYiKdPn8LBwQFOTk6i+zs4OMDBwUFru1wu1/pP0LshsOGy8LWiAiCzBbLzASc7oIUfUJ7+74iVE9FG5SQNlZM0VE76WbpsyuTAscDAQFy5cgX5+fmC7RcvXgQAnXOeb9y4gezsbEyYMAEVK1ZU/hw/fhxXrlxBxYoVMX36dJPkMdQX2NQfeLMx9zvUF/BzBw4OAZaGc7/93E1yKkIIIaVUmaxJR0REYNWqVdi6dSsGDBig3L5u3TooFAq0bNlSdL+QkBDExcVpbZ84cSLS09Oxdu1a+Pr6miyfob7cjzo/dwrOhBBCOGUySHfr1g3h4eEYO3YsMjIyULduXcTExODAgQNYv369co70qFGjsG7dOiQlJaFGjRrw8PBA+/bttY7n4eGB/Px80fcIIYQQcymTQRoAtm3bhpkzZ2LOnDlIS0uDv78/YmJi8OabbyrTFBQUoKCgAIwxC+aUEEIIEVdmg7SrqyuWLVuGZcuW6UwTFRWFqKgog8c6evSo6TJGCCGESFQmB44RQgghZQEFaSuUnA5sucz9JoT/PqRkWDonhJCSVmabu0ur5HSg83rVXGmailW+qX8fPOTAJz6WzhEhpCRRTdrKnL7LXZAB7vfpu5bND7Es9e9DTr7+tISQsoeCtJVp4cPVoIGXq45RzalcU/8+OFK7FyHlDv23tzL8qmOn73IXaGrqLt/Uvw9NqwB//2npHBFCShIFaStEq44Rdfz3IS+v/D0ZjZDyjpq7CSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQIkAPeCHEetBiJoQQJXrACyHWhWrShBAlesALIdaFgrSVoibH0qcs/M3oAS+EWBdq7rZCRW1yTE6nB3NYSllpJqYHvBBiXYpdk167di2ysrJMkRfyUlGaHPkgMTmW+12aa3OlUVlqJvZzB15vRAGaEGtQ7CA9evRoXL58WbDtypUrxT1suVaUJseyFCRKI2omJoSYQ7GDNGNM8LqgoAABAQE4e/asYHtKSgrWrl1b3NOVC3yT49Jw6c2mFCQsqyh/M0IIMcQsA8c0AzcA3L17F6NHjzbH6coUfvARYFyTY1kLEqVxEBY1ExNCTI0GjlmR4g4+8nMvGwGitAzCooF6hBBzoyBtRcT6lcvjxb80lENpuZHg0Q0FIaWTSYL0o0ePTHGYco/vV+Yv/OW1X7k0lENpuJHglbYbCl3oRoOURyYJ0r1794abmxuCg4MRGBgImUyGlJQUNGnSBDKZzBSnKBfK6xzVlAzV71pepaMcSsONBE/zhmJ5PDCuuXWWqy5l5UaDEGMVO0jv3bsXZ8+exdmzZ5GQkIDff/8dABAREQFHR0cEBAQgODgYLi4uxc5seVBW+pWlSk4H+m4CPvHhfu8erCoDay6H0nAjwVO/oQCAjZeAnf8YH+gsWZMtTS0XhJhSsYN0t27d0K1bN+XrtLQ0ZcDmg/fq1avBGKNaNdFy+i6Q8/Lim1PKLr7WfiPB428olsdzARowPtBZuiZbmlouCDElkw8c8/T0RKdOndCpUyfltoyMDCQkJODcuXOmPh0p5Vr4AI4vv4WOdPE1Gz93rol75z9FC3QlUZPVV1O3ppYL6hsnJanYQXrZsmXo378/fH19daZxc3NDWFgYwsLCins6Usb4uQM7BgB//8n9poue+RQn0BlbkzU2kEmpqVtDy4WlWxRI+VPsxUw+++wz1KhRAy1btsQXX3yBpKQkU+SLlCO+bsLfxHyKuuCKMYvlFGUd+dKyrG1pyScpO4odpO/fv49Dhw6hefPm+Oabb1C/fn2EhIRgwYIFWmt6E/MpjSt0lWel8e8lNcAXJZCVlmVtS0s+SdlR7OZuGxsbZVP2//73Pxw/fhxbtmzBTz/9hLlz56J+/fro378/+vfvjyZNmpgiz0QDNcGVLtb69zJVX2tRBnlZU5+zPqUln6TsMPna3a+++iq+/vpr3Lp1CydOnEDv3r2xadMmNGvWDLVr18bUqVNNfcpyj5rgSoapar/W+Pcq7qNO1cumqOvIl5a1z0tLPknZYJYHbPBatGiBzz//HNeuXUNCQgIGDx6MvXv3mvOU5RI1wZmfKZ/XbY1/r+LcOIiVTVkMZKWxi4KUfmZbu7uwsBA5OTlwdnYGAISEhCAkJASffvqpuU5ZblETnPGMbdo15RSkovy9ktOB08mAvGinNEisiVpqGZWHhUastYuClH0mC9I5OTnYuHEj9u7di+PHj+Phw4dgjMHBwQGNGjVChw4dMHjwYAQHB5vqlIQUSVEuuKZeTMOY6UR8flkBsKSGavlUU9K8cQCkl5GlFxoxVV+6vuOUhxsRYp2KHaSzs7OxZMkSLFu2DOnp6fD390fHjh1RuXJlODo6Ii0tDTdu3MCqVauwdOlStG7dGkuWLEGrVq1MkX8CusvnmbPmZ8nWCj6/ji8X7Eu4Z/ogDQhvHLZcll5Gliwbse9+NWfj9j99l5v+N3yn7v9Dlr4RIeVXsYN0vXr14OLiglmzZmHw4MGoUqWKaDrGGOLi4rB27VrlSPDRo0cX9/QEdJcPGHejUtQLrqUW0+Dzywq4100VJXdOqWUkpWzMsVKX2He/Tz1p+6p/Z+Q2QF6h8DjqeaQuJWIpxQ7Sn3zyCYYNGwZbW1u96WQyGTp06IAOHTpg/vz5SE5OLu6pyUumvssvjcseGnOjon7B9XVTDZIypvm5JMtHmd9kADeKt+iL1LybOiiZorVHLO/F+e6rf2fyClWBWtdxrGHFM1L+FDtIjxw50uh9ateujdq1axf31OQlU15QzdV0XpzAJmVfqRdr9WO18DH+s1qqa8HPHajmD+y7IS29WJlp5v3zjsBfKUCEPxAqsqqvKYNScVt7dJW72Hc/L0/aMTW/M1F9uP7+0nRzSsq+YgfptWvXIjIyEq6urqbIDykiU11Qjb2YSgmgxQlsUveVcqOieawPWxofOIoSbEq65q2rzDTzPv437t8bLwGb+osHalMpbmuPvnIv6ne/tDZhl8aWLlJ0xQ7So0ePRuPGjdGiRQvltitXrqBhw4bFPTSxAGMuplIDaHFqUUXdV+xCpnksyIwPHEV50IR6v+f6CPMGQ0B3mWk+V1rd9qvmzVdxA6K5Bm6VtiZsGiRa/hQ7SDPGBK8LCgoQEBCA+Ph4vPLKK8rtKSkpiI2NxYgRI4p7SmJGxlxM9QVQzWblol5gjWnG5i9eDraATMY9n1r9QqZ5rG51uR9jAoexwUaz33PIduDIUPNeWHWVGZ939edK8yL8zZcfXnECYmmt9ZoaDRItf8yymIlm4AaAu3fvYvTo0RSkSwGpF1NdwUDsbl9zDu6Wy/qbpdUvxlL2Vb945RaotqtfyMQGjLXw4VbGkkI9X+r76Gt+bOEjHDmcV2j+C6u+gObnLnyutK0M+Lqz+Wv3plDaar3mQFPByh+zrThGyj5dwUDsbp9fIlIzgEf1AVKeqFbS0jdASOw9/ny+bqqLl2ZNWnP1LFMOGDPU/OjnzjVxD9muf+SwroFeRa056gtopqqVUt9oyaMWhfKHgjQpFrFgoO9uXzOAD9kO2DLVSloJD3Q352nuu/8a8PUp7dG5vm7Ahf8AyLjmbMB8A8akND+G+nJN3LourMnpQPh67qbC0Q6IHaKdZ2MX6TCkuPOaqW/UcqhFoXwxSZB+9OiRKQ5DzMAStR19d/vqAZxvBrZVW0mrhZ/uAK8Z/CETBkh++ox68OD7nNXTPckVnxNrqNlaLF9Smx/1XVj3X+cCNMD93n8d8HIybpEOXXlPTuduZvgbFlONqjdn3yjV0K0L/T0syyRBunfv3nBzc0NwcDACAwMhk8mQkpKCJk2aQCaTmeIUpAgsWdvRFZQ0+4WH7xSupGWoP1Wzf/rrk8IAKRY81AOpgy2w9rxq8YqoPoab0sX6xtWn/xS7+VFzCAcr+ih79RHkyelAp19UffRfneRq6VLyaCgIm6tvlGro1oX+HpZX7CC9d+9enD17FmfPnkVCQgJ+//13AEBERAQcHR0REBCA4OBguLi4FDuzRD/NO15rHQmqHsCj+gC7rgB4rlpJSyzA6xq0JRYgNYOHeiBNzQY++5NLl1fI1b4B/U3pDrbA5FZcTVRskFlxmx+71QOWnuSCqYMt99qYRTp0jSA/fVc4iC7HiO+AoSBsrr5RU35n+e4TqgEWnbVeQ8qTYgfpbt26oVu3bsrXaWlpyoDNB+/Vq1eDMUa1ajMSu+O19pGgyemqmnTLGrqf7qTvbl4zQPLBY/91CGqo6oPPNGvfAPfbwVYVKNWb0nMLuMD+9UnT1CTUH+rAN9Efeku8li41oNrJgPyXnzevkLvJ6FZP9ZkArpYtdUlRfUFY1w2TKZjyO9t3E/A0j2qAxWHt15DywOQDxzw9PdGpUyd06tRJuS0jIwMJCQk4d+6cqU9HXtLVzPthKACmqp2VBGOfRmXo6U5in43fruscfCD++pR2UNcVfPh7SJkMCK6ivfCHKWoS6jccPD6IFCfg2cgguCn54gQQXJUL/tEXgZ/OccF7+E79AUvzbyfWv/3VKe056KZiyhp6DtUAi41Gk1teiYzudnNzQ1hYGMLCwkridKVKUQZlnEzhVohSX3NZ847X101jAJWeQUemVJSnURl6upOhz2bsoCaxGurpu8LBWykZL2vkGkFJ6vQpXdTzxpN648HbeZUbYMen2X8deFEoTKPe7F3XU/8TntQ/h65yVR+BrplvU0/hMtSCILW8He2AnDyqARYXjSa3LJqCZUFFGZRxMgUYsJX7t/qay5p3vObuS9J1oTSm1svn+cC/AO7pPrahz7b/OvBOU9X+vm6Gn2ikqYXPy4t6PvebP/c7zbgbHH0jp42pWYotzSnlxgNQ9Z/PjANktmppNAeevcQvnCK1yVLfd0Z9BLp6vosa/Io6IMmY/XYMoD5pUvpRkLagogTS7Ve1X/O1ac07XnP1Jem7UBpb6wW4ZSo/8eH6EFf25ppkdS1Ywv/mAyoALDnONU/zI5qH7uAClJ1MNXpbCn6hPP63vr5XsWZrKX9DzdHtfJ+0lO9Cwj3Voi/qaYKrip9L/WZDSpOl3mCucSPQsx7w8atFD35FvYnUtx//92r68pH2vm7i3SfmQlOViDlQkLagogzKiPAXrrusueay+oVCylOhivK+vgul5sCtC3oWJ+GPpd7MvP2qdi2Z71/myW2ANxoDGy5yr/OZqml3/zXVQKl8xi1qwgdvfZ9VfSR0boH2Qin6mtR5Uudcq99MqQ8iM/RdaKoA/r6hOpevG7dEaupz7bQAMKmlcYPQ9AVzzRHoxQnQQNFaOwDd/2fUb5o85NxNX0mypqlKdLNQtpgtSP/999/w8/ODh4eHuU5R6hVlUEaoL9fErdknDYhfKHQNRjLU/6ivGVestqy5njYfWO1thOfVHF3MNzMD3O8If9W60k52AJh2MMwrBDYliq+JDc0JBDL9S3ryZW9ooRR984TtZMDbTYGBAarjSrlgiy6Rqud5xr5uwN8AFoYBvhVVLQ6OdsLR3QBXNkUZh6ArmPu5i49ALwp+VL/mXHWp+TO0FK1ms3xJsJapStZ0s0BMw2xBukmTJli1ahVGjhxprlOUCUUZlBHqK/5ABGMuFLrSSmnG1Wyy1WyeVj+25oAmvl9V/fPvGAD8/Sf3u5aX8CJ8L1M8//kMGPuKatSyeq3qq5OqvmWxFcf4fnJDDwERm6qlnu+oPqo1udee54K0MX8HzXQpGdJGePfxB3ZeEwalsc2An86qAt/6CN3Tpix90dac1635nTBE7P+M+k2To9muarpZy1Qla7lZIKZjtq+z2JOwiHkZc6HQlXb/NWGABlR9m+r4C+WWy+JTv/hj29tw3Zn6mjX5GqLmYibJ6dp98Dy5DTAokPvRDD6xIjUt9aVI+b5tzXzzDwHh8aO7tWrnL6VkiI+a1vd3MNcjPAcFAO1rSG9hMfbCbcoamjkCmvqNY9Mq3E2fMYp7E2MtU5Ws5WaBmA71SZcyhvo6pSxAITZiGgBWnuH6HTUVFmpv44ldFPha5i9/A7E3gBcFXHPs8BDjPqdmjV7d9FeFNfvkdGGTu+YocvVa7/Cd3GspFzO+yf+rk8A6jeZoXRdEsbLdcpm7ORi6Q9Wve+itol/Yxc7Bt2js/EcYRE1RuzJlDc1cAY3/u+flcTd9UpnqBsTcU5Wk3EhYy80CMZ0yG6SzsrIwa9Ys/Prrr0hLS4O/vz+mTZuGN9980+C+Dx8+xNSpU7Fnzx48f/4cwcHBWLBgATp27FgCOddNyvrSYhcKfY9/1NXEre7FyxWs3mlmONj7uXPTxPiAyMtnwA9ngKjzwovgyRRuWdCWGufUHJjV0gc4dVf1etFxIL+QG9l84T/D06A0a70pGVyg5mueYhczzQdfDNoGFDBVgNV3QVSfurT0BBeY1fuNcwu4RUbqehav9sbvJ9aioWsMQVFqV6auoZk7oAHSlwUtDU3ExtxIlETZkpJTZoN0v379EB8fj8WLF6N+/fqIjo7GwIEDUVhYiEGDBuncLzc3Fx07dsTTp0+xbNkyVK5cGcuXL0fXrl1x6NAhtGvXrgQ/BYcPjKnPNUY+GxiBzNO8CC2PB8Y1F69p6aRnAJbmaGXNAK1O/SLIz/l2lHHLgsbfA1rX4NJp9jG+UhU4e184UOyz4/qPr05ssJuumqeSRo9NgVqA5W9aNC+I6qO1+ePz8jWOp96fXtwBPvqCqClqV4aOYU193jypy4KWhibi0nAjQcyjTAbpffv2ITY2VhmYASAsLAy3b9/GlClTMGDAANja2oruu3r1aiQmJuKvv/5Cq1atlPsGBwdj6tSpOHXqVIl9DkAYGB3tVGsxSxmBzNNcQGPjJWFgUn/fwRZ4vSGwIVG1v72N7gFYmufbf108QItNt/lFo01yz79ckOYv+FF9uClcS08APyRw+VAf0S1G10VWM8iILYji5SScu6w+7cgWQIH6AUX6qTWfRqUvn7Yy6auA8XnW9zxpQ0HUFLUrXceQ0sJjiDmCvNRlQUtDE3FpuJEg5lEmg/T27dvh6uqKyMhIwfYRI0Zg0KBBOHXqFFq3bq1z3wYNGigDNADY2dlhyJAhmDFjBu7evQsfn5L7H6I5tWTGa1wwERuBLDYVClBdhJbHq+ZYq1+4xAKYepD+qJX2ACydFwqmvcnRDhgRDFR0Uj3T+GTKy7nUaloqtC/4H4aq5i+/KOQW0fBzB1afEz7hSW4DTGmt/5nJmkFGvabON0nz+Cdf/dyXC9q+bsCwncJR4zxlS0e2cNSy+vOyC9TKxUbjNZ9/TZplsd9AT42l+kTFbnjUv5eGWgmM7RM2x7Kg1t5EXBpuJIh5lMkgnZiYiIYNG8LOTvjxgoKClO/rCtKJiYlo06aN1nZ+30uXLukM0rm5ucjNzVW+zsjg5pbk5eUhT9dzBg1oWoVbnIEPDp1rCuca73+TW4lK4Qa8s0uVbscAYbpqzsC7TbglOPk0Tatwg2xSMrhjNFVw6TTPGVgJ2HKRe58/H59W82MFVlI9MAPgApItA9aeVeX/Ziowage3epZcBjjIuIPk5efhdDK3lrejjPstK1DlBQAOXQcc7IDRQUDmC+DXyy9PxIBAb/E8ianmrPosadlcjVk93ygEvjoO2NkAq3pxZXJA5LOnZHDNqjn5gL0t4C4Hcl+W2w89gHsZXACeelgjAxo18cNJ3PGaKlR/N82ySEjJgxwo8nepONQ/p+b3S/P7IisQ5vt0MlDNX/exNT+nvvQpGUCfTVwZO9gBOwdoz73ny2dr/zycf6j7u1raVHMG+ryc+26Kz8KXkyW+T6WJpctHxsw0V8rGxgY//fSTReZJ169fH7Vr18aBAwcE2+/fvw+FQoHPPvsM06dPF93X3t4eI0eOxIoVKwTbT5w4gdatWwua0DXNmzcP8+fP19oeHR0NZ2c9bZWEEEKs0vPnzzFo0CCkp6fDzU3is15NqEzWpAHofXa1oedaF3Xf6dOnY9KkScrXGRkZqF69OsLCwuDlZbpFhMVqNYBqG0+sRq1p51XugQ28yaHAiCaq86w6C2y9onq/Sx3gzQCuhqhe61PPky72tsCuN4HYJOFUr1aKPETIYxEeHg65XI6UDC4NZEB4bVU/saHjr+0DNFd7kpZ6C4Gh5yjzaRVuQOIDYNlpbvQ4b2EYt4iI5jHj7wEjdgo/44sCVU167F7t2mf8PWDMbuHx1fHn0vwMVZzyEBurKqfi0FczBoR51PdZ9B1fStmrl7vmd0rM2nPC747695WXl2e6cirLqJykSU1Ntej5zRak4+Li0KBBA3MdXi8vLy/Rgk1LSwPAPfPaHPs6ODjAwcFBa7tcLjfpf4KEB9yoVYDrc0t4wC3EsXuwsN+Zf0/fQwZa+AHMRtUf+8VpgNlyD6zQHJ0MADuvcz+AsP9QPU+6TArl8tLVjjvPi5fnPH0fiPAD1ifKwWzlCK4CfBXPnfur09w5anlxn099SpPchpsv/W+a+AIe3TbqHsyk/m8/d8DWDoAt4OvBDV4LUqhGqTvZceV0/7n2Me89A3LU2qKUU7byuFXBNP9Otby44+9/C/j8OLDnmrCM+HPxX5daXqq/H9/qJvZ9Mnbgldh3iD9PcjowdJdqYFtOHvc5V/ZWTVlTTyt2XvV86yL2NzI0eK5rA+DL06q55l0bqMpKkynKqTww9fWprLF02ZgtSFtiqhIvMDAQMTExyM/PF/RLX7zIPZEhICBA7758OnVS9i0p+hbRGNdcuPa1lAEzk1sBn71coelFAfdvQ6OTAe1VtvhapBj1wVZ+7sBHoaopVPwgqqUnuYCnfm5+qpmXM3eOd5pyNxB8AP3ihGrN6+R03Qt4qA9mspMBtjaqUfJRfbSXNg315ZbWVJ9HbXB1NY3Pb2vDjUh/UcgFFM2/xaGbqn/b23JlElxVtWyp1CBSlMU49K2/fvqu8G/Pr9KmOWUNKN4iIFJmC4h9tqKuIU7rWpPSqEw2d0dERGDVqlXYunUrBgwYoNy+bt06KBQKtGypuXSGcN/33nsPp06dUqbLz8/H+vXr0bJlSygUCp37lhRDi2gYWs4SENYoutXVftKU+uhkRztu1THNdbjlNsLAo7kyWVhN4Hgyt5/myIdu9YAlf2nPHdY8t72taqESRzvuyU58GoDLMx+w+fe71RN5YAZTfb58BuQXqPbXfPIWHyQ1gxL/MBA+L/zyovxNQupz4dztDWr3epq9JOpP/wKAzrW5AK15s2AoiCSnc60nxs6hVf8Oaa6/rr4aG78OeEqGeBkZe15jl0WVsnyrVDTXmJRGZTJId+vWDeHh4Rg7diwyMjJQt25dxMTE4MCBA1i/fr1yjvSoUaOwbt06JCUloUYNbhWNkSNHYvny5YiMjMTixYtRuXJlfP/99/jnn39w6NAhS34sAV0ri/EXXX6Rk69Pal/sxWoUUX2AwdtUQdPBVjX9iH/QxdcngZNqq34tDRfWXDUD7v0sVWDPLdC+KNrIIDply8FWFdTzCwA+9ufkc0HQQWOKu7JZ9uX7X58Sf2DGFye0Wwec7LSfvCU2j5o/Dp+vQqZa4lO9OZ2fV60pRyMoaM5d33MN+C1J2txpnthKcfqm4mniv0OaLQQpGdo3gcnp4gHVmLm7Yt87Q9OKDAVyY5qvaa4xKY3KZJAGgG3btmHmzJmYM2eOclnQmJgYwbKgBQUFKCgoEDwMxMHBAYcPH8bUqVPxwQcf4Pnz5wgJCcH+/fst2oRviK6FNMQu9rqeCqUeZAsKAUUF1bOYxfqn1QNeCx/VQiu8q4+F6dUHBJ2+q10zlwF4rxng4aCqkYq1uGsGQc2meV01rvURqlq3vS1Xe30riPuMYsFC7KEcyjnbBeLnU+860Myj2Cpg6mMI1FsQpAQRzZXi3mzM3XBo1sb5tGLrue+/BjzJFbYQ8M+nVn8il67WG2Pm7halVqyv1cjY5mtLzjWmvnBSVCJLKBhn7dq1yMrKMkVeTMrV1RXLli3D/fv3kZubiwsXLmit2x0VFQXGGGrWrCnYXqVKFaxbtw6pqanIzs7GiRMn0KlTpxLMvfE0H//HL44hdrHXHD3r68alUV9QI58BH8Vyi47oWjZU/Th+7lzNO7iK7jyqP5KQD+rqGIA6Fbnmas314MJqaqcHuG3rI7hFXvj3HWy5hUWS04VpQ32BI0O5tDJwtdfhO1V92eoBw88d+LyjamWw4Tu5z+v08raWX/0NEJZxt7raj0q0herRkSdTgI8Pcb/5MQT8MZ3suHRLw43rV+b3Hddcu1l6/zUukE2O5X7zZZKcDnT6hbsZ+uEM1zLwXjOupeCz48K0PM0yEtvGP+hE/Tz8a838aj4dTH0/Q+cFdN9s6qPrWObE30xo/g0IkcLomvStW7cEQW306NFo3LgxWrRoodx25coVNGzY0CQZJNJoNuXx/aRid+6az+9NyeAC2NJw4MODqoFcp+5ya2vPaSM+kEz9OLpq2zzNZliAq3XeTgd2vKxJqj8SU2YDQTW6lQ/wSXvtEdGTW3F5V1TgXj/JBtZe4GqzYk39fu7cim25an3SugYsTY5VlUV2PrdEqWYTumbtyM+de1TmigQgJpELfnI7Ln/8WuUAV3ve1F93LV4MX94pGdzIaV01Q0FfvI6lY9VbBQCuZSAt23C56KNZsxUbkCeW36IO6NJc3z31OXcsfcunWgL1hZPiMCpI/+9//8Pnn3+OO3fuKLdproVSUFCAgIAAxMfH45VXXlFuT0lJQWxsLEaMGFHMLBMxxjTlidWkk9OBjw9zQUmzq3jhn+Lbr6epaqGate2e9YCgKlzNml9Wk79gO9pxNTa+P3dNb+DReW7uLd9HqjmHOLgq1y/+W5Jqm9wG+PsBN9f748PSmvoB7QFgugYsad6ULD3B1ZQ1m4E1H5MJAL9e4gI0oOqPjr8nPN5XJ4Evw6UtSZmczs1r/sSH+717sGo/sYGD6jcS6kt08ts0uyccdfTNG0P9WeS6BuSJ1WL3Xy9aEFMOkrzO/W348QiGlk8tadQXTopDcpBu3749nj17hmPHjhlMK7aI2d27dzF69GgK0mYkdf1hsZq0ejOp5l+Pr01qbv8hAYi6oBr5rH4heiuIOybfr60+OEl9VHN2/sulM6G6WUjNFj7WEeAeR6k58CuvkKtVq9espfbr8l9R/rdmn6HmwC5AfPCbWC1QbApTCx/u8/H9zwDXUtF5vbSao+ZjMw09MEL9PbGbNz93bioTPwuAX/O8qH22yencKHyeg620oJ+czgVYnq6bJl3EWkYS7nHfJ2tB626T4pAcpM+cOYNJkyahdu3a5swPMYHkdCA6EbiTzg2OyisUPt1J1529+rbPOwJ/pQD1PVU1aTHqI58/bAlAJlwIhZ8WFVxV2DSpXpNuqgD+vsHlj1/cwsEWQCEXqPlmW/XAp1mr51+rN/Xzg70A4YVx/zXVRT23QPuBEPz+UX20n1WtGUB0jQLna6m2MmD0ywalUF+uifurk6pnY0sdxb30hGpGnUMRApnY8f3cucdtSklriOaUMr4bwlBw0mx2n9TS+PNrfp/575MmSw7eKmq5qqPBZ+WT5CB94cIFjBo1Cp07d8bBgwfNmScikdh/2pMpwqlUYitaqfcN+rqpalN8UIIMaFKN++m8XjtAaz628ffb3GAoPqB+GCqsNX92XBX8LjwAwLigzd80VHMG/gZXA+L3yy3Q/cQvfqWxT/5Q5YGBq32PCOFq74oK4v2cYjU+9TnU6vOu+f261dN9cdR1w8PPiy5g3MCsqPOqRVK+DBfmzdC0Ij6Q8Q8Aeb+Z9V2kNctBfeEafXnV2q+e8efWrKny3yd1pX0hk9Kef1J0koN0nTp1cPToUXz77bda7z169MikmSKG6XqG75Dt4guE8LLzucFXfH8xP9cXUI1Wzi0AvvwLaFhJeyCYnQx4uynw/RnVtp3/Co8Ppt1UnJ3P3QDw87fVLzT8cpdNFdoXen3Nto0rC2ul+YzL19rzXK1e14ApzRpft7qqfIn1aesbDSzWlLnlsvYa4+p5MHZaER/I2Mu/U3gd8bwYw9S1sqI26ZqqKVj9ZkDsoUWlffBWac8/KTqjR3ePHz9ea1vv3r3h5uaG4OBgBAYGQiaTISUlBU2aNDH4MAsiLjn95fOWGVe7kDr9xNBSnoCqH1dzxLZgtG/hy1qvhnwGpOfoXjZUbsPlt1s94Trb+kYa83zduNq2+lKc6jRrZaG+3HrgfM2Xl53PnUushitW41MOQLoGPM1VPata6iAfzXz5ummXj+axdNUwdc0lPjiEe4Qjbhh+WIgh5qqVFbVJ1xRNwYaU9sFbpT3/pOiMDtKa9u7di7Nnz+Ls2bNISEjA77//DoBbXtPR0REBAQEIDg6Gi4tLsTNbXiSnA+HrVbWxpSe5QT7qFzL1/7T8vODgKqptdjKgTQ2gTXVg0XEuYNjKhE3XeYXCAVrqNWl1Db2BpDTV4iMbErnfdjJu1TB+u51MNR8Y4NbZ7lbX8EhjXkqG9lKchvpqh+/U/hwAVxa6BkzpqrnxtWlHO2DGq+I3R4ao50luw01ryyuUXkvUty57NX9gn0Zfq2aNWEoN2ZS1stLST1raB2+V9vyToit2kO7WrRu6deumfJ2WlqYM2HzwXr16NRhjVKuWSLNJVmxUsXrt76tT3LxgXfOjw+uo+p/Vm7cd7YAlHYGDN7i0AwO47XOOAnG3VOdqX4Or2WqupJXPgLFNgYqO0FnjlzLSmKfeJ52dz63GNa65tGCj2cTPz/3WNWBKX8tETj73QI+iDqBSX1Qmr1A4ZcsQYy7GUuYliwVuUy21Wdr6SUuixm5OpT3/pGiKHaQ1eXp6olOnToIVujIyMpCQkIBz586Z+nRlkvo8XkD8CUrAy+knzqp02flccNIMCur/udWn3VR14Rbs4AdJDQzg0n3SHujws6o2OCiQ21fzIRwAUNGBqzFLpe9Co94nDXDTldRr1PqCjeaIcb681Ncz17W4C2C65kRTHEfqxVizRqzrQSFigdQUS22as5+0tNTQCTE3kwdpMW5ubggLC0NYWFhJnK7U83PnVq3S1yfNK0pQ8HLmgpZ6X252Pnc+LyfuPZuXjR78b/7CHpMIrDzDje52sC3aaFxdfN2017M2JtgA2k2/Yg+gEAs8phzAJPZwCnMEHM2/vdQHhfA3AVL7xE353ZOitNXQCTGnYgfpgIAAfPrpp4iIiJCU/v79+1i0aBEUCgWmTZtW3NOXWX7u0mqo6kFBc14w/wAFfrEKQPwhHADXX80P8lJ/j59HzAfvtee5AC234dbqNsXFU3O5S7FnYksNNur/FltzXF/g0RW4jA2y6scxZ8ARuyEwtEyooUBqTOA1Vz8pjWQmRKXYQfqNN97A0KFDUbFiRQwePBjt27fHK6+8Am9vb8hkMmRnZyMpKQknT57Ezp078dtvv6F58+Z49913TZF/AtUFTLN/Ur3/+auT3Eho9f5S9fnOBQwoKFC9p/4sabHgnVeovXJZUQYxiS13CagWRlGfhmVsrU1s1TBja3zFDbKaAWf/deO6BwwRu0kRG7sgNZAWJb2pA6jYjQI1f5PyqthBes6cORgzZgy++eYb/PTTT/j8888hk8kgk8kgl8vx4sULANxSoW3atMHGjRvRr1+/YmecCIn1T6qP0s7Jh2D+spMdN7BKfYAYP/rb3gZ4zQ9o4MUNCuMfG6kevNWfNKXZtKxvEBMgvOCqD5LLydde/Ut9UQx9/ahi2zVbGfT1SUstV2NrdZrjC/j1v0sy0BgbSC09QEmsG4P/bjnYqua2U7Am5YFJ+qSrVauGzz//HAsWLMCpU6dw4sQJ3Lt3D9nZ2fD29oa/vz/at28PX19fU5yOiNCsfbT2Fa4Tzfcfq6+edS9TGKS/7gxcecytyX3kFvczNFB4nqXhwH9ZqhHl/JOmpAxi0hXM+Uc7OtpBa/Uv9aAoFjzUjyd2AS9uwDFFv2unWqqV38RG6hNt6n839XXfcwt0P+GMkLLIpAPH5HI5XnvtNbz22mumPCzRQ70WqV774PumeZNDhYGL/72pv2rxkFBf4K9Dwv1+vih8nVeoPaKcH+DG1xh1DWICtIN5SgbwQw/uKVg/9AB8PVRrZas/bEFXbVn9ePwF/KuT3MA7U1zAi9PvKjZwTf2RipYOMKZoQi6JZmixbgvqqyblRYmM7iamxw8KU3/4w8EhwulXUtZEDvXlfngR/sIaOKBqBhd7GId6n7WDLbfetr4nKmnWTH3dgHd2cX3SY/cCK3uLP6FKV7+w2AWcbzY3Vd9vUWvjmgPXetYDYm+oHqloyZqgKQa0ldQobP67pPl9p1W3SHlAQboUEquhiTUNF6UGGOoLzGkjfHjF1521V83ij536XNVnnVvAjQLX19SsmS/NPmn1vnS+aZj/fPo+p+YCLH8/MH9t1VAtUvOGJKiyqtnb0jVBU4ygLslR2H7u3BO71LtrqBZNygMK0qWQ2NQisZpFUWuAo17hHl6x/SrXty22rCV/7JMp+p/fLBbINPOl3ietq5ncUL/wX3eEr/dcAw7fNF/tTkotUmwAlPoDRopTEyxuM7Mp+totsZ60pQe1EVLSKEiXQprrdmsOltJ3AdecO63rghfqyz3ukV9D3NFOu59Xc53qqD7a70sJZDsGAH//yf2u5SXeAqCvVeD0Xe31xgHz1u6k1iI1g4op5hWrP3e7qM3MppjjbK550oQQFbME6eTkZCgUCtjZ2Qn+TUzD0HQkXYExOR3o+AvwQm3utL4BVvuva0+PUu/n1VynWnPetNRA5uvGPf+Xf7qTrmZysX2T07mpYA62qn5xmcy0/ZZiNz1FXf/aFDVBzTXOi3ojYoq8UM2WEPMyeeQsKChArVq1cObMGQQFBaFWrVqIj4/HK6+8YupTlWu6Lo76AuP+66oADXCBTO8Fnul/bShQGdscuvMq0MJP+kVf/YZE/clVgOlqd7pueop6o2QKms/dpgFUhJRdZqneMsbAXg7N5X+TkqEvMD7JFqaV2+i/wHerxz0mk6+hao4QN9TcKbU5lK+Bz4wDCm2AzrWBt4KEo87F6HtylamCor6bnqLcKJkCv8Y5NTMTUvZRG3QZoyswJqcDay+o0tlC+OxnXcc69Jb+pT4NNXdKaQ5NuAfIX/77RQE36GvPNW4Ot75AXRIDl4pyjpLIV2lpZqblPAkpHgrSZRB/MVR/2Ib6VCcA+Pg1wzVVfl/1QG+OZtymCuDvG9rbt1/Vn0dzDVzSDCy6HmBS0vkqbehpVoQUHwXpUsCY2oiuRU40a3f8mtjGMFczLj9wbFQTYPlZ1fYIf8P7mrpGqS+wGBNwSktN15zoaVaEFB8FaSt3MkX13GdDwUHfIievN+KmSPFLgFpqbq2+G44PQ4G2tYTLlJY0XYGlpAJOWWoetsQ8akLKGgrSViw5XRWgAcPBQd8iJ/yc5ux8brEQS8ytldL8qblMaUnTFVikBpziBNmy1jxclpr9y9LNEyldKEhbsdN3VQEaMDwaW98iJ+pPErLU3NrS0PypK7BICTimfva0NZaPscpCs39Zu3kipQsFaSumHnTlNtJGY+sKJNbQ9GjOPJiypqMrsBgKOKZ49rSl/0ZEW1m8eSKlBwVpK2ZMc6F6kFJ/ElZRjmUu5hyNbQ01neIGWWv4GxFtdPNELMnkQdrW1hZxcXFo0KCB4N+kaKQ0F0oNUqZqeixOrdUczZ/WUtMxRZAtC83DZQ3dPBFLMktNul27dqL/JuZhbJAy1eAmvgnekgO9AOuq6VCQLZvo70oshZq7ywBjgpQpBzflFXKjz48MtewFjGo6hJCyioJ0GSAWpHTVlk0xuIl/fjTA/baGgTRU07EMmppEiHnZmOpAL168wNtvv42LFy+a6pDECH7u3IAxPkB3Xg9MjuV+J6er0vG1bqDog5vWR3CBuqjHIGWDvu8ZIcQ0TFaTtre3x/r16zFkyBBTHZIUkaEnNxW3aTjUl2viphpU+WYtA/YIKctMVpMGgJCQENy8edOUhyR6JKdzi5Ro1mAM1ZbVa91FZYpjkNKtuK0yYnR9pwkpr0zaJz1r1ixMnjwZkZGRcHZ2NuWhiQZ9A8BoIBUpCab+nlnLfHdCrIlJa9KLFy9GVlYWXnnlFezYsQO5ubmmPDxRI9bUqI5qukVHtTnpTPk9M/SdJqQ8MmlN2tHREbm5ubh79y769+8PBwcHBAUFoVmzZsqfgIAAU56y3LKmucFlCdXmLIe+04RoM2mQPnjwIAAgOTkZCQkJOHv2LBISErBlyxZ8//33kMlkKCgoMOUpyy1q0jYPGgxlOfSdJkSbWeZJ+/n5wc/PDxEREcptd+/eRUJCgjlOV27R3GDTo9qcZdF3mhChElvMxMfHBz4+dMUj1o1qc4QQa0IrjhGiwZS1OVqRixBSHBSkCTETGoRGCCkuk07BIoSo0JQiQkhxUZAmxEzMsSIXIaR8oeZuQsyEBqERQoqLgjQhZkRTigghxUHN3YQQQoiVoiBNCCGEWCkK0oQQQoiVoiBtBeipS4QQQsTQwDELowUvCCGE6EI1aQujBS8IIYToQkHawmjBC0IIIbpQc7eF0YIXhBBCdKEgbQVowQtCCCFiqLmbEEIIsVIUpAkhhBArVWaDdFZWFiZOnAiFQgFHR0eEhIRg48aNRTrWrFmzIJPJEBAQYOJcEkIIIbqV2T7pfv36IT4+HosXL0b9+vURHR2NgQMHorCwEIMGDZJ8nPPnz+PLL79ElSpVzJhbQgghRFuZDNL79u1DbGysMjADQFhYGG7fvo0pU6ZgwIABsLW1NXic/Px8jBgxAu+88w4uXLiAx48fmzvrhBBCiFKZbO7evn07XF1dERkZKdg+YsQI3Lt3D6dOnZJ0nMWLFyMtLQ0LFy40RzYJIYQQvcpkTToxMRENGzaEnZ3w4wUFBSnfb926td5jXL58GQsWLMC2bdvg6uoq6by5ubnIzc1Vvs7IyAAA5OXlIS8vz5iPUK7wZUNlpB+VkzRUTtJQOUlj6fIpk0E6NTUVtWvX1tru6empfF+fwsJCjBw5Ev369UP37t0ln3fRokWYP3++1va4uDg4OztLPk55FRsba+kslApUTtJQOUlD5aTf8+fPLXp+qw/SR48eRVhYmKS0586dQ0hICABAJpPpTKfvPQD46quvcO3aNezatUtyPgFg+vTpmDRpkvJ1RkYGqlevjrCwMHh5eWmlT8kAEu4BTRWAr5tRpypT8vLyEBsbi/DwcMjlcktnx2pROUlD5SQNlZM0hip15mb1QbpBgwZYtWqVpLR+fn4AAC8vL9GCTUtLA6CqUYtJTk7GnDlzsHjxYtjb2+Pp06cAuEFkhYWFePr0KRwcHODk5KS1r4ODAxwcHLS2y+Vyrf8EyelAt4309Ct1YuVEtFE5SUPlJA2Vk36WLhurD9LVqlXD6NGjjdonMDAQMTExyM/PF/RLX7x4EQD0zne+ceMGsrOzMWHCBEyYMEHr/YoVK2LChAn45ptvjMqTJrGnX5X3IE0IIUTI6oN0UURERGDVqlXYunUrBgwYoNy+bt06KBQKtGzZUue+ISEhiIuL09o+ceJEpKenY+3atfD19S12HvmnX/E1aXr6FSGEEE1lMkh369YN4eHhGDt2LDIyMlC3bl3ExMTgwIEDWL9+vWCO9KhRo7Bu3TokJSWhRo0a8PDwQPv27bWO6eHhgfz8fNH3ioKefkUIIcSQMhmkAWDbtm2YOXMm5syZg7S0NPj7+yMmJgZvvvmmIF1BQQEKCgrAGCvxPNLTrwghhOhTZoO0q6srli1bhmXLlulNFxUVhaioKIPHO3r0qGkyRgghhEhUJlccI4QQQsoCCtKEEEKIlaIgTQghhFgpCtKEEEKIlaIgTQghhFgpCtKEEEKIlaIgTQghhFgpCtKEEEKIlaIgTQghhFgpCtKEkFJl9+7d8Pf3R1ZWlsG0n3/+OTp37lwCuSLEPChIE0JKlU2bNqFq1apwdXU1mLZFixaIjY3FjRs3SiBn0mVlZWHixIlQKBRwdHRESEgINm7caPL9paY7evQoZDKZ6M/JkyeL9VlJ8ZTZtbsJIWVPfn4+9u3bhzlz5khK36ZNG3h6emLHjh2YNGmSmXMnXb9+/RAfH4/Fixejfv36iI6OxsCBA1FYWIhBgwaZbH996SIjI7WO+9lnnyEsLEywLSAgoPgfmBQdI2aTnp7OALDHjx9bOitW7cWLF2zHjh3sxYsXls6KVaNyYuzw4cMMAEtKStKZhi+nzMxMxhhjb731Fmvbtm1JZdGgvXv3MgAsOjpasD08PJwpFAqWn59vkv0NpcvOzlZ+n+Li4hgAtnnzZhN8wrLl8ePHDABLT0+3yPmpuZsQUmrs3LkTgYGBqF27tmB7586d0aZNG+zduxdt2rTBG2+8gfHjxwMA+vTpg+PHjyM1NdUSWdayfft2uLq6atVkR4wYgXv37uHUqVMm2d9QutOnT5vg0xBzoyBNCCk1du3ahT59+mhtP3fuHG7cuIFJkyZh9OjRmD17Nt577z0AQJcuXWBnZ4c9e/YYPD5jDPn5+ZJ+iioxMRENGzaEnZ2wtzEoKEj5vin2N5Tu0qVLWsceN24c7Ozs4Obmhi5duuDPP/804pMRc6AgTQgpFS5cuIBbt25pBenbt2/j8ePHcHZ2xpkzZzBs2DAEBgYiODgYAPds+Y4dO2Lnzp0Gz3Hs2DHI5XJJP7du3SrS50hNTYWnp6fWdn6boRq/1P2NOY+7uzsmTJiAlStXIi4uDsuWLcOdO3fQvn17/PbbbxI/GTEHGjhGCCkV4uLi4OrqiqZNmwq2JyQkAAAWL16MChUqIC8vT2vfsLAwLFy40OA5mjZtivj4eEn5USgUet/XrG3b2tpCJpMBgPK3GH3vSUmj/p7UdE2aNEGTJk2Ur9u0aYOIiAgEBgZi6tSp6NKli8E8EfOgIE0IKRWaN2+OrKws/PPPP/D391duP3v2LBwcHNC9e3ed+8bHx6N58+YGz+Hq6oqQkBBJ+dFsRlZ369Yt1KpVS7AtLi4O7du3h5eXl2htOS0tDQBEa7/qpO5f3PN4eHigZ8+eWLFiBbKzs+Hk5KQ3PTEPau4mhJQKrVq1QuXKlbWarRMSEhAUFKQziLx48QIHDhwQ7cvWZKrmboVCgfj4eMEP3wIQGBiIK1euaNW0L168CMDwlCep+xtK17hxYwOlwfXRA9Jq98Q8KEgTQkoFGxsb9OzZEzt27BBsP3v2LJo1a6ZzvyNHjiAjIwO9e/c2eA6+uVvKj77mbnt7ezRr1kzwU6FCBQBAREQEsrKysHXrVsE+69atg0KhQMuWLfXmUer+htK1aNFC73mePHmCPXv2ICQkBI6OjnrTEvOh5m5CSKnRt29frF27Fv/99x+qVq2KlJQUPHz4UG+Q3rlzJ5o2bYrq1asbPH6FChX0HssUunXrhvDwcIwdOxYZGRmoW7cuYmJicODAAaxfvx62trbKtMeOHUPHjh0xZ84c5QIuUvc35jyDBg2Cn58fmjVrBm9vb1y7dg1Lly7FgwcPEBUVZdbyIPpRkCaElBqdOnWCk5MTdu/ejTFjxigHjekKrIwx7N69G++8805JZtOgbdu2YebMmZgzZw7S0tLg7++PmJgYvPnmm4J0jDEUFBSgsLCwSPvrS6c+wC4oKAibNm3CihUrkJWVBU9PT7z22mv45ZdfJPXlE/ORMb7TgZhcRkYG3N3d8fjxY3h5eVk6O1YrLy8P+/btQ/fu3SGXyy2dHatF5cSJiIhAXl6eznnP6uV0/vx5tGjRAhcuXFDODyYc+j5Jk5qaCm9vb6Snp8PNza3Ez0990oSQUqVPnz44fPiwpKdg7dy5EzVr1qQATUotCtKEkFJl+PDhyM7OlvQUrAULFuDmzZslkCtCzIOCNCGEEGKlKEgTQgghVoqCNCGEEGKlKEgTQgghVoqCNCGEEGKlKEgTQgghVoqCNCGEEGKlKEgTQgghVoqCNCGEEGKlKEgTQgghVoqCNCGEEGKlKEgTQgghVoqCNCGEEGKlKEgTQgghVoqCNCGEmEhWVhYmTpwIhUIBR0dHhISEYOPGjZL2zczMxNSpU9G5c2dUqlQJMpkM8+bN00p35MgRjBw5Ev7+/nBxcYGPjw/69OmDhIQErbRHjx6FTCYT/Tl16lRxPy4pARSkCSHERPr164d169Zh7ty52L9/P5o3b46BAwciOjra4L6pqan48ccfkZubi759++pM98MPP+DWrVuYMGEC9u3bh2XLluHhw4cIDQ3FkSNHRPf57LPPcOLECcFP48aNi/oxSQmys3QGCCGkLNi3bx9iY2MRHR2NgQMHAgDCwsJw+/ZtTJkyBQMGDICtra3O/WvUqIEnT55AJpPh8ePH+Omnn0TTLV++HJUrVxZs69q1K+rWrYvPPvsMHTp00NqnXr16CA0NFWzLy8sz9iMSC6CaNCGk1MjPz8cXX3yBoKAgODk5aTXh+vn5WSxv27dvh6urKyIjIwXbR4wYgXv37hlsXuY/gyGaARoAXF1d0ahRI9y5c8e4TBOrR0GaEFJqjBo1CtOnT0enTp2wfft2REVFKYPWhx9+iPnz5xt1PMYY8vPzJf0YkpiYiIYNG8LOTthAGRQUpHzfXNLT03H27FmdTdjjxo2DnZ0d3Nzc0KVLF/z5559mywsxLWruJoSUCtHR0fj555/x448/YsyYMcrtjDGMGDECHTp0QM+ePY1qxj127BjCwsIkpb158yZq1qyp8/3U1FTUrl1ba7unp6fyfXMZN24cnj17hpkzZwq2u7u7Y8KECWjfvj28vLxw/fp1fPHFF2jfvj127txptvwQ06EgTQgpFVauXImGDRsKAjQANGzYEACQlpYmul9+fr5W7ZbXtGlTxMfHSzq/QqEwmEZfc7WUpuyimD17NjZs2IDvvvsOTZs2FbzXpEkTNGnSRPm6TZs2iIiIQGBgIKZPn45PP/3ULHkipkNBmhBi9Z48eYI//vgDkyZN0novJSUFAODr6wsA+O6777B582bExMRg7969iImJQc+ePUWP6+rqipCQEEl50BXoeV5eXqK1Zf7mga9Rm9L8+fOxYMECLFy4EO+//76kfTw8PNCzZ0+sWLECubm5Js8TMS3qkyaEWL07d+6AMSZam92xYwe8vb3Rtm1bAMDFixdx6dIljB49Gunp6ejSpYvO4x47dgxyuVzSz61bt/TmMTAwEFeuXNHqv7548SIAICAgwMhPrd/8+fMxb948zJs3DzNmzDBqX8YYAPPV7onpUE2aEGL1PDw8AABXrlwRbP/rr78QExODuXPnKmu6Fy9exIABA9CuXTvY2NjAxkZ3XcSUzd0RERFYtWoVtm7digEDBii3r1u3DgqFAi1btpR0Hik+/fRTzJs3D7NmzcLcuXON2vfJkyfYs2cPgoODYW9vb7I8EfOgIE0IsXp+fn5o27YtoqKiUKtWLbRo0QKnT5/GZ599hs6dOysHTBUWFuLKlSt47733JB23QoUKaNasmUny2K1bN4SHh2Ps2LHIyMhA3bp1ERMTgwMHDmD9+vXKOdLHjh1Dx44dMWfOHMyZM0dwjP379+PZs2fIzMwEAFy+fBlbtmwBAHTv3h3Ozs5YunQp5syZg65du6JHjx44efKk4Bjq86EHDRoEPz8/NGvWDN7e3rh27RqWLl2KBw8e4KeffqLm7tKAEbNJT09nANjjx48tnRWr9uLFC7Zjxw724sULS2fFqpX3cnrw4AEbOHAg8/DwYPb29iwgIIB9+eWXLC8vT5nmn3/+YVWqVLFYOWVmZrLx48ezqlWrMnt7exYUFMRiYmIEaeLi4hgANnfuXK39a9SowQCI/ty8eZMxxli7du10ptG8pC9atIiFhIQwd3d3ZmtryypVqsQiIiLY6dOny/33SarHjx8zACw9Pd0i56eaNCGkVKhcubLB5TX//vtv5bxkS3B1dcWyZcuwbNkynWnat2+v7BPWZKjfG+DW45Zq2rRpmDZtmuh7tOJY6UADxwghZcbFixcRGBho6WwQYjJUkyaElBnz589HXl4e9u3bZ+msEGISVJMmhBBCrBQFaUIIIcRKUZAmhBBCrBQFaUIIIcRKUZAmhBBCrBQFaUIIIcRKUZAmhBBCrFSZDdJZWVmYOHEiFAoFHB0dERISgo0bN0rePy4uDuHh4ahcuTJcXV0RFBSEb7/9FgUFBWbMNSGEEKJSZhcz6devH+Lj47F48WLUr18f0dHRGDhwIAoLCzFo0CC9+x46dAhdunRB27ZtsWrVKri4uGDXrl2YMGECkpKS9C75RwghhJhKmQzS+/btQ2xsrDIwA0BYWBhu376NKVOmYMCAAcon0oiJioqCXC7Hnj174OLiAgDo1KkT/vnnH0RFRVGQJoQQUiLKZHP39u3b4erqisjISMH2ESNG4N69ezh16pTe/eVyOezt7eHk5CTY7uHhAUdHR5PnlxBCCBFTJmvSiYmJaNiwofIh8Dz+6TiJiYlo3bq1zv3fffddxMTEYPz48ZgxYwacnZ2xe/dubN++HYsWLdK5X25uruD5rBkZGQC4p83QE2d048uGykg/KidpqJykoXKSxtLlUyaDdGpqKmrXrq213dPTU/m+Pi1btsSRI0cQGRmJ5cuXAwBsbW2xaNEiTJ48Wed+ixYtwvz587W2x8XFwdnZ2ZiPUC7FxsZaOgulApWTNFRO0lA56ff8+XOLnt/qg/TRo0cRFhYmKe25c+cQEhICAJDJZDrT6XsPABISEhAREYGWLVti5cqVcHFxwZEjRzBr1izk5ORg9uzZovtNnz4dkyZNUr7OyMhA9erVERYWBi8vL0mfoTzKy8tDbGwswsPDIZfLLZ0dq0XlJA2VkzRUTtIYqtSZm9UH6QYNGmDVqlWS0vr5+QEAvLy8RAs2LS0NgKpGrcu4ceNQpUoVbN++XTnALCwsDDY2Npg3bx4GDx4sWlN3cHCAg4OD1na5XE7/CSSgcpKmPJZTVFQURowYoXxta2uLSpUqoV27dvj0009Rr149QfqZM2eib9++qFWrFpKSkrRuzH///Xe0a9cOALB27VoMHz5c+d6pU6ewePFiJCQk4MGDB/Dw8EDt2rXRunVrLF261Hwf0kIMfZ9u3LiBjz76CEeOHEF+fj5atWqFzz//HK+88orBY//555+IiorCuXPnkJiYiBcvXuDmzZuoWbOmIN2///6LH3/8EXFxcUhKSoKNjQ0aNmyIDz/8EK+//rogreZ3Qd39+/dRtWpVwx/aCJb+v2b1QbpatWoYPXq0UfsEBgYiJiYG+fn5gn7pixcvAgACAgL07n/+/HkMHDhQawR48+bNUVhYiCtXrogGaUKIea1duxb+/v7IycnB8ePHsXDhQsTFxeHq1auoWLGiIG2FChVw8+ZNHDlyBB07dhS8t2bNGri5uSnHjfD27t2L3r17o3379liyZAmqVauG+/fv48yZM9i4cWOZDNL6PHr0CG3atEHFihWxZs0aODo6YtGiRWjfvj3i4+PRoEEDvfsfPnwYhw4dQpMmTeDm5oajR4+Kpjt48CD27t2Lt956C82bN0d+fj42bdqEyMhIzJ8/H3PmzNHah/8uqCuTLZasDNq3bx8DwDZu3CjY3rVrV6ZQKFh+fr7e/WvVqsUCAgK00s2YMYMBYOfPn5eUj/T0dAaAPX782LgPUM68ePGC7dixg7148cLSWbFq5bmc1q5dywCw+Ph4wfb58+czAGzNmjXKbS9evGCNGzdmjRo1YqGhoWzQoEGCfTIyMpizszMbM2YMA8DWrl2rfK9t27asTp06LC8vTysPBQUFpv1QGp4/f84KCwtF33v27JnJzyfl+zRlyhQml8vZrVu3lNvS09OZt7c3e+ONNwyeQ73MvvjiCwaA3bx5Uyvdo0ePRD97jx49mLOzM8vJyVFu0/VdMJfHjx8zACw9Pb1EzqepTE7B6tatG8LDwzF27FisWrUKcXFxePvtt3HgwAEsWbJEUEMeNWoU7OzscPv2beW2Dz/8EImJiejVqxd27tyJ2NhYTJs2DUuWLEGnTp0QHBxsiY9FCNHQrFkzAMCDBw9E3x85ciS2bduGp0+fKrfxKw+++eabWulTU1Ph7e2tNTMEAGxspF0uz5w5g969e8PT0xOOjo5o0qQJfv31V0GaqKgoyGQyHDx4ECNHjkSlSpXg7OyM3NxctG/fHgEBAfj999/RunVrODs7Y+TIkZLObWrbt29Hhw4dUKNGDeU2Nzc39OvXD7t370Z+fr7e/aWWmbe3t+hYoRYtWuD58+fKrsryqEwGaQDYtm0b3nrrLcyZMwddu3bFqVOnEBMTg8GDBwvSFRQUoKCgAIwx5bYPPvgAW7duRWZmJkaPHo2IiAjs2bMHc+fOxY4dO0r4kxBCdLl58yYAoH79+qLvv/nmm7C1tUVMTIxy2+rVq/H666/Dzc1NK32rVq1w6tQpjB8/HqdOnTJ6+k1cXBxeffVVPH36FCtWrMDOnTsREhKCAQMGICoqSiv9yJEjIZfL8csvv2DLli3K/s/79+9jyJAhGDRoEPbt24f33ntPch4KCwuRn59v8MfQEsfZ2dlISkpSTl1VFxQUhOzsbNy4cUNyvooiLi4OlSpVQuXKlbXe69mzJ2xtbeHp6Yl+/fohMTHRrHmxGIvU38sJau6Wpjw34xqjPJcT38R58uRJlpeXxzIzM9mBAwdY1apVWdu2bQXN0+rN3YwxNmzYMNasWTPGGGOXLl1iANjRo0dZfHy8VnP348eP2WuvvcYAMABMLpez1q1bs0WLFrHMzEyD+fT392dNmjTRai7v2bMnq1atmrL5l/88Q4cO1TpGu3btGAB2+PBho8uJ/7x8/vX9tG3bVu/36e7duwwAW7RokdZ70dHRDAD766+/JOdLX3O3mFWrVjEAbNmyZYLt+/fvZzNnzmS7d+9mx44dY//73/+Yr68vc3FxkdwVaQxLN3db/cAxQgjhhYaGCl43bNgQO3fuFG2e5o0cORLt2rXDxYsXERUVhTp16qBt27ZISEjQSuvl5YU//vgDZ86cweHDh3HmzBkcPXoU06dPx8qVKxEfHw9vb2/R81y/fh1Xr17Fl19+CQCCpuDu3btjz549+Oeff9CwYUPl9v79+4seq2LFiujQoYPugtBj3rx5eP/99w2mc3R0RFJSksF0xZnOWlT79+/HuHHj8Prrr+ODDz4QvNe1a1d07dpV+bpt27bo0aMHAgMDMWfOHOzcudMsebIUCtKEkFLj559/RsOGDZGZmYlNmzZh5cqVGDhwIPbv369zn7Zt26JevXpYuXIlfv31V0ycONFgcGnWrJmyvzsvLw8ff/wxvv76ayxZsgRLliwR3YfvF//oo4/w0UcfiaZ5/Pix4HW1atVE0+naLoWfnx98fX0NpsvPz9cbpCtWrAiZTFas6axF8dtvv6Ffv34IDw/Hhg0bJN0I1KxZE6+99hpOnjxp8vxYGgVpQkip0bBhQ2XwDAsLQ0FBAX766Sds2bJFaz6tuhEjRmDWrFmQyWQYNmyYUeeUy+WYO3cuvv76a739nnwNe/r06ejXr59oGs0pS7oCUHFqqCNHjsS6desMpmvbtq1g8SVNTk5OqFu3rnLqqrqLFy/CycnJ5FNRf/vtN/Tt2xft2rXD1q1bYW9vL3lfxpjkgWqlCQVpQkiptWTJEmzduhVz5sxBv379dF6khw0bhlOnTqFhw4bw8fHRebz79++L1mKvXLkCAFAoFDr3bdCgAerVq4cLFy7gs88+M/KTmI4pm7sjIiLwzTff4M6dO6hevToAIDMzE9u2bUPv3r31djMY6+DBg+jbty9ee+017NixQ3RhKF1u3ryJ48ePo1OnTibLj7WgIE0IKbUqVqyI6dOnY+rUqYiOjsaQIUNE0ykUCkkzM7p06QJfX1/06tUL/v7+KCwsxPnz57F06VK4urpiwoQJevdfuXIlunXrhi5dumD48OHw8fFBWloarly5grNnz2Lz5s1F+Zi4ffs26tSpg2HDhmH16tV609asWVNrRS8xeXl5BoP0Rx99hF9++QU9evTAJ598AgcHByxevBg5OTmYN2+eIG3dunUBcH3zvEePHuHYsWMAVItJ7d+/H5UqVVKuGAdwK5P17dsXVatWxYwZM3D+/HnBsRs1aqQcjd+pUye0bdsWQUFBcHNzw8WLF7FkyRLIZDJ8+umnBj93aUNBmhBSqn3wwQf43//+h08++UT5/PiimjVrFnbu3Imvv/4a9+/fR25uLqpVq4ZOnTph+vTpgkFfYsLCwnD69GksXLgQEydOxJMnT+Dl5YVGjRrhjTfeKHK+GGPK6aIlqVKlSvjjjz/w0UcfYdiwYcplQY8ePaq12pfYnOlLly5pPTKYn07Wrl075Qpkhw4dQnZ2Nm7duiU6YC4uLg7t27cHwK0ouWnTJnz55ZfIzs5G5cqV0aFDB8yePVvnVLzSTMaY2gRhYlIZGRlwd3fH48ePy+ZydSaSl5eHffv2oXv37hZfJ9eaUTlJQ+UkDZWTNPwCN+np6aJz682t7PWyE0IIIWUEBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqwUBWlCCCHESlGQJoQQQqyUnaUzUJYxxgAAmZmZkMvlFs6N9crLy8Pz58+RkZFB5aQHlZM0VE7SUDlJk5mZCUB1PS9pFKTNKDU1FQBQq1YtC+eEEEJIcaSmpsLd3b3Ez0tB2ow8PT0BAMnJyRb545YWGRkZqF69Ou7cuQM3NzdLZ8dqUTlJQ+UkDZWTNOnp6fDz81Nez0saBWkzsrHhuvzd3d3pP4EEbm5uVE4SUDlJQ+UkDZWTNPz1vMTPa5GzEkIIIcQgCtKEEEKIlaIgbUYODg6YO3cuHBwcLJ0Vq0blJA2VkzRUTtJQOUlj6XKSMUuNKyeEEEKIXlSTJoQQQqwUBWlCCCHESlGQJoQQQqwUBWkzyMrKwsSJE6FQKODo6IiQkBBs3LjR0tkqEUeOHMHIkSPh7+8PFxcX+Pj4oE+fPkhISNBKe/bsWXTq1Amurq7w8PBAv379cOPGDdHjfvfdd/D394eDgwNq1aqF+fPnIy8vz9wfp8T89NNPkMlkcHV11XqvvJfTn3/+ie7du6NixYpwcnJCvXr18OmnnwrSlPcyOnfuHPr27QuFQgFnZ2f4+/vjk08+wfPnzwXpylM5ZWZmYurUqejcuTMqVaoEmUyGefPmiaY1R7k8fPgQw4cPh7e3N5ydndGqVSscPnzY+A/CiMmFh4czDw8PtmLFCnbkyBE2evRoBoBt2LDB0lkzu9dff52FhYWx77//nh09epRt3ryZhYaGMjs7O3b48GFluitXrrAKFSqwNm3asL1797KtW7eyxo0bM4VCwR4+fCg45oIFC5hMJmPTp09ncXFxbMmSJcze3p6NGTOmpD+eWaSkpDB3d3emUCiYi4uL4L3yXk4bNmxgNjY27M0332S7du1iR44cYatWrWLz589XpinvZXTp0iXm6OjIgoOD2aZNm9jhw4fZ3Llzma2tLevdu7cyXXkrp5s3bzJ3d3fWtm1b5TV47ty5WunMUS45OTksICCA+fr6svXr17ODBw+yPn36MDs7O3b06FGjPgcFaRPbu3cvA8Cio6MF28PDw5lCoWD5+fkWylnJePDggda2zMxMVqVKFdaxY0fltsjISObt7c3S09OV227dusXkcjmbOnWqctvjx4+Zo6Mje/vttwXHXLhwIZPJZOzSpUtm+BQlq2fPnqxXr15s2LBhWkG6PJdTSkoKc3FxYWPHjtWbrjyXEWOMzZw5kwFg169fF2x/++23GQCWlpbGGCt/5VRYWMgKCwsZY4w9evRIZ5A2R7ksX76cAWB//fWXclteXh5r1KgRa9GihVGfg4K0iY0ePZq5urqyvLw8wfbo6GgGgB0/ftxCObOssLAwVr9+fcYY92V1cnJi77zzjla6zp07s3r16ilfr1+/ngFgJ06cEKS7d+8eA8AWLlxo3oyb2S+//MIqVKjA7ty5oxWky3s5zZs3jwFgt27d0pmmvJcRY6pyevTokWD71KlTmY2NDcvKyir35aQrSJurXDp16sQaNGigdczPPvuMAWApKSmS80590iaWmJiIhg0bws5OuCx6UFCQ8v3yJj09HWfPnkXjxo0BAElJScjOzlaWibqgoCBcv34dOTk5AFTlFRgYKEhXrVo1eHt7l+ryfPjwISZOnIjFixfD19dX6/3yXk6///47PD09cfXqVYSEhMDOzg6VK1fGu+++i4yMDABURgAwbNgweHh4YOzYsbhx4wYyMzOxZ88erFy5EuPGjYOLiwuVkw7mKpfExESdxwSAS5cuSc4jBWkTS01NFX1aCr+Nf3xleTJu3Dg8e/YMM2fOBKAqA13lxBjDkydPlGkdHBzg4uIimrY0l+d7772HBg0aYOzYsaLvl/dyunv3Lp4/f47IyEgMGDAAhw4dwpQpU/Dzzz+je/fuYIyV+zICgJo1a+LEiRNITExEnTp14Obmhl69emHYsGFYtmwZAPou6WKucjFlHKCnYJmBTCYr0ntl0ezZs7FhwwZ89913aNq0qeA9qeVUFstz69at2L17N86dO2fwM5TXciosLEROTg7mzp2LadOmAQDat28Pe3t7TJw4EYcPH4azszOA8ltGAHDr1i306tULVapUwZYtW1CpUiWcOnUKCxYsQFZWFlavXq1MW57LSR9zlIupypBq0ibm5eUlepeUlpYGQPyOrayaP38+FixYgIULF+L9999Xbvfy8gIgfjeZlpYGmUwGDw8PZdqcnBytqSR82tJYnllZWRg3bhw++OADKBQKPH36FE+fPsWLFy8AAE+fPsWzZ8/KfTnxn79Lly6C7d26dQPATZsp72UEANOmTUNGRgZ+++039O/fH23btsWUKVPwzTffYM2aNTh27BiVkw7mKhdTxgEK0iYWGBiIK1euID8/X7D94sWLAICAgABLZKvEzZ8/H/PmzcO8efMwY8YMwXt16tSBk5OTskzUXbx4EXXr1oWjoyMAVf+PZtr//vsPjx8/LpXl+fjxYzx48ABLly5FxYoVlT8xMTF49uwZKlasiMGDB5f7chLr0wMA9vJxAzY2NuW+jADg/PnzaNSokVYzbPPmzQFA2Qxe3stJjLnKJTAwUOcxASPjgOQhZkSSffv2MQBs48aNgu1du3YtF1OwGGPsk08+YQDYrFmzdKZ54403WOXKlVlGRoZy2+3bt5m9vT37+OOPldtSU1OZo6Mje/fddwX7L1q0qFRNB1GXnZ3N4uLitH66dOnCHB0dWVxcHLt48SJjrHyX02+//SY6mvirr75iANgff/zBGCvfZcQYN3OiUqVKLDMzU7D9xx9/ZADYjh07GGPlu5z0TcEyR7l8//33DAA7efKkclteXh5r3Lgxa9mypVF5pyBtBuHh4axixYrsxx9/ZEeOHGFjxoxhANj69estnTWz+/LLLxkA1rVrV3bixAmtH96VK1eYq6sra9u2Ldu3bx/btm0bCwgI0LuAwIwZM9jRo0fZF198wRwcHErVwgpSiM2TLu/l1KtXL+bg4MA+/fRTFhsbyxYtWsQcHR1Zz549lWnKexnt3LmTyWQyFhoaqlzMZOHChczV1ZU1atSI5ebmMsbKZznt27ePbd68ma1Zs4YBYJGRkWzz5s1s8+bN7NmzZ4wx85RLTk4Oa9y4MatevTrbsGEDi42NZREREbSYibXIzMxk48ePZ1WrVmX29vYsKCiIxcTEWDpbJaJdu3YMgM4fdWfOnGEdO3Zkzs7OzM3NjfXt21drQQbesmXLWP369Zm9vT3z8/Njc+fOZS9evCiJj1RixII0Y+W7nJ4/f84+/vhjVr16dWZnZ8f8/PzY9OnTWU5OjiBdeS4jxhg7cuQI69y5M6tatSpzcnJi9evXZ5MnT2aPHz8WpCtv5VSjRg2d16KbN28q05mjXP777z82dOhQ5unpyRwdHVloaCiLjY01+jPQ86QJIYQQK0UDxwghhBArRUGaEEIIsVIUpAkhhBArRUGaEEIIsVIUpAkhhBArRUGaEEIIsVIUpAkB8Ndff2HevHl4+vSppbNCiumPP/6Ag4MDbt++LdheWFiI9evXo0uXLqhcuTLkcjk8PDwQGhqKL7/8Eo8fPzbqPHl5eahSpQpCQ0N1piksLISfn59yidPDhw/D1dUVd+/eNf6DkXKJgjQh4IL0/PnzKUiXcowxTJw4EWPGjEGNGjWU27Ozs9G1a1cMHToUnp6e+Pbbb3H48GGsX78eHTp0wBdffIGIiAijziWXy/HWW2/h1KlTuHz5smiaQ4cO4c6dOxg1ahQAoGPHjmjRooXWevaE6GT08ieElEFffPGF1ipE+jx//ty8GbJihYWFVvv5+bXzr169Ktj+9ttvMwAsOjpadL9nz56xH3/80ejzXb58mQFgkydPFn1/wIABzN7eXrDy15YtW5itrS1LTk42+nyk/KEgTcq9uXPnii4bGBcXxxjjlhbs0aMH27p1KwsJCWEODg7s448/Zjdv3mQA2Nq1a7WOCZHF/P/99182cOBAVqlSJWZvb8/8/f3Z//73P0l5LCwsZMuXL2fBwcHM0dGReXh4sP79+7OkpCRBunbt2rHGjRuz06dPs9dee405OTmxWrVqsUWLFrGCggJB2vT0dDZ58mRWs2ZNJpfLmUKhYBMmTGBZWVlan2XcuHHshx9+YP7+/kwul7MffviBMcbYH3/8wUJDQ5mDgwNTKBRs1qxZbNWqVYIbnpEjR7KKFSsq10pWFxYWxho1aiSpDKTo1asXa968uWDbvXv3mJ2dHevRo4dRx5Ja5q1atWKVK1dmeXl5gu1Pnjxhjo6OLDIyUrA9NzeXubu7s9mzZxuVH1I+UZAm5d6dO3fYBx98wACwbdu2KR8Gkp6ezhjjgnS1atVY7dq12Zo1a1hcXBw7ffq0UUH60qVLzN3dnQUGBrKff/6ZHTx4kE2ePJnZ2NiwefPmGczjmDFjmFwuZ5MnT2YHDhxg0dHRzN/fn1WpUoX9999/ynTt2rVjXl5erF69emzFihUsNjaWvffeewwAW7dunTLds2fPWEhICPP29mZfffUVO3ToEFu2bBlzd3dnHTp0YIWFhYLP4uPjw4KCglh0dDQ7cuQIS0xMZBcuXGCOjo4sKCiIbdy4ke3atYt1796d1axZUxCkL1y4wACwVatWCT7TpUuXGAC2fPlyKX8mg3Jzc5mTkxObOnWqYPuGDRsYALZy5Uqjjie1zH/66SfB06Z4y5cvZwDYgQMHtI7drVs39sorrxiVH1I+UZAmhOlv7q5RowaztbVl//zzj2C7MUG6S5cuzNfXVxn4ee+//z5zdHRkaWlpOvN24sQJBoAtXbpUsP3OnTtaQYl/wMmpU6cEaRs1asS6dOmifL1o0SJmY2PD4uPjBem2bNnCALB9+/YJPou7u7tWHiMjI5mLiwt79OiRcltBQQFr1KiRVlm2a9eOhYSECPYfO3Ysc3Nz03rEYlGdOnVK9DGxixcv1hks8/LyBD88Y8o8MzOTubq6st69ewvSNm3alFWvXl2rBYMxxmbOnMlsbGy0Wi0I0UQDxwiRICgoCPXr1y/Svjk5OTh8+DAiIiLg7OyM/Px85U/37t2Rk5ODkydP6tx/z549kMlkGDJkiGDfqlWrIjg4GEePHhWkr1q1Klq0aKGVf/XRznv27EFAQABCQkIEx+zSpQtkMpnWMTt06ICKFSsKth07dgwdOnSAt7e3cpuNjQ3eeOMNrc8wYcIEnD9/HsePHwcAZGRk4JdffsGwYcPg6uqqt/ykunfvHgCgcuXKktKfP38ecrlc8MOP8DamzF1dXfHGG29g3759ePDgAQAgMTERCQkJGD58OGxstC+zlStXRmFhIf77779ifmpS1lGQJkSCatWqFXnf1NRU5Ofn47vvvtMKCt27dwcAvdN/Hjx4AMYYqlSporX/yZMntfb18vLSOoaDgwOys7MFx/z777+1jlehQgUwxrSOKfb5U1NTUaVKFa3tYtv69OmDmjVrYvny5QCAqKgoPHv2DOPGjdP5uY3Ffz5HR0fBdj8/PwDQmpLVoEEDxMfHIz4+HmPGjBG8Z2yZjxo1Cvn5+fjll18AAGvWrIFMJsOIESNE88rnUf1vQogYO0tngJDSQCaTaW3jL7S5ubmC7ampqYLXFStWhK2tLd566y2dQalWrVo6z+3t7Q2ZTKac/6tJbJsh3t7ecHJywpo1a3S+r07s83t5eSlrjurEaoc2NjYYN24cZsyYgaVLl+L7779Hx44d0aBBA6Pzrguf57S0NMH29u3bw87ODrt27cLbb7+t3O7k5IRmzZoB4GrOmscypsxbt26Nhg0bYu3atZgwYYJyapeuvyufR81yJkQTBWlCoLroGlOzqVKlChwdHfH3338Ltu/cuVPw2tnZGWFhYTh37hyCgoJgb29vVN569uyJxYsX4+7du6JNyUXRs2dPfPbZZ/Dy8tJ7g6BPu3btsG/fPjx+/FgZbAoLC7F582bR9KNHj8a8efMwePBg/PPPP/j888+LnH8xDRs2BAAkJSUJtlerVg0jR47Ejz/+iI0bN+LNN980eKyilPnIkSMxZcoUzJo1C48ePcLIkSN1pr1x4wa8vLxEWx0IUUdBmhAAgYGBAIBly5Zh2LBhkMvlaNCgASpUqKBzH77Pcs2aNahTpw6Cg4Nx+vRpREdHa6VdtmwZXnvtNbRp0wZjx45FzZo1kZmZievXr2P37t04cuSIzvO8+uqrePvttzFixAicOXMGbdu2hYuLC+7fv48///wTgYGBGDt2rFGfd+LEidi6dSvatm2LDz/8EEFBQSgsLERycjIOHjyIyZMno2XLlnqPMXPmTOzevRsdO3bEzJkz4eTkhBUrVuDZs2cAoNUX6+HhgaFDh+KHH35AjRo10KtXL6PybIivry9q166NkydPYvz48YL3vvnmG9y8eRODBw/Grl270KdPHygUCjx//hxXr17Fxo0b4ejoCLlcDqBoZT506FDMmDEDX3zxBTw8PNCvXz+deT158iTatWsn2kJBiICFB64RYjWmT5/OFAoFs7GxEZ0nLSY9PZ2NHj2aValShbm4uLBevXqxW7duic6TvnnzJhs5ciTz8fFhcrmcVapUibVu3ZotWLBAUv7WrFnDWrZsyVxcXJiTkxOrU6cOGzp0KDtz5owyDT9PWtOwYcNYjRo1BNuysrLYrFmzWIMGDZi9vb1yitiHH34omGKEl/Okxfzxxx+sZcuWzMHBgVWtWpVNmTKFff755wwAe/r0qVb6o0ePMgBs8eLFkj6zsWbPns0qVqzIcnJytN4rKChgP//8MwsPD2fe3t7Mzs6Oubu7sxYtWrDZs2ezlJQUrX2klLm6iIgIBoC99957OvN4/fp1BoBt3bq16B+UlBsyxhiz4D0CIaSM6dy5M27duoV///1X673Jkyfjhx9+wJ07d0QHuBXXvXv3UKtWLfz8888YMGCAyY9vCrNnz8bPP/+MpKQk2NlRYybRj74hhJAimzRpEpo0aYLq1asjLS0NGzZsQGxsLFavXi1Id/LkSfz777/4/vvv8c4775glQAOAQqHAxIkTsXDhQkRGRopOf7Kkp0+fYvny5fjuu+8oQBNJ6FtCCCmygoICzJkzB//99x9kMhkaNWqEX375BUOGDBGka9WqFZydndGzZ08sWLDArHmaNWsWnJ2dcffuXVSvXt2s5zLWzZs3MX36dAwaNMjSWSGlBDV3E0IIIVbKutqCCCGEEKJEQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECtFQZoQQgixUhSkCSGEECv1f+RTyPxwEWccAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# make plots\n", "matplotlib.rcParams.update({'font.size':12})\n", "fig, ax = plt.subplots(1,1)\n", "plt.gcf().subplots_adjust(bottom=0.15)\n", "plt.gcf().subplots_adjust(left=0.15)\n", "ax.set_xlim((0.,1000.))\n", "ax.set_ylim((0.,1000.))\n", "x0,x1 = ax.get_xlim()\n", "y0,y1 = ax.get_ylim()\n", "ax.set_aspect(abs(x1-x0)/abs(y1-y0)) # make square plot\n", "xtick_spacing = 200\n", "ytick_spacing = 200\n", "ax.yaxis.set_major_locator(ticker.MultipleLocator(xtick_spacing))\n", "ax.yaxis.set_major_locator(ticker.MultipleLocator(ytick_spacing))\n", "plt.xlabel('true energy (GeV)', labelpad=3)\n", "plt.ylabel('reconstructed energy (GeV)', labelpad=3)\n", "plt.scatter(y_test, y_pred, s=3, color='dodgerblue', marker='o', label=model_label)\n", "plt.figtext(0.3, 0.81, model_label)\n", "plt.figtext(0.3, 0.73, f'$R^2$ = {R2:.3f}')\n", "plt.show()\n", "plt.savefig(\"scatterplot.pdf\", format='pdf')\n", "\n", "fig, ax = plt.subplots(1,1)\n", "ax.set_xlim((0.,1000.))\n", "ax.set_ylim((-1,1.))\n", "x0,x1 = ax.get_xlim()\n", "y0,y1 = ax.get_ylim()\n", "ax.set_aspect(abs(x1-x0)/abs(y1-y0)) # make square plot\n", "xtick_spacing = 200\n", "ytick_spacing = 0.4\n", "plt.grid(True)\n", "ax.yaxis.set_major_locator(ticker.MultipleLocator(xtick_spacing))\n", "ax.yaxis.set_major_locator(ticker.MultipleLocator(ytick_spacing))\n", "plt.xlabel('true energy (GeV)', labelpad=3)\n", "plt.ylabel(r'$r = (\\hat{E}-E)/E$', labelpad=3)\n", "residual = (y_test - y_pred)/y_pred\n", "mean_resid = residual.mean()\n", "sigma_resid = residual.std()\n", "rootMSE = np.sqrt(sigma_resid**2 + mean_resid**2)\n", "plt.scatter(y_test, residual, s=3, color='dodgerblue', marker='o')\n", "plt.figtext(0.25, 0.81, model_label)\n", "plt.figtext(0.55, 0.28, f'$\\\\langle r \\\\rangle$ = {mean_resid:.3f}')\n", "plt.figtext(0.55, 0.22, f'$\\sigma_r$ = {sigma_resid:.3f}')\n", "plt.figtext(0.55, 0.16, f'RMS err. = {rootMSE:.3f}')\n", "plt.show()\n", "plt.savefig(\"resolution.pdf\", format='pdf')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.7.12" } }, "nbformat": 4, "nbformat_minor": 4 }