Skip to content

Instantly share code, notes, and snippets.

@RyotaBannai
Last active January 12, 2019 09:22
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save RyotaBannai/0e94870c7eb9e7305e446fc08a9090d6 to your computer and use it in GitHub Desktop.
Save RyotaBannai/0e94870c7eb9e7305e446fc08a9090d6 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib\n",
"import matplotlib.pyplot as plt\n",
"import pandas as pd\n",
"from statsmodels.formula.api import ols\n",
"plt.style.use('ggplot')"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"114674"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"path = '/Users/ryota/data/Kaggle/falldeteciton.csv'\n",
"df = pd.read_csv(path); df.size"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>dfb_Intercept</th>\n",
" <th>dfb_TIME</th>\n",
" <th>cooks_d</th>\n",
" <th>dffits</th>\n",
" <th>dffits_internal</th>\n",
" <th>hat_diag</th>\n",
" <th>standard_resid</th>\n",
" <th>student_resid</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>0.030580</td>\n",
" <td>-0.024454</td>\n",
" <td>0.000498</td>\n",
" <td>0.031543</td>\n",
" <td>0.031553</td>\n",
" <td>0.002506</td>\n",
" <td>0.629467</td>\n",
" <td>0.629277</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>0.017298</td>\n",
" <td>-0.014205</td>\n",
" <td>0.000156</td>\n",
" <td>0.017648</td>\n",
" <td>0.017655</td>\n",
" <td>0.002840</td>\n",
" <td>0.330821</td>\n",
" <td>0.330674</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" dfb_Intercept dfb_TIME cooks_d dffits dffits_internal hat_diag \\\n",
"0 0.030580 -0.024454 0.000498 0.031543 0.031553 0.002506 \n",
"1 0.017298 -0.014205 0.000156 0.017648 0.017655 0.002840 \n",
"\n",
" standard_resid student_resid \n",
"0 0.629467 0.629277 \n",
"1 0.330821 0.330674 "
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"N = 1000\n",
"p = 1\n",
"model = ols('HR ~ TIME', data = df[:N], hasconst=True) # allow constant, default is True \n",
"m = model.fit()\n",
"infl = m.get_influence()\n",
"sm_fr = infl.summary_frame(); sm_fr[:2]"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"t = 3*(p+1)/N\n",
"mask=sm_fr.cooks_d > t"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, 'HR')"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfsAAAF6CAYAAAAeZ/GvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzs3XmMXNl96Pfv3aqqq7qqt+qd7IXd3JfhaFaNpaeRRpI9I+NJz0++thM4tqXnCRC/l4cYQWwDQRAEAZ4fkDzHfwTGE6wEUhxHvk+OLT1HeF60y5oZaUazcR2STTZ737ura7vryR+32OwmOWQ3p9nr7wMQrHvr1K1zqpv81Tn3nN/RlFIIIYQQYu/St7sCQgghhHi0JNgLIYQQe5wEeyGEEGKPk2AvhBBC7HES7IUQQog9ToK9EEIIscdJsBdCCCH2OAn2QgghxB4nwV4IIYTY4yTYCyGEEHucud0V2ESS91cIIcR+o62n0F4K9oyPj293FTZVPp9ndnZ2u6ux5aTd+8t+bTfs37ZLuzdHV1fXusvKML4QQgixx0mwF0IIIfY4CfZCCCHEHifBXgghhNjjJNgLIYQQe5wEeyGEEGKPk2AvhBBC7HFbts7etu3/BvgXxMlv3gV+C+gEvgY0Az8Dft1xHM+27STwVeAJYA74FcdxbmxVXYUQQoi9ZEt69rZtdwP/NfCk4zinAAP4VeDfAn/kOM5hYAH4Yu0lXwQWHMcZBP6oVk4IIYQQD2Erh/FNoM62bRNIAxPAJ4Cv157/CvC52uPP1o6pPf+CbdvrSgkohBBCiLW2JNg7jjMG/C/ATeIgvwS8ASw6jhPUio0C3bXH3cBI7bVBrXzLVtRVCCGE2Gu25J69bdtNxL31fmAR+A/Ai/coemszm3v14u/a6Ma27ZeBlwEcxyGfz29KfXcK0zT3XJvWQ9q9v+zXdsP+bft+bLfnwfz89rV7qybofRK47jjODIBt2/8v8BzQaNu2Weu9HwBu7WQzChwERmvD/g3A/J0XdRznS8CXaodqr22sIJtF7C/S7v1nv7Z9v7V7fl7j3LkEmUyOM2dmsazNue5GNsLZqmB/E3jWtu00UAFeAF4Hvgt8nnhG/m8A36iV/2bt+JXa899xHEe2sBVCCLFrKAXXrhkMDZmk04pnn1X4/vbUZavu2b9GPNHuZ8TL7nTiHvnvAb9r2/ZV4nvyX6695MtAS+387wK/vxX1FEIIsX9os7MkfvADkn/3d1ivvIJWKm3atSsVeP31BENDJl1dIc8+69HQsGmX3zBNqT3TYVayn/3eIO3eX/Zru2H/tn0ntFubnSX5wx+iMpn4RBSheR7VT38aUqkPdO3paZ3z5y2UguPHfTo7I+CR7We/rpVqkkFPCCHEvmNduHA70APoOso0MS9ffuhrhiFcvGjy1ltWbdjeWwn0223LMugJIYQQO4XmunefNE30hxzKLxY13n3XYnlZo68vZHAwQN9B3WkJ9kIIIfadqL4efWmJNRHZdYl6ezd8rbExg0uXTAwDPvQhn3x+Z/TmV9tB3zuEEEKIrRGcOYPmeRDU8rq5LiQSBAMD676G78M771icP2/S2Bjx7LPujgz0ID17IYQQ+5DKZKj+/M9jXrqEViqhenvjQG+uLywuLsbD9tWqxuHDAX19IdoOTuouwV4IIcT+lEwSPPbYhl6iFNy4YXD1qkkqpXjqKY/Gxp2/qk2CvRBCCLEO1SqcP28xN6fT0RFx/Li/adnwHjUJ9kIIIcQDzM7qnDtnEYZw8mRAd3e43VXaEAn2QgghxPuIIrh61eTGDYNsVnH6tE99/c4ftr+TBHshhBAbpo+MYF26hFYuozIZ/JMniTo7t7tam6pcjifhLS1pHDwYcuRIgGFsd60ejiy9E0IIsSHazAyJ118HpVB1dRBFJF59FW1xcburtmkmJnRefTVBuaxx9qzP8eO7N9CDBHshhBAbZF26tDbVLKDq6jAvXtymGm2eIIBz50zefdcim1U8+6xLW9vOXDu/ETKML4QQYmPutU+rrqNt1/6tm6RQ0HjnHYtyWePQoYCBgZ29dn4jpGcvhBBiQ6KWlrsDvusStbdvT4U2wfCwwU9+kiCKNJ580mNwcO8EepCevRBCiA0KTpzAmJpCq1RQqRRatYrKZgkGB7e7ahvmefHa+ZkZndbWiJMnfRKJ7a7V5pNgL4QQYmMsC/eTn0QfH0efmSFqb49n4u+yrvD8vMa5cwk8D44dC+jp2V1r5zdCgr0QQoiN03WiAweIDhzY7ppsmFJw7ZrB0JBJJqN45hmfbHb3rZ3fCAn2Qggh9o1KBc6dS7CwoNHVFXLsWLDevW92tX3QRCGEELtKFKHfvIkxPo6qqyM8dixez/8BTU3pXLhgoRScPu3T2bn7l9StlwR7IYTYCZRCn5pCn5ggamoi6ukBfR8umFKKxI9+hL64iEqlYGkJa3iY6kc/imppeahLhiG8957JyIhBQ0Oc8jad3tvD9neSYC+EENstiuIANz8fJ6cZGUFduoT7iU+wJ6eG34c+ORl/Dul0fMIwiNJprLffxvvEJzZ8vWIxTnm7vKzR1xcyOBjsy+9Q+7DJQgixs+jXr8c92XQaNC3u0SqF9c472121LXdr6P5OerG44WuNjRm89loCz9P40Id8jhzZn4EepGcvhBDbzpiaigP8mpMG+sLC9lTog6hWsd58E31pCWWaBIcPE/X2rvvlYWMj+tgY3PF53PX53Ifvw8WLFpOTOi0tEadO+SST6375niTBXgghtlsyCYXCXffo1QeZJl4uY50/j7a8jKqvJzh58q589psuDEl++9tgGHH63Cgi8dZb+FFE2N+/rktEfX1w+XK8t2zt89DKZbyzZ9f1+sXFOOWt62ocPhzQ17eDMuGp7ZsnsE8HNIQQYucIjh1Dq1TWnNPKZYIjRx7ugtUqqX/4B/S5ObQgQJ+fJ/n3f49WKm1Cbd+ffuPGmiANtQ1yrlxZ/0UMA/cTnyBqbo531TNNvCefjL8E3IdSMDRk8NOfJtA0eOopj/7+HRDolcJ8+21Sf/M3aF//Oonvf/+R/xzuRXr2QgixzVQmg/uRj2C98w56sYiqq8M7e5aou/uhrmddvIiyrNtBV9dRqRTmu+/iP/vsJtZ8LWNx8Z4TCrVqdWMXSqXwn3563cWrVTh3zmJ+XqejI+L4cR/L2thbPirmW29hjI6ikkmoq0NbXCTxve/hvvjilq62kGAvhBA7gGptxXvhhU25llYocNfm67qOXi5vyvXfT9jdjTEycnsmfU1UX//I3nN2VufcOYswhJMnA7q7d1DK2yiKV1asnm+g62jVKvr4+JZmH5RhfCGE2GNULhdvzL5aFD3SoAvEOfJbWm7fkggCtFIJf5332zf0XhFcvmzys59ZJJOKZ5/1dlagh7iS0d2Je5RloW/xUP6W9Oxt2z4K/MWqU4eA/wH4au18H3ADsB3HWbBtWwP+GHgJKAO/6TjOz7airkIIsdv5x4+TGh1FAZgmhCGa5xGcOvVo31jT8D7yEfSxsXjouq4unnewCdnvViuX40l4hYLGwYMhR4/u0CV1pklUX48Whms2CdJcl/Dgwa2tyla8ieM4l4GzALZtG8AY8FfA7wPfdhznD23b/v3a8e8BLwKHa3+eAf6k9rcQQogHSaWofupTWBcvohUKqEwmno2/yUH3njTtkW6QMz6uc+mShabB2bM+bW07O+Wt/8QTJH/wA5Sux1+6ikWCI0fuutXxqG3Hd6EXgGuO4wwDnwW+Ujv/FeBztcefBb7qOI5yHOdVoNG27c6tr6oQQuxSqRT+44/jfexj+E8+uTWB/hEKAjh3zuTcOYtsVvHss+6OD/QAqqmJ6ksvEZw4gertxf3Upx79CMs9bEew/1Xg/6k9bnccZwKg9ndb7Xw3MLLqNaO1c0IIIfaZQkHj1VcTTEwYHDoU8OST3mbfGXi0LItwYADOnkVls9tShS2djW/bdgL4p8AfPKDovVZG3pWNwLbtl4GXARzHIZ/Pf+A67iSmae65Nq2HtHt/2a/thv3b9o20e2gILl3SqK+Hj35U8ZB74ewI2/nz3uqldy8CP3McZ6p2PGXbdqfjOBO1Yfrp2vlRYPXshQPA+J0XcxznS8CXaodqdnb2EVV7e+TzefZam9ZD2r2/7Nd2w/5t+3ra7Xlw/rzFzIxOa2vEiRM+SsFu/rg2++fd1dW17rJbHex/jdtD+ADfBH4D+MPa399Ydf5f2rb9NeKJeUu3hvuFEELsbfPz8dp5z4NjxwJ6enbYkrpdaMuCvW3baeBTwH+56vQfAo5t218EbgK/XDv/LeJld1eJl9791lbVUwghxPZQCq5dMxkaMshkFM8845PN7q995x+VLQv2juOUgZY7zs0Rz86/s6wCfmeLqiaEEGKbVSpw7lyChQWNrq6Q48eDu5IAiocn6XKFEEJsq6kpnQsXLJSC06d9Ojt3/pK63UaCvRBCiG0RhnHK29FRg4YGxenTPum0DNs/ChLshRBCbLliMU55Wyxq9PWFDA7u0JS3e4QEeyGEEFtqeBheey2BacITT/i0tMiw/aMmwV4IIcSW8H24eNGiXNZobIw4dconmdzuWu0PEuyFEEI8couL8bC962o89ZSiocFfvRGceMTkDokQQohHRikYGjL46U8TaBo89ZTH4CAS6LeY9OyFEEI8EtUqnDtnMT+v09ERcfy4j2Vtd632Jwn2QgghNt3sbJzyNgzh5MmA7m5JebudJNgLIYTYNFEEV66YDA8bZLOKM2d8MhlZO7/dJNgLIYTYFOVyPAmvUNDo6Qk5ckTWzu8UEuyFEEJ8YOPjOpcuWWganD3r09Yma+d3Egn2QgghHloQwKVLJuPjBk1NilOnPOrqtrtW4k4S7IUQQjyUQiEetq9UNAYGAg4dCmVJ3Q4lwV4IIcSGDQ8bXLlikkjAE094NDfLJLydTIK9EEKIdfO8eO387KxOW1vEiRM+icR210o8iAR7IYQQ6zI/H6+d9zw4diygp0fWzu8WEuyFEELcl1Jw7ZrJ0JBBJqN45hmfbFaG7XcTCfZCCCHeV6UC584lWFjQ6O4OOXYswDC2u1ZioyTYCyGEuKepKZ0LFyyUgtOnfTo7Ze38biXBXgghxBphCJcvm4yOGjQ0xClv6+pk2H43k2AvhBBiRbEYr50vFjX6+0MGBiTl7V4gwV4IIQQAo6MGly6ZWBY88YRPS4sM2+8VEuyFEGKf8324cMFiakqnpSXi1CmfZHK7ayU2kwR7IYTYxxYWNN5918J1NY4cCejtlZS3e5EEeyGE2IeUguvXDa5dM0mlFE8/7dHQIJPw9ioJ9kIIsc9Uq3HK2/l5nWzDMgMDZerrs4AsoN+rJNgLIcQm0SoVCENUff12V+V9zczonD9vUan6kDhHZC1z9SZcG9EY7O2jtbllu6soHgEJ9kII8UG5Lokf/xh9cRENiOrq8J55BtXUtN01WxFFcOWKyfCwQTarSOUuYCU8NO32TLyrwzdoamjElBR5e86WBXvbthuBPwVOAQr4AnAZ+AugD7gB2I7jLNi2rQF/DLwElIHfdBznZ1tVVyGE2IjEK6+guS4qk+HWXe/EP/4j7ksvsRMWqZfL8dr5QkGjpyfk0IDL6+8W1wR6AE2D+cV52lpat6mm4lHZyt/CPwb+k+M4x4DHgIvA7wPfdhznMPDt2jHAi8Dh2p+XgT/ZwnoKIcT6uS76wsJdQV0LAvTp6W2q1G3j4zqvvpqgUtE4e9aPc9vrwPvMuNe17f9yIjbflvxUbdvOAf8E+DKA4zie4ziLwGeBr9SKfQX4XO3xZ4GvOo6jHMd5FWi0bbtzK+oqhBAbou4zgz3avqQ0QQDnzpmcO2eRzUZ8+MMubW1xfQzDoCGbI7yzfppGc+POufUgNs9WDeMfAmaA/9O27ceAN4B/DbQ7jjMB4DjOhG3bbbXy3cDIqteP1s5NbFF9hRBifVIpolwOLQhYs0Bd14na27elSoVCPGxfqWgMDAQcOnT32vmj/QNcunqFQqmIAlKJBCcHj6DvgNsOYvNtVbA3gQ8B/8pxnNds2/5jbg/Z38u9Bpju+vps2/bLxMP8OI5DPp/fjLruGKZp7rk2rYe0e3/ZE+3+zGfQvvMdKJfjY8tCffzj1D8g2D+Ktg8NwaVLGvX18NGPKlruM7m+o70dz/OIVEQqmdrUetzPnviZP4TtbPdWBftRYNRxnNdqx18nDvZTtm131nr1ncD0qvIHV73+ADB+50Udx/kS8KXaoZqdnX0kld8u+Xyevdam9ZB27y97pt0f/jDawgJaGBK1tMT38B/Qrs1su+fFa+dnZ3Xa2iJOnPBR6oFVWFFcLm5KPdZjz/zMN2iz293V1bXuslsyXuM4ziQwYtv20dqpF4ALwDeB36id+w3gG7XH3wT+C9u2Ndu2nwWWbg33CyHEjqRpqOZmotbWLZ+BPz+v8+qrSebndY4fDzh71ieR2NIqiB1uK9fZ/yvg/7ZtOwEMAb9F/GXDsW37i8BN4JdrZb9FvOzuKvHSu9/awnoKIbaDUuhjYxjj46j6eoKBAWQ3lvuLIhgaMhkaMshkFM8845PNSspbcbctC/aO47wFPHmPp164R1kF/M4jr5QQYmdQisQPf4i+sICqq4OZGcyrV3Gffx6Vy2137XakSiXewGZxUaO7O4yX1EkuHPE+ZNqlEGLb6ePj6PPzcaCHeIJbMon11lvbW7EdamoqXjtfLGqcOeNz8qQEenF/ki5XCLHtjPFxVDq99qSmoRe3btLYbhCGcPmyyeioQUOD4swZn7o6GbYXDybBXgix7VR9PUxNceesMiWzzFYUi/Ha+WJRo78/ZGAg2AmZeMUuIb8qQohtFwwOoim1JuOcVqngHz68jbXaOUZHDV59NYHvazzxhM/hwxLoxcbIr4sQYvtZFu4nPkFUXx9notN1vMcfJ+rt3e6abSvfh7fftrhwwaSpKeLZZ11aWrYvBa/YvWQYXwixI6hMBv8jH8Hf7opskFKKMAzRdX1TU80uLMSz7V1X48iRgN7eu1PeCrFeEuyFEOIhzS7McX10FN/3MQyDztY2erq6P9A1lYLr1w2uXTNJpRRPP+3R0CCT8MQHI8FeCCEeQrlS5vL1IeqSKYxa8p/x6Uks06Sz7eE2wKlW45S38/M6nZ0Rx4/7mPK/tNgE8mskhBAPYXRqilRibYa/hJVgcm7moYL9zIzO+fMWYQgnT/p0d8u9ebF5JNgLIcRDCMMA7R430aNwY0PuUQRXrpgMDxtks/Ha+UxGhu3F5pLZ+EII8RBaW/JUPXfNuTAMadxAet9SCX7ykwTDwwY9PSHPPONJoBePhPTshRDiIbQ0NJJvbGZ2YQ7DMAnDkHQqRV/3gXW9fnxcZ2xMo1LROHvWp61Nhu3FoyPBXgghHoKmaRzpP8SBjk4Wlwtk6tI0ZLMPfF0QwKVLJuPjBv39cPq0Syq1BRUW+5oEeyHE1lMKrVAAy7p/uSiKyyWTtzfJeZBKBc11Udks990dRin0mzfRKhXCwUFuTXsPw5Cq62IYBkEYUHVdNE0jU5dG13V83yeVTBJG0crjfFMTS8vLTM7OkLAsUskk6dTt+kZRRMWtUi6aXL5cT6WiMTAQ8PTTirm5W9VRVKpVNE2jTqK/2GQS7IUQW0qfmsJ6/XU014UgQOvthbNn79q7Xh8dJfHWW1DbDCfq6sL78Iff/wuC72O98grGrehpmvgnTxIeOnRXUW18nNy/+TcYMzNESkE6TfELX+DG4UFGJydYKhQYnhpjabFAGIUYukEykSTf3ER3WweziwskE0kas1nGpiaYmJmmWCpTKpdIJVP0dh/gxOBhTh4+iuu5vDc8zPiYxdhIlkza5Bd/oZ22VmMlSU6xVOLi0FU830MDUskkJwaPkrrjM8HzQNeR9Xhio+Q3RgixdXyfxKuvQhRhXr2KVq2iDQ+Tfvttyl/84u2eeLVK4gc/IPHaaxjz8wBE2SyqWsX/9KfXXjOKsN58E+vHP0afmYF0mmBgAJXJYL31FlFra9zLX1U+++/+HVqlQtTYGJ9TCuvf/3tGXv4XmE1NzC7OUyyWWC4VSdfVEYQBlUIFTdeoeh7pVIpiucTw2AgTM9PMLS2CitA1nUK5SNWLRwMipai6iomRVqZnIVAT5NLjfP8Ni8eOHqe5uRmlFBeuvYdpmNQlU7XqKC5eu8LjJ04BoC0tYf3kJ/EugLpO2N6O/9RT9x+5EGIVmY0vhNgy+vg4KoqwLlwAQKVSkEphjI1h/fSnK+XM4WESP/whermMymRQmQxaEJD52tfWbJYDYL77LvrERBwI6+tB17EuXYIgQKXTGFevrq3DxATG5OTa3rGmUXZd8hcusLRcoOq6uJ6Hoet4vocXBCvD9gtLC+i6xuLSEvNLiyyXixBFuL5PqCI0TadUKTM2NcE7FyY593YTS4smRvIKPX0zpNMWvu8zv7TEpaFrFIpFfD9YU0dN0yhXK1RdF8KQxA9/iBZF8WdRV4c+N4f1+uub/NMRe5kEeyHE1lEKfX4epe5YXmYYGOPjK4Fcm5/HKBTW9lx1HSoV9JGRNS81x8bu3hrXMNAnJ+Nr3fleSvHQi9tU/Md1PcIoRKmIKIjiIA+oSIGCUGnMTrdw6WIThhFxoH+YltblNZeyTJPZ+VmUuvcs/Fsr+PXxcbQwZE1ifMuKv7CE4cO2ROwzEuyFEFsm6u5G833W7M8aRUTZbHyuFpjDzs67A1kUQV1dfK9/tSDuFaumppXH6Dqa76OVywSDg2sv09VF1Np6uyyAUqQTCWaPHaOhPksykSRpWYRRRMJKkDBNdF3Hsiyaco34YUCoIhpyOUzLQCnQDQNQBIHJ4tQgXuUAnV0Bp84sYiU8NE2vNSMiXVcXD/NHivp0BuuOe/BKKVLJJKlkEs3zUPfaYOeOLYGFuB8J9kKIrWNZVH/+59EqlfhPuYwyDILBQaKGhpWefNTTQ9DTEwd810XzPDAMwu5uwv7+NZcMm5shigj6+uLh/moVbWmJqLER/8wZ1J1JbnSd5d/9XVQyibawgDY/H88F+MIX6D56lCAKyTc1UZ/OkElnADANg8ZclsaGLAM9vWRSGdLJFM25RnKZHOm6FDoapeUG5iaOYWg5Thyv8uILrRwb6CdTV0e5WsH1PUzDoLutA6UU2fp6TNPk+ODheF6AW6XiVlFKcXzwSNy+7u5429/VlIq/ID1oNYMQNdpdw2m7lxofH9/uOmyqfD7P7Ozsdldjy0m79z7jyhUSr72GSiTItrRQqFbxnn8elcncLnP5Muk/+7N4dxjDQLMsyv/8nxOcPbv2YtUqye9/H61SiYfvy2WC3l78J5+8/6x1peL7+dUq4ZEjK6sBgjCkUq1iGgZ+4Mf3zdHIZjLouo7n+9SlUkzOTPPejWuoSHFjbJJLF3WqpVYaGkIePxsx0N/B6SPHMAyDKIq4evMG41NT1CWThCrCMkyef+4jlGqrDZRSlCsV0CBTl177eV27hvXOOyjDiO/dWxbexz6Gqq/flJ/HVttPv+urbXa7u7q64PYdn/uSYL+DyT+I/WXftdt10cfHae7qYtay1g7trypjvvMOhCHBY4/B+621Vwp9ehqKRaKurvcvt8mCIGBkfImLFzOUijrtHQX6+n0ac/Xk6u9OsOP5PvNLiyQti4Zsjra2tvX/zGufF8kkUUfHvT+vXWLf/a7XbGewl6V3QojtkUwS9fdDPg/v9x9gMknw1FMPvpamEbW3Q/vDbS37sCYnk1x7r4tMCp550qelJXHf8gnLoiPf+nBvduvzEuIhSLAXQogN8n24cMFiakqnpSXi1Cn/zpxAQuwoEuyFEGIDFhY03n3XwnU1jhwJ6O0NucdOt0LsKBLshRDrpxT6yAjG2BikUgTHjq0/Z/0OV65UGJ2aJAxD2vN5mhsa1zyvFFy/bnDtmkkqpXj6aY+Ghj0z50nscRLshRDrlvjHf0Sfm4sD/OIiyRs3cD/6UVQ+v91VW2N6bo7x6UmCMCCTzjBwsJfEfZapzczP8d6NIVKJJJqmcenaVdpaWhjsje+RV6vw7rsWCws6nZ0Rx4/7kp5e7Cq7dzqnEGJL6TMz6LOzt3vyhhHnn3/77e2t2B0mZ2a4dvMGURTnqi+VSrx98TzR+ySgUUpxY2yEumQKrTYen0ommZqbpeq6zMzovPpqkkJB5+RJn9OnJdCL3WfLfmVt274BLAMhEDiO86Rt283AXwB9wA3AdhxnwbZtDfhj4CWgDPym4zg/26q6CiHupo+Px7ns7zy/vHyP0ttnbGqC5Kr0ubqu4wc+03MzdLTePVvfDwJ838dIrt1URsfgjTddSoUc2azizBmfTEaG7cXutNU9+487jnPWcZwna8e/D3zbcZzDwLdrxwAvAodrf14G/mSL6ymEuEPU2Hh3qlpA7bBp6P6d2eYA0zBZLlfuWd40DIw7do+rVgzOn2tkdrqenp6QZ57xJNCLXW27h/E/C3yl9vgrwOdWnf+q4zjKcZxXgUbbtju3o4JCiFh08GAc2FflrNcqFYIjR7axVndL3WP0wfM88o2N9ygd9/zb8614vg/AzHSKt99qQNfSPPdhk2PHgt2cv0YIYGuDvQL+zrbtN2zbfrl2rt1xnAmA2t9ttfPdwOqtrUZr54QQ20XXcT/+caKWlnjnOF3HO3uWcGBgu2u2xsDBHqquS1i7R+96HtlslsZcw/u+pq/7IJ2tXVy7mmPoWpbOjhS/8kuttLbKRjNib9jKaSY/5zjOuG3bbcDf27Z96T5l77Vq9a4xtNqXhpcBHMchv8NmBH9QpmnuuTath7R7h+ve3O/dm93ufD5PV2cXN0Zv4vo+rc0tdORbVybf3cviIhTmW+npgBc+qjh8mC1ZO79rfuabTNq9De+9VW/kOM547e9p27b/CngamLJtu9NxnInaMP10rfgocHDVyw8AdyW+dxznS8CXaodqr+ValvzR+8t+b7dSisXCEq7v09zQeN+lcuvRkrs9bD83N3fPMkrB8LDB1asmiYTi9GmfpibF+xTfdPv9Z77fPKLc+OuyJcHetu0MoDuOs1x7/GngfwK+CfwG8Ie1v79Re8k3gX9p2/ZqgEbaAAAgAElEQVTXgGeApVvD/UKIvcfzPd65dBHP9zF0naGbNzjY2c3BzvX/Z7ZRrgvnz1vMzuq0tUWcPOmv2TE2DEPK1SqWaZLaYZMQhdiorerZtwN/Zdv2rff8c8dx/pNt2z8FHNu2vwjcBH65Vv5bxMvurhIvvfutLaqnEGIbvDc0hKZpK0HVwuLmxBj5pmbq7jXhzveZmJkiDCM6W9vuWeZ+5uZ0zp2z8H04fjzg4MFwzfPj01PcHB+t3ffXyNXXc2Lg8JpZ++YPfkD6L/8SY2ICVakQDQzgHz2K//TTBGfPYrz1FumvfQ19fh6Vy1H+/OcJnnsONA19ZATz8mW0dBoriiAISLz5Jlq1StTeTuUzn0Ft8aY+Ym+TLW53MBnq2l/2c7v/5h/+loS1dse4MIpobW6hr/vAmvPzi4tcun4VyzDRNA3X9+jrPkh3e8cD3yuK4No1k+vXDTKZeO18Nrv2/8BypcKbF86t+QIRhCFNuQYO98UZ9awf/YjMn/wJuu+jzc6i+T6EIe6HP4xqbyfo7CT1938PqVS8FW0UQbFI6Xd+h/DAARJvvolKp8nlcpR/9COst98mOH4cTDPeZcc0Kf32b6Oamx/2Y93R9vPv+nZtcSsLSoQQ207T7v6vKIoirDvWvyuluDZyg1QiiWEY6LpOXTLF8PgowT3W169WqWi8/nqC69cNurtDnn3WuyvQA0zMTq9JygPxWvyFwtLKcfI//kdIJNBmZ9ELBbQgQNN1zGvX0FyX5N/+bVzw1po9XYd0muS3voV18SIqnY7P+z7m9euQTKIv1a5vWWjVKtabb963PUJshAR7IcTGBAFapRLPbnsIYRjieh6rRxXzTU0E4dpgHUUR7a1ta865nofn3x3UNTQWi7cz+QVBgOd7eL5HEIZMTur8449N5hcCjh4rcfy4h1IBnu/jeh5VLy57Pyuz86MIfWkJY3oavVKBKELzPCiX0ctloJZV8M7PxzTjLwaVVcl9PC/+Y5rxJIJVdlpmQrG7SYZnIcT6KIX15psYIyMQRahkEv/MGaIDtWF2z8M8fx59YQF1a0e85ua497u0RNjayrXFBWbm54iiEMu06DtwkHw+T/+BHsJIMb+wQKQikokEJwePYN7RszdN855L4iKlSJoWQRhy6doVpuZmmZiZRkU6QXWAmakUqTqXTOM1pt/w4GcaCdPE830Wl5dJ16VozDbS29VFb/cBJmemqUuuGsYPAppvJeXRdTDNeC3wrSCtaWhhiKolzQ87O1cC/4pqlbC/P85E6Hnxt4e6OlQmg7awcNdmQmHHg29LCLFeEuyFEOtinj+PPja2ZkvbxE9+gtvSgrIskt/5Tnxv2jTRlpfjY10HXUcZBkszMwR1SRJnzgDxtPf3blyn92APuq5zpK+fqKeXsPZF4J51MAyaG5pYWi5g1QJrFEWkU0my9fWcv3KZYqXC1OwsKshw80ae2fkCnZ0TaKkxctlWpudKlKsVPN8nYVlYpkWxVCaTSjM9P0sQhvR1H2RkYpwwCtGWizQrxaFDgyv1cD/6Uczr14kSCXTXRVUqkEgQtbQQ5XJ4zz9P+s//HK1YhHQaKhVUIkHl859H5fMkv/tdVCIBuk5w7BiJN96Ij5VCLxYJDxzAP3v20f0wxb4jwV4IsS7G6CjcsQRN1dVhXL6MymbjiWW3ntc09NlZtEKBoBa0FqOQhrk5wukZqm2tACQti5vjo+QbmoA4da3+gNy0R/r6uT46wtziAlEU0ZDNcbi3jyAMWSouU6lWWJjLMjOVxzAimtuGSNX7BGFIEMVlCsvLRCoiDEIaGxrI1NURhCHLxRItjU0YhsHTvf2oP/sz0uPjJAC+/V3c557De/FFot5eyl/4AtaPf4wxOUlkWdDSgt/fT/XXfg3V0sLS2bOk/vqvMa9cIeztpfrP/tnKhLvqL/wC5uXLqESC6ksvUfmlXyLx4x+jLy7iPv88/oc+BB8wz4AQq0mwF0KsixZFqFog1mdnMSYn4/vV1SphZ+ddXwT0QgFWTZpTKIJUHemJiZVgr2kaQbB22duD6LrOQE8vAz29a877gU/gw5X3mpgcC8g1VjlwcI7xmRJRpKPQ8D2PcqUcf6GIINIjoiiiVKnQEkVESqFXqhg/+CG5n/wEc3wccjmC7m4wDFI/+hHhkSMEJ05gTEzgffrT8XC8UmiVCu4nP4mqr48rlM1S/fVfv3cjUimCxx6DfJ6oNjvb/dzn7l1WiE0gE/SEEOsS5vMQBBijo5jXrsUT9UolqFbRx8agWl37Ak1bsyNeJlVHFIZEqzaDd32P7s7N2eOquJzg2ntdhF4jrR2z9PbPYFoRCStBQy6HoetUPRfLsNB1HcM0Vpb7+X5AMmGRTiXhwnnaUimM+fl4CN7zMMfGAIjq60m88goql8N9/nlUJhPvE5DJ4H7847cDvRA7jPTshRDr4j/+OMnvfje+V51MxhPOurtRLS1oy8uoKIqH8i0Looiovh6Vy628vj3fSrC4yPXOdpTvEUWKrrZ2mnINH2jtsVIwNGQwNGTS09WFnjpH2dOYXagCikM9vRRLJeoz9UxMTZFIWuihgWUYaMByqUwykSBUirxucriunsSqLyToOlQqt3f7q41uqKYmvI985KHrLcRWkmAvhFgfy8L7yEfQJydRpolqaFi5r6xME//MGfRSCX1uDpVK4T333JpRAD2VouOzn6OuqYmKV6Uxm7srkc5GVavw7rsWCws6nZ0Rx4+Drp+gUIx3vgNFXSpFNlPPcqnE7Pw8F4feozHXQBRG+EGAHwZkM/Uc6eunpVgivVxCAUF3N9bwMNQmzqEUeqlE+bnnPvhnKcQWk2AvhFg3lUqhmprWDM9D7X5+aytBf/+a88GJEwTHjsU9/kQCNI0skOWDD3fPzOicP28RhnDypE93963taHUaV40o3NKQzdKQzVKfSXN9dIRQCzFNk85cjmOHBuP7+LkGqCWzCZ56CqNYRJ+cjK/qulReeIGor+8D112IrSbBXojd7NZ98kQCrVqNl2+Zj/Cfta7jHz2Kdf58fL8a0CqVeDh/1f3qKIpWlrbpun7X5L31CsMQPwhIWFacGtfz0HWToWtJbt40yOXineoymfUn+GnPt9La3ELFjTe5WTO6YJp4TzxB4o03QCm8n/s5lO8TnDpFeOjQQ7dDiO0mwV6IXUgrl7FeeQW9UECfmUGfnyfs60Ol0wQ9PfFM70e0IXt45Agql4uH58MQf2CA6NChlefHpiYZnZwgCEMMXaeztY3eO/LbP4hSiuujN5menSNUEVW3QhhFGGS5ca0ZQ2vkmSebOXo04gEr9e5J13Uydel7PhcdOEC1owN9agosiyif56HeRIgdRIK9ELuQ9aMfoUXxsLUxMYFKpdCHhwkefxzz5k1IpwmOHHlk7x91dODdI8PbYmGJ4bFRUsnkStKb8elJUskk7fnWlXJ+4DM6MUHZrdJTKVOXSK7Jljc+Ncn03CyJRIIgDBgenWFhvh7l9ZNKmnT1jpGsX0LXe++qw6YwTaLu7kdzbSG2gXxdFWKX0ZaX0Usl0HWM8fH4/rmmoVeraJUKKpXCuHlzW+o2Pj111yYyyUSSqbnbs+2rrssb595ldmEe13UZnRznrQvnCMLb6+2n5+dWhten5xaZGOtgaqKDIJrj1JkF8vmI+cXFrWmUEHuABHshdpswvL3JShStHa5ffX4bVKou80uLFMslFLfvo6/e9Ob62AiWaa7sDW+ZFlEUMTY1sVImiuLyxaLJ5QttFJYytHcs0D84TSIRty1S29NGIXYjGcYXYpdRDQ0rs+Gj9vY47WoyibKsOG+97xNuwxD0tZvDTM3NMrswj2kYmKZBf/dBIqXIN7WslKtWq2h3zCcwTZNiqbRy3JDLcf5ClYmxBurTVXLN50lnPJoa4gQ8YRSRq797xv1qQRhSLJVIJZOkZGKd2Ock2Aux22ga3jPPkHzlFZRlEeVy6LOzBEePopfLhC0tBKdPb2mVlpaXmZqbob2lhSgMWCou4/sR18dHODl4lAOr7u9bloXrumsCfhhFKwHZdWFpZoDpyQnSmUX6B5ZZXLYIAgPLtKi4LulUisHevvetz+jEBCOT4yu9/1x9lpODRx6Yd1+IvUqCvRC7kMrnqX7mM+jT0/Dcc0T19fHWstlsnOzmEYmi6J4Bc3puhqSVQNM0uto74qVt1SqGaXJycO1Ewb6uA7xz+SLJRFw+UoogCDjY2cXcnM65cxa+Dy9+soOm5gyVaoVc9hSaprG0vEzSSpC9T1raYqnE8MTomi1qK9Uq10dvMtDTt2mfhRC7iQR7IXYrXSda1WOOauveH4X5xUVujN6kUC6BUgwc7ONgV9fK86ZpEimFUeutW5aFZVlr7tXfUp/JcProMYbHxvB8j/p0mu7WDoZvpLl+3SCTUXzoQz7ZrALSZNK3l8jlm5ofWNeJmWlSibXD9qZhsFBYesjWC7H7SbAXQtxXuVLhwtXLzMzPU65WABgeG+PsiVM8duw4AN3tHUzMTGOsCrKe79PV1n7Pa2Yz9Zw6chSAdDrP975XYHFRo7s75NixgFWr8DYs3oRO3TUvAB5N3gEhdgO5gSWEuK/RqUkWC8u4vk/CSpCwEqTr6rhw5fLKpLqEleD4wGGUUlQ9l/mlJZZLReYWF7g+cpMwvPc2tpOTOj/4gUappHHmjM/Jkx8s0AN0tXXget6ac0EY0NLQ+MEuLMQuJsFeCHFfQeBTqpQx7rhXrxsGEzPTK8dNuQaeOHWG3u4D1CUTtDQ2oZRiZmGedy9fXDOkH4Zw/rzJO+9YZLPw7LMeHR2bs5QuXVfHQE8fnu9TcV2qnkuuPrfhLH5C7CUyjC+EuK98UxNBGJDQbyfLiVREuq7unpP1JqenqUvVrRybhkHFjXv7LY2NLC9rvPOORamkcehQyNNPK+bn15/bfj06Wltpa2mh6rpYlollWpt6fSF2Gwn2Qoj7am3O09HaxvjUFImERRiGWJZFY7aBzta2NWVvbYBzZxY9yzRZLhUpL7dw+bKJZcGTT/o0Nz9cbvv10HWddF3dgwsKsQ9IsBdC3JemaXz8med488J5xqYnSSWS5Orr6T9w8K5gqus6CevuXnS5EjByo5VKySSfjzh50pcN5ITYQhLshRAPpOs6T5w6zZngGL4fkEom3zdBTW/3AS5du0oYhnFvvpyitDjIwY4cR48G9PSEH2xDPqXQp6cxRkZQmQzBwADcMZIghFhLgr0QYt0s03rg/e+mXANe4HF9ZAS3dIDlxXYaG1xOnS7Q0f7Bu/PWa6/FO/2l0zA1hXnlCu7HPvZIkwkJsdt9oLtltm2fsW37P2xWZYQQu9vN8TH+/G/+mp+9+x5jwwcYH83Q0uLx+IeWmVkc+sDX12dmMCYnUZlMvKDeslDJJNabb25C7YXYux7Ys7dtOw38AXAWuAL8j0Ae+F+BTwFfeYT1E0Ksw8TMNBPTUwRhWFt61rsmXeyDLJeKXB8diWevmybpVIpSpUIQhpSrFZKJBIZukEmnGezpZWm5yNj0BL4fkEwksEyLoZFhLl+/RmGxjsWZPkIVMDA4TrrBB62fcjW453uXymXOvXeZUqXM9NwMheUiuqGTy2Q5PjhI/4GelQQ5+shIvNnPKsbICMbERLwvQGsr/tmzcI95A0LsZ+sZxv/fgceBvwVeBE4Dx4iD/G87jjN7n9euYdu2AbwOjDmO84u2bfcDXwOagZ8Bv+44jmfbdhL4KvAEMAf8iuM4N9bdKiH2uCiKuDkxxmJhiYVCAc/zaG1uwTQMXNflnUsXePLUYyvbyN5Pxa3y7uVLpJJJTMNgfnGRtycnOHSwh6pbZXZ+Hk3XGOztp1qt8qPXXyOZrCOdSmEaBjfGRllYWiKKFHPT7cxOZdDNAoePLqOZVSDPcmmZptzdSW2CIOCn77yF67mMT09yc2KcMAhpbmhkWSvx5sXzhGHE4b5+IN7xj7Exbs3uM4eG0Obm4h3/TBN9Zobk97+P+8lPburnLcRut55h/J8HPu04zu8BLwEvAP+Z4zj//UYCfc2/Bi6uOv63wB85jnMYWAC+WDv/RWDBcZxB4I9q5YQQNefeu8TU7CxRpJiZn2NucYHZhQWAWi9YY3JVwpv7GZ2YWDODfm5pgUxdHXML8ywWlrAsCw2NhcISmqYxs7CI67oAhGFIqVTCdU0uX8hTWGiivbNCvvMKml6JN7mJQsJI3bVMD2BsehLDMPCDgMWlpXhZn2myXCph6Dqu6zI1O7OSgS/s7Y2H76MIogh9dhYtighv7RFgmmjLy2hzcx/g0xVi71lPsK93HGcawHGcUaDoOM4PN/pGtm0fAD4D/GntWAM+AXy9VuQrwOdqjz/L7dsDXwdeqJUXYt8rFJdZLpexzHhgLg6QFvNLCyji5DSmYVB2q+u6nuf7a2bWR7cCaxQS1B4bhoHnxQFe09RKOtogCllcqGfsRi+ep9M3MEu+bZqGbD26pqGUwjDine8675Env1qNbxsEgU8Qhmi1/PW3tqZVSuGHwUo9ME28j38cVV+PXqlAFBH09BCt2pRH6Tpasbi+D1OIfWI9w/imbdsfZ9UuEnceO47znXVc538D/jsgWztuARYdx7l1I28U6K497gZGatcObNteqpVfM5Jg2/bLwMu1cuTz+XVUY/cwTXPPtWk9pN33FkURCqj4Hs1NjVimRRRFNDY2EAYRvu9Tn6nHMAwqbpWjA4M0NzatuUYYRWiwJrj3e1VujI6QsOLla41NTXhevBtdueqCioN7V0cnufosdak6OtrbSafqmbiaZm4mTUNjlROnXVJ19SwWQgrFIvnmZrrbO3jp+RfWZNRb0yYNLl27Smu+lan5OfwgwDANTMMknU6TCBN0trXT1dl5e2ObfB56e0EptG98A8w7/hszDNSJE7ALEurI7/r+sp3tXk+wnwb+j1XHc3ccK+DQ/S5g2/YvAtOO47xh2/bztdP36qmrdTy3wnGcLwFfuvX87OxG7yrsbPl8nr3WpvWQdq8VRRFXh68zv7RIpBSmrjM5N0vVdQmCgDCMqFQr1GcyFEtFXM+jMdeAW6lyZe5qvG+8rvPe0DVK1Qq6ptPc0Mhgbx+6rpMyE3hVl0KhEG9yYyW5OTKCYRgsFYtMzkxRl0oxMzvL1Nw0pm7wxjvDlBYG6ek8RK5xEl97D99voVj28MKAMAyZmZmhXCrzl+UKz5z9ECiFZVnUp29vxaspRba+npGxUXKZeiampygtl+hqbWdhcZ6WxibyjU3Mvc+wvN7XR+KnP0WlUqBp6JUK/uHDBKUS1Dbp2cnkd31/2ex2d60a0XqQBwZ7x3H6Pkhlan4O+Ke2bb8EpIAccU+/0bZts9a7PwCM18qPAgeBUdu2TaABmN+Eegix67x34zqF5cJKz7tcrTA0cpO2phYsy8I0FH7o09XWQRhGBGHA0Mgw71y6SGtTM8lkgonpKfq6D67s875YWOK9G9c5dmgAXdc5c+wEswtzLBYKqCjCMAwmZ2aYmJnC0HTKlQo3x8dJWAmCajfFpWYy6SXMzP/HmROd9Hc/ycLSEn7g8971a2iahjIMlksl3jj/LpeHrvGhk6fRDZ10KsWpI8ewTAtN0zh7/CR1iSTzC/P0dR1guVKisFykvSVP/8EeUvdJtRcdOIDb3Ixx5QpaEOD196OaH7znvRD7zXqW3v2Qe/SqV3Mc55884Pk/IF6+R61n/986jvOf19bof554Rv5vAN+oveSbteNXas9/x3Gczd0pQ4hdIAxDFgqLJK3bGeLmFhfJNzZRn8kQRRG6rtHd3oGmxbPbTd1gqdZLn5ydoT2fx/U8puZmV/aXN02ThaXF2ut1dF2nraWVxlwj3/r+twmjCNdzSSaT6Fo8OS8KDBYWDuG7DTS1lGhqGcFKJHBdDz/w6WhtZXJ2huVSkbq6NPWJZHwd38PzPEqVMm0tecIw4r0b1zk5eASIJxTmaz34h6HSaYLHHvvgH7YQe9h6hvH/dNVjjXgp3n+1Se//e8DXbNv+n4E3gS/Xzn8Z+L9s275K3KP/1U16PyF2lUhFqGjt99woCuM173VpmlZljRseH6G36yDj01NYtdn1pmEwOz+PaRiUymUAlpaXWVou4Ic+hxf7ya/qCU/PzYIC16viBQG6ruN6LpVShuX5fsCkoeUm+XaXIAzx/ThYLywtUZ/OUK1WiZQiqtW56roYaARRSKlSAeL5AoXlZZRSt+/DCyEeqfUM469JmmPb9r+789xGOI7zPeB7tcdDwNP3KFMFfvlh30OIvcIyLVKp5JqxtVx9PculIvWZDIXlZZbLRQzdwKr1/nVdB6Xie9i6jqHrVF0X04SJ6WmWlguYVvxP/70bQ7i+R3d7x8prc/VZwlCha1AolSksdLA424SRqNLUeoVMWgFplIJkMkXSSqwE7YZcDqUUqUT8ZUPTIAI0Q6cpl1tpg6ZpEuiF2EKPaHNJIcRmOdJ7CC/w8XwfPwhIJVL0dB7g+thNRqYmWC6WKJZLVF2XquvS3NBI1fMolIosFQvUZzI0NzaiGzpzSwtoukYYhhxo7yCVTDI6OUEQBERRRFtLnrpUitbmFtyqzvz0UcqFdjK5eVo7LmOYFUzTIPADdF3jmTOPE0QB9Zk0QRhSl0xy6GAfSlN4QbykLwhCOvNtNNSCfRCGNEseeyG2lGyEI8QOV5/J8NSpx5hdmCcIAlpbWlhcLhCpiCAISCVTZNJ1BH7A6PQkdckUruexWFhC1zUmTJPHjp2krbmFNy+eJ51KkqvPoWkaQRhwfXSYqlslYSVoyGYZ7OtnclLDLeRIWkVaW0fINVXw/UYM3SSRMDF0g6fPPEZjrpHHjh0nmUpRqVRoaWziiVNneOP8OUbGR1CazmPHc6ST6TgRj6bRmM0x2Nu/3R+rEPvKeibofeLO1zzkOnshxEMyDIP2fOvK8fziIo3Z28Pinu+jAZ35NhaWFujMtzJwsJe6VBLTMIlURHNjE025HKlaznyF4sbYKFGkyNSl43vpxQpX3iuD/xinj1/l6LEKiUQ3mgamaVFYLvD4idPksll83yeZSNxes79qgt1zjz+Be/J0PKSfTKKUqt1KMB64a54QYvOtp2f/5TuON7zOXgixuVLJJMvFIkEYMDI5ge/7AAR+yEBPD+l0ek35qutiGgZNuUYKpSKWaVJ1XSqVCm0teXRdp1wyuHqlmcJyxCc/ptPpVrBMHU2Ll75FUURrcwvNjXGOe/MBefeTq/aY1zSNutT6N+YBwPcxz53DmJtDGQbB4cNEBw5s7BpC/P/s3VmsnGl62Pf/++21L6fOvvFwbbJ7epnxzEgjS1YkBZEcA75xKoABwwkE6Ma5CJwLG7lJLpObGAYSOBCgCzkI4BwoMaQY8aI1jmRJs/TMNMlubofk2bfat29/31x8xTMkm71Mq5vdJN9fo9GsOl9VfV8VTz/1Ls/zaMCn26Cn59s07StmeX6Bg+Njdg73AfCDgHHgU8oXuL+/yxuXrjxxvIHBaDJGKsl4MkYYBo5lM1tvMN+Y5fjIY2e7iGkqLl3pcv78HGl6mVv37xGEAQoo5QtcvXjp+VygUrh/8icQx1m9eylxvv994jAkvXDh+ZyDpr1E9Jq9pr2AbMvmwuo6D/Z2eLi3B0gatQZ5L8f2wR5Lc/NEcYRMJYV8njSV/MF/+FPynsdMrY6UkkqxhJQmd26X6XVcKtWI8xcHmKbKqu4JwTfeeJMgDDEMcVbU53kwDg8Rk8kT7WxVPo91544O9pr2Gehgr2lfYePJhIcHe0RRhOM4nFtaoTCdok9k1ikun8tlefR+NnKvlsvcvHObmXodQwgOTo8JooiFmVmG4zHdwYC1pSV2D3yG7YuctnqsrfWYmx8jFbx27tITaXEfV8Hui2J0Oqhn9KQXUZR1vDN0IpGm/TR0sNe0ryg/DPjx7ffP8tiDILv9ztU3yHkeJ+02SmTd6QajIVJJOv0+ruNwZeMCF9bWCaKQ+zs79AdDkmn7WFsYXL8eoeIVVpdN/u7fWSJM2hjGDLO1OuYnrMU/D+nCAta9e6inmtyoXE4Hek37DPRvjaY9B0mSsHNwwO37W5x0WkgpP/ExO/v7TxSsEULg2g7bB3tAtulurt6gNxhkz6fANAySOCEIA4Io5MHONsNJ1u41jELi2OTh1jwHB0XK1RE//3OKmbrF0tw8C43Zr0SgB1CNBrLRQEyr7iElYjwmfuONL/fENO0FpUf2mvYZpI/1ef8kfhjw41sfYJDVpO/0exyfnPL65StPtJp9WhiFH6oyJ4QgmvaStywT13OZqVRJlcRAoIQgCAMG4xHv3rxOEIaMJ2PiJKXbdfD78xhGyPq5Fm+9dY5atfjZ34QvWPSd72Bsb2MdHKBsm+S111CPVeHTNO3T08Fe034KcRJz+/4Ww/EYhaKcL3LlwoWPzR2/v7uDbZoMxkN6gwECwbhQZK7TfiJ3/mmFfJ5Ov4/52BeCVEoqpRIAS3MLbO/vsTi/QG84II4jDCFYX1jkh7fen27MSxiOA9onDQy5TGMmobHwkHqjxGsbX/GNboaB3Ngg2tAJQZr2V6WDvab9FG7evUOSJGc55EEUcuPuHd65+jqQFbvZPtgjiiMc22F9aQXf9zlqnWb57Wb2K3fSbnHrgfmxwX5taZl2t0uiFJZpkqQpKMXa0jKQpd9du3CZH35wk2qpTM5zWZydY//4CCklURSRxC7to0uEgUutfsTX3qww3zhPtVTisHXC+pLOW9e0V4EO9pr2KQVhyNifkHN/UhzGMAx8f4IfBCRpwq37d/FcD9uyUUpx6/5dpJQMhoOzRjWPHjcYDs9azD6Lbdl8/fWvsXd8yHjiU8/nWJlfxHps09rVi5cwLZNWt4ttmqRJtlmvXqkS+bP02rNY5pD5jQO83IDXzr99tvQwebQermnaS08He037lFKZPtF97hGlsgTq+BwAACAASURBVJ/tHh7gOk+mqbmOS7ffI0qSs2CvlELKlEqpQhhFH1tZzrIszi2vfuTPhRBc2bjA6sKEbr9PPpcnTg1u3Sow6pcoFHwqC9sIAmzbPgv0cZKwUNLr35r2qtDBXtM+pZzrYdsf/pWxLJO8lyNJkmduqKuWK0RJwu37W7R7HeIkoVap0qjWcJ6RS/5Z5HN58rk8/b7gZO8CDibzi9uUqqckqcV4YrC2mH1piJME27ZY+JglBE3TXi462Gvap2QYBpfOnefW1j2ymC5QSnFl40LWB75U4qTVemKaPUkS5hqz3Nm+z8ifYJgmedtmPBnz3es/5sr5i6wv/9XXzZWC7W2Tu3ctGrUZfv7nd5BUmAQupmFM29Y2sCyTuZkGi7NzH5sJoGnay0UHe+2VlaQpveEA2zQpF0sfGpU/S61c4Ztvvk2r2wGgUaufNYRZXVii0+9l1e5sJ9uk5zjMzzSIohjLMFCGiUBQLpRAwPXbH7A8v/DEF4SfVhjCjRs27bbB/Lzk2rUExWts7+8x9ic4ts3a4jLFQuEzv4amaS82Hey1V9Jx+5T7OztAtobu2DavX77yxOa7pwVhyP7xEUmS0KjVCeOIe9sPKOYLLMzOYZkmb7/2OqedNoPRkHJxgdn6DJMgIEljXNelmC+cvebY99k5PODG3VtsrKxRLH36nHcpJdv7e9x7MOL2LY/hKOTatZTElLR6CyzOznHpnE5Z0zQto4O99spJkoT7OztPtGBVSnHn/hZvTVPontYd9Png3h0c20Epxfdv/Jic67G+vEJ/OOTg+Ih3Xn8D27KZb8w+kVJXyOXIeXmGo8nZfcPxiLE/oVatcv32bX70/vtcvXgJx3F4/eLlT1zLv/twmxs3E44Oq7T625TrD9k+8UnFAp1BD0OIj03re9GIXg9zbw9ZLCJXV+ErUulP014UOthrrwSlFP3hkESmhEGA8YyNdGPfJ0nTZ/Zpv7+zjTcd9T/Y2+Ok3SZJY8b+mEZthkatzsO9vWeOpoUQfOftb/B7f/jviJIEJVNG/oRGtU4URcxUsv7w7X6XudoMt7bu8uZr1/DDgN2DfcIoopDPs7a4jGVZjEaSP/szQRSWsL0DSuZNLFuQpilj38cPQ27bWy9NsLfffRdrZwfpulhxjLpxg+iXfgk1bQikadon08Fee+lNfJ+bd28TJTGGMBiMhwgEjVr9ieOEEB/6EgDZ2n4YRXiuSxzH7B3tY5kmhhBESbbub5oG9lOj8d5gwMP9XXqDPsPxmLevvc5Jp8XRSYuL6xuUC0WiODo7PoxihCEY+RMGwyE3791BCOj0+0xCn+9d/zEL1a9xcjSP78OV1/ocdh7g90BM/0nShKJb4Pj09It5M58z0W5jbm8jp/sNlGmClNjvvkv01//6l3x2mvbi0MFee+ndfrCFYRhn6/ENu87Ne3eoVapnpWjTNKVSKj1zh7ppGGf56a1eF8e2iZMEqRSOZWGZFr3BgKX5xbPHTHyfH9x8j1tb99g/PmQSBliGycLsPG9cvsJoPOS41aLV6VAs5CkXS8w1sha0/cGA6+o2Mk3ZPTpASslxq8PxUQORtHntoqK+eJdafZHWwEIqhQEoFDnXQyn1lWlo81dlbW9/eARvGBjd7pdzQpr2gtLBXnupxUmMHwRP9GQXQrC6uMhwNMJ1XQRQKZe5/Fit+PFkws7hAUkSZ4G4PsNxu5UVwymWOWqdYFk27vR5wzhmfXH57PH3dh5w885tDlsnDEdDwCBRMcetU9rdNo1qDWEYKCT90YgwjgmiiJlKBcdxeLi3S2/YZ7Y2w97xkIOdDeLYpTazx8XLkiDKc9ppU69VOem08MOQWrmMQpGmKRfW1p7XW/yFUrkcJAk8NWvydOtbTdM+nv6N0V5qQhjwjIw6z/G4vHGBYj6PaRhPpL71hwNu3L2N57gIIThunWLZFisLC8RxTLffp1Iq0R+NebC7jVQKmUr+13/xvyFTievaBFHEaafDcDxCAIaRjbS7gz7FnMfecUi5VCJOEuI4RsqEnOuCEKwuLnFwcsLR6Qnb2wb37s0i8ak17pIvSVo9j/WlZVzHpZDLkcQJQhikMsW2bQSCOE74ix+9i2NZrC4tM1ufeU7v+OcruXAB684d1GPBXgQByaVLX+JZadqLRwd77aVmmSaVYolJEJxN2SulMA2DWrn8zGn7B/t7T6TgWZbFcDRGScXK4hKnvQ6nrRZKSpIkZTgeZzXyw6wl7TgUtLod1LQZTbYXIEUJMISBlBDGIRiCarHM/EyD03ab+dk5lufmsS2benmWP9+PGQ5yGHaL+YVDDCPBtou0uh0WZhrMzTRYWVhkbWmZo9Yp3UEfFLS7HZzHMg3uPryPZVnUypUv/g3/vDkO4S/8Ava772IMh2DbJBcvkly9+mWfmaa9UHSw1156V85f5NbW3awtrZK4rsfrlz66l3wURU/syD9pt2n3OhTzBR7u72ab9TwP23bo9nv0hwMUkjCWWJaJbVqkcUKcZiNulCJVKUqB6RjZaN8wiEOT06jNcDyiXCxlTXa8HIOBzdbdGcp5j0j9JQvzE+JUUSqUEcKgmMvT6nf55lvvAFllv6W5eZbm5rl9f4vCU2vcnuuxf3T4YgZ7QNVqRL/8y1mZwE9R+EjTtA/TwV576VmmyRuXXyNJEqSSOI91n3sWx7aRUgIQxzHtXgfTMCkW8rR7PYQQ7B8fUqtUSWWKH4bE0/azj8roShSmMDBMgzCOkUqhpERJk1RJLCUIkwh3mrefpAmO5bCz43F4UCBKujQW7lBVgjh1sWIT13GplsvkPI+NpdVnpgjGSfzMSoBJmn4eb+WXSwd6TfvMdLDXXhmWZTHxJ+wfH+E6LnMzjWcGzPXlFa7f/oC8l2MwGiEQGKZBtVTh8PSU03aLKE6IoojTTpsgDDAME0MYxHGCVCmWaeF6HkpKhCGI4gRMhW3bGIaJlClJnODYNtVSGZk6DLtXubU/wcvtkq885Pz6JSZBjiiOUUpRyOdZmpsnCEOWFxaeeY3VcoX9o8Mn0gBTKanoDnea9kp7LsG+2Wx6wL8H3Olr/s7m5uZ/12w2N4B/AdSBd4G/t7m5GTWbTRf458A3gDbwn29ubj58Hueqvbwe7u9mgd52SNOU7f09vnb5tSdqxg/HI+7v7hDHCQ/au8RxjB/6lJ0KhycnhFEICEqFwnRELs829wlAGGBIgVIyy733XAwEg/GIOE6y40Q29S4ME8dyIJ1nd7vKxjeXufzzJwRpCz+cpzscsLa4zPb+LgoY+2OCMKRaqVD9iCn5pbl5Wt0Ok8DHc1yiOMK27c+l2Y6maS+u59X2KgR+aXNz8y3gbeBXm83mzwD/I/BPNjc3LwFd4Nenx/860N3c3LwI/JPpcZr2mflBwP7xETnXwzCyAjiObXNv+8HZMUmacuPO7SwVr1TCNi32j47YOzpi7+iAh/u77B0c4rku5UKJhdk5CjmPhcYcF1bXWFlcplqqsDy/SKVUZrZep1auIJE4tk0+51HKFwCBaVjk3Bzj/iqHe4ucW5lh6dx9uuNb7B0fZjMGfoBSiovrG8zWZyjmirx+6TLXLlz6yKY9hmHw1mvXuLxxgVKxxMbKGu9cfeOZMxiapr06nsvIfnNzUwGj6U17+q8Cfgn4u9P7fxv474F/Bvzt6Z8Bfgf4n5vNppg+j6b91FrdDrb1ZK62EIJJGCClxDAMTjutswp6u0cHHJ2e4jgOY39CfzBAGIJEJvSHfX7u69+ikM/THw4ZjofM1hu4tkOrWGAS+CzMzpHzPEbjMXEc45VdDtvHWLbNTK1GpxPRPryA6zT45jdLzMy2GE0EhmFQK2ezCKdBmwvnzmEYBsV8gcvnLnyq6XghBI1qjUa19oW8l5qmvXie25p9s9k0gR8AF4H/BdgCepubm8n0kD3gUVWSZWAXYHNzM2k2m31gBmg9r/PVXi6u6yBl+qEGKoYwzkbJUZxtbuv2+5y0O0x8nyiOSNKEuZkGSZpSyOUJwoCdg11s20UpRbvTZeL7eK6HZVpUS2XKhRJJmjBbn8EwTQbDAUEYEoQhkT9POL7ChbVVfvZbDr/yC9f4nX/zrygVS8TEWKZFvVql3evQGw6xTIvZ+gyri4vPujRN07RP9NyC/ebmZgq83Ww2q8C/BJ6VKPto5P6sOcoPjeqbzeZvAL8xfX4ajcbndLZfDZZlvXTX9Gl8Edddr9dp9/uY05r2kAX3Wq3G/b0dgigiTmJuPdjCNE0OT48IgpBUZmvvUiocx2G2PkOqEg6PT9lYW2PGFJiWwfFpm1RKGvUK59fP8Us/+9f54N4dBqMhB9/7C8I44tLGFXbulwiCAvVaxLd/JuBXf+mbNGp1Lp2/wPHpKaaVfRmp12q888YbzDfmuHrh0ktT/vZZXtW/5/DqXru+7i/htZ/3C25ubvaazeafAD8DVJvNpjUd3a8AB9PD9oBVYK/ZbFpABeg847l+E/jN6U3Var1cA/9Go8HLdk2fxhd13eeWlrn38AFjf4IwTEr5AgeHh+S8rIDO9v4ek4nPYDQgTbNqdJNBn0K+wPG0PG4UhdQrNRq1Osuz81y/c4uT0zZSpqRpQpLE/PjmTbqdNhur64S+z8rCEne32uxuLSFjm+WVY86tJzRq6wipaLfb2IbJpXMbtNptAFzHYTwaU1zO0X3J68C/qn/P4dW9dn3dn4+lpaVPfexz2aDXbDZnpyN6ms1mDvgV4APgj4G/Mz3s7wO/O/3z701vM/35H+n1eu2vKud6fO3KVX7m7W/w7TffxjDEWc38KI5p9boU8jks22Fhdg4BFAoFBuMhvcEAlGS2PstwPCGKQgD2jg6JkxgAqRT7x8cMRwO6/T5SSg5OT9l+YNI/fY1Cvshrb/R5840clXLpLJcf4OLaOeIkxjRNXMchjELmZhpPZApomqZ9Vs9rN/4i8MfNZvM94HvA729ubv4r4B8B/7DZbN4jW5P/renxvwXMTO//h8A/fk7nqb1C0jRFCMFJu83WzkMGwyGnnQ790QBTGFQrFTzHme6kzxEmCSedU8rFPHGa0B8OSdOUNE3xA58oytrVZhX2XFqdMTeu57lz1yYVR9Tm3qdUzL4YeJ5HuVQ6O5dCPs/PfeNbrMwvUC2X+dqVq1xc3/hS3hdN014+QqmXZsCsDg4OPvmoF4ie6vrpBWHI4ekxpmGxMDuL81S3tMftHx1yb/sBe8dHOLZDp5dttIvSmImf1dLvDwekUmIZJo16HdMwCOMIgcBxHE5ap3QG/WyTiWHgmDb1aoV68QIPH1SR0iRXesjaqswq6UnJ8sICF9fP887V1584P/15v3pe1WvX1/35mE7jf6rSkrqCnvbSODg+5uH+Dva0BO3e0QGvnb9IvVp95vGLc/O8d+cWgqw5TqlQQCrFTK7Gj26/TymXx7Ed4iQm53kcnBxjCAPHtsnlcpimSW88olIsocgq1QVByP5uCVnbwDTbrG0c4XmKSrnKbG2GSRjwjdffYG1x+SNz5TVN0z5vOthrL4UkSdg+2MN7vFudabK1+5Ba5a2zwCqlZBL4WKaF57q8dv4iD/Z28QOfnJe1vW31uswez3BuZYXdo0P6/QGn7TZ+4ON5HkEUEqcJURTiWllLWQEEoUX35BxpXOT1CyOU8xDTElTLdQzDJOe5lIsF6pXqTxXoJ76PQpH3cvoLgqZpn4kO9toLL0kS/vAv/ozrtz/Atmxm6zNcPrfBcDTioHVKkqRUKxXKhSK7hwekaQoCSoUiF9c32D7YJU5iHtzb4e7D+6SpBAXdQQ+FYhIE+GFAqhRj38c0TITI+t57rsco9AlGNQadVSxTUKlvUZ9XSFkk7+WwbZs4SVBK4dgOxfyHN91Fccz3rr/HzVvvY5smuVwe0zQ4abUQQtAbDojihDcvX+Hta6/jOi7tXpd2v0ccRYRxjGNbzDdmaVTrH9nRT9O0V5MO9toL79/+6Z9wa+sevcEA0zTpdLsct1vMz8wgyEb4/UGf77/3I65euHS2Th5GEfd3t0kSSa/f5+adW1kevmFimRbD0QjDMEimu+0fkTLNmtMAo5HPZLiBihcwzB5e5S6WB3fuK371F36R7rCP7wfEaUK5WOLqxQ+Xuh0Mh/zeH/8+vf6A4XjAcDymkC9QzOdJ0pRWt4NAYBoG+0cH3N1+yMrCAgLB/skxJ61TEIJGrUYxX2BteYW3XrumS+RqmnZGB3vthTYcj/lg6y6One2aV0qBKdjd36dRq1IulTFNk3a/hzAM+sMhtUrWRMY0DB7s7jDfmGU0meC6LrblIES2ho+ASRhgGyYgsEyLJE1QQJwkkBbpnq6jVIF8YZt85QDHsSnm8kRJzA9uvse33/46pjC4tHGeRq3+zGv4Dz/8PqPxGMsyieIEz/XwA5/BaIjnuoRhSG46QxDFMSftFsPJkPmZWYIwQBjZQkIYxRhGwGQyZnt/jwtr68/vg9A07StNB3vthdbt9YjjBNd2qZTKDMdjkiQmTmOUhJW5n7SCNRAMxyMKuRwIgW1ZRHGEZZqEUYQQAiUlUZqQpglSKdI0xTVthIBUphhCIJUiGM3hD9eoVQrUZrcJ0zZ5r4JlWSipKBYLlIslFmfnOL+6/rFr7d1BnzhJcJ1sY6EQAtO0GIzHKKkQhkE6zck3TINJMME0DUb+hDCKMI1sBB9EAeVikeF4wnA8+sjX0zTt1aODvfZCK5dLiOn6tCEElWIRyEbeq4uLCCGQUtIb9Hn3/RsopTAEzFRrzDZmWZpboNVpI4SgPxgQJykKhSEMQGEYJkEckaZptuM+tQiGF0iiOjMzMevnduiPushE4gcBpmnieR71Sg3HsbFM6xM31VmmNX09MA0TNa0MbVsmhmkQRwn2tIyuABzbQSCwTTNbZkgTUGAJizRNcRwbU6/Za5r2GP1/BO2FVimWWF9aIggDpJJIJQmigEvr53Bsh1RKHh7scuv+FsIQgEKiOO606Q/6KJly1GoRxhFScRZos+dSyDSZBmtBGpfpn7xO5JdZWOrytbcClEiolEoIIUilIopjTMOgUa1hGRb1SoW7Dx/wwdZdTjotnlXX4rWNC9i2RZwkFKbr9Gmasr64QjFXyGr5T8+tVChQLVepVarMVOvkPY80TUnS9Kzanuu4LC98+jKamqa9/PTIXnuhCSH4W7/4H/OnP/hLdg4PAbi8fp5vv/11Cvk8O4cHfP/Gj8h73jRoK5IkBSkxhMHByTFzMw1u3b9HIecRRBFKQZLECMPAEIJ6ucrxUYmgv0ixkLC+0eHrb64xmoyZr88wCQJKhRLH7RNMIxvJH52e8J1vfJN3P7iJZzvYts3W9jatTpdrFy89cQ1vXHmNVEk+2LrL/vEh9UqFxdk5NlbW8TyX3cMD7m4/IE1TGrUaVy9cZH15jb3DA+I4RqFI0hTXsVlZWOLi+gb16b4ETdM00MFeewmUikW+/dbXqZS2SGXK8sIi1XIZwzA4v7JKpVgmChMGo1E2gpeSKIqYhAG1cpXdowM6vR6pzLotm4aB6bokcUoaW5wenicOCriFQ2bmjxiFCVs7Wd57mlNUiiWKuQKObdPqdugO+9CF//uPfh/HcVhbWKRSLrMyv0h30GcwGlIu/qRUrhCCt6++zq/8/N94ZnWtN688q0EkzM+8el3DNE37bPQ0vvbC2zs64N33b2JbFoVcnv3jY/799/6CwWiIaZqUC0WEKQjCgDiOCcOQJElI05ROv8toMsa0DKRS0816IUhFFFXpnlxjMrEozjxgbvEQ2xI4lo0QBtVyhYk/5rTT5rTbzgr2hAF5L3f2OmmacHB8RBTF7B0f4to2nX7vy37LNE17xeiRvfaVNfEnHLVamKbB4uz8h+rcSyn54P49/r/v/gWpTMnn88RRTH88JAojfnDjOvOzs/zs2+/w3es/pjcaksQxSZoCMJyMn3g+AYSAUoLTk3niYBHXC8lX3wMjJAgNFJD3crTa2Zp/tVQmTMbkPY8oilBSYVomMpWEUYTjOkipCOOIVKYEUUgxl38+b6CmadqUDvbaV9LO4T67hwe4tjNtHXvE1QuXqJV/shZ99+EDfnDjPYaTMUopTrtdBsMB8zOzhEmEMATHrVN+9w9+nzRJKBeKjCcT0mn52WwF/ycUIJMco94FLFGjVDkgV9whlSlhJDGEwLKsrHytUriOQ5oqivkcnu1gWRZ+GJCkCalMEUZW6lYIgzRNznbQ16u15/12apr2itPT+NpXThRF7B0dknM9DMPAMk08x2Vr+yEAk8BnOB5z8+5tojhiOBox8X26/S5BGHLSbSMw8BwXx8rW0REKz3Yo5vJY1qPvuE+mxAXjBqPOmyBd8pX3KdcOMEwD0zLxHAcAmUoSmbXGRYASklKhiOO6fPONtyjk8rhOdt6u7ZL3clhCEIQRlVKRyxsXiOInK/JpmqZ90fTIXvvKafe6Z3nnj5sEPt9970dZ9Tql2NrdJk1TyoUiA388bUejmPgTqqUytvWTaX/btkkSSaoUj8bzZ2l20mTcWyf063i5PpXZXUwjQggbpbL2tsC0oI7EUiagME0TlGCuMUucJliWRa1aZe/wENM0EYAwBBc3NjCEQblY4v17d0FAtVzh6vmLuoa9pmnPhQ722leO63qoacW4xx2ennBuaYWxP2E4HDIJJvhByOLsPFXTQKYpJ20f27RAKfwwwDBMlucXaPe7FHMehiEYT8bESYJlGPhBjknvAjJ1qdSOqNSzPQJRBFESo6TCtgT2NGffMU1c18WzHVzXIeflcB2HvJnjW2++nZXGXdvAD7Ppe8/1kDKl0+0+sQN/PB7zYH+XC6u6pK2maV88Hey1L02apuwfHzPysw1uK/OL2ei4XJ5ubJNnI98wirBMi52DfQ5OT5AyRSnFaDLmpHOKYZj4YYhlWkgBJ51sd7xlWVRKZQbDIf3RIJsVAJSCQX8Of7CMacU0FrdYmHcYjy2iJEFKRSrTrLCOn5AkWU18UxjZsoFpYJk2c3XBj27dwLMcHuxtM1OtMVOpsTg3z1HrlO39PTqDHtVCiR/cuI5lW6wtLhKGMQenx+Rcj4XGrB7ha5r2hdLBXvtSJGnKj96/QSoltmUxGo85arV45+o1hBC8eeUqdx8+YDDtPNeo1Rj5Y370/vt4jkMK2JaDKQS9wYBysUQcxViWDShc22I08QmjiNN2Cz8Isk1zgEwtxv2LRH6Rej3g9TcUpjnH2tIyH9y7Oy2LmzCaTBiORmdBv1Qo4ochhpGV4PUTn73jI1ASx3Fweg6tThfHsShs5cEw8GybdqfN3uEBjVqNvJfj7sP7rC4sMVOtsXt4wEm79ZG59JqmaZ8HHey1L8X+0SFyGugha0MrpWR7f4/lpWVsy+baxctPPObm1h0MIRj5Y/wgJE4iojQlDENsyyKMI1zHIUkSwijEdlx8f4JpmEgls251QYVRfwOkTXV2n1p9yEztGm9cusL68jKdfp/do33SOAUhKORzgMB1bEqFMqedFqbrImVKEEWgIE6zGQSlVLYT30+J4oS3XrvGaDyetsw1GU18ZCpxbZvjdou1pWVsyyIIQ1p7u8zdvo1zeIgqFkneeAOV1yl6mqZ9PvTcofalGI5Hj+2KzxiGwSQIPvIxG0trSCUZT0fsUZzgWBaObWMYBoYhMIQgThL8MMyK2iQSKSWphHF/hUH7EoaRUJu7Sa7YIlWSMAqRSjLxfcr5/LTTXLbd79FMQZoqxsEEENMpfgkq622Pym4LBEmaImVKGEUEYYgfBFiWdbZ2Pw58pFQIoaiWygDkkgTr3/xrRLeLSBKMTgf3D/4AfP+L+wA0TXul6JG99qVwHZcgDJ9Yq1ZKnRXOUUpx1Dpl7+ggK4E7u8Dq0tJZjvpwMiaHS6vXIYxjZJISy5SJ72OYJkqpLBALRRhaDE7PE8cFcoUWxeoepgWQrbsrpZgEAatLK3z/+ntZJ7kkwTRMkjRrn5vIFIOsjaxpmCgps+dHgMj6yUdJhAwlQggc18EyTWzbJgxDwjjCEIKcl6OQz+O4ztmsRv72baxiGQwDMR5jHB+DELj/+l9j+D50uyRvvUXy1/4aKpd7vh+UpmkvBR3stS/F2tIy795sY1s2o/GYJE1wHJcrGxcIgoDf+4N/x/29XcI4IAxDTMNifWWFjeVVhuMRSRwzGI0I4og4jnk8c/1R7/c0TQn9GfrtFZSSlOr3cHM9FKBk1kq2XCpiGiadXpcf3HwPPwxxbJskic9G7LYFoBhMXzdSUdYkB4FhgFCKMAhJkVimiVAKJSX3dh5SKZYYTcZIpci5LoZh0O73uHb+IpZlZUsZvs/Myhrcv49z9y7StjFv3ya3vU26uAieh/v97xN/97v4zSZydfVL+MQ0TXuR6WCvfSlcx+HS+nn+3+/9OWEUYVsW9WqNwWjE+w/u0e51SWVCt9cjjGMMkU2R39/ZZn52jjhJGYxGSPXhFD2BQCqDZPIa/qBKoTBkbuGAOA0wzSJKSYqFAuVikfnGHDnPI04S/uy732UcTnBsB0OY2K5NvVqj0+9RLVdQEnpWHz8MUDKbhfBcDyGgO+hjGzaeY1Ov1inkcnT6A1BDZqp1lJKYloUAVheXKU3T8AqFAueuvY552kLs7yNzOUSrhdVqgWFg+j5ppQJSYt6/j/O97xEsL4Peva9p2k9BB3vtuRgMhzw82CMIQ1zH4dzyCntHB5xbXp22ns28+/51do+O6HQ7tHpt/CBESkmSJPSHAxzbxg9DBsPBM3vDA8Sxx7BzAUMVacx3qNZPMEyBqwrkXQ+pFK7jUi6W2T86Yq4xQxzH2K6DlUSMJ2Mc2yZOE8Ikmn7RSLAMi7yXY6ZSJUVRzheZb8ygFGwfHlCdBnCFmi5RCAr5Aktz80yCCUIISvkiruuyXq3yrVTBYIS8eAnj+g0QApRCDAYQRZDLQZq148UwEEmCeXSEGI9RycQHtgAAIABJREFUpdIzr13TNO1ZdLDXvnDjyYQbd2/huR6WaZKmKe/dvkUcx3iuy2mnRZKkpFIyGA6I4gClFOOJTxSFWSrcdGo+jGP8MPsC8GRl+4w/mmPcX8E0Uyqzd8iXY0zTI0lS/ChAyuwx3cGA0XiMZVnEcUK716NWLmNbFpFpYZgGpALXdhjIEdFwjOvYSJWl5Tmui2EIhGFioH5STlcpesM+aSqRKiWKQm4/3KJRrQGC4WjMRcNifes+5voGCIF1/z7x5cuo27dRvo+qVknn5zF7vewLwPRLAID0PNT0tTRN0z4tHey1L9zu0QGu4z5xX8512T3cRyqJUAIloNvvEUQRszMz9PtZAZxUyidG/sDZfY+P7GVqMeqdIwqqOF6Pmdl9ymUXzykQJwmJTJFpOk3xS/FcBz8KcJSDIQSmaTAJfJbmF7m9dRchspa3ge+Tcx3iOMu1d2ybOE6I4hhzWrO/Wqng2A4PD/aJ4wikQqaSmUo2fZ9ONw4W8wXiNEUcHTHz1tfBzMrwqmIRczCAa9dIJhPEuXOI0Qij18s25AkBYUi6sEBy9Sq4T76XmqZpn0QHe+0LF8XxhwK2YRj0RwOOTk4ZjLNGNnEcYdlZ5TxzWvJWwYen65XCMAwerdbHYZFh5zxSWhQqu+SKx+TzJTzHxTRN+qMBj5reJGmKbdkoFKZhIAwDpSSe6+E5LnnPwzQtEplgCAPHdTBSiedkjW8EZKN2KUEJpEyol8pc3bhIoVDk3ZvvYVoW6yurLM7Ns3d4wNifEIYhjmOzXqsxn0jU07X/0xT1zW8i330XM46JfvZns9H9wQFGGJKsrRH+2q8Rf+tbX8RHpGnaS04He+0LVy4WOWm3saYjWchG561Ol7E/YTyZkCQJUiniOOK4dUocxUimQf2pOvmKLGgrBf5wiWC0hGFF1GbvkMtHOHYRIbJjhGGytrRMd9DHFCapkkRRhGma1MoVzq+uk/NcRuMxcZownkxYaMwy9ic4tk2UJCQypF6tsjy/wPL8Ils7DzENk8sb57NrSVPG/pi/+Qv/EWsLi2dd7fwgwHM9SsUiOddjeX4B4hh5eIT11AY7IQRUq8Q/93Ponniapn3enkuwbzabq8A/BxYACfzm5ubmP202m3Xg/wDOAQ+B5ubmZrfZbArgnwJ/E5gA/8Xm5ua7z+Nctc9XkiSkqWT/5BDbtJifmUVKydifZMVrHk3HC4FpmMRxTH84xDQEOdvFj7J690maPPG8aeIw7G4gZI2FeZ/F1T5B6DH2FQpFzstnxXLSlOriEq7rMRwNs854uRx5L0+pWGRuZoZysYRhCN68co07D7YYTiY83N3Gsm1kKjlqnZKmKTPTnfli2rXuEdM06Q0HAKwuLHHj3h1yrkvO83AdGz8IWF1Yyt4PQzCzsIiVpj/ZUZ8kpPU6lEoQhs/ng9E07ZXyvPJ3EuC/2dzcvAr8DPAPms3mNeAfA3+4ubl5CfjD6W2AXwMuTf/9DeCfPafz1D5HcRLzg5vXafe6rC1kJXB3jw+Zb8xyYW2dfM4ljCOklAiRbW4TQiBlShTH2I5NuVjAtq0nlgFCv0rv5BoyLlCuP6A2d0CSRuRcD8e2qZZKFPJ51pZWWFteoV6p8J23v87V8xeplysU80XevHKVb735NquLS6wuLvGN19/MutflctiWxdrySlY8B8XszAwLs/N4rodSkmq5ymy9/tTVZudXKZe5snEeIbLr31hZ4+qFSxiGQSpTGrU6G7/2nyIXFxFJgohj5Ows0Xe+8xw/GU3TXjXPZWS/ubl5CBxO/zxsNpsfAMvA3wZ+cXrYbwN/Avyj6f3/fHNzUwF/0Ww2q81mc3H6PNoLYufgAMMwMI2sU1x/OASleO/2+9SrNVKpUOmjAjgSmaZkFeyzgC9lSrVUATFGCEHgx4TjdYLRLJ43oTKzheslIIo4loVEUi6VmavXMQxj2hwnmm6sc7h07jwbq+tUyxUun9t45jmvLi5x3G6Tcz02VlZRSjEcjSgVS5imweriMu1e94kvH0mSUK9Wz243anUatae/DDwp/vrX9XS9pmnPzXNfs282m+eAd4C/BOYfBfDNzc3DZrM5Nz1sGdh97GF70/t0sH+BjP0JpmFw/fYH7Bwe0Bv0GU3GCLL8cz8MmPgTkjR94nGpzG4HYUTgRDiWyWhoM+m9Rhx5FCsneKU9lEoJI1BqSJqmlEtFBsMRpXyeYqHA/EwDx3E4Oj0hiCIEikqpzMW1j+4h79gO71x7nYd7u/hBAAgkWZqflDKryheFTHx7mhUnqJbLnNd96TVN+wp7rsG+2WwWgf8T+K83NzcHzWbzow4Vz7jvQ0nVzWbzN8im+dnc3KTRaHxep/qVYFnWC31Nc+1Z7tzfYvf4CD/wieKYOE1RUhINE1zbIe/l8KOQOP7JONcQBsZ0g12r38EVGxSMr1Gah1j8kJQj4rMZgYQQhWmGKFmkXq1iOw6/+O3vYJomY3/Cd77xTbqDPp1eF9MwafV7XDl/AdM0SWW2Yc+xbczHNhCuLq8AcP32B5TGJYzHRvKlUolyqcjGyjq2ZeF8TnnvL/rn/Vm9qtcNr+616+v+El77eb1Qs9m0yQL9/765ufl/Te8+fjQ932w2F4GT6f17wOMFwFeAg6efc3Nz8zeB35zeVK1W64s5+S9Jo9Hgq35Nw9GI/mhIqVCk8lRVt0qxyI07twjDgCiO8QM/qxCHyKbsTYmS8kPf7BSKVCqkNEkmlwnlInNzQ/7Wf7LADz/wuPMgxbGzXvaRGfPou6HrONRrNVzLZmtnmyjKdtH/2z/5I/L5PPVKNtU+GAw4PD6iUauzd3RImiYYhslCY5ZzK0/WnT85zTbnPS0IAhbqs4ynOf9PpxZ+Fi/C5/1FeFWvG17da9fX/flYWlr61Mc+r934Avgt4IPNzc3/6bEf/R7w94H/Yfrf333s/v+q2Wz+C+DbQF+v13+1SCl5f+sug9EIyzRJkn1K+QKvX75y1sku53oszc9y0j7NHiQElmGQpBIpFXGSINME+aE0ekUcFRh2zuOYRWYXjlhYTpCqzOrCIp1+D8/1GI/HEIDjOpgYzM00qJUrHLdPma3VcR2HVEqOTluMJhPKhSKWZWGaJiftNp1el0qpfNZ97rB1gue5LDTmzs4l53oMRsMnuvOlUmJIyQ/fv44fhpimxeLsLGuLy1/sm65pmvYZPa+R/c8Bfw+43mw2fzS9778lC/KbzWbz14Ed4D+b/uz/IUu7u0eWevdfPqfz1D6lw5NjxpPxWZlYyzSZhAG7RwesL60gpWRnf4/RaEK718ex7Wmp3Kw1rCDrBZ8qhSnEWcBXCvzRApPBMqYZU5u7w8JigbxXZWtvhySOMYWBZZpnMwnCyHrMFwuFabEbhTs9rzRNUUphmSadfo+5mWwKrT/sM1N9chOdazuctFpPBPtzKyu8+/4NTMA0jOm0f9anvpjP400rAx4cH2EaZpZLr2ma9hXzvHbj/ynPXocH+OVnHK+Af/CFnpT2mTwqcNMZ9LEt++z+NE0Zjkf4gc9CY44//vM/44cf3OC008YPfAajAYZhECcJlmlimCakKbZpooQAKZGpxbC7QRxWcHNdao09qpWsO11n2GfRnaNaruCHIUmSECcJtmMTxzGzM3WEAtu1WH1shG3b2Ug+S+lTT1xHIZ//8PU91UXPsR2+fu1r7B4e4Ac+uVwOmaT0RsMPHXfcOtXBXtO0ryRdQU/7VJIk4fbD+wyGWZDrDno0qnUsy2IwGnFwfIiYptn99r/c5IOte/QGA1ASYRgIIYiTnxTGyXsefhBiGQLTslDxDP3eGgYGjfk9HO+EfC7H0tw8tmlhGgaLs3PMNxqsLixxa+seI3/MUqlM3suxtriEZVvMVOtIKZn4EwzDQCCYn2nwcH+P5fkFlFKEUcjljQv4QfBEVb8kSZifnf3QtTu2zYXHdvDfvr+F+YwWs4+yCDRN075qdLDXPpX3790hiuOz6fFirsD93R0undvg6PQY23aI45hEpWztbNMb9LMiOQjSaSncRyzTIgwjQDEJE8Yns/ijBWxrQrF+DzuXknc9HNvCMExq1Sq1SgXDMDCEQT6X41tvvc1wPObNK1epVSoEUYhtWdiWTZIkXL97Owv4wsBzXX7hm98GshmIy+fOUyoWubf9gJN2e1rFT1IpV84q3X2cuZkZPtjq4j3WkEZKSbFQ/HzfdE3TtM+JDvbaJwrCkOFkTM71zu4rFgqUS0X6wyF+FJF3DUrFAt+/cZ3eYHC2Vi4M44lADxBEWUnYNHEZdi6SxAW8wgmF8h7CkMjUYBKGuK6LYRjM1mbojQYctU7Z2t1mtjaDaRqUCkWq5fLZCP641cLzPOrlCu9cfZ3RZEwYRVSKJSzrw3/VL65vsLq4zHAyJu965HO5T/V+1CpVZmo12t1O9uUiTbEt62Pz9zVN075MOthrnyhNU1AQxRFhFJHzctkGuWKZyxvnse1sRP2DGz8mjuOsJ7xhEqcJ6qkmNo+Ekzqj3jqgKNXv4eZ6Zz9zLBvDNAnCiMFowPdu/Ih6uUp32KeYK5CmaVbQxhD0R0N6gz6Hp6fY07x527J487WrFPMFivnCx16b6zhnsxU/jSsbFxjNLdDudynksrQ+4xlT+5qmaV8FOthrn8h1HA5Oj4miKGskA8zUalQrFWaqNYqFIjsH+6jpsZZlkfM8pD8hfSrYS2kw7q8RThpYzpBS7QGmFX3oNW3LnLbGNUiSBGEIPMcjTVNs0+Lcyiqu43Br6x5JmlDMF4iiiMlkQqvX4d7uNqsLiyzONlhbXH3myP6vqlgoUCx8/JcJTdO0rwId7LVPtH2wT71S5bh1ijHtw35wckypWOC92+8jU0lv2Kc/HBDHMYYQWbe7pwJ9EucYti+Qph750gG50gHPqkUTyxQrzerZA5imwWg8wg8iUpmyMDtLq9PGD0O6gz6e7dAbDoiShNNuO3vtJKFQKGAIg42VVX72nW/w+qUrX/h7pWma9lWkg732ibqDPrVyhXKhSH80RErJ/d1dvv/eeziOg2kY7B4d4DoOSimSNGESBk88hz+aY9xfwTASKo1b2O7oI19PSYlCsji7iB/6FHJ5pMrS4vI5D9dxuPPwPvVqlZlKld2jQ3qDPkmaItOs/G2UJORyeTzP5rTb5gc3rjNTrbEwO/eRr6tpmvay0sFe+0SPpu5N08QwTN774CannTapUuQ8DxTkvRxHpyeEUbau/4hMLUa9c0RBFcfrUaw+xDCTj3opANYWl7Esi0I+x7mVFcb+BCklYS6PaZoMRiNAIJVgdmaG2w+2MAwDfzLCmE77W6ZBksQgPMIwwrYt7m4/0MFe07RXkg722idq1Gps7Tzkg3t32T85ZuJPkDLF83JEUUS72yVK4g91r4vDIsPOeaS0KFR2yBVPPuIVnhQlEYZpIqUk7+VwbJs4SbFti+W5BU67bTzXpZjPpumr5QqtbgeBQSrVNN9/mj+vFIaZLT08vaygaZr2qtDBXvtEg9GQP3v3e7R7fcbBBJVKhICR72dT508FUaXAHy4xGS5hmgHV/7+9O4+RI7sPO/6tq6vvnvviMUMOyeVyuavVytKuLTiRZGcjO7LlIHZZCmLLhhI5gRQnQBDYDhJISBxEAZI4SuIIUGxFkpFYfnDiWI7l2Jas05IleaXd1S6XXB7Lc+6j7+6680cVR0NyhuQuhzPDnt8HGLC7urrq/aqa/et69Y7h05iZ9l3vb3F5hWw2i6ZpPHfmJfrLFXJZm0eOPETH7RIEAYurKwxU+sjncvRXKniBT1+5wmq9BsR4QYChG0QxDFf68IOA6QMHt/jICCHEg0GSvbitl185z5/8+ZfpuB66FpO3s3i+j+t7eOumpb0uDDI0Vg8ReCXs/BKFymV0/dVdUfthgBkGXJ2dRdOhVCjieh6f+9pXmBgZpZDLMz40wsLKMmgwPDBI1+2iaTqWaaDFMR23S8ayKeTzFHI5Hjp0iP13MWCOEEL0Ikn2e1TH7XLx6lW6bpdMJsPUxP5bxoq/OjvL82fOABBGyTSwnu+uTVOr6/oNV/Vup4/m6hSgUeq/gJ1fuevy6Jq2NviOoSXbNQyDMAxZra5y4sgxWp0OuqEzPTmFZZpEUUSz3eKJRx7lh3/gB1mt12i12xQLBeIootXpJK33h4bJZe9uwBwhhOhFkuz3IM/3ePbUi+m98IArszO8cOY0jxw7zujgIIVcnnwux+zSPIVcFg0Ny7AI/TCd+S1EJ0ZHJyIijnRa9f10WyOYVovSwAUM073r8hi6nvxwSMfOj6KIMAjQrXR5HOH6HraVIQzCtSlpdV0nY2XIZ3Pous5gXz+Dff1r2x3s33B3Qgix50iy34OuzM5gmSadbpdLM1exTJNqo87n//zLDA8NcXB8gv5yBc/36a/04Yc+YRSyWq8SRhEayZzzuVyeWjWgtnqY0M+TK86RL19D0+I7lgHAMgxs28b3feyMTRi2iOIYK2MR+CF+4JPNZCgXS2QskyAMbhnS1jQNGblOCCHuQJL9HtTpdtF1nfnlJTJWho7bpd3pYJoGGjGe59HqtJlbXGSluoqGRqvVRjcM/CBA1zRKhSKrq2Vqi6NEWkh58GUy2fqm+9RgLSlbpkWxWGB8cJhKuQJxnAx7W6vTctvYVoZ2p0MYBui6ycTIKL4fEIQBfeXy2ja7nsvBsX2S7IUQ4g4k2e8BcRwzv7zE0soyXc9lfmmJjtuh1mhQyOWZX1yk3qyDphFFIdV6g1w2S61RZ7lWpdPp0HE7hGGIaZiEkc7VyyN02n1Ydo1y/0UM49bGerqmY5lmMludrpHJZKgUivRVKowPjWDZNgdGx7h07RojQ0N4ns9KrcrlmWsYuk6lNMBApY9qo042Y/PkY6+n0W5zZW6GqX0HODp5iOGBwR04okII8WCRZL8HnLlwntV6jU63w+zCPGEc4QcBtVqNyzMz1Bt1/CAkn8ty4coV+krJTHIr9SrdTgfDsui4yT143yskfefDDPnyVXLFuQ2HvIVkxLsgDMmaBuViGU2DcqlMq9Ph3JVLjA+PUqvVsEwLPwxwXRc/CJjcvx/f9xkdHGZxdYXHj59Yq77vK1fwfI/RoWFJ9EIIcZck2fe4dqfNcq1K1s5weeYamUwyB7uOS7PTYWZ+FmJA12l3Wui6geu6tLudtdbxWreb9J1vjtOuT6AbHpXh01iZ1h33bxjJ/PONVpP+coWO2yUMAvK5PFfnZtZmixvs72d2cQHLNCgUhugvVQijkJxts1yv3nCvPmNlWF5ZYWxo+L4cMyGE6DWS7HtcvdnA1HWiMCIMI3RdJ44jZhbmqNaSaWUjYqLAR9M0wijC8z3idXPQB6FFc/UQvlvGzq1Q6LuEroeb7XKNaRj0l8u88bHXc21+DsuysHQTz/dottuARtd1cX0PTdPZPzoGwOjICN2OS39fhZVqleCm/vxxHN+XWeyEEKJXyTdmjyvkCoRxhKmb6EZS397qdJhfXqTV7RLHMXGcDDELEEVh0uAtjgmiCK9bprF6GGKdYt9FsoWlu963YZhUShU0TUPXdHzfQ7O0te50YRhSazQwLRPP9+kYBnEUcezQNPV60tgvjmOWV2/sr+96LscPH9miIySEEL1PmjH3uFKxSKlQJAhDBir9+IHP8uoqoGEaBno6uU0Uhvi+TxQl99ljdFq1A9SXj6XV9qdeVaIHsEyT0aHhdFrcmIFyP4VcjiAMsUxrbUAfQ9PJZjL0lUpEUZROdJN277OzPHHiMcIoxPWSGodjh6ZvGQBICCHE5uTKfg945MixdApai1zWptlusbCSJO4oCAmikDAMiQHTMPFdk+WVwwR+nmxhgUL5KtpdDnmradraLYCMadJ1u1TrdQb7+tk/No7r+zRaLSDGNEyCwKdULFIuFjFNk9efOEkcRRi6TsayOH74CIV8nmkmiaJIutkJIcRrIMm+x0RRxOLKMtV6nVw2x8ToKKZhMDmxn8mJ/UAyHe2psy8ThCFh2mIekr7w7WaFxupBIKY0cA47V31V+9dIEn65WOJ1Dz+CoetkbZt3vPWvkclkmFmYZ9/ICI12m5dfuUA+l1ubwa6vVMbzfZ44+RhafOvAPJLohRDitZFk30OiKOL5My8l1eNWhlqjzszCHI8//AhZ215bb2F5iSiO0TTWWtxHkU6rdhCvPYSd7VDoP4+md1/V/k3DwDBMSoUihw8cYGRwiH2jY4RRxMLKEocPTDK1bz+Q/OiYPjDFhauXyKY9BIIwaX0/NDDA8vLy1hwUIYQQkux7ycLyIl3XJWNlADBNk67r8uxLL/D4wyfXEv7s4gK2laHruURhRODlaawcJgxt8uVZiqU5dFPDMmxc30fXdcIgwDQt/MDH0JIW/RFJVX02m6WQyydX6Lk8g/0D5LM5ioXkvrquaQTpuPfrjQ0Pk8lkuDY3QxiFDJcHODixD22zjvtCCCFeE0n2PcAPfGbm5zlz4TyZjIVpGniez9ziAl3XJdZiwihmsK+fiZERqrUa1WaDKArx2uM0VsdAD6gMncHOtojjGEPPYGcy6JpOuVjCMAzKpRLtdosoimi02uSyNuMjozx2/ASz8/OUiyXQoNFqYtsZyoUSAF3P48jk1IZlH6hUGKhU7jrWOJ2+Vtf0G2orhBBCbE6S/QOu2Wrx9WefwfV8mq0mV+au0XHdJJG7PjGQtW1WqlXmlhbpdj3anTauG1FdOojX7SOTrVLsu4huBJiGRRCFgIbn+ei6xsjQEK7n0e12sW2bfaPjeL5Po9mgv9JHGAYcPXSYSrGEF/jMLy2iazpzi4vousb05BSVUvlOodxRo9XkzIXzuL4PxORzeR45coyMZd3ztoUQopdJsn+ANZpN/uALf0IYhGiazsXZK8wvLmFZFmEYrv0ZhsHM/Byu7xEDXrdEbekwUWRSqFymUFwgJmlYF8URpUKBYi5Pq93GtCxq9QbZbIZao0FfuUS1XsfOWDx0aJpcNsfgwACvf/iRter3lzMZLs1cxTQNsnaWWr2O67nYmdd+JR5FES+eO0vGNMmlV/RhEHD6/FkeO35iC46mEEL0Lkn2D5j5pUVOnT/HuUuv8NKFsyyvrBJFSbc5XdOSkeXaybj0kHaFg+SeeQztxj7ajXEMs0v/0EvYOR9TtzEMnUqhhJ21GezrxzJNVqpVGu0WK7UVct0c+UKejGVj6DrFfJHxkVHCKGT/2Nhaoq81GiyurjLYN7BW5jiOOXf5Eo8cOfaa467W60RRyPqPrK7rNNtt/MDHMuXqXgghNrMtyd5xnI8D7wAWlFIn02UDwO8AU8BFwFFKrTqOowEfAX4UaAM/p5T69naUc7c7dfYsf/ilz9FoNbl8bYYgurHRWxTHeDcNLXu9K1wYZGisHiLwStj5JQqVy+h6RBRq+HFMEGrUjTbDWZu+coVGq4nr+3i+jxcExHRpdtpExYC+cgnTNHF9j1w2x9C6xL6wvEg2k7mxDJpGs3XncfRvJyZGY+OGexv00hNCCLHOdnVc/gTw9puW/TLweaXUUeDz6XOAHwGOpn/vAz66TWXctfzAZ3ZxgT/80uewMxnqzcYtiX4zMdBpV6gunCD085T6L1Dqv4ieDpKTdMHT0DQIPI84gpValdmFOTqdDkHgo2sahpaMuNfxPDpuF9PQGRsa5nXHT9zQ/z2TyRBGtw7AYxrGPR2D/vKtjfjiOCafzco9eyGEuINtSfZKqS8DKzctfifwyfTxJ4GfWLf8U0qpWCn1F0Cf4zjj21HO3abrdjl/+RLffP5ZvvGdZ1it1Vit1/E8767eH0c6zdVJGitHMEyXvpFT2PkbT4OmacRRhAZYVgZNh6RxXoDre8m+NA0vCPB8nzAKKBWLvPkNb+LI5KFbkvjEyChheOMPEc/3GRkcuocjkVTZH58+gh8EdLpdOt0uMXB8+ug9bVcIIfaCnbxnP6qUmgVQSs06jjOSLt8HXFm33tV02ew2l2/HeL7Pi+deJibiuRdPkctmqTbqtDpt6q0mnc6dB7sJ/CyNlWnCIEeuOEe+fA1Nu7W++3rFeC6bjGSXs23KhQKr2SztbptsNkccx8mY+UFAIVdganw/+8c2/v1lmRYnjx3n/OVLdLpddMNgYmSUA+MT93JIgOTq/o2Pvo5mu41h6BRyMj6+EELcjd3YQG+jG7Mb3pV1HOd9JFX9KKUYGrq3q8fd4pnvPkchn2dxdZmhwQEWlpeZW5rHMk18t0u08eFY020N06odQNNCyoMvk8nW115bP3Y9JNXruWwOwzQwTZNysczo8AiD/f28ePYM1VodTU9mqhuo9FPI5ziwb/9tj/XQ0BCHJqde81j2pmnedvsjm77yYLtT3L1qr8YNezd2iXsH9r0je03MO44znl7VjwML6fKrwIF16+0HZjbagFLqY8DH0qfx0tKrm5VtNwqCgGtzs2QzNs1Wi+dOvciV2Rla3c4d3xtFBq3qJG5nAMuuURm8hK77ROt+G+iaBmnL+TCK0HUDO2PT7rbJZ/NoGszMz/L6hx8hm7G5OHON5ZUVDF3DNA3KxRLFXI5rMzPYNzXE2ypDQ0P0wrl8tSTuvWevxi5xb42JibuvMd3JZP8Z4D3Ah9N/f3/d8g84jvNp4Emgdr26fy+IiSFO7tefu/gKs4sLd5Xofa9AY2WaKLTIl6+SL86h6UliNzSNMIowjSSxR3GcVp/ETI7vI5PJYFomBTtPvd0kaLb44je+jh8GZCyLwf4+TMMka2fQNI3LMzNcnZulVCjx8JGj99z4TgghxP21XV3vfht4CzDkOM5V4IMkSV45jvNe4DLwU+nqnyXpdneOpOvdz29HGXcLy7TI53K8ePYMS6srbNCw/QZxDJ3mOO36RDrv/GmszPe6uenpVbxlmui6Ts62MQ0TL/CYnNjPW578AYIwZGZhjna7jR8GhFZI4AfUGk30gka7M3f9AAATyUlEQVSr3WF0aJCFpSWOHT5CPpcDoOu5nL30Cg8fPnLfjocQQoh7ty3JXin17k1e+qEN1o2B99/fEu0ey9Uq5y+9wqkLZ2k0mpimycTIKH/5/HOsNusbTiBzXRhaNFcP4btl7NwKhb5L6HqIoetJ3/owJIoistks/eUKo0PDLCwu0HZdpienePz4ozTbbTzfY7VaQ9M1Vms1NE3DD3x0Q8cwTMZHhslYNvvGJnA9d23/hq5TrdWI0+57Qgghdqfd2ECvZzVbLS5eu5oMHWvbhGHAhcuXefb0i1RrNbzAR9d0vvHcdwjCAF3bvHGb1y3TWD0MsU6x7yLZQnIfyND1pMreNDEtCy0GLdaIwohOt42VselLh7B9/vR32T+2n4mREeqtJgsrSxDHZO0s7U6HmJhSoYCmGcnMd1F0Q+M+IYQQD4btGlRnz2u12zz/8ku4noumadTqNZ554bucu/wKrucRRjGWadHpdgnCAA2I41vr8ONYo1U7QH35WFptf2ot0UMyhryh6ViGQd7O0l+pkLFM/MCn1mhimSamYeAHIY1Wm2vzs+hpTYBhGARhSExM1rbxfX9tNLx9o2N4vkcpnckOkgZ+lVJJruqFEGKXkyv7bXJ5dgbbyqwlxpVajcvXrjKzMA+aRhRFaJqGFyTD3W50/RwGNo2VaQI/T7awQKFy5Ya+85Zhgq5BnAxCU8zn0XWDjutCHBGnzfLjtIFeEIb4gU/H7dJfKuN6HkEQ4LrJoD39lQqarjPY108YhRydPAQkU8wSQ7FQ4NjU4ft52IQQQmwBSfbbIAxD5hbn6XoulWKZMAx59tQLrNRrRFGEFwRrE9dsptseoFWdBGLKA+ew87Ukuccxsa7RXyqTy2YBsE0L3bTouh00NCzLxDYt4rSPfRTF5HM5gjBE13VsK0M2a3No/wFGO0PsGxnFtExiNI5OHaJcKGBnbOxMhiiKaHU6af/87DYcPSGEEPdKkv191Gg1OXvxFc5euki308ayMqxUa1TrVVqdDp7vEYThbRN9FOlp3/lBzEyDUv8rGKaHrukYhkE2m2WgUuHJ1z3BU697guXaKn/0pT9jfnmZjGkRxSG2maGvXE4a3ekaxUIRXdcwg4DhvgEsy6K/0sfs4gLDA4MM9PcTpd0AJoZHbhgYR9d1SoXCfT92Qgghto4k+y0WBAHXFuaYX1xkYWUlubrWwDBNLs1ewzJMzl+5RLPdwtCN297vDrw8jZXDhKFNvjRDrjRzfTwcdF1H15Nx7R+ePsrbnnozWdumUipxYGKCOIqoNhvYls1DU9NUShX6y2VOv3KeTqfLwEAfDx2appgvUG82GezrZ6BvAE2DKAopFopMHzz4mkbAE0IIsbtIst9CrufynVMvomsal2eu4YcB80uLjA2NsFKr4nZdrlZnabXaxMSEaWO4jXSao7Rq+9D1gMrQGSy7mVSjhyGWlcEwDOyMxbGpaZ448ShZ2wag0WoxMjDE5Pj+G7fndnn8xEl+4A1vvO/HQQghxO4iyX4LXbh6Bcs0k37qYYBpmERxzOziAjExrW6bIAiwbRvXczesvo9Ck2Z1Cq/bRyZbpdh3Ed0I0NEIwwjLshgbGibWoFIo8X0nH8NfN8tcJmOx0fQCGhqmIadbCCH2Ivn230KdTmetWt4yTKI4wtQNZhfnIYZao47nBViWiWnoeMH3kr2ORtct0lo5TBgZFCqXyRUXsEwTO5MnCpPBcfLZHLquY5om+yfGQdMopCPaAeSzOQq5HEEQrFXBB2FIpVSSed+FEGKPkmS/RcIwpN5s0u50qJRKeIHPV5/5JvVm85Z1A/d7V+IaEMXQau7Dbe7Dtn2Gx69gWl2y9hCFfA4dnUKhAHFMkI5xPz48ysHRCYhj9o/eON3syaMP8fLFC9QaDTRNY6BSYfrg1H0+AkIIIXYrSfZboNPt8tzpU8RxxGqtyvnLF/n2Sy/g+/6d3xxn8RrH0KN+yv1LHD3aYXLfMab2H+Abz32Ha3MzRHGM3WoyOjjMxMgwoJHPZamUyhw9dBjTvPE0mqbJiSPH7k+wQgghHjiS7LfA2UuvYJkmmVKZKI75+rN/SXSbMe2vczt9RJ3jZDI2A2OzZPOr2HY/x6YO8ZVnvkmtXiOXTjvb7nS5Nj9LuVTkycdeD0DHc7FMOYVCCCFuT/pVbYH2unv1ruehaTrabca1jyOdVnWKdu0hcvmY4bGXyReq5LM2hVye1Vqd5dVVYjQ0jXTku5goDFlcWaXRaqFpGq7r0mq3tylKIYQQDyq5LNwChvG9xB6GITk7S7vTJojCW9YN/CyNlWm0qED/4CpvflOFRrcfLYZYg4MT+/EDjyiOuT5L/PXOebEGcRQShsl2NU0jCG/dhxBCCLGeJPt7dO7SKzz/0kss1VbJ2TZjw6MUCwVWaqu3rNttDdOsHkDXQ44cXebkw0OYhkmj00TXdcqlMgfGxnA9j2IuT9dziWPQdA0N0nWKlIpFAAzDoJw+FkIIITYjyf4eXJuf46vf/kvKpRIx0Oq0eOncy9i2TV+5QrPdptVuEUYGrdUpjHiMpx4f5Sd/Yorj0wcAuHjtKucuX8TzXIb6B/ECH9u2eecPPc0ffeWLLCwvARpo0Ffq49ihw4RhcnX/0KHDMsKdEEKIO5Jkv4Gl1RVeOHOaWrPB8MAg+WyWpWqVcqnAY0cfJptOAPPC2dNkM8nIdZVSiUqpRBzHeIHPEycfw+26+F6BuWtj5CYrPPaoxdNvnWL9CLlHJqc4MjmF5/us1KpkM/batLE/+OT386df/RKLK8scnTrEkQNTNDpt4ihioNKHYRgblF4IIYS4kST7dYIw5PNf+ypf+ubXqDUbeJ5H1/MwDYNCoUDezvKFr32Nv/G2H+bk0Yfoui76TWPba4DvB5i6SbU5wMJcH6YVMDk9x+h4P5sNhZ9JR8ZbL5/P8+Ynbhze1k6HxRVCCCHultQBr3P6/FmeO/Mi7W43HXA2uVceRhGBHxBHMa7v8vmvfQXXc5kYGcPboC99pTjI+bODzM/2Ua60OTh9mbERi1KxtM0RCSGEEJLsb3Dp2jVq9TqGrhNGEVEUEqdN4f0gICLGD0I6rseVuVkee+hhioU87W6HKI7odDsU7Cmi1vfjdXIMj88xMnGNUt6mUCgwObH/9gUQQggh7gOpxl/HtCz0tf7x36tv1wB93ex0OjG2lcHOZPjxtz3NK1evcGVuDr81SeQdpHgwYmzfLI2OhW4YTAwNM9DXL43phBBC7AhJ9uscP3yIMxfOceHyJUzDIDRNtCAgimOsjI2u6VimSblUZt/oGJB0hxsfmmJx5iierzE5GXLsWICuDwPDt9+hEEIIsQ0k2a8zOXGAt77p++m6LlfnZgDIZDLohkHGtLAzFuMjw/ytp9+xdpU+M6Nz+rSFpsHjj/uMjNw6ba0QQgixkyTZ3+TR4w/zyLGHaLZa+L5PLptFM3Q818M0TQr5PABBAKdPm8zMGPT3x5w86bFuplkhhBBi15Bkv4FkpLobW87n7Oza43pd4/nnLTodjenpgMOHw0271AkhhBA7TZL9q3TpksHZsyaZDLzhDR4DA/Gd3ySEEELsIEn2d8nz4IUXLJaWdEZGIk6c8MlkdrpUQgghxJ1Jsr8LKys6L7xg4Xlw/HjAwYMy05wQQogHx65N9o7jvB34CGAAv6GU+vB2lyGO4fx5kwsXDAqFmCef9CmVpNpeCCHEg2VXjvLiOI4B/DrwI8AJ4N2O45zYzjJ0Ohrf+laGCxcM9u0LeeopTxK9EEKIB9JuvbJ/E3BOKXUBwHGcTwPvBE5tx86XlnS++12LOIZHH/UZH5e+80IIIR5cuzXZ7wOurHt+FXhyu3ZuWTHFYszJkz65nFzNCyGEeLDt1mS/Ua/1W7Ku4zjvA94HoJRiaGhoS3Y+NATT01uyqXtimuaWxfQgkbj3lr0aN+zd2CXuHdj3juz1zq4CB9Y93w/M3LySUupjwMfSp/HS0tI2FG37DA0N0Wsx3Q2Je2/Zq3HD3o1d4t4aExMTd73ubk323wKOOo5zCLgGvAv42ztbJCGEEOLBtCtb4yulAuADwB8DLyWL1Is7WyohhBDiwbRbr+xRSn0W+OxOl0MIIYR40O3KK3shhBBCbB1J9kIIIUSPk2QvhBBC9DhJ9kIIIUSPk2QvhBBC9DhJ9kIIIUSPk2QvhBBC9DhJ9kIIIUSPk2QvhBBC9DgtjntmCteeCUQIIYS4SxvNEnuLXrqy13rtz3GcZ3a6DBK3xC1xS+wS966O+670UrIXQgghxAYk2QshhBA9TpL97vaxnS7ADpG495a9Gjfs3dgl7m3WSw30hBBCCLEBubIXQgghepy50wXYaxzHuQg0gBAIlFLf5zjOAPA7wBRwEXCUUquO42jAR4AfBdrAzymlvp1u5z3AP083+6tKqU9uZxx3w3GcjwPvABaUUifTZVsWq+M4bwA+AeSAzwL/SCm141VVm8T9IeDvAYvpav9MKfXZ9LVfAd5L8pn4RaXUH6fL305yTAzgN5RSH06XHwI+DQwA3wZ+RinlbU90G3Mc5wDwKWAMiICPKaU+skfO92axf4jePudZ4MuATZJLflcp9cHNyuo4jk1ynN4ALAM/rZS6mG7rVR2PnXSbuD8B/FWglq76c0qpZ3fLZ12u7HfGW5VSjyulvi99/svA55VSR4HPp88BfgQ4mv69D/gorCXMDwJPAm8CPug4Tv82lv9ufQJ4+03LtjLWj6brXn/fzfvaKZ9g47L8WnreH1/3pX8CeBfwSPqe/+o4juE4jgH8OslxOQG8O10X4N+m2zoKrJJ8Se60APgnSqmHgaeA96fl3Qvne7PYobfPuQu8TSn1OuBx4O2O4zzF5mV9L7CqlDoC/Fq63ms9Hjtps7gB/um68/1sumxXfNYl2e8O7wSuX5l/EviJdcs/pZSKlVJ/AfQ5jjMO/HXgT5VSK0qpVeBP2T1ffGuUUl8GVm5avCWxpq+VlVJfT3/xfmrdtnbUJnFv5p3Ap5VSrlLqFeAcyX/8NwHnlFIX0iu4TwPvTK8S3gb8bvr+9cdwxyilZq9frSilGsBLwD72xvneLPbN9Mo5j5VSzfSplf7FbF7W9Z+F3wV+KI3tVR2P+xzWHd0m7s3sis+6JPvtFwN/4jjOM47jvC9dNqqUmoXkiwMYSZfvA66se+/VdNlmyx8EWxXrvvTxzct3sw84jvO84zgfX/cL/tXGPQhUlVLBTct3DcdxpoDXA99gj53vm2KHHj/n6RX4s8ACSbI6z+ZlXYsvfb1GEtsD9z13c9xKqevn+1+n5/vX0tsWsEs+65Lst9+blVJPkFTtvN9xnL9ym3U3Gh0pvs3yB9mrjfVBOwYfBaZJqv1mgX+fLu+puB3HKQL/C/jHSqn6bVbtqbhhw9h7/pwrpUKl1OPAfpIr8Yc3WO16WXs2bsdxTgK/AhwH3kjSXuGX0tV3RdyS7LeZUmom/XcB+D2S/yDzadUN6b8L6epXgQPr3r4fmLnN8gfBVsV6NX188/JdSSk1n35BRMB/Iznv8OrjXiKpBjRvWr7jHMexSJLd/1BK/e908Z443xvFvhfO+XVKqSrwRZI2C5uVdS2+9PUKye2uB/Z7bl3cb09v58RKKRf477z2831fPuuS7LeR4zgFx3FK1x8DTwMvAJ8B3pOu9h7g99PHnwF+1nEcLW0AUkurQv8YeNpxnP60avDpdNmDYEtiTV9rOI7zVHrf72fXbWvXuZ7wUn+T5LxDEve7HMex01bMR4FvAt8CjjqOc8hxnAxJA6bPpPfwvgD8ZPr+9cdwx6Tn4DeBl5RS/2HdSz1/vjeLfQ+c82HHcfrSxzngh0naK2xW1vWfhZ8E/iyN7VUdj/sf2e1tEvfpdT9qNZJ77OvP945/1qXr3fYaBX7PcRxIjv3/VEr9P8dxvgUox3HeC1wGfipd/7Mk3TXOkXTZ+HkApdSK4zj/iuQ/A8C/VErdbYOwbeM4zm8DbwGGHMe5StLy9MNsXaz/gO91T/mj9G/HbRL3WxzHeZykOu4i8AsASqkXHcdRwCmSVt3vV0qF6XY+QPKFYAAfV0q9mO7il4BPO47zq8B3SBLNTnsz8DPAd9N7mQD/jD1wvtk89nf3+DkfBz6ZtprXAaWU+r+O45xi47L+JvBbjuOcI7mifxe85uOxkzaL+88cxxkmqYZ/Fvj76fq74rMuI+gJIYQQPU6q8YUQQogeJ8leCCGE6HGS7IUQQogeJ8leCCGE6HGS7IUQQogeJ8leCCGE6HHSz14IAYDjOM11T/Mks3uF6fNfIBns5IhS6u+k68ckI+Ltuz4Wejoy2gwwrJTS0mVfJBlZLVi3/S8opX7s/kUjhFhPkr0QAgClVPH6Y8dxLgJ/Vyn1uXXLPrTB26ok8zz8Qfr8R0mmNR2+ab0PKKV+YyvLK4S4e1KNL4S4F79FMpzndT9LMiWnEGIXkSt7IcS9+D/AP7w+Vjjwg8CHgF/dsRIJIW4hyV4IcS+6JFX4P00yJvhn0mU3+0+O4/y7dc//s1LqX2xD+YQQSLIXQty7TwH/hiTZ/9Im6/yi3LMXYufIPXshxL36CslMYKPAV3e4LEKIDUiyF0Lck3RO8h8Dfjx9LITYZaQaXwhxz+5invH/4jjOf1z3/IxS6g33s0xCiO+R+eyFEEKIHifV+EIIIUSPk2QvhBBC9DhJ9kIIIUSPk2QvhBBC9DhJ9kIIIUSPk2QvhBBC9DhJ9kIIIUSPk2QvhBBC9DhJ9kIIIUSP+//vSW4EE8WapwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 576x432 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"c=['#415952', '#f35134', '#243AB5']\n",
"tm=np.array(df[:N].sort_values(by=['TIME'])['TIME'].tolist())\n",
"inx=sm_fr[mask].index.tolist()\n",
"carray=np.array([c[0]]*N)\n",
"carray[inx]='red'\n",
"\n",
"fig, ax = plt.subplots(figsize=(8,6))\n",
"plt.scatter(df['TIME'][:N], df['HR'][:N], c=carray, alpha=0.3)\n",
"plt.plot(tm, (lambda x: x*m.params[1] + m.params[0])(tm), c='b', alpha=0.3)\n",
"plt.xlabel('{}'.format(model.exog_names[1]))\n",
"plt.ylabel('{}'.format(model.endog_names))"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"t = 2*np.sqrt((p+1)/(N-p-1))\n",
"mask=sm_fr.dffits > t"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, 'HR')"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfsAAAF6CAYAAAAeZ/GvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzs3XmMnOl94Pfve1VV19VX9U32wW4ezWs4M9RwJI/WksaWPLKxMrLr106CjWwrVoDsbhYxgtgGgiAIAqwXSNbxH7sGBCuBlM1GfteJYyVxNvZKliVLM6O5h/fRTTb7vqq76656jyd/vMVmN8khu8lmn78PQLDqred963mqm/zV+xy/R1NKIYQQQoj9S9/pCgghhBDi+ZJgL4QQQuxzEuyFEEKIfU6CvRBCCLHPSbAXQggh9jkJ9kIIIcQ+J8FeCCGE2Ock2AshhBD7nAR7IYQQYp+TYC+EEELsc+ZOV2ALSd5fIYQQB422kUL7KdgzNTW101XYUplMhoWFhZ2uxraTdh8sB7XdcHDbLu3eGt3d3RsuK934QgghxD4nwV4IIYTY5yTYCyGEEPucBHshhBBin5NgL4QQQuxzEuyFEEKIfU6CvRBCCLHPbds6e9u2/3PgPyZMfnMR+A2gC/gO0AK8D/wDx3Fqtm1HgW8DLwOLwK86jnNnu+oqhBBC7Cfbcmdv23YP8J8B5x3HOQ0YwK8B/wz4A8dxjgJLwNfqp3wNWHIcZwj4g3o5IYQQQjyF7ezGN4EG27ZNIA5MA18A/rT++reAX64//kr9OfXXX7dte0MpAYUQQgix3rYEe8dxJoH/HrhLGORXgPeAZcdxvHqxCaCn/rgHGK+f69XLt25HXYUQQoj9ZlvG7G3bbia8Wx8AloF/A7zxiKL3NrN51F38Qxvd2Lb9deDrAI7jkMlktqS+u4VpmvuuTRsh7T5YDmq74eC2/SC2u1aDbHbn2r1dE/R+DrjtOM48gG3b/wfwGaDJtm2zfvd+CLi3k80EcBiYqHf7NwLZBy/qOM43gG/Un6r9trGCbBZxsEi7D56D2vaD1u5sVuPSpQiJRJqzZxewrK257mY2wtmuYH8XeNW27ThQBl4H3gX+Gvj7hDPyvwr8eb38d+vP36y//n3HcWQLWyGEEHuGUjAyYjA6ahKPK159VeG6O1OX7Rqzf5twot37hMvudMI78t8Bftu27VuEY/LfrJ/yTaC1fvy3gd/djnoKIYQ4OLSFBSI//CHRv/xLrDffRCsWt+za5TK8+26E0VGT7m6fV1+t0di4ZZffNE2pfXPDrGQ/+/1B2n2wHNR2w8Ft+25ot7awQPRHP0IlEuGBIECr1ah88YsQiz3TtefmdC5ftlAKhodduroC4LntZ7+hlWqSQU8IIcSBY125cj/QA+g6yjQxr19/6mv6Ply9avLhh1a92762Guh32rZl0BNCCCF2C61affigaaI/ZVd+oaBx8aJFPq/R3+8zNOSh76LbaQn2QgghDpwgmURfWWFdRK5WCfr6Nn2tyUmDa9dMDANeesklk9kdd/Nr7aLvHUIIIcT28M6eRavVwKvndatWIRLBGxzc8DVcFz7+2OLyZZOmpoBXX63uykAPcmcvhBDiAFKJBJUvfQnz2jW0YhHV1xcGenNjYXF5Oey2r1Q0jh716O/30XZxUncJ9kIIIQ6maBTvhRc2dYpScOeOwa1bJrGY4lOfqtHUtPtXtUmwF0IIITagUoHLly0WF3U6OwOGh90ty4b3vEmwF0IIIZ5gYUHn0iUL34dTpzx6evydrtKmSLAXQgghPkEQwK1bJnfuGKRSijNnXJLJ3d9t/yAJ9kIIITZNHx/HunYNrVRCJRK4p04RdHXtdLW2VKkUTsJbWdE4fNjn2DEPw9jpWj0dWXonhBBiU7T5eSLvvgtKoRoaIAiIvPUW2vLyTldty0xP67z1VoRSSePcOZfh4b0b6EGCvRBCiE2yrl1bn2oWUA0NmFev7lCNto7nwaVLJhcvWqRSildfrdLevjvXzm+GdOMLIYTYnEft06rraDu1f+sWyeU0Pv7YolTSOHLEY3Bwd6+d3wy5sxdCCLEpQWvrwwG/WiXo6NiZCm2BsTGDn/40QhBonD9fY2ho/wR6kDt7IYQQm+SdPIkxO4tWLqNiMbRKBZVK4Q0N7XTVNq1WC9fOz8/rtLUFnDrlEonsdK22ngR7IYQQm2NZVH/u59CnptDn5wk6OsKZ+HvsVjib1bh0KUKtBidOePT27q2185shwV4IIcTm6TrBoUMEhw7tdE02TSkYGTEYHTVJJBQXLrikUntv7fxmSLAXQghxYJTLcOlShKUlje5unxMnvI3ufbOnHYAmCiGE2FOCAP3uXYypKVRDA/6JE+F6/mc0O6tz5YqFUnDmjEtX195fUrdREuyFEGI3UAp9dhZ9epqguZmgtxf0A7hgSikif/u36MvLqFgMVlawxsaofPazqNbWp7qk78ONGybj4waNjWHK23h8f3fbP0iCvRBC7LQgCANcNhsmpxkfR127RvULX2BfTg1/DH1mJvwc4vHwgGEQxONYH31E7Qtf2PT1CoUw5W0+r9Hf7zM05B3I71AHsMlCCLG76Ldvh3ey8ThoWnhHqxTWxx/vdNW23b2u+wfphcKmrzU5afD22xFqNY2XXnI5duxgBnqQO3shhNhxxuxsGODXHTTQl5Z2pkLPolLB+uAD9JUVlGniHT1K0Ne34dP9pib0yUl44PN46PN5DNeFq1ctZmZ0WlsDTp92iUY3fPq+JMFeCCF2WjQKudxDY/TqWaaJl0pYly+j5fOoZBLv1KmH8tlvOd8n+r3vgWGE6XODgMiHH+IGAf7AwIYuEfT3w/Xr4d6y9c9DK5WonTu3ofOXl8OUt9WqxtGjHv39uycTnlI7N0/ggHZoCCHE7uGdOIFWLq87ppVKeMeOPd0FKxVi/+7foS8uonkeejZL9K/+Cq1Y3ILafjL9zp11QRrqG+TcvLnxixgG1S98gaClJdxVzzSpnT8ffgl4DKVgdNTgnXciaBp86lM1BgZ2PtArpRgdH+Ptjz/ghz99i4+vX6VcrWx7PeTOXgghdphKJKi+9hrWxx+jFwqohgZq584R9PQ81fWsq1dRlnU/6Oo6KhbDvHgR99VXt7Dm6xnLy4+cUKhVNhncYjHcV17ZcPFKBS5dsshmdTo7A4aHXSxrc2/5vIyO32VhKYtlmkQjEVbKZS5eu8r5My+gb+MEAgn2QgixC6i2Nmqvv74l19JyOR7afF3X0UulLbn+J/F7ejDGx+/PpK8Lksnn9p4LCzqXLln4Ppw65dHTs3tS3gZBwHx2gYh1/wuQpmn4QUB2eYlMy9MtJXwa0o0vhBD7jEqnw43Z1wqC5xp0gTBHfmvr/SEJz0MrFnE3ON6+qfcK4Pp1k/fft4hGFa++WttVgR7CLvzgEeP0hq5TqdW2tS7bcmdv2/Zx4E/WHDoC/NfAt+vH+4E7gO04zpJt2xrwh8CXgRLw647jvL8ddRVCiL3OHR4mNjGBAjBN8H20Wg3v9Onn+8aaRu2119AnJzEmJlANDeG8gy3IfrdWqRROwsvlNA4f9jl+fHcuqTMMg4ZojCAI0NZMHnB9j0xzy7bWZVuCveM414FzALZtG8Ak8GfA7wLfcxzn923b/t36898B3gCO1v9cAP6o/rcQQognicWo/PzPY129ipbLoRKJcDb+FgfdR9K057pBztSUzrVrFpoG5865tLfv7pS3R/sGuHjjGnq9+75crdDT0Ulsm9cC7sR3odeBEcdxxoCvAN+qH/8W8Mv1x18Bvu04jnIc5y2gybbtru2vqhBC7FGxGO6LL1L72Z/FPX9+ewL9c+R5cOmSyaVLFqmU4tVXq7s+0AMkEwk+dfYch7t76My08eLwafp7Dm97PXYi2P8a8L/VH3c4jjMNUP+7vX68Bxhfc85E/ZgQQogDJpfTeOutCNPTBkeOeJw/X9vqkYHnyjQMuts7GOzrJ75DFd/W2fi2bUeAvwv83hOKPmpl5EOzHGzb/jrwdQDHcchkMs9cx93ENM1916aNkHYfLAe13XBw276Zdo+OwrVrGskkfPaziqfcC2dX2Mmf93YvvXsDeN9xnNn681nbtrscx5mud9PP1Y9PAGv7OQ4BUw9ezHGcbwDfqD9VCwsLz6naOyOTybDf2rQR0u6D5aC2Gw5u2zfS7loNLl+2mJ/XaWsLOHnSRSnYyx/XVv+8u7u7N1x2u4P9v8/9LnyA7wJfBX6//vefrzn+j2zb/g7hxLyVe939Qggh9rdsNlw7X6vBiRMevb27a0ndXrRtwd627Tjw88B/subw7wOObdtfA+4Cv1I//heEy+5uES69+43tqqcQQoidoRSMjJiMjhokEooLF1xSqYO17/zzsm3B3nGcEtD6wLFFwtn5D5ZVwD/cpqoJIYTYYeUyXLoUYWlJo7vbZ3jYeygJoHh6ki5XCCHEjpqd1blyxUIpOHPGpatr9y+p22sk2AshhNgRvh+mvJ2YMGhsVJw54xKPS7f98yDBXgghxLYrFMKUt4WCRn+/z9DQ7kx5u19IsBdCCLGtxsbg7bcjmCa8/LJLa6t02z9vEuyFEEJsC9eFq1ctSiWNpqaA06ddtjlF/IElwV4IIcRzt7wcdttXqxqf+pSisdFFe1SuVPFcyAiJEEKI50YpGB01eOedCJoGn/pUjaEhJNBvM7mzF0II8VxUKnDpkkU2q9PZGTA87GJZO12rg0mCvRBCiC23sBCmvPV9OHXKo6dHUt7uJAn2QgghtkwQwM2bJmNjBqmU4uxZl0RC1s7vNAn2QgghtkSpFE7Cy+U0ent9jh2TtfO7hQR7IYQQz2xqSufaNQtNg3PnXNrbZe38biLBXgghxFPzPLh2zWRqyqC5WXH6dI2Ghp2ulXiQBHshhBBPJZcLu+3LZY3BQY8jR3xZUrdLSbAXQgixaWNjBjdvmkQi8PLLNVpaZBLebibBXgghxIbVauHa+YUFnfb2gJMnXSKRna6VeBIJ9kIIITYkmw3XztdqcOKER2+vrJ3fKyTYCyGEeCylYGTEZHTUIJFQXLjgkkpJt/1eIsFeCCHEJyqX4dKlCEtLGj09PidOeBjGTtdKbJYEeyGEEI80O6tz5YqFUnDmjEtXl6yd36sk2AshhFjH9+H6dZOJCYPGxjDlbUODdNvvZRLshRBCrCoUwrXzhYLGwIDP4KCkvN0PJNgLIYQAYGLC4No1E8uCl192aW2Vbvv9QoK9EEIccK4LV65YzM7qtLYGnD7tEo3udK3EVpJgL4QQB9jSksbFixbVqsaxYx59fZLydj+SYC+EEAeQUnD7tsHIiEkspnjllRqNjTIJb7+SYC+EEAdMpRKmvM1mdVKNeQYHSySTKUAW0O9XEuyFEGKLaOUy+D4qmdzpqnyi+Xmdy5ctyhUXIpcIrDy37sLIuMZQXz9tLa07XUXxHEiwF0KIZ1WtEvnJT9CXl9GAoKGB2oULqObmna7ZqiCAmzdNxsYMUilFLH0FK1JD0+7PxLs1dofmxiZMSZG372xbsLdtuwn4Y+A0oIDfBK4DfwL0A3cA23GcJdu2NeAPgS8DJeDXHcd5f7vqKoQQmxF58020ahWVSHBv1Dvy4x9T/fKX2Q2L1EulcO18LqfR2+tzZLDKuxcL6wI9gKZBdjlLe2vbDtVUPC/b+Vv4h8C/dRznBPACcBX4XeB7juMcBb5Xfw7wBnC0/ufrwB9tYz2FEGLjqlX0paWHgrrmeehzcztUqfumpnTeeitCuaxx7pwb5rbXgU+Yca9rO//lRGy9bfmp2radBv4O8E0Ax3FqjuMsA18BvlUv9i3gl+uPvwJ823Ec5TjOW0CTbdtd21FXIYTYFPWYGezBziWl8Ty4dMnk0iWLVCrg05+u0t4e1scwDBpTafwH66dptDTtnqEHsXW2qxv/CDAP/M+2bb8AvAf8E6DDcZxpAMdxpm3bbq+X7wHG15w/UT82vU31FUKIjYnFCNJpNM9j3QJ1XSfo6NiRKuVyYbd9uawxOOhx5MjDa+ePDwxy7dZNcsUCCohFIpwaOoa+C4YdxNbbrmBvAi8B/9hxnLdt2/5D7nfZP8qjOpge+vps2/bXCbv5cRyHTCazFXXdNUzT3Hdt2ghp98GyL9r9i7+I9v3vQ6kUPrcs1Oc/T/IJwf55tH10FK5d00gm4bOfVbQ+ZnJ9Z0cHtVqNQAXEorEtrcfj7Iuf+VPYyXZvV7CfACYcx3m7/vxPCYP9rG3bXfW7+i5gbk35w2vOPwRMPXhRx3G+AXyj/lQtLCw8l8rvlEwmw35r00ZIuw+WfdPuT38abWkJzfcJWlvDMfwntGsr216rhWvnFxZ02tsDTp50UeqJVVhVyBe2pB4bsW9+5pu01e3u7u7ecNlt6a9xHGcGGLdt+3j90OvAFeC7wFfrx74K/Hn98XeB/8i2bc227VeBlXvd/UIIsStpGqqlhaCtbdtn4GezOm+9FSWb1Rke9jh3ziUS2dYqiF1uO9fZ/2Pgf7VtOwKMAr9B+GXDsW37a8Bd4FfqZf+CcNndLcKld7+xjfUUQuwEpdAnJzGmplDJJN7gILIby+MFAYyOmoyOGiQSigsXXFIpSXkrHrZtwd5xnA+B84946fVHlFXAP3zulRJC7A5KEfnRj9CXllANDTA/j3nrFtXPfQ6VTu907XalcjncwGZ5WaOnxw+X1EkuHPEJZNqlEGLH6VNT6NlsGOghnOAWjWJ9+OHOVmyXmp0N184XChpnz7qcOiWBXjyepMsVQuw4Y2oKFY+vP6hp6IXtmzS2F/g+XL9uMjFh0NioOHvWpaFBuu3Fk0mwF0LsOJVMwuwsD84qUzLLbFWhEK6dLxQ0BgZ8Bge93ZCJV+wR8qsihNhx3tAQmlLrMs5p5TLu0aM7WKvdY2LC4K23Iriuxssvuxw9KoFebI78ugghdp5lUf3CFwiSyTATna5Te/FFgr6+na7ZjnJd+OgjiytXTJqbA159tUpr686l4BV7l3TjCyF2BZVI4L72Gu5OV2SzlAqjsmlu6fr6paVwtn21qnHsmEdf38Mpb4XYKAn2QgjxlPSJCayLF9EqFbAsvCNH8E6efKZrKgW3bxuMjJjEYopXXqnR2CiT8MSzkWAvhBBPQcvliPz0p6hkEpVIAGDeuoWKRvEHB5/qmpVKmPI2m9Xp6goYHnYx5X9psQXk10gIIZ6CeePGQ8sFVSyGefv2UwX7+Xmdy5ctfB9OnXLp6ZGxebF1JNgLIcTTcN1Hj9Fvcg/7IICbN03GxgxSqXDtfCIh3fZia8lsfCGEeAp+by/avS1t7/E8/La2DV+jWISf/jTC2JhBb6/PhQs1CfTiuZA7eyGEeApBdzf+oUMY4+Moy0JzXVQ6jXfmzIbOn5rSmZzUKJc1zp1zaW+Xbnvx/EiwF0KIp6FpuOfP4x0/jj43R9DYiMpknnia58G1ayZTUwYDA3DmTJVYbBvqKw40CfZCiO2nFFouB5b1+HJBEJaLRu9vkvMk5TJatYpKpXjs7jBKod+9i1Yu4w8NsTrt3fPQCgX8SIQSikq1iqZpJBri6LqO67rEolH8IAgfx+P4vYdZyefxF+aJWBaxaJR47H59gyCgXK1QKphcv56kXNYYHPR45RXF4uK96ijKlQqaptEg0V9sMQn2Qohtpc/OYr37Llq1GgbWvj44d+6hvev1iQkiH34I9c1wgu5uap/+9Cd/QXBdrDffxLgXPU0T99Qp/CNHHiqqTU2R/qf/FGN+nkApiMcp/OZvojo6sK5dY7FY4P3lLDfdGnMRC8MwiUaiZFqa6WnvZGF5iWgkSlMqxeTsNNPzcxSKJYqlIrFojL6eQ5wcOsqpo8ep1qrcGBtjatJicjxFIm7yS7/QQXubsZokp1AscnX0FjW3hgbEolFODh0n9sBn4nkemqZhyBZ3YpMk2Ashto/rEnnrLQgCzFu30CoVtLEx4h99ROlrX7t/J16pEPnhD4m8/TZGNgtAkEqhKhXcL35x/TWDAOuDD7B+8hP0+XmIx/EGB1GJBNaHHxK0tYV3+WvKp/75P0crlwmamsJjSpH8l/+S6htvUGtv53K2zJTrEiwv09zUSDYC5VwZTdeo1GrEYzEKpSJjk+NMz8+xuLIMKkDXdHKlApVa2BsQKEWlqpgeb2NuATw1TTo+xd+8Z/HC8WFaWlpQSnFl5AamYdIQjdWro7g6cpMXT54GoFgqcf3OKJX6nX9TYyPH+4+gS4J8sUHymyKE2Db61BQqCLCuXAHCdenEYhiTk1jvvLNazhwbI/KjH6GXSqhEApVIoHkeie9856GlbebFi+jT0+F2uMkk6DrWtWvgeah4HOPWrfV1mJ7GmJlhXbaa+na6xp07zJVLFFyPmlsFw8AslKh53mq3/dLKErqusbyyQnZlmXypAEFA1XXxVYCm6RTLJSZnp/n4ygyXPmpmZdnEiN6kt3+eeNzCdV2yKytcGx0hVyjgut66OmqaRqlSplKtEgQBl25eQwUBsWiUaCRCvlDgxp3bW/zTEfuZBHshxPZRCj2bRakHlpcZBsbU1Gog17JZjFxu/Zi7rkO5jD4+vu5Uc3Ly4a1xDQN9Zia81oPvpRSfuLhtzReJR86NV+GfarWGH/goFRB4QRjkARUoUOArjYW5Vq5dbcYwAg4NjNHall93Kcs0WcguoNSjZ+HfS4O/uLyEUgptTWJ80zBYXlkh2OSafnFwSbAXQmyboKcH7cFkNEFAkEqFx+qB2e/qAt9/4OQAGhrCsf61vPCuWDU3rz5G19FcF61UwhsaWn+Z7m6Ctrb7ZQGUIkgk8AcGaGuIk7BMolYEAh8vESdimui6jmVZNKebcH0PXwU0ptOYloFSoBsGoPA8k+XZIWrlQ3R1e5w+u4wVqaFper0ZAfGGhrCbP1Ak4wmsB3LiKqWIRaPEolFcz1s9d107VBDONxBiAyTYCyG2j2VR+dKX0Mrl8E+phDIMvKEhgsbG1Tv5oLcXr7c3DPjVKlqtBoaB39ODPzCw7pJ+SwsEAV5/f9jdX6mgrawQNDXhnj2LSqfX10HXyf/2b6OiUbSlJbRsFlWpUPit38J95RUi1SonIzG6DBO9sZGlWBTTMGhKp2hqTDHY20ciliAejdGSbiKdSBNviKGjUcw3sjh9AkNLc3K4whuvt3FicIBEQwOlSpmqW8M0DHraO1FKkUomMU2T4aGjeL5HuVqhXK2glGJ46BgAmeYW/Ae++CiliDc0YMpEPbFB2kPdaXuXmpqa2uk6bKlMJsPCwsJOV2PbSbv3P+PmTSJvv42KREi1tpKrVKh97nOrG8oAGNevE/9X/yrcHcYw0CyL0t/7e3jnzq2/WKVC9G/+Bq1cDrvvSyW8vj7c8+d57C4ySoXj+ZUK/rFj91cDuC5aPk8QjZJXAZVqFdBIJRLouk7NdWmIxZiZn+PGnRFUoLgzOcO1qzqVYhuNjT4vngsYHOjkzLETGIZBEATcunuHqdlZGqJRfBVgGSaf+8xrFOurDZRSlMpl0CDRsD7n/tTcLHcm7qLrBoEKMA2TM8dO7Nklegfpd32trW53d3c33B/xeSwJ9ruY/IM4WA5cu6tV9KkpWrq7WbCsR+eZr1YxP/4YfB/vhRfgk9baK4U+NweFAkF39yeX22Ke5zE+tcLVqwmKBZ2Ozhz9Ay5N6STpZOqh8jXXJbuyTNSyaEylaW9v3/DP3PVcFpeXsUyL5nR6T8/EP3C/63U7Gexl6Z0QYmdEowQDA5DJwCf9BxiN4n3qU0++lqYRdHRAR8fW1vEJZmaijNzoJhGDC+ddWlsjjy0fsSw6MxvPnb+WZT79uUJIsBdCiE1yXbhyxWJ2Vqe1NeD0affBnEBC7CoS7IUQYhOWljQuXrSoVjWOHfPo6/PRNtSRKsTOkWAvhNg4pdDHxzEmJyEWwztxYuM563e5UrnMxOwMvu/TkcnQ0ti07nWl4PZtg5ERk1hM8corNRob982cJ7HPSbAXQmxY5Mc/Rl9cDAP88jLRO3eofvazG9rtbTvNLS4yNTeD53sk4gkGD/cRecymO/PZRW7cGSUWiaJpGtdGbtHe2spQX7jMr1KBixctlpZ0uroChofdx070F2K32bvTOYUQ20qfn0dfWLh/J28YYf75jz7a2Yo9YGZ+npG7dwiCMFd9sVjko6uXPzHbnFKKO5PjNERjq1nqYtEos4sLVKpV5ud13norSi6nc+qUy5kzEujF3rNtv7K2bd8B8oAPeI7jnLdtuwX4E6AfuAPYjuMs2batAX8IfBkoAb/uOM7721VXIcTD9KmpMJf9g8fz+UeU3jmTs9NE16TP1XUd13OZW5yns+3h2fqu5+G6LkZ0fYIaHYP3PqhSzKVJpRRnz7okEtJtL/am7b6z/7zjOOccxzlff/67wPccxzkKfK/+HOAN4Gj9z9eBP9rmegohHhA0NT2cqhZQu2wauut5Dx0zDZN8qfzI8qZhPLRlbKVscPlSEwtzSXp7fS5cqEmgF3vaTnfjfwX4Vv3xt4BfXnP8247jKMdx3gKabNvu2okKCiFCweHDYWBfk7pVK5fxjh3bwVo9LPaI3odarUamqekRpcM7/45MGzXXBWB+LsZHHzaia3E+82mTEye8R+b7EWIv2c5fYQX8pW3b79m2/fX6sQ7HcaYB6n+314/3AGu3tpqoHxNC7BRdp/r5zxO0toY7x+k6tXPn8AcHd7pm6wwe7qVSreLXx+irtRqpVIqmdOMnntPfc5iutm5GbqUZHUnR1RnjV/+9NtraZFc5sT9s5zSTn3EcZ8q27Xbgr2zbvvaYso9atfpQH1r9S8PXARzHIbPLZgQ/K9M0912bNkLavcv1bO337q1udyaToburmzsTd6m6Lm0trXRm2tZtEfug5WXIZdvo7YTXP6s4epRtWTu/Z37mW0zavQPvvV1v5DjOVP3vOdu2/wx4BZi1bbvLcZzpejf9XL34BHB4zemHgIcS3zuO8w3gG/Wnar/lWpb80QfLgW+3Uuizs1AuE3R1wTNu8tKavt9tv7i4+MgySsHYmMGtWyaRiOLMGZfmZsUnFN9yB/5nfsA8p9z4G7Itwd627QSgO46Trz/+IvDfAt8FvgqSN00QAAAgAElEQVT8fv3vP6+f8l3gH9m2/R3gArByr7tfCLEPlctEf/ADtGoVZRhoH36INzyMd+LEc3vLahUuX7ZYWNBpbw84dcpl7VJ83/cpVSpYpklsl01CFGKztuvOvgP4M9u2773nv3Yc59/atv0O4Ni2/TXgLvAr9fJ/Qbjs7hbh0rvf2KZ6CiF2QOSdd0DTUPFwa1cVjWJeuYJ/6BAqmXyofM11mZ6fxfcDutraN73V6+KizqVLFq4Lw8Mehw+v3y9+am6Wu1MT9XF/jXQyycnBo+tm7X98/So/ef8dFpaXcGsu3R2dHOro4tjAAIO9/YzeHeMHP32TfLFIvCHGay+9wsmjx9A0jfnsIhOzMyTicWq1GoHvcWvsLq7n0pxOc/7suYcy+AnxLGSL211MuroOloPc7sI3v/lw2l3Pw+/txTt9et3h7PIy127fwjJMNE2j6tbo7zlMT0fnE98rCGBkxOT2bYNEIlw7n0qt/z+wVC7zwZVL675AeL5Pc7qRo/1hRr3LN67zf/3grwiCgJV8Hs/zCQg4eeQoTelGWhqb+ODqRSwzgq5pBCqgUqvxS597nUxLK6N37xCNREmn03x4+RK3x8c43N2DoRv4vo+u67zxs58nlXj4i85+cJB/13dqi1tZUCKE2HkPrHMHwiV+D3SfK6UYGb9DLBLFMAx0XachGmNsagLvEevr1yqXNd59N8Lt2wY9PT6vvlp7KNADTC/MrUvKA+Fa/KXcyurztz56n4hhsZLPUypX8AMfHY2p+Vlc1+W9Sx+hFOj1WX66phO1LH568SMmpqeIRsJ2eZ7H7MI8pmlRqucBMAyDmutya+zOEz82ITZKgr0QYnM8D61cDme3bdH5fk8P1GrrimlBgNffv+5YtVaj5j4c1DU0lgv3M/l5nkfNrVFza3i+z8yMzo9/YpJd8jh+osjwcA2lPGquS7VWo1ILyz7Ovdn5QRBQLJdZLuSoui5KBfieT7XmUq3WQNMoVys82Gtq6AblcpnKmsREnu/heR6mYeB67rr3KlUenQRIiKchGZ6FEBujFNYHH2CMj0MQoKJR3LNnCQ4dCl+v1TAvX0ZfWkLd2xGvpQVtYQF9ZYWgrQ1jZARzfBx8Pzz/zBnIZHBfeAHL98Pd9IIAFY9Tfe01eGDzGtM0H7kkLlCKqGnh+T7XRm4yu7jA9PwcKtDxKoPMz8aINVRJNI0w914N3teImCY112U5nyfeEKMp1URfdzd9PYeYmZ+jIbqmG9/zaKkn5dF1HdPQUQos3cDXXDRNw/e8+pi+ormxmWptfbbBmudyuK2NZCKB53lomkYkEiUai1IslkinU/c/aqC1qXlLfmxCgAR7IcQGmZcvo09Orhtbj/z0p1RbW1GWRfT73w8HxU0TLZ8Pn+s66DrKMDDGxtA8D+/kyfvnv/suDA2BruOeP4/74ovgeQ9136/WwTBoaWxmJZ/Dqu9GEwQB8ViUVDLJ5ZvXKZTLzC4soLwEd+9kWMjm6OqaRotNkk61MbdYpFQpU3NdIpaFZVoUiiUSsThz2QU836e/5zDj01P4gU+lUsE0Lbrb788JOH3sBNPz8/WhBINqrYZlmqQSSeINcV4/Psz33/4J5WqFWCRKzXUxDZPXXn6FxmSKj69fwTRMdE3jUHsXt+7exjJMFIpKpUKmuYUjh/ue009SHEQS7IUQG2JMTDw8ht7QgHH9OiqVAte9/7qmoS8soOVyeOfOAaCXSmjlMtriIqq1dfV8rl6FI0fqb2I8evx+jWP9A9yeGGdxeYkgCGhMpTna14/n+6wU8pQrZZYWU8zPZjCMgJb2UWJJF8/38YKwTC6fJ6h3vzc1NpJoaMDzffKFIq1NzRiGwamjx/irH/8N2ZUchq4zNjXOycFjnD/zAm2tGb702c9x9dYNsivLGLpOOpGio62dz7/6adLJFAO9/fzkvfeYnp+mraWVz7x0nnQyvHt/+fQLTM5ME4s38Mq5c/zMy+e5MnKLUrnI4WPDDPUPYD7hcxBiMyTYCyE2RAsCVD1JvL6wgDEzA0GAVqngd3U99EVAz+XCu/R7fD/8cjA/j1cP9mha+CVhE3RdZ7C3j8He9Xe+rufiuXDzRjMzkx7ppgqHDi8yNV8kCHQUGm6tRqlcQtd1CCDQg9Ux+NYgIFCKWs3l4vWrjM9MsLi0QryhgUxTC7quc+nmdbo7uujt7GZpZZmXTp1B0zSUUtRcl3PDp1Zn8cdjMX7uZ37mkW2IWBYDh3vXzc7+mZfOP7KsEFtBJugJITbEz2TA8zAmJjBHRsKJdsUiVCrok5NQqaw/QdPW7YinmprA81BrNoPXSiU4enRL6lfIRxi50Y1fa6Ktc4G+gXlMKyBiRWhMpzF0nUqtimVY6LqOYRpErHDWvet6RCMWDbEoY1MTpBNJ8oUysUgU3/dZWM4C0BCLcXXkJol4nLPHholGoyiliEajnD0+vOn1/kJsF7mzF0JsiPvii0T/+q8xb98miEahUsHv6UG1tqLl86ggCO/SLQuCgCCZRKXTq+d7AwNYCwsE6TSUy+Fs+8FBaG+HZ1h7rBSMjhqMjpr0dnejxy5RqmksLFUAxZHePgrFIslEkunZWSJRC903sAwDDcgXS0QjEXylaIjFiFpRrDUTAzU0ajWXoL6xjla/RUomEpw+evyp6y3EdpJgL4TYGMui9tpr6DMzKNNENTauzpZXpol79ix6sYi+uIiKxah95jPregFULEbp138dlUigFQoEbW3wYCKdTapU4OJFi6Ulna6ugOFh0PWT5Aq99SVuYQBPJZLki0UWslmujt6gKd1I4Ae4nofre6QSSY71D+CrgNt37wKQaW5idmEBy7QAhQIqlTInByXAi71Hgr0QYsNULIZqbl7XPQ/18fy2NryBgXXHvZMnw/z2rguRyOpiddXS8sx1mZ/XuXzZwvfh1CmXnp5729HqNK3pUbinMZWiMZUimYhze2IcX/MxTZOudJoTR4bQdR3f91eD/bH+QSrVGosrSwC4vsu5k6fpPIC7tYm9T4K9EHvZvXHySAStUkFFImA+x3/Wuo57/DjW5cuoRAIArVwOu/PX5LAPgmB1aZuu65+4lO5JfN/H9TwilhWmxq3V0HWT0ZEod+8apNPhTnWJxMYT/HRk2mhraaVcDTe5uTduD2H2uqP9A9y8cxuA4cGj+L5PX3cPne0dRB5Y9y/EXiHBXog9SCuVsN58Ez2XQ5+fR89m8fv7UfE4Xm8v3gsvPLcN2f1jx1DpdNg97/u4g4ME95bOAZOzM0zMTOP5Poau09XWTl/PoU29h1KK2xN3mVtYxFcBlWoZPwgwSHFnpAVDa+LC+RaOHw/Qn2Kasa7rJBrij3wt09xCc7qRpXwOUzdIJ5PhFxYh9jAJ9kLsQdbf/i1afcKYMT2NisXQx8bwXnwR8+5diMfxjh17bu8fdHZS63x445nl3ApjkxPEotHVpDdTczPEolE6Mm2r5VzPZWJ6mlK1Qm+5REMkum5d+dTsDHOLC0QiETzfY2xinqVsElUbIBY16e6bJJpcQdefT+IZwzDISAY7sY/I11Uh9hgtn0cvFkHXMaamwvFzTUOvVNDKZVQshlEfd95uU3OzD20iE41EmV28P9u+Uq3y3qWLLCxlqVarTMxM8eGVS3j+/W1m57KLq93rc4vLTE92MjvdiRcscvrsEplMQHZ5eXsaJcQ+IMFeiL3G9+9vIhME67vr1x7fAeVKlezKMoVSsT5//V617j++PTmOZZqre8NbpkUQBEzOTq+WCYKwfKFgcv1KO7mVBB2dSwwMzRGJhG0L1M60UYi9SLrxhdhjVGPj6mz4oKMD8/p1VDSKsqww/azrhrvIbbORu2PMLi6wsJTFNAxM02Cg5zCBUmSaW1fLVSoVtAfmE5imSaFYXH3emE5z+UqF6clGkvEK6ZbLxBM1mhu7APCDgHTy4Rn3a3m+T6FYJBaNEnvKCYJC7BcS7IXYazSN2oULRN98E2VZBOk0+sIC3vHj6KUSfmsr3pkz21qllXye2cV5OlpbCXyPlUIe1w24PTXOqaHjHFozvm9ZFtVqdV3A94NgNSBXq7AyP8jczDTxxDIDg3mW8xaeZ2CZFuVqlXgsxlBf/yfWZ2J6mvGZqdW7/3QyxamhYzLRThxYEuyF2INUJkPlF38RfW4OPvMZgmQy3Fo2lQqT3TwnQRA8MmDOLc4TtSJomkZ3R2e4tK1SwTBNTg2tnyjY332Ij69fJRoJywdK4Xkeh7u6WVzUuXTJwnXhjZ/rpLklQblSJp06jaZprOTzRK0IqTXL/B5UKBYZm55Yt0VtuVLh9sRdBnv7t+yzEGIvkWAvxF6l6wRr7piD+rr35yG7vMydibvkSkVQisHD/Rzu7l593TRNAqUw6nfrlmVhWda6sfp7kokEZ46fYGxykppbIxmP09PWydidOLdvGyQSipdeckmlFBAnEb+/RC7T/ORkPNPzc8Qi67vtTcNgKbfylK0XYu+TYC+EeKxSucyVW9eZz2YpVcoAjE1Ocu7kaV44MQxAT0cn0/NzGGuCbM116W7veOQ1U4kkp4+FaWfj8Qw/+EGO5WWNnh6fEye8J+1y+1iaFk4IfHBeADyfvANC7AUygCWEeKyJ2RmWc3mqrkvEihCxIsQbGrhy8/rqpLqIFWF48ChKKSq1KtmVFfLFAovLS9wev4u/ZlndWjMzOj/8oUaxqHH2rMupU88W6AG62zup1mrrjnm+R2tj07NdWIg9TIK9EOKxPM+lWC5hPDBWrxsG0/Nzq8+b0428fPosfT2HaIhGaG1qRinF/FKWi9evruvS9324fNnk448tUil49dUanZ1bs5Qu3tDAYG8/NdelXK1SqVVJJ9ObzuInxH4i3fhCiMfKNDfj+R4R/X6ynEAFxBsaHjlZb2ZujobY/d3sTMOgXA3v9lubmsjnNT7+2KJY1DhyxOeVVxTZ7MZz229EZ1sb7a2tVKpVLMus71wnxMElwV4I8VhtLRk629qZmp0lErHwfR/LsmhKNdLV1r6u7L0NcB7MomeZJvligVK+levXTSwLzp93aWl5utz2G6HrOvFn3EJXiP1Cgr0Q4rE0TePzFz7DB1cuMzk3QywSJZ1MMnDo8EPBVNf1R+4MVyp7jN9po1w0yWQCTp1yn3YjPCHEU5BgL4R4Il3Xefn0Gc56J3Bdj1g0+okJavp6DnFt5Ba+74d386UYxeUhDnemOX7co7fXf6YN+ZRSLOdyzGcXiUWjdLd3YD7PbX2F2AfkX4gQYsMs03ri+HdzupGaV+P2+DjV4iHyyx00NVY5fSZHZ8ez385fHx0hu7JMNBJhObfC5OwMZ48Pr1uPL4RY75lGy2zbPmvb9r/ZqsoIIfa2u1OT/Ov/+//k/Ys3mBw7xNREgtbWGi++lGd+efSZr7+Sz5FdWSYWjaJpGqZpErEsRsbHtqD2QuxfT7yzt207DvwecA64Cfw3QAb4H4CfB771HOsnhNgAY3QUc2QEzXUJ0mncF19EbSKjXr5Y4PbEeDh73TSJx2IUy2U836dUKRONRDB0g0Q8zlBvHyv5ApNz07iuRzQSwTItRsfHuH57hNxyA8vz/fjKY3BoinijC9oApYr3yPculkpcunGdYrnE3OI8uXwB3dBJJ1IMDw0xcKh3NUHOfDa7bvKfQjGfzbK4skSlWqUpleJIbz/msy7WF2Kf2Ug3/r8AXgT+P+AN4AxwgjDI/5bjOAuPOXcd27YN4F1g0nGcX7JtewD4DtACvA/8A8dxarZtR4FvAy8Di8CvOo5zZ8OtEmK/833MK1cw5ubQZmfRKhWCQ4dQpolWKBD567+m+gu/ABsYyy5XK1y8fo1YNIppGGSXl/loZpojh3upVCssZLNousZQ3wCVSoW/ffdtotEG4rEYpmFwZ3KCpZUVgkCxONfBwmwC3cxx9HgezawAGfLFPM3ph5PaeJ7HOx9/SLVWZWpuhrvTU/ieT0tjE3mtyAdXL+P7AUf7B4BwDf3iUharPglwdn6eXKGAZZiYhsFKocDF61d58eTpLf24hdjrNtKN/yXgi47j/A7wZeB14D9wHOe/2kygr/snwNU1z/8Z8AeO4xwFloCv1Y9/DVhyHGcI+IN6OSFEXeRHP8IYGwuD/t27mJOTGJOT4Yu6Dkph3r69oWtNTE+vm0G/uLJEoh5Ul3MrWJaFhsZSbgVN05hfWqZarQLg+z7FYpFq1eT6lQy5pWY6uspkum6i6eVwk5vAxw/UQ8v0ACbnZjAMA9fzWF5ZCZf1mSb5YhFD16lWq8wuzK9m4OtozUA9HW4QBKwU8gQqoLm++Y+h65QqFXKF/LN8vELsOxsJ9knHceYAHMeZAAqO4/xos29k2/Yh4BeBP64/14AvAH9aL/It4Jfrj7/C/eGBPwVer5cX4sDTFhfRl5ag3p2teR4qGsWYnr5fKBJBy+U2dL2a666bWR/UA6sf+Hj1x4ZhUKuFAV7T1Go6Wi/wWV5KMnmnj1pNp39wgUz7HI2pJLqmoZTCMMKd77oekSe/UgmHDTzPxfN9tHr++ntb0yqlcH1vXT3OnjhJNBqlWqsRBAHtLa20NjWvXtPQNSrVyobaLsRBsZFufNO27c+zZheJB587jvP9DVznfwT+SyBVf94KLDuOc28gbwLoqT/uAcbr1/Zs216pl1/Xk2Db9teBr9fLkclkNlCNvcM0zX3Xpo2Qdn+CIAClYHERraUFotFwy9nWVjTXhWqVWDwedt0XCqizZ+GB6/lBgAbrgvtArcKdiXEiVvjloam5mVot3I2uVKmCCoN7d2cX6WSKhlgDnR0dxGNJpm/FWZyP09hU4eSZKrGGJMs5n1yhQKalhZ6OTr78udfXZdRb1yQNro3coi3Txmx2EdfzMEwD0zCJx+NE/Ahd7R10d3Wt29jmcM8hlFL85L13MB4Yny9XKwwODD6U2Gc3kt/1g2Un272RYD8H/E9rni8+8FwBRx53Adu2fwmYcxznPdu2P1c//Kg7dbWB11Y5jvMN4Bv3Xl9Y2Oyowu6WyWTYb23aCGn3A3wf6/33MaamQCmUYVCcnWGacB/4iOvSMzVNMpHALZXQSiX8jg7KQH7kVrhvvK5zY3SEYqWMrum0NDYx1NePruvEzAi1SpVcLhducmNFuTs+jlEfA5+Zn6UhFmN+YYHZxTlM3eC9j8coLg3R23WEdNMMrnYD122lUKpR8z1832d+fp5SscT/Xipz4dxLoBSWZZGM3584qClFKplkfHKCdCLJ9NwsxXyR7rYOlpaztDY1k2lqZnFx8ZGfWWtjEzdujxKxLDRNo1qr0t3RST6XYy905Mvv+sGy1e3uXrPN9JM8Mdg7jtP/LJWp+xng79q2/WUgBqQJ7/SbbNs263f3h4CpevkJ4DAwYdu2CTQC2S2ohxB7jvXuu+gLC6h6trpyucTyjWuYh/vQGhpQhslEawuZo8cIDJ07rU3crhbJ/+X/S1tzC9FohOm5Wfp7Dq/u876cW+HGnducODKIruucPXGShaVFlnM5VBBgGAYz8/NMz89iaDqlcpm7U1NErAhepYfCSguJ+Apm4v/h7MkuBnrOs7Syguu53Lg9gqZpKMMgXyzy3uWLXB8d4aVTZ9ANnXgsxuljJ7DMMECfGz5FQyRKdilLf/ch8uUiuXyBjtYMA4d7iT0m1V6muYVkPMHU3Ay+H9CZOUIqmdyWn4sQe8lGlt79iEfcVa/lOM7fecLrv0e4fI/6nf1/4TjOf1hfo//3CWfkfxX48/op360/f7P++vcdx9nanTKE2As8D2NmZjXQAyyurOAe6qXS3obhugSGSeHUSe7GooSdYorlerf8zMI8HZkM1VqN2cWF1f3lTdNkaWU5HAbQdXRdp721jaZ0E3/xN9/DDwKqtSrRaBRdCyfnBZ7B0tIR3Gojza1FmlvHsSIRqtUarufS2dbGzMI8+WKBhoY4yUg0vI5bo1arUSyXaG/N4PsBN+7c5tTQMSBMx5up38E/jVg0ypHDfc/6SQuxr22kG/+P1zzWCJfi/adb9P6/A3zHtu3/DvgA+Gb9+DeB/8W27VuEd/S/tkXvJ8Te4vvhWP0aQeCDYVBraaG0phtvdmqcvu7DTM3Nri5NMw2DhWwW0zAolkoArOTzrORzuL7L0eUBMi0tq9eYW1wABdVahZrnoes61VqVcjFBPjsAmDS23iXTUcXzfVw3DNZLKysk4wkqlQqBUgRB+N28Uq1ioOEFPsVyGQjnC+TyeZRS68bhhRDPz0a68dclzbFt+58/eGwzHMf5AfCD+uNR4JVHlKkAv/K07yHEvhGNopLJcGJeXTqZZGF5hXJrK7l8nnypgKEbWPUJdnp96R2ahq7rGLpOpVrFNGF6bo6VfA7TCv/p37gzStWt0dPRuXpuOpnC9xW6BrliidxSJ8sLzRiRCs1tN0nEFRBHKYhGY0StyGrQbkynUUoRi4RfNjQNAkAzdJrT6dU2aJomgV6IbfScNpcUQmyV2vnzaNUqWqUCtRpNukn19Cmuz0wyPjtNvlCkUCpSqVapVKu0NDZRqdXIFQusFHIkEwlamprQDZ3FlSU0XcP3fQ51dBKLRpmYmcbzvHAZW2uGhliMtpZWqhWd7NxxSrkOEuksbZ3XMcwypmnguR66rnHh7It4gUcyEcfzfRqiUY4c7kdpipoXLunzPJ+uTDuN9WDv+T4t9XXxQojtIRvhCLHLqeZmKm+8gT4xge66eIcOkS4Xab9xHc/ziEVjJOINeK7HxNwMDdEY1VqN5dwKuq4xbZq8cOIU7S2tfHD1MvFYlHQyjaZpeL7H7YkxKtUKEStCYyrFUP8AMzMa1VyaqFWgrW2cdHMZ123C0E0iERNDN3jl7As0pZt44cQw0ViMcrlMa1MzL58+y3uXLzE+NY7SdF4YThOPxsNEPJpGUyrNUN/ATn+sQhwoG5mg94UHz3nKdfZCiKdlmgT9/dwbvc9OT9KUut8tXnNdNKAr087SyhJdmTYGD/fREItiGiaBCmhpaqY5nSYWjQFhXvk7kxMEgSLREA/H0gtlbt4ogfsCZ4ZvcfxEmUikB00D07TI5XO8ePIM6VQK13WJRiL31+yvmWD3mRdfpnrqTNilH42ilKoPJRhP3DVPCLH1NnJn/80Hnm96nb0QYmvFolHyhQKe7zE+M43rugB4rs9gby/xB7Z7rVSrmIZBc7qJXLHA/8/efcXalaUHfv/vvE/ONwdexiKrukK3qrsVWqNpybAsDDAv42NggMHYEKCX8YMxfpiBX+xH+8WDAWyMIUAPkmFAvpA9kDwwPKNouVvqVN3VRRbzJW+OJ6ed1/LDPnWLLLKiiiyG9SsQxXPu3meHQ/LbK32fZZr4QYDneczU6ui6zmRscPdOlcFQ8Bt/T2c+8LBMHU1Ll74JIWhUa1TLaY77Tys282BSG03TyLju57rGOEnY2ttlMBygGwaLs7PUK7XP9RmKoqQ+ywQ91d+mKM+Yxdk59o+O2D5I8+F7vs/Y9yhkc9zb2+G1C5ce2l5HZzQZI6RgPBmj6Tq2adGo1pmtNzg6dNneymMYkguXupw9O0OSXOTmvbv4gY8ECtkcl89feCrXJ6Xk6q3rxHGCYRgkScLtzU3CKD5dPqgoymenxuwV5TlkmRbnlle5v7vN5u4uIKhX6mTdDFv7uyzMzBJGISIR5LJZkkTwZ3/zPbKuS61SRQhBKV9ACIPbt4r0Og6lcsjZ8wMMQ6ZZ9zSNb7z2On4QoOvaaTrdp6HT7+OHIc4Dx3Rtm72jAxXsFeULUMFeUZ5h48mEzf1dwjDEtm3OLCyRm3bRxyKtFJfNZNJ19F7aci8Xi7x/+xa1ahVd09g/OcIPQ+ZqDYbjMd3BgJWFBXb2PYbt85y0eqys9JiZHSMkvHLmwkPL4j4pg92TMhyPMI1H/3mK4+Q0EZCiKJ+dCvaK8ozyAp+f37p+uo7d99PXb11+jYzrctxuI7W0Ot1gWuq10+/j2DaX1s5xbmUVPwy4t71NfzAknpaPtTSdq1dDZLTE8qLBP/5HCwRxG12v0ahUHyks81WoFEvsHx1ifORBw7EsFegV5QtQf2sU5SmI45jt/X1u3dvguNNCfCQr3uNs7+09lLBG0zQcy2ZrfxdIJ93NVOv0BoP082Razz2OYvzAxw8D7m9vMZyMAAjCgCgy2NyYZX8/T7E84ju/LKlVTRZmZpmrN56JQA9QKhQoFQqEUVpKV0qJFwScWVr+is9MUZ5PqmWvKF9EPK3MbH76XyEv8Pn5zRvopDnpO/0eR8cnvHrx0ie2UoMweCTLnKZphNNa8qZp4LgOtVKZRAp0NKSW1nIfjEf89P2r+EHAeDImihO6XRuvP4uuB6yeafHGG2eolJ/dojGXz13guNOm3e1gGiaXz50nm8l++o6KojxCBXtF+TyCAPtHP0LvdEBKRLVK+K1vwSeMa9/b2cYyDAbjIb3BAA2NcS7PTKfNbL3xsfvlslk6/T7GAw8EiRCUCgUAFmbm2NrbZX52jt5wQBSF6JrG6tw8P7t5fToxL2Y49mkf19HFIvVaTH1uk2q9wCtr5768+/IE6LrOXL3B3CfcI0VRPhsV7BXlc7C//320IDitQqdNJtjf+x7hr/86AJ1ej639XcIoxLZsVheW8DyPw9ZJur59OunsuN3i5n3jE4P9ysIi7W6XWEpMwyBOEpCSlYVFIF1+d+XcRX52433KhSIZ12G+McPe0SFCCMIwJI4c2ocXCHyHSvWQr71eYrZ+lnKhwEHrmNWFpSd8xxRFeRaoYK8on5E2HmP0+4hc7sM3dR1jMEAbjRhocPPeHVzHxTItpJTcvHcHIQSD4eC0UE26W1r57ZNmllumxddf/Rq7RweMJx7VbIal2XnMB4YOLp+/gGEatLpdLMMgidPJetVSmdBr0Gs3MI0hs2v7uJkBr5x983RcfjKtQqcoyotPBXtF+aziGPmYidBa1ZUAACAASURBVHVSSogidtonOPZHZo/bDt1+jzCOT4O9lBIhEkqFEkEYfmJmOdM0ObP48ZPSNE3j0to5lucmdPt9spksUaJz82aOUb9ALudRmttCw8eyrNNAH8Uxcw+k21UU5cWmgr2ifEayUEA+LjBbFrJUIj46eOyEunKxRBjH3Lq3QbvXIYpjKqUy9XIF2/py8sRnM1mymSz9vsbx7jlsDGbntyiUT4gTk/FEZ2U+fWiI4hjLMtVYuKK8RFSwV5TPSteJfuEXsH/wAySApqFJSfjNb4KuUywUOG61Hupmj+OYmXqD21v3GHkTdMMga1mMJ2N+dPXnXDp7ntXFv/u4uZSwtWVw545JvVLjO9/ZRlBi4jsYuj4tW1vHNA1manXmGzNqvbqivERUsFdeXlGEfnwMjoOo1eAjrfLHEbOz+L/1W+i76Vp3sbQE09b58twCnX4vzXZn2ekkPdtmtlYnDCNMXUfqBhoaxVwBNLh66waLs3MPPSB8XkEA165ZtNs6s7OCK1diJK+wtbfL2JtgWxYr84vkH5xroCjKS0UFe+WlpG9uYr/7bvpCSqTrEv7qryI/ISD6QcDe0SFxHFOvVAmikMHOFvlsjrnGDKZh8OYrr3LSaTMYDSnm52hUa0x8nziJcByHfDY3PaRk7HlsH+xz7c5N1pZWyBc++5p3IQRbe7vcvT/i1k2X4SjgypWE2BC0enPMN2a4cEbVsFIUJaWCvfLyCUPsd99FPlgGVgisH/6Q8Lvffewu3UGfG3dvY1s2Ukp+cu3nZByX1cUl+sMh+0eHvPXqa1imxWy98dCSulwmQ8bNMhxNTt8bjkeMvQmVcpmrt27x7vXrXD5/Adu2efX8xU8dy7+zucW192MOD8q0+lsUq5tsHXsk2hydQQ9d0z5xWd/zZjQe0+p1yTgOjWpNDUEoyuekgr3ycpAS/eQEogjGY+RH08LqOvpgkP78MYH23vYWrpNOzru/u8txu02cRIy9MfVKjXqlyubu7mNb05qm8UtvfoM/+fP/QBjHSJEw8ibUy1XCMKRWSuvDt/tdZio1bm7c4fVXruAFPjv7ewRhSC6bZWV+EdM0GY0E3/++RhgUsNx9Csb7mJZGkiSMPQ8vCLhlbbwwwf7u1n2OO21s00pr3O/v8vqlK19JgR5FeV6pYK+88LThEPuv/xotDNNJdZ0OaBpicfEjG2rwmBZjnCQEYYjrOERRxO7hHqZhoGsaYZzQGw4wDB3rIw8JvcGAzb0deoM+w/GYN6+8ynGnxeFxi/OraxRz+dPc7wBBGKHpGiNvwmA45P27t9G0tNzrJPD48dWfM1f+GseHs3geXHqlz0HnPl4PtOl/cRKTd3IcnZw8kXv5tA1GQ447bdzpkkZb15FScnd7k9cuXPqKz05Rnh8q2CsvPOuHPwTTRE6DsXRdnL/5G4LZ2Q9z28cxSaMBjykEY+j66fr0Vq+LbVlEcYyQEts0MQ2T3mDAwuz86T4Tz+Od99/j5sZd9o4OmAQ+pm4w15jltYuXGI2HHLVatDod8rksxXyBmXpagrY/GHBV3kIkCTuH+wghOGp1ODqso8VtXjkvqc7foVKdpzUwEVKiAxJJxnGRUj4zBW3+ro7brYdq2kPaUzIaj7+iM1KU55MK9sqLLQjQh8OHJ97pOuGlS2jdLmSzoGkk9TrR22+fbjKeTNg+2CeOozQQV2sctVtpMpx8kcPWMaZp4Uy7koMoYnX+w56Cu9v3ef/2LQ5axwxHQ0AnlhFHrRPa3Tb1cgVN15EI+qMRQRThhyG1Ugnbttnc3aE37NOo1Ng9GrK/vUYUOVRqu5y/KPDDLCedNtVKmeNOCy8IqBSLSCRJknBuZeVp3eEnyrYckqT3yGoFw1Bj9oryeahgr7zYdP3xS+pyOcJvfhNZqaTj9/aHrcf+cMC1O7dwbQdN0zhqnWBaJktzc0RRRLffp1Qo0B+Nub+zhZASkQj+lz/8XxGJwHEs/DDkpNNhOB6hAbqetrS7gz75jMvuUUCxUCCKY6IoQoiYjOOAprE8v8D+8TGHJ8dsbencvdtA4FGp3yFbELR6LqsLizi2Qy6TIY5iNE0nEQmWZaGhEUUxP3j3p9imyfLCIo1q7Snd8C/XwswMe0cHD/1DFUYRC7OzX9k5KcrzSAV75cVmWSS1Gvpw+GGXvZRI00TMzT12jP7+3i4Z58NMeaZpMhyNkUKyNL/ASa/DSauFFII4ThiOx+i6jhekJWnHgUar20FOi9FomoauJUgNdE1HCAiiAHSNcr7IbK3OSbvNbGOGxZlZLNOiWmzwt3sRw0EG3WoxO3eArsdYVp5Wt8Ncrc5Mrc7S3DwrC4sctk7oDvogod3tYD/w8HJn8x6maVIplp747f6ymabJ1y6+wt3tTTzfRzcMFmZmWZlf/PSdFUU5pYK98sKLvv1t7B/8AL3dTgN9Lkf4ne88NtADhGGI+cCY93G7TbvXIZ/Nsbm3k07Wc10sy6bb79EfDpAIgkhgmgaWYZJEMVGStriRkkQmSAmGraetfV0nCgxOwjbD8YhivsDYm5BxMwwGFht3ahSzLqH8IXOzE6JEUsgV0TSdfCZLq9/l7TfeAtKiOgszsyzMzHLr3ga57MM1313HZe/w4LkM9gD5XI43L7+KlPKRdMSKonw2KtgrLz7LSoN7GEKSwLQ87cexLQsxLXgTRRHtXgdDN8jnsrR7PTRNY+/ogEqpTCISvCAgmpafTWORhkBiaDq6oRNEEUJKpBBIYZBIgSk1gjjEma7bj5MY27TZ3nY52M8Rxl3qc7cpS40ocTAjA8d2KBeLZFyXtYXlhx5IPhDF0WMDYpwkX8ad/EqpQK8oX5wK9srLw7aZeBOOdrdxbIeZWv2xAXN1cYmrt26QdTMMRiM0NHRDp1wocXBywkm7RRjFhGHISaeNH/jouoGu6URRjJAJpmHiuC5SCDRdI4xiMCSWZaHrBkIkxFGMbVmUC0VEYjPsXubm3gQ3s0O2tMnZ1QtM/AxhFCGlJJfNsjAzix8ELM7NPfYSy8USe4cHDy0DTISgpCrcKcpL7akE+2az6QJ/DTjTY/7R+vr6f9tsNteAPwSqwE+Bf7K+vh42m00H+APgG0Ab+M/W19c3n8a5Ki+uzb0d9o4OcSybJEnY2tvlaxdfeShn/HA84t7ONlEUc7+9QxRFeIFH0S5xcHxMEAaARiGXm7bIxelMcQ3QdNCFhpQiXXvvOuhoDMYjoihOt9PSrndNN7BNG5JZdrbKrL29yMXvHOMnLbxglu5wwMr8Ilt7O0hg7I3xg4ByqUT5Y7rkF2ZmaXU7THwP13YIoxDLsr6UYjuKojy/ntb6lQD47vr6+hvAm8BvNpvNbwP/A/Cv1tfXLwBd4Len2/820F1fXz8P/KvpdoryhXm+z97RIRnHRdfTBDi2ZXF36/7pNnGScO32LTSgVChgGSZ7h4fsHh6ye7jP5t4Ou/sHuI5DMVdgrjFDLuMyV5/h3PIKS/OLlAslFmfnKRWKNKpVKsUSAoFtWWQzLoVsDtAwdJOMk2HcX+Zgd54zSzUWztyjO77J7tFB2mPg+UgpOb+6RqNaI5/J8+qFi1w5d+Fju7R1XeeNV65wce0chXyBtaUV3rr82mN7MBRFeXk8lZb9+vq6BEbTl9b0lwS+C/zj6fu/D/x3wL8B/uH09wB/BPxPzWZTm36OonxurW4Hy3w4w52maUwCHyEEuq5z0mmhT4PozuE+hycn2LbN2JvQHwzQdI1YxPSHfX75698kl83SHw4Zjoc0qnUcy6aVzzHxPeYaM2Rcl9F4TBRFuEWHg/YRpmVRq1TodELaB+dw7Dpvv12g1mgxmmjouk6lmPYinPhtzp05g67r5LM5Lp4595m64zVNo16uUC9Xnsi9VBTl+fPUxuybzaYBvAOcB/5nYAPora+vx9NNdoEP1tMsAjsA6+vrcbPZ7AM1oPW0zld5sTiOjRDJIxnydE0/bSWHUTq5rdvvc9zuMPE8wigkTmJmanXiJCGXyeIHPtv7O1iWg5SSdqfLxPNwHRfTMCkXihRzBeIkTou2GAaD4QA/CPCDgNCbJRhf4tzKMr/4TZvf+NUr/NH/8+8o5AtERJiGSbVcpt3r0BsOMQ2TRrXG8vz84y5NURTlUz21YL++vp4AbzabzTLwb4HLj9nsg5b74/ooH2nVN5vN3wF+Z/r51Ov1L+lsnw2mab5w1/RZPInrrlartPt9jGlOe0iDe6VS4d7uNn4YEsURN+9vYBgGByeH+H5AItKxdyEktm3TqNZIZMzB0QlrKyvUDA3D1Dk6aZMIQb1a4uzqGb77i7/Cjbu3GYyG7P/4BwRRyIW1S2zfK+D7OaqVkG992+c3v/s29UqVC2fPcXRygmGmDyPVSoW3XnuN2foMl89deGHS3z7Oy/rnHF7ea1fX/RUc+2kfcH19vddsNv8K+DZQbjab5rR1vwTsTzfbBZaB3WazaQIloPOYz/pd4HenL2Wr9WI1/Ov1Oi/aNX0WT+q6zywscnfzPmNvgqYbFLI59g8OyLhpAp2tvV0mE4/BaECSpNnoJoM+uWyOo2l63DAMqJYq1CtVFhuzXL19k+OTNkIkJElMHEf8/P336XbarC2vEngeS3ML3Nlos7OxgIgsFpeOOLMaU6+soglJu93G0g0unFmj1W4D4Ng249GY/GKGbrf7pd+LZ8nL+uccXt5rV9f95VhYWPjM2z6VCXrNZrMxbdHTbDYzwG8AN4C/BP7RdLN/Cvzx9Pd/Mn3N9Od/ocbrlb+rjOPytUuX+fab3+Bbr7+JrmunZVLDKKLV65LLZjAtm7nGDBqQy+UYjIf0BgOQgka1wXA8IQwDAHYPD4jiCAAhJXtHRwxHA7r9PkII9k9O2Lpv0D95hVw2zyuv9Xn9tQylYuF0LT/A+ZUzRHGEYRg4tk0QBszU6g+tFFAURfmintZs/HngL5vN5nvAj4E/XV9f/3fAvwD+ebPZvEs6Jv970+1/D6hN3//nwL98SuepvESSJEHTNI7bbTa2NxkMh5x0OvRHAwxNp1wq4dr2dCZ9hiCOOe6cUMxniZKY/nBIkiQkSYLne4RhWq42zbDn0OqMuXY1y+07Fol2SGXmOoV8+mDgui7FQuH0XHLZLL/8jW+yNDtHuVjka5cuc3517Su5L4qivHg0KV+YBrPc39//9K2eI6qr6/Pzg4CDkyMM3WSu0cD+SI35B+0dHnB36z67R4fYlk2nl060C5OIiedj6Dr94YBECEzdoF6tYug6QRSioWHbNsetEzqDfjrJRNexDYtquUQ1f47N+2WEMMgUNllZFmkmPSFYnJvj/OpZ3rr86kPnp77vl8/Leu3qur8c0278z5RaUmXQU14Y+0dHbO5tY01T0O4e7vPK2fNUy+XHbj8/M8t7t2+iAVJKCrkcQkpqmQrv3rpOIZPFtmyiOCLjuuwfH6FrOrZlkclkMAyD3nhEKV9Akmaq8/2AvZ0CorKGYbRZWTvEdSWlYplGpcYk8PnGq6+xMr+o0r8qivLUqGCvvBDiOGZrfxf3wWp1hsHGziaV0hungVUIwcT3MA0T13F45ex57u/u4PkeGdfl4to5Wr0ujaMaZ5aW2Dk8oN8fcNJu4/keruvihwFREhOGAY6ZlpTVAD8w6R6fIYnyvHpuhLQ3MUyNcrGKrhtkXIdiPke1VP5cgX7ieUgkWTejHhAURflCVLBXnntxHPPnP/g+V2/dwDItGtUaF8+sMRyN2G+dEMcJ5VKJYi7PzsE+SZKABoVcnvOra2zt7xDFEffvbnNn8x5JIkBCd9BDIpn4Pl7gk0jJ2PMwdANNS+veu47LKPDwRxUGnWVMQ6NU3aA6KxEiT9bNYFkWURwjpcS2bPLZRyfdhVHEj6++x/s3r2MZBplMFsPQOW610DSN3nBAGMW8fvESb155Fcd2aPe6tPs9ojAkiCJsy2S23qBerqJ/TEU/RVFeTirYK8+9f/+9v+Lmxl16gwGGYdDpdjlqt5it1dBIW/j9QZ+fvPcul89dOB0nD8KQeztbxLGg1+/z/u2b6Tp83cA0TIajEbquE09n239AiCQtTgOMRh6T4RoymkM3erilO5gu3L4n+c1f/TW6wz6e5xMlMcV8gcvnH011OxgO+ZO//FN6/QHD8YDheEwumyOfzRInCa1uBw0NQ9fZO9znztYmS3NzaGjsHR9x3DoBTaNeqZDP5lhZXOKNV66oFLmKopxSwV55rg3HY25s3MG20lnzUkowNHb29qhXyhQLRQzDoN3voek6/eGQSiktImPoOvd3tpmtNxhNJjiOg2XaaFo6ho8Gk8DH0g1AwzRM4iRGAlEcQ5Kne7KKlDmyuS2ypX1s2yKfyRLGEe+8/x7fevPrGJrOhbWz1CvVx17D3/zsJ4zGY0zTIIxiXMfF8z0GoyGu4xAEAZlpD0EYRRy3WwwnQ2ZrDfzAR9PTgYQgjNB1n8lkzNbeLudWVp/eF6EoyjNNBXvludbt9YiiGMdyKBWKDMdj4jgiSiKkgKWZD0vB6mgMxyNymQxoGpZpEkYhpmEQhCGapiGFIExikiRGSEmSJDiGhaZBIhJ0TUNIiT+awRuuUCnlqDS2CJI2WbeEaZpIIcnncxTzBeYbM5xdXv3EsfbuoE8Uxzh2OrFQ0zQMw2QwHiOFRNN1kumafN3QmfgTDENn5E0IwhBDT1vwfuhTzOcZjicMx6OPPZ6iKC8fFeyV51qxWECbjk/rmkYpnwfSlvfy/DyapiGEoDfo89Pr15BSomtQK1do1BsszMzR6rTRNI3+YEAUJ0gkuqYDEl038KOQJEnSGfeJiT88RxxWqdUiVs9s0x91EbHA830Mw8B1XaqlCrZtYRrmp06qMw1zejwwdAM5zQxtmQa6oROFMdY0ja4G2JaNhoZlGOkwQxKDBFMzSZIE27Yw1Ji9oigPUP8iKM+1Ur7A6sICfuAjpEBIgR/6XFg9g23ZJEKwub/DzXsbaLoGSASSo06b/qCPFAmHrRZBFCIkp4E2/SyJSOJpsNZIoiL941cJvSJzC12+9oaP1GJKhQKappEISRhFGLpOvVzB1E2qpRJ3Nu9zY+MOx50Wj8tr8craOSzLJIpjctNx+iRJWJ1fIp/Jpbn8p+dWyOUoF8tUSmVq5SpZ1yVJEuIkOc2259gOi3OfPY2moigvPtWyV55rmqbxD37tP+J77/yQ7YMDAC6unuVbb36dXDbL9sE+P7n2LlnXnQZtSRwnIAS6prN/fMRMrc7Ne3fJZVz8MERKiOMITdfRNY1qsczRYQG/P08+F7O61uHrr68wmoyZrdaY+D6FXIGj9jGGnrbkD0+O+aVvvM1Pb7yPa9lYlsXG1hatTpcr5y88dA2vXXqFRApubNxh7+iAaqnEfGOGtaVVXNdh52CfO1v3SZKEeqXC5XPnWV1cYfdgnyiKkEjiJMGxLZbmFji/ukZ1Oi9BURQFVLBXXgCFfJ5vvfF1SoUNEpGwODdPuVhE13XOLi1TyhcJg5jBaJS24IUgDEMmgU+lWGbncJ9Or0ci0mrLhq5jOA5xlJBEJicHZ4n8HE7ugNrsIaMgZmM7XfeeZCSlfIF8JodtWbS6HbrDPnTh//qLP8W2bVbm5ikViyzNztMd9BmMhhTzH6bK1TSNNy+/ym985+89NrvW65ceVyASZmsvX9UwRVG+GNWNrzz3dg/3+en197FMk1wmy97REX/94x8wGA0xDINiLo9maPiBTxRFBEFAHMckSUKn32U0GWOYOkLK6WS9AIQkDMt0j68wmZjka/eZmT/AMjVs00LTdMrFEhNvzEmnzUm3nSbsCXyybub0OEkSs390SBhG7B4d4FgWnX7vq75liqK8ZFTLXnlmTbwJh60WhqEz35h9JM+9EIIb9+7y//3oByQiIZvNEoUR/fGQMAh559pVZhsNfvHNt/jR1Z/TGw2Jo4g4SQAYTsYPfZ4GBICUGifHs0T+PI4bkC2/B3qAH+hIIOtmaLXTMf9yoUgQj8m6LmEYIoXEMA1EIgjCENuxEUISRCGJSPDDgHwm+3RuoKIoypQK9sozaftgj52DfRzLnpaOPeTyuQtUih+ORd/ZvM87195jOBkjpeSk22UwHDBbaxDEIZqucdQ64Y//7E9J4phiLs94MiGZpp9NR/A/JAERZxj1zmFqFQqlfTL5bRKREIQCXdMwTTNNXysljm2TJJJ8NoNr2ZimiRf4xElMIhI0PU11q2k6SRKfzqCvlitP+3YqivKSU934yjMnDEN2Dw/IOC66rmMaBq7tsLG1CcDE9xiOx7x/5xZhFDIcjZh4Ht1+Fz8IOO620dBxbQfbTMfR0SSuZZPPZDHND55xH14S54/rjDqvg3DIlq5TrOyjGzqGaeDaNgAiEcQiLY2LBlITFHJ5bMfh7dfeIJfJ4tjpeTuWQ9bNYGoafhBSKuS5uHaOMHo4I5+iKMqTplr2yjOn3euerjt/0MT3+NF776bZ66RkY2eLJEko5vIMvPG0HI1k4k0oF4pY5ofd/pZlEceCREo+aM+fLrMTBuPeKoFXxc30KTV2MPQQTbOQMi1vC0wT6ghMaQASwzBAaszUG0RJjGmaVMpldg8OMAwDDdB0jfNra+iaTjFf4PrdO6BBuVji8tnzKoe9oihPhQr2yjPHcVzkNGPcgw5OjjmzsMTYmzAcDpn4Ezw/YL4xS9nQEUnCcdvDMkyQEi/w0XWDxdk52v0u+YyLrmuMJ2OiOMbUdTw/w6R3DpE4lCqHlKrpHIEwhDCOkEJimRrWdM2+bRg4joNr2TiOTcbN4Ng2WSPDN19/M02Nu7KGF6Td967jIkRCp9t9aAb+eDzm/t4O55ZVSltFUZ48FeyVr0ySJOwdHTHy0gluS7Pzaeu4WJxObBOnLd8gDDENk+39PfZPjhEiQUrJaDLmuHOCrht4QYBpmAgNjjvp7HjTNCkVigyGQ/qjQdorAEgJg/4M3mARw4yoz28wN2szHpuEcYwQkkQkaWIdLyaO05z4hqanwwaGjmlYzFQ13r15Dde0ub+7Ra1coVaqMD8zy2HrhK29XTqDHuVcgXeuXcW0TFbm5wmCiP2TIzKOy1y9oVr4iqI8USrYK1+JOEl49/o1EiGwTJPReMxhq8Vbl6+gaRqvX7rMnc37DKaV5+qVCiNvzLvXr+PaNglgmTaGptEbDCjmC0RhhGlagMSxTEYTjyAMOWm38Hw/nTQHiMRk3D9P6OWpVn1efU1iGDOsLCxy4+6daVrcmNFkwnA0Og36hVweLwjQ9TQFrxd77B4dghTYto3ds2l1uti2SW4jC7qOa1m0O212D/apVypk3Qx3Nu+xPLdArVxh52Cf43brY9fSK4qifBlUsFe+EnuHB4hpoIe0DK0Qgq29XRYXFrFMiyvnLz60z/sbt9E1jZE3xvMDojgkTBKCIMAyTYIoxLFt4jgmCAMs28HzJhi6gZAirVbnlxj110BYlBt7VKpDapUrvHbhEquLi3T6fXYO90iiBDSNXDYDaDi2RSFX5KTTwnAchEjwwxAkREnagyClTGfiewlhFPPGK1cYjcfTkrkGo4mHSASOZXHUbrGysIhlmvhBwM7BPkfdNkfHJ2Qch9XFJVzHefpfjKIoLyTVd6h8JYbj0QOz4lO6rjPx/Y/dZ21hBSEF42mLPYxibNPEtix0XUfXNXRNI4pjvCBIk9rEAiEEiYBxf4lB+wK6HlOZeZ9MvkUiBUEYIKRg4nkUs9lppbl0ut8HPQVJIhn7E0CbdvELkGlte2T6WkMjThKESAjCED8I8Hwf0zRPx+7HvocQEk2TlAvF9MKk5MfXfs5wPEaIhOF4xM+uXyOMwid2/xVFebmolr3ylXBsBz8IHhqrllKeJs6RUnLYOmH3cD9NgduYY3lh4XSN+nAyJoNDq9chiCJEnBCJhInnoRsGUso0EGuSIDAZnJwlinJkci3y5V0MEyAdd5dSMvF9lheW+MnV99JKcnGMoRvESVo+NxYJOmkZWUM3kEKkn48GWlpPPoxDRCDQNA3bsTENA8uyCIKAIArRNY2MmyGXzWI79mmvxv7x8WnBGy/w6Q0GAPzovZ8RhBEjb8K5pRUunDmLM10CqCiK8nmoYK98JVYWFvnp+20s02I0HhMnMbbtcGntHL7v8yd/9h+4t7tDEPkEQYChm6wuLbG2uMxwPCKOIgajEX4UEkURD65c/6D2e5IkBF6NfnsJKQWF6l2cTA8JSJGWki0W8hi6QafX5Z3338MLAmzLIo6j0xa7ZQJIBtPjhjJMi+SgoeugSUngByQITMNAkxIpBHe3NynlC4wmY4SUZBwHXddp93tcOXse0zQRQhAlEbXKAvvHh+zs7WGaJtsH+5y0W1TLZRzL5s79e9zZvMev/MK3aFRrX8E3pijK80wFe+Ur4dg2F1bP8v/++G8JwhDLNKmWKwxGI67fv0u71yURMd1ejyCK0LW0i/ze9hazjRmiOGEwGiHko0v0NDSE1Iknr+ANyuRyQ2bm9okSH8PII6Ugn8tRzOeZrc+QcV2iOOb7P/oR42CCbdnomoHlWFTLFTr9HuViCSmgZ/bxAh8p0l4I13HRNOgO+li6hWtbVMtVcpkMnf4A5JBauYqUAsM00YDl+UUK02V4uVyOV89fYDAa0+50sS2b/nDIYDhE0zTCMCLnZtE1wf7xMbfu3aNWrqjZ+4qifC4q2CtPxWA4ZHN/Fz8IcGybM4tL7B7uc2ZxeVp6NvXT61fZOTyk0+3Q6rXx/AAhBHEc0x8OsC0LLwgYDAePrQ0PEEUuw845dJmnPtuhXD1GNzQcmSPruAgpcWyHYr7I3uEhM/UaURRhOTZmHDKejLEtiyiJCeJw+qARY+omWTdDrVQmQVLM5pmt15AStg72KU8DuEROhyg0ctkcCzOzTPwJmqZRyOZxHIfZWo1KqYSQkmqtwfb+O6cJ/TzPS0vWTtf2CyS6rpPECb1Biz4RMAAAIABJREFUHz8MyLqZJ/6dKYry4lDBXnnixpMJ1+7cxHVcTMMgSRLeu3WTKIpwHYeTTos4TkiEYDAcEEY+UkrGE48wDNKlcNOu+SCK8IL0AeDhzPYpbzTDuL+EYSSUGrfJFiMMwyWOE7zQR4h0n+5gwGg8xjRNoiim3etRKRaxTJPQMNENHRINx7IZiBHhcIxjWwiZLsuzHQdd19B0Ax35YTpdKekN+ySJQMiEMAy4tblBvVwBNIajMaVCgSiOsK10n8PjYxbn5jloHeOJCblclkpQYOR56KTPAFJK0CWWbZ2O9SuKonxW6l8N5YnbOdzHsR9eRpZxHHYO9hBSoEkNqUG338MPQxq1Gv1+mgAnEeKhlj9w+t6DLXuRmIx6Zwj9Mrbbo9bYo1h0cO0cURwTiwSRJNMlfgmuY+OFPra00TUNw9CZ+B4Ls/Pc2riDpqUlb33PI+PYRFG61t62LKIoJowijGnO/nKphG3ZbO7vEUUhCIlIBLVS2n2fTCcO5rM5IiFodTtcOX/xtCs+47pMvAmrC4uMx2Nm6w1+FviMfQ/TstDQiJKISqnEytz8Q2mAFUVRPgsV7JUnLoyiRwK2ruv0RwMOj08YjNNCNlEUYlpp5jxjmvJWwqPd9TLt1v5gtD4K8gw7ZxHCJFfaIZM/Ipst4NoOhmHQHw34oI88ThIs00IiMXQdTdeRUuA6Lq7tkHVdDMMkFjG6pmM7NnoicO208I0GaatdCJAaQsRUC0Uur50nl8vz0/ffwzBNVpeWmZ+ZZfdgn7E3IQgCbNtitlFnNJ6gf+R+CCG4dPYcN+7eIU4SLp+9QLlQotPrEsYRjWqNt19/k0tnzj6R70hRlBebCvbKE1fM5zlutzEN4/S9RAhanS5jb8J4MiGOY4SURFHIUeuEKIxOx6rFR/LkS9KgLSV4wwX80QK6GVJp3CaTDbGtPJqWbqPpBisLi3QHfQzNIJGCMAwxDINKscTZ5VUyrsNoPCZKYsaTCXP1BmNvgm1ZhHFMLAKq5TKLs3Mszs6zsb2JoRtcXEsDb5IkjL0xv/Wrf5+VufnTqnae7+M6LoV8nozjsjg7R5wk3Lh7G/2BewGAppHP5Xn1I4mEFEVRvgxPJdg3m81l4A+AOUAAv7u+vv6vm81mFfjfgTPAJtBcX1/vNptNDfjXwG8BE+A/X19f/+nTOFflyxXHMUki2Ds+wDJMZmsNhBCMvUmavOaD7nhNw9ANoiiiPxxi6BoZy8EL03z3cRI/9LlJbDPsrqGJCnOzHvPLffzAZexJJJKMm02T5SQJ5fkFHMdlOBqmlfEyGbJulkI+z0ytRjFfQNc1Xr90hdv3NxhOJmzubGFaFiIRHLZOSJKE2nRmvjatWvcBwzDoDdO18ctzC1y7e5uM45BxXRzbwvN9lucWpltL5hozaaGfaTd+kiQUcjmymQyT8fipfC+Korxcntb6nRj4r9fX1y8D3wb+WbPZvAL8S+DP19fXLwB/Pn0N8J8AF6a/fgf4N0/pPJUvURRHvPP+Vdq9LitzaQrcnaMDZusNzq2sks04BFGIEAJNSye3aZqGEAlhFGHZFsV8DssyHxoGCLwyveMriChHsXqfysw+cRKScVxsy6JcKJDLZllZWGJlcYlqqcQvvfl1Lp89T7VYIp/N8/qly3zz9TdZnl9geX6Bb7z6elq9LpPBMk1WFpfS5DlIGrUac41ZXMdFSkG5WKZRrX7katPzKxWLXFo7i6al17+2tMLlcxfS2fQioV6p8uu/+CtUy2XiJCFOYkqFwiOpgRVFUb5MT6Vlv76+fgAcTH8/bDabN4BF4B8Cvzbd7PeBvwL+xfT9P1hfX5fAD5rNZrnZbM5PP0d5Tmzv76PrOoaeVorrD4cgJe/duk61XCEREpl8kABHIJKENIN9GvCFSCgXSqCN0TQN34sIxqv4owauO6FU28BxY9Dy2KaJQFAsFJmpVtF1fVocJ5xOrLO5cOYsa8urlIslLp5Ze+w5L88vcNRuk3Fc1paWkVIyHI0o5AsYhs7y/CLtXvehh484jqmWy6ev65Uq9cpHHwYedn718cdXFEV5Ep76mH2z2TwDvAX8EJj9IICvr68fNJvNmelmi8DOA7vtTt9Twf45MvYmGLrO1Vs32D7YpzfoM5qM0UjXn3uBz8SbECfJQ/slIn3tByG+HWKbBqOhxaT3ClHoki8d4xZ2kTIhCEHKIUmSUCzkGQxHFLJZ8rkcs7U6tm1zeHKMH4ZoSEqFIudXPr6GvG3ZvHXlVTZ3d/B8H9AQpMv8hBBpVr4wYOJZaBpomka5WOSsqkuvKMoz7KkG+2azmQf+D+C/Wl9fHzSbzY/bVHvMe48sqm42m79D2s3P+vo69Xr9yzrVZ4Jpms/1Nc20G9y+t8HO0SGe7xFGEVGSIIUgHMY4lk3WzeCFAVH0YcJbXdPRpxPsWv0OjrZGTv8ahVmItJ+RcEh02iMQEyAxjAAp8lTLZSzb5te+9UsYhsHYm/BL33ib7qBPp9fF0A1a/R6Xzp7DMAwSkU7Ysy0L44FJc8uLSwBcvXWDwrjw0Oz5QqFAsZBnbWkVyzSxv6R89c/79/1FvazXDS/vtavr/gqO/bQO1Gw2LdJA/7+tr6//n9O3jz7onm82m/PA8fT9XWD5gd2XgP2Pfub6+vrvAr87fSlbrdaTOfmvSL1e51m/puFoRH80pJDLUyoUHvpZKZ/n2u2bBIFPGEV4vgdSpulskwRhCKQQjzzZSSSJkAhhEE8uEoh5ZmaG/IP/eI6f3XC5fT/BttJa9qER8cGzoWPbVCsVHNNiY3uLMExn0f/7v/oLstks1VLa1T4YDDg4OqReqbJ7eECSxOi6wVy9wZml5YfO5fgknZz3Ub7vM1dtMJ6u+f/o0sIv4nn4vp+El/W64eW9dnXdX46FhYVP32jqac3G14DfA26sr6//jw/86E+Afwr899P///ED7/+XzWbzD4FvAX01Xv9sEUJwfeMOg9EI0zCI4z0K2RyvXrz0YbIYx2VhtsFx+yTdSdMwdZ04EQghieIYkcSIR5bRS6Iwx7BzFtvI05g7ZG4xRsgiy3PzdPo9XMdlPB6DD7ZjY6AzU6tTKZY4ap/QqFRx7DTd7OFJi9FkQjGXxzRNDMPguN2m0+tSKhRPM9IdtI5xXYe5+szpuWQcl8Fo+FAu+kQIdCH42fWreEGAYZjMNxqszC8+2ZuuKIryBT2tlv0vA/8EuNpsNt+dvvffkAb59Waz+dvANvCfTn/2f5Muu7tLuvTuv3hK56l8RgfHR4wn49M0saZhMAl8dg73WV1YQgjB9t4uo9GEdq+PbVnTVLlpaViNtBZ8IiWGpp0GfCnBG80xGSxiGBGVmdvMzefIumU2dreJowhD0zEN47QnQdPTGvP5XG6a7EaeloJNkgQpJaZh0On3mKmlXWj9YZ9a+eFJdI5lc9xqPRTszywt8dPr1zAAQ9en3f5pnfp8Nos7zQy4f3SIoRsszs49wbuuKIryxTyt2fjf4/Hj8AC//pjtJfDPnuhJKV/IBwluOoP+Q2lbkyRhOB7h+R5z9Rn+8m+/z89uXOOk08bzPQajAbquE8UxpmGkSWWSBMswkJoGQiASk2F3jSgo4WS6VOq7lEtpdbrOsM+8M0O5WMILAuI4JopjLNsiiiIatSqaBMsxWX6ghW1ZaUs+XdInH7qOXDb76PV9pIqebdl8/crX2DnYx/M9MpkMIk7ojYaPbHfUOlHBXlGUZ5LKoKd8JnEcc2vzHoNhGuS6gx71chXTNBmMRuwfHaBNl9n9/r9d58bGXXqDAUiBputomkYUf5gYJ+u6eH6AqWsYpomMavR7K+jo1Gd3sd1jspkMCzOzWIaJoevMN2aYrddZnlvg5sZdRt6YhUKRrJthZX4B0zKplasIIZh4E3RdR0NjtlZnc2+Xxdk5pJQEYcDFtXN4vv9QVr84jpltNB65dtuyOPfADP5b9zYwHlNi9oNVBIqiKM8aFeyVz+T63duEUXTaPZ7P5Li3s82FM2scnhxhWTZRFBHLhI3tLXqDfpokB41kmgr3A6ZhEgQhIJkEMePjBt5oDsuckK/excokZB0X2zLRdYNKuUylVELXdXRNJ5vJ8M033mQ4HvP6pctUSiX8MMAyTSzTIo5jrt65lQZ8Tcd1HH717W8BaQ/ExTNnKeTz3N26z3G7Pc3iJygVSw9kuvt4M7UaNza6uM6HxX2EEORz+S/3piuKonxJVLBXPpUfBAwnYzKOe/pePpejWMjTHw7xwpCso1PI5/jJtav0BoPTsXJN1x8K9AB+GACQxA7DznniKIebOyZX3EXTBSLRmQQBjuOg6zqNSo3eaMBh64SNnS0alRqGoVPI5SkXi6ct+KNWC9d1qRZLvHX5VUaTMUEYUsoXMB9TFvb86hrL84sMJ2Oyjks289lqxFdKZWqVCu1uJ324SBIs0/zE9fuKoihfJRXslU+VJAlICKOQIAzJuJl0gly+yMW1s1hW2qJ+59rPiaIorQmvG0RJnOaAf4xgUmXUWwUkhepdnEzv9Ge2aaEbBn4QMhgN+PG1d6kWy3SHffKZHEmSpAltdI3+aEhv0Ofg5ARrum7eMk1ef+Uy+WyOfDb3idfm2PZpb8XncWntHKOZOdr9LrlMuqxPf0zXvqIoyrNABXvlUzm2zf7JEWEYpoVkgFqlQrlUolaukM/l2d7fQ063NU2TjOsivAnJR4K9EDrj/grBpI5pDylU7mOY4SPHtExjWhpXJ45jNF3DtV2SJMEyTM4sLePYNjc37hInMflsjjAMmUwmtHod7u5ssTw3z3yjzsr88mNb9n9X+VyOfO6THyYURVGeBSrYK59qa3+PaqnMUesEXUtbr/vHRxTyOd67dR2RCHrDPv3hgCiK0DUtrXb3kUAfRxmG7XMkiUu2sE+msM/jctFEIsFM0nz2AIahMxqP8PyQRCTMNRq0Om28IKA76ONaNr3hgDCOOem202PHMblcDl3TWVta5hff+gavXrj0xO+VoijKs0gFe+VTdQd9KsUSxVye/miIEIJ7Ozv85L33sG0bQ9fZOdzHsW2klMRJzCTwH/oMbzTDuL+ErseU6jexnNHHHk8KgUQw35jHCzxymSxCpsvishkXx7a5vXmParlMrVRm5/CA3qBPnCSIJE1/G8YxmUwW17U46bZ559pVauUKc42Zjz2uoijKi0oFe+VTfdB1bxgGum7w3o33Oem0SaQk47ogIetmODw5JgjTcf0PiMRk1DtD6Jex3R758ia6EX/coQBYmV/ENE1y2QxnlpYYexOEEASZLIZhMBiNAA0hNRq1Grfub6DrOt5khD7t9jcNnTiOQHMJghDLMrmzdV8Fe0VRXkoq2Cufql6psLG9yY27d9g7PmLiTRAiwXUzhGFIu9sljKNHqtdFQZ5h5yxCmORK22Tyxx9zhIeFcYhuGAghyLoZbMsiihMsy2RxZo6Tbpv/v707j5Ejuw87/q2rq++eezjDY4YcksvlclerlaVdW3Aiyc5GdmTLQeyyFMSWDSVyAilOgCCwHSSQkDiIAiRxlMQRoNiKJCOx/ODEsRzLsS1ZpyVL8kq7q10uuTyWx3Duo+/uuvNHFUdDcobkLoczw57fBxiwu7q66v2qmv3revWOrG1TzCfV9H3lCkurK2johFGc9vdP+8/HMbqR3Hq4+baCEELsFZLsxR3Vmw3+/NvfYrlao9VtE4cRmgbNTiepOr8picYxdBrjtBvjGEaXvuEzmJn2Xe9vcXmFbDaLpmk8d/Yl+ssVclmbR44+RMftEgQBi6srDFT6yOdy9FcqeIFPX7nCar0GxHhBgKEbRDEMV/rwg4Cpg4e2+MgIIcSDQZK9uK2XX7nAn/z5l+m4HroWk7ezeL6P63t466alvS4MMjRWDxN4Jez8EoXKFXT91V1R+2GAGQZMz86i6VAqFHE9j8997SuMj4xSyOUZGxphYWUZNBgeGKTrdtE0Hcs00OKYjtslY9kU8nkKuRwPHT7MgbsYMEcIIXqRJPs9quN2uTQ9TdftkslkmBw/cMtY8dOzszx/9iwAYZRMA+v57to0tbqu33BV73b6aK5OAhql/ovY+ZW7Lo+uaWuD7xhasl3DMAjDkNXqKiePHqfV6aAbOlMTk1imSRRFNNstnnjkUX74B36Q1XqNVrtNsVAgjiJanU7Sen9omFz27gbMEUKIXiTJfg/yfI9nT7+Y3gsPuDo7wwtnz/DI8ROMDg5SyOXJ53LMLs1TyGXR0LAMi9AP05nfQnRidHQiIuJIp1U/QLc1gmm1KA1cxDDduy6PoevJD4d07PwoigiDAN1Kl8cRru9hWxnCIFybklbXdTJWhnw2h67rDPb1M9jXv7bdwf4NdyeEEHuOJPs96OrsDJZp0ul2uTwzjWWaVBt1Pv/nX2Z4aIhDY+P0lyt4vk9/pQ8/9AmjkNV6lTCK0EjmnM/l8tSqAbXVI4R+nlxxjnz5GpoW37EMAJZhYNs2vu9jZ2zCsEUUx1gZi8AP8QOfbCZDuVgiY5kEYXDLkLamacjIdUIIcQeS7PegTreLruvMLy+RsTJ03C7tTgfTNNCI8TyPVqfN3OIiK9VVNDRarTa6YeAHAbqmUSoUWV0tU1scJdJCyoMvk8nWN92nBmtJ2TItisUCY4PDVMoViONk2NtanZbbxrYytDsdwjBA103GR0bx/YAgDOgrl9e22fVcDu3bL8leCCHuQJL9HhDHMfPLSyytLNP1XOaXlui4HWqNBoVcnvnFRerNOmgaURRSrTfIZbPUGnWWa1U6nQ4dt0MYhpiGSRjpTF8ZodPuw7JrlPsvYRi3NtbTNR3LNJPZ6nSNTCZDpVCkr1JhbGgEy7Y5OLqPy9euMTI0hOf5rNSqXJm5hqHrVEoDDFT6qDbqZDM2Tz72ehrtNlfnZpjcf5BjE4cZHhjcgSMqhBAPFkn2e8DZixdYrdfodDvMLswTxhF+EFCr1bgyM0O9UccPQvK5LBevXqWvlMwkt1Kv0u10MCyLjpvcg/e9QtJ3PsyQL0+TK85tOOQtJCPeBWFI1jQoF8toGpRLZVqdDuevXmZseJRarYZlWvhhgOu6+EHAxIED+L7P6OAwi6srPH7i5Fr1fV+5gud7jA4NS6IXQoi7JMm+x7U7bZZrVbJ2hisz18hkkjnYdVyanQ4z87MQA7pOu9NC1w1c16Xd7ay1jte63aTvfHOMdn0c3fCoDJ/ByrTuuH/DSOafb7Sa9JcrdNwuYRCQz+WZnptZmy1usL+f2cUFLNOgUBiiv1QhjEJyts1yvXrDvfqMlWF5ZYV9Q8P35ZgJIUSvkWTf4+rNBqauE4URYRih6zpxHDGzMEe1lkwrGxETBT6aphFGEZ7vEa+bgz4ILZqrh/HdMnZuhULfZXQ93GyXa0zDoL9c5o2PvZ5r83NYloWlm3i+R7PdBjS6rovre2iazoHRfQCMjozQ7bj091VYqVYJburPH8fxfZnFTgghepV8Y/a4Qq5AGEeYuoluJPXtrU6H+eVFWt0ucRwTx8kQswBRFCYN3uKYIIrwumUaq0cg1in2XSJbWLrrfRuGSaVUQdM0dE3H9z00S1vrTheGIbVGA9My8XyfjmEQRxHHD09RryeN/eI4Znn1xv76rudy4sjRLTpCQgjR+6QZc48rFYuUCkWCMGSg0o8f+CyvrgIapmGgp5PbRGGI7/tEUXKfPUanVTtIffl4Wm1/+lUlegDLNBkdGk6nxY0ZKPdTyOUIwhDLtNYG9DE0nWwmQ1+pRBRF6UQ3afc+O8sTJx8jjEJcL6lxOH546pYBgIQQQmxOruz3gEeOHk+noLXIZW2a7RYLK0nijoKQIAoJw5AYMA0T3zVZXjlC4OfJFhYolKfR7nLIW03T1m4BZEyTrtulWq8z2NfPgX1juL5Po9UCYkzDJAh8SsUi5WIR0zR5/clTxFGEoetkLIsTR45SyOeZYoIoiqSbnRBCvAaS7HtMFEUsrixTrdfJZXOMj45iGgYT4weYGD8AJNPRnj73MkEYEqYt5iHpC99uVmisHgJiSgPnsXPVV7V/jSThl4slXvfwIxi6Tta2ecdb/xqZTIaZhXn2j4zQaLd5+ZWL5HO5tRns+kplPN/niVOPocW3DswjiV4IIV4bSfY9JIoinj/7UlI9bmWoNerMLMzx+MOPkLXttfUWlpeI4hhNY63FfRTptGqH8NpD2NkOhf4LaHr3Ve3fNAwMw6RUKHLk4EFGBofYP7qPMIpYWFniyMEJJvcfAJIfHVMHJ7k4fZls2kMgCJPW90MDAywvL2/NQRFCCCHJvpcsLC/SdV0yVgYA0zTpui7PvvQCjz98ai3hzy4uYFsZup5LFEYEXp7GyhHC0CZfnqVYmkM3NSzDxvV9dF0nDAJM08IPfAwtadEfkVTVZ7NZCrl8coWeyzPYP0A+m6NYSO6r65pGkI57v96+4WEymQzX5mYIo5Dh8gCHxvejbdZxXwghxGsiyb4H+IHPzPw8Zy9eIJOxME0Dz/OZW1yg67rEWkwYxQz29TM+MkK1VqPabBBFIV57jMbqPtADKkNnsbMt4jjG0DPYmQy6plMuljAMg3KpRLvdIooiGq02uazN2Mgoj504yez8POViCTRotJrYdoZyoQRA1/M4OjG5YdkHKhUGKpW7jjVOp6/VNf2G2gohhBCbk2T/gGu2Wnz92WdwPZ9mq8nVuWt0XDdJ5K5PDGRtm5VqlbmlRbpdj3anjetGVJcO4XX7yGSrFPsuoRsBpmERRCGg4Xk+uq4xMjSE63l0u11s22b/6Bie79NoNuiv9BGGAccOH6FSLOEFPvNLi+iaztziIrquMTUxSaVUvlMod9RoNTl78QKu7wMx+VyeR44eJ2NZ97xtIYToZZLsH2CNZpM/+MKfEAYhmqZzafYq84tLWJZFGIZrf4ZhMDM/h+t7xIDXLVFbOkIUmRQqVygUF4hJGtZFcUSpUKCYy9NqtzEti1q9QTabodZo0FcuUa3XsTMWDx2eIpfNMTgwwOsffmSt+v3lTIbLM9OYpkHWzlKr13E9Fzvz2q/EoyjixfPnyJgmufSKPgwCzlw4x2MnTm7B0RRCiN4lyf4BM7+0yOkL5zl/+RVeuniO5ZVVoijpNqdrWjKyXDsZlx7SrnCQ3DOPod3YT7sxhmF26R96CTvnY+o2hqFTKZSwszaDff1YpslKtUqj3WKltkKumyNfyJOxbAxdp5gvMjYyShiFHNi3by3R1xoNFldXGewbWCtzHMecv3KZR44ef81xV+t1oihk/UdW13Wa7TZ+4GOZcnUvhBCb2ZZk7zjOx4F3AAtKqVPpsgHgd4BJ4BLgKKVWHcfRgI8APwq0gZ9TSn17O8q5250+d44//NLnaLSaXLk2QxDd2OgtimO8m4aWvd4VLgwyNFYPE3gl7PwShcoVdD0iCjX8OCYINepGm+GsTV+5QqPVxPV9PN/HCwJiujQ7baJiQF+5hGmauL5HLptjaF1iX1heJJvJ3FgGTaPZuvM4+rcTE6OxccO9DXrpCSGEWGe7Oi5/Anj7Tct+Gfi8UuoY8Pn0OcCPAMfSv/cBH92mMu5afuAzu7jAH37pc9iZDPVm45ZEv5kY6LQrVBdOEvp5Sv0XKfVfQk8HyUm64GloGgSeRxzBSq3K7MIcnU6HIPDRNQ1DS0bc63geHbeLaejsGxrmdSdO3tD/PZPJEEa3DsBjGsY9HYP+8q2N+OI4Jp/Nyj17IYS4g21J9kqpLwMrNy1+J/DJ9PEngZ9Yt/xTSqlYKfUXQJ/jOGPbUc7dput2uXDlMt98/lm+8Z1nWK3VWK3X8Tzvrt4fRzrN1QkaK0cxTJe+kdPY+RtPg6ZpxFGEBlhWBk2HpHFegOt7yb40DS8I8HyfMAooFYu8+Q1v4ujE4VuS+PjIKGF44w8Rz/cZGRy6hyORVNmfmDqKHwR0ul063S4xcGLq2D1tVwgh9oKdvGc/qpSaBVBKzTqOM5Iu3w9cXbfedLpsdpvLt2M83+fF8y8TE/Hci6fJZbNUG3VanTb1VpNO586D3QR+lsbKFGGQI1ecI1++hqbdWt99vWI8l01GssvZNuVCgdVslna3TTabI47jZMz8IKCQKzA5doAD+zb+/WWZFqeOn+DClct0ul10w2B8ZJSDY+P3ckiA5Or+jY++jma7jWHoFHIyPr4QQtyN3dhAb6MbsxvelXUc530kVf0opRgaurerx93ime8+RyGfZ3F1maHBARaWl5lbmscyTXy3S7Tx4VjTbQ3Tqh1E00LKgy+TydbXXls/dj0k1eu5bA7DNDBNk3KxzOjwCIP9/bx47izVWh1NT2aqG6j0U8jnOLj/wG2P9dDQEIcnJl/zWPamad52+yObvvJgu1PcvWqvxg17N3aJewf2vSN7Tcw7jjOWXtWPAQvp8mng4Lr1DgAzG21AKfUx4GPp03hp6dXNyrYbBUHAtblZshmbZqvFc6df5OrsDK1u547vjSKDVnUCtzOAZdeoDF5G132idb8NdE2DtOV8GEXouoGdsWl32+SzeTQNZuZnef3Dj5DN2FyaucbyygqGrmGaBuViiWIux7WZGeybGuJtlaGhIXrhXL5aEvfes1djl7i3xvj43deY7mSy/wzwHuDD6b+/v275BxzH+TTwJFC7Xt2/F8TEECf3689feoXZxYW7SvS+V6CxMkUUWuTL0+SLc2h6ktgNTSOMIkwjSexRHKfVJzETY/vJZDKYlknBzlNvNwmaLb74ja/jhwEZy2Kwvw/TMMnaGTRN48rMDNNzs5QKJR4+euyeG98JIYS4v7ar691vA28BhhzHmQY+SJLkleM47wWuAD+Vrv5Zkm5350m63v38dpRxt7BMi3wux4vnzrK0usIGDdtvEMfQaY7Rro+n886fwcp8r5ubnl7FW6aJruvkbBvTMPECj4nxA7zlyR8gCENmFuZot9v4YUBohQR+QK1zSd05AAATzklEQVTRRC9otNodRocGWVha4viRo+RzOQC6nsu5y6/w8JGj9+14CCGEuHfbkuyVUu/e5KUf2mDdGHj//S3R7rFcrXLh8iucvniORqOJaZqMj4zyl88/x2qzvuEEMteFoUVz9TC+W8bOrVDou4yuhxi6nvStD0OiKCKbzdJfrjA6NMzC4gJt12VqYpLHTzxKs93G8z1WqzU0XWO1VkPTNPzARzd0DMNkbGSYjGWzf984rueu7d/Qdaq1GnHafU8IIcTutBsb6PWsZqvFpWvTydCxtk0YBly8coVnz7xItVbDC3x0Tecbz32HIAzQtc0bt3ndMo3VIxDrFPsukS0k94EMXU+q7E0T07LQYtBijSiM6HTbWBmbvnQI2+fPfJcD+w4wPjJCvdVkYWUJ4pisnaXd6RATUyoU0DQjmfkuim5o3CeEEOLBsF2D6ux5rXab519+Cddz0TSNWr3GMy98l/NXXsH1PMIoxjItOt0uQRigAXF8ax1+HGu0agepLx9Pq+1PryV6SMaQNzQdyzDI21n6KxUylokf+NQaTSzTxDQM/CCk0WpzbX4WPa0JMAyDIAyJicnaNr7vr42Gt390H57vUUpnsoOkgV+lVJKreiGE2OXkyn6bXJmdwbYya4lxpVbjyrVpZhbmQdOIoghN0/CCZLjbja6fw8CmsTJF4OfJFhYoVK7e0HfeMkzQNYiTQWiK+Ty6btBxXYgj4rRZfpw20AvCED/w6bhd+ktlXM8jCAJcNxm0p79SQdN1Bvv6CaOQYxOHgWSKWWIoFgocnzxyPw+bEEKILSDJfhuEYcjc4jxdz6VSLBOGIc+efoGVeo0oivCCYG3ims102wO0qhNATHngPHa+liT3OCbWNfpLZXLZLAC2aaGbFl23g4aGZZnYpkWc9rGPoph8LkcQhui6jm1lyGZtDh84yGhniP0jo5iWSYzGscnDlAsF7IyNnckQRRGtTiftn5/dhqMnhBDiXkmyv48arSbnLr3CucuX6HbaWFaGlWqNar1Kq9PB8z2CMLxtoo8iPe07P4iZaVDqfwXD9NA1HcMwyGazDFQqPPm6J3jqdU+wXFvlj770Z8wvL5MxLaI4xDYz9JXLSaM7XaNYKKLrGmYQMNw3gGVZ9Ff6mF1cYHhgkIH+fqK0G8D48MgNA+Pouk6pULjvx04IIcTWkWS/xYIg4NrCHPOLiyysrCRX1xoYpsnl2WtYhsmFq5dptlsYunHb+92Bl6excoQwtMmXZsiVZq6Ph4Ou6+h6Mq79w1PHeNtTbyZr21RKJQ6OjxNHEdVmA9uyeWhyikqpQn+5zJlXLtDpdBkY6OOhw1MU8wXqzSaDff0M9A2gaRBFIcVCkalDh17TCHhCCCF2F0n2W8j1XL5z+kV0TePKzDX8MGB+aZF9QyOs1Kq4XZfp6iytVpuYmDBtDLeRTnOUVm0/uh5QGTqLZTeTavQwxLIyGIaBnbE4PjnFEycfJWvbADRaLUYGhpgYO3Dj9twuj588xQ+84Y33/TgIIYTYXSTZb6GL01exTDPppx4GmIZJFMfMLi4QE9PqtgmCANu2cT13w+r7KDRpVifxun1kslWKfZfQjQAdjTCMsCyLfUPDxBpUCiW+79Rj+OtmmctkLDaaXkBDwzTkdAshxF4k3/5bqNPprFXLW4ZJFEeYusHs4jzEUGvU8bwAyzIxDR0v+F6y19HoukVaK0cII4NC5Qq54gKWaWJn8kRhMjhOPptD13VM0+TA+BhoGoV0RDuAfDZHIZcjCIK1KvggDKmUSjLvuxBC7FGS7LdIGIbUm03anQ6VUgkv8PnqM9+k3mzesm7gfu9KXAOiGFrN/bjN/di2z/DYVUyrS9YeopDPoaNTKBQgjgnSMe7Hhkc5NDoOccyB0Runmz117CFevnSRWqOBpmkMVCpMHZq8z0dACCHEbiXJfgt0ul2eO3OaOI5YrVW5cOUS337pBXzfv/Ob4yxe4zh61E+5f4ljxzpM7D/O5IGDfOO573BtboYojrFbTUYHhxkfGQY08rkslVKZY4ePYJo3nkbTNDl59Pj9CVYIIcQDR5L9Fjh3+RUs0yRTKhPFMV9/9i+JbjOm/XVup4+oc4JMxmZg3yzZ/Cq23c/xycN85ZlvUqvXyKXTzrY7Xa7Nz1IuFXnysdcD0PFcLFNOoRBCiNuTflVboL3uXr3reWiajnabce3jSKdVnaRde4hcPmZ438vkC1XyWZtCLs9qrc7y6ioxGppGOvJdTBSGLK6s0mi10DQN13VptdvbFKUQQogHlVwWbgHD+F5iD8OQnJ2l3WkTROEt6wZ+lsbKFFpUoH9wlTe/qUKj248WQ6zBofED+IFHFMdcnyX+eue8WIM4CgnDZLuaphGEt+5DCCGEWE+S/T06f/kVnn/pJZZqq+Rsm33DoxQLBVZqq7es220N06weRNdDjh5b5tTDQ5iGSaPTRNd1yqUyB/ftw/U8irk8Xc8ljkHTNTRI1ylSKhYBMAyDcvpYCCGE2Iwk+3twbX6Or377LymXSsRAq9PipfMvY9s2feUKzXabVrtFGBm0Vicx4n089fgoP/kTk5yYOgjApWvTnL9yCc9zGeofxAt8bNvmnT/0NH/0lS+ysLwEaKBBX6mP44ePEIbJ1f1Dh4/ICHdCCCHuSJL9BpZWV3jh7BlqzQbDA4Pks1mWqlXKpQKPHXuYbDoBzAvnzpDNJCPXVUolKqUScRzjBT5PnHoMt+viewXmru0jN1HhsUctnn7rJOtHyD06McnRiUk832elViWbsdemjf3BJ7+fP/3ql1hcWebY5GGOHpyk0WkTRxEDlT4Mw9ig9EIIIcSNJNmvE4Qhn//aV/nSN79GrdnA8zy6nodpGBQKBfJ2li987Wv8jbf9MKeOPUTXddFvGtteA3w/wNRNqs0BFub6MK2Aiak5Rsf62Wwo/Ew6Mt56+XyeNz9x4/C2djosrhBCCHG3pA54nTMXzvHc2Rdpd7vpgLPJvfIwigj8gDiKcX2Xz3/tK7iey/jIPrwN+tJXioNcODfI/Gwf5UqbQ1NX2DdiUSqWtjkiIYQQQpL9DS5fu0atXsfQdcIoIopC4rQpvB8ERMT4QUjH9bg6N8tjDz1MsZCn3e0QxRGdboeCPUnU+n68To7hsTlGxq9RytsUCgUmxg/cvgBCCCHEfSDV+OuYloW+1j/+e/XtGqCvm51OJ8a2MtiZDD/+tqd5ZfoqV+fm8FsTRN4hioci9u2fpdGx0A2D8aFhBvr6pTGdEEKIHSHJfp0TRw5z9uJ5Ll65jGkYhKaJFgREcYyVsdE1Hcs0KZfK7B/dByTd4caGJlmcOYbna0xMhBw/HqDrw8Dw7XcohBBCbANJ9utMjB/krW/6frquy/TcDACZTAbdMMiYFnbGYmxkmL/19DvWrtJnZnTOnLHQNHj8cZ+RkVunrRVCCCF2kiT7mzx64mEeOf4QzVYL3/fJZbNoho7nepimSSGfByAI4MwZk5kZg/7+mFOnPNbNNCuEEELsGpLsN5CMVHdjy/mcnV17XK9rPP+8RaejMTUVcORIuGmXOiGEEGKnSbJ/lS5fNjh3ziSTgTe8wWNgIL7zm4QQQogdJMn+LnkevPCCxdKSzshIxMmTPpnMTpdKCCGEuDNJ9ndhZUXnhRcsPA9OnAg4dEhmmhNCCPHg2LXJ3nGctwMfAQzgN5RSH97uMsQxXLhgcvGiQaEQ8+STPqWSVNsLIYR4sOzKUV4cxzGAXwd+BDgJvNtxnJPbWYZOR+Nb38pw8aLB/v0hTz3lSaIXQgjxQNqtV/ZvAs4rpS4COI7zaeCdwOnt2PnSks53v2sRx/Dooz5jY9J3XgghxINrtyb7/cDVdc+ngSe3a+eWFVMsxpw65ZPLydW8EEKIB9tuTfYb9Vq/Jes6jvM+4H0ASimGhoa2ZOdDQzA1tSWbuiemaW5ZTA8SiXtv2atxw96NXeLegX3vyF7vbBo4uO75AWDm5pWUUh8DPpY+jZeWlrahaNtnaGiIXovpbkjce8tejRv2buwS99YYHx+/63V3a7L/FnDMcZzDwDXgXcDf3tkiCSGEEA+mXdkaXykVAB8A/hh4KVmkXtzZUgkhhBAPpt16ZY9S6rPAZ3e6HEIIIcSDblde2QshhBBi60iyF0IIIXqcJHshhBCix0myF0IIIXqcJHshhBCix0myF0IIIXqcJHshhBCix0myF0IIIXqcJHshhBCix2lx3DNTuPZMIEIIIcRd2miW2Fv00pW91mt/juM8s9NlkLglbolbYpe4d3Xcd6WXkr0QQgghNiDJXgghhOhxkux3t4/tdAF2iMS9t+zVuGHvxi5xb7NeaqAnhBBCiA3Ilb0QQgjR48ydLsBe4zjOJaABhECglPo+x3EGgN8BJoFLgKOUWnUcRwM+Avwo0AZ+Tin17XQ77wH+ebrZX1VKfXI747gbjuN8HHgHsKCUOpUu27JYHcd5A/AJIAd8FvhHSqkdr6raJO4PAX8PWExX+2dKqc+mr/0K8F6Sz8QvKqX+OF3+dpJjYgC/oZT6cLr8MPBpYAD4NvAzSilve6LbmOM4B4FPAfuACPiYUuoje+R8bxb7h+jtc54FvgzYJLnkd5VSH9ysrI7j2CTH6Q3AMvDTSqlL6bZe1fHYSbeJ+xPAXwVq6ao/p5R6drd81uXKfme8VSn1uFLq+9Lnvwx8Xil1DPh8+hzgR4Bj6d/7gI/CWsL8IPAk8Cbgg47j9G9j+e/WJ4C337RsK2P9aLru9ffdvK+d8gk2Lsuvpef98XVf+ieBdwGPpO/5r47jGI7jGMCvkxyXk8C703UB/m26rWPAKsmX5E4LgH+ilHoYeAp4f1revXC+N4sdevucu8DblFKvAx4H3u44zlNsXtb3AqtKqaPAr6XrvdbjsZM2ixvgn64738+my3bFZ12S/e7wTuD6lfkngZ9Yt/xTSqlYKfUXQJ/jOGPAXwf+VCm1opRaBf6U3fPFt0Yp9WVg5abFWxJr+lpZKfX19Bfvp9Zta0dtEvdm3gl8WinlKqVeAc6T/Md/E3BeKXUxvYL7NPDO9CrhbcDvpu9ffwx3jFJq9vrVilKqAbwE7GdvnO/NYt9Mr5zzWCnVTJ9a6V/M5mVd/1n4XeCH0the1fG4z2Hd0W3i3syu+KxLst9+MfAnjuM84zjO+9Jlo0qpWUi+OICRdPl+4Oq6906nyzZb/iDYqlj3p49vXr6bfcBxnOcdx/n4ul/wrzbuQaCqlApuWr5rOI4zCbwe+AZ77HzfFDv0+DlPr8CfBRZIktUFNi/rWnzp6zWS2B6477mb41ZKXT/f/zo937+W3raAXfJZl2S//d6slHqCpGrn/Y7j/JXbrLvR6EjxbZY/yF5trA/aMfgoMEVS7TcL/Pt0eU/F7ThOEfhfwD9WStVvs2pPxQ0bxt7z51wpFSqlHgcOkFyJP7zBatfL2rNxO45zCvgV4ATwRpL2Cr+Urr4r4pZkv82UUjPpvwvA75H8B5lPq25I/11IV58GDq57+wFg5jbLHwRbFet0+vjm5buSUmo+/YKIgP9Gct7h1ce9RFINaN60fMc5jmORJLv/oZT63+niPXG+N4p9L5zz65RSVeCLJG0WNivrWnzp6xWS210P7Pfcurjfnt7OiZVSLvDfee3n+7581iXZbyPHcQqO45SuPwaeBl4APgO8J13tPcDvp48/A/ys4zha2gCkllaF/jHwtOM4/WnV4NPpsgfBlsSavtZwHOep9L7fz67b1q5zPeGl/ibJeYck7nc5jmOnrZiPAd8EvgUccxznsOM4GZIGTJ9J7+F9AfjJ9P3rj+GOSc/BbwIvKaX+w7qXev58bxb7Hjjnw47j9KWPc8APk7RX2Kys6z8LPwn8WRrbqzoe9z+y29sk7jPrftRqJPfY15/vHf+sS9e77TUK/J7jOJAc+/+plPp/juN8C1CO47wXuAL8VLr+Z0m6a5wn6bLx8wBKqRXHcf4VyX8GgH+plLrbBmHbxnGc3wbeAgw5jjNN0vL0w2xdrP+A73VP+aP0b8dtEvdbHMd5nKQ67hLwCwBKqRcdx1HAaZJW3e9XSoXpdj5A8oVgAB9XSr2Y7uKXgE87jvOrwHdIEs1OezPwM8B303uZAP+MPXC+2Tz2d/f4OR8DPpm2mtcBpZT6v47jnGbjsv4m8FuO45wnuaJ/F7zm47GTNov7zxzHGSaphn8W+Pvp+rvisy4j6AkhhBA9TqrxhRBCiB4nyV4IIYTocZLshRBCiB4nyV4IIYTocZLshRBCiB4nyV4IIYTocdLPXggBgOM4zXVP8ySze4Xp818gGezkqFLq76TrxyQj4u2/PhZ6OjLaDDCslNLSZV8kGVktWLf9Lyilfuz+RSOEWE+SvRACAKVU8fpjx3EuAX9XKfW5dcs+tMHbqiTzPPxB+vxHSaY1Hb5pvQ8opX5jK8srhLh7Uo0vhLgXv0UynOd1P0syJacQYheRK3shxL34P8A/vD5WOPCDwIeAX92xEgkhbiHJXghxL7okVfg/TTIm+GfSZTf7T47j/Lt1z/+zUupfbEP5hBBIshdC3LtPAf+GJNn/0ibr/KLcsxdi58g9eyHEvfoKyUxgo8BXd7gsQogNSLIXQtyTdE7yHwN+PH0shNhlpBpfCHHP7mKe8f/iOM5/XPf8rFLqDfezTEKI75H57IUQQogeJ9X4QgghRI+TZC+EEEL0OEn2QgghRI+TZC+EEEL0OEn2QgghRI+TZC+EEEL0OEn2QgghRI+TZC+EEEL0OEn2QgghRI/7/xbwbZeLSEukAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 576x432 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"tm=np.array(df[:N].sort_values(by=['TIME'])['TIME'].tolist())\n",
"inx=sm_fr[mask].index.tolist()\n",
"carray=np.array([c[0]]*N)\n",
"carray[inx]='red'\n",
"\n",
"fig, ax = plt.subplots(figsize=(8,6))\n",
"plt.scatter(df['TIME'][:N], df['HR'][:N], c=carray, alpha=0.3)\n",
"plt.plot(tm, (lambda x: x*m.params[1] + m.params[0])(tm), c='b', alpha=0.3)\n",
"plt.xlabel('{}'.format(model.exog_names[1]))\n",
"plt.ylabel('{}'.format(model.endog_names))"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"hide_input": false,
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.6"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": false
},
"varInspector": {
"cols": {
"lenName": 16,
"lenType": 16,
"lenVar": 40
},
"kernels_config": {
"python": {
"delete_cmd_postfix": "",
"delete_cmd_prefix": "del ",
"library": "var_list.py",
"varRefreshCmd": "print(var_dic_list())"
},
"r": {
"delete_cmd_postfix": ") ",
"delete_cmd_prefix": "rm(",
"library": "var_list.r",
"varRefreshCmd": "cat(var_dic_list()) "
}
},
"types_to_exclude": [
"module",
"function",
"builtin_function_or_method",
"instance",
"_Feature"
],
"window_display": false
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment