car_rental/car_rental_2.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Car Rental 2\n",
    "\n",
    "## Objective and Prerequisites\n",
    "\n",
    "In this example, which is an extension of Vehicle Rental Optimization I, you’ll see how mathematical optimization can be used to figure out in which locations a car rental company should expand repair capacity. We’ll guide you through the process of creating a mixed integer linear programming model of this problem using the Gurobi Python API and then show you how you can find an optimal solution with the Gurobi Optimizer.\n",
    "\n",
    "This model is example 26 from the fifth edition of Model Building in Mathematical Programming by H. Paul Williams on pages 287 and 342-343.\n",
    "\n",
    "This example is at the intermediate level, where we assume that you know Python and the Gurobi Python API and that you have some knowledge of building mathematical optimization models.\n",
    "\n",
    "**Download the Repository** <br /> \n",
    "You can download the repository containing this and other examples by clicking [here](https://github.com/Gurobi/modeling-examples/archive/master.zip). "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Problem Description\n",
    "\n",
    "A small car rental company\n",
    "(which rents only one type of car) \n",
    "has depots in Glasgow, Manchester, Birmingham and Plymouth. There is an estimated demand for each day of the week except Sunday (when the company is closed). These estimates are given in the following table. It is not necessary to meet all demand.\n",
    "\n",
    "![weeklyDemand](weeklyDemand.PNG)\n",
    "\n",
    "\n",
    "Cars can be rented for one, two or three days and returned to either the\n",
    "depot from which they were rented or another depot at the start of the next morning. For\n",
    "example, a 2-day rental on Thursday means that the car has to be returned on\n",
    "Saturday morning; a 3-day rental on Friday means that the car has to be returned\n",
    "on Tuesday morning. A 1-day rental on Saturday means that the car has to be\n",
    "returned on Monday morning, and a 2-day rental on Tuesday morning.\n",
    "\n",
    "The rental period is independent of the origin and destination. From historical\n",
    "data, the company knows the distribution of rental periods: 55% of cars are rented\n",
    "for one day, 20% for two days and 25% for three days. The current estimates\n",
    "of percentages of cars rented from each depot and returned to a given depot\n",
    "(independent of the day) are given in the following table.\n",
    "\n",
    "![FromToPct](FromToPct.PNG)\n",
    "\n",
    "The company's marginal cost\n",
    "of renting out a car (‘wear and tear’, administration, etc.) is estimated as follows:\n",
    "\n",
    "| Days rented | Marginal cost |\n",
    "| --- | --- |\n",
    "| 1-day | $\\$ 20$ |\n",
    "| 2-day | $\\$ 25$ |\n",
    "| 3-day | $\\$ 30$ |\n",
    "\n",
    "The ‘opportunity cost’ (interest on capital, storage, servicing, etc.) of owning\n",
    "a car is $\\$ 15$ per week.\n",
    "\n",
    "It is possible to transfer undamaged cars from one depot to another depot,\n",
    "irrespective of distance. Cars cannot be rented out during the day in which they\n",
    "are transferred. The costs in USD, per car, of transfer are given in the following table.\n",
    "\n",
    "![FromToCst](FromToCst.PNG)\n",
    "\n",
    "Ten percent of cars returned by customers are damaged. When this happens,\n",
    "the customer is charged an excess of $\\$ 100$  (irrespective of the amount of damage\n",
    "that the company completely covers by its insurance). In addition, the car has\n",
    "to be transferred to a repair depot, where it will be repaired the following day.\n",
    "The cost of transferring a damaged car is the same as transferring an undamaged\n",
    "one (except when the repair depot is the current depot, \n",
    "in which case the cost would be $\\$0$). \n",
    "The transfer of a damaged car takes a day, unless it is already at a repair depot.\n",
    "Having arrived at a repair depot, all types of repair (or replacement) take a day.\n",
    "Only two of the depots have repair capacity. The (cars/day) capacity at each repair depot is as follows:\n",
    "\n",
    "| Repair depot | Capacity |\n",
    "| --- | --- |\n",
    "| Manchester | 12 |\n",
    "| Birmingham | 20 |\n",
    "\n",
    "Having been repaired, the car is available for rent at the depot the next day\n",
    "or may be transferred to another depot (taking a day). Thus, a car that is returned\n",
    "damaged on a Wednesday morning is transferred to a repair depot (if not the\n",
    "current depot) on Wednesday, repaired on Thursday and is available for rent\n",
    "at the repair depot on Friday morning.\n",
    "The rental price depends on the number of days for which the car is rented\n",
    "and whether it is returned to the same depot or not. The prices (in USD) are given in the following table.\n",
    "\n",
    "![RentalPrice](RentalPrice.PNG)\n",
    "\n",
    "We assume the following at the beginning of each day:\n",
    "1. Customers return cars that are due that day.\n",
    "2. Damaged cars are sent to the repair depot.\n",
    "3. Cars that were transferred from other depots arrive.\n",
    "4. Transfers are sent out.\n",
    "5. Cars are rented out.\n",
    "6. If it is a repair depot, then the repaired cars are available for rental.\n",
    "\n",
    "The car rental company wants to determine in which locations it should expand  repair capacity. The weekly fixed costs, given below, include interest payments on the necessary\n",
    "loans for expansion. The options are as follows:\n",
    "\n",
    "1. Expand repair capacity at Birmingham by 5 cars per day at a fixed cost\n",
    "per week of $\\$18,000$.\n",
    "2. Further expand repair capacity at Birmingham by 5 cars per day at a fixed\n",
    "cost per week of $\\$8,000$.\n",
    "3. Expand repair capacity at Manchester by 5 cars per day at a fixed cost per\n",
    "week of $\\$20,000$.\n",
    "4. Further expand repair capacity at Manchester by 5 cars per day at a fixed\n",
    "cost per week of $\\$5,000$.\n",
    "5. Create repair capacity at Plymouth of 5 cars per day at a fixed cost per\n",
    "week of $\\$19,000$.\n",
    "\n",
    "If any of these options is chosen, it must be carried out in its entirety; that\n",
    "is, there can be no partial expansion. Also, a further expansion at a depot can be\n",
    "carried out only if the first expansion is also carried out -so, for example, option\n",
    "(2) at Birmingham cannot be chosen unless option (1) is also chosen. If option (2)\n",
    "is chosen, thereby also choosing option (1), these count as two options. Similar\n",
    "stipulations apply regarding the expansions at Manchester. At most three of the\n",
    "options can be carried out.\n",
    "\n",
    "The goal is to determine in which locations to expand  repair capacity, the numbers of cars the car rental company should own, and where they should  be located at the start of each week day in order to maximize weekly profit. The company wants a ‘steady state’ solution in which the same expected number of cars will be located at the same depot on the same day of subsequent weeks."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Model Formulation\n",
    "\n",
    "$d,d2 \\in \\text{Depots}=\\{\\text{Glasgow}, \\text{Manchester}, \\text{Birmingham},  \\text{Plymouth}\\}$\n",
    "\n",
    "$\\text{NRD}=\\{\\text{Glasgow}\\}$: Depots without repair capacity.\n",
    "\n",
    "$\\text{RD}=\\{\\text{Manchester}, \\text{Birmingham}, \\text{Plymouth}\\}$: Depots with repair capacity.\n",
    "\n",
    "$t \\in \\text{Days}=\\{\\text{Monday},\\text{Tuesday},\\text{Wednesday},\\text{Thursday},\\text{Friday},\\text{Saturday}\\}$\n",
    "\n",
    "$r \\in \\text{RentDays}=\\{1,2,3\\}$: Number of days rented.\n",
    "\n",
    "### Parameters\n",
    "\n",
    "$\\text{demand}_{d,t} \\in \\mathbb{R}^+$: Estimated rental demand at depot $d$ on day $t$.\n",
    "\n",
    "$\\text{pctDepot}_{d,d2} \\in \\mathbb{R}^+$: Proportion of cars rented at depot $d$ to be returned to depot $d2$.\n",
    "\n",
    "$\\text{cstTransfer}_{d,d2} \\in \\mathbb{R}^+$: Transfer cost of a car from depot $d$ to depot $d2$.\n",
    "\n",
    "$\\text{pctRent}_{r} \\in \\mathbb{R}^+$: Proportion of cars rented for $r$ days.\n",
    "\n",
    "$\\text{capRepair}_{d} \\in \\mathbb{R}^+$: Repair capacity of depot $d$.\n",
    "\n",
    "$\\text{cstSameDepot}_{r} \\in \\mathbb{R}^+$: Rental cost for $r$ days with return to same depot.\n",
    "\n",
    "$\\text{cstOtherDepot}_{r} \\in \\mathbb{R}^+$: Rental cost for $r$ days with return to other depot.\n",
    "\n",
    "$\\text{marginalCost}_{r} \\in \\mathbb{R}^+$: Marginal cost to company of  $r$ days rental of a car.\n",
    "\n",
    "$\\text{pctUndamaged } \\in [0,1]$: Percent of cars returned by customers that are undamaged.\n",
    "\n",
    "$\\text{pctDamaged  } \\in [0,1]$: Percent of cars returned by customers that are damaged.\n",
    "\n",
    "$\\text{cstOwn} \\in \\mathbb{R}^+$: Cost of owning a car.\n",
    "\n",
    "$\\text{damagedFee} = 10$: Damaged car fee. Ten percent of the cars are damaged and the fee for a damaged car is $\\$100$.\n",
    "\n",
    "$\\text{cstExpCapB1} \\in \\mathbb{R}^+$: Fixed cost per week of expanding capacity at Birmingham.\n",
    "\n",
    "$\\text{cstExpCapB2} \\in \\mathbb{R}^+$: Fixed cost per week of further expanding capacity at Birmingham.\n",
    "\n",
    "$\\text{cstExpCapM1} \\in \\mathbb{R}^+$: Fixed cost per week of expanding capacity at Manchester.\n",
    "\n",
    "$\\text{cstExpCapM2} \\in \\mathbb{R}^+$: Fixed cost per week of further expanding capacity at Manchester.\n",
    "\n",
    "$\\text{cstExpCapP} \\in \\mathbb{R}^+$: Fixed cost per week of expanding capacity at Plymouth.\n",
    "\n",
    "### Decision Variables\n",
    "\n",
    "$\\text{xOwned} \\in \\mathbb{R}^+$: Total number of cars owned.\n",
    "\n",
    "$\\text{xUndamaged}_{d,t} \\in \\mathbb{R}^+$: Number of undamaged cars available at depot $d$ at the beginning of day $t$.\n",
    "\n",
    "$\\text{xDamaged}_{d,t} \\in \\mathbb{R}^+$: Number of damaged cars available at depot $d$ at the beginning of day $t$.\n",
    "\n",
    "$\\text{xRented}_{d,t} \\in \\mathbb{R}^+$: number of cars rented  out from depot $d$ at the beginning of day $t$.\n",
    "\n",
    "$\\text{xUDleft}_{d,t} \\in \\mathbb{R}^+$: Number of undamaged cars available at depot $d$ at the beginning of day $t$.\n",
    "\n",
    "$\\text{xDleft}_{d,t} \\in \\mathbb{R}^+$: Number of damaged cars left in depo $d$ at the end of day  $t$.\n",
    "\n",
    "$\\text{xUDtransfer}_{d,d2,t} \\in \\mathbb{R}^+$: Number of undamaged cars  at depot $d$ at the beginning of day $t$, to be transferred to depot $d2$. \n",
    "\n",
    "$\\text{xDtransfer}_{d,d2,t} \\in \\mathbb{R}^+$: Number of damaged cars at depot $d$ at the beginning of day $t$, to be transferred to depot $d2$. \n",
    "\n",
    "$\\text{xRepaired}_{d,t} \\in \\mathbb{R}^+$: Number of damaged cars to be repaired at depot $d$ during day $t$.\n",
    "\n",
    "$\\text{yExpandCapB1} \\in \\{0,1 \\}$: This binary variable is equal to 1, if the Birmingham repair capacity is expanded by 5 cars per day.\n",
    "\n",
    "$\\text{yExpandCapB2} \\in \\{0,1 \\}$: This binary variable is equal to 1, if the Birmingham repair capacity is expanded by a further 5 cars per day.\n",
    "\n",
    "$\\text{yExpandCapM1} \\in \\{0,1 \\}$: This binary variable is equal to 1, if the Manchester repair capacity is expanded by 5 cars per day.\n",
    "\n",
    "$\\text{yExpandCapM2} \\in \\{0,1 \\}$: This binary variable is equal to 1, if the Manchester repair capacity is expanded by a further 5 cars per day.\n",
    "\n",
    "$\\text{yExpandCapP} \\in \\{0,1 \\}$: This binary variable is equal to 1, if the Plymouth repair capacity is expanded by 5 cars per day."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Objective function\n",
    "The objective is to maximize profit.\n",
    "\n",
    "\\begin{equation}\n",
    "\\sum_{d \\in \\text{Depots}}\n",
    "\\sum_{t \\in \\text{Days}}\n",
    "\\sum_{r \\in \\text{RentDays}} \n",
    "\\text{pctDepot}_{d,d}*\\text{pctRent}_{r}*(\\text{cstSameDepot}_{r} - \\text{marginalCost}_{r} + \\text{damagedFee})*\\text{xRented}_{d,t}\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "+ \\sum_{d \\in \\text{Depots}} \\sum_{d2 \\in \\text{Depots}}\n",
    "\\sum_{t \\in \\text{Days}}\n",
    "\\sum_{r \\in \\text{RentDays}} \n",
    "\\text{pctDepot}_{d,d2}*\\text{pctRent}_{r}*(\\text{cstOtherDepot}_{r} - \\text{marginalCost}_{r} + \\text{damagedFee})*\\text{xRented}_{d,t}\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "- \\sum_{d \\in \\text{Depots}} \\sum_{d2 \\in \\text{Depots}}\n",
    "\\sum_{t \\in \\text{Days}} \\text{cstTransfer}_{d,d2}*(\\text{xUDtransfer}_{d,d2,t} + \\text{xDtransfer}_{d,d2,t} )\n",
    "- \\text{cstOwn}*\\text{xOwned}\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "-(\\text{cstExpCapB1}*\\text{yExpandCapB1} + \\text{cstExpCapB2}*\\text{yExpandCapB2} + \\text{cstExpCapM1}*\\text{yExpandCapM1} \n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "+ \\text{cstExpCapM2}*\\text{yExpandCapM2} + \\text{cstExpCapP}*\\text{yExpandCapP} )\n",
    "\\end{equation}\n",
    "\n",
    "### Constraints\n",
    "\n",
    "**Undamaged cars at a non-repair depot** <br />\n",
    "Number of undamaged cars available at a non-repair depot $d$ at the beginning of day $t$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\sum_{d2 \\in \\text{Depots}} \n",
    "\\sum_{r \\in \\text{RentDays}} \\text{pctUndamaged}*\\text{pctDepot}_{d2,d}*\\text{pctRent}_{r}*\\text{xRented}_{d2,(t-r)mod(6)}\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "+ \\sum_{d2 \\in \\text{Depots}} \\text{xUDtransfer}_{d2,d,(t-1)mod(6)} + \\text{xUDleft}_{d,(t-1)mod(6)} = \\text{xUndamaged}_{d,t} \n",
    "\\quad \\forall d \\in NRD, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "Demand of undamaged cars at the non-repair depot $d$ during day $t$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{xUndamaged}_{d,t} = \\text{xRented}_{d,t} + \n",
    "\\sum_{d2 \\in \\text{Depots}} \\text{xUDtransfer}_{d,d2,t} + \\text{xUDleft}_{d,t}\n",
    "\\quad \\forall d \\in NRD, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "**Undamaged cars at a repair depot** <br />\n",
    "Number of undamaged cars available at a repair depot $d$ at the beginning of day $t$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\sum_{d2 \\in \\text{Depots}} \n",
    "\\sum_{r \\in \\text{RentDays}} \\text{pctUndamaged}*\\text{pctDepot}_{d2,d}*\\text{pctRent}_{r}*\\text{xRented}_{d2,(t-r)mod(6)}\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "+ \\sum_{d2 \\in \\text{Depots}} \\text{xUDtransfer}_{d2,d,(t-1)mod(6)} \n",
    "+ \\text{xRepaired}_{d,(t-1)mod(6)}  + \\text{xUDleft}_{d,(t-1)mod(6)} = \\text{xUndamaged}_{d,t} \n",
    "\\quad \\forall d \\in NRD, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "Demand of undamaged cars at the repair depot $d$ during day $t$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{xUndamaged}_{d,t} = \\text{xRented}_{d,t} + \n",
    "\\sum_{d2 \\in \\text{Depots}} \\text{xUDtransfer}_{d,d2,t} + \\text{xUDleft}_{d,t}\n",
    "\\quad \\forall d \\in NRD, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "**Damaged cars at a non-repair depot** <br />\n",
    "Number of damaged cars available at a non-repair depot $d$ at the beginning of day $t$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\sum_{d2 \\in \\text{Depots}} \n",
    "\\sum_{r \\in \\text{RentDays}} \\text{pctDamaged}*\\text{pctDepot}_{d2,d}*\\text{pctRent}_{r}*\\text{xRented}_{d2,(t-r)mod(6)}\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "+ \\text{xDleft}_{d,(t-1)mod(6)} = \\text{xDamaged}_{d,t} \\quad \\forall d \\in NRD, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "Demand of undamaged cars at the non-repair depot $d$ during day $t$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{xDamaged}_{d,t} = \n",
    "\\sum_{d2 \\in \\text{Depots} \\cap RD} \\text{xDtransfer}_{d,d2,t} + \\text{xDleft}_{d,t}\n",
    "\\quad \\forall d \\in NRD, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "**Damaged cars at a repair depot** <br />\n",
    "Number of damaged cars available at a repair depot $d$ at the beginning of day $t$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\sum_{d2 \\in \\text{Depots}} \n",
    "\\sum_{r \\in \\text{RentDays}} \\text{pctDamaged}*\\text{pctDepot}_{d2,d}*\\text{pctRent}_{r}*\\text{xRented}_{d2,(t-r)mod(6)}\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "+ \\sum_{d2 \\in \\text{Depots}} \\text{xDtransfer}_{d2,d,(t-1)mod(6)} \n",
    "+ \\text{xDleft}_{d,(t-1)mod(6)} = \\text{xdamaged}_{d,t} \n",
    "\\quad \\forall d \\in RD, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "Demand of undamaged cars at the non-repair depot $d$ during day $t$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{xDamaged}_{d,t} = \\text{xRepaired}_{d,t} +\n",
    "\\sum_{d2 \\in \\text{Depots} \\cap RD} \\text{xDtransfer}_{d,d2,t} + \\text{xDleft}_{d,t}\n",
    "\\quad \\forall d \\in RD, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "**Repair Capacity** <br />\n",
    "Repair capacity of depot $d$ for each day $t$, including expansion capacity options. \n",
    "\n",
    "\\begin{equation}\n",
    "\\text{xRepaired}_{Birmingham,t} \\leq \\text{capRepair}_{Birmingham}\n",
    "+ 5*\\text{yExpandCapB1} + 5*\\text{yExpandCapB2}\n",
    "\\quad \\forall t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{xRepaired}_{Manchester,t} \\leq \\text{capRepair}_{Manchester}\n",
    "+ 5*\\text{yExpandCapM1} + 5*\\text{yExpandCapM2}\n",
    "\\quad \\forall t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{xRepaired}_{Plymouth,t} \\leq \\text{capRepair}_{Plymouth}\n",
    "+ 5*\\text{yExpandCapP}\n",
    "\\quad \\forall t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "**Depot Demand** <br />\n",
    "Demand at depot $d$ for each day $t$. \n",
    "\n",
    "\\begin{equation}\n",
    "\\text{xRented}_{d,t} \\leq \\text{demand}_{d,t}\n",
    "\\quad \\forall d \\in Depots, t \\in Days\n",
    "\\end{equation}\n",
    "\n",
    "**Number of cars** <br />\n",
    "Total number of cars owned equals number of cars rented out from all depots on Monday for 3 days, plus those on Tuesday for 2 or 3 days, plus all damaged and undamaged cars in depots at the beginning of Wednesday. \n",
    "Rationale: Let’s choose a day (Wednesday), count the cars undamaged and damaged that were returned to the depots and that are available on Wednesday morning. Let’s count the cars that have been rented and have not been returned: Cars rented on Monday for 3 days, and cars rented on Tuesday for 2 or 3 days.\n",
    "\n",
    "\\begin{equation}\n",
    "\\sum_{d \\in \\text{Depots}} (0.25*\\text{xRented}_{d,0} + 0.45*\\text{xRented}_{d,1} + \\text{xUndamaged}_{d,2}  + \\text{xdamaged}_{d,2} ) = \\text{xOwned}\n",
    "\\end{equation}\n",
    "\n",
    "**Capacity expansion constraints** <br />\n",
    "Additional expansion at a depot can be carried out only if the first expansion is also carried out.\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{yExpandCapB1} \\geq \\text{yExpandCapB2}\n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{yExpandCapM1} \\geq \\text{yExpandCapM2}\n",
    "\\end{equation}\n",
    "\n",
    "**Maximum number of expansions** <br />\n",
    "At most three of the options can be carried out.\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{yExpandCapB1} + \\text{yExpandCapB2} +\n",
    "\\text{yExpandCapM1} + \n",
    "\\end{equation}\n",
    "\n",
    "\\begin{equation}\n",
    "\\text{yExpandCapM2} +\n",
    "\\text{yExpandCapP} \\leq 3\n",
    "\\end{equation}\n",
    "\n",
    "\n",
    "---\n",
    "## Python Implementation\n",
    "\n",
    "We import the Gurobi Python Module and other Python libraries."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "%pip install gurobipy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "from itertools import product\n",
    "\n",
    "import gurobipy as gp\n",
    "from gurobipy import GRB\n",
    "\n",
    "# tested with Python 3.7.0 & Gurobi 9.0"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input data\n",
    "\n",
    "We define all the input data for the model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# list of depots and working days of a week\n",
    "\n",
    "depots = ['Glasgow','Manchester','Birmingham','Plymouth']\n",
    "NRD = ['Glasgow'] # Non-repair depot\n",
    "RD =['Manchester','Birmingham','Plymouth'] # Repair depot\n",
    "\n",
    "days = [0,1,2,3,4,5] # Monday = 0, Tuesday = 1, ...  Saturday = 5\n",
    "rentDays = [1,2,3]\n",
    "\n",
    "d2w, demand = gp.multidict({\n",
    "    ('Glasgow',0): 100,\n",
    "    ('Glasgow',1): 150,\n",
    "    ('Glasgow',2): 135,\n",
    "    ('Glasgow',3): 83,\n",
    "    ('Glasgow',4): 120,\n",
    "    ('Glasgow',5): 230,\n",
    "    ('Manchester',0): 250,\n",
    "    ('Manchester',1): 143,\n",
    "    ('Manchester',2): 80,\n",
    "    ('Manchester',3): 225,\n",
    "    ('Manchester',4): 210,\n",
    "    ('Manchester',5): 98,\n",
    "    ('Birmingham',0): 95,\n",
    "    ('Birmingham',1): 195,\n",
    "    ('Birmingham',2): 242,\n",
    "    ('Birmingham',3): 111,\n",
    "    ('Birmingham',4): 70,\n",
    "    ('Birmingham',5): 124,\n",
    "    ('Plymouth',0): 160,\n",
    "    ('Plymouth',1): 99,\n",
    "    ('Plymouth',2): 55,\n",
    "    ('Plymouth',3): 96,\n",
    "    ('Plymouth',4): 115,\n",
    "    ('Plymouth',5): 80\n",
    "})\n",
    "\n",
    "#repairCap\n",
    "depots, capacity = gp.multidict({\n",
    "    ('Glasgow'): 0,\n",
    "    ('Manchester'): 12,\n",
    "    ('Birmingham'): 20,\n",
    "    ('Plymouth'): 0\n",
    "})\n",
    "\n",
    "# Create a dictionary to capture \n",
    "# pctRent: percentage of cars rented for r days \n",
    "# cstMarginal: marginal cost for renting a car for r days\n",
    "# prcSameD: price of renting a car r days and returning to same depot\n",
    "# prcOtherD: price of renting a car r days and returning to another depot\n",
    "rentDays, pctRent, costMarginal, priceSameD, priceOtherD = gp.multidict({\n",
    "    (1): [0.55,20,50,70],\n",
    "    (2): [0.20,25,70,100],\n",
    "    (3): [0.25,30,120,150]\n",
    "})\n",
    "\n",
    "# Cost of owing a car per week.\n",
    "cstOwn = 15\n",
    "\n",
    "# damaged car fee\n",
    "damagedFee = 10\n",
    "\n",
    "# Create a dictionary to capture the proportion of cars rented at depot d to be returned to depot d2 \n",
    "d2d, pctFromToD = gp.multidict({\n",
    "    ('Glasgow','Glasgow'): 0.6,\n",
    "    ('Glasgow','Manchester'): 0.2,\n",
    "    ('Glasgow','Birmingham'): 0.1,\n",
    "    ('Glasgow','Plymouth'): 0.1,\n",
    "    ('Manchester','Glasgow'): 0.15,\n",
    "    ('Manchester','Manchester'): 0.55,\n",
    "    ('Manchester','Birmingham'): 0.25,\n",
    "    ('Manchester','Plymouth'): 0.05,\n",
    "    ('Birmingham','Glasgow'): 0.15,\n",
    "    ('Birmingham','Manchester'): 0.2,\n",
    "    ('Birmingham','Birmingham'): 0.54,\n",
    "    ('Birmingham','Plymouth'): 0.11,\n",
    "    ('Plymouth','Glasgow'): 0.08,\n",
    "    ('Plymouth','Manchester'): 0.12,\n",
    "    ('Plymouth','Birmingham'): 0.27,\n",
    "    ('Plymouth','Plymouth'): 0.53\n",
    "})\n",
    "\n",
    "# Create a dictionary to capture the transfer costs  of cars\n",
    "d2d, cstFromToD = gp.multidict({\n",
    "    ('Glasgow','Glasgow'): 0.001,\n",
    "    ('Glasgow','Manchester'): 20,\n",
    "    ('Glasgow','Birmingham'): 30,\n",
    "    ('Glasgow','Plymouth'): 50,\n",
    "    ('Manchester','Glasgow'): 20,\n",
    "    ('Manchester','Manchester'): 0.001,\n",
    "    ('Manchester','Birmingham'): 15,\n",
    "    ('Manchester','Plymouth'): 35,\n",
    "    ('Birmingham','Glasgow'): 30,\n",
    "    ('Birmingham','Manchester'): 15,\n",
    "    ('Birmingham','Birmingham'): 0.001,\n",
    "    ('Birmingham','Plymouth'): 25,\n",
    "    ('Plymouth','Glasgow'): 50,\n",
    "    ('Plymouth','Manchester'): 35,\n",
    "    ('Plymouth','Birmingham'): 25,\n",
    "    ('Plymouth','Plymouth'): 0.001\n",
    "})\n",
    "\n",
    "# Proportion of undamaged and damaged cars returned\n",
    "pctUndamaged = 0.9\n",
    "pctDamaged = 0.1\n",
    "\n",
    "# Expansion of capacity binary costs\n",
    "cstExpCapB1 = 18000\n",
    "cstExpCapB2 = 8000\n",
    "cstExpCapM1 = 20000\n",
    "cstExpCapM2 = 5000\n",
    "cstExpCapP = 19000\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Preprocessing\n",
    "We prepare the data structures to build the mixed integer linear programming model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Build a list of tuples (depot, depot2) such that d != d2\n",
    "list_d2notd = []\n",
    "\n",
    "for d,d2 in d2d:\n",
    "    if (d != d2):\n",
    "        tp = d,d2\n",
    "        list_d2notd.append(tp)\n",
    "\n",
    "d2notd = gp.tuplelist(list_d2notd)\n",
    "\n",
    "# Build a list of tuples (depot, depot2, day)\n",
    "list_dd2t = []\n",
    "\n",
    "for d,d2 in d2notd:\n",
    "    for t in days:\n",
    "        tp = d,d2,t \n",
    "        list_dd2t.append(tp)\n",
    "                    \n",
    "dd2t = gp.tuplelist(list_dd2t)\n",
    "\n",
    "# Build a list of tuples (depot, rent_day)\n",
    "list_dr = []\n",
    "\n",
    "for d in depots:\n",
    "    for r in rentDays:\n",
    "        tp = d,r\n",
    "        list_dr.append(tp)\n",
    "        \n",
    "dr = gp.tuplelist(list_dr) \n",
    "\n",
    "# Build a list of tuples (depot, day, rent_days )\n",
    "list_dtr = []\n",
    "\n",
    "for d in depots:\n",
    "    for t in days:\n",
    "            for r in rentDays:\n",
    "                tp = d,t,r\n",
    "                list_dtr.append(tp)\n",
    "                \n",
    "dtr = gp.tuplelist(list_dtr) \n",
    "\n",
    "# Build a list of tuples (depot, depot2, day, rent_days)\n",
    "list_dd2tr = []\n",
    "\n",
    "for d,d2 in d2notd:\n",
    "    for t in days:\n",
    "        for r in rentDays:\n",
    "            tp = d,d2,t,r\n",
    "            list_dd2tr.append(tp)\n",
    "                    \n",
    "                    \n",
    "dd2tr = gp.tuplelist(list_dd2tr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Model Deployment\n",
    "We create a model and the variables. The main decision variables are to choose the depots were repair capacity is \n",
    "worthwhile to expand, the number of cars to own\n",
    "and where should they be located at the start of each day of a week to maximize weekly profits."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Using license file c:\\gurobi\\gurobi.lic\n"
     ]
    }
   ],
   "source": [
    "model = gp.Model('RentalCar1')\n",
    "\n",
    "# Number of cars owned\n",
    "n = model.addVar(name=\"cars\")\n",
    "\n",
    "# Number of undamaged cars\n",
    "nu = model.addVars(d2w, name=\"UDcars\")\n",
    "\n",
    "# Number of damaged cars\n",
    "nd = model.addVars(d2w, name=\"Dcars\")\n",
    "\n",
    "# Number of cars hired (rented) cannot exceed their demand\n",
    "tr = model.addVars(d2w, ub=demand, name=\"Hcars\")\n",
    "#for d,t in d2w:\n",
    "    #tr[d,t].lb = 1\n",
    "\n",
    "# End inventory of undamaged cars\n",
    "eu = model.addVars(d2w, name=\"EUDcars\")\n",
    "\n",
    "# End inventory of damaged cars\n",
    "ed = model.addVars(d2w, name=\"EDcars\")\n",
    "\n",
    "# Number of undamaged cars transferred\n",
    "tu = model.addVars(dd2t, name=\"TUDcars\")\n",
    "\n",
    "# Number of damaged cars transferred\n",
    "td = model.addVars(dd2t, name=\"TDcars\")\n",
    "\n",
    "# Number of damaged cars repaired\n",
    "rp = model.addVars(d2w, name=\"RPcars\")\n",
    "\n",
    "# Number of damaged cars repaired cannot exceed depot capacity\n",
    "for d,t in d2w:\n",
    "    if d == 'Glasgow':\n",
    "        rp[d,t].ub = capacity[d] #repair capacity of Glasgow is zero\n",
    "        \n",
    "# Expansion of capacity binary variables\n",
    "zB1 = model.addVar(vtype=GRB.BINARY, name='expCapB1') # Expansion at Birmingham\n",
    "\n",
    "zB2 = model.addVar(vtype=GRB.BINARY, name='expCapB2') # Extra expansion at Birmingham\n",
    "\n",
    "zM1 = model.addVar(vtype=GRB.BINARY, name='expCapM1') # Expansion at Manchester\n",
    "\n",
    "zM2 = model.addVar(vtype=GRB.BINARY, name='expCapM2') # Extra expansion at Manchester\n",
    "\n",
    "zP = model.addVar(vtype=GRB.BINARY, name='expCapP') # Expansion at Plymouth"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Constraints\n",
    "The number of undamaged cars available at a non-repair depot d at the beginning of  day t should be equal to the demand of undamaged cars at the non-repair depot d during day t."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Undamaged cars into a non-repair depot constraints (left hand side of balance equation -availability)\n",
    "\n",
    "UDcarsNRD_L = model.addConstrs((gp.quicksum(pctUndamaged*pctFromToD[d2,d]*pctRent[r]*tr[d2,(t-r)%6 ] for d2,r in dr ) \n",
    "                              + gp.quicksum(tu.select('*',d,(t-1)%6)  ) \n",
    "                              + eu[d,(t-1)%6 ] == nu[d,t] for d in NRD for t in days ), \n",
    "                             name=\"UDcarsNRD_L\")\n",
    "\n",
    "# Undamaged cars out of a non-repair depot constraints (right hand side of balance equation -requirements)\n",
    "\n",
    "UDcarsNRD_R = model.addConstrs((tr[d,t] \n",
    "                                + gp.quicksum(tu.select(d,'*',t )) \n",
    "                                + eu[d,t] == nu[d,t] for d in NRD for t in days ), name='UDcarsNRD_R' )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The number of undamaged cars available at a repair depot d at the beginning of  day t should be equal to the demand of undamaged cars at the repair depot d during day t."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Undamaged cars into a repair depot constraints (left hand side of balance equation -availability)\n",
    "\n",
    "UDcarsRD_L = model.addConstrs((gp.quicksum(pctUndamaged*pctFromToD[d2,d]*pctRent[r]*tr[d2,(t-r)%6 ] for d2,r in dr ) \n",
    "                              + gp.quicksum(tu.select('*',d,(t-1)%6)  ) + rp[d, (t-1)%6 ]\n",
    "                              + eu[d,(t-1)%6 ] == nu[d,t] for d in RD for t in days ), \n",
    "                             name=\"UDcarsRD_L\")\n",
    "\n",
    "# Undamaged cars out of a repair depot constraints (right hand side of balance equation -requirements)\n",
    "\n",
    "UDcarsRD_R = model.addConstrs((tr[d,t] \n",
    "                                + gp.quicksum(tu.select(d,'*',t ) ) \n",
    "                                + eu[d,t] == nu[d,t] for d in RD for t in days ), name='UDcarsRD_R' )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The number of damaged cars available at a non-repair depot d at the beginning of  day t should be equal to the demand of damaged cars at the non-repair depot d during day t."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Damaged cars into a non-repair depot constraints (left hand side of balance equation -availability)\n",
    "\n",
    "DcarsNRD_L = model.addConstrs((gp.quicksum(pctDamaged*pctFromToD[d2,d]*pctRent[r]*tr[d2,(t-r)%6 ] for d2,r in dr ) \n",
    "                              + ed[d,(t-1)%6 ] == nd[d,t] for d in NRD for t in days ), \n",
    "                             name=\"DcarsNRD_L\")\n",
    "\n",
    "# Damaged cars out of a non-repair depot constraints (right hand side of balance equation -requirements)\n",
    "\n",
    "DcarsNRD_R = model.addConstrs(( gp.quicksum(td[d,d2,t] for d2 in RD ) \n",
    "                                + ed[d,t] == nd[d,t] for d in NRD for t in days ), name='DcarsNRD_R' )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The number of damaged cars available at a repair depot d at the beginning of  day t should be equal to the demand of damaged cars at the repair depot d during day t."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Damaged cars into a repair depot constraints (left hand side of balance equation -availability)\n",
    "\n",
    "DcarsRD_L = model.addConstrs((gp.quicksum(pctDamaged*pctFromToD[d2,d]*pctRent[r]*tr[d2,(t-r)%6 ] for d2,r in dr )\n",
    "                              + gp.quicksum(td[d2,d,(t-1)%6 ] for d2, dd in d2notd if (dd == d)) \n",
    "                              + ed[d,(t-1)%6 ] == nd[d,t] for d in RD for t in days ), \n",
    "                             name=\"DcarsRD_L\")\n",
    "\n",
    "# Damaged cars out of a repair depot constraints (right hand side of balance equation -requirements)\n",
    "\n",
    "DcarsND_R = model.addConstrs((rp[d,t] + gp.quicksum(td[d,d2,t ] for d2 in NRD ) \n",
    "                                + ed[d,t] == nd[d,t] for d in RD for t in days ), name='DcarsND_R' )"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Total number of cars equals the number of cars rented out from all depots on Monday for 3 days, plus those on Tuesday for 2 or 3 days, plus all damaged and undamaged cars in depots at the beginning of Wednesday."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Total number of cars owned constraint\n",
    "# Note: 25% of cars are rented for 3 days, and 20% + 25% = 45% of the cars are rented for 2-days or 3-days\n",
    "\n",
    "carsConstr = model.addConstr((gp.quicksum(0.25*tr[d,0] + 0.45*tr[d,1] + nu[d,2] + nd[d,2] for d in depots ) \n",
    "                              == n ),name='carsConstr')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Repair capacity expansion constraints."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Depot capacity constraints\n",
    "\n",
    "RepairCapB = model.addConstrs((rp['Birmingham',t] <= capacity['Birmingham'] \n",
    "                               + 5*(zB1 + zB2) for t in days ), name='RepairCapB')\n",
    "\n",
    "RepairCapM = model.addConstrs((rp['Manchester',t] <= capacity['Manchester'] \n",
    "                               + 5*(zM1 + zM2) for t in days ), name='RepairCapM')\n",
    "\n",
    "RepairCapP = model.addConstrs((rp['Plymouth',t] <= capacity['Plymouth'] \n",
    "                               + 5*zP for t in days ), name='RepairCapP')\n",
    "\n",
    "# Additional expansion constraints\n",
    "\n",
    "expandB1B2 = model.addConstr((zB1 >= zB2 ), name='expandB1B2')\n",
    "\n",
    "expandM1M2 = model.addConstr((zM1 >= zM2 ), name='expandM1M2')\n",
    "\n",
    "# At most three expansions are possible\n",
    "\n",
    "atMost3 = model.addConstr((zB1 + zB2 + zM1 + zM2 + zP <= 3), name='atMost3')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The objective function is to maximize profit."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Maximize profit objective function.\n",
    "\n",
    "model.setObjective((\n",
    "    gp.quicksum(pctFromToD[d,d]*pctRent[r]*(priceSameD[r] - costMarginal[r] + damagedFee)*tr[d,t] for d,t,r in dtr )\n",
    "    + gp.quicksum(pctFromToD[d,d2]*pctRent[r]*(priceOtherD[r]-costMarginal[r]+damagedFee)*tr[d,t] for d,d2,t,r in dd2tr)\n",
    "    - gp.quicksum(cstFromToD[d,d2]*tu[d,d2,t] for d,d2,t in dd2t) \n",
    "    - gp.quicksum(cstFromToD[d,d2]*td[d,d2,t] for d,d2,t in dd2t) \n",
    "    -(cstExpCapB1*zB1 + cstExpCapB2*zB2 + cstExpCapM1*zM1 + cstExpCapM2*zM2 + cstExpCapP*zP) - cstOwn*n ), GRB.MAXIMIZE)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Gurobi Optimizer version 9.1.0 build v9.1.0rc0 (win64)\n",
      "Thread count: 4 physical cores, 8 logical processors, using up to 8 threads\n",
      "Optimize a model with 118 rows, 294 columns and 1136 nonzeros\n",
      "Model fingerprint: 0xe441c90b\n",
      "Variable types: 289 continuous, 5 integer (5 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e-03, 5e+00]\n",
      "  Objective range  [2e+01, 2e+04]\n",
      "  Bounds range     [1e+00, 3e+02]\n",
      "  RHS range        [3e+00, 2e+01]\n",
      "Found heuristic solution: objective -0.0000000\n",
      "Presolve removed 49 rows and 91 columns\n",
      "Presolve time: 0.00s\n",
      "Presolved: 69 rows, 203 columns, 987 nonzeros\n",
      "Variable types: 198 continuous, 5 integer (5 binary)\n",
      "\n",
      "Root relaxation: objective 1.357349e+05, 58 iterations, 0.00 seconds\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "     0     0 135734.925    0    2   -0.00000 135734.925      -     -    0s\n",
      "H    0     0                    130734.92534 135734.925  3.82%     -    0s\n",
      "H    0     0                    131725.36880 135734.925  3.04%     -    0s\n",
      "H    0     0                    132341.46511 135734.925  2.56%     -    0s\n",
      "\n",
      "Explored 1 nodes (58 simplex iterations) in 0.03 seconds\n",
      "Thread count was 8 (of 8 available processors)\n",
      "\n",
      "Solution count 4: 132341 131725 130735 -0 \n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 1.323414651079e+05, best bound 1.323414651079e+05, gap 0.0000%\n"
     ]
    }
   ],
   "source": [
    "# Verify model formulation\n",
    "\n",
    "model.write('CarRental2.lp')\n",
    "\n",
    "# Run optimization engine\n",
    "\n",
    "model.optimize()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## Analysis"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The optimal number of cars to be owned is: 983\n",
      "The optimal profit is: $132,341.47\n"
     ]
    }
   ],
   "source": [
    "# Output report\n",
    "\n",
    "# Total number of cars owned\n",
    "print(f\"The optimal number of cars to be owned is: {round(n.x)}\")\n",
    "\n",
    "# Optimal profit\n",
    "print(f\"The optimal profit is: {'${:,.2f}'.format(round(model.objVal,2))}\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The depot at Manchester expanded its capacity to 5 more cars per day\n",
      "The depot at Manchester further expanded 5 more cars per day of additional capacity\n",
      "The depot at Plymouth expanded its capacity to 5 more cars per day\n"
     ]
    }
   ],
   "source": [
    "# Repair capacity expansion plan\n",
    "\n",
    "if zB1.x > 0.5:\n",
    "    print(f\"The depot at Birmingham expanded its capacity to 5 more cars per day\")\n",
    "    if zB2.x > 0.5:\n",
    "        print(f\"The depot at Birmingham further expanded 5 more cars per day of additional capacity\")\n",
    "\n",
    "if zM1.x > 0.5:\n",
    "    print(f\"The depot at Manchester expanded its capacity to 5 more cars per day\")\n",
    "    if zM2.x > 0.5:\n",
    "        print(f\"The depot at Manchester further expanded 5 more cars per day of additional capacity\")\n",
    "\n",
    "if zP.x > 0.5:\n",
    "    print(f\"The depot at Plymouth expanded its capacity to 5 more cars per day\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "_________________________________________________________________________________\n",
      "Estimated number of undamaged cars in depot at the beginning of each day: \n",
      "_________________________________________________________________________________\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Day</th>\n",
       "      <th>Glasgow</th>\n",
       "      <th>Manchester</th>\n",
       "      <th>Birmingham</th>\n",
       "      <th>Plymouth</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Monday</td>\n",
       "      <td>95</td>\n",
       "      <td>227</td>\n",
       "      <td>242</td>\n",
       "      <td>66</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Tuesday</td>\n",
       "      <td>99</td>\n",
       "      <td>162</td>\n",
       "      <td>282</td>\n",
       "      <td>65</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Wednesday</td>\n",
       "      <td>100</td>\n",
       "      <td>168</td>\n",
       "      <td>242</td>\n",
       "      <td>67</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Thursday</td>\n",
       "      <td>101</td>\n",
       "      <td>233</td>\n",
       "      <td>168</td>\n",
       "      <td>80</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Friday</td>\n",
       "      <td>117</td>\n",
       "      <td>226</td>\n",
       "      <td>162</td>\n",
       "      <td>71</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Saturday</td>\n",
       "      <td>103</td>\n",
       "      <td>182</td>\n",
       "      <td>232</td>\n",
       "      <td>67</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Day Glasgow Manchester Birmingham Plymouth\n",
       "     Monday      95        227        242       66\n",
       "    Tuesday      99        162        282       65\n",
       "  Wednesday     100        168        242       67\n",
       "   Thursday     101        233        168       80\n",
       "     Friday     117        226        162       71\n",
       "   Saturday     103        182        232       67"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Create a list to translate the number label of each day to the actual name of the day\n",
    "dayname = ['Monday','Tuesday','Wednesday','Thursday','Friday','Saturday']\n",
    "\n",
    "# Number of undamaged cars in depot at the beginning of each day.\n",
    "print(\"\\n\\n_________________________________________________________________________________\")\n",
    "print(f\"Estimated number of undamaged cars in depot at the beginning of each day: \")\n",
    "print(\"_________________________________________________________________________________\")\n",
    "\n",
    "undamaged_cars = pd.DataFrame(columns=['Day','Glasgow','Manchester','Birmingham','Plymouth'])\n",
    "for t in days:\n",
    "    undamaged_cars = undamaged_cars.append({\"Day\": dayname[t], \"Glasgow\": round(nu['Glasgow',t].x), \"Manchester\":  round(nu['Manchester',t].x), 'Birmingham':  round(nu['Birmingham',t].x),'Plymouth': round(nu['Plymouth',t].x)   }, ignore_index=True) \n",
    "undamaged_cars.index=[''] * len(undamaged_cars)\n",
    "undamaged_cars"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "_________________________________________________________________________________\n",
      "Estimated number of damaged cars in depot at the beginning of each day: \n",
      "_________________________________________________________________________________\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Day</th>\n",
       "      <th>Glasgow</th>\n",
       "      <th>Manchester</th>\n",
       "      <th>Birmingham</th>\n",
       "      <th>Plymouth</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Monday</td>\n",
       "      <td>12</td>\n",
       "      <td>22</td>\n",
       "      <td>20</td>\n",
       "      <td>7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Tuesday</td>\n",
       "      <td>11</td>\n",
       "      <td>24</td>\n",
       "      <td>20</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Wednesday</td>\n",
       "      <td>11</td>\n",
       "      <td>22</td>\n",
       "      <td>20</td>\n",
       "      <td>7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Thursday</td>\n",
       "      <td>11</td>\n",
       "      <td>25</td>\n",
       "      <td>24</td>\n",
       "      <td>9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Friday</td>\n",
       "      <td>19</td>\n",
       "      <td>22</td>\n",
       "      <td>20</td>\n",
       "      <td>11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Saturday</td>\n",
       "      <td>17</td>\n",
       "      <td>22</td>\n",
       "      <td>20</td>\n",
       "      <td>13</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Day Glasgow Manchester Birmingham Plymouth\n",
       "     Monday      12         22         20        7\n",
       "    Tuesday      11         24         20        8\n",
       "  Wednesday      11         22         20        7\n",
       "   Thursday      11         25         24        9\n",
       "     Friday      19         22         20       11\n",
       "   Saturday      17         22         20       13"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Number of Damaged cars in depot at the beginning of each day.\n",
    "print(\"_________________________________________________________________________________\")\n",
    "print(f\"Estimated number of damaged cars in depot at the beginning of each day: \")\n",
    "print(\"_________________________________________________________________________________\")\n",
    "\n",
    "damaged_cars = pd.DataFrame(columns=['Day','Glasgow','Manchester','Birmingham','Plymouth'])\n",
    "for t in days:\n",
    "    damaged_cars = damaged_cars.append({\"Day\": dayname[t], \"Glasgow\": round(nd['Glasgow',t].x), \"Manchester\":  round(nd['Manchester',t].x), 'Birmingham':  round(nd['Birmingham',t].x),'Plymouth': round(nd['Plymouth',t].x)   }, ignore_index=True) \n",
    "damaged_cars.index=[''] * len(damaged_cars)\n",
    "damaged_cars"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "_________________________________________________________________________________\n",
      "Estimated number of undamaged cars rented out from each depot and day: \n",
      "_________________________________________________________________________________\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Day</th>\n",
       "      <th>Glasgow</th>\n",
       "      <th>Manchester</th>\n",
       "      <th>Birmingham</th>\n",
       "      <th>Plymouth</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Monday</td>\n",
       "      <td>85</td>\n",
       "      <td>204</td>\n",
       "      <td>86</td>\n",
       "      <td>59</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Tuesday</td>\n",
       "      <td>89</td>\n",
       "      <td>129</td>\n",
       "      <td>176</td>\n",
       "      <td>58</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Wednesday</td>\n",
       "      <td>90</td>\n",
       "      <td>72</td>\n",
       "      <td>218</td>\n",
       "      <td>50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Thursday</td>\n",
       "      <td>75</td>\n",
       "      <td>203</td>\n",
       "      <td>100</td>\n",
       "      <td>72</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Friday</td>\n",
       "      <td>106</td>\n",
       "      <td>189</td>\n",
       "      <td>63</td>\n",
       "      <td>64</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <td>Saturday</td>\n",
       "      <td>93</td>\n",
       "      <td>88</td>\n",
       "      <td>112</td>\n",
       "      <td>61</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Day Glasgow Manchester Birmingham Plymouth\n",
       "     Monday      85        204         86       59\n",
       "    Tuesday      89        129        176       58\n",
       "  Wednesday      90         72        218       50\n",
       "   Thursday      75        203        100       72\n",
       "     Friday     106        189         63       64\n",
       "   Saturday      93         88        112       61"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Undamaged car rented out from each depot and day.\n",
    "print(\"_________________________________________________________________________________\")\n",
    "print(f\"Estimated number of undamaged cars rented out from each depot and day: \")\n",
    "print(\"_________________________________________________________________________________\")\n",
    "\n",
    "rentedOut = {}\n",
    "\n",
    "for d in depots:\n",
    "    for t in days:\n",
    "        count = 0\n",
    "        for d2 in depots:\n",
    "            for r in rentDays:\n",
    "                #print(f\"Depot {d}, day {t}: cars rented out {tr[d,t].x}\")\n",
    "                count += pctUndamaged*pctFromToD[d,d2]*pctRent[r]*tr[d,t].x\n",
    "        rentedOut[d,t] = round(count)\n",
    "    \n",
    "\n",
    "#print(rentedOut)\n",
    "\n",
    "rentout_cars = pd.DataFrame(columns=['Day','Glasgow','Manchester','Birmingham','Plymouth'])\n",
    "for t in days:\n",
    "    rentout_cars = rentout_cars.append({\"Day\": dayname[t], \"Glasgow\": round(rentedOut['Glasgow',t]), \"Manchester\":  round(rentedOut['Manchester',t]), 'Birmingham':  round(rentedOut['Birmingham',t]),'Plymouth': round(rentedOut['Plymouth',t])   }, ignore_index=True) \n",
    "rentout_cars.index=[''] * len(rentout_cars)\n",
    "rentout_cars"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "---\n",
    "## References\n",
    "\n",
    "H. Paul Williams, Model Building in Mathematical Programming, fifth edition.\n",
    "\n",
    "Copyright © 2020 Gurobi Optimization, LLC"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}