{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import mlatom as ml\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# prepare H2 geometries with bond lengths ranging from 0.5 to 5.0 Å\n",
    "xyz = np.zeros((451, 2, 3))\n",
    "xyz[:, 1, 2] = np.arange(0.5, 5.01, 0.01)\n",
    "z = np.ones((451, 2)).astype(int)\n",
    "molDB = ml.molecular_database.from_numpy(coordinates=xyz, species=z)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# calculate HF energies\n",
    "hf = ml.models.methods(method='HF/STO-3G', program='PySCF')\n",
    "hf.predict(molecular_database=molDB, calculate_energy=True)\n",
    "molDB.add_scalar_properties(molDB.get_properties('energy'), 'HF_energy') # save HF energy with a new name"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# calculate CISD energies\n",
    "cisd = ml.models.methods(method='CISD/cc-pVDZ', program='PySCF')\n",
    "cisd.predict(molecular_database=molDB, calculate_energy=True)\n",
    "molDB.add_scalar_properties(molDB.get_properties('energy'), 'CISD_energy')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAGyCAYAAADnH8C6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAABpF0lEQVR4nO3dd3gU1R7G8e+m9w2QBiQQekeRGkSKIEWl2BFpoqAC4hXbxYoKolhQrwoqCljBAlJEBBFQpEkJSAu9k4QQSEgvO/ePhdUYCElIMinv53n2yU7b+W0Ssi9nzpxjMQzDQERERKSCcjK7ABEREREzKQyJiIhIhaYwJCIiIhWawpCIiIhUaApDIiIiUqEpDImIiEiFpjAkIiIiFZqL2QWUdjabjRMnTuDr64vFYjG7HBEREckHwzA4d+4c1apVw8kp77YfhaHLOHHiBGFhYWaXISIiIoVw9OhRQkND89xHYegyfH19Afs308/Pz+RqREREJD8SExMJCwtzfI7nRWHoMi5cGvPz81MYEhERKWPy08VFHahFRESkQlMYEhERkQpNYUhEREQqNPUZKiLZ2dlkZmaaXYaUQa6urjg7O5tdhohIhaUwdIUMwyA6OpqzZ8+aXYqUYf7+/oSEhGgsKxEREygMXaELQSgoKAgvLy99mEmBGIZBSkoKsbGxAFStWtXkikREKh6FoSuQnZ3tCEJVqlQxuxwpozw9PQGIjY0lKChIl8xEREqYOlBfgQt9hLy8vEyuRMq6C79D6ncmIlLyFIaKgC6NyZXS75CIiHkUhkRERKRCUxgSERGRCq3MhKEzZ84waNAgrFYrVquVQYMG5et29l27dtGnTx+sViu+vr60a9eOI0eOFH/BIiIiUiaUmTA0YMAAIiMjWbJkCUuWLCEyMpJBgwblecz+/fvp0KEDDRs2ZOXKlWzdupXnnnsODw+PEqq69Bo6dCj9+vXLtX7lypVYLBbOnj3reP7vx7PPPpvrmKpVq2IYBrGxsTzwwAPUqFEDd3d3QkJC6NGjB2vXrr3k6/3zMXPmTMB+p96UKVNo3rw5Hh4e+Pv706tXL/7444/LvrfTp0/Ts2dPqlWrhru7O2FhYYwePZrExMQc+xmGwccff0xERAR+fn74+PjQpEkTHnnkEfbt21fo762IiORPVraN+OQMEtPMvXmkTNxav2vXLpYsWcK6deto27YtgONDLCoqigYNGlz0uGeeeYYbb7yRyZMnO9bVrl27RGouT6KiovDz83Ms+/j45Ni+YMEC+vTpg8Vi4bbbbiMzM5NZs2ZRu3ZtYmJiWL58OfHx8dxwww2cPHnScdwjjzxCYmIiM2bMcKyzWq0YhkH//v355ZdfeP311+natSuJiYm8//77dO7cmW+//faiQe4CJycn+vbty4QJEwgMDGTfvn2MGjWK+Ph4vvrqK8AehAYMGMAPP/zA008/zZQpUwgKCuLgwYMsW7aMCRMmOIKZiIjkLSPLRkJqJgmpGZxNybQ/UjM5m5JBQmomZ1Ls6xNSL2yzL59LywLg0W71eaRbPdPqLxNhaO3atVitVkcQAmjXrh1Wq5U1a9ZcNAzZbDZ+/PFHnnzySXr06MGWLVuoVasW48aNy/ODND09nfT0dMfyv1sTLscwDFIzswt0TFHwdHUutjuSgoKC8Pf3v+T2BQsW8M4773D27FlWr17NypUr6dSpEwA1a9akTZs2jn1DQkL+rtnTk/T09BzrAObMmcN3333HggUL6N27t2P9Rx99xOnTp7n//vu54YYb8Pb2vmg9lSpV4qGHHnIs16xZk5EjR/L666/nOMfs2bOZP38+ffr0cayvXbs2Xbt2xTCMy3xXRETKp7TMbE4nZ3AmOSPX1/iUDBL+EWbsjwySM67scy8lM6uIqi+cMhGGoqOjCQoKyrU+KCiI6Ojoix4TGxtLUlISr776KhMmTOC1115jyZIl3HrrraxYscLxYf1vkyZN4sUXXyx0ramZ2TR+/udCH19YO1/qgZdbyf84d+zYQXR0NF27dsXFxQUfHx9++OEH2rVrh7u7e6Fe86uvvqJ+/fo5gtAFjz32GHPnzmXZsmV5htp/OnHiBHPnzs3xM//6669p0KBBjiD0T7rVXUTKA5vNICE1k/iUDOKTcz8cQSclg9NJ9q8phQw2Fgv4ebji7+WKv6crVi83Kv3jub/n+W1erlg93fD3cqWSlxt+Hi64OJvba8fUMDR+/PjLBo8///wTuPiHk2EYl/zQstlsAPTt25dHH30UgKuvvpo1a9Ywbdq0S4ahcePGMXbsWMdyYmIiYWFhl38zZdCiRYtyXfLKzs79jyA0NDTH8uHDhx0jbs+fP58ePXo4+mHNnDmT4cOHM23aNK655ho6depE//79ad68eb7r2rNnD40aNbrotgvr9+zZc9nXufvuu5k/fz6pqan07t2b6dOn5zjHv1sU//Of/zj28ff359ixY/muWUSkpGTbDE4npxN3LoNTSenEnUsnLimdUxe+Jtm3xSWlcyYlA1shGrpdnS1U8nKjsnfOR6ULAcfLDev5MHMh5Ph6uOLsVDb/I2lqGBo9ejT9+/fPc5/w8HC2bdtGTExMrm2nTp0iODj4oscFBATg4uJC48aNc6xv1KgRq1evvuT53N3dC92iAfbLVTtf6lHo46/kvAXVpUsXpk6dmmPd+vXrGThwYI51v//+O76+vo7lSpUqOZ7Pnz+fkSNHOpZvu+02brrpJn7//XfWrl3LkiVLmDx5MtOnT2fo0KEFrvFSLoTgXr168fvvvwP2y2E7duxw7DNlyhReeOEFoqKiePrppxk7diwffPBBrte44JlnnmH06NHMnTuXV155pchqFRG5HMMwSEzNIjoxjdhzaX8Hm3PpxCVl5Ag7p5MzKOiVfF93Fyr72MNMFW83Knn//bWytxuVvdyo7PP3Ol93lwrVQm5qGAoICCAgIOCy+0VERJCQkMCGDRsc/U/Wr19PQkIC7du3v+gxbm5utG7dmqioqBzr9+zZQ82aNa+8+EuwWCymXK4qDG9vb+rWrZtj3cVaQ2rVqnXRPkPR0dFs3ryZm266Kcd6Dw8PbrjhBm644Qaef/557r//fl544YV8h6H69euzc+fOi27btWsXAPXq2TvaTZ8+ndTUVABcXV1z7BsSEkJISAgNGzakSpUqXHfddTz33HNUrVqVevXqsXv37hz7BwYGEhgYeNFLsiIihZWWmU1MYhoxien2sJOYRnRCGjHn0olJSCPmXBoxiWmkZdry/ZoWC1TxdiPAx50AH3cCfd0J8HH7x3P7o8r5AOTmUmZuHjdFmfjUbtSoET179mT48OF8+OGHAIwYMYKbb745x6WOhg0bMmnSJG655RYAnnjiCe666y46duxIly5dWLJkCQsXLmTlypVmvI1yZ8GCBURERFw20DZu3Jgffvgh36/bv39/BgwYwMKFC3P1G3rzzTepUqUKN9xwAwDVq1fP12te6BB9oXP83XffzYABA5g/fz59+/bNd20iIv+UmpHNiYRUTpy1P46fTePk2dTzoccefhJS83/buL+XK0G+9kATmCPouBNwYZ2vvSXH7H425UmZCEMAX375JWPGjKF79+4A9OnTh/feey/HPlFRUSQkJDiWb7nlFqZNm8akSZMYM2YMDRo04Pvvv6dDhw4lWnt5tWDBghxB4vTp09xxxx0MGzaM5s2b4+vry8aNG5k8eXKBAkf//v359ttvGTJkSK5b6xcsWMC33357yTvJABYvXkxMTAytW7fGx8eHnTt38uSTT3LttdcSHh7uOMfcuXPp378/48aNo0ePHgQHB3P48GHmzJmjmeNFBJvNIC45nRNn0/4Rdi4EnzSOn00lPjkjX6/l4epEiJ8HQX4ehPh5EOznTrCfh+Nh3+aORyG6PMiVKzNhqHLlynzxxRd57nOx26GHDRvGsGHDiqusCis5OZnly5czZcoUxzofHx/atm3LlClT2L9/P5mZmYSFhTF8+HCefvrpfL+2xWLhm2++4Z133mHKlCmMGjUKd3d3IiIiWLFixWXDrKenJx9//DGPPvoo6enphIWFceutt/Lf//43xznmzJnDxx9/zIwZM5g8eTKZmZmEhobStWtX3nrrrYJ/U0SkzElIzeRofApH4lM4fNr+9Wh8CkfPpHDybBoZ2Ze/dOXt5kz1Sp5U8z//sJ4PONa/w46fR8Xqg1PWWAwNqJKnxMRErFYrCQkJOQYeBEhLS+PgwYPUqlWrwo1qPXfuXJ599tlL9u2RgqnIv0sixSkr28bJhLRcYefI+cflLmFZLBDs60E1fw+q+XtSvZIn1f09qWa1B5/q/p74eSrolEZ5fX7/W5lpGZLSxcfHh9dee83sMkREMAyD08kZHDiVzMG4JA7EJZ9/nszh08lkZuf9f/4AH3dqVPakRmUvalT2Iuz8o7q/JyFWD1zVN6fcUxiSQrnQd0tEpKRkZNk4EJfEvtgkDp5KtoeeuGQOnkoiMe3SIxi7uTgRVunvsFOjirfjeWglT7zd9VFY0ek3QERESpXMbBuH4pLZE5PEnphz7I09x56YJA7GJZN9iREELRao7u9JrQBv6gT6UCvAm9qB3tQK8Kaa1ROnMjoYoJQMhSERETGFYRicSEhj54lEdp1MtAefmCQOxCVd8tKWr7sLdYN9qB3gQ+1Ab2oHeFM70IeaVbx0J5YUmsKQiIgUu8xsG/tPJbHzRCI7TySy40QiO08mXrIDs7ebM3WDfakf5EP9YF/qh/hSP9iHED8PdVaWIqcwJCIiRSojy0ZU9Dkij51l+7EEdp5MJCrmHBlZuW9Td3GyUDfIh8ZV/RyBp36wry5tSYlSGBIRkUKz2QwOnU5m67GzbD2awNZjZ9lxIvGiwcfH3YXGVf1oXO38o6of9YJ9cHfR5S0xl8KQiIjk29mUDDYdPsPmI2fYdiyBrUfPXvROLqunK81DrVwV6k+T8+EnrJKXWnukVFIYEhGRizIMg2NnUvnzUDx/HjrDpsPx7IlJyrWfu4sTTatbaR5q5eowf64K9admFS/17ZEyQyNJVWDR0dE8/PDD1K5dG3d3d8LCwujduzfLly8HIDw8nLffftux/5YtW7j55psJCgrCw8OD8PBw7rrrLuLi4gA4dOgQFovF8fD19aVJkyaMGjWKvXv3XrSGmTNn0q5du2J/rxfTrFkz7r///otu+/rrr3F1dSUmJoaVK1c63pOTkxNWq5UWLVrw5JNPcvLkSccx/37//37UqlWrpN6aSKHYbAY7TyQy84+DjPpqM+0mLee6ySsY+81Wvt5wxBGEagd6c0fLUCbe0pRFD3dg+4s9+P6h9rzQuwl9r65OeIC3gpCUKWoZqqAOHTrEtddei7+/P5MnT6Z58+ZkZmby888/M2rUKHbv3p1j/9jYWLp160bv3r35+eef8ff35+DBgyxYsICUlJQc+/7yyy80adKElJQU/vrrL9555x2uuuoqFi5cSNeuXXPs++/JXkvSfffdx/PPP8+7776Ll5dXjm2ffvopN998M8HBwezatQuwTwTs5+dHYmIimzdvZvLkyXzyySesXLmSZs2aERYWliMcXbBx40b69evHqFGjSuR9ieSXYRgcOp3CH/viWLv/NGsPnM418aiLk4VmoVZah1emZc1KtKpZiSo+7iZVLFJMDMlTQkKCARgJCQm5tqWmpho7d+40UlNTTajsyvTq1cuoXr26kZSUlGvbmTNnDMMwjJo1axpTpkwxDMMw5s2bZ7i4uBiZmZmXfM2DBw8agLFly5Yc67Ozs43OnTsbNWvWNLKyshzrU1NTDW9vb2P79u2GYRhGWlqa8cQTTxihoaGGm5ubUbduXWP69OmO/bdv327ceOONhq+vr+Hj42N06NDB2Ldv3yXr6dSpkzFq1Chj1KhRhtVqNSpXrmw888wzhs1mMwzDMOLi4gw3Nzdj5syZOY47fPiw4eTkZCxcuNAwDMNYsWKFATi+LxekpKQYDRo0MK699tpL1hAdHW2EhoYa99xzzyX3ufC9KKu/S1K2nDibYny38agxdk6k0e6VX4yaTy3K8Wj03E/GoE/WG/9bvsdYuz/OSEnPuvyLipRCeX1+/5tahoqaYUBmyuX3K2quXvYhWPMhPj6eJUuWMHHiRLy9vXNt9/f3z7UuJCSErKws5s2bx+23316gJnAnJyceeeQRbrnlFjZt2kSbNm0AWL58OSEhITRp0gSAwYMHs3btWt59912uuuoqDh486LgEd/z4cTp27Ejnzp359ddf8fPz448//iAr69JD8APMmjWL++67j/Xr17Nx40ZGjBhBzZo1GT58OFWqVKFv377MmDGDIUOGOI6ZMWMGwcHB9OrVK8/X9vT05MEHH+TRRx8lNjaWoKCgHNszMzO57bbbCAkJYfr06fn+fokUpfSsbP48eIaVUbGs3HOKfbE5+/y4OTvRooY/19YNoH2dKlwV5q+5uKTCURgqapkp8Eq1kj/v0yfALXewuZh9+/ZhGAYNGzbM98u3a9eOp59+mgEDBvDggw/Spk0brr/+egYPHkxwcPBlj79wrkOHDjnC0Pz58x2XyPbs2cM333zDsmXL6NatGwC1a9d2HP/+++9jtVqZPXs2rq6uANSvX/+y5w0LC2PKlClYLBYaNGjAX3/9xZQpUxg+fDgAw4YN48Ybb+TAgQPUrl0bwzCYOXMmQ4cOxdn58rf7/vN9/TsMjR49mn379rFx40bNRC8l6tiZFFZGnWJlVCxr9p8mJSPbsc3JAs1C/WlfpwrX1gmgZc1KeLrp1nap2BSGKiDDsA9zX9AOjhMnTmTs2LH8+uuvrFu3jmnTpvHKK6/w22+/0axZswKd0zAMFi5cyOzZswGIjIzE2dmZTp06XfT4yMhIrrvuOkcQ+qfff/89RyvOhx9+yD333APYQ9w/32dERARvvvkm2dnZODs70717d0JDQ5kxYwYvv/wyv/76K4cOHeLee+/N1/fkUt/LadOmMXPmTFasWEFoaGi+XkuksGw2gy1Hz7B0RwzLd8fmav0J8nWnU/1AujQM4tq6AVg9c/87EqnIFIaKmquXvZXGjPPmU7169bBYLOzatYt+/foV6DRVqlThjjvu4I477mDSpEm0aNGCN954g1mzZuV53IVOyBfuqNqwYQMZGRl06NABsF9yykte21u1akVkZKRjOT8tVRc4OTkxdOhQZs6cyYsvvsiMGTPo2LEj9erVy9fxF95XeHi4Y93q1asZM2YMH3zwAe3bt893LSIFkZ6Vzdr9p/l5RwzLdsYQl5Tu2ObsZOGaGv50bhBE5waBNK7qp7u7RPKgMFTULJZ8X64yS+XKlenRowfvv/8+Y8aMydVv6OzZsxftN/Rvbm5u1KlTh+Tk5Dz3s9lsvPvuu9SqVYsWLVoA9ktkN910k+NSVLNmzbDZbKxatcpxmeyfmjdvzqxZs8jMzMzVOuTp6UndunUveu5169blWq5Xr16OS2D33nsvEyZMYO7cucydO5dp06Zd9r0DpKam8tFHH9GxY0cCAwMBOHr0KLfddhsjRoy45G37IoWVlJ7Fit2xLN0Zw4rdsSSl/91nztfdhS4Ng7ihcTAd6wVi9VLrj0h+KQxVUBdaLdq0acNLL71E8+bNycrKYtmyZUydOtXR4nHBokWLmD17Nv3796d+/fqOy1yLFy9mxowZOfY9ffo00dHRpKSksH37dt5++202bNjAjz/+6AghCxYs4MUXX3QcEx4ezpAhQxg2bJijA/Xhw4eJjY3lzjvvZPTo0fzvf/+jf//+jBs3DqvVyrp162jTpg0NGjS45Ps8evQoY8eO5YEHHmDz5s3873//480338yxT61atbj++usZMWIErq6u3H777Rd9rdjYWNLS0jh37hybNm1i8uTJxMXFMXfuXADS0tK45ZZbqF69Ov/973+Jjo7O9RohISF5/FREckvLzGbF7lgWbD3Br7tjSf/HNBdBvu7c0DiYHk1CaFe7Cm4u6vgsUhgKQxVUrVq12Lx5MxMnTuSxxx7j5MmTBAYG0rJlS6ZOnZpr/8aNG+Pl5cVjjz3G0aNHcXd3p169ekyfPp1Bgwbl2PdCy46Xlxc1a9akS5cufPTRR47Wm/3797Nv3z569OiR47ipU6fy9NNPM3LkSE6fPk2NGjV4+umnAfvluV9//ZUnnniCTp064ezszNVXX821116b5/scPHgwqamptGnTBmdnZx5++GFGjBiRa7/77ruP5cuXM2LEiFxjDl3QoEEDLBYLPj4+1K5dm+7duzN27FhHwFm/fj2bNm0C7B23L+ZCHyORvGRm21i9N44FW0+wdEc0yf/oAF07wJvuTULo0SSYq0L9Nb2FSBGwGPrrnKfExESsVisJCQn4+fnl2JaWlsbBgwepVauW7hYqgLfeeotffvmFxYsXF+t5OnfuzNVXX51jFO3SSr9LYrMZbDgUz/zIE/y0/SRnUzId26r7e9L7qmr0vqqq+v+I5FNen9//ppYhKXGhoaGMGzfO7DJESoVjZ1L4ftNxvt98jCPxf49RFuDjzs3Nq9L7qmpcU8NfAUikGCkMSYm78847zS5BxFSpGdks2XGSbzceY83+0471Pu4u3NgshL5XV6dd7So46xKYSIlQGJJya+XKlWaXIJLD9uMJfLn+CAu3nshxJ1j7OlW4o1UoPZqE4OWmP8siJU3/6kREilFaZjaLtp3ki3WHiTx61rE+rLInt18Txm0tqxNaKf/jhIlI0VMYKgLqgy5XSr9D5c+BU0l8tf4I3246RkKqvTO0q7OFnk2rMqBNDdrWqqw7wURKCYWhK3Bh8L+UlJTLjqAskpeUFHvH2YtNNyJlh81msGrvKT5dfZDf98Y51lf392RA2xrc2SqMQF93EysUkYtRGLoCzs7O+Pv7ExsbC9jH1dEdH1IQhmGQkpJCbGws/v7++ZocVkqftMxs5m05zierDzrmBbNYoEuDIAa2q0Gn+kHqDC1SiikMXaELA+5dCEQiheHv76/RqcugU+fS+XzdYb5Yd5j45AzAfkfYXa3DGNo+nLDK6gskUhYoDF0hi8VC1apVCQoKIjMz8/IHiPyLq6urWoTKmINxyUxbuZ95W46TkW2fHqO6vyf3XhvOXa3D8PXQ5U6RskRhqIg4OzvrA02knIuKPsf7K/axaNsJbOf7vF9Tw5/7OtSmR5NgXJw1N5hIWaQwJCJyGduOneW9X/exdGeMY931DYMY1aUOLWtWNrEyESkKCkMiIpew8VA87/66j9/2nALsnaJ7NQ1hZOe6NK1uNbk6ESkqCkMiIv/y17EE3lgaxarzIcjZyULfq6oxsksd6gb5mlydiBQ1hSERkfP2xJzjraV7WLIjGgAXJwu3twxlZOe61KiiO8NEyiuFIRGp8A6fTubtX/byQ+RxDMN+Oazf1dX5T7d61KzibXZ5IlLMFIZEpMI6m5LBu8v38fm6Q2Rm228P69kkhLHd61M/WJfDRCoKhSERqXDSs7L5fO1h3l2+l8Q0++zx19UL4IkeDWge6m9ucSJS4hSGRKTCMAyDxX9F89qS3RyJt88H1zDEl6dvbETH+oEmVyciZlEYEpEKYceJBF6Yv4ONh88AEOjrzuPd63N7yzDNGyZSwSkMiUi5lpCSyRtLo/hy/WFsBni6OjOiY21GdKyNt7v+BIqIwpCIlFM2m8G3m47y2pIoxySqva+qxtM3NqSq1dPk6kSkNFEYEpFyZ9uxszw3fwdbj54FoF6QDy/2bUL7OgHmFiYipZLCkIiUG8npWby5dA8z1xzEZoCPuwv/6VaPIe3DcdUkqiJyCQpDIlIurIiK5dl52zl+NhWAvldX45kbGxHk52FyZSJS2ikMiUiZFpeUzsuLdjI/8gQA1f09eeXWZnTSrfIikk8KQyJSJhmGwdzNx3n5x52cTcnEyQLDrq3F2O718XLTnzYRyT/9xRCRMif2XBpPz93OL7tiAGhU1Y/Xbmum0aNFpFDKTI/CM2fOMGjQIKxWK1arlUGDBnH27Nk8j7FYLBd9vP766yVTtIgUucV/naTHlN/4ZVcMrs4WnujRgAWjr1UQEpFCKzMtQwMGDODYsWMsWbIEgBEjRjBo0CAWLlx4yWNOnjyZY/mnn37ivvvu47bbbivWWkWk6J1NyeCFBTscfYMaVfVjyl1X0TDEz+TKRKSsKxNhaNeuXSxZsoR169bRtm1bAD7++GMiIiKIioqiQYMGFz0uJCQkx/L8+fPp0qULtWvXLvaaRaTorIyK5cnvthF7Lh1nJwsjO9fh4evr4eZSZhq3RaQUKxNhaO3atVitVkcQAmjXrh1Wq5U1a9ZcMgz9U0xMDD/++COzZs3Kc7/09HTS09Mdy4mJiYUvXESuSHpWNq/9FMWnfxwEoHagN2/deTVXh/mbW5iIlCtlIgxFR0cTFBSUa31QUBDR0dH5eo1Zs2bh6+vLrbfemud+kyZN4sUXXyxUnSJSdA6cSuLhr7ew44T9PyRD24fz314N8XB1NrkyESlvTG1jHj9+/CU7OV94bNy4EbB3hv43wzAuuv5iPv30U+655x48PPIegG3cuHEkJCQ4HkePHi34GxORQjMMg283HuXm/61mx4lEKnu78cmQVozv00RBSESKhaktQ6NHj6Z///557hMeHs62bduIiYnJte3UqVMEBwdf9jy///47UVFRzJkz57L7uru74+7uftn9RKTonUvL5Nkftjs6SUfUrsLb/a8mWKNIi0gxMjUMBQQEEBBw+YkTIyIiSEhIYMOGDbRp0waA9evXk5CQQPv27S97/CeffELLli256qqrrrhmESkeO08k8tCXmzh8OgVnJwtjb6jPg53q4OyUv9ZfEZHCKhO3YjRq1IiePXsyfPhw1q1bx7p16xg+fDg333xzjs7TDRs2ZN68eTmOTUxM5Ntvv+X+++8v6bJFJJ++33SMWz74g8OnU6ju78k3D7RjVJe6CkIiUiLKRAdqgC+//JIxY8bQvXt3APr06cN7772XY5+oqCgSEhJyrJs9ezaGYXD33XeXWK0ikj/pWdm8tHAnX64/AkCn+oG80/9q/L3cTK5MRCoSi2EYhtlFlGaJiYlYrVYSEhLw89PgbiJF5fjZVEZ+uZmtR89iscCY6+vxSNd6OKk1SESKQEE+v8tMy5CIlB+r98YxZvYW4pMzsHq68nb/q+nSIPfwGSIiJUFhSERKjGEYzPjjEBN+3InNgKbV/Zh6T0vCKnuZXZqIVGAKQyJSIjKybDw/fzuz/7SP3XV7y1Am9GuqsYNExHQKQyJS7OKTM3joi02sPxiPkwWevrER93Wole9BU0VEipPCkIgUqz0x57hv1p8cjU/Fx92F/93dgi4N1T9IREoPhSERKTa/7o5hzNeRJKVnUaOyF58MaUW9YF+zyxIRyUFhSESKxefrDvPC/O3YDGhbqzJTB7aksrfGDxKR0kdhSESKlGEYTP45iqkr9wNwV6swXu7XFDeXMjHgvYhUQApDIlJkMrJsPPX9NuZtOQ7Ao93qM6ZrXXWUFpFSTWFIRIpEYlomD32xiT/2ncbZycKkW5txZ6sws8sSEbkshSERuWLRCWkMnbGB3dHn8HJzZurAlnSqH2h2WSIi+aIwJCJX5MCpJAZ9soHjZ1MJ8HFn5r2taVrdanZZIiL5pjAkIoW280Qigz9dT1xSBrUDvJk1rI2m1hCRMkdhSEQKZfORMwz9dAOJaVk0rurHZ/e1IcDH3eyyREQKTGFIRArsj31xDP9sIykZ2bSsWYlPh7bG6ulqdlkiIoWiMCQiBbJsZwyjvtpMRpaN6+oF8OGglni56U+JiJRd+gsmIvk2P/I4Y7/ZSrbNoHvjYP43oAXuLpp1XkTKNoUhEcmX7zcd4/HvtmIYcGuL6ky+vTkuzhpVWkTKPoUhEbmsfwahAW1rMKFvU5ycNKq0iJQP+m+diORp7mYFIREp3xSGROSS5m05xmPf2oPQ3W0UhESkfFIYEpGL+mHLcR775u8gNLGfgpCIlE8KQyKSyw9bjjP2m0hsBtzdJkxBSETKNYUhEcnhp79OOoJQ/9ZhTOzXTEFIRMo1hSERcVi15xRjZm/BZsAdLUN55RYFIREp/xSGRASAjYfieeDzjWRmG9zUrCqv3tZcQUhEKgSFIRFhx4kE7p35J2mZNjrVD2TKXVfjrCAkIhWEwpBIBXfgVBKDP9nAubQsWodXYtrAlri56E+DiFQc+osnUoEdP5vKwOnrOZ2cQZNqfnwytDWebpprTEQqFoUhkQoqPjmDQdPXcyIhjTqB3nw2rA1+Hq5mlyUiUuIUhkQqoNSMbO6b9ScH4pKp7u/J5/e1pYqPu9lliYiYQmFIpILJthmMmb2FLUfOYvV0ZdawNlTz9zS7LBER0ygMiVQghmEwfsEOlu2Mwc3FielDWlE3yMfsskRETKUwJFKBTFt1gM/XHcZigXfuuprW4ZXNLklExHQKQyIVxA9bjvPakt0APHdTY3o1q2pyRSIipYPCkEgFsGZfHE98txWA4dfVYliHWiZXJCJSeigMiZRze2PO8cDnm+zTbDSvyrhejcwuSUSkVFEYEinH4pMzuG/WRs6lZ9EmvDJv3nGV5hsTEfkXhSGRciojy8aDX2ziSHwKYZU9mTaoJR6uGl1aROTfFIZEyiHDMHjuh+1sOBiPj7sLnwxpTWVvN7PLEhEplRSGRMqhT1YfZM7GozhZ4H8DWlA/2NfskkRESi2FIZFyZsXuWF5ZvAuAZ25qTJcGQSZXJCJSuikMiZQje2LO8fDXW7AZ0L91GMOuDTe7JBGRUk9hSKScOJOcwX2z/iQpPYu2tSrzUt+mWCy6c0xE5HIUhkTKgQuTrx6NT6VGZS+mDWyJm4v+eYuI5If+WoqUA28ti+L3vXF4ujrz4aCWVNKdYyIi+eZS0AMOHTrE77//zqFDh0hJSSEwMJAWLVoQERGBh4dHcdQoInlYsj2a91fsB+DV25rRqKqfyRWJiJQt+Q5DX331Fe+++y4bNmwgKCiI6tWr4+npSXx8PPv378fDw4N77rmHp556ipo1axZnzSJy3r7YJB7/1j7n2H0datH36uomVyQiUvbkKwxdc801ODk5MXToUL755htq1KiRY3t6ejpr165l9uzZtGrVig8++IA77rijWAoWEbtzaZk88PlGR4fp//ZqaHZJIiJlUr76DL388sts3LiR0aNH5wpCAO7u7nTu3Jlp06axa9cuwsPDi7pOzpw5w6BBg7BarVitVgYNGsTZs2fzPCYpKYnRo0cTGhqKp6cnjRo1YurUqUVem0hJMwyDx7/dyv5TyYT4efDegGtwdVYXQBGRwshXy9BNN92U7xcMCAggICCg0AVdyoABAzh27BhLliwBYMSIEQwaNIiFCxde8phHH32UFStW8MUXXxAeHs7SpUsZOXIk1apVo2/fvkVeo0hJmbpqPz/viMHN2YmpA68h0Nfd7JJERMqsQv1Xcv/+/Tz77LPcfffdxMbGArBkyRJ27NhRpMVdsGvXLpYsWcL06dOJiIggIiKCjz/+mEWLFhEVFXXJ49auXcuQIUPo3Lkz4eHhjBgxgquuuoqNGzcWS50iJWHN/jje+Nn+ez++TxNa1KhkckUiImVbgcPQqlWraNasGevXr2fu3LkkJSUBsG3bNl544YUiLxDsocZqtdK2bVvHunbt2mG1WlmzZs0lj+vQoQMLFizg+PHjGIbBihUr2LNnDz169LjkMenp6SQmJuZ4iJQWp86l88jsSGwG3NEylLvbhJldkohImVfgMPTf//6XCRMmsGzZMtzc/h7LpEuXLqxdu7ZIi7sgOjqaoKDc8ysFBQURHR19yePeffddGjduTGhoKG5ubvTs2ZMPPviADh06XPKYSZMmOfolWa1WwsL0YSOlQ7bN4NE5kZw6l079YB+NMC0iUkQKHIb++usvbrnlllzrAwMDOX36dIFea/z48VgsljwfFy5pXeyPvmEYeX4YvPvuu6xbt44FCxawadMm3nzzTUaOHMkvv/xyyWPGjRtHQkKC43H06NECvSeR4vLBin2s3mcfWPH9Adfg6eZsdkkiIuVCgQdd9Pf35+TJk9SqVSvH+i1btlC9esHGOBk9ejT9+/fPc5/w8HC2bdtGTExMrm2nTp0iODj4oselpqby9NNPM2/ePEcH8ObNmxMZGckbb7xBt27dLnqcu7s77u7qjCqly7oDp5nyyx4AXu7XlHrBviZXJCJSfhQ4DA0YMICnnnqKb7/9FovFgs1m448//uDxxx9n8ODBBXqt/N55FhERQUJCAhs2bKBNmzYArF+/noSEBNq3b3/RYzIzM8nMzMTJKWfjl7OzMzabrUB1ipgpLimdMednor+9ZSi3tww1uyQRkXKlwJfJJk6cSI0aNahevTpJSUk0btyYjh070r59e5599tniqJFGjRrRs2dPhg8fzrp161i3bh3Dhw/n5ptvpkGDBo79GjZsyLx58wDw8/OjU6dOPPHEE6xcuZKDBw8yc+ZMPvvss4te5hMpjWzn+wnFnkunXpAPL/VtYnZJIiLljsUwDKMwB+7fv58tW7Zgs9lo0aIF9erVK+racoiPj2fMmDEsWLAAgD59+vDee+/h7+/v2MdisTBjxgyGDh0K2Dtejxs3jqVLlxIfH0/NmjUZMWIEjz76aL47niYmJmK1WklISMDPT3M+Scl6f8U+Xv85Cg9XJxaM7kB9XR4TEcmXgnx+FzoMZWRkcPDgQerUqYOLS4GvtpUZCkNilk2H47nzw3Vk2wwm396cO1vpzkYRkfwqyOd3gS+TpaSkcN999+Hl5UWTJk04cuQIAGPGjOHVV18tXMUiksO5tEwemR1Jts2g39XVuEP9hEREik2Bw9C4cePYunUrK1euxMPDw7G+W7duzJkzp0iLE6moXpi/g2NnUgmt5MnL/TSekIhIcSrw9a0ffviBOXPm0K5duxx/oBs3bsz+/fuLtDiRimjB1hPM3XIcJwu8fdfV+Hq4ml2SiEi5VuCWoVOnTl10NOjk5GT971XkCh0/m8oz8/4CYPT19WgVXtnkikREyr8Ch6HWrVvz448/OpYvBKCPP/6YiIiIoqtMpIK5MN3GubQsWtTwZ8z1dc0uSUSkQijwZbJJkybRs2dPdu7cSVZWFu+88w47duxg7dq1rFq1qjhqFKkQpq3az4aD8Xi7OfP2XVfj4lzg/6uIiEghFPivbfv27VmzZg0pKSnUqVOHpUuXEhwczNq1a2nZsmVx1ChS7m07dpYpy+zTbYzv04SaVbxNrkhEpOIoUMtQZmYmI0aM4LnnnmPWrFnFVZNIhZKSkcUjsyPJshnc1KyqptsQESlhBWoZcnV1dUx3ISJF47WfdnMwLpkQPw8m3qLb6EVESlqBL5Pdcsst/PDDD8VQikjFs2ZfHLPWHgZg8u3N8fdyM7kiEZGKp8AdqOvWrcvLL7/MmjVraNmyJd7eOfs2jBkzpsiKEynPzqVl8sR32wC4p20NOtYPNLkiEZGKqcBzk9WqVevSL2axcODAgSsuqjTR3GRSXMbN3cbXG44SVtmTJY90xNu9/M7xJyJS0gry+V3gv74HDx4sdGEiYrciKpavNxwF4PXbr1IQEhExUYH7DL300kukpKTkWp+amspLL71UJEWJlGcJKZn893v75bF7rw2nXe0qJlckIlKxFfgymbOzMydPnsw1Jcfp06cJCgoiOzu7SAs0my6TSVEbOyeSuVuOUyvAm8VjrsPTzdnskkREyp2CfH4XuGXIMIyL3vq7detWKlfWPEoieVm6I9oxCesbd1ylICQiUgrku6NCpUqVsFgsWCwW6tevnyMQZWdnk5SUxIMPPlgsRYqUB2eSM3j6/CSswzvWpmXNSiZXJCIiUIAw9Pbbb2MYBsOGDePFF1/EarU6trm5uREeHq6JWkXy8PKincQlZVAvyIdHu9U3uxwRETkv32FoyJAhZGVlAdCtWzdCQzVlgEh+rYyKZe6W41gs8NrtzfFw1eUxEZHSokB9hlxcXBg5cmS56yQtUpyS0rN4Zt52AO5tX4traujymIhIaVLgDtRt27Zly5YtxVGLSLn0xs9RHD+bSmglTx7voctjIiKlTYFHehs5ciSPPfYYx44du+h0HM2bNy+y4kTKuk2H45m19hAAk25thpebBlcUESltCjzOkJNT7sYki8XiuOW+vF1C0zhDUljpWdnc9O5q9sUmcXvLUN644yqzSxIRqTA0HYdIKfD+r/vYF5tEgI87z97UyOxyRETkEgochmrWrFkcdYiUK7tOJvLByv0AvNS3Cf5ebiZXJCIil1LoDgw7d+7kyJEjZGRk5Fjfp0+fKy5KpCzLthn89/ttZNkMujcOplfTELNLEhGRPBQ4DB04cIBbbrmFv/76y9FXCHCMSF3e+gyJFNRnaw+x9VgCvh4uvNyv6UWnrxERkdKjwLfWP/LII9SqVYuYmBi8vLzYsWMHv/32G61atWLlypXFUKJI2RGdkMabS/cA8FTPhgT7eZhckYiIXE6BW4bWrl3Lr7/+SmBgIE5OTjg5OdGhQwcmTZrEmDFjNAaRVGgvLdpBUnoWLWr4M6BNDbPLERGRfChwy1B2djY+Pj4ABAQEcOLECcDesToqKqpoqxMpQ1bsjmXxX9E4O1mY2K8ZTk66PCYiUhYUuGWoadOmbNu2jdq1a9O2bVsmT56Mm5sbH330EbVr1y6OGkVKvdSMbJ6bb59yY9i14TSupjGpRETKigKHoWeffZbk5GQAJkyYwM0338x1111HlSpVmDNnTpEXKFIW/O/XvRw7k0o1qwf/0Yz0IiJlSoHDUI8ePRzPa9euzc6dO4mPj6dSpUq6a0YqpD0x5/jotwMAjO/TBG93TbkhIlKWFMlf7cqVKxfFy4iUOTabwbPztpNlM+jWKJjuTTSmkIhIWVPgMJScnMyrr77K8uXLiY2NxWaz5dh+4MCBIitOpLT7bvMxNhyKx9PVmRf7NjG7HBERKYQCh6H777+fVatWMWjQIKpWrapLY1JhnUnOYNLiXQA8ekM9qvt7mlyRiIgURoHD0E8//cSPP/7ItddeWxz1iJQZry+N4kxKJg1DfLn32lpmlyMiIoVU4HGGKlWqpD5CUuH9dSyBrzccAeClvk1xdS7wPyURESklCvwX/OWXX+b5558nJSWlOOoRKfVsNoPnF2zHMKDf1dVoU0v/ORARKcvydZmsRYsWOfoG7du3j+DgYMLDw3F1dc2x7+bNm4u2QpFS5vvNx9hy5Czebs6Mu7GR2eWIiMgVylcY6tevXzGXIVI2JKRm8tqS3QA80q2eJmIVESkH8hWGXnjhheKuQ6RMePuXPcQlZVAn0Juh7dVpWkSkPCiWXp+GYRTHy4qYKir6HJ+tPQzYR5p2c1GnaRGR8iBff80bNWrEV199RUZGRp777d27l4ceeojXXnutSIoTKS0Mw+CFBdvJthn0bBLCdfUCzS5JRESKSL4uk73//vs89dRTjBo1iu7du9OqVSuqVauGh4cHZ86cYefOnaxevZqdO3cyevRoRo4cWdx1i5SoRdtOsu5APO4uTjx7szpNi4iUJ/kKQ9dffz1//vkna9asYc6cOXz11VccOnSI1NRUAgICaNGiBYMHD2bgwIH4+/sXc8kiJSs5PYuJP9pHmh7VpS6hlbxMrkhERIpSgUagbt++Pe3bty+uWkRKpfdW7CM6MY0alb0Y0bG22eWIiEgRUw9QkTwcjEtm+u/2yYefv7kxHq7OJlckIiJFrcyEoTNnzjBo0CCsVitWq5VBgwZx9uzZPI+JiYlh6NChVKtWDS8vL3r27MnevXtLpmApFyb+uIvMbINO9QPp2ijI7HJERKQYlJkwNGDAACIjI1myZAlLliwhMjKSQYMGXXJ/wzDo168fBw4cYP78+WzZsoWaNWvSrVs3kpOTS7ByKav+2BfHL7ticHay8NzNjXOMwi4iIuVHgWetN8OuXbtYsmQJ69ato23btgB8/PHHREREEBUVRYMGDXIds3fvXtatW8f27dtp0qQJAB988AFBQUF8/fXX3H///SX6HqRsybYZvLxoJwCD2tWkbpCPyRWJiEhxKRMtQ2vXrsVqtTqCEEC7du2wWq2sWbPmosekp6cD4OHx93QJzs7OuLm5sXr16kueKz09ncTExBwPqXjm/HmU3dHnsHq68kjXemaXIyIixajAYahz58589tlnpKamFkc9FxUdHU1QUO7+GkFBQURHR1/0mIYNG1KzZk3GjRvHmTNnyMjI4NVXXyU6OpqTJ09e8lyTJk1y9EuyWq2EhYUV2fuQsiExLZM3l0YB8J9u9ajk7WZyRSIiUpwKHIZatmzJk08+SUhICMOHD2fdunWFPvn48eOxWCx5PjZu3Ahw0f4ahmFcsh+Hq6sr33//PXv27KFy5cp4eXmxcuVKevXqhbPzpe8IGjduHAkJCY7H0aNHC/3+pGx6/9d9nE7OoHagNwPb1TS7HBERKWYF7jP05ptvMnnyZBYtWsSMGTPo2LEjdevWZdiwYQwaNIjg4OB8v9bo0aPp379/nvuEh4ezbds2YmJicm07depUnudr2bIlkZGRJCQkkJGRQWBgIG3btqVVq1aXPMbd3R13d/d8vwcpXw7FJfPpHwcBeO6mxrg6l4krySIicgUsxhXOqnrq1Ck+/PBDJk6cSHZ2NjfeeCNjxozh+uuvL6oa2bVrF40bN2b9+vW0adMGgPXr19OuXTt279590Q7UF7N3714aNmzITz/9RPfu3fN1TGJiIlarlYSEBPz8/Ar9HqRseODzjfy8I4aO9QOZdW9r3UEmIlJGFeTz+4r+27thwwaef/553njjDYKCghg3bhxBQUH07t2bxx9//EpeOodGjRrRs2dPx2W5devWMXz4cG6++eYcQahhw4bMmzfPsfztt9+ycuVKx+31N9xwA/369ct3EJKKZc3+OH7eYb+V/tmbGikIiYhUEAW+TBYbG8vnn3/OjBkz2Lt3L71792b27Nn06NHD8eFx55130q9fP954440iK/TLL79kzJgxjiDTp08f3nvvvRz7REVFkZCQ4Fg+efIkY8eOJSYmhqpVqzJ48GCee+65IqtJyg/7rfT2+ccGtKlB/WBfkysSEZGSUuDLZG5ubtSpU4dhw4YxdOhQAgMDc+2TmJhI3759WbFiRZEVahZdJqsYZm84wn/n/oWvhwurnuhCZd1BJiJSphXk87vALUPLly/nuuuuy3MfPz+/chGEpGI4l5bJG0v3APBI13oKQiIiFUyB+wxdLgiJlDUfrjpAXFI6tQK8GRwRbnY5IiJSwgrcMtSiRYuLdiy1WCx4eHhQt25dhg4dSpcuXYqkQJHiFJ2QxvTV9lnpn+rZEDcX3UovIlLRFPgvf8+ePTlw4ADe3t506dKFzp074+Pjw/79+2ndujUnT56kW7duzJ8/vzjqFSlSby6NIi3TRqualejRJP9jZImISPlR4JahuLg4HnvssVx3ZU2YMIHDhw+zdOlSXnjhBV5++WX69u1bZIWKFLVdJxP5bvMxAJ7WrfQiIhVWgVuGvvnmG+6+++5c6/v3788333wDwN13301UVNSVVydSjF79aTeGATc1q8o1NSqZXY6IiJikwGHIw8PjojPFr1mzxjFDvM1m05QWUqqt3hvHqj2ncHW28GTP/I1gLiIi5VOBL5M9/PDDPPjgg2zatInWre3TFWzYsIHp06fz9NNPA/Dzzz/TokWLIi9WpCjYbAavLLYPsDiwXU1qVvE2uSIRETFToeYm+/LLL3nvvfccl8IaNGjAww8/zIABAwBITU113F1W1mnQxfLn+03HeOzbrfi6u7DqSQ2wKCJSHhXboItZWVlMnDiRYcOGcc8991xyP09Pz4K8rEiJScvM5s2l9hA/sktdBSERESlYnyEXFxdef/11srOzi6sekWL16R8HOZGQRjWrB/deG252OSIiUgoUuAN1t27dWLlyZTGUIlK84pMzmLpiPwCP92iAh6uzyRWJiEhpUOAO1L169WLcuHFs376dli1b4u2ds/Npnz59iqw4kaL07vK9nEvPokk1P/pdXd3sckREpJQocAdqJ6dLNyZZLJZydwlNHajLh0NxyXR7axVZNoMv72/LtXUDzC5JRESKUbHOWm+z2QpdmIhZJv+8myybQecGgQpCIiKSwxXNSpmWllZUdYgUm81HzrD4r2icLDCuVyOzyxERkVKmwGEoOzubl19+merVq+Pj48OBA/YZv5977jk++eSTIi9Q5EoYhsHkJbsBuO2aUBqE+JpckYiIlDYFDkMTJ05k5syZTJ48GTe3v8doadasGdOnTy/S4kSu1O9741h3IB43Fyf+c0N9s8sREZFSqMBh6LPPPuOjjz7innvuwdn571uTmzdvzu7du4u0OJErYbMZTP7Z/js5qF1NqvtrMFAREcmtwGHo+PHj1K1bN9d6m81GZmZmkRQlUhQWbz/J9uOJ+Li7MLJzHbPLERGRUqrAYahJkyb8/vvvudZ/++23mpxVSo3MbBtvLt0DwPDralPFx93kikREpLQq8K31L7zwAoMGDeL48ePYbDbmzp1LVFQUn332GYsWLSqOGkUK7LtNxzgYl0wVbzfuu66W2eWIiEgpVuCWod69ezNnzhwWL16MxWLh+eefZ9euXSxcuJAbbrihOGoUKZC0zGze/sXeKjSqS1183Auc+UVEpAIp1KdEjx496NGjR1HXIlIkPlt7iJjEdKr7e3JPuxpmlyMiIqVcof/LnJGRQWxsbK4RqWvU0IePmCchNZP3z0/G+p9u9XB30WSsIiKStwKHob179zJs2DDWrFmTY71hGOVybjIpWz7+7QAJqZnUC/Lh1mtCzS5HRETKgAKHoaFDh+Li4sKiRYuoWrUqFoulOOoSKbBT59L5ZPVBAB7v0QBnJ/1uiojI5RU4DEVGRrJp0yYaNmxYHPWIFNp7v+4lNTObq8P86d442OxyRESkjCjw3WSNGzcmLi6uOGoRKbSj8Sl8teEIAE/2bKAWSxERybcCh6HXXnuNJ598kpUrV3L69GkSExNzPETMMGXZHjKzDa6rF0D7OgFmlyMiImVIgS+TdevWDYCuXbvmWK8O1GKW3dGJzIs8DsCTPXT5VkRECqbAYWjFihXFUYdIob3xcxSGATc1q0qzUKvZ5YiISBlT4DDUqVOn4qhDpFA2HY7nl12xODtZGNu9vtnliIhIGVTgPkMAv//+OwMHDqR9+/YcP26/PPH555+zevXqIi1OJC+GYfDaT1EA3NEylDqBPiZXJCIiZVGBw9D3339Pjx498PT0ZPPmzaSnpwNw7tw5XnnllSIvUORSVu45xYZD8bi5OPFIt3pmlyMiImVUgcPQhAkTmDZtGh9//DGurq6O9e3bt2fz5s1FWpzIpdhsBq8vsbcKDYmoSVWrp8kViYhIWVXgMBQVFUXHjh1zrffz8+Ps2bNFUZPIZS366yQ7Tybi6+7CyM51zS5HRETKsAKHoapVq7Jv375c61evXk3t2rWLpCiRvGRm23hzqb1VaETH2lTydjO5IhERKcsKHIYeeOABHnnkEdavX4/FYuHEiRN8+eWXPP7444wcObI4ahTJ4ZuNRzl8OoUAHzeGdahldjkiIlLGFfjW+ieffJKEhAS6dOlCWloaHTt2xN3dnccff5zRo0cXR40iDqkZ2bzzy14ARnepi7d7gX+FRUREcrAYhmEU5sCUlBR27tyJzWajcePG+PiUz9uaExMTsVqtJCQk4OfnZ3Y5Fd7Ulft5bcluQit5svyxTri7OJtdkoiIlEIF+fwu9H+rvby8aNWqVWEPFymwhJRMpq6091d7tFt9BSERESkShRp0UcQMH/62n8S0LOoH+9CvRXWzyxERkXJCYUjKhNjEND794yAAT/RoiLOTxeSKRESkvFAYkjLhf7/uIy3TxjU1/OnWKMjsckREpBxRGJJS7/DpZL7ecASAJ3s2xGJRq5CIiBQdhSEp9d5atocsm0HH+oG0q13F7HJERKScURiSUm3niUQWbD0BwJM9GphcjYiIlEcKQ1KqvbE0CsOAm5tXpWl1q9nliIhIOVRmwtDEiRNp3749Xl5e+Pv75+sYwzAYP3481apVw9PTk86dO7Njx47iLVSKzJ+H4vl1dyzOThYe665WIRERKR5lZi6DjIwM7rjjDiIiIvjkk0/ydczkyZN56623mDlzJvXr12fChAnccMMNREVF4evrW8wVy5UwDIPXftoNwJ2twqgV4G1yRSKljGGcf2SDYbM/bBeeZ5/f9u915/czDOD88Rdei389v+hX8thmODZffp/8nuOSb75Qm/I+7nKTMRT22FJ03JUeW5yq1IPgxqadvsyEoRdffBGAmTNn5mt/wzB4++23eeaZZ7j11lsBmDVrFsHBwXz11Vc88MADxVWqFIEVUbFsPHwGdxcnHulaz+xypDwyDMhKh8wU+9estL+/ZmfkXHZ8/fe6dLBlQXYm2DLPP8+yP8/O/Me2C+v+tS3Hsdn/Ci//CDW51p1/LlJedBgLwS+YdvoyE4YK6uDBg0RHR9O9e3fHOnd3dzp16sSaNWsuGYbS09NJT093LCcmJhZ7rZKTzWYweUkUAEPbhxNi9TC5IilVbDZIT4DUs5CWAGlnzz8/v5x+DjKSISPp/Nfkv5fTk3IuG9nmvpcSYQEnZ7A4gcUZLBb7ulxfz+8Leexj+fs1L7WP42Uut09+znGZ95Xn5ssNwXElxxfjuYu17mI+95Xwr1F8r50P5TYMRUdHAxAcHJxjfXBwMIcPH77kcZMmTXK0Qok5Fmw9we7oc/h6uPBQ5zpmlyPFLTMNkk9Bciwkx9mfJ/3jeXIspMT/HXbSEiny5nyLE7h4gos7uHjk/6uzGzi7gJMrOLvaQ4fjuevf25xczq9z+de2fx1ruRBanP4RYC627PyPdZZLHHdhm8blErkcU8PQ+PHjLxs8/vzzzyuaEPbfA/QZhpHnoH3jxo1j7NixjuXExETCwsIKfX4pmIwsG28us7cKPdCxNv5ebiZXJFckIxkSjkHicUg4Dokn7M8T//E8LaFwr+3qBR5W8PAHT/+/v7r7gpsPuHn/42sez1097WFERCosU8PQ6NGj6d+/f577hIeHF+q1Q0JCAHsLUdWqVR3rY2Njc7UW/ZO7uzvu7u6FOqdcudl/HuFofCqBvu4M61DL7HLkcgzD3mpz5iDEHzz/9cDfz5Ni8vc6Tq7gEwTeAeAdCN7/eO4TBJ6VcwYeD6u9dUZEpAiYGoYCAgIICAgolteuVasWISEhLFu2jBYtWgD2O9JWrVrFa6+9ViznlCuTnJ7Fu8v3AjCmaz283MrtVdyyKekUxO6EU7shdtffX9PO5n2cuxWs1cGv2vlH6D+eVwffYHvI0eUcETFJmfm0OXLkCPHx8Rw5coTs7GwiIyMBqFu3Lj4+PgA0bNiQSZMmccstt2CxWPjPf/7DK6+8Qr169ahXrx6vvPIKXl5eDBgwwMR3Ipfy6eqDxCVlULOKF/1b69KkaQzD3qpzIhJObLE/YndCyulLH+NbDSrXsj8qnf9aubb9uad/SVUuIlIoZSYMPf/888yaNcuxfKG1Z8WKFXTu3BmAqKgoEhL+7n/w5JNPkpqaysiRIzlz5gxt27Zl6dKlGmOoFIpPzuDD3w4A8Fj3Brg6l5nxQMu+pFNwdB0c33Q+/EReorXHApXCIagRBDaEoMYQ1BCq1LX3uxERKaMshnHZkaYqtMTERKxWKwkJCfj5+ZldTrk1YdFOpq8+SOOqfix6uANOTrpkUiwMA84cgiNr7Y/Da+H03tz7ObtBcFOo1sL+CGkGAfXBzavESxYRKYyCfH6XmZYhKb+On03ls7X24Q6e7NlAQaiopcTDgZWwfznsX2G/g+vfghpDaGuofo09/AQ2AhfdySciFYPCkJju7WV7yMi20a52ZTrVDzS7nLLPlg3HNsK+ZbD/Vzi+mRzj8ji52ANPjQio2R7C2oJXZdPKFRExm8KQmGpvzDm+33wMgCd7NsxzDCjJQ1YGHPwNdi+E3YvtAxX+U1BjqHO9/VEjQpe7RET+QWFITPX6z1HYDOjRJJhralQyu5yyJTMN9v4MOxfA3qWQ/o+pY9ytUK8b1OkKdbrYb2MXEZGLUhgS02w+coalO2NwssATPRqYXU7ZYLPB4T9g2xx7CEr/x+jNPsHQ8CZoeDOEX6c+PyIi+aQwJKYwDIPXftoNwO0tQ6kbpOEO8hR/ADZ/Btu+ydkB2i8Umt4CjfpA9VbgpCEJREQKSmFITLFqzynWH4zHzcWJ/3Srb3Y5pVNWBkQthk0z7HeDXeBuhSZ9ofldUKO9ApCIyBVSGJISZ7MZTF5in4x1SERNqvlrwL4cEk/Aho9hyxf/6AhtgbpdocUgqN8TXD1MLVFEpDxRGJISt3DbCXaeTMTX3YWRneuaXU7pcWILrP0AdswFW5Z9nU+wPQBdMxgq1TS3PhGRckphSEpURpaNN5fuAWBEx9pU8q7gnXxtNtizBNa+D4dX/72+Zgdo+wA06AXOrubVJyJSASgMSYn6esMRjsSnEODjzrAOtcwuxzw2G+yaD6teh9gd9nVOLtD0Nmg3EqpdbWp5IiIVicKQlJjEtEzeWW6fB+uRbvXwdq+Av362bNg+F35/A07Z76bDzRda32dvCdJ4QCIiJa4CfhqJWaat3E98cga1A73p3zrM7HJKlmHAzvnw64S/J0Z1t0K7h6Ddg+CpASdFRMyiMCQl4sTZVD5ZfRCA//ZsiKtzBbod/NBqWPY8HN9kX/asBBGjoM0I8LCaW5uIiCgMScl4c+ke0rNstAmvzA2Ng80up2TE7IBfxtunygBw9Yb2oyFiNHj4mVqaiIj8TWFIit3OE4nM3WKfjPXpmxqV/8lYU+Lh15dh00wwbPaO0S2HQscnwbeCBEERkTJEYchkWdk2XMrxJSPDMHhl8S4MA3pfVY2rw/zNLqn42LLto0X/OgFSz9jXNe4LXV+AKnXMrU1ERC5JYcgkM/44yIerDnBX6zAevaH8Tkexas8pVu+Lw83ZiSfL82Ssh9fCT09A9F/25aAmcONkCO9gbl0iInJZCkMmcbJYiE5M46/jCZffuYzKthlMWmy/fXxwRE3CKnuZXFExSD1r7xy9eZZ92cMKXZ6FVsPAWf+8RETKAv21NkmzUPtdRNuOncUwjHLZj+b7TceIijmHn4cLo68vh9Nu7FwAix+HpBj78jVD7JfEvKuYW5eIiBSIwpBJGlf1w8XJQlxSBicS0qheziYrTcnI4s1l9slYH76+Hv5e5WjajcQTsPgJ2L3IvlylHvR5F2q2N7cuEREpFIUhk3i4OtMgxJcdJxL569jZcheGPvrtADGJ6YRW8mRw+3IywahhwNbZ8NOTkJ5ov0usw6Nw3eOaRV5EpAwrv7cxlQHNz18q23qsfPUbOnE2lWmr9gPw314NcXdxNrmiIpB8Gr4ZBD88aA9C1VvCA7/B9c8qCImIlHEKQyZqHuoPwF/lLAy9+tNu0jLtAyze1Kyq2eVcuT0/wwftYNdCe2vQ9c/BsKUQ3MTsykREpAjoMpmJmlUvf52oNx6KZ8HWE1gs8HzvxmX7PWUkw89P2wdPBAhsCLd+BFWvMrUsEREpWgpDJmoQ4ou7ixOJaVkcOp1CrQBvs0u6IjabwYsLdwJwV6swmlYvw/NuxeyEb4dCXBRgsc8ldv1zuiQmIlIO6TKZiVydnWhczT5H1bZjZ80tpgh8t/kYfx1PwNfdhce6l9EBFg0DNs2Cj7vYg5BPCAyeDz0mKgiJiJRTCkMma36+9WTr0bLdb+hcWiaTl5y/lb5rXQJ93U2uqBDSEuH7+2HhGMhKgzpd4cHVULuT2ZWJiEgxUhgyWYsalQDYeDje5EquzPsr9hOXlE6tAG+Gtq9ldjkFF7MDPuoE278DizN0Gw/3fAc+gWZXJiIixUx9hkzWtnZlALYfT+BcWia+Hq4mV1Rwh08n8+nqgwA8c2Mj3FzKWMb+6ztY8DBkpoBfKNz+KdRoa3ZVIiJSQsrYp1b5U9XqSY3KXtgM2Hj4jNnlFJhhGIxfsIOMbBvX1Quga6Mgs0vKv+ws+PkZ+P4+exCq3QUe/F1BSESkglEYKgXa1rK3Dm04WPYulf28I4YVUadwdbYwvk+TsnMrfXIcfN4P1r5nX+7wKAz8Hrwqm1qWiIiUPIWhUqDN+TC0/sBpkyspmOT0LF5auAOABzrWoU6gj8kV5dPJbfBhJzj0O7h6w52f2fsIOZWDkbJFRKTA1GeoFGhX2z7L+bZjCaRkZOHlVjZ+LO/+upcTCWmEVvJkVJcyMiv97sX2O8Yyk6FyHej/FQQ1NLsqERExkVqGSoHQSp5Us3qQZTPYfPis2eXky56Yc3zyu73T9It9muDpVspbVQwD/ngXZg+wB6HanWH4rwpCIiKiMFQaWCwW2p5vHVpXBi6VGYbBcz9sJ8tm0K1RMF0bBZtdUt6yMux3iy17DjCg1TD7bfOe/mZXJiIipYDCUCkRUccehn7be8rkSi7v+83HWX8wHg9XJ17o3djscvKWEg9f3ApbPgeLE/R8DW56C5zL3hAGIiJSPMpG55QKoHMD++B+244lEJuYRpBf6Zz6IfZcGi8vss8/NqZrPcIqe5lcUR7iD8AXt0P8fnDzsY8fVL+H2VWJiEgpo5Yhs6Sfg0OrIXYXAEG+HlwVap+aY0VUrJmV5Wn8gh0kpGbStLofI66rbXY5l3ZiC3zS3R6ErDXgvqUKQiIiclEKQ2ZZMQlm3gQbZzhWXeh7s3xX6QxDS7afZPFf0Tg7WXjttua4OJfSX599v8CMmyD5FIQ0g/uXQXATs6sSEZFSqpR+mlUA1VrYv57Y7Fh1fUP76M2r98WRlpltRlWXlJCSybM/2McUerBTbZpUs5pc0SVsnQ1f3WW/Y6xWJxi6GHxDzK5KRERKMYUhs1S/xv41+i/IzgSgSTU/gv3cScnIZn0pG4365R93EpeUTp1Abx6+vp7Z5eRmGLD6bZj3ANiyoNkd9jvGPPzMrkxEREo5hSGzVK4N7lbISnP0G7JYLFzf0H6pbNnOaDOry2FlVCzfbTqGxQKTb2+Oh2spG1PIZoMl4+CXF+zLEaPhlo/Axc3cukREpExQGDKLxQLVrrY//8elsp5N7Zd0Fv8VTUaWzYTCcopPzuCJ77YBMCQinJY1S9ncXdlZ8MODsH6qfbn7ROgxEZz0qy0iIvmjTwwzXbhUdmKLY9W1daoQ6OtOfHIGq/aYO+aQYRiMm7uNU+fSqRvkw397lbLRmrPS4dshsG0OWJzh1o+h/WizqxIRkTJGYchMFzpRH/+7ZcjF2Ym+V1UDYN6WY2ZU5fDtpmP8vCMGV2cLb991dem6PJaRbO8ovXsROLtD/y+h+Z1mVyUiImWQwpCZqp1vGYrdCZlpjtW3XFMdgF92xZKQmmlGZRyKS+bFBfa7x8be0ICm1UvR3WOpZ+HzW+HACvus8/d8Aw16mV2ViIiUUQpDZrKGgneg/e6nmO2O1Y2r+tEg2JeMLBuL/zpZ4mWlZWYz8svNJGdk06ZWZUZ0LEWDKybHwazecHQdeFhh8A/2SVdFREQKqcyEoYkTJ9K+fXu8vLzw9/fP1zFz586lR48eBAQEYLFYiIyMLNYaC8xiueilMovFwq3nW4e+WHcYwzBKtKzxC3aw82QiVbzdeLd/C5ydLCV6/ktKPAEzekH0NvAKgCGLIKyN2VWJiEgZV2bCUEZGBnfccQcPPfRQvo9JTk7m2muv5dVXXy3Gyq7QhUtlxzfmWH1nqzA8XZ3ZcSKRNftLbib7bzceZfafR7FY4J3+LQixlpI50s4chk97Qtwe8KsOw5ZA1eZmVyUiIuVAmZmo9cUXXwRg5syZ+T5m0KBBABw6dKgYKioiNdravx5em2N1JW837mwVyqy1h/notwNcWzeg2EvZcSKB5+bbL9c92q0+HeoV/znz5cwhmHkzJByFSrVg8HyoVNPsqkREpJwoMy1DJSU9PZ3ExMQcj2IV2sZ+W3jCETh7JMemYR1q4WSBVXtOERV9rljLiE5I476ZG0nLtNGpfiCju9Qt1vPlW/xB+zxjCUehSl24d7GCkIiIFCmFoX+ZNGkSVqvV8QgLCyveE7r7/N1v6F+tQzWreDsGYfzfr3uLrYTk9CyGzfyT6MQ06gb58O7dLXAqDf2E4g/YJ7NNPAZV6tn7CPlVM7sqEREpZ0wNQ+PHj8diseT52Lhx4+VfqAiNGzeOhIQEx+Po0aPFf9Ka7e1fD6/OtWlUl7pYLLBo20k2FMN8Zdk2g4e/3sLOk4kE+LgxY2hrrJ6uRX6eAju9394ilHgcAurD0EXgV9XsqkREpBwytc/Q6NGj6d+/f577hIeHl0wx57m7u+Pu7l6i56TmtbDmXTi8JtemJtWs9G8dxtcbjvLiwh0sGN2hyO7ustkMnvxuG7/ujsXdxYmPB7cirLJXkbz2FTm9395H6NwJCGgAQxaCb7DZVYmISDllahgKCAggIKCUdNI1U422gAVO74NzMbk++B/v3oBF206y40Qis/88wj1tr7zPjM1m8PS8v/h+8zGcnSy8e3cLWtSodMWve8Xi9sGsm+HcSQhsaA9CPkFmVyUiIuVYmekzdOTIESIjIzly5AjZ2dlERkYSGRlJUlKSY5+GDRsyb948x3J8fDyRkZHs3LkTgKioKCIjI4mOLj0zwgPgWQmCm9qfH8ndOlTFx51Hu9UHYOKPuzgYl3xFp8vIsjH2m0hm/3kUJwu8fdfV9GgSckWvWSTi9tr7CJ07CYGNFIRERKRElJkw9Pzzz9OiRQteeOEFkpKSaNGiBS1atMjRpygqKoqEhATH8oIFC2jRogU33XQTAP3796dFixZMmzatxOu/rAv9hg7l7jcEMKR9OG1rVSYlI5sHP99EUnpWoU6TmJbJfbP+5IfIEzg7WXjrzqvpfVUp6JR8ao89CCVFQ1BjBSERESkxFqOkhzcuYxITE7FarSQkJODn51d8J9r9I8weAJXCYUykfXTqfzmZkEqf9/7g1Ll0rq1bhemDW+Pplv/JU7cfT2DUV5s5fDoFT1dnPhh4DV0alILAcSrK3kcoORaCmsCQBeCty6ciIlJ4Bfn8LjMtQ+VerU7g5GofYPD0/ovuUtXqyfTBrfB2c+aPfae5Z/o6Ys+lXXTff0rLzGbKsj3c+sEaDp9Oobq/J3MeaFc6glDs7r+DUHBTe4uQgpCIiJQghaHSwt0Hwq+1P9+79JK7XRXmz8xhbfDzcGHzkbN0n/IbX6w7TFpmdq59E1Iy+WT1Qbq+uYp3lu8lI9tGt0bB/DimA81D/YvpjRRA7C57Z+nkWAhudj4IVTG7KhERqWB0mewySuwyGcDa9+Hnp6F2F/ts7HnYfyqJ0V9tYddJ+wjZvh4utK1Vmer+nmRk29h58hy7TiSSkW0DoJrVg2dvbkyvpiFYLnIJrsTF7LTPPp8SByHNYPAC8KpsdlUiIlJOFOTzW2HoMko0DMXthfdagbMbPHnQ3lqUh8xsG1+us89ddiLh4pfLGob4MrBdTW69pjpebqVkKrqYHeeD0GkIaW6fa0xBSEREilBBPr9LyaejAPa5tyqF2/sNHfodGvTKc3dXZyeGXluLwRHhRB47y7ajZzmVlI6zkxP1gnxoVt1KzSpepaMl6ILo7fBZH3sQqnoVDPpBQUhEREylMFSaWCxQ9wb482OI+umyYegCJycL19SoxDWlYdDEvET/BbP6QGq8fT62QfPsYyyJiIiYSB2oS5uG9jGR2LUQsjPNraUo5QhC19hbhBSERESkFFAYKm1qdQTvIHtoOLDK7GqKxslt9j5CqfFQveX5FiF/s6sSEREBFIZKHydnaNLP/nz7d6aWUiRObj0fhM5A9VYKQiIiUuooDJVGTW+zf921CDIvP6hiqXUi0n5pLO0shLaGQXPBw2p2VSIiIjkoDJVGoW3ALxQyzsG+ZWZXUzgnttjvGks7a38/AxWERESkdFIYKo2cnKDprfbnkV+bW0thHN8En/WFtAQIawsDvwePYh6jSUREpJAUhkqrFgPtX/f8BGePmltLQRz64/ylMQUhEREpGxSGSqvABhB+HRg22DTT7GryZ98v8MVtkJFkvytu4Fxw9zW7KhERkTwpDJVmbYbbv26eBVnp5tZyObsWwlf9ISsV6vWAAd9edjoRERGR0kBhqDRrcBP4VoXkU7BzgdnVXNrWOfDNELBlQuN+cNcX4OphdlUiIiL5ojBUmjm7QKth9ud/vAM2m7n1XMz6D2HeA2Bkw9X3wO2fgoub2VWJiIjkm8JQadf6fnDzhZi/YPdCs6v5m80Gy16An54EDGjzAPR5zz5opIiISBmiMFTaeVWGiJH25ysmlY7WoawM+OFB+ONt+/L1z0Gv1+xDAoiIiJQx+vQqC9qNBHcrnNoFO+aaW0taInx1J2ybAxZn6PsBdHwcLBZz6xIRESkkhaGywNMf2o+2P1/2PKSfM6eO0/thejc4sAJcvWDAHGhxjzm1iIiIFBGFobIiYjT414TE4/DrhJI//77l8HEXiIuy3+E29Eeod0PJ1yEiIlLEFIbKCjcvuHmK/fn6D+HYppI5r2HA2vfhy9vto0qHtoYRK6H6NSVzfhERkWKmMFSW1O0Kze4EDPh+GKTEF+/5Us/AN4Ph56ftI2FfPdDeIuQbUrznFRERKUEKQ2VNr9fsl8vOHILv7wdbdvGc5/BamNoBdi0AJxfo+Sr0fQ9c3IvnfCIiIiZRGCprvCpD/y/BxRP2L7eP81OUt9unJ8HPz8DMGyHxGFSuDfctg3YP6Y4xEREplxSGyqKQZtDvfcACf06HhQ8XTQtR1E/wfltY+579sthVd8MDv6l/kIiIlGsuZhcghdT0NsjOhB8egi1fQOIJ+5g/flUL/lrHNsHy8XDwN/uyfw246S3dLSYiIhWCwlBZdlV/cPWEuSNg/68wNQK6PAMtBtrX5yUrA/b+DOumweHV9nVOrvbRrjv91373moiISAVgMQzDMLuI0iwxMRGr1UpCQgJ+fn5ml3Nxp/bA3OFwMtK+7BUAjftC7U5QuQ54VYGsVEiOg+htcGQ97PkZ0hPs+1uc7cGq83/trUIiIiJlXEE+vxWGLqNMhCGwXzLbNBP+eBcSjuTvGJ8QuPpu+2Sw1tBiLU9ERKQkFeTzW5fJygtnV2gzHFoOhX2/wP4VcGQNJJ6ElNP26TM8K0FQQ6jWAup2g9A2mlxVREQqPIWh8sbZFRr0sj8uMAzdFi8iInIJahaoCBSERERELklhSERERCo0hSERERGp0BSGREREpEJTGBIREZEKTWFIREREKjSFIREREanQFIZERESkQlMYEhERkQpNYUhEREQqNIUhERERqdAUhkRERKRCUxgSERGRCk1hSERERCo0F7MLKO0MwwAgMTHR5EpEREQkvy58bl/4HM+LwtBlnDt3DoCwsDCTKxEREZGCOnfuHFarNc99LEZ+IlMFZrPZOHHiBL6+vlgsFrPLKZUSExMJCwvj6NGj+Pn5mV1OhaefR+min0fpop9H6VKcPw/DMDh37hzVqlXDySnvXkFqGboMJycnQkNDzS6jTPDz89Mfl1JEP4/SRT+P0kU/j9KluH4el2sRukAdqEVERKRCUxgSERGRCk1hSK6Yu7s7L7zwAu7u7maXIujnUdro51G66OdRupSWn4c6UIuIiEiFppYhERERqdAUhkRERKRCUxgSERGRCk1hSERERCo0hSEptN9++43evXtTrVo1LBYLP/zwg9klVViTJk2idevW+Pr6EhQURL9+/YiKijK7rApt6tSpNG/e3DGYXEREBD/99JPZZcl5kyZNwmKx8J///MfsUiqk8ePHY7FYcjxCQkJMq0dhSAotOTmZq666ivfee8/sUiq8VatWMWrUKNatW8eyZcvIysqie/fuJCcnm11ahRUaGsqrr77Kxo0b2bhxI9dffz19+/Zlx44dZpdW4f3555989NFHNG/e3OxSKrQmTZpw8uRJx+Ovv/4yrRZNxyGF1qtXL3r16mV2GQIsWbIkx/KMGTMICgpi06ZNdOzY0aSqKrbevXvnWJ44cSJTp05l3bp1NGnSxKSqJCkpiXvuuYePP/6YCRMmmF1Ohebi4mJqa9A/qWVIpBxKSEgAoHLlyiZXIgDZ2dnMnj2b5ORkIiIizC6nQhs1ahQ33XQT3bp1M7uUCm/v3r1Uq1aNWrVq0b9/fw4cOGBaLWoZEilnDMNg7NixdOjQgaZNm5pdToX2119/ERERQVpaGj4+PsybN4/GjRubXVaFNXv2bDZv3syff/5pdikVXtu2bfnss8+oX78+MTExTJgwgfbt27Njxw6qVKlS4vUoDImUM6NHj2bbtm2sXr3a7FIqvAYNGhAZGcnZs2f5/vvvGTJkCKtWrVIgMsHRo0d55JFHWLp0KR4eHmaXU+H9s4tFs2bNiIiIoE6dOsyaNYuxY8eWeD0KQyLlyMMPP8yCBQv47bffCA0NNbucCs/NzY26desC0KpVK/7880/eeecdPvzwQ5Mrq3g2bdpEbGwsLVu2dKzLzs7mt99+47333iM9PR1nZ2cTK6zYvL29adasGXv37jXl/ApDIuWAYRg8/PDDzJs3j5UrV1KrVi2zS5KLMAyD9PR0s8uokLp27ZrrbqV7772Xhg0b8tRTTykImSw9PZ1du3Zx3XXXmXJ+hSEptKSkJPbt2+dYPnjwIJGRkVSuXJkaNWqYWFnFM2rUKL766ivmz5+Pr68v0dHRAFitVjw9PU2urmJ6+umn6dWrF2FhYZw7d47Zs2ezcuXKXHf+Scnw9fXN1YfO29ubKlWqqG+dCR5//HF69+5NjRo1iI2NZcKECSQmJjJkyBBT6lEYkkLbuHEjXbp0cSxfuM47ZMgQZs6caVJVFdPUqVMB6Ny5c471M2bMYOjQoSVfkBATE8OgQYM4efIkVquV5s2bs2TJEm644QazSxMx3bFjx7j77ruJi4sjMDCQdu3asW7dOmrWrGlKPRbDMAxTziwiIiJSCmicIREREanQFIZERESkQlMYEhERkQpNYUhEREQqNIUhERERqdAUhkRERKRCUxgSERGRCk1hSERERCo0hSERkXIuJiaGrl27EhoaqkliRS5CI1CLiJRzTz31FJmZmdx2223cdttt7N27F19fX7PLEik11DIkIiWic+fO/Oc//zHl3EOHDqVfv36X3D5z5kz8/f1LrJ7LCQ8P5+233y7wcadPnyYoKIhDhw7lWG+1WqlWrRr169cnICAADw+PHNtvv/123nrrrSuoWKRsUxgSqWAuFQxWrlyJxWLh7Nmzlzw2r9BgsVj44YcfiqTGiqKoQ9ikSZPo3bs34eHhOdaPHTuWV199lZCQEO69915cXV1zbH/++eeZOHEiiYmJRVaLSFmiMCQiUg6kpqbyySefcP/99+fatmfPHlJSUhg9ejR//PFHru3NmzcnPDycL7/8siRKFSl1FIZEpMRkZWUxevRo/P39qVKlCs8++yz/7LZ45swZBg8eTKVKlfDy8qJXr17s3bvXsf1CS8rPP/9Mo0aN8PHxoWfPnpw8edKxT3Z2NmPHjnWc48knn6QwXSMXLlxIy5Yt8fDwoHbt2rz44otkZWU5tlssFqZPn84tt9yCl5cX9erVY8GCBTleY8GCBdSrVw9PT0+6dOnCrFmzHK1vK1eu5N577yUhIQGLxYLFYmH8+PGOY1NSUhg2bBi+vr7UqFGDjz76KM96f/rpJ1xcXIiIiMi17dNPP6Vv37489NBDLFq0iFOnTuXap0+fPnz99dcF/C6JlA8KQyJSYmbNmoWLiwvr16/n3XffZcqUKUyfPt2xfejQoWzcuJEFCxawdu1aDMPgxhtvJDMz07FPSkoKb7zxBp9//jm//fYbR44c4fHHH3dsf/PNN/n000/55JNPWL16NfHx8cybN69Adf78888MHDiQMWPGsHPnTj788ENmzpzJxIkTc+z34osvcuedd7Jt2zZuvPFG7rnnHuLj4wE4dOgQt99+O/369SMyMpIHHniAZ555xnFs+/btefvtt/Hz8+PkyZOcPHky1/to1aoVW7ZsYeTIkTz00EPs3r37kjX/9ttvtGrVKtf6zMxMvvzySwYOHEjDhg1p3rz5RVuA2rRpw4YNG0hPTy/Q90qkXDBEpEIZMmSI4ezsbHh7e+d4eHh4GIBx5syZSx47Y8YMA8h1rLe3twEY8+bNu+SxnTp1Mho1amTYbDbHuqeeespo1KiRYRiGsWfPHgMw/vjjD8f2uLg4w9PT0/jmm29ynH/fvn2Ofd5//30jODjYsVy1alXj1VdfdSxnZmYaoaGhRt++ffN8X1ar1bF83XXXGa+88kqOfT7//HOjatWqjmXAePbZZx3LSUlJhsViMX766SfHe2vatGmO13jmmWdyfI//fd4LatasaQwcONCxbLPZjKCgIGPq1KmXfA99+/Y1hg0blmv9999/bwQGBhqZmZmGYRjGlClTjObNm+fab+vWrQZgHDp06JLnECmvXMyLYSJili5dujB16tQc69avX8/AgQMdyz4+Po7nAwcOZNq0aQD4+vqyefPmXK9Zr169y563Xbt2WCwWx3JERARvvvkm2dnZ7Nq1CxcXF9q2bevYXqVKFRo0aMCuXbsc67y8vKhTp45juWrVqsTGxgKQkJDAyZMnc1wqcnFxoVWrVgW6VLZp0yb+/PPPHC1B2dnZpKWlkZKSgpeXF2Dva3OBt7c3vr6+jlqioqJo3bp1jtdt06ZNvmv452tbLBZCQkIcr30xqampue4SA5gxYwZ33XUXLi72P/d33303jz/+OJs3b+aaa65x7Ofp6QnYW95EKhqFIZEKyNvbm7p16+ZYd+zYsRzLkZGRjud+fn6O505OTrmOLQqXCiuGYeQIUP++E8pisRSqT1BebDYbL774Irfeemuubf8MHBerxWazAbnrvrAuv/J67YsJCAjgzJkzOdZFR0ezZMkSbDZbjvCbnZ3NjBkzcoShC5f3AgMD812jSHmhMCQiF1UcgWfdunW5luvVq4ezszONGzcmKyuL9evX0759e8A+bs6ePXto1KhRvl7farVStWpV1q1bR8eOHQF7p+1Nmzbl+OC/nGuuuYaoqKgr+h40bNiQxYsX51i3cePGHMtubm5kZ2cX+hz/1KJFC7744osc6z7//HNq1qyZa8iD5cuX89JLL/HGG2/g7u4OwPbt2wkNDSUgIKBI6hEpS9SBWkRKzNGjRxk7dixRUVF8/fXX/O9//+ORRx4B7JfZ+vbty/Dhw1m9ejVbt25l4MCBVK9enb59++b7HI888givvvoq8+bNY/fu3YwcOTLPsZMu5vnnn+ezzz5j/Pjx7Nixg127djFnzhyeffbZfL/GAw88wO7du3nqqafYs2cP33zzDTNnzgRwtBiFh4eTlJTE8uXLiYuLu6JLVD169GDHjh05WodmzJjB7bffTtOmTXM8hg0bRnJyco67337//Xe6d+9e6POLlGUKQyJSYgYPHkxqaipt2rRh1KhRPPzww4wYMcKxfcaMGbRs2ZKbb76ZiIgIDMNg8eLFuS4Z5eWxxx5j8ODBDB06lIiICHx9fbnlllsKVGePHj1YtGgRy5Yto3Xr1rRr14633nqLmjVr5vs1atWqxXfffcfcuXNp3rw5U6dOddxNdqE1pn379jz44IPcddddBAYGMnny5ALV+U/NmjWjVatWfPPNN4C9D9iuXbu47bbbcu3r6+tL9+7d+fTTTwFIS0tj3rx5DB8+vNDnFynLNDeZiEgJmThxItOmTePo0aPF8vqLFy/m8ccfZ/v27Tg55f//uu+//z7z589n6dKlxVKXSGmnPkMiIsXkgw8+oHXr1lSpUoU//viD119/ndGjRxfb+W688Ub27t3L8ePHCQsLy/dxrq6u/O9//yu2ukRKO7UMiYgUk0cffZQ5c+YQHx9PjRo1GDRoEOPGjXPc5i4ipYPCkIiIiFRo6kAtIiIiFZrCkIiIiFRoCkMiIiJSoSkMiYiISIWmMCQiIiIVmsKQiIiIVGgKQyIiIlKhKQyJiIhIhfZ/TKWCrcFn89MAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plot the energies\n",
    "plt.plot(xyz[:, 1, 2], molDB.get_properties('HF_energy'), label='HF/STO-3G')\n",
    "plt.plot(xyz[:, 1, 2], molDB.get_properties('CISD_energy'), label='CISD/cc-pVDZ')\n",
    "plt.legend()\n",
    "plt.xlabel('H-H bond length (Å)')\n",
    "plt.ylabel('energy (hartree)')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# train ANI model with HF energies\n",
    "! rm ANI-HF.pt &> /dev/null\n",
    "ani = ml.models.ani(model_file='ANI-HF.pt', verbose=False)\n",
    "ani.train(molecular_database=molDB, property_to_learn='HF_energy')\n",
    "# predict with trained ANI model\n",
    "ani.predict(molecular_database=molDB, property_to_predict='ANI_HF_energy')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Sequential(\n",
       "  (0): AEVComputer()\n",
       "  (1): ANIModel(\n",
       "    (H): Sequential(\n",
       "      (0): Linear(in_features=48, out_features=160, bias=True)\n",
       "      (1): CELU(alpha=0.1)\n",
       "      (2): Linear(in_features=160, out_features=128, bias=True)\n",
       "      (3): CELU(alpha=0.1)\n",
       "      (4): Linear(in_features=128, out_features=96, bias=True)\n",
       "      (5): CELU(alpha=0.1)\n",
       "      (6): Linear(in_features=96, out_features=1, bias=True)\n",
       "    )\n",
       "  )\n",
       ")"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# show the model structure\n",
    "ani.model\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# fix some of the layers\n",
    "! cp ANI-HF.pt ANI-HF-TL.pt\n",
    "ani = ml.models.ani(model_file='ANI-HF-TL.pt', verbose=False)\n",
    "ani.fix_layers([[0, 6]])\n",
    "# transfer leaning with every 40th of the data  \n",
    "step = 40\n",
    "val = molDB[::step][::10]\n",
    "sub = ml.molecular_database([mol for mol in molDB[::step] if mol not in val])\n",
    "ani.energy_shifter.self_energies = None # let the model recalculate the self atomic energies\n",
    "ani.train(molecular_database=sub, validation_molecular_database=val, property_to_learn='CISD_energy', hyperparameters={'learning_rate': 0.0001}, reset_optimizer=True)\n",
    "# predict with TL model\n",
    "ani.predict(molecular_database=molDB, property_to_predict='ANI_TL_energy')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# train ANI model with CISD energies directly\n",
    "! rm ANI_CISD.pt &> /dev/null\n",
    "ani_cisd = ml.models.ani(model_file='ANI_CISD.pt', verbose=False)\n",
    "ani_cisd.train(molecular_database=sub, validation_molecular_database=val, property_to_learn='CISD_energy')\n",
    "# predict with trained ANI model\n",
    "ani_cisd.predict(molecular_database=molDB, property_to_predict='ANI_CISD_energy')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAGyCAYAAADnH8C6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAAClfElEQVR4nOzdd3QU5dfA8e9syaZvegECSSihF+lFinQUFBBFEEUUVAQUOyqIhSJi+dkbUhRfsSBNBamKVCmRFkILPYX0nmyZ94/ASkhIgSSbcj/n7El25pmZu0vYuftURVVVFSGEEEKIGkpj7wCEEEIIIexJkiEhhBBC1GiSDAkhhBCiRpNkSAghhBA1miRDQgghhKjRJBkSQgghRI0myZAQQgghajSdvQOo7KxWKxcvXsTNzQ1FUewdjhBCCCFKQFVV0tLSqFWrFhpN0XU/kgwV4+LFiwQFBdk7DCGEEELcgHPnzlGnTp0iy0gyVAw3Nzcg7810d3e3czRCCCGEKInU1FSCgoJs9/GiSDJUjCtNY+7u7pIMCSGEEFVMSbq4SAdqIYQQQtRokgwJIYQQokaTZEgIIYQQNZr0GSojFosFk8lk7zBEFefg4FDsEFAhhBBlS5Khm6SqKjExMSQnJ9s7FFENaDQaQkJCcHBwsHcoQghRY0gydJOuJEJ+fn44OzvLxIzihl2Z4DM6Opq6devK35IQQlQQSYZugsVisSVC3t7e9g5HVAO+vr5cvHgRs9mMXq+3dzhCCFEjSOeEm3Clj5Czs7OdIxHVxZXmMYvFYudIhBCi5pBkqAxIc4YoK/K3JIQQFU+SISGEEELUaJIMCSGEEKJGqzLJUFJSEmPGjMFoNGI0GhkzZkyJhrNHREQwZMgQjEYjbm5udOrUibNnz5Z/wEIIIYSoEqpMMjRq1CjCw8NZu3Yta9euJTw8nDFjxhR5zMmTJ+nWrRuNGzdmy5Yt/Pvvv0yfPh1HR8cKirryGjt2LHfddVeB7Vu2bEFRFJKTk22/X/t45ZVXChwTGBiIqqrExcXx6KOPUrduXQwGAwEBAfTv358dO3Zc93xXPxYtWgTkdSB+7733aNmyJY6Ojnh4eDBw4EC2bdtW7GtLSEhgwIAB1KpVC4PBQFBQEJMmTSI1NTVfOVVV+fLLL+ncuTPu7u64urrSrFkznnzySU6cOHHD760QQoiSMVmsxKVlkpiZadc4qsTQ+oiICNauXcvOnTvp2LEjgO0mFhkZSVhYWKHHvfzyywwaNIh58+bZtoWGhhZ5rZycHHJycmzPr72B1kSRkZG4u7vbnru6uubbv2rVKoYMGYKiKAwfPhyTycTixYsJDQ0lNjaWjRs3kpiYSN++fYmOjrYd9+STT5KamsrChQtt24xGI6qqMnLkSDZs2MDbb79N7969SU1N5eOPP6Znz578+OOPhSZyV2g0Gu68807efPNNfH19OXHiBE888QSJiYl89913QF4iNGrUKFasWMFLL73Ee++9h5+fH1FRUaxfv54333zTlpgJIYQoXK7ZwsW0ZC6mJhKTnkBcRjIJmckkZaWQnJtKWm4q6aZUssxpZFszyLWmYyYTC7mo5IJiQtFYaO06gm+Gz7Db66gSydCOHTswGo22RAigU6dOGI1Gtm/fXmgyZLVa+fXXX3n++efp378/+/fvJyQkhGnTphV5I50zZw6vvfbaDceqqipZpoofFu2k15bbSCQ/Pz88PDyuu3/VqlX873//Izk5mb///pstW7bQo0cPAOrVq0eHDh1sZQMCAv6L2cmJnJycfNsAli1bxk8//cSqVasYPHiwbfsXX3xBQkICjzzyCH379sXFxaXQeDw9PXn88cdtz+vVq8fEiRN5++23813j+++/Z+XKlQwZMsS2PTQ0lN69e6OqajHvihBCVA85JgsxaSlcSEsgJi2JSxlJxGddTmhyUkjJTSXDlEamOY1sSzq5agZmMrAqmaDJQlFK+HmpANr8T20xWLLL8iWVWpVIhmJiYvDz8yuw3c/Pj5iYmEKPiYuLIz09nblz5/Lmm2/y1ltvsXbtWoYNG8bmzZttN+trTZs2jaefftr2PDU1laCgoBLHmmWy0HTGuhKXLytHXu+Ps0PF/3MePnyYmJgYevfujU6nw9XVlRUrVtCpUycMBsMNnfO7776jUaNG+RKhK5555hmWL1/O+vXri0xqr3bx4kWWL1+e79/8//7v/wgLC8uXCF1NhrgLIaqaHHMOMelJnE+J52JaIjFpiVzKTCIhK5nknBRSLyc1WZZ0cqzp1yQ01pJd5DrJDFYdiuqCDmf0igsGjSvOOldcdO646t0xGtzxcDTi6WjEx9mIl5MrXs6ueDu74OXsgrvB9dorVSi7JkMzZ84sthbmn3/+AQq/Oamqet2bltWa9w975513MnXqVABat27N9u3b+eyzz66bDBkMhhu+iVc1a9asKdDkVdhkf3Xq1Mn3/MyZM7YZt1euXEn//v1t/bAWLVrE+PHj+eyzz7jlllvo0aMHI0eOpGXLliWO69ixYzRp0qTQfVe2Hzt2rNjz3HfffaxcuZKsrCwGDx7MV199le8a19YoPvXUU7YyHh4enD9/vsQxCyFEWTJZTMRkxBOVFMOZ5FjOp8YSl5FAUnYKKbkpVzU9pWNS0y8nNbklO7kGW4/hfHdQVYtidUGLMw5XJzR6d9z07ng4uuclM04e+Ll44u/qRaCbJ7WN3jjrq3ZfXLsmQ5MmTWLkyJFFlgkODubAgQPExsYW2Hfp0iX8/f0LPc7HxwedTkfTpk3zbW/SpAl///33jQddDCe9liOv9y+38xd13dLq1asXn376ab5tu3bt4v7778+3bevWrbi5udmee3p62n5fuXIlEydOtD0fPnw4t99+O1u3bmXHjh2sXbuWefPm8dVXXzF27NhSx3g9V5LggQMHsnXrViCvOezw4cO2Mu+99x6vvvoqkZGRvPTSSzz99NN88sknBc5xxcsvv8ykSZNYvnw5s2fPLrNYhRACwGK1cD41nhMJFzmdHM3FtDhiMi6RkJVAcm4i6aYksqwpmEhB1ZSwQ/E1w6BUVQGrExrVGT2uGDSuOGndcNG74e6QVzvj7WTEx9kTfxdP/F09qeXuRaCbFy56pxpbK27XZMjHxwcfH59iy3Xu3JmUlBR2795t63+ya9cuUlJS6NKlS6HHODg40L59eyIjI/NtP3bsGPXq1bv54K9DURS7NFfdCBcXFxo0aJBvW2G1ISEhIYX2GYqJiWHfvn3cfvvt+bY7OjrSt29f+vbty4wZM3jkkUd49dVXS5wMNWrUiCNHjhS6LyIiAoCGDRsC8NVXX5GVlQVQYC2vgIAAAgICaNy4Md7e3tx6661Mnz6dwMBAGjZsyNGjR/OV9/X1xdfXt9AmWSGEuJ5cs5lTibEcS7jA6aQYzqfFEJNxicTsS6TkJpJpTcJEMlZNavHNUVclN6qqAYsLGqs7BsWIk9aIq94dN70Ro8GIt5MHPs6eBLh6EeDqSZCHD7XdPXExyLqGpVUl7tpNmjRhwIABjB8/ns8//xyACRMmcMcdd+Rr6mjcuDFz5sxh6NChADz33HPce++9dO/enV69erF27VpWr17Nli1b7PEyqp1Vq1bRuXPnYhPapk2bsmLFihKfd+TIkYwaNYrVq1cX6Df0zjvv4O3tTd++fQGoXbt2ic55pUP0lZGC9913H6NGjWLlypXceeedJY5NCFFzmMxmTiTGcDz+IlHJ0ZxPiyE24xIJ2ZdIzU0k05pIrpqCqk27fpJzVafhK3UuqtkVnWrEoPHAReuJ0cETL4MXvs4+BLr5UcfdjxDPAII9ffBwMqDR1MzamopUJZIhgKVLlzJlyhT69esHwJAhQ/joo4/ylYmMjCQlJcX2fOjQoXz22WfMmTOHKVOmEBYWxs8//0y3bt0qNPbqatWqVfkSiYSEBEaMGMG4ceNo2bIlbm5u7Nmzh3nz5pUq4Rg5ciQ//vgjDz74YIGh9atWreLHH3+87kgygN9++43Y2Fjat2+Pq6srR44c4fnnn6dr164EBwfbrrF8+XJGjhzJtGnT6N+/P/7+/pw5c4Zly5ah1Za+2VEIUTWoqkpydjLHEs9yLP4CZ1KiOZcay6XMSyTlxJNuTiJHTcKqSbv+SKlrkhxVVVCsrjhgxFnrhbveC29HX/xcfKnj5k9dYwChXrVo4BWAq6Fq96+pjqpMMuTl5cW3335bZJnChkOPGzeOcePGlVdYNVZGRgYbN27kvffes21zdXWlY8eOvPfee5w8eRKTyURQUBDjx4/npZdeKvG5FUXhhx9+4H//+x/vvfceTzzxBAaDgc6dO7N58+Zik1knJye+/PJLpk6dSk5ODkFBQQwbNowXX3wx3zWWLVvGl19+ycKFC5k3bx4mk4k6derQu3dv3n333dK/KUKISsNqtXL00gXCY05wNP40USlniM28QLI5mmw1DlWTdf2Dr+pc/F+S45EvyfF3yavBCfaoRQPvQBp4BWLQSfNUVaWoMqFKkVJTUzEajaSkpOSbeBAgOzubqKgoQkJCatys1suXL+eVV165bt8ecWNq8t+UEKWVlJlNeHQUB2NPcSLxNOfSzxGffZF0SwwmTTyKxlTk8VazG1qLEUeNJ256L7wMPgS4+lHHzZ9Qr1o08qlFI59AnPQOFfSKRFkq6v59rSpTMyQqF1dXV9566y17hyGEqOYSMjLYc/4U/8Ye51jiaS6knych5yJZahyqLhFFuWY6EAXQ/Vero7F44azxw8uhFoHOdQgx1iXMJ4QW/qEEe3ngeAMjcUX1I8mQuCFX+m4JIcTNSspKZ/f54xyIOcHxxDOcTz9LQs5FMtU4VG1ywX47l+9cCoCqRa/64qYNwM+pNnXdgmjkFUIL/xBaB4bi7FAz5o0TN0eSISGEEOXOYrUQlXyOHeciOBB7jBPJp4jNOkeGJRartpA1IK8egWU1YMAXD30gtVzqUN8jmGZ+IbQJbECwRy20GqndETdHkiEhhBBlKiU7nY2n9rHz/AGOJh0hNjuKTDUGFHPBwlfyGIszBvxsCU+oRz1a+NWnXZ0G1DX619jJAEXFkGRICCHEDYtNTWX9yf3svPAvx5MjiMs9iUkTU7BpSwHVqkMx++KmqU2gcz0aeIbQ0r8BnYIaEurtKwmPsBtJhoQQQpTIyfhLrDu+j70xhziVeoxEUxQW3TWJj/Zy05bZiKsSTG3nhjTxakK7Wo3pULc+AW7OkvSISkeSISGEEAWcTLzIuuN7+Cf6IKdSjpFkPo2qS8xfSJ+X+CgWdzy0oQS7hdHGvwW3hdxCy8A6kvSIKkOSISGEqOHScrL4/dge/jqzl4ikg8SbjmPVJuUvdPluobV44+0QTH1jGG0DmtOnflvqe9Wq+KCFKEOSDAkhRA2iqipH48/w67Gd7L64n6j0I2Qp5/LP16PNm6NHa/HHVx9KQ8/GdKrdgn4NbyHQzct+wQtRTjTFFxHVVUxMDJMnTyY0NBSDwUBQUBCDBw9m48aNAAQHB/P+++/byu/fv5877rgDPz8/HB0dCQ4O5t577yU+Ph6A06dPoyiK7eHm5kazZs144oknOH78eKExLFq0iE6dOpX7ay1MixYteOSRRwrd93//93/o9XpiY2PZsmWL7TVpNBqMRiNt2rTh+eefJzo62nbMta//2kdISEhFvTQhbFRV5Z8LR3lpw+f0WzqO1ou6cs9vg1l8YhYRmb+RrTmdlwhZXPFQW9PefTSTm8xn7V1/8u/D69nwwOd8OngqD97SRxIhUW1JzVANdfr0abp27YqHhwfz5s2jZcuWmEwm1q1bxxNPPMHRo0fzlY+Li6NPnz4MHjyYdevW4eHhQVRUFKtWrSIzMzNf2Q0bNtCsWTMyMzM5ePAg//vf/2jVqhWrV6+md+/e+cpeu9hrRXr44YeZMWMGH3zwAc7Ozvn2ff3119xxxx34+/sTEREB5C0E7O7uTmpqKvv27WPevHksWLCALVu20KJFC4KCgvIlR1fs2bOHu+66iyeeeKJCXpcQkZcu8MPhTey4uJML2QewapP/26kBVdWiN9ehlmMYt/i3YkDDjnQKaohWK9+PRQ2liiKlpKSogJqSklJgX1ZWlnrkyBE1KyvLDpHdnIEDB6q1a9dW09PTC+xLSkpSVVVV69Wrp7733nuqqqrqL7/8oup0OtVkMl33nFFRUSqg7t+/P992i8Wi9uzZU61Xr55qNptt27OyslQXFxf10KFDqqqqanZ2tvrcc8+pderUUR0cHNQGDRqoX331la38oUOH1EGDBqlubm6qq6ur2q1bN/XEiRPXjadHjx7qE088oT7xxBOq0WhUvby81Jdfflm1Wq2qqqpqfHy86uDgoC5atCjfcWfOnFE1Go26evVqVVVVdfPmzSpge1+uyMzMVMPCwtSuXbteN4aYmBi1Tp066ujRo69b5mpV+W9K2M/55AT13b+Xq3f937Nq6wV91OaLmud7NPu6ldruq2HqyGVvqJ/uWK+eS0q2d8hClLui7t/XkpqhsqaqYMosvlxZ0ztDCUduJCYmsnbtWmbNmoWLi0uB/R4eHgW2BQQEYDab+eWXX7j77rtLNUpEo9Hw5JNPMnToUPbu3UuHDh0A2LhxIwEBATRr1gyABx54gB07dvDBBx/QqlUroqKibE1wFy5coHv37vTs2ZNNmzbh7u7Otm3bMJsLmcTtKosXL+bhhx9m165d7NmzhwkTJlCvXj3Gjx+Pt7c3d955JwsXLuTBBx+0HbNw4UL8/f0ZOHBgked2cnLiscceY+rUqcTFxeHn55dvv8lkYvjw4QQEBPDVV1+V+P0SojhZpmx+idjObyf+4mjyvstNXZeHt1/u7+NgCSLUtQ231evC3c274eda9EKVQtRkkgyVNVMmzLbDyIqXLoJDwcSmMCdOnEBVVRo3blzi03fq1ImXXnqJUaNG8dhjj9GhQwduu+02HnjgAfz9/Ys9/sq1Tp8+bUuGVq5caWsiO3bsGD/88APr16+nT58+AISGhtqO//jjjzEajXz//ffo9XoAGjVqVOx1g4KCeO+991AUhbCwMA4ePMh7773H+PHjARg3bhyDBg3i1KlThIaGoqoqixYtYuzYsWi1xU/xf/XrujYZmjRpEidOnGDPnj2yAr24aYdiT7H43zXsjtlBouUoaC5/Ebg8r4/W7EeQU0u61enCPc17EOLlV+T5hBD/kQbiGkhV875BlnYOkFmzZhETE8Nnn31G06ZN+eyzz2jcuDEHDx4s9TVVVWX16tUMGTIEgPDwcLRaLT169Cj0+PDwcG699VZbInS1rVu34urqanssXbrUtq9Tp075Xmfnzp05fvw4FkveyJl+/fpRp04dFi5cCMCmTZs4ffo0Dz30UEnekuu+l5999hmLFi1i+fLl1KlTp0TnEuJqVquV34/t4YHlr9P26wHct/ZO1kZ/SaJ6KC8RMrvhq3TmjsCn+KbfasIf3sjqUe/xQvcRkggJUUpSM1TW9M55tTT2uG4JNWzYEEVRiIiI4K677irVZby9vRkxYgQjRoxgzpw5tGnThvnz57N48eIij7vSCfnKiKrdu3eTm5tLt27dgLwmp6IUtb9du3aEh4fbnpekpuoKjUbD2LFjWbRoEa+99hoLFy6ke/fuNGzYsETHX3ldwcHBtm1///03U6ZM4ZNPPqFLly4ljkWIHLOJZQf/YsWxdZzM2PnfXD9aUFUNjuaGtPDqzJCGPbm9SWscdLJAqRBlQZKhsqYoJW6ushcvLy/69+/Pxx9/zJQpUwr0G0pOTi6039C1HBwcqF+/PhkZGUWWs1qtfPDBB4SEhNCmTRsgr4ns9ttvtzVFtWjRAqvVyp9//mlrJrtay5YtWbx4MSaTqUDtkJOTEw0aNCj02jt37izwvGHDhvmawB566CHefPNNli9fzvLly/nss8+Kfe0AWVlZfPHFF3Tv3h1fX18Azp07x/Dhw5kwYcJ1h+0LcbXU7CwW7f+D30+t53zOHtBe/v+kBdWqx0vTgltr9WJsm4E09Cl5oi+EKDlJhmqoK7UWHTp04PXXX6dly5aYzWbWr1/Pp59+aqvxuGLNmjV8//33jBw5kkaNGtmauX777TdbE9MVCQkJxMTEkJmZyaFDh3j//ffZvXs3v/76qy0JWbVqFa+99prtmODgYB588EHGjRtn60B95swZ4uLiuOeee5g0aRIffvghI0eOZNq0aRiNRnbu3EmHDh0ICwu77us8d+4cTz/9NI8++ij79u3jww8/5J133slXJiQkhNtuu40JEyag1+u5++67Cz1XXFwc2dnZpKWlsXfvXubNm0d8fDzLly8HIDs7m6FDh1K7dm1efPFFYmJiCpwjICCgiH8VUVOkZKezYO/v/HbqD2LM+1E0OXk7tIDFmVoObekb3JuH2vTD28XNrrEKURNIMlRDhYSEsG/fPmbNmsUzzzxDdHQ0vr6+tG3blk8//bRA+aZNm+Ls7MwzzzzDuXPnMBgMNGzYkK+++ooxY8bkK3ulZsfZ2Zl69erRq1cvvvjiC1vtzcmTJzlx4gT9+/fPd9ynn37KSy+9xMSJE0lISKBu3bq89NJLQF7z3KZNm3juuefo0aMHWq2W1q1b07Vr1yJf5wMPPEBWVhYdOnRAq9UyefJkJkyYUKDcww8/zMaNG5kwYUKBOYeuCAsLQ1EUXF1dCQ0NpV+/fjz99NO2BGfXrl3s3bsXyOu4XZgrfYxEzZOUlcLX+1fze9QGYk3/gpLXAVrRgGIxEuLUkbsaDmBkqx446R3sHK0QNYuiyqdzkVJTUzEajaSkpODunn9oanZ2NlFRUYSEhMhooVJ499132bBhA7/99lu5Xqdnz560bt063yzalZ38TVUvOeYclvy7lp8iV3Ixdx9cveSFyYdGrp25p+lAhjfrgq4EoxeFECVX1P37WlIzJCpcnTp1mDZtmr3DEKJcqKrKltP/8Nm+ZUSkbkXVZOXtUIDcABq6dOHeZoMY1qwdeukALUSlIMmQqHD33HOPvUMQoswlZKby3o7/449zK8hSzudt1IBqNlLP0I2RTe/k3pYdcdDJjCZCVDaSDIlqa8uWLfYOQdQA647v5+M9S4jK3gqaHFBAteowWtszuP5gHuvQFw9nafIUojKTZEgIIUopNSeTd7f9xK+nl5OtPZm3UQOKyY92XoOY2nkULQID7RukEKLEJBkSQogS2nX2GO/sXExE+sa8+YAuT4boq2nL6CYjGXtLH3Sy8rsQVY4kQ0IIUQSr1cpXe9ez5PA3JHMgb0FULShmD27xGsjzXR+gqZ8suSJEVSbJkBBCFCI5K4M5f33HH+d/xKyLBiVvQJi72pwRje7l8fa3YyhkrTwhRNUjyZAQQlzlcOwZZm9dwIHUP/KawnSgWh1o4NST5zs/Qpd6TewdohCijEkyJIQQwIaT+5m341MumnehKNbLTWFedPO7k5d7jKW2u5e9QxRClBNJhoQQNZaqqvx85C8+2PsFSeoBIG+tZSdLI+5uOJInO90lTWFC1AAy7KGG2759O1qtlgEDBuTbfvr0aRRFwc/Pj7S0tHz7WrduzcyZM23Pe/bsyVNPPVXkdYKDgwtdFmPmzJm0bt0633NFUQo8NmzYUNqXJsR1WVUrX+5dRefFQ3ltzySS1AOoqoI3HXmt3ZfsHvczz986QhIhIWoIqRmq4b7++msmT57MV199xdmzZ6lbt26+/WlpacyfPz/fCvPlrVmzZgWSHy8vaaIQN89itfDp7hUsiviCHM1F2wSJtXXdmdb1UXqGNrV3iEIIO5BkqAbLyMjghx9+4J9//iEmJoZFixYxY8aMfGUmT57Mu+++yxNPPIGfn1+FxKXT6WwrwQtRFlRVZcG+NXx+4BOyNefzlsmwGGjg2J8Z3R/lljp1iz+JEKLakmSojKmqSpY5q8Kv66RzQlGUUh2zbNkywsLCCAsL4/7772fy5MlMnz4933nuu+8+1q9fz+uvv85HH31U1mELUa5UVWXpgXV8uP9jMpXTtiSoifPtzOr9BI18KybBF0JUbpIMlbEscxYdv+tY4dfdNWoXznrnUh2zYMEC7r//fgAGDBhAeno6GzdupE+fPrYyiqIwd+5cBg8ezNSpU6lfv/4Nx/jCCy/wyiuv5NuWm5tL06b5myYOHjyIq6ur7XnTpk3ZvXv3DV9X1Ey/Ru5g9s63SeX45eYwBxoYBjB34GQa+0nNoxDiP5IM1VCRkZHs3r2b5cuXA3lNU/feey9ff/11vmQIoH///nTr1o3p06fz3XffFXne2bNnM3v2bNvzI0eO2PohPffcc4wdOzZf+Q8++IC//vor37awsDBWrVple24wGEr9+kTNte/iCZ7fNJdYyy4AVKueevo+zOozhda1ZaZoIURBkgyVMSedE7tG7bLLdUtjwYIFmM1mateubdumqip6vZ6kpKQC5efOnUvnzp157rnnijzvY489xj333GN7XqtWLdvvPj4+NGjQIF/5wjpGOzg4FCgnRHEupCbwzB/vcCj9NxTFgqoq+CvdeKPXM3QJvvEaTSFE9SfJUBlTFKXUzVUVzWw2s2TJEt555x369euXb9/w4cNZunQpd9xxR77tHTp0YNiwYbz44otFntvLy0tGfokKlZmbzcsbP2dD9HegzURRwMXSlJc6P8eQJu3sHZ4QogqQZKgGWrNmDUlJSTz88MMYjcZ8++6++24WLFhQIBkCmDVrFs2aNUOnkz8bUTl8sms1Xxx+F4s2HrSgNQfyUOPJTOlyR6kHFAghai6ZdLEGWrBgAX369CmQCEFezVB4eDiJiYkF9jVq1Ihx48aRnZ1dEWEKcV27z0XSffH9fHr0pbxEyOxGX99J7BjzK092HSyJkBCiVBRVVVV7B1ESSUlJTJkyxdaxdsiQIXz44Yd4eHhc95jrfSDOmzev2L4vV6SmpmI0GklJScHd3T3fvuzsbKKioggJCcHR0bFkL0SIIsjfVNGSstKZ8vt89qeuuNwvSEtjp9v5cMBzBBo97B2eEKISKer+fa0q094xatQozp8/z9q1awGYMGECY8aMYfXq1dc9Jjo6Ot/z33//nYcffpjhw4eXa6xCiLKlqirvbv+RJZEfYtUm5/ULsjZldvdXuK1+C3uHJ4So4qpEMhQREcHatWvZuXMnHTvmzeHz5Zdf0rlzZyIjIwkLCyv0uGtnMV65ciW9evUiNDS03GMWQpSN8ItRTF4/nWT+ta0kP6rBZJ6/dRgajbT0CyFuXpVIhnbs2IHRaLQlQgCdOnXCaDSyffv26yZDV4uNjeXXX39l8eLFRZbLyckhJyfH9jw1NfXGAxdC3LBcs5ln133CprjFKJpcVFVLC5ehfDjoGXxcXIs/gRBClFCVSIZiYmIKXRfLz8+PmJiYEp1j8eLFuLm5MWzYsCLLzZkzp0IXJRVCFPTH8X28tHUGOdozKBpwtNTnrR5vSJOYEKJc2LWOeebMmSiKUuRjz549QOGdoVVVLfGoka+//prRo0cX2yl12rRppKSk2B7nzp0r/QsTQtyQlOxMRv00g6e3PUSO9gyq1ZF+fhPZ8eDPkggJIcqNXWuGJk2axMiRI4ssExwczIEDB4iNjS2w79KlS/j7+xd7na1btxIZGcmyZcuKLWswGGT5ByHsYPnhHby+cwYWXQyKAp6045OBr9E8QFaUF0KUL7smQz4+Pvj4+BRbrnPnzqSkpLB79246dOgAwK5du0hJSaFLly7FHr9gwQLatm1Lq1atbjpmIUTZSs/J5tHVc/g3fQWKzgoWN8Y0nMpz3e6W+YKEEBWiSgzFaNKkCQMGDGD8+PHs3LmTnTt3Mn78eO644458nacbN27ML7/8ku/Y1NRUfvzxRx555JGKDlsIUYzfjv1D96V3ciBjOYpixU/TiTVDV/D8rSMkERJCVJgq0YEaYOnSpUyZMsW2ltaQIUP46KOP8pWJjIwkJSUl37bvv/8eVVW57777KixWIUTRsk25PP7rPP5J/hFFawWLK6MbTOXF7vcUf7AQQpSxKjMDtb3IDNSiItWEv6ltZw7z5KbnyNHkDU7woh1f3j6bRj6Bdo5MCFGdlGYG6irRTCbKz/bt29FqtQwYMCDf9tOnT6MoCn5+fqSlpeXb17p1a2bOnGl73rNnT5566qkC5y7JaMHTp08zc+ZMWrduXQ6vTlQmqqoyfdOXPLZpTF4iZHFmeJ0X2fLA15IICSHsSpKhGu7rr79m8uTJ/P3335w9e7bA/rS0NObPn39D53722WeJjo62PerUqcPrr7+eb1tQUNDNvgRRBUSnxdN36ThWnPsANCacLU34dsCPzOw9WvoGCSHsTpKhGiwjI4MffviBxx9/nDvuuINFixYVKDN58mTeffdd4uLiSn1+V1dXAgICbA+tVoubm1uBbaJ6++HQZgb8OJRYyx5UVUt7twfZ+sB3tKolQ+aFEJVDlelAXVWoqoqalVXh11WcnEr9DXvZsmWEhYURFhbG/fffz+TJk5k+fXq+89x3332sX7+e119/vUCHdSGKYrFamPjrPLYl/B+KVkUx+fFy+zcZ2aqzvUMTQoh8JBkqY2pWFpG3tK3w64bt24vi7FyqYxYsWMD9998PwIABA0hPT2fjxo306dPHVkZRFObOncvgwYOZOnUq9evXL9O4RfV0ITWe+1Y8SZJ6AEUBX7rxzYi51DYa7R2aEEIUIM1kNVRkZCS7d++2zQCu0+m49957+frrrwuU7d+/P926dWP69OnFnnf27Nm4urraHoX1QxLV22/HdjHop2EkqQdQrTr6+T7Jxgc+kURICFFpSc1QGVOcnAjbt9cu1y2NBQsWYDabqV27tm2bqqro9XqSkpIKlJ87dy6dO3fmueeeK/K8jz32GPfc899cMbVq1SpVXKLqUlWVVzZ9wcpzn6JoLShmH17v9BZDm3Wwd2hCCFEkSYbKmKIopW6uqmhms5klS5bwzjvv2CaxvGL48OEsXbqUO+64I9/2Dh06MGzYMF588cUiz+3l5YWXl1eZxywqt8zcHEYuf56onE0oCrhbbmHpsHcJ9vS2d2hCCFEsSYZqoDVr1pCUlMTDDz+M8Zqmi7vvvpsFCxYUSIYAZs2aRbNmzdDpyv7PJisri/Dw8HzbXF1dadCgQZlfS5St4/ExjF71OFnaE6iqQhvX0Xx917PodTJSUAhRNUifoRpowYIF9OnTp0AiBHk1Q+Hh4SQmJhbY16hRI8aNG0d2dnaZx3Ts2DHatGmT7yHryVV+vx7dy/CV95ClPQFWRx5p+Abf3P2CJEJCiCpFluMohizHISpSVfqbenvrTyw+MQdFk4vW7Mv7t31Az5Dm9g5LCCGA0i3HIc1kQohSsVqtPLzyLfakfoeiAVdrE34Y/glBHj72Dk0IIW6IJENCiBJLzcpm2I9TiVX/BqCBYQDfDX8TJ73BzpEJIcSNk2RICFEiJ+LjGLnycXJ0x1BVDYMCH2de/8fsHZYQQtw0SYaEEMXaGhXJpE1PYNXFgtXA063eYNwtA+0dlhBClAlJhoQQRVoavpW5+54HXToaiwcf9PqIHiGt7B2WEEKUGUmGhBDXNefPH1h6ag6K1ozBGsR3Q76gkU8de4clhBBlSpIhIUQBVquVx9e8x7bExSgaFU9a8su9n+LtXPTwVCGEqIokGRJC5JNrNnP3j88TlbseRYH6hn78cPdcHHR6e4cmhBDlQpIhIYRNWk4WQ5ZNJF7dg6oq9PB5iI9ufwpFUewdmhBClBtJhoQQAFzKSGHIjw+TrkSiWrWMDp3GtB732jssIYQod7I2WQ23fft2tFotAwYMyLf99OnTKIqCn58faWlp+fa1bt2amTNn2p737NmTp556qthr7d+/nxEjRuDv74+joyONGjVi/PjxHDt2LN81r16w9eeff6Zjx44YjUbc3Nxo1qwZzzzzjG3/okWLUBQFRVHQarV4enrSsWNHXn/9dVJSUkr/htRQp5NjGPjDfZcTIQNPNp8niZAQosaQZKiG+/rrr5k8eTJ///03Z8+eLbA/LS2N+fPn3/R11qxZQ6dOncjJyWHp0qVERETwzTffYDQamT59eqHHbNiwgZEjR3L33Xeze/du9u7dy6xZs8jNzc1Xzt3dnejoaM6fP8/27duZMGECS5YsoXXr1ly8ePGmY6/uDsScZOjyUeRozoHFldfafcT49v3sHZYQQlQYaSarwTIyMvjhhx/4559/iImJYdGiRcyYMSNfmcmTJ/Puu+/yxBNP4Ofnd0PXyczM5KGHHmLQoEH88ssvtu0hISF07NiR5OTkQo9bs2YN3bp147nnnrNta9SoEXfddVe+coqiEBAQAEBgYCBNmjRh8ODBNGvWjOeff55vv/32huKuCbadPcTEDY9i1aaC2Yv3e3xC7wbN7B2WEEJUKKkZKmOqqmLKsVT4Q1XVUse6bNkywsLCCAsL4/7772fhwoUFznPffffRoEEDXn/99Rt+T9atW0d8fDzPP/98ofs9PDwK3R4QEMDhw4c5dOhQqa/p5+fH6NGjWbVqFRaLpdTH1wSbo8J5fMMjWLWpaEy1WNhvsSRCQogaSWqGypg518oXT/5Z4ded8L8e6A3aUh2zYMEC7r//fgAGDBhAeno6GzdupE+fPrYyiqIwd+5cBg8ezNSpU6lfv36pYzt+/DgAjRs3LtVxkydPZuvWrbRo0YJ69erRqVMn+vXrx+jRozEYil8YtHHjxqSlpZGQkHDDtVrV1fqTe3n6r4mgzURnqsv/DVlAY78Ae4clhBB2ITVDNVRkZCS7d+9m5MiRAOh0Ou69916+/vrrAmX79+9Pt27drtu352qzZ8/G1dXV9jh79uwN1VoBuLi48Ouvv3LixAleeeUVXF1deeaZZ+jQoQOZmZnFHn/lujIsPL/fj+3k6b8eA00melMwPw9bLImQEKJGk5qhMqZz0DDhfz3sct3SWLBgAWazmdq1a9u2qaqKXq8nKSmpQPm5c+fSuXPnfP13CvPYY49xzz332J7XqlWLRo0aAXD06FE6d+5cqjgB6tevT/369XnkkUd4+eWXadSoEcuWLeOhhx4q8riIiAjc3d3x9vYu9TWrq9+O7eKFbZNAk4ODqT7Lh31NPS8ve4clhBB2JclQGVMUpdTNVRXNbDazZMkS3nnnHfr1yz9qaPjw4SxdupQ77rgj3/YOHTowbNgwXnzxxSLP7eXlhdc1N9d+/frh4+PDvHnz8nWgviI5Ofm6/YauFRwcjLOzMxkZGUWWi4uL47vvvuOuu+5Co5EKUIANJ/bzwrYplxOhhqy4+2uCSvi+CyFEdSbJUA20Zs0akpKSePjhhzEajfn23X333SxYsKBAMgQwa9YsmjVrhk5Xuj8bFxcXvvrqK0aMGMGQIUOYMmUKDRo0ID4+nh9++IGzZ8/y/fffFzhu5syZZGZmMmjQIOrVq0dycjIffPABJpOJvn372sqpqkpMTAyqqpKcnMyOHTuYPXs2RqORuXPnlirW6mrLyUNMvdxHyMEcwoq7F0giJIQQl8lX5hpowYIF9OnTp0AiBHk1Q+Hh4SQmJhbY16hRI8aNG0d2dnapr3nnnXeyfft29Ho9o0aNonHjxtx3332kpKTw5ptvFnpMjx49OHXqFA888ACNGzdm4MCBxMTE8McffxAWFmYrl5qaSmBgILVr16Zz5858/vnnPPjgg+zfv5/AwMBSx1rdbDt9jMlbHgNtOnpzEMuHLSTIw9PeYQkhRKWhqDfau7WGSE1NxWg0kpKSgrt7/hW7s7OziYqKIiQkBEdHRztFKKqTsv6b2nf+NGPXjUXVJaCzBLD8rm8J8fIvg0iFEKJyK+r+fS2pGRKimjoaF81D6x5B1SWgtfiwbMgiSYSEEKIQkgwJUQ2dToznvtXjsOpi0Vg8WDxwAY18ahd/oBBC1ECSDAlRzUSnpjD8l3GYdedRLK583vdLWgWG2jssIYSotCQZEqIaic9IZ8iPD5OriwKLM+/3+JROQaWb+VsIIWoaSYaEqCbSsrMZvGw82bpIsBqY0/UDbqvf2t5hCSFEpSfJkBDVQK7ZzOBlE0nXHkK16nml3XzuCOto77CEEKJKkGRIiCrOarUy7IdnSOAfVKuWp1rM4t4WPe0dlhBCVBmSDAlRxY1d8SZnTJtQVYUx9V/ikXYD7R2SEEJUKZIMCVGFPbP2E/an/QhAP//HeaH7PcUcIYQQ4lpVJhlKSkpizJgxGI1GjEYjY8aMITk5uchj0tPTmTRpEnXq1MHJyYkmTZrw6aefVkzAosROnz6NoiiEh4dX2DUXLVpU4sVhr+jZsydPPfVUucRzI+ZtXca6mM8AaOV6D+8OfNzOEQkhRNVUZZKhUaNGER4eztq1a1m7di3h4eGMGTOmyGOmTp3K2rVr+fbbb4mIiGDq1KlMnjyZlStXVlDUJWO1Wjh3+AAR2/7k3OEDWK2WcruWoihFPsaOHWsrt2LFinKLo7S2bNmCoijFJsAlde+993Ls2LFSHbN8+XLeeOONMrn+zVqyfz1LTs5BUVTq6vqwZOjL9g5JCCGqrCqxan1ERARr165l586ddOyYN0Lmyy+/pHPnzkRGRuZbtPNqO3bs4MEHH6Rnz54ATJgwgc8//5w9e/Zw5513VlT4RTq+azubFn1BemK8bZurlw+3jZ1Aw45dyvx60dHRtt+XLVvGjBkziIyMtG1zcnIq82tWpNzcXBwcHIot5+TkVOrX6uXldaNhlakNJ8J5O3waisaCp9qO5ffMQ6OpMt9rhBCi0qkSn6A7duzAaDTaEiGATp06YTQa2b59+3WP69atG6tWreLChQuoqsrmzZs5duwY/fv3v+4xOTk5pKam5nuUl+O7trPq3dn5EiGA9MR4Vr07m+O7rv/ablRAQIDtYTQaURSlwLYbMXPmTOrWrYvBYKBWrVpMmTLFtq+wWiYPDw8WLVqUb9vRo0fp0qULjo6ONGvWjC1btgB5zWi9evUCwNPTM18NVs+ePZk0aRJPP/00Pj4+9O3bF4B3332XFi1a4OLiQlBQEBMnTiQ9Pd12rWubyWbOnEnr1q355ptvCA4Oxmg0MnLkSNLS0mxlrm0mCw4OZvbs2YwbNw43Nzfq1q3LF198ke81bd++ndatW+Po6Ei7du1YsWLFTTUJHo+/wNN/TQFNDk6WRqy+9xMMev0NnUsIIUSeKpEMxcTE4OfnV2C7n58fMTEx1z3ugw8+oGnTptSpUwcHBwcGDBjAJ598Qrdu3a57zJw5c2z9koxGI0FBQWXyGq5ltVrYtOiLIstsXvxFuTaZlZWffvqJ9957j88//5zjx4+zYsUKWrRoUerzPPfcczzzzDPs37+fLl26MGTIEBISEggKCuLnn38GIDIykujoaP73v//Zjlu8eDE6nY5t27bx+eefA6DRaPjggw84dOgQixcvZtOmTTz//PNFXv/kyZOsWLGCNWvWsGbNGv7880/mzp1b5DHvvPMO7dq1Y//+/UycOJHHH3+co0ePApCWlsbgwYNp0aIF+/bt44033uCFF14o9ftyRXJWOqNWPYqqTUJj9uOHuz7DWMVr8oQQojKwazI0c+bMYvuw7NmzB8irXbiWqqqFbr/igw8+YOfOnaxatYq9e/fyzjvvMHHiRDZs2HDdY6ZNm0ZKSortce7cuZt/oYW4EHG4QI3QtdIS4rkQcbhcrl+Wzp49S0BAAH369KFu3bp06NCB8ePHl/o8kyZNYvjw4baO7kajkQULFqDVam1NVH5+fgVqsBo0aMC8efMICwujceO8pSeeeuopevXqRUhICLfddhtvvPEGP/zwQ5HXt1qtLFq0iObNm3PrrbcyZswYNm7cWOQxgwYNYuLEiTRo0IAXXngBHx8fW43W0qVLURSFL7/8kqZNmzJw4ECee+65Ur8vABarhWE/TSZbewYsznx020cEe/ne0LmEEELkZ9c+Q5MmTWLkyJFFlgkODubAgQPExsYW2Hfp0iX8/f0LPS4rK4uXXnqJX375hdtvvx2Ali1bEh4ezvz58+nTp0+hxxkMBgwGQylfSemlJyeVaTl7GjFiBO+//z6hoaEMGDCAQYMGMXjwYHS60v15de7c2fa7TqejXbt2REREFHtcu3btCmzbvHkzs2fP5siRI6SmpmI2m8nOziYjIwMXF5dCzxMcHIybm5vteWBgIHFxcUVeu2XLlrbfrzQ5XjkmMjKSli1b4ujoaCvToUOHYl9PYcb8MpNL1j2oVi3PtprDrSFNbug8QgghCrJrzZCPjw+NGzcu8uHo6Ejnzp1JSUlh9+7dtmN37dpFSkoKXboU3snYZDJhMpkKdCzVarVYrdZyfV0l4erhWabl7CkoKIjIyEg+/vhjnJycmDhxIt27d8dkMgF5SYKqqvmOubKvOEXV/F1xbXJz5swZBg0aRPPmzfn555/Zu3cvH3/8cbHX1V/T90ZRlGL/Voo6prCay2vfh5KYvvFrDqavAODOOk8xtu1tpT6HEEKI66sSfYaaNGnCgAEDGD9+PDt37mTnzp2MHz+eO+64I99IssaNG/PLL78A4O7uTo8ePXjuuefYsmULUVFRLFq0iCVLljB06FB7vRSb2k2a4erlU2QZN28fajdpVkER3RwnJyeGDBnCBx98wJYtW9ixYwcHDx4EwNfXN98otuPHj5OZmVngHDt37rT9bjab2bt3r63Z68oIMYul+D5Ue/bswWw2884779CpUycaNWrExYsXb+r13YjGjRtz4MABcnJy8sVWGt/s38wv5/L6RzV3vptZfceWZYhCCCGoIskQ5PW/aNGiBf369aNfv360bNmSb775Jl+ZyMhIUlJSbM+///572rdvz+jRo2natClz585l1qxZPPbYYxUdfgEajZbbxk4oskyvByeg0WgrKKKCoqKiCA8Pz/e4ekTWFYsWLWLBggUcOnSIU6dO8c033+Dk5ES9evUAuO222/joo4/Yt28fe/bs4bHHHitQowLw8ccf88svv3D06FGeeOIJkpKSGDduHAD16tVDURTWrFnDpUuXCo3jivr162M2m/nwww9t8Xz22Wdl9K6U3KhRo7BarUyYMIGIiAjWrVvH/PnzgZLVeO05f5J5+19CUaz40IFvh00v75CFEKJGqhLzDEHeHC/ffvttkWWubYIICAhg4cKF5RnWTWnYsQtDnn6pwDxDbt4+9HqwfOYZKo2nn366wLbNmzfb5m26wsPDg7lz5/L0009jsVho0aIFq1evxtvbG8gbcfXQQw/RvXt3atWqxf/+9z/27t1b4Nxz587lrbfeYv/+/dSvX5+VK1fi45NXe1a7dm1ee+01XnzxRR566CEeeOCBAkPzr2jdujXvvvsub731FtOmTaN79+7MmTOHBx544ObekFJyd3dn9erVPP7447Ru3ZoWLVowY8YMRo0ala8fUWGSMjOYsG4S6NJxsASxfOT/0GqrzHcXIYSoUhT1Rjox1CCpqakYjUZSUlJwd3fPty87O5uoqChCQkKKvbkVxWq15I0uS07C1cOT2k2a2bVGSJSfpUuX8tBDD5GSklLopI/Z2dmcOnWKmQcWEJHzB1hcWTrw/2gZGFzxwQohRBVW1P37WqWuGTp9+jRbt27l9OnTZGZm4uvrS5s2bejcufNNJQQ1mUajJahZy+ILiipnyZIlhIaGUrt2bf79919eeOEF7rnnniJnv07ISiVZPYiqanm57RxJhIQQopyVOBn67rvv+OCDD9i9ezd+fn7Url0bJycnEhMTOXnyJI6OjowePZoXXnjB1ldEiJouJiaGGTNmEBMTQ2BgICNGjGDWrFnXLR+fkUKuNQOA22s9zn2telZQpEIIUXOVKBm65ZZb0Gg0jB07lh9++IG6devm25+Tk8OOHTv4/vvvadeuHZ988gkjRowol4CFqEqef/75Yme+viI1O5Ok3EsA1NF34a1+j5ZnaEIIIS4rUZ+hX3/91TZxYXHi4+OJioqiffv2Nx1cZVARfYaEyDWbOJ50Cqs5l0vnk2jTuDWebkW3cQshhLi+Mu8zVNJECPImUrwyAkgIUTyrqnIy+SwoZlC1+Dh54aR3sHdYQghRY9zQWN2TJ0/yyiuvcN9999mWHli7di2HD1f+dbSEqGxOJ13ESjag4O8SgF4rIwmFEKIilToZ+vPPP2nRogW7du1i+fLltsnvDhw4wKuvvlrmAQpRncWmJ5FlTQbA08Efd4OzfQMSQogaqNTJ0Isvvsibb77J+vXrbUskAPTq1YsdO3aUaXBCVGcZudnEZ+ctU2JQ3Knl7m3niIQQomYqdTJ08ODBQtf28vX1JSEhoUyCEqK6M1stnEk9C6goGAjxrG3vkIQQosYqdTLk4eGRb9HNK/bv30/t2vKBfiNUq0r2yWQyw+PIPpmMaq15k4KfPn0aRVEIDw8HYMuWLSiKQnJy8nWPWbRoER4eHjd97bI6T0mpqsqppPOomAANwe5BaDWy1IYQQthLqT+BR40axQsvvEBMTAyKomC1Wtm2bRvPPvtsha/9VB1kHYon5q3dxH95kMTvI4n/8iAxb+0m61B88QffAEVRinyMHTvWVm7FihXlEkNJdOnShejoaIxGY5meNzg4mPfffz/ftnvvvZdjx46V6XWKciEtHpOa19fOz6kWzg6GCru2EEKIgkq9HMesWbMYO3YstWvXRlVVmjZtisViYdSoUbzyyivlEWO1lXUonoRvIwpst6TkkvBtBN73N8GpedlOU3B1rd6yZcuYMWMGkZGRtm1FLRNRkRwcHAgICKiQazk5OVXY607OziDl8sSKLlovfF3KNtkTQghReqWuGdLr9SxdupRjx47xww8/8O2333L06FG++eYbtDIkuMRUq0ry6pNFlklefarMm8wCAgJsD6PRiKIoBbaV1n333cfIkSPzbTOZTPj4+LBw4UIgb+qFbt264eHhgbe3N3fccQcnT17/9RfWTLZo0SLq1q2Ls7MzQ4cOLdBH7eTJk9x55534+/vj6upK+/bt2bBhg21/z549OXPmDFOnTrXVhF0577XNZJ9++in169fHwcGBsLAwvvnmm3z7FUXhq6++YujQoTg7O9OwYUNWrVpV5Ptkspi5kHYeUNHiRD2Pikn2hBBCFO2GOyoEBQXRokULhg0bRsOGDcsyphohJyoFS0pukWUsKTnkRKVUUEQ3bvTo0axatco2zQLAunXryMjIYPjw4QBkZGTw9NNP888//7Bx40Y0Gg1Dhw7FarWW6Bq7du1i3LhxTJw4kfDwcHr16sWbb76Zr0x6ejqDBg1iw4YN7N+/n/79+zN48GDOnj0LwPLly6lTpw6vv/460dHRhfZ9A/jll1948skneeaZZzh06BCPPvooDz30EJs3b85X7rXXXuOee+7hwIEDDBo0iNGjR5OYmFjoOfP6CV3Im1gRLSEeQbZkTAghhH2VOhnKzMzk4YcfxtnZmWbNmtluNFOmTGHu3LllHmB1ZU0rOhEqbTl76t+/Py4uLvzyyy+2bd999x2DBw+2TYE+fPhwW+LcunVrFixYwMGDBzly5EiJrvG///2P/v378+KLL9KoUSOmTJlC//7985Vp1aoVjz76KC1atKBhw4a8+eabhIaG2mpsvLy80Gq1uLm52WrCCjN//nzGjh3LxIkTadSoEU8//TTDhg1j/vz5+cqNHTuW++67jwYNGjB79mwyMjLYvXt3oee8kJqAmbxkMdC5NgadvkSvWwghRPkrdTI0bdo0/v33X7Zs2ZJvPa4+ffqwbNmyMg2uOtO4lWy5hZKWsye9Xs+IESNYunQpkFcLtHLlSkaPHm0rc/LkSUaNGkVoaCju7u6EhIQA2JLp4kRERNC5c+d82659npGRwfPPP0/Tpk3x8PDA1dWVo0ePlvgaV1+ra9eu+bZ17dqViIj8/btatmxp+93FxQU3NzfbjOxXS83OIsWUt91F54mXs1up4hFCCFG+St2BesWKFSxbtoxOnTrlq+Zv2rRpkX1ARH6GECNao0ORTWVaowFDSNXoYDt69Gh69OhBXFwc69evx9HRkYEDB9r2Dx48mKCgIL788ktq1aqF1WqlefPm5OaWrOarBOsJ89xzz7Fu3Trmz59PgwYNcHJy4u677y7xNa52bROWqqoFtun1+Wt3royuvJrZauV82nlQVDQYqGuUfkJCCFHZlLpm6NKlS/j5+RXYnpGRIX0gSkHRKHgMrl9kGY/BoSiaqvGedunShaCgIJYtW8bSpUsZMWKEbYbyhIQEIiIieOWVV+jduzdNmjQhKSmpVOdv2rQpO3fuzLft2udbt25l7NixDB06lBYtWhAQEMDp06fzlXFwcMBisRR5rSZNmvD333/n27Z9+3aaNGlSqpgBopIuoiq5gIZgYxAaReYTEkKIyqbUNUPt27fn119/ZfLkycB/36C//PLLAs0WomhOzX3wvr8JyatP5qsh0hoNeAwOLfNh9aUVFRVlmwTxigYNGuDq6lqgrKIojBo1is8++4xjx47l62zs6emJt7c3X3zxBYGBgZw9e5YXX3yxVLFMmTKFLl26MG/ePO666y7++OMP1q5dWyC25cuXM3jwYBRFYfr06QVqaoKDg/nrr78YOXIkBoMBH5+C7/Fzzz3HPffcwy233ELv3r1ZvXo1y5cvzzcyrSRi0pLIVfM6wPs6BeCkl/mEhBCiUlJLadu2baqbm5v62GOPqY6OjuqTTz6p9unTR3VxcVH37NlT2tNVeikpKSqgpqSkFNiXlZWlHjlyRM3Kyrqpa1gtVjXrRJKasT9WzTqRpFot1ps6X0ktXLhQNRqNhe4DCn1s3rz5uuc7fPiwCqj16tVTrdb8r2H9+vVqkyZNVIPBoLZs2VLdsmWLCqi//PKLqqqqGhUVpQLq/v37VVVV1c2bN6uAmpSUZDvHggUL1Dp16qhOTk7q4MGD1fnz5+eLPyoqSu3Vq5fq5OSkBgUFqR999JHao0cP9cknn7SV2bFjh9qyZUvVYDCoV/78C3sfPvnkEzU0NFTV6/Vqo0aN1CVLlhR4f67EfoXRaFQXLlyoqqqqpudkq4cuHVEPXTqknkg4d9337Fpl9TclhBA1XVH372spqlqCzhjXOHToEG+//TZ79+7FarVyyy238MILL9CiRYuyydAqkdTUVIxGIykpKbaRUVdkZ2cTFRVFSEhIvs7komazWK1EJkShKtko6AnzalDi5Tbkb0oIIcpGUffva5WqmcxkMjFhwgSmT5/O4sWLbypIIaqrM8mxqEo2oFBP1h0TQohKr1Sf0nq9Pt9cMkKI/BIy08iy5k286Gnww8WhcixvIoQQ4vpK/ZV16NChdl3AU4jKKtdsJibjAgB6xZVAV287RySEEKIkSj2arEGDBrzxxhts376dtm3b4uLikm//lClTyiw4IaoKVVWJSj4PigXQEeJRW6aaEEKIKqLUydBXX32Fh4cHe/fuZe/evfn2KYpSI5OhG+iDLqqZ6LREzGQAUMulNnptqf9rAfK3JIQQ9lDqT+yoqKjyiKNKujIDcWZmJk5O0jekpsrIzSEpNxYAZ60nnk4F52EqqSuzZWu12jKJTQghRPFKnQy9/vrrPPvsszg7O+fbnpWVxdtvv82MGTPKLLjKTqvV4uHhYVuPytnZWZpGahirqnI6+RxWLCjo8XfxIDs7+8bOZbVy6dIlnJ2d0elurGZJCCFE6ZV6niGtVkt0dHSBJTkSEhLw8/MrdqmDqqa4eQpUVSUmJobk5OSKD07YXUJmKjnWdEDB29Hnplej12g0hISE2JYyEUIIcWPKbZ4hKHzBSoB///0XLy+v0p6uylMUhcDAQPz8/DCZTPYOR1SgzScP8W7kqyiKlV7+o5naostNn9PBwQGNzEskhBAVqsTJkKenJ4qioCgKjRo1ypcQWSwW0tPTeeyxx8olyKpAq9VKP48aJCU7k5d2voRZF4230poXeo6RJEYIIaqoEidD77//PqqqMm7cOF577TWMRqNtn4ODA8HBwbJQq6gxxq+ag1kXDRZXvr7zbUmEhBCiCitxMvTggw9iNpsB6NOnD3Xq1Cm3oISozH44sI0jmStRFBgX9hyh3gH2DkkIIcRNKNXXWZ1Ox8SJE6tdJ2khSiopK4PZ/8xEUVQCtJ2Z2nWYvUMSQghxk0pdt9+xY0f2799fHrEIUek9snIOFl0MisWVrwfPsnc4QgghykCpR5NNnDiRZ555hvPnzxe6HEfLli3LLDghKpOl+/8mMnsVigKPNHmOIKOvvUMSQghRBko9z1BhHUUVRbENua9uTWilmadAVF8JGRnc9n93YdXHUFvXmbWjv7B3SEIIIYpQrvMMyXIcoiZ6eOUcrPq85rEFQ6V5TAghqpNSJ0P16tUrjziEqLQW79nKidzVKAo81vx5artL85gQQlQnN7wA0pEjRzh79qxtYckrhgwZctNBCVFZXErL4J3wN1D0VoIcOjOxw1B7hySEEKKMlToZOnXqFEOHDuXgwYO2vkKAbUbq8uozlJSUxJQpU1i1ahWQl3R9+OGHeHh4XPeY2NhYXnjhBf744w+Sk5Pp3r07H374IQ0bNiyXGEX1M37VHFR9NIrVla/vmG3vcIQQQpSDUg+tf/LJJwkJCSE2NhZnZ2cOHz7MX3/9Rbt27diyZUs5hJhn1KhRhIeHs3btWtauXUt4eDhjxoy5bnlVVbnrrrs4deoUK1euZP/+/dSrV48+ffqQkZFRbnGK6uPbfds4YVoNwMTmzxPg5mPniIQQQpSHUo8m8/HxYdOmTbRs2RKj0cju3bsJCwtj06ZNPPPMM+UyB1FERARNmzZl586ddOzYEYCdO3fSuXNnjh49SlhYWIFjjh07RlhYGIcOHaJZs2ZAXq2Vn58fb731Fo888kiJri2jyWqm5Mwsenx3J1Z9NHX0nfh91Jf2DkkIIUQplOb+XeqaIYvFgqurK5CXGF28eBHI61gdGRl5A+EWb8eOHRiNRlsiBNCpUyeMRiPbt28v9JicnBwAHB0dbdu0Wi0ODg78/fff171WTk4Oqamp+R6i5hm/ei5WfTSKxZUvb5fRY0IIUZ2VOhlq3rw5Bw4cAPJmo543bx7btm3j9ddfJzQ0tMwDBIiJicHPz6/Adj8/P2JiYgo9pnHjxtSrV49p06aRlJREbm4uc+fOJSYmhujo6Otea86cORiNRtsjKCiozF6HqBp+OriTiKwVAIxr8gx1jAX/9oQQQlQfpU6GXnnlFaxWKwBvvvkmZ86c4dZbb+W3337jgw8+KNW5Zs6ciaIoRT727NkD/NdB+2pXJnosjF6v5+eff+bYsWN4eXnh7OzMli1bGDhwIFqt9roxTZs2jZSUFNvj3LlzpXpNompLy87mzd2voihWArQdeKrz3fYOSQghRDkr9Wiy/v37234PDQ3lyJEjJCYm4unped3E5HomTZrEyJEjiywTHBzMgQMHiI2NLbDv0qVL+Pv7X/fYtm3bEh4eTkpKCrm5ufj6+tKxY0fatWt33WMMBgMGg6HkL0JUK4+veQeL7iJYXPhyyBx7hyOEEKIC3PA8Q1fz8vK6oeN8fHzw8Sl+hE7nzp1JSUlh9+7ddOjQAYBdu3aRkpJCly5dij3eaDQCcPz4cfbs2cMbb7xxQ/GK6m39sYOEp/2EooFRDaYQ7CnNY0IIUROUejRZRkYGc+fOZePGjcTFxdmazK44depUmQZ4xcCBA7l48SKff/45ABMmTKBevXqsXr3aVqZx48bMmTOHoUPzJsb78ccf8fX1pW7duhw8eJAnn3yStm3b8vPPP5f4ujKarGbINVvosmgEOfrjeCot+HPM0lLXdAohhKg8ynVtskceeYQ///yTMWPGEBgYWGE3jKVLlzJlyhT69esH5E26+NFHH+UrExkZSUpKiu15dHQ0Tz/9NLGxsQQGBvLAAw8wffr0ColXVC3P/P4VOfrjYNXz0cA3JRESQogapNQ1Qx4eHvz666907dq1vGKqVKRmqPrbe/4sD/4xAkWbSb+Ah3mn/1P2DkkIIcRNKtd5hjw9PW+4j5AQlY2qqjz5x+so2kyc1CDm9nnC3iEJIYSoYKVOht544w1mzJhBZmZmecQjRIWau2UVKdpdqKrCnO6vo9fq7R2SEEKIClaiPkNt2rTJ14fixIkT+Pv7ExwcjF6f/+axb9++so1QiHJyPjmF706+C3pobbyd3qHXn3JBCCFE9VWiZOiuu+4q5zCEqHiPrZkD+kS0Vk8+GviSvcMRQghhJyVKhl599dXyjkOICvV/4Ts5bf4dRYEprV/Aw9HN3iEJIYSwk1L3GSqJUg5QE6JCZeTk8taeN1AUK7X1HRnX5nZ7hySEEMKOSpQMNWnShO+++47c3Nwiyx0/fpzHH3+ct956q0yCE6I8TP71Yyz6s2B15NOBr9s7HCGEEHZWomayjz/+mBdeeIEnnniCfv360a5dO2rVqoWjoyNJSUkcOXKEv//+myNHjjBp0iQmTpxY3nELcUO2RZ1gd8q3KBoYFjyBEM9a9g5JCCGEnZVq0sXt27ezbNky/vrrL06fPk1WVhY+Pj60adOG/v37c//99+Ph4VGO4VY8mXSx+rBYrHRd9AAZun9xpQHbHvgZjVIuLcVCCCHsrNyW4+jSpUuJFkYVojKauWkZGbp/UVUt79z2hiRCQgghgHLqQC1EZXMqIZ4VZ/PWsuvkNZwudZvbOSIhhBCVhSRDokZ47Lc3QZeKzurHBwOfs3c4QgghKhFJhkS1t2DPZi5aNgHwYrtXcNY72jkiIYQQlYkkQ6JaS8vO5sN/56AoKiGGntzbope9QxJCCFHJSDIkqrXJv7+HRRcNFhc+HSQzqQshhCio1MlQz549WbJkCVlZWeURjxBlZsfZo+xJ/gGAe0ImUdvdx84RCSGEqIxKnQy1bduW559/noCAAMaPH8/OnTvLIy4hborVauXZTa+haMy4Wpvyco9R9g5JCCFEJVXqZOidd97hwoULLFmyhEuXLtG9e3eaNm3K/PnziY2NLY8YhSi1d7f/TKpyCNWqZV6vmWg00iIshBCicDd0h9Bqtdx5552sWLGCCxcuMGrUKKZPn05QUBB33XUXmzZtKus4hSixxMx0lhz7AICWbndxa3ATO0ckhBCiMrupr8u7d+9mxowZzJ8/Hz8/P6ZNm4afnx+DBw/m2WefLasYhSiVib/NRdUmo5i9+GjQM/YORwghRCVXquU4AOLi4vjmm29YuHAhx48fZ/DgwXz//ff0798fRVEAuOeee7jrrruYP39+mQcsRFG2RB3iUPpqFAXGNn4KL2c3e4ckhBCikit1MlSnTh3q16/PuHHjGDt2LL6+vgXKdOjQgfbt25dJgEKUlNVqZdqW11A0VoxqK6Z2ucveIQkhhKgCSp0Mbdy4kVtvvbXIMu7u7mzevPmGgxLiRsz+6/9I1xxFtep4p8+rtppKIYQQoiil7jNUXCIkhD3Epifzw6lPAGhrHEHHoIZ2jkgIIURVUeqaoTZt2hT6jVtRFBwdHWnQoAFjx46lVy9Z9kBUnIm/zkHVpqIx+/LhoKn2DkcIIUQVUuqaoQEDBnDq1ClcXFzo1asXPXv2xNXVlZMnT9K+fXuio6Pp06cPK1euLI94hShg3fH9RGatBeDR5s/i7uhk54iEEEJUJaWuGYqPj+eZZ55h+vTp+ba/+eabnDlzhj/++INXX32VN954gzvvvLPMAhWiMBarhel/53Wa9qYdEzsOsndINYJqVYk7k0bs6VQyU3JQtAoevk4ENvTA3VuSUSFE1aKoqqqW5gCj0cjevXtp0KBBvu0nTpygbdu2pKSkcPToUdq3b09aWlqZBmsPqampGI1GUlJScHd3t3c44hrTNy5kxfl3Ua0OfDvgZ1oHBts7pGrNYrFy+K+LhG84S1pCdqFlvGq50LJXHcI6BqBz0FZwhEIIkac09+9S1ww5Ojqyffv2AsnQ9u3bcXR0BPKGOBsMhtKeWohSOZeSwIozX4AWOnuNlESonCXHZrL2i4MkXMgAQO+opXZDD9x8nLCarcSfTyfuTBqJFzPYsjSS3Wui6DKsAY06+MvIPiFEpVbqZGjy5Mk89thj7N27l/bt26MoCrt37+arr77ipZdeAmDdunW0adOmzIMV4mpP/PYmaNPRmgN4f+AUe4dTrZ0+GM/6BYfJzbbg6Kqn4+AQGncOLFDzk5NpImJ7NAc2nSctMZsNC49w5O+L9B7bRJrPhBCVVqmbyQCWLl3KRx99RGRkJABhYWFMnjyZUaPyVgbPysqyjS6r6qSZrHJaGbGLl3eNR1FUpjZ/j3Ft+9g7pGrrVPgl1n5+EFWFwPpG+k9ojoux6Jpfs8nCvxvPsee305hzrTg46eg5OoyG7fwrKOrrU60qOVlmcrPMWK0qWp0Gg7MOB8dSfzcUQlRipbl/lyoZMpvNzJo1i3HjxhEUFHTTgVYFkgxVPiazmc5L7iJHewY/pTMbH/jC3iFVW+cjk1jz4b9YzFYadfTntjFN0OpKPgg15VIW678+TGxUKgCNuwTSfWQj9BXYlyg5NpNzEYlEn0gm4WIGKXFZWMzWAuX0jlrcfZzwreuGX103/EPc8Q1yQ9FIE58QVVG5JUMArq6uHDp0iODg4JuJscqQZKjyeX7d5/we8xGqxZEfbv+Fpv517B1StZQcm8kPc/7BlG0hpJUPAyY0R6Mt/drOFouVPb+eZu/vp1HVvA7W/R9pjlctl3KIOk9ulpkj2y4SuSuG+HPphZbR6jVoNAoWsxWrpfCPQSd3B+o19yakpQ/1mnuXKhEUQthXuXag7tOnD1u2bGHs2LE3Gp8QNywqMZbfLywALfT0u18SoXJiMVn5Y8FhTNkWAhsY6fdIsxtKhAC0Wg0dh4RSJ8yTPxYcJvFiBj/O/Yeeo8II6xRYpnGnJ+VwYNM5Dm+9QG62BQCNRiGwoZHajTzxC3bHw88ZV09DvsQmN9tMRnIOSTGZXDqbRtyZNKJPJJOVmsvR7dEc3R6No4uehh38adIlEN8gWQBYiOqk1DVDn3/+OTNnzmT06NG0bdsWF5f83+6GDBlSpgHam9QMVS6DvpvEOdOf6Mx12P7ASpz0DvYOqVra+sMxDmw6j6OLnntf6YCrZ9mMDs1MzWX914c5fzQJgMadA+h2TyMMTjfXXyfhYjrh689ybHesrZbHM8CZlrcF0eAWPxxd9aU+p8Vs5eKJZM4cSOD43lgyU3Jt+2o19KBV7yCCW/qgkWY0ISqlcm0m02iu/+1QURQsFktpTlfpSTJUefxw8C/e2PcEANNaf8yoVt3tHFH1dD4yiZXv7Qfg9idaEtzCp0zPb7Wq7P39NLvXRIEKLkYHuo8MI6S1T6mG4KuqyoXIJMI3nOPMoQTb9loNPWjTty71mnuXWX8fq8XKuaNJHN0ezan9l7Ba8z423X0cadkriCZdA6UDthCVTLkmQzWNJEOVQ47ZROclQzBpz1NL2511939s75CqJbPJwvdv7CYlLotm3WvTc1RYuV3r4vFkNi2JIOVSFgD+Ie60vz2EoKZeRda2ZKXnErkzhsNbL5Icm5m3UYH6rX1p3a8uASHGcosZ8priDv55nsNbL5CTYQbA4KyjeY/atOwVhLO71FYKURlUWDKUnZ1dLYbPF0WSocrhqd8+YOOlL8HixPIhq2joE2DvkKqlXatPsefX0zi7OzDqtU433XxVHHOuhT2/nebfjecwm/JGeDkbHQhp6YNvXTdcPR1Bgay0XJKiM7l4PInYqFSufGrpDVrCOgbQqncQHv7O5RrrtUy5FiJ3xhC+4SwpcXkJnVanoXHnAFr3qVvh8Qgh8ivXZMhisTB79mw+++wzYmNjOXbsGKGhoUyfPp3g4GAefvjhmwq+spFkyP4iL13g7jV3gSabvn4TeXfg4/YOqVpKjsvk/17bhdWi0n98cxq09auwa2ek5LB/3VmO7oq21bYUxbeuG81urUXD9v52b56yWlWi/r3E/j/O2qYQQIHQ1r60qYCaKiFE4cp1NNmsWbNYvHgx8+bNY/z48bbtLVq04L333qt2yZCwv8nrXgdNNg6WesztO774A8QN2fnLSawWlbpNvah/i2+FXtvFaKDbPQ3pPKw+Zw8nEH0yhfhzaWRnmFFVFYOzDk9/F3yCXKnbzBs3r8pTI63RKNRv40doa1+iT6Sw/48znD6YwKn9lzi1/5Kts3W9Ft5ob3BEXmlZrRYuRBwmPTkJVw9PajdphkZz43M7qWYzOSdOkHP8OLnnzmFJSMSafnntSZ0OrYcHOl9fDA0aYmjYAJ2fX5H9v8o6vrJW2eNTrSo5USlY03LRuDlgCDFWqvmwSvv+qbm5qBYLGif7zVJf6mRoyZIlfPHFF/Tu3ZvHHnvMtr1ly5YcPXq0TIMT4pvwjURbtqOqCi91eBkHnXRSLQ/RJ5I5uf8SigJdhjew21piWp2GkFa+hLSq2GSstK73YV+roQe1GnqQeDGD/RvOcmxXDBePJ3PxeDJObnoatvenUfsA/OqV32SOx3dtZ9OiL0hPjLdtc/Xy4baxE2jYsUuJz2NJTydj61bSNm0m/a+/sKaklPhYracnzh074tKpEy5du+Bw1SS9ZRVfeans8WUdiid59UksV41u1Bod8BhcH6fmZTvY4UaU5P1TVZWcY8fJ2PoX6du2kbU/HP8XX8Bz5Eh7hV36ZjInJyeOHj1KvXr1cHNz499//yU0NJQjR47QoUMH0tMLn+DsZs2aNYtff/2V8PBwHBwcSE5OLvYYVVV57bXX+OKLL0hKSqJjx458/PHHNGvWrMTXlWYy+8ky5dDlmzswa2Oop+/DmlHv2TukaklVVX6et5fYqFSadqtFr/sb2zukSq00N8v0pBwObD7H0R3RZKWZbNud3PTUbepNnSae+NV1xyPAuUyG6B/ftZ1V786+7v4hT79U5A3ddOECaX/+SfqmzWTs2gWm/2LWuLpiaNQIh9AQdD4+aN3cQNGg5uZiSU7GdPEiOSdOkHvmDFwzqtjQpAnu/fsRV9uf379dcMPxlbebff/KW9aheBK+jbjufu/7m9g1ISru/es78C58T54lfetWzDEx+fYZhw+j1qxZZRpPuTaTNWvWjK1bt1KvXr1823/88cdyXZw1NzeXESNG0LlzZxYsuP5/pqvNmzePd999l0WLFtGoUSPefPNN+vbtS2RkJG5uMmlaZff02g8xa2PA4sonQ162dzjV1umDCcRGpaIzaOkwOMTe4QCVt5nieh/26YnxrHp3doGbpaungS7DGtDpzlDOHkkkcmcMZw4lkJVmInJXDJG78m4IOoMWr0AXjD6OuHk74eppwOCiw9FZj8FFj8FZh6OLHgcn3XWTJqvVwqZFRS9Ns3nxF4Q0b4n1UjymmFjMsTHknjtPblQUWQcPYjp7Nl95h3r1cO3dG7fet+HUujWKtvh/A2tODtmHD5OxYwcZO3aQtT+cnIgI4iIi2NykHuh1cJ28b/PiL6jfvmO5/1urqoolPh5TXBxqVhbWrCwsubls/K7oe0tFxVcY1aqSvPpkkWWSV5/CsWnZTSlRGsX+/akqW1f+SK+IsyiAYjDg3LEDrt1uxaVzJxwaNKiwWAtT6mTo1VdfZcyYMVy4cAGr1cry5cuJjIxkyZIlrFmzpjxiBOC1114DYNGiRSUqr6oq77//Pi+//DLDhg0DYPHixfj7+/Pdd9/x6KOPlleoogwcij3D1vj/Q9HAHXUeoa6n/at/qyNVVdnzaxQALXvWKXYB1opQGZopVFUFkwlrrgk1Jws1KxNLZhqbvvyoyOM2ffk/gg0J6Aw6FCxgNYHFhMZqJthqIfgWC5YWFqJjHTh7zomYOEcuJTpizoG406nEnU4tLjIMOhMOmhwclGwc1Cz0ahaO1nRMOedIT4sv8ui0hHh2dr8V74zswgso4FTbCbeGLrg2dMHgrQf1Tzj4JxxQQbUCKnnD+a75eXmfRlVxRsVZq+LbDcxtjaSfzOVEFGQ7FH3LSUuI58Lc2wjyvNyJvkDDxTXPC23YyL/NnGkl+5KF7DgL2XFWchItmNJU1GumxEtwcSSjQe3i45vVlSBjTiGXvflYr316tZzchlhSphYZnyUlh5y3h+PocLzIcoW7uVl2LqQ4kp5YxChfRSHbQY+pmRP16+bgXMuERrcdErbDGqD9I9DtqZuK4WaUOhkaPHgwy5YtY/bs2SiKwowZM7jllltYvXo1ffv2LY8Yb0hUVBQxMTH069fPts1gMNCjRw+2b99+3WQoJyeHnJz//tBTU4v7cBLlYcq611A0uTha6vNG77H2DqfaOns4kbgzaegcNLTuY//Fl0tT82LNycGScAlL3EUs8TFY4uOwpiZhzUjFmpGOmpmBNSsTNSsTa3Y21uxc1JxcrDkmVJMZq8mCaragmq15D4t61aPgN+sEF0fSi7lZpqdlsOvhGXnJhkZFo1XR6PIeivaqn1qVulqVYK0KOsgwBJDmEEim1od0jS/ZGg9yFRdyFde8nxoXLFpHQCHH7EAODkD+2m2L1QU4Vex7nKPXoeis6J0t6Jzyfjq4m3E0mnDyzUWrv3xTTL38uEk6wMMHXPS+cNG/2PLxkXHUqh2N1qF0N2dVBXO2hpxkPdlJerIT8x6mzOvc5hQVncGKRp/3b2J1KVmn/At7M/GqFY2zTy5KBS5VZ7UEl6xcugm0Z4svWIZUFZIuBgDFT3niUDcVV+MlyLhmR3bJ+6SVhxvqjdq/f3/69+9f1rGUqZjL7ZH+/vn/8/n7+3PmzJnrHjdnzhxbLZSwj6/++Z1L6j+oqoaZXaajK0HVvCg9VVX553KtUPPutXFys9NkgaqKNTEa0/kTbPy8iJoXVWX926/DmTNYslRUc3k1BRRyXkUlx1CyO1+O/vLHqlXBalWwmoounycVJ1JxIhLv65SwKlrMOmdMemesbh5Y3T2wuBgxO7iS7WAkUbFwsQRXchk1kvodvNE7aAAFFOWan5prtpH3U7le+aL2/Xdd11PnYMHPxcaXcciRY7sC0ft5YQiuhUNQADqjG1qjGxpnR7BYUS0WrFk5mBOSMSWkYIq5RM7pC1hSr73D5tHX8sOpYT0cG9bDMTQIfaAfel9PlKsGZOhPneWfz5cWG5/5nJazR33Qehnx6NMV44DuGOoE/Pc+5XPN80IHJlxbpvDras5b4JdCaqSuLXfHdKhzo5+ZpZsBPjvyFKl/7iLtz12kpqdDCVq6XO96CxoU0hzvWnyiXJ5ueGhObm4ucXFxWK3WfNvr1q1b4nPMnDmz2MTjn3/+oV27djcUI1BgVIyqqkWOlJk2bRpPP/207XlqaipBQfb/xlxTpOdk89HBt0ELDR37c3vjtvYOqdo6fzRvAkOtXkPrviX/f1sSqqpiTUrCdPIQ5rORmKPPYLkUgyUhAXNSMpaUdMxp2VgyzJiz8mpiim2mUBSytDpirY54my838ygqWgcrWkfQOipoDVoUgx6NQY/G4IDiaEDj6IjGyQnF0RGNkyOKoxMaRycURycUR+e8nwZnFCcXFCdXNJd/Kk6uKM6uKI4uKA6OuB6NJHz2zGJfe4PPP6dOvVCsWVlYs7Lzaqeys7FmZuX9npODNTsbNTsHa3ZWvp+qxYyid0DR61H0ejROjmi9vdFdfmi9vdH5+qIxFGzOtFotfPnEw/maFwu+h27sOdyKf0/oqNfMm9A2vgS38MGhnCfXBKgdasH1lz+LjM9Z50CAbwCmjNOY4hIxxSXC7kMlv4hWi0NIMI6NwnBs1izv0bQJ2hIMfqldtxOuP64rMj4XVzeC+w0kY9MmLIkpJPzwGwk//IZzu3Z4jLgbt379ym14uCFQRbtpd75RZNfSGg0Y2reHcuozpJrNZP37L2mbNpG2dh2mCxds+7ydnXHS6MiyXn+eMDdvH2rfehdUgv5/1yr1/4Djx48zbtw4tm/fnm/7lSSjNGuTTZo0iZHFDKULDg4ubYgABATkZeoxMTEEBv63MnZcXFyB2qKrGQwGDIV80IiKMeX397BoL4HFnY8HTbN3ONWW1Wph6/ebseTGEtw8GCe30t8MLYkJ5B7eTW5EOLmnj2OOjsF0KRFTUgamFDNq8XMnXpb3wZ3jULIPSPeHR1C/Uye0voFovPxQHN1BW/qFWEurdos2uHr5FHmzdPP2IahFKzQaLVoPj3KP6WoajZbbxk4ocjRPgw7DSb7kRHpSDif3X+Lk/ktodRrqNvOi/i1+hLQsv8SoJPH1mfIsDTp2wZKcTM7x42QfO4bp7FnMSUlYkpKxZmTkdeLWadE4OqHz80Mf4I8uMBDHRo1wqF+/0ESxrOLrPWEydTp2Qc3NJX3rVpJ//In0v/4ic88eMvfsQfPmLNzvuB2Pu+/GqRSjlktC0Sh4DK5f5Ggyj8GhZd552nzpEulb/yZ9619kbNuO9aquI4qTE269euLWfwCu3W9Fd2B/ke9frwcnVIqBEIUp9dD6rl27otPpePHFFwkMDCxQy9KqVasyDfBaixYt4qmnnip2aL2qqtSqVYupU6fy/PPPA3m1WX5+frz11lsl7kAtQ+srzj/nT/DQ+hEoGjMjgl5gxm332zukaun4ru2s/+ozslITbduu10FZzc0l9/Aecg5sJ/d4BLlnz5F7MYHc+Ews1+mDezWtwYrOVYPOzRGtuys6T3e0Xl7ofPzQ+tdCFxCEtlYI2tr1uXDuPD++8VKx57xnxmyCmrUs9esuC5V96DUU3gHdzduHXg/m/fuqqsqls2mcupwM2dZ3gwpJjIqLz95KG58pNpaUX34h+aefMZ0/b9tuaNIEj7uHY7zjDrTGspuFvPB5hgx4DA4tk2H1pgsXyNy7l8y9+8jcu4fcE/lHsGmNRly6dcOtXz9cu99aoCasMv37lutyHC4uLuzdu5fGjSt2LpKzZ8+SmJjIqlWrePvtt9m6dSsADRo0wNXVFYDGjRszZ84chg4dCsBbb73FnDlzWLhwIQ0bNmT27Nls2bKlVEPrJRmqOLcuHkMy4bhYw9j+4A9oNBXYO7GGKHYekK7tqZ10kZxjJ8k+d4nc+BxU6/W/aeqcLDh4O6D3MaIP8EMfWAt93WD0IU3QhTZH4xNU4irxkjTzuHn78MhHC+z67bIyfdhfT0mnJlBVlYQLGZzYG8vJfRWXGFXWqROuuJH4VKuVzF27SP7xJ9LWr0e9PEeTYjDg1rcv7gP649K1a5k0o5XVDNSmuDhyIiLIPnKE7CMRZB06hDk6ukA5x2bNcO3RHZdbb8WpZctip1ioLP++5ZoMtW/fnvfee49u3brdVJClNXbsWBYvXlxg++bNm+nZsyeQ1z9o4cKFjB07Fvhv0sXPP/8836SLzZs3L/F1JRmqGO9vX8GC49NRVS0f3voNveq3sHdI1U6xyYaq4mgy2+YBuUKjs+LgrcPg74FDUCAOIfVxCGuBQ/OOaPzrQxkmrVWh5gUqz4d9WSouMQpq4kmdxl7UauSBT23XSrX8Q2VjTkoidfUakn/6iZxjx2zbFScnXLt1w61Pb5w7dkQfUP4LTquqijkujtxTp8g5eSrvZ9Qpck6cwHKpkM8CnQ7Hpk1xbtsW53ZtcbrlFnSenuUeZ3ko12Ro06ZNvPLKK8yePZsWLVqg1+dvq69uCYMkQ+UvKTOdHt/fjqpNpKnznSwb8aa9Q6o2zImJZB+JIOffnZzds51NmcWPRullSiYoJADHpi0x3NIVfdPOKIaKW4G9KtS8VHdFJUYABmcdAfWN+NR2xbu2K56BLrh5GXBw0t30Ui5Wq4rFbMVisub/abZiMamoqorOQYve8N9D56Cx2xIyRVFVlexDh0hds4a09RswXcw/3k9fuzZObW/J6+/UoAGG+vXR+fujcSjZyE5VVbFmZGJJTsaSEI/p4sW8x4WLV/1+AWtG4aPs0GhwCAnBsWlT28OpRXM0zhX3/708lWsydKXp4nqjtErTgboqkGSo/I36aSYHM35GMXuw+b7f8HaW2cFvhDUjg8x9+8k6cIDs/bvIPhKBOTHNtv+ihyvh9YofvjpoynM06dqjPEMtVnWseamqriRGZw8ncOFYMtEnkjHlFP45rzdocXZ3wMFJh4OjFq1ec3lORvXyvIwqZltyo2IxWS7//C/xsVpLP/mfolFwdNHh6OqAk6seJzc9Tq4OOF79u5seZ7fL21z1aMph0dwCiZzZimpV0Wg1aDSQe+I4WVs3k/33X+RERMA1o7Gv0BqNaH180Dg5oXF0RHFwQDWbUU2mvJ/Z2XkJUHKyrTmuSFotDkFBOISGYggNwSG0PobQEAyNGlWbxKcw5bocx+bNm284MCGu9VfUYQ6krUDRwJhGT0oiVAqqxULWv/+S8fc2MnZsI+vAQbBc++Gqone14Ohpwa+WN5Sg07Orh/2rxDUard06SYv8FEXBp44rPnVcuaV/PawWK5fOphN7OpWEi+kkXkgnOTaL7AwTphwLKZeyyvDioNNp0Oo1aHQadDoNKGDOtWDKsWDOzft7V60qWWkmstJMJJXw1AZnHU5uDji66NA5aPNqmxw0aB20aHUauJLAXf5ptVgx51gx5Vow51gw5VowZVtsz825JU3kmqMNbIljAy0OGjMOpnQcsxIwJJ5DH3MSp/RoXNJi0KYUvfRGvrfJwQGtlxf6wED0tWqhr1378s9a6GvVwiEoCKWEtU01VamToR497PuNUVQfVquVF7a8hqKx4K624Jmuw+wdUqWnqipZe/aQ+vtaUtf+hiUxOd9+vbMZJ99cnHxUHMMaYujQE21YD6jdDv1FE7z4KKjXX0zZzduH2k3KdkiwqF40Wg3+Ie74h+T/pm3KtZCRlENmai652WZMORYsJuvl+ReVy3MzKmgvJzdanQbd5Z952/Lv0+o1aDRKkc1fqlXFlGshN8tCdkYuWWkmstNNZKXn/Z6VbiI7LZesdBNZl39mZ5hAhZxMMzmZJZ7/4YZodRoUrYJqUbFYrLYVLyxmKxkp1suTMLvkPYx1wdjVdqy7u4KnUcXDxYyHqxlvT3B206PR6/KSHw8PdJ6eaD08UJycKmUzYVVyQ8MDtm7dyueff86pU6f48ccfqV27Nt988w0hISEV3rFaVF1z/lpGuiYC1arj3d4zZfRYEUyxsaQs/4XkH/8P08U423aN3oprYA7O/jm4NK6Fwy19oUEfqNcFHPJXf+9bexC9cy9MGauve53KPA+IqNz0Dlo8/J3x8K+4ZhdFo+DgqMPBUYerZ8nmF7JaVXIyTJeTpVxyMsyYTf/VNJlNl5M4jVIgkdM7aNE7atE7aNFd6bPk8F+/Je1VyZ1GWzCRs1pVzLkWsjPykrbsDBOZqbmkxmeTGp9FanwWybGZZKWZSE1VSU2FM+i4cqt2ctPjW9cN3yA3fJ3d8HVww83JURKhMlDqZOjnn39mzJgxjB49mn379tnW8UpLS2P27Nn89ttvZR6kqH5i0pL5/tTHoIW2xuF0rNvI3iFVOqqqkrljB4lffET6rv22b5UanRW3oGzc65lx6doFpentUL83eFx/pvT48+mcCr+E1tCQbvc8zT8rl0gHZVEjaTQKTm4Ol5efcanwa19J3ty9rz/EPjM1l4Tz6cRfSCfhQjrx59JIjM5Lks4eTuTs4f/mCHN00eNb1zUvSarrjm9dV9x9pKaotErdgbpNmzZMnTqVBx54ADc3N/79919CQ0MJDw9nwIABtjXBqgvpQF0+hi97kWPZv6Ix+/DX6F8xOlbfTnylZc3NJXXZQhIXLSbnwn89IJx8cvBoZMW9Tw80Le+Ehn3BULI+Vuu+PMSJvXE0aOtH//HNpYOyEFWMOddCwoUMLp1N5dLZNOLOppF4MQOrpeAt3MFJdzlBykuOfIPccPd2yuvQXoOUawfqyMhIunfvXmC7u7t7sbNCCwHwe+Q+IrN+R1HgsebPSCJ0mTn6LMkfvU7i7zuwZOZ1DFW0VjwamPAc2BnDbQ9AaE/QlW65gcSLGZzYl9e01m5QMCAdlIWoanQO2gJ9tSwmKwkX07l0Ns32SLiQQW6WmQuRyVyITP7vBAq4uDvg5u2Em7cjbl6Ol0fa6XG8PPLO1plcr0Gn16LRFd1n61qqqmK15I2ou/LTYrZiNef1mbKaC2678rtnoDM+dew3gKbUyVBgYCAnTpwosGbY33//TWhoaFnFJaopi8XKjG1voGit+CjteLzjHfYOye5ytq8h8bP3Sdl7HtWS98Gjc7Lg2SkAz/sfQttuBBhcb/j8u9ecAhVC2/jiXfvGzyOEqFy0eg1+9dzxq3dVgmSxkhSdQdyZNOIv1yAlXMzAnGMhIyWXjJRcYk6llOj8ipJ3Da1eg3KdFe2vTYBuVNsB9apWMvToo4/y5JNP8vXXX6MoChcvXmTHjh08++yzzJgxozxiFNXIKxsXka09gWrV8/6AV+0dTrkprhlKzc0i4//eJvH/fibj9JU1hhQc/TR43XUb7mOfR/G6fh+gkrp0Lo2T+y6BAh3uCLnp8wkhKjetVoNPHbe8xOLy4DRVVcnOMJGWkJ3XWTshi/SknLxO3OmXR9ld7tBtNv036k1VyetUnlv4fEjFUvhvtKBOQaPN+5nXwfyqbXoFN2/HsnkDblCpk6Hnn3+elJQUevXqRXZ2Nt27d8dgMPDss88yadKk8ohRVBMnE+JYc/4L0EI3n/toFRhs75DKRWEzKF9ZCLV+w3qkfvoqiWv+JifpyjctFbem3niNG4/ToAdQynBU3e7VUQA0bOsntUJC1FCKouDk6oCTq0O+WqTCqKqK1axiNlkwm/ISIYv5+smQonBNYqNBq1XQ6PKmRqgqSt2B+orMzEyOHDmC1WqladOmtsVSqxvpQF12+n/zFBetG9FbAth2/xqcHErX96UqKG5trXYXLuIXnzcpnaIDjx7N8Zr8Mg6NW5d5LLGnU/lp7h4UBe57tSOeARU7ckYIIeypXDtQX+Hs7Ey7du1u9HBRwyzZu5ULlk0oCrzQ/qVqmQhZrRY2Lfri+gVUlUO+vvTNOY/3XX3xmDgdrad3ucWze/UpABp1DJBESAghinDDyZAQJZWancU7+2eh6FXqOnTl3ha97B1SubgQcfj6K8IDKArZDnocv/oG71a3lG8sx5I4ezgRRaPQ/vbgcr2WEEJUdZIM2cm60+tYeWIl3et0Z2TjkfYOp1xNXPM/rPoLKFZnPhv0ur3DKVNWq8r5o4mc2BvH2QN7SnRMRnpa8YVuMqa/fzwOQLNba2H0lakLhBCiKJIM2cn6YxFsjd5KXIparZOhjSeOEJ62DEUDI+s/QZDRz94hlQnVqnJ0Zwy7V58iPSlvFnaLqWSTFpb3QqhHt0cTfy4dBycdHQbLCDIhhCiOJEN2YrDWAVXhTPoJe4dSbswWKy/+OQNFZ8ZDacq0W8fYO6QykRybyfqvDxN3Jq+GR29Kxy9uLx5JRznkbyVHd/3RYHonD7zLcemRnEwTO1fmrXbd/vZgnFxlpWohhChOzZqbuxLp5tGc4QefwZjqREZuhr3DKRcv/vE12bpIsOr4X583q8VaOaf2X+KHWbuIO5OG1pxF/ZO/0HX7y7SK/5ZbuifQ7+FRRZ9AcyvL3tzD+cikosvdoG0/nSArzYSHvzMtetYpl2sIIUR1IzVDdqKLtOKbEUSf4w+wKeIAg1t1tndIZepI7AXWRn+JooUefvdzS62G9g7ppu376V92bEgAwJh8guZHvsbLMw7v4YE4j/sYJbQ77sAQ79AC8wy5efvQsu8oTv3rTlJMJivf30+7gcG0vz0YjbZsvpOcOZxAxPZoUKDX/Y3RFlFDJYQQ4j+SDNlJj3sasWdfJO45XhxdcYk7WqrVouYE8ibtevz3V1G0mThY6vBO/8n2DummmGLj+Pvt1RxJz+t/U+f8FtpkL8V3kANOo96CFnfnzTx2WcOOXajfvmOhM1C3HWRh6w/HiNgWzZ7fTnPhWBL9H2mOi8fNTTWQnW5iy7dHAWjZqw61Gnrc1PmEEKImkWTIThycdJxoFUmrf27BOdqVQ39eqDbNGnP//IVEZReqqjDr1tcw6KpmvxVLcjIJX33F/o3nOVFvCABh8Svo2vAnnAY/Ad2eAofC5++53kKoeoOW28Y0IaiJF1u+PUr0iRSWzdpNv0eaUyfsxjpWWyxW1n55iPSkHNx9neh0Z/0bOo8QQtRUUo9uR80a12dHvZUA7PjlJOlJ2XaO6OadjI/ju5PvANDKfTADGnawc0SlZ83MJP6zzznRtx+HVh+0JUJtlf+jzz3JOM3YAbe9fN1EqCQatvNnxLT2eNd2JSvNxKr397Pn99Oo1tJPCL/txxNciExCb9Ay6LEW6A0lG9UmhBAijyRDdjQ4rDOHArYS7RqFKcfC1mXH7R3STVFVlUd+nQG6VHQWPz67/SV7h1QqqsVC8k8/caJffy69/z4J2gCONr4fgFaem+k49V64Zwl41iuT63n4O3P3C21p3CUQVYVdK0+x5qN/SUssWVKsqiq7Vp3i4JbzAPR5qKmsPyaEEDdAkiE7auhdB43qztbQZaConAq/RNS/l+wd1g17e+svxLMNVVWY2fl13AxVZwmIzH/+IWrECKJfmY4lPh6zlwdHWj2MqmhpGBRH11dfRGnUp8yvq3PQ0vuBJvQak9fh+eyRRP7v9V0c2Hy+yMURc7PNbFwcwZ7fTgPQ9e4GhLb2LfP4hBCiJpA+Q3bm7xBGtGY354KiCTpbi20/n6Buc2+0ZTTCqKJEJV7im+PvgA5aug3mziZd7R1SiZji4oidPYe0tWsB0Bg0eDTJYIvPM+SY3fCtpeO250agOJRv01PTrrUICDWyaUkEsVGpbF12jPANZ2nTty71b/HD2T2v31VOpokTe+PY8/tp0hNzUDQKt97TsNr0NxNCCHuQZMjO2vi1IvribrZ4rWd84iOkxGVxdHs0zW6tbe/QSuWRNa9cbh7z57PbX7Z3OMVSVZXU1auJmTUba0oKKAoeDXPwbZrIVtPjxGeE4uSmZ+Ck9ujKORG6wivQhWHPteXwXxfY89tp0hKy+ev7Y/z1/TGc3B3QaBQyUnLgcrciNy9Hej/YhNo32PFaCCFEHkmG7Gxok+78dvFL0h2OENojgCO/nuOfNVE06hiAvoJuwjdr/tblxKnbUVWFGZ1ew92xcq+FZU5KIvrlV0jftAkAxwBHAludxdHTzHHnBzhyqg8o0PfhZrh5OVZobBqNQouedWjcJZAjWy8SuSuGS2fTyErNtZXxDHShaddAmnWvXWX+RoQQojKTZMjO2tVqjkZ1xKrNZp/LWXy9nEhLzObglvPc0q9sOuqWp+OXYlh8bD7ooIXrEIY2rdzNY1mHDnN+ymTMF6NBp8W3RRbeDS+i6A2ktn+TLWubAhbaDqhHUGMvu8Wpd9DSqncQrXoHkZ1hIi0hG1VVcfV0tDWZCSGEKBtVq2NKNaTT6Kjj1ByAzee30/6OvIn9wjecw2yy2DO0YlmtVsb9+iLo0tBZAvhicOVuHkv++WfOjBqF+WI0em8nQnpH4xMWj1LnFizj/+KPf9uTm20hsL6RDndUngVOHV30+NZ1w6+euyRCQghRDiQZqgRurdMJgDMZ+wlt54urp4Gs1Fwid8bYObKivbj+a5KVvaiqhjndZuNmcLJ3SIVSVZVLH3xA9MuvoObm4hqsJaTHKRw9LXDrs/DwH+zeoSM2KhWDs44+45qW2RIZQgghKj/5xLcT1aqSfTKZzPA47nXtjkZVUJ1OsfNMDK16BwGwf/1ZrDcwCV9F2HY6gt8ufgpAN+/7GdCovZ0jKpxqtRI7azbxn+TF6tM8gzodz6H1DoQHV0Pv6ZyNTGXfurMA9BrTGHfvypnUCSGEKB/SZ8gOsg7Fk7z6JJaUvE6xemCxbjaf+i/jh8MbeW/AaPb8dpqUuCyiwi9R/xY/+wZ8jSxTDk9uegFFm4uLtREfDpxq75AKpZrNRL8ynZQVKwDwb5uMV8NMaHwHDPkQnL3ISMlhw8IjADTvUZv6bSrXey2EEKL8Sc1QBcs6FE/CtxG2ROgKb7M7r1wYj3LqNHqD1jZvzL+bztkjzCKNX/kWOdoosDjyWf/56HWVL6e25uZyYerUvERIgVqdkvBqYoXb34V7vwVnL1SryoaFR8hKM+Fd25Wudzewd9hCCCHsQJKhCqRaVZJXnyx0n4KCCoxLbEP42Tia3VobRaMQfSKFhAvpFRtoERbu3UR4+k8A3Bs6lda1Kk9H4yusmZmcf+xx0tZvQNGo1OmaiLFtHRi/Gdo/bFthft8fZzh/NAmdg4Z+jzRDp5dh6kIIURNJMlSBcqJSCtQIXU2Dgp/Zk3VbNuHqaSC0lQ8Ah/66UFEhFulYfDTvhk9HUVRq67rxSs9R9g6pAEtqKmcfGkvG9u0oOitB3RNw6zcAJmwG/6a2cjGnUti1KgqAW+9thFdg1Vk6RAghRNmSZKgCWdOunwhdLTrhCKqq0qxH3izUkbtiyM02l2doxco1m3lw9VO2Waa/GzrPrvEUxhwfz5lR95L170E0eiv1bkvGZezrMGIRGNxs5bIzTKz76hCqVaVhOz+adAm0X9BCCCHsrvJ19qjGNG4lmyPmkiGCbafO0DWsHh7+ziTHZnJsdyzNu9tviY5HV79FuuYIqlXP/B7z8XJ2K/6gcmS1WrgQcZj05CRcPTzxM3pxfsx95EbHo3W0UHeQFsfHVkHdjvmOU1WVTUsiSE/MwejrRM/RjVEuN5sJIYSomSQZqkCGECNao0ORTWXxulQOuxzjy/2/0K3+VJrdWottP50gYttFuyVDX+/9g39SlqEocFfQZHrXb22XOK44vms7mxZ9QXpivG2bo9lC08wsgpzN1B1TH4fx34BrwVXcD2w+T9S/8Wh0Cv3HN8fBSf4LCCFETSfNZBVI0Sh4DK5fZJl9zZKwKir7kn4jK9dMWMcANBqFuDNpJEZnVFCk/zkQfZb3/n0VRVGppe3Om30eqvAYrnZ813ZWvTs7XyIEkK3VsK+eP6bRPXCYsqbQRCjuTCrbfz4BQNfhDfGta9/aLSGEEJWDJEMVzKm5D973N0FrzN9kpjUa8L6/CUOHDgarARxi+GD7rzi5OVC3uTcAkTujKzTW5KwMHvr9CdCmo7fU4vthb1Xo9a9ltVrYtOiLwncqCigKW0+asRbS7JWTZWbdl4ewWlRCW/vSoqf9mhyFEEJULpIM2YFTcx8CXuiAzyATXvp5+Hi+T8AL7XFq7oOno5EWxr4A/HzyOwCadM7r4Bu5M6bCZqS2Wq2M+OlpcrWnweLEZ30/wNPZtUKufT0XIg4XqBG6VlpCPBciDufbpqoqW749Smp8Nm5ejvQaI/2EhBBC/EeSITtRNAqOHdrhrPsbx6wNKOn/1fpM6zoeVVXI0h1hzdG91GvhjaOLnoyUXM5HJFZIfI+teYcY63ZUVcPUVm/SISisQq5blPTkpBsqd+jPC5zYG4dGo9BvfDMcXfTlEZ4QQogqqsokQ7NmzaJLly44Ozvj4eFRomOWL19O//798fHxQVEUwsPDyzXGUnN0B/+8Fes5u9O2uYV/KP7avLW+3v/nc7Q6DQ3b+wN5w+zL21t/LWd74jcA9PYdz7i2/cr9miXh6u5esnIenrbfY06l8PePxwHoPKw+ASHGcolNCCFE1VVlkqHc3FxGjBjB448/XuJjMjIy6Nq1K3Pnzi3HyG5S3c55P89sz7f5qXaPARBr3cHa43to1CEvGYr6Nx5zrqXcwvnp4D98c2I2iqISaujD+4OeKLdrlUpOOk5fPINjrhnU6zcVunn7ULtJMwAyU3NZ+0VeP6H6t/jaFsAVQgghrlZlkqHXXnuNqVOn0qJFixIfM2bMGGbMmEGfPn3KMbKbFNwt7+epzfk2D27SHh86ATBr+3z8Q9xx9TJgyrFw9kj5NJX9dSqS13Y/haLNwZ0wlt39VqXoW6OmRnPpsR5cWh9N04vxtuU0CtPrwQloNFqsFit/LDhERnIOngHO3PZAk0rxWoQQQlQ+VSYZqig5OTmkpqbme5Sr0B6gaCHhBCSdybfr9e7PoqpakjnM/x3cQIPLq9ef2BtX5mHsOXeGSZseA10qDtZa/HL3lzjqSjZJZHlSY48Q90hv4ndkAtD8zv4Mfnoarl4++cq5efsw5OmXaNixCwC7Vp3iQmQyeoOWAY+2wMFR5hMSQghROLlDXGPOnDm89tprFXdBRyPUaQ/ndsLJjdBunG3XrSFh1N3Wm3OWP5i/dy6/dPs/wjecI+pAXlOZzqFsFhY9cDGGcWsfRXWIR2vxZtmdX+Pn4ln8geVMPb6J6KcnkHI8r8Oz/5MT8Hp8Kr5Ag/ad8s1AXbtJMzSavPfj1P5L7Ft3FoDbHmgi644JIYQokl1rhmbOnImiKEU+9uzZU6ExTZs2jZSUFNvj3Llz5X/RBr3zfp7YWGDXe/1fQDW7Y9LEMfPwR7h5OWLOsXDmcEKZXHrf+XPc/+tYVIcLaKxuLBn0FQ287T8Hj/rPEi488XBeIqRA4MyX8Hp8qm2/RqMlqFlLmnTtQVCzlrZEKDk2kw2LjwDQqk8QDdr62SV+IYQQVYdda4YmTZrEyJEjiywTHBxcMcFcZjAYMBgMFXpN6veGzbMg6i+wmED739DvMF8/RgQ/xU/nX+ef5BX0Ce0LiXlNZfXb3NyN/o9jETzz1yRwiENjdeOLvp/TMqDBzb6am6OqqBtncWHuV6SddwKtQu35b+M+8PZiDzXlWPj984OYsi0ENjDSeWjRs30LIYQQYOdkyMfHBx8fn+ILVne1WoOTJ2Qlwfk9UK9zvt0zbrubDQvXkazdxcKsjxnKBE4fiMeUa0F/g01lH25bx+eRr6LoM9BaPVk04Cv+v707D4+iytc4/q3OnpCFbJCQkAQIgkA07DBuuICoCI6oqCiIozMoipfLHWZGB8ELwr2O+4I7IOI2ShQZN1wQYQQTMOwkYU2ALIRshOzddf9oac1lSyBJd9Lv53nymOqqOvXrtI/9eurUORdGdW+CN3MO6mowU6Zw6LWvOHrAD8PTQswLL9HuskvPeKppmnz39k6KDh3DP9ibEff0xsNDQ+JEROTMWs23RXZ2Nunp6WRnZ2O1WklPTyc9PZ3y8nLHMT169CAlJcWxXVRURHp6Otu322+bZGRkkJ6eTl5e88/V0ygWD+h6uf333SfeKjMMg1evnQO14eQHbuOoTyl1NTb2b2n8rbLK2hpu+2Aur2T9GcPjGP5mHClj3nF+EKosxnzrBg69/iVl2X7gYaHT8y80KAgBbFl1gKzUfCwW+wKsAcEt3LsnIiKtVqsJQzNnziQ5OZlHH32U8vJykpOTSU5OrjemKCMjg9LSUsf28uXLSU5O5tpr7bdYxo0bR3JyMi+//HKL139Gib9MbLjzXyfd3bNDRx4b/CSmzZddYakAZKQ1bq2yf+3cwEVvjWVL5XsYho14n4v4+tb3SWgfc06ln7Pi/ZivDyf3n1sp2+8PHhZinn2WwGHDGnR67q4S1v7TvgDr0Bu7Ed0tpBmLFRGRtsYwzdPMYCeUlZURHBxMaWkpQQ2cAfmsVBbDE4lgq4X7UyHi5D01L/y4gmVpr3LjtoeosdTwu+m9GNjl9JMJpufu4eHvnmV/zXcYhgk2X8Z1eZC/XTLe+XPvHNyIufQW8lZVU7InACwWOj31JEFXX92g04+VVvPB46lUlNbQrX8kw+/u5fz3JCIiTteY7289Wu8q/NpDl8tg10rY/glc+l8nPWzKkOvw9oSSzFICaoOZmfIIIVEDuW/AWAbEhePn7YFpmuQUH+WjHd/y6e4VFNhSMQwbhgHhxkBeuGYmvTrEtez7O5mdn2F+eDf567wcQSj6f/6nwUHIarXx5WtbqSitITQ6gGHjtQCriIg0nsKQK+k15pcw9PEpwxDAvQOu490NaylKr6ZLWVd+CH+ZyauXYKuKxWL6YHgVgXcehqUWsE/Y7G/twQN972P8hQ279dTs1r+C+dkM8n8OpHhXABgGUY/PJXjUdQ1u4sdlu8ndVYqXrwcjNbGiiIicJX17uJLzrgGLJ+RvhcJdEH7qx9yHXtSDFemb6Fk2iPV8SY1nGZZ2O+sdY7EF09V/MJMuuIVrz+vvGr0mNit89Qjmjy9RsCmI4sx2AET992OEjBnT4GayUvPZ9I19DqgrJ55PSAf/5qhWRETcgMKQE5mmWT+g+IdCwiWw+1t779Al0095bsx57fHy8aC2ElZcsoJDfnvYVLCTo9WVhPt2IDnqPHqEdcNiuNAY+ZoKWHYP5o4VHN4SSNFOexDqOOtRQsaObXAzRw6W8+2SHQD0vTqOLhdGNEu5IiLiHlzom9K9VPz8M/tuGUf5mrX1d/S6wf7PTe+ednV2Dy8LnXuFArB/cxH9OvZjUtLtTB3wB27vM4rzw7u7VhAqPwyLr4OdKyjcHsyR7YEAdHj4YdqfYeLN36qprOPzV7ZQV2Mjpkd7Bl3fpbkqFhERN+FC35bu5egXX1K1eTP58+dh1tX9uqPXDeAVYF+4df/aUzcAJFxg7xHZu6mwOUs9d4cz4fUr4OAGCjMjKNxiXysscsYMQu8Y36imVr+XSWlBJe1CfRj+h15YLC5w609ERFo1hSEnCb9vMh4hIdTs2k3x0qW/7vAJhD432n/fsPi0bcT1DsOwGBQdOkbp4cpmrPYc7Pke3rgSs3g/hXtiObzRvtRIxLRphN01sVFNZaXmk7E+D8OA4Xf3xq+ddzMULCIi7kZhyEk8goOJeOghAAqeeprqXbt+3dl3ov2f2z+xzz90Cr4BXkQnBgOwd9PhZqr0HGxcAm//HrOylII953H4JysAEVMfJPzeexrV1NGiKla9kwFAv2viieoa3OTlioiIe1IYcqKQm28i4KKLMKurOTj9v7DV1Nh3dOoLHXqDtRo2vX/aNhKSXPBWmc0GX8+C5VMw6+rI25tMUepRACL/MoPwyZMb2ZzJ1wu3U1NZR4eEIAZcE9/0NYuIiNtSGHIiw2Ihet7jeLRvT/XOneT/9xxM07RPDNRvov2g9QvAWnfKNhIusC90m7u7lKpjtS1Q9RnUVsKHd8Gap7HWGhzYPoCSn/LBMOj42GzCJk5sdJNbvjvAoawSvHw8uGrS+Vi0AKuIiDQhfas4mWdEBNHz54FhUPLPf1K4YIF9x4W3gX8YFO+DbSmnPD8o3I/Q6ABMm0n2tsYv3Nqkygtg0XWw/WNqKnzY/9MFlG89iOHtTaennqT9zTc3usmywkrWLd8D2NcdC47QfEIiItK0FIZcQLtLL6XD3x8BoPC55yl+913wDoDBv9xOWvOU/dbTKcT3sfcO7dvsxFtleVvhtSvgYBrlR0LZ930c1TkFeEZEEPf2EoJGjmx0k6Zp8v27GdRVW4nqFkyvi6KboXAREXF3CkMuIvS22wj70x8ByJv9GEVL3oYB94B3IBRsh8wvTnlufJI9DO3fVoTVeurQ1Gy2LoM3rsIszqYgI5aclb5YS8vx7dWL+A//iV9S0lk1u3vjYbK3FWHxNOzrjukxehERaQYKQy4kYupUQu+eBED+3LkceWcZDLjbvvPb/z7l2KEOCUH4BXpRU1lH7q7SlirXvrTGykfhw7uoKapm/9ouHPnZ/sRYyK3jiFv6Nl4dOpxV07U1VtZ+lAVAvxFxtO8Y0GRli4iI/JbCkAsxDIPI6dMJv89+e6zgiSco3BFiX9G+YDtsWHjS8ywWg7jeYQDsa6mnyiqKYOlNmGueoSgjgD1fRVN5sApLYCCdnnmGqEcfxeLre9bNp6/MpryomnahPiSPiGvCwkVEROpTGHIxhmEQ8eCDREx9EIDDL71KQeGl9pU5vp1jDyEncfxW2d4thfYn0ppT7iZ47XKqNnzP/m8jyP85GLPWiv+gQSSkpBB09Yhzar68uIqNX+wHYOjvu+Hl7dEUVYuIiJyUwpCLCp88mcj/si/UeuTTn8jbGo9ZUQIr/37S42N7hmLxNCg7XElxXsU5X99ms5KzbTM71n5PzrbN2GxW+1pp617G+uJV5H1dxN4vI6k87IXF35+Osx6l88I38Y7pdM7XTl2xl7paG1HdgunWL/Kc2xMRETkdrVrvwsLuvhtLUBB5j86iZFsNdSWhdDLfxtLtKug1pt6x3r6exHRvT/b2IvZtLiQ06uzH2GSt/zffLnqV8qJfb7m1ax/KsK4VhKevpXBLe6w19hwdeNVVRM6Y0SQhCKAkv4IdP+YB9l4hw9CgaRERaV7qGXJx7W+6iZjnnsXw9qb8oC85q8KwfvgglB444djjt8r2bTn7cUNZ6//N8qcerxeEAMqLj/BpaiWbdnXCWmPBJzGRzosWEvP8c00WhAB++nQPps0kvk8YHbtoyQ0REWl+CkOtQOCVVxL7+mtY2rWj4rAP+/7lTfWLN0FV/SfHjoehvN2lVJbXNPo6NpuVbxe9eoq99h6a7TERRPz9YRJSlhEweHCjr3E6hQfKyUorAGDQ6C5N2raIiMipKAy1EgEDBxK35C08I8OpKfNi37tFHH1sDNRWOY4JDPUlLKYdpgn7tzZ+NuqDO7ad0CNUj2FQ5elB5YVJGJ5Nf4d145f2QdPd+kUSHhPY5O2LiIicjMJQK+LbsycJHy3Dr/d52GotHFiWx6HbL8FakOM4JiHp7GejLsnMaNBx5SXFjW77TEoPV7IrLR+AvlfrUXoREWk5CkOtjGdEBHHvfEDY2BGASenmo+weMZzi157BrKtzLM2Rvb0Ia92ZZ6O2lpZS8tEy9t10A6VzHm9QDe1C2p/LWzip9JXZmCZ07hVKRKx6hUREpOXoabJWyPD2JnLOM7S7+ANyZz5KTSnkPfkKhW+8RfCtd+AXkEzlMSuHMkuIPT+03rmm1Up1VhYVqWkc/eZrKlJT4ZclPEINEz9rLZUeXqe8dmBYOJ169mrS91NRVsOOf+cC0FcTLIqISAtTGGrF/EfcTJeBv6N49gQKv8umrqSSIwteJaT7bVRG/45NT7yNEbQHDAOzqoraQ4eoyd6PWVVdrx2f4FqC4qoIvvoyvGJGsfyNJae85rAJ92KxNO0kiJu+zcFaZ6NDQhDRiSFN2raIiMiZKAy1ckb7ToQ+8zUhmz/h6GuzKN1cQkTxZnKjf0e+GUXZlwv4/zP1WDxt+IXVEhBVRWA3P7wvHgf9J0FYVxKB64NjT5hnKDAsnGET7iVx0NAmrb+6so6tq+zTBPQdEad5hUREpMUpDLURlqTRBD93PcH719Lhp/fY9m0tVb5h+PYPJMSWg+EBXn5WvIIteCd0w+hyCSReCXEXgad3vbYSBw2l64BB9qfLSoppF9KeTj17NXmPEMCOtYeoqbLSvqO/Y/C3iIhIS1IYaksMA+Ivwif+ImIOb2L/liMUjXmJLgOt4OULvsEQFAMeZ/7YLRYPYnslNWu5NpvJll96hS68sjOGRb1CIiLS8vQ0WRvleMR+N9B5EERdAO3jGxSEWsr+LYWUFVbhE+BJ4sAOzi5HRETclMJQG3X8Efv8fWVUlDV+NuqWsPk7e69Qr4uitTK9iIg4jcJQGxUQ4kNE50Awz22tsuZy5GA5B3YWYxjQ+9IYZ5cjIiJuTGGoDYs/h9mom9vmX8YKdbkwgsBQXydXIyIi7kxhqA07Pm4oZ0cRdbVWJ1fzq6pjtWSuywMg6XL1ComIiHMpDLVh4bHtCAjxoa7GxoGdTb+e2Nna+WMudbU2wmLaEdUtxNnliIiIm1MYasMMw/j1VtmWxq9i3xxM02TbD4cA6HNpJ02yKCIiTqcw1MbF9wkD7I+xm6bp5GrgUFYJJfkVePp4kDhAj9OLiIjzKQy1cTE92uPpbaG8uJrCnHJnl+PoFeo+oAPevq4z55GIiLgvhaE2ztPLg9ie9pXr9zr5qbLK8hp2/1wAQK+Lo51ai4iIyHEKQ24g4QL7uKE9Px92ah0Z6/Kw1ZlEdA4kMi7IqbWIiIgcpzDkBhIuiMBiMThysJyi3GNOqeG3A6fVKyQiIq5EYcgN+AZ4EXu+/VbZrrR8p9RwfOC0lwZOi4iIi2k1YWju3LkMHToUf39/QkJCznh8bW0tM2bMoE+fPgQEBBAdHc2dd97JoUOHmr9YF5TYPxKAXRsKnPJU2fFeocSBGjgtIiKupdWEoZqaGm666SYmT57coOMrKirYuHEjf//739m4cSPLli0jMzOT66+/vpkrdU0JF0Tg4WmhOK+CIwdb9lZZvYHTF+kWmYiIuJZW87/os2fPBmDRokUNOj44OJiVK1fWe+35559n4MCBZGdn07lz56Yu0aV5+3nSuVcoezcVsistn/CYdi12bQ2cFhERV9ZqeoaaQmlpKYZhnPY2W3V1NWVlZfV+2orE/vaxOllp+S12q0wDp0VExNW5TRiqqqriL3/5C7fddhtBQafunZg3bx7BwcGOn9jY2BassnnFJ4Xj6W2hrLCKw9lHW+SaGjgtIiKuzqlhaNasWRiGcdqftLS0c75ObW0t48aNw2az8dJLL5322L/+9a+UlpY6fnJycs75+q7Cy8eD+D72OYeyUlvmqTINnBYREVfn1G+nKVOmMG7cuNMeEx8ff07XqK2t5eabb2bv3r18++23p+0VAvDx8cHHx+ecrunKEgd0YNeGAjJ+ymfwDV3x8Gi+PPzbgdO9L+7UbNcRERE5F04NQ+Hh4YSHhzdb+8eDUFZWFt999x1hYWHNdq3WIq53GH6BXlSW1ZC9rYiEpOb7++/8968DpyM6BzbbdURERM5FqxkzlJ2dTXp6OtnZ2VitVtLT00lPT6e8/NfFR3v06EFKSgoAdXV1jB07lrS0NJYuXYrVaiUvL4+8vDxqamqc9TaczsPTQvdBHQHYsbb55lwyTZNtaw4CGjgtIiKurdUM4pg5cyaLFy92bCcnJwPw3XffcdlllwGQkZFBaWkpAAcOHGD58uUAXHjhhfXa+u057qjn0Cg2fZ3D/i1HqCirwT/Iu8mvcTCzhNKCSrx8NXBaRERcW6sJQ4sWLTrjHEO/fVw8Pj7eKTMttwZh0e3okBBE/t4yMtblkTy86edc2vaDvVeo+8COGjgtIiIurdXcJpOm1XNoFABbfziIzda0obGirIY9Px8GNOO0iIi4PoUhN9V9YEd8/D0pO1zJ/q1HmrTt7WsOYbOaRMZp4LSIiLg+hSE35eXjwfm/s/fabPqm6eZSstba2LLqAABJl7edCStFRKTtUhhyY32GxWBYDA5mFHPkYPmZT2iArLR8KspqCAjxoVu/yCZpU0REpDkpDLmxwFBfulxon2co/evsc27PNE1HO0nDYvDw1L9eIiLi+vRt5eYuvMr+JFnG+nxK8ivOqa2c7UUcOXgMTx8PztfAaRERaSUUhtxcx4Rg4vqEYdpMUv+196zbMU2T9cv3ANDrd9H4Bng1VYkiIiLNSmFIGDSqCwCZqfkcOXR2Y4f2/HyYgv1H8fTxoO/VcU1ZnoiISLNSGBIiOgfSJTkCTPhx2e5GT1Zps9ocvUIXXhHbLDNai4iINBeFIQFg8OguWDwN9m89wq4NBY06d/uaQxTnVeAb4OUYgyQiItJaKAwJAO07BtB/ZDwAP7yfSVV5bYPOKy+u5seU3QAMuC4BHz8tvSEiIq2LwpA49B0RR2h0AJVHa1n1TsYZb5fZbCZfL9pGTZWVyPggel/aqYUqFRERaToKQ+Lg4Wlh2B09sFgMdm8sIO2zfac9/t/LdnEwowRPHw+unNgTi8VomUJFRESakMKQ1NMxIZiLb0kE4KdP95L6r70n9BCZpsm6T3az6Wv7Mh7Dxp9H+44BLV6riIhIU9AADzlB70tjqCyv5adP9/LTp3s5mFlsv4UWFUBR7jE2frGfg5klAAz9fTe6D+jo3IJFRETOgcKQnNSAaxPwD/Lmh/ezOJhRwsGMknr7Pb0sXHJrd3oO1UzTIiLSuikMySn1urgTsT1D2fhVNvs2F1JRWk1AiA9xfcJJvqozwRF+zi5RRETknBlmY2fYczNlZWUEBwdTWlpKUFCQs8txKtM0MQwNkhYREdfXmO9vDaCWBlMQEhGRtkhhSERERNyawpCIiIi4NYUhERERcWsKQyIiIuLWFIZERETErSkMiYiIiFtTGBIRERG3pjAkIiIibk1hSERERNyawpCIiIi4NYUhERERcWsKQyIiIuLWFIZERETErXk6uwBXZ5omAGVlZU6uRERERBrq+Pf28e/x01EYOoOjR48CEBsb6+RKREREpLGOHj1KcHDwaY8xzIZEJjdms9k4dOgQgYGBGIbh7HJcUllZGbGxseTk5BAUFOTsctyePg/Xos/DtejzcC3N+XmYpsnRo0eJjo7GYjn9qCD1DJ2BxWIhJibG2WW0CkFBQfqPiwvR5+Fa9Hm4Fn0erqW5Po8z9QgdpwHUIiIi4tYUhkRERMStKQzJOfPx8eHRRx/Fx8fH2aUI+jxcjT4P16LPw7W4yuehAdQiIiLi1tQzJCIiIm5NYUhERETcmsKQiIiIuDWFIREREXFrCkNy1lavXs2oUaOIjo7GMAw+/vhjZ5fktubNm8eAAQMIDAwkMjKSMWPGkJGR4eyy3NqCBQtISkpyTCY3ZMgQPv/8c2eXJb+YN28ehmHw0EMPObsUtzRr1iwMw6j307FjR6fVozAkZ+3YsWNccMEFvPDCC84uxe19//333H///axbt46VK1dSV1fH8OHDOXbsmLNLc1sxMTHMnz+ftLQ00tLSuPzyyxk9ejTbtm1zdmluLzU1lVdffZWkpCRnl+LWevXqRW5uruNny5YtTqtFy3HIWRs5ciQjR450dhkCfPHFF/W2Fy5cSGRkJBs2bOCSSy5xUlXubdSoUfW2586dy4IFC1i3bh29evVyUlVSXl7O7bffzmuvvcacOXOcXY5b8/T0dGpv0G+pZ0ikDSotLQUgNDTUyZUIgNVq5b333uPYsWMMGTLE2eW4tfvvv59rr72WK6+80tmluL2srCyio6NJSEhg3Lhx7Nmzx2m1qGdIpI0xTZNp06Zx0UUX0bt3b2eX49a2bNnCkCFDqKqqol27dqSkpHD++ec7uyy39d5777Fx40ZSU1OdXYrbGzRoEG+99Rbdu3cnPz+fOXPmMHToULZt20ZYWFiL16MwJNLGTJkyhc2bN7NmzRpnl+L2zjvvPNLT0ykpKeGjjz5iwoQJfP/99wpETpCTk8PUqVP56quv8PX1dXY5bu+3Qyz69OnDkCFD6Nq1K4sXL2batGktXo/CkEgb8sADD7B8+XJWr15NTEyMs8txe97e3nTr1g2A/v37k5qayrPPPssrr7zi5Mrcz4YNGygoKKBfv36O16xWK6tXr+aFF16guroaDw8PJ1bo3gICAujTpw9ZWVlOub7CkEgbYJomDzzwACkpKaxatYqEhARnlyQnYZom1dXVzi7DLV1xxRUnPK1011130aNHD2bMmKEg5GTV1dXs2LGDiy++2CnXVxiSs1ZeXs6uXbsc23v37iU9PZ3Q0FA6d+7sxMrcz/33388777zDJ598QmBgIHl5eQAEBwfj5+fn5Orc09/+9jdGjhxJbGwsR48e5b333mPVqlUnPPknLSMwMPCEMXQBAQGEhYVpbJ0TTJ8+nVGjRtG5c2cKCgqYM2cOZWVlTJgwwSn1KAzJWUtLS2PYsGGO7eP3eSdMmMCiRYucVJV7WrBgAQCXXXZZvdcXLlzIxIkTW74gIT8/nzvuuIPc3FyCg4NJSkriiy++4KqrrnJ2aSJOd+DAAW699VYKCwuJiIhg8ODBrFu3jri4OKfUY5imaTrlyiIiIiIuQPMMiYiIiFtTGBIRERG3pjAkIiIibk1hSERERNyawpCIiIi4NYUhERERcWsKQyIiIuLWFIZERETErSkMiYi0cfn5+VxxxRXExMRokViRk9AM1CIibdyMGTOora3lxhtv5MYbbyQrK4vAwEBnlyXiMtQzJCIt4rLLLuOhhx5yyrUnTpzImDFjTrl/0aJFhISEtFg9ZxIfH88zzzzT6POOHDlCZGQk+/btq/d6cHAw0dHRdO/enfDwcHx9fevtHzt2LE899dQ5VCzSuikMibiZUwWDVatWYRgGJSUlpzz3dKHBMAw+/vjjJqnRXTR1CJs3bx6jRo0iPj6+3uvTpk1j/vz5dOzYkbvuugsvL696+2fOnMncuXMpKytrslpEWhOFIRGRNqCyspI33niDP/zhDyfsy8zMpKKigilTprB27doT9iclJREfH8/SpUtbolQRl6MwJCItpq6ujilTphASEkJYWBiPPPIIvx22WFxczJ133kn79u3x9/dn5MiRZGVlOfYf70n58ssv6dmzJ+3atePqq68mNzfXcYzVamXatGmOa/z5z3/mbIZGfvrpp/Tr1w9fX1+6dOnC7Nmzqaurc+w3DIPXX3+dG264AX9/fxITE1m+fHm9NpYvX05iYiJ+fn4MGzaMxYsXO3rfVq1axV133UVpaSmGYWAYBrNmzXKcW1FRwaRJkwgMDKRz5868+uqrp633888/x9PTkyFDhpyw780332T06NFMnjyZFStWcPjw4ROOuf7663n33Xcb+VcSaRsUhkSkxSxevBhPT0/Wr1/Pc889x9NPP83rr7/u2D9x4kTS0tJYvnw5P/74I6Zpcs0111BbW+s4pqKign/84x8sWbKE1atXk52dzfTp0x37n3zySd58803eeOMN1qxZQ1FRESkpKY2q88svv2T8+PE8+OCDbN++nVdeeYVFixYxd+7cesfNnj2bm2++mc2bN3PNNddw++23U1RUBMC+ffsYO3YsY8aMIT09nT/+8Y88/PDDjnOHDh3KM888Q1BQELm5ueTm5p7wPvr378/PP//Mfffdx+TJk9m5c+cpa169ejX9+/c/4fXa2lqWLl3K+PHj6dGjB0lJSSftARo4cCA//fQT1dXVjfpbibQJpoi4lQkTJpgeHh5mQEBAvR9fX18TMIuLi0957sKFC03ghHMDAgJMwExJSTnluZdeeqnZs2dP02azOV6bMWOG2bNnT9M0TTMzM9MEzLVr1zr2FxYWmn5+fuYHH3xQ7/q7du1yHPPiiy+aHTp0cGxHRUWZ8+fPd2zX1taaMTEx5ujRo0/7voKDgx3bF198sfn444/XO2bJkiVmVFSUYxswH3nkEcd2eXm5aRiG+fnnnzveW+/eveu18fDDD9f7G///6x4XFxdnjh8/3rFts9nMyMhIc8GCBad8D6NHjzYnTZp0wusfffSRGRERYdbW1pqmaZpPP/20mZSUdMJxmzZtMgFz3759p7yGSFvl6bwYJiLOMmzYMBYsWFDvtfXr1zN+/HjHdrt27Ry/jx8/npdffhmAwMBANm7ceEKbiYmJZ7zu4MGDMQzDsT1kyBCefPJJrFYrO3bswNPTk0GDBjn2h4WFcd5557Fjxw7Ha/7+/nTt2tWxHRUVRUFBAQClpaXk5ubWu1Xk6elJ//79G3WrbMOGDaSmptbrCbJarVRVVVFRUYG/vz9gH2tzXEBAAIGBgY5aMjIyGDBgQL12Bw4c2OAaftu2YRh07NjR0fbJVFZWnvCUGMDChQu55ZZb8PS0/+f+1ltvZfr06WzcuJG+ffs6jvPz8wPsPW8i7kZhSMQNBQQE0K1bt3qvHThwoN52enq64/egoCDH7xaL5YRzm8KpwoppmvUC1P9/EsowjLMaE3Q6NpuN2bNn8/vf//6Efb8NHCerxWazASfWffy1hjpd2ycTHh5OcXFxvdfy8vL44osvsNls9cKv1Wpl4cKF9cLQ8dt7ERERDa5RpK1QGBKRk2qOwLNu3boTthMTE/Hw8OD888+nrq6O9evXM3ToUMA+b05mZiY9e/ZsUPvBwcFERUWxbt06LrnkEsA+aHvDhg31vvjPpG/fvmRkZJzT36BHjx589tln9V5LS0urt+3t7Y3Vaj3ra/xWcnIyb7/9dr3XlixZQlxc3AlTHnzzzTc89thj/OMf/8DHxweArVu3EhMTQ3h4eJPUI9KaaAC1iLSYnJwcpk2bRkZGBu+++y7PP/88U6dOBey32UaPHs0999zDmjVr2LRpE+PHj6dTp06MHj26wdeYOnUq8+fPJyUlhZ07d3Lfffeddu6kk5k5cyZvvfUWs2bNYtu2bezYsYP333+fRx55pMFt/PGPf2Tnzp3MmDGDzMxMPvjgAxYtWgTg6DGKj4+nvLycb775hsLCwnO6RTVixAi2bdtWr3do4cKFjB07lt69e9f7mTRpEseOHav39NsPP/zA8OHDz/r6Iq2ZwpCItJg777yTyspKBg4cyP33388DDzzAvffe69i/cOFC+vXrx3XXXceQIUMwTZPPPvvshFtGp/Of//mf3HnnnUycOJEhQ4YQGBjIDTfc0Kg6R4wYwYoVK1i5ciUDBgxg8ODBPPXUU8TFxTW4jYSEBD788EOWLVtGUlISCxYscDxNdrw3ZujQofzpT3/illtuISIigv/93/9tVJ2/1adPH/r3788HH3wA2MeA7dixgxtvvPGEYwMDAxk+fDhvvvkmAFVVVaSkpHDPPfec9fVFWjOtTSYi0kLmzp3Lyy+/TE5OTrO0/9lnnzF9+nS2bt2KxdLw/9d98cUX+eSTT/jqq6+apS4RV6cxQyIizeSll15iwIABhIWFsXbtWp544gmmTJnSbNe75ppryMrK4uDBg8TGxjb4PC8vL55//vlmq0vE1alnSESkmfzHf/wH77//PkVFRXTu3Jk77riDv/71r47H3EXENSgMiYiIiFvTAGoRERFxawpDIiIi4tYUhkRERMStKQyJiIiIW1MYEhEREbemMCQiIiJuTWFIRERE3JrCkIiIiLi1/wNy2eLT3dd+JgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plot the energies\n",
    "plt.plot(xyz[:, 1, 2], molDB.get_properties('HF_energy'), label='HF/STO-3G')\n",
    "plt.plot(xyz[:, 1, 2], molDB.get_properties('CISD_energy'), label='CISD/cc-pVDZ')\n",
    "plt.plot(xyz[:, 1, 2], molDB.get_properties('ANI_HF_energy'), label='ANI-HF')\n",
    "plt.plot(xyz[:, 1, 2], molDB.get_properties('ANI_TL_energy'), label='ANI-TL')\n",
    "plt.plot(xyz[:, 1, 2], molDB.get_properties('ANI_CISD_energy'), label='ANI-CISD')\n",
    "plt.plot(sub.xyz_coordinates[:, 1, 2], sub.get_properties('CISD_energy'), 'o', label='TL subtraining')\n",
    "plt.plot(val.xyz_coordinates[:, 1, 2], val.get_properties('CISD_energy'), 'o', label='TL validation')\n",
    "plt.legend()\n",
    "plt.xlabel('H-H bond length (Å)')\n",
    "plt.ylabel('energy (hartree)')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "base",
   "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.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
