test_haze_detection_v1.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%load_ext autoreload\n",
    "%autoreload 2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "solve_haze_detection 2x1, 2 pixels: solve_haze_detection 2x1, 2 pixels: 0.0 seconds, 0.0 seconds CPU\n"
     ]
    }
   ],
   "source": [
    "import common\n",
    "common.regression_test()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Low resolution haze detection v1\n",
    "\n",
    "Results look good, but it requires subsampling and there's a 5000x performance cliff between 600 and 726 pixels, for the sample image\n",
    "![Performance Chart](https://docs.google.com/spreadsheets/d/1VZqHtatv2zCw0UsI8FsD_WXfeI_6b4PR6LFLl-2Dnt8/pubchart?oid=801690291&format=image)\n",
    "\n",
    "Source: https://docs.google.com/spreadsheets/d/1VZqHtatv2zCw0UsI8FsD_WXfeI_6b4PR6LFLl-2Dnt8/edit?gid=0#gid=0\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(22, 33, 4)\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGMCAYAAAD5g/maAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA9/ElEQVR4nO3deXhU9d3//9eZJZMASVAIWQRZlR1UlBg3UFJDaqmg5VZuexGQ4lULbTVupV8LuFyNS916Q6W9eyO11qLcVVCLVFmtZbEsuZX7Vgo0ELaERZOQQJLJzPn94Y+575EEMp+TOCfh+biuua7MzHnPeZ/P+Zwz75w553ws27ZtAQAAuJgn3gkAAACcDQULAABwPQoWAADgehQsAADA9ShYAACA61GwAAAA16NgAQAArkfBAgAAXI+CBQAAuB4FCwB8DaZMmaJOnTrFOw2gzaJgAVxs0aJFsixLmzdvjncqkqQTJ05o7ty5Wrt2bbxTAXCOoWAB0GwnTpzQI488QsEC4GtHwQIAAFyPggVoY06dC3HgwAGNHz9enTp1Ulpamu6//36FQqHIdHv27JFlWfrFL36h5557Tj179lRSUpJGjRql7du3R33m6NGjNXr06Ebn1atXr8jnpaWlSZIeeeQRWZYly7I0d+5cSVIwGNRnn32mQ4cOtdgySFJNTY3uu+8+9ejRQ4FAQP3799cvfvELNWeg+Z07d+rWW29VRkaGEhMT1b17d91+++2qrKyMTGNZlmbOnKklS5Zo0KBBSkpKUk5Ojj755BNJ0q9//Wv169dPiYmJGj16tPbs2XPafJYsWaIRI0YoKSlJXbt21Xe/+10dOHDgrPkVFxcrLS1No0ePVnV1tSTpwIEDuvPOO5Wenq5AIKDBgwdr4cKFZ/0soL3zxTsBALELhULKy8tTdna2fvGLX2jlypV65pln1LdvX919991R07788ss6fvy4ZsyYodraWr3wwgu64YYb9Mknnyg9Pb3Z80xLS9OLL76ou+++WxMmTNAtt9wiSRo2bJikL79oBw4cqIKCAi1atKhFlsG2bX3729/WmjVrNG3aNF1yySX6y1/+ogceeEAHDhzQc8891+Tn19fXKy8vT3V1dfrhD3+ojIwMHThwQO+8844qKiqUmpoamfavf/2r3nrrLc2YMUOSVFRUpG9961t68MEH9atf/Uo/+MEP9MUXX+ipp57SnXfeqdWrV0diFy1apKlTp+qKK65QUVGRysvL9cILL+hvf/ubtm3bps6dOzea39///nfl5eXp8ssv17Jly5SUlKTy8nJdeeWVkSIqLS1N7777rqZNm6aqqirdc889Z21XoN2yAbjWSy+9ZEuy//73v0deKygosCXZjz76aNS0l156qT1ixIjI85KSEluSnZSUZO/fvz/y+qZNm2xJ9r333ht5bdSoUfaoUaNOm39BQYHds2fPyPMjR47Ykuw5c+acNu2p+RUUFJx1uZq7DEuXLrUl2Y8//njUdN/5zndsy7LsXbt2NTmPbdu22ZLsJUuWnDEXSXYgELBLSkoir/3617+2JdkZGRl2VVVV5PVZs2bZkiLT1tfX2926dbOHDBlinzx5MjLdO++8Y0uyZ8+eHbXMHTt2tG3btj/88EM7JSXFvummm+za2trINNOmTbMzMzPto0ePRuV4++2326mpqfaJEyfOuCxAe8ZPQkAb9f3vfz/q+bXXXqt//vOfp003fvx4XXDBBZHnI0eOVHZ2tpYvX96i+fTq1Uu2bTfr6MopZ1uG5cuXy+v16kc/+lHUdPfdd59s29a7777b5GefOoLyl7/8RSdOnDhjHmPGjIn89CVJ2dnZkqRbb71VycnJp71+KsfNmzfr8OHD+sEPfqDExMTIdDfddJMGDBigP//5z6fNa82aNcrLy9OYMWP0xhtvKBAISPryaNKf/vQnjRs3TrZt6+jRo5FHXl6eKisrtXXr1jMuB9CeUbAAbVBiYmLkfJJTzjvvPH3xxRenTXvRRRed9trFF1/c6LkYX6fmLMPevXuVlZUVVTRI0sCBAyPvN6V3794qLCzUb3/7W3Xt2lV5eXmaP39+1Pkrp1x44YVRz08VOz169Gj09VM5npp///79T/vMAQMGnJZfbW2tbrrpJl166aV6/fXXlZCQEHnvyJEjqqio0G9+8xulpaVFPaZOnSpJOnz4cJPLC7R3nMMCtEFer7dFP8+yrEZPYv3qCbAtqaWXoTHPPPOMpkyZomXLlum9997Tj370IxUVFWnjxo3q3r37WXNp6vXG2qo5AoGAvvnNb2rZsmVasWKFvvWtb0XeC4fDkqTvfve7KigoaDT+1PlCwLmIggVo53bu3Hnaa//4xz+ifgI577zzGv056atHCCzLavH8zqRnz55auXKljh8/HnWU5bPPPou8fzZDhw7V0KFD9fDDD2v9+vW6+uqrtWDBAj3++OMtkp8k7dixQzfccEPUezt27DgtP8uy9Ic//EE333yzJk6cqHfffTdydVZaWpqSk5MVCoWUm5vrODegveEnIaCdW7p0adQlth999JE2bdqk/Pz8yGt9+/bVZ599piNHjkRe+6//+i/97W9/i/qsDh06SJIqKipOm08slzU31ze/+U2FQiHNmzcv6vXnnntOlmVFLcNXVVVVqaGhIeq1oUOHyuPxqK6urkXyu/zyy9WtWzctWLAg6jPfffddffrpp7rppptOi0lISNAbb7yhK664QuPGjdNHH30k6cujObfeeqv+9Kc/nXbZuaSodQOcizjCArRz/fr10zXXXKO7775bdXV1ev7559WlSxc9+OCDkWnuvPNOPfvss8rLy9O0adN0+PBhLViwQIMHD1ZVVVVkuqSkJA0aNEivvfaaLr74Yp1//vkaMmSIhgwZEvNlzc0xbtw4XX/99fp//+//ac+ePRo+fLjee+89LVu2TPfcc4/69u3bZOzq1as1c+ZMTZw4URdffLEaGhr0+9//PlIYtAS/368nn3xSU6dO1ahRozRp0qTIZc29evXSvffe22hcUlKS3nnnHd1www3Kz8/XunXrNGTIED3xxBNas2aNsrOzNX36dA0aNEiff/65tm7dqpUrV+rzzz9vkbyBtoiCBWjnJk+eLI/Ho+eff16HDx/WyJEjNW/ePGVmZkamGThwoF5++WXNnj1bhYWFGjRokH7/+9/r1VdfPe02/L/97W/1wx/+UPfee6/q6+s1Z84cDRkypFVy93g8euuttzR79my99tpreumll9SrVy89/fTTuu+++84YO3z4cOXl5entt9/WgQMH1KFDBw0fPlzvvvuurrzyyhbLccqUKerQoYOeeOIJPfTQQ+rYsaMmTJigJ598ssl7sEhSSkqK/vKXv+i6667TN77xDf31r39Vv3799NFHH+nRRx/VG2+8oV/96lfq0qWLBg8erCeffLLFcgbaIss2PXsMgKvt2bNHvXv31tNPP637778/3ukAgCOcwwIAAFyPggUAALgeBQsAAHA9zmEBAACuxxEWAADgehQsAADA9drFfVjC4bAOHjyo5OTkr/3W4QAAwIxt2zp+/LiysrLk8Zz5GEq7KFgOHjx42qiqAACgbdi3b1/UgKSNaRcFy6lB0W7KvVZ+X+yLZDs4KFNfbz4miccy/0Wu5kStcWyHQMLZJ2qCk/F1fQnm860LmZ8bXlNTbxxbWW0e63XQrzp26mAcm+A171fBr4y9EwuPx3x34vP7jWNlObtuwPKYr6iwg/Zy8ot8uCFoHOvzm8/XtpzsARycgeBg1HDbDhvHJndKMo5NSjTv019UnTSODTvY74QdbEqm1+80NDRo7QfrogY3bUq7KFhO/Qzk9/nk93+9BYsdNt+QnBQsPoPC7BSTou4URwWLwbo5JezgS8nnM99heb3m69fr4OdJJ+vX56BgcfLV76RgcdIn41mwmPcOyXLwBR5ysKb8vngVLA5iHWxLTgoWv4NCOsFBrN9nXgg7K1jM+5XTC46bczpHq510O3/+fPXq1UuJiYnKzs6OjEjalCVLlmjAgAFKTEzU0KFDtXz58tZKDQAAtDGtUrC89tprKiws1Jw5c7R169bIIGSHDx9udPr169dr0qRJmjZtmrZt26bx48dr/PjxjQ6xDgAAzj2tUrA8++yzmj59uqZOnapBgwZpwYIF6tChgxYuXNjo9C+88ILGjh2rBx54QAMHDtRjjz2myy67TPPmzWuN9AAAQBvT4gVLfX29tmzZotzc3P+dicej3NxcbdiwodGYDRs2RE0vSXl5eU1OX1dXp6qqqqgHAABov1q8YDl69KhCoZDS09OjXk9PT1dZWVmjMWVlZTFNX1RUpNTU1MiDS5oBAGjf2uSdbmfNmqXKysrIY9++ffFOCQAAtKIWv6y5a9eu8nq9Ki8vj3q9vLxcGRkZjcZkZGTENH0gEFAgEGiZhAEAgOu1+BGWhIQEjRgxQqtWrYq8Fg6HtWrVKuXk5DQak5OTEzW9JL3//vtNTg8AAM4trXLjuMLCQhUUFOjyyy/XyJEj9fzzz6umpkZTp06VJE2ePFkXXHCBioqKJEk//vGPNWrUKD3zzDO66aabtHjxYm3evFm/+c1vWiM9AADQxrRKwXLbbbfpyJEjmj17tsrKynTJJZdoxYoVkRNrS0tLowY5uuqqq/Tqq6/q4Ycf1k9/+lNddNFFWrp0qYYMGdIa6QEAgDbGsp3eT9cFqqqqlJqaqvFjr//ab81fXxefsYSqHYwl1DFeYwk5mG9dg3k3rXYwllDFcfP16+TW/J2Sz62xhJzcAj2ut+Z30F6Obs3vYCwhf9zGEnIQGzJvZ0djCTnYDjs4GEvo88p4jSX09d+aP9jQoJWrV6myslIpKSlnnLZdjCV0yqHy4/J5Y18kSw4GP7PMNwaFHezsHHRKr4Me7Xewcy87Wm0cWxdyMt6Lgy80B4Mudkk73zi2U4dE41gn3991QfMCTQ76ht9vHutgU/iSg/ayfOZFeDDoZBwy8y//pCTznJ20VW2teZF1os78H7QkB/8o2Q720baDnpmSbF7shELm30m27aCINpxtMNj8ebbJy5oBAMC5hYIFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA61GwAAAA16NgAQAArkfBAgAAXI+CBQAAuB4FCwAAcD0KFgAA4HoULAAAwPUoWAAAgOv54p1AS+qUGJDPF/siebzmzWDb5kN5V9fUGMfWmY96rhMN5kPT26F649iUlI7GsV6D9XqK32u+vH6veU1vOcjZydD0TvhkPqy9JfNtIcFBO0u2g1hJloN17GC2dtg8b9tBW3s98VnesIN9ZWWN+X6nIWQcqg4dA8axtbXmO2mPx7ylLctBrJNtyTDncAzdkSMsAADA9ShYAACA61GwAAAA16NgAQAArkfBAgAAXI+CBQAAuB4FCwAAcD0KFgAA4HoULAAAwPUoWAAAgOtRsAAAANdr8YKlqKhIV1xxhZKTk9WtWzeNHz9eO3bsOGPMokWLZFlW1CMxMbGlUwMAAG1Uixcs69at04wZM7Rx40a9//77CgaDuvHGG1VzloH+UlJSdOjQochj7969LZ0aAABoo1p8tOYVK1ZEPV+0aJG6deumLVu26LrrrmsyzrIsZWRktHQ6AACgHWjxguWrKisrJUnnn3/+Gaerrq5Wz549FQ6Hddlll+nnP/+5Bg8e3Oi0dXV1qqurizyvqqqSJI0c3kuJgYSYc6w6fuajP2fSEDQfQvx4jfnPXjv2fe5gvieMY7skm3eZfj1SjGNTO3Y0jpX5qPayfF7j2P1Hzdu5ouakcazHMs/Zwaj2sjzmB2zDYfOVZBtHfsly9AnmscGQ+b6jPhgyjk1IMI+Vbb68DQ7Wcaek2PfrpyT4zbeHhpD5BmE5aCuPx0Gs+eLKdtCfLcMfbMLh5s+zVU+6DYfDuueee3T11VdryJAhTU7Xv39/LVy4UMuWLdMrr7yicDisq666Svv37290+qKiIqWmpkYePXr0aK1FAAAALtCqBcuMGTO0fft2LV68+IzT5eTkaPLkybrkkks0atQovfHGG0pLS9Ovf/3rRqefNWuWKisrI499+/a1RvoAAMAlWu0noZkzZ+qdd97RBx98oO7du8cU6/f7demll2rXrl2Nvh8IBBQIBFoiTQAA0Aa0+BEW27Y1c+ZMvfnmm1q9erV69+4d82eEQiF98sknyszMbOn0AABAG9TiR1hmzJihV199VcuWLVNycrLKysokSampqUpKSpIkTZ48WRdccIGKiookSY8++qiuvPJK9evXTxUVFXr66ae1d+9efe9732vp9AAAQBvU4gXLiy++KEkaPXp01OsvvfSSpkyZIkkqLS2V5/9cUfDFF19o+vTpKisr03nnnacRI0Zo/fr1GjRoUEunBwAA2qAWL1jsZlzKtXbt2qjnzz33nJ577rmWTgUAALQTjCUEAABcj4IFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA61GwAAAA12u1sYTiIaVTkhIDsQ9D/umOxscsao6uXVKMY0MyH+Y9GDQfmt7nYLh1v8d8uPVwvXnOieeZ52xb5jl7HNT0loP1a8nBsPYOhoh3MFvJQTs7yFi2o2jJasa9o5riseLzP5/lbJGN2Q76h4Ndh87r3Mk4NsFvvo7CtnlsqCFsHOtx0K0sB9uhk83f9PiHFcM2xBEWAADgehQsAADA9ShYAACA61GwAAAA16NgAQAArkfBAgAAXI+CBQAAuB4FCwAAcD0KFgAA4HoULAAAwPUoWAAAgOtRsAAAANejYAEAAK5HwQIAAFzPF+8EWtKePfuV4PfHHJfgTTCep8dB7OHPK4xjg0HjUHl85mPTB0Mh89gG8/mGHYxr74lh+PKvCoYctJWDlRRqaDCOtfzmm7UdNp9vqCFsHOtoYHsH61eK339twQYH21LYvK2DDvpW2HawDTvYd8jB8ipsvoZt27ytGkLmOTsIlcdjvi1ZHvO2Mp1rLNsBR1gAAIDrUbAAAADXo2ABAACu1+IFy9y5c2VZVtRjwIABZ4xZsmSJBgwYoMTERA0dOlTLly9v6bQAAEAb1ipHWAYPHqxDhw5FHh9++GGT065fv16TJk3StGnTtG3bNo0fP17jx4/X9u3bWyM1AADQBrVKweLz+ZSRkRF5dO3atclpX3jhBY0dO1YPPPCABg4cqMcee0yXXXaZ5s2b1xqpAQCANqhVCpadO3cqKytLffr00R133KHS0tImp92wYYNyc3OjXsvLy9OGDRuajKmrq1NVVVXUAwAAtF8tXrBkZ2dr0aJFWrFihV588UWVlJTo2muv1fHjxxudvqysTOnp6VGvpaenq6ysrMl5FBUVKTU1NfLo0aNHiy4DAABwlxYvWPLz8zVx4kQNGzZMeXl5Wr58uSoqKvT666+32DxmzZqlysrKyGPfvn0t9tkAAMB9Wv1Ot507d9bFF1+sXbt2Nfp+RkaGysvLo14rLy9XRkZGk58ZCAQUCARaNE8AAOBerX4flurqau3evVuZmZmNvp+Tk6NVq1ZFvfb+++8rJyentVMDAABtRIsXLPfff7/WrVunPXv2aP369ZowYYK8Xq8mTZokSZo8ebJmzZoVmf7HP/6xVqxYoWeeeUafffaZ5s6dq82bN2vmzJktnRoAAGijWvwnof3792vSpEk6duyY0tLSdM0112jjxo1KS0uTJJWWlsrzfwZYuuqqq/Tqq6/q4Ycf1k9/+lNddNFFWrp0qYYMGdLSqQEAgDaqxQuWxYsXn/H9tWvXnvbaxIkTNXHixJZOBQAAtBOtftLt1+mT3Ufk88a+SCfrzYcQD5UeNY4NOhr23HzVWQ6GtT9abRyqo58dMY7d+Gn52Sdqivlo6wrbToa1N5+x1+s1jrVlnrNlm/eNpCS/cazPbx5rG0d+KcFjvp4sB32rNmieuW2bz7ghWG8cW1sfNI51sj0kBpKMYxtC5n067KB3nQyat1V90Pw7ye8x/15xEKpQ2KytGhqav6wMfggAAFyPggUAALgeBQsAAHA9ChYAAOB6FCwAAMD1KFgAAIDrUbAAAADXo2ABAACuR8ECAABcj4IFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA6/ninUBL8vp98vpiX6Sgg2HAbXmNY0MOhj33esyHiO+UFDCOrQ2aD7de42Bo+lCD+XwdjGov2zafr8fBfM17lVTfEDaOdTI0fSd/kvl8/eZLHA6bL68kWZb5inIQKo/HPO9Q2LxfyjJva49lvg03ONgQT9aZ7+98PvP5Jhh8n5xiOeiXHst8O7Q95jnbtnnOSYlm8w02NL8/coQFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA61GwAAAA16NgAQAArkfBAgAAXI+CBQAAuB4FCwAAcD0KFgAA4HotXrD06tVLlmWd9pgxY0aj0y9atOi0aRMTE1s6LQAA0Ia1+OCHf//736MG9du+fbu+8Y1vaOLEiU3GpKSkaMeOHZHnTgYjAwAA7U+LFyxpaWlRz5944gn17dtXo0aNajLGsixlZGS0dCoAAKCdaPGC5f+qr6/XK6+8osLCwjMeNamurlbPnj0VDod12WWX6ec//7kGDx7c5PR1dXWqq6uLPK+qqpIkhUNhha3Yh8cO+M1/GautD519oqZ4HAx7nmCec2JignFsbbDWONYf8BvHBgLmbRV2MGR60MEQ8T6vg83Lto1DPd7mD9f+VX6veTsnd+zgYL7mOYcdrCNJOlFr3qfrQ+bryXZwJNmyHMw3bB6b6DffdwQdHDivdbCObAdt1VBfbxxbHzTvlx6feTsnp6Yax4Yt8++VkOG2EI4hrlVPul26dKkqKio0ZcqUJqfp37+/Fi5cqGXLlumVV15ROBzWVVddpf379zcZU1RUpNTU1MijR48erZA9AABwC8u2HfwrdxZ5eXlKSEjQ22+/3eyYYDCogQMHatKkSXrssccanaaxIyw9evTQ1SOukM8X+3+1taFgzDGRWAdHWIIN5k2f5OAIS0on8/+EK447+G/UNv8XyyOOsDSX6X86krMjLBekmf9n11aPsAQdHK1wsuu1HSyz32O+7/A6OFoRbDAOleVxciSLIyzNZcm8rRIM9x0NDQ1auW6NKisrlZKScsZpW+0nob1792rlypV64403Yorz+/269NJLtWvXrianCQQCCgQCTlMEAABtRKv9JPTSSy+pW7duuummm2KKC4VC+uSTT5SZmdlKmQEAgLamVQqWcDisl156SQUFBaf9RDN58mTNmjUr8vzRRx/Ve++9p3/+85/aunWrvvvd72rv3r363ve+1xqpAQCANqhVfhJauXKlSktLdeedd572XmlpqTz/5zfUL774QtOnT1dZWZnOO+88jRgxQuvXr9egQYNaIzUAANAGtepJt1+XqqoqpaamctJtM3HSbfNx0m3zcdJtbDjptvk46bb52vNJt4wlBAAAXI+CBQAAuB4FCwAAcD0KFgAA4HoULAAAwPUoWAAAgOtRsAAAANdrtbGE4qG2vkE+g3tQ+BP8xvMM+c2b0OCWMRHJ5imrwcENEfwOljfgM0/a0e2CLPN7iwQbzO+z47HM7y0iB/dDkIN71sg2X94GB/d/kWU+X1vO7sNieR3cl8TBevI4WE1h27xv2Q3m239iUqJ5rOXg/2MH27CT9XvgULlx7BdVJ4xju3Y1v1eW18HyOrmnUX3Q7H5mDTH0R46wAAAA16NgAQAArkfBAgAAXI+CBQAAuB4FCwAAcD0KFgAA4HoULAAAwPUoWAAAgOtRsAAAANejYAEAAK5HwQIAAFyPggUAALgeBQsAAHA9ChYAAOB6vngn0JLCshRS7MOQOxki3g6ZD9UecjCUd3XIPOcTtWbDgEuSz2de4/r8CcaxobD58PLmrSzJNm9nr22es2UeqrCDnIPBeuPYUG2tcazPb77AIY/58kqSZZv3aTsYMo5NSko0jj3hYD1ZBvvIUzwO+qWcbA8eB+vIQeyJWvP1Gwo5aKwGB9twnfl3UoNtvrdMSfQaxXms5i8rR1gAAIDrUbAAAADXo2ABAACuR8ECAABcL+aC5YMPPtC4ceOUlZUly7K0dOnSqPdt29bs2bOVmZmppKQk5ebmaufOnWf93Pnz56tXr15KTExUdna2Pvroo1hTAwAA7VTMBUtNTY2GDx+u+fPnN/r+U089pV/+8pdasGCBNm3apI4dOyovL0+1Z7iC4LXXXlNhYaHmzJmjrVu3avjw4crLy9Phw4djTQ8AALRDMRcs+fn5evzxxzVhwoTT3rNtW88//7wefvhh3XzzzRo2bJhefvllHTx48LQjMf/Xs88+q+nTp2vq1KkaNGiQFixYoA4dOmjhwoWNTl9XV6eqqqqoBwAAaL9a9ByWkpISlZWVKTc3N/JaamqqsrOztWHDhkZj6uvrtWXLlqgYj8ej3NzcJmOKioqUmpoaefTo0aMlFwMAALhMixYsZWVlkqT09PSo19PT0yPvfdXRo0cVCoViipk1a5YqKysjj3379rVA9gAAwK3a5J1uA4GAAoFAvNMAAABfkxY9wpKRkSFJKi8vj3q9vLw88t5Xde3aVV6vN6YYAABwbmnRgqV3797KyMjQqlWrIq9VVVVp06ZNysnJaTQmISFBI0aMiIoJh8NatWpVkzEAAODcEvNPQtXV1dq1a1fkeUlJiYqLi3X++efrwgsv1D333KPHH39cF110kXr37q2f/exnysrK0vjx4yMxY8aM0YQJEzRz5kxJUmFhoQoKCnT55Zdr5MiRev7551VTU6OpU6c6X0IAANDmxVywbN68Wddff33keWFhoSSpoKBAixYt0oMPPqiamhrdddddqqio0DXXXKMVK1YoMfF/RybdvXu3jh49Gnl+22236ciRI5o9e7bKysp0ySWXaMWKFaediAsAAM5Nlm07GIveJaqqqpSamqrJ3x6lBH/s5xHbDobUDpmP5K2Dnzd9M72zqTphPrx85fGTxrGDe3Y2ju2SYn6Ot9fBEPGW5eSXT/Ph5e2weaxs85wbQubzDYbNtwWfxzKO9ThYR7b5bCVJwZD5B3x6oM44NmSZz/fC8xOMY/2Wef/onJJkHGs5+KbxOOhbIQfz3faPI8axX9SY76MH9+xiHNvtPPN1FDTf/NUxyaxP1geD+vXr76myslIpKSlnnJaxhAAAgOtRsAAAANejYAEAAK5HwQIAAFyPggUAALgeBQsAAHA9ChYAAOB6FCwAAMD1KFgAAIDrUbAAAADXo2ABAACuR8ECAABcj4IFAAC4HgULAABwPQoWAADger54J9CS9pdXyefzxhwXCASM53nk8+PGsbLNQz0Gy3lK2DKf7+Evao1jK6vDxrE+B6W112PeVj5fyDg2kGA+304J5ptmgs9JrHnOtoN+5Zd5sMdyMGNJdTLvl14H8w42GIcqHDYPTjTf3Skl0Xx7cNJBnPxn7XGw8+iSYh7bYJtvSx0crKO0ZPNYJ9twKFxvFOezm9+XOcICAABcj4IFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA61GwAAAA16NgAQAArkfBAgAAXI+CBQAAuF7MBcsHH3ygcePGKSsrS5ZlaenSpZH3gsGgHnroIQ0dOlQdO3ZUVlaWJk+erIMHD57xM+fOnSvLsqIeAwYMiHlhAABA+xRzwVJTU6Phw4dr/vz5p7134sQJbd26VT/72c+0detWvfHGG9qxY4e+/e1vn/VzBw8erEOHDkUeH374YaypAQCAdirmUdLy8/OVn5/f6Hupqal6//33o16bN2+eRo4cqdLSUl144YVNJ+LzKSMjI9Z0AADAOaDVz2GprKyUZVnq3LnzGafbuXOnsrKy1KdPH91xxx0qLS1tctq6ujpVVVVFPQAAQPtlPg59M9TW1uqhhx7SpEmTlJKS0uR02dnZWrRokfr3769Dhw7pkUce0bXXXqvt27crOfn0sbKLior0yCOPnPb6SfnlM1ikE7V2zDGnlFcEjWNTO/qNY/122Di2+kStcazdYD68vM+XYB7rNR+qPRAwr8vNl1by+8zHavd7zOccDjV/uPavsi3zbcHrNW9nj2XeVh6Zx0rSyXrzbami2nw9WQ624ZIG82XukppoHBvymn9leB1sw3LQP8IO/i//4qT5fE+GzOe7o8x8H308aP69Uh8033cEw2bbQkND8+fZakdYgsGg/uVf/kW2bevFF18847T5+fmaOHGihg0bpry8PC1fvlwVFRV6/fXXG51+1qxZqqysjDz27dvXGosAAABcolWOsJwqVvbu3avVq1ef8ehKYzp37qyLL75Yu3btavT9QCCgQCDQEqkCAIA2oMWPsJwqVnbu3KmVK1eqS5cuMX9GdXW1du/erczMzJZODwAAtEExFyzV1dUqLi5WcXGxJKmkpETFxcUqLS1VMBjUd77zHW3evFl/+MMfFAqFVFZWprKyMtXX10c+Y8yYMZo3b17k+f33369169Zpz549Wr9+vSZMmCCv16tJkyY5X0IAANDmxfyT0ObNm3X99ddHnhcWFkqSCgoKNHfuXL311luSpEsuuSQqbs2aNRo9erQkaffu3Tp69Gjkvf3792vSpEk6duyY0tLSdM0112jjxo1KS0uLNT0AANAOxVywjB49Wrbd9JUEZ3rvlD179kQ9X7x4caxpAACAcwhjCQEAANejYAEAAK5HwQIAAFyPggUAALgeBQsAAHA9ChYAAOB6FCwAAMD1KFgAAIDrtcrgh/FSe6JBXu/Zb1z3VbEMb316bNA41qME41ifx3zVdUw0HzgyNdF86PLOyR2MYy3bbOhySUpJMl/ek7Xm891fddI4ttYyH9ZesW8CEZaDf2H8/rB5rNe8P3vkoK0khczTlm2b7ztCDtZTSOZ9uvyLWuPYk/XmbZ0UMN93eB10TMtv3rfqG7zGscGg+XdDTdi8Ux44dsI4NtSMG782zSznUKj5+1iOsAAAANejYAEAAK5HwQIAAFyPggUAALgeBQsAAHA9ChYAAOB6FCwAAMD1KFgAAIDrUbAAAADXo2ABAACuR8ECAABcj4IFAAC4HgULAABwPQoWAADgeubjbrvQ8ZMn5fXGPhy4g5G85fGaN2HNSfNh3mvrjEPl8yYYx9bVNRjH1lrmSXtiGIL8q5IdlOU9UjoYx+4qrzCObZD5sPZyMkS8g7ayHMzW7zeP9Xot82BJwbB5nw462HlYHvPGrg+Z5xx2kHPVyXrjWFnmfdrvM28rj8e8Y9bUmS9vfYP5Pitsm/dp28F3gyXz+VqGoaFQ89cPR1gAAIDrUbAAAADXo2ABAACuF3PB8sEHH2jcuHHKysqSZVlaunRp1PtTpkyRZVlRj7Fjx571c+fPn69evXopMTFR2dnZ+uijj2JNDQAAtFMxFyw1NTUaPny45s+f3+Q0Y8eO1aFDhyKPP/7xj2f8zNdee02FhYWaM2eOtm7dquHDhysvL0+HDx+ONT0AANAOxXyJS35+vvLz8884TSAQUEZGRrM/89lnn9X06dM1depUSdKCBQv05z//WQsXLtRPfvKTWFMEAADtTKucw7J27Vp169ZN/fv31913361jx441OW19fb22bNmi3Nzc/03K41Fubq42bNjQaExdXZ2qqqqiHgAAoP1q8YJl7Nixevnll7Vq1So9+eSTWrdunfLz8xVq4j4aR48eVSgUUnp6etTr6enpKisrazSmqKhIqampkUePHj1aejEAAICLtPiN426//fbI30OHDtWwYcPUt29frV27VmPGjGmRecyaNUuFhYWR51VVVRQtAAC0Y61+WXOfPn3UtWtX7dq1q9H3u3btKq/Xq/Ly8qjXy8vLmzwPJhAIKCUlJeoBAADar1YvWPbv369jx44pMzOz0fcTEhI0YsQIrVq1KvJaOBzWqlWrlJOT09rpAQCANiDmgqW6ulrFxcUqLi6WJJWUlKi4uFilpaWqrq7WAw88oI0bN2rPnj1atWqVbr75ZvXr1095eXmRzxgzZozmzZsXeV5YWKh///d/1+9+9zt9+umnuvvuu1VTUxO5aggAAJzbYj6HZfPmzbr++usjz0+dS1JQUKAXX3xRH3/8sX73u9+poqJCWVlZuvHGG/XYY48pEAhEYnbv3q2jR49Gnt922206cuSIZs+erbKyMl1yySVasWLFaSfiAgCAc5Nl206GdnWHqqoqpaamql+/gV/7aM319eajpiZ6zZve6+THPAejNfvC5iOQpiYGzj5RE5yM1tyns/l8Oyebj9a8fKf5jQ/PvdGazddRPEdrrg+axzoZrdnvM+8fTkZrDjgYVjulQ5JxrN9nfn2Ix2/eVocOHz37RE04WRc0jvV5zZfXn2AeG5/RmkPa8WmxKisrz3o+aotfJRRPiR6zL/Kwg+8Gn2XehPX15uOA1wYdFDveWuPYLAcnOJ+sNd+APQ42pHrbwdD0pluhpICDHYffQcHS1C0EmiO5k/mXykV9za/U++fexm9h0BzHqqqNYyU5O5PPMg/2Ovivw8n/mYGAedER8JnHNgTN+6VlmW8PPr+TL2HzdWSZL65CclD9h8xjPQ7maxoby/6KwQ8BAIDrUbAAAADXo2ABAACuR8ECAABcj4IFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA61GwAAAA16NgAQAArkfBAgAAXI+CBQAAuB4FCwAAcD1fvBNoSZ0TffJ5Yx+GvCbYYD7TgPnQ5VVB89naDoamDzgYbj3gdzDsuYOh2j0yzznoNc+5piFsHNvR7+D/gdi7cYTl4P+Qjknmu4SOCc0fJv6rMtM6GMdWHq8wjpWkxETzeZ8M1xrHBnzmfdrvc9BBZL49+Dzm69ifkGA+X8t8O/R5zJfXwW5WHgeryOM17xuJDvbvYdu8nb2G++iQ1fw4jrAAAADXo2ABAACuR8ECAABcj4IFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA61GwAAAA16NgAQAArhdzwfLBBx9o3LhxysrKkmVZWrp0adT7lmU1+nj66aeb/My5c+eeNv2AAQNiXhgAANA+xVyw1NTUaPjw4Zo/f36j7x86dCjqsXDhQlmWpVtvvfWMnzt48OCouA8//DDW1AAAQDsV80hn+fn5ys/Pb/L9jIyMqOfLli3T9ddfrz59+pw5EZ/vtFgAAACplc9hKS8v15///GdNmzbtrNPu3LlTWVlZ6tOnj+644w6VlpY2OW1dXZ2qqqqiHgAAoP0yH0u+GX73u98pOTlZt9xyyxmny87O1qJFi9S/f38dOnRIjzzyiK699lpt375dycnJp01fVFSkRx555LTXPX7LbFhu8xG1ZVkOaj7bPDbgNx+7PJBgPt+aUNA41usgZ9s2HyL+WL15zl/UNxjHJnVIMI5tCIeMY0Nh8w5dUV1rHLvl4z3GsQ0OcrYNh7U/JdRg3tYy75byec23Q7/HfMYej/l8PeabsBIcfNt4Zb6OkpPMl/egg/mGQuaxfst8/XZ1sN+prXOQs2F/bgg1f/tt1SMsCxcu1B133KHExMQzTpefn6+JEydq2LBhysvL0/Lly1VRUaHXX3+90elnzZqlysrKyGPfvn2tkT4AAHCJVjvC8te//lU7duzQa6+9FnNs586ddfHFF2vXrl2Nvh8IBBQIBJymCAAA2ohWO8LyH//xHxoxYoSGDx8ec2x1dbV2796tzMzMVsgMAAC0NTEXLNXV1SouLlZxcbEkqaSkRMXFxVEnyVZVVWnJkiX63ve+1+hnjBkzRvPmzYs8v//++7Vu3Trt2bNH69ev14QJE+T1ejVp0qRY0wMAAO1QzD8Jbd68Wddff33keWFhoSSpoKBAixYtkiQtXrxYtm03WXDs3r1bR48ejTzfv3+/Jk2apGPHjiktLU3XXHONNm7cqLS0tFjTAwAA7ZBlO7n8wiWqqqqUmpqq60YMk88b+2nsx+vMryKxLPPT5isrzK/I8MTpKqFEv/kVGV6fec4OLopQgsmVY6fm6+AKlJCDvhGvq4Tqg+Z9I+xgvk6uEqp3ECtJCX7zqyrqHew7Oib5jWMTHFyt4+wqIfPtISFg3s5eB5djJXfqaBz7WWm5cezx43XGsYk+83WU2dV8eWsd9Ge/wXevJDWEQtpQ/F+qrKxUSkrKGadlLCEAAOB6FCwAAMD1KFgAAIDrUbAAAADXo2ABAACuR8ECAABcj4IFAAC4HgULAABwvVYb/DAeGsJhyYr9xkbBkIMbXjmItQ1vtCNJHsu81vTKwfIGndyky/zmT+a3UZPCYQc3Q3OQs+wG81AHN50LObjpnJPbSJrfUkwK+MyjAx4nc5Zkm98sKyHJvG/5vObbUjhsvqJq681vaObzm39lnKg7aRzbMdH8pnPBSuNQnagxb6twg/l2WOfg5n7VdebrqC5o3q+qKk4YxYVj2F9xhAUAALgeBQsAAHA9ChYAAOB6FCwAAMD1KFgAAIDrUbAAAADXo2ABAACuR8ECAABcj4IFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA61GwAAAA1zMfh9qFgsEG2QbDrnts8+HpO3jMh+O2fOaxdQ0NxrGKYTjvr0pOMu8yPst8vpLXONK2zdvZYzlYvw76hkdO+pVxqGS+KSgUMs85bDvoz3bYPFZy0NKSPOYN5rHN+7THNv9fMxw2by+fwf71FI9lvrxW0Dzn49U1xrE+B+s3qVPAOLauwXxfmZqaYBzbrUuacWzxJzuN4mLZb3CEBQAAuB4FCwAAcD0KFgAA4HoxFSxFRUW64oorlJycrG7dumn8+PHasWNH1DS1tbWaMWOGunTpok6dOunWW29VeXn5GT/Xtm3Nnj1bmZmZSkpKUm5urnbuNPs9DAAAtD8xFSzr1q3TjBkztHHjRr3//vsKBoO68cYbVVPzvyc13XvvvXr77be1ZMkSrVu3TgcPHtQtt9xyxs996qmn9Mtf/lILFizQpk2b1LFjR+Xl5am2ttZsqQAAQLti2Q4uoThy5Ii6deumdevW6brrrlNlZaXS0tL06quv6jvf+Y4k6bPPPtPAgQO1YcMGXXnllad9hm3bysrK0n333af7779fklRZWan09HQtWrRIt99++1nzqKqqUmpqqq4YMkA+b+xnotc5uEDB5+Aag5MNTq4SMo/1O7hax9FVQo6uSXNylZCDS1+cXCXkYLYeB5frOJmvHFzJwVVCMYY66NOesPmv+XXBeuNYn99vHOtx0DF9XvPYmjrzNVx5ss441u8zX79OrhLq2b2bcWy3LqnGseZXCYX0j3/uVGVlpVJSUs44raNzWCorKyVJ559/viRpy5YtCgaDys3NjUwzYMAAXXjhhdqwYUOjn1FSUqKysrKomNTUVGVnZzcZU1dXp6qqqqgHAABov4wLlnA4rHvuuUdXX321hgwZIkkqKytTQkKCOnfuHDVtenq6ysrKGv2cU6+np6c3O6aoqEipqamRR48ePUwXAwAAtAHGBcuMGTO0fft2LV68uCXzaZZZs2apsrIy8ti3b9/XngMAAPj6GBUsM2fO1DvvvKM1a9aoe/fukdczMjJUX1+vioqKqOnLy8uVkZHR6Gedev2rVxKdKSYQCCglJSXqAQAA2q+YChbbtjVz5ky9+eabWr16tXr37h31/ogRI+T3+7Vq1arIazt27FBpaalycnIa/czevXsrIyMjKqaqqkqbNm1qMgYAAJxbYipYZsyYoVdeeUWvvvqqkpOTVVZWprKyMp08eVLSlyfLTps2TYWFhVqzZo22bNmiqVOnKicnJ+oKoQEDBujNN9+UJFmWpXvuuUePP/643nrrLX3yySeaPHmysrKyNH78+JZbUgAA0GbFdKHpiy++KEkaPXp01OsvvfSSpkyZIkl67rnn5PF4dOutt6qurk55eXn61a9+FTX9jh07IlcYSdKDDz6ompoa3XXXXaqoqNA111yjFStWKDEx0WCRAABAe+PoPixuwX1YYsN9WGLAfViajfuwxBjKfViaP1/uw9Js7fk+LI6+QtymPuRRyOA84gS/+cYQcLAhJXrMm//g0ePGsRW15hvDsZqgcWzYwReL5eCLVA6+/B0VDs6+Co05qc/CTmIdLK5XzooOJ5yspbCDaI+DFeUk1kmRZXvMC0snBYvloHM52e/INt9XJnjNi8rEQIJxrNc230eXHzC/4rZLwKyYbQg1v50Y/BAAALgeBQsAAHA9ChYAAOB6FCwAAMD1KFgAAIDrUbAAAADXo2ABAACuR8ECAABcj4IFAAC4HgULAABwPQoWAADgehQsAADA9ShYAACA67WL0Zpt+8uRPEMhs5E1GxyMIup1MJqv7WDE1XDYfBTRcNh89NKwg1F1nY3W7GB0WgfryEksozU3n3UOjtbsaEU5iXWyz3IwV8VttGYn68jB/s7BKgo52L83NJiPqN3Q4GC+pt+//3+c3Yz11C4KluPHj0uStv9jR5wzAQAgfnaVxDsDM8ePH1dqauoZp7Hs5pQ1LhcOh3Xw4EElJyc3+p94VVWVevTooX379iklJSUOGbYdtFXz0VbNR1vFhvZqPtqq+dzYVrZt6/jx48rKypLHc+azVNrFERaPx6Pu3bufdbqUlBTXrCS3o62aj7ZqPtoqNrRX89FWzee2tjrbkZVTOOkWAAC4HgULAABwvXOiYAkEApozZ44CgUC8U3E92qr5aKvmo61iQ3s1H23VfG29rdrFSbcAAKB9OyeOsAAAgLaNggUAALgeBQsAAHA9ChYAAOB6FCwAAMD12n3BMn/+fPXq1UuJiYnKzs7WRx99FO+UXGnu3LmyLCvqMWDAgHin5QoffPCBxo0bp6ysLFmWpaVLl0a9b9u2Zs+erczMTCUlJSk3N1c7d+6MT7Jxdra2mjJlymn9bOzYsfFJNs6Kiop0xRVXKDk5Wd26ddP48eO1Y0f0eGi1tbWaMWOGunTpok6dOunWW29VeXl5nDKOn+a01ejRo0/rW9///vfjlHH8vPjiixo2bFjkbrY5OTl69913I++35T7VrguW1157TYWFhZozZ462bt2q4cOHKy8vT4cPH453aq40ePBgHTp0KPL48MMP452SK9TU1Gj48OGaP39+o+8/9dRT+uUvf6kFCxZo06ZN6tixo/Ly8lRbW/s1Zxp/Z2srSRo7dmxUP/vjH//4NWboHuvWrdOMGTO0ceNGvf/++woGg7rxxhtVU1MTmebee+/V22+/rSVLlmjdunU6ePCgbrnlljhmHR/NaStJmj59elTfeuqpp+KUcfx0795dTzzxhLZs2aLNmzfrhhtu0M0336z//u//ltTG+5Tdjo0cOdKeMWNG5HkoFLKzsrLsoqKiOGblTnPmzLGHDx8e7zRcT5L95ptvRp6Hw2E7IyPDfvrppyOvVVRU2IFAwP7jH/8Yhwzd46ttZdu2XVBQYN98881xycftDh8+bEuy161bZ9v2l/3I7/fbS5YsiUzz6aef2pLsDRs2xCtNV/hqW9m2bY8aNcr+8Y9/HL+kXOy8886zf/vb37b5PtVuj7DU19dry5Ytys3Njbzm8XiUm5urDRs2xDEz99q5c6eysrLUp08f3XHHHSotLY13Sq5XUlKisrKyqH6Wmpqq7Oxs+lkT1q5dq27duql///66++67dezYsXin5AqVlZWSpPPPP1+StGXLFgWDwai+NWDAAF144YXnfN/6alud8oc//EFdu3bVkCFDNGvWLJ04cSIe6blGKBTS4sWLVVNTo5ycnDbfp9rFaM2NOXr0qEKhkNLT06NeT09P12effRanrNwrOztbixYtUv/+/XXo0CE98sgjuvbaa7V9+3YlJyfHOz3XKisrk6RG+9mp9/C/xo4dq1tuuUW9e/fW7t279dOf/lT5+fnasGGDvF5vvNOLm3A4rHvuuUdXX321hgwZIunLvpWQkKDOnTtHTXuu963G2kqS/vVf/1U9e/ZUVlaWPv74Yz300EPasWOH3njjjThmGx+ffPKJcnJyVFtbq06dOunNN9/UoEGDVFxc3Kb7VLstWBCb/Pz8yN/Dhg1Tdna2evbsqddff13Tpk2LY2ZoT26//fbI30OHDtWwYcPUt29frV27VmPGjIljZvE1Y8YMbd++nfPGmqGptrrrrrsifw8dOlSZmZkaM2aMdu/erb59+37dacZV//79VVxcrMrKSv3nf/6nCgoKtG7dunin5Vi7/Umoa9eu8nq9p539XF5eroyMjDhl1XZ07txZF198sXbt2hXvVFztVF+in5np06ePunbtek73s5kzZ+qdd97RmjVr1L1798jrGRkZqq+vV0VFRdT053LfaqqtGpOdnS1J52TfSkhIUL9+/TRixAgVFRVp+PDheuGFF9p8n2q3BUtCQoJGjBihVatWRV4Lh8NatWqVcnJy4phZ21BdXa3du3crMzMz3qm4Wu/evZWRkRHVz6qqqrRp0yb6WTPs379fx44dOyf7mW3bmjlzpt58802tXr1avXv3jnp/xIgR8vv9UX1rx44dKi0tPef61tnaqjHFxcWSdE72ra8Kh8Oqq6tr+30q3mf9tqbFixfbgUDAXrRokf0///M/9l133WV37tzZLisri3dqrnPffffZa9eutUtKSuy//e1vdm5urt21a1f78OHD8U4t7o4fP25v27bN3rZtmy3JfvbZZ+1t27bZe/futW3btp944gm7c+fO9rJly+yPP/7Yvvnmm+3evXvbJ0+ejHPmX78ztdXx48ft+++/396wYYNdUlJir1y50r7sssvsiy66yK6trY136l+7u+++205NTbXXrl1rHzp0KPI4ceJEZJrvf//79oUXXmivXr3a3rx5s52Tk2Pn5OTEMev4OFtb7dq1y3700UftzZs32yUlJfayZcvsPn362Nddd12cM//6/eQnP7HXrVtnl5SU2B9//LH9k5/8xLYsy37vvfds227bfapdFyy2bdv/9m//Zl944YV2QkKCPXLkSHvjxo3xTsmVbrvtNjszM9NOSEiwL7jgAvu2226zd+3aFe+0XGHNmjW2pNMeBQUFtm1/eWnzz372Mzs9Pd0OBAL2mDFj7B07dsQ36Tg5U1udOHHCvvHGG+20tDTb7/fbPXv2tKdPn37O/gPRWDtJsl966aXINCdPnrR/8IMf2Oedd57doUMHe8KECfahQ4fil3ScnK2tSktL7euuu84+//zz7UAgYPfr189+4IEH7MrKyvgmHgd33nmn3bNnTzshIcFOS0uzx4wZEylWbLtt9ynLtm376zueAwAAELt2ew4LAABoPyhYAACA61GwAAAA16NgAQAArkfBAgAAXI+CBQAAuB4FCwAAcD0KFgAA4HoULAAAwPUoWAAAgOtRsAAAANf7/wB5pUwl6f0kqAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGMCAYAAAD5g/maAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA/FElEQVR4nO3deXhU5d3/8c+ZSWYSsrIEQmRXWQSBQmsaBUFJDakPigtVagsoYrXQiri0+HMBtcaligsU2qeP4lKL0kfBWqQiCmhZLGiqVOUhCASEsGl2ss2c3x8tU0eSkLlP4pyE9+u6znUxZ+7vnPvcc5+TL2fOzNeybdsWAACAi3mi3QEAAIATIWEBAACuR8ICAABcj4QFAAC4HgkLAABwPRIWAADgeiQsAADA9UhYAACA65GwAAAA1yNhAeAqvXr10pQpU5rc9r/+679atkPfEMuyNGPGjGh3A3AtEhYgyhYvXizLsrR58+Zod0WSVFlZqTlz5mjNmjXR7ook6eOPP9acOXO0a9euaHcFQBTFRLsDANylsrJSc+fOlSSNHj36G9/+tm3b5PH85/9SH3/8sebOnavRo0erV69e33h/ALgDCQsAV/H7/dHuAgAX4iMhwIWmTJmixMREff755xo/frwSExOVlpamW265RYFAINRu165dsixLv/71rzVv3jz17NlT8fHxGjVqlLZu3Rr2mqNHj673ismUKVNCVy527dqltLQ0SdLcuXNlWZYsy9KcOXMkSbW1tfr000+1f//+Rvv/6quvyrIsffjhh6F1//u//yvLsnTppZeGtR0wYICuuOKK0OOv3sOyePFiTZgwQZJ03nnnhfrz9Y+r3n33XZ111lmKi4tTnz599Oyzzzbav2OWLFmi4cOHKykpScnJyTrzzDP1+OOPh54/9nHdu+++q5///OdKS0tTamqqfvKTn6impkbFxcWaNGmS2rdvr/bt2+u2226Tbdth26ioqNDNN9+s7t27y+/3q1+/fvr1r399XLv63HffffJ4PHryySdD615//XWNHDlSCQkJSkpK0oUXXqh//vOfTdpfoDUjYQFcKhAIKCcnRx07dtSvf/1rjRo1So888oh+97vfHdf22Wef1RNPPKHp06dr9uzZ2rp1q84//3wdOHAgom2mpaVp4cKFkqRLLrlEzz33nJ577rlQkvH5559rwIABmj17dqOvM2LECFmWpXXr1oXWvfPOO/J4PHr33XdD6w4dOqRPP/1U5557br2vc+655+rnP/+5JOn2228P9WfAgAGhNgUFBbr88sv1ve99T4888ojat2+vKVOmnPCP+KpVqzRx4kS1b99eDz74oB544AGNHj1af/vb345r+7Of/Uzbt2/X3LlzddFFF+l3v/ud7rzzTo0bN06BQED333+/RowYoYcffljPPfdcKM62bV100UWaN2+exo4dq0cffVT9+vXTrbfeqlmzZjXavzvuuEN33XWXfvvb3+pnP/uZJOm5557ThRdeqMTERD344IO688479fHHH2vEiBHc44O2zwYQVU8//bQtyf773/8eWjd58mRbkn3PPfeEtf3Wt75lDx8+PPR4586dtiQ7Pj7e3rt3b2j9pk2bbEn2TTfdFFo3atQoe9SoUcdtf/LkyXbPnj1Djw8dOmRLsu++++7j2h7b3uTJk0+4XwMHDrR/8IMfhB4PGzbMnjBhgi3J/uSTT2zbtu2XX37ZlmT/4x//CLXr2bNn2OsvXbrUlmS//fbbx22jZ8+etiR73bp1oXUHDx60/X6/ffPNNzfavxtvvNFOTk626+rqGmxz7L3Jycmxg8FgaH1WVpZtWZZ9/fXXh9bV1dXZ3bp1CxvjZcuW2ZLs++67L+x1L7/8ctuyLLugoCC0TpI9ffp027Zt++abb7Y9Ho+9ePHi0PNlZWV2amqqPW3atLDXKioqslNSUo5bD7Q1XGEBXOz6668Pezxy5Eh99tlnx7UbP368TjnllNDjs846S5mZmVqxYkWz9qdXr16ybVuLFy8+YduRI0fqnXfekSSVlZXpH//4h6677jp16tQptP6dd95RamqqBg0aZNynM844QyNHjgw9TktLU79+/eodp69KTU1VRUWFVq1adcJtTJ06VZZlhR5nZmbKtm1NnTo1tM7r9erb3/522HZXrFghr9cbukp0zM033yzbtvX666+HrbdtWzNmzNDjjz+u559/XpMnTw49t2rVKhUXF2vixIk6fPhwaPF6vcrMzNTbb799wv0AWjMSFsCl4uLiQveTHNO+fXt9+eWXx7U9/fTTj1vXt2/fqH5MMHLkSO3fv18FBQVav369LMtSVlZWWCLzzjvv6Jxzzgn7VlCkevTocdy6hsbpq37605+qb9++ys3NVbdu3XTNNddo5cqVTdpGSkqKJKl79+7Hrf/qdnfv3q2MjAwlJSWFtTv2kdbu3bvD1j/77LNasGCBnnzySU2cODHsue3bt0uSzj//fKWlpYUtb7zxhg4ePNjo/gKtHd8SAlzK6/U26+tZllXvjZ5fvYm3OY0YMUKStG7dOn322WcaNmyYEhISNHLkSD3xxBMqLy/XBx98oF/96leOttPQONW3r1/VuXNn5efn669//atef/11vf7663r66ac1adIkPfPMM03aRn3rT7TdxpxzzjnKz8/X/Pnz9YMf/EAdOnQIPRcMBiX96z6W9PT042JjYjido23jCgvQBhz73/dX/d///V/Y75a0b99excXFx7X7+v/yv/rRhxM9evRQjx499M477+idd94JfWxz7rnnateuXVq6dKkCgUCDN9w2d3/q4/P5NG7cOP3mN7/Rjh079JOf/ETPPvusCgoKmuX1e/bsqX379qmsrCxs/aeffhp6/qtOO+00vfHGG9q3b5/Gjh0bFnfqqadK+leilZ2dfdwSjd/MAb5JJCxAG7Bs2TJ9/vnnocfvvfeeNm3apNzc3NC6U089VZ9++qkOHToUWvePf/zjuG/FtGvXTpLqTW6a+rXmY0aOHKm33npL7733XihhGTp0qJKSkvTAAw8oPj5ew4cPb/Q1EhISGuyPE0eOHAl77PF4NHjwYElSdXV1s2zj+9//vgKBgObPnx+2ft68ebIsK+z9OWbw4MFasWKFPvnkE40bN05Hjx6VJOXk5Cg5OVn333+/amtrj4v76vsKtEVcQwTagNNOO00jRozQDTfcoOrqaj322GPq2LGjbrvttlCba665Ro8++qhycnI0depUHTx4UIsWLdLAgQNVWloaahcfH68zzjhDL774ovr27asOHTpo0KBBGjRoUOhrzZMnT27yjbd/+MMfZFlW6CMir9ers88+W3/96181evRo+Xy+Rl9j6NCh8nq9evDBB1VSUiK/36/zzz9fnTt3Nhusf7v22mv1xRdf6Pzzz1e3bt20e/duPfnkkxo6dGjY16adGDdunM477zz9v//3/7Rr1y4NGTJEb7zxhpYvX66ZM2eGrpp83Xe/+10tX75c3//+93X55Zdr2bJlSk5O1sKFC/XjH/9Yw4YN05VXXqm0tDQVFhbqL3/5i84555zjEiOgLeEKC9AGTJo0ST/72c80f/58/epXv9LAgQP11ltvqWvXrqE2AwYM0LPPPquSkhLNmjVLr776qp577jkNGzbsuNf7/e9/r1NOOUU33XSTJk6cqD/96U9G/Tp2VaV///7q2LHjceu/+u2ehqSnp2vRokU6ePCgpk6dqokTJ+rjjz826s9X/ehHP1JcXJx+85vf6Kc//ameeeYZXXHFFXr99dcd3QT8VR6PR6+++qpmzpyp1157TTNnztTHH3+shx9+WI8++mijseeff75eeuklvfHGG/rxj3+sYDCoH/7wh1q9erVOOeUUPfzww7rxxhu1ZMkSDR06VFdffXWz9BlwK8t2cocYgKjatWuXevfurYcffli33HJLtLsDAC2GKywAAMD1SFgAAIDrkbAAAADX4x4WAADgelxhAQAArkfCAgAAXK9N/HBcMBjUvn37lJSU1KI/4w0AAJqPbdsqKytTRkbGCX//qE0kLPv27TuuaioAAGgd9uzZo27dujXapk0kLMdKt2eP/K5ZxVIHV2UCgTrjWFvmVXIry2uMYzukJJ24UQPqqiqNY33+OOPYo0Hz96j0qPl7VOEg1uPgfnafP9Y41ivzsbKdzGcHt+/HxJqfipx+a8DJ9w7qAkHjWK+DX9MNBM3PHU6uQsfGmM/LGAf76+Q9qg0eX3epqZKTzM9ZSQ6O4S9Ky41jPZ7GS1005lhFcBN20OzcUVdXp79t3hz6O96YNpGwHDsAY2JiFPsNJyweB59ABR38YYmJMT9hGY3Rv1kx3qhst85BwuJgs/Ka7648tvnBH+M177TXwXx28sffUcLi4E2KZsJiy8F77OAPuJOkw0msk/cpWgmLHTSPdXLOinWQhDsZZ4/HPNZZwmIcKqlp87LFbrpdsGCBevXqpbi4OGVmZuq9995rtP3SpUvVv39/xcXF6cwzz9SKFStaqmsAAKCVaZGE5cUXX9SsWbN099136/3339eQIUOUk5OjgwcP1tt+/fr1mjhxoqZOnaoPPvhA48eP1/jx47V169aW6B4AAGhlWiRhefTRRzVt2jRdffXVOuOMM7Ro0SK1a9dOTz31VL3tH3/8cY0dO1a33nqrBgwYoHvvvVfDhg2jVDoAAJDUAglLTU2NtmzZouzs7P9sxONRdna2NmzYUG/Mhg0bwtpLUk5OToPtq6urVVpaGrYAAIC2q9kTlsOHDysQCKhLly5h67t06aKioqJ6Y4qKiiJqn5eXp5SUlNDCV5oBAGjbWuUv3c6ePVslJSWhZc+ePdHuEgAAaEHN/rXmTp06yev16sCBA2HrDxw4oPT09Hpj0tPTI2rv9/vl9/ubp8MAAMD1mv0Ki8/n0/Dhw7V69erQumAwqNWrVysrK6vemKysrLD2krRq1aoG2wMAgJNLi/xw3KxZszR58mR9+9vf1llnnaXHHntMFRUVuvrqqyVJkyZN0imnnKK8vDxJ0o033qhRo0bpkUce0YUXXqglS5Zo8+bN+t3vftcS3QMAAK1MiyQsV1xxhQ4dOqS77rpLRUVFGjp0qFauXBm6sbawsDCsyNHZZ5+tF154QXfccYduv/12nX766Vq2bJkGDRrUEt0DAACtjGU7+c1jlygtLVVKSorGnjfiG/9pfie1hIKOaglVG8d2TEk2jq2rqjCO9fnjjWOd1BIqcVAPqLzSSS0h89+q9vvN64E4+mn+uijVEvK1zp/mr62Lzk/zR62WUOzJVUsoxUktoTjzsTpScnLVElq7caNKSkqUnNz436Y2UUvomKKiMnmdFH8x4HFyunRwB1GMgzpEFR7zZMeS+fge/rLMOLbGa36TdZ1lPtAe23x/OzahmFdDYn3mfQ7aDpLoWgdj5WBOeh3UXTH/0/0vzm7kM9/n2lrzAqZ1DnbaH2f+R9hJclh11Py8U1VlPlZ+v4O5VWN+LAUc1ANKjDd/j2prHRTzdFCHSAHDI8luelyr/FozAAA4uZCwAAAA1yNhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HokLAAAwPVIWAAAgOuRsAAAANcjYQEAAK5HwgIAAFyPhAUAALgeCQsAAHA9B7Wk3ScQCMqk3LvHY14ivs5jnvPVOagRH2uZb9dTEzSOra01L7duW+bTzRvjNY5tF2deqj3WNi/VHu/3Gcda5rurYMB8XgVinPwfxvw48njNt+s1f4v+xbzbjgQD5hsO2uaxMQ7G2nJwPBx1EFtcWmEcGx9nfhy2i29vHFtWaX6ujHHwl9lyMKF9lnlswHBeBe2mx3GFBQAAuB4JCwAAcD0SFgAA4HokLAAAwPVIWAAAgOuRsAAAANcjYQEAAK5HwgIAAFyPhAUAALgeCQsAAHA9EhYAAOB6zZ6w5OXl6Tvf+Y6SkpLUuXNnjR8/Xtu2bWs0ZvHixbIsK2yJc1D/BQAAtC3NnrCsXbtW06dP18aNG7Vq1SrV1tbqggsuUEVF48WrkpOTtX///tCye/fu5u4aAABopZq9WvPKlSvDHi9evFidO3fWli1bdO655zYYZ1mW0tPTm7s7AACgDWj2hOXrSkpKJEkdOnRotF15ebl69uypYDCoYcOG6f7779fAgQPrbVtdXa3q6urQ49LSUklS9/QUxRrU5bbtYMQxx9Sah6ouYF5u/YvicuPY8krzUu0pCeYf1aWndTSOjfH7jWN9PvM+1wUCxrHVNeaxlV+Z35EKmk8ryeM1DvV6W/x0Uv92Pc4uFNc6eI+dzI+ExGTj2PITXLFulIPzji0HJzwH59l28T7j2JjYWOPYoIOxqgma72+sbRnHes1D5XHwoYvpSNkRjHGL3nQbDAY1c+ZMnXPOORo0aFCD7fr166ennnpKy5cv1/PPP69gMKizzz5be/furbd9Xl6eUlJSQkv37t1bahcAAIALtGjCMn36dG3dulVLlixptF1WVpYmTZqkoUOHatSoUXr55ZeVlpam3/72t/W2nz17tkpKSkLLnj17WqL7AADAJVrsGu6MGTP02muvad26derWrVtEsbGxsfrWt76lgoKCep/3+/3yO/iIAAAAtC7NfoXFtm3NmDFDr7zyit566y317t074tcIBAL66KOP1LVr1+buHgAAaIWa/QrL9OnT9cILL2j58uVKSkpSUVGRJCklJUXx8fGSpEmTJumUU05RXl6eJOmee+7Rd7/7XZ122mkqLi7Www8/rN27d+vaa69t7u4BAIBWqNkTloULF0qSRo8eHbb+6aef1pQpUyRJhYWF8nzlzv4vv/xS06ZNU1FRkdq3b6/hw4dr/fr1OuOMM5q7ewAAoBVq9oTFtk/8FaU1a9aEPZ43b57mzZvX3F0BAABtBLWEAACA65GwAAAA1yNhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HrRqQffQhLbxSnWoJT4kS8Om28zKcE49ujROuPYQNC8rL3XQbn1dj7zHDfBZ173vFOHJONY8yLvUq0dbxx7uPSocazXQWl6j+WgRLyDeeX3+4xjTY7b5oiVpPg4r3mw+ZRWnD/OOLampsY4trS4zDi2rLTcONZrnfg3uhqS1iHFONbjdfBnrgm/K9YQJ8ewZL5dyzKfz0EHE9pjGcZ6mn6+4goLAABwPRIWAADgeiQsAADA9UhYAACA65GwAAAA1yNhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HokLAAAwPVIWAAAgOuRsAAAANcjYQEAAK7noO62+xytrlVdIPKy3F6vz3ibMTHmsRUVlcax1Qb7eUycg3e9OmhefrzWSelyr3lu7feZv0eqc7DdeAexicnGscFAwDi2orLCODa+XbxxbKeUROPY9E7msZKU1M5rHFtdXWsca9vm86O2tsY4dq/5YaiSklLj2No687HyxvjNYz3m58pA0PxYqq01319PwPxNsr3m+2s7mBuWZRZcV1fX5LZcYQEAAK5HwgIAAFyPhAUAALhesycsc+bMkWVZYUv//v0bjVm6dKn69++vuLg4nXnmmVqxYkVzdwsAALRiLXKFZeDAgdq/f39oeffddxtsu379ek2cOFFTp07VBx98oPHjx2v8+PHaunVrS3QNAAC0Qi2SsMTExCg9PT20dOrUqcG2jz/+uMaOHatbb71VAwYM0L333qthw4Zp/vz5LdE1AADQCrVIwrJ9+3ZlZGSoT58+uuqqq1RYWNhg2w0bNig7OztsXU5OjjZs2NBgTHV1tUpLS8MWAADQdjV7wpKZmanFixdr5cqVWrhwoXbu3KmRI0eqrKys3vZFRUXq0qVL2LouXbqoqKiowW3k5eUpJSUltHTv3r1Z9wEAALhLsycsubm5mjBhggYPHqycnBytWLFCxcXFeumll5ptG7Nnz1ZJSUlo2bNnT7O9NgAAcJ8W/6Xb1NRU9e3bVwUFBfU+n56ergMHDoStO3DggNLT0xt8Tb/fL7/f/JcPAQBA69Liv8NSXl6uHTt2qGvXrvU+n5WVpdWrV4etW7VqlbKyslq6awAAoJVo9oTllltu0dq1a7Vr1y6tX79el1xyibxeryZOnChJmjRpkmbPnh1qf+ONN2rlypV65JFH9Omnn2rOnDnavHmzZsyY0dxdAwAArVSzfyS0d+9eTZw4UUeOHFFaWppGjBihjRs3Ki0tTZJUWFgoj+c/edLZZ5+tF154QXfccYduv/12nX766Vq2bJkGDRrU3F0DAACtVLMnLEuWLGn0+TVr1hy3bsKECZowYUJzdwUAALQRLX7T7Tdpx+6D8nojLxXvoKK2Yg4fNY6trTMvXe6kNH1VddPLeX/dkaD5aAVUYhxbVmW+3bh2icax5VXmJeJj4+ONYxOTzftsOfig1+9grIIe8/eoNmB+LChQYx4rqbY6aBxbWW5+/NfUmZ9+q2rNx+vQl+a/W7XvsPkxXOPgvJPmt41j7YD5MRyQ+XYrqquMY+tqzMcqNtZ8XgUdHIfVNWbjHIhgmxQ/BAAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HokLAAAwPVIWAAAgOuRsAAAANcjYQEAAK5HwgIAAFyPhAUAALgeCQsAAHA9EhYAAOB6JCwAAMD1zOtQu1BV0JJXkZe4D3oijzkm3jYvP15jm5e193lijWMT4sxjq4LmfT4aMM+PrWrjUMUEzcvLy2P+/sZ6zPc3NtZnHFvjoKx9MGi+v3EO9red33x/HUxJSZLH9hrHtvP7jWNty/y8E/Saj1e7xETz2IQE41h5zY/Do1XmsZbMJ4jXbz43PA7eX4/XfLs1gYBxrOrMY2O9ZumEJ4K/2VxhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HokLAAAwPVIWAAAgOuRsAAAANcjYQEAAK5HwgIAAFyPhAUAALhesycsvXr1kmVZxy3Tp0+vt/3ixYuPaxsXF9fc3QIAAK1Ysxc//Pvf/67AV4ovbd26Vd/73vc0YcKEBmOSk5O1bdu20GPLQdEoAADQ9jR7wpKWlhb2+IEHHtCpp56qUaNGNRhjWZbS09ObuysAAKCNaPaE5atqamr0/PPPa9asWY1eNSkvL1fPnj0VDAY1bNgw3X///Ro4cGCD7aurq1VdXR16XFpaKkmK9XnlNSjLbdvmV3RqA+Zlzz0yLyEe5zfvc3y7WOPYmsoa49igg3GurjEfZ9vBOKcktzOOTfSZj3NcjINPa20H72/QvLx8Qrz5OPvNu6yjXzkXmLDrzPsdH2fe8fg48+3WHDXfZ4/M3+OUJPPjIWAfNY794kiJ+XYdnKOTEs3/RAYD5uN8tNo81htj3mevbX7e8XnNtmtF8HehRW+6XbZsmYqLizVlypQG2/Tr109PPfWUli9frueff17BYFBnn3229u7d22BMXl6eUlJSQkv37t1boPcAAMAtLNu27ZZ68ZycHPl8Pv35z39uckxtba0GDBigiRMn6t577623TX1XWLp3765+fc/8xq+wBIPm2btlm/8PKzHePItOTIg3ji1xcIXFE2u+3Vif+Y3YPgexTq6wpCSaxyalJBnHVtWaz8maqirj2PYJPvPYJL9xrGXXGcdKUpzBOeMYJ1dYgpb5dkscXGHZt++geeyBYuPYL0q5wtJU0bvCYhyqOMPtBgIB5W/LV0lJiZKTkxtt22IfCe3evVtvvvmmXn755YjiYmNj9a1vfUsFBQUNtvH7/fL7zU9wAACgdWmxj4Sefvppde7cWRdeeGFEcYFAQB999JG6du3aQj0DAACtTYskLMFgUE8//bQmT56smK9dJpo0aZJmz54denzPPffojTfe0Geffab3339fP/rRj7R7925de+21LdE1AADQCrXIR0JvvvmmCgsLdc011xz3XGFhoTye/+RJX375paZNm6aioiK1b99ew4cP1/r163XGGWe0RNcAAEAr1CIJywUXXKCG7uVds2ZN2ON58+Zp3rx5LdENAADQRlBLCAAAuB4JCwAAcD0SFgAA4HokLAAAwPVIWAAAgOuRsAAAANcjYQEAAK7XYrWEoiFYVyfLoJajL9Z8GOo85gXfnJSdbOc3L7pWV2dewDDGvE6kHNSYlOUgNtZnnpcnxJsXToxvZ17vyhfrpM/mRRetJPP9jYsJGsd6LfPYmBhndcW8lvlYWw4mZpyD97i43LygX0W1ebHIhIQE41g7YP4ep7XPMI6tM9+sCnZ9bhzrYErL5+C8Ewiab9gj8/lcY1hkMhBBkUiusAAAANcjYQEAAK5HwgIAAFyPhAUAALgeCQsAAHA9EhYAAOB6JCwAAMD1SFgAAIDrkbAAAADXI2EBAACuR8ICAABcj4QFAAC4HgkLAABwPRIWAADgejHR7kBzCgYCsuzI46xY82GwbIMN/ltdwLzMe7mD8vLVdebb9cX6jGMt2zw/Dnq9xrGB2hrj2IrqKuNYX4LfOLadml5y/ev8DsaqnT/WODYx3ny71TXm87muzvwYlKQYJ8e/x3zb7eLNx/qU2BTj2FiP+f7GBSuNYzsmdzaO7dA+yTg2EDAf5wef2m8c27mr+f4O6G0eu2rDJ8axFRVHjWNTEszGORBo+jHEFRYAAOB6JCwAAMD1SFgAAIDrkbAAAADXizhhWbduncaNG6eMjAxZlqVly5aFPW/btu666y517dpV8fHxys7O1vbt20/4ugsWLFCvXr0UFxenzMxMvffee5F2DQAAtFERJywVFRUaMmSIFixYUO/zDz30kJ544gktWrRImzZtUkJCgnJyclRV1fC3LV588UXNmjVLd999t95//30NGTJEOTk5OnjwYKTdAwAAbVDECUtubq7uu+8+XXLJJcc9Z9u2HnvsMd1xxx26+OKLNXjwYD377LPat2/fcVdivurRRx/VtGnTdPXVV+uMM87QokWL1K5dOz311FP1tq+urlZpaWnYAgAA2q5mvYdl586dKioqUnZ2dmhdSkqKMjMztWHDhnpjampqtGXLlrAYj8ej7OzsBmPy8vKUkpISWrp3796cuwEAAFymWROWoqIiSVKXLl3C1nfp0iX03NcdPnxYgUAgopjZs2erpKQktOzZs6cZeg8AANyqVf7Srd/vl99v/iuiAACgdWnWKyzp6emSpAMHDoStP3DgQOi5r+vUqZO8Xm9EMQAA4OTSrAlL7969lZ6ertWrV4fWlZaWatOmTcrKyqo3xufzafjw4WExwWBQq1evbjAGAACcXCL+SKi8vFwFBQWhxzt37lR+fr46dOigHj16aObMmbrvvvt0+umnq3fv3rrzzjuVkZGh8ePHh2LGjBmjSy65RDNmzJAkzZo1S5MnT9a3v/1tnXXWWXrsscdUUVGhq6++2vkeAgCAVi/ihGXz5s0677zzQo9nzZolSZo8ebIWL16s2267TRUVFbruuutUXFysESNGaOXKlYqLiwvF7NixQ4cPHw49vuKKK3To0CHdddddKioq0tChQ7Vy5crjbsQFAAAnJ8u2bWe12V2gtLRUKSkpumhUpmJjIr+P2OuxjLcdsIPGsQeLq41jq6rrjGPLjtYYx57evb1xbIcEn3Gs3xd34kYNsLxe41hfXLxxbFxCO+PYBAc3lXst809628WZlYiXpIT2icaxgaD5MVhbXWscK0l1AfNjuLzS/Djs1cV8bqV3TjWOrSwzP+/s/+yfxrFWXYVxbHyc+fHfs9dpxrFPvLTGOLZr967GsecMOdU49u//3G0cW1Cw0zg2MdFsPtfW1enlNzeppKREycnJjballhAAAHA9EhYAAOB6JCwAAMD1SFgAAIDrkbAAAADXI2EBAACuR8ICAABcj4QFAAC4HgkLAABwPRIWAADgeiQsAADA9UhYAACA65GwAAAA1yNhAQAArkfCAgAAXC8m2h1oTkWHSxTj9UYcF+f3GW+zurraONa2zIffa7CfIZZ5aFlljXFsXU2tcWxcrPl2fbGx5rE+8z4nVJnPjdiUBONYy2P+/5DaKvN5FeM3n89JyYnGsXat+dyQJK/HfJ/tBAfjFaw0jvXVmh/ENXbAOLZw3wHj2KqKYuPYjK6djGPjD5v3OSHO/P1Na59iHBtj2caxQ/tmGMfaZZ8bx9YEgmZxtU2P4woLAABwPRIWAADgeiQsAADA9UhYAACA65GwAAAA1yNhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HoRJyzr1q3TuHHjlJGRIcuytGzZstBztbW1+sUvfqEzzzxTCQkJysjI0KRJk7Rv375GX3POnDmyLCts6d+/f8Q7AwAA2qaIE5aKigoNGTJECxYsOO65yspKvf/++7rzzjv1/vvv6+WXX9a2bdt00UUXnfB1Bw4cqP3794eWd999N9KuAQCANiriamW5ubnKzc2t97mUlBStWrUqbN38+fN11llnqbCwUD169Gi4IzExSk9Pj7Q7AADgJNDi97CUlJTIsiylpqY22m779u3KyMhQnz59dNVVV6mwsLDBttXV1SotLQ1bAABA22VeD74Jqqqq9Itf/EITJ05UcnJyg+0yMzO1ePFi9evXT/v379fcuXM1cuRIbd26VUlJSce1z8vL09y5c49bH/AkSJ7Id6miNuKQkENHzJOl5GTzEvGWzEvEV1RWG8cesM3Lrft8scaxcbHm5dYTk3zGsVat+XZjqs3fo6IK8/fI46A0vc9n/n+Y0+MTjWMV4zcOPVplVtb+GI9VZxyb6Dc/Hmory41jy7zmfba95sdhWpcOxrE+byfj2I6pCcaxR2vMx6qq2nxuVVYeNY6tqTPv8yEH/3//21Ynfa4xigsEm36ebLErLLW1tfrBD34g27a1cOHCRtvm5uZqwoQJGjx4sHJycrRixQoVFxfrpZdeqrf97NmzVVJSElr27NnTErsAAABcokWusBxLVnbv3q233nqr0asr9UlNTVXfvn1VUFBQ7/N+v19+v/n/yAAAQOvS7FdYjiUr27dv15tvvqmOHTtG/Brl5eXasWOHunbt2tzdAwAArVDECUt5ebny8/OVn58vSdq5c6fy8/NVWFio2tpaXX755dq8ebP+8Ic/KBAIqKioSEVFRaqp+c/nW2PGjNH8+fNDj2+55RatXbtWu3bt0vr163XJJZfI6/Vq4sSJzvcQAAC0ehF/JLR582add955ocezZs2SJE2ePFlz5szRq6++KkkaOnRoWNzbb7+t0aNHS5J27Nihw4cPh57bu3evJk6cqCNHjigtLU0jRozQxo0blZaWFmn3AABAGxRxwjJ69GjZdsPfQmjsuWN27doV9njJkiWRdgMAAJxEqCUEAABcj4QFAAC4HgkLAABwPRIWAADgeiQsAADA9UhYAACA65GwAAAA1yNhAQAArtcixQ+jJRgMyLIij6upqTXeZlXAvPx4kmWeL8Z4DHb03+JifcaxKe3Mi04mJ5rHeh2k1p0c9Nl2MFaHKiqNY4urq4xj/T6vcazHaz5WxV8WG8cmxJmPs+3gOJKkutqml7f/uvJq8/fY8pn3u/xonXFsYko749jOHVKNY4+Wfmkcm9Y+0Tg26DHfX7/3M+PYsi9LjGNLjnYyjt2974hxbHmt+d/COsvsvBPJX1CusAAAANcjYQEAAK5HwgIAAFyPhAUAALgeCQsAAHA9EhYAAOB6JCwAAMD1SFgAAIDrkbAAAADXI2EBAACuR8ICAABcj4QFAAC4HgkLAABwPRIWAADgejHR7kBzKqk4Kq8n8hLXgTrzUu1er1lJbUmqqKoxjo2VZRzri/UZx9ZUm/e5zhNJIfGvbTdoG8f6PbHGsSnt4oxjd5ZXGsdW1JjPyRgHc7Kqyny7tTXVxrEJ7czfo3YJCcaxkhS0zccr4OAU6rPMj+FK80NJVq157P7D5nP64L4vjWPj28Ubx57ap4NxrN9n/v564s37fLjC/Fjad9h8nKvrzM+zAQWM4oLBpsdxhQUAALgeCQsAAHA9EhYAAOB6EScs69at07hx45SRkSHLsrRs2bKw56dMmSLLssKWsWPHnvB1FyxYoF69eikuLk6ZmZl67733Iu0aAABooyJOWCoqKjRkyBAtWLCgwTZjx47V/v37Q8sf//jHRl/zxRdf1KxZs3T33Xfr/fff15AhQ5STk6ODBw9G2j0AANAGRXwLdG5urnJzcxtt4/f7lZ6e3uTXfPTRRzVt2jRdffXVkqRFixbpL3/5i5566in98pe/jLSLAACgjWmRe1jWrFmjzp07q1+/frrhhht05MiRBtvW1NRoy5Ytys7O/k+nPB5lZ2drw4YN9cZUV1ertLQ0bAEAAG1XsycsY8eO1bPPPqvVq1frwQcf1Nq1a5Wbm6tAoP7vWh8+fFiBQEBdunQJW9+lSxcVFRXVG5OXl6eUlJTQ0r179+beDQAA4CLN/sNxV155ZejfZ555pgYPHqxTTz1Va9as0ZgxY5plG7Nnz9asWbNCj0tLS0laAABow1r8a819+vRRp06dVFBQUO/znTp1ktfr1YEDB8LWHzhwoMH7YPx+v5KTk8MWAADQdrV4wrJ3714dOXJEXbt2rfd5n8+n4cOHa/Xq1aF1wWBQq1evVlZWVkt3DwAAtAIRJyzl5eXKz89Xfn6+JGnnzp3Kz89XYWGhysvLdeutt2rjxo3atWuXVq9erYsvvlinnXaacnJyQq8xZswYzZ8/P/R41qxZ+u///m8988wz+uSTT3TDDTeooqIi9K0hAABwcov4HpbNmzfrvPPOCz0+di/J5MmTtXDhQn344Yd65plnVFxcrIyMDF1wwQW699575ff7QzE7duzQ4cOHQ4+vuOIKHTp0SHfddZeKioo0dOhQrVy58rgbcQEAwMkp4oRl9OjRsu2GKzr+9a9/PeFr7Nq167h1M2bM0IwZMyLtDgAAOAk0+7eEosnvseX1Rl4eu1bmZd593ljj2JqaOuPYiqB5GfC4OPP9TXVQMr2yyrzPQeNIqaLGvFR7YlWVcWx1rXmvLQe3lwUC5tsNBszfo4SEROPYI6UVxrFlNWZl7Y+xPOanwfZJCcaxpTXmY11Za37uCDjY38PFlcaxFdXm+1tcar6/hfsOGceWOziWir8sN45VXDvj0GoHcyMoB8eS6dtrN32bFD8EAACuR8ICAABcj4QFAAC4HgkLAABwPRIWAADgeiQsAADA9UhYAACA65GwAAAA1yNhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HokLAAAwPXM64y7UPt2sYrxeiOOK62qNd6mpci3d4xdZ14G3BNjnmvGOYhN8hmH6qh5qIKWZRwb8JiXtT/q4D3yWg5iY8znlddj/v76fUHj2Hif+XtUY34I6kjxYfNgSSnJycaxsQ62GxtrPl4K1hiHlpaWG8d+fvAL49ienVKMYxNSOxnHHio1P/NU1ATMt1v0pXHs4cMVxrG1teYHU4xlvr+WZXbOiuSMwxUWAADgeiQsAADA9UhYAACA65GwAAAA1yNhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HokLAAAwPUiTljWrVuncePGKSMjQ5ZladmyZWHPW5ZV7/Lwww83+Jpz5sw5rn3//v0j3hkAANA2RZywVFRUaMiQIVqwYEG9z+/fvz9seeqpp2RZli677LJGX3fgwIFhce+++26kXQMAAG1UxMUPc3NzlZub2+Dz6enpYY+XL1+u8847T3369Gm8IzExx8UCAABILXwPy4EDB/SXv/xFU6dOPWHb7du3KyMjQ3369NFVV12lwsLCBttWV1ertLQ0bAEAAG1XxFdYIvHMM88oKSlJl156aaPtMjMztXjxYvXr10/79+/X3LlzNXLkSG3dulVJSUnHtc/Ly9PcuXOPW2/FeGV5Iy9x7fXURRxzjO0xLxFv27ZxrC/WrJS3JMV4zftcaV59XDE+8+lmOUitvwiaj3NpZaVxbEK7WOPYQDCSouvhfA7GuS5oXpr+s8L9xrGJ7doZx3pjzI8FSaqpqTaOraw0n1sJCeb7bAfNj+GyCvPzndcfbxybfkoX49jk5ATj2Div+ckjUGN+HH5ZWm4cG1Nh3ufkBL9xbGo789gYr9l5JxAIaGcT27boFZannnpKV111leLi4hptl5ubqwkTJmjw4MHKycnRihUrVFxcrJdeeqne9rNnz1ZJSUlo2bNnT0t0HwAAuESLXWF55513tG3bNr344osRx6ampqpv374qKCio93m/3y+/3zwTBAAArUuLXWH5n//5Hw0fPlxDhgyJOLa8vFw7duxQ165dW6BnAACgtYk4YSkvL1d+fr7y8/MlSTt37lR+fn7YTbKlpaVaunSprr322npfY8yYMZo/f37o8S233KK1a9dq165dWr9+vS655BJ5vV5NnDgx0u4BAIA2KOKPhDZv3qzzzjsv9HjWrFmSpMmTJ2vx4sWSpCVLlsi27QYTjh07dujw4cOhx3v37tXEiRN15MgRpaWlacSIEdq4caPS0tIi7R4AAGiDIk5YRo8efcJvt1x33XW67rrrGnx+165dYY+XLFkSaTcAAMBJhFpCAADA9UhYAACA65GwAAAA1yNhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HotVvwwGqprAwoYVAOvCwSMt1lXZ15+POig7LnHMo+1rMZ/+K8x5VVVxrHeGPPpFusgtsbjYJw95mNlyTKODZzgxxkbYwdrjWM9lvl89nq8xrGVFeXGsT5/rHGsJFkyP/4VqDEOtYPmY330qPl2A0HzuVVbaz63CosOGcfG1Jmfd1IDPuPYL78oNY6tqqk2jo2PNe9zaZX53CivNB9n0z8rwQiOA66wAAAA1yNhAQAArkfCAgAAXI+EBQAAuB4JCwAAcD0SFgAA4HokLAAAwPVIWAAAgOuRsAAAANcjYQEAAK5HwgIAAFyPhAUAALgeCQsAAHA9EhYAAOB6MdHuQHMKBmwFFHmNa6/XfBi8xpFSjYOq9rUB89L0qjEvLx8X62CsLONQKWg+WHawzny7lnmnLYO5eIyTeVXtYLuWbb5ljxy8Rw5iKyuOGsdKcrBlyQ442Gfzt0lOZojHY37u8FrmsSVflhnHfvqx+TFcVmH+Hn1x1Dw2IT7ROLa2uto4NlBlPlbxsT7j2KOGfbYjOBC4wgIAAFyPhAUAALgeCQsAAHC9iBKWvLw8fec731FSUpI6d+6s8ePHa9u2bWFtqqqqNH36dHXs2FGJiYm67LLLdODAgUZf17Zt3XXXXeratavi4+OVnZ2t7du3R743AACgTYooYVm7dq2mT5+ujRs3atWqVaqtrdUFF1ygioqKUJubbrpJf/7zn7V06VKtXbtW+/bt06WXXtro6z700EN64okntGjRIm3atEkJCQnKyclRVVWV2V4BAIA2xbIjuUX3aw4dOqTOnTtr7dq1Ovfcc1VSUqK0tDS98MILuvzyyyVJn376qQYMGKANGzbou9/97nGvYdu2MjIydPPNN+uWW26RJJWUlKhLly5avHixrrzyyhP2o7S0VCkpKRrW/wx5vZHfPR8IOvjGjQNl1eZ3oDu40V+xMdH5llBMrINvNjj4to7t5Ns6DrbrcbBdJ18gqYnat4TMP2F28i0hZ6PFt4Qi2qqDbwklxfuNY2ProvMtoX1Ha41jvTHm37hx8i0heczPWV6v+TFs+i2hYDCofUWFKikpUXJycqNtHd3DUlJSIknq0KGDJGnLli2qra1VdnZ2qE3//v3Vo0cPbdiwod7X2Llzp4qKisJiUlJSlJmZ2WBMdXW1SktLwxYAANB2GScswWBQM2fO1DnnnKNBgwZJkoqKiuTz+ZSamhrWtkuXLioqKqr3dY6t79KlS5Nj8vLylJKSElq6d+9uuhsAAKAVME5Ypk+frq1bt2rJkiXN2Z8mmT17tkpKSkLLnj17vvE+AACAb45RwjJjxgy99tprevvtt9WtW7fQ+vT0dNXU1Ki4uDis/YEDB5Senl7vax1b//VvEjUW4/f7lZycHLYAAIC2K6KExbZtzZgxQ6+88oreeust9e7dO+z54cOHKzY2VqtXrw6t27ZtmwoLC5WVlVXva/bu3Vvp6elhMaWlpdq0aVODMQAA4OQSUcIyffp0Pf/883rhhReUlJSkoqIiFRUV6ejRf9XySElJ0dSpUzVr1iy9/fbb2rJli66++mplZWWFfUOof//+euWVVyRJlmVp5syZuu+++/Tqq6/qo48+0qRJk5SRkaHx48c3354CAIBWK6LvqC5cuFCSNHr06LD1Tz/9tKZMmSJJmjdvnjwejy677DJVV1crJydHv/nNb8Lab9u2LfQNI0m67bbbVFFRoeuuu07FxcUaMWKEVq5cqbi4OINdAgAAbY2j32FxC36HJTL8DkvT8TssTcfvsEQYy++wNBm/wxKBNvw7LOZ/fVyoJhAwOozjfLHG24yxzM867RwcwPuPmJdqL6swP/gtq8Y41sn52Zb5Qejs75l5sMfByV0O9rfOQZJl2Q5Odg5iLZmPle0g9l/x5n/8nbxPDma0LEeT2ny8LEdzutI4MmiZv0d1AQfnOwfvks98s4rxmP9p9jpISH0x5vsbExNvFBeIIOmn+CEAAHA9EhYAAOB6JCwAAMD1SFgAAIDrkbAAAADXI2EBAACuR8ICAABcj4QFAAC4HgkLAABwPRIWAADgeiQsAADA9UhYAACA65GwAAAA12sT1Zrtf9dpj6Tq41fVBRzkbQ6qNTupBBoMmpdMDwadVGs1Dj3pqjUrStWag1Gq1uwoNqrVmp3VTf7mI1trtWZzTua0o/Odo+PQwTnaOFKOzjuBgPm8ChqOVeDff8uO/R1vTJtIWMrKyiRJn3xWEOWeAACASJWVlSklJaXRNpbdlLTG5YLBoPbt26ekpCRZ9WTipaWl6t69u/bs2aPk5OQo9LD1YKyajrFqOsYqMoxX0zFWTefGsbJtW2VlZcrIyJDH0/inHW3iCovH41G3bt1O2C45Odk1b5LbMVZNx1g1HWMVGcar6RirpnPbWJ3oysox3HQLAABcj4QFAAC43kmRsPj9ft19993y+/3R7orrMVZNx1g1HWMVGcar6RirpmvtY9UmbroFAABt20lxhQUAALRuJCwAAMD1SFgAAIDrkbAAAADXI2EBAACu1+YTlgULFqhXr16Ki4tTZmam3nvvvWh3yZXmzJkjy7LClv79+0e7W66wbt06jRs3ThkZGbIsS8uWLQt73rZt3XXXXeratavi4+OVnZ2t7du3R6ezUXaisZoyZcpx82zs2LHR6WyU5eXl6Tvf+Y6SkpLUuXNnjR8/Xtu2bQtrU1VVpenTp6tjx45KTEzUZZddpgMHDkSpx9HTlLEaPXr0cXPr+uuvj1KPo2fhwoUaPHhw6Ndss7Ky9Prrr4eeb81zqk0nLC+++KJmzZqlu+++W++//76GDBminJwcHTx4MNpdc6WBAwdq//79oeXdd9+NdpdcoaKiQkOGDNGCBQvqff6hhx7SE088oUWLFmnTpk1KSEhQTk6OqqqqvuGeRt+JxkqSxo4dGzbP/vjHP36DPXSPtWvXavr06dq4caNWrVql2tpaXXDBBaqoqAi1uemmm/TnP/9ZS5cu1dq1a7Vv3z5deumlUex1dDRlrCRp2rRpYXProYceilKPo6dbt2564IEHtGXLFm3evFnnn3++Lr74Yv3zn/+U1MrnlN2GnXXWWfb06dNDjwOBgJ2RkWHn5eVFsVfudPfdd9tDhgyJdjdcT5L9yiuvhB4Hg0E7PT3dfvjhh0PriouLbb/fb//xj3+MQg/d4+tjZdu2PXnyZPviiy+OSn/c7uDBg7Yke+3atbZt/2sexcbG2kuXLg21+eSTT2xJ9oYNG6LVTVf4+ljZtm2PGjXKvvHGG6PXKRdr3769/fvf/77Vz6k2e4WlpqZGW7ZsUXZ2dmidx+NRdna2NmzYEMWeudf27duVkZGhPn366KqrrlJhYWG0u+R6O3fuVFFRUdg8S0lJUWZmJvOsAWvWrFHnzp3Vr18/3XDDDTpy5Ei0u+QKJSUlkqQOHTpIkrZs2aLa2tqwudW/f3/16NHjpJ9bXx+rY/7whz+oU6dOGjRokGbPnq3KyspodM81AoGAlixZooqKCmVlZbX6OdUmqjXX5/DhwwoEAurSpUvY+i5duujTTz+NUq/cKzMzU4sXL1a/fv20f/9+zZ07VyNHjtTWrVuVlJQU7e65VlFRkSTVO8+OPYf/GDt2rC699FL17t1bO3bs0O23367c3Fxt2LBBXq832t2LmmAwqJkzZ+qcc87RoEGDJP1rbvl8PqWmpoa1PdnnVn1jJUk//OEP1bNnT2VkZOjDDz/UL37xC23btk0vv/xyFHsbHR999JGysrJUVVWlxMREvfLKKzrjjDOUn5/fqudUm01YEJnc3NzQvwcPHqzMzEz17NlTL730kqZOnRrFnqEtufLKK0P/PvPMMzV48GCdeuqpWrNmjcaMGRPFnkXX9OnTtXXrVu4ba4KGxuq6664L/fvMM89U165dNWbMGO3YsUOnnnrqN93NqOrXr5/y8/NVUlKiP/3pT5o8ebLWrl0b7W451mY/EurUqZO8Xu9xdz8fOHBA6enpUepV65Gamqq+ffuqoKAg2l1xtWNziXlmpk+fPurUqdNJPc9mzJih1157TW+//ba6desWWp+enq6amhoVFxeHtT+Z51ZDY1WfzMxMSTop55bP59Npp52m4cOHKy8vT0OGDNHjjz/e6udUm01YfD6fhg8frtWrV4fWBYNBrV69WllZWVHsWetQXl6uHTt2qGvXrtHuiqv17t1b6enpYfOstLRUmzZtYp41wd69e3XkyJGTcp7Ztq0ZM2bolVde0VtvvaXevXuHPT98+HDFxsaGza1t27apsLDwpJtbJxqr+uTn50vSSTm3vi4YDKq6urr1z6lo3/XbkpYsWWL7/X578eLF9scff2xfd911dmpqql1UVBTtrrnOzTffbK9Zs8beuXOn/be//c3Ozs62O3XqZB88eDDaXYu6srIy+4MPPrA/+OADW5L96KOP2h988IG9e/du27Zt+4EHHrBTU1Pt5cuX2x9++KF98cUX271797aPHj0a5Z5/8xobq7KyMvuWW26xN2zYYO/cudN+88037WHDhtmnn366XVVVFe2uf+NuuOEGOyUlxV6zZo29f//+0FJZWRlqc/3119s9evSw33rrLXvz5s12VlaWnZWVFcVeR8eJxqqgoMC+55577M2bN9s7d+60ly9fbvfp08c+99xzo9zzb94vf/lLe+3atfbOnTvtDz/80P7lL39pW5Zlv/HGG7Ztt+451aYTFtu27SeffNLu0aOH7fP57LPOOsveuHFjtLvkSldccYXdtWtX2+fz2aeccop9xRVX2AUFBdHuliu8/fbbtqTjlsmTJ9u2/a+vNt955512ly5dbL/fb48ZM8betm1bdDsdJY2NVWVlpX3BBRfYaWlpdmxsrN2zZ0972rRpJ+1/IOobJ0n2008/HWpz9OhR+6c//andvn17u127dvYll1xi79+/P3qdjpITjVVhYaF97rnn2h06dLD9fr992mmn2bfeeqtdUlIS3Y5HwTXXXGP37NnT9vl8dlpamj1mzJhQsmLbrXtOWbZt29/c9RwAAIDItdl7WAAAQNtBwgIAAFyPhAUAALgeCQsAAHA9EhYAAOB6JCwAAMD1SFgAAIDrkbAAAADXI2EBAACuR8ICAABcj4QFAAC43v8Hf8UTvkBRnsEAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "solve_haze_detection 33x22, 726 pixels: solve_haze_detection 33x22, 726 pixels: 4.1 seconds, 257.9 seconds CPU\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Clipping input data to the valid range for imshow with RGB data ([0..1] for floats or [0..255] for integers). Got range [1.0112319538263253e-31..95.76781607305696].\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGMCAYAAAD5g/maAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA9uklEQVR4nO3df3wU1b3/8ffsZrOJmAQUyA9Ffolg+WVFjaggltRAWxW1CpQWsCithVaLWMVeBaptqN5ataC291axUn/2ClaLtIoItYAWkCJtpUADga8EhGsSEkiyyZ7vH162XZNA9pyEnYTX8/HYxyM7O5+dz5w5O/vJ7MwczxhjBAAA4GOBZCcAAABwLBQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfI+CBQAA+B4FCwAA8D0KFgAA4HsULADQwnbs2CHP8/Sf//mfyU4FaDcoWIA26q9//au++tWv6rTTTlM4HFZeXp4mTJigv/71r07v+6Mf/UhLlixpmSSPYfXq1ZozZ47Kysqc3qe2tlYPP/ywPvvZzyozM1MdO3ZU//79NXXqVH3wwQctkyyApKJgAdqgl156Seeee66WL1+uG264QY8++qimTJmiFStW6Nxzz9XixYut3/t4Fyxz5851LliuvfZa3XbbbRowYIDmzZunuXPnavjw4Xrttde0du3alkkWQFKlJDsBAInZvn27vva1r6lXr15atWqVunTpEnvtlltu0bBhw/S1r31NmzZtUq9evZKY6fHx5z//Wa+++qp++MMf6q677op7bf78+c7FEAB/4AgL0MY88MADOnTokH7xi1/EFSuS1LlzZ/385z9XVVWV7r///tj0yZMnq0ePHg3ea86cOfI8L/bc8zxVVVXpqaeekud58jxPkydPjpv3gw8+0PXXX6/MzEydeuqpuuWWW1RdXR17jyPnbyxcuLDB8jzP05w5c2Lvd/vtt0uSevbsGVvejh07JEn79+/XBx98oEOHDh21PbZv3y5Juvjiixu8FgwGdeqppzZY33/84x/66le/qqysLHXp0kV33323jDHatWuXrrrqKmVmZionJ0c/+clPGrznvn37NGXKFGVnZystLU2DBw/WU089ddQcJckYo6lTpyo1NVUvvfRSbPqiRYs0ZMgQpaen65RTTtG4ceO0a9euY74fcKKhYAHamFdeeUU9evTQsGHDGn19+PDh6tGjh373u98l/N5PP/20wuGwhg0bpqefflpPP/20vvGNb8TNc/3116u6ulpFRUX6whe+oEceeURTp05NeFnXXHONxo8fL0n66U9/GlvekSJs/vz5Ovvss/Xuu+8e9X26d+8uSfr1r3+turq6Zi177NixikajmjdvnvLz83XffffpoYce0uc//3mddtpp+vGPf6wzzzxTM2fO1KpVq2Jxhw8f1ogRI/T0009rwoQJeuCBB5SVlaXJkyfr4YcfbnJ59fX1mjx5sn71q19p8eLFuuaaayRJP/zhDzVx4kT16dNHDz74oG699VYtX75cw4cP58gQ8GkGQJtRVlZmJJmrrrrqqPNdeeWVRpKpqKgwxhgzadIk07179wbzzZ4923x6N9ChQwczadKkJue98sor46Z/61vfMpLMX/7yF2OMMcXFxUaSefLJJxu8hyQze/bs2PMHHnjASDLFxcVNLm/FihVHXddoNGouvfRSI8lkZ2eb8ePHmwULFpidO3c2+Z5Tp06NTaurqzOnn3668TzPzJs3Lzb9448/Nunp6XFt8dBDDxlJZtGiRbFptbW1ZujQoebkk0+OtfeRNnjggQdMJBIxY8eONenp6eb3v/99LG7Hjh0mGAyaH/7wh3E5vv/++yYlJaXBdOBExxEWoA05ePCgJCkjI+Oo8x15vaKiosVzmDZtWtzzb3/725KkpUuXtuhy5syZI2OMRowYcdT5PM/T73//e913333q1KmTnn32WU2bNk3du3fX2LFjGz1SceONN8b+DgaDOu+882SM0ZQpU2LTO3bsqL59++qf//xnbNrSpUuVk5MTOzIkSaFQSN/5zndUWVmplStXxi2ntrZW1113nV599VUtXbpUl19+eey1l156SdFoVNdff732798fe+Tk5KhPnz5asWJFc5sKOCFw0i3QhhwpRI4ULk1pbmFjo0+fPnHPe/furUAgEDv3JBnC4bC+//3v6/vf/7727NmjlStX6uGHH9YLL7ygUCikRYsWxc1/xhlnxD3PyspSWlqaOnfu3GD6gQMHYs937typPn36KBCI/1/v7LPPjr3+74qKilRZWanXXnutQeG1detWGWMatOcRoVDo2CsOnEAoWIA2JCsrS7m5udq0adNR59u0aZNOO+00ZWZmSlLcibX/rr6+3jmnT793ay6rOXJzczVu3Dhde+216t+/v1544QUtXLhQKSn/2t0Fg8EGcY1Nkz45WdZWYWGhli1bpvvvv18jRoxQWlpa7LVoNCrP8/Taa681uuyTTz7ZerlAe8RPQkAb86UvfUnFxcV6++23G339j3/8o3bs2KEvfelLsWmdOnVq9KeRTx8RkJouOI7YunVr3PNt27YpGo3GrkLq1KmTJDVYns2yXIRCIQ0aNEiRSET79+9vkffs3r27tm7dqmg0Gjf9yM3pjpwAfMSFF16oJUuWaPXq1bruuuviTgru3bu3jDHq2bOnCgoKGjwuvPDCFskZaC8oWIA25vbbb1d6erq+8Y1vxP1cIUn/+7//q29+85s66aSTYpcMS598OZaXl8cdmdmzZ0+jN5jr0KHDUa9QWbBgQdzzn/3sZ5Kk0aNHS5IyMzPVuXPnuKtrJOnRRx9tdFlSw+JGav5lzVu3blVJSUmD6WVlZVqzZo06derU4PJvW1/4whdUWlqq559/Pjatrq5OP/vZz3TyySfr0ksvbRBTUFCg5557TsuWLdPXvva1WLFzzTXXKBgMau7cuQ2O4hhjGmxb4ETHT0JAG9OnTx899dRTmjBhggYOHKgpU6aoZ8+e2rFjh375y19q//79evbZZ9W7d+9YzLhx43THHXfo6quv1ne+8x0dOnRIjz32mM466yxt2LAh7v2HDBmiN954Qw8++KDy8vLUs2dP5efnx14vLi7WlVdeqVGjRmnNmjVatGiRvvKVr2jw4MGxeW688UbNmzdPN954o8477zytWrVK//jHPxqsy5AhQyRJ3//+9zVu3DiFQiFdccUV6tChg+bPn6+5c+dqxYoVRz3x9i9/+Yu+8pWvaPTo0Ro2bJhOOeUU/b//9//01FNP6cMPP9RDDz3U5M89iZo6dap+/vOfa/LkyVq/fr169Oih3/zmN/rTn/6khx56qMlzhsaMGaMnn3xSEydOVGZmpn7+85+rd+/euu+++zRr1izt2LFDY8aMUUZGhoqLi7V48WJNnTpVM2fObJG8gXYhiVcoAXCwadMmM378eJObm2tCoZDJyckx48ePN++//36j8//hD38wAwYMMKmpqaZv375m0aJFjV7W/MEHH5jhw4eb9PR0Iyl2We+Ref/2t7+ZL3/5yyYjI8N06tTJTJ8+3Rw+fDjuPQ4dOmSmTJlisrKyTEZGhrn++uvNvn37GlzWbIwx9957rznttNNMIBCIu8S5uZc1792718ybN89ceumlJjc316SkpJhOnTqZz33uc+Y3v/lN3LxH3vOjjz6Kmz5p0iTToUOHBu996aWXmv79+zdY3g033GA6d+5sUlNTzcCBAxtcwv3vlzX/u0cffdRIMjNnzoxN+5//+R9zySWXmA4dOpgOHTqYfv36mWnTppktW7Ycdb2BE41njMMZZQBOGHPmzNHcuXP10UcfNbiaBgBaG+ewAAAA36NgAQAAvkfBAgAAfI9zWAAAgO9xhAUAAPgeBQsAAPC9dnHjuGg0qg8//FAZGRmteqtvAADQcowxOnjwoPLy8hoMKvpp7aJg+fDDD9WtW7dkpwEAACzs2rVLp59++lHnaRcFy5HbYXfMzLI8wmJ/VMblgI6R/fnOnkPOKQGXpO1zNg45Rx1+vYw6LNfpjHSn89kdYh2W67lsX5fVTeKBUZe82+I6u+Tssr9L1sFvl00UdDhpwm03m6SO5bbDs4syRh9XlDc5rMW/axcFy5EixfM8BdpUwWLPpWA51mG3o4o6fKE57bHsc3YplJLxAU5mrJesnJP6U26SLpRsg79eu/zk7raJk7ONXHIOOFQs5lOjgSe4ZIfY49/OR9a0OX2r1U66XbBggXr06KG0tDTl5+fr3XffPer8L774ovr166e0tDQNHDhQS5cuba3UAABAG9MqBcvzzz+vGTNmaPbs2dqwYYMGDx6swsJC7du3r9H5V69erfHjx2vKlCl67733NGbMGI0ZM0abN29ujfQAAEAb0yo3jsvPz9f555+v+fPnS/rkKp5u3brp29/+tu68884G848dO1ZVVVV69dVXY9MuvPBCnXPOOXr88cePubyKigplZWWpU1bHNvaTUJLOYXH5YTZJPwm5nMNS73IOi9MvJCfaOSxt8ychl7xPvHNY2t5PQsk6hyXYFn8SSsI+K2qM/re8TOXl5crMzDzqvC1+hKW2tlbr169XQUHBvxYSCKigoEBr1qxpNGbNmjVx80tSYWFhk/PX1NSooqIi7gEAANqvFi9Y9u/fr/r6emVnZ8dNz87OVmlpaaMxpaWlCc1fVFSkrKys2INLmgEAaN/a5J1uZ82apfLy8thj165dyU4JAAC0oha/rLlz584KBoPau3dv3PS9e/cqJyen0ZicnJyE5g+HwwqHwy2TMAAA8L0WP8KSmpqqIUOGaPny5bFp0WhUy5cv19ChQxuNGTp0aNz8kvT66683OT8AADixtMqN42bMmKFJkybpvPPO0wUXXKCHHnpIVVVVuuGGGyRJEydO1GmnnaaioiJJ0i233KJLL71UP/nJT/TFL35Rzz33nNatW6df/OIXrZEeAABoY1qlYBk7dqw++ugj3XPPPSotLdU555yjZcuWxU6sLSkpibvb6kUXXaRnnnlG//Ef/6G77rpLffr00ZIlSzRgwIDWSA8AALQxrXIfluON+7AkhvuwNB/3YUlksdyHJSHchyUB3Iel+drvfVjaxVhCR0Qi9VYfKKcBvVyKHYcOHXT68rfnMi6PcViwcRlJzLP/ALvsoANesgoW+4Z2GkvI6cu7DQ6sI8cRnxy2U7L+zXRaX4cv4ahLEe7yT2U0SQOnuoydlqSPsIxdWyWSb5u8rBkAAJxYKFgAAIDvUbAAAADfo2ABAAC+R8ECAAB8j4IFAAD4HgULAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfI+CBQAA+F5KshNoSZFInTyLscRtYmKxAfuaL2Ac6sWAy0DgDsPay76toi7r68ALOGxfh74R8OzX1/OSNka8daTn0Dfk8hl0Wa4kubS1A+PweTDGPmeHUEWj9vuOOvvFyjgs12Wf5QXtY+vtU3bt0Q7LtV+ycdvxNAtHWAAAgO9RsAAAAN+jYAEAAL5HwQIAAHyPggUAAPgeBQsAAPA9ChYAAOB7FCwAAMD3KFgAAIDvUbAAAADfo2ABAAC+1+IFS1FRkc4//3xlZGSoa9euGjNmjLZs2XLUmIULF8rzvLhHWlpaS6cGAADaqBYvWFauXKlp06Zp7dq1ev311xWJRHT55ZerqqrqqHGZmZnas2dP7LFz586WTg0AALRRLT5a87Jly+KeL1y4UF27dtX69es1fPjwJuM8z1NOTk5LpwMAANqBFi9YPq28vFySdMoppxx1vsrKSnXv3l3RaFTnnnuufvSjH6l///6NzltTU6OamprY84qKCklSfb2xGqU+6DCEuMuI2i7DvCdrePlAwP6gXCAQdFiwfawXtO/mLsOtu/Qrz6ljuRw4tV9uwHPpG/Zt5dInJce2dth1BB3yrquvt4+N1FnHRhxio9GofazNjv3/GJeN5NA1og7BDhk7xbqwXdtEvo9a9aTbaDSqW2+9VRdffLEGDBjQ5Hx9+/bVE088oZdfflmLFi1SNBrVRRddpN27dzc6f1FRkbKysmKPbt26tdYqAAAAH/CMy7/qx3DzzTfrtdde09tvv63TTz+92XGRSERnn322xo8fr3vvvbfB640dYenWrZtSAmnyLCrxYNC+bvOcjjjYx6Y4/EcadPpv1uEoiedwQI8jLM1n7P/75ghLwsHWTrQjLJE6+1iX9XU5whJ0ObLj0DeSdYTFJdb+CItR2cGPVV5erszMzKPO22o/CU2fPl2vvvqqVq1alVCxIkmhUEif/exntW3btkZfD4fDCofDLZEmAABoA1r8JyFjjKZPn67FixfrzTffVM+ePRN+j/r6er3//vvKzc1t6fQAAEAb1OJHWKZNm6ZnnnlGL7/8sjIyMlRaWipJysrKUnp6uiRp4sSJOu2001RUVCRJ+sEPfqALL7xQZ555psrKyvTAAw9o586duvHGG1s6PQAA0Aa1eMHy2GOPSZJGjBgRN/3JJ5/U5MmTJUklJSVxvzt//PHHuummm1RaWqpOnTppyJAhWr16tT7zmc+0dHoAAKANatWTbo+XiooKZWVlcdJtM3HSbQKxnHTbbJx0mxhOum0+Trr1f+zxOOmWsYQAAIDvUbAAAADfo2ABAAC+R8ECAAB8j4IFAAD4HgULAADwPQoWAADge602llAyBALecb8Pi8s9IJyGAXe62N7l3iL26xtMse9uwWDIOtYLutw7xuH+Lw7t7BAqmahLsHVkikM7h0L2sakOsZIUSnH4DDtsJ5d9h8u9VCorD1nHVlUdto6NRh3upWLst7HTjcZcblPmEJusm6MlY7mJ3AqOIywAAMD3KFgAAIDvUbAAAADfo2ABAAC+R8ECAAB8j4IFAAD4HgULAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfI+CBQAA+F5KshNoScGUFHkW470HAvZjxAdcxpdPFoecjUNsIGhfH4dC9l01JRSyjlXAfrmBFIecU4LWsU7D2puodWzIYX07pNlvo5NPcti+ksIh+35pnNraPramptZ+uVH7bVxdbb9cyX7f4bKbddpG9otV1KGdk8ZlhS2DE9k+HGEBAAC+R8ECAAB8j4IFAAD4XosXLHPmzJHneXGPfv36HTXmxRdfVL9+/ZSWlqaBAwdq6dKlLZ0WAABow1rlCEv//v21Z8+e2OPtt99uct7Vq1dr/PjxmjJlit577z2NGTNGY8aM0ebNm1sjNQAA0Aa1SsGSkpKinJyc2KNz585Nzvvwww9r1KhRuv3223X22Wfr3nvv1bnnnqv58+e3RmoAAKANapWCZevWrcrLy1OvXr00YcIElZSUNDnvmjVrVFBQEDetsLBQa9asaTKmpqZGFRUVcQ8AANB+tXjBkp+fr4ULF2rZsmV67LHHVFxcrGHDhungwYONzl9aWqrs7Oy4adnZ2SotLW1yGUVFRcrKyoo9unXr1qLrAAAA/KXFC5bRo0fruuuu06BBg1RYWKilS5eqrKxML7zwQostY9asWSovL489du3a1WLvDQAA/KfV73TbsWNHnXXWWdq2bVujr+fk5Gjv3r1x0/bu3aucnJwm3zMcDiscDrdongAAwL9a/T4slZWV2r59u3Jzcxt9fejQoVq+fHnctNdff11Dhw5t7dQAAEAb0eIFy8yZM7Vy5Urt2LFDq1ev1tVXX61gMKjx48dLkiZOnKhZs2bF5r/lllu0bNky/eQnP9EHH3ygOXPmaN26dZo+fXpLpwYAANqoFv9JaPfu3Ro/frwOHDigLl266JJLLtHatWvVpUsXSVJJSYkCgX/VSRdddJGeeeYZ/cd//Ifuuusu9enTR0uWLNGAAQNaOjUAANBGecZlKEufqKioUFZWltLTMqxGaw62xdGaHRYbDNqPBJziMCJvaqr9eUepIftYRmtOJJTRmhPRFkdrPvC/9reBcIk9dPiwdWx91L6toknaRtFovXVs0iRptObKwwdVXl6uzMzMo87b6ifdHk/RaNSqYPEcfhnzHIodt1rHPtilRHWJDTi0VSjVpVBKtY71gg5Fh8NyU1Ptv4RtPgMtIeiw2JPC9gVaRrrbbiyUYp+4y5dSXb39h6m+3r6wDDj8wyKHz7DxHApDz76dnYpKh2/wqEPxbxwKNKf1dSoM7dY3kXwZ/BAAAPgeBQsAAPA9ChYAAOB7FCwAAMD3KFgAAIDvUbAAAADfo2ABAAC+R8ECAAB8j4IFAAD4HgULAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA33Mbl91n6uqNPJvRz+1HTJdnPxq3y0jtCjgkHbBqpCPsYx1GLneLdcg5FLSv6VND9h+vtHCqdaxLv3IZmj4g+9jUUNA6NuiwjSQpFHRoMIfYgBe1jq1z6FupqSHr2FCqfb9MidRbx0ZrHXYADv3SJdZz2M8ah+VGo/b9yiXWWO6kE9nncIQFAAD4HgULAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfI+CBQAA+B4FCwAA8D0KFgAA4HstXrD06NFDnuc1eEybNq3R+RcuXNhg3rS0tJZOCwAAtGEtPvjhn//8Z9XX/2uQq82bN+vzn/+8rrvuuiZjMjMztWXLlthzl0GjAABA+9PiBUuXLl3ins+bN0+9e/fWpZde2mSM53nKyclp6VQAAEA70eIFy7+rra3VokWLNGPGjKMeNamsrFT37t0VjUZ17rnn6kc/+pH69+/f5Pw1NTWqqamJPa+oqJAkeQHvuB+dcRkGXLLP1WU93drIZdhz+9hInf3Q9MGgfWw4nGodGwrYt3MoaB8bcFhuIkO9N1iu7IemDzr8OG2M/XI/iQ9ax4YcEvdC9rG1Dp+HUMh+fdPCIevYuro661gTtV9fOfRpOewrTdQ+NiL7tqp3aCuXbwbb7xUjqb6Zm6hVT7pdsmSJysrKNHny5Cbn6du3r5544gm9/PLLWrRokaLRqC666CLt3r27yZiioiJlZWXFHt26dWuF7AEAgF94xuVfqmMoLCxUamqqXnnllWbHRCIRnX322Ro/frzuvffeRudp7AhLt27dlBLqYFXlBQP2dVvQ4b9Zt1j7/5KCLv8VOrRVKDVsHZsatj8RO80h9qQO6daxHdIdck63b6u2eIQl7HC0IZzidlQ1nOJyhMV+2c39r7IxlYdrjj1TE/Z/fNA69oBDbGXVYevYw4errWMjEfujFfLsN5KJ2n8eIg5HoyKRiHWsqXc4WmnZVMYYRUytysvLlZmZedR5W+0noZ07d+qNN97QSy+9lFBcKBTSZz/7WW3btq3JecLhsMJh+x06AABoW1rtJ6Enn3xSXbt21Re/+MWE4urr6/X+++8rNze3lTIDAABtTasULNFoVE8++aQmTZqklJT4gzgTJ07UrFmzYs9/8IMf6A9/+IP++c9/asOGDfrqV7+qnTt36sYbb2yN1AAAQBvUKj8JvfHGGyopKdHXv/71Bq+VlJQo8G/nQXz88ce66aabVFpaqk6dOmnIkCFavXq1PvOZz7RGagAAoA1q1ZNuj5eKigplZWVx0m1zYznpttk46TaB5XLSbUI46bb5OOk2gdh2fNItYwkBAADfo2ABAAC+R8ECAAB8j4IFAAD4HgULAADwPQoWAADgexQsAADA91ptLKFkCHh2o4F7DtfauwzI7TbquUvOLtfaO9x3ot5+2HOX2KixX1+n4dYd7ofici+V1JDDPXpc7g3k2bezw+1MlOJwbyDJ7b5EQYfEg7Y3rpCUmmKfc1qqff/IOCnVOjbk0D8y0+2/qmoj9vuOKof73RiH/U6kzn4bHXLoV3Wyv/9LNGp7IxZJzdxEHGEBAAC+R8ECAAB8j4IFAAD4HgULAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfI+CBQAA+B4FCwAA8D0KFgAA4HsULAAAwPfsx+z2oaAneRajvQdsgv6PQ6jk2Q8DbqL2Q5dHZZ90MGhf43oOw557sl9fOQzzXh+1H5q+3mG5LkPTBzz7oenDIfvYUErIOtatb7h8CKVg0D4+JcUhNmC/ziGX5Tq0dUbIPlZR+/4RCNivb02t/Wdp974y61jPYV8ph22056My69iDlYetY6OW30nGGNU2czfLERYAAOB7FCwAAMD3KFgAAIDvUbAAAADfS7hgWbVqla644grl5eXJ8zwtWbIk7nVjjO655x7l5uYqPT1dBQUF2rp16zHfd8GCBerRo4fS0tKUn5+vd999N9HUAABAO5VwwVJVVaXBgwdrwYIFjb5+//3365FHHtHjjz+ud955Rx06dFBhYaGqq6ubfM/nn39eM2bM0OzZs7VhwwYNHjxYhYWF2rdvX6LpAQCAdsgzxlhfP+V5nhYvXqwxY8ZI+uToSl5enm677TbNnDlTklReXq7s7GwtXLhQ48aNa/R98vPzdf7552v+/PmSPrk8qlu3bvr2t7+tO++8s8H8NTU1qqmpiT2vqKhQt27dlB7uIM/iOuNAwOFSXZfrmh1CU1wuxXa4RDAYtL8SPuhw2WtqOGwdm5aWbh17UoeTrGM7nOwQm26/vulh+3ZOS7XfvqGUZF3y7nZZc8jhsubUJF3WXF9vf7n9wYP2l64ePlRpHatonXUolzU3X1u8rLmq+rDKy8uVmZl51Hlb9ByW4uJilZaWqqCgIDYtKytL+fn5WrNmTaMxtbW1Wr9+fVxMIBBQQUFBkzFFRUXKysqKPbp169aSqwEAAHymRQuW0tJSSVJ2dnbc9Ozs7Nhrn7Z//37V19cnFDNr1iyVl5fHHrt27WqB7AEAgF+1yTvdhsNhhR1+IgAAAG1Lix5hycnJkSTt3bs3bvrevXtjr31a586dFQwGE4oBAAAnlhYtWHr27KmcnBwtX748Nq2iokLvvPOOhg4d2mhMamqqhgwZEhcTjUa1fPnyJmMAAMCJJeGfhCorK7Vt27bY8+LiYm3cuFGnnHKKzjjjDN16662677771KdPH/Xs2VN333238vLyYlcSSdLIkSN19dVXa/r06ZKkGTNmaNKkSTrvvPN0wQUX6KGHHlJVVZVuuOEG9zUEAABtXsIFy7p163TZZZfFns+YMUOSNGnSJC1cuFDf+973VFVVpalTp6qsrEyXXHKJli1bprS0tFjM9u3btX///tjzsWPH6qOPPtI999yj0tJSnXPOOVq2bFmDE3EBAMCJyek+LH5RUVHxySXOHTKs7ovici8Vl9uwuDW8S872saGQ/XnaKQ73cAmFUu1jU+1jww73cElPTzv2TE3FptmfVB52uJdKetihrdLtY53uheT4SXK4fYxOCtl/ljLCyfkMV1XZ32vj4wP7jz1TEyK19ssNOvSPQND+s7S3zP6+MyGHC0OCKUHr2A/3lVvHfnTgY+vYSMTuPjtRY3Sg/ODxvw8LAABAa6BgAQAAvkfBAgAAfI+CBQAA+B4FCwAA8D0KFgAA4HsULAAAwPcoWAAAgO9RsAAAAN+jYAEAAL5HwQIAAHyPggUAAPgeBQsAAPA9ChYAAOB7FCwAAMD3UpKdQEsyxljFeS2cR7OX6yVryQ4s2/j/gh1C6+1jo/ax0UitdWwkYL99Aw7tHK0L2sfW27eV5/DvTzgt1T7Y8WPkObxBStB+O6UE7NvafgtLNV6ddWyktsY6trq62jo2NWT/VRU0DtvX4TN8UjhkHRtMsV/fjhknWcdWHay0jj0UjVrFRaPN/wxxhAUAAPgeBQsAAPA9ChYAAOB7FCwAAMD3KFgAAIDvUbAAAADfo2ABAAC+R8ECAAB8j4IFAAD4HgULAADwvYQLllWrVumKK65QXl6ePM/TkiVLYq9FIhHdcccdGjhwoDp06KC8vDxNnDhRH3744VHfc86cOfI8L+7Rr1+/hFcGAAC0TwkXLFVVVRo8eLAWLFjQ4LVDhw5pw4YNuvvuu7Vhwwa99NJL2rJli6688spjvm///v21Z8+e2OPtt99ONDUAANBOJTzC0ujRozV69OhGX8vKytLrr78eN23+/Pm64IILVFJSojPOOKPpRFJSlJOTk2g6AADgBNDq57CUl5fL8zx17NjxqPNt3bpVeXl56tWrlyZMmKCSkpIm562pqVFFRUXcAwAAtF/2Y1g3Q3V1te644w6NHz9emZmZTc6Xn5+vhQsXqm/fvtqzZ4/mzp2rYcOGafPmzcrIyGgwf1FRkebOndtguqeA1VDxxn6EeEWN3ZDakuQl6ZRnz2F96xMYCrzBcuvt2yrg2ccGg/XWsfX19rG1tbXWsYkMuf5pddGgdaxL5zip3n53kiL7nIMBtw9SWor9OoeDie9vjgjKvm+FHLZTesg+NrNDyGG5adaxKSn227jOvpl1KGDfVkGHfVbY4Zs5xeHzEKmzX25NjV1bRRP4Am61r8xIJKLrr79exhg99thjR5139OjRuu666zRo0CAVFhZq6dKlKisr0wsvvNDo/LNmzVJ5eXnssWvXrtZYBQAA4BOtcoTlSLGyc+dOvfnmm0c9utKYjh076qyzztK2bdsafT0cDiscDrdEqgAAoA1o8SMsR4qVrVu36o033tCpp56a8HtUVlZq+/btys3Nben0AABAG5RwwVJZWamNGzdq48aNkqTi4mJt3LhRJSUlikQi+vKXv6x169bp17/+terr61VaWqrS0tK43/NHjhyp+fPnx57PnDlTK1eu1I4dO7R69WpdffXVCgaDGj9+vPsaAgCANi/hn4TWrVunyy67LPZ8xowZkqRJkyZpzpw5+u1vfytJOuecc+LiVqxYoREjRkiStm/frv3798de2717t8aPH68DBw6oS5cuuuSSS7R27Vp16dIl0fQAAEA7lHDBMmLECJmjnNV7tNeO2LFjR9zz5557LtE0AADACYSxhAAAgO9RsAAAAN+jYAEAAL5HwQIAAHyPggUAAPgeBQsAAPA9ChYAAOB7FCwAAMD3WmXww6Tx/u+RoObc7K4pUftQOYxcLs9+VHsnnk0Dx2JdODSW/SjvkkPfqK+3H6u9PlpvHWsUtI5NceiUtdX2///UpdrHBj233ZiJ2i+7vs5h3+HQp6MB+04dTrFf38wOqdaxtUH7Ph1ODVnHRqL2e56qwxHr2KjD57+u3v4zHInUHnumJtRH7fuV7dZN5FPAERYAAOB7FCwAAMD3KFgAAIDvUbAAAADfo2ABAAC+R8ECAAB8j4IFAAD4HgULAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfM9tXHafiUaNPIuRxI2xH+bdPlKKOgS7VJqeTSMd4dBWMvZDl5t6+5yNZz+sfTRgv9x6+8XKOPSs+qj90PQu2zclaB8bTrVv52AgzTpWkoJeyDo24tA/ahw+hlFjH+w5fA4jEfvY6po669hg0H6PlxJMtY4NOGzfiMOuMlJr31aHamqtY+ui9ts3arnvSOT7lyMsAADA9yhYAACA71GwAAAA30u4YFm1apWuuOIK5eXlyfM8LVmyJO71yZMny/O8uMeoUaOO+b4LFixQjx49lJaWpvz8fL377ruJpgYAANqphAuWqqoqDR48WAsWLGhynlGjRmnPnj2xx7PPPnvU93z++ec1Y8YMzZ49Wxs2bNDgwYNVWFioffv2JZoeAABohxK+Smj06NEaPXr0UecJh8PKyclp9ns++OCDuummm3TDDTdIkh5//HH97ne/0xNPPKE777wz0RQBAEA70yrnsLz11lvq2rWr+vbtq5tvvlkHDhxoct7a2lqtX79eBQUF/0oqEFBBQYHWrFnTaExNTY0qKiriHgAAoP1q8YJl1KhR+tWvfqXly5frxz/+sVauXKnRo0ervombUuzfv1/19fXKzs6Om56dna3S0tJGY4qKipSVlRV7dOvWraVXAwAA+EiL3zhu3Lhxsb8HDhyoQYMGqXfv3nrrrbc0cuTIFlnGrFmzNGPGjNjziooKihYAANqxVr+suVevXurcubO2bdvW6OudO3dWMBjU3r1746bv3bu3yfNgwuGwMjMz4x4AAKD9avWCZffu3Tpw4IByc3MbfT01NVVDhgzR8uXLY9Oi0aiWL1+uoUOHtnZ6AACgDUi4YKmsrNTGjRu1ceNGSVJxcbE2btyokpISVVZW6vbbb9fatWu1Y8cOLV++XFdddZXOPPNMFRYWxt5j5MiRmj9/fuz5jBkz9F//9V966qmn9Pe//10333yzqqqqYlcNAQCAE1vC57CsW7dOl112Wez5kXNJJk2apMcee0ybNm3SU089pbKyMuXl5enyyy/Xvffeq3A4HIvZvn279u/fH3s+duxYffTRR7rnnntUWlqqc845R8uWLWtwIi4AADgxecZlqGKfqKioUFZWljLSM6xGIk7WaM2ew0igAYcRlwOey8inLrEuI/Laj0CcmmI/Gm8wZH9eusNgzU6jNQdT7NsqPWw/sm3myeFjz9SEUzt2sI7tcJLbaM3hkH3/SA/Z9+k0h0seQgH7/uEZ+555sNz+FhKHD1VZx3Y4yb5fBh1Ga9790UHr2Ihnv4GjDvu70o/st9Hu0jLr2EPVEas4Y4wOHjqo8vLyY56P2uJXCSWT533ySJRL0eF0EpDDgp0KJYdh7V2WbKL2C66X/bDndfX2Q7W7tLP9UuWwtlIw6tKxHL4IPfusU1Ls+4btsPZHRML2X+Amal/s1NU5FPCefc5e1KFgOVRrHVtf69A/gvaxJmCf86Fa+0/x4SZu5dEc9Q7fLIdrHLaRQ87G2G2jRA4YMPghAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfI+CBQAA+B4FCwAA8D0KFgAA4HsULAAAwPcoWAAAgO9RsAAAAN+jYAEAAL6XkuwEWlIg4MnzLIZstx+5XHIb2d6a/cD0kk0THRFwWXCSGiuR4cs/rT5qP9y6/cD0ctvAUfv1rau3j62ttU/6YGWNdWy9Q86SdFKa/Zaqrwtbx6YE7dvLM/b9Mlpnv741VbXWsWkp9v8f1zl8VVVXR6xjKw/bt1X5Ifvl1tbZfylVHbbfRpE6+35lLPc7ieyfOcICAAB8j4IFAAD4HgULAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfI+CBQAA+F7CBcuqVat0xRVXKC8vT57nacmSJXGve57X6OOBBx5o8j3nzJnTYP5+/folvDIAAKB9Srhgqaqq0uDBg7VgwYJGX9+zZ0/c44knnpDnebr22muP+r79+/ePi3v77bcTTQ0AALRTCY8oNXr0aI0ePbrJ13NycuKev/zyy7rsssvUq1evoyeSktIgFgAAQGrlc1j27t2r3/3ud5oyZcox5926davy8vLUq1cvTZgwQSUlJU3OW1NTo4qKirgHAABov+zH7G6Gp556ShkZGbrmmmuOOl9+fr4WLlyovn37as+ePZo7d66GDRumzZs3KyMjo8H8RUVFmjt3boPpR85/SVTAsx+e3mVg+0SG1W5ZDutr7Gtcm20TW67sY+sc2tlz2kYusfbr6zltX/th7SN1ddaxtbUO6+vw+ZXc1jkatY9NCQatY11yrovYbyfV2bd1aihkv9io/X4nErEO1eFq+7aqqKy2jq2J1Ccltt6hPzvsspqtVY+wPPHEE5owYYLS0tKOOt/o0aN13XXXadCgQSosLNTSpUtVVlamF154odH5Z82apfLy8thj165drZE+AADwiVY7wvLHP/5RW7Zs0fPPP59wbMeOHXXWWWdp27Ztjb4eDocVDoddUwQAAG1Eqx1h+eUvf6khQ4Zo8ODBCcdWVlZq+/btys3NbYXMAABAW5NwwVJZWamNGzdq48aNkqTi4mJt3Lgx7iTZiooKvfjii7rxxhsbfY+RI0dq/vz5seczZ87UypUrtWPHDq1evVpXX321gsGgxo8fn2h6AACgHUr4J6F169bpsssuiz2fMWOGJGnSpElauHChJOm5556TMabJgmP79u3av39/7Pnu3bs1fvx4HThwQF26dNEll1yitWvXqkuXLommBwAA2iHPJO9SlRZTUVGhrKwsdcrIsroSxTicGe3Sei5N73LFTTBo/0tgMGB/ZUPAIWcFXK4icYgN2LeVcblKKFnbN2i/XJerXtLD9leQhB1iJSnNIT49nGod2zavEqqxDu14kn1bneTQzocP2+e8Y2+Zdez+isPWscm6SsjlSj/brzNjjKqqD6m8vFyZmZlHnZexhAAAgO9RsAAAAN+jYAEAAL5HwQIAAHyPggUAAPgeBQsAAPA9ChYAAOB7FCwAAMD3Wm3ww2RIxi3wXBZpXG5o5jCWt0s7uQw/HnW4CZvnMuq5Q1nuuQy37sJz6ln2oS6r69BWhxxiIw43ypKk2lr7m2VVV9daxwZdbhwXdbljpX1be1GHm845bOPa6oh1bF2Nfc5Vh+xvOne4xj7nWpcbxznc/K2u3n65AcsdbSI3UOUICwAA8D0KFgAA4HsULAAAwPcoWAAAgO9RsAAAAN+jYAEAAL5HwQIAAHyPggUAAPgeBQsAAPA9ChYAAOB7FCwAAMD3KFgAAIDvUbAAAADfo2ABAAC+l5LsBHzB85IR6jLKuxwGl1c0geG8W5IXtV+u5znEOqyuw+Z1DLZnP7i85Dl0aNvh5T+JdRjWPhCxjpUkL+Cyzg77DutIt+CAwwciIPudVkXQvn8EHT7DkYh93zpYY/9pqo3YJ11bZ9/ONQ451zt8KQWDdnEmge8jjrAAAADfo2ABAAC+R8ECAAB8L6GCpaioSOeff74yMjLUtWtXjRkzRlu2bImbp7q6WtOmTdOpp56qk08+Wddee6327t171Pc1xuiee+5Rbm6u0tPTVVBQoK1btya+NgAAoF1KqGBZuXKlpk2bprVr1+r1119XJBLR5Zdfrqqqqtg83/3ud/XKK6/oxRdf1MqVK/Xhhx/qmmuuOer73n///XrkkUf0+OOP65133lGHDh1UWFio6upqu7UCAADtimcSOUX3Uz766CN17dpVK1eu1PDhw1VeXq4uXbromWee0Ze//GVJ0gcffKCzzz5ba9as0YUXXtjgPYwxysvL02233aaZM2dKksrLy5Wdna2FCxdq3Lhxx8yjoqJCWVlZ6nhyluVVDvZncztdreNwlZDLFQYuVza5XVFl/wukS85OsfahSbtKyKVPJu8qIYflOlzlI3GVUEKxDlcJpZ5gVwlVR+zbqrbOYbk1tdaxblcJ2V0mZIxRNBpReXm5MjMzjzqv0zks5eXlkqRTTjlFkrR+/XpFIhEVFBTE5unXr5/OOOMMrVmzptH3KC4uVmlpaVxMVlaW8vPzm4ypqalRRUVF3AMAALRf1gVLNBrVrbfeqosvvlgDBgyQJJWWlio1NVUdO3aMmzc7O1ulpaWNvs+R6dnZ2c2OKSoqUlZWVuzRrVs329UAAABtgHXBMm3aNG3evFnPPfdcS+bTLLNmzVJ5eXnssWvXruOeAwAAOH6sCpbp06fr1Vdf1YoVK3T66afHpufk5Ki2tlZlZWVx8+/du1c5OTmNvteR6Z++kuhoMeFwWJmZmXEPAADQfiVUsBhjNH36dC1evFhvvvmmevbsGff6kCFDFAqFtHz58ti0LVu2qKSkREOHDm30PXv27KmcnJy4mIqKCr3zzjtNxgAAgBNLQgXLtGnTtGjRIj3zzDPKyMhQaWmpSktLdfjwYUmfnCw7ZcoUzZgxQytWrND69et1ww03aOjQoXFXCPXr10+LFy+W9MlVCbfeeqvuu+8+/fa3v9X777+viRMnKi8vT2PGjGm5NQUAAG1WQoMfPvbYY5KkESNGxE1/8sknNXnyZEnST3/6UwUCAV177bWqqalRYWGhHn300bj5t2zZErvCSJK+973vqaqqSlOnTlVZWZkuueQSLVu2TGlpaRarBAAA2hun+7D4BfdhSTCW+7A0P9Y+lPuwJBTLfViOVzD3YWk+7sPSfMfjPiwJHWHxP7te7bKTduHw+VW9Q7UTdfjwey71rcOOMlm8NpizcfkSNS5f3vZ9MmkFuCS5FNIOmSerb3kORYfLNj7stpWtI132lZF6+21U5/APqcs/s1GHz7AXsC8JTMDuc2SMkaKRZs3L4IcAAMD3KFgAAIDvUbAAAADfo2ABAAC+R8ECAAB8j4IFAAD4HgULAADwPQoWAADgexQsAADA9yhYAACA71GwAAAA36NgAQAAvkfBAgAAfK9djNZs/m8EYeMykrDtspMU7LKubW/84WRqe62VtD7pslwH7st1GWXaYWTcJLWY23LtY10H1bYVddlXOsVahzru35P0IbbMOZHv73ZRsBw8eFCSVF51MMmZAADQRrnUOnVuiz548KCysrKOOo9nknFYooVFo1F9+OGHysjIkOc1rOMrKirUrVs37dq1S5mZmUnIsO2grZqPtmo+2ioxtFfz0VbN58e2Msbo4MGDysvLUyBw9LNU2sURlkAgoNNPP/2Y82VmZvpmI/kdbdV8tFXz0VaJob2aj7ZqPr+11bGOrBzBSbcAAMD3KFgAAIDvnRAFSzgc1uzZsxUOh5Odiu/RVs1HWzUfbZUY2qv5aKvma+tt1S5OugUAAO3bCXGEBQAAtG0ULAAAwPcoWAAAgO9RsAAAAN+jYAEAAL7X7guWBQsWqEePHkpLS1N+fr7efffdZKfkS3PmzJHneXGPfv36JTstX1i1apWuuOIK5eXlyfM8LVmyJO51Y4zuuece5ebmKj09XQUFBdq6dWtykk2yY7XV5MmTG/SzUaNGJSfZJCsqKtL555+vjIwMde3aVWPGjNGWLVvi5qmurta0adN06qmn6uSTT9a1116rvXv3Jinj5GlOW40YMaJB3/rmN7+ZpIyT57HHHtOgQYNid7MdOnSoXnvttdjrbblPteuC5fnnn9eMGTM0e/ZsbdiwQYMHD1ZhYaH27duX7NR8qX///tqzZ0/s8fbbbyc7JV+oqqrS4MGDtWDBgkZfv//++/XII4/o8ccf1zvvvKMOHTqosLBQ1dXVxznT5DtWW0nSqFGj4vrZs88+exwz9I+VK1dq2rRpWrt2rV5//XVFIhFdfvnlqqqqis3z3e9+V6+88opefPFFrVy5Uh9++KGuueaaJGadHM1pK0m66aab4vrW/fffn6SMk+f000/XvHnztH79eq1bt06f+9zndNVVV+mvf/2rpDbep0w7dsEFF5hp06bFntfX15u8vDxTVFSUxKz8afbs2Wbw4MHJTsP3JJnFixfHnkejUZOTk2MeeOCB2LSysjITDofNs88+m4QM/ePTbWWMMZMmTTJXXXVVUvLxu3379hlJZuXKlcaYT/pRKBQyL774Ymyev//970aSWbNmTbLS9IVPt5Uxxlx66aXmlltuSV5SPtapUyfz3//9322+T7XbIyy1tbVav369CgoKYtMCgYAKCgq0Zs2aJGbmX1u3blVeXp569eqlCRMmqKSkJNkp+V5xcbFKS0vj+llWVpby8/PpZ01466231LVrV/Xt21c333yzDhw4kOyUfKG8vFySdMopp0iS1q9fr0gkEte3+vXrpzPOOOOE71ufbqsjfv3rX6tz584aMGCAZs2apUOHDiUjPd+or6/Xc889p6qqKg0dOrTN96l2MVpzY/bv36/6+nplZ2fHTc/OztYHH3yQpKz8Kz8/XwsXLlTfvn21Z88ezZ07V8OGDdPmzZuVkZGR7PR8q7S0VJIa7WdHXsO/jBo1Stdcc4169uyp7du366677tLo0aO1Zs0aBYPBZKeXNNFoVLfeeqsuvvhiDRgwQNInfSs1NVUdO3aMm/dE71uNtZUkfeUrX1H37t2Vl5enTZs26Y477tCWLVv00ksvJTHb5Hj//fc1dOhQVVdX6+STT9bixYv1mc98Rhs3bmzTfardFixIzOjRo2N/Dxo0SPn5+erevbteeOEFTZkyJYmZoT0ZN25c7O+BAwdq0KBB6t27t9566y2NHDkyiZkl17Rp07R582bOG2uGptpq6tSpsb8HDhyo3NxcjRw5Utu3b1fv3r2Pd5pJ1bdvX23cuFHl5eX6zW9+o0mTJmnlypXJTstZu/1JqHPnzgoGgw3Oft67d69ycnKSlFXb0bFjR5111lnatm1bslPxtSN9iX5mp1evXurcufMJ3c+mT5+uV199VStWrNDpp58em56Tk6Pa2lqVlZXFzX8i962m2qox+fn5knRC9q3U1FSdeeaZGjJkiIqKijR48GA9/PDDbb5PtduCJTU1VUOGDNHy5ctj06LRqJYvX66hQ4cmMbO2obKyUtu3b1dubm6yU/G1nj17KicnJ66fVVRU6J133qGfNcPu3bt14MCBE7KfGWM0ffp0LV68WG+++aZ69uwZ9/qQIUMUCoXi+taWLVtUUlJywvWtY7VVYzZu3ChJJ2Tf+rRoNKqampq236eSfdZva3ruuedMOBw2CxcuNH/729/M1KlTTceOHU1paWmyU/Od2267zbz11lumuLjY/OlPfzIFBQWmc+fOZt++fclOLekOHjxo3nvvPfPee+8ZSebBBx807733ntm5c6cxxph58+aZjh07mpdfftls2rTJXHXVVaZnz57m8OHDSc78+DtaWx08eNDMnDnTrFmzxhQXF5s33njDnHvuuaZPnz6muro62akfdzfffLPJysoyb731ltmzZ0/scejQodg83/zmN80ZZ5xh3nzzTbNu3TozdOhQM3To0CRmnRzHaqtt27aZH/zgB2bdunWmuLjYvPzyy6ZXr15m+PDhSc78+LvzzjvNypUrTXFxsdm0aZO58847jed55g9/+IMxpm33qXZdsBhjzM9+9jNzxhlnmNTUVHPBBReYtWvXJjslXxo7dqzJzc01qamp5rTTTjNjx44127ZtS3ZavrBixQojqcFj0qRJxphPLm2+++67TXZ2tgmHw2bkyJFmy5YtyU06SY7WVocOHTKXX3656dKliwmFQqZ79+7mpptuOmH/gWisnSSZJ598MjbP4cOHzbe+9S3TqVMnc9JJJ5mrr77a7NmzJ3lJJ8mx2qqkpMQMHz7cnHLKKSYcDpszzzzT3H777aa8vDy5iSfB17/+ddO9e3eTmppqunTpYkaOHBkrVoxp233KM8aY43c8BwAAIHHt9hwWAADQflCwAAAA36NgAQAAvkfBAgAAfI+CBQAA+B4FCwAA8D0KFgAA4HsULAAAwPcoWAAAgO9RsAAAAN+jYAEAAL73/wGsTWlN/Jt31gAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "def low_resolution_haze_detection_v1():\n",
    "    subsample = 12\n",
    "    url0 = \"https://thumbnails-v2.createlab.org/thumbnail?root=https%3A%2F%2Fbreathecam.org%2F%23v%3D4654%2C2127%2C4915%2C2322%2Cpts%26t%3D1015.36%26ps%3D0%26bt%3D20240519140924%26et%3D20240519140924%26startDwell%3D0%26endDwell%3D0%26d%3D2024-05-19%26s%3Dclairton4%26fps%3D9&width=400&height=300&format=png&fps=9&tileFormat=mp4&startDwell=0&endDwell=0&fromScreenshot&minimalUI\"\n",
    "    im0 = common.read_image_from_url(url0, subsample=subsample) / 255.0\n",
    "    print(im0.shape)\n",
    "    plt.title(\"Input: no smoke\")\n",
    "    plt.imshow(im0 * 255.0)\n",
    "    plt.show()\n",
    "    url1 = \"https://thumbnails-v2.createlab.org/thumbnail?root=https%3A%2F%2Fbreathecam.org%2F%23v%3D4654%2C2127%2C4915%2C2322%2Cpts%26t%3D984.02%26ps%3D0%26bt%3D20240519135036%26et%3D20240519135036%26startDwell%3D0%26endDwell%3D0%26d%3D2024-05-19%26s%3Dclairton4%26fps%3D9&width=400&height=300&format=png&fps=9&tileFormat=mp4&startDwell=0&endDwell=0&fromScreenshot&minimalUI\"\n",
    "    im1 = common.read_image_from_url(url1, subsample=subsample) / 255.0\n",
    "    plt.title(\"Input: with smoke\")\n",
    "    plt.imshow(im1 * 255.0)\n",
    "    plt.show()\n",
    "    haze_image = common.solve_haze_detection_v1(im0, im1)\n",
    "    plt.title(\"Output: Smoke\")\n",
    "    plt.imshow(haze_image * 255.0)\n",
    "    plt.show()\n",
    "\n",
    "low_resolution_haze_detection_v1()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "venv",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}