observation_hours.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Analyse des données Vigilo pour identifier les heures où ont lieu les observations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                               address approved categorie  \\\n",
      "0                        Avenue de Lodève, Montpellier        1         2   \n",
      "1                        Avenue de Lodève, Montpellier        1         4   \n",
      "2                        Avenue de Lodève, Montpellier        1         2   \n",
      "3                       Rue Joseph Cugnot, Montpellier        1         3   \n",
      "4                          Route de Mende, Montpellier        1         2   \n",
      "5                                   golf de castelnau         1         2   \n",
      "6                         Avenue de Nîmes, Montpellier        1         2   \n",
      "7    Avenue de l'École d'Agriculture Gabriel Buchet...        1         2   \n",
      "8                            Rue du Truel, Montpellier        1         4   \n",
      "9             Rue du Faubourg Figuerolles, Montpellier        1         2   \n",
      "10                           Avenue Lepic, Montpellier        1         3   \n",
      "11                    Boulevard Renouvier, Montpellier        1         3   \n",
      "12                         Cours Gambetta, Montpellier        1         6   \n",
      "13                           Rue du Truel, Montpellier        1         3   \n",
      "14                           Rue du Truel, Montpellier        1         4   \n",
      "15                           Rue du Truel, Montpellier        1         4   \n",
      "16      Avenue de la Justice de Castelnau, Montpellier        1         6   \n",
      "17                         Cours Gambetta, Montpellier        1         2   \n",
      "18                         Cours Gambetta, Montpellier        1         2   \n",
      "19                     Avenue des Moulins, Montpellier        1         2   \n",
      "20                    Boulevard Renouvier, Montpellier        1         2   \n",
      "21                       Rue de la Galéra, Montpellier        1         2   \n",
      "22                       Avenue de Lodève, Montpellier        1         2   \n",
      "23                     Rue du Plan d'Agde, Montpellier        1         5   \n",
      "24                       Place Saint-Roch, Montpellier        1         5   \n",
      "25                    Place de la Comédie, Montpellier        1         5   \n",
      "26                             Rue Campan, Montpellier        1         2   \n",
      "27                       Rue Paul Rimbaud, Montpellier        1         2   \n",
      "28                Rue de la Vieille Poste, Montpellier        1         4   \n",
      "29           Rond-Point Benjamin Franklin, Montpellier        1         6   \n",
      "..                                                 ...      ...       ...   \n",
      "570                      Rue de Barcelone, Montpellier        1         2   \n",
      "571                      Avenue de Lodève, Montpellier        1         2   \n",
      "572                         Avenue de Fès, Montpellier        1         3   \n",
      "573                      Rue de Barcelone, Montpellier        1         2   \n",
      "574                          Rue du Truel, Montpellier        1         2   \n",
      "575                      Avenue de Lodève, Montpellier        1         2   \n",
      "576                      Avenue de Lodève, Montpellier        1         2   \n",
      "577                      Avenue de Lodève, Montpellier        1         2   \n",
      "578                      Rue de Barcelone, Montpellier        1         2   \n",
      "579          Rue du Faubourg Saint-Jaumes, Montpellier        1         2   \n",
      "580                       poste rondelet, rue Rondelet        1         5   \n",
      "581                       Place de Phocée, Montpellier        1         2   \n",
      "582                         Rue de Raguse, Montpellier        1         2   \n",
      "583                       Rue Ray Charles, Montpellier        1         3   \n",
      "584                    Placette Le Moulin, Montpellier        1         3   \n",
      "585        Avenue du Martinet, 34830 Clapiers, Francia        1         3   \n",
      "586               Rue de la Vieille Poste, Montpellier        1         3   \n",
      "587                       Avenue du Pirée, Montpellier        1         4   \n",
      "588                       Avenue du Pirée, Montpellier        1         4   \n",
      "589                Avenue Henri Becquerel, Montpellier        1         3   \n",
      "590                    Avenue Jean Mermoz, Montpellier        1         6   \n",
      "591                   Place de la Comédie, Montpellier        1         5   \n",
      "592                    Rue Georges Méliès, Montpellier        1         3   \n",
      "593                  Pont Jean Zuccarelli, Montpellier        1         3   \n",
      "594                       Avenue du Pirée, Montpellier        1         3   \n",
      "595              Rue du Comté de Melgueil, Montpellier        1         3   \n",
      "596               Boulevard de Strasbourg, Montpellier        1         3   \n",
      "597              Boulevard de Strasbourg, Montpellier         1         3   \n",
      "598                      Avenue de Lodève, Montpellier        1         3   \n",
      "599                      Avenue Jean Prat, Montpellier        1         5   \n",
      "\n",
      "                                               comment     coordinates_lat  \\\n",
      "0           Stationnement habituel devant le snack...    43.60912577551626   \n",
      "1            Piste cyclable non réparée après travaux.            43.61044   \n",
      "2                                                                 43.60921   \n",
      "3                                trottoir obligatoire   43.595030651498796   \n",
      "4    Marquage temporaire insuffisant, plots nécessa...   43.63599459547765   \n",
      "5       GCUM comme d'habitude sur la PC devant le golf            43.62383   \n",
      "6                                             gcum ups   43.61562995419824   \n",
      "7                             Ça dépasse ? Vraiment ?    43.61274158756225   \n",
      "8                                                       43.627542181092146   \n",
      "9                                Ah ! c'est ma place !  43.606361395818475   \n",
      "10                                                                43.60084   \n",
      "11                                                                43.60167   \n",
      "12      pas de signalisation sur la chaussée pour velo            43.60646   \n",
      "13    piste cyclable droite termine sur passage piéton  43.629117437641405   \n",
      "14                                                      43.627966395360104   \n",
      "15   haie sur le trottoir -> piétons sur piste cycl...   43.62756693509621   \n",
      "16                                    peinture effacée   43.62690269912045   \n",
      "17                              gène un autre cycliste   43.60600403457251   \n",
      "18                             Gambeta style... Zeeman   43.60595402325053   \n",
      "19                                                      43.620590028190655   \n",
      "20                                    Every morning...  43.605787237618955   \n",
      "21            voiture bloquant trottoir danger pietons  43.637660137338784   \n",
      "22                         poubelle sur piste cyclable   43.61297705083104   \n",
      "23                                                                43.60816   \n",
      "24                                                                43.60843   \n",
      "25                      au lieu des grilles du parking              43.609   \n",
      "26      passage piéton toujours utilisé comme parking             43.60747   \n",
      "27              voiture sur piste cyclable et trottoir   43.61797083971064   \n",
      "28                                     marquage effacé   43.61833298432741   \n",
      "29   piste cyclable traverse rue passante sans marq...   43.61845386053034   \n",
      "..                                                 ...                 ...   \n",
      "570                      Piste toujours obligatoire !     43.6037506980055   \n",
      "571                Camionnettes sur une piste cyclable    43.6099860343091   \n",
      "572  piste cyclable double sense sur trottoir + inv...   43.63306231432475   \n",
      "573            Piste cyclable encore obligatoire !!!!   43.603774170889984   \n",
      "574                                     tous les jours   43.63032400644474   \n",
      "575                                               GCUM   43.60996157128582   \n",
      "576                          Camion sur piste cyclable   43.60993113090111   \n",
      "577                          Camion sur piste cyclable   43.60999971994542   \n",
      "578                      Piste cyclable obligatoire!!            43.603801   \n",
      "579                                          mal garé    43.61438545045207   \n",
      "580           il y a de la place pour quelques arceaux   43.60396428464999   \n",
      "581                                      File parking   43.600311020638046   \n",
      "582                           Parking sur le trottoir    43.60063364932313   \n",
      "583                                      Coupé en deux   43.59892875839891   \n",
      "584                       est-ce vraiment nécessaire ?  43.609900420323775   \n",
      "585                          bande cyclable dangereuse   43.63900151824729   \n",
      "586  Stop qui croise passage piéton et piste cyclable    43.61126831598584   \n",
      "587                                      nids de poule   43.60858126078484   \n",
      "588                                                      43.60488009204526   \n",
      "589                                                     43.618068981440864   \n",
      "590                                                      43.61248582088645   \n",
      "591         Difficile de garer un vélo sur la comédie    43.60855277659723   \n",
      "592                       Trottoir ? piste cyclable ?   43.603388417630384   \n",
      "593                       piste cyclable non connectée   43.60072945843522   \n",
      "594                       piste cyclable non connectée   43.60032157603305   \n",
      "595                       Piste cyclable ou trottoir?   43.603608623587334   \n",
      "596                                                     43.603740049960415   \n",
      "597                       piste cyclable non connectée    43.6038261567388   \n",
      "598                                      la Continuité  43.610655044541424   \n",
      "599  Seulement 10 arceaux non sécurisés pour 450 go...  43.590615337302076   \n",
      "\n",
      "        coordinates_lon  group        time                           token  \n",
      "0    3.8665290176868434      0  1554551040                        K04X8OOZ  \n",
      "1               3.85726      0  1554550560                        DZWQB2UL  \n",
      "2               3.86622      0  1554547380                        1FN1RHVJ  \n",
      "3     3.876265436410904      0  1554500400                        AQ249KG3  \n",
      "4     3.868088722229004      0  1554479470                        E08DBKLC  \n",
      "5               3.91871      0  1554478260                        Z4GPXLO7  \n",
      "6     3.883943259716034      0  1554477420                        I44XLBKL  \n",
      "7    3.8595851138234143      0  1554476880                        VBVBVCL3  \n",
      "8     3.869011402130127      0  1554464520                        RUDVWHTB  \n",
      "9    3.8641616329550743      0  1554463189                        9M6GGBWC  \n",
      "10              3.86767      0  1554461460                        UTDVYJ6O  \n",
      "11              3.86818      0  1554461400                        092SFUKP  \n",
      "12               3.8722      0  1554461280                        N2BE1DF7  \n",
      "13   3.8678030669689174      0  1554453360                        SI75BT6B  \n",
      "14   3.8686265051364903      0  1554453300                        4XUXXEUX  \n",
      "15    3.868977203965187      0  1554453240                        VUZYF8J3  \n",
      "16   3.8699260354042053      0  1554453180                        GNKEBNC7  \n",
      "17    3.872769512236118      0  1554447720                        YQQVS5YL  \n",
      "18   3.8728583604097366      0  1554447240                        EW36RUWM  \n",
      "19   3.8333574682474136      0  1554447060                        EUBY8K2O  \n",
      "20    3.867492601275444      0  1554446580                        Y8QL7KC9  \n",
      "21   3.8373623415827747      0  1554444900                        5UNYMG1Y  \n",
      "22   3.8437459990382195      0  1554443700                        2DJL1EFU  \n",
      "23              3.87591      0  1554402660                        6XHBAWZW  \n",
      "24              3.87574      0  1554402480                        PER935GN  \n",
      "25              3.88023      0  1554401700                        RK0W92UZ  \n",
      "26              3.88296      0  1554398040                        FBEB7A0V  \n",
      "27    3.839509449899197      0  1554397680                        996A0FKW  \n",
      "28   3.9126596599817276      0  1554394560                        CUDBWD2X  \n",
      "29   3.9126137271523476      0  1554394380                        E90PGQD1  \n",
      "..                  ...    ...         ...                             ...  \n",
      "570  3.8878363720566567      0  1547816225  0urUpvj0r7OCYtgyTYHyVe2Bs9dhhZ  \n",
      "571  3.8614859338127427      0  1547813650  PkLhQbkDIeGWpUYmtQ4cEi9zz2Bc2w  \n",
      "572   3.833777904510498      0  1547807400                        9XDNQTI1  \n",
      "573   3.887851017098436      0  1547747043  xmqS5TAAL0FAu7F0gkHdsEltpofjVJ  \n",
      "574  3.8667565276177784      0  1547713500                        PL0OUTR2  \n",
      "575  3.8614299273099073      0  1547710560                        EULQUIVK  \n",
      "576  3.8620907077295286      0  1547652699  n1W704hR5hQU1tMwidc6vT4EaMAmmT  \n",
      "577  3.8615543961457055      0  1547651846  gZNBPKDK5ZQi75vdrz2LJh0nQSbSll  \n",
      "578            3.887844      0  1547643193  chNij7UxBA4pBOMTpyXhRKm5ev0udh  \n",
      "579  3.8701996868441135      0  1547635860                        NMFD7PBG  \n",
      "580   3.875076947699778      0  1547546400                        82YC5HZM  \n",
      "581  3.8988551326987686      0  1547298517  ynD9Ozt9uNB4Q3Zbmp1smRKwfhbrrF  \n",
      "582  3.9004433155059948      0  1547298289  ttueacArvsHCA81zNHD4tN8KvVj63l  \n",
      "583    3.90309290376607      0  1547294320  rOVRRvqt2pghAy9PReqWjdqRYIKgcQ  \n",
      "584   3.896498811879505      0  1546784280                        7XWGGKKT  \n",
      "585   3.847742805512162      0  1546707000                        KHP7HYEM  \n",
      "586   3.900763918444552      0  1544801820                        A6KL9NK4  \n",
      "587   3.896885812282549      0  1544773500                        2M42KQ4H  \n",
      "588   3.896570101961383      0  1544773440                        16S026U6  \n",
      "589  3.9132821559905873      0  1544687520                        J4MSEZNQ  \n",
      "590   3.889216726476503      0  1544614800                        258M9SU9  \n",
      "591  3.8792932033538685      0  1544366340                        XYKWVDQQ  \n",
      "592  3.9141734631966285      0  1543519260                        ZCYBC2SQ  \n",
      "593  3.8969062844047064      0  1542043620                        A0T6SG5Y  \n",
      "594   3.896162509918226      0  1542043560                        TW7NMG5B  \n",
      "595  3.8908052444457875      0  1541934720                        0DLKNYCX  \n",
      "596  3.8887864351272716      0  1541840700                        O6J81W9Y  \n",
      "597   3.889809250831604      0  1541840640                        Q18A635V  \n",
      "598  3.8494745841036515      0  1541336400                        WDHP90GZ  \n",
      "599   3.852594271302223      0  1537100520                        OPZJJTVI  \n",
      "\n",
      "[600 rows x 9 columns]\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "\n",
    "import requests\n",
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "\n",
    "# Getting data\n",
    "resp = requests.get('https://vigilo.jesuisundesdeux.org/get_issues.php')\n",
    "df = pd.DataFrame(resp.json())\n",
    "\n",
    "print(df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Separating week days and weekends\n",
    "df.time = pd.to_datetime(df.time, unit='s')\n",
    "df['day_type'] = 'week'\n",
    "df.loc[df.time.apply(lambda x: x.weekday() > 4), 'day_type'] = 'weekend'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3Xd4VFX6wPHvmcwkkx5SSAUSEhJI6B1pBkVBUUQEsaJYFgW77rqroqL7W7Gt6+raO0sQUYoiukoXkRKkJIRQExLSQ3oymXZ+f9yIgHRmMpPkfJ5nnpA7d849d4B557T3CCkliqIoiuJudK6ugKIoiqKcjApQiqIoiltSAUpRFEVxSypAKYqiKG5JBShFURTFLakApSiKorglFaAUpZkJIZYLIaa6uh6K4u6EWgeltCVCiGHAi0AKYAOygAellJtdWjFFUf5A7+oKKEpzEUIEAN8A9wALAE9gONDoynopinJyqotPaUsSAaSUaVJKm5SyQUr5PynlDgAhxDQhRJYQokII8b0QotNvLxRCSCHEvUKIvUKIGiHEc0KIeCHEBiFEtRBigRDCs+ncdkKIb4QQpU1lfSOEiDmmrNVCiDub/nybEOInIcTLTeceFEKMPebcQCHEB0KIQiHEYSHE80IIj+Z6wxTFlVSAUtqSPYBNCPGJEGKsEKLdb08IIa4B/gZcC4QB64C0E14/BugHDAb+DLwL3AR0ALoDNzSdpwM+AjoBHYEG4I3T1GsQkA2EonU/fiCEEE3PfQJYgQSgD3AZcOe53riitEQqQClthpSyGhgGSOA9oFQIsVQIEQ78CfiHlDJLSmkF/g/ofWwrCpgjpayWUmYCGcD/pJQHpJRVwHK0AIKUslxK+aWUsl5KWQP8HRh5mqrlSinfk1La0AJSJBDeVK+xaGNkdVLKEuCfwBSHvSmK4sbUGJTSpkgps4DbAIQQXYG5wGtorZ1/CSFeOeZ0AUQDuU2/Fx/zXMNJfo9oKtcHLZCMAX5rpfkLITyagtCJio6pX31T48kPCAYMQOHvDSp0QN5Z37CitGCqBaW0WVLK3cDHaN1zecCfpJRBxzy8pZQ/n0fRjwBJwCApZQAwoum4OPVLTioPbQJH6DF1CpBSppxHnRSlxVEBSmkzhBBdhRCP/DZhQQjRAW3c6BfgbeCvQoiUpucChRCTzvNS/mgtqkohRDDw9PkUIqUsBP4HvCKECBBC6JomZpyuu1BRWg0VoJS2pAZtQsJGIUQdWmDKAB6RUi4C5gDzhRDVTcfHnrKk03sN8AbKmq7x3QXU+Va06fC7gApgIdoYlaK0emqhrqIoiuKWVAtKURRFcUsqQCmKoihuSQUoRVEUxS2pAKUoiqK4pXNaqBsaGipjY2OdVBVFURSlLUhPTy+TUoad6bxzClCxsbFs2bLl/GulKIqitHlCiNwzn6W6+BRFURQ3pQKUoiiK4pZUstg2SkrJ7qIaSmoaqaw3U22yEhFgpGdMIOEBRldXT1EURQWotmZ3UTVLthXw9bYC8isbTnpOmJ8XqV3DuGNYZ5Ii/Ju5horSclksFvLz8zGZTK6uilswGo3ExMRgMBjO6/UqQLURJTUmnvsmi6+3F6AT0CM6kCt7RhIV5I2vlx5vgwelNY0cLKtlX0kti7cVsGBLPiMTw5iRmsDAuGBX34KiuL38/Hz8/f2JjY3lmC1S2iQpJeXl5eTn5xMXF3deZagA1crZ7ZL5m/N4YXkW9WYbE/tGc1lyBAHef/xGE+zrebTFVGOy8GNWCf/LLGLyOxu4ZXAnHh/bFV8v9U9GUU7FZDKp4NRECEFISAilpaXnXYb6tGnFTBYbD32+jeUZRSRHBnDHsDiigrzP6rX+RgMT+kRzRY8IFmzOY+4vuazKLuHlSb0Y3DnEyTVXlJZLBaffXeh7oWbxtVJVDRZu/WATyzOKuGlQR568sttZB6djeek9uGVILLPGJWO12bnxvV/4dEOOw+urKIpyIhWgWqGiKhOT3v6ZrYcquG9UAuN6Rl3wN5mukQH849qe9OnYjllLMnn260xsdrVVi6K4o+LiYm688UY6d+5Mv379GDJkCIsWLQJg06ZNjBgxgqSkJLp27cqdd95JfX09zzzzDC+//PJx5cTGxlJWVgZoraFbbrnl6HNWq5WwsDDGjRvntPtQAaqVqaq3cNP7v5B3pIG/jOnKRfGhDivbaPDg4UsTGds9go/W53D3Z1toMNscVr6iKBdOSsk111zDiBEjOHDgAOnp6cyfP5/8/HyKi4uZNGkSc+bMITs7m6ysLMaMGUNNTc0Zy/X19SUjI4OGBm327w8//EB0dLRT70UFqFak0Wrj7s+2kFtez2OXJ9E9OtDh19DpBLcOieX2i2JZmVXCnZ9uVkFKUdzIypUr8fT0ZPr06UePderUifvuu48333yTqVOnMmTIEEBrFV133XWEh4efVdljx45l2bJlAKSlpXHDDTc4/gaOoSZJtBJSSv6ycAcbDx5hZmoC3SIDnHq9y1IiMBo8eHvNfu78dDPv3zoAb08Pp15TUVqSZ7/OZFdBtUPLTI4K4OmrUk57TmZmJn379j3pcxkZGUydOvW8rz9lyhRmz57NuHHj2LFjB9OmTWPdunXnXd6ZqBZUK/HPH/eyeFsB1w/owNAEx3Xrnc6IxDCmj4zn533l3PnpZkwW1ZJSFHczY8YMevXqxYABA0573qnGqY893rNnT3JyckhLS+OKK65waD1PRrWgWoF1e0t5fcVeRiaGMb5XVLNee0SiljH/7TX7mTFvK+/c3A+9h/reoyhnauk4S0pKCl9++eXR3998803Kysro378/Y8aMIT09nfHjx//hdSEhIRQWFh53rKamhqCgoOOOXX311Tz66KOsXr2a8vJy59xEE/VJ0sKV1Tby0OfbiGnnze1DXbNAcERiGLcNjWVFVgl/W7QTKdXsPkVxlVGjRmEymXjrrbeOHquvrwdg5syZfPLJJ2zcuPHoc3PnzqWoqIgRI0awdOnSoxMmvvrqK3r16oWHx/Fd99OmTWPWrFn06NHD6feiWlAtmN0ueWTBdqobrDx2eVe89K4bA7osOYKqBgsLtuQT6ufFn8d0dVldFKUtE0KwePFiHnroIV588UXCwsLw9fVlzpw5hIeHM3/+fB599FFKSkrQ6XSMGDGCa6+9loiICGbOnMmwYcMQQtC+fXvef//9P5QfExPDAw880Dz3ci7fdvv37y/VhoXu4/11B3h+WRbThsYyOjnC1dVBSskHPx1kxe4Snr4qmduHnl/+LUVpqbKysujWrZurq+FWTvaeCCHSpZT9z/Ra1YJqofaV1DDnu90MiG3Hpd3OboqoswkhmDY0jmqThdlf7yIqyJvLU1wfOBVFaZnUGFQLZLdL/vLlTox6D+4Y1tmtcn/pdIIZqQkktPfj/rRf+fVQhaurpChKC6VaUC3QfzcdIj23gukj4wk8SVbyZiUlxuqD+JdtxcNcg5A2hLTzSmIQc7YZ+NPHv/DFjBF0CvF1bT0VRWlxVIBqYQqrGnhheRY9ogMZ0aV51jv9gZQEFP9C+J55BBb9gqHx5FNN/wuYpZ7sNxOoT70Vnz6TwK9989ZVUZQWSwWoFkRKyVOLM7DaJHcMi2v+rj0pCSpYS/TONwgoTcfiGUhdSA/qgrpR3y4Rm2cgEh0IHXpTOcaaXOpKctAX7cRnxd+QK59ExI+CITOgcyq4UdekoijuRwWoFuTHrBJ+zCrhpkEdCQ8wNuu19aZyEtY/RruC1ZiNIRR0vZ3KqJFID8+Tnm/2i8bsFw2RF7EuDO7fms/DIb8wpmAN4rMJENUHhj0MXceBTg2FKoryR+qToYUwWWzM/iaTmHbejOnevDPjAgp/ptc3VxBYtJ7CxFvYN/SfVHQYfcrgdKLhUTC8Wwz3lF3HS9FvwJD7oLYYFtwC746EfT+CWtyrKA61aNEihBDs3r37jOe+9tprRxfzuhMVoFqID346SN6RBm4Z3Al9c7U4pCR6xxsk/3gLdp2BAwNnc6TTWKTu3Bve13aGK2PhPzsln1lSYfxbMPwRqC2BuRPh06vhcLrj70FR2qi0tDSGDRvG/Pnzz3ju6QKUzea6HJsqQLUARVUm3li1jwGx7egZE3TmFziClHTa+g86bn+VqsihHBj0PI3+nc67OCHgTykwKByeXm/ix0N2bRzqmrdg4J+gcAe8NwoW3Apl+xx4I4rS9tTW1rJ+/Xo++OCDowFq9erVx20uOHPmTD7++GNef/11CgoKSE1NJTU1FQA/Pz9mzZrFoEGD2LBhA+np6YwcOZJ+/fpx+eWXH83Z9/rrr5OcnEzPnj2ZMmWKw+9DjUG1AHO+243VZuemQecfIM6JtBO3+Vkisj+jvMNlFCXdCuLCv8t46ODPfeHxDXDfigbmX6WjV3sDdLsKEi6BzEWQuRiyvoG+t8LFj4O/WuirtFDLH4einY4tM6IHjH3hjKctXryYMWPGkJiYSHBwMFu3bj3luffffz+vvvoqq1atIjRUmxlcV1dH9+7dmT17NhaLhZEjR7JkyRLCwsL4/PPPeeKJJ/jwww954YUXOHjwIF5eXlRWVjrsNn+jWlBubuuhChb9epgre0Q1z8QIKem88Skisj+jrNOVFCVNdUhw+o1RD7MGQIAnTFtez4HKpu4Dgw/0vgmufReSxsKvn8G/esOK2WCqctj1ncbaCA2VYLO4uiaKQlpa2tEWzZQpU0hLSzun13t4eDBx4kQAsrOzycjIYPTo0fTu3Zvnn3+e/Px8QNt+46abbmLu3Lno9Y5v76gWlBuTUvL3ZVkE+RgY37t5ttGI2flvwvemURo3npL4yU6ZCh5shNmD4C/rJTd9U88X432J8W8Kgt7tYNB06DYetv8X1r0CWz7UxqsG3AWG5p29eJSUUHkI8jdD6W4o36c9qvKhsRbsxwQmvRGMQRCaAOHdITwFOg2FkHjX1F1xjbNo6ThDeXk5K1euJCMjAyEENpsNIQRXX301drv96Hkmk+mUZRiNxqNZzKWUpKSksGHDhj+ct2zZMtauXcvSpUt57rnnyMzMdGigUgHKjX2fWUx6bgV3Do/DaHB+pvKQg1/TYftrVESOcFpw+k2MH8weDH/bILnpmzq+GO9Le59jWmoBkTD8UUieAL9+Cv97En55C1KfgJ6TwcPJGTTsdijZBTnr4OA6yN8IdWXac8ID/MPBPwo6DtFaf54+oPMEawOY68FcowWv9I/B2vRBENIFul4ByeMhup9z66+0WQsXLuTWW2/lnXfeOXps5MiRAOzatYvGxkZMJhMrVqxg2LBhAPj7+1NTU3O0i+9YSUlJlJaWsmHDBoYMGYLFYmHPnj1069aNvLw8UlNTGTZsGPPmzaO2tvYP+0ddCBWg3JTFZmfOd7uJbufNxYnOz77gV/orCT8/Rl1QVwqT72iWRbTxgfDMQHjyF8nN39STdpUPId4ndCeGxMOlz0LhNtj6CSy5F358GnrfCH2nOq5V8ltAyl0POT9B7k9Qf0R7zj8SInpCaFcIS4SgTmcfIKUdqguhYCvkbYQNb8D6f0FUXxh4N6RMcF2rUGmV0tLSePzxx487NnHiRObNm8fkyZPp2bMnXbp0oU+fPkefv/vuuxk7diyRkZGsWrXquNd6enqycOFC7r//fqqqqrBarTz44IMkJiZy8803U1VVhZSShx56yKHBCdR2G27rsw05PLUkk8cuS6Jvp3ZOvZZnXQE9vh2PFHoODnwWm2eAU693ou1l8Owm6OCvY944H9r7nmLMS0o4vAX2fAf5W0DatMCRcAnEXwIdBoLe6+wu2lgDJbshfxPkrIdDP0NDU2Jbv3CtWy6ip/ZwZHomcy0cWAO7l0HVIfANhWGPQP9pKlC1Amq7jT+6kO02VIByQ7WNVka8uIqIAC+evDLZuSmN7FZSfrgJ3/KdHBj0HGbfaOdd6zR2lsGzm6G9j2DeOF+i/c8wMaO+HPavgsNNY0J2m9b1FtQRQhO1nwZvbTxIpwdTpdZFV1faNHaU93tZ/lFNAak7hPdonnyBUkLRdtj5BRRuh4BobdZirxvBQ3VstFQqQP2R2g+qlXl37QGO1Jl5eHSi0/Ptxez8NwElm8nvPsNlwQmgRyg8PxhmbZRMWlrHp1f4kNDuNONuPiHQ4zrtYa6Hoh1QvheqDkPZHq1FZDWDzQxIbZzIGAhe/hDcWVuDFdRJ6yL0dUHSXSEgsrf2KNyudV8uvQ82vQtXvqq1BhWljVMBys2U1zby/roDDIoLJj7Mz6nXCijaQMyON6iIGklV5FCnXutsdG0H/zcYnt4kuXZxHW+N9mFozFn8E/X0gY6DtceJpNTGgXTOn2Ry3iJ7wRWvaONemz+AD0ZDn5vh0tngG+Lq2innSErpVnu0udK59NCdjFoH5Wb+s3o/JouNSf07OPU6elM5XX56CLNvJIVdpzr1WuciIQheHQbBXnDrt/XM22W+sAKFcO/g9BshIHa4llkj5VrYngZv9ocdC1SewhbEaDRSXl5+wR/MrYGUkvLycozG8x9bVS0oN1JQ2cBnG3IZ3iWM6CBv511ISjr/8iT6xgoODJyN9HCvwflwH3hpKLywFf62zsTOMhuzhhjxNrSBb6UGb23CRPwo2PBv+OouLUiN+ycEOfdLi3LhYmJiyM/Pp7S01NVVcQtGo5GYmJjzfr0KUG7k9RV7sUvJxL7n/xd6NkJyviEk73uKutxwQfn1nMnHAE8PgE+zIS3LwuZCG69f6k1ySAtoDTlCu1gY8yLs/kbLqvGfQXDJMzDgTrU9iRszGAzExcW5uhqthgpQbuJAaS1fbMlndHI4Yf5nOVX6PBgaSonb9DT1gQmUd7rSaddxBA8d3N4NeofCq7/aGf9VHQ8P8OKOHp54erSB1pTOQ1vU23EwbHgTlj8GGQvh6n9DWJKra3d+TNXaOrPiDG0yS9lebUamtGtdmR4GCO2izcQMS4KESyHQuV/YFPelppm7iZnztvLjrmL+eX1vgnzObp+lcyYliWvupd3hlewf9H/ahoItRFUjvLEDfi6C+CAdzw41MuxsJlC0FlLCgZXaJAqrCUY8BkMfBL2T/q04UnWBNqa29wctVZTdqh33ba9Nr/dpB+i0cTirSTu/Kh9sjdp50QMgZTz0nAJ+YS67DcVx1DqoFiSzoIorX/+Ja3pHc/0A540zhOR8Q+K6+ynqcgPlsVc57TrOtLkY3smEwjq4pJMH9/c10qt9G+n2A20x8ab3IGcttE+Gq9+AGDdMmyQlHFwDm96H7G+1RdXB8RDdV8uiEdpFW6N2ytfbofowHNqgLaQ+sl8bnxtwF1x0vwpULZwKUC3I7R9tYnNOBa9d3xtfL+e0Cjwaq+iz5BIsXu04OPBZh2Yob25mGyw6AF/th1oLjIjx4J4+XgyO9Gg703vzNsLGt7R0TIPugVFPgKevq2ulZXTfnqYF0SP7tXVnCaMhcayWX/F8VebBzgVa0NN7absyj3j07DOHKG5FBagWYnPOESa9vYEbBnTg6t7O63KL2/gk4Xvmc2Dw85j8Y512neZUb4FlubBoP1SZoVOAYHJXTyYmGog4Vbqk1sRcpy3wzf5Wy5xx1b+02X+uULANNr+vZcawmrTxo6QrIXYYeDiwG7IqH7bN01qQYV3hmv+oxLstkApQLYCUksnvbGBvSS2vXd8bL71zuqr8SrfR/buJHOk4hqKkW5xyDVcyWWF9IfyQBzvLQQC92+sYHWvg0k56urTTte6WVVGGloS2Ol+bVJH6RPNMorA0aJtMbn4fDqdrXXZxIyHpCudvLZK/RbvnhiPaWFzqEypFVAuiAlQLsDq7hNs+2sztQ2O5LNlJO8farfT8djye9cXsu+gl7Honrq9yAwV1sOYwbCyGvU0bfIZ6CwZFejAoUs+gKA+6tNOha20By2bWWi+7FmubJ/aYDCP/7JxAUb4f0j+CX+dqY2KBHbSgFJ8Knmef/URKicUOBh3n9wXCXAdbPoC9/9MC46SPwSf43MtRmp0KUG7Obpdc9cZPlNY08sqkXug9nNMlFZH1EXFbniOv5wNUhw9yyjXcVVkDbCnRWlWZR6C0QTse5AUDI/Va0IrS0y1Yh4eulQQsUxVkfKmtn7KZtZyD/adpuxRfyB5adeWQ+ZW2aDh/kzYFvsNgrRsvosdpt2epNEl+LrCSWWZjV7mdPUdsVJsl9RawSTB6QHtfQbiPjq7BOgZHaV8kQk/ceuVU9v2oTcMPiIIb0rTEv4pbUwHKzS3dXsD9ab9y78XxDO/inBlJhvpi+iy5lPrABA71+XOz7PHkrqSE4gbIKP/9UVivPRfgCQMiPRgcqWdwlJ7kkFYQsBoqtG1J9n6vZXH3CdG2JEm4RAta/uGnf725Xktie2C19ji8RZseHtQJOl+sjXX5nDpP4JEGO9/st/LdQQsbC23YJHgI6OAPsf7ae270AC8PbaLLkUbtC8WBamhomoXeI0zHDV09uTrBgJ/nGf4+SnfD6v8Diwmu/0xbP6W4LRWg3JjFZueSV9YgBPzfhB5O627qsnYmwXk/sm/IHCw+Z/hAaoPKGrTW1c6mgHW4Tjse6AmpHQ1cGqtnZAc9/mf6cHRndps2XpO7TpvIYGrq9zQGQrs4aNcJ9N7aNHBph9oSbTuSmsKmAoQ2JTyytzbhoV3cab/o7Cq38fFOM4v3WTDboIMfDImAQRHQOQA8zzDMarXDvirYUQZrCiCnGnz0MCHRwMw+XkT6naZVVV8OK57VZvxNfB9Srjm390ppNipAubG5v+Ty5OIMHrs8ib4dnbMZYWDBWpJX3EZx/CTKOk9wyjVam3KTFqy2lsDmEqg2a+Mjl3bSc00XAxd31OPVkjNYSDscOaBNqqgphNoiqC1uWjjbtFDWy1/rKguI1mYGhvcArzOPK+0otfHSJhPr8m0YPWBUDFwZC7EXsPellJBdCd/lwurDoBNwW3dP7untRZDxFH8P5lpYMVtrUV31L+h76/lXQHEaFaDcVIPZxogXVxHi58mscc7ZjFBnNdHr6zEIu5X9Q/6B1F3A2EMbZZOw+wj8VAhrC6CyURu7mpTkyS0pnnQMaAPT2M/Cvgobr2xuZPlBKwGecG08jOkI/g5OcFFcD3OzYVW+VvaTQ4xMSjKc/P+PxaR19xVshbEvwqA/ObYyygVTAcpNvblqHy99n80zV6WQFOHvlGvEbP8nHXb8m5x+T1AXrAaML5TNDr+WwQ+HtFRLUkJqRw/u6tXGFgcfo9YseT29kQ93mjF4wITO2sPHyd+FcqrhPzu1SS8jO3jwjxHeRJ2s289mgTVzIO8XGPca9L/duRVTzokKUG6ost7M8BdXkdjen0cvd846FWP1AXp9PZbq9gM43GOmU67RlpU1wHeHYHmu1qrqF65jRh8vUjvq20SgklKydL+Vv28wUVIvuawjTO2qtS6bi13Cshz4OEtLKPzCCG+uSjhJZLRZYNXzcHgrTHgbek1pvkoqp3W2AUr1UzSjt9bsp9ZkZbKz8u1JSdzGp5E6A8WJNzvnGm1cqDfcnAQfXgL39IC8ajvTvmtgwuI6NhZYXV09pzpcY+f25fU8sKKBQE/Jq8PggV7NG5xAG4u6Kg7evBg6+sF9Kxp4Zr0Js+2EL9seBrj4bxDZExbfoy0qVloUFaCaSVGViY/X5zAsIZSOwT5OuUZIztcEFa2nOOF6rF5BTrmGovHygHGx8N4ouL8X5Ffbuf7reu74ro69FTZXV8+h7FLyWaaZ0Qtq+aXAxt0p8MowSHLO/J6zFuEDL1wE13SGjzPMTF5aR2Gt/fiT9F6Q+pSWFumru+HgOtdUVjkvKkA1k3+t2IPNLrmun3P2tvEwVxO75XnqAzpTEXOJU66h/JFeB5d3hHdSta6uDYdtXP5FHY+vaaC4zn7mAtzcgUob1y+t56mfTCQGaa2W8Z21NU3uQK+Du1Lgr/0g+4id8YvqyCg74QuCwQijngK/CJh/IxTvck1llXOmAlQz2F9ay4LN+VzSLZz2Ac7ZXr3DtlcxmI5Q2G1ai85U3lIZ9TC5C7w/Cq6KhYXZFkam1fLKZhN1lrMf53UXVrvk7W2NjF1YR1a5jQd7wfODtVaLOxoWBS8PbcpvuaSOFbmW40/w8odLnwGdHuZO1PacUtye+iRrBq/+bw8GvWBCH+dkK/ct30HEnrkc6TAaU0Bnp1xDOTuBXnB3d3g7FQaGw7+3mkmdX8uXe8zYz2FCkittK7Fx1Vd1vLCxkX5h8NbFMLqj+yciiQ2AV4dBtC/c9X0Dn2Wajz/Br70WpEwVMPc6LS2U4tZUgHKyXw9VsGxnIVd0jyTQ2wlzcO02Ov/yJFbPAEriJzm+fOW8RPrCX/pp3+rbeUoeWWVi/KI6thS570SKGrPkmfUmJiyqo6TOzt/6wd/6Q7BzGv1OEWzUxqUGtIenfjLx762NHDdTObizNnGibDd8fgtYzacuTHE5FaCcSErJc9/sIsjHwFW9opxyjfC9/8XvSAZFiTdjN7hp/0sb1i0YXh4Gj/SBwho71y2p574f68mvcZ/xKSkl3x2wcOnntXySYebKWHj7Yhga5f6tppMx6uGJ/nBJDLyyuZHnNjQe33qN6qPtyntwDSyZoS1sU9yS2kDFib7ZUcjWQ5XcPbwzRoPj93oyNJTS8deXqQ3uTnX4EIeXrziGTmipfy6KgIX74ct9Vr7PqeW2FE/u7XOatD3NIL/GztPrG1iRa6NzgHvMznMEDx082Bt8DfDhTjM1ZskLI4y/JwGOv0RLovvrZxAYrXX9KW5HBSgnMVlsvLB8N51CfBiZ6Jxs5Z22/B2dzURh19tb5lfdNsao19ZQXdZBS9vz3g4zabvNTO/lxdTunmfO2O1AlSbJm7828nGGGQ8B05Lhmjjtg7210Am4OwX8DDAv24JNwksjjwlSPSZDXSn89E9tT6sBd7i2wsofqADlJB+tz+FwZQNPXNENnRO2bggoXE9YzlJKOl+L2TfS4eUrztPeBx7uo+Wt+ywbXtrcyDvbG5na3ZPbunsScrb7IJ2HSpNk7i4z725vpMYMl3TQgmZYK93HUgi4KUkLVnOzLUgJL1/cFKSEgEH3QH0ZfPuoFqQSL3N1lZVjqFTdTimZAAAgAElEQVRHTlBa08jFL6+ia3iAU1Ia6awmen4zFp2tkf2DX0B6ODgzp9Ks9lTCgr2woUjbI2l8FwPXdzXQp73j8vzlVtn5KMPMgt1m6q0wIBxu63ph2cZbmvl7tC8E4xP0vJLqjf63L46WBvjucS3D+7TvILKXayvaBpxtqiPVgnKCF5bvptFi56ZBHZ1SfsyO1/CuyeVgvydUcGoFEoPgyQFwqAYWHYAley18vttClyAdV3cxcHEHPSmh575NfUm9nW8PWFm6z8zWYjseAi6OhgnxENeGAtNvpiRqjaZPd1uRNPDqb0HK4A2XzNJaUfMmw50rINA5C+qVc6NaUA62JecI1729gat7RXHDQMcHKN/ynfRYPoHKqJEUJN/l8PIV16u3wrrD8L882F2hHQv1FgyK9KBLOw8Sg3V09NfhbQCjh0AnoLxBUtpgJ79Gsq3ExtZiKwertP/bcQEwIkqbqBHaSrvyzsWCvfDJbhgXr+e1Uce0pCpyYPlftE0cp30PxjYYxZuJakG5gNVm56klGYT6eTplUa6wW4jf8BesnoEUdbnR4eUr7sFHD5d30h4VjdoGiltKJFuLrXx7wMqZvlIGeULXdjAyCgaFQ0fn7OrSYk3uoo1JfZRlRcoGXhvljcFDQLtYuPivsOJp+GIq3LhASziruIwKUA70342HyCqs4cFLujhlWnlk5nv4VuzmUK+HsBt8HV6+4n7aeWkTGS5pSoBvsmpb05fUQ6MdGm3a9hNBXtq5IUYINapJnWdyXYL2Hn24ywocE6SiesPgGfDz67DsYbjqdfVmupAKUA5SVtvIy//Lpkd0IAPjgh1evndFNh12/Iuq9gOpaT/A4eUrLYNRD/GB2kO5MBPjQQAfnBikulwGtcWw9VOtVTX8ERfXtO1SAcpBnv16Fw1mG1MvinX4xnXCZqbL+oew6X20ZLCKojjEtfHazz8Eqd43a0FqxWwI6gQ9rnNpPdsqFaAc4MddxXy9vYBJ/WKIDnL8KHTMjtebuvYeweapBm4VxZGubWpJvX9ikLroAS3bxOJ7ICAKOl3k6qq2Oa1o3bhrVJssPLF4Jx2DfbjaCfn2/Eq3Ep35NhVRF1PTvp/Dy1cURZt6f2cyLDtg5YGVDVhsUpsgkfqElgV9/o1QttfV1WxzVIC6QC8s301pTSN3j+iM3sF5YnSWOhLWP4LFK4SiJLWFu6I4029B6ttjg5SXP1zyjJZQdu5EqClydTXbFBWgLsCG/eXM23iIsd0jiQ/zc2zhUtJ541MYa/I4nDIdu15lKlcUZztpkPKP0Bby1hbDp9dA/RFXV7PNUAHqPFU1WHh4wTYiAoxM6u/4Vedh+78g7OBiSjtfS31wN4eXryjKyZ00SIUmwqgnoXyflm3CXOfqarYJKkCdByklTy7OoLjaxIzUBLz0jl3z5F2RTdymZ6gNTqG08zUOLVtRlDM7aZCK7A0jHoPD6fD5zWBtdHU1Wz0VoM7D4m2H+Xp7Adf160BCe8d27eks9SSunYndw8jh7jNAqL8iRXGFY4PUvT82YLJKbSbfkJmwfyUsuFXtyOtk6tPvHOUdqefJxRl0jfBnvKNn7Uk7CT8/hnf1AQ53vxerV5Bjy1cU5ZxMiIfp3eGHHCvTltdTZ5HaQt5B98Ke77SUSDaLq6vZaqkAdQ4arTZmztuKlHDvxQkO3+cpZse/CTm0nOIuN1AX0t2hZSuKcn6uioOHe8PGQhs3fVNHpUlC1ytg4J8g+1tYeLsKUk6iAtQ5eHpJJtvzq5g+Ip4wfy+Hlh2c+y0ddvyLisgRlHe60qFlK4pyYS7pAH/tB5lldiYtqaOg1g7droIBd0HW17DgFrCYXF3NVkcFqLM0b+Mh5m/O45reUQxwcK493/IMEtY/Sn1gFwqT71DJKRXFDQ2JhGcHweFaO9csqiOr3AbJ47VdebOXw3+vg8YaV1ezVVEB6ixsPVTBrCUZ9IoJZFK/Dg4t21idQ9eV07AZfDnU62GkTqX3VxR31SsUXhoKdrtk0pI61udboeuVMOwRyP0ZPrlarZNyIBWgziDvSD1/+jSdYF9PZqZ2cei4k2d9Ed1+vAWdzUxu38exeakU1Yri7mID4OVh2rYmU7+t57+7zBCfChf/DYoz4IPRcOSAq6vZKqgAdRpH6sxM/XAT9RYrj16WhJ/Rcbl19Y0VdPvxVgyNR8jt+2fMvo7f4FBRFOcI84YXh0LvMHhinYlZPzVgiR4Io2dDbQm8dwnkbnB1NVs8FaBOocFs446PN5NXUc+jlyXRIdhxqYa04HQbxppcDvV6BFNAZ4eVrShK8/A1wKyBMKEzfJppYeq39ZQHJMMVL4HBGz69GrZ/7upqtmgqQJ2E2WpnxrytbM+v5L7ULnSNcNwWF4b6ElK+n4JP5W7yej5IfXCyw8pWFKV5eQi4MwUe7A2bi2xcsbCOLfXhcMXLEJoEi+6Gbx9TC3rPkwpQJzBZbEz/bAsrd5dw20VxDp2x51WbT/fvJ+NVm0dunz9TG9bHYWUriuI6ozvAK8NAJyTXL63nvd2e2EfPhuRrYNO78NFYqMxzdTVbHBWgjtFgtnHnJ1tYmV3KHcPiGJ0c7rCyfcszSPluEvrGI+T2+xv1wSkOK1tRFNeLD4R/DYdB4fD3Xxq59TszRcnT4OK/QskueGc4ZC52dTVbFBWgmlSbLNz20SbW7ytj+sjOXNrNccEp5ODXpHw/GSGt5PR7kobABIeVrSiK+/A1wN/6w8wesKXQxmULalliGQhX/hN8w7TUSF/eqaainyUVoICcsjomvLmeLbkVzEhNYGRie8cUbLfR4deXSPzpAUz+nTgw6Hka/Ts6pmxFUdySEDA2Fl4fAVG+8MCKBu76pR0Fw+dA75sg8yt4awjsWqJthKickpDn8Ab1799fbtmyxYnVaX4/7y/jnrlbsdklD13aheQox6xF8qrNJ2H9IwSUbOZI9CiKut6G1DlumrqiKO7PZofFB+C/e7QJFQ8P8OK2yEPoN7wOFQchfhSMfQlC21avihAiXUrZ/4zntdUAZbNL3l93gJe+zyY80MhjlyURHmC88IKlJOzAV8RtegakncKuU6mKHK7SFylKG1ZcD2/thM0l0KWdjicG6hlZ/z1i2zywmWHQn2DYw+Ab4uqqNgsVoE7jcGUDj3y+jV8OHmFgbDB/GtkZH88Lb90Yqw8Su/k52hWspq5dNw6nTMfiHXbhFVYUpcWTEjYUwUdZUFAHQ6M9eLJ3A91y5sKBlWDw0faaGjIDjI5b2uKOVIA6Cbtd8kV6Hs8vy8Jqs3PrkFhGJoYhLrB1o7PUEbPzDSKzPkAKA6XxEynvOEZtNqgoyh9Y7LA8B+bthRozpHb04LHEEpLz0rR8fl4B0P92GDQdAhy855ybUAHqBOm5FTyzNJOdh6voGuHP9JHxF9ylp7PUEbFnLlGZ72FoPEJF1AhKEqaojQYVRTmjOgt8kwOL9kONBYZEeXB/bB6DyhejO/Sz9gU35Vrod5u2k28rGiZQAapJVmE1b6zax7IdhQT7enLjwI5cFB9yQa0mQ0Mp7fd+TmTWhxjMldSE9KQ0/jo1fVxRlHPWYIVvc2BpDpQ1QIyfYHpCBddal+OTuwIs9RDcGfrcAt0nQrtOrq7yBWvTAUpKydZDFby1ej8/ZpXgbfDg8pQIxveOwmjwOM9C7QQUbyR8TxrBed+hs1upCe1DaecJKjApinLBbHb4pQi+zoGd5aATMCrKwj1BW+hVvRJ9aaZ2YlQfLUNF0hUQ2qVFtqzaZIA6Umfmq635zN+cx76SWvy89IzpHsHlKRH4eZ3HJAi7Df+yrYTkLic4dzleDcVYDX5URg6nIuYSzL6ts39YURTXOlwLK/NhVT4UN4BeB1eHl3GT7xa6123Aq2KvdmJADHS5FOJGat2A/hGurfhZajMBKu9IPT9mFfO/zGI2HTyCTUq6tPcjNak9Q+JDzq3FJCVetYcILPqFwMJ1BBaux2Cuwq4zUBvSi+rwQVS3H4D08HTeDSmKojSxS9hdobWsNhRps/8A+vmVc1PgTgbJ7URW70RnrdeeaBcLHQZDRHcI7w4RPcA31GX1P5VWGaAsNjv7SmrZmV/FxoNH2HignPzKBgA6tPOmX6d2DIkPpePZbI0h7XjV5uFTkY1PZTZ+ZdvxL/sVQ2OFdi2vYGqDu1Mb0pPasN7Y9Y7bbkNRFOVcSQn5tbCtDLaVwo5yqLeCHitDjLmM8d1DP7LpZN6Ht+WYVEp+4VqgCu8OIfEQ1EkbxwqIAQ/XJA9osQHKarNTWttIQaWJ/Ip6csrqOVhWy96SWvYU12CxafUNMOrpGhFAt0h/enUIIjLQ+7hydJZ6PBtKMDSU4NlQiqGhBK/6IrxqcjE2PTxspqPnm3yjaQiMpyGwC/VBSTT6RrfIvl1FUdoGm4RDNVoLa3cF5FRDXi002iCYarrqDjHQ6xB9DYeIl4eIsObjIa1HXy+FBwREI9p1gsAO4BcGvu3Br72WN9Cvvfa7MRD0ju01ckmAKqoyUVRtotFio9FiQ9aWYLZasVgsmK02rMf8uc5kprbBTF1jI7UNFmpNjdSZzNSZzOilFS9hwRMLnlgJNdoJ8xZE+EJ7b0GoNwQabBgsNXhYatCba/Aw1+BhqUZvrsFgKsfDWveH+tl1esze4Vi822P2icDkG0OjfwcafWOw6x2QRUJRFMWF7FLLWnGoBnJrtIBVUAflJqhssBEuyukgSukgSo7+jNWVEinKCaYKA7aTlmvReWHR+2M1+GP1DMDu6Y/d0x+pN/7+MHgjPbzBYAS9j/bTYER4eCI89Og8DKDzQOr0RPa+/KwClEPbdx+uP8i7aw8A4IWZbONt517IyQK1HahrejSRCOwGH+wGX6TBV/uzpy9W3/Y0evXF5h2Mzag9rN4h2IzB2D39/7B4VgAqNCmK0lr4tYP4kxy32SVHTKGU1nWhtN5Oab2drWY7axolNWZJjcmGNNdiMFfiY60iwFZJoKzEx15PgKjH31xPgKjDnwYCRAF+mDAKM16Y8caMETN6YXfovZxTC0oIUQNkO7QGbVsoUObqSrQi6v10LPV+OpZ6P3/XSUp5xjxw59qCyj6bZplydoQQW9T76Tjq/XQs9X46lno/z51KFqcoiqK4JRWgFEVRFLd0rgHqXafUou1S76djqffTsdT76Vjq/TxH5zRJQlEURVGai+riUxRFUdySClCKoiiKW1IBSlEURXFLKkApiqIobkkFKEVRFMUtqQClKA4ihLhYCJHfTNeSQgi1lbPSqqkApbR6Qoi/CiG+PeHY3lMcm9K8tVMU5VRUgFLagrXAUCGEB4AQIgIwAH1POJbQdK6iKG5ABSilLdiMFpB6N/0+AliFlpn/2GP7pZQFQoiuQogfhBBHhBDZQojJvxUkhPASQrwshDgkhCgWQrwthDh+t8zfz71fCLFLCBHT9Ps4IcQ2IUSlEOJnIUTPY87NEUI8KoTYIYSoEkJ8LoQwHvP8Y0KIQiFEgRBimiPfHEVxVypAKa2elNIMbEQLQjT9XAf8dMKxtUIIX+AHYB7QHrgB+I8QIqXpvDlAIlpgSwCigVknXlMI8RRwGzBSSpkvhOgLfAj8CQgB3gGWCiG8jnnZZGAMEAf0bHo9QogxwKPAaKALcOl5vxmK0oKoAKW0FWv4PRgNRwtQ6044tgYYB+RIKT+SUlqllFuBL4HrhBACuAt4SEp5REpZA/wfcOy4lRBCvApcDqRKKUubjt8FvCOl3CiltEkpPwEagcHHvPZ1KWWBlPII8DW/t+4mAx9JKTOklHXAMw55RxTFzTl0R11FcWNrgRlCiHZAmJRyrxCiGPik6Vj3pnMmAYOEEJXHvFYPfAaEAT5AuharAG1TZo9jzg0C7gaul1JWHXO8EzBVCHHfMcc8gahjfi865s/1xzwXBaQf81zu2d2yorRsKkApbcUGIBAteKwHkFJWCyEKmo4VSCkPCiHygDVSytEnFiCE0AENQIqU8vAprlMB3AwsEEJMkFKubzqeB/xdSvn386h7IdDhmN87nkcZitLiqC4+pU2QUjYAW4CH0br2fvNT07HfZu99AyQKIW4RQhiaHgOEEN2klHbgPeCfQoj2AEKIaCHE5SdcazVwE7BICDGo6fB7wHQhxCCh8RVCXCmE8D+L6i8AbhNCJAshfICnz+c9UJSWRgUopS1Zgzbx4adjjq1rOrYWoGlc6TK0caUCtG63OcBvkxn+AuwDfhFCVAM/AkknXkhK+QNwO9pEiH5Syi1o41BvoLWy9tE0CeJMpJTLgdeAlU2vW3m2N6woLZnaD0pRFEVxS6oFpSiKorglFaAURVEUt6QClKIoiuKWVIBSFEVR3JIKUIqiKIpbOqeFuqGhoTI2NtZJVVEURVHagvT09DIpZdiZzjunABUbG8uWLVvOv1aKoihKmyeEOKt0XaqLT1EURXFLKhef0mI1Wm1kF9VQ22jFZLFhtkpiQ32ID/PD4KG+eylKS6cClNKilNU2svjXw6zbW8amg0dosNj+cI6nh47ECD8uS45gyoAOtA8wnqQkRXE8i8VCfn4+JpPJ1VVxC0ajkZiYGAwGw3m9/pxSHfXv31+qMSjFFUprGnl37X7m/nKIBouNqCAj3aMCSY4MwN+ox1PvgU5AYZWJnPI69pbUkl1Ug14nuCwlnOkj4+kZE+Tq21BauYMHD+Lv709ISAjHbMnSJkkpKS8vp6amhri4uOOeE0KkSyn7n6kM1YJS3JrNLvnPqn28uWofZpudofGhjO8TTXTQSXdZp3OYH0MTQgEorGpgRVYJa/eUsnxnETcP7sSjlycR6H1+3+YU5UxMJhOxsbFtPjgBCCEICQmhtLT0zCefggpQitvKr6jnwfnb2JJbwaC4YK4f0IHIwJMHppOJDPTm5sGduLZvNF+k5/Pfjbkszyhk9vjuXNEj0ok1V9oyFZx+d6HvhRpJVtzS95lFjP3XOjILqrn34ngevDTxnILTsXw89UwdEsvz1/Qg0NvAvf/dyv99m4XVZndwrRVFcSQVoBS389+NuUz/LJ32/l7849oeDO9yxvV8ZyUu1JdnrkrhsuRw3l17gKkfbaKizuyQshXFnRQXF3PjjTfSuXNn+vXrx5AhQ1i0aBEAmzZtYsSIESQlJdG1a1fuvPNO6uvreeaZZ3j55ZePKyc2NpaysjJAaw3dcsstR5+zWq2EhYUxbtw4p92HClCK25BS8uaqfTyxKIPeHYN4alwy4Q6egaf30HH70DjuHt6ZjQeOMP7N9eRX1Dv0GoriSlJKrrnmGkaMGMGBAwdIT09n/vz55OfnU1xczKRJk5gzZw7Z2dlkZWUxZswYampqzliur68vGRkZNDQ0APDDDz8QHR3t1HtRAUpxC1JK/rF8Ny99n82whFAeHp2Il97DaddL7dqep8YlU17byOS3N5BTVue0aylKc1q5ciWenp5Mnz796LFOnTpx33338eabbzJ16lSGDBkCaK2i6667jvDw8LMqe+zYsSxbtgyAtLQ0brjhBsffwDHUJAnFLfxn9X7eXXuAy5LDmXpRLLpmGGhODPfniSuT+cfyLCa9s4G0uwaT0N7P6ddV2oZnv85kV0G1Q8tMjgrg6atSTntOZmYmffv2PelzGRkZTJ069byvP2XKFGbPns24cePYsWMH06ZNY926dedd3pmoFpTicgs25/HS99kMTQhttuD0m7hQX566MhmL1c7172zgoGpJKa3MjBkz6NWrFwMGDDjteaeacXfs8Z49e5KTk0NaWhpXXHGFQ+t5MqoFpbjUiqxi/vrVTnrGBDJ9ROdmDU6/6RDsw1Pjknn260xu+WAjX91zkco+oVywM7V0nCUlJYUvv/zy6O9vvvkmZWVl9O/fnzFjxpCens748eP/8LqQkBAKCwuPO1ZTU0NQ0PEL3K+++moeffRRVq9eTXl5uXNuoolqQSkus7uompnzfiU21IeHLk1E78L8eVFB3vx5TFfKahuZ+tEmqk0Wl9VFUS7EqFGjMJlMvPXWW0eP1ddrE4FmzpzJJ598wsaNG48+N3fuXIqKihgxYgRLly49OmHiq6++olevXnh4HD8WPG3aNGbNmkWPHj2cfi8qQCkuUdVg4e5P0zEadDx6WRJGg/MmRJyt+DA/Hro0kT3Ftdz1yRYarX/M86co7k4IweLFi1mzZg1xcXEMHDiQqVOnMmfOHMLDw5k/fz6PPvooSUlJdOvWjXXr1hEQEEDPnj2ZOXMmw4YNo3fv3rz99tu8//77fyg/JiaGBx54oHnuReXiU5qb3S6589MtrN1TylPjkkkM93d1lY6zfl8Zb6zax3X9Ynjpup4qM4By1rKysujWrZurq+FWTvaeqFx8itv698p9rNxdwu1DY90uOAEMTQilsKqBhen5JIX7c9eIzq6ukqK0SaqLT2lWP+8r47Uf9zC8Syiju53d2gtXuLZvDAPjgvnH8ixWZZe4ujqK0iapAKU0m4o6Mw8t2EZkkJE7hsW5ddeZTgjuGRlPx2Af7pv3K/tLa11dJUVpc1SAUpqFlJK/LtpJea2ZmaldnJolwlGMBg8euSwJnYB75qbTYFaTJhSlOakApTSLL7bk811GEZP7dyAu1NfV1TlroX5ezEhNYG9xLU8tyXB1dRSlTVEBSnG63PI6nv46k+5RAVzZs+Xtw9QzJogJfaJZmJ7Pgi15rq6OorQZKkApTmW3Sx77Ygc6YPrIeJdkinCEiX1jSIkK4KnFGWQVOja/mqI4w6JFixBCsHv37jOe+9prrx1dzOtOVIBSnGruxlw25RzhpsGdCPHzcnV1zptOJ5iZmoCPpwf3zE2nRmWaUNxcWloaw4YNY/78+Wc893QBymZz3dirClCK0+Qdqecfy3fTKyaQixMds+mgKwX5eHLfqC4cOlLP41/u5FwWuStKc6qtrWX9+vV88MEHRwPU6tWrj9tccObMmXz88ce8/vrrFBQUkJqaSmpqKgB+fn7MmjWLQYMGsWHDBtLT0xk5ciT9+vXj8ssvP5qz7/XXXyc5OZmePXsyZcoUh9+HWqirOIWUkse/3AES7hze2a2nlJ+LbpEBXN+/A2mb8xj0SzC3Dol1dZUUd7X8cSja6dgyI3rA2BfOeNrixYsZM2YMiYmJBAcHs3Xr1lOee//99/Pqq6+yatUqQkNDAairq6N79+7Mnj0bi8XCyJEjWbJkCWFhYXz++ec88cQTfPjhh7zwwgscPHgQLy8vKisrHXabv1EtKMUpFmzJY/3+cm4c1JHQFty1dzLjekXRt2MQs7/exfY8x/+nVJQLlZaWdrRFM2XKFNLS0s7p9R4eHkycOBGA7OxsMjIyGD16NL179+b5558nPz8f0LbfuOmmm5g7dy56vePbO6oFpThcaU0jf1+WRXKkP5d0be/q6jictog3gb8u2sGMeVtZdv9wAr0Nrq6W4m7OoqXjDOXl5axcuZKMjAyEENhsNoQQXH311djt9qPnmUymU5ZhNBqPZjGXUpKSksKGDRv+cN6yZctYu3YtS5cu5bnnniMzM9OhgUq1oBSHe37ZLurNNqYNaz1deyfyM+q5f1QXCqtM/PWrHWo8SnEbCxcu5NZbbyU3N5ecnBzy8vKIi4sDYNeuXTQ2NlJVVcWKFSuOvsbf3//oNhsnSkpKorS09GiAslgsZGZmYrfbycvLIzU1lRdffJHKykpqax2bcUUFKMWh1uwpZcm2Asb3jiI6yNvV1XGqLuH+XN+/A9/uLGLuxkOuro6iAFr33oQJE447NnHiRObNm8fkyZOPdsv16dPn6PN33303Y8eOPTpJ4lienp4sXLiQv/zlL/Tq1YvevXvz888/Y7PZuPnmm+nRowd9+vThoYce+sPmhhdKbbehOEyD2cbof67BLiUvXNsTgws3IGwudil56ftsdhVUs2jGRaREBbq6SooLqe02/uhCttto/Z8gSrP598q95Fc0cMewzm0iOMHvSWX9jHpmzvuV2karq6ukKK1G2/gUUZxuX0kt7649wPAuoSRHBri6Os0qwNvAjNQEcsvreHKRWh+lKI6iZvEpF0xKydNLMvDU67hxYEdXV8clkiMDmNg3hi/S87koIZTJ/Tu4ukqOVVsKJZlQvAvKsrXfG45AXRlIGxh8wdMHvAIgNBHCkiCsK0T2AoPR1bVvVlLKVjs56Fxd6Jc1FaCUC7ZsZyHr95dz+0WxBPl4uro6LnNN72iyCquZtSSDPh2C6OKGuwWftboyOLAaDq6FA2ugMuf357wCwCdY++nXHoQH2BrB2gh1ByBnHVibpjAbvKFzKiReDoljwd99N6l0BKPRSHl5OSEhIW0+SEkpKS8vx2g8/y8oapKEckFqG62Menk1vl56nh/fHZ2ubf+nrKg389evdtLe34ulM4fh7en++14dVV0AWV/DriVwaANIO3j6QngPCE+Gdp2hXSfwbnf6cqQdakugIgcKfoXDW6C2GHR6SL4GBk2HmP7QCj/ALRYL+fn5p11j1JYYjUZiYmIwGI5fJ3i2kyRUC0q5IP/6cQ8lNY08NyqhzQcngHY+ntx7cTwvLN/NM0szmXNdT1dX6fTMdbBrKWyfBwfXARKCOkHP6yFmAATHg+4cg6zQgX+E9ug4GKSEylzY9wNkL4eMhRDVFy55CuJHOeW2/r+9O4+Purr3P/46s89k3yE7ayBsAQTZERXcrntttUWt2mrr0nq1drG3XltttT9sbbXVn2213rZexQ2FqnW3LqgQwh4gEEzIvk4ymX3m+z33j2/AjSUhE7Kd5+Mxj+9kthy+TOY933PO93MGitVqPXTOkdJ3KqCU47ansYvH3q9iWVEm4zP73p1linixBZqxBtsxaUFM0RBCD4Ewo5usSJMVzRpP1J5MxJ6KZkswPgwHmem5yZxfksPq0hrmj0vjgpk5A92kL6srg9LHYOfzRkgljIYZl0LhEkiO8Qprf4AAAB7WSURBVPiZEJBSCHO+DSUrofIt4/f+/UKYsAKW3wWZk2L7O5VhQQWUclyklPzshR04bWYundu7DzRLoJWEljJcHXtwdezF1bEHu68Oc7R369HoJgth12hCcTmE4nMJJI3DnzwJf0oRYWfWgHYhfWV2LrsbPdy+ZjvTc5MYmxE/YG05JOQ1jl5KH4OGrWBxQOFiGH8aZE45MfvL6oRJ5xjBtGsdbF8NDy+Aed+FU3824iZUKEenxqCU47Jmcy3/uXor31o8htMmHX3g2xQNktTwHsn1/yax8WNcnspD94WcWYTicoi4sojYk4naU4haE5BmO7rJhjRbQeoIPYpJj2DSgpjDXZgjXizhTqzBVqzBVuPIK/Rp4daoLRFf8iT8KcbFm16CP2lC77ur+qDNG+Ina7aTm+JkzfULcVgHaDyqcYcRSttWQ9hrHM1MPAvGnmKMMQ2kYCds/jtU/MuY/XfhI5Aza2DbpPS7no5BqYBSeq0zEOG037xDssvGz8+bcthVcoUWJqX2TdKqXyal7i3M0QCa2YE/uQh/ymR8KZMIxeejW2L3jdkc8WL31uDoOmBsvTXYvTWYNWPAOmqNx5s+k66MWXRlzMKbXmJ0E/ajzQfc/L9X93DxrFzuu2T6iZvZFfTAjueg7O9QvwnMNihcZARTxqTBN0GhrgzWPwABNyz9ESy5DUyDr/tWiQ01SULpN/e/XkGbN8wty4u+FE4OTxWZ+1aTue9ZrKE2IrYkOrMW4Mmagy+l2JjJ1U80azz+lMn4Uz5TVkXq2PxNODv34eqswNmxl6SG9xFIJAJ/chFdWXPwZM7FkzmHiCu21ddn5qdw8awcniurpSQ/mcvnFcT09T9HSjjwkXFEsvN5iAQgudAY+xl3KtgH8bT3nFlw3h/g44fhnV8ZM/8u+jM4Y1vbTRla1BGU0is76jo57w/vc9rkLK5e+OlspfjWreRs/yOptW8ghYmu9Fm4c0/FmzZ90E1kMEX9ODsru8e/duPs3HfoKCuQUEBX5lw8WUZgheLz+ny0oUvJfa/uYXtdJ6uvm8/sgmNM0+6tjgOw43kjmNr2GeM8hUuMcZ70iYPvaOlopDRm+m18BJLz4dInIVPVthtuVBefEnO6LvnK/1/PvhYvv7mkhHi7hYSmDeRue4DkxvVErfG0563AnXMaUUeMP4T7kx7F0VVNnHs3ro7duDoqsESMpQdCziy6usPKkzWXQNL44wpcbyjKf72wHSlh3U2LyErsY9dm+yew52XYsQbqNhq3ZRYboVSwaOhPNmjaCf++F7QwfOVxmLhioFukxJAKKCXmnt5Yww+f28Z3lo7jzKxO8jf9mtS6N4nYk2nLPxt37mnolmGwxIbUsfvqcLl3d4fWHqyhdgAitmS6MufgTZ+OL2Uy/tTiHs8YrG7zcee6nUzITODp6+b37iTegBtqNsL+t6HiVWjvnmiSNt4IpMJFxnlHw4mvFd6+2wjj8x6Emd8Y6BYpMaICSompDn+YZfe9w/j4MA9m/ZOsfavRzQ5aC8+jLf9MpHkYlziSEmugmbiO3bjce4xp8f6GQ3dH7Mn4UorxpxbjS55EKCGfYHw+EWfGl4KrtLqd375WwdnTR/OHy2Z+edKElEYVhpZd0LzbqH9XswFadhv3m20wahrknGRcEkf3979+YEX88M49RkWK0+6ARbcMrS5L5bBUQCkx9dPntxLd9A9+4VqNLeqlPfd0WsZeiGYbWZXLDzJFA9i9B3B2VePoqsbeVY3DW4NJjxx6jG62E3akE3WkEbGnoNkS0c02qtwRdjcHmJXjYmqG1ahb52uFrnroajLq2h1kT+guvjrZ6MLLmGicvzSSaBFY/3ujNuDJ34Uz71EhNcSpWXxKzJRv28gFW65njrUCX1wRNZOvNiYPjGC6xUkguYhActFnboxi8zdhCzRjC7ZgDbRgCXdiDntwdB3ApAUQepSZMspkq4a30Yq30068024UXk0ZYxwVxaUbEwSS8o26dyP9w9hsNY6c7InGLD8tDGffp6ahjwAqoJQj06JoHzzA+Ld+hd9k55NJ1+LPXTLoZuUNGiYL4fgcwvHHLm0U0eDODbCjDf4838myfOsxnzOiCZMxXd5shdJHQY/Af/xehdQwp/53lcNrqYDHzsD81s95UyvhpUmr8OedosIpRqxm+OlJUJAI178eYHOTWon3mISAWd+EaV+Dsr/B2htB1we6VUo/Up82yudJaVQfeGQxWsseboneyN9Sb6YkL2mgWzbsuKzw85Mh2QZXveJnn1sb6CYNfkLAzJUw4zLY8gS8cpvxnlWGJRVQyqeCHnj+27D2RmT6RG5NuI9XWMB100b4GEg/SrHDXfOMP8RL16mQ6hEhYMbXYcpFsPEv8OYvBrpFSj9RAaUY6srgkcVG/baZl/PyuDt4oS6BlUWQMQxObRrMRsfBr+aDpksuXednrwqpYxMCZl8FE8+E938L798/0C1S+oEKqJFOSvjwIXh0hVHp+ox78BR9lTvXRxifBOcWDnQDR4a8BLjnYEit9VPRrkLqmIQwpp0XLoE37oSNjw50i5QYUwE1kvnb4cmvwas/gZzZcO6DkDWFVRuCtAUlN04Hs3qHnDAHQ0pKyVde9PFxvZo4cUwmMyy+xVj996VbYdszA90iJYbUx89I1bAVHlkC+96CudfBsp+CPYGypij/KI9wbiFMUIWkT7i8BLhvESTaYOVLfl6qjBz7SSOdyQJLfwyjpsKa62D3ywPdIiVGVECNRFufgkeXGxUMzvo1TD4XhCCiSW5/N0iaA1aqFbgHTJYLVi00viDc+EaAh7eE6E3FlxHJYjdW5E0bB898E/b/e6BbpMSACqiRRIvAKz8yvmWmT4Rz7je23R7ZGmZ3u853poJLncI9oBJtcPc8WJQNv/44xLWvBegMqZA6KqsLTrvTKJr71GVG/T5lSFO1+EYKbzM8fSUcWA/F58Psqz+3/Pmedo1znvOxYBT8aPYAtlP5HClh7SfwaDlkxwseWu5iWsYALR0fI5ouqfNKqjp1qj06NR6d1qCkLaDTHpQEoxDWJBEdrCZwWAQOM6Q4TGS6BBkuQV6CiXHJJsanmEhxfOF7tr8NXvkh6FG45nXjqEoZVFSxWOVTtaWweqUxKWLBTTD2lM/dHdUlF6zxUevReegUSLIPRCOVo9nVDveWgTsI355h4+bZdhyWwX9+mi4l+9w6ZU0a21o0yts09rTrBD4z/8NqMs4HS7JDohVsZuM2qwmiEsIahDToCkN7CNwh0D/zsZXpEpRkmpmZZWZ2lpmSTDM2bz3864dgT4JrXhv+Vd+HGBVQimHT4/DybeBMhWW3Q+rYLz3kj5tDrNoQ4sezYXH2iW+i0jNdYeNI6vUayE8U/HKxk8W5g6svVpeS8jad9XVR1tdpbGqK0hU27ou3wthEGJNolHjKiTPOAUu1964eriahJQA1XVDjhcpOqOiAep9xv8sC87MtXJhaxdn7/htTaiFc9YpReFcZFFRAjXTRkNHNselxyJ4FS24zlm74gorurr25WXD7Md8uymCwtRUe3AYNPliQbebWOXZmjxqYoJJSUtmhs75OY319lA/ro3R2rxaSFw9TUmFyKkxOgey4/i3M7gnDzjYoa4EtrUZgLTRt53HbKlqTpiIvX0N2emr/NUDpMRVQI5mnHlZfDnWlMO0SKFn5ufGmg0Ka0bXX0KXzx1MgWXXtDRkhDV6ugmcroSMES3LNrJxiY1meBau5/1JASsknnTqljdqhUGr2G58hWU6Yng4z0o1t2gAvW1XrhfUNYKn7iJ+GH+QtfSaP593NBbMLOGfa6N6taKzElAqokap6PTx9hVEVYuHNULDwiA+956Mgj2wNc8ccOHmYrRY+UgSjsK4K1u43xmcynIKLJlo5vcBCSaa5z2EV0iQ7WjRKGzVKmzQ2NWq0B43PjBR7dyClGdtRrsG7dJWp8nUm7/8r60yncpP/GhKdVi6dk8/KkwvIT3MNdPNGHBVQI42UsOFP8OrtEJ9lnHibnH/Eh39YH+Xr6/yckQ83zTiB7VT6habDxmZjfGpDkzGJIM4K80ZbmJJuYlyKmXHJJjKcgkSbwGEBIQSaLglr4A5JGn06jT7JAY9ORbvOHrfGXrdOuLvqUnac0VVXnGpc8uIHbyAdTkbls2Tuf55tBVfxy/DX2FjVjpSwbFIm31xQyKLx6ZhMQ+gfNISpgBpJwn546RbY+iTkzYVFt4It7ogP7wxJznzGiwnJA0vAMbjG2ZU+6grDtjbY0mJs673wxVWTDk4AjB7hzz/dAfkJUJBgjCEVpxpHTEOalIze/VdSa9+gavbt7Cy4gjd2NfHW7mY6AxHGpMdx9aIxXDI7F4dVdf/1JxVQI0VbJTx9OTSVG2vkzLj0qIsKSim54Y0Ar34SZdVCKFITm4a9sGZMGKj1GhMJfFHwd1dQsnZP6Y63QJrTCKZMlzHjbliSOrnbHySp6WP2LvwNrWMvJKLpfLS/jVd3NlLZ4iMjwc61i8fy9ZPzibOrb2/9QQXUSLD7JaMqBAIW32oUfD2Gv+0Ic8cHQa6aDF8Z3/9NVJTBRugR8jevIs69i92nPEJH7qmA8eVtZ72HF7bUsbPeQ5LTyjWLxnDl/EKSXMM1sQeGCqjhTIvCW3fBB7+D9Amw9CcQn3nMp21v0bjoBR8l6XDHXFDd7cpIZYoGKNz0S+y+Onad+lc8o+Z97v69TV28sKWOsgMdxNnNXDm/kGuXjCXZZRugFg8vKqCGK28zPHsVVL0PE8+CudeC+djf7jpDknOe8xKMSh5cYtR6U5SRzBz2UFh6N9ZQG+Wn/wNvxswvPaa6zceLW+r5aH8b8Q4L1y4ey1WLxhCvuv76RAXUcFT9ITxzJQTcMP8GGHdaj56mS8m1rwZ4+0CUexcYA96KooAl5Kaw9C7MER/ly5/Alzb1sI870O7nmdIaSqvdpLis3LBsPCvnFajJFMeppwGlqpkPBboG7/waHj/bmABx9n09DieA35WGeKM6yreKVTgpymdF7SlUz7od3Wxn8htX4HLvPuzj8lNd3LqiiLvOn0puiou7X9rF0lVv88TH1US0L86RVGJFHUENdp218Py3jRNwxy4zlri29fzEwpf3R7j+9QDL8+D7M4bWeSuKcqJY/U0UbroboWuUL38Cf+rkoz6+vL6T1aU1VDR5yUtxcsuKiZw3IwezGtjtEdXFNxyUr4W1N4EWMoJp3Km9evquNmNSRGH3UuKqN0JRjswIqV8i9Ajly/+BP3XKUR8vpWRzTQfPlNZQ1eZnfGY8t51RxIriLIT6JnhUKqCGsrDfqAix6a/GgoKLfwCJvSsz3ujTuXCNj4gmuX8xpA5wXTRFGQoOhpRJD1N++t/xpU075nN0KdnwSTvPlNZQ3xlkZl4yPzprEvPGpp2AFg9NKqCGqoZt8Ny3oHUPTL3YKPTag1l6n+UJSS5Z66PGo3PvAhiX1E9tVZRhyBpoobD0bsxRP7tP/QuerJN79DxNl/y7ooXnympp94VZOjGD284oYmqO+gP8IhVQQ000DO+ugvd/C/ZEWPSfkP3laa/HEtIk33zZz8YGjZ+fDCUZ/dBWRRnmLME2CsruxRZooWLJH3Dnnd7j54ajOq+VN/Lilnq8oSjnzsjm1uUTKUw/cvmxkUYF1FBSVwYvXg/Nu4xxpjnfPuzaTcei6ZLvvxXgn5VRfjATluX2Q1sVZYQwhz3kb16Fs6uKffN/Teu4i3r1fF8oyj+31fPKjkY0XXLp3Dy+d+oEMhNVf7sKqKEgEoR37oH1Dxirfc6/EXLnHNdLabrktneCPL83wtXFcPG4GLdVUUYgUzRA3tb7iW/fQc3071M7/Xu9ngrr9odZs7mOt3Y3YzULrl44huuWjiPJOXLLJ6mAGuxqNsIL34W2vTB+Ocy5Bmzxx/VSupT86N9BntkT4fIiuHRijNuqKCOY0CNkl/+F5Ib3aBlzPpXz70Wae1/avckT5JlNNXywr40kp5XrTxnHlQsKR+TJviqgBit/O7z5C2Mp9rh046ipB0Vej0TTJf/1XpAnd0e4bCKsLIpdUxVF6SYl6Z+8SFbl03gyTmLP0oeIOtOP66Wq2nys3ljDlpoOshLt3Hz6RC6ZnYvFPHLqJqiAGmx0zZg2/uZdEOqCSecYM/R6cdLtF4U0yS1vBXhpf5SvjocrJqkTcRWlPyU2fkjOzkeI2FOoWPrQYev39VR5g4fVGw9Q0eRlTHocP1hRxFlTR42IRRNVQA0m1R/CK7dB43YYNd0o8JpS2KeX9IYl173m54M6jasnw8Vq6QxFOSEcnirytv0OS6idqpP+i6aJlx/3N0MpJZsOuHl6Yw017gBFoxK4+bQJnDFleAeVCqjBwNMAr98B25+GuAw46WooWNTnw5xmn841r/opb9W5aQYsz4tRexVF6RFTxEvujodJaN1Ma8E57D/5bjT78Z/vpOuSD/e38fzmWuo7gkwalcDNp09gRfHwDCoVUAPJ3w4f/B42PGKs3TTlIph2CVj7Pr10S7PGta/68YQkP5wFJ4+KQXsVRek9qZNetY7MymcJO9KpXLiKztGL+vSSui5Zv7+NNWW11HcGmTw6ge+fNpEVxVnDKqhUQA2EoAc+ehg+fBBCXhizxBhnShwdk5d/riLMT94NkmKHn82BMYkxeVlFUfrA4dlPzo6HcPjqaZj0TQ6U3Ipu7dtJuZouWV/ZyvNldTR6gozPjOc7S8dx3oxsbJahP5lCBdSJFPbDxj/D+7+DQDvkz4eSb/R5nOkgf0Tyi/VBntodYUY6/GgWJPV+lquiKP1EaCGy9j5FWs2rhFyj+WTunbjzlvf5dbXurr91W+s50O5ndJKDby0ey6Vz8ogbwosmqoA6EUJeKPsbfHC/sdJtzmzjiCl9Qsx+xfYWje+96aeqU3LJeGMa+QiajaooQ4qzYw/Zux7D4a2hPfd0qk76GaGEvg8SSynZUtPBum317GroIslp5coFhVwxv4D0+KH3bVUFVH/ythjjSxv+DMEOyJoGM78BWYdfjfN4BKOSh7eE+OPmMMl2uLUEph/faReKopxIepS0A6+QWfkcoNM48XLqpt1A1BGb1UIrmrpYt7We0mo3VrPg3OnZXLGgkJK85Ji8/omgAqo/NGyFj/8E258BLQz584yK4xmTYvpr1tdFuf3dAFUeySk58J2pkGCL6a9QFKWfWYLtZFY+S3L9u2hWF/XF19I46Qo0W2wGj+s6Arxe3sS7FS0EIhrTcpK4ckEh/zF99KCvTqECKlYiAWPhwNLHoOYjsDiMgq6Tz4Wk2M7vrurUWbUhyEv7o4yOgxumwUxVjVxRhjS7t5bMfU+R2FKGZomjsWglDZOvJuKMzR93IKzx3r4WXi9votYdINll5asn5XHxrFyKRvW+6PSJoAKqL6SEhi2w+QnYthpCHmPBwIlnw4TTj7tm3pG0+HUeLAvxv+URLCa4cBxcMh7sg/tLkKIoveDwVJFetZbEpo+RJittBWfTWLQSb/rMmJSAkVJS3uDhtZ1NbDrgRtMlU7ITuWhWLufNyCYjYfCMVamAOh5tlbD9WePE2rZ9YLZBwQKYsMIYXxKxnZ1Q1anzp60hnq2IoOlwRj5cNlGtfqsow5nN10DagX+R1Pg+5mgAX8okmsd9lbbCc2J2VOUJRFhf2cb7+1qobPFhFoKlRRlcMDOHZUUZJDgGtpK6CqiekNJYg2nXOti1Fpp2GLdnTYOxS42qD/bYHi3pUvJercaTu8K8VhXFLODUXGN5jJzY/ipFUQYxUzRAUuN6UmrfxNlVhRQmOkctoK3gHNy5p8YsrGrdft7b28r7+1pp94WxmgULxqWzvDiL5cVZZA3A+lQqoI7E1wr734H9b0PlO+CpBQRkFhuTHgoXGWWJYmyfW+OflVGe2ROmzitJtMGKfDh/jDpiUpSRzu6tJalxPYmNH2IPNCEReNOn4849nc5RC/GmTQVT38570nVJRVMXpdVuNlW7afQEAZiem8SK4iwWjk9nak4S1hNwHosKqIMCHVC36dNQatxu3G6Lh9HTIXsW5J1sLBgYQ5ou2d6q886BKC/vj1Dh1hHAjHQ4swDmZcEgn2ijKMqJJiUObzUJLWXEt2zG5akEQLPE0Zl1Ml2Zs/GlTcObOq1Ptf+klNS6A2yqdrPpgJt9zV4AXDYzcwpTmTc2jXljU/stsEZmQIX90LjNWEK9vszYthv/wZgskDEZsktgdAmkjQdT7BJC0yUVbp3SRo0NDVHeq43SEQIBFKfComxYOBrS1NGSoig9ZA53Ete+i7j2ncS5y7H7Gw7dF4zPw5s2HV/aVLypUwkkTTC6BY9jwkVnIMKuBg/lDR52NXiodQcAI7Cm5yYxNTuJKTmJTM1OYmxGPOY+1gUc3gHlbzcmMbTuNbZte43rrXtBasZj4tIhbYIRRGkTIHMyWJ0x+fVdYUl1p85et8bONp2drRo7WjW6wsb9qQ6YmQ6zMoxp4qoskaIosWCKeHF6qnB69uPw7MfZVYUt0Hzofs3iIpBQSDBxDMGEQgKJYwgmFBB2jSLizECae3ZCZWcgwu7uwPqk1Ud1m5+wpgPgtJqZNCqBSaMTGZPuojAtjjHpceSlunp8/tXQDChdg2CnMU7UVW8sV3Fo2wCeOnBXQcD96XNMZkgYDQnZkDoG0iYapYZcx3fWtqZL2oOSloCk1S9pCejUdUmqPDpVncalPfjpPrOZoDARxiVBcQoUp0GWUy0cqCjKiWEOe3B0VWP31WPzN2ILNHVvmxFS/9xjw440ws4swq5RhF1ZRJyZRO0pROwpRB3dW3sKUXsquuXT7h5Nl9R3BKhq8/FJq4+qNh917gCeYPTQYwSQneygIC2OUUkOshIdjEp0kJVoJ7P7emqcDYfV3OOA6t2omxaG9k+MINGjxkUevN691SIQDULEb5zkGvEbXW8Hr0cCEPEZtwU7jbDxt0PQbfx8GLo1nqgzDc2RSiRzDuH4HMJxOYTicgi7MtGFBV2CJiEUlYTaIdwSIaRBSIOwJruvS0JR8IQl3oikKyzxho2fu8KS9oCkLSjRD5PZGU4YHQdzMiEnzrieEw+5cao2nqIoA0ezJeJLm4Yvbdrnbhd6FGugGVugGUvQjTXUjiXUjjXYjrOzksTmjVgiXUd8Xd1kQbe40CxxaNZ4Zljj0axxaBYXWko8emYcAeHEE7XQETbhDgnaQ4L2NnA3mKgKCXZLC2EshKSVKBY0TJjNPR9a6V1ANe2EB0p69ZTPMduMSgwWu7G1x4M9EdLHG1t7AjgS2etzccdGCw0ylSaZQiDogCPux1D3peesJoizCuKsApcVXFZBmtPE2GRBqkOQ4jSR4ui+7jCR5hQ4LOqQSFGUISYhlTCTCB/pfl3DHOnCHPJgCnswH7p0YYr4MEUDmKLB7m0AS7AVk2ZcN0UCmKJ+xMFhlS86yqlWPf007VUXnxCiC9jT4ycox5IOtA50I4YRtT9jS+3P2FL781MFUspjns/T24n1e3rSb6j0jBCiVO3P2FH7M7bU/owttT97T42eKIqiKIOSCihFURRlUOptQP2pX1oxcqn9GVtqf8aW2p+xpfZnL/VqkoSiKIqinCiqi09RFEUZlFRAKYqiKINSjwJKCHGmEGKPEGKfEOLH/d2o4U4IUSWE2C6E2CKEGATFDYceIcRjQohmIcSOz9yWKoR4XQixt3sb2xL1w9gR9uedQoi67vfpFiHE2QPZxqFCCJEnhHhbCLFLCLFTCPH97tvV+7OXjhlQQggz8EfgLKAYuEwIUdzfDRsBlkkpS9R5EcftceDML9z2Y+BNKeUE4M3un5WeeZwv70+A+7vfpyVSypdPcJuGqihwq5RyMjAPuKH7M1O9P3upJ0dQc4F9Usr9Usow8BRwfv82S1GOTkr5LtD+hZvPB/6n+/r/ABec0EYNYUfYn8pxkFI2SCnLuq93AbuAHNT7s9d6ElA5QM1nfq7tvk05fhJ4TQixSQhx7UA3ZhjJklI2gPEhAWQOcHuGgxuFENu6uwBVl1QvCSEKgZnAx6j3Z6/1JKAOV9dPzU3vm4VSylkY3aY3CCGWDHSDFOUwHgbGASVAA/CbgW3O0CKEiAeeA26WUnoGuj1DUU8CqhbI+8zPuUB9/zRnZJBS1ndvm4E1GN2oSt81CSFGA3Rvm4/xeOUopJRNUkpNSqkDf0a9T3tMCGHFCKcnpJTPd9+s3p+91JOA2ghMEEKMEULYgEuBtf3brOFLCBEnhEg4eB1YAew4+rOUHloLXNl9/UrgxQFsy5B38MO024Wo92mPCCEE8CiwS0r528/cpd6fvdSjShLd00t/B5iBx6SUv+zvhg1XQoixGEdNYFST/1+1P3tPCPEkcArGEgZNwH8DLwBPA/nAAeASKaUa+O+BI+zPUzC69yRQBVx3cAxFOTIhxCLgPWA7cHBJ29sxxqHU+7MXVKkjRVEUZVBSlSQURVGUQUkFlKIoijIoqYBSFEVRBiUVUIqiKMqgpAJKURRFGZRUQCmKoiiDkgooRVEUZVD6P1NsKc7FHrL8AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, (ax1, ax2) = plt.subplots(2, sharex=True)\n",
    "\n",
    "sns.kdeplot(df[(df.categorie == '2') & (df.day_type == 'week')].time.apply(lambda x: x.hour),\n",
    "            shade=True, legend='GCUM', label='GCUM', ax=ax1)\n",
    "sns.kdeplot(df[(df.categorie != '2') & (df.day_type == 'week')].time.apply(lambda x: x.hour),\n",
    "            shade=True, label='Autres', ax=ax1)\n",
    "# sns.kdeplot(df[df.day_type == 'week'].time.apply(lambda x: x.hour),\n",
    "#             shade=True, label='Total', ax=ax1)\n",
    "sns.kdeplot(df[(df.categorie == '2') & (df.day_type == 'weekend')].time.apply(lambda x: x.hour),\n",
    "            shade=True, legend='GCUM', label='GCUM', ax=ax2)\n",
    "sns.kdeplot(df[(df.categorie != '2') & (df.day_type == 'weekend')].time.apply(lambda x: x.hour),\n",
    "            shade=True, label='Autres', ax=ax2)\n",
    "# sns.kdeplot(df[df.day_type == 'weekend'].time.apply(lambda x: x.hour),\n",
    "#             shade=True, label='Total', ax=ax2)\n",
    "\n",
    "ax1.set_title('Semaine')\n",
    "ax2.set_title('Weekend')\n",
    "\n",
    "ax1.set_xbound(0, 24)\n",
    "\n",
    "ax1.get_yaxis().set_visible(False)\n",
    "ax2.get_yaxis().set_visible(False)\n",
    "\n",
    "plt.tight_layout()\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}