Skip to content

Instantly share code, notes, and snippets.

@mcohen01
Created May 13, 2020 17:34
Show Gist options
  • Select an option

  • Save mcohen01/0606ee4307ab20bda2535a122838d8a1 to your computer and use it in GitHub Desktop.

Select an option

Save mcohen01/0606ee4307ab20bda2535a122838d8a1 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Daily New Cases"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import math\n",
"import os\n",
"\n",
"import numpy as np\n",
"import pandas as pd\n",
"import scipy.stats as stats\n",
"import matplotlib.pyplot as plt\n",
"import matplotlib.dates as mdates\n",
"import us as us_states\n",
"\n",
"plt.style.use('bmh')"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<style>.container { width:90% !important; }</style>"
],
"text/plain": [
"<IPython.core.display.HTML object>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from IPython.core.display import display, HTML\n",
"display(HTML(\"<style>.container { width:90% !important; }</style>\"))"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/Users/michael/.pyenv/versions/3.7.5/lib/python3.7/site-packages/ipykernel_launcher.py:23: SettingWithCopyWarning: \n",
"A value is trying to be set on a copy of a slice from a DataFrame.\n",
"Try using .loc[row_indexer,col_indexer] = value instead\n",
"\n",
"See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n"
]
}
],
"source": [
"def good_states(x):\n",
" non_states = ['Evacuee', 'Islands', 'Recovered', 'Princess', 'Guam', 'Samoa']\n",
" return x is not None and x.split(' ')[-1] not in non_states\n",
"\n",
"\n",
"def lookup_state(x):\n",
" s = us_states.states.lookup(x.split(',')[-1].strip())\n",
" if s is not None:\n",
" return s.name\n",
" return None\n",
"\n",
"\n",
"def parse_daily_reports(path):\n",
" states = None\n",
" for f in sorted(os.listdir(path)):\n",
" if f[-3:] == 'csv':\n",
" df = pd.read_csv(path + f)\n",
" df['Date'] = pd.to_datetime(f[0:-4])\n",
" try:\n",
" us = df[df.Country_Region == 'US']\n",
" except:\n",
" us = df[df['Country/Region'] == 'US']\n",
" us['Province_State'] = us['Province/State'].apply(lookup_state)\n",
"\n",
" idx = us['Province_State'].apply(good_states)\n",
" us = us[idx]\n",
"\n",
" if states is None:\n",
" states = us\n",
" else:\n",
" states = pd.concat([states, us])\n",
" return states\n",
"\n",
"\n",
"all_states_df = parse_daily_reports('./csse_covid_19_data/csse_covid_19_daily_reports/')\n",
"\n",
"def daily_new_cases(state):\n",
" if state:\n",
" df = all_states_df[all_states_df.Province_State == state]\n",
" else:\n",
" df = all_states_df\n",
" df = df[df.Date > '2020-03-01']\n",
" diff = df.groupby(by='Date').Confirmed.sum().diff().apply(lambda x: max(x, 0))\n",
" tail = diff.rolling(10).mean().tail(20)\n",
" slope = stats.linregress(range(len(tail)), tail).slope\n",
" return slope, diff, state\n",
"\n",
"\n",
"def plot_state(diff, ax):\n",
" diff.rolling(10).mean().plot(ax=ax, color='red')\n",
" ax.bar(diff.index, diff.values)\n",
" ax.xaxis_date()\n",
" ax.xaxis.set_major_formatter(mdates.DateFormatter('%m-%d'))\n",
" plt.gcf().autofmt_xdate();"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Daily New Cases in US"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAuAAAAG0CAYAAABzKZJ1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdeXxU1fnH8c+ZLZlkkgBJBIoKLoAiKhVQ6/JDpW6IIuCCWy1urcXWVtzrVq241Fqtu3Xf9xUVRUVcKqggbiiKSwRLYhJIMpNttvP7Y4YYMECAmTuZzPf9es2Lyb0z9z4n93LyzJlzn2ustYiIiIiIiDNcmQ5ARERERCSXKAEXEREREXGQEnAREREREQcpARcRERERcZAScBERERERBykBFxERERFxkBJwEZEuwhhzjzFmepq2faMx5o10bFtERNaPEnARkY1gjHnDGHNjB8t/a4wJrefmTgeOXde208UYc6gx5l1jTJ0xJmSM+cIYc0e79RvSJowxexljrDGmLLURi4hkJyXgIiJdhLW23lpbl4l9G2NGA48DzwO7Ar8EzgJMJuIREenOlICLiDhg5fQSY8zpxpgfjDErjDF3G2MKVn/NyufAKGBKcvTYGmMGJNcNMca8YIwJGmN+NMY8bIzp0247bmPMNcl9rDDGXAe41xHiwcBca+00a+0X1tqvrLXPW2tPTG5zL+BuoLBdPJck1x1rjHm/XTyPG2P6JdcNAGYl91GdfN89yXXGGHO2MeZrY0yzMeYTY8yx7YMyxlxkjKkwxrQaYyqNMfet569eRKTLUQIuIuKcPYGhwK+BI4HxJKaddOR04F0SSW/f5GOJMaYv8CbwKbBzclsB4FljzMo+fSpwMvA74Fckku9j1hFbJbCNMWbHNaz/L/BnoKldPNck1/mAi4EdgbFAGfBwct0SYGLy+XbJ961s89+BE4EpwBDgCuA2Y8xBAMaYicCZwB+Agcltv7eOdoiIdHmeTAcgIpJDGoDfW2tjwOfGmMeB0SQSz1VYa+uNMWGgyVpbuXK5MeZU4CNr7Tntlv0GWA6MIJGg/hm42lr7WHL96cD+64jtBhIfEBYYY5YCc4FXgQestSFrbdgYU58I7ad4krHe1e7Hb5Ixfm6M2dRau9QYszy57kdrbU0ypkLgDGA/a+1byfXfGmN2JpGQvwD0B5YBr1hrI8D3wAfraIeISJenEXAREecsTCbfK/0P2GQ9tzEc+L/kRZKh5EWRS5LrtjLGlJAYZX535RustXESCfUaWWsbrbUHAVsDfwPqSHww+MwY03tt7zXG7GSMeTY5VSTIT0ny5mt52xAgH5ixWltOBbZKvubx5Gu+NcbcaYw53BiTt7ZYRESygUbARUQ2TgNQ0sHyHkD9assiq/1sWf+BEBeJ0eEzO1hXtQHbWzUga78GvgbuMMZcDnxJIim+pKPXJ0eyXyYxWn4c8COJKShvkZiasiYr4zyYxMh2e5FkLEuMMYNJfEvwa+CfwMXGmF2stY3r3TgRkS5CCbiIyMZZBIwxxhhrrW23fKfkuo0R5ucXT84HjgAqktMyfsYYs4xEJZPXkz8bEvPFl63n/r8jMec7sJZ4tiGRcJ9vrf02ub8JHbSD1d67EGgF+ltrX19TANbaFhIfOF4wxlxJYq767sAr69kWEZEuQwm4iMjGuQU4DbjBGPMfoAUYAxwFHLKR2/4O2DlZSSREYp73TSQusHzUGHMVUA1sSSIpn2qtDQLXA+cZY74EPiFxEWNf1pKAJyuaFAAvAhUkRvD/RCL5fq5dPPnGmH2BD0kk59+TSKRPM8bcBGwLXLba5itIjPYfZIx5Hmi21gaNMdcA1yQ/ILyZ3NeuQNxae7sx5rck/k7NTbb/SBKj41918vcnItIlaQ64iMhGsNZ+A/wfiSodr5C4CHIScLi19qWN3Pw1JEaPF5JItDe31v6PxAhwHJgBfEYiKW9NPiAxVeNu4A4SyasLeHAd+5oNbAHcC3xOYlrJAOAQa+2bybb+F7iVRIWTauBsa201cDxwaDLOi0lcXNnGWvtDcvnlJKbJrLy50IUkpracmWzHTBIVU75Nrq8jUSXlLRJVXyYCE1aOtIuIZCuz6jemIiIiIiKSThoBFxERERFxkBJwEREREREHKQEXEREREXGQEnAREREREQcpARcRERERcZAScBERERERBykBFxERERFxkBJwEREREREHKQEXEREREXGQEnAREREREQcpARcRERERcZAScBERERERBykBFxERERFxkBJwEREREREHKQEXEREREXGQEnAREREREQcpARcRERERcZAScBERERERBykBFxERERFxkBJwEREREREHKQEXEREREXGQEnAREREREQcpARcRERERcZAn0wE4adasWTY/Pz/TYYhINxeNRvF4cqp7FZEMUF/TtTU1NdWMHj26vKN1OXXUPB4P22yzTabDEJFurqamhrKyskyHISLdnPqarm3+/PkVa1qnKSgiIinmdrszHYKI5AD1NdkrpxLwWCyW6RBEJAc0NDRkOgQRyQHqa7JXTiXgXq830yGISA4oL+9wyp+ISEqpr8leOZWAR6PRTIcgIjlg+fLlmQ5BRHKA+prslVMJuIiIE6y1mQ5BRHKA+prslVMJuEr1iIgT9LWwiDhBfU32yqkEPBKJZDoEEckBVVVVmQ5BRHKA+prslVMJuMr1iIgTAoFApkMQkRygviZ75VQCLiIiIiKSaTmVgKsOuIg4IRQKZToEEckB6muyV04l4KoDLiJO6N27d6ZDEJEcoL4me+VUAq464CLihOrq6kyHICI5QH1N9sqpBFxExAnGmEyHICI5QH1N9sqpBFx1wEXECb169cp0CCKSA9TXZK+cSsBVB1xEnKCvhUXECeprsldODQmrDriIOKG4uDjTIYh0K/vd8eFa179y0i8diqRrUV+TvXJqBFxExAkqeSoiTlBfk71yKgHXiSoiTmhsbMx0CCKSA9TXZK+cSsBVB1xEnNCnT59MhyAiOUB9TfbKqQRcF2GKiBMqKyszHYKI5AD1NdkrpxJw1csUESfo2zYRcYL6muzlSAJujMk3xrxnjPnIGPOZMeZvyeX3GGO+NcYsSD6GJZcbY8y/jTGLjTEfG2N2aret440xXyUfx7dbPtwY80nyPf82HWTbqoIiIk4oKSnJdAgikgPU12Qvp0bAW4F9rLU7AsOAA4wxuybXnWWtHZZ8LEguOxAYmHycAtwCYIzpBVwM7ALsDFxsjOmZfM8twMnt3nfA6kHoVvQi4oSamppMhyAiOUB9TfZyJAG3CaHkj97kw67lLeOA+5LvmwP0MMb0BfYHZlprl1trVwAzSSTzfYFia+0ca60F7gMOXX2jGgEXESdoVEpEnKC+Jns5diMeY4wbmAdsDdxkrZ1rjDkVuNwYcxHwGnCutbYV6Acsaff2pclla1u+tIPlq0jk5iIi6RUOhzMdQtbTjVdE1k19TfZyLAG31saAYcaYHsDTxpihwHlAJeADbgfOAS5NVwy1tbXsvvvueDweYrEYEyZMYMqUKVRWVlJYWIjb7aahoYHy8nKWL1+OtZby8nKqqqoIBAIAhEIhevfuTXV1NcYYevXqRXV1NcXFxcRiMRobG+nTpw+VlZV4vV5KSkqoqamhpKSEcDhMc3Nz23qfz0dRURG1tbX07NmT5uZmWlpa2tbn5+fj9/tZsWIFpaWlBINBwuFw23q/34/P56O+vp6ysjLq6+uJRCJt69UmtUltykybotEora2t3apNTh+nUWVhvgy66euPU+SxzKvzMLxHlNqwi2DUUFFRkXVt6o7Hyak27VEaJmoNX4XcDC2OsrTZhdcFvfPizKvzUFFRkXVtSsVxCoVCFBQUdKs2dafjtDYmE6PCyRHvJmvtNe2W7QWcaa0da4y5DXjDWvtwct0iYK+VD2vt75LLbwPeSD5mWWu3SS4/qv3rVnrnnXfskCFD0to2EZHW1lby8vIyHUZW0wi4tKfzoWPqa7q2+fPnzxs9evSIjtY5VQWlPDnyjTHGD+wLfJGcu02yYsmhwKfJtzwH/CZZDWVXoN5auwx4GdjPGNMzefHlfsDLyXUNxphdk9v6DfDs6nGoDriIOEG1eUXECeprspdTU1D6Avcm54G7gMestdONMa8bY8oBAywAfp98/YvAGGAx0ARMBrDWLjfGXAa8n3zdpdba5cnnfwDuAfzAS8nHKlyunCp7LiIZ4vP5Mh2CiOQA9TXZy5EE3Fr7MfCz74estfus4fUWmLKGdXcBd3Ww/ANg6NriUAIuIk4oKirKdAgikgPU12SvnMpIVQdcRJxQW1ub6RBEJAeor8leOZWAezyOFX0RkRzWs2fPdb9IRGQjqa/JXjmVgMfj8UyHICI5YF3lp0REUkF9TfZSAi4ikmItLS2ZDkFEcoD6muyVU3MyvF5vpkMQkRzQp0+fTIcgIt3Q6vXQA544oejytp9ztR56NsqpEXDVARcRJ6g2r4g4YXgPFZfIVjmVgKsMoYg4IT8/P9MhiEgOqIuYTIcgGyinMlIl4CLiBL/fn+kQRCQH1IaV12SrnDpyqgMuIk5YsWJFpkMQkRywVWEs0yHIBsqpBFx1wEXECaWlpZkOQURywJdBd6ZDkA2UUwm4yhCKiBOCwWCmQxCRHNDXr7wmW+XUkLAScBFxQjgcTvs+Vi9H1p5KkYnkhiKPzXQIsoFyagRcdcBFxAmqAy4iTphXl1PjqN1KTiXgqgMuIk5QHXARcYLqgGevnErAVYZQRJygMoQi4gSVIcxeOfXdhTEqWC8i6efz+TIdQkasbV46aG66pF+uXRsRjCqvyVY59dEpFlO9TBFJv/r6+kyHICI5YECB8ppslVMJuOqAi4gTysrKMh2CiOSAz4PKa7JVTh05jYCLiBPq6+spLCzMdBiSQ3Jt6oUk9C+I8WNrTo2ldhs5lYBbq3qZIpJ+qri0dkoWRVKjwK28Jlvl1Mcm1QEXESeoDriIOEF1wLNXTiXgGpUSESeoDriIOEF1wLNXTiXgbrc70yGISA7Q/G8RcUKV5n9nLX13ISKSYhvzYV+1tKU9nQ+yNpF4piOQDZVTH51UBUVEnNDQ0JDpEEQkB2zqVwaerXIqAddFmCLihPLy8kyHICI54NMGTWTIVjmVgEejulhBRNJv+fLlmQ5BRHLAwIC+2c9W+ugkIpJiuueA80x9Pe733sOWlhLffHNsaSkYk+mw1kr10GVjeYz6mmyVUwm4bkUvIk7QFBTnuCoqyLv1VvIefBATCrUttwUFxDfbjPjmmxPbfPO25/HNNyc2ZAjk52cwapHU+ERTULJWTh051QEXESdUVVXRv3//TIfRrbnfe4/8m2/GO306Jp64EC26004QieD6/ntc9fW4Fy3CvWgRq1/9E+/bl6ZrriFy4IHOBy6SQsNKosyu8WU6DNkAOZWAqw64iDghEAhkOoRuyR2LMeqztznmzccp/v5zAKzXS+sRR9B66qnEtt++7bWmvh7XkiWJZHzlY8kS3AsX4v72WwLHHEN43DiarrwS27t3ppokslGWteTUpXzdSk4l4CIikn0KWpo45IOXmPTWU/RbkbjLaLxHD1onT6b1pJOwffv+7D22pIRYSQmxoUNXXRGLkXf77finTcP37LN43niD5ksvJXzssV1+zriIdB859dFJdcBFxAmhdnORZeNss/RLHvvnZKY+dzP9VlTyfWk/rjr0T9R/8gktF17YYfK9Vm43raeeSsN//0vk17/GVV9P4emnEzj0UFxff52eRoikSd981QHPVjmVgKsOuIg4obemNKTEPh/P5j+3/Jne9TUs3HQwU4+/lMPPupsndhsHhYUbte34ZpsRevRRQrffTry0FO9bb1G8557kXXcd6HohyRIL6jWRIVvlVAKuOuAi4oTq6upMh5DdrOWE1x7gqgcuJT/SyrMjD+DEP1zPm9vtTtyVwmt5jCFy2GE0zJlD66RJmJYWCi69lKLRo3HPn5+6/YikyfbFymuyVU4l4CIiTjCaS7zhWloo+N3vOPXlu4kbw3UH/Y6/H3YmUU/6vsG0paU03XwzwSefJNa/P55PP6Vov/3wn3kmZsWKtO1XZGNFrfqabOXIdxfGmHzgTSAvuc8nrLUXG2O2AB4BSoF5wHHW2rAxJg+4DxgO1AJHWmu/S27rPOBEIAb8yVr7cnL5AcD1gBu4w1p75epxqA64iDihV69emQ4hK/UKLqfokEPwfPABjT4/Fx59Pm8N2c2x/Uf33puGt9/Gf9VV5N1yC/l33YXvmWdovuACwscdBxmopOWLhNnmh6/YvuIztv9+IT1D9bR6fYQ9Plq9Pgo+7At5eZyxuCGxzOMj5A/wyrC9qS3SedjdfRVSdbds5VRG2grsY60NGWO8wNvGmJeAM4B/WWsfMcbcSiKxviX57wpr7dbGmEnAVcCRxpghwCRgO+AXwKvGmEHJfdwE7AssBd43xjxnrV3YPgjVARcRJ1RXV6sO+HraetnXXHv3BXjqfiS26aacdPhFLO67lfOBFBbSfOmltB51FAXnnov3rbcoPOMM8u67j6arriI2cmRad1/aUMsOFZ/hv/ApPO+9x6z5C/DF1vK366PEP0ettvjUGXfxyB7juX/UkWmLtSvzRcJsUfUd3/QZQMTTfetkDy1WHfBs5UgCbhP3ZV5ZFsCbfFhgH+Do5PJ7gUtIJODjks8BngBuNInvdMcBj1hrW4FvjTGLgZ2Tr1tsrf0GwBjzSPK1qyTgqgMuIk4oLi7ucLluPd6xPRf+l78/dDkF4RaiI0cSuv9+Fj/3Q0Zjim+7LaFnnsH7zDMUXHghngULKN5/f1qPPprmiy/GbuzdTq1lk7pqtl72DVtXfsPAZd+wfcXCtjKLbXEYw+I+W/Bx/yF83H87lvXsjS8awRcNkxdp5a+79YXWVm5742vyIq34ookR8z0/n8PkWQ8zcc7zeGN/oeWUUzb6wtVssGXldxz63guMmTeTkuYg9f4iXt1hFC8O3w/ssG5XanJps2YSZyvH5mQYY9wkpplsTWK0+mugzlq78gqCpUC/5PN+wBIAa23UGFNPYppKP2BOu822f8+S1ZbvkoZmiIisk0qedpK1HPvm4/zxxdtxWctLvxzNrs/en7xNfGYTcCBxkeb48dTvuy/5//oX+TfeSN5DD+GdPp2W886j9cQToTNTG4NB3J9/nrgJ0MrHZ5/xQn39z14ayivg0823Zcfxo4nuvDMHf+al0b/mGzuddXTig9vDsVU/3G33/edMmXEnIxd/CJddRt5tt9EydSqtxx8Pvu41YpoXbsH30EPk3Xcfj773Xtvy6qJSyoO1TJw7nYlzpxN76V+EjzyS8BFHEB8wIHMBp5BX+XfWciwBt9bGgGHGmB7A08A2Tu17pZqaGnbffXc8Hg+xWIwJEyYwZcoUKisrKSwsxO1209DQQHl5OcuXL8daS3l5OVVVVW13tguFQvTu3Zvq6mqMMfTq1Yvq6mqKi4uJxWI0NjbSp08fKisr8Xq9lJSUUFNTQ0lJCeFwmObm5rb1Pp+PoqIiamtr6dmzJ83NzbS0tLStz8/Px+/3s2LFCkpLSwkGg4TD4bb1fr8fn89HfX09ZWVl1NfXE4lE2tarTWqT2pSZNkWjUZqbm3/Wpm2KovTOizOvzsPwHlGaYoaKJjfbFkWpq6sjHA4zqizctj4YNSxrdjGoKMbXjW5KfXEqKiro06cPo8rC1EUMtWEXWxXG+DLopq8/TpHH0trampHjNKosvEqbvmtyU+SxlPoSba6oqMDn89E3L8agohg7P/YAR8x4EIA7DvwtXxx2BFvU1+NP/h7at2nl76Q27CIYNVRUVDh77p1xBsFRo+j/r39RMHs2Beedh/e++6jbYQc84TB54TDRujp8kQiEQtDYiLe1FUIh3GuoC19fUET15lvwbd8t+LrPFjRuM4g3irdkswCcMqI3ffr0YUTDAqpao0TisKk/zqcNHgYGYniM5ZOGxO80EAgwKBClb36cBfUeti+OEu21NRdufTXjl7zP8TMfIf/jjyk45xx8N97IDyedROjQQynp1Ssr/j/17t2bPUrDRK3hq5CbocVRlja72OqHrxnz7ovs+O4sClsaAWjKK2DBr/bi6V0OomaLrdnux2/oN2sWYz58jZJvv8V/5ZX4r7ySlpEjWTF2LDV77035wIFZ0+8NLY62/X8a3iNK3/wY9WHT1kdUVVWpL+9CbVobk5gd4ixjzEVAM3AO0Cc5yv0r4BJr7f7GmJeTz981xniASqAcOBfAWntFcjsv89NUlUustfsnl5/X/nUrvfPOO3bIkCFpb5+I5LbW1lby8vJ+trwzU1DW9prOvi5T01nWJ/Yj3nmas569kajLxUWTzmfmsL1/9prObMtx1uKdMQP/+efjrqjo3Ft8PmKDBhHbbjti226b+HfIEPadvmyNUyJSej6cOAzvCy/gv/xy3IsWARAbNIimadOI7rNPp9qQae3bt/cnb3L8G4+w3ZJFbcuiI0bQ+pvfcEDDlrT4/D97vyse49WBQXyPPopv+nRMUxMANi+P1uOPp/myyyAL7hWy+nEOeOKEoj8Ng+fyVLauaP78+fNGjx49oqN1TlVBKQci1to6Y4yfxMWSVwGzgMNIVEI5Hng2+Zbnkj+/m1z/urXWGmOeAx4yxlxL4iLMgcB7gAEGJquq/EDiQs2Vc8vb6CJMEXFCZWWlLsJci30XzGLqczcBcPlhU9uS7w3h+AcRY4gceCCRvfbC9+STmMZGbEEBtrAQGwjAyueFhdiCAggEEss7mqpiKn++LB2MITJ2LJEDD8T3+OPkX3kl7i+/JHDEETTdcAPho1a/hLOLspbfv3w3J76e+NYkmF/IizvtyzO7jOHmCw8DoGUN50Pc5Sa6995E996bplAI3/Tp+B59FM+bb5J/++24vvuOxrvugoICx5qTCsN76CLMbOXUFJS+wL3JeeAu4DFr7XRjzELgEWPM34EPgTuTr78TuD95keVyEgk11trPjDGPkbi4MgpMSU5twRhzGvAyiTKEd1lrP1s9CNXmFREn6K67a+aZNYu/PXolLmv595iTmT7igEyHtGH8fsLHHpvpKNaP20140iTCEyaQf/XV+K+9lsIpUzDBIK2nnJLp6NbKE41w4RPXMGb+q8SMi38fdApP7nowrb789d9YIJD4PUyahPv99wkcdRS+V17BNX48oUcewfbsmfoGpElTTHlNtnKqCsrHwM+GIpJVS3buYHkLcPgatnU5cHkHy18EXlxbHKqCIiJOKCkpyXQI62Vjp8Z0dqTZPX8+gd/8BhOL8uCeh+VsibyM8/loueACbM+eFFx4IQXnnosJBmk544wuWSXE1NVxwx3nMOKbj2jy5XPeMRfx321TU2chNnIkwRdfpGjiRDzvv0/RmDEEn3gC26/fut/cBVQ0Ka/JVjl1Zxrdil5EnFBTU0NhDpR8Wx/9f1xC4MipmMZGXtzp11x/0O+6ZLKXS1qnTMEGAhSccQb+yy/HNDTQfMklXeq4uJYsIXD44Yz45kuqi0r5y+S/s2jTQet+43qIDxpEw4wZFB1+OO4vvqD4gAMIPvEE8cGDU7qfddmQD7nbFkX5sVVTULJRTiXgGgEXESdk2wh4upXXV/PvO8/BtaKWyK9/zaWjz8K6VD+tKwgffzw2EKDw1FPJv+EGTDBI0z/+kZG7fq7OvWABgaOOwlVVxde9B3D6CdOo6tk7Lfuy/foRfPFFAkcdhWfuXIrGjCH0yCNpv/HSRrGW5VVBNl/eRHFTkOKmIL7HvsLU1WFWrEg86usxzc1Ed9+d8MEHY/v0yXTUkpRTCXgmKr6ISO4Jh8OZDqHLKGoK8u87z+MXK6qIjhhB6O67iT38ZabDknYiEycSKioi8NvfknfPPRAK0XTTTRmtCuJ55RUCJ5yAaWoi8n//x0mjpxJaSz30VLA9ehB88kkKTzoJ34wZFB16KKF77iG6775p3e+GGLH4Q85+5t9s8eP3q664u+PX+557Dv+55xLdZRci48YRHjs2a6bZdFc5lYDH4/FMhyAiOWBd9V9zRV64hWvvuYCtK7/lm002p9cjj+TE3RizUXS//Qg99hiBo48m74knMI2NNN55Z/KmSM7y3XUXBWefjYnHaZ00iabrriN038/qKqRHQQGN992H/fOfyXvoIQJHH52oFDNpkjP7X4eeoRX8efqtjJn/KgAt+X6qC3vSUFBEg7+Y4dtvTrxnT2xJCbZnz8QFpbEY3hkz8L72Gt45c/DOmUPBeecR3XlnwoccQviQQ7CbbprhluWenErAVZlARJzQR1/z4o5FueLByxj23adUlZTzp5Ou4v5evTIdlqxFdI89CD79NIHDD8f30kuYSZMIPfAABNI78tymuRn/lVeSf8MNiR/POouWc891fk66x0PTDTcQ32QT/NddR+Ef/oCprqb1j390No724nF8DzzA4/+4kJLmIC0eH3eNPpanRx9GHT/dc+BnF003JlfstR0Fu57KHp/P4eLgR3hffRXPe+/hee89Ci64gOjw4bScdhqRceMcbljuyqkEXHXARcQJOV8H3FrOffp69vx8DnUFRZx20lVU9dgkY+F02Zv6dEGx4cMJTp9O0YQJeN98k6Jx42idPJnoiBHEBw2CVM/dj0S4+Py72G/B6+z16TvktzYRdbmZNvEMni8/gFcydUGoMbRcdBG2vJyCv/6Vgosvxv3llzRdfTX4f36jn3TasvI7ig46H8/cuQC8O2gEV40/nR9Kf8GosjCzazq3nab8Al755T6cedJUCIXwzpyJ79ln8c6ciWfePApPPJFQr15E99wzja2RlXIqAXfpoh8RcYDPl9tVCcbPnc6h771IizePv0yexne9c/jDSBaKDxmSuCBx/Hg8H36I58PEB5h4SQmx4cOJjhyZeIwYAcXFG7CDOJ65c/E++SS+Z5/l37W1basWbjqYG8aczAdbd40PRa2nnkp8k00o/OMfyXvwQdyffELjPfcQHzAg7fvOC7dw0msPcOzsx/DEY8Q32YQLR5/MK8P2bvtWIBjdwA8ogQCR8eOJjB8PjY34r7iC/JtvpvDkk2mYPRvbOz0Xu8pPlICLiKdMpLQAACAASURBVKRYUVFRpkPImO2+/5yznr0RgMsnnsGn/YdkOCLndYcR9/iWWxJ89dXE3SI/+ADP++/jWrYM1+uv4339dQCsMcQHDSI6ciTxzTbDFhVhA4HEvysfK38uLsb93Xf4nnwS31NP4frhh7Z9fbvJ5rw8bB9e2XFvlpR3vbnIkYkTCQ4eTOHxx+P5+GOK9t6bpltvJbL//unZobXstug9zn7mBvotXwZA6+TJNF90Ea88/s0qL13WnIK8prCQ5r/9DffHH+N9+20KTzmF0FNPdYlKON1ZTiXgqgMuIk6ora0l4NS82S6kR6iOq+7/G95YlEd3O5QZO/060yHJRrDl5bSedhqtyZ/N0qV43n+/LSF3f/wx7kWLcC9atN7bjm26KZEJEwhPnMgRcyJdqvZ4R2JDhxJ8/XUKTj0V34wZBI46iuapUxNz1FOUqPpbm9l/wetMfPc5tvnfYgC+6rMFV0z8C/+8rOO7rg4qirGsNQX7d7tpvP12ivfaC+9bb5F/1VW0nH/+xm9X1iinEnCPJ6eaKyIZ0jOLbmWdMtEo0x76O73rq/m4/xCuG/v7TEckKWY33ZTIppsmpi0AtLbi/vhjPPPnY2pqMMHgzx+hUNtzW1hI+OCDCU+YQGznnX+aTz537d8YdBW2pITGBx4gdv315F9+Of5//hPPvHk0/uc/2NLSDd6ua+FC8u6+mxcfeJhAaxMAdQXF3LvXJB7ecyIx95pzl68bUzdKbfv0ofH22wmMH0/+P/9JdNddie6zT8q2L6vKqYxUZQhFxAnNzc0Ub8jc2CyWP20aIxd/SG2gB+ceexFRj6pOdXt5ecRGjlzjzWrWOBXnU+DTj7JiKs7PuFy0/OUvRHfaicKTT8b7xhsU77UXoXvuITZ8eOe309qK77nnyLvrrraLKwEWDBjKU7uO5bXtRxH2rvtaklJfnKXNqUvCo//3f7Sccw7+K6+k8Pe/p+GNN7C/+EXKtt8ZG3JH0GykBFxEJMVaWloyHYKjvNOn47/uOqIuF+cfcyHVJeWZDkkkraKjRtEwaxaByZPxfPABRWPG0HTFFYQnT4Z4HBobMU1NiUdz8yo/e+bOxffQQ7iSF5/aQIDWI49kcvEuLO671XrF0cOb+hsMtkydimfOHLxvvEHhyScTevZZ0AyClMup36jqgIuIE7pKHXAnLgbcvHoJhZeeBsCNB57M/K2GbfQ2RbKB7deP4PTp+C+8kPz//IfCM8+k4PzzMZ28E250++1pnTyZ8MSJUFTE4nX8f+3IvLo0pHFuN4233ZaYD/7uu+RPm0bLRRelfj85LqcScNUBFxEn5EodcH9rM1ffdwkmFCJ8yCE8uPvhmQ5JxFk+H81XXUV05EgKzjwTV0MD1hgoLMQWFGD9figoSDwvLMT6/dhf/ILWo45KTFnZyItPh/eIMrsm9WVPbXk5jf/5D4Fx4xLfbu26K9H99kv5fnJZTiXgKkMoIk7Iz8Dtux1nLRc8cQ1bVX1HbOBAGm+4AR5dnOmoRDIicthh1I8bB7EY5OU5VtWlLpK+/UR3243mv/6VgksvpfDUUxP1wXXL+pTJqYxUCbiIOMHv8J3yMmHSO0+x30dv0OjzE7rvPsjh2uciAHi9kJ/vaEnF2nB685rWP/2JyL774lqxgsAJJ0Anp9fIuuVURqo64CLihBUrVmQ6hLQa9u3HnD79NgAuPeIs4oMHZzgikdy0VWEsvTtwuWi85Rbi/frh+eAD/Jdemt795ZCcmoKiOuAi4oTSjagJ3NX1XV7JFQ9chice44H/O5zXdxjFuZkOSmQ9dIc7la70ZTD9d6u0vXoRuvNOisaOJf/mm4n36UPrlCld/uZJXV1OZaQqQygiTggGg93yTpgDqiq48Y6zKQsuZ96WO3LjgSdnOiTJkO6UxGazvv54au6EuQ6xnXem6corE5VeLroI99df03T11YlpN7JBlICLiKRYuBvOkxy89EtuuPNcejbWM3+L7Zl6/KXEUnQLbpHuLJ0fVoo8qa8DvibhE07A9uxJ4ZQp5N17L67vvqPx7ruxPXo4FkN3klMJuOqAi4gTukod8FTxvPsut95+JoGWRv47eCRnH3cJrb4cqPQi0sWlpQ74WkTGjye42WYEjj0W7+zZFO2/P6FHHiG+xRaOxkE8juuLL/D+97943n0XEwoRGziQ2ODBxAYPJj54MLakxNmY1lNOJeCqAy4iTuhOdcA9r75K4PjjMS3NzNxhFBdNOk+3mRfpItJVB3xtYiNGEJw5k8KjjsKzcCFF++5L4333Ed1tt7Tt0x2LMviHxeTd8Daed9/FM2cOrrq6VV7jnTlzlZ/jffsSGzQokZRvsw2xIUOIjRzZZeau51QCrjKEIuKE7lKGcJ+PZxP46xWYSIRnRx7AtIlnEHfl9rQTzX2WriTdZQjXJL7ZZgRffJHAySfjnTmTwPjxNF1/PeFJk1K2j+KmBibMeZ4RXy9gh+8W4o+0rBpDv35EdtuN6K9+hS0txf3ll7i+/BL3okWJ58uW4Vq2DO/s2W3vCR96KI233go+Zz+0dCSnEnDTRT71iEj35usCnfvGOvj9Gfz1iX9ibJyWU0/l7wMmdpmRI9kw+vDQ/QSjGfw/WVxM6KGH8F9wAfm33UbhH/6Aa/FiWs4/HzZiwNMdizFh7vOc8so99GgKti2vKNuUPgfuTXS33Zj0fQmVPXsnVkSBKqBkUxi5D4yEVybvgGvJkkRC/sUXuBctwvf88/ieeQYTChG65x4oKNio5m+snErAY7E018sUEQHq6+vpkcUXJk16+0mmPnczAM3nnEPL2WfDnQsyHJWIrG5AQYyKpgx+K+V203zFFcQGDqTgnHPwX3st7q++ovGWWzYowfXMmsWD101lq6rvAHh/q2E8+atD+HCL7Vle1KvtQ2LlOj5M4nYTHzCA+IABRPfbD4DWk08mcPjheF99lcDhhxN6+GEoLl7vGFMlp+ZkqA64iDihrKws0yFsGGs5aeZ9bcn3tQf/gZZzztHIt0gX9Xmwa+Q14RNOIPTII9iiInzPP0/J0KEUTJmCd8YMaGlZ5/tdX39N4dFHUzRxIltVfcfSXn056zd/4w+nXMNrO4xieVGvjY4xNmwYwRdeIN63L95336Vo3DhMTc1Gb3dD5VQCrhFwEXFCfX19pkPYIH+YcSe/m3kvMePi0sPP5OE9J2Y6JBFZi/4FXSeviY4eTcPLLxPdfntcdXXkPfwwgaOPpsegQRSedBLeZ5+FxsZV39TQgP/CCynebTd8M2ZgAwFuOPAkjpx6F28M3SPlH/7jgwYRfOklYltuieejjyg66CDMDz+kdB+d1TU+OjnEWufqZYpI7srGiktHvPM0k2c9TNTl5oKj/8prO4zKdEgisg4F7q6V18S32Ybg7Nm4knOuvdOn4/n4Y3xPPYXvqaewfj+R0aOJjB0LTU34p03DVVODNYbWY46h+YILuO/5/6U3xs03J/jCCwQmTkxUcRkzhtDTTxPfcsu07nd1OZWAqw64iKyvtV24tqaL1rKtDvjen7zJ1OduAuDyw85Q8i2SJZyuA95Z8cGDaRk8mJYzz8T13Xd4n38e3/PP4/ngA3zTp+ObPr3ttZFdd6V52jRiw4Yll6Q3AQewvXsTev55AkceieeDDygaM4bgU08RHzIk7fteKaemoGTjqJSIZJ/KyspMh9Bp7jlzuOzhabis5eb9T2D6iAMyHZKIdNLwHtFMh7BO8QEDaP3jHwm+8gp1n3xC05VXEtljD6LDhhG64w5CL7zQLvl2ju3Zk+BTTxEZNQrXjz9SNHYs7vffd2z/OZWAu3XbZBFxQGFhYaZD6BTXokUEjj6avGiEJ3cZy937HJ3pkERkPVS1ZlcaZ/v1o/WUUwg99xzB118nMmFCZi/yDgQIPfww4TFjcNXVUTRhAp52dcPTKbuOnIhIFsiGD/ulDbUEjjgCV10db277K/5x6J9U7UQky0TimY6gG8jPp/Gee2g98khMYyOB447DtWRJ2nebUwm4qqCIiBMaGhoyHcJaFbY0cv1d5+NesoTo8OGcf8wFxLLgQ4OIrGpTvzLwlPB4aLrpJsJjx2JCIQr++EeIp/d3m1MJuC7CFBEnlJeXZzqENfJEI1x5/98Y/L/FxLbcktDDD9Pqy890WCKyAT5t6JoXYWYll4uma68lXlaG9803ybvzzvTuLq1b72Ki0a5/sYKIZL/ly5dnOoSOWctfn/wnu341j9pAD0KPP47N1psGiQgDA/pmP5VsWRlN//wnAP6//Q3X11+nbV/66CQikmJd9Z4Dp758F2PnzaTZm88Zky/nui22yHRIIrIRPKZr9jVdxYaUkY0cfDCthx9O3uOPUzhlCsEXXoA0TNHLqQRct6IXESd0xSkoE959nhNef4ioy8W5x17Ews22yXRIIrKRPtmIKShrS05hzQlqLmi+6iq8b7+N5733yLv5Zlr/+MeU78ORKSjGmM2MMbOMMQuNMZ8ZY05PLr/EGPODMWZB8jGm3XvOM8YsNsYsMsbs3275Aclli40x57ZbvoUxZm5y+aPGGN/qcagOuIg4oaqqKtMhrGLsBzM455nrAbhiwl/477a7ZDgiEUmFYSWaWpsOtkcPGq+7DiBxt84vvkj5PpwaEo4CU621840xRcA8Y8zM5Lp/WWuvaf9iY8wQYBKwHfAL4FVjzKDk6puAfYGlwPvGmOestQuBq5LbesQYcytwInBL++1mQ2kwEcl+gUAg0yG0GT9nOuc/9S8AbjrgBJ7becw63pG7NuTr6u5Iv4fssawlpy7lc1R0331pPe448u6/PzEVZcYMSGExD0cScGvtMmBZ8nnQGPM50G8tbxkHPGKtbQW+NcYsBnZOrltsrf0GwBjzCDAuub19gJV3kbgXuITVEnARkXRYPWEZFIjyZej7tp8zlbQc8c7TnPXsjQBcP+YUHtjryIzEISKSjZouuwzPrFl4PvyQ/Ouuo+Wss1K2bccnRRtjBgC/BOYCuwOnGWN+A3xAYpR8BYnkfE67ty3lp4R9yWrLdwFKgTprbbSD17epqalh9913x+PxEIvFmDBhAlOmTKGyspLCwkLcbjcNDQ2Ul5ezfPlyrLWUl5dTVVXVNqIVCoXo3bs31dXVGGPo1asX1dXVFBcXE4vFaGxspE+fPlRWVuL1eikpKaGmpoaSkhLC4TDNzc1t630+H0VFRdTW1tKzZ0+am5tpaWlpW5+fn4/f72fFihWUlpYSDAYJh8Nt6/1+Pz6fj/r6esrKyqivrycSibStV5vUJrVp49sU8MQZ3iNKVauLSDxRd/fTBg8DAzGWLFlCeXk5o8rCbSNRffPj5Lksm+TFiVrDVyE3FRUVFBcXs01RlN55cebVeRjeI0pTzFDR5Gbboih1dXWEw2FGlYXb1gejhmXNLgYVxfi60U2pL05FRQV9+vRhVFmYuoihNuxiq8IYXwbd9PXHKfJYWltbiV19dVvy/cxRv+eVXScysiBCgdu2bb+mpga3282osnBbmzzG8kmDh2ElUZa1uKitrSUUClHijbN9cbStTUOLoyxtduF1sUpM7dv0XZObIo+l1Jdoc0VFBT6fj755sVXa1MP7U0xVVVX4/X5GlYVXadPK9bVhF8GooaKigrKyMkb2XLVN7Y9TU1MTy5cvZ4/S8CptWnmcFtR72L44ytKlS+nVqxejysJtbVr9ODU2NlJTU0P/gtgqbVr9ODU0NLBjSWSVNq1+nJYtW0Y4HG47t1a2aUBBjM+DHvoXxFb5nXZ07q08ThUVFQQCAQYFoqu0qf1xWrFiBbFYbJVzq6Pj1H6fHZ17PbyJc6uyspIdSyIdnnsrt//jjz/SvyC2SptWP04rVqygoaGBUl+8w3MPEudW79692aM03OG5t/I4VVRU4PV62SQv3uG5N7xHlGXLllFUVMSosnCH597K41RRUZHSfq/EG1/juRe1hqamJqqrq9mqMNrhuVfRlOhDSkpKGFocXaVNffNjBCOm7ThVVVXR0tLSdm511Ee0P84dnXsrj1NFRUVG+vKdekQ6PPdWxpSfn8+m/tgaz73asIu6ujrq6+vZJC/e4bkXiSfOrc60aflFF7HVKaeQ/49/8P2OO+IdObLTbVprPuzk1frGmAAwG7jcWvuUMaY3UANY4DKgr7X2BGPMjcAca+0DyffdCbyU3MwB1tqTksuPI5GAX5J8/dbJ5ZsBL1lrh7bf/zvvvGOHDBmS7maKSDfSma/jV39NiTdOfcS1ztd1Zlsb8ro3617DP20aANMm/Jmndz14g7fldOzdYVvZHHsqt5XNsadSOmN3oq9xmtPnQ2f4zzmH/P/8h+iQIQRfew3y8jr1vvnz588bPXr0iI7WOTZ5yBjjBZ4EHrTWPgVgra2y1sastXHgP/w0zeQHYLN2b980uWxNy2uBHsYYz2rLV6E64CLihO2LM9TXWMvvX74b/7RpWGO49PCz1ph8i0j2y1hfk2OaL7qI2JZb4lm4kPx//CMl23SqCooB7gQ+t9Ze225533YvGw98mnz+HDDJGJNnjNkCGAi8B7wPDExWPPGRuFDzOZsYxp8FHJZ8//HAs+lsk4jImkStcX6n1vLHF2/nxNcewLrdNN52G8+PPMD5OETEMRnpa3JRYSGNN96INYb8667D/cEHG71Jp0bAdweOA/ZZreTg1caYT4wxHwN7A38BsNZ+BjwGLARmAFOSI+VR4DTgZeBz4LHkawHOAc5IXrBZSiLhX4XqgIuIE74KOVxxyVqmPncTv5n9GFGXm8Y77iBy2GHrfp+IZDXH+5ocFtt1V1pPOw0Tj1M4ZQqEQhu1PaeqoLwNdPQx7cW1vOdy4PIOlr/Y0fuSlVF2Xn15e6oDLiJOGFocZXbNz25FkBYmHuecp69n4tzphN1ezjv2Ii4YN86RfYtIZjnZ1wg0n3ce3ldewb1oEcV7703jbbcR22mnDdpWpxJwY8xRwAJr7efGmMEk5mvHgFOttamvTp4mqgMuIk5Y2uxcbd7fznqYiXOn0+Lxcfbxf+PdwWsdhxCRFMr0hYxO9jUC5OcTuv9+Cn/7WzwLF1J0wAG0nHMOLX/+83rfrr6zI+B/B3ZLPr+GxHzsEHAzifrbIiKS5HXob2KgOcRxsx8F4NzjLlbyLZJjnOprUiHTH1ZSJb711gRffRX/ZZeRf8st+C+/HO+rr9J4663E+/fv9HY6e+jKrbVVxph8YA/gr8ClwLD1Dz1zYrFYpkMQkRzQOy/uyH4mvf0URS2NvL/1L3ln210d2aeIdB1O9TWymvx8mi+/nOCTTxLv0wfP3LkU77knvkcfhU6W9+5sAl5tjNkaOBB4P3mHynw6ntfdZXlTeAtREZE1mVeX/strClsaOertJwH4z69/k/b9iUjX40RfI2sW3XtvGt5+m/DYsZhQiMJTT6XwpJMwdXXrfG9nj9xlwDwS875X3sv418BHGxRxhugiTBFxwvAe6b8w6vD/PkNxc4h5W+7Ih1vukNZ9iUjX5ERf0xndZXrJhrC9etF4771EHnyQgvPOw/f003jee4/GW26BgoI1vq9TI+DW2nuAvsCm1tqZycVzSNThzhqJcuQiIunVFEtzXxMMcsybjwNwx6+PS+++RKTLSntfI51jDOFjj6Vh9myiw4fj+uEHAuuoRrU+0/f9wERjzNnJnz04VMYwVVQFRUScUNGU3r4m7+676dEUZMGAoXywVVZdiiMiKZTuvkbWT3zLLQm++CLNZ50F6xj07VQCbowZBSwCjgEuTC4eCNyyMYE6TbeiFxEnbFuUxr6msZH8G28EkqPf+mZPJGelta+RDeP10nLeeQRfeGGtL+vsCPh1wJHW2gOAlUd7Luu48U1XoxFwEXHCd2kclcq75x5cNTV8svm2zB04PG37EZGuL519jWyc2C67rHV9ZxPwAdba15LPV9ZXCZNlU1BsJ0vDiIhsjCJPevqavHAL+TfcAMAdozX6LZLr0tXXSPp1NgFfaIzZf7VlvwY+SXE8aRWPq16miKRfqS89fc34917A9eOPRIcN47/bZNUXkCKSBunqayT9OjuCPRWYbox5AfAbY24DDgbWfolnF6M64CLihHTU5vVFwvxmVuKuly1nnQU/aPRbJNepDnj26mwZwjnADsBnwF3At8DO1tr30xhbyqkOuIg4YXiP1F8YNe79FykP1hIdOpTIAQekfPsikn3S0deIMzr90cla+z/gagBjjB/Iuu89XK71qbooIrJhgtHUjk57o2GOn/UIkBz91txvESH1fY04p7NlCK8xxuycfH4QsBxYYYw5OJ3BpZoScBFxwrLm1PY1Yz94md711SzuswWRgw5K6bZFJHuluq8R53T2yB0DfJp8fhFwLHAIMC0dQaWL6oCLiBMGFcVSti1PNMLk1x8C4M7Rx4IGEkQkKZV9jTirsz15gbW2yRhTCmxprX3SWvsq0D+NsaWcx6OLFUQk/b5uTF1t3jHzZ9K37ke+2WRzXt9+z5RtV0SyXyr7GnFWZxPwL40xxwCnATMBjDFlQHO6AksHlSEUESekqjSYOxblhOTo912jjyXu0h9bEfmJyhBmr84OCf8BuJ7EzXdOTC7bH3glHUGlixJwEXFCD29qbo5x4Iev0m/5MirKNmXmjnulZJsi0jn73fHhGte9ctIvHYxkzVLV14jzOpWAJ8sN7rbasgeBB9MRVLqoDriIOCEVtXndsRiTX1s5+n2MRr9F5GecqAOeDR9EslGnj5wxxgcMBsqAtro31trX0xBXWqgOuIg4YXiPKLNrfBu1jX0+mc3mtT+wtFdfXh42OkWRiUh3koq+RjKjUwm4MWYP4HEgDygGGoAiYAmwZdqiSzGVIRQRJ9RFNrI2r7VtlU/u22sSMbdGv0Xk5za6r5GM6WxG+i/gamttLyCY/Pcy4Oa0RZYGSsBFxAm14Y3ra/b4fA4DK7+luqiU6SP2S1FUItLdbGxfI5nT2SM3iMRFmO1dCfwlteGkl+qAi4gTtirciNq81jL59cTlNQ+MOpyIR18vi0jHNqqvkYzq7BzwehJTT+qAZcaYIUAtEEhXYOmgOuAi4oQvgxs+ZWSnbz5ih+8/p66giKd3GZvCqESku9mYvkZSb/ULVq/cac2v7ewI+FPAmOTzu4BZwDzgifWOLoNUhlBEnNDXv+F9zcq534/uPoHmPH+qQhKRbmhj+hrJrM6WIfxzu+fXGGPmkLgI8+V0BZYOSsBFxAlFng2rzeueP59dv5pHo8/Po7uPT3FUItLdbGhfI5m31hFwY4zfGDN09eXW2rdJVEDJqsmJqgMuIk7Y0Nq8+dddB8ATux1CsKAolSGJSDfkRB1wSY91TUE5m5/ufLm6ycBZqQ0nvVQHXEScMLzH+l/w7friC3zTp9Pq8fLwnoelISoR6W42pK+RrmFdCfiRwDVrWHctcFRqw0kvlSEUESdsSGmw/OsThaaeG3kgtUW9Uh2SiHRDKkOYvdZ15PpZa3/oaEVyeb/Uh5Q+xqhgvYikXzC6fn3NL5Yvw/fEE1i3m/tHHZmmqESku1nfvka6jnVNHmo0xmxmrV2y+gpjzOZAU3rCSo9YTPUyRSRh9XJRq3vlpF9u8LYHFMSoaOp8ebDjZj+GicVonTSJZb36bPB+RSS3rG9fI13HukbAXwSmrWHdZcALqQ0nvVQHXESc8Hmw831NaUMtB7//EtYYWk4/PY1RiUh3sz59jXQt6zpyFwDvGmM+IlELfBnQFxhP4sY8u6U3vNTSCLiIOKF/QYwfWzs3N/Pot54gLxohPHYs8cGD4a21j8yLiKy0Pn2NdC1rTcCttZXGmJ2AqcABQCmJO2A+D1xrrV2R/hBTx1rVyxSR9Ctwd66vKW5qYOK7zwPQ8pe/pDMkEemGOtvXSNezzu8ukkn2BclHVlMdcBFxQmdr8x7xzjMUhpt5d9AItvnlhs85F5HcpDrg2cuR7y2MMZsZY2YZYxYaYz4zxpyeXN7LGDPTGPNV8t+eyeXGGPNvY8xiY8zHyVH4lds6Pvn6r4wxx7dbPtwY80nyPf82HZQ8UR1wEXFCZ2rz+lubOfKdpwC4Z++j0x2SiHRDqgOevZyaOBQFplprhwC7AlOMMUOAc4HXrLUDgdeSPwMcCAxMPk4BboFEwg5cDOwC7AxcvDJpT77m5HbvO2D1INxuXSksIulX1Yk5mePnTqdHU5CP+m/H/C13cCAqEeluOtPXSNfkyJGz1i6z1s5PPg8Cn5OoIT4OuDf5snuBQ5PPxwH32YQ5QA9jTF9gf2CmtXZ5cmrMTOCA5Lpia+0cm5jofV+7bYmIOCoSX/t6bzTMsbMfB+DufY4G3aNARDbAuvoa6boc/+hkjBkA/BKYC/S21i5LrqoEeief9wPa1x5fmly2tuVLO1i+ClVBEREnbOpf+1/Fg99/mfJgLV/23ZJ3ttnFoahEpLtZV18jXVenZu8bY2qBN4HZyccCuwElRYwxAeBJ4M/W2ob207SttdYYk9bLeevr69l9993xeDzEYjEmTJjAlClTqKyspLCwELfbTUNDA+Xl5SxfvhxrLeXl5VRVVREIBAAIhUL07t2b6upqjDH06tWL6upqiouLicViNDY20qdPHyorK/F6vZSUlFBTU0NJSQnhcJjm5ua29T6fj6KiImpra+nZsyfNzc20tLS0rc/Pz8fv97NixQpKS0sJBoOEw+G29X6/H5/PR319PWVlZdTX1xOJRNrWq01qk9q05jbt1CNCkccyr87D8B5RasMuglHDgIIYnwc9/O9//yMSiRDwxBneI0pVq4tIPPEH79MGDwMDMZYsWUJ5eTmjysIsa0mMZ/TNj1PR5GKP0jBRa/gq5KaiooLi4mK2KYrSz9XCSa8/AMDsQyaxSb5l26IIdXV1hMNhRpWF22IKRg3Lml0MKorxdaObUl+ciooK+vTpw6iyMHURQ23YxVaFMb4Muunrj1PksbS26K3vPQAAIABJREFUtlJZWcnQ4ugqbepfEKPA/VOba2pqcLvdjCoLt7XJYyyfNHgYVhJlWYuL2tpaQqEQJd442xdH29o0tDjK0mYXXherxNQUM1Q0udm2KMp3TW6KPJZSX5x5dR4qKirw+Xz0zYut0qYe3p9iqqqqwu/3M6osvEqbVj9OFRUVlJWVMbJnZJU2tT9OTU1NLF++nD1Kw6u0aeVxWlDvYfviKEuXLqVXr16MKgu3tal3Xrxtm00xQ2NjIzU1NfQviK3SptWPU0NDAzuWRFZp0+rHadmyZYTD4bZza/Vzr39BbJXfaUfn3srjVFFRQSAQYFAgukqb2h+nFStWEIvFVjm3OjpO7ffZ0bnXw/vTubVjSaTDc2/l9n/88Uf6F8TWeO5VtbpYsWIFDQ0NlPriHZ57kDi3evfuvcr/p/bn3srjVFFRgdfrZZO8eIfn3vAeUZYtW0ZRURGjysIdnnsrj1NFRcU6+4jGxsb/Z+++w+OorsaPf4+0WpVVL5ZtbGQwmEBMNTVADC9gDIFQQktCKIGE3kliorwkoSckhBJCCfiFJPRO6A4/MDEdTDPN4MSLZcvqvW27vz9mVt6VdmUZa3c02vN5Hj2WZ2Znzsy9e/ZqduYMHR0dccc0WY4oyYkk7XshI/T29tLU1MRMXyhh3/P3WjmkpKSE2cWhuH0yhrj3U0NDA/39/YN9K1GOiG3nRH0v2k5+vx+fz8dMXyhh39uhJERLSwsAcysDCftebI6IxpQsR0T71k6lwYR9LxpTXl4e0/LDSfteSyCL9vZ2Ojo6mJQbSdj3ghErprH+fJpbGYjbJ0he/ENGM44WkZnAt4G59k8psARYbIz5w3pXYK0jB3gKeN4Yc5097XNgH2NMvX0ZycvGmK1E5Db79/til4v+GGNOs6ffBrxs/7xkjPmGPf37sctFLV682Gy3nV5rqZQa/ZMwR1ou2TK7lwd5ozUn4XLHLnmUi5+8meVTNuf4827DZGWNensbG1eq1uXm2FOxLjfHPpbrcnPsY7kup3LNxmwvFesar204loZu75qdzLv77bffzomWHdUZcGPMCmAF8H8iMgs4ETgb65rs9Q7A7YokdwKfRgfftiftdV1j//tEzPSzReR+rBsuO+xB+vPAVTE3Xs4DLjHGtIpIp4jsjnVpywnATaPZN6WUu6Q7oX4dniRf5uUG+jn5/90LwG3zThocfCul1NeRLNeo0RvtYH6sjfYSlDOwzoDvCazBuhzl+1hnwUdjT+BHwEci8r497ZdYA+8HReQUwA8cY897BjgY+BLoBU4GsAfalwNv28tdZoxptX8/E7gLyAeetX/i6KPolVLp8FFn4lxz9OtPUtHdxifTtuKVbVz1IGGl1DiULNeo8W+0LXcz1hnwy4GnYm6cHBVjzBIg2W3++yVY3gBnJVnXQmBhgunvALNHikPrgCul0mGHkhCLm73xE7u6OOHl+wG4bd6JWvlEKbXREuYa5QqjHYBPw7r2+9vAefb13K9gXQP+j1QFN9a0DrhSmcGprxSjojdbxcr761/J7+ngw5pteG2rXVO6faVUZkiUa5Q7jKrljDFrjDH3GWPOwLqU5GGsy0XuHvmVSimlpKOD3Jus21JunXeynv1WSqkMN6oBuIhcICJP2OUInwA2BS4Gtk5lcGNN64ArpdJhSl58bd7cv/yFrI4O3tl8e97eYnzcKKqUcr+huUa5x2gvQdkeeBy4wBjznxTGk1I5OcnrMSql1Fh5v2Ndai3p6SDvllsAuPVAPfutlBo7sblGuctoL0E5yRjzf8BKu163K4VCIadDUEplgG2L1+Wa4195COnuJrjvvnyw2bYORqWUmmhic41yl9FeglIqIvcC/VilARGR74rIFakMTiml3ChkrLPcZd1tHLvkMQD6fvlLJ0NSSk1A0Vyj3Ge0t8/eCnQANUDAnvY6cGwqgkoVrQOulEqHL7qtiksnvnQ/+cF+AgceSHjOHIejUkpNNNFco9xntCPS/YCpxpigiPXYJWNMk4hMSl1oY0/rgCul0mF2cYiPV3TyvdefBKD/kkscjkgpNRHNLtY64Oky1k9hHu0Z8A6gMnaCiGwKbNADeZymdcCVUulQ15fFSS/dR14owP+bvTfh7bZzOiSl1ARU16d1wN1qtC13B/CIiOwLZInIHlg1wG9NWWRKKeVSk9sbOeLNp4mIcPu8E50ORyk1QeXo+Nu1Rtt0vwMewHokfQ7Wo+CfAG5IUVwpoXXAlVLpcNzz9+ANB3lh+31ZMXkzp8NRSk1Q1blaB9ytRnUNuDHGYA22XTXgHkrrgCulUm2TljXsvGQRYcnijv1PcDocpdQE9m67FpdwqxFbTkS+vb4VGGNeGbtwUktvwlTKGWN988p4duq//oYnEuapOfPwT5rudDhKqQlsTqnehOlW6/vT6Z4k0w1QCvgA19zZKPoEOqVUCm238mMOeXcRoWwPd+z/I6fDUUpNcL1hHde41YgDcGPMsNM3dunBWuBkXHYTplZBUWr8GuksOYz/M+XZ4RALHv0TAI/8z9GsrpjqcERKqYnO36vjGrca9f2z9tMwrwaWA0XAdsaYs1IWWQroo+iVUqny/SWPsuXa/1JXPoX3vueqZ5QppVxq6yId17jVeq/eFxEfcBFwPrAI2M0Y83mqA0sFPQOulEqFyW0N/PSFuwH4/RHnsiZU4HBESqlMsFLPgLvWiGfAReRi4L/AHGAfY8yxbh18A1jFXJRSamxd/MSfyQ/2s2i7fXh9q10p8miuUUqlnuYa91rfGfDfA61AOfDnRDcxGmPWWyllvIhEtF6mUmps5Tz9NHM/eY3uPB/XHXoGABVezTVKqdTTXONe6xuAn5yWKNJE64ArpcZUdzcFv/gFAH858Mc0l1QCWptXKZUemmvca31VUO5OVyDpoHXAlVJjKf+aa8has4ZPpm3FI3scOjhda/MqpdJBc417jboKykSQlZVRu6uUSqFZa74k97bbMFlZXH3k+USy1t0M1RXS2rxKqdTTXONeGTUi1QG4UmosZEXCLHj0eiQcZuDUU/ls2qy4+fV9mmuUUqmnuca9MqrltA64UmosHP7mM2z71adEpkyh75e/HDZ/VlHYgaiUUplGc417jWoALiIVqQ4kHTwevVlBKbVxKrpaOfvZvwLQe9VVUFw8bJkVPVqbVymVeppr3Gu0Z8C/EpEnROQoEXHt1f5ahlAptbHO/+etFPX3sOQbuxH87ncTLqOlwZRS6aC5xr1GOwCfAbwI/AJYKyK3i8heKYsqRXQArpTaGLsuf4f5779If04u1x5+DiR4NgJAaY4+HEMplXqaa9xrVANwY0yTMeZGY8wuwB5AI/B3EfmPiFwmIjUpjXKMaB1wpdTX5Q0G+MVjNwJwx/4/Yk35lKTLam1epVQ6aK5xr69zE+Zk+6cYWAFsArwnIgvGMrBU0DrgSqmv65jXHmfTltX8p7qGe/Y+asRl55TqDd9KqdTTXONeo/rTSUS+CRwP/ADoAe4GtjfG1NnzLwc+BK5JUZxjQssQKqW+Dl9/Dye+dC8A1x9yOiHPyN+mtQe1Nq9SKvU017jXaL+7eAW4DzjaGPPW0JnGmJUicv2YRpYCOgBXSn0dP3zlIUp7u1i62ba8PmuX9S7fEtBco5RKPc017jXaAfhkY8yI128YYy4dg3hSSuuAK6U2VGl3Oz945WEA/jL/lKQ3Xsaa6QtT16flwZRSqaW5xr2SDsBF5MdD/p9wOWPMwjGOKWW0DrhSakOd9NJ9+AJ9LPnGbnyw2bajes3yLv1AVEqlnuYa9xppRPqjUbzeAK4ZgGsZQqXUhpDVqznq9ScAuOXAk0f9uin5EeoH9INRKZVammvcK+kA3BizbzoDSQcdgCulNkT+tdeSGwqyaLt9WL7JlqN+XZFHa/MqpVJPc417jXQJihhjjP170qv8jTGuGdVqHXCl1GhlrViB9557CGVlceuBJ23Qa7U2r1IqHTTXuNdIt892xPweAoJDfqLT1ktEFopIo4gsi5n2GxFZLSLv2z8Hx8y7RES+FJHPReTAmOnz7WlfxtYdF5HNRORNe/oDIuJNFIfWAVdKjVb+Ndcg4TBPzTmQr6qmb9BrtTavUiodNNe410gD8G/G/L4ZsPmQn+i00bgLmJ9g+p+MMTvYP88AiMg2wHH29ucDfxGRbBHJBm4GDgK2Ab5vLwvwO3tdWwBtwCmJgtAyhEqp0dhyzQq8jzyC8Xq5Y//R3A4TT0uDKaXSQXONe410DfiqmN/9G7MRY8wrIjJjlIsfBtxvjBkA/isiXwK72vO+NMb8B0BE7gcOE5FPgf/BekgQWA8J+g1wy9AVJ6vkopRSsc543rq3fODkk2koq97g13eFNNcopVJPc417jfpPJxH5roj8UUTuFpG/RX82cvtni8iH9iUqZfa0TYBVMcvU2dOSTa8A2o0xoSHThwmHwxsZrlJqottu5cfs/ekbGJ+P/gsv/FrrmFGguUYplXqaa9xrtI+i/zVwOnA/cDRwG9YZ5wc2Ytu3AJdjlTK8HPgj8OMRX7GR2tvb2XPPPfF4PITDYY488kjOOuss1q5di8/nIzs7m87OTqqqqmhtbcUYQ1VVFQ0NDRQWFgLQ3d1NdXU1TU1NiAjl5eU0NTVRXFxMOBymp6eHyZMns3btWnJycigpKaG5uZmSkhICgQB9fX2D871eL0VFRbS0tFBWVkZfXx/9/f2D8/Py8sjPz6etrY2Kigq6uroIBAKD8/Pz8/F6vXR0dFBZWUlHRwfBYHBwvu6T7tN42afdy4N4xPBRp4cdSkLU91t/+0/Ji9Df309TUxO7lwf5ojub2cUh6vqyyMmC6twI77Z7mFMaYs2aNZSUlDC3MsDK3myKPIYK77r5XSGhu7ublpYWpuWHqfBGKM0xg/Pbg0JLIIuZvjDd3d3sVBqkyLNufksgi64g/PLFvwLQfPzxrO7tpdATYU5piIaBLIIRmJYfYVmnhy0Lw6xatYqqqirmVgbi9qmuL4u9KgKEjPBFdzZ+v5/i4mK+URSK26fesODvzWbrohDt7e0EAgHmVgbi9qm+L4tZRWFW9GRT4Y3g9/uZPHkycysDcfu0vCubKfkRijyGgYEB1q5dy+xiax0zCsJ82uWhpiBMQfa6fW5ubiY7O5u5lYHBfRraTi0tLXR3d1OSE2Hb4tDgPsW2U2xMsfs0tJ38fj9er5cpueG4fYptp4aGBvLz85lbGYjbp7h2Cgl+v5/Kykp2KQvG7VNsO/X29tLa2speFYGEfe/9Dg/bFoeoq6ujvLycuZWBhH2vNyz09PTQ3NxMTUE4Yd+LtlNnZyfblwST9r3lXdnU19cTCAQG+1Z0n2LbKfaYJup70Xby+/0UFhYyqzAUt0+x7dTW1kY4HI7rW4naKXabifpeac66vrV9STBh34uuv7GxkZqCcNK+1zCQRVtbG52dnVR4Iwn7Hlh9q7q6Ou79lChH+P1+cnJymJQbSdj35pSGqK+vp6ioiLmVgYR9L9pOfr+fioqKxDnCbqeenh46OjrijmmyHFGSE0na90JG6O3tpampiZm+UMK+5++1ckhJSQmzi0Nx+wTEvZ8aGhro7+8f7FuJckRsOyfqe9F28vv9+Hw+ZvpCCfveDiUhWlpaAJhbGUjY92JzRDSmZDki2rd2Kg0m7HvRmPLy8piWH07a91oCWbS3t9PR0cGk3EjCvheMWDGt7/Pp/Q4Pq1atQkSo8EY2+PMJkhf/GO3tsz8GDjDGLBORk40xF4jIfcCvRvn6YYwxDdHfReSvwFP2f1cDsXc8TbOnkWR6C1AqIh77LHjs8nHKysp49dVXh02vqamJWwagoKAg4fyKigorkOnTE86vrKwcNs3n8424zegAq7i4OOH86PTockPnl5aWDtuO7pPu03japzda1yWhxc3r7pFe3g21eXlMnz6dN55vHjb/s651r6mdOnXY/KHrLCwspLCwkLpFrXFPh4t9TV1fNrWFhSxtHx7TbsvfYebnH9GRX0TOggXUlJTQHWqNe/2KHuvfltYsLrPbYug+7VIWZEnLumnRY/JZlydun6IaB7zU2sdycfPqYfOjdX7r+rIH1zV0n2KXq83NpaamhmWLWgHw92bb21n3pefiZi+1dp+JrqulNX4+WP2uoqKCjmBr3D7Fbj9RTI0Dw9spulz9wNq4fYqLqbo6bl2xNY4TbfPttsR9a0UP1BYUUFBQwJLnmobNX95t/bukxcul06YNmz+0nXw+Hz6fD799TBNts34gm9riYj7oSBxTdF+nTJkCMKxvxbbTlQmOaWzfG3oclnd/FbdPsduvtd+Hi5tXDVtnbDslaseh7RTtWx/Yx2Fo3xvc5qRJ+HtXJ+17YOWHsrIyWgKtCftebEzJ+l60naLLNQ40JOx7i5u91NrHPfr6ZDkiuq5EOQKsdqq1+8Pi5oZh84fmiI5gc9K+B1bOr6mpYYV9TBPliGhMyzrXDdsWN3vZpSxI/UD24PGvtt8/Q/tWbDslaudE7RRdbkXPqoR9b3Gzl1r7c2lx81dx+5Qo/qExDW2naN9aah+HRDkiGlNd35qkfQ+sz6bS0lIaB1oT9r3YdY30+QTrPnNbAk0b/Pl00KbJy0SO9hKUUmNMtIJJQERyjDFvAXNH+fphRGRKzH+PAKLrfxI4TkRyRWQzYEvgLeBtYEu74okX60bNJ+1SiS8BR9mvPxF4ItE27aqKSik1nDGc+eydANy973GYkpKvvaqCbM01SqnU01zjXqM9A75CRL5pjPkYa6B8hoi0YVUcWS/7bPk+QKWI1AG/BvYRkR2wLkFZCZwGYIz5WEQeBD7BKnV4ljEmbK/nbOB5IBtYaMcD8AvgfhG5AngPuDNRHFoHXCmVzL7LlrDN6uU0F5Xz4LcOH9WjgJPR2rxKqXTQXONeo225X2Hd7AhwCXAPUAicOZoXG2O+n2BywkGyvfyVwJUJpj8DPJNg+n9YVyklKa0DrpRKJCsS5nS78smd+x3PgDdvo9Y3pzSU8FIZpZQaS5pr3GtUA/BojW779zeBLVIWUQplZ2evfyGlVMY5aOm/2LzxK1aXT+HxXQ9e/wvWo2FAa/MqpVJPc417jbYKyjbA3kA50Ar82xjzSSoDU0qpdCjp6eCcZ24H4PYDTiDk2fhL1YKRjV6FUkqtl+Ya9xpxAC7Wk2vuxLqxsQ5Yg1Vje6qI/B34sXHRnY1aB1wpNdTPH7+Riu52lm62Hc/uuP+YrHNafmSwaoBSSqWK5hr3Wt93Fz/Funlyd2NMjTFmD2PMpsAeWGfET0txfGNKb8JUSsXa78PFzPvgZXq9eVx2zM8wWWPzdW5sqTCllEoVzTXutb5Pmx8B5xpj3o6daP//fHu+a4RCofUvpJTKCOVdrSx49HoAbvzOaayumDpm696yUL9tU0qlnuYa91rfAHwbYHGSeYvt+Uop5S7GsOCxGyjt7eStLXbi0d0OGdPVe8Q1V+YppVxMc417rW8Anm2M6Uo0w57uqttvPR79qkYpBTmPPMK+y5bQnVvA5UdfPGaXnkR9pF8LK6XSQHONe62v5XJEZF9AvubrxxWtA66UquxopuCqnwNw/SGns7asesy3sUOJ1uZVSqWe5hr3Wt8AuhFYuJ75rqF1wJXKcMbwy0euI6u9nVe32pUnxqDmdyL1/a76clAp5VKaa9xrxAG4MWZGmuJQSqmUO+Td59n7szeJFBdz5VEXgiT7ck8ppZRKnYz600nrgCuVuarbG7noyb8A0HfNNTSVVKVsW1Py9OkYSqnU01zjXhk1ANc64EplKGOoffiPFPb3sHibbxE49tiUbu79DlfdHqOUcinNNe6VUQNwrQOuVGY64s2n2WP5O7QXFHH19y5I+aUn2xZrrlFKpZ7mGvfKqAG4UirzTGldy3lP3QrA7w8/l5ai8pRvM2T02nKlVOpprnGvjBqAax1wpTKLRCL870PX4gv08a9tv82i7fdNy3a/6NaKS0qp1NNc414ZNQDXOuBKZZaTX7qXXVa8T6uvlN8dcV7aqp7M1q+FlVJpoLnGvTJqAK51wJXKHLsuf4fTXriLiAi/Oe4XtBeWpm3bdX0ZlVqVUg7RXONe2nJKqQmnur2RK+67iixjuHO/43l9q13Tuv0czaxKqTTQXONeGdV0WgdcqQwQCHDVPy6jrKeD12ftzB37/yjtIVTnam1epVTqaa5xr4wagGsdcKUmvvxLL2W7rz6loaSK//3+L4lkpf/Ss3fb9YZvpVTqaa5xr4wagOtNmEpNbDmPPELe7bcTzPaw4PhL6fCVOBLHnFK9MUoplXqaa9wro/50kjRVQFAqU8y7470R579w6o5pigRmNPjx/fZ8AP50yBksq9kmbdseqjesuUYplXqaa9wrowbgWgVFqYkpf6CP3/39N0hPD4Ejj+Sh3Q5zNB5/r+YapVTqaa5xr4y6BEUfRa/UBGQMtQ//kc0bvyI8axY911+ftnrfyWxdpLlGKZV6mmvcS8+AK6Vc7ejXnuDAD16i15tH8O67obDQ6ZBYqWellFJpoLnGvTLqDLgxxukQlFJjaLb/Ey546hYArjjqYiJbbeVwRJYij+YapVTqaa5xr4wagEciWi9TqYmitLuda/5xGTnhEPfveQSLdtjX6ZAGVXg11yilUk9zjXtl1ABc64ArNXFc/MSfqe5o4sNNt+aG75zmdDhxtDavUiodNNe4V0YNwLUOuFITw3YrP+bAD16i3+PlVz/4FSHP+PrjWmvzKqXSQXONe2XUADwrK6N2V6kJSSIRLvznzQDc8+2jqS+f7HBEw3WFtDavUir1NNe4V0aNSHUArpT7zX//Rb656nOai8q5e9/vOx1OQvV9mmuUUqmnuca9MqrltA64Ui7X08NZz94BwM0HnUJfbr7DASU2qyjsdAhKqQyguca9MurqfY8no3ZXqQkn76abKOto5tNNtuTpneY5HU5SK3q0Nq9SKvU017hXRp0B1zKESrmXrF5N3k03AfCnQ8/EjONLyrQ0mFIqHTTXuNf4/QRLAR2AK+Ve+ZdfjvT18a9tv817m2/ndDgjKs3Rh2MopVJPc417ZdQAXOuAK+VO23z1GbkPPojxernpOz91Opz10tq8Sql00FzjXmkZgIvIQhFpFJFlMdPKRWSRiHxh/1tmTxcRuVFEvhSRD0Vkp5jXnGgv/4WInBgzfY6IfGS/5kYRSViXR+uAK+VCxnCRXXZw4IwzWFM+xeGA1k9r8yql0kFzjXul6wz4XcD8IdMWAC8aY7YEXrT/D3AQsKX981PgFrAG7MCvgd2AXYFfRwft9jI/iXnd0G0BWoZQKTc64IOX2M7/CZGqKvouuMDpcEalPai1eZVSqae5xr3SMiI1xrwCtA6ZfBhwt/373cDhMdP/ZixvAKUiMgU4EFhkjGk1xrQBi4D59rxiY8wbxhgD/C1mXXF0AK6Uu+QGBzjnmb8C0FdbC8XFDkc0Oi0BzTVKqdTTXONeTl48VG2Mqbd/XwtU279vAqyKWa7OnjbS9LoE04fROuBKjd68O95LOu+FU3dMSww/eOVhprQ3snzKTKp++MO0bHMszPSFqevT8mBKqdTSXONe4+LqfWOMEZGU38rb3t7OnnvuicfjIRwOc+SRR3LWWWexdu1afD4f2dnZdHZ2UlVVRWtrK8YYqqqqaGhooLCwEIDu7m6qq6tpampCRCgvL6epqYni4mLC4TA9PT1MnjyZtWvXkpOTQ0lJCc3NzZSUlBAIBOjr6xuc7/V6KSoqoqWlhbKyMvr6+ujv7x+cn5eXR35+Pm1tbVRUVNDV1UUgEBicn5+fj9frpaOjg8rKSjo6OggGg4PzdZ90nzZmn6blh5npC7O8K5sp+RGKPIZ32z3MKQ3R2NiI1+tlbmWAT7s81BSEKcheN79hIItgBPx+P1VVVexeHsQjho86PexQEqK+3zprMyUvQn9/P01NTexeHuSL7mxmF4eo68uisquFU16+F4AXj/8J325ooKSkhLmVAVb2ZlPkMVR4I4Pb7AoJ3d3dtLS0MC0/TIU3QmnOupjag0JLIIuZvjDd3d3sVBqM26eWQBZdIWFGQZhPuzysWbOGYDBIoScSt0/T8iMs6/SwZWGYVatWUVVVxdzKQNw+rekT9qoIEDLCF93Z+P1+iouL+UZRiOrcdTH3hgV/bzZbF4Vob28nEAgwtzIQt0/1fVnMKgqzoiebCm8Ev9/P5MmTmVsZiNun2HYaGBhg7dq1zC4Oxe3T0HZqbm4mOzubuZWBwX0a2k4tLS10d3dTkhNh2+LQ4D5F2ykni7iYYvdpaDv5/X68Xi9TcsNx+xTbTg0NDeTn5zO3MpCw70Xbye/3U1lZyS5lwYR9b1p+hN7eXlpbW9mrIpCw773f4WHb4hB1dXWUl5cztzIwuE9D26mnp4fm5mZqCsIJ+160nTo7O9m+JJi07y3vyqa+vp5AIDDYt4b2vZqCcNwxTdT3ou3k9/spLCxkVmEobp9i26mtrY1wOBzXtxK1U+w2E/W90px1fWv7kmDCvhebI2oKwkn7XsNAFm1tbXR2dlLhjSTse2D1rerq6rj3U2zfi7aT3+8nJyeHSbmRhH1vTmmI+vp6ioqKmFsZSNj3ou3k9/upqKgYMUf09PTQ0dERd0yT5YiSnEjSvhcyQm9vL01NTcz0hRL2PX+vlUNKSkqYXRyK26csTNz7qaGhgf7+/sG+lShHxLZzor4XbSe/34/P52OmL5Sw7+1QEqKlpQWAuZWBhH0vNkdEY0qWI6J9a6fSYMK+F41pfZ9PLYEs2tvb6ejoYFJuZKM+n97v8LBq1SpEhApvJGHfm1MaYs2aNQk/nyB58Q8nB+ANIjLFGFNvX0bSaE9fDUyPWW6aPW01sM+Q6S/b06clWH6Y8vJyXn311WHTa2pqBn8vK7MuKy8oKEg4v6KiAoDp06cnnF9ZWTlsms/nG3Gb0QFWccxWzJM7AAAgAElEQVTX67Hzo9Ojyw2dX1paOmw7uk+6Txu7T3V9awbPrNQPrDvDsrjZS+2kSYO/AzQOZMXNH7quN1pzEs5f3g21eXlMnz6dN55vjpv/g4f/Ru5APy99c08WVu3CcVOnDnv90HUWFhZSWFhI3aLWuLNCsa+p68umtrCQpe2JY/L3Wq+bam+vO9QaN39Fj/VvS2sWl9ltMXSfdioNsqRl+HH4rMvDZ13Dt9k44KXW7h+Lm1cPmx89/nV92YPrGrpPscvV5uZSU1PDskWtcfs0tJ1q7T4TXVdL6/B2rKiooKKigo5ga9w+JWrnofs0VHS5+oG1cfsUF1N1ddy6hva9oet6uy1xO67ogdqCAgoKCljyXNOw+cu7rX+XtHi5dNq0YfOHtpPP58Pn8+FfFH815dB2qi0u5oOOxDFF93XKFOtG4qF9K7adrkxwTGP73tDjsLz7q7h9it1+rZ1bFjevGrbO2HZK1I5D2ynatz6wj8NIOcLfuzpp3wMr55WVldESaE3Y92JjStb3ou0UXa5xoCFh31vc7KXWPu7R1yfLEdF1jZQjau3+sLi5Ydj8oTmiI9ictO+BlfNrampYYR/TRDkiGtOyznXDtsXNXmuwOpA9ePyr7ffP0L4V206J2jlRO0WXW9GzKmHfW9zspdb+XFrc/FXcPiWKf2hMQ9sp2reW2schUY4YzecTWJ+3paWlNA60btTnE6z7zG0JNCXNEbVJPp8O2jT5uWUnB+BPAicC19j/PhEz/WwRuR/rhssOe5D+PHBVzI2X84BLjDGtItIpIrsDbwInADcl2qDWAVdq5EtLIH2Xl4xkq7rlHPLu8wSzPdz4ndOcDmeDFXm0Nq9SKvU017hXWgbgInIf1tnrShGpw6pmcg3woIicAviBY+zFnwEOBr4EeoGTAeyB9uXA2/ZylxljoqcizsSqtJIPPGv/DKN1wJVyAWO48J+3kGUM9+55BHWVCW/pGNe0Nq9SKh0017hXWlrOGPP9JLP2S7CsAc5Ksp6FwMIE098BZq8vDq0DrtT4d9hbz7DTfz+kzVfCnfsd73Q4X8uc0lDCS2WUUmosaa5xr4yqX6NlCJUa36o6mjj/6dsA+MN3z6Y7v3A9rxiftDSYUiodNNe4V0a1XJIHZCqlxgNjWPDoDRT29/DK1nvwwg77Oh3R19YV0lyjlEo9zTXulVED8HA47HQISqkk5r3/Et/+9HW683xcc+R54OI/mGcUaK5RSqWe5hr3yqgBuMejNysoNR5JczMXP/lnAG74zmk0lVQ5HNHG+bRLc41SKvU017hXRg3A9Qy4UuNTwSWXUNbTwdtb7Mjjux7sdDgbrUbPSiml0kBzjXtl1ADcKrCilBpPcp57Du8jj9CXk8eV37vQ1ZeeRBVka65RSqWe5hr3yqgBuNYBV2p8KezrpuCiiwC4Zf6PWV0x1eGIxobW5lVKpYPmGvfKqAG41gFXanw57+nbyKqvJ7TLLjyw5+FOhzNm5pSGnA5BKZUBNNe4V0YNwLOzs50OQSll2+WLpRz+1jMYr5eeG28kkjVx3p8NAxmVWpVSDtFc417ackqptMsL9FH7yHUA9P/sZ0S22srhiMZWMOJ0BEqpTKC5xr0y6uIhrYKiJrp5d7w34vwXTt0xTZGM7MznFrJJaz2fT92CSeee63Q4Y25afoQVPU5HoZSa6DTXuFdGnQHXmzCVct52Kz/m2FcfI5SVxWVHXwwT8H25rDOjzm0opRyiuca9MmoAHgrpzQpKOckbDPCrh/9AljH8fe6xLN9kS6dDSoktC/XbNqVU6mmuca+MGoArpZx10kv3slnjV6ysms4d+5/gdDgp4xGtzauUSj3NNe6VUQNwfRS9Us6paVzFiS/dD8CVR11IIMfrcESp85F+LayUSgPNNe6VUQNwrQOulEOM4RePXY83HOTxXQ7i/c22czqilNqhRC93U0qlnuYa98qoAbjWAVfKGQe99y92WfE+7QXF/PngnzgdTsrV92dUalVKOURzjXtpyymlUkra2zn/n7cCcMMhp9HhK3E4IqWUUspZGTUA1zrgSqVf/mWXUd7TztLNtuOpOQc6HU5aTMnTp2MopVJPc417ZdQAXOuAK5Ve2W+9Re5ddxHKyuaaI88DEadDSov3O/TGKKVU6mmuca+MGoBrHXCl0ic7HKbgoosA+PvcY/hv9QxnA0qjbYs11yilUk9zjXtl1ABcKZU+x776KJ6PPya86abcud/xToeTViGTGWf6lVLO0lzjXhk1ANc64EqlR3V7I6e9cBcAvddey4A3z9mA0uyLbq24pJRKPc017pVRA3CtA65Uelz05M0UBPoJHHoooQMOcDqctJutXwsrpdJAc417ZdQpYa0DrlTq7f3Ja+y7bAk93nyCV13ldDiOqOvLqHMbSimHaK5xL205pdSYyQv08bPHbwLg1vknYzbZxOGInJGjmVUplQaaa9wro5pO64ArlVqn/uvvTGlv5LOpW/DQHoc7HY5jqnO1Nq9SKvU017hXRl2ConXAlZvNu+O9pPNeOHXHNEaS2Mz6//DDVx4mIsLVR15AOIMv+Xq3PaNSq1LKIZpr3CujzoDrTZhKpUZeoI9LH/oDnkiYR3Y/lE82/YbTITlqTqneGKWUSj3NNe6VUQNwyZCn8CmVTp5QkN//7bdsU/c59aWT+Mv8U5wOyXG9Yc01SqnU01zjXhk1ANcqKEqNraxImN8+cA17LH+bNl8J55z6e7rzC50Oy3H+Xs01SqnU01zjXhk1ANdH0Ss1hozh54/dyLwPXqY7t4BzTrkG/6TpTkc1LmxdpLlGKZV6mmvcK6Ou3tcz4CrdxvuNkxsj78or+d6bTzHgyeGik67g82mznA5p3FipZ6WUUmmguca9MuoMuDHG6RCUmhBy//xn8q+7jlBWFpccfylLZ27vdEjjSpFHc41SKvU017hXRp0Bj0S0XqZSG+vQt5+j4KFrAbjs6J/z722+5XBE40+FV3ONUir1NNe4V0YNwLUOuFIbZ59lS6h9+I8A9F59Nc9m7+pwROOT1uZVSqWD5hr3crzlRGQl0AWEgZAxZmcRKQceAGYAK4FjjDFtYtURvAE4GOgFTjLGLLXXcyLwK3u1Vxhj7h66La0Drsajka4Th/FzrfguXyzlynuuINtE6PvZzxg47TRYT+yZak5piMXNXqfDUEpNcJpr3Gu8XAO+rzFmB2PMzvb/FwAvGmO2BF60/w9wELCl/fNT4BYAe8D+a2A3YFfg1yJSNnQjWVnjZXeVcpdtvvqMP9z9v3jDQR741uH0L1iw/hdlsK6Q1uZVSqWe5hr3Gq8j0sOA6Bnsu4HDY6b/zVjeAEpFZApwILDIGNNqjGkDFgHzh65UB+BKbbhZq7/ghoULKAj08+yO+/HH754F+lCrEdX3aa5RSqWe5hr3Gg8tZ4AXRORdEfmpPa3aGFNv/74WqLZ/3wRYFfPaOntasulxtA64UhtmzpfvcdutF1La28W/t96d3x7zc4z+Ibtes4rCToeglMoAmmvcy/FrwIG9jDGrRWQSsEhEPoudaYwxIjImdXba29vZc8898Xg8hMNhjjzySM466yzWrl2Lz+cjOzubzs5OqqqqaG1txRhDVVUVDQ0NFBZaT/fr7u6murqapqYmRITy8nKampooLi4mHA7T09PD5MmTWbt2LTk5OZSUlNDc3ExJSQmBQIC+vr7B+V6vl6KiIlpaWigrK6Ovr4/+/v7B+Xl5eeTn59PW1kZFRQVdXV0EAoHB+fn5+Xi9Xjo6OqisrKSjo4NgMDg4X/fJ+X2aWxmgJZBFV0iYURDm0y4PNQVhCrINAwMDrF27lm8UhQhGYFp+hGWdHrYsDOMRw0edHnYoCdHS0gLA3MoA73d42LY4RMgIX3RnM7s4RF1fFs3NzfT09FDoiTCnNERvWPD3ZrN1UYiVvdkUeQwV3ggDAwPMrQzQFRLq+7KYVRRmRU82Fd4IpTmGd9s9+P1+8vLyOOrzl7no7qvxhEK8vfNcnvrJheT3ZDGnNEBjYyNer5e5lYG4fXq33cOc0hANA1kEI+D3+6mqqmL38mDcPtX3W4P4KXkR+vv7aWpqYvfyYNw+5WRBdW5kcJ1r1qyhpKSEuZWBuH2Kzu8KCd3d3bS0tDAtPxy3T3NKQ7QHhZZAFjN9Ybq7u9mpNEiRZ938oe20Zs0agsHg4DGN7lNsO61atYqqqirmVgbi9qmhX9irIjDYTn6/n+LiYr5RFIrbp9h2am9vJxAIMLcyELdPQ9vJ7/cP9q3YfVrelc2U/AhFnnV9a3ZxKGHfi66/ubmZ7Oxs5lYGEva9+v4sWlpa6O7upiQnkrDv5WQRF1OyvhftW16vlym54YR9b05piIaGBvLz85lbGYjbp6Ht5Pf7qaysZJeyYMK+Ny0/Qm9vL62trexVEUjY96Lvp7q6OsrLy5lbGUjY93rDQk9PD83NzdQUhBP2vWg7dXZ2sn1JMGnfW96VTX19PYFAYLBvJcoRscc0Ud+LtpPf76ewsJBZhaG4fYptp7a2NsLhcFzfStROsdtMliOifWv7kmDCvhddf2NjIzUF4aR9r2Egi7a2Njo7O6nwRhL2PbD6VnV1ddz7KVGO8Pv95OTkMCk3krDvzSkNUV9fT1FREXMrAwn7XrSd/H4/FRUVI+aInp4eOjo64o5pshxRkhNJ2vdCRujt7aWpqYmZvlDCvufvtXJISUkJs4tDcfvkERP3fmpoaKC/v3+wbyXKEbHtnOzzKXpMfT4fM32hjfp8iuaI0Xw+rV27lp1Kg+v9fJqWH07a91oCWbS3t9PR0cGk3MhGfT693+Fh1apViAgV3sgGfz5B8uIfjg/AjTGr7X8bReQxrGu4G0RkijGm3r7EpNFefDUQ+6i9afa01cA+Q6a/PHRb5eXlvPrqq8NiqKmpGfy9rMy6dLygoCDh/IqKCgCmT5+ecH5lZeWwaT6fb8RtRgeNxcXFCedHp0eXGzq/tLR02HZ0n8bHPsXeHBN9ZHDjgPXmrs3Npaamhs8WtQKwosdarqU1K+71tfa+LG7+CoAlLd64+dH9qayspDvUGrfNxoH4m3Nqc3Pj5tcPWDHV9a17mENNTQ3eu+7iZwuvIMsYHvjW4fzxu2dhOrLWxTRpUtz2o/s0dJ+jx+eN1pyE85d3Q21eHtOnT+eN55uHzf+sK+Y4TJ06bP7QdRYWFlJYWEjdota4fYp9TV1fNrWFhSxtTxxTtJ2m2tsbekxj2+kyu38N3aftS4Jx7RQ9Dp91eeL2KapxwEut3ecXN68eNj+2naLrGrpPsctF+9Yyu28N7XvR19fa74Pouob2PbDeSxUVFXQEWxP2vdj9G6nvxS5XP7A2Yd9b3Oyltro6bl3R5ZJt8+22xO24ogdqCwooKChgyXNNw+Yv77b+XdLi5dJp04bNH9pOPp8Pn8+H3z6mibZZP5BNbXExH3Qkjim6r1OmTAGG963YdroywTFNlCOix2F591dx+xS7/Vo7Xy5uXjVsnbHtlKgdh7ZTtG99YB+HoX1vcJuTJuHvXZ2074GVx8vKymgJtCbse7ExJet70XaKLtc40JCw7y1u9lJrH/fo65PliOi6RsoRtXZ/WNzcMGz+0BzREWxO2vfA+hyrqalhhX1ME+WIaEzLOtcN2xY3e9m+JEj9QPbg8a+23z9D+1ZsOyVq50TtFF1uRc+qjfp8iq5rNJ9PNTU1LLWPQ7LPJ2vamqR9D6wxRGlpKY0DrRv1+QTrxhEtgaYN/nw6aNPk548d/S5ZRHwiUhT9HZgHLAOeBE60FzsReML+/UngBLHsDnTYl6o8D8wTkTL75st59rQ4WgdcqfUwhrzf/x7fhReSZQy3HHgyfzjsbL3sZAOV5ujDMZRSqae5xr2cPgNeDTxmVRfEA9xrjHlORN4GHhSRUwA/cIy9/DNYJQi/xCpDeDKAMaZVRC4H3raXu8wYE3+aAq0DrtRIsiJhLnryZvJfewKTlcXVh5/HY7sf4nRYrqS1eZVS6aC5xr0cbTljzH+AYc+wNsa0APslmG6As5KsayGwcKTtaR1wpRLLCQX47f2/44APX8bk5tJz++081jDN6bBcS2vzKqXSQXONe2XU98pahlCp4Qr6e7l+YS0HfPgy3Xk+uh96iOChhzodlqu1B7VMo1Iq9TTXuFdGjUh1AK5UvLLuNm657SJ2/XIpLYVlnHb6dYT22svpsFyvJaC5RimVeppr3CujWk7rgCu1ziYta7jj5vPYZvVyVlVM5ZSzbmT51C2cDmtCmOnT2rxKqdTTXONeGXX1vseTUburVFLbrPqMPy2spbynnc+mbsF5p1xNa1G502FNGMu7ste/kFJKbSTNNe6VUWfAtQyhUvCtz97k1lsvorynnddn7cxpp/9JB99jbEq+5hqlVOpprnGvjDolrANwlem899zDH+/6FZ5IhKfmHMCV37uIkEfLc461Io/W5lVKpZ7mGvfKqAG41gFXGcsYTnnxH/heuAuAhf/zA2458Mcgegd9KmhtXqVUOmiuca+MugRF64CrTJQdDnPJo3/i9Bfuwojwu8PP5Zb5p+jgO4XmlOoN30qp1NNc414Z9aeTliFUmSY30M9V91zBtz99nX6Pl9DCO3h47SZOhzXhaWkwpVQ6aK5xr4xqOdEzfiqDlPR0cMvtF/PtT1+nI7+Is356LcFD9NHy6dAV0lyjlEo9zTXulVED8HBY62WqzLD52pXcefM5bPvVp9SXTuLUM2/gwxmznQ4rY8wo0FyjlEo9zTXulVGXoGgdcJUJDn73BS559HryggMsnzKT8358Fc0llU6HlVE+7dJco5RKPc017pVRLadnwNWE1tdH7UN/4PC3nwXg6Z0O4Jojz6Pfm+9wYJmnpiBM40BGfcGolHKA5hr3yqgBuDFaL1NNTFkrVuA76SQO//hj+j1erj38HJ7c5SCtdOKQgmzNNUqp1NNc414ZNQDXOuBqItrvw8UUX3Yd0t2Nv3ITLjn+13wxdabTYWU0rc2rlEoHzTXulVEtp3XA1USSEwpw3lO3cexrjwMQOOwwTpxzCj15PocjU3NKQyxu9jodhlJqgtNc414ZdeFQdna20yEoNSamtK7lr7ecz7GvPU4w20Pv735Hz8KFOvgeJxr0mkylVBpornEvbTml3MQY5i/9F3+/4XS+uepz1pRVc+qZNzDwk5/o9d7jSDDidARKqUyguca9MuoSFK2Cotxsq7rlXPzkzeywchkA/956d35z7C/oLCh2ODI11LT8CCt6nI5CKTXRaa5xr4wagOtNmMqNpKmJ/Cuu4G//+AdZxtBSWMrNB53KU3MOxGTpl1jj0bLOjEqtSimHaK5xr4xquVAo5HQISo1adjjEMa8+TvEV95DV2UkoK5t79z6SO/Y7np78QqfDUyPYsjBMS6v+caSUSi3NNe6VUQNwpdxit+XvcNGTN7NZ41cABPfbj+/v+CP8k6Y7HJkaDY9obV6lVOpprnGvjBqA66Po1Xg3rXk15z91K3M/eQ2Aryo2ofzPfyA4bx7+O993ODo1Wh/p18JKqTTQXONeGdVyWgdcjVebNq3ipJfu46Cli/BEIvR487lz/+O5f68jefrA3ZwOT22gHUq0Nq9SKvU017hXRg3AtQ64Go15d7yXdN4Lp+44ptvK+uQTrrjnCvb/cDHZJkIoK4snd57PX+b/mJbiijHdlkqf+n69JlMplXqaa9wrowbgSo0XW6/6HN+PrsP79NMcCASzPTw55yDu3vc4VldMdTo8pZRSSqVQRg3AtQ64ctp2K5dxyov/4Fufvw2AycvjgTkH8Y+5x9BQOsnh6NRYmZIXYXm301EopSY6zTXulVEDcK0Drpyyrf9jznz2Tnb+zwcA9HrzyDrtJ/SfeSZ//Ocah6NTY+39joxKrUoph2iuca+MajmtA67GykjXicO6a8U3aVnDWc/ewQEfLgagK8/HA3sewf17HclD5+1jL60D8Ilm2+IQS1r0xiilVGpprnGvjBqAK5Uu0t5O3h/+wEO33U5OOER/Ti737H0Uf597jD5EJwOEjDgdglIqA2iuca+MGoBrHfDMNtqz1hvDEwpy9OtPUnzVvWS1txMR4ak587jlwB/TWFq10etX7vBFt1ZcUkqlnuYa98qoEanWAVcpYwz7Lvs35zzzV6a3WJeUBPfem5N2+iHLN9nS4eBUus0u1tq8SqnU01zjXhk1ANc64GqsSCRCeU87k9samNpazzGvPcEOK5cB8N9Jm1J1w+8IzpvHcn16ZUaq69PavEqp1NNc414ZNQBXE9eYX15iDFWdzcxoXMXk9gamtDUwub2Rwqf6yFq1iiVf1eENx3+j0uor5fZ5J/L4rt/h2QN33tBdUBNIjn4mKqXSQHONe2XUAFzrgKuhssNhprbWs1njV+Te8ArZy5fzf0veY0bTKgr7e5K/DmjzlbC2dBJrSyfx2bRZPLDnEfTk+dIXvBq3qnMjfNbldBRKqYlOc417ZdQAXOuAjy+jPWv9dR8NL5EIJb2dVHW24FnUzGFvvklVZwuVnS1M6mxmSlsDmzYNP5M92/63vaAIf9Wm1FVMtQbaZdWce9weRKZP55B/NdPvzR/djqqM8257RqVWpZRDNNe414RqORGZD9yAdYLyDmPMNbHz9SZMd/KEgvgGeikY6KNgoJfivm6Kezsp7usmd+A1pL2dn73+JcW9XRT3dVHU101FVytVnS3khO3a79fDr5Ksv6Gkiv9O2pQd9p1DeKutuHC58N/qTWnzlYLEl3g6cz9rwN//SvKz40rNKdUbo5RSqae5xr0mzABcRLKBm4EDgDrgbRF50hjzSXSZ9vZ2p8JzhdGcad7Qs9bVbQ1Ma1lDQaCP/IF+CgJ9XLBjOdLTwzmvraAg0E/+QB8FAWtwnT/Qh2+gj/xAHyVXB5Dubl4f6Q+nh6x/jkkyuyO/iObiCjbdZjOe6/LSVFxBc3EFjcWVNJZU4a+aTm9eQVzsS9ezj0qtzzsvPQvbHuZ0GEqpCU5zzfjW2tpamWzeRLp8f1fgS2PMf4wxAeB+IK5XjnYA/txzz6VtmfG8rvZPXtvoZQ57+1luvf1irrvrf7nyviupfeQ6Cn71K/KvvpoTFj/IUa8/yXeWLmLfZUvY7YulbPfVp8xsWMnUtgay2tqQYJBQVjYd+UWsLZ3EF6VVvDdjNq9svQdPzTmA/tNOo2/BAq497GwuPW4B5598JaeceSP7HXkBe135DPv/9nGOu+hOuh95hMuO+Tm3zD+Fh751OItn78XrXS2Dg+9UH4exXle6tzde1zVeY3/v5WfHbF1uPg7aH5xZl5tjH8t1uTn20S43mlwzXmN387pGu73Ozs6kDwCZSAPwTYBVMf+vs6cNMsaMakXjZQA+7473Bn/Ovf6euP8nWmak5b5OXO2fvLHRy6yq2IR3N9+eF0om8cL2+/D4LgfRf/rp9F10ETcddCq/P/wcfnPMz/n5j37N0Vvtyiln3shxF/yVwxb8g/YvvqCtvp49rn6e/X/7OIf+8j6+NWtXfnrmDVx08hX89tgF9F19Nf0//zkP7nkEz+50AK9uvTsfzvgmH9YtZyAnd6P3b6yOw1ivy82xj+W6xmvs+aPMrBP9OGh/cGZdbo59LNfl5thHu9xocs14jd3N6xrt9kYiox2UjncichQw3xhzqv3/HwG7GWPOji7zz3/+s7+xsXGwFEpxcXFTeXl589B1dXZ2lhQXF3eMtL2xWkbXNb63N17X5ebYx3Jd4zX21tbWykS5JVVxjdfjoP3BmXW5OfaxXJebYx/tcqPJNeM1djeva7TbGxgY2Orggw8uSjRvwlwDDqwGpsf8f5o9bdChhx6al9aIlFJKKaWUGmIiXYLyNrCliGwmIl7gOOBJh2NSSimllFIqzoQ5A26MCYnI2cDzWGUIFxpjPnY4LKWUUkoppeJMpDPgGGOeMcbMMsbMNMZc6XQ844WIzBeRz0XkSxFZYE/7HxFZKiLLRORuEZkwf4wlIiILRaRRRJYNmX6OiHwmIh+LyO+dii9dRGS6iLwkIp/Y+3yePf1yEflQRN4XkRdEZKrTsaaSiOSJyFsi8oF9HH5rT99MRN603ysP2N+mTWhJ8oOIyJUislxEPhWRc52OM5US5YdMe0/AiPmhXEQWicgX9r9lTseaSiMchwfs/vC+iKwUkfedjjXVEuWHmHk3iki3U7GlS5L8cLTdNyIisvPXWrExZsL8ACdiPciw2P6/OB2T0z9Y3wasADYHvMAHwDZYFWNm2ctcBpzidKwpPg7fBnYClsVM2xf4F5Br/3+S03Gm4ThMAXayfy8Cltv9oThmmXOBW52ONcXHQYBC+/cc4E1gd+BB4Dh7+q3AGUlev5nT+zBGxyFZfjgZ+BuQZS83od8bSfJDRr0n7P1Mlh9+Dyywpy8Afud0rE4chyHL/BG4NA2x7Ah4HToOCfODPW9n4O9At9PtlYbjkCg/bA1cDXwEHG5Py9qQ9U6IM+AicoSIfAB8D/g5cAWAsY9IhktUH/17QMAYs9xeZpE9bcIyxrwCtA6ZfAZwjTFmwF6mMe2BpZkxpt4Ys9T+vQv4FNjEGNMZs5gPmNDvHWOJnrnJsX8M8D/Aw/b0u4HDY18nIgeISA9wp4gkre/qIsmen3AGcJkxJgIT/72RKD9k2nsCkucHrD5xt73YsPfFRDPCcQCsb4iwnv92X6piEJG5ItKO9YDB7VK1nfVImB/sBx9eizXemvCG5gcRORDrWMwEioE/2ctFNmS9rh+Ai8g04LvAOcaY7wLXAz4R2cnZyMaNRPXRJwOemK9NjiK+gkymmAXsbV9ysFhEdnE6oHQSkRlYZ1fetP9/pYisAn4IXOpcZOkhItn2V8iNWH+ErgDajTEhe5G4ZwmISDFWwj0XaAf2FxG359Bkz0+YCRwrIu+IyLMisqUj0Tks094TsYbkh2pjTL09azsTEMMAAAuPSURBVC1Q7VBYaTc0T9r2BhqMMV+kaJt5QAXWAPd14FsiUp6Kba1HsvxwNvBkTJ/IGCJSiHWlxQJjzDHAf4FVInLKhq7LlR8eIlImIgeLSIExpg7rq6Al9uzVwAygwan4XMBgVYn5k4i8BXQB4ZFfMiF5gHKsSw9+Bjxon9mY8Owk8ghwfvRMnzGm1hgzHbgHK8FOaMaYsDFmB6ySpbsC30iwWJada/Lt4/SYMeZOrLMfPwBq0hdxWuUC/caYnYG/AgsdjscRmfaeiEqUH6Lsb5Yn/LcBMOJx+D5jfPbbvs7+BBEpM8b0Ay8YY24HnsG63GOHsdzeRigAjgZucjqQdLFPvmQD2N+cPgU8H3MCpg7rBM4Gcd0AXEROBJYBZwJ/EZH5xphlxpiIfTACQAjIzZTB1HokrI9ujHndGLO3MWZX4BWsa9wyTR3wqH05wltABKh0OKaUE5EcrA+Ve4wxjyZY5B4m+CVJsYwx7cBLwB5Aqay7IfknWNc+ngncIiIHGWMa7Nc8iPWH69Ei4oPBr6XdJtnzE+qAaN94DOe+Ah8vMuY9kSQ/NIjIFHv+FKxvjSa0ZHnSzg9HAg+M4baOAd6z1/sHETkzeomcMeZFrPfjXBHZzF4+XWO3RPlhBbAF8KWIrAQKROTLNMWTdiJyDtY3noM3ohtjPjfGRGIuOZmJdRnjBn0OuG4ADmwLHGOMOQRYDJwoInvC4PU3s4Ei+5olIyIF4NoPx7GQsD66iEwCEJFc4BdYN5xlmsexbsRERGZh3WSy3qcXupn9PrgT+NQYc13M9NhLDA4DPkt3bOkkIlUiUmr/ng8cgHWd50tYl2SBdWnbn2NyzQnRXGO7EZgPTLL/78aBeLLnJwy+N4C5ZOAf6Jn2noDk+QGrT5xo/34i8ES6Y0unEY4DwP7AZ/a372Nlc6xLGg7Husb+aBE5KGb+vViD3s3t/6cr1yTKD48bYyYbY2YYY2YAvcaYLVIchyNEZCZWfj8TOIjEpbtzgSpjzCL7/yX2a9ffNmN9t2iqf7A6xMH27+XAecDfY+b/EDjfPmj3k6SKQSb9AAdjfYCuAGrtaddiDTg+x/p6zfE4U3wM7gPqgSDW2YRTsAbc/8D6RmUp8D9Ox5mG47AX1tfHHwLv2z8HY53pWWZP/yfWjZmOx5vC47Ad1hmnD+39vtSevjnwFvAl1k0337WnD8s19vRf2XnmZaybFh3ft69xLBLlh1Lgaaw7/F8Htnc6zhQfg0T5IaPeE/ZxSJYfKoAXgS+wKkeVOx2rE8fBnncXcPoYb+9x4Psx/z8JWDJkmePtXLMMuC2Nx2JYfhgyf0JXQbE/E+4DerG+JY/mhyPs3wNAD9aVBM8BN4123WJvYNwTkWxjTFhEzgDmGWOOsKdvBVwC3GeMeV5ErsG6g/9LrIfx3Oxc1Eopt1lPrvkl8A9jn+0QkYeA/YDLjTF/cixopZTriEiWsS6fPQG4wBizoz09F2uw/ZAx5l572oNYVZouN8bc4FjQGUZExBhj7OvAlwInGWOWxEw/CesemaXAncaYW0a77nF5CYqIHCMiu4pIkf1/sT8QBeuGhJB9zRRYZ6k6Wbcv22AdjL118K2UGomInC4i37V/94B1c+YIuaYD+2tIsaoIfQTU6OBbKTUSETlXRH4qVgk/wLps1h6E/w1oEZHTo4tjffNUZr92E+AdrOcP6OB7jCVqG3t6dJCdY6ybcG/Dem4KZt3Z62LgGmDuhgy+gfF1BlxE5gKXA/1YpW8asb7W7RMRj7HLg4nIEVhf/e5jjOkSkfuAR4wxD9uVUXqd2gellDuIyObAv4GwMWZTe1o2Vl4cKdc8aIx5LHr2yqn4lVLjn33fyNVYf7wvMMZ8bE8XrAe3hGOWexSYY4ypE5E/A+8YY+5yJvKJbz1tI4nyu4gswXog0RZYz5x5LdqGG2pcnAEXqx5vPnACcKMxZh7W9XdFWBVNMMaERGSGiFwE/D+s6zRvE5F/YZUd/MheTgffSqmEYiqcgFWq9C/AGhGJ3mhlRpFrPrEX1MG3Uiohe1zjxXpy6CfGmEONMR/bl5dgLGE71/weq8747cBlIvIKsBvWNehqjI2ybSJ221wdU33GA1QBB2JdbvLvrzv4BofPgNs7cw2Qh3WX77vGfiqhiCwENsWq0LECqw7m/cCvjTE326/dAtjV/vpGKaUSisk1OcA/jTH/EpEdsO4X+V+sp70V28vOxyo7d6nmGqXUhhgyrlmIVbltM+AFrGpLpVjVlt60p/8N66bKP9rfwFUBexljHk6werURvmbb/MUYc71YpR9PASqMMdeMSTxODcDtU/w3Y10/8yzWXb+PY/0FeBbWgXkD+CbWJSm/BEqNMa3R1xsn/3pQSrlCglxzMtY3bM8DPzXGLLC/7p2HlXzPBoqNMS3R12uuUUqtz5Bc8zxwLNa3aNEHfb2BVVFld6DJGPMzESk11rMIBm/KdCT4CW5j28Zex5i2T6KahulShPVkpwPtayubgUOBQ4wx10cXEpE5WNdgzjDG/Mf+CzGiH4hKqVEammtasGpcLwCaRWQLYCesRyx/YowJYt0QpblGKbUhhuaaVmAXrAfaXGGMWQYgIu/8//bu5sWqOo7j+PsTCAWaEklJhNK2wDYVBdUqsCyKQCKhRwitf8Bq5aYSoejBhRTtoqQWka2iRWVq4KJAQTdRA0YthGnGHhR8+Lb4HeEyzEwyOufe8b5fMFyYe3+Xc+Y3fO/nnvN7ALYlWVdVEwO1xvC9eC61b+py98/QxoB3M0onaFe+AQ7QZvk+kOTGgZeuBiar6peu3Tk/ECVdrDlqzVHams9baLcfd9LW+X55oJ21RtJFm6XW7AN+pQ15GNzkbTVwsqomunbWmkU2in0zzCvg0LY43pBkTVX9keQwbYOMNUmuBjYDTwLvgbeCJS3YzFpzjLat8nO0MeEF0NUda42khZpZa44A64EbkpwDXqDlmnfBWtOzkeqbYa+Csp/2zeNZgKr6kTYeB9ofYT1t96n3u+f9J5W0EDNrzSHa1sKTVVXdKkxU1a7u0VojaSHmyjVXAU/R5rc9VFUfdM9ba/ozUn0z1Cvg3TeQL4AdSX6mbTN/hrb04M5LWd5Fki6Yo9acpm0jTFWdGubxSboyzJNrzgBvO857eEatb0ZiI54kDwKbgHuAXReuQknS5WStkdQHa83oGpW+GYkADpBkGd0mGMM+FklXLmuNpD5Ya0bXKPTNyARwSZIkaRwMexKmJEmSNFYM4JIkSVKPDOCSJElSjwzgkiRJUo8M4JIkSVKPDOCSJElSjwzgkjRGkkwkOZXkryRTSQ4m2Zrkfz8PkqxLUkmGuouyJC11BnBJGj+PVNUKYC2wA9gGfDjcQ5Kk8WEAl6QxVVXTVbUXeAJ4JsltSTYm+SnJySTHk2wfaLKve5xK8neSuwGSPJ/kWJI/k3yVZG3PpyJJS4oBXJLGXFUdAn4D7gX+AZ4GVgEbgReTPNa99L7ucVVVLa+qH5I8CrwKPA6sBr4HPunz+CVpqTGAS5IAfgeuq6pvq+pIVZ2vqsO0MH3/PO22Am9U1bGqOgu8DtzuVXBJmpsBXJIEcBMwmeSuJN8kOZFkmhawr5+n3VrgnW5C5xQwCaR7P0nSLAzgkjTmktxBC8z7gY+BvcDNVbUS2E0L1AA1S/PjwJaqWjXwc01VHezj2CVpKTKAS9KYSnJtkoeBPcBHVXUEWAFMVtXpJHcCmweanADOA7cM/G438EqSW7v3XJlkUz9nIElLk2u5StL4+TLJWVqYPgq8RQvSAC8BbybZBXwHfEqbkElV/ZvkNeBAkmXAhqr6PMlyYE837nsa+Br4rNczkqQlJFWz3VGUJEmStBgcgiJJkiT1yAAuSZIk9cgALkmSJPXIAC5JkiT1yAAuSZIk9cgALkmSJPXIAC5JkiT1yAAuSZIk9cgALkmSJPXoP863PD3F5pJwAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 864x504 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(figsize=(12, 7))\n",
"\n",
"t = daily_new_cases(None)\n",
"plot_state(t[1], ax)\n",
"ax.set_ylabel('Daily New Cases');\n",
"fig.suptitle('United States', y=.95, fontsize='14');"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Daily New Cases by State\n",
"\n",
"States are \"sorted\" by how well they're containing the virus. Those states whose most recent 20 measurements of their 10 day moving average of daily new cases are in steepest decline are shown first. States whose 10 day moving average of daily new cases is trending up are shown last."
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment