{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Task 0 — Datasaurus Warm‑Up (Starter)\n", "\n", "**Goal:** Show why we must *always* visualize by comparing groups with nearly identical summary statistics but very different shapes when plotted.\n", "\n", "This starter uses `datasaurus_task0.csv` (long format) with four groups: `dino`, `star`, `circle`, `bullseye`.\n", "\n", "**What to do:**\n", "1. Load the CSV (as strings first).\n", "2. Compute basic stats per group (mean, std, correlation).\n", "3. Generate a SweetViz report (optional but recommended).\n", "4. Plot x vs y for each group (facet grid).\n", "5. Write *2–4 sentences* reflecting on why summary stats were misleading.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 0) (Optional) Install packages in this environment\n", "Uncomment if needed." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "tags": [ "setup" ] }, "outputs": [], "source": [ "!pip install -q numpy pandas seaborn matplotlib sweetviz dtale\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 1) Load the dataset\n", "Load as strings first (safer), then coerce numeric columns." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "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", "
datasetxy
0dino61.0562093838706650.000786041836115
1dino57.9149519364229659.204465996007286
2dino61.4702319605998743.113610600617946
3dino57.8003536701023647.24321395851151
4dino53.5871245928257750.05299250969186
\n", "
" ], "text/plain": [ " dataset x y\n", "0 dino 61.05620938387066 50.000786041836115\n", "1 dino 57.91495193642296 59.204465996007286\n", "2 dino 61.47023196059987 43.113610600617946\n", "3 dino 57.80035367010236 47.24321395851151\n", "4 dino 53.58712459282577 50.05299250969186" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import pandas as pd\n", "csv_path = 'datasaurus_task0.csv'\n", "df_raw = pd.read_csv(csv_path, dtype=str)\n", "df_raw.head()\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Coerce numeric columns and quick sanity checks" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "RangeIndex: 568 entries, 0 to 567\n", "Data columns (total 3 columns):\n", " # Column Non-Null Count Dtype \n", "--- ------ -------------- ----- \n", " 0 dataset 568 non-null str \n", " 1 x 568 non-null float64\n", " 2 y 568 non-null float64\n", "dtypes: float64(2), str(1)\n", "memory usage: 13.4 KB\n" ] } ], "source": [ "df = df_raw.copy()\n", "for c in ['x','y']:\n", " df[c] = pd.to_numeric(df[c], errors='coerce')\n", "df.info()\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 2) Basic summary stats by group (fill in)\n", "Compute mean, std for x & y by `dataset`, and the correlation within each group." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "( x y\n", " dataset \n", " bullseye 54.502310 47.851315\n", " circle 53.569122 47.547658\n", " dino 53.833282 48.328789\n", " star 53.580462 47.571224,\n", " x y\n", " dataset \n", " bullseye 3.693820 4.779339\n", " circle 3.595391 5.104593\n", " dino 4.054927 4.994526\n", " star 4.175326 4.905460,\n", " dataset\n", " bullseye -0.018665\n", " circle 0.053956\n", " dino -0.049008\n", " star 0.037045\n", " dtype: float64)" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# TODO: groupby summaries\n", "g = df.groupby('dataset')\n", "means = g[['x','y']].mean()\n", "stds = g[['x','y']].std()\n", "corr = g.apply(lambda d: d[['x','y']].corr().iloc[0,1])\n", "means, stds, corr\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 3) (Optional) SweetViz profile\n", "Generate a quick report to observe that top-level stats look very similar." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "d:\\Projects\\43679_InteractiveVis\\VI_Lab_01_EDA\\.venv\\Lib\\site-packages\\tqdm\\auto.py:21: TqdmWarning: IProgress not found. Please update jupyter and ipywidgets. See https://ipywidgets.readthedocs.io/en/stable/user_install.html\n", " from .autonotebook import tqdm as notebook_tqdm\n", "Done! Use 'show' commands to display/save. |██████████| [100%] 00:00 -> (00:00 left)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Report task0_sweetviz_report.html was generated! NOTEBOOK/COLAB USERS: the web browser MAY not pop up, regardless, the report IS saved in your notebook/colab files.\n", "Wrote task0_sweetviz_report.html\n" ] } ], "source": [ "import sweetviz as sv\n", "report = sv.analyze(df)\n", "report.show_html('task0_sweetviz_report.html')\n", "print('Wrote task0_sweetviz_report.html')\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 4) Visualize — scatter by group (facet)\n", "Create a facet grid with one subplot per dataset and compare shapes." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxAAAAMQCAYAAAC3+YP9AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAizZJREFUeJzt3QmcVNWZ8P+nN5ZuFtkaFeLCpgKCkoBrjAFGowGMGiFBYkQdBDLxD44BI4jCBPUlLjiYUTsgGBIj8GpAo77OgGPUiAoGlc2GFkRBoFmUBlqg6a7/5xxSWNVdy62qu5x77+/7+WBbS3fd2s45z3mec25eJBKJCAAAAABYkG/lTgAAAABAAAEAAAAgI2QgAAAAAFhGAAEAAADAMgIIAAAAAJYRQAAAAACwjAACAAAAgGUEEAAAAAAsI4AAAADIEeflRZgQQAAee/755+WMM86QrVu3yqxZs/T/AwD8Y+PGjfLTn/7U68MAXFPo3kMBSOe6666T7373u7xQAOAj/+///T9ZtWqV14cBuIYAAjDIiSeeqP8BAACYigACcFFdXZ088cQTsnDhQvnyyy/loosukr59+x6/XZUwPfbYY1JeXq4v/+xnP5NTTjlFTj31VHnmmWdkz5490qNHD7nrrrukV69ex39v9erVMnPmTFmzZo3U1NRIv3795N///d+la9euvL8AYAPVvv72t7/VP1Vb3rt3bxk3bpy8+eabut1WVAnqv/3bv8kvf/lL2bt3r27TX3/9ddm1a5cUFxfr9v7Xv/61dOzY8Xgb3759ezly5Ii88cYbcu6558rcuXN5v2A81kAALlKdz+9+9zv58Y9/rDucE044QR566KGUv/Pqq6/KsmXLZPLkyfLwww/L7t27dedUW1urb3/nnXeO197ed9998pvf/Ea2b98uP/nJT+STTz5x5XkBQJAdOHBAbrnlFmnVqpUOCh555BH5+uuv5eabb5ZrrrlGt+nKggULdCmqWlB96623yt///ne54447ZM6cOTqwWL58udxzzz1xf/uVV16RkpISefzxx/VjAH5ABgJwSVVVlcyfP19GjhypOxJFrXeorKzUM1jJHD16VHc+zZo105cPHjwoEydOlPXr10vPnj11AKIyFGVlZVJQUKDvc/HFF8u//Mu/yH/+53/Ko48+6tIzBIBgqqio0FnjG264Qfr06aOv69Spkw4YVLsbLT0955xz9M+dO3dK06ZNdVv9ne98R1933nnnyWeffaZ/J1ZRUZFMnTpVGjVq5PrzArJFAAG45IMPPtDlRd///vfjrr/iiitSBhBdunQ5HjwoKt2tqNmv6upqXb6kApJo8KC0aNFCP87f/vY3R54LAISJKgdt3bq1jB49Wn7wgx/oyR9VgvqrX/0q4f1VO/2HP/xBZyLUDntbtmyRTZs2yT/+8Q9drhRLBSIED/AbAgjAJfv27dM/VQo8Vrt27VL+nprFipWff6zyUNXg7t+/X3dQbdu2bfB76jp1OwAgN6rE6E9/+pMuM1IlRyqL0KRJE7nqqqt0eWkiL7zwgi47VSWlqlz1rLPO0r+T6G8DfkMAAbgkGjiohdBqxinqq6++yvpvNm/eXPLy8vS6iPrUoj3VaQEAcqfabbWOTa0/++ijj2TJkiXy5z//WW90Ud/KlSt1+ZJaJK3WSUQzxzNmzJD333+ftwO+xyJqwCVqdw01+6T2C4/1v//7v1n/TbWrh1oHoWbEoouqFZV5UDt/fPvb387pmAEAx87zcP755+uJGVUuqtrze++9V5eLfvHFF8czw1HqnBAqS6w2vIgGD6qNfvvtt/X/q9sAPyMDAbhEpanHjh2rt1tVZUmqM1JrFHIJIBS1Xaua4Ro1apQMHz5cr7NQC6pVne0vfvEL244fAMJKLZxWg37Vpqq2VrXnauJGTdZcdtll+kzUyl//+le9vWt0m+1p06bJtddeq0tYVQnUxx9/rK9X69di17YBfkMGAnCR2tZPncNBzWaNGTNGn+9BpblzccEFF+h9ww8dOiS333673H333XrGS51rolu3brYdOwCEVWlpqcyePVuXjU6aNEm35WvXrtVbuqrJIBVEnH322XLnnXfqXfPUjktTpkzRmYh//dd/lQceeEBOPvnk4+eLoIwJfpcXUSswAQAAAMACMhAAAAAALCOAAAAAAGAZAQQAAAAAywggAAAAAFhGAAEAAADAMgIIAAAAAJaFLoAYMWKE/gcAMB9tNgCYJ3Rnot6+fbvXhwAAsIg2GwDME7oMBAAAAIDsEUAAAAAAsIwAAgAAAIBlBBAAAAAALCOAAAAAAGAZAQQAAAAAywggAAAAAFhGAAEAAADAMgIIAAAAAJYRQAAAAACwjAACAAAAgGUEEAAAAAAsI4AAAAAAYBkBBAAAAADLCCAAAAAAWEYAAQAAAMAyAggAAAAAlhFAAAAAALCMAAIAAACAZQQQAAAAACwjgAAAAABgGQEEAAAAAMsIIAAAAABYRgABAAAAwDICCAAAAACWEUAAAAAAsIwAAgAAAIC/AojFixfLlVdeKWeffbb88Ic/lFdeeeX4bVu3bpVbb71V+vTpIxdffLHMnDlTamtrPT1eAAAAIKw8DyCWLFkikyZNkuuvv15eeuklGTRokNx+++2yatUqqampkZtvvlnf79lnn5V7771X/vznP8vvfvc7rw8bAAAACKVCLx88EonIo48+KjfccIMOIJQxY8bIypUr5b333pNt27bJF198IQsXLpSWLVtKt27dZM+ePTJjxgwZPXq0NGrUyMvDBwAALlq7abcsWLpR9h88LM1LGsmwgd2kR6e2vAdAmAKIzZs36yBh8ODBcdfPmTNH/1QZhx49eujgIer888+XAwcOyPr166V3796uHzMAAPAmeJgx/33ZW3Xo+HVbtu+XCT/7NkEEEKYSJhVAKNXV1bpU6YILLpDrrrtOXnvtNX39jh075MQTT4z7ndLSUv1z+/btHhwxAADwgso8xAYPirq8cOlG3hAgTAGEyiQoEydO1GsfnnrqKbnoootk7Nixsnz5cjl06FCDMqXGjRvrn4cPH/bkmAEAgPtU2VIiVUmuBxDQEqaioiL9U2Ufrr76av3/Z511lqxbt07mzp0rTZo0kSNHjsT9TjRwKC4u9uCIAQCAF5qXHJtArK9FkusBBDQD0b59e/1TLY6O1aVLF719qypfqqysjLstejn6uwAAIPiGDewqrVs0ibtOXR46sKtnxwSElacBhFogXVJSIh9++GHc9Rs2bJBTTjlF+vbtq7MR0VIn5Z133tG/c+aZZ3pwxAAAwAtqtyW1YLrPGaXSpWNL/ZMF1EAIS5hUidItt9yiz+ugMgq9evXS54L4+9//LvPmzZNzzjlHnzhu3Lhxcscdd+isxMMPPyw33XQTW7gCABDCIGLqKLZtBUIdQChqwXTTpk3lkUcekZ07d0rnzp1l1qxZct555+nbZ8+eLVOnTpWhQ4fq7VyHDx+ufwcAACBbnFMCyF5eRJ3NLUQGDBigfy5btszrQwF8jw4YTqPNhlvnlFDrKSiJAnywBgKA/zvgVeWVUrF1n6wq36Uvq+sBwGScUwLIDQEEgKzQAQPwK84pAeSGAAJAVuiAAfgV55QAckMAASArdMAA/IpzSgA+34UJgH874C3bqxosQuSkTgD8ck6JhUs3StXBw/ps1qrtUtcDSI8AAkBW6IAB+BnnlACyRwABIGt0wAAAhA9rIAAAAABYRgbChzh5FwAAALxCABGAs2du2b6fs2cCAADAFQQQATp519RR7B4BAIBpqBxA0BBA+Awn7wIAwD+oHEAQsYjaZzh5FwAAwagcAPyKAMJnOHsmAAD+QeUAgogSppCevMuLekxqQAEAYUPlAIKIACKEJ+/yoh7T1BpQghoAgNOVA1u2V8X1f61bNNGTf4BfEUCEkBc7Odn1mHYO+E0NagAAwWFX5QBgEgKIEPKiHtOOx7R7wM+WuP5CtghAWCsHANOwiDqEvKjHtOMx7d7JgoVt/hENHleVV0rF1n2yqnyXvqyuBwAA7iKACCEvdnKy4zHtHvAnC2q27TooU8reZnBqELZBBADAHJQwhZAX9Zh2PKbdmZNEC9uUrw8f1TPcrIcwB9kiAADMQQARUl7UY+b6mHbvZBEb1Kz/dK8OHNxcWA7r2AYRAABzEEAg1JmTaFAz/pHXdW29kwvLWQScPbZBBADAHAQQ8BWnMidOz3CzZWxu2AYRAABzEEAALsxws2Vs7tgGEQAAMxBAAC7McLMIGADChbJVBBkBBODQDHds57Fz79eun3sDAOANylaT94fNSxrJsIHdOBO3zxFAAC51Hnl5IpGIe+feAAB4g7LVbxBMBRMnkgNc6jxU8NCiuJF06dhS+pxRqkumnDz3BgDAG5StfoMTgQYTGQjAxc6jtHVTeWT8pbzmAIxhQnmJCcdgJ85d8w2CqWAigAAcQOcBwA9MKC8x4RjsxrlrvkF/GEyUMAEOdR5qjUMs1jwAMI0J5SUmHINTO/upctWwl63SHwYTGQjAAZz4DIAfmFBeYsIxOIFz17izTTq8QQABOITOA4DpTCgvMeEY4Cz6w+ChhAkAgJAyobzEhGMAkBkyEAiM6C4elXsOSlX1EWlR0khKWxf7fjcPAAhyeYkJxwAgM3mRSOyprYJvwIAB+ueyZcu8PhQ4vItH7ExWWBevAX5Hmw07BG2bWMBrlDAhEBLt4hGU3TwAALlPMK0qr5SKrftkVfkufVldDyA7lDAhEJLt4hGU3TzgPGYogWBKtU3s1FFkIYBsEEAgEJLt4uGn3TwYwHr3WgXxRFYAgr1NLBDqAGLnzp1yySWXNLj+/vvvl2uuuUYmT54sixYtirutQ4cO8tprr7l4lPDjWT/9tJsHA1hvXytmKIHgYptYIIABxMcffyyNGzeWpUuXSl5e3vHrmzdvrn+Wl5fL6NGjZcSIEcdvKygo8ORYYa7YXTx27jko+6trpHlJkbRvXeKL3TwYwHr7WjFDCYRrgskPE0uAyTwPIDZs2CCnnXaalJaWNrhNbRBVUVEho0aNknbt2nlyfPAPP5+ohgGst68VM5RAcLFNLCWyCGAAoTIMnTt3TnjbZ599JtXV1dKpUyfXjwtwEwNYb18rZiiBYPPzBFOuKJFFILdxVRmIvXv3yvXXXy8XXnih/PSnP5U33njj+G3K/PnzpX///jJw4ECZNm2a7N+/3+OjBuzFmVi9fa2iM5R9ziiVLh1b6p8soAYQhOBh+twVScs+AV9mII4ePSqbNm2SLl26yJ133inNmjWTl156SZcszZ07VwcQ+fn5urzpiSee0BmJGTNmyMaNG+Xpp5/WtwFBQIrd+9cqzDOUAIKbedhffSTh7exCBd8GEIWFhfLuu+/qRdFNmhybUezZs6cOEObMmSNlZWUyfPhwadWqlb6tW7duei3E0KFDZfXq1dK7d28vDx+wVboBLNu8Wn+t3MJ7AsCPJ1j1y/bmMJfnayBKSkoaXNe1a1d56623dIYhGjzE3qbs2LGDAAKhQQ2reXhPEHQEyP5+X1KdYJVdqODrAEJlGoYNGyaPP/64nHfeecevX7NmjS5rmjBhglRWVsq8efOO36YyD4q6HQgLtnk1D+8JgsyLAJmAxd73JemGE8WNEt6f1x+Z8HQRgdp9Se2wpBZGr1y5Uj755BN9ArkPPvhAxowZI5dffrksX75cHnvsMb3+4W9/+5vcddddMmjQoKQ7NwFBxDav5uE9QVgDZCcHxqvKK6Vi6z5ZVb5LX1bXI7v3JdmGE3eN7JsweOD1h28CCFWipBZH9+rVS8aNGydXX321fPjhh3oBtVrvMGDAAJk5c6YsW7ZMBg8eLJMmTZLLLrtM7rvvPi8PG3Ad27yah/cEQeZ2gOx2wBKG9yWT3eV4/eG7NRBt27bVWYdkrrjiCv0PcIuJaVzOU2Ae3hMEmdsBMhk9Z94XqxtO8PrDdwEEYBJTF8YGaZtXEwO0sL8nQDYBsp3fZTJ63k5c8PojUwQQgE8WxpqydWkQA7QwvydANgGy3d9lMnreTlzw+iNTBBBADNK44Q3QAFgPkO3+LpPRy+y1sru95PVHpggggBikcZ1FgAYEgxPfZTJ63uL1RyYIIIAYpHGdRYAGBENQvstBWZMVlucFcxBAwHhuNoSkcZ1FgAYEQxC+y0FbkxX05wWzEEDAaF40hKRxnUOABgRDEL7LQV2TFdTnBbMQQMBoNITBzCBNHXWB14cFIOSTLUFdkxXU5wWzEEDA6BIlGkJ/I5UOwFR5eXkJr89Pcr2TOKcG/IYAAkYPMIOyUC+syCABMFXSMCHPXxMt9YOPft3b+359CsxHAAGjB5h+XKjH7hffIIMEwFR1kUji6+sSX2/iREuy4OO6AV1kxbpK365PgfkIIGD0ANNvC/Uo2YlHBgmAqUxpn3KZaEkWfKjgwc9rzZiIMx8BBDz98ltpwP20UM+ukp2gNJ5+zCABCAdT2qdcApkgZnmZiPMHAgh4+uU3pQG3ix2NeZAaT79lkACEhyntUy79oClZFDuxds4fCCDg6ZfflAbcLnY05kFrPP2UQXJaUDJLQFC+Fya0T7n0g0GbhAtqViWICCDg+ZffhAbcLnY05k42ngxgvROkzBJgF74XufWDQZuEC2pWJYgIIODqlz/oA1g7GnOnGk866sSviVufx6BllgA78L3IXZAm4YKaVQkiAgi49uXv2700FDOwuTbmTjWedNTeBlSk5QHvvxdBn8QKgiBmVYKIAAKuffkZwHrbeDKAjef255G0vP8w2HSem98LsrD++f4ELasSRAQQcET0yx9tQGYvWSPbdh1IeF8WRiV//ezEANbbgIq0vL8w2HSHm98LJrHcw/cn+Agg4GoDkggLo8zqqMMy6+p2QEVa3l8YbKZmVzvh5veCLKx7+P4EHwEEXG1A6vPbwig/D66tdNRhmjXyIiNAWt4/gjTYtLvdsrudcOt7QRbWPUH6/iAxAgg4pnLvwYTX5+eJFBbmS+OiQrluQBffDEyDMLhO11GHadaIjADCMNh0ot3yaztBGaF7gvL9QXIEEHCs09q5tzrhbXURkSM1dXKk5ogsWlYhp5/c0hcDcL92mpkI26wRGQEEfbDpRLvl13aCSQP3BOX7g+QIIOBYp3W0NpL2fn4agPu108wEs0ZAsAabTrRbfm4nmDRw73UOwvcHyRFAwNVOy88DcD93mlYxawQEa7DpRLtFO4GwfH+QHAEEXO20/DwAN6nTdGoxN7NGQLA40W7RTgDIi0Qi6etMAmTAgAH657Jly7w+lEBLtHCveXGR5OWpjENNXEcWXcznhx2O1DF6nZJN9tqe1KZE6iKRjF47P7zmCDfa7GC0W4BV9Ev+QAABVzstJVFHlmhQHBtc5HIMQRsgTylbLqvKK1Pex8pr59RrDtiJAAIID/ol/6CEKUdBHKA6Xf+Y6DondgoJwrar2a4vsfLaObmrFN8LAECmwrDbYVAQQOQgqAPUoOwU4seGyMrA2+r6knSvnVO7SvG9AOAkJiiCKwy7HQYFAUQO/DhANbVxdmKnEL81RFYH3okWRWbz2jm1q1RQvxfwHwaaweOnCYr6n79+3dvLe+sqqVgI+W6HQUEAkQO/DVBNbpyd2CnEbw2R1YF3/R1Q8vPyZMfegw0Wp6d77ZzaVSqI3wv4j58GmvDXBIWVwDTR5++DDbskdtsaNz+PfgmmTdrtEKkRQOTAbwNUkxtnJ7YF9FtDlMnAu/76kmx2WXFqK8Ygfi/gPyYMNGE/rycorAamiT5/9fe8dOvz6Kdgmi2C/YMAIgd+G6Ca3jjbfdIZvzVEuQy8s33tnDjRTxC/F/AfrweacIbXExRWA1OrJ1N14/PodTCdafaDE9D5AwFEiAaofmic7W6c/NQQBWXgHcTvBfzHj20ZzG8nrQamVje7cOPz6GUw7afsBzJDAJEjPw1Q/dA4h7lxCtLAO2jfC/iP39oy+KOdtBqYJvr8qROpxpYxufV59DKY9jr7AecQQMCoxjnsjRMDbyCcbRnMbiejme7KvQelsCBPjtZGUgYCiT5/fbuXyop1la5/HjMNpu1ccE0pYXARQMDXg1i7Gie/7FABIJhtGcyVKNNdWJAv7Vs1lfZtSpIGAok+f4Mu7uzKMWcbTNud1aeUMLg8DyB27twpl1xySYPr77//frnmmmtk/fr1Mn36dFmzZo20bt1abrzxRrnhhhs8OdYg8+sA2o7GKUhlUACCxa9tc5AkynQfra3TwcPUURdIkIJpu7P6lBIGl+cBxMcffyyNGzeWpUuXSp4qEPyn5s2by5dffikjR46U/v37y9SpU+WDDz7QP0tKSuTaa6/19LiDxMoA2tROzI7GyekyKFNfOwBmY3LDDGEqw7H7uVJKGFyeBxAbNmyQ0047TUpLSxvc9vTTT0tRUZFMmzZNCgsLpXPnzrJlyxYpKysjgLBRugG0yZ1YNo1T/QF95Z5qWxvM2L+vYuKde6vjTvJmymtHYAOYLWhrvPwqTGU4TjxXSgmDyfMAory8XAcGiaxcuVL69eung4eo888/X5588knZvXu3tG1LA+rGjIMJnZhdW7Umq2W1q8FM9PfrS/TauT2YNzkodEImry+BFez+TNndNq/ZtFumlL1NNtMlYSrDCdNzRQAyEK1atZLrr79eNm/eLKeeeqqMGTNGr4vYsWOHdOvWLe7+0UzF9u3bCSBcmnHwOn1r52A3WS2rlV01sv37iezcczDl81tdsUduHtLdsQV3JgSFbsnk8xO2wArZcetzkqxtPlJTJ6vKd/HZdCmgC1MZTpieK3wcQBw9elQ2bdokXbp0kTvvvFOaNWsmL730kowaNUrmzp0rhw4dkkaNGsX9jlovoRw+HLzaQ1NnHLxO39o52E0WDLVvXSztW5fk3GBaPfvozi+/1h2eeoxkQc2cF9bK6Se3dKTh9joodFMmn58wBVbInlufE9U2V3z+peyv/qYE0unH9CM3Ajovt45NFxTZHTxRcgTjAwhVmvTuu+9KQUGBNGnSRF/Xs2dP2bhxo8yZM0dfd+TIkbjfiQYOxcXFnhxzEKWbcfA6pWnnYDdZMKSCBzt207B69lEVIEQ7/mTPT2VEMhkcZNKJeB0UuimTz0+YAitkz83PSczeIq49pt8EMfC3GhRlcj8280CgSpjUjkr1de3aVd566y058cQTpbKyMu626OX27du7doxhkGrGweuUpp2DXaeDoUR/Pz9PpC7m7KP1O/5UQYfVwUGmM3DZvg7qcZ56Ya18XnlARCLSsbSZ3Dykp6XPglcdWCafnzAFVsieW58T9X2J3YDBjcfMhVff8SAG/laDIiv3ozQTTki8etQlKtPQp08fnYWIpc75oMqa+vbtK++//77U1tYev+2dd96R008/Xdq0aePBEYc9wLhAHhl/qf7pZj2kGuyqwW2sbAf90WCozxml0qVjS/3T7jR3/b/fpeMJCe+7c+/XMv6R1+VA9WEpUFFGDoODVJ2IXa+D6oTum/eebPj8K/n68FH5+nCtbPx8n0yf+56+LZVoB7aqvFIqtu7T9dvqcrrfc/vzY+dnDcHl1uckXUmkSZ9NL7/jQQz8rQZFVu6Xaf8AGJ+BULsvderUSW/Tqs7voBZTL1y4UJ/v4bnnntNBwuzZs2XSpElyyy23yEcffSTz5s3T90V42J0Bcbq+s/7fTzT7o8oS9lcf0f+U4sYFcrimTmpjUhWZDA6ymYHL9HVINhuq6rPTlQp4WWKQyefH62wb/MGtz0mygXHTxoVy1mmtjfpsOvUdt5LV8LrMNhfJnp/VoCjZ/fLz8mRK2XL9d7ftUhnjYGVoEPIAIj8/X5544gl56KGHZNy4cVJVVSXdu3fXC6ijuy+pAEKdifrqq6+Wdu3ayYQJE/T/w1xOpLH9vKir/mBDZR6igUNU9eFa6XbKCdKsaaOsBiRuzMClmg1N1xF5XWKQyefHz581BKtEJ9nA2MRdwZz4jlstvfFr4J/q+VkNihLdr3lxkezYe1CqPv8q5eP7OUMD73m+BkKdy+H+++9PenuvXr1kwYIFrh4TsketZfpBqSpbqh9AKAera6Sk6bFdxyKSYNFECm7MwKVaq5GuIwpiiQHCxYu2zU8DYye+45lkNfwY+Kd+fhdYeu8TfUYOVB/Rpaap+CVDA3N5HkAgWIK4G4ZbHa3a2nXb7oNZDU7cGGioIOWTrV82KGNSs13pOiI/lxgAXrZtfhkYO/Ed9ypz6VamKd3zs/re17+fmqRKVvrWoV2J0YEo/IMAArbyulTFDxJ1tMdOZFeX0+DEjbUdd93YT556ca18vvNYTe23SpvJTUN6pO2I/DSTCiRC2yaOfMdTDdbdzlyqY1Hn39m0bV/cejSnMk1OPb9kf1etm7Fju3JAIYAICFP2eKZUJbuOdufeg7Jt1zfZB1MDL3XsD/1/3wv0TCqQCG2b/d/xRGVhn2z9Sp+XJxKJSF5ens5wxp5Ir35Ww66+L9GxOJ1pciozS8YXbiCACACT1h0EreFyKjCr39Gq3TISBRCsEQDMELS2zdSyMFUiWXXwm/r9FiVF0u1bJ0hdJNIgq5EuAMmkzU50LLHWf7pHP56dfapTmVm/ZXxNmQBFZgggAsCkdQd+a7hMCcwYnABmC1LbZop057mIBhRdOjZKWHpjJQCx2manOxZ13hvVH+TS/icbKDvRT/sl42vSBCgyQwARAKbV5vql4TIpMPN6cBLt2CrV1n8Ha6RFcZGUtilmJgio9z21+t1nVjW3nd2s9GVWAhCrbbaVY8ml/WegbP4EKDJDABEA1OY6I1nn5EQq28vAK1HHpraZVTtCMRME2POd4rtkLfOaSLJSzlwDkGyOpf7fshooMlD2xwQorCOACADKX5yRrHOyI5Xthlw6tnQzQcyuAskxWMwu85qfnyfbdx9MuWjazgCkflvWpmVjaXtCE6mriyQ84Wf9v5VJoOjGQNmP7TIToP5FABEATpW/+LExslOqzsn0FKsdHVuqGTdmV4HkmFXNPvOq2herfVm075s+d0XCwX66ACRRWxY907eS6LbYv5VJoOj0QFk9l/vmvRd3nh61oFxtvW1Hv+3UeIAJUP8igAgIu8tf3BgkZtIgeRHMRDune8qWy+Ga+HM0KGrrVVPZ0bEl6+CYXQVSY1bVvb5M3b9966YJAwh14rRUfVauZ4LOJFB0eqD81AtrG5zkU11W5+3JduttN8YDXq//Q/YIIJCQ04PETBqkbBqvTAOOZPdX/xoVFcrhmoad0/56jbVJcu3YUnVwzK4iLLKduGBW1V2pTpyW6v3K9UzQmQSKTg+UP688kPj6f5700+TxQFA2XgkbAggk5PQgMZMGKdPGK9OAI9391Y5EiWa31AmOTJVtx6ayKiowUs+tfZuShB0cs6sIg1xmXZlVdVe2AVuubVmmj2vXQDlRYCvyzZmzY9XVNcyeZ4pJIyRCAIGEnB4kZtIgZdp4ZRpwpLt/aZsSvSNRfWqAbSonOzZmVxEGuc66MqvqnmwDtlzbMi8CxWSBbZsWTWRrgpOR1hyty3nXQCaNkAgBBDwZJGbSIGXaeGUacKS7vx8HzE52bMyuIgyYdfWXbAI2O9qyTB831/V8yQLbtt86QfLzROrqJSLU5VxLjfzYB8J5BBDwZJCYSYOUaeOVLOBQ2/KNf+T1Bo12ugDFrwNmJ2dAmV1NLuy7lwUFs67h4GRZkR3r+awGtnWRiJzUtkS2JchC5Fp67Kc+kPbXPQQQ8GwAarVByrTxShRw5OUdOzladC1DbKNtJUBhwAwr2OI2OJh1tYYBm/XvvR2LkVMFtuq2RAGEHaXHfugDTW1/1wZ0UokAAp7JpEHK9L6xAUeiEwLFNtp+mF0JagMUNGxxGxx+aBe8ZuqAzdTvvR1lcekC2zCXGpnY/q4N8HeEAAKBFBtwqLKlRLsoxTbaJs+uuLGNLexB3XywmNwu+HnA5mX7lO6xszk2q997O8ri0gW2cWf2Vql3EZm9ZE0o+gET298FBgY1ngUQhw8flsaN7dmJB3CD32uZnd7GFvbx+2cNcHrA5mX7lO6xsz02q997u8riUgW20dvC2A+Y2P7uNzCosUt+pr9w0UUXyT333CMfffSRM0cE5Eg1nFPKluvMw5Syt6Vf91LdSMfyU1rXzm1s7fLXtz6R4Xe/Itfe+aIMv/tlfRnHBgh+/qwBTg/Y3Gifsn3sbI/N6vc+mj3oc0apdOnYUv90akDv5evsFRPb3+YGBjWeZSBuuukmWbJkiSxcuFBOP/10ueaaa+Sqq66Sdu3aOXOEQAaSzbpcN6CLrFhX6ctaZru2sVUnibODChbKFq+RyD+3CzxSU6cvK4Mu7ixhRt08wiSbGXUvZ2TTPXa2x5bppiC5lq5YKbMK8sy3n9rfYQHeAjfjAGLs2LH63z/+8Q/5y1/+Ik8++aTMnDlTLrzwQrn22mulf//+UlRk7hl6EWzJZl1U8DB11AVigkxrbO3bxrY65xMKKc+8uuF48BClLv/51Q0pA4iwrMugbh5hkc2AzYsZ2Wjbs23XgZSPncuxufW9t1qaFOSZbz+1vz0MDGo8X0Tdp08f/e/uu++WN998U+bNmyfjxo2TFi1a6KzEiBEjpEOHDvYeLRwVhAGe6bMu2dSlZrON7eqK3XK0ti7u+qO1EZk6+10587RWOb23h2uOJrz+UJLrFbvqcYPwGQXCPGBze0Y2UdsTK/ax/TBbbHVNnB+eS1j0MCyoMWIXpu3bt+typldeeUXKy8t1SdOll14qb7zxhjzzzDNy//33y5VXXmnf0cIxQVlwZfqsS7Y7MmS6jW37Vk1l2+6GJUtfHz4qq8p35fTeNi4qlCM1DXe1alJU6OhOFEH5jAJh5uaMrGozps9dkXAXvqaNC+Ws01rHPbZJs8XJJkusTpJl81yYoIGjAcSBAwfk1VdflcWLF8v7778vTZo0kR/84Ad6YbXKSCgTJ06UW2+9Ve677z4CCJ8IylZjyWZd+nYv1QurM5m5dqIxdStDUtqmJGEAYcd7O/zybnFrIBS1W+BPL+/m6PMOymcUCDs3ZmSjEw6JggelQ7uShGWtJswWp5ossTJJVr/vuuWqnpb6OyZo4GgAoXZhOnLkiPTu3VumTZumA4Ti4uIG9zv77LNl3bp1mf55eMT00h+rEs26qOBh0bKKjM+j4ERj6laGJFEgZdd7G13noNY8qLIllXlQwUOq9Q92PO+gfEYBOC/RhIOJWelMJ0vSlSal2kjkvXWVSSfEmKCB4wHE9ddfLz/+8Y+lU6dOKe83cuRIGTNmTMYHBG+YXvqTifozSCrzkOnMtVONqVt1qbGB1PpP9+rSJTvfWxUsZLLjkh3PO0ifUQC5S5UlTjbhYKXt8bqUJ9VkSbrSpGR915wX1sWti6s/IcYEDRwPICZMmGDpfiUlJRkfDLwT5AVX2TSMlUm2PE03252u43GzxjbVCYXcfm/teN5B/owCJvN6QJ1NljjphENxI0+yz5lIN1mSqswqWX9Xf1ON+hNiTNDA1UXUCA6TFo/ZLdOGUXUgasvTTH4nk44nUePvZAdtynuba22xKc8DCBMTBtTZZImTTTikO24TSnlymSxJ1t+lmxBjggaZIoCAUYvHnBhUZ9owqsdUW57WV1iQn7IBz7bjcaODNu29DfvzAPwiXbvmVXYiXWY52wkHE0p5cpksSdTfFRbkJezTYifEmKBBpgggYKRcB9XRTk2VIlUdrJFGRXnSvLhIWpQ0kvatS1I2xsk6ELU1aqrHzrbjMWHGCwAybde8zE5YySxnM+FgSilPtpMlVjcSSTSJxgQNMkEAASPlMqhOdeKgosKCrM+U2r5NiSMdjwkzXgCQabvm5eSHUyU3QSjlSRQInH5yS8o/YSsCCBgp2aB6Z5LFzVa377PSuWXbgWT7e8k6aLUOQwVD1PgD8Eqqdm32kjWeTX7kUnKTquwqqKU8ZBdgNwIIGCmXQXWq7fv039iTOgjJtgPJ9veSnbNhf3WNzqR4vVgRQHilatdSZSfcWBuRzaDYStkVg20gPQIIuM5Kx6IG1asrdjfYek4tBEuXQUi3C4UamKc7jlzqTzP9vWgHPX3uigZnTWUtBACvJWvXkmUnVM29F2sjrPQtrDkD7EEAAVcb70y2OlWLlrftPphxevxY8LEr4a4TSvOSIuO2JtTPt3XTBgGEF2shTNzzHYB/shNeDNKttul+XnPmVNtMm49sEEDA1cY7k46ltE1JwgAi3aJk9Xg3D+khZX9ZI4lCCLULk4mzUCbs/mFaYAWEiR8HcomyE/sP5rY2IpvXwWqbbkI7a1LbTJuPbOVn/ZsINdXoTClbLuMfeV2mlL0tc15Ym7Txjr2/+mm1Y1GZBJUOj2V1N4xBF3eWUVf31PtfJ/p9E2ehcnm+dknVCQNwTnQgt6q8Uiq27pNV5bv05WRtpslyGaRn+zpYbdNNaGdNaptp85EtMhCwZcaiID9+oJ5ur3ArHUuuu2GoICLZ1nUmzkKZsPuHiYEVEAYmZkW92Ap11sIPsnodrLbpJrSzJrXNtPnIFgEEbOnoausiGe0VHitVx5LrbhiZLv7zehbK690/TAysgCBIV5YTpIFctoP0v771iWzblXiXPCtr36y26V63sya1zbT5CEQAsXnzZrnmmmvk7rvv1j+VyZMny6JFi+Lu16FDB3nttdc8Okok6+hUFiI2kEi3V3ijonzp2amtJ7M/fp2FcpqpgRXgZ1bqzIM2kMtmkP7MqxuS3mZl7VuQ23ROnAfTGBNA1NTUyB133CHV1dVx15eXl8vo0aNlxIgRx68rKCjw4AgRlayj69yxpTRr2shyuZAKHqaOusCzRYZ+nIVy+rUKeicMmFqe5HXwbsIC7sM1R5PeZuV1iLbp0eeiJq/8shg9HafaZtp8+D6AmDVrljRr1izuukgkIhUVFTJq1Chp166dZ8eGeMk6upsG90jYmNnVMbJbhDuvVdADK8BtVsqTvBzImdK2Ni4qlCM1DbeyLm5SaPk4THkuTnCqbabNh293YVqxYoUsWLBAHnjggbjrP/vsM52R6NSpk2fHhoaiHV2fM0qlS8eW+meqxjnT+yfDbhG8VoAfZbLAV2VlHxl/qf7p1oDXlLZ1+OXdJK/efhzq8s+uONN3zwUIOs8zEFVVVTJhwgS91uGkk06Ku23DhmP1kPPnz5c33nhD8vPz5ZJLLpHx48dL8+bNPTpiZDNjYccMR5AWGTqN1wowh9flSX5pL9TOecqfX90gh2qOSpOiQvnp5d2OX++n5xIEJpS1wVyeBxD33nuvnHvuuTJ48OAGt6kAQgUNpaWl8sQTT+iMxIwZM2Tjxo3y9NNP69sQHkFbZOiX14pOBJBA15mb1LaqYCGTgMHk5+JnQS4FQwACiMWLF8vKlSvlxRdfTHj7mDFjZPjw4dKqVSt9uVu3bnotxNChQ2X16tXSu3dvl48YXjJ9Fs+kATfrTgCzmFxn7pe21annYkKbbZognZcEAQwgnnvuOdmzZ49ceumlcdffc8898vLLL8vs2bOPBw9RXbseawR27NhBABHyWbz8fxbLmrTThhezNsk6PztmPOlEgOAzPUPi5HNhpj33UjACsHDyNIB48MEH5dCh+Aj3sssuk9tuu02GDBmi10ZUVlbKvHnzjt+uMg9Kly5dXD9eeC92mz4T06tuD7jTvQ6sOwHg9wxJLs/lm8Ft4okmJklyKwUztS9GwAOI9u3bJ7y+TZs2+rbLL79cxo4dK4899pgOKNSJ5qZNmyaDBg2Szp2zr5GE/5na6CebtVn/6R4Z/8jrtmdKnH4dqCcG4NdZZiuDWxZd51YKZmpfjBAsok5lwIABMnPmTCkrK5Pf//73eucltdh63LhxXh8aPOZUo59rJ5lswP314Vqp2LrP9tkZpzu/RJ1IYUG+9O1easvfB2A+v84yWxncMkmSWykYAVh4GRdAqDNPx7riiiv0PyCWE42+HZ1kogF3fX7KEKjnfd2ALjLnhbVytDairztaWyeLllXI6Se3NHrwAMCeyRK/zjJbGdwGaQG5F2VtBGDhxT6o8CXV6KtGPlaujb4dJyCqf9K8po0LHc8Q2P061PfeusrjwUMUJ2YC/Cs6WbKqvFJnRleV79KX1fWJWJ1lVr8/pWy5LtecUvZ20r/nluhGGw2uz8+z/USnYeVGHwQzGZeBALzaNcSuVGzsrI3qTFUn7WSGwOndU0hRA8GSaUbByixzthlcJ9dWxE97fKP66xrdNsc+pjrzN8K9gxcyQwCBUOy0YYUTqVg30uNO755CihoIlkwnBay0Y9mUOTm9tiISSRxCbN9TLVt3HXTkMcMoSDt4wToCCPieXZ2QE4N9u2dnvNgJhRphIFgynRSw0o5lk6n0ahe52rrEJZkMggHrCCDge3Z1Qk6lYu2anfFqJxRS1ECwZDMpkK4dyyZT6cUucgX5eQ0CCDsfEwgLAgj4np2dkMmpWC93QjH5dQHCxI4spBOTAtkEJW7sIlf/eR6oPiIbPv/KsceMvj+Vew9K1cEaaVFcJKVtin1x3owwnBsE9iGAgO+FpUafxcxAuNmZhbR7UiCboMTtNWLq9XvqhbUNshB2PWai92d/9RHZtvtg4NZZ+PXcILAPAQR8Lyw1+mEJlACIL8/HkGlQ4mZ5ZKIBrwokOndoKTcN6WHLYyZ6f0x8n8LwWYTzCCDge2Gp0Q9LoAQgPFlIt8ojEw14VRaiWXEj2/qKZO9PEN6nMHwWkRkCCARCGGr0wxIoAUiMLKTZA95k70+qbLFf1xHwWQQBBOAjYQiU3OTXzhvhFPYsZC7fVzcGvInen6jCgjzp2700MOsIwv5ZBBkIACHl584b4RTmLGSu31e3FmxH35/PdlbJ7q++eayjtRFZtKxCTj+55fHj9fM6gjB/FnEMGQgAoeTnzhvhzVyFNQuZ6/fVrQFv9P2ZUrY8LoBIdLx+X0cQ1s8ijiGAABBKfu+84Q0yV/79vro54LVyvKwjgJ8RQAAIJTpvZIPMlTf89n21crysIzAP6+KsI4AAEMrGmM4b2SBz5Q2/fV+tHC/rCMxCdjEzBBCAjxu7OS+sla2V+/XljqXN5WabTogUhsaYzhthmAkPCr99X60eL+sIzEF2MTMEEIAPqcH59Lnvyf7qmuPXbfz8K7lv3nty1439jOhU/dAY03kj6DPhQctgTh11gS1/x402kvbFX8guZoYAAvAh3RnGBA9RVQdrjBmg0xgjiPw2E+5XdmUwTc+EwhxkFzNDAAE4zInZr2SDc5N2EaIxRlAxs2xeBjNZO+uHTCjMQHYxMwQQCC030tpOzX4lG5ybVItNYwzAjQxmqnY2DJlQkzer8BOyi5khgEAouZXWdmr2Sw3OKz7/skEZU4uSImNqsWmMAbiRwUzVzgY9E/rXtz6ROS+sk6O1dcevo0Qre2QXrSOAQCi5ldZ2avZLNXKTRvaTp15YK59XHtDXfat9M7lpsDm7MCk0xgCymUXv17295cXqqdrZW67q6dqid7czAdGd+I7WRuKup0QLbiCAQCi5ldZ2cvZLdUwPjfuemI70OoBsMsLXDegiK9ZVpl2snqqddSsT6sVibRWs1A8egliiBTMRQCCU3Eprh30dADugAMg2I6yCByvbtqZrZ93IhHqxWDvVZhpBKdGCufK9PgDAC6rDUR1MLCcG9tHZrz5nlEqXji31zzBtH5iqUwUAOzLCJrSzXizWTjYRVliQH5pJKniHDARCyc0FvmFeBxCGHVAAeJ8R9rqd9WKxdqLMS2FBntw8pHtoJqngHQIIhJbXHU4YBH0HFAC5C0KppxfPgZ3u4CUCCACOCcLAAICzgjAQ9uo5MBEGr+RFIpHES/gDasCAAfrnsmXLvD4UIDQLqU0eGLBLlNlos53DZx9AtshAACnQwQZ7hoxdohBWfPYB5IIAAoFh92CfDjb4vNh6ETABn30AuSCAQCAG/04M9ulgg49dohBWfv/skx0GvEUAAc/ZMfh3YrDv9w4W6bFLFMLK689+LgEA2WHAe5xIDoE42ZgTg32vO1j494SCaoAzpWy5jH/kdZlS9ra+DITxZJqpAoBV5ZVSsXWfrCrfpS9b/Z5wgkrAe2Qg4LlsB/+xM1g791bbPthnC9Lgc2LrRWZHEcbPfiYZhVwzxmSHAe8RQMBz2cz0Jxqk5eWJxG5KnOtsWhD2Jof7u0SxdgZh++xnGjTnGgCQHQa8RwABz2Uz059okKaChxbFjaS0dVPbBvsmb0EKMzE7irDJNGjONQAgOwx4jwACnstmpj/ZIE0FD4+Mv1TCip1JvMfsKMIm06A51wCA7DDgPaMCiM2bN8s111wjd999t/6prF+/XqZPny5r1qyR1q1by4033ig33HCD14cKm8XO9H8zCF6TtJaWQVpD1N6bgdlRhE2m7bEdAQDZYcBbxgQQNTU1cscdd0h19TeLYb/88ksZOXKk9O/fX6ZOnSoffPCB/llSUiLXXnutp8cLbwfBDNIaovbeDMyOImwZx2zaYwIAf77XgHEBxKxZs6RZs2Zx1y1cuFCKiopk2rRpUlhYKJ07d5YtW7ZIWVkZAUTIB8EM0hqi9t4cDI4Qpowj7XF43mvAqABixYoVsmDBAlm8eLFceuk39esrV66Ufv366eAh6vzzz5cnn3xSdu/eLW3b8sUJ8yCYQVp2ZQTMZAH+Z1rGkfY4PO81YEQAUVVVJRMmTJDJkyfLSSedFHfbjh07pFu3bnHXlZaW6p/bt28ngAgg1jY4W0aggofpc9+T/dU1x6+r+PwrmTSyHzNZgI+QcUwsiBMkvNcwkednor733nvl3HPPlcGDBze47dChQ9KoUaO46xo3Pjabevhw9mcYhrm8PDuq30XLCPqcUSpdOrbUP+unuOe8sDYueFDU5adeWOvBEQPIFpMt9p/h2lS81zCRpxkIVbKkypRefPHFhLc3adJEjhw5EnddNHAoLi525RjhLmppc3/9UqW0t1buT3j955UHcnxkAG5iIwn3Sn28zmrwXsNEngYQzz33nOzZsydu3YNyzz33yMsvvywnnniiVFZWxt0Wvdy+fXtXjxXhrqX1ugOxT57XBwDABky2uFPqY8ICZt5rmMjTAOLBBx/UZUqxLrvsMrnttttkyJAhsmTJEnn22WeltrZWCgoK9O3vvPOOnH766dKmTRuPjhphY0IHYpdvlTaTDZ9/1fD69vE7oAEwn4mTLUEr9TFlATPvNUzj6RoIlUU49dRT4/4pKjhQt6lzPRw4cEAmTZokFRUV8vzzz8u8efPk1ltv9fKwYdDAfkrZchn/yOsypextx+pcU3UgfnPTkB7SoqQo7jp1+abBPTw7JgAwdQ0dC5gBQ3dhSkUFErNnz9Znor766qulXbt2escm9f8INzezAkHqQNRrc9eN/XI6AywAmMiJUh8WMAM+CSDKy8vjLvfq1UufIwLwKq0ctA6EVDiAoLK7fWMBM+CTAAIwLStABwIA4WQ1qxGcjTYAawgg4EtuZgXYAQMAwitdViNIG20AVhFAwJfczgpQ9gMAMHmnJsBNBBDwJbICAAATBGmjDcAqAgj4FlkBAIDXgrbRBmD8eSAAAAD8zInzTwCmIwMBAACQJUpqEUYEEAAAADmgpBZhQwkTAAAAAMsIIAAAAABYRgkTfIMzfQIAbSoA7xFAwBc40ycA0KYCMAMlTPD9mT4BALSpANxDAAFf4EyfAECbCsAMBBDwBc70CQC0qQDMwBoI+OZMn1u2V8WVMXGmz/BhIT1gD9pUALkggIAvcKZPsJAeoE0FYAYCCPgGZ/oMt1QL6aeOauvZcQF+RZsKIFusgQDgCyykBwDADAQQAHyBhfQAAJiBAAKAbxZ9qoXzsVhIDwCA+1gDAcAXWEgPAIAZCCAA+AaLPgEobOkMeIsAAoFGJwMAwcKWzoD3WAOBwHcyq8orpWLrPllVvktfVtcDAIK3pTMAdxBAILDoZAAgeNjSGfAeAQQCi04GAIKHLZ0B7xFAILDoZAAgeNjSGfAei6gR6E5my/aquFpZzhsAAP7Gls6A9wggEFh0MgAQTGzpDHiLAAKBRicDAABgL9ZAAAAAALCMAAIAAACAZQQQAAAAACwjgAAAAABgGQEEAAAAAMsIIAAAAABYRgABAAAAwDICCAAAAACWEUAAAAAA8M+ZqPfs2SMPPPCAvPnmm3L48GHp27evTJw4UTp37qxvnzx5sixatCjudzp06CCvvfaaR0cMAAAQHGs37ZYFSzfK/oOHpXlJIxk2sJv06NTW68OCwTwPIH7xi19IXV2dlJWVSUlJiTz66KNy4403yn//939L06ZNpby8XEaPHi0jRow4/jsFBQWeHjMAAEBQgocZ89+XvVWHjl+3Zft+mfCzbxNEwMwSpn379ulswm9+8xvp1auXzjqMHTtWKisrZePGjRKJRKSiokJ69uwp7dq1O/6vdevWXh42AABAIKjMQ2zwoKjLC5du9OyYYD5PMxAtW7aUhx566PjlvXv3yrx58+TEE0+ULl26yGeffSbV1dXSqVMnLw8TAAAgkFTZUiJVSa4HjChhirr77rtl4cKF0qhRI3n88celuLhYNmzYoG+bP3++vPHGG5Kfny+XXHKJjB8/Xpo3b+71IQMAAPha85LGCa9vkeR6wKgA4uc//7kMGzZM/vSnP+l1Ec8884wOIFTQUFpaKk888YTOSMyYMUOXNz399NP6NgAAAGS3KHrYwK6yZXtVXBlT6xZNZOjArrykSCovohYaGEQtqB40aJD07t1bpk+frtdJtGrV6vjtH374oQwdOlRnK9R9MjVgwAD9c9myZbYeNwDAfrTZQO6LolVAkGpRtPodteZBlS2pzIMKHtiFCcZmINSah+XLl8vll18uhYXHDkVlFdT6B7WQWv1/bPCgdO16LCLesWNHVgEEAABAGBdFTx2VOIBQwUKy24BEPK0B2r17t9x+++06iIiqqamRdevW6R2ZJkyYoLd0jbV69Wr9UwUZAAAA+AaLohH4AKJbt256UbTaxnXFihV6zcOdd94pVVVVOnBQmQkVXDz22GN6/cPf/vY3ueuuu3SJU/REcwAAAEi9KHrn3q9l/COvy5Syt3XJEuDrRdQPP/yw3spV7ay0f/9++c53vqMXUp988sn638yZM/VJ5n7/+9/rnZcGDx4s48aN8/qwAQAAjJNoUXRensj+6iP6n8KJ4hC4RdROY0EewrrTBuBHtNlA5mIXRavMQzRwiNXnjFKZOuoCXl74MwMBwL6dNphVAgDELopWZUuJAghOFIdccCIFIIA7bQAAoHCiODiBAALwKXbaAABYWROhzgMRixPFIVeUMAE+xawSAMBKOZM6iRwnioOdCCCAAO20wawSAKA+ThQHuxFAAD7FrBIAAPACAQTgY8wqAQAAt7GIGgAAAIBlBBAAAAAALCOAAAAAAGAZAQQAAAAAywggAAAAAFhGAAEAAADAMgIIAAAAAJYRQAAAAACwjAACAAAAgGUEEAAAAAAsI4AAAAAAYBkBBAAAAADLCCAAAAAAWEYAAQAAAMAyAggAAAAAlhFAAAAAALCMAAIAAACAZQQQAAAAACwjgAAAAABgGQEEAAAAAMsIIAAAAABYRgABAAAAwDICCAAAAACW5UUikYiEyNlnny21tbVy0kkneX0oABBoqp394x//mNPfoM0GAPPa7dBlIBo3biyFhYVeHwYAwALabAAwT+gyEAAAAACyF7oMBAAAAIDsEUAAAAAAsIwAAgAAAIBlBBAAAAAALCOAAAAAAGAZAQQAAAAAywggAAAAAFhGAAEAAADAMgIIAAAAAJYRQAAAAACwjAACAAAAgGUEEAAAAAAsI4AAAAAAYBkBBAAAAADLCCAAAAAAEEAAAAAAsB8ZCMBDd955p/Tv3z+nv7F161Y544wz5Pnnn7ftuAAgSFQ7q9rbXMyaNUu3tXb+TcCvCr0+ACDMxo4dKzfccIPXhwEAAGAZAQTgoVNOOYXXHwAA+AolTIDDIpGIzJs3T6644grp1auX/Mu//IvMmTNHX1+/hEn9/3333Sc///nP9X0nTZqkr6+srJSJEyfKBRdcIOeee66MGDFCVq1alfQxv/jiC7n99tulX79+0rt3b/331q1bx3sNILRqamrkN7/5jfTt21e+853v6DZ17969+raf/exn+l+sd999V5csqZ9W/PWvf5UhQ4botvv888+XO+64Q3bu3Bl3n0WLFskPf/hD6dmzp1x66aW6LKq2tlbf9vrrr+vHe+utt+J+Z+XKlfr6999/X1/+6quvZMqUKXLhhRfK2WefLUOHDpXly5fn9NoAmSIDAThsxowZ8vTTT8vIkSPloosuktWrV8uDDz4oR48eTXj/P/3pT/q+//qv/yolJSVy8OBB+elPf6o7mV/96lfSvn17eeqpp+Smm26Sv/zlL1JYGP81Vh3iT37yE2natKncfffd+qd6/Ouvv17+7//9v9K5c2fecwCh88orr+gJlQceeEC3k6odrqiokIULF+b8t9XgfsKECbosVQUoO3bskN/+9rfy7//+7/LHP/5R3+fJJ5+URx55RE8A/frXv5b169frAGL79u164ui73/2ulJaWypIlS+Tiiy8+/rcXL14sp512mnz729+Ww4cP6wmh3bt3y/jx4/X9n3vuObnllltk9uzZepIJcAMBBOCgqqoq+cMf/qA7DDX4V9Ss0a5du2TFihXStm3bBr9z8skn65mrKNX5bNu2TQcLZ511lr6uT58+8qMf/Uj/jfodhgoW1AzVn//8Z+nQoYO+7pJLLpErr7xSHn30UfnP//xP3nMAodOqVSud/S0uLj5++Re/+IW88cYbtgQQTZo0kVGjRkmjRo30dSeccIKeMFLZ5gMHDsh//dd/ybBhw2Ty5Mn6dhUkqPuoy2rSqGvXrnL11VfL/Pnz9cSRmkA6dOiQDnzU31VUcPHxxx/roEcFQ9H2XWVPVECkggnADZQwAQ764IMPdKbhsssui7tedRhqtiiRaJAQ2zF17Ngx7nqVVXj11Vfluuuua/D7KpWt7qsyFeqx1b/8/Hzdybz99tu2PTcA8JPvfe97x4OHaMmoyuCqiZhcqazD119/LYMGDZKHHnpIlx2pAOHf/u3fJC8vT5ecqmBAPWa0XVb/oiWsf//73/XPa6+9Vqqrq+V//ud/9GX1U11WE0bR9r1du3bSo0eP439DZae///3vy5o1a2Tfvn05PxfACjIQgINUJkBp3bq15d+J7eCif6NNmzYZPeaWLVt0B5OI6uRUAAIAYaIG3rHUxIrKQqhMca7U2rSysjK93m3u3Ln6/1WGefTo0To7EO0LopmE+tQ6N+XUU0/Va9dU2ZIKGtRPlbVWE0KK+jsqg52sfVe3tWzZMufnA6RDAAE4qEWLFvqnqrft1KlT3CLnzz77TC/qS6d58+b6XA/1/eMf/9AdRePGjRvcX3VAqh43kWh6HQDCJDqIj1Iz919++aWeoFGTLtHFzFFq5j8Tag2D+qcmad555x1dvqoWbatSo2hfoMqM1HqG+mLLWVUW4q677pJPPvlEZxzU78S27+r3Y6+LpbLVgBsoYQIcpHbjKCoqkv/93/+Nu14tgla7JBUUFKT9G2q3kM8//1w2btx4/Dq1kO6Xv/ylXhRdnwoeNm/eLKeffrreoSP6T9XOqvtbeUwACBpVJhS7eYUqA1WXzzvvPGnWrJle+BwruuuRFf/n//wfPfBX6x1UhleVFKldnqITRiqIUH2B2pUptl1WJVQPP/xw3CTR5Zdfrv/Gvffeq9dBDBw4MK59V4uuVdAT+3fUc1NlsbTvcAsZCMBBqnRJnShOpbXVzL9q/D/88EO9wFllCNQuHOlcc801elHdmDFj5LbbbtMpdzWzpbIXw4cPb3D/G2+8UQcL6qfaqUnd/+WXX9aL7tTOHwAQRqq8R028qJKiTz/9VA/c1c54aiMKNYB/7bXX5P7779frEtQaBlU+ZJXatlWVLqmtudVWrqp9VgN6tUha3aZ+qp2S1EYWakG1ClpUMKEuqzUSZ5555vG/pYIHtdXrggUL9A58sVlj1R+ojTXUomtVHnXSSSfptW2///3v9WYdKkgB3EAAAThM7b6kZoueffZZ3aGoFLPaXlVttao6m3TUzJjqMNR2sP/xH/8hdXV1cs455+gg4lvf+laD8iZVK6seSy3kUzNYKluhUt7Tp0+XH//4xw4+UwAwl5pw2b9/v955SQ3KBw8erNtnNYBX2QNVVqp2u1Ptp1oUrXasUwN4qwu0VVmRyi5HF06rbVdVO62CB2XcuHF6HcYzzzyj+wJVgqqCF5WNVqVJsdQ5IlQAoQKG+mvk1Fbfqn1X28Sq56N221PbxaoJI8AteRGVbwMAAIAR7rnnHp2tziQLAriJDAQAAIABVMZi06ZNuuRUZRgAUxFAAAAAGECtvXjzzTf12abVOSUAU1HCBAAAAMAytnEFAAAAYBkBBAAAAADLCCAAAAAAWBa6AEKdaEX9AwCYjzYbAMwTul2Y1CngAQD+QJsNAOYJXQYCAAAAQPYIIAAAAABYRgABAAAAwDICCAAAAACWEUAAAAAAsIwAAgAAAIBlBBAAAAAALCOAAAAAAGAZAQQAAAAAywggAAAAAFhGAAEAAADAskLrdwUAAICd1m7aLQuWbpT9Bw9L85JGMmxgN+nRqS0vMoxGAAEAAOBR8DBj/vuyt+rQ8eu2bN8vE372bYIIGI0SJgAAAA+ozENs8KCoywuXbuT9gNEIIAAAADygypYSqUpyPWAKAggAAAAPNC9pnPD6FkmuB0xBAAEAAOCBYQO7SusWTeKuU5eHDuzK+wGjsYgaAADAA2q3JbVgWq15UGVLKvOgggd2YYLpCCAAAAA8ooKFqaPYthX+QgkTAAAAAMsIIAAAAABYRgABAAAAwDICCAAAAACWEUAAAAAAsIwAAgAAAIBlBBAAAAAALCOAAAAAAOCvAGLx4sVy5ZVXytlnny0//OEP5ZVXXjl+29atW+XWW2+VPn36yMUXXywzZ86U2tpaT48XAAAACCvPA4glS5bIpEmT5Prrr5eXXnpJBg0aJLfffrusWrVKampq5Oabb9b3e/bZZ+Xee++VP//5z/K73/3O68MGAAAAQqnQywePRCLy6KOPyg033KADCGXMmDGycuVKee+992Tbtm3yxRdfyMKFC6Vly5bSrVs32bNnj8yYMUNGjx4tjRo18vLwAQAAgNDxNAOxefNmHSQMHjw47vo5c+bosiUVSPTo0UMHD1Hnn3++HDhwQNavX+/BEQMAAADhVuh1AKFUV1frUqV169ZJx44ddRaif//+smPHDjnxxBPjfqe0tFT/3L59u/Tu3duT4wYysXbTblmwdKPsP3hYmpc0kmEDu0mPTm15EQEAgC95moFQmQRl4sSJeu3DU089JRdddJGMHTtWli9fLocOHWpQptS4cWP98/Dhw54cM5Bp8DBj/vuyqrxSKrbuk1Xlu/RldT0AAIAfeZqBKCoq0j9V9uHqq6/W/3/WWWfpTMTcuXOlSZMmcuTIkbjfiQYOxcXFHhwxkBmVedhbdSjuOnV54dKNMnUUWQgA8AMyyYBBGYj27dvrn2pxdKwuXbro7VtV+VJlZWXcbdHL0d8FTKbKlhKpSnI9AMAsZJIBwwIItUC6pKREPvzww7jrN2zYIKeccor07dtXZyOipU7KO++8o3/nzDPP9OCIgcw0LzlWcldfiyTXAwD8k0kGwsrTAEKVKN1yyy36vA5//etf5bPPPpPHH39c/v73v8vIkSNl4MCB0q5dOxk3bpx8/PHHsnTpUnn44YflpptuYgtX+MKwgV2ldYsmcdepy0MHdvXsmAAA1pFJBgxbA6GoBdNNmzaVRx55RHbu3CmdO3eWWbNmyXnnnadvnz17tkydOlWGDh2qt3MdPny4/h3AD9RuSxN+9m09U6XKllTmQQUP7MIEAP5AJhloKC+izuYWIgMGDNA/ly1b5vWhAADSoM2GKWsgYsuYVCZZTQ4xGYSw8jwDASAz7AYCAO4hkww0RAAB+HwmbMv2/cyEAYDDQQRbbwOGLKIGkBl2AwEAAF4jgAB8hN1AAACA1wggAB9hNxAAAOA11kAAPlow3a97e9myvarBbiBenFeCxdwAAIQTAQTgswXT1w3oIivWVXp6XgkWcwMAEF4EEIDPFkyr4GHqqAvE1MXc7FQCAECwsQYCMJTJC6ZNPjYAAOAsAgjAUCYvmDb52AAAgLMIIABDDRvYVS+QjuXVgmk/HRsAAHAWayAAQ6mF0RN+9m29rsDLBdN+OzYAAOAsAgjAYGpAbuqiZJOPDUA8tl0GYCcCCAAAAoxtlwHYjTUQAAAEWKptlwEgGwQQAAAEGNsuA7AbAQQAAAHGtssA7EYAAQBAgLHtMgC7sYgaAIAAY9tlAHYjgAAAIODYdhmAnShhAgAAAGAZGQgAAOAYTmIHBA8BBAAAcAQnsQOCiQACAAC4fhK7qaPa8qrD19Zu2q0/4+pcK81LGsmwgd30eqMwIIAAAACO4CR2CKq1m3bLjPnvxwXIW7bvlwk/+3YogggCCCBAwjwbAsA8nMQOQbUg5Nk1AgggIMI+GwLAzJPYbdleFdcutW7RRIYO7OrpcQG52n/wcMLrq5JcHzRs4woExFMvrE06GwIAXp7Ers8ZpdKlY0v9k0kNBEHzksYJr2+R5PqgIQOBwPJrOU82x61+55Nt+0I9GwLAG3996xN55tUNcrjmqDQuKpDhl58hgy7ufPx2TmKHIBoW8uwaAQQCya/lPNketwo4ausioZ4NAeBN8FC2eI1E/tn8HKmp05eV2CACCGp2beHSjXqiTvW1KngweYxhJwIIBJJfFzdle9zJajEL8vNCMxsCwH0q8xANHqLU5T+/uoEAAoHXo1Nbo8cUTmINBALJr4ubsj3uZLWYnTu0DM1sCAD3qbKlRA4luR5AMBBAIJD8urgp2+NWtZiq9jKWunzTkB62Hh8AxGpclLiQoUmS6wEEAwEEAinZgNr0cp5sj5udTgB4Yfjl3SQvL/46dfmnl3fjDQECLC8SqV+9GGwDBgzQP5ctW+b1ocCFBcl+XNzk1+MGnECb7Y+F1GrNgypbUpkHFTwEZQF1ul3x/LrbH5ArAgiEAo084E8EEDBpVzyVEY7uipfudiDIKGFC4EUb+VXllVKxdZ+sKt+lL6vrAQDIdFc8K7cDQUYAgcCjkQcA2L0rnl93+wPsQACBwKORBwDYvSueX3f7A+zAPmsIPBp5AEA2u+Jt2V7VYI1DdFe8dLfXx1o8BInnAcTOnTvlkksuaXD9/fffL9dcc41MnjxZFi1aFHdbhw4d5LXXXnPxKOFnmTbyAABEt8dOtiteuttjJVpwvWX7fhZcw7c8DyA+/vhjady4sSxdulTyYjaTbt68uf5ZXl4uo0ePlhEjRhy/raCgwJNjhT9l0sgDABDbf0wd1Tbj2+tnGw5U1yRdcJ3q7wOm8jyA2LBhg5x22mlSWlra4DZ1ioqKigoZNWqUtGvXzpPjQzg6AQAA7JAo21CQX+9se//Egmv4lecBhMowdO6c+IQzn332mVRXV0unTp1cPy4AAIBMg4fpc1fI/uojcdfX1iU+Zy8LruFXRmQgWrVqJddff71s3rxZTj31VBkzZoxeF6FuU+bPny9vvPGG5Ofn6+vHjx9/vMQJMFEQFssF4TkAgNuZh/rBQ2wWIjaQYC0e/MzTAOLo0aOyadMm6dKli9x5553SrFkzeemll3TJ0ty5c3UAoYIGVd70xBNP6IzEjBkzZOPGjfL000/r2wDTBGGxXBCeAwB4fc6hWJ07tpRmTRuxFg+B4GkAUVhYKO+++65eFN2kSRN9Xc+ePXWAMGfOHCkrK5Phw4frDIXSrVs3vRZi6NChsnr1aundu7eXhw9kfOI6v6zDCMJzAAATzjkUzTbcNLiHpxMwZJURqBKmkpKSBtd17dpV3nrrLZ1hiAYPsbcpO3bsIIAIET81fEE4cV0QngMAGHHOoeJGnmdvySrDbp7WAKlMQ58+fXQWItaaNWt0WdOECRPkxhtvjLtNZR4UdTvCIdrwrSqvlIqt+2RV+S59WV1voiCcuC4IzwEA3D7nkMo0xFKX7xrZVwcPqs+aUrZcxj/yukwpe9vVPixVVhnwXQZC7b6kdliaNm2aTJ06VWcbFi5cKB988IE899xz8vnnn8vYsWPlsccekyFDhuhF1uq+gwYNSrpzE4LHb+U0QThxXRCeAwBrWd3oOZjU1ummZ3j9es4hrzMAZJURqABClSipxdEPPfSQjBs3TqqqqqR79+56AbVa76D+zZw5U6+F+P3vf693Xho8eLC+L8JBNboff7rXV+U0QThxXRCeA4DEEg1mY4V5w4Rcy2WTnXPI64kwssoI3BqItm3byv3335/09iuuuEL/Q3g7ua8PH/VdOU0QTlwXhOcAIPPdgkzO8DrJySyB1xkAssqwG/ugwpedHOU0AGD/bkGmZ3j9uk7A6wxANKvc54xS6dKxpf4Z1iwTApKBADLt5Jo2LqThAwCbB7N+yfA6xcksgQkZALLKsBMBBHzXyZ11WmtmTQDAxsFsrLBmeJ3MErCuDEFDAAFjmTBjAwBBU38wm5+fJxIRqYtEQr1hgtN9DhkABEleRO3bFiIDBgzQP5ctW+b1ocDiojY/7QTkpxPeAX5Amw03+a3PAbxCBgJGs3PGxunBvdf7fAMAgpUlYFIKpiKAgK9ZbVwTDe5XV+yW9q2LpbR1sS3BhNf7fAMAgiPdpBTBBbxEAAHfymTGP9Hg/mhtRLbtOqj/2ZEpSLWDR6YNPR0DAIRbqkmpoQOFjDc8RQAB38pkxj/dvud2ZAqS7eCRn5eXUUNPKRQAZCeTyRfTJ2pSTUol6/+eenGtPPT/fc+lI0SYcSI5hGLPbiv7nue617fawUPt2BErejmTkxM5eTIjAAiq6OTLqvJKqdi6T1aV79KX1fW53NfEbWWT9X+fbN1n1HNAcBFAwGiqIZxStlzGP/K6TCl7O65hzGTP7kSDeyu/Z8eZPtXWiJkELE6ezAgAgiqTyRc/TNQkm5RSO0Ml6/9q6yJGPQcEFyVMMFa6Up5M9uyO3fd8556DsvPLr+VobV3a37NjB49MT06U6v6mp9wBwCuZTL74YaIm1cnnhg0U+WjjLh0wmPwcEFwEEPDtGodMz+wZO7h3c6/vTE9OlOz+fbuXsmgOAJLIZLLGybNOu7GtrLq+U4eWsvHzr4x/DggmAggYy8oMUbZ7dru513c2gU6i+7NNLAA3eJnpzOWxM5mscfqs0264eUiPBpNKfnsO8C8CCBjLLzNETgQsie6//+CahPfdufdgzscHAF7vApfrY6v7XDegizzz6gY5XHNUGhcV6MuJfjfTiZ3YYzSljDTb5wDYgQACxgrCDJEbAdXOvdW6U6PTAJArLzOduT62agcXLauQ/dVH9OUjNXX68uknt0waRGTynEzcYtu0M2cjPNiFCcZKtqtRWAfKKqAqLGj4lVUnxLO660aqXa0AwMvFxbk+ttM7K/lh5ybALWQgYDRmV+Jfi/atmsq23Qez6mBNnD0DYBYvS0dzfWyngx+3giuTyqSAZAggAA9k20GUtilJGEBY6WCtlAfQcQHBl+p77mXpaK6P7XTw40ZwxUQP/IIAAnBZLh1ELh1sutkzOi4g+NJ9z71cmJvrYzsd/LgRXLHbHvyCAAJwecZPLXreX12T1ULBXDrYdLNndFxA8Fn5nntZOprLYzsd/GT697PJ6PrhBHeAQgABeDDjl0sHkW0Hm272jI4LCL6gf8+dDn6s/v1sM7pWyqQoNYUJCCDgurA1folm/LxYpJhu9sy0826E7XMCuMG073lQZZvRTTfRQ6kpTEEAAVeFsfFLNuMXy61Fiqlmz0w670YYPyeAG0z6ngdZtpmedBM9dpSaMjkDOxBAwFVhrLNPNuMXVViQl/RsqWE9q2kYPydA2L7nQZZLpifVRE+uJWhMzsAuBBBwVdDrb63O+NU/EdyKdZUy6OLO4jVTzrsRxs8JELbvuV1MnFF3KtOTawkakzOwCwEEXBXG+lvVkakMw5wX1upgIREGxvHC+DkBkDlTZ9SdyvTkGpgkm5zZubfh+YWAVAggELj6WxNno95bV5k0eFAYGMejThuA32fUncj05BqYJJucUduLq77T674S/kEAgUDV35o6G5VqITULGBuiThtALm2raVldOye2cglM1OTM6ordcrS2Lu56NcFlQtAF/yCAQKDqb63ORrmdpUhaklPcyPPgxlRBq9MGEM5yR5MmttTjtW/VVLbtPmh80AWz5Xt9AIDbs1HRxnxVeaVUbN0nq8p36cvqeqeoWR+VaYilLt81si/BAwDY3LaatC1tqoktL5S2KTE+6IL5CCAQutkoLxrzaElOnzNKpUvHlvonmQcACH7balqZlR+CLpiPEiaEbvGtV42530pyTFyMDgB+a1tNK7NijRnsQACBQLHSMJrWmJvIpJpdAPAzE3eVMz3ogvkIIOB7iWbKp466wFeNuWlM3hoRAPyEGX8EEQEEQjdT7lVj7qeSINNqdgEgyBNbYexn4G8EEAjlTLnb6Vu7SoLc6hxMLfOicwRgOq9KQCk9hZvYhQm+5peZcjt2fnJz+1kTd+nwYvtdAPDLtq1uPK5qb6eULZfxj7wuU8repv0NMTIQ8DVTZ8qdCHTcXJdgYs0u6zKA4AtCltGriS2nH5cMB4wKIHbu3CmXXHJJg+vvv/9+ueaaa2T9+vUyffp0WbNmjbRu3VpuvPFGueGGGzw5VpjHLwuiswl06neklXuqXe2UTNulwy/ZJgDhHqB6NbGV7HHz8/N01iDXoIxJHBgVQHz88cfSuHFjWbp0qeTl5R2/vnnz5vLll1/KyJEjpX///jJ16lT54IMP9M+SkhK59tprPT1umMHLmfJMZsoyDXQSdaSFBfk5ByF+nM3zW7YJCLts252gDFC9mthK9LgtSopk++6Dsr/6q5yDMiZxYFQAsWHDBjnttNOktLS0wW1PP/20FBUVybRp06SwsFA6d+4sW7ZskbKyMgIIF/hl8OnFTHmmM2WZBjqJOtKjtXVSWJAnR2sjljqlv771icx5YW3c/f04m+e3bBMQZrlkEYIyQI22948t+lC+2HVAVAv89eEa2fzFPkfb3kT9zP7qI7Lx82+Ch1yCMquTOH4ZO8DnAUR5ebkODBJZuXKl9OvXTwcPUeeff748+eSTsnv3bmnblg+kU4KSSrYq0wYv2UzZUy+slZLiRgn/TiaBTrKOtH3rYmnfuiRtEKKez5wX1sUFD3bO5rnVQdR/nOsGdJEV6yqNWZcBwL4sQpCyjCpY2KaCh382wV8frpWyxWv0/w+6OPGYxw71+xm12NmuoMzKJE7Yxg5hZkQGolWrVnL99dfL5s2b5dRTT5UxY8bodRE7duyQbt26xd0/mqnYvn17qAMIpwdwQUklW5FNg5dsgP/Jtn1SW5f7jH9sOV+skiZFei/x6Ps/e8mahO+/uk1lLJyYzXOrg6AjAvwnlyyCm1nGXPpQK7/7zKsbjgcPUeryn1/d4GgA4WRQZiWTHqaxQ9h5GkAcPXpUNm3aJF26dJE777xTmjVrJi+99JKMGjVK5s6dK4cOHZJGjRrF/Y5aL6EcPuyvlKad3BhYBSWVbEU2DV6yRjk2eLDyd5LJS3GDlfc/2ftnx2yeWx0EHRHgP7kMWN1a05ZLH2r1dw/XHE34+4eSXO8Uu4OydJn0MI0dws7TAEKVJr377rtSUFAgTZoc23O+Z8+esnHjRpkzZ46+7siRI3G/Ew0ciouLJazcGFgFKZWcTjYNXqJGuSA/r0EAke7vJFNXf+oqen1dxNL7n+z9U2socp3Nc6uDoCMC/CfXAasba9py6UOt/m7jokI5UhM/flGaFBUGeqORMI0dws7zEia1o1J9Xbt2lbfeektOPPFEqaysjLsterl9+/YSVm4MrIKwYNVqijqbBi9Ro3yg+ohsqLdYLd3fSSbVMSV7n9X10edcufdggwXXahenm4d0z7njcKuDoCMC/CdR29i3e+k/2+LEJZd+6kOt/u7wy7vpNQ+xc0GqMvWnl8eXZQdto5EgjB3ggwBCZRqGDRsmjz/+uJx33nnHr1fnfFBlTWeddZY8++yzUltbq7MUyjvvvCOnn366tGnTRsIqk4FVtnWeJp5IzKkUdbYNXv1GOdFjZttw9uteKh9sqGzQ+aiO+L118UF1VH5eXsKtX9u3airt25TY9v651UHQEQH+FNs2mriWKVkfum3XQX125VT9pNX+N7rOQa15UGVLKvOgggc31z94IduxAzs3+U9eJJKkVsIFdXV1MnToUPn666/1+R3UYuqFCxfKM888I88995wOEq644gp9HohbbrlFPvroI7n33nv1fa+++uqsHnPAgAH657Jly8Svkg1U6zfIVu8XROqkOavKGw60+5xRqhch16deKzuCJbv+TqrjV38z0fvatmWThBmQZM85F3Y9T1MeB+YKQpsdZpm2xW5I1DfGStVPhrlfdQqvqT95moHIz8+XJ554Qh566CEZN26cVFVVSffu3fUC6ujuS7Nnz9ZnolYBQ7t27WTChAlZBw9hi/CDsgg1m5mJbFLU0Ug6cvz/vEsVpzr+ZO+/2pEp2e/4NSVu2tmwAfh/LVNsG7r+073y9eGjDfrJ6XNXyKSRfRv0NX7PzpsoKGOVsPF8DYTaivX+++9PenuvXr1kwYIFrh6TH1gZWJnYcGfq2MnQ1sVtSWol/Z1pmZdfUuzR40/0/rNmAIBpEz6qtDLbciEnRdtQdZ6Eiq37GtyuTsA26fG3pVOHlnLzkB5xx+jkxEYYS3mCMFYJo3yvDwDO8fuA8tjJ0NY2OJ9BdGYiXf28SivHSlann2r2wyuZHH8uvwMAdolOxqiSJTUoX1W+S7bvOSgtSooa3FfN+qvb1f3V75nWTypqVz11Fme3jjHR6+f16+MGv49VwsrzDASc4/dFqMdOhhbJamYikzSzibMf2aTJ/ZZaD+NMGxBkiSZj9lfXSLdvnSBdOjZKWi7kZalKon6yvkyOMZd2zW+lPHa14X4fq4QVAUSA+W1AWV+uJ0OzmmY2dfYjmzS5X9YMmFg2BsCZNlud10YtmE5WLmTCZI1a86DKlpKxcozp2rX6A+5+3dvrXfWilyv3VGf92H5uw/0+VgnrpBwBRMD5ZUCZaWpZbWdqF2Y/3Oe3mTYAkvNkjMmTNWrBdKqdmawcY6p2behAafD3P9iwK26rbrXtdraP7fZg1O423M9jlbBOyrEGAsZSA3t1MrREViQ5F0I2orMfalvBLh1b6p9+nglXDYnaOlHN9qlFiibWz5pYNgYgN+nWYSW6XbXxO/ce9LytivYDqtyqID++37FaTpOqXUs04K6/ib5a71e/z/OilMfKWgzacLN4sZaTDASMpRr09q2L9W4dTg80gzL74ZfSIFNnIgFkL10pSuztO/cclJ1ffq0HzaqNV/+8bqvU4z407ntZn38mVbtmtc9SfV771iWelvJYyS7QhpvFi4COAAJGK21dkjCAYKDZkOr0EtXxmlgaRNkYEEzpJmOit6ss6bbdB3Nuq5yo+65/Ju1jf39N2r/fr3uprK7YFbf5RzSDoP6GFSp48OoEe5kMRmnDzeJFQEcAAaPRSGWWeUi2CNC00iAWzQHhZseMabYZV6tBRyZ/X9130bKKuOBBrWm4bkAXfd9hA6XBTkPqFBmxZUym7DxkZTBKG24WL8ZKBBAwGo1U9iln0zM2QSkbA+DNjGk2C3kzCQoy+fuJ7qvKs9R6vUEXd07Yl6nNQNTtpu08ZHUwShse7rESAQSMRyOV25a3psxqAfAnJ8qE7JgxzSaLkUlQkMnft3LfRH2ZCi5M26qTiTt/6uHypBwBBHwnUQOqZNOoOtEYe3GCtKSzecWNLC1K5KRuANzcmCHTQWqiNiqbLEYmQUEmf9+NGnQ3N8lg4g7pEEDAVxI1oBWff6VrSasO1mRcB2t3Y5zsb6o62NgTBtkdVCSbzbMaPPhh5yYAwTpni9VBaqp2NdMsRiYD/UTtauy2s7HtuBs16Jw/BybhPBDwlUQN6P7qmrjgwer+x07sm5zsbz75lzUp99T28lwWXuwfDcAfTNjvP1kbpdYPZNrupTtXRbJ2tUPbEr0oWi2SVjsD1m/H3TifkAnvBRBFBgKBqfXPtFF1ojG2enxObK2a7daD2bwOlDwB4ZBsxn7n3mrdDriRpUzVRmVaapNp6VQm2846XfbDuRdgEgII5MzNwWSyBjSb2lMnGuNMjs+pWaNMS5IyfR0oeQLCI1FpTjTzq9oZN0od7W6rsxnomzD7z7bmMAkBBHIKAtweTCZqQIsbF8jXh2slZjttvSZCbZGX6d/KtWY1WWebyM69Xzsyg5dpnWymr0Oyv//Ui2ulpGkjVxePA3BWdMbey5NUJmqjCvLz9PG4lQUxYfa/fvYkX3V0IjJ7SfpMM2A3AgjkFAS4vagrUfpZdSIbP/8q7n7q5DzR/bfd3Dc59m+u/3SvfH34aNL7quN2YgYv05myTF+HZH//k637pLbumzCOhdhAMKi2oH3rpglPVOnGDHy0jXrqhbXyybZj7Yz6p9p9t7IgViZa3NpiVfWtZIKzQ/mtfQggQiz6Rfo4wUDXahDgRVq3fvp5/COvZ30MTtSspmrg63Mi2MpmpiyT1yHZ348NHrJ5bjTsgLm8noFXbVRJcaOc25lcHj/VRIsdA/pM2kB2ZMocQZfHAcThw4elcWPzzmqLzFgZ3FoZgHvdqaQ6BlUipIILr1K79TsctXNHooyE3cGW03WyycoJ6nfsmTw3pxt2ghPA//X3Xq9DSDXRkuuAPtM2MJPXgvbPnvcIOQYQF110kfzwhz+Ua6+9Vnr16pXpr8PgL1I2QYAJnUqiY1CloSrdHk25e1VOE9vhqF081FauTgdbTp9F1GoZWSbPzcmGnVknIBhnJzZhwiqZdAN61Q7NeWGtbK3cry93LG0uNw/pcfz1y7QNtPpa0P5Zf4/gcABx0003yZIlS2ThwoVy+umnyzXXXCNXXXWVtGvXLtM/BQ+l227UahBgQqdS/xhU5sGrxX6mBFtObydY/+8n6qQyeW5ONuzMOgHBODuxCRNWyaQa0Kv2cfrc9/TOVVFqwuW+ee/JXTf2069rpm2g1dci1aYXD/1/35Nc+C2zYXIAGooAYuzYsfrfP/7xD/nLX/4iTz75pMycOVMuvPBCnZXo37+/FBUVOXO0cPyL1LRxoZx1WuuMggAnOxWrDVTsMaiyJa8W+3kRbJnQiOf63Jxs2Jl1AoLBhAmrZFIN6HX7HBM8RKkToEYntjJtA62+Fqk2vchlBys/ZjZMDkBDtYi6T58++t/dd98tb775psybN0/GjRsnLVq00FmJESNGSIcOHew9Wjj+RTLpy59tA2XyLIPdwZZJjXguz83Jhj3Z50FtgajKyvwyewbA2yxIqsmaVAN6dULPZKITW9m0gVZei1SbXuSSlfdjZtfkADR0uzBt375dlzO98sorUl5erkuaLr30UnnjjTfkmWeekfvvv1+uvPJK+44WofoiZdtAZdoQmzCDny0/NuJufx4TfR6aFxfJjr0HpSpm3Ybps2eAH/m5fc10sibZgD7VCUajE1tOtYGq/fto466cNrkIUmbX6zK8IMk4gDhw4IC8+uqrsnjxYnn//felSZMm8oMf/EDuuecenZFQJk6cKLfeeqvcd999BBAGM/2LlG0DlUlDbNIMfpgacTc/j4k+Dweqj8iGeou+/Rh4ASbze/tq12SNGsRXfP5lgzKmFiVFcRNbTm0r3qlDy5w2ufBbph8G78J05MgR6d27t0ybNk0HCMXFxQ3ud/bZZ8u6devsOk6EUC4NlNWG2O8z+DTi7p87BIA1JrSv6TIg9W/v1729vLeussH9c5msUb8/aWQ/fSK8zysP6Ou+1b6Z3DT4m12YnKR2e8plkwsnyk6DkpkKs4wDiOuvv15+/OMfS6dOnVLeb+TIkTJmzJhcjg0hl20DlUnD5PcZfBaFZYfAC3Ce1+1rugxIots/2LBLIjHVPtH759pmqMd7aFxuux5ly4nyqFz+ZpAyU2GWcQAxYcIES/crKSnJ5ngQQNnONGTTQGXaMOXSKZgwg+KHtSwmIvAC/BOoZ9vWpsuAJLo9NniIvb/f2wynyqOif/Ob92iN52fRNqFvDoOcFlEDTs80JGr0UjUOmTZMuWQ5rDwvNxoy09eymIjAC3CeHYPuXNradBmQdOdDir0/bUbu75EbmSmyG+4hgICj7J5pSNc4ZNowZdspWHleNGRmI/ACzM+Q5tLWtmnZJGUGJNXuSInuT5uR/XvkVgmpCetuwoIAAo6ye6YhXeOQTcOUTadg5XmZ1pCR1gXgdtuR66A7l7a27QlNdMYjWQYkUYYkLy++jMlPZUpeceos2m4cC7JHAAFH2T3TkK5xcKtO1crzstqQpdoFRHVmeZIndZGIpY47WUefaIau4vOv5KQ2JZb/NoDw8TqTmktbW1cXSZkBUT+vG9BFnnl1gxyuOSqNiwrk+9/uKFsrD7KmzOb3KJaT5WBskOEeAgg4yu4BfbrGwa06VSvPy0pDZmUXkFipOu5UHX2iGTq1J/n+ak6kBviRWxlFrzOpuba1qTIg6jVctKxC9lcf0ZeP1NTJWx9uZzcgB96j+pwqB/N6sfvaEC3gJoCAo+we0FtpHNyoU7XyvPp1L5UPNlTGBQMqo9C3e+nxy1Z2AbHacafq6K0sFqROFPAHN7MCXpeEWGlrsx00eh0cBUX0PXrqxbXy+c5j57lo09KbE8p5udh9bci2pyWAgOPsHNCbtBNGuuelypDqBwPq8op1lTLo4s4Z7QJipeNO1dFbXSxInShgPjcHvqlm992abU3X1mbbL3gdHAXN7q8OydeHj+r/3/j5Pj2Y9mLw7NVi9wUhC0gJIOA70cYh2nnNXpJ+32kvWOmcrA7srdSVpurohyaYocvkbwMwh5sD32Sz+yqTatJsazaDRurl7RO2wXMiYQtICSDgS7mkCt2aNbPSOVnZBSRWqrR8qjR+/Rm6/Lw82bFXLRSssfS3AZjDzYFvstn9IAwYva6XD5KwDZ4TCVtASgABX8q283KzRtHqeo36nbOa2VNlTtGBvuQd200kXVo+XRq//gydei1MKAUDYPbAN9HsvjrjsN8HjCaVxPpd2AbPiYQtICWAQKhmO6wGHnZkKax2Tok65+gaCSfT+KaeFClMu1gAfh34BmXAaGo76Lf2N2yDZ1O/l24igIAvZdt5WQk87MxS+K1z8nrwHrZdLIBsed22ZDtg9LqNCSsrr3su7W/YBs+mfi9DG0Bs3rxZrrnmGrn77rv1T2Xy5MmyaNGiuPt16NBBXnvtNY+OEibItvOyEngEobbXr+tKwvraA36TzYCRCQJvWH3dc21/wzR4hkEBRE1Njdxxxx1SXV0dd315ebmMHj1aRowYcfy6goICD44QJsl2tsNK4BHWxWAmrCsJ62sP+FGmA0YmCLxh9XWn/YUvA4hZs2ZJs2bN4q6LRCJSUVEho0aNknbt2nl2bDArFV3/cW+5qqflx7USeASltte0dSVWhPW1B8LQDzBAzf61y4XV1532F74LIFasWCELFiyQxYsXy6WXXnr8+s8++0xnJDp16uTp8cGcVLQdj5tu1iysi8GcXFdiVVhfe8CPMm2PGaBm/9rlwurrTvuLTOSLx6qqqmTChAl6rcNJJ50Ud9uGDRv0z/nz50v//v1l4MCBMm3aNNm/f79HRwsrs85+f9xolqLPGaXSpWNL/TMMi3hV56EG67HsWldiVVhfe8CPMm2Ps21jgsjNPtTK6x7NhjRtVCDNi4ukQ7sS2l+YnYG499575dxzz5XBgwc3uE0FEPn5+VJaWipPPPGEzkjMmDFDNm7cKE8//bS+Dd6UEnmVinbqcRO9hlNHXSBhksu6korPv5T91d+clE51QNkOCliIB/hDpu2x2zv1mLzjk5t9aLrXPVE2pKiwIJS7KMEnAYQqWVq5cqW8+OKLCW8fM2aMDB8+XFq1aqUvd+vWTa+FGDp0qKxevVp69+4tYedVKZFXqWgnHpedQXIfvKvz3aW6DCB4smmP3ZogSNauXzegi7y3rtLzoMLtPjTV687idmTD0yn85557Tvbs2aPXPagshPqn3HPPPXLLLbfoDEM0eIjq2vXYrOaOHTs8OWbTeFVK5FUq2onH9eo1DAr1+lUd/Cb7oKjLvH6A99RAekrZchn/yOsypextfdkuJpckJWvX57ywVlaVV0rF1n2yqnyXDjLsfE38+NqxuB2+y0A8+OCDcuhQ/Bf8sssuk9tuu02GDBmi10ZUVlbKvHnzjt+uMg9Kly5dXD9eE3n1xffqpDFOPG5YG0+70vthff0A0zmdXTX55GHJ2qWjtREjzjNj0mvH4nb4LoBo3759wuvbtGmjb7v88stl7Nix8thjj+mAQp1oTi2iHjRokHTu3Nn14zWRl198r2rV7X7cMDaedg4swvj6AX6YBHCjNMXUNUvJ2iWTJjtMee3YfQnZMHoV8oABA2TmzJmybNkyvch60qRJOkNx3333eX1oxjApDepXYXwN7SzbCuPrB5gyCZCqHCfM2cFE7VJhQeIhT9gnO9j9Dr7chak+debpWFdccYX+B/PToH4VxtfQzoFFGF8/wGtWsgthzg4mapf6di+VRcsqOM+MwdkQ+IdxAQQyxxffudfQ5G0Ac2H3wILPIGDeJEDYS1MStUunn9ySyQ6PpepXg9rnBhEBBEIrXUMV5O1dwz6wAMIwCRCm7KDVgaeb28h6ORD2+vFTHVeyflUJap8bRAQQCCUre4Tv3Pu17K8+YsSOHXYL08ACCCKrkwBhyA6aNtnj9fF4/fjZlt6p/bGcXvQP+xBAIJSS7xG+To7W1qX83aAsQAzDwAIIKiYB0rfn0+eukPatm+Y8A5/pbL7XJ2bz+vHtXn8XlD43aAggEErJ9whPHTyEZQEiAPN5PQlgSplM5d6DCa9XGeRoFjnbs1BnM5vv9e5XXj9+tqV38WfoiL8N5iGAgJEdgkl7hMdinQAA5F4mE+1rKvcclKrqI9KipJGUti7OuM9Rf2fn3uq091PH+eRf1sRdZ+V4s5nN93r3K68fP5fSO9bm+QcBRADYNeg3uW7S7ueaqBErLMhrcJZSpUWx6tiask4AAGwok/nrW5/InBfWxrW3+6trZNuugxn3Oeo4ErXbVlg53mxm873epMLrx8+l9I61ef5BAOFzdg76Ta6btPu5ZrJHuEkBFACYINsyGdWOH1trlnjQn2mfk+w4rEp3vNnM5nu9PsXrx8+l9M7rsjxYRwDhc3YO+k2um3QiwEm3R3h+fp5EIiKzl6wJdDkXALhVJnMsY2DfRhXJjiNZRjnT4812Nt/rgbDXj4/gI4DwOTsH/SbXTboV4EQbXT+UcwGAV7IdWFvJGGTS5yQ7DrVgesW6St0/JNqSWyksyLcUCJg8m++lsKyZRGIEED5n56Df5LpJtwOcXLIdNKoAgi7bgXW6DSwy7XNSHcegizvr+ySaEFIZipuHdLc04PVyNt/U/oRJNhBA+Jydg37TZ1rcDHByqe8Na+bC1I4OgDOyGVgnaseVtic0kVPat8iqz0l3HKb3bX7sT0xfMwnnEUD4nN0No8l1k252ArnU94axUTW5owNgDq8G8271bXZOpJjcn5i+ZhLOI4AIAJMH/X59rnbX9wa9UTW5owNglqD2WXZPpJjcn5i+ZhLOI4BAaKWaKbK7vjfojarJHR0A+HEiJT8vL+H1h47UypSy5Z6Wi5q+ZhLOI4BAKFmZKUo3S5YoAMmkUQ3SmoGwBk5AmAWpDTNxIiXZJrRf7DogWysPeFou6td1JbAPAQRCKdeZolQBiJVG1bQ1A7kOBJiNAsLFpDbMlEDG7omUiDoRUQJ19a72qlw0qKVosIYAAr6VS6eR60xR6gDkgrSNqklrBuwYCDAbBYSLKW2YSYGM3RMp6ba8jUW5KNxGAAFfzhLl2mnkOlOUawBi0poBuwYCzEYB4WFKG2ZKIOPEREqigESd/C7RmbwpF4XbCCACkDb1A7tniXLtNHKdKco1ADFpzYApAwEA/mFKG2Za+2XnREqigKRv91JZtKyCxcvwHAFEANKmfmD3LFGunUauM0W5BiAmrRkwZSAAwD+8aMMSTdoFvf1KFJCcfnJLFi/DcwQQAUib+oHds0R2dBq5zBTlGoCYtGbApGAGgD+43YYlm7S7bkCX0LVflIvCBAQQAUmbml62ZfcskQmD3lwbcVM6AZOCGQD+4WYblmzSbsW6StovwAMEEFkKetrU7rItuwf8Tg96w7a+xZRgBgAynbQzpf0KW7+BcCOAyJIJM+B+KttyYsDvVKfB+hYAMIvpk3b0GwgbAogshbHsw46Fy5kO+L2Y0WF9CwCYxfRJO/qN4CGjlBoBRA5MSZsGdQbIqxmdoK1voREE4HemT9oFrd8IOzJK6RFAhEwug0m3Z4C8mtGxK1AyYeBOIwggKEyetDO9xAqZIaOUHgFEiOQ6mHR7BsirGR07AiVTBu40ggDgPNNLrJAZMkrpEUCEiB2DSTdngLya0bEjUDJl4J6sEVyzabdMKXubXUIABIpXmV/TS6yQGTJK6RFAhIjfImovZ3RyDZRMea2TNYJHaupkVfkuzp4OIDC8zvwm6jdMKGVF5sgopZdv4T4ICL9F1NEZnT5nlEqXji31T7dLgPz+WqtGUAVdyUSzIgDgd6kyv14GNKvKK6Vi6z49aaMuq+thNj+PP9xCBiJE/BhRm7xozg+vdWxaXZUtqcyDXzJQAILF6dl4UzK/ppWyIlzjD7cQQIQINZrhfK2jjeCUsuV6JswvGSgAwQkI3CgvsiPza2eQY1pAA9iJACJkiKjD+1qbkhUB4H+ZBgRuzMbn2sbZHeSYUsoKOIE1EEBIUNMJwKv1Bm7Mxufaxtm9hiLRGjQmbRAUZCCAEDEtKwLAnzINCNyajc+ljbM7yHGylJXdneA1AggAnqETBPwp04DADyWUTgQ5TkzaeL1dLaBQwgTAE2xxCPhXpuU5fiih9EvJkWnb1SKcjAogNm/eLOeee648//zzx69bv369jBgxQs455xzp37+//OEPf/D0GAHYg04Q8K9MAgI1WaB2gZu9ZI1EJCK3XNVTpo66wKjgwS9BjsLuTjCBMSVMNTU1cscdd0h1dfXx67788ksZOXKkDhymTp0qH3zwgf5ZUlIi1157rafHC2SLsp1j6AQBf7NSnuO3chs/rBNjdyeYwJgAYtasWdKsWbO46xYuXChFRUUybdo0KSwslM6dO8uWLVukrKyMAAK+5LfO1El0gkDwcTI1+/lhPQmCz4gSphUrVsiCBQvkgQceiLt+5cqV0q9fPx08RJ1//vny6aefyu7dnAoe/kPZjv/qjROJlmSMf+R1mVL2tr4MoCEyjeEttUKweZ6BqKqqkgkTJsjkyZPlpJNOirttx44d0q1bt7jrSktL9c/t27dL27Z8WeAvdKZmnq07E2SRAOvINIa31ArB5nkAce+99+qF04MHD25w26FDh6RRo0Zx1zVufGw7tcOHORU8/IfO1P+dICUZgHWU2wDB5GkAsXjxYl2m9OKLLya8vUmTJnLkyJG466KBQ3FxsSvHCNi5SNoPnSmLvFMjiwQEP9MIwOAA4rnnnpM9e/bIpZdeGnf9PffcIy+//LKceOKJUllZGXdb9HL79u1dPVbAjvIW0ztTynPSI4sEBD/TmC0mYBAWngYQDz74oC5TinXZZZfJbbfdJkOGDJElS5bIs88+K7W1tVJQUKBvf+edd+T000+XNm3aeHTUQG7lLSZ3ppTnpOeHLBIA9zEBgzDxdBcmlUU49dRT4/4pKjhQt6lzPRw4cEAmTZokFRUV+gRz8+bNk1tvvdXLwwYCW97i9+N3AzugAEiEXfYQJp4vok5FBRKzZ8+W6dOny9VXXy3t2rXTOzap/wdM5PfyFr8fv1tMziIB8AYTMAgT4wKI8vLyuMu9evXS54gA/MDv5S1+P34A8AoTMAgT4wIImIdFYdaZvkg66McPAKZPwNCnIgjyIpFIREJkwIAB+ueyZcu8PhTfLgpTDWIYznpJIw94jzYbfuoP1P1STcCEuU9FsJCBQEph3ZWH3TQAAJn2B+nWR1npU5m8gh8QQCClsC4K82PgRKcDAGb3B+n6VCav4BeebuMK84V1UZjfAqdop7OqvFIqtu6TVeW79GV1PQDAjP4gXZ/KVrDwCwIIpF0UpuozY4VhVx6/BU50OgBgfn+Qrk/12+QVwosSJoRuVx4rpT5+286UTgcAnGFnf5CuT/Xb5BXCiwACoTppltX6Ur8FTnQ6AOAMu/uDVH2q3yavEF4EEAiVTBbD+SlwotMBAOe41R/4bfIK4UUAgVAJaqkPnQ4ABIOfJq8QXgQQCJUgl/rQ6QAAADcQQCBUKPUBgPDg/DiAMwggECqU+gBAOHBSNsA5BBAIHUp9ACD47DyDNIB4BBAAACBwgrpphl0o70IuCCAAAEDgBHnTjFxR3oVc5ef8FwAAAAzcNEOdhC0WJ2VLX94FWEEGAgAABA6bZiRHeRdyRQABAAACiU0zEqO8C7mihAkAACBEKO9CrshAAAAAhAjlXcgVAQQAAEDIUN6FXFDCBAAAAMAyAggAAAAAllHCBAAe44ywAAA/IYAAAA9xRlgAgN9QwgQAHuKMsAAAvyGAAAAPcUZYAIDfUMIE+AA18sHFGWEBAH5DAAEYjhr54J8Rdsv2Ktlbdej4da1bNJGhA7t6elwAACRDCRNgOGrkw3FG2D5nlEqXji31T3VZXQ8AgInIQACGo0Y++DgjLADAT8hAAIajRh4AAJiEAALwQY28qomPRY08AADwCiVMgE9q5Bcu3ShVBw9Li5LGeoEtNfIAAMALBBCAD1AjDwAATEEJEwAAAADLCCAAAAAAWEYAAQAAAMAyAggAAAAA/llEvWfPHnnggQfkzTfflMOHD0vfvn1l4sSJ0rlzZ3375MmTZdGiRXG/06FDB3nttdc8OmIAAAAgvDwPIH7xi19IXV2dlJWVSUlJiTz66KNy4403yn//939L06ZNpby8XEaPHi0jRow4/jsFBQWeHjMAAAAQVp6WMO3bt09nE37zm99Ir169dNZh7NixUllZKRs3bpRIJCIVFRXSs2dPadeu3fF/rVu39vKwAQAAgNDyNAPRsmVLeeihh45f3rt3r8ybN09OPPFE6dKli3z22WdSXV0tnTp18vIwAQAAAJhSwhR19913y8KFC6VRo0by+OOPS3FxsWzYsEHfNn/+fHnjjTckPz9fLrnkEhk/frw0b97c60MGAAAAQseYXZh+/vOfy3PPPSeDBg3S6yLWrl2rAwgVNJSWlsoTTzwhd955p7z11lu6zEmtmwAAAAAQ0gyEKllSpk+fLh9++KH88Y9/1P8/fPhwadWqlb6tW7dueg3E0KFDZfXq1dK7d2+PjxoAAAAIF08zEGrNw0svvSRHjx795oDy83UwoRZSq/+PBg9RXbt21T937Njh+vECAAAAYedpALF79265/fbbZfny5cevq6mpkXXr1ukdmSZMmKC3dI2lMg+xGQsAAAAAIQkgVEmSWhSttnFdsWKFXvOg1jlUVVXpwOHyyy/XwcVjjz2md2T629/+JnfddZdeJxE90RwAAACAEK2BePjhh/VWrmpnpf3798t3vvMd+dOf/iQnn3yy/jdz5kx9krnf//73euelwYMHy7hx47w+bAAAACCU8iLqbG0hMmDAAP1z2bJlXh8KACAN2mwAMI8x27gCAAAAMB8BBAAAAAD/rIEAnLJ2025ZsHSj7D94WJqXNJJhA7tJj05tecEBAAByQACBwAYPM+a/L3urDh2/bsv2/TLhZ98miAAAAMgBJUwIJJV5iA0eFHV54dKNnh0TAABAEBBAIJBU2VIiVUmuBwAAgDUEEAik5iWNE17fIsn1AAAAsIYAAoE0bGBXad2iSdx16vLQgV09OyYAAIAgYBE1AknttqQWTKs1D6psSWUeVPDALkwAAAC5IYBAYKlgYeootm0FAACwEyVMAAAAACwjgAAAAABgGQEEAAAAAMtYAwEkOZO1OhmdOp9E85JGMmxgNxZgAwAAEEAAiYOHGfPfjzuT9Zbt+/WuTuziBAAAwo4SJqAelXmIDR4UdVltCQsAABB2BBBAPapsKRF1PgkAAICwI4AA6mle0jjha6JORgcAABB2BBBAPcMGdpXWLZrEXacuqzNZAwAAhB27MAH1qIXSasG0WvOgypZU5kEFDyygBgAAIIAAElLBwtRRbXl1AAAA6qGECQAAAIBlBBAAAAAALCOAAAAAAGAZAQQAAAAAywggAAAAAFhGAAEAAADAMgIIAAAAAJYRQAAAAACwjAACAAAAgGV5kUgkIiFy9tlnS21trZx00kleHwoABJpqZ//4xz/m9DdoswHAvHY7dBmIxo0bS2FhodeHAQCwgDYbAMwTugwEAAAAgOyFLgMBAAAAIHsEEAAAAAAsI4AAAAAAYBkBBAAAAADLCCAAAAAAWEYAAQAAAMAyAggAAAAAlhFAAAAAALCMAAIAAACAZQQQAAAAACwjgAAAAABgGQGEhzZv3iznnnuuPP/888evmzx5spxxxhlx//r37y9BtXPnzgbPV/2Lvibr16+XESNGyDnnnKNfhz/84Q8SZOlej7B9PhYvXixXXnmlnH322fLDH/5QXnnlleO3bd26VW699Vbp06ePXHzxxTJz5kypra2VIEv1ejz++OMJPzuwD202bXZ9tNmZtVNha7cXB7jNLvT6AMKqpqZG7rjjDqmuro67vry8XEaPHq0HzVEFBQUSVB9//LE0btxYli5dKnl5ecevb968uXz55ZcycuRIPUCeOnWqfPDBB/pnSUmJXHvttRK21yNsn48lS5bIpEmT5K677pLvfve78tJLL8ntt98uJ554ovTs2VNuvvlmOe200+TZZ5+Vzz77TN83Pz9fbrvtNgnb66EmItRn46qrrpJf/epXXh9qINFmH0ObHY82Ox7tdnjabAIIj8yaNUuaNWsWd10kEpGKigoZNWqUtGvXTsJgw4YNehBYWlra4Lann35aioqKZNq0aVJYWCidO3eWLVu2SFlZWWADiFSvR5g+H+q5Pvroo3LDDTfI9ddfr68bM2aMrFy5Ut577z3Ztm2bfPHFF7Jw4UJp2bKldOvWTfbs2SMzZszQAVajRo0kTK+H6ozUZ2fo0KGB/2x4hTb7GNrseLTZ36DdDlebTQmTB1asWCELFiyQBx54IO56NYuqMhKdOnWSsFARuAoMElFftH79+ungIer888+XTz/9VHbv3i1hez3C9PlQpSIqSBg8eHDc9XPmzNHpb/XZ6NGjhw4eYj8bBw4c0GVvYXs9jhw5or8XYfhseIE2+xu02fFos79Bux2uNpsAwmVVVVUyYcIEXct+0kknxd2molFl/vz5umxn4MCBevZ9//79ElTqOe/du1dH6BdeeKH89Kc/lTfeeEPftmPHDp3qixWdmd++fbuE7fUI0+dDNb6KCphUqdIFF1wg1113nbz22muh/Gykez1UZkrVEb/66qty+eWXy6WXXqrT4pWVlR4fuf/RZsejzbb+eoSpzVZot8PVZhNAuOzee+/Vqav6UWm0sVE13Gog9MQTT8idd94pb731lowdO1bq6uokaI4ePSqbNm2Sffv2yS9/+UtdmqQWS6sSneXLl8uhQ4calKKo9QHK4cOHJWyvR5g+HyqToEycOFEGDRokTz31lFx00UX6uYbxs5Hu9YgOVJo2barT5tOnT9efJZU+V68Vskeb/Q3a7Hi02Zm1U2Fqtw+EoM1mDYTLq/FV6cWLL76Y8HZVHzd8+HBp1aqVvqzqulVtnKqRW716tfTu3VuCRJUmvfvuu3oRcJMmTfR1anHsxo0bdZpPXafSfLGijUxxcbEETbrXQwUUYfl8qLUvipq5ufrqq/X/n3XWWbJu3TqZO3du6D4b6V4P9dm45JJLpHXr1sd/p2vXrvo6NeOldgFB5miz49FmZ/Z6hKnNVmi3w9Vmk4Fw0XPPPacXeqpUlcpCqH/KPffcI7fccoueXY42NLEfqGjJRhCpHZWiDW/sc1Zb46kSlfrpvOjl9u3bS9hejzB9PqLvr+pwY3Xp0kVvAxi2z0a610OJ7YgUlak64YQTAvfZcBNtdkO02dZfjzC12QrtdrjabAIIFz344IPy8ssv61mt6D9FbTup0ldqbcSNN94Y9ztqliL6oQsaNUuj9oJWMzix1qxZo59v37595f3334/bI/qdd96R008/Xdq0aSNhez3C9PlQC6RVx/zhhx/GXa/Svqeccor+bKiZnGiaOPrZUL9z5plnStCkez0eeeQRXUerdv6IUp2U2go5aJ8NN9Fmx6PNzuz1CFObrdBuh6zNjsBT3bp1izz33HP6/5cuXaovz5o1K7Jly5bI66+/Hunfv3/k9ttvD+S7VFtbG7n22msjV155ZWTFihWRioqKyH333Rfp2bNnpLy8PLJ79+5I3759IxMnToxs3LhRv05nn3125Pnnn4+E8fUI2+fjd7/7XeTcc8+NvPjii/r5/td//VfkzDPPjLzzzjuRQ4cORQYOHBi5+eabI+vXr4/8z//8T6Rfv376tQmqVK/H6tWrIz169IhMmTIlsmnTpsh7770X+dGPfhT5yU9+Eqmrq/P60AOFNps2O4o2uyHa7fC02QQQBnVGyssvv6w/RL169YpcdNFFkQceeEAPloJq165dkTvvvFM/VxUcDBs2TA+eoz788MPI0KFD9SD6+9//fmT+/PmRIEv3eoTt8/HUU0/pIEk1tEOGDNGBQtSnn34aGTlypH6dLr744sjMmTN1hx5kqV6Pt99+W39ezjnnHB1M/frXv4589dVXnh5vENFm02bHos1uiHY7HG12nvqP11kQAAAAAP7AGggAAAAAlhFAAAAAALCMAAIAAACAZQQQAAAAACwjgAAAAABgGQEEAAAAAMsIIAAAAABYRgABAAAAwDICCAAAAACWEUAAAAAAsIwAAgAAAIBlBBCAoZYtWyZnnHGGzJo16/h1n3zyifTq1UvuuusuT48NABCPNhthkheJRCJeHwSAxH71q1/JK6+8IkuWLJFTTz1Vhg0bJlVVVbJ48WIpKSnhZQMAg9BmIywIIACD7du3TwYNGiSdOnWS888/X2cjnnnmGTnnnHO8PjQAQD202QgLAgjAcK+//rrceuutkp+fL2PHjpVf/vKXXh8SACAJ2myEAWsgAMNdeOGFUlpaKnV1dfL973/f68MBAKRAm40wIIAADPfYY4/JV199pcuYJk+eLDU1NV4fEgAgCdpshAEBBGCwjz76SGbPni1jxoyR3/72t7JhwwZ5/PHHvT4sAEACtNkIC9ZAAIY6cuSI/OhHP5KCggJ5/vnnpaioSO677z7505/+JIsWLZLu3bt7fYgAgH+izUaYkIEADDVz5kzZvHmz/Md//IcOHpRx48bp9RATJ07UnRUAwAy02QgTMhAAAAAALCMDAQAAAMAyAggAAAAAlhFAAAAAALCMAAIAAACAZQQQAAAAACwjgAAAAABgGQEEAAAAAMsIIAAAAABYRgABAAAAwDICCAAAAACWEUAAAAAAsIwAAgAAAIBY9f8Dq6vnQ4BrSoIAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import seaborn as sns\n", "import matplotlib.pyplot as plt\n", "sns.set_theme(style='white', context='notebook')\n", "g = sns.FacetGrid(df, col='dataset', col_wrap=2, height=4, sharex=True, sharey=True)\n", "g.map_dataframe(sns.scatterplot, x='x', y='y', s=20, edgecolor=None)\n", "g.set_titles('{col_name}')\n", "for ax in g.axes.flatten():\n", " ax.set_xlabel('x'); ax.set_ylabel('y')\n", "plt.tight_layout()\n", "plt.show()\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 5) Reflection (write here)\n", "**Prompt:** If the per-group mean/variance/correlation were similar, why do the plots look different?\n", "- Which shapes do you observe?\n", "- What does this imply for relying solely on `.describe()` or correlation before plotting?\n" ] } ], "metadata": { "kernelspec": { "display_name": ".venv", "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.11.9" } }, "nbformat": 4, "nbformat_minor": 5 }