Graph-SCP.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "2dc7b190",
   "metadata": {},
   "source": [
    "### Imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "a5a72020",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/zohairshafi/miniforge3/lib/python3.10/site-packages/tqdm/auto.py:22: TqdmWarning: IProgress not found. Please update jupyter and ipywidgets. See https://ipywidgets.readthedocs.io/en/stable/user_install.html\n",
      "  from .autonotebook import tqdm as notebook_tqdm\n"
     ]
    }
   ],
   "source": [
    "%run utils.py"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "dc649651",
   "metadata": {},
   "source": [
    "### Data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "13b92f38",
   "metadata": {},
   "outputs": [],
   "source": [
    "with open('data/training_instances_all.pkl', 'rb') as file: \n",
    "    train_inst = pkl.load(file)\n",
    "\n",
    "with open('data/test_instance_type_1.pkl', 'rb') as file: \n",
    "    test_1 = pkl.load(file)\n",
    "    \n",
    "with open('data/test_equal_instances.pkl', 'rb') as file: \n",
    "    test_2 = pkl.load(file)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5ef3a092",
   "metadata": {},
   "source": [
    "### Train Models"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "75511e24",
   "metadata": {},
   "outputs": [],
   "source": [
    "f1_list, auc_list, aup_list, auc_list, ratio_list = train_set_cover_unsup(train_instances = train_inst,\n",
    "                                                        epochs = 200,\n",
    "                                                        patience = 100,\n",
    "                                                        cost_type = 'equal',\n",
    "                                                        model_path = None,\n",
    "                                                        model_name = './models/set_cover_sage', \n",
    "                                                        m_type = 'sage',\n",
    "                                                        feats = ['costs', 'cover', 'ppr', 'bi_sf', 'hyp_feat'])\n",
    "\n",
    "f1_list, auc_list, aup_list, auc_list, ratio_list = train_set_cover_unsup(train_instances = train_inst,\n",
    "                                                        epochs = 200,\n",
    "                                                        patience = 100,\n",
    "                                                        cost_type = 'custom',\n",
    "                                                        model_path = './models/set_cover_sage.pkl',\n",
    "                                                        model_name = './models/set_cover_sage', \n",
    "                                                        m_type = 'sage',\n",
    "                                                        feats = ['costs', 'cover', 'ppr', 'bi_sf', 'hyp_feat'])\n",
    "f1_list, auc_list, aup_list, auc_list, ratio_list = train_set_cover_unsup(train_instances = train_inst,\n",
    "                                                        epochs = 200,\n",
    "                                                        patience = 100,\n",
    "                                                        cost_type = 'custom',\n",
    "                                                        model_path = './models/set_cover_sage.pkl',\n",
    "                                                        model_name = './models/set_cover_sage', \n",
    "                                                        m_type = 'sage',\n",
    "                                                        feats = ['costs', 'cover', 'ppr', 'bi_sf', 'hyp_feat'])\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7986263a",
   "metadata": {},
   "source": [
    "### Custom SCP Instances"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "0845af58",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 14 rows, 5 columns and 28 nonzeros\n",
      "Model fingerprint: 0x6777a76e\n",
      "Variable types: 0 continuous, 5 integer (5 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+00, 1e+00]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Found heuristic solution: objective 2.0000000\n",
      "Variable types: 0 continuous, 5 integer (5 binary)\n",
      "\n",
      "Root relaxation: cutoff, 3 iterations, 0.00 seconds (0.00 work units)\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "     0     0     cutoff    0         2.00000    2.00000  0.00%     -    0s\n",
      "\n",
      "Explored 1 nodes (3 simplex iterations) in 0.01 seconds (0.00 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 1: 2 \n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 2.000000000000e+00, best bound 2.000000000000e+00, gap 0.0000%\n",
      "\n",
      "User-callback calls 39, time in user-callback 0.00 sec\n",
      "20\n",
      "Set parameter Presolve to value 0                         \n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 14 rows, 0 columns and 0 nonzeros\n",
      "Model fingerprint: 0x47b2dfcb\n",
      "Coefficient statistics:\n",
      "  Matrix range     [0e+00, 0e+00]\n",
      "  Objective range  [0e+00, 0e+00]\n",
      "  Bounds range     [0e+00, 0e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Iteration    Objective       Primal Inf.    Dual Inf.      Time\n",
      "       0    0.0000000e+00   1.400000e+01   0.000000e+00      0s\n",
      "\n",
      "Solved in 0 iterations and 0.00 seconds (0.00 work units)\n",
      "Infeasible model\n",
      "\n",
      "User-callback calls 19, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 14 rows, 0 columns and 0 nonzeros\n",
      "Model fingerprint: 0x47b2dfcb\n",
      "Coefficient statistics:\n",
      "  Matrix range     [0e+00, 0e+00]\n",
      "  Objective range  [0e+00, 0e+00]\n",
      "  Bounds range     [0e+00, 0e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Iteration    Objective       Primal Inf.    Dual Inf.      Time\n",
      "       0    0.0000000e+00   1.400000e+01   0.000000e+00      0s\n",
      "\n",
      "Solved in 0 iterations and 0.00 seconds (0.00 work units)\n",
      "Infeasible model\n",
      "\n",
      "User-callback calls 19, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 14 rows, 0 columns and 0 nonzeros\n",
      "Model fingerprint: 0x47b2dfcb\n",
      "Coefficient statistics:\n",
      "  Matrix range     [0e+00, 0e+00]\n",
      "  Objective range  [0e+00, 0e+00]\n",
      "  Bounds range     [0e+00, 0e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Iteration    Objective       Primal Inf.    Dual Inf.      Time\n",
      "       0    0.0000000e+00   1.400000e+01   0.000000e+00      0s\n",
      "\n",
      "Solved in 0 iterations and 0.00 seconds (0.00 work units)\n",
      "Infeasible model\n",
      "\n",
      "User-callback calls 19, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 14 rows, 1 columns and 8 nonzeros\n",
      "Model fingerprint: 0x73afaa05\n",
      "Variable types: 0 continuous, 1 integer (1 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+00, 1e+00]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R0 by 1.000000000\n",
      "\n",
      "\n",
      "Explored 0 nodes (0 simplex iterations) in 0.00 seconds (0.00 work units)\n",
      "Thread count was 1 (of 10 available processors)\n",
      "\n",
      "Solution count 0\n",
      "\n",
      "Model is infeasible\n",
      "Best objective -, best bound -, gap -\n",
      "\n",
      "User-callback calls 27, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 14 rows, 3 columns and 22 nonzeros\n",
      "Model fingerprint: 0x64909387\n",
      "Variable types: 0 continuous, 3 integer (3 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+00, 1e+00]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R0 by 1.000000000\n",
      "\n",
      "Found heuristic solution: objective 2.0000000\n",
      "Variable types: 0 continuous, 3 integer (3 binary)\n",
      "\n",
      "Root relaxation: cutoff, 0 iterations, 0.00 seconds (0.00 work units)\n",
      "\n",
      "Explored 1 nodes (0 simplex iterations) in 0.00 seconds (0.00 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 1: 2 \n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 2.000000000000e+00, best bound 2.000000000000e+00, gap 0.0000%\n",
      "\n",
      "User-callback calls 38, time in user-callback 0.00 sec\n",
      "#########################################\n",
      "Graph-SCP Time:  0.03359508514404297\n",
      "Random Time: 0 0 0\n",
      "Gurobi Time: 0.0124\n",
      "Greedy Time: 1.3113021850585938e-05\n",
      "\n",
      "Optimal Objective:  2\n",
      "Gurobi Objective:  2\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  1.0\n",
      "Gurobi Node Count :  1.0\n",
      "\n",
      "Num Times Gurobi Run :  5\n",
      "Problem Size :  3 / 5\n",
      "#########################################\n",
      "################################\n",
      "#############0/1################\n",
      "################################\n"
     ]
    }
   ],
   "source": [
    "# Simple test instance \n",
    "\n",
    "inst = {}\n",
    "inst['universe'] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14}\n",
    "inst['subsets'] = {\n",
    "                   0 : {1, 2}, \n",
    "                   1 : {3, 4, 5, 6}, \n",
    "                   2 : {7, 8, 9, 10, 11, 12, 13, 14}, \n",
    "                   4 : {1, 3, 4, 7, 8, 9, 10}, \n",
    "                   5 : {2, 5, 6, 11, 12, 13, 14},\n",
    "                  }\n",
    "inst['solutions'] = [['4', '5']]\n",
    "\n",
    "results, inst = test_set_cover(test_instances = [inst],\n",
    "                                 num_runs = 10,\n",
    "                                 cost_type = 'equal',\n",
    "                                 model_path = './models/set_cover_sage_unsup.pkl',\n",
    "                                 m_type = 'sage',\n",
    "                                 feats = ['costs', 'cover', 'ppr', 'bi_sf', 'hyp_feat'],\n",
    "                                 init_sol_perc = 55,\n",
    "                                 perc_start = 55,\n",
    "                                 dec_size = 10, \n",
    "                                 obj_threshold = 1, \n",
    "                                 time_limit = 3600, \n",
    "                                 run_random = False, \n",
    "                                 gnn_solver = 'gurobi', \n",
    "                                 neurons = 512)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cdbcd1d2",
   "metadata": {},
   "source": [
    "### Larger SCP instances using SCIP"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "f0cf5907",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "651\n",
      "#########################################                 \n",
      "Graph-SCP Time:  2.3824403285980225\n",
      "Random Time: 0 0 0\n",
      "SCIP Time: 414.2942190170288\n",
      "Greedy Time: 0.015257120132446289\n",
      "\n",
      "Optimal Objective:  8\n",
      "SCIP Objective:  8\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  158603\n",
      "SCIP Node Count :  2646\n",
      "\n",
      "Num Times Gurobi Run :  1\n",
      "Problem Size :  65 / 350\n",
      "#########################################\n",
      "################################\n",
      "#############0/5################\n",
      "################################\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "651\n",
      "#########################################                 \n",
      "Graph-SCP Time:  2.6388721466064453\n",
      "Random Time: 0 0 0\n",
      "SCIP Time: 412.30747389793396\n",
      "Greedy Time: 0.016705989837646484\n",
      "\n",
      "Optimal Objective:  8\n",
      "SCIP Objective:  8\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  161099\n",
      "SCIP Node Count :  3697\n",
      "\n",
      "Num Times Gurobi Run :  1\n",
      "Problem Size :  65 / 350\n",
      "#########################################\n",
      "################################\n",
      "#############1/5################\n",
      "################################\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "651\n",
      "#########################################                 \n",
      "Graph-SCP Time:  3.1719789505004883\n",
      "Random Time: 0 0 0\n",
      "SCIP Time: 415.5094940662384\n",
      "Greedy Time: 0.015203714370727539\n",
      "\n",
      "Optimal Objective:  8\n",
      "SCIP Objective:  8\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  174071\n",
      "SCIP Node Count :  4368\n",
      "\n",
      "Num Times Gurobi Run :  1\n",
      "Problem Size :  65 / 350\n",
      "#########################################\n",
      "################################\n",
      "#############2/5################\n",
      "################################\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "651\n",
      "#########################################                 \n",
      "Graph-SCP Time:  4.110252141952515\n",
      "Random Time: 0 0 0\n",
      "SCIP Time: 250.72438406944275\n",
      "Greedy Time: 0.0168759822845459\n",
      "\n",
      "Optimal Objective:  8\n",
      "SCIP Objective:  8\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  108853\n",
      "SCIP Node Count :  6243\n",
      "\n",
      "Num Times Gurobi Run :  1\n",
      "Problem Size :  65 / 350\n",
      "#########################################\n",
      "################################\n",
      "#############3/5################\n",
      "################################\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "651\n",
      "#########################################                 \n",
      "Graph-SCP Time:  3.5928940773010254\n",
      "Random Time: 0 0 0\n",
      "SCIP Time: 476.8254590034485\n",
      "Greedy Time: 0.015661001205444336\n",
      "\n",
      "Optimal Objective:  8\n",
      "SCIP Objective:  8\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  197195\n",
      "SCIP Node Count :  5345\n",
      "\n",
      "Num Times Gurobi Run :  1\n",
      "Problem Size :  65 / 350\n",
      "#########################################\n",
      "################################\n",
      "#############4/5################\n",
      "################################\n"
     ]
    }
   ],
   "source": [
    "results, inst = test_set_cover(test_instances = test_2[-5:],\n",
    "                                 num_runs = 10,\n",
    "                                 cost_type = 'equal',\n",
    "                                 model_path = './models/set_cover_sage_unsup.pkl',\n",
    "                                 m_type = 'sage',\n",
    "                                 feats = ['costs', 'cover', 'ppr', 'bi_sf', 'hyp_feat'],\n",
    "                                 init_sol_perc = 90,\n",
    "                                 perc_start = 90,\n",
    "                                 dec_size = 10, \n",
    "                                 obj_threshold = 1, \n",
    "                                 time_limit = 3600, \n",
    "                                 run_random = False, \n",
    "                                 gnn_solver = 'scip',\n",
    "                                 neurons = 512)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "26afcef0",
   "metadata": {},
   "source": [
    "### Larger SCP instances using Gurobi "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "3dd5fc30",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 250 columns and 9693 nonzeros\n",
      "Model fingerprint: 0x7e35746b\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Found heuristic solution: objective 2054.0000000\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Found heuristic solution: objective 1520.0000000\n",
      "\n",
      "Root relaxation: objective 4.629387e+02, 214 iterations, 0.01 seconds (0.01 work units)\n",
      "pressed CTRL-C 1 times (5 times for forcing termination)\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "     0     0  462.93873    0   48 1520.00000  462.93873  69.5%     -    0s\n",
      "H    0     0                    1387.0000000  462.93873  66.6%     -    0s\n",
      "H    0     0                    1226.0000000  462.93873  62.2%     -    0s\n",
      "H    0     0                     975.0000000  462.93873  52.5%     -    0s\n",
      "     0     0  480.54333    0   57  975.00000  480.54333  50.7%     -    0s\n",
      "H    0     0                     884.0000000  480.54333  45.6%     -    0s\n",
      "     0     0  503.26150    0   57  884.00000  503.26150  43.1%     -    0s\n",
      "     0     0  503.26150    0   58  884.00000  503.26150  43.1%     -    0s\n",
      "     0     0  503.52028    0   53  884.00000  503.52028  43.0%     -    0s\n",
      "     0     0  505.37418    0   52  884.00000  505.37418  42.8%     -    0s\n",
      "     0     0  505.37418    0   52  884.00000  505.37418  42.8%     -    0s\n",
      "     0     0  505.37418    0   51  884.00000  505.37418  42.8%     -    0s\n",
      "     0     0  505.37418    0   51  884.00000  505.37418  42.8%     -    0s\n",
      "     0     0  505.63917    0   53  884.00000  505.63917  42.8%     -    0s\n",
      "     0     0  505.63917    0   53  884.00000  505.63917  42.8%     -    0s\n",
      "     0     0  506.58394    0   53  884.00000  506.58394  42.7%     -    0s\n",
      "     0     0  507.41442    0   52  884.00000  507.41442  42.6%     -    0s\n",
      "     0     0  507.41442    0   52  884.00000  507.41442  42.6%     -    0s\n",
      "     0     0  508.94672    0   51  884.00000  508.94672  42.4%     -    0s\n",
      "     0     0  508.94672    0   51  884.00000  508.94672  42.4%     -    0s\n",
      "     0     0  508.94672    0   51  884.00000  508.94672  42.4%     -    0s\n",
      "     0     0  508.94672    0   50  884.00000  508.94672  42.4%     -    0s\n",
      "     0     0  509.38589    0   50  884.00000  509.38589  42.4%     -    0s\n",
      "     0     0  531.25129    0   55  884.00000  531.25129  39.9%     -    0s\n",
      "     0     2  531.25129    0   55  884.00000  531.25129  39.9%     -    0s\n",
      "H  359   370                     867.0000000  531.25129  38.7%  28.8    0s\n",
      "H  568   515                     839.0000000  531.25129  36.7%  27.3    0s\n",
      "H  572   515                     827.0000000  531.25129  35.8%  27.2    0s\n",
      "H  576   515                     818.0000000  531.25129  35.1%  27.4    0s\n",
      "H 1737  1256                     806.0000000  570.88181  29.2%  26.0    0s\n",
      "* 2449  1471              26     783.0000000  570.88181  27.1%  29.1    1s\n",
      "* 4780  1700              25     771.0000000  570.88181  26.0%  28.6    1s\n",
      "*10135  1615              23     753.0000000  699.50198  7.10%  27.1    2s\n",
      "\n",
      "Cutting planes:\n",
      "  Gomory: 1\n",
      "  MIR: 15\n",
      "  Zero half: 4\n",
      "\n",
      "Explored 12074 nodes (310637 simplex iterations) in 2.25 seconds (4.83 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 10: 753 771 783 ... 975\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 7.530000000000e+02, best bound 7.530000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 25003, time in user-callback 0.03 sec\n",
      "502\n",
      "Set parameter Presolve to value 0                         \n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 3 columns and 241 nonzeros\n",
      "Model fingerprint: 0x330d88ff\n",
      "Variable types: 0 continuous, 3 integer (3 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 1e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R2 by 1.000000000\n",
      "\n",
      "\n",
      "Explored 0 nodes (0 simplex iterations) in 0.00 seconds (0.00 work units)\n",
      "Thread count was 1 (of 10 available processors)\n",
      "\n",
      "Solution count 0\n",
      "\n",
      "Model is infeasible\n",
      "Best objective -, best bound -, gap -\n",
      "\n",
      "User-callback calls 27, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 6 columns and 475 nonzeros\n",
      "Model fingerprint: 0x78d47e83\n",
      "Variable types: 0 continuous, 6 integer (6 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 1e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R2 by 1.000000000\n",
      "\n",
      "\n",
      "Explored 0 nodes (0 simplex iterations) in 0.00 seconds (0.00 work units)\n",
      "Thread count was 1 (of 10 available processors)\n",
      "\n",
      "Solution count 0\n",
      "\n",
      "Model is infeasible\n",
      "Best objective -, best bound -, gap -\n",
      "\n",
      "User-callback calls 27, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 75 columns and 4654 nonzeros\n",
      "Model fingerprint: 0x6b89bd63\n",
      "Variable types: 0 continuous, 75 integer (75 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R2 by 1.000000000\n",
      "\n",
      "Found heuristic solution: objective 1624.0000000\n",
      "Variable types: 0 continuous, 75 integer (75 binary)\n",
      "Found heuristic solution: objective 1276.0000000\n",
      "\n",
      "Root relaxation: objective 4.678607e+02, 148 iterations, 0.00 seconds (0.00 work units)\n",
      "Another try with MIP start\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "H    0     0                     923.0000000  467.86065  49.3%     -    0s\n",
      "     0     0  467.86065    0   46  923.00000  467.86065  49.3%     -    0s\n",
      "H    0     0                     820.0000000  467.86065  42.9%     -    0s\n",
      "     0     0  487.93180    0   50  820.00000  487.93180  40.5%     -    0s\n",
      "     0     0  518.76371    0   51  820.00000  518.76371  36.7%     -    0s\n",
      "     0     0  518.77334    0   50  820.00000  518.77334  36.7%     -    0s\n",
      "     0     0  521.13677    0   48  820.00000  521.13677  36.4%     -    0s\n",
      "     0     0  524.32844    0   48  820.00000  524.32844  36.1%     -    0s\n",
      "     0     0  528.35740    0   49  820.00000  528.35740  35.6%     -    0s\n",
      "     0     0  533.96105    0   50  820.00000  533.96105  34.9%     -    0s\n",
      "     0     0  537.22214    0   49  820.00000  537.22214  34.5%     -    0s\n",
      "     0     0  537.62873    0   49  820.00000  537.62873  34.4%     -    0s\n",
      "     0     0  561.63895    0   51  820.00000  561.63895  31.5%     -    0s\n",
      "     0     0  562.58887    0   50  820.00000  562.58887  31.4%     -    0s\n",
      "     0     0  563.57381    0   51  820.00000  563.57381  31.3%     -    0s\n",
      "     0     0  565.13554    0   46  820.00000  565.13554  31.1%     -    0s\n",
      "     0     0  566.23786    0   47  820.00000  566.23786  30.9%     -    0s\n",
      "     0     0  569.53739    0   48  820.00000  569.53739  30.5%     -    0s\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "     0     0  569.59179    0   47  820.00000  569.59179  30.5%     -    0s\n",
      "     0     0  570.33685    0   48  820.00000  570.33685  30.4%     -    0s\n",
      "     0     0  570.64369    0   49  820.00000  570.64369  30.4%     -    0s\n",
      "     0     0  572.59887    0   49  820.00000  572.59887  30.2%     -    0s\n",
      "     0     0  572.69718    0   48  820.00000  572.69718  30.2%     -    0s\n",
      "     0     0  572.69718    0   49  820.00000  572.69718  30.2%     -    0s\n",
      "     0     0  590.56441    0   44  820.00000  590.56441  28.0%     -    0s\n",
      "     0     0  590.56441    0   44  820.00000  590.56441  28.0%     -    0s\n",
      "     0     2  590.56441    0   44  820.00000  590.56441  28.0%     -    0s\n",
      "*  123   121              13     815.0000000  590.56441  27.5%  29.2    0s\n",
      "H  370   241                     783.0000000  590.56441  24.6%  23.2    0s\n",
      "* 3612  1073              25     771.0000000  679.10434  11.9%  19.2    0s\n",
      "* 5272   364               7     753.0000000  726.98562  3.45%  19.0    0s\n",
      "\n",
      "Cutting planes:\n",
      "  Gomory: 3\n",
      "  MIR: 48\n",
      "  Zero half: 4\n",
      "\n",
      "Explored 5740 nodes (105706 simplex iterations) in 0.61 seconds (0.95 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 8: 753 771 783 ... 1624\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 7.530000000000e+02, best bound 7.530000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 11814, time in user-callback 0.02 sec\n",
      "#########################################\n",
      "Graph-SCP Time:  0.6466751098632812\n",
      "Random Time: 0 0 0\n",
      "Gurobi Time: 2.262\n",
      "Greedy Time: 0.0016622543334960938\n",
      "\n",
      "Optimal Objective:  753\n",
      "Gurobi Objective:  753\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  12074.0\n",
      "Gurobi Node Count :  5740.0\n",
      "\n",
      "Num Times Gurobi Run :  3\n",
      "Problem Size :  75 / 250\n",
      "#########################################\n",
      "################################\n",
      "#############0/5################\n",
      "################################\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 250 columns and 10469 nonzeros\n",
      "Model fingerprint: 0xed2a6546\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Found heuristic solution: objective 1533.0000000\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Found heuristic solution: objective 1443.0000000\n",
      "\n",
      "Root relaxation: objective 4.445766e+02, 190 iterations, 0.00 seconds (0.01 work units)\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "     0     0  444.57665    0   55 1443.00000  444.57665  69.2%     -    0s\n",
      "H    0     0                    1380.0000000  444.57665  67.8%     -    0s\n",
      "H    0     0                    1256.0000000  444.57665  64.6%     -    0s\n",
      "H    0     0                     896.0000000  444.57665  50.4%     -    0s\n",
      "     0     0  455.38658    0   60  896.00000  455.38658  49.2%     -    0s\n",
      "H    0     0                     848.0000000  455.38658  46.3%     -    0s\n",
      "H    0     0                     799.0000000  506.56572  36.6%     -    0s\n",
      "     0     0  506.56572    0   61  799.00000  506.56572  36.6%     -    0s\n",
      "     0     0  506.56572    0   57  799.00000  506.56572  36.6%     -    0s\n",
      "     0     0  506.56572    0   56  799.00000  506.56572  36.6%     -    0s\n",
      "     0     0  506.56572    0   57  799.00000  506.56572  36.6%     -    0s\n",
      "     0     0  506.56572    0   56  799.00000  506.56572  36.6%     -    0s\n",
      "     0     0  506.56572    0   57  799.00000  506.56572  36.6%     -    0s\n",
      "     0     0  506.56572    0   58  799.00000  506.56572  36.6%     -    0s\n",
      "     0     0  506.56572    0   58  799.00000  506.56572  36.6%     -    0s\n",
      "     0     0  529.73340    0   64  799.00000  529.73340  33.7%     -    0s\n",
      "     0     2  529.73340    0   64  799.00000  529.73340  33.7%     -    0s\n",
      "H   88    89                     795.0000000  529.73340  33.4%  41.9    0s\n",
      "* 1148   841              19     793.0000000  529.73340  33.2%  32.0    0s\n",
      "H 1625  1120                     783.0000000  532.23879  32.0%  31.9    1s\n",
      "H 1658  1089                     779.0000000  532.23879  31.7%  32.3    1s\n",
      "H 1975  1183                     771.0000000  532.23879  31.0%  33.5    1s\n",
      "*15515  1274              32     770.0000000  723.98010  5.98%  28.0    3s\n",
      "\n",
      "Cutting planes:\n",
      "  Gomory: 1\n",
      "  MIR: 23\n",
      "\n",
      "Explored 17505 nodes (474458 simplex iterations) in 3.13 seconds (7.17 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 10: 770 771 779 ... 1256\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 7.700000000000e+02, best bound 7.700000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 36256, time in user-callback 0.04 sec\n",
      "502\n",
      "Set parameter Presolve to value 0                         \n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 2 columns and 160 nonzeros\n",
      "Model fingerprint: 0xafedb600\n",
      "Variable types: 0 continuous, 2 integer (2 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 1e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R0 by 1.000000000\n",
      "\n",
      "\n",
      "Explored 0 nodes (0 simplex iterations) in 0.01 seconds (0.00 work units)\n",
      "Thread count was 1 (of 10 available processors)\n",
      "\n",
      "Solution count 0\n",
      "\n",
      "Model is infeasible\n",
      "Best objective -, best bound -, gap -\n",
      "\n",
      "User-callback calls 27, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 10 columns and 742 nonzeros\n",
      "Model fingerprint: 0xb7a4275f\n",
      "Variable types: 0 continuous, 10 integer (10 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 1e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R0 by 1.000000000\n",
      "\n",
      "\n",
      "Explored 0 nodes (0 simplex iterations) in 0.01 seconds (0.00 work units)\n",
      "Thread count was 1 (of 10 available processors)\n",
      "\n",
      "Solution count 0\n",
      "\n",
      "Model is infeasible\n",
      "Best objective -, best bound -, gap -\n",
      "\n",
      "User-callback calls 27, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 76 columns and 4922 nonzeros\n",
      "Model fingerprint: 0x1e17695b\n",
      "Variable types: 0 continuous, 76 integer (76 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R0 by 1.000000000\n",
      "\n",
      "Found heuristic solution: objective 1350.0000000\n",
      "Variable types: 0 continuous, 76 integer (76 binary)\n",
      "Found heuristic solution: objective 1320.0000000\n",
      "\n",
      "Root relaxation: objective 4.513184e+02, 138 iterations, 0.00 seconds (0.00 work units)\n",
      "Another try with MIP start\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "H    0     0                     882.0000000  451.31840  48.8%     -    0s\n",
      "     0     0  451.31840    0   43  882.00000  451.31840  48.8%     -    0s\n",
      "H    0     0                     857.0000000  451.31840  47.3%     -    0s\n",
      "H    0     0                     784.0000000  451.31840  42.4%     -    0s\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "     0     0  470.45963    0   47  784.00000  470.45963  40.0%     -    0s\n",
      "     0     0  491.62071    0   45  784.00000  491.62071  37.3%     -    0s\n",
      "     0     0  493.06918    0   47  784.00000  493.06918  37.1%     -    0s\n",
      "     0     0  495.40723    0   43  784.00000  495.40723  36.8%     -    0s\n",
      "     0     0  495.40723    0   44  784.00000  495.40723  36.8%     -    0s\n",
      "     0     0  501.26787    0   43  784.00000  501.26787  36.1%     -    0s\n",
      "     0     0  501.26787    0   43  784.00000  501.26787  36.1%     -    0s\n",
      "     0     0  501.26787    0   45  784.00000  501.26787  36.1%     -    0s\n",
      "     0     0  523.94716    0   46  784.00000  523.94716  33.2%     -    0s\n",
      "     0     0  528.36163    0   44  784.00000  528.36163  32.6%     -    0s\n",
      "     0     0  528.36163    0   46  784.00000  528.36163  32.6%     -    0s\n",
      "     0     0  530.73456    0   47  784.00000  530.73456  32.3%     -    0s\n",
      "     0     0  531.83581    0   47  784.00000  531.83581  32.2%     -    0s\n",
      "     0     0  534.56433    0   48  784.00000  534.56433  31.8%     -    0s\n",
      "     0     0  569.28549    0   49  784.00000  569.28549  27.4%     -    0s\n",
      "     0     2  569.35011    0   49  784.00000  569.35011  27.4%     -    0s\n",
      "H 1002   608                     778.0000000  569.35011  26.8%  22.8    0s\n",
      "* 2085  1179              31     776.0000000  603.91514  22.2%  21.6    0s\n",
      "H 3476  1566                     774.0000000  634.24274  18.1%  20.9    0s\n",
      "H 3502  1503                     772.0000000  634.24274  17.8%  21.0    0s\n",
      "* 7293  1058              27     770.0000000  676.38309  12.2%  21.9    1s\n",
      "\n",
      "Cutting planes:\n",
      "  MIR: 23\n",
      "  Zero half: 1\n",
      "\n",
      "Explored 10251 nodes (219957 simplex iterations) in 1.21 seconds (2.08 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 10: 770 772 774 ... 1350\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 7.700000000000e+02, best bound 7.700000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 21055, time in user-callback 0.02 sec\n",
      "#########################################\n",
      "Graph-SCP Time:  1.2801659107208252\n",
      "Random Time: 0 0 0\n",
      "Gurobi Time: 3.1414\n",
      "Greedy Time: 0.0017549991607666016\n",
      "\n",
      "Optimal Objective:  770\n",
      "Gurobi Objective:  770\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  17505.0\n",
      "Gurobi Node Count :  10251.0\n",
      "\n",
      "Num Times Gurobi Run :  3\n",
      "Problem Size :  76 / 250\n",
      "#########################################\n",
      "################################\n",
      "#############1/5################\n",
      "################################\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 250 columns and 9788 nonzeros\n",
      "Model fingerprint: 0x23165593\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Found heuristic solution: objective 1884.0000000\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Found heuristic solution: objective 1641.0000000\n",
      "\n",
      "Root relaxation: objective 4.794783e+02, 213 iterations, 0.01 seconds (0.01 work units)\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "     0     0  479.47832    0   67 1641.00000  479.47832  70.8%     -    0s\n",
      "H    0     0                    1637.0000000  479.47832  70.7%     -    0s\n",
      "H    0     0                    1270.0000000  479.47832  62.2%     -    0s\n",
      "H    0     0                     930.0000000  479.47832  48.4%     -    0s\n",
      "     0     0  497.20891    0   72  930.00000  497.20891  46.5%     -    0s\n",
      "     0     0  532.83323    0   62  930.00000  532.83323  42.7%     -    0s\n",
      "     0     0  533.28078    0   66  930.00000  533.28078  42.7%     -    0s\n",
      "     0     0  533.28078    0   67  930.00000  533.28078  42.7%     -    0s\n",
      "     0     0  533.28078    0   67  930.00000  533.28078  42.7%     -    0s\n",
      "     0     0  533.28078    0   65  930.00000  533.28078  42.7%     -    0s\n",
      "     0     0  533.28078    0   66  930.00000  533.28078  42.7%     -    0s\n",
      "     0     0  533.28349    0   66  930.00000  533.28349  42.7%     -    0s\n",
      "     0     0  534.88838    0   65  930.00000  534.88838  42.5%     -    0s\n",
      "     0     0  534.88838    0   65  930.00000  534.88838  42.5%     -    0s\n",
      "     0     0  534.88838    0   67  930.00000  534.88838  42.5%     -    0s\n",
      "     0     0  534.88838    0   65  930.00000  534.88838  42.5%     -    0s\n",
      "     0     0  534.88838    0   66  930.00000  534.88838  42.5%     -    0s\n",
      "     0     0  563.22433    0   65  930.00000  563.22433  39.4%     -    0s\n",
      "     0     2  563.22433    0   65  930.00000  563.22433  39.4%     -    0s\n",
      "H  384   344                     890.0000000  563.22433  36.7%  29.2    0s\n",
      "H  386   344                     889.0000000  563.22433  36.6%  29.2    0s\n",
      "H 1669  1223                     875.0000000  563.22433  35.6%  29.1    0s\n",
      "*27147  9913              20     871.0000000  752.61197  13.6%  29.9    3s\n",
      "*28837  8722              27     849.0000000  756.14716  10.9%  29.8    3s\n",
      "\n",
      "Cutting planes:\n",
      "  Gomory: 3\n",
      "  MIR: 17\n",
      "  Zero half: 1\n",
      "\n",
      "Explored 42675 nodes (1209326 simplex iterations) in 4.85 seconds (10.36 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 10: 849 871 875 ... 1884\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 8.490000000000e+02, best bound 8.490000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 86348, time in user-callback 0.08 sec\n",
      "502\n",
      "Set parameter Presolve to value 0                         \n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 0 columns and 0 nonzeros\n",
      "Model fingerprint: 0xdb7c14f9\n",
      "Coefficient statistics:\n",
      "  Matrix range     [0e+00, 0e+00]\n",
      "  Objective range  [0e+00, 0e+00]\n",
      "  Bounds range     [0e+00, 0e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Iteration    Objective       Primal Inf.    Dual Inf.      Time\n",
      "       0    0.0000000e+00   2.000000e+02   0.000000e+00      0s\n",
      "\n",
      "Solved in 0 iterations and 0.02 seconds (0.00 work units)\n",
      "Infeasible model\n",
      "\n",
      "User-callback calls 19, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 0 columns and 0 nonzeros\n",
      "Model fingerprint: 0xdb7c14f9\n",
      "Coefficient statistics:\n",
      "  Matrix range     [0e+00, 0e+00]\n",
      "  Objective range  [0e+00, 0e+00]\n",
      "  Bounds range     [0e+00, 0e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Iteration    Objective       Primal Inf.    Dual Inf.      Time\n",
      "       0    0.0000000e+00   2.000000e+02   0.000000e+00      0s\n",
      "\n",
      "Solved in 0 iterations and 0.01 seconds (0.00 work units)\n",
      "Infeasible model\n",
      "\n",
      "User-callback calls 19, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 75 columns and 4569 nonzeros\n",
      "Model fingerprint: 0x00550351\n",
      "Variable types: 0 continuous, 75 integer (75 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R0 by 1.000000000\n",
      "\n",
      "Found heuristic solution: objective 1650.0000000\n",
      "Variable types: 0 continuous, 75 integer (75 binary)\n",
      "Found heuristic solution: objective 1255.0000000\n",
      "\n",
      "Root relaxation: objective 4.875064e+02, 171 iterations, 0.00 seconds (0.01 work units)\n",
      "Another try with MIP start\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "H    0     0                    1004.0000000  487.50638  51.4%     -    0s\n",
      "     0     0  487.50638    0   57 1004.00000  487.50638  51.4%     -    0s\n",
      "H    0     0                     988.0000000  487.50638  50.7%     -    0s\n",
      "H    0     0                     979.0000000  487.50638  50.2%     -    0s\n",
      "H    0     0                     884.0000000  487.50638  44.9%     -    0s\n",
      "     0     0  503.10700    0   57  884.00000  503.10700  43.1%     -    0s\n",
      "     0     0  504.98088    0   55  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   53  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   57  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   56  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   55  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   56  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   56  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   56  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   55  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   57  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  504.98088    0   57  884.00000  504.98088  42.9%     -    0s\n",
      "     0     0  511.00667    0   54  884.00000  511.00667  42.2%     -    0s\n",
      "     0     0  528.41551    0   57  884.00000  528.41551  40.2%     -    0s\n",
      "H    0     0                     867.0000000  528.41551  39.1%     -    0s\n",
      "     0     0  573.51666    0   56  867.00000  573.51666  33.9%     -    0s\n",
      "     0     0  576.42445    0   54  867.00000  576.42445  33.5%     -    0s\n",
      "     0     0  576.42516    0   58  867.00000  576.42516  33.5%     -    0s\n",
      "     0     0  576.42516    0   56  867.00000  576.42516  33.5%     -    0s\n",
      "     0     0  580.02753    0   58  867.00000  580.02753  33.1%     -    0s\n",
      "     0     0  580.02753    0   57  867.00000  580.02753  33.1%     -    0s\n",
      "     0     0  580.02753    0   58  867.00000  580.02753  33.1%     -    0s\n",
      "     0     0  580.17299    0   57  867.00000  580.17299  33.1%     -    0s\n",
      "     0     0  580.17299    0   57  867.00000  580.17299  33.1%     -    0s\n",
      "     0     0  583.35569    0   56  867.00000  583.35569  32.7%     -    0s\n",
      "     0     0  585.69132    0   56  867.00000  585.69132  32.4%     -    0s\n",
      "     0     0  585.69132    0   57  867.00000  585.69132  32.4%     -    0s\n",
      "     0     0  587.20585    0   56  867.00000  587.20585  32.3%     -    0s\n",
      "     0     0  587.20585    0   55  867.00000  587.20585  32.3%     -    0s\n",
      "     0     0  587.20585    0   56  867.00000  587.20585  32.3%     -    0s\n",
      "     0     0  587.52506    0   57  867.00000  587.52506  32.2%     -    0s\n",
      "     0     0  587.52506    0   55  867.00000  587.52506  32.2%     -    0s\n",
      "     0     0  587.52506    0   57  867.00000  587.52506  32.2%     -    0s\n",
      "     0     0  590.22281    0   57  867.00000  590.22281  31.9%     -    0s\n",
      "     0     0  591.47838    0   57  867.00000  591.47838  31.8%     -    0s\n",
      "     0     0  614.84799    0   57  867.00000  614.84799  29.1%     -    0s\n",
      "     0     2  614.84799    0   57  867.00000  614.84799  29.1%     -    0s\n",
      "* 6906  1988              22     866.0000000  708.72634  18.2%  22.8    1s\n",
      "* 9903  2070              18     852.0000000  745.22868  12.5%  21.8    1s\n",
      "*15292    25              18     851.0000000  823.37988  3.25%  20.8    1s\n",
      "*15690     1              24     849.0000000  840.02124  1.06%  20.6    1s\n",
      "\n",
      "Explored 15722 nodes (323634 simplex iterations) in 1.51 seconds (2.54 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 10: 849 851 852 ... 1255\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 8.490000000000e+02, best bound 8.490000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 32027, time in user-callback 0.03 sec\n",
      "#########################################\n",
      "Graph-SCP Time:  1.585658073425293\n",
      "Random Time: 0 0 0\n",
      "Gurobi Time: 4.8521\n",
      "Greedy Time: 0.0016999244689941406\n",
      "\n",
      "Optimal Objective:  849\n",
      "Gurobi Objective:  849\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  42675.0\n",
      "Gurobi Node Count :  15722.0\n",
      "\n",
      "Num Times Gurobi Run :  3\n",
      "Problem Size :  75 / 250\n",
      "#########################################\n",
      "################################\n",
      "#############2/5################\n",
      "################################\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 250 columns and 9782 nonzeros\n",
      "Model fingerprint: 0x95211dce\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Found heuristic solution: objective 1531.0000000\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Found heuristic solution: objective 1376.0000000\n",
      "\n",
      "Root relaxation: objective 4.704933e+02, 195 iterations, 0.00 seconds (0.01 work units)\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "     0     0  470.49325    0   57 1376.00000  470.49325  65.8%     -    0s\n",
      "H    0     0                    1307.0000000  470.49325  64.0%     -    0s\n",
      "H    0     0                    1271.0000000  470.49325  63.0%     -    0s\n",
      "H    0     0                     839.0000000  470.49325  43.9%     -    0s\n",
      "     0     0  487.61455    0   69  839.00000  487.61455  41.9%     -    0s\n",
      "H    0     0                     793.0000000  487.61455  38.5%     -    0s\n",
      "     0     0  490.73924    0   67  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  490.73924    0   62  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  490.73924    0   62  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  490.73924    0   61  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  490.73924    0   62  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  490.73924    0   60  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  490.73924    0   60  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  490.73924    0   61  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  490.73924    0   59  793.00000  490.73924  38.1%     -    0s\n",
      "     0     0  500.60823    0   60  793.00000  500.60823  36.9%     -    0s\n",
      "     0     0  504.35987    0   64  793.00000  504.35987  36.4%     -    0s\n",
      "     0     0  505.66406    0   65  793.00000  505.66406  36.2%     -    0s\n",
      "     0     0  505.97731    0   63  793.00000  505.97731  36.2%     -    0s\n",
      "     0     0  505.97731    0   65  793.00000  505.97731  36.2%     -    0s\n",
      "     0     0  506.93261    0   68  793.00000  506.93261  36.1%     -    0s\n",
      "     0     0  506.93261    0   68  793.00000  506.93261  36.1%     -    0s\n",
      "     0     0  506.93261    0   68  793.00000  506.93261  36.1%     -    0s\n",
      "     0     0  506.93261    0   68  793.00000  506.93261  36.1%     -    0s\n",
      "     0     0  506.93261    0   67  793.00000  506.93261  36.1%     -    0s\n",
      "     0     0  506.93261    0   66  793.00000  506.93261  36.1%     -    0s\n",
      "     0     0  506.93261    0   66  793.00000  506.93261  36.1%     -    0s\n",
      "     0     0  546.83349    0   63  793.00000  546.83349  31.0%     -    0s\n",
      "     0     2  546.83349    0   63  793.00000  546.83349  31.0%     -    0s\n",
      "\n",
      "Cutting planes:\n",
      "  Gomory: 2\n",
      "  MIR: 21\n",
      "  Zero half: 1\n",
      "\n",
      "Explored 14138 nodes (416457 simplex iterations) in 3.09 seconds (6.67 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 6: 793 839 1271 ... 1531\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 7.930000000000e+02, best bound 7.930000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 29355, time in user-callback 0.03 sec\n",
      "502\n",
      "Set parameter Presolve to value 0                         \n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 0 columns and 0 nonzeros\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Model fingerprint: 0xdb7c14f9\n",
      "Coefficient statistics:\n",
      "  Matrix range     [0e+00, 0e+00]\n",
      "  Objective range  [0e+00, 0e+00]\n",
      "  Bounds range     [0e+00, 0e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Iteration    Objective       Primal Inf.    Dual Inf.      Time\n",
      "       0    0.0000000e+00   2.000000e+02   0.000000e+00      0s\n",
      "\n",
      "Solved in 0 iterations and 0.01 seconds (0.00 work units)\n",
      "Infeasible model\n",
      "\n",
      "User-callback calls 19, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 3 columns and 235 nonzeros\n",
      "Model fingerprint: 0xce0790ae\n",
      "Variable types: 0 continuous, 3 integer (3 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 1e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R0 by 1.000000000\n",
      "\n",
      "\n",
      "Explored 0 nodes (0 simplex iterations) in 0.01 seconds (0.00 work units)\n",
      "Thread count was 1 (of 10 available processors)\n",
      "\n",
      "Solution count 0\n",
      "\n",
      "Model is infeasible\n",
      "Best objective -, best bound -, gap -\n",
      "\n",
      "User-callback calls 27, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 75 columns and 4913 nonzeros\n",
      "Model fingerprint: 0x3b6d935e\n",
      "Variable types: 0 continuous, 75 integer (75 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R0 by 1.000000000\n",
      "\n",
      "Found heuristic solution: objective 1515.0000000\n",
      "Variable types: 0 continuous, 75 integer (75 binary)\n",
      "Found heuristic solution: objective 1307.0000000\n",
      "\n",
      "Root relaxation: objective 4.783460e+02, 148 iterations, 0.01 seconds (0.00 work units)\n",
      "Another try with MIP start\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "H    0     0                     839.0000000  478.34600  43.0%     -    0s\n",
      "     0     0  478.34600    0   45  839.00000  478.34600  43.0%     -    0s\n",
      "     0     0  499.66990    0   47  839.00000  499.66990  40.4%     -    0s\n",
      "     0     0  524.29373    0   47  839.00000  524.29373  37.5%     -    0s\n",
      "     0     0  524.29373    0   47  839.00000  524.29373  37.5%     -    0s\n",
      "     0     0  530.92898    0   46  839.00000  530.92898  36.7%     -    0s\n",
      "     0     0  531.30430    0   46  839.00000  531.30430  36.7%     -    0s\n",
      "     0     0  540.58108    0   48  839.00000  540.58108  35.6%     -    0s\n",
      "     0     0  541.56574    0   49  839.00000  541.56574  35.5%     -    0s\n",
      "H    0     0                     811.0000000  591.75244  27.0%     -    0s\n",
      "     0     0  591.75244    0   50  811.00000  591.75244  27.0%     -    0s\n",
      "     0     0  593.89317    0   48  811.00000  593.89317  26.8%     -    0s\n",
      "     0     0  595.92932    0   49  811.00000  595.92932  26.5%     -    0s\n",
      "     0     0  597.83247    0   47  811.00000  597.83247  26.3%     -    0s\n",
      "     0     0  599.11624    0   50  811.00000  599.11624  26.1%     -    0s\n",
      "     0     0  599.20903    0   48  811.00000  599.20903  26.1%     -    0s\n",
      "     0     0  599.21832    0   46  811.00000  599.21832  26.1%     -    0s\n",
      "     0     0  599.21832    0   46  811.00000  599.21832  26.1%     -    0s\n",
      "     0     0  599.43617    0   47  811.00000  599.43617  26.1%     -    0s\n",
      "     0     0  599.95185    0   48  811.00000  599.95185  26.0%     -    0s\n",
      "     0     0  602.87188    0   50  811.00000  602.87188  25.7%     -    0s\n",
      "     0     0  603.37592    0   49  811.00000  603.37592  25.6%     -    0s\n",
      "     0     0  617.06032    0   48  811.00000  617.06032  23.9%     -    0s\n",
      "     0     2  617.06032    0   48  811.00000  617.06032  23.9%     -    0s\n",
      "*  323   293              32     793.0000000  617.06032  22.2%  25.1    0s\n",
      "\n",
      "Cutting planes:\n",
      "  Gomory: 3\n",
      "  MIR: 49\n",
      "  Zero half: 2\n",
      "\n",
      "Explored 4935 nodes (98747 simplex iterations) in 0.58 seconds (0.82 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 5: 793 811 839 ... 1515\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 7.930000000000e+02, best bound 7.930000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 10187, time in user-callback 0.02 sec\n",
      "#########################################\n",
      "Graph-SCP Time:  0.6619358062744141\n",
      "Random Time: 0 0 0\n",
      "Gurobi Time: 3.0899\n",
      "Greedy Time: 0.0017387866973876953\n",
      "\n",
      "Optimal Objective:  793\n",
      "Gurobi Objective:  793\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  14138.0\n",
      "Gurobi Node Count :  4935.0\n",
      "\n",
      "Num Times Gurobi Run :  3\n",
      "Problem Size :  75 / 250\n",
      "#########################################\n",
      "################################\n",
      "#############3/5################\n",
      "################################\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 250 columns and 10617 nonzeros\n",
      "Model fingerprint: 0xb252d6e5\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "Found heuristic solution: objective 1630.0000000\n",
      "Variable types: 0 continuous, 250 integer (250 binary)\n",
      "Found heuristic solution: objective 1265.0000000\n",
      "\n",
      "Root relaxation: objective 4.464195e+02, 217 iterations, 0.01 seconds (0.01 work units)\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "     0     0  446.41953    0   59 1265.00000  446.41953  64.7%     -    0s\n",
      "H    0     0                    1180.0000000  446.41953  62.2%     -    0s\n",
      "H    0     0                    1066.0000000  446.41953  58.1%     -    0s\n",
      "H    0     0                    1033.0000000  446.41953  56.8%     -    0s\n",
      "H    0     0                     884.0000000  446.41953  49.5%     -    0s\n",
      "     0     0  469.35270    0   67  884.00000  469.35270  46.9%     -    0s\n",
      "H    0     0                     815.0000000  469.35270  42.4%     -    0s\n",
      "     0     0  487.25097    0   60  815.00000  487.25097  40.2%     -    0s\n",
      "     0     0  487.25097    0   59  815.00000  487.25097  40.2%     -    0s\n",
      "     0     0  487.43007    0   60  815.00000  487.43007  40.2%     -    0s\n",
      "     0     0  487.50008    0   61  815.00000  487.50008  40.2%     -    0s\n",
      "     0     0  487.50008    0   59  815.00000  487.50008  40.2%     -    0s\n",
      "     0     0  487.50008    0   58  815.00000  487.50008  40.2%     -    0s\n",
      "     0     0  487.50008    0   59  815.00000  487.50008  40.2%     -    0s\n",
      "     0     0  488.30840    0   58  815.00000  488.30840  40.1%     -    0s\n",
      "     0     0  517.22116    0   61  815.00000  517.22116  36.5%     -    0s\n",
      "     0     0  521.28001    0   61  815.00000  521.28001  36.0%     -    0s\n",
      "     0     0  521.36328    0   58  815.00000  521.36328  36.0%     -    0s\n",
      "     0     0  525.01570    0   62  815.00000  525.01570  35.6%     -    0s\n",
      "     0     0  530.65225    0   63  815.00000  530.65225  34.9%     -    0s\n",
      "     0     0  530.65225    0   62  815.00000  530.65225  34.9%     -    0s\n",
      "     0     0  533.49767    0   61  815.00000  533.49767  34.5%     -    0s\n",
      "     0     0  533.49767    0   64  815.00000  533.49767  34.5%     -    0s\n",
      "     0     0  533.83331    0   63  815.00000  533.83331  34.5%     -    0s\n",
      "     0     0  534.59173    0   60  815.00000  534.59173  34.4%     -    0s\n",
      "     0     0  534.59173    0   63  815.00000  534.59173  34.4%     -    0s\n",
      "     0     0  536.57325    0   61  815.00000  536.57325  34.2%     -    0s\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "     0     0  536.57325    0   63  815.00000  536.57325  34.2%     -    0s\n",
      "     0     0  564.68629    0   61  815.00000  564.68629  30.7%     -    0s\n",
      "     0     2  564.68629    0   61  815.00000  564.68629  30.7%     -    0s\n",
      "*   92    78              10     791.0000000  564.68629  28.6%  40.1    0s\n",
      "*  171   142              15     789.0000000  564.68629  28.4%  34.0    0s\n",
      "* 5313  2203              19     781.0000000  622.64577  20.3%  34.2    2s\n",
      "\n",
      "Cutting planes:\n",
      "  Gomory: 2\n",
      "  MIR: 14\n",
      "  Zero half: 2\n",
      "\n",
      "Explored 23287 nodes (676005 simplex iterations) in 3.90 seconds (7.98 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 10: 781 789 791 ... 1630\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 7.810000000000e+02, best bound 7.810000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 47830, time in user-callback 0.05 sec\n",
      "502\n",
      "Set parameter Presolve to value 0                         \n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 1 columns and 80 nonzeros\n",
      "Model fingerprint: 0x0395ae4d\n",
      "Variable types: 0 continuous, 1 integer (1 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 1e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R1 by 1.000000000\n",
      "\n",
      "\n",
      "Explored 0 nodes (0 simplex iterations) in 0.00 seconds (0.00 work units)\n",
      "Thread count was 1 (of 10 available processors)\n",
      "\n",
      "Solution count 0\n",
      "\n",
      "Model is infeasible\n",
      "Best objective -, best bound -, gap -\n",
      "\n",
      "User-callback calls 27, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 6 columns and 478 nonzeros\n",
      "Model fingerprint: 0x6c84f796\n",
      "Variable types: 0 continuous, 6 integer (6 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 1e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R1 by 1.000000000\n",
      "\n",
      "\n",
      "Explored 0 nodes (0 simplex iterations) in 0.00 seconds (0.00 work units)\n",
      "Thread count was 1 (of 10 available processors)\n",
      "\n",
      "Solution count 0\n",
      "\n",
      "Model is infeasible\n",
      "Best objective -, best bound -, gap -\n",
      "\n",
      "User-callback calls 27, time in user-callback 0.00 sec\n",
      "Set parameter Presolve to value 0\n",
      "Set parameter TimeLimit to value 3600\n",
      "Gurobi Optimizer version 10.0.0 build v10.0.0rc2 (mac64[arm])\n",
      "\n",
      "CPU model: Apple M1 Max\n",
      "Thread count: 10 physical cores, 10 logical processors, using up to 10 threads\n",
      "\n",
      "Optimize a model with 200 rows, 76 columns and 5055 nonzeros\n",
      "Model fingerprint: 0xcdc41e79\n",
      "Variable types: 0 continuous, 76 integer (76 binary)\n",
      "Coefficient statistics:\n",
      "  Matrix range     [1e+00, 1e+00]\n",
      "  Objective range  [1e+02, 2e+02]\n",
      "  Bounds range     [1e+00, 1e+00]\n",
      "  RHS range        [1e+00, 1e+00]\n",
      "\n",
      "User MIP start did not produce a new incumbent solution\n",
      "User MIP start violates constraint R1 by 1.000000000\n",
      "\n",
      "Found heuristic solution: objective 1616.0000000\n",
      "Variable types: 0 continuous, 76 integer (76 binary)\n",
      "Found heuristic solution: objective 1369.0000000\n",
      "\n",
      "Root relaxation: objective 4.570024e+02, 159 iterations, 0.00 seconds (0.01 work units)\n",
      "Another try with MIP start\n",
      "\n",
      "    Nodes    |    Current Node    |     Objective Bounds      |     Work\n",
      " Expl Unexpl |  Obj  Depth IntInf | Incumbent    BestBd   Gap | It/Node Time\n",
      "\n",
      "H    0     0                     920.0000000  457.00245  50.3%     -    0s\n",
      "     0     0  457.00245    0   46  920.00000  457.00245  50.3%     -    0s\n",
      "H    0     0                     916.0000000  457.00245  50.1%     -    0s\n",
      "H    0     0                     815.0000000  457.00245  43.9%     -    0s\n",
      "     0     0  479.26723    0   51  815.00000  479.26723  41.2%     -    0s\n",
      "     0     0  519.06130    0   51  815.00000  519.06130  36.3%     -    0s\n",
      "     0     0  519.32223    0   50  815.00000  519.32223  36.3%     -    0s\n",
      "     0     0  526.14406    0   51  815.00000  526.14406  35.4%     -    0s\n",
      "     0     0  526.14406    0   51  815.00000  526.14406  35.4%     -    0s\n",
      "     0     0  534.84742    0   50  815.00000  534.84742  34.4%     -    0s\n",
      "     0     0  559.76133    0   50  815.00000  559.76133  31.3%     -    0s\n",
      "     0     0  560.69370    0   52  815.00000  560.69370  31.2%     -    0s\n",
      "     0     0  561.11330    0   49  815.00000  561.11330  31.2%     -    0s\n",
      "     0     0  561.11330    0   51  815.00000  561.11330  31.2%     -    0s\n",
      "     0     0  561.16980    0   51  815.00000  561.16980  31.1%     -    0s\n",
      "     0     0  561.16980    0   51  815.00000  561.16980  31.1%     -    0s\n",
      "     0     0  561.74564    0   49  815.00000  561.74564  31.1%     -    0s\n",
      "     0     0  562.25570    0   51  815.00000  562.25570  31.0%     -    0s\n",
      "     0     0  562.25570    0   52  815.00000  562.25570  31.0%     -    0s\n",
      "     0     0  562.51217    0   52  815.00000  562.51217  31.0%     -    0s\n",
      "     0     0  562.51217    0   51  815.00000  562.51217  31.0%     -    0s\n",
      "     0     0  562.96760    0   51  815.00000  562.96760  30.9%     -    0s\n",
      "     0     0  562.96760    0   50  815.00000  562.96760  30.9%     -    0s\n",
      "     0     0  579.41877    0   53  815.00000  579.41877  28.9%     -    0s\n",
      "H    0     0                     800.0000000  579.41877  27.6%     -    0s\n",
      "     0     2  579.41877    0   53  800.00000  579.41877  27.6%     -    0s\n",
      "H 1007   674                     789.0000000  579.41877  26.6%  21.7    0s\n",
      "* 9594  1399              26     787.0000000  700.47325  11.0%  21.3    1s\n",
      "*11449   788              24     781.0000000  731.67839  6.32%  21.2    1s\n",
      "\n",
      "Cutting planes:\n",
      "  Gomory: 1\n",
      "  MIR: 21\n",
      "  Zero half: 2\n",
      "\n",
      "Explored 12634 nodes (262779 simplex iterations) in 1.45 seconds (2.41 work units)\n",
      "Thread count was 10 (of 10 available processors)\n",
      "\n",
      "Solution count 9: 781 787 789 ... 1616\n",
      "\n",
      "Optimal solution found (tolerance 1.00e-04)\n",
      "Best objective 7.810000000000e+02, best bound 7.810000000000e+02, gap 0.0000%\n",
      "\n",
      "User-callback calls 25903, time in user-callback 0.03 sec\n",
      "#########################################\n",
      "Graph-SCP Time:  1.486314058303833\n",
      "Random Time: 0 0 0\n",
      "Gurobi Time: 3.9056\n",
      "Greedy Time: 0.0016629695892333984\n",
      "\n",
      "Optimal Objective:  781\n",
      "Gurobi Objective:  781\n",
      "Random Objective:  0 0 0\n",
      "\n",
      "Gurobi Node Count :  23287.0\n",
      "Gurobi Node Count :  12634.0\n",
      "\n",
      "Num Times Gurobi Run :  3\n",
      "Problem Size :  76 / 250\n",
      "#########################################\n",
      "################################\n",
      "#############4/5################\n",
      "################################\n"
     ]
    }
   ],
   "source": [
    "results, inst = test_set_cover(test_instances = test_1[-5:],\n",
    "                                 num_runs = 10,\n",
    "                                 cost_type = 'custom',\n",
    "                                 model_path = './models/set_cover_sage_unsup.pkl',\n",
    "                                 m_type = 'sage',\n",
    "                                 feats = ['costs', 'cover', 'ppr', 'bi_sf', 'hyp_feat'],\n",
    "                                 init_sol_perc = 90,\n",
    "                                 perc_start = 85,\n",
    "                                 dec_size = 25, \n",
    "                                 obj_threshold = 1, \n",
    "                                 time_limit = 3600, \n",
    "                                 run_random = False, \n",
    "                                 gnn_solver = 'gurobi', \n",
    "                                 neurons = 512)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e44e5bd1",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "49e39d46",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.8"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}