{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Card-Heining-Kline event study example" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# Add PyTwoWay to system path (do not run this)\n", "# import sys\n", "# sys.path.append('../../..')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Import the PyTwoWay package\n", "\n", "Make sure to install it using `pip install pytwoway`." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import pytwoway as tw\n", "import bipartitepandas as bpd" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Get your data ready\n", "\n", "For this notebook, we simulate data." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "scrolled": false }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
ijytlkalphapsi
00810.2839380240.000000-0.114185
1081-1.6135041240.000000-0.114185
201751.0574882280.0000000.908458
30450.8943133220.000000-0.604585
40901.1223514240.000000-0.114185
...........................
4999599992-2.978308000-0.967422-1.335178
4999699992-2.532244100-0.967422-1.335178
49997999914-2.428863200-0.967422-1.335178
49998999914-0.928106300-0.967422-1.335178
49999999921-3.030404401-0.967422-0.908458
\n", "

50000 rows × 8 columns

\n", "
" ], "text/plain": [ " i j y t l k alpha psi\n", "0 0 81 0.283938 0 2 4 0.000000 -0.114185\n", "1 0 81 -1.613504 1 2 4 0.000000 -0.114185\n", "2 0 175 1.057488 2 2 8 0.000000 0.908458\n", "3 0 45 0.894313 3 2 2 0.000000 -0.604585\n", "4 0 90 1.122351 4 2 4 0.000000 -0.114185\n", "... ... ... ... .. .. .. ... ...\n", "49995 9999 2 -2.978308 0 0 0 -0.967422 -1.335178\n", "49996 9999 2 -2.532244 1 0 0 -0.967422 -1.335178\n", "49997 9999 14 -2.428863 2 0 0 -0.967422 -1.335178\n", "49998 9999 14 -0.928106 3 0 0 -0.967422 -1.335178\n", "49999 9999 21 -3.030404 4 0 1 -0.967422 -0.908458\n", "\n", "[50000 rows x 8 columns]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "df = bpd.SimBipartite().simulate()\n", "display(df)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Prepare data\n", "\n", "This is exactly how you should prepare real data prior to generating the CHK plot.\n", "\n", "- First, we convert the data into a `BipartitePandas DataFrame`\n", "\n", "- Second, we clean the data (e.g. drop NaN observations, make sure firm and worker ids are contiguous, etc.)\n", "\n", "- Third, we cluster firms by the quartile of their mean income, to generate firm classes (columns `g1` and `g2`). Alternatively, manually set the columns `g1` and `g2` to pre-estimated clusters (but make sure to [add them correctly!](https://tlamadon.github.io/bipartitepandas/notebooks/custom_columns.html#Adding-custom-columns-to-an-instantiated-DataFrame)).\n", "\n", "- Fourth, we convert the data into extended event study format\n", "\n", "Further details on `BipartitePandas` can be found in the package documentation, available [here](https://tlamadon.github.io/bipartitepandas/).\n", "\n", "
\n", "\n", "Note\n", "\n", "In general, the CHK event study is generated by clustering the data using quartiles of firm-level mean income.\n", "\n", "
" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "checking required columns and datatypes\n", "sorting rows\n", "dropping NaN observations\n", "generating 'm' column\n", "keeping highest paying job for i-t (worker-year) duplicates (how='max')\n", "dropping workers who leave a firm then return to it (how=False)\n", "making 'i' ids contiguous\n", "making 'j' ids contiguous\n", "computing largest connected set (how=None)\n", "sorting columns\n", "resetting index\n" ] } ], "source": [ "measures = bpd.measures.Moments(measures='mean')\n", "grouping = bpd.grouping.Quantiles(n_quantiles=4)\n", "cluster_params = bpd.cluster_params(\n", " {\n", " 'measures': measures,\n", " 'grouping': grouping\n", " }\n", ")\n", "\n", "bdf = bpd.BipartiteDataFrame(\n", " i=df['i'], j=df['j'], y=df['y'], t=df['t']\n", ") \\\n", " .clean() \\\n", " .cluster(cluster_params) \\\n", " .to_extendedeventstudy(transition_col='j', periods_pre=2, periods_post=2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Creating CHK event study plot\n", "\n", "Once the data is cleaned, clustered, and in extended event study form, we can generate the event study plots." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAABdw0lEQVR4nO2dd3gc1dW437Or3rvc5N4NGLDpxTIQIBBCCRBq4oSSSkISEkjIFyD5JR8p3wf5SKMlpJMQSsAUY4OFDRhiG2xwlZuMqyTbktXLas/vjxlJK3nXKruSVtJ5n2eemblz5twzc/bumXvnzr2iqhiGYRhGtOEZaAMMwzAMIxgWoAzDMIyoxAKUYRiGEZVYgDIMwzCiEgtQhmEYRlRiAcowDMOISgZlgBKRRSJSLiLfD3G8UESO64G+00XkQxFpEJExIWRSReSp3trs6jjCbhGZJyJvicgbIrJURArc9MtFZKOINHTS8UkReVdElovI9W5ajoj8LRzbIkEf+OU77rW+JSIPiYgEkRknIo+EYfM5IlLkLu+LyOpOx+8Tka0B+zki8g8ReV1EXg1IG/D730of+OEaEXlTRJaJyEIRSQsiE1b5EJE0EXnb9cN/RORcNz1U+fCKyC9EZIl7zkw3/RkRSe6tHZEm0r4IOK/D77LTsXDLRKqIrBCRShG5ISA9QUT+6v73/FVEEtz0p1zfvSsiC9y0iJSJQRmggJuAbx/leCHQE6evB04D3jmKzK3A33ugMxjB7F6hqmeo6jzgz8DX3PRlwAnA7lZBEfEAPwfOA84Bvioi6ap6AKgWkdlh2hcukfbLs6p6iqqeAeTjXHNn7gB+1wOdHVDV11W1UFULcfzb9icrIvnA1E6nPAj8UFXPUdXzXR3Rcv9bibQfnlHVM1X1bOA94MYgMuGWjxrgbNcP1wD3u+mhysetQLGqnuf6b4Ob/hRwcxh2RJpI+yLU7zKQsMoEUA9cjvNbD2QBsElVzwI2u/sA31PV04F5wPdFJCFSZWJQBihV3R3qmIhk4dy4u90nK6+I3OJG93dF5PNB9B1W1Zousv0UsMjNY4GIPOc+ra0TkbN6a7eqNgXspgEfuOkHVbWhk3gOUK6q1araDBQDp7jHXgKu6o4dfUUf+GVLwG4j4Aui+mxVfc/N4173ye55EVkjItN7eAnXAYFPff8F/HfANXiBY4BvuU/0Xw6QHfD730of+CHwN5qM80DXmbDKh6r6VbXVv4HlIGj5wLnX49xa1a9EJM5Nfxm4sqv8+otI+8Klw+8yCGGVCVX1qer+IIfmAQvd7Rfc/cBy2gS0AK2jP4RdJmLCOTkaUdVDIvIEsFVV/yIiucBXgZNckZUi8oKqlvdQdbaq1nbK6woROR34JrBcRL4EfDrIuVeo6qFQikXkYuA+nAJ40VFsKAdyRGQ0UA2ciVMgwQlWoX7QA044fhGRecBInFplYHouzn0IpFxVrxeR63CepO8QkftwC1MATa01IFfXccBhVf3I3Z8CpKjqB9LespgHHAt8FtgIvC4iS1V1I1F+/1vprR9E5Cbgdpyn6/s5krDLh/u7/gdO7aDtXoYoH6OBfao6X0R+4cr/TlUr3RpG1NMbX4T4XRJwPGJlIgjZQIW7XQlkdTr+XeBJVW1098MuE0MuQAVhIvBh65OYiHwITMD5sw+H1ncVH+E4DlX9LfDbnipS1ReBF0XkauAnwNUh5FREbsVp6qgFPgT29tz0qKBbfnEDx/3AJdq9cbkC/fIxAFW9pxvn3QD8JWD/XuAHnWQqgL2quta1rQgnYG3shv5opVt+UNXHgcdF5Ds4TVbf6UJvj8uHqu4BzhSR8UAR7tN6iPJxCHjFPfUV4IquLzXq6Y4v7uXI32VX9LZMdOYQkOFup7v7AIjIZ3CaKq/thd6QDNUA1UT7te0AjgtoAjjWTespB0UkOeApMfDPUgB6U4Ny22tbm/IqgbqjGaGqy4BzRCQVeJb292ZTgXXduZABpEd+EZHJwO+BT7lt2h1Q1XL3PnRIDlTh6jnq06L7bu8KYG7A8YnAr92n1JEi8n+q+jUR2S4iBaq6C5gDPOPKD4b730pP/dD5N5oURGdY5QOoDXjyrsKtBRylfBTh+GtrwBoRyQBKg+QRrfT0vyro77L1YKTKRAjewKnBrnHXb7i6LsVpHv+kqvoD5MMvE6o66BbgUZx28K3Ac0GOnwS8BfwL5z3bF4H/uMstQeSnAktwnpCXA18KInMHzh8lOO3G33e3xwBFvbUbp8r9BrAU50lwnJt+lmtTnbu+wk3/mSu7GDgpQPfDwPFDzC8LXV1F7nJxEJlfAXPc7XuBG9ztM4Enumn3OcDTRzm+NWD7eNdfbwP/L5rufx/64b8CfPAMkBHp8oET7Je5v+3lwLldlI9M4LkAm5Ld9GuB2wfaB33li1C/y07pkSgTLwDbcQLM79y0RJyOMMvddYKbXgOsCviNjI5UmRBXkdEF4nSt/YOqfmqgbemMiOQAD6lqRKvXgwG3OehuVb1lAG0Ytve/lWgpHyLyLHCjdt3pacgylMqEBSjDMAwjKolIN3MR+b2IlIlI0PZGEckUkWdF5ANxPsI7JhL5GoZhGEOXSH0H9QRw4VGOfw9Yo6rHAZ8BfhmhfA3DMIwhSkQClDo9y0J+5wPMBF53ZTcB4wfLtwqGYRjGwNBf3czX4nQlXS4iJwPjcHr3HNEd1P3O51aAxMTEOQUFBQD4/X48nshU+CKmy9+AX2PweCNzG/vyGouLiw+oam53zw/lh762s/eKhp8vzA/h6bIy0UN1A+GLCHanHA+sC3EsDfgDTv/5PwMr6Ub3wzlz5mgrS5cu1UgRMV2L5+nS5x+OjC7t22sEVmkvfRvoh762s9cMQ1+YH8LTZWWiZwyEL/qlBqWqVcDnAMT5wmwHTh97wzAMwwhKvwwWKyIZAV9H3wwsc4OWYRiGYQQlIjUoEfk7zrDxOSKyG7gHiAVQ1d8BM4A/iojifFV9UyTyNQzDMIYuEQlQ2sXXwqq6gqPPX2IYhmEYHRiU80EZhmEYQx8LUIZhGEZUYgHKMAzDiEoiNRbfhSKyWUS2ishdQY4/4E43vEZEikWkMhL5GoZhGH1PQ3MLLf7+H1g87E4SIuIFfo0zU+NunGmKn1fVDa0yqvqNAPnbgBPCzdcwDMPoWz46WMfjb27nn6t209TUzLj3ipiQk8yEnGTG5yQzMSeZCbnJ5Kcm4PEcOQV9uESiF9/JOBNnbQcQkSeBS4ENIeSvxemGbhiGYUQha3ZV8uiy7by8bh9ej/DJ2aPJ//X/suOWr7HjQC1vbTtAQ3P75LmJsV7GZScxMdcNXtnJ7nYKmUmxuDMA95iw54MSkSuBC1X1Znf/RuAUVf1qENlxOFOUj1HVlhD62sa7ys/Pn/Pkk08CUFNTQ0pKSli2thIxXdXF1PizSEnPCV8XfXuN8+fPX62qc49ySgdC+aGv7ew1w9AX5ofwdFmZ6IhflbXlLby8o5niCj+JMXBOQSznjYshM8HDqEceYe+tt7bJVjQopXXK/lo/pbV+9tcppbV+yuuVloCwkhQDI5I95CcLI5I85Cd7OOXlp5DPt3+dFNIXwcY/6skCXAk8FrB/I/CrELJ34syyaOOOBcHGHQuDYegL80N4uqxMONQ3+fRv7+7U+b9YquPuXKin//dr+vjy7Vrd0Nwru5p8Lbq9vEZf31iqjy3frnc/+4Fe/+g7evp/v6bj7lyo4+5cqD+af1OHc0L5IhJNfHuAgoD9MW5aMK4BvhKBPA3DMIwwOFTbxF/e2cmfVpRwoKaJY0an8X/XnsBFx4wgxtv7/nOxXk/be6r5nY41NLdQcrCW5OndG0woEgFqJTBFRCbgBKZrgOs6C4nIdCATWBGBPA3DMIxesPNgLY+/uYN/rtpFQ7Of+dNyueXsiZw2MbvX74q6S0Ksl+kj0ii57Gjz27YTdoBSVZ+IfBVYBHiB36vqehH5IU617XlX9BrgSbc6ZxiGcVT8fqWy0U+LX/H2QQ+x4cZ7H1Xw6LLtvLJ+P7EeD5edMIqbz5rI1PzUfrelZMECxndDLlJj8b0EvNQp7Qed9u+NRF6GYQxdmnx+3tl+kEXr97N4Qyll1Y18e9nLjMpIpCAziYKsRMZkJlGQlcSYTCctJyWuz5/8Byt+hVfX7+fR5dtZWVJBWkIMX5o3iQWnjycvLWGgzeuS/ppR1zAMIyi1jT7eKC7n1fX7eW1TGdUNPhJjvRROy2Xi26/Qcsmn2FVRx+5DdSxaX8Wh2qYO5yfGep1glZVEQWZrAGsPZOmJsQN0ZQNHQ3MLT++ezUMleex/ezVjMhO555KZXD23gOT4wfO3P3gsNQxjyHCotoklG0t5df1+lm05QJPPT2ZSLBfOGsH5s0Zw1pQcEmK9FI2pprBweodzaxt97K6oZ9ehOidwtW3Xs3LHIaobfR3k0xJiGJOZRPaOzRx30ulkJccxVDlU28SfVzgdHw7Wns+ElCZ+dd0JXDgrvI4PA4UFKMMw+oXdFXW8ur6URev3s7LkEH6F0RmJXHfyWC6YNYKTxmd26080OT6GaSNSmTbiyHcnqsrh+uYOAWzXoXp2V9Sxq8ZPUpy3Ly5twNlxoJbH39zOv1bvpqHZz7nT87gl+ZfUx5/N/ONGDbR5vSasACUiWcA/gPFACXC1qlYEkRsLPIbTHV2Bi1S1JJy8DcOIblSV4tIaFq3fz6sb9rNujzOJ9tT8FL4yfzIXzBrBrFFpEX1/JCJkJMWRkRTHMaPTOxwreeM3JMR+KWJ5RQOrdzodHxZtcDo+XHHiaG4+awKT81JhyW6K6gfawvAItwZ1F/Caqt7vDhJ7F87HuJ35E/BjVV0sIimAP4iMYRiDHL/C1qo4Vry0kUXr91NysA6AE8dmcNfHp3PBrBFMyEkeENu623Ms2mnxK0s2lvLIsu2s3llBemIsXymczGdOH0deavR3fOgJ4QaoS3Gmegf4I1BEpwAlIjOBGFVdDKCqNWHmaRhGFNHk87Oitefd+1+ivCmFGM8OTpuUzc1nTeT8mfmDosdYtNPQ3MK/Vu/m8Td3sONALQVZidz3yVlcNXcMSXFD821NuFeVr6r73O39QH4QmalApYg8A0wAlgB3aYix+AzDiH5qG30UbS7n1Q37ed3teZcU56UwYw8FGbl8+dprhmXvub7gYE0jf1qxkz+/s5NDtU3MLsjgN9efyAWzRgz578O6HCxWRJYAI4Icuhv4o6pmBMhWqGpmp/OvBB7HmWLjI5x3Vi+p6uMh8rOBMftAlw2M2cXpTcruaj97avzsrfFT29hMfFwsMQJeD3hFiPGA19130gWvEJAuAfKt6UKcv55JuYN/sNiqJmVNmY/VpS2sP9iCzw+psXB8Xgwn5nuZle0lrm7LsCsTv1/XSObGdcSceCzp8UJGvJAeL6THCalx0uMg0mrn/lo/i0qaeXOPj2Y/nJDn5cLxsUzN9HTvvd0Q+H8KazRzEdkMFKrqPhEZCRSp6rROMqcCP1XVee7+jcCpqtrlmHxz587VVatWAVBUVERhYWGvbQ0kYrqWFFJUfx2Fl9wavi769hpFpEeFMZBAPwTTHQ797Yuqhma2lFZTXFrD5v3VFLvbB2oa22TSEmJIrD0MmVk0tyjNLX58reteTNp2cd12fv1/t7XtR8oX/eGH3RV1LHJ73q0K6Hl3/qx8Lpg1grnjOvW8G2ZlQlW54MFl7Nu5n+qEI/+8PQJZyfHkprpLSjx5ac66NS3PXafExyAiPPrsa6ysTmfxxlJivR4+deIYbjpzApPzehgchoAvwm3iex74LHC/u/53EJmVQIaI5KpqOXAOsCqInGFEjLomH1vLnCC0paw9GO073NAmkxznZUp+KudMz2VqfipT852uy3mp8ez83OcY//MnjtCrqvj8TrBqblF87ro1eDnp7QGtuUVp+skrEb++F9bu5TfPf8SfS1aSnRJHTko82Snx5Ljbzn4cmUlxPXqCD+x5t2j9ftbvdXreTctP5avzJ3N+H/S8G8yICK9+Yx4lCxYw4tHHKa9upLymkbIqZ11e3bo0UF7dyNbSasprGmluOfJBJyHWQ1pCLGXVjWQktXDb/MnceNp4clPjB+DKooNwA9T9wD9F5CZgJ3A1gIjMBb6oqjeraouI3AG8Js6vejXwaJj5GgYAjb4WtpfXUrx/Bosr0/jLH1dRXFrNroo6WhsH4mI8TMlL4dSJ2W4QSmFKXiqjMxJDzgIaqseXiBDrFWJ78NFj0Y2f6vmFdUF8jIfc3SXsGzeBdXsPc7CmKWjtznmCj2sLWNnJ7cEr113npMSz+VALb7s973YerEMETijI4Ltuz7vxA9TzbrBQsmAB42O9zmgWWUlHlVVVKuuaOwSwMjeAHaxp4viffp8r33luyHZ86Alh3QFVPQicGyR9FXBzwP5i4Lhw8jKGN74WPyUH6ygurXZrRc665GAdLX4FPoFXlEl5tRw3Jp0r54xxa0UpjMtOHnIvk8+fNYKpSdsZ/3VnyEu/X6lqaOZATRMHapw/OmfdSHlNEwdrGjlQ08jaikoOVDdS23RkH6VY7w5Om5TDrWdP5GMzrOddXyEiZCbHkZkcF3Sg1pKXcyw4udhdMHrE5b95i5ptu/n99v+Q4z6F56TEk5Pa3rSUkxJPVnLPmpZa8fuVXRV1FJfWuO+HnEC0vbyWphbn8zmPwLjsZKbmp3DRsSOdWtGO2/hILuC8SyPT3j4YCKzleTztH6h2511FfVOLE8BqmzhQ3ci+X/6SS395H2kJ1vNuoBkq32tFgkERoBqaW2juxctpI/IcMyqdLcX1VNY1sbW0mgM1TW2BI5DWpqXs5CODV05KHDmp8eQkx/NBuY/iZdvYvL+GLWXVbCmtob65/el+dEYi00akUjgtj6n5KUzNT2VyXgoJsZ2GrCk/yN5B/tV8f5IY17E5quja8y04GVHHoAhQL6/bxzderSPn7cWMTE9kZHoCozKc9ciMREa56/zU+EE5IGJPUFWq6n1tbdaBL2TLqhra9vdW1PL6iQ3kR7iZ5keXHUNRxgEKC89st6fB16FZ6UBNIweqnaal1mam9z+q5EBNI3VBmpZgE/lp8UzNT+W6U8a2BaIp+amkDKKRlw3DiCyDovRPH5HG56vXUX/yx9lb2UDJwVpWbDt4xKjFHoG81ARGZiQwKr1jABvhBrXclL7pEaOqqDoDDfrdbb/7lt4f5JgGyLSmH6z3s3ZXJWXVR748bUuraaTJd2SNJT7G09Z9dVJuCgUJDfRHRysRIT0xlvTEWCbldi1f1+TjQHUT5W7ganjwF5z9y5+QkTR0R5g2DKN3hDtY7FXAvcAM4GS3c0Qwua8DtwACPKqqD/Yknxkj0zj7qlMoLOzYz6K6oZl9hxvYW1nPvsMN7KusZ+/hBvYdrmfjvipe21RKQ3PHP/MYj5ASC4krXgsIJACKvzVw4LwLUWgLJn4FdWVoDSz6Lfwq6Fsv9uRyjs4bb3XYzUqOa/t2YmJOMrlprd9SJLR/U5EaT6r7DUUrRUVFUTkuV1JcDGOzYxib7TYtXf8JC06GYQQl3BrUOuAK4OFQAiJyDE5wOhloAl4RkYWqujXMvElNiCU1ITbklMWt3Tn3Hq5nX6UTuPYebmDdlp2MHJmDIHg8AIJHQAQ8IohjNyI4MgHHCJQp+Qsf+WYxYeqJbfKtx1q7L3dIa9XZth+QjwiH//gE02/7ovPxXlo82cnxxMUM7SZLwzCMUITbzXwj0NVHezOAd1W1zpV9Ayeo/SycvLtDYHfOWaPah94vSthPYeHs8DNY8iZF9WMpPH9a17LdoKhhHoUzgw1naBiGMfwIa6ijNiUiRcAdwZr4RGQGzggTpwH1wGvAKlW9rbOsK9823hUwDdjsbucAB8I2dvjqGqeq3XhL5HAUPwTTHQ7DUVekfDFYrjdadVmZiB5dwX3hvKwPveCMPr4uyHJpgEwRMPcoOm7CGUFiGfBb4MGu8g2iY1VPzzFdkdE1WOwcDroGg43DQddgsXOw6+qyiU9Vz+tKphs6HscZ0RwR+QmwO1ydhmEYxtCmX7qZi0ieqpa5U79fAZzaH/kahmEYg5ewuoiJyOUishvn/dKLIrLITR8lIi8FiD4tIhuAF4CvqGplL7J7JBxbTVdEiVY7h4OuvtJruqJHt+lyiUgnCcMwDMOINPaRjWEYhhGVWIAyDMMwohILUIZhGEZUYgHKMAzDiEosQBmGYRhRiQUowzAMIyqxAGUYhmFEJRagDMMwjKjEApRhGIYRlViAMgzDMKISC1CGYRhGVGIByjAMw4hKBnWAEpETROQtEVkmIq+LyMQgMpe503x0V+flIrJRRBqOIjNORHo9sq+InCMiRe7yvoisdtOvEZE33etZKCJpbvp49/reEpHvuWlxIvK8iESFD/vIFw+KyDvuclcImVNF5L4w7L4uwBcbReRpN/1ed7/1mNdNfzUgrV5EjhWRHBH5W29tiCR95IdvuvreEpE/iUhsEJmwyoSr4wIReU1ElorIN92000XkQxFpEJExAbLDskwEnPdHEVkS4li4ZWKea/cbri8K3HSviPxCRJa4v/+ZAefEisgWEfm+uz9TRB7srQ1t9NUMk/2xACOAVHf7IuDPQWSeAM7sgc5sIAHYehSZh4ATI3QN3wHucrfjAtJ/iDM1CcCTwFnu9hJgurv9XQJmNh6Cvpjirj3ACmBSEJlngbwIXcNvgGvc7XuBG44iOxJYF7D/MDB7iPoh8Hf5J+DiIDJhlQmcKcAXBublpqcDKTizdo8JSB+WZcI95zjg38CSEMfDKhOd/P154Ofu9peAW0Oc81Wc6ZS+H5D2ApAbzj2MiieN3qKq+1W12t1tBHyBx90IfyHwkIg85abdIyIrRORdEbk4iM6Dqhqy9uRytqq+5+q7V0T+6j65rRGR6T28jOuAv7l5NwWkJwPr3e3jVXW5u/0iMM/dfgm4qof59Ql95Ist7trv6mvppDMVpwCUuftPiMijIvKiW+vK6679bq3g4zgFv5XvuDXarwU55VqcP8lWosIXfeSHJldOcB4WtgbJOtwycTFwCHheRF4RkVlu3odVtSaI/LAsEy7/Bfwk2IFIlIlO/0NpwAfu9lXAOLdW9SsRiXPzSMEpO093UrUEuLSr/I5Gv8yo29eISDLw/4CbAtNVdYOIvAI8pqpvisjxwFnA6ThPZv8RkZfdP8Du5pULVHdKLlfV60XkOuBm4A63ij2vk1yTqp4foOs44LCqfhSQdhNwO1AP3O8mBz5IVOI8mQEUA8d21/b+oC98ISLXA9tVtaTToelA57T1qnqL2+xzNfArEXkYmNZJ7iNV/UzA/seBZapa7+4/BNyHU5t+QUTeD/hDBLgeuDJgvxjnaTMqiLQfRORuYAGwBdjV6VgkysQoYDJQCMzAmdDujKNc4rAsEyJSiHONpSGyi0iZcIPjfTgB6iL3+Ghgn6rOF5Ff4Pzefwd8G3jQPR5IMU4A7jWDPkC5T77/AH6qqhu6EJ8GvKNO/bNSRMpwmhbKwjRjtbv+CPgYgKre043zbgD+Epigqo8Dj4vId3Ac/x0g8M8iHedJM+roC1+IyHnA54BLumlGoC8mAajqF7px3g3Ao607qnrQ3awXkWeAucBy16YZQL2q7uimTf1KX/hBVX8sIj8BfoUTqH7Thd6elolDwOvu0/vabjzpD9cycRdwDZDRAzN6XCZU9UWcWdKvxqmtXY1zj19xRV4BrhCRfOAEVb1HRBb0wKZuMagDlDgvQ/8CPKeqz4UQa6L9OouBW9yminQgDzjQkzxVtdytRndIDjTLte2oNSjX9itw/vharychoHmxEkhyt9eKyOmq+jbOk/7tbvpUYF1P7O8r+sIXInIK8CPg4wE1m0A2AeM7pQXzRVdPi2nAHJyC35p3hqpWuvYV4rwraOVG4K+d9EWFL/rIDwmq2qCqKiKHgbrA4xEqE0XAA65cAXC4i0sddmXCvccjcJqWE4FZInK3qv44QF/YZSLI/1Crv4tw/q+2BqyPBXLdmuBoIF5E1qrqC0TAF4M6QOH8wV8M5IvIDcCHqnpbJ5mFwA9FZKOqfkFE3sZ54e4BvhWkKeMs4B5glDi9ZH6jqs900rlcROao6mpC0I0aVCGwVlUrA9K+LSLnutuHaG8y+i5OrSoOeFlVN7rpFwFPdZFPfxFxXwCPu+vnnDLLtwLvuapWi8gBEclX1VBNHt2pQV2J8ycSmP+DIjINp0AXqepL0PYe5krgtE46LgJ+20U+/UFf+OF/3HdCre+fgv22wyoTqrpZnJ5hy4BY4OsAIjIVp7Y2G/i7iPxNVX/LMCwT7vus48HpxYjTNBgYnCJVJm4QkRtxaqmNQKvsz4A/iMgXcf6fblTVWpx3Tbg1qDFucAI4j07Nmj1FnNqk0RPcH8fdqnrLANoQB/wLuKwn79CGGiJyGk4N6wcDaEMO8JCqXjtQNgw0ViaihygpEzOBL6jq18PSYwHKMAzDiEYi0s1cRH4vImUiErS9UUQyReRZEflARP4jIsdEIl/DMAxj6BKp76Ce4OjdCb8HrFHV44DPAL+MUL6GYRjGECUiAUpVl3H0bp4zgddd2U3AeLd7omEYhmEEpb968a3F6dGyXEROBsYBYwjysZmI3ArcCpCYmDinoKAAAL/fj8cTmQpfxHT5G/BrDB5vZG5jX15jcXHxAVXN7e75ofzQ13b2XtHw84X5ITxdViZ6qG4gfKGRG3dqPAFjk3U6lgb8AVgD/BlYiTNUyVF1zpkzR1tZunSpRoqI6Vo8T5c+/3BkdGnfXiOwSnvp20A/9LWdvWYY+sL8EJ4uKxM9YyB80S81KFWtwhkNoPU7kh3A9v7I2zAMwxic9MtgsSKS4X6jAM64XMvcoGUYhmEYQYlIDUpE/o4zMkKOiOzG+dI8FkBVf4cz+OMfRURxRugO6+tiwzAMY+gTkQClXXxBr6orcMZlMgzDMIxuMdjH4jMMwzD6kOqGZlZsO8iHpT7SP6ogPy2B3NR4Yr19/4bIApRhGIbRgT2V9by2sZTFG0p5Z/tBmlucIfEeev9tAEQgOzme/LR4RqQlkJeWQH5aPPlpCe6+s52VFIfHI722I1LvoC7EGR3CizPC7v2djj8AzHd3k3CmI86IRN6GYRhGePj9yod7DjtBaWMZG/c5fdgm5iTzuTMmcM70PKp/8iO8d3yX/YcbKa1qoKy6gf2HG9h3uIG1uys5UNN0hN5Yr5CX6gas1ARGpDvbMRUtFHbDrrADlIh4gV/jTEq2G1gpIs9rwORcqvqNAPnbgBPCzdcwDMPoPQ3NLby19QBLNpby2sYyyqob8QjMHZ/F9y6azrkz8pmUm9ImX3TDZRRODz0AUJPPT3mNG7yqnOBVWt2638jW8hre2naA6gYf19ZsgsvPDamrlUjUoE4GtqrqdgAReRJnHvpQs0deS/D5ZAzDMIw+pLy6kdc3lbJ4Qxlvbi2nodlPSnwM86bmcu6MPOZPyyMzOa5rRUGIi/EwOiOR0RmJR5Wra/Kx/QvPdktnJALUaGBXwP5u4JRggiIyDpiAOy6fYRiG0Xeowub91Sxx3yet3V2JKozOSOTTcws4d0Y+p0zMIj7G2282JcXFcOCzN3RLNuz5oETkSuBCVb3Z3b8ROEVVvxpE9k6cGRc7zyoZKNM23lV+fv6cJ598EoCamhpSUlJCndYjIqarupgafxYp6Tnh66Jvr3H+/PmrVXXuUU7pQCg/9LWdvWYY+sL8EJ6uoVomfH6luMLP+7sP8N6hJA42Or3tJqR7OCHPywl5MYxJEdxZqvvNrqPpCumLYOMf9WTBmfp6UcD+d4HvhpB9Hzi9u7pt3LHI6cLGHesRg8EX5ofwdA2lMlFZ26TPvb9bv/LX1XrMD17RcXcu1KnffVY/ef9T+rd3d2rp4foBsau7ukL5IhJNfCuBKSIyAdgDXANc11lIRKYDmcCKCORpGMYQ5VBtE9vKa9heXsM7W5so9mwjMS6G5DgvSXFekuJiSI73khjrruO8JMfFkBjrDatL82Cj5EAtSzaWsmRjKStLKmjxKzkp8Vx07EjOm5nPmbuu592mqyk8eexAm9prwg5QquoTka8Ci3C6mf9eVdeLyA9xouLzrug1wJNutDQMYxjja/Gzq6KebWU1bD9Qw7ayWraV17CtvIaKuuYOss9u3dRtvYmx3o5By10nxXmprmjghFOaSU+MjfTl9AstfuX9jypY7Pa621pWA8D0Eal8cd5EzpuRz+wxGe1Beq9vAK2NDJEa6ugl4KVOaT/otH9vJPIyDGPwUNXQzPbyWraV1bi1IicQlRysbfv4EyAnJY6JuSlceMxIJuUmMyk3hUm5KTR96yuMePhR6pp81DW2UNfU4my769rGFuqaW6hr9HU61nH7QE0jVbsP08PXLgNObaOP5VvKWbyhjKWbyzhU20SMRzh1YjY3nDKWc2fkU5CVNNBm9hk2koRhGGHh9yt7D9ezrbyWbR+dyPKqdB5+5B22lddQVt3YJhfjEcZmJzEpN8X9xiaZSXkpTMpJIT0peK2maMFnmRwfQ0p8DKSGZ2fJggWkJXw6PCVBePiNbWzb0cTOuBLiYzzEx3pIiPESH+shPsZLgruOj/GQEOt1Zdx1jOeIzgqHGvz8+Z2dLNlQyoptB2lq8ZOeGMv8abmcNzOfs6fmkpYwOGuBPcUClGEY3aK+qcVpjguoEW0rr2XHgRoamv2u1Lkkef1MH93C2VNz3ZqQE4jGZiX1y/htoShZsIDxfaD3fxcX0+jzw+b1vTq/NVDFx3qJ9Qh7DzcA6xifncRnThvHeTPzmTsuk5gBvHcDhQUowxiEvLX1AL9f18ibNRtISXBqGCnxMSG3k+NjuhUcVJWy6sa24BPYNLensr5NTgQKMpOYlJvMGZOymdgaiDbdyIe+y5j/yUv68vKjivX3XcCWm29hxK9+S4OvhcZmP40+Pw3NLR3Wjb4WGpqddWOzP6hso6+F3Dde5qp7bmNSbnKPu4IPNSxAGcYgZHdFHet2VrOy9CPqm1u6dU5CrKdD8EqOiyHVDWKlpQ08sO5NtpXXUtPY/nI9Oc7LxNwUThqfyTW5BU6TXG4K47KTSIgN8nHn9jqke+YMGWK8Hko/eyMzejkCQ2eK8iqZnBeZ740GO2EFKBHJAv4BjAdKgKtVtSKI3FjgMaAAUOAiVS0JJ2/DiDZUlbqmFg7VNnVYDhzq3sCYPeHTJ43llF//gPFPPIGvxU9tUwu1jT5qGn1UN/jatmsafdQ0OOvaRh/V7n7r9r7DDdQ0+mgsrWLKjEw+deLotiA0KTeF/LT4Yf8Ubwwc4dag7gJeU9X7ReQud//OIHJ/An6sqotFJAXwB5ExBgHffmotJXsaeL50TdtL4MA29LbttpfEATId5L3ExXg43KgcrGlERPAICIJ4QACPCNKaJk6zkkek47E+/PP0tfipqGumoq6JgzVNzrq2iYqA4NP5WJPvyJ/2x+t2wBVdD4zZU1rfqcR4PaQnesLqPl2yYAHjb3oiUqYZRkQIN0BdCm0Ph38EiugUoERkJhCjqosBVLWmp5nUN7VQ06T4WvzD8kVhNLG7op7S3RXsa4xz2szb2s7DeOZYuiQsm0TAw7cAwfvOyyC0BTuPG8ScQEd7IHTXtMm0Bz8RoerQYapfeTlknqnxMWSlxJGZFMfI9ARmjkojOzmOzOQ4spLjyEqKIyvFWdd9u3sDYw4kfdWBwDDCIayx+ESkUt15ncR5lK3QTvM8ichlwM1AE85AsUuAu1Q1aEt1sPGu3trTzKMfOnONxHshKUZIjHXXMUJSDCTGirsPSbHt6YHbiTFCQgzU1dbauGNdcLRxx0Y98gh7b721g7yq4lPw+aG5BZr9SrMfZ2lp3XbX7vGmFkh9800qzzgTBdD2qrVf2/UqzqCXCh23A9MaD9GkCcTEJQWcr8HPC6ZDO+abtXkDnuNnkRonHZaUWEiNE2J6MGKBjcXXewZLmTBfhKcrlC+6rEGJyBJgRJBDdwfuqKqKSLBoFwOchTMH1Ec476wWAI8Hy09VHwEeAZg7d64WFhYyurSauLf+j9grrqK6wUdVfTPVDT6qG511WYOPqmpnu6nl6E/yIpDa0sSEcbFMzElmQqclOb4Hlcol91JUfx2FhVd2/5yjUFRURGFhYVToCuaHNt0QOTvHxkZG15JCxxeX3Nq1bDcoKoqPel9E0++ljeFaJqLIzjaGgC+6/DdW1fNCHRORUhEZqar7RGQkUBZEbDewRtvni3oOOJUQASoYU/JTOenqMygsnNqlbENzixO8Gpqpctdt+/XOumThIg7Fj+Y/Ow7x7Pt7OpyfnxbP+OxkJua2Bq0UJuQkMzYribgYa140DMPoL8J9B/U88Fngfnf97yAyK4EMEclV1XLgHGBVmPmGJCHWS0Ksl9zU+JAyRXH7KCx0pqxqaG6h5GAtO8pr2X6glpIDtew4UMur60s5WNs+hbFHoCArqa2mNTEnmfEHx3HA48Xv12E1SKVhGEZ/EG6Auh/4p4jcBOwErgYQkbnAF1X1ZlVtEZE7gNfc91SrgUfDzDdiJMR6mT4ijekj0o44driumR0HnS/lWwPYjgO1/GfHIeqaWnAvl7vXvsL4bLfGldsewCbkJJOVHNdv3XRVlepGH2VVDWw42MLpPr/V+gzDGLSEFaBU9SBwRP9ZVV2F0zGidX8xcFw4eQ0E6UmxHJ+UwfEFGR3SVZXy6ka2L7qZVyovJG70HLaX17KlrJrXNpV2GAQzNSEm4F1XChNy3dpXTrIzvlg3qWvyUVrVSGlVA6VVDZS1blc3uvsNlFY1dvho8xOFDYzNHroDSRqGMbSxkSR6gYiQl5ZAXuZuGhLqKLxoRtsxX4ufPZX1Tm2rvNZpPjxQy8qSCv69di+BnSbzUuOd2lZuMuOzkzm018f2N3dQWh0QgNxgVN145ND5ibFe8tPiyUtL4NgxGZyXGk9+WgJ5afHorx8gN/XC/rgdhmEYfYIFqAgT4/UwLjuZcdnJzJ/W8VhDcws7D9ax40BNWwA74n3XBxuIi/GQnxZPfmoC00ekcfZUJ/Dkp8WTl5rQFpRS42NCNh8WXf8JEuOCDEVjGIYxSLAA1Y8kxHqZNiKVaSOOnDfgcH0za776TWY/9D+kJ8ba8DKGYQx7LEBFCemJseiNnyIjKTIDThqGYQx2wuriJSJXich6EfG7PfdCyX1dRNa5sreHk6dhGIYxPAi3D/I64ApgWSgBETkGuAU4GZgNfEJEJoeZr2EYhjHECStAqepGVd3chdgM4F1VrVNVH/AGTlAzDMMwjJCENVhsmxKRIuAO9/unzsdm4IwwcRpQD7wGrFLV20LoahuQEZgGtAbAHOBA2MYOX13jVDW3uycfxQ/BdIfDcNQVKV8MluuNVl1WJqJHV3BfqOpRF5zRx9cFWS4NkCkC5h5Fx004I0gsA34LPNhVvkF0rOrpOaYrMroGi53DQddgsHE46Bosdg52XWENFttdVPVx3MFhReQnOAPIGoZhGEZI+qWbuYjkqWqZO/X7FTijmRuGYRhGSMLtZn65iOzGeb/0oogsctNHichLAaJPi8gG4AXgK6pa2YvsHgnHVtMVUaLVzuGgq6/0mq7o0W26XCLSScIwDMMwIo3NxWAYhmFEJRagDMMwjKjEApRhGIYRlViAMgzDMKISC1CGYRhGVGIByjAMw4hKLEAZhmEYUYkFKMMwDCMqsQBlGIZhRCUWoAzDMIyoxAKUYRiGEZUM6gAlImki8raIFInIf0Tk3CAyhSJyXA90ni4iH4pIg4iMCSGTKiJP9YXdIvIZEXlNRJaKyHVu2oMi8o673OWmxYnI8yISFT7sI198R0TeFZG3ROQhEZEgMuNEJKxBLEXkgoB7/k03LejvQETGi8jrrk3fc9Oixhd95IdrRORNEVkmIgtFJC2ITLhlYqZrc5GIrBCRgwHHgpWJLNeW5YG/DRF5RkSSe2tHpOgLPwScd5+IbA1xLKzyICLnBPjhfRFZ7aaHKg9fFpHiQHsiWh76agKv/lhwAmyMuz0RWBlE5l7ghh7oTAdScCZhHBNC5lvAFZG2G5gF/Al3EN8A+SkB560AJrn73yVg4sgh6IspAdv/BM4NIvMQcGIYducAC4G47vwOgCeBs9ztJcD0aPJFH/khLmD7hzgzEkS0THTSdTXwO3c7VJm4H7jR3f49cKG7fS3w9aHoB/ecfODvwNYQx8MqD510fQe4y90OVR7ygdjO9kSqPAz4E184qKpfVX3ubhrwQeBxEckCFgB3u08EXhG5xX0qf1dEPh9E52FVreki608BrVOLLBCR59wnt3UiclYYdl8J1AGvisizrU8qqrql9TzAB7S48i8BV3WVX3/QR77YErDbiHPtnTlbVd9z87hXRP7qPr2tEZHp3TD9YuAQ8LyIvCIis9y8Q/0OjlfV5e72i8A8dzsqfNFHfmgK2E0G1gfJOqwy0YkbgL+420HLBM59X+huv0C7H152zxlQ+sIPLv8F/PdRsg63PARyHfA393qClgdVLVXV5iDnRqQ89MuEhX2JiIwG/gFMBTo4VVUPicgTONH9LyKSC3wVOMkVWSkiL6hqeQ+zzVbV2k55XSEipwPfBJaLyJeATwc59wrXrmB2jwKygPNx/jh/AVwTcK3XA9tVtcRNKgaO7aHtfUZf+UJE5gEjgWWd0nOB6k7i5ap6vdsUdDNwh4jcR/sfWCtNqno+zj2fDBQCM3DmqTnjKJcZ+FBXCYxwt6PGF33hBxG5CbgdqMepvXQm7DLh5pMNTAfeco+FKhNZOPcfd53l5lkpIvlBb0w/E2k/iMgUIEVVP5AjW7sjVR5adR0HHFbVj3p84Q4RKQ+DPkCp6h7gTBEZj1P9XHgU8YnAh61PhCLyITAB6GmA6sxqd/0RkO3a9Vvgtz20+xBOU4CKM/njT1vlReQ84HPAJWHa2mf0hS/cgnI/cIm6bQddEOiLj7l23XMU+UPA664da0Ukrwv9/oDtdPf8qKIv/KCqjwOPi8h3gG/jNP8cjR6XCZdPA08F+DpUmajAuf+VDB8/3Av8oIdm9LQ8tBJYix0wBnWAEpF4VW10d6s48ukBoIn269wBHCcice7+sW5aTzkoIskBT4yBf5ytL2tDPi0CtSHsLgIuBx4D5gDbXF2nAD8CPq6q9QG6pgLremF/xOkLX4jIZJz3C59S1QOdlalquYikdk4OVOHqOdoTYxHwgCtXABwOdY0ua0XkdFV9G/g4Tq0CosQXfeSHBFVtcHcrgaQgOsMqE601KOB6nCf9VooIUiaAN4CLcJqgLgKecfPIAEqD5NGv9NF/00Tg127taaSI/J+qfq31YITKA+J0brgCmHv0qzwqkSkP4b7EGsgF5we7DFgKLCf4S/STcJoL/oXTPPNF4D/ucksQ+ak4L78rXJ1fCiJzB86fJjjtyN93t8cARb21G+cH9ABOoVxG+wv4de5S5C5ztP1F5GUD7Yc+9MVCYGvAdV8cROZXAffjXtyXzsCZwBPdtP07ru0rgDOO9jvA+ZNY6l7H9wN0RIUv+sgP/xXgg2eAjEiXiYB7u6pTWqgykY3zDnA58GvA46ZfC9w+FP3Q6dxQnSQiUR7OAZ7ulBaqPFzlpte569MjWR5syvdeIE432z+o6qcG0IY4nB/2Zep0nhiWuM0nd6vqLQNow7D3RTSUCdeOZ3F693XV0WlIMtTKgwUowzAMIyqJSDdzEfm9iJSJSNA2RxHJdLuIfiDOR2vHRCJfwzAMY+gSqe+gngAuPMrx7wFrVPU44DPALyOUr2EYhjFEiUiAUtVlHL2b50zgdVd2EzA+Wr5VMAzDMKKT/hpJYi1Ot0VE5GRgHE7vHsMwDMMISn99B3U/8EsRWQN8CLxP+3A9HRCRW4FbARITE+cUFBQA4Pf78XgiE08jpsvfgF9j8Hgjcxv78hqLi4sPqGpud88P5Ye+trP3ioafL8wP4emyMtFDdQPhiwj2+x8PrOuGnAAlQFpXsnPmzNFWli5dqpEiYroWz9Olzz8cGV3at9dIp+9LerIE+qGv7ew1w9AX5ofwdFmZ6BkD4Yt+qUG5X3fXqTOMx83AMlWt6o+8DcMwjN5xuK6ZjfurKC6tZsfuZmK3HmB0RiIjMxKIj/H2ef4RCVAi8necwTZzRGQ3cA/OEOyo6u9wBuH8o4gozkjIN0UiX8MwDCN8WvxKycFaNu2rZuO+Kjbuq2LT/mr2VNZ3kPvDuncBEIG81HjGZCYxOiORMZmJjM5MdLeTGJOZSEJs+AEsIgFKVa/t4vgKnKEyDMMwjAHkcH0zm9wA1BqMNpdW09DsDPrg9QiTcpOZMy6TG04dx4yRqUwbkUrJN78F3/sReyrr2V1Rx56KenZX1LNmVyUvfbgPn7/joA85KXFtAWt0phvEMpxA1uDr3gARg3qwWMMwDCM4LSrsr4/hxQ/2sWl/lRuMOtaKMpNimTEyjetOdgLRjJFpTM5LCVr72Xzj1RROyg6el18pq25gd0U9eyrq24LY7op6Nu6rYvHGUpp87aMeXVG7hQvPm9/lNViAMgzDGORUNTS3Nc9t2l/Fhn3VFO/5GvX+OHjvPbweYWKOUyu6/tSxzBiZxsyRaeSlxhNsbqme4vUII9MTGZmeyEnjjzzu9ysHahvbal2xDy3pll4LUIZhGIMEv1/ZeajOCUT7nEC0cV9Vh1pRRlIsM0akce3oD5D4SVz+8UtC1or6C49HyEtNIC81gRPGZlJ0w6XdOi9SnSQuxBm+yAs8pqr3dzr+ANBan0sC8lQ1IxJ5G4ZhDEWqG5oD3hM56837q6lvdj4hba0VndhaKxqRxoyRaeSnubWiJXdRVD+SY0anD/CV9J6wA5SIeHHmY/kYsBtnquLnVXVDq4yqfiNA/jbghHDzNQzDiAaWFZezpsxH0/r9tPiVFlVn7Vd8fsXfulbF1+Ku/e0ybYsqzT4/q4sb+P67r7O74sha0TUnFzBjZBozRqQxJX9ga0X9QSRqUCfjTJ61HUBEngQuBTaEkL8Wpxu6YRjGoOcLf17t1GreW921cAi8HnEWEfLqKzj+xGlce/JYZo5MY/rIVEakJUTkXdFgI+z5oETkSuBCVb3Z3b8ROEVVvxpEdhzwDjBGVbsc6ig/P3/Ok08+CUBNTQ0pKSlh2dpKxHRVF1PjzyIlPSd8XfTtNc6fP3+1qnZ7CudQfuhrO3vNMPRFODaqKk0tUOdT6nxQWV1HclIiHgGPiLtuX7zifPviEcEDeD3OQJ7SQU6GpR+2V7aQ+9yzHLri8rZ75+10/1rvT+d0rzhD6wQGn1GPPMLeW28N/4KHgC/6u5PENcC/QgUnAFV9BHgEYO7cuVpYWAhAUVERrdvhEjFdS+6lqP46CguvDF8X0XWNofwQCd2BmC+6Jpgvmnx+Xn7tDcbOmkt1g4+qhmaqG3xUu+uq+maqGnxtae3H22VbOny3IkBDWNcpAl5GAB5iYxqI8Qget2bgESGmddsDXpG2WoMnYLu1FuFx11WH68nOTkJo/5OXtgDq6BIR97iElAFh/75GfnXzmSTH9+5vL1SZKASKMrxcEqnfC1iZcIlEgNoDFATsj3HTgnEN8JUI5GkYw5o/rSjh/y2tg6VvhJRJjY8hLTGW1IQYUhNiyE9LYEpeDKkJTlr7sVhqfvMQOV/7On5VWvy471H8tPidnmMt2v4+pcXf8T1K63G/X/Ft/wslzbMYPfE4R4/f7+qiw7uYtnc1Le3nBr67cd7f+GnwQVV9M6qKX0FR/H7wq6LqrAO3FTfNT9s5fnfdcrgGX4vNID6YiESAWglMEZEJOIHpGuC6zkIiMh3IBFZEIE/DGNacNimbr1e9x/hbPkdqfGyHQJSWGEtKXAweT/ffWRRdewGFs0aEb9iS5RTVF1B48czwddH6pH1GRHSVLFhAetKnIqLL6B/CDlCq6hORrwKLcLqZ/15V14vID3FGqH3eFb0GeFLDfellRB2qSl1TC1UNzVTV+9x1e7NSa3NTa1qgTNz+Pfxj92oKspIoyExkTFYSBREcy2uoMmtUOuVXn0XhCTatWncpWbCA8QNthNEjIjUW30vAS53SftBp/95I5GUMLN/85xo2ltTzPx++GRCIOr/POJLEWC9piU5zUlpCDFnJcYzLTqZ07zY2l1bz2qayDkOhgDMYZWvgKshyglZBZhIFWUmMTE8gxttf820ahjEQDLmRJFSVRp+fmkYftY0+d91CbaOPajdt865mxh+oZXxO8kCbO+gorWrAu7+MnGOzmZibTFpCLGmJMe46Nsi+E5TiYoIHk6KRhyksLMTvV8prGtl1qI5dFXXsOlTftr2ypILn1+4lMAY6Q6skuAErkYJDp1LtTSR15yEKMpPISYnvUROXYRjRx6AIUJv2V/Hi9iZWNW6mpi3otK9rG1s6pHceVTcYT6wvYmJOMvOn53HO9DxOGp8V8k/UaOevN59KyYLfMf5zn42oXo9HyE9LID8tgbnjs4443tziZ//hho4BrKKOXYfqWLq5nPLqswB4ZIvzijM+xsPothpXe82rdT89MXZYfldiGIOJQRGg1u2p4qniZjxbtpISH0NKfAzJ7pISH0NeaoK77W1LT02IITmuXSY53uucmxDD1m9+m22f/xavby7nz+/s5PE3d5Ac5+XMKTmcMz2P+dPyyEtLGOjLjloGoi0/1utxAkxWUtDjDYvO45nK6xg5+xNtgas1iK3ZVcnh+uYO8qnxMe77Lmc6gA5BLCvR+U7I56e5xY+vRWn2u+sWN83vHPf5FV+Ln2b3mM8fsO2uKypaKOyHe2QYQ41BEaAumT2SYx66h2l/eCwiT72bb7ySBWdMYMEZE6hr8vHW1oO8vqmMos1lLFpfCsAxo9M4Z1oe86fnMXtMhjUXRTkJXh+jknwUTs8LevxwfTO73ZrX7tYAVlHPjgO1LNtS3jYXTgcWvRwR2y6q2w6XnxsRXYYxnAgrQIlIFvAPYDxQAlytqhVB5MYCj+F8L6XARapa0t184mO87F9wI9P7oEkmKS6Gj83M52Mz81FVNu2v5vVNZSzdVMavlm7l/17fSnZyHPOm5XLO9DzOmpJLemJsxO0w+pb0xFjSE9OZNerIgTNVlQM1TW01r90V9ZT/6xlyr7yCGI8Q4/UQ53XWMR4hLsZDjMdDjFeI8zrrGI+HuBjplO7IV3zruf6/YMMYAoRbg7oLeE1V7xeRu9z9O4PI/Qn4saouFpEUIMjj6sAjIs5AjCPT+Mr8yVTUNrFsSzmvbyrj9U1lPPPeHrweYe64TM6Znsc5NdmovbYa9IgIuanx5KbGc+LYTACK5GQKCydHRP+2z14fET2GMdwIN0BdCm3N63/EGaWjQ4ASkZlAjKouBlDVmjDz7Dcyk+O49PjRXHr8aFr8yvsfVbQFq/9+eRP/zefJifdxka5j/vQ8TpuYbd/uGIZhRIiwBosVkcrWeZ3EeTlU0XmeJxG5DLgZaAImAEuAuwb7YLEH6/18sHsv7x1KY3OVl6YWiPPAjGwvs3OdJTuxZ9WrwTAwZl/b2WuGwMCYoRgsZQIYln7oazt7zVDwhaoedcEJKOuCLJcClZ1kK4KcfyVwGJiIU2N7Gripq3xVlTlz5mgrS5cu1UgRMV2L5+nS5x/W+iafFm0u03v+vU7P+unrOu7OhTruzoV6/v++of/90kZ9d/tBbfa19J9dQXThjOrR5T0PtgT6oa/t7DWuLyLFYPCF+SE8XVYmesZA+KLLJj5VPS/UMREpFZGRqrpPREYCZUHEdgNrtH2+qOeAU4HHu8p7sJAQ62Xe1FzmTc3lnktmsq28lqVuU+Bjy7fzuze2kZ4Yy9lTczlnei7zpuaRlRzXpza1Oti+9TEMY7AS7juo54HPAve7638HkVkJZIhIrqqWA+cAq8LMN2oRESbnpTA5L4Vbzp5IVUMzb2450NaN/YW1exGBEwoynG+upucxc2RatwJJk8/PodomDtQ0cqCmkYM1TRysddYHapz0tv3KWpaf2MiIdPueyzCMwUm4Aep+4J8ichOwE7gaQETmAl9U1ZtVtUVE7gBec99TrQYeDTPfQUNaQiwXHTuSi44did+vfLjncFuw+sWrxfzi1WJGpCUwf3ouKfXN7FpRwoEOgad9XdXgC5pHXIyH3JR4clLiyEtNYMaINGLeWEOs94J+vlrDMIzIEVaAUtWDwBFfIKrqKpyOEa37i4HjwslrKODxCLMLMphdkME3PjaV8upGijaXsXRzGQvX7qO60QcfrkcEMpPiyE6OIyclnpmj0shJiSc7OY5sNxAFrpPjvEfUwIpyK8hOiR+gKzUMwwifQTGSxFAlNzWeq+YWcNXcAppb/Lx969eZ+eAvyEyKtZG6DcMY9liAihJivR70xk+Rm2q1HsMwDAB7TDcMwzCikrAClIhcJSLrRcTvdowIJfd1EVnnyt4eTp6GYRjG8CDcGtQ64ApgWSgBETkGuAU4GZgNfEJEIjPImWEYhjFkCStAqepGVd3chdgM4F1VrVNVH/AGTlAzDMMwjJCENRZfmxKRIuAOt3t552MzcD7gPQ2oB17DGdbithC62sa7AqYBrQEwBzgQtrHDV9c4Vc3t7slH8UMw3eEwHHVFyheD5XqjVZeViejRFdwXwcY/Clw4ylh8ATJFwNyj6LgJ5wPdZcBvgQe7yjeIjl6Pm2W6IrtEq53DQddgsHE46Bosdg52XWGNxdddVPVx3LH3ROQnOOPzGYZhGEZI+uU7KBHJU9Uyd2bdK3AGizUMwzCMkITbzfxyEdmN837pRRFZ5KaPEpGXAkSfFpENwAvAV1S1shfZPRKOraYrokSrncNBV1/pNV3Ro9t0uUSkk4RhGIZhRBobScIwDMOISixAGYZhGFGJBSjDMAwjKrEAZRiGYUQlFqAMwzCMqMQClGEYhhGVWIAyDMMwohILUIZhGEZUYgHKMAzDiEosQBmGYRhRiQUowzAMIyqxAGUYhmFEJUMiQInIVBFpFpEzgxy7zJ3mo7u6LheRjSLScBSZcSLS65F9RSRVRFaISKWI3BCQfrqIfCgiDSIyJiD95yLyhoj8R0R+7qbFicjzIhJVPoywLx4UkXfc5a4QMqeKyH1h2DtPRN5y7+9SESlw07NEZKGILBeRh0REAs6JFZEtIvJ9d3+miDzYWxv6ggj74Zsissy9T38SkdggMuGWiXNEpMhd3heR1W76d0TkXTfvNj+IyCfd9OUicr2bliMif+utDZEkkvc/4Lw/isiSEMfCLQeTRGS1iNQE2hyqDIrIRBF5QUReF5E/BaRfICKvuWXpm27aF0XkU70yrK9mmOzPBfgzzsy/ZwY59kSw9KPoygYSgK1HkXkIODEMe2OAEcC9wA0B6elACs4MxWMC0uMCtt8AZrnb3yVgZuNoWCLsiynu2gOsACYFkXkWyAvD3sB7+3ng5+72/cCN7vbvgQsD5L6KM3XM9wPSXgByB/r+95EfAu/Rn4CLg8iEVSY66foOcFfgb8Dd/idwrvt72AykArHubyPdlXkYmD2U7r97znHAv4ElIY6HWw6SgKzOtoUqg8BLwMhOOnKAhYG/Fzc9ASjqjV1R9fTdG0TkFGA/QWbpFZGZwIXAQyLylJt2j1t7eVdELu58jqoeVNWQtSeXs1X1PVffvSLyV7c2s0ZEpndls6r6VHV/kPTDqloTJL3JzSsWqAH2uodeAq7qKr/+og98scVd+wEf0NJJZypOUChz958QkUdF5EX3iS+vK5tb761LGvCBuz0Pp7CBE3zmuXmkAB8Hnu6kaglwaVf59Qd94IfW35/g/FFtDZJtWGWiE9cBf3Pz3hKQ3ojzO8gBylW1WlWbgWLgFFdmwMtEpO+/y38BPwmRXyTKQZ2qHgqSfkQZFJFxOAHtl27LQ2vt6GLgEPC8iLwiIrPccxuAChGZ3JUdnRn0AQq4G+dp9whUdQPwCnCbql4lIscDZwGnAxcAD0gPm8hEJBeo7pRcrqqfBH4G3OzK3RfQZNG6vNqTvDrl+xCwHeeHf9hNLgaO7a3OPqBPfOE24WxX1ZJOh6YDndPWq+rFwPPA1e75DwfxRWCzxMUisgr4Ms5TIjhPk5XudqW7D/Bt4MEgZkaTLyLuBxG5G+cas4BdnY5FrEyIyHHAYVX9qFP6PGAksAwoB3JEZLSIpAFn0u6faPBDRO+/iBTiXFdpiPwiUg6ORqcyOAo4AfgGcAlwn4hkuumTgU8Cd9JxUsJe+aVfpnzvK9ynjVWqelDaXxEcjWnAO+rUOytFpAznaawsTFNWu+uPgI8BqOo9YersgKreJiLfwHlyvxDnSTFq6CtfiMh5wOdwCkJ3CPTFJABV/cLRTlDVF3FmhL4a5yn1aqACp8m10l0fEpF84ARVvUdEFnTTnn6lr/ygqj8WkZ8AvwIWAL/pQm9vy8QNwF8CE9ygdT9wiWsnInIrTjNaLfAh7a0KA0of3f+7gGuAjB6Y0uNyEIogZfAQ8KGq7nGPrwGmuOmvuzXutd2puXXFoA5QwPFAoYicjhOdp4vIp1V1Z4BME+3XWQzc4jZVpAN5wIGeZKiq5W6VukNywHbrS9z7cJuFAm1R1fN7kp+rK0FVG1TVJyK1QJ17aCqwrqf6+ojjibAv3KaSHwEfV9X6IHluAsZ3Sgvmi4dx/ggC+UhVP9N6b920Strv7RvARThNTRcBz7jXlSsirwCjgXgRWauqLxA9vjieyPuh9fenInKY9nsERK5MuDWHK4C5AXlPxnkH+ClVbbNLVZcB57j5Pgu84x4aaD8cTwTvv3t9I4AngURglojcrao/DtAXdjkIdTEhyuBWIMm1rR6YCezEadl5wD2vgPaWHnD88liofELS25dq0bYQ4sUjTnVzGfCwu/9DnB/zf3CeyDrLn4XzPqHOXV8RROZXwBx3+17cjg44TQ1PdNPeF3Ca7NYBv3PTprp5VgDLgS+56U/jdJx4C7g/QMd3gcsG+t73oS/WuUuRu8wJIvMckN85X5wn8Xu7YevNOMFoKU7Tyzg3PRt40fXDrwFPp/MWcGQniV6/pI5yP/zavf/LcIJFbBCZSJSJc4CnO6UtxPlDbP0NXOym/8z12WLgpAD5h4HjB/reR/L+B5w3ntCdJMItB2k4/z17gZXAfW560DIInAe8CbyL+z/lpn/HvbYVwBluWgLwRm/uobgKjB4gIuOBu1X1lgG0IQ74F06A8g+UHQONiJyG83T3gwG0YSbwBVX9+kDZMNBESZnIAR5S1WsHyoaBIhrKQShE5Is47yQ7dyzq+lwLUIZhGEY0EpFefCLyexEpE5Ggbb8ikikiz4rIB+J8bHpMJPI1DMMwhi6R6mb+BE7PslB8D1ijqscBnwF+GaF8DcMwjCFKRAKUOj1qjvjIK4CZwOuu7CZgvNtl1zAMwzCC0l/dzNfidB9dLiInA+OAMQT58Mz9vuFWgMTExDkFBQUA+P1+PJ7IVPgipsvfgF9j8Hgjcxv78hqLi4sPqGpud88P5Ye+trP3ioafL8wP4emyMtFDdQPgi4h1knB78SxU1SPeL7lfe/8S5+vjD3G+fL5FVdccTefcuXN11apVABQVFVFYWBgRWyOma0khRfXXUXjJreHrom+vUURWq+rc0GeEJtAPwXSHg/miZ1iZiJyuaCwTDc0tLFm6jJNPO40Wv+JrUXx+xdfip7lFafErzX6/s27x42tNa3HT/EqL35H1rfs5G5tOYcLM09vPa3FkfK3yLa68X91jR9FZU8kJ08YxJjOJgsxExmQmMSYzkYRYb4+vs7u+6JcalKpW4XyJ3Dqe1w6cb4AMwzCGLL4WPxV1zRysbeRgTRMHapx1+76b5u7XNbnDTb7+WgRyv8BZbd9wxJEYjxDjFWI8HnftbHs9QqxXiPF62mS8Hmd7d5WfNW+W0NTS8auWnJR4xmQmMiYzkYKsJHfbWY/O6F0Aa7Oz12f2ABHJAOrUGQLjZmCZG7QMwzAGLarK3/+zi1XbmiiqWt8WgJyg00RFXRPBGqm8HiErOY6clHhyUuIYn51Edko8WclxNDz1JHmfudENGkKstz1weN2AEusJTOso0xZo3r6adxsvY97HP9MecDyOfDeHYepAUVERZ589j7LqRnZX1LG7or5tvauijg/3HGbR+v00t3S84NxUJ4AVZLYHL19VS4hcOhKRACUifwcKcQZw3A3cgzMMPqr6O2AG8EcRUWA9cFMk8jUGJ7WNPt7ZfpC3tx3E+/IKJh13MgVZSQNt1rCgxa8crGmktKqR0qoGSqsbKK1qpHhbI7sTdjIpN4XJeSnkpMT16k9suCEi/PjFDdQ2tZC6azc5KfFkJ8cxKTeFkyfEke0GoOzkeLJT4tq20xNj8XiC398iOZPCU8eFb1xCDenqJyMpLnxdLh6PMCI9gRHpCcwdf+TxFr9SVt3QHrwOOcFrd0U9a3ZV8tKH+/D5lWtqNsEnz+0yv4gEqK6+3FbVFTjD+BjDEFVl0/5qlhWX80ZxOatKKmhq8RMX46Ep9Tge+dlSThibwSXHjeLi40aSn5Yw0CYPOvx+paKuyQk81Q2UVTW0B6GqRsqqGyitaqC8uhF/pyd6EYj1t/BKSftnjGkJMUzKS2FybgqT8lKYlJvCpNxkxmYlEeMdCpMgRI6ldxRS+bUvM/WJ3w+0KQOO1yOMTE9kZHoiJ43POuJ4i18prWpg3zdf6Ja+wT5YrBGlVNQ2sXzrAZYVl7OsuJyy6kYApo9IZcEZ4zl7Si5zx2fy/hdv4/0bv8kLa/fxw4Ub+NGLGzhlQhaXzB7Fx48ZSVZy5J7+hhJLN5fxq/cb+OWGtyhzA1DnphWArOQ48lLjyU9LYPqIVPLTEshLSyDfTctPSyAnJY6Pbvo8CQ/+lm3lNWwtq2FbeQ3bymp5o7icp1a3T2sU6xXGZyc7ASsvmclu8JqYm0JK/PD8O8lLS2DDgs/YE3g38HqEURmJFH/mmm7JD5tfVJPPz7q9h1lVcojip99nzMy5TM5LGWizhgy+Fj9rd1fyRvEB3igu54PdlahCemIsZ07JYd7UXM6eksuI9I61o8Ybr+bLhZP5cuFktpZV88LafbzwwV7ufnYdP/j3es6cnMMls0dx/qx80hKOmGl82FJZ18TB3YcYMyOTUyYmO8HGDTp5aQnkp8WTmxpPfEz3XlB/tGABhRmJjMpI5KwpHXv7Hq5vZnt5DdvKa9uCV3FZNYs3ltISUB0bkZbA5JiriI9P56MVJW3NhXmp8dZcaPSKIRugqhqaeW9nBatKKlhZcog1uypp9Dm9T2ITJvDMA29w5ZwxfP28qYzOSBxgawcneyvrnRrSlnLe3HKAqgYfHoHZBRl8/dwpnD01l9ljMvCGaGvvzOS8VL7xsVRuP28KG/ZVOcFq7V7ueGotcc94KJyWyyWzR3HujDyS4obsT7dbXH7CGE745auMv/mJPs8rPTGWE8ZmcsLYzA7pTT4/Hx2qZWtZrVvjqmHbtnhWHU7itX+vb5NLiY9hUm5yQFNhCpPzkhmXnUysNRcaRyFSnSQuxPnOyQs8pqr3dzr+ADDf3U3CmZYgIxJ5t7LvcD0rSypYVXKIlSUVbNpfhapTpTxmVBo3nDqOk8ZnMmdcFnu+dhsvfOp2/rxiJ8+t2cuNp47jK/MnW3NSF/j9yptbD/D3TY38+L032FLmzE4/Ii2BC48ZwbypeZwxOTvsl7IiwqxR6cwalc6dF07j/V2VLFy7j4Uf7OXVDaUkxno5b2Y+lxw3knnTcrtdSxhqlCxYcMQkQP1JXIyHyXmpTM4LmApqye0srbuOmYWfCWgqrGFreQ1vbz3IM+/taRON8Qhjs5MCglZKWyCz2rIBEQhQIuLFmS/mY8BuYKWIPK/O1MYAqOo3AuRvw/lgt9f4/cqWshpWlhxqC0h7Kp25tJLjvJw4LpPbz53KSeMzOX5sxhFP2+s/ex3/VTiTz585gQcXF/OHt3bwj5W7uOWsidx01oRh25beFSLwrafWUlnVyClTMrh6bgFnT81lan5KnzXhiAgnjs3kxLGZ3H3xDFaWHOKFtXt56UOndpWaEMMFs0ZwiX88PutbERWI0PZ+64zJOR2O1TT62N7pPde28hqKNpd1eIeWmxrP5NwUYhobWXToQ0TAIyCIu+383jwi7cdEENx1gLzHSWRnSRMnn+4b9rXvwUQkPHUysFVVtwOIyJPApcCRX4c5XIvTDb3bNDS3UPXXF/gNY1jl1pKqGnwA5KXGc9L4LG4+awInjc9i+ojUbvcyGp2RyM+vms0X5k3kF4uKeWBJMX9aUcJX5k/m+lPHDtsn81CICH/6/MnId7/B9J8+3u/5ez3CqROzOXViNvd+chZvbzvI82v2smjdfv7VeBWpMS180v8hl8wexcnjs0J24zUGjpT4GI4bk8FxYzI6pDe3+Nl1qK7De65t5TXs+qiaTVWl7gR24FdFcR5SFSAwzZVROEK+9Vuke5r9RLDXtdHHhD3UkYhcCVyoqje7+zcCp6jqV4PIjsOZOXKMqgb9UitwvKv8/Pw5Tz75JK991MyfNzQBMCpFmJrhZUqmh6mZXnISe/7RWU1NDSkpR3aQ2F7Zwr+2NLHhoJ/sBOGyybGcMTqm7WntCKqLqfFnkZKeE/x4DwllVyR0zZ8/v0fDugTzQ3/Y2RuaWpQPd+3irfJ01lXG0NQCGfHCySO8nDIyhonpnoj9RnpDX/ki2vwARLxMjHrkEfbeGplhk0Y+8gh7b7ml7bcwlMsEMDT+nyIwrfGVOO+dWvdvBH4VQvZOnBkvu6V7zpw5qqq6t7JOF086WQ/VNGokWLp06VGPLy8u10seWq7j7lyo5/1Pkb784T71+/1HCi6ep0uffzgiNnXHrnB0Aau0lz5u9UN/2NlrXF/UNjbr82v26C1/XKlTvveSjrtzoZ5x/2t6/8sbdf2ew8H92Jd2BdEVKV9Esx8ixWDwQ1/b2WuGgC8i0cS3BygI2B/jpgXjGuArPc1gZHoijWfOILOfOjGcOSWHMyafwSvr9vPzVzfzxb+sZnZBBndeMI3TJ0fmacToG5LiYrhk9igumT2Kw/XNLN5Qygtr9/LIsu38tmgbk3KT245PyrXPDAwjmolEgFoJTBGRCTiB6Rrgus5CIjIdyARW9CaT/u6xJCJ8/NiRfGxmPs+8t4cHlhRz3WPvctaUHL59wbQj2tCN6CM9MZYr54zhyjljOFjTyCvr9/PC2r388rUtPLhkCzNHpnHJ7FF84riRNtSSYUQhYQcoVfWJyFeBRTjdzH+vqutF5Ic41bbnXdFrgCfd6tygIcbr4eqTCvjk8aP4yzs7+fXSrXzyV29x0bEj+FZKZuTmJDb6lOyUeK4/ZRzXnzKO0qoGXvzA+SD4p69s4qevbLKhlgwjConUWHwvAS91SvtBp/17I5HXQJEQ6+Xmsyby6ZMKeGz5Dh5bvp1FzZ/njLx6ph2uZ2S6few7WMhPS+DzZ07g82dOYNehOhZ+4HRZDxxqKaG5kZcPfNA2944vYA4dX9s8PYFz6gST8ePzK4X7P4jYnEaGMZywDwJ6SGpCLN/42FRuPG0cv/7j/fx59xzm/byIz542ji8V2se+g42CrCS+VDiJLxVOYmtZDQs/cL6xKt9bQ9zhspDz5Xjd/Vivh4RY97g7pUHg9AcxHuHkvy8EvjnQl2oYgw4LUL0kJyWee6YtZVbeGN5hLo+/uYO//2cXt549kZvOnECyfew76Jicl8Lt503l9vOmOu88f/5ERPSWvDA2InoMY7hhb1DCJCehhV9cNZtFt5/NGZOz+d/FxZz9s6X84a0dNPq6NymXEX2ULFgQlboMYzgRVoASkSwRWSwiW9x1Zgi5sSLyqohsFJENIjI+nHyjkSn5qTx841ye/fLpTM1P5b4XNnDOL97g6dW7O4z4bBiGYXSPcGtQdwGvqeoU4DV3Pxh/An6uqjNwhkYqCzPfqOWEsZn87ZZT+PNNJ5OVHMe3nlrLx3+5jEXr9zPIOjAahmEMKOEGqEuBP7rbfwQu6ywgIjOBGFVdDKCqNapaF2a+UY2IcNaUXJ7/6hn85voT8bUoX/jzai7/zdus2HZwoM0zDMMYFIQ1Fp+IVKo7bYY4A1xVaKdpNETkMuBmoAmYACwB7tIejMUHg3usqxa/8uZeH//e2syhBuWYbC9XTo1lfHrHwWjDtavep+yr9bOvxk/Nm+8x/8qTifPauGO9wcbi6yXD0A99bWevGQq+CDb+UeCCE1DWBVkuBSo7yVYEOf9K4DAwEafX4NPATV3lq0Nw3LH6Jp8+8sY2nX3fIh1350L98l9W69ay6h7Z5ff7dV9lvS4vLtcn3tqh//Xch3rdoyv0lB8v0XF3LmxbJn7737qltF03Nu5YjxgMY8CZH8LTZWWiZ0TlWHyqel6oYyJSKiIjVXWfiIwk+Lul3cAabZ+O4zngVKD/52sYYBJivdxy9kQ+fXIBjy3bzmNv7uCV9fu5eu4YvnbulA6yTT4/Ow86c+U40w+0z1pa29Re+UyNj2FiXgqnT87uMPHb2JGZxP2sub8v0TAMI2KE+7HO88Bngfvd9b+DyKwEMkQkV1XLgXOAVWHmO6hJS4jlm+dP48bTxvPrpVv567s7efq9Pcyv3s5fdqawrbyGjw7Vdej9Nyo9gUl5KVw1t8CZddQNRLmp8UGnkii58foBnW3VMAwjXMINUPcD/xSRm4CdwNUAIjIX+KKq3qyqLSJyB/Ca+55qNfBomPkOCXJT47n3k7O46cwJPLhkC4tW1DD6UB0zRqbyieNGttWIJuYm9/jD34GeDtwwDCNcwgpQqnoQODdI+iqcjhGt+4uB48LJayhTkJXE/1w9m9teeoDxDzwx0OYYhmFEBTaSRBRhIw4YhmG0YwHKMAzDiEosQBmGYRhRiQUowzAMIyoJd7DYq0RkvYj43Z57oeS+LiLrXNnbw8nTMAzDGB6EW4NaB1wBLAslICLHALfgDBI7G/iEiEwOM1/DMAxjiBNWgFLVjaq6uQuxGcC7qlqnqj7gDZygZhiGYRghCWuw2DYlIkXAHe73T52PzcAZYeI0oB5nWo5VqnpbCF1tAzIC04DWAJgDHAjb2OGra5yq5nb35KP4IZjucBiOuiLli8FyvdGqy8pE9OgK7otgA/QFLhxlsNgAmSJg7lF03IQzgsQy4LfAg13lG0RHrwd2NF2RXaLVzuGgazDYOBx0DRY7B7uusAaL7S6q+jju4LAi8hOcAWQNwzAMIyThjsXXLUQkT1XLRGQszvunU/sjX8MwDGPwEm4388tFZDfO+6UXRWSRmz5KRF4KEH1aRDYALwBfUdXKXmT3SDi2mq6IEq12DgddfaXXdEWPbtPlEpFOEoZhGIYRaWwkCcMwDCMqsQBlGIZhRCWDKkCJyM9FZJOIfCAiz4pIRhi6ujVMUxc6LhSRzSKyVUTuCsOW34tImYis662OAF0FIrJURDa41/f1cHWGyCdqfBEpP7i6IuKL4egHV4eViSjwxZApE331XUBfLMD5QIy7/VPgp2HomoHzsV0RR/mG6yjne4FtwEQgDlgLzOylLWcDJwLrInCPRgInutupQHFv7RoMvoikHyLpi+Hmh0j7wsqElQlVHVw1KFV9VZ3hkgDeAcaEoas7wzQdjZOBraq6XVWbgCeBS3tpyzLgUBi2BOrap6rvudvVwEZgdCR0d8onWnwRMT+4tkTEF8PQD2BlIlp8MWTKxKAKUJ34PPDyAOY/GtgVsL+bPvjRh4OIjAdOAN7t46wG0hfmh3asTHTBMPHFkPFDv3yo2xNEZAkwIsihu1X1367M3YAP+Gu4uoYqIpICPA3crqpVvdRhvggT80P0YL6IDnrih6gLUNrF0EoisgD4BHCuuo2ZvdUVJnuAgoD9MW7agCMisTg/gL+q6jO91TNIfGF+iA4/gPkiWnwxZPwwqJr4RORC4DvAJ1W1boDNWQlMEZEJIhIHXAM8P8A2ISKCM+7hRlX93z7MJ1p8YX6IDj+A+SJafDF0/BBOr4z+XoCtOG2ra9zld2HouhynbbYRKAUW9ULHRTg9UbbhVMt7a8vfgX1As2vTTWHoOhNQ4IOA+3TRUPZFpPwQSV8MRz9E0hdWJqxMqKoNdWQYhmFEJ4Oqic8wDMMYPliAMgzDMKISC1CGYRhGVGIByjAMw4hKLEAZhmEYUYkFKMMwDCMqsQBlGIZhRCX/HyvqHuXvUtXxAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "tw.diagnostics.plot_extendedeventstudy(bdf, periods_pre=2, periods_post=2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
\n", "\n", "Warning\n", "\n", "Be careful not to include too many clusters!\n", "\n", "
" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "checking required columns and datatypes\n", "sorting rows\n", "dropping NaN observations\n", "generating 'm' column\n", "keeping highest paying job for i-t (worker-year) duplicates (how='max')\n", "dropping workers who leave a firm then return to it (how=False)\n", "making 'i' ids contiguous\n", "making 'j' ids contiguous\n", "computing largest connected set (how=None)\n", "sorting columns\n", "resetting index\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbUAAAEYCAYAAADbKGjtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAACSVUlEQVR4nOydd5gcxbX2f6dndjbNrlY5BwTI5LQLGAxoBSYaXWycyJYDXGxznS8GA0bAB8Y2OAHXZGMwNjYYY8AyQaBFIoMwUYAQCATKabV5J/T5/qia3dZoZnZmZ1aIpd7n6aerq09X1dun4unqKlFVHBwcHBwcBgO8DzsBDg4ODg4OpYJr1BwcHBwcBg1co+bg4ODgMGjgGjUHBwcHh0ED16g5ODg4OAwauEbNwcHBwWHQYKs1aiLyoIisFZHzs9xvFJE9CgjvQBF5RUS6RGRCFpkaEbmzv2m2YXSKSJM9vm79horIQyLymIg8ISJ7WH4dIvK+iDwtIudY2YiI3CsiMwrkd7aIPGPDv0pEJIPMZBG5vghu00TkScvtCRHZM+3+RSKyJHD9gojEReRdETnZ+o0QkT9bd6E6PEFEHheR+SJyv4jUZpApSocisrflNl9EHhWRqdk4Wh3GRGSpfSedIrJ7imOh/DLFkeFesTqsEZGnRKRZRE4J+M8K8GgSkfHW/1Z7vUFEfBG5LZBHvcDzheoyY3wZ0lqMLncJhP+UiKy3/icF/F8Xkb9b/3dEZGOp6hz7zI9FZK6N69AM94vlWBsok8+KyGGBtK4M8Ky3+bVdRP4vS1iF6vAHgfCXisiVpeZnw9iiLcihw8tEpKGgCFR1qxzABGAWcH6W+7OBUwoIbwgQBZqACVlkfggcX2S6l2TwOwu40Lobgb9afj8Gzsd0Fp4Ctrcy5wJ/KZDfjgH334DDMshcBexTBLcwINZ9KHBn4N5om+Yl9toD3gbOBH5q+Q2x964D9uyHDiMB98XAt0utQ2AMUGPdxwC35eDYk0eBscCrAdnrgN8Xwi9THAOkwzHp7z5XWbP3pwJrA/n4XOC4wP1CdZkzvlLoMi2sLwHXZvD/P+AE624Ans+Wrn5wPBq4bCA52nIWDujoOetuBG5Mk50AfB94rRT80p6dA3xyIHRI321BUIcTgL8XEv5WG6mp6gfZ7onIMAzJ82xLHRKR08WMVJ4Rka9lCG+Tqrb1Ee3ngQdtHLNE5B4RuVtEXhWRg/NM+hgxI7K7RWSK9XsdSI0qhgJrLL/VNm0+kACSVuZxYGaB/N4KXHbb8NJxiKq+YPnNFpHbbY/7RRHZqS9iqppQm3Msn5cDty8Afha4HgGsBLoAH1gM7G/vzQFOpXAdxgKX1cBrGZJZlA5VdZWqttrL9Pe4Gce0PHoicEfger71y5tfpjgyoBQ6XJXl9mliRsKXBEdhFnsBi+jNo3OAL9p0FFwe84gPSlMeUzgF+FPQQ0TKMA3PPwFU9XlM+SxLf7ifHL8EVIjII2JGuENKzVFVfVVN5dH0MnmkiCwQY7mptPl1I1ArIlUl4Jd6dhSwnao+XWp+lmOutiBdhx8AE9L55UI4X8GBhKpuEJFbML3ZP4nISMxoaF8r8pyI3KeqawsMeriqtqfFdbyIHAj8AFggIt8Evpzh2eNVdQMwRVXXiciRwE3AYcBC4GIReRWoAw4KPijGNPeOqr5rvZ4H2oBLC+UnItMxo4b5af4jgdY08bWqerKInAR8A/iRiFwETE+Ti6nqETaceuBqYBJwvPXbEYiq6svSa/Vci2nY6oByy/nf9t5i4GvALRSoQzEm3e8BncDl6fwpjQ4RkWrg/wEpE3ImjkGcDHwhcP0CsIoCdNhXHKXSYRb8E7jNuv9g+dwWuH8KZrSdwmJgd+h3eewrPiidLocDOwFPpMkcDcxX1c6A3zpMvt0M/eQ4DlivqoeJyFmY0e05peYoxnT7V2AaplyBqXN2VNUuEbkU+BFwib3XDHwC+E+R/FI4AWMdyoSS6DAHMunw3XR+ubBNNGoZMBV4JdWTF5FXgO0wFWsxWGjPy4DhAKr6e4xZKSNUdZ09Pygi11jvszFD4l+JyAHANcBn7L3tMaaCmTnSkRc/Mfbwy4GZgRFVvvwOt+m+MNcDqroQOEBE9sM0bvthzBY/TZNTETkDuBZjInkFWJEj6Lw4qupNwE0icjbwv5h3my/HvHRoe39/BX6uqous9xYcAxgJdKrq0hzpyIdfrjiyoWAdZoKqbky5ReQO4EhsIyMidZgG7F6Meacv9Mk1V3w5ULAuLb6MMZWnl4lTgBv6eDYb8tHnBuAB634A+F0e4fanzlkOHCTGMtQE3B+wNgDcTu7RfyYUUqeejHmX+aC/OsyGYnQIbFuNWoze9CwF9hCRiL3e3foVivUiUh3oWQQLgQDk6lXYNHWqatI2MOsCz6bca4Bh1j0VM5LbI62nMQ3Ty8+bn4jsANwMfD7VsAahqmtFpCbdOwO/rL18EalQ1S7r1wx0BHhcY0cXY0Xkd6r6HVWdLyK/sPcPBFLmiWnAqxSowwzxZzIxFKvDZoyZ6h5VvSdwLyNHe28vTMURxDSM+bWQPJorjpLoMANn7DN1qtpsLw8F3gzc/hJwdwZ+rwbDpzBd5oovhaJ0Gejln4wZxQbjrwXqMaOMIEbQW1bTUWid04T5TjfXnjNN/ik2v7arare9bsGO5EVkiKpusv7p77eOzO+74DpVRKZh+rBvpd8rBb9cI7UcOpxCZn6ZoSX4aJvPgWl9X8NkhHsy3N8XY064CzMSOBN41h6nZ5CfhslcG4EFwDczyPwI0yhA4MMkpnfalEea98MMeecDjwF7Wv9xwCOYTP4MZmR2A+Z7UzemEDUB9Vb+XJuWQvjdb99Vkz0+k0Hm6kAcs7EfhTGmwVvy4PcZy22ePfbKILMk4H4JaMc0fk0B/+swDUGhOrwgwO9uoG4AdPgFjOk3Fc9V2TjSm0djwL/SZK7DVKZ588v2HkupQyt7H/AOplG61vpdavPmE5jOUVlAfr7l0FMebR79bBHlMWt8pdKllZ0KPJ/B/2vAlWl+f8Y0DKWqc8qBWzFl5SFgzADk13p6y+QC7AQx4FuYzxjzsWXF5tc3Mfm7aH72mYuB/8mRvlLoMGNbkEWHE4C78wm355lChD9qB+ZDa0EzZwYgDRGMmccbgLCnADd8yPxGAH8Z5DocMI7biA4HLI9+mLoELgP23crvclBz/CjwS03ndnBwcHBw+Mij6Cn9IjJRROaJyCIReU1EvptB5mQReVnMz9JPStpPvg4ODg4ODqVA0SM1ERkLjFXVF+xH74UY2/yigMyBwOuqulFEjgZmq+r+WYJ0cHBwcHDoF4qe/aiqKzGzwlDVVhF5HRiP+bEzJfNk4JGnyW8asYODg4ODQ0EoyTc1ETkK+C1mdlAdMElVWwL3fw3MsJfjgVpVLc8S1hnAGQCVlZX148ePx/NyW0l9388t43fhaxgvlLsN7zMcK7NkyZJ1qjoyp2AWpPObOHFi3vFmlSkxP8/zWLx4cck4Oh3mkbZtSIcDwg+cDvMJI03G6bAfKMHslBBmTcDdMKsuvAvskkV2BuZ/rdvzCbu+vl7nzZunfaFPmYen67x7rys+HCtDhinF/Tnq6+sLijcrSsxPVUvK0ekwj7RtozosGT9Vp8N8wkiT+SjpsCue0PmL1+gvHnhdz775IZ33xmr9YGOH+r6fNa5S6TB4lOLn6/1so3Yl5ifXcuA4AuZH6Fkd40bMKhS3lCBeBwcHB4cPEcubO2l6cw3z3ljLk2+voyOWJOQJSV/565vPARAtD7PDqCjTRkeZNrrGumtSA52Soyjzo5hFM+cBO2O+q+0FHAvsr6pnWZkZmB9MdwCWY5Zm+byqpq9okAqzZ9g8evTo+htvvJFoNJozHW1tbbllWhfT5g8jOmSLJeAKC8fKzJw5c6GqFrYdgkU6vzvuuCPveLPKlJhfNBplxowZJePodJhH2rYhHQ4IP8jKMeErazuUFe0+q9t92rtiVJRH8AQ8ETwxS1WYa3NMrezmK1/ctnWY8JXVHcryNp/lrT7L23w+aPPZ1BojWWbXWNbeJTl6amKFnTtX88PPTt3mdJjYtJhXW0awuCPKy2sTfNBmUj28QthzVIg9RoTYeXiI4X/4IwuPP9Vwb/NZ0eazvE1pifW2N/9vzTxO/cOV/dZhNhQ7UjsHY07cA/MD5xuYVTVeFpEzAVT1WhF5CrMKRxVGd+ez5TI9WPnrgesBGhoaNBqN0tjYmDMRTU1NuWXmzqap8yQaG7+QXSafcKxMMUjn19jYmHe8WWVKzK8vmb7gdNiPtGXh5/vK6tYulq3v4L0NHTy/9HV2Hz6FqkiY6vIQ1eXhXnckTFUkxMJnnixKhwPCD2h/4FLu3HgS8UmfZMmaNt5e086StW28t76DeDJ9taV4zrh+tvHx/MhkQSk5Jn3l/Q0dvPnQn3i0+VDiHUN4a3Urb69t6+HlCUweXs2eU6LUPLuAodOPJLXIdc9S1wJiryrvXUhjY9aF9LcqvzWtXTz25lqa3lzL/NeH05ooJ+wl2G+7YXzlkFHM2Gkk24+M9vABaAoJ/50hro3tMRavbuWtNW3sdt0/+80vF4pt1L6BWTapDbOc1I6YJXJQ1WvB/MeGWeR3pb3/pqruU2S8Dh9DqJpCsaqli9UtXaxp6e5xr27pZk1rF82bOql79XFCnhASIeQJ4ZAQ8jxCAiHPY9+aZN+RbUV0xZMsaxvOfzZV8M7jS1m2oYNlGzp4b30772/sJJbwN5P/25uZdujpxZT4Jo44bCBTnB2qyvr2GEvWtJmGa609r2ljxabvGaGXXiLkCZOHV7H9yCiH7zKaHUZG2X5UlKkjq1lz5hlMvPFGfB98VZKq+L6S9FNu2PTdv5c03bc9/R7Llsfpfm0VNeVhqu0RLe/tMIjAuk6fR99YzeLVbSxe1criNa28tbqN7oQPfA6ACa0b+cToGho/MYpPjImy4yhjcqsoCwHQNLGNxsZdcqanqfKQkvIrBElfeemDZpreWMO8N9fyynKz5OTo2nI+M+pNRtZO4IwTTqSmYosdffrE0OoI+08dzv5Th9PUdRLcUtTaxRlRrPmxDbMy8+uYtdBmYdYiO0lVX7MyYzFbp3QALwLrgcM18B9bWpg9w+ZRo8fUX37VDQyrraYyDGUem/UGUhho05Wq0p2EzoSS7Orgi58rndlj1s9uIx7rZli0nKoyobpMqA4LVWVQXSZEQtKTturqahIK3QmTlu4kdCWUrrYP2BSvgfIaupLQnVA67bnLynQnoTOWAC+Er5BUzNnXza5ry3wuOLB05sfho8fXH33OtQypilBdZnhFy4SqsFBtOVaXCe3t7cTDVWzsUjZ2K83dPhu7lOZupblL2dgZozkWIqFb6j9aBnXlQl25B5rA88Km8ktxtOeUe+akJD8+rXQ6nPLtmynDp64yTG1EqIkIQ8qF2gjURITaiOG4srmDNipY2+GzpkNZ22nOzd2bl8GKEIyq8hhZJYys9BhVJfbwkFgHkcpqupPao9cuq2fjB2Mfb2LXkw8vmenqtMtu5d9LY5SFPMIehAXCnhh3zyG0dSVY3SW0BwZZ5SEYW+0xtloYW7aBYRUVbDeqjtFVQtjLuOXPh2BCHlNfMevGPp8Jick/KQwtF8ZHPcbXmPME7wOGVNQwfNi2Z0I+7bJbefq1dYSGDyUSgvKQ9JzLw/bagxUtMV5vFtriZgS5Q53HHiND7DEyxKQaD2l7a5v5DJANfY7URGQuZmfddJwHJFS12+4tdCNmksjfVPU1EbkYM7PlXmCliFyB2VxyIWn/sQURHDZvt/MeetFCIbV4fMgTqiMhairKqC4P2V5UmLZNXYwcEUUCw3cReq9Xt7I2OZRR42tTnBArA1i38P6KLiprK2ntStDenaCtO0FbV4K2WIJU2//13TL+iZA30s0C97/r09YtmDV0t0Qk7FFbUUZnt9Cd7CThZ+qEpDJYbxgiUB00S5WHINHG8KFD7cjFVCrm7PVcdzavKan5cbud99AX1wsdaxIkM6a9J8WYLdV6URUJMaa2glHDytm143G6QxNp2Ht/xtRWMLq2nNG1FYysKe/pAcOHY348eKfxLF62EsqjvNsWY92abttzz8SxGxEYU1vBpGFV7Dm1iknDqpi08krWefvzueNOZmhVWcbOW978JpeV1Px41IF7sq7zRUaMHkMs4RNPKt0Jn3jSJ5bwiSWNuyzRynH7TGCHkVF2GGVGXmNrK/BSjdfcRmNiPfb43On/EHQ457zDWPy9H1F38aW0dfeW//buZI+7K5Gk+t47OfDHZ7HjqBqGVKWNVFL8toHPAOn8Dth7F1a/fD8ydgc6Ykk6Y0ma40k6O5J0xBJ0xJJ0J3yG+DGOqJ9K406jOGTHEdRVRTYPeO5F28xngGzos1FT1U9nuyciq0VkrKrOEbNhpqjqpfa59H2k3sesBn4QZjXvPhEtD3PC7hEmTt2xp4ExGSxJW3ec9u4krV0J1ncqnRs6ehoeRdHgB9i2kbRrGetWtvR8mE2NUI3bPJOM+YwuT1BTEWbskAqi5WGiFcYEkXJ7a9/OJ+l546ULj+Dfc5vYvX4/NnXGMx4tnQnWrFrBJ6ZONmaRSIiq8nBvo/XK91kU/zSNR3yp5/tKZVloi4rRZLRP5kxPqTPa8OoIVxxWzfTp02nrTmzOraPXveTttzlw710YXVPBqNoKxtj334O555rCdNB2JU1fKfDLL+5JU9NGGhvNXrGqSnssyfq2bta1xVjX1s2mjjj+DdfQ8P/OZ8LQys0aYgDmvkFT5z4Mq45kiOHDxaE7jcZbVU5jY+7V7Uz+2n0rpaq0GFVTwaJTv8Bu44fklGuqPISGKcNyymyL+NzeExi6aScaG/fLKpP0lWVf+yrbffmWrZewAUCf5sc+RmoHYcyJIcxW93NUNeMmjyLyLDAEODfbzEcr1zNsxux2up7s+yGlkGvPpIGQqdYS/PSJ4fdmidJW6ncwuYQcnQ7zS9s2ocMB5FdqGafDLPg46DArivnJDfM9bRHGbvQoMMz6NwA3BuR2wMyK/GE/4ujz57xtTebD4DjY+X0cOA52fh8HjoOd37bOUbX4n6/3tedJqhrc2v157M60YmxgP8VsXnhlkfE5ODg4ODhkRbFbz1wN1AAPi8iLIpKaxj9OROZYmU8BpwKHWpkXReSYIuN1cHBwcHDYAkWN1FR1hyz+K4BjrPtxAv8X9gPXfwRlCkUp4h3s/La2TKFwOtz2ZAqF0+G2J1Mw3M7XDg4ODg6DBkXvfO3g4ODg4LCtwDVqDg4ODg6DBq5Rc3BwcHAYNHCNmoODg4PDoIFr1BwcHBwcBg1co+bg4ODgMGjgGjUHBwcHh0ED16g5ODg4OAwauEbNwcHBwWHQwDVqDg4ODg6DBq5Rc3BwcHAYNNgqjZqI7C0iT4jIfBF5VESmZpD5rIhMKiDMz4nI6yLSlUNmsoj0e9FMEZkmIk+KSJNN/55p9y8SkSUBfu+JSJeIPBGQGSEif+4Hv9+IyNP2OCeLzCdF5KIi+E236X5MROaJyETrv8W7tRwXiUiLiDSLyPeC/Ky7UI4/sHniCRG5VUTKMsgMiA5F5EARecXqa0JAh2+JSIeIvBsIIyIi99r3kje/wPN/tJvtZrpXrA63F5GFItImIgcF/Bts3nlMROaISM0A6XC2zStN9ghlkClWhzUi8pRN8ykB/xNE5HGbh+4XkVrrP0VEnrU8l9o65xNWh17g+UK5hkTkChGZa7nuMgBcIyLyNxFZICLPiMihWbgeZPPrKnsuRZ36i4AeV4rI/wwAv2zlsUJEbre8bxeRCut/s4hMKCiSgdikLcNmcGOAGus+Brgtg8wtwEEFhDkcqACW5JC5CtiniHSH6V30+VDgzsC90cBfgCUpfsBY4GvAyrRwrgP+WSC/He3ZA54Cts8g8w9gVBH8IgH314BfZnu3lttiy3MmsBYYEuC3Zz90GIz/VuAzW0uHmF3Yo0ATMCGgw9GWX0taOOcCjxTCzz63h9X93Cz3i9VhFTAs/d0DdwHTrXs28M0B0uFs4JQ+ZEqhwzHpcaXln4uBb1v3HcBxludczA7Qt1kdHhd4plCu3wTOGGCuxwJ/sO4pwDNZuJ5j+R1g81DRdWras68A4wZIl5nK45nABdb9U+BM6z4A+HUhcWyVkZqqrlLVVnvZDSSC922P5yjgKhG50/pdaHtnz4jIZzKEuV5Vs47SLA5R1RdseLNtD+BeMXu67ZRHuhNq3yxQC7wcuH0B8LMgP1VdidnhO33rg1eAwwrk95Y9+5j3lQzeF5EaYKSqrrHXt4jIDSLyL9tDH5UHv1jgsodflnebBNZYPXYAm4D97b05mAJfqA5jVk4wjfeSDMkcEB2q6iZVbQvIpXS42vJLxxvAgYXws7gAuCzTjRLpsENVN2S49RpQZ91DgTUMgA4tzrajiO9kuV8KHa7K4B/Mv9WWM8BeqvpPy/NfwM6YMjQH+KJNR8F1jn12shirxtUiEik1V+BtoNyWiZTeMnF92vJ7GtMZKbpODTy7D7BazRZiJeWXo06dDtxv3ffZayy/w/oKNz2SrXZglPEUsEuuXgWwF6aHJZiCuRjwsoSZcaQGjAQeD1zPBn5j3ScBV1j3RZjeevB4KPBcvU3zcmB/67cjcEt6/JbfC8BTaWnZGVjWT34nA3/M4L8v8Ke09/c96/4JcJZ1X5eB362B5z4DPG/TsEO2d2vT+gawA/Cs5XNCgN8/+6ND4DzgLUyFU7W1dBi41wRMyJBH30uTqwQ2FMIPaAQuxfS4txiplUqH6eXHXu8BvAe8CjyJ7SGXWoeYUb3Y9zMXOHggdBh49pQ0v69jOo3PYjoIAIsD9//bvoddbBpf6m+dA7wZ0MkV2NFEKbkC5ZiR15vAamC/PrhWA63AJ0tYp/4K+OpWrlMfAqZY93bAgwH5Z4ChmdKa6Shqk9BCIOZ7yV+Bn6vqoj7EP4HpiSjQLCJrgBHYXksRWGjPy4DDAVT1wlwPqOpC4AAR2Q+z0/d+GGX+NCgX4Hct8KU+0pEXPxH5NPBVjKkoHwT5bW/T/9+5HlDVfwH/EpEvYUYUGdOuqioi38T0nJZiGu9MPbkU8uKoqpeKyGWYdzsL+L9c6aV0OtwCwTyKqbRyIR9+5wAn0DtiygcF6zALrgU+p6oviMi5wPdV9Zel1qGqrrfOThG5G2gAFvSRtoJ1mA2qehNwk4icDfwvcDbgQ48+vwP8S1UXiUhljqDy0ecG4AHrfgA4Po8kFsr1K8D7qvo5EZkC3A3sk4mriJyHya8rMY1gLuRb54QwptvZeXDrD79s5XEDveVkiL3uF7ZKoybm4+yfgHtU9Z4sYrFAehYDp9sh+BBgFLCukDhVda0172zmHUyWTdtF9A51e9KiqkeISIX2muGa6TVLTQWuMcljrIj8DvMt5h5Mbyi9YZhm0583PxHZH7gEOFpVOzNQfAMzAuiL33WYDB3EMlU9LQe/LWB1eCamov4rpjf5dIDfq5ieXCEcK1S1yzaYm9LjH2AdZuLXk0dFJL1Rmwasz5efTfcYzPedSmBXETlPVS8NhFm0DjNxCTybSs8aYIcB0mGdqjZbmUbM6KCXTAl0mJXglrqtsu6XRORTmAZNMd+Bghx7wi+EK2bE0YAxk6fOpeYa1NtGzHezTFyrsfkV+K6qbmRL9KdOPQx4XlVb0m8McHl8DDPf4kV7fsyGJRgLTiZ+GbG1RmrHY8xco8XMXnpFVdNn1twPXCwir6vqf4vIk5ghqgf8UM23pR6IyMHAhcA4MTPL/k9V704Lc4GI1NueQUb00as4TER+TO/3rO/ZZw4IpGMJMB9TmPcBrgQiNk3/rapvY5T0x0L4ATfZ8z228fxhkIeqtorIOhEZreY7UDZ+uXr5p4jIqZiebTfGVJPx3dp0fg44AvgF8LD22vmPAX4PTCqQ45Uisiu939My6WJAdCgi0yyvPTETfl7H5NHdROSXQNRy/6mqPmk5/i1ffmq+d+xl45oC3JjWoJVEh2Jm/N2NMa/tKiJz7Ps4B/ibmBmsPnAKphyWWoe/EZFPYCq0JlWdkyGZxeoQEbkP2BXoEJGDVPVMzGgl9b1lA2ayE5gJIffSa/b/vYi8ghmV3hkItqA6B/PO/iAiZ9r4Th0Arn8C/iIij2Ea6fOsfzrXe4AfYj6FDBWRq4qtUy1OsWnIhgEpj5j682YRWQB8gLFQAXwSM0Erf+Rrp/woHphe8A0fchpGAH8ZoLAPAC4erPy2IR1GMJVkxm8QTocfbx1+mFyBmwl8E94K8W3z/FJTKx0cHBwcHD7yKHpKv4hMFDPFdZGIvCYi380gc7KIvCzmZ9cnJe0nZgcHBwcHh1Kg6JGaiIwFxqqZZVWDmQ3zWQ3McBSRA4HXVXWjiBwNzFbV/bME6eDg4ODg0C8UPVFEzQ/HK627VUReB8YDiwIyTwYeeRqzgoODg4ODg0NJUZJvaiJyFPBbzI+DdcAkDUwJFZFfAzPs5XigVlXLs4R1BmZZGyorK+vHjx+P5+W2kvq+n1vG78LXMF4odxveZzhWZsmSJetUdWROwSxI5zdx4sS8480qU2J+nuexePHiknF0OswjbduQDgeEHzgd5hNGmozTYT9QgtkpIczSLrthfuZ8lwwrhljZGcAq4PZ8wq6vr9d58+ZpX+hT5uHpOu/e64oPx8pg/uMo+t3V19cXFG9WlJifqpaUo9NhHmnbRnVYMn6qTof5hJEm81HR4cb2bn38rbV6bdMS/f5f/6Mn/fYBvfCfr+oVD76h1zYt0duffk//+eJyffT11frc0vX6+spN+sDDj5ZMh8GjFP+p7WcbtSsx/zeUY/5I32zVEBHZA7gR86/ILSWI18HBYRDB95UNHTG6k4NvRnY86fPBxk7eW9/OkysSdL26kqpImKpIiMpIiOqAuyqy1RZ6KhiqsLHbY+6i1by2ooXXVmzitRUtLG/uXR9idG058ViSlzd8QFt3As2izl9taBqQNBZlfhSRYcA8zE+OKzE/mx6LWc/rLCszA7MUyg6Ytb62Az6vW/4onQqzZ9g8evTo+htvvJFoNJozHW1tbbllWhfT5g8jOmREceFYmZkzZy5U1YacglmQzu+OO+7IO96sMiXmF41GmTFjRsk4Oh3mkbZtSIcDwg/wWxazvGs4HeE61ncp6zp91nUq61PnLiVhfweuCsPQCqGuXBha4TG0XBhaYY9yoYZOvvDZbUuHsebFvNc5nLbQEFa3K2s6fNZ0KKs7fNZ3KX4BVe2uXWv4389u96Hq0FdlTYeyrMXnvZ4jRmvc7C4kwOhqYXKNx+Raj0m1ISbXetREpCccX5XuJHQllM4EdCS0xz39/tvZ//Zb+q3DbCi2S3AOxpy4B+bnxjcwq9S/bP+6R1WvFZGngHGYP+QVOB+zCsIWUNXrgesBGhoaNBqN0tjYmDMRTU1NuWXmzqap8yQaG79QXDhWphik82tsbMw73qwyJebXl0xfcDosLG1JX9nw71/wQPsXGLprI+vaulnXFmN9Wzfr22KsbzfX69q62dQuDKtRqiNhohVhasrNOVree71p07ucd2zutBTCb+j2e9HU9Bx7T9iRmoow1ZEw1eVh4y43cZeHPR6Z18TU3ffjg+YOlm/s5IONnSxv7uSDjR0sb+5k5cYxJDSEqSIMRkQjjB8apWFcJROGVjJ2SAXL//xXuo84ltUtXaxq6WbJpi7WrOjarFH4fxufpBikc5y8274sfO4Z9qv/JJGQRyTsEQl5hDzBrugDwL/nzmPEjnuzbEMH765v57115rxsQwcrN423Ut0A1FaE2W5ElAMmVDN5eBWTh1czZXgVXZf8lGGXXEZHLEFHLGkP4+6MJWmPJfD+vpDGxq+mJ7vf/N6Uidwx7xlG7DiMIZVl1FaWUVtRxpDKsp7rIZVlvLj8ZRL+SBataGHRyhbaus3i/2UhYcdRNRw14kXKKyfy2SOOZuextVSXZ25C8ioPY8Jw+y395pgNxTZq3wDagTbgP5glW24G05iB+Y8NszDrSnv/TVXdp8h4Bw3eWNXC+60+b61uJeQJYc8jFBLCnuCJOYdCQmdCaetOIECqjAli3MkQcR9iCR8RrIxsLiuCqhJL+MSSPvGETzxp3LGETzypxJM+77UkMye0n2iPJfj7WzH+E1+MJ0LIA88TQiKELEdP4K334rw1/x2Ttp5Dbdp84u8fyfvxodz55xd60h5P6ubyCaWlrYOK55tI+krSV3zd/Jz0ldN22mIfy6Jw61Pv8uZ7cd578t2AbgCrAzA9uf+8HaOp5TXWBRqr9W0xNnTEUP22EXzxmZ5wQ54wvDrC8Gg5I6IRthtRzcZ1q6kbXkdbd4LWrgSrWrpoW5ugrStBa3eCWMJnUo3Xs7ZSKTD39dXc8lqMW157MatM2BOSvqIPzevxE4HRNRWMH1rJPpOGMr52Hh2hTzBj+pGMr6tkfF0llZEtddEU/ySNjbtt5pf0lXVt3aza1MWqli5GXnlvyfgpMOOKJnPx2OYrMolgGrmQhwi0dCVg7uM990dEy5k8vIoDth/OlOa/0hHehaMPn8nk4VXUVWXamQaaTv0Cu4yrzZmmpvAhxVDaAqNqy5mQbCXuCSuau3h9ZSstXXFauxJbyFZF3mfnsbUcv894dh1Xy67jhrDj6Cjl4RDMvYCmzpNomDKspOkrJYo1P7Zhtp54HbN1wCzMatEnqeprVmYsZnPCDsxileuBwzXLSv3BYXPd6An1I772e8KeUOYJZR6EPexZetyiSSrKwvbayJWFsM9BOL4OX6qIVlX3PFPmiZGRXrn2jk7C5ZXEfCXuYyp6H+JJevw+EY3xrRNLZ/YY9o2b6NgyX31oGBtv5Wczx5TM/Dhk9MT6obP+D0X6eGpzhARCXko/ECaBJx5l4ZDRfeq+B2HbWIY98JMJIuEwntDTYIpgr81qA/XD4iXVYfVXbyLfz0CVYaiNCEPKhZqIUGuPGl1LJFzJ6LohPf5VZYZDEH2ZkOK+0tzSzsi60pkfb739L0Rv+RPvfvlkupJKVwI6A2aklF8yHmPskHJGVHqMqBSGVZhOWQ+2YRPyd371J6Lzmlh38HQSPiRseU/4kFBz7StEvTgT6ioYVSWMqvKoDA8cv61hQvZV6Ygbs2BHXJnw9zvxvvrlLfJdD7YhHWZDnyM1u6jrmAy3zgMSqtotImdhJoGUA39T1ddE5GLMzJZ7gZViVj33MT9nb/YfWxDBYfMue+yth0yKMGbceLoTfs8oI5ZIBtw+6zZswi+vIpbw6e45kkYu6aM6xIYeyxRlkC1B00gmfHXXjH8i5I10s8AvTmng/V//jmHf+nbPSCLpKwlfSfp+j/utJW+z/fZTez66KuajraLoWzfwTnxPtttpX+vXe48eN7z37rtM22EqZSEhEvIoC3uUhTzK7TkS8mj79S9pbOxr15XCOF5xVJTp06fjK5uPmlRRH5KqPPHEE8yYfjBlIaHM8/C8tEI1t9GYH2eekTPuD8P8+PD5h/PeWWcx8SqzGHz6+0cBgeb/+RbTbrk5c6A9/HLvZrI1TMjp/I44bAZNIeGkQWxC/slJn6ZpXHjQfAbot4l8SOgjo8Ns6LNRU9VPZ7snIqtFZKyqzhGRVzEjv0vtcz9NE38fs7L0QZhN3/pEVSTECTtFaGzcNaeceYGZh+uqSvzhw3ik/UT2O/y0noYw1Uh2J5I97nVXXsHkn5xDedijPByiomzzc3nYY/78x/JJet44dKfRNJ14OI17jsvNMbmMxkO2z3wz8QxNndvTeOiOucNoWkFj4w65ZU4+Nuf9/kJE7Ogrcw+wJiJEs9jnt3UMrY7w0ldOYq9o7g7PolmnMW0rpcnB4eOKPs2PfYzUDsKYE0PAicAcVT07SzjPYvbxOTfbzEcr1zNsxuwhtZ6+91IbsZVlqrUEP31i+L1ZorSV+h1MLiFHp8P80rZN6HAA+ZVaxukwCz4OOsyKYn5yw3xPWwR0Ao8Cw6x/A2b/qJTcDhi73g/7EUefP+dtazIfBsfBzu/jwHGw8/s4cBzs/LZ1jqrF/3y9rz1PUtW1KU9VfR4zMzK1c+lPgWtV9coi43NwcHBwcMiKYreeuRqz3fjDIvKiiKSm8Y8TkdQOuJ/C7BB7qJV5UUSOKTJeBwcHBweHLVDUSE1VM846UNUVmO3hUdXHocD53Jvj+o+gTKEoRbyDnd/WlikUTofbnkyhcDrc9mQKhtv52sHBwcFh0KDona8dHBwcHBy2FbhGzcHBwcFh0MA1ag4ODg4OgwauUXNwcHBwGDRwjZqDg4ODw6CBa9QcHBwcHAYNXKPm4ODg4DBo4Bo1BwcHB4dBA9eoOTg4ODgMGrhGzcHBwcFh0MA1ag4ODg4OgwZbpVETkVoReVJEmkTkWRE5LINMo4jsUUCYB4rIKyLSJSITssjUiMidRaR7bxF5QkTmi8ijIjLV+p8tIs/Ye1cF+L0iIm0i8rKInGxlR4jIn/vBLz2OLRaFFpHJItLvRUFF5EgReVpEHhOROSIy3Po3pPnXWI7v2/e9KaXDFD/rLpTjCSLyuH2/94tIbQaZAdFh4P5FIrIkoMNVItIiIm+KyPeDHAvllx5HlnsDpcMfWM5PiMitIlJmOb5u+bWKyGVBftZdqA5nichSW7abRGR8BpmBKocVInK7iCyw55FWh8ssxzdE5JdWNiIi94rIjEJ1KCI/FpG5lt+hpeZnw3hQRNaKyPkBv0YRWRl4t/XWf5iIPGA5fiC2ThWRu0WkOu35QnT5g0BcS0Vki63CSqBLEZGrReQpEXlORE60/kNF5CGbj59IpVtELhORhoIiGYhN2jJsBucBYeueCjyXQWY2cEoBYQ4BokATMCGLzA+B44tI9xigxrqPAW6z7h0DMn8DPg2UYXaX3R14HngKGGJlrgN+XyC/9DgOyyBzFbBPEfwmAeXW/S3gEuu+C5ge0Ms3rQ4nWP09EdSh5bdnP3QYCbgvBr69tXRor0cDfwGWpPIoEEnlUeCtwLMF6zA9jiz3B0qHwXd7K/AZzA71b2K2i5oGtKXl0f7ocBZwfh8yA1UOzwQusO6f2nyarsPHgF2tzLlWF4XwOxq4bCD52TAmpL9LoJHAZssB/8uB0yzXm+1zzwEnAt8NyBWky7Q45gCfHABd7gbMs+4a4G3rPgu4MMD7r4H38vdC4tgqIzVV9VU1YS9rgZeD90VkGEYx59leQkhEThczUnlGRL6WIcxNqtrWR9SfBx60ccwSkXtsb+ZVETk4j3SvUtVWe9kNJKz/WwGxbiAODAXWYiqOl4DFwP5WZj4mwxXCLz2ORLoMcIiqvmD5zba91XvF7Fm3Ux78lqlqd4Y4XgPqrHsosMbq8APAtxyDOpyD2TNvVoEcY4HLahtvOgZEhxYXAD+zcr6qJmyaajE7ui8DOqxswTpMjyMLBkSHqXcrIoJpsJdgdqpfa99HFbCS3jzaLx1anCZmxH2JiGSqUwZKh9OB+637Psy7DOrwFUzDvcLKPA7MLJDfl4AKEXlERG4TkSGl5mc5fpDl1pFiRqJXiUhlgPd9tk69DzgYUx7/DXzBpqPgOjUFERkFbKeqTw8A1xVATETKMI3aBuv/OkZnYOsc6HkvE0SkKo+wDYrpXRTYQo/HZKo1wLEZ7s/G9iqAkZiGIWKPl4CRWcJtIvtI7c2AexZwj3UfCNxl3d+0YaQfwwLPVmNGXrukhT8dmIvZL248pgJcB3wReBs4wcrtDLzRT349caT5jwQeT3t/v7Huk4ArrPuiDNweSgtrNPAfYJS93gN4D3gVeJLeUfZ4TI8wFtSh5ffP/ugQ+Dqm8nk2y/0B0SGwI3CLdS9Jy6NdQDNwURrHgnSYKY6tpUPrdx5mtDkH04gJpnF7FtMBW8nmebRgHWIqoJA9bgVO3Yo6fAiYYt3bYSrblA47LMebAM/KVAKrCuT3YEAPZwGXDyC/WWw+UqsBKqz7UnpHpW/SW+e8AnRiyyOwuAR16ncI5P1ScrXpvhpYimkLjrPP1GEsQK8CH6T0au/9Fdg7U3oyHUVtEloIVHU5cJCITLEE788hPhV4RXt7m69gMu3aIpOx0J6XYXqtqOrvMWaljLA9ir8CP1fVRQH/PTBmgJlq3vxyETnK+v0RU+BWZAgS8uSXIY5C+B1u+V2Y6wEx37HuAs5U1TXW+1rgc6r6goicC3wf+KWqLheRLwJ/wmTMonWoqjcBN4nI2cD/AmcXwLEYHc7GmKyCaQnm0ceAz4jIX4N6L5DfFnHkgVLpEFW9VMx3s6uBWar6f7aHPhtTqR1I9jwKeXBU1Y2BdNwBHAncVgDHYnS4gV6LwhBgQ4Z6ZgRwFKZhL5ifjeMB634A+F0f3PrFLxO0d3QKcDu9I/6NGLPxchG5ADiCEpVHi5OBU/JMZkFcReQITGO8A0ZnC0TkAUy5/7uq/kpEDgCuwZjMC8ZWadREpFx7TSQtQGsGsVggPUuBPUQkYq93t36FYr2IVKtqu70ONgxi0/ZN4MsZnj0e01v/E6Y3ck/PgyI7YGzZn1fVdSl+qjpfRI4FFmC+96WG79MwveK8+aXHkZ44VV0rIjXp3hn4XYQZ7QURU9UjrDnjH8ClqvpM2rOpONcAO6TpMMHmOpyG6WEVpEMRqVDVLnvZjBlNpGNAdIgp5NcY6xxjReQaTE+8DJNHWzA94M4Ax4J0mCGO36nqd1I3B1KHqXerqioim4AOq8P5wKEiMgljYg3m0f7osE5Vm+3loZhRRDoGSoePYb6xvWjPT6S401vPtNNrQp6GGakVosMmoAFjLWnAjHRLyk9VN2TwR0SGqOomexl8t48B/4UZGR+D6UAfLCJ1wOpAEAXXqSIyDVDd/PNHybhamY2qmhSRVsyoMcSWdc6wwHNTyJyvMiPfIV0xB1CP+SYxD1PhZ5r0sC8mU96F+QZwJsZM8ixwegb5aZiMttGG+c0MMj/CNAoQGNpjPj425ZHuL2Bs8k32uMr634/J3Cn/71p+yzAFcCGwbyCc6zC9n0L4pcfxmQwyVwP1uqWp4SCs2asPfj/CZKRUHOdZ/+mYyq4JeBQYZ3W4xPKLW47bB/jt1Q8dXhCI+26gbmvpME1mieW3wOanZsz3ve8Vo8P0OLL4D5QOr7HX8zGdozLL8X3LcSPwrTR+/dHhpcAz9pmbgbKtWA4rMRM/FtjzAZbvGqvDVwmYCzETRX5UIL9yTOMxD9N4jCk1Pyt7g81zS+g16X0LM+lsPoHygRkRLQA2YUx1C4DDMN98g3m2IF3aZy4G/qeP/FaMLkPALRgT8XPAd6z/OOARq99ngMZAuHfn8w574ihE+KN2YD48FjRzZgDSMAL4ywCFPQW4YbDyczr86PPbhnQYAe7Ffl8bbPxsOv4BRAc4jq3KFbiMwAAhn0Psgw4ODg4ODh95FD2lX0Qmisg8EVkkIq+JyHczyJws5ofkV8T8HLlnsfE6ODg4ODiko+iRmoiMBcaqmSlXg/nW8lndfKbggcDrqrpRRI4GZqvq/lmCdHBwcHBw6BeKnv2oqisxs8JQ1VYReR0zZXNRQObJwCNPYz7+OTg4ODg4lBQl+aZm/8/6LWamUB0wSVVbAvd/Dcywl+OBWlUtzxLWGcAZAJWVlfXjx4/H83JbSX3fzy3jd+FrGC+Uuw3vMxwrs2TJknWqOjKnYBak85s4cWLe8WaVKTE/z/NYvHhxyTg6HeaRtm1IhwPCD5wO8wkjTcbpsB8oweyUEGb1jN2AF4B3SVt5IyA7A/OfyO35hF1fX6/z5s3TvtCnzMPTdd691xUfjpUBnk9Pa3+O+vr6guLNihLzU9WScnQ6zCNt26gOS8ZP1ekwnzDSZLZ1HSaTvi5b367z3litNy54R39y98t60g1P6dG/mKNf+8Oz+q3bF+oP/vqinvePl/WS+17TKx58Q69+9C29ccE7etecR0qmw+BRip+v97ON2pWYHyTLgeMImB+hZ3WMGzErGNxSgngdHBwctll0xpKsauli1aYuVrd09bjffLeLf6z6D6rmz2VTGYNizwqhjhiNjR82g160dydYuq6dt1ftzCPNNdz55xd4e00bS9e1053we+SGVJax3Yhq4klY1dJFVzxJV9ynK56k0x5qjYO/WT9vQNJalPnRLpo5D7Nu3ErMz5vHAvur6llWZgbmB9MdgOWYpVk+r6p3ZwmzZ9g8evTo+htvvJFoNJozHW1tbbllWhfT5g8jOmREceFYmZkzZy5U1cK2Q7BI53fHHXfkHW9WmRLzi0ajzJgxo2QcnQ7zSFsOfglfifsQT8KmtnaqqqpsZWjup0pwqpKMdXYwaUT/dVgqfr4qXQnoTCgdCejYtIwNiVo0UktHXHv8O+NKR0LptLKJZJKyUAhP6DlExLiBkAeHjklw1kkfjg6rqqtpi8HGbp+NXWqObmVjazPruitoSYbZ2GW4paMiBNVhJRTySO0lJYBIajkOw3Fspc+360unwytv/DMfbOigtrqSMg/CHpSFhDIPexj3mk3ttGglK9t9Vrb7rGr3WdmubOjqbScEZVSVx5hqj7HVwthq64561JQZXWV7l6pKQiGWhCl/uInd/nJ7v3WYDcWO1M7BmBP3wPzc+AZmMdiXReRMAFW9VkSewvwxXoUpd+dj/pDfAqp6PXA9QENDg0ajURr76LI0NTXllpk7m6bOk2hs/EJx4ViZYpDOzxu3K8teeZn9d9+DirIQlWWhzc7lYQ/Pk9xpKzG/vmT6QpDjnnvX65pEJdtN25uykGcPIRLyCFt3WcjjiQXzOfiQ6T0V2BYokKOq0p3w6Y77dCeSdCdMb7E74dP22gsl49fQ0KBPtI9i2fvdTJw4Cs8TW0HZSthWxgBLlr/HqLEjbe+1N01dcZ+u9cK6RIiySgL+SboSPkk/fSWizi3SFMTU+CYevfLYkvFbX7MDNz/3CkOHVuZ8bs36Trxyj9auBC2dcdpiCTbvM6c+nXT3+ERCHrWVYWoqyqipCDOqPMym5o0MqRtK0ld8VZK+klRI+j5xH3xfCZXH+80vnWN9Q4P+9Fmlu8ujuhpCIkZvnhDyzHUqT76/1qMl3kk8uflgwBMYUaZUR5LsMmkcY2orGDOkouc82p6j5eGtUg7TdfiujOGGRUsx1XMuSI9MTXmYqSNraJwUZftRUbYfWc3UJWfxnhzJ4cedUXT6m8oE/nJ7foQKQLGN2jcwa6u1YVYI3xGzVA6qei2Y/9iA7TEjuR0xqzzvk0/grV0JrvpPF3eueIEyT3oqwbDnEbaVYdgTVn4Q4/2K9xhaVcawqghDqyMMrYpQV1VGRVmoSIoDi2//+QVauxLw3FNZZcrDHmHxqX5ibk9DV17mUR4OUVHmUd58PC06jHva/xO4b89hr8fv7ffjvPvEUlPZJ0xlH0u5beXftqG7pGaP5o4YFz/dBU8/3rfwQ71rzqZ6rp6tYNDvAx6hZx/oafh6ZDyxvV2hsztGcu6/iQVMIun4aXP2d90f3LXwAzq7E3grl6FqRigpc5IfuA57ULXyA6OPMo+KsOm4VJR5REMxCCUZN67W+ns991Jy5WUeLbf+kZFf+yqCIZ/i3dPbF+i69uqS8kv4Pt1J6Iwnc8p5ApOGVfU0ULUVYWorjbumooyaRefwVrKRxsO/2COTqXyaCvGTOeMqtnOZjvrJQ1m5qpvhI2vRVEPqG92lGlaAacM89txxCmNqyzdrrEZGywnPO9R0vGZ+rqRpKwW++qnt2OueGxl2znnEkj7dtlPVnfBtHWCuV7z3Dkd9am92GBllZE35lh3MVetZkbtP9aGj2EYtAozF7IWzDLPo5gmYbTNSSGA2luvALDxaKyK7aOaVzzcbNg8dPYEVrQmWt60i6WN6a0qg59br948lr2ZMYHkIouHvUxX2GPLGA0QjEC0TaiJCtT1Hy4RoBNraO3ntrkfoSBhzSEccYxaJW/NJXDlodKZtzfJHmlmAy/cOE73n36w85lhiPnQnIZ5UYkmI+ebcnYSO7iSEksT9BPEkxGPQ2aW0JCEen0i3H+a9N1caM1XAXLWFcfm13tfuCZubH0IwLOIXXWEEOY4cPZZfrZnH6iOPIulDwuov0eOGhCqdXTEikUiPWW2Lc6yTmFYQziijhqcqfkKpqghR5oUoCxleEY8ed5kHO92/smT8Ro8ezR2HRGhrixGNVuR8zphkgpN+lZ6twVrLaPMhOqQl06NmF7tuaPtMPdH2d3LHM/OgonSYzm90+zt8b7ck0Wgs53NtbUmi0bQtDpOYbm87JKs+xSi/hmWvPd9HOG19pr+tra+tFHMjneNxo5tpq04QjWZ5/z3xJohWrjKrn66DTevMAoxvAnSeRJs/LK+0l0ImF9L5LX7xGfS4Q+l+/xWgd/+Z9NW0JwzrJvb+qyx6P0vAJeY4IOhrJglm0eBXMxzHAc1W5hjMBBCld0HVi4H/CoRzBaZhuxc4PJ9ZLPnOnHv4kUd19aZOfWNliz65ZJ3OeXmF/unpd/XqR9/Si+97Tb//65/rzMvv0v+6+nE95BeP6m4XPqCTf3x/XseUc+7X3S98QD91+SN69G/m6xV3PPyRmXXl+752x5Pa0hnTta1deve/H9ENbd3a3h3XeCKZM55ScnQz5/JI28dw5ly/w1Gnw0KOj4MOg0efIzVV/XS2eyKyWkTGquocEXkVM/HkUvtc+j5S7wP3YFYff4YSIuwJo2orGFWbpac891vWLPD5Hq940qe5I87Gjhgb22Ns7Iix4arfscM5/9tj46+tCFMdCeN5vUPwUps9BhIiQiQsRMIeNcDQCo+h1ZE+n3NwcHD4qKLP2Y8iMhcYk+HWeZgGaj3mX7UTgTmqmnGTRxF5FrMp3LmaZeajlesZNgOfsOFvsZ9YGkZsZZlqLcFPnxh+b5YobaV+B5NLyNHpML+0bRM6HEB+pZZxOsyCj4MOs6KYYR5mX59FmOlYj2K3JsdspndjQG4HzJSaH/Yjjj6Hp9uazIfBcbDz+zhwHOz8Pg4cBzu/bZ2javE/X+9rz5NUNbjF+/OYmZGImT7zU+BaVb2yyPgcHBwcHByyotitZ67GTKB5WEReFJHUNP5xIpKan/0p4FTMFvIv2uOYIuN1cHBwcHDYAkWN1FR1hyz+KzAzIlHVxwHJJJcnrv8IyhSKUsQ72PltbZlC4XS47ckUCqfDbU+mYLidrx0cHBwcBg2K3vnawcHBwcFhW4Fr1BwcHBwcBg1co+bg4ODgMGjgGjUHBwcHh0ED16g5ODg4OAwauEbNwcHBwWHQwDVqDg4ODg6DBq5Rc3BwcHAYNHCNmoODg4PDoIFr1BwcHBwcBg1co+bg4ODgMGiwVRs1EZkmInEROSjDvc+KyKQCwvqciLwuIl05ZCaLSL8XzRSRGhF5SkSaReSUgP8sEVkqIk32GG/9vyEivoi8JCK3i0hYRCIicq9NbyH8fiMiT9vjnCwynxSRi4rgd5qIPCsi80XkDhEpt/5TRORREXlCRH4SkD9XRF62HK+yfruIyG+su1Ad/sDG/YSI3CoiZRlkitXhkfYdPiYic0RkeNr9P9qNcBGRYSLyjOX3H5tXh6Y4FsovUxwZ7g2IDgP3m0TkxsD1LBF5wXL8mvU7SkR+aN2F6nC2LYepshDKIDMgOrT6ul9EFojIVXabK0TkGhFREXkuEMYIEflzP/iFROQKEZlr+e0yAPyy1TNTrV6bRGSeiEyw/kcHyuFDafWM1w+Ovwjob6WI/M9W5JitLr014LdRRGYGOeaMbCA2acuxKdxtwFzgoAz3bsnknyOs4UAFsCSHzFXAPkWkN4zZ9Xs2cErAfxZwfgb59cDjmB3BbwGOtv7nAo8UyG9He/aAp4DtM8j8AxhVBL+pQMi6fwF83brvAA627rnATpgtht4C/mS5LMPsWgtwHzCyHzqMBNy3Ap8ZAB1OAsqt+1vAJYF7ewD/BOZmyKOnAw8E/O8D/lIIv2xxbA0d2utjbbpvtNdDgReA220+fTkgOx8o64cONysbWWQGRIfA5cCp1n0zcJR13wU8ATyXFs51VheF8PsmcMYA88tWz1wBfMW6ZwE/t+7ngbttPp3D5vXMcYXqMC0trwDjtiLHWWSoSwP3y4GlAf2fCxyXK66tNlITkf2BVcAHGe7tAhwFXCUid1q/C23L/oyIfCb9GVVdr6pZR2kWh6jqCza82WJGT/eK2dNtp77SrKoJVV2V5fZpIvK4iFxie0f7AysxDRvAECC1ceobwIEF8nvLnn0gASSD90WkBhipqmvs9S0icoOI/Mv2akflwe8dVU2F223jAdhLVRdY97+A6ZjdzVsxW7SvBGJA3MrMBc6kcB3GrJxgGu8lGZJZrA6XqWp3Bo4AFwCXpS7S8ugRmAY8hVcw2ynlzS9THEEMpA5tb/bbwDUB8f2AxcAKzLuukt6R3ULMvocF6dDibFsWvpPl/kDpcDpwv3XfB0y3OlyKyaPpeAU4rEB+XwQm25HS1SISGQB+2eqZ14A66x4KrLHu1UALJp9W01vPzMFsztwfHSIi+wCr1WwdtrU4QlpdmnbvM8AjAf3PwegkK4rd+boQnAd8Fdhi92tVXSQiD2B6lI+LyF7AwZiGYAjwrIj821bweUFERmIq4SDWqurJInISRvk/EmP6mZ4mF1PVI3IE/09Mjx7gD8DJmBd9AaYXvDfwtJodwAEewDQK/1MoPxE5GXhHVd9Nu7UTkO73mqqeLsZk+CXgahG5DvhEmtwyVT0tEMdOmIJwsPUKZqxmYIyqJkSkAviyTfPvU40SpqI8yvIsSIcich6mt/YW8H7avZLpUERGA2cBR9rrRpvu1YFnUnn0VzbdXw3cWwB8lgJ0mCWOIAZSh1/B9OaDHb/hmLz5bUw5bAWGYRqAxZZHoTq8CrgIYzW5T0T+E+gQDagObdqbrbvZXqd0eB1b4hFgA4WVw/HASlWdISJXAF8Drh0IfhkwF3hQRL6OGbHsZ/2HAvVAJbAwUM8sxoxq+1uXnoKpvzbDh1CX3ha4fwomj6WwGNg9R3hbp1GzvYPnVXW96ZT3iU9gGgUFmkVkDTCC3p5Kf7HQnpcBhwOo6oWFBqKqG1NuEbkD+G+MSeCXmJ7j74BTROSLqnpnhiDy4icin8YU0Jl5Ji3Ib3ub1v/O9YAYO/0fgRMCI99ghh8CbBCR04FaYApwI3C0iPxKVZdnCTovjqp6qYhchtlFfRbwfwVwzEuHIlKLMUmdmRoVAecAJ2B7wml5dBymcWnLEWw+/DaLI08UrUPb+TgZ08gFv19vD2wIlMMoppLPhj45qmrKMtEpIncDDZgOQL4ci9HhRkz+bLbnoZS+ntmAaSSw5+PzCLck9Qzwc4xp7m4RORG4TETm2HTvBVwC1OWoZyD/uiaEMV3OzjNtA1WXHolt1ESkDtOANRUS5tYaqe0FNIrIgZhE7iQiX1bV9wIysUB6FgOni8mZQ4BRGLNX3lDVtWLMO5t5B9ypj8oF9y5EpE5Vm+3loTasRmAsphKZBDyL6TkCTMOYJfPmZ80ol2Ds5Z0ZkvEGpoHpi1/WXr6IjAD+jqko3g7cf0lEDlTVJ4Gjge8BZwBVmJ7V7phvbNOA5fb8KrBvgRwrVLVLVVVENgEdm5EpgQ5FpBLz3epSVX3Gytdg7Pt3YHq7uwLfBcI2j34KeE9EJgfy6DSb/rz4ZYpDRM5T1UsDaRwoHW6HaUjvx+TBsSLyDSBi0/EgsCdGn2OA9yy/OcDO+XK0aatT1WYr04j5ntNLZoB0aPEYxiT8Z3vuoLee2QuoLEaHFk2YhnpJ4FxSfmSHBNKzBqPLvTCj4puBXWxcqXeSKof9qUsPw3QIWtJvbOW69M3A7S8Bd9sGOYUUx+zQfn746+9Blo+YwH9hPlZfZ68vBp7GNA4zM8gfjBmed9jz8RlkrgbqrXs29gMldiJHnum9D3jHvshrrd+lmIz0BCZzlVn/L2Ls2y8C9wJR7f24eWmB/F61R5M96jPI3AOMTn+vmCH77Dy4XY2xy6fiSE0UmQrMs/zOD8j/zKZ5LfCntHc0qh86vMbGOz/4HkupQ+BHmEKc4nhe2v0pbD5RJIoxtUzPkA9OLYRftji2hg4D9xuxE0Xs9deAJzEm0eCkkgWYRq9QHd6Cmcj0NHB5jjSWXIcYc+q/bNqvATzrf5bltwlTN2xv/a/DdF4K4TfU6qgJY86t3or1zK6YCT1NmLK4W6CeeRbz/fdxNq9nPluoDq3MrdnuDTDHjHWpvTcf2D0tjHOBz+aMJ5/EfFQPTGVyw4echgimgfMGIOwDgIs/ZH67AL8d5DocMI7biA6PAn44yHU4AvjLIOY3YPXMR42jWGEHBwcHB4ePPIqe0i8iE8VMd10kIq+JyHczyJws5mfBV0TkSRHZs9h4HRwcHBwc0lH0SE1ExgJjVfUF+zFxIcbmuSggcyDwuqpuFJGjMd8K9i8qYgcHBwcHhzQUPftRVVdif3RU1VYReR3zb8eigMyTgUeeBiYUG6+Dg4ODg0M6SvJNTUSOAn6L+UGwDpikgamhIvJrYIa9HA/Uqmp5ejhW9gzM9HEqKyvrx48fj+fltpL6vp9bxu/C1zBeKHcb3mc4VmbJkiXrVHVkTsEsSOc3ceLEvOPNKlNifp7nsXjx4pJxdDrMI23bkA4HhB84HeYTRprMR02HCiR8UN+nLOyR62/BYnWYFSWYkRIC3gZ2w6wr9y6wSxbZGZhpqLfnE3Z9fb3OmzdP+0KfMg9P13n3Xld8OFYG8z9H0e+uvr6+oHizosT8VLWkHJ0O80jbNqrDkvFTdTpU1c5YQu9/6FFNJv284tkWdfjIP6/T99a167w3VuvNj7+jF9zzip58w9N64M8e0Snn3K+Tf9x77HLBv/Xgnz+qn73mcf36Lc/pj+96SX/xwOt604J39K45j5RMh8GjFD9f72cbtSsxa+WVY/5MXxQUEpE9MCtRrCDtB00HBwcHAN/v6QBv8+iMJdnQEWNje4yNHTE2rNqJ5zqq+c/Di2nuiLGhI27O7TGaO+JsaI/RGTfLdMqjc6gpD1NbWcaQwFFbUcaQqjLi6+M0fojcVJV1bTGWrmtn6bo23lnXztK17Sx976ss7RhO4sl5PbI15WGmjqxm3ylDmTpyIpOGVbHy+hvgy6eyrq2b9W3drG+Psby5k5c/aGZ9e4ykr/x2/bwcKeg/ijI/isgwzE+6O2O+q+2FWRl8f1U9y8rMwPy4twNm9YntgM+r6t1ZwuwZNo8ePbr+xhtvJBqN5kxHW1tbbpnWxbT5w4gOGVFcOFZm5syZC1W1IadgFqTzu+OOO/KON6tMiflFo1FmzJhRMo6DXYfX3vJnujs7qKutpswDL8sSTfnq0FclqZD0CZyNX1t7B9VVVTYdZimH1BnMuauzgzFD+6/DdH6/vOHPvLW2g7JIOT7g+2YdNV97D1WlozuGhCPEkxD3lbhvz0mI+5CIt9OVLMP3yognU/c3l0na6igk4Ik5hzwIiRD2ev2+NNXnR6eWSodj6vf8/s0kk0lCoVDPu029zxQUaO1O0pn0aI2Z9GdDVRhqIkJ1mVATEaJlQjQC0TLBT8RIehE64kp7QumIQ3tc6YgrHQnjnhRVLjiwdDr8/EW38vB7MSrCHuUhqAiLOYeE8jCUh4SKELR2xVkXC7Gq3aczsPR3WGBUtTAm0sHwijATh1UzptpjTJVHTQTSlyXLldd9NZyn3HITu//l9n7rMBuKHamdgzEn7oH5Me4NzAKqL4vImQCqeq2IPAWMwyzLo8D5mL/zt4CqXg9cD7DrHnvro6sjTK0a3/PSTAGWzTLd8g3vcNi0XZk8rIrxQyspC6XZhOfOpqnzJBobv5CTTFNTE42NjX3KFIMgv4aGBh39iX14/cnnmLrdTmzqjNNij02Bo607Qaw9zCcmD2VkTTkjouWMrClnZLScETXljFx4GQtjXygZv75kCuG48+576f0flDF2/Eh8X03l7Jsj4Su+Pa9Z08X4sXWEPSEc8igLCWHPnkNCeMVDvJ+MMq1uImFPKAt5hENCmWfO4ZBHmScsWvka08bvSCzpk0gq8aRvD+NOJH1GZl2usnB+DQ0NevaCbmJJIbXKVyTkUV7mUVEWojxszpGQR2tbiPJKIWHTk/B705iMj6TbL8OnHT9nP1Mwa2Nnx9T4Jh698tiS8VtXMYFbF7+FWX0pd9pE4lSEQ5SXeT3vodxel8eW4XsjGD1mDOVh6x8OyIQ9ImGPpUvfZcKkSfbdmPcUTyqJpE/SV+K+Mrwq15KVhXGsb2jQoUPr2Lixmbq6ut71nzR10hQ7Kto3seOksQytjjC0KsLQqrIe97CXv8GixDEc89mvEU6vgwLIp4zNfXReUeUwXYefn74niaYXGTZqHJ2xBO2xJB2xBB2xJG3dSdZ0JWjvTiJJj10mDuPgEdVMHVHNdiOjTB1Rzbi6SkKewNxGU5fOPCNn/HnVNRGBv2yxfnLRKLZROw74oXXPwizncguYxiwlpKrfAL4hIsuAh1X1SPJARzzJ3W/F4a3FfcreuuhZwPTkxtVVMnl4FZOGVTFpWDWTV09jjVdGS1ec2oot9qH8UPGla5+itTsBzy3s8RMxQ/ohVcYkES0P09ytNL25tmfovjm+S1iU0YseNY1ctJyRNZHeRs+e32/1ee7dDT2NZktnnJauRMAdp725myLbtM3QGU8y7/0EkVUfEPKEkIg5e4InpsEKidDR4bMytmGziixV4Sd8Jel/ygS47I2+I33pxZy3z9wz4xylfuOi43bl1UVvMnG7qXTFk3TFfboT9hxP0p0w1+uT7YwdHSXsbd4gl4U8wh/cyYrkTkydtnevX1oDH/aEDTfdyIjTT8dXBTUVrqZGS9Ydv+WGkvKbdeAU9r3jt0z+5S/xPCFsdZfSp+dByBOefHwBh81o3KLX3oO5F9gKMfeawE1NK2hszL2bSbGdyyAEuOOMA2xFfECf8TY2ZvnN9u0NfNDp52zQ8kXYy2tB5rxx6E6j8VaV09iYc4F7y++j/bdVsY3aZMxGdtthGrMxmFXJT0oJiMhEzLpiEzFT+R/YIpQA0obNXHegUh2tRgHt6TkR6EXBuk3ttEslazp81nYqazpirFzbzcvL1tMaA9P2wuyXHiJaBiOrPEZVCqOqPEZW2XOl0NnRwe33PUpHQs0Rh85Er1mgI6HsPTROMUjnd+5uIcbMeYDOzx5LdZlQVSZUhoMmrCSQpG1ckmi0El9DtMWhpVvZ1K1siimbWtawPlZFB3E2dXTzZjM82+3TGtt81VEAnnhqizRVhKCqTKgKw/Byv+gKI53j1Qco0WjuhqStzScazV4ZaMtiNiWHUVkznISmzHJpZjqF0XfdyYYvfZGQQNhWtmFrwgpbk1Z7e3tJ+Y3teIeaEd1E9X1TorKUqrbhCaLR9B08LCaOo80PEy0PbAOmmJ3DAmagqmP3Jdqaadu5QDzHfrIoHabze+m5J2n73NEkXn4253NdHe089thj2QU6T6LNH9Zn2tra2vKSKQbpHJuamvKON6tMifmVUocl4QfblA6zoq+ZJJgFQV/NcBwHtGGK8DGYLRp8ehcbvRizsOZYYB9M49eBWTU64+zITLN2ip2x09IZ01fv/qL+8g9/1Gublui5d7+sp9z4tB7880d16rn/2mymTq5jx5/M0fpLHtaf3f7wR2bWVTyR1NUtnfra8k3a9OYa/cMZP9H5i9foS+9v1KVr23R9W7fGEsmM8ZSSo5s5l0fa3OzH/MNRp8NCjo+DDoNHnyM1Vf10tnsishyzc+8cEfkR8Bu1W2uo6k8DoitF5GDM7MjRpP2cPZCoqShj15o1rA130Th9+83uxZM+K5u7eG9DO8s2dPDW4sU07LkrtRVl1FSYmUkpd0WZ+YBcSrPHQCMc8hhVU8GomgoA9MTDOXjH0v4S4uDg4LAtodjZj78E1qvq5XYlkbdVNeMXahF5GvgN8AvMFgpb7Ntj5XqGzZg9pNbT915qI7ayTLWW4KdPDL83S5S2Ur+DySXk6HSYX9q2CR0OIL9SyzgdZsHHQYdZ0ddQjtzmx+GYLdLXYTexs880sPkeTlMw/6ctJMO+Z33E3+fwdFuTKfQoRbyDnd/HgeNg5/dx4DjY+W3rHFWLND8CiMhtmKn6h6lqh33meeAbAbHlwCvAg5rl/zQHBwcHB4diUdTsR7vm49mYXYI7ssgIcBNmlf5fFROfg4ODg4NDLhT7Q8XVQA3wsIi8KCLXAojIOBGZY2U+BZwKHGplXhSRYwqI4/qPoEyhKEW8g53f1pYpFE6H255MoXA63PZkCobb+drBwcHBYdCg+F/fHRwcHBwcthG4Rs3BwcHBYdDANWoODg4ODoMGrlFzcHBwcBg0cI2ag4ODg8OggWvUHBwcHBwGDVyj5uDg4OAwaOAaNQcHBweHQQPXqDk4ODg4DBq4Rs3BwcHBYdDANWoODg4ODoMGW71RE5FOEWmyx9cz3G8UkT0KCO9AEXlFRLpEZEIWmRoRubOINIuIXC0iT4nIcyJyYiCtKy2XpIg8b93/EpGGInmdLSLPiMgTInKV3e0gXWayiPR7UVD7Xp4SkWYROaWvuEUkbt/zxpTuRGSEiPy5nxxPEJHHRWS+iNwvIrVZ0thv3dkwHhSRtSJyfsCvSkTusvr6h4jUBfLm+yKyyLp3sfJ3i8hRhfALxHWRiCzJcq9YHR4pIk+LyGMiMkdEhlv/z4nI6yLSlSYfE5FWEdkkIg9Yv2J0OEtElgbK9PgMMsWWv71tXpwvIo+KyNTAvR+LyFwb96FWhy9Yfm+KyP9auYiI3CsiMwrVYXocW5Ofvb9Z/rFlcJM9brJ+xejwBwH9LRWRKz9sjtavTETeSpVbEdlFRH7TZ2QDsUlbHxvDLenj/mzglALCGwJEgSZgQhaZH1Lg5qRpz+8GzLPuGswO3wCN2M1Qg7yACcDfi+S1Y8D9N8x+dekyVwH7FMErDIxJT1u2uIGlwFRgblo41wF79oNjJOC+GPh2qXUX0Mcs4PyA3/eAc6z7y8ClwBLgm8AZGcI4Efh3Ifzsc6OBv2TL9yXQ4SSg3Lq/BVxi3cOBivR4gTgQte4mYOcidbjZe80iU2z5GwPUWPcxwG3WfTRwWZrsEuA5YJK9/hcwzbrPtboohN8WcWwtftnyD9AJbAdEgBeAocXoMC0tc4BPftgcrf9ZwH1p5fY+YGSuuD4M8+MY26u8W0SmBG+IyDBMITnP9hpCInK6mFHDMyLytfTAVHWTqrb1EefngQdtHLNE5B4b/6sicnAeaV4BxESkDNOobQjcO1JEFgCTRGSBiNyNaSwmiEhVEbzeClx2A4kM6TpEVV+wccwWkdttb/RFEdmpL1KqmlDVVQXEPRL4M7BHmu7mYLYXKpRjLHBZDbyWIZnF6g5V/SCD9zTgeet+FpiBKXiXAT8QM4K9WkQiVuYpK5M3P4sLgJ/lSF6xOlymqt32skdXqrpeVbsyPCKYraL+gcnLzda/Xzq0OE3MiPsSEclUpxSlQ1Vdpaqt6RyBLwEVIvKIiNwmIkMwOtwV+I3No89jOp8AjwMzC+SXKY6txQ/S8o+IlGMas1uAO4AXgf3s7WJ0iIiMArZT1ac/TI42jiimQ/H3tKDmAsf1FdlWPYAR9nwk8EiG+7OxvQxMJfoSRokR687YSpN7pPZmwD0LuMe6DwTusu5v2jDSj2GYiuBqzEhlDXCcfaYGqLDuX1nlHAk8AvwV2LsEvKZbRUqa/0jg8bTwf2PdJwFXWPdFGTg9lO2d54obGAFMwfQOHwnI7Qz8sz8cga9jdkV/Nsv9onSX9mywx/fNwDv6NvCG5fcmcI3V4RXAmYFn1hfCD9gRuMW6txiplViHo4H/AKPS/NN7v2cCq+zxXrE6BIYCIXvcCpw6gDqsxnQudrHXDwbe0VnA5VaHT1qejwJP0Dsir7S8C+G3RRxbkd8W+QcYBzwRqEPfAU4sthxa2e8AFw1EHVoIx0CeP5wty+3RwG8zpTF1FLXzdX+gquvs+UERuaYP8anAK2p79CLyCmbYvbbIZCy052UYMw2q+nvg95mEReQIYDywA8bcuUBEHtDengfAjcDPVPUSy2tdjvjz4iXGLn45MFOtRgvgdbjldWEez22BTHGr6jrbg9qAadxyIS+OqnoTcJOInA38L2Yn9VwoSHc5cBPwKxGZhylkKyy/DcCvgQcwvcfjszyfD7/ZwE8LTFfBOhTzLfIuTAO8JodcDfAj4BNAG9AiIvup6rNZHumTo6puDIR/B6aivS1XeumHDq2V5K/Az1V1kfXegNET9vw7q8PTgSuB/YG7MZaWfvHLFEcf3ErJbzZb5p8NmIYhVYcOY3PLUToKqUNPBk7J4J8JA8ZRREZjBgQXisisPNPTg63aqNkKsVNVk7bSzFTxxwLpWooxdaVMQLtbv0KxXkSqVbXdXgcbiNQkiG9ivq2k43grs9GmuxXT4wmJyBBV3WR5HQa8GeA1BdPr7xcvEdkBuBn4fKojEISqrrWV1GbeGXhdhBlxBRFT1SMycM0ad0p3ViTK5rqbBrzaD44V2msiawaqMiSnKN2pasYCbwv5Wfb5M4C1IhLC9CyPt/waMN9oEJE6oLUQfpgK5Rox82zGisjvVPU7gTQUrUMRqQT+AVyqqs9k4hpAJUZHbRgTXQdmpAX912Gdqjbby0PZPM+nUGz5awb+hBkd3BO414TR0Vx7fk9EQqr6mu0k3YCxtPw7wHFVIfwyxJFpws9A8dsi/wA/AdpFZBJmVBrGWDlS/ArWoU3fNEB1808PHxbHe4GRYiYyjQfKReQlVb0vwDE7cg3jSn1gbL//AeYDjwF7ZpDZF2MyuAszO/NMjNKeBU7PID8Nk+E2AguAb2aQ+RGmgobAcBYzgaApj3SHMDbsxzEfob9j/b+Fsdn/xyrtCcvrcODuInndjylATfb4TAaZq4F6655Nr8nhIOyQPg9u92FMGK8C12aL2+rufWATZrLB08D2Vv46YK9+cLwgEMfdQF2pdWdlb8B8r1tCr9lkFxvvI5iR2QFWj09iGrTnbJqqrfyJmN5/3vzS0pBtokhROrTvZ13gPZ5n/Q/GlIsOez4+oMMWq8d/AqEidXgp8Ix95magbADK3xcwDXGK41XWvxxj8pwHPAQcZXX4NqY8PgscEwjnXJuWQvilxzFma/HLlH+sDt+w+tsEXBiQ6ZcO7bMXA//TRz7bKhzT/HriCtRXo3LGlU+l8FE/gFrSZiMOcHyXAftuhXimADd8yO92BPCXwaK7HOn4B3bWoNPhR0+HGOvKvYA3SPl9HHS4C318T1PVngkADg4ODg4OH3kUPaVfRCaKyDwxP6u+JiLfzSBzsoi8LOYn6SdFZM9i43VwcHBwcEhH0SM1ERkLjFXVF+xH74XAZ7V3dgsiciDwuqpuFJGjgdmqun9RETs4ODg4OKSh6NmPqroSWGndrSLyOmbGyqKAzJOBR57GfFx0cHBwcHAoKUryTU1EjgJ+i5kpVIdZoqYlcP/XmNUYwDR4tapaXnTEDg4ODg4OAZTC/BgCFmOWLrkV81/IMUHzY0B2BmaNr0dU9eQs4Z0BnAFQWVlZP378eDwv96c/3/dzy/hd+BrGC+UemPYZjpVZsmTJOlUdmVMwC9L5TZw4Me94s8qUmJ/neSxevLhkHJ0O80jbNqTDAeEHH0sdKqC6+dn3fUKeh4j5wcvbYqnyD1eH4nkkffAVkqokld5rP0lCBV/F3t/8h7VCMKoSPljafx1mQykatQMwS5ooZkmZcgBVTV/Law/MtOhNwI9V9eG+wm5oaNArrriCxsbGnHJNTU25ZeY20tR5Eo0zzyguHCszY8aMharakFMwDzQ0NOjzzz+fd7xZZUrMr7GxEREpGcePuw4TSZ/OeJInn3icw2c04mWqxbZRHTY0NOgj85/iwaYF7Lvf/viq+L7aCk7NtZrrhQsXUl9fbytqIbWvRMotz3yN57uOYf8ZX7T3BE9ARDZ75smnnqZ+332JJ5Wkr8STvj0rCd8n4Sut777GcUcdWhId1jc06EU338tLr7zGDtM+QSzh053wiSV9YonAkfR59/3lDB0xiq54ku6ET1c8SVfcnGMt79GcqIWyKrrjSboSJt35IBLyqCjzqIyEqCwLMSbSzR3fPapkOvzR/93Nnx5bRM2QOhK+TyypJJI+8aRPIqnE7Lm5vYuuZJY0hj1GhJqJhCuYMmE8w6ojjIiWM7w6wvBoOcOjEYZXRxhaFeGpp59mt70a6Igl6Igl6YglaO9O9ly3x5JMjH/Alz5zWEl0GERRjZpdomUeZs2xlZgf/44F9lfV1GoNMzA/mO4ALMcs0fJ5Vb07S5g9PYxRo8fUX3HNDQyJVhP26Dm8tF1Y2traiEaj2RPaupg2fxjRISNy8ukzHCszc+bMfisiyG/06NH1J/6/W+nsjlFRHsETCIvhFxIIeaYXFxIo87sZW1dJXblHTSTtHZSYXzQaLarSD3IcMXpc/WkXXUekfHNrc3q1nox3M6S6gqowVIaFqjKh0rorw+C1vbXN6vCwn/yRtu44CUJ0J6EroT3nriTE/c2fDwmUefYIiTkTQyRERVmYshCUeUI4JeOJ8RPQZJzqighlHoQ9+2yavJfoYrex/ddhOr8vX3Ir/3w73p9XNWC4aG0Ts26+okQ6HFNfMevGrLIh6a17ykQpD3s9uosEdBjxWxGJUFVRQVkIIim9eUIkZHQT8YThj85l9YzDiCWhOwmxpBJLQszXnush4QQn7lo6HX798tt4Ynk3kXDI1i3Sw8ucbZ3jxxkejVATEWpTR7lQExEqQiDbUDnMhmInipyDWXZmD8zPjW8AXcDLInImgKpeKyJPYRbirMKM6M7HrNawBVT1euB6gO123kMveE4wiyL0IuSJyVAhj0jYI+R7TBtfyaRhVb3HcHOuqSiDubNNL7jxCznJ5NsLLgZBfg0NDTp/pbChXVD6qjQE82oN/5HRckbXljOypoLRLWvpCNXwyXFTGVVTwajackbXVjCsKrLZqKAvfklfmfto3++gEI6TP7G7/n2pYFbu6Ytfd9a70dAIIqEQI4YKNRVl1FSEA+cwtda9+L23mFQ7ka64T3ciSXfc9Lq7E6Zn3R33qY9m6Yr2g19DQ4OuipfT3Z1k5LBahpWHqYqEqI6EqSoPUV0epjoSprIsxOIlS5gwabvetCSSPaOC7uVNrEhMIjpsZI9fayJJd8zvlYkn6YwLvubOK1PjmzjrymNLxu+Mz+zPLhf/ilH/fTqeiOl0eWaUFbxe/atfMeoH3zdmNgVftdfspoq+fCGvxg5il/pP41s/sHKBZzbeeANjvnmmaahDHmFPCIfElnvPnC+8p9/8MnG8/XuHsPa8c9nuyl/21CvlYY9IyMu/DOU72h4X5vRiLRd9IJ3fpbMOL94iBDD3om2mLs2GYhu14zD77IBZzqQJs5wUqnptSkhVvwF8Q0SWAQ+r6pH5BB4tD/Pfm56n9rSvEEv4xJNqz+botu633ltOS2ecf72ykuaOzQv80KoyJoVPoSIyhOcibzBpWBUTbcM3dkgloUymoK2IhReYzHbIIdNJ+Nbc4vskk0rCN+aWWMLn1XMuIPTds1nT2sXqli7WtHSzurWbDzZ28ML6aWyIV/GPZa9sFnbYE0bWlDOqtoLRNeV0bOrmzuUv0B5L0N6doK07ZRYwpoHOeJLRVcKRny4dvxHRci4+rIqDDjqo1zODcWDeggXsvs9+tHYlaO2Kp50TtC7+C29170z1iCit3XE2tMd4b30HrV1xWroSxBKB4dDrbwBmlFtRFqI87FEeDlFe5lERDrFLRWkXHHjkh422EH8qp1yTv4zGxh0z35z7fVshfj53GE1NHHzIdNvQmcYx5e6yDfi6Sy7CrE1bGuw6bghrT5hB4965Jy03nXQUjTuNzi6wejFVnQ007jEudzhth9K4Zx8yXzkZbsk+uioUnxhTw8pTv8D4usqShenw4aBY82MX8C7GpNgEfNJen6Sqr1mZiZgJJBMxi1feqKpZuzLpw+Ybb7wxb5MZQHtcWdfps6ZDWdvhs6ZTWdvaxurOCBtjHskA3ZDAiEphZJXHqEoBP04oXIav9HwETX0s9X3wgU+OSPCjU0tnurrjjjsKMgtmROtiNiaGkawYRnO30tyl5pw6upTmbp+2mE9lmUdFSKgIQ0VYKA+Zc0XqrDGO2rF05sf+6DAbx1xmj7ivdMaho6OduppqyjyydlhKbX4slQ63FRPygPCDbfozgNNhnmnbhnSYFXmstzUXs9ht+nEcZoHKMGYn0w2Yej+1oOrFwH8BY4F9MHtTdWBmSu6Sz1pf9fX1Om/ePO0Lfco8PF3n3XudxhNJXba+XR9/a63++Zn39GdzXtdv/WmhHvu7BbrH7Ad12k/u190ufED3vOhBrb/kId3v0of1wJ89ogf9/BFt/OU8PezKJv3FXx5W4Pl80p8Pv7zS35eM5VdUGGkypeRYSh0WHY6VcTrcyvxUnQ7zCSNNxumw8KNP86OqZjVGichyzIZzc0TkR5gNDi+1zwX3yFkpZnfUP2E2M9zs5+ythXDIY6I1P2YyFH2YdmAHBwcHh+JRrPnxl8B6Vb3criTytqpm/EItIk8DvwF+AeymgZ+z0+R6hs2YzQzXk3vDTTArVG9NmWotwf8xGH5vlihtpX4Hk0vI0ekwv7RtEzocQH6llnE6zIKPgw6zoq+hHLnNj8Mx+1GtA9Zgt+3GbKZ3YyCMKZjdZxdiNm7MeyhJHsPTbU2m0KMU8Q52fh8HjoOd38eB42Dnt61zVC3S/AggIrdhpuofpqod9pnngW8ExJYDrwAPapb/0xwcHBwcHIpFUVP67ZqPZwPTUw1aBhkBbsKs0v+rYuJzcHBwcHDIhWL3U7saqAEeFpEXReRaABEZJyJzrMyngFOBQ63MiyJyTAFxXP8RlCkUpYh3sPPb2jKFwulw25MpFE6H255MwXA7Xzs4ODg4DBoUvfO1g4ODg4PDtgLXqDk4ODg4DBq4Rs3BwcHBYdDANWoODg4ODoMGrlFzcHBwcBg0cI2ag4ODg8OggWvUHBwcHBwGDVyj5uDg4OAwaOAaNQcHBweHQQPXqDk4ODg4DBq4Rs3BwcHBYdBgqzdqIlIvIg+JyDwR+UWG+40iskcB4R0oIq+ISJeITMgiUyMidxaR5loReVJEmkTkWRE5zPoPtVweE5EnUukWkZtF5OkScjxbRJ6xcVxldz5Il5ksIv1eIFREIiLyNxFZYOM61PpPFZH5lvtzlus8EXk5/X33g9cJIvK4Df9+EanNIFOU7mwYD4rIWhE5P+A33b7PFJ+JNm++JyKbRGSFiJxjZSMicq+IzCiEXyCui0RkSZZ7xeqtRkSeEpFmETkl4D9aRB6w3P4oIuXW/1Mi8oGIbBCRZTYPp/h5/dDhLBFZavNHk4iMz5LGYsrf3lZX80XkURGZGrj3YxGZa+M+1OrwXfs+VorZyLgoHabHMQD8stUv24vIQhFpE5GDrF+9fQfrROT9YH0gIneLyFEF6u8HAd0tFZErtyK/Lcqg9R9m64MFGfhV9xnhQGzSlmNTuAjwMFCTQ2Y2cEoBYQ4BokATMCGLzA8pcHPStOc9IGzdU4HnrPss4ELrbgT+ajnOB/5ZQo47Btx/w+xdly5zFbBPERyPBf5g3VOAZ6z7CuArlterwK+t/wEpdxG8IgH3xcC3S607G8YEYBZwfpa4vwZcafPmXgGdPwVsb6/PBf5SCD/73Gj73JIs94vVWxgYk/7uMbvMn2DdPwZOtzpcChyXIZxzMRv/FqrDzd5rFpliy9+YVJ0BHAPcZt1HA5cFdWp1OCzg9xiwa391mB7HAPHLVr9UAcOAW4CDAvx+BZxqZW4GjrLuE4F/F5pHA+mYA3xyK/JLL4O/tO7Lc/D7bl/xbe2R2gFAG/Bn29s4OHhTRIZhCsl5tlUPicjpYkYOz4jI19IDVNVNqtrWR7yfBx60ccwSkXtsq/9qehoyQVV9VU3Yy1rgZet+3V4DDMXs/n0AsB44xHIoBce3ApfdQCJdBjhEVV+wccwWkdttz/RFEdmpL47A20C57RWluAC8BtRZXj5whIg8iqlMDyuSVyxwWW3jSkdRurPxfNBH3LVAFyZvXmL5fQrznpNW5nFgZiH8LC4AfpYjeUXpTVUTqroqw61pwPPW/SwwA6PDOuBXIrJRRP4QkJ8DnEyBOrQ4TcyI+xIRyVSnFFv+Vqlqq70M5v8vARUi8oiYzYo/jdHhH2390mivV1j5/uhwszhEZMgA8MtYv6hqh6puCIim6s9ZwOk27PuA6fb+Uxg9F6o/RGQUsJ2qPr0V+aWXwVS9Oh2437qD/P4NfKGv+PrV8hbRYp8IvGMJjAfewG5/E5CZje1pACOBlzA9lIh1j8wSdhPZR2pvBtyzgHus+0DgLuv+pg0j/Rhm74/HFIo1wLHWrw54AjOC+QAzwklx/DtwRIk5TgfmZghvJPB4Wvi/se6TgCus+6IM/B6y98qBfwBvAquB/az/RGAR8B4Qw4x6Urp7BhhaDC/g65hd0Z/Ncr9o3QWePT8t7M9gKv7FwPfYPG+uAP4YkK0EVhXCD9gRuMW6txiplUJvmd69vf45cJZ1/xJ4AJM3Ffii5dhCby+40nIoSIeYDlDIHrdie9gDpMNqTMW9i71+MPCOzgL+GdDhzUAcs0GxV4QO0+O4fGvVL4Ewb8GM1FJ1y1v0lsEZwHUB2fWF8As89x3goq1Zf2Yogzuk4qN3W7R0foszpTF4FLXzdT+wAXhSVVuAFhFZh3nxa7LITwVeUduii8grwHbA2iLTsdCelwHDAVT198Dvsz2gqsuBg0RkCkZZ92N2/f67qv5KRA4ArgF+BzwJlNl0loSjGDv55cBMtdotgOPhlsOFOeS/Aryvqp+zHO8G9sFUjucD7fZ8rqp+2+quMkd4efFS1ZuAm0TkbOB/Me80X1556S4bVPVfwL9E5EvAt7F5U0T2w5h+LsjxeD78ZgM/LTBZheotGy4DrhaR4zGV2QpM+esC/qGqCRFZjen9P5AljD45qurGlFtE7gCOBG7rI20F61BEyjDm/Z+r6iLrvSGQ9geA0+itX74mItOAccBRmJFowfwyxPG7Prj1i1+W+iUdGzB1yw6Y8rgOmGT9M6GQ+vNk4JQM/plQMn5pZfAyzMh4I+azUrM9Z+OXEVu7UXsGY94JYyrEUZieRRCxQLqWAnuISMRe7279CsV6EalW1XZ7HWwUUh8hvwl8OcOzxwPtqtptr1uA1sCz66x7Dcb+/QxwiY1jOSXgKCI7YHqen1fVdaRBVdeKSE26dwaOF9E7lO9Ji6oekcZlI2ZH8yDHlzGN83Ab1yigO1ip9YNXhap22ctmTEOSjqJ0p5ubb3LFvRyYJiIHApdazssDj0zD9PILyZtTgWuMRZexIvI7Vf1O6maJ9JYRqroJs+M8InIZ5lvMM0AHsL+IpPT5SoDfqxSuwzpVbbaXh2J62ekotvw1A3/CjBDuCdxrAhow1osG4AVgHzGTCTxMHn3fck5xLFSH6XFkmvAzUPVLOlJ1S5N9bhRmBPd3G0edfbag+tM2/qqbf+YYcH4ZymBKT49hvp3+2Z7vDvBbnSWNvehrKFfqA1PQHsco6HMZ7u+LMendhcmYZ2JMU88Cp2eQn4bJcBuBBcA3M8j8CNMgQMAMhTGlNeWR5nrM5I95No7DrP844BFMJnsGaNTeofy6EnK8H1OYmuzxmQwyVwP11j2bXhPEQVgTWB8cq4F7bYZ6DviS9d/V6qsJYyJ4wfL6MdZUVgSvCwKc7gbqSq07K3sD5nvdEnpNJ9+wXOdheuCTbd5sxxSuV226Uu/0XJuWvPmlpSHbRJGi9GZl78OYpV4FrrV+h1pujwA/Cch+D1NWWjDlxgvw+2w/dHipzQ9PYDpeZQNQ/r6A+ZaUyitXWf9yjMlzHvAQZkLJqZiy14L55n15IJyCdZgpjq1Yv9RaHa3AlMmLLL+nrQ4XYaxDKR2eiJnwVFAexUzS+p8c6RsofluUQes/HPiXlU3n970+48un0HzUD5s5/r4V47sM2Hcrc5wC3LAV47uZLN8wP8q6y5KGCKbB9z7qetva/D4mOvzQ+dl0/AOIftz5pT7GOTg4ODg4fORR9JR+MT+tzhORRSLymoh8N4PMyWJ+1n1FzE94exYbr4ODg4ODQzqKHqmJyFhgrKq+YD96LwQ+q70zlLAf319X1Y0icjQwW1X3LypiBwcHBweHNBQ9+1FVVwIrrbtVRF7H/JOwKCDzZOCRpzEfGB0cHBwcHEqKknxTE5GjgN9iZgrVAZPU/CuSuv9rzE90YBq8WlUtLzpiBwcHBweHAEphfgxhpnofh5n6Ogw4Jmh+DMjOwKy99oiqnpwlvDOAMwAqKyvrx48fj+fl/vTn+35uGb8LX8N4odwD0z7DsTJLlixZp6ojcwpmQTq/iRMn5h1vVpkS8/M8j8WLF5eMYy4dKpD0Ien7RMIeW6zU3JMwp8M+w0iT6a8OB4QfbLM6rKisrB83fiJJ30c8D1WTL1XVnuk5+74inhCsNnucfgIfj9RKYcGaVQMORQmJIAKePYSU2/iL+pSFSqfD0eMnEkv4IGLSEuCkqfQp+Koogt/jxrot/0ITkgGe5Tq8At5f2n8dZkMpGrUDMP9PKGZJmXIAVf1ZmtwemCmZm4Afq+rDfYW91z71+t0Lf8Hue+5NLOHTnfSJJQKHvV72zlt8+fBPsv3IKCEvQ7U4t5GmzpNonHlGzviamppobGzsU2bGjBkLVbWhr/T3hYaGBn3++ed74vV9ZV17Nyubu1i5qZMV9rymtZvu5jXMqN+Z7UdG2WFUlLqqSG9AJebX2NiIiJSE49771Ov/XPgLxk7dheXNnazc1MWK5k5zbOpibWt3j6wnMK6ukinDq5k8vIrtRlQzeXg1U4ZXMfHFz/J07MvbpA6/e9VdvLhoMcPHTqSlM05LV5xNnXFaOhM97vbuBGFRhtdUUltRRm1lmCGVZdZdxpAVt7Fad6N+3+nG394z7jCVZSFEZKvrsKGhQb/92zu59pHXKa+oRFF8W7P5qrZSVHyFZDxGXU0VVZEQVWVhKiMhqiKhnnPVyrtY7e/ELrvtZ/3DvffLQlRFzDPPPvsse9c3kEgqcd8356RPPGncCd+na/nrfO6ow0qiw/r6Bl1/+EXFBgOAoIQ8D08Ez4OQiHULIU/wBLq6Y8TUM41MFuwcW8+/f3VayXR4/MW3cf38d7KnWyAS8giLT111JdHyMNXlIarLw9RUhKmOhI175Z9Z7e/K7nt/ippy4xe1R3V5iGhFmKpImEeaFrDr3g1s6ozT3BHf7Jw69qvewEkzS6PDIIr6piZmEdtrgZ0x39VOwKz2vn9AZgbmB9MdMCs07E3vahWZwuzpYdSNnsBFT3XBU0/1mZabXp1PRQi2G+KxfV2I7es8pg4JMaRcoPMk2vxhNDU15Qyjra0tL5liEOQ3evRoTv7dg6xui/O/j81hY5eSSOtjlHlQGxE2dfs88O4rPf41ERhX7TG22mNs+HSGVlSyds6jDK80hagYfn3J9IUgxyGjJ3LxU13w1AsAREIwvEIYViHsPMTjoNFlDKsQhj3WxOL9prO6I8bytd288O46OgLLNgvfpy6ijHn934yu8hhd7TGqShhdZc6RkBTEsVT8Ro8ezUX3vUbCh8jb71BVJlSFMecyYWwEdqgWKkJh2rtixCRORzzGho3KB2uUjgS0x5Xu5EEA3L70pYxxhgSqyiAiSmj+nJ4RrWAqJLEXAowo9zH/KJeG37r332ZitU9ZuBsR2TzOgLs75uN7XXTHoKVTWZuE7qTSnYBYUulONhDzPf75/ht9J+LxBTlvX7ym//xMmjfneOauEfx4NzVVFZR5EAkJEQ/KQlDmCRF77uxopyZabUZVlnvPaKv1Ldp1GNEhI3LG3dbWTTRaScJXuhLQlbTnhNKVVDoTMP6hl4oqh+n8dvZWcv3K+9nwpS9S5kHIM3zCAmHPpl+EtrY2olEPMybzMctnBgOeTJsfIdq11Cy4hllhupvNl00Kxdv5YJFZTSuEMd8NA7N4YBlQC21tHQwEip0ocg5m2Zk9MD83voGh+rKInAmgqteKyFOY1TeqMCO687FLn6RDVa8HrgfYY6999McbH2fCD39AJORRXuaZc9gjkjpCHq9+90es/daPefH9jfxnWTP/freVpG9ahwlDK9k70kK0qpov7r0nu4ytpaIslJFMvr3gYhDk19DQoB3hKF6ohQOnjmHskErG1VUwdkglY4dUMK6ukqFVZYgIj86bx/Z77Mfba9t4e007S9a08fbaNl5c20ZTR50J/I1OKso8po4wo7ntR0bZflQ1O4yKMmV4NU8/sWALfr5vesPxpBJP+Cx44ok+30EhHHffax89Z+Pj7HXhOYwbUkmd5ZOOpgllfDst3uaOGO+u7+C99e0sffYanmndm67Kcby8voMNH8Q2kx07pIJJw6robO1iyNBKfFWSvhlVJFNuez56XEXJ+DU0NOiDP/k0Lzz7JIcfOiPnc7nyV/yhQ3mg7SR2O+TErKO9ls4477y/nNGjRhtzEEETmR01Kfht64rSYTq//z3h0wWNELNibiOPdpzE/kd+jY5Yks5Yko54go5Yko7uJB2xBJ3xJOuuvZZx3/4W4ZBHOCSUeR5lISEcsmfPQ8/7R7/5ZeI4+9TDi+c49yJjMWnMvZB8XvGMCZdUh5876lCaKjxOKFqHs0vHsci6NBuKbdSOw+y1A2b5lCbMitKo6rUpIVX9BvANEVkGPKyqR+YTeCTsMf7EI5k+LbfJVU/9PF+on8AX6s2kys5YkldXbOI/yzby4vvNPL94PCtX1/KX/3uSspCwy7gh7D2xjr0n1bHXxDomDavaoqLtiidZ09LNyk2drGrpYuWmLlZt6mKin8yUhH7jH9/6lM0Ae+eU80SYPNyY4w4NbEiiqmz49zHcteHL1Ox4KG+vbWPJmjZeWLaR+15e0WP798SM+MKPP0wsYRuxpE/C33xoOC7Rxn9lXVGwcJRbHe46bkjBz9ZVRdirKsJeE+tg3ZM0dU6hceYXAdjUEee9De2m0VvX3tP4tcYUrytByBNj+vGgzJqDUn5lXix3xAVieLScskxm7wJQ5vnUlPlsNyL3HohNTev7zCsDVVmUAp5AtTVbZUPTCYfRuPvYnOE0feUUuOWmUifPYRCgqG9qItIFvItZ+bkJ+KS9PklVX7MyEzETSCZiFni9UVWzfhhJGzbX33jjjUSj0ZzpMEPmHDKti/mgcxir/Dre2eTzdnOSpS0+Mds+1ZTBdnUh/ESC1mSIDV0+rRnqvcow/M+qJ/n2TZf12w6czu+OO+7oO/19cWxdTJu/pdmjO6msavdZ2a6sbPNZ0xajMlJGyNr6w54xPYQFQp65nrBgHjufckRR35wGSoeZOBYcjpWZOXPmR0KHBYWRJtNfHQ4IP3A6zCeMNBmnw34gj/W25mIWSk0/jsMsMhrGrKS8AWOEPc8+dzHwX8BYzBYmV2AWil2M3Q+pr6O+vl7nzZunfaFPmYen67x7r9vMK55I6mvLN+mfnn5Xf/S3F/WIXz2mB/2/OfrVPzyr5979sv5u7mL963PLdP7iNfrW6hZt6YypqurSr3xFgefzSX8+/PJKf18yGfgVHEaaTCk5DpQO+xWOlXE63Mr8VJ0O8wkjTcbpsPCjT/Ojqn462z0RWY7ZdG6OiPwIs2r7pfa54D5SK+0OqX/CbG+/2c/ZHwbCIY9dxtWyy7haTt5/MpCyA++b87l3Z82CP/5xK6TQwcHBwaFQFGt+/CWwXlUvtyuJvK2qx2aRfRr4DfALYDcN/JydJtczbAY+gZlUs8UeYmkYsZVlqrUE/8dg+L1ZorSV+h1MLiFHp8P80rZN6HAA+ZVaxukwCz4OOsyKvoZy5DY/Dsfs17QOu0mmfaYB8+0sFcYUzJ5ACzEbN+Y9lCSP4em2JlPoUYp4Bzu/jwPHwc7v48BxsPPb1jmqFml+BBCR2zBT9Q9T1Q77zPOYDeBSWI7ZYfdBVc04ld/BwcHBwaFYFPvz9VHA2cD0VIOWQUaAmzCr9P+qmPgcHBwcHBxyodj91K7GrA7ysIi8KCLXAojIOBGZY2U+hdmC/FAr86KIHFNAHNd/BGUKRSniHez8trZMoXA63PZkCoXT4bYnUzDcztcODg4ODoMGRe987eDg4ODgsK3ANWoODg4ODoMGrlFzcHBwcBg0cI2ag4ODg8OggWvUHBwcHBwGDVyj5uDg4OAwaOAaNQcHBweHQQPXqDk4ODg4DBq4Rs3BwcHBYdDANWoODg4ODoMGrlFzcHBwcBg02KqNmojsIiJN9nhKRNZnkGkUkT0KCPNAEXlFRLpEZEIWmRoRubPItD8oImtF5Py0tK4McKq3HN8TkYUl5Hi2iDwjIk+IyFV254N0mcki0u8FQkVkmog8aXk8ISJ7Wv/RIvKAiMwTkT+KyJ5WZqXlWAp+J4jI4yIyX0TuF5HaDDJF6VAMrrY6eU5ETgz4X5XiYtPxlIh0iEhDMbwCz10kIkuy3CtWbzU2vc0ickrAf4tykVb+NolIl/UfISJ/tu5CdTdLRJYGwh2fJY3F6G5vmyfni8ijIjLV+qfnm31tGl4WkTYRiYvI7SISFpGIiNwrIjMK1aGI/FhE5tqwDx0AfrWBsvesiBxm/YeKyEMi8pjlv4fV4Qci0iIirSLSbmVT/Lx+6PAHAf0tFZErS83RhtEZiOfr1u80y3m+iNwhIuXW/077Tp4RkVnWryef5sRAbNKW54Z1XwKuzeA/GzilgHCGAFGgCZiQReaHFLg5aYYwJgCzgPMDfo0ENkNNk/17CTnuGHD/DbN3XbrMVcA+RfAL07vA9aHAndb9G+AE6/4xcLp1HwDcXyJ+kYD7YuDbpdYhsBswz7prMLu0AxwF3GTdpwGXW73dBvy9GF72mdHAX4AlWe6XQm9j0tOWq1wAxwLPY7aDSvldB+zZD91tViayyBSruzFAjXUfA9yWK99Ybt8CrgVuAY62/udaXRTC72jgsgHm5wFh654KPGfdZwEXWncj8Ffr3tGevwysArYP8DuuP/k0kJY5wCdLzdGGsUUZsHxD1v0L4OtpHCuAJUBFMJ/miufDND+eAvwp6CEiwzCF5DzbmodE5HTbWj8jIl9LD0RVN6lqWx9xfR540MYxS0TuEZG7ReRVETk4n8Sq6gdZbh0pIgtsb78yIDsB+EqJOL4VuOwGEhnScYiqvmDjmG17qPeK2epnpzz4JdTmGqAWeNm6p2EqCYBngRnW/TRwcIn4xQKX1cBrGZJYrA5XADERKcM0ahus/3RM4wxwn70+BVN4JohIVX95WVwA/CxHukqht1UZ/DOWCxHxgG8DcUxlkcIczBZR/eF4mh0xXWLDT0dRulPVVaraai978n+OfPMapmPyJ0zjvtb6Pw7MLJDfl4AKEXlERG4TkSEDwM9X1VSZDpa91+01wFBgjZVP1QcnA6uBpL2eY/1mFcgRm/ZRwHaq+nSpOVqMETPqvFtEplgu76hqKv1B3aY4xiy/VN00B/hizliKaXmLaLGHA4uxI4O0e7OxvQxgJPASELHHS8DILGE2kX2k9mbAPQu4x7oPBO6y7m/aMNKPYWnPBkdqNfT2IC4FLgjc+wfwXok5Tgfmpodpw3g8LfzfWPdJwBXWfVEGfg8FnqsHnsLsVL6/9fs5cJZ1/xJ4IKDDTmBoKfgBX8fsjv5slvtF6RAQzP5/SzGVw3H2meuBRusWTL5cbN1/BfbuLy9gR+AWzd5LLYne0tOWq1wAXwW+B7xPwMoA7Az8sx8chwIhe9wKnDqA5a8akz93yZVvgM9hKsc3gb8FZCsxI5tC+D0Y0MNZwOUDwQ8Yj2l01wDHWr864AngVeADYEognuHASuCPafxeKlSHgee/A1w0UHUoMMKejwQeSQt/J+A5bH0a8D8vmCZsPs2UxtRR1M7XReDLGPOW9iE3FXhFbY9MRF4BtqO359VfLLTnZZjMgar+Hvh9IYFob+8R4HY275FPAh4uFUdrI78cmJlHmLA5x8Ntei/M9YCqLgQOEJH9MA3AfsBlwNUicjymQKyw4l+md7STC3nxU9WbgJtE5GzgfzE7qufLr08disgRmIpjB0zvfYGIPGA51FmxIRhz3l9UVWXLT5eF8poN/LQPHrl45aW3fCEiFZie/D2YyqYv9MlRVTcGwr8DU2Hd1ke4BZc/O8L+K/BzVV0UiD9TvrkO+D9V/Y6IXCsiX1TVTN+D8tHhBuAB634A+F0f3PrFT1WXAwfZEUwTxnpwNsYE/isROQC4BviMfeSnQBdwZh9pKaQOPRljpcgH/eG4zp4fFJFrUv5ivvn+EfOZoyvgfxqwB3BinmkC+NAatZOBb2S5F6M3XUuBPUQkYq93t36FYr2IVKtqu70ONgoCICLfxFTU6TheVTNW3iIyRFU32ctDMT3DFKZRIo4isgNwM/D5VMYIQlXXikhNuncwCBvORZjR3mZpUdUjRKQikKGagQ4b9iaMWQoRuQx42MqcDLQHK7Ui+KXHXZUhzKJ0aGU2qmpSRFoxvdYQ8BimZ38P5ntNhF6T6hQ212mheXMqcI1tHMeKyO9U9Tupm6XQWwauubAdpgG/BDNSGyoi31DVGzH59dVCOYpInao228v0MpBCsbprxujkHlW9JxB3tnxTg6kkwVTew6x7GmakVogOm4AGjIWkgc1NtqXi166q3fa6BWgNPJsq72tSPERkf8yI+9Oq2hkIq186tGFOA1Q3/9RRSo4xoNOWvz1SvERkBGb+wZmq+nYgPcdhrBX/pap+Bo7ZkWsYNxAHpqA/n+P+vpgh912YD6hnYkwLz2InKaTJT8NkuI3AAuCbGWR+hGkQIGBCxHz3asoz3TdgbPVL6B16fwvzvWk+cDdQp71D8o0l5Hi/jbfJHp/JIHM1UG/ds+k1PxyENYH1we8zlsc8e+xl/Q+1148APwno8HWsqawE/C4IcOt5j6XUIaYBuwVj4nkO+I719zA94AXAo8B/AuHeXQyvtGezTRQpSm9W9j7gHUxhvzZXubC6e560SU6Y0c1e/dDdpcAz9pmbgbIB0N0XgLZAHrkqW76x/N626Z0P3AtErfy5Ni2F8CvHmFXnAQ8BYwaAXz29ZW8BdiIYMA5T7prsO260/m9iTP8p7vUBfp8tVIf22YuB/8mRxmI57gf8x/J8DDvZA5P/PwhwSU0UacPk05T/+GA+zRlXPoXmo35gPrb+fSvGdxmw71bmOAW4YSvGdzNZvmE6HW67esuShhEYk+ug0F2WNEQwDZzn+H1kOeaVT1NTuB0cHBwcHD7yKHpKv4hMFPNj7iIReU1EvptB5mQxP0S+IuaHuj2LjdfBwcHBwSEdRY/URGQsMFZVX7AfvRcCn9XADCURORDzo+dGETkamK2q+xcVsYODg4ODQxqKnv2oqisx/0ugqq0i8jpm6nRw2u2TgUeexnxcdHBwcHBwKClK8k1NRI4CfouZKVQHTFLVlsD9X9O7EsV4oFZVy4uO2MHBwcHBIYBSmB9DmBUYjsNMfR0GHBM0PwZkZ2DWXntEVU/OEt4ZwBkAlZWV9ePHj8fztvz0F/ehK6l0J4CubiLV5VSGhbJMXwn9LnwN44VyD0x9388YV7rMkiVL1qnqyJyCWZDOb+LEiVvEm/Ah5ivxJMR8SPhKWKAiLFRk4lhifp7nsXjx4pJwLK+sqq8dOQ5fhaRCUtWc/c1/dEkh7EHEEyIhKPMgErJ8t1EdVlRW1teMmkDSV1SFpCq+QlLB1y05CoZX2BPCnuFb5kGYGJ6ECG0DOkzPoyPGTqA9rihiuRmOqWOL54GQgOdBSISQmOsQCUSEcChMyKPHv78cS6lDhozH7+OZvuABIubocQOeSMANqOJ5hrgE5Ai4BaWqrHQ6rB45gZZYfnW9pLsl6K+A9PhvIWsdmsr/Stb3OiHewjvr1/Rbh9lQikbtAMwyPopZUqYcQFV/lia3B2bpqE3Aj1X1YfpAQ0ODXnHFFezzyU/x4rJmXli2kReWNfPiso20dJml0uqqyqheuZzlQ0YBMCJazoHbD7fHCCYOq0QemUFT50k0zjwjZ3xNTU00Njb2KTNjxoyFqtqQUzAPNDQ06C9vn8MDT75EMjqKxatbeWt1G53xZI/M2CEVTBpWxevLN/ZkyqFVZdRPHsa+U4ay73bD2O2Nz/Nk9wkl49fY2IiIlITjlJ12V/ns5QyrjjAiWs7wqDmn3COj5YyoibD8179lwxe/ypurW3hjVSvvrmvvqTAjYY8dKlZSV1nD9IZ9+P/tnWtwlNUZx3/PJtncNndJwk2uioqKmuiIdTRexqID1Y76QaZUq4x1Om21re3YYh3asaPWjuM4TgcFL6it+KG11optpZLiDOJgrSAWEgJiIATMDchmk2yy+/TDeRNy2Wt2EzbL+c2c2X3f93kv//3v7nnP2T3PWVBZwDmVhVQUZjMy68ep8LD7pkcRDVBelE9JvpuSvCyK84Y+uvHkZFL33Et0LbmFg+0+Gtt9HOzw0dM3/CNfUZjNmaV5zCzNY2ZJ3uDzM0vzKC/IZsuWf0+oh9XV1fqdJzfw+/frOKMwj+LcLIry3BTnZlGclzVsuXXtWmT5t2nz+mnx9tLm9dPqPLZ19dIXGP1dIwKleSffG2WebNpbvqKisiLidS3KaeOum69PmodLH1lPS3MTZ8+bQ67bRW5WBtlZGeQ6JScrA3emi8bHHqfw/h/T3RfA5zelpy+Az9+Pr+EN9vnPpbjyrOHr/YHB+G5/gG5/PyFeimEs9LfyzlN3Js3Dl9/6Fxu3bKfqogvIdRtNee5Mo8+dMbjug2jvr001cX+XBoOKry+At6efzp4+Onv78fb0U/zYL1j08rqkeDiUhCo1MUleN2PycTVjBm8uxeQN/L4Tcw1mgN18TE7BOZhBfH8Oc8zBO4zSiulV5y9/gkZ3sXN/ANM9wvySDOYXu5hfnEFFnjB97Vp2fGslu9sD7G4LsLs9yLFeo6ssRziv0MucwiwunlFESU74O0Cv14vH44mo2ev1smzZsjEbMVRfRUVFVenKF/D1Q1G2MN0jTPe4Bss0j4v8LPOl3dnppcuVx96OAPUdQeo7Ahz1GY1uV5BZHuWcKTksKHExrziD3MzRt8Cx6vN4PAl96Q/XWFn13Nq1FBXEdt4B/AHlsDfIIW+QQ51BDh3z0tiVzYm+k7rys2CGx8WMAtfgY9DfTXZO7rDWkmkZ6uByeWYPK25Lnod/fP11fF1dMb+2A6gqx/1Ki09paT9MU08Bx4O5tPiCtHQrHT06rKWX6YKCLMUlLhRzF6wMtAjNQhCoyFUeuWLsHo7Ut2HDhrjeO6FQVXzHGmjuLSXgLua4XznRq5zwm9LpV473msf+QBCJ0lK7Y16QB1ckz8NkaKSzHm+wFE/RGVGPkZ+fP9hj0a+mdyagSn/QPJ/+xgb0nuUp5WG8Gsf7uzQcif5R5CFM2pkLMYP/9mDyke0UkfsAVHWNiHyIGR2fh/kcPozJADAKVX0ek2SWBedfpHMONnDrfcu55MwSFs0soiAna9Q+tSLcPuTuQlXZ19LFh/ta2bqvjQ/rfHzwlZtXGrqZOyV/sBV3+dwySvPdJ48T411+IgzVV11dreu/dyX9q37KopciT6lVW1vLshHX1tLZy8cH2tm+5Xneb7+Qd/b38baaLo7zphVSPauUS2eXcumcEsoLckbp6/YHaPX20uLtpbXTPO49WM/qpcPPk6jGogJPzK2LsDh3iIuuvYs9RzqpO3KCuqOd7DnSybYjnXT5BxK2C+YtGJ5fHvsoZi2hGKnv2muuiasVHJIQd8C9/QGaOro52NFtWnbtPnY1NDK1shKXDOnWEhlcdgl0tR2Oei3x6KupqUlcH8CmX1GbuZyaZbcmdhyS/zlMisZNq42HNbeN/RgDMR5XinqYRI0JehiORCu1mzHz7IBJnVKLSUWEqq4ZCFLVlcBKEWnEJPn9eiwHL8jJZHX+F8y+7qy4LkpEmF/uYX65hxWLZxN87ye82nonfTOuZuu+Nt78pInXtjUCcO7UQq6YV8biuWU0nQiwq+m4068tg/3jTg8yItDVl9zB6gunFVF75/Ix7TulIJsbL5jKjUc3c1X3VKpvuJv/Nnaw/UAH279oZ8P2Rl7eegCAWWV5lGb08uzurbR6e2n1+vH2jp7BpjzgY3UCesabknw3i+eVsXhe2eC6YFBpOtZN/dFODj79DNN++AMyM4QMl4ssl5DhksHlTJegD//l1AmIg+zMDOZO8TB3ysk73trao9TURB7mWVs7Kj2oxXLakGj3Yw9wANOlWAtc7iwvV9XPnZiZmD+QzMTkZVunqmE7ZEc2m9etW5f0JnN/UDlwIuh0VQbY2xGkL8ZfiX/Uto0H1j06Kbo9+oPKlyeC1HcE2dsRoNnbT1FOBoVuoShbKHIeC4c8P/uVF2m5d2USux/Hx8MxH4fkdyFPdNfVeHchnw5dV9bDMV5bCnkYlhjybW3CJEodWW7GJJ3MxGQ3b8d06a9y9vs18A1gKnAJ8DtM5vd6hsyHFKlUVVXp5s2bNRpRY967Wjf/9bmwm7v9/frR/jZ99d6f6993Neu7nw2Uw7px52F9Z+dh/duOw/r2jib9sqhCiZCQOZ5SVVUV2/VHi4miL6ZjjIhJpsaJ8DDm4zgx1sMJ1qdqPYzlGCNirIfxl6jdj6p6fbhtItKEmXBuo4g8iMna/htnv6HzSDU7s6O+hpneftjg7FNNTlYGl80pxXfHDdQsrIwYe+CWJbB+fcQYi8VisZwaEu1+fBJoU9XHnUwi+1R1aZjYbcDTwG+B83XI4OwRcYPNZmAB0MbJOYXCccYEx+RrEsbHYPTVJenakv0azEqiRuthbNeWEh6Oo75kx1gPw3A6eBiWaE05Inc/lmHm+2nFmcTO2aea4XM1zcbMmPwfzKSbMTcliaF5mmox8ZZknDfd9Z0OGtNd3+mgMd31pbpG1QS7HwFE5FXMX/WvU9WB2ZI/Zvisz03AZ8A/NMz4NIvFYrFYEiWhv/Q7OR9/Blw9UKGFiBHgBUyW/qcSOZ/FYrFYLJFIdD61Z4EC4D0R+VRE1gCIyDQR2ejEfA1YAVzrxHwqIjfFcY7Io5JTMyZeknHedNc30THxYj1MvZh4sR6mXkzc2JmvLRaLxZI2JDzztcVisVgsqYKt1CwWi8WSNkyKSk1EnhSRPSKyU0TeFJHiEDG3i8jnIhIUkeoh65eISJ2INIjIQ2GO/6KIfCUiu8Jsnykim0Xkf8457k+auJPniKgxnD5nW0SN0fQ5MeOq0XpoPUx1DxPR52yzHqbA5zDpYwTGowA3AJnO8yeAJ0LEnIsZZFgLVDvrMoB9mJyTbmAHIVJ0AVdhUnntCnP+qcAlzvMC4kj1lSyNofTFqjGavonQaD20Hqa6h2PVZz1MHQ9VdXK01FT1n6o6kFJ+GzAjRMxuVa0bsfoyoEFV96uqH9iAGTQ+ct8tmNyV4c7frKqfOM87gd2YVF9JI5rGMPogBo3R9Dkx46rRemg9THUPE9AH1sOBfU/553BSVGojuBt4N8bY6cDBIcuHSPAFFJHZwMVAYpNyRSbdNaa7Pkh/jVbfcNJd46TRl+h8aklDRDYBobIJr1LVt5yYBqAcmDOiT3cwZpyv0QP8CXhAw+SujLJ/RI3O9guAXIZrnBB9zjWOWaP10HqY6h5OBn3ONaS1xkQ/h5FImUpNo6fjugtoAS7UMNlLQtCEmcdtgBnOurgRkSyMCX/QMab6iqYRM4vBdxmSciwGUkaj9RCwHoYi3fVB+mtMGX1RSfRHuYkowBLMVDVTYoit5eSPm5nAfswkpgM/bi4Ms99swv+4KZiJTp8+1RoZ/QN1TBoj6ZsIjdZD62GqezhWfdbD1PFQVSdNpdaA6c/91ClrQsR8E9PP2wscxSRPBjOBaT3mnzurwhz/daAZ6HOOcc+I7VcCCuwccg03TaTGcPpi0RhN30RotB5aD1Pdw0T0WQ9Tw0NVtWmyLBaLxZI+TMZ/P1osFovFEhJbqVksFoslbbCVmsVisVjSBlupWSwWiyVtsJWaxWKxWNIGW6lZLBaLJW2wlZrFYrFY0ob/A8NHIFZvKwdwAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "measures = bpd.measures.Moments(measures='mean')\n", "grouping = bpd.grouping.Quantiles(n_quantiles=7)\n", "cluster_params = bpd.cluster_params(\n", " {\n", " 'measures': measures,\n", " 'grouping': grouping\n", " }\n", ")\n", "\n", "bdf = bpd.BipartiteDataFrame(\n", " i=df['i'], j=df['j'], y=df['y'], t=df['t']\n", ") \\\n", " .clean() \\\n", " .cluster(cluster_params) \\\n", " .to_extendedeventstudy(transition_col='j', periods_pre=2, periods_post=2)\n", "\n", "tw.diagnostics.plot_extendedeventstudy(bdf, periods_pre=2, periods_post=2)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.10.4 ('tw-env')", "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.10.4" }, "vscode": { "interpreter": { "hash": "cc42643cd328f586448453625167b9fc1f3293a5a6342212b4e72707ccae656b" } } }, "nbformat": 4, "nbformat_minor": 4 }