{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Synthetic 1D LECA Introduction" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Import LECA & Standard Python Modules" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from LECA import prep, fit, analyze # LECA Modules\n", "import pandas as pd # Pandas DataFrames\n", "import numpy as np # Numpy for standard math operations mostly\n", "random_state=0 # (int) random_state for reproduceability, set to None if this is not desired" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Import Interactive widgets, plots and display" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Magic command for clean inline graphics (inline) or interactive (notebook)\n", "%matplotlib inline\n", "## Imports for interactive widgets\n", "from ipywidgets import interact, widgets, interact_manual\n", "import matplotlib.pyplot as plt\n", "from IPython.display import display\n", "import itertools\n", "from uncertainties import unumpy, ufloat\n", "from importlib import reload" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Define Synthetic experimental setup" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# True underlying function f(x) -> y\n", "def f(x):\n", " \"\"\"1D Polynomial function for our toy data.\"\"\"\n", " return np.ravel(np.exp(-(x - 2)**2) + np.exp(-(x - 6)**2/10) + 1/ (x**2 + 1))\n", "\n", "\n", "# Synthetic experimental setup with normal measurement error sigma\n", "def run_experiment(x, sigma=0.1, n_resample=3): ## take dataframe -> output dataframe of experimental results\n", " result = pd.DataFrame(columns=['x','y'])\n", " for _ in range(n_resample):\n", " noisy_y = [y + np.random.normal(0, sigma) for y in f(x)]\n", " single_sample = pd.concat([x, pd.Series(noisy_y, name='y')],axis=1)\n", " result = pd.concat([result, single_sample],axis=0)\n", " return result.reset_index(drop=True)\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Initialize DataFrame of experimental measurements" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Combined 45 datapoints down to 15\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkoAAAGGCAYAAACE4a7LAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB2GUlEQVR4nO3deXiTZbr48e+bpEm6pDtdoC1LwbbstOAGbqCi4AI6OuKMiuP5jTOO4D7ijstRh3GcA66jx310HEXhiMuoIIjgihQBaVlbutI2XdJ0SdIk7++P2khtU7qkTdLen+vqdZG8253QvrnzLPejqKqqIoQQQgghOtD4OwAhhBBCiEAliZIQQgghhBeSKAkhhBBCeCGJkhBCCCGEF5IoCSGEEEJ4IYmSEEIIIYQXkigJIYQQQnghiZIQQgghhBc6fwcgWrndbsrKyjCZTCiK4u9whBBCiEFLVVWsVivDhw9Ho+m6zUgSpQBRVlZGamqqv8MQQgghhozi4mJSUlK63EcSpQBhMpmA1v+0yMhIP0cjhBBCDF719fWkpqZ6Pnu7IolSgGjrbouMjJRESQghhBgA3RnqIoO5hRBCCCG8kERJCCGEEMILSZSEEEIIIbyQMUpCCCGCgtvtxuFw+DsMEST0ev0xp/53hyRKQgghAp7D4aCgoAC32+3vUESQ0Gg0jB49Gr1e36fzSKIkhBAioKmqSnl5OVqtltTUVJ+0EojBra2Ic3l5OWlpaX0q5CyJkhBCiIDmdDppampi+PDhhIWF+TscESSGDRtGWVkZTqeTkJCQXp9H0nIhhBABzeVyAfS5C0UMLW2/L22/P70liZIQQoigIOtgip7w1e+LJEpCCCGEEF5IoiSE8Amr1Up5eTlWq9XfoQgh+mj58uVMnTrV32Fw+umnc+ONN/o1BkmUhBB95nA4WLJkCYsWLeLrr7/2dzhCBIwjR45www03MHbsWIxGI4mJicycOZNnnnmGpqYmf4fXa5s2bUJRFOrq6gLyfL4ks96EED6RnJxMXFwc2dnZ/g5FiIBw6NAhZs6cSXR0NA8//DCTJk3CYDCwa9cunnvuOUaMGMEFF1zQ6bEtLS19mqkVKBwOR9APwpcWJSGETxgMBkwmEyaTyd+hCBEQrrvuOnQ6Hdu2bePSSy8lKyuLMWPGcOGFF/LBBx9w/vnne/ZVFIVnnnmGCy64gPDwcP77v/8bgGeeeYb09HT0ej0ZGRm89tprnmMKCwtRFIUdO3Z4nqurq0NRFDZt2gT83FKzYcMGpk+fTlhYGCeffDJ79+5tF+ujjz5KYmIiJpOJa665BpvN5vV1FRYWcsYZZwAQExODoigsXrwYaO0qu/7667nxxhuJj49n7ty5x4yzq/NBa02kP//5z8TGxpKUlMTy5cu7+1/gE5IoCSGECEoOp9vrj9Pl7va+Ld3ctyeqq6v55JNP+NOf/kR4eHin+/xyVtby5ctZuHAhu3bt4ne/+x1r1qzhhhtu4JZbbmH37t1ce+21XH311WzcuLFHsQDcdddd/O1vf2Pbtm3odDp+97vfeba99dZbLF++nIcffpht27aRnJzM008/7fVcqampvPPOOwDs3buX8vJyVq5c6dn+yiuvoNfr2bp1K88+++wxY+vO+cLDw/nmm29YsWIFDzzwAJ9++mmP34Pekq43IYQQQempjQe8bhsdH86CaSM8j5/bfJAWl9rpvikxoVwyPdXz+MWtBTQ7Otbeuems47od24EDB1BVlYyMjHbPx8fHe1pr/vSnP/GXv/zFs+3yyy/n6quv9jxetGgRixcv5rrrrgPg5ptv5uuvv+axxx7ztMB013//939z2mmnAbBs2TLmz5+PzWbDaDTyP//zP1xzzTVcc801ADz00EOsX7/ea6uSVqslNjYWgISEBKKjo9ttHzduHCtWrPA8Liws7DK2Y51v8uTJ3HfffZ5zP/nkk2zYsIGzzjqrW6+9r6RFSQjRJbdb5f92lLKvwoqqdv5BI4Tonm+//ZYdO3YwYcIE7HZ7u23Tp09v9zgvL4+ZM2e2e27mzJnk5eX1+LqTJ0/2/Ds5ORmAyspKz3VOOOGEdvufdNJJPb5Gm5ycnF4f25mjY4fW+NtiHwjSoiSEaKesrpkfiuuYOyEJjUZhT3k9h6oaOVTVyJhh4ZyekUBUaPAPMhXB709njPW6TfOLWoO/PzXd676/rEv4u5mj+xIWAGPHjkVRlA5jgcaMGQNAaGhoh2O8ddF507bm3dFfYFpaWjrd9+iB4W1dfv21wPAvX0dP4uzMLwe1K4oyoIsjS4uSEAJovYl9sb+Kt7YVk3/Eyo6SOgAyk0ycOCYOrUbhUFUjr31VyN4jUitJ+J9ep/H6o9Nqur1vSDf37Ym4uDjOOussnnzySRobG3v1+rKysti6dWu757Zu3cr48eOB1rXMAMrLyz3bjx4w3ZPrfPPNN+2eO1aZj54sD9KdOH213Eh/kBYlIQQABeZGthXWAjBheCTjkyMB0Gk1nJQex3GJEWzIr6S0tpn1eRWkxoYSppdbiBDePP3008ycOZPp06ezfPlyJk+ejEaj4bvvviM/P/+YXVS33XYbl156KdOmTePMM89k3bp1vPvuu6xfvx5obZU68cQTefTRRxk9ejSVlZXcfffdPY7zhhtuYPHixUyfPp2ZM2fy+uuv8+OPP3pavzozcuRIFEXh/fffZ968eYSGhhIREdHpvt2JsyfnG2jSoiSEwO1W2XLADEDOyBjOnpCEMUTbbp+4CAOX5KSQEGnA4XTzbUGNP0IVImikp6eTm5vLmWeeyR133MGUKVOYPn06TzzxBLfeeisPPvhgl8cvWLCAlStX8thjjzFhwgT+8Y9/8NJLL3H66ad79nnxxRdxOp3k5ORw44038tBDD/U4zl//+tfcc889/PnPfyYnJ4fDhw/zxz/+sctjRowYwf3338+yZctITEzk+uuv73L/Y8XZ0/MNJEWV0ZkBob6+nqioKCwWC5GRkf4ORwwxu0stfLqnAmOIlqtnjuqQJB2tqLqJd7aXEKbXcvXM0eh1GhwOBw8//DAAd955Z9AXmBOBxWazUVBQwOjRozEajf4ORwSJrn5vevKZK+3mQgxxDqebrw5WA3D86NgukySAtLgw5mQlMC7B1ONxG0IIEWwkURJiiKu3taDTKkSGhjAlJapbx0xOie7foIQQIkBIoiTEEBcfYeDKk0ZhaW7pMFPoWFRV5Ui9jbjQrluhhBAiWEmiJIRAq1GIDe/ZuCJVVVm7o5RCcxPnT0rop8iEEMK/ZICBEENUbaODH4rrcLl7N59DURSiw1qTqy0HzNTW1lFRUUFhYeGAFoMTQoj+JC1KQgxRP5TUkVtUR5XVzpnjE3t1jhNGx/LZV9v5ZNN/OLz9cxSXA6vVyvjx41m4cCFZWVk+jloIIQaWtCgJMQS53aqnuvaYYT1bNuFohw/uJ/ejf1F6IA8lNJK4uDji4uLIzc1l1apVvVqTSgghAokkSkIMQUU1TTQ5XITqtYyM612i5Ha7WbNmDe7mehJHpuM2mFAUDZGRkYwfPx6z2czatWulG04EFKvVSnl5OVarLMMjukcSJSGGoPwj9QBkJJrQ/nL10G4qKioiPz+fjPRR6HVa3KqCg9bZb4qikJKSQl5eHkVFRT6LW4i+cDgcLFmyhEWLFh1zLbNANWrUKP7nf/6n2/tv2rQJRVGoq6vrt5gGO0mUhBhiHE43B6taF+nMSDL1+jxWqxWbzUZERAQxYa2re9vUn4c9hoeHY7PZ5Ju7CCjJycnk5OSQnZ3dr9dRFKXLn+XLl/fqvN999x2///3vu73/ySefTHl5OVFR3auR1lttCZmiKGg0GqKiopg2bRp//vOf2y2G212KorB27VrfB9oLMphbiCHmkLkBh9NNVGgIyVG9Xw7CZDJhNBppbGwkKdJIrKYJPT+v/N3Y2IjRaMRk6n0yJoQvud1umpubsdvtVFdXExMTg0bTP+0FRycH//73v7n33nvZu3ev57mjF3xVVRWXy4VOd+yP5GHDhvUoDr1eT1JSUo+O6Yu9e/cSGRlJfX0927dvZ8WKFbzwwgts2rSJSZMmDVgcviQtSkIMMVVWOwCZySYUpXfdbgBpaWlkZmZSXFyMMUSDQXHRdjpVVSkpKSErK4u0tDRfhC1En+Tl5fHXv/6VjRs3snnzZh544AEeffTRfptwkJSU5PmJiopCURTP4/z8fEwmEx999BE5OTkYDAa2bNnCwYMHufDCC0lMTCQiIoIZM2awfv36duf9Zdeboij87//+LwsXLiQsLIxx48bx3nvvebb/suvt5ZdfJjo6mo8//pisrCwiIiI455xz2iV2TqeTpUuXEh0dTVxcHLfffjtXXXUVCxYsOObrTkhIICkpieOOO47LLruMrVu3MmzYsHaL7H733XecddZZxMfHExUVxWmnncb27dvbvUaAhQsXoiiK53F33p/+IImSEEPMKeOG8buZo5nSx2VINBoNCxcuJD4+nry8PGw2G263G4vFwp49e4iPj2fBggX99o1diO7Ky8tj1apV7Nixg7CwsICZnbls2TJPsjZ58mQaGhqYN28eGzZsIDc3l3POOYfzzz//mOP87r//fi699FJ27tzJvHnz+M1vfkNNTY3X/Zuamnjsscd47bXX2Lx5M0VFRdx6662e7X/5y194/fXXeemll9i6dSv19fW97gYLDQ3lD3/4A1u3bqWyshJo7ba/6qqr2LJlC19//TXjxo1j3rx5nm767777DoCXXnqJ8vJyz+Pevj99JXcwIYagqLAQwg1973nPyspi6dKlTJ06laamZkrM9WzbV8yUqdNYunSp1FESftc2O9NsNpOVlYXBYECjCYzZmQ888ABnnXUW6enpxMbGMmXKFK699lomTpzIuHHjePDBB0lPT2/XQtSZxYsXs2jRIsaOHcvDDz9MQ0MD3377rdf9W1paePbZZ5k+fTrZ2dlcf/31bNiwwbP9iSee4I477mDhwoVkZmby5JNPEh0d3evXmZmZCUBhYSEAs2fP5re//S2ZmZlkZWXx3HPP0dTUxOeffw783L0YHR1NUlKS53Fv35++kjFKQgwhthYXxhDfrsuWlZXFbbfdhtlcTX5zBDkzT+dXZ88gM7l/B48K0R1tszNTU1M7dDX/cnZmWxfPQJk+fXq7xw0NDSxfvpwPPviA8vJynE4nzc3Nx2wxmTx5suff4eHhREZGelpvOhMWFkZ6errncXJysmd/i8VCRUUFxx9/vGe7VqslJyen18mkqrZW/297/ysqKrj77rvZtGkTlZWVuFwumpqajvk6e/v+9JUkSkIMEXVNDl796jCj4sM5b1Iyml6WBeiMRqMhJiaakZEmYhNHsL+yURIlERDaZmeGh3deLyw8PJzS0lK/zM78ZUy33norn376KY899hhjx44lNDSUX/3qVzgcji7PExIS0u6xoihdJjWd7d+WzPSHtq7NtkT0qquuorq6mpUrVzJy5EgMBgMnnXTSMV9nb9+fvpJESYghIv+IFZdbxely+zRJOlqMphmAQnMjdqcLg863rVdC9NTRszM7S5YCaXbm1q1bWbx4MQsXLgRaW1DauqsGSlRUFImJiXz33XeceuqpALhcLrZv387UqVN7fL7m5maee+45Tj31VE8X2tatW3n66aeZN28eAMXFxZjN5nbHhYSE4HK52j3nr/dHxigJMUTsr2j9xtyX2knHYsRJbLgep1vlYGVjv11HiO46enbmL1tNAm125rhx43j33XfZsWMHP/zwA5dffrlfxk4tWbKERx55hP/7v/9j79693HDDDdTW1nZrlmxlZSVHjhxh//79vPnmm8ycOROz2cwzzzzj2WfcuHG89tpr5OXl8c033/Cb3/yG0NDQducZNWoUGzZs4MiRI9TW1nqO88f7I4mSEEOApbkFc4MDRYH0YRHHPqCXFAXG/rR2XIFZEiXhf8E0O/Pxxx8nJiaGk08+mfPPP5+5c+f2e2HMztx+++0sWrSIK6+8kpNOOomIiAjmzp2L0XjsumsZGRkMHz6cnJwcHn30Uc4880x2797N+PHjPfu88MIL1NbWkp2dzRVXXMHSpUtJSEhod56//e1vfPrpp6SmpjJt2jTAf++PovZnx6Totvr6eqKiorBYLERGRvo7HDHI7CypY0NeJSOiQ7l0RqrPz+9wOHj44YcBuPpPN/PujiMYQ7Rce+qYfuvmE0OHzWajoKCA0aNHd+vDujN5eXmsXr2a1atX43Q6mTZtGhMmTGDBggUyO/MY3G43WVlZXHrppTz44IP+Dqfbuvq96clnrv9T6G7au3cvTzzxBIsXL2bSpEnodDoUReGhhx7q8bncbjdffvkl9957L7NmzSIuLo6QkBDi4+M566yzeP31170ObHv55ZePWZr+P//5T19frhA+1da6Myq+dwvg9kSiycAwk4HjEiNokQVxRYBom515xhlncOqpp3Lvvfdy++23S5LUicOHD/P888+zb98+du3axR//+EcKCgq4/PLL/R2aXwTNYO5nnnmGlStX+uRchw4dYubMmQDExsYyffp0YmJiOHToEOvXr2f9+vW8+eabvPPOO+j1+k7PkZ6ezqxZszrdNmLECJ/EKYQvOF1uimuaABgVH9bv19NoFH574sh+v44QPaXRaAgNDUWr1RIXFxcQ3W2BSKPR8PLLL3PrrbeiqioTJ05k/fr1QzapDJpEaeLEidx6661MmzaN7OxsHn74YV577bVenUtRFGbPns1tt93GWWedhVb788yczz//nPnz5/P+++/z6KOPcu+993Z6jlmzZvHyyy/36vpCDCS3CjPHxlNusTEswuDvcITwG71ez9y5c9m0aRO7du3i9NNP93dIASk1NZWtW7f6O4yAETSJ0n/913+1e9yXbwLp6entqpAe7bTTTmPZsmXcc889vPrqq14TJSGChV6nYVpaDNMG+Lput8qRehvDTAZCtPLNXQSGnJwcMjIy2i1KK0RXgiZRGkhtI+yLi4v9HIkQwev1b4swW+0snDZiQMZGCdEdJpMpIGomieAhiVIn9u/fD7SWdffmwIED3H333VRWVhIREcHEiRO54IILiI+PH6gwhTgmq62Fw9VNjIwLw2QMOfYBPpQUacRstXO4pkkSJeETMklb9ISvfl8kUfqFpqYmVq1aBcDFF1/sdb+tW7d26MM1Go0sX76c22+/vV9jFKK7DlY1sjG/kpSYUC6Z7vuyAF1Jiw1jd6mFop8GkgvRW23jSB0OR4fChEJ407a0ydHjkHtDEqVfuO666ygoKGD48OHceeedHbYnJSVx1113ccEFFzBmzBgMBoOndMFrr73GsmXLcLlcnR57NLvdjt1u9zyur6/3+WsRovCnsgCj/dCikxrb+oFmttpptDsJN8jtRvSOTqcjLCyMqqoqQkJCZLaaOCa3201VVRVhYWHodH2798id6ygPPvggr7zyCkajkbfeeou4uLgO+5xzzjmcc8457Z6bPn06r7zyClOmTOGWW27hgQce4JprriExMdHrtR555BHuv/9+n78GIdq0HFUWYGTcwCdKYXodCZEGKuvtFNU0kZUshVRF7yiKQnJyMgUFBRw+fNjf4YggodFoSEtL69bSK12RROknjz/+OPfeey8Gg4E1a9Z46iz1xA033MAjjzyC2Wzmk08+4YorrvC67x133MHNN9/seVxfX09q6sB2jYjBrbS2GadbxWTUER/ReT0wX7Lb7TgcDqxWq+dLRlpsGJX1dg5XS6Ik+kav1zNu3Lh+XyleDB56vd4nrY+SKAFPPPEEt9xyC3q9nnfeeadDi1F3abVaxo0bh9lspqSkpMt9DQYDBoPUtBH9p6D6p2rcceF9/kZ1LN7q04yMDWdbYS3FNU2oqtrvcYjBTaPR9HoJEyF6a8gnSk899RRLly71JEnz58/v0/mqq6sBZPqp8LvCAVy2BDqvTzM82shJ6XGMjOv/iuBCCNEfhnSi9Oyzz3L99dd7kqTzzjuvT+fbvn07+/btA+D444/3RYhC9Eq9rQVLcwtajeIZVN3fOqtPo9NqOHFMx7F+QggRLAb11IEnn3ySzMxMrrzyyg7bnn/+ea677roeJUlNTU089dRTWK3WDts2b97sKScwa9YsSZSEX0UaQ7j21HQWThuBQde3qbFCCDGUKWqQVPDavn071113nefxwYMHMZvNpKSktFuEds2aNZ5CkcuXL+f+++/ntNNOY9OmTZ59duzYQXZ2NqqqkpmZyQknnOD1ukev51ZXV0dMTAwGg4Fp06aRlpaG0+lk37597N69G4BJkybx8ccfd1mssjP19fVERUVhsViIjJRBr2LwcLlVDlQ2UGZp5rRxw9BoZJySEMK/evKZGzRdb/X19XzzzTcdni8pKWk3cPro2kTe1NXVeSp25ufnk5+f73XfoxOlsLAw7rnnHrZt20Z+fj4//vgjzc3NxMTEcOaZZ3LJJZewePFi9Pr+n2EkRLBQgA35Fdhb3IxPjiQxUgbjCiGCR9C0KA120qIkfKWy3samfVWkD4sgZ2SMv8MBYG1uKQXmRk49bljAxCSEGLp68pk7qMcoCTHUWK1Wtu8rpqCijtK6Zn+H4zEipnVAeVkAxSSEEN0RNF1vQoiuORwOlixZwu6GME4663zOnpTi75A8hkf/nChJPSUhRDCRREmIQSQxKZlK13CSk5NJiw2c2kWJJgM6jUKTw0VtUwux4TKOTwgRHKTrTYhBxKk3EaI3EGMKC6hkRKfVkBjVOoi7tFa634QQwUMSJSEGEavauixOakxYwHVvpfzU/VbdeOyZqUIIESik602IQaThp0QpJWZgqnH3xOTUaCalRGEyhvg7FCGE6DZJlIQYJFRVJURxo0UlLQATpQiD3G6EEMFH7lxCDBKKojBaW4OqQoRR/rSFEMIX5G4qxCATYEOT2imqbmJ7US3xEQZmjYv3dzhCCHFMMphbiEGiweb0dwjHZHe6KDA3UmBu8HcoQgjRLZIoCTEIWG0tvPjlYfKcCbgDeFGitsKT5gYHthaXn6MRQohjk0RJiEGguKYZt9tNc10VVZUVFBYW4na7/R1WB+EGHTFhrbPeAmmJFSGE8EbGKAkxCGzZtoON777Lge1bOOhqxmq1Mn78eBYuXEhWVpa/w2tnREwYtU0WyuqaSR8W4e9whBCiS9KiJESQ27NnD6+/8A9KD+RhCjMQFxdHXFwcubm5rFq1iry8PH+H2M7w6NYK3bJArhAiGEiiJEQQc7vdvP7v1dTVVpM8aiwmQwgajYbIyEjGjx+P2Wxm7dq1AdUNlxLdugZdRb2dFlfgxCWEEJ2RREmIIFZUVMQPu/cQHZ9MZGhIu9IAiqKQkpJCXl4eRUVF/gvyFyJDW8cpDY8OpckhA7qFEIFNxigJEcSsViu19Q2EDYsl0hhC1S+2h4eHU1paitVq9Ut8nVEUhatOHhVwa9EJIURnJFESIoiZTCYSYyPR6F2e2WRHa2xsxGg0YjKZ/BCdd5IkCSGChXS9CRHE0tLSmD5lIkpjDcaQ9n/OqqpSUlJCVlYWaWlpfoqwa80OF6oawIWfhBBDniRKQgQxjUbDwoULiY+PJy8vD5vNhtvtxmKxsGfPHuLj41mwYAEaTWD9qauqyr++LeLZzw9S3ejwdzhCCOFVYN09hfARq9VKeXl5QI3N6Q+7SixEJo/iT9cvYerUqTQ3N1NdXU1NTQ3Z2dksXbo04OooQWvXW4i29fYjZQKEEIFMxiiJQcfhcLBkyRIKCwu56667OOuss/wdUr+wtbj4LL8St6pyzSnjuO2226iursZut3PTTTcxduzYgGtJOtrwKCPFNU2U1dmYnOLvaIQQonOSKIlBKTk5mbi4OLKzs/0dSr8pqW3GrarEhIUQaQzB4VCJjo4GYNSoUQGdJMHP675Ji5IQIpAF9p1UiF4yGAyYTKaAm+3lS8U1TQCMjAv3cyS9kxRlRFHA0txCo93p73CEEKJTkigJEaQOVzcCkBob5udIescYoiUuwgBAuUValYQQgUkSJSGCUL2thdqmFjSKQkpMqL/D6bXhUa3rvpXW2fwciRBCdE7GKAkRhIqqW7vdkqIMGEO0fo6m90bHh+Nyq6QGcbI3VFmtVhoaGoiIiBjwLm5/XlsMPdKiJEQQKv1pAHSwdru1GTMsgrMnJDFmWIS/QxE9YLPZuPrqqzn//PMHfNHltlmtixYt4uuvvx6w64qhS1qUhAhCZ2UlMjklijC9/AmLgZWXl8fq1avJz8/H4XDw/vvvU1xczMKFCwesZtdQmNUqAofcZYUIQhqNQnLU4OiucrtVqhrsACRGGv0cjehKXl4eq1atorKyEpPJREhICElJSeTm5lJcXDxgBU4NBoNnZqsQ/U263oQYROx2O1arNagqku8oqeONb4r4+lC1v0MRXXC73axZswaz2UxWVhYGgwGNRkNkZCTjx4/HbDYPeDecEANBEiUhgsy720v4dE8FVltLu+f1ej1z587FZDKxa9cuP0XXc8Oj2gpP2mSB3B4ayKV6ioqKyM/PJzU1FUVR2m1TFIWUlBTy8vIoKirq91iEGEiSKAkRRCzNLRyubmJPWb1nrbSj5eTkcO2115KTk+OH6HpnmMlAiFbB1uKiRhbI7baBHtRstVqx2WyEh3de4DQ8PBybzRZUrZlCdIeMURIiiByqagAgOdrYaVmAYKxGrtUoJEYaKaltptxi8xShFMc2kIOaTSYTRqORxsbGTpOlxsZGjEZj0P3+CXEs0qIkhjyny83Xh6p574cyz7IggarA3FqNO31YcC5b4k3bum+lsu5bjwzkUj1paWlkZmZSXFzcoYtUVVVKSkrIysoiLS2t32MRYiBJi5IY0srqmlmfV0F1Q2uXz8HKBsYmRHDe5OQO4zD8ze50UVLbmkiMiR9cdYdkgdzAp9FoWLhwIcXFxeTl5WGz2dDr9VgsFsrLy4mPj2fBggUBvxizED0liZIYkhxON18eNLOjuA5VhTC9lpFx4eyrsBIXoQ+4JAlaq3G73CoxYSHEhOv9HY5PJf+0QG5dUwsNdicRBrk1BaKsrCyWLl3K6tWr2bt3Lw6Hg4qKCrKzs1mwYMGA1VESYiDJ3UgMSfsqrOQW1QEwfngkp44bRqhey0lj4jDqf/5G7HC6URQ6HTg90A5WtXa7DcYq1sYQLWdkJBAbrseo8/97LbzLysrirrvuYty4cXz++eecdtppXHrppdKSJAatoPnN3rt3L0888QSLFy9m0qRJ6HQ6FEXhoYce6tN5169fz7x584iPjyc0NJTMzEzuuusuGhoaujzuwIEDLF68mJSUFAwGAykpKSxevJhDhw71KR4xMMYnR5KRZGLhtBHMnZBEqL51YHRUWAgGXeu/K602/vVtEVv2m/0ZqkdMWAgxYSGMjh9c45PaTEmNJjU2DF0AJKWiaxqNhvnz53Pvvfcyf/58SZLEoBY0LUrPPPMMK1eu9Ok5//73v3PzzTejKAqnnHIKiYmJfPHFFzz88MO88847bNmyhfj4+A7Hbd26lbPPPpumpiYmTJjArFmz2L17N6+88gqrV69m/fr1nHjiiT6NVfiWRqMwb1Jyl/s02Vunq9c0OkiLCyPdzy05J4yJ44QxcVJrSAQEk8lEaFg4R/827quw8tXBalpcbtyqilsFl1vFoNMQaQxhxuhYT6KvqmpAdnEL8UtB8zVg4sSJ3Hrrrbz++uvk5eVxxRVX9Ol8ubm53HLLLWi1Wj744AM+//xz3nrrLQ4ePMicOXPYu3cvf/jDHzoc19TUxKWXXkpTUxN33HEHu3fv5s0332T37t3ccccdNDY2cumll9LcLINSA1GBuRGHs3uVg0fFh5M9Mgag0wKP/jJYP1xUVWV/hZVNeyuxtbj8HU7Ac7vd1NXVUVFRQWFh4YBVxG5yOPmxzMK6H8r4x+ZDHK7+eaaoRoGaRgdWm5NGu4tmhwuH043V5qS0rhnXUTEW1zTz4pYCNuRVsL/Cit0p/+ciMAVNi9J//dd/tXvc16beRx55BFVVufrqqzn33HM9z4eFhfHCCy8wZswY3nnnHfLz88nMzPRsf/nllykrK+O4447r0O330EMP8c4777Bv3z5effVVrr322j7FKHyrttHBezvKCNVr+O2JI7u1oOzM9DhKapuorLfz8Y8VXDRtBBrNwCcqxTVNJEcZB3W3lKIobDlgpq6phZFx4YO2i9EX2ham3bhxI06nE6vVyvjx4/ttYVq3W+VgVQM7iusorWvm6EZNS/PPRUJTYsL4VU4Kep0GRQGtoqBRFOxON/W2FpKOWp+wqKYJS3MLO0ss7CyxoNdpyEo2MWlENMNMUktLBI7Be9ftgsPh4IMPPgDg8ssv77B95MiRzJw5E4A1a9a029b2+LLLLuuQrGk0Gn79618D8O677/o8btE3Xx2qxq2qJEYau5UkAei0Gs6dmEyIVqG4ponvi2r7OcqO6pocrP6+hOe/KMDpGtzraHnqKdVKi6w3bQvT7tixg7CwMOLi4oiLiyM3N5dVq1aRl5fn0+tZmlt46ctC3t9ZTklta5KUEGngxDFx/OaENLLTYjz7GkO0pMaGkRhpJMFkJC7CQEy4nqQoI8clmtrNZpwxOoYLpw5nalo00WEhOJxufii28M+vD/PWd8U02p0+fR1C9FbQtCj50r59+2hqam0unj59eqf7TJ8+nS+++ILc3Nx2z7c97uq4o/cTgaGy3sbeI1YUBU5O7zjurCux4XpOz0jg0z0VfHmgmrHDIgZ0ev6+itaJBQkmw6BuUQIYER3KnrJ6qafkxS8Xpq2ubl1IODIykujoaPbs2cPatWvJyMjw2QDrSKMOrQKhei2TU6KYOCKKSGNIn89r0GkZMyyCMcMiUI9TKalt5oeSOg5WNmK1OwntpPK8EP4wJBOlgoICAKKjo71WtE1NTW23L7SuddR2Y/JWfbbtuKqqKq+l/sXA23qwdeZaZpKpV836E4ZHcqCyNWEZyCFCqqryY5kFgMzkwb80xIifWpSO1NtwutyDPjHsqZ4sTDtq1Kgen19VVfZWWNlZYuGiaSPQaTUoisL5U4YTGRrSb2UyFEUhNTaM1NgwrLYW6m1OTxe3y62yMb+SqWnRxMvyNsIPhmSi1LZoY1dJTERE6wyn+vr6Dsd1dWzbcW3HetvPbrdjt9vb7Sv6R3FNE4XmJjSKwolj4np1DkVRmD85ecDrKZXUNlPX1IJep2FcwuBPlKLDQgg3aGm0uzhSbyMlJszfIQWU7ixMW1pa2quFaSvqbWzaW0lZnQ2AnaUWT7faQK6/ZzKGYDqqxWpXqYVdpRZ2l1mYlhbDjNTB/3cgAot8XfOTRx55hKioKM9PW0uU8L0dxXUATBwRSXRY77vMupskWa1WysvLfbKK+o9lrQl0RqIJ/RAoxKgoioxT6sLRC9N2pjcL07rcKl8drObNb4spq7Oh12mYOTaeySOifBV2n4yMDWNsQgSqCtsP1/LmtlKa1L53/QnRXYP/ztuJtpuIt5sN4Ck4GRkZ2eG4ro49ulDl0cf+0h133IHFYvH8FBcXdy94cUxHT5suKCjA5XajURQmpfjmxt/kcLIhr4I9ZR1bAR0OB0uWLGHRokV8/fXXfbqO3eniQGVrsjVhhPffpcGmrfutptFxjD2HHl8vTFvb6OCtbcV8/dNEh4wkE1edPIrjR8cGTLdnTLie86cM54Kpwwk3aKlpdLDPNYwjbhMut9QUE/1vSHa9tfXd19XVYbVaO/321Za4HN3PbzKZiI2NpaamhqKiIqZMmeL1uPj4+C679gwGAwaD9Lf7mrdp0+fOv4BhPuo+yCtvHcNxqKqRsQkRHVp6kpOTiYuLIzs7u0/XKa5posWlEhehJynS2KdzBZPMpEjGDIsg0jgkb09d8vXCtJv3V3HEYsMQomFOZiIZSYHbrZU+LILhUaF8vLuUL1Qod5lYn1/J+VOlNV70r8D4yjDAMjIyCAtrHfuwbdu2Tvdpe/6XH3Ztj3t6nOh/XU2bfu7Zp8jPz/fJdaakRBEVGkKD3cn3hzuWCzAYDJhMph51f3RmbIKJK08ayRkZCYO2yGRnQvVaokJDhtRr7om2hWmnTp1Kc3Mz1dXV1NTUkJ2dzdKlS3tUR+mMzATGJkTw2xNHBnSS1CZUr+XcCYmM1NaiV1zMGBlz7IOE6KMhmSjp9Xrmz58PwBtvvNFh++HDh/nyyy8BWLhwYbttbY/ffPPNDpVw3W43//73vwG46KKLfB638O6X06YNBgMoGoxhEYwfPx6z2czatWt9Ur1Yp9VwyrjWEgPfH67p14rdcREGUmNlQLNoLysri9tuu40zzjiDU089lXvvvZfbb7/9mElSg93JrhKL53GkMaR1RpsPpvsPFEVRiNU0k6WtIPaoMh21A9RV68sxiCI4DOpE6cknnyQzM5Mrr7yyw7Zly5ahKAovvfQS//nPfzzPNzU1cc011+Byubj44ovbVeUGWLx4McOHD2ffvn3cc8897bbdc8897Nu3j5SUlE6vKfpPZ9Om7ejILa7jYFVju2nTvjA2IYIR0aG0uFS2HvD9orlDfexFuaWZ934oY/2eCn+HErA0Gg3R0dEkJiYyatSoY3a3HbHYePPbItbnVbD3SPB/yB9dIL+4polXvzrM5n1VuPvxb8eXYxBF8AiaQQDbt2/nuuuu8zw+ePAgAP/4xz94//33Pc+vWbOG5OTWxU7NZjN79+4lKSmpw/mys7P529/+xs0338y8efM47bTTSEhI4IsvvqC8vJyMjAyeffbZDseFhYXx1ltvcfbZZ/Pwww/z3nvvMXHiRHbv3s3u3bsJDw/n7bffJjQ0tMOxov90Nm26WQ0hCjDoNISHG3s9bboziqJwWsYw/vVtEXnlViaOiCIlJswzkNxut1NYWMjYsWN7XPjP5VZ5aWsBw6NDOSMjgVD90Cu853KrHKxsIEyvZY46tLoe+8O+Cisf7z6C09065i0xcnCNj6yot+FWVb4/XEuD3cncCUlo+2mpIV+NQRTBI2gSpfr6er755psOz5eUlFBSUuJ5fHRtomO56aabmDRpEn/729/49ttvaWxsJC0tjTvuuIM77rjD6xiTmTNn8sMPP/Dggw+yfv163nnnHYYNG8aVV17JvffeS3p6es9foOiTo6dNh4eH41bB8dOvd3yEgcbGhh5Pmz6WxEgjk0ZEsbPEwlcHq5kYftgn628VmBuw2pyU1DZhGAIlATqTFGlEp1FocrioaXQMaB2fwWZnSR2f5VeiqjBmWDjnTEzCoBtcyff0UbGYjCF8/OMR9h6x0uJyM29S/9Q9a5uI48t7iQhsQZMonX766R2mwx7L8uXLWb58eZf7nHnmmZx55pk9jmfs2LG88sorPT5O9I+2adO5ublkZGRgIwRVhTC9FmOIhkP7S8jOzu72tOnumjk2Hp1WQ5StglWrnqKyspKwsDBCQkI8A8mLi4u7PchWVVV2FLeOIRmfHOWXBXgDgU6rYXh0KEU1TRTXNkui1AuqqvJdYa2na3hKahSnH5cwaH+nMpJMGHQa3t9ZxqGqRtbmlnLB1OGDLikUA29ofl0Vg07btOn4+Hjy8vKot7XgdrvQu23s2bOnx9Omu8sYouWUsXF8+P577QaSazQaIiMjezyQ/EBlA8U1Teg0vqv7FKzaBrEX1TT5OZLgVG6xeZKkE0bHckbG4E2S2oyKD2dhdgp6nYaS2mbe+b4Uu9Pl77BEkJNESQwabdOmJ0ycTGOTjcaaSrA19GradE+0DSRPSUmh0dH+pvzL9be60uJy8/m+KgByRsUQFRo8M5H6Q9pPiVJJbVO/DtAdrIZHh3JyehynHjeMk8fGD5lxXiOiQ7kkJ4VQvZa4CD36ACmcKYJX0HS9CdEdWVlZXHDVH/nKHIKupYGHbl3cqwHVPWG1WmlubqbY6qbeVo9D1aJXfk6Yurv+1neFNVhtTiJDQ5gxKrbf4g0WCSYDhhAN9hY3FVYbyVEyQeKX7HY7DocDq9VKXFwcqqridKuesTkn9HJtw0DWnQkTCZFGFh2fhsmgGzIJoug/kiiJQSczKZKsOC0aIrs1bbqvTCYToaGhNLTYAB0W1UgcP3cXdWf9LZdbJb+8NZE67bj4AV98NxBpNAqj4sJpsDuHfLmEzuj1eubOncumTZvYtWsXp512Gpv2VlFptXHh1BEYQwbf2Bxvlfc7mzBxdIus262SW1zHlJSogFmaRQQPSZTEoKPXaYjVDNyCqm0Dyb/7/nv0Eck4VQ21hOJ2q2g0retvHWsguVaj8NsTR7L3iJX0YREDFnugO3dikrQIdCEnJ4eMjAzCw8P5Yr+ZHcV1KAqU1jUPut+jtsr7vZkwsSG/kt2lFkpqmzhv8vB+Kx0gBidJrYXoo7aB5EkJCSh1pbgcTdhcsKvwCD/++GO3B5LrdRompURJYnAUeS+6ZjKZSE5OZnelw7OczplZiYMuSeqs8n5PJkxkJpnQaRQOVTXy8Y9HejyDWgxtkiiJQSW3qJZth2uxqwPb7dA2kPz46dnom2toqq2g9EgVsSMzu/ymW9fkYE9Zvdy4j8HW4qLR7vR3GAHpu8IavimoAVrXbps4YvDNluys8n6b7kyYSI0N4/wprS1Je49Y+WK/76vpi8FLut7EoLKjuI5qqw2bqsOgDOy04Lb1t6qrq6lsVhh18nwSh6cyfOSYDvu21kuqY+sBMy0ulQJzI/MnJw9ovMHiq4PVfFNQTc7IGE4ZN8zf4QSUPWX1bPnpQ//U4+KZmhrt34D6SWeV94/WnQkTo+LDOXtCIh/tOsL3h2uJDgthckp0P0UsBhNJlMSgYWlqoa6pBY0CEcrALJD5S23rb0VHw7xZk0mNNxEV1jqo1NLcgqWpBZNRx/q8CkpqW8dRpcaGMeunRXZFR9FhrcVDpZ5Sew6nmy0HWstJTB8VQ87IwTtT8peV93+pOxMmoHWih6WphS8PVvNZfiUmYwij4ztPvoRoI4mSGDQO1zQCkBxlpF7xf1fW1NRo9PqfVzf/sczCN4dqPI9DtAqnjBvGZBmX1KW2wpNVVjvNDteQXPuuM3qdhl/lpLK71MKssYM70f5l5f2jqWr3Jky0OX50LHXNLeyvsMpsStEtkiiJQeNwdWuLQ1psGHv9HEtntIpCqF5Ls8PFiJhQ5o5P8rQ2Ce8iDDriI/SYGxyU1DYxLnFor7GlqqonsY4N13PqcYO/O7JtwkRxcTF5eXnYbDb0ej0Wi4Xy8vIeVd5XFIUzsxLJTothmEmWxhHHJomSGBTcbpXi2tZEKTXWv4UJf1kEsM0JY+KYMSqWphYX4XqttCL1QEpsGOYGB0U1QztRsrW4eHd7KbPGxpMWF+bvcAZU24SJ1atXc+DAAaxWKzU1NWRnZ7NgwYIeVd7XapR2SVKj3YlBp5EaS6JTkiiJQaHCasPe4sYYoiXBjwuo/rII4Omnn95uu0ajEGGQP7ueSosNY0dR3ZAep+Ryq7y/s5yKehsb8iu48qRRQ64e0NETJux2OzfddFOfK++X1TXzwc5yRseHc+b4RB9GKwYLuWOLQaGuqYUQrUJqbKjfF/5sKwIYETG4atn4U0pMKBpFoa6pBUtzy5BbB09VVTbtraS4pgm9TjOkiya2TZgAfFJ53+F00+hwsqvUQlKUcVCWVxB9I4mSGBSykiMZlxCB3ekGOi86N1BMJtMxZ9+InjHotMwYHUOkMQSDbuh1j+wormNniQVFgXMmJsnYGh8aFR/OyenxbD1g5rP8SuIi9LKuoGhn6N1xxKCl02oIl26tQevk9HgmjogalGuYdaXA3Mjn+1rLAJwyLn7QVd0OBDNGxTA2IQKXW+WDneVS3FS0I4mSCHpumeIrBqnaRgcf7ipHVWHC8Eiy02L8HdKgpCgKZ09IJDZcj9Xm5MNd5XJfER6SKImg9/m+Kl75spB9Fd6r8orBwdLUwveHaympHRqDuk1GHeMSIhgRE8qcrESZKdmPDDot508Zjl6noaS2mdziOn+HJAKE9FOIoHe4upHaphY08iEy6O0oqWP74VqykiNJiRn80+N1Wg1njU/E6VaH7ODtgRQbrmdOVgKHqhqZOCLS3+GIACEtSiKoWZpbPElSSowMwBzs0oe1LjdRYG4c1F0jRdVNntenKAohUt9nwGQmRXLuxCQMuqE1Fk54J399IqgV/1RXJynKMOQG+Q5Fw6NCCdVrsbW4KK1r9nc4QOuCreXl5V0uyNoTByqtvLO9hPd+KJMlNvykrYtTVVUOVFpRVfl/GMokURJBzbOw7BDohhGtBTvH/LSI6YGqBj9HAw6HgyVLlrBo0SK+/vrrPp+vusHOxz9WAK2LAUt3m399/OMR1v1QzrbDtf4ORfiRJEoiaKmq6hnUO0K63YaM9ITW6fEHKxsC4pt+cnIyOTk5ZGdn9+k8thYX634ow+F0kxobxqnjBv8aboGubRzclweqKQuQFkwx8CRREkGrvtmJ1eZEoyhSIG4ISYsNI0SrYLU5qbLa/R0OBoOhz0VGVVXl0z0V1Da1YDLqmDcpye8V5kVrSYbMJBNuVeXjH4/gcPq3mK3wD0mURNBSUZk4IopxiRHoh2C15qEqRKthZFw4Oo2CucHh73B84vvDtRyobECrUThv8nDC9DIhORAoisIZmQmYjDrqmlr44oDZ3yEJP5C/RhG0osP0nCWLWA5Jp2UMY+6EpEGRIDc5nHxTUAPAGRkJJEUZ/RyROJoxRMvcCUm8s72EH8us1LmNRGts/g5LDKDgv8sIIYacSGPIoEiSAML0Oi7JSeGE0bFSuydApcaGkTOytSp6iTsatyrdokOJtCiJoGRrcVHX1EKCySBjOYa4Fpc76OsMJUQaSYiUlqRAdtKYOMz1zVRoa9Ao/p9EIAZOcN9dxJB1qKqRf31bxDvbSzrdbrfbsVqtPqttIwJPuaWZ174+zLtefgcC3VcHqzlikS6cYKHTapg/KYlwZXCMixPdJ4mSCEptZQE6G8+h1+uZO3cuJpOJXbt2DXRoYoBEGHRUN9gpq7NR0xhcH14/lln4+lA1b28rpkFWqg8abreburo6Kioq+CF/P7WN/p91KfqfdL2JoNRWldnbel85OTlkZGQQERExkGGJAWQyhjA6PpxDVY3sKrVw2nHBUXeo0mrjs7xKAKaPiiXCILfhYJCXl8fq1avZuHEjTS1utu0vY2xGBvf8aTHjx4/3d3iiH0mLkgg6VlsLdU0tKAoMj+58XIfJZCI5OblPtW1E4Js0IgqAPWX1OF2BX+PG1uLi/R/KcbpVRseHc+KYWH+HJLohLy+PVatWsWPHDsLCwhgWF0t4ZAx7du3kvkceIy8vz98hin4kiZIIOm3LliSYjLJw5RA3Ki4ck1FHfWMTX+0pDOgxaepPRQstzS1EhoZwzsQkz5piInC53W7WrFmD2WwmKysLg8GAXqswLnUYiWljOXC4jDffege3O/ATddE7fUqUxo4dy1/+8hcqKyt9FY8Qx1Ra29btJtW4hzqNRiEzIYwPP/qIGx/4m0/WW+uJo8esFBYWdvlh+U1BDYeqGtFpFM6fnCyLOAeJoqIi8vPzSU1NbZfYJpoMRIWFEBmfxMZvczl8+LAfoxT9qU+J0qFDh7jzzjtJTU3l0ksvZf369b6KSwiv2gZyS6IkAMYnR2KKiCBhVBajMycN2HXz8vL461//ysaNG9m8eTMPPPAAjz76aKfdMKqqepZbmZ2VIKUAgojVasVmsxEeHt7ueUVRSB8WgTEsnPqGJr7bX+anCEV/61OidNdddzF8+HBaWlpYvXo1c+fOZezYsaxYsUJamUS/OXN8IielxzE8WhIlAWF6DVFNJcTX78fVZBmQLpBfjlmJi4sjLi6O3NxcVq1a1SFZUhSF8yYns2DaCCYMj+r3+ITvmEwmjEYjjY2NHbYZQ7QMM6ro9AbyzA4sTS1+iFD0tz4lSg8++CCHDx/mvffeY/78+Wg0Gg4dOsQdd9xBamoqv/71r6WVSfhcSkwYJ46Jk64L4WnV2f35B2zfsp4HH3zQa6uOr3Q2ZkWj0RAZGcn48eMxm82sXbsWt9uN0+VGVVuLEyqKwuj48GOcXQSatLQ0MjMzKS4u9vxftlFVFVtdJZmZmZx1/EQijDKDcTDq82BujUbDeeedx7p16zh8+DD3338/aWlptLS08Pbbb0srkxCiX/S0VcdXvI1ZgdZkKCUlhby8PA4fPswHu8r5+McjtATBjDzROY1Gw8KFC4mPjycvLw+bzYbb7cZisbBnzx6GDRvGnX+8kpPHDkMrqwQMSj6d9TZ8+HDuueceDh06xEcffcTChQvR6XQ+bWV6++23Of3004mJiSE8PJwpU6awYsUKWlp61uQ5atQoFEU55s8DDzzQ7rhNmzYd85hnn322169PdO3rQ9Xsr7DKB88Q11mrjqposTh1hMSltmvV8TVvY1bahIeHY7PZ2JJXzKGqRvZXNFDbFFwFMUV7WVlZLF26lKlTp9Lc3Ex1dTU1NTVkZ2ezdOnSdnWUXG4VW4vLj9EKX+uXdkJFUZg7dy5z586lrKyMyy+/nM2bN3vGMq1evZqxY8dyyy23cM0116DVdq8L5cYbb2TlypXodDpmz55NREQEn332Gbfffjvr1q3jk08+ITS0e+NWfvWrX2E2mzvdVlNTw7p16wA444wzOt0nMTGRc845p9NtGRkZ3YpB9Eyj3clXB6tRFPjDaelIz9vQ5a1Vp8xiQ1E0jEpIJi8vj6KiIkaNGuXTax89ZqWzZKmxsRG7qmVftZPYpJ8Gb5tk8Hawy8rK4rbbbqO6uhq73c5NN93E2LFj0Wh+bm+otNr4+McKYsJCOG/ycD9GK3yp3zpUi4qKeOGFF3jxxRcpK2udDaAoClOnTmXXrl3s37+fP/7xjzz//PN8+OGHDBvWdVXdtWvXsnLlSiIiIvj888/Jzs4GwGw2M3v2bLZs2cI999zDY4891q34utpvxYoVrFu3juOOO45TTjml030yMzN5+eWXu3Ut4RtlP1XjjoswyPikIa6zVh2tohIbrqe2yUm1HTQ2W7/UVWobs5Kbm9vhS5Gqquw7VIgmfgzRCcOZMSpWBm8PIhqNhujoaKC1V+LoJKlNTYMDs9XOgcoGxibIygCDgU+73lwuF2vXrmXevHmkp6fz0EMPUVpaSmxsLLfccgv79u3j+++/p7i4mHvvvZfw8HC2b9/OHXfcccxzP/zwwwAsW7bMkyQBxMfH8/TTTwPw5JNPYrFY+vw6XnzxRQB+97vf9flcwnfaCk2myGy3Ic/bTKTUmFA0ClTV1ONUdP1Smb2rMSs7du7GqoYx4eQzyUiOZObYOJ9fXwTuotcJJiPTR8UAsDG/UrrgBgmfJEqFhYXcddddpKWlcfHFF/Of//wHl8vFSSedxGuvvUZpaSl//etfSU9PB1q7rZYvX85nn32Gqqp89NFHXZ6/tLSU7777DoDLL7+8w/ZZs2aRmpqK3W7nww8/7NNr2bp1K3v37kWn03HVVVf16VzCt0rqpNCkaOVtJpIxREtCpIE68xF0sSmkpKT2y/U7G7NSXV2NPimdky/8LVMnTWDuBKm83R8CfdHrE0bHEhMWQoPdydYDnQ/vEMGlT11vq1ev5rnnnvMkPKqqYjKZ+O1vf8sf//hHJk6c2OXx06dPJykpiSNHjnS5X25uLgCxsbGMHj3a67mKi4vJzc1l0aJFvXtB/NyaNG/ePJKSkrzuV1FRwQMPPEBpaSlGo5HMzEzmz59PWlpar68tvGt2uDD/VLBvhCRKQ15bq05xcbGnVUev12OxWLCWlxIVHceY6Wew54iVySnR/RJDZ2NWDLHD2XqwhgumDidEKytE9ZdAXvRap9UwJyuR1d+XsLPEQlZypNR8C3J9SpQuvfRSz7+nTp3KH//4Ry6//HKvs0E6o9frj7lPQUEBQJdJSGpqart9e6OxsZG33noLgGuuuabLffPz87nvvvvaPafT6ViyZAkrVqxAp5N6Gr5U6hmfpCdML++t+LlVZ/Xq1Rw4cACr1UpNTQ0zpk/nuOPPoMgdy7cFNUwcHoWmn6Zt/3LMil6vJy0uXFqS+pnJZAroBa9TY8OYMDySH8vq2ZBXweUnjJTSAUGsT584RqORX//61/zxj3/k+OOP79U5CgsLj7lPWz90VwlY2zeL+vr6XsUB8NZbb9HQ0EBSUhLz5s3rdJ+oqChuvPFGFi5cyHHHHUdkZCQHDx7kpZde4sknn+Tvf/87DQ0NPPfcc11ey263Y7fbPY/7EvdQUGm1ATBCvpmJo3ibiaSisPWAmalp0f2WJAE4XW4KXTEkaH4eKyNJkgA4ZdwwDpkbMYRosbW4CDfIF7xg1af/ubKyMs+3qcHghRdeAODKK6/02iI0bdo0pk2b1u65SZMm8fjjjzNr1iwuvvhinn/+ea677jqmTp3q9VqPPPII999/v89iH+xOTo9vnT2kHntfMbR4m4l06nFdz6TtK6fLzYe7K6h1h9Ko6nG55ZdT/CxUr+WyGalEhYZI8hzk+tSJPlBJUlsTa2dr7bRpaGgAIDIyslfX2LdvH1u3bgV6P9vtoosu8iRHbXWYvLnjjjuwWCyen+Li4l5dcyiJCg0hKizE32GIAHSsWVD5R+p9WqTU5Vb5YFc5hdVNaBSVkdpa6VoRHUSH6SVJGgSCYrRhW8G4rpKJtm29LS7XNoh71qxZfSoYmZWVBUBJSUmX+xkMBiIjI9v9CCF67lizoL4+VM1Hu47w4a5yn7T6NNqdvLO9hENVjYRoFcZoqolQpPK28M7udLFpbyV7jwRWOQPRPUHRadrW1VVdXU1BQUGnM9+2bdsG0K7GUne5XC5effVV4NiDuI+luroaIKAHGgabbw5Vc6TextTUaEbGyaKioqOuZkGlxISi0ygcqmrk0z0VzJ2Q2Otv+WV1zXyws5wGuxO9TsPcrAT+uVWSJNG13aUWcovq2Ku3MjIuTArmBpmgaFFKSUlhxowZALzxxhsdtm/ZsoXi4mIMBoPXQdhd+fDDDykvL8dkMnHJJZf0Os7S0lK++OILgF4PbhcdHTI3cqiqkUa7FG8TnTOZTCQnJ3f6BSUlJox5k5PRKAp55fVs3m/usAp8d/1QXEeD3UlchJ7Lj08jLTasr6GLIWBqagyx4XqaHC6+PCi1lYJNUCRKAHfeeScAjz76KNu3b/c8X11dzXXXXQfA9ddfT1TUz8sFrFmzhszMTObMmdPludu63S677LJjljZYuXJlp2vE7dy5k/PPP5/m5mbS09O58MILu/fCRJfsTheV9a2zA1NiZcab6J30YRGcOT4BgO2Ha/nyYHW3xyw5j9pvdlYCOSNjuGxGGjHhxy5tIgSAVqMwO7P1929niYUjFpufIxI9ERRdbwALFixg6dKlrFq1ihNPPJE5c+YQHh7Ohg0bqKurY+bMmTz44IPtjrFYLOzduxebzfsvZWVlJR988AHQvW63++67j1tuuYWpU6cyevRoNBoNBw8eJDc3F7fbTVpaGuvWrcNgMPTtBQsAjlhsuFWVyNAQIo0ykFv03oThUdhaXGzeZ+bbghpKapv49QzvtdnKLc18W1BDk8PFZTNaF9816LT9PptODE6psWFkJZvIK7fyWX4ll81I7dfSFcJ3giZRgtbWnJkzZ/LUU0/x5Zdf0tLSQnp6OsuWLeOmm27qVvHKX3rttddoaWlhwoQJnHDCCcfc/6677mLr1q38+OOPfPrppzQ2NhIZGcnJJ5/MhRdeyLXXXivjk3yobX03qZ8kfCFnZCwGnZZvCmrISv55AkVxTRMHqhpwON04nG4a7E7Pt35FgaoGOwkmo7/CFoNEW22linobO0stTE2N9ndIohsUtbed9cKn6uvriYqKwmKxyAy4o7z1XTGldc2cNT6RiSNkFXbhG22z39qm9OcW1bJpb1W7fTSKQlayiemjYon10s3mcDi47777cDgc3HnnncTFySK4g5nD4fAs0H7nnXf26sv5D8V1fJZficmo4+qZo6WshJ/05DM3qFqUxNDS4nJzpL71W70shCt86ZcfTklRRo4fHYtep0Gv1aDXaRgRE3rM7t620gSbNm1i165dnH766f0YtRgMJo2IotHuZEpqtCRJQUISJRGwmhwukqOMNNidRIXK+CTRf5KjQkmO6l0yHsgLtIrAo9EonDw23t9hiB6QREkErKjQEC6ZnorbrUp1WxGwAn2BVhHYimuaGB4dKq1LASxoygOIoUtmhgghBqMNeRWs/r6E7UW1/g5FdEESJRGQnC43zQ4pMCmEGLySolpnUn5zqJp6W4ufoxHeSKIkAlJRTRP/2HyQ934o83coQgjRL8YnRzIiOpQWl8rmfVXHPkD4hSRKIiAV1zajqhAqayIJIQYpRVE4IzMBjaKwv6KBQnOjv0MSnZBESQSkktomAFJl2RIhxCA2zGRgSmprjbiNeyvbLZkjAoMkSiLgNDtcVFlb13dLjZFFR4UQg9tJ6XFEGHTUNbXw/WEZ2B1oJFESAae0rglVhbgIPeEGqWAhhBjcDDotpxwXT2RoCPEmWSc00MinkAg4xTWt67tJa5IQYqjISDSRPiyCEK20XwQa+R8RAadYxicJIYYYRVHaJUmyDGvgkERJBBRVVclOiyEr2cSIaGlREkIMLW63ys6SOv79XbEM7A4QkiiJgKIoChNHRHHOxGRC9VIaQAgxtLS43XxzqIZyi43tRXXttlmtVsrLy7Farf4JboiSREkIIYQIEG0DuwG+LajG0txasdvhcLBkyRIWLVrE119/7c8QhxxJlERA2V1qocpql/55IURAstvtWK3Wfm3VyUg0kRLTsWJ3cnIyOTk5ZGdn99u1RUeSKImAYWlu4dM9FbzxTREO6ZsXQgQYvV7P3LlzMZlM7Nq1q9+uc3TF7gOVP1fsNhgMmEwmTCZTv11bdCTlAUTAaKvGnRhpwKCT8UlCiMCTk5NDRkYGERER/Xqd+AgDU9Oi2X64lk17K7k0O7lfrye8k0RJBIyi6rayADLbTQgRmAayReeE0bHsPVJPXXMLpXW2Abmm6EgSJREQ3G6VwzWtidLIOEmUhBDCGKLl7PFJhOq1xBhlpIy/SKIkAkKl1U6zw4VepyE5SgpNCiEEwKj4cKB11pvwD0mUREA4XN06WDEtNgytRvFzNEIIEXhsqg4Xcn8caJIoiX5jtVppaGggIiLimH36xbWt67tJt5sQQnRUYG5kr2sYely43FI+ZSBJp6foFz0tjnbh1OFclD2C9GH9O5NECCGC0fAoI1pUbKqOHcUWf4czpEiiJPpNT4qjhWg1jIwLJ9wgjZxCCPFLhhAtyZrWBOnbwlqsthY/RzR0SKIk+o0URxNCCN+JVZqJUBy0uNxs3mf2dzhDhiRKwq9UVeXd7SVs3ldFs8Pl73CEECJgKQqM0NahKLCvwuqpPSf6lyRKwq9qGh0crm7ih+I6dFqZzSGEEF0JU5xMHhEFwMa9lTKwewBIoiT8qvCnb0QjYkIJ0cqvoxBCHMuJo2OICg1h/PBIf4cyJMjIWeFXbfWTRsaF+zkSIYQIDoYQLVedPEpqzg0Q+Qov/KbF5ab0p/pJo6R+khBCdNvRSZJbut/6lSRKwm9Ka5txulVMRh2x4Xp/hyOEEAHL7XZTV1dHRUUFhYWFuN1uAIprmvjnN4c9rfPC96TrTfhN4VHdbooiTchCCNGZvLw8Vq9ezcaNG3E6nVitVsaPH8/ChQs5oomjusHBxvxKfnviSHQy1tPn5B0VfmMM0RKm10q3mxBCeJGXl8eqVavYsWMHYWFhxMXFERcXR25uLqtWrSLGXkm4QUttUwvfH671d7iDkiRKwm9OHBPH708dwxhZtkQIITpwu92sWbMGs9lMVlYWBoMBjUZDZGQk48ePx2w289EH65iZHgfAtwU1WJqkYrevSaIk/EpRFJm5IYQQnSgqKiI/P5/U1NQOwxMURSElJYW8vDzCHLWkxITidKts2lfpp2gHL0mUhF9UWe2oqszUEEIIb6xWKzabjfDwzsunhIeHY7PZaGhoYHZmAhpF4VBVIwerGgY40sFNEiUx4BrsTl7/5jD/+0UBDqfb3+EIIURAMplMGI1GGhs7n9HW2NiI0WjEZDIRF2Ege2Q0AAcqJVHypaBLlN5++21OP/10YmJiCA8PZ8qUKaxYsYKWlp71y7788ssoitLlz3/+8x+vx1dUVHD99dczevRoDAYDiYmJXHLJJWzfvr2vL3HQK6hqRFUhwqhDrwu6X0EhhBgQaWlpZGZmUlxc3KEFXlVVSkpKyMrKIi0tDYATRscxf3IyZ49P9Ee4g1ZQlQe48cYbWblyJTqdjtmzZxMREcFnn33G7bffzrp16/jkk08IDQ3t0TnT09OZNWtWp9tGjBjR6fP79u3jlFNOobKykjFjxrBgwQIKCgpYvXo1a9eu5a233mLhwoU9fn1DRVuzcLoM4hZCCK80Gg0LFy6kuLiYvLw8bDYber0ei8VCeXk58fHxLFiwAI2m9QunXqfhuESTn6MefIImUVq7di0rV64kIiKCzz//nOzsbADMZjOzZ89my5Yt3HPPPTz22GM9Ou+sWbN4+eWXu72/qqpcdtllVFZWcsUVV/DSSy+h1WoBeO6557j22mu58sor2b9/P0lJST2KZSiwO10U1bSu75Y+TJYtEUKIrmRlZbF06VJWr17NgQMHsFqt1NTUkJ2dzYIFC8jKyur0OFuLiz3l9UxLjZY6dX0UNP0eDz/8MADLli3zJEkA8fHxPP300wA8+eSTWCyWfo3jo48+Ijc3l+joaJ5++mlPkgTw+9//njlz5tDQ0MDKlSv7NY5gdbi6CZdbJSYsRKpxCyFEN2RlZXHbbbdxxhlncOqpp3Lvvfdy++23e02SXG6VN74p4vO9VeQfsQ5wtINPUCRKpaWlfPfddwBcfvnlHbbPmjWL1NRU7HY7H374Yb/GsmbNGgAuuOACIiI6dh21xffuu+/2axzB6uBPgwzTEyLkW44QQnSTRqMhOjqaxMRERo0a5elu64xWozBxRBQAm/dVYWtx9fn6VquV8vJyrNahl3gFRaKUm5sLQGxsLKNHj+50n+nTp7fbt7sOHDjA3Xffze9//3tuvvlmXnzxRcxm8zFjabuetzj279/vdabCUOVyqxT8tGyJjE8SQoj+kzMyhrgIPU0OF1sPeP9M6w6Hw8GSJUtYtGgRX3/9tY8iDB5BMUapoKAAwDOyvzOpqant9u2urVu3snXr1nbPGY1Gli9fzu23397jWNriUFWVwsJCJkyY0KN4BjONAhdOHcFhcyNJkUZ/hyOEEIOWVqNwRkYCq78vYVephfHDI0mO6tlkp6MlJycTFxfXbujLUBEULUptTX3eim4Bnm6w+vr6bp0zKSmJu+66i2+++Yaqqirq6+v57rvvuPLKK7Hb7SxbtswzLqonsRzdHddVLHa7nfr6+nY/g52iKIyIDuXksfFopBq3EEL0q9TYMLKSI1FV2JBXidvd+yK/BoMBk8mEyTT0ZtUFRaLUH8455xweeughjj/+eOLj4zGZTEyfPp1XXnnFM3PugQceoKKiol+u/8gjjxAVFeX5aWuJEkIIIXzl1OPiMYZoqbLa2VFS5+9wglJQJEptGWxXY34aGloHCUdGRvb5ejfccAPx8fHY7XY++eSTHsXSFsexYrnjjjuwWCyen+Li4j7HHcgqrTY25FVQ/FNpACGEEP0vTK/jlHHxjB8eSVZS3z8fh6KgGKM0atQogC6TibZtbfv2hVarZdy4cZjNZkpKSjrEUlNTQ1FRUZdxKIrCyJEjvV7DYDBgMBj6HGuw2FNWz84SC3anm9TYMH+HI4QQQ8bEEVGeWXCi54KiRWnatGkAVFdXex2svW3bNgCfDTSrrq4G6NAf23b+tut5i2PcuHGdlg8YKtxuN3V1dVRUVFBQUEB+eWt9q4ykode/LYQQgUJVVRrsTn+HEVSCIlFKSUlhxowZALzxxhsdtm/ZsoXi4mIMBgPz5s3r8/W2b9/Ovn37ADj++OPbbWtbmuS9997rtPutLb6LLrqoz3EEq7y8PP7617+yceNGNm/ezLK772Pdq89QW1bAqDipxi2EEP7Q7HDx3g9lvPHNYZ/UVhoqgiJRArjzzjsBePTRR9stPFtdXc11110HwPXXX09U1M/Ni2vWrCEzM5M5c+a0O1dTUxNPPfVUp4WzNm/ezMUXXwy0FrL8ZaJ07rnnMm3aNOrq6rjuuutwuX7+ZXvuuefYsGEDERER3HDDDX18xcEpLy+PVatWsWPHDsLCwoiLiwODiZIDe8j96A327c33d4hCCDEkhWgVahsdNNr7XltpKAmKMUoACxYsYOnSpaxatYoTTzyROXPmEB4ezoYNG6irq2PmzJk8+OCD7Y6xWCzs3bsXm83W7nmHw8H111/PLbfcwrRp00hLS8PpdLJv3z52794NwKRJk3jrrbc6xKEoCv/617845ZRTePXVV9myZQszZsygoKCAb7/9Fp1Ox6uvvjok13lzu92sWbMGs9lMVlYW1dXVqCrYNAYS08bitpaydu1aMjIyuqwqK4QQwvd0Wg1zshJZ/X0JO0ssZCZHMiK697WVhoqg+rRauXIl//73vznppJP48ssv+fDDD0lJSeHRRx/ls88+IzS0e//hYWFh3HPPPcyePZuKigo++ugj3nvvPSoqKjjzzDP5xz/+wbZt20hOTu70+IyMDHbu3Mmf/vQnXC4Xa9asoaCggIsuuohvvvnG0z031BQVFZGfn09qaqpneRI7OlxuFWOIlnGjR5KXl+d1ILwQQoj+lRob5hnYvX5PBU6X288RBb6gaVFqc+mll3LppZd2a9/FixezePHiDs/r9XoeeOCBPsWRlJTEk08+yZNPPtmn8wwmVqsVm83WrhinChh0WuIjDERE6CkrKxuSawUJIUSgOGVcPAXmBmoaHXxbWMPJ6fH+DimgBVWLkghsJpMJo9HYbpB7qOJkSkokI2JCaWxsxGg0DsnKrkIIESiMIVrOyEgA4LuCWqqsdj9HFNgkURI+k5aWRmZmJsXFxajqz6XyFUVBo0BJSQlZWVldrtknhBCi/41NiCA9IQKTUYfTLd1vXQm6rjcRuDQaDQsXLqS4uJi8vDyszS2E6bVYLBbKy8uJj49nwYIFMpBbCCF6wW6343A4sFqtrTOK+0BRFM7KSkSnVQjRyj25K/LuCJ/Kyspi6dKljJ84mcomFwXVjVSZq8nOzmbp0qVkZWX5O0QhhAg6er2euXPnYjKZ2LVrl0/OGarXtkuSju4JED+TFiXhc1lZWcz7zbV8bQ4hpKWB+29dzNixY6UlSQgh+iAnJ4eMjAyfr/qgqio7iusoqmniginDPbOWRStJlITPqapKfmUjYdHxpGpDGDVqlCRJQgjRRyaTqV8mw9TbnGw9YKbFpbK7tJ5JKbIu3NHk00v4XJnFRnWDA42iEq00+TscIYQQXYgKDeGkn0oEbN5fRb2txc8RBRZJlITP7SqpAyBGaUanSJ+3EEIEummp0QyPNuJwulm/p0LGKx1FEiXhU80OF/srGgCI03RcNFgIIUTg0WgUzhqfhE6jcLi6iR/L6v0dUsCQREn41OGaRpxulQSTgTCk+VYIIYJFbLiek9Jbyw5s3l+FVbrgABnMLXwsMymSYREGGprtHPrU39EIIYToiey0GPZXNlBZb6eszkZGUoi/Q/I7SZSEz8VFGDDpZXqpEEIEG41GYe6EJFpcbhIjjf4OJyBIoiR8xuF0o9dJb64QQgSz2HC9v0MIKPKpJnyi0e7k+S8O8eGucpwuWTdICCEGgyqrnY17qxjKk+CkRUn4xJ7yehxON/XNLei0Ghwuf0ckhBCiLxxON6u/L6HR5qBKDSdBGZozmaVFSfSZ262yq8QCIBVdhRBikNDrNMwaG4/b7eZAjYPiI2YKCwtxu4dWr4G0KIk+O1jVgKW5BWOIluMSfV9eXwghhH9o60vZ8f4r5H2xmXynA5u1lgkTJrBw4cIhs8i5tCgNclarlfLycqxWa7+cX1VVviusBWBKalS7laiFEEIEr7y8PJ544gkay/ZjCAvHGJuMQxdObm4uq1atIi8vz98hDgj5VBvEHA4HS5YsYdGiRXz99df9co3immYq6m2EaBWmpcb0yzWEEEIMLLfbzZo1azCbzUyaOIF4o4Ki0WBxhZAyehxms5m1a9cOiW44SZQGueTkZHJycsjOzu6X8+cWt7YmTRgeRahe2y/XEEIIMbCKiorIz88nNTUVRVEwKk7ClNZK3VUNdlJSUsjLy6OoqMjPkfY/GaM0yBkMBgwGAyZT/4wdmjshiZ0lFjKSZGySEEIMFlarFZvNRnh4uOc5k2JjZFwYw6PDcLlclJaW9tuwjkAiLUqiT4whWo4fHUtUqJS5F0KIwcJkMmE0Gmls/LkkgEaBpEgjiqLQ2NiI0Wjsty/hgUQSJdErTpcbdShXIBNCiEEsLS2NzMxMiouLO9zrVVWlqKgY47A0ouKT/RThwJFESfTK5/uqePO7Ykrrmv0dihBCCB/TaDQsXLiQ+Ph48vLysNlsuN1uLBYLe/bswa4LJ3b8TP6zpwKXe3B/aZZESfRYo93JnrJ6jlhsuAf5H4gQQgxVWVlZLF26lKlTp9Lc3Ex1dTU1NTVkZ2dz77JbSEsfxxGLja0HzP4OtV/JYG7RY9sO1+J0qyRHGUmJCfV3OEIIIfpJVlYWt912G9XV1djtdm666SbGjh2LRqMhtrKBdT+U8f3hWlJiQhkzLMLf4fYLaVESPWJpbuGH4joAThwTh6Io/g1ICCFEv9JoNERHR5OYmMioUaPQaFpTh7EJEUxNiwbgkz0VWG0tfoyy/0iiJHrkq4PVuNwqqbFhjIwL83c4Qggh/OiUsfEkRBpodrj4aPeRQTkcQxIl0W1VVjv5R+oBmDU2/pitSXa7HavVOiTqbAghxFCk02qYPykZvU5DdYOD2iaHv0PyOUmUBjG3201dXR0VFRU+WfH5xzILqgrHJZpIijJ2ua9er2fu3LmYTCZ27drVp+sKIYQIXNFhes6bnMxvTkwjLsLg73B8TgZzD1J5eXmsXr2ajRs34nQ6sVqtjB8/vk8rPp923DCSo0JJMHXvDyEnJ4eMjAwiIgbnAD8hhBCtRsaFt3usquqgGcMqidIglJeXx6pVq6isrCQsLIyQkBDi4uLIzc2luLiYpUuX9ipZUhSlR0uVmEymIVG1VQghxM8OVjWwo6iOC6cOR6fte8eV1WqloaGBiIgIv3ymSNfbIHP0is9ZWVkYDAY0Gg2RkZGMHz++Vys+VzfYcTgH/wrRQggh+sbudLF+TwVFNU1s3FvV5/M5HA6WLFnCokWL+Prrr30QYc9JojTI/HLF56MpitLjFZ+dLjfv/VDGy18WUG6RKtxCCCG8M+i0nDMxCUWB3aUWdpda+nzO5ORkcnJyyM7O9kGEPSeJ0iDT2YrPRwsPD8dms3V7Jtq3hTXUNbWgoBAbrvdlqEIIIQahkXHhnDQmDoCN+ZUcsdj6dD6DweDXoRySKA0yna34fLSerPhc0+hgW2EtAKdlDMOg0/o0ViGEEIPT8aNjGTMsHKdbZd0PZTTYnf4OqdckURpkjrXic0lJCVlZWaSlpXV5HlVV+Sy/EpdbZXR8OOMSZOaaEEKI7lEUhXMmJhEXoafB7uSDnWU4XcE51lUSpUHmWCs+x8fHs2DBAk8Jem/yj1gprmlCp1E4IyNh0EzzFEIIMTAMOi3nTx6OIUTDMJMhaD9Hgi5Revvttzn99NOJiYkhPDycKVOmsGLFClpaerbGTG5uLo888ghz5swhMTGRkJAQYmJiOOWUU3jqqae8nm/Tpk0oitLlz7PPPuuLl9prXa343J3SALYWF5v3tc5WOGFMHFFhIQMRthBCiEEmJlzPFSeOZHZmIlpNcCZKQVVH6cYbb2TlypXodDpmz55NREQEn332Gbfffjvr1q3jk08+ITT02KvZO51Oz+j5iIgIZsyYQWJiIiUlJXz11Vds2bKFV199lY8//pjo6OhOz5GYmMg555zT6baMjIxev0Zf6WrF52NR1dbBeFVWGzkjYwYgWiGEEIOVyfjzl22XW6Wm0cGwbhYuDgRBkyitXbuWlStXEhERweeff+5JdMxmM7Nnz2bLli3cc889PPbYY906X05ODrfffjsXXHABBsPP/2G7du1i7ty5fPvtt9x88828+OKLnR6fmZnJyy+/3OfX1Z/aVnwG2q34fCyh+tbpnQ6nO2i/AQghhAgsdqeL938o50i9jV/PSCU+SJY7CZqut4cffhiAZcuWtaulEB8fz9NPPw3Ak08+icVy7JoNOp2Obdu2cckll7RLkgAmTZrEihUrAHjzzTd73KUXzFpc7nYDwPW6oPn1EEIIEeC0ioJLVXE43azNLaUxSGbCBcUnYWlpKd999x0Al19+eYfts2bNIjU1Fbvdzocfftjn602bNg2A5uZmzGZzn88XDFRV5aPdR/i/HWU0OYLjl1cIIUTw0Gk1nD95ODFhIVhtTv5vRxktQTATLigSpdzcXABiY2MZPXp0p/tMnz693b59sX//fgD0ej2xsbGd7lNRUcEDDzzAtddeyw033MAzzzzT7WrXgWh3aT0HKxsoqmmiwSaJkhBCCN8L1WtZMG0EoXotFfU2Ptp9BLdbPfaBfhQUY5QKCgoAuqz9k5qa2m7f3lJV1dP1dt5553XommuTn5/Pfffd1+45nU7HkiVLWLFiBTpdULy1QGthyc/3VQIwc2wcCZFGP0ckhBBisIoO03P+lOG8+30JBysb2Ly/itOOGxaw5QOCokWpbbkNb8tyQOvsNYD6+vo+Xev+++/nq6++IiIigkcffbTD9qioKG688UY+//xzysvLaWxsZOfOndx0000oisLf//53rrvuumNex263U19f3+7HH5wuNx/tLqfFpTIyLozsNJnlJoQQon+NiA7l7AlJAPxYVo81gMcrBU+zxwB49dVXeeCBB9BoNLz44ouMGzeuwz7Tpk3zjGFqM2nSJB5//HFmzZrFxRdfzPPPP891113H1KlTvV7rkUce4f777/f1S+ixL/abqay3E6rXcvaEpIDN6IUQQgwuGUkm7E4XSVFGIo2d1+tzu93U1dVht9spLCzsdpkbXwqKRKltXTJv65cBNDQ0ABAZGdmra7z99tv87ne/A+D555/nkksu6fE5LrroIqZOncqOHTtYt25dl4nSHXfcwc033+x5XF9f7+k+HCi7Sy3sKK4D4KzxiUQYguLXQQghxCAxOSW63eMWl5sQbWsilJeXx+rVq9m4cSNOpxOr1cr48eNZuHDhMQsn+1JQfDKOGjUKgOLiYq/7tG1r27cn3n33XS6//HLcbjf/+Mc/PAlTb2RlZbFjxw5KSkq63M9gMHgd/zRQEiINRIaGMHF4JOnDZC03IYQQ/lNa18yHO8s5d1IS1iOHWbVqFZWVlYSFhRESEkJcXBy5ubkUFxd3a5UJXwmKMUptXV3V1dVeB2tv27YNoF2Npe5Yu3Ytl112GS6Xi2eeeYb/9//+X59ira6uBn5uBfM3u92O1Wr1jPM6WoLJyG9OSOP40Z3P7BNCCCEGyo6iOhrsTtbmlvDav97GbDaTlZWFwWBAo9EQGRnJ+PHjMZvNrF27Frd7YEoLBEWilJKSwowZMwB44403OmzfsmULxcXFGAwG5s2b1+3zrlu3jksvvRSn08kzzzzDtdde26c4S0tL+eKLLwA4/vjj+3QuX9Dr9cydOxeTycSuXbsAcDjdHLHYPPsYQ7QyLkkIIYTfnT0hkRHRoVSUlvDZN7nEJSR3+HxSFIWUlBTy8vIGrCRPUCRKAHfeeScAjz76KNu3b/c8X11d7Zlldv311xMVFeXZtmbNGjIzM5kzZ06H83344Yf86le/wul08uyzz3Y7SVq5cmWnRSh37tzJ+eefT3NzM+np6Vx44YU9en39JScnh2uvvZacnBzcbpX//HiEt7YVk3/EP7PshBBCBJ+ueid8JUSr4YKpw4nQOrE12yiqd9HscHXYLzw8HJvN1q+xHC0oxigBLFiwgKVLl7Jq1SpOPPFE5syZQ3h4OBs2bKCuro6ZM2fy4IMPtjvGYrGwd+9ebDZbu+crKyu56KKLcDgcpKSk8OWXX/Lll192et3HHnuM+Ph4z+P77ruPW265halTpzJ69Gg0Gg0HDx4kNzcXt9tNWloa69at8/v4ozYmkwmTyYSqqnyyp4KDlQ1oNYrXGQZCCCHE0dp6JzZt2sSuXbs4/fTT++1axhAtC2aks+6fYTQ2NpF3RINTVdApPxelbGxsxGg0DtgQl6BJlKC1NWfmzJk89dRTfPnll7S0tJCens6yZcu46aab0Ov13TpPU1MTdrsdgJKSEl555RWv+y5fvrxdonTXXXexdetWfvzxRz799FMaGxuJjIzk5JNP5sILL+Taa68NmPFJbVRVZdO+KvaU1aNRFOZNSmZ4dKi/wxJCCBEkcnJyyMjI8NQs7E/j0kdz1knT+GDjl+hDQ2lUDUQprQ0eqqpSUlJCdnZ2l0WofUlRj14FVfhNfX09UVFRWCyWXpc48ObLA2a+KagB4JyJSWQl+/b8QgghhC/l5eXxt7//D/sKS7FbzBgMeqZOnUp5eTnx8fF9nvXWk8/coBmjJHpnW2GNJ0manZkgSZIQQoiAl5WVxS033chZp5yAzdZMdXU1NTU1ZGdnD2hpAAiyrjfRM6qqUtPoAGDWuHimpEb7NyAhhBCim7Kysrjtttuorq7Gbrdz0003SWVu4VuKonDW+ETGJkQwRgpKCiGECDIajYbo6GigtaD0QCdJIF1vg56iKJIkCSGEEL0kiZIQQgghhBeSKAkhhBBCeCGJkhBCCCGEF5IoCSGEEEJ4IYmSEEIIIYQXkigJIYQQQnghiZIQQgghhBeSKAkhhBBCeCGJkhBCCCGEF5IoCSGEEEJ4IYmSEEIIIQKW3W7HarVitVr9cn1JlIQQQggRkPR6PXPnzsVkMrFr1y6/xKDzy1WFEEIIIbohJyeHjIwMIiL8s8C7JEpCCCGECFgmkwmTyeS360vXmxBCCCGEF5IoCSGEEEJ4IYmSEEIIIYQXkigJIYQQQnghiZIQQgghhBeSKAkhhBBCeCGJkhBCCCGEF5IoCSGEEEJ4IYmSEEIIIYQXkigJIYQQQnghS5gECFVVAaivr/dzJEIIIcTg1vZZ2/bZ2xVJlAKE1WoFIDU11c+RCCGEEEOD1WolKiqqy30UtTvplOh3brebsrIyTCYTiqL47Lz19fWkpqZSXFxMZGSkz847WMj70zV5f45N3qOuyfvTNXl/jq0/3iNVVbFarQwfPhyNputRSNKiFCA0Gg0pKSn9dv7IyEj5I+yCvD9dk/fn2OQ96pq8P12T9+fYfP0eHaslqY0M5hZCCCGE8EISJSGEEEIILyRRGuQMBgP33XcfBoPB36EEJHl/uibvz7HJe9Q1eX+6Ju/Psfn7PZLB3EIIIYQQXkiLkhBCCCGEF5IoCSGEEEJ4IYnSEFNZWcmrr77K5Zdfzrhx4zAajYSFhZGZmcnSpUspLCz0d4gD4u233+b0008nJiaG8PBwpkyZwooVK2hpafF3aH7V0tLChg0buO2225gxYwbR0dGEhISQlJTEBRdcwAcffODvEAPSn//8ZxRFQVEUHnroIX+HExAcDgerVq1i1qxZxMbGYjQaSUlJ4dxzz+Xf//63v8Pzu6KiIq6//noyMjIIDQ3FaDQyevRorrrqKn744Qd/h9fv9u7dyxNPPMHixYuZNGkSOp2u238/69evZ968ecTHxxMaGkpmZiZ33XUXDQ0N/ROsKoaU3/zmNyqgajQadfLkyeoll1yizps3Tx02bJgKqOHh4eonn3zi7zD71Q033KACqk6nU88++2z1oosuUqOjo1VAnTVrltrU1OTvEP3m008/VQEVUJOSktT58+erl156qTpx4kTP87///e9Vt9vt71ADxtatW1WNRqMqiqIC6oMPPujvkPyuuLhYHT9+vAqo8fHx6nnnnaf++te/Vk8++WQ1LCxMvfjii/0dol99/fXXqslkUgF1xIgR6gUXXKAuXLhQHT16tOfe9NZbb/k7zH7Vdh/+5c+x/n4ef/xxFVAVRVFPPfVU9ZJLLlGTkpJUQM3IyFCrqqp8HqskSkPMkiVL1Pvvv18tKSlp97zValUvu+wyFVBjY2PVmpoaP0XYv9asWaMCakREhPr99997nq+qqlInTZqkAuott9zixwj9a8OGDerFF1+sbt68ucO2N998U9VqtSqgvvLKK36ILvA0Njaq48aNU0eMGKEuWLBAEiVVVZuamtTMzEwVUJcvX646HI522xsbG9Xc3Fz/BBcgJk+e7PnScfT743K51LvvvlsF1OjoaLW5udmPUfav559/Xr311lvV119/Xc3Ly1OvuOKKY/79bN++XVUURdVqteqHH37oeb6xsVGdM2eOCvRLEi6JkvBobGz0fMt57bXX/B1Ov5gxY4YKqA899FCHbV988YUKqAaDQa2rq/NDdIHvmmuuUQF1zpw5/g4lICxdulQF1A8++EC96qqrJFFSVfWee+7xJAGiI7PZ7Gk9qays7LDd6XSqoaGhKqBu377dDxH6R3f+fi655BIVUP/rv/6rw7bCwkJVo9GogJqXl+fT2GSMkvAICwsjIyMDgOLiYj9H43ulpaV89913AFx++eUdts+aNYvU1FTsdjsffvjhQIcXFKZNmwYMzt+Pntq0aRNPPPEEV155JfPmzfN3OAGhpaWFZ555BoDbbrvNz9EEpp7UAoqPj+/HSIKLw+HwjJHs7P49cuRIZs6cCcCaNWt8em1JlIRHS0uLZzB3cnKyf4PpB7m5uQDExsYyevToTveZPn16u31Fe/v37wcG5+9HTzQ0NPC73/2OxMRE/ud//sff4QSM7du3YzabGT58OGPHjmXXrl3cf//9XHvttSxbtowPPvgAt9vt7zD9KiIiglNOOQWAu+++u90EErfbzfLly2lububcc88lNTXVX2EGnH379tHU1AT8fJ/+pf66f8uiuMLjhRdewGw2ExoayrnnnuvvcHyuoKAAgLS0NK/7tN2Y2vYVPzty5Agvv/wyABdffLF/g/GzW2+9lYKCAtasWUNMTIy/wwkYO3fuBCAlJYVly5axYsUK1KNqGv/lL39h2rRprF27tsu/w8Hu+eefZ968eTz33HN88MEHTJ8+Ha1WS25uLqWlpVxxxRU8+eST/g4zoLTdk6OjozGZTJ3u01/3b2lREgDs2rXL01R+zz33kJiY6OeIfM9qtQIQHh7udZ+IiAgA6uvrBySmYOF0Ovntb3+LxWJh0qRJXHvttf4OyW8++eQT/vGPf3DZZZexYMECf4cTUKqrq4HWb/R/+ctfuO6669i7dy8Wi4VPP/2U4447jtzcXObPnz+kS3FkZGTw1VdfcfbZZ1NaWsr//d//8e6771JQUMDYsWM5/fTTiYyM9HeYAcWf929pUQoif/7zn3nvvfd6fNz//u//MmvWLK/bS0pKOP/882loaOCCCy5g2bJlfQlTDEJ/+MMf2LBhA3FxcaxevRq9Xu/vkPzCYrFwzTXXMGzYMJ544gl/hxNw2lqPWlpaWLRoUbtWkTPPPJNPP/2UjIwMdu/ezZtvvskVV1zhr1D9auvWrVx00UXodDreeOMNZs+ejV6vZ+vWrdx8881cc801bN26lRdeeMHfoQokUQoqZWVl7N27t8fHdVWE68iRI8yZM4fDhw8zd+5c3nrrLRRF6UuYAautubaxsdHrPm3vlXyb+9kNN9zACy+8QExMjKdVYKi68cYbKSkp4d///rcMtO3E0V0inbU6pqWlMX/+fN555x3Wr18/JBOluro6Fi5ciNls5quvvuKEE07wbDvvvPMYP348kyZN4sUXX+S3v/0tZ5xxhh+jDRz+vH9LohRE/vnPf/LPf/7TZ+errKxk9uzZ7Nu3jzPPPJO1a9cO6hWsR40aBXQ9Y6ttW9u+Q90tt9zCqlWriI6O5pNPPvHMehuq1qxZg06n4+mnn+bpp59uty0/Px9oHeu3fv16kpKSePPNN/0Rpt+MGTOm0393tk95efmAxBRoPvjgA6qqqkhPT2+XJLUZM2YMJ5xwAhs3bmT9+vWSKP2k7Z5cV1eH1WrtdJxSf92/JVEaoqqqqpg9ezZ5eXnMmTOH9957D6PR6O+w+lXbh3x1dTUFBQWdznzbtm0bANnZ2QMaWyD685//zOOPP05UVBSffPKJ15kmQ43T6eTzzz/3ur2wsJDCwkJGjhw5gFEFhuzsbBRFQVVVzGZzp7O2zGYz8PN4kqGmqKgI6LrVIyoqCoCampoBiSkYZGRkEBYWRlNTE9u2bes0geyv+7cM5h6CzGYzs2fP5scff2TOnDmsW7eO0NBQf4fV71JSUpgxYwYAb7zxRoftW7Zsobi4GIPBMOTr4ixbtoy//vWvREVF8emnn3ret6Gurq4OtbVQb4efq666CoAHH3wQVVWHzLqJR0tKSvKMh1y/fn2H7S0tLZ4k8/jjjx/Q2ALFiBEjgNYWSIvF0mF7S0sL27dvB/BaxmQo0uv1zJ8/H+j8/n348GG+/PJLABYuXOjbi/u0fKUIeNXV1Z7y+WeeeeaQW9fM2xImZrNZljD5yV133eVZQuHbb7/1dzhBQypzt1q/fr0KqDExMepXX33leb6lpUVdsmSJCqgmk0k9cuSIH6P0n8rKSjU8PFwF1EsuuUS1Wq2ebXa7Xf3Tn/6kAmpISIh68OBBP0Y6sLrz9/P99997ljD56KOPPM/39xImiqoeVeRCDHoXXXQRa9asQVEULrnkEq8tSQsWLBi0U59vuOEGVq1aRUhICHPmzCE8PJwNGzZQV1fHzJkz+fTTT4dEC1tn3nvvPS688EKgtXjbhAkTOt0vPj6exx57bCBDC3iLFy/mlVde4cEHH+Tuu+/2dzh+9dBDD3HPPfeg0+k4/vjjSUpKYvv27RQWFhIaGsrbb7/taR0Yiv75z39y9dVX43Q6GTZsGDNmzCAkJIRt27ZRWlqKRqPhqaee4g9/+IO/Q+0327dv57rrrvM8PnjwIGazmZSUFE+rG7SOCzy6wO3f//53br75ZhRF4bTTTiMhIYEvvviC8vJyMjIy2LJli+8nWvg89RIB7bTTTut0xeZf/tx3333+DrVf/fvf/1ZPPfVUNTIyUg0NDVUnTpyoPvroo6rdbvd3aH710ksvdev3Y+TIkf4ONeBIi1J7H3/8sXruueeqsbGxakhIiJqamqouXrzY5+twBasdO3aoixcvVseMGaMaDAZVr9erI0eOVH/zm9+o33zzjb/D63cbN27s1r2moKCgw7Gffvqpes4556ixsbGqwWBQx40bp95xxx1qfX19v8QqLUpCCCGEEF7IYG4hhBBCCC8kURJCCCGE8EISJSGEEEIILyRREkIIIYTwQhIlIYQQQggvJFESQgghhPBCEiUhhBBCCC8kURJCCCGE8EISJSGEEEIILyRREkIIIYTwQhIlIYQQQggvJFESQgghhPBCEiUhhBBCCC8kURJCCCGE8EISJSGE6Ka//OUvKIqCXq/n22+/7XSfDz/8EI1Gg6IovP766wMcoRDC1xRVVVV/ByGEEMFAVVXOPvts1q9fz5gxY9ixYwcmk8mzvby8nClTplBVVcWVV17JK6+84sdohRC+IImSEEL0wJEjR5gyZQqVlZX85je/4Z///CfQPokaO3Ysubm5RERE+DlaIURfSdebEEL0QFJSEi+//LKna62t1egvf/kL69evJyQkhH/961+SJAkxSEiLkhBC9MItt9zC448/TkREBM888wy/+93vaGlp4a9//Su33nqrv8MTQviIJEpCCNELDoeDk046ie3bt3ueO/vss/nPf/6Doih+jEwI4UuSKAkhRC/t3r2bSZMmARAVFUV+fj5JSUl+jkoI4UsyRkkIIXrpueee8/y7vr6eHTt2+C8YIUS/kBYlIYTohffff5/zzz8fgMmTJ7Nz504SEhLYuXMniYmJfo5OCOEr0qIkhBA9VF5eztVXXw3A1VdfzebNmxk1ahSVlZVcddVVyPdPIQYPSZSEEKIH3G43V1xxBWazmXHjxvHEE08QFRXFG2+8gU6n4+OPP+bxxx/3d5hCCB+RREkIIXpgxYoVbNiwwVMvKTw8HICTTjqJ++67D4A777yz3Ww4IUTwkjFKQgjRTd9++y2zZs3yWi/J7XYzZ84cNm3axHHHHcf27ds9iZQQIjhJoiSEEN1gtVqZOnUqhw4d4qyzzuLjjz/utF5SSUkJU6ZMoaamhsWLF/PSSy/5IVohhK9IoiSEEEII4YWMURJCCCGE8EISJSGEEEIILyRREkIIIYTwQhIlIYQQQggvJFESQgghhPBCEiUhhBBCCC8kURJCCCGE8EISJSGEEEIILyRREkIIIYTwQhIlIYQQQggvJFESQgghhPBCEiUhhBBCCC8kURJCCCGE8EISJSGEEEIIL/4/PpeJjjPAlSkAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "training_data_size = 15\n", "sigma = 0.1 # Measurement noise\n", "n_resample = 3 # How many times a given x is tested for objective function y (generating statistics)\n", "\n", "# Generate training dataset\n", "features, objective_functions = 'x', 'y'\n", "starting_data = pd.DataFrame({features: np.linspace(-2,10,training_data_size)}) # sample points for initial model training\n", "df = run_experiment(starting_data[features], sigma, n_resample)\n", "# Use prep.combine_cut to combine repeated measurements\n", "combined_df = prep.combine_cut(df, objective_functions, features)\n", "\n", "%matplotlib inline\n", " \n", "# Initialize plot figure\n", "fig = plt.figure(figsize=(6,4))\n", "ax = fig.add_subplot(111)\n", "\n", "\n", "n_samples = 100\n", "X = np.linspace(-2, 10, n_samples)\n", "ax.errorbar(X,f(X),label='Ground truth', fmt='--', alpha=0.5,capsize=1)\n", "ax.errorbar(combined_df['x'],combined_df['y'],yerr=1.96*(combined_df['y_std_mean']),xerr=None,label='Training Data', fmt='o', c='black', alpha=0.5,capsize=1)\n", "\n", "ax.set_ylabel('y', fontsize=18)\n", "ax.set_xlabel('x', fontsize=18)\n", "ax.tick_params(labelsize=16)\n", "ax.legend()\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Start workflow and run autoML" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Small dataset (n_data < 500), using exclusively GPR models, can manually include other models with WorkFlow.add_regr('', '')\n", "Best performing model for objective y: iMatern_alpha\n", "CV average model score on test (unseen) data MSE: 0.1218300770445662, R2: 0.5670326395256643\n", "Model performance on validation set:\n", "No validation set defined\n" ] } ], "source": [ "wf = fit.WorkFlow(combined_df,\n", " [features],\n", " objective_functions,\n", " random_state=random_state,\n", " polynomial_degree=3,\n", " validation_holdout=0)\n", "wf.autoML()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Use activeLearn to find expected points of interest" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'Query By Committee (maximal uncert):'" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
xyy_std
07.0467520.8544030.114471
\n", "
" ], "text/plain": [ " x y y_std\n", "0 7.046752 0.854403 0.114471" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "'Expected Improvement'" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
xyy_std
02.0355871.2623290.074863
\n", "
" ], "text/plain": [ " x y y_std\n", "0 2.035587 1.262329 0.074863" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkoAAAGGCAYAAACE4a7LAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAADLeUlEQVR4nOzdeXxddZ34/9fZ736zL93SJV1S0pUiICiLKKiooOOGo6LgfMcFFXdcZnB0YMYZZ8TxNy4zCqOOy4CCooCKgILstIWWLnRLs683ufu9Z/39cbI0bdImTdIk7ef5eNxHm5uTcz83vb33fT6f9+f9ljzP8xAEQRAEQRCOIc/2AARBEARBEOYqESgJgiAIgiCMQwRKgiAIgiAI4xCBkiAIgiAIwjhEoCQIgiAIgjAOESgJgiAIgiCMQwRKgiAIgiAI4xCBkiAIgiAIwjjU2R7A6cZ1Xdrb24lGo0iSNNvDEQRBEAThKJ7nkU6nWbBgAbJ8/DkjEShNs/b2dhYvXjzbwxAEQRAE4QRaWlpYtGjRcY8RgdI0i0ajgP/Lj8ViszwaQRAEQRCOlkqlWLx48fBn9vGIQGmaDS23xWIxESgJgiAIwhw2kRQZkcwtCIIgCIIwDhEoCYIgCIIgjEMESoIgCIIgCOMQgZIgCIIgCMI4RKAkCIIgCIIwjnkTKO3du5f/+I//4Nprr2XdunWoqookSXz1q189qfPdfPPNSJJ03NuePXum+VkIgiAIgjCfzJvyAN/+9re57bbbpv28GzZsYOPGjWN+Lx6PT/vjCYIgCIIwf8ybQKmxsZFPfepTbNq0ic2bN3PLLbfwox/9aMrnveqqq7j55punPkBBEARBEE478yZQuv7660d9faLeLIIgCIIgCFMlog1BEARBEIRxzJsZpZmydetWPve5z5FIJIjH42zatIk3vOENE+r/IgjzTjENZhb0MBjiNS4IgnAiZ3ygdO+993LvvfeOui8ej/PNb36T97znPbM0KkGYAbYJ930aBprhFZ+G+ktme0SCIAhz3hm79LZixQpuueUWtm3bRiKRIJFI8Nhjj3HllVeSTCZ573vfy//+7/+e8DzFYpFUKjXqJghzVqQGajdB7YbZHokgCMK8cMYGSu9+97u56aab2LhxI6WlpZSWlnLBBRdw7733csMNNwBw4403Yprmcc9z6623Eo/Hh2+LFy8+FcMXhJOjGmBE/Nt4imlId/p/CoIgnOHO2EDpeG6++WYURaGnp4ennnrquMfedNNNJJPJ4VtLS8spGqUgzICh5bm73g8tz872aARBEGbdGZ+jNJaysjKqqqro6OigtbX1uMcahoFhGKdoZIJwCkRqIFgulucEQRAQgdKYHMchmUwCiN1vwplHNUaW6ARBEM5wYultDL/+9a/J5XJIksSWLVtmeziCIAiCIMyS0zpQ+ta3vsWaNWuO2ebf3NzMj3/8YwqFwjE/c8899wxXAX/Xu95FTU3NKRmrIAiCIAhzz7xZetu6dSsf+tCHhr8+cOAAAN/97nf5zW9+M3z/3XffTW1tLQC9vb3s3bv3mGAnkUjw7ne/mw9+8INs2rSJhQsXks/n2bVrF/v27QPgkksu4dvf/vZMPy1BEARBEOaweRMopVKpMXegtba2jkq4LhaLJzzX4sWL+exnP8szzzzD/v372bp1K6ZpUlFRwZVXXsk111zD29/+dtFPThAEQRDOcJLned5sD+J0kkqliMfjJJNJYrHYbA9HEEbYJjz6df/vr/gkqPrJHSMIgjDPTeazWkyZCIIgCIIgjEMESoIgCIIgCOMQgZIgCIIgCMI4RKAkCIIgCIIwjnmz600QhBMopsHMgh4GQ1SUFwRBmA5iRkkQTgcz3Mw2nU7T0dFBOp2e9nMLgiDMZSJQEoTTRaQGajdNezNb0zS54YYbeOc738mTTz45recWBEGY68TSmyCcLmawmW1tbS3l5eVs3rx52s8tCIIwl4lASRCEEzIMA8MwiEZF7pMgCGcWESgJwmnOclxsxwPHQS3mwc5jZwYgVA6AqkhoiliFFwRBGIsIlAThNJfImnSnigS6tlKx63eoToFiJkf/mndSqFhHVcygOhaY7WEKgiDMSSJQEoTTXFlYJyblUZ75JVKxB1fRCPVupWRnD9alN6OGVo+edXJcAGzLAdfBNB0c10ORpdl9IoIgCLNABEqCcJrTFBnNTuI5OfLhatRiEtkIIaVaUZ//HyhbTiK6ig5jBY4SwOhIoLgFCi3teKEKLNMkbzlEDPF2IQjCmUe88wnCmSBYhhcoQXJtzFA1WrAUgqV+Ycpcgnj2CSKlPTjxJTgDzyA7BdQ938Ld+G6KVQ0ENWX0+URxS0EQzhAig1MQTneuC5IM69+BqwZR7BxeSR1c+kV45adhzetRypahV9Vj7LoTxcnj6mG0Qj+BnT8j6OVGL7vNcHFLQRCEuUTMKAnC6a7zBdj/ICx5OdkFL0eyCxgX3wSRCv/7tev9W38TUmEAVwmgmoMVuPP9/u1okRoIlk97cUtBEIS5RswoCcLpzC7CoT+DY4Hn4ikGrhEfe7lscHlONZPIZhr6XoJA3F+iO9pQYcsZKG4pCIIwl4hASRBOZ4cf93OJQmWw4ARVtQMxWP8ObKMEyXVA1mHpBSIHSRCEM5pYehOE01W+H1qf8f++4lKQleMfD1C7nuyCl2MlDxFfsBLM3MyOURAEYY4TM0qCcLo68DC4DpQtg/L6Cf+YpxiY8WV4WhAGmnHTXQwMDNDV1UVTUxOu687goAVBEOYWMaMkCKej/sPQsxckCVa8yv9zEjzFgPJ6dm97krt+9QUefvh5bNsmnU6zds1qrl6UoKGuaoYGLwiCMHeIQEkQ5rHxKmpLA50oLkgLN6BGKk/q3LvTUb71y8fpTuYIhUJomkF5eTnbtm2n5elWPvrml9Pwiul8NoIgCHOPCJQEYR4b6uPm2UW07gwAVlcaSVuDUldNeWmM6pM4r+t63PPws/TmoWH9ZhI7DuEhE4vFKFnbwK4/7+Wev+xi9TWuWL8XBOG0JgIlQZjHysI6sYCGa2kMBDUASiojyJoBRFCVk+vP1tqTZM/eXhavuwApHMbj8PD3JEliUWWc3Ye7aW5uZmn9qul4KoIgCHOSuBgUhHlMU2SCukJQV1BliUChh2Cxd/g+TTm5/+LZfJFCoUA4HB7z++GATsG0SaczUxm+IAjCnCcCJUE4TcjFFIHEbnjmvyDdNaVzhYMGgUCAbDYLnkuIHDEvieW4WI5LKl9E11S0YIi86ZA3HSxH7IYTBOH0IwIlQTgddLxAtOVB9HQzdD4P6fYpnW5RZZw1q1fT0tKCZ2apoo+Yl2IglaEvU2Rva4LK8jgFvYT93Rn2d2dIZM1pejKCIAhzh8hREoT5rpCC525HMTN4io6khWH7T6FshV9t+yTIssRVb3oTre0d7D7QglWQiWkOaqqVtrRLdUmYN13QwMrq2GA+FCedDyUIgjCXiRklQZjv8gmkVCueouNqYShbOn4z20loaFjDRz/6UTZu3EhXXqY1kSPZ3cw5Z2/iA284l4YlVcO5UFPJhxIEQZjLxIySIMx3wTI8x0JyTFyjBFIdEKkau5ntJDU0NPDpT3+aRF8Pa4rbueyCc6i75E30PfvLqY9bEARhHhCXgIIw33kOUvkKPEVD8ly8UDlsvOakl92OJssy8ZIyAtX1LK0pRe7cPi3nFQRBmA/EjJIgzHeFJF7FKvI9bZixOoyLb4JIxbQ/TAdVeJKEPNCMZOXwtNC0P4YgCMJcIwIlQZgPimkws6CHwYiO/l7pUjjnbyh0d4GsHvv9QTu6d3Cw5XE8PP70wI24inLMMbLjcFHLS0hILO/ewYaFW4a/VyAAZSvw7DxSYifedD4/QRCEOUoESoIw19km3PdpGGiGV3wa6i859hhJ8oOk4/hL819I9+zAkOBnvY/Srxwb6miuRKVXjud5dLY8PipQAnAb3oisKLjtt0zpKQmCIMwXIkdJEOaDSA3UboLaDaPvL6TAdSZ0indtfB9/UUpIexINuFiudezNs/A8D13Wede6dx17ksFgTHKKyMWkP9MlCIJwGhMzSoIwH6iGfzMio+/ffS9ku2HlFSc8RUSPsLrxOv57+/+QIAXYcHTpIw9Mz6Cx8mXgBfC8MRbYOl4g3PYYqpmGR2TY/B5YsPFkn5kgCMKcJgIlQZivzCwkW8DzIFR+3EOTeYum3gwryq8kLz+H7QTQpAFsuRNPygIgeWFUt4YU9Zjqep7d30lFRTkLotrIiYppeP4n6NkOPFlFSrbC9p9A2fJp22UnCIIwl4ilN0GYr3pf8oOkaA0E4uMelsxb7GxL0jFQoCQY5vWrL8KjgEoFursCyQsjeWF0dwVQSXjxuVSGNMLZw3QmC+zsSFH0BhO/8wmkYgpHi4GkgGJMS3FLQRCEuUoESoIwX/W85P9ZuXrcQzzP43BflmzRprYkiCTByxddSkQrwaWA4pWj2/Xo9koUr4KQprFs5dV0L7qcQM1qauNBckWHtKvjefjFLQMlSNjgOZBs9gtbTkNxS0EQhLlIBEqCMEdZjkvedMhbDpbj+l9bjn9fLouTaPIPrBg/UMqaDomMSWlIJ1O06UoWyRdl3rDqXSheKXgKmrcMnaXgKbxx1TXknTDddpi0qwNQEtIooGEj+6UH1r8DR48juRYeMjReLZbdBEE4bYlASRDmqETWZH93hn1dafZ3Z0b9vXXfC+SLJoQrIDx+fpLtuNiuR8F26EwWANBVmTetvpygJuGRBhw8HIKaxdUNVxDQ/LeF3rSJ7bpoiozrSbhDmd+168ksvJBCyUqkpReAIYIkQRBOXyJQEoQ5qiysU18VYUVlhFhQIxbUWFEZob4qQp3bQlBXjrvsBqAqMjnTpq0/D0DYUKiKGiyrKOO9m9+IK6dx5F5cqZe/3ng1S8pKqIgY6KqM43n0pItYjotjm2QzadJpvxyApxjY4Vp/J17i4Iz/LgRBEGaLCJQEYY7SFJmgrhDUFVRZQpWl4a/1+otQll4AVWuPe468aeO4Hqm8TTyoURY2kCR/ZuitZ72NiFqDTRcBPc+lS98AgCRJlIY1JCBbdOhKW1xy/hZKIyF27NgxfG7HKPH/0t/kJ5ULgiCchkR5AEGYj6LV/m2Q53lkHRnblQgWbaKqTk+6SFeqSE08iCbL2I5L0XbQZJmC5ZDKy7xj/VX85/ab+NiFf0t5OErHQH74mLAuczhRIGIoXHreJrSXrSUSGanj5OhRvDVXIlWu8iuDC4IgnIZEoCQI81wyb9HUleZgKoLjSbQeHkBSsgQ0hYihsqIyzNoFMZp6M3QmC2Q8m9KiRm1JgC9e+h7WLJB4/6b3Y9naqGNKQmEWlgYoCxugBqgtLwPAtYr+A0uyP6OlGbP47AVBEGaWCJQEYT5xbdjzIJSvhPJ6kkWHnW1J0rkCIcVBkzzylkNbb5GQrvDy+nKqYgEAGhfEKVgujuvRuDBONKAiSRIfedlH/HNrxx6jyhIHe7Mk8xbJvEU8qB1ncIIgCKcfkaMkCPNJfxN0vAD7/4AHIzWSYgEM2SNhaxQsl8qoQUCXyZr2cBsSSZII6SrRgEbEUIdzlY509DEhQ6Ui4s8YtfXnsR139A+4Dhx+Arb/1G/ee6RiGtKdoh+cIAjzmgiUBGE+6dnr/1mxmqzlDtdIAug1NTK2X0G7KmZQVxamL22SNSfWNHc81TEDQ5NxXI+OwRIDwyQZ2rf5AVyyZeR+24T7Pg13vR9anp3S4wuCIMwmESgJwjyQyRfp6ktiduzy76hcNVwjSVdlBvIWGccPkmpiBrGAhqbI2K537CzQJEmSxKLSIAADOX8JLuvIJC2FjOnglS71Dzy6TECkBmo3Qe2GKT2+IAjCbBI5SoIwCyzHxXbG31KvKhKa4l/HmKbJ5793P262j3+6/jWsWLsBYotQLRdVlkgXLPqyJrg2FUqeiGwCISzH/76qTP16KKSrVER1mnpz/GlfP04qgjeYOL7QraTOdggkDh31JAz/ZkTGPqkgCMI8IAIlQZgFiaxJd6qI63kc7MkCsLwyjDyYN1QVM6geTMIGqC6NsHqhTU1Njd+yRJYJ6369o63NA8QCGnJVA5HsDn9mZ+Fm+nMmtSUBwroyLWMOqgrNiSwDmSIRR6bWMAnrCi35KoJpiyq6CeQHIFgyLY8nCIIwF4hASRBmQVlYJxbQcF0P0/aXxuorI8iyHyipyuhEa0NTqIu6GIYOlasAf0nM0BQ0WaY/axKuWM5A5UKC8VKSyTxhQ2VJWXjMpO3J8jyP5v4cQU1BjRikPZmiK1OuKYSCJaT0Sox8D0rvAezqjeA4qINLfrbl+EnfjJ4pEwRBmA9EoCQIs0BTZDQFXNcjoPkzPkFdGQ6UjmbILqYr46kBKKkD/PpJlu2xoioCnsfezjRZVyfmGdSWGCwpC0/bdv6h5rq18SDJLKSBPlOjxvOQAbV8OfnWThIte+iW6vHsIlp3BgCrK42k+jvnjp4pEwRBmOtEoCQI84AH7MlGed3at6LJCpbjDvdvW1oRoipiIMvSMTWSpsuRieMVYY1OycP2JNJFm1LNwClZht32NIFAgPqqCK6lMTAYpJVURpAHi1IePVMmCIIw14lASRDmOKlrB1dUdBBSbNSn/wO2XEuLVo/jegR1mZpYAM/zE66BcWskTYWqyKiyhGm76JJEVLHptzUyBZvSCOSMKnrX/g01K6oI6gqupJAZnB0L6gqyNj15UoIgCKeaSBYQhLmskEJ5/idEVJuUrSHl+sg980Py6QEkCRaVhqY9KBpLWFcoi+j05/yikhHVzz/KWy6W49KftyiLh6YtcVyYPpbjkjedcW/WFMtHCMLpTswoCcJclk8g9R8kqlgUPI2cUY2d6kUpDlBVXTWc3zTTJEmirjxMumDTkcojeaBLHkXb5UBPmgUloZHE8WIGZNHqZK6Y7A5LQRBGE4GSIMxhSSmGVDAxNJmCF6S7ownCVXiBEsrC+ikdSzyo0bgwTlOXx0FHwfQkXNMhHg7QuDBOXPPgqe9BPgEv+3+ndGzC+Ca7w3IqJlIfDJhwDTFBmAtEoCQIc1Qyb7GntYcqfRGH+3fTnMvhVKik615HhW3MSpPaeFCjcUGUWLifou0wUKESioUwVBlUBSQJPA/6D5/ScQnjm+wOy6mYyOwVIGa4hHlFBErC6aOYBjMLehiM6GyPZko8z+NwX5am5x7i7ocO8afnLLotlUVdGdb2PcnZToyS0HoaF8RPSY7SkSTVwKu7kNKOJ9DNw2SlMpJ5y/8QLlsG2V6/95sw4yZT4f1UmOjs1ama4RKE6SACJeH0MNSEdaAZXvFpqL9k7ONOQTA1HR9eWdNh2/M7+f3//YSBRA+WUYIeiFBeUc2BXTvobm9DU65jWcV5RIxT/984X7GOYslyqivKyeb8HnDVsQCULYeWZ6D/kD+zdIqDuDPNXMs/mujs1ama4RKE6SACJeH0EamBYPn4TVgnGkxN0WQ/vGQzg2xnoShBMAaAadk8+sA9pPr7WFpbRktTHsWDitI4VZUV7H9pD3/+w/284aItMAuBkqtHcPUIsViU9nwa03bJmTah+GKQVSQzg2Tn8bTQKR/bmeRU5h+dSnNtpkw4s4lASTh9TKQJ64mCqWkwqQ8vx6T28S+hZ1rhss9B/aUAdLS10rHveapLI5hqBNsz0d0CSBKSJFFZvYCWQ/voaGulbHX9jD2XE5FliXhQYyBnMZCzCJUEoWQJ9L6EWhzAEoHSjDqV+UejTGBmdqwLgImaazNlwplNBErCmeUUdLSf7IeXHarGCZQRqhkJ3pxiDhOdTHgJuhFAdTuRGLnCljQDybVwirkZex4TFQ/5gVIyb1EbDyCVLYPel1DMJBYLZnt4wnSbyMzsOBcAE3W6zpQJ85MIlARhlnmKjqfoo4K3SCSKEY6SlsPoShB5MEiyHZe8bWObBcpiYWKxyV2pz4SooaLIErbjkSnaRMvr8Wo6sHq7Z3towgSNWuoqpsHKgjYyW3TMUtcEZmbHugA42nizTrM2UyYIYxCLvIIwB+klVayoX4WZ6qUiouFKCrakYzkeYUOhONDD+sazWLJkyWwPFUmSiIf8MgUDOQtCZbDqCpxA2SyPTJioRNZkf3eGfe19dP/8o2R+/B7adz7K/u4M+7szJLLm6B8YmpU9zsysp+i4Wnj8YwZnnRY/9BFoe3Yan40gTC8xoyQIc0y6YFE8+ARXnNdApq+dRNsh5HwCTVUJS0V6W1pZUFPFVVddhSxPz7XO0IyC63oULAeAvOmMWuo4XvJsSVAjkTFJFSxcd/wkXGFuGl7qsjSSJQuwY5XUrDkXOewHOaoijcw6OQ7qYNsT23LAdYaPmWyC9URmncYjEr6FU0UESsK8kE6nyWQyRCIRotH5XSPpeGzHpbU3RWnHk1wQcFn6t+/hzt/8kb3bniCfy5NODnD22Wdz1VVX0dDQMG2PO5Q8CwwvdRzszQ5//0TJs2FDRVMlLNsjXbSJKi5Kvg+lmPCXcjRj2sYqTL/hpS5JIaMZgEEwEh/VzLgrVaA7VcSzi2jdGQCsrjSS6v/bnkyC9VjLzhMlEr6FU0UESsKcZ5omN9xwA01NTXzhC1/g1a9+9WwPaca0DxRQBprQJJtIvJKGs1/Bp1ZvoOW5BzAth89/8QusXLN2UjNJE5ktGppRGM9EkmfjQY3etEkyZxFNPkes+fdIrg1//DKccx0s2DjhMQszYyo70Y6cdRoYrAhfUhlBHgyCT3WCtUj4Fk4VESgJ80JtbS3l5eVs3rx5tocyrTzPI+vI2K5EsT/HQFEinjxALKAhV64GSUKWZeJh/8q4rq5u0sttE50tmmp/3ZKgTm/aJJNK4L14N5IHnqwhDTTD9p/4xSgDs598fsaa4k60UbNOg8FIUFdGzTrNBJHwLcw2ESgJ84JhGBiGcVotuyXzFk1daQ6mIhQcmeKeXuJhg1dmDqLpMlSsnJbHmY7ZookI6gqGJuOm+7Gz/bh6DMXOgaJDvt+/iUBpVk0lJ2hWTDG4E4TpIAIlQZgFybzFzrYk6VyBkOyQdxR0TcEdaCabyRAqLyUYXzwtjzV05X0qlAQ1eoxS8mqUgGeD50CqFUqXQrD01AziNDBTicpTyQmaLfMuuBNOOyJQEoRTbKjhbbZoUxsL0OQqmJ5EiabQGOjGcjy6tMXUSRLzbREhHtLo0qP0LLua2o7tGGYaFANWXS5mkyZhthKVj1wKDhZtoqp+ypsuHy1temTzLrrpEZ/VkQhnKhEoCcLRZrhxbtZ0SGRMSkM6Bdui3/b/G1ZENLRcHlVTaNWWUGE6s9LwdioMVSGoy2TLG2mveiU1qReILlkPumhlMhmzkah85FKw40m0Hh6gLFqk0NNCLtlHVVUVa9eunfbHPR7TNPn89+6npXuAm2su4zWvfd0pfXxBABEoCcJop6Bxru242K6HIkN7ogBAWHGJBjRSy6/EW3gRmbyEPVir5nimWv9oJsSDOplsng47imEswkBDH2hBqjulw5jXZiJR+XizRaOWghUHTfLobj7A/9z3Wx564Dfkk71s2rCeSy65hDdd+XoqpuVZTkx1aYTSaJBNmzae1M+LekvCVIlASZjXZqII3kw3zlUVGUWGlv48tuuhSR7l2kjl46IcQlFt1CPGbVo2pu2QTqcJhEdmuaZa/2gmyBIc6M3Slg7TrdSxoOxSoqXLqctbxIPjJ5ULU3O8rf9jzRZVxC3qysPEAuqopeB+2eNwWze/eeAH9CcSRErKKCsro76+nm3btnG46RBvaTRYuqDqlCzP6ZqKrqnH3chxvOcu6i0JUyUCJWFeG3oTPF4RvLKwPrlgaoYb54Z1BZDoShapDMtU6iayBJJdBFWnP2dSWxIYPA50XefiTfU8vrOJHTt3cmnNSKPZU7WjbaKSeYs9nWkyRZuA7KKqKmrJYjrTFmkzSePCuAiWJimbSZPLZklHJeLxcfK8jrM7bKzZorCu0JkskC7Y1JWF6BgooCsy/TmTnqLCH57eRV9GZcWq1QxknsdFJlZSRqy0nK3bn+eHCYc3vnYpFUcEXLP273qCnXGi3pIwVSJQEuYHuwiO5ecP6eXDd0+kCN5EgqmJXFFO1xR+qmAT0BRCuoLruXiAbOcoe+E79AcWE175ZpaUhUddpW9cVs6q6iCLz1o16lynckfbiRyZpL68IswB2SPnKv5zDep0JPM0J7I0LojPeoLwbJrM68g0TW750mdob23mli//Ha95zWvG/bmxdocd+W9SEzVodWVynoSVt3Akhd0dKQ71ZMhbDmVhA1ybw90pDrX1UrpsPcmCje0puJJEpmiTtz0i5dV0HtqGneoZFXDNZhB8vJ1xot6SMFUiUBLmPtuEfX+AQhJaXgNrRz4sJlIEb7oqCk/HFH7RdmjtzxExVM5ZVkqhUOTgiwp6PokVDFISVFmyqGT0B46iw8rXUNnxBNF8M7B0QuM91Y5MUlexAbA8CTuforx/O+HcAM3KG8jOwyT16TTZ11FldTUlpWVs2jR+sVXP88gQwCaIQYCo5yFJElnToTddRJKgKZFnwPJfV27RQVZkIoZKMmciyxKyDGFdRbVyuLZFMBjC88BGwvUk2pMFXA8ioRA9pk2uYM6ZIHg+lj0Q5o8z991KmF/0CGhBqF0/7iGSU0SyC8f0FpuuisJTncJ3XY/mvhyuCyFDYXlFGNfSicYyhMw2SuON6Cs3Io1xVZ6vWEexZDmVtTWTGvOpNJSkrqsyOBIhxSXnyGQsiUW9O/BcD6+QwnZKZnuos2qyryNdN9D18Yutjpd/tKQ0RG/W5HBfjlhQA9dDxSOoupRHNAw9gCzBQFYjGtTIFG2qwwqJqELUkDGwMAIGkgSq5+I4Lrqmkkyl8BSdpBQlWbAojeiUhnT60uYZHwQLpyfxihbmB1kBlPGvGDteINz+OLJThEduhc3vnvbeYiczhX9kkmlbXqVguSiyxJKyEJIkIUkSEQqE3CS6KiONU43b1SO4emRGyhVM1Il22HmAKkuYtosuQUh2yDkyaUfBCtUgpduJ5ltQlaWjTzzD5RjmmulcChor/yikybzUmWZna5KaeBAJ8HCpjuo4gSJIUBrUkVWVguVgaApLK8Ic6s3SkcoTLy1l8YJKDrW1UVNXj46DLjtEAxqGptB9uJNFtdWUlpXSkzbJWBLlYR3b9Sa0U1MQ5hsRKAnzXyEFL/wMtZjE0UJIuT6/t1jJElA0fzZqcDlgvFmnmeDZRWKPfhk500Hz+R8jVbkFWZZYUh4alcukFPv946O1SEcUZZxrW/9PtMOuMqpTFtHpTBaoDisEFQcsjaLtkY0sRk20UGm3DyepA6ekHMPp6ujCpQnJI2UrpFImiqySKVpkiiYrqiLkTJtYQKX/qFhsaONAbTxA2FBp6vI46Gms27CR7hc66G45gJdPYssy6WQ/nQMJqqsqed26ckoNB1eSKFouTX1ZVEUeXj6ci4UrBeFkzZtAae/evfz+97/nueee47nnnmP37t04jsNXvvIVvvjFL570eR988EH+7d/+jaeffppsNktdXR1vectbuOmmm4hExHr3bBoKFEzLwfH85Ne85eCYR+1WyyeQMt1ITgFJC0G0FnK90LsfDvzRn40yopDtI3bofiQ8eESGze8ZnnWa7jIDQ8shh5wV5NTlFAdqKSHDhsUlxyxNqLluP3iLLRx1/1zb+j+RHXZBXSVdsOlI5ZE8CV3yKNou++1K1ikylU7nsdXGZ7gcw+nqyJywrGnSUjBwkYh6HoYisbwqjCLBiqrI8GyR5Epokh94J7N5woY6vHEgHtRoXBAlGstwViTMBZd/jD/8/nc88sgjtLW10dvZypoN5/CGK19PRfsf0dw+jIhNrxPgUG+RyqhO20COVEGjL5mleUBDsfN07G+ltKJ6dnfGCcIUzJtA6dvf/ja33XbbtJ7z3//93/nEJz6BJEm84hWvoLq6mkcffZRbbrmFX/ziFzz22GNUVJzK0mrCkYYCBTOXJpH16wzt70oPF3keDhSCZXhWFskxkawcpDsgUgWq5s8kuQ6ku5AOPoJi+vkVUu++UR3tJ7sz7kQ1a4aWQ3RVIuUYGMEgOdOmI1mgPGKMfGC0PkOw70Uk14I9v4FY7XDwNte2/k9kh108KNO4MO7PTDgKRU/CMx3KShdQZcUIeCakO/3nOWSGyzHMd5Jj+q+PYhoCI7OgQzlhecumJ1nARUKVPKqjOrFICNeF3kyRkK6M+jdxPAlMh9qSMEvKRgcvkiQRVlxQoH5DI2dv3sQb3/hGuru7CcXLsaK1KB3bWNj7Jww3D9tSFBa/mWWVZ1Ea0hnI2TzXNEDpwAs0DjxE0M3B7gSHFr+ZnYUNojyEMC/Nm0CpsbGRT33qU2zatInNmzdzyy238KMf/eikz7dt2zY++clPoigK9957L6997WsByOVyvPGNb+SPf/wjf/u3f8tdd901XU9BmKShQKGQ9TBU/xN6eWWEQNj/QB0OFBwLYgvxlJ0gK3ihcqSN1/gBR816/wOmexd07cTOplCcAlgFyHYPd7Sf1M6449RtGbUdO2awx9ZxJYmoobC4Mk53qjCyO6iYhp134Sk6rhpGNbOjgre5tPV/Mo6cmcg7MqnaKOFgAM2rg8QB6D80OlASxuVaBdw9D2Dn0xx+egsrX/VuZNmf2VQVmWTBwrJdNAkiikOFZhENaEhIWI6DKkuog7vbhv5NbFeipq6EaDg05nLYkcvTcqSCxsZG/xuOTartRayWO8nZeZJKnJJMD+t2fg3lrCsxMiFe6s6yfCDFiv5HkHNtWHqMYHGAle2/5vlIHc0JdVbLQxzvAkcQxjNvAqXrr79+1NdDbxYn69Zbb8XzPN73vvcNB0kAoVCI73//+yxfvpxf/OIX7NmzhzVr1kzpsYSTMxQoSJbCUJ5rUFcI6EdFD4f/ghSpIl++DjO6mPKLb4LI4EygrECwBKob8UrqoPcwLgG/QWu4arij/WR3xo1Xt+XI5ZDeTJ6C6ze2rYkFUCRp9O6gfAKpmMYKVoCsYcRqIds7HLzNZ0MzE2HFRQ1rmJ5ENryEWD7hlzsQTmj37t3c+fOf8ejeXmzHIf/fv+GsZ9q4+uqrWb16Db3pAjKQytusqjRAt0b9/NGFS4+cLYoY6tjBypGbIh7+R1j1Gj/Hb6AZUu3Esr24xVaUUICQ5BGpWoze1oKUacUsxgjk89QpSYLmAHlJJu9pFLVK4sUBKqQU0q7HyVsNhKpWgBo6pTmDJypMKQjjmTeB0nQyTZPf/va3AFxzzTXHfL+uro4LLriARx99lLvvvpubbrrpVA9RmKhsrz9bBJjxZbjaODunAjFY/w7cpmeQnSJe6VJ/1ukkA5Lx6rYMLYckCybJvP/BVaFb/pZ5/IBseHdQsAwvUIJiZnD0CKQGlwwHg7fTRSyg0Zv36C85i1j9ebM9nHlh9+7dfPOb36SzswPdCBBSVWLllWzdupWmw828+d3Xs3DpSmpLgoR0hYJtHzf/aEIGN0XIdg7JNpGan4CO52HZK/2lUYBwOQTjGOmDOHoEvdCDVLMeNl5DwdFpb0tSrlsErSRex4sU5SCRTAepUAWKWyTUvxf5pWY4pEC+n9ih34EEPCzB2e+d9p2qRxvvAkf0gxOO54wMlF566SVyuRwAW7ZsGfOYLVu28Oijj7Jt27ZTOTRhspoeA8/DK1+Jm9l9/GNr15Nd8HIku4AxNOtUSIGVg+j01CdSFZlUwcIcXA6p0CzCijP8fctx/eUQWYKXfge1G3CbnkaxciNLhvN8NulosYBKb94iXXRxXQ9ZlmamR98cMB0fuK7rcvfdd9Pe2U3VwuUcONSEjUTG1SlduJwD+/Zy/2/u5QMfuZF1C+PYrkdT18AJ84+O/6COvymiMICjR9GcPr9cg2NBySJYsBlK6iBUBi3P4v76k/5rNlyJtPndULsBuWhjZhKkDBXprPcR6vsqcSdHwVjB4YVXkS5GWVB1LgvCA5A6jNTy9EjOYOuzYOfhtV+b0df/eBc4oh+ccDxnZKB06NAhAEpKSsYt4rZ48eJRxwpzkGNBptv/e93L4fAJAiXAUww8xfBnnZJtsOP/QDHgnOtGrpqnIFuwkPCXQ9ZUB3BVZ9T3h5dDMk3Qtx8kmWzNeUiuNRK8nWaCuoKm2li2R7poEzdk+vt66TIDU24rM9dMxwduc3MzO17cRaC0iqzpIHsOEn4w2ZUqokUr2P/SXrR8H2GjBGDC+UfHKKah+Sno3Qvr3zE8w2nrUfRIFcQXwbq3jQ5exrrgwO9hOFQeQi9vpKv8IjQ3j/6y68nZQbr6C3jxzaxcXUUg14zXuZNiNoPkFJFlDcXMIs3SsrPoBycczxkZKKXTaQDC4fC4xwyVBkilUsc9V7FYpFgsDn99ouOFaaRocM71kGyGSPXkfz5cAWoA8gN+i5SGK6c0nO50ge60SU08iK7K5Exn7OWQ0hDS7sf9H6rdiJd8HA9O62KLsYBGX8Yk132IePO9lGsRopuvn3JbmblmOj5wU6kUfckMFQsriRgyEi4OCgXbJaireF6AYraPYn6kVMSE8o84IlE72QaJff7SmuuX37D6mvHWvhXn0DPITgErXIPb+E48OYzquKNmwkZdcBwxhrry8Eh5CCmAphpIWhjV8SgLa1THgrQk8iSkEFG1ioLdQk4up9zTQSpDk6LEAQ4/AYE4hBdhyA5hxfaDuvDM/B8R/eCE4zkjA6XpdOutt/LlL395todx5pJlKF0KVvGEhx5DNWDNlbD9f6FzB1SshJKlYx6aTqfJZDJEIhFChnZMMb2+rElX0h/Diqowa9XY+Msh+VZItYOswqItsOfxk3/+80Qs6AdKA1KcGttEsRME3SyubkyprcxcMx0fuIoR8hPerSIYQRwUbFnHA3RFRpIdioruHzcZHS8QaXsUtdAPv3ze3xEarYH4IvorzqFDqkFSLKoW+rNFPQ034Mjl0J2Z8AxfPKiNWYpgYVmYJaUhMqZDc1+O57styuKXc5a6i7ibxo5v4tCCN+L2ujR6vcSbHgXXQcmneHN1CwDqY1+DLdfOeB6TIBztjAyUhpbbstnsuMdkMv5yQCx2/Gngm266iU984hPDX6dSqeFlO2F6uK5LMlvAtBwOHz7MyjVrkfsP+TkTyhRfwiWLYfG50Pwk7L0fNr37mENM0+SGG26gqamJGz/zeZau3TSqr5au51FkiYihUhUzqIr6HyjjLofsfcI/ce0Gf0fRaWqk/k+GcLgMWQabAGaoGiPXBf1NULF6toc559QuXMSiZfW07HsRPbAMW/Z3CQZUmbChsP9wB3UrG6lduGjiJy2kYOsdaJl2PyfIzPq//5f9DVSvJe56hB0PHBM1FgWiRBbVDu9QnMwM35HlIY5+7ceCHge7M+Qth3x5I8/EXk2tkiR29geIhMrpSOZpSUvElpyP1PQYcvPjLAjkSVkqUqp9VPkMQThVzshAaenSpQAMDAyQTqfHzFNqaWkZdex4DMPAMGZ4W+sZbGiL9F92NOE4Lv/w1X/krJXLuHpZnobVK/2lN3WK282XvsKv75Pp8ROsPW+45cmQ2tpaouXV6DUr6EiN9NVyXZdd7SlCusLLlpeOuuoeczkk2Qr9h0GSYcm5Uxv3XKboZBZdRLjjCeh+EWn5K4kFNAZyFpnQ4sFA6ZAIlMagayoXvfp13N3dzt49e7FME10Bt5jjYMsh4qXlvPLVr0XXJvH2nU8g2QVso9R/aS86F5yCn9QsSWiK5NfsshUYXGLTNAXUkRm+4ZY6loPtDlbKH8yhgtGJ6uMtBWZNB8f1WFkdYSCdZ4Agh5wAa+QwBvhlNXI22WXnE4lUQdNj5GwVVQZyfeD5SeciUBJOpfm3rWQarF69mlDIn7Z+9tlnxzxm6P7NmzefsnEJow1tkd7+/PMEDZXSaJDysjK2PfEw3/zl4+zudaceJIE/K9XwRpAVpMQB1ELvMYfouoEXKsdTDGpjAQzZI+8qJAsOlVEDVZEo2i6eN/6OJ8DPvQCoafRzME5j+Yp1JBreBbUbAYarjA8EBmdCEof8oPQMJpsZ1FyXn38zKKwrbFh/Fi97/buoW7UWq5gnl+wnOZCgYd1GXvu2a9m8sXF0z7zxdL3on3uwFAWyjK3HwUz7ZSgmUYoikTXZ353hUG92eOfmod4s+7sz7O/ODFfPP56h8hmlIZ1FJQFUycNGonWgQNF2RpfPiNbila8k7agUHNkv0FpI+kvWgnAKnZGvOF3Xef3rX8+dd97JT37yEy65ZHQjzsOHD/P4437eyNVXXz0bQzzjDW2R7u3tpWHNGjr3+WUaYgGF0poguw73cM/WLla/3p1y8VHAr1+07CK8bC/OwJ+Qi8lRRfBsZAqolIT8D/uMrdBraUSBaFBlUSA4UkjSOM5/q6UX+LNVS86f+pjnOFeP4OqR4YTfaEBFkiBnVGNJKpqVH9m1eCYapwCiJPlByJLlqyitXoDZcwDXzPOuD3+CaPVSokH9xPWRrALs+x107YKyZbD+7cN1xBT75EpRDLfUcUzUsH+BEquKTGp5TlVkVFnCtF10VabWKNJVNHA8j9aBPOVhbbiaOEYMp/FtpB+8l5Bi41asRG58s78TbwYdXb37yPzE8XZJC6e303pG6Vvf+hZr1qzhPe95zzHf+9znPockSdx+++088MADw/fncjmuu+46HMfhLW95i6jKPUMsxyVvOuPeDhxqYs+ePSxevHjUB4KUakGSJBbVLWP3gcM0NzdP36CWnAvhKsIdTxLpeAIeuRXat+O6LomBAXq7Ouhqa6atP0ev5QdMsYBKdTQw+kr4eGILYN1f+fVozjDyYB4XskIuNNgAeODMLr9hh6rJV6yHIwogdqUKeJ5EfXWEdQtihKIlhCoWEq9YyMKy0Lj90iSn6Af4Pbvh2e/7QZIk+4GF5w1v68/Ung8X3zTppGhNkQnqCkFNGUxalwlqin+frkyo7tVQCYH+nD/7pEhQYxQJajKeC3s6MgR1ZXi2zKtexwO9tTzQW4t9+b/AureMnGyg2Z+hnc5ZycHgdfFDH4G2Z4fzE9/5znfy5JNPTt/jCPPKvJlR2rp1Kx/60IeGvz5w4AAA3/3ud/nNb34zfP/dd99Nba3fR6q3t5e9e/dSU3NsMcHNmzfz9a9/nU984hO87nWv46KLLqKqqopHH32Ujo4OVq9ezXe+850ZflZnrhPVm+nrTlAoFAiHw3iODYAhO36dFUUlXLOStkMtw6UepkUhBS/8FM8qkFcjqL0H2PPrb/KLjgU89vBDpC3obTnAkmUraKxVOWtxBZVR/2p6uJDkPCySeCrFghrpgs1AvIF4VZ2fkM9Dsz2sWXN0AcRk3qI75e+eXF0dJeBplBfb8CSZzUtilJeO0ydtsPWIWkjAvc9D7XqoaoCGN0B8MCh1xt7WfyodU0JgsHxGWUhjX1+RoKbgAamCPRwMFl2FoquMHrNdhF2/gmIG+g4gmVkkz56WVihHV++ura2lvLxcpGGcweZNoJRKpXjqqaeOub+1tZXW1tbhr4+saXQiN954I+vWrePrX/86Tz/9NNlsliVLlnDTTTdx0003iWnWGXSiejOBYhmBQIBsNksoYKBIHtV6AVwbYrVkTZdAIDCt/0ap/m7sRC+9bhmhXD/7OlJ87y/b6A6vRguGicQi6OE4h1/aRWpPHwuvuICl+OM9uq/WMfb93i9HsOS8MzoRNRrw33JS4eWYNVFUzzrBT5y+PM8bVWZC8xRaEn7HgPKITmlYp5AtonkmeMepj1RIwbYfomc78SQJycpCuh1e9SW/VtgcM2YJActlw+I4siTheRItiRxuSZCA52FJOp4kkynaGCHP/x0our8JY/8fkFqeJN70IJ4ShEcU2PzuKZUQGBW8eiMbdsTnwZlr3gRKF1988YkTZY9y8803c/PNNx/3mMsuu4zLLrtsCiMTTsaJ6s2sWLaUNWvWsG3bNlbVRlkZShNVLUi24pXX09rayubNm1myZMm0jCeZt9jdr1BLmIibBCQeeL6D3oEC0QUlSBkI4FFVXoJSHqd5WzOPPvsiq95kkc46x+2rJdkF6HzBXwapaYRAbFI7iE4nmiITMhRyRYd0waL0DN0wmsxbNHWlh8tMHD7UT8GRqYga1MQD1MYnUZU8n0AqpnBlzc9/W7jZb0liZudkoATjlxAAaO3PM5Cz2NOZJp/N0GcsxJVktjanqMn5M1LxoOYHQ0YM71cfRraLeJ6H1LtPlBAQpt28CZSEM4ssy1x99dW0HDrAnmf/TMCyMR0NbIP2rU9QsfZCrrrqqmlJ5PY8j8N9WVJekLKz3obS9xIdPf282CdTXl2Lkm9DR0WVZJaUBkkVLEpKS9nX2suhw82sO2vNcftqadl2pJIIlK/wc5QYWXqUHJPKweDwUG8WT/FnWOZrC4+JiAZUP1DKpClNtaHke3GCc/MDfSYk8xY725Kkc4NlJvBIFm16si5Z02Z1dWTijWzB39EWqsBTDRwtgm4V5kVz5fFKCCwuC5EzUxxoyZDKZnEkmaCTIWTIdCYLpAv2SJ6WLCNFa3AGupCdAlhZ6H0JZqAViuW4OKYz7vdP14sbQQRKwhzW0NDAR697J7/4xjb+vL2Tgm0hxTQ21wW56v3X0NDQMKnzHb3UEVV1JEkiazokMialIZ1c+CwOllzKnkQLbaGDrC7RCTgDLKabPsoxNIUqTUaJwsF0gWVxicYF4+SN4Bdd1PI9wKpRO92mYwfRfBULaHQli1g9B/G6foeeOkxR1qYlv2SuGwrKs0Wb2liAftkjYWl4nh8cqwq0JfOUhvUTB0tdu8ApwoJNsP4dOE3PIDvFY3a0zbfZS8/zKNouuiZRGlJJZoukikUSXe3ULV9FV7pIcyLr/78LluEFy/BkBUcKoRXTftmNGSi9kciaDBRE49wzkQiUhDmtYePL+MxbzmFzuJOOnMLFL1/JsmX1yBvPmdR5jl7qaD08QEXcoq48jOd5FGwHT/JIZgoMOBHscBV6qIdCpIKo6SFLPVTTA3YRSdWxzCJRQ6a6rOS4H2ha+jCSlcMLViCVjCwTDi09Hq/A3+kqoCkYmkwxWof5UieBgZfQc53+LsMp5pfMdUcG5eCQsRVStkIUqI4FUGXpmDITpmVj2g7pdJrAUK+zlqdh/x/9pbZw1biNamH+zV4O/Y7ob+MP997Ni1ufwXEcbvvXf6bxrEYuuOy16MoK/3cUiA2XPZCdIl71WUjnfxiCJdM+rrKwTkVMNM49E4lASZjbAjGk+ALOqnBZ6sjU1dUhb37XpKbVj1nqkDzCukJnskBfpkgsqNKZLBDQFDTJQ8FjeVWUxobV7N31IsG6ekwOE6YIdgFP0ejoTbN+Re3xc6SaHiPc9ZzfxqP9Ob8B6WkcBExGLKDRl7X8HUuOhSspSLm+OZ1fMjQzM56JzMwMFVzUVRmraNE3WGaiPKwRMVQc1xtVZkLXdS7eVM/jO5vYsXMnl1bXwoE/Qssz/gkXbvGXc21z3B1t82320nZcDuzbywP/dzuJnm5CmgLhCEa0hG3bt9He1soVb38fG4cq2x8ZJL72iyNBoudB106oOsvvCTlFmiKj64ponHsGEoGSMLdZBZAk9ucibE+WsvzCz6CWLZjwj4+11AHgeOB6sK8rS0VEIxZQSRdtFpUZ2IEiSHDp5a+nq7OTvfv2ky9EKNUjlBah/eBuymMhrjh39fg5UoUUbPsRkmPiamFUx57TQcCpFguqDBT7cTwJRdGRB3czku2dkfyS6XCikhZHz8wcXbgQRgouFm2H7nQRDwjI3uAM09hlJjbU11K/sJy6jRtg96/9JTeAFZf4fQpPsEQ332YvZQmeeOh+Er19LF+5ikTHIRxsQuEIkWic9kP7ePLhB3jrZSMtgMYMEg8+DM1PQe8+WPumWXgmwuli7ixMC8JY+g8heS4pS6fDDB23/ovkmMhW1q+tMujIpQ4Pj4yt0Fow6EgWyJsOsaCK5XpsWFxCfVWUjOlQ9CRcDxbUreA1b72WNY3rKeYLdCTSJBIJNm/axAcub2RN5WBejef5xe/23g+tz/kPnE8gSRKOFsEOVkK81g8A8v0z/RubF0K6ihQqoxBeiGNbuGYWO9GMFyyZs0nIZWGd+qoI9ZURFpUGWVQapL4y4t9XFaEsfEQ7naMKFw4ZKrh4qDdL0XaRgQp9pPVHf86kPKqPKjMR1aE2aBHdf89IEcmGN/ilJiaT9D1P9HW20dVykHhlzfCytoJDRFeQJYlASRVNB/bR29F6/BPFFoGsQM9e2HEXOGduKQphasSMkjC39e4DoNs8fg6FJ2t0V78SveMZjNadRFe/EkmShpc6TNulI50nUZSQXBPFKRKPBokYQZJ5i/KIQUU0MLq2i+lw7tnreeMrNnGLl6dYLHLjjTey3NqHc9//QrMCv+7yd7PpfsFAwhX+9uyhJFM1gCcBqY55sRPpVPKMKM/XvJmG5mcIWTlSOYf02jdS5QWZi13wTlTS4mhHFy4Ef6dXZcTgBdtjoFCkUrWR8ShYDsls/tgyE0OFJItJpOQOWHQ2nP9hf2aS+ZeoPRGZTAbVs4lGI2QKNh4yEi5Dv+ZwJESmK8H+tl6WLVs2/okqV/lV8Hf+AhIHYeddSFbOXwo/AzYOCNNHBErC3OW6kPArsPea47+pJfMWTb0Ztll1eCXV1JllVLQlqSsPU7AculMFZFlCUxTssgYqM7uoow0lUkPBcoaXOiKGOmZtF8uyKCkpAWBpbTny/bfg2UU8RUdKHPT7la28DGo3QPVZ/qCOSDJVrBP01rKL/tVuMQPq/G5tMvzB7fof/jD4wX1EwqumyCTzFq2JHAeUFVjhV3CWs5v80ito0lbS3ZY8tk1HMe3XBdLDs1ZVerKOrroN4LgeyYLFiqoI6axM9jAkLA1Mh9qS8OgyE4UUvPAzZDuPFShD10P+UnRo/iZqT0Q0GiUaDhGWbWxDxZUUPFQsx6M0bBCRLbxQEFcN0JzIsTAij7mbFfADyvXvgB3/h9T2LPGmh3HV0JQ3Doy1rCqcvkSgJMy6cd90Um1gFfDUAAP22DWKhhO18xaBUARNiRGOBDncl+WlzjQLSoIEdYXejMnKqjDSonokdxFSZTlwbEXt8Wq7DMsnkGQFO1iBbOcgvsSf3l//dn9m6UjH2Yk0TNWh/lXQ9Bh0vwjLXnHSv8e5YOiDGxiedTnYmx3+flXMoCpqcLgvi+N6VEcD9ATreEmvoqL+9dSqOh3J/Mj2b0kC24T7Pu0vb77i01B/yZiPPR+0D+SxbI+ysM7ZC0O07h8dlI96vaU7kPIJHD0Csua3JTkqh2vWErUnENxLjjk4e5OZ1OzNkiVLhovNrqxfQdBOgSRRGzcIGDq7Ww6wacMGFixaTPtAgV0tOXKpCN5Ru1mHA86SxdDwJry7/xbZygEyUrb35HMGx2lmLJy+RKAkzK7jven0+ctuXslSPI59wx+VqF0SJG+5WI5Lf87EdiCRtVAViYbaKN3pIrmigykH8NQQBSlAMjnGUsc4isUipmmScTRKAiUAWKEqdCPiL6mFK8f8uQn11lqwESpW+rMl89zwB/c4VGWkblVZ2CDh2gzIOklkhsLI0pB+zBZ5IjUQLPdn7eapgZzJQM5CkmBxWRBVcsYPyu0iHHgYL9eHYqZx9OiYy7ezkqg9keBe0cksuohwxxP+MZFXTvj0w8VmW1rYs2cPtlVEUxQK2QyHDh6koqKCd/zVW4hFDF7Y20Pbnk7k53qID/Sy57cPs+wdG0hvrB09K6moSNFq3IF2HCMO8QVT2jgw1rKqcPoSgZIw68Z901l8HoQrcb2xP3hHJWp7Hv05k0zBRpGDyJLEwpIAuqawqCxMbUmIpt4MnckCGc+mtKhRWxI4bkXtIbquc/nll/PII4/wwt4mXnmCJbVJ540Y0XmznHQiwx/cx5Et2sNb5KOGSgeQtRU8PCS7gKYYo7bIA36fPNUYtYw1n5i2S9tAHoCqqEFIV3Gtcao82yb28/+Hl+mDkhU4yV5kM4sVKMVtfAeeHEZ13JnLPZrIUvAEgvt8xTqKJcuprD22KfmQ8WadGhoa+OhHP8qdP/8Ze7Y9QcYx6Usk2Lx5M1dddRVr1qxhR1sS8/AAu//tKZauK2K6AbqfbKPniU647XJiQXVkVjJYhhepwdVCuIo25ZzBsZZVhdOXCJSEWTfum44egpp1kE2P+XNDidqaKtGTKpIp2ACEDYWKiIEqy/RmitiOS0lIp3FBnILl4rgejQvjRAPjNBkdw9lnn83q1auJRCIQ0I+7pHY65o1Mp6Et8qbtEjIUZDwcD4L7fkNJoZX2+mtQ5dioLfLzzZFV4AMFi568iev6yd+V0eMsQzk27PwFhZ4mMo7CwKrrKXMMJLtAT8MNOHI5dGdm7jU00aXgCQT3rh7B1SPjH3eCWaeGhgY+9clP0PLcA5iWw+e/+AVWrlmLLPsNchMZk/b/20U4VQQ8LFklG9YJZRwO/XQnCxurR2Ylx8oZPOsqP8lb1DYTTkAESsK8pSoyigxtiTyFwUq55RGd2rg/o3Rkojb4+Uch3X/Jj9uJfRzRaHS4e7hrFY+7pDbfCvydakNb5DuTBarDChHVIWWrFIpFcB3szp2Ur7p41Bb5+eToKvDK3l6QVRaUBllVc5w+bq4DL94N/U0EgwHUs95GPFyF2hYFokQW1Z6a19ApXAo+0ayTLMvEw35AWFdXN1y3bOgiKb0vgWo6lJkZEnoEW5PIhFUC+xLHzkoemTN44Y1+TapC0i+xMI+XdIWZJwIlYW7a/6C/5b5m3biHhDQZ2/FoGyhQEdEpj+jDgRAcm6g9VRNdUptvBf5ONUnyO8CnCzYdqTya5OJ60GYsJ5RooVzZT3np6ybXGHaOOLoKvOtJDNgu6aJJUFcoWC7GWK8Dz4Vdv4K+/SCrKOvfhlJa5yeyn+rX0BSWgie663HICWedxjE0KxldWUYhkUH3bMrMDIc8cAwV1pYjSxwzKzl8gROuhKoGaH4S9j7gv9ccvRlDEAaJQEmYe6yCX7jRc6FiFXBsDpHn+QFSJKARMhQkQJGl4TfoZN6acKL2RIkltekTD2o0Loz7das8mbyj0K0uZq1mUGMUCNi9wMQrsE/VZNuTjLVT8+gq8AnJo83UCSoyyyqCeHijd/MBklNEsgsw0AL9Tf4OysY3Q2ndzDzRGS5FMZFdj9Pxf2RoVnLZOxvpfboFJAnds4lkLTJRlarXrcRyPELaEYFSoQDPPw+Njf7Xyy8GMwOdO/2ZvI3vguDUfifT0eZGmHtEoCRMuym/WfQf8oOkUDmEyvAyKSxJx5P83AQ96NI2UGAgZxENqFy0upJEpsi2wwOTTtSeDLGkNr3iQW24btWiQBF7YRmasZpA8ZBfgTp26gKlSbUnGWen5tENb/tsDduTUGWJyqiO5Xijd/MNFpOUnSI8rcCa10HJkpmb2TgFpSgmsutxsiTP8d8PihkYbAo8PCu5sRbv66+m+Kc2nFSW2i21LHzberQlcWJBjeZEnrrywbIL99wDTz8F4cElRUmC1a/zz9vfBDv+D9a94ySe9YjJtrkR5gcRKAnTbspvFn37/T/LV5DMW+zryNBnLMSVZJ5rTiJ3FIkENKIBlcVlIeJBjQWxwEknak+UWFKbfkN1qwKBIn2GQia+GqvjAFr3Llhx6bQ0M52IoQ/4iXaGH2un5pENbzO5Ihnbf11UxwwUWQaOaHg7WExSsbLYRsxvCPzS7+DSL87sE53h/KOJ7HqcFEXHjC1FS7dA9y4oqx3+1tCsZET1OJhbguNJ1L/sEipLwpSHDRJZk3TBpjmRY0lZCO+XvyArG9iH2kcKUw7N4G37sV84duddSGYWybP9Iqf65EoHTPZ1JMwPIlASpt2U3ixcF/r8atypyFJ2tiXpTxVRPRPJtcmaDr3JHCFD4aLVlcMzRlNJ1BZmnyL5H3xJqY5cu0bczMLAYSg7TouKaTTZ9iRj7dQcypvJFm160n7/trhqExw836iGt7kEUt9+JDuHpEdOXUPgeViKwgpVYQdKx0y4PnJW8ujCndGAyuG+HKm8za6DXVhb93PovHqcfo3WlzqpqCofLExpwPq3wdYfQaqNSOcTSK4Nj9yKtH5yM0yTfR0J84MIlIRpN6U3i1QbWHk81aDJKiVbNKmJG8ieTV6JYLkeVYOzUQM5k0UlQREUnSZKQhpJ06avYgvRihhypHq2hzQsnU6TyWSIRCKEA/qYx4R1hdKwxvbWASKqhCF7lKr28PdHbS7o7oZCEtku+NWiRS/AcY0EpWMHeGNW03/pJaJXXcUS02FXySL2hypJlG4g5Hp4RZU9f/WfLCvvIk0XjcU+4nf9DNa8Hh79NyTHxNEjSLk+pBd+hoFOEdEX7kwmssqEuWVw2a0QW0oiaw/me0BeiWDLOhJQEw+wsCQ4nO8x3STHRLayfu6CcMpEDBVdlclWnk2ybL1fR2saWI5L3nTGvVlHbiEfg2ma3HDDDbzzne/kySefHPc4SZKIGCoKMn3ZIjHVwgUKlkPHkVXg+5vg4CNQ0+hXiUbCC5XDeL0AhcmrroZ4nOiBl6Cri/a8wVb7ZfQqZdio9BQqeLqlkaYBg+aqJXjV1YCHZGZGWsaEKpDyAwQpnNQQZDODmuvyl/CEeU3MKAlzi+eAomKVLMfO+vke/VkbW/YDppqoQVhXcVzv2DopUzCyrVkhUfsKop1PUmzbMZzwKnarnBqlYY2uZJFEzqQ0PPbMzWRNR4JtbW0t5eXlbN68edxjcqZN3nJZURXBtjR6m/zdfKMa3toJePGX4Dp4yy4iZYeRnOL4vQCFkxOPw6OPkv2Hf6Twv/eS7q3EC0kUlKHXlISHxKHSi1j4zb8mGwgT8SS8QAmKOXiB1PQYXs068pxE8rXoB3daEYGSMLfUXwbLLkY2bdTDKQqWQ1/W335vOHlCxhj5HtPgyG3NcvUGUuX1DKhh3G7/TVPsVjk1SoI6Xcki+WyWXHY3ipXDXfQy1MGA2LYcvzAjYwSvxTSYWT9R+YhlmulIsDUMA8MwBouOHruj03E9mhM5PA8WlQZZEAlzeNdReTOFAXjh535tpJIlsPp1eId24qmB6ckbmuGt/3PZmK1QVBX7czdhh2owv9FByPXfRwqKTj6oYmQ9UoUQNrJ/wRUaqd6t5XshWAWBGNIY/94TIfrBnT5EoCTMPYpKOODXSdnTmcLzPGTPRXfzw4dMdzHJ0duaj+3fJHarnBq6KhMNqOQTadJ7fgOSQre0HHUwYLW60kiq/0E4Kni1Tbjv0zDQDK/4NNRfMnzOyebMjVUj6UTaB/JYtj8DuqAkiOSYx+bNHH7cD+QildD4Fhij0fNJOwVb/+es47RCURUZ9cBBoqpJwTKI2XlSapBiUEVyHapXlo2+4Bqq3m3lKV9Yj1RIchZb2cZZkx6W6Ad3+hCBkjB3mCNXg5IksSAe5IXWJKmMiebmkVyHQjZLoihPezHJad/WLJy00pBOOl+BHSynSkoTk9tJDO5uLKmMIA++Ro4JXiM1ECyfWjuKk1gy6c+aDOQsJAkWlwX9wqdjpc6tvNyvubX4XNACYBVPfpxjOYWtR+aa8VqhhHEoe/A+lgXr6M1WEnKKuIMBaiGkUvPm1ZRH9VEXXMPVuze8A2/rD4mSZRWHwDu5mSVh/hNJF8Lc4Nrw1Hfh6f8Cy585ypkOKyojLCyLUNDKSXghelv2UFsSoHFhfFqLSQpzRyyooigyhdI1uJ6H0b/Xv+qXJYK6MnzTFHkkUdtysCQNSw2Rl4MTTtQeix2qJl+xHsZZMhlqeJu0FBJZk9b+HABVUWNUCx3JKSIXBkaSeRUVVr565hK2jShEa+bd9v/p4OoR7FD1Mc9deuIJ6joOsbQKtty0Ga06Tky1qJX7aAxtI2D3UhbSx77gCsRwG96IB9TQg9Sx/ZQ8F2HuETNKwpygFJNIXsCvwKsFyRZtknmLiKGyrqGWjldewHNPP4EmyaNaQAinH0mSKAlp9JeuJt/9BNpAM5Jj+ssYRxnKLfPsItqJlueY2LLa8ZZMUnmLtv4cB1MRbFei+GI34aDBisoIldEjtpAPVt3W8r3w6064+HOwYNNUfi3CyTjrLOL//I80vvtaIimTg9EXcTyJFZ98E+7v/0hk3Qr6cxalYXtUkDvEjS3m6YEyKouHWfSXu1hevlzsTDwDTSlQqq+v5wMf+ADve9/7qKqqmq4xCWcgtdgPVEP5CjzPoyPpzyqVRXRChsrGc85n1dr1nF1fO7+CpDM4wXaixkrELQvr9GVKyBo1RN0u1EIvVvjYlibDidqWxsCJluemuBOp6Cns7EhRMG1CikPG8/Oc+rImpWGTVMH2ZzkHq26r+T7wHKREEzz933DFreJD9lSrrISPfYw40BgpjBSmXFVDZOP1NPXlyBYdmnpzLK8Mow3OFtquRPPzO7n/vt/yi4dfJGgP8FDXTtYe/i/edOXrEfsTzyxTWno7ePAgn//851m8eDFve9vbePDBB6drXMKZxPNQigP+38tX0p+zyJsusgzVg1fp4UiUyuqawV1H88RQgq0R8ZNMhWMNJuK6WnjU7yig+ctr+bIG8raLmu1CLiaPqUmjKfLwUtx4y3NHOtGy2ng8D9KuQa7oUBsL4CGRdlQMVeasBXFM26U5kcXzPMgnkBIH/TYYkgI1jeCYftVtRtd1sgfLXEymrpNwcoYKU8Y1v9eeLMvUlYcJ6jKO67GzLcn21iTPpyLctzfLP/zLbfz5iWfRg2H08jrKF69i27Zt/Me3/j/2HW4b8/UonJ6mFCh94QtfYMGCBViWxV133cXll19OfX09X/va1+ju7p6uMQqnOdnKILkWnhrAiS6kM+kXeKuOBaZt+/+sWbARtrzP/1MYU75iHYmGd0HtxlH3l4Y0CiUrMZNdSImDeIf+TO6Bm/Hat530Y3mK7gdlk9yJZCNTQKUkpGE5Lj1Ff/YqHlSJGCqlIX2kAGqmBy8/4Fd4VkN+OYNg6XDV7UTWZH93hkO92eHg7lBvlv3dGfZ3Z0hkzZN+fsLkKLLE0vIwpu2ypyPNI48dpvvpw/zpzj9y8C97CIdqcPUorqQSi8VYu3YtvS0v8dDvHyDU9jg8ciu0b5/tpyHMsCl9Cn3lK1/h8OHD/PrXv+b1r389sixz8OBBbrrpJhYvXszb3/52McsknJCa60ayCxCpoTtr4bgehiZTPk0FB2eEXfSXik5UvfsMTrA9niNnVXJSkIxWQV4KjZpVKQnp5ItF+jImz7Tb3PdSgZZD++l77A6SA4lTNlbXdUkMDNDb1UF3ewtt/TkcJDTJo2LwNaopMrbr4fYegIMPI9Wc5VfdluRjqm6XhXXqqyKsqIpQFtYpC+usqIpQP3grO/p1P9HX2jxx5L99wfJvszmjpsgSSB5d+3rZ8fUn2deapL2/H31Ap+P3BynkbExPwfM8JLvAIruJl9oGaO63/GbG23/iL7cKp60pJ3PLssyVV17JlVdeSXt7O9///vf5wQ9+wOHDh7nzzju56667WLZsGX/zN3/DtddeK3KZhNE6XiDUvRXZKeDu/yNZdRmUr6M2Hpi7uUhncs2aaXJkgc+h2kYHe7PD36+KGQQ0hf7eLmLFNEkpTCRgE6ipx8wkaGpqZuXq6IzvfNy9ezd33XUXjz38EBnL49/bDrKwbjnra1W2LCsdfo1ajovuFjBe+jXg4C27eNyq28OlKGwFBmdMNU0BdYz6FKfha20i//YzVdx1rHy4rOmQztukf/0SuB5ZxyHnSgTQ8YBCRwpjRTmW46G6RcIBlVZHIp9N+bOE+f6Zb2YszKppXddYsGABX/rSlzh48CD3338/V199NaqqilkmYWyFFDz/U2xJo6CWkMkXKdl7JzE5TzQwx7f+iyW1KRmaVRnvVhrSONyXxTZK8Iw4di5DIpUlc/h51FApSSIjOUEzZPfuPXzzm99k+/bthENBYuVVKMEYTXtf5LGHH6GtvWv42P6cSWlpHH3t66C83q+6rQZwjfjUZxNPs9faif7tj5lRmy7j5MPZjuvnib2UIJQsYqg6sqpRUGzwwM6ZeIDref6OXDlGQJWIaB60PQeBuGhmfJqbkfIAkiRx+eWXc/nll9Pe3s4111zDn//85+Fcprvuuov6+no++clPct1116EootLfmSjV343d30uvVEseA8kppyTVS9hJATUn/PlZZUTFctoUnKjAZ6Zok8iYZAYS/MvjRZqe7CHg5nlkb4pFGxewpbwHJRAja/qJudPNdT3u+dWv6O3tpaGhgc7eflRPIhaLUlEao+P5h3j02RdZ9SaLdNYZKYAaPAuq1vqVwqfLafZam83irmMVplQVGVWWiK4so/BUhkWGSkW8lP19STS5EjXkN+OWJQlP0Wi1S1m3sIy6mASS7Bf5nMRs0lBfyfGIvpJzz4z9azQ3N/P3f//3nHvuuTz66KOAH0Bt3LgRRVHYt28fH/zgBznvvPPo6emZqWEIc1Qyb7GrXyFNmIibRHFNYlY3BSXG/rRGMm/N9hCFWWQ7Lgf27eXOH/4X+1p6aZdrscO1hCsWsWvXbu684zsc2Ld32poiH621J8mevXtZvHgxRdvFQkGRPBaXhigL65SUlrKvtZfDB/aypu8PNFbII8uAc3XJWBizMGVY99slLXtnI5IkEfEsNi6oI6BF6KUHudxAdi3y2TS7du2iYtEKXvXqy8guvgiWXgi5BPTum/AYhpL5X+pK8+D2Azz87A72t3SIZP45bFoDJcdxuOeee3jd617HihUr+OpXv0pbWxtlZWV88pOf5KWXXuK5556jpaWFv/u7vyMcDrN161Zuuumm6RyGMMd5nsfhvixpR0OuXoslaQTdDE6gDG/jNaS84LQuq8y15FHhxGQJnnjofhK9faxavRpZM2hzS1HDJaxeWEK6bR9PPvwA47RrO2lDVbc70zaZXB7NCJAq2ACouJSENKqiBkvCDhEnyYquB1jstBA/dJ9ocTFPSZJEXXmYpRtr2fL1V6NVxVhdHubl515Fw6svwPYsMoluEokEmzdv5oYPf5iVdQuxYovx6i7wT3LoTxP+9x9eeizTOW/vP3Hh85+h3twz80uPwkmbljnrpqYm/uu//os77riDzs7O4Q+4l7/85Xzwgx/krW99K7o+8o9fXV3NzTffzJVXXsnLXvYy7r///ukYhjBPZE2HRMak2u1BzXbSIVdwMLyOBWe/i2hpNaWWM7zVejqWVWYzeVQ4OX2dbXS1HCReWYMsy2iuiSUHSAYXo+cPUht26D74In2dbcSWLZuWx0zmLZq60hxMRehwLdJWF7ubeygvjaPgouE3cJMyXUj9TZTaaSra/4RU/Up42f8TM0nzWDyo0bgwTkT1OJhbguNJXLPlbaSKFr/4n/9CNbN84RMfYuXKleBY9Bwa/MFlrwRVgyXnT/jff7hJs6SQDVdTDJYRWbwJeZoafAvTb0qfQnfddRff+973eOihh/A8D8/ziEaj/PVf/zUf/OAHaWxsPO7Pb9myhZqaGjo7O6cyDGGeGUqejOSaKdou/VoVRb2UULQEGNlqPV3LKkPVm8dzTHNVYdZlMhlUzyYajZAp2ICE4prk1RgJSojpCVKpVjKp6dmWncxb7GxLks4VCCkOK6rCVC1ext5dO1EDQRRc/3PQLuJ1vEBHb4pNCw2WxDwws6DM8c0HwgnFgxqNC6Ij1buXleJ6Eo+UluBRQrC8FlmWRzc8llWov+ykH/N47XKEuWNKgdLb3va24b9v3LiRD37wg1xzzTWEwxPvXn3kTJNwZhhKnpT7D1KwXNJqOXHVRmJkq7UqS9NWbHI2k0eFkxONRomGQ4RlG9tQcSUFD784o11SR7Grj6DsEC20AZOrsn20oaXgbNGmNhagX/LodIKce/Hl9Pd00Xn4ANmCRUyXSfZ20L7vIBVhhTesK0FetGWk6rbYHj4jTLOIZVqk02kCAePEPzAFQ9W7USBiqNieREQqkvYMejMm8ZBFRB1pcxIs2kTVI5rqdu6EYAmEKmd0nMKpNaVPokAgwHvf+16efPJJtm7dygc+8IFJBUngL9s5jnPiA4U54ch8n7FuE8n3CesKFVqe3EAXjiSRU0uIKCOvgf6cSXlUJyymos9YS5YsYc2aNfR0tlMR0QjaKUJOimXlQapKorTmDBas3kTtxldN+bGGloJLQzqe59FraRRcibrlK3n///sg69ZvIJcv0JtIkGg/xNmLI3z45RHqF1eDlR1VdVuYXrquc8FFryIcibBjx45ZGYMhOQQlf3PJ7s4UW5sHeD4V4YV0hGcPD7CjLelvPmnfBrvvhV2/Ais/K2MVZsaUZpTa29spKSmZpqEI88FQvo/reRzs8fN8lleGkQevqCaS7yNJEkvposmFXrkcTZVxcSlYDslsfmSrtcj5OGPJsszVV19NS0sLe/bswbaKaIpCIZuhrbOLskUruOAtf0NrymJ5xdRmpYeWgj08WgfyZBw/QK+OGURq1rJk6XISpkLM7uezN1zP8vZfYe7+HZJn44XKkY6ouj35BxdNk09k7bqNLF2+krPra2dtDCEswGN/V5aiaVIrOwRkj7Cu0JkskC7YNFavJB58yp9d3Pd7P7lbvIedFqYUKIkg6cwz3K3d9TBtf/aovjKCPLj96Oh8H9nMINtZKEoQHPkwkQcOEg9p9IeXI3c0kbA0MB1qS8IsKQvPeMVlYe5raGjgox/9KHf+/Gfs2fYEGcekL5Fgy9ln89rXvwGtYjG5okN7f5YaOTX+csgJKLJEpmjRly2i4qHgUaHbhIZmNIsp4qWlVMsaS1euRV65hkw6j2QXjqm6PSmnYdXtmRCORAlHoqesIfao6t36yGO6Lliuv/MxZeuE9CIBTSEU1OlI5mlO2TQ2vAFp24+ReveiZdpw9KjfOFeb2SVDYWbNSMFJ4fQ1vGPD9YZ3kAV1ZThQGsUxqX38S+iZVrjsc1B/qX+/55FLtKErMvVrGlmY3eknT9aVEA2HxEySMKyhoYFPffITtDz3AKbl8PkvfoGVa9YiyzKpgkVLWzvFF+6hycyyIxnERqX18AAVcYu68tEB99DW/yODKdNx6UwW8DxI5mzqSjVKAkWG4v1A706Mvb9lBQUKDAb6koSnGHiKMT1VtytWgj65lIXT3VBRRtf1KFj+snzedEZdkM1IUcbB6t3hjif84HXReYDfFLk/Z7KyKkJPMkvBlUjYKkMLrsNNkSuqiSy/GLb9iEj7Y7haxG+cu/ndp01l9TORCJSEGWWHqnECZYRqRhJus6ZDx5r3o+c6WFq1mIEjkidFkCQcTZZlgrqKIkuUlZUhy/4HZCygEQhFSSTTSIV+FuSS2KFKdClPZ1Lxl0MWxokHtVFb/x1PovXwAJqWQ1NlwrrKgtIgQV3Bc208x0R3shhNf0Tu3oEiw0K5j4PSDCRrn2ZVt6fLbJb0GKt6t4uE7UI0oCG7BkkgZatkBwPuUTt1q9bgDTT7M1NSESnb4zfOLVsuEv7nKREoCTNqrO2vvZkiSDLhqqVoYzUCFYQj6LrOxZvqeXxnEzt27uTSmgXAyAxRe9UraGz6H6IDT+LmIsg70sRXv4ND1NOcyLKkNMTO9tTw1n8ZSBdsunqLhHSF9YvibFpcSt5y6Nr1GErvn4ia3Rjpx7Cr16Ge9SYOvvjEhMY6PBNiOdiuX08ubzrInj8rItpTTMxslvRw9QiuHvEDWNtP4pbxUGUwbZewoRJTbVK2SnfaJBAI4LjeyE7dXD9SoBRXDSDh+rOFonHuvCYCJeGUKlgOqZwFkkRFxADs2R6SMEcdGXSctbyGpbVlLGjcQN70g46C7e9Wq11YR3hvJ6pnUjTz0PUiau528ms/xQvpIId6suRMh5qIQtJRSNoKYdejKmZgOy6SDJoioVtZogd/jGk348oyhmwiKwWspecAT05ozEMzIZJjUjm4THSoN4un+B+4orjpxMy1kh4qLuVhg96cSXVYoUy1KTgKjufRlS4gAwtKg/5OXa8ML1wBso6tRdCsPESqxM7IeUwESsIplUj0Uv7ij1EqVhJY+AbcKcZJs5bLIMy4UUFHOEBJOEBPUaanOwOAofnLHRFSqME4eSWG7NlIrkNZ4gVCmVZa5aUk8zZlYZ2+rEXa9t/yDFWmpjSI5zFSBT7ThdSzG1mSkGUFedEWP6k33z/hMQ/PhDgm6mArilhVBBT/76K46fwkSbCkLEjeydORyiN5EhW6Sb/j0dKXZ0FJYGSnbiAG69+B3fQMslOc+s5IYdaJQEk4ZWzHJdexl2gxRYmbmJats6I9yenrREFHwXZo78+Tk+K4oQoU9RCOHCGgSjiBCoI1y6nIajiux2Kln7xeAuSIeVmqg0UC6S6UVAvJ0gv83JJoLV7JYrzsAHagHM3KTXomYHgmxFZgMEDXNAXEEvO8FwtqNC7UaeryOOgoOJ5ESJfRdYXyiIFy5IaW2vVkF7x8ZGekosLhJ6Du/Mk/cDHtV3/XwyKfbZaIQEk4ZfqyJnryMJoiY1StnJZzivYkp68TBR0BTaYsotOZdInWv5lQ+y5UJ0++8mwGVv4VOcIsrVQpZDPU7b8Ho5jA638C1SmgP/QXiqVrKAQqCYVWoSrlYMTgwk9h//qTo2cCxIeTMOiYNierKkiaMJCzaUnkWFkVGe4oMLwzUlbh2R/4NbOCJVDVMPEHtE2479Mw0Ayv+DTUXzIzT0w4LhEoCaeE63r0pfKUpZv9+jRly6flvHMtl0E4dYa6vqcLNgedetSyiwm5OfR115HwIoQNldXVMTqau8hbDqXdW/HMPhxZRzHTGP37OFz/CuIlJSNV4I+eCYhUgGlSLBYxTZN0Ok1pTPTlOpMd3eYkEtLJWxmKlktrf56lFUeVetCCsPBsOPw4vPQAxBZObhkuUgPBcqidWqse4eSJQEk4JfrzFnK6FdWzMEKlEK058Q8Js8JxHCzLmu1hjLBNkEP+3wsFUEfa5BgSrKowaOszadEryHoQdzRqYjK1cYOA7FJRs4TD2TcQyjQDEZBVvPJ6XNvGrWygqjROsegv37qWiWmUgwEFT0cuFAB47Wtfy9NPP82LL77Ilk0bsdTI4HCKyI439phdDVwLUgkIlczkb2jWaJqGopzZVyqyLLGkLMT+7gzpgk1Pukj50av9Sy+ExEFId8Le+2D92497zqHcSxwHVdJA1bDlIJhi9+RsEIGSMOM8D3ozJoFkEyFdQSpfLkr7z0Ge59HZ2cnAwMBsD2U0z4PYuf7fm1vHfO3onsfSVRvxAEXOIOdy9Oagd/D7wUCI/nUfwLNNPEn2k/0lhaCu09vROnyc53m4NRcCkGppHa7rVV5ezhVXXIEkSTS1tI55zDGqXuUvt7R1gDrxhPD5pqSkhJqamjO6BlpAU6iNB2gfKNCVKhCU1NHFTcM6UsMb4NnbIXEI2rZCdeO45xvKvfTsItrg5gWrK42k+hW+Re7lqSUCJWHGDNW56bdUzJxFXbqJoKZA2YrZHpowhqEgqaqqilBoDlVI91zIDoYy4QqQjr2S9lwXO9MDgBqpRJLHOMbM4aQ6kfCQ9BBSoARJD07+PBM4BgDX8ccuySCffrMunueRy+Xo7u4GoLZ29nqxzQXlEYNM0aZ9oMAfuwawUxG8weKmfqX4OPEVl8C+P8DBhyA6/u9ruFWUpTEwWF2+pDKCPNgKReRenloiUBJmxJGVkNsLBsH2JGG1mlJNIVi2bLaHJxzFcZzhIKm8vHy2hzOa54Kl+jNLuu7vIDrmEBfb9O9XA4GxAxxdx3Zy4Lmo8Vok5dhNAMc7j+t54IHnuShF/xjFMJCGAjeJ4ebQZ4pg0A80u7u7qaqqOmOW4Ub1gzuij1ssoPFEXy+Hd3WiPNtFPJlgz28fZtk7NpDeWEvjgvXEy/b7s0rN49fmGm4VJSlkBnfTBXUFWSRkzgqxyClMu2TeYmdbkh1Pt9LxTBt9Tzfxwve2sr23ga01byNpj/6gkxwT2cr6bzrCrBjKSQqFQrM8kjFIsr/zTJLBzk/xXJI/uyNN/gPHcT0KtoNpu0iSfyrTdinYDgXbwXHHyFU6Awy9ZuZUXttMGuwH52phvx/cIM/zaOnPkTs0wO5/f4pMv0XalOl5so2nP/Y7mrZ30Nyfw1v9Wqh7Oay6HMkpIheTfgkAYc4SgZIwrTzP43BflqbtHTz7yT+Q7DdxizbZ5zp47uODbxaJLJ43+KEyzpuOMDvmzHLb0bSQv+ymzV4gp8gSAVXB0BQUWUaRZQxNIaD6N2WsxtBngNl6zViOS950yJsOBcu/DX2dNx0sxz3xSU5SvmIdiYZ3Qe3G4fuypl8pvv3/dmFkLcAjpYZwAA849NOdfnFTKQzLL4KuXYTbHyfS8YTfOLd9+4yNV5gasfQmTKuhN4tDP92JLUNRUpFwqZV6SEqlHPrpThY2VvuVkA3/5TdWE0pBGEVWgNlddpAlCST8JcAhEmJjwiyZzWKzo/rBDbIdF9v1SO9LYOQsNNfGklVyEZ3QgE1qX2KkcW4hBS/8DNnOAxJStls0zp3DRKAkTKsj3yyKg29elXKKZSsPkMmHeWZfaOTNYtBYbzog2pMIgjC+uVZsVlVkVFkiurKMQiJD3M7Rq8dwVIliRKNmZdkRjXMTSIUBJMdCdk1wHNE4dw4TnzICMHoae6zbRKexh94swivLMIN+HF6r9gGQzMeIrSwfebM4gUTWZH93hoO9WQKaQkBTONibZX93hv3dGRJZ8+SfsDBzCgX42tf8P0+BlpYW3n/dddQ1nE24ahlLly3jYx/7GH19faOOu+TSS9FLF6GXLkJWFKqrq3nrW9/K4cOHRx1nmib/ett/snHTJkKhEBUVFVxwwQXcfvvtZ04ezjygKTJBXRn3dqovosK6QllEZ9k7G5EkCRWPEisLsoQZVFn49rMoj+p+cdNgGV6gBAnAcyCxH2R5Uu1ypus9WzgxMaMkACPT2K7ncbDHn75eXhke3sUz0WnsoTeLqr9q4OA/9aB6LuVakqwUpTdZyrJPNI68WZzAXLtiFCboF7+Az34WFi2Ca66Z0Yc6ePAg559/PqtWreJH//0tli5Zwt7mHj7z2c9y//338+STT1JWVjZ8/HXvvYa/v+lTKJEqmlta+PjHP85f//Vf8+ijjwJ+kPT6t7yLF3bu4h++8lUuvPBCYrEYTz75JP/6r//Kpk2b2Lhh/Yw+p8nwPA/HcVBV8VY+24YrxW+sxfv6qyn+qQ05laV2bQXlr1tDcGmcRSWhYxrnqnnLz70Llg33MZyI6XrPFk5MzCgJgB+U1FdFqK+MsKg0yKLSIPWVEf++qghl4Yn9B5YkibqyEPHlpTTceC4V5S7hgEVoYYylf/82lm6oHemyfQJz7YpRmKA77/T/vOuuGX+oD3/4w+i6zu8eeIBXXnA+SxYv5LWvfS0PPvggbW1tfOELXxh1fCgYpKa6itraWs477zw+8pGPsHXr1uHvf+O223j08ad44Fc/58Mf+hAbN25k+fLlXHPNNTz11FOsXDnUo9Dzyxa4I1ftd9xxB0uWLCEUCnH11Vfz9a9/nZKSkuHvX3vttVx11VWjxvPxj3+ciy++ePhr13W59dZbWbZsGcFgkA0bNnDXEb/HRx55BEmSuP/++zn77LMxDIMf//jHyLLMs88+O+rc3/jGN6irq8N1xczCqRIPajQujLPuZYsoO2cJZZeu5bIvXkTDhlrqyiOkC/bIwYPtclJLXoXXcKXfa3D/Hyb8WEe+Zy+JOCwz0tTHvEm/ZwsnJj5pBGB0UDK0zHWyQYkHLK2IsHpdNasaAwSXVVL1uvM46+UraVwYJx4cf5ZImOeyWXjgAf/v993nfz1DEokEv/vd7/jQhz40XM9nSE1NDe9617v4+c9/PrLDcoyf/7//+z/OPffc4ft+8pOf8KqLX8Gm9cdWTdY0jXA4PGa5gqeeeorrrruOj3zkI2zfvp1LLrmEr371q5N+Trfeeis//OEP+c53vsOLL77IjTfeyF//9V/zpz/9adRxn/vc5/inf/ondu/ezRvf+EYuu+wybr/99lHH3H777Vx77bXI4xXEFGbEUOPcDbEM66MZXra0lEvWVBINqCTzFv1HpAx4ioEbKIWz3uxvCujcCT0vjTpfJl+kK5EmnR5dQmD4PVtxWPb0zdQ/+lGCPdvEheQMEPO1wrTrzZhEDJWlS0tRth1GU1LodcsILYjP3e3nwsl56SW46qqRgMg0YbBvGsUi1Nf7RSIBwmG45x5YtWpaHnrfvn14nkdDw9jd2BsaGujv76enp4eqqioAvvP9H/KDH/0Uz4NcLseqVav43e9+N+qcrzz/nSd+cC0EamC4Svhtt93GFVdcwWc+8xkAVq1axeOPP84DQ0HjBBSLRW655RYefPBBzj//fACWL1/OY489xne/+10uuuii4WP/4R/+gVe/+tXDX19//fX87d/+Lf/2b/+GYRhs3bqVHTt28Ktf/WrCjy9Mn6Mb58qaRlXUoCtVpD2ZJ2QojLpcjC2Exef6RSj3/wHKV4CsYJomn//e/bR0D3BzzWW85rWvG/Px7FA1TqCMUI1onDsTRMgpTKucaZM3HSQJKpK7iCdeIJJrJrj/PqSO52d7eMJ0q66GeByam/1bZ+fo73d2jnwvHoea6S8BMd6M0RBdH1mCeOdbr+aZP/+O7du28dhjj1FfX89rXvMakqkUruvheR6e51cAcD1/x6Xren5V7iPJCijacGuS3bt3j5qZAoaDnYnav38/uVyOV7/61UQikeHbD3/4Qw4cODDq2C1btoz6+qqrrkJRFO6++27AXwa85JJLWLp06aTGIMycyqhByFBwXWhJ5HFdl6wjk7QUMkUbb+mFUNPoN8w9ouVNdWmE9Stq2bRp47jn9hTdr0VnRE7BMznziEBJmFa9aX9auVQpoOy6E08NYQfKkayCXyekkJrlEQrTKh6HRx+Fv/u7warXR72lKIp//9//vX9cbPq2PtfX1yNJErt37x7z+7t376aysnJUnlA8FqV++TLq6+u54IIL+P73v8++ffv42c9+RsF2WLlyJS/tPzAjVbdlWT4mqDtyF10m41em/+1vf8v27duHb7t27RqVpwT4S4BH0HWd97znPdx+++2YpslPfvIT3v/+9095zII/05dOH7v0NVmSJLG4NIQsQ0+6yKP7EzyfivBCOsKzhwfY0ZElufQKv7DqEXRNJRI0iEaj45xZmGkiUBKmjWm7pAr+G3+FkkUqJLEDJTiBUojXjtQJEU4vqgpf/jJ8/eujkpsBvz7Mv/0b3Hyzf9w0Ki8v57JXv5r//M//JJvLjZoJam/v4H//939573vfe9xzDPUmKxYKBFSFa665hj8+8igv7HjxmKrblmWRHSfnqqGhgaeeemrUfU8+ObqXV2VlJR0dHaPu2759+/Df165di2EYNDc3U19fP+q2ePHiE/4+rr/+eh588EH+8z//E9u2efOb33zCnzkTTKWCt67rXH755USjUXbs2DGpxx2rNZOuykQDKge6M7zUnUGWPMo0i7Cu0JkssLMtSTI/GDwnW6fc2kSUEJgeIkdJmDZ92SKeB2FDwYhW4AZKUMwMjh6BVAdEqiZVJ0SYZ7Zv92eQHGfkPkXx758h37jtm1z0igu54orX8g83fZyldUvY8dg2Pn/T56hfuZLPf/FLo47P5fN0dnWj5iS6e3r4yle+QiAQ4PLLL0eWJW78+Me47957ePUb385XvvIVLnzFK4hGozz77LP88z//M9///vfZuHHjMeP46Ec/ygUXXMC//uu/8qY3vYnf/e53x+QnXXrppfzLv/wLP/zhDzn//PP58Y9/zM6dO9m0aRMA0WiUT33qU9x44424rsuFF15IMpnkL3/5C7FY7IRBX0NDA+eddx6f/exnef/7339MgvuZaqoVvM8++2xWr15NJDKJZa3B1kzhjif81kyRVwL+MvFAzvJTEyIGaVslqjgENIVQUKcjmac5kaWRQ0j7foccqsbfHnNyuZ2ihMD0EDNKwrRwXW+4AGRF1PCvhGILQVJQrBxeqBw2XiOqzp6uTBN++Us/SNI0+Ju/8f90HP9+c2aKg65ZvYqnnnqaFSuWc837P0j9+nO56o1vYPWq1fzlsb8Qj41ervj+//yEJWs2s2DhQi655BJ6e3u57777WL16NQCGYfCHe37GZz72Qb77ve9x3nnncc455/DNb36Tj370ozQ2HrsbDuC8887jv/7rv7jtttvYsGEDv//97/niF7846pjLL7+cL33pS3zmM5/hnHPOIZ1O8573vGfUMV/5ylf40pe+xK233kpDQwNXXHEFv/3tb1m2bNmEfh/XXXcdpmmKZbcjDG+jH+d2om300WiU2traSS99Ha8f3IqqCKosYXsSCWskrbs0pPv94EILQFKQBw6zIpihTCue1OzSdJV9OdOJGSVhWvTnTFzXn1qOBTRo2YVkZSnGlmCFF2BcfBNEKk58ImF+evxxyGRg9Wq/jtK6dfCRj8Bf/ZW/M+7xx+GIekHTRZYkli9fxv/8z//4HyTFDH//z7fxb9/8/9i5cwfnnXfe8LEPP/QQdspf+lJjtUjjbJs3DJ3P3fhhPvelfwBl4m+R73//+0cFKHfccccxx3z5y1/my1/+8rjnkCSJj33sY3zsYx8b8/sXX3zxcZPX29raWLduHeecc86Ex3260xQZbRbaBB6vH1xQU6iOGfQDGUchXbCJR3Q0RfZbPBmlsPxieOYHvLaygwFbQ33sa7DlWliwccJjGHrurusNz6YFdWW4DZQwMWJGSThpsplBzXVBMU3f0GxSRPcTRXr2AGCHanCN+DF93ITTzFlnwTe+AVu3+kES+H9u2+bff9ZZMz8GLQThCr78lX/km9/8Jk8++eTkiy2OUSNpPshkMuzcuZNvfetb3HDDDbM9HGEcQy2eTNslqCmUqH4Byt6Mien4OUPDLZ4qV8FAC4biEJAdpGyf2BAzS8SMkjAh6XSaTCZDJBLxp6Adk9rHv4SeaSX1ys9SjG1Blv2pY5ItUMzgqQEcIz7bQxdOhcpKGGsWJBQa+/6ZICuAf9X8vve97+TPc1SNpPngIx/5CD/96U+56qqrxLLbHDbU4qkzWaA67AdKBVfB8Ty6UkUUCRaUBv0WTwP9SME4OVtFkz08SRaNc2fJ/HknEGaNaZrccMMNvPOd7xy1k8cOVZOvWE9fxM/vKAvr/pRuz17/gLIV8+rDRjgDeJ6/M887zkzTUTWSTta1117LwMDAlM4xUXfccQfFYpGf//znwzv5hLlnqB9c2FDpSBUoehJlmonteLQm8li2N9LiKViGF67E8iRkPKTe3aAFxYaYWSA+xYQJqa2t5eyzz2bz5s3D93mKTk6OkHENAMrDxqhlNyrXzMZQBWGU4cKRSNhaFE+S8az8+MUkBWEGDfWDq40FyDkKKVslbMhURnXKowbD6UOBGE7j2+gxDVwkvPLVsOV9YjZpFoilN2FCDMPAMEaKnnmeR9aRaSsYqKbDgoiKrsow4C+7oepQWjfLoxYEcFxvuF6MpAVxVANHkvFsv4yBpsjIikhuFU6doX5w0VgG25WoWVlBypLoz1q09udZWaWgKjJe9Toe6F1ASLH55Of+GbVswUk/pmxmkO0sFCUIimBrMkSgJExaMm/R1JVmXypCaz5ArCONomqUhHXijulXlo3WgCxeXsLsU2QJRRpajhpjWUrESGc8y3GxHX+GsWD5AXTedIZ3h6mKNO1NZo/uBxcJ6YNFMV1a+/MsrfCrrxddhaKrjGyIsQpgZo6p4H1cR+SUctnnoP7SaX0upzvxSXYGGHoTGM9k3gRSeYu9PXnSuQKOKxFTbUrCGgN5i51tSRoXLiH+sg+AY4PrnPiEgjDDZEkSwZBwXFMtSjkdZFlicVmI/d0Z0gWb3kyRyNGv20w37LgTJAW2vN+fuZ+g4zXOnc7PiNORCJTOANNVndXzoDmRJ1u0qYkZ9LgysgRVEYN4JDhSVXZBHElRjxsozcYVnCAIwljKwrpf/20c6ilamg1oCrXxAO0DBTqTBRaEXCxJx5NkMkUbIxpF8jwo9MPBh2HV5RM+t6foeIo+ZuNcUcH7+ESgdAYYehNwXQ/T9nM16isjo4KSibCR6csWKQ0HyRQKOEgoeEQC/hVYpZQkMWCTrYgQMY7/0poLV3CCIAgwe0Upx1IeMcgUbdoHCjzckqDJimPZNvc9+SLrG9eyYulriO79BbRthYqVEF045cecrs+I05W4ZD8DaIpMUFcI6goBzb8NfR3UlQnP3LhI2C5oqkR/zm/cGFMdJCTwPKoO/YqFu/8bZ6D9hOeaalsBQQB/C74kScfcrrjiiuFjli5dyje+8Y3ZG6QgTFLUUHn82e388Hv/H89ve57d257jJ9/9d277+r/w6xd6yVYOLp/tuQ+sqRdGna7PiNOVmFESJkzGQ5UhkTExHQ8Zj+hgZVkt2w5mBllWUWJVJzzXXLqCE07eC10v8N9b//u4bTWGSJLE9ZuvZ331+mkdwxVXXMHtt98+6j7DMKb1MWaCaZrourggOJNJjonkWv5OYc1/zXqex5+f2c4f7vofBvq60QJBArJBWXk5bft38dPb24n8v7/hjaFSpHw/HPjjLD+L058IlIQJU3EpC+m80JEhpkvEVWe45ofR/xIFy0Fd0EA4MPc/pITp8ejhR/mPp/8DRVKQj1Nc1PVcHM9hTcWaaQ+UDMOgpqbmpH/ecRw+/elP84Mf/ABFUbjuuuvo6uoimUxyzz33AP6s1Mc//nE+/vGPD//cxo0bueqqq7j55psBGBgY4FOf+hS/+tWvKBaLbNmyhX//939nwwb/6v/mm2/mnnvu4SMf+Qj/+I//yOHDh7njjju48cYbaW9vHxXcXXXVVUSjUX70ox+d9PMS5jhFJ7PoIsIdT0D3ixB5JQDpgsVvfv0rrEyS+pWrSfX34gKhUITlq0rZt2cPv/nNfVz62euJ7v45UvculHwfTrB8dp/PaezMnk8TJkWSoCSkockyiZyJLju4HhRMG7NjF5oiU7Z0g19VVjgjvGfDeygxSnA8B8u1xr05nkOJUcJ7Nrxntod8jK9//evccccd/OAHP+Cxxx4jkUhw9913T/o8b33rW+nu7ub+++/nueeeY/PmzbzqVa8ikUgMH7N//35+8Ytf8Mtf/pLt27fz1re+Fcdx+PWvfz18THd3N7/97W9FK5IzQL5iHYmGd0HtxuH7mpqaOHxwP7ULFxLSFWTPxQMyRRuQqF24kKYD+ziUsGDJ+XgVq3HVAHIx6TeGFqadCJSESSnYLiuqItRXhCm4CglLw062UqIWqSqNEa1dNdtDFE6hqBHlsxd+1s9TOw4Jic9d+Dki+rE7bqbqN7/5DZFIZNTtlltumfDPf+Mb3+Cmm27izW9+Mw0NDXznO98hHp9cj8LHHnuMp59+mjvvvJMtW7awcuVK/vVf/5WSkhLuuuuu4eNM0+SHP/whmzZtYv369QSDQa655ppRS4c//vGPWbJkCRdffPGkxiDMP64ewQ5Vj2oans9lsYoF9EAIJAnVM5HwsD3ImTa6EcQyi+RzWVj2SihdSrjzaSIdT8Ajt0L79tl7QqcpsfQmTFjRUyjaLrGgxtkLy2jd61eVXWR0EowaSNWrQREvqTPNh8/5MP/82D8zUBwY95i4EefDL/vwjDz+JZdcwre//e1R95WVlU3oZ5PJJB0dHZx77rnD96mqypYtWyaUdzXk+eefJ5PJUF4+evkjn89z4MCB4a/r6uqorKwcdcwHPvABzjnnHNra2li4cCF33HHHcJK6cOapKishGgkxkEpTURJBwkN1TSSgaLsk0imi4SBVZSX+DNILP0e28zhaBCndAdt/AmXLRauTaTTvZpTuvPNOLr74YkpLSwmHw2zYsIGvfe1rWJY1qfPccccdY+6WOfL2wAMPzNCzmJ/ynl9npDJioCoyYcUlrtoEUgf8GYWqhlkeoTAbTjSrNJOzSQDhcJj6+vpRt4kGShMly/IxgdOR7zmZTIba2lq2b98+6rZ3714+/elPjxrr0TZt2sSGDRv44Q9/yHPPPceLL77ItddeO63jF6bOclzypjNYPdu/DX2dN53hNjlTVVdXx6Z1Z9Hf3UE6b+EhI+OiKRLZgk13ZzvrzlpLXV0d5BNIhQEcLYxqpiDbC5kuyPdPy1gE37y6/P/4xz/ObbfdhqqqXHrppUQiER566CE++9nPcu+99/L73/+eYDA4qXOuWLGCCy+8cMzvLVw49foUp4uip2AjI0sS5REDHNP/hiTB+nfCQBOULpvVMQqz53izSjM5mzRV8fj/396dhzdVpo0f/56sTZPuLW2BFgqUUpSlLCqCwgAKbgzq64KKbL46P0FgRBlcUFRex30fndGRxW1wdxxxYRlwQUQRUFBkbSmUltItbdLsOb8/QiOlDbTQNm25P9eVS3rOycmdmKZ3nvM89x1DamoqGzZs4PzzA5NpvV5vcI5RjaSkJAoLC4M/V1ZWkpubG/x5wIABFBUVodPp6Nq1a6PjuOmmm3jmmWcoKChg9OjRpKWlnfyTEs2ipWq/aTQaJlz9Pxw4sJ99+/fi9PjQ6vQ47TYqi4vplJrM0FEX4VdBY4pHjYhF67GDzwPOikC17lMYTQrVE+50rt7dZhKljz76iGeffRaLxcKXX34Z/BArKSlh5MiRfPPNN8yfP58nnniiUecdNmwYS5YsaYaI25ea0aQEsx6tRqlddNsUC9Fn13u/+pa/ivanZlTp7tV3o/L7h2lzjyYBuFwuioqKam3T6XQkJjasF9asWbN45JFHyMzMpFevXjz11FNUVFTUOmbkyJEsWbKEyy67jNjYWO677z602t/rW4wePZohQ4Ywfvx4HnvsMXr27MnBgwdZvnw5l19+OYMGDTpuDNdddx133HEHr7zyCq+99lrDnrhoUS1ZvTs7O5s7b/8z7yz7F8t++xG3X0XvsTNm+BAGnH8hnTMyyS+rJiMxCvpeiz/vB7SoqMZolPjuUPgzdBve+Ac+Tk+407l6d5tJlGomZ86bN6/WN73ExERefPFFzjvvPF544QXmz5/f6ImY4vgqHR68aFBQSWhMIcgQy19F+1TfqFJLjCZ9/vnnpKam1tqWlZXFb7/91qD7z5kzh8LCQiZNmoRGo2Hq1KlcfvnlWK3W4DF33XUXubm5XHrppcTExPDQQw/VGlFSFIVPP/2Ue+65hylTpnD48GFSUlI4//zzSU5OPmEMMTExXHnllSxfvpzx48c37ImLFtXStd+ys7O58445HNj0BW6Pj7vvv5fMXr1x+1R2F9uwu3wcqnTRIbUv9o7nonidJJx1NeR+CfnrA/OUYhs/MhmqJ9zpXL27TYyTFRQU8MMPPwCBb17HGjZsGGlpabhcLj799NOWDq/dK7YFhptNihfdUUOrGncVEWW/weEdIe9b3/JX0T4dO1epJUaTlixZgqqqdW5HJ0l5eXm16h8dS6fT8cwzz2C1WikvL+fJJ5+sM5E6OjqaZcuWYbVayc/PZ9KkSWzZsiVYQwkgKiqK5557joKCAtxuN/n5+bzxxhvBy2gLFixgy5YtIeMoKCjg+uuvbxPFMkXL0Gg0xJgjSIo106VLFzQaDRF6LZ3jAlNMDle5sFa7sWGiXBOPLbEfasqZgcac2/8DHmejH1PVGvDrzcGecFVVVRQWFuKstp+21bvbxDPbvHkzEFjJkpFR/zyYmqHtmmMbavfu3dx7773cfPPN3H777SxatIiSkpJTC7gdsTo8OD1+FCCC2hPmddVF6OyFUPRzyPvXt/xVtF/TB08nxhgY0W3Nc5Nak/Lycj788EPWrl3L9OnyeokTi400kGAxYHN5+XJXKT9ao/i5ysLGfRVsizwLhy4KnFbYvfKUHsftdnPbbbcxYcIEvvvuuyaKvu1pE5feaoa409PTQx5T863t6OHwhli3bh3r1q2rtS0iIoIFCxbwl7/8pZGRtg9HT+Y7XBX4Zh2heIJVuAEo/InI4s1ofK5AF+uUM6Fj/7DEK1qPmlGlu1bf1eyjSe1FTk4O5eXlPProo2RlZYU7HNFGmA1aDpRXU1blwuDVkh7hxGzQUmhTcUafz5nWLzGl9Dnlx0lNTSUhIaHWlJfTTZtIlKqqAtVG61taW8NiCXwgV1ZWNuicKSkp3HPPPYwbN45u3bphNBrZsWMHzz//PK+//jrz5s3D5/Nx9913H/c8LpcLl8sV/Lmhj99aqV4X0V8/gMZWSOGwWVTHD0KjgOno0SRnJWx8FY3XiV9nROfzSO0OETTz7JlYDBam5rTdytItucAjLy+vxR5LtA+qqrKvrBqLUYteiaBCgXKvnkS9lkiTgUJrB3ZnXMeZsfEnKAV7YkajEaPRSFTU6XtVoE0kSs1h7NixtTqMQ+Dy3dKlS+nXrx9z5szhwQcfZNq0acedjPnXv/6VBx54oLnDDakpl2xaHR7yDlWR6+uO15iBs6IjZp+Nnh0ia48mOcpQKvJRtQZUvQViOgXqdzjKJVESROojmXHWjHCHIUS7oKg+UP2BlcPmQLJid/sos7lJtETg82qoAGw+LRUOD/FRBuIiDZTavdjdPixGXeDLraoGyrmIRmsTc5RqMlm73R7yGJvNBgQmXZ6qWbNmkZiYiMvlYsWKFcc99q677sJqtQZv+/fvP+XHb4wyu5vdxTZ2Hqri821FfL6tiJ2HqthdbGN3sY0yu7tB57E6PGwrsFJY6cSkVzDp9WgMEZTY3ByyOnGpRy330JlQvS4Unxuf1giVhWCKC9yEEEI0Da0Bd3RXVI0ein8Nbvb6/Hj9KgadBpNeS5wuMOJfYnNjd3vRazV4/Spenx+KtsKPi9DbC0M9ijiBNjGiVFPA7XhJSM2+kyn2diytVktmZiYlJSUcOHDguMfWDEuGS1Ms2VRVlX2lduwuL6nREZRrVMo8Oow6DalxJuwOFwVlNiI9leTl5dHDUIKS3Bt/ZREanwc1MgGl/3UymiSEOK3VjPD7/SpOT6DYnMPtq/V53NjVYZ7IDngj4iD19+X6Oq0GnUbB7fVjUCBG78NzZNyjqNJJotmATqMEVin7fSg+N4bKXFRFCbQ9aaaadu21KGWbSJRycnIAKC0tJTc3t96Vbxs3bgRosglnpaWlAK3+umxNbQ+/Xw1WizUZtMFfzIaoGcaNizQAPhw+LS6/QoQCZQW5rFj+Meu/XI3O68Bmq6J3lw78sW8yqSnn4jNYMI64CywNK+4nhBDtVXNU71a1BlStodbKYbNBS7zFQJHVSbI58DiJOg/Veg0uP+wqtpGTFoPZoIXUfqg7V2CoysdQmQ9rFsLAyc2y+Ka9FqVsE4lS586dGTx4MD/88ANvvfUW99xzT63933zzDfv378doNHLxxRef8uNt2rSJnTt3AnDWWWed8vlau6OHcVWflzJP4G1RfjCPz999jdKSwxhMZmIMZhISEti8cz/5h21c30dPz/QYWfovhBC0XPVuRVHokmCmyumlsNKB4lfQKyqxJh07SlzoNRoU5cjIjasKKg+g+H2oGh1K8fZmW3zTXotStpkxsJrVZ4888gibNm0Kbi8tLeXWW28FYMaMGbWqcn/44Yf06tWLUaNG1TpXdXU1f/vb34Kr6Y721VdfceWVVwKBQpanQ6J09DBuRbUHj6qg+H1s+noFFWWldOneE6MxAq1GQ3R0NL1796aktITPv9+J3193mLWlmkcKIURrotdqahVhPPbWlJedYkx6zuwUQ2p0BNU+LWUePU6vSt/OMWSmWNAoCgUVjsDiG7cdT2QyKFrwe6F83yk1ztW4beiqDwWSsKMc/fzbU1HKNjGiBDB+/HhmzpzJc889xznnnMOoUaMwm82sXr2aiooKhg4dykMPPVTrPlarlR07duB01q5O6na7mTFjBnPmzCEnJ4f09HS8Xi87d+5k27ZtAPTp04d33nmnxZ5fONUM4+aX2XG7A5MCfVUl7Nuzm+SOnXB5VXT4MOKGygKU6I507tSZXVt3c+CwlWPXBLZU80ghhDidxZj0nNkxiqhoG16/QkqXWKLMkVS5vOwrqabc7iEiwkJ8RCyK34NPZ0LvqgqsflNOMmk5Tj+49qpNpXfPPvssb7/9NkOGDOHbb7/l008/pXPnzjzyyCP897//xWQyNeg8kZGRzJ8/n5EjR3Lo0CE+++wzPv74Yw4dOsTo0aP5xz/+wcaNG+v0j2qvaoZxnW4/h6tcaFHxux1UO5z4NEb0OgWD4iNGsaGU50HpbszmSJxuH3aHq8754s0GenSwhLzFN6ZfnBDHMXnyZBRFQVEU9Ho9GRkZzJ07t86Xo5pjFEVBp9ORnp7O7bffXqsG2pIlS1AUhezs7DqP8+6776IoSr2LRRwOB/Hx8cGVssfq2rUrzzzzzCk/VyHqoygKZq2fGH2gFICiKERH6EmNDXwZLXTqsWX9D25tJB50eCOTUPteC9EdT/oxvZHJOBL7wjH94NqrNjOiVOPqq6/m6quvbtCxkydPZvLkyXW2GwwGHnzwwSaOrG3z+1U6xUWixYtX66dKZ0arN6D1u+hgiSVf8WHGDiRAVAp2ezURBi1mU93VEy3dPFKc3saOHcvixYvxeDz8+OOPTJo0CUVRePTRR2sdt3jxYsaOHYvH4+Gnn35iypQpmM3mWiPRZrOZ4uJi1q9fz5AhQ4LbX3311ZCdAd5//33OOOMMVFXlo48+4pprrmmeJ9rM3G43BoN8iWkvEi1GnB4f+8sc/Neais80FpPqwJxzLXFJHeni8hNjOrkP6t8nmJ8elffb1IiSaB5en59CqxOLUceQjHgGxlQxsouOc3LOxFFejEGnwYIdDSroI1GN0RwoOEBm5yQ6J8Wc+AGEaEZGo5GUlBTS0tIYP348o0ePZuXKuj2uYmNjg8ddeuml/PGPf6w13xECDXKvu+46Fi1aFNx24MAB1q5dW29DbggkUTfccAM33HADr7766ik/n+LiYi677DJMJhMZGRm8+eabtUal8vLyUBSlVoPdiooKFEVh7dq1wW3btm3joosuwmKxkJyczMSJE2v1sRwxYgQzZsxg9uzZJCYmMmbMGKZOncqll15aKx6Px0OHDh2a5LmJxnO7Xdhttnrn1J5IlFFHQXk1h6pcWP2R6CPMRFjiKbI62VZgxWp3gq24GaJuXyRREhRanfj8KhF6DUlRRsxaP3FGlWuuvJykpCS2//orGmc5Pr9KBVH8un07iYmJjD07q1FlCIRobtu2bePbb7894cjIzp07+e9//8vZZ59dZ9/UqVN55513qK6uBgKX5MaOHVtvhf49e/awfv364Ej3119/zb59+07pOUyePJn9+/ezZs0a3nvvPV588UWKixv3x6yiooKRI0eSk5PDxo0b+fzzzzl06FCd0filS5diMBhYt24df//737npppv4/PPPKSz8vTjhJ598QnV1dZsdKWttGrPYxWAwMHT4KMwWC1u3bm3U49S0OTEbtaTGRKDRQKnHiFGvITXGhNNupeq7xaib3ww00BUhtblLb6Jp2VxeKqoDE7g7xZlQFF9wX3Z2L2bOnMl7r/2Dr/d8R3mVipqmMGDAAMZdcjGJuR+EK2zRSuWvy+erhV9RvLWYDn06cP6955M+NHQz66bwySefYLFY8Hq9uFwuNBoNL7zwQp3jJkyYgFarDR536aWXctddd9U5Licnh27duvHee+8xceJElixZwlNPPcXevXvrHLto0SIuuugi4uICVenHjBnD4sWLWbBgwUk9l507d/LZZ5/x/fffM3jwYCAwYlXfvKnjeeGFF8jJyeHhhx+uFWtaWho7d+6kZ8+eAGRmZvLYY4/Vum9WVhavv/46c+fOBQKXLK+66qpgP01xahq72KV3n/507ZbJwB6NmzN7dJsTxaelAnD6FYqr3KTGGYm2RGHf78GjVmP49WPof/0pP7f2ShKlVqy5q5z6/SoF5Q4A4i0GIg06/B5frWOys7P5yxUDyHRsYq8rnv/58wJ69OgBPg+Hc0/6oUU7lL8un6UjlqKqKqpPxVZkI3dVLpPWTmrWZOkPf/gDL730Ena7naeffhqdThcs8XG0p59+mtGjR+Pz+di9eze33347EydOZNmyZXWOnTp1KosXLyY9PR273c7FF19cJ/ny+XwsXbqUZ599Nrjthhtu4I477uC+++5Do2n87+b27dvR6XQMHDgwuK1Xr17ExsY26jw//fQTa9asqTe52bNnTzBROvpxatx00028/PLLzJ07N7jY5b///W/jnogIqbG1lsyWKMyWqEYXPz66Ph6KhiSDm0NuA1VOL3qbk7hII4WdxpBe/m+wHoC8ryBtyIlPfBqSRKkVa+4qp4dtLtxePzqtQkqo8/j9KIZIOsRGkkcvunbtikajwe+r/3Bx+vpq4VfBJAkI/Fcb2H7DZzc02+OazeZA8k5g1KRfv368+uqrTJs2rdZxKSkpweOysrKoqqpiwoQJLFy4MLi9xvXXX8/cuXNZsGABEydORKer+1H5xRdfUFBQUOeSlM/nY/Xq1VxwwQVN+TSDahIwVf39S5TH46l1jM1m47LLLqszoR2otZrXbDbX2X/jjTcyb9481q9fz7fffktGRgbnnXdeU4V/2muOxS6Kz43i9wQa5x5pT3JsmxOT1k+i3oMLqKj24vWr6E1xqEkXwa5PIP87sKQ0bWDthMxRasWCy+yTLHSOM9E5zkSPpJNfZq+qKk6bFUf5QcrKSimuDCyh7hhjQhtqrpFGg5o9jm8ZhBOpfSRCK95aHEySaqg+leKtLTdZVKPRcPfdd3PvvfficDiOe6xWG/hrVd9x8fHxjBs3ji+//JKpU6fWe/9XX32Va6+9li1bttS6XXvttSc98blXr154vV5+/PHH4LYdO3ZQUVER/DkpKQmg1jyioyd2Q6CV0y+//ELXrl3p0aNHrVt9ydHREhISGD9+PIsXL2bJkiVMmTLlpJ6LaCFaA7bOw/HrzVD8S3BzTX288urfG6NbdD4SzIHRrNzD1Rh1Wkydzgy0M1FV+G05GrcNjctap5hkUzl6jlZ9t9ZYkFhGlFqxpujjVsPq8JB3qIKiNS9CdTm7s67A1PEMenSwEBMZehi4hlfeKuIEOvTpgK3IVitZUrQKHfp0aNE4rrrqKu68807+9re/cccddwS3V1RUUFRUhN/vZ9euXTz44IP07Nkz5PyfJUuW8OKLL5KQkFBn3+HDh/nPf/7Dxx9/zJlnnllr34033sjll19OWVkZ8fHxABQUFNRJZrp06RKc21QjKyuLsWPHcsstt/DSSy+h0+mYPXt2rRpxJpOJc845h0ceeYSMjAyKi4u59957a51n+vTpvPLKK0yYMIG5c+cSHx/P7t27WbZsGf/85z+DSWIoN910E5deeik+n49JkyYd91gRfo7EPrhiu5GU+vuIUKg2J2a9lkKbm0iDFkUTmMtk6TEarAdQCn8iKv9LVF0krP0rSt9rmzzWttgPTkaUTgNWh4dtBVYKK50YI6MhJg1i0iixuSm1u7E6PPXfsSIf7KUtG6xos86/9/xAUccjcywUbaDA4/D5w1s0Dp1Ox4wZM3jsscew23+fJDtlyhRSU1Pp3LkzEyZM4IwzzuCzzz6r97IaBBKS+pIkgNdeew2z2VynPRLAqFGjMJlMvPHGG8FtTzzxBDk5ObVuy5cvr/fcixcvpmPHjgwfPpwrrriCm2++mQ4daiebixYtwuv1MnDgQGbPns3ChQtr7e/YsSPr1q3D5/Nx4YUX0qdPH2bPnk1sbGyD5k6NHj2a1NRUxowZQ8eOJ1+YULQMv8GCNzK5Tt/N+tqc2N0+zugYRf/0GMwGHftK7Tj9GugxGrU8DzR6fPpIlOpSlJ+XYaRuEdWGCNXmpKmvlLQEGSZo51RVZV+pHbvLS2p0BCVaLeUeM+aISDrHmHF6fOSX2TmzYwyKohx9R9i1AlyVkDUufE9AtBnpQ9OZtHZSrVVvw+cPJ+3ctGZ7zCVLltS7fd68ecybNy/489HzeUIJVaC2xuzZs5k9ezYAc+bMYc6cOfUeZzAYKC//vY9WXl7eCR/7aCkpKXzyySe1ts2fP7/Wz9nZ2Xz77be1th37HDMzM/ngg9ArU4+uuXQsu91OeXl5nXleou0J1eZEVSG31E61y0duiZ3uOtCZO+C0VeNTDegjU9A5DmMiGhd1Cwsf13HanDTllZKWIolSO1ezRDQu0gD4KPUY8ANGnYbYSD0uj5bSKndg+NX4+9tB67aiOPRgiIS4LmGLX7Qt6UPTm3Xitmhefr+fkpISnnzySWJjYxk3Tr4khUPNime/X8V5ZCWyw+0LJhONXfFc0+YELcE2J4oCXeIj2Vtix+Xx86tNS0dfJK5qJ9VaI/6CX8GSQoVqhpPIYbyRyfgi4olsB21OJFFq545eIlpa6aDaF/jl6hBlRCHwy+b1q3iPmkCn+FwYKvZCTBSkDwFt6xsKFUI0vfz8fDIyMujcuTNLliwJeVlSNK+Waiyu02rommDmpwMV7ChR2R81hr667SR4D2F0mzlEDw75o4jVuE98smO0pzYn8lvQztUsES2xuYKFJRP1Hoy6QMLk8fnRaRR0Nd9OCn/GUvA1enshqFWQeWG4QhdCHNHYy3cnq2vXrg26TCmaV2NrLZ0KvVZBAZxeH6YOfdhUNJpOlJCWGIXZ7+NMdrPd3/20fl/IZO52zmzQEmnQBlcXxOq8WHS/F0Eqr3aTEGXAbNCCsxJ+XobWWYaq0YNGA7990mzLRIUQQtSl12owGbQhb6dSaPhYdrcPm9NLZrIFRQGrGsl+bSfKu1+OT1WIsO5BW/wL23ftxe9vfUv3W4KMKLVzHp8KioJJr6Xa7aGD1odfBafHh9XuwGzUkR5vDkzkdpSh2IpQVD8oWkg6AxzlgVs96ityJoQQou2omZ6RZDaiU31UAlVeLev3lrNtVQE/fbMTl/c3ikoq6TfwbP546SUkNnNMzd2VorEkUWoj7LYqqu12qqIUYmKiG3Qfnz+w4s2k13Jmp2hUr4e8XVp8qgJuH6mxZtLjzcSYjgzxmuJRdSYUvxef3gTuKrB0ANPvtV48Pj8+tw+/X0tZ6nlEFX2Hq2ArZAQq97b0G1gIIcTJO7qCt8WoI1Hv4ef8UtZ+9gPOqnL0phg6GTx0oIzNGzeQn5fH9X119ExLaraYWlutJUmU2gC3283D8+dy8EA+Dz9wHxdeeOJ5Q6qqsr+sGqcn0KIkKyUGrd9NzDFLRGuVBIiIhrNuwV24A63PgRqZgNL/ulq1OcrsbiqcgUmGmuR+VCb0oEJnxl9sA1pnsTAhhBD1q6ngXWR1kmzWEqnx8NtPW6iqUOiUkcmunT5ilf3EmXXExyXwy4ESPt9gpUen5htXqpmj5feruL2By309kiy1Vv21JEmU2oik5GRi4+LJyRlQ735VVal2e/H5VWwuLzaXhyqnN7AENCESg06D31N3iWgdqX2xdxqK4nWSMOIusCSC+/cVD/FmA4nRNZMM665maOk3sBBCiJN3bAXvQ4etFBUdIjVjAA6PHx8aDpFI74RMNDEd6axa2PXzTg4WHCDZVXVK0y40bhsarx1cCph+v1LS2motSaLURhgMRgwGY70dpK0OD3klNrbnH0L1OjlcHoNfayQlxkTv1GgiDQ3433zoF0CBuG6oWiOq1linyisE3sAGQxN3dBRtXmubUyCEqO14c0prKnjnHVL5ze6lyuWjQ6QZi1ZHEV7cGHAYErEAkaoNj/UQ7n0b8az+PwyDb0TpmNP4gI5TlLK1kU+uNi7YnqTcRucdi+iy9Tmqi3dzuMpNqc3VsEJhHifsWgm//hsO/9bsMYv2p8zuZnexjZ2Hqvh8WxGfbyti56Eqdhfb2F1so8ze+DosrVHXrl155plngj8risJHH310SudsinMIcVwhGuceraaC94AkL8mRCjF6H+nxJoyKF4AqlxdrlY2y3Rsx4Mali6S8pIjSb5ZgrSg7qbC8kck4EvtCKy9KKSNKbdix7UkOR8RxkCRM8V3oFmPGr6r1tyc5Vv634HFAZAIk9my5JyDajZo5BT6fH2u1B6/fT2p0BJaIwCXe9npJtrCwsE5j21AWLFjARx99VKc5bmPOIcTJqq9x7rEURSEzOZre6fHsLDxIXGwsBsUPamCyd3lpKY7CInKSjZwZXY03Ih63rYy8vHwys6J+XxjUQG2lKKWMKLVhR7cncXr9FHrMeLURRJgiSYoyEhdpCLYnCclRDgc2Bv7dYxRo5LKaaDy9VoPb52dPiY29h23sLbGz7aCV3YdtuH3+VnXZze1uutGtlJQUjMZTK43RFOcQ4kRCNc49lkajcNHZvUhMSGT79u04nU60fg8VViu5+w4QHZfIxdlm9HiILt2MwRyDFQv5ZfZ2W5Sy9Xx6iUarqX/h8HgpqHDgQ0GvqKREG9Eo9bcnqWPPGvD7ID4D4ru1XPCiXQleAq5wEmHQEhdpwGzUUWR1sq3AitXhabbHHjFiBDNmzGDGjBnExMSQmJjI/Pnzgx/aXbt25aGHHuLGG28kOjqam2++GYBvvvmG8847D5PJRFpaGjNnzsRu/71NRHFxMZdddhkmk4mMjAzefPPNOo997GWzAwcOMGHCBOLj4zGbzQwaNIgNGzawZMkSHnjgAX766acjfbaUYEPfY8+xdetWRo4ciclkIiEhgZtvvhmbzRbcP3nyZMaPH88TTzxBamoqCQkJTJ8+HY+n+V5jcXrpmZbEbTOm079/fxwOByVlZdis5fTuP4hRN8witlM3UP2oihZvRCLRUTEn/lLehsmltzZMoyiUV7vx+VX0Cpi1fhL1bnSaEO1JjlWRD4d3gKJA95GB/wpxFL//xN8QVVUlr8RGlcNDSnRE8MPSqNWQHB1BYYWDvBLbiS8BH3EyK1uWLl3KtGnT+P7779m4cSM333wz6enp/O///i8ATzzxBPfddx/3338/AHv27GHs2LEsXLiQRYsWcfjw4WCytXjxYiCQkBw8eJA1a9ag1+uZOXMmxcXFIWOw2WwMHz6cTp068fHHH5OSksKmTZvw+/1cc801bNu2jc8//5xVq1YBEBMTU+ccdrudMWPGMGTIEH744QeKi4u56aabmDFjRjCxAlizZg2pqamsWbOG3bt3c80119C/f//g8xXiVGVn9+LOnlmUlpZS5fRy3h9voG/vTA5bq/nxuwqiPaX07aAlApXEgtVYE0Yc/0t5GyaJUhugqipetPhRsLm8GI0G3D4/RVYHOo1CRbWHzCQjGGpfUiivdpMaGxFoT1L3pLB7deDfqf0ChSWFOIrfr/LLwcoTHlft9rL9YCURBi12t4+Cckdwn0ZRcHl9FFmdOD3+Bq3APKNjdKOTpbS0NJ5++mkURSErK4utW7fy9NNPBxOHkSNHMmfOnODxN910E9dffz2zZ88GIDMzk+eee47hw4fz0ksvkZ+fz2effcb333/P4MGDAXj11VfJzs4OGcNbb73F4cOH+eGHH4iPjwegR48ewf0WiwWdTkdKSug5Im+99RZOp5PXXnsNs9kMwAsvvMBll13Go48+SnJyMgBxcXG88MILaLVaevXqxSWXXMLq1aslURJNSqPREBsbi1nV0CktDZ9fITnaSKXeQKmmE9uTcsguW42h9FcSlGh0Pdpnb1C59NbKWR0ethZYOVDpJa/Ezqfrt7J+Twlb9lfg8qp0ioskI8mM0+PH5VeC7UkKrce0JzmWokC3ERCbBl3PO2EcLpeLqqoqqqqk75uozedX8apqyHlIeo0Gr6ria8Do1Mk655xzar3PhwwZwq5du/D5AqNbgwYNqnX8Tz/9xJIlS7BYLMHbmDFj8Pv95Obmsn37dnQ6HQMHDgzep1evXsTGxoaMYcuWLeTk5ASTpJOxfft2+vXrF0ySAIYOHYrf72fHjh3BbWeccQZa7e9fgFJTU4872iXEqdDhJ8FspLzajYJCkt5NlM6H3ZzO3thzOewxYkzNqv9LeTsgI0qtmNXh4ZMvf+CLTz5i+7cr8Xk9LCnLJzo1gyEjL2LM0IHkpMdS7faRd6iCvb7jtCepT3xG4HYCBoOBMWPGsHbtWrZu3cqIESOa7kmKVkujUTij44nb5dhcXpweH2ajDuNRyVL3RDOKRsHp8RHnCtRpsRhP/JHTHEXljk48IHCZ7JZbbmHmzJl1jk1PT2fnzp2NfgyTyXTS8TWWXl/791pRlNO2YenpqqZ2md+v4vQEvhA43L5a1aubahGFokB6vAmHz0FhpQNFVYjTeXDpNezSZmLp1J1BUUkNurTeWKGKUrYkSZRaKVVVWfPdJv61+B84KsswRphAH09EVCzFudtZ+W4RWSlm+qWdHax/EXW89iRHKD4XiqcaGllRdeDAgWRlZWGxtO5lnKJpNSRpiYrQkRhlDLRAiI4I9mNSNAoaRcHq8JAaG0FURIhq8E1gw4YNtX7+7rvvyMzMrDXqcrQBAwbw66+/1ro0drRevXrh9Xr58ccfg5feduzYQUVFRcgY+vbtyz//+U/KysrqHVUyGAzBEa5QsrOzWbJkCXa7PZjcrVu3Do1GQ1ZW1nHvK04vNf3QgGD16r0lvy9GaOp2UtEmPWd2MpB3SA1+KY/Xa8lONaHTanB5VA6UV9OJYlRFj92nwetXMLm8ROkMJ/e730qKUsqlt1aqyunhk4//jaOygq7de6IYIlEUMJtM9MrujeqoYvl//kOVM7DSRVEC7Uli9L7Q7UkKf8ZS8A2xuZ/Cf2ZCweYGxxMVFUVqamq9lcHF6a2mBYLZqKOwwoHL6wt+yz3hJeAmkp+fz+23386OHTv417/+xfPPP8+sWbNCHv+Xv/yFb7/9lhkzZrBlyxZ27drFv//9b2bMmAFAVlYWY8eO5ZZbbmHDhg38+OOP3HTTTccdNZowYQIpKSmMHz+edevWsXfvXt5//33Wr18PBFbf5ebmsmXLFkpKSnC5XHXOcf311xMREcGkSZPYtm0ba9as4bbbbmPixInB+UlCQKB2WY8OlpC3eLOhyR+z5kt5v2gbfaNsDOoSy7DMRLJSAn8X7Ad3UPzNUgq/fYvfylR2l3nYvPsAW09h5WuFEkuethtV5vCtypZEqZXKy8tj397dpHbqhFZR0VcfQm8vJNpfjlGvJbVTJ/L27CIvL69hJ3RWws//QucoQdVoUSoPwk9vBbYLcYpqWiCkxkbgdPsod7gDhVBjIzizU0yjC9E11o033ojD4eCss85i+vTpzJo1K1gGoD59+/blyy+/ZOfOnZx33nnk5ORw33330bFjx+AxixcvpmPHjgwfPpwrrriCm2++mQ4dQi96MBgMrFixgg4dOnDxxRfTp08fHnnkkeCo1pVXXsnYsWP5wx/+QFJSEv/617/qnCMyMpIvvviCsrIyBg8ezP/8z/8watQoXnjhhVN4dUR7pNdqMBm0IW/NVbusvi/l8WYD6QmRWHWJ7LdpsR/4hbNKPmBA5Wqytj+HPe/HkyoT4na7uevV1dz87Gds2PJrszyfhpBLb62Uo9qOx+XEEBEJiooOHygqWnMCAAajCY/bhaPafoIz1ZywDKUsF1BB0UJKX3BUBApORoTnuq9oXwLfNmNwevz4/Cpndopp1sttR9Pr9TzzzDO89NJLdfaF+jIxePBgVqxYEfKcKSkpfPLJJ7W2TZw4sdbPxxbY69KlC++991695zMajfXuO/Ycffr04b///W/IuI4uE1Dj6LYqQpyKWj3hDA2/ghAdoUPVR7I77jxGH1qB1luNX6ND7yon8+DH/GTpQn6ZrsFlQmokx1mIizKRk9P/JJ5N05ARpVaqQ3wsUZZIKiqPrDJTFFA0oA18M6+oqiLKbKJDfGzDTuj1oDorUHxufHozuCrBFBe4CXGKPD4/DrcPp8cfLHaqVRScnsB2TzutryJEu9KAnnCh2N0+HG4f3ZItuI3xuJQIVK8X1etC67KSqK0+qaKUBr0Oi6n+hvAtRUaUWqkuXbqQ0+cMvl7/A0ZjV1Q0KPjx+vxUezyUFxdx/rln0aVLlxOfzOOEPatQks/AV1UCKKiRCSj9rwuOJgVXUHh8eI8s43a4fWjUwJtaur+L42npiaVCiObRkJ5w9anpFGGK7oA2tjOeigPo/E60Fbl4YrqgRMTi9ZygU0QrJYlSK6XRaJhw9f9w4MB+9u3fi9PjQ6vTU2GtxFZeQre0VK696ko0mhMkL6oKO5aDowI1sSeVHhOK34NxxF1gSQweVvOHTvG5STqy0im3xI6qDVxTlj904nhqmuKG0pxNcdeuXdts5xbidOM3WPAbLIGecN6GzynSaTXoNApOjRlrzysxFW3H7wKv1sjBmIHssymYjITuFNGKSaLUimVnZ3Pn7X/mnWX/YtlvP+L2q+g9dv4w7BzGjx9/3CrBtcSkQVku9BqHWrIYFeo0Rgz+ofO50R1ZLRHdwQLawL/ba/d30TT0Wg369llrTgjRAGaDlniLgSKrE0PCmRxKGI7e70DXaywHDN0prnSTGqsJfHknMDevSUoItABJlFq57Oxs7rxjDgc2fYHb4+Pu++8ls1fvE48k1VAUSDsLks8MTOIOIfiHzqsF1QM+D3q/A1qwiJ4QQoi2qaZMSJXTGyhKqUSg1xkxJ2ShrfYRbdKRaDZwcO+vaJO6YbM72VtpwacqHNhXQWKMhy4JJyiSHCaSKLUBGo2GGHPgsleXLl0aliR5nIEkSXekqKQhEjx167bUoTNAj1GQ901gMl/GidubCCGEEDVlQo4uSonbR9dEM51jIvD98m8cBdvYmj+Ew1E9SdL4iNComA1aiqxOqpzeFikn0liSKLVHqgo7PoHqUjjj8sY3vO3YHxIzwWA+4aFCCCFEjZCdIgA1sQP7D0DyoS8x6LUUug0kGrwk6LVEmgwUWh3kl9kbXUKguUmi1B4VbISSXaDRgr9xSzGBwPwlo1TgFifJVQVueyDRlveREKedmqKUaKnVKcLeaRgV+w6TxM9YitfgcPmwe80UFh8mOr4DsZH6YAkBi1GHqqp4FAOqosHm8mKMVMOSQEmi1I4oPhdaRynsWgF6U+ASWnRquMMSpxOvGz69Eyry4bw7occfwh2REKKV8PpVipPPJ8oACftWca71G+yKBfcvReR3uRJbXG+8fhW3x4fVr7Kr0EapsRN+RcOm/EpSqpWwzGNqe+v0RP0Kf8ZcsI7oA2tQcr8KJEodB4Q7KnE6sqRAag6k9gt3JK3aggUL6N+/f7jDYMSIEcyePTvcYYhmUFMINlAMNnCr+TkchWB1Wg06rYayxMHoqwsxqi6isBHnKqDr/g9xV1updHjYWmDl652HOWh1oVPdRPhsRBo1FFmdJ9UK5VRJotQeOCvhpzfRVx9CVY5cbrPuD1wCEaKl6YxgtARuLaCoqIhZs2bRo0cPIiIiSE5OZujQobz00ktUV1e3SAzNYe3atSiKQkVFRas8n2j9yuxudhfb2FtiJ0KvJUKvZW+Jnd3FNnYX2yizu1s0npoSAtWVJfgMsTj0saDRokalEqupJtJfSYLFQKnNxcEKJw63l3K7m7LySsoOHSQ5yojd5SW/zF6n9U9zkktv7YGjDKVkFygKqkYHnQYEkifp4ybaub179zJ06FBiY2N5+OGH6dOnD0ajka1bt/Lyyy/TqVMnxo0bV+99PR4Pen3rWl1zMtxuNwZD03eKF21fcxWCrdUPTm9s+P2OlBDYXplIlWJBr+hxaaOIcNuo0scRn5BMemocm/IrcJXs4sOVy9my5Wd8Ph+upx+nZ6/eDB19EQZt9+A8ppYgI0rtgSkeNaEHqCpeYyw4rdLHTYSP1xVouNwCI5q33norOp2OjRs3cvXVV5OdnU23bt344x//yPLly7nsssuCxyqKwksvvcS4ceMwm8383//9HwAvvfQS3bt3x2AwkJWVxeuvvx68T15eHoqisGXLluC2iooKFEUJVgSvGalZvXo1gwYNIjIyknPPPZcdO3bUivWRRx4hOTmZqKgopk2bhtPpDPm88vLy+MMfAvO74uLiUBSFyZMnA4FLZTNmzGD27NkkJiYyZsyYE8Z5vPMB+P1+5s6dS3x8PCkpKSxYsKCh/wtEK6bXajAZtCFvJ9WW6hT6wUFgVVx2RmdcZ1yNQ2PB5/Pi0Mfi7nMtvTsYiSr+nkP7dvPVR29QtPc3IgwGoqOjiEuIZ/u2n3h36Svs2bWjRVuhSKLURrg9XmwOF1VVR/3xqRl6jIiG/jfgMaeg8XtRIxPgqD5uQpwSrzv0zeetfWzBj7D3S8j9Ev77EOz/4ahjPSc+byOUlpayYsUKpk+fjtlcfymLY1fILFiwgMsvv5ytW7cydepUPvzwQ2bNmsWcOXPYtm0bt9xyC1OmTGHNmjWNigXgnnvu4cknn2Tjxo3odDqmTp0a3PfOO++wYMECHn74YTZu3EhqaiovvvhiyHOlpaXx/vvvA7Bjxw4KCwt59tlng/uXLl2KwWBg3bp1/P3vfz9hbA05n9lsZsOGDTz22GM8+OCDrFy5stGvgTg9OBL7UJZ9PaT2P6n7x5j0dO8zBE36Wfg6n03MRffRo/dgone+j37f1/z07xcpKykhMyuLSL2KEQ8WSxTdevairLSU79Z8jqYFF7/Jpbc2wGAwMCKnB99uy2Prtm2MTOkY+KOz7QNI6A6dB0FqX+wdz0XxOuv0cZOGt+KUfP1k6H0J3aHv1YF/OythxT1gzQ+0vtm3Hgp/hozzA/OWYtMh5/rf7/vdi+Bx1D7fH+5qcFi7d+9GVVWysrJqbU9MTAyO1kyfPp1HH300uO+6665jypQpwZ8nTJjA5MmTufXWWwG4/fbb+e6773jiiSeCIzAN9X//938MHz4cgHnz5nHJJZfgdDqJiIjgmWeeYdq0aUybNg2AhQsXsmrVqpCjSlqtlvj4eAA6dOhAbGxsrf2ZmZk89thjwZ/z8vKOG9uJzte3b1/uv//+4LlfeOEFVq9ezQUXXNCg5y7aL5fLhdvtpqqqioSEBOCYfnAnSVEUIg16MOixRMeh6I3Q5VwKv3qbor2/khqXioJa5z4xSckU5e+htKiA6IyMU3puDSV/HduIfj1SuXHMQAYOGHAkSXofyvbC3rXBSxyq1ojfGFPnzVszoS+3xI5Oo6DTKOSGcUKfaKccZYGK8FpDoIaXMQq8TvA6TnzfJvT999+zZcsWzjjjDFyu2tXoBw0aVOvn7du3M3To0Frbhg4dyvbt2xv9uH379g3+OzU1UJajuLg4+Dhnn312reOHDBnS6MeoMXDgwJO+b32Ojh0C8dfELk5fBoOBMWPGEBUVxdatW5v/ATsNxJY6BJ9fJV5jRzm8A53iJ1LjwedxYnN6iYqyoFO92Gy25o/nCBlRaiMsJiMWk5GoyIgjSVIuaPWBb/PGqOO2J5GGt+KUnDcn9D7lqO9apvhAX8F96wPvyeiOYEmC8+cGfj62UNw5t55SWD169EBRlDpzgbp16xYIp54+haEu0YVS0y7o6BU2Hk/9S5OPnhhec8nP72+eeRTHPo/GxFmfYye1K4rSbLGLtmXgwIFkZWVhsbTMKtaorjmYU3qA9zDR7iIGmQ/hUrVEFW7Al9IHnz4K1RxJVFTLFbOVEaU2QvG50DjL4Ke3aidJseknvG9wQp9ee6T5rQaT/hQn9InTh84Q+qY96rtWRDT0vz5QkdtTHUiSciaCOeHIsfoTn7cREhISuOCCC3jhhRew2+0n9dSys7NZt25drW3r1q2jd+/eACQlJQFQWFgY3H/0hOnGPM6GDRtqbfvuu++Oe5+alWw+34mr6zckzsacT4gaUVFRpKamtlhikp6eTq9+gyl0RxLhKiVScaFVvVg0buIqf+NwQT7Z2dmkp5/4b19TkRGltqDwZ8wH16GvPoxSsRU6D4ShsxqUJAnRolL7BeYkeR2BkSRzQrM+3IsvvsjQoUMZNGgQCxYsoG/fvmg0Gn744Qd+++23E16iuvPOO7n66qvJyclh9OjR/Oc//+GDDz5g1apVQGBU6pxzzuGRRx4hIyOD4uJi7r333kbHOWvWLCZPnsygQYMYOnQob775Jr/88ktw9Ks+Xbp0QVEUPvnkEy6++GJMJlPIb/UNibMx5xMiXDQaDZdffjn7d//Kr79sx+TR4NcoVPkMHMwtIjG7O+PHj29Yc/imiqnFHkmcHGcl/LwMncsKqIE5H+5qiIgFalde9fpVvH41rJVXhUBnDLw/W6DPW/fu3dm8eTOjR4/mrrvuol+/fgwaNIjnn3+eO+64g4ceeui49x8/fjzPPvssTzzxBGeccQb/+Mc/WLx4MSNGjAges2jRIrxeLwMHDmT27NksXLiw0XFec801zJ8/n7lz5zJw4ED27dvH//t//++49+nUqRMPPPAA8+bNIzk5mRkzZhz3+BPF2djzCREu2dnZzJw5i5yeaVQ4/RyscFNWXsaAnmnMnDmL7OzsFo1HUVuyvOVpoLKykpiYGKxWK9HRTbA8vzwPdfVD2Iv3ovU6iUgfhOJzwsh7Ia4LhyqdFFe6UHxukra8AMDh/jNQj8w/6hBtJDk6InAurxvWPByYDH7eHDDHn3p8ol1wOp3k5uaSkZFBRETEyZ/I6/59ldx5cxp9OU20PU323hFhE1wZ7Vf5tbASgN6p0WiOrME/emW036/yy8HAMWd0/P2Yo/k9Lg5/+jAASRffjaaeopQNOcax+xv+s/A6VJ+HgedfRLeLZqDp3DStuRrzt1ouvbViHp8fry4GrSEajduGV2/B66wESxJeXTQ6n79xE7V1hkCj3LxvAoXCMs4Lx9MS7Z3XFUjGXTbQSTIuRGtXZndTXBlYEBSh1wKwt+T3eX+1vnC3ICW1Hz+504nUevnjlQvRxHds8RhAEqVWLfDm1WDqNI7U/B/ReauxKtGUdx6P06qhg+omOToCvRbwauFIxq/Xa0Gnrf+kHftDYmZgwq0QTU2ScSHanOZodXKybU6O5fJrcfm1LXIpPxRJlFqx4Js34Rx0h0eA14l52B2YIwMTZE9qWb8xKqxvOHEakGRciDYlsBq6CU94pM2JuXB94AuT5fwmPHnLk0SpFQu+eb1aMJrAaEJviQ09WiREayDJuBCnPUdiH1yx3UhKTQl3KKdMEiUhhBBCNKmmaHPSWkh5ACGEEEKIECRREkIIIYQIQS69CSGaVFVVFTabDYvF0qL9mIQQojnIiJIQosm43W5uu+02JkyYcMJeZkKItuPoLhBOT+B2unSBkBElIUSTSk1NJSEhgQEDmqaCrghtxIgR9O/fn2eeeSbcoYh2rrUWpWwJMqIkhGgyfr8fh8NBdXU1paWl+P3N/y1z8uTJKIrCn/70pzr7pk+fjqIoTJ48udnjEKI9izcb6NHBEvIWb258uyLF50bjsQeKUobg9nixOVxUVVWdSvinRBKl9sTrCrzhjvOmE6K5bN++nccff5w1a9bw1Vdf8eCDD/LII4+wffv2Zn/stLQ0li1bhsPhCG5zOp289dZbpKenN/vjnwq32x3uEIQ4Ib1Wg8mgDXmr6QXXYEeKUvr15kBRynoYDAZG5PTAYjKyddu2JngWJ0cSpfaipnWE0RLyTSdEc9m+fTvPPfccW7ZsITIykoSEBBISEti8eTPPPfdcsydLAwYMIC0tjQ8++CC47YMPPiA9PZ2cnJzgNr/fz1//+lcyMjIwmUz069eP9957L7jf5/Mxbdq04P6srCyeffbZWo+1du1azjrrLMxmM7GxsQwdOpR9+/YBgdGt8ePH1zp+9uzZjBgxIvjziBEjmDFjBrNnzyYxMZExY8YAsG3bNi666CIsFgvJyclMnDiRkpKS4P3sdjs33ngjFouF1NRUnnzyyVN+3YQIJ0diH8qyr4fU/iGP6dcjlRvHDGRgGC/lS6LUnnTsD4OmBP4rRAvx+/18+OGHlJSUkJ2djdFoRKPREB0dTe/evSkpKeGjjz5q9stwU6dOZfHixcGfFy1axJQpU2od89e//pXXXnuNv//97/zyyy/8+c9/5oYbbuDLL78MPpfOnTvz7rvv8uuvv3Lfffdx991388477wDg9XoZP348w4cP5+eff2b9+vXcfPPNKErj2gktXboUg8HAunXr+Pvf/05FRQUjR44kJyeHjRs38vnnn3Po0CGuvvrq4H3uvPNOvvzyS/7973+zYsUK1q5dy6ZNm0725RIi7PwGC97I5OMWpbSYjCTHR4V1Ba1M5m5PpHWECIP8/Hx+++030tLS6iQMiqLQuXNntm/fTn5+Pl27dm22OG644Qbuuuuu4OjOunXrWLZsGWvXrgXA5XLx8MMPs2rVKoYMGQJAt27d+Oabb/jHP/7B8OHD0ev1PPDAA8FzZmRksH79et555x2uvvpqKisrsVqtXHrppXTv3h2A7OzsRseamZnJY489Fvx54cKF5OTk8PDDDwe3LVq0iLS0NHbu3EnHjh159dVXeeONNxg1ahQQSLY6d+7c6McWQjSOJEpCiFNSVVWF0+nEbK6/Ca7ZbKagoKDZJ2MmJSVxySWXsGTJElRV5ZJLLiExMTG4f/fu3VRXV3PBBRfUup/b7a51ee5vf/sbixYtIj8/H4fDgdvtpn///gDEx8czefJkxowZwwUXXMDo0aO5+uqrSU1NbVSsAwcOrPXzTz/9xJo1a7BYLHWO3bNnTzCOs88+O7g9Pj6erKysRj2uEM3J4/Pj9an4/SpOjw8Ah9uHRhP4AqXTKo2fy9QKSKIkhDglUVFRREREYLfb602W7HY7ERERLTJ0PnXqVGbMmAEEEp6j2WyBRQ7Lly+nU6dOtfYZjUYAli1bxh133MGTTz7JkCFDiIqK4vHHH2fDhg3BYxcvXszMmTP5/PPPefvtt7n33ntZuXIl55xzDhqNBlVVa53b4/HUifPY18lms3HZZZfx6KOP1jk2NTWV3bt3N/QlECJs2msJAUmUhBCnJD09nV69erF58+Y6IxyqqnLgwAEGDBjQIqvPxo4di9vtRlGU4CTpGr1798ZoNJKfn8/w4cPrvf+6des499xzufXWW4Pb9uzZU+e4nJwccnJyuOuuuxgyZAhvvfUW55xzDklJSWw7ZnXOli1b0Ov1x417wIABvP/++3Tt2hWdru7Hcvfu3dHr9WzYsCH4OpaXl7Nz586Qz0WIlhZvNhAdEfq9rtMqbXLUSRIlIcQp0Wg0XH755ezfv5/t27fjdDoxGAxYrVYKCwtJTExk/PjxaDTN/+Gn1WqDK+y0Wm2tfVFRUdxxxx38+c9/xu/3M2zYMKxWK+vWrSM6OppJkyaRmZnJa6+9xhdffEFGRgavv/46P/zwAxkZGQDk5uby8ssvM27cODp27MiOHTvYtWsXN954IwAjR47k8ccf57XXXmPIkCG88cYbbNu2rdalvfpMnz6dV155hQkTJjB37lzi4+PZvXs3y5Yt45///CcWi4Vp06Zx5513kpCQQIcOHbjnnnta5DUVoqH0Wg167fGPOVTpPOGoU7zZEEimPD68/sAIrcPtQ6MGEquWTqYkURJCnLLs7GxmzpzJe++9x+7du6mqqqKsrIwBAwYwfvz4k5rwfLKio6ND7nvooYdISkrir3/9K3v37iU2NpYBAwZw9913A3DLLbewefNmrrnmGhRFYcKECdx666189tlnAERGRvLbb7+xdOlSSktLSU1NZfr06dxyyy0AjBkzhvnz5zN37lycTidTp07lxhtvZOvWrceNuWPHjqxbt46//OUvXHjhhbhcLrp06cLYsWODydDjjz8evEQXFRXFnDlzsFqtTfGSCdFiGjLqVHMJT/G56aB6UPwe9hUewhcRD7T8JTxFPfaCujgllZWVxMTEYLVaj/uB3SheN3x9pGbKeXMCNZOEaEJOp5Pc3FwyMjKIiDj5DyCn08m8efNwuVz8+c9/pkePHjLq0c411XtHiBo1l+cANPu+QbPvG/xdhuHvMgxomhGlxvytlhGltsLrAp8nUHVbFx/uaISol0ajwWQyodVqSUhIkCRJCNFotS7hpQ+AlCwwmMFwgut6zUQ+xdoCqbot2giDwcCYMWOIioo64eUmIYQ4IWMURKWEtUagjCi1FR37Q2JmIKsWohUbOHAgWVlZ9dYEEkKItkYSpbZCqm6LNiIqKrztBoQQoinJpTchhBBCiBAkURJCBMkiWNFY8p4R7Z0kSkKIYOXo6urqMEci2pqa98yJqo8L0VbJHCUhBFqtltjYWIqLi4FAYUVFUcIclWjNVFWlurqa4uJiYmNj61RCF6K9kERJCAFASkoKQDBZEqIhYmNjg+8dIdojSZSEEAAoikJqaiodOnSot+O9EMfS6/UykiTaPUmUhBC1aLVa+eMnhBBHyGRuIYQQQogQJFESQgghhAhBEiUhhBBCiBBkjlITqym+VllZGeZIhBBCCFGfmr/RDSmYKolSE6uqqgIgLS0tzJEIIYQQ4niqqqqIiYk57jGKKvXnm5Tf7+fgwYNERUU1acG+yspK0tLS2L9/P9HR0U123vZIXquGkdep4eS1ahh5nRpOXquGaa7XSVVVqqqq6NixIxrN8WchyYhSE9NoNHTu3LnZzh8dHS2/VA0kr1XDyOvUcPJaNYy8Tg0nr1XDNMfrdKKRpBoymVsIIYQQIgRJlIQQQgghQpBEqY0wGo3cf//9GI3GcIfS6slr1TDyOjWcvFYNI69Tw8lr1TCt4XWSydxCCCGEECHIiJIQQgghRAiSKAkhhBBChCCJUhtVXFzMa6+9xnXXXUdmZiYRERFERkbSq1cvZs6cSV5eXrhDbHHvvvsuI0aMIC4uDrPZTL9+/XjsscfweDzhDq1V8Hg8rF69mjvvvJPBgwcTGxuLXq8nJSWFcePGsXz58nCH2KrNnTsXRVFQFIWFCxeGO5xWxe1289xzzzFs2DDi4+OJiIigc+fOXHTRRbz99tvhDq/VyM/PZ8aMGWRlZWEymYiIiCAjI4NJkybx008/hTu8FrNjxw6ef/55Jk+eTJ8+fdDpdA3+vVq1ahUXX3wxiYmJmEwmevXqxT333IPNZmu+gFXRJl1//fUqoGo0GrVv377qVVddpV588cVqUlKSCqhms1ldsWJFuMNsMbNmzVIBVafTqRdeeKF6xRVXqLGxsSqgDhs2TK2urg53iGG3cuVKFVABNSUlRb3kkkvUq6++Wj3zzDOD22+++WbV7/eHO9RWZ926dapGo1EVRVEB9aGHHgp3SK3G/v371d69e6uAmpiYqF566aXqNddco5577rlqZGSkeuWVV4Y7xFbhu+++U6OiolRA7dSpkzpu3Dj18ssvVzMyMoKfXe+88064w2wRNZ/Xx95O9Hv11FNPqYCqKIp6/vnnq1dddZWakpKiAmpWVpZ6+PDhZolXEqU26rbbblMfeOAB9cCBA7W2V1VVqddee60KqPHx8WpZWVmYImw5H374oQqoFotF/fHHH4PbDx8+rPbp00cF1Dlz5oQxwtZh9erV6pVXXql+9dVXdfYtW7ZM1Wq1KqAuXbo0DNG1Xna7Xc3MzFQ7deqkjh8/XhKlo1RXV6u9evVSAXXBggWq2+2utd9ut6ubN28OT3CtTN++fYNfRo5+nXw+n3rvvfeqgBobG6s6HI4wRtkyXnnlFfWOO+5Q33zzTXX79u3qxIkTT/h7tWnTJlVRFFWr1aqffvppcLvdbldHjRqlAs2WlEui1A7Z7fbgN5fXX3893OE0u8GDB6uAunDhwjr7vv76axVQjUajWlFREYbo2o5p06apgDpq1Khwh9KqzJw5UwXU5cuXq5MmTZJE6Sjz588P/vEXoZWUlARHTYqLi+vs93q9qslkUgF106ZNYYgwvBrye3XVVVepgHrTTTfV2ZeXl6dqNBoVULdv397k8ckcpXYoMjKSrKwsAPbv3x/maJpXQUEBP/zwAwDXXXddnf3Dhg0jLS0Nl8vFp59+2tLhtSk5OTlA+3/PNMbatWt5/vnnufHGG7n44ovDHU6r4vF4eOmllwC48847wxxN69aYGkCJiYnNGEnb5Ha7g3Mo6/uc79KlC0OHDgXgww8/bPLHl0SpHfJ4PMHJ3KmpqeENpplt3rwZgPj4eDIyMuo9ZtCgQbWOFfXbtWsX0P7fMw1ls9mYOnUqycnJPPPMM+EOp9XZtGkTJSUldOzYkR49erB161YeeOABbrnlFubNm8fy5cvx+/3hDrNVsFgsnHfeeQDce++9tRaY+P1+FixYgMPh4KKLLiItLS1cYbZaO3fupLq6Gvj98/xYzfk5L01x26FXX32VkpISTCYTF110UbjDaVa5ubkApKenhzym5oOn5lhRV1FREUuWLAHgyiuvDG8wrcQdd9xBbm4uH374IXFxceEOp9X5+eefAejcuTPz5s3jscceQz2qfvGjjz5KTk4OH3300XF/P08Xr7zyChdffDEvv/wyy5cvZ9CgQWi1WjZv3kxBQQETJ07khRdeCHeYrVLNZ3dsbCxRUVH1HtOcn/MyotTObN26NTgMPn/+fJKTk8McUfOqqqoCwGw2hzzGYrEAUFlZ2SIxtTVer5cbbrgBq9VKnz59uOWWW8IdUtitWLGCf/zjH1x77bWMHz8+3OG0SqWlpUDgG/yjjz7Krbfeyo4dO7BaraxcuZKePXuyefNmLrnkEinRAWRlZbF+/XouvPBCCgoK+Pe//80HH3xAbm4uPXr0YMSIEURHR4c7zFYp3J/zMqIUBnPnzuXjjz9u9P3++c9/MmzYsJD7Dxw4wGWXXYbNZmPcuHHMmzfvVMIUp4k//elPrF69moSEBN577z0MBkO4Qworq9XKtGnTSEpK4vnnnw93OK1WzeiRx+NhwoQJtUZDRo8ezcqVK8nKymLbtm0sW7aMiRMnhivUVmHdunVcccUV6HQ63nrrLUaOHInBYGDdunXcfvvtTJs2jXXr1vHqq6+GO1RxDEmUwuDgwYPs2LGj0fc7XkGtoqIiRo0axb59+xgzZgzvvPMOiqKcSphtQs0wrN1uD3lMzesm39bqmjVrFq+++ipxcXHBUYDT3ezZszlw4ABvv/22TKw9jqMvgdQ3Cpmens4ll1zC+++/z6pVq07rRKmiooLLL7+ckpIS1q9fz9lnnx3cd+mll9K7d2/69OnDokWLuOGGG/jDH/4Qxmhbn3B/zkuiFAZvvPEGb7zxRpOdr7i4mJEjR7Jz505Gjx7NRx99dNp0pO7atStw/JVaNftqjhUBc+bM4bnnniM2NpYVK1YEV72d7j788EN0Oh0vvvgiL774Yq19v/32GxCYB7hq1SpSUlJYtmxZOMIMu27dutX77/qOKSwsbJGYWqvly5dz+PBhunfvXitJqtGtWzfOPvts1qxZw6pVqyRROkbNZ3dFRQVVVVX1zlNqzs95SZTauMOHDzNy5Ei2b9/OqFGj+Pjjj4mIiAh3WC2m5o97aWkpubm59a5827hxIwADBgxo0dhas7lz5/LUU08RExPDihUrQq4kOV15vV6+/PLLkPvz8vLIy8ujS5cuLRhV6zJgwAAURUFVVUpKSupdrVVSUgL8Pn/kdJWfnw8cf7QjJiYGgLKyshaJqS3JysoiMjKS6upqNm7cWG8i2Zyf8zKZuw0rKSlh5MiR/PLLL4waNYr//Oc/mEymcIfVojp37szgwYMBeOutt+rs/+abb9i/fz9Go1Hq4Bwxb948Hn/8cWJiYli5cmXw9RMBFRUVqIFivHVukyZNAuChhx5CVdXTsqdijZSUlOCcyVWrVtXZ7/F4gsnmWWed1aKxtTadOnUCAiOSVqu1zn6Px8OmTZsAQpY5OZ0ZDAYuueQSoP7P+X379vHtt98CcPnllzd9AE1ewlK0iNLS0mBJ/NGjR5/WvcxCtTApKSmRFibHuOeee4KtEr7//vtwh9PmSGXu2latWqUCalxcnLp+/frgdo/Ho952220qoEZFRalFRUVhjDL8iouLVbPZrALqVVddpVZVVQX3uVwudfr06Sqg6vV6dc+ePWGMNDwa8nv1448/BluYfPbZZ8HtLdHCRFHVowpfiDbjiiuu4MMPP0RRFK666qqQI0njx48/LZY3z5o1i+eeew69Xs+oUaMwm82sXr2aiooKhg4dysqVK0+70bZjffzxx/zxj38EAsXZzjjjjHqPS0xM5IknnmjJ0NqMyZMns3TpUh566CHuvffecIfTKixcuJD58+ej0+k466yzSElJYdOmTeTl5WEymXj33XeDowGnszfeeIMpU6bg9XpJSkpi8ODB6PV6Nm7cSEFBARqNhr/97W/86U9/CneozW7Tpk3ceuutwZ/37NlDSUkJnTt3Do6+QWC+4NEFcJ9++mluv/12FEVh+PDhdOjQga+//prCwkKysrL45ptvmmcBRrOkX6LZDR8+vN7uy8fe7r///nCH2mLefvtt9fzzz1ejo6NVk8mknnnmmeojjzyiulyucIfWKixevLhB75kuXbqEO9RWS0aU6vfFF1+oF110kRofH6/q9Xo1LS1NnTx5crP03WrLtmzZok6ePFnt1q2bajQaVYPBoHbp0kW9/vrr1Q0bNoQ7vBazZs2aBn0W5ebm1rnvypUr1bFjx6rx8fGq0WhUMzMz1bvuukutrKxstnhlREkIIYQQIgSZzC2EEEIIEYIkSkIIIYQQIUiiJIQQQggRgiRKQgghhBAhSKIkhBBCCBGCJEpCCCGEECFIoiSEEEIIEYIkSkIIIYQQIUiiJIQQQggRgiRKQgghhBAhSKIkhBBCCBGCJEpCCCGEECFIoiSEEEIIEYIkSkIIIYQQIUiiJIQQR3n00UdRFAWDwcD3339f7zGffvopGo0GRVF48803WzhCIURLUlRVVcMdhBBCtBaqqnLhhReyatUqunXrxpYtW4iKigruLywspF+/fhw+fJgbb7yRpUuXhjFaIURzk0RJCCGOUVRURL9+/SguLub666/njTfeAGonUT169GDz5s1YLJYwRyuEaE5y6U0IIY6RkpLCkiVLgpfWakaNHn30UVatWoVer+df//qXJElCnAZkREkIIUKYM2cOTz31FBaLhZdeeompU6fi8Xh4/PHHueOOO8IdnhCiBUiiJIQQIbjdboYMGcKmTZuC2y688EI+//xzFEUJY2RCiJYiiZIQQhzHtm3b6NOnDwAxMTH89ttvpKSkhDkqIURLkTlKQghxHC+//HLw35WVlWzZsiV8wQghWpyMKAkhRAiffPIJl112GQB9+/bl559/pkOHDvz8888kJyeHOTohREuQESUhhKhHYWEhU6ZMAWDKlCl89dVXdO3aleLiYiZNmoR8xxTi9CCJkhBCHMPv9zNx4kRKSkrIzMzk+eefJyYmhrfeegudTscXX3zBU089Fe4whRAtQBIlIYQ4xmOPPcbq1auD9ZLMZjMAQ4YM4f777wfg7rvvrrUaTgjRPskcJSGEOMr333/PsGHDQtZL8vv9jBo1irVr19KzZ082bdoUTKSEEO2PJEpCCHFEVVUV/fv3Z+/evVxwwQV88cUX9dZLOnDgAP369aOsrIzJkyezePHiMEQrhGgJkigJIYQQQoQgc5SEEEIIIUKQREkIIYQQIgRJlIQQQgghQpBESQghhBAiBEmUhBBCCCFCkERJCCGEECIESZSEEEIIIUKQREkIIYQQIgRJlIQQQgghQpBESQghhBAiBEmUhBBCCCFCkERJCCGEECIESZSEEEIIIUKQREkIIYQQIoT/D1mauvaLfp7LAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "query_ei = wf.optimize('EI', n_restarts_optimizer=10)\n", "query_qbc = wf.optimize('max_uncert', n_restarts_optimizer=10)\n", "display('Query By Committee (maximal uncert):', query_qbc, 'Expected Improvement', query_ei)\n", " \n", "# Initialize plot figure\n", "fig = plt.figure(figsize=(6,4))\n", "ax = fig.add_subplot(111)\n", "\n", "\n", "n_samples = 50\n", "X = np.linspace(-2, 10, n_samples)\n", "x_input = pd.DataFrame({'x':X})\n", "pred = wf.predict(x_input, X_scaled=False, min_max=True, return_std=True)\n", "rbmal_queries, _ = wf.RBMAL(pool=x_input, batch_size=5)\n", "ax.errorbar(X,pred['y'],yerr=1.96*(pred['y_std']),xerr=None,label='prediction', fmt='o-', alpha=0.2,capsize=3)\n", "ax.errorbar(X,f(X),yerr=1.96*(sigma)/np.sqrt(3),xerr=None,label='Ground truth', fmt='.--', alpha=0.5,capsize=1)\n", "ax.errorbar(combined_df['x'],combined_df['y'],yerr=1.96*(combined_df['y_std_mean']),xerr=None,label='Measured', fmt='o', c='black', alpha=0.5,capsize=1)\n", "ax.scatter(query_qbc['x'],query_qbc['y'],label='QBC query', marker='*', alpha=1, c='r',s=60)\n", "ax.scatter(query_ei['x'],query_ei['y'],label='EI query', marker='v', alpha=1, c='g',s=60)\n", "ax.scatter(x_input.loc[rbmal_queries]['x'],pred.loc[rbmal_queries]['y'],label='RBMAL query', marker='.', alpha=1, c='purple',s=60)\n", "\n", "\n", "ax.set_ylabel('y', fontsize=18)\n", "ax.set_xlabel('x', fontsize=18)\n", "ax.tick_params(labelsize=16)\n", "ax.legend()\n", "plt.tight_layout()\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.18" } }, "nbformat": 4, "nbformat_minor": 4 }