Maximizing%2Bthe%2Bprofit%2Bof%2Ban%2Boil%2Bcompany%2BLocal.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "# Maximizing the profit of an oil company\n",
    "\n",
    "This tutorial includes everything you need to set up the decision optimization engines and build mathematical programming models.\n",
    "\n",
    "\n",
    "When you finish this tutorial, you'll have a foundational knowledge of _Prescriptive Analytics_.\n",
    "\n",
    ">This notebook is part of [Prescriptive Analytics for Python](https://rawgit.com/IBMDecisionOptimization/docplex-doc/master/docs/index.html).\n",
    "\n",
    ">It requires a valid subscription to **Decision Optimization on Cloud** or a **local installation of CPLEX Optimizers**. \n",
    "\n",
    "Discover us [here](https://developer.ibm.com/docloud).\n",
    "\n",
    "\n",
    "Table of contents:\n",
    "\n",
    "-  [Describe the business problem](#Describe-the-business-problem)\n",
    "*  [How decision optimization (prescriptive analytics) can help](#How--decision-optimization-can-help)\n",
    "*  [Use decision optimization](#Use-decision-optimization)\n",
    "    *  [Step 1: Download the library](#Step-1:-Download-the-library)\n",
    "    *  [Step 2: Set up the engines](#Step-2:-Set-up-the-prescriptive-engine)\n",
    "    -  [Step 3: Model the data](#Step-3:-Model-the-data)\n",
    "    *  [Step 4: Prepare the data](#Step-4:-Prepare-the-data)\n",
    "    -  [Step 5: Set up the prescriptive model](#Step-5:-Set-up-the-prescriptive-model)\n",
    "        * [Define the decision variables](#Define-the-decision-variables)\n",
    "        * [Express the business constraints](#Express-the-business-constraints)\n",
    "        * [Express the objective](#Express-the-objective)\n",
    "        * [Solve with the Decision Optimization solve service](#Solve-with-the-Decision-Optimization-solve-service)\n",
    "    *  [Step 6: Investigate the solution and run an example analysis](#Step-6:-Investigate-the-solution-and-then-run-an-example-analysis)\n",
    "*  [Summary](#Summary)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "## Describe the business problem\n",
    "\n",
    "* An oil company manufactures different types of gasoline and diesel. Each type of gasoline is produced by blending different types of crude oils that must be purchased. The company must decide how much crude oil to buy in order to maximize its profit while respecting processing capacities and quality levels as well as satisfying customer demand.\n",
    "\n",
    "* Blending problems are a typical industry application of Linear Programming (LP). LP represents real life problems mathematically using an objective function to represent the goal that is to be minimized or maximized, together with a set of linear constraints which define the conditions to be satisfied and the limitations of the real life problem. The function and constraints are expressed in terms of decision variables and the solution, obtained from optimization engines such as IBM® ILOG® CPLEX®, provides the best values for these variables so that the objective function is optimized.\n",
    "\n",
    "* The oil-blending problem consists of calculating different blends of gasoline according to specific quality criteria.\n",
    "\n",
    "* Three types of gasoline are manufactured: super, regular, and diesel. \n",
    "\n",
    "* Each type of gasoline is produced by blending three types of crude oil: crude1, crude2, and crude3. \n",
    "\n",
    "* The gasoline must satisfy some quality criteria with respect to their lead content and their octane ratings, thus constraining the possible blendings.\n",
    "\n",
    "* The company must also satisfy its customer demand, which is 3,000 barrels a day of super, 2,000 of regular, and 1,000 of diesel. \n",
    "\n",
    "* The company can purchase 5,000 barrels of each type of crude oil per day and can process at most 14,000 barrels a day.\n",
    "\n",
    "* In addition, the company has the option of advertising a gasoline, in which case the demand for this type of gasoline increases by ten barrels for every dollar spent. \n",
    "\n",
    "* Finally, it costs four dollars to transform a barrel of oil into a barrel of gasoline."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "## How  decision optimization can help\n",
    "* Prescriptive analytics technology recommends actions based on desired outcomes, taking into account specific scenarios, resources, and knowledge of past and current events. This insight can help your organization make better decisions and have greater control of business outcomes.  \n",
    "\n",
    "* Prescriptive analytics is the next step on the path to insight-based actions. It creates value through synergy with predictive analytics, which analyzes data to predict future outcomes.  \n",
    "\n",
    "* Prescriptive analytics takes that insight to the next level by suggesting the optimal way to handle that future situation. Organizations that can act fast in dynamic conditions and make superior decisions in uncertain environments gain a strong competitive advantage.  \n",
    "<br/>\n",
    "\n",
    "+ For example:\n",
    "    + Automate complex decisions and trade-offs to better manage limited resources.\n",
    "    + Take advantage of a future opportunity or mitigate a future risk.\n",
    "    + Proactively update recommendations based on changing events.\n",
    "    + Meet operational goals, increase customer loyalty, prevent threats and fraud, and optimize business processes.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "## Use decision optimization"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "### Step 1: Download the library\n",
    "\n",
    "Run the following code to install the Decision Optimization CPLEX Modeling library.  The *DOcplex* library contains the two modeling packages, Mathematical Programming (docplex.mp) and Constraint Programming (docplex.cp)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "import sys\n",
    "try:\n",
    "    import docplex.mp\n",
    "except:\n",
    "    if hasattr(sys, 'real_prefix'):\n",
    "        #we are in a virtual env.\n",
    "        !pip install docplex \n",
    "    else:\n",
    "        !pip install --user docplex"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "### Step 2: Set up the prescriptive engine\n",
    "\n",
    "* Subscribe to the [Decision Optimization on Cloud solve service](https://developer.ibm.com/docloud).\n",
    "* Get the service URL and your personal API key and enter your credentials here:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "url = None\n",
    "key = None"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "### Step 3: Model the data\n",
    "\n",
    "* For each type of crude oil, there are capacities of what can be bought, the buying price, the octane level, and the lead level. \n",
    "* For each type of gasoline or diesel, there is customer demand, selling prices, and octane and lead levels.\n",
    "* There is a maximum level of production imposed by the factory's limit as well as a fixed production cost. \n",
    "* There are inventory costs for each type of final product and blending proportions. All of these have actual values in the model.\n",
    "\n",
    "* The maginal production cost and maximum production are assumed to be identical for all oil types.\n",
    "\n",
    "\n",
    "Input data comes as *NumPy* arrays with two dimensions. [NumPy](http://www.numpy.org/) is the fundamental package for scientific computing with Python."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "The first dimension of the *NumPy* array is the number of gasoline types; \n",
    "and for each gasoline type, we have a *NumPy* array containing capacity, price, octane and lead level, in that order."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of gasoline types = 3\n",
      "Number of crude types = 3\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "\n",
    "gas_names = [\"super\", \"regular\", \"diesel\"]\n",
    "\n",
    "gas_data = np.array([[3000, 70, 10, 1], [2000, 60, 8, 2], [1000, 50, 6, 1]])\n",
    "\n",
    "oil_names = [\"crude1\", \"crude2\", \"crude3\"]\n",
    "\n",
    "oil_data = np.array([[5000, 45, 12, 0.5], [5000, 35, 6, 2], [5000, 25, 8, 3]])\n",
    "\n",
    "nb_gas  = len(gas_names)\n",
    "nb_oils = len(oil_names)\n",
    "range_gas = range(nb_gas)\n",
    "range_oil = range(nb_oils)\n",
    "print(\"Number of gasoline types = {0}\".format(nb_gas))\n",
    "print(\"Number of crude types = {0}\".format(nb_oils))\n",
    "\n",
    "# global data\n",
    "production_cost = 4\n",
    "production_max = 14000\n",
    "# each $1 spent on advertising increases demand by 10.\n",
    "advert_return = 10"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "### Step 4: Prepare the data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "[Pandas](http://pandas.pydata.org/) is another Python library that we use to store data. *pandas* contains data structures and data analysis tools for the Python programming language."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "gaspd = pd.DataFrame([(gas_names[i],int(gas_data[i][0]),int(gas_data[i][1]),int(gas_data[i][2]),int(gas_data[i][3])) \n",
    "                      for i in range_gas])\n",
    "oilpd = pd.DataFrame([(oil_names[i],int(oil_data[i][0]),int(oil_data[i][1]),int(oil_data[i][2]),oil_data[i][3]) \n",
    "                      for i in range_oil])\n",
    "gaspd.columns = ['name','demand','price','octane','lead']\n",
    "oilpd.columns= ['name','capacity','price','octane','lead']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Use basic HTML and a stylesheet to format the data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       "body {\n",
       "    margin: 0;\n",
       "    font-family: Helvetica;\n",
       "}\n",
       "table.dataframe {\n",
       "    border-collapse: collapse;\n",
       "    border: none;\n",
       "}\n",
       "table.dataframe tr {\n",
       "    border: none;\n",
       "}\n",
       "table.dataframe td, table.dataframe th {\n",
       "    margin: 0;\n",
       "    border: 1px solid white;\n",
       "    padding-left: 0.25em;\n",
       "    padding-right: 0.25em;\n",
       "}\n",
       "table.dataframe th:not(:empty) {\n",
       "    background-color: #fec;\n",
       "    text-align: left;\n",
       "    font-weight: normal;\n",
       "}\n",
       "table.dataframe tr:nth-child(2) th:empty {\n",
       "    border-left: none;\n",
       "    border-right: 1px dashed #888;\n",
       "}\n",
       "table.dataframe td {\n",
       "    border: 2px solid #ccf;\n",
       "    background-color: #f4f4ff;\n",
       "}\n",
       "    table.dataframe thead th:first-child {\n",
       "        display: none;\n",
       "    }\n",
       "    table.dataframe tbody th {\n",
       "        display: none;\n",
       "    }\n",
       "</style>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "CSS = \"\"\"\n",
    "body {\n",
    "    margin: 0;\n",
    "    font-family: Helvetica;\n",
    "}\n",
    "table.dataframe {\n",
    "    border-collapse: collapse;\n",
    "    border: none;\n",
    "}\n",
    "table.dataframe tr {\n",
    "    border: none;\n",
    "}\n",
    "table.dataframe td, table.dataframe th {\n",
    "    margin: 0;\n",
    "    border: 1px solid white;\n",
    "    padding-left: 0.25em;\n",
    "    padding-right: 0.25em;\n",
    "}\n",
    "table.dataframe th:not(:empty) {\n",
    "    background-color: #fec;\n",
    "    text-align: left;\n",
    "    font-weight: normal;\n",
    "}\n",
    "table.dataframe tr:nth-child(2) th:empty {\n",
    "    border-left: none;\n",
    "    border-right: 1px dashed #888;\n",
    "}\n",
    "table.dataframe td {\n",
    "    border: 2px solid #ccf;\n",
    "    background-color: #f4f4ff;\n",
    "}\n",
    "    table.dataframe thead th:first-child {\n",
    "        display: none;\n",
    "    }\n",
    "    table.dataframe tbody th {\n",
    "        display: none;\n",
    "    }\n",
    "\"\"\"\n",
    "\n",
    "from IPython.core.display import HTML\n",
    "HTML('<style>{}</style>'.format(CSS))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Let's display the data we just prepared."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Gas data:\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>name</th>\n",
       "      <th>demand</th>\n",
       "      <th>price</th>\n",
       "      <th>octane</th>\n",
       "      <th>lead</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>super</td>\n",
       "      <td>3000</td>\n",
       "      <td>70</td>\n",
       "      <td>10</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>regular</td>\n",
       "      <td>2000</td>\n",
       "      <td>60</td>\n",
       "      <td>8</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>diesel</td>\n",
       "      <td>1000</td>\n",
       "      <td>50</td>\n",
       "      <td>6</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      name  demand  price  octane  lead\n",
       "0    super    3000     70      10     1\n",
       "1  regular    2000     60       8     2\n",
       "2   diesel    1000     50       6     1"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Oil data:\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>name</th>\n",
       "      <th>capacity</th>\n",
       "      <th>price</th>\n",
       "      <th>octane</th>\n",
       "      <th>lead</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>crude1</td>\n",
       "      <td>5000</td>\n",
       "      <td>45</td>\n",
       "      <td>12</td>\n",
       "      <td>0.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>crude2</td>\n",
       "      <td>5000</td>\n",
       "      <td>35</td>\n",
       "      <td>6</td>\n",
       "      <td>2.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>crude3</td>\n",
       "      <td>5000</td>\n",
       "      <td>25</td>\n",
       "      <td>8</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     name  capacity  price  octane  lead\n",
       "0  crude1      5000     45      12   0.5\n",
       "1  crude2      5000     35       6   2.0\n",
       "2  crude3      5000     25       8   3.0"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from IPython.display import display\n",
    "print(\"Gas data:\")\n",
    "display(gaspd)\n",
    "print(\"Oil data:\")\n",
    "display(oilpd)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "### Step 5: Set up the prescriptive model\n",
    "\n",
    "#### Create the DOcplex model\n",
    "A model is needed to store all the variables and constraints needed to formulate the business problem and submit the problem to the solve service."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "from docplex.mp.model import Model\n",
    "\n",
    "mdl = Model(name=\"oil_blending\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "#### Define the decision variables\n",
    "\n",
    "For each combination of oil and gas, we have to decide the quantity of oil to use to produce a gasoline. A decision variable will be needed to represent that amount. \n",
    "A matrix of continuous variables, indexed by the set of oils and the set of gasolines needs to be created.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "blends = mdl.continuous_var_matrix(keys1=nb_oils, keys2=nb_gas, lb=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "We also have to decide how much should be spent in advertising for each time of gasoline. To do so, we will create a list of continuous variables, indexed by the gasolines."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "adverts = mdl.continuous_var_list(nb_gas, lb=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "#### Express the business constraints\n",
    "\n",
    "The business constraints are the following:\n",
    "\n",
    "* The demand for each gasoline type must be satisfied. The total demand includes the initial demand as stored in the data,plus a variable demand caused by the advertising. This increase is assumed to be proportional to the advertising cost.\n",
    "- The capacity constraint on each oil type must also be satisfied.\n",
    "- For each gasoline type, the octane level must be above a minimum level, and the lead level must be below a maximum level."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "##### Demand\n",
    "\n",
    "+ For each gasoline type, the total quantity produced must equal the raw demand plus the demand increase created by the advertising."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Model: oil_blending\n",
      " - number of variables: 12\n",
      "   - binary=0, integer=0, continuous=12\n",
      " - number of constraints: 3\n",
      "   - linear=3\n",
      " - parameters: defaults\n"
     ]
    }
   ],
   "source": [
    "# gasoline demand is numpy array field #0\n",
    "mdl.add_constraints(mdl.sum(blends[o, g] for o in range(nb_oils)) == gas_data[g][0] + advert_return * adverts[g]\n",
    "                   for g in range(nb_gas))\n",
    "mdl.print_information()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "##### Maximum capacity\n",
    "\n",
    "+ For each type of oil, the total quantity used in all types of gasolines must not exceed the maximum capacity for this oil."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Model: oil_blending\n",
      " - number of variables: 12\n",
      "   - binary=0, integer=0, continuous=12\n",
      " - number of constraints: 6\n",
      "   - linear=6\n",
      " - parameters: defaults\n"
     ]
    }
   ],
   "source": [
    "mdl.add_constraints(mdl.sum(blends[o,g] for g in range_gas) <= oil_data[o][0]\n",
    "                    for o in range_oil)\n",
    "mdl.print_information()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "##### Octane and Lead levels\n",
    "+ For each gasoline type, the octane level must be above a minimum level, and the lead level must be below a maximum level."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Model: oil_blending\n",
      " - number of variables: 12\n",
      "   - binary=0, integer=0, continuous=12\n",
      " - number of constraints: 12\n",
      "   - linear=12\n",
      " - parameters: defaults\n"
     ]
    }
   ],
   "source": [
    "# minimum octane level\n",
    "# octane is numpy array field #2\n",
    "mdl.add_constraints(mdl.sum(blends[o,g]*(oil_data[o][2] - gas_data[g][2]) for o in range_oil) >= 0\n",
    "                   for g in range_gas)\n",
    "    \n",
    "# maximum lead  level\n",
    "# lead level is numpy array field #3 \n",
    "mdl.add_constraints(mdl.sum(blends[o,g]*(oil_data[o][3] - gas_data[g][3]) for o in range_oil) <= 0\n",
    "                   for g in range_gas)\n",
    "\n",
    "mdl.print_information()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "##### Maximum total production\n",
    "\n",
    "+ The total production must not exceed the maximum (here 14000)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Model: oil_blending\n",
      " - number of variables: 12\n",
      "   - binary=0, integer=0, continuous=12\n",
      " - number of constraints: 13\n",
      "   - linear=13\n",
      " - parameters: defaults\n"
     ]
    }
   ],
   "source": [
    "# -- maximum global production\n",
    "mdl.add_constraint(mdl.sum(blends) <= production_max)\n",
    "\n",
    "mdl.print_information()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "#### Express the objective\n",
    "\n",
    "* The objective or goal is to maximize profit, which is made from sales of the final products minus total costs. The costs consist of the purchase cost of the crude oils, production costs, and inventory costs.\n",
    "- The model maximizes the net revenue, that is revenue minus oil cost and production cost, to which we subtract the total advertising cost.\n",
    "- To define business objective, let's define a few KPIs :\n",
    "    * Total advertising cost\n",
    "    - Total Oil cost\n",
    "    - Total production cost\n",
    "    - Total revenue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "# KPIs\n",
    "total_advert_cost = mdl.sum(adverts)\n",
    "mdl.add_kpi(total_advert_cost, \"Total advertising cost\")\n",
    "total_oil_cost = mdl.sum(blends[o,g] * oil_data[o][1] for o in range_oil for g in range_gas)\n",
    "mdl.add_kpi(total_oil_cost, \"Total Oil cost\")\n",
    "total_production_cost = production_cost * mdl.sum(blends)\n",
    "mdl.add_kpi(total_production_cost, \"Total production cost\")\n",
    "total_revenue = mdl.sum(blends[o,g] * gas_data[g][1] for g in range(nb_gas) for o in range(nb_oils))\n",
    "mdl.add_kpi(total_revenue, \"Total revenue\")\n",
    "\n",
    "# finally the objective\n",
    "mdl.maximize(total_revenue - total_oil_cost - total_production_cost - total_advert_cost)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "#### Solve with the Decision Optimization solve service\n",
    "\n",
    "If url and key are None, the Modeling layer will look for a local runtime, otherwise will use the credentials.\n",
    "Look at the documentation for a good understanding of the various solving/generation modes.\n",
    "\n",
    "If you're using a Community Edition of CPLEX runtimes, depending on the size of the problem, the solve stage may fail and will need a paying subscription or product installation.\n",
    "\n",
    "We display the objective and KPI values after the solve by calling the method report() on the model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "* model solved with objective = 287750.000\n",
      "*  KPI: Total advertising cost = 750.000\n",
      "*  KPI: Total Oil cost         = 487500.000\n",
      "*  KPI: Total production cost  = 54000.000\n",
      "*  KPI: Total revenue          = 830000.000\n"
     ]
    }
   ],
   "source": [
    "assert mdl.solve(url=url, key=key), \"Solve failed\"\n",
    "mdl.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "### Step 6: Investigate the solution and then run an example analysis\n",
    "\n",
    "#### Displaying the solution\n",
    "First, get the KPIs values and store them in a *pandas* DataFrame."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": true,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "all_kpis = [(kp.name, kp.compute()) for kp in mdl.iter_kpis()]\n",
    "kpis_bd = pd.DataFrame(all_kpis, columns=['kpi', 'value'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "blend_values = [ [ blends[o,g].solution_value for g in range_gas] for o in range_oil]\n",
    "total_gas_prods = [sum(blend_values[o][g] for o in range_oil) for g in range_gas]\n",
    "\n",
    "prods = list(zip(gas_names, total_gas_prods))\n",
    "prods_bd = pd.DataFrame(prods)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Let's display some KPIs in pie charts using the Python package [*matplotlib*](http://matplotlib.org/)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAicAAAFeCAYAAABelb9JAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3XecVNX9//HX+85SpAjSUWnSRIpi7yWKDbvYfrF3oiaK\nxhIjJtZgL7EEQew9sYavsaAhURMVY4loVBSxIYiISN+Z8/vjnIHLMLs7uzu7d3b383w85rE79557\n72futM+cduWcwxhjjDGmVERJB2CMMcYYE2fJiTHGGGNKiiUnxhhjjCkplpwYY4wxpqRYcmKMMcaY\nkmLJiTHGGGNKiiUnxhhjjCkplpwYY4wxpqRYcmKMMcaYkmLJiTElSFIvSRlJR8eW/U5SJsm4mjJJ\nx4bnpGcCx75L0mf1fVxjkmLJSZFIOiZ8cG2as3xtSW9IWixp97Ds4lA2e/tJ0vuSLpXUtqp9Fjnu\nl3NiWSzpHUm/kqS6Om5TJKmDpKslfShpiaR5kp6VNLKCTXKvLeHyLKvseDtL+oukbyQtk/StpKck\nHVjjB1H58QaF13bRv7xj74XsbYmk/0m6WVKXYh+vAtU6/9UlqXs4f8MqOHZiiamktiG2tyUtDJ8T\n70m6UlL3OjrmaEnH1MW+TekrSzqARma1D66QaDwPDAYOcM49l1P2VGAR0AbYHbgQ2AXYvqJ91gEH\nfAGcDwjoBPw/4Prw/0V1fPwmQdJA4EWgI3AnMA1oD/wceFrS1c6587LlnXOfS1oLWFHD4/0OGAt8\nBNwOfB6OvTfwmKSfO+ceqvkjymsj4GLgJWBWkfcN/rV6ETATaIl/n4wG9pI0xDm3tA6OWZ/WxZ+/\nz4B3c9adSEI/JiVtALwArA88CvwJ/7ocCpwAHAhsWAeH/gUwF7i7DvZtSpwlJ3VEUhvgOWAYcGBO\nYpL1Z+fc9+H/8ZIeAw6UtJVz7t/1FSuwwDn3YPaOpD8BHwJnSBrr7OqQtSKpDHgMaAfs4Jx7M7bu\neuAB4BxJbzrnHs2uc84tr+HxRuETk0eAnzvn0rHV10oaATSryb6rOjR1n0w/65x7K/x/p6TvgbOA\n/YGH8wYltXLOLa7juIqhwprK8BymK1pfVySlgL8AnYGdnHOv5ay/EDgv37bG1IY169QBSa2BvwGb\nAAc5554tcNMp4W+fSvbdVdIkSV9IWirpa0lPxKvSQ1PSQElr1yR+59wy4A2gLbBalbmk5pJ+L+nj\ncPxZksZJah4r856kF/LELklfSXo4Z9mZkv4bqupnS7pdUvucbWeGJontJP07lJ0h6aiccnn7ZaiC\n/gKS9pI0NTSt/SjpGUkbVfOUVWUUvvbsynhiAhASv1OABcDvYnGt0eekGi4F5gEn5CQm2WM+75yb\nHDtWZ0kTw7lfEqru1ziupMMlvRnO0wJJ70o6I6w7Bp8MAWSbCtOSdqxB/NUxBf+l3ifEkX2ed5R0\nq6Rv8TWD2ccwXNL/hfgXSnpB0lZ5HutGkqaE5osvwpfwGp+X4Vhj8yyfKenOnGXtJF0v6bPw3vlC\n0t3yzX07Aa/jk7u7Yufv6LDtGn1OJLWSdG14Dy6Vby48u4IYb5K0f3hvLg3vtz0KOL+j8D+wLstN\nTACccz8551arXZV0SHidLJY0V9K9ktbNKVPp51h4rIOBnbWqKW8KpsmwmpPiawM8C2wGHOyc+79q\nbNsv/J1XSZm/AIOAm/BV9V2AEUBPVlWlHwhMAo4F7qnG8eP64D8of8gukCTgaWBbfNXuh/iq3bOA\n/sBBoejDwMWSujjn5sT2uQPQHYg3J4wHjsY3ddwYjnsGsImk7WJfri4c41FgInAXcDwwSb7G4YNY\nuXy/3tdYHhKbu/DP17lAK3wzwT8kDXfOFatpYt9w7HvzrXTO/SjpSeBoSRs45z6t6YEk9QMGAhOc\nc4sKKN8SeBnoC9yMbzI5BP8F2c45d3MoNwJfw/M8/lyBfx1uG7abin9NngFchn9tAGSfl7qS+57J\nPse3AnOA3wOtASQNDnEuAP4AlOMTw5cl7eiceyOU64o/JxFwBbAYOBmoTrNR7mutNfBP/HMzEfgP\nvtl0P3xzyXR8bdcl+PfWP8Kmr8b2l/u6fhrYKezvbWAP4GpJ6zrncpOUHfDvz1uBhcAv8c17vWK1\nt/nsF457X9UP2SeH+Pfyv/FNxV2BM4Ftw3vqx1C0qs+xXwF/DLFehk9Avy0kBtNIOOfsVoQbcAy+\nw9pn+A+xfSspezG+irY/vh9AL/yH32LgK6BlbJ9pYNNwv104xpgCYkkDRxcQ90vA+yGOjsAA4Kpw\nnCdzyh6Jb2veJmf5yeF4W4f7/cP2v8gpdwvwI9Ai3N8+lDssp9yIsPzw2LLPwjG2jS3rBCwBrso9\nt5Wck57hfmvge+C2nHKdgfnA7UV8bbwFfF9FmTNDfCPD/V7hHBwdK5P3seXsZ9+w3S8LjO1X4bjx\nc50CXsF/ibcOy64v4DEcHPa1Yx29v9L4PlkdgfWAw/B9En4Cuue8D18GlLOPx8PrpVdsWbfwOF+K\nLbs+HGuz2LKO4XWx8jUUlmeAsXni/Qy4M3b/92Hb/Sp5jJvlPuexdZOAT2P39w9lz88p9zA+6eqT\nE+MSoHds2VDyvEfzHHdaVc97rGwZMBufKDWPLd87HOvicL/Qz7H3gCnFfi3ZrWHcrFmn+Lrgk5Mv\nqign4H/4D9fPgNuAj/FfThX9QlsCLMdXdbavoAzOubudcynnXKG1JoNCHHPxv3jPAZ4EjsspNwr/\nS/gjSR2zN3yCI/wXB865j/EfUIetfLBShP/l9qTzzUbZ/f0AvJizv//gv3B2yTn+dOdc9pckzrnv\n8OdwgwIfZ9zu+A/Jh3KO7fC/+nKPXRtt8b8AK5NdX6OmuJjs9lUdL2svYLaLdY51vrbqJnwt4E5h\n8Q9AmwKbAuqK8J2K5+LfXw/gk90DnHPfxMo54A7n3MqahvD6GwE87pz7fGVB52aH/ewg308M/Dn5\nl3NuWqzcPOD+WsR+EPCOc+6pWuwjbi98EnJzzvLr8DU+e+Usf945NzN7xzn3Hv7cVfXeWZvCX0ub\n4z//bnWx/lLONyF+CGRHpRX0OWaaNmvWKS6Hr0W4AfibpO3DF3VFZQ/Cv/FXAF865yqdx8A5t1zS\necA1wLeS/gU8A9zjnKtNledn+NEAKXz1/oX4GoTcJKk/vlf+3HzhsXr/lIeByyV1D18cu+CreOOd\nFvvjR6zEm34q2h/kHwEyH1gnz/Kq9MN/2b1UwbEXVLRh+KLrnLP4e+dcRSNrFgK9q4gnO4S80C+C\nimSrzdtWWmqVXvikONcH+PPTK9y/Fd/cM1nS1/jO3o845/5WkyBDc1K7+LICXsMOP4LjY/wX87fO\nuf9VUHZmzv3O+Ga7j/KUzT7WHuH/XsC/8pSr6FiF6IvvFF0svYCv3ZpNdx/E1sfl+7FUyHvnRyrp\nA5cnJkf+c/whsB3U6eeYaUQsOSm+D4A98V96z4d+E19VUPYfrvL23jU4526U9BRwAL6N+RLgAkm7\nOOfeqWHMi5xz2S/pFyS9im+KuALf3JAV4atazyL/yIL4B+DDwJX4L7SbgEPxv77jX2YRvh35/1Ww\nv9wkqKLRCvFtKxotksq5H4WyR5K/Lbu8gv2A/xL7LGyfHaGyC74/Qz4fABtLWt8592UFZTYOf6dX\nctxCZPt6DC2wfEFz2Tjn5kraBP+a2yvcjpN0t3Mut4atEIfhmypWHoI1n6N83nCrRutUZknO/erO\n2ZPvdVSdfRTyWGqjolgqev0X8t7J50N8/6/1KvkcK3RfK9XR55hpRCw5qQPOuWmS9gcm4xOUHUK1\ncLH2/xm+Xfx6SX2Bd4Cz8R1Li7H/9yTdB5wi6ZrYF+oMYFgskalsHzMlvQ4cJukWfCfdx3NqF2YA\nuwKvxpp6ams++BFLblXnO1iz5mIG/sN0rnOuuqMAZgO75Syr7AP1GeAI/PNzRe5K+flw9gM+cLXo\nDAu+SU3S/4D9Jf3KVT2Edib5E5lB4W+8CaQc+Gu4Iek24GRJl4a4qzOM+FnWPId1aQ6+T9fAPOsG\nsWq+H/CPeUCecvm2nY+v/VtJUjN8x++4GcCQKmKszvmbCfxMUuuc2pPsSLPP19ykRp7Gv3aPBMYV\nEJPw5+nlnHUDc2Mq4HPMpjBowqzPSR0JX+BH4Jsuno21Z9eYpLUktchZ/Bm+KaBFrFythhIHVwHN\ngTGxZY8A60s6KU9sLSW1yln8MLA1fqKmTqw5D8Uj+AQ531DMlKR2ucsLkE06Vg5hDSMlchO3v+Gr\nrH8jPw9J7vE7VXQA59wy59yUnFuFzUD46vzpwPmSNss5jvCTpLXHd5oshovx53ui/DwVq5E0Qqtm\npZ0MdJMU7x+Uwo+6WQj8PSzrkOc474W/2dfeIvy5r7IfgXPu29xzWNhDqxnnXAbfFLW/Vh923xX/\nPp3qnPspLJ4MbC1p81i5zqFcrhnEXmvBqaxZc/JnfO3Z/pWEmU0yCumHMRn/3jk9Z/lZ+M6m1Rkl\nWJnH8M/zhZK2zl0pP3PsZeHum/gk8NSQoGXL7IVPAJ8J9wv6HMOfD+uT0kRZzUlxrVat6Zx7InyR\nTwSekbRHDWoI4vscgO88+gj+y64c32+lC/BgrFythxI75z6QNBk4Mfwyno8fCnsocJukXfAjOlL4\nD55D8J1M41Xuj+Dbla/GD/V8MecYU+UnfDs/NBk8h+9/MwDfWfaX+CGH1fEcvm/KnZKuxn9QH4f/\n0OwRO/ZCSaPx5+ctSQ/hm5F64jvu/TMcv9accyvkJ0Z7AfinpEn4D/L2+Cat4cA1LjYBWy2P94ik\nocBvgOGSHmTVDLF7Aj8LxwU/lPsU/NDhzVk1lHgb4FexX+UTQoIyBfgSXxN1OvC2WzWM+21888F5\noaPjMuDF0HG5GAptNqio3G/xtTWvSLoVH+vJ+CT83Fi5q4Cj8P3GbsTXuJyEP4e5U8tPAG6Xn0Dx\neXzz3O6s2SR5Nf41/Wh4/qfhn499gVNCB9UZ+KbPUyX9hP9y/le8A2+Wc+4p+Xk/LpefwTU7lHhf\n4Pqq+q8VyjlXLumg8Nimhs+eV/Dv08H419H3wG9D2fPwQ4mnhtddN/z76FN8Xzwo/HNsWjgXFwKf\nAHMKqbU1jUSSQ4Ua042cYb8568aEdU/ia6uyQ4k7VGefQAd8/4338b/6v8fPg3BQBdsVOpT4nQrW\n7Rj2Mza2LIUfzfMu/kP7O/zkURcCbfLs4x9hHxUOzcXXrLyOH6HzA/6D9gqga6zMp+QMbY7F/2LO\nsk3CeVmC/0X2S3KGEuc8xsnhXC7Cd+abCAyvg9dIR/yX1P/CuZuHr8HZO0/ZXrnPYXjdlFfjeDvj\nk7tv8InCbOAJwnDlWLlO+C/Zb8M5exs4KqfMgfhf49/EzustQJeccsfjO6wup4jDiit7f1WnHD55\nmIzv8LwQ/6W7ZZ5yg/GJ2CJ8snsBPsnNHUqs8Fr9Nuzvr/gOpJ8CE3P22R4/l8+scA4/D6+1dWJl\n9sHXVCyLP//4HxszcvbXCp/8f4HvvP4hcFaex5IGbsyzfI0YKzmva4fX39vhcS7CN8Ncmuc1MAqf\nfC9m1fTz3WPrC/0c6wI8hf9MSGPDipvUTeFFYIwxxhhTEqzPiTHGGGNKiiUnxhhjjCkplpwYY4wx\npqRYcmKMMcaYkmLJiTHGGGNKiiUnxhhjjCkplpwYY4wxpqRYcmKMMcaYkmLJiTHGGGNKiiUnxhhj\njCkplpwYY4wxpqRYcmKMMcaYkmLJiTHGGGNKiiUnxhhjjCkplpwYY4wxpqRYcmKMMcaYkmLJiTHG\nGGNKiiUnxhhjjCkplpwYY4wxpqRYcmKMMcaYkmLJiTHGGGNKiiUnxhhjjCkplpwYY4wxpqRYcmKM\nMcaYkmLJiTHGGGNKiiUnxhhjjCkplpwYY4wxpqRYcmKMMcaYkmLJiTHGGGNKiiUnxhhjjCkplpwY\nY4wxpqRYcmKMMXVM0oOSHijSvq6U9Gpd7NuYUmHJiTGm0ZOUkZQOf3NvaUljC9zPa5KuqKMYU5LO\nlfRfSUskfS/paUlb5hS9FNi7LmKoLkkDwzkckHQspnEpSzoAY4ypB91i/x8O/B4YACgs+6neI4qR\nJOBxYGvgHOBlYB3gV8A/JO3vnHsWwDm3GFicUKi5BLikgzCNj9WcGGMaPefcnOwNWOAXubmx5YsB\nJO0m6U1JSyV9JemSkDgg6UFgK+C8WI1LF0nNJU2SNFPSYkkfSBpdzRCPBkYChzvn7nHOzXLOveOc\nOx54HpgoqVmI40pJr1Vn55J2kjRV0iJJ8yT9VVKrsK6lpFslzQk1Ni9L2iS2bUdJD0maG3t8R0hq\nAUwPxT4M52NyNR+3MXlZzYkxDZyk1viagW5Ad6AT0Crc1or9jf8PsDTcluT8XYr/ZT4H+Ar4Gvja\nObeofh5RMiT1Ap4GbgH+HzAYmICvVbkKOAXoD7wCXA4+6Qlf8jOAG4H5wI7A7ZK+cM49U+DhjwDe\ndc5NybPuWuBFYBfgubCs4NoKSVuF7W4BRodtdwVSociNwB4hhq+B3wJ/k9TXOfcTMA7oDYwAvsfX\nOKWcc8sk7QBMBbbDn4NlhcZlTGUsOTGmxIVf7t2BoeHWB+hJin5kWB9os/oGOMrIhBs0W3kTzYlo\nFmpMy3GsIMMKHCtwlEO4iXIiylevWVVKPyFmk2EWji+Bz4H/httHzrnyOjwN9eEM4EPn3Dnh/keS\negO/Aa5yzv0oaQWwKNTAACubWS6L7edeSTsChwKFJicDgH9VsO6DWJnnKihTmXOBvzvnxsSWTQeQ\n1B44HhjlnHsxLDsemAUcg09oegDTnHNvh21nxfbzXfg7L35OjKktS06MKSHhV/hgfBIyDLEJEZuQ\noR0AZaRpT4YOlNEO0Q5YG5+etA1/WyJEilW/jCs8XKVllgELw+1HYCFtWEg/FtKPBZTzPbAkfIaI\nFSrTJ6R5C3gPn7C8B3zhnGsofRIGAa/mLHsF6CCpk3PuuzzbACDpTOAooCfQEmgOVKvppRLZfjE1\nPY+bAOMrWNcf37y/8nGHGpFp+PMBcCvwYKiBeR74i3PujRrGYkxBLDkxJkGS1sE3A+xMxG7ARmT7\ngq3DCrrTjC5A13BrT4qoyqSjOFqEW6e8a/1nxyJ8488cmjGHQcymP3M4nBUhxog5kl4E/h5u/yvh\nZCVf584qEwNJx+JrTn4FvIlP5y4CBlbj2B/hn/t8BsXK1MSSStZV9PhWngvn3JOhyWsksBswVdLV\nzrmCRjgZUxOWnBhTj0I1+g7ALiEZGQKItVnBBjSjJz4J6Qw0p1mCoRamNb6Rqc/KJWU4fJfTb4FZ\ndOFTDmE2h+MQEfMkTcGPRvk7ML2EkpXp+L4YcdvhmyzmhfvLWbO2aVvgJefcxOwCSf2qeeyH8J1e\nd802r8ScDXyDP1818S7+cY3Ls+4jIANsDzwBvoMssCnwVLaQc24ucBdwl6TX8U1dY/Hno/IaOGNq\nwJITY+pY+KI6mIhDgOFARBvK6UsZvfFdDddpAIlIoQS0Dzdfd1DGMuALYCYd+YyD+IZRZBARcyU9\nCvwF3y8iyX4rNwO/kHQt8Cd809qFwB9iZWYC20jqASwOScvHwEGSfoZ/lCeEbadTuLuBUcADks4F\nXgI64GtjdgX2d84tr+Hjuhx4S9L1+A6+mbDPu51zP0iaAFwv6Sd8EnQhvtbkHgBJl+ObqKbj09G9\nYo/tG3wD4F6SvgOWOucW1jBOY1ay5MSYOiBpED4hORwYTIoM/REDUEhGylZWqDcFLYB+4QYpluO/\nxj+mM+9zMgv5BRELJP0Fn6i84JxbWp8hOuc+l7QPvobhF8A8fMJyTazYH4BJwIdAS0ndQ5mhwJ/x\nXYrvB+7A17oUemwnaX98LcmvgdvxzTGvANs7596sxeN6X9Ke+CTlFPxIrNeAO0ORMUAaeADfa+nf\nwB6x0Vnl+NFKvfANeX8HTgr7XirpLOAC/Hl7nhKZIM40bCqdGlVjGjZJQ4BDiDicDAMoI82GRGyE\n6IfvImnW5PC/vz8A/ks58ylDLMbxFPAY8HQtag2MMQ2QJSfG1EIYXXMYEaeRYTOakWYQKTYC+kIj\naqypHw6Yi09U3qecOZQR8T0ZJgB3OOc+STZAY0x9sOTEmBoIzTanIo7H0Ya+ZNiciP5YY2kxzQGm\nAf8hzXJSiCk4bgeetNoUYxovS06MKVCYrvvAUEuyPWtRzmaUsSm+66KpOyuA94E3SfMlqVCbMh5f\nm/JpwtEZY4rMkhNjqiCpLXAaEb8mQwd6kWYLUmyI1ZIk4Vt8bcrbpFlOhO9Ae6VzblqygRljisWS\nE2MqEOYkOYOIs4G2bErEVvg5SEzyluNn8Pgn5fxAGeJFHJcDL5fQ3CnGmBqw5MSYHJI6AmciziJi\nLTYnYlsIE8ibUpPGd6CdurID7atkuAg/MZp9wBnTAFlyYkwgqQswBvFLUrRgi5CUtE06MlMQh5/v\n9CXSzCYVkpQLnXMvJxyZMaaaLDkxTV7o6Hom4mJSNGdrUmyDnwvTNDwOP2frlJCkiKdwnG3DkI1p\nOCw5MU2WJAH7EnETjp5sidgJaJV0ZKYoHP7ayM9Rzk+A43rgMufcj8kGZoypiiUnpkmStBHiJhy7\n0ocMexHRJemoTJ1YDrwK/IMMjvlkOB+Y5JxLJxyZMaYClpyYJkXSOsDvgNNpT4a9KGMANKnr3DRV\nC4AXcLyHiHiXDGc456YmHZYxZk2WnJgmITThHEnETUS0ZWdSbI3NU9IUfQH8H2m+JgXcBZzlnPsh\n2aCMMXGWnJhGT1JXxHgc+zEExx7IRuA0cRngbXySkuY7MhzrnHs26bCMMZ4lJ6ZRk3QIEeNpQVv2\nDRfkMyZrAfAkGT4lAiYCZzvnFiQclTFNniUnplGS1AlxC45DGYRjH2RDg01eDngLeJY0aeaQ4Rjn\n3PNJh2VMUxYlHUChJD0o6YGk48iStIekjKTmdXycKyW9WpfHaGwk7UfEhzTnYA4GDrXExFRCwGbA\naaToSVfgOUl/CtdUMsYkoODkJHwRp8Pf3Fta0tgC9/OapCtqHnJJKVq1k6QW4VzunrPqUmDvYh0n\nCfX1nEtqKWk88CR9WYfTSTEUG4ljCtMeOIaIkUAZJxDxtqTBSYdlTFNUnbEK3WL/Hw78HlYbhPlT\nsYKqL5KaOedWJB1HkPcr1Dm3GFhcz7E0OJL6EvE4YjB7A5sSWVJiqk3AFsAGpHiInnzHm5KOd849\nmHRoxjQlBdecOOfmZG/4bmTOOTc3tnwxgKTdJL0paamkryRdEoZxIulBYCvgvFiNSxdJzSVNkjRT\n0mJJH0gaXZ0HIukUSd9IOkTSDElLJD0jqVuszJXhV/xoSZ8B88PylpJulTQnbPeypE1y9r+/pI9D\nfH8DeuSsv1LSaznLzpP0QZ44p4fz86Wka8Kqz/A1Mc+GczM9lP9DfL+SIkmXhnO7NJzrn8XWDwzb\n7ytpqqRFkt6StFkV56+DpImSvg2P8W1JI2LrDw9xLwvn94yc7c+U9EmIabak+8LyvM95ZbFUl6Sf\nA2+zNoM4kYjNsNoSUzsdgZMoYygtgAck3VTXTbjGmFWK2udEUi/gaeBlYBhwOnAa8OtQ5BR817M/\n4mtiuodkpwyYARwADAKuAK6RtE81Q2gPnAUcCmwfjnFvTpnBwJ7AfvgvTYAb8U0nRwCbAl8Bf5PU\nJjyuvsAjwMPAxsD9wOV5jp+vmWflMklnAdcAN4U4DsQnJeB/rwlfK9UtxJ/dPr7f84DR+HM7DJgK\n/FVSz5zjXgZcEuKdFWLOS1IEPA9sAhwGbARchB9wiaRtgfuASSHuK4CrJB0a1m8PjAPOBfrjz2+2\nn0xFz3mtSUpJugy4TynacDJldC/Gno0BmgMHIvYGIk4n4h+S1ks6LGOagmJPQXUG8KFz7pxw/yNJ\nvYHfAFc5536UtAJYFP+CCrUul8X2c6+kHfFJxjPVOH5z4GTn3H8BJJ0A/EfSkOwyfEJ2pHNuYSjT\nHjgeGOWcezEsOx7/hX4McAs+wXrXOffbsI+PJQ0HflloYKH26ALgCufc7WHxDOCN8P/c8PeHKr68\nzwYudc49Hu6PkbQr/tz/OlbuSufcC+HYvwfelNTTOTcrzz73AYYA/WPrZ8bWjwGecc5dHe5/ImlY\nON4j+FqkBcBk59xS/DRXbwNU9JzXlqT2ETzoYM8jgXvTwOPAz4t1BGPwPxe2BLojHmZTFvOupFHO\nuZeSDs2YxqzYo3UGseoXc9YrQAdJnSrbMDQLTJM0V9JC4GggtzagKktiSQjOuXeAJSGurBnZxCTo\njz8Pr8a2WwZMi223IfDvnGO9RvWsD3QCplRzu5UkdQY6kP8cxx+jA96L3f8G/zFbUXPKxsCnFSQu\nhH2/kueYG4b/J+OTq5mS7gpNQC0qeyy1IalvGbzVGkb8H3APcD6gj1nzzBhTDD2A0ZTRk/bAi5J+\nkXRIxjRmxU5OxJpNG9nW/wpHtkg6Fl9zciuwG/7L8gF8TUh1VHSM+PJFBcYXfyz5HleuDGv2dGgW\n+39JFdsXopBYs+IdfbPrKnq+q4qtsueVMGnVMOAoYA6+yestSUUfwCtpizJ4vSf0eAtSe4TllwI7\nAanngC+LfVRjgNbA0URsjYBbJF2e7U9njCmuYicn04HtcpZtB8xzzs0L95cDqZwy2wIvOecmOufe\ncc59CvSrwfFbSRqSvSNpY6Al8EHFm/ARPrHI9vFAUkt835PpYdF0YOuc7bbJuT8X1ujxMDz7j3Pu\nO2A2sGsFcZTjE4Dcc7NSaBaZF481Fktlj7Eq7wIb5Om3kjU9zzG3ix/TOZd2zj3vnDsXf+42BHYI\nq/M959Umaa8Ipg6Hdv+GsvgLpAzfvtQFSN0NLK3t0YzJI8L3qPID/n8D3CmpWWWbGGOqr9jJyc3A\nQEnXShog6WDgQuDqWJmZwDaSekjqGJZ9HJb9TFJ/SX8Ahtbg+MuB2yVtLmkL4A580vN+RRuEC35N\nAK6XH2k0GLgTnyjcE4rdCgyTdFmI72h859m4KUAPSb+S1FfSmcAuOWV+B1wgP1qor6TNstXD4fLt\nXwK7SeoRfX+mAAAgAElEQVQqqV0FIV8D/FbSQeEcX4cf0n1zrEy1fs05554D3gQel7SLpN6S9o6N\nAroGGCnpXEn9JJ0InEh4XiUdKOkXkoaFTtHHAWn88wr5n/NqkXSc4JmR0PxlSOVrI+wMPAloBf4Z\nNaaubAscBEQcg3g623neGFMcRU1OnHOf4ztX7gi8gx8FczP+yy3rD/jmmg+BOWFY6c34fgt/xvca\naIlPLKprfjjmI8Df8U0MRxWw3Rjgr/impDfwNSC7O+cWhcc1A98597DwuI7GJ10rOefeBX6F77D6\nFn5Uyw05ZcbjR9v8CngfeALoHStyJn4U0Swq7tNyNT5ZuinEsgOwt3Pui/ih8mxXVbPUfvgalEdC\nbJcTkhzn3L/wXU2PAf6L79j7a+fco2Hb+fhzMyVsexS+g/GMsD7fc14QeRcBd54C0V8galVJ+S3w\nJ4fv8JmKMXVlGPBzRBm7ETG12EPkjWnKGs21dSSdAlzsnFs36VhMcYQhzn8ERl+Oz4gKrRI6AbgL\nyByE/xIxpq58A9xLOUv5kgy7hmZpY0wtNJhr65imRVIkGC8YPRHfuF+dtqpb8O2CZY/je+kYU1e6\nAydSRjvWJ+JfkvonHZIxDZ0lJ6bkSErJX77++Hvwk9BUV0t8q05rB9EEfHdjY+pKB+AEymjPOqGJ\np0/SIRnTkDWaZh3TOIQakzuBo+8H5fY6rq7ngT0A15OaZTnGVMePwJ2U8yOzybBtTl8wY0yBrObE\nlIzQx+R2ipSYAIwgTD08i1pMf2dMgdYGjqWMNnQLNSjWB86YGrDkxJSEMJnVjYKTJhUpMck6H9gX\nSE3FXzDAmLrUHjiOMlqxPhF/l9Q16ZCMaWgsOTGl4iLg9Nvx45WLKcJf/bEnkHoA+KnIBzAm1zr4\nBKUlvYl4qarLdxhjVmfJiUmcpGOA318KnFxHx2gHPAU0S4PuIFxv2Zg61BHfxNOCAUQ8XxeXczCm\nsbLkxCRK0gjBxBPJmdWuDgzBz33iFgCPVl7WmKLoAhxNioihiHtDvypjTBXsjWISI2lYCp7YHXQb\n1Zxzv4YOA84C9AHwej0c0JjuwChSOA4k9M82xlTOhhKbREhavwzeHAyd/gmp+rwwyQr8RY/+BaRP\nYc3LNRpTF14FngPgOOfcXYnGYkyJs+TE1DtJbcvgX11hwBtQlkRuMBvYGJjXHNLn4K/8Y0xdcsDT\nON4iDezqnJuadEjGlCpr1jH1Sr7R/Z4WMPC5hBITgG7A4+CvY31nQkGYpkXA3oheCPGUpH5Jh2RM\nqbLkxNS38zNwwAOQ2ijhQLYlXDZ6Nv6a2MbUtTLgMFK0pzURz0paJ+mQjClFlpyYeiPpKODy3wD7\nJR1McBrwcyB6HZiecDCmaWgFHEkZZfRGTAoTEBpjYqzPiakXknoQRdPJZNq0l3jVOQYlHVSwGNgC\n+EhQ/kv8BFrG1LUPgIcBOM05d2uywRhTWqzmxNQ5Sc1IpR6lffuW/O53/NC2LUNSKSYkHVjQCj9B\n21oOojuAdMIBmaZhED4rFjdIGpZ0OMaUEktOTH24Aue25JJLythpJ5gwgUz//pwkcSilMVlrX+BB\nILMYuD/hYEzTsTvQiYiIx2wGWWNWseTE1ClJOwPncNJJYvBgv7BzZ7jpJjjoIB4FNgDmJBfiSiOB\nsQCfAjbI09SHZsChpIjoS+ifbYyxPiemDklqSyo1nUGDunPDDSlSqTULTZkC48bRPJ3miXSaveo/\nzNVkgL2BF4D0cUCvZOMpGW8AbwI/hPudgZ2A/uH+JODzWHkBmwH7VLLPJ4C3c5b1A44M/5fj29v+\nB7TBZ48bxMr+E1gIib9oiuEt/GOFw51zDycbjDHJs+TE1BlJd9C8+fFMmhSx7roVF5w5E377W/jm\nG87JZLi63iLM73tgE+DrMkiPwXdKaeo+wiccHcL9t/Eznp6KT1Tuwl/o7mf4ycbA1wq0qGSfTwCL\ngANi25QBLcP//8YnRIcCHwOvAL8O6+YD9wGn0Dgm0HPAYzimsxjHUOfcZ0mHZEySrFnH1AlJI4ET\nOf30yhMTgN69Yfx42GEHrgE2x4+gSUoH4EkgVQ7YFYy9Afhako7htis+KfgyVqYZ0Bpfy9GGyhOT\nrFTONi1j674DBuKTny3wiUz2hfFXYASNIzEBn/jti2hLCyLutuHFpqmz5MQUnaSOpFJ3scUWGfap\nrF4/plUruPhiOO00pkUR3aJojRr/+jQcn5cwnzCVrFkpA7yHv0hRj9jy94CrgFvx7WIrCtjXTOBq\n4GbgGVbPSrsCs8J+PgHa4mux3sHXsGxYi8dQiloC+1NGhh2AE5MOx5gkWbOOKTpF0UOstdYo7r47\nRadO1d/Be+/B2LHoxx+5MZPhjOKHWLBfALcDbj9g0wQDKQXfAhPxfUGaAwezqs/JNKA9PoH4Fnge\nWA9/GeiK/Bdf27IOvi3txbDfE/E1CWngWXyTTmtgT6ATMB44Ft/k8z6+qmv/cOzG4Ekcb7MYx0Dn\n3FdJh2NMEiw5MUUlaQ/gWS64AHbfveY7+v57+P3v4b33GOkcT5FMNd9yYHvgP4LyU/G/5puqNLAA\nWIqfQGwacBy+2SXXZ8A9QHUmtZsP3AgcA/SpoMwT+KtIt8cnMyfhO8bOofJEqCFZAtxMOUuY7DJu\n/6TDMSYJ1qxjikZSC1Kp2xg6NM2IEbXbWYcOcN11cPjh/BVYn9W7N9SX5vhWnbUdRHdSWFNFY5XC\n11Ksi+9z0g3faTWf9fCdPL+vxv7XwTfbVLTNp8BcYEt8c1B/fM3LYFYfKdTQrQWMpAzHfpJK5UoP\nxtQrS05MMY0hk+nNWWelKEZ/vlQKTj4ZLr2Ub1q2pE8qlUj3j/WAPwMsw49KMZ7DN/HkMxvfNNOm\nGvtbgK81yLdNOf7ijPuG/WZY1VE5TePrtLwR0JcMEbfZ5GymKbLkxBSFpJ5E0VhGjRJ9KqqTr6Ht\nt4cJEyjv0YODoohfFHfvBdkZ39eTr4DnEgggaS/iayd+wPcpeQFfezEMX9Pxd+DrsP5DfHVTL1Zv\nBrs5rAPfXvYcvjrsB3ytyEP4kUD98hz/ZfyIoW7hfk9809K3wOvhfmMiYCQRohtwUdLhGFPfrM+J\nKQpF0WOsvfb+3H9/Ga3r6Ife0qW+qef55xmCn/Zi7bo5Ul4OP+XG40D6CPww16biSXw/kp/wQ4S7\n4jvjbICv8fgLvsllOdAOf92YHVh9OPHv8R1XN8E3jz2Er2FZiu/M2g/YBd/5NW4O/gJ5p+KbccA/\nGZPxI4Q64jvndqDx+TvwEmlgsHPuf0mHY0x9seTE1JqkEcBzXHgh7LZb3R7MOXjmGbjxRlo5x5RM\nhq3q9oir+Qk/aOfTCNJnUr/ZkWl6VgA3Us5PPOmcG5V0OMbUF2vWMbUiKUUqdRNDh6bZddf6OCDs\nuy/88Y8sXmcdto4i/lD3R12pDfA00CIDsgnaTF1rBuxKGXCwpC2SDseY+mLJiamtw0mnN2T06OJ0\ngi3UhhvCxIkwfDgXALtRcd/MYhuInzndLcRfytiYurQx0JFyxFU2c6xpKiw5MTUmqRmp1BVsu22G\nQYPqP4B27WDcODj2WF4EukvU1wVJDgTOB/Qx/hozxtSVCBhBGY6d8ZP2G9PoWXJiauM40umenHBC\ncq+jVAqOOQbGjeO7Vq3on0pxf06R2/E/PtuF27b4iUcrUg5cgu+fuRZ+Kvu/5ZS5P9xSgLKjTrLm\n40emLKvZQzJmDQOB9UkTcZUk+9w2jZ69yE2NSGpJKnUJu+7q2GCDqjeoa1tuCRMmkO7dmyMljo2t\n6gGMw09oOg1/4dz98SNR87kQf12dW0KZU/A1Je+E9fPwE5NeB/wfftRnNAk/6gRWXZSukAvfGVMI\nASNIkWFj4JCkwzGmrtloHVMjks5CupZ77xXrrZd0OKssXw5//CM8/TT98BOY5hth2hG4Bj/7eq71\n8BNLnBpbNgo/eek9wBv45ObrsG53/DQgmU7AjsB04PDaPxRj1nA/GWYwiwwDnHNNeb5i08hZzYmp\nNkltSKUuYu+9SysxAWjeHMaMgfPO45OyMtaNIl6Orc7gp9dYDGxTwS6WsWalx1r4S7iAnzV9Mb4m\n5XtgBv4SMnyHrzXZuwiPw5h8diMiQ2/g/yUdijF1yZITUxPHkcm058gjk46jYnvuCbfdxrLOndkl\nijgFP89XC/yVhh8HNqxg0z3wTTaf4Of6eh4/x9g3YX174G7gKGBr/AVyr8cnLVqGv1ru7cCt+FoU\nY4qlK9CPDBHn2sgd05hZcmKqJcxrcg477wzdulVZPlH9+sGECbDVVozHz37+T2A0cDSrZlLPdSM+\n0dgQn8z8Ejge3/k1a3/gXeAjfBPQy/jOtkMApgLb4aeTfRJYVMwHZZq8bYnIsBH+8ovGNEqWnJjq\n2p90uieHHdYwfrW1aQOXXQYnn8xbEntKHIkfvXNjBZt0wteULMZfTuYD/IzqFV0xaDlwGjABuJrQ\nQXYyvoqlI/56PMYUSx+gC+WIXycdijF1xZITUz2p1LkMGZJmYAO6sEwUwRFHwLXX8kPbtgxJpfiS\nqkf6Nge642cQ/zNwQAXlLsF3M9kYX+veFnBL8L1n41fPNaYYBGxHGY7dJQ1OOhxj6oIlJ6ZgkrYm\nnd6Kww5LVV26hEyYAO+9B927w9ixZNq350N8R9YMvonnN7Hir+P7pHwG/APYC9/3JN/P1PeBR/EJ\nCqxqCjoAYBb+onUl1mfYNAKDgdaUA2OSDsWYulCWdACmAZHOpmvXcrbZpmG9bubPhyuvhHnzoHVr\n6NMHBg9m6tSp9JHo4dxqb4SlwG/xyUkbYCR+uvp81/g7BbgBP5oHoCVwF77TbQtgeQbcHHx1ijHF\nUgZsQxkvcrSkC51zs5MOyZhisnlOTEEk9QA+55e/FAcemHQ4xTFlCowbR/N0mifSafYq8u4X4GeX\nnZWC9Fn4TMeYYlkCXEuacq50zl2UdDjGFJM165hCHUPz5o499kg6juL52c9g/HiWd+vG3lHE2UXe\nfTvgKaBZ2q5gbOqAv7ZCiohTJTWs2kxjqmDJiamSJJFKncjOO0e0apV0OMXVqxeMHw877sh1wOb4\nUTrFMgTfzOMW4DunGFNMw4EMnfAX5jam0bDkxBRiB9LpXuxV7IaPEtGqFYwdC6edxrQoolsU8XYR\nd38YcBagD/C9bY0plu5AR8rx/bqNaTQsOTGFOI4uXcoZNizpOOqOBKNGwQ03sHDttdk0iiqcB6Um\nxuGvhpyazKqpZo2pLQHDKUMcLClfn21jGiRLTkylJLUhig5n5Mgyoibwchk6FCZOxA0dypkS+1Cc\nriLNgMfwc7KlJuFnbjOmGIYCjub461Ma0yg0gW8bU0uHkMm0bFQdYavSoQNcey0cfjh/BdbHT1lS\nW93w86ewHLizCDs0BnzP6z5kUN6LbBvTIFlyYioXRUczfHiGrl2TjqR+pVJw8slw6aV807IlfVMp\n/lyE3W6LnxeF2cDkIuzQGIBNiHBsL6lX0qEYUwyWnJgKSepIJrMjO+/cdF8n228PEyZQ3qMHoyRG\nF2GXpwE/B6LXsasWm+LYECgjDZTwpcKNKVzT/dIxhRgJRGy3XdJxJGu99eC222C33bgdPzz4x1rs\nTsB4wvfJo8D8IsRomrYWwCBSRPw86VCMKQZLTkzFpAPZcMM0HTsmHUnyWraECy6AMWN4P5WiWxTx\nWi121wo/QdtaDqI7gHRxwjRN2EAgw6Awm7MxDZolJyYvSa2Q9mLHHRvWRf7qkgT77gu33MKSDh3Y\nNoq4sha76ws8CGQWA/cXJ0TThPUFhIOiX4nBmHpnyYmpyG5kMi2afJNOPgMH+isdDx/Ob4BdgfIa\n7mokMBbgU2BqkeIzTdNawPpk8CPgjWnQLDkxFTmA9dYrp2fPpOMoTe3awbhxcOyxTAG6Scyo4a4u\nBvYAUlOAz4sWoWmKBpACRkhqkXQoxtSGJSdmDeFaOgewww52MbHKpFJwzDEwbhzzWrViYCrFfTXY\nTQQ8AKwLpO6luBf3MU3LAMDREtgx6VCMqQ1LTkw+G5FOr8PmmycdR8Ow5ZYwYQLpPn04SuIYwFVz\nFx2AJ4FUOWBXMDY11QVoQzm+xdCYBsuSE5PPzqRSjo02SjqOhqNbN7jlFthnH+7B/4D9vpq7GI7P\nS5hPmErWmGoSMJAyIvZLOhRjasOSE5PPTgwYkGGttZKOo2Fp3hzGjIHzzuOTsjLWjSKmVHMXRwOj\nAb0HvFX8EE0T0B/I0EfSBkmHYkxNWXJiVhP6m+zKppvaEOKa2nNPuO02lnXuzK5RxEXV3PwGYHOg\n7GlgTvHDM43cqj7sWycYhTG1YsmJybUh6XQHNt446Tgatn79/HDjrbbiMmA7YFmBmzbHt+qs7SCa\nCKyoqyBNo9QKaM8KLDkxDZglJybXzkSRY8iQpONo+Nq0gcsug5NP5lWJbhLvF7jpeuAvNLgMuKuu\nAjSNVk+aEWGTFJkGy5ITk2tH629SRFEERxwB117LD23bMiyVYnyBm+4MXAXwFfBcXQVoGqX1AMdQ\nm+/ENFSWnJjVlZVtw5Ah1t+k2IYPhwkTyAwYwCkSoyhstPAYYBSQehX4X51GaBqT9QFHM/wgMGMa\nHEtOzEqS2lFe3ov+/ZMOpXHq3BluugkOPpg/A30kZlexiYBJwAZA6mFqdzlk03R0BSIywFZJh2JM\nTVhyYuL8ryxLTupOWRmcdhqMHcusZs3oGUVMrmKTNsDTQIsMyCZoM4UoA7pbcmIaLktOTNwmlJVl\n7Ho69WCXXWD8eFZ0787IKOLsKooPBO4D3EL8pYyNqUoPyojYIekwjKkJS05M3DB69cqQsi4n9aJX\nLxg/HnbckeuAzaj8sjoHAucD+hh4tT4CNA1aNyDD+pJaJx2KMdVlyYlZpaxsOP3728X+6lOrVjB2\nLJx+Om9FEV2jqNKJYS8FdgJSzwFf1k+IpoHquPK/fglGYUyNWHJiAJAUkU4Pok+fpENpeiQ4+GC4\n4QZ+ateOzaOIGysoWgY8gr++W+puYGm9RWkamg4r/7NOZKbBseTEZK2Lcy3o0SPpOJquoUNhwgTc\n0KGcKTGS/H1fO+OvYKwVwIR6jdA0JK2A5qSx5MQ0QJacmKzeAHTvnmwUTV2HDnDttXDEEUwG1pOY\nlafYFsCtAN/hMxVjcgnoiMOSE9MAWXJisnoD0LVrslEYSKXgpJPg0kuZ3bIlfaOIx/IUOwk4Hoj+\nA7xbvyGaBqITZURsmHQYxlSXJScmqzdt25bbtPUlZPvt4Y47KO/Zk0MkTs1T5BZgKFD2ODCvfsMz\nDYDvdzIg4SiMqTZLTkxWb7p1SzoGk2u99eC222C33fgTMJjVJ4ltiW/Vae0gmgCUJxGkKVkdgAwd\nJa2ddCjGVIclJ8aLog1Yd10bRlyKWraECy6AMWOYnkrRLYp4Lba6F/Ao4JYA9yQToilRq0bsbJBg\nFMZUmyUnxouivlZzUsIk2HdfuOUWlnTowLZRxOWx1SOAywBmAVMSidCUojYr/+ucYBTGVJslJ8ZL\np7vT2T6/St7AgTBhAmy6Kb8FfsaqlpzzgX2B1FRgRlIBmpKyqgtZx0pKGVNyLDkxSFoL55qxtjVL\nNwjt2sG4cXDssbwEdJOYgX8z3wv0BFIPAD8lGaQpCS0A4Yg38BjTAFhyYgDaA9CmTRXFTMmIIjjm\nGBg3jnmtWjEgleI+oB3wFNAsbVcwNvi5TlqSxmpOTANjyYkBWAew5KQh2nJLmDCBTJ8+HCVxNLAR\ncBfgFuB7ypqmbS0clpyYBsaSEwPZ5KRt24TDMDXSrRvccgvssw/34ie1GAGcBegD4PVEozNJa02E\nJSemgbHkxIA16zR8zZvDmDFw/vnMKCtj3ShiL2BbIDUZ+Cbh+ExyWpECOiUdhjHVYcmJAas5aTz2\n2ANuu41lnTuzexSxMf4nc2oSsDzh2EwyWgERXZIOo7okPSjpgTrY72uSrij2fqtx/G8knVzsso2N\nJScGoD2plKN586TjMMXQr58fbrzVVtyKH7CRXg7cmXBcSXoduAE/GcwdwFeVlJ0DPBzK/w74V54y\n7wLXAeOAv+Wsmw/cDCyrVcTF40fstKvuZpIyktLhb+4tLWlsgftJNBlIiqRTJOWrsxwC3F3gbqpT\ntiTV9Pm35MQAlJFKOaSk4zDF0qYNXH45nHwyX0h+ErfZwOSkA0vAf/EJxM7AKUA3/JjrRRWUX4Ef\neLsbkK8ycTF+SNQewFH4ROWj2Pq/4jv9tKh96EXhP+VrMvtzN6B7+HsmsADoGlt+TXECbHwkleHH\nSrncdc65ec65glLX6pRtbCw5MQBCWuNNZBo4CY44Aq67zjfZpVLodWB60oHVs9eAzYFN8POk7gM0\nA/5TQfn18MnFECCVZ/18/EWNBgPr4q/n/V1Y9y4+DSil6wDXMDlxzs3J3vCJiXPOzY0tXwwgaTdJ\nb0paKukrSZdI/peOpAeBrYDzYjUuXSQ1lzRJ0kxJiyV9IGl0deIL+3lY0peSFkl6W9JBOWXaSHpA\n0k+SvpB0es766yS9lGffH0o6O3Z/dFi2RNL7kk6MrRsYHtvBkv4haQlwMHAr0DX2uM8N5VdrqpF0\nuaRZ4fzNknRVbN3KspJahH0dLenp8Jg/lLRnTuwHS/oknNe/STohbFdh1bikDpImSvo2bPe2pBGx\n9YdLmi5pmaQZks7I2f7McMylkmZLui8sz/v8VxRHnCUnBux10Lhtsolv5hkwACcRPYL/gm0K0vjO\nwH1iy4S/0syXNdxnB3ztymx8LcrX+PqEJcDLwN6hnMPPM5MJcWRv5eG2ItyWh9uy2G1puC0Jt8Xh\ntih2+yncFobbj7HbgnD7IRyjZjUnVZLUC3g6PPJhwOnAacCvQ5FTgLeAPxJqXEKyU4afx/gAYBBw\nBXCNpH2qcfi1gFeAvfCp5F3AQ5KGxcrcBGyBf1b2wqemg2Pr7wd2kLRu7DFtDvQDHgj3TwDOA87B\np51jgaskHZITz5X4hr4NQ1znAXPxr47u+Ma+1Ug6En+OjgvHPJiqfz78DpiEvyD5S8ADktqE/Q0A\nHgqxD8M3CV1KnhqcWAwR8Dw+fT8MPxvBRYRZkiRtC9wXjjkY/1xdJenQsH778LjPBfoDewKvht1X\n9PxXyS70ZsDXnCQdg6lLnTvDTTfBn/5E5rHH/MfkWlTykRVUtL46ywupk6vucaqzPoPvQxJ/iWcn\np7ukin044Nlwy3V77P97Y/9fV0VMSVDo9F58ZwAfOufOCfc/ktQb+A1wlXPuR0krgEXxL6VQ63JZ\nbD/3StoROBR4ppADO+c+xycfWTdIGgmMAt6VtA6+4W1/59xUAElHAV/E9jFN0sfAEcC1YfERwN+d\nc9n+Ir8DTnfOZeP6XNImwKmsPpPQVbEySPoRyDjn5lbyMHrge0C95JzL4FPmN6p46OOdc38Jx/gN\nPgHYFJgKjAb+45zL9gf6RNJwYEwl+9sHn9z1d87NCstmxtaPAZ5xzl0d2+cwfAL6SHgMC4DJzrml\n+PP7NkBFz38hLDkxYM06TUNZGZx2Gmy0EbzwQtLR1I/ly+HNN2HI0NVHo82cCQsXwtChlW8/bRqs\nuy507155uQUL4PPPYfBg+M9/YMAAaNYM3n0XNt3U/5+UL76Ar79eUUd7H8SqX8lZrwAdJHVyzn2X\nZxvANwXgk4ee+Iay5rDaBbcrFfp1XAQchG9gax5u2e7O/fG1witn+nHOzZX0ac6uHgB+DlwbmqMO\nB34bjrEOvqHvPq3+Ay6FrzuLm1Zo7DEP4WuaPpX0LL7H0l9DolKR97L/OOfmS1oOK0djDWTNmY2q\nmuloY+DTWGKSaxC+ViruFeD48P9k/PmaGR7Ds8Djte0rY8mJAYis5qQJ2WUXf2sKysthzz3hsMNg\nu+1WLf/DH2DRIrj00sq3P+IIGDkSDj644jIrVsBJJ8H11/vLCpxzDtwcavBHj+b/t3fncVbWdf/H\nX59zHbbcAsMlH+Jy58ItZpniXmZ1a2bdomlmaZZpkprlvuSGpUYumJi4guLyUzE1vMEAFRVRjBAR\nZZV9EwYQZoZl5lzX5/fH94xzGLbZrzMz7+fjcR7OOec657wHinnP9/u9vhe9esERRzT8e6mvgQPh\nqaea6kpLm1r0WfWPyZamEs4hjJxcAownTExdT/jhWlt/BM7Pv8cUwmTXAEJBqVWOvCeBm8xsf0IR\n6QI8n3+uqtGeDXxQ43W5gq+dzS+x3ix3n21mXwH+h7AE+yFgipkd5+6by12zaDrVU/Nb+vvYnLVb\neX6L7+nuq/IjKccRVmv9GbjezHq6e53/TKporYGApnWktcpmwyjGhAnVj7mH+wccsPnX1cXjj8Ph\nh4dTuJME4rj6uVwuPJam8Pnx1g6rp4+Bo2o8dhSw3N2X5+9XsPHS4iMJUxmPuPsH7j6LsOaiLo4E\nhrj7s+7+ITCPMFpSZTphAu/wqgfyizH3LnyT/GePA34OnEkYuVidf3o+Ybnz3u4+q8ZtcyMNVTb1\nfW/E3de5+z/d/XeEkvIt6lbSCk0FetZ47NCtvGYSsLeZddvM8x8DR9d47ChCIQTA3WN3H+nuVxKm\nmPYHjsk/Xas/h5o0ciIAFeRyaifSOp12Whgp2Xdf6N4dnnsO1q8PIyoAt94a1uScd164n8uFaZ+q\nr0tKYOZM6NQJdtttw/eePRveeAMeeijc79YtjJ4MGwadO4cplf1TPnWnacvJvcBvzexO4AHCIs3r\ngNsLjpkDHGFmuwNr8qVlBnCKmR1HKADn5l9bl3PJZgAnmFlPwtLgK6F6bU1+ymMwcLeZlRKWB9/G\npnegeQq4jHDtzHML3sPN7GbgVjNbA4wiTEH1BDq6+335Qzf17+ccwvTW0eRHdvJrMj6XX2ybI6wz\nWUeYXiqlYF1MHd0PXGhmfQgroQ4lFC7YzAiSu48ws/HAC2Z2OTCbsCh2nbu/Rjhl/K382Ub/IJyU\n/wVzmtIAABshSURBVGvgF/nvoRdhwe8YwtqTXoT/vc0o+HOo+fe/VSonArCa9esjkiT8wyrSmnz7\n22FNyMCBsHJlGOHo2xe+GK7aQEkJRAW/2JWUwPnn8/lo4rPPhttBB4XTsgvddVdYx9Mhv6lJ+/Zw\n1VXQr18oNpdcAjumfFmbNWvArEmmddx9bv4Mm78AvwWWEwpL4R4otxPO9JgKdDSzqjNXDiRMn+QI\nUysPsfEozJbcSFiMOYrwA/3vbLyTzyWE0jSM8IPzdtjk4uBnCEuZywjrPgq/x/vyi1svLThmEtUL\naGHT+5m8bmaDgBfzn3kN0LfGsasIC0v3IxScScCJBdMhNd93i0vO3X26mZ2R/5zLgbcIZ9fc4+5b\nWnf0I8Lf2bOEPYWn51+Pu79rZj8j/Hn3IazpucLdqxYDryTsg9OHUNymAT9290/yz2/091+bxbG2\n+WktaSvyp4Q9w8svwzbbpB1HRBrTrbfCa6+N9VyuLj/4pZUws1uA0929vlNFqdDIiUDYGSH8hqVy\nItK6lJVBHNdqKF1aPgsbzY0ljGgcSxg9ui3NTPWhciIQhkTD2Qtdu6YcRUQaVWlpTJg+kLahO2EK\nqTMwl7AJW4u71IDKiUDhyImItC6lpQkqJ22Gu19I2DulRdPqRwGVE5HWq7wcVE6khVE5Eai60srq\n1Vs5TERanPJyI5xGK9JiqJwIQCmZzFpKNrvTtIi0RBUVsHZtlurrJou0CCongrs7mcxCPv007Sgi\n0piWfr6dxNw0Y4jUlcqJBHE8u+AfMhFpDap/4VA5kRZF5UQC93ksWZLb+oEi0mJUl5MFacYQqSuV\nE6kyXyMnIq3M0qWQzZY09PL1Is1N5USqzGf16iwVFWnnEJHGEkZONKUjLY7KiVQJl/9etizlGCLS\naBYvTsjlPtn6gSLFReVEqkwFqi8VLyIt34IFMeGS9SItisqJVFlIFJUye3baOUSkMZSVwfLl7YAP\n044iUlcqJwLk9zqBySonIq1E9f+XJ6UZQ6Q+VE6kWhxPZObMyrRjiEgj+OQTMMtRNWUr0oKonEih\nD1mwIEul+olIizdrFkTRDHfXKXjS4qicSKHJJImxQPs1ibR4M2bkyOX+k3YMkfpQOZFCk4EwHCwi\nLVeSwKxZGbTeRFoolRP5nLuvJJudx5QpaUcRkYZYtAgqKlROpMVSOZEN5XKvM3GirrEj0pJ9+CGA\nA+NSTiJSLyonUtNbzJ6dpbw87RwiUl+TJkE2+5G7f5Z2FJH6UDmRmsbgDh99lHYOEamvCRMqyeVe\nSzuGSH2pnEhN04milUyenHYOEamPZctg6dJ2wJtpRxGpL5UT2YC7O0kymokT47SziEg9TPp8Dexb\nacYQaQiVE9mY+1tMmWJUaO8mkRYnrDf5xN2Xph1FpL5UTmRTRpPLZbTuRKQF+s9/cuRyr6YdQ6Qh\nVE5kUyYSRSWM01mIIi3KwoWwcGEWeCXtKCINoXIiG3F3J46H8vbb2u9EpCUZOxbMKoGRaUcRaQiV\nE9mcYSxYkGXx4rRziEhtjRkTY/aqu5elHUWkIVROZHNGYJbj7bfTziEitbF6NUyenCFJXkw7ikhD\nqZzIJrn7asxe5a23dEqxSEswbhwkiQFD044i0lAqJ7J5SfICkydnWLUq7SQisjVvv+1E0fvuvijt\nKCINpXIiW/ISSQJvaqNJkaJWUQHvvpsQxy+kHUWkMaicyGa5+xIymVG88oqmdkSK2ZgxsH59BDyT\ndhSRxqByIluWJIP4+OOIhQvTTiIim/PKKwlR9J67T087ikhjUDmRrXmRTGYto0alnUNENmX5chg/\n3ojjR9OOItJYVE5ki9x9DUnyDMOH53BPO46I1DRyJEAOTelIK6JyIrXxOJ9+mmXy5LRziEghdxg+\nPAe84O6fpR1HpLGonEhtvEEULWbEiLRziEih6dNh3rws7oPSjiLSmFROZKvcPSGOH2bEiIQy7Yot\nUjSGD4coKkHX0pFWRuVEamsAlZUJw4alnUNEAEpLYfjwhDi+z911kU5pVVROpFbcfRHuzzBkSI5Y\n256IpO7ll6GyMgHuTzuKSGNTOZG66MeyZVneeSftHCJtWy4HQ4bkcH/C3T9NO45IY1M5kVpz9/FE\n0TiefVZDJyJpevNNWLEiC/RLO4pIU1A5kbqJ47v48MOImTPTTiLSNrnDM8/EZDKj3f2DtOOINAWV\nE6mrfxBFi3nuOe3IJpKGjz6C6dMjkuTOtKOINBWVE6kTd88Rx39h5Eh0vR2RFDz9dEIUzQJ06py0\nWionUh8PksmU8PjjGj0RaU7Tp8PYsRni+GZ3T9KOI9JUVE6kztx9LXHch5EjYd68tOOItB2DBiVE\n0WzgqbSjiDQllROpr4fJZD7V6IlIM5k2Dd55J0Mc36BN16S1UzmRenH3dcTxzbz6qjFnTtpxRFq/\nBx+MiaIZwNNpRxFpaion0hCPEkWLeOwxjZ6INKX334cJEyLi+Gp31z5D0uqpnEi9uXsFcXwjo0cb\nU6akHUekdXKHBx6IiaIJwAtpxxFpDion0lADiaLJ3HNPTKKTB0Qa3ahRMG1aRBxf7u4apZQ2QeVE\nGsTdY+L4QqZNixgxIu04Iq1LeTncd18Os+fc/fW044g0F5UTaTB3fxOzZxkwIEd5edpxRFqPwYOh\ntDSH+2VpRxFpTion0jjcL6e0NGbw4LSTiLQOc+fCc885SXKLu89PO45Ic1I5kUbh7vNJkj8xZIgz\nX/+OijSIO9xzT4LZPEDX0JE2x7S+ShqLmXUiiqZz0EFf5o47MpilHan5lJTAgw/CuHGwfj3sthtc\ndRXsu2/1MY8+CsOGQVkZ9OgBf/hDOG5zHnss3Ap16waDBlXfv+8+GDECOnWCX/8avvvd6udefz0s\npvzznxvlW5Rm9OabcOONACe5+/+lHUekuWXTDiCth7uvNbPzmTBhGMOGwQ9+kHak5lFWBhdfDAcf\nDH37wg47hIsibrdd9TFPPw0vvghXXw277gqPPAJXXhmKRrt2m3/vvfaCO+8Mv0kDRFH1c2PHhgJy\nxx0wf3747J49YfvtQ6ZHH4W77mqSb1maUGkp9OuXI5MZ4XGsYiJtkqZ1pFG5+3BgIP37x3z6adpx\nmsdTT8FOO8EVV8B++8Euu8A3vhFKSJXnn4ezzoIjjwyF45prwmjLmDFbfu8ogi9+ETp3Drftt69+\nbt48OOgg2GcfOO44+MIXYPHi8NwDD8DJJ0PXro3//UrTuvdeZ9WqdSTJb9KOIpIWlRNpCpdSWVnC\nHXcktIVpw7FjQym56SY45RQ4/3z4v4JfeBcvhhUrwshKlW22ge7d4eOPt/zeCxbAaafBz34WpmeW\nLq1+7r/+K1yltqwsXHeloiJME334IcycGbJIyzJmDIwcaSTJhe6+IO04ImnRtI40Onf/zMx+yfjx\nbWN6Z/Fi+Oc/4fTT4ec/hylT4N57oX17+N73QjExCyMfhTp3Ds9tTvfuYd1Kt26wfHmYArrkkjBd\n06kTHHpoWGNywQXQoUMYjenYEfr1C9NHL74IL7wQRl4uvRT23LMp/xSkoVatgr/+NYfZK7jrtDdp\n01ROpEm4+3AzG0j//mdzyCERO++cdqSm4x5GTn71q3D/K1+BOXPgpZdCOdnS67akZ8/qr/faC/bf\nH844A0aPhu9/Pzz+i1+EW5VBg8KUUhTBE0/AwIFhZOe228JUjxSvfv2csrJy3M/TTrDS1mlaR5pS\nmN75y1+SVr21fZcuYXSj0B57VE/BdOkSisjKlRse89ln4bna2nZb2H33sNh2U+bNg9deCyVp4sSw\nHmX77eHYY2HGDFi7tvafJc3r9ddh9GgjSS5w9yVpxxFJm8qJNBl3/4w4/jnvv288+WTacZpOjx5s\ntLfL/Pl8Plq0666hhEyYUP18eXmY/jnggNp/ztq1sGgR7Ljjpp+/807o3TtM7SQJ5HLh8VwuTCu1\n5oLYki1ZAnfcEWP2PPBM2nFEioHKiTQpdx8F3MLAgc7776cdp2mcdlooGk8+GUY1Ro0KC2J79ao+\n5tRTwzTL2LEwa1aYZunaFY46qvqYSy8N60SqDBgAH3wQfnhNngzXXx+ma447buMMQ4eGNSxHHBHu\n9+gB778fFtw+91wYydlmm6b5/qX+Kirg+utj1q9fpOkckWrahE2anJlFZDKj2G67o3n00WydpjJa\ninffhYceCuVkl13C4tgTT9zwmEGD4OWXw9k1Bx4Iv//9hpuwnXkmHH989RqSW26BSZNg9eqwqLVH\nj7DRWuEpyhCmiy68EPr333CaaPBgGDIklJZrrgnrYqS43H03DB1aifvh7j5h6y8QaRtUTqRZmNnO\nRNFkevTozJ13RhtsJibSFlXv3vsbd38w7TgixUTTOtIs3P1T4vg0Pvggs9GW7CJtzZw58Ne/xpg9\nATyUdhyRYqNyIs3G3UcD1zN4MLz9dtpxRNKxdi1cf32OOJ6B+wVaZyKyMU3rSLMyswxmz9Ou3Q/p\n3z9in33SjiTSfOIYrrsu4d//XkeSfMPdp6YdSaQYaeREmpW7J7j/nDj+kKuuyrFsWdqRRJqHe1i0\n/N57kCSnqJiIbJ7KiTQ7dy8njk9k9eplXH11rM3BpE0YMiScKu7e293/lXYckWKmciKpcPfFxPEJ\nzJmznltuSYjjtCOJNJ233oK//x2gr87MEdk6lRNJjbtPIkl+zDvvGAMGpB1HpGlMmQK33JJgNgS4\nJu04Ii2Byomkyt2HAxcxZAiteot7aZvmz4err86RJONxP9vddQ0BkVpQOZHUufvfgZt4+GF44YW0\n44g0joUL4fe/z1FePos4PsndtbhKpJZUTqRY9AHu4m9/g1deSTuLSMMsWRKKyapV84njY91dp6WJ\n1IH2OZGiYWYGDADO57rr4LvfTTuSSN0tWwYXX5yjpGQRcXyUuy9IO5JIS6OREyka+Z0yewMDufVW\n5/XX044kUjfLl4cRk5KSpcTxN1VMROpHIydSdMwswmwQ8DOuucb43vfSjiSydStWhGKyaNHy/IjJ\nJ2lHEmmpVE6kKJlZBDwAnMtFF8Gpp6YdSWTzFi6Eyy7LUVKygjg+xt2npx1JpCVTOZGilV+Dcjtw\nJWedBb/8JZilHUtkQzNmwBVX5Cgrm0scf8fd56YdSaSlUzmRomdmVwJ/4Uc/gt/9DqIo7UgiwcSJ\ncM01MZWVk4jj43VWjkjjUDmRFsHMzgUe4thj4dprjXbt0o4kbd2bb0KfPgnuo0mSk929NO1IIq2F\nyom0GGbWC7NnOfBAo0+fiB12SDuStFUvvgh/+5sDz+V3fl2fdiSR1kTlRFoUMzuGKHqJHXfcjttu\ny7L33mlHkrakshLuvReGDgX4G3Cpu+uqlSKNTOVEWhwz24Moepko+m/++McMxxyTdiRpC1asgBtu\niPn4Y8f9And/JO1IIq2Vyom0SGa2DWaP4X4q55wDZ50FGe0pKE1k2jS49tocq1Z9Rhz/yN3fSTuS\nSGumciItVv5U4z8CfTjmGOfqq40vfCHtWNLajBwJffsmuL9PHP+vuy9MO5JIa6dyIi2emZ1MJvMU\nO+3UjptuyrLffmlHktagogLuvz8sfoVBQG93X5duKJG2QeVEWgUz+wpR9CzuX+O884zTT9c0j9Tf\n3Llw440x8+YluF8CDHD9YynSbFROpNUws/bALcCVHHxwzLXXRuy4Y9qxpCVxD2fi9O+fkCSfEMen\nuvuHaccSaWtUTqTVMbPvEUVP0alTZ669NuKII9KOJC3BihVhbcm4cRnCdZ0udfc1accSaYtUTqRV\nMrOuZDKPkSTf54QTnN69je23TzuWFCN3eOMNuOuuHGvWrCaOz3H3oWnHEmnLVE6k1cqfzXMemcyd\nbLttRy65JMu3v62LB0q1JUvg7rsT3nsvQybzEknyG3f/NO1YIm2dyom0ema2K2b9cT+Fww5L+MMf\nMuy8c9qxJE1xDM8/D488khDHy4jjC9z9xbRjiUigciJthpn9L1H0AFHUlfPOy9Crl65w3BZNmwZ9\n+8bMmpUB7gWud/fVaccSkWoqJ9KmmNn2wG1Ab/bYI+bCC7McemjasaQ5rFwJjz8OL73kZDKTieNz\n3f3faccSkY2pnEibZGY9iaJ7iOPDOfTQhN69M+y1V9qxpCmsWwdDhsCTT8ZUVKwlSW4C7nH3XNrR\nRGTTVE6kzcovmO1FFN1JkuzBiScav/oVdOmSdjRpDHEM//oXPPxwjs8+A/f+wJ/cfXna0URky1RO\npM3Lb97Wmyi6mSjajp/+NMOpp8J226UdTerDHcaNg/vvzzFvXhazZ3C/1t1npR1NRGpH5UQkz8w6\nA9dhdjEdOkT06hVx2mnQuXPa0aQ2kgTGjIHBg2NmzozIZMaQJJe5+3tpRxORulE5EanBzHYGLiWT\nuZgo6sAPf5jhjDOga9e0o8mm5HIwahQ88USOhQuzZDJvkSR/AkbqejgiLZPKichmmFkX4GIymUsx\n247jjzd+8hPo1i3taAJhoeuwYfD00zlKSrKYDcX9Vnd/N+1oItIwKiciW2Fm2wEXEEVXEcc78tWv\nxpx8csTRR0O7dmnHa3vmzg0X5xs+PGbtWgP+H+63ufvktKOJSONQORGpJTPrAJxCFF1EHB/JDjvk\n+OEPs5x0EtpxtolVVITr3/zznzGTJ0dE0Uri+CFggLvPTjueiDQulRORejCzHkBvMplzcO/EYYc5\nJ5yQ4YgjoH37tOO1Du4wZw688goMG5ajrKxqPcl9wIvuvj7tiCLSNFRORBrAzLYFziSKfkscH0TH\njjHHHhvxne/A17+u7fHrY84cGD0aXn01x4IFWaJoFXH8CPCgu09LOZ2INAOVE5FGYmb7E4rK2cTx\nHmy7bY5vfSvLt74FBx2kEZUtmTs3FJLXXgt7k2Qy5STJP4BnCWfdaJREpA1RORFpZPmdZw8GfkwU\n/ZQ43oP27WO+/nWjZ88MhxwCu+8OZmlHTU9ZGUycCOPHw7hxOZYsyZLJrCkoJCNUSETaLpUTkSaU\nLypfBY4nk/k+7kfjnuVLX8px+OFZDjkEvvrV1r/RWy4HU6eGMvLeezHTpmVIEiObnUsuNxx4BfiX\nu69LO6qIpE/lRKQZmdk2wLHA8WSzPyCX2xuArl0rOeCAdnTvDt27wz77QMeOaUatvySBhQtDGZk6\nFT7+OMfMmRlyuQyZTCnuI3AfSZiu0ZbyIrIRlRORFJlZN+AIoCdRdCRJcjDu7clknD33jNl33yy7\n787nt912K669VUpLYd68cJs7F6ZPT5g61Vm7NqwEDiMjbwP/Bt4Bxrt7nGZkESl+KiciRcTM2gE9\ngJ7AYWSzXyNJ9iVJtskfADvtVMmee2b58peNHXcMU0JduoRb587hls02LEiShOKxciV89ln1f1es\ngEWLYMGCHAsXQnl51Qc52exi4ng84Vo27xGKyMqGBRGRtkjlRBqFmT0NuLufmXaW1ia/bmVnYL/P\nb2b7E0X74L4TcbzDRi/aZpsc7ds7HTo47dsbHTpAhw4ZOnQwOnTI4A4VFQkVFf75rbISKiuddeuM\n0tIsSWI1glQSRStw/4Q4ng58AswEpgLT3X1NU/9ZiEjb0MBfr6RYmFkCOLCpU0AcuNnd+9Tifd4B\nXnf3axs5otRT/uJ1S/K3N2o+b2btga7ALvnbzpSX70R5eSegE9Cx4FZ1PwEqgPUFt6r75cDS/O3T\nz792X+2VlfptRkSanMpJ67FLwddnADcD+1JdVsqaPdEWmFmG8HNXP+wayN0rgIX5m4hIi5dJO4A0\nDndfWnUDVoWHfFnB42sAzOy7ZjbezNaZ2UIz65OfNqiamjkMuMrMEjOLzWwnM2tvZgPNbI6ZrTGz\nKWbWuy75zOw3ZrbYzE4xs6nAOsJv+5hZbzObamZrzewjM/t1wesmmNkNNd7ry2aWM7Nv5O93NLN+\n+e+n1MzGmNmRm/jsH+Q/Z7WZDTWzHQuOecfMbq3xOcPN7O8F97f4OSIi0jhUTtoQM9sDGAqMJuy9\ncRFwIXBF/pDfABOA/oSRmF3zZSdLWF9wMtAduBW4w8xOqmOELwK/A84CDgRWmtm5wFXA5cD+wA1A\nXzM7Lf+aJ4Gf1XifM4GZ7v6f/P2HgIOAU/Lf18vAiPyZMIWffSHwE8KpvPsBt9cxf20+R0REGkjT\nOm3LxcBUd788f3+6me0JXAv0dffVZlYJlOdLCQD5UZc/FbzPYDP7JnA64Qd0bbUHfu3uM6seMLOb\ngIvcvep95prZ14ALgOeAp4HbzewbBWXkp8AT+dfvQygcu7j7ivzzt5vZD4CzC3K3B85198X5192f\n//OolTp8joiINJDKSdvSHRhb47G3gS5m9iV3L9ncC83s94QRj26EBZXtCftW1EVZjWLSGdgNeMI2\n3Mo9Iiz+xN0XmdmbhNGT/+SvX/M1oGpk5cD88XNswzdpTziTpMqKqmKStxjYqQ7Za/s5IiLSQCon\nbYsRztyp+RibeLz6ALNzCCMDlwDjgVLgesLUSF2U17i/Xf6/ZwMf1HguV/D1k8AtZnYZoaS8V7Cz\n6LaEM0wOYuMzlUoLvq6s8Zyz4bRmsonXF+52VtvPERGRBlI5aVs+Br5T47GjgOXuvjx/v4IwQlDo\nSMLpxY9UPWBmX2mEPPOBEmBvd39hC8cNIayD+Q5hSufugucmAB2ALgXTPvWxDNi16k5+M7T/BqY3\n8ueIiMhWaEFs23IvsJ+Z3Wlm+5rZqcB1wF8LjpkDHGFmuxeczTIj/9hxZraPmd1OmOZokPxpxDcD\nN+TP2NnHzA40s3PN7MKC41YDw/I5uxGuWlv13GTgH8DTZvYjM9vTzA4zs+vMrGYR25LXgJPN7H/y\nU0cPEvYEaezPERGRrVA5aUPcfS5wEvBNwjTKPYTCckfBYbcT1lFMBZaa2U75Y4YBzxPWrHQknLnS\nGJnuI5w1dD4wiVASzgRm1zj0ScIZMqPcfVmN584kFJZ++dxDCNMvC+oQZQDwDPBUPsNENl5T0xif\nIyIiW6Ht60VERKSoaOREREREiorKiYiIiBQVlRMREREpKionIiIiUlRUTkRERKSoqJyIiIhIUVE5\nERERkaKiciIiIiJFReVEREREiorKiYiIiBQVlRMREREpKionIiIiUlRUTkRERKSoqJyIiIhIUVE5\nERERkaKiciIiIiJFReVEREREiorKiYiIiBQVlRMREREpKionIiIiUlRUTkRERKSoqJyIiIhIUVE5\nERERkaKiciIiIiJFReVEREREiorKiYiIiBQVlRMREREpKionIiIiUlRUTkRERKSoqJyIiIhIUfn/\njS5nv52CUx8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fdbad631590>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "def display_pie(pie_values, pie_labels, colors=None,title=''):\n",
    "    plt.axis(\"equal\")\n",
    "    plt.pie(pie_values, labels=pie_labels, colors=colors, autopct=\"%1.1f%%\")\n",
    "    plt.title(title)\n",
    "    plt.show()\n",
    "                                                           \n",
    "display_pie( [kpnv[1] for kpnv in all_kpis], [kpnv[0] for kpnv in all_kpis],title='KPIs: Revenue - Oil Cost - Production Cost')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "##### Production"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeQAAAFeCAYAAABU/2zqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xd4FOXax/HvPbPpCQFCQpUuAlIUBBSxoYCAXexif+29\nosfeyzkqIBZU7L17jl2xN0RQ6aGG0AmhQ9jszvP+MRNYQgJJSDKb3ftzXXslmZ1y7252fzvPPPOM\nGGNQSimllL8svwtQSimllAayUkopFRU0kJVSSqkooIGslFJKRQENZKWUUioKaCArpZRSUUADWSml\nlIoCGshKKaVUFNBAVkoppaKABrKKGSJyp4g4paYtEJFxftUUTURkkIg4ItLb71p2RUQu9mrN8WHb\nb4rIjNrerlIayKrCRKS1iDwhIrNEZKN3m+ZN6+p3fYDxbpGcMqbVKBGZ74XJzm5hETmrkuu9TUSG\n7mZ5u3wuROSiUrVuFpGZIvKYiDTaze1XVFmvZbURkT1E5A4R6VzOtp0ypitVowJ+F6DqBhE5CngT\nKAZeA/7G/dDqCJwAXCwibYwx+f5VWaa9qP0P16uA9Ii/hwKnAlcDqyKm/1LJ9d4OPAt8slvVVYwB\nRgCLgRTgYOBKYJCIdDfGFNdCDTWpJXAHMAOYXuq+MwGp9YpU3NNAVrskIm2BN4D5wOHGmBWl7r8R\nuIwo3KvwIziMMR9H/i0iTXED+SNjzMLarmc3fGKMKQmr50VkPXAJMAT4qKwFRCTVGLOptgrcDeUG\nrjEmXJuFKFVCm6xVRdwEpALnlg5jAON6whizuGSaiHQVkRdEZK7X5LlURJ4XkYaRy4pIuog87jXz\nFonIchH5UkT2KTXfSSIyUUQ2ichKEXlFRJrtqvDSx5BF5GyvGbaviDwqIitEZIOIvC8iWWUsP1hE\nfvDmWSci/yunmbPKRCRDREaJyCLvOZguIldG3J/kHRu3cFsiSpqSn/Tubysiz4hIbsTz84aItKjO\nOoHxuEHWxttuSS0HiMhYEVkJzI6ou5eIfOU9b+tE5AsR6VnG4+/uPcebRCTP+4InpeZJ8rZ1YxnL\nLyt5LiKmNfSe0zzvOc0TkXEiUk9EBgE/4LYCvBlxCOFkb9kdjiHv6jUqVePDIjJM3MM5RSLyj4j0\nr9QzreKS7iGrihgKzDHGTKzEMgNwP7jHAcuAvYGLgM7AARHzPYPb5D0at/kwCzgQ6AT8BSAi53jr\n+R23GbUxbvNvXxHZ1xizbid1lHcccjRQCNwJtAauAZ4ATiuZQUSGAy8CnwM34n4puQT40dvubu/t\niogFfAbsj/tcTMV9vh8XkcbGmH8BQdxm1JdxQ/EFb/Fc7+cBwL7Aq7hNzO2AS4EeItKlGlsJ2ns/\nS5rdS57bZ4EluE3qyd7j2hf4DigA7vPmuwT4QUT6GmP+9uZrAXyLeyjkXu+xXgxsqERd273GIlIP\n93BAa+A53MMrOcBxQBPv73uA23Bf89+8RX+NWJ+JWF9FXqNIR+C2iDwJbAKuBd4TkZbGmPWVeFwq\n3hhj9Ka3cm9ABm5T9Htl3JeJG6Alt+SI+5LKmP8UIAwcGDFtNTBqJ9sP4Ab6X0BixPQhXl13REy7\nAwiXWn4+MC7i77O95T4vNd9/cMMgw/s7DTewnyo1X7ZX89OVeA6v8x53y3KeEwe4utT0j7x6mkdM\nKwaeLGMdZT3XB3vrPTFi2iCvjt67qPcib76+3uvaAjjDez7WAI0i5nOAL8tYx2e4odosYloLYCPw\nWcS0p4AQ0CViWmNgvVdDTslj9LZ1YxnbWhr5vAAPecsO3MljPNBb38ll3PcGML2yr1FEjRtLvW69\nvOnn1cR7VG+xc9Mma7Ur9byfZe2xfAesjLhdWnKHMWZLye9eU14W7h6uAD0i1rEG6O0dZy3Lfrh7\nN08aY4IR6/8UmIm7p1JZBhhbatqPgA208v4eiPuF400RySq5ecv+DhxWhe2WZTCwGXi61PRHcb+M\nDNrVCko91wneYYHpuHtnPcpdcOcE+An3dV2Iu3deABxrjCmI3DzuXuO2BUUSgf7A28aYJRF1LgLe\nBvqLSJI3eTDwvTFmasR8y4G3qlg3uC0uvxtjvtyNdUSq7Gv0iYk4fGOM+QPYArStpnpUjNIma7Ur\nJU1s6WXcdyHuHnRj3ObSrUSkAW5z8Cm4gVrC4AZdiRtxm4XzReRP4FPgZWPMfO/+Vt4yuexoJu6e\nTlWU7g2+2vvZwPvZHjeUvi1jWQOsreJ2S2sF5BtjikpNnxFx/06JSCrwL9y9/6ZsO/5a+rmuDANc\nAOTh7sEuNcbMLmfeBaX+bgokUPZrNgP3c6eZiCwA9sA9JFDarMqXvFUbyn7dqqqyr1FZZxqsYdv/\nllJl0kBWO2WMWSciS4EuZdz3B4CItGLHXqvv4B5zexj3mN0G3E5JXxDRmdAY846I/AAcj7tXej1w\nk4gcb4z5ooz1VpeyetJKxPYs3FA6E1hexryhaqqjOh7fWOAk3D22CcA63NrfZ/c6bv5utvWy3pnN\npf6u7GMq6zh/6XXs7Jxke+tCIjXx/1LZdZbXS1tPpVI7pYGsKuIT4HwR2c9UoGOXiNTHbbK8zRhz\nX8T09mXN7zVRPg08Le7AE5Nx9/i+wN37Etzzib8rteheuHtw1SXyQ3+ut92Vxpjx1biN0hYAvUQk\nKbLpGbdTG2z/+MoLpROAscaYm0smiEg62w431LYluMe79yrjvk64X2aWGGOMiOQDHcqYb7tljTFB\nEdkE1I+c7rUONIqYz4jIfMr4AllKZQYdWUDFXyOlqkyPIauKeBh3L2iclD2UYen/o3A506+hVO9V\nr0fsVt7xySW4HWQAJgIrcE/3SYhYdjDuB+L/KvdQKuwL3D3NW0Rkhy+uUn0jVn2KO/DGxaWmX4Mb\nXJHNuRspFUieMGU/177wjvV/AwyTiFPTRKQ57p78NxHB9ilwiESM9Ob1JzipjFXPw+2sFumyMuZ7\nD+jjnd5Uno3ez7Kez9Iq8xopVWW6h6x2yRgzR0ROB14HZolIyUhdJeekno4bCou8+dd7zdA3eh18\nFuM2R7dh+2a7DGCRiLzLtmbtAbgdua711hUSkZtwT3v6QUTewD115UrcD+jHq/CQyms63DrdewyX\n4HZmmiQib+J2cGqJ25HsJ6+G3fUu8DPwbxHpwLZTagYDD0R2igL+BAZ7578uxz0V7U/cFowLRGQz\n7nHbfrjH1tfs7DHWsFtwn6NfROQpb7slgTYiYr4HcE8R+kZERrLttKe5QLdS63wO91SjN3GPEffE\nDejSx/MfwD0E8rGIPI/bQ78R7mlPZxpjcnGPUW8ELheRYtwOcD97Hc9Kq8xrpFTV+d3NW29154Yb\nqE+w7cNsAzDNm9a11LxNcT/IVuGeLvMGbuevMG5TNrgdfx4EJuGGxzrv9wvL2PYw3L3lTbjB+BLQ\ntNQ8dwChUtPmAc9H/H22V0OPUvMd4k0/uNT0g3H3kAq9x5wLPA/sW4nnrdzTnrz703G/WCwCinA7\nC11RxnydcQe02Oit70lvegPcjnErvOfxY++1WgKMiVi+sqc9dd6d+XADs6SlYR3unmSPMubrDnzv\nvbZ5wA24obz1tCdvPgt4xHuc63FPO9qj9OP05s0CxnjP6Wbc09/GAvUi5jne+//d4m3rZG/6G8C0\nyr5GuK06YeChMh7jDjXqTW+lb2JMrY67r5RSSqky6DFkpZRSKgpoICullFJRQANZKaWUigIayEop\npVQU0EBWSimlooAGslJKKRUFNJCVUkqpKKCBrJRSSkUBDWSllFIqCmggK6WUUlFAA1kppZSKAhrI\nSimlVBTQQFZKKaWigAayUkopFQU0kJVSSqkooIGslFJKRQENZKWUUioKaCArpZRSUUADWSmllIoC\nGshKKaVUFNBAVkoppaKABrJSSikVBTSQlVJKqSiggayUUkpFAQ1kpZRSKgpoICullFJRQANZKaWU\nigIayEoppVQU0EBWSimlooAGslJKKRUFNJCVUkqpKKCBrJRSSkUBDWSllFIqCmggK6WUUlFAA1kp\npZSKAgG/C1BKbU9EBGgANAOaez+zgGQgpdTPZCBFIMlAMbC5jNsm7+c6YAmw2LsVGmNM7T0ypdTO\naCAr5QMRSQc6A128n3sEbFqK0EKExsaQEDl/SpKEkxLESUoQkhKFpAQhJVEkOUkkOdGyEhNEwmHD\npiLjbA46zuYtxhQFDUVBw5ZiI1uCRoLFxo5MX0sIJibIcschP+ywEDekZwPTgenGmFW19HQopQDR\nL8hK1RwRCQCdgK5AF4EuAZt9i8O0ABCgaSO7uGVOgt24oW3lNLBp3CBA4wY2OQ1scuoHyK5vk5gg\nu11LcchQsDbM8sIQywrDLF/t/r58dZilq0LhJatCzrJV4YBjEICAzYZQmPnAF8Bk75ZrjAnvdjFK\nqR1oICtVjUQkEegFHGwJhyIc5DikAGRlWsWdWibae7VMtPZskUiHPRJo1yyBlKTo6coRLDYsWFbM\nnMXF3PPSKlZtdiCFYtZ5e+xCEcLfOHwFfA38aowJ+lq0UjFCA1mp3SAiycD+wCGWxaEYDnAMSSmJ\nEu65V5L06Zxs9dgzmQ4tE8hMs/0ut8Icx9D9vIVsaW3gDNwj0MuApcBiDHMJU0QAoQj4FsOXuAE9\nTY9LK1U1egxZqUoSkUxgqAgnWsIQx5CcniKh3p2S7d6dkqVXx2Q6tUq0A/buNzP7ZdHKEFuKDezh\nTUgB2ng3EBwCLAPmkcxcBpHHIBwsLFaKyOe44fy1MWaJLw9AqTpIA1mpChCRHOBYy2KYCP2NIdC5\nVWLoyD6pgUP2SaVDi4SAZdXdAC5tVr7XCt2+nBks3L7fzYB+WBQDC4G5ZDOH01jBcACxZRIOLwJv\nGmNW1nTdStVlGshKlUNEGgGn2xanAAeIwH4dkpyBvdPsAful0jQrELPvn9z8Ymwbwk0quEAC0M67\nDSTABmAeMJV9mM1IDI+JJZ9heAn4nzGmqIZKV6rOitkPFKWqQkRsYKAI54twrCXYB3VLMQN7pUr/\nnqk0zLDrzoHg3TBzYRCTiMGiarv96UA3oBsWG4Gp2PzFkSzlKIQNIvI68Arwsx5zVsqlnbqUAkSk\nGXBBwOaSUJgm7ZsnhE7pnxE4+sA0GmbERQZv5/BrFpEvIbikmle8EvgH+IsQ6wlgkY/DC8CLxpj5\n1bw1peoUDWQVt7wRsQ4V4XIMxyUkCMf1S7NO6Z9BlzaJuHfHn6KgQ/fzFmK6AifU0EYc3GPOfwNT\nCVOMBbwNPGiM+auGtqpUVNMmaxV3RMQCjrItbg879GzdJBA6c2A967h+6WSkRs85wX6Zu7gYY4BW\nNbgRC2jt3QZj8xfwEyeyjlPEkq8w3A98r83ZKp5oIKu44R0fPilgc1soTOfu7ZPClx1fn35dkwPx\nujdclln5xe4vHWppg4lAb6AnAaYBP3IYKxmAxZ8ich/wkTHGqaVqlPKNBrKKed7oWcMDNreGwrTe\nv3Oyc8lx9enVMTn+Dg5XQG5+kEAChDJqecM2bkewrgSYA/zIPizkfSzmisj9wKs6KpiKZRrIKmZ5\n40hfYNvcHg7T9LB9U80lx2XSpU0UjVUZhWbmBQkl+ViAAHsCe2KTD/xEG2bxPBb3i8idwLM6nraK\nRdqpS8UkERlo24x0wnQc2jfNXHpsprRvkeh3WXVCn4sWsrqhA+f5XUmElcCPGP5BsJiOw+XGmG/9\nLkup6qR7yCqmiEhHS3gUGNy9XVL41rMa0qVNkh4grqDC9WFWb3Bgb78rKSUbOAGhD/Ape7GY8WLJ\nBxiu09OlVKzQQFYxQUSygDtEuLRxQ9vccmZDBvZKtbWzVuXkLvQO0bb1t45yNQcuwGYK8CVHs5Gj\nROQR4AFjzAafq1Nqt2ggqzrNO058qW1xT2JA0i47ob599qAMkhL1MHFV5OYXIwImWgMZ3GPM3YCO\nBPgJ+JkRGP5PRK7H7filPbJVnaSBrOosEelsW7zsOPQ88ZB0rj6pAY0yteP07piVH8ROgFCC35VU\nQCLQH+iBxZc0YjovYXGliFxqjJngd3lKVZYGsqpzvL3iG0S4u3l2gEcuacS+eyb7XVZMmJEXJJTm\ndxWVVB84GWEB8CndWcFv3vnLdxtjiv0tTqmK017Wqk4RkS62xcuOYZ8LhtaTK06sT7I2T1cLxzF0\nP28hW1obOMPvaqooDPwEfIdB+AuH04wxs/wuS6mK0E8yVSeISIKI/EuEyS0bB7q+fWdTueG0hhrG\n1WjRyhBbig3s4Xclu8EGDgEuQMikK8LfInKpaO8+VQfop5mKeiKyl23xhyXcc+HRmYGP728W6N7e\nz5ErYtOsfK+HdXt/66gWzYFLCNCTJGAMwmciUtGrOyvlCz2GrKKaiJxsWbzYIjuQ8Ojl2dK1rQZx\nTcnNL8a2IBwrsZUIHIU7JveHHE4RM0TkPGPMBz5XplSZdA9ZRSURSRSRkcBbR/ZOTf7w/mYBDeOa\nNWthEJNI7H0qdAAuI8Ce1APeF5FxIlLbI3UrtUu6h6yijojsYVu8B+x3y/CGnDkgQw8B1oLpC4I4\nmX5XUUPSgFOxmAx8ytk4HCQig40xc/wuTakSsfZdWNVxIjLQtvg7q56975t3NJXhA+uhYVzztgQd\nFq0MQY7fldQgAXoAF2ORSWuEiSJyiN9lKVVCA1lFBRGxROR24PMD9k7O/O+D2nGrNs1dUoxjgNZ+\nV1ILGgH/R4BWZABfi8j5fpekFGggqyggIkkivArcdeWJ9eW5GxtbDTJ0xK3aNHOhN37Gnv7WUWtS\ngeFY9CQAPCci/xYR/adTvtJjyMpXIlLftvhQhIP+fWk2Q/ava8NExYbc/CCBAITq+V1JLbJxe2Fn\nA59zLUInETnVGLPe58pUnNI9ZOUbr/PWL8mJ0u+lm5tYGsb+mZkXJBSPo48KsD9wBkKAQVj8JiKt\n/C5LxScNZOULEelqW/zRqL695zt3N7V7dYrHNIgeMxYGoYHfVfhoT+D/sMmgAxaTRKSv3yWp+KOB\nrGqdiPS3LH5t3yIh6727mwbaN0/0u6S4tnp9mNXrHWjmdyU+ywEuJEBzMhG+FZEj/C6pNngXa1FR\nQANZ1SoRGSbCFwfsnZzyxu1NAzkN9LPAb7klQ2a28beOqJAGnIVNOwIIn9REKLvvAflHRDaJSIGI\nfCkiqSLyrYg8WmreD0RkXMTf80XkVhF5XUQ2iMgiEbm01DKZIvKciKwQkbUi8rWIdIu4/w4RmSwi\n54vIPKCouh+jqhoNZFVrROQEgbeGHpBmj72+sZWeov9+0WBWfjEiQDu/K4kSCcApWBGhPKC6Vu2N\np/068BzQEfdSGO/jHs2uqOuBycA+wIPASBE5POL+d4EsYBDumdeTcE/vqh8xT3vgBOB4bz0qCuju\niaoVInKsCG8P7pMqD1/cSAK2DvYRLWYtDGInYEIJlQqF2FYSym8RYC7/E5GjjDFfVcOam+L27/7A\nGJPvTZsGVGYAnJ+NMY94vz8hIgcC1wDfiEg/YD8gJ+Ja0DeKyPHAMNwvAuA+wuHGmMLdejSqWuku\niqpxInK0CO8O6pVq/fvSbA3jKDM9L2hCqRrGO0jAHW7T3VP+XzXtKf8NfANMFZG3ReSCUnuuFfFr\nGX938n7vBmQAhSKyvuSGO+RLZBtInoZx9NFAVjVKRIaI8P6A/VLt/1ymYRxtHMcwZ1GxkO13JVEq\nQOlQHrg7qzPGOMaYgcCRuHvGVwAzRaQ14LBj03VCRVft/UwHluAGc/eI217AIxHzb6xC+aqGaSCr\nGiMigyzhw/49UqzHLs+WhICGcbRZVBBiS7GBFn5XEsW2D+X/7m4oAxhjfjXG3AXsCxQDxwErcZu0\nAXc4WaBLGYvvX8bfM73fJwFNgLAxZl6pm+4RRzkNZFUjROQwS/j44H1S7JFX5lgaxtEpN96GzKyq\nHfeUq9R8LSK9ReRmEekpInsAJ+KOrj0DGA8MdVuVZC/gKaCs5uwDReR6EdlTRC7DPTb8OIAx5mvc\nJuwPRWSAiLQSkb4icq+I9KhKzar2aKcuVe1EpJNt8XHvTsmBJ67KsRI1jKNWbn4Q24ZwE78rqQNK\nQvlNYC4ficgBxpi/K7mWdcDBwFVAPSAPuNYY84V3PnA34CUgBDyGG9Kl/Qe349adwFrgGi+ISwwB\n7gPG4Q4Mugz4AVheyVpVLRNjzK7nUqqCRKRxwGZi6yYJTd6+s2kgPVUbYaLZVaNW8MXfm3BG+F1J\nHbIFGEeYlRTg0NMYs7i2Ni0i84HHjDGjamubqvbop6WqNiKSYlv8r16q1eS5GxtrGNcB0+YHceLp\nghLVIQk4A5s0srD4XEQy/C5JxQb9xFTVQkQsS3jZtujx7I2NA80a6dGQaLcl6LBoZQga+11JHVQP\nOJMANp0Q3qnF4Se1STOGaSCr6nKfMQx77Ipsq2vbJL9rURUwd0kxjgH02kZV0xg4BRsYCDxcG5s0\nxrTV5urYpYGsdpuIXACMuOmMBgzYTy+hWFfMyvd6WHfwt446rT0wCAGuEZGz/S5H1W0ayGq3uKdU\n8PRpR2Rw7mA9GFmX5OYHCQRwm19V1fUB9sUgPCsiffwuR9VdGsiqykQk27Z4f5/2Sdw6vGFlxuJV\nUWBmXpCQHl3YfQIMRWiGhcV/RSTeL2SpqkgDWVWJiNiWxevpKVajkVdm2zrwR90zY2EQGvpdRYxw\nz1G2SaUBFu+IiO13Saru0UBWVXWLcTji8Suz7SYNtUd1XbN6fZjCdQ7ovlz1yQCGEcChL+4lEpWq\nFA1kVWkichBw12UnZHJglxS/y1FVkJsfdH9p428dMac1cCAA94rIvr7WouocDWRVKSLS0LZ5q0eH\nJOfS4yp71TgVLWblFyMCtPW7khh0GJCDYPGmiOg3VlVhGsiqwkRELOG5lETJefSybFsvpVh35eYH\nsRMwJPpdSQwKAMOwcU+KesDnalQdooGsKmO4Yzj+gYsa2ToSV902PS9IKHWHa++q6pIDDMACrqrq\nlaFU/NFAVhUiIjm2xaij+qaZQb108I+6zHEMcxYVuxf9UzWnD9AaB4tXRUT7s6td0kBWFSLCyNRk\nK/3W4Q11r6qOW1wQoihoYA+/K4lxFnA8FgGyEJ4WPVFf7YIGstolERlqDKfedlZDu2E9Pb2yrsst\nGTKzvb91xIVM4BhsDCcBZ/hdjopuGshqp0Qkw7Z59sCuyc6x/bSpOhbMWhjEtoCmflcSJ7oAXTEI\nT4qIHihQ5dJAVrtyf8CSxvecn2Vpi1tsmJUfxCSi7/7aNAghQCpwl9+lqOilb0lVLhHpC1x23Sn1\nrRbZCX6Xo6rJtAVBHL2gRO1KBw7DBi4Rka5+l6OikwayKpOIBAI2z3dpk+gMH6Sf3rEiWGxYtCLk\nXstX1a7eQAPCCCO1g5cqiwayKs/ZoTAd7zovy7Yt/eyIFXMWB3EM0MrvSuJQADiSAIbDgGP8LkdF\nHw1ktQMRSbMt7h96QJrp2lavzxdLtvaw3tPfOuJWB6AtDhYjRUTfXGo7GsiqLFcD2deeXF93jWNM\nbn6QQAD3dBxV+wQ4EgtDS+Aqv8tR0UUDWW1HRHIsi1uGD6one+RoR65YM2NhkJDul/krB+iNINwp\nIk38LkdFDw1kVdrtyYmSdMmxugsVi2bkBUEHcfTfIUAiicB9fpeioocGstpKRDqIcPHlx9e3G2To\niFyxZs2GMIXrHB0QJBqkAodjA+eKSHe/y1HRQQNZbWUJD+TUt83wgRl+l6JqwNYOXXoN5OjQE8gk\nDNzsdykqOmggKwBEpKtjOOGak+oHkhL13yIWzcoPIoIGcrSwgQMJACeLiL4qSgNZuQSuz6lvh44+\nMN3vUlQNyc0PYidgSPS7ErXVPkAKYeB6v0tR/tNAVojIHghnnH9UvUBCQM90ilXTFwQJpaIvcDRJ\nBPYngHC+iOT4XY7yV8DvAlRUuDo1STjpUD12HKscxzB7UTG0rMWN/gjMAAqABNzrLx8BlFzvaDPw\nLTAXWIfb0akjcBiQXM46w8B4YDaw2puvrbfekn/fEPAxMAt3DOmhbN9M/xOwHhi8m4+vuvQCfsQm\nxJXArX6Xo/yje8hxTkQaWBaXnDWonp2eov8OsWpxQYiioIEWtbjRPKAP8H/AWbhh+grg9S1jPbAB\nGARcChwHzMEN0/IUA8uAQ4GLgVNwA/+NiHn+BJYCF+B2nHov4r7VwGTg8Co/quqXCuyHjXCliOi3\n4jimn8DqEktI0gtIxLatPazb1+JGzwS6A9m4F7M4DlgLLPHuzwFOxh1OsgHQBugP5AJOOetMBoYD\nnYEs3C8YQ7x1rvXmKQD28rbbC9gIbPLu+wQYAFF3HP0AQEjH/fqi4pQGchwTkWTb4rphh6RbjTL1\nvONYlpsfxLaAZj4WUYQ7dGTKLuZJonKfTCXrLWnmbgwsxN2bnoPblJ0K/I17kK5jpaquHZlAN8Di\nRhGJtq8LqpZoIMe34Y5Dw/OG6qhcsW5WfhCTiH/veAN8jnsMu7yuSxuBH3CbmSsqBHwNdMUNcoB9\ncUN5DO7x4pNxj1d/h3vc+BtgFPAqbrN5tDgQwaExcLrfpSh/aCDHsYDNJYfsk+K0bqJjVse6aQuC\nOH4elfgEWAkMK+f+LcDruGF9aAXXGQbext07Hhox3fb+vhq3AXgP4Atgf9zjz7OAS4DmwKeVeAw1\nLRvogIPFDX6XovyhgRynRKRrKMy+ww5N1/+BGBcsNuQvD5W/Z1rTPsHtFX0O23pCR9qC29krGbeT\nVkX+I8PAO7i9s4ezbe+4LPNwvwz0BhbgXnoyAdgbt+NZNNkPC4fOOpxmfNIP4/h1bmaaFTp031S/\n61A1bM7iII4BWvmw8U9w90jPBuqXcX9JGCcAp1GxEzFLwng1bu/tnR2TDuHuBR+NuyftsK3DWJjy\nO4/5pR2QQgj3kak4o4Ech0Qk0bY45/iD0wOJOhBIzNvaw7pDLW/4f8AU4ETcXs0bvFvJaU9bgJe9\nv4/B7ZxVMk9kUI4GZnq/O7jN1EuBE7y/S5YJl1HDd7iPu+Qihy1xz41eDkygds/Lrggb6EYAi7NE\nRMeJiDM4GS/kAAAgAElEQVT6gsenoWGHBicerMNkxoPc/CCBAIRqu+/eRNy90hdLTT8Wd8jIpWw7\nBWqU99N4y1zFtj3qQtywBreJOtf7/elSy5wNtI7Yzgrc8L04Ylpn3GbrF3BPmzqxko+pNnQDfqcR\n7tnSX/hcjapFYozxuwZVy2xL/tuxVeKRH97XTL+QxYHzHlzGT/OK0K5CdYQBRhOikLeMMWf6XY6q\nPdpkHWdEpIkxDDnp0HQN4zgxIy/oDryh6gYB9iWAMExH7oovGsjx50zbhqP6pvldh6oFazeGWbXO\n8XdAEFV5XQFDEu6RchUnNJDjjG1x0qH7pEhmmo7MFQ9yF3o9qNr4W4eqpPpAK8IIZ/tdiqo9Gshx\nREQahx16HdEzVbtWx4mZ+UFEcE+nUXVLd2wMh4pIbV4SRPlIAzm+DBWQQ/bRc4/jRW5+EDsBE3UX\nU1C71hmwMcCpfpeiaocGchwR4Zhu7RPDWXohibgxfUGQUCraIlIXJeOexiVRc+VmVcM0kOOEiKQI\nDBqwnx48jhfGGGYvKoZGfleiqqwdFnCQiOxsPDIVIzSQ40d/x5Dcv4e+r+PF4oIQRUHjXkRB1U3t\nAEMCcKDfpaiap4EcP45u1sgOtWumV3aKF1uHzNzT3zrUbsihZGzrAX6XomqeBnIcEBEJ2Bw/sFda\nQEQPJ8aLWflBbAs9B7kuE6A9ASw9jhwPNJDjQ8dQmJyDu2lzdTyZtTCISUTf5XVdO8Chq4hob4AY\np2/V+HCgCOyz584uGqtizfQFQZx6flehdlvbrb8d7mMVqhZoIMeHvnu2SAilp+jLHS+CxYaFy0Pu\nMUhVt9UDsihGjyPHPP2EjgMJNof26pisF5OII3OXFOMYoJXflahq0Z4ELAaLdgKJaRrIMU5EsovD\ntOmhzdVxJTc/6P6iPaxjQ1vAoRnQ3u9SVM3RQI59fQF6dEj2uw5Vi3LzgwQCuBcpUHXftpaO/X2s\nQtUwDeTY17dRphVq1kgH6IonMxcGCWmjSOxIBjIoxr0wo4pRGsgxzrY4uFfHZD3/OM5MXxCEBn5X\noapVUwII3fwuQ9UcDeQYJiK2Y+i5rx4/jitrN4ZZtc7RAUFiTWMEYV+/y1A1RwM5trUxhoQ9W+i1\n9+JJ7kJvyMzWvpahqlsO4JAjIg39LkXVDA3k2NYRQMevji+z8oOIO+SiiiWNt/6mx5FjlAZybOuU\nnCjhxg21Q1c8yc0PYidg0IaR2JIFWBg0kGOWBnJs69imaYLRDl3xZXpekFAK+qLHGhvIIgTasStW\naSDHsIBNlz1bJOgIXXHEGONedlEvQxCbmpKApR27YpUGcowSETGGTnr8OL4sLghRFDTQwu9KVI3I\nAQx7i4h+dscgfVFjV3bYIaNtcw3keJKb7/Ww1iEzY1NjwJCCjlIekzSQY5fbw7qpBnI8mZUfxLbQ\nc5BjVebW3/QVjkEayLGrFUCLHD2EHE9y84sxieg7O1albf2t8U7mUnWUvm1jV9PUJAknJ+pLHE+m\nzd+Ck+F3FarGpACCQa90HZP00zp2Ncuubzt+F6FqTzBkWLgipPtOscwCUgihr3JM0kCOXU2bNAzo\niCBxZN6SYhwHaOl3JapGuc3WuoccgzSQY5Rt0Sy7vq2vbxzJzQ+6v3Twtw5VwzKw0T3kmKQf2DHK\nsshukKEvbzyZtTBIIADU97sSVaPSsbBo6ncZqvrpJ3aMMoZGDTK0xTqezFwYJKxX2ox9aYBoIMci\nPScmBomIiJBZP12/b8WTGXlBTAO/q1A1Lg1wyPa7DFX99BM7NqUZQ0ADOX6s3RimYK2D7jfFgXTA\nkC4i2h4SY/QTOzYlAiQk6AV/4sXWITPb+FuHqgWpW3/TS4jEGA3k2BQACFgayPEiNz/oXm+xnd+V\nqBq37UCjjosbYzSQY1MCgJ6FHD9y84uxEzFoI2bs2/aprX2AYowGcmzyAln3kOPFtAVbCKWgL3g8\n0ECOWRrIscltstY95LhgjGH2omI9ohgvtn3t0kCOMfqCxibdQ44jSwrCbN5iYAHwiN/VqBq3bYT6\nFB+rUDVAAzk2eXvIGsjxwLIgq0ECq1YXw0ZIa5hG616t0Qbs2LR57Wbm/TrP7zJUDdBAjk3aqSuO\nNM0K8POopjz98TrG/m8dGws3sqFgA0NuGUKHQzsgoskcSxZMXMDjAx8H2OJ3Lap66THk2BQGCOnF\nF+OGZVlcelx9Jo1twVXDMlk5fRFPnfgUjw96nNzvczHG+F2iqibG2fpa6js8xmggx6b1ABs36/s1\n3liWxWXHN2DS2OZceWImK6ct4snjn2Tk4JHk/qDBHAuc8Nb3ddjPOlT100COTW4gF2kgxyvLsrj8\nBDeYLz8+kxVT8nnyuCcZNWQUs3+c7Xd5ajds2bC1pXq9n3Wo6qeBHJu8PWTdG4p3lmVx5TA3mC87\nLpNl/yxkzLFj3GD+SYO5Ltq8bnPJr2v9rENVPw3k2LRFIKx7yKqEZVlcdVIDJo9tziXHZrL0rwWM\nOWYMo4aOMnN+nuN3eaoSitYXgXv8eIPPpahqpoEcg4wxxrLYtEGPIatSLMvimpMbMPnZFlx8TCZL\nJufJE0c/weijRjP317l+l6cqoGhdEZZtbTTaISDmaCDHKEvYsLFI36+qbJZlce0pDfjr2eZceHQm\niyfNZ/TQ0Yw+erSZ95ue4xrNitYVIZbo8eMYpIEcq4R12sta7YplWVx/agMmj23B/x1dj8V/zpdR\nQ0bxxDFPoMEcnTav2wyix49jkQZyjDKGtes1kFUF2bbFDac2ZNIzLbjgqHosmjifUUNGMebYMcz/\nfb7f5akIReuKMI5Z7XcdqvppIMeoUJiFSwvCmsiqUgIBixtPa8ikZ5pz/pB6LJwwl5GDR/Lk8U8y\nf4IGczRYs2SN44ScxX7XoaqfDp0Zu/LyVxaH0S9dqgoCAYubzmjIdafU55E31/D6+DmMPHIkHQ7t\nwJCbh7hjZdeyu/e5m9X5O+4Y9rugHyc+dOJOl5303iReufAVug7tynkvn7d1+vjR4/l2zLcA9L+i\nP4dddtjW+xZMXMB7N73HtV9fG1XDjxbmFYaBhX7XoaqfBnLsyluxOhxwHINlRc+HiapbAgGLm89s\nyA2n1ufhN1bzxrdzeHzQ4+x12F4MvnkwrfdrXWu1XDf+ushRqlg6fSlPnfgU+xy3z06XK8wv5OM7\nPqZt37bbTV8yfQmfP/Q5F751IcYxjD11LB37d6Rpp6aEQ2Hevf5dThl5SlSFseM4rF221kYDOSbp\n3lPsyguFkYK1Orqe2n2BgMUtw7OYPLYFZw3KYMEvs3l84OM8Pexp8v7Mq5Ua0hqmkZGdsfU29fOp\nNGrTiHYHtCt3GcdxePWiVxl882CyWmZtd9+K3BU027sZ7Q9sz54H7Umzzs1YMXsFAONHjadd33bs\n0X2PGn1MlbWxYCPh4rCFBnJM0kCOXXkASwpCftehYkggYHHrWVn8ObYFwwdmMO/nXB4b8BhPn/Q0\nCyfVXkaEi8P8+e6f9Dmzz07n++KhL0hvlE6fM3acr2nnpqycu5I1i9dQmF/Iynkr3WnzVvLHW38w\n5F9Daqr8Klu9aGuTvQZyDNIm69iVB7C4IMw+e/pdioo1iQGL287O4qYzGvDga6t5+/tcHj1iJp2O\n6MTgmwfTct+WNbr9f/73D0Xriuh9au9y55n32zx+f/13bvzxxjLvb9yhMUNvHcqTJzwJwNG3H01O\n+xyePP5JjrnzGGZ8PYMvHv4CO9Hm+PuP3+meeG3RQI5tGsgxyhizJmDLxiUFoTS/a1GxKzFgcfvZ\nWYw4owEPvFrIOz/M4tGvZ9BpQCeG3DyEPfapmSbf31/7nU5HdKJe43pl3r9lwxZeu+Q1Tnn8FFLr\np5a7nr7n9KXvOX23rff130nOSKbVfq24v/f9XPftdaxetJqXL3iZ2/+6HTvB34uMr160GrFki3HM\nKl8LUTVCAzmGWRb5iwtCHf2uQ8W+xIDFHec04uYzHe5/pZB3v5/Ff76aQeeBnRl88+BqPRa7etFq\ncr/P5fxXzy93noIFBRTmF/Lc6c9tveRkyXWEr2t8HbdMuIWsVtsfU96wagNf/vtLrvzkSvL+zCOn\nfQ6NWjeiUetGhIvDrJizgqadmlbb46iKFXNWYNnW/FA4pMPwxSAN5BhWHGL63MXFewL+fq1XcSMx\nYHHnuW4w3/dKIe99N5PpX05n7yP3ZvCIwbTo1mK3t/Hbq7+RkZNBpwGdyp2ncYfG3PTTTdtN++S+\nT9iyYQsnPHgC9ZvX32GZD//1IYdddhiZTTPJm5RHOLStQ6QTcrYGup8WT10cDheHJ/ldh6oZGsix\n7a9pC7YcY4yJqlM3VOxLSrC4+7xG3OIF8/vjZzDt82l0GdyFI0ccSYuuVQtmYwwT3phA79N6Y1nb\n90l97dLXyGyayVG3HUUgMUCTjk22uz+lXgoiQpO9tp8OMOvbWayct5Iznz4TgFY9WrFi9gpmfD2D\n1YtWYwUsctrnVKnm6mKMYem0pQD/+FqIqjEayLFt8obNJrCkIEzzbH2pVe1LTrS45/xG/Gu4w70v\nF/LBN9OZ+tlUugzpwuARg2nepXml1pf7XS5rFq+h9+k7duZas3gNUoVz7ouLinlvxHuc88I5W6dl\nNs3khAdP4I0r3iCQFOCMJ88gkOTve6gwv5Dg5qANTPG1EFVjRK/gFbtEpDmwaMw12QzYT/t2Kf8V\nBR3ufqmQj37ZRHHQoevQrgweMZhmezfzu7SoN/XzqTx3+nMALY0x+X7Xo6qfnocc25YEbFZPXxD0\nuw6lAHeP+f7/a8TEp1sw7JA0Zn41jYcPephxZ41jyfQlfpcX1ZZOX4plWxuARX7XomqGBnIMM8YY\nx2Hi9AVBbQZRUSUlyeL+C7P546kWnHhwGjO+nMrD/R5m3NnjWDp9qd/lRaUl05eAMMVos2bM0kCO\ncY5h0pR5QR2uS0Wl1GSLBy5yg/n4g9KY/vlUHur3EC+c8wJLZ2gwR8qbmFfshJyJ1bEuEflWRB71\nfp8vIldWx3oruO1a3V5dooEc+yYXrA0nFK7XMa1V9EpNtnjo4mz+eLoFx/VLY9pnU3io30O8eO6L\nGszA+pXrKVxYmAD8UgOr3w8YWwPrVZWkgRz7JgFMmbvF7zqU2qW0ZIuHL3GD+di+aUz99B83mM97\nkWUzl/ldnm8W/LGg5NdqD2RjzCpjTFF1r1dVngZy7JsTsFnx23T/32/9r17EXmcu2OF290vuKIDB\nYsOdL6yi90UL2ff8PK54fAWrKnG1qtueL2CvMxfw8ufrtk4LhgzXP7mSHhfkMej6RfwydfN2yzz7\n37Xc+7KOQhht0pItHrnUDeZj+qYx9ZN/ePDAB3npgpdYNiv+gnn+7/OxAtZyoNK9q0UkVUReFpH1\nIrJYRK4tdf92Tcgikikiz4nIChFZKyJfi0i3iPu7ich4EVnn3f+HiPSIuL+fiPwgIptEJE9ERopI\n+eOXqq00kGOcMcaEwnz+4z+bfT+O/N49TfllzB5bby+MaIwAg/u4p2Td90oh3/+1iSeuzuG125qw\nYk2Yy0euqNC6v5q4kX/mBmncYPtByd78Zj0z8oK8c3dTTj4sg+vGrNx6X/6KYt79fj3Xntyg2h6j\nql5pyRb/vjSbCU+14OgD0pjy3795sO+DvPx/L7M8d7nf5dWa2T/MDjkhZ3wVO3T9GzgIOBoYCBwK\n9NzJ/O8CWcAgoAduK9s3IlIyvNlruF8Menr3PwgUA4hIO+Az4B2gC3AKcCAwugp1xx0N5PjwTW5+\nccDv48gNMmyyMrfdvp28iZaNA/TqmMyGTQ7vfr+em89sSO9OyXRuncQDFzZiUu4W/t5Fc/uywhD3\nvlzIfy5rRMDefmCIeUuK6d8jlXbNEjlzQAaF6x1We8/DnS8UcuNpDUlN1rdBtEtPsfjPZdlMeLIF\nQ/dP5Z+P/+LBAx7k5QtfZvns2A7movVFLJqyyAa+r+yyIpIGnAdcZ4z5zhgzDTibcobTFZF+uMeU\nTzbGTDbGzDXG3AisAYZ5s7UEvjbGzPbuf88YUzJYyQjgVWPMaGPMPGPMb8DVwNkikljZ+uONfhLF\nh28Afpvmf7N1ieKQ4eOfNzLs0AwApszfQjgMffdO2TpP22YJNMuy+Wt2+YFsjOHGpwu44KhM2jff\n8f3esVUif84qYkvQ4Yd/NpNT36ZBhs1HP20gOVE4vKe2pNUl6akWj12ew4QnWzC4Twr/fPQXD+7/\nIK9c9Aor5lSsNaWumT9hPsYxQhUCGWgHJAATSiYYY1YDs8qZvxuQARR6TdzrRWQ90NpbF8CjwPMi\n8pWI3CQibSOW7w6cU2rZz7372lSh/rii4ynGAWPM4oSAzPl12ub2Q/aPjhG7vpq4iQ2bHI4/yK1n\n1dowCQEhPXX774hZmTYr15S/Z//Mx2tJsIXhA8u+DN+wQ9KZtTDI4BuX0LCexairslm7Mczo99fw\n6r+a8Njbq/n0t420bJzAAxdmkdNA3xJ1QXqqxeNX5LB+o8Otzxfw1YeTmfTuJHqc1IMjbziS7HbZ\nfpdYbWaOn4kVsFY4Iae8EN2ZkiajijZ1pwNLgEMili2xBsAYc5eIvAYMBYYAd4nIKcaYj7zlnwFG\nlrG8XsN5F/TTJ06Ewnz+w9+bL8L9tuy7975fz8HdU8iuv/N/QWOgvOtiTJ2/hVe+XM+H95V/SbyA\nLdxxzvaX2RvxTAFnDarH9Lwg4ydt4r8PNmPsf9dyz0uFjL7a3wsIqMrJSLMYeWUO6zaGufW5VXz1\nwWQmvTOJnif3ZNANg8huW7eD2RjDlE+mhJyQ82EVjx/PAULA/sB7ACLSAOgAfFfG/JOAJkDYGFNu\ngBpj5uCG7kgReR04F/jIW35vY8z8KtQa97TJOn58s3RVOCF/RbHfdbCkIMQvU4s4+bCMrdMaZdoU\nhwwbNjnbzVu4LkyjzLKvHvnnrC0UrgtzyJWL6HzWAjqftYDFBSEeeK2Q/leXPbrgr9M2M3dxkOED\nM/h9ehEH75NCcqLFkD5pTJgZPU36qnLqpdmMuiqH38c0Z+B+KUx+bxL3976fVy99lZXzVu56BVFq\nxewVFC4sDAD/rcryxpiNwPPAIyJymIh0AV4Aymx2MsZ8DfwKfCgiA0SklYj0FZF7RaSHiCSLyGgR\nOUREWorIgUAvYLq3ioeAA7x5uotIexE5VkS0U1cF6B5y/PgOcH6eWmSd2t/fneR3v19PVqbNIfts\nO17cpU0Stg2/TNvMwF5uM/b8pcUsWRVmnz2TylzPcQelcWCX5O2mnfvgco47KJ0TD0nfYf5gseHu\nFwt57IpsRATHGEIhd/c7GDI4zg6LqDqmXprN6KtzWLMhxK3PFvLNu5P48+0/2e/k/Rh0wyAatWnk\nd4mVMu2LaYglQeOY8buxmhuANOBjYD3wH6Ae25qxS+95DwHuA8YB2cAy4AdgOW6QZwEvAY2BAtw9\n7zsBjDFTROQQb/kfcJut5wJvRaxfh/4shwZynDDGrAnY8uNnv23sd2r/jLJ3OWunDj74YQMnHJyO\nFXGpvPRUi5MOzeDB11aTmWaRlmJx70uF9OiQRPd22wJ50PWLuOHUhhyxXyqZaTaZads/lISAkJ1p\n07rJjl86Rr+/hsP2TaFjS7fzV48OyTzyxmpOODidV79aT48OZQe/qnvqpwd44ho3mP/1bCHj3/2T\niW9PpPdpvRlw3QAata4bwTzlsylh4CtjzKaqrsPbSz7bu5X4T8T9bcuY/2rvVpbTd7G9P4Ejd3J/\n2/Lui3cayHEk7PDm79OLDi5cF6ZhPX8y+ZepRSxdFebEg3fcg73lzIZYViFXjFxJcchwULeUHY7/\n5i0LsX5z5XdlZy8K8uUfG/no/m2X+TuydyoTZhRxxj3LaNM0wH8ur9vHG9WO6qcHGHNNDqvXh/jX\ns6v49q0/mPDmBHqf2puB1w8kq1XWrlfik42rN7JgwgLLOOZjv2tRtUOvhxxHRCRbYNmd5za0Tjui\n7F7JSsWywnUhbnl2Fd//swXHGPqc1ocB1w8gq2X0BfPEdyby6kWvArQwxiz2ux5V8zSQ44xty9c9\nOyQd+tptTX1rtlbKb4XrQtw8dhU/TPGC+fQ+DLguuoJ57KljnZnjZ04OF4f387sWVTs0kOOMiJwv\n8NyPT7TQc25V3Fu1NsTNz67ixylbMAb6nOEGc8M9Gvpa14aCDdzW6TZjwuZKY8wTvhajao2e9hR/\nPkAIf/FHlfuIKBUzsjIDjL2+MT+NasZBXZOY8Prv3NvjXt665i1WL1rtW12T3p+EcYwDvOlbEarW\n6R5yHLIt+ax7+6QBb92pzdZKRVq5xm3K/nmau8e8/1n7M+CaATRoUbsXIPn3Yf8OLZ6y+DMn7BxT\nqxtWvtJAjkMicjbw4ncjW9CskTZbK1XayjUhRjyzip+nbwEDB5x9AEdcfUStBPPy3OU8sP8DACcZ\nY96t8Q2qqKGBHIdEJMOyWH7xMZkpV5+klx5UqjzLV4e4+ZkCfpkeBIG+Z/fliKuPoH7z+rteuIo+\nufcTvhn5zQYn7GQbY3T4uDiigRynRGRMg3Trwh/H7BFIDJQzWLRSCnCDecTTBfw6MwjAgeceyOFX\nHU79ZtUbzE7Y4a5ud4XWLl37gjHmwmpduYp6GshxSkT2BqY+elkjjuq74yAdSqkdLVsVYsTYAn6b\nEUQsoe+5fTniqiPIbJpZLev/55N/GDd8HEBvY8wf1bJSVWdoIMexgC0/dG2b1Pftu7Rzl1KVsbTA\nDebfZ7rBXLLHvLvBPPro0eH5v8//M1wc7lNNpao6RAM5jonIMOCdj+5vRqdWiX6Xo1Sds6QgxIhn\nCpgwywvm87xgblL5YF4yfQkP93sY4FRjzFu7ml/FHg3kOCYiCbbFohMPSc+594K6Mdi+UtFo8cpi\nRoxdxR+ztiC2Rb/z+nH4VYdTr3HFh6h986o3mfDGhBVOyGlhjPH/Oqmq1mkgxzkRuT0xQe74ZUwL\nq16atlwrtTsWryzmpqcLmDg76Abz+f04/MpdB/PGwo3c3vl2JxwM326Mua+WylVRRgM5zolIUxHy\nbzq9gX3ekOrpmKJUvFu0opibningz9lBLNum3wVuMGfkZJQ5/9ePf80n935SbBzT3BizspbLVVFC\nA1lhibzYoJ51xvcjWwSSEnU0VaWqS/4Kd4950hw3mA/6v4Pof2V/MrK3BXNxUTF3db8rtGHlhleM\nMef5WK7ymQayQkQ6ADPvPLehnK6XZVSq2uUtDzLi6VVMnhvECtgcfOHBHHb5YWRkZ/DTuJ9494Z3\nDYZOxphZfteq/KOBrACwLHkju749bPzjLXSgEKVqSN7yIDc9vYq/SoL5ooOZ+PbE8Lrl694xjjnN\n7/qUvzSQo5iIzAceM8aMqoVt7Q1Muef8LDmlf9nHuZRS1WPB0iA3PbOKv+cFccIGoIsxZprfdSl/\n6QFDBYAxZpoIb49+b00oWKxf0pSqSa2bJvLyLY1pkCYh4GMNYwUayJUmIgl+11BRImKJSIXbn43h\njhVrwvab49fXZFlKKeCt8RtYtc6xgRv8rkVFBw3kXRCRb0VktIg8JiIrgc9FJFNEnhORFSKyVkS+\nFpFupZa7VUSWe/c/KyIPiMjkUut9tNQyH4jIuJ3Uco2I/CMiG0RkoYiMEZG0iPvPFpHVInK0iEwD\nioA9KvpYjTGzBF4a88Ga0OYtTkUXU0pV0qYihzEfrAkJvGiMyfW7HhUdNJAr5ixgC9AXuBh4B8gC\nBgE9gEnA1yJSH0BEzgBuwf3m2xNYCFwC7G5bcBi4Atjbq+kw4KFS86QCNwLne/OtqMwGDNy9ZoPD\nuE/X7WapSqnyjPt0HWs2Ohi42+9aVPTQQK6YOcaYEcaY2UAO0As42Rgz2Rgz1xhzI7AWGObNfznw\nrDHmZWPMHGPMPcCU3S3CGDPKGPO9MSbPGPMdcBtwcqnZAsAlxpjfjDGzK3s9VWPMfGN49KkP1ziL\nV4Z2t2SlVCmLVhbz1EdrHGP4jzFmgd/1qOihgVwxEyN+7w5kAIUisr7kBrQG2nrz7AWUvnTahN0t\nQkSO8JrHF4nIOuAVIEtEUiJmCxpjpu7mpu4NO6x64LVC7d2lVDV74NXVxnEoAO71uxYVXTSQK2Zj\nxO/pwBKgG244l9z2Av4dMV/pMCvducopY1q5HcZEpBXwX+Av4ATcpvLLylhuc3nrqChjzPqww9Vf\n/rFJfp6626tTSnl+nrKZryZukrDD1caYDX7Xo6KLBnLlTQKaAGFjzLxSt0JvnllA71LL7Vfq75VA\n05I/RMQCuuxkuz0ByxhzvTFmgjFmDtB8tx7Jzr1hW/xy5wurQsGQ7igrtbuCIcOdL64K2RY/AW/6\nXY+KPhrIlWSM+Rr4FfhQRAaISCsR6Ssi94pID2+20cAFInKWiLQXkVtx96gjk208MFREhojIXsBT\nQP2dbHoOEBCRK0WkjYgMBy6q9gfoMcaYsMOlectC9itfaAcvpXbXK1+sI29ZyA47XGZ0RCZVBg3k\nXSvrjTME+AEYh7s3/DrQElgOYIx5HbgfeAT4E2gFvIh7GlKJccBL3u07YC5uSJe5bWPMP8C1uD2o\npwCnASN243HtkjHmb2DMyHfXhFes1g5eSlXV8tUhRr67Jgw84b2XldqBDp1ZS0TkS2CpMeZsv2up\nDBFpYFvMPbJPWv3HLs/WQa6VqoJrRq80n0/YuCbs0M4Ys9rvelR00j3kGiAiKd4gHp1FpOP/t3fv\nQVKVZx7Hf8/79tyHgREG0FCQCLgxRhQVYlzvSIJRd5cEQXPBW9j1ktVERlNqJZvdKIqI0WTFjRBc\nXSsKosYoMKaIiAYVMOEiLCoDkQm3GWboufR093Sf9332D8YUaynhNpzT3b9P1fmDKqbqmX/mW+c9\n73mPiPw7gDHYe5ecU1Q17jxuWfhWp9St6PzbP0BE/8/iFZ1Y+HanOI+bGWPaH94h9wARKcXeHdEj\nAaoSsSIAAAx0SURBVJRg77L2T1X1xVAHO0QiIkbwXEWZ+YfF9x9n+1fHwh6JKCc0xQNcfNt215nW\nF71iAp8d0/4wyHRARKTGGmw8++Sy6sdu628O4ohsooKkqpgyo8kvfze1x3mcqKrNYc9E0cYlazog\nqrrbeVy9bG3KzF/K1yeJ/pZ5SxN4fW3KOI+rGWM6EAwyHTBVfVmAOXc/ucc3NGXDHocoshoas7jn\nyT0Oe4/QXRj2PJQbGGQ6KArc6rxur53V7Jzn4w6ij3NeUTur2QVedwCYGvY8lDsYZDooqtoROHxr\nzaYu89hLbWGPQxQ5j/22DWvqu4xz+Kaq8uPidMAYZDpoqvoGgGkPPduqf1jHs66JPvLGuhQeerYV\nAO5R1T+EPQ/lFu6ypkMiItYYLCovkTG/mXacHdz/U7+LQVQQGpqy+Kc7d7hkly7xHpeoqgt7Jsot\nvEOmQ6Kqzntckc7othtmNgXJtA97JKLQJNMe1z/Q5FIZ3eY9rmSM6VAwyHTIVDUeOFy6eUc2uHN2\nM888oIKkqrjjsWbdsiObcQ6X8jQuOlQMMh0WVV3vPSYvejspcxfxq1BUeOYuasfiFUnxismquj7s\neSh3Mch02FT1WQDT7386rsvXc5MXFY7l61O4/+m4ArhPVReEPQ/lNm7qoiOie5NXXXmJXDD/J8fa\nYYOKwx6JqEfVb8tg4k92umSXLvUe4/jcmA4X75DpiOje5DUx1aXvX3VvY7Cjmd9Ppvy1fXeAydMa\ng1SXvu89LmeM6UhgkOmI6f5U49g97W7n5Gm7gj0d/BtF+WdPu8NV03YF8YTb6TwuUtXWsGei/MAg\n0xGlqjucx5jtu4O26+5rdIkUX4ei/JFIeVw7vdFtbw7anMOFqroz7JkofzDIdMSp6ibnMXZjQ6br\nxp81+UyW+xQo92WyihsfbPLvNWTSzmOsqtaHPRPlFwaZeoSqrvYel6z837SrnbVb+SEKymXOK6Y+\nsltXbEy77lO4Voc9E+UfBpl6jKq+5hWT6lYm8W9zW+AZZcpB3it+/KsWvLIqCVVMVNVlYc9E+YlB\nph6lqi8AuGb+0oTe/l/NGjhGmXJH4BS3Pdqsz76WUABXq+pvwp6J8hffQ6ajQkQmiuDXY04vNw99\nr0aKiyTskYj2K5NV3PKLJv/qn1Kqiiu7D8Ah6jEMMh01InKpETz/5S+W2lk/6G/KSrhAQ9GU6vK4\n4cEm/9aGtFPFeFVdGPZMlP8YZDqqRGSMMXh55PCSotm1A2xlOaNM0ZJIenx3RqNbU9+V7d7A9WrY\nM1FhYJDpqBORs6zBK58fXFz2+B0DbJ9KG/ZIRACA1oTDNfc2uvcaMinn8RVVfSvsmahwMMgUChE5\nzRosGTKwqGrO7f3toJqisEeiAveXpiymzGgKtu7KdjiPMXy1iY42BplCIyInWou6ilLzmUd/0N+O\nOrE07JGoQK3amMYNP2sKkmm/PXAYp6rvhT0TFR4GmUIlIn2twfMAzvmP6/rK5ef3CnskKjDzl3bg\nx3NbVIDXncc3VLUl7JmoMDHIFDoRKQLwcwDXXzWuF374zWMQs3wtinpW4BTTf70HT9R1AMCjAG5R\n1WzIY1EBY5ApMkTkRgF+8eUvluLnN9eYqgpu9qKe0dbpcPPDu/3bG9JQ4F9VdVbYMxExyBQpIjLG\nGjz/mZpY+S9r+8eGHlcc9kiUZzbvyOBfZjQF25uDpPMYz9eaKCoYZIocERkWs1hkjAz90eRjzMQL\nKiHCJWw6PKqKeUsTuPuJFu8V9YHDJfxiE0UJg0yRJCIVAjyowD9fdEa53vPdvlLdi0vYdGjiHQ53\nzWnRJe8kBcAvAdyqqsmw5yLaF4NMkSYi463B49W9bMWDN/WLnXlSWdgjUY55a0MKUx9pDuIdrtN5\nXNP9wROiyGGQKfJEZJA1eMp7nDflst64ZUIfFMW4hE37lwkUDy+IY85L7TAGy5zHt1R1e9hzEX0a\nBplygohYALUiuOcLQ4ox86Yae/xxPN2LPtmWHVnc+shut3FrBqq4C8ADqurCnotofxhkyikiMsoa\nzBfB4JvG9zFTLuuNYt4tU7dMoJj9UhseeaHVA9gaOExS1VVhz0V0IBhkyjkiUgbgRyK4/bMDYzpt\nSr/Y6X/HYzcL3R/fT+PO2c3Bh7sCUcX9AH6qqqmw5yI6UAwy5SwRGWENfuU8zphwXqXWXlEtx1Rx\nJ3ah2dPuMOOZuD63LCHW4B3nca2qvhv2XEQHi0GmnNb9bHmKNZheViIVtVdU20kX9oI1XMbOd84r\nnvl9Bx54Ju7SGe10Hj8EMJvPiilXMciUF0SkRoB7FbjuhEFFrvbKanveKWU8UCQPqSqWrU3hgafj\n7oNtWSvAHAXuUNXmsGcjOhwMMuUVETnTGsx0HmeNHF7iaq+otqM+z+fL+WLVxjRmPBN3a+q7rDV4\n03ncqqorwp6L6EhgkCnvyN7b4q9ag+nOY8TZJ5f6qZOqzUmfKwl7NDpE6//chZnz4n75u2ljDdZ2\nL0//TvkHjPIIg0x5S0QMgK/HLO4LHIZ+dXS5fv/yPsIPVuSOzTsyeGh+q76yKikxi/rA4Q4AzzHE\nlI8YZMp7IhID8J2Yxd3O49hxo8txzcVVcupwLmVH1epNafz34natW5mENdgZONwF4ClVDcKejain\nMMhUMESkBMB1MYvawOFzJx9f7K79Wm/7lVHlPIozArKB4pVVnXh8Ubt7d0vGxiz+HDjMADBXVbvC\nno+opzHIVHC6l7K/Zg2mOo/za/rY4OpxVbGJF1aidwXfYz7a2jod5r2awBN17cHuVhezBkudx0wA\ni1XVhz0f0dHCIFNBE5ERAL4vgu8UxcRMOK/STDi/Eid9tpivTPUgVcWGDzNY8FoCC5YlfBCo84r/\nAfCwqq4Lez6iMDDIRABEZACA663B95xHvyEDYsHXz62MXfb3FRhUw49YHCl/acri5Tc78fzriWBr\nYxCzBs3O4z8BPKqqTWHPRxQmBploH90bwMYA+LYRTPCK0pHDS9z4cyrtuC+Vo08ll7QPVmvCoW5F\nEi+8kXCrN3VZI0h7xQIATwH4PTdqEe3FIBN9ChGpBPCPxmCy9xhrDXD+qWW46IxyOWdEGfpXx8Ie\nMbKa4gHeWJfCkneS+tqaFJyHGoMl3uNJAC+qaiLsGYmihkEmOgAiMhDAJGvxbedwOgAZPqgouGBk\neezcU8owcnhJQe/UzgaK1Zu68PraFJauTgabtmVjANRa/NE5PAVgnqruCntOoihjkIkOkojUABgL\n4GJrcInzqC4rEXf2yWXm3FPKZPSJpRgyIAaTxx+48F6xtTHAyo1pLFub0uXrUj6VUWsN4s7jZQB1\n2HuSFs+XJjpADDLRYeh+heo0AONiFpcGDqMAmIpScacMLZERw0rMiONLMGJocU4vcTfFA6zbnMG6\nLV1YW9/l127u0mRaLQAfs1gZOCwEsBjAar6qRHRoGGSiI0hE+gAYBWC0CL5kDc4KHPoCQL/eNnva\nCSVFXxhSjMEDYxgyoAhDBsRQFaF3n9s6HRoaA2xtzKJhV4ANH2Z09aauoLnNFQFAzKLFebypihUA\nVgJYpaqt4U5NlB8YZKIe1P2hi0EARgMYbQ3OBHCq86j66P9UlkkwZGARjj+2KDZ4QAyDBxShb5VB\n70qL3hUGvSsMqirMYX3j2XlFe6dH20dXwqGl3aOhMYuGxgBbdmaDrbuySKT0r7fx1qAdwBrn8Tb2\nxnclgG08R5qoZzDIRCEQkWoAQ7uvYQCGWYMTRHBC4NDvk36molSCqgrj+1Ra6VUmxloRIxBrBdZA\nnIc6p/AKdU61I6W+NeG0vdObzrR+4np5zKJZFR84jw8A1HdfmwFsVtV4z/z2RPRJGGSiiBGRMgB9\nARyzz1X9sX/3AmD3uQwAD8Dtc3UA2LPPFf/Yv1tUNXW0fi8i2j8GmYiIKAJM2AMQERERg0xERBQJ\nDDIREVEEMMhEREQRwCATERFFAINMREQUAQwyERFRBDDIREREEcAgExERRQCDTEREFAEMMhERUQQw\nyERERBHAIBMREUUAg0xERBQBDDIREVEEMMhEREQRwCATERFFAINMREQUAQwyERFRBDDIREREEcAg\nExERRQCDTEREFAEMMhERUQQwyERERBHAIBMREUUAg0xERBQBDDIREVEEMMhEREQRwCATERFFAINM\nREQUAQwyERFRBDDIREREEcAgExERRQCDTEREFAEMMhERUQQwyERERBHAIBMREUUAg0xERBQBDDIR\nEVEE/B8RIcrtskuSOQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fdbad6b3650>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display_pie(total_gas_prods, gas_names, colors=[\"green\", \"goldenrod\", \"lightGreen\"],title='Gasoline Total Production')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "We see that the most produced gasoline type is by far regular.\n",
    "\n",
    "Now, let's plot the breakdown of oil blend quantities per gasoline type.\n",
    "We are using a multiple bar chart diagram, displaying all blend values for each couple of oil and gasoline type."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": true,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "sblends = [(gas_names[n], oil_names[o], round(blends[o,n].solution_value)) for n in range_gas for o in range_oil]\n",
    "\n",
    "blends_bd = pd.DataFrame(sblends)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABSkAAAHpCAYAAABuqLm0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3X+YVnWd//Hne9QQSpC8HQxjQlERMU3I1LY1w6TEcllt\nt0bZhUWBJGXFtjATZenKyFYxNn+blvHDvoHKqigpqJW60sKqaKBp6pCAOLgMygyp8Pn+cd9D99zM\nMMMwcAbn+biuuYbzOe/zOe9zz3WNc738nHMipYQkSZIkSZIkZaUs6wYkSZIkSZIkdWyGlJIkSZIk\nSZIyZUgpSZIkSZIkKVOGlJIkSZIkSZIyZUgpSZIkSZIkKVOGlJIkSZIkSZIyZUgpSZIkSZIkKVOG\nlJIkSZIkSZIyZUgpSZIkSZIkKVOGlJIkSe8zEbE5Ii7L6NyPRMTCLM7dnmzPzyAiXomIW3d2T5Ik\nSe2ZIaUkSVI7FxHDC6FX8dfrEbEwIr6YdX8lUkuKCsFc8fW8HRFPRsQ/7ewG20pEnBoRlzexO1H0\nWUTECRFxeUR0baR2My383CRJkt6v9sy6AUmSJLVIAiYCrwAB9ABGAPMi4ksppXnZtdYqCfhf4D/I\nX89HgHOBn0fEB1JKP82yuRYaAowF/r2RfZ2B94q2Pw1cBtwGrC+p7Us+qJQkSeqwDCklSZJ2Hw+k\nlJbUbxRuEX4dqAR2t5AS4LWU0qz6jYj4OfAnYDzQbkPKiOiSUqolH642KqX0Tulh26h9t616kyRJ\n2l15u7ckSdJuKqW0Dqij4Yq9RkVEz4i4NSJWR8TGiHg2IkaW1Hy2cOv1P0TEdyNiRUTURcRDEdGn\nkTlHR8SLEVEbEf8dEZ/ZweupBpYDjZ2rW0RcExFVhf7/GBHfjogoqvlYof+LIuLCwi3ltYXnZPYv\nme/jEXFbRLxUuMZVEfHTiPhwSd2kwpz9ImJmRLwJ/DYibiO/irL++ZObI2JT0XFbnklZuCX8ysKu\n+tvcN0VERWH/Vs+kjIiDIuJXEbE2IjZExBMRMaSkZrt+XpIkSe2ZKyklSZJ2H90iYj/yq/LKgXHA\nB4FfbOugiCgHngQ2AdOAauBU4JaI+FBKaVrJIRcXan8EdAMmANOBE4rmPAe4AfgdMBU4GPgv4E2g\nqjUXFxF7AB8F/q9kvDPwG6AncD2wgvzt0z8ADgAuKplqOPAh4CfA3sC/Agsi4uMppTcKNacABwG3\nAquB/sAY4Iji6+Svz4r8FfAC8B3yn///Fvr5PHA221gpCdwJHAZ8rdDL2sJ4fS8NnkdZ+Hk9Uej9\nx+Q/0+HAPRFxRkppbsn8zf68JEmS2jtDSkmSpN1DAAtKxjYCI1NKzb1N+4rC8Z8orL4EuCkiZgKT\nIuLGlNJfiuo7AUenlDYBRMQ64JqIOCKl9IeI2BP4PrAEGJRSeq9Q9wfgZloeUu5VCF0hHzZOIP+s\nzZ+U1H2TfKD4iZTSnwpjN0fEKuDfIuKqlNJrRfV9gENSSqsLfc0nH9JOAP6tUHNtSunq4pNExJPA\nzIj4m5TSYyU9PJVSGlZS/wLw+eJb1huTUloaEUvIh5RzU0rNfT7fAfYHPpNSeqJwrluAZ4CrgdKQ\ncps/r2bOJUmS1C54u7ckSdLuIQHnkV+5V79672HgpxExtJljzwDuAfaIiP3qv4Bfk195N6Ck/tb6\nwKvgt+RDzoML258kv5LzhvqAsuDnwDpa7gvkVxO+ASwFhpF/scy3S+q+UuihpqT/BeT/p/uJJfV3\n1QeUACml35MPKYcUjW0JZSOiU2G+JwvXWfp5JPKrRneVU4FF9QElQEppA3AT0Dsijiipb+7nJUmS\n1O65klKSJGn38fuSF+fcQX41408i4t6SwLC+Zn9gX2A0+duZSyXygWOxFSXb9bdfdy98/1jhuBcb\nTJTSexHxcguvBeC/ge+S/5v0SODSwjlKXzpzKPBx/np7dHP9v9hI3Qvkw04AIqI7MAn4asnxiXxw\nW2p7rmtHfYz8Z1NqWdH+4hWSzf28JEmS2j1DSkmSpN1USilFxCPkn015KH8NsYrV3zkznfxKx8Y8\nU7K9qdGqvz53sf572kZNS1SnlB4u/PvBiHgeuJf8cxuvKaorAx4EftjE/C+04Fylx/0KOJ78C22e\nBt4unGc+jd9tVNeCc2SluZ+XJElSu2dIKUmStHur/3vuQ03sfwN4C9ijBc+ubKlXyAdghwGP1g8W\nnlXZG3iqNZOmlOZFxKPAJYXnZNYHgy8BHyoKNJtzaBNjrxb63BcYBExMKX2/qP9DtrflnVT7KtC3\nkfF+RfslSZLeV3wmpSRJ0m6qEAp+gfzt0Y2toiSltBmYA5wZEf0bmSPXilP/D/nw8+uFHur9C/lb\ny3fED4EcMKpo7P8BJ0TE4NLiiOhWeCt4saER0bOo5lPAccC8wlD9ysPSv4XHs31h4obC/F1bWkvL\nPp95wKci4rj6gYj4IPlb9l/2ZTiSJOn9yJWUkiRJu4cAhkRE/Wq6cvIvz+kD/CCl9PY2jr0YOAl4\nMiJuJv88ww8DA8mvKNyuoLLw7MlLyb9M5uGI+CX5t2//C/lVj62WUnogIp4FLoqIawsvhPkRcDpw\nb0T8DFgMfBA4ivxLgXoDbxZN8yLwu4i4Htib/O3jbxTmIaX0VkT8Bvh2RHwAeA0YXLiG7blFenGh\n/j8LbxDflFL6ZTO1VxSeJfou8F9Fq0WLTQEqgQciYlrh2kaQfxblGdvRnyRJ0m7DkFKSJGn3kIB/\nL9reCCwHvp5SurmR2i0rAlNKawqrCS8D/p78W8LXAs+x9Zu0m1pJ2GA8pXRzRJQB3yL/XMelwJeB\n721jjiZ7LPEf5N/yfTZwe0qpLiJOBC4B/gH4J2A9+WdRXgbUlBx/O7AZuJB8mPskcEFK6fWimkrg\nP4Gx5MPD+cAXgZUt7B/gTmAa8LVCrwHUh5SlP4P/KQS7Xye/+rWMfCha1Ujtmog4gfyq0vPJB63P\nAF9KKT1Q0kOLfl6SJEntXaTk3y+SJEna/UXEx8i/hfvfUkpXZ92PJEmSWs5nUkqSJEmSJEnKlCGl\nJEmSJEmSpEwZUkqSJOn9ZFvPupQkSVI75TMpJUmSJEmSJGXKlZSSJEmSJEmSMmVIKUmSJEmSJClT\nhpSSJEmSJEmSMmVIKUmSJEmSJClThpSSJEmSJEmSMmVIKUmSJEmSJClThpSSJEmSJEmSMmVIKUmS\nJEmSJClThpSSJEmSJEmSMmVIKUmSJEmSJClThpSSJEmSJEmSMmVIKUmSJEmSJClThpSSJEmSJEmS\nMmVIKUmSJEmSJClThpSSJEmSJEmSMmVIKUmSJEmSJClThpSSJEmSJEmSMtXuQsqI+E5EbI6Iq4vG\nHimM1X9tiojrSo7rFRH3RcSGiFgdEVdGRFlJzUkRsTgiNkbECxExfFddlyRJkiRJkqTG7Zl1A8Ui\n4lhgFPB0ya4E3ARMBKIwVlt0XBkwD1gJHA/0BH4BvANcWqjpDdwLXAecBXweuCUiVqaUHtwpFyRJ\nkiRJkiSpWe1mJWVEfAiYDpwLrGukpDal9EZKaU3h6+2ifV8ADgfOTiktTSnNJx9ofiMi6oPY84A/\npZS+nVJ6PqV0LTAbGL/TLkqSJEmSJElSs9pNSAlcC9yTUlrYxP6zI+KNiFgaEVdEROeifccDS1NK\n1UVj84FuQP+imodK5pwPnNAGvUuSJEmSJElqpXZxu3dEfA34BPDJJkpmAK+Sv537KOBK4DDgK4X9\nBwCvlxzzetG+p7dR0zUiOqWU/tJIX/uRX6X5CrCx5VckSZIkSZIkCdgb6A3MTymtbaoo85AyIj4K\nXAOcklJ6t7GalNItRZvPRcRqYEFEHJRSermZU6Rtnb6Zmi+QD0glSZIkSZIktd7ZwMymdmYeUgID\ngf2BxRFRHxruAZwYEecDnVJKpSHik4XvhwAvA6uBY0tqehS+ry763qOkphxYn1J6p4neXgGYPn06\n/fr1a9nVqMXGjx/P1KlTs25DklTC38+S1D75+1mStDtatmwZw4YNg0LO1pT2EFI+BHy8ZOxnwDJg\nSiMBJcAx5Fc/ripsPwFcEhG5oudSDgZqCvPU15xaMs/gwnhTNgL069ePAQMGNH8l2i7dunXzc5Wk\ndsjfz5LUPvn7WZK0m9vmoxQzDylTShuAPxSPRcQGYG1KaVlEHAycBcwD1gJHA1cDj6aUni0c8uvC\nHL+IiAnAR4DvAT8puoX8BuD8iPghcCtwMvlnWg7ZmdcnSZIkSZIkadva09u9ixWvnnwH+Dz5N3Ev\nA34E/Ao4fUtxSpuBLwGbgMeB28mvxry8qOYV4LTCXE8B44FzUkqlb/yWJEmSJEmStAtlvpKyMSml\nQUX//jNwUguOWUE+qNxWzaPkn4EpSZIkSZIkqZ1orysp1QFUVlZm3YIkqRH+fpak9snfz5Kk97N2\nuZJSHYN/ZElS++TvZ0lqn/z9LGl3V1VVRXV1dfOF2q3kcjkqKip2eB5DSkmSJEmSJO1UVVVV9OvX\nj9ra2qxbURvr0qULy5Yt2+Gg0pBSkiRJkiRJO1V1dTW1tbVMnz6dfv36Zd2O2siyZcsYNmwY1dXV\nhpSSJEmSJEnaPfTr148BAwZk3YbaIV+cI0mSJEmSJClThpSSJEmSJEmSMmVIKUmSJEmSJClThpSS\nJEmSJEmSMmVIKUmSJEmSJClTvt1bkiRJkiRJmaqqqqK6ujrTHnK5HBUVFZn2sCPKysqYNGkSl112\nWdattIohpSRJkiRJkjJTVVXF4Yf3pa5uY6Z9dO68N8uXP79bB5Wtdf311/Pwww/z5JNPsmLFCkaM\nGMGtt966S3swpJQkSZIkSVJmqqurqavbyBlnQC6XVQ9w550bqa6u7pAh5ZVXXsnbb7/Npz71KVav\nXp1JD4aUkiRJkiRJylwuBz17Zt3FrlVbW0uXLl2yboPf/OY39OrVC4B99tknkx58cY4kSZIkSZLU\nBlauXMk555zDgQceyN57783BBx/M2LFjee+99/jZz35GWVkZv/nNbxg7diw9evTYEgyOGDGCgw46\naKv5Jk2aRFlZw/junXfeYfz48ZSXl9O1a1eGDh3Ka6+91mQ/I0eO5IADDmDvvffmyCOPbPQ27vo+\nsuRKSkmSJEmSJGkHrVq1imOPPZb169czZswY+vbty2uvvcbs2bOpra0lIgAYO3Ys5eXlXH755dTW\n1gIQEVv2F2ts/JxzzmHmzJmcffbZnHDCCSxcuJDTTjttq7o1a9Zw3HHHscceezBu3DhyuRz3338/\n5557Lm+//Tbjxo3bSZ9E6xhSSpIkSZIkSTvo4osvZs2aNSxatIhjjjlmy/ikSZMa1OVyORYsWNBo\nKNmcZ555hhkzZnD++eczbdo0AM477zyGDRvG0qVLG9RecsklpJR46qmn2HfffQEYPXo0Z511FpMm\nTWLMmDF06tRpu3vYWQwpJUnajVRVVVFdXZ11G20il8t1yIeSS5Ik6f0npcTcuXM5/fTTGwSUpSKC\nUaNGtSqgBJg3bx4RwQUXXNBg/MILL2TmzJkNxu68806++tWvsmnTJtauXbtlfPDgwfzyl79kyZIl\nnHDCCa3qY2cwpJQkaTdRVVXF4f0Op662LutW2kTnLp1Zvmy5QaUkSZJ2e2+88Qbr16+nf//+zdb2\n7t271ed59dVXKSsro0+fPg3G+/btu1U/69at46abbuLGG2/cap6IYM2aNa3uY2cwpJQkaTdRXV1N\nXW0dw24cRo/DemTdzg55/YXXmT5mOtXV1YaUkiRJ2u2llFpc27lz563GmlpZuWnTpladZ/PmzQAM\nGzaM4cOHN1pz1FFHtWiuXcWQUpKk3UyPw3rQ6+js374nSZIkKa/+TdvPPvtsq47v3r0769at22r8\nlVdeabDdu3dvNm/ezEsvvcShhx66ZXz58uUN6vbff3/22WcfNm3axKBBg1rV065W1nyJJEmSJEmS\npKZEBEOHDuWee+5hyZIl2318nz59qKmpaRByrlq1irvvvrtB3amnnkpKactLc+pdc801DVZjlpWV\nceaZZzJnzhyee+65rc7XHp9z70pKSZIkSZIkZS7L3Kwtzn3FFVfw4IMPcuKJJzJ69Gj69evHypUr\nmT17No899hjQ9O3alZWVTJgwgaFDhzJu3Dg2bNjADTfcQN++fRuEnkcffTSVlZVcd911rFu3jk9/\n+tMsWLCAl156aau5p0yZwiOPPMJxxx3HqFGjOOKII3jzzTdZvHgxCxcubBBU3nvvvTz99NOklHj3\n3Xd5+umn+f73vw/A3/3d33HkkUfu+AfUDENKSZIkSZIkZSaXy9G5897ceefGTPvo3Hlvcrlcq4/v\n2bMnTz75JBMnTmTmzJmsX7+eAw88kCFDhtClSxeg6WdPdu/enbvvvpuLLrqICRMmcNBBBzFlyhRe\neOGFrVZm3nbbbZSXlzNjxgzmzp3LySefzH333UevXr0azF9eXs6iRYuYPHkyd911F9dffz377bcf\n/fv358orr2ww55w5c7j99tu3bD/11FM89dRTAPTq1WuXhJSxPQ/27GgiYgCwePHixQwYMCDrdiRJ\nHdySJUsYOHAg33z4m7v9MylXPL2Cqz53Ff43VpIkqWOo/1u2qb//qqqqMr8FOZfL+VLH7dTcz7W4\nBhiYUmryXnhXUkqSJEmSJClTFRUVBoQdnC/OkSRJkiRJkpQpQ0pJkiRJkiRJmTKklCRJkiRJkpQp\nQ0pJkiRJkiRJmTKklCRJkiRJkpQpQ0pJkiRJkiRJmWp3IWVEfCciNkfE1UVjnSLi2oiojoi3ImJ2\nRJSXHNcrIu6LiA0RsToiroyIspKakyJicURsjIgXImL4rrouSZIkSZIkSY1rVyFlRBwLjAKeLtl1\nDXAacCZwItATmFN0XBkwD9gTOB4YDowAJhfV9AbuBRYARwM/Bm6JiFN2xrVIkiRJkiRJapl2E1JG\nxIeA6cC5wLqi8a7ASGB8SunRlNL/Av8C/E1EfKpQ9gXgcODslNLSlNJ8YCLwjYjYs1BzHvCnlNK3\nU0rPp5SuBWYD43fF9UmSJEmSJElqXLsJKYFrgXtSSgtLxj9JfoXkgvqBlNLzQBVwQmHoeGBpSqm6\n6Lj5QDegf1HNQyVzzy+aQ5IkSZIkSVIG2kVIGRFfAz4BfKeR3T2Ad1JK60vGXwcOKPz7gMJ26X5a\nUNM1Ijq1pm9JkiRJkiSpPSgrK2Py5MnNF7ZTezZfsnNFxEfJP3PylJTSu9tzKJBaULetmmhBDePH\nj6dbt24NxiorK6msrGzB6SVJkiRJkrQtVVVVVFdXN1+4E+VyOSoqKjLtIQt//vOf+elPf8q8efP4\n4x//yB577MGRRx7JpZdeysknn7xdc82aNYtZs2Y1GKupqWnRsZmHlMBAYH9gcUTUh4Z7ACdGxPnA\nF4FOEdG1ZDVlOX9dGbkaOLZk3h5F++q/9yipKQfWp5Te2VaDU6dOZcCAAS29HkmSJEmSJLVQVVUV\nh/c7nLraukz76NylM8uXLe9wQeXcuXP50Y9+xNChQxkxYgTvvfcet99+O6eccgq33XYbw4cPb/Fc\njS3qW7JkCQMHDmz22PYQUj4EfLxk7GfAMmAK8BrwLnAycBdARBwGVACPF+qfAC6JiFzRcykHAzWF\neeprTi05z+DCuCRJkiRJkjJQXV1NXW0dw24cRo/DSteX7Rqvv/A608dMp7q6epeGlLW1tXTp0mWX\nna8xgwYNoqqqig9/+MNbxsaMGcMnPvEJLrvssu0KKXdE5iFlSmkD8IfisYjYAKxNKS0rbP8UuDoi\n/g94C5gGPJZS+n3hkF8X5vhFREwAPgJ8D/hJ0S3kNwDnR8QPgVvJh55fAYbszOuTJEmSJElS83oc\n1oNeR/fKuo0dsnLlSiZOnMgDDzzA2rVr6dmzJ1/84heZNm0a06dPZ+TIkTzyyCPccccdzJkzh/fe\ne4+1a9cyYsQIHn30UV5++eUG802aNInJkyezefPmLWPvvPMOEyZMYMaMGWzcuJFBgwZx7bXXNtnP\npZdeyrx581i3bh2HHHIIF110ESNHjtxS069fv62O+8AHPsCQIUOYOnUqGzZs4IMf/GAbfUJNyzyk\nbELpMyLHA5uA2UAn4AHgG1uKU9ocEV8Crie/unID+dWYlxfVvBIRpwFXA+OAPwPnpJRK3/gtSZIk\nSZIkbZdVq1Zx7LHHsn79esaMGUPfvn157bXXmD17NrW1tdQ/5XDs2LGUl5dz+eWXU1tbC0BE8Nen\nIP5VY+PnnHMOM2fO5Oyzz+aEE05g4cKFnHbaaVvVrVmzhuOOO4499tiDcePGkcvluP/++zn33HN5\n++23GTduXLPX06VLl1220rNdhpQppUEl238BLih8NXXMCuBLzcz7KPlnYEqSJEmSJElt5uKLL2bN\nmjUsWrSIY445Zsv4pEmTGtTlcjkWLFjQaCjZnGeeeYYZM2Zw/vnnM23aNADOO+88hg0bxtKlSxvU\nXnLJJaSUeOqpp9h3330BGD16NGeddRaTJk1izJgxdOrUqdHzvPjii9x111189atfbVWfrVG2S84i\nSZIkSZIkvU+llJg7dy6nn356g4CyVEQwatSoVgd/8+bNIyK44IKG6/guvPBCUmp4Y/Kdd97Jl7/8\nZTZt2sTatWu3fA0ePJiamhqWLFnS6Dnq6ur4h3/4B7p06cIVV1zRqj5bo12upJQkSZIkSZJ2F2+8\n8Qbr16+nf//+zdb27t271ed59dVXKSsro0+fPg3G+/btu1U/69at46abbuLGG2/cap6IYM2aNVuN\nb968ma997WssX76cBx54gI985COt7nV7GVJKkiRJkiRJO6B0FeO2dO7ceauxplZWbtq0qVXnqX/R\nzrBhw5p8O/dRRx211di5557Lfffdx8yZM/nsZz/bonO1FUNKSZIkSZIkaQeUl5fTtWtXnn322VYd\n3717d9atW7fV+CuvvNJgu3fv3mzevJmXXnqJQw89dMv48uXLG9Ttv//+7LPPPmzatIlBgxq8+qVJ\n3/rWt/j5z3/Oj3/8Y/7xH/9x+y9iB/lMSkmSJEmSJGkHRARDhw7lnnvuafJZj9vSp08fampqGoSc\nq1at4u67725Qd+qpp5JS2vLSnHrXXHNNg9WYZWVlnHnmmcyZM4fnnntuq/NVV1c32P7Rj37EVVdd\nxXe/+13OP//87e6/LbiSUpIkSZIkSZl7/YXXd+tzX3HFFTz44IOceOKJjB49mn79+rFy5Upmz57N\nY489BjR9u3ZlZSUTJkxg6NChjBs3jg0bNnDDDTfQt2/fBqHn0UcfTWVlJddddx3r1q3j05/+NAsW\nLOCll17aau4pU6bwyCOPcNxxxzFq1CiOOOII3nzzTRYvXszChQu3BJV33XUXEyZM4LDDDqNv377M\nmDGjwTyDBw9m//333+HPpzmGlJIkSZIkScpMLpejc5fOTB8zPdM+OnfpTC6Xa/XxPXv25Mknn2Ti\nxInMnDmT9evXc+CBBzJkyBC6dOkCNP3sye7du3P33Xdz0UUXMWHCBA466CCmTJnCCy+8sNXKzNtu\nu43y8nJmzJjB3LlzOfnkk7nvvvvo1atXg/nLy8tZtGgRkydP5q677uL6669nv/32o3///lx55ZVb\n6p555hkigj/+8Y/88z//81a9Pfzww7skpIztebBnRxMRA4DFixcvZsCAAVm3I0nq4JYsWcLAgQP5\n5sPfpNfRvbJuZ4eseHoFV33uKvxvrCRJUsdQ/7dsU3//VVVVbXUL8q6Wy+WoqKjItIfdTXM/1+Ia\nYGBKqcl74V1JKUmSJEmSpExVVFQYEHZwvjhHkiRJkiRJUqYMKSVJkiRJkiRlypBSkiRJkiRJUqYM\nKSVJkiRJkiRlypBSkiRJkiRJUqYMKSVJkiRJkiRlypBSkiRJkiRJUqYMKSVJkiRJkiRlypBSkiRJ\nkiRJUqYMKSVJkiRJkiRlypBSkiRJkiRJ2s2VlZUxefLkrNtotT2zbkCSJEmSJEkdW1VVFdXV1Zn2\nkMvlqKioyLSHLGzcuJFvfOMbLFq0iBUrVrBp0yb69OnDyJEjGTt2LHvuuWviQ0NKSZIkSZIkZaaq\nqorDDz+curq6TPvo3Lkzy5cv73BBZV1dHcuWLeO0006jd+/elJWV8fjjjzN+/HgWLVrE9OnTd0kf\nhpSSJEmSJEnKTHV1NXV1dfztqNF069kzkx5qVq7ktzffRHV19S4NKWtra+nSpcsuO19junfvzuOP\nP95gbPTo0XTt2pVrr72Wq666ih49euz0PgwpJUmSJEmSlLluPXuy38d6Z93GDlm5ciUTJ07kgQce\nYO3atfTs2ZMvfvGLTJs2jenTpzNy5EgeeeQR7rjjDubMmcN7773H2rVrGTFiBI8++igvv/xyg/km\nTZrE5MmT2bx585axd955hwkTJjBjxgw2btzIoEGDuPbaa5vs59JLL2XevHmsW7eOQw45hIsuuoiR\nI0c2ey0f+9jHAKipqTGklCRJkiRJknYHq1at4thjj2X9+vWMGTOGvn378tprrzF79mxqa2uJCADG\njh1LeXk5l19+ObW1tQBExJb9xRobP+ecc5g5cyZnn302J5xwAgsXLuS0007bqm7NmjUcd9xx7LHH\nHowbN45cLsf999/Pueeey9tvv824ceMa1L/77rusX7+euro6fv/733PVVVfRu3dvDjnkkLb8mJpk\nSClJkiRJkiTtoIsvvpg1a9awaNEijjnmmC3jkyZNalCXy+VYsGBBo6Fkc5555hlmzJjB+eefz7Rp\n0wA477zzGDZsGEuXLm1Qe8kll5BS4qmnnmLfffcF8rdxn3XWWUyaNIkxY8bQqVOnLfV33nknlZWV\nW7aPPfZYbr31VsrKyra7z9bYNWeRJEmSJEmS3qdSSsydO5fTTz+9QUBZKiIYNWpUqwJKgHnz5hER\nXHDBBQ3GL7zwQlJKDcbuvPNOvvzlL7Np0ybWrl275Wvw4MHU1NSwZMmSBvWDBg3ioYceYvbs2Zx3\n3nnstddevP32263qszVcSSlJkiRJkiTtgDfeeIP169fTv3//Zmt79+7d6vO8+uqrlJWV0adPnwbj\nffv23arhmnsmAAAgAElEQVSfdevWcdNNN3HjjTduNU9EsGbNmgZj+++/P4MGDQLgjDPO4Ac/+AGn\nnHIKL774IuXl5a3uuaUMKSVJkiRJkqQdULqKcVs6d+681VhTKys3bdrUqvPUv2hn2LBhDB8+vNGa\no446aptzfOUrX+G73/0uc+fOZdSoUS06744wpJQkSZIkSZJ2QHl5OV27duXZZ59t1fHdu3dn3bp1\nW42/8sorDbZ79+7N5s2beemllzj00EO3jC9fvrxB3f77788+++zDpk2btqyO3F51dXVA/u3eu4LP\npJQkSZIkSZJ2QEQwdOhQ7rnnnq2e9dgSffr0oaampkHIuWrVKu6+++4GdaeeeioppS0vzal3zTXX\nNFiNWVZWxplnnsmcOXN47rnntjpfdXX1ln+vXbu20Z5uvvlmIoJPfvKT2309reFKSkmSJEmSJGkH\nXXHFFTz44IOceOKJjB49mn79+rFy5Upmz57NY489BjR9u3ZlZSUTJkxg6NChjBs3jg0bNnDDDTfQ\nt2/fBqHn0UcfTWVlJddddx3r1q3j05/+NAsWLOCll17aau4pU6bwyCOPcNxxxzFq1CiOOOII3nzz\nTRYvXszChQu3BJXTp0/nhhtuYOjQoRx88MG89dZbzJ8/n4ceeojTTz+dk046aed8YCUyDykj4uvA\neUDvwtBzwOSU0gOF/Y8AJxYdkoAbU0pji+boBdwAnAS8BdwOXJxS2lxUcxJwFdAfqAK+n1L6+c64\nJkmSJEmSJG2fmpUrd+tz9+zZkyeffJKJEycyc+ZM1q9fz4EHHsiQIUPo0qUL0PSzJ7t3787dd9/N\nRRddxIQJEzjooIOYMmUKL7zwwlYrM2+77TbKy8uZMWMGc+fO5eSTT+a+++6jV69eDeYvLy9n0aJF\nTJ48mbvuuovrr7+e/fbbj/79+3PllVduqfvMZz7DE088wR133MHq1avZa6+96Nu3L1OnTuX888/f\n4c+lpTIPKYEVwATgxcL2CGBuRHwipbSMfCh5EzARqP+ka+sPjogyYB6wEjge6An8AngHuLRQ0xu4\nF7gOOAv4PHBLRKxMKT248y5NkiRJkiRJ25LL5ejcuTO/vfmmTPvo3LkzuVxuh+b46Ec/ym233dbo\nvuHDhzf5EhuAk08+maeffnqr8csvv7zB9gc+8AGmTp3K1KlTG4yXvmQH8p/ttGnTtro9vNjAgQO5\n4447mty/q2QeUqaU7isZujQiziMfOC4rjNWmlN5oYoovAIcDn0spVQNLI2IiMCUiJqWU3iO/UvNP\nKaVvF455PiI+A4wHDCklSZIkSZIyUlFRwfLlyxs8JzELuVyOioqKTHvoyDIPKYsVVkX+I9AFeLxo\n19kR8U/AauAe4HsppbrCvuOBpYWAst584Hryt3Y/Xah5qOR084GpSJIkSZIkKVMVFRUGhB1cuwgp\nI+JI4Algb/LPlPz7lNLzhd0zgFfJ3859FHAlcBjwlcL+A4DXS6Z8vWjf09uo6RoRnVJKf2m7q5Ek\nSZIkSZK0PdpFSAksB44G9gXOBG6PiBNTSstTSrcU1T0XEauBBRFxUErp5WbmbfyVSXnRghoAxo8f\nT7du3RqMVVZWUllZ2dyhkiRJkiRJUocwa9YsZs2a1WCspqamRce2i5Cy8NzIPxU2l0TEp4B/Jf8s\nyVJPFr4fArxM/hbwY0tqehS+ry763qOkphxYn1J6p7n+pk6dyoABA5orkyRJkiRJkjqsxhb1LVmy\nhIEDBzZ7bNnOamoHlQGdmth3DPnVj6sK208AH4+I4tcvDQZq+OuLd54ATi6ZZ3BhXJIkSZIkSVKG\nMl9JGRHfB+4HVgD7AGcDnwUGR8TBwFnAPGAt+VvCrwYeTSk9W5ji18AfgF9ExATgI8D3gJ+klN4t\n1NwAnB8RPwRuJR9YfgUYsvOvUJIkSZIkSdK2ZB5Skr8N+3by4WIN8AwwOKW0MCI+Cnye/K3fHyQf\nZP4K+H79wSmlzRHxJfJv834c2AD8DLi8qOaViDiNfMA5DvgzcE5KqfSN35IkSZIkSZJ2scxDypTS\nudvY92fgpBbMsQL4UjM1jwLN3wAvSZIkSZKknWLZsmXNF2m30ZY/z8xDSkmSJEmSJL2/5XI5unTp\nwrBhw7JuRW2sS5cu5HK55gubYUgpSZIkSZKknaqiooJly5ZRXV2ddStqY7lcjoqKih2ex5BSkiRJ\nkiRJO11FRUWbhFl6fyrLugFJkiRJkiRJHZshpSRJkiRJkqRMGVJKkiRJkiRJypQhpSRJkiRJkqRM\nGVJKkiRJkiRJypQhpSRJkiRJkqRMGVJKkiRJkiRJypQhpSRJkiRJkqRMGVJKkiRJkiRJypQhpSRJ\nkiRJkqRMGVJKkiRJkiRJypQhpSRJkiRJkqRMGVJKkiRJkiRJypQhpSRJkiRJkqRMGVJKkiRJkiRJ\nypQhpSRJkiRJkqRMGVJKkiRJkiRJypQhpSRJkiRJkqRMGVJKkiRJkiRJytSeWTcgSZK0LVVVVVRX\nV2fdRpvJ5XJUVFRk3YYkSZLUrhhSSpKkdquqqorDDz+curq6rFtpM507d2b58uUGlZIkSVIRQ0pJ\nktRuVVdXU1dXx9+OGk23nj2zbmeH1axcyW9vvonq6mpDSkmSJKmIIaUkSWr3uvXsyX4f6511G5Ik\nSZJ2El+cI0mSJEmSJClThpSSJEmSJEmSMmVIKUmSJEmSJClThpSSJEmSJEmSMmVIKUmSJEmSJClT\nhpSSJEmSJEmSMpV5SBkRX4+IpyOipvD1eER8sWh/p4i4NiKqI+KtiJgdEeUlc/SKiPsiYkNErI6I\nKyOirKTmpIhYHBEbI+KFiBi+q65RkiRJkiRJUtMyDymBFcAEYGDhayEwNyL6FfZfA5wGnAmcCPQE\n5tQfXAgj5wF7AscDw4ERwOSimt7AvcAC4Gjgx8AtEXHKTrsqSZIkSZIkSS2yZ9YNpJTuKxm6NCLO\nA46PiNeAkcDXUkqPAkTEvwDLIuJTKaVFwBeAw4HPpZSqgaURMRGYEhGTUkrvAecBf0opfbtwjucj\n4jPAeODBnX6RkiRJkiRJkprUHlZSbhERZRHxNaAL8AT5lZV7kl8BCUBK6XmgCjihMHQ8sLQQUNab\nD3QD+hfVPFRyuvlFc0iSJEmSJEnKSLsIKSPiyIh4C/gLcB3w9yml5cABwDsppfUlh7xe2Efh++uN\n7KcFNV0jolMbXIIkSZIkSZKkVsr8du+C5eSfFbkv+WdP3h4RJ26jPoDUgnm3VRMtqAFg/PjxdOvW\nrcFYZWUllZWVLWhBkiRJkiRJev+bNWsWs2bNajBWU1PTomPbRUhZeG7knwqbSyLiU8C/Av8P+EBE\ndC1ZTVnOX1dGrgaOLZmyR9G++u89SmrKgfUppXea62/q1KkMGDCgRdciSZIkSZIkdUSNLepbsmQJ\nAwcObPbYdnG7dyPKgE7AYuA94OT6HRFxGFABPF4YegL4eETkio4fDNQAy4pqTqahwYVxSZIkSZIk\nSRnKfCVlRHwfuB9YAewDnA18FhicUlofET8Fro6I/wPeAqYBj6WUfl+Y4tfAH4BfRMQE4CPA94Cf\npJTeLdTcAJwfET8EbiUfWH4FGLIrrlGSJEmSJElS0zIPKcnfhn07+XCxBniGfEC5sLB/PLAJmE1+\ndeUDwDfqD04pbY6ILwHXk19duQH4GXB5Uc0rEXEacDUwDvgzcE5KqfSN35IkSZIkSZJ2scxDypTS\nuc3s/wtwQeGrqZoVwJeamedRoPkb4CVJkiRJkiTtUu31mZSSJEmSJEmSOghDSkmSJEmSJEmZMqSU\nJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmS\nJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClD\nSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmS\nJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZ\nMqSUJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZMqSUJEmSJEmSlClDSkmS\nJEmSJEmZMqSUJEmSJEmSlClDSkmSJEmSJEmZyjykjIjvRMSiiFgfEa9HxF0RcVhJzSMRsbnoa1NE\nXFdS0ysi7ouIDRGxOiKujIiykpqTImJxRGyMiBciYviuuEZJkiRJkiRJTcs8pAT+FvhP4Djg88Be\nwK8jonNRTQJuAnoABwAfAb5dv7MQRs4D9gSOB4YDI4DJRTW9gXuBBcDRwI+BWyLilJ1yVZIkSZIk\nSZJaZM+sG0gpDSnejogRwBpgIPC7ol21KaU3mpjmC8DhwOdSStXA0oiYCEyJiEkppfeA84A/pZTq\nw83nI+IzwHjgwTa7IEmSJEmSJEnbpT2spCy1L/mVk2+WjJ8dEW9ExNKIuKJkpeXxwNJCQFlvPtAN\n6F9U81DJnPOBE9qudUmSJEmSJEnbK/OVlMUiIoBrgN+llP5QtGsG8CqwEjgKuBI4DPhKYf8BwOsl\n071etO/pbdR0jYhOKaW/tNV1SJIkSZIkSWq5dhVSAtcBRwB/UzyYUrqlaPO5iFgNLIiIg1JKLzcz\nZ9rGvmhBDePHj6dbt24NxiorK6msrGzm1JIkSZIkSVLHMGvWLGbNmtVgrKampkXHtpuQMiJ+AgwB\n/jaltKqZ8icL3w8BXgZWA8eW1PQofF9d9L1HSU05sD6l9M62TjZ16lQGDBjQTEuSJEmSJElSx9XY\nor4lS5YwcODAZo9tF8+kLASUf0f+xTdVLTjkGPKrH+vDzCeAj0dErqhmMFADLCuqOblknsGFcUmS\nJEmSJEkZyTykjIjrgLOBs4ANEdGj8LV3Yf/BEXFpRAyIiI9FxOnAz4FHU0rPFqb5NfAH4BcRcVRE\nfAH4HvCTlNK7hZobgD4R8cOI6BsRY8k/0/LqXXe1kiRJkiRJkkplHlICXwe6Ao+QfzFO/dc/Fva/\nA3ye/Ju4lwE/An4FnF4/QUppM/AlYBPwOHA78DPg8qKaV4DTCnM9BYwHzkkplb7xW5IkSZIkSdIu\nlPkzKVNK2wxKU0p/Bk5qwTwryAeV26p5FGj+JnhJkiRJkiRJu0x7WEkpSZIkSZIkqQNr8UrKiGjx\nsxtTShe1rh1JkiRJkiRJHc323O59TMn2QGAP4PnC9mHknwm5uA36kiRJkiRJktRBtDikTCl9rv7f\nEXER8BYwPKX0f4Wx7sBtwG/buklJkiRJkiRJ71+tfSblN4Hv1AeUAIV/X1rYJ0mSJEmSJEkt0tqQ\nsiuwfyPj+wP7tL4dSZIkSZIkSR1Na0PKu4DbIuKMiPhoRBwYEWcCPwXubLv2JEmSJEmSJL3fbc+L\nc4p9HfgPYCawV2HsPfIh5bfaoC9JkiRJkiRJHUSrQsqUUi0wNiK+BfQBAngxpbShLZuTJEmSJEmS\n9P7X2pWUABRCyWfaqBdJkiRJkiRJHVCrQsqI+CBwMXAyUE7Jsy1TSgfveGuSJEmSJEmSOoLWrqS8\nBfgs8AtgFZDarCNJkiRJkiRJHUprQ8pTgdNSSo+1ZTOSJEmSJEmSOp6y5ksa9X/Am23ZiCRJkiRJ\nkqSOqbUh5URgckR0actmJEmSJEmSJHU8rb3d+5tAH+D1iHgFeLd4Z0ppwA72JUmSJEmSJKmDaG1I\neXebdiFJkiRJkiSpw2pVSJlS+ve2bkSSJEmSJElSx9TaZ1ISEftGxLkR8YOI+HBhbEBEHNh27UmS\nJEmSJEl6v2vVSsqIOAp4CKgBegM3k3/b9xlABfDPbdSfJEmSJEmSpPe51q6kvBr4WUrpUGBj0fg8\n4MQd7kqSJEmSJElSh9HakPJY4MZGxl8DDmh9O5IkSZIkSZI6mtaGlH8BujYyfhjwRuvbkSRJkiRJ\nktTRtDak/C/gsojYq7CdIqIC+CEwp006kyRJkiRJktQhtDak/CbwIWAN0Bl4FHgReAv4btu0JkmS\nJEmSJKkjaNXbvVNKNcApEfEZ4CjygeWSlNJDbdmcJEmSJEmSpPe/VoWU9VJKvwN+10a9SJIkSZIk\nSeqAWhxSRsS4ltamlKa1rh1JklqmqqqK6urqrNtoE7lcjoqKiqzbkCRJkqTMbM9KyvEtrEuAIaUk\naaepqqri8MP7Ule3MetW2kTnznuzfPnzBpWSJEmSOqwWh5QppYMaG4+IKOxPbdVUe7Ns2bKsW2gT\nrtSR9H5RXV1NXd1GzjgDcrmsu9kx1dVw550bqa6u9ne0JEmSpA6r1c+kjIhzyK+uPLSw/UfgmpTS\nLW3UW7sxbNiwrFtoE67UkfR+k8tBz55ZdyFJkiRJ2lGtCikjYjJwEfCfwBOF4ROAqRFRkVK6rI36\naxc+9zk49NCsu9gxrtSRJEmSJElSe9XalZTnAaNSSrOKxv4rIp4hH1y2OKSMiO8Afw8cDtQBjwMT\nUkovFNV0Aq4Gvgp0AuYDY1NKa4pqegE3ACcBbwG3AxenlDYX1ZwEXAX0B6qA76eUft5cj927u1JH\nkiRJkiRJ2lnKWnncXsD/NDK+mO0PPv+WfLB5HPD5wty/jojORTXXAKcBZwInAj2BOfU7I6IMmFc4\n9/HAcGAEMLmopjdwL7AAOBr4MXBLRJyynf1KkiRJkiRJakOtXUn5C/KrKS8qGR8NzNieiVJKQ4q3\nI2IEsAYYCPwuIroCI4GvpZQeLdT8C7AsIj6VUloEfIH8SszPpZSqgaURMRGYEhGTUkrvFfr9U0rp\n24VTPR8RnyH/XM0Ht6dnSZIkSZIkSW2nxSFlRFxdtJmAcyNiMPDfhbHjgV7kb7PeEfsW5n+zsD2w\n0OeCLSdP6fmIqCL/HMxFhXMvLQSU9eYD15O/tfvpQs1DJeeaD0zdwX4lSZIkSZIk7YDtWUl5TMn2\n4sL3PoXvbxS++re2mYgI8rd2/y6l9IfC8AHAOyml9SXlrxf21de83sj++n1Pb6Oma0R0Sin9pbV9\nS5IkSZIkSf+/vTuPsqyq7wX+/SkqgwGMLY28gDO2URpjq2EQh4eBFzUmDitaSjBREzWKCs8hec9E\njHl5UaNGcQgRVFCoF4Mm6jMBA1FJhASkUUEGo0FaQFpKm26UJkw7f5xTevumhyq6uk9V1+ez1l11\n79n7nPs7vXqdPv2tvc/mrptxSNlae8q2LKT3gSQ/n+QJM+hb6UZcbsnm+tQM+uS885JLL91w2wEH\ndC8AAAAAIJmcnMzk5OQG29auXTujfe/qMynnXFW9L8nTkhzWWrtupOn6JPesqt3HRlPulZ+OjLw+\nyePGDrl0pG3659KxPnslWddau3VztR1ySLJ8+czOAwAAAAAWo4mJiUxMTGywbeXKlVmxYsUW972r\nq3vPqT6g/NV0C9+sGmu+KMntSQ4f6b9/kv2SnNdvOj/JAVW1ZGS/I5KsTXL5SJ/Ds6Ej+u0AAAAA\nwEAGH0lZVR9IMpHkmUl+XFXTox3XttZuaa2tq6qTk7yrqtYkuSnJe5N8ubV2Yd/380kuS/Kxqnpj\nkvsneWuS97XWbuv7/EWSV1XV25J8OF1g+dx0ozcBAAAAgIHMh5GUL0+ye5IvJrlu5PXrI32OTfL/\nk5wx0u85042ttTuTPCPJHelGV56a5KNJ3jzS5ztJnp7kqUm+2h/zJa218RW/AQAAAIDtaPCRlK21\nLQal/crbx/SvTfX5brqgcnPH+VKSLU+CBwAAAAC2m/kwkhIAAAAAWMSElAAAAADAoISUAAAAAMCg\nhJQAAAAAwKCElAAAAADAoISUAAAAAMCghJQAAAAAwKCElAAAAADAoISUAAAAAMCghJQAAAAAwKCE\nlAAAAADAoISUAAAAAMCghJQAAAAAwKCElAAAAADAoISUAAAAAMCghJQAAAAAwKCElAAAAADAoISU\nAAAAAMCghJQAAAAAwKCElAAAAADAoISUAAAAAMCghJQAAAAAwKCElAAAAADAoISUAAAAAMCghJQA\nAAAAwKCElAAAAADAoISUAAAAAMCghJQAAAAAwKCElAAAAADAoISUAAAAAMCghJQAAAAAwKCElAAA\nAADAoOZFSFlVh1XVZ6rq2qq6s6qeOdb+kX776Ovvxvrcp6pOq6q1VbWmqk6qqt3G+iyvqnOran1V\nXV1Vr98e5wcAAAAAbNq8CCmT7Jbkq0lemaRtos/fJ1maZO/+NTHWfnqSRyQ5PMnTkzwxyYnTjVX1\nM0nOSnJVksckeX2S46vqpXN2FgAAAADArO00dAFJ0lo7M8mZSVJVtYlu/9Fau2FjDVW1LMmRSVa0\n1i7utx2T5HNV9brW2vVJjkpyjyQvaa3dnuTyqvqFJMclOWlOTwgAAAAAmLH5MpJyJp5cVaur6oqq\n+kBV/exI28FJ1kwHlL2z043K/MX+80FJzu0DymlnJXl4Ve2xTSsHAAAAADZpoYSUf5/k6CT/Pckb\nkjwpyd+NjLrcO8n3R3dord2R5Id923Sf1WPHXT3SBgAAAAAMYF5M996S1tonRj5+o6ouSfLtJE9O\n8oXN7FrZ9DMup9uzhT4577zk0ks33HbAAd0LAAAAAEgmJyczOTm5wba1a9fOaN8FEVKOa61dVVVT\nSR6aLqS8Psleo32q6u5J7tO3pf+5dOxQ0/uMj7DcwCGHJMuXb23VAAAAALDjmpiYyMTEhmtdr1y5\nMitWrNjivgtluvcGqurnktw3yff6Tecn2bNfCGfa4elGSl4w0ueJfXg57YgkV7bWZhbpAgAAAABz\nbl6ElFW1W1UdWFWP7jc9uP+8b9/29qr6xap6QFUdnuRvk3wz3cI3aa1d0b//UFU9rqoOTXJCksl+\nZe8kOT3JrUk+XFU/X1XPS/LqJO/cjqcKAAAAAIyZL9O9H5tu2nbrX9PB4SlJfjfJ8nQL5+yZ5Lp0\ngeQfttZuGznGC5K8L92q3ncmOSPJa6YbW2vrqurIvs9XkkwlOb61dvK2Oy0AAAAAYEvmRUjZWvtS\nNj+q83/M4Bg3JjlqC30uSbcyOAAAAAAwT8yL6d4AAAAAwOIlpAQAAAAABiWkBAAAAAAGJaQEAAAA\nAAYlpAQAAAAABiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAABiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAA\nBiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAABiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAABiWkBAAAAAAG\nJaQEAAAAAAYlpAQAAAAABiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAABiWkBAAAAAAGJaQEAAAAAAYl\npAQAAAAABiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAABiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAABjUv\nQsqqOqyqPlNV11bVnVX1zI30+aOquq6qbq6qf6iqh46136eqTquqtVW1pqpOqqrdxvosr6pzq2p9\nVV1dVa/f1ucGAAAAAGzevAgpk+yW5KtJXpmkjTdW1RuTvCrJy5I8PsmPk5xVVfcc6XZ6kkckOTzJ\n05M8McmJI8f4mSRnJbkqyWOSvD7J8VX10m1wPgAAAADADO00dAFJ0lo7M8mZSVJVtZEur0ny1tba\nZ/s+RydZneTXknyiqh6R5MgkK1prF/d9jknyuap6XWvt+iRHJblHkpe01m5PcnlV/UKS45KctE1P\nEAAAAADYpPkyknKTqupBSfZOcs70ttbauiT/muTgftNBSdZMB5S9s9ONyvzFkT7n9gHltLOSPLyq\n9thG5QMAAAAAWzDvQ8p0AWVLN3Jy1Oq+bbrP90cbW2t3JPnhWJ+NHSMjfQAAAACA7WxeTPe+iyob\neX7lLPtMTy3f7HHOOy+59NINtx1wQPcCAAAAAJLJyclMTk5usG3t2rUz2nchhJTXpwsTl2bDkZB7\nJbl4pM9eoztV1d2T3Kdvm+6zdOzY0/uMj7DcwCGHJMuXz7puAAAAAFg0JiYmMjExscG2lStXZsWK\nFVvcd95P926tXZUuYDx8eltV7Z7uWZPn9ZvOT7JnvxDOtMPThZsXjPR5Yh9eTjsiyZWttZlFugAA\nAADAnJsXIWVV7VZVB1bVo/tND+4/79t//vMkb6qqX6mqA5KcmuSaJJ9OktbaFekWwflQVT2uqg5N\nckKSyX5l7yQ5PcmtST5cVT9fVc9L8uok79wuJwkAAAAAbNR8me792CRfSPdsyJafBoenJHlxa+3t\nVbVrkhOT7Jnkn5L8cmvt1pFjvCDJ+9Kt6n1nkjOSvGa6sbW2rqqO7Pt8JclUkuNbaydvyxMDAAAA\nADZvXoSUrbUvZQujOltrxyc5fjPtNyY5agvHuCTJk2ZfIQAAAACwrcyL6d4AAAAAwOIlpAQAAAAA\nBiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAABiWkBAAAAAAGJaQEAAAAAAYlpAQAAAAABiWkBAAAAAAG\nJaQEAAAAAAYlpAQAAAAABrXT0AUAAAA7jlWrVmVqamroMubEkiVLst9++w1dBgAsCkJKAABgTqxa\ntSrLlj0869ffMnQpc2KXXXbOFVdcKagEgO1ASMl/sSP99jvxG3AAgO1lamoq69ffkmc/O1myZOhq\nts7UVPKpT92Sqakp95IAsB0IKdnAqlWrsuwRy7L+5vVDlzJndtl1l1xx+RVuLgEAtpMlS5J99hm6\nCgBgIRFSsoGpqamsv3l9jjrxqCzdf+nQ5Wy11d9cnY+/7ON+Aw4AAAAwjwkp2ail+y/NvgfuO3QZ\nAAAAACwCdxu6AAAAAABgcRNSAgAAAACDElICAAAAAIMSUgIAAAAAgxJSAgAAAACDElICAAAAAIMS\nUgIAAAAAgxJSAgAAAACDElICAAAAAIMSUgIAAAAAgxJSAgAAAACDElICAAAAAIMSUgIAAAAAgxJS\nAq+P6EYAABEHSURBVAAAAACDElICAAAAAINaECFlVb25qu4ce1020n6vqnp/VU1V1U1VdUZV7TV2\njH2r6nNV9eOqur6q3l5VC+L8AQAAAGBHttPQBczCpUkOT1L959tH2v48yS8neU6SdUnen+STSQ5L\nkj6M/Lsk1yU5KMk+ST6W5NYkb9oOtQMAAAAAm7CQQsrbW2s3jG+sqt2TvDjJ81trX+q3/VaSy6vq\n8a21C5IcmWRZkqe01qaSXFJVf5DkT6vq+Nba7ePHBQAAAAC2j4U03flhVXVtVX27qj5eVfv221ek\nC1vPme7YWrsyyaokB/ebDkpySR9QTjsryR5JHrntSwcAAAAANmWhhJT/kuQ3042IfHmSByU5t6p2\nS7J3kltba+vG9lndt6X/uXoj7RnpAwAAAAAMYEFM926tnTXy8dKquiDJ1Ul+Pcktm9itkrSZHH5L\nHc47L7n00g23HXBA9wIAAAAAksnJyUxOTm6wbe3atTPad0GElONaa2ur6ptJHprk7CT3rKrdx0ZT\n7pWfjpa8Psnjxg6ztP85PsLyvzjkkGT58q0sGgAAAAB2YBMTE5mYmNhg28qVK7NixYot7rtQpntv\noKruneQh6VbrvijdSt+Hj7Tvn2S/JOf1m85PckBVLRk5zBFJ1ia5bHvUDAAAAABs3IIYSVlV70jy\n2XRTvP9bkrekCyb/X2ttXVWdnORdVbUmyU1J3pvky621C/tDfD5dGPmxqnpjkvsneWuS97XWbtu+\nZwMAAAAAjFoQIWWSn0tyepL7JrkhyT8nOai19oO+/dgkdyQ5I8m9kpyZ5JXTO7fW7qyqZyT5YLrR\nlT9O8tEkb95O9QMAAAAAm7AgQsrW2sQW2v8jyTH9a1N9vpvkGXNcGgAAAACwlRbkMykBAAAAgB2H\nkBIAAAAAGJSQEgAAAAAYlJASAAAAABiUkBIAAAAAGJSQEgAAAAAYlJASAAAAABiUkBIAAAAAGJSQ\nEgAAAAAYlJASAAAAABiUkBIAAAAAGJSQEgAAAAAYlJASAAAAABiUkBIAAAAAGJSQEgAAAAAYlJAS\nAAAAABiUkBIAAAAAGJSQEgAAAAAYlJASAAAAABjUTkMXAAAAwMKzatWqTE1NDV3GnFmyZEn222+/\nocsAWLSElAAAAMzKqlWrsmzZsqxfv37oUubMLrvskiuuuEJQCTAQISUAAACzMjU1lfXr1+ew3/6d\n7LHPPkOXs9XWXndd/ulDf5mpqSkhJcBAhJQAAADcJXvss0/u+4AHDl0GADsAC+cAAAAAAIMSUgIA\nAAAAgzLdG3pWJwQAAAAYhpASYnVCAAAAgCEJKSFWJwQAAAAYkpASRlidEAAAAGD7s3AOAAAAADAo\nISUAAAAAMCghJQAAAAAwqEUXUlbVK6vqqqpaX1X/UlWPG7qmxeqiT140dAkAbITrM8D8NDk5OXQJ\nALDNLKqQsqqel+SdSd6c5BeSfC3JWVW1ZNDCFqmVn1w5dAkAbITrM8D8JKQEYEe2qELKJMcmObG1\ndmpr7YokL09yc5IXD1sWAAAAACxeiyakrKp7JFmR5Jzpba21luTsJAcPVRcAAAAALHaLJqRMsiTJ\n3ZOsHtu+Osne278cAAAAACBJdhq6gHmgkrRNtO2cJKtWbb9itpU1a7qfl19++Wb7TbdfdvZlWf3N\n8Tx3bt143Y35yl9/ZZt+xw9W/SDJzM/72q9/PWu/d902rWl7+NENU0m2fN6wUE3/3f63f0umpgYu\nZivN9Po82mdbX6Ndn7cd12d2dIv1+pwkN9xwQ6a28Ulfc801Oe2007bpdyTJkiVLcr/73W+zfVyf\nAZipkWvrzpvrV92M5x1fP9375iTPaa19ZmT7R5Ps0Vp71kb2eUGSbX8XAAAAAAA7the21k7fVOOi\nGUnZWrutqi5KcniSzyRJVVX/+b2b2O2sJC9M8p0kt2yHMgEAAABgR7Jzkgemy9k2adGMpEySqvr1\nJKckeVmSC9Kt9v3cJMtaazcMWRsAAAAALFaLZiRlkrTWPlFVS5L8UZKlSb6a5EgBJQAAAAAMZ1GN\npAQAAAAA5p+7DV0AAAAAALC4CSkBgBmrqquq6tVD1wGw0FXVF6rqXf377XptdS0HYD5aVM+kBAAA\nmIcem+THQxcBAEMykpIFr6qE7cCiVlX3GLqGmaqqu1VVDV0HwHzSWvtBa+2WoesAgCEJKZm1qnpu\nVX29qm6uqqmq+nxV7To6ZWWk799U1YdHPl9VVW+qqtOr6kdVdU1V/e7YPntU1UlV9f2qWltVZ1fV\n8pH2N1fVxVX1kqr69yRu6IBFpb/enlBV766qG5KcuaVrZ7/fm6pqdd/+oar6v1V18dhxN3sd30gt\nx/b/JvyoqlZV1furareR9hdV1Zqq+pWq+ka6a/a+c/aHAbAA9PfKp1bVTVV1bVUdN9a+wfTrGdwP\nL6+qf6yqdX37hVX1mJH2J1TVuf39+tVV9Z6q2nX7nC0A3DVCSmalqvZOcnqSk5IsS/KkJJ9KMptR\nMa9LcnGSRyf50yTvqarDR9rPSHLfJEcmeUySlUnOrqo9R/o8NMmzkzyrPw7AYnN0kv9IckiSlyf5\n62zm2llVL0zyv5K8PsmKJKuSvCJJ28o67khyTJJH9jU9JcnbxvrsmuQNSV7S9/v+Vn4nwELzZ0kO\nS/IrSY5I8uR01+JN2dj98Dkj98OnJfluf4zHpLunvi1JquohSf4+3b8Lj0ryvCSHJjlhLk8IAOaa\nabLM1v2T3D3J37TWvttv+0aSzGL23pdba+/o37+vqg5Ncmy6G68npHsmz16ttdv6Pm+oqmcleW66\ncDRJ7pHkN1prP9yqswFYuL7VWvu9JOmvo4/L5q+dr0ryodbaqX37W6vqiCS7ZSu01t478vHqqvqD\nJB/sv2/aTkle0Vq7dGu+C2Ah6keXvzjJC1prX+y3vSjJNZvoP5P74f2SvL219m99+7dHDvF7ST7e\nWpsOJf+9ql6b5ItV9YrW2q1zd3YAMHeMpGS2vpbknCSXVtUnquqlYyMcZ+L8jXx+RP9+eZKfSfLD\nfjrMTVV1U5IHJnnIyD5XCyiBRe4rI+8PzKavnQ/u+zw8yYVjx7hga4uoqqf20xCvqap1ST6W5L5V\ntctIt1sFlMAi9pB0v2D/yTW3tbYmyZWb6D+T++F3JTm5qv6hqt5YVQ8e2f/AJL85tu+ZfduD5uys\nAGCOGUnJrLTW7kxyRFUdnG6qyjFJ/riqDkpyZ/7rtO+ZLuYwPd3w3kmuSzeNfPxYN468t/ohsNiN\nXgdneu0cn9o93ndW1/GqekCSzyZ5f7qp5D9MN53xpH6/9X3X9Rs9AMDiMH1dnenjNbZ4TW+tvaWq\nTkvy9CRPS/KWqnpea+3T/f4nJnnPRvZfNfvyAWD7EFJyl7TWzk9yflW9NcnVSX4tyQ3ppoMn6VZw\nTfccnH8c2/2gjXy+on+/MsneSe5orbmJApiZmVw7r0zy+HTPMZv22LE+M72OT1uR5G6ttdeN7PP8\n2ZUOsMP7VpLb093zfjJJquo+SfZP8sWN9J/R/XBr7Vvpgsj3VNXpSX4ryaf7/R/ZWrtqDs8BALY5\n072Zlap6fFX9flWtqKp9kzwnyZIkl6f7T+zTq+ppVfXwdM8k29hU8EOr6nVV9bCqemW6Z+v8eZK0\n1s5ON/37b6vql6rqAVV1SFX98eiKhQD81AyvnSckeWlVHV1VD62qN6WbUjg6smem1/Fp30qyU1W9\nuqoeVFW/keRlc36CAAtYa+3HSU5O8o6qekpVPSrJR9ItPLax/pu9plfVzlV1QlU9qar2G3ku8WX9\nId6W5OC+z4H9Nf9Xq8rCOQDMa0ZSMlvrkjwxyWuS7J5uFOVxrbWzqmqndP/hPSXdb4vfnY2Pvnln\nutE7xydZm+TY/mZs2tOS/J8kH05yvyTXJzk3yeptcD4AC9HGpgxu9trZWju9qh6U5B1Jdk7yiSQf\nTfcf22kfzpav4z/57tba16vquHQrd/9J/32/l+TUADDq9ekWKvtMkpvS3Q/vnp9eU8ev65u7pt+R\nbuXvU5IsTTKVboTm8UnSWrukqp7U739uuinf307yVyPHn+nUcwDYbqo1/z6x/VTVVUnePbYaLAAD\nqKrPJ/lea+1FQ9cCAAAsbkZSAsAi0K+2/fIkZ6VbIGciyeFJnjpkXQAAAImQku3P0F2AYbR00wf/\nd5J7pVtI59mttS8MWhUAAEBM9wYAAAAABmZ1bwAAAABgUEJKAAAAAGBQQkoAAAAAYFBCSgAAAABg\nUEJKAAAAAGBQQkoAAAAAYFBCSgAAFryq+khVfWrk8xeq6l1D1gQAwMztNHQBAACwDTwryW1DFwEA\nwMwIKQEA2OG01m4cugYAAGbOdG8AALZaVd27qk6rqh9V1bVV9drRKddV9cKqurCq1lXV9/q+9xvZ\nf89+2/er6uaqurKqXjTS/qiqOqdvm6qqE6tqt83Us8F076q6qqp+v6pO7mu4uqp+e2yfn6uqv6qq\nNf13/G1VPWCk/clV9a/9Oa6pqn+qqn3n6s8QAGAxE1ICADAX3p3k4CTPSPJLSQ5L8piR9nskeVOS\n5Ul+NckDknx0pP2PkyxLcmT/8xVJppKkqnZJcmaSHyRZkeS5SZ6a5IRZ1nhckguTPDrJB5J8sKr2\n779jpyRnJVmb5ND+dVOSM6tqp6q6e5K/SfKFJI9KclCSv0zSZlkDAAAbYbo3AABbparuneToJM9v\nrX2x3/ZbSa6b7tNa++jILt+pqtcm+deq2rW1dnOSfZNc3Fq7uO+zaqT/UUl2TnJ0a+2WJJdX1auS\nfLaq3thau2GGpX6utfYX/fu3VdWxSZ6c5JtJnp+kWmu/M3JeL0mypu9zUZLd+2N8p+9y5Qy/FwCA\nLTCSEgCArfXgdL/8vnB6Q2ttXUZCvKpaUVWf6adZr0vyxb5pv/7nB5NMVNXFVfW2qjp45PjLknyt\nDyinfTndvezDZ1HnJWOfr0+yV/9+eZKHVdVN0690IzfvleQhrbU1SU5J8vn+PF5dVXvP4rsBANgM\nISUAAFur+p/jU58rSapq13TTtW9M8oIkj023+naS3DNJWmtnpgss353k/knOqaq3jxxnU9OqZzPd\neny175af3g/fO8lX0oWVB4689k9yel/ji9NN8/5ykuclubKqHj+L7wcAYBOElAAAbK1vJ7k9yU8C\nu6raPcnD+o/Lktw3ye+31r7cWvtmkqXjB2mt/aC1dmpr7egkr00yPfX6siSP7p9NOe0JSe5IN1V7\nLqzs672htfbvY6+bRmr8Wmvtba21Q5N8I13oCgDAVhJSAgCwVVprP0o3FfrP+hWwH5nk5HQhYkv3\nfMlbk7y6qh5UVc9Mt4jOT1TVW6rqmVX1kH7/Z6QLJ5PktCS3JDmlqh5ZVU9J8t4kp87ieZRbclq6\nhXo+XVVPqKoH9ufynqrap//8J1V1UFXtV1VHpAs1L9v8YQEAmAkhJQAAc+HYJOcl+WySzyf55yRX\nJLmltTaV5EXpVuX+RpI3JPmfY/vfmuRPknwt3fMqb08ykSSttfXpVv3+2SQXJPlEkn9Icsxm6hmf\nBr6xaeE/2dZ/xxPTBaqfTBc+fijdMynXJbk53YjQM9I9a/MvkpzQWvvLzdQAAMAMVWuzeYwPAABs\nWf8cymuTHNda+8jQ9QAAML/tNHQBAAAsfFX16HQjDS9IsmeSP0w3UvHTQ9YFAMDCIKQEAGCuvC7d\nati3JrkoyRNaaz8ctiQAABYC070BAAAAgEFZOAcAAAAAGJSQEgAAAAAYlJASAAAAABiUkBIAAAAA\nGJSQEgAAAAAYlJASAAAAABiUkBIAAAAAGJSQEgAAAAAY1H8CB4hSuaouk6UAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fdba6f75110>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "f, barplot = plt.subplots(1, figsize=(16,5))\n",
    "\n",
    "bar_width = 0.1\n",
    "offset = 0.12\n",
    "rho = 0.7\n",
    "\n",
    "# position of left-bar boundaries\n",
    "bar_l = [o for o in range_oil]\n",
    "\n",
    "mbar_w = 3*bar_width+2*max(0, offset-bar_width)\n",
    "\n",
    "tick_pos = [b*rho + mbar_w/2.0 for b in bar_l]\n",
    "\n",
    "colors = ['olive', 'lightgreen', 'cadetblue']\n",
    "\n",
    "for i in range_oil:\n",
    "    barplot.bar([b*rho + (i*offset) for b in bar_l], \n",
    "                blend_values[i], width=bar_width, color=colors[i], label=oil_names[i])\n",
    "\n",
    "plt.xticks(tick_pos, gas_names)\n",
    "barplot.set_xlabel(\"gasolines\")\n",
    "barplot.set_ylabel(\"blend\")\n",
    "plt.legend(loc=\"upper right\")\n",
    "plt.title('Blend Repartition\\n')\n",
    "                                                                          \n",
    "\n",
    "# Set a buffer around the edge\n",
    "plt.xlim([0, max(tick_pos)+mbar_w +0.5])\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "Notice the missing bar for (crude2, diesel) which is expected since blend[crude2, diesel] is zero in the solution.\n",
    "\n",
    "We can check the solution value of blends for *crude2* and *diesel*, remembering that crude2 has offset 1 and diesel has offset 2.\n",
    "Note how the decision variable is automatically converted to a float here. This would raise an exception if called before submitting a solve, as no solution value would be present."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "collapsed": false,
    "deletable": true,
    "editable": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "* value of blend[crude2, diesel] is 0\n"
     ]
    }
   ],
   "source": [
    "print(\"* value of blend[crude2, diesel] is %g\" % blends[1,2])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "## Summary\n",
    "\n",
    "\n",
    "You learned how to set up and use IBM Decision Optimization CPLEX Modeling for Python to formulate a Mathematical Programming model and solve it with IBM Decision Optimization on Cloud."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true,
    "deletable": true,
    "editable": true
   },
   "source": [
    "## References\n",
    "* [CPLEX Modeling for Python documentation](https://rawgit.com/IBMDecisionOptimization/docplex-doc/master/docs/index.html)\n",
    "* [Decision Optimization on Cloud](https://developer.ibm.com/docloud/)\n",
    "* Need help with DOcplex or to report a bug? Please go [here](https://developer.ibm.com/answers/smartspace/docloud).\n",
    "* Contact us at dofeedback@wwpdl.vnet.ibm.com."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<div class=\"alert alert-block alert-info\"> Note: To save resources and get the best performance please use the code below to stop the kernel before exiting your notebook.</div>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true,
    "deletable": true,
    "editable": true
   },
   "outputs": [],
   "source": [
    "%%javascript\n",
    "Jupyter.notebook.session.delete();"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "deletable": true,
    "editable": true
   },
   "source": [
    "<hr>\n",
    "Copyright &copy; IBM Corp. 2017. Released as licensed Sample Materials."
   ]
  }
 ],
 "metadata": {
  "celltoolbar": "Dashboard",
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}