Linear_programming-Analytics-Python.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# >> Linear programming:\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One of the Optimization topics is Linear Programming. In this category of optimization problems, both the cost function and all the restrictions are linear.\n",
    "\n",
    "Linear programming (LP) is one of the simplest ways to perform optimization. We can solve some very complex optimization problems by making a few simplifying assumptions. As an analyst one is bound to come across applications and problems to be solved by Linear Programming.\n",
    "\n",
    "It is a simple technique where we depict complex relationships through linear functions and then find the optimum points. The real relationships might be much more complex – but we can simplify them to linear relationships.\n",
    "\n",
    "It could be a breakthrough method while solving optimization problems. One can use it to find Minimum cost expenditure, maximum profit gained, shortest path to travel, etc.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Simple Example:\n",
    "> Here I am providing a simple example to illustrate solving a simple LP optimization problem.\n",
    "\n",
    ">> Find the maximal and minimal value of **z = 3x + 4y** subject to the following **constraints**:\n",
    "\n",
    ">>> x + 2y <= 14,\n",
    "\n",
    ">>> 3x - y >= 0, \n",
    "\n",
    ">>> x - y <= 2\n",
    "\n",
    ">> Here we need to find **optimum (Max/Min)** solution for the **Objective Function: z = 3x + 4y**. First of all, we will convert constraints equations into more simpler way:\n",
    "![alt txt](eqn.png \"Simplified equation\")\n",
    "\n",
    "\n",
    ">> Now we will plot graph for these constraint equations.\n",
    "\n",
    "![alt txt](graph.png \"Graph system\")\n",
    "\n",
    ">> The soltuion lies in the bounded (feasible) region. And the optimal solution is  @ either of these corners.\n",
    "\n",
    ">> Its tabular repsentation to find the corner points.\n",
    "\n",
    "![alt txt](Table.png \"Table\")\n",
    "\n",
    ">>> Solution for these corner points:\n",
    " \n",
    ">>>> (2, 6):      z = 3(2)   + 4(6)   =   6 + 24 =   30 \n",
    "\n",
    ">>>> (6, 4):      z = 3(6)   + 4(4)   = 18 + 16 =   34 \n",
    "\n",
    ">>>> (–1, –3):  z = 3(–1) + 4(–3) = –3 – 12 = –15\n",
    "\n",
    ">>>>> Then the maximum of **z = 34** occurs at (6, 4), \n",
    "\n",
    ">>>>> and the minimum of **z = –15** occurs at (–1, –3)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# >> Common terminologies used in Linear Programming\n",
    "\n",
    "> **Decision Variables**: The decision variables are the variables which will decide my output. They represent my ultimate solution. To solve any problem, we first need to identify the decision variables. For the above example, the total number of units for A and B denoted by X & Y respectively are my decision variables.\n",
    "\n",
    "> **Objective Function**: It is defined as the objective of making decisions. In the above example, the company wishes to increase the total profit represented by Z. So, profit is my objective function.\n",
    "\n",
    "> **Constraints**: The constraints are the restrictions or limitations on the decision variables. They usually limit the value of the decision variables. In the above example, the limit on the availability of resources Milk and Choco are my constraints.\n",
    "\n",
    "> **Non-negativity restriction**: For all linear programs, the decision variables should always take non-negative values. Which means the values for decision variables should be greater than or equal to 0."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# >> Solving LP problems:"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "> There are many ways to solve a LP problem; like:- solving it throgh converting it to Standard Form, Graphical Method.\n",
    "Out of which, Graphical Method is most popular method.\n",
    "\n",
    "> Additionally, we have multiple libraries for python that can help us soving a LP problem.\n",
    "viz.\n",
    "\n",
    ">> 1.) scipy.optimize.minimize\n",
    "\n",
    ">> 2.) scipy.optimize.linprog\n",
    "\n",
    ">> 3.) PuLP\n",
    "\n",
    ">> 4.) Analyzing graphs through plotting with Matplotlib\n",
    "\n",
    ">> ..... And there must be few other methods also."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "# >> Circumstances in which we apply LP technique to solve our daily analytics operation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "we use LP technique to solve the optimiziation cases:\n",
    "  >> When we want to get best outcome (such as **maximum profit** or **lowest cost**) in a mathematical model; considering whose **requirements are represented by \"linear relationships\" **.\n",
    "  \n",
    "  >> When its **objective function** can be represented through a ** <u>linear function </u>**.\n",
    "    and subject to subject to **<u>linear equality</u>** and **<u>linear inequality</u> constraints**.\n",
    "    \n",
    "  >> ->It can be applied in following application indaily analytics operations:\n",
    "  \n",
    "  >>> 1. Manufacturing industries use linear programming for <u>analyzing their supply chain operations</u>. Their motive is to **maximize efficiency** with **minimum operation cost**.\n",
    "  \n",
    "  >>> 2. Linear programming is also used in **organized retail** for <u>shelf space optimization</u>. Since the number of products in the market have increased in leaps and bounds, it is important to understand what does the customer want.\n",
    "  \n",
    "  >>> 3. Optimization is also used for <u>optimizing Delivery Routes</u>. This is an extension of the popular traveling salesman problem.\n",
    "   >>> 4. Optimizations is also used in <u>Machine Learning</u>. **Supervised Learning** works on the fundamental of linear programming.\n",
    "   >>> 5. There are many more applications of linear programming in **real-world** like applied by **Shareholders, Sports, Stock Markets**, etc.  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Question:"
   ]
  },
  {
   "cell_type": "raw",
   "metadata": {},
   "source": [
    "> A car company produces 2 models, model A and model B. Long-term projections indicate an expected demand of at least 100 model A cars and 80 model B cars each day. Because of limitations on production capacity, no more than 200 model A cars and 170 model B cars can be made daily. To satisfy a shipping contract, a total of at least 200 cars much be shipped each day. If each model A car sold results in a $2000 loss, but each model B car produces a $5000 profit, how many of each type should be made daily to maximize net profits?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    ">  car models - **A** & **B**\n",
    "\n",
    "> *demand: A >= 100 ; B >= 80\n",
    "\n",
    "> *produce: A <= 200 ; B <= 170\n",
    "\n",
    ">> So, we can conclude that: **200 >= A >= 100**   and  **170 >= B >= 80**\n",
    "\n",
    "> *Shipping: A + B >= 200\n",
    "\n",
    "> *objective function: **( z = 5000B - 2000A )**\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# visualization of the problem:\n",
    "> # The solution lies in the feasible region( green region) ------"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEWCAYAAAAzcgPFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd4FOX2wPHvoVcDIsYIKJEgEBQCRMACxCsoNlCxYEH4\ngaKoVwG9iO2KvXdEBaWoCGIDK9iIV70CFzCCNImASkSQskBESsL5/fHOxiWmk2R2yPk8zz7ZndmZ\nOfNmds6Ud95XVBVjjDEm2lTyOwBjjDEmL5agjDHGRCVLUMYYY6KSJShjjDFRyRKUMcaYqGQJyhhj\nTFSKygQlIhNF5N4CxmeKyFFlsNyymu8QEVnvzb9Bac+/gOXeKiIvltfyTMVU2O9GRNaISPdizO9e\nEdkoIr8V4bs58xaRUSLyalGXU8RYymSfEER+7E+iMkEVRlXrqOqq/ZmHiKSKyBWlPd88llMVeBw4\n1Zv/ptKcf8RyUkRkbeQwVb1fVa/Ib5qg8NZtr7ezCL/6R4yvLiLjRWSbiPwmIsNzTZ8kIgtEZIf3\nNynX+GHedNu8+VSPGLdGRP70lrlFRD4QkSb7sS79vRi2ichaEXlYRKpEjD9YRN4RkT9E5CcRuSTX\n9KeIyHJvXWaLyJER40REHhKRTd7rIRGRXNOLiKwSkaUlXYfcIn83hR1cFkZEjgBuBBJV9bDSirGk\nymKf4Ie89nfFVZz9SWkdLAQyQQVMLFADWOJ3IAH3q7ezCL8mRYwbBTQHjgROBkaISE8AEakGzABe\nBeoDk4AZ3nBE5DRgJHCKN/1RwF25ln22qtYB4oD1wDP7sR61gKHAIUAnb7k3RYx/FtiN224uBZ4T\nkdZerIcAbwN3AAcD84HXI6YdDJwDtAXaAGcDV+VaflfgUOAoETluP9ajrBwBbFLVDX4HYv4SeRBV\nrlTVlxfQCkgFQridd6+IcROB54FPgO3AF8CREeMVSPDeVwceBX7G7TyeB2pGfLc3kAZsA34EegL3\nAdnATiATGB05X9yO4zegcsR8zgUWee8r4XZqPwKbgGnAwXms49HAH958M4HPgabe5yoR30sFrvDe\nDwC+8tZpC7AaOD3iuwcDE4BfvfHTgdrAn8BebzmZwOG4HferEdP28so65C2zVcS4Nbgd5SJgK27H\nV6OI/8vvIpab6a1fSiluKynA2gLG/4o7Qw1/vhuY6r0/FcgAJGL8z0BP7/1rwP0R4/4B/JarXLpH\nfD4D+KEU12048J73vjYuOR0dMf5l4EHv/WDgvxHjwv/3lt7n/wKDI8YPBObkWt54YDIu0Y0uIK7/\nC8flfV4JvBHx+RcgKdfvZjCwx1uHzIj1KtK2BXTPtR1PLOJ22917X6TtvSTrFrFfehb4ALdfmgs0\ni5juVGCFt45jcPutK/Ip38rArbh9yHZgAdDEG3cC8D9vPv8DTsi1r7gH+Nqb7mPgEG9cDdyB2CZv\nnf+HO9DJb3/3lLeu27zld4lYzijgTW9+24ArIsuXv/Zj/XG/p43Abd64nt42sMdb3ncR+7ZVXtyr\ngUsL/X2U1g+tmD/KqkC69w+qhtspbAdaRGwI23FHe9W9gvwqYvrIjeYJ4F3cjrsu8B7wgDeuo/dP\n7oFLKo3468ecmnvjyTXfH4EeEePeAEZ6728A5gCNvfheAKbks67hf2SVvD7njsX7J+4BrsRtxENw\nO2Dxxn+A+4HX98qxmzc8hVw78VwbVDhZ9vCmG+H9D6pF/NDn4RLbwcAy4OqIeYWAk4rwvx0MLAcO\nymPcEd588ntdks88U3Ab/Hrchv0EUNsbV98rz9iI7/cBFnvvhwEf5Zrfe8CN3vvvgIsixjXw5tcg\njx1gLdwZ2Mul+FuYzl8JqB2wI9f4G/lrR/8U8Fyu8YuBPt77rUCniHEdgO0Rn2vhdjZneGW0Mfz/\nzyOuo7z/SSVvm/gpvH1547YAlfL43UwE7s01rwK3rTz+12sjPhdlu/1bgipouv1ct024/UoVXKIP\nHwgd4pXted64G3C/4/wS1L+8/10LQHBnvQ288tkC9PPmc7H3Obw9puL2TUcDNb3P4e3nKty2XQu3\n7+iA9zsk7/3dZd4yq+C2s9/wDhy8styDOyOv5C0rsnybemUzzhvXFtjFXwcBOd/1Ptf2yie8j48D\nWhf2+/DrEl9noA6uYHer6ufA+7h/RtgHqvofVd0F3AYcn/vav3d9fTAwTFU3q+p24H6gr/eVQcB4\nVf1EVfeqaoaqLi9ijFPC8YhIXdyPeoo37mrc0cJaL75RwPmleBr8k6qOU9Vs3A4xDogVkTjgdNyP\ne4uq7lHVL4o4z4twZfqJqu7BnaHVxB2thT2tqr+q6mbchp5zr0ZV66nqVwUtQEROAu7FnQ1vyz1e\nVX/25pPf67V8Zr3ciyUOdzDTAXdfD9x2BG7nHLYNd7ASHh85rrDx4bjrRgybLiIh/jrYeSSfOItF\nRAYCybj/RTiW3OW2v+tSJ+I+1Hm4ncjHuAOdqsCZecWm7r7Ldly5dwVmAb+KSEugG/Clqu4t0oo6\n+W5bhSjKdlus6fZz3d5R1XmqmoVLUOH1OANYoqpve+Oexu3w83MFcLuqrlDnO3X3p88EVqrqK6qa\npapTcNv/2RHTTlDVH1T1T9zVm3AMe3AJJ0FVs1V1QV6/wzBVfVVVN3nLeQx3sN0i4ivfqOp0b9/5\nZz6zuUtV/1TV73AHe20LWOe9wDEiUlNV16lqobc9/EpQhwO/5NoIfsKd4YT9En6jqpnAZm+6SA1x\nRwsLRCTk7URmesMBmuCONkriNeA874b5ecBCVf3JG3ck8E7EMpfhTqFjS7is3HI2bFXd4b2tg1uf\nzaq6pQTzDB8phue7F1fGkWUe+YPawV87/0J5Bw/TgP6q+kMJ4suXqv6mqku9H8pq3NFwH290pvf3\noIhJYnA7oPD4yHGFjY/x/m6PGHaOqtbDXUK5DvhCRP52A19EukRU4ijwxyci5wAP4C7fbixhrEVZ\nl0z1Dllxl2OmeTukncBb3rD8fIE7o+nqvU/F7cC7eZ+Lo6TbVlG225JMV9J1y289DmfffZYC+1Ra\nyiW/fdM+cXty7xvzi+EVXLKdKiK/ehVwquYXgIjcJCLLRGSrtx+LwZ0Jhv2Sz6SRivR/VdU/cAcN\nVwPrvMpGLQubuV8J6legiYhELv8I3L2CsJyzJRGpgzv1/TXXfDbirlm3jjgKj1F3QxtcATfLJwbN\nZ7gbqboUt2GcDlyCS1hhv+B2LJFH/zVUNSOveeXyh/e3VsSwotZW+gU4WETq5RVyIdP+ikusQM7Z\nZxP2LfMSEZGauEtVT6rqRwV874hcNfFyvy4t4iIVb9v1kvU69j1ya8tflVKWAG1y1WZrk2t87mnX\nax61Lb2j0rdxByMn5TH+S/2rEkfr/IL3KnCMw1W+WBwx6gegiog0L2BdcmIVkdq47bugdVnifbcx\n7uzzMq/G4m/A+cAZXuWLvIR34l28919Q+E68sO2wuEq63RY2XUnWrSDrcJf8I5fXOP+v57tv2idu\nT+59Y568Kyp3qWoi7gzzLODy8OjI74pIF9yB3oVAfe8AbCvuciN5TVNMf5tWVWepag/clZDluN9A\ngfxKUHNx2XaEiFQVkRTcKezUiO+cISInebWt7sHd7N0no3tHReOAJ0TkUAARaeTVzAJ4Cfg/r2pu\nJW9cOGuvx11vLshruGvJXXH3oMKeB+4LV/EVkYYi0rsoK66qv+M2tstEpLJ3mSe/JJp72nXAR8AY\nEanvlV3XiPVpICIx+Uw+DTjTK4uquGvOu3A31vfXeGC5qj5cSPw/67418XK/Juc1nYicLCJHelWk\nmwAP4Wrmhb0M3O6VSSvc/buJ3rhUXEK5Xlx19OtxP57PI6YdJCKJIlIfV0NuInnwlt8bd99rWSFl\nkicR+Qfu0lAfVZ0XOc47ynwbuFtEanuXTHvhjowB3sFdIukjIjWAO3E3oMOXrV8GhnvbeSPc/zi8\nLv1wCbAF7pJQEu4+xlr2vbQe6QtcrciaqroW+BJ3A7wB8G0+0xTld1UcJd1uC5uuJOtWkA+AY0Xk\nHO9S/7UUfOD5InCPiDT3tqs24p6R/BA4WkQuEZEqInIRkIi7BVIg73dyrIhUxl3e3YO7rAZ//7/U\nBbKA33EHRf/m72fn+2M90DR8EiIisSLS2zuo2oU72y/0ErEvCUpVd+MS0um4s6AxwOW67/2h13A/\nwM24ew6X5TO7m3E3P+eIyDbgU7zrqN4O4P9wN9W34tUG9KZ7CnffaIuIPJ3PvKfgjqg+j7gME572\nXeBjEdmOqzDRqWhrD7gd6L9wN1xbU7wk0Q+34S0HNuCqLOOV3RRglXfpcZ/Loaq6AleGz+DK/Gzc\nEfzuoizUO8Ppks/ovsC5uc6G8vtuSbTDldEf3t9FwPUR4+/EXS75CZeQHlbVmZCzrZ2DO5IM4Sqh\nnBNeb+97DwOzvelXe/OL9J6IZOJ+9PfhLmOW9LGBO3CXUj6MKKvIs85rcPdKNuB+A0PCy/IObvp4\nMWzB3azvGzHtC7j7O4u91/veMHCX8sZ4l0tzXriDrTwv83mXajNxO2+8+xmrgK/V3R/Ny0tAorcN\nTi9qoeSnpNttYdOVcN0KWt5G4ALctrQJl1Tm43bGeXkcl0Q/xm1XL+GS5Sbcmc+N3nxGAGfl2v/k\n5zBczbttuAOoL/jr4Cb3/m4W7nbID7jtfidFu6RXVOED+k0ishCXa4bjzhA34/arQwqbSbhmWGB4\nGTkbV+38Z7/jMcaY3Lz91FpcVerZfscTVEF8UPcYXLYvtBkUY4wpLyJymojUE1ex6lbc/Zw5PocV\naIFKUCLSB3cp5uaiXpoyxphycjzuUnP4kuI5BVTPNkUQuEt8xhhjKoZAnUEZY4ypOPxpALAcHXLI\nIdq0adMSTfvHH39Qu3bt0g2onAQ5drD4/bJi0wqys7NJPDTR71BKLLBlvwJX9omV/Q4FgAULFmxU\n1YaFf7PsHPAJqmnTpsyfP79E06amppKSklK6AZWTIMcOFr9fUiamEAqFmD+0ZL+ZaBDYsk/Blf38\nvJ7DL38ikrtFi3Jnl/iMMcZEJUtQxhhjopIlKGOMMVHJEpQxxpioZAnKGGNMVPItQYlIDRGZJyLf\nicgSEbnLG36wiHwiIiu9v/UjprlFRNJFZEVEi+XGGGMOQH6eQe0C/qGqbXFN//cUkc7ASOAzVW0O\nfOZ9RkQScS03t8Y1iz/Ga1beGGPMAci3BKVOuDfUqt5Lgd64bs7x/p7jve8NTFXVXV6vqum47gbK\nxDPPwNdfN8BagjLGGH/4+qCudwa0AEgAnlXVuSIS63XMB67F8nA36o3Yt2XgteTT7bOIDAYGA8TG\nxpKamlqsuLKz4cknk1m16lgmT97GoEGr6NAhVKx5+C0zM7PY6x1NLH5/hEIhsrOzAxl7WHDLPinw\nZV/afE1QXsdgSeK6MH9HRI7JNV5FpNjnMKo6FhgLkJycrCV5qnz5crj11hW8/noLbropiZNPhnvv\nhRNOKPasfBHUp+nDLH5/1FtTj1AoFMjYwwJb9vUIfNmXtqioxaeqIVw3Gj2B9SISB+D93eB9LQNo\nEjFZY29YmahaFc48cx0//ABPPQVLlsCJJ8JZZ0FaWlkt1RhjTJiftfgaemdOiEhNoAeuG/N3+asL\n6v7ADO/9u0BfEakuIvFAc2BeWcdZowZcfz2sWgUPPAD//S+0awcXXujOsowxxpQNP8+g4oDZIrII\n+B/wiaq+DzwI9BCRlUB37zOqugSYBiwFZgLXepcIy0Xt2jBypEtUt98OH30ErVvDgAGwenV5RWGM\nMRWHn7X4FqlqO1Vto6rHqOrd3vBNqnqKqjZX1e6qujlimvtUtZmqtlDVj/yIu149uOcel6iGDoWp\nU6FFC7jmGvj1Vz8iMsaYA1NU3IMKooYN4bHH4McfYdAgGDcOmjWDm26CjRv9js4YY4LPEtR+atQI\nnnvOdTZ24YXwxBMQHw///jds3ep3dMYYE1yWoErJUUfBpEnw/ffQs6e7DBgfDw8+CH/84Xd0xhgT\nPJagSlmrVvDGG7BwIRx/PNxyi7v09/TTsGuX39EZY0xwWIIqI+3awQcfwNdfu6R1ww3QvDm8+CLs\n2eN3dMYYE/0sQZWxE06Azz+HTz6BuDi48kpITITXXoO9e/2OzhhjopclqHIgAt27w5w5MGMG1KoF\nl14KbdvC9OlYg7TGGJMHS1DlSAR69YJvv4UpU2D3bjj3XOjYEWbNskRljDGRLEH5oFIl6NvXte/3\n0kuwYYOr+detG3z5pd/RGWNMdLAE5aMqVWDgQPjhB9f/1MqV0LWrS1YLFvgdnTHG+MsSVBSoXh2u\nu861SvHQQ/C//0FyMvTp486yjDGmIrIEFUVq1YIRI1w7f3fe6Wr+HXss9OvnkpcxxlQklqCiUEwM\njBrlWkm/6SZ46y1o2RKuugrWrvU7OmOMKR+WoKJYgwbw8MPu7Omqq2DCBEhIgGHDXMUKY4w5kFmC\nCoC4OBg92lWmuOQS12zSUUfBbbfBli1+R2eMMWXDElSANG0K48fD0qWu6/n773eJ6r77IDPT7+iM\nMaZ0+dnlexMRmS0iS0VkiYjc4A0fJSIZIpLmvc6ImOYWEUkXkRUicppfsfutRQvXUWJaGnTp4nr4\nPeoo19XHzp1+R2eMMaXDzzOoLOBGVU0EOgPXikiiN+4JVU3yXh8CeOP6Aq2BnsAYEansR+DRom1b\nePdd14RS27YwfLi7RzV2LGRlid/hGWPMfvGzy/d1qrrQe78dWAY0KmCS3sBUVd2lqquBdKBj2Uca\n/Tp1clXSP/8cjjjCVajo378jr74K2dl+R2eMMSVTxe8AAESkKdAOmAucCPxTRC4H5uPOsrbgktec\niMnWkk9CE5HBwGCA2NhYUlNTSxRXZmZmiaf1g4i7HzV37sGMG3ck/frV5Pbb/2DgwNV06bIRCdBJ\nVdDKPregxh8KhcjOzg5k7GHBLfukwJd9qVNVX19AHWABcJ73ORaojDu7uw8Y7w0fDVwWMd1LwPmF\nzb9Dhw5aUrNnzy7xtH777LPZOm2aasuWqqDaoYPqRx+p7t3rd2RFE+SyVw1u/N0mdNO2T7T1O4z9\nEtiy76batu0Wv8PIAcxXn/ODr7X4RKQq8BYwWVXfBlDV9aqarap7gXH8dRkvA2gSMXljb5jJQ6VK\ncMEFrgv6SZNg82Y4/XTX1t9//uN3dMYYUzg/a/EJ7ixomao+HjE8LuJr5wLfe+/fBfqKSHURiQea\nA/PKK96gqlwZLr8cli+H555zzSh16wanneba/DPGmGjl5xnUiUA/4B+5qpQ/LCKLRWQRcDIwDEBV\nlwDTgKXATOBaVbUqAEVUrRpcfTWkp8Ojj8LCha4fqnPPhcWL/Y7OGGP+zrdKEqr6FZDXbfsPC5jm\nPtx9KVNCNWvCjTfC4MHw5JPw2GOuivrFF8Ndd7lq6sYYEw2sJYkKqm5duOMOd8lv5EjX9XzLlnDl\nlfDzz35HZ4wxlqAqvIMPdk0mrVrl+qR6+WVo3hxuuAHWr/c7OmNMRWYJygAQG+su+aWnQ//+8Oyz\nrvmkW25xNQCNMaa8WYIy+2jSxDWVtGwZnHOO6+E3Ph7uvhu2b/c7OmNMRWIJyuSpeXOYPBkWLYJT\nTnE9/MbHu0oVf/7pd3TGmIrAEpQp0DHHwNtvw7x50KGD6+E3IcE9U7V7t9/RGWMOZJagTJEcdxzM\nmgVffOHuTV1zjev2Y9IkyMryOzpjzIHIEpQplnBTSR995LqkHzAAjj0Wpk2DvXv9js4YcyCxBGWK\nTQR69nRNJb31lmv376KL3CXADz4A15avMcbsH0tQpsRE4LzzXEWKV16BbdtcV/QnngizZ5d8vikp\nKQwdOrT0AjXGBJIlKLPfKleGyy5zDdK+8IJrieIf/4Du3V1vv8YYUxKWoEypqVrVtfGXng5PPOHO\nrI4/Hnr1gu++8zs6Y0zQWIIypa5GDRg61DWfdN998OWXkJQEffvCihV+R2eMCQpLUKbM1KkDt97q\nEtVtt8H770NiIgwcCD/95Hd0xphoZwnKlLn69eHee12iuuEGeO0111LFddfBunV+R2eMiVaWoEy5\nOfRQePxxd49q4EBXoaJZMxgxAjZt8js6Y0y08bPL9yYiMltElorIEhG5wRt+sIh8IiIrvb/1I6a5\nRUTSRWSFiJzmV+xm/zRuDM8/72r9nX++6+E3Ph5GjXJV1Y0xBvw9g8oCblTVRKAzcK2IJAIjgc9U\ntTnwmfcZb1xfoDXQExgjIpV9idyUimbNXP9T338Pp57qevSNj4eff+5LdnYNv8MzxvjMtwSlqutU\ndaH3fjuwDGgE9AYmeV+bBJzjve8NTFXVXaq6GkgHOpZv1KYsJCbCm2/C/PnQqROsXn01y5e/x+jR\nsGuX39EZY/xSxe8AAESkKdAOmAvEqmr41vlvQKz3vhEQ+djnWm9YXvMbDAwGiI2NJTU1tURxZWZm\nlnhavwU19hEj4IcfJrJu3bX8858NuffenfTrt4aePddTuXJw2lAKavmHQiGys7MDGXtYcMs+KfBl\nX+pU1dcXUAdYAJznfQ7lGr/F+zsauCxi+EvA+YXNv0OHDlpSs2fPLvG0fgty7N26ddM2bdrqxx+r\nHnecKqg2b6762muq2dl+R1c0QS3/bhO6adsn2vodxn4JbNl3U23bdovfYeQA5qvP+cHXWnwiUhV4\nC5isqm97g9eLSJw3Pg7Y4A3PAJpETN7YG2YOQCLQowfMnQvTp0P16nDJJe6B3xkzrEFaYyoCP2vx\nCe4saJmqPh4x6l2gv/e+PzAjYnhfEakuIvFAc2BeecVr/CECvXu7ppJeew127nRd0XfuDJ98YonK\nmAOZn2dQJwL9gH+ISJr3OgN4EOghIiuB7t5nVHUJMA1YCswErlXVbH9CN+WtUiW4+GJYuhRefNE9\n4HvqqXDyyfD1135HZ4wpC37W4vtKVUVV26hqkvf6UFU3qeopqtpcVbur6uaIae5T1Waq2kJVP/Ir\nduOfKlVg0CBYuRKeeso9S3XSSXDGGbBwod/RGWNKk7UkYQKpenW4/nr48Ud48EHXrUeHDu7B36VL\n/Y7OGFMaLEGZQKtdG26+GVavhn//G2bNcl3Q9+/v2v4zxgSXJShzQIiJcS1RrFoFw4bBtGnQogVc\nfTVkWF1PYwLJEpQ5oDRs6Nr2+/FHuPJKGD/eNal0443w++9+R2eMKQ5LUOaAdPjhMGaM6yCxb194\n8kk46ii44w4IhfyOzhhTFJagzAEtPh4mTnQN0p5+uuuXKj4eHngAMjP9js4YUxBLUKZCaNXK3Zda\nuBBOPNH19NusmauqvnOn39EZY/IieoA/ip9YN1End5hcomm3hrYSUy+mlCMqH0GO/bu078jKyqJD\ncocyW8a2ba7mXygE1avBkU3hsMNcyxWlIajlf2XSlWRlZTHh+wl+h1JigS37tASysrKZkLza71AA\naP9F+wWqmuxnDFHRmrkx5e2gg6BtWwhtcYnqhx/gl5+haVPX8y+llKiCSFHWhNb4HUaJ1cquFcj4\nd2Y1portkvdxwJdGrRa1aJfarkTTpqam0i6lZNP6LcixD0sZRigUIi01rVyWl6Lw/vuuAsV330Hr\nynDPPa7Nv5KeUQW1/OtMrEPG7xmc1/A8v0MpsUePfpSbfrjJ7zCKb8Jsjok5psT7q1IXBQdpdg/K\nVHgicPbZ7v7U669DVhacdx507Oge/D3Ar4IbE7UsQRnjqVQJLrzQ1fibMAE2boSePaFbN/jyS7+j\nM6bisQRlTC5VqsCAAe4ZqmefhfR06NrVJav58/2OzpiKwxKUMfmoVg2uucYlqEceccnpuOPc5b8l\nS/yOzpgDnyUoYwpRqxbcdJNr5++uu+Czz1yDtJdd5pKXMaZs+N3l+3gR2SAi30cMGyUiGbk6MQyP\nu0VE0kVkhYic5k/UpqI66CDXYvqqVTBiBLz9NrRsCYMHwy+/+B2dMQcev8+gJgI98xj+RGQnhgAi\nkgj0BVp704wRkcrlFqkxngYNXB9Uq1a5S4CTJkFCAgwdCuvX+x2dMQcOXxOUqv4H2FzoF53ewFRV\n3aWqq4F0oGOZBWdMIQ47DJ5+2j3k268fjB7tGqS99VbYvv2Af8TQmDIXrb+if4rI5cB84EZV3QI0\nAuZEfGetN+xvRGQwMBggNjaW1NTUEgWRmZlZ4mn9FuTYQ6EQ2dnZgYr/ssugW7eaTJzYlAceiKV2\n7U5cdNFq+vRZS61a2X6HV2ShUIhqUo1Hj37U71BKrHH1xoGMf0ytZuheArXdlzXf2+ITkabA+6p6\njPc5FtgIKHAPEKeqA0VkNDBHVV/1vvcS8JGqvlnQ/JOTk3V+CesGp6amkpKSUqJp/Rbk2FNSUlxL\nEmnl05JEaVu0CK65ZiNff30IhxwCt9wCQ4ZAzZp+R1a4lIkpZPyeQfqO4Nb+CHpLEou/O8TvSAAQ\nEd/b4vP7HtTfqOp6Vc1W1b3AOP66jJcBNIn4amNvmDFRpU0buPfe75k7F9q1c50lNm8Ozz8Pu3f7\nHZ0xwRF1CUpE4iI+nguEa/i9C/QVkeoiEg80B+aVd3zGFFXHjvDxx5Ca6hqhHTLE1fp7+WXIDs5V\nP2N843c18ynAN0ALEVkrIoOAh0VksYgsAk4GhgGo6hJgGrAUmAlcq6r2MzdRL9xU0ocfQr160L+/\ne47qzTdh716/ozMmevldi+9iVY1T1aqq2lhVX1LVfqp6rKq2UdVeqrou4vv3qWozVW2hqh/5Gbsx\nxSHievSdP98lJoALLoDkZJe4rEFaY/4u6i7xGXMgq1QJ+vSBxYvd81OhEJx5JnTpAl984Xd0xkQX\nS1DG+KByZbj8cli+HJ57znWamJICp54K8+zOqjGAJShjfFWtGlx9tWvT77HH4NtvoVMn11niokV+\nR2eMvyxBGRMFataE4cNd80n33ONq/iUlwcUXu5YqjKmILEEZE0Xq1oXbb3eJauRIePddSEyEK66A\nn3/2OzpjypclKGOi0MEHw/33u0R13XXwyivuYd/rr4fffvM7OmPKR7ESlIjEeK2KG2PKQWwsPPmk\nu0fVvz9HXG/+AAAgAElEQVSMGeMapB05EjZt8js6Y8pWoQlKRD4TkYNEpD6QBrwiIo+UfWjGmLAm\nTWDsWFfr77zz4OGHXaK66y7Yts3v6IwpG0U5gzpYVbcB5wGvqmoHwDoLNMYHCQnw6quuht8pp8Co\nUS5RPfII7Njhd3TGlK6iJKgqItIQuAB4r4zjMcYUwTHHuB59//c/1xrFiBEueY0ZYw3SmgNHURLU\nfcAXwM+qOk9EjgJWl21YxpiiSE6GmTNdKxQJCXDttXD00TBhAmRl+R2dMfun0ASlqlNVNVFVB3uf\nV6lq77IPzRhTVF27uiQ1cyYccggMHOjOsl5/3RqkNcGVb4ISkRu9v0+IyOO5X+UXojGmKETgtNPc\nZb+334YqVaBvX2jfHt57zxqkNcFT0BnUj97f74ElebyMMVFIBM49F777zlWoyMyEXr3ghBPg88/9\njs6Yoss3QanqdO/tq143GDkv4J3yCc8YU1KVK8Oll8KyZa6K+tq1rubfKafAN9/4HZ0xhStKJYm5\nInJc+IOI9MZ1MmiMCYCqVeHKK2HlSvfQ7/ffu7Ops86CtDS/ozMmf0VJUAOAF0TkARGZBPwT6FEa\nCxeR8SKyQUS+jxh2sIh8IiIrvb/1I8bdIiLpIrJCROxZLGOKoUYNuOEG+PFH14zS119Du3Zw4YXu\nAWBjok1RavGlAXfiEtOpwJWqWlrNVk4EeuYaNhL4TFWbA595n/GaWOoLtPamGSMilUspDmMqjDp1\n4JZbXB9Ut9/uevRt3RoGDICdO/2Ozpi/FKWpoxeAEUAScAXwkYhcVRoLV9X/AJtzDe4NTPLeTwLO\niRg+VVV3qepqIB3oWBpxGFMR1avnuvZYvRqGDoWpU11niRvW1YNtcX6HZwxVivCdlcDVqqpAuoj8\nF3gSeKGMYopV1XXe+9+AWO99I2BOxPfWesP+RkQGA4MBYmNjSU1NLVEgmZmZJZ7Wb0GOPRQKkZ2d\nHdj4IXjlf/bZ0LlzNYbM3cWmTbWpOvlnTjg1jX+cPY/aB/3pd3jF0rh6Yx49+lG/wyi2MbWaoXsJ\n1HZT1gpNUKr6aK7PW4D+ZRbRvstSESn20xuqOhYYC5CcnKwpKSklWn5qaiolndZvQY69Xr16hEKh\nwMYPwS3/Z/+ozoYN2zju4iq8+moHvvy0FXs7PQYnPAY1gtEq7aNHP8pNP9zkdxjFt6MDx8QcE8jt\npqwU5RJfMxGZKiKLROSH8KsMY1ovInHesuOADd7wDKBJxPcae8OMMaWoWrW9TJoEixcL551dC/7z\nbyo9/RN8eTPsruV3eKYCKUotvonABECA04FpwOtlGNO7/HWG1h+YETG8r4hUF5F4oDkwrwzjMKZC\nS0yEN94QFiyA07rFwGcPUumZNTDnn5BVze/wTAVQlARVS1VnAajqj6p6Oy5R7TcRmYJ7pqqFiKwV\nkUHAg0APEVkJdPc+o6pLcMlxKTATuFZVs0sjDmNM/tq3hw8/FL76Crq0bwAzn6bS6FWwYBBkW0Va\nU3aKkqB2iUgl4EcRuVpEzgbqlsbCVfViVY1T1aqq2thrqWKTqp6iqs1Vtbuqbo74/n2q2kxVW6jq\nR6URgzGmaE48EWbPrsQnn0D75nHw3otUGvMDLLoY9orf4ZkDUFES1DCgNnA9cCKuqvnAsgzKGBOd\nRKB7d5g3txIzZkCruCPh7deo9MJiWNYbrEFaU4qK8qDuXFXd7j2cO0xVe6vq1+UQmzEmSom4BmgX\nfVeZKVOgad2W8Pp0Kr34P0jvYYnKlIqinEFF+rhMojDGBFKlSq5LjxXLKvPSSxBXqR28+jEy6Qv4\n6US/wzMBV5QHdSMF70LzihVQwucKkkIh97h9AAU59ifT0sjKyirx/y0aBLb8k9KoU4Kyr4K77j+g\nKayrDmt+UvZMqAI12yD11yDVtpdBsHmLr/kcx/0ZvF3VDb9BrY2rIeV8v0OJGvkmKBH5ELhGVddE\nDB5f5hEZYwKrUiVo1AgOixMyMuDnn+uT9Wt9GhwC8fFQuxweowqFQnQ9tGvZL6iUxWyp5w7MTI6C\nzqAmAB97LZg/rKp7VPWZcoqr9LRoASVsOiQtoK0BQLBjH5qSQigUIi3ATb4EtvwnppAZCsHo1P2a\nTWXgCCBmKzzxBDz+OGTOh8sugzvvhGbNSiPYvAW27FNwZR8t2734fxZaUIeFbwDtgYOA+SJyk4gM\nD7/KLUJjTGDFxMCoUbBqFdx0E7zxBrRsCVdd5TpQNKYghVWS2A38AVTHPfsU+TLGmCI55BB4+GHX\nF9VVV8GECZCQAMOHw4YNhU9vKqZ8E5SI9ATSgFpAe1W9U1XvCr/KLUJjzAHj8MNh9Gj44Qe45BJ4\n6ik46ijXL1Uo5Hd0JtoUdAZ1G3CBqo5U1R3lFZAx5sDXtCmMHw9Ll8KZZ8J997lKFPffD5mZfkdn\nokVB96C6eO3fGWNMmWjRAl5/Hb79Fk46CW67zVWgePJJ693XFP9BXWOMKXVJSfDee/DNN3DssTBs\nGDRvDmPHwp49fkdn/GIJyhgTNTp3hk8/hc8/hyZNXIWKVq1g8mTItr4LKhxLUMaYqHPyyfD11/D+\n+1Cnjnt+qm1bePttUGvnr8KwBGWMiUoirgLFwoXuPlV2NvTpA8cdBzNnWqKqCCxBGWOiWqVKcOGF\nsHixe35q0yY4/XTo1g2+/NLv6ExZitoEJSJrRGSxiKSJyHxv2MEi8omIrPT+1vc7TmNM+ahSBQYM\ncO0/jxkD6enQtSv07Anz5/sdnSkLUZugPCerapKqJnufRwKfqWpz4DPvszGmAqlWDYYMca1SPPqo\nS07HHQfnnQfff+93dKY0RXuCyq03MMl7Pwk4x8dYjDE+qlkTbrzRtfN3113w2WfQpg3ce28r0tP9\njs6UhuL2B1WeFPhURLKBF1R1LBCrquu88b8BsXlNKCKDgcEAsbGxpJawdeDMzMwST+u3IMceCoXI\nzs4ObPwQ3PIPatl37QpJSVV4/fUmvPVWI1q0UE4/fR2XX/4Thx66y+/wiiQUSgpk2Zcl0SitCiMi\njVQ1Q0QOBT4B/gm8q6r1Ir6zRVULvA+VnJys80t4gTo1qM32E+zYU8LdbaSl+R1KiQW1/FMmemU/\nNLhl//bb/+WLL07g+efd5yFD4JZbIDbPw9nokZKCt91HR0eXIrIg4vaKL6L2Ep+qZnh/NwDvAB2B\n9SISB+D9tXaQjTH7OPjg3Tz1FKxcCZdf7hqnPeoouPVW2LLF7+hMcURlghKR2iJSN/weOBX4HngX\n6O99rT8ww58IjTHR7ogjYNw4WLYMzjkHHnzQNUh7zz2wvfx6oDf7ISoTFO7e0lci8h0wD/hAVWcC\nDwI9RGQl0N37bIwx+Wre3DWV9N13roWKf//bnVE99hj8+aff0ZmCRGWCUtVVqtrWe7VW1fu84ZtU\n9RRVba6q3VV1s9+xGmOC4dhj4Z13YO5caN/e9fCbkADPPw+7d/sdnclLVCYoY4wpKx07wqxZkJrq\nLvkNGeK6oX/5ZWuQNtpYgjLGVEjhppI+/BDq14f+/eGYY+CNN2DvXr+jM2AJyhhTgYm4dv3mz4c3\n33SfL7wQkpPhgw+sQVq/WYIyxlR4Iq6l9MWL3aW+rVvhrLNcL7/23Kx/LEEZY4yncmXo1w+WL3eV\nJ376ydX869HDVa4w5csSlDHG5FK1quvNd+VKePxxV0W9c2fo1QsWLfI7uorDEpQxxuSjZk0YNsw1\nSHvvvfCf/7ieffv2dd1+mLJlCcoYYwpRpw7cdhusXu2aTHr/fUhMhIED3WVAUzYsQRljTBHVrw/3\n3efOqK6/Hl57zbVUcd11sG5d4dOb4rEEZYwxxXToofDEE+4e1f/9H7zwAjRrBiNGuC7pTemwBGWM\nMSXUpIlLTsuXu2rqjz7qWqcYNQq2bfM7uuCzBGWMMfupWTN45RX3HFWPHq6H3/h4ePhh2LHD7+iC\nyxKUMcaUktat4a23XMsUHTvCzTe75DV6NOwKRse+UcUSlDHGlLIOHeCjj1y19KOPhn/+E1q0gPHj\nISvL7+iCwxKUMcaUkS5dXFNJs2ZBw4YwaJA7y5o61RqkLYrAJSgR6SkiK0QkXURG+h2PMcYURARO\nPRXmzXP9UVWrBhdfDO3awbvvWoO0BQlUghKRysCzwOlAInCxiCT6G5UxxhROxHU9n5bmevjdsQN6\n94bjj4fPPvM7uuhUxe8AiqkjkK6qqwBEZCrQG1haFgsbOnQo9erVK4tZl7lQKBTY2NPS0sjKyiIl\nJcXvUEosqOWflmRlX17i4ipTqVJP0tIup3v3WCpX3k7VqutISTmnVOafegA0wx60BNUI+CXi81qg\nU+4vichgYDBAbGxsif9R2dnZhEKhEk3rtyDHnpWVhaoGNn4IbvlnZWWhWNmXl5o1J3P00W+waVMf\nfv11ONnZrVi48B4OO2w0tWrtX2N/lqCilKqOBcYCJCcna0mPBp955pnAHkmmpqYGNvaUlBRCoRBp\naWl+h1JiQS3/lIle2Y+2si9vXbrAjz/+yc6dXVi5sgvnnw933w2tWvkdmX8CdQ8KyACaRHxu7A0z\nxphAq1wZDj10F6tWwR13wMyZrgv6AQNcI7UVUdAS1P+A5iISLyLVgL7Auz7HZIwxpaZePXfmtGqV\n6+pj6lT3DNU118Cvv/odXfkKVIJS1SzgOmAWsAyYpqpL/I3KGGNKX8OGrm2/H390z0+NG+dapbjp\nJti40e/oykegEhSAqn6oqkerajNVvc/veIwxpiw1agTPPec6SLzwQteK+j33+B1V+TggK0kYY8yB\n5qijYNIkGDnS9UtVEVTIBLVnzx7Wrl3Lzp07C/xeTEwMy5YtK6eoSleQY7/zzjtR1QLjr1GjBo0b\nN6Zq1arlGJkx/qtItfoqZIJau3YtdevWpWnTpohIvt/bvn07devWLcfISk+QY69UqRJZWVm0yueX\nqKps2rSJtWvXEh8fX87RGWPKS+DuQZWGnTt30qBBgwKTk4leIkKDBg0KPQM2xgRbhUxQgCWngLP/\nnzEHvgqboIwxxkQ3S1A+mj59OiLC8uXLizxNSkoKLVq0ICkpiVatWjF27NhSi+eTTz6hQ4cOHHvs\nsXTo0IHPP/88Z9yCBQs49thjSUhI4Prrr0e9PgJ27drFRRddREJCAp06dWLNmjWlFo8xpmKzBOWj\nKVOmcNJJJzFlypS/jZs4cSKjRo3Kc7rJkyeTlpbG119/zc0338zu3bsB2LJly37Fc8ghh/Dee++x\nePFiJk2aRL9+/XLGDRkyhHHjxrFy5UpWrlzJzJkzAXjppZeoX78+6enpDBs2jJtvvnm/YjDGmLAK\nWYsv0tChrn+WvGRn16Ry5eLPMykJnnyy4O9kZmby1VdfMXv2bM4++2zuuuuuYi8nMzOT2rVrU9kL\n8pxzziEmJoYrrriCLl26FHt+7dq1y3nfunVr/vzzT3bt2sXmzZvZtm0bnTt3BuDyyy9n+vTpnH76\n6cyYMSMnkZ5//vlcd911qKrdIzLG7Dc7g/LJjBkz6NmzJ0cffTQNGjRgwYIFRZ720ksvpU2bNrRo\n0YI77rgjJ0GlpqYyfPhw3nzzTZKTk7n11ltJT0/Pme6RRx4hKSnpb6/rr7/+b8t46623aN++PdWr\nVycjI4PGjRvnjGvcuDEZGa6N3oyMDJo0ce33VqlShZiYGDZt2lSiMjHGmEgV/gyqoDOd7dv/LLNn\niaZMmcINN9wAQN++fZkyZQpNmzbllFNOAWDz5s3s3r2b6dOnA/DKK69w7LHHAu4SX3JyMr///jsn\nnHACPXv25Mgjj0RESElJISUlhYyMDMaMGUPLli15/fXX6dOnD//617/417/+VWhsS5Ys4eabb+bj\njz8uk3U3xpiiqPAJyg+bN2/m888/Z/HixYgI2dnZiAiPPPJITh9IEydOZM2aNfnehwJo2LAh7du3\nZ+7cuRx55JEA/Pnnn7zzzjuMGzeO7du389RTT9GjRw/AnUFNnjz5b/Pp2rUrTz/9NOAeYj733HN5\n+eWXadasGQCNGjVi7dq1Od9fu3YtjRo1yhn3yy+/0LhxY7Kysti6dSsNGjTY/0IyxlR4donPB2++\n+Sb9+vXjp59+Ys2aNfzyyy/Ex8fz5ZdfFms+O3bs4Ntvv81JJCNGjCAxMZH//ve/3HPPPcyfP59r\nr72Wgw46CIB//etfpKWl/e0VTk6hUIgzzzyTBx98kBNPPDFnOXFxcRx00EHMmTMHVeXll1+md+/e\nAPTq1YtJkyblrNc//vEPu/9kjCkVlqB8MGXKFM4999x9hvXp0yfP2nx5ufTSS0lKSqJDhw4MGDCA\nDh06AK4K+rJlyxg9ejRt27YtdlyjR48mPT2du+++O+f+1IYNGwAYM2YMV1xxBQkJCTRr1ozTTz8d\ngEGDBrFp0yYSEhJ4/PHHefDBB4u9XGOMyYtd4vPB7Nmz/zYsd0WFAQMG5DltampqvvM944wz9ics\nbr/9dm6//fY8xyUnJ/P999//bXiNGjV444039mu5xhiTFzuDMsYYE5WiLkGJyCgRyRCRNO91RsS4\nW0QkXURWiMhpfsZpjDGmbEXrJb4nVPXRyAEikgj0BVoDhwOfisjRqprtR4DGGGPKVtSdQRWgNzBV\nVXep6mogHejoc0zGGGPKSLSeQf1TRC4H5gM3quoWoBEwJ+I7a71hfyMig4HBALGxsX+rWBATE8P2\n7dsLDSI7O7tI34tGQY49KysLoND4d+7cWWClET9lZmZGbWwFCYVCZGdnBzL2sOCWfVLgy760+ZKg\nRORT4LA8Rt0GPAfcA6j39zFgYHHmr6pjgbEAycnJmpKSss/4ZcuWFamFiCD3Shvk2KtUqUJWVlah\n8deoUWOf9gOjSWpqKrm3uyCot6YeoVAokLGHBbbs6xH4si9tvlziU9XuqnpMHq8ZqrpeVbNVdS8w\njr8u42UATSJm09gbFji//PILJ598MomJibRu3ZqnnnoKcC1M9OjRg+bNm9OjR499Wid/4IEHSEhI\noEWLFsyaNavUYxo6dCiNGjVi7969pTrf0aNHk5CQgIiwcePGnOGqyvXXX09CQgJt2rRh4cKFOeO+\n/PJLzjrrLBISEuy5KmMqsKi7ByUicREfzwXCD9+8C/QVkeoiEg80B+aVd3yloUqVKjz22GMsXbqU\nOXPm8Oyzz7J06VIefPBBTjnlFFauXMkpp5ySs3NeunQpU6dOZcmSJcycOZNrrrmG7GxXN2R/u9gA\n2Lt3L++88w5NmjThiy++yBm+devW/U5YJ554Ip9++mlOU0xhH330UU7XHWPHjmXIkCGAuzR59913\n89xzz7F06VKmTJnC0qVL9ysGY0wwReM9qIdFJAl3iW8NcBWAqi4RkWnAUiALuLZUavAV0N9Gzexs\nyqK/jbi4OOLiXB6uW7curVq1IiMjgxkzZuRcf+7fvz8pKSk89NBDzJgxg759+1K9enXi4+NJSEhg\n3rx5HH/88ft0sXHGGWdQpUrx/6Wpqam0bt2aiy66iClTpnDyyScD8NVXXzF06FAuvfRSBg4cyBFH\nHFHseed3CW7GjBlcfvnliAidO3cmFAqxbt061qxZwxFHHEGTJk2oVq0affv2ZcaMGSQmJhZ72caY\nYIu6MyhV7aeqx6pqG1XtparrIsbdp6rNVLWFqn7kZ5ylZc2aNXz77bd06tSJ9evX5ySuww47jPXr\n1wP7dmkB+3Z3EdnFRqtWrUrUxcaUKVO4+OKLOffcc/nggw/Ys2cPAGeeeSbffPMNMTEx9OrVi549\ne/LGG2/kdJC4ffv2POedlJRU6FlPfuuUkZGRUwa519UYU7FE4xlU+SrgTOfPMq5okJmZSZ8+fXjy\nySdzGnQNE5EiNboa2cXGtm3beOihh2jZsiUTJ07ksssuK7SLjd27d/Phhx/y+OOPU7duXTp16sSs\nWbM466yzANfL7rBhwxg2bBjffPMNAwcO5J577mHRokXUrVs3p/V1Y4wpbZagfLJnzx769OnDpZde\nynnnnQe4KvHr1q0jLi6OdevWceihhwJ/dWkRFtndBfzVxcb48eMJhUI89dRTOZfpCutiY9asWYRC\noZy+pnbs2EHNmjVzEhS4e2ATJkxg+vTpdOvWjSuvvBJwZ1D59dz72muvFXhZLr912rNnD+vWrfvb\ncGNMxWMJygeqyqBBg2jVqhXDhw/PGR7uumLkyJFMmjRpny4tLrnkEoYPH86vv/7KypUr6djRVW4c\nMWIEb7zxBmeeeSaPPPJIzj2f8DNEhZ1BTZkyhRdffJGLL74YgD/++IP4+Hh27NjB8uXLueaaa6hU\nqRKDBg3i22+/pU6dOjnT7s8ZVK9evRg9ejR9+/Zl7ty5xMTEEBcXR8OGDfnpp59Yu3YtzZs3Z+rU\nqbz22mslWoYxJtgsQfng66+/zukhNykpCYD777+fkSNHcuGFF/LSSy9x5JFHMm3aNABat27NhRde\nSGJiIlWqVOHZZ5/N6eY9JSWFu+++mxo1ahQ7jh07djBz5kyef/75nGG1a9fmpJNO4r333qNNmzZM\nmDCBVq1alXhdn376aR5++GF+++032rRpwxlnnMGLL77IGWecwYcffkhCQgK1atViwoQJgKvheMcd\nd3DVVVdRuXJlBg4cSOvWrUu8fGNMcImq+h1DmUpOTtb58+fvM2zZsmVF2ukG+WHXIMe+YsUKsrKy\nCk1MRf0/+iGoD4umTEwhFAqRNjS49xYDW/Yp7kHdtLR6focCgIgsUNVkP2OIulp8xhhjDFiCMsYY\nE6UsQRljjIlKlqCMMcZEJUtQxhhjopIlKGOMMVHJEpRPBg4cyKGHHsoxxxyTM+yiiy7KacuuadOm\nOc9IQdl1t5GamkpMTAxJSUm0adOG7t27s2HDhlKb/2effUb79u1JSkripJNOymknsKDuNowxBixB\n+WbAgAHMnDlzn2Gvv/46aWlppKWl0adPn5wmkMq6u40uXbqQlpbGokWLOO6443j22WeB0uluY8iQ\nIUyePJm0tDQuueQS7r33XiD/7jaMMSaswrcksXLoSjLTMvMcl52dndNiQ3HUSapD8yebF/idrl27\nsmbNmjzHqSrTpk3j888/Byjz7jYil7t9+3YSEhKA0uluQ0TYtm0b4BLe4YcfnrNOeXW3EdmSuTGm\nYqvwCSoaffnll8TGxtK8uUtyGRkZdO7cOWd87u42vvjiC8aPH8+NN97IBRdcwMCBA4mNjQUKbyw2\nvLykpCQ2bdpE7dq1uf/++wHX3UanTp145ZVX6NWrF4cddhiDBg2id+/eVKtWrUiNxYabNapZsyYH\nHXQQc+bMyVmnvLrbsARljAmr8AmqoDMdv5oLCvfPVBT7290GuEt877//PgAPPfQQI0aMyGmfb3+7\n23jiiSf48MMP6dSpE4888gjDhw/nxRdfLNK6GWMqNl8SlIhcAIwCWgEdVXV+xLhbgEFANnC9qs7y\nhncAJgI1gQ+BG/QAbEgwKyuLt99+mwULFuQMK8vuNnLr1asXffr02WdYSbvbaNiwId999x2dOnUC\nXCWQnj17FmmdjDHGrzOo74HzgBciB4pIItAXaA0cDnwqIkd7Xbs/B1wJzMUlqJ7AAdGrbqRPP/2U\nli1b0rhx45xhZdndRm5fffUVzZo1A2DhwoX71d1GVlYWW7du5YcffuDoo4/mk08+yWncNb/uNowx\nJsyXBKWqy4C8eoztDUxV1V3AahFJBzqKyBrgIFWd4033MnAOAU5QF198MampqWzcuJHGjRtz1113\nMWjQIKZOnfq3y3tl1d1GWPgelKoSExOTcwmuZs2a+9XdRpUqVRg3bhx9+vShUqVK1K9fn/HjxwPk\n292GMcaE+drdhoikAjeFL/GJyGhgjqq+6n1+CZeE1gAPqmp3b3gX4GZVPSuf+Q4GBgPExsZ2mDp1\n6j7jY2JicmqqFaSktfiiQZBjh6LFn56eztatW8spouLJzMzc52wzKIamDSU7O5tnOjzjdyglFtiy\nH5rkyv6ZxX6HAsDJJ5/se3cbZXYGJSKfAoflMeo2VZ1RVssFUNWxwFhw/UHl7htm2bJlRar8EOQ+\nlYIcOxQt/ho1auRc0ow2Qe2TqN6aeoRCoUDGHhbYsq9H4Mu+tJVZggqf7RRTBtAk4nNjb1iG9z73\ncGOMMQeoaGtJ4l2gr4hUF5F4oDkwT1XXAdtEpLO4G1eXA/t1FnYAVgCsUOz/Z8yBz5cEJSLnisha\n4HjgAxGZBaCqS4BpwFJgJnCtV4MP4BrgRSAd+JH9qCBRo0YNNm3aZDu5gFJVNm3atF8VQ4wx0c+v\nWnzvAO/kM+4+4L48hs8Hjvn7FMXXuHFj1q5dy++//17g93bu3BnYnWCQY4fC469Ro8Y+VfGNMQee\nCtmSRNWqVYmPjy/0e6mpqVF7E74wQY4dgh+/MWb/Rds9KGOMMQawBGWMMSZKWYIyxhgTlXxtSaI8\niMjvwE8lnPwQYGMphlOeghw7WPx+CnLsEOz4oyn2I1W1oZ8BHPAJan+IyHy/m/ooqSDHDha/n4Ic\nOwQ7/iDHXhbsEp8xxpioZAnKGGNMVLIEVbCxfgewH4IcO1j8fgpy7BDs+IMce6mze1DGGGOikp1B\nGWOMiUqWoIwxxkQlS1CAiFwgIktEZK+IJOcad4uIpIvIChE5LWJ4BxFZ7I17WvLov94PIjJKRDJE\nJM17nRExLs91iSYi0tOLL11ERvodT1GIyBpvW0gTkXDv0AeLyCcistL7W9/vOMNEZLyIbBCR7yOG\n5RtvNG03+cQeiG1eRJqIyGwRWertb27whgei7H2hqhX+BbQCWgCpQHLE8ETgO6A6EI/r5qOyN24e\n0BkQXNcfp/u9Hl5co4Cb8hie77pEywuo7MV1FFDNizfR77iKEPca4JBcwx4GRnrvRwIP+R1nRGxd\ngfbA94XFG23bTT6xB2KbB+KA9t77usAPXoyBKHs/XnYGBajqMlVdkceo3sBUVd2lqqtxfVF1FJE4\n4LPA/jMAAARoSURBVCBVnaNuS3oZOKccQy6JPNfF55hy6wikq+oqVd0NTMXFHUS9gUne+0lE0fah\nqv8BNucanF+8UbXd5BN7fqIt9nWqutB7vx1YBjQiIGXvB0tQBWsE/BLxea03rJH3PvfwaPFPEVnk\nXQ4JXy7Ib12iSRBizIsCn4rIAhEZ7A2LVdcTNMBvQKw/oRVZfvEG5X8SqG1eRJoC7YC5BL/sy0yF\nSVAi8qmIfJ/HK3BH6IWsy3O4S2RJwDrgMV+DrRhOUtUk4HTgWhHpGjnSO8sOzPMcQYuXgG3zIlIH\neAsYqqrbIscFsOzLVIXpsFBVu5dgsgygScTnxt6wDO997uHloqjrIiLjgPe9j/mtSzQJQox/o6oZ\n3t8NIvIO7jLMehGJU9V13iXhDb4GWbj84o36/4mqrg+/j/ZtXkSq4pLTZFV92xsc2LIvaxXmDKqE\n3gX6ikh1EYkHmgPzvNPxbSLS2au9dzkww89Aw7wNPOxcIFzbKc91Ke/4CvE/oLmIxItINaAvLu6o\nJSK1RaRu+D1wKq7M3wX6e1/rT5RsHwXIL96o326Css17+4qXgGWq+njEqMCWfZnzu5ZGNLxwG/Va\nYBewHpgVMe42XO2ZFUTU1AOScT+EH4H/b+9+XmyM4jiOvz9jY6EwFhqRH9mIqbEyC2WFZuQ/EAsG\nK5GZQVOykYUpGYOUHaL8AzSksNGk4U5WmmgWVpSV0sjX4pzJnWkGdee65zafVz11n3vP89zvuZ36\ndk7PPd9h8q4cjT6AO8A4UCEN8La/9aWkA+gmPd00AQw0Op5/iHcT6Umrt8C76ZiBVcBT4D3wBGht\ndKxVMd8nLYVN5XF/+E/xljRu5om9KcY8sJO0fFcB3uSju1l++0Yc3urIzMyK5CU+MzMrkhOUmZkV\nyQnKzMyK5ARlZmZFcoIyM7MiOUGZ1SDvUP1BUms+X5nPN9R43+OSDi5EjGbNyo+Zm9VIUj+wOSKO\nSroFfIyIS42Oy6zZeQZlVrsrQKekk6Q/Yw7ObiBpv6RXksbyXoqr8/tXJZ3Pr/dKei6pJdc46s3v\nn8g1hCqSHvzHfpk11KLZi8+sXiJiSlIf8AjYExFTczR7CXRGREg6AvQDp4FzwKikF8AQ0B0RPzWz\n/uVZYGNEfJe0oq6dMSuIZ1BmC6OLtAXPtnk+Xws8ljQO9AFbASLiG9ADjADDETExx7UV4J6kA8CP\nhQ7crFROUGY1ktQB7CZVWD4lqU3SxekS5LnZNVICageOAUurbtEOfAHWzPMV+4DrpEqyo5K88mGL\nghOUWQ3yDtU3SbV9JoHLwGBEDERER6Q6UQDL+V0q4VDV9etJS33bgS5JO2bdvwVYFxHPgDP5Psvq\n2SezUjhBmdWmB5iMiJF8fgPYImnXrHYXgIeSXgOfYUb5hd6I+ETamfu2pOrZ1RLgbl4aHAOGIuJr\n3XpjVhA/Zm5mZkXyDMrMzIrkBGVmZkVygjIzsyI5QZmZWZGcoMzMrEhOUGZmViQnKDMzK9IvIrzR\n5Al1KRkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7efd3a8c4b70>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " calculating for point: A = 200.000000 and B = 170.000000\n",
      "solution for z =  450000.0\n",
      "\n",
      " calculating for point: A = 100.000000 and B = 170.000000\n",
      "solution for z =  650000.0\n",
      "\n",
      " calculating for point: A = 100.000000 and B = 100.000000\n",
      "solution for z =  300000.0\n",
      "\n",
      " calculating for point: A = 120.000000 and B = 80.000000\n",
      "solution for z =  160000.0\n",
      "\n",
      " calculating for point: A = 200.000000 and B = 80.000000\n",
      "solution for z =  0.0\n",
      "\n",
      " the maximum profit z = $650000.000000 @ A = 100.000000 and B = 170.000000\n"
     ]
    }
   ],
   "source": [
    "import matplotlib.pyplot \n",
    "from matplotlib.pyplot import *\n",
    "import numpy\n",
    "from numpy import arange\n",
    "\n",
    "figure()\n",
    "A = arange(-100, 220, 10)\n",
    "B = arange(-100, 220, 10)\n",
    "\n",
    "# constraint equation\n",
    "\n",
    "B1 = 200.0 - A\n",
    "\n",
    "\n",
    "\n",
    "xlim = (-100, 220)\n",
    "ylim= (-100, 220)\n",
    "hlines(0, -100, 220, color = 'k')\n",
    "vlines(0, -100, 220, color = 'k')\n",
    "grid(True)\n",
    "\n",
    "xlabel('X-axis')\n",
    "ylabel('Y-axis')\n",
    "\n",
    "#Plotting graph\n",
    "\n",
    "plot(A, B1, color = 'b')\n",
    "axhline(y = 80, color = 'r')\n",
    "axhline(y = 170, color = 'm')\n",
    "axvline(x = 200, color = 'b')\n",
    "axvline(x = 100, color = 'g')\n",
    "\n",
    "title('objective function: z = 5000B - 2000A with following constarints')\n",
    "legend(['A+B>=200','200>=A>=100','170>=B>=80'])\n",
    "\n",
    "# get the co-ordinates of intersection points by mere visualisation\n",
    "A = [200.0,100.0, 100.0, 120.0,200.0]\n",
    "B = [170.0, 170.0, 100.0, 80.0, 80.0]\n",
    "fill(A,B,'g+')\n",
    "\n",
    "\n",
    "show()\n",
    "\n",
    "#Getting Solution\n",
    "checker = 0\n",
    "for i,j in zip(A,B):\n",
    "    print('\\n calculating for point: A = {0:f} and B = {1:f}' .format(i,j))\n",
    "    print('solution for z = ', 5000*j-2000*i)\n",
    "    if(checker <= (5000*j-2000*i)):\n",
    "       checker = (5000*j-2000*i)\n",
    "       X,Y = i,j\n",
    "       \n",
    "print('\\n the maximum profit z = ${0:f} @ A = {1:f} and B = {2:f}' .format(checker,X,Y))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Solving using library PULP:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "printing status of the LP problem:  Optimal\n",
      "Value of model A car:  100.0\n",
      "Value of model B car:  170.0\n",
      "the optimal solution or say maximum profit: $ 650000.0\n"
     ]
    }
   ],
   "source": [
    "from pulp import *\n",
    "# declare your variables\n",
    "A = LpVariable(\"A\", 100, 200)   # 100 <= A <= 200\n",
    "B = LpVariable(\"B\", 80, 170) # 80 <= B <= 170\n",
    " \n",
    "# defines the problem: optimization - Maximization\n",
    "prob = LpProblem(\"problem\", LpMaximize)\n",
    " \n",
    "# defines the constraints\n",
    "prob += A + B >=200 \n",
    "prob += A<=200\n",
    "prob += A>=100\n",
    "prob += B>=80\n",
    "prob += B<=170\n",
    " \n",
    "# defines the objective function to maximize\n",
    "prob += 5000*B-2000*A\n",
    " \n",
    "# solve the problem\n",
    "status = prob.solve()\n",
    "print('printing status of the LP problem: ', LpStatus[status])\n",
    " \n",
    "# print the results A = 100, B = 170\n",
    "print('Value of model A car: ', value(A))\n",
    "print('Value of model B car: ', value(B))\n",
    "print('the optimal solution or say maximum profit: $', value(prob.objective))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Solving using Scipy Library:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Positive directional derivative for linesearch    (Exit mode 8)\n",
      "            Current function value: -650000.0\n",
      "            Iterations: 7\n",
      "            Function evaluations: 12\n",
      "            Gradient evaluations: 3\n",
      "About solution from Optimize.minnimize: \n",
      "      fun: -650000.0\n",
      "     jac: array([ 2000., -5000.])\n",
      " message: 'Positive directional derivative for linesearch'\n",
      "    nfev: 12\n",
      "     nit: 7\n",
      "    njev: 3\n",
      "  status: 8\n",
      " success: False\n",
      "       x: array([ 100.,  170.])\n",
      "\n",
      "\n",
      "\n",
      "Printing the Results: \n",
      "\n",
      "\n",
      "\n",
      "maximum Profit Gained: $ 650000.0\n",
      "\n",
      " Number of Model A Cars:  100.0\n",
      "\n",
      " Number of Model B Cars:  170.0\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "from scipy.optimize import minimize\n",
    "\n",
    "#defining our Objective function\n",
    "def objective(x):\n",
    "    A = x[0]\n",
    "    B = x[1]\n",
    "    #since we are finding Max result using Minimizer. So returned the negative value.\n",
    "    return -(5000*B-2000*A)\n",
    "\n",
    "def constraint(x):\n",
    "    return x[0]+x[1]-200\n",
    "\n",
    "#boundations\n",
    "\n",
    "b1 = [100, 200]     #bounds on A: 200 >= A >= 100\n",
    "b2 = [80, 170]      #bounds on B: 170 >= B >= 80\n",
    "\n",
    "bnds = (b1,b2)\n",
    "\n",
    "#defining Inequality type of constraint.\n",
    "con = {'type': 'ineq', 'fun': constraint}\n",
    "\n",
    "# passing any guess value for A and B into minimize fuinction.\n",
    "Xguess = [100,100]\n",
    "\n",
    "#Predicting Maximization values for the objective function by seeding guess value Xguess; boundations included..\n",
    "sol = minimize(objective, Xguess, method = 'SLSQP', options={'disp':True}, bounds = bnds, constraints = con)\n",
    "\n",
    "print('About solution from Optimize.minnimize: \\n' ,sol)\n",
    "\n",
    "\n",
    "#assigning values to desired results.\n",
    "profit = -sol.fun   #converting back to positive answer\n",
    "CarA = sol.x[0]     # number of car A models\n",
    "CarB = sol.x[1]     # number of car B models\n",
    "\n",
    "###################################################################\n",
    "\n",
    "print ('\\n\\n\\nPrinting the Results: \\n\\n\\n')\n",
    "print('maximum Profit Gained: $', profit)\n",
    "print('\\n Number of Model A Cars: ',CarA)\n",
    "print('\\n Number of Model B Cars: ',CarB)\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}