Single Level Reformulation.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "from docplex.mp.model import Model\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import random\n",
    "import copy\n",
    "import os\n",
    "import math\n",
    "from decimal import *\n",
    "# from docplex.mp.environment import Environment\n",
    "# Environment().print_information()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "def read_instance():\n",
    "    os.chdir(\"D:\\AA MASC\\C Codes\\FINAL\\INSTANCES\")\n",
    "    with open(\"test_10k_10ta_3.txt\", 'r' ) as file:         # Hardcoded file name \n",
    "        lines = file.readlines()\n",
    "        text_file = file.read()\n",
    "        file.close()\n",
    "        \n",
    "    i=1\n",
    "    TF = []\n",
    "    Rdash = []\n",
    "    RR = []\n",
    "\n",
    "    for line in lines:\n",
    "        if i==1:\n",
    "            J = int (line.split()[2])      # No. of Products \n",
    "            i+=1\n",
    "#             print(J)\n",
    "            continue\n",
    "\n",
    "        if i==2:\n",
    "            K = int (line.split()[2] )     # No. of Customers \n",
    "            i+=1\n",
    "#             print(K)\n",
    "            continue\n",
    "\n",
    "        if i==3:\n",
    "            i+=1\n",
    "            continue \n",
    "\n",
    "        if i == 4:                       # i = 4 , denotes 4th line \n",
    "    #         print(line)\n",
    "            for iter in range(0,K):\n",
    "#                 print(iter)\n",
    "#                 print(line.split()[iter+3 ])\n",
    "                TF.append( int ( line.split()[ iter+3 ] ) ) \n",
    "#                 print(TF)\n",
    "            i+=1\n",
    "            continue\n",
    "\n",
    "        if ( 4 < i < 5+ J):\n",
    "            for iter in range(0,K):\n",
    "#                 print('iter' , iter)\n",
    "                if i == 5:\n",
    "#                     print(line.split()[iter+3])\n",
    "                    Rdash.append( int ( line.split()[ iter+3  ] ) ) \n",
    "                else:\n",
    "#                     print(line.split()[iter])\n",
    "                    Rdash.append( int ( line.split()[ iter  ] ) ) \n",
    "            i+=1\n",
    "            RR.append(Rdash)\n",
    "            Rdash=[]        \n",
    "            continue\n",
    "        \n",
    "#         print(RR)\n",
    "        NestList = []\n",
    "        R_K =[] \n",
    "        BUDGETS = []\n",
    "        \n",
    "        for k in range(0,K):                   # For each customer \n",
    "            for j in range(0,J):   \n",
    "#                 print( TF[j])\n",
    "                NestList.append( TF[k] -  RR[j][k]) \n",
    "            R_K.append(NestList)\n",
    "            BUDGETS.append(max(NestList))\n",
    "            NestList=[]\n",
    "        \n",
    "        \n",
    "         \n",
    "        return K,J , R_K, BUDGETS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "K,J , R_K, BUDGETS = read_instance()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "# S_K = [[0, 20, 7, 9, 0, 12, 0, 10, 15, 19, 17, 13, 11, 16, 8, 14, 6, 18, 0, 0], [0, 7, 14, 15, 16, 18, 9, 10, 6, 17, 20, 19, 8, 0, 11, 0, 12, 0, 0, 13], [7, 11, 17, 6, 0, 0, 12, 14, 9, 10, 19, 18, 0, 13, 15, 16, 0, 8, 20, 0], [15, 19, 0, 12, 0, 14, 10, 16, 0, 0, 8, 20, 7, 13, 11, 0, 17, 6, 18, 9], [20, 7, 0, 17, 0, 0, 0, 0, 8, 9, 16, 10, 15, 11, 19, 6, 18, 14, 12, 13], [13, 0, 10, 7, 14, 18, 0, 0, 0, 16, 0, 11, 12, 8, 20, 15, 6, 17, 19, 9], [14, 0, 12, 10, 13, 17, 0, 7, 9, 0, 16, 18, 15, 8, 19, 0, 20, 0, 11, 6], [0, 11, 17, 19, 7, 0, 8, 20, 12, 0, 15, 18, 6, 0, 13, 10, 14, 0, 9, 16], [16, 0, 18, 11, 9, 0, 12, 17, 20, 19, 6, 15, 8, 10, 13, 14, 0, 0, 7, 0], [16, 18, 13, 0, 8, 11, 0, 17, 15, 6, 20, 19, 7, 10, 12, 0, 9, 0, 14, 0], [17, 0, 13, 8, 0, 15, 0, 20, 10, 16, 7, 9, 0, 0, 12, 19, 6, 14, 18, 11], [0, 0, 18, 11, 19, 9, 6, 12, 7, 0, 0, 20, 10, 0, 8, 16, 15, 13, 17, 14], [15, 0, 19, 8, 18, 20, 0, 0, 0, 12, 16, 17, 7, 14, 0, 10, 11, 6, 13, 9], [16, 0, 7, 15, 20, 6, 0, 11, 14, 8, 12, 0, 10, 0, 13, 0, 18, 17, 19, 9], [18, 15, 20, 16, 8, 0, 0, 7, 0, 0, 11, 17, 10, 14, 6, 9, 12, 19, 0, 13], [0, 17, 0, 15, 14, 9, 18, 6, 11, 20, 0, 12, 0, 0, 16, 13, 10, 19, 7, 8], [15, 11, 0, 17, 0, 20, 8, 18, 0, 7, 0, 6, 0, 10, 19, 12, 14, 16, 13, 9], [8, 10, 16, 17, 0, 19, 6, 11, 12, 14, 18, 7, 13, 0, 15, 0, 9, 20, 0, 0], [19, 0, 0, 13, 6, 0, 11, 0, 17, 7, 15, 20, 16, 0, 18, 9, 10, 12, 14, 8], [17, 14, 0, 10, 0, 9, 6, 0, 7, 16, 8, 0, 20, 15, 13, 18, 12, 11, 19, 0], [0, 10, 0, 17, 13, 11, 14, 20, 8, 19, 7, 0, 9, 18, 15, 0, 16, 0, 6, 12], [19, 17, 14, 16, 0, 0, 7, 11, 9, 6, 20, 0, 12, 15, 8, 18, 13, 10, 0, 0], [12, 6, 0, 7, 16, 17, 15, 18, 19, 11, 20, 0, 8, 10, 13, 0, 14, 0, 9, 0], [0, 17, 19, 11, 18, 13, 8, 0, 15, 0, 6, 7, 9, 10, 12, 14, 16, 0, 20, 0], [18, 11, 19, 20, 0, 15, 8, 0, 0, 16, 6, 10, 14, 0, 17, 7, 13, 0, 12, 9], [7, 16, 15, 6, 17, 8, 13, 12, 19, 14, 20, 18, 0, 9, 0, 0, 11, 10, 0, 0], [0, 6, 7, 0, 0, 10, 11, 12, 0, 14, 20, 17, 0, 8, 9, 16, 19, 18, 15, 13], [0, 14, 19, 13, 18, 6, 12, 17, 0, 0, 16, 11, 7, 15, 8, 10, 20, 0, 9, 0], [0, 7, 19, 0, 13, 8, 9, 11, 0, 18, 0, 10, 15, 17, 16, 20, 6, 12, 0, 14], [0, 13, 0, 9, 0, 19, 20, 7, 17, 14, 0, 12, 18, 0, 15, 6, 11, 10, 8, 16], [0, 0, 9, 7, 17, 14, 0, 16, 18, 0, 8, 15, 6, 20, 12, 19, 10, 11, 0, 13], [13, 10, 0, 15, 16, 8, 7, 6, 12, 14, 0, 0, 20, 19, 11, 18, 0, 0, 9, 17], [12, 17, 0, 19, 7, 0, 16, 20, 18, 14, 0, 0, 0, 8, 10, 13, 15, 11, 6, 9], [16, 8, 6, 15, 18, 13, 0, 20, 0, 12, 0, 9, 7, 17, 19, 11, 14, 0, 10, 0], [15, 17, 7, 9, 11, 13, 0, 6, 0, 14, 8, 0, 16, 20, 19, 12, 0, 18, 10, 0], [18, 9, 0, 13, 0, 14, 16, 0, 7, 17, 8, 20, 6, 10, 11, 15, 12, 0, 19, 0], [0, 8, 10, 9, 13, 0, 16, 0, 0, 18, 7, 14, 15, 12, 20, 6, 0, 19, 17, 11], [16, 7, 0, 11, 13, 15, 0, 9, 0, 14, 19, 12, 0, 18, 0, 20, 17, 8, 10, 6], [17, 12, 9, 14, 0, 19, 8, 0, 10, 0, 7, 16, 13, 11, 0, 0, 20, 18, 6, 15], [6, 18, 13, 17, 20, 0, 0, 7, 11, 14, 19, 0, 8, 0, 16, 15, 9, 10, 12, 0], [11, 9, 7, 0, 6, 0, 0, 12, 0, 0, 8, 19, 20, 14, 15, 17, 13, 10, 16, 18], [20, 0, 11, 6, 14, 17, 7, 13, 8, 9, 12, 0, 18, 10, 19, 0, 16, 0, 0, 15], [0, 15, 0, 17, 8, 13, 0, 18, 19, 16, 0, 7, 12, 14, 11, 10, 9, 0, 20, 6], [0, 7, 20, 11, 19, 0, 17, 0, 14, 15, 8, 12, 0, 6, 16, 10, 18, 0, 13, 9], [7, 14, 16, 18, 19, 9, 10, 0, 8, 12, 0, 6, 11, 15, 20, 13, 0, 17, 0, 0], [12, 8, 18, 13, 6, 11, 0, 10, 0, 9, 14, 19, 7, 16, 0, 17, 20, 0, 0, 15], [0, 18, 19, 7, 13, 12, 0, 6, 16, 15, 8, 11, 9, 17, 20, 10, 0, 0, 14, 0], [0, 20, 0, 8, 13, 7, 12, 15, 18, 9, 6, 0, 19, 10, 17, 11, 16, 0, 0, 14], [9, 19, 8, 14, 6, 18, 11, 7, 12, 0, 15, 17, 0, 0, 10, 0, 20, 0, 16, 13], [7, 0, 20, 11, 10, 12, 13, 14, 17, 15, 9, 8, 18, 6, 0, 0, 0, 19, 16, 0], [8, 0, 15, 14, 7, 11, 9, 16, 10, 6, 0, 13, 17, 20, 19, 18, 0, 0, 12, 0], [0, 8, 13, 17, 20, 0, 10, 7, 14, 6, 19, 15, 9, 12, 0, 11, 0, 16, 18, 0], [17, 11, 0, 9, 18, 10, 0, 19, 6, 7, 0, 0, 13, 20, 0, 16, 14, 12, 15, 8], [12, 0, 15, 20, 0, 0, 8, 10, 19, 9, 0, 7, 6, 14, 18, 11, 0, 13, 16, 17], [8, 19, 12, 11, 10, 0, 14, 13, 17, 15, 0, 0, 6, 20, 9, 0, 7, 18, 16, 0], [12, 9, 18, 10, 0, 14, 0, 13, 0, 0, 20, 17, 7, 0, 6, 16, 8, 11, 19, 15], [8, 0, 15, 19, 0, 0, 7, 20, 6, 11, 18, 0, 0, 9, 17, 13, 14, 12, 10, 16], [0, 11, 8, 0, 19, 17, 20, 15, 13, 14, 16, 10, 0, 0, 6, 9, 12, 18, 0, 7], [19, 20, 9, 0, 8, 10, 13, 0, 12, 6, 11, 7, 0, 0, 17, 16, 14, 0, 15, 18], [0, 18, 6, 13, 19, 15, 20, 0, 0, 17, 7, 12, 0, 10, 0, 11, 16, 8, 14, 9], [16, 19, 10, 0, 6, 9, 15, 7, 8, 18, 13, 14, 0, 0, 0, 0, 12, 20, 17, 11], [0, 0, 6, 12, 0, 8, 13, 18, 9, 11, 14, 17, 7, 15, 20, 0, 16, 10, 0, 19], [17, 14, 6, 16, 18, 13, 0, 12, 15, 19, 0, 7, 8, 9, 0, 0, 0, 20, 10, 11], [15, 10, 0, 11, 0, 0, 12, 8, 0, 9, 13, 17, 0, 20, 19, 7, 18, 14, 6, 16], [14, 20, 10, 16, 9, 13, 0, 0, 7, 18, 8, 0, 11, 15, 12, 6, 0, 17, 0, 19], [0, 0, 20, 9, 10, 11, 17, 15, 18, 7, 16, 6, 8, 0, 14, 19, 0, 12, 0, 13], [7, 18, 8, 6, 9, 16, 12, 0, 20, 11, 15, 19, 0, 13, 0, 10, 0, 0, 14, 17], [17, 0, 16, 10, 0, 11, 18, 0, 0, 8, 9, 19, 13, 15, 12, 14, 20, 7, 0, 6], [10, 11, 20, 0, 12, 13, 18, 9, 0, 17, 14, 0, 0, 7, 19, 8, 16, 0, 6, 15], [8, 9, 11, 18, 17, 16, 20, 0, 0, 10, 7, 6, 12, 15, 0, 14, 19, 0, 0, 13], [17, 0, 13, 6, 7, 0, 11, 15, 14, 8, 19, 0, 10, 16, 0, 9, 20, 12, 18, 0], [17, 12, 0, 19, 20, 14, 16, 6, 9, 18, 15, 13, 0, 11, 10, 0, 0, 7, 8, 0], [0, 0, 17, 6, 0, 11, 0, 12, 20, 8, 10, 14, 15, 16, 0, 18, 19, 13, 7, 9], [8, 18, 7, 0, 19, 16, 0, 20, 0, 11, 6, 13, 10, 0, 0, 9, 17, 14, 12, 15], [11, 14, 19, 7, 0, 16, 13, 12, 10, 6, 0, 20, 0, 18, 0, 15, 17, 0, 8, 9], [16, 7, 0, 19, 12, 10, 17, 15, 6, 14, 8, 0, 0, 11, 9, 13, 18, 0, 0, 20], [17, 0, 0, 12, 20, 14, 10, 0, 0, 7, 8, 0, 6, 16, 9, 19, 15, 11, 13, 18], [13, 20, 19, 9, 8, 7, 17, 0, 10, 12, 11, 0, 6, 18, 0, 0, 15, 0, 16, 14], [19, 9, 12, 0, 8, 0, 18, 0, 17, 16, 14, 6, 0, 11, 13, 7, 15, 10, 20, 0], [14, 18, 10, 8, 0, 19, 11, 20, 15, 12, 7, 6, 0, 0, 9, 0, 16, 17, 0, 13], [17, 0, 0, 0, 14, 0, 19, 11, 18, 9, 0, 15, 16, 12, 20, 6, 7, 10, 13, 8], [0, 7, 0, 8, 0, 9, 11, 19, 0, 15, 6, 17, 16, 20, 14, 10, 12, 18, 13, 0], [0, 0, 7, 0, 10, 16, 6, 8, 0, 14, 18, 15, 13, 20, 9, 12, 19, 17, 0, 11], [17, 20, 7, 14, 9, 0, 16, 12, 0, 0, 8, 11, 10, 19, 0, 15, 13, 18, 0, 6], [0, 18, 19, 7, 10, 12, 0, 11, 0, 13, 6, 20, 17, 14, 8, 15, 16, 0, 0, 9], [0, 0, 7, 17, 13, 8, 12, 14, 16, 0, 0, 18, 20, 19, 10, 6, 11, 9, 0, 15], [6, 13, 9, 10, 12, 7, 15, 20, 0, 18, 11, 8, 0, 0, 19, 14, 16, 0, 17, 0], [8, 20, 0, 17, 7, 16, 13, 0, 6, 0, 11, 15, 12, 14, 10, 0, 0, 9, 19, 18], [14, 0, 0, 9, 11, 17, 16, 19, 20, 0, 13, 7, 8, 0, 6, 0, 15, 10, 18, 12], [17, 0, 0, 16, 19, 6, 11, 13, 8, 0, 15, 0, 7, 14, 10, 12, 18, 20, 9, 0]]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "NUMBER OF THREADS BEING USED , DISPLAYING BEFORE SOLVE  0\n",
      "Version identifier: 12.10.0.0 | 2019-11-26 | 843d4de2ae\n",
      "CPXPARAM_Read_DataCheck                          1\n",
      "CPXPARAM_MIP_Strategy_File                       3\n",
      "CPXPARAM_TimeLimit                               14400\n",
      "CPXPARAM_MIP_Tolerances_MIPGap                   1.0000000000000001e-17\n",
      "Found incumbent of value 0.000000 after 0.00 sec. (0.06 ticks)\n",
      "Tried aggregator 1 time.\n",
      "MIP Presolve eliminated 355 rows and 118 columns.\n",
      "MIP Presolve modified 2043 coefficients.\n",
      "Reduced MIP has 1023 rows, 292 columns, and 3197 nonzeros.\n",
      "Reduced MIP has 141 binaries, 10 generals, 0 SOSs, and 0 indicators.\n",
      "Presolve time = 0.00 sec. (2.94 ticks)\n",
      "Probing time = 0.00 sec. (2.75 ticks)\n",
      "Tried aggregator 1 time.\n",
      "Detecting symmetries...\n",
      "MIP Presolve eliminated 124 rows and 48 columns.\n",
      "MIP Presolve modified 520 coefficients.\n",
      "Reduced MIP has 899 rows, 244 columns, and 2865 nonzeros.\n",
      "Reduced MIP has 121 binaries, 87 generals, 0 SOSs, and 0 indicators.\n",
      "Presolve time = 0.02 sec. (3.54 ticks)\n",
      "Probing time = 0.02 sec. (2.23 ticks)\n",
      "Clique table members: 763.\n",
      "MIP emphasis: balance optimality and feasibility.\n",
      "MIP search method: dynamic search.\n",
      "Parallel mode: deterministic, using up to 8 threads.\n",
      "Root relaxation solution time = 0.01 sec. (5.96 ticks)\n",
      "\n",
      "        Nodes                                         Cuts/\n",
      "   Node  Left     Objective  IInf  Best Integer    Best Bound    ItCnt     Gap\n",
      "\n",
      "*     0+    0                            0.0000     1647.0000              --- \n",
      "      0     0      231.0000    49        0.0000      231.0000      429     --- \n",
      "      0     0      231.0000    49        0.0000       Cuts: 8      447     --- \n",
      "      0     0      231.0000    49        0.0000      Cuts: 66      488     --- \n",
      "      0     0      231.0000    49        0.0000      Cuts: 23      523     --- \n",
      "*     0+    0                          225.0000      231.0000             2.67%\n",
      "*     0+    0                          227.0000      231.0000             1.76%\n",
      "      0     2      231.0000    17      227.0000      231.0000      523    1.76%\n",
      "Elapsed time = 0.28 sec. (72.73 ticks, tree = 0.02 MB, solutions = 3)\n",
      "*     4+    1                          228.0000      231.0000             1.32%\n",
      "*   654     8      integral     0      229.0000      230.0000    14187    0.44%\n",
      "\n",
      "Clique cuts applied:  4\n",
      "Implied bound cuts applied:  338\n",
      "Mixed integer rounding cuts applied:  3\n",
      "Zero-half cuts applied:  2\n",
      "Gomory fractional cuts applied:  4\n",
      "\n",
      "Root node processing (before b&c):\n",
      "  Real time             =    0.28 sec. (72.41 ticks)\n",
      "Parallel b&c, 8 threads:\n",
      "  Real time             =    0.25 sec. (95.83 ticks)\n",
      "  Sync time (average)   =    0.06 sec.\n",
      "  Wait time (average)   =    0.00 sec.\n",
      "                          ------------\n",
      "Total (root+branch&cut) =    0.53 sec. (168.25 ticks)\n",
      "Solve Details in order mip relative gap, best bound 0.0 229.0 integer optimal solution\n",
      "\n",
      " Gap is  0.0\n",
      "time limit =  14400.0\n",
      "threads =  0\n",
      "threads   0\n",
      "\n",
      " Upper Level Objective \n",
      " 229.0 229.0\n",
      "\n",
      " Lower Level Objective \n",
      " 63.4\n",
      "\n",
      " Assignment of products \n",
      " [0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0]\n",
      "Prices of products  [27.0, 19.0, 27.0, 22.0, 20.0, 24.0, 25.0, 26.0, 24.0, 23.0]\n",
      "Numbered Prod Purchased  [3, 8, 4, 7, 2, 5, 9, 5, 6, 4]\n"
     ]
    }
   ],
   "source": [
    "K,J , R_K, BUDGETS = read_instance()\n",
    "Cust_Ass, Opt_Prices  ,upper_objective ,upper_objective_kpi, lower_objective_kpi , Q_aux  = cplex_linearized ( K,J , R_K, BUDGETS , S_K)\n",
    "\n",
    "print('\\n Upper Level Objective \\n',  upper_objective , upper_objective_kpi  )\n",
    "print('\\n Lower Level Objective \\n',lower_objective_kpi )\n",
    "print('\\n Assignment of products \\n' , Cust_Ass  )\n",
    "print('Prices of products ',    ( [float(Decimal(\"%.2f\" % price)) for price in Opt_Prices] ) )\n",
    "print('Numbered Prod Purchased ',  prod_number(Cust_Ass , K , J ))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 112,
   "metadata": {},
   "outputs": [],
   "source": [
    "# len(Opt_Prices)  135807 97135     2071.0814   723     1961.0000     2078.3045  9811353    5.98%"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 95,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[4, 15, 9, 19, 16, 11, 19, 15, 14, 5, 6, 13, 3, 5, 15, 19, 8, 15, 17, 7], [6, 15, 15, 2, 17, 6, 5, 16, 16, 8, 5, 18, 20, 8, 10, 4, 14, 4, 4, 18], [27, 15, 23, 19, 23, 25, 13, 11, 19, 17, 16, 13, 21, 12, 16, 11, 18, 16, 27, 16], [5, 17, 21, 6, 14, 11, 5, 11, 4, 22, 4, 15, 11, 6, 22, 5, 11, 19, 12, 13], [19, 4, 9, 3, 9, 18, 5, 17, 11, 15, 8, 12, 13, 14, 15, 12, 9, 10, 18, 4], [9, 17, 10, 11, 20, 21, 19, 13, 26, 26, 24, 25, 19, 18, 12, 21, 21, 20, 23, 21], [13, 17, 20, 25, 8, 9, 11, 14, 22, 26, 10, 17, 13, 24, 20, 18, 13, 18, 14, 13], [19, 25, 20, 12, 11, 17, 18, 27, 26, 23, 14, 26, 9, 24, 15, 15, 19, 9, 10, 21], [20, 16, 15, 18, 9, 10, 18, 6, 20, 15, 16, 19, 10, 6, 23, 17, 16, 7, 20, 19], [23, 26, 19, 9, 21, 17, 15, 11, 15, 19, 18, 10, 27, 24, 27, 14, 11, 18, 19, 23], [21, 13, 11, 9, 17, 14, 21, 24, 23, 21, 17, 16, 17, 27, 8, 9, 22, 12, 10, 17], [7, 4, 4, 18, 16, 9, 15, 3, 13, 17, 13, 11, 13, 13, 6, 10, 13, 9, 8, 10], [21, 5, 8, 5, 17, 16, 12, 20, 14, 13, 14, 7, 7, 6, 21, 8, 9, 17, 17, 22], [14, 6, 19, 9, 16, 17, 11, 7, 18, 22, 5, 18, 13, 18, 19, 20, 11, 19, 5, 7], [10, 19, 14, 22, 21, 18, 13, 6, 16, 7, 8, 20, 13, 12, 18, 15, 6, 14, 15, 24], [16, 16, 13, 16, 23, 10, 8, 11, 20, 19, 24, 19, 12, 22, 12, 20, 15, 13, 12, 9], [22, 9, 8, 6, 21, 9, 14, 11, 4, 14, 15, 22, 18, 5, 13, 22, 18, 18, 9, 16], [13, 3, 5, 5, 1, 9, 13, 11, 14, 10, 9, 9, 18, 3, 16, 16, 13, 14, 17, 14], [17, 22, 8, 8, 8, 6, 8, 13, 9, 6, 21, 19, 18, 17, 18, 5, 14, 7, 12, 10], [11, 20, 23, 19, 17, 7, 7, 22, 11, 7, 20, 20, 20, 23, 13, 15, 10, 16, 23, 18]] \n",
      " [[17, 20, 18, 7, 6, 10, 8, 15, 16, 9, 13, 0, 19, 0, 12, 0, 14, 11, 0, 0], [18, 14, 11, 15, 6, 17, 12, 16, 7, 10, 19, 0, 9, 8, 0, 0, 13, 20, 0, 0], [16, 0, 20, 0, 8, 15, 13, 0, 6, 9, 12, 19, 7, 10, 0, 18, 17, 14, 11, 0], [13, 0, 17, 8, 0, 16, 6, 18, 12, 15, 0, 9, 14, 7, 10, 11, 19, 0, 20, 0], [6, 15, 0, 20, 19, 0, 18, 0, 11, 13, 14, 9, 12, 8, 17, 16, 7, 0, 0, 10], [12, 9, 14, 0, 11, 0, 0, 20, 16, 0, 7, 13, 17, 15, 8, 6, 18, 0, 19, 10], [20, 9, 0, 17, 16, 19, 15, 10, 7, 13, 0, 8, 14, 0, 6, 11, 18, 0, 0, 12], [17, 9, 11, 18, 0, 12, 14, 7, 0, 16, 6, 10, 0, 15, 0, 20, 13, 19, 0, 8], [6, 0, 17, 19, 12, 14, 0, 11, 7, 0, 15, 9, 13, 16, 0, 8, 0, 18, 20, 10], [13, 8, 11, 17, 0, 6, 15, 12, 7, 0, 16, 14, 20, 18, 9, 0, 0, 19, 10, 0], [10, 20, 0, 0, 14, 11, 6, 8, 15, 9, 0, 0, 7, 12, 0, 17, 13, 16, 18, 19], [0, 18, 11, 7, 12, 19, 17, 0, 6, 15, 0, 16, 14, 13, 10, 20, 9, 8, 0, 0], [9, 0, 18, 17, 0, 7, 10, 0, 8, 16, 14, 11, 0, 13, 19, 20, 15, 6, 0, 12], [20, 18, 16, 8, 0, 11, 15, 7, 14, 0, 6, 9, 19, 17, 13, 0, 12, 0, 0, 10], [18, 0, 15, 6, 19, 9, 14, 0, 11, 0, 16, 0, 17, 20, 7, 13, 12, 10, 8, 0], [13, 14, 16, 0, 0, 17, 19, 11, 6, 0, 8, 18, 9, 0, 20, 12, 15, 7, 10, 0], [15, 19, 11, 20, 0, 9, 0, 18, 0, 13, 0, 12, 7, 8, 17, 14, 6, 0, 16, 10], [0, 15, 10, 11, 18, 20, 0, 13, 14, 16, 9, 0, 0, 12, 17, 7, 6, 0, 8, 19], [0, 17, 20, 18, 8, 9, 10, 7, 11, 0, 14, 0, 15, 16, 0, 12, 0, 19, 13, 6], [7, 15, 20, 0, 9, 13, 10, 0, 0, 14, 6, 18, 12, 0, 0, 16, 11, 19, 17, 8]] \n",
      "\n",
      " [19, 20, 27, 22, 19, 26, 26, 27, 23, 27, 27, 18, 22, 22, 24, 24, 22, 18, 22, 23]\n"
     ]
    }
   ],
   "source": [
    "print(R_K,'\\n', S_K, '\\n\\n' , BUDGETS )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "# -------------CPLEX Reformulation ------------------\n",
    "def cplex_linearized( K,J ,  R_K, BUDGETS, S_K  , **kwargs):  \n",
    "     \n",
    "    #iterations +=1 \n",
    "    #     ----- Result Variable\n",
    "#     K = len(BUDGETS)\n",
    "#     J = len(BUDGETS)\n",
    "    Cust_Ass =[]              #Empty list to store decision variable x_i_j\n",
    "    Prices =[]\n",
    "    Q  =  []\n",
    "    lower_objective = 0\n",
    "    max_b = max(BUDGETS)\n",
    "    mdl = Model(name='lowerlevel' )\n",
    "    L = min(BUDGETS) \n",
    "    U = max(BUDGETS) \n",
    "    pi_min = min(min(R_K))\n",
    "    pi_max = max(BUDGETS)\n",
    "    epsilon=0.01\n",
    "    weight = 0.8\n",
    "    big_M = 1000\n",
    "   \n",
    "    J_range = range(0,J )\n",
    "    K_range = range(0,K)\n",
    "    # ------ decision variables ---\n",
    "    x  = mdl.binary_var_matrix(K,J ,\"cust_assign\", key_format =\"_%s\")\n",
    "    p  = mdl.continuous_var_matrix(K,J, name=\"price\"  , key_format =\"_%s\")\n",
    "    pi = mdl.integer_var_list( J , name=\"pi\" , key_format =\"_%s\")\n",
    "    q  = mdl.binary_var_matrix(K,J ,\"aux_q\", key_format =\"_%s\")\n",
    "    t  = mdl.continuous_var_matrix(K,J ,name=\"aux_t\", key_format =\"_%s\")\n",
    "\n",
    "    # --------constraints 1---\n",
    "    for k in K_range:  \n",
    "        mdl.add_constraint( mdl.sum(  x[k,j] for j in J_range if S_K[k][j] > 0 )  <= 1   )\n",
    "    \n",
    "    # --------constraints 2---\n",
    "#      Binary constraints \n",
    "\n",
    "    # --------- constraints3 ---\n",
    "    for k in K_range:\n",
    "        for j in J_range:\n",
    "            for i in J_range:\n",
    "                if (i !=j  and  S_K[k][j] > 0 and S_K[k][i] > 0    ):\n",
    "                    mdl.add_constraint(  ( weight*S_K[k][j] +  (1-weight) *R_K[k][j] ) *x[k,j ] - ( 1-weight)*p[k,j]  >= ( weight*S_K[k][i] + (1-weight)*R_K[k][i] )  *x[k,j] - (1-weight) * pi[i] - big_M * q[k,i]  )  \n",
    "   \n",
    "    # --------- constraints4 ---\n",
    "    for k in K_range:\n",
    "        for j in J_range:\n",
    "            if S_K[k][j] > 0:\n",
    "                mdl.add_constraint( p[k,j]  <= pi[j] )      \n",
    "            \n",
    "    # --------- constraints5 ---     \n",
    "    for k in K_range:\n",
    "        for j in J_range:\n",
    "            if S_K[k][j] > 0:\n",
    "                mdl.add_constraint( p[k,j]  <= BUDGETS[k] * x[k,j] )\n",
    "    \n",
    "     # --------- constraint redundant ---   FOR PROOF THIS CONSTRAINT BECOMES REDUNDANT   \n",
    "#     for k in R:\n",
    "#         if S_K[k][j] > 0:\n",
    "#               mdl.add_constraint( mdl.sum(   p[k,j]   for j in R  )  <= BUDGETS[k]    )\n",
    "                 \n",
    "    # --------- constraints6 ---\n",
    "    for k in K_range:\n",
    "        for j in J_range:\n",
    "            if S_K[k][j] > 0:\n",
    "                mdl.add_constraint( p[k,j]  >= pi[j] - max_b * (1-  x[k,j] )  )\n",
    "        \n",
    "#     # --------- constraints7 ------------  Auxiliary q  \n",
    "    for k in range(0,K):\n",
    "        for i in range(0,J):\n",
    "            if S_K[k][i] > 0:\n",
    "                mdl.add_constraint( BUDGETS[k]  >=  pi[i] -   t[k,i] + L * q[k,i]  )\n",
    "    \n",
    "    # --------- constraints8 ------------      \n",
    "    for k in range(0,K):\n",
    "        for i in range(0,J):\n",
    "            if S_K[k][i] > 0:\n",
    "                mdl.add_constraint(  BUDGETS[k]   <=  t[k,i] - epsilon * q[k,i]  + U*(1- q[k,i] )    )  \n",
    "                \n",
    "    # --------- constraints9 a------------      <=   pi_max * q[k,i] \n",
    "    for k in range(0,K):\n",
    "        for i in range(0,J):\n",
    "            if S_K[k][i] > 0:\n",
    "                mdl.add_constraint( pi_min* q[k,i] <= t[k,i]   )  \n",
    "    \n",
    "    # --------- constraints9 b ------------     \n",
    "    for k in range(0,K):\n",
    "        for i in range(0,J):\n",
    "            if S_K[k][i] > 0:\n",
    "                mdl.add_constraint(   t[k,i]   <=   pi_max * q[k,i]  )  \n",
    "    \n",
    "#     --------- constraints10_a ------------      <=   -pi_min * (1- q[k,i])\n",
    "    for k in range(0,K):\n",
    "        for i in range(0,J):\n",
    "            if S_K[k][i] > 0:\n",
    "                mdl.add_constraint( -pi_max* ( 1- q[k,i]) <= t[k,i] - pi[i]   )   \n",
    "    \n",
    "   # --------- constraints10 b ------------      \n",
    "    for k in range(0,K):\n",
    "        for i in range(0,J):\n",
    "            if S_K[k][i] > 0:\n",
    "                mdl.add_constraint(   t[k,i] - pi[i] <=   -pi_min * (1- q[k,i])   )   \n",
    "                \n",
    "    # --------- constraints11 ------------\n",
    "    for k in K_range:\n",
    "        for j in J_range :\n",
    "            if S_K[k][j] > 0:\n",
    "                mdl.add_constraint( p[k,j] >= 0   )   \n",
    "                \n",
    "    # --------- constraints11 ------------                                     \n",
    "    for j in J_range:\n",
    "        if S_K[k][j] > 0:\n",
    "            mdl.add_constraint( pi[j]    >= 0   )               \n",
    "    \n",
    "    # ------------ objective ---         \n",
    "    mdl.maximize( mdl.sum(   p[k,j]   for k in range(0,K) for j in range(0,J) if S_K[k][j] > 0  ) )        \n",
    "    \n",
    "    # ------------ KPI  ---   \n",
    "    mdl.add_kpi(  mdl.sum( x[k,j ]*pi[j]      for k in K_range for j in J_range )  , publish_name = 'UpperObjective')  \n",
    "\n",
    "    # ------------ KPI Lower Level ---   \n",
    "    mdl.add_kpi( mdl.sum(   weight*  x[k,j ] *S_K[k][j] +  x[k,j] *(1-weight)* ( R_K[k][j] - pi[j]  ) for k in range(0,K) for j in range(0,J)) , publish_name = 'LowerObjective') \n",
    "    \n",
    "    \n",
    "    #------------SOLVE SETTINGS --------------------------------\n",
    "    mdl.parameters.mip.tolerances.mipgap.set(0.00000000000000001)          #Assigns Gaps \n",
    "    mdl.set_time_limit(14400)                    #Setting Time Limit \n",
    "    mdl.parameters.mip.strategy.file = 3         # Tranferring Node FIle , mmeory full issue \n",
    "#     model.context.cplex_parameters.threads = 4   # Number of Threads . NOT USEDDD \n",
    "\n",
    "    print('NUMBER OF THREADS BEING USED , DISPLAYING BEFORE SOLVE ' , mdl.context.cplex_parameters.threads.get()   )\n",
    "    msol = mdl.solve(log_output=True)                 #log_output=True              # msol is a data structure of the class Solve Solution\n",
    "    \n",
    "    # ----- Retrieving values of decision variables Assignment \n",
    "    for i in range(0,K):\n",
    "            for j in range(0,J):\n",
    "                    ld = msol.get_value(x[i,j ])\n",
    "                    Cust_Ass.append(ld) \n",
    "    # ----- Retrieving values of decision variables PRICES\n",
    "    for j in range(0,J):\n",
    "        ld = msol.get_value(pi[ j ])\n",
    "        Prices.append(ld) \n",
    "    \n",
    "    # ----- Retrieving values of q\n",
    "    for k in range(0,K):\n",
    "            for i in range(0,J):\n",
    "                    ld = msol.get_value( q[k,i] )\n",
    "                    Q.append(ld) \n",
    "    \n",
    "    #----------------Solve module , solve detals ----------------\n",
    "    print('Solve Details in order mip relative gap, best bound' ,msol.solve_details.mip_relative_gap , msol.solve_details.best_bound , msol.solve_details.status )\n",
    "    print('\\n Gap is ', msol.solve_details.gap)\n",
    "    \n",
    "    # Some Prints \n",
    "    mdl.export_as_lp(basename=\"cplex_%s\", path=\"D:\\AA MASC\\C Codes\\FINAL\" )  #Printing the constraints \n",
    "    print(\"time limit = \",mdl.get_time_limit()) #Prints time limit \n",
    "    print(\"threads = \",mdl.parameters.threads.get()) #Prints thread\n",
    "    print('threads  ' , mdl.context.cplex_parameters.threads.get()   ) \n",
    "#     time_lim =  mdl.get_time_limit()           #getting the time limit \n",
    "#     mdl.report()                  #Prints KPIs  , value of obj\n",
    "#     mdl.print_solution() \n",
    "    upper_objective = msol.get_objective_value()                # Retrieves obj value \n",
    "    upper_objective_kpi = msol.kpi_value_by_name('UpperObjective')\n",
    "    lower_objective_kpi = msol.kpi_value_by_name('LowerObjective')\n",
    "    \n",
    "    return Cust_Ass, Prices  ,upper_objective ,upper_objective_kpi, lower_objective_kpi , Q"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  READ AT SOLVE \n",
    "#  K,J , R_K, BUDGETS = read_instance()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "#-----------DONT RUN  (RUN ONLY ONCE )----------------------------\n",
    "def Rank_Generator(K,J):\n",
    "    Ranked_List = [ [0 for x in range(J)] for x in range(K) ] \n",
    "    Bool_list =  [ [True for x in range(J)] for x in range(K) ] \n",
    "    \n",
    "    for k in range(0,K):\n",
    "#         print('Customer  is ' ,k )\n",
    "        Rank_value = J\n",
    "        iterator = 0\n",
    "        Nested_Rank = []\n",
    "        preferred_size_not_obtained =  math.ceil (0.75 * J ) \n",
    "#         print(preferred_size_not_obtained)\n",
    "        while preferred_size_not_obtained:\n",
    "            random_index =  random.randint(0, J-1 )\n",
    "#             print(random_index)\n",
    "            if Bool_list[k][random_index]:\n",
    "                Ranked_List[k][random_index] = Rank_value\n",
    "#                 print('Value updated at position ', Ranked_List[k])\n",
    "                Bool_list[k][random_index] = False\n",
    "                Rank_value -= 1 \n",
    "                preferred_size_not_obtained -=1\n",
    "#                 print('Updated size',preferred_size_not_obtained)\n",
    "                 \n",
    "\n",
    "    return Ranked_List\n",
    "\n",
    "Ranked_List = Rank_Generator(K,J)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "S_K = Ranked_List"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "# print(S_K)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [],
   "source": [
    "# COUNT NUMBER OF ZEROES S\n",
    "# def ranked_list_check(Ranked_List)\n",
    "# for indv_list in Ranked_List:\n",
    "#     print(indv_list.count(0))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "# ----------Funtion to Compute , Assignment \n",
    "def prod_number(cust_list, no_cust, no_prod  ):\n",
    "    numbered_list = []\n",
    "    iterator = 0\n",
    "    for i in range(no_cust):\n",
    "        flag = 0\n",
    "        j = 0\n",
    "        for j in range(no_prod ):  \n",
    "#             print(int(str(i) + str(j)))\n",
    "#             if ( cust_list[int(str(i) + str(j))]  >  0 ) :  \n",
    "            if ( cust_list[iterator]  >  0 ) :  \n",
    "                numbered_list.append(j+1)\n",
    "                flag=1\n",
    "            j += 1\n",
    "            iterator += 1\n",
    "        if flag ==0:\n",
    "            numbered_list.append(0)\n",
    "            \n",
    "    return numbered_list\n",
    "\n",
    "# This list contains actual index i.e. Product1, Prod 2 NOT python indices from 0,\n",
    "# A value of 0 indicates customer did not made any purchase "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Numbered Prod Purchased by checking   \t  [3, 8, 4, 7, 2, 5, 9, 5, 6, 4]\n",
      "Result Obtained from Reformulation   \t [3, 8, 4, 7, 2, 5, 9, 5, 6, 4]\n",
      "Lower Level is Optimal \n"
     ]
    }
   ],
   "source": [
    "# ----------Funtion to Check Optimality of Lower Level \n",
    "# -------------CPLEX LOWER LEVEL ------------------\n",
    "def master_lower_level( K, J , PRICES, R_K, BUDGETS, S_K  , **kwargs):  \n",
    "   \n",
    "     \n",
    "    #iterations +=1 \n",
    "    #     ----- Result Variable (Variables to store) ----------------\n",
    "    Cust_Ass =[]\n",
    "    lower_objective = 0\n",
    "    \n",
    "    mdl = Model(name='lowerlevel' )\n",
    "\n",
    "    weight =0.8          # Weight \n",
    "   \n",
    "    R = range(0,K)\n",
    "    J_Range = range(0,J )\n",
    "    # ------ decision variables ---\n",
    "    x  = mdl.binary_var_matrix(K,J ,\"cust_assign\", key_format =\"_%s\")\n",
    "     \n",
    "    # --------constraints 1---\n",
    "    for i in R:  \n",
    "        mdl.add_constraint( mdl.sum(  x[i,j] for j in J_Range )  <= 1   )\n",
    "\n",
    "    # --------- constraints2 ---\n",
    "    for i in R:\n",
    "        mdl.add_constraint( mdl.sum(   x[i,j]* PRICES[j] for j in J_Range  )  <= BUDGETS[ i]    )\n",
    "      \n",
    "    # ------------ objective ---         \n",
    "    mdl.maximize( mdl.sum(   weight*  x[i,j ] *S_K[i][j] +  x[i,j ] *(1-weight)* (R_K[i][j] - PRICES[j]  )     for i in R for j in J_Range) )      \n",
    "    \n",
    "    # ------------ KPI  ---   \n",
    "    mdl.add_kpi(  mdl.sum( x[i,j ]*PRICES[j]      for i in R for j in J_Range )  , publish_name = 'UpperObjective')  \n",
    "\n",
    " \n",
    "    msol = mdl.solve()                 #log_output=True              # msol is a data structure of the class Solve Solution\n",
    "    # ----- Retrieving values of decision variables \n",
    "    for i in range(0,K):\n",
    "            for j in J_Range:\n",
    "                    ld = msol.get_value(x[i,j ])\n",
    "                    Cust_Ass.append(ld) \n",
    "                    \n",
    "    lower_objective = msol.get_objective_value()                # Retrieves obj value \n",
    "    upper_objective = msol.kpi_value_by_name('UpperObjective')\n",
    "    \n",
    "    return Cust_Ass , lower_objective ,upper_objective\n",
    "\n",
    "## VALIDATING PRICES FROM CPLEX \n",
    "soln = master_lower_level(K, J ,  Opt_Prices, R_K, BUDGETS, S_K  ) \n",
    "print('Numbered Prod Purchased by checking   \\t ',  prod_number(soln[0], K , J ))\n",
    "print('Result Obtained from Reformulation   \\t', prod_number(Cust_Ass,K,J) )\n",
    "if (soln[0] == Cust_Ass):\n",
    "#     print(soln[0] , Cust_Ass)\n",
    "    print('Lower Level is Optimal ')\n",
    "else:\n",
    "#     print(soln[0] , Cust_Ass)\n",
    "    print('Not Optimal')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "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>0</th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>7</th>\n",
       "      <th>8</th>\n",
       "      <th>9</th>\n",
       "      <th>Budgets</th>\n",
       "      <th>Selection</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>7.2</td>\n",
       "      <td>5.8</td>\n",
       "      <td>7.4</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>5.8</td>\n",
       "      <td>2.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>0.8</td>\n",
       "      <td>1.8</td>\n",
       "      <td>27.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-1.0</td>\n",
       "      <td>5.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>5.2</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>6.4</td>\n",
       "      <td>8.0</td>\n",
       "      <td>0.6</td>\n",
       "      <td>3.6</td>\n",
       "      <td>26.0</td>\n",
       "      <td>8.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>-1.0</td>\n",
       "      <td>0.2</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>8.0</td>\n",
       "      <td>6.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>4.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>5.2</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>3.2</td>\n",
       "      <td>5.8</td>\n",
       "      <td>6.0</td>\n",
       "      <td>3.6</td>\n",
       "      <td>6.0</td>\n",
       "      <td>1.2</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>2.4</td>\n",
       "      <td>27.0</td>\n",
       "      <td>7.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>-1.0</td>\n",
       "      <td>7.8</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>19.0</td>\n",
       "      <td>2.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>3.8</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>5.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>-1.0</td>\n",
       "      <td>4.6</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>4.2</td>\n",
       "      <td>4.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>6.4</td>\n",
       "      <td>4.8</td>\n",
       "      <td>24.0</td>\n",
       "      <td>9.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>3.8</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>5.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>1.6</td>\n",
       "      <td>5.6</td>\n",
       "      <td>5.8</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>0.6</td>\n",
       "      <td>5.8</td>\n",
       "      <td>24.0</td>\n",
       "      <td>6.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>-1.0</td>\n",
       "      <td>2.2</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>6.4</td>\n",
       "      <td>6.2</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>-1.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>4.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Prices</th>\n",
       "      <td>27.0</td>\n",
       "      <td>19.0</td>\n",
       "      <td>27.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>25.0</td>\n",
       "      <td>26.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>23.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           0     1     2     3     4     5     6     7     8     9  Budgets  \\\n",
       "0        7.2   5.8   7.4  -1.0  -1.0   5.8   2.0   5.0   0.8   1.8     27.0   \n",
       "1       -1.0   5.0  -1.0  -1.0   5.2  -1.0   6.4   8.0   0.6   3.6     26.0   \n",
       "2       -1.0   0.2  -1.0   8.0   6.0  -1.0  -1.0  -1.0  -1.0  -1.0     22.0   \n",
       "3        5.2  -1.0   3.2   5.8   6.0   3.6   6.0   1.2  -1.0   2.4     27.0   \n",
       "4       -1.0   7.8  -1.0  -1.0  -1.0  -1.0  -1.0  -1.0  -1.0  -1.0     19.0   \n",
       "5       -1.0  -1.0  -1.0  -1.0   3.8  -1.0  -1.0  -1.0  -1.0  -1.0     20.0   \n",
       "6       -1.0   4.6  -1.0   4.2   4.0  -1.0  -1.0  -1.0   6.4   4.8     24.0   \n",
       "7       -1.0  -1.0  -1.0  -1.0   3.8  -1.0  -1.0  -1.0  -1.0  -1.0     20.0   \n",
       "8       -1.0  -1.0  -1.0   1.6   5.6   5.8  -1.0  -1.0   0.6   5.8     24.0   \n",
       "9       -1.0   2.2  -1.0   6.4   6.2  -1.0  -1.0  -1.0  -1.0  -1.0     22.0   \n",
       "Prices  27.0  19.0  27.0  22.0  20.0  24.0  25.0  26.0  24.0  23.0      0.0   \n",
       "\n",
       "        Selection  \n",
       "0             3.0  \n",
       "1             8.0  \n",
       "2             4.0  \n",
       "3             7.0  \n",
       "4             2.0  \n",
       "5             5.0  \n",
       "6             9.0  \n",
       "7             5.0  \n",
       "8             6.0  \n",
       "9             4.0  \n",
       "Prices        0.0  "
      ]
     },
     "execution_count": 35,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# --------Visualizing and Comparing the Result ---------\n",
    "\n",
    "def obj_evaluation(K,J, PRICES, BUDGETS, R_K, S_K):\n",
    "    Obj_list = []\n",
    "    weight = 0.8\n",
    "    for i in range(K):         # For each customer \n",
    "#         print(\"Customer is %d\"%(i+1))\n",
    "        Nest_O = []\n",
    "        for j in range(J):\n",
    "#             print(\"Product is %d\"%(j+1))\n",
    "#             print(\"Reservation: %d and Rank: %d\"%( R_K[i][j], S_K[i][j]))                   or BUDGETS[i] < PRICES[j]\n",
    "            if (S_K[i][j] == 0 or PRICES[j]>BUDGETS[i] ):\n",
    "#                 print(\"Customer %d , is unable to purchase Prod:%d as Budget is %d and Price is %d\"%( i+1,j+1, BUDGETS[i], PRICES[j] ))\n",
    "                obj = -1 \n",
    "            else:\n",
    "#                 obj = round( weight*S_K[i][j] + (1-weight)* ( R_K[i][j] - PRICES[j]), 2)  \n",
    "                obj = weight*S_K[i][j] + (1-weight)* ( R_K[i][j] - PRICES[j])      \n",
    "            Nest_O.append(obj)\n",
    "        Obj_list.append(Nest_O)        \n",
    "    return Obj_list\n",
    "\n",
    "new_list = obj_evaluation(K,J, Opt_Prices ,BUDGETS,  R_K, S_K)\n",
    "\n",
    "columns=[\"Cust1\", \"Cust2\",\"Cust3\",\"Cust4\",\"Cust5\",\"Cust6\",\"Cust7\", \"Cust8\", \"Cust9\", \"Cust10\"]\n",
    "products = [\"Prod1\",\"Prod2\",\"Prod3\",\"Prod4\",\"Prod5\",\"Prod6\",\"Prod7\",\"Prod8\",\"Prod9\",\"Prod10\" ]\n",
    "# , columns=products, index= columns  \n",
    "df =pd.DataFrame(new_list   )\n",
    "df['Budgets'] = BUDGETS\n",
    "df['Selection'] = prod_number(Cust_Ass, K , J )\n",
    "pric =[]\n",
    "pric = copy.deepcopy(Opt_Prices)\n",
    "pric.append(0)\n",
    "pric.append(0)\n",
    "df.loc['Prices'] = pric\n",
    "\n",
    "df_RK = pd.DataFrame(R_K, columns=products, index= columns  )\n",
    "df_RK['Budgets'] = BUDGETS\n",
    "df_SK = pd.DataFrame(S_K, columns=products ,index= columns )\n",
    "df_concat = pd.concat([df_RK,df_SK] , axis = 1)\n",
    "\n",
    "df "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "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>Prod1</th>\n",
       "      <th>Prod2</th>\n",
       "      <th>Prod3</th>\n",
       "      <th>Prod4</th>\n",
       "      <th>Prod5</th>\n",
       "      <th>Prod6</th>\n",
       "      <th>Prod7</th>\n",
       "      <th>Prod8</th>\n",
       "      <th>Prod9</th>\n",
       "      <th>Prod10</th>\n",
       "      <th>...</th>\n",
       "      <th>Prod1</th>\n",
       "      <th>Prod2</th>\n",
       "      <th>Prod3</th>\n",
       "      <th>Prod4</th>\n",
       "      <th>Prod5</th>\n",
       "      <th>Prod6</th>\n",
       "      <th>Prod7</th>\n",
       "      <th>Prod8</th>\n",
       "      <th>Prod9</th>\n",
       "      <th>Prod10</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Cust1</th>\n",
       "      <td>27</td>\n",
       "      <td>24</td>\n",
       "      <td>24</td>\n",
       "      <td>19</td>\n",
       "      <td>21</td>\n",
       "      <td>25</td>\n",
       "      <td>15</td>\n",
       "      <td>19</td>\n",
       "      <td>16</td>\n",
       "      <td>16</td>\n",
       "      <td>...</td>\n",
       "      <td>9</td>\n",
       "      <td>6</td>\n",
       "      <td>10</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7</td>\n",
       "      <td>5</td>\n",
       "      <td>8</td>\n",
       "      <td>3</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust2</th>\n",
       "      <td>22</td>\n",
       "      <td>20</td>\n",
       "      <td>22</td>\n",
       "      <td>15</td>\n",
       "      <td>14</td>\n",
       "      <td>25</td>\n",
       "      <td>21</td>\n",
       "      <td>26</td>\n",
       "      <td>15</td>\n",
       "      <td>25</td>\n",
       "      <td>...</td>\n",
       "      <td>5</td>\n",
       "      <td>6</td>\n",
       "      <td>7</td>\n",
       "      <td>0</td>\n",
       "      <td>8</td>\n",
       "      <td>0</td>\n",
       "      <td>9</td>\n",
       "      <td>10</td>\n",
       "      <td>3</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust3</th>\n",
       "      <td>13</td>\n",
       "      <td>8</td>\n",
       "      <td>18</td>\n",
       "      <td>22</td>\n",
       "      <td>22</td>\n",
       "      <td>17</td>\n",
       "      <td>21</td>\n",
       "      <td>8</td>\n",
       "      <td>16</td>\n",
       "      <td>7</td>\n",
       "      <td>...</td>\n",
       "      <td>6</td>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>7</td>\n",
       "      <td>9</td>\n",
       "      <td>8</td>\n",
       "      <td>4</td>\n",
       "      <td>5</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust4</th>\n",
       "      <td>21</td>\n",
       "      <td>18</td>\n",
       "      <td>19</td>\n",
       "      <td>11</td>\n",
       "      <td>14</td>\n",
       "      <td>22</td>\n",
       "      <td>27</td>\n",
       "      <td>20</td>\n",
       "      <td>13</td>\n",
       "      <td>19</td>\n",
       "      <td>...</td>\n",
       "      <td>8</td>\n",
       "      <td>0</td>\n",
       "      <td>6</td>\n",
       "      <td>10</td>\n",
       "      <td>9</td>\n",
       "      <td>5</td>\n",
       "      <td>7</td>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust5</th>\n",
       "      <td>6</td>\n",
       "      <td>18</td>\n",
       "      <td>5</td>\n",
       "      <td>2</td>\n",
       "      <td>18</td>\n",
       "      <td>4</td>\n",
       "      <td>18</td>\n",
       "      <td>19</td>\n",
       "      <td>14</td>\n",
       "      <td>16</td>\n",
       "      <td>...</td>\n",
       "      <td>4</td>\n",
       "      <td>10</td>\n",
       "      <td>0</td>\n",
       "      <td>8</td>\n",
       "      <td>0</td>\n",
       "      <td>7</td>\n",
       "      <td>5</td>\n",
       "      <td>9</td>\n",
       "      <td>6</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust6</th>\n",
       "      <td>4</td>\n",
       "      <td>3</td>\n",
       "      <td>7</td>\n",
       "      <td>13</td>\n",
       "      <td>11</td>\n",
       "      <td>14</td>\n",
       "      <td>12</td>\n",
       "      <td>15</td>\n",
       "      <td>8</td>\n",
       "      <td>20</td>\n",
       "      <td>...</td>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>5</td>\n",
       "      <td>6</td>\n",
       "      <td>7</td>\n",
       "      <td>8</td>\n",
       "      <td>9</td>\n",
       "      <td>4</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust7</th>\n",
       "      <td>24</td>\n",
       "      <td>6</td>\n",
       "      <td>17</td>\n",
       "      <td>19</td>\n",
       "      <td>12</td>\n",
       "      <td>10</td>\n",
       "      <td>18</td>\n",
       "      <td>15</td>\n",
       "      <td>16</td>\n",
       "      <td>15</td>\n",
       "      <td>...</td>\n",
       "      <td>3</td>\n",
       "      <td>9</td>\n",
       "      <td>4</td>\n",
       "      <td>6</td>\n",
       "      <td>7</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>5</td>\n",
       "      <td>10</td>\n",
       "      <td>8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust8</th>\n",
       "      <td>8</td>\n",
       "      <td>5</td>\n",
       "      <td>20</td>\n",
       "      <td>17</td>\n",
       "      <td>19</td>\n",
       "      <td>16</td>\n",
       "      <td>3</td>\n",
       "      <td>6</td>\n",
       "      <td>20</td>\n",
       "      <td>18</td>\n",
       "      <td>...</td>\n",
       "      <td>10</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4</td>\n",
       "      <td>5</td>\n",
       "      <td>9</td>\n",
       "      <td>7</td>\n",
       "      <td>8</td>\n",
       "      <td>6</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust9</th>\n",
       "      <td>14</td>\n",
       "      <td>15</td>\n",
       "      <td>24</td>\n",
       "      <td>10</td>\n",
       "      <td>20</td>\n",
       "      <td>21</td>\n",
       "      <td>23</td>\n",
       "      <td>21</td>\n",
       "      <td>11</td>\n",
       "      <td>16</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>5</td>\n",
       "      <td>7</td>\n",
       "      <td>8</td>\n",
       "      <td>3</td>\n",
       "      <td>6</td>\n",
       "      <td>4</td>\n",
       "      <td>9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust10</th>\n",
       "      <td>18</td>\n",
       "      <td>18</td>\n",
       "      <td>12</td>\n",
       "      <td>22</td>\n",
       "      <td>15</td>\n",
       "      <td>14</td>\n",
       "      <td>5</td>\n",
       "      <td>6</td>\n",
       "      <td>8</td>\n",
       "      <td>17</td>\n",
       "      <td>...</td>\n",
       "      <td>0</td>\n",
       "      <td>3</td>\n",
       "      <td>0</td>\n",
       "      <td>8</td>\n",
       "      <td>9</td>\n",
       "      <td>4</td>\n",
       "      <td>10</td>\n",
       "      <td>6</td>\n",
       "      <td>7</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>10 rows × 21 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        Prod1  Prod2  Prod3  Prod4  Prod5  Prod6  Prod7  Prod8  Prod9  Prod10  \\\n",
       "Cust1      27     24     24     19     21     25     15     19     16      16   \n",
       "Cust2      22     20     22     15     14     25     21     26     15      25   \n",
       "Cust3      13      8     18     22     22     17     21      8     16       7   \n",
       "Cust4      21     18     19     11     14     22     27     20     13      19   \n",
       "Cust5       6     18      5      2     18      4     18     19     14      16   \n",
       "Cust6       4      3      7     13     11     14     12     15      8      20   \n",
       "Cust7      24      6     17     19     12     10     18     15     16      15   \n",
       "Cust8       8      5     20     17     19     16      3      6     20      18   \n",
       "Cust9      14     15     24     10     20     21     23     21     11      16   \n",
       "Cust10     18     18     12     22     15     14      5      6      8      17   \n",
       "\n",
       "        ...  Prod1  Prod2  Prod3  Prod4  Prod5  Prod6  Prod7  Prod8  Prod9  \\\n",
       "Cust1   ...      9      6     10      0      0      7      5      8      3   \n",
       "Cust2   ...      5      6      7      0      8      0      9     10      3   \n",
       "Cust3   ...      6      3      0     10      7      9      8      4      5   \n",
       "Cust4   ...      8      0      6     10      9      5      7      3      0   \n",
       "Cust5   ...      4     10      0      8      0      7      5      9      6   \n",
       "Cust6   ...      3      0      5      6      7      8      9      4      0   \n",
       "Cust7   ...      3      9      4      6      7      0      0      5     10   \n",
       "Cust8   ...     10      0      0      4      5      9      7      8      6   \n",
       "Cust9   ...      0      0     10      5      7      8      3      6      4   \n",
       "Cust10  ...      0      3      0      8      9      4     10      6      7   \n",
       "\n",
       "        Prod10  \n",
       "Cust1        4  \n",
       "Cust2        4  \n",
       "Cust3        0  \n",
       "Cust4        4  \n",
       "Cust5        3  \n",
       "Cust6       10  \n",
       "Cust7        8  \n",
       "Cust8        3  \n",
       "Cust9        9  \n",
       "Cust10       5  \n",
       "\n",
       "[10 rows x 21 columns]"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_concat"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "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>Prod1</th>\n",
       "      <th>Prod2</th>\n",
       "      <th>Prod3</th>\n",
       "      <th>Prod4</th>\n",
       "      <th>Prod5</th>\n",
       "      <th>Prod6</th>\n",
       "      <th>Prod7</th>\n",
       "      <th>Prod8</th>\n",
       "      <th>Prod9</th>\n",
       "      <th>Prod10</th>\n",
       "      <th>Budgets</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Cust1</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>27.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust2</th>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>26.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust3</th>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>22.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust4</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>27.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust5</th>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>19.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust6</th>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>20.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust7</th>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>24.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust8</th>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>20.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust9</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>24.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Cust10</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>22.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Prices</th>\n",
       "      <td>27.0</td>\n",
       "      <td>19.0</td>\n",
       "      <td>27.0</td>\n",
       "      <td>22.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>25.0</td>\n",
       "      <td>26.0</td>\n",
       "      <td>24.0</td>\n",
       "      <td>23.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Prod1  Prod2  Prod3  Prod4  Prod5  Prod6  Prod7  Prod8  Prod9  Prod10  \\\n",
       "Cust1     0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0     0.0   \n",
       "Cust2     1.0    0.0    1.0    0.0    0.0    0.0    0.0    0.0    0.0     0.0   \n",
       "Cust3     1.0    0.0    0.0    0.0    0.0    1.0    1.0    1.0    1.0     0.0   \n",
       "Cust4     0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0    0.0     0.0   \n",
       "Cust5     1.0    0.0    0.0    1.0    0.0    1.0    1.0    1.0    1.0     1.0   \n",
       "Cust6     1.0    0.0    1.0    1.0    0.0    1.0    1.0    1.0    0.0     1.0   \n",
       "Cust7     1.0    0.0    1.0    0.0    0.0    0.0    0.0    1.0    0.0     0.0   \n",
       "Cust8     1.0    0.0    0.0    1.0    0.0    1.0    1.0    1.0    1.0     1.0   \n",
       "Cust9     0.0    0.0    1.0    0.0    0.0    0.0    1.0    1.0    0.0     0.0   \n",
       "Cust10    0.0    0.0    0.0    0.0    0.0    1.0    1.0    1.0    1.0     1.0   \n",
       "Prices   27.0   19.0   27.0   22.0   20.0   24.0   25.0   26.0   24.0    23.0   \n",
       "\n",
       "        Budgets  \n",
       "Cust1      27.0  \n",
       "Cust2      26.0  \n",
       "Cust3      22.0  \n",
       "Cust4      27.0  \n",
       "Cust5      19.0  \n",
       "Cust6      20.0  \n",
       "Cust7      24.0  \n",
       "Cust8      20.0  \n",
       "Cust9      24.0  \n",
       "Cust10     22.0  \n",
       "Prices      0.0  "
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# --------Visualizing Aux Q  ---------\n",
    "def q_aux_nest( Q_aux , K,J , PRICES, BUDGETS):\n",
    "    List_q1 = []\n",
    "    for i in range(K):         # For each customer \n",
    "        Nest_q= []\n",
    "        for j in range(J):\n",
    "            val = Q_aux[int(str(i)+str(j))]\n",
    "            Nest_q.append(val)\n",
    "        List_q1.append(Nest_q)\n",
    "    return List_q1\n",
    "\n",
    "Q_aux_nest =  q_aux_nest ( Q_aux, K,J ,  Opt_Prices ,BUDGETS )\n",
    "columns=[\"Cust1\", \"Cust2\",\"Cust3\",\"Cust4\",\"Cust5\",\"Cust6\",\"Cust7\", \"Cust8\", \"Cust9\", \"Cust10\"]\n",
    "products = [\"Prod1\",\"Prod2\",\"Prod3\",\"Prod4\",\"Prod5\",\"Prod6\",\"Prod7\",\"Prod8\",\"Prod9\",\"Prod10\" ]\n",
    "df_aux =pd.DataFrame(  Q_aux_nest, columns=products, index= columns   )\n",
    "df_aux['Budgets'] = BUDGETS\n",
    "pric =[]\n",
    "pric = copy.deepcopy(Opt_Prices)\n",
    "pric.append(0)\n",
    "df_aux.loc['Prices'] = pric\n",
    "df_aux"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "###########SOLUTIONS OF INSTNACE SET 4 "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  INSTANCE 1 \n",
    "\n",
    "167525 114661     2040.0605   612     1935.0000     2094.0000 34757959    8.22%\n",
    " 167967 115030     2089.3367   870     1935.0000     2094.0000 34862143    8.22%\n",
    " 168363 115510     1950.0670   669     1935.0000     2094.0000 35009090    8.22%\n",
    "\n",
    "GUB cover cuts applied:  15\n",
    "Clique cuts applied:  1562\n",
    "Cover cuts applied:  1250\n",
    "Implied bound cuts applied:  3679\n",
    "Flow cuts applied:  387\n",
    "Mixed integer rounding cuts applied:  2296\n",
    "Zero-half cuts applied:  45\n",
    "Lift and project cuts applied:  1\n",
    "Gomory fractional cuts applied:  53\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =   10.09 sec. (4695.74 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14390.25 sec. (2099875.92 ticks)\n",
    "  Sync time (average)   = 1540.60 sec.\n",
    "  Wait time (average)   =    0.09 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14400.34 sec. (2104571.66 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.08217054263565396 2093.9999999999986 time limit exceeded\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1935.0 1935.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1168.9999999999995\n",
    "\n",
    "Prices of products  [26.0, 27.0, 23.0, 24.0, 27.0, 25.0, 20.0, 29.0, 24.0, 20.0, 17.0, 22.0, 27.0, 24.0, 24.0, 23.0, 28.0, 27.0, 28.0, 26.0]\n",
    "Numbered Prod Purchased  [18, 18, 11, 4, 10, 12, 10, 7, 7, 8, 4, 7, 9, 9, 12, 11, 14, 14, 11, 3, 7, 7, 7, 12, 12, 15, 3, 10, 11, 4, 6, 3, 13, 4, 12, 11, 7, 7, 11, 7, 16, 0, 11, 4, 4, 18, 15, 11, 10, 10, 1, 14, 9, 11, 12, 10, 11, 7, 10, 7, 19, 4, 3, 11, 5, 6, 7, 6, 12, 12, 20, 11, 11, 11, 15, 7, 7, 11, 10, 10, 4, 16, 10, 13, 16, 2, 17, 10, 11, 11]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  INSTANCE 2   "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  INSTANCE 3 , SOLVED\n",
    "\n",
    "Elapsed time = 13940.88 sec. (1181341.24 ticks, tree = 23648.24 MB, solutions = 64)\n",
    "Nodefile size = 21575.19 MB (19919.69 MB after compression)\n",
    " 205878 154913     1982.5329   434     1961.0000     2077.9770 14754464    5.97%\n",
    " 206458 155093     2038.0007   739     1961.0000     2077.9767 14782063    5.97%\n",
    " 207180 155604     1992.8674   492     1961.0000     2077.9752 14830572    5.97%\n",
    " 208074 155772     1992.3905   479     1961.0000     2077.9735 14849839    5.96%\n",
    " 209096 156222     1970.7508   475     1961.0000     2077.9726 14877072    5.96%\n",
    " 209628 156806     2009.5564   541     1961.0000     2077.9702 14920889    5.96%\n",
    " 209958 158286     2015.7485   397     1961.0000     2077.9677 15037273    5.96%\n",
    " 210566 158455     1966.1964   307     1961.0000     2077.9660 15046757    5.96%\n",
    " 211492 158543     2074.9803   772     1961.0000     2077.9660 15065010    5.96%\n",
    " 212536 159210     2061.2869   544     1961.0000     2077.9560 15127402    5.96%\n",
    "Elapsed time = 14383.34 sec. (1219547.83 ticks, tree = 24408.26 MB, solutions = 64)\n",
    "Nodefile size = 22332.84 MB (20614.55 MB after compression)\n",
    " 213354 159772     1970.2249   381     1961.0000     2077.9546 15192710    5.96%\n",
    "\n",
    "GUB cover cuts applied:  4\n",
    "Clique cuts applied:  1145\n",
    "Cover cuts applied:  957\n",
    "Implied bound cuts applied:  3858\n",
    "Flow cuts applied:  184\n",
    "Mixed integer rounding cuts applied:  824\n",
    "Zero-half cuts applied:  39\n",
    "Lift and project cuts applied:  4\n",
    "Gomory fractional cuts applied:  77\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =    7.64 sec. (3055.95 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14404.36 sec. (1222658.42 ticks)\n",
    "  Sync time (average)   = 4452.59 sec.\n",
    "  Wait time (average)   =    1.62 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14412.00 sec. (1225714.37 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.059635380667564276 2077.9449814890995 time limit exceeded\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1961.0 1961.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1176.8000000000006\n",
    "\n",
    " Assignment of products \n",
    " [0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0]\n",
    "Prices of products  [22.0, 26.0, 22.0, 21.0, 25.0, 25.0, 23.0, 27.0, 26.0, 28.0, 19.0, 23.0, 21.0, 27.0, 25.0, 28.0, 22.0, 23.0, 27.0, 17.0]\n",
    "Numbered Prod Purchased  [4, 11, 4, 20, 1, 10, 18, 20, 17, 20, 13, 20, 15, 19, 1, 4, 19, 12, 5, 17, 16, 11, 15, 20, 18, 18, 4, 13, 12, 11, 3, 13, 11, 11, 20, 9, 18, 7, 13, 6, 13, 4, 18, 13, 15, 11, 6, 1, 5, 11, 4, 7, 17, 13, 7, 4, 13, 13, 11, 17, 11, 14, 17, 3, 2, 5, 1, 6, 4, 13, 4, 3, 20, 9, 1, 4, 4, 20, 11, 12, 9, 18, 11, 12, 10, 1, 11, 20, 11, 11]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  INSTANCE 4 \n",
    "\n",
    "Elapsed time = 14300.70 sec. (2066464.35 ticks, tree = 14154.38 MB, solutions = 68)\n",
    "Nodefile size = 12084.22 MB (11052.11 MB after compression)\n",
    " 204829 130108     2073.4382   557     1969.0000     2083.0000 28528084    5.79%\n",
    " 205275 130299        cutoff           1969.0000     2083.0000 28571907    5.79%\n",
    " 205806 130426     2078.0006   785     1969.0000     2083.0000 28591965    5.79%\n",
    "*205880+130154                         1970.0000     2083.0000             5.74%\n",
    " 206434 130480     2074.2557   777     1970.0000     2083.0000 28600201    5.74%\n",
    "\n",
    "GUB cover cuts applied:  9\n",
    "Clique cuts applied:  1486\n",
    "Cover cuts applied:  1127\n",
    "Implied bound cuts applied:  3936\n",
    "Flow cuts applied:  208\n",
    "Mixed integer rounding cuts applied:  1445\n",
    "Zero-half cuts applied:  58\n",
    "Lift and project cuts applied:  2\n",
    "Gomory fractional cuts applied:  76\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =    9.56 sec. (4204.05 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14403.58 sec. (2079059.75 ticks)\n",
    "  Sync time (average)   = 3608.73 sec.\n",
    "  Wait time (average)   =    1.59 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14413.14 sec. (2083263.80 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.05736040609136695 2082.9999999999986 time limit exceeded\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1970.0 1970.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1084.4\n",
    "\n",
    " Assignment of products \n",
    " [0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]\n",
    "Prices of products  [26.0, 23.0, 25.0, 28.0, 26.0, 25.0, 24.0, 29.0, 24.0, 23.0, 24.0, 24.0, 24.0, 27.0, 26.0, 20.0, 18.0, 26.0, 22.0, 27.0]\n",
    "Numbered Prod Purchased  [2, 4, 2, 16, 17, 8, 19, 6, 11, 7, 17, 17, 12, 16, 13, 19, 17, 16, 17, 5, 19, 10, 17, 7, 7, 11, 17, 19, 2, 12, 10, 2, 17, 16, 13, 3, 8, 17, 19, 16, 17, 17, 15, 16, 19, 18, 17, 6, 17, 19, 11, 5, 16, 16, 1, 17, 9, 17, 16, 11, 19, 17, 16, 16, 3, 16, 19, 16, 16, 19, 10, 13, 16, 16, 17, 17, 2, 16, 19, 16, 2, 12, 14, 16, 10, 2, 6, 10, 17, 5]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  INSTANCE 5 \n",
    "\n",
    "Elapsed time = 13711.64 sec. (1183557.07 ticks, tree = 25551.63 MB, solutions = 60)\n",
    "Nodefile size = 23464.71 MB (22528.43 MB after compression)\n",
    " 177137 126188     1984.3528   438     1941.0000     2072.0000 13706881    6.75%\n",
    " 177830 126822     2049.8593   540     1941.0000     2072.0000 13812338    6.75%\n",
    " 178605 126901     1976.0438   423     1941.0000     2072.0000 13821232    6.75%\n",
    " 179149 127858     2026.5937   485     1941.0000     2072.0000 13904923    6.75%\n",
    " 179341 127525     2039.4818   475     1941.0000     2072.0000 13879839    6.75%\n",
    " 179548 128170     2062.6005   724     1941.0000     2072.0000 13943286    6.75%\n",
    " 179741 128615     1960.9460   522     1941.0000     2072.0000 14021294    6.75%\n",
    " 180021 128767     1966.4347   537     1941.0000     2072.0000 14057588    6.75%\n",
    " 180515 128946     2050.0234   665     1941.0000     2072.0000 14093337    6.75%\n",
    " 181020 128932     2053.2550   741     1941.0000     2072.0000 14091961    6.75%\n",
    "Elapsed time = 14135.36 sec. (1221826.24 ticks, tree = 26142.20 MB, solutions = 60)\n",
    "Nodefile size = 24047.09 MB (23088.50 MB after compression)\n",
    " 181788 129859     1955.2850   441     1941.0000     2072.0000 14202436    6.75%\n",
    " 182212 129538     2020.6030   492     1941.0000     2072.0000 14165185    6.75%\n",
    " 182467 130297     2071.0000   834     1941.0000     2072.0000 14259166    6.75%\n",
    " 182644 130905     1958.0580   438     1941.0000     2072.0000 14330413    6.75%\n",
    " 182864 131138     1959.8505   414     1941.0000     2072.0000 14377331    6.75%\n",
    " 183137 131207     1993.8379   433     1941.0000     2072.0000 14389344    6.75%\n",
    " 183562 131358     2041.3856   657     1941.0000     2072.0000 14424244    6.75%\n",
    " 184275 131697     1959.8351   537     1941.0000     2072.0000 14461639    6.75%\n",
    "\n",
    "GUB cover cuts applied:  11\n",
    "Clique cuts applied:  1494\n",
    "Cover cuts applied:  1134\n",
    "Implied bound cuts applied:  3468\n",
    "Flow cuts applied:  226\n",
    "Mixed integer rounding cuts applied:  1527\n",
    "Zero-half cuts applied:  36\n",
    "Lift and project cuts applied:  5\n",
    "Gomory fractional cuts applied:  121\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =   10.27 sec. (3598.95 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14407.42 sec. (1252411.04 ticks)\n",
    "  Sync time (average)   = 4644.13 sec.\n",
    "  Wait time (average)   =    2.43 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14417.69 sec. (1256009.99 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.067490984028842 2071.999999999989 time limit exceeded\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1941.0 1941.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1126.6000000000004\n",
    "\n",
    " Assignment of products \n",
    " [0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0]\n",
    "Prices of products  [24.0, 27.0, 18.0, 25.0, 22.0, 20.0, 26.0, 23.0, 28.0, 25.0, 24.0, 23.0, 23.0, 26.0, 21.0, 23.0, 26.0, 27.0, 26.0, 27.0]\n",
    "Numbered Prod Purchased  [3, 15, 3, 13, 3, 4, 12, 9, 3, 10, 15, 14, 3, 11, 15, 3, 3, 3, 5, 5, 6, 5, 3, 10, 3, 20, 3, 17, 3, 8, 3, 3, 6, 3, 8, 14, 2, 20, 13, 1, 14, 11, 7, 3, 5, 15, 3, 8, 3, 8, 3, 12, 19, 6, 1, 11, 6, 1, 15, 5, 20, 15, 3, 3, 3, 15, 2, 5, 6, 5, 3, 10, 1, 19, 6, 4, 19, 6, 3, 6, 8, 15, 0, 3, 13, 18, 5, 6, 12, 16]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "## INSTANCE 6 \n",
    "\n",
    "Elapsed time = 14046.41 sec. (2482851.77 ticks, tree = 1592.20 MB, solutions = 83)\n",
    " 279615 191907     2026.3842   529     1925.0000     2073.0000 45176113    7.69%\n",
    " 279841 192542     2057.8969   558     1925.0000     2073.0000 45279362    7.69%\n",
    " 280034 192742     2042.3405   428     1925.0000     2073.0000 45344348    7.69%\n",
    " 280218 192721     1933.9883   314     1925.0000     2073.0000 45335571    7.69%\n",
    " 280556 192774     2064.4423   789     1925.0000     2073.0000 45353131    7.69%\n",
    " 280985 192994     1957.7239   531     1925.0000     2073.0000 45422918    7.69%\n",
    " 281565 193232     1951.9137   293     1925.0000     2073.0000 45504940    7.69%\n",
    " 282082 193639     1950.4174   365     1925.0000     2073.0000 45610221    7.69%\n",
    " 282560 193569     2031.4364   581     1925.0000     2073.0000 45583224    7.69%\n",
    " 283227 194479     2061.2188   547     1925.0000     2073.0000 45753299    7.69%\n",
    "Elapsed time = 14244.00 sec. (2521097.80 ticks, tree = 1620.52 MB, solutions = 83)\n",
    " 284004 193958     2005.2286   645     1925.0000     2073.0000 45667505    7.69%\n",
    " 284545 194632     2060.0182   779     1925.0000     2073.0000 45773397    7.69%\n",
    " 285009 194788     2060.2778   864     1925.0000     2073.0000 45851476    7.69%\n",
    " 285456 195016     2073.0000   590     1925.0000     2073.0000 45900740    7.69%\n",
    " 285704 196693        cutoff           1925.0000     2073.0000 46212630    7.69%\n",
    " 285922 196512     2012.9173   444     1925.0000     2073.0000 46160209    7.69%\n",
    " 286229 196550     2012.2121   504     1925.0000     2073.0000 46166469    7.69%\n",
    " 286489 196947     1928.6427   572     1925.0000     2073.0000 46282430    7.69%\n",
    "\n",
    "GUB cover cuts applied:  10\n",
    "Clique cuts applied:  1338\n",
    "Cover cuts applied:  1176\n",
    "Implied bound cuts applied:  4640\n",
    "Flow cuts applied:  241\n",
    "Mixed integer rounding cuts applied:  1152\n",
    "Zero-half cuts applied:  53\n",
    "Lift and project cuts applied:  2\n",
    "Gomory fractional cuts applied:  102\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =   10.48 sec. (2929.03 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14389.94 sec. (2552432.29 ticks)\n",
    "  Sync time (average)   = 1269.72 sec.\n",
    "  Wait time (average)   =    0.19 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14400.42 sec. (2555361.32 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.07688311688311265 2072.9999999999995 time limit exceeded\n",
    "\n",
    " Gap is  0.07688311688311265\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1925.0 1925.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1195.9999999999998\n",
    "\n",
    " Assignment of products \n",
    " [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]\n",
    "Prices of products  [27.0, 27.0, 21.0, 19.0, 24.0, 26.0, 24.0, 23.0, 24.0, 24.0, 18.0, 26.0, 29.0, 22.0, 26.0, 24.0, 20.0, 27.0, 23.0, 28.0]\n",
    "Numbered Prod Purchased  [19, 4, 19, 19, 3, 11, 19, 14, 11, 5, 17, 14, 19, 11, 11, 4, 11, 20, 18, 17, 15, 7, 16, 3, 10, 16, 17, 11, 8, 6, 11, 18, 17, 17, 5, 11, 19, 15, 5, 19, 9, 3, 13, 3, 11, 14, 9, 16, 4, 11, 11, 14, 10, 17, 5, 17, 19, 7, 4, 11, 4, 4, 11, 17, 17, 17, 12, 17, 11, 17, 4, 12, 7, 8, 11, 11, 17, 19, 8, 1, 11, 17, 14, 3, 19, 0, 10, 3, 4, 9]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#### INSTANCE 7\n",
    "\n",
    "Elapsed time = 13446.56 sec. (1372457.53 ticks, tree = 26171.01 MB, solutions = 55)\n",
    "Nodefile size = 24095.46 MB (22326.51 MB after compression)\n",
    " 205484 159887     1997.4786   602     1933.0000     2053.8998 16296477    6.25%\n",
    " 205841 161090     2044.7132   729     1933.0000     2053.8975 16425060    6.25%\n",
    " 206534 161137     1965.2788   333     1933.0000     2053.8969 16421897    6.25%\n",
    " 207406 161292     1978.8497   426     1933.0000     2053.8954 16443700    6.25%\n",
    " 208413 162231        cutoff           1933.0000     2053.8945 16539862    6.25%\n",
    " 209297 163113     2052.1551   856     1933.0000     2053.8945 16611076    6.25%\n",
    " 210164 163717     1937.6183   323     1933.0000     2053.8935 16662665    6.25%\n",
    " 210832 164366     2030.9035   536     1933.0000     2053.8935 16713473    6.25%\n",
    " 211159 165340        cutoff           1933.0000     2053.8929 16837390    6.25%\n",
    " 211326 165264     2037.6094   704     1933.0000     2053.8929 16822529    6.25%\n",
    "Elapsed time = 14043.72 sec. (1410704.56 ticks, tree = 27271.43 MB, solutions = 55)\n",
    "Nodefile size = 25197.99 MB (23353.31 MB after compression)\n",
    " 211541 165695     2046.9999   559     1933.0000     2053.8929 16915321    6.25%\n",
    " 211755 165726     2050.9357   649     1933.0000     2053.8929 16919174    6.25%\n",
    " 211948 165771     2049.0524   697     1933.0000     2053.8929 16930243    6.25%\n",
    " 212201 166113     1962.2274   384     1933.0000     2053.8929 16976729    6.25%\n",
    " 212868 166293     1952.8334   327     1933.0000     2053.8929 17003740    6.25%\n",
    " 213727 166441     1954.3046   773     1933.0000     2053.8929 17018762    6.25%\n",
    "\n",
    "GUB cover cuts applied:  2\n",
    "Clique cuts applied:  1016\n",
    "Cover cuts applied:  1084\n",
    "Implied bound cuts applied:  4064\n",
    "Flow cuts applied:  169\n",
    "Mixed integer rounding cuts applied:  1249\n",
    "Zero-half cuts applied:  46\n",
    "Lift and project cuts applied:  2\n",
    "Gomory fractional cuts applied:  116\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =   10.36 sec. (4128.70 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14428.73 sec. (1433339.18 ticks)\n",
    "  Sync time (average)   = 4263.73 sec.\n",
    "  Wait time (average)   =    1.56 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14439.09 sec. (1437467.88 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.06254157120685519 2053.8928571428573 time limit exceeded\n",
    "\n",
    " Gap is  0.06254157120685519\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1933.0 1933.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1179.0000000000014\n",
    "\n",
    " Assignment of products \n",
    " [0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0]\n",
    "Prices of products  [25.0, 26.0, 21.0, 23.0, 26.0, 18.0, 17.0, 21.0, 24.0, 24.0, 27.0, 27.0, 22.0, 23.0, 27.0, 22.0, 20.0, 25.0, 25.0, 25.0]\n",
    "Numbered Prod Purchased  [7, 6, 16, 3, 19, 7, 3, 7, 10, 5, 3, 4, 20, 8, 3, 17, 10, 3, 9, 8, 8, 14, 2, 17, 16, 13, 13, 3, 19, 14, 1, 11, 13, 1, 19, 4, 8, 6, 7, 4, 4, 3, 8, 20, 9, 7, 3, 7, 3, 17, 9, 1, 15, 5, 6, 7, 16, 7, 3, 12, 16, 3, 3, 16, 11, 7, 15, 7, 8, 9, 7, 3, 13, 7, 3, 1, 13, 17, 19, 16, 6, 6, 7, 6, 10, 6, 6, 6, 7, 16]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "######### INSTANCE 8 \n",
    "\n",
    "\n",
    "Elapsed time = 14031.16 sec. (2371713.78 ticks, tree = 12098.53 MB, solutions = 54)\n",
    "Nodefile size = 10030.20 MB (9056.87 MB after compression)\n",
    " 252188 166305     2012.5049   546     1945.0000     2061.0000 36332652    5.96%\n",
    " 252467 166334     2049.6139   722     1945.0000     2061.0000 36336397    5.96%\n",
    " 252800 166923     2000.4163   413     1945.0000     2061.0000 36411961    5.96%\n",
    " 253365 167155     1992.1041   425     1945.0000     2061.0000 36441411    5.96%\n",
    " 253787 167496     1968.2923   432     1945.0000     2061.0000 36506210    5.96%\n",
    " 254558 167189     2016.8915   580     1945.0000     2061.0000 36440461    5.96%\n",
    " 255143 167836     1975.6164   439     1945.0000     2061.0000 36538993    5.96%\n",
    " 255558 168430     1981.1319   344     1945.0000     2061.0000 36618616    5.96%\n",
    " 255953 167867     2061.0000   922     1945.0000     2061.0000 36539907    5.96%\n",
    " 256179 169529     2043.6449   698     1945.0000     2061.0000 36770152    5.96%\n",
    "Elapsed time = 14328.91 sec. (2410000.82 ticks, tree = 12853.51 MB, solutions = 54)\n",
    "Nodefile size = 10791.13 MB (9776.01 MB after compression)\n",
    " 256423 169029     1948.2595   462     1945.0000     2061.0000 36701919    5.96%\n",
    " 256642 169859     1977.8857   568     1945.0000     2061.0000 36818042    5.96%\n",
    " 256839 169793     2051.6079   817     1945.0000     2061.0000 36810798    5.96%\n",
    " 256958 169933     2058.9076  1047     1945.0000     2061.0000 36833775    5.96%\n",
    "\n",
    "GUB cover cuts applied:  12\n",
    "Clique cuts applied:  1419\n",
    "Cover cuts applied:  1290\n",
    "Implied bound cuts applied:  4468\n",
    "Flow cuts applied:  251\n",
    "Mixed integer rounding cuts applied:  1372\n",
    "Zero-half cuts applied:  37\n",
    "Lift and project cuts applied:  2\n",
    "Gomory fractional cuts applied:  77\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =    7.39 sec. (3461.08 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14398.89 sec. (2424365.83 ticks)\n",
    "  Sync time (average)   = 1934.07 sec.\n",
    "  Wait time (average)   =    0.83 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14406.28 sec. (2427826.91 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.05964010282775903 2060.9999999999973 time limit exceeded\n",
    "\n",
    " Gap is  0.05964010282775903\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1945.0 1945.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1175.1999999999996\n",
    "\n",
    " Assignment of products \n",
    " [0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0]\n",
    "Prices of products  [25.0, 23.0, 26.0, 24.0, 23.0, 17.0, 24.0, 24.0, 29.0, 24.0, 19.0, 25.0, 26.0, 23.0, 26.0, 24.0, 27.0, 22.0, 20.0, 26.0]\n",
    "Numbered Prod Purchased  [7, 18, 16, 19, 7, 19, 14, 11, 14, 14, 11, 6, 20, 19, 5, 4, 11, 13, 8, 8, 19, 6, 2, 17, 8, 13, 1, 5, 19, 14, 1, 11, 19, 11, 19, 14, 11, 8, 7, 10, 4, 11, 4, 2, 8, 7, 9, 0, 14, 11, 9, 11, 6, 6, 6, 19, 12, 18, 6, 12, 16, 11, 2, 15, 6, 19, 2, 15, 11, 5, 7, 3, 4, 18, 2, 6, 18, 6, 19, 6, 4, 16, 19, 6, 10, 6, 11, 11, 7, 11]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  INSTANCE 9\n",
    "\n",
    " \n",
    "Elapsed time = 13990.64 sec. (1950093.57 ticks, tree = 29367.85 MB, solutions = 72)\n",
    "Nodefile size = 27252.34 MB (26092.54 MB after compression)\n",
    " 201515 148988     2039.1680   825     1969.0000     2070.0000 18737165    5.13%\n",
    " 201796 149634     1983.3146   623     1969.0000     2070.0000 18785546    5.13%\n",
    " 202024 150460     2049.6348   671     1969.0000     2070.0000 18906544    5.13%\n",
    " 202221 150772     2031.9252   544     1969.0000     2070.0000 18960071    5.13%\n",
    " 202423 150692     1987.1851   466     1969.0000     2070.0000 18925715    5.13%\n",
    " 202581 150733     2012.3643   787     1969.0000     2070.0000 18956307    5.13%\n",
    " 202898 151068     1993.4021   408     1969.0000     2070.0000 19012737    5.13%\n",
    " 203438 151222     2028.1705   543     1969.0000     2070.0000 19025666    5.13%\n",
    " 203981 151294     2040.8000   669     1969.0000     2070.0000 19052201    5.13%\n",
    "*204386+151676                         1971.0000     2070.0000             5.02%\n",
    " 204472 151959     2032.9642   561     1971.0000     2070.0000 19129994    5.02%\n",
    "Elapsed time = 14247.00 sec. (1988380.38 ticks, tree = 30072.90 MB, solutions = 73)\n",
    "Nodefile size = 27988.24 MB (26793.15 MB after compression)\n",
    " 205004 151650     2011.9799   545     1971.0000     2070.0000 19211753    5.02%\n",
    " 205424 151798     1997.9346   603     1971.0000     2070.0000 19236083    5.02%\n",
    " 205884 151823     2026.2492   633     1971.0000     2070.0000 19239750    5.02%\n",
    " 206168 152477     2016.4992   566     1971.0000     2070.0000 19318089    5.02%\n",
    " 206285 152546     1998.0268   637     1971.0000     2070.0000 19325610    5.02%\n",
    "\n",
    "GUB cover cuts applied:  7\n",
    "Clique cuts applied:  1635\n",
    "Cover cuts applied:  1076\n",
    "Implied bound cuts applied:  3821\n",
    "Flow cuts applied:  274\n",
    "Mixed integer rounding cuts applied:  1726\n",
    "Zero-half cuts applied:  40\n",
    "Lift and project cuts applied:  4\n",
    "Gomory fractional cuts applied:  102\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =    8.03 sec. (3665.57 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14403.00 sec. (2007680.31 ticks)\n",
    "  Sync time (average)   = 4319.69 sec.\n",
    "  Wait time (average)   =    1.91 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14411.03 sec. (2011345.88 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.05022831050227986 2069.9999999999986 time limit exceeded\n",
    "\n",
    " Gap is  0.05022831050227986\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1971.0 1971.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1163.4\n",
    "\n",
    " Assignment of products \n",
    " [0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0]\n",
    "Prices of products  [24.0, 25.0, 28.0, 23.0, 21.0, 27.0, 22.0, 24.0, 23.0, 19.0, 29.0, 26.0, 17.0, 26.0, 27.0, 23.0, 27.0, 23.0, 21.0, 22.0]\n",
    "Numbered Prod Purchased  [10, 13, 19, 19, 1, 10, 13, 20, 9, 19, 16, 13, 12, 5, 3, 18, 6, 10, 12, 19, 8, 4, 10, 19, 4, 10, 11, 5, 10, 7, 14, 13, 4, 10, 15, 19, 18, 16, 20, 5, 20, 1, 13, 3, 13, 12, 15, 13, 2, 19, 14, 5, 5, 19, 2, 19, 10, 10, 2, 7, 18, 15, 5, 20, 18, 13, 20, 7, 10, 5, 17, 10, 16, 5, 10, 20, 5, 10, 19, 8, 7, 8, 17, 18, 13, 13, 8, 2, 19, 13]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#  INSTANCE 10\n",
    "\n",
    "Elapsed time = 14227.05 sec. (2415121.25 ticks, tree = 11619.77 MB, solutions = 52)\n",
    "Nodefile size = 9555.19 MB (8452.62 MB after compression)\n",
    "*246387+166846                         1949.0000     2089.9975             7.23%\n",
    " 246457 167074     2068.0646   611     1949.0000     2089.9975 38895379    7.23%\n",
    " 246730 166921     2013.9832   332     1949.0000     2089.9975 38843475    7.23%\n",
    " 246963 167035     2062.7656   604     1949.0000     2089.9975 38874664    7.23%\n",
    " 247572 164520     2043.0005   475     1949.0000     2089.9975 38969470    7.23%\n",
    " 248004 164229     1996.4378   611     1949.0000     2089.9975 38922161    7.23%\n",
    " 248211 164347        cutoff           1949.0000     2089.9975 38937100    7.23%\n",
    "\n",
    "GUB cover cuts applied:  14\n",
    "Clique cuts applied:  1402\n",
    "Cover cuts applied:  1407\n",
    "Implied bound cuts applied:  4093\n",
    "Flow cuts applied:  248\n",
    "Mixed integer rounding cuts applied:  1412\n",
    "Zero-half cuts applied:  57\n",
    "Lift and project cuts applied:  3\n",
    "Gomory fractional cuts applied:  109\n",
    "\n",
    "Root node processing (before b&c):\n",
    "  Real time             =    5.25 sec. (2825.29 ticks)\n",
    "Parallel b&c, 8 threads:\n",
    "  Real time             = 14405.89 sec. (2437825.52 ticks)\n",
    "  Sync time (average)   = 1784.05 sec.\n",
    "  Wait time (average)   =    0.43 sec.\n",
    "                          ------------\n",
    "Total (root+branch&cut) = 14411.14 sec. (2440650.81 ticks)\n",
    "Solve Details in order mip relative gap, best bound 0.07234322981726599 2089.9969549138586 time limit exceeded\n",
    "\n",
    " Gap is  0.07234322981726599\n",
    "time limit =  14400.0\n",
    "threads =  0\n",
    "threads   0\n",
    "\n",
    " Upper Level Objective \n",
    " 1949.0 1949.0\n",
    "\n",
    " Lower Level Objective \n",
    " 1106.3999999999996\n",
    "\n",
    " Assignment of products \n",
    " [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.0, 0, 0, 0, 0]\n",
    "Prices of products  [24.0, 25.0, 28.0, 28.0, 28.0, 28.0, 26.0, 24.0, 21.0, 22.0, 18.0, 25.0, 27.0, 26.0, 20.0, 23.0, 28.0, 29.0, 27.0, 28.0]\n",
    "Numbered Prod Purchased  [11, 11, 11, 9, 15, 10, 7, 7, 8, 8, 11, 10, 2, 9, 15, 9, 10, 2, 15, 10, 13, 12, 15, 7, 9, 8, 8, 15, 14, 15, 11, 5, 3, 0, 11, 10, 13, 11, 15, 11, 10, 16, 2, 0, 9, 10, 6, 7, 16, 19, 15, 16, 12, 15, 10, 11, 15, 15, 14, 11, 11, 9, 11, 15, 11, 1, 11, 20, 14, 18, 11, 7, 19, 1, 16, 10, 10, 8, 10, 4, 10, 9, 14, 10, 15, 1, 10, 9, 0, 16]"
   ]
  }
 ],
 "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.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}