{ "cells": [ { "cell_type": "markdown", "id": "f56b531d", "metadata": {}, "source": [ "# Grid-Based Stop Detection" ] }, { "cell_type": "markdown", "id": "9cfdb26e", "metadata": {}, "source": [ "The stop detection algorithms implemented in `nomad` support different combinations of input formats that are common in commercial datasets, detecting default names when possible\n", "- timestamps in `datetime64[ns, tz]` or as unix seconds in integers\n", "- geographic coordinates (`lon`, `lat`) which use the Haversine distance or projected coordinates (`x`, `y`) using meters and euclidean distance.\n", "- Alternatively, if locations are only given through a spatial index like H3 or geohash, there is a **grid_based** clustering algorithm requiring no coordinates. \n", "\n", "The algorithms work with the same call, provided there is at least a pair of coordinates (or a location/spatial index) as well as at least a temporal column." ] }, { "cell_type": "code", "execution_count": 1, "id": "24b50a14", "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib\n", "import matplotlib.pyplot as plt\n", "\n", "\n", "# Imports\n", "import nomad.io.base as loader\n", "import geopandas as gpd\n", "from shapely.geometry import box\n", "from nomad.stop_detection.viz import plot_stops_barcode, plot_time_barcode, plot_hexagons, plot_pings\n", "import nomad.stop_detection.grid_based as GRID_BASED\n", "import nomad.filters as filters \n", "\n", "# Load data\n", "import nomad.data as data_folder\n", "from pathlib import Path\n", "data_dir = Path(data_folder.__file__).parent\n", "city = gpd.read_parquet(data_dir / 'garden-city-buildings-mercator.parquet')\n", "outer_box = box(*city.total_bounds)\n", "\n", "filepath_root = 'gc_data_long/'\n", "tc = {\"user_id\": \"gc_identifier\", \"x\": \"dev_x\", \"y\": \"dev_y\", \"timestamp\": \"unix_ts\"}\n", "\n", "users = ['admiring_brattain']\n", "traj = loader.sample_from_file(filepath_root, format='parquet', users=users, filters=('date','==', '2024-01-01'), traj_cols=tc)\n", "\n", "# Grid-based - data is in Web Mercator (EPSG:3857) projected coordinates\n", "traj['h3_cell'] = filters.to_tessellation(traj, index=\"h3\", res=11, traj_cols=tc, data_crs='EPSG:3857')\n", "stops_gb = GRID_BASED.grid_based(traj, time_thresh=240, complete_output=True, location_id='h3_cell', traj_cols=tc)" ] }, { "cell_type": "code", "execution_count": 2, "id": "4fea8a03", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmkAAAKbCAYAAACuKx3AAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACp70lEQVR4nOzdd3wU950//tfO7GxRWfWGekEUCQSiF9G7sAFjigvGvaQ47XL55XLJfXPO5S6XSy+2E8d23HClmF4FiN4FAoR6QRVp1bVldmd+f6wkJNAMoC3aXd7PxwOklT77mc+MRrsvzXyKQhRFEYQQQgghxK0wQ90AQgghhBByNwpphBBCCCFuiEIaIYQQQogbopBGCCGEEOKGKKQRQgghhLghCmmEEEIIIW6IQhohhBBCiBuikEYIIYQQ4oYopBFCCCGEuCEKaYQQQgghbohCGiGEEEKIG6KQRgghhBDihiikEUIIIYS4IQpphBBCCCFuiEIaIYQQQogbopBGCCGEEOKGKKQRQgghhLghCmmEEEIIIW6IQhohhBBCiBuikEYIIYQQ4oYopBFCCCGEuCEKaYQQQgghbohCGiGEEEKIG6KQRgghhBDihiikEUIIIYS4IQpphBBCCCFuiEIaIYQQQogbopBGCCGEEOKGKKQRQgghhLghCmmEEEIIIW6IQhohhBBCiBuikEYIIYQQ4oYopBFCCCGEuCEKaYQQQgghbkg51A0gxBMUFhbCbDY7rX6VSoXU1FSn1U8I8Qz0WkP68viQtnPnTuj1eqfVHxwcjOzsbOTm5qK8vNzjt+MK3rgvDOPci8499Tv7fHYVbzwHaF/cjyv2x9XvAfRac/96jpk38/iQptfr0dDQ4PTttLW1edV2XMGb9sVVXHU+u4o3nQO0L+7LFfvjbcfM215rvBX1SSOEEEIIcUMU0gghhBBC3BCFNEIIIYQQN0QhjRBCCCHEDVFII4QQQghxQxTSCCGEEELcEIU0QgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtwQhTRCCCGEEDdEIY0QQgghxA1RSCOEEEIIcUMU0gghhBBC3BCFNEIIIYQQN0QhjRBCCCHEDVFII4QQQghxQxTSCCGEEELcEIU0QgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtyQcqgbQPqLiooa6ibYLTg4eKib4LG85dhFREQMdRMczhW/mzzPo7293Wn19/xcdDodwsPDnbYdV/HG88xVvOG15mH4+VNIczOZmZlD3QSHYVl2qJvgMXqOVXZ29hC3xLEYxnsu1nvT72ZWVhaysrKGuhkO403nmbN562uNt6KQ5mZyc3NRXl4+1M2wS3BwMLKzs2G1Woe6KR6j51jl5uaira1tiFtjP51Oh6ysLAiCMNRNcRhn/272/N5cuHABtbW1TtsOwzBe83PxxvPM2bzptabn5+/NKKS5mba2NjQ0NAx1M8gQKS8v94qff3h4uNe9eLrqd7O2thYFBQVO34438MbzzFW84bXmYfj50zViQgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtwQhTRCCCGEEDdEIY0QQgghxA1RSCOEEEIIcUMU0gghhBBC3BCFNEIIIYQQN0QhjRBCCCHEDVFII4QQQghxQxTSCCGEEELcEIU0QgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtwQhTRCCCGEEDdEIY0QQgghxA1RSCOEEEIIcUMU0gghhBBC3BCFNEIIIYQQN0QhjRBCCCHEDSmHugH2Cg4Odmr9ERERAACdTofw8HCP344reOO+uIqzz2dX8cZzgF4D3I8rjpmrf/6u4g2vNa4+ZkNBIYqiONSNIMTdFRQUwGKxOK1+lUqF1NRUp9VPCPEMzn6t8fX1RWJiotPqJ47l0VfSCgsLwTDOv2OrUCjgiizrqu24gjfti8ViceqLJgCYzWaXnc+u4k3nAL0GuC9XHDNX/Vxc8VpjtVoBALm5uWhra3PqtpxNp9MhKytrqJvhVB4d0sxm81A3gRCHofOZEOIq5eXlaGhoGOpm2CU8PNzrQ5r3/NlOCCGEEOJFKKQRQgghhLghCmmEEEIIIW6IQhohhBBCiBuikEYIIYQQ4oYopBFCCCGEuCEKaYQQQgghbohCGiGEEEKIG6KQRgghhBDihiikEUIIIYS4IQpphBBCCCFuiEIaIYQQQogbopBGCCGEEOKGKKQRQgghhLghCmmEEEIIIW6IQhohhBBCiBuikEYIIYQQ4oYopBFCCCGEuCEKaYQQQgghbohCGiGEEEKIG6KQRgghhBDihiikEUIIIYS4IeVQN8AeKpUKDOM9OVMQBJjNZqduw1XHjGVZWK1Wp2/HFbxpXwDXnGeAd/1+etM54E37AnjX/rjqd5N4Do8OaampqUPdBIcrLCx02i+pSqXyymNGHpwzzzOAzjVCBsvZv5vEs3h0SAOAnTt3Qq/XD3Uz7BYcHIzs7GynXnnoqdvZx6xnX3Jzc1FeXu607biCN+0L4JrzDHDdueYK3nQOeNO+AN61P6763SSexeNDml6vR0NDw1A3w6O46pi1tbV5zc/Gm/bFlbzp99ObzgFv2hfA+/aHkB4U2QkhhBBC3BCFNEIIIYQQN0QhjRBCCCHEDVFII4QQQghxQxTSCCGEEELcEIU0QgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtwQhTRCCCGEEDdEIY0QQgghxA1RSCOEEEIIcUMU0gghhBBC3BCFNEIIIYQQN0QhjRBCCCHEDVFII4QQQghxQxTSCCGEEELcEIU0QgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtyQcqgbQLxXVFTUUDfBbsHBwUPdBI+WkJDg8cfQ09s/EJ1Oh/Dw8KFuht0iIiIAeMf+9OyLq3jDee3qYzYUKKQRp8nMzBzqJjgMy7JD3QSP0nO8srKyhrgljuNN50BWVpZX/Wy8aX8UCoVT6+85j7Ozs526HeIYFNKI0+Tm5qK8vHyom2GX4OBgZGdnw2q1DnVTPErP8dq5cyf0ev0Qt8Y+3ngOXLhwAbW1tU6rX6fTISsry+nbAQCGYSAIglO34Qo9x0wURadup+c8zs3NRVtbm1O35Ww9x8ybUUgjTtPW1oaGhoahbgYZQnq9ns4BN1RbW4uCggKn1R8eHo6srCynb8eb9BwzVykvL/f4301XH7OhQAMHCCGEEELcEIU0QgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtwQhTRCCCGEEDdEIY0QQgghxA1RSCOEEEIIcUMU0gghhBBC3BCFNEIIIYQQN0QhjRBCCCHEDVFII4QQQghxQxTSCCGEEELcEIU0QgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtwQhTRCCCGEEDdEIY0QQgghxA1RSCOEEEIIcUMU0gghhBBC3BCFNEIIIYQQN6Qc6gbYKyEhAcHBwUPdDLu5ch9SUlKcuj2dTgcAiIqKcto2evA8j/b2dqfVHxERAcC2T+Hh4U7bjqv07I+reMPvpjedA67aF286Zq7i6t9N4hkUoiiKQ90IcltJSQkMBoNT6tZqtUhOTnZK3cSzFBQUwGKxOK1+lUqF1NRUp9VPiLdy5nsAAGg0GqSkpODDDz9EQ0OD07bjCuHh4diwYcNQN8OpPP5K2s6dO6HX64e6GXYLDg5GdnY2nJmZe+r2tmNWVVUFk8nktO0oFAqn/lxczWKxODWgAYDZbEZhYSEYxjt6VHjTOeCqfXHFdtRqNWJjY5Gbm4u2tjanbsvZdDodsrKyvOY8I47h8SFNr9d7/F8DruZtx8xkMsFoNA51M8gdzGbzUDeBPCTKy8s9/jUtPDwcWVlZQ90M4ma8489cQgghhBAvQyGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtwQhTRCCCGEEDdEIY0QQgghxA1RSCOEEEIIcUMU0gghhBBC3BCFNEIIIYQQN0QhjRBCCCHEDVFII4QQQghxQxTSCCGEEELcEIU0QgghhBA3RCGNEEIIIcQNUUgjhBBCCHFDFNIIIYQQQtwQhTRCCCGEEDdEIY0QQgghxA1RSCOEEEIIcUMU0gghhBBC3JByqBtACCGE3IsgCDCbzTCbzeB5HmazGQzDwGQyoby8HA0NDQgMDERgYOBQN5UQh6GQRggh5L5ZrdZ+YaknMPUNT32/NmA5Cw+zyfax9+u8GRaeh9nMw8zbyll4HiaT7flWqwUAwKlUUHX/U6tV8PP3h9FoAMMwqDxShfXr1kOj0QzxUSLEMTw+pCUkJCA4OHiom2E3b9gHb6VSqcAwzu8ZwLIsrFar12yn58qHM9HPxsZsNuPWrVt3haW+j00m011fN/O2f7zJFoz6lePNd9dnMkMQbO1TqdW9YUml4qBS3w5OnEoFtVoNJaeEWq2CSm17zPlx8FUFQKVWQ61SgVNxUGvU4Djb89UqW1mOU93x+HY5hULRu98Mw0CtVqO1tRVWiwW//Z/fIz8/H8uWLXPYz8RVIiIiXLo9b3jPcfUxGwoeH9KysrKGugkO1fcFiAw9lUqF1NTUoW6GxyosLHRaUKOfjY0gCNj47EZcL7iOyMgIW3hSd4cljrOFpJ6wo1JBrVFDreOgU/v1Bque76t7gxcHtUbTG77UKvXtst313s9rlclkhCCITtv3niYEBAQAAL7//30PLz39KiZMmIDRo0c7bbvO5Oz3AJZlAQDZ2dlO3Q5xDI8PaTt37oRerx/qZtgtODgY2dnZEEXnvaCRB9dzlcbZ51nPzz83Nxfl5eVesx1nXuWin41NXl4eGhrr8dFX7zv5qqIIs9kEs9l0z5KsUtkdnBQAnPea1vNy2dHRAavVClbJYvX6VfjOd7+D17/9ukuusjqKTqdDVlaW098Deq7U5ubmoq2tzWnbiYqKQmZmplO303PMvJnHhzS9Xo+Ghoahbgbxcq46z9ra2rxqO67wMP9srFYrDuUcxPd+/B1YBQFWQXBy69yT2WQCb7H1WZuzcDZ2bN2Jr776CmPGjBnilt2/8PBwlwaOnsEWnrwdVx+zoeA5f2YQQgjp58qVK4iJj8H4ieOHuilug2VZvPbdV3Hu/Fl0dXUNdXMIsQuFNEII8UBmsxl5ly9h44sbqC/rHVJHDsfsBbNw8uTJoW4KIXahkEYIIR7o0qVLSM9Ix8i0kUPdFLf09HNPoba+FjU1NUPdFEIGjUIaIYR4GIPBgKvX8vH0808OdVPcli5Ah+de2Yhjx3MhPKR99Yjno5BGCCEe5vz585g+axriE+OHuilube7COQiLCMOlS5eGuimEDAqFNEII8SDt7e0oLCrE+mfWDXVT3B7DMHj1O68g73IeOjo6hro5hDwwCmmEEOJBzp49i4XLFiAi0vtnW3eExOQELFy2AMdPHB/qphDywCikEUKIh9Dr9aisqsDapx4f6qZ4lCc3rkdTcxMqKyuHuimEPBAKaYQQ4iHOnDmNFWseRUBgwFA3xaP4+PrgpW+8gOMnjsPSPektIZ6AQhohhHiAuro6NOobsWL1o0PdFI80c84MxCbE4OLFi0PdFELuG4U0QgjxAGfOnsbap9fCx9dnqJvikRQKBV55/WXkX73i1DUrCXEkCmmEEOLmKisrYTQZsSR70VA3xaPFxsUge1U2jh07NtRNIeS+ePwC64QQ4s1EUcTpM6fx9PNPQaVWDXVz7kkQBBgNRrS1tqG9vR1GgxFmsxlmMw+TyQTezMNsNsFsMsNkNnc/Nt9+bDLDZDJ1P8f2NQvPY80Ta5A2ZrTd7Vvz5GocPnAYpaWlSEpKcsAeE+I8FNIIIcSNFRcXQ61VYfb8WQ6pr7S4DEajEWaTLQTxPA+zqftjdzAym3jwvLk3XJm7g5LJZOoud7ssz/PgTbcfW61WAIBGo4FSpQTHcVCpOHBqle2jioNKrbJ9VHV/VHNQcraPKhUHH38NdJw/1BoVOI7DrfpG/OZ/fou33vsrVCr7gqpWq8Ur334Jf/r1XxEXFwelkt4Gifuis5MQQtyUIAg4e+4svvmDV8GyrN31ffT+J/h689dIHJ4AlVoFJaeEWq3uDVNqjS00cSoOKn8OPuoAcJwSqt6AZfu+Wt2nnEp1R8jq+Z4KoiACEO1utyiKKLpajC8++RJPPWv/UliTp01G6uj9OHfuHKZOnWp3fYQ4C4U0QghxU9euXUN4ZBgmT5tsd1211bXYv2c//rTpdwgODXJA66QpFAqH1/fstzfg37/x/zBv4VxERUfZXd9L33wR3335exgxYgSCgpx7PAgZLBo4QAghbshiseDipQvY+NIGh4Sezz7+HAsenef0gOYsUTGRWLxqIf76x7cgivZfnYsaFonVTzyGY8dyHVIfIc5AIY0QQtxQXl4eUkakYMy4MXbXVVpchrxLl7F87VIHtGzorHzyEVRWVuDksVOOqe/xFeAFC4qLix1SHyGORiGNEELcjNFoxOX8y3jmxacdUt+mDzdh+fpl8PP3c0h9Q0WlVuGZbz2Nt/78NowGo0Pqe+27r+Dk6ZMwm80OaCEhjkUhjRBC3MzFixcxYXImkocn211X/uWrqKioxOIVCxzQsqE3Ydp4xKXE4pMPPnVIfeMnjENG5hicOXPGIfUR4kgU0gghxI20tbXhesF1PP2c/aMYRVHEJx9swmPPrIBao3ZA69zDxm8+jV07dqGi3DELpr/w6vMoLi1GY2OjQ+ojxFE8fnRnQkICgoODh7oZdvOGffBmzj7PeurW6XQIDw932nYiIiJcuh1XcPbvjquP2fHjx7Fg6TzEJ8TZXefpE2fQ0dmOOUtmOXzEpSwnbyssMhTL1y3DX373F/z+zd/ZXV9oeCjWb1iLPV/vw4pHV7j2WBEiw+NDWlZW1lA3waHoxcG99MxN5arzLCsryyXbctV2LBaL0+oWBAEAkJ2d7bRt9OWKY1ZdXY28y5fwt4/eQnBIiF11CYKAzz75AuteWAO1RuOgFj4gBRwxTdqAstcsRe6+Yzi0PwfjJ2TYX9/KZTiw5yAKCgowatQoB7SQEPt5fEjbuXMn9Hr9UDfDbsHBwcjOzqah4G6mZ/Z0Z59nPT//qqoqmEwmp20HsP0h4IrzzGKxODWkmc1mFBYWgmGc32vDVcfsd7//HRZmL4QuQGf3ebB/9wFACUyZNbl7UlkACtu+CFYBogPSk9FgxOE9R1BZWoW4pFjMWTIbGq0GCijAsIzTAhoAcJwSz77+DP76qzfx+7/+Biq1fbdzOY7Dq995Bf/1779EYmIiNEMVbAnpw+NDml6vR0NDw1A3g3g5V51nJpMJRqP9o9YeFt40Iq+8vBwXLl7AP374t96rhINlNpvxwXsf4vnvPQOG6XM5S1QACkAUBYeEzpxdh7H/64NQa9UouHwDoiBiyapF3bc7nR+ex2SmY8TY4dj04WfY+OIzdteXPjYNk6dPxunTpzF79mwHtJAQ+9DAAUIIcQOfff4ZVq9fBX+dv9117dy2C6GRIRg70f451uRUlVVBrVUjLjEWaq0aVWU3nbq9gTz92lM4fOgIykrKHVLfsy8/g4qqCtTX1zukPkLsQSGNEEKG2PXr11FRWY4Vqx+1u66uzi589vHnWPfiGqf3cY1NjIXJYEJlWRVMBhNiE2Ocur2BBIcGYdWGR/HOW/+w+wokAAQFB+Hp559E7rFch9RHiD0opBFCyBASRRGffvYpnti4Hhqt/f2gtnyxFcPTU5A6OsUBrZM3Z8ksLHx0PhJTErDw0fmYs2SW07c5kMUrF8LEG5Fz4LBj6steBD+dH65eveqQ+ggZLApphBAyhC5cuIC29lYszl5kd10tzS3Y8tVWrHt+tQNadm8arQZLVi3CS99/HktWLXJIyBwMlmXx7Osb8MkHm9De1u6Q+l777is4d+Ecurq6HNBCQgaHQhohhAwRQRDw6Wef4pkXN4DjOLvr++zjLzBxRiZiElx/23GojUhLxfhp4/DJB5scU9+oVMyal4WTJ086pD5CBoNCGiGEDJFjx46BUysxa57986/V1dZj7+69WP3MKge0zDM98eIanDpxGoUFRQ6pb8PzT6G2rgY1NTUOqY+QB0UhjRBChgDP8/jyqy/x7MsbHTLX28f//ARzls5CWGSoA1rnmXSBOqx9YbXDBhHoAnTY+MozOHb8GA0iIEOCQhohhAyBgwcPIiIqDBOnTLC7rvKyCpw4fgIrn3zEAS3zbPOWzQGjVGDvrn0OqW/+onkICQtGXl6eQ+oj5EFQSCOEEBczGAzYsm0Lnn1lo0OmyfjgHx9g6WOLoQvUOaB1no1hGDz7+gZ8/skXaGlucUh9r333VVzKu4jOzk77G0jIA6CQRgghLrZ7926MSh+J0emj7a7r+tXruHr1GrLXLHVAy7xD8ogkTJ07BR++97FD6ktKScSCpQtw/Phxh9RHyP2ikEYIIS7U1taGnbt3YuOLG+yuSxRFvPf3f2LFk8uh9aG1Jvta+9xjyLuUh6tXrjmkvic3rkejvhGVlZUOqY+Q+0EhjRBCXGjb19swdcYUJCQl2F3X+TMXUFtbi4WPzLe/YV7Gz98PT7y0Fv94613wPG93fb5+vnjxm8/j+InjsFqtDmghIfdGIY0QQlyksbEROTk5ePr5p+yuSxAEvPfO+3j82VXgVPbPseaNshbOgG+AD3Z9vdsx9c2ZiZj4aFy4cMEh9RFyLxTSCCHERb766kssXDYfkVERdtd17PAxmC1mZC2c4YCWeSeFQoFnX9+ALV9sReOtJofU9+p3XkH+1Stoa2tzQAsJkUchjRBCXKC6uhpnzp3F+g3r7K7LYrHgn+9+iHXPr3bIHGveLC4xFrOXZuGf7/zTIfXFxsUge1U2DSIgLkG/3YQQ4gKfff4ZVqx5FIFBgXbXtW/XfvgG+GDC9Ez7G/YQWP3MKhQWFuHi+UsOqW/Nk6vR3tmGsrIyh9RHiBQKaYQQ4mQlJSW4UViAx9astLsuo9GEjz/4BOtfXOuQOdYeBlofDZ5+bT3effs98Gb7BxFotVq8/O2XcOLkcVgsFge0kJCBUUgjhBAn2/TpJqzbsBa+fr5217V983bEJsdgdMZIB7Ts4TFl1mSERoVgy5dbHVPf9MlIGZmCc+fOOaQ+QgZCIY0QQpzoypUrqGuoQ/aKZXbX1d7Wji8+/RLrnn/cAS17uCgUCjz3+jPYuW0X6mrrHVLfy998EdcLrqGlpcX+BhIyAApphBDiJKIoYtOnm7DhhaegUqnsru/LTzcjfeJoJA5PsL9xD6GomEgsWrUA7/3tPYiiaH990VFYtW4lco/lOqQ+Qu5EIY0QQpzk9OnT4K1mzF80z+66mhqbsOPrHVj73GoHtOzhteKJ5bh58ybOnjzrkPoeW7sKvMWM4uJih9RHSF8U0gghxAmsVis++/wzbHzpGbAsa3d9n3zwKabPm4LI6EgHtO7hpdaoseGbT+G9d/4Jo9Fod30qtQqvfudlnDx9Emaz2QEtJOQ25VA3gLheQkICgoODh7oZdnPlPjj7mLlyX1QqFc2t5QL79++HLkiHmbNnPMAoTBGCcPdts+qb1Th86DB+/e5/O6ZxCgWcOS7U5aNOH3B/Jk7PxOHdR7Hl86145vmn7d78lGmTMW7COFy/fh0LFy4cVB0REfZPcPwgnP2a01O/M7fj6mM2FCikPUR6/prPysoa4pY4ljPfEFx9zJz95qZSqZCamurUbRDAZDJh29fb8IOffBcazYMtfG4yGe8Kah+++zEWPDoPwaFB9jWs+/RyxJU99yACUAxqf557fSP+7ZWfYtWalYhPiLe7JT/8yQ/w2rPfxIwZM5CUlDToepw9pUfPscrOznbaNkRR7F0v1ZnbeRhQSHuI9CwKvHPnTuj1+iFujf2Cg4ORnZ3t1A67rjpmrtgXAL1X0LzlHHBXx48fR0CIDvGJcdA33d9yRKxSiYCAANiS1O3zoKSoBOfPncfvPvxf+xtmyzQQrAJEOO9cU0ABhnXN1VpREIBBXBcMiwjFsjVL8Ltf/R7/87v/tvsPpKDgQKxauwL/9pN/wxv/+cag6rNYLA4NaYIgwGKxwGw29/5jGAYGgwGnTp1CS0sLeJ4Hz/OwWCy9n/c8NplMMFt4mHkzzLwJJt4Mnrd0P+bB82bwFgt4c/dHnoel+7kAsGHNU079A1en03ndRYc7UUh7COn1ejQ0NAx1MzyKtx0zb9sfd2I0GnHsxDH88rdvgHfAG+7773yA5euXwc/fzwGtsxFFwbl/ECgUcG2X58Hty/I1S/GjfceQs/8w5iyYbXcrVq5Zgf27D2D//v2YNWtW79cFQegNPyaTqffznuB05+d9y5nNZpjMZhjNRpjMptv/TCaYeR5Gk8kWmnqeb7YFKEt3PSIABaMAyynBciyUKiXUWg06DJ0QFCIYJQOFkoGC6/54x+eMkgWjZsD4sWBYBgzHQqFUglGqwSpZcD1llLbvMUoGDMvCoO/ARx9ugkajgb+/v93HdiDh4eEU0gghhNy/ixcvYuKUCUgenmx3XZcvXUFJSTG+9R+vOKBl5E6cisOzrz+Dt/7n75g0daLdkw1zHIdVa1fiZz/9T4S9H2YLXGbblSURIhilEqySAcuxYJRKMBwDRZ+Ao1AqAJYBwzGAUgEoGShYBRRcdxDiWLD+bG8YYjkOjFINNcdCy7FglWxvOVsZ29cUzO2reoyCgVqjRlNTEyy8826t+kUFQjcuAjm5h/Hoskecth1vRyGNEEIcpLOzE9cLruP3b//G7rpEUcR7f38fjz2zEmqN2gGtIwNJGzcKgiigrLQc6WPT7K7vg/c/gmZ4IBIen2ALVT2hie0flh4G0bNTce2vx1BRUYH4ePv7/T2MaIgXIYQ4yLlz5zB30WwMixlmd12nT5xGS2sz5iyZde/CZNAO7zmKyIhIpI0ZbXdd58+ex+XrVzFmzTT4hPhDo/MB56MGyykfuoAGAEo1h6iFqdh7ZD+tcTpIFNIIIcQBWltbUVpWirVPrbG7LkEQ8P47H2DNc6uhVNIND2cxm8zY/MFWPPfysw4ZWf3zn/4C0bNT4RPsuP6Dni40bRj4QAXOnac1TgeDQhohhDjAuXPnkL1yKUJCQ+yuS9+kR/XNaqSNG+WAlhEpe7buQ0JCIjLGj7W7roN7D6G0ugIjloyzv2FeRKFQIHbpaJy4eAqtra1D3RyPQyGNEEIcgGEZdHR0OKSu0LBQLFg0H5ve+cIh9ZG7dbZ3YvumXXj2xY0Oqe+X//W/iFuUBrXvg82L9zDwCfVH4KRoHDx6aKib4nEopBFCiANMnTIVRw7m4sb1QofU9+zLz+LSyUu4ke+Y+tyV0WDE7s178fb/vYPdm/fCaLB/qab7sf2LXcickInhI1LsruvLTV/hVkcTUueOcUDL3JMoihAsVlgMZpjbjTA2d6LrVjs6alvQVqWHuUP+5xY9czgqmqpRUlLiohZ7B+rsQAghDqDVajExcyLe/MPb+M1f/tfuWf0DAnR49sVn8e4fPsAv3/q5F60S0F/O7iPYt+0A1Fo1Cq7cAAAsfWyxU7fZ3NSC/VsP4o9v/d7uugRBwG9/90ckLh0DpYazv3H30BOWBN4Kq8UKwSJA4Ps+tkIw9/mcv6Os2QoFgPbWdljMPKx899ct1u7PrRAtAqy8FaJVuF2HVQBgm3NNwfaZU41lAAZgBQUmvDYPLDdwrGBVSgxbMgL7dh/Ai3Fx4DjnHytvQCGNEEIcJC0tDTcKb2Dfrv1Y+sgSu+tbtGwh9u3ehz1b9iP7cfvrc0eVpVVQa9WIS4xFZVkVKstuOn2bWz7ehjlzZyM6Jtruuv7x1nswMGYkZcmPDm2r1qO57NbtgMVbYeEt3Z9bYO35aBEgmC22UNUTrrrDU0+YAgAwitsT0XaHpt451lhbkEL3vGsKVmH7mtL2PUbJQlQD8FFAwaq6A5cCyjvDl1IBhZK1zdXWs50BRqmKoojmnUW4eawI8XOl+1EGp0ai6cJNnD5zGjNnzLTnsD80KKQRQoiDMAyDrJlZ+PAfH2PqjCkICrZvrU2GYfDN734DP/rBjzFtzhT71+50Q3FJsSi4cgOVZVUwGUyIS4xx6vbqaupx/MBJvP3em4Ouw2gwYt/u/SguKsG777yP4U9NAquUvtLZWtmEY3/cBZ/RYQDbZ5LankCk7glGmt6vcawCCo7tX45lZMPS/VIwDKwWi8NWnVAoFNBlxeHmF9cRPjYW2pCBR7cqFArELRmN0387idGjRjt9kXdvQCGNEEIcKCIiAvGx8Xj/7x/gez/6jt31JQ9PxvwF8/Hx25vw7Z98wwEtdC9zl9qWY6osu4m4xJjex87yxfubkb18GULDQgddx77d+7Fjy06UlpTByJuhFKUDmiiKKNh5Af4zY6Ebb//8ee5KGaCBZmw4indfRvpT0ySnNNEE+SJkWhz2Hz6AtavWOGTqE29GAwcIIcTBpkyZgjPHzyD/8lWH1PfM80/j2sUCXLmQ75D63IlGq8HSxxbjlR+8gKWPLYZG67zRkeXFFbhyJh9rnnzcrnrKSsqg5JSovFkFXWwwOhqkp5bQF9Wh5VYz/MdE2rVNT+CfGYX2xlY0Xa+VLTdsejJqO2+hsNC7B8U4AoU0QghxMI1Gg8mTpuCtP7wNnuftrs/Xzxcvv/YS3v/jh+DtWG/RNpJyD/72m39gz5Z9LhtJ6S4+e/dLrF77GPx19i34nZiciLyLlwEVA7WPGv6RgQOWEwUR13ddhO+kYbb+YF5OoWTgPysOJXuvwGqSPk8ZJYvopaNw4NghmM1mF7bQ83j/WUMIIUNg5MiRsJit2LVtt0Pqm7NgNkJDwrDzi8HXl7P7CPZu2Y+y4nLs//ogDu856pC2eYJreQWoKKrEitWP2l1X5qTxaGxuQlRmAmKnp2LYpKQBy9XnV6HLYIDvyDC7t+kpNPGBQJgGlUcLZMsFJYeDjfXB8VMnXNMwD0UhjRBCnEChUGDmjJnY9MFnaGxodEh93/jOa9jx2S7cqhtcfZWlVdBoNYhLjIVaq0aVC0ZSugNRFPHZP77Ak8884ZDbqb/6xa8RlhmHSS/PR9yMVChVd3fvFqwCbuy+CL8pwx66dTt1M2NRfaEcXQ1tsuViF4/GhWuX0Nho/++Ht6KQRgghThIaGoqUpBT84633HFJffEIclmUvxQd//XhQz49LioXRYOwdSRnr5JGU7uL8yYtob+7Akmz751+rKKvAoSNHMOLRTNly1WdLYGYFaJPcdwSjKIoQLQIEowXWTjMsrUbwTV0wN3TAVNMOY1UrDGXN6CpqQmfBLXTk16M9rw6d1xpkR4Yq/dXwyYxE0a482XJqnRZhWQnYe2gfBEFwxi56PBrdSQghTjR58mR8+vmnuHj+EsZPGGd3fU88sx4vP/sqLpy6hMypD1bf3KWzoWBsHehjE2IwZ8ksu9vj7gRBwOf/+BIbn9/gkMXqf/ZvP0fYhDgExUiPDrWaLSjcfxn+8+Lua/SiKIgQrQJEqwjRYoVosYUn0Sr0/2iReGwRIPR+tPaWEXiJ5/T5HADA2OZRA6sAmNvzqvX9J/aUYRTALSNCGYXsbVy/jAg03riGW1duInxsrGS5yClJuJ53HNevX0daWto9j9XDhkIaIcRj8TyP/Px8NDY2IjQ0FOnp6W43k7lKpcK0KdPw1h/exp/f+SM4lX3t0/po8eq3XsHf/vQ3pI8fDZVadd/PtY2kXOLQObLc3fGDJ8ApVciam2V3XVfy8nH20gXM+o9VsuUqTxQCgSpoYgIky4iCiLotV2GoaOn5CqDoP/ns7X/M7YDE3g5L4p1fYxWAhumugx3guX0mtmVgm7CWVTzw7Vihvgv6o+XQJgaBUQ8cIxQsA11WHEr25yN4eASU2oHPU4ZlELNsNA59cRjJycnQaGjt074opBFCPFZ+fj7y8vLAcRyqq6sBAOPHjx/iVt0tJSUFBQXXseWLrVj71Bq765ueNQ17du7Ftk07sObZxxzQQu/Em3l8+f4WfPu73wbD2N+75//9+39i2Ixk+IXqpLfZZULRoXwEPiK/JmjnjUYYmjrBPJoAdK8E4NI5w+wI6UyED6wharScrELwnETJcuoYHbqifVGeU4CUZWMlywXEh6ApJQC5J45h4bwFg26XN6I+aYQQj9XY2AiO4xAaGgqO49y2A7JCocDMmVn46tMtqKutd0h9r337FezZvA+1N+sc0ELvdHBnDiIiIjFhsnz/sftx/OgJFJQVY8RS+T8Cyo5chzLaD6rwgWfdBwDRKkB/vByKjBAoNErbKgIeNqkrMy4U7dcaYG7olC2nmx6LuvxKdNS2yJaLWTASl4uvoq6Ozue+KKQRQjxWaGgoeJ5HY2MjeJ5HaOjgZ5F3tqCgIIwaMQp/+/PfHXKrcVjMMDz2+Cr8888fPjS3Lh+EocuIrR9vx3MvbXRIAHrj579E7NwR0Ab4SpYxtRlQdrwA/lPk1wRtv1IPq0oBxTDputydQqsERgWi6VAJREH6/GN9VfCdGIWinZdky6n8NAifk4R9h/fTIII+KKQRQjxWeno6MjIyEB4ejoyMDKSnpw91k2RNnDgRhdeKcObkGYfU9/gTq9Fw8xZO556VLWebxHYv3v6/d7B7856HYhLb3Zv3YtSokRiVJr3g9/3auW0XqhvrkLooQ7Zc8cErUCUHgQvSSpYRzFY0n6wExgR73NWzOzHDA2E28ui82iBbzndMBIxmM+ovVsiWi5yQgBZFF67kX3FkMz2ax/dJS0hI8IpFWl25D95wvADbGomu4uxj5sp98SYcx7llHzQpSqUS06dNx9/+9A4yMjPs7iStVqvxje98A7/7ze+RMXEstD4D15ez+wj2bTsAtVaDgisFEAQBS1YusmvbclwfPvpvr621Hbu+3IPf/uk3dvdFEwQB//ur3yBhaTo0PtLhq7OxHVXnShH6pPwfCm0Xa2AN5MCG+djVLnegYBRQZIai+XgFtMnBYH0GHhSjYBTwnx2H0l3XEDIyCpyvWrJcbPZo5Hx8FMNThsPHx/OPkb08PqRlZdk/YsedOPPFrecScnZ2ttO2MRQslsEvk3Mvrj5mztwX4h6SkpJQUFCALz75Chuef8ru+iZOmYBRI0di80db8dTL6wcsU1laBbVWjbikWFSVVuFmeQ1YB0xH4Q4UA4Sw7Z/uwMysmRgxMtXu+t//+wdotXRi+uKJA05a2+PygRPwSQ8Dp5MO3lYDj5ZzN6GY5T0LrTOhWlgjtWg5XoGQhdKDJdSR/uASdCg7eA2pj0r/YeUfHQTf0aE4cuwoli5a4owmexSP/y3duXMn9Hr9UDfDbsHBwcjOznZq3xKz2YzCwkKHjHJyFxaLxanBxpXHzNn7QtzHzJkz8dXmLzFnwWzExtk/oezL33wZ33jxm5i9aCZiEu6uLy4pFgVXbqCytAomgxExCdGwWq12b1eKAgowrCteZ2x/1Fqt1t7Rik0NTTi4/RB+++ffQN/UZFftFosF//eb3yFm2Ui0tksvot5R14qqy2UIe3oMrDK/w80nKyBEaMEGqgYeXemhtz+ZjFB07KmEX3oE1FHS66L6T49Bw8f5iBwfD12s9N2JmHkjcO2vuRhbPQbR0fL9+7ydx4c0vV6Phgb5++HkNlrM9sHRMSOOptPpkJ42Bm//8W9449c/t/sKenhEGNY9sRbv/uED/PS3P76rvrlLZwMAKstuIj45FrMWzoTozM7ZCgVc2uVZFHv/wP3yg62YPW82gkODwdv5R8/f//IOzBoBkePiYZFZ2L5472VoMsKhULOSf2hb2k1ou1wHZqH3rfKgULNAWhCaDpUi6omxkvOusRoOflOjUbQrD5kvzR7wKigAcD4qRC4Yjr2H9+PZJ57xqgsLD+rh3XNCCBlCmZmZqCqvwrHDxx1S38o1K9DZ2oVjB+6uzzaJ7WK88oMXsfSxJfe1fqXRYMSeLfvw99/+A3u27LNrsEH/gQt7nTZwoaayBmePnsNja+Unm70fRqMRH378CRKWpsuGhNaKJrTU6OGXId+vtPX0TYhxvlD4uddky47CJAeAt1rQkSc/hYbPqDCYYUHNmTLZcuFjY9Gl5XHp0iUHttLzUEgjhJAhwLIsZkyfiXf++i66Orvsro/jOHzjO6/h47c/Q2e7/NxV92I0GPGHN/6Mv/3mHziecwp7tuzD4T1HB11fz8CFsuJy7Nt2ADm7j9jVPimfv/cVlj6yFAGB0jP9368//eYvYMO0CB8tfbtNFEWUHbwK34lRYDhWshyvN6Djxi0wad4xaGsgCoUCiswwNJ+qhLVT+u6DglFANzseFUcLYG6XDusKRoHYpaNx5OwxtLe3O6PJHoFCGiGEDJG4uDgEBwVj0wefOqS+jPFjkZk5Hp+995Vd9RzecxSnj56F0WhCq74NLc2tqCq7Oej6egcuJMZCrVWj0o66pJTcKEXB5UI8ssr+QT5trW3Y8vXXSFiaDoVC+m2yuageXW2d8B0tvYYlALScrARSAqDQeHwPI1lMsAZijC+aj5bLllOF+UI1PAil+/Jly/lGBsA/IwKHjzkn1HsCCmmEEDKEZkyfgX07D6C0RP72z/168bUXcOrQaZQVlQ+6jqqyKgQE6qBWq2AymdDW3IbYxMH3pYpLioXJYEJlWRVMBhPi7KhLymf/+BKr1qx0yLQN//tf/wefhCCEJEdKlhEFEaUHr8J3yjDbGpgSzPUd6KpsATMi0O52eQJmTAg6K5phrJIeaAEA/lOi0VRWj5ayW7LlYmanoqimFJWVlY5spsegkEYIIUPIz88PGWMz8Nffvz2omdaNBiO+3rwdf/j1H/H15u3w8fHBUxufxLt/+OegZ26PTYxFQHAAAoMCoNGoMWXWJMxZMmtQdQG2gQuLVixA4vBELFqxoHcgg6NcuXAVdVX1WLzU/rnf6mvrcSAnB4lL5ec7u5V/E7xohTYlRLZc8/EKYEQgFCrp26HeRKFigbEh0OeUQrRKn3+MWgnf6TEo3pUHwSJdTqnhELUoFXsO73soR79797VXQgjxAOPGjcPX27dhz469mL943gM9d9/u/dixZSc0Wg2uXr4GAFi+Mhv79uxHzu4jmJ8994Hb0xPIqspuIjYxBnOWzLqvwQZSegYuOIMoivjsnc+x9ok14FSD65RvNBpxaF8OKssrcejAYQSMDEdAtHT4EixWlOVch19WtORIRgAwVrXC2NgJZmLcoNrlqZh4f/BlbWi/WAvdROk+fT6pIdBfa0TNqWLEzJSe0y40LRr6i9U4f+E8pkye4owmuy26kkYIIUOMYRgsWbwU7739Ptpa2x7ouWUlZdBoNUhISoBGq0F5aTlYlsW3vvsNfPbOl2hrffBO1xqtBktWLcJL338eS1YtsiugOdvpo2dgNlgwa97gJzY/tC8He3fuw+VLV3CjoBB+wwJly9edr4Doy0IdJz1AQRRFNB+rAEYFQaF8uN5qFQoFmMwwtJy5CUubSbacblYcKo4XwtgiPXhGoVAgdsloHL9wCq2t8rdRvc3DdeYQQoibiouLw5xZc/He3/75QM9LTE6E0WBEeWk5jAYjEpISAACj0kZh+sxp2PT3z53QWvdgtVrx+T++wlPPPAGWHfztxMrySmg0GhQXlUAd4gvBKD3Rr9VkQUXuDfhNjZad385Qooe5ywwmUTfodnkyRYAaYoI/mo/K97XkQnygHh2C0r3y63X6hPkjcOIwHDxyyJHNdHsU0gghxE28/vrrOJZzDAXXCu77OYuWLsTyVdlISU3G8lXZWLR0Ye/3nnv5OVw4fhGFV4sc3lZHzqM2WEf2HkWgLgBTpk+2q564hDjU19WjvqEB2hA/aMOlZ82vPl0CJlwjO7O+KIhoPl4JRXqQ7O1Qb8ekBaOrpg2G8hbZcv4To9Fc3YTmonrZctFZw1Ghv4mSkhIHttK9UUgjhBA3ERQUhKeefBp/+e2b971sk0arwaOPPYLX/+XbePSxR/rdmgwI0OHZFzfi3T/80+HLQB3ecxT7vz6IsuIK7P/6oF3zqA2G2WTGF+99heGpw/HmH97Crq93w2gcXFCct2gu6uoboEsORfS0ZESMG7gPGd9pQtXJIvhPlR+d2nX9FiwKAYoYv0G1x1soOAaKjO5BBDKDAxgVC78ZsSjanQeBlz5PWZUSwxaPxP6jB8HzvDOa7HYopBFCiBtZuHAhRCuwc9tuh9S3OHsRtCot9m074JD6elSV9Z/7zJ551AZj77b9YMAi7+JllJaUY+/OfTi0L2dQdZ05cRYtpjaMf3Uuhk1OAiuxkHrV8SJwCQHgQmzTfAi8Fe2XaqE/WIr2S7UQeCtEi4Dmk5VQpIcMeDtUtAiwFjbDcqYO1sJm2fDiDRQxfrBoGbSdq5Ytp00OgqjjUHVc/qpv8IhICOEcTp8948hmui0KaYQQ4kYYhsHzzz2Pj979GM36ZrvqMhqM2LF1J3S6APz9N++irlp+yZ4HEZvYf+4ze+ZRe1Cd7Z3YvmkXUoanQKvVID4hDhqNBpUVVYOq7/9+/TtEz02FSquWLGNq7ULNhTL4Tx52ux1XG9CRV4euG41oPliGpl2FaLtQA6svCyZy4PnahNJWCEUtEPUmCEUtEEq9uyO8QqEAMz4UrRdqwLcYZMvpsuJw83QxDE0dsuXilozG6ctn0djY6IwmuxUKaYQQ4maSk5MxbcpU/OPN9+yqp2d6jmZ9MxgFg//58a8d1ELbNB0LH52PxJQELHx0/n3Po+aIdTy3f7EL4zPHY+KUCTAYjKgor4TRaERcfOwD17Xtq+1o7GxG7IwRsuUqjtyAZkQIlLrbt5P5JgOsnTwEAw+Bt8BQpEfz0XIo0qWn7xBbTADLQBGoBlgGYov0EkreQuGvgpisQ/PhcskF6AFAGaiBZmw4SvZckS2nCfZFyLQ4bN25TbacN6B50gghxA2tXbsOP/jhD3D50hWMHTdmUHX0nZ7DYrHg8qUruHrxGkaOkZ6Tqi+jwYjDe46iqqwKsYmx/eZL65mm40H1rOOp1qpRcOUGADzQHGot+hbs33oQf3zr9wgNCwPHcbh65SqiY6Mxb9GDzQknCAL+8uc3EbNwFJQStzgBoOtWO+rzb0I7Phz6g6XgQrTwTQsHF6KFYLBA4K1gOBaCVYCgYaEMkZ6yRBGohnjLYAtrVgGKQNUDtdlTMaOCYNhbBUNJM3xSpNcw9c+Mwq1N+dAX1CJk1DDJcsOmJ+PG307i4MGDWLBggTOa7BboShohhLghPz8/PLn+Sfz1d28OupN03+k5LLwF2SuW4b0//BM8f38ztztjcIC963hu/uhrzJ4zC9Ex0dBqNVi1ZiW+8Z1XsezRpdBoHmw+t4/e+xhdLI/oSUmy5cpzroMJUcNwoxH8rU50XqlH59UG+KaFQ5sSBIZjoVArYe3iwaTIL+7OJAWAGR4IRbAGzPBAMEn2LwbvCRRKBorxIdAfLoVglh4coFAy0M2KQ/GeK7CapM9TRskietko/OyX/w9dXdJzrHk6CmmEEOKmZs2aBa3aB9u++npQz79zeo7v/ei7CA0Jw+6v9tzX850xOOB+1/G0Wq0wdBnR1tKGxvpG1N6sw9WL13D8wAk8ufEJu9thsVjw7nsfIG7xaNk51lrKG1F/pRLmpi5Y2s1QBmmg4FjwTQYwHIuQZakImp8ERqsEE+sHVuZWJ2ALIWxqEJSTI8CmPlwT3TLD/GAN4NB6Rv480sQHAqFqVObekC0XPDwC5jAGf3nzrw5spXuh252EEOKmFAoFnn/uefz8jZ9j9rxZCAsPe6Dn90zP0YNhGHz3h9/Bay98A9PmTkVouHygiEmMwbXLN1BWWA5DlwGBwQGoqaoFb+ZhNplhNvPgeR5mo9n20WQGb+bBm3lYrQJMJhN4Ew+z2QzebIHZZIKhywhDlxF1N+uh0WpwaMdh7N92qPu5ZvC8BSajCaJg62ukUqugUqmgUqvAcRye3PAkQsNCH/xg3uHN378FMUCJyLHSSzaJooiCr85CZBVQCCIszQYYSgClvwpciNZ2TDkW2sQg6I9XQLk41iGhS7QIEEpaIbaYoAhUg0kO8Jowx4wLRdv+m/AbFdY7SnYguqw4VH92DREZcfAJk56TLnl5Bt770wd4bOUqJCcnO6PJQ4pCGiGEuLG4uDjMmT0Hf/vzO/jJf/7Y7voSEuOxaMlC/Or/+z+EhIf0Bi7ezMNs7g5UJlv4MhpNaGtphdnMw8dHi31bDuLonuNQqVTgOCU4VXeA6n6sUquh4jioNWoEBQVBq/CBzo/rDlpqqFTdny+1ha6eshyn6g1j6u4w1hPK5Gb1l9J3Lc64hDjMWzS3363Qro4ufPblV0h6egIY5nb4sZotqL9Uia6GNviE66AJ8EFXYwe4YX7ggrUwwHYlzHdMBHzTwnuf13KqCkjUQeEzuLVD7ySUtEIoagaUDMRbtlt57Iggh9Q91BS+HJAagKacUkSsTpP8+Sr91fDJjEDRrjyMfWaGZDltkC9CZ8bjp2/8Bz5+78NBnS/ujEIaIYS4ucdXP45/+eEPcO70eUycMsHu+p7YsA6JyYlgGAU4lS0MiaKIi+cuob62HnEJsZi3cC78dH7gOA5KpfKB3vw4pRLBISEwmUwQBNfPA9azFqdGo8H1q7bVG5Y9urT3+7/51e+gjvZDWGpUv+fVX6pE3blysCoWbRV6GNu7oB0VAoveAEuLEUo/FXzHRMB/3O3nmRs70VmiB7PUcYuoiy0mQGkbASq2mGyPvQgzIgimfVXoutEI35HSV4f9MiLRWHANjfnVCBsjPcVLwqyRuPjHQ9i1axeys7Od0eQh4x3XTwkhxItptVo89eTTePP3b8Fstn/KBpVKhWkzp2LK9CnInDgeYzLScbPyJi6cvYCmxiacPHYK586ch4+Pz6CvZg2lnrU44waYP62psQm79uxFwtIxUCj6vwV2NbSBVbHwCdeBN5ph7DAiaEEyfMdEgAvzvesKGgC0nKgChuugUA9+7dA7KQLVgEWwhTOLYHvsRRSsAorMUOiPlkOQGRygYBn4z4pF8b4rsBilB88wShYJj47Bz3/1Btra2pzR5CFDIY0QQpyE53lcvHgR+/fvx8WLF+1aymbatGkICQ7FF5986cAW3iYXbDxNXEIcjEYjKgeYP+1Xb/wa/iPCEBR/9xUcn3AdrGYrOuva0FHTAp+0MLBqJfzHRSF4fhL8x0WB4W6HMVNNOww1rWBSHXsrkkkOADM8CIogDZjhQWCSvW8EKBPhAyFEjZaT8ueZJiYA7DBfVORcly0XkhIJxGrwhz//0ZHNHHIU0gghxEny8/ORl5eHhoYG5OXlIT8/f9B19Qwi2PzZVtRW1zqsjUajEbu+3o3CgiKUl5WjrKRs0BPDuot5i+ZicfYiJKYkYnH2ot7506oqqnDkxHEkLE4b8HkR4+IQOTEBgAhliBa6GbF3LfvUQxRFNB+rAEYGQcE59q1UoWTAjgiCckok2BHeOwKUyQhF+7UGmBs6ZcvpZsSi9kolOmpbZMslZWdg07bPcf26fKDzJB7fJy0hIQHBwdIT43kKV+2DSqXq11HWWViWdfiCzgMRBMEht3/kuOqYuWJfiGs1NjaC4ziEhoaisbHR7mVsoqKisHjhYrz5x7/h5//zs0HfhmSVShgNRhzcdwj7du9HdVUNhsUMA9N9+y/70aWYv2geOOXdbxE9z6soq0B8YjzmL5rXb1F3AFByjulAP1gajaZfH7Qev/z5/yA4Ixq6yIGvfLEqJSLHx6PyeCGCs4fDWKRH55V6KDgW5hrbbbSe/mjGihaYWg1gpj3YiFtym8JHCYwKhD6nFJFrx0DBSAwi8FPDd2IUSnZfRuZLc3vLcXecZxqdFuFzkvCTn/8UX37yuUtet53N40NaVlbWUDfBoZzZ90OlUiE19f5mGvckhYWFTgs3rj5mztwX4nqhoaGorq5GY2MjeJ5HaKj9U0esXLkSP/zXH+LksVOYnjXtAZ9tm9YiICAAh/bl4OC+HFSWV6GlpQXxifEYMSoVKanDseH5DZI1bPliKw7uy4FWq0FxUSn8/P2xas3Kgbc2yCV7jAYj9u3ej7KSMiQmJ2LR0oV3BcEHdS3/Gi5cuYwJP5RfJaHmbBkQqIImJgBdN5qg4FhwIT7gm7rAN9nWnhSF7qtoaUFQsPa/ZosWAUJpK8QWs+um3HCTfoZMahBMFVXoKmiEbmyUZLmAzBg0FuWjo7AJCTJLeMVNT8WlC4ewectmPL76cWc02aU8PqTt3LkTer1+qJtht+DgYGRnZzt1HbKevyqcfcwSEhKQlZXl9O30HDNn/rXkqmPmin0hrpeeng7AdkUtNDS097E91Go1ntnwDN7649vInDj+gcKLIIgwmYwAFCgsKIJKpUJCUgIunr+EkqISxCcmIDY+BiaT9GjCnufFxsehvLQcRTeKBiwviuKgX8961hzVaDW4evkaAPSb720wfvnG/yJ8ajx8gvwG/L7VbEHt2TIUfn0RflOjIfBWcCFamGvawDd1Qex+DABdRU3gLVYw8dLzdz2InkXXXTrlhpusealQAIrMUDQdLYUmKQisVvoqrG9WLPK2nIRvYhBUvmqwDAs///4/T4ZlkLgiA7/87a8wf958BAV59tQlHh/S9Ho9GhoahroZHsXZx6zn1q03/Wy8aV+I63Ach/Hjxzu83okTJyLncA42ffApnnvl2Qd6riCIAEQkJMUjPy8fnIpDWFgYYhNisHT5EixcskB22oye55WVlMFoMCI+Md7h02z0XXO0vLQc5aXldtV3+sQZFFWUYvK6JZJl6i9VovzANYhKBUwVLejUqXtHcvJNht71OkWrgObjFVCMDXLYnY++i65745Qb98KEamGN0KI5txwhC6QnpFVF+IKN98eNnRcx/JFxUHJK+OHu0B2UEIa6FH/83+9/g//6+S+c2XSn8/iQRgghD6ONz2zEj//tx5i/ZP6gOvkvWroQAFBeWo7F2Yvu+5Zi3+clJCX0PnakxOREXL18DeWl5TAajEhISrCrvl/99/8halYy1P5ayTLtVXoYWzqhTQ+DubYDbWerAQC+aeH9RnS259XBqlKAifK1q0199Vt03Qun3LgfzNgQdOytgl9aONRR0lco/afFoP6TfESOj0NQYrhkueSlY7Hl919j7WNrkJGR4YwmuwSFNEII8UDh4eFYnr0cf/3dm/jv3/3XA1/VuXPJKGc/70E4Mgju3bUftfoGTM66eyBBX12NHVD4cDDXdsCiN0AZpEXnlXoAtwcLCGYrWk5XQTEl3KH9h3sWWe/bJ+1ho9AogbQgNB0qRdQTYyUHEbBaDr5TolG0Mw+TvjFfsj6VnwaR84fjJ//5U2z7fIvs+qzujEIaIYR4qOXLl+Poj47iyMGjmLNgtlO35YzO/HL1OiIICoKA3//2T4heMAKcRiXdBn0nOpraETAjDh2X66AM0kKbHARLs7F3sAAAtOfVQghQgQ2TviI3GD2LrrtLZ/6hwiQHgC9vR8flun6rOtzJd3QYmgoaUX26FKHLpUfXxk5JwaULOdiydYvHDiKgkEYIIR6K4zg8/9zzePMvf8WkqRPh6+e4W3B3up/O/IMJcs4YJNDjy083o83SiZFTU2TLlR++Du2YUARMiQGjZtF5pR6WZmO/wQJWI4/Wc9VQzB7mkLaRuykUCihGB6P1fI1sSFMwCqhiA9Ba1XTPOkWzAJ2/zpHNdCkaSkYIIR5szJgxGJE6Eh+997FTt9O3M79GqxmwM39P4CopKsWOLTuxb/d+h9Q7GBaLBW+/9XfELh4FdoD53np01rWisbgO/pm28OWbFj7gMlBtZ2sgRvqAeQj7i7mSWNUB3+EhsmUEowWG/AbEz5SeigMAqi+UIj4oGgsWLHBkE12KQhohhHi4p596Ggd2H0RpcanTtpGYnAijwSjbmX8wget+6h2M995+H2atiGGZibLlyg5dgzYjAozGFuQYjr1rGShLhxntV+rApHn+xOnuTGwxQVHXBd2kaNly7RdrEZocCb9hgZJlrLwVdYeK8ePv/8ijpzai252EEOLhQkJCsHLFKvzlt2/i13/+lVPelO6nM/9gRmUOdpAAz/PobO8Ez/Ooqa5GW2s7Oru6YOgyoKOjE//88GMkrh8veyzaq5vRWq1H6LwxstvqvNYAMUwDxm9oV1LwdkK+HoETomXnSrN2mmHMv4W0l+bK1lV1qghjk9IwbdqDTvjsXiikEUKIF1iyZAmO5h7FgT0HsGiZ/Kz6gzFQZ36jwYgd23aiuLAEYeGhSBubhlFjRqGyrAIxcdHo6urCR+99DEOXAQaDAUajEV1dBphNJhgNJphMtn88z8NsNuPo4Vz842/vwmKxgud5WK1WWC0WWKxWWK1WCN0frVYBgAiFQgElpwSrZMEqleBUtn8KhoHRbIRfpPwoSaWGgyiIEAx8v2k27uSTEoKWc9UQO3goKKg5hdBoANNilu2LBgDtZ2sQmR4HbcjAkxIDgMXIo+FoKf749kdOXcXHFSikEUKIF1AqlXj+uefx29//BiFhobBarOjq6kJnRye6ugwwdhnQZTDA0GWAyWiEwWiE0WCEyWiC2WyG0Wi0hSXeDJ7nwfMWWC0W8BZLd1CywGoVYLVYYBUECFYrujq6YDKabR2+GcDX3xcBwQHgVByUnBLnLp8Hp+KgUnNQqVXg1H0+D1UhUK0Dp1JCpVZDqVJCpeJ6y3O9n6ttn6s5qDWq7jIqWz0aNbQaDZqammCxWPodj9/88I8o3p6HjGellw7UhvghZmIi6nMrEZwtvfwbF6xFQEYUWi81gp0pHyLIgxNFEeJlPYKmxIBRSYdlS4sR5qJmxH1zgmx9FbkFmDVxJsaOHevoprochTRCCPESI0eOhEqtwivPvwY/nU93WOoOPT3hSMVBpVH1PlbpOPiqA8CpQsCpbF/nVBw4Tgm1RtX9/D7P6wldKhW2f7ADtZX1iE4chtrKOsQmRWP9N9e6bH8ZmaskL//kOfzgyZ+g/tpNRIyOkSwXmzUC9X+9CUNpM7RJtiWEBN6KzqsNvSsNaFNDwGhYiOVt4JWAcmKE89fWfIiIdV1gzQL80iJky7WfqUb0pCSo/KVHDJs7jGg6WYkfbvqjo5s5JCikEUKIl7hx4wY6jZ14c/cfodE6fxRifGo8KourUVtZB5PBjKh497nK5B/oj0fXL8WOrQcQOjwSLKeE1WxB/aVKdDW0wSdch4hxcWBVSqQsGYuCPZegjtWB4Vh0Xm1A55V6KDgW5po2GCtaYG0zgQvSwnRVD2uAGspRNIjAEURRhHhFj6DpcbLB13yrE5aqdsSsnCpbX/mh63hk4XIkJ0svL+VJ6E8BQgjxAqIo4sNPPsT8tXNdEtAAYOqCychaNh2xSdHIWjYdUxdMdsl279eypxfDX6VF2aGrAGzrc9adK0dnXSvqzpWj/lIlACB4RCQCwgPRcb4GgG2tTgXHggvxsQW1hk4oOBaahECwWg5CUeuQ7ZO3ESvbwbEsfFJDZct1nKpG3IxUKLXSkxIbmjvQcqEG3/nmtx3dzCFDIY0QQrzApUuX0NBcj2kLXReU1Bo1Zi/PwvpvrsXs5VlQa9xrDjGGYfDKvz2HmiPF6LjViq6GNrAqFj7hOrAqFoaGdgC2SVRTloyF8Uoj+BbbLU6Rt4Jv6oLIW6EK94XIW2FpMoAL9oGi3QyxzTzEe+f5RKsI8WozgmbESy4DBQCm6jaIehOiJstPp1K+/xo2rHkSw4Z5z4TDdLuTEEI8nCAI+PjTj7H4yQXgVM4ffWgymnDqwBnUVNRiWHwUpi6YPGQBzWQw4fjuk6iuqEFkbAQmzs6Eqs8SUCnpyRg3IR1FWy8iYmws2ir16Gpog9VshTbcv98tUF1UEFoPlSNo+XAA6NcnzVDYBL7JAJ8QLXy6eLRevAVm1jCPHz04lITSVqh1GmgSAyXLiKKIjlPViJ81EiwnHVk66lthuNGMV37/ihNaOnQopBFCiIc7efIkjKIBE2aNd8n2Th04g9xdJ6DWqlBytQwAMHv57VGUrgxxx/edxIl9p8CpOZReLUNhXhF0wbp+ge2FHz+D76/7N4gTYhE5MQGGhnZow/0RMS6u9xYoq2JhNVlg1XfCVNF611QQfR+LFgGdhY2wVnVAEefvlP3ydqJFAK43I+iRUbJB11jWAtYkImJ8nGx95fuu4tVnX0JwsHf1FaSQRgghHozneXzy2cfIfnGxy2ZWr6mohVqrwrD4KNRU1KK2orbf9+8V4qRYrVbwZgt4Mw8Lb/vY88/CW2A29f+6lbfi2O4TaNG3YnTmSFw7X4C6qnokjIxHeUEFAGD64qnw8fPB488/ii8+2o5JP1zSb7H1vrdAuxraEBgfitbjN6GJD5ScDkKhZBAyNwn1ewohRvlCwVHPoQclFLZAG6WDOlp6XU1RENFxuhqp89LBsNLHuKWyEag2YuMzG53R1CFFIY0QQjzYoZxD0AZrkDZptMPqtFqs4HkLeJMZPG+BpTcsWcDzFgiCiIbqW2hq0IM3mqEL1OHQtsO9Aev0wbNovtUMX50vWvWt+PqDHci/cA1mo62+O8MXb+ZhMVsgCAJEiFBySihVHJQcC6VKCbb7o5JjwXJK2+cq2+cdQica6m6BP8HDZDLBT+eHyNgI1FXVo66qvnefFqyeh0NfH0XpgXyMWJ7Z+3WfcF2/W6BhY2PBajm0n6lGwEzpqzea+ED4xASgK78J7Pgwhx37h4FosgKFrQhcky5brquwERpOhZBR8n3MKvdew/dfex2+vr6ObKZboJBGCCEeymg04ovNn+OpH66TvWXUWNuIXZ/thaHL0Bukev/1C2MW8DwPURABoDsQdQcklW1mf07NAQwgciI6DZ3wDfKFkTPhRkWxbfZ/FQtVqAp8owWd5i4otAyGZcQgYVqiLXT1hixb+OJU3O0Q1r2N++7npQB4E4+TW07j8Ps5GD1hJLraDairqofZaEZkbP95t1799+fx82//L6ImJEIXZZsTLWKcLYj1vQUaMjIK5948CH5kKLhQH8nNB81OhOHDixATdVDQwuv3TShohm9SEFRh0qFKtAjoPFOL0Y9MkB1UcKuwBj6dSqxd47r5+VyJQhohhHio3bt3Iyo5AslpSZJlBEHAW798B9pwH6RMSgGn7g5F3WFJFEWUnC9BS10LwuLDMHbeGGj8NA8Wlu5gNppx5XA+GisbERoXijFz0vt15ncklUaFmWumg1MpUX3xJqYtmoL6mw29fdL6ihseh8kzxuPq1gvIfHUuFAoGrEqJYZP7Hz9WpUR81khUHS1HyKqRksdB6adC4OQYNF+oBzM3mgYR3AexiwdK2xD49DjZch1XG+AX7I/AZOmrlKIo4ubeArzx+r9DpXLO+TXUKKQRQogH6ujowNadW/HqG8/Lljt35AKsjIBHvrt8wD5r5/dcQPG5YnBqDo03G6H2VWPCkswBarp/Ko3K7joe1Ni5Y3Dyy5MQLSLGz8y4a5Rnj43/8hS+v/bHqL5QhpgJ0hOeDpuchLpLFegqaITvKOmg4D8uCh3XboEvbwebKN2/itgI15rhPzocygDpVQMEsxVd5+uQ8cQ02eDbdL0GkeoQZGdnO6OpboF6OxJCiAfasm0LUjNTEJ0YLVmGN/PY/skuTF83XXJQQWPlLXBqDmFxYbagVtnorCY71bUT16HUcLhyOh/Hd5/EuSMXBiyn0Wqw/tXHUb7jCswGk2R9jJLB8OwMdJ6shmC0SJZTsAxC5iUBV5psfa0eAqJFgLWwGZYzdbAWNttGat7P89rNUNzsRMBk6WW6AKDjUh2C40LhHyM9UlMUBNQeLMb/971/BctKr/fp6Tz+SlpCQoJXDLn1hn0ghLhGU1MT9h3Yi+/99luy5Y7vPQWfUB8kjI3v/40+FydC48Jws6Aat6pugTfxCI0L7ff9QRMdUMf9UgCNVY0ICA+AaBXQ1KBHQ/UtKJUDv8XNWzEbB7ccRtm+fKQ9PkWy2pCUSISnDkPr6WoEzkmQLKeJCYBfYgg68pvATgi3d2/cnlDaCqGoBWAZiLcMAAA2Nejez8vXI2BcFFhf6VuT1i4ehrwGjHphtmxdDXlVCNcEY86cOQ/SdI/j8SEtK+vew7o9CfVpIITcy5ebv8T4ORkIjZJeSsdoMGH3F3uR/f1l/V5X7uyEPXbeGCgUCtyqbERYXCjGzE13zOuQAr0DEJyNZVlEJESgprAG/qH+KD1fiqDgIISEhEg+5ye//yG+ve6HsM4fh/AU6dGDE5+Yjf2/+AKWRgPUkdJzooXOS0bXP85BSDSBCZG+lecNxBYTwDJQBKohtpggttx79QVBbwRzywjdcvmRmh3naxA2ahh8wqSPtcBb0XCkFC+tec7r3zM9PqTt3LkTer1+qJtht+DgYGRnZ0MUXfnnJyHE09TW1iL3RC5+9Jfvy5bL2XYEEcMjMGx4nzfF7vczq9Xae6WLVbIYtzCj33OtFjtv2ylswQkKuOSKmiAIGD1rNARBQGNVI/xDdMi/kI+FzfMkb/PqQnWYOX8Kzrx7CHEzR6CroR1+kQGImpAAVnX7rZHjOKQtn4grORcR8vho6ZGGahZB0+PQdKEa4vxBDiLwkMChCFRDvGWwhTWrAEXgvTvti/l6BE6KBqOWjh2WdhOMBU1If3W+bF1158oRHRiJpCTpATPewuNDml6vR0NDw1A3gxBCXGLTZ5swY/lUBARLd1Jvb+3Awe05ePxnqwcuIAImgxn5R/LRWNWI0NhQpM9Oh0rjmCWlFH3ulwpWARazBRbeAgtvtX3e/djKW2AxdX+9ewoQi9kCa5/HPXOo2eZUs4I388iYMxYpfTr9i4IITqXE+EXjere56aef4tieE5g8b+KAbTQbzUgYEY/9Xx3CtdoWBCaGQl9aD7PZfNdoz4QZI1F45Ao68uvglx4xYH0A4Jsejvb8OphLWsGmBA7+ALo5JikAACC2mKEIVPU+liLUd4HtsMA/I0q2XMeZagwbnwBNoPS0JxYjj1vHyvHUo+sevOEeyONDGiGEPCxKS0txKf8i/u1b/ypbbv+XB5A4IQGhMdK3Q/OP5OPc9rMAw6DobDFuVdxC0vik7vBkhcXEw2KxhSqr2QK+O1jxfUKWhbeAN3WHLb5v+LL2BjJRECDCdsWOVSnBKhmwKhZM98S0jJIBy9keMxxre6y6/TnDsWA1LFh/JRhODaXJgh1v7sIrv3sRWp12wH1jWAZzn5uDHb/ZibRJo+Drf/d8XOeOXMDZnPPQBfqjqaEZlmFWcFqud9H1O+sb8ch4XPhnLrRJwWB9Bg6zCkaB4HnJqNt8FWKMHxQa73yLVSiZ++qDBtimyRDz9QicGguFUnqsIt/UBXNZK2K+NVm2vtqTJUiJTkJkZOQDtdlTeecZRAghXujjTz/G3NWz4eM3cDgBAP2tZhw7cApP/fcTsnVdyclHWX4FfIP9YDFZ0NbagZr6W2CVLBiVLSBBAbTebIaxwwTfEB+EjYwCF6CCUqWBmmNt4Uppm2uM7fNYqVKC6Z6Ytudr93v7z2LiUX6iBK03WxAQE4iE6clQqvuHovbaVhzZlIslryySrGfY8GFIyEzAnk/3Y/VLK+/6fl1VPVQaFZJHJ6Gl6SJaSxsQNDwS2vCB+0LpYkMQMToGLSerEDhf+jabOtIP/iPD0H65Cexk6atuDwuxuhNKK+A7Wn5ARfvpasROHQ6Vr/SkwOYOI/RnbmLF2mcc3Uy3RSGNEEI8QH5+PspvluHJnzwuW273p3sxavZIBIRJ34IydhpRXVyD0OHh0AT5wmrikTRnBFLmjuhXrjjnBhpLGqDyU8PYZgTno7qrzEAUgK1/1SD6o5WfKEHJ4UIo1Uo0Ftu6sqTMHdmvTNrKcch5YxfGzktHdKr0FCRZT8zEBz/8EBPnZCJ+eP8lniJjI1BeUAElp0R0QhRuVtTCPz6odwWCgSTMH4WzfzkIU0071MOkO7YHTI9F5z8vQbhlABMmHai9Xc9VtKCsRNlVA0x1HRDquzBsrXwfs5pjxRgzfPRDNRsCzZNGCCFuThRFfLTpIyxcP9+2LJOE2so6nD95EVNWyN8yOrvjHGImxmPEsjEIigtC0pwRSJh+98SubTebwao5BEQHglVzaKtuHvQ+WEw8inMKcP7DUyjOKYDFxA9YrvVmC5RqJQKiA6FUK9Fa3XJXGbW/BqmPjMHed/ZDsErP0aX112La2mnY+t72u8pNnJ2JqQsnIyY5GiueewSLH5+P5hv1UCilwwTno0bS/DS0Ha2UHbnKajgEzYyHeOGWy0a4uiOhrA0qrQraZJn5zkQRHSdvIm7miLuumPZlbO5EW14dpk2e5oymui0KaYQQ4ubOnj2LVkMLpswfuBN8jx0f78K4xePgo5PueN3Z0olzey9g9KoMpMwdgcynpyJl7ggoBxh1p4sJgtXEo7W6BVYTD130/fVDGkjPFbKWKj1KDhei/ETJgOUCYgJhMVnQWt0Ci8mCgOjAAcslzkiBSbTg4r5LsttNn5MGQSni5IHT/b6u0qgwffFUPPbiCkxfPBVPfmcdlF0iKk8UydYXMS4OPloNOi7XyZbzHRUGta/aNp/YQ0i0CsC1ZgTNiJe/ilbVCkWnBVETE2TrqzlShIlpmfD3l76C6Y0opBFCiBuzWq345LOPseSphbIzq5cXVqDwWjEmZk+Qre/4VycROSEWfhH3XsIoYXoykuaMuOtqm8VkQXHODVz48BSKc27AYpKekb/H/Vwh69lm8pxUBMYFI3lO6oBX+ABbJ/0x6ybgyOe56GzplNwuw9gGERzYnIO25rsHBfRQKpV49vtPomrvNRjapOtTMAqkLMtA17k6WDuk5wdTMAoEz00CrrfY1qu8w2Bn7fcUQnErNCG+UMdKn2eiIKLjVDUS54wGo5Q+tzvrWtFV2IzJk+SvEHsj6pNGCCFu7OjRo4AGyJg2pt/XTUYTTh04g5qKWgyLi8L54xcwccVEqLTSc1a11Lfg8pErmPPvS+5r20q1csA+aOUnSlB6+AZYNYeGwnrUXa2BT4AWupggJExPBjfAVbmAmEA0Fjfc8wqZUs3d1QdNSlB8CCInxODwJ0ew7LWlkuUiEiOQMj0FuzbtwfpvrJEsN37mOCQnx6NkRx7Sn5wuWc4vMgBRGXFoOl6FoMXS63+qwnyhSwtHW14T2Gn9RyMOOGv/CO/oayXyVuBGC4JWpckOGDGU6MGBRdgY+WWianKKMG38ZGg03j1J8EDoShohhLgpnufx6ZefYtmGxXdNynrqwBnk7jqBm6XV2PXpXhRdLUHGgrGy9eV+fhyxM5OgDZK+HXo/+vZVM7YaUHO+As1VzSg9fEPyNub9XiF7UKMeGYuivFJUXb8Js5HHhb0Xse+d/biw9yLMxttXsKY/Pg2F+cUouVYmW98rP30BLfl1aCyqBQBYzRaUHb2Owm0XUXOmFFaz7aph/JxREGo7Yaxqla0vYGosGL0JQl1Xv6/3nbUfLHNfs/Z7CuFGC3xiA6GK8JMsI1oFdJ6pQeI8mQmCAbRVNsFS04XM8ZnOaKrbo5BGCCFuau++vQiODsSIjOF3fa+mohZqrQpRcZGor65HWFI4lJz0zZGGigYUXihG6uLRdrerb181c7sJKn/tPQcX9Fwhm/D0VKTMHSnbSfx+9AxEuLL5IgISg7H3nX24fPAy8vbnoaG8AXn785B/JL+3vMZXgxlPzsDW97fDwkvfng0OC8LiVfNQvOUirBYras+Xo/zYdXTUtqDuXDnqL1Xa9kfDIXnRGLQfrZS9VcmolQiZnQjx4i1bP61uikA1YBUeaNZ+TyAaLUBxKwKnSY+SBYDO67eg9dUiOFV6vjNRFFFzsAhZk6aD4xwz0bKnoZBGCCFuqKurC19t/QrLNiwe8JbRsPgomAxmXD13HTxvweiZo2TrO/JpLpLmp8rOQ9XjXn3O+vZVGzYhDhqd2iGDCx5E70CESj3aqlvQ2tKBKzlXwGk4hMaGgtNwaKxq7PecUTNGQqXjcHzPSdm6H3tpBTRWFhW5BeioawWr4uAboQOrYvtNdhuaHg2/QD90XJIfRKAdHgJNoA+EG7cDLJMUAGZ4IBTBGjDDA+85a7+nEK43w294KLhg6alHBN6KrnN1SFogfzu0uagByjYRY9LHSJbxdhTSCCHEDe3YsR0J6XF3ze/VY+qCyZixeCqabjVhysrJGDtP+o3sZsFN3CypQXJ3/7J7hbCePmdN5U24/MU55PxqT79yPX3VMp+eiqkvZSFl/ijZqTycod9ABA2HoOQQ1JTVwdBmQGNVI3gjj9DY/isuKBQKzNk4B4e2H0FLY4tk3QzD4MUfbcTNAzeg9OVgNfPorG+D1WztN9mtQqFAytIMdF2sg6XNKFmfQtE9iKCwDWKH7RZsz6z9yskRYFODZGfj9xRiBw9FRQcCpsbKluu83ICAyCAExEuviCEKImoPFWLO1CzZATPezvPPCkIIcWM8z+PixYvYv38/Ll68CJ4feH6wvlpbW7Fjzw4sfVJ6Rn21Rg2VWoW49DgsfW0JVJqBb5eJooicT44gZcno3luM5SdKUHzwOirPVuDCBydx6u+5/YJaT58z0SrAoO9Ci0x/s76BTWoqD2foP1UHj/BRUUhZNApQKBCeEI6MhRlIn51+1/NCY0Mxes5obP9ot2z9aRNHYVT6cDSX1CNh5ij4RQUicmLCXZPd+oT5I3piEtpyq2Tr44K1CMiIgnCpUbacJxOu6uGfHgGlv/TVWsFoQdelOiTMk7/y23i1Gv6CFqmpqY5upkehkEYIIU6Un5+PvLw8NDQ0IC8vD/n5+fd8zuYtm5E2bTQiY6WXFTKbzNj5+R7MWD9d9pZR6cVSNDe1ImHG7StcrTebYWwzwdxmgMVkQc35yn4BrG+fM1GEQyazdbT+AxFsV/DSVmTAyJsxfNJwZC4eL7lg/NRVU1BRWonCPPk50V7+t2fRXtIEha8SqSvGY9jkJLCqu0NoXNYIoMkEQ5n88dFNjgbbboFQ3XH/O+ohxBYTFHVd0E2SXgECANov1iI0ORJ+UYGSZQSLgLqcEsydPvuuATMPm4d77wkhxMkaGxvBcRxCQ0PBcRwaG+WvpNTX1+Pg4QNYvH6BbLmjO48hKDYIsaOkby0JgoCcT44idXk62D7zUAXEBMHcboDFbAGrUkLlr+4XwHr7nMUGwSfEFwzLDNjfzGKy4Ma+azjw8x3Y/587ULjv2n3NmeYIvQMRNty+gqfUcBi1OgP73z8oOzhApVUh6+ksbPnndvBm6SubuiAdVj29HFc+PgGLWbo+Vq1EypKxaM+tkh9EwLEImZMI8VKT982Llq9HQGY0WK10B39rpxnG/FuInyM/xUr9xQpE+AQjISHBwa30PBTSCCHEiUJDQ8HzPBobG8HzPEJDpfvhAMAnn36CyYsnIjhMugN+V3sX9m4+gOlrpefyAoDrxwtgFnjE3jGbe8L0ZAybEA+lmoNap4EmQNsvgPXcwpzzoyUYu2YCghNDBuxvVn6iBFe3XkT99VrUX61F/tZLklNwuEp0Zhy4IDXOfH1WttzwSSnwj/DHkR3HZMs9ujEbAWoflB6UvwIaPDISARGBaD9fI1tOkxQEbbgfhOvuc1XSXkKjAUyzGf7jo2TLtZ+tQWR6HLQh0lNzWE0W3DpahjkzZsteIX5YUEgjhBAnSk9PR0ZGBsLDw5GRkYH09Lv7SfUoKSnBqbMnsWD1PNk6D27NQUx6DCISwiXLWHkrjn6ei5Erxt41D5VSrcTUl7KQ+cxUxE1OQMq8kQN2+L9Xf7O2m82w8gJUfmqo/dWw8tYhvyWqUNhWIjiz+xxa6ltky815ZjZy9xxHU71eshzDMPjBf7+O6pxCdNySnhNNoVAgeclYGC/fAt9ikC0XPCcRipI2iO2ePzeaKIoQr+gROCUGjEq6g7+lxQhzUTPiZt89OXJftWfKkBARi2HDhjm6qR6JQhohhDgRx3EYP348Fi5ciPHjx8vO9/SHP/0BWStmwE/nK1mmVd+Gw7tyMX2N/ELTeYcug9WpEJE28JudIzr862KCwHIMzB0mmNpNYDnWZVNwyPGL0CF+VgoO/PMQRFF6gfOgqCCMXTQG2/65Q7bc6PEjMW5CGoq3XYIoyizoHuyL2GnDbQuwy9SnDNAgYEI0hAu3ZMt5ArGuC6xJgF+6dP9JAGg/U43oSYlQ+UuvGsB3mdB0sgKzpmU5upkei0IaIYS4gcrKSpzPO4c5j86SLbfn830YPn04giKlw5DZYMaxr45j1IqMe94yGsw6nD0SpicjbeV4RIyOQkRaFNJXjhvwilzPxLPnPzyF4pwCGNsM/R5bTPce8fqgUpeMRl1VA4rPyd9+nfzoZNTV1uPq2euy5V7+9+fRVdGMuivyozhjpqeAabfCUCJ9dQ4AdBOGQWkUIVZ57iCCnqtoQdPjZKcQMd/qhKWqHTHT756Uua/aE6UYmZB6zy4BDxNau5MQQoaYKIo4cGg/5q+dA41WDUEY+GrNrdpGnD5yFs/87wbZ+s7tPg+/uECEJIfdc9t91+FsLG4AgAHX6xyIUq3EiEWjMWKR/CoGPRPPKtVKNBY3oO5qDTpvdfQ+tm3z/tbrvF+sSom0NeNx4INDiB8TLznSU6lSYvYzs7D9vZ1IHZsiOZWJj58PHn/2UXy5aQdCU6PASZRjOBbDl43FtW3noYkLvOsWoMBb0Xm1AXyTAdrYALRdbIQiyhcKzvOumYiVHeBYFj6p8qGq41Q14makQimzrqypzYDm89VY/cRzjm6mR6OQRgghQ6yiogLtpjbMXDIdgszttB2f7MKYBWPgFyzd8bqrzYDTO85g6nfn3te2+67D2Vrd4pQ+ZX0nnm2tbkFzeRN8gn17H7dWtzh8mwAQNTYGlcdLcWrLKcx6QvoWWtL4JOTnXMWhrYexZL1tbjqz0YxzRy6gsaYJyaOSMHLiCCxcMx+Hvj6KsgNXkbp8vGR9QSkRCBwWjFtbr4ML8QEXooVvWjgYjkXn1QZ0XqmHgmMh8lZwWg6Wq01gx907UN+LKIqAteef0Pu5eMdjWAXb6FKhf3nRKkJhFQEBgCDYPreiu5xwu66e55kFBC0fIbv2pqm6DaLehKjJibJtrzlShHEjxiIgwDtWXnAUCmmEEDKEBEHAyTMnsOiZBeBUHEwm04DlqsuqkX/hGp777bOy9Z36+jRC0yIRcJ99w3QxQWgsbnDqsk4BMYG927CYLAhKCEHnrY7exwHRgQ7fZo/0NZnI/e99SJs1GiHRIZLl5myYjU9+sgnjZ45DREw4zh25gFP7z0Drq0VV8U10tHdgysJJeOUnz+ON7/4akRMToJO55ewfFYia06UQ4niYa9psXxsXBb7JAAXHggvxAd/UBXWID6zlzbByDKBQ3A5YAroDU98gJfYLTHeFMFEEoICCVUChZKBgGdvHOz5nlCzAovvrrO1rGgYKjr37ub3PU9xdn4oF6yPdx1IURXScqkb8rJFgZdaV7WpsR8e1W5j2zMp7/jwfNhTSCCFkCBUVFYH1YTBuxljZcts+3IEJyzOh8ZPueN3W1I6LBy5h1o8X3/f2e/qQtVU3Qxcd5JRlnXrqbK1ugV+4PwTeCoO+E4JFQMKMZKcuJeUb6ofEBSOw/72DWPeTNXf10TMbeeQfyUdjVSNCEkKw+d1tePWnL6Kuqh4qjW0B+8baJtRW2dbnTBgZj0nTx+HalgvIfHUuFIqBb1OaO4zwC9fB2GaGMkwLvsk24pML0cJc0wa+qQsib4V6mD98UkNgrGwFq1J2hyS2TxiSD1x3PWYUsle2AEDBMLBaLE4ftGAsawFrEhExXn6x9ZqcIkweOxE+Pj5ObY8nopD2EEpJSUFwcLDT6o+KinLJdnQ6ndPqvlNEhPzIJXs58zi5w/bIwCwWC85fOo8Vry2DWi29lE7x1RJUlFZh7rfkp+Y4/uVxDJscD99Q6duhPawmC8pPlKD1ZjMCYoKQsWaiXUs6We6oL2FGMpQq21WWnolnAaA4pwDlx0ugVCthMVnAcmzvclXOkrJgJI6c3oOCkzcwanr/vm/5R/KRtz8PnIaD2WBGW0M78k5cRmRsBMoLKlBbWQfBIiAqNrL3Oc/+8Gl8f+2PUX2hHDETkgbcpk+4DuoALQwVXbDoDfAdabud6ZtmmzaFbzL0uw3qOzwUClGB9su14Ju6wIX4wK/7e06hUMCZs5CJgojO09UYuShDsp8fALRXN8N6swtLv7sUGo30HyADcfbrsjugkPYQ6fkLcto0+aH7juKq7Thz8d2eY7ZokfQaio7k7IWEezqkZ2dnO3U75P588cUXiE6OwsyFM3rPtTuvboiiiG0f7MDkxyaDkwkzTdV6XDtVgLk/W3bvDSsUd3TmvwUASJk3+M77A9Y3wGCAO/unOas/Wl8sxyJ9/QQc+iAHSeMSofa5HYgbqxrBaTiExoaisaoRw0YOw45P9uD1X7xm+36tHskjEzFy4u3BFBqtButefgwfvfMVwkdHQ6W9O2D3rPGp1mlxq7AW2hTbH0YMx8J/3MCTvnbk16P9ch0UHAtjte0WqW68c+YLc/ZrTce1Bvj7+WFE1t3z9PVVvfkqfv7jn2HD0/KDYR5WFNIeIj0v/jt37oReLz883BMEBwcjOzsbVqvVadvoOWa5ubloa2tz2nZ0Oh2ysrKcui8AYDabUVhY+NCvh+cOjEYjfvvH3+Kpf10Hs9k2qakoineFtKtnr0Hf0ozls5fL1nf0s6OInz0cGp323hsXbet39g9Lzd19mgbn7vpaBix3Z/80Z/ZH6yt8ZCQCkkNx7MsTmP/M7UEVobGhqC6oRmNVI3gjj9FZo1FXVIcj23PxyDPZ0Gg0CAgIQFNTEyyW29OTzH4kCwe2HEbJ3isYtXLiXdtjVcreoNZe24Km7YUIe3y07JUx861OKDgWqlAfmBu7em+ROpooMXrYYfVbBHScrsak9Vm95/ZAmorqoGjkMSFzAoqLiwe1rZSUlME20yNQSHsI6fV6NDQ0DHUzPEp5eblTj1l4eDiyslwzgaPciyZxnS1btiAyKQxJoxIkp9wQBAFbP9yOaY9P7bf25p3qSupQdrUS839+v1dIxd4BAy19BgzcT0Trua3ZdrMZuhhbHzalWnlXfVLhq2//tIDoQKf2R7tT2upxOPKLPRgzOx3h8bbbj+mzbStANFY1IjQ2FOmz05E6eTg++tHHmDArE0kjpUclvvazF/GzV/8LLRMTERhz96CE+kuVqDtXDpW/Bh036tByrAJckPauW509VGG+MFa3wtxo66/GhdxH4HZD7fn1CAwLQNDwCMnRyqIoomrfdXxj/QsQBAFGo9HFrfQMFNIIIcTFOjo6sHXnVrz6xvOy5c4duQArIyB1aqpsucObjiJ50UhwMvNQ3WmwAwak5lXrW19AtK1P2kD69k9zNW2gD1KWjsb+9w7giZ+tA8MwUGk4ZC7uP52GSsNh0qpJ2Pr+dnz3l9/q/XrPtBx1VfWIjI3AxNmZmD5nMi5uuYDMb86/6wp1V0MbWBULn3AdTK0GtJ+uhmqYPxgV22/EZw+/9AiIgtAb4u41i787EsxWdJ2vw9gX5stOpNyQfxNBCn/MmiU/efPDjkIaIYS42JZtW5CamYLoxGjJMryZx/ZPdmHWc7Nkb09XXKlA3c0GzH9x8gO1oWdZqL6krpL1JTWvWt/6FIBtOgk3XPEoaU4qjp4sw9UjVzFm7hjJcuMXjcP13Os4m3Meix5bAMAWmk/sPYX2lg6caDuFwrwirP/G47iw4aeoPluC2Cn9Z9T3CdehrVKProY2qP01MDR1wNplhjoqBHzT3bczGY51Wh+0ByUKIsTu+dR6P1r6fK3n67zQr5y5rgNhSZEIjJee4FYQBNQeKMS/vvA9p/eN83QU0gghxIWampqw78BefO+335Itd3zvKfiE+iBhbLxkGVEUkfPJEaRmp4N1wCjA+1l9YKB51e4Md4nTk6GUmOF/qDEsg/T1E3D477lImZgCrf/AtxQZlsHcZ+di++93YvoC2yCouqp6tLd0oKO1A2YTjytnriI1Yzie/vY6/OMPHyM8PRZq39sjFHv6pBka2qEN94cuIQRl+6/B5NcBhSje1+3MgcOSbeJZkbfaPlru+P4d5dA9ca1o6Z5PrTtswSpA4G9PbtvveYLYexwYjgWjZMBwSjBKFizHglWyYPp+5JRQciqwUX5IyJK/Ulp7rgyxgcMwefKD/WHxMKKQRgghLvTl5i8xfk4GQqOkrzQYDSbs/mIvsr+/TPaWUeGZInR0dWHiFPnZ3O/X/aw+MNBt0jvDnQL2jRR1ttCUcISkRSD38+NY9MICyXLRqcMQNz4em9/biuynlyAyNgIn2k7BbOKhUnPw0/mhrqoej724Avu+OIiSXXkYvWZK7/NZlRLDJt+eosNqtqCtshkttU3gArXgK1qhL2/rDkwiGFEBgbf2hiWhZ1UAAIyyOyjdGY6ULBiVEqySAcspbQGKU4PhWCh9bI9tIcsWru58vi189X3cp5ySveecaw/KyltQf7gEb/zgP+65riyhkEYIIS5TW1uL3BO5+NFfvi9bLmfbEUQMj8Cw4dK3vgSrgMObjmLEI2PAsI4ZrXs/qw8MdJv0znDX6oSlpRxt9MpxOPzGLoyZnYaolIGnxACAmU/MwMc/+gTpU0Zj4uxMFOYV4cqZq/DT+cHX3wdt+jZsfmcbJs2dgC/f2wb9pEQEJ4QPWBerUmLs8zOhv14LMIp+wYvTqBAcGgyrKAAs+gcn1vFhaajcPFmM9MRRSE9PH+qmeAQKaYQQ4iKbPtuEGcunIiBYeiLm9tYOHNyeg8d/tlq2ritH8iGqgGHjYh3WvsEOJrgz3N3vklRDSaPTYMQjY7Dv3QPY8IunJPv9+eh8MG3tNGz+xzZ88z9fwZOvr+sdPNCmb0NjbRM62jphNpqROioZxVsuYOJ3FoJhBr79rFRzCB939wz8Sk6JwJAQmIwm2fVbPRlvMOPWsXL8+P/971A3xWPQZEmEEOICpaWluJR/EXNXzJYtt//LA0ickIDQGOnboRazBblfHMfIFRkOvWXUc5Us8+mpSJk74r5XIEiYnoykOSMQFBeEpDkjXDqthj0SZqbAJFqQd+CybLkx89JhYaw4feAsVBoVpi+eisdeXAFdsA4aXw0iYyOg0qiQPDoRbKeAyhNFLtoDz1KVW4jp46YgOdkzzg93QFfSCCHEBT7+9GPMXT0bPn7SncX1t5px7MApPPXfT8jWdWHvRWgifBE+MlK2nKvceQvUU27MKRgFxqybgNy/HEXqlFT4Bgy8diTDMMh6ciY2/fRTlFwrhcVsga/OF51tnTB2GVFXVQ+z0Yxh8VF49ntP4a//8y4ixsZCq/N18R65L1O7AfozN7Hh1z8a6qZ4FLqSRgghTnb16lWU3yzDrGUzZMvt/nQvRs0eiYCwAMkyxk4jTmw7hVGPyC/ITu5PUEIIIjJjcPjjI7LlGipuQYCIE3tP4/TBc8g/cw2NtU0IjQzBsIQoTF04GRNnZyIzaxySk+NQsiPPRXvgGapybmBR1nwMG+YeU4x4CgpphBDiRKIo4qNNH2Lh+vmya2/WVtbh/MmLmLJCflqCMzvOIWh4KIIS7p7hngzOqEfHojivFDev35Qsc6vqFiKSwmE0GAGIUHIsNL4a6IJ1eOzFFZi+eGrvQuKv/PQFtFytQ2NJnYv2wL116TvQdrke69esG+qmeBy63UkIIU507tw5tHS1YMr8u9d37GvHx7swbvE4+OgGvuUGAJ0tnTi7+xxm/utCRzfzoabyVWPkyrHY994BbPzlhgGX4AqLDUN1QTX8Q/2hr9bDYrbAbDQjMvbuVQGCw4KweOU8HNx8AkHfWyy7pJcriYIIwWKFYLHCylsh8N2fW7o/560QLAKsvAWCRYDA9/mexQqBF27Ps2YRIfI9c631TBlitX3s8zxr9/aeXvMUQkLoD4sHRSGNEEKcRBAEfPzpR1jy1ELZmdXLCytQeK0Yz730rGx9JzafRGRmLPwjpEeHOtL9rEDgLeKmJqLqRCku7L6ASY9Muuv76XPSIUJEfWk9rhy6Av8g/95bnAN57KUVOH7gNCpyC5A0N63f926HJaH3IwOA7RRh6DDAyvN3B6ee0MNbge4JZ2G1BSXRYpuwVrSKveUF3haYbIGrO5RZrVBAAaVSCU7FgeM4qFRqcBwHtUoFlUoFrUoFFaeGRqWCWq2GWqWGWq2GVqeBurusqrtsz+d9P6rVasmvkwfnnb9thBDiBo4ePQpogIxp0ssPiaKIbR9sx8QVE6GSWXuzpb4Fl4/mY86/L3VGUwd0PysQeAuFQoEx6yfi+O8OYeSMkfAP9u/3/b5rfE5cPhFfvfEVxkxJ773FeSeGYbDxe0/ij//xNrou3+oXyESrbYoNVqmEkmXBKpXQajRgYmJgtVqhZJXwU6uh4jio1b62oKRSQ61RQ6OzfX5nKBooPPV8fmc5uWXGiHuhkEYIIU7A8zw+/XITHvvWCtk3xRuXClFX24DF/7JEtr7cL44jdmYStEHSt0Md7X5WIPAmAdGB8AnzQ9XVKozOGi1ZTqXhbMs13WM+s9amNoxJHY35c+dDqVSCZVkolUoolUowDNPvvAgPD8eGDRtQXFwMo9HosH0ino3iNCGEOMHefXsRNCwQIzKGS5YRBAHbPtyBKY9PgZKT/pu5oaIBhReKkbpIOjg4gy4mCFYTL7sCgTdpKKiDpc2M1CmpsuVObT6NyXMmQhckfdvZ0GnA4c1HMTtrNkJCQhAQEAA/Pz9oNJrekEbIvdCVNEIIcbCuri5s3vYVnvvpBtnJZvNOXEanqQujZ46Sre/Ip7lImp8Kla9r+/UMdgUCd3K//epEUUTB11cwc/U0KFXSb42NVY0oP1+Ox/9vpex2c3ccx7CIaERE3D2wgJD7RSGNEEIcbNfuXYgbHYv44Xcv/9PX/q2HkLk8U/aqimAVUJ5fgRmLpRcCd5aB1un0NPfbr6720k0ojALSZ8uvKXnii5PIWjIdfjIT1bY1t+PYjlNYu2qtfY0nDz263koIIQ6mUAC8iYcoirLl0iaMRuHJQtlyDMtg2oqpuLb10j3rI3fr26+OVXMD9qsTrAJubL+C2etmyi5WX1NUg4biBsxcOn3A75uNZpzYewp/+rc34aPygZ+fn8P2gzycKKQRQoiDPbL8UTTfbMHlU1dkyy16bD466ztQdLZYttzUFZNh0RtRfaHSkc18KNxPv7rK02Xw1WqRMilFsh5RFHHisxOYt2I21NqBbzufO3IBR3ccQ3lBBUxGE/Lz8x22H+ThRCGNEEIcTKVS4YVnX8TWd3bAaDBJluPUHNa9/DiOfHAUZoNZspySU2Lx8wtxbfMl8EbeGU32Wvda/N3KW1G08ypmr8+Sve1cfrkcnY1dmDLv7jnUetRV1aOptgnRUdHQarVobGx02H6QhxOFNEIIcYIJEyYgNXEE9n1+QLbc6AmjkJKaiFNbTsmWSxyXiPjUGNzYRVdnHkRPv7rMp6ciZe6IuwYNlB0tQviwUMSlS/cfFAQBJz47iUWr58uOwlVr1WhubIVarQbP8wgNDXXYfpCHE4U0Qghxkuc3Po8ze8+htlJ+DcfVL67C1ZxraLwpf+Vl/sb5uHmiDK01LQ5s5cOLN/Ao3nsds9dnyY7CvXGqEIyVwbgZ8ova15bVYUzaGERFRSEjIwPp6fKDEAi5F48f3ZmQkIDg4OChbobdXLkPKSkpXnHMdDrXLI3jSt60dArLsrBarR6/DcB2JcVslr4dKSUsLAyrHnkMX729Fd/8xSuSQSA4LAhLVi9EznuH8fhPV0uW04X6Y8Zj05D36TnM+N582WBxL4N/5tDX7qgtlhwoQOLoeESlREmWEawCTn1xCqufWQGVSnpFiNLrZWgob8R3vvkd2XJSXDlVh0qlcsk8ba76/fRmHh/SsrKyhroJDmXPi+791j1t2jSnbWMoyK2J6Cn8/W1L0MTGxg5xS4iUwsLCQQW15cuXI+doDs4duYBJcyZIlpu3Yi7OHD6H68cLkCYz2/2kZRNx5Ug+Ks+UIX5q0gO3B4Bt+Kk3GcT+GNsMKD9ShOf+6xnZBdDz9+cjNDgE0+ZNlXx9FkURH2z9BD/98U+xZs2aB25LX858DwBsAS01VX6yXuI+PD6k7dy5E3q9fqibYbfg4GBkZ2c7dYh9T93edsy84S81juMAeM/PJiEhAVlZWU7dH1dsA7h9ng32ygPHcXjp+Zfw2z//BmkTR8PHTztwOZUSG779JP70xptIGBsPja9m4AoVwMJn5+OrP2xD1JhoqHwe5KqNwnbJSUTPf06icO3FtEHsT+Huqxg5eQQCowIlX0N4E49TX53Gc99/Bm1tbZJ1XT13DQ2Vt6BWq7Fz584HakcPnU6HrKys+3oPMBqNyMnJQUVFBeLj4zF37lxoNBLnyx16zmNn/9644vez53fTm3l8SNPr9WhoaBjqZngUOmbuy1t+Nj230525P67YhqOMGTMGY0dnYPcne7H65ZWS5UZkpGJE+nAc/+IE5m2cK1kuZlQMUsYk4vrXlzF2/cT7bocCInpSmnMjWs92XOXB9qezsQM3T5Xjhf99FqIg/cwLuy9iWGwUohOjJNfTFAQBX7+3Exlp41FUVPSA7b4tPDz8vu8M5eTkYO/evdBoNCgoKAAALF269IG25+zfG0/6/XRnNHCAEEJcYOOGjbh05DKqSm7Kllv+9DIUnyhGfZn8G9ucp2aj7nwVmiuaHNnMh8KNnVcwft5Y6EKl+7UaO4y4uPMilq5bJFtX3okrsHYJGD5ceo1WR6uoqIBGo0FcXBw0Gg0qK2n+PG9FIY0QQlwgODgYax5bi6/e3gZBECTL6YL8sWDVXOS8nyNbzjfQFzMfn4Ern56XvRpE+muracGt/DpMWTFFttzZHecwYuxwRMVHSpax8Bbs23QQ0yZPd+mC6fHx8TAajaisrITRaERcnPzyY8RzUUgjhBAXWbpkKRRGBqcOnJUtN23hFDC8AvmHr8qWG7cgAyqwKD8uv2IBue36tsuYvGwitP4D9w0EgHZ9O/IP5mPR6vmSZcxGMz763adoqmmCXq8Hz7tukuG5c+di8eLFSExMxOLFizF3rvStceLZKKQRQoiLsCyLl194Gbs/3IuOtk7JcgzLYOVzj+Dk5ydhaDfIllv0/ELc+DofpvaB+0yR25pKbqG9sgUTlmTKlju99Qwyp49DSGSIZJmT+0/jxJ6T8PXxw+XLl126BJRGo8HSpUvxyiuvYOnSpfc9aIB4HgpphBDiQiNGjMCUzKnY+eFu2XLxw+MwZkIajn16XLbcsOFRGDV5BK5tzXNkM72OKIq4vu0ypq+cCpVWekRsc20zik8WY97KObL1nT10Hlq1D2JjY8FxHC0BRZyCQhohhLjYU08+hasnC1BWUC5bbsn6hSg7X4aaohrZcrPWz0Rjfh2aSm45sJXepf5aLaytJoydN0a23MmvTmH6winQBflLluls70RNaS0CAwPR2NhIS0ARp6GQRgghLhYQEICn1j+Fr97eJjvPn6+/L5asta1EIFilBxFo/bWYvT4LVz49L1vuYSUKIgq2XUbWmpmya2/WlzXg5pWbmJU9U7a+w1tzMXrkaEyaNAnh4eG0BBRxGgpphBAyBObPmw8fxhfHdp+ULTdpzgRolBrkHZC/nZk+Ow1+Pj4oPVzoyGZ6hZvnKqBSKDFy+gjZcie+OIE5j8yC1ld6UEFLYwtO7T2L6dNmYPz48Vi4cCHGjx/fOyE1IY5EIY0QQoYAwzB4+YWXsf/Tg2jVS89mzzAMVj33CM5sPouOZpnBBgyDRc/NR/HuazC0dDmjyR5JsFhRuCMfs9fNlJ0mo+paFZqrmjFt4WTZ+g58mYORw0ciICDA0U0l5C4U0gghZIgkJSVh1rTZ2P5P+aWEohOHYfz0DORuypUtF54QjjFZabi6+ZIDW+nZyo+XIDBYh6Tx0uuciqKI45+dwIJVc6FSSw8qaKi+hcvHr2LyRPkgR4ijUEgjhJAhtH7depRcKkPhZfklhRY9Ph81+dWoulYlW27mmhloKWpEQ0GdI5vpkSwmC4p2XcPs9bNkFy4vOV8Ca6cFE2fJT82x77ODGDt6LHx8fBzdVEIGRCGNEEKGkJ+fHzY+9Sw2v70NFt4iWU7rq0X2k0uR8/5hWHnpwQZqHzXmPz0H+Z+dh2CRLvcwKDl0A9HJUYgZGS1ZRhAEnPz8FBY9vgCskpUsV1VyE8WXSpE5Xj7IuYLRaMTu3bvx1ltvYffu3ZLrihLPRyGNEEKGWFZWFkJ14Ti49bBsuXEzxiJQF4ALey7Ilhs5fSSCQwJRdKDAga30LOYOE0oP3sCsdfKLll8/dh0alRrpk0dLlhFFEbs/2o8JGROgVqsd3dQH1rPAenl5Ofbu3YucnJyhbhJxEgpphBAyxBQKBV57+TXkfHkEzY0tsuVWPvsIzm+/gLZG6cEGCoUCC5+dj9L9N9DV1OGEFru/on3XMHx8MsLjwyTLWHgLTn91BkvXLZIdVFCSX4qGigaMHTvWGU19YLTA+sODQhohhLiB+Ph4rFu9Djve3yVbLiImHFPmTMLRj47KlguJCUHmggzkf3HRkc30CIbmLlQeK8XMx6fLlss7cBlh4aFISU+WLCMIAnZ/vA+TM6dAqZSeY82VaIH1hweFNEIIcROvvvIaKq5VoTBPfhDB/FVzcKu0EWWXymTLTVs5FZ03W1F3pdqBrXR/N3blY0xWGgIjAiXLmAwmnPv6HJasXSg7qODq2evo0hsxatQoJ7R0cGiB9YcHhTRCCHETPj4+WDx/Cba+s0N2EIFaq8ajTy/D4Q+OwGKWLqfSqrDgmXnI//wCrDLlvEl7fRvqLtzE1JVTZctd2H0RSamJiBseK1nGarFiz8f7MXXSVNnboa5GC6w/PNznrCOEEIK0tDT4Kv1weJv8nGhpk0YjMjICZ3ecky2XMikFkTHhKNx7zZHNdFs3vr6CCYsz4RfkK1mms7ULeXvysHjNAtm6zh+9CNaqRHKy9O1QQpyJQhohhLgRhUKB2VmzcXTbMTTV62XLPbohG3l78tBc1yJbbv7GeSg/XIT2eunBBt6gubwJ+sJbmJQ9Qbbc2a/PIn3iaETEhEuW4c089n92CNOnTJe9HUqIM1FII4QQNxMUFIQxo8Zi27s7IIqiZLnQqBDMXDQNRz44IlsuKDIQU5ZNwpUvLsiW83QFX1/GtEcnQ+Mrffuv9VYrrh+5jgWr5snWdWLvaQT6BCE2Vvp2KCHORiGNEELc0ORJk1FbXI+rZ6/Llpvz6Cy01bah+GyxfH2PTALfZED1RfkVCzxVQ0EdDPWdGL94vGy5U5tPY/KciQgKC5QsY+g04PDmo5g2ZZqDW0nIg3GP8cR2SEhIQHBwsFO3oVKpYDabnboNZ+8DsU9ERIRT69fpdC7Zjqv07I8zz+ueup39u+Pqn0nf/VmRvQK7/rkDozJHQq0ZeE1JpVKJx55bgc//8RUSxiZApR24HKfmsOi5hdj+9i5EjIoEp+Gc0n4b194eFEURBdsuY+bqaVCqpN/WGqsaUXa+DI//30rZ+nJ3HMewiGiv+X0knsvjQ1pWlvxs0p6G+j64F4vFNiJu0aJFLtmeq7bjKtnZ2V6xDeD2ueAsgiAAuHt/2jvbcWbfOax87hHJ585cOB15Jy7j7PZzmPWk9GtiyoRkxKfG4Mauq0hfLX/FyZPUXKqCwiQifXa6bLkTX5xE1pIZ8NNJDypoa27HsR2nsHbVWkc3k5AH5vEhbefOndDrpTvX2ishIQFZWVlO305wcDCys7O9ur+IJ7JYLCgoKHDJJJYKhcKrfv4sy8Jqde7aka46ZhaLxekhzWw2o7Cw8K6pHtatWYdvf+9bGDdzLKLiIiWf/9jzK/CL13+FUVkjERoTKllu/sZ5eOdf3kPc1ETohgU6qvlDRrAKKNiejwXrZoNhpXvw1BTVoKG4Ac++9pRsfTlbjiAlMQVBQUGObiohD8zjQ5per0dDQ4PT6u+59eDs7RD35Yo3aEIADNitQqfTYcXylfjyrS34xhsvS15tDwwNxKLH5uPQuzlY8++PS5bThegwfeVUXP7sPGZ8d57HX72vPFUGPx8tUialSJYRRREnPjuBeY/OhlorvfZmU70e5w9dxJNr5YMcIa5CAwcIIcTNPbL8ERiajDifK7/E07wVc2Btt+D6cfmF1Sctmwh0WVF1Rn7FAndn5a0o3JmPOetnyU42W365Ap23OjFl/iTZ+vZ/fhBpI9Lg7+/v6KYSMigU0gghxM1xHIeXnn8Z29/dBUOnQbKcklNi3SuP49gnx2DqMkmWYzkWS15YiOtb8mDucu6gKGcqPVqIiJgwxKVLr10pCAJOfHYCi1YvgJKTvnlUW1GH62cKMXGifJAjxJUopBFCiAcYM2YMxo7KwO5N+2TLpY4djlFjR+L4Fydky8WlxSFlbBKub89zZDOdQhRFWC1W8AYexjYDOps60FrdjJJ9BZj75GzZW7ZFp4vAWBUYN3Os7Db2fLIf48aMoyWWiFvx+D5phBDysNi4YSO+98PvYvK8iYhJipYs99hzK/CL1/8HabPTEJEgPav+3Kfn4O//8i7ipiUjKO7eU5mIogjBIkCwWGE1W2Hlbf8GfmyBlRdsj3krrLztsWCxQuAFCLwVgkWA1Wyx1Wm2wmqx9pYX+Nv1WXkrIIqAQgElp4RSpYSSUyJ9xmhEpUTBahl4gIrVYsXJL05hxVPLZW+Hll0vR9WNm5j75Px7HgNCXIlCGiGEeIjg4GCseWwtvnp7K779369JBo+AYB2WrV2MQ+8ewrr/t1aynF+QH7JWz8Dxvx+HLjqgNzxZu0OW0B2yrN0hS+AFACIUCgZKFQuWs4UlpVoJJcd2hyeu+yMLTmX7XKNSQsn5gPPjukMW2xu22N7QxfY+5tRKsMrucqru7yuVYJRM/6tm9xjzkP//t3fv4VHV977HP3Nbk5BkooMQiQYCcpOtqGC59DTeigjiZburXDZ4KRR7etT2VE+lT1s3Vrtre9r91GefR+zu2UjOUy0icikFhQoC0i0KFvBSK2gBARGDjTAhIZnbOn+EGZKQzEwuM/Ob4f16Hh5m1vzW7/dba31nzSczmZXNf1Gpz6eLRw3rsI1t23r5uT/qysu/JI8nndeOAzqPkAYAOWTypMnauHmj3nz1LY2fMKbDdldN+Yq2bnhT7236i0Zed2mH7UbdcIWKzy2WHY3G36FqDkuuM+57LI9cHpdcblc6Nq1HhZpC2r5iu2bdNz3hx6Ef7NytY0cCuuSaxNdYA7KBkAYAOcTlcuneOffqp//2U40c+w8qKmn/wqwul0szvjVVT/3rbzT4ysHq5Stst53T5dTw8R2/05Srdq7bpQv6l2vgxZUdtolGo1r73Csad+U4uVzmB0+cffjiAADkmGHDhmnsFWO1+rdrE7YbOLxSV4wZqT8t+VOGZmaGxhON2rlmpyZNvT5hu7dff1eRhqiGDBmSoZkBnUNIA4AcNPOfZ+ovr7+v/bs/Ttjulrtu0v639uvwh4czNLPs2776LQ29dIjKK/t12CYcCuuPizdo/JgvJ/xSAZBNVCYA5KDS0lLNnD5Ly369Mv53P9tTUlqsW2fdpFef2ahopON2+aKutk7vrn9XN9ye+Jua2zf+WYXuQg0YMCBDMwM6j5AGADnqq9d9VYWOXvqvtVsTths/cayKrF7atd78a6J115srt2nUly9X7/N7d9im6WST1r+wUePHfDnn/ywW8htfHACAHOV0OjV3zr368ROP6rLxl8p3rq/DdtP+++369/kLNGzsUBWd0/6XDTIpds21cDiscDD2L6JIOKxQU1iR0KlloYjCwZb3m9uFT923w1E1NQab2wQj+uSvn+ih//3thGP/19o31OecviovL8/Q1gJdQ0gDgBx20UUX6arxV+sP/+8lzfyf0zts139whcZc/SW99txrmnzf5Pjy2NX8I6fCUDwItbkfCUUUago1twudDlax+5F4mGq+xlo41LrPtuuFW1yg1mO55fF4mv+3mq+l5vF65LE8p5Z7ZHlP3y6yCuSxLFnFlnoVFaqxqVFOl1Nuj1tTbruhw7AqSfV19Xpt5Z90203/1KPHAUgHQhoA5Ljp06brOw99Wx/95W8a/A8Xddjupn+erJ/c/4T+7/3/2So8SZLD6TgdkCzP6X8t7lve5scty5LH8qjAKpLldcvjs+Lt3Ja71TqxPpvXPR3EYstdbleXP3J0Op3yer36+9//rnA4nNI6m36/RQMuGKDzzjuvS2MCmURIA4AcV1xcrLtn3qMlv35e/+tX3+nwYrO9igv1w/8zT/V1DW3CmPusuE7Ysb8f15vrtmv612ZkeypASvjiAADkgaqqKvUu6q1Nf9iSsF1RSZH6lveRv8+5KiktVkGh96wIaJK0/sWNGjZ4uEpLS7M9FSAlhDQAyAMOh0PfmD1Xry7dpC8+P5bt6Rjn6OGjeudP72nMlR3/KS3ANIQ0AMgT/fv31/XXTtTKhX/I9lSMs+75DRo5YqR69eqV7akAKSOkAUAeueP2O/TJ7k/1/o4Psj0VYxza+4k+2rVXo64Yle2pAJ3CFwcAII8UFhbq63d+XdW/WaQh/36RPJYn21PKiGgkqlAwpFAorHAwpHAoolAwpHAorJeeXafRl42W1+vN9jSBTiGkAUCeGTdunF559RW9umKTbpiW+I+M96TYNdfCoXBzYAqGFAqGFQ6FFGwKtVkeC1ThVvebr8UWVqjpVMg6dZ22eF+nglewKaRIOCI7auvYF8cUCoflkORyueV2x/655HK5VFLs08j/NjJj+wHoKYQ0AMgzDodDc2fP1cM/fFgjrrxYvYp7xd9VahmSYmHnzOWnA1Ko6VRIirdrDkqhVuHpdH/RqC2HwyHL8sjjsWRZ1qnbnubbp5Z5PJa8p/5vXl6oYsvbfD224lib5vW8Xm98/Zb/l5SUaPjw4Vq+fLlqa2vlcnX9mmuAiQhpKaqsrJTf709b/+nsO5tjpVNZWZkk5c1HGC6XS5FIJNvT6DHRaFTBYDCtY1iWJaeTX61tz8CBAzXta9P0H48sPBV2mkOR5fXK43HLsryyrNNByXvqfrFV1Hy75HQgiv1rG5haLm95P1Nhyev1qnfv3iovL1dBQUHax0un2PksU9L9OhDrP53jZHqfZYPDtm0725Pojt/+9reqqalJW/8jR47U9ddn7uOCv/3tbzp58mRa+i4qKtLAgQPT0jfQnj179qQtqFmWpaFDh6albyBb0vkaIPE6kGt4Jy2J2AvMmjVrVFtbm7Zx/H6/pkyZonRm5ti7NFu2bFEgEEjbOP369dOoUaMyNk66j00mVFZWqqqqKi+2RTpdz+l8lyvWd77ss0zItzqTpJKSEnk86ftyRKbOZz6fT1VVVWl9DZAy9zogpf/Tgdg+y2eEtBTV1tam9R27TNq/f39GtiVT4+TDsYl9JJAP25Jp7LPU5WOdZWo70n0+69u3b0YDR6bOz+mU6X2WDfwyBwAAgIEIaQAAAAYipAEAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABiKkAQAAGIiQBgAAYCBCGgAAgIEIaQAAAAYipAEAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABnJnewK5orKyUn6/P239p7NvIJ+l+7mZTyoqKrI9BQCdQEhLwrIsSVJVVVVGxnM4HBkZB8h1LpdLUuaemwCQaYS0JILBoCRpzZo1qq2tTds4fr9fU6ZMkW3baRsDyCeRSERS+p+b+aSyspJQC+QQQlqKamtrVVNTk+1pAGiD52bq+FgYyC18cQAAAMBAhDQAAAADEdIAAAAMREgDAAAwECENAADAQIQ0AAAAAxHSAAAADERIAwAAMBAhDQAAwECENAAAAAMR0gAAAAxESAMAADAQIQ0AAMBAhDQAAAADEdIAAAAMREgDAAAwECENAADAQIQ0AAAAAxHSAAAADERIAwAAMBAhDQAAwECENAAAAAO5sz2B7vL7/RnpP93jlJWVpbX/TI7l8/kyOk66j00mZKrOMiWT9Zwv+ywT8q3OMiEfXwOk/KiBTO+zbHDYtm1nexI47YMPPlA4HE5L3263W8OHD09L30B70lnPlmVp6NChaekbyJZ0PmckqaioSAMHDkxb/+hZOR3S9uzZI6cz/Z/YOhwOZWI3hcPhtD45peag5nan/w3UTO0zl8ulSCSS9nEyIVP7LFMyUc+WZWXkHJBP8q3OMiGfXgMKCgo0ePBgbdmyRYFAIK1jpZvP51NVVVW2p5FWOf1xZzAYzPYUck4mTgJApnAOALpm//79qqmpyfY0uqVv3755H9L4ERQAAMBAhDQAAAADEdIAAAAMREgDAAAwECENAADAQIQ0AAAAAxHSAAAADERIAwAAMBAhDQAAwECENAAAAAMR0gAAAAxESAMAADAQIQ0AAMBAhDQAAAADEdIAAAAMREgDAAAwECENAADAQIQ0AAAAAxHSAAAADERIAwAAMBAhDQAAwEDubE+gOyzLktOZPznT5XIpEolkexo9IlPbkk/7LBqNKhgMpn2cTD1v8unYZEo+7bN8Ogfk03Hxer3ZngI6IadD2tChQ7M9BaBH7dmzJ61BzbIsnjcAkCNyOqRJ0po1a1RbW5vtaXRbZWWlqqqq8mJ7MrUt+bTP/H6/pkyZkvZ3uGL9c2zMk0/7LJ/OAfl0XKTT24PckPMhrba2VjU1NdmeRrf5/X5J+bE9mdqWfNpnmcaxMU8+7bN8Ogfk03GRTm8PckP+/EIXAABAHiGkAQAAGIiQBgAAYCBCGgAAgIEIaQAAAAYipAEAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABiKkAQAAGIiQBgAAYCBCGgAAgIEIaQAAAAYipAEAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYyJ3tCSB/+f3+jPSf7nEyoaysLNtTAHpcZWVlWp+fFRUVaR8nE2NkUmx78mFbzobzJiENPc6yLEnSlClTMjJepsbJBIfDke0pAN0WOwdUVVVlZLxMjJOpbcmUfDpv5jNCGnpcMBiUJG3ZskWBQCBt4/Tr10+jRo1K+ziZ4PP5VFVVJdu2sz0VoNti54A1a9aotrY2beNUVlaqqqoqreNkYoxMim1PPp038xkhDWmzf/9+1dTU5M046dS3b9+8P9ng7FNbW5vW52bsI7t0jpOJMTIptj2cN3MDXxwAAAAwECENAADAQIQ0AAAAAxHSAAAADERIAwAAMBAhDQAAwECENAAAAAMR0gAAAAxESAMAADAQIQ0AAMBAhDQAAAADEdIAAAAMREgDAAAwECENAADAQIQ0AAAAAxHSAAAADERIAwAAMBAhDQAAwECENAAAAAMR0gAAAAxESAMAADAQIQ0AAMBA7mxPoLv8fn+2p9AjYtuRD9uTqW3Jp31WVlaW0fE4NubJp32WT+eAfDouUn5tT6bPm9ngsG3bzvYkAAAA0BofdwIAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABiKkAQAAGIiQBgAAYCBCGgAAgIEIaQAAAAYipAEAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABiKkAQAAGIiQBgAAYCBCGgAAgIEIaQAAAAYipAEAABjInUqjaDSqw4cPq6SkRA6HI91zAgAAyEu2bauurk7l5eVyOhO/V5ZSSDt8+LAqKip6ZHIAAABnu4MHD+rCCy9M2CalkFZSUhLv0OfzdX9mWRKOhnQsWCuHHDqb3w+0JdmydY7ll9vp6bBdY6RR+wMfyiGHnI70fDIetaOyZavSN0QFroIe7/9EqE7bal6XU065emAbInZUUUU1pu+XVewp6YEZdk1tY61W7f+93A63XE5X1ubRkUg0orAd1i2Vt8pf4E/a/rP6o3rqnWfkcbrldqZ0WkooHA0rFA3rvpGzVVbUp9v9meJg4Ige2fwreVweeVwdP3c7EoqEFIqE9PjV31WF7/wuzWFf7Sf61oofy3J5ZLmtLvVhsmA4qGAkpKdvm6+B/guyPZ20+tvRA5r2nw/JcnlU4PEmbNsYalIwEtKSb/ybLurTP0MzzE+BQEAVFRXxbJVISmfD2EecPp8v50NaJBiSUw450hQ6coFtRxWVLZ/lSxjSrIilYhXL5XClNaRF7Ih8Pl9aQpoz5FDRyV7yODxy9cCLfyQaVsgOyefzZTWkha2wCosL5XV55UlwDLMlFA2pKdLUfM4oSH7OaHA1yltUoF7uAlmu7r/wByNBNYQbVeIrka8od89ZbZWoXp5eXhVZvVTgTvyi2p7GcJPqgw3N+6WL5/KScEDuQktF3iIVJnlhz0UnQ02KNtV3ax/liuLGErm8HhUVFquXlfj86wp6FDl5QsUl+b9fMiWVXx87e5MKAACAwQhpAAAABiKkAQAAGIiQBgAAYCBCGgAAgIEIaQAAAAYipAEAABiIkAYAAGCgvAxpixcv7nDZ8iUrOlxv2ZLlSW+nqivrdEVXx0m0H1r6ycNPnLFs9Ytr2r3dGbH1Vr/4ktYse7ndNosXL9YDDzzQ7vL2brdn6fNLtXHlZknS+hWvJmybyuMbVm6M95fKnBLVYrJlyWxbs01bV7/R7mOvr349pWUdPR673XZZ2z6qH6tOOL9UxLb9/fXvpNS+rbfX/bnVv5i2/bXcx+3VVUfzWrx48RnHtL3H2juGyeo3VS3XObBld6fXb6nl+l2txZrt++K3D25tPZ/Y/YNbd5/xWHck6qvlmJ3po6P2Lbevpc6cexKt2502iaRS120F3v0kLW1b6sy2p7oP2mvf3f1norMupK1YurLD9Za/sCLp7VR1ZZ2u6Oo4ifZDS3/8/fozlq1+8aV2b3dGbL01y17WS8vWtttm8eLFWrp0abvL27vdnqVLXtTG358KaSuThLAUHn/195vi/aUyp7SGtJe26401HYW0rSkt6+jx2O22y9r28ebajoPYtpe2JxwvJrbtf13/bkrt23r7jzta/Ytp21/LfdxeXXU0r+6EtGT1myrTQtrR7fvjt/MxpLXcvpZMD2mp1HVbgfcOp6VtS4S0rsvLkAYAAJDrCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABiKkAQAAGIiQBgAAYCBCGgAAgIEIaQAAAAYipAEAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABiKkAQAAGIiQBgAAYCBCGgAAgIEIaQAAAAYipAEAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABiKkAQAAGIiQBgAAYCB3tieQDjNmzOhw2W13/GOH6/3T1NuS3k5VV9bpiq6Ok2g/tDTx1glnLLvp9hvbvd0ZsfWmfG2yooq222bGjBkaMGBAu8vbu92eO6bdrg+PfyBJmvCP1yVsm8rjUTuiiB1JeU6JajHZsmTG3PgluZ2edh/78k3jU1rW0eOx2+0ta2nspDEJ55eK2LZfPOHSlNq3ddnEUe0ub9tfy318xx13pDyvtsemM8e6vXG6cqxbrtO/alin12+p5fpdrcU+X6qM364Y33o+sfttl3dXov5SHbOjubbVcvta6sy5J9G63WmTSCp13ZbvkvK0tG2pM9ue6j5or31395+JHLZt28kaBQIBlZaW6vjx4/L5fJmYV1qEoyF9Efy7nHLI4Th730S07aiisnWu1bvDF3lJaow0am9gt1wOl5xp2l9RO6qIHdEg3zAVuAp6vP8ToTpt/ew1eRweuZzd/5kkEg0rZIc0vuwqFXtKemCGXVPbWKtle1+U1+WVJ8ExzJZQNKSmSJO+Nuh2+Qv8Sdsfqa/Rr3b+h3q5C2S5rG6PH4wE1RBu1Hev+KbOL+rb7f5McSDwqb634ecqsnqpwO3t9PqN4SbVBxv0i6/OU39fvy7NYW/tId3zwg9U4i1SoafzczDdyVCT6prqVT31pxrkvzDb00mrD2s+1i0L7pOvsFi9rMTn34ZgowInT2jV/3hKQ/qe+cMzUteZTHX2JhUAAACDEdIAAAAMREgDAAAwECENAADAQIQ0AAAAAxHSAAAADERIAwAAMBAhDQAAwEApXd0zdr3bQCCQ1smkWzgaUl2wTg455Mj2ZLLIlmTLlsvyJL2Y7YnACTnkSOvFbG3ZCiigoCvY4/2fCNWpvq5BTjnl6oFtiNhRRRVVoDCgqCfpdaDTJtAY0MkTJxVyhORyurI2j45EohGF7bACgYDcweSnmbr6OjXVNyrqDMvdAxcdDkfDCkXDqgvUqVek5y+SnC11gTqFGppU3xRV0NXU6fVDkZBCkZDqAnUKqKjLcwifDKo+aCvk7vnnbLYFw0GFY/vInduvecmcqKtTpCmk+vAJRTyhhG0bQ02KREI6UVenQEF+75d0i2WpFP6WQGp/ceDQoUOqqKjo/swAAACggwcP6sILE/9Vi5RCWjQa1eHDh1VSUiKHw7z3oAKBgCoqKnTw4MGc/rNV6B7qABJ1gGbUASQz68C2bdXV1am8vFxOZ+JPeFL6XMHpdCZNeybw+XzGHARkD3UAiTpAM+oAknl1UFpamlI7vjgAAABgIEIaAACAgfIipHm9Xs2fP19erzfbU0EWUQeQqAM0ow4g5X4dpPTFAQAAAGRWXryTBgAAkG8IaQAAAAYipAEAABiIkAYAAGCgrIW01157TTfffLPKy8vlcDi0cuXKM9osX75cEydOVO/eveVwOLRr166U+z906JAsy9Ill1zS7uMOh0MOh0NvvPFGq+VNTU3x8TZt2tSJLUJXJKuDUCikefPm6dJLL1VRUZHKy8t111136fDhwwn7veeee+LH2LIsDR48WI899pjC4bAkadOmTXI4HDr33HPV2NjYat3t27fH10XmPPXUU6qsrFRBQYHGjh2rbdu2tXr8N7/5ja655hr5fD45HA4dO3YsaZ/UQe5JVAe1tbV64IEHNGzYMBUWFqp///769re/rePHjyfs85prrokfy4KCAo0YMUILFiyIP15dXS2Hw6GLL774jHWXLl0qh8OhysrKHttGJJfsfPDNb35TF110kQoLC9WnTx/deuut+uCDDxL2mYt1kLWQVl9fr8suu0xPPfVUwjZf+cpX9POf/7zT/VdXV2vq1KkKBAJ68803221TUVGhRYsWtVq2YsUKFRcXd3o8dE2yOmhoaNCOHTv0yCOPaMeOHVq+fLl2796tW265JWnfkyZN0qeffqoPP/xQDz30kB599FH94he/aNWmpKREK1asaLVs4cKF6t+/f9c3Cp22ZMkSPfjgg5o/f7527Nihyy67TDfccINqamribRoaGjRp0iT94Ac/6FTf1EHuSFYHhw8f1uHDh/XLX/5S7733nqqrq7V27VrNmTMnad9z587Vp59+qvfff19Tp07Vfffdp8WLF8cfLyoqUk1NjbZu3dpqPeog81I5H4wePVqLFi3SX//6V61bt062bWvixImKRCIJ+865OrANIMlesWJFh4/v27fPlmTv3Lkzpf6i0ag9aNAge+3atfa8efPsuXPntjvmj370I9vn89kNDQ3x5ddff739yCOP2JLsjRs3dnJL0B3J6iBm27ZttiT7448/7rDN3Xffbd96662tll1//fX2uHHjbNu27Y0bN8ZrYMKECfE2DQ0NdmlpabwGkBljxoyx77vvvvj9SCRil5eX20888cQZbWPH7osvvkjaL3WQWzpTBzEvvPCCbVmWHQqFOmxz9dVX29/5zndaLRsyZIg9ffp027Zte9GiRXZpaal9//3329/4xjfibQ4ePGh7vV77+9//vj1gwICubRQ6rSt18Pbbb9uS7I8++qjDNrlYB3n5O2kbN25UQ0ODJkyYoFmzZun5559XfX39Ge1Gjx6tyspKLVu2TJJ04MABvfbaa7rzzjszPWV0wvHjx+VwOHTOOed0ar3CwkIFg8FWy+68805t2bJFBw4ckCQtW7ZMlZWVGjVqVE9NF0kEg0H9+c9/1oQJE+LLnE6nJkyYcMZPsz2BOjBTV+vg+PHj8vl8crtT+lPUce3VwezZs/XCCy+ooaFBUvMnMpMmTVJZWVmn+kbXdaUO6uvrtWjRIg0cOFAVFRWdGs/0OsjLkLZw4UJNnz5dLpdLl1xyiQYNGqSlS5e223b27Nl65plnJDUfiBtvvFF9+vTJ5HTRCY2NjZo3b55mzJiR8h/LtW1b69ev17p163Tddde1eqxv376aPHmyqqurJUnPPPOMZs+e3dPTRgKff/65IpHIGSfAsrIyHTlypMfGoQ7M1pU6+Pzzz/X444/r3nvvTXmcSCSiZ599Vu+8884ZdXDFFVdo0KBBevHFF2Xbtqqrq6mDDOtMHSxYsEDFxcUqLi7Wyy+/rFdeeUWWZaU0Tq7UQd6FtGPHjmn58uWaNWtWfNmsWbO0cOHCdtvPmjVLW7du1d69e3lCGi4UCmnq1KmybVtPP/100varV69WcXGxCgoKNHnyZE2bNk2PPvroGe1mz56t6upq7d27V1u3btXMmTPTMHtkC3WQnwKBgKZMmaIRI0a0ezzbir2gFxYWau7cufrud7+rb33rW2e0mz17thYtWqTNmzervr5eN954Yxpmj54wc+ZM7dy5U5s3b9bQoUM1derUM74A1Fau1UHn3h/OAb/73e/U2NiosWPHxpfZtq1oNKo9e/Zo6NChrdr37t1bN910k+bMmaPGxkZNnjxZdXV1mZ42kogFtI8//livvvpqSu+iXXvttXr66adlWZbKy8s7/Dhk8uTJuvfeezVnzhzdfPPN6t27d09PHwmcd955crlc+uyzz1ot/+yzz3T++ed3u3/qIDd0pg7q6uo0adKk+Bc+PB5P0v5nzpypH/7whyosLFS/fv3kdLb/HsXMmTP18MMP69FHH9Wdd97Z6Y9R0T2dqYPS0lKVlpZqyJAhGjdunM4991ytWLFCM2bM6LD/XKuDvHsnbeHChXrooYe0a9eu+L+3335bVVVV8Y8125o9e7Y2bdqku+66Sy6XK8MzRjKxgPbhhx9q/fr1Kb94FhUVafDgwerfv3/CJ5jb7dZdd92lTZs28U5qFliWpdGjR2vDhg3xZdFoVBs2bND48eO73T91kBtSrYNAIKCJEyfKsiytWrVKBQUFKfVfWlqqwYMH64ILLujwhVmS/H6/brnlFm3evJk6yIKung9s25Zt22pqakrYf67VQdZ+RDhx4oQ++uij+P19+/Zp165d8vv98a+51tbW6sCBA/FrYu3evVuSdP7557f7E/auXbu0Y8cOPffccxo+fHirx2bMmKHHHntMP/nJT844UU+aNElHjx5N+Xec0HOS1UEoFNLtt9+uHTt2aPXq1YpEIvHfS/D7/Sn//kEyjz/+uL73ve/x7kmWPPjgg7r77rt15ZVXasyYMXryySdVX1+vr3/96/E2R44c0ZEjR+L18u6776qkpET9+/eX3+/vkXlQB9mVrA5iAa2hoUHPPvusAoGAAoGAJKlPnz499kN2dXW1FixYQB1kSbI62Lt3r5YsWaKJEyeqT58+OnTokH72s5+psLCwRz+WNKEOshbS3nrrLV177bXx+w8++KAk6e67747/8u6qVatanaSnT58uSZo/f367v4OwcOFCjRgx4oyAJkm33Xab7r//fr300ktnXGPL4XDovPPO6+4moQuS1cEnn3yiVatWSZIuv/zyVutu3LhR11xzTY/Mw7IsaiCLpk2bpqNHj+pf/uVfdOTIEV1++eVau3Ztq18e/vWvf60f//jH8ftXXXWVJGnRokW65557emQe1EF2JauDHTt2xK97OXjw4Fbr7tu3r8cuNFpYWKjCwsIe6Qudl6wOCgoKtGXLFj355JP64osvVFZWpquuukqvv/66+vbt22PzMKEOHLZt21mdAQAAAM6Qd7+TBgAAkA8IaQAAAAYipAEAABiIkAYAAGAgQhoAAICBCGkAAAAGIqQBAAAYiJAGAABgIEIaAACAgQhpAAAABiKkAQAAGIiQBgAAYKD/Dzy7yCRhNcarAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, (ax_map, ax_barcode) = plt.subplots(2, 1, figsize=(6,6.5),\n", " gridspec_kw={'height_ratios':[10,1]})\n", "\n", "gpd.GeoDataFrame(geometry=[outer_box], crs='EPSG:3857').plot(ax=ax_map, color='#d3d3d3')\n", "city.plot(ax=ax_map, edgecolor='white', linewidth=1, color='#8c8c8c')\n", "\n", "plot_hexagons(stops_gb, ax=ax_map, color='cluster', cmap='Greens', location_id='h3_cell', data_crs='EPSG:3857')\n", "plot_pings(traj, ax=ax_map, s=6, color='black', alpha=0.5, traj_cols=tc)\n", "ax_map.set_axis_off()\n", "\n", "plot_time_barcode(traj['unix_ts'], ax=ax_barcode, set_xlim=True)\n", "plot_stops_barcode(stops_gb, ax=ax_barcode, cmap='Greens', set_xlim=False, timestamp='unix_ts')\n", "\n", "plt.tight_layout(pad=0.1)\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.12.3" } }, "nbformat": 4, "nbformat_minor": 5 }