Instantly share code, notes, and snippets.
Created
September 23, 2019 10:18
-
Star
0
(0)
You must be signed in to star a gist -
Fork
0
(0)
You must be signed in to fork a gist
-
Save mra1385/e370eb927ca62b6062cdb7eeb3809ed6 to your computer and use it in GitHub Desktop.
BikeshareSystemAge
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 1, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# import modules\n", | |
"import matplotlib.pyplot as plt\n", | |
"import matplotlib.ticker as tkr\n", | |
"import matplotlib.dates as mdates\n", | |
"import matplotlib.cbook as cbook\n", | |
"import numpy as np\n", | |
"import pandas as pd\n", | |
"import seaborn as sns\n", | |
"from datetime import datetime, timedelta\n", | |
"import os as os\n", | |
"%matplotlib inline\n", | |
"\n", | |
"# use fivethirtyeight plot style\n", | |
"plt.style.use('fivethirtyeight')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 2, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# load pickle\n", | |
"bike_share = pd.read_pickle('/Users/mazar/Desktop/personal/chart-it/posts/Bikeshare/pt 4/bs_pickle_2010-2018')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 3, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# summary bike stats\n", | |
"bikes = bike_share.reset_index().groupby('Bike').agg({'Minutes': ['median', 'sum', 'count'], 'Start Date': 'min', 'End Date': 'max' })\n", | |
"stations = bike_share.reset_index().groupby('Starting Station').agg({'Minutes': ['median', 'sum', 'count'], 'Start Date': 'min', 'End Date': 'max' })" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 4, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_bikes = bike_share.reset_index().set_index('Bike')\n", | |
"bike_share_stations = bike_share.reset_index().set_index('Starting Station')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 5, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# add start of operations for bike/stations\n", | |
"bike_share_bikes['Bike COD'] = bikes['Start Date']['min']\n", | |
"bike_share_stations['Station COD'] = stations['Start Date']['min']" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 6, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# add latest usage date for bike/station\n", | |
"bike_share_bikes['Bike Latest Usage'] = bikes['End Date']['max']\n", | |
"bike_share_stations['Station Latest Usage'] = stations['End Date']['max']" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 7, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# filter out bikes and station that have not operated during last 6 months (assumed to be decomissioned)\n", | |
"bike_share_bikes = bike_share_bikes[bike_share_bikes['Bike Latest Usage'] > '2019-02-28']\n", | |
"bike_share_stations = bike_share_stations[bike_share_stations['Station Latest Usage'] > '2019-02-28']" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 8, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# add bike/station life \n", | |
"bike_share_bikes['Bike Life'] = (bike_share_bikes['End Date'] - bike_share_bikes['Bike COD']) / np.timedelta64(1, 'Y')\n", | |
"bike_share_stations['Station Life'] = (bike_share_stations['End Date'] - bike_share_stations['Station COD']) / np.timedelta64(1, 'Y')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 9, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_bikes = bike_share_bikes.reset_index().set_index('Start Date')\n", | |
"bike_share_stations = bike_share_stations.reset_index().set_index('Start Date')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 10, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"name": "stderr", | |
"output_type": "stream", | |
"text": [ | |
"C:\\Users\\mazar\\Anaconda3\\lib\\site-packages\\ipykernel_launcher.py:2: FutureWarning: how in .resample() is deprecated\n", | |
"the new syntax is .resample(...)..apply(<func>)\n", | |
" \n" | |
] | |
}, | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAq4AAAGICAYAAACA4mKwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xd4m9XZBvD7aFqSbXnImU7IciYhYSeQQBIyGWWUj1FGoIUy20JLP6BQWsoqFEr7sSllUygjEEYCScggISEBEgLZJnvHe2nrPd8fTmK9r+RYtqVX6/5dVy7i15J1ZBvn9qPnPEdIKUFERERElOoMyV4AEREREVEsGFyJiIiIKC0wuBIRERFRWmBwJSIiIqK0wOBKRERERGnBlOwFRFNXV8dRB0RERERpzul0inh+PFZciYiIiCgt6BZchRB/FkJIzZ99ej0+EREREaU3vSuuGwF0D/szXOfH11V5eXmyl5BUfP58/tmMz5/PP1tl83MH+PwT/fz17nENSilZZSUiIiKidtO74tpPCLFbCLFVCPGWEKKfzo9PRERERGlKSKnPBn4hxDQAeQA2AOgC4G4AgwEMk1JWhd82fKpAtpfciYiIiNJJWVnZ4b/He6qAbq0CUsrZ4W8LIb4CsAXAdAB/b+1+4U8+3ZSXl6f1+juLz5/Pn8+fzz9bZfPzz+bnDvD5J7rgmLRxWFLKRgBrAWTvV5eIiIiIYpa04CqEyEFzq8DeZK2BiIiIiNKHnnNcHxVCnC6E6CuEOBnAuwAcAF7Raw1ERERElL70HIdVCuBNAC4AFQC+AjBKSrldxzUQERERUZrSc3PWJXo9FhERERFlHr0PICAiIiJKqlqfggV7vLAaBaaU5sBoiOvEJkogBlciIiLKCmuqA3hhfSPe3uKBO9g8Mv7aIQ78bVRBkldGsWJwJSIioowVUCQ+3u7B8+ubsGy/P+L9L29swj3H5yPPnLRBS9QODK5ERESkq/3uEJ5d1whFAr8cmoueDmNCHqfGp+DCOZX4tjLQ6m0CCrBojw9nH2VLyBoovhhciYiISFe3LK3F7J1eAMCKCj9mTXNBiPj2mbqDCi6eW3XE0HrIvF1eBtc0wbo4ERER6aYhoBwOrQCwbL8f5XXBuD5GQJG4akE1VlREtgaM7WbBrcNzVdfm7fZBShnXNVBiMLgSERGRbjbURIbUBXt8cfv4ipT41ZIazNml/piju1qw7Lwu+GhaCW4fmQ+bsaXCu6sphA218Q3PlBgMrkRERKSbdTWRL90vjGNw/dM39Xhrs0d1bXiRGW9NLMaQQjMAIMckMLa7RXWbebu8oNTH4EpERES6WRsluC7Z50NA6fxL9U/80IAn1jSqrvXJM+LdScVwWtSRZ2LPHNXbc3fHLzxT4jC4EhERkW6iVVwbAhLfRulHbY+XNjThj9/Uq66V5Bjw/mQXutojpxZMKlUH12X7fWgIKJ1aAyUegysRERHpQkqJdVF6XIHO9bk+s7YRty6rVV3LMwu8O7kYffOjD1Dqm29C//yWQBtQgC/i2LJAicHgSkRERLrY71FQ7Yte1exon+s/f2jAnSvqVNcsBuD1CcUYUWxp5V7NtO0C83azzzXVMbgSERGRLqK1CRzyTYUfdf7YX6qXUuLh7+rxJ017gNUIvDahGKf3sLb5MbTtAnN3cSxWqmNwJSIiorjxhlp/35GCa0gCS/bGVnWVUuL+lfV4aFWD6rrNKPDWGcWY0iunlXuqndrNipyw9tddTSFsjPNMWYovBlciIiLqtGpvCBM+OoBxy2z45aLqqJVLbX9rrkl9WtbCGIJrSJG4Y3kdHvtePT3AYRJ4Z3IxxveMLbQCgM0kMLabujI7l2OxUhqDKxEREXXak2sbsbIygBAE3t7iwYfbIwOgtuJ6+UC76u22+lw9QYmrFlbjufVNquv5ZoEZk4sxplvb7QFaEzXtAvN2cYNWKmNwJSIiok5bfkA9zurTnergGlIkNtSqg+svh+Qi7AArlNcFsasx+kv1Vd4QzvusEh9pAnGBRWDmVBdO7tr+0ApEH4vVyLFYKYvBlYiIiDpFSon1mjaAhXu8qnaBrQ1BVf+rK8eAfvkmnFii3vkfbSzWtoYgJn9SERGOu9sN+HhaCY51HXl6wJH0yzehX15Lo6tfAb6IsdeW9MfgSkRERJ1S4Y0cc7XXrag2Oq3VBNuhB49fHafZ/a9tF1hZ4cekjyuwuV6962tIgQlzzyrB0UXmTq+f7QLpg8GViIiIOkVbbT1kQdgxqtr+1qGFzQcDjI8SXJWDldpPtntw9qeVqPCqQ/HYbhbMPrMEpbnRDxdor4ixWLu9HIuVohhciYiIqFPW10Yfc7VwT0s/amRwba6UHl9iQb65pdG1yqfgh+oAnlnbiMvnV8MdVAfI/+lnw7uTXSiwxi/CjNGMxdrZGEI5x2KlJAZXIiIi6pQNrcxnXbLPD3+oOXhqg+uwg8HVZBAY011ddb12UQ3uXFEHbc3z1uG5eO60QliNAvFkMwmM0mzuWl/L4JqKGFyJiIioU1oLeU1Bia8r/HAHFWwJ61EVAAYXtLzMr20X2KSpdhoF8PjoAvzpBCcMIr6h9ZAB+eq2g52tTDeg5IpPcwgRERFlJSllq60CQPOUALtJqKqnffKMcJhbamfa4Bou1yTw8viiiA1U8dYr16h6e1fTEY4AS0EhReLdrR7sd4cwoWdOXDatpSIGVyIiIuqwvW4F9f7WNzIt2O1Fnzx1KDzU33pI/3wTSh3GiLDYw27Afye5MFyHEFbq0ATXxvQJrv6QxPQF1Zh9cHbuPd/U45SuFlw3NBdn9s6B2ZCYKnUysFWAiIiIOkxbbe1hVU8AWFUVwNJ96vmr2uAqhMC03uqK6vAiM+ad3UWX0ApECa5pUnENKBJXL2wJrYcs3e/H9AXVGPHOPjy6ugH1/sw4VIHBlYiIiDpsvWbT1QkFCoaG9a8qEnhvq1t1m2GFkWH0D8fm46QSC3KMwM8G2DHrTBd6aMJkImlHa+1Mg4prQJH4+cJqfLIj8njdQ/a4FTz+fUPERrd0xVYBIiIi6jDtxqz+dgWlLivWhV33aTLgoRmu4QqtBsw5uwRBRcKUhJe2u9kMMArg4BAEVPkUuIMK7KbUrPEFFIlfLKyOOAI33yxQH1DH1EsG2OG0pObzaK/MeBZERESUFNpRWP3sCsb3aH0jldXYfMxqa5IRWgHAaBARFd7dKdouEFQkrl1Ugw81obV/vhHLL+iKped1wdWD7LCbmj+X1w5xJGOZCcHgSkRERB2iSIkNmoprP7vEKV0taK3AN8hpTlo4bUu6bNC6ZWktPtjmUV3rm2fER1NL0N1uxNBCMx4/pRDrLuqGl8YVYnBB5kwYYHAlIiKiDtnZGEJT2MlWTotAiUXCYTbgpC6WqPeJ1iaQKnppguvOFKy4rqr04/Vydc9wnzwjPpoa2RNcYDXg/L52PZeXcAyuRERE1CHaauuQAjMOnQ/QWruAdqJAKilNg1muT6xpVL3dO7c5tGo3l2UqBlciIiLqEO1EgSGFrZ+GdUhKB1eHOvylWqvA9oZgRIvAY6ML0CtLQivA4EpEREQdpJ3hGt5LOaLYjEJrZC9ragfX1K64PrOuEUrYwIChBSZM7Nn6qWOZiMGViIiIOmR9jaZVICyUGg0Cp3dXtwsUWAS621M3ekS0CjQGW7ml/mp9Cl7bpO5tvenoXAiRmhvdEiV1v3uIiIgoZYUUiU11mlaBAvVL1tp2gaGF5pQOWtqK6253CIpMjdH9L25sUm2E62Yz4MJ+mbXxKhYMrkRERNRu2xpC8Ia9kl5sNaDEpg5+5/axqSqs0wel9jzRfIsB+ZaWYO0LAZXe5B+V6gtJPLdOvSnr+qG5sBpT95eARMmebl4iIiKKG21/65AoY64KrAZ8NNWFGVs9OLrIjKm9Wj+YIFWUOoxY529pEdjVGEIXm35Hz0bz9mY39ntaAnSuSeCqFP8lIFFYcSUiIqJ2izYKK5oBTjP+d2Q+zuxtgyGF2wQOSbVZroqUeFIzAuuKgXYUWLMzwmXnsyYiIqJOiRyFlbrTAtpDOw812ZMF5u3yYWNdyy8JRgHcMCw3iStKLgZXIiIiarfIUViZ0X0YeexrcicLPLGmQfX2+X1t6J1Fc1u1GFyJiIioXQKKRHld66Ow0lkqzXJdWeHH4n1+1bWbs7jaCjC4EhERUTttqQ8iELbZvpvNgMIM6blMpWNf71tZr3p7bDcLRrosSVpNasiM7zIiIiLSjfbggcEZUm0ForUKJCe4LtjtxYI9PtW1W47JS8paUgmDKxEREbVLxCisDOlvBYDudiMMYcMPKrwKPEF9DyFQpMSfv42stk7okV3Hu0bD4EpERETtkqkTBQDAZBDoYVdXXffo3C7w/lYPVlepP8f3nuBM6VPH9MLgSkRERDHbUh/El5oNQ63NcE1X2naBnU36TRbwh2REb+v5fWw4riS7e1sPYXAlIiKimJTXBXDW7ApU+Vp2ZlkMwKAMahUAIjdo7dSxz/XljU3Y1tDyeCYB3H1cvm6Pn+oy6zuNiIiIEmJDbQDnflqpOnoUAK4bmot8S2bVwZI1Eqver+CR1eq5rVcNcqC/k3HtEH4miIiI6IjW1QTwk08rUelVh9ZL+tvw5+MzrxqYrOD65NpG1efYYRL4/QhOEgiXWb8iERERUVz9UB3A2bMjQ+sVZXY8NaYQRkPmbRiKmOWqQ6vAfncIT61pVF276ehcdNVsFMt2SQuuQog/CCGkEOLJZK2BiIiIWrd8vw/nzK5AtU8dWn8+yIF/nlqQkaEVAHo51C9I79Jhc9Y939ShKWzslivHkPWnZEWTlOAqhBgF4FoA3yfj8YmIiOjI5u7y4rzPqlDrV88w/eUQBx4b7YQhg0czRTs9S8rEzXJ9o7wJ/93sUV37/Yi8jOsdjgfdPyNCCCeANwD8AkCN3o9PRERER/b2ZjcunVcFT0gd1m4alouHT878eaJOiwH55pbn6AsholUiXtbXBHDbsjrVtUFOE64e5EjI46U7kcjfIKI+oBD/BbBNSnm7EGIhgDVSypvDb1NXV3d4UeXl5bquj4iIKJu9tceEx7ZEzgz9ZW8/rukVRIZn1sMuWZmDze6W+t4rI7wYmhff8OoJAdNX52Br2ONYDRIvj/BigEPffBZPZWVlh//udDrj+h2j61QBIcS1AAYAuCLW+4Q/+XRTXl6e1uvvLD5/Pn8+fz7/bJWOz19KiQdXNeCxLepxTALAI6OcuHZIbP2W6fjco+m/rRKb3b7Db4uiHijrY2vzfu15/jcursFWt1t17bFTCjGtLH2rrYkuOOoWXIUQgwA8CGCslNLf1u2JiIhIPw+sasCjmhmiJgE8d1ohftrPnqRVJU+pwwSgJbjGeyTWG+VN+M+P6tB6SX8bLhuQfZ/r9tCz4joagAvAmrDeGCOA04QQ1wNwSCl9rd2ZiIiIEuPJNZGh1W4SeG1CEc7omZOkVSVX5Aat+E0W+K7SH7Wv9bHRBRnfP9xZegbXDwB8o7n2EoByNFdiWYUlIiLS2WubmnD31/Wqa06LwLuTXDixS2Sva7aIOISgk7NcmwIKPtzuxRvlTViyTx15bEaBl8YXwWHmFIG26BZcpZS1AGrDrwkhmgBUSynX6LUOIiIiajZzmwe/War6pxkOk8B7k104oSR7QysQv9OztjUE8fj3DZix1YOGQPQNV38b7cTQQnOHPn624ZGvREREWejz3V5cs6gaSliWshiA/5xRlPWhFYg+y7W9NtQGcOasyogDHMJNH2hnX2s7JDW4SinHJfPxiYiIstGKAz5cMb8agbA8ZRTAi+OKcHqP7Oxp1epuN8IgcDjYH/Ao8AYlckyx9aBKKfG7ZbWthtaTSiz4+WAHLu5vY19rO7DiSkRElEW2NQRx6bxquIPql62fHFOIs49qe9xTtjAbBLrbjNjtbqm07nGH0C8/tug0Y6sHX2p6WbvaDLikvx2XldkxsICtAR3B4EpERJQl6vwKLp5bhSpNFfDhk524lC9XRyjNVQfXnY2xBdfGgIK7v1ZPDTijpxX/nVgMk4HV1c7g9jUiIiKdBRSJGVvceHptI/a54zsftDVBReLqBdXYWKce6/S/I/Nw3dDYDhfINpEbtGIbifXY6gbsdbf8cmA2NP9ywNDaeay4EhER6ajer+Dy+dX4Ym/z6PK3N7sx7+yShIeaO5fXYf4e9bj0C/vZcOfIvIQ+bjrTBtcdMYzE2u4ReHJto+razcNyMcDJ1oB4YMWViIhIJwc8IZzzaeXh0AoA31UFsPxAYkeZP7+uEf/a0KS6dmKJGU+eWsiNQUfQSzNZ4PVNbjQEWp8QIKXEo5stqk1vPewG/G4EfzmIFwZXIiIiHWxrCGLqJxVYXRWIeN/ivYk7OHLuLi/uWKHutyx1GPHGGcUx75DPVhNLc2AM+xTtdofw0Kr6Vm8/a4cXX9Wqw+79JzqRy4MF4oafSSIiogRbUx3AlE8qsKUh+kvNi/clJriuOODD1QvUs1pzTQL/nViMLjZj63ckAECfPBNuGqbu/312XRNWV0VWyD1BiTs1vyCM6WbB+X05qSGeGFyJiIgSaOk+H86cXYH9ntZfYv76gB+eYPRTlTrqu0o/LpxThcawj2sQwL/HFWFYEfstY3X7yDxVy4AigVuW1iIU9tuANyhx4+IaVQ+sUQCPjCpgK0acMbgSERElyJydXlwwpxL1fnUovbCfTbXxx68AK+LY57qmOoDz51SiXnPE6IMnOTGlFw8YaA+H2YBHRxWorq2qDODfB3uGq70hnD+nEu9v86hu88shDh7jmgAMrkRERAnw/lY3fvZ5Fbya7oDrhjjw/GmFOK27VXV9SZzaBTbVBnDeZ5Wo8alD6x0j83A9x151yJReOfjJUerAf9/Keny5z4dJn1Rg2X71Lx2lDiPuODZfzyVmDQZXIiKiOJu5z4hfLKqB9tX/u4/Lx19PdsIgBMYmILhurQ/i3M8qUelVtyXcMjwXt3PsVaf89eQC5JlbXvZvCEicNbsSm+vVv5n0sSn4eJoLTgsjViLws0pERBRHT69txP0/WlUbogDg0VFO3DYi73DP49huFtX7v6nwwx1svQ+2LbubmkdthQ++B5orvH86Pp+9lp3Uw2HE3ccduYp6ajcL/j3Ciz55HJOfKAyuREREcfLwd/X4g2ZnuVEAz51WiGuGqF+mL801oW9eS59rQAGW7+9Yn2udX8H/zKnEriZ19e+qgXb89WQnQ2ucXDPYgWNd0ftWL+pvw4zJLsRwIix1AoMrERFRHLy4oQkPrWpQXbMYgFfGF+Hi/vao99G2C3RkLJY/JHHl/Gqsq1UfR3pJfxv+fgp3tceT0SDw+OgCaA85+/2IPDw3thBWIz/XicbgSkRE1Elf7ffh9uW1qmt2k8Dbk4px9lGtz/Ec000TXNt5EIGUEr/6sgaLNPc7s3cOnhxTCANDa9yNdFnwj1MKYDMKlOQY8PSYAtx1HFsx9MKCNhERUSfsdYcwfUG16phPm0Hi/ckunNzV2vodEVlxXVkZQGNAifmkpQdWNeC/m9VjmE4sMeOF0wth0pYFKW6uHOjA5WV2hCRg5udZV6y4EhERdZAvJHHl/KqIwwX+PNDfZmgFgO52IwaENUWGJPBVjH2ur2xswqOr1a0J/fKMeHNiMewm/vOeaAYhGFqTgN/ZREREHXT7V7X4uiKguvbbY3IxwRX9aNdoxnZXTxeIpV1g5jYPfrtM3ZpQbDXg3ckuuHJ4lCtlLgZXIiKiDnh5YxNe3uRWXZvY04q72jl4fqy2z/UIG7SklHh0dQOmL6hGKGzcls0o8NbEYvTjlnbKcPwOJyIiaqeVFX78/it1xbNPnhEvnF4EYztfPj5VE1y/qwqgzq9EDLD3BJs3Yr27Rd3TKgD86/RCnNhFXbklykSsuBIREbXTw6sbVJuxHCaBNyYUo8Da/n9Wu9qNGORsqSMpEli2X1113esO4azZFRGh1SSAJ8YUHHFyAVEmYXAlIiJqBykllmuC5ROnFmBYUfTB9LGIOP51b/MGrZAi8dlOL8746ABWVqp7aQssAjOmuHB5maPDj0uUbtgqQERE1A573Apq/S0Npg6TwHl9O1fxHNvdihc2NB1+e95uL3JXCbxe7o44DQsABjlNeJM9rZSF+B1PRETUDutq1JXPoYWmTg/6H9NN3Z+6oTaIv37XEPW2k3pa8cK4oogeWKJswO96IiKidlhbrQ6uwwo73iJwSHGOEUML264l3TjMgbcmFjO0UtZixZWIiKgdIiuunQ+uQPNYrHU1wYjrNmNzK8LPBzk4OYCyHoMrERFRO6zRBNfObMoKd80QB17a2AT/wWkFw4vMmD7Qjgv72Ts0rYAoEzG4EhERxcgfkthUq66KxqNVAADKnGasuKArvtznw7BCM0YUmyE62TtLlGkYXImIiGJUXhdEMOzEqp52Y1yroX3yTOiTx3+aiVrD1x6IiIhitDaiTYAhk0hPDK5EREQxStTGLCKKDYMrERFRjBIxCouIYsfgSkREWUdKie8q/ZixxY2GgBLz/bTjquI1UYCIYsPmHCIiyireoMSty2rx5o9uAMBxLjPmnFUCk+HIO/hrfAp2u1uOXzUbgDIn/xkl0hMrrkRElDX2NIVw1uyKw6EVAFZWBrBwj6/N+2o3Zg10mmBuI+wSUXwxuBIRUVZYvt+HcR8dwLeVgYj3LdnXdnBdx/5WoqTjaxxERJTxXtnYhNu+qkVr7ayL97a/4sr+ViL9MbgSEVHGqvSG8IcVdXh7s+eIt/uuKoB6v4J8S+svRGpHYbHiSqQ/tgoQEVHGUaTEq5uacMJ7+6OG1t8dk4sB+S21m5AElh/wH/HjaScKcIYrkf5YcSUiooyyoTaAW5fWYtn+yCBqNwk8PaYQ5/W1ocYn8WN9SxhdvNeHSaU5UT/m9oYQmsLOei20CnS3s/ZDpDf+X0dERBkhoEg8sLIeY2ceiBpahxaaMOesEpzX1wYAGNPNonr/kTZoaftbhxaaIQQnChDpjRVXIiJKe9sagrhmUTW+qYicGGAzCtx5bB5uGJarGl81prtVdbsj9blGbMximwBRUrDiSkREae29LW6cNvNA1NA6udSKr87vgl8Pz4uYudrFZsSgsAMEFImolVogcmPW0ZwoQJQUDK5ERJSWGgMKblpSg18sqkF9QKre19VmwCvji/DficU4Kq/1FxfHaqqurbULrK3mxiyiVMDgSkREaWdzXRDjPqzAG+XuiPdN6ZWDped1wbl9bG32oY7p1nZwdQcVbA7bxCUADC5gpx1RMvD/PCIiSiubagP4yaeV2OdRnyZgNQL3neDEtUMcMW+cOlWzQWt1VQB1fgXOsD7XjbVBhNdz++QZkWtm3YcoGfh/HhERpY31NQGcNTsytA50mjDv7C745dDcdu32L7EZMaRA2+eqrrqu4VGvRCmDwZWIiNLCmuoAzp5diQqvOrRe1M+GBeeUYHgHN0xp2wUW71Vv0Io4MYsbs4iShsGViIhS3uoqP875tAJVPnVovXqQHc+eVghHJ166147F0va5ruWJWUQpg8GViIhS2tcH/PjJp5Wo8aknB1w72IG/jy6AoZMHAWj7XL+vCqD2YEAOKRJrNa0CRzO4EiWNbsFVCHGTEOJ7IUT9wT/LhBBn6fX4RESUXkKKxOPfN2DarArU+dWh9YahDjwyyhmX06tcOUYMDetzlQCW7vfBG5SYvqBaVeW1GQX65Bk7/ZhE1DF6ThXYBeB2AOVoDszTAXwghDheSvm9jusgIqIUt7sphOu+qMaSfZEHAvz66Fzce0J+XI9cHdPdinW1LS0Bs3Z48dTaRnypefxJpVYYDTzqlShZdAuuUsqZmkt3CSFuADAaAIMrEREBAGZu8+A3X9agVlNlBYDbRuThrmPz4hpageYNWs+vbzr89utR5sOWOoy470RnXB+XiNpHSBn5gyHhDyqEEcD/AHgVwPFSyh/C319XV3d4UeXl5TqvjoiI9KZIYEWtATP2mbCgKrKm4jRJ3F3mx7jiUEIevzYATFpub/X9fe0KnhjmQ1er/v9mEqWbsrKyw393Op1x/S1T1wMIhBDDASwDkAOgEcD52tCqFf7k0015eXlar7+z+Pz5/Pn8+fzbsqMxiDfK3Xij3I1dTdFD6endrXj2tEJ0tye2t3TYpv0REwQA4KQSC/47qRiF1ti3hWTz1z+bnzvA55/ogqPeJ2dtBDASQAGAnwJ4RQgxTkq5Rud1EBFREu13h3DrslrM3uFFazVMswG457h83HR0bqcnB8RiTDdrRHCd1NOKl8cXdWrcFhHFj67BVUrpB/DjwTe/EUKcCOBWAL/Qcx1ERJQ8G2sDuHBuFXY2tv6y/9BCE54eU4iRLkurt4m3SwfY8cKGJoQOJumL+tvw1JhCmLkZiyhl6F1x1TIAsLZ5KyIiyghf7vPhss+rom68shkFftInB1cOdOCUrpa4b8Bqy0iXBa9PKMLHO7w4vbsVF/az6VLpJaLY6RZchRB/BfAJgJ0A8gD8DMA4AJzlSkSUBWZsceP6xTXwqw+/wtFFZvx8kAM/7WeD05Lcl+Sn9bZhWm9bUtdARK3Ts+LaDcDrB/9bh+YRWNOklJ/puAYiItKZlBJPrmnEH7+pj3jf9IF2PDa6ACa+HE9EMdBzjutVej0WERGlBikl/rCiDs+sa4p4393H5eN3x+Tq3hJAROkr2T2uRESUoUKKxK3LavHqJvUwf5MAnhxTiEsGtD43lYgoGgZXIiKKu6AErl9cg3e2eFTX880Cr00owuk9cpK0MiJKZwyuREQUV76QxJ0bLFhYpQ6tJTkGzJjiwvAic5JWRkTpjsGViIjixhOUuGJ+FRZqjm3tYTdg5lQXypwMrUTUcQyuREQUF+6ggovnVmHxPr/q+lG5Rsyc6kKfPP6TQ0Sdw58iRETUaf6QxPT51RGhtcxpwgdTXOjpMCZpZUSUSRhciYioU0KKxHVf1GDubp/q+tDC5tDaxcbQSkTxweBKREQdJmXzyKv3t6k3YpXZFXw81YWiHIZWIoofBlciIuoQKSXu/ro+Yk5rvzwjnhjiZmglorhL7qHQRESUtv62ugGmam65AAAgAElEQVRPrW1UXSt1GPHBVBeKLUlaFBFlNAZXIiJqt5c2NOHBVQ2qa64cA96fUozeuXwxj4gSo8PBVQjRRwjB36mJiLLM57u9uO2rWtW1fIvAjMnFnNNKRAkVU3AVQtwnhLjy4N+FEGI2gC0A9gkhTkzkAomIKHWsqwng6gXVCMmWa3aTwNsTi3EM+wOIKMFirbheCaD84N+nADgRwBgA/wHwcALWRUREKeaAJ4SL51WhPtCSWgWAf51WiFFdrclbGBFljVgbkboC2HXw72cCeFtKuVQIUQFgRUJWRkREKcMdVHDpvCrsbAyprt9/khNnHWVL0qqIKNvEWnGtBtDr4N+nAJgfdn/OOyEiymCKlLhhcQ2+rQyorv98kAM3DnUkaVVElI1irbjOAPCGEGIjABeATw9eHwngx0QsjIiIkq95VmsdZm7zqq5P6GHFI6OcEEIkaWVElI1iDa63AtgB4CgAd0spDw3uKwXwXCIWRkREySWlxF1f1+HptU2q60MKTHhpfBFMBoZWItJXm8FVCGEGcC+AZ6WUO8LfJ6V8LFELIyKi5JFS4o7ldXhuvTq0luQY8N9JxXBaOAaciPTX5k8eKWUAwK9juS0REaU/RUr8/qvI0FpoFXh3Mg8YIKLkiTWMzgMwLoHrICKiFKBIid8tq8ULG9ShtchqwIdTSzCCs1qJKIli/bV5NoCHhBBHA/gWgOonmpTyw3gvjIiI9KVIiVuW1uLVTW7VdVeOATOnuDCsiKdiEVFyxRpcnzn4399GeZ8ER2IRUQar9Sn40zd1OOBR8JvhuRk5bF8erLRqQ2sXmwEfTnVhcAFDKxElX6zBlT+xiChr3bSkBp/saB4HNXeXF6+ML8qooftSSty+vA4vbVSH1m4HQ+tAhlYiShEx9bhKKUNH+pPoRRIRJcvW+uDh0AoAQQlctbAas3d4kriq+Gme01qP5zUbsbrZDPh4GkMrEaWWmLeGCiGcaD41qzcAVXe+lPLBOK+LiCglvFbeFHEtoADTF1TjtQnFmNIrJwmrig8pJe79th5PrW1UXS/JMWDmVBcGOBlaiSi1xBRchRAnonmDlgKgCMBeAN0AeAHsBMDgSkQZJ6BIvFHujvo+vwJcMb8Kb5xRjEml6RleH/quAf/4QR1ai6zNoXUQK61ElIJiHYf1GIC3AHQF4AFwOporr98C+EtilkZElFyf7vRiv0c5/LZ25r5fAS6fX4VPtnsgpdR5dR0npcT9K+vxyHcNqusFFoEPphRjaCFDKxGlpliD6zEAnpDNP5lDAKxSyr0A/hcMrkSUoV7dqG4TuKzMjsdHF6iu+ULAZfOrcfx7+3H/ynqsrwnoucR2k1LiDyvq8OhqdWjNtwh8MMWFYzinlYhSWKzBNYDmsVcAsB/N1VYAqAdQGu9FEREl287GIObt9qmuTR/owNWDHXhstDPi9lsaQnh0dQNGf3AAp7y/H2/uNiGopFYVNqQ0z2l9Zp06kOeZBd6b5MJIF0MrEaW2WDdnrQJwAoBNABYB+IsQwgXgCgDfJ2htRERJ81q5G+Gxc0Sx+XCw+8XgXCgS+P1XdVHvu642iHW1FhjyG3D3cfk6rLZtQUXixsU1eHuLehqC0yIwY7ILx5cwtBJR6ou14no3miuth/5eD+BfaN6gdV0C1kVElDQhReINzSD+6QMdqrevHZKLj6a6cFbvnIje10Pe+tGdEr2vvpDEVQuqI0KrK8eAj6eVMLQSUdqIqeIqpVwR9vcDACYlbEVEREk2b7cPu90tI6rtJoEL+0UeODC2uxVju1tR61PwyQ4PZmz1YP5u3+FK7a6mELY3htAnL+bJg3EXVCSuXliNWWGzaAGgh92AD6ZwTisRpZdYK64AACHESCHET4UQ9oNvW4UQ7foYRESp7pVN6h7Q8/vakN9aWRVAgdWAy8oceG+yC2O6qauXX+z1tXKvxJOyuadVG1qPyjVi1pklDK1ElHZiCp1CiC5CiC8BrATwNppbBADgSQB/T9DaiIh0t9cdwmc71UFv+kB7zPcf292qenvJvuQF1/tX1uN1zRzaMqcJs84sSWoVmIioo2Ktlv4dQDWALgDCfwq+jebTtIiIMsIb5W6EwtpShxaYcGI7ekDHdNME172+pPS5PruuEY99rz5coFeuETOnuNDTYdR9PURE8RBrcJ0E4A9SykrN9R/RMhqLiCitBRWJ1zRtAlcOckAIEfPHOL7EApux5fZ73Aq21IeOcI/4e3eLG3csV088KLYaMGNyMXowtBJRGos1uNrQfLyrlquV60REaefFDU3Y3tgSMq1G4OL+sbcJNN9H4OSu6grtYh3bBebv9uKGxTWqaw6TwDuTilHmZE8rEaW3WIPrEgBXhr0tD27K+j2ABXFfFRGRzg54Qrh/Vb3q2oX97Ci0tn//aUS7gE7BdeEeLy6fX41Ayym1MAngtQlFOI4jr4goA8Tanf+/ABYKIU4AYAHwCIBhaK64npqgtRER6ebP39Sj3t/Si5prEh0+PGCsZrLA4oN9ru1pOWivz3Z6ceWCKvg0XQnPjC3EhJ45CXtcIiI9HbGUIIQYBQBSyjUAjkHzVIEFAJwAPgRwrJSyPNGLJCJKpOX7ffjPj+rd93ccm4fu9o71gx5XYoHN0BKC93sUlNcFO7XGI/louweXz48MrQ+d5MT/tLPVgYgolbVVcf1CCPEwgHullHsA3KXDmoiIdBNSJG7THN06pMCE64bmdvhjmg0CI/IVfFXbEnwX7/MlZG7qu1vcuO6LGtUkBAB48CQnbhjW8edARJSK2mremgbgCgArhBBDdVgPEZGuXtzYhB+qA6prj4wqgNnQuZf1j3eqy59L9vo79fGieaO8Cdcuigytj4124kaGViLKQEcMrlLKzwEMB7AKwDdCiN/qsioiIh1UeEK4b6V2Q5Yt4hCBjjihQFG9vWRf/Oa5hhSJ+7+tx01LahH+EQ0CeGpMAX4xmKGViDJTm9tlpZQNUspfoLny+ogQolEIUR/+J/HLJCKKvz9/G7kh674TnXH52INzFeSZW6q2FV4FG2rVfa47GoOYuc2DCk/sc16rvCFcOLcKj37foLpuFMC/TivEZWWOzi2ciCiFxTRV4OA0gfsBlAN4FEDidhkQpbmQIrFgjw9bq4zo3U/CakzcTnLquLd+dOMNzXGot3diQ5aWSQCju1owZ1fLKKzFe30YUtjc5/rZTi+uWlANT0ii2GrAzKkuHF105B7Ybyv8mL6gGrua1EHXbABeHFeEc46yxWXtRESpqq2pAiYhxH0AlgL4DM1TBP4tpXwl/I8uKyVKE3/5th4Xzq3C79dbcc2i6qQc90lH9szaRlyvGdI/uMCE6zuxISua1ua5Lt/vOxxaAaDKp+CSeVU40ErlVUqJFzc0YeqsiojQ2tVmwAdTXAytRJQV2moV+BrNBw9Mk1LeIqXkKVlER9AYUPDs+pbz4T/a7sWivfqdmkRHJqXEAyvrcecK9RQBgwAeG935DVla2l7ZJfv8WFcTwMXzqg6H1kN2NYVw+efV8AbV171BiZuW1OK3y2pVBwsAzRXdRT/pglO7db4nl4goHbQVXNcAOObgJi0iasOCPb6IWZp/W90Q/cakq5Ai8btldRFfD7MBePH0ooSEv2OKzMi3tIThap+CabMqUOuPXoVfUeHHr76sOVyl39UYxLTZFREzZgHgxmEOfDjVhW5xam0gIkoHR+xxlVJeoddCiDLB7B2RL0p8uc+Ppft8OIVVsaTxhySuX1yDGVs9qusOk8DrE4owPkEnSxkNAqd0teLTnS3fF3Wa0NrTbsRud8tvO+9s8WCg04RTullx1YJqVHjVZVaHSeDJMQU4vy8PFiCi7NP+Q7g7SAhxpxDi64OTCCqEEB8JIY7W6/GJEi2kSHy2M3o3DauuyXXPN3URobXQKjBzqithofWQMZrjX8P9fJADX57XBWVOdQ3hgVUNOPfTyojQOiDfhM/PKWFoJaKspVtwBTAOwNMATgEwAc2TCeYJIYp0XANRwnxd4UeVT4n6vgV7fPimIv4D6Klt2xqCeGF9k+paD7sBs88swQklrYfKeGltJuy5fXLwt1FOFFgN+O/EYhRa1f21mlZXTOmVg8/PKcHgBJy+RUSULnQLrlLKKVLKl6SUa6SUP6B5LmwJgFP1WgNRIkVrEwjHqmty/HVVvSoEljqM+PQs/QLg8CIzCizqUDq2mwXPn1YE48HNYP3yTXh1fDFMrewN+/2IPLx5RhGcFj1rDUREqSeZPwXzDj5+TVs3JEoHszVtAud1VY87/mynF6urWHXV08baAN7eom4R+OPx+eidG9MI67gwCKEas3Wsy4w3ziiOmO87trsVfz+lQHUt92AP7l3H5cMgOA+YiEgka8akEOJtAGUATpBSqvZh19XVHV5UeXm53ksjarftHoELv22Zo2kSEnNP9uDXa634oaFl1/f44iAeGcLwqpc71lvweVVLSO1rV/DmsV7ofSaElMDyWgPqgwLji0MwH6Fk8O5eE/6z24SeORK/7edHXzvnABNReikrKzv8d6fTGdefuPqVHcIIIf4OYAyAMdrQqhX+5NNNeXl5Wq+/s7Lp+X/6QwOAltOPx3bPQa7Jgz+e3AUXzas6fH1BlQkBVw8MLWx5mTqkyMMvGWeSZH/9V1f58XlVheran092YXAffQb1a5//wBjvd2cZcGdilqSrZH/9ky2bn382P3eAzz/RBUfdg6sQ4nEAlwAYL6XcovfjEyXCLE2bwLRezTvVJ5VaMbLYjO+qAoffd/WCarhsBhzwKNjvDqEhIDHSZcZr44tQquNL2JnuwZX1qrdHFJvxk6MSO0GAiIgSS9ceVyHEPwH8DMAEKeUGPR+bKFGqvCEsP6B++X9q7+aAJITAbSPyVO/bWBfEl/v8KK8Loj4gIQGsqgzg+sU1PB42TlYc8OGzXeoTy+4+Lh+CfaJERGlNzzmuTwG4GsClAGqEEN0O/onv4eBEOpuzywclLG8eXWRWbf45s3cOhha2XUldss+PWW1MJqDY3Petuto6qosFE3vyAAgionSnZ8X1RjRPEvgcwN6wP7fpuAaiuJu9Q71r/VCbwCEGIfDQSc5WRx2F++PXdfCHWHXtjEV7fFi8T10Bv4vVViKijKBbQ52Ukv9qUMbxBiU+361+SfrM3pF9lKf3yMGCn3TB8v0+OMwGdLUZ0NVmRGNAwbTZlYcrtlsaQnhhQxNuHMYXIjrqAU1v67ge1lYPASAiovTCnSBEnbBknw9NYdPtu9sNGFEcfbD98CIzhhdFvu+qgQ68uLHlZKdHvqvHpQPsKLRy2Hx7bWsIYoXmhLI/HpefpNUQEVG88V9Gok7QHjowtVdOuwfF33lsHvLMLfep9Uv8dVX9Ee5BrfmhOqB6++QuFhyvw7GuRESkDwZXog6q9SmYFdHf2v4ZoSU2I353jHrywL83NKG8LtDKPTLHyxubcM7sCjz8XX1cJiqsq1F/zo5ppfpNRETpicGVqJ2klHh7sxsnztiPvW7l8HW7SeC0DvZSXj80F71zW07YCkrgnq8zu+o6b5cXtyytxeJ9fjy0qkHVLtFR2uA6rJDBlYgokzC4ErXDptoAzv2sCr/8ogYVXkX1vsmlOciJZXRAFDkmgXtPUPdizt7pxaI9vlbukf6eWNOoevu5dU2drrqurwmq3o5lDBkREaUPBleiGAQViQdX1ePUmQfwxd7IMDkg34Q/Hd+5TUDn9bHhJE0/5n0r6zr1MVPVptoAFmk+j5vqglihOcihPbxBic316uA6hBVXIqKMwuBK1IY6v4KL51Xhke8aEFAXWWExAHeMzMOSc7ugb37nqntCCDx4slN17ZuKAOr9Siv3SF8vbIjeFvBaubvDH3NjXQDhI3B75xqRZ+aPOCKiTMKf6pQVFCnx2qYmXLOoGh9s9bR9h4O2NQQx5ZOKiFmtADC+hxXLzuuKO47N73CLgNYJJRb0yzOqru1sDMXlY6eKhoCCt36MHlDf3+pBg/a3gxiti2gTYLWViCjTMLhSxpNS4s7ldfjVl7V4d4sHVy2sjpgGEM3y/T5M/LgCG2rVgairzYAXTy/EjMnF6O+Mfw9laa76Y+5qyqzg+vZmN+oD0XtZm4IS77fjF4tw2o1Z7G8lIso8DK6U8R5d3YDn1qtfmr57xZGPVn1nsxs/+awSlZoNWMe7zPjiJ11wQT97wo4Q7ZWrrrjuagq2csv0I6XEC5qvhStH/WPotU0dmy4QGVxZcSUiyjQMrpTR/r2hEQ+saoi4vqUhhJdaGb/0/LpGXPtFDXyaQuf5fWz4eFoJutqNUe8XL6WOzG0V+HK/H+vDKthGAbwyvkh1m68rAthQ2/4ZtgyuRESZj8GVMtb7W924bVnru/If/q4BtT51RXX5fh/uXBF5n9+PyMO/xxXCFqde1iOJrLhmTnDVVlvP7J2DU7tZMbabeprCa5vat0mrxqeoZuqaDUBZAto4iIgouRhcKSPN3+3FL7+oQXgzgM0o4AgLntU+BX//vqUaW+UN4ecLa1Q70y0G4PnTCnHXcfntPsq1o3plaMV1T1MIH29X969eMzgXAHDlQIfq+ls/uo/YyqGlrbaWOU0wG/T5ehERkX4YXCnjrKzw4/L51arRVUYBvDy+CLcMz1Xd9tl1jdjeEIQiJW5YXIPdbnVIfHl8ES7qb9dj2YeVOjSbszIkuL68qQnBsCw6yGnCad2bK61nH2WD09ISNKt8Cmbv9Mb8sXliFhFRdmBwpYzSGFBw9cJquIPqat3TYwsxpVcObjo6Fz3sLd/2fgW4b2U9/u+HRszZpR559aujc3Fmb5su6w7XU1Nx3esJIaB07kSpZPOHJF7R9BRfM8RxeIObzSRwUT/1Lwjt2aTF/lYiouzA4EoZ5S/f1mO7pkL515OduPhg1dRuMuDu49QnXL27xYP7Vtarrp1UYsE9nTwJq6NyTAJdbC3/ayqy+WX2dPb+Ng/2e1pK4LkmcfhrcsjlA9Vvf77bh12NsU1U0M5wHcJRWEREGYnBlTLG0n0+PK/Z/HPtYAeuH6puD7i4vx3Di9QVufB2ykKrwIvjCpPaIxnR55rGwXXxXh9uXVqrunbJADvyLeofPyOKLTgm7OsiAfynlYMKwkkpsZ4VVyKirMDgShnBHVRw85Ia1bU+eUb8+YTIqqnRIHD/ia1XU58dWxRxCIDeSrWTBdK0z3XRHh8umlulat0wiOY2gWiu1FRd3/rRDSmP3CaxsymkOtAg3ywigj8REWUGBlfKCA+sbMCWBnW4e+LUQjhaOav+9B45mFJqjbj+m6NzMaVXTkLW2B69NBu0dsb4knkqWV5rwMXzKuHRTAe49/h8DC6IXhG9sJ8d4YXYLQ0hbKo78nPX9rcOKTQn7HAIIiJKLgZXSnsrDvjw9NpG1bVfDHZgbPfIYBru3hOdCM+1o7pYcHeS+lq1IiquadYqMH+3F79bZ4VXs+y/nJCPXw3Pa/V+BVZDxNft0zamC6zX9LfyqFcioszF4EppzRuUuHlJrWpea6kjeouA1uACM96ZVIwzelpx5UA73ppYnDKzP7UvdadTcP18txeXfl4Fn6L+XN5/Yj5+fYTQesg0TcV79o4jB1dOFCAiyh4sTVBaqvEpWLTHhzd/bIp4Kfn/Ti1AXistAlrjeuRgXI/ktwZoaSuuyT6EYG11AO9scWOvO4SfD3Lg5K7Rq9nldQFMn18dcVzuQyc5ccOw3Kj30ZraKwe3fdVyetnyA35UekNw5UTvW13L4EpElDUYXCktSCnxXVUAs3d6MX+3FysrA4g22vSKMjsm9Ey9INpevTWbw3Y1hSCl1LV3syGgYMYWD17d1IRvK1vC4btbPHhtQlHEjNvGgIIr5lejUTND95GTnfjl0NhCKwCU5ppwTJEZ31c3P6YE8NlOLy4ri9zQFVAkyuu0rQIMrkREmYrBlVKWlBJraoJ4f6sbM7Z6sK3hyFXHHnYD7j/JqdPqEqvA0nw8bdPBEOgOSlT7FBS3UnWMp2pvCH/+th7vbfEcfvxwIQlcvbAa7012YUy35sqrlBK3LK3Fhlp1iHzgpPaF1kOm9s45HFyB5j7XaMG1vC6oOiGtu92AQis7oIiIMhWDK6UcT1Di2XWNePNHd5s7yg8pthrwwulFcFoyI7QIIdAr16gKgjsbQwkPrg0BBVNmVUZUMbV8IeDSeVX4aKoLI10W/Gt9E97d4lHdZlpJEDcOjT72qi1n9srBI981HH57/m4fvEGJHJO64sz+ViKi7MLgSillnzuEyz6vUr003ZphhSZM6JmDM3paMaqLNSLUpLtShya4NoUw0pW4x5NS4pYva1sNrSeUmPFNRcvXpSEgceHcKtx3ohN3fV2nuu3QAhPuHODucGvDiGIzutsN2OtuLqc2BSUW7/NhUqm6DYTBlYgouzC4Usr4vsqPS+dVY7c7ekuAxQBMKs3BWb1zMKFnDrrZM3vIfKl2skCCN2i9vNGN97aqq6Y97AZcVubAZWV29Mkz4Z8/NOBP37Qcj1vpVXDDYvXBD/lmgVcnFEEeUB+j2x5CCEztlYOXNracnDV7hzdKcGV/KxFRNmFwpZQwa4cH1y6qieipNAlgQk8rzu9rx5m9czKmFSAWvTQbtHY2Je4Qgu+r/LhjhfpY1mGFJsw7uwtsYZXs3wzPQ5VXwf+tadR+iMOeGluIAU4zyg90bk3TetlUwfXTnR48Jp2qKm5kxZU/0oiIMhl/ylNSSSnx5JpG3PNNPbTbgCb2tOLZ0wpbHYOU6fQ69rXer+DqheoRVrkmgZfHF6lC6yH3npCPGp+C18rdEe/79dG5OOcoW8T1jjituxV2kzh8XOwet4LVVQGMdFkANPfj7gj7nBgFMMjJiisRUSbLnvIVpZygInHr0lr8MUpovW6IA29NLM7a0ApEHkKwMwGHEByaBrC5Xv2x/3FqAcpaCYFCCDx+SgHOOUr9sv0pXS24J44nj+WYBCb0UM+LnR12itZ6TbW1f74p4/qciYhIjcGVksIblLhqQTVe3qSu2hkF8OgoJx4eVQBTipxilSx6VFxf2tg8aizc1YPsuLCf/Yj3MxkEXji9CFcNtCP3YMB8ZXxR3L9mU3tHP0VrR2MQv12m3hA2hG0CREQZjz/pSXf1fgWXfV6Fxfv8quv5FoFXxhVhfAYcIBAPPexGGAQOH7RQ4VXgCcqoL993xJJ9PtyxXN3XenSRGQ+eVBDT/a1GgX+cWoh/nFoYl/VEM6U0BwI4XJH/vjqAdza7ceeKOlR6FdVtjzvYQkBERJmLFVfSVYUnhHM+rYwIraUOI+aeVcLQGsZkEOihmZywO04btH6oDuBn86rgD8t+uSaBl8cVxi0Yx0OJzYiTuqgD6bVf1ESE1uFFZkwf2LGZsURElD4YXEk32xqCmPJJBVZXqXsTBxeY8NlZJRhUwI01Wr207QJx6HPd1hDEhXMqUR9Qdxb/89QCDEjBzU3Teh35l5lzjsrB7DNdKOCJWUREGY8/6UkXqyr9mPpJBbZojm09scSMWdNc6OnI3k1YR6Kd5bqjk32uFZ4QLvisEvs96orln47Px0/b6GtNFm2fa7jbR+bhlfFFyDXzRxkRUTZgjysl3Ptb3bhxcS08IXWF74yeVrw6vggOho5WxbPiWu9XcOHcqohfHq4f6sAtw3M7/HETbZDThL55RmwNW7fdJPDM2EKc2yc+o7eIiCg9MDFQwkgJ/HVVPa5eWBMRWn/a14Y3zyhmaG1DqUNzCEEHK66+kMQV86sj2jT+p58ND57k7PDRrHoQQuDOY1vGbJU6jPj0TBdDKxFRFmLFlRLCE5S4a6MFcysbIt5387Bc/OXEfBhSOCylioiKa2P7N2cFFImrF1Zj0V6f6voZPa14akxhWnwdLupvx0CnCdsbQ5hcmpNSG8iIiEg/DK4Ud/vdIVz6eRVWVqq/vcwG4O+jC3AFd3/HTNvj2t5DCIKKxLWLajBrh1d1/XiXGa+ML4LFmD4BcKTLgpGuZK+CiIiSicGV4mprfRDnz6nENk0fZZHVgNcmFOHUbtZW7knRaA8h2N0UgiJlTFXSkCJx45IafLBNfcBAmdOEtycVc0MTERGlHf7LRXHzQ3UAU2ZVRITWwQUmzD+nhKG1A/LMBhRYWkJqQEHERIBoFClx67JavL1ZHVr75hkxc4oLxVl8lC4REaUvBleKi6X7fDhrdgUOaELVxJ5WzDmrBH3yWNzvqNJc9eeuraNfpZS4/as6vKo5TrdXrhEfTnWhB0ePERFRmmJwpU6bvcODC+ZUot6vnhwwrSSINycWI9/Cb7PO6OXQjsQ68gatB1Y24F8bmlTXetgN+HCKC71y+QsEERGlL/4rRh2mSIl/b2jCHcvroJl2hRuGOnBVYQXMhvTZ/JOqtH2uRxqJNXeXF49+r57k0MVmwMypLvTN5//uRESU3vgvGXXI5rogfr20Bl/u80e874/H5eO3x+Tixx8rkrCyzNM7xskCe90hXP9FjepasdWAmVNcKEvBo1yJiIjai8GV2iWoSDy9thEPrqqHV5OfDKJ53NVVgzjuKp5iqbiGFInrvqhBla+lx9gogNfPKMKQQoZWIiLKDAyuFLO11QHc/GUNVlUGIt6XaxJ4ikdwJoS2LzXasa+P/9CILzQHDNx5bD5Gd+UkByIiyhwMrhSTj7d7cM2i6ogqK9A8OeDvpxSgNzf+JETEIQSa07OW7ffhoVX1qmundbfi1uG5CV8bERGRnpg0qE0vbWjC776qhaLZgFVoFXjopAJc3N+W0mfdp7suNgMsBsB/sAugzi9R71eQbzGgxqfg2kU1qs1xrhwDnjutEEZujCMiogzD4EqtklLir9814OHvGiLed14fGx4Z5UQXG2eCJppBCPR0GLE17GCH497bD0UCnqCERzPS4Zmxhehu59eFiIgyD4MrRRVUJG5bVouXNUPsDQJ4fAQD6EsAACAASURBVHQBpnMDlq5KNcG10hv99Kybh+ViUmmOXssiIiLSla6T4YUQpwkhPhRC7BZCSCHEVXo+PsXGE5S4ckF1RGjNMQJvTChiaE2CWCYDHOsy457j83VYDRERUXLofaRRLoA1AH4DwNPGbSkJan0KLphTiVk7vKrrhVaBmVNcmNabUwOS4Yahueid2/rL/8e7zHh1fBEsRva1EhFR5tK1VUBKOQvALAAQQrys52NT2/a5Q/jpnEqsrVHvWi91GPHe5GIMKuA80GTpm2/Cdxd2xc7GEIQAzAYBs6H5vxaDgM3EwEpERJlPSCnbvlUiHliIRgA3Sylf1r6vrq7u8KLKy8v1XFbW2ukRuHmNFXt86iJ8f7uC/xvmQxdrcr5PiIiIKL2UlZUd/rvT6YxrZSXlN2eFP/l0U15enhbrX13lx/Vzq3DAp97wM6qLBW9NLEaBtWMdJeny/BOFz5/Pn8+fzz8bZfNzB/j8E11wTPngSom1ZJ8PP5tXhfqAuqI6pdSKl8YXwW7Suw2aiIiIKDoG1yy2YLcXP/u8OmIO6CX9bXhiTCHMHGBPREREKYTltCw1Z6cXl3xeFRFabx6Wi6fHMrQSERFR6tG14iqEyAUw4OCbBgC9hRAjAVRLKXfouZZs9sl2D65aWI2AZob9n4/Pxy3H5CVnUURERERt0LviegKAVQf/2ADce/Dvf9F5HVnrg60eTF8QGVofPtnJ0EpEREQpTe85rgsB8DXoJHlnsxvXLa6Bopls9fjoAlw9mKdhERERUWrj5qws8eaPbty4uAbhmVUAeGJMAS4vY2glIiKi1MfgmgXeKG/CzUtqVaHVIIBnxhbi4v72pK2LiIiIqD0YXDPca5ua8Osv1aHVKIB/nVaIC/oxtBIREVH6YHDNYK9sbMJvltaqrpkE8O9xRTi3jy1JqyIiIiLqGAbXDPXihib8dllkaH1pfBHOOYqhlYiIiNIPDyDIQE+tbYwIrWYD8OoEhlYiIiJKX6y4ZhApJe79th7/+KFRdd1yMLRO7cXQSkREROmLwTVDBBWJW5bW4vVyt+q6xQC8PqEYk3vlJGllRERERPHB4JoBPEGJny+sxuydXtX1XJPA62cUYVwPhlYiIiJKfwyuaa7Wp+DSz6uwbL9fdb3YasC7k4txrMuSpJURERERxReDaxrb5w7hgjmVWFcTVF3vlWvE+5OLMcBpTtLKiIiIiOKPwTVNbakP4vzPKrG9MaS6PrTAhPemuNDdbkzSyoiIiIgSg8E1DX1f5ceFc6twwKOoro/qYsFbE4tRYOWUMyIiIso8DK5p5st9Plw6rwr1Aam6PqVXDl4aVwi7iaGViIiIMhNTThqZtcODn86pjAitl/S34fUJRQytRERElNGYdNLEf8qbcMX8anjVLa24cZgDT48thNkgkrMwIiIiIp2wVSANPLGmAX/8uj7i+j3H5+PW4bkQgqGViIiIMh+DawqTUuIv39bjcc0RrgYBPD66ANMHOZK0MiIiIiL9MbimqJAiceuyWry6KfII13+dXoRz+9iStDIiIiKi5GBwTUG+kMQ1i6rx0fbII1zfOKMYp/ewJmllRERERMnD4JpifCGJK+ZXYc4un+o6j3AlIiKibMfgmkJaC62lDiPen1KMMh7hSkRERFmMwTVFtBZaBzpNmDG5GKW5/FIRERFRdmMaSgGthdZBThM+nOpCV7sxSSsjIiIiSh08gCDJGFqJiIiIYsOKaxI1BBRMn1+N+XsYWomIiIjawuCaJPvdIVw0rwqrqwKq6wytRERERNExuCbB5rogLphTie2NIdV1hlYiIiKi1jG46uzbCj8umluFKp+iun6sy4y3JxajxMbQSkRERBQNg6uOPtvpxdULq+EOStX1iT2teHl8EXLN3CtHRERE1BoGV528u8WN676oQUidWfGzAXb889QCmA0iOQsjIiIiShMs8englY1NuHZRZGi97Zg8PDWGoZWIiIgoFqy4JtiTaxpw99f1qmsCwN9GOXHNkNzkLIqIiIgoDTG4JoiUEv/aYcLzO9Sh1SSA508rxAX97ElaGREREVF6YnBNACkl/vh1PZ7fYVFdtxqBl8cVYVpvW5JWRkRERJS+GFzjTEqJP6yowzPrmlTXHSaB/5xRhNN75CRpZURERETpjcE1zh5Y2RARWp0WgXcmFeOkLtYkrYqIiIgo/TG4xtGjqxvw6PcNqmvFVgPen1KMY4otrdyLiIiIiGLB4BonT61txP0r1Rux8k0SM6e6cHSROUmrIiIiIsocnOMaBy9uaMJdK+pU1/LMAv83zMfQSkRERBQnDK6d9NqmJvx2Wa3qmt0k8PakYgzLU5K0KiIiIvr/9s483sqq6uPfHyAos4KKKIhT4jxrTnkNnHPKCs05zWjQykrLnHqdGgzFNHtfozAtHHLKNOdAcSzTnGcRREBAQFQckPX+sfaBfQ/n3Hu53HPOPbC+n8/53OfZ41p772c/a0/PDZY9wnBtJWbGiKfmcuKDjY3WLh1hzJBV2HH1OIgVBEEQBEHQlsQe11awwIwfPzqH/3u+8dcDVugAV+3eJz55FQRBEARBUAHCcF1CPvrUGH7/LG6aMK+Re5eOMGq3VdhzQBitQRAEQRAElSAM1yVgzscLOPzemYyf+nEj956dxV+G9GGXfrE9IAiCIAiCoFKE4dpCnp/1CV8b+w7Pz57fyH2Nrh346x592SS+HhAEQRAEQVBRwnBtBjPjjy9+wGmPzebDTxv7bdCrEzfs2YeB3aMYgyAIgiAIKk1YXE0w66MFfGf8LG6b+OFiftutugLXDu3DKit2rIFkQRAEQRAEyx9huJbhgSkfMfz+WUz+4NPF/IattxIX7dSbrp3ia2JBEARBEATVIgzXDDPjvrc+4pKn32PclI8W8+/eSfx6p94MW69rDaQLgiAIgiBYvqm64SrpW8CPgDWAZ4HvmdkD1ZYj55MFxo2vz+OSp+fy7Kz5JcNs3XcFRu22Cuv0DFs/CIIgCIKgFlTVCpM0DBgJfAsYn/7+Q9LGZjaxmrLknPLIbP744gdl/b+3WXdO26onnTuqilIFQRAEQRAEOdXepHkyMNrMrjCz583sRGAK8M0qy9GIIzfoVtJ9l36duW2fvpy9ba8wWoMgCIIgCGpM1WZcJXUGtgEuLPK6C9ipWnKUYutVO7NLv86Mn/oxHQQHrL0SJ23ana1X7VxLsYIgCIIgCIIMmVl1MpL6A5OB3czs/sz9TOBwM9uw4DZnzpyFQr388stVke+RWR0Y905HDu8/n7VWqk6ZBEEQBEEQLGtssMEGC6979erVpkvWtThpVGwVqoTbQnLlK8kGwJFtnObLL79cNfnbI6F/6B/6h/7LK8uz/suz7hD6V3rCsZp7XGcAnwL9itxXA6ZVUY4gCIIgCIKgDqma4WpmHwOPA3sUee0BPFQtOYIgCIIgCIL6pNpbBUYAV0l6DHgQGA70B35XZTmCIAiCIAiCOqOqhquZXSupD3A6/g8IngH2NbM3qilHEARBEARBUH9U/XCWmf0W+G218w2CIAiCIAjqm2r/A4IgCIIgCIIgaBVhuAZBEARBEAR1QRiuQRAEQRAEQV0QhmsQBEEQBEFQF4ThGgRBEARBENQFYbgGQRAEQRAEdYHMrNYyLMacOXPan1BBEARBEATBEtGrVy+1ZXox4xoEQRAEQRDUBWG4BkEQBEEQBHVBu9wqEARBEARBEATFxIxrEARBEARBUBeE4RoEQRAEQRDUBWG4NoGkn0j6l6R3JU2XdKukTYvCSNLZkt6SNE/SWEmbFIX5qaQHJb0vqeTeDEkjJf1b0oeSJlRQrRZTLf0lbSFpjKRJKY0XJf1IUs3aZxV1X1XSnSmNj1IZXCapV6V1bIpqtv0sbF9JkyWZpL6V0KulVPnZtxK/4ZXUrzmqXf+SjpD0ZOr/Zkj6U6V0awlVfP6PKVP/Jmm7SutZjiq3/+0k3SNplqTZku6VtH0l9WuOKus/RNJDkuZKmiLpF5I6VVK/pmgL3SUNkjRK0mvJ/zVJF0haqSidgSn999Nzf4mkzs3JGIZr0zQAvwV2Aj4PzAfukbRKFuYU4AfAicB2wNvA3ZJ6ZGG6ADcCFzeRVwfgSqCmHXYRDVRH/22A6cCRwCbAWcCZwI/bSpFW0EB1dF8A3ATsD3wGOAYYAlzRRnq0lgaq1/YL/BF4cqklbxsaqK7+XwfWyH5XLr0KS0UDVdJf0knAr4ALgU2B3YFb2kqRVtJAdfS/lsb1vgZwNfA68O820qU1NFAF/SV1B+4A3kp57QhMAe4sSqfaNFAd/TcHbgfuArYCDgUOAH7ehrosKQ0sve6DgY7AN/F3+onAUcDIQgKSOgK3AT2AXYHDgC8Bv25WQjOLXwt/QHfgU2D/dC/8IftpFmYlYC7wjRLxv+RF3mQePwQm1FrXWumfhf0l8Hitda6R7icBU2qtczX1B74L3It3lAb0rbXO1dI/6fulWutYC/2B3sD7wB611rFW9V8UriswGzit1jpXqf63Te1/ncxtneS2ba31roL+5wNPFLntD8wDetRa77bQPQvzLWBmdr8PPnEzIHM7AvgQ6NmUTDHjumT0wGdGZ6X7dYB++GgJADObB9yPj1aWNaqpf88sn/ZAVXSX1B/4IjCu1ZJWhorpL2kr4FR8RL6gLYStAJWu/5FpqexfkoarhttkylAp/ffEZ2ZWl/ScfKvITZLWbRux24xq9X1fAbrhqw/tiUrp/yK+2nacpC6SuuCrDxOBZ9tA7raiUvp3wQ21nHnAivhKZHugrXQvfqfvCDxvZpMytzvxMmlS9/bWObZ3RuJLmQ+n+37p77SicNMyv2WJqugvaWt8yfzy1qZRASqqu3yP7wfAZHzkemwr5awUFdFfUjdgDHCimU1eWiErSCXr/0xgGDAUuAZfKjutdWJWjErpvy7+HjodOBk4GFgB+Kekrq2Wtu2pVt9/AvB3M5uyFGlUgorob2Zz8aXpYcAH6TcMn4GftxTytjWVqv87gR0kHSmpk6Q18f4AfNtIe2CpdZc0EF9N/m3m3K9EGjPw2d0myzAM1xYiaQSwC3CImX1a5F286Vol3OqaaukvaUN838vFZnZDa9Joa6qk+/eBrYGD8Jd5S/aEVoUK638J8GB7qetSVLr+zewcMxtvZk+a2a+BnwE/arXAbUyF9e+AG6onmdkdZvYYcDiwGr5kWnOq2Pdtgs9C1Xp/eyMqqX86rPMH3Cj6LLAz8ARwSxrU1pxK6m9md+EG3aX4zOtL+J5XcAOuprSF7pJWxw30u4GLmkmjOXcgDNcWIekifOPw583stcxravpbPDpYjcVHEnVLtfSXNBgYC1xjZrU8mLWQauluZlPN7AUzuwX4BnCCpAGtkbktqYL+Q4BjJM2XNB/f5wowVdJ5rZG5LanRs/8o0DN1+DWlCvoXZhafKziY2Rz8sM7AJZO27aly/Z8ATMIPK7ULqqD/V4H1gGPN7F9m9khyG4jPvteUatS/mY3A93oPBPqy6GDi60sscBvSFrpL6gf8E3gGONLSRtYsneI0+uJbh5oswzBcm0HSSPxB+ryZvVDk/Tpe+Htk4VfET8g9VDUhK0i19Je0MW60Xm9m318amduKGtZ94bnsspTpLBVV0n9PYAtgy/Q7Prk34LOxNaOG9b8lPvsyeynTWSqqpP+D6e+GWTrd8WXSN1ohdptRzfpPcY8E/mBm7WKfd5X074rPruU6L0huNbVPqln/5ryVtkcchg9g/tNa2ZeWttBd0hr4O/154DAzm1+UzsPARpLWytz2AD4CHm9SwFqfWGvPP+Ay4F38pHO/7Nc9C3NqCvNF/FMu1+CzBT2yMAPxl9EP8Qey8JLO01k/uY1I8QthOi/r+uOfy5iW4ub59FsOdP8CcHSKPwjYD599enh5aftF+TbQDr4qUMX63x8/jLIpPvN0PDAHGLk86J/C3IzPyOwMbAxcD0wAui4P+qdwR+BLwwNrWe81aP+D8UHa5cBG+LvgqvQMrLWs65/C/AjYLOl+BvAxcFA96w70x7c9jAUGFKXTMYXpCDwN3Id/CmwofsbjN83KWOsHpD3/UkMr9Ts7CyPgbHzJ60P8NPimRemMLpNOQxZmbJkwg5Z1/VP8knktB7oPxUees/HTpC8BvwBWXl7aflH4BtqH4Vqt+t8b39M3F/8s1NP4p8E6LQ/6pzA9gN8D7+Cnjm8F1lte9E/hxgG311LnGtb/HsB4vA+chS8t77Qc6X8fi/r/R4B96l13/HB1uXQGZeEGAn/HD+XNBH4DdGlORqXIQRAEQRAEQdCuiT2uQRAEQRAEQV0QhmsQBEEQBEFQF4ThGgRBEARBENQFYbgGQRAEQRAEdUEYrkEQBEEQBEFdEIZrEARBEARBUBeE4RoEQRAEQRDUBWG4BkEQBEEQBHVBGK5BEARBEARBXRCGaxAEQRAEQVAXhOEaBEEQBEEQ1AVhuAZBEARBEAR1QRiuQRAEQRAEQV0QhmsQBEEQBEFQF4ThGgRBEARBENQFYbgGQRAEQRAEdUEYrkEQBEEQBEFdEIZrEARBEARBUBeE4RoEQRAEQRDUBWG4BkEQBEEQBHVBGK5BEARBEARBXRCGaxAEQRAEQVAXhOEaBEEQBEEQ1AVhuAZBEARBEAR1QRiuQUkkjZV0fJXzPE3S75cwzsaS/t2KvH4n6Yx03SDpzSVNI1gyJL0pqaHKef5e0mlN+J8s6dxm0rha0hFtL11QaSR1kfSepP6tjD9c0j3pekVJJmmtMmGPk3RrS8K2Upa9Jb3SVukF1UfSYEnzay1HvROGaxWQNEHSvNSBvrc0HWm9Iml0cwaCmZ1vZsen8INSx9+pmaTPAS7M8pkg6WNJfYvyfzKlNyjlNdzMzmmNLkXp5nU7NenZvUS4s1P+2y9tnu0RSX+Q9LUq5zle0jFNhTGz483s/BR+qKQJRUF+Bxwrqc9SyHGCpBclzU1t4O+SurU2vaDtMLOPzKy7mb0FIOkaSadXKK9RZrZ/JdJe3pB0v6TPVSDditV/WxMDlfKE4Vo99k8daPe8I81pgZFWl0jqWKF01wB2B24u8nodOCwLtxmwUiVkSOxvZt2BLYGtgJ8UySngSOAd4OhKCFCpMm5h3gL2Av5Rpfw6SGqTvsvMPgDuwuunNbIMAX4GfMXMegCbAH9tC9mK8mkznYOgPb9rJPUCNgIeqrUsQfskOsIaks0qHidpInBfcj9A0rOSZqcl+42yOBMk/UjSU5LelzRK0uqS/pFmfO6RtPISyHBgmo18V9KrkvbOvNeW9GBK9658FlPS9Wl2aU4aHW+S+Y2WdLmk2yW9DxwHHA6ckmYmby0jy9mSrk6396e/s1OcHUtE2QP4j5l9WOR+FXBUdn808KeivMrOAEs6SdJzhWU+SV9IZTRb0kOSNi8Vz8ymAnfiBmzOrkB/4LvAoZI6Z3n9t2gm3pSW05ewjHcvocfakm6WNF3SDEkjk3sHSWdKekPS2ymtnslv/STDMfKl/XckfV3SDpKeTmUwsiirrYBpZjYlpfENSS+kdvOMpC2ysFundOZIGiOpS4rTJ+kyXdIsSbdKWjPTZbykcyQ9DLwPjAF2BH6Xyu3iMnV5dWpXvYBbgYFZWa+Wgo0F9isVvwVsBzxoZv8FMLOZZjbazN5P+fdOMkxPz+5PJCn5nStpdCbr+pKsCZ0HpnIaLWlKKqcbsvAHpPY0O8XdtCUKSDpP0kXpuoukDyUVZqm7p/tekv4s6bvJfe3UTk5I94OTjpL3GQcm94YUbs90v7fKbO2R1E3SJZImpfYxTlKn9LtB0rSk2z8lbZjFu0bSb5L7XEn3FtqOsiV7SScBhwBnpPq/PoU5U9LrWXttri0clOpyeiq7Qn0u3FZQQrfdk147p/tNJd2X6vB5SQdlYQ/Mnp9JSe48rdNS3pMlHZ65H5zq/135s31a5jdY0nz5szwJuD257yrp0VSu/ynIV0aHLTKZp0r6QXJfSdJlqU2+KelXklZIfntLekXS6fI+aLKkfZOOr0qaWUgnYy/gn2Y2X9LOkp5IOk2VdEFK915JXy+S76WUX0dJl6YympPKZMMm6n+ApFuSfK9JGp6l+XN5u782xXlS0jqSzkrhJ0harO/N4neSNDLp+Qr+zsr9877yFaVVK/kK0E3AulrUX/VJ5fFo0ustSRepHQ9CKoaZxa/CP2ACMLSE+yDAcKOqGz4r+Bn8JbUHsAJwCvAK0DlL6xFgdWBN4G3gP7jx0AU3fs/K8ngK+GoZubYH5qS8OqT0Bie/scCrSZ6V0v3Ps7hfA3qkPC8Gnsz8Rqd0d07prpjczm2mnM4Gri4qm05NhP8VcFmpsgZexEftHYFJwNopvUGZjOem6wbgzXR9RirPVdP91qmMd0hpHZ3y6FJct8BawNPAyCKZRgHXpfqcCXyxjD4nAC8APVtTxkVpdQKewbdRFNrWzlk+LwHrpPRvAf6Y/NZP5XRpyndfYB7eia6adJxZSCvFOR04J10flsp7G0B4+xmQ/N7E224/oE+S4fjktypwcJKzJ3Aj8Ncsj/GprDdK5dgpuR3TTJu6Gjg7XQ8FJpR5Dt5uJo0jyvg1pPI5C9ip0C4y/78kXXoA6+LP8tHJ71xgdBZ2fcCa0fnOlObKQGfgcynsdsC09Ldjajuvsqjf+F/gkjI67Ak8ka4/l+I9mPk9nrWbm9L1USncnzO/G9L1+cBF6frMFO68zO/XZeQYhc9+90s67Jr+dsKfu+54X3I58EgW7xpgNj6QWRHf/nFP8lsRb89rZWFPL8p3GLAG/hwdCcwF+ia/4SXSuhPojT8/rxXaRpmwawEHABOBrZJfT2AKPpjvmOrsHWD95D8T2D5d98ni7Q18Avw0tYeDk6zdk/8QfMa/A95vvQPsnfwGJ3l+D3TFn7NBKa+hKc6+wHRg5RJ1s3Ly+w7eL/QEtkt+vwQeAPri76V/AT8tkvnUVI8nAlPxyYVu+HvrQ2DNor6t8Iw8AXw5XfcAdsja37gszg4p3Y7AgcDDScYOqUxWK1X/KfzTSb7OeH81Edgt+f8c+ACfGOgEXIuv6P0w0+f5JvqO76X0++N93HhgfuZ/AN6OlOphHrBJVnavlOirCs/4enh/MrypPnBZ/NVcgOXhh7983sM719nAzcl9UOpM1s3CngFcl913ACYDDVlah2f+NwCXZ/cnFtJvgVz/S3rBlPAbW/SAfwu4o0zY3kmPXul+NPCnojCjaXvD9QoyYzorn6G4MXVBevjvTp1MU4brZGBE6lh6ZeldTjLKMrcXs46tULdzU/r3Ar2zsF2Bd4GDsjK/pYQuu+AG8mdaW8ZF4XcldeQl/MYBJ2T3mwAfpbZWMFxXz/znAIdk97cA38nuHwZ2TNf3At8uI9ObwKHZ/Qjg0jJhtwWmZ/fjgTOLwrSV4boR8HEzaZQ0XJP/fsDfUznNxQdUHXDjYn5ep8C3WWTctMRwPTO7H5DS61VChivIBqzJ7VWyAUYT8ndL9d8bf25OBd5Kbfc8YEQKtyFu6Ag3gE4AJia/PwMnpeu98JUQgHuA44Hx6f5B4IASMqyAGzgbtkDefsAC0mANN0byclwlteFVaYHhWiL9F4C90nUpY7QhC3sycFsTYU/FDZ3BWZyjgbuL8rwSODVdTwOOBXoUhdk7tbEOmdu7wJZl9PgdcEG6Lhiu/TP/s4AriuKMA4aVSOtY4OEy+UwGPp/dHwi8UErmVCcGbJGFf5ZFBrZwo75gaD6GG+p9SrTZd4GB6f7SrJ3um9LcHlBRvGLDdTfg5aIwPyO9U3HD9dbM78ukZ6BIn5XKlM1DZH0UbqjOLxU2+d8BfCMru1fKhU1hfgyMae6ZWdZ+sVWgehxkZr3T76Aiv0nZdX/gjcKNmS1I/mtmYaZl1/NK3C92OKgMA/CXWzmmZtcfFNJNSzE/T0s97+LGG/iIu0Cu02JIOjxbAmnt3shZ+Ci8FFcBXwWOoWibQBl64y/iC8xsTua+NvCDtJQ2W9JsvNzyw3UHme9vbMBfEHk5HIwbG7en+z8D+0hatRBA0gB8RvZoM3spuS1tGQ/AjbRPS/g1amPpujPeCQNgZi1qY2lJa13g0Szf1rSpbvIvAExM+t5HY12h+TZ1RtamLm0qbBE98AFlqzCz28zsC/is1BeBr+Mv+tXwmZHisl5zsUTKk+s8AJhR1D4LrA2cWtRO12hJXubbGp7AZ1s/hw9aH8FnMHfDjRnM7EXgY2AzfGD0N2CGpPXycLhxuklq45viRtm6qa1sg8/OFbMGPrh8rdgjLbdemJZx38UNS+GzkQUWlpOZvYMPJlt0AFa+VeuprNzWZ/G2l5PXyRvN5HMyPhh/IXNbG/hcUV0dgpcBwEHpfmJamt82izs9vRMK5M/QzvLtFdMlzcH7vlyPBdb4bMXawBFFcmxbRp+Sz7Uk4QOJptp4LvO89LfcO2tbfPXr7XR/NLA58FJaIt8LFrbZG4HD07aEYXifD77XfhQ+STBN0m9V4sBsVgaDisrg5KRTgWJZp1uyGjN9uskPfxb6n8eTe38Wby8LkW/veUy+JWs28HmaaHvyr+j8Q75t5l18RaOptrpMEoZr+8Cy67fwhwlY2DEMwEe1bc0kfLlhSfkqPqoeCvTCZ0fBXyYFrChOo3sz+7MtOqi2T4k8iuOX4il8aWfxyGZv4DMd++IdXHPMAr4A/FGN93lNwpc5e2e/rmY2pkSe4/CZ0Asz56PxTnmipKnA9fjs0mHg+8Pww2UXm1luwLemjHMm4XuUSx3aatTGgIG4QTK9ifTKsTc+e1R4MbW2TZ2CL5ltb2Y98Q68mOba1DlZm/pOC+IX2Aj475IKvFjiZgvM7G7c8NsUn0H/lMXLuvAsv4/PahbIX5YLk82uJwF9lfYjFzEJ+FmJdnpdC8Ufh7e1zYDH0/0+LG5o3g8cis8MT03hjkt6PI17vAc8q/gu3QAABnhJREFUCXwf397yCT6w+QE+EzerRP5T8AHeuiX8jsW3LOyOPwuDk3v+LAwoXEhaBX/mppRIq1EbkPQZ4Df4oHUVM+uNL7+qRNzF8sLrc7GDthkH48bh8MxtEnBXUV11N7PvAZjZw2kgtDq+dWKxvqYM1+FL2QPMrBfeFzXVX0wCfl8kRzczu6hE2iWf62TATaV8G19S9gVuy9J/3syG4YPAS4AbteiMwJXAEXgfNM3MnijIZGYjzGwr3OjdAj9fAKXL4IWiMuhhZgcvqeBmdk/W/2yTnKeweHsBfLCOvw/OwWeYe+MD9kKdleqvrsC3sq2X+sn/oem2ukwShmv74zpgP0lD0kjyB/gyXiVOWI7CPwU0RH5gZ01Jg5uN5TNUH+FLJl3xfWvNMY3SL6VyTMeXA5uKczd+2GfFMv7H4UtY77ckQzMbi+87u0nSDsn5CmC4/HCS0szgfpLKzfReDOwhaUv5AZEhuEG8ZfptAfyCRV8X+APecf6yKJ3WlHHOwynu+ZK6yg9QFAzyMcDJ8sOBPfDl4DFFMzktZT8WzSaDLyGfImmrVF4bpBnl5uiBzx7NSjNzZ7YgzpK2qWm44Vdcd7vRyi8iyA/EfEXSyknfz+KzkY8kg+2veB10l7QObswVDiA+CewmPxzSG1/2K4uZTcKX3i+TH/paQYs+GfR/wLclbZfk6C5pf7X8s1zj8Bm6p8xsPm58nwC8lGYw83DfYdHs6th0/0BR+ykXbhwlSGX1J2Ck/LBpR0m7pIFXD3wf5Ex8ibjUocoD0zPaJfn/M5u1yyluM93xfmY60CEZmOuXkjHjVPlhtUFJp2ubCDsR7wNO06LPxd0MbCVpWKrDzpI+K+kzqX85NA1OPsG3npRaNWlEmuDoDsw0sw8l7YQvazfFlcCXU//fMfURQySVGkDdDKwv6ZtJ3p6Stkt+Y4Cz5IeHVsOX9q8ukUZLaNSfSDpKUp+0cjQHN+YK7Wxs0vk8slW1VJbbyg8tvY8PygtlWFz/41Oc78kP8nWStLmkrVspfzHXAd+XtIb8cPMpmd9K+CTG28ACSQfgq3YFpgGrqfFscQ9gjpm9Jz+s2+iA2vJCGK7tjLQcdwQ+CzAD2B//3NLHrUlP/nWCw0v5mdlj+GzGRXinMI7GI+dy/Alf8pgMPIcvKzbHKGBj+XJM8eerSsn2Ad4hPZjifLZEmGn4CPXAMmm8amZL9M8J0ozZscDfJG2T4n8d30M1C5+NOaaJ+NPx8jkDP+jxpJndZWZTCz985mBz+anvQ4GD1fjLArvSijKWL5XfmuSYjxvMG+GzChOBL6WgV+Av2wfwpdm5LJqRaDHyzzMNwQ+rFPQfgxvm1+J70G7El9CbYwQ+mzYTH6S1xJC8GDgstY8RzQU2s2fwPeETUpzV5DPee9Oy7SSlmI3vbXwF1/dK4HwzKxgz38JfnK/jz9eVWV534Ifensb38v2tBfkV/hHCS/iL7cSk26PAN/E92bOS/8J/miDfhtHU9onx+ACp8DWPp3DD6f6icOPwl2fB/QHceGguXPF94bu6M7I4J+HL0U/g7eAcfDZpFG5YTsXLanwJ+a/G9yPOwNv80SXCgBv426X6v8bM/oPvBf03Pju2TrpuitvwGfp/4zNmTRppZvYa/pz8j6Qj04zzXng/MwWfsT0XN2LAD9a9gffJRzWhS56H4e3wQklzcQPp+hbIdQi+p3NGyvO7JLtA/vWKi1PYWfgh3kNxQ+tFfF8++CDzOXxf6ZP4VpHigXizyLeWDKRx+X8BeDHpdAH+2bn5mc5X4Xv0/5LF6Y3PNs/G+7c38D4XFq//T/BZ3p1SuOn4M9TS7XbNcSn+jDyLrzosXAExsxn4Ia9b8fZ+EI0nAf6L9wlvJHlXwQe+x0t6D7iMpgdNyyyFDcZBUJdI2hg3Bra3aMxVJc3qXGhmO9ValtYi6fv4FySa+u9aV+MHE1s7ixRUEEnXAM+YWZP/4CRo30g6Cv9Cy1HNBl4U5wTcmB1aOcmC9sby9/2vYJnCzJ7DPw8SVJ8F+GxN3VJmP18QBNVnBr7S2CLSNphv4jOxwXJEGK5BELQKM2vJFpFlgRvxrQBBEFQIM7u9+VBO2g86Bl9ab/P/VBe0b2KrQBAEQRAEQVAXxOGsIAiCIAiCoC4IwzUIgiAIgiCoC8JwDYIgCIIgCOqCMFyDIAiCIAiCuiAM1yAIgiAIgqAu+H8TP1r/vCdmaAAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<Figure size 720x432 with 1 Axes>" | |
] | |
}, | |
"metadata": { | |
"needs_background": "light" | |
}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"# mean bike age\n", | |
"ax = bike_share_bikes['Bike Life'].resample('M', how= lambda x: x.median()).plot(figsize=(10,6), fontsize=14, legend=False) \n", | |
"# ax = bike_share_bikes['Bike Life'].resample('M', how= lambda x: x.quantile(.75)).plot(figsize=(6,4), fontsize=14, label='25th Percentile', color='r')\n", | |
"\n", | |
"# ax.legend(fontsize=14, loc='best').set_zorder(20)\n", | |
"ax.set_ylabel('Years', fontsize=14)\n", | |
"ax.set_xlabel('')\n", | |
"\n", | |
"\n", | |
"ax.text(0.08, -0.15, \"From: chart-it (MikeRAzar.com/chart-it) | Source: www.capitalbikeshare.com/system-data\", \n", | |
" fontsize=12, transform=ax.transAxes) \n", | |
"\n", | |
"plt.savefig('/Users/mazar/Desktop/personal/chart-it/posts/Bikeshare/pt 4/bike_age_years.png', dpi=200, orientation='landscape',\n", | |
" bbox_inches='tight', pad_inches=0.1)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 11, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_stations_rides = bike_share_stations.groupby(['Starting Station', lambda x: x.month, lambda x: x.year]).\\\n", | |
"agg({'Minutes': 'count'})" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 12, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_bikes_mins = bike_share_bikes.groupby(['Bike', lambda x: x.month, lambda x: x.year]).\\\n", | |
"agg({'Minutes': 'sum'})" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 13, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_stations_rides.reset_index(inplace=True)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 14, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_bikes_mins.reset_index(inplace=True)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 15, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_bikes_mins['Date'] = bike_share_bikes_mins.apply(lambda x: datetime(int(x['level_2']), int(x['level_1']), 15), axis=1)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 16, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_stations_rides['Date'] = bike_share_stations_rides.apply(lambda x: datetime(int(x['level_2']), int(x['level_1']), 15), axis=1)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 17, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_stations_rides.set_index(['Starting Station', 'Date'], inplace=True)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 18, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_bikes_mins.set_index(['Bike', 'Date'], inplace=True)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 19, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"name": "stderr", | |
"output_type": "stream", | |
"text": [ | |
"C:\\Users\\mazar\\Anaconda3\\lib\\site-packages\\ipykernel_launcher.py:1: FutureWarning: sortlevel is deprecated, use sort_index(level= ...)\n", | |
" \"\"\"Entry point for launching an IPython kernel.\n" | |
] | |
} | |
], | |
"source": [ | |
"bike_share_stations_rides.sortlevel(['Starting Station','Date'], inplace=True)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 20, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"name": "stderr", | |
"output_type": "stream", | |
"text": [ | |
"C:\\Users\\mazar\\Anaconda3\\lib\\site-packages\\ipykernel_launcher.py:1: FutureWarning: sortlevel is deprecated, use sort_index(level= ...)\n", | |
" \"\"\"Entry point for launching an IPython kernel.\n" | |
] | |
} | |
], | |
"source": [ | |
"bike_share_bikes_mins.sortlevel(['Bike','Date'], inplace=True)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 21, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_bikes_mins = bike_share_bikes_mins.groupby(level=0).cumsum()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 22, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"bike_share_stations_rides = bike_share_stations_rides.groupby(level=0).cumsum()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 23, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAr4AAAFoCAYAAACxCPjaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xl4lNXZx/Hvmclk3xP2RQQBEVHcAQVXFJfWWq3aWpe6tHWpWrVaq63aWlvX2mq1arW2tr7VurTugAiibIoKAgJG2ffseyaznPePCWSeyUICk8xM8vtc11xkznPmyTljTO6c3Oc+xlqLiIiIiEhP54r1AEREREREuoMCXxERERHpFRT4ioiIiEivoMBXRERERHoFBb4iIiIi0iso8BURERGRXkGBr4iIiIj0Ct0a+BpjphhjXjPGbDbGWGPMJbvpP6ypX+RjWkS/Y40xnxhjGowxa4wxP27lXlcZY9Y29fnEGDM5ytMTERERkTjW3Su+mcBy4DqgvhOvmwYMCHu8t/OCMWZf4C1gPnAI8DvgEWPM2WF9zgP+CNzT1Gc+8LYxZujeTEZEREREEoeJ1cltxpga4Bpr7bPt9BkGrAWOsNYubqPPvcC3rbUjw9r+Coy11k5ser4I+Nxae0VYnyLgJWvtrXs/GxERERGJd0mxHkAHvWKMSQWKgD9Ya18KuzYRmBHRfzpwsTHGAxjgMOCBiD4zgEmRn6iyslJnOIuIiIgkuJycHBPZFu+b22qAm4BzgdOAWcALxpjvh/XpD2yPeN12QkF9YdPD3Uaf/l0wZhERERGJQ3G94mutLQEeDGtabIwpBG4G/hneNeKlJqzdtNNHq7siIiIivUS8r/i2ZhEwMuz5Nlqu3PYF/EApUAIE2ugTuQqc8IqKimI9hJjS/DX/3qo3zx00f81f8+/NOjP/RAx8xwNbw54vAE6K6DMVWGyt9VlrG4FPmtoi+8zvslGKiIiISFzp1lQHY0wmsF/TUxcw1BgzHiiz1m5opf/FgA/4DAgC3wCuBm4J6/YX4BpjzMPAE8DRwCXAd8P6PAQ8Z4z5CJgH/BgY2PRaEREREekFujvH93Bgdtjzu5oefwcuMcbcCdxhrQ3fhXc7sA+hdIUvgUuttbvye621a40xpwF/AK4EtgDXWmtfDuvzgjGmoOleAwjVEj7NWrs++lMUERERkXjUrYGvtXYOzZvNWrMvMDOs/98JBcW7u+/7wKG76fMY8FiHBioiIiIiPU7cVHUwxhjgBODEWI9FRERERHqeuAl8begIuSGxHoeIiIiI9EyJWNVBRERERKTTFPiKiIiISK+gwFdEREREegUFviIiIiLSKyjwFREREZFeQYGviIiIiPQKCnxFREREuoAvaFmw3cu2ukCshyJN4qaOr4iIiEhPsbS0kUtml7G2OkBGkuFvx+Vz8pDUWA+r19OKr4iIiEgU/ePLWk5+s5i11aGV3lq/5ZI5ZSwpaYzxyESBr4iIiEgU1PstV39YzrXzKvBGZDfU+S3nvVvK+mp/bAYngAJfERERkb22tsrPyW8W86+iujb7bK8Pcu7MUiq8wW4cWWwEgpZ1cRjkK/AVERER2QuLixs59vUdLCvzOdrTkwzTIvJ6V1f6ueC9UrwB251D7Fbrqv2c8U4Jp75VTHmcBfkKfEVERET2UFlDgIvfK6Oq0RnIjsxJYtYZfXj+xHzOGpbmuDZvWyNXf1hO0Pas4Ndayz++rOWY/+5gwfZGttYFuWlBRayH5aDAV0RERGQPWGu5Zl4FmyPKlZ05LJX3vtGHMXkeXMbw+OQ8JvZLdvR5aU09d39a1Z3D7VLF9QG+N6uMa+dVUONvDuhfXlvPK2vaTv/obgp8RURERPbAUytreWtDg6Ptxwdk8Oxx+WR5mkOs1CTD8ycWMDLHWUX2oc9reOHr+AkK98T2ugD3L6liwqs7eHtjQ4vrUwelMKl/SgxG1jrV8RURERHppGVlPn65uNLRdmihh18fnoMxpkX/vBQX/5lawNQ3iiluaM57vXZeOcOzkjiib3KL18Qray0fbGvkmVW1vLG+Hn8rGRvpSYa7j8jhB6PTW30/YkWBr4iIiEgn1PqCXDqnzFGyLMtjeOa4fJLdbQd5w7KS+L+TCjj97eJdr/UG4IL3SnnvjD4MzoyvsMxay+pKP+urA2yo8bOhJvTvslIfa6rbPo3u8D4enpicz4ic+JoPKPAVERER6ZSbF1VSVOks1fXwpFyGZe0+rDq8TzKPHJ3HD+eW72rbUR/ke7PKePu0QjI88ZGFWtoQ4Btvl/BFRcdLkqW5DTcenMX14zJJcsXPKm+4+Hh3RURERBLAi1/XtajV+/2R6Zw9PL3D9zh3RDo3HJTpaPu8zMdVcVTp4aHPazoc9I7KSeL3R+Ww8rz+3HRwVtwGvaAVXxEREZEOmbW5gWs+LHe0jcpJ4t6jcjp9r9sPzWZlud+xIex/6xq4d0k1tx6Svddj3VvTW9moFi7JwDf2SePS/TM4pn9yXOXxtkeBr4iIiMhuzNvm5fuzymgMO48hxQ1PH5e/R+kJLmN48tg8Tnmj2LGyeu+Sak4enMphfWK32W19tZ+vqprHlGTg/P3SGZrpZkhmEkMz3RyQ5yEvJfESBxT4ioiIiLRjcXEj580spT7itLWHJuYyLt+zx/fN8rh4/qQCTny9mNKwE87u/rSKV08p3OP77q33Nnsdz4/sm8yjx+TFaDTRlXihuoiIiEg3+by0kbNnlDgOZQC496gcLhiZsdf3H5aVxBNTnEHl7C1ePtjqbeMVXW/WZmeaw0mDU9vomXgU+IqIiIi0YnWFj7Oml1IZcRzxnYdl86MDMtt4VeedOCiFY/o7Uxvu/rQKG4ONbr6gZW5E0H3CwPg5gGJvKfAVERERibC8zMeZ75Q4UhAAfnZwFtcflBXVz2WM4ZeHOje0LdrRyMxN3b/qu7i4kSpfc8BdmOrioII9T+eINwp8RURERMLM3erltLeK2VbvDHqvGZvJLw6JbtC701H9Ujh5sHNl9e5Pq7q9vNmsTS1Xe10JUrGhIxT4ioiIiDR5eU0dZ88ocax6Aly2fwa/OSK7S8t23Rax6vt5mY/X17dfVqwtNb4gK8p8zNjYwOoKX4dfN2uL8/OdMKjn5PeCqjqIiIiIAPDI8mp++XFVi/Yfjcngd0fldHmt2oMLkjlzWCr/W9ccfN7zaRVnDE3FvZtDIV5fX88b6+tZVx1gbbWfHWGr1Qa4f0IOl49pPy+5tCHAkhJnkHzCoJ6T3wta8RUREZFerqQhwM8WVrQa9P768Gx+f1ROt/25/xeHZBMe466u9PPimvp2X/PmdjcXvlfGC1/Xs2hHoyPoBbDArR9Vsqys/ZXf2Vu8hK9zj8v30DfN3ckZxDcFviIiItLr+IOWdzbWc+F7pez/7208tbLWcd3jgqem5HHtuKxuPZVsdK6Hc4enOdp+/1kVjYHWc3231AZ4YM3uD7vwBeHHc8vwtnEfgFkR9XtP7GGrvaDAV0RERHqJal+Q2Zsb+NXHlYx9cRvnv1vG6+sbiCjRS5bH8NLUAr4zIj0m4/z5IdkkhcXa62sCPLGypkU/ay0/nV9OTaBlYJ5kYGimc7V2Rbmf+5a0XNXeea/3Nvfs/F5Qjq+IiIj0UFvrAiza3sjCHV4WbG9kWZmP4G6KJPRPc/Hi1AIOKojdkcHDspK4eHQGT69qXoW+Y3EVI3OSmDakeTX431/XMz2iCsMt47P47n7pDM5wk+Qy/PD9MkeqxB+W1TBtSBpH9HXOb3m5n+1hKRIZSYYJfWP3HnQVBb4iIiKSMJaX+Xh1bR2bagMUprrpl+bCVrk5OKOBjCQXS0sbWbQj9NhQE+jwfUfmJHHBfulcMjqD3JTY/0H8poOzeOGrul0nxgUtXDqnnNenuTmsTzLb6gL8fFGF4zUT+yVzy/gsRz7yfRNy+WCbl611wV33ufKDcuae2Yf0pOZ5Rq72Th6QQrK755Qx20mBr4iIiMS10oYAL62p519FdXze6gatFPiytNP3zUwyfHt4Ghfsl86RfZO7NZd3dwaku3liSh4Xzi7btUpd57ec924pM0/vw20fVzpOlEt1w6NH57XYhJeb4uKRo/M4Z2bz+/NVlZ+7Fldx74TcXW29Ib8XujnH1xgzxRjzmjFmszHGGmMu2U3/44wx/zPGbDXG1BljPjfGXNpKH9vKY/+IfmcbY74wxnib/j2rC6YoIiIiUVJU6ePi2aXs/8I2bllU2UbQ2zmjcpK4aFQ6T03JY/X5/fnT0Xkc1S8lroLenU7fJ437J+Q42koagpz4RjFvbXCu0N52aDYjclpfzzxpcCqXjHLmKz+xspZLZpfx+vp6yr1BFm6PDHx7Xn4vdP+KbyawHPhH02N3JgHLgPuArcApwJPGmAZr7fMRfccCZWHPi3d+YIyZCLwA3AG8Anwb+I8x5mhr7aI9nIuIiIh0kfe3NHDhe2UtDpLoDI8Lxhd4OKpvChP7JXNUv2QKUxOrPNdl+2eyqSbAH5Y1b24rizhGeVxWgKsOaL9G72+OzGH2Fi/rw9I//ruunv+uqyfFDY1htxyW5WZ4ds9MCujWWVlr3wLeAjDGPNuB/vdEND1ujDkeOBuIDHx3WGtL2rjV9cBsa+1vm57/tuk+1wPf7eDwRUREpBs8X1TLtfMqWlRbgNBhDCcMSmHakFRqfJZtdQG+Lq6k1p1OpTfI0KwkJvQNBbmHFCSTmhR/K7md9avDstlcF+DFr1vW801xwy9HNu72gIssj4vHJudxxtslRL6t3ohU6J662guJmeObDWxqpX2xMSYF+AK421o7O+zaROCRiP7TgWu6ZogiIiLSWdZafr+kmnuXVLe4NiLbzQUjMzhvRDqDMpyrtkVFxYwcuU93DbPbGWN49Og8ttcFeX+rMyXh1vHZ7JtW16H7HN0/hX+fVMBNCyvY2M7GvxMG9sz8XgBj7Z7/CWGvPrExNcA11tpnO/GaM4BXgaOttR81tY0Gjgc+BpKBC4EfA8dZa+c29WkELrfW/iPsXhcBT1lrHf91Kysrd70hRUVFezY5ERER6RRfEH77VTJv7mi5JnftsEa+P8hPHKbhdqsaP/xwWSpFtaEtWgdnB/jLOC+dXdS2Fj6vdjGj2M27JUmU+ZpvUJgc5JXDGkjUA9tGjhy56+OcnJwW70zCrPgaY44mlN5w7c6gF8BauxpYHdZ1gTFmGHATMDesPTLCN620OYS/eYmiqKgoIccdLZq/5t9b59+b5w6af6LPv8Ib5KLZZczd4VzNTHHDE5Pz+da+aW28MiTR598Zc/YL8szqWrBw2ZgM0pNcezT/UcA5hE6w+2Crlzc2NFDtC3LN2MyY1jDeE52Zf0IEvsaYYwjlBv/KWvt4B16yCDg/7Pk2oH9En77A9uiMUERERPbEmio/571bSlGl39Gen+Li+RPzmdCv5/7ZfU9keFz85MCsqN0vyWU4flAqx/fgvN5wsa/QvBvGmCnA28Bd1tqHO/iy8YSqQOy0AJga0WcqMH/vRygiIiJ74sNtXk58Y0eLoHffLDczT++joFeirltXfI0xmcB+TU9dwFBjzHigzFq7oZX+xwFvAo8B/zLG7Fy1DVhri5v6XA+sA1YQyvH9PvAtQpUfdvojMNcYcyuhHOGzCOUFHxPN+YmIiEjHPPdlLTcsqMDnrMzFkX2Sef6k/IQrOyaJobtXfA8HPmt6pAF3NX38awBjzJ3GmPC820uAdEL5ulvDHh+H9UkGHgA+Bz4gFMyebq19ZWcHa+18QqkPFzf1uwg4TzV8RUREulfQWn71cSU/mdcy6D13eBqvTStU0Ctdprvr+M4htKmsLfsCM8P6X0Io+G3vnvcROuBid5/7JeClDgxTREREusidi6v40/KaFu23H5rNjQdlxuUJatJzxM3mNhP6Sj8BODHWYxEREZHoe2N9fYugN81t+MuUPM4c1n7lBpFoiJvA14YKCg+J9ThEREQk+tZV+7nqw3JHW980Fy+cVMAhhYlVPksSV9xXdRAREZHE1uC3XDy7jKrG5m08SQb+eUK+gl7pVgp8RUREpEv94qNKlpb6HG2/PiKHI/uqXJl0LwW+IiIi0mX+83Vd6KSxMN/cJ5UrD8iI0YikN1PgKyIiIl1idYWP6+dXONqGZ7l55Jg8VW+QmFDgKyIiIlFX2RjkovfKqPU35/WmuOHZ4/PJSVb4IbGhrzwRERGJKn/QcumcMlZHHEV8/4RcDirQZjaJHQW+IiIiElW3fVTJrM1eR9v5I9K4cGR6jEYkEqLAV0RERKLmmVW1PLHSuZnt8D4e/jBJeb0Sewp8RUREJCrmbGngZwudm9kGZ7j51wkFpCUp6JXYU+ArIiIie62o0sfFs8sINO9lIyPJ8PyJ+fRLd8duYCJhFPiKiIjIXilrCHDezFIqw05mM8CTU/K0mU3iigJfERER2WMNfsv3ZpWxpjrgaL/z8GxO3yctRqMSaZ0CXxEREdkjQWu58oNyFu5odLR/d790rj0wM0ajEmmbAl8RERHZI3ctruLVdfWOtkn9knl4Uq4qOEhcUuArIiIinfb0qhr+uLzG0TYqJ4nnTywgxa2gV+KTAl8RERHplHc21vOzhZWOtj6pLl6cWkBuikILiV/66hQREZEO+6ykkUvnlBMMK1uWnmR4cWoBw7KSYjcwkQ5Q4CsiIiIdsrS0kbOml1Dnb456XQaePjaPQwpVtkzinwJfERER2a3PSxv51vQSKsJq9QLce1QOpw5V2TJJDAp8RUREhLKGALM3N7CtLtDi2vIyH9+aXkq51xn03nBQJleMUdkySRxKxhEREenFyhoC/HFZDU+urKU+YDHA5AEpnDM8jW/uk8bm2gBnvlNCmTfoeN21B2byy0OzYzNokT2kwFdERKQXqvYFeXxFDY8ur6HK17ySa4G5W73M3erlxgUVpLoN1T7nSu/VYzO56/Bs1eqVhKPAV0REpBep8wf52+o6HlpaTWnEKm4kXxB8QWfQe+UBGdx9hIJeSUwKfEVERHqBCm+Qp1fV8vgXNZQ0tB7wZiQZav221WsAPxqTwT1H5ijolYSlwFdERKQH21Ef4PEVNTy9qtaR0hAuL8Xw03FZXD4mg821AV5aU8/La+r5qsq/q88VYzL4/VEKeiWxKfAVERHpgZaWNvLkylpeWlOHt2WhBgAykwxXHZjJ1WMzyUkOFXoamePi1kM8/Hx8FktLfSza0ch+OUmcMDBFQa8kPAW+IiIiPYQvaHl9XT1Prqxl4Y7GNvtleQyX7Z/BNQdmUpjqbrWPMYbxhcmM18EU0oMo8BUREUlwFT54cGk1f11Vw9a6tjesFaS4uGpsJpftn0Fuikr5S++jwFdERCRBrany89iKGv75ZRoNwao2+w3JdHP12EwuGpVOepICXum9FPiKiIgkmEXbvfxpeQ1vbWggtF2t9dzbYwek8MMxGUwbkorbpfxcEQW+IiIiCeLjHY3c81kVs7d42+yTnmQ4f0Q6V4zJYEyepxtHJxL/FPiKiIjEuSUloYB3xqa2A94B6S5+NCaTS0Yrf1ekLQp8RURE4pC1lnnbG3lsRSiloS3j8j2cU1jDlRP2JdmtdAaR9ijwFRER6YRA0PK31bUsL/Nx/n7pTOiXEtX7VzYGeeGrOp5ZXcuqCn+b/Y7un8zNB2cxZUAKX31VqaBXpAO6NfA1xkwBbgIOAwYCP7DWPrub14wDHgWOBMqAJ4DfWGttWJ+zgd8AI4Cvgdusta+GXTfAHcAPgTxgEXC1tXZF1CYnIiI9nrWW6+dX8FxRHQD/KKrj4Um5XDQqo9X+Fd4gL62po6gyFMDuPP/BNH3swuA24HaByxi21wV4ZW19u8cGH9U3mV8cks2UAck6UEKkk7p7xTcTWA78o+nRLmNMNjATmAscAYwGngVqgQeb+kwEXiAU2L4CfBv4jzHmaGvtoqZb3QzcCFwCrAZ+Bcw0xoy21lZHaW4iItLDPfh5za6gFyBo4dp5FZQ2BLl+XKYjEJ25qYGffFjOtvq26+p2xqGFHm47NFsnqInshW4NfK21bwFvARhjnu3ASy4A0oGLrbX1wHJjzBjgBmPMQ02rvtcDs621v216zW+NMcc3tX+3abX3euD31tqXmz73xcAO4HuEVpBFRETa9cLXddz9aeu1cu/6pIqShiC/OSKbOr/l9o8qefbLulb7dobbwOlDU7ls/wymDFDAK7K34j3HdyLwQVPQu9N0QmkNw4C1TX0eiXjddOCapo/3BfoDM3ZetNbWG2PmApNQ4CsiIrsxd6uXaz4sb7fPn1fUsKHGz7IyH+uqA3v1+Qaku7h4VAYXjcpgYEbrRwqLSOeZsFTZ7v3ExtQA17SX42uMmQFsstZeGtY2FFgPTLLWLjDGNAKXW2v/EdbnIuApa22KMWYSMA/Yx1q7IazPM8Aga+0p4Z+zsrJy1xtSVFS0t9MUEZE4U+uHv23y8HWti1S3Jc0F6W5Lujv076BUy+jMIINTLS4Da+oMly1NpSbQvNrqMZYr9/Hx9EYPtYH2V2FdWL4zwM/AVOfP26AFCwQsBJv+BRiVYTk6L4AOWBPpvJEjR+76OCcnp8X/nPG+4gsQGZmbVtpb6xPZ1pE+DuFvXqIoKipKyHFHi+av+ffW+ffmuUPH51/SEOCKGaUsKfXttm9mkuHAfA8bavzUBJx5uo9Pyeec4emcXdrIOTNKKW5oPY93v+wk/jIlj8P7JHdsIntI//01f82/Y/OP998ntxFKUwjXt+nf7bvpE36d3fQREZEebmONn1PfKulQ0AtQ47cs3NHIljpnUHvHYdmcMzwdgIMLknnntD4MzWyZjvDDMRnMPbNPlwe9ItJx8R74LgAmG2NSw9qmAluAdWF9pka8biowv+njtYSC3119mu43OayPiIj0YKsrfEx7s2RXWbE99YPR6Vw/LtPRNiIniemn92HKgFA935E5SfzvlALum5BLuvIVROJKd9fxzQT2a3rqAoYaY8YDZeH5t2GeJ1Sm7FljzN3AKODnwF1hdXz/CMw1xtwKvAqcBRwPHANgrbXGmIeB24wxq4AvgduBmqb7i4hID7CkpJE3NzRQ6w8yIN3NwHQ3AzPc1Pstl79fTpnXuXJ73MAULhqZTo3fUuOz1PqClHmDrCj383lpIxWNzmy4UwancP+E3FYrKwxId/PatELKvUHydFywSNzq7hzfw4HZYc/vanr8HbjEGHMncIe11gBYayuNMVOBPwOLgXJC9Xsf2nkDa+18Y8z5wN1N9/oaOC+shi/AfUBa0312HmBxsmr4iogkvnJvkHu+8vDfbcXtb9wI861haTwxJY+UNk47s9aysTbA56U+Vlf46Zfu4vwR6SS52t/IpqBXJL51dx3fOTRvTmvNvoQOrAh/zTJgym7u+xLwUjvXLXBn00NERHoAay3//rqeX35cSUmDp8Ov+8HodB6YkIu7nSDWGMPQzCSGZiZxxj7RGK2IxIO4qerQdNDECcCJsR6LiIjEt1UVPm5cUMG8bY2det2NB2Vy+6HZOghCpJeKm8C3aVV2SKzHISIi8WtHfYD7llTzt9W1u+rehtsvO4lzR6RRXB9kS12ArXUBttQGSHIZbjwoix/sn9H9gxaRuBE3ga+IiEhbanxBHl1ew6PLa6jxt4x4k43lpvHZXDcuq828XRERBb4iIhK3/EHLc1/W8bslVeyob/2QiJMGpXB1/wqOP2hwN49ORBKNAl8REYlLS0oauW5+BUvbOHBin0w3vz4ih2/uk8pXX5V38+hEJBEp8BURkbhS4wvyu8+qefyLGoKt5PHmp7i4eXwWl47OIFlpDSLSCQp8RUQkbkzf2MCNCyrYVBtocS3NbbhqbAbXjssiJ1n1ckWk8xT4iohIzNX6gvx0QQUvfl3f6vWzhqVx95E5DMpwd/PIRKQnUeArIiIxta7azwWzSllR7m9xbUimmwcn5HLykNQYjExEehoFviIiEjOzNzfwgzllVDQ6k3ldBq46IJNbD8kiw6O0BhGJDgW+IiLS7ay1/Gl5DXd9UtViA9sBuUk8NjmP8YXJsRmciPRYCnxFRKRbeQOWqz4o5+W1LfN5zxqWxqPH5GqVV0S6hAJfERHpNtZafvJhy6DXAHccls114zIxRiXKRKRrKPAVEZFuc9/Sal5c4wx6c5INzxyXz4mDtIFNRLqWAl8REekW//m6jt99Vu1oG5Ht5qWpheybrR9HItL1lEQlIiJdbuF2L1d/6DxWOD/FpaBXRLqVAl8REelSa6v8XDCrjMZgc1uyC/51Yr6CXhHpVgp8RUSky1R4g5z3biml3qCj/dFj8pjYLyVGoxKR3kqBr4iIdImgtVz+fhlfVjpPZLtlfBbnjkiP0ahEpDdT4CsiIl3iyZW1vLvZ62g7Z3gaPx+fFaMRiUhvp8BXRESibmW5jzsWVzrajuyTzKNH56lOr4jEjAJfERGJKm/AcsXccryB5rbsZMPTx+WRmqSgV0RiR4GviIhE1W8/rWJ5mc/R9uCEXIZkqoKDiMSWAl8REYmauVu9PLK8xtF2zvA0vqPNbCISBxT4iohIVFR4g1w5txwb1jY4w80DE3JjNiYRkXAKfEVEJCpuWljB5rrmxF4DPDY5j9wU/agRkfig70YiIrLXXl5Tx0tr6h1tPzkwkykDdEiFiMQPBb4iIrJXttcFuGlhhaPtwHwPtx2aHaMRiYi0ToGviIjsMWst182voNzbnNmb7IKnpuSR4lbpMhGJL50KfI0xLmOMK+x5f2PM5caYo6M/NBERiXf/91Ud72xscLTddmg2Y/I8MRqRiEjbOrvi+ybwEwBjTCawGLgfmGOMuSjKYxMRkTi2uTbAzz9qeTrbNWMzYzQiEZH2dTbwPQx4r+njbwNVQF/gCuCmKI5LRETimLWWa+eVU9XYnOKQ5jY8NjkXt0spDiISnzob+GYBO3cwnAy8aq31EQqGR0RzYCIiEr/+8WUdszZ7HW2/Oiyb/XKU4iAi8auzge8G4GhjTAZwCjCzqT0fqIvmwEREJP5Ya3lnYz23RaQ4TOqXzI8OyIjRqEREOqazB6c/BDwH1ADrgblN7VOAZVEcl4iIxBFrLbM2e7nnsyo+LfE5rmWB/6gWAAAgAElEQVQkGR6bnIfLKMVBROJbpwJfa+0TxpjFwFBgprU22HTpa+CX0R6ciIjElrWW97d6uefTaj4qbmy1z6+PyGZYVmfXUUREul+HUx2MMR5jzCKgxlr7qrW2Zuc1a+2b1tp5XTFAY8w6Y4xt5fFm0/U7W7m2LeIepqnfFmNMvTFmjjFmbFeMV0Skp7DW8vNFlXxremmbQe+VB2Rw6WilOIhIYujwr+jWWp8xZl/A7rZzdB0BuMOeDwA+AV4Ma1sNHBf2PIDTzcCNwCVNfX8FzDTGjLbWVkd5vCIiUbNgu5cb51ewtjrA+EIPxw1M4bgBKRzWJ5mkLq6e8PxXdTyxsrbVa6cPTeXnh2QzLl+b2UQkcXT2b1N/J1S67GddMJZWWWuLw58bYy4jVEbtP2HNfmutY5U3rL8Brgd+b619uantYmAH8D3gia4Yt4jI3rDW8pcvavnlx5X4m5YbFmxvZMH2Rn73WTXZHsPR/VO4cFQ6pw1Ni/rnL6r0cfPCyhbtpwxJ5dbxWYwvTI765xQR6WqdDXwzgAuMMVMJrbo6lgKstddGa2CtaQpiLwP+aa0NryIx3BizGWgEFgG/sNauabq2L9AfmBE2znpjzFxgEgp8RSTO1PiCXDuvglfW1rfZp8pneXtjA29vbOCc4Wlc1Td6n98bsFw2p5xaf/Mf+NKTDC+fXMDEfinR+0QiIt3MWNvxzAVjzOx2Lltr7Ql7P6R2P//JwHTgEGvtkqa2UwnVF15F6DCN24H9gbHW2lJjzCRgHrCPtXZD2L2eAQZZa08J/xyVlZW73pCioqKunI6ISAvr6gw3r0phbV3nqk32Swly16hGDssJ7r7zbvxhjYfntzhTGG7fz8uZ/SOzyERE4svIkSN3fZyTk9MiH6yzVR2Oj8KY9sYVwMc7g14Aa+3b4R2MMQuBNcDFhMqv7eoacS/TSptD+JuXKIqKihJy3NGi+Wv+iTb/0oYAK8r9fFHu44tyH6+urafa5/zWlO0x3D8xl8aAZc4WL+9v9VLS4Axwt3tdXLkslevGZfKLQ7JJdu9Z/u+MjQ08v6XU0XbWsDRuPGYgJo7LlSXif/to0vw1f82/Y/NPmPozxpi+wJnA1e31s9bWGGNWADvfgZ25v/2BjWFd+wLboz1OEZHdCVrL75dU89yXtWyta3+F9oC8JJ47voAROaFv1xeOyiBoLZ8U+7h+fjkryv27+lrg4WU1vLvZy6WjM5g2JJWBGe427tzStroAV35Q7mgbkunmD5Ny4zroFRHpqE4FvsaY19q7bq395t4Np10/ALzAv9vrZIxJJZTqsDMtYy2h4Hcq8HFYn8l04yY9EZGdHlxazX1Ldl9Q5tzhafxhUi4ZHmfag8sYjuibzHvf6MtvPqni0RU1juvLy3zcsKCCGxbA+AIP04akMmVAClW+IFtqg2ypDbCp1k9JQ5Bqn6WqMfRvmTdIXVher9vAX6fkkZvS2UM+RUTiU2dXfEsjnnuAg4EhwCtRGVErmja1XQ78O7L8mDHmAeB1Qscp9yV0kEYGoQoUWGutMeZh4DZjzCrgS0J5wDXA8101ZhGR1sza3MA9n7Uf9GYkGe44LJsrxmS0u9Ka4jbcfWQOUwencMXsYnY0tgxQl5T6WFLq4/cdCLQj3XpINkdpM5uI9CCdzfH9QWvtxpgHga6sh3scsB9wQSvXBgP/BxQCxcBCYIK1dn1Yn/uANODPQB6hyg8nq4aviHSn9dV+Ln+/zLG5IMnA/nkeDshLYmyehwPyPBzRJ7lTq6zHDkzl+UMaeGxHQbuVIDrjmP7J/HRcZlTuJSISL6KV4/sE8CFwZ5Tu52CtnU1oM1pr187vwOstobHdGdWBiYh0UIPfctHsMsq9zWGvy8DLJxdw7MDUvb5/jgeeOS6fm8f7eGtDA29vqGdxsa/TJw4Z4Oj+yTxzXD7uLj4gQ0Sku0Ur8B0dpfuIiPRIP1tYwdJSn6Ptl4dmRyXoDbd/rof9cz3ccFAWO+oDTN/YwDsbG1hfE6Aw1cWgDDcD090MynDTL81FTrKL7GQXWR6z69+uPhFORCRWOru57U+RTYSOED4VeCZagxIR6Un+8WUtzxXVOdpOH5rK9V2cStA3zc2FozK4cFRGl34eEZFE0dkV33ERz4OE8mp/igJfEZEWPi1u5KYFFY62EdluHpucpxJhIiLdLNEOsBARSRhFlT7OfbeUxrBSvelJhudOKCAnWSXCRES62x7l+DbVwd2PUL30r621DVEdlYhIgttU4+es6aUtTlj709G5HJDnaeNVIiLSlTq15GCM8Rhj7gfKgaXAMqDcGHOfMUbfyUVEgJKGAGfNKGVTbcDRft2BmZwzPD1GoxIRkc6u+N4LfBf4MaHyZRA6Ae13hILom6I3NBGRxFPZGOTsGaUUVfod7ReOTOfOw7NjNCoREYHOB77fAy611r4V1va1MaYY+CsKfEWkF6v3W777bmmLsmVnDkvl4Um52swmIhJjnd1dkQN83Ur710Du3g9HRCQxWWu57P0y5m9vdLQfPzCFJ6foMAgRkXjQ2cB3KXBtK+3XAUv2fjgiIonpzQ0NvLXBuc/3iD4enjshnxS3gl4RkXjQ2VSHm4G3jDFTgQWEqjpMBAYSOsRCRKTXsdby4OfVjrYD8pL4z9RCMj0qWyYiEi869R3ZWjsXGAX8B8gEsps+Hm2t/bC914qI9FSzt3j5rMSZ1/vUlHxyUxT0iojEkw6t+BpjhkY0PRF5H2PMUGvthugMS0QkcTyw1Lnae9rQVMbmq8KjiEi86WiqwzpCaQ1tMU3X3Xs7IBHpmO11AaZvaqDSGyTFbUhNMlQWu1mVXE9hqosj+iSTpA1VXW7Bdm+LDW03HZQVo9GIiEh7Ohr4HhH2sQHeJ1TabFPURyQi7arwBvnT8mr+8kUtdf7I30dT4MsyACb0TeaFqYl1NO6nxY3csbiSJJfhN0fkcGACrJo+GLHae/zAFA7tkxyj0YiISHs6FPhaaz8Jf26MCQLLrLVrumRUItJCnT/IE1/U8vCyaiob2/sDTMjCHY2cO7OUl04uICsBNlgtL/Nx5vQSqn2huf1obhkfntk3rmvfLilp5N3NXkfbjQdrtVdEJF51tqqDiHSjqsYgHxc3smBbI/8sqmVbfbBTr1+0o5Hz3y3lP1MLSE+K3+B3a12A82aW7gp6AVaU+1lXHWDf7Pj9NhVZyWFC32SO7qfVXhGReBW/P1FEeiFrLe9u9vLupgYWbG9kebmPYDuLu33TXJw5LI2gDZ0aVlJRxXpfKqvDjsudt62R780q498nFpCaFH+rp9W+IOfOLGVzXaDFtblbvXEb+K4s9/H6emfd3hsPzorrFWoRkd5ub36i7P5vrSLSYRtq/Fw7r4I5W7y77ZudbLjuwCx+fEAGGWFpDEVFJfTdZwjfml7iKK81Z4uXC98r5Z8nFuw6TCEQtDQGIS2GwbA/aLlsThnLynytXp+71cvFozO6eVQd84eI1d6DCzycNCglRqMREZGO6Gg5s9cimlKBp4wxdeGN1tpvRmtgIr2FtZZ/fFnH7R9XOv7U35o0t+FHB2Rw3bgs8tqoEZuT7OKVkwv55jsljoBy5mYvR726HQOUe4NUNlosMC7fw0/HZXLmsLRuPVbXWsvPF1UyY1Pbgf4H27xYa+NuFXV1hY+X1tY72m44SKu9IiLxrqMrvqURz/8Z7YGI9GSVjUE+L/WR6jbkpRjyUlzkJLvYVhfg2nkVvNfOKu/onCQm9EtmYr8UThqcQmHq7qsG5qW4+O8pBZzxdgkrK5rTHtZVt0wnWFbm49L3yxm9pJqbDs7i2/t2TwD85xU1/HVVraPt4AIPX1f6qWmqVrGjPsjqSj/758ZPdYftdQHOe7fUkYIyOieJb+yTGrtBiYhIh3S0qsMPunogIj3VvG1eLnqvjFJvy41pLkOrObyT+yfzowMymdAvuUOBbmsKUt38b1ohp79dQlFYzm9bVlf6uWJuOfcuqeZn47M4d3hal61gLtru5VeLqxxtgzPcvHBSAdfOK3esAs/d4o2bwLeyMcg5M0tb/AJx08FZuLTaKyIS9+J3m7dID/DSmjrOml7SatALLYPe9CTDAxNy+N+0Qs7YJ22Pg96d+qa5eW1aIWNyO57O/1WVnx/NLefWjyr36nO3pdYX5MoPyh1zz/IYXjipgP7pbiYPcObJzt26+5zn7tDgt1wwq7RFPvJ5I9I4Z3hajEYlIiKdEZ/bpUUSnLWWPy6r4c5PqnbfucnR/ZP58zF5DMuK7v+WA9LdzPlmXz4racQA+aku8lJc5Ca7WFft54Gl1fxnTT2BiCD86VW13HZodtRrAN/5SRVrIlZMnz42f9cRv1MiAt8Pt3kJBG235h9HCgQtP5xbxofbnCe0nTw4hUePyVNur4hIglDgKxJl/qDl5oWVPLO6tsW1cfkeqn3BXZvLAAamu7h2XBY/HJPRZX8uT3EbJvRrWXFgvxwPf5mSz83jQwHwC1/X7QqAfUGYvdnLN4dFbzXz/S0NPLXS+b5ctn8GJw9pzo8dl+8hN9lQ0fT+VDRalpX5GF/Y8fq4/qClpCHIjvoAxQ1BdtQHqfMHGZfv4fA+yZ16n621/GxhJa9FlC47oo+Hvx2Xj0fHQouIJAwFviJRVOMLctmcMqZHVCrwuODRY/I4b0T6rrZA0FLts+Qkm5ivGA7PTuKxyXnkJBse/6I5MJ2+qSFqgW9lY5CrP6xwtA3LcnPX4dmONpcxTB6Q4qiR+8FW724DX2stM4vdXLxiOyvL/W3WWxyU7uYbw1L51rA0juy7+yD4wc9rWvwSMzoniRdOKnCUkhMRkfin79oiUdLgt5z/bmmLoDc72fDyyYWOoBfA7TLkprhiHvSGmzbEWZlg5qYGgjY6Jbtv+6iSTbXNKQ4GeOyYPDJbCR4n9+9cnu+mGj/nzyrjF6tT+KKdoBdgc12Av3xRy7S3Shj74jae+KIG28Yc393UwG8/daarDEp38/LJBeTvZf61iIh0PwW+IlEQCFquaCUHdHCGm+mn9WmRtxqvJvZLIcvTHIjvqA+ypKT1wyU64+0N9fyzyFH2m2sOzGRS/9bflykDne3ztzfia6X8RSBoeeKLGia8uoPpGxtaXN+drXVBbllUyW0fV7YIftdV+7n8/TJHEJ2bbHj5lAIGZ+qPZSIiiUiBr8hestZy08KKFsfXjsv3MPOMPozJi49SXB2R7DYcHxF0vrOp8wFluLKGANfNd6Y4jM5J4rZDstt4Reh637Tmb0+1fsunxc5fKlZV+Jj2VjG3LKrcVfc3XH6Ki/1zk5jcP5lzhqfxjX1SSW/jlLrHVtTy80XNwW+D33Lx7LJdecYQKj339+ML4qa0moiIdJ6WLUT20r1Lqvnbaudq5uicJF6bVtjm6Wrx7JQhqY6NXNM3NvCLdoLU9lhruX5+BTvqm8u5uQ38ZUoeqe0clWyMYcqAFF5a03w62tytXo5q2qC3qsLHyW8UU9XKSXdH9U3m4Um5rf7CUecPMnOTl/+tq+ftDQ3Uh5WyeGJlLQEL903I4WcLK1ha6lzp/tWh2Rw7MDFW7kVEpHWJ91NZJI48s6qW3y+pdrTtzAFNxKAXYOrgVMJD0qWlPrbWtTzxrSP+9VVdi2oINx6cxSEdqNAQmR6yM8+3qjHI92eVtQh6szyGm0c08vZphW2usqcnuThzWBrPHJfPjDP6kB/x3+ivq2o59a0SnotIyzhjaCrXjcvc7ZhFRCS+JeZPZpE48Nq6em5c4PwTfk/IAe2b5ubQQmfgOHMP0h3WVPm5ZaHzEIyDCzzcdFBWh14fGfh+VNxInT/IVR+U81WV8yS604amsvCsfnxngL/DpcrG5Xt4fVohhanOb4OLdjhTKkZku/nzZNXqFRHpCRT4iuyBZWU+fjjXufEpzW14cWrPyAE9JaK6wzud3DjmC1queL+M2rDc2/Qkw1+PzSPZ3bEAcp9MN0MymysneANw6Zxy3tjgHMv39kvnXyfkMyij81UWxuZ7eOPUQkc+cbj0JMNzJxSQk6xvlSIiPYG+m4t0UoU3yEXvldIQ9td/t4G/HZ/HkX17Rg5oZOA7Z4uXhlY2kLXl3iXVfBJRDeKeI3MYmdPxXwp25vmGiwzADy7w8ODE3L1ajd0/18Mb0wrp30rw+8dJuRyQQJsTRUSkfQp8RTrBWstVH5azNuLI3T9MymXakOidcBZrB+V7GJDe/O2hzm+Zt739Wro7zd/m5aHPnXnPpw1N5eJR6W28om3tlYHLSzH84/h80trZJNdRo3I9vHlqHwaGzfmqsRl8Z0TnxywiIvErcRMRRWLgT8treCviT+0/GJ3ORaMyYjSirmGM4eTBqfz9y+ZNXu9sbODEQantvCq0Gv7DueWEl9ztl+bikaP3bFU28iCLXeMDnj42n32yovctbEROEvO/1Y+X19YxIN3NqUPan6uIiCSeuF/xNcbcaYyxEY9tYddNU58txph6Y8wcY8zYiHvkGWOeM8ZUNj2eM8bkdv9sJJF9sNXLXZ84T/E6pNDD747smV9KkekO0zc2tHnC2U43L6xwnM4G8PjkPAr28JSzgRluRua0DG5vOzSbE3YThO+J3BQXl+2fyWlD07SZTUSkB4r7wLfJamBA2GNc2LWbgRuBnwBHADuAmcaY8K3jzwOHAqcC05o+fq7rhy2JKhBxStjWugCXzilzrGTmJhuePS6/3Xq0iezYASmkhMWrG2oCrKrwt9l/xsYGXgyruwuhdIG9DVAjD9SYNiSVGw5SaTEREem8REl18Ftrt0U2mtCSzPXA7621Lze1XUwo+P0e8IQxZgyhYPcYa+38pj4/Aj4wxoy21q7urklI/Ju1uYHr51ewsSZAmtuQnWzI8rio8QUpbmg+hMEAT0X5T+3xJsPjYkr/FGZubs7tnbGpodUaubW+IDcudJZ2G5uXxK8OzdnrcVw3Lot3NjawoSbAxH7J/GVyXodLlomIiIRLlBXf4caYzcaYtcaYfxtjhje17wv0B2bs7GitrQfmApOamiYCNcD8sPvNA2rD+kgvZ63lL1/U8J2ZpWysCf2pvj5g2V4f5KsqP9vCTh4D+Nn4LKYO7vk5oCd3sKzZvUuqd71vEDre98/HtH86W0cNynCz+Nv9WHFuf96YVkhugh4MIiIisWd2l7MXa8aYU4EsYBXQF7gd2B8YC4wmFMTuY63dEPaaZ4BB1tpTjDG/AC631g6PuO8a4Clr7e/C2ysrK3e9IUVFRV0zKYkr/iA8sMbDy9s6VrZqQm6Ah8d66WA52oS2pcFw5uLmahUuLPeOaeS4guYg98saw0VLUgmEnff23YE+bhjuLGcmIiLS1UaOHLnr45ycnBY/qeP+77TW2rfDnxtjFgJrgIuBhTu7RbzMRLS1Ft1H9mkh/M1LFEVFRQk57mjp7PwrvEEumVPGnG0dK9V1QF4S/5rWf483a3W1aP/3Hwkc8NV2vmjK7Q1iuGVVCg9MyOXS/TMIBC1XvlVMgOYgd3CGm/tPGECmp/tXZnvz139vnjto/pq/5q/5d2z+cR/4RrLW1hhjVhD6mfzfpub+wMawbn2B7U0fbwP6GmOMbVrebsoN7hPWR3qhNVV+zn+3lC8rnRu2sjyGvx6bz6T+yVQ3Wqp8Qaoag7iM4cA8T4/dzNaWWw7J5uLZZbueBy3csKCCLbUB+qW7WFzsXNm9b0JOTIJeERGR3Um4wNcYk0oo1WE2sJZQYDsV+Djs+mTgZ00vWQBkEsr13ZnnOxHIwJn3K73I9roAp75VzPaI3N0hmW7+fWIBY/NDaQ9ZHhhIfK7udpczh6Xxx0m5/HRBhaOqxQOfVxP5K8AZQ1M5bWjPOchDRER6lrhfljHGPGCMOdYYs68x5ijgJUJB69+bVnAfBn5ujPm2MeZA4FlCm9meB7DWrgTeIVThYYIxZiLwBPCGKjr0TsGm09cig94j+yTz3hl9dgW90uzi0Rk8f2I+aRGJzeG5Qlkew70TemZNYxER6RniPvAFBgP/R6iW7yuAF5hgrV3fdP0+4CHgz8BiQnV+T7bWhp+ZegGwlFD1h+lNH1/YLaOXuPPkylpmbXbm9H5neBqvTSukT1rvXt1tz7Qhabx+aiH5bVRVuP3QbAZl6P0TEZH4FfepDtba83dz3QJ3Nj3a6lMGfD+qA5OE9EW5jzsWVzradtaGdbt6V+7unji8TzIzTi/k7BmlrA8rX3ZIoYfL9+9ZxzaLiEjPkwgrviJR0eC3XP5+Gd6wE3WzPYYnpijo7Yz9cjzMOL0Px/RPBkJVHJ7UeygiIgkg7ld8RaLl159W8kW5s4LDgxNzGZqp/w06q1+6mzdO7cPWugCFqS48CnpFRCQB6Ce+9DjWWioaLdW+IDU+S63PsrLCx2Mrah39zh2exndGpMdolD3DgHTl9IqISOJQ4Cs9RmlDgGc2JvHap9vZXBdot++QTDf3T1QFAhERkd5Ega8khM21AT7e0YjHBfvlJDEsK4mUptJaX1b4ePyLGv7vqzoaAslA+0Gvy8CTU/LISVaKu4iISG+iwFfiWlVjkAeWVvP4FzX4wsruugwMzXRTkOLikxJf2zdoxY0HZTGxX0qURyoiIiLxToGvxKWgtbzwdT13Lq5scdBE6Dqsqw6wrrr11d0UNxSmuMn0GDI9hgyPi2yP4diBKSq7JSIi0ksp8JW482lxIzcvqmBxcedWcgH6pQS5ZlwuF47KILeNgxZERESkd1LgK3Hln0W1/OTDCsdRuDsNSnczIieJryv9LTavHVbo4eqxmYzxb2bMqKzuGayIiIgkFAW+EjeWl/m4YX7LoDfFDT85MIufjsskwxNaxa31BVlTHWBzrZ+hmUmMyU3CGENRUfePW0RERBKDAl+JC/V+yxXvl9EYkc57xtBU7j4yh2FZzi/VDI+LcfkuxuV7unGUIiIiksgU+EpcuGNxJSsrnKeqPXpMLt8fqY1oIiIiEh3a/SMxN2NjA0+udJ6qdv6INAW9IiIiElUKfCWmiusDXP1huaNtn0w3903QqWoiIiISXQp8JWastVzzYTnFDc2JvW4DTx2bR7ZOVRMREZEoU3QhMfP0qlqmb/I62m4en8WRfXWqmoiIiESfAl+JiaJKH7d/XOloO6pvMjcepBq8IiIi0jUU+Eq3CwQtV31QTkPYGRRZHsMTU/JIcpnYDUxERER6NAW+0u0eWV7DxxHHEd97VMtavSIiIiLRpMBXutUX5T7u+azK0XbqkFS+u196jEYkIiIivYUCX+k2vqDlx3PLHaez5ae4eHhSLsYoxUFERES6lgJf6TYPLK3m8zJnisODE3Pol+6O0YhERESkN1HgK91iSUkjDy6tdrR9e980ztpXKQ4iIiLSPRT4SpfzBixXflCO3za39U1z8cCEnNgNSkRERHodBb7S5R5eVs3KCr+zbVIu+alKcRAREZHuo8BXutTGGj8Pf17jaPvufumcNjQtRiMSERGR3kqBr3SpX31cRX2gOcehIMXF745UioOIiIh0PwW+0mU+3Obl1XX1jrY7Ds8mN0VfdiIiItL9FIFIl/AHLbcsrHC0HVzg4QIdVCEiIiIxosBXusTfv6xlRblzQ9u9R+XgdumgChEREYkNBb4SdeXeIHd/6jyW+NwRaUzolxKjEYmIiIgo8JUucM+nVZR7mze0ZSQZ7jpcG9pEREQkthT4SlStKPPx9OpaR9tNB2cxQMcSi4iISIwp8JWosdZy60eVBMNOaNs3y81VYzNjNygRERGRJgp8JWre2+Jl7lavo+2eI3NIcWtDm4iIiMSeAl+JiqC13LnYuaHt+IEpTBuSGqMRiYiIiDjFfeBrjLnVGPOxMabKGFNsjHndGHNgRJ9njTE24rEwok+KMeYRY0yJMabWGPOaMWZw986m53p5TT3LynyOtjsPz8YYrfaKiIhIfIj7wBc4DngMmAScAPiBd40x+RH93gUGhD1Oi7j+MHA28F1gMpANvGGM0a6rvdQYsPz2M+dq79n7pnFwQXKMRiQiIiLSUlKsB7A71tpTwp8bYy4EKoGjgdfDLnmttdtau4cxJge4DPiBtXZm2H3WAycB07tg6L3Gs6trWVcd2PU8ycBth2bHcEQiIiIiLRlr7e57xRFjzABgCzDZWvthU9uzwLeARqACeB+4zVq7o+n6CcAsoK+1tjjsXiuAl6y1d+xsq6ys3PWGFBUVdfl8El1dAM5anEaZrzml4ZwBPm4Z4WvnVSIiIiLRN3LkyF0f5+TktMi3jPsV31b8EVgCLAhrewd4BVgLDAPuBt4zxhxmrfUC/YEAUBJxr+1N11oV/uYliqKiom4d971LqijzVe96np5k+O2UIfSLUd3e7p5/vNH8e+/8e/PcQfPX/DV/zb9j80+owNcY8xBwDHCMtXbX39attf8O67bMGPMJoTSG0wkFxG3eEkisJe84Ulwf4JFlNY62q8ZmxizoFREREWlPImxuA8AY8wdCG9NOsNauaa+vtXYLsAnYGf5vA9xAYUTXvoRWfWUPPLC0mhp/8+8N+Skurj1Qh1WIiIhIfEqIwNcY80fge4SC3lUd6F8IDAK2NjV9AviAqWF9BgNjgPlRH3Av8FWlj2daOZo4OzkhvqRERESkF4r7VAdjzJ+BCwltXis3xuzMya2x1tYYYzKBO4GXCQW6w4DfATuAVwGstZXGmKeB+40xO4BS4CHgc0Jl0KQTrLXcvLASX7C5bUimm8v2z4jdoERERER2I+4DX+Cqpn9nRbTfRSjgDQDjgIuAXELB72zgXGttdVj/nxKqAfwCkNZ0v4vCc4WlY15f38B7W5xHE99xWLaOJhYREZG4FveBr7W23WjKWlsPnNJen6Z+DcBPmh6yh2p9QX7xUaWj7Zj+yZy9b1qMRiQiIiLSMUrIlE556PNqNtU6D6u4f0KujiYWERGRuKfAVzrsq0ofjyx3li/78Y4VGPQAACAASURBVAGZjMnzxGhEIiIiIh2nwFc6xFrLLYsqaQzb0NY/zcUth2TFblAiIiIinaDAVzrkjQ0NzNrs3NB295E5ZHn0JSQiIiKJQVGL7FadP8iti5wb2o7WhjYRERFJMAp8Zbf+tKzGsaHNbeABbWgTERGR/2/vvMPtKqr+//nmhvQGoYeELqFIBwVEgiGKKE1UQPorIhZQQeV9UYqCXQNRsPwQBQRpNkRRmhINYKGEIl0IBEhCEtJJIcn8/lhzkn3PPbfk5p5zbnK+n+c5z917ZvbMWjOzZ6+9ZmbfNQwbvqZNZixaxmXe0GaMMcaYtQAbvqZNxj46j/lL04rz9Xr34JxdvaHNGGOMMWseNnxNq7w8fylXPrWgWdhZOw9gUC93G2OMMcasediCMa3y7UfmsbjwD52H9Wvi1JED6ieQMcYYY8xqYMPXVOS5OW9y3bNvNAs7Z7eB9OnpDW3GGGOMWTOx4Wsq8rWH5rFs5dJeth7UxEe26Vc/gYwxxhhjVhMbvqYFj8xcwm8nLWwW9qXdBtGzh729xhhjjFlzseFrWnDxg3Obnb91vXU4wv+swhhjjDFrODZ8TTPum7qYO8v+NfH5ewyih/9ZhTHGGGPWcGz4mhXMXrycz943u1nYPhv14qBhveskkTHGGGNM12HD1wDw5vLESX99nWfmLG0Wfv4eg/yviY0xxhizVmDD15BS4uz7ZzN+SvMlDkdv3Zd9NrK31xhjjDFrBzZ8DT94fD7XPNP8m717b9CLcfuuWyeJjDHGGGO6Hhu+Dc7vJy3kggeaf8Vh8wFNXDd6Pf+zCmOMMcasVdjwbWAenrGEj/9tFoX/U8GgXuLGMUPZoG9T3eQyxhhjjKkGNnwblH9OW8yH7pzJwsK/Z2sSXHPgeowcsk4dJTPGGGOMqQ496y2AqT2/fHYBn71vNkuWNw8fu88QRm3apz5CGWOMMcZUGRu+DcSy5YmvPjiXcY/PbxF3xk4DOGm7/nWQyhhjjDGmNtjwbRDmvbmcj42fxZ8nL2oR9/mdB3Lu7gPrIJUxxhhjTO2w4buW8/qiZfz0qQX85IkFzFzcfG1D7ya4bL91+dDW/eoknTHGGGNM7bDhu5Yyef5SLv/PfH7xzBssWJpaxG/UtwfXjR7Knhv0qoN0xhhjjDG1x4bvWsS0N5Zx67QmJk6eyW0vLaKCvQvALkPX4ZejhzKsvz9ZZowxxpjGwYbvGszS5YmHZizhjpcXc+fLi3hk5ptAb6DlOl6AQeuIU7fvz9k7D6T/Ov6SnTHGGGMaCxu+axgvz1/KX15dzN2vLOKeVxczZ0krbt0CG/ftwSd3jK82DO5lg9cYY4wxjYkN3zWARUsTP3t6Adc8s4CnZi/t8HXbDu7JGTsN4Oit+9G7yf9+2BhjjDGNjQ3fbsybyxPXPvMG33lkLq++sbz9C4AdBizj0G2GMGazPuy+/jr0kA1eY4wxxhiw4dstWbY8cdPzC/nWxLlMmreszbTr9hYHbtqHMZv1YfSw3sx5+Xm23XZQjSQ1xhhjjFlzsOHbzXhg+hLOmDCLJ1tZ0tAk2GuDXrxrWG9GD+vDrkPXoanHSq/unFoJaowxxhizhmHDt5uwZFni2xPnMfaxeSyvsF+tdxOcOnIAn33rADbo68+QGWOMMcasKjZ8uwH/ef1NTv/7LB57/c0WcT0FJ76lP5/fZSCb+ru7xhhjjDGdpqG+bSXpk5JekLRI0oOS9q+nPMtTYtxj8zjw1tcqGr0f3rovDxy1EWP3HWKj1xhjjDFmNWkYj6+ko4FxwCeBCfnvnyTtkFJ6qS4yAfdNXcySsg82bNy3Bz94x7qM2axPPcQyxhhjjFkraSSP71nAVSmlK1JKT6aUzgCmAJ+ol0CS+P5+67Je75XNcNSWfbn/yI1s9BpjjDHGdDFKqf3//LWmI6kX8AZwbErp5kL45cBOKaUDSmFz5sxZUSHPPvtsTeS7a0YT33yuF+dsvYQxG7T9+TJjjDHGGFOZbbfddsXx4MGDW/wzg0ZZ6rA+0ARMKwufBhzU2kXFyqsm224Lx+6+nCG9V98B/+yzz9ZM7u6I9bf+jap/I+sO1t/6W3/r3zH9G2mpA0C5e1sVwupCVxi9xhhjjDGmdRrF2poBLAM2LgvfkJZeYGOMMcYYsxbSEIZvSmkJ8CAwpixqDHBf7SUyxhhjjDG1plHW+AKMBX4h6V/AvcDpwKbAj+sqlTHGGGOMqQkNY/imlG6UNBT4MrAJ8DhwSErpxfpKZowxxhhjakHDGL4AKaUfAj+stxzGGGOMMab2NMQaX2OMMcYYY2z4GmOMMcaYhsCGrzHGGGOMaQhs+BpjjDHGmIZAKXWLf1zWbZgzZ44rxBhjjDFmDWfw4MEqD7PH1xhjjDHGNAQ2fI0xxhhjTEPgpQ7GGGOMMaYhsMfXGGOMMcY0BDZ8jTHGGGNMQ2DDt4uQ9H+S/i1prqTpkm6VtFNZGkm6UNKrkhZKukfSjmVpviTpXkkLJFVchyJpnKQHJC2SNKmKanWIWukuaRdJ10uanPN4WtIXJNW1H9dQ/w0k3Z7zWJzr4XJJg6utY1vUsu8X0q4v6RVJSdL61dCro9T43k8VfqdXU7/2qHX7Szpe0sQ8/s2QdE21dOsINbz/T26l/ZOkvaqtZ2vUuP/vJekuSbMkzZZ0t6S9q6lfe9RY/9GS7pM0T9IUSd+S1LOa+rVHV+gvaQtJV0p6Psc/L+kbkvqW5TMi578g3/vfl9RrVWW24dt1jAJ+COwLvAtYCtwlab1Cmi8CZwNnAHsBrwF3ShpYSNMb+A1waRtl9QCuBuo64BcYRW103wOYDpwA7AhcAJwP/G9XKdJJRlEb/ZcDvwUOBd4CnAyMBq7oIj06yyhq1/dL/ByYuNqSdw2jqK3+HwM2KfyuXn0VVotR1Eh/SWcC3wG+C+wEHAjc0lWKdJJR1Eb/G2ne7psA1wIvAA90kS6dYRQ10F/SAODPwKu5rH2AKcDtZfnUmlHURv+dgduAO4DdgGOAw4BvdqEunWEUq6//SKAJ+ATxbD8DOBEYV8pAUhPwR2AgsD9wLPBB4HurLHFKyb8q/IABwDLg0Hwu4ib9UiFNX2Ae8PEK138wmqfNMj4PTKq3rvXQvZD228CD9da5jvqfCUypt8611B/4DHA3McgmYP1661wr/bO+H6y3jvXQHxgCLADG1FvHerV/Wbp+wGzg3HrrXKP23zP3/y0LYVvmsD3rrXcN9P868HBZ2KHAQmBgvfXuKv0LaT4JzCycv5dw/gwvhB0PLAIGrYqM9vhWj4GEZ3ZWPt8S2Jh4WwMgpbQQ+BvxprQ2UUvdBxXK6S7URH9JmwIfAMZ3WtLqUDX9Je0GnEN4A5Z3hbBVoNrtPy5P8/1b0umq81KfClRL/3cTXqGNJD2hWOryW0lbdY3YXUatxr8PA/2J2Y/uRLX0f5qY8fuopN6SehOzHy8B/+kCubuKaunfmzDyiiwE+hCzod2FrtK//Nm+D/BkSmlyIex2ol5WSf/uNmCuTYwjpmLvz+cb57/TytJNK8StLdREd0m7E9P9P+psHlWiqvor1jm/AbxCvDWf0kk5q0VV9JfUH7geOCOl9MrqCllFqtn+5wNHAwcBNxDTfOd2TsyqUS39tyKeWV8GzgKOBNYB/iqpX6el7XpqNfafBvwhpTRlNfKoBlXRP6U0j5hWPxp4I/+OJmYAFq6GvF1Ntdr/duBtkk6Q1FPSMGI8gFj20l1Ybf0ljSBmtH9YCN64Qh4zCO/yKt1HNnyrgKSxwDuAo1JKy8qiyxetq0LYGkutdJe0HbHe59KU0q87k0c1qJH+nwN2B44gjIGOrImtCVXW//vAvd2pvcupdvunlC5KKU1IKU1MKX0P+ArwhU4L3MVUWf8ehKF7ZkrpzymlfwHHARsSU751p4bj346EB6ze6/ubUU3980annxEG1duB/YCHgVvyS3Hdqab+KaU7CGPwMsLz+wyx5hfC+Ks7XaG/pI0II/9O4JJ28mgvvCI2fLsYSZcQi67flVJ6vhA1Nf8tfzPZkJZvMWsktdJd0kjgHuCGlFK9N7atoFb6p5SmppSeSindAnwcOE3S8M7I3JXUQP/RwMmSlkpaSqzzBZgq6WudkbkrqdO9/09gUH5Y1JUa6F/ybD5RCkgpzSE2O41YNWm7nhq3/2nAZGKzV7egBvp/BNgaOCWl9O+U0j9y2AjC+19XatH+KaWxxFr3EcD6rNzY+cIqC9zFdIX+kjYG/go8DpyQ8kLeQj7leaxPLH9apXq04duFSBpH3IjvSik9VRb9AtFwYwrp+xC7E++rmZBVola6S9qBMHpvTil9bnVk7krq2Pale7j3auazWtRI/3cDuwC75t+pOXwU4Q2uG3Vs/10J78/s1cxntaiR/vfmv9sV8hlATPO+2Amxu4xatn++9gTgZymlbrHOvUb69yM8e0Wdl+ewen/Ssmbtn4JX8/KOY4kXoIc6K3tX0BX6S9qEeLY/CRybUlpals/9wPaSNiuEjQEWAw+uksD13P23Nv2Ay4G5xE7zjQu/AYU05+Q0HyA+xXMD4a0YWEgzgniYfZ64oUsP+WI+2+Swsfn6Uppea7PuxGdOpuVri+Vs3AhtD7wfOClfvwXwPsL7dX8j6F+h3FF0g6861LD9DyU28+xEeL5OBeYA4xpB/5zmd4Q3aD9gB+BmYBLQrxH0z+mOJ6a2R9Sz3evQ/0cSL3k/ArYnnge/yPfAZmu7/jnNF4C3Zt3PA5YAR6zp7Q9sSizduAcYXpZPU07TBDwG/IX4nNtBxD6XH6yyzPW+adaWX+6olX4XFtIIuJCYsltE7MbfqSyfq1rJZ1QhzT2tpNlibdY9X1+xrEZo+3yj30949xbmgeJbwLqNoH+FckfRPQzfWrX/wcSaxnnEZ70eIz7t1rMR9M9pBgI/BV4ndnzfCmzdKPrndOOB2+qpcx3bfwwwgRgDZxHT4vs2kP5/YeX4/w/gvWtD+xOb1FvLZ4tCuhHAH4iNjTOBHwC9V1Vm5cyMMcYYY4xZq/EaX2OMMcYY0xDY8DXGGGOMMQ2BDV9jjDHGGNMQ2PA1xhhjjDENgQ1fY4wxxhjTENjwNcYYY4wxDYENX2OMMcYY0xDY8DXGGGOMMQ2BDV9jjDHGGNMQ2PA1xhhjjDENgQ1fY4wxxhjTENjwNcYYY4wxDYENX2OMMcYY0xDY8DXGGGOMMQ2BDV9jjDHGGNMQ2PA1xhhjjDENgQ1fY4wxxhjTENjwNcYYY4wxDYENX2OMMcYY0xDY8DXdDkn3SDq1xmWeK+mnq3jNDpIe6ERZP5Z0Xj4eJenlVc3DrBqSXpY0qsZl/lTSuW3EnyXp4nbyuFbS8V0vnak2knpLmi9p005ef7qku/JxH0lJ0matpP2opFs7kraTshws6bmuys/UHkkjJS2ttxzdARu+azmSJklamAfg+aszEK+pSLqqPQMjpfT1lNKpOf0W+cHRs52sLwK+WyhnkqQlktYvK39izm+LXNbpKaWLOqNLWb7Ftp2a9RxQId2Fufy9V7fM7oikn0n6nxqXOUHSyW2lSSmdmlL6ek5/kKRJZUl+DJwiaehqyHGapKclzct94A+S+nc2P9N1pJQWp5QGpJReBZB0g6QvV6msK1NKh1Yj70ZD0t8kvbMK+Vat/buatf1Fx4ZvY3BoHoAHFAfiIh0w8tZIJDVVKd9NgAOB35VFvQAcW0j3VqBvNWTIHJpSGgDsCuwG/F+ZnAJOAF4HTqqGANWq4w6WLeA9wJ9qVF4PSV0ybqaU3gDuINqnM7KMBr4CfDilNBDYEfhVV8hWVk6X6WxMd37WSBoMbA/cV29ZTPXwYNagFLyaH5X0EvCXHH6YpP9Imp2XHGxfuGaSpC9IelTSAklXStpI0p+yx+kuSeuuggyHZ2/oXEn/lXRwIXpzSffmfO8oelEl3Zy9W3Py2/mOhbirJP1I0m2SFgAfBY4Dvpg9o7e2IsuFkq7Np3/Lf2fna/apcMkY4KGU0qKy8F8AJxbOTwKuKSurVQ+0pDMlPVGappT0/lxHsyXdJ2nnStellKYCtxMGcJH9gU2BzwDHSOpVKOuRspmApLwcYBXr+MAKemwu6XeSpkuaIWlcDu8h6XxJL0p6Lec1KMdtk2U4WbE04XVJH5P0NkmP5ToYV1bUbsC0lNKUnMfHJT2V+83jknYppN095zNH0vWSeudrhmZdpkuaJelWScMKukyQdJGk+4EFwPXAPsCPc71d2kpbXpv71WDgVmBEoa43zMnuAd5X6foOsBdwb0rpEYCU0syU0lUppQW5/CFZhun53v0/ScpxF0u6qiDrNpJSGzqPyPV0laQpuZ5+XUh/WO5Ps/O1O3VEAUlfk3RJPu4taZGkkpd8QD4fLOk6SZ/J4ZvnfnJaPh+ZdZRizDg8h4/K6d6dzw9WK0uTJPWX9H1Jk3P/GC+pZ/79WtK0rNtfJW1XuO4GST/I4fMk3V3qOyosOZB0JnAUcF5u/5tzmvMlvVDor+31hSNyW07PdVdqzxXLIirodmDWa798vpOkv+Q2fFLSEYW0hxfun8lZ7mJe5+ayX5F0XCH8yNz+cxX39rmFuJGSliru5cnAbTl8f0n/zPX6UEm+VnTYpSDzVEln5/C+ki7PffJlSd+RtE6OO1jSc5K+rBiDXpF0SNbxv5JmlvIp8B7grymlpZL2k/Rw1mmqpG/kfO+W9LEy+Z7J5TVJuizX0ZxcJ9u10f7DJd2S5Xte0umFPL+p6Pc35msmStpS0gU5/SRJLcbewvU9JY3Lej5HPLOK8cWx8jnlWTPFDNRvga20crwamuvjn1mvVyVdom78EtMmKSX/1uIfMAk4qEL4FkAijLL+hFfyLcRDbgywDvBF4DmgVyGvfwAbAcOA14CHCOOjN2E8X1Ao41HgI63ItTcwJ5fVI+c3MsfdA/w3y9M3n3+zcO3/AANzmZcCEwtxV+V898v59slhF7dTTxcC15bVTc820n8HuLxSXQNPE16DJmAysHnOb4uCjBfn41HAy/n4vFyfG+Tz3XMdvy3ndVIuo3d52wKbAY8B48pkuhK4KbfnTOADrehzGvAUMKgzdVyWV0/gcWIZSKlv7Vco5xlgy5z/LcDPc9w2uZ4uy+UeAiwkBuENso4zS3nla74MXJSPj831vQcgov8Mz3EvE313Y2BoluHUHLcBcGSWcxDwG+BXhTIm5LrePtdjzxx2cjt96lrgwnx8EDCplfvgtXbyOL6VuFG5fi4A9i31i0L8L7MuA4GtiHv5pBx3MXBVIe02QGpH59tznusCvYB35rR7AdPy36bcd/7LynHjJ8D3W9Hh3cDD+fid+bp7C3EPFvrNb/PxiTnddYW4X+fjrwOX5OPzc7qvFeK+14ocVxLe942zDvvnvz2J+24AMZb8CPhH4bobgNnEi1AfYvnKXTmuD9GfNyuk/XJZuUcDmxD30QnAPGD9HHd6hbxuB4YQ98/zpb7RStrNgMOAl4DdctwgYArhDGjKbfY6sE2OnwnsnY+HFq47GHgT+FLuD0dmWQfk+NHEjEMPYtx6HTg4x43M8vwU6EfcZ1vksg7K1xwCTAfWrdA26+a4TxPjwiBgrxz3beDvwPrEc+nfwJfKZD4nt+MZwFTCOdGfeG4tAoaVjW2le+Rh4EP5eCDwtkL/G1+45m053ybgcOD+LGOPXCcbVmr/nP6xLF8vYrx6CTggx38TeINwLPQEbiRmFD9f0OfJNsaOz+b8NyXGuAnA0kL8YUQ/Um6HhcCOhbp7rsJYVbrHtybGk9PbGgO766/uAvhX5QaOh9d8YnCeDfwuh2+RB6OtCmnPA24qnPcAXgFGFfI6rhD/a+BHhfMzSvl3QK6fkB9QFeLuKRsgPgn8uZW0Q7Ieg/P5VcA1ZWmuousN3ysoGOOF+jmIMMa+kQePO/Mg1Zbh+wowNg9Mgwv5/Yhs1BXCni4MjKW2nZfzvxsYUkjbD5gLHFGo81sq6PIOwsB+S2fruCz9/uQHQYW48cBphfMdgcW5r5UM340K8XOAowrntwCfLpzfD+yTj+8GPtWKTC8DxxTOxwKXtZJ2T2B64XwCcH5Zmq4yfLcHlrSTR0XDN8e/D/hDrqd5xAtZD8I4WVpsU+BTrDSOOmL4nl84H57zG1xBhisovPDmsP9SeEFpQ/7+uf2HEPfNOcCrue9+DRib021HGEoiDKjTgJdy3HXAmfn4PcRMDMBdwKnAhHx+L3BYBRnWIQyk7Tog78bAcvLLHmHMFOtxvdyHN6ADhm+F/J8C3pOPKxmzowppzwL+2EbacwhDaWThmpOAO8vKvBo4Jx9PA04BBpalOTj3sR6FsLnArq3o8WPgG/m4ZPhuWoi/ALii7JrxwNEV8joFuL+Vcl4B3lU4Pxx4qpLMuU0SsEsh/X9YaaCLeCkoGar/Igz9oRX67FxgRD6/rNBPD8l57g2o7Lpyw/cA4NmyNF8hP1MJw/fWQtyHyPdAmT59W6mb+yiMUYShu7RS2hz/Z+Djhbp7rrW0Oc3/Ate3d890x5+XOjQGR6SUhuTfEWVxkwvHmwIvlk5SSstz/LBCmmmF44UVzltsrmqF4cTDsTWmFo7fKOWbp5K+maeq5hLGH8Qbf4miTi2QdFxhCqeza0NnEV6ASvwC+AhwMmXLHFphCPEg/0ZKaU4hfHPg7DwVOFvSbKLeipsTj0ixvnMU8YAp1sORhLFyWz6/DnivpA1KCSQNJzzCJ6WUnslhq1vHwwkjb1mFuGZ9LB/3IgZxAFJKHepjeUpuK+CfhXI706f6K77A8FLW9y801xXa71PnFfrUZW2lLWMg8ULaKVJKf0wpvZ/win0A+BhhKGxIeGbK63pYi0xap6jzcGBGWf8ssTlwTlk/3aQjZaVYlvEw4e19J/HS+w/Cg3oAYQyRUnoaWAK8lXix+j0wQ9LWxXSEcbtj7uM7EUbdVrmv7EF4B8vZhHg5fb48Ik8XfzdPQ88lDFMR3tASK+oppfQ68TLaoQ3EiqVmjxbqbRta9r0ixTZ5sZ1yziJe5p8qhG0OvLOsrY4i6gDgiHz+Ul5asGfh2un5mVCieA/tp1geMl3SHGLsK+qxPDXfW7I5cHyZHHu2ok/F+1qSiBeRtvp4UeaF+W9rz6w9idm31/L5ScDOwDN5iv89sKLP/gY4Li+rOJoY8yH2GlxJOBmmSfqhKmw4LtTBFmV1cFbWqUS5rNNTtjoL+vRXbJ4tjT8P5vBNadlfVqBYnvQvxZKy2cC7aKPvKb5i9CfFsp+5xIxKW32122LD16TC8avEzQisGFiGE2/VXc1kYrpkVfkI8VZ/EDCY8M5CPIxKpLJrmp2nlK5LKzf6vbdCGeXXV+JRYmqq5cUpvUh4Wg4hBsj2mAW8H/i5mq9zm0xM0w4p/PqllK6vUOZ4whP73ULwScSg/pKkqcDNhHfrWIj1ccTmvEtTSsUXgM7UcZHJxBrtSpvemvUxYARh0ExvI7/WOJjwXpUebJ3tU18kpvz2TikNIh4A5bTXpy4q9KlPd+D6EtsDj6yqwC0yT2l5SulOwnDcifDgL6NlXZfu5QWEV7VE8WG7ItvC8WRgfeX12GVMBr5SoZ/e1EHxxxN97a3Ag/n8vbQ0VP8GHEN4pqfmdB/NejxGRMwHJgKfI5bnvEm8GJ1NeAJnVSh/CvGCuFWFuFOIJRcHEvfCyBxevBeGlw4krUfcc1Mq5NWsD0h6C/AD4qV3vZTSEGL6WBWubVEW0Z4tNioXOJIwLk8vhE0G7ihrqwEppc8CpJTuzy9SGxFLP1qMNa1wEzEVPzylNJgYi9oaLyYDPy2To39K6ZIKeVe8r7MBOJXW+/iqcgjwx0L+T6aUjiZeIr8P/EYr90hcDRxPjEHTUkoPl2RKKY1NKe1GGM27EPsroHIdPFVWBwNTSkeuquAppbsK488eOXgKLfsLEC/7xPPgIsLDPYR44S+1WaXx6gpiKd7WeZz8Km331W6LDV9T5CbgfZJG5zfZs4lpyGrscL2S+JTTaMWGp2GSRrZ7VXjIFhNTPv2IdXvtMY3KD7XWmE5MZ7Z1zZ3EZqk+rcR/lJiCW9CRAlNK9xDr7n4r6W05+ArgdMXmLmXP5PskteZpvhQYI2lXxQab0YRBvWv+7QJ8i5Vfd/gZMfB+uyyfztRxkfvztV+X1E+xAaVk0F8PnKXYXDmQmM6+vsyT1FHex0pvNsQU+Bcl7Zbra9vs0W6PgYT3alb2DJ7fgWtWtU9NIwzH8rY7gE5+kUKxoejDktbN+r6d8Ib+Ixt8vyLaYICkLQljsLSBcyJwgGJzzRBi2rJVUkqTiaUDlys2za2jlZ98+n/ApyTtleUYIOlQdfyzauMJD+GjKaWlhPF+GvBM9qAW032ald7de/L538v6T2vpxlOBXFfXAOMUm3WbJL0jv7gNJNaBziSmuCttSj0836O9c/xfC17DIuV9ZgAxzkwHemQDdZtKMhY4R7HZb4us041tpH2JGAPO1crP/f0O2E3S0bkNe0l6u6S35PHlmPxy8yaxdKbSrE0zsoNkADAzpbRI0r7EtHxbXA18KI//TXmMGC2p0gvY74BtJH0iyztI0l457nrgAsXmqw2JpQnXVsijIzQbTySdKGlonrmaQxiDpX52T9b5axRm9XJd7qnY9LWAeKkv1WF5+0/I13xWsRGyp6SdJe3eSfnLuQn4nKRNFJvDv1iI60s4QV4Dlks6jJg1LDEN2FDNvdUDgTkppfmKzc7NNvitSdjwNSvI04nHE16IGcChxOeylnQmP8XXIY6rFJdS+hfhTbmEGFTG0/zNvTWuIaZsXgGeIKZF2+NKYAfFdFL558cqyfYGMaDdm695e4U004g35MNbyeO/KaVV+ucW2WN3CvB7SXvk6z9GrCGbRXiDTm7j+ulE/ZxHbJSZmFK6I6U0tfQjPBc7K3bdHwMcqeZfdtifTtSxYqr/1izHUsLgA3ISCgAAAqpJREFU3p7warwEfDAnvYJ4WP+dmFqex0qPSIdRfF5rNLHZp6T/9YRhfyOxBu83xBKA9hhLePNmEi95HTFELwWOzf1jbHuJU0qPE2viJ+VrNlR43A+mY8thKjGbWNv5HKHv1cDXU0olY+iTxIP3BeL+urpQ1p+JTYOPEWsZf9+B8kr/SOMZ4sF4Rtbtn8AniDXps3L8in+6oVhG0tbyjwnEC1bpayqPEobX38rSjScevqXwvxPGR3vpys9L31WeUbjmTGI6/WGiH1xEeLOuJAzTqURdTagg/7XEeswZRJ8/qUIaiBeEvXL735BSeohYC/sA4Z3bMh+3xR+JGYIHCI9dm0ZeSul54j75qqQTssf7PcQ4M4XwGF9MGEEQGxNfJMbkE9vQpVhGIvrhdyXNIwysmzsg11HEmtYZuczPkG0SxddDLs1pZxGboI8hDLWniX0JEC+pTxDraicSS13KX+TbRbE0ZgTN6//9wNNZp28Qnw1cWtD5F8QehV8WrhlCeLtnE+Pbi8SYCy3b/03Cy7xvTjeduIc6ulywPS4j7pH/ELMeK2ZgUkoziE1ytxL9/QiaOxEeIcaEF7O86xEvzqdKmg9cTtsvXd2a0iJpY8wqImkHwpjYO/lGqinZq/TdlNK+9Zals0j6HPEFj7b+u9u1xMbOznqxTBWRdAPweEqpzX+QY7o3kk4kvpBzYruJV15zGmEMH1Q9yUw1WDO/wWZMNyCl9ATxeRdTe5YT3qI1llbWMxpjas8MYqazQ+RlPJ8gPMFmDcOGrzFmjSOl1JElLmsDvyGWMhhjqkRK6bb2UwV5Pez1xNKALv9Piab6eKmDMcYYY4xpCLy5zRhjjDHGNAQ2fI0xxhhjTENgw9cYY4wxxjQENnyNMcYYY0xDYMPXGGOMMcY0BP8f3khks0ku8ooAAAAASUVORK5CYII=\n", | |
"text/plain": [ | |
"<Figure size 720x432 with 1 Axes>" | |
] | |
}, | |
"metadata": { | |
"needs_background": "light" | |
}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"# bike age by hours ridden\n", | |
"ax = bike_share_bikes_mins.groupby(level=1).median()['Minutes'].plot(figsize=(10,6), fontsize=14, legend=False)\n", | |
"# ax = bike_share_bikes_mins.groupby(level=1).quantile(0.75)['Minutes'].plot(figsize=(10,6), fontsize=14, color='r', label='25th Percentile')\n", | |
"ax.set_xlabel('')\n", | |
"ax.set_ylabel('Hours', fontsize=14)\n", | |
"\n", | |
"# formatter function to format y-axis labels to have a comma at the thousands place \n", | |
"def func(x, pos): \n", | |
" s = '{:0,.0f}'.format(int(x) / 60)\n", | |
" return s\n", | |
"\n", | |
"years = mdates.YearLocator() # every year\n", | |
"months = mdates.MonthLocator() # every month\n", | |
"yearsFmt = mdates.DateFormatter('%Y')\n", | |
"\n", | |
"ax.xaxis.set_major_locator(years)\n", | |
"ax.xaxis.set_major_formatter(yearsFmt)\n", | |
"ax.xaxis.set_minor_locator(months)\n", | |
"\n", | |
"plt.xticks(rotation=0)\n", | |
" \n", | |
"y_format = tkr.FuncFormatter(func) # make formatter\n", | |
"ax.yaxis.set_major_formatter(y_format)\n", | |
"ax.yaxis.set_ticks(np.arange(0, 105000, 15000))\n", | |
"# plt.legend(fontsize=14, loc='best')\n", | |
"\n", | |
"ax.text(0.08, -0.15, \"From: chart-it (MikeRAzar.com/chart-it) | Source: www.capitalbikeshare.com/system-data\", \n", | |
" fontsize=12, transform=ax.transAxes) \n", | |
"\n", | |
"plt.savefig('/Users/mazar/Desktop/personal/chart-it/posts/Bikeshare/pt 4/bike_age_hours.png', dpi=200, orientation='landscape',\n", | |
" bbox_inches='tight', pad_inches=0.1)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 24, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"25080173" | |
] | |
}, | |
"execution_count": 24, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"len(bike_share)" | |
] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python 3", | |
"language": "python", | |
"name": "python3" | |
}, | |
"language_info": { | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 3 | |
}, | |
"file_extension": ".py", | |
"mimetype": "text/x-python", | |
"name": "python", | |
"nbconvert_exporter": "python", | |
"pygments_lexer": "ipython3", | |
"version": "3.6.7" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 1 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment