UCB.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": "UCB.ipynb",
      "provenance": [],
      "authorship_tag": "ABX9TyOsE3PpETA+ETSojmPGWuGN",
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    },
    "language_info": {
      "name": "python"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/PRAN20/Upper-Confidence-Bound/blob/main/UCB.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "**Importing Libraries**"
      ],
      "metadata": {
        "id": "UnFcUvaWbReI"
      }
    },
    {
      "cell_type": "code",
      "execution_count": 1,
      "metadata": {
        "id": "AoNWd6-2bKxA"
      },
      "outputs": [],
      "source": [
        "import matplotlib.pyplot as plt\n",
        "import numpy as np\n",
        "import pandas as ps"
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "**Importing DataSet**"
      ],
      "metadata": {
        "id": "qEuR9CTHblS1"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "dataset = ps.read_csv('Ads_CTR_Optimisation.csv')"
      ],
      "metadata": {
        "id": "vx-FAiQZbnqE"
      },
      "execution_count": 2,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "source": [
        "**Performing UCB Operation**"
      ],
      "metadata": {
        "id": "6VGBE4OXzBXl"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "import math\n",
        "N = 10000\n",
        "d = 10\n",
        "ads_selected = []\n",
        "numbers_of_selections = [0] * d\n",
        "sums_of_rewards = [0] * d\n",
        "total_reward = 0\n",
        "for n in range(0, N):\n",
        "  ad = 0\n",
        "  max_upper_bound = 0\n",
        "  for i in range(0, d):\n",
        "    if (numbers_of_selections[i] > 0):\n",
        "      average_reward = sums_of_rewards[i] / numbers_of_selections[i]\n",
        "      delta_i = math.sqrt(3/2 * math.log(n + 1) / numbers_of_selections[i])\n",
        "      upper_bound = average_reward + delta_i\n",
        "    else:\n",
        "      upper_bound = 1e400\n",
        "    if (upper_bound > max_upper_bound):\n",
        "      max_upper_bound = upper_bound\n",
        "      ad = i\n",
        "  ads_selected.append(ad)\n",
        "  numbers_of_selections[ad] = numbers_of_selections[ad] + 1\n",
        "  reward = dataset.values[n, ad]\n",
        "  sums_of_rewards[ad] = sums_of_rewards[ad] + reward\n",
        "  total_reward = total_reward + reward"
      ],
      "metadata": {
        "id": "hlo7FzhcdV3s"
      },
      "execution_count": 3,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "source": [
        "**Visualising the Result**"
      ],
      "metadata": {
        "id": "QjH4-CvDzGZQ"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "plt.hist(ads_selected)\n",
        "plt.title('Histogram of ads selections')\n",
        "plt.xlabel('Ads')\n",
        "plt.ylabel('Number of times each ad was selected')\n",
        "plt.show()"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 295
        },
        "id": "gjn0z4TKdZG0",
        "outputId": "b684866e-877d-4fc3-9fe7-5b0f1be997ae"
      },
      "execution_count": 4,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEWCAYAAACXGLsWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de7hVVb3/8fdHLooKgoI+CiikZD87eYsU02MmhrcUj+fkT0tF40SdzLRjJ7HHsjJLK81LZVmSaKYSeSEzjR+p5SkvoOYFNFBBQBCSq5Ao8v39MceSyW6vvSaw12Xv9Xk9z3r2nGPevmttWN89xxhzDEUEZmZmbdmi3gGYmVnjc7IwM7OKnCzMzKwiJwszM6vIycLMzCpysjAzs4qcLKzdSHpW0mH1jqOeJP2bpLmSXpe032ae6zBJ89ortgLXmy3piCqct+n/XXQGThZWSGtfJJLOkPRQaT0i3hsRD1Q4zyBJIalrlUKtt+8Bn4uIbSPiiXoHU2uSbpD0zXxZkX8X1vicLKxTaYAktBvwbJ1jMGt3ThbWbvJ3H5IOkDRV0gpJr0q6Iu32x/RzWaqqOUjSFpIulDRH0iJJN0raLnfe09O21yR9pcV1viZpoqRfSFoBnJGu/RdJyyQtkPQDSd1z5wtJn5U0U9JKSRdL2l3Sn1O8E/L7t3iPrcYqaUtJrwNdgL9KeqHM8VelaqoVkqZJ+tfcth7pL/OlkqYDH2hx7PmS5qeYn5c0vMw1jpE0Pe03X9IXc9s+KunJ9Nn8WdLebbzPsZJeSJ/7BEnb57Yfko5flt7PGZLGAJ8AvpR+t79J++Z/X1tKulLSK+l1paQt07bDJM2TdF76bBdIOrPI+7IaiAi//Kr4AmYDR7QoOwN4qLV9gL8Ap6XlbYFhaXkQEEDX3HGfBGYB70r73g7clLbtBbwOHAJ0J6vmeSt3na+l9RPI/vjpAbwfGAZ0TdebAZybu14AdwG9gPcCa4Ap6frbAdOBUWU+h7Kx5s69Rxuf46nADim284CFwFZp26XAn4DtgYHAM8C8tG1PYC6wS+5z3L3MNRYA/5qW+wD7p+X9gEXAgWRJbVT6nW3Zyu/vHOBhYACwJfAT4Ja0bTdgJXAK0C29n33TthuAb5b7twN8I513R6Af8Gfg4rTtMGBt2qcbcAywGujT1vvyq0bfAfUOwK+O8Ur/4V8HluVeqymfLP4IfB3o2+I8g/jnZDEF+GxufU+yBNAV+GrpSypt2xp4kw2TxR8rxH4ucEduPYCDc+vTgPNz65cDV5Y5V9lYc+cumyxaOd9SYJ+0/CJwVG7bGNYniz3SF/0RQLcK53wZ+DTQq0X5taUv5lzZ88CHWvn9zQCG5/bbOfc7uSD/ebY4X6Vk8QJwTG7bkcDstHwY8I8W/zYWsf4PjVbfl1+1ebkayjbGCRHRu/QCPtvGvqOBdwPPSXpM0kfb2HcXYE5ufQ7Zl9JOadvc0oaIWA281uL4ufkVSe+WdLekhalq6ltA3xbHvJpb/kcr69tuQqwVSfqipBmSlktaRnYnU4ptg/eav05EzCJLel8DFkm6VdIuZS7z72R/lc+R9KCkg1L5bsB5qepoWbr+wHTdlnYD7sjtNwN4O73PgWRf+puitc8vf/3XImJtbn01638X5d6X1YCThVVFRMyMiFPIqhsuAyZK2obsL++WXiH7cirZlaw64lWyqocBpQ2SepBVe2xwuRbr1wLPAUMiohfwZUCb/m4Kx9qm1D7xJeAksqqV3sDyXGwLyL6I8+d+R0T8MiIOSdcPss/1n0TEYxExkuyzvxOYkDbNBS7JJ/yI2DoibmnlNHOBo1vsu1VEzE/bdi/zNisNY93a5/dKhWMqvS+rAScLqwpJp0rqFxHryKqsANYBi9PPd+V2vwX4gqTBkrYluxO4Lf2FORE4TtIHU6Pz16j8xd8TWAG8Luk9wH+11/uqEGslPckSy2Kgq6SvkrWblEwALpDUR9IA4OzSBkl7Sjo8NQa/QXb3s67lBSR1l/QJSdtFxFtkn0Npv58Cn5F0oDLbSDpWUs9WYv0xcImk3dJ5+0kambbdDBwh6SRJXSXtIGnftO1VNvzdtnQLcGE6X1+yasZftLF/kfdlNeBkYdVyFPBs6iF0FXByRPwjVSNdAvxvquIYBowDbiJr53iJ7MvwbICIeDYt30r2l/frZPXYa9q49heBj5M1wv4UuK0d31fZWAu4D7gX+BtZ9csbbFjt9PVU/hLw+3Sdki3JGsD/TtYoviNZ20FrTgNmpyq4z5D1UCIipgKfAn5A1lYyi6yTQmuuAiYBv5e0kqxR+sB0npfJqoPOA5YATwL7pOOuB/ZKv9s7WznvN4GpwFPA08DjqayIVt+X1YYiPPmRdRzpr/llZFVML9U7HrNm4TsLa3iSjpO0dWrz+B7ZX6Sz6xuVWXNxsrCOYCRZI+grwBCyKi3fEpvVkKuhzMysIt9ZmJlZRfUedK0q+vbtG4MGDap3GGZmHcq0adP+HhH9WtvWKZPFoEGDmDp1ar3DMDPrUCTNKbfN1VBmZlaRk4WZmVXkZGFmZhU5WZiZWUVOFmZmVpGThZmZVeRkYWZmFTlZmJlZRU4WZmZWUad8gtuskkFjf1u3a8++9Ni6XdtsU/nOwszMKnKyMDOzispWQ0m6Big72UVEfL4qEZmZWcNp685iKjAN2ArYH5iZXvsC3asfmpmZNYqydxYRMR5A0n8Bh0TE2rT+Y+BPtQnPzMwaQZE2iz5Ar9z6tqnMzMyaRJFkcSnwhKQbJI0HHge+VeTkknpLmijpOUkzJB0kaXtJkyXNTD/7pH0l6WpJsyQ9JWn/3HlGpf1nShq1KW/UzMw2XcVkERE/Bw4E7gBuBw4qVVEVcBVwb0S8B9gHmAGMBaZExBBgSloHOBoYkl5jgGsBJG0PXJRiOAC4qJRgzMysNiomC0kCjgD2iYi7gO6SDihw3HbAocD1ABHxZkQsA0YCpWQzHjghLY8EbozMw0BvSTsDRwKTI2JJRCwFJgNHbcybNDOzzVOkGupHwEHAKWl9JfDDAscNBhYDP5f0hKSfSdoG2CkiFqR9FgI7peX+wNzc8fNSWbnyDUgaI2mqpKmLFy8uEJ6ZmRVVJFkcGBFnAW8ApL/ui3Sd7UrW5fbaiNgPWMX6KifSuYI2nuXYGBFxXUQMjYih/fr1a49TmplZUiRZvCWpC+lLXVI/YF2B4+YB8yLikbQ+kSx5vJqql0g/F6Xt84GBueMHpLJy5WZmViNFksXVZI3bO0q6BHgI+HalgyJiITBX0p6paDgwHZgElHo0jQLuSsuTgNNTr6hhwPJUXXUfMEJSn9SwPSKVmZlZjVQcdTYibpY0jezLXsAJETGj4PnPBm6W1B14ETiTLEFNkDQamAOclPa9BzgGmAWsTvsSEUskXQw8lvb7RkQsKXh9MzNrBxWThaSbIuI04LlWytoUEU8CQ1vZNLyVfQM4q8x5xgHjKl3PzMyqo0g11HvzK6n94v3VCcfMzBpR2WQh6QJJK4G9Ja1Ir5VkDdJ3lTvOzMw6n7LJIiK+HRE9ge9GRK/06hkRO0TEBTWM0czM6qxINdSj6Wls4J3xnk5o6wAzM+tciiSLiyJieWklDdlxUfVCMjOzRlMkWbS2T8VeVGZm1nkUSRZTJV0haff0uoJsBj0zM2sSRZLF2cCbwG3ArWRjRLX6PISZmXVORZ7gXgWMlbRNWjYzsyZTZD6LD0qaTjZxEZL2kfSjqkdmZmYNo0g11PfJJiB6DSAi/ko2qZGZmTWJIsmCiJjboujtKsRiZmYNqkgX2LmSPgiEpG7AOaQqKTMzaw5F7iw+Q9b7qT/ZpEP74t5QZmZNpUhvqL8Dn6hBLGZm1qDKJgtJ19DG/NgR8fmqRGRmZg2nrTuLqTWLwszMGlrZZBER4/PrkraOiNXVD8nMzBpNkYfyDkoP5T2X1v1QnplZkynSG+pK/FCemVlT80N5ZmZWkR/KMzOzivxQnpmZVeSH8szMrKIivaG+I6mXpG6SpkhaLOnUIieXNFvS05KelDQ1lW0vabKkmelnn1QuSVdLmiXpKUn7584zKu0/U9KoTX2zZma2aYpUQ42IiBXAR4HZwB7A/2zENT4cEftGxNC0PhaYEhFDgClpHeBoYEh6jQGuhSy5ABcBBwIHABeVEoyZmdVGkWRRqqo6FvhVRCzfzGuOBEoP/I0HTsiV3xiZh4HeknYm67Y7OSKWRMRSYDJw1GbGYGZmG6FIsrhb0nPA+4EpkvqRzcNdRAC/lzRN0phUtlNELEjLC4Gd0nJ/IN9Fd14qK1e+AUljJE2VNHXx4sUFwzMzsyKKNHCPlfQdYHlEvC1pNdldQBGHRMR8STsCk1PSyZ87JJUdrHBjRMR1wHUAQ4cObZdzmplZpuhDeUsi4u20vCoiFhY8bn76uQi4g6zN4dVUvUT6uSjtPh8YmDt8QCorV25mZjVSKFlsCknbSOpZWgZGAM8Ak4BSj6ZRwF1peRJweuoVNYzsTmYBcB8wQlKf1LA9IpWZmVmNFHmCe1PtBNwhqXSdX0bEvZIeAyZIGg3MAU5K+98DHAPMAlYDZ0J2VyPpYuCxtN83ImJJFeM2M7MWKiYLSQcDT0bEqvR8xf7AVRExp63jIuJFYJ9Wyl8DhrdSHpR5MjwixgHjKsVqZmbVUaQa6lpgtaR9gPOAF4AbqxqVmZk1lCLJYm36q38k8IOI+CHQs7phmZlZIynSZrFS0gXAqcChkrYAulU3LDMzayRF7iz+L7AGGJ26zA4AvlvVqMzMrKEUeShvIXBFbv1l3GZhZtZUiow6O0zSY5Jel/SmpLclbe74UGZm1oEUqYb6AXAKMBPoAfwn8KNqBmVmZo2l6HAfs4AuEfF2RPwcj/pqZtZUivSGWi2pO/BkGlBwAVUcJsTMzBpPkS/909J+nwNWkQ3q9+/VDMrMzBpLkTuLPYBFaba8r1c5HjMza0BF7ixOB/4q6WFJ35V0nKc1NTNrLkWesxgFIGkX4D+AHwK7FDnWzMw6hyKjzp4K/CvwPuDvZF1p/1TluMzMrIEUuTu4kmyk2R8D90fE7KpGZGZmDadim0VE9AU+CWwFXCLpUUk3VT0yMzNrGEWG++gF7ArsBgwCtgPWVTcsMzNrJEWqoR7KvX4QEfOqG5KZmTWaIr2h9q5FIGZm1rg8bIeZmVXkZGFmZhU5WZiZWUVl2ywkXQNEue0R8fmqRGRmZg2nrTuLqcA0sucr9ieb/GgmsC/QvfqhmZlZoyibLCJifESMB/YGDouIayLiGmA4WcIoRFIXSU9IujutD5b0iKRZkm5Lc2Ugacu0PittH5Q7xwWp/HlJR27aWzUzs01VpM2iD9Art75tKivqHGBGbv0y4PsRsQewFBidykcDS1P599N+SNoLOBl4L9kMfT+S1GUjrm9mZpupSLK4FHhC0g2SxgOPA98qcnJJA4BjgZ+ldQGHAxPTLuOBE9LyyLRO2j487T8SuDUi1kTES8As4IAi1zczs/ZR5KG8n0v6HXBgKjo/IhYWPP+VwJeAnml9B2BZRKxN6/OA/mm5PzA3XXOtpOVp//7Aw7lz5o95h6QxwBiAXXfdtWB4ZmZWRNGus2vI5t5eCrxb0qGVDpD0UbIZ9qZtRnyFRcR1ETE0Iob269evFpc0M2saReaz+E+ydocBwJPAMOAvZNVJbTkYOF7SMWQ9qnoBVwG9JXVNdxcDgPlp//lk83vPk9SVbMDC13LlJfljzMysBorcWZwDfACYExEfBvYDllU6KCIuiIgBETGIrIH6DxHxCeB+shn3AEYBd6XlSWmdtP0PERGp/OTUW2owMAR4tMibMzOz9lFk1Nk3IuINSUjaMiKek7TnZlzzfOBWSd8EngCuT+XXAzdJmgUsIUswRMSzkiYA04G1wFkR8fZmXN/MzDZSkWQxT1Jv4E5gsqSlwJyNuUhEPAA8kJZfpJXeTBHxBvCxMsdfAlyyMdc0M7P2U6Q31L+lxa9Jup+sLeHeqkZlZmYNpcidxTsi4sFqBWJmZo3Lo86amVlFThZmZlaRk4WZmVVUMVlIOlHSTEnLJa2QtFLSiloEZ2ZmjaFIA/d3gOMiYkbFPc3MrFMqUg31qhOFmVlza2ta1RPT4lRJt5E9lLemtD0ibq9ybGZm1iDaqoY6Lre8GhiRWw/AycLMrEmUTRYRcWYtAzEzs8ZVpDfU+DQ2VGm9j6Rx1Q3LzMwaSZEG7r0j4p0hySNiKdkw5WZm1iSKJIstJPUprUjano0cU8rMzDq2Il/6lwN/kfQrQGQTE3m4cDOzJlJkiPIbJU0DPpyKToyI6dUNy8zMGkmh6qQ0W91isrm0kbRrRLxc1cjMzKxhFOkNdbykmcBLwIPAbOB3VY7LzMwaSJEG7ouBYcDfImIwMBx4uKpRmZlZQymSLN6KiNfIekVtERH3A0OrHJeZmTWQIm0WyyRtC/wJuFnSImBVdcMyM7NGUuTOYiTZ2FDnAvcCL7DhuFFmZtbJFek6u0rSbsCQiBgvaWugS/VDMzOzRlGkN9SngInAT1JRf7Lhyisdt5WkRyX9VdKzkr6eygdLekTSLEm3SeqeyrdM67PS9kG5c12Qyp+XdOTGv00zM9scRaqhzgIOBlYARMRMYMcCx60BDo+IfYB9gaMkDQMuA74fEXsAS4HRaf/RwNJU/v20H5L2Ak4G3gscBfxIku9szMxqqEiyWBMRb5ZWJHUlm8+iTZF5Pa12S68ADie7UwEYD5yQlkemddL24ZKUym+NiDUR8RIwCzigQNxmZtZOiiSLByV9Gegh6SPAr4DfFDm5pC6SngQWAZPJGseXRcTatMs8smot0s+5AGn7cmCHfHkrx5iZWQ0USRZjgcXA08CngXuAC4ucPCLejoh9gQFkdwPv2cQ4K5I0RtJUSVMXL15crcuYmTWlIr2h1gE/Ta9NEhHLJN0PHAT0ltQ13T0MAOan3eYDA4F5qaprO+C1XHlJ/pj8Na4DrgMYOnRoxWoyMzMrrsidxSaR1K80w56kHsBHgBnA/WTDnAOMAu5Ky5PSOmn7HyIiUvnJqbfUYGAI8Gi14jYzs39WzUmMdgbGp55LWwATIuJuSdOBWyV9E3gCuD7tfz1wk6RZwBKyHlClEW8nANOBtcBZEfF2FeM2M7MWNipZSNoC2DYiVlTaNyKeopXpVyPiRVrpzRQRbwAfK3OuS/CES2ZmdVPkobxfSuolaRvgGWC6pP+pfmhmZtYoirRZ7JXuJE4gm8diMHBaVaMyM7OGUiRZdJPUjSxZTIqItyjwUJ6ZmXUeRZLFT8hmx9sG+GMaVLBim4WZmXUeRZ6zuBq4Olc0R9KHqxeSmZk1miIN3DtJul7S79L6Xqx/HsLMzJpAkWqoG4D7gF3S+t/IJkIyM7MmUSRZ9I2ICcA6eGeQPz8UZ2bWRIoki1WSdiD1gEpzUiyvalRmZtZQijzB/d9k4zPtLul/gX6sH9vJzMyaQJHeUI9L+hCwJyDg+fSshZmZNYmKySINBHgMMCjtP0ISEXFFlWMzM7MGUaQa6jfAG2STH62rbjhmZtaIiiSLARGxd9UjMTOzhlWkN9TvJI2oeiRmZtawitxZPAzckeayeIuskTsioldVIzMzs4ZRJFlcQTZ39tNpmlMzM2syRaqh5gLPOFGYmTWvIncWLwIPpIEE15QK3XXWzKx5FEkWL6VX9/QyM7MmU+QJ7q/XIhAzM2tcZZOFpCsj4lxJv6GVaVQj4viqRmZmZg2jrTuLm9LP79UiEDMza1xlk0VETEuL+0bEVfltks4BHqxmYGZm1jiKdJ1tbQrVMyodJGmgpPslTZf0bEowSNpe0mRJM9PPPqlckq6WNEvSU5L2z51rVNp/piRP6WpmVmNttVmcAnwcGCxpUm5TT2BJgXOvBc5LQ5z3BKZJmkyWaKZExKWSxgJjgfOBo4Eh6XUgcC1woKTtgYuAoWRtJ9MkTYqIpRv3Vs3MbFO11WbxZ2AB0Be4PFe+Eniq0okjYkE6nohYKWkG0B8YCRyWdhsPPECWLEYCN6aH/x6W1FvSzmnfyRGxBCAlnKOAWwq9QzMz22xttVnMAeaQDfWxWSQNAvYDHgF2SokEYCGwU1ruT/a0eMm8VFauvOU1xgBjAHbdddfNDdnMzHKKtFlsFknbAr8Gzo2IFflt6S6iXYYRiYjrImJoRAzt169fe5zSzMySqiYLSd3IEsXNEXF7Kn41VS+Rfi5K5fOBgbnDB6SycuVmZlYjZZOFpCnp52WbcmJJAq4HZrQYR2oS63tYjQLuypWfnnpFDQOWp+qq+8imcu2Tek6NSGVmZlYjbTVw7yzpg8Dxkm4lm8fiHRHxeIVzHwycBjwt6clU9mXgUmCCpNFkbSInpW33kM31PQtYDZyZrrNE0sXAY2m/b5Qau83MrDbaShZfBb5CVu3TcoTZAA5v68QR8RAtEkzO8Fb2D+CsMucaB4xr63pmZlY9bfWGmghMlPSViLi4hjGZmVmDKTLq7MWSjgcOTUUPRMTd1Q3LzMwaScXeUJK+DZwDTE+vcyR9q9qBmZlZ4ygy+dGxZIMJrgOQNB54gqyx2szMmkDR5yx655a3q0YgZmbWuIrcWXwbeELS/WS9mw4lG/zPzMyaRJEG7lskPQB8IBWdHxELqxqVmZk1lCJ3FqURZCdV3NHMzDqlqg8kaGZmHZ+ThZmZVdRmspDURdJztQrGzMwaU5vJIiLeBp6X5NmEzMyaWJEG7j7As5IeBVaVCiPi+KpFZWZmDaVIsvhK1aMwM7OGVuQ5iwcl7QYMiYj/J2lroEv1QzMzs0ZRZCDBTwETgZ+kov7AndUMyszMGkuRrrNnkc16twIgImYCO1YzKDMzayxFksWaiHiztCKpK9lMeWZm1iSKJIsHJX0Z6CHpI8CvgN9UNywzM2skRZLFWGAx8DTwaeAe4MJqBmVmZo2lSG+odWnCo0fIqp+ejwhXQ5mZNZGKyULSscCPgRfI5rMYLOnTEfG7agdnZp3DoLG/rct1Z196bF2u2xkVeSjvcuDDETELQNLuwG8BJwszsyZRpM1iZSlRJC8CK6sUj5mZNaCyyULSiZJOBKZKukfSGZJGkfWEeqzSiSWNk7RI0jO5su0lTZY0M/3sk8ol6WpJsyQ9JWn/3DGj0v4z0/XNzKzG2rqzOC69tgJeBT4EHEbWM6pHgXPfABzVomwsMCUihgBTWD+X99HAkPQaA1wLWXIBLgIOBA4ALiolGDMzq52ybRYRcebmnDgi/ihpUIvikWQJB2A88ABwfiq/MfWyelhSb0k7p30nR8QSAEmTyRLQLZsTm5mZbZwivaEGA2cDg/L7b+IQ5Tul+bwBFgI7peX+wNzcfvNSWbny1uIcQ3ZXwq67evoNM7P2VKQ31J3A9WRtFeva68IREZLa7XmNiLgOuA5g6NChfg7EzKwdFUkWb0TE1e10vVcl7RwRC1I106JUPh8YmNtvQCqbz/pqq1L5A+0Ui5mZFVSk6+xVki6SdJCk/UuvTbzeJKDUo2kUcFeu/PTUK2oYsDxVV90HjJDUJzVsj0hlZmZWQ0XuLN4HnAYczvpqqEjrZUm6heyuoK+keWS9mi4FJkgaDcwBTkq73wMcA8wCVgNnAkTEEkkXs76r7jdKjd1mZlY7RZLFx4B35YcpLyIiTimzaXgr+wbZvBmtnWccMG5jrm1mZu2rSDXUM0DvagdiZmaNq8idRW/gOUmPAWtKhZvYddbMzDqgIsnioqpHYWZmDa3IfBYP1iIQMzNrXEWe4F7J+jm3uwPdgFUR0auagdWTx943M9tQkTuLnqVlSSIbx2lYNYMyM7PGUqQ31DsicydwZJXiMTOzBlSkGurE3OoWwFDgjapFZGZmDadIb6jjcstrgdlkVVFmZtYkirRZbNa8FmZm1vGVTRaSvtrGcRERF1chHjMza0Bt3VmsaqVsG2A0sAPgZGFm1iTamlb18tKypJ7AOWSjwd4KXF7uODMz63zabLOQtD3w38AnyObM3j8iltYiMDMzaxxttVl8FziRbKrS90XE6zWLqkn5yXEza1RtPZR3HrALcCHwiqQV6bVS0orahGdmZo2grTaLjXq622xT1Otuysw2TpGH8sysE3Bits3hZGFmnZbbAduPq5rMzKwiJwszM6vI1VDmumyzdlbP/1PVqgJzsjCrMSdn64hcDWVmZhV1mGQh6ShJz0uaJWlsveMxM2smHSJZSOoC/BA4GtgLOEXSXvWNysyseXSIZAEcAMyKiBcj4k2ykW89W5+ZWY10lAbu/sDc3Po84MD8DpLGAGPS6uuSnt+M6/UF/r4Zx3cm/iw25M9jPX8WG2qIz0OXbdbhu5Xb0FGSRUURcR3ZCLmbTdLUiBjaHufq6PxZbMifx3r+LDbU2T+PjlINNR8YmFsfkMrMzKwGOkqyeAwYImmwpO7AycCkOsdkZtY0OkQ1VESslfQ54D6gCzAuIp6t4iXbpTqrk/BnsSF/Huv5s9hQp/48FBH1jsHMzBpcR6mGMjOzOnKyMDOzipwscjykyHqSBkq6X9J0Sc9KOqfeMdWbpC6SnpB0d71jqTdJvSVNlPScpBmSDqp3TPUk6Qvp/8kzkm6RtFW9Y2pvThaJhxT5J2uB8yJiL2AYcFaTfx4A5wAz6h1Eg7gKuDci3gPsQxN/LpL6A58HhkbEv5B1wjm5vlG1PyeL9TykSE5ELIiIx9PySrIvg/71jap+JA0AjgV+Vu9Y6k3SdsChwPUAEfFmRCyrb1R11xXoIakrsDXwSp3jaXdOFuu1NqRI03455kkaBOwHPFLfSOrqSuBLwLp6B9IABgOLgZ+narmfSdqm3kHVS0TMB74HvAwsAJZHxO/rG1X7c7KwNknaFvg1cG5ErKh3PPUg6aPAooiYVu9YGkRXYH/g2ojYD1gFNG0bn6Q+ZLUQg4FdgG0knVrfqNqfk8V6HlKkBUndyBLFzRFxe73jqaODgeMlzSarnjxc0i/qG1JdzQPmRUTpTnMiWfJoVkcAL0XE4oh4C7gd+GCdY2p3ThbreUiRHEkiq5OeERFX1DueeoqICyJiQEQMIvt38YeI6HR/ORYVEQuBuZL2TEXDgel1DKneXgaGSdo6/buKPD8AAAGxSURBVL8ZTids8O8Qw33UQh2GFGl0BwOnAU9LejKVfTki7qljTNY4zgZuTn9YvQicWed46iYiHpE0EXicrBfhE3TCoT883IeZmVXkaigzM6vIycLMzCpysjAzs4qcLMzMrCInCzMzq8jJwqwKJJ0gKSS9p8z2ByQNrXVcZpvKycKsOk4BHko/zTo8JwuzdpbG0zoEGE0aqlpSD0m3prkf7gB6pPIukm5I8yA8LekL9YvcrDw/wW3W/kaSzfXwN0mvSXo/8CFgdUT8H0l7kz3tC7Av0D/Ng4Ck3vUJ2axtvrMwa3+nkA04SPp5Ctn8D78AiIingKfS9heBd0m6RtJRQFOO7GuNz3cWZu1I0vbA4cD7JAXZOGNBNl7QP4mIpZL2AY4EPgOcBHyyRuGaFeY7C7P29R/ATRGxW0QMioiBwEvANODjAJL+Bdg7LfcFtoiIXwMX0txDfVsD852FWfs6BbisRdmvyWYa7CFpBtnw1aWJlPqTzThX+sPtgppEabaRPOqsmZlV5GooMzOryMnCzMwqcrIwM7OKnCzMzKwiJwszM6vIycLMzCpysjAzs4r+P0F3+9tvTHiiAAAAAElFTkSuQmCC\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ]
    }
  ]
}